1887
Volume 53, Issue 2
  • ISSN: 0812-3985
  • E-ISSN: 1834-7533

Abstract

As the key point of full-waveform inversion (FWI), numerical solution of wave equation is required to simulate the forward propagated and adjoint wavefield. When traditional regular grids are used for forward modelling, scattering artifacts may occur due to the stepped approximation of layer interfaces and rugged topography. Unstructured mesh or irregular grids method can achieve certain geometric flexibility, however, its algorithm is quite complex. Mesh-free FWI can effectively reduce the scattered artifacts under regular grids and avoid the extra computation in the process of irregular grids generation. For the implementation of mesh-free FWI method, an algorithm with fast generation of node distributions is used to discretize the velocity model, radial-basis function generated finite difference is used to realise seismic wave propagation numerical simulation, and Limited-memory BFGS (L-BFGS) algorithm is used for iteration. The mesh-free FWI method we proposed achieves flexibility of simulation region and abundant wavefield information. It reduces the storage required for FWI and illustrates the applicability of high-precision reconstruction of underground velocity in the case of rugged topography, which can provide more accurate velocity information for oil and gas exploration under complex geological conditions.

Loading

Article metrics loading...

/content/journals/10.1080/08123985.2021.1906643
2022-03-04
2026-01-22
Loading full text...

Full text loading...

References

  1. Bollig, E.F., N.Flyer, and G.Erlebacher. 2012. Solution to PDEs using radial basis function finite-differences (RBF-FD) on multiple GPUs. Journal of Computational Physics231: 7133–51.
    [Google Scholar]
  2. Daniel, A., and N.A.Petersson. 2009. A stable finite difference method for the elastic wave equation on complex geometries with free surfaces. Communications in Computational Physics5: 84–107.
    [Google Scholar]
  3. Fasshauer, G.E.2007. Meshfree approximation methods with Matlab. Singapore: World Scientific.
  4. Fornberg, B., and N.Flyer. 2015. Fast generation of 2-D node distributions for mesh-free PDE discretizations. Computers and Mathematics with Applications69: 531–44.
    [Google Scholar]
  5. Lailly, P., and J.B.Bednar. 1983. The seismic inverse problem as a sequence of before stack migrations. In Conference on inverse scattering: Theory and application, 206–20. Philadelphia, PA: SIAM.
  6. Li, S.Z., C.Y.Sun, and L.Tian. 2018. Unsplit complex frequency shifted PML implementation in mess-free discretization. 80th EAGE conference and exhibition.
  7. Lines, L.R., R.Slawinski, and R.P.Bording. 1999. A recipe for stability of finite-difference wave-equation computations. Geophysics64: 967–9.
    [Google Scholar]
  8. Liu, D.C., and J.Nocedal. 1989. On the limited memory BFGS method for large scale optimization. Mathematical Programming45: 503–28.
    [Google Scholar]
  9. Martin, B., B.Fornberg, and A.St-Cyr. 2015. Seismic modeling with radial-basis-function-generated finite differences. Geophysics80: T137–T146.
    [Google Scholar]
  10. Nocedal, J., and S.J.Wright. 2006. Numerical optimization. 2nd ed. Berlin: Springer.
  11. Plessix, R.E.2006. A review of the adjoint-state method for computing the gradient of a functional with geophysical applications. Geophysical Journal International167: 495–503.
    [Google Scholar]
  12. Qu, Y.M., J.P.Huang, and Z.C.Li. 2015. Undulating free surface elastic wave prestack inverse time migration by hierarchical coordinate transformation method. Chinese Journal of Geophysics (in Chinese)58: 2896–911.
    [Google Scholar]
  13. Qu, Y., Z.Li, J.Huang, and J.Li. 2018. Multi-scale full waveform inversion for areas with irregular surface topography in an auxiliary coordinate system. Exploration Geophysics49: 68–80.
    [Google Scholar]
  14. Sun, C.Y., S.Z.Li, and N.Xu. 2019. PML and CFS-PML boundary conditions for a mesh-free finite difference solution of the elastic wave equation. Applied Geophysics16: 438–54.
    [Google Scholar]
  15. Takekawa, J., and H.Mikada. 2016. An absorbing boundary condition for acoustic-wave propagation using a mesh-free method. Geophysics81: T145–T154.
    [Google Scholar]
  16. Takekawa, J., H.Mikada, and N.Imamura. 2015. A mesh-free method with arbitrary-order accuracy for acoustic wave propagation. Computers and Geosciences78: 15–25.
    [Google Scholar]
  17. Takekawa, J., and H.A.Mikada. 2018. A mesh-free finite-difference method for elastic wave propagation in the frequency-domain. Computers and Geosciences118: 65–78.
    [Google Scholar]
  18. Tarantola, A.1984. Inversion of seismic reflection data in the acoustic approximation. Geophysics49: 1259–66.
    [Google Scholar]
  19. Thompson, J.F., Z.U.A.Warsi, and C.W.Mastin. 1985. Numerical grid generation: Foundations and applications. Amsterdam: North-Holland.
  20. Vigh, D., E.Starr, and J.Kapoor. 2009. Developing Earth models with full waveform inversion. The Leading Edge28: 432–5.
    [Google Scholar]
  21. Virieux, J., and S.Operto. 2009. An overview of full-waveform inversion in exploration geophysics. Geophysics74: WCC1–WCC26.
    [Google Scholar]
  22. Wang, Y., and L.G.Dong. 2014. Selection strategy of the initial matrix for L-BFGS method in time domain full waveform inversion. Geophysical Prospecting for Petroleum53: 545–55.
    [Google Scholar]
/content/journals/10.1080/08123985.2021.1906643
Loading
/content/journals/10.1080/08123985.2021.1906643
Loading

Data & Media loading...

  • Article Type: Research Article
Keyword(s): finite difference; Full waveform; inversion; time-domain; wave propagation

Most Cited This Month Most Cited RSS feed

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error