Abstract
Seismic inversion is one of the key techniques used for reservoir characterization. Depth-domain seismic inversion eliminates the cumulative errors associated with depth-to-time and time-to-depth conversions, thus providing geologists and reservoir engineers with an intuitive basis for geological interpretation. The method has received increasing attention in the field of reservoir characterization. Extracting accurate depth-domain seismic wavelets is a prerequisite for successful depth-domain seismic inversion. However, the depth-domain wavelet is velocity-dependent and exhibits significant non-stationarity, which leads to the failure of seismic wavelet estimation methods based on the stationary convolutional model. To this end, we propose a modified wavenumber-domain unscaled S-transform (MWUST) method to accomplish accurate estimation of depth-domain seismic wavelets. The proposed method enhances the accuracy of wavenumber components by removing the linear wavenumber-dependent term from the S-transform. Furthermore, it introduces slope and intercept parameters to improve the depth resolution at low wavenumbers, thereby yielding a more reliable depth–wavenumber spectrum. Subsequently, the relationship between the non-stationary depth-domain seismic wavelet and the depth–wavenumber spectrum is established, allowing for the accurate extraction of non-stationary wavelets under the assumption that the depth-domain reflectivity is a random sequence. Synthetic and real data applications have been used to verify the effectiveness of the proposed method.
1 Introduction
With continuous advancements in oil and gas exploration, the demands for accurate reservoir characterization are becoming increasingly stringent. Seismic inversion is a key technique used for obtaining elastic parameters for reservoir characterization (; ; ; ; ). Many depth migration methods have been developed and are widely used to generate depth-domain seismic images, leading to a growing need for depth-domain seismic inversion of subsurface elastic properties for reservoir characterization (; ; ; ; ; ). Depth-domain seismic inversion avoids the cumulative errors associated with depth-to-time and time-to-depth conversions, providing geologically significant results for geologists and reservoir engineers. However, depth-domain seismic inversion is different from time-domain seismic inversion, in that methods applicable in the time domain cannot be directly applied to the depth domain. The primary difference is that depth-domain seismic data exhibit strong spectral variability. This is because the depth-domain seismic wavelet depends on velocity, i.e., higher velocity results in longer waveforms, and vice versa. This makes it challenging to assume stationarity in depth-domain seismic data, which can lead to erroneous inversion results. In order to obtain reliable results, the accurate extraction of the depth-domain seismic wavelet is a prerequisite ().
Typically, seismic wavelets can be estimated by statistical and deterministic approaches (). Statistical methods estimate seismic wavelets using only seismic data (; ; ; ; ; ; ), whereas deterministic methods rely on a combination of seismic and well-log data to estimate the wavelets (; ; ; ; ; ). Given the limited well data in the field region, the time–frequency decomposition method, which is a statistical approach, serves as a suitable choice for estimating non-stationary seismic wavelets. extended the S-transform (ST)-based time–frequency decomposition method from the time–frequency domain to the depth–wavenumber domain, thereby completing the estimation of non-stationary seismic wavelets in the depth domain. successfully estimated the depth-domain seismic wavelet using the ST and used it for depth-domain pre-stack seismic inversion, thereby obtaining reliable depth-domain elastic parameters including P- and S-wave velocity and density. developed a depth–wavenumber decomposition technique that applies the ST to depth-domain pre-stack angle gathers, generating depth and angle-variant wavelets for inversion. introduced a generalized unscaled ST for the spectral decomposition of depth-domain seismic data, estimating seismic wavelets by utilizing local spectral information at each depth sampling point.
The spectral decomposition method entails performing depth–wavenumber spectral decomposition on depth-domain seismic data and estimating depth-variant wavelets from the resulting depth–wavenumber spectrum. As a hybrid of the short-time Fourier transform and the continuous wavelet transform, the ST combines the advantages of both, offering adaptive resolution adjustment and lossless reversible transformation. The generalized S-transform (GST) was later developed for specific applications. However, the linear frequency-dependent term in the ST or GST causes its frequency distribution to deviate from the true values. This leads to incorrect estimation of seismic wavelets, which affects the accuracy of the subsequent depth-domain seismic inversion. developed an unscaled ST (UST) algorithm to reduce the bias in frequency components by removing the linear frequency-dependent term of the ST. derived the frequency domain form of the UST that preserves the signal amplitude. However, the removal of the linear frequency-dependent term reduces their temporal resolution at low frequencies. introduced slope and intercept parameters into the generalized ST to achieve the desired time–frequency resolution. further improved the temporal resolution of the time–frequency spectrum by modifying the basis in the study by , while ensuring the preservation of frequency characteristics. Existing wavelet extraction methods based on spectral decomposition are constrained by the limitations of the spectral decomposition algorithm, which reduces the accuracy of wavelet extraction. Inspired by this, we combine the advantages of the aforementioned time–frequency decomposition methods to re-derive a modified frequency-domain unscaled ST algorithm and extend it to the wavenumber domain (MWUST), thereby achieving an accurate estimation of non-stationary seismic wavelets in the depth domain. The effectiveness of the proposed method is verified through a series of numerical experiments and by comparison with traditional ST ()- and UST ()-based depth-domain wavelet extraction methods.
2 Methods
Depth-domain seismic inversion requires a reliable depth-domain seismic wavelet as a prerequisite. In order to obtain a reliable depth-domain seismic wavelet, a reliable depth–wavenumber decomposition method is required. In this section, we first derive a modified frequency-domain unscaled ST algorithm and extend it to the wavenumber domain. Then, we describe in detail the estimation of the depth-domain seismic wavelet based on the depth–wavenumber spectrum obtained through the proposed method.
2.1 The modified wavenumber-domain unscaled S transform
For a time-domain signal , the time–frequency spectrum based on the ST is expressed as Equation 1 ()where is the Gaussian window of the specific form , is a translation factor to control the position of the Gaussian window on the time axis , represents the frequency, and is imaginary units.
Then, the depth–wavenumber spectrum of a depth-domain signal based on the ST can be expressed as Equation 2where represents the depth-domain Gaussian window of the specific form , is a translation factor to control the position of the Gaussian window on the depth axis , and represents the wavenumber.
To reduce the bias of the frequency components in the time–frequency spectrum based on the ST, the time-domain UST removes the linear frequency-dependent term of . Then, the depth-domain UST is derived as
Then, the expression to the right of the equal sign in Equation 3 in the wavenumber domain (i.e., WUST) can be written aswhere is the Fourier transform of . denotes the translating wavenumber ().
However, the UST sacrifices time/depth resolution at low frequencies/low wavenumbers in order to obtain a reliable frequency/wavenumber distribution. To overcome this shortcoming, slope () and intercept () parameters are introduced to obtain the desired depth–wavenumber resolution. The equation of the modified UST in the depth domain is
Based on the convolution theorem, the expression to the right of the equal sign in Equation 5 in the wavenumber domain (i.e., MWUST) can be written aswhere is the Fourier transform of . The constant average of the signal is put into zero wavenumber, which ensures the feasibility of the inverse transform. By integrating and inverse Fourier-transforming along different axes, we can reconstruct the original signal. It is worth noting that Equation 6 degenerates to Equation 4 when and .
There are two parameters (i.e., A and B) in the proposed method that need to be determined in advance. For different tasks and frequencies, the values of A and B may vary, which requires manual adjustment. One can use trial-and-error methods to obtain relatively optimal results. However, before that, we can use the full-window spatial or wavenumber width at half-maximum (FWHM) to roughly estimate parameters A and B (; ). For the depth and wavenumber domains, there are different expressions for A and B.
For the depth domain, A and B have the following expressions (Equation 7) ():where and represent specified spatial (or depth) FWHM resolutions at wavenumbers and , respectively.
For the wavenumber domain, A and B have the following expressions (Equation 8):where and represent specified wavenumber FWHM resolutions at wavenumbers and , respectively.
Figure 1 shows relationships between FWHM amplitudes and wavenumber for different A and B expressions. In Figure 1A, the width of the MWUST is narrower than that of the WUST in the low-wavenumber region when (i.e., orange-, green-, and gray-solid curves). It indicates that the MWUST has higher resolution in the low-wavenumber region, and the difference decreases as the wavenumber increases. When , shown as a purple solid curve, the width of the MWUST is broader than that of the WUST in the high-wavenumber region, indicated by larger depth FWHM values. The wavenumber FWHM values increase linearly as the wavenumber increases, as shown in Figure 1B. From Figure 1, it can also be noticed that A plays a greater role for FWHM values than B, indicating that A has a great influence on the depth–wavenumber resolution. B controls the starting depth resolution.
FIGURE 1
2.2 Depth-domain seismic wavelet estimation based on the depth–wavenumber spectrum
The depth-domain non-stationary convolution model is expressed aswhere represents depth-domain synthetic seismic records, represents the non-stationary depth-domain seismic wavelet, and represents the depth-domain reflectivity series. Equation 9 is transformed into the wavenumber domain with the following expression (Equation 10):where is the 2D Fourier transformation of . and are the Fourier transformation of and , respectively. If is the random reflectivity, we have ()where is a constant value, which can be determined from the well-side seismic traces. We then perform the inverse Fourier transform to Equation 11 to reproduce the original signal, as follows (Equation 12):
Combined with the original signal reconstruction method, the following expression is obtained:where represents the depth-domain seismic wavelet and represents the results after inverse Fourier transform along the depth of the depth–wavenumber spectrum (e.g., ). Then, the estimation of the depth-domain seismic wavelet can be achieved based on the depth–wavenumber spectrum calculated by the method proposed in Section 2.1.
3 Application
To validate the proposed approach, we apply it to synthetic experiments and real data. We first illustrate the differences between depth-domain and time-domain wavelets under different velocity model conditions and further apply the proposed method to the single-wavelet model. A three-wavelet model is then designed based on the single-wavelet model to verify the validity of the proposed method. Finally, well-log data are used for testing, and the extracted wavelet is used for depth-domain inversion to obtain satisfactory results.
3.1 Synthetic data test
We first experiment with the proposed method using a simple three-wavelet model, referencing the model used in . Figure 2A shows a 40-Hz zero-phase time-domain Ricker wavelet, and its corresponding amplitude spectrum is shown in Figure 2B. Taking this wavelet as the source, the corresponding depth-domain wavelets as they propagate through the different velocity strata are shown in Figure 2C, and their corresponding amplitude spectrums are shown in Figure 2D. It can be seen that as the velocity changes, the depth-domain wavelet is significantly stretched or compressed, i.e., the depth-domain wavelet is velocity-dependent. Since subsurface velocities are spatially variable, the interpretation and inversion of depth-domain seismic data can result in unreliable results if a constant wavelet is utilized.
FIGURE 2
Figure 3 shows the single-wavelet (i.e., sky-blue line shown in Figure 2C) depth–wavenumber spectrum experiments for different methods. It can be seen that WUST (Figure 3C) obtains a more accurate wavenumber spectrum than ST (Figure 3B) by removing the linear wavenumber-dependent term, but it performs poorly in the low-wavenumber region. By introducing parameters A and B, MWUST ensures that the wavenumber spectrum is accurate while performing well in the low-wavenumber region (Figures 3D–F). Figures 3D–F show that A plays a greater role in the depth–wavenumber spectrum than B, which is consistent with the conclusion reached in Figure 1. The amplitude spectra are extracted from the red dashed lines in Figure 3 to describe the wavenumber distribution in detail, as shown in Figure 4. Here, 1/3 and 10 are used for A and B in MWUST, respectively. It can be seen that the central frequency of the amplitude spectrum obtained based on the ST deviates from the reference frequency. However, MWUST matches the FT curve well, which further validates the advantages and effectiveness of the proposed method (i.e., MWUST) over ST and WUST.
FIGURE 3
FIGURE 4
Subsequently, a simple model in which the three reflection coefficients (Figure 5A) are related to the wavelet velocities mentioned in Figure 2C is used to perform wavelet estimation. Figures 6A–C show the depth–wavenumber spectrum of ST, WUST, and MWUST, where the black and red asterisks indicate the reference primary wavenumber of the signal at the position of the reflection coefficient and the primary wavenumber corresponding to the depth–wavenumber spectrum calculated by the different methods, respectively. It can be seen that the wavenumber distributions of WUST and MWUST match the reference values more closely than those of ST, but the depth–wavenumber spectrum of MWUST has a higher depth resolution in the low-wavenumber region (white arrows).
FIGURE 5
FIGURE 6
Figures 7A–C show the depth-variant wavelets extracted from the depth–wavenumber spectrum obtained using different methods. Based on the known reflection coefficients, the depth-domain wavelets extracted by the different methods are convolved with them to obtain reconstructed seismic records, as shown by the red line in Figures 5B–D. The seismic records reconstructed by the proposed method (i.e., MWUST) are best matched to the reference seismogram (black line in Figure 5D). Figure 8 shows the normalized errors of the reconstructed seismic records of the wavelet extracted by the different methods. The white dashed line indicates the exact solution, i.e., the closer the focus is to the white dashed line, the more accurate the reconstructed result is. As expected, the more accurate depth–wavenumber spectrum obtained by MWUST resulted in more accurate extracted depth-domain wavelets, thus minimizing the error between the reconstructed seismogram and the reference seismogram.
FIGURE 7
FIGURE 8
3.2 Field data example
A more realistic example is further used to test the effectiveness of the method. The well–seismic ties between velocity curves and depth-domain seismic data (black lines) are displayed in Figure 9. Figures 10A–C show the depth–wavenumber spectrum of depth-domain seismic data using ST, WUST, and MWUST, respectively. The depth-variant wavelets extracted using the different methods are shown in Figures 11A–C. The reconstructed seismograms (red lines) obtained by combining known velocity well-log data and depth-domain wavelets extracted using the different methods are shown in Figures 9B–D. In this test, the density is assumed to be a constant density that does not vary with depth. The seismic records reconstructed (red lines) by the three methods roughly match the reference seismic records (black lines). However, the reconstructed seismic records of the proposed method are more consistent with the reference seismic records in some aspects compared to the reconstructed records of the other two methods (e.g., green arrows).
FIGURE 9
FIGURE 10
FIGURE 11
We perform further inversion tests to demonstrate the impact of the extracted wavelet on subsequent seismic inversion and interpretation. The algorithm used here is a Bayesian-based inversion method, the details of which can be found in . In order to ensure a fair comparison of the inversion results, we perform the inversion by changing only the input wavelet, leaving the other inversion parameters unchanged. Figure 12 shows the inversion results obtained from the wavelet extracted by the different methods (i.e., ST, WUST, and MWUST) as input. It can be seen that the results obtained from the inversion of the wavelet extracted based on the ST show significant oscillations compared to those obtained from the inversion of the wavelets extracted by the other two methods. Although small differences are exhibited in the reconstructed seismic records (Figure 9B), the impact of the small errors on the inversion results is significant. The accuracy of the inversion results (green arrows in Figure 12) of the proposed method is greater due to the extraction of a more accurate depth-domain wavelet.
FIGURE 12
4 Discussion
Many depth migration techniques have been developed to create depth-domain seismic datasets, which are increasingly being utilized for oil and gas exploration and demonstrate greater advantages over time-domain datasets. This places higher requirements on direct deep-domain processing and interpretation techniques, including depth-domain seismic inversion. To obtain reliable results from depth-domain seismic inversion, it is necessary to have an accurate depth-domain wavelet. The wavelet has two key parameters, dominant frequency and phase, mostly obtained from the depth–wavenumber spectrum for the depth-domain seismic data. To this end, we developed a workflow for extracting depth-variant wavelets to accommodate the potential non-stationarity of seismic data in the depth domain. The method improves the precision of wavenumber components by removing the linear term related to the wavenumber in the S-transform. Additionally, it introduces slope and intercept parameters to enhance the depth resolution at low wavenumbers, resulting in a more reliable depth–wavenumber spectrum. Then, on the premise that the reflectivity in the depth domain is a random sequence, the non-stationary wavelet can be accurately extracted based on the obtained depth–wavenumber spectrum. However, the proposed method is not limited to the extraction of non-stationary wavelets in the depth domain. Depending on the difference in the input spectrum, the proposed method can be easily extended to the extraction of seismic wavelets from both stationary and non-stationary signals in the time domain.
examined the impact of varying angles on seismic wavelets and extended the spectral decomposition-based wavelet extraction method to pre-stack seismic data. In this study, we focus on achieving a more precise spectrum by re-deriving the expression for the depth–wavenumber spectrum to ensure the accuracy of the subsequent wavelet extraction. This improved method is tested and validated using post-stack seismic data. Indeed, seismic gathers exhibit angle dependence in actual exploration processes, and it is crucial to consider the effects of this angle dependence. In addition, attenuation and dispersion effects of seismic signals during propagation through subsurface media are not considered. For future research, the proposed depth-domain wavelet estimation method is suggested to be improved by incorporating the depth-variant effects of seismic attenuation and dispersion. It is hoped to obtain a wavelet matrix that can more accurately characterize seismic wave propagation and improve the accuracy of subsequent reservoir prediction and interpretation.
5 Conclusion
An accurate depth-domain wavelet plays a crucial role in the reliable interpretation and inversion of depth-domain seismic data. In this paper, we propose a MWUST method to accomplish accurate depth-domain seismic wavelet estimation. The method obtains more accurate wavenumber distribution with guaranteed depth and wavenumber resolution, thus providing a prerequisite for reliable depth-domain wavelet estimation. Numerical experiments demonstrate the advantages of the proposed method over ST- and WUST-based wavelet extraction methods. We further apply depth-domain wavelets extracted using different methods to perform inversion tests and show that even small errors can lead to unreliable inversion results. This illustrates the importance of extracting an accurate wavelet for subsequent seismic inversion and interpretation, as well as the effectiveness of the proposed method. Additionally, the proposed method can be easily extended to the wavelet extraction of arbitrary non-stationary seismic data, in addition to its application to depth-domain seismic wavelet extraction.
Statements
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Author contributions
YX: conceptualization, methodology, writing–original draft, and writing–review and editing. JZ: software, supervision, and writing–review and editing. GZ: investigation, validation, and writing–review and editing. DZ: supervision and writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was financially supported by the National Natural Science Foundation of China (42204108), the National Key Laboratory of Petroleum Resources and Engineering, China University of Petroleum, Beijing (PRE/open-2305), the Natural Science Foundation of Sichuan (2023NSFSC0768), and the Software Development Project for Concrete Defect Measurement and Processing at Multi-Layer Interfaces (2024-ZH004-N).
Conflict of interest
Authors GZ and DZ were employed by Sichuan Water Development Investigation, Design & Research Co., Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
BulandA.OmreH. J. (2003). Bayesian wavelet estimation from seismic and well data. Geophysics68 (6), 2000–2009. 10.1190/1.1635053
2
CaiR.SunC.YaoZ. A.LiS. (2024). A depth-variant seismic wavelet extraction method for basis pursuit inversion with an impedance trend constraint. Geophysics89 (3), R275–R286. 10.1190/geo2023-0255.1
3
DanielsenV.KarlssonT. (1984). Extraction of signatures from seismic and well data. First Break2 (4), 15–21. 10.3997/1365-2397.1984008
4
De MacedoI. A.De FigueiredoJ. J. S. J. (2020). On the seismic wavelet estimative and reflectivity recovering based on linear inversion: well-to-seismic tie on a real data set from Viking Graben, North Sea. Geophysics85 (5), D157–D165. 10.1190/geo2019-0183.1
5
GeorgeN. V.SahuS. S.MansinhaL.TiampoK.PandaG. (2009). Time localised band filtering using modified S-transform. IEEE68, 42–46. 10.1109/icsps.2009.63
6
GunningJ.GlinskyM. E. (2006). Wavelet extractor: a Bayesian well-tie and wavelet extraction program. Comput. Geosci.32 (5), 681–695. 10.1016/j.cageo.2005.10.001
7
KeX.ShiY.FuX.SongL.JingH.YangJ.et al (2023). The nth power Fourier spectrum analysis for the generalized seismic wavelets. IEEE Trans. Geoscience Remote Sens.61, 1–10. 10.1109/tgrs.2023.3243184
8
LaakeA.ZhangR.DengZ. (2019). Depth domain pre-stack seismic inversion with depth and angle variant wavelets. OnePetro, 1–5.
9
LazearG. D. J. (1993). Mixed-phase wavelet estimation using fourth-order cumulants. Geophysics58 (7), 1042–1051. 10.1190/1.1443480
10
LiD.CastagnaJ.GoloshubinG. (2016). Investigation of generalized S-transform analysis windows for time-frequency analysis of seismic reflection data. Geophysics81 (3), V235–V247. 10.1190/geo2015-0551.1
11
LiM.LuJ.StovasA.WangY. (2022). Joint PP, PS1, and PS2 AVA inversion of HTI media. IEEE Trans. Geoscience Remote Sens.60, 1–13. 10.1109/tgrs.2022.3170240
12
LiS.PengG.WangJ.TianT. (2020). Advantages and applications of seismic data interpretation in depth domain. OnePetro, 1–5.
13
LinesL. R.TreitelS. J. (1985). Wavelets, well logs and Wiener filters. First Break3 (8), 9–14. 10.3997/1365-2397.1985014
14
PaxtonA.KholiqA.ZhangL.ShadlowJ.ShadrinaM. (2022). Large-scale pre-stack seismic depth-domain inversion: a case study from the northern carnarvon basin. Offshore Australia: European Association of Geoscientists and Engineers, 1–5.
15
SacchiM. D.UlrychT. J. J. (2000). Nonminimum-phase wavelet estimation using higher order statistics. Lead. Edge19 (1), 80–83. 10.1190/1.1438466
16
SenguptaM.ZhangH.JervisM. (2018). Direct depth domain elastic inversion. OnePetro, 206–510.
17
SenguptaM.ZhangH.ZhaoY.JervisM.GranaD. (2021). Direct depth-domain Bayesian amplitude-variation-with-offset inversion. 86(5): M167–M176. 10.1190/geo2020-0219.1
18
ShadlowJ.ChristiansenD.Al-HouliM.PaxtonA.WilsonT. (2022). Getting the most out of a large data set: a case study for a large 3D seismic interpretation project in the Carnarvon Basin, NW Australia. Lead. Edge41 (12), 857–863. 10.1190/tle41120857.1
19
StockwellR. G.MansinhaL.LoweR. (1996). Localization of the complex spectrum: the S transform. IEEE Trans. signal Process.44 (4), 998–1001. 10.1109/78.492555
20
TangY.ChenS.WangG.LiX. (2024). Depth-domain direct inversion based on the generalized unscaled S-transform method. Oslo, Norway: European Association of Geoscientists and Engineers, 1–5.
21
Van Der BaanM. (2008). Time-varying wavelet estimation and deconvolution by kurtosis maximization. Geophysics73 (2), V11–V18. 10.1190/1.2831936
22
WaldenA. T.WhiteR. E. (1998). Seismic wavelet estimation: a frequency domain solution to a geophysical noisy input-output problem. IEEE Trans. Geoscience Remote Sens.36 (1), 287–297. 10.1109/36.655337
23
WangB. (2016). An amplitude preserving S-transform for seismic data attenuation compensation. IEEE Signal Process. Lett.23 (9), 1155–1159. 10.1109/lsp.2016.2586445
24
WangN.ShiY.NiJ.FangJ.YuB. (2024). Enhanced seismic attenuation compensation: integrating attention mechanisms with residual learning in neural networks. IEEE Trans. Geoscience Remote Sens.62, 1–11. 10.1109/tgrs.2024.3445130
25
WangN.XingG.ZhuT.ZhouH.ShiY. (2022a). Propagating seismic waves in VTI attenuating media using fractional viscoelastic wave equation. J. Geophys. Res. Solid Earth127 (4), e2021JB023280. 10.1029/2021jb023280
26
WangP.CuiY.-A.ChenX.PanX.LiuJ. (2022b). Analysis and application of the sparse prior in probabilistic prediction of elastic parameters. IEEE Trans. Geoscience Remote Sens.60, 1–9. 10.1109/tgrs.2022.3181175
27
WuL.CastagnaJ. (2017). S-transform and Fourier transform frequency spectra of broadband seismic signals. Geophysics82 (5), O71–O81. 10.1190/geo2016-0679.1
28
YilmazÖ. (2001). Seismic data analysis: processing, inversion, and interpretation of seismic data. Tulsa: Society of exploration geophysicists.
29
ZhangJ.ChenX.JiangW.LiuY.XuH. (2022a). Estimation of the depth-domain seismic wavelet based on velocity substitution and a generalized seismic wavelet model. Geophysics87 (2), R213–R222. 10.1190/geo2020-0745.1
30
ZhangJ.LiJ.ChenX.LiY. (2021a). Geological structure-guided hybrid MCMC and Bayesian linearized inversion methodology. J. Petroleum Sci. Eng.199, 108296. 10.1016/j.petrol.2020.108296
31
ZhangJ.LiJ.ChenX.LiY.HuangG.ChenY. (2021b). Robust deep learning seismic inversion with a priori initial model constraint. Geophys. J. Int.225 (3), 2001–2019. 10.1093/gji/ggab074
32
ZhangJ.LiJ.WangS.GengW.ChengW. (2022b). An improved unscale S transform in frequency domain. IEEE Geosci. Remote Sens. Lett.20, 1–5. 10.1109/lgrs.2022.3232593
33
ZhangJ.ZhaoX.LiJ.ChenX. (2022c). Structure-oriented prestack waveform inversion. Geophysics87 (3), M73–M85. 10.1190/geo2021-0452.1
34
ZhangR.DengZ. (2018). A depth variant seismic wavelet extraction method for inversion of poststack depth-domain seismic data. Geophysics83 (6), R569–R579. 10.1190/geo2017-0816.1
35
ZhangR.DengZ. (2023). Depth-domain angle and depth variant seismic wavelets extraction for prestack seismic inversion. Geophysics88 (1), R1–R10. 10.1190/geo2021-0647.1
36
ZhouL.LiuX.LiJ.LiaoJ. (2021). Robust AVO inversion for the fluid factor and shear modulus. Geophysics86 (4), R471–R483. 10.1190/geo2020-0234.1
Summary
Keywords
depth-domain seismic data, non-stationary signal, time–frequency decomposition, modified S-transform, seismic wavelet extraction
Citation
Xue Y, Zhang J, Zhang G and Zhao D (2024) Estimation of the depth-variant seismic wavelet based on the modified unscaled S-transform. Front. Earth Sci. 12:1480487. doi: 10.3389/feart.2024.1480487
Received
14 August 2024
Accepted
16 October 2024
Published
29 October 2024
Volume
12 - 2024
Edited by
Xingye Liu, Chengdu University of Technology, China
Reviewed by
Ning Wang, Northeast Petroleum University, China
Guangtan Huang, Chinese Academy of Sciences (CAS), China
Updates
Copyright
© 2024 Xue, Zhang, Zhang and Zhao.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Gan Zhang, 404085229@qq.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.