Comptes Rendus
Parameter sensitivity analysis of a nonlinear least-squares optimization-based anelastic full waveform inversion method
Comptes Rendus. Mécanique, Inverse problems, Volume 338 (2010) no. 7-8, pp. 364-376.

In a recent article, we described a seismic inversion method for determining the crustal velocity and attenuation of basins in earthquake-prone regions. We formulated the problem as a constrained nonlinear least-squares optimization problem in which the constraints are the equations that describe the forward wave propagation. Here, we conduct a parametric study to investigate the influence of parameters such as the form of the regularization function, receiver density, preconditioning, noise level of the data, and the multilevel continuation technique on the cost and quality of the inversion. We use the same 2D Los Angeles example as in our earlier study.

Nous avons décrit dans un article récent une méthode d'inversion sismique, formulée comme un problème de moindres carrés non linéaires assujetti à des contraintes correspondant aux équations décrivant le problème direct de propagation d'ondes, permettant de déterminer la célérité et l'atténuation crustales de bassins situés dans des zones à forte seismicité. Dans cet article, une étude paramétrique est menée afin d'évaluer l'influence sur le coût et la qualité de l'inversion de paramètres tels que la forme de la fonctionnelle de régularisation, la densité de capteurs, le préconditionnement, le niveau de bruit entachant les données et la méthode de raffinement progressif de l'espace des inconnues. L'exemple utilisé (modèle 2D du bassin de Los Angeles) est le même que dans l'article précédent.

Published online:
DOI: 10.1016/j.crme.2010.07.002
Keywords: Waves, Waveform inversion, Adjoint formulation, Intrinsic attenuation
Mots-clés : Ondes, Inversion par forme d'onde, Formulation adjointe, Atténuation intrinsèque

Aysegul Askan 1; Volkan Akcelik 2; Jacobo Bielak 2; Omar Ghattas 3

