Regularisation parameter adaptive selection and its application in the prestack AVO inversionGuangtan Huang 1 2 Jingye Li 1 2 Cong Luo 1 2 Xiaohong Chen 1 2 3
1 State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum, Beijing 102249, China.
2 National Engineering Laboratory for Offshore Oil Exploration, China University of Petroleum, Beijing 102249, China.
3 Corresponding author. Email: firstname.lastname@example.org
Exploration Geophysics - https://doi.org/10.1071/EG16100
Submitted: 29 June 2016 Accepted: 8 February 2017 Published online: 21 March 2017
The numerical solution of the inverse problem is usually obtained by solving a set of linear algebraic equations, while the system of equations may suffer from ill-posedness due to insufficient data. Regularisation is a technique for making the estimation problems well posed by adding indirect constraints on the estimated model, but the regularisation parameter selection is difficult. In geophysics, without explicit calculation methods and quantitative evaluation criteria, it is usually based on the experience of the inversion engineers to try to achieve the best inversion results by continuously modifying the regularisation parameter. For prestack amplitude variation with offset (AVO) inversion for real seismic data, fixed regularisation parameters cannot satisfy the optimisation conditions in the seismic data with different signal-to-noise (SNR) in one area. Besides, fixed regularisation parameter may cause that the model constraint misfit is too large or too small compared to data misfit, which may guide the inversion to generate undesirable results. Therefore, adaptive selection of regularisation parameter according to the seismic data can help guarantee a good inversion result. Based on the traditional L-curve criterion, we derive the theoretical formula of the adaptive computation of the regularisation parameter, which can be applied to any norm constraint. We proposed the application of this selection scheme in prestack AVO inversion. Model tests show that the improved L-curve method has better stability than its main competitor, the generalised cross-validation (GCV) method. Prestack AVO inversion on logging data and real seismic data demonstrate that the proposed method can improve the accuracy of the inversion, and it is more immune to strong noise.
Key words: generalised cross-validation, L-curve, prestack inversion, regularisation parameter, sparse constraint condition.
ReferencesAki, K., and Richards, P. G., 1980, Quantitative seismology: theory and methods: W. H. Freeman.
Blumensath, T., and Davies, M. E., 2008, Iterative thresholding for sparse approximations: The Journal of Fourier Analysis and Applications, 14, 629–654
| Iterative thresholding for sparse approximations:CrossRef |
Daubechies, I., Defrise, M., and De Mol, C., 2004, An iterative thresholding algorithm for linear inverse problems with a sparsity constraint: Communications on Pure and Applied Mathematics, 57, 1413–1457
| An iterative thresholding algorithm for linear inverse problems with a sparsity constraint:CrossRef |
Fatti, J. L., Smith, G. C., Vail, P. J., Strauss, P. J., and Levitt, P. R., 1994, Detection of gas in sandstone reservoirs using AVO analysis: a 3-D seismic case history using the Geostack technique: Geophysics, 59, 1362–1376
| Detection of gas in sandstone reservoirs using AVO analysis: a 3-D seismic case history using the Geostack technique:CrossRef |
Gardner, G., Gardner, L. W., and Gregory, A. R., 1974, Formation velocity and density: the diagnostic basis for stratigraphic traps: Geophysics, 39, 770–780
| Formation velocity and density: the diagnostic basis for stratigraphic traps:CrossRef |
Golub, G. H., Heath, M, and Wahba, G., 1979, Generalized cross-validation as a method for choosing a good ridge parameter: Technometrics, 21, 215–223
| Generalized cross-validation as a method for choosing a good ridge parameter:CrossRef |
Guitton, A., Ayeni, G., and Diaz, E., 2012, Constrained full-waveform inversion by model reparameterization: Geophysics, 77, R117–R127
| Constrained full-waveform inversion by model reparameterization:CrossRef |
Haber, E., 1997, Numerical strategies for the solution of inverse problems: Ph.D. thesis, University of British Columbia.
Hansen, P. C., and O’Leary, D. P., 1993, The use of the L-curve in the regularization of discrete ill-posed problems: SIAM Journal on Scientific Computing, 14, 1487–1503
| The use of the L-curve in the regularization of discrete ill-posed problems:CrossRef |
Hansen, P. C., 2000, The L-curve and its use in the numerical treatment of inverse problems, in P. Johnston, ed., Computational inverse problems in electrocardiology: WIT Press, 119–142.
Huang, G, Chen, X, Li, J, Luo, C, and Wang, B, 2017, Application of an adaptive acquisition regularization parameter based on an improved GCV criterion in pre-stack AVO inversion: Journal of Geophysics and Engineering, 14, 100–112
| Application of an adaptive acquisition regularization parameter based on an improved GCV criterion in pre-stack AVO inversion:CrossRef |
Li, X. C., Song, B., and Gan, L. Z., 2015, The application of improved iterative algorithm to regularization model of image restoration: Acta Electronica Sinica., 43, 1152–1159
Loris, I., Douma, H., Nolet, G., Daubechies, I, and Regone, C, 2010, Nonlinear regularization techniques for seismic tomography: Journal of Computational Physics, 229, 890–905
| Nonlinear regularization techniques for seismic tomography:CrossRef | 1:CAS:528:DC%2BD1MXhsFSgsbnF&md5=5d4fe8f4206b9979dfcda7f8b127065bCAS |
Phillips, D. L., 1962, A technique for the numerical solution of certain integral equations of the first kind: Journal of the Association for Computing Machinery, 9, 84–97
| A technique for the numerical solution of certain integral equations of the first kind:CrossRef |
Sacchi, M. D., and Ulrych, T. J., 1995, High-resolution velocity gathers and offset space reconstruction: Geophysics, 60, 1169–1177
| High-resolution velocity gathers and offset space reconstruction:CrossRef |
Tikhonov, A. N., 1943, On the stability of inverse problems: Comptes Rendus (Doklady) de l’Academie des Sciences de l’URSS, 39, 176–179
Tikhonov, A. N., 1963, Solution of incorrectly formulated problems and the regularization method: Soviet Mathematics – Doklady, 5, 1035–1038
Tikhonov, A. N., and Arsenin, V. Y., 1978, Solution of ill-posed problems: Mathematics of Computation, 32, 1320–1322
| Solution of ill-posed problems:CrossRef |
Twomey, S., 1963, On the numerical solution of Fredholm integral equations of the first kind by the inversion of the linear system produced by quadrature: Journal of the Association for Computing Machinery, 10, 97–101
| On the numerical solution of Fredholm integral equations of the first kind by the inversion of the linear system produced by quadrature:CrossRef |
Wang, Z. J., 2003, Research on the regularization solutions of ill-posed problems in geodesy: Ph.D. thesis, Institute of Surveying and Geophysics, Chinese Academy of Sciences.
Xiang, Y., Yu, P., and Chen, X., 2013, An improved adaptive regularized parameter selection in magneto telluric inversion: Journal of Tongji University (Natural Science), 41, 1429–1434
Xue, Z., and Zhu, H., 2015, Full waveform inversion with sparsity constraint in seislet domain: 85th Annual International Meeting, SEG, Expanded Abstracts, 1382–1387.
Xue, Z., Chen, Y., Fomel, S., and Sun, J., 2016, Seismic imaging of incomplete data and simultaneous-source data using least-squares reverse time migration with shaping regularization: Geophysics, 81, S11–S20
| Seismic imaging of incomplete data and simultaneous-source data using least-squares reverse time migration with shaping regularization:CrossRef |
Zhang, F., Dai, R., and Liu, H., 2014, Seismic inversion based on L1-norm misfit function and total variation regularization: Journal of Applied Geophysics, 109, 111–118
| Seismic inversion based on L1-norm misfit function and total variation regularization:CrossRef |
Zwartjes, P., and Gisolf, A., 2007, Fourier reconstruction with sparse inversion: Geophysical Prospecting, 55, 199–221
| Fourier reconstruction with sparse inversion:CrossRef |