基于岩石物理学参数计算地层压力的方法技术

一种基于岩石物理学参数计算地层压力的方法，利用地层的纵波速度VP、横波速度VS、密度ρ和孔隙度φ测井数据，并按如下步骤：利用Gassmann方程推导出干岩体积模量Kdry和岩石基质物质的体积模量Kma的比值；根据V-G-H模型计算矿物成分的有效体积模量Kma，Wood模型计算孔隙流体的体积模量Kf，求取有效应力系数α和Skempton系数B；根据Skempton系数B和多孔介质有效应力并且根据有效应力系数α和Skempton系数B估算地层压力。本方法适用于混合岩性和单一岩性，不依赖于传统的泥岩或页岩的欠压实理论来进行地层的压力估算，计算得到的地层压力更逼近于地层的真实压力。

【技术实现步骤摘要】

本专利技术涉及勘探地球物理资料处理
，特别是涉及到岩石物理学中一种基于岩石物理学参数计算地层压力的方法。
技术介绍
目前地层压力(地层孔隙流体压力)预测的应用主要分为两个方面，一是利用地震速度对地层压力进行钻前预测，为钻井工程服务，主要是对地层可能存在的异常高压进行预测，防止钻井事故发生，钻井泥浆比重一般介于地层的正常压力和地层破裂压力之间。二是利用测井资料对已钻井进行地层压力估算，进一步预测地层的应力、破裂压力和坍塌压力，以便对储层进行压裂，为石油勘探开发服务。除此之外，一些地质学家们期望精确的地层压力预测能为油气聚集成藏研究提供一些有利的证据，以便进行油气成藏的分析。精确的地层压力预测一直是石油勘探开发中所面临的一大难题，由于地层压力的成因复杂，不同地区的压力特征也有所不同。以泥岩或页岩(以下简称泥页岩)欠压实和等效深度法为理论基础的地层压力预测方法长期以来在地层压力预测研究中占主导地位。关于地层异常高压力的成因和表现方面，一般认为泥页岩的欠压实作用，由于孔隙中的流体不能排出，孔隙体积维持压实前的状态(比正常压实状态下的孔隙体积大)，此时泥页岩纵波速度会比正常压实情况下的纵波速度低，在这种情况下泥页岩能导致地层产生异常高压；而在实际钻井中，地层的异常高压都表现在可渗透性地层(如砂岩)发育层段。根据实际资料分析，即使是欠压实的泥页岩的孔隙度也是非常低的，孔喉半径十分的微小，毛细管压力非常的大，孔隙中的流体以束缚水的状态存在，同时泥页岩的渗透率是非常低的(除泥页岩裂缝外)，其内部和外部的流体无法通过泥页岩，所以在钻井过程中钻遇到泥页岩层段地层一般不会表现出异常高压。在上述这种情况下，一直存在争论，有的人认为此时泥页岩不存在高压，而有的人认为此时泥页岩存在高压，只是泥页岩中的流体太少，渗透率极低，所以表现不出高压。关于此争论目前还没有权威来定论谁是正确的，谁是不正确的。李传亮(2009)指出泥岩地层是否存在异常高压从来没有被证实过，因为其地层压力无法实测；并且指出等效深度法并不等效。Shaker(2002)认为目前地层压力的术语可能被混淆了，预测的地层压力是那些相对不渗透岩石(页岩和泥岩)的地层压力，而实际测量的地层压力却是储层特性岩石(砂岩)的地层压力。在很多情况下，砂岩孔隙中的压力与页岩夹层的压力没有直接的关系，地层孔隙压力的预测存在疑点和误区。Zoback(2011)认为有两种情况需要根据地球物理数据预测地层压力：第一种情况是钻前通过地震反射数据预测孔隙压力，在可能出现的高压地区进行钻井安全设计需要孔隙压力剖面；第二种情况是估算页岩中的孔隙压力，由于渗透性较差难以直接测量，即使在钻后也需要检测。在这两种情况下，在某些区域应用效果较好的技术手段有可能在其他区域变得不适用。并且指出，在存在较大挤压应力或者沉积和初始压实后孔隙压力增加较大的情况下，不能利用简单压实曲线对孔隙压力进行预测。一个以泥页岩欠压实或等效深度法为理论基础的地层压力预测经验公式很难揭示地层压力的本质，这一直是地层压力预测技术研究的一个瓶颈，由此可见地层压力的研究目前还存在一些疑点需要解决，寻找地层压力的本质才是最根本的途径。根据石油勘探实践发现，地层异常高压力与油气聚集有着非常密切的关系，简单来讲地层孔隙中的流体受到挤压或膨胀作用，如果流体不能发生转移，那么就意味着会产生异常高压，这是一个遵循石油地质学油气运移理论的原理。从岩石物理学的角度来讲，如果砂岩储层中存在异常高压的流体，那么岩石物理学参数(如泊松比、体积模量、剪切模量等岩石物理参数)与正常压力状态下的岩石物理参数有什么区别呢？一直以来，大家公认为纵波速度(如声波测井速度、VSP测井速度和地震速度)是预测地层压力的最佳参数，即VP=K+4/3μρ]]>其中VP是岩石的纵波速度，K是岩石的体积模量，μ是岩石的剪切模量，ρ是岩石的密度。可以看出，地层压力其实是地层岩石物理学参数的一个综合表现结果，这就提出一个问题，我们在预测地层压力时，实质上在预测什么，要想获得更加精确的地层压力预测结果，有没有更好的途径。
技术实现思路
本专利技术的目的是针对目前地层压力预测的难点和疑点，如何从地层压力的力学特性来展开岩石物理学方面的研究将会是一个好的出发点，地层压力的研究要建立在地质背景的分析和岩石力学理论的基础上，结合地球物理勘探方法来对地层压力进行表征，提出了一种基于岩石物理学参数计算地层压力的方法。该方法根据Skempton系数B和多孔介质有效应力原理推导出一种新的地层压力计算方法，通过求取有效应力系数α和Skempton系数B估算地层压力。本专利技术的目的通过如下技术措施来实现。基于岩石物理学参数计算地层压力(地层孔隙流体压力)的方法，利用地层的纵波速度VP、横波速度VS、密度ρ和孔隙度φ测井数据，并按如下步骤：步骤1，利用Gassmann方程推导出干岩体积模量Kdry和岩石基质物质(颗粒)的体积模量Kma的比值的表达式；步骤2，根据V-G-H模型计算矿物(颗粒)成分的有效体积模量Kma，Wood模型计算孔隙流体的体积模量Kf，求取之值，由此求取有效应力系数α和Skempton系数B；步骤3，根据Skempton系数B和多孔介质有效应力原理推导出一种新的地层压力计算方法，并且根据有效应力系数α和Skempton系数B估算地层压力。上述方案进一步包括：在步骤1中，利用Gassmann方程推导出干岩体积模量Kdry和岩石基质物质(颗粒)的体积模量Kma的比值的表达式：KdryKma=Kwet/Kma[1+φ(Kma/Kf-1)]-1Kwet/Kma+1+φ(Kma/Kf-1)-2---(1-1)]]>式中：Kdry为干岩体积模量，Kwet为饱和岩石的体积模量(湿岩)，Kma为岩石基质物质(颗粒)的体积模量，Kf为孔隙流体的体积模量，φ为岩石的孔隙度。在步骤1中，具体推导步骤如下：由Gassmann方程得ρVP2=Kdry+43μdry+(1-Kdry/Kma)2(1-φ-Kdry/Kma)/Kma+φ/Kf---(1-2)]]>同时由波动方程得VP=Kwet+4/3μdryρ---(1-3)]]>VS=μdryρ---(1-4)]]>根据(1—2)、(1—3)两式得到岩石在不排水情况下的体积模量Kwet=Kdry+(1-Kdry/Kma)2本文档来自技高网...

