基于Hampel三截尾函数的储层弹性和物性参数同步反演  

Simultaneous inversion on elastic and physical properties of reservoir based on Hampel's three-part redescending function

作　　者：李志勇[1] 张家树[2] 蔡涵鹏[1] 胡光岷[1] 

机构地区：[1]电子科技大学通信与信息工程学院,四川成都611731 [2]西南交通大学信号与信息处理四川省重点实验室,四川成都610031

出　　处：《石油物探》2017年第2期261-272,共12页Geophysical Prospecting For Petroleum

基　　金：国家自然科学基金联合基金项目(U1562218)资助~~

摘　　要：常规的地震反演方法通常假设噪声信号服从高斯分布,以此构建出的目标函数在特定地区反演效果好,但其普适性较小。假设地震噪声信号和岩石物理模型误差服从高斯分布且带有加性脉冲噪声,引入Hampel三截尾函数构造能够同时压制高斯噪声和脉冲噪声的反演目标函数,利用牛顿法求解目标函数得到储层弹性和物性参数反演的迭代公式,并通过阈值的选择适应噪声分布不同的地区。该算法利用Hampel三截尾函数构建加权矩阵,能在一定程度上自适应调节地震数据、岩石物理约束和先验信息之间的权重,消除地震和岩石物理噪声的统计性误差造成的反演不稳定性,最终得到稳定的储层弹性和物性参数反演结果。通过理论模型测试和实际资料应用显示了算法的稳定性和可靠性。The discrepancies between geophysical measurements and forward modeling data are commonly modeled as Gaussian er-rors by conventional seismic inversion methods. However, misfit functions constructed by this distribution always have good inver-sion result in specified regions but have weak adaptability. We assume Gaussian distribution with additional impulse distribution for both seismic noise data and elastic property deviations from the rock-physics model.In this case, Hampel’s threepart redescending estimate function is included in the misfit function, which could perform better to suppress both Gaussian and impulse errors as a robust measure.Minimizing the misfit function classically leads to the Newton？s optimization method, with which elastic and phys-ical properties of reservoir are obtained.This approach can be applied to the areas with different distribution of errors by choosing the degree of threshold.The relative contributions of the rock-physics constraints and prior information are adaptively adjusted, which can increase the stability of the algorithm and ensure a better inversion Both synthetic and real seismic data examples show that the method can obtain stable and reliable inversion results.

关 键 词：AVO反演 岩石物理反演 联合反演 高斯噪声 脉冲噪声 Hampel三截尾函数 

分 类 号：P631[天文地球—地质矿产勘探]