CN104573150A - 一种针对富含有机质页岩的各向异性岩石物理建模方法 - Google Patents

Abstract

本发明提供了一种针对富含有机质页岩的各向异性岩石物理建模方法，属于油气地球物理勘探领域。本方法包括：第一步，输入组成岩石矿物的比例和模量，计算得到多种矿物混合的岩石基质的弹性模量；第二步，向所述岩石基质中依次加入不同类型的孔隙得到干岩石，计算得到含孔隙干岩石的骨架模量；第三步，向所述干岩石中加入干酪根得到含有机质干岩石，计算得到含干酪根的干岩石刚度矩阵；第四步，向所述含有机质干岩石中加入流体，计算得到整个有机页岩的弹性模量。

Description

一种针对富含有机质页岩的各向异性岩石物理建模方法

技术领域

本发明属于油气地球物理勘探领域，勘探对象为页岩气，具体涉及一种针对富含有机质页岩的各向异性岩石物理建模方法。

背景技术

岩石物理的主要任务在于如何将地下岩石的物性参数与弹性参数(如速度、密度)联系起来。岩石的矿物成分，孔隙度、含流体性甚至孔隙形状都会影响岩石的弹性模量，为了研究和模拟地下复杂的岩石与模量间的相互关系，人们发展了KT理论(Kuster and Toksoz，1974)，自相容模型(0’Connell andBudiansky，1974)，等效微分介质理论(Cleary，1980；Zimmerman，1991a)。一些经典的岩石物理建模方法依赖于建立速度与孔隙度、泥质含量和其他一些物性参数之间的关系，这种方法在砂岩储层和灰岩储层中取得了巨大的成功，如针对砂泥岩的Han模型，Xu-White模型，Gassmann模型以及针对碳酸盐岩的Xu-Panye模型等。

对于致密砂岩储层，往往具有较低的孔隙度和渗透率并发育有较强的裂缝，目前国外已有的技术包括Hudson模型，Eshelby-Cheng裂缝模型。然而对于富含有机质的页岩，这些模型都显得苍白无力。页岩中复杂的岩性(常包括砂岩、石灰石、泥质-硅质混合物)和强非均匀性给建模过程带来了巨大的麻烦。此外，页岩中所富含的有机质以及应力非均匀性还会带来较强的各向异性。因此，在针对有机页岩的岩石物理建模过程中不得不考虑这些因素的影响。目前，我国尚不具备一套完整的针对有机质页岩的各向异性岩石物理建模方法与流程。

发明内容

本发明的目的在于解决上述现有技术中存在的难题，提供一种针对富含有机质页岩的各向异性岩石物理建模方法，针对富含有机质的页岩，用于页岩气勘探领域。意在建立一种富含有机质的各向异性页岩岩石物理建模方法和流程，帮助更好的了解和刻画页岩岩石物性对于弹性参数和力学参数的影响，指导页岩气勘探甜点预测。

本发明是通过以下技术方案实现的：

一种针对富含有机质页岩的各向异性岩石物理建模方法，包括：

第一步，输入组成岩石矿物的比例和模量，计算得到多种矿物混合的岩石基质的弹性模量；

第二步，向所述岩石基质中依次加入不同类型的孔隙得到干岩石，计算得到含孔隙干岩石的骨架模量；

第三步，向所述干岩石中加入干酪根得到含有机质干岩石，计算得到含干酪根的干岩石刚度矩阵；

第四步，向所述含有机质干岩石中加入流体，计算得到整个有机页岩的弹性模量。

所述第一步是这样实现的：

由公式(1)计算多种矿物混合的岩石基质弹性模量：

$M_{\mathrm{VRH}}=\frac{M_v+M_R}{2},---\left(1\right)$

其中

$M_V=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{f}_i{M}_i,\frac{1}{M_R}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}\frac{f_i}{M_i}$

f_i和M_i分别代表第i种矿物的比例和模量；N表示N种矿物。

所述第二步是这样实现的：

基于第一步得到的多种矿物混合的岩石基质的弹性模量，利用各向同性等效微分介质理论向岩石基质中依次加入不同形状的孔隙，计算得到含孔隙干岩石的骨架模量：等效体积模量K^*和剪切模量μ^*；

不同的形状是通过孔隙的宽长比来刻画的；

计算得到含孔隙干岩石的模量是这样实现的：

