Layered low-frequency extrapolation with deep learning in full-waveform inversion

2021 ◽  
Author(s):  
Yao Liu ◽  
Baodi Liu ◽  
Jianping Huang ◽  
Jun Wang ◽  
Honglong Chen ◽  
...  
Geophysics ◽  
2020 ◽  
Vol 85 (6) ◽  
pp. A37-A43
Author(s):  
Jinwei Fang ◽  
Hui Zhou ◽  
Yunyue Elita Li ◽  
Qingchen Zhang ◽  
Lingqian Wang ◽  
...  

The lack of low-frequency signals in seismic data makes the full-waveform inversion (FWI) procedure easily fall into local minima leading to unreliable results. To reconstruct the missing low-frequency signals more accurately and effectively, we have developed a data-driven low-frequency recovery method based on deep learning from high-frequency signals. In our method, we develop the idea of using a basic data patch of seismic data to build a local data-driven mapping in low-frequency recovery. Energy balancing and data patches are used to prepare high- and low-frequency data for training a convolutional neural network (CNN) to establish the relationship between the high- and low-frequency data pairs. The trained CNN then can be used to predict low-frequency data from high-frequency data. Our CNN was trained on the Marmousi model and tested on the overthrust model, as well as field data. The synthetic experimental results reveal that the predicted low-frequency data match the true low-frequency data very well in the time and frequency domains, and the field results show the successfully extended low-frequency spectra. Furthermore, two FWI tests using the predicted data demonstrate that our approach can reliably recover the low-frequency data.


2021 ◽  
Author(s):  
Senlin Yang ◽  
Peng Jiang ◽  
Yuxiao Ren ◽  
Xinji Xu

<p>The seismic full waveform inversion (FWI), as one of important ways to obtain the seismic wave velocity, has made rapid development in the last decade. In response to problems of cycle-skipping artifacts, dependence on the initial model, and low-frequency information in FWI, researchers have made many improvements, such as multi-scale envelope inversion and low-frequency extension. Recently, deep learning has been also adopted seismic data processing and interpretation, because of its strong nonlinear mapping ability. However, these works depend on labels used for training heavily, especially for the velocity model in the inversion, which prevents them from real application. Referring to these studies, this work combines low-frequency extension commonly as well as multiscale inversion with deep learning, and proposes a multi-scale FWI gradient optimization method based on CNN. CNN we designed is trained to predict the inversion gradient corresponding to the low-frequency band data in FWI, so that multi-scale gradient optimization can be directly used in multi-scale inversion, expanding the low-frequency information in the actual data and reducing the calculation in FWI. With a specially designed dataset, CNN is trained to optimize the gradients computed from the high-frequency band data by predicting the gradients corresponding to the low-frequency band data and the gradients corresponding to the mid-frequency band data, respectively. The predicted gradients are used in different stages of the multi-scale inversion. The low-frequency gradients are used to invert the initial structural construction so as not to rely on a good initial model, and the high-frequency gradients are used to improve the accuracy of the inversion results. In this way, low-frequency expansion and multiscale inversion can be achieved simultaneously. Our method achieves good results on the initial model for a given uniform wave velocity, effectively alleviating the reliance on the initial model in FWI. This study provides a new idea of combining deep learning and full waveform inversion, which will be effectively used in seismic data processing.</p>


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 ◽  
2021 ◽  
pp. 1-54
Author(s):  
Milad Bader ◽  
Robert G. Clapp ◽  
Biondo Biondi

Low-frequency data below 5 Hz are essential to the convergence of full-waveform inversion towards a useful solution. They help build the velocity model low wavenumbers and reduce the risk of cycle-skipping. In marine environments, low-frequency data are characterized by a low signal-to-noise ratio and can lead to erroneous models when inverted, especially if the noise contains coherent components. Often field data are high-pass filtered before any processing step, sacrificing weak but essential signal for full-waveform inversion. We propose to denoise the low-frequency data using prediction-error filters that we estimate from a high-frequency component with a high signal-to-noise ratio. The constructed filter captures the multi-dimensional spectrum of the high-frequency signal. We expand the filter's axes in the time-space domain to compress its spectrum towards the low frequencies and wavenumbers. The expanded filter becomes a predictor of the target low-frequency signal, and we incorporate it in a minimization scheme to attenuate noise. To account for data non-stationarity while retaining the simplicity of stationary filters, we divide the data into non-overlapping patches and linearly interpolate stationary filters at each data sample. We apply our method to synthetic stationary and non-stationary data, and we show it improves the full-waveform inversion results initialized at 2.5 Hz using the Marmousi model. We also demonstrate that the denoising attenuates non-stationary shear energy recorded by the vertical component of ocean-bottom nodes.


