1887
Volume 63, Issue 5
  • E-ISSN: 1365-2478

Abstract

ABSTRACT

Prediction of elastic full wavefields is required for reverse time migration, full waveform inversion, borehole seismology, seismic modelling, etc. We propose a novel algorithm to solve the Navier wave equation, which is based on multi‐block methodology for high‐order finite‐difference schemes on curvilinear grids. In the current implementation, the blocks are subhorizontal layers. Smooth anisotropic heterogeneous media in each layer can have strong discontinuities at the interfaces. A curvilinear adaptive hexahedral grid in blocks is generated by mapping the original 3D physical domain onto a parametric cube with horizontal layers and interfaces. These interfaces correspond to the main curvilinear physical contrast interfaces of a subhorizontally layered formation. The top boundary of the parametric cube handles the land surface with smooth topography. Free‐surface and solid–solid transmission boundary conditions at interfaces are approximated with the second‐order accuracy. Smooth media in the layers are approximated up to sixth‐order spatial schemes. All expected properties of the developed algorithm are demonstrated in numerical tests using corresponding parallel message passing interface code.

Loading

Article metrics loading...

/content/journals/10.1111/1365-2478.12231
2015-03-12
2024-03-28
Loading full text...

Full text loading...

References

  1. AlpertB., GreengardL. and HagstromT.2000. Rapid evaluation of nonreflecting boundary kernels for time‐domain wave propagation. SIAM Journal on Numerical Analysis37, 1138–1164.
    [Google Scholar]
  2. AppelöD. and PeterssonN.A.2009. A stable finite difference method for the elastic wave equation on complex geometries with free surfaces. Communication in Computational Physics5(1), 84–107.
    [Google Scholar]
  3. ChararaM., VershininA., SabitovD. and PekarG.2011. SEM wave propagation in complex media with tetrahedral to hexahedral mesh. 73rd EAGE Conference & Exhibition, Vienna, Austria, Extended Abstracts.
  4. ChapmanC.2004. Fundamentals of Seismic Wave Propagation, pp. 608. Cambridge University Press.
    [Google Scholar]
  5. ClaytonR. and EngquistB.1977. Absorbing boundary conditions for acoustic and elastic wave equations. Bulletin of the Seismological Society of America6, 1529–1540.
    [Google Scholar]
  6. CourantR. and HilbertD.1989. Methods of Mathematical Physics: Partial Differential Equations Vol. 2, pp. 879. Wiley‐Interscience, New York.
    [Google Scholar]
  7. De la PuenteJ., KaeserM., DumbserM. and IgelH.2007. An arbitrary high‐order discontinuous Galerkin method for elastic waves on unstructured meshes–IV anisotropy. Geophysical Journal International169, 1210–1228.
    [Google Scholar]
  8. DovgilovichL.2013. High‐order FD method on curvilinear grids for anisotropic elastodynamic simulations. 75th EAGE Conference & Exhibition, Extended Abstracts.
  9. DovgilovichL.E. and SofronovI.L.2013. High‐order finite‐difference method of calculating wavefields in anisotropic media. Seismic Technologies2, 4–31 (in Russian).
    [Google Scholar]
  10. DovgilovichL.E. and SofronovI.L.2014a. Anisotropic elastodynamic simulations by FDM and SEM on curvilinear nonmatched grids. SEG Technical Program, Expanded Abstracts, pp. 3361–3365.
  11. DovgilovichL.E. and SofronovI.L.2014b. High‐accuracy finite‐difference schemes for solving elastodynamic problems in curvilinear coordinates within multiblock approach. Applied Numerical Mathematics.
    [Google Scholar]
  12. GroteM.J. and KellerJ.B.1995. Exact non‐reflecting boundary conditions for the time dependent wave equation. SIAM Journal of Applied Mathematics55, 280–297.
    [Google Scholar]
  13. GroteM.J. and KellerJ.B.2000. Exact nonreflecting boundary condition for elastic waves, SIAM Journal of Applied Mathematics60(3), 803–819.
    [Google Scholar]
  14. HagstromT. and HagstromG.2012. Grid stabilization of high‐order one‐sided differencing II: second‐order wave equations. Journal of Computational Physics231(23), 7907–7931.
    [Google Scholar]
  15. HairerE., NorsettS.P. and WannerG.1993. Solving Ordinary Differential Equations I: Nonstiff Problems. Springer Science & Business, Berlin.
    [Google Scholar]
  16. HirschC.2007. Numerical Computation of Internal and External Flows: The Fundamentals of Computational Fluid Dynamics Vol. 1, 2nd edn. Elsevier Science.
    [Google Scholar]
  17. KomatitschD. and VilotteJ.P.1998. The spectral element method: an efficient tool to simulate the seismic response of 2D and 3D geological structures. Bulletin of the Seismological Society of America88, 368–392.
    [Google Scholar]
  18. KosloffD. and CarcioneJ.M.2010. Two‐dimensional simulation of Rayleigh waves with staggered sine/cosine transforms and variable grid spacing. Geophysics75(4), T133–T140.
    [Google Scholar]
  19. KostinV.I., LisitsaV.V., ReshetovaG.V. and TcheverdaV.A.2013. Finite difference simulation of elastic waves propagation through 3D heterogeneous multiscale media based on locally refined grids. Numerical Analysis and Applications6(1), 40–48.
    [Google Scholar]
  20. KozdonJ.E., DunhamE.M. and NordströmJ.2013. Simulation of dynamic earthquake ruptures in complex geometries using high‐order finite difference methods. Journal of Scientific Computing55, 92–124.
    [Google Scholar]
  21. LambH.1904. On the propagation of tremors over the surface of an elastic solid. Philosophical Transactions of the Royal Society A203, l–42.
    [Google Scholar]
  22. LanH. and ZhangZ.2011. Three‐dimensional wave‐field simulation in heterogeneous transversely isotropic medium with irregular free surface. Bulletin of the Seismological Society of America101(3), 1354–1370.
    [Google Scholar]
  23. LeeS.H., JennyP. and TchelepiH.A.2002. A finite‐volume method with hexahedral multiblock grids for modeling flow in porous media. Computational Geosciences6, 353–379.
    [Google Scholar]
  24. MattssonK.2012. Summation by parts operators for finite difference approximations of second‐derivatives with variable coefficients. Journal of Scientific Computing51(3), 650–682.
    [Google Scholar]
  25. MattssonK. and AlmquistM.2013. A solution to the stability issues with block norm summation by parts operators. Journal of Computational Physics253, 418–442.
    [Google Scholar]
  26. PateraA.T.1984. A spectral element method for fluid dynamics: Laminar flow in a channel expansion. Journal of Computational Physics54, 468–488.
    [Google Scholar]
  27. SofronovI.L.1993. Conditions for complete transparency on the sphere for the three‐dimensional wave equation. Doklady Mathematics46(2), 397–401.
    [Google Scholar]
  28. SofronovI.L1998. Artificial boundary conditions of absolute transparency for two‐ and three‐dimensional external time‐dependent scattering problems. European Journal of Applied Mathematics9(6), 561–588.
    [Google Scholar]
  29. SofronovI.L.2009. Differential part of transparent boundary conditions for certain hyperbolic systems of second‐order equations. Doklady Mathematics79(3), 412–414.
    [Google Scholar]
  30. SofronovI.L. and ZaitsevN.A.2010. Numerical generation of transparent boundary conditions on the side surface of a vertical transverse isotropic layer. Journal of Computational and Applied Mathematics234(6), 1732–1738.
    [Google Scholar]
  31. SofronovI.L., ZaitsevN.A., DaryinA., VoskoboinikovaO. and DovgilovichL.2012. Multi‐block FD method for 3D geophysical simulation with explicit representation of sub‐horizontal interfaces. 74th EAGE Conference & Exhibition, Extended Abstracts.
  32. SofronovI.L., DovgilovichL. and KrasnovN.2014. Application of transparent boundary conditions to high‐order finite‐difference schemes for the wave equation in waveguides. Applied Numerical Mathematics.
    [Google Scholar]
  33. StrandB.1994. Summation by parts for finite difference approximations for d/dx. Journal of Computational Physics110, 47–67.
    [Google Scholar]
  34. TarrassI., GiraudL. and ThoreP.2011. New curvilinear scheme for elastic wave propagating in presence of curved topography. Geophysical Prospecting59, 889–906.
    [Google Scholar]
  35. TessmerE., KosloffD. and BehleA.1992. Elastic wave propagation simulation in the presence of surface topography. Geophysical Journal International108(2), 621–632.
    [Google Scholar]
  36. ThomsenL.1986. Weak elastic anisotropy. Geophysics51, 1954–1966.
    [Google Scholar]
  37. ThompsonJ.F. and WeatherillN.P.1999. Fundamental concepts and approaches. In: Handbook of Grid Generation (eds J.F.Thompson , B.K.Soni and N.P.Weatherill ), pp. I1–I30. CRC Press, New York.
    [Google Scholar]
  38. TrompJ., KomatitschD. and LiuQ.2008. Spectral‐element and adjoint methods in eismology. Communications in Computational Physics3, 1–32.
    [Google Scholar]
  39. WheelerM.F., ArbogastT., BryantS., EatonJ., LinQ. and PeszynskaM.1999. A parallel/multidomain approach for reservoir simulation. Proceedings of 15th SPE Symposium on Reservoir Simulation, Houston, TX, 1999, pp. 51–61.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/journals/10.1111/1365-2478.12231
Loading
/content/journals/10.1111/1365-2478.12231
Loading

Data & Media loading...

  • Article Type: Research Article
Keyword(s): Anisotropy; Elastodynamics; Full waveform; Numerical modelling

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