每加入一种孔隙，计算一次含孔隙干岩石的模量，具体如下：加入第一种孔隙，计算得到第一次含孔隙干岩石的模量，然后在此基础上加入第二种孔隙，计算得到第二次含孔隙干岩石的模量，以此类推，直到加入最后一种孔隙，然后计算得到最后的含孔隙干岩石的模量。

所述第三步是这样实现的：

基于第二步得到的含孔隙干岩石的模量，利用各向异性等效微分介质理论向干岩石中加入干酪根，计算得到含干酪根的干岩石刚度矩阵。

所述第四步是这样实现的：

利用各向异性Gassmann方程向含有机质干岩石中加入流体，计算得到整个有机页岩的弹性模量。

与现有技术相比，本发明的有益效果是：本发明利用介质平均理论、微分等效介质理论、各向异性等效介质理论及各向异性流体替代方程成功建立了页岩岩石物理建模流程，考虑了不同矿物含量、孔隙度、孔隙形状及干酪根有机质对于页岩岩石弹性模量的影响。可用于分析不同的岩石物性参数对弹性参数及地震响应的影响，有效指导和支持页岩气勘探的甜点圈闭工作。

附图说明

图1有机质页岩岩石物理建模流程示意图

图2-1预测速度随干酪根含量的变化图

图2-2各向异性参数ε随干酪根含量的变化图

图2-3各向异性参数γ随干酪根含量的变化图

图2-4各向异性参数δ随干酪根含量的变化图

图3-1岩石杨氏模量随粘土含量和TOC的变化

图3-2泊松比随粘土含量和TOC的变化

图4-1利用测井及岩芯资料来进行的岩石物理标定中的伽马值

图4-2利用测井及岩芯资料来进行的岩石物理标定中的井径

图4-3利用测井及岩芯资料来进行的岩石物理标定中的纵波速度

图4-4利用测井及岩芯资料来进行的岩石物理标定中的横波速度

图5-1基于页岩气各向异性岩石物理模型对矿物含量的反演结果中的石英含量

图5-2基于页岩气各向异性岩石物理模型对矿物含量的反演结果中的粘土含量

图5-3基于页岩气各向异性岩石物理模型对矿物含量的反演结果中的长石含量

图5-4基于页岩气各向异性岩石物理模型对矿物含量的反演结果中的干酪根含量

具体实施方式

下面结合附图对本发明作进一步详细描述：

发明的内容为针对有机质页岩的各向异性岩石物理建模方法，即包括多种矿物混合的岩石基质模量计算、含孔隙干岩石骨架模量计算、干酪根有机质的植入和含流体岩石模量计算一整套建模流程。页岩岩石物理之所以区别于常规油气藏，一方面是因为岩石中发育的复杂孔隙，另一方面是由于应力和有机质带来的各向异性，在该建模方法中，分别考虑了不同孔隙和有机质对岩石模量带来的影响。为了更好的说明整个建模方法及流程思路，下面结合图1来展示说明。

第一步，由Reuss-Voigt-Hill平均((1)式)计算多种矿物混合的岩石基质弹性模量：

$M_{\mathrm{VRH}}=\frac{M_v+M_R}{2},---\left(1\right)$

其中

$M_V=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{f}_i{M}_i,\frac{1}{M_R}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}\frac{f_i}{M_i}$

f_i和M_i分别代表第i种矿物的比例和模量；N表示N种矿物；

对于陆相页岩，这些矿物通常包括：粘土、石英、长石等；而对于海相页岩，通常包括：粘土、白云石、石灰石等。各种矿物的参数可以通过X射线衍射分析得到。

第二步，利用各向同性等效微分介质理论(DEM)向岩石基质中依次加入不同形状的孔隙(加入的孔隙依实际地区的情况而定，并不是一成不变的。例如对于砂泥岩储层可以加入10％的孔隙，而对于碳酸盐岩储层，只能加入2％左右的孔隙；同样，孔隙的形状也是依实际情况而定的，该模型可用于不同地区的岩石物理建模)，不同的形状由孔隙宽长比来刻画(如刚性孔隙宽长比一般设为0.2-0.5，软孔隙或微裂隙宽长比一般小于0.1)，计算得到含孔隙干岩石的模量(等效体积模量K^*和剪切模量μ^*)。每加入一种孔隙都需要计算一次，但计算的结果是以前一次为基础的。例如加入第一种孔隙后，可以得到模量1，然后在此基础上加入第二种孔隙得到模量2，依次加入不同形状的孔隙。

