[发明专利]一种致密砂岩储层岩石弹性模量计算方法

权利要求书

1.一种致密砂岩储层岩石弹性模量计算方法，其特征在于，包括：

获取测井参数和致密砂岩储层的岩石物性参数；

根据所述测井参数和所述岩石物性参数，结合预先建立的孔隙裂缝衰减模型计算高频下的岩石弹性模量；

根据所述高频下的岩石弹性模量，结合预先建立的自洽模型，计算岩石的高频基质弹性模量；

利用所述岩石的高频基质弹性模量，结合所述孔隙裂缝衰减模型，计算任意频率下的饱和岩石弹性模量。

2.根据权利要求1所述的方法，其特征在于，所述岩石的高频基质弹性模量，结合所述孔隙裂缝衰减模型，计算任意频率下的饱和岩石弹性模量，之后还包括：

根据所述任意频率下的饱和岩石弹性模量，结合开始模拟退火寻优的方法，计算岩石的声波速度。

3.根据权利要求1所述的方法，其特征在于，所述根据所述测井参数和所述岩石物性参数，结合预先建立的孔隙裂缝衰减模型计算高频下的岩石弹性模量包括：

根据所述测井参数与岩石物性参数，利用wood公式和patchy模型计算岩石流体体积模量；

根据所述岩石物性参数与测井参数，利用VRH模型计算岩石基质弹性模量；

根据所述岩石流体体积模量和所述岩石基质弹性模量，结合所述孔隙裂缝衰减模型计算高频下的岩石弹性模量；

其中，所述孔隙裂缝衰减模型的计算公式为：

K1=K-ϵ{4(3λ+2μ)(λ+2μ)μ(λ+μ)[1-3A(ω)]-4πrA(ω)}-φP{3λ+2μ4μ[λ+2μ3λ+2μ+B(ω)]-3B(ω)}]]>和μ1=μ-φP15μ(λ+2μ)9λ+14μ-16ϵμ(λ+2μ)45(1+KC)(3λ+μ)(KC+11+iωτ)-3215ϵμ(λ+2μ)3λ+4μ;]]>

其中，K、λ和μ分别代表岩石基质体积模量、拉梅常数和岩石基质剪切模量，ε为裂缝密度，φ_P为刚性孔隙度，ω为角频率，r为软孔隙纵横比，K_f为岩石流体体积模量，

A(ω)=(1+iωτ)λ+2μλ+μ[16ϵ27φ(1+Kp)+λ+μ3λ+2μ]+iωτ[13(1+Kc)-γ′]1+iωτ+(1+iωγτ)16ϵ(1+Kc)(λ+2μ)9φ(1+Kp)(λ+μ)]]>

KP=4μ3Kf,γ=3π(λ+μ)(1+KP)4(λ+2μ)(1+KC),γ′=γλ+2μ(3λ+2μ)(1+KP),]]>

B(ω)=13(1+Kc)+9(1+Kp)16λ+μ(3λ+2μ)(1+Kc)+iωγτ1+iωτ[γ′-13(1+Kc)]9(1+Kp)(λ+μ)16(1+Kc)(λ+2μ)+1+iωγτ1+iωτ]]>

τ为时间常数且η为流体粘滞性，α为柔性孔隙半径，υ为岩石基质泊松比，κ为渗透率，ζ为特征矿物颗粒大小，φ为总孔隙度，K1为高频下岩石弹性模量的体积模量，μ1为高频下岩石弹性模量的剪切模量。