Recovery of Low Frequency Components of Seismic Data via Sparsity Driven Inversion

Author(s):  
L. Han ◽  
L.G. Han ◽  
Z.F. Tong ◽  
Y.F. Ye
2019 ◽  
Vol 16 (6) ◽  
pp. 1017-1031 ◽  
Author(s):  
Yong Hu ◽  
Liguo Han ◽  
Rushan Wu ◽  
Yongzhong Xu

Abstract Full Waveform Inversion (FWI) is based on the least squares algorithm to minimize the difference between the synthetic and observed data, which is a promising technique for high-resolution velocity inversion. However, the FWI method is characterized by strong model dependence, because the ultra-low-frequency components in the field seismic data are usually not available. In this work, to reduce the model dependence of the FWI method, we introduce a Weighted Local Correlation-phase based FWI method (WLCFWI), which emphasizes the correlation phase between the synthetic and observed data in the time-frequency domain. The local correlation-phase misfit function combines the advantages of phase and normalized correlation function, and has an enormous potential for reducing the model dependence and improving FWI results. Besides, in the correlation-phase misfit function, the amplitude information is treated as a weighting factor, which emphasizes the phase similarity between synthetic and observed data. Numerical examples and the analysis of the misfit function show that the WLCFWI method has a strong ability to reduce model dependence, even if the seismic data are devoid of low-frequency components and contain strong Gaussian noise.


Geophysics ◽  
2019 ◽  
Vol 84 (6) ◽  
pp. R989-R1001 ◽  
Author(s):  
Oleg Ovcharenko ◽  
Vladimir Kazei ◽  
Mahesh Kalita ◽  
Daniel Peter ◽  
Tariq Alkhalifah

Low-frequency seismic data are crucial for convergence of full-waveform inversion (FWI) to reliable subsurface properties. However, it is challenging to acquire field data with an appropriate signal-to-noise ratio in the low-frequency part of the spectrum. We have extrapolated low-frequency data from the respective higher frequency components of the seismic wavefield by using deep learning. Through wavenumber analysis, we find that extrapolation per shot gather has broader applicability than per-trace extrapolation. We numerically simulate marine seismic surveys for random subsurface models and train a deep convolutional neural network to derive a mapping between high and low frequencies. The trained network is then tested on sections from the BP and SEAM Phase I benchmark models. Our results indicate that we are able to recover 0.25 Hz data from the 2 to 4.5 Hz frequencies. We also determine that the extrapolated data are accurate enough for FWI application.


2020 ◽  
Author(s):  
Hyunggu Jun ◽  
Hyeong-Tae Jou ◽  
Han-Joon Kim ◽  
Sang Hoon Lee

<p>Imaging the subsurface structure through seismic data needs various information and one of the most important information is the subsurface P-wave velocity. The P-wave velocity structure mainly influences on the location of the reflectors during the subsurface imaging, thus many algorithms has been developed to invert the accurate P-wave velocity such as conventional velocity analysis, traveltime tomography, migration velocity analysis (MVA) and full waveform inversion (FWI). Among those methods, conventional velocity analysis and MVA can be widely applied to the seismic data but generate the velocity with low resolution. On the other hands, the traveltime tomography and FWI can invert relatively accurate velocity structure, but they essentially need long offset seismic data containing sufficiently low frequency components. Recently, the stochastic method such as Markov chain Monte Carlo (McMC) inversion was applied to invert the accurate P-wave velocity with the seismic data without long offset or low frequency components. This method uses global optimization instead of local optimization and poststack seismic data instead of prestack seismic data. Therefore, it can avoid the problem of the local minima and limitation of the offset. However, the accuracy of the poststack seismic section directly affects the McMC inversion result. In this study, we tried to overcome the dependency of the McMC inversion on the poststack seismic section and iterative workflow was applied to the McMC inversion to invert the accurate P-wave velocity from the simple background velocity and inaccurate poststack seismic section. The numerical test showed that the suggested method could successfully invert the subsurface P-wave velocity.</p>


Geophysics ◽  
2017 ◽  
Vol 82 (3) ◽  
pp. M43-M53 ◽  
Author(s):  
Zhaoyun Zong ◽  
Kun Li ◽  
Xingyao Yin ◽  
Ming Zhu ◽  
Jiayuan Du ◽  
...  

Seismic amplitude variation with offset (AVO) inversion is well-known as a popular and pragmatic tool used for the prediction of elastic parameters in the geosciences. Low frequencies missing from conventional seismic data are conventionally recovered from other geophysical information, such as well-log data, for estimating the absolute rock properties, which results in biased inversion results in cases of complex heterogeneous geologic targets or plays with sparse well-log data, such as marine or deep stratum. Broadband seismic data bring new opportunities to estimate the low-frequency components of the elastic parameters without well-log data. We have developed a novel AVO inversion approach with the Bayesian inference for broadband seismic data. The low-frequency components of the elastic parameters are initially estimated with the proposed broadband AVO inversion approach with the Bayesian inference in the complex frequency domain because seismic inversion in the complex frequency domain is helpful to recover the long-wavelength structures of the elastic models. Gaussian and Cauchy probability distribution density functions are used for the likelihood function and the prior information of model parameters, respectively. The maximum a posteriori probability solution is resolved to estimate the low-frequency components of the elastic parameters in the complex frequency domain. Furthermore, with those low-frequency components as initial models and constraints, the conventional AVO inversion method with the Bayesian inference in the time domain is further implemented to estimate the final absolute elastic parameters. Synthetic and field data examples demonstrate that the proposed AVO inversion in the complex frequency domain is able to predict the low-frequency components of elastic parameters well, and that those low-frequency components set a good foundation for the final estimation of the absolute elastic parameters.


Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. R385-R400
Author(s):  
Luca Bianchin ◽  
Emanuele Forte ◽  
Michele Pipan