向基质中加入矿物或孔隙计算等效体积模量K^*和剪切模量μ^*的DEM表达式为(Berryman，1992b)：

$(1-y)\frac{d}{\mathrm{dy}}\left[K^{*}\left(y\right)\right]=(K_2-K^{*})P^{*2}\left(y\right),$

$(1-y)\frac{d}{\mathrm{dy}}\left[{\mathrm{\μ}}^{*}\left(y\right)\right]=({\mathrm{\μ}}_2-{\mathrm{\μ}}^{*})Q^{*2}\left(y\right)---\left(2\right)$

初始条件为：K^*(0)=K₁，μ^*(0)=μ₁，K₁和μ₁为基质(相1)的体积和剪切模量；K₂和μ₂为矿物(相2)的体积模量和剪切模量，y为所加入相2的体积比例。对于流体和孔隙空腔，y等于孔隙度。参数P和Q用以刻画任意形状的包裹物(相2)对岩石的影响，具体如下：

$P=\frac{1}{3}{T}_{\mathrm{iijj}}$

$Q=\frac{1}{5}(T_{\mathrm{ijij}}-\frac{1}{3}{T}_{\mathrm{iijj}})---\left(A1\right)$

Berryman(1980b)给出了计算P和Q所需的相关标量：

T_iijj=3F₁／F₂；

$T_{\mathrm{ijij}}-\frac{1}{3}{T}_{\mathrm{iijj}}=\frac{2}{F_3}+\frac{1}{F_4}+\frac{F_4{F}_5+F_6{F}_7-F_8{F}_9}{F_2{F}_4}---\left(A2\right)$

其中：

$F_1=1+A[\frac{3}{2}(f+\mathrm{\θ})-R(\frac{3}{2}f+\frac{5}{2}\mathrm{\θ}-\frac{4}{3})]$

$F_2=1+A[1+\frac{3}{2}(f+\mathrm{\θ})-\frac{R}{2}(3f+5\mathrm{\θ})]+B(3-4R)+\frac{A}{2}(A+3B)(3-4R)[f+\mathrm{\θ}-R(f-\mathrm{\θ}+{2\mathrm{\θ}}^2)]$

$F_3=1+A[1-(f+\frac{3}{2}\mathrm{\θ})+R(f+\mathrm{\θ})]$

$F_4=1+\frac{A}{4}[f+3\mathrm{\θ}-R(f-\mathrm{\θ})]$

$F_5=A[-f+R(f+\mathrm{\θ}-\frac{4}{3})]+\mathrm{B\θ}(3-4R)$

F₆=1+A[1+f-R(f+θ)]+B(1-θ)(3-4R)

$F_7=2+\frac{A}{4}[3f+9\mathrm{\θ}-R(3f+5\mathrm{\θ})]+\mathrm{B\θ}(3-4R)$

$F_8=A[1-2R+\frac{f}{2}(R-1)+\frac{\mathrm{\θ}}{2}(5R-3)]+B(1-\mathrm{\θ})(3-4R)$

F₉=A[(R-1)f-Rθ]+Bθ(3-4R)   (A3)

其中A，B和R为：

$A=\frac{{\mathrm{\μ}}_i}{{\mathrm{\μ}}_m}-1$

$B=\frac{1}{3}(\frac{K_i}{K_m}-\frac{{\mathrm{\μ}}_i}{{\mathrm{\μ}}_m})$

$R=\frac{1-{2v}_m}{2(1-v_m)}---\left(A4\right)$

函数θ为：

