PP- & PS-wave prestack nonlinear inversion based on adaptive MCMC algorithm
Zhang Guangzhi1,2, Pan Xinpeng1, Sun Changlu1, Yin Xingyao1,2
1. School of Geosciences, China University of Petroleum(East China), Qingdao, Shandong 266580, China;
2. Laboratory of Marine Mineral Resources, Qingdao National Laboratory of Marine Science and Technology, Qingdao, Shandong 266071, China
Abstract:Through amplitude versus offset (AVO) inversion, we can obtain elastic parameters related to reservoir lithology and fluid content. However, the solutions of prestack inversion with PP data only are quite unstable, so the joint application of PP and PS seismic data could further reduce the non-uniqueness of the inversion, and greatly improve its stability and accuracy of the inversion. In this article, we develop a new method of joint prestack nonlinear inversion of PP and PS waves based on the adaptive Markov chain Monte Carlo (MCMC) algorithm, and it can be used for inverting three parameters including P- and S-wave velocities and density via the exact Zoeppritz equation. Based on Bayesian framework, logging-constrained priori information is introduced in the inversion process. Then we sample posterior probability density and apply the adaptive MCMC algorithm to obtain the posterior probability density of AVO three parameters by making statistical analysis of these random samples. Finally, we analyze the uncertainty of AVO inversion results of three parameters. Tests on synthetic data show that all three parameters are well retrieved, and the proposed method is quite stable, accurate, and anti-noise, which demonstrates its reliability and effectiveness.
Hou Dongjia,Liu Yang,Hu Guoqing et al. Pre-stack multiwave joint inversion for elastic moduli based on Bayesian theory. Chinese Journal of Geophysics,2014,57(4):1251-1264.
Zhang Guangzhi,Du Bingyi,Li Haishan et al. The method of joint pre-stack inversion of PP and P-SV waves in shale gas reservoirs.Chinese Journal of Geo-physics,2014,57(12):4141-4149.
[11]
Zhi L X,Chen S Q and Li X Y. Joint AVO inversion of PP and PS waves using exact Zoeppritz equation. SEG Technical Program Expanded Abstract,2013,32:457-461.
[12]
Larsen J A. AVO Inversion by Simultaneous P-P and P-S Inversion. University of Calgary,1999.
[13]
Chen T S,Wei X C and Jia L H. Nonlinear simultaneous PP and PS AVO inversion based on Zoeppritz equations. Beijing 2009 International Geophysical Conference and Exposition,2009,24-27.
[14]
Chen T S and Wei X C. Zoeppritz-based joint AVO inversion of PP and PS waves in ray parameter domain. SEG Technical Program Expanded Abstract,2012,31:1-5.
[15]
Liu H X,Li J Y,Chen X H et al. Joint AVO inversion of PP and PS reflections using the reflectivity method. 76th EAGE Conference & Exhibition,2014.
Zhang Fengqi,Wei Fuji,Wang Yanchun et al.Genera-lized linear AVO inversion with the priori constraint of trivariate Cauchy distribution based on Zoeppritz equation. Chinese Journal of Geophysics,2013,56(6):2098-2115.
Huang Handong,Wang Yanchao,Guo Fei et al. High precision pre-stack inversion algorithm based on Zoeppritz equations. OGP,2013,48(5):740-749.
[18]
Liu Hongxing,Li Jingye, Chen Xiaohong et al. Joint AVO inversion of PP and PS reflections using the reflectivity method. 76th EAGE Conference & Exhibition,2014.
[19]
Hastings W K. Monte Carlo sampling methods using Markov chains and their applications. Biometrika,1970,57(1):97-109.
[20]
Metropolis N,Rosenbluth A W,Rosenbluth M N et al. Equations of state calculations by fast computing machines. The Journal of Chemical Physics,1953,21(6):1087-1092.
Zhang Guangzhi,Wang Danyang,Yin Xingyao et al. Study on pre-stack seismic inversion using Markov chain Monte Carlo. Chinese Journal of Geophysics,2011,54(11):2926-2932.
Zhang Guangzhi,Wang Danyang and Yin Xingyao. Seismic parameter estimated using Markov chain Monte Carlo method. OGP,2011,46(4):605-609.
[26]
Bachrach R,Sengupta M,Salama A et al. Reconstruction of the layer anisotropic elastic parameter and high resolution fracture characterization from P-wave data:A case study using seismic inversion and Bayesian rock physics parameter estimation. Geophysical Prospecting,2009,57(2):253-262.
[27]
Grana D and Rossa E D. Probabilistic petrophysical properties estimation integrating statistical rock phy-sics with seismic inversion. Geophysics,2010,75(3):O21-O37.
[28]
Ulvmoen M and Omre H. Improved resolution in Bayesian lithology/fluid inversion from pre-stack seismic data and well observations:Part 1-Methodology. Geophysics,2010,75(2):R21-R35
[29]
Rimstad K,Avseth P and Omre H. Hierarchical Bayesian lithology prediction:A North Sea case study.Geophysics,2012,77(2):B39-B85.
[30]
Haario H,Saksman E and Tamminen J. Adaptive proposal distributions for random walk Metropolis algorithm.Computational Statistics,1999,14(3):375-396.
[31]
Haario H,Saksman E and Tamminen J. An adaptive Metropolis algorithm. Bernoulli,2001,7(2):223-242.
[32]
Gelman A,Roberts G and Gilks W. Efficient Metropolis jumping rules.Bayesian Statistic,1996,5:599-607.
[33]
Downton J E and Ursenbach C. Linearized amplitude variation with offset (AVO) inversion with supercri-tical angles. Geophysics,2006,71(5):E49-E55.
[34]
Skopintseva L,Ayzenberg M,Landrø M et al. Long-offset AVO inversion of PP reflections from plane interfaces using effective reflection coefficients. Geophysics,2011,76(6):C65-C79.
[35]
Aki K and Richards P G. Quantitative Seismology:Theory and Method. W H Freeman and Co,1980,123-192.