Frequency-domain elastic full-waveform multiscale inversion method based on dual-level parallelism

2015 ◽  
Vol 12 (4) ◽  
pp. 545-554 ◽  
Author(s):  
Yuan-Yuan Li ◽  
Zhen-Chun Li ◽  
Kai Zhang ◽  
Xuan Zhang
Geophysics ◽  
2006 ◽  
Vol 71 (6) ◽  
pp. R91-R100 ◽  
Author(s):  
Kun Xu ◽  
Stewart A. Greenhalgh ◽  
MiaoYue Wang

In this paper, we investigate several source-independent methods of nonlinear full-waveform inversion of multicomponent elastic-wave data. This includes iterative estimation of source signature (IES), standard trace normalization (STN), and average trace normalization (ATN) inversion methods. All are based on the finite-element method in the frequency domain. One synthetic elastic crosshole model is used to compare the recovered images with all these methods as well as the known source signature (KSS) inversion method. The numerical experiments show that the IES method is superior to both STN and ATN methods in two-component, elastic-wave inversion in the frequency domain when the source signature is unknown. The STN and ATN methods have limitations associated with near-zero amplitudes (or polarity reversals) in traces from one of the components, which destroy the energy balance in the normalized traces and cause a loss of frequency information. But the ATN method is somewhat superior to the STN method in suppressing random noise and improving stability, as the developed formulas and the numerical experiments show. We suggest the IES method as a practical procedure for multicomponent seismic inversion.


2021 ◽  
Vol 11 (14) ◽  
pp. 6330
Author(s):  
Kai Wang ◽  
Meiyan Guo ◽  
Qingxia Xiao ◽  
Chuanyi Ma ◽  
Lingli Zhang ◽  
...  

Ahead geological prospecting, which can estimate adverse geology ahead of the tunnel face, is necessary in the process of tunnel construction. Due to its long detection range and good recognition effect on the interface, the seismic method is widely used in tunnel ahead prospecting. However, the observation space in tunnels is quite narrow compared to ground seismic prospecting, which leads to some problems in the acquisition of wave velocity, including: the velocity of the direct wave is used to replace the wave velocity of the forward rock approximately; the arrival time information of seismic waves is the main factor in time-travel inversion or the tomography method, which is sufficient to provide a simple model rather than deal with complex geological conditions. In view of the above problems, the frequency domain full waveform inversion method in ground prospecting is introduced to tunnel seismic prospecting. In addition, the optimized difference format is given according to the particularity of the tunnel environment. In this method, the kinematics and dynamics of the seismic wavefield are fully used to obtain more accurate wave velocity results. Simultaneously, forward modeling and inversion simulations on tunnel samples with typical adverse geological bodies are given here, which verified the validity and reliability of the proposed method.


Geophysics ◽  
2011 ◽  
Vol 76 (4) ◽  
pp. R109-R124 ◽  
Author(s):  
Hafedh Ben-Hadj-Ali ◽  
Stéphane Operto ◽  
Jean Virieux

Three-dimensional full waveform inversion (FWI) still suffers from prohibitively high computational costs that arise because of the seismic modeling for multiple sources that is performed at each nonlinear iteration of FWI. Building supershots by assembling several sources allows mitigation of the number of simulations per FWI iteration, although it adds crosstalk artifacts because of interference between the individual sources of the supershots. These artifacts themselves can be reduced by encoding each individual source with a random phase shift during assembling of the sources. The source encoding method is applied to an efficient frequency-domain FWI, in which a limited number of discrete frequencies or coarsely sampled frequency groups are inverted successively following a multiscale approach. Random codes can be regenerated at each FWI iteration or for each frequency of a group during each FWI iteration, to favor the destructive summation of crosstalk artifacts over FWI iterations. Either a limited number of sources (partial assembling) or the total number of sources (full assembling) can be combined into supershots. Wide-aperture acquisition geometries such as land or marine node acquisitions are considered, to allow one to stack a large number of shots in the full computational domain and to test different partial assembling strategies involving sources that are close to or distant from each other. Two-dimensional case studies show that partial-source assembling of distant shots has a limited sensitivity to noise, for a computational saving that is roughly proportional to the number of shots assembled into the supershots. On the other hand, full assembling is more sensitive to noise, and it requires successive inversions of finely sampled frequency groups with a large number of FWI iterations. In contrast, refining the shot interval to improve the fold degrades the models when full assembling is applied to noisy data. Preliminary 3D application of the method leads to the same conclusions that 2D case studies do, with regard to the footprint of crosstalk noise in the imaging.


Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. R339-R347
Author(s):  
Ramzi Djebbi ◽  
Tariq Alkhalifah

Full-waveform inversion (FWI) using the scattering integral (SI) approach is an explicit formulation of the inversion optimization problem. The inversion procedure is straightforward, and the dependence of the data residuals on the model parameters is clear. However, the biggest limitation associated with this approach is the huge computational cost in conventional exploration seismology applications. Modeling from each of the source and receiver locations is required to compute the update at every iteration, and that is prohibitively expensive, especially for 3D problems. To deal with this issue, we have developed a hybrid implementation of frequency-domain FWI, in which forward modeling is combined with ray tracing to compute the update. We use the sensitivity kernels computed from dynamic ray tracing to build the gradient. The data residual is still computed using finite-difference wavefield modeling. With ray theory, the Green’s function can be approximated using a coarser grid compared to wave-equation modeling. Therefore, the memory requirements, as well as the computational cost, are reduced significantly. Considering that in transmission FWI long-to-intermediate wavelengths are updated during the early iterations, we obtain accurate inverted models. The inversion scheme captured the anomaly embedded in the homogeneous background medium. For more complex models, the hybrid inversion method helps in improving the initial model with little cost compared with conventional SI inversion approaches. The accuracy of the inversion results shows the effectiveness of the hybrid approach for 3D realistic problems.


Symmetry ◽  
2021 ◽  
Vol 13 (2) ◽  
pp. 260
Author(s):  
Meng Suo ◽  
Dong Zhang ◽  
Yan Yang

Inspired by the large number of applications for symmetric nonlinear equations, an improved full waveform inversion algorithm is proposed in this paper in order to quantitatively measure the bone density and realize the early diagnosis of osteoporosis. The isotropic elastic wave equation is used to simulate ultrasonic propagation between bone and soft tissue, and the Gauss–Newton algorithm based on symmetric nonlinear equations is applied to solve the optimal solution in the inversion. In addition, the authors use several strategies including the frequency-grid multiscale method, the envelope inversion and the new joint velocity–density inversion to improve the result of conventional full-waveform inversion method. The effects of various inversion settings are also tested to find a balanced way of keeping good accuracy and high computational efficiency. Numerical inversion experiments showed that the improved full waveform inversion (FWI) method proposed in this paper shows superior inversion results as it can detect small velocity–density changes in bones, and the relative error of the numerical model is within 10%. This method can also avoid interference from small amounts of noise and satisfy the high precision requirements for quantitative ultrasound measurements of bone.


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.


Sign in / Sign up

Export Citation Format

Share Document