scholarly journals Frequency Domain Full Waveform Inversion Method of Acquiring Rock Wave Velocity in Front of Tunnels

2021 ◽  
Vol 11 (14) ◽  
pp. 6330
Author(s):  
Kai Wang ◽  
Meiyan Guo ◽  
Qingxia Xiao ◽  
Chuanyi Ma ◽  
Lingli Zhang ◽  
...  
Keyword(s):  
Frequency Domain ◽  
Wave Velocity ◽  
Waveform Inversion ◽  
Inversion Method ◽  
Full Waveform Inversion ◽  
Tunnel Face ◽  
Detection Range ◽  
Seismic Prospecting ◽  
Full Waveform ◽  
Geological Conditions

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.

2D frequency-domain elastic full-waveform inversion using the block-diagonal pseudo-Hessian approximation

Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. R247-R259 ◽  
Author(s):  
Yuwei Wang ◽  
Liangguo Dong ◽  
Yuzhu Liu ◽  
Jizhong Yang
Keyword(s):  
Frequency Domain ◽  
Wave Velocity ◽  
Hessian Matrix ◽  
Waveform Inversion ◽  
Full Waveform Inversion ◽  
S Wave ◽  
Trade Off ◽  
Full Waveform ◽  
Trade Offs ◽  
Block Diagonal

Elastic full-waveform inversion (EFWI) of multicomponent seismic data is a powerful tool for estimating the subsurface elastic parameters with high accuracy. However, the trade-offs between multiple parameters increase the nonlinearity of EFWI. Although the conventional diagonal-approximate Hessian matrix describes the illumination and limited bandwidth effects, it ignores the trade-off effects and decreases the convergence rate of EFWI. We have developed a block-diagonal pseudo-Hessian operator for 2D frequency-domain EFWI to take into account the approximate trade-offs among the P-wave (compressional-wave) velocity, S-wave (shear-wave) velocity, and density without extra computational costs on forward simulations. The Hessian matrix tends toward a block-diagonal matrix as the frequency grows to infinity; thus, the proposed block-diagonal pseudo-Hessian matrix is more accurate at higher frequencies. The inverse of the block-diagonal pseudo-Hessian matrix is used as a preconditioner for the nonlinear conjugate-gradient method to simultaneously reconstruct P- and S-wave velocities and density. This approach effectively mitigates the crosstalk artifacts by correcting the gradients from the trade-off effects and produces more rapid inversion convergence, which becomes more significant at higher frequencies. Synthetic experiments on an inclusion model and the elastic Marmousi2 model demonstrate its feasibility and validity in EFWI.

An efficient frequency-domain full waveform inversion method using simultaneous encoded sources

Geophysics ◽  
2011 ◽  
Vol 76 (4) ◽  
pp. R109-R124 ◽  
Author(s):  
Hafedh Ben-Hadj-Ali ◽  
Stéphane Operto ◽  
Jean Virieux
Keyword(s):  
Frequency Domain ◽  
Case Studies ◽  
Waveform Inversion ◽  
Random Phase ◽  
Inversion Method ◽  
Multiscale Approach ◽  
Full Waveform Inversion ◽  
Computational Domain ◽  
Multiple Sources ◽  
Full Waveform

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.

Seismic imaging of complex onshore structures by 2D elastic frequency-domain full-waveform inversion

Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. WCC105-WCC118 ◽  
Author(s):  
Romain Brossier ◽  
Stéphane Operto ◽  
Jean Virieux
Keyword(s):  
Free Surface ◽  
Surface Waves ◽  
Frequency Domain ◽  
Wave Velocity ◽  
Waveform Inversion ◽  
Surface Effects ◽  
Full Waveform Inversion ◽  
S Wave ◽  
Velocity Models ◽  
Full Waveform

