1887
24th International Geophysical Conference and Exhibition – Geophysics and Geology Together for Discovery
  • ISSN: 2202-0586
  • E-ISSN:

Abstract

The key computational kernels of most advanced 3D exploration seismic imaging and inversion algorithms involve calculating solutions of the 3D acoustic wave equation, most commonly with a finite-difference time-domain (FDTD) methodology. While well suited for regularly sampled rectilinear computational domains, FDTD methods seemingly have limited applicability in scenarios involving irregular 3D domain boundaries and mesh interiors best described by non-Cartesian geometry (e.g., surface topography). Using coordinate mappings and differential geometry, I specify a FDTD approach for generating numerical solutions to the acoustic wave equation that is applicable to generalized 3D coordinate systems and (hexahedral) structured meshes. I validate the method on different computational meshes and demonstrate the viability of the modelling approach for 3D non-Cartesian imaging and inversion scenarios.

Loading

Article metrics loading...

/content/journals/10.1071/ASEG2015ab023
2015-12-01
2026-01-17
Loading full text...

Full text loading...

References

  1. Appelo, D., and Petersson, N.A., 2009, A stable finite-difference method for the elastic wave equation on complex geometries with free surfaces: Communications in Computational Physics, 5, 84-107.
  2. Cockburn, B., Karniadakis, G., and Shu, C., 2000, Discontinuous Galerkin methods: Theory, computational and applications: Springer-Verlag.
  3. Fornberg, B., 1988, Generation of finite difference formulas on arbitrarily spaced grids: Mathematics of Computation, 51, 699-706.
  4. Guggenheimer, H., 1977, Differential geometry: Dover Publications, Inc., New York.
  5. Hestholm, S., 1999, Three-dimensional finite difference viscoelastic wave modelling including surface topography: Geophysical Journal International, 139, 852-878.
  6. Komatitsch, D., and Vilotte, J.-P., 1998, The spectral element method: An efficient tool to simulate the response of 2D and 3D geological structures: Bulletin of the Seismological Society of America, 88, 368-392.
  7. Marfurt, K., 1984, Accuracy of finite-difference and finite-element modelling of the scalar and elastic wave equations: Geophysics, 49, 533-549.
  8. Ohminato, T., and Chouet, B.A., 1997, A free-surface boundary condition for including 3D topography in the finite-difference method: Bulletin of the Seismological Society of America, 87, 494-515.
  9. Priolo, E., Carcione, J. and Seriani, G., 1994, Numerical simulation of interface waves by high-order spectral modeling techniques: Journal of the Acoustical Society of America, 95, 681 -693.
  10. Robertsson, J., Blanch, J.O., Nihei, K., and Troemp, J. (Editors), 2012, Numerical Modeling of Seismic Wave Propagation: Gridded Two-way Wave-equation Methods: SEG, Geophysics Reprint Series 28.
  11. Sava, P., and Fomel, S., 2005, Riemannian wavefield extrapolation: Geophysics, 70, no. 3, T45-T56.
  12. Shragge, J., 2008, Riemannian wavefield extrapolation: Non-orthogonal coordinate systems: Geophysics, 73, no. 2, T11-T21.
  13. Shragge, J., 2014, Reverse time migration from topography: Geophysics, 79, no. 4, T1 -T12.
  14. Synge, J.L., and Schild, A., 1978, Tensor Calculus: Dover Publications, Inc., New York.
/content/journals/10.1071/ASEG2015ab023
Loading
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