$\mathrm{\&theta;}=\left\{\begin{array}{cc}\frac{\mathrm{\&alpha;}}{{({\mathrm{\&alpha;}}^2-1)}^{\frac{3}{2}}}[\mathrm{\&alpha;}{({\mathrm{\&alpha;}}^2-1)}^{\frac{1}{2}}-\mathrm{arccosh}\mathrm{\&alpha;}]& (\mathrm{\&alpha;}>1)\\ \frac{\mathrm{\&alpha;}}{{(1-{\mathrm{\&alpha;}}^2)}^{\frac{3}{2}}}[\arccos \mathrm{\&alpha;}-\mathrm{\&alpha;}{(1-{\mathrm{\&alpha;}}^2)}^{\frac{1}{2}}]& (\mathrm{\&alpha;}<1)\end{array}\right.---\left(A5\right)$

对扁圆球体α<1，扁长球体α>1；

第三步，利用各向异性等效微分介质理论(Hornby，1994)向干岩石中加入干酪根有机质，表达式为：

$\frac{d}{\mathrm{dy}}\left(c_{\mathrm{ijkl}}^{\mathrm{DEM}}\left(y\right)\right)=\frac{1}{1-y}(c_{\mathrm{ijkl}}^{*}-c_{\mathrm{ijkl}}^{\mathrm{DEM}}\left(y\right))*{[I+G(c_{\mathrm{ijkl}}^{*}-c_{\mathrm{ijkl}}^{\mathrm{DEM}}\left(y\right))]}^{-1}---\left(3\right)$

式中c为各向异性介质中的刚度系数，c^DEM和c^*分别表示基质刚度矩阵和包裹物(这里指有机质)的刚度矩阵(即含干酪根的干岩石刚度矩阵，是通过求解公式(3)得到该矩阵的)。I为单位矩阵，G为刻画不同形状有机质对整个岩石模量影响的参数，具体如下：

对于TI介质，用以下参数来表示刚度系数：

d=c₁₁，e=(c₁₁-c₁₂)／2，f=c₄₄，g=c₁₃+c₄₄，h=c₃₃   (B1)张量包含了12个非零元素，分别表达为：

$\begin{array}{c}{\stackrel{\&OverBar;}{G}}_{1111}={\stackrel{\&OverBar;}{G}}_{2222}\\ =(\mathrm{\&pi;}/2){\&Integral;}_0^1\mathrm{\&Delta;}(1-x^2)\left\{\right[(1-x^2)+h{\mathrm{\&rho;}}^2{x}^2\left]\right[(3e+d)(1-x^2)+4f{\mathrm{\&rho;}}^2{x}^2]\\ -g^2{\mathrm{\&rho;}}^2{x}^2(1-x^2)\}\mathrm{dx}\end{array}---\left(B2\right)$

${\stackrel{\&OverBar;}{G}}_{3333}=4\mathrm{\&pi;}{\&Integral;}_0^1\mathrm{\&Delta;}{\mathrm{\&rho;}}^2{x}^2[d(1-x^2)+f{\mathrm{\&rho;}}^2{x}^2][e(1-x^2)+f{\mathrm{\&rho;}}^2{x}^2]\mathrm{dx}---\left(B3\right)$

$\begin{array}{c}{\stackrel{\&OverBar;}{G}}_{1122}={\stackrel{\&OverBar;}{G}}_{2211}\\ =(\mathrm{\&pi;}/2){\&Integral;}_0^1\mathrm{\&Delta;}(1-x^2)\left\{\right[f(1-x^2)+h{\mathrm{\&rho;}}^2{x}^2\left]\right[(e+3d)(1-x^2)+4f{\mathrm{\&rho;}}^2{x}^2]\\ -{3g}^2{\mathrm{\&rho;}}^2{x}^2(1-x^2)\}\mathrm{dx}\end{array}---\left(B4\right)$

$\begin{array}{c}{\stackrel{\&OverBar;}{G}}_{1133}={\stackrel{\&OverBar;}{G}}_{2233}\\ =2\mathrm{\&pi;}{\&Integral;}_0^1\mathrm{\&Delta;}{\mathrm{\&rho;}}^2{x}^2\left\{\right[(d+e)(1-x^2)+2f{\mathrm{\&rho;}}^2{x}^2\left]\right[f(1-x^2)+{\mathrm{h\&rho;}}^2{x}^2]\\ -g^2{\mathrm{\&rho;}}^2{x}^2(1-x^2)\}\mathrm{dx}\end{array}---\left(B5\right)$

$\begin{array}{c}{\stackrel{\&OverBar;}{G}}_{3311}={\stackrel{\&OverBar;}{G}}_{3322}\\ =2\mathrm{\&pi;}{\&Integral;}_0^1\mathrm{\&Delta;}(1-x^2)[d(1-x^2)+f{\mathrm{\&rho;}}^2{x}^2]\left[e(1-x^2+f{\mathrm{\&rho;}}^2{x}^2)\right]\mathrm{dx}\end{array},-,-,-,\left(B6\right)$

${\stackrel{\&OverBar;}{G}}_{1212}=(\mathrm{\&pi;}/2){\&Integral;}_0^1\mathrm{\&Delta;}{(1-x^2)}^2\{g^2{\mathrm{\&rho;}}^2{x}^2-(d-e)[f(1-x^2)+h{\mathrm{\&rho;}}^2{x}^2\left]\right\}\mathrm{dx}---\left(B7\right)$

$\begin{array}{c}{\stackrel{\&OverBar;}{G}}_{1313}={\stackrel{\&OverBar;}{G}}_{2323}\\ =-2\mathrm{\&pi;}{\&Integral;}_0^1\mathrm{\&Delta;g}{\mathrm{\&rho;}}^2{x}^2(1-x^2)[e(1-x^2)+f{\mathrm{\&rho;}}^2{x}^2]\mathrm{dx}\end{array}---\left(B8\right)$

上式中的Δ表达为：

Δ^-1=[e(1-x²)+fρ²x²]{[d(1-x²)+fρ²x²][f(1-x²)+hρ²x²]-g²ρ²x²(1-x²)}   (B9)ρ为宽长比的倒数，即ρ=1／α。

这里可以将干岩石作为背景介质，也可将干酪根作为背景介质，具体依赖于不同的地质条件和岩芯资料情况。干酪根的形状同样由宽长比决定，由于干酪根有机质在页岩岩石中通常会呈一定形状排练，因此会带来较强的VTI各向异性，计算得到的富含有机质干岩石模量以刚度系数来表达。

第四步，利用各向异性Gassmann方程向含有机质干岩石中加入流体(可以是水、气、油或任意比例的混合)。利用的表达式为(Gassmann，1951)：

$c_{\mathrm{ijkl}}^{\mathrm{sat}}=c_{\mathrm{ijkl}}^{\mathrm{dry}}+\frac{(K_0{\mathrm{\&delta;}}_{\mathrm{ij}}-c_{\mathrm{ijaa}}^{\mathrm{dry}}/3)(K_0{\mathrm{\&delta;}}_{\mathrm{kl}}-c_{\mathrm{bbkl}}^{\mathrm{dry}}/3)}{(K_0/K_{\mathrm{fl}})\mathrm{\&phi;}(K_0-K_{\mathrm{fl}})+(K_0-c_{\mathrm{ccdd}}^{\mathrm{dry}}/9)}---\left(4\right)$

(公式(4)中哪里用到了第三步)

为干岩石刚度系数(就是指第二步得到的含孔隙干岩石的骨架模量)，为含流体岩石的刚度系数(就是指第四步最终要求的整个有机页岩的弹性模量)，K₀为矿物体积模量(可以从实验室中测量得到)，K_fl为流体体积模量(可以从实验室中测量得到)，φ为孔隙度(孔隙度的测量通常来自于测井解释或岩芯资料)。

${\mathrm{\&delta;}}_{\mathrm{ij}}=\left\{\begin{array}{cc}1,& i=j,\\ 0& i\&NotEqual;j,\end{array}\right.---\left(5\right)$

由此，可最终计算得到富含有机质页岩的刚度系数。对于VTI介质，相应的速度及各向异性参数表达为(Thomsen，1986)：

垂向速度：

$V_{P0}=\sqrt{C_{33}/\mathrm{\&rho;}},$

$V_{S0}=\sqrt{C_{44}/\mathrm{\&rho;}},---\left(6\right)$

各向异性参数：

$\mathrm{\&epsiv;}=\frac{C_{11}-C_{33}}{{2C}_{33}},$

$\mathrm{\&gamma;}=\frac{C_{66}-C_{44}}{2C_{44}},$

${\mathrm{\&delta;}}^{*}=\frac{1}{{2C}_{33}^2}[2{(C_{13}+C_{44})}^2-(C_{33}-C_{44})C_{11}+C_{33}-{2C}_{44})],---\left(7\right)$

岩石物理模型并不是一种具体的物理模型，它侧重以一种数学的方式来描述物理模型，因此最后的结果即得到岩石的刚度矩阵(即整个有机页岩的弹性模量)，由公式(4)来表达，公式4的结果是以前面三步为基础的。

具体实施时，第一步输入组成岩石矿物的含量，输出的是岩石基质的模量；第二步以岩石基质模量为输入，加入不同的类型孔隙得到干岩石骨架模量；第三步在干岩石的基础上加入干酪根的影响，并考虑干酪根对介质带来的影响，输出含干酪根的干岩石刚度矩阵；第四部，以第三步结果做为输入，在岩石中加入流体，最终得到整个有机页岩的弹性模量。

为验证有机页岩各向异性岩石物理模型的有效性和应用性，首先利用该建模方法分析了岩石弹性参数、各向异性参数及脆性参数随粘土含量及有机质(干酪根)含量的变化情况如图2-1至图2-4所示、如图3-1至图3-2所示，得出的结论与实验室测量结果有较好的一致性。接下来，将该模型应用于油田中的页岩气勘探生产，通过岩芯资料标定验证了该模型的准确性，如图4-1至图4-4所示。图5-1至图5-4为利用页岩岩石物理模型对各种矿物含量进行的反演结果。

非常规油气藏，如页岩气，正日益成为我国重点勘探和开发的目标。由于对该类型的油气藏缺少足够的地球物理认识，因此，如何利用地球物理来对非常规油气藏进行定量刻画目前面临着巨大的挑战。了解有机质、定向分布的应力所带来的各向异性及地球物理响应，对于非常规储层地球物理来说至关重要。该发明提出了一种新的包含有机质的页岩各向异性岩石物理模型，建立了矿物成分、孔隙度、TOC(有机碳总量)及含流体性与岩石弹性参数之间的定量关系。该模型通过了基于岩芯和测井资料的标定，充分证明了其具有较好的应用前景。

上述技术方案只是本发明的一种实施方式，对于本领域内的技术人员而言，在本发明公开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本发明上述具体实施方式所描述的方法，因此前面描述的方式只是优选的，而并不具有限制性的意义。

Claims (5)

1.一种针对富含有机质页岩的各向异性岩石物理建模方法，其特征在于：所述方法包括：

第一步，输入组成岩石矿物的比例和模量，计算得到多种矿物混合的岩石基质的弹性模量；

第二步，向所述岩石基质中依次加入不同类型的孔隙得到干岩石，计算得到含孔隙干岩石的骨架模量；

第三步，向所述干岩石中加入干酪根得到含有机质干岩石，计算得到含干酪根的干岩石刚度矩阵；

第四步，向所述含有机质干岩石中加入流体，计算得到整个有机页岩的弹性模量。

2.根据权利要求1所述的针对富含有机质页岩的各向异性岩石物理建模方法，其特征在于：所述第一步是这样实现的：

由公式(1)计算多种矿物混合的岩石基质弹性模量：

$M_{\mathrm{VRH}}=\frac{M_v+M_R}{2},---\left(1\right)$

其中

$M_V=\underset{i=1}{\overset{N}{\mathrm{\&Sigma;}}}{f}_i{M}_i,\frac{1}{M_R}=\underset{i=1}{\overset{N}{\mathrm{\&Sigma;}}}\frac{f_i}{M_i}$

f_i和M_i分别代表第i种矿物的比例和模量；N表示N种矿物。

3.根据权利要求2所述的针对富含有机质页岩的各向异性岩石物理建模方法，其特征在于：所述第二步是这样实现的：

基于第一步得到的多种矿物混合的岩石基质的弹性模量，利用各向同性等效微分介质理论向岩石基质中依次加入不同形状的孔隙，计算得到含孔隙干岩石的模量：等效体积模量K^*和剪切模量μ^*；

不同的形状是通过孔隙的宽长比来刻画的；

计算得到含孔隙干岩石的模量是这样实现的：

每加入一种孔隙，计算一次含孔隙干岩石的模量，具体如下：加入第一种孔隙，计算得到第一次含孔隙干岩石的模量，然后在此基础上加入第二种孔隙，计算得到第二次含孔隙干岩石的模量，以此类推，直到加入最后一种孔隙，然后计算得到最后的含孔隙干岩石的模量。

4.根据权利要求3所述的针对富含有机质页岩的各向异性岩石物理建模方法，其特征在于：所述第三步是这样实现的：

基于第二步得到的含孔隙干岩石的模量，利用各向异性等效微分介质理论向干岩石中加入干酪根，计算得到含干酪根的干岩石刚度矩阵。

5.根据权利要求4所述的针对富含有机质页岩的各向异性岩石物理建模方法，其特征在于：所述第四步是这样实现的：

利用各向异性Gassmann方程向含有机质干岩石中加入流体，计算得到整个有机页岩的弹性模量。