1 Department of Civil Engineering and Earthquake Engineering Research Center, Middle East Technical University, Ankara, 06531, Turkey
2 Computational Seismology Laboratory, Department of Civil and Environmental Engineering, Carnegie Mellon University, Pittsburgh, PA 15213, USA
3 Jackson School of Geosciences, Department of Mechanical Engineering, and Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX 78712, USA
@article{CRMECA_2010__338_7-8_364_0,
     author = {Aysegul Askan and Volkan Akcelik and Jacobo Bielak and Omar Ghattas},
     title = {Parameter sensitivity analysis of a nonlinear least-squares optimization-based anelastic full waveform inversion method},
     journal = {Comptes Rendus. M\'ecanique},
     pages = {364--376},
     publisher = {Elsevier},
     volume = {338},
     number = {7-8},
     year = {2010},
     doi = {10.1016/j.crme.2010.07.002},
     language = {en},
}
TY  - JOUR
AU  - Aysegul Askan
AU  - Volkan Akcelik
AU  - Jacobo Bielak
AU  - Omar Ghattas
TI  - Parameter sensitivity analysis of a nonlinear least-squares optimization-based anelastic full waveform inversion method
JO  - Comptes Rendus. Mécanique
PY  - 2010
SP  - 364
EP  - 376
VL  - 338
IS  - 7-8
PB  - Elsevier
DO  - 10.1016/j.crme.2010.07.002
LA  - en
ID  - CRMECA_2010__338_7-8_364_0
ER  - 
%0 Journal Article
%A Aysegul Askan
%A Volkan Akcelik
%A Jacobo Bielak
%A Omar Ghattas
%T Parameter sensitivity analysis of a nonlinear least-squares optimization-based anelastic full waveform inversion method
%J Comptes Rendus. Mécanique
%D 2010
%P 364-376
%V 338
%N 7-8
%I Elsevier
%R 10.1016/j.crme.2010.07.002
%G en
%F CRMECA_2010__338_7-8_364_0
Aysegul Askan; Volkan Akcelik; Jacobo Bielak; Omar Ghattas. Parameter sensitivity analysis of a nonlinear least-squares optimization-based anelastic full waveform inversion method. Comptes Rendus. Mécanique, Inverse problems, Volume 338 (2010) no. 7-8, pp. 364-376. doi : 10.1016/j.crme.2010.07.002. https://comptes-rendus.academie-sciences.fr/mecanique/articles/10.1016/j.crme.2010.07.002/

[1] A. Tarantola Inversion of seismic reflection data in the acoustic approximation, Geophysics, Volume 49 (1984), pp. 1259-1266

[2] A. Tarantola Inverse Problem Theory and Model Parameter Estimation, SIAM, Philadelphia, PA, 2005

[3] A. Tarantola; B. Valette Generalized nonlinear inverse problems solved using the least squares criterion, Rev. Geophys. Space Phys., Volume 20 (1982), pp. 219-232

[4] J. Tromp; C. Tape; Q. Liu Seismic tomography, adjoint methods, time reversal and banana-doughnut kernels, Geophys. J. Int., Volume 160 (2005), pp. 195-216

[5] P. Chen; L. Zhao; T.H. Jordan Full 3D tomography for crustal structure of the Los Angeles region, Bull. Seism. Soc. Am., Volume 97 (2007), pp. 1094-1120

[6] R.G. Pratt; M.H. Worthington Inverse theory applied to multi-source cross-hole tomography. Part I: Acoustic wave-equation method, Geophys. Prospect., Volume 38 (1990), pp. 287-310

[7] L. Sirgue; R.G. Pratt Efficient waveform inversion and imaging: a strategy for selecting temporal frequencies, Geophysics, Volume 69 (2004), pp. 231-248

[8] S. Operto; C. Ravaut; L. Improta; J. Virieux; A. Herrero; P. Dell'Aversana Quantitative imaging of complex structures from dense wide-aperture seismic data by multiscale traveltime and waveform inversions: a case study, Geophys. Prospect., Volume 52 (2004), pp. 625-651

[9] C. Ravaut; S. Operto; L. Improta; J. Virieux; A. Herrero; P. Dell'Aversana Multiscale imaging of complex structures from multifold wide-aperture seismic data by frequency-domain full-waveform tomography: application to a thrust belt, Geophys. J. Int., Volume 159 (2004), pp. 1032-1056

[10] A. Tarantola Theoretical background for the inversion of seismic waveforms including elasticity and attenuation, Pure Appl. Geophys., Volume 128 (1988), pp. 365-399

[11] G.J. Rix; C.G. Lai; A.W. Spang In situ measurements of damping ratio using surface waves, J. Geotech. Geoenviron. Engng., Volume 126 (2000), pp. 472-480

[12] G. Hicks; R.G. Pratt Reflection waveform inversion using local descent methods: estimating attenuation and velocity over a gas-sand deposit, Geophysics, Volume 66 (2001), pp. 598-612

[13] A. Askan; V. Akcelik; J. Bielak; O. Ghattas Full waveform inversion for seismic velocity and anelastic losses in heterogeneous structures, Bull. Seism. Soc. Am., Volume 97 (2007) no. 6, pp. 1990-2008

[14] H. Magistrale; S. Day; R.W. Clayton; R. Graves The SCEC Southern California reference three-dimensional seismic velocity model, Version 2, Bull. Seism. Soc. Am., Volume 90 (2000), p. S65-S76

[15] A. Askan, Full waveform inversion for seismic velocity and anelastic losses in heterogeneous structures, PhD thesis, Carnegie Mellon University, Pittsburgh, PA, 2006.

[16] S.M. Day; C.R. Bradley Memory efficient simulation of anelastic wave propagation, Bull. Seism. Soc. Am., Volume 91 (2001), pp. 520-531

[17] P. Moczo; J.O.A. Robertsson; L. Eisner The finite-difference time-domain method for modelling of seismic wave propagation (Ru-Shan Wu; Valerie Maupin, eds.), Advances in Wave Propagation in Heterogeneous Earth, Advances in Geophysics, vol. 48, Academic Press, New York, MA, 2007, pp. 421-516

[18] G. Biros; O. Ghattas Parallel Lagrange–Newton–Krylov–Schur methods for PDE-constrained optimization. Part I: The Krylov–Schur solver, SIAM J. Sci. Comput., Volume 27 (2005) no. 2, pp. 687-713

[19] G. Biros; O. Ghattas Parallel Lagrange–Newton–Krylov–Schur methods for PDE-constrained optimization. Part II: The Lagrange–Newton solver, and its application to optimal control of steady viscous flows, SIAM J. Sci. Comput., Volume 27 (2005) no. 2, pp. 714-739

[20] V. Akcelik; G. Biros; O. Ghattas; J. Hill; D. Keyes; B. van Bloemen Waanders Parallel algorithms for PDE-constrained optimization (M. Heroux; P. Raghaven; H. Simon, eds.), Frontiers of Parallel Computing, SIAM, 2006

[21] V. Akcelik, G. Biros, O. Ghattas, Parallel multiscale Gauss–Newton–Krylov methods for inverse wave propagation, in: Proceedings of ACM/IEEE SC2002, 2002.

[22] V. Akcelik, J. Bielak, G. Biros, I. Epanomeritakis, A. Fernandez, O. Ghattas, E.J. Kim, J. Lopez, D. O'Hallaron, T. Tu, J. Urbanic, High-resolution forward and inverse earthquake modeling on terascale computers, in: Proceedings of ACM/IEEE SC2003, 2003.

[23] J. Nocedal; S.J. Wright Numerical Optimization, Springer, New York, 1999

[24] R. Acar; R.C. Vogel Analysis of bounded variation penalty methods for ill-posed problems, Inverse Problems, Volume 10 (1994), pp. 1217-1229

[25] C. Vogel Computational Methods for Inverse Problems, SIAM, 2002

[26] P. Craven; G. Wahba Smoothing noisy data with spline functions, Numer. Math., Volume 31 (1979), pp. 377-403

[27] S.C. Eisenstat; H.F. Walker Choosing the forcing terms in an inexact Newton method, SIAM J. Sci. Comput., Volume 17 (1996), pp. 16-32

[28] C. Bunks; F.M. Saleck; S. Zaleski; G. Chavent Multiscale seismic waveform inversion, Geophysics, Volume 50 (1995), pp. 1457-1473

[29] J.L. Morales; J. Nocedal Automatic preconditioning by limited memory quasi-Newton updating, Siam J. Optim., Volume 10 (2000) no. 4, pp. 1079-1096

[30] I. Epanomeritakis; V. Akcelik; O. Ghattas; J. Bielak A Newton–CG method for large-scale three-dimensional elastic full-waveform seismic inversion, Inverse Problems, Volume 24 (2008), p. 034015 (26 pp)

Cited by Sources:

Comments - Policy