2021 ◽  
Author(s):  
J. Pruessmann ◽  
G. Eisenberg Klein ◽  
E. Schuenemann ◽  
E. Verschuur ◽  
S. Qu

Geophysics ◽  
2020 ◽  
Vol 85 (6) ◽  
pp. R553-R563
Author(s):  
Sagar Singh ◽  
Ilya Tsvankin ◽  
Ehsan Zabihi Naeini

The nonlinearity of full-waveform inversion (FWI) and parameter trade-offs can prevent convergence toward the actual model, especially for elastic anisotropic media. The problems with parameter updating become particularly severe if ultra-low-frequency seismic data are unavailable, and the initial model is not sufficiently accurate. We introduce a robust way to constrain the inversion workflow using borehole information obtained from well logs. These constraints are included in the form of rock-physics relationships for different geologic facies (e.g., shale, sand, salt, and limestone). We develop a multiscale FWI algorithm for transversely isotropic media with a vertical symmetry axis (VTI media) that incorporates facies information through a regularization term in the objective function. That term is updated during the inversion by using the models obtained at the previous inversion stage. To account for lateral heterogeneity between sparse borehole locations, we use an image-guided smoothing algorithm. Numerical testing for structurally complex anisotropic media demonstrates that the facies-based constraints may ensure the convergence of the objective function towards the global minimum in the absence of ultra-low-frequency data and for simple (even 1D) initial models. We test the algorithm on clean data and on surface records contaminated by Gaussian noise. The algorithm also produces a high-resolution facies model, which should be instrumental in reservoir characterization.


Geophysics ◽  
2021 ◽  
Vol 86 (1) ◽  
pp. R31-R44
Author(s):  
Bin Liu ◽  
Senlin Yang ◽  
Yuxiao Ren ◽  
Xinji Xu ◽  
Peng Jiang ◽  
...  

Velocity model inversion is one of the most important tasks in seismic exploration. Full-waveform inversion (FWI) can obtain the highest resolution in traditional velocity inversion methods, but it heavily depends on initial models and is computationally expensive. In recent years, a large number of deep-learning (DL)-based velocity model inversion methods have been proposed. One critical component in those DL-based methods is a large training set containing different velocity models. We have developed a method to construct a realistic structural model for the DL network. Our compressional-wave velocity model building method for creating dense-layer/fault/salt body models can automatically construct a large number of models without much human effort, which is very meaningful for DL networks. Moreover, to improve the inversion result on these realistic structural models, instead of only using the common-shot gather, we also extract features from the common-receiver gather as well. Through a large number of realistic structural models, reasonable data acquisition methods, and appropriate network setups, a more generalized result can be obtained through our proposed inversion framework, which has been demonstrated to be effective on the independent testing data set. The results of dense-layer models, fault models, and salt body models that we compared and analyzed demonstrate the reliability of our method and also provide practical guidelines for choosing optimal inversion strategies in realistic situations.


Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. R409-R423
Author(s):  
Polina Zheglova ◽  
Alison Malcolm

Vector-acoustic full-waveform inversion (VAFWI) directly inverts vector-acoustic (VA) data, which consist of pressure and particle displacement components, at the cost of conventional acoustic full-waveform inversion (FWI). VA data contain information about the direction of arrival of the recorded seismic waves. In VAFWI, this directional information is taken into account by introducing an appropriate data weighting. With this weighting, in the geometry of a marine seismic experiment, the VAFWI adjoint calculation approximates inverse wavefield extrapolation, resulting in the natural separation of up- and downgoing recorded waves. If the free-surface effects are modeled during the inversion, the wave separation leads to (1) suppression of surface-related artifacts, (2) constructive interference of receiver ghosts with their primaries leading to preservation of the low-frequency content in the adjoint fields, and (3) compensation for insufficient spatial wavefield sampling on the receiver side. The horizontal displacement component helps interpolate the missing data. Synthetic examples demonstrate that for undersampled data, VAFWI consistently recovers the subsurface properties with higher resolution and fewer artifacts than conventional FWI.


Sign in / Sign up

Export Citation Format

Share Document