1887
Volume 72, Issue 2
  • E-ISSN: 1365-2478

Abstract

Abstract

In this paper, we present a frequency‐domain volume integral method to model the microseismic wavefield in heterogeneous anisotropic‐elastic media. The elastic wave equation is written as an integral equation of the Lippmann–Schwinger type, and the seismic source is represented as a general moment tensor. The actual medium is split into a background medium and a scattered medium. The background part of the displacement field is computed analytically, but the scattered part requires a numerical solution. The existing matrix‐based implementation of the integral equation is computationally inefficient to model the wavefield in three‐dimensional earth. An integral equation for the particle displacement is, hence, formulated in a matrix‐free manner through the application of the Fourier transform. The biconjugate gradient stabilized method is used to iteratively obtain the solution of this equation. The integral equation method is naturally target oriented, and it is not necessary to fully discretize the model. This is very helpful in the microseismic wavefield computation at receivers in the borehole in many cases; say, for example, we want to focus only on the fluid injection zone in the reservoir–overburden system and not on the whole subsurface region. Additionally, the integral equation system matrix has a low condition number. This provides us flexibility in the selection of the grid size, especially at low frequencies for given wave velocities. Considering all these factors, we apply the numerical scheme to three different models in order of increasing geological complexity. We obtain the elastic displacement fields corresponding to the different types of moment tensor sources, which prove the utility of this method in microseismic. The generated synthetic data are intended to be used in inversion for the microseismic source and model parameters.

Loading

Article metrics loading...

/content/journals/10.1111/1365-2478.13416
2024-01-30
2025-07-14
Loading full text...

Full text loading...

/deliver/fulltext/gpr/72/2/gpr13416.html?itemId=/content/journals/10.1111/1365-2478.13416&mimeType=html&fmt=ahah

