An optimal frequency-domain finite-difference operator with a flexible stencil and its application in discontinuous-grid modeling

Geophysics ◽  
2021 ◽  
pp. 1-41
Author(s):  
Na Fan ◽  
Xiao-Bi Xie ◽  
Lian-Feng Zhao ◽  
Xin-Gong Tang ◽  
Zhen-Xing Yao
Keyword(s):  
Finite Difference ◽  
Frequency Domain ◽  
Waveform Inversion ◽  
Uniform Grid ◽  
Optimal Method ◽  
Velocity Models ◽  
Full Waveform ◽  
Rotationally Symmetric ◽  
Grid Points ◽  
Coarse To Fine

We develop an optimal method to determine expansion parameters for flexible stencils in 2D scalar wave finite-difference frequency-domain (FDFD) simulation. The proposed stencil only requires the involved grid points to be paired and rotationally symmetric around the central point. We apply this method to the transition zone in discontinuous-grid modeling, where the key issue is designing particular FDFD stencils to correctly propagate the wavefield passing through the discontinuous interface. The proposed method can work in FDFD discontinuous-grid with arbitrary integer coarse-to-fine gird spacing ratios. Numerical examples are presented to demonstrate how to apply this optimal method for the discontinuous-grid FDFD schemes with spacing ratios 3 and 5. The synthetic wavefields are highly consistent to those calculated using the conventional dense uniform grid, while the memory requirement and computational costs are greatly reduced. For velocity models with large contrasts, the proposed discontinuous-grid FDFD method can significantly improve the computational efficiency in forward modeling, imaging and full waveform inversion.

3D finite-difference frequency-domain modeling of visco-acoustic wave propagation using a massively parallel direct solver: A feasibility study

Geophysics ◽  
2007 ◽  
Vol 72 (5) ◽  
pp. SM195-SM211 ◽  
Author(s):  
Stéphane Operto ◽  
Jean Virieux ◽  
Patrick Amestoy ◽  
Jean-Yves L’Excellent ◽  
Luc Giraud ◽  
...  
Keyword(s):  
Wave Propagation ◽  
Finite Difference ◽  
Frequency Domain ◽  
Waveform Inversion ◽  
Staggered Grid ◽  
Full Waveform Inversion ◽  
Domain Modeling ◽  
Acoustic Wave Propagation ◽  
Full Waveform ◽  
Grid Points

We present a finite-difference frequency-domain method for 3D visco-acoustic wave propagation modeling. In the frequency domain, the underlying numerical problem is the resolution of a large sparse system of linear equations whose right-hand side term is the source. This system is solved with a massively parallel direct solver. We first present an optimal 3D finite-difference stencil for frequency-domain modeling. The method is based on a parsimonious staggered-grid method. Differential operators are discretized with second-order accurate staggered-grid stencils on different rotated coordinate systems to mitigate numerical anisotropy. An antilumped mass strategy is implemented to minimize numerical dispersion. The stencil incorporates 27 grid points and spans two grid intervals. Dispersion analysis showsthat four grid points per wavelength provide accurate simulations in the 3D domain. To assess the feasibility of the method for frequency-domain full-waveform inversion, we computed simulations in the 3D SEG/EAGE overthrust model for frequencies 5, 7, and [Formula: see text]. Results confirm the huge memory requirement of the factorization (several hundred Figabytes) but also the CPU efficiency of the resolution phase (few seconds per shot). Heuristic scalability analysis suggests that the memory complexity of the factorization is [Formula: see text] for a [Formula: see text] grid. Our method may provide a suitable tool to perform frequency-domain full-waveform inversion using a large distributed-memory platform. Further investigation is still necessary to assess more quantitatively the respective merits and drawbacks of time- and frequency-domain modeling of wave propagation to perform 3D full-waveform inversion.

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.

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.

A finite-difference iterative solver of the Helmholtz equation for frequency-domain seismic wave modeling and full waveform inversion

Geophysics ◽  
2020 ◽  
pp. 1-43
Author(s):  
Xingguo Huang ◽  
Stewart Greenhalgh
Keyword(s):  
Finite Difference ◽  
Helmholtz Equation ◽  
Frequency Domain ◽  
Iteration Method ◽  
Krylov Subspace ◽  
Waveform Inversion ◽  
Full Waveform Inversion ◽  
Iterative Solver ◽  
Subspace Iteration ◽  
Full Waveform

We present a finite difference iterative solver of the Helmholtz equation for seismic modeling and inversion in the frequency-domain. The iterative solver involves the shifted Laplacian operator and two-level pre-conditioners. It is based on the application of the pre-conditioners to the Krylov subspace stabilized biconjugate gradient method. A critical factor for the iterative solver is the introduction of a new pre-conditioner into the Krylov subspace iteration method to solve the linear system resulting from the discretization of the Helmholtz equation. This new pre-conditioner is based upon a reformulation of an integral equation-based convergent Born series for the Lippmann-Schwinger equation to an equivalent differential equation. We demonstrate that the proposed iterative solver combined with the novel pre-conditioner when incorporated with the finite difference method accelerates the convergence of the Krylov subspace iteration method for frequency-domain seismic wave modeling. A comparison of a direct solver, a one-level Krylov subspace iterative solver and the proposed two-level iterative solver verified the accuracy and accelerated convergence of the new scheme. Extensive tests in full waveform inversion demonstrate the solver applicability to full waveform inversion applications.

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