scholarly journals Frequency Domain Wave Equation Inversion and Its Application on the Heterogeneous Reservoir Model Data

2016 ◽  
Vol 6 (1) ◽  
pp. 55
Author(s):  
Zhiyuan Tang
Keyword(s):  
Wave Equation ◽  
Velocity Field ◽  
Frequency Domain ◽  
Waveform Inversion ◽  
Gradient Methods ◽  
Carbonate Reservoir ◽  
Physical Models ◽  
Full Waveform Inversion ◽  
Matrix Formalism ◽  
Full Waveform

Seismic full waveform inversion seeks to make use of the full information based on full wave field modeling to extract quantitative information from seismograms. Its serious nonlinearity and high dependence on initial velocity model often results in unsatisfactory inversion results in paleo-karsts carbonate reservoir characterized by strong heterogeneity. The paper presents an improved strategy of multi-scale inversion to establish velocity field model of waveform tomography. the forward wave equation algorithm was derived in frequency domain, and then the Matrix formalism for the iterative inverse methods is derived by gradient methods to speed up calculation and to avoid convergence to local minimum value. After massive amount of frequencies tests, the appropriate bandwidth are extracted, and the velocity field calculated at low frequency is used as the input of the high frequency. After the iteration, the accurate velocity field is inverted. Finally, frequency domain wave equation full waveform inversion in mathematical and physical models is conducted in order to verify the inverse program. The method of selecting the inverse frequencies is proved to be effective.

On efficient frequency-domain full-waveform inversion with extended search space

Geophysics ◽  
2020 ◽  
pp. 1-61
Author(s):  
Hossein S. Aghamiry ◽  
Ali Gholami ◽  
Stéphane Operto
Keyword(s):  
Wave Equation ◽  
Frequency Domain ◽  
Computational Efficiency ◽  
Waveform Inversion ◽  
Search Space ◽  
Full Waveform Inversion ◽  
Lu Decomposition ◽  
Multiple Sources ◽  
Full Waveform ◽  
Convergence Point

Efficient frequency-domain Full-Waveform Inversion (FD-FWI) of wide-aperture data is designed by limiting inversion to few frequencies and by solving the Helmholtz equation with a direct solver to process multiple sources efficiently. Some variants of FD-FWI, which process the wave-equation as a weak constraint, were proposed to increase the computational efficiency or extend the search space. Among them, the contrast-source reconstruction inversion (CSRI) reparametrizes FD-FWI in terms of contrast sources (CS) and contrasts and update them in an alternating mode.This reparametrization allows for one lower-upper (LU) decomposition of the Helmholtz operator to be performed per frequency inversion hence further improving the computational efficiency of FD-FWI.On the other hand, Iteratively-refined Wavefield Reconstruction Inversion (IR-WRI) relies on the alternating-direction method of multipliers (ADMM) to extend the search space by matching the data from the early iterations via an aggressive relaxation of the wave-equation while satisfying it at the convergence point thanks to the defect correction performed by the Lagrange multipliers. In contrast to CSRI, IR-WRI requires to redo one LU decomposition when the subsurface model is updated.In both methods, the CSs or the wavefields are computed by solving in a least-squares sense an overdetermined linear system gathering an observation equation and a wave-equation.A drawback of CSRI is that CSs are estimated approximately with one iteration of a conjugate gradient method, while the wavefields are reconstructed exactly by IR-WRI with a Gauss-Newton method. We combine the benefits of CSRI and IR-WRI to decrease the number of LU decomposition during IR-WRI with a fixed-point (FP) algorithm while preserving its search space extension capability. Application on the 2D complex Marmousi and the BP salt models shows that our FP-based IR-WRI manages to reconstruct these models as accurately as the classical IR-WRI while reducing the number of LU factorizations considerably.

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

scholarly journals Velocity model building by 3D frequency-domain, full-waveform inversion of wide-aperture seismic data

Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE101-VE117 ◽  
Author(s):  
Hafedh Ben-Hadj-Ali ◽  
Stéphane Operto ◽  
Jean Virieux
Keyword(s):  
High Resolution ◽  
Frequency Domain ◽  
Distributed Memory ◽  
Waveform Inversion ◽  
Steepest Descent Method ◽  
Frequency Domain Method ◽  
Full Waveform Inversion ◽  
Direct Solver ◽  
Velocity Models ◽  
Full Waveform

We assessed 3D frequency-domain (FD) acoustic full-waveform inversion (FWI) data as a tool to develop high-resolution velocity models from low-frequency global-offset data. The inverse problem was posed as a classic least-squares optimization problem solved with a steepest-descent method. Inversion was applied to a few discrete frequencies, allowing management of a limited subset of the 3D data volume. The forward problem was solved with a finite-difference frequency-domain method based on a massively parallel direct solver, allowing efficient multiple-shot simulations. The inversion code was fully parallelized for distributed-memory platforms, taking advantage of a domain decomposition of the modeled wavefields performed by the direct solver. After validation on simple synthetic tests, FWI was applied to two targets (channel and thrust system) of the 3D SEG/EAGE overthrust model, corresponding to 3D domains of [Formula: see text] and [Formula: see text], respectively. The maximum inverted frequencies are 15 and [Formula: see text] for the two applications. A maximum of 30 dual-core biprocessor nodes with [Formula: see text] of shared memory per node were used for the second target. The main structures were imaged successfully at a resolution scale consistent with the inverted frequencies. Our study confirms the feasibility of 3D frequency-domain FWI of global-offset data on large distributed-memory platforms to develop high-resolution velocity models. These high-velocity models may provide accurate macromodels for wave-equation prestack depth migration.

Sign in / Sign up

Export Citation Format

Share Document