裂缝储层岩石物理建模方法及系统与流程

1.本发明涉及地球物理勘探领域，更具体地，涉及一种裂缝储层岩石物理建模方法及系统。


背景技术：

2.岩石物理建模技术包含等效介质理论、孔隙弹性理论，等效介质理论分为不考虑几何形态影响的等效介质理论(通常所指的界限模型)和考虑几何形态的岩石物理模型。不考虑几何细节影响的界限模型的发展较为提前，包括voight和reuss提出的v-r界限，hashin与shtrikman给出的h-s界限，hill求取平均数给出v-r-h模量。后续又发展出了考虑几何细节影响的各向同性等效介质模型sca，微分等效介质模型dem模型和backus平均理论。裂缝等效介质理论硏究的发展是建立在材料科学基础上的，hudson提出单套扁圆币状裂缝理论用于研究ti介质schoenberg提出了广泛适用的线性滑动理论。现有技术只能通过岩石物理模型构建地震参数(振幅、相位信息等)和弹性参数(泊松比，杨氏模量等)之间的联系，不能直接推广到裂缝信息(裂缝密度，裂缝内流体类型等)，并且假设为各向同性介质建模。因此，有必要开发一种裂缝储层岩石物理建模方法及系统。
3.公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。


技术实现要素：

4.本发明提出了一种裂缝储层岩石物理建模方法及系统，其能够通过建立双层介质模型研究页岩储层和碳酸盐岩储层特性，通过改变裂缝储层相关油藏参数，分析不同情况下的avo响应规律，对裂缝型油气藏进行定性描述得出变化规律，反射系数方位各向异性特点对裂缝方位预测有重要意义。
5.根据本发明的一方面，提出了一种裂缝储层岩石物理建模方法。所述方法可以包括：通过各向同性基质建模，建立背景介质岩石体；根据所述hudson模型与所述schoenberg模型，在所述背景介质岩石体嵌入一组垂直发育的裂缝，构建ti介质模型；针对所述ti介质模型的等效岩性进行brown-korringa各向异性流体替换，混入不同油气水比例的孔隙，获得最终的裂缝储层岩石物理模型。
6.优选地，所述hudson模型的弹性系数矩阵为：
7.c＝c0+c1+c2ꢀꢀꢀ
(1)
8.其中，c0为个各向同性背景介质的弹性刚度系数矩阵，c1，c2分别为一、二阶修正分量，一阶修正量为：
[0009][0010]
二阶修正量为：
[0011][0012]
其中，
[0013]
优选地，所述schoenberg模型的刚度矩阵为：
[0014][0015]
其中，δn为法向柔度，δt为切向柔度，
[0016]
优选地，brown-korringa模型为：
[0017]
[0018]
其中，s
(dry)
、s
(sat)
、s0分别为干岩石、流体饱和岩石、岩石组成矿物成分的柔度表示形式，b
fl
、b0分别为岩石孔隙内含有的流体、岩石自身矿物质的可挤压性，f为岩石孔隙度。
[0019]
优选地，还包括：针对所述最终的裂缝储层岩石物理模型进行地震相应分析，确定裂缝参数随着地震响应参数与裂缝储层参数的变化规律与响应特征。
[0020]
根据本发明的另一方面，提出了一种裂缝储层岩石物理建模系统，其特征在于，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：通过各向同性基质建模，建立背景介质岩石体；根据所述hudson模型与所述schoenberg模型，在所述背景介质岩石体嵌入一组垂直发育的裂缝，构建ti介质模型；针对所述ti介质模型的等效岩性进行brown-korringa各向异性流体替换，混入不同油气水比例的孔隙，获得最终的裂缝储层岩石物理模型。
[0021]
优选地，所述hudson模型的弹性系数矩阵为：
[0022]
c＝c0+c1+c2ꢀꢀꢀ
(1)
[0023]
其中，c0为个各向同性背景介质的弹性刚度系数矩阵，c1，c2分别为一、二阶修正分量，一阶修正量为：
[0024][0025]
二阶修正量为：
[0026][0027]
其中，
[0028]
优选地，所述schoenberg模型的刚度矩阵为：
[0029][0030]
其中，δn为法向柔度，δt为切向柔度，
[0031]
优选地，brown-korringa模型为：
[0032][0033]
其中，s
(dry)
、s
(sat)
、s0分别为干岩石、流体饱和岩石、岩石组成矿物成分的柔度表示形式，b
fl
、b0分别为岩石孔隙内含有的流体、岩石自身矿物质的可挤压性，f为岩石孔隙度。
[0034]
优选地，还包括：针对所述最终的裂缝储层岩石物理模型进行地震相应分析，确定裂缝参数随着地震响应参数与裂缝储层参数的变化规律与响应特征。
[0035]
其有益效果在于：相比现有技术只能通过岩石物理模型构建地震参数(振幅、相位信息等)和弹性参数(泊松比，杨氏模量等)之间的联系的不足，本发明能直接推广到裂缝信息(裂缝密度，裂缝内流体类型等)，裂缝储层建模更能够直接的建立裂缝储层参数信息和地震参数之间的直接联系，为后续直接反演裂缝参数信息提供理论依据。
[0036]
本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
[0037]
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
[0038]
图1示出了根据本发明的裂缝储层岩石物理建模方法的步骤的流程图。
[0039]
图2示出了根据本发明的一个实施例的建立背景介质岩石体的流程示意图。
[0040]
图3a和图3b分别示出了根据本发明的一个实施例的schoenberg线性滑动理论模型的裂缝参数与法向、切向分量之间的变化关系曲线的示意图。
[0041]
图4a和图4b分别示出了根据本发明的一个实施例的油、气、水对裂缝纵横比与裂缝密度的敏感度的示意图。
[0042]
图5a、图5b和图5c分别示出了根据本发明的一个实施例的水油混合、水气混合、油气混合的岩石物理参数随裂缝纵横比的变化示意图。
[0043]
图6a、图6b和图6c分别示出了根据本发明的一个实施例的ti介质相速度随各向异
性参数ε、裂缝密度与裂缝纵横比的变化示意图。
[0044]
图7a、图7b、图7c、图7d、图7e和图7f分别示出了根据本发明的一个实施例的各向异性参数对反射系数的影响示意图。
[0045]
图8示出了根据本发明的一个实施例的ti储层反射系数的三维立体图。
具体实施方式
[0046]
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
[0047]
图1示出了根据本发明的裂缝储层岩石物理建模方法的步骤的流程图。
[0048]
在该实施例中，根据本发明的裂缝储层岩石物理建模方法可以包括：步骤101，通过各向同性基质建模，建立背景介质岩石体；步骤102，根据hudson模型与schoenberg模型，在背景介质岩石体嵌入一组垂直发育的裂缝，构建ti介质模型；步骤103，针对ti介质模型的等效岩性进行brown-korringa各向异性流体替换，混入不同油气水比例的孔隙，获得最终的裂缝储层岩石物理模型。
[0049]
在一个示例中，hudson模型的弹性系数矩阵为：
[0050]
c＝c0+c1+c2ꢀꢀꢀ
(1)
[0051]
其中，c0为个各向同性背景介质的弹性刚度系数矩阵，c1，c2分别为一、二阶修正分量，一阶修正量为：
[0052][0053]
二阶修正量为：
[0054][0055]
其中，在一个示例中，schoenberg模型的刚度矩阵为：
[0056][0057]
其中，δn为法向柔度，δt为切向柔度，
[0058]
在一个示例中，brown-korringa模型为：
[0059][0060]
其中，s
(dry)
、s
(sat)
、s0分别为干岩石、流体饱和岩石、岩石组成矿物成分的柔度表示形式，b
fl
、b0分别为岩石孔隙内含有的流体、岩石自身矿物质的可挤压性，f为岩石孔隙度。
[0061]
在一个示例中，还包括：针对最终的裂缝储层岩石物理模型进行地震相应分析，确定裂缝参数随着地震响应参数与裂缝储层参数的变化规律与响应特征。
[0062]
具体地，根据本发明的裂缝储层岩石物理建模方法可以包括：
[0063]
通过各向同性基质建模，建立背景介质岩石体：选定背景基质的岩性，假定背景介质是各向同性介质；运用sca模型迭代计算，针对页岩储层的矿物成分比例特点，设定加入不同的矿物成分黏土、干酪根等；结合微分等效介质模型dem模型混入不同孔隙度大小的孔隙；通过设定不同水油比例的流体，采用gassmann方程进行流体替换；计算得到背景基质等效岩石体的弹性模量信息。
[0064]
sca模型能够考虑到孔隙形状这一影响因素，能够适用于孔隙度较大的岩石弹性模量的求取。该模型能够充分考虑到内含物之间的相互作用，适用于求取大孔隙度的岩石
弹性模量。berryman给出了更为通用的n相介质的自适应模型迭代公式如下：
[0065][0066]
其中，p、q为形状因子，k
sc
与μ
sc
方程是耦合的，需要通过多次迭代的方法来求解等效体的弹性模量值。
[0067]
微分有效介质模型dem中，各相内含物加入的顺序会导致有效模量的改变。berryman建立了关于体积模量k与剪切模量μ的相互耦合的微分方程如下：
[0068][0069]
其中要满足k
*
(0)＝k1，μ
*
(0)＝μ1等式中符号k1，μ1分别为1(原始基质)的体积模量和剪切模量；其中k2，μ2分别为2(外来内含物)的体积模量和剪切模量，y的含义为2的体积含量。p和q为形状因子分量，上标*2代表的是内含物2在基质k*和μ*中的影响因子。
[0070]
gassmann方程应用于干岩石到流体饱和岩石之间流体替换环节，可以求得流体饱和岩石的等效弹性参数，并且建立了于纵横波速度之间的联系，方便直接应用和换算：
[0071][0072]
由以上公式，可以计算流体饱和情况下的岩石体等效弹性模量，只需要将上述公式经过转化表达改写，再结合岩石密度与孔隙度关系式便可以得到纵横比速度求取表达式如下：
[0073][0074]
根据hudson模型与schoenberg模型，分析背景介质岩石体的裂缝参数与弹性参数之间的相关关系、裂缝参数的变化特点以及裂缝参数对含油气的敏感程度；根据裂缝参数的变化特点，根据hudson模型与schoenberg模型，在背景介质岩石体嵌入一组垂直发育的裂缝，构建ti介质模型。
[0075]
hudson模型弹性系数矩阵为公式(1)，其一阶修正量为公式(2)，二阶修正量为公式(3)。
[0076]
schoenberg模型中假设裂缝面为位移间断，但是界面两侧应力保持连续。线性滑动理论的思想是背景介质和裂缝柔度叠加，表达式为：
[0077]
s＝s0+δs
ꢀꢀꢀ
(10)
[0078]
求逆得到s-1
即为弹性刚度矩阵，整合有c＝c
0-c0δsc0形式。为简化表达，引入无量纲裂缝法向、切向柔度δn，δt表征裂缝，取值范围0-1，则schoenberg刚度矩阵为公式(4)。
[0079]
根据裂缝参数对含油气的敏感程度，针对ti介质模型的等效岩性进行brown-korringa各向异性流体替换，混入不同油气水比例的孔隙，获得最终的裂缝储层岩石物理模型。brown-korringa模型公式(5)被认为是各向异性流体替换公式，可以用来表示各向异性条件下流体饱和岩石的有效弹性参数同各向异性干岩石有效弹性参数两者之间的必然关系。
[0080]
针对最终的裂缝储层岩石物理模型进行地震相应分析，确定裂缝参数随着地震响应参数与裂缝储层参数的变化规律与响应特征。
[0081]
本方法通过建立双层介质模型研究页岩储层和碳酸盐岩储层特性，通过改变裂缝储层相关油藏参数，分析不同情况下的avo响应规律，对裂缝型油气藏进行定性描述得出变化规律，反射系数方位各向异性特点对裂缝方位预测有重要意义。
[0082]
应用示例
[0083]
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
[0084]
图2示出了根据本发明的一个实施例的建立背景介质岩石体的流程示意图。建模流程首先是一个岩石体基质，向其中加入不同岩性组分的颗粒基质成分；然后进行不同水油比例流体替换填充孔隙内结构；再向岩石体中嵌入平行裂缝组；再成分考虑到ti介质特点，进行各向异性流体替换嵌入水油混合的流体。
[0085]
图3a和图3b分别示出了根据本发明的一个实施例的schoenberg线性滑动理论模型的裂缝参数与法向、切向分量之间的变化关系曲线的示意图。设定裂缝内填充物为黏土矿物，裂缝参数随裂缝形状呈非线性曲线变化，随裂缝密度呈线性变化。
[0086]
图4a和图4b分别示出了根据本发明的一个实施例的油、气、水对裂缝纵横比与裂缝密度的敏感度的示意图。可以看出裂缝法向分量δn在不同流体类型填充的情况下显示为不同的数值，说明在油气水不同流体填充情况下，δn对流体类型有较好的分辨能力。而且随裂缝纵横比以及裂缝密度的变化程度不同。根据这些规律，可以通过裂缝的法向和切向分量来估测流体信息。
[0087]
图5a、图5b和图5c分别示出了根据本发明的一个实施例的水油混合、水气混合、油气混合的岩石物理参数随裂缝纵横比的变化示意图，其中，箭头指向裂缝纵横比asp增大方向0.001、0.003、0.005、0.007、0.01。可以看出裂缝法向分量δn随含水饱和度或含油饱和度的增大而呈现非线性减小变化，气水混合裂缝的减小其实更明显，油水裂缝次之，最后是油气裂缝。说明裂缝内流体不同对裂缝纵横比影响程度较明显，可以运用这一特点很好的区分并指示流体类型。
[0088]
图6a、图6b和图6c分别示出了根据本发明的一个实施例的ti介质相速度随各向异性参数ε、裂缝密度与裂缝纵横比的变化示意图。图6a可以看出ti各向异性介质中的纵波相速度形状为一个椭圆而非圆形，说明具有椭圆各向异性特点，随着ε的增大速度椭圆形状更加明显，即各向异性特点增强。图6b是对裂缝密度影响的研究，裂缝纵横比选定为0.01，裂
缝密度取值依次增大变化，具体数值为0.01、0.03、0.05、0.1、0.2，图中极坐标中的箭头方向代表裂缝密度取值增大的方向，由图可知随着裂缝密度的增大，相速度椭圆各向异性特点逐渐增强。图6c为对裂缝纵横比的影响，裂缝密度为0.01，裂缝纵横比取值从0.001变化到0.02，间隔为0.005变化的四个数值，由图可知随着裂缝纵横比的增大，相速度椭圆各向异性特点逐渐增强，并且从数值大小上分析来看，裂缝密度比裂缝纵横比对相速度各向异性的影响程度更大。
[0089]
图7a、图7b、图7c、图7d、图7e和图7f分别示出了根据本发明的一个实施例的各向异性参数对反射系数的影响示意图，其中，箭头指向参数取值增大方向，由上到下依次为变化ε、δ、γ取值情况，研究纵波反射系数分别随入射角和方位角的变化影响情况，可以得出反射系数随各向异性参数变化敏感度规律为：对γ变化最敏感，其次为δ，最弱为ε。根据这个规律，应用在avaz反演中，可以参考各向异性参数的敏感度情况，在反射系数近似式简化改写提供理论依据。
[0090]
图8示出了根据本发明的一个实施例的ti储层反射系数的三维立体图。反射系数最大值处，能量速度最强，此时波传播方向与裂缝发育方向平行一致，根据反射系数随方位变化情况，预测ti介质裂缝储层裂缝发育方向。
[0091]
综上所述，本发明通过建立双层介质模型研究页岩储层和碳酸盐岩储层特性，通过改变裂缝储层相关油藏参数，分析不同情况下的avo响应规律，对裂缝型油气藏进行定性描述得出变化规律，反射系数方位各向异性特点对裂缝方位预测有重要意义。
[0092]
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
[0093]
根据本发明的实施例，提供了一种裂缝储层岩石物理建模系统，其特征在于，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：通过各向同性基质建模，建立背景介质岩石体；根据hudson模型与schoenberg模型，在背景介质岩石体嵌入一组垂直发育的裂缝，构建ti介质模型；针对ti介质模型的等效岩性进行brown-korringa各向异性流体替换，混入不同油气水比例的孔隙，获得最终的裂缝储层岩石物理模型。
[0094]
在一个示例中，hudson模型的弹性系数矩阵为：
[0095]
c＝c0+c1+c2ꢀꢀꢀ
(1)
[0096]
其中，c0为个各向同性背景介质的弹性刚度系数矩阵，c1，c2分别为一、二阶修正分量，一阶修正量为：
[0097][0098]
二阶修正量为：
[0099][0100]
其中，
[0101]
在一个示例中，schoenberg模型的刚度矩阵为：
[0102][0103]
其中，δn为法向柔度，δt为切向柔度，
[0104]
在一个示例中，brown-korringa模型为：
[0105]
[0106]
其中，s
(dry)
、s
(sat)
、s0分别为干岩石、流体饱和岩石、岩石组成矿物成分的柔度表示形式，b
fl
、b0分别为岩石孔隙内含有的流体、岩石自身矿物质的可挤压性，f为岩石孔隙度。
[0107]
在一个示例中，还包括：针对最终的裂缝储层岩石物理模型进行地震相应分析，确定裂缝参数随着地震响应参数与裂缝储层参数的变化规律与响应特征。
[0108]
本系统通过建立双层介质模型研究页岩储层和碳酸盐岩储层特性，通过改变裂缝储层相关油藏参数，分析不同情况下的avo响应规律，对裂缝型油气藏进行定性描述得出变化规律，反射系数方位各向异性特点对裂缝方位预测有重要意义。
[0109]
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
[0110]
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。