[发明专利]一种基于分频AVO反演的储层含气性半定量识别方法

摘要

本发明提供了一种基于分频AVO反演的储层含气性半定量识别方法，属于地震勘探领域。本方法包括(1)利用三维地震数据，分别计算纵波速度衰减梯度和横波速度衰减梯度；分别绘制纵波速度衰减梯度与生产井产能之间的关系图和横波速度衰减梯度与生产井产能之间的关系图，开展线性拟合，计算两个衰减梯度与产能之间的相关系数R2；纵波速度衰减梯度与横波速度衰减梯度之间的权重值x的初始值选择为1，即x0=1；(2)沿各井提取纵波速度衰减梯度和横波速度衰减梯度，计算各井生产层段反演值的加权和gi,0,其中i表示井号；(3)以井产能Y为纵坐标、以G为横坐标，进行最小二乘线性拟合。

主权项

一种基于分频AVO反演的储层含气性半定量识别方法，其特征在于：所述方法包括：(1)利用三维地震数据，分别计算纵波速度衰减梯度和横波速度衰减梯度；分别绘制纵波速度衰减梯度与生产井产能之间的关系图和横波速度衰减梯度与生产井产能之间的关系图，开展线性拟合，计算两个衰减梯度与产能之间的相关系数R2；纵波速度衰减梯度与横波速度衰减梯度之间的权重值x的初始值选择为1,即x0＝1；利用地球物理反演技术对公式(1)开展反演，求得纵波速度衰减梯度和横波速度衰减梯度：R(θ,f)≈A(θ)ΔVpVp(f0)+(f-f0)×A(θ)ddf(ΔVpVp)(f)+B(θ)ΔVsVs(f0)+(f-f0)B(θ)ddf(ΔVsVs)(f)---(1)]]>其中，R(θ,f)表示地震反射系数；θ,f分别是地震反射系数的入射角与频率；A(θ),B(θ)分别是只与入射角相关的系数函数；Vp,Vs,ΔVp,ΔVs分别表示地震纵波速度、横波速度、地震纵波速度变化、横波速度变化；分别表示是频率f的函数；(2)沿各井提取纵波速度衰减梯度和横波速度衰减梯度，计算各井生产层段反演值的加权和gi,0，其中i表示井号；所述步骤(2)中是利用公式(5)计算各井生产层段反演值的加权和gi,0：gi,0=aix0+bi,ai=ddf(ΔVpVp)(f),bi=ddf(ΔVsVs)(f)---(5),]]>每改变一次xk，利用公式(5)得到一个对应的gi,k；(3)以井产能Y为纵坐标、以G为横坐标，进行最小二乘线性拟合，形成如下拟合公式：Y＝CG+E   (2)其中C、E为常数；G是gi,0的大写表示；(4)计算各井点产能与步骤(3)得到的拟合线之间的距离和d0；所述步骤(4)所用的计算公式为：d0=Σi=1l|Cgi-yi+E|C2+1=Σi=1l|C(aix0+bi)-yi+E|C2+1---(3)]]>其中d0为距离和、yi为第i井的产能，C、E为最小二乘线性拟合常数；每改变一次xk，利用公式(3)得到一个对应的dk；(5)令xk+1＝xk+t，其中t为步长，k表示迭代次数；(6)重复步骤(2)至步骤(5)，直到dk>dk‑1，转入步骤(7)；(7)输出gi,k，gi,k＝aixk+bi (4)；根据反演结果数值进行含气性识别，反演结果数值越高则含气量越大。