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

权利要求书

1.一种基于分频AVO反演的储层含气性半定量识别方法，其特征在于：所述方法包括：

(1)利用三维地震数据，分别计算纵波速度衰减梯度和横波速度衰减梯度；分别绘制纵波速度衰减梯度与生产井产能之间的关系图和横波速度衰减梯度与生产井产能之间的关系图，开展线性拟合，计算两个衰减梯度与产能之间的相关系数R²；纵波速度衰减梯度与横波速度衰减梯度之间的权重值x的初始值选择为1,即x₀＝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(θ)分别是只与入射角相关的系数函数；V_p,V_s,ΔV_p,ΔV_s分别表示地震纵波速度、横波速度、地震纵波速度变化、横波速度变化；分别表示是频率f的函数；

(2)沿各井提取纵波速度衰减梯度和横波速度衰减梯度，计算各井生产层段反演值的加权和g_i,0，其中i表示井号；所述步骤(2)中是利用公式(5)计算各井生产层段反演值的加权和g_i,0：

gi,0=aix0+bi,ai=ddf(ΔVpVp)(f),bi=ddf(ΔVsVs)(f)---(5),]]>

每改变一次x_k，利用公式(5)得到一个对应的g_i,k；

(3)以井产能Y为纵坐标、以G为横坐标，进行最小二乘线性拟合，形成如下拟合公式：

Y＝CG+E (2)

其中C、E为常数；G是g_i,0的大写表示；

(4)计算各井点产能与步骤(3)得到的拟合线之间的距离和d₀；

所述步骤(4)所用的计算公式为：

d0=Σi=1l|Cgi-yi+E|C2+1=Σi=1l|C(aix0+bi)-yi+E|C2+1---(3)]]>

其中d₀为距离和、y_i为第i井的产能，C、E为最小二乘线性拟合常数；每改变一次x_k，利用公式(3)得到一个对应的d_k；

(5)令x_k+1＝x_k+t，其中t为步长，k表示迭代次数；

(6)重复步骤(2)至步骤(5)，直到d_k>d_k-1，转入步骤(7)；

(7)输出g_i,k，g_i,k＝a_ix_k+b_i (4)；

根据反演结果数值进行含气性识别，反演结果数值越高则含气量越大。