Low-frequency components of reflection seismic data are of paramount importance for acoustic impedance (AI) inversion, but they typically suffer from a poor signal-to-noise ratio. The estimation of the low frequencies of AI can benefit from the combination of a harmonic reconstruction method (based on autoregressive [AR] models) and a seismic-derived interval velocity field. We have developed the construction of a convex cost function that accounts for the velocity field, together with geologic a priori information on AI and its uncertainty, during the AR reconstruction of the low frequencies. The minimization of this function allows one to reconstruct sensible estimates of low-frequency components of the subsurface reflectivity, which lead to an estimation of AI model via a recursive formulation. In particular, the method is suited for an initial and computationally inexpensive assessment of the absolute value of AI even when no well-log data are available. We first tested the method on layered synthetic models, then we analyzed its applicability and limitations on a real marine seismic data set that included tomographic velocity information. Despite a strong trace-to-trace variability in the results, which could partially be mitigated by multitrace inversion, the method demonstrates its capability to highlight lateral variations of AI that cannot be detected when the low frequencies only come from well-log information.


Geophysics ◽  
2013 ◽  
Vol 78 (2) ◽  
pp. WA15-WA20 ◽  
Author(s):  
Lasse Amundsen ◽  
Hongbo Zhou

We evaluated a solution to seismic deghosting that deghosts the low-frequency components of the seismic pressure data. In an approximation that neglected the dependence on wavenumbers, the low-frequency deghosted pressure field was computed trace-by-trace as the sum of the pressure field and its scaled temporally integrated and temporally differentiated fields. We gave simple numerical examples that demonstrated the concept. The method was found to deghost data up to a frequency that is typically half of the second notch frequency. On the low-frequency side, the deghosting method was limited by the signal-to-noise ratio. The low-frequency deghosting technique can be appropriate to apply to the part of seismic data that have penetrated and reflected beneath complex and attenuating overburdens such as basalt, salt, and chalk.


2022 ◽  
Vol 41 (1) ◽  
pp. 34-39
Author(s):  
Vincent Durussel ◽  
Dongren Bai ◽  
Amin Baharvand Ahmadi ◽  
Scott Downie ◽  
Keith Millis

The depth of penetration and multidimensional characteristics of seismic waves make them an essential tool for subsurface exploration. However, their band-limited nature can make it difficult to integrate them with other types of ground measurements. Consequently, far offsets and very low-frequency components are key factors in maximizing the information jointly inverted from all recorded data. This explains why extending seismic bandwidth and available offsets has become a major industry focus. Although this requirement generally increases the complexity of acquisition and has an impact on its cost, improvements have been clearly and widely demonstrated on marine data. Onshore seismic data have generally followed the same trend but face different challenges, making it more difficult to maximize the benefits, especially for full-waveform inversion (FWI). This paper describes a new dense survey acquired in 2020 in the Permian Basin and aims to objectively assess the quality and benefits brought by a richer low end of the spectrum and far offsets. For this purpose, we considered several aspects, from acquisition design and field data to FWI imaging and quantitative interpretation.


Author(s):  
G. Y. Fan ◽  
J. M. Cowley

It is well known that the structure information on the specimen is not always faithfully transferred through the electron microscope. Firstly, the spatial frequency spectrum is modulated by the transfer function (TF) at the focal plane. Secondly, the spectrum suffers high frequency cut-off by the aperture (or effectively damping terms such as chromatic aberration). While these do not have essential effect on imaging crystal periodicity as long as the low order Bragg spots are inside the aperture, although the contrast may be reversed, they may change the appearance of images of amorphous materials completely. Because the spectrum of amorphous materials is continuous, modulation of it emphasizes some components while weakening others. Especially the cut-off of high frequency components, which contribute to amorphous image just as strongly as low frequency components can have a fundamental effect. This can be illustrated through computer simulation. Imaging of a whitenoise object with an electron microscope without TF limitation gives Fig. 1a, which is obtained by Fourier transformation of a constant amplitude combined with random phases generated by computer.


Author(s):  
В. М. Мойсишин ◽  
M. V. Lyskanych ◽  
R. A. Zhovniruk ◽  
Ye. P. Majkovych

The purpose of the proposed article is to establish the causes of oscillations of drilling tool and the basic laws of the distribution of the total energy of the process of changing the axial dynamic force over frequencies of spectrum. Variable factors during experiments on the classical plan were the rigidity of drilling tool and the hardness of the rock. According to the results of research, the main power of the process of change of axial dynamic force during drilling of three roller cone bits is in the frequency range 0-32 Hz in which three harmonic frequency components are allocated which correspond to the theoretical values of low-frequency and gear oscillations of the chisel and proper oscillations of the bit. The experimental values of frequencies of harmonic components of energy and normalized spectrum as well as the magnitude of the dispersion of the axial dynamic force and its normalized values at these frequencies are presented. It has been found that with decreasing rigidity of the drilling tool maximum energy of axial dynamic force moves from the low-frequency oscillation region to the tooth oscillation area, intensifying the process of rock destruction and, at the same time, protecting the tool from the harmful effects of the vibrations of the bit. Reducing the rigidity of the drilling tool protects the bit from the harmful effects of the vibrations generated by the stand. The energy reductions in these fluctuations range from 47 to 77%.


Sign in / Sign up

Export Citation Format

Share Document