[发明专利]基于自适应基质矿物等效弹性模量反演的横波速度估算方法

权利要求书

1.基于自适应基质矿物等效弹性模量反演的横波速度估算方法，其特征在于采用以下步骤：

(1)以单个测井深度点为研究对象，根据该深度点的测井资料(纵波速度、密度和孔隙度)，采用砂泥岩的纵波速度与横波速度的经验关系式(公式1)估算初始横波速度，设定该测井深度点的干岩石的泊松比的取值范围(一般介于0.0～0.4)、初始值(一般等于取值范围的下限值或上限值)和增量(一般介于0.01～0.05之间)，并计算基质矿物的等效体积模量的取值范围(公式2)、初始值(一般等于取值范围的下限值或上限值)和增量(一般介于0.0～1.0GPa之间)；

V_s＝0.804V_p-0.856    (1)

下限值：ρ(Vp2-43Vs2)/(1-φ)[3/(1-φ)],]]>上限值：[ρ(Vp2-43Vs2)+aφ]/(1-φ)[3/(1-φ)]---(2)]]>在公式1和公式2中，V_p、V_s和ρ代表流体饱和岩石的纵波速度、横波速度和密度，φ代表孔隙度，a为常数(一般介于100～400)；

(2)分别采用两种不同的流体项计算方法：基于Gassmann方程(公式3)和Gassmann-Boit-Geertsman方程(公式4)的方法和基于Russell等人2003年提出的流体识别因子(公式5)的方法，同时计算流体项，并设定上述两种流体项之差的绝对值为反演目标函数，采用寻找全局最优解的方法反演最优的干岩石的泊松比和基质矿物的等效体积模量，当目标函数的值＜设定的阀值(根据计算精度和计算效率选择，取值一般介于0～0.5之间)，则停止计算，输出当前深度点的岩石基质矿物的等效体积模量和干岩石的泊松比；

f=β2M=(1-KdryK0)2φKfl+1-φK0-KdryK02,]]>β=1-KdryK0,]]>1M=β-φK0+φKfl---(3)]]>

(Y-1)β2+[Yφ(K0Kfl-1)-Y+NK0]β-φ(Y-NK0)(K0Kfl-1)=0,]]>Y=3(1-σdry)(1+σdry),]]>N=ρVP2---(4)]]>

ρf=(Zp2-cZs2),]]>c=(VpVs)dry2=2(1-σdry)1-2σdry,]]>Z_p＝ρV_p，Z_s＝ρV_s    (5)

在公式3、公式4和公式5中，f代表流体项，V_p、V_s和ρ代表流体饱和岩石的纵波速度、横波速度和密度，K_dy和σ_dry分别代表干岩石的体积模量和泊松比，K₀代表基质矿物的等效体积模量，K_fl代表流体的体积模量，φ代表孔隙度；

(3)基于步骤2反演得到的干岩石的泊松比和基质矿物的等效体积模量，采用Nur(1995)基于临界孔隙度模型的线性平均公式(公式6)和Gassmann方程(公式7)计算基质矿物的等效剪切模量(μ₀)；

μ0=μdry(1-φφc),]]>μdry=43Kdry(S-1),]]>S=3(1-σdry)1+σdry---(6)]]>

Kdry=Ksat(φK0Kfl+1-φ)-K0φK0Kfl+KsatK0-1-φ,]]>Ksat=ρ(Vp2-43Vs2)---(7)]]>

在公式6中，μ_dry和σ_dry代表干岩石的剪切模量和泊松比，φ_c代表临界孔隙度(一般取值0.40)；

(4)基于步骤2和步骤3得到的基质矿物的等效体积模量和等效剪切模量，采用简化的Xu-White模型(Keys和Xu，2002)估算横波速度，计算公式见公式8和公式9；

Vs=μd/ρ,]]>μ_d＝μ₀(1-φ)^q，q=15Σl=s,shυlFiijj(αl)---(8)]]>

Fiijj(α)=2F3+1F4+F4F5+F6F7-F8F9F2F4---(9-1)]]>

F2=1+A[1+32(g+γ)-R2(3g+5γ)]+B(3-4R)+A2(A+3B)(3-4R)[g+γ-R(g-γ+2υ2)]---(9-2)]]>

F3=1+A2[R(2-γ)-1+α2α2g(R-1)]---(9-3)]]>

F4=1+A4[3γ+g-R(g-γ)]---(9-4)]]>

F5=A[R(g+γ-43)-g]+Bγ(3-4R)---(9-5)]]>

F₆＝1+A[1+g-R(γ+g)]+B(1-γ)(3-4R)    (9-6)

F7=2+A4[9γ+3g-R(5γ+3g)]+Bγ(3-4R)---(9-7)]]>

F8=A[1-2R+g2(R-1)+υ2(5R-3)]+B(1-γ)(3-4R)---(9-8)]]>

F₉＝A[g(R-1)-Rγ]+Bγ(3-4R)           (9-9)

A＝-1，B＝0，R=3μ03K0+4μ0,]]>g=α21-α2(3γ-2),]]>γ=α(1-α2)3/2[cos-1(α)-α1-α2]---(9-10)]]>

在公式8和公式9中，υ_l和α_l分别代表砂岩和泥岩占岩石基质的体积百分比(根据测井泥质含量计算)和孔隙扁率(砂岩孔隙的扁率取值0.10，泥岩孔隙的扁率取值0.04，与常规的Xu-White模型的取值一致)，F_iiij(α_l)代表孔隙扁率为α_l时的Eshelby张量；

(5)选取下一个测井深度点，重复步骤1～5的计算，直至目标深度段内的所有测井深度点计算完成，则停止计算；

(6)输出该目标深度段内每个测井深度点估算的横波速度值。