【技术保护点】
基于岩石物理学参数计算地层压力的方法，其特征在于利用地层的纵波速度VP、横波速度VS、密度ρ和孔隙度φ测井数据，并按如下步骤：步骤1，利用Gassmann方程推导出干岩体积模量Kdry和岩石基质物质的体积模量Kma的比值的表达式；步骤2，根据V‑G‑H模型计算矿物成分的有效体积模量Kma，Wood模型计算孔隙流体的体积模量Kf，求取之值，由此求取有效应力系数α和Skempton系数B；步骤3，根据Skempton系数B和多孔介质有效应力原理推导出地层压力计算方法，并且根据有效应力系数α和Skempton系数B估算地层压力。

【技术特征摘要】
1.基于岩石物理学参数计算地层压力的方法，其特征在于利用地层的纵波
速度VP、横波速度VS、密度ρ和孔隙度φ测井数据，并按如下步骤：
步骤1，利用Gassmann方程推导出干岩体积模量Kdry和岩石基质物质的体积
模量Kma的比值的表达式；
步骤2，根据V-G-H模型计算矿物成分的有效体积模量Kma，Wood模型计算
孔隙流体的体积模量Kf，求取之值，由此求取有效应力系数α和Skempton
系数B；
步骤3，根据Skempton系数B和多孔介质有效应力原理推导出地层压力计
算方法，并且根据有效应力系数α和Skempton系数B估算地层压力。
2.根据权利要求1所述的基于岩石物理学参数计算地层压力的方法，其特
征在于，在步骤1中，利用Gassmann方程推导出干岩体积模量Kdry和岩石基质
物质的体积模量Kma的比值的表达式：
KdryKma=Kwet/Kma[1+φ(Kma/Kf-1)]-1Kwet/Kma+1+φ(Kma/Kf-1)-2]]>式中：Kdry为干岩体积模量，Kwet为饱和岩石的体积模量，Kma为岩石基质物
质的体积模量，Kf为孔隙流体的体积模量，φ为岩石的孔隙度。
3.根据权利要求1或2所述的基于岩石物理学参数计算地层压力的方法，
其特征在于，在步骤2中，已知饱和岩石的体积模量Kwet,根据V-G-H模型...

【专利技术属性】
技术研发人员：王鸿升，贺洋洋，王兴谋，闫昭岷，冯德永，宫红波，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司胜利油田分公司物探研究院，
类型：发明
国别省市：北京;11