基于岩石孔隙结构参数的地层孔隙流体压力的计算方法与流程

1.本发明涉及岩石物理学技术领域，具体为一种基于岩石孔隙结构参数的地层孔隙流体压力的计算方法。


背景技术：

2.地层孔隙流体压力定义为任一多孔隙体系，当其承受外加负荷时，由孔隙中流体所支撑的那部分负载称为孔隙流体压力。在正常情况下，当孔隙与外界连通，孔隙中的流体可以自由排出时，孔隙流体压力等于静水压力；反之，当孔隙与外界相互阻断时，构成异常地层压力，异常地层压力分布是含油气沉积盆地普遍存在的现象,有效预测地层压力对于了解油气富集、运移状况和保证钻井安全具有重要意义。
3.现有计算地层孔隙流体压力的方法可概括分为三类：第一类为建立泥岩或页岩的正常压实趋势线判断地层是否超压，再根据等效深度法和有效应力原理估算地层孔隙流体压力；第二类利用泥岩或页岩的正常压实趋势线方程建立泥岩或页岩的有效应力与速度之间的关系，再根据有效应力原理估算地层孔隙流体压力；第三类为直接建立速度与有效应力之间的关系，不依赖于泥岩或页岩的正常压实趋势线方程，根据有效应力原理估算地层孔隙流体压力。上述三类估算地层孔隙流体压力的方法本质均为建立有效应力与速度的关系，再根据有效应力原理估算地层孔隙流体压力，仅适用于单一岩性(同种岩性)泥岩、页岩或砂岩，对于岩性复杂的地层，尤其是混合岩性的，利用上述三类方法进行地层孔隙流体压力估算时存在诸多的不确定性，导致估算结果的可靠度低，不利于了解油气富集和运移状况，钻井存在安全隐患。


技术实现要素：

