基于Metropolis抽样的弹性阻抗随机反演

笔者提出了基于Metropolis抽样的弹性阻抗随机反演方法,该方法是一种基于蒙特卡洛的非线性反演方法,能够有效地融合测井资料中的高频信息,提高反演结果的分辨率。应用贝叶斯理论框架,首先通过快速傅里叶滑动平均模拟算法(fast Fourier transform-moving average,FFT-MA)和逐渐变形算法(gradual deformation method,GDM)得到基于地质统计学的先验信息,然后用Metropolis抽样算法对后验概率密度进行抽样,得到反演问题的解。其中FFT-MA模拟作为一种高效的频率域模拟方法,结合GDM后,在保持模拟空间结构不变的前提下,可以连续修改储层模型,直至满足实际观测地震记录。数值试验表明:FFT-MA模拟有效地提高了计算效率,融入GDM更新算法,可以保证反演结果有效地收敛,并且反演结果与理论模型吻合较好,具有较高的分辨率,该方法采用两步法反演弹性参数,运算速度较快。     

基于分形先验信息的非线性反演

基于分形先验信息的非线性反演方法能综合利用测井数据和地震信息,在贝叶斯框架下,通过分形高斯噪音算法得到基于分形理论的先验信息,然后根据地震资料构建似然函数,最终利用基于快速模拟退火算法(Very Fast Simulated Annealing,VFSA)改进的粒子群优化(Particle Swarm Optimization,PSO)算法(VFSA-PSO)实现后验概率密度的抽样。与确定性反演结果相比,该方法能够有效地融合测井资料中的高频信息,提高反演结果的分辨率,并且目标函数的建立融合了确定性反演中的低频约束,从而得到宽频带的反演结果。数值模拟试验表明:基于分形先验的非线性反演结果与理论模型吻合较好,实际资料的应用效果也证明了该反演方法的有效性。     

基于保幅拉东变换的多次波衰减

为在去除多次波时有效保护地震一次反射波数据的AVO现象,给后续反演、解释提供准确的地震数据,本文提出了一种基于保幅拉东变换的多次波衰减方法,该方法是对常规抛物拉东变换的修改,把常规的稀疏拉东变换在拉东域分成两部分：一部分用于模拟零偏移距处的反射波能量,增加的另一部分用于模拟反射波振幅的AVO特性.该方法不仅考虑了反射波同相轴的形状,还考虑了反射波同相轴振幅幅度的变化,从而可把反射波信息进行有效转换,进而有利于多次波的消除,更好地恢复有效波的能量.在把地震数据由时间域转换到拉东域时,本文采用了IRLS算法实现保幅拉东算子的反演.模型数据和实际地震道集的试算分析表明,与常规拉东变换相比,保幅拉东变换在去除多次波的同时可有效保护一次反射波的AVO现象.     

各向异性碳酸盐岩储层精细横波速度估算方法

横波速度测井曲线是叠前地震反演必不可少的资料,在储层预测中起着非常重要的作用,但在研究中比较缺乏.考虑到地下储层的复杂多样性,为了准确求取地下储层的横波速度曲线,通过对碳酸盐岩储层物性的分析及其裂缝对储层的各向异性分析,根据其特征,首次建立了各向异性碳酸盐岩储层的岩石物理理论模型,并依据此模型对各向异性碳酸盐岩储层的横波速度进行了估算.在此基础上,引入反演的思想,利用测井约束反演及叠前地震道集的弹性波形反演等方法,对基于岩石物理理论模型的初步估算横波速度进行了修正,确保了估算结果的准确性,得到了与地下储层更加吻合的精细横波速度;通过实际资料试算,对估算横波与实测横波进行了对比分析.结果表明,该估算方法可成功模拟各向异性碳酸盐岩储层的横波速度,并在实际工区得到很好的应用.