基于稀疏约束贝叶斯估计的相对波阻抗反演
|
摘要
地震道积分是一种利用地震资料进行的无约束反演技术,可以方便地得到地层的相对波阻抗,但采用稀疏反演的反射系数递推反演相对波阻抗,横向连续性差、可用性低。给出了一种基于稀疏约束贝叶斯估计的地震相对波阻抗反演算法,即在反射系数Cauchy概率分布稀疏约束下,基于贝叶斯最大后验概率估计,采用预条件共轭梯度法估计地震子波和反射系数,从而得到高分辨率的反褶积结果,进而得到高分辨率的相对波阻抗剖面。利用理论模型和实际数据对算法进行了验证,结果表明:基于稀疏约束贝叶斯估计的相对波阻抗反演方法可行;与直接法稀疏反演相比,预条件共轭梯度法稀疏反演精度高、收敛快、数值计算稳定。
Seismic trace integration is a non-constraint inversion technique. However,the relative acoustics impedance recursed by reflectivity obtained from sparse inversion has bad lateral continuity and worse adaptability.Based on Cauchy sparseness constraint Bayesian estimation,an inversion algorithm for seismic relative acoustic impedance was proposed.Assuming the Cauchy sparseness distribution of reflectivity,based on Bayesian maximum posterior probability estimation,pre-condition conjugate gradient inversion algorithm was adopted to realize the simultaneous estimation of reflectivity and wavelet.Then high-resolution deconvolution result was obtained,and high-resolution relative acoustic impedance profile was achieved.The inversion algorithm was tested by both theoretical model and real 2-D seismic data.Test results indicate that, relative acoustic impedance inversion based on Cauchy sparseness constraint Bayesian estimation is feasible.Compared to direct sparseness inversion algorithm,pre-condition conjugate gradient sparseness inversion algorithm is stable,with higher accuracy and faster converging speed.
Seismic trace integration is a non-constraint inversion technique. However,the relative acoustics impedance recursed by reflectivity obtained from sparse inversion has bad lateral continuity and worse adaptability.Based on Cauchy sparseness constraint Bayesian estimation,an inversion algorithm for seismic relative acoustic impedance was proposed.Assuming the Cauchy sparseness distribution of reflectivity,based on Bayesian maximum posterior probability estimation,pre-condition conjugate gradient inversion algorithm was adopted to realize the simultaneous estimation of reflectivity and wavelet.Then high-resolution deconvolution result was obtained,and high-resolution relative acoustic impedance profile was achieved.The inversion algorithm was tested by both theoretical model and real 2-D seismic data.Test results indicate that, relative acoustic impedance inversion based on Cauchy sparseness constraint Bayesian estimation is feasible.Compared to direct sparseness inversion algorithm,pre-condition conjugate gradient sparseness inversion algorithm is stable,with higher accuracy and faster converging speed.
引文
[1]张永刚.地震波阻抗反演技术的现状和发展[J].石油物探,2002,41(4):385-389
[2]黄文松,年静波,刘喜武.地震波阻抗反演的预条件共轭梯度法[J].新疆石油地质,2006,27(4): 425-428
[3]刘胜,邱其祥,程增庆,等.利用相对波阻抗识别煤层缺失变薄带[J].物探与化探,2008,32(6):682-684
[4]姚建阳.用地震道积分方法提高地层的识别能力[J].石油物探,1990,29(1):40-49
[5]崔炯成,吴媚,李维新,等.高分辨率非线性地震波阻抗反演方法和应用[J].石油物探,2007,46(1):1-12
[6]王振国,陈小宏,王克斌,等.多参数优化的模拟退火法波阻抗反演[J].石油物探,2007,46(2):120-125
[7]陈双全,王尚旭,季敏,等.地震波阻抗反演的蚁群算法实现[J].石油物探,2005,44(6):551-553
[8]李晶,陈裕明,唐湘蓉.用退火遗传算法进行地震波阻抗反演[J].石油物探,2004,43(3):234-237
[9]刘喜武,刘洪.地震盲反褶积综述[J].地球物理学进展,2003,18(2):203 209
[10]张繁昌,刘杰,印兴耀,等.修正柯西约束地震盲反褶积方法[J].石油地球物理勘探,2008,43(4): 391-396
[11]刘喜武,宁俊瑞,张改兰.Cauchy稀疏约束Bayesian估计地震盲反褶积框架与算法研究[J].石油物探, 2009,48(5):459-464
[12]Sacchi M D.Statistical and transform method in geophysical signal processing[D].Edmonton,Canada: Department of Physics,University of Alberta,2002
[13]刘喜武,刘洪,李幼铭.反射系数和子波同时迭代反演的预条件共轭梯度法[J].物探化探计算技术, 2006,21(3):211-215
[14]孟小红,吴何珍,刘国峰.盲源反褶积方法与应用研究[J].石油地球物理勘探,2005,40(6):642-645
[15]Claerbout J F.Fundamentals of geophysical data processing-with application to petroleum prospecting [M].New York:Blackwell Scientific Publications, 1985:130-143
[16]刘喜武,刘洪.实现稀疏反褶积的预条件共轭梯度法[J].物探化探计算技术,2003,25(3):215-219
[17]朱海龙译,Cheng Q著.对反射地震信号同时进行子波估计和反褶积[J].石油物探译丛,1996(6):41-50
[18]赵宪生,梁桂容,罗运先.子波反演与波阻抗反演[J].物探化探计算技术,1999,21(3):207-211
[19]牟永光,陈小宏,李国发,等.地震数据处理方法[M].北京:石油工业出版社,2007:227-232
[2]黄文松,年静波,刘喜武.地震波阻抗反演的预条件共轭梯度法[J].新疆石油地质,2006,27(4): 425-428
[3]刘胜,邱其祥,程增庆,等.利用相对波阻抗识别煤层缺失变薄带[J].物探与化探,2008,32(6):682-684
[4]姚建阳.用地震道积分方法提高地层的识别能力[J].石油物探,1990,29(1):40-49
[5]崔炯成,吴媚,李维新,等.高分辨率非线性地震波阻抗反演方法和应用[J].石油物探,2007,46(1):1-12
[6]王振国,陈小宏,王克斌,等.多参数优化的模拟退火法波阻抗反演[J].石油物探,2007,46(2):120-125
[7]陈双全,王尚旭,季敏,等.地震波阻抗反演的蚁群算法实现[J].石油物探,2005,44(6):551-553
[8]李晶,陈裕明,唐湘蓉.用退火遗传算法进行地震波阻抗反演[J].石油物探,2004,43(3):234-237
[9]刘喜武,刘洪.地震盲反褶积综述[J].地球物理学进展,2003,18(2):203 209
[10]张繁昌,刘杰,印兴耀,等.修正柯西约束地震盲反褶积方法[J].石油地球物理勘探,2008,43(4): 391-396
[11]刘喜武,宁俊瑞,张改兰.Cauchy稀疏约束Bayesian估计地震盲反褶积框架与算法研究[J].石油物探, 2009,48(5):459-464
[12]Sacchi M D.Statistical and transform method in geophysical signal processing[D].Edmonton,Canada: Department of Physics,University of Alberta,2002
[13]刘喜武,刘洪,李幼铭.反射系数和子波同时迭代反演的预条件共轭梯度法[J].物探化探计算技术, 2006,21(3):211-215
[14]孟小红,吴何珍,刘国峰.盲源反褶积方法与应用研究[J].石油地球物理勘探,2005,40(6):642-645
[15]Claerbout J F.Fundamentals of geophysical data processing-with application to petroleum prospecting [M].New York:Blackwell Scientific Publications, 1985:130-143
[16]刘喜武,刘洪.实现稀疏反褶积的预条件共轭梯度法[J].物探化探计算技术,2003,25(3):215-219
[17]朱海龙译,Cheng Q著.对反射地震信号同时进行子波估计和反褶积[J].石油物探译丛,1996(6):41-50
[18]赵宪生,梁桂容,罗运先.子波反演与波阻抗反演[J].物探化探计算技术,1999,21(3):207-211
[19]牟永光,陈小宏,李国发,等.地震数据处理方法[M].北京:石油工业出版社,2007:227-232
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心