Memoryless quasi-Newton (MLQN) method for 2D acoustic full waveform inversion
Cai Liu 1 2 Fengxia Gao 1 Xuan Feng 1 Yang Liu 1 Qianci Ren 1
+ Author Affiliations
- Author Affiliations
1 College of GeoExploration Science and Technology, Jilin University, No. 938, Ximinzhu Street, Changchun, Jilin 130061, China.
2 Corresponding author. Email: liucai@jlu.edu.cn
Exploration Geophysics 46(2) 168-177 https://doi.org/10.1071/EG13090
Submitted: 20 June 2013 Accepted: 13 March 2014 Published: 22 April 2014
Abstract
Full waveform inversion (FWI) is an efficient way to solve parameter reconstruction problems, such as velocity, density, and viscosity coefficient. In this study, we apply a memoryless quasi-Newton (MLQN) method in FWI to invert velocity from surface seismic data for the first time. This method can attain acceptable results with low computational cost and small memory storage requirements. To ensure that the inverted velocity is maintained between the lower and upper boundaries of the velocity model, a nonlinear transformation is added to velocity as a priori information. To test the efficiency of the MLQN method in FWI, two synthetic models, a modified Marmousi model and a modified overthrust model, are examined from the surface seismic data with and without white Gaussian noise. For comparison, the conjugate gradient (CG) method is carried out for the same velocity models with the same parameters. We compare the inverted velocities by the two methods based on the aspects of memory storage requirements, computation time for each iteration, and error. By keeping the memory storage requirements and computation time in each iteration similar, the reconstructed velocity models obtained using the MLQN method are closer to the true velocity models than those obtained using the CG method. Our numerical tests show that the MLQN method is feasible and reliable in FWI.
Key words: CG method, full waveform inversion, memoryless quasi-Newton method, velocity.
References
Anagaw, A. Y., and Sacchi, M. D., 2012, Full waveform inversion with simultaneous sources using the full Newton Method: 82nd Annual International Meeting, SEG, Expanded Abstracts, 1–5.Bae, H. S., Pyun, S., Chung, W., Kang, S. G., and Shin, C., 2012, Frequency-domain acoustic-elastic coupled waveform inversion using the Gauss-Newton conjugate gradient method: Geophysical Prospecting, 60, 413–432
| Frequency-domain acoustic-elastic coupled waveform inversion using the Gauss-Newton conjugate gradient method:Crossref | GoogleScholarGoogle Scholar |
Brossier, R., 2011, Two-dimensional frequency-domain visco-elastic full waveform inversion: parallel algorithms, optimization and performance: Computers & Geosciences, 37, 444–455
| Two-dimensional frequency-domain visco-elastic full waveform inversion: parallel algorithms, optimization and performance:Crossref | GoogleScholarGoogle Scholar |
Brossier, R., Operto, S., and Virieux, J., 2009, Seismic imaging of complex onshore structures by 2D elastic frequency-domain full-waveform inversion: Geophysics, 74, WCC105–WCC118
| Seismic imaging of complex onshore structures by 2D elastic frequency-domain full-waveform inversion:Crossref | GoogleScholarGoogle Scholar |
Gauthier, O., Virieux, J., and Tarantola, A., 1986, Two-dimensional nonlinear inversion of seismic waveforms: numerical results: Geophysics, 51, 1387–1403
| Two-dimensional nonlinear inversion of seismic waveforms: numerical results:Crossref | GoogleScholarGoogle Scholar |
Greenhalgh, S. A., and Zhou, B., 2003, Surface seismic imaging by multi-frequency amplitude inversion: Exploration Geophysics, 34, 217–224
| Surface seismic imaging by multi-frequency amplitude inversion:Crossref | GoogleScholarGoogle Scholar |
Greenhalgh, S. A., Zhou, B., and Green, A., 2006, Solutions, algorithms and inter-relations for local minimization search geophysical inversion: Journal of Geophysics and Engineering, 3, 101–113
| Solutions, algorithms and inter-relations for local minimization search geophysical inversion:Crossref | GoogleScholarGoogle Scholar |
Habashy, T. M., and Abubakar, A., 2004, A general framework for constraint minimization for the inversion of electromagnetic measurements: Progress in Electromagnetics Research, 46, 265–312
| A general framework for constraint minimization for the inversion of electromagnetic measurements:Crossref | GoogleScholarGoogle Scholar |
Hu, W., Abubakar, A., and Habashy, T. M., 2009, Simultaneous multifrequency inversion of full-waveform seismic data: Geophysics, 74, R1–R14
| Simultaneous multifrequency inversion of full-waveform seismic data:Crossref | GoogleScholarGoogle Scholar |
Hustedt, B., Operto, S., and Virieux, J., 2004, Mixed-grid and staggered-grid finite-difference methods for frequency-domain acoustic wave modeling: Geophysical Journal International, 157, 1269–1296
| Mixed-grid and staggered-grid finite-difference methods for frequency-domain acoustic wave modeling:Crossref | GoogleScholarGoogle Scholar |
Liu, C., Gao, F. X., Feng, X., Liu, Y., and Liu, Y. W., 2013, Incorporating attenuation effects into frequency-domain full waveform inversion from zero-offset VSP data from Stokes equation: Journal of Geophysics and Engineering, 10, 035004
| Incorporating attenuation effects into frequency-domain full waveform inversion from zero-offset VSP data from Stokes equation:Crossref | GoogleScholarGoogle Scholar |
Ma, Y., and Hale, R., 2012, Quasi-Newton full-waveform inversion with a projected Hessian matrix: Geophysics, 77, R207–R216
| Quasi-Newton full-waveform inversion with a projected Hessian matrix:Crossref | GoogleScholarGoogle Scholar |
Malinowski, M., Operto, S., and Ribodetti, A., 2011, High-resolution seismic attenuation imaging from wide-aperture onshore data by visco-acoustic frequency-domain full-waveform inversion: Geophysical Journal International, 186, 1179–1204
| High-resolution seismic attenuation imaging from wide-aperture onshore data by visco-acoustic frequency-domain full-waveform inversion:Crossref | GoogleScholarGoogle Scholar |
Métivier, L., Brossier R., Virieux J., Operto S., 2012. The truncated Newton method for full waveform inversion: 82nd Annual International Meeting, SEG, Expanded Abstracts, 1–5 .
Mora, P., 1988, Elastic wave-field inversion of reflection and transmission data: Geophysics, 53, 750–759
| Elastic wave-field inversion of reflection and transmission data:Crossref | GoogleScholarGoogle Scholar |
Nocedal, J., and Wright, J. S., 2006, Numerical optimization (2nd edition): Springer-Verlag.
Pratt, R. G., 1999, Seismic waveform inversion in the frequency domain, Part 1: theory and verification in a physical scale model: Geophysics, 64, 888–901
| Seismic waveform inversion in the frequency domain, Part 1: theory and verification in a physical scale model:Crossref | GoogleScholarGoogle Scholar |
Pratt, R. G., Shin, C., and Hicks, J. G., 1998, Gauss-Newton and full Newton methods in frequency-space seismic waveform inversion: Geophysical Journal International, 133, 341–362
| Gauss-Newton and full Newton methods in frequency-space seismic waveform inversion:Crossref | GoogleScholarGoogle Scholar |
Ravaut, C., Operto, S., Improta, L., Virieux, J., Herrero, A., and Dell’Aversana, P., 2004, Multiscale imaging of complex structures from multifold wide-aperture seismic data by frequency-domain full-waveform tomography: application to a thrust belt: Geophysical Journal International, 159, 1032–1056
| Multiscale imaging of complex structures from multifold wide-aperture seismic data by frequency-domain full-waveform tomography: application to a thrust belt:Crossref | GoogleScholarGoogle Scholar |
Shanno, F. D., 1978, Conjugate gradient methods with inexact searches: Mathematics of Operations Research, 3, 244–256
| Conjugate gradient methods with inexact searches:Crossref | GoogleScholarGoogle Scholar |
Shi, Y. M., Zhao, W. Z., and Cao, H., 2007, Nonlinear process control of wave-equation inversion and its application in the detection of gas: Geophysics, 72, R9–R18
| Nonlinear process control of wave-equation inversion and its application in the detection of gas:Crossref | GoogleScholarGoogle Scholar |
Shin, C., Yoon, K., Marfurt, J. K., Park, K., Yang, D., Lim, Y. H., Chung, S., and Shin, S., 2001, Efficient calculation of a partial-derivative wavefield using reciprocity for seismic imaging and inversion: Geophysics, 66, 1856–1863
| Efficient calculation of a partial-derivative wavefield using reciprocity for seismic imaging and inversion:Crossref | GoogleScholarGoogle Scholar |
Sourbier, F., Operto, S., and Virieux, J., 2009, FWT2D: a massively parallel program for frequency-domain full-waveform tomography of wide-aperture seismic data - Part 1: algorithm: Computers & Geosciences, 35, 487–495
| FWT2D: a massively parallel program for frequency-domain full-waveform tomography of wide-aperture seismic data - Part 1: algorithm:Crossref | GoogleScholarGoogle Scholar |
Tarantola, A., 1984, Inversion of seismic reflection data in the acoustic approximation: Geophysics, 49, 1259–1266
| Inversion of seismic reflection data in the acoustic approximation:Crossref | GoogleScholarGoogle Scholar |
Vigh, D., Starr, E. W., and Kapoor, J., 2009, Developing earth models with full waveform inversion: The Leading Edge, 28, 432–435
| Developing earth models with full waveform inversion:Crossref | GoogleScholarGoogle Scholar |
Xu, C. X., Chen, Z. P., and Li, N. C., 2002, Modern optimization methods: Science Press (Beijing).
Yu, J. Y., 1990, A memoryless quasi-newton method for solving unconstrained optimizaiton problems: Computation Mathmatics, 3, 259–269
Zhou, B., and Greenhalgh, A. S., 2003, Crosshole seismic inversion with normalized full-waveform amplitude data: Geophysics, 68, 1320–1330
| Crosshole seismic inversion with normalized full-waveform amplitude data:Crossref | GoogleScholarGoogle Scholar |