4.本发明意在提供一种基于岩石孔隙结构参数的地层孔隙流体压力的计算方法，以解决现有技术计算地层孔隙流体压力的方法不适用复杂岩性地层的问题。
5.为了实现上述目的，本发明提供如下技术方案：
6.基于岩石孔隙结构参数的地层孔隙流体压力的计算方法，包括以下步骤：
7.s101：利用常规测井资料，先计算出各测井曲线点下的饱和岩石体积模量，接着计算各孔隙不同流体的流体体积模量，再计算孔隙混合流体的流体体积模量；
8.s102：利用常规测井资料，基于s101计算饱和岩石体积模量和流体体积饱和模量，并根据自适应基质体积模量反演理论，计算各个测井点处的岩石基质体积模量；
9.s103：基于s101得到的所述孔隙混合流体的流体体积模量、s102得到的所述饱和岩石体积模量和所述岩石基质体积模量，确定表示孔隙结构的综合孔隙结构参数；
10.s104：计算考虑孔隙结构参数影响的岩石biot系数，基于实际地质资料及其测井资料，分析其孔隙结构特征，得到可调节参数c的取值范围；
11.s105：利用所计算的考虑孔隙结构参数影响的biot，基于有效应力理论及aplin表达式计算得出地层孔隙流体压力。
12.进一步地，所述常规测井资料包括密度测井曲线、孔隙度测井曲线、含水饱和度测井曲线、纵波速度测井曲线和横波速度测井曲线。
13.进一步地，利用测井资料所得的饱和岩石纵波速度、饱和岩石横波速度和饱和岩石密度计算s101中所述饱和岩石体积模量的表达式为：
[0014][0015]
上述表达式中，k
sat
为饱和岩石体积模量，v
p
为饱和岩石的纵波速度，vs为饱和岩石的横波速度，ρ
sat
为饱和岩石密度。
[0016]
进一步地，计算s101中各孔隙中不同流体的体积模量的表达式为：
[0017][0018]
上述表达式中，kw、ko、kg分别代表水、油、气的体积模量，h表示岩石的埋藏深度，单位为km。
[0019]
进一步地，计算s101中孔隙混合流体的体积模量表达式为：
[0020][0021]
上述表达式中，k
fl
为孔隙混合流体的体积模量，si分别表示含水、含气、含油饱和度，ki对应表示为水、油、气的体积模量。
[0022]
进一步地，确定s103中综合孔隙结构参数的表达式为：
[0023][0024]
上述表达式中，k
sat
为饱和岩石体积模量，k
fl
为孔隙混合流体体积模量，k0为岩石基质体积模量，为孔隙度，c为可调节参数，c的取值范围在0～1之间，根据岩性不同改变c的取值。
[0025]
进一步地，联系biot系数和孔隙结构参数，计算s104中考虑孔隙结构参数biot系数的表达式为：
[0026][0027]
上述表达式中，β为biot系数，为孔隙度，c为可调参数，c的取值范围在0～1之间，根据岩性不同改变c的取值。
[0028]
进一步地，s105中有效应力理论的表达式为：
[0029]
σ＝p
d-β
·
p
p
[0030]
上述表达式中，σ为有效应力，pd为上覆地层压力，β为考虑岩石孔隙度结构参数的biot系数，p
p
为地层孔隙流体压力。
[0031]
进一步地，计算s105中aplin的表达式为：
[0032]
v＝v
0-aln(σ/σ0)
[0033]
上述表达式中，为空隙率，a为压实系数，σ为有效应力。
[0034]
进一步地，计算s105中地层孔隙流体压力的表达式为：
[0035][0036]
上述表达式中，p
p
为地层孔隙流体压力，pd为上覆岩层压力，σ0初始有效压力，v为空隙率，a为地层压实系数，β为考虑岩石孔隙度结构参数的biot系数。
[0037]
由以上本技术提供的技术方案可见，本技术先利用常规测井资料，计算饱和岩石的体积模量k
sat
和孔隙中混合流体的体积模量k
fl
；然后根据自适应反演计算岩石基质的体积模量k0；基于饱和岩石体积模量k
sat
和岩石基质体积模量k0计算出表征岩石孔隙结构的综合孔隙结构参数s；联系综合孔隙结构参数s与biot系数的关系，得出考虑孔隙结构参数影响的biot系数β，确定可调节参数c的取值范围；最后基于有效应力理论及aplin表达式计算出考虑岩石孔隙结构的地层孔隙流体压力p
p
。本技术中计算地层孔隙流体压力不用建立正常压实趋势线，摆脱由于趋势线预测不准确导致地层孔隙流体压力计算不准确的问题，且本技术中引入孔隙结构参数，综合考虑了复杂岩石的矿物组成成分、孔隙流体、胶结物及孔隙结构参数(包括扁率、固结系数、临界孔隙度等)对于地层孔隙流体压力的影响，使地层孔隙流体压力的计算更加精确，能够适用于复杂岩性地层的地层孔隙流体压力的计算。
附图说明
[0038]
图1为本发明基于岩石孔隙结构参数的地层孔隙流体压力的计算方法的流程图；
[0039]
图2为本发明基于岩石孔隙结构参数的地层孔隙流体压力的计算方法中s101和s102所使用的常规测井曲线图；
[0040]
图3为本发明基于岩石孔隙结构参数的地层孔隙流体压力的计算方法中s102饱和岩石及孔隙混合流体体积模量曲线图；
[0041]
图4为本发明基于岩石孔隙结构参数的地层孔隙流体压力的计算方法中s102自适应反演基质体积模量流程图；
[0042]
图5为本发明基于岩石孔隙结构参数的地层孔隙流体压力的计算方法中s103基质体积模量曲线图；
[0043]
图6为本发明基于岩石孔隙结构参数的地层孔隙流体压力的计算方法中s104综合孔隙结构参数及biot系数曲线图；
[0044]
图7为本发明基于岩石孔隙结构参数的地层孔隙流体压力的计算方法中s105地层孔隙流体压力及压力系数图。
具体实施方式
[0045]
下面结合附图和实施方式对本发明作进一步的详细说明：
[0046]
参考图1至图7，基于岩石孔隙结构参数的地层孔隙流体压力的计算方法，包括以下步骤：
[0047]
s101：利用常规测井资料，先计算出各测井曲线点下的饱和岩石体积模量，然后计算各孔隙流体的流体体积模量，再计算孔隙混合流体的流体体积模量。
[0048]
常规测井资料如图2所示，包括纵波速度测井曲线、横波速度测井曲线、密度测井曲线、孔隙度测井曲线、含水饱和度测井曲线。纵波速度表示为v
p
，单位为米/秒,横波速度表示为vs,，单位为米/秒，密度表示为ρ
sat
，单位为克/立方厘米，孔隙度表示为含水饱和度表示为sw。
[0049]
其中，计算各测井点下的饱和岩石体积模量的表达式如下：
[0050][0051]
式中，ρ
sat
为饱和岩石的密度，v
p
为由测井资料所得到的纵波速度，vs为由测井资料所得到的横波速度，k
sat
为计算得到的饱和岩石体积模量；
[0052]
计算各测井点下的孔隙中不同流体的体积模量的表达式如下：
[0053][0054]
式中，kw、ko、kg分别代表水、油、气的体积模量，h表示岩石的埋藏深度，单位为km；
[0055]
将不同的流体混合之后，计算孔隙中混合流体的体积模量的表达式如下：
[0056][0057]
式中,k
fl
为孔隙混合流体的体积模量，si分别表示含水、含气、含油饱和度，ki对应表示为水、油、气的体积模量；
[0058]
s102：利用常规测井资料，基于s101计算饱和岩石体积模量和流体体积饱和模量，根据自适应基质体积模量反演理论，计算各个测井点处的岩石基质体积模量。
[0059]
在自适应基质体积模量反演理论中，首先需要输入测井曲线v
p
，vs，ρ
sat
，sw，参考图4，判断如若不满足此条件需进行下一测井点的计算，如若满足此条件，可根据s101中所得的饱和岩石体积模量、剪切模量和流体体积模量确定c值和k0的取值范围及初始值，流体因子之差绝对值的初始值b。
[0060]
确定k0的取值范围需引入biot系数表达式k
dry
＝k0(1-β)及岩石基质与岩石骨加参数的关系式，干岩石骨架的统一表达式为参考表1为干岩石骨架统一模型及对应p和q值
[0061]
表1常见干岩石骨架模型及对应p和q值
[0062][0063]
表1中p、q为可调节的孔隙结构参数，m为与干岩石泊松比有关的函数，α为岩石孔隙扁率，c
cond
为固结系数，为临界孔隙度。
[0064]
由于岩石物理学中有关于体积模量的相对关系：
[0065]
k0＞k
sat
＞k
dry
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2-1)
[0066]
由干岩石骨架统一表达式可知在k
sat
》k
dry
两端分别乘可得则有
[0067][0068]
根据式(2-2)，通过调节孔隙结构参数p、q来确定对应深度点下岩石基质的体积模量。将初始值设置为岩石基质体积模量的下限。
[0069]
基于russell流体因子表达式，可得流体因子f1的表达式为：
[0070][0071]
式中，z
p
＝ρ
satvp
表示岩石纵波阻抗，zs＝ρ
satvs
，c为中间变量系数，定义为：
[0072][0073]
式中，σ
dry
为干岩样泊松比，取值范围为0.10～0.40，对应的中间变量c的取值范围在2.25～6.0之间。
[0074]
基于gassmann方程得到流体因子f2的表达式为：
[0075]
[0076]
由于上式的计算中没有干岩石骨架的体积模量，其计算式为k
dry
＝k0(1-β),β为中间变量可根据gassmann-biot-geertsaman方程来获取，其方程可表达为：
[0077][0078]
式中，为纵波模量，为中间参数，σ
dry
为干岩样的泊松比，与russell流体因子中求取c系数的参数一致。由上述改变的岩石基质体积模量和干岩石泊松比，按各自的增量，求取式(2-6)中的β，以此获取k
dry
的值。从而式(2-5)变为gassmann流体因子f2关于岩石基质体积模量和干岩石泊松比的关系式。
[0079]
依据图4自适应基质体积模量反演流程图，将russell流体因子f1与gassmann方程中的流体识别因子f2的绝对值之差b作为反演的目标函数，当改变岩石基质体积模量和干岩石泊松比计算b，判断b是否小于前一次的b值，如果小于，则更新b。不断循环计算和更新b值，直至岩石基质体积模量和干岩石泊松比达到上限值，取得在此范围内最小的b值，这个b值所对应的基质体积模量即为对应测井深度点下的自适应基质体积模量。
[0080]
s103：基于s101得到的孔隙混合流体的流体体积模量、s102得到的饱和岩石体积模量和岩石基质体积模量，确定表示孔隙结构的综合孔隙结构参数。
[0081]
基于gassmann方程代入干岩石骨架的同一表达式可得含孔隙结构参数p、q的孔隙度表达式：
[0082][0083]
令s＝(p q),可得表征孔隙结构参数的综合孔隙结构参数：
[0084][0085]
s104：计算考虑孔隙结构参数影响的岩石biot系数，基于实际地质资料及其测井资料，分析其孔隙结构特征，给出孔隙结构参数的取值范围。
[0086]
由于biot的理论计算表达式为：
[0087][0088]
又有干岩石统一骨架模型为：
[0089][0090]
将(4-2)式代入(4-1)时可得：
[0091][0092]
将s＝(p q),代入式(4-3)，整理得考虑岩石孔隙结构参数的biot系数：
[0093]
[0094]
确定c值的取值范围，由可知，c为孔隙结构参数p、q的函数，其物理意义为孔隙结构参数q在综合孔隙结构参数s＝(p q)中所占有的比重，取值范围在0～1之间。根据表1可知通常将其近似视为0或1，对于岩石孔隙结构参数的选取十分有限，精确度受到一定的影响。加入参数c后，不再近似于1，如表2所示，都有了具体的数值，提高对于岩石骨架模型的计算精度。对于各类岩石具有适用性，可根据研究区域岩性、沉积相的不同，选取不同的c值进行综合孔隙结构参数的研究计算。
[0095]
表2常见干岩石骨架模型及c值
[0096][0097]
s105：利用所计算的考虑孔隙结构参数影响的biot，基于有效应力理论及aplin表达式计算得出地层孔隙流体压力。
[0098]
有效应力理论的表达式如下：
[0099]
σ＝p
d-β
·
p
p
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5-1)
[0100]
式中，σ为有效应力，pd为上覆地层压力，β为虑岩石孔隙度结构参数的biot系数，p
p
为地层孔隙流体压力。
[0101]
上覆地层压力可利用密度积分得出：
[0102][0103]
式中，z为相对应的深度,ρ为饱和岩石密度。
[0104]
aplin表达式如下：
[0105]
v＝v
0-aln(σ/σ0)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5-3)
[0106]
式中，为空隙率，a为压实系数，σ为有效应力。
[0107]
联系有效应力理论及aplin表达式计算基于岩石孔隙结构参数的地层孔隙流体压力的表达式如下：
[0108][0109]
式中，p
p
为地层孔隙流体压力，pd为上覆岩层压力，σ0初始有效压力，v为空隙率，a为地层压实系数，β为考虑岩石孔隙度结构参数的biot系数。
[0110]
当地层孔隙流体压力为正常压力时，此时地层孔隙流体压力为
[0111][0112]
式中，z为相对应的深度,ρf为孔隙流体密度。异常地层孔隙流体压力与正常孔隙流体压力的比值称为压力系数参考图7，在储层区其压力系数明显偏大，可有效区分储层区和非储层区。
[0113]
以上所述的仅是本发明的实施例，方案中公知的具体技术方案或特性等常识在此未作过多描述。应当指出，对于本领域的技术人员来说，在不脱离本发明技术方案的前提下，还可以作出若干变形和改进，这些也应该视为本发明的保护范围，这些都不会影响本发明实施的效果和专利的实用性。本技术要求的保护范围应当以其权利要求的内容为准，说明书中的具体实施方式等记载可以用于解释权利要求的内容。