Cauchy稀疏约束Bayesian估计地震盲反褶积框架与算法研究

Cauchy稀疏约束Bayesian估计地震盲反褶积框架与算法研究

以二阶统计学方法为基础,从Canadas等提出的非最小相位子波和非白噪反射系数地震盲反褶积框架出发,给出了Cauchy稀疏约束Bayesian估计地震盲反褶积框架.基于反射系数与子波相互独立(或弱相关)的假设,分别构建了反射系数和子波最优估计方程,并采用预条件共轭梯度法迭代反演实现反射系数和子波的同时估计.在方法具体实现时,以传统脉冲反褶积结果作为迭代初值,通过迭代得到反射系数和任意相位子波;然后再对子波进行最小相位化,通过反演得到反子波;最后将反子波与地震道褶积,得到反褶积结果.利用理论模型和实际数据对算法进行了试算,结果表明,给出的地震盲反褶积理论框架是正确的;与直接(共轭梯度求解正则方程)稀疏同时迭代反演法的对比显示,预条件共轭梯度算法稳定,精度高,收敛快.

作 者： 刘喜武 宁俊瑞 张改兰 Liu Xiwu Ning Junrui Zhang Gailan   作者单位： 刘喜武,Liu Xiwu(中国石油化工股份有限公司石油勘探开发研究院,北京100083;中国石油大学(北京)资源与信息学院,北京102249;中国海洋大学海洋地球科学学院,山东青岛266100)

宁俊瑞,张改兰,Ning Junrui,Zhang Gailan(中国石油化工股份有限公司石油勘探开发研究院,北京,100083) 

刊 名： 石油物探  ISTIC PKU 英文刊名： GEOPHYSICAL PROSPECTING FOR PETROLEUM  年，卷(期)： 2009 48(5)  分类号： P631.4  关键词： 地震盲反褶积   子波估计   Cauchy稀疏约束   Bayesian估计   预条件共轭梯度法   seismic blind deconvolution   wavelet estimation   Cauchy sparse constraint   Bayesian estimation   pre-conditional conjugate gradient algorithm