基于贝叶斯理论的AVO三参数波形反演
|
摘要
在实际的AVO反演问题中,叠前数据体中的噪声或其他因素严重影响了AVO反演问题的适定性,而采用先验地质信息作为AVO反演问题的约束条件是解决AVO反演问题不适定的一种可行方法.文中的似然函数采用了ιp范数的解,并用Cauchy分布表示先验模型参数的分布.以此为基础,在反演中建立了测井数据的参数协方差矩阵对反演过程进行约束,并采用了共轭梯度算法实现多参数非线性的反演过程.同时,为了提高反演精度,避免动校正拉伸及依赖于炮检距的调谐效应对参数估计的影响,反演采用动校前地震数据进行参数估计.从应用效果分析来看,即使叠前道集的信噪比不高,反演的结果也能较好地与实际情况相匹配,为识别储层流体性质提供了新的手段.
In the actual AVO inversion, the noise and other factors heavily influence the well-posed property of the inversion problem. So AVO inversion constrained by a priori geological information is a practical method to the solution of ill-posed inversion problem. In this paper the solution for the ιp norm is used to describe the likelihood function. The probability of the prior model parameter is independently Cauchy distribution. Based on these theories the impedance reflectivity attribute parameter covariance matrix from well log is established to constrain inversion, and the conjugate gradient algorithm is used to implement multi-parameter non-linear inversion. In order to improve the inversion accuracy and avoid NMO stretch and offset dependent tuning which influences the estimation of the attribute parameters, the attribution parameter inversion is based on synthetic gather before NMO. Tests on synthetic data show that all inverted parameters were almost perfectly retrieved compared to the actual values, even if the signal/noise ratio is not higher. Therefore it can provide a new tool for identifying fluid content in reservoir pores.
In the actual AVO inversion, the noise and other factors heavily influence the well-posed property of the inversion problem. So AVO inversion constrained by a priori geological information is a practical method to the solution of ill-posed inversion problem. In this paper the solution for the ιp norm is used to describe the likelihood function. The probability of the prior model parameter is independently Cauchy distribution. Based on these theories the impedance reflectivity attribute parameter covariance matrix from well log is established to constrain inversion, and the conjugate gradient algorithm is used to implement multi-parameter non-linear inversion. In order to improve the inversion accuracy and avoid NMO stretch and offset dependent tuning which influences the estimation of the attribute parameters, the attribution parameter inversion is based on synthetic gather before NMO. Tests on synthetic data show that all inverted parameters were almost perfectly retrieved compared to the actual values, even if the signal/noise ratio is not higher. Therefore it can provide a new tool for identifying fluid content in reservoir pores.
引文
[1]Simmons J L,Backus M M.Waveform-based AVO inversion and AVOprediction error.Geophysics,1996,61(6):1575～1588
[2]Buland A,Omre H.Bayesianlinearized AVOinversion.Geophysics,2003,68(1):185～198
[3]牛滨华,文鹏飞.基于BSR的AVO正演估算水合物含量方法的研究.地球物理学报,2006,49(1):143～152Niu B H,Wen P F.Estimation of hydrate concentration based on AVO modeling of BSR.Chinese J.Geophys.(in Chinese),2006,49(1):143～152
[4]李景叶,陈小宏,郝振江.多波时移地震AVO反演研究.地球物理学报,2005,48(4):902～908Li J Y,Chen X H,Hao ZJ.Astudy on multiple time-lapse seismic AVOinversion.Chinese J.Geophys.(in Chinese),2005,48(4):902～908
[5]Debski W,Tarantola A.Information on elastic parameters obtained fromthe amplitudes of reflected waves.Geophysics,1995,60(5):1426～1436
[6]Westergren R.Mathematics Handbook for Science and Engineering.Springer-Verlag,1999
[7]Sacchi MD,Ulrych TD.High-resolution velocity gathers and offset space reconstruction.Geophysics,1995,60(6):1169～1177
[8]Dong W.Fluid line distortion due to migration stretch.66th Ann Internat Mtg:Soc of Expl.Geophys.,1996.1345～1348
[9]Dong W.AVO detectability against tuning and stretching artifacts.68thAnn Internat Mtg:Soc of Expl.Geophys.,1998.236～239
[10]Shuey R T.Asimplification of the Zoeppritz equations.Geophysics,1985,50(4):609～614
[11]Claerbout J F.Earth Soundings Analysis:Processing Versus Inversion.Blackwell Scientific Publications,1992
[12]DowntonJ,Lines L.High-resolution AVOanalysis before NMO.73rd Ann Internat Mtg:Soc of Expl.Geophys.,2003.219～222
[13]GidlowP M,Smith GC,Vail PJ.Hydrocarbon detection usingfluid factor traces,a case study:Howuseful is AVOanalysis?Joint SEG/EAEG Summer Research Workshop,Technical Program and Abstracts,1992.78～89
[14]Aki K I,Richards P G.Quantitative Seismology:Theory and Methods.Cambridge:WHFreeman and Co,1980.144～154
[15]Tarantala A.Inversion Problem Theory and Methods for Model Parameter Estimation.Societyfor Industry and Applied Mathematics,Philadelphia,2004
[16]Gardner L W,Gregory A R.Formation velocity and density-The diagnostic basics for stratigraphic traps.Geophysics,1974,39(6):770～780
[17]Smith G C,Gidlow P M.Weighted stacking for rock property estimation and detection of gas.Geophysical Prospecting,1987,35:993～1014
[18]Greenberg M L,Castagna J P.Shear-wave velocity estimation in porous rocks:theoretical formulation preliminary verification and applications.Geophys.Prosp.,1992,40(2):195～210
[2]Buland A,Omre H.Bayesianlinearized AVOinversion.Geophysics,2003,68(1):185～198
[3]牛滨华,文鹏飞.基于BSR的AVO正演估算水合物含量方法的研究.地球物理学报,2006,49(1):143～152Niu B H,Wen P F.Estimation of hydrate concentration based on AVO modeling of BSR.Chinese J.Geophys.(in Chinese),2006,49(1):143～152
[4]李景叶,陈小宏,郝振江.多波时移地震AVO反演研究.地球物理学报,2005,48(4):902～908Li J Y,Chen X H,Hao ZJ.Astudy on multiple time-lapse seismic AVOinversion.Chinese J.Geophys.(in Chinese),2005,48(4):902～908
[5]Debski W,Tarantola A.Information on elastic parameters obtained fromthe amplitudes of reflected waves.Geophysics,1995,60(5):1426～1436
[6]Westergren R.Mathematics Handbook for Science and Engineering.Springer-Verlag,1999
[7]Sacchi MD,Ulrych TD.High-resolution velocity gathers and offset space reconstruction.Geophysics,1995,60(6):1169～1177
[8]Dong W.Fluid line distortion due to migration stretch.66th Ann Internat Mtg:Soc of Expl.Geophys.,1996.1345～1348
[9]Dong W.AVO detectability against tuning and stretching artifacts.68thAnn Internat Mtg:Soc of Expl.Geophys.,1998.236～239
[10]Shuey R T.Asimplification of the Zoeppritz equations.Geophysics,1985,50(4):609～614
[11]Claerbout J F.Earth Soundings Analysis:Processing Versus Inversion.Blackwell Scientific Publications,1992
[12]DowntonJ,Lines L.High-resolution AVOanalysis before NMO.73rd Ann Internat Mtg:Soc of Expl.Geophys.,2003.219～222
[13]GidlowP M,Smith GC,Vail PJ.Hydrocarbon detection usingfluid factor traces,a case study:Howuseful is AVOanalysis?Joint SEG/EAEG Summer Research Workshop,Technical Program and Abstracts,1992.78～89
[14]Aki K I,Richards P G.Quantitative Seismology:Theory and Methods.Cambridge:WHFreeman and Co,1980.144～154
[15]Tarantala A.Inversion Problem Theory and Methods for Model Parameter Estimation.Societyfor Industry and Applied Mathematics,Philadelphia,2004
[16]Gardner L W,Gregory A R.Formation velocity and density-The diagnostic basics for stratigraphic traps.Geophysics,1974,39(6):770～780
[17]Smith G C,Gidlow P M.Weighted stacking for rock property estimation and detection of gas.Geophysical Prospecting,1987,35:993～1014
[18]Greenberg M L,Castagna J P.Shear-wave velocity estimation in porous rocks:theoretical formulation preliminary verification and applications.Geophys.Prosp.,1992,40(2):195～210
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心