References

  1. Aki, K. & Richards, P.G. (2002) Quantitative seismology. Mill Valley, CA: University Science Books.
    [Google Scholar]
  2. Alles, E.J. & van Dongen, K.W.A. (2011) Perfectly matched layers for frequency‐domain integral equation acoustic scattering problems. IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, 58(5), 1077–1086.
    [Google Scholar]
  3. Auld, B. (1990) Acoustic fields and waves in solids. Malabar, FL: Krieger.
    [Google Scholar]
  4. Barrett, R., Berry, M., Chan, T.F., Demmel, J., Donato, J.M., Dongarra, J., Eijkhout, V., Pozo, R., Romine, C. & Vorst, H.V. (1994) Templates for the solution of linear systems: building blocks for iterative methods. Philadelphia, PA: Society for Industrial and Applied Mathematics.
    [Google Scholar]
  5. Červený, V. & Ravindra, R. (1971) Theory of seismic head waves. Toronto: University of Toronto Press.
    [Google Scholar]
  6. Červený, V. & Hron, F. (1980) The ray series method and dynamical ray tracing system for three‐dimensional inhomogeneous media. Bulletin of the Seismological Society of America, 70, 47–77.
    [Google Scholar]
  7. Červený, V. (1985) Ray synthetic seismograms for complex two‐and three‐dimensional structures. Journal of Geophysics, 58, 2–26.
    [Google Scholar]
  8. Červený, V. (2001) Seismic ray theory. Cambridge, UK: Cambridge University Press.
    [Google Scholar]
  9. Gibson, R.L. & Ben‐Menahem, A. (1991) Elastic wave scattering by anisotropic obstacles: application to fractured volumes. Journal of Geophysical Research, 96(B12), 19905–19924.
    [Google Scholar]
  10. Graves, R.W. (1996) Simulating seismic wave propagation in 3D elastic media using staggered‐grid finite differences. Bulletin of the Seismological Society of America, 86(4), 1091–1106.
    [Google Scholar]
  11. Grechka, V. (2020) Moment tensors of double‐couple microseismic sources in anisotropic formations. Geophysics, 85(1), KS1–KS11.
    [Google Scholar]
  12. Huang, X., Jakobsen, M., Eikrem, K.S. & Nævdal, G. (2019) Target‐oriented inversion of time‐lapse seismic waveform data. Communication in Computational Physics, 28(1), 249–275.
    [Google Scholar]
  13. Hudson, J.A. & Heritage, J.R. (1982) The use of the Born approximation in seismic scattering problems. Geophysical Journal of the Royal Astronomical Society, 66, 221–240.
    [Google Scholar]
  14. Jakobsen, M. & Johansen, T.A. (2000) Anisotropic approximations for mudrocks: a seismic laboratory study. Geophysics, 65(6), 1711–1725.
    [Google Scholar]
  15. Jakobsen, M. & Ursin, B. (2015) Full waveform inversion in the frequency domain using direct iterative T‐matrix methods. Journal of Geophysics and Engineering, 12(3), 400–418.
    [Google Scholar]
  16. Jakobsen, M. & Tveit, S. (2018) Distorted Born iterative T‐matrix method for inversion of CSEM data in anisotropic media. Geophysical Journal International, 214, 1524–1537.
    [Google Scholar]
  17. Jakobsen, M. & Wu, R.S. (2018) Accelerating the T‐matrix approach to seismic full‐waveform inversion by domain decomposition. Geophysical Prospecting, 66, 1039–1059.
    [Google Scholar]
  18. Jakobsen, M., Psencík, I., Iversen, E. & Ursin, B. (2020) Transition operator approach to seismic full‐waveform inversion in arbitrary anisotropic elastic media. Communications in Computational Physics, 28(1), 297–327.
    [Google Scholar]
  19. Jost, M. & Herrmann, R. (1989) A student's guide to review of moment tensors. Seismological Research Letter, 60(2), 37–57.
    [Google Scholar]
  20. Lei, L., Tan, J., Zhang, D., Malkoti, A., Abakumov, I. & Xie, Y. (2021) FDwave3D: a MATLAB solver for the 3D anisotropic wave equation using the finite‐difference method. Computational Geosciences, 25, 1565–1578.
    [Google Scholar]
  21. Lecture notes on Fourier transform in N dimensions . Departments of Radiology and Medical Physics, University of Wisconsin‐Madison.
  22. Madariaga, R. (2015) Seismic source theory. Amsterdam: Elsevier B.V.
    [Google Scholar]
  23. Malovichko, M., Khokhlov, N., Yavich, Zhdanov, M,. (2017) Approximate solutions of acoustic 3D integral equation and their application to seismic modeling and full‐waveform inversion. Journal of Computational Physics, 346, 318–339.
    [Google Scholar]
  24. Malovichko, M., Khokhlov, N., Yavich, Zhdanov, M,. (2018) Acoustic 3D modeling by the method of integral equations. Computers and Geosciences, 111, 223–234.
    [Google Scholar]
  25. Miles, J.W. (1960) Scattering of elastic waves by small inhomogeneities. Geophysics, 25, 642–648.
    [Google Scholar]
  26. Monchiet, V. & Bonnet, G. (2012) A polarization‐based FFT iterative scheme for computing the effective properties of elastic composites with arbitrary contrast. International Journal for Numerical Methods in Engineering, 89(11), 1419–1436.
    [Google Scholar]
  27. Nocedal, J. & Wright, S.J. (2000) Numerical optimization. Berlin: Springer Science+Business Media.
    [Google Scholar]
  28. Osnabrugge, G., Leedumrongwatthanakun, S., Vellekoop, I.M. (2016) A convergent Born series for solving the inhomogeneous Helmholtz equation in arbitrarily large media. Journal of Computational Physics, 322, 113–124.
    [Google Scholar]
  29. Pujol, J. (2003) Elastic wave propagation and generation in seismology. Cambridge, UK: Cambridge University Press.
    [Google Scholar]
  30. Rice, J.R. (1998) Notes on elastodynamics, Green's function, and response to transformation strain and crack or fault sources. Course on Earth and Planetary Sciences 263. Cambridge, MA: Harvard University.
  31. Saad, Y. (2003) Iterative methods for sparse linear systems. Philadelphia, PA: Society for Industrial and Applied Mathematics.
    [Google Scholar]
  32. Saputera, D.H., Jakobsen, M., van Dongen, K.W.A., Jahani, N., Eikrem, K.S. & Alyaev, S., (2023) 3D induction log modelling with integral equation method and domain decomposition preconditioning. arXiv. [Preprint] available from arXiv:2306.17537.
  33. Sei, A. & Symes, W. (1995) Dispersion analysis of numerical wave propagation and its computational consequences. Journal of Scientific Computing, 10(1), 1–27.
    [Google Scholar]
  34. Shewchuk, J.R. (1994) An introduction to the Conjugate Gradient method without the agonizing pain. Pittsburgh, PA: School of Computer Science, Carnegie Mellon University.
    [Google Scholar]
  35. Shi, P., Angus, D., Nowacki, A., Yuan, S., Wang, Y. (2018) Microseismic full waveform modeling in anisotropic media with moment tensor implementation. Surveys in Geophysics, 39, 567–611.
    [Google Scholar]
  36. Snieder, R. (2002) General theory of elastic wave scattering. In: Scattering and Inverse scattering in pure and applied science. Amsterdam: Elsevier. pp. 528‐542.
    [Google Scholar]
  37. Stein, S. & Wysession, M. (2003) An introduction to seismology, earthquakes, and earth structure. Oxford, UK: Blackwell.
    [Google Scholar]
  38. Stierle, E., Vavryčuk, V., Kwiatek, G., Charalampidou, E.M., Bohnho, M. (2016) Seismic moment tensors of acoustic emissions recorded during laboratory rock deformation experiments: sensitivity to attenuation and anisotropy.. Geophysical Journal International, 205(1), 38–50.
    [Google Scholar]
  39. Tong, M.S. & Chew, W.C. (2009) Multilevel fast multipole algorithm for elastic wave scattering by large three‐dimensional objects. Journal of Computational Physics, 228(3), 921–932.
    [Google Scholar]
  40. Touhei, T. (2011) A fast volume integral equation method for elastic wave propagation in a half space. International Journal of Solids and Structures, 48, 22–23.
    [Google Scholar]
  41. Van der Vorst, H.A. (1992) Bi‐CGSTAB: a fast and smoothly converging variant of Bi‐CG for the solution of nonsymmetric linear systems. SIAM Journal on Scientific and Statistical Computing, 13(2), 631–644.
    [Google Scholar]
  42. Vavryčuk, V. (2005) Focal mechanisms in anisotropic media. Geophysical Journal International, 161(2), 334–346.
    [Google Scholar]
  43. Vavryčuk, V. (2015) Moment tensor decompositions revisited. Journal of Seismology, 19, 231–252.
    [Google Scholar]
  44. Virieux, J. (1986) P‐Sv wave propagation in heterogeneous media: velocity‐stress finite‐difference method. Geophysics, 51(4), 889–901.
    [Google Scholar]
  45. Worthington, M.H. (2008) Interpreting seismic anisotropy in fractured reservoirs. First Break, 26, 57–63.
    [Google Scholar]
  46. Wu, R.S. & Ben‐Menahem, A. (1985) The elastodynamic near field. Geophysical Journal of the Royal Astronomical Society, 81, 609–622.
    [Google Scholar]
  47. Xiang, K., Eikrem, K.S., Jakobsen, M. & Nævdal, G. (2021) Homotopy scattering series for seismic forward modelling with variable density and velocity. Geophysical Prospecting, 70(1), 3–18.
    [Google Scholar]
  48. Ying, L. (2015) Sparsifying preconditioner for the Lippmann‐Schwinger equation. SIAM Multiscale Model Simulation, 13(2), 644–660.
    [Google Scholar]
  49. Zhdanov, M.S. & Fang, S. (1997) Quasi‐linear series in three‐dimensional electromagnetic modeling. Radio Science, 32(6), 2167–2188.
    [Google Scholar]
/content/journals/10.1111/1365-2478.13416
Loading
/content/journals/10.1111/1365-2478.13416
Loading

Data & Media loading...

  • Article Type: Research Article
Keyword(s): anisotropy; elastics; full waveform; modelling; seismics; wave

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