Abstract
We interpret uncertainty in a model for seismic wave propagation by treating the model parameters as random variables, and apply the Multilevel Monte Carlo method to reduce the cost of approximating expected values of selected, physically relevant, quantities of interest (QoI) with respect to the random variables. Targeting source inversion problems, where the source of an earthquake is inferred from ground motion recordings on the Earth’s surface, we consider two QoIs that measure the discrepancies between computed seismic signals and given reference signals: one QoI, \({\mathscr {Q}}_E\), is defined in terms of the \(L^2\)-misfit, which is directly related to maximum likelihood estimates of the source parameters; the other, \({\mathscr {Q}}_W\), is based on the quadratic Wasserstein distance between probability distributions, and represents one possible choice in a class of such misfit functions that have become increasingly popular to solve seismic inversion in recent years. We simulate seismic wave propagation, including seismic attenuation, using a publicly available code in widespread use, based on the spectral element method. Using random coefficients and deterministic initial and boundary data, we present benchmark numerical experiments with synthetic data in a two-dimensional physical domain and a one-dimensional velocity model where the assumed parameter uncertainty is motivated by realistic Earth models. Here, the computational cost of the standard Monte Carlo method was reduced by up to \(97\%\) for \({\mathscr {Q}}_E\), and up to \(78\%\) for \({\mathscr {Q}}_W\), using a relevant range of tolerances. Shifting to three-dimensional domains is straight-forward and will further increase the relative computational work reduction.
Similar content being viewed by others
References
Aki, K., Richards, P.G.: Quantitative Seismology: Theory and Methods, 2nd edn. University Science Books, Sausalito (2002)
Albaric, J., Perrot, J., Déverchère, J., Deschamps, A., Gall, B.L., Ferdinand, R., Petit, C., Tiberi, C., Sue, C., Songo, M.: Contrasted seismogenic and rheological behaviours from shallow and deep earthquake sequences in the North Tanzanian Divergence. East Africa. J. Afr. Earth Sci. 58(5), 799–811 (2010). https://doi.org/10.1016/j.jafrearsci.2009.09.005. (Active Volcanism and Continental Rifting in Africa)
Berenger, J.P.: A perfectly matched layer for the absorption of electromagnetic waves. J. Comput. Phys. 114(2), 185–200 (1994). https://doi.org/10.1006/jcph.1994.1159
Bissiri, P.G., Holmes, C.C., Walker, S.G.: A general framework for updating belief distributions. J. R. Stat. Soc. Ser. B Stat. Methodol. 78(5), 1103–1130 (2016). https://doi.org/10.1111/rssb.12158
Blanc, E., Komatitsch, D., Chaljub, E., Lombard, B., Xie, Z.: Highly accurate stability-preserving optimization of the zener viscoelastic model, with application to wave propagation in the presence of strong attenuation. Geophys. J. Int. 205(1), 427–439 (2016). https://doi.org/10.1093/gji/ggw024
Brown, T.S., Du, S., Eruslu, H., Sayas, F.-J.: Analysis for models of viscoelastic wave propagation. Appl. Math. Nonlinear Sci. 3(1), 55–96 (2018). https://doi.org/10.21042/AMNS.2018.1.00006
Carcione, J.M.: Wave Fields in Real Media; Wave Propagation in Anisotropic, Anelastic, Porous and Electromagnetic Media, 3rd edn. Elsevier Science, Amsterdam (2014)
Carcione, J.M., Kosloff, D., Kosloff, R.: Wave propagation simulation in a linear viscoelastic medium. Geophys. J. 95(3), 597–611 (1988). https://doi.org/10.1111/j.1365-246X.1988.tb06706.x
Collier, N., Haji-Ali, A.L., Nobile, F., von Schwerin, E., Tempone, R.: A continuation multilevel Monte Carlo algorithm. BIT Numer. Math. 55(2), 399–432 (2015). https://doi.org/10.1007/s10543-014-0511-3
Cruz-Jiménez, H., Li, G., Mai, P.M., Hoteit, I., Knio, O.M.: Bayesian inference of earthquake rupture models using polynomial chaos expansion. Geosci. Model Dev. 11(7), 3071–3088 (2018). https://doi.org/10.5194/gmd-11-3071-2018
Dahlen, F.A., Tromp, J.: Theoretical Global Seismology, 2nd edn. Princeton University Press, Princeton (1998)
Emmerich, H., Korn, M.: Incorporation of attenuation into time-domain computations of seismic wave fields. Geophysics 52(9), 1252–1264 (1987). https://doi.org/10.1190/1.1442386
Engquist, B., Hamfeldt, B.: Application of the Wasserstein metric to seismic signals. Commun. Math. Sci. 12(5), 979–988 (2014). https://doi.org/10.4310/CMS.2014.v12.n5.a7
Engquist, B., Froese, B.D., Yang, Y.: Optimal transport for seismic full waveform inversion. Commun. Math. Sci. 14(8), 2309–2330 (2016). https://doi.org/10.4310/CMS.2016.v14.n8.a9
Giles, M.B.: Multilevel Monte Carlo path simulation. Oper. Res. 56(3), 607–617 (2008). https://doi.org/10.1287/opre.1070.0496
Giles, M.B.: Multilevel Monte Carlo methods. Acta Numer. 24, 259–328 (2015). https://doi.org/10.1017/S096249291500001X
Haji-Ali, A.L., Nobile, F., von Schwerin, E., Tempone, R.: Optimization of mesh hierarchies in multilevel Monte Carlo samplers. Stoch. Partial Differ. Equ. Anal. Comput. 4(1), 76–112 (2016). https://doi.org/10.1007/s40072-015-0049-7
Heinrich, S.: Multilevel Monte Carlo Methods. In: Margenov, S., WasÌniewski, J., Yalamov, P. (eds.) Large-Scale Scientific Computing: Third International Conference, LSSC 2001, Sozopol, Bulgaria, June 6–10, 2001. Lecture Notes in Computer Science, vol. 2179, pp. 58–67. Springer (2001). https://doi.org/10.1007/3-540-45346-6_5
Hoel, H., Krumscheid, S.: Central limit theorems for multilevel Monte Carlo methods (2018). J Complex. (2019). https://doi.org/10.1016/j.jco.2019.05.001
Hughes, T.J.R.: The Finite Element Method: Linear Static and Dynamic Finite Element Analysis. Dover Civil and Mechanical Engineering, 2nd edn. Dover Publications, Mineola (2000)
Komatitsch, D., Tromp, J.: Introduction to the spectral element method for three-dimensional seismic wave propagation. Geophys. J. Int. 139(3), 806–822 (1999). https://doi.org/10.1046/j.1365-246x.1999.00967.x
Komatitsch, D., Tromp, J.: A perfectly matched layer absorbing boundary condition for the second-order seismic wave equation. Geophys. J. Int. 154(1), 146–153 (2003). https://doi.org/10.1046/j.1365-246X.2003.01950.x
Komatitsch, D., Vilotte, J.P.: The spectral element method: an efficient tool to simulate the seismic response of 2d and 3d geological structures. Bull. Seismol. Soc. Am. 88(2), 368–392 (1998). https://doi.org/10.1190/1.1820185
Komatitsch, D., Vilotte, J.P., Cristini, P., Labarta, J., Le Goff, N., Le Loher, P., Liu, Q., Martin, R., Matzen, R., Morency, C., Peter, D., Tape, C., Tromp, J., Xie, Z.: Specfem2d v7.0.0 [software] (2012). GitHub repository. https://github.com/geodynamics/specfem2d. Accessed 25 Oct 2017
Kortas, S.: decimate v0.9.5 [software] (2018). GitHub repository. https://github.com/KAUST-KSL/decimate. Accessed 2 Oct 2018
Lions, J.L., Magenes, E.: Non-homogeneous Boundary Value Problems and Applications, vol. I. Springer, Berlin (1972a)
Lions, J.L., Magenes, E.: Non-homogeneous Boundary Value Problems and Applications, vol. II. Springer, Berlin (1972b)
Liu, H.P., Anderson, D.L., Kanamori, H.: Velocity dispersion due to anelasticity: implications for seismology and mantle composition. Geophys. J. R. Astron. Soc. 47, 41–58 (1976). https://doi.org/10.1111/j.1365-246x.1976.tb01261.x
Liu, J., Chauris, H., Calandra, H.: The normalized integration method—an alternative to full waveform inversion? In: 25th Symposium on the Application of Geophpysics to Engineering and Environmental Problems, Best of 2011 EAGE/NSGD. United States (2012). https://doi.org/10.4133/1.4721678
Moczo, P., Robertsson, J.O., Eisner, L.: The finite-difference time-domain method for modeling of seismic wave propagation. Adv. Geophys. 48, 421–516 (2007). https://doi.org/10.1016/S0065-2687(06)48008-0
Motamed, M., Appelö, D.: Wasserstein metric-driven Bayesian inversion with application to wave propagation problems (2018). arXiv:1807.09682 [math.NA]
Motamed, M., Nobile, F., Tempone, R.: A stochastic collocation method for the second order wave equation with a discontinuous random speed. Numer. Mathe. 123(3), 493–536 (2013). https://doi.org/10.1007/s00211-012-0493-5
Métivier, L., Allain, A., Brossier, R., Mérigot, Q., Oudet, E., Virieux, J.: Optimal transport for mitigating cycle skipping in full-waveform inversion: a graph-space transform approach. Geophysics 83(5), R515–R540 (2018). https://doi.org/10.1190/geo2017-0807.1
Roecker, S., Ebinger, C., Tiberi, C., Mulibo, G., Ferdinand-Wambura, R., Mtelela, K., Kianji, G., Muzuka, A., Gautier, S., Albaric, J., Peyrat, S.: Subsurface images of the Eastern Rift, Africa, from the joint inversion of body waves, surface waves and gravity: investigating the role of fluids in early-stage continental rifting. Geophys. J. Int. 210(2), 931–950 (2017). https://doi.org/10.1093/gji/ggx220
Tromp, J., Komatitsch, D., Liu, Q.: Spectral-element and adjoint methods in seismology. Commun. Comput. Phys. 3(1), 1–32 (2008)
Virieux, J.: P-SV wave propagation in heterogeneous media: velocity–stress finite-difference method. Geophysics 51(4), 889–901 (1986). https://doi.org/10.1190/1.1442147
Xie, Z., Komatitsch, D., Martin, R., Matzen, R.: Improved forward wave propagation and adjoint-based sensitivity kernel calculations using a numerically stable finite-element PML. Geophys. J. Int. 198(3), 1714–1747 (2014). https://doi.org/10.1093/gji/ggu219
Xie, Z., Matzen, R., Cristini, P., Komatitsch, D., Martin, R.: A perfectly matched layer for fluid-solid problems: application to ocean-acoustics simulations with solid ocean bottoms. J. Acoust. Soc. Am. 140(1), 165–175 (2016). https://doi.org/10.1121/1.4954736
Yang, Y., Engquist, B.: Analysis of optimal transport and related misfit functions in full-waveform inversion. Geophysics 83(1), A7–A12 (2018). https://doi.org/10.1190/geo2017-0264.1
Yang, Y., Engquist, B., Sun, J., Hamfeldt, B.F.: Application of optimal transport and the quadratic Wasserstein metric to full-waveform inversion. Geophysics 83(1), R43–R62 (2018). https://doi.org/10.1190/geo2016-0663.1
Zhang, C.H., Xie, Z., Komatitsch, D., Cristini, P., Matzen, R.: Revisit the 1/L problem in rheological models for time-domain seismic-wave propagation. In: Society of Exploration Geophysicists, 3955–3959 (2016). https://doi.org/10.1190/segam2016-13964889.1
Acknowledgements
The authors are grateful to Prof. Martin Mai and Dr. Olaf Zielke, Dr. Luis F.R. Espath and Dr. Håkon Hoel, Prof. Mohammad Motamed and Prof. Daniel Appelö, and Prof. Jesper Oppelstrup for valuable discussions and comments. We are grateful for the support provided by Dr. Samuel Kortas, Computational Scientist, High Performance Computing, KAUST. In particular, we are using the job-scheduler extension decimate 0.9.5 (Kortas 2018), developed by Dr. Kortas. We would also like to acknowledge the use of the open source software package SPECFEM2D (Komatitsch et al. 2012), provided by Computational Infrastructure for Geodynamics (http://geodynamics.org) which is funded by the National Science Foundation under awards EAR-0949446 and EAR-1550901. M. Ballesio, J. Beck, A. Pandey, E. von Schwerin, and R. Tempone are members of the KAUST SRI Center for Uncertainty Quantification in Computational Science and Engineering.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This work is supported by the KAUST Office of Sponsored Research (OSR) under Award No. URF/1/2584-01-01 in the KAUST Competitive Research Grants Program-Round 4 (CRG2015) and the Alexander von Humboldt Foundation. For computer time, it used the resources of the Supercomputing Laboratory at KAUST, under the development Project k1275.
Rights and permissions
About this article
Cite this article
Ballesio, M., Beck, J., Pandey, A. et al. Multilevel Monte Carlo acceleration of seismic wave propagation under uncertainty. Int J Geomath 10, 22 (2019). https://doi.org/10.1007/s13137-019-0135-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s13137-019-0135-5
Keywords
- Multilevel Monte Carlo
- Propagation of uncertainty
- Seismic wave propagation
- Partial differential equations with random data