Least-squares seismic inversion with stochastic conjugate gradient method

Abstract

With the development of computational power, there has been an increased focus on data-fitting related seismic inversion techniques for high fidelity seismic velocity model and image, such as full-waveform inversion and least squares migration. However, though more advanced than conventional methods, these data fitting methods can be very expensive in terms of computational cost. Recently, various techniques to optimize these data-fitting seismic inversion problems have been implemented to cater for the industrial need for much improved efficiency. In this study, we propose a general stochastic conjugate gradient method for these data-fitting related inverse problems. We first prescribe the basic theory of our method and then give synthetic examples. Our numerical experiments illustrate the potential of this method for large-size seismic inversion application.

Keywords

least-squares seismic inversion / stochastic conjugate gradient method / data fitting / Kirchhoff migration

References

[1] Claerbout J. F. Basic Earth Imaging, 2010, 67-69.

[2] Dai W., Wang X., Schuster G. T. Least-Squares Migration of Multisource Data with a Deblurring Filter. Geophysics, 2011, 76(5): R135-R146.

[3] Fletcher R., Reeves C. Function Minimization by Conjugate Gradients. The Computer Journal, 1964, 7(2): 149-154.

[4] Godwin J., Sava P. Blended Source Imaging by Amplitude Encoding. SEG Technical Program Expanded Abstracts, 2010, 3125-3129.

[5] Huang W., Ma H. D., Vigh D., . Velocity Model Building with Long-Offset and Full-Azimuth Data: A Case History for Full-Waveform Inversion. SEG Technical Program Expanded Abstracts, 2013, 32: 4750-4754.

[6] Jiang H., Wilford P. A Stochastic Conjugate Gradient Method for the Approximation of Functions. Journal of Computational and Applied Mathematics, 2012, 236(9): 2529-2544.

[7] Krebs J. R., Anderson J. E., Hinkley D., . Fast Full-Wavefield Seismic Inversion Using Encoded Sources. Geophysics, 2009, 74(6): WCC177-WCC188.

[8] Nemeth T., Wu C. J., Schuster G. T. Least-Squares Migration of Incomplete Reflection Data. Geophysics, 1999, 64(1): 208-221.

[9] Saad Y. Iterative Methods for Sparse Linear Systems. SIAM, 2003, 105-127.

[10] Schraudolph N. N., Graepel T. Combining Conjugate Direction Methods with Stochastic Approximation of Gradients. Proc. 9th Intl. Workshop Artificial Intelligence and Statistics (AIstats), Society for Artificial Intelligence and Statistic, 2003, 7-13.

[11] Suh S. Y., Yeh A., Wang B., . Cluster Programming for Reverse Time Migration. The Leading Edge, 2010, 29(1): 94-97.

[12] Tarantola A. Inversion of Seismic Reflection Data in the Acoustic Approximation. Geophysics, 1984, 49: 1259-1266.

[13] Tibshirani R. Regression Shrinkage and Selection via the LASSO. J. Royal Statistical Society, 1996, 58(1): 267-288.

[14] Leeuwen T. V., Aravkin A. Y., Herrmann F. J. Seismic Waveform Inversion by Stochastic Optimization. International Journal of Geophysics, 2011, 1-18.

[15] Vigh D., Starr E. W. 3D Prestack Plane-Wave, Full-Waveform Inversion. Geophysics, 2008, 73(5): VE135-VE144.