Quantitative imaging of the elastic properties of the subsurface at depth is essential for civil engineering applications and oil- and gas-reservoir characterization. A realistic synthetic example provides for an assessment of the potential and limits of 2D elastic full-waveform inversion (FWI) of wide-aperture seismic data for recovering high-resolution P- and S-wave velocity models of complex onshore structures. FWI of land data is challenging because of the increased nonlinearity introduced by free-surface effects such as the propagation of surface waves in the heterogeneous near-surface. Moreover, the short wavelengths of the shear wavefield require an accurate S-wave velocity starting model if low frequencies are unavailable in the data. We evaluated different multiscale strategies with the aim of mitigating the nonlinearities. Massively parallel full-waveform inversion was implemented in the frequency domain. The numerical optimization relies on a limited-memory quasi-Newton algorithm thatoutperforms the more classic preconditioned conjugate-gradient algorithm. The forward problem is based upon a discontinuous Galerkin (DG) method on triangular mesh, which allows accurate modeling of free-surface effects. Sequential inversions of increasing frequencies define the most natural level of hierarchy in multiscale imaging. In the case of land data involving surface waves, the regularization introduced by hierarchical frequency inversions is not enough for adequate convergence of the inversion. A second level of hierarchy implemented with complex-valued frequencies is necessary and provides convergence of the inversion toward acceptable P- and S-wave velocity models. Among the possible strategies for sampling frequencies in the inversion, successive inversions of slightly overlapping frequency groups is the most reliable when compared to the more standard sequential inversion of single frequencies. This suggests that simultaneous inversion of multiple frequencies is critical when considering complex wave phenomena.

scholarly journals Hybrid frequency-domain full-waveform inversion using ray+Born sensitivity kernels

Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. R339-R347
Author(s):  
Ramzi Djebbi ◽  
Tariq Alkhalifah
Keyword(s):  
Frequency Domain ◽  
Ray Tracing ◽  
Waveform Inversion ◽  
Hybrid Approach ◽  
Computational Cost ◽  
Equation Modeling ◽  
Inversion Method ◽  
Model Parameters ◽  
Full Waveform Inversion ◽  
Full Waveform

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.

scholarly journals Application of an Improved Ultrasound Full-Waveform Inversion in Bone Quantitative Measurement

Symmetry ◽  
2021 ◽  
Vol 13 (2) ◽  
pp. 260
Author(s):  
Meng Suo ◽  
Dong Zhang ◽  
Yan Yang
Keyword(s):  
Nonlinear Equations ◽  
Waveform Inversion ◽  
Optimal Solution ◽  
Inversion Method ◽  
Full Waveform Inversion ◽  
Inversion Algorithm ◽  
Elastic Wave Equation ◽  
Full Waveform ◽  
Ultrasonic Propagation ◽  
Joint Velocity

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.

scholarly journals Multi-scale time-frequency domain full waveform inversion with a weighted local correlation-phase misfit function

2019 ◽  
Vol 16 (6) ◽  
pp. 1017-1031 ◽  
Author(s):  
Yong Hu ◽  
Liguo Han ◽  
Rushan Wu ◽  
Yongzhong Xu
Keyword(s):  
Frequency Domain ◽  
Seismic Data ◽  
Waveform Inversion ◽  
Low Frequency ◽  
Full Waveform Inversion ◽  
Local Correlation ◽  
Time Frequency ◽  
Model Dependence ◽  
Full Waveform ◽  
Frequency Components

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.

Application of the variable projection scheme for frequency-domain full-waveform inversion

Geophysics ◽  
2013 ◽  
Vol 78 (6) ◽  
pp. R249-R257 ◽  
Author(s):  
Maokun Li ◽  
James Rickett ◽  
Aria Abubakar
Keyword(s):  
Frequency Domain ◽  
Waveform Inversion ◽  
Simulated Data ◽  
Calibration Procedure ◽  
Full Waveform Inversion ◽  
Inversion Algorithm ◽  
Unknown Parameters ◽  
Variable Projection ◽  
Full Waveform ◽  
Data Calibration

We found a data calibration scheme for frequency-domain full-waveform inversion (FWI). The scheme is based on the variable projection technique. With this scheme, the FWI algorithm can incorporate the data calibration procedure into the inversion process without introducing additional unknown parameters. The calibration variable for each frequency is computed using a minimum norm solution between the measured and simulated data. This process is directly included in the data misfit cost function. Therefore, the inversion algorithm becomes source independent. Moreover, because all the data points are considered in the calibration process, this scheme increases the robustness of the algorithm. Numerical tests determined that the FWI algorithm can reconstruct velocity distributions accurately without the source waveform information.

Acoustic multiscale full waveform inversion of frequency domain based on FCG

2015 ◽  
Author(s):  
Changlu Sun* ◽  
Guangzhi Zhang ◽  
Xinpeng Pan ◽  
Xingyao Yin
Keyword(s):  
Frequency Domain ◽  
Waveform Inversion ◽  
Full Waveform Inversion ◽  
Full Waveform
Sign in / Sign up

Export Citation Format

Share Document