Monte Carlo-Markov Chain stochastic inversion constrained by seismic waveform
ZHOU Shuangshuang1,2, YIN Xingyao1,2, PEI Song1,2, YANG Yaming1,2
1. School of Geosciences, China University of Petroleum(East China), Qingdao, Shandong 266580, China; 2. Laboratory for Marine Mineral Resources, Qingdao National Laboratory for Marine Science and Technology, Qingdao, Shandong 266071, China
Abstract:The resolution of seismic stochastic inversion based on geostatistics and logging data is higher than that of conventional deterministic inversion, so the former is quickly and widely used, but it is difficult to improve calculation efficiency and eliminate randomness. This paper proposes a Monte Carlo-Markov Chain (MCMC) stochastic inversion method based on the constraint of seismic waveform. By making full use of the geophysical mapping relationship between seismic data and parameters to be inversed, and a correlation coefficient to guide pseudo ordinary Kriging interpolation to well data according to the similarity of known seismic waveforms, an initial model is established; then the posterior probability density distribution is constructed under the constraints of seismic data and logging data on the Bayesian framework, and the initial model which can indicate seismic waveforms is randomly simulated multiple times by using the Metropolis-Hastings sampling algorithm. The posterior mean value is the optimal solution to the model parameters. This method effectively improves inversion stability and lateral continuity, reduces randomness, effectively weakens the impact of seismic noises on inversion results, and greatly accelerates the convergence of the Markov chain, which effectively improves computing efficiency and estimation accuracy. Applications on model and real data have proved the MCMC stochastic inversion method constrained by seismic waveforms has good noise resistance, can effectively improve inversion accuracy, and are advantageous in identifying thin reservoirs within a tuning scale. It improves both vertical resolution and ho-rizontal resolution.
印兴耀,张世鑫,张繁昌,等.利用基于Russell近似的弹性波阻抗反演进行储层描述和流体识别[J].石油地球物理勘探,2010,45(3):373-380.YIN Xingyao, ZHANG Shixin, ZHANG Fanchang,et al. Utilizing Russell approximation-based elastic wave impedance inversion to conduct reservoir description and fluid identification[J].Oil Geophysical Prospecting, 2010, 45(3):373-380.
[2]
印兴耀,曹丹平,王保丽,等. 基于叠前地震反演的流体识别方法研究进展[J].石油地球物理勘探, 2014, 49(1):22-34.YIN Xingyao, CAO Danping, WANG Baoli, et al. Research progress of fluid discrimination with pre-stack seismic inversion[J].Oil Geophysical Prospecting, 2014, 49(1):22-34.
[3]
印兴耀,宗兆云,吴国忱.岩石物理驱动下地震流体识别研究[J].中国科学:地球科学,2015,45(1):8-21.YIN Xingyao, ZONG Zhaoyun, WU Guochen. Research on seismic fluid identification driven by rock physics[J].Science China:Earth Sciences,2015,45(1):8-21.
[4]
Li K, Yin X Y, Liu J, et al. An improved stochastic inversion for joint estimation of seismic impedance and lithofacies[J].Journal of Geophysics and Engineering,2019,16(1):62-76.
[5]
高君,毕建军,赵海山,等.地震波形指示反演薄储层预测技术及其应用[J].地球物理学进展,2017,32(1):142-145.GAO Jun,BI Jianjun,ZHAO Haishan,et al.Seismic waveform inversion technology and application of thinner reservoir prediction[J].Progress in Geophy-sics, 2017,32(1):142-145.
[6]
Yin X Y,Li K,Zong Z Y.Resolution enhancement of robust Bayesian pre-stack inversion in the frequency domain[J].Journal of Geophysics and Engineering,2016,13(5):646-656.
[7]
Li K,Yin X Y,Zong Z Y.Bayesian seismic multi-scale inversion in complex Laplace mixed domains[J].Petroleum Science,2017,14(4):694-710.
[8]
Liu X X,Yin X Y,Li H S.Optimal variable-grid finite-difference modeling for porous media[J].Journal of Geophysics and Engineering,2014,doi:10.1088/1742-2132/11/6/065011.
[9]
王宇,韩立国,周家雄,等. L1-L2范数联合约束稀疏脉冲反演的应用[J].地球科学——中国地质大学学报, 2009,34(5):835-840.WANG Yu,HAN Liguo,ZHOU Jiaxiong,et al.Application of combined norm constrained sparse spike inverse[J].Earth Sciences:Journal of China University of Geosciences, 2009,34(5):835-840.
[10]
郭朝斌,杨小波,陈红岳,等.约束稀疏脉冲反演在储层预测中的应用[J].石油物探, 2006,45(4):397-400.GUO Chaobin,YANG Xiaobo,CHEN Hongyue,et al. Constrained sparse pulse inversion research in north of Haitongji depression[J].Geophysical Prospecting for Petroleum, 2006, 45(4):397-400.
[11]
赵晨,张广智,张佳佳,等.基于Metropolis优化的叠前全局迭代地质统计学反演方法[J].地球物理学报,2020,63(8):3116-3130.ZHAO Chen,ZHANG Guangzhi,ZHANG Jiajia,et al.Prestack global iteration geostatistical inversion method based on Metropolis sampling algorithm[J].Chinese Journal of Geophysics,2020,63(8):3116-3130.
[12]
孙瑞莹.先验信息构建与地震随机反演方法研究[D].山东青岛:中国石油大学(华东), 2015.SUN Ruiying.Study of the Priori Information Construction and Seismic Stochastic Inversion Method[D].China University of Petroleum (East China),Qingdao,Shandong, 2015.
[13]
杨锴,艾迪飞,耿建华.测井、井间地震与地面地震数据联合约束下的地质统计学随机建模方法研究[J].地球物理学报,2012,55(8):2695-2704.YANG Kai,Ai Difei,Geng Jianhua.A new geostatistical inversion and reservoir modeling technique constrained by well-log,crosshole and surface seismic data[J].Chinese Journal of Geophysics,2012,55(8):2695-2704.
[14]
王小丹,印兴耀,金惠,等.叠前地震随机反演方法及实际资料应用[J].地球物理学进展,2018,33(6):2471-2476.WANG Xiaodan,YIN Xingyao,JIN Hui,et al.Pre-stack seismic stochastic inversion and application on real data[J].Progress in Geophysics, 2018, 33(6):2471-2476.
李宁.基于模拟退火的地质统计学反演方法研究[D].山东青岛:中国石油大学(华东), 2013.LI Ning.The Study of Geostatistics Inversion Based on Simulated Annealing Method[D].China University of Petroleum (East China),Qingdao,Shandong,2013.
[18]
印兴耀,贺维胜,黄旭日.贝叶斯-序贯高斯模拟方法[J].石油大学学报(自然科学版),2005,29(5):28-32.YIN Xingyao,HE Weisheng,HUANG Xuri.Bayesian sequential Gaussian simulation method[J].Journal of the University of Petroleum,China(Edition of Natural Science),2005,29(5):28-32.
[19]
Yin X Y,Sun R Y,Wang B L,et al.Simultaneous inversion of petrophysical parameters based on a geostatistical priori information[J].Applied Geophysics,2014,11(3):311-320.
[20]
王保丽,印兴耀,丁龙翔,等.基于FFT-MA谱模拟的快速随机反演方法研究[J].地球物理学报,2015,58(2):664-673.WANG Baoli,YIN Xingyao,DING Longxiang,et al.Study of fast stochastic inversion based on FFT-MA spectrum simulation[J].Chinese Journal of Geophy-sics,2015,58(2):664-673.
[21]
潘新朋.优化MCMC方法在地震反演中的应用研究[D].山东青岛:中国石油大学(华东),2016.PAN Xinpeng.The Application Study in Seismic Inversion Based on Optimized Markov Chain Monte Carlo Method[D].China University of Petroleum (East China), Qingdao,Shandong,2016.
[22]
王丹阳.基于MCMC方法的叠前反演方法研究[D].山东青岛:中国石油大学(华东),2012.WANG Danyang. Pre-stack Seismic Inversion Based on Markov Chain Monte Carlo Method[D].China University of Petroleum (East China),Qingdao,Shandong,2012.
[23]
张广智,王丹阳,印兴耀.利用MCMC方法估算地震参数[J].石油地球物理勘探,2011,46(4):605-609.ZHANG Guangzhi,WANG Danyang,YIN Xingyao. Seismic parameter estimation using Markov Chain Monte Carlo Method[J].Oil Geophysical Prospecting,2011, 46(4):605-609.
[24]
Hong T C,Sen M K. A new MCMC algorithm for seismic waveform inversion and corresponding uncertainty analysis[J].Geophysical Journal International,2009, 177(2):14-32.
[25]
Metropolis N,Rosenbluth A W,Rosenbluth M N,et al. Equation of state calculations by fast computing machines[J].The Journal of Chemical Physics,1953,21(6):1087-1092.
[26]
Hastings W K.Monte Carlo sampling methods using Markov chains and their applications[J].Biometrika,1970, 57(1):97-109.
[27]
Smith A F M,Roberts G O.Bayesian computation via the Gibbs sampler and related Markov Chain Monte Carlo methods[J].Journal of the Royal Statistical Society,1993, 55(1):3-23.
[28]
张广智,潘新朋,孙昌路,等.纵横波联合叠前自适应MCMC反演方法[J].石油地球物理勘探,2016,51(5):938-946.ZHANG Guangzhi,PAN Xinpeng,SUN Changlu,et al.PP- & PS-wave prestack nonlinear inversion based on adaptive MCMC algorithm[J].Oil Geophysical Prospecting,2016, 51(5):938-946.
[29]
张广智,王丹阳,印兴耀.基于MCMC的叠前地震反演方法研究[J].地球物理学报,2011,54(11):2926-2932.ZHANG Guangzhi,WANG Danyang,YIN Xingyao.Study on prestack seismic inversion using Markov Chain Monte Carlo[J].Chinese Journal of Geophy-sics,2011,54(11):2926-2932.
[30]
Pan X P,Zhang G Z,Chen H Z,et al.McMC-based AVAZ direct inversion for fracture weaknesses[J].Journal of Applied Geophysics,2017,138(7):50-61.
[31]
Pan X P,Zhang G Z,Chen H Z,et al. McMC-based nonlinear EIVAZ inversion driven by rock physics[J].Journal of Geophysics and Engineering,2017,14(2):368-379.
[32]
Pan X P,Zhang G Z,Zhang J J,et al. Zoeppritz-based AVO inversion using an improved Markov chain Monte Carlo method[J].Petroleum Science,2017,14(1):75-83.
[33]
李坤,印兴耀,宗兆云.岩石物理驱动的相约束叠前地震概率化反演方法[J].中国科学:地球科学,2020,50(6):832-854.LI Kun,YIN Xingyao,ZONG Zhaoyun.Facies-constrained prestack seismic probabilistic inversion dri-ven by rock physics[J].Science China:Earth Sciences,2020,50(6):832-854.
[34]
章雄,张本健,梁虹,等.波形指示叠前地震反演方法在致密含油薄砂层预测中的应用[J].物探与化探,2018,42(3):545-554.ZHANG Xiong,ZHANG Benjian,LIANG Hong,et al.The application of pre-stack inversion based on seismic waveform indicator to the prediction of compact and thin oil-bearing sand layer[J].Geophysical and Geochemical Exploration,2018,42(3):545-554.