本发明属于油气地球物理勘探领域，尤其涉及一种基于指数约束的地层孔隙压力预测方法及系统。
背景技术：
目前，业内常用的现有技术是这样的：
在油气地球物理勘探中，地层孔隙压力为油气的分布、运移、储集提供了重要信息，它不仅是确定钻井液密度和井深结构的依据，还是决定钻井成败的重要因素。目前常规的地层孔隙压力预测方法可以分为两类：一类是利用测井资料进行地层孔隙压力预测，另一类是运用地震层速度进行地层孔隙压力预测。其中利用测井资料的预测方法是公认的较理想的地层孔隙压力预测手段。
目前，较常用的地层孔隙压力预测方法主要有Eaton模型方法、Fillippone公式方法，如下所示：
Eaton模型方法
Eaton BA(1976.Graphical methodpredicts geoopressure worldwide[J].World Oil,183:100-104)所提出的模型方法是基于正常压实趋势来分析速度场偏差,再根据模拟井建立与孔隙压力数据直接相关的速度扰动经验关系式。
该模型描述了地层孔隙压力Pf与上覆地层压力Pov、静水压力Pw、正常压实速度Vn及地层速度Vi之间的关系，通过测井资料拟合得到参数C，最后将该模型应用于目标段从而得到目标段的地层孔隙压力Pf，但人为建立正常压实趋势线易产生误差，影响预测结果。
Fillippone公式方法
云美厚(地震地层压力预测[J].石油地球物理勘探,1996,31(4)：575～586)所提出的Fillippone公式及其改进在国内广泛用。该方法的主要原理如下：
通过给定或搜索目标范围内的最大速度(即骨架速度)和最小速度(即孔隙流体速度)，结合统计的思维对地层孔隙压力进行预测，其公式如(2)所示：
式中，Pf为地层孔隙压力，单位：MPa；h表示上覆地层深度，单位：m；ρ为上覆地层平均密度，单位：g/cm3；g为重力加速度；vmax为最大层速度，单位：m/s；vmin为最小层速度，单位：m/s；vi为第i层的层速度，单位：m/s，该公式由统计的方法得到，具有较强的区域局限性，适用范围较窄。
综上所述，现有技术存在的问题是：
(1)Eaton模型方法人为建立正常压实趋势线易产生误差，影响预测结果，降低地层压力的预测效率。
(2)Fillippone公式方法由于是从某些特定地区中通过统计的相关方法确定的，其区域性特征较强，综合的区域适应性较差。
(3)在现有技术中，大多都是将Eaton模型方法和Fillippone公式分开讨论，其方法中的优点并未很好的较全面展示。
解决上述技术问题的难度和意义：
本发明为油气勘探开发领域内的地层压力预测提供基于指数约束的地层孔隙压力预测方法，构建了通过指数来约束的地层孔隙压力预测公式(指数约束公式)，将Eaton模型方法和Fillippone公式的假设前提进行统一，减少了油气勘探中人为因素的影响，拓宽应用范围，提高地层压力的预测效率，更好的作用于油气勘探。
技术实现要素：
针对现有技术存在的问题，本发明提供了一种基于指数约束的地层孔隙压力预测方法及系统。
本发明是这样实现的，一种基于指数约束的地层孔隙压力预测方法，所述基于指数约束的地层孔隙压力预测方法为：
基于地震层速度、密度、时深关系和C系数，利用指数约束公式计算获得地层孔隙压力数据体，进行地层孔隙压力的预测。
该方法拓宽了压力预测的应用范围，减少了人为因素的影响，提高了地层孔隙压力的预测效率。其中，指数约束公式如下所示：
式中，vmax为岩石骨架速度，vmin为孔隙流体速度，vi为地层纵波速度，Pf为地层孔隙压力，Pov为上覆地层压力，Pw为静水压力，C为测井拟合参数。
进一步，所述的基于指数约束的地层孔隙压力预测方法，具体包括：
1)根据测井所得数据，选取其中的纵波速度vi，岩石骨架速度vmax，孔隙流体速度vmin，密度ρ、深度h和地层孔隙压力Pf，根据上覆地层压力的计算公式计算得到上覆地层压力，最后再通过指数约束公式计算获得C测井拟合系数；其中，h表示上覆地层深度，ρ(h)为深度为h时所对应的密度，g为重力加速度；
2)采用基于三维叠后地震资料和叠加速度数据反演得到的地震层速度、密度；
3)利用三维叠后地震资料与测井资料得到对应的时深关系；
4)将步骤1)中得到的C系数代入到公式
中，结合步骤2)中三维地震数据与测井资料反演得到的速度、密度和步骤3)得到的时深关系利用指数约束公式计算出地层压力数据体，进行地层孔隙压力预测。
本发明的另一目的在于提供一种实现所述基于指数约束的地层孔隙压力预测方法的计算机程序。
本发明的另一目的在于提供一种实现所述基于指数约束的地层孔隙压力预测方法的信息数据处理终端。
本发明的另一目的在于提供一种计算机可读存储介质，包括指令，当其在计算机上运行时，使得计算机执行所述的基于指数约束的地层孔隙压力预测方法。
本发明的另一目的在于提供一种基于指数约束的地层孔隙压力预测系统，包括：
C测井拟合系数获取模块，用于根据测井所得数据，选取其中的纵波速度vi，岩石骨架速度vmax，孔隙流体速度vmin，密度ρ、深度h和地层孔隙压力Pf，根据上覆地层压力的计算公式计算得到上覆地层压力，最后再通过指数约束公式计算获得C测井拟合系数；其中，h表示上覆地层深度，ρ(h)为深度为h时所对应的密度，g为重力加速度；
地震层速度、密度参数获取模块，用于通过基于三维叠后地震资料和叠加速度数据反演得到的地震层速度、密度参数；
时深关系获取模块，用于通过三维叠后地震资料与测井资料得到对应的时深关系；
地层孔隙压力预测模块，用于通过C测井拟合系数获取模块得到的C系数代入到公式
中，结合地震层速度、密度参数获取模块获得的的速度、密度参数和时深关系获取模块得到的时深关系，利用指数约束公式计算出地层压力数据体，进行地层孔隙压力预测。
本发明的另一目的在于提供一种搭载有所述基于指数约束的地层孔隙压力预测系统的信息数据处理终端。
综上所述，本发明的优点及积极效果为：
本发明克服现有技术的不足，提供了基于指数约束的地层孔隙压力预测方法，该方法统一了Fillippone公式和Eaton模型方法的假设前提，通过指数约束的方式构建了一种新的地层孔隙压力预测公式(指数约束公式)来进行地层孔隙压力预测，省去了建立正常压实趋势线的步骤，解决了Fillippone公式方法中区域适应性较差的问题，兼具有Fillippone公式方法和Eaton模型方法的优点，拓宽了应用范围，提高了地层孔隙压力预测结果的准确性。
相比于现有技术的缺点和不足，本发明具有以下有益效果：
(1)本发明提供的地层压力预测方法，省去建立正常压实趋势线的步骤，减少了人为因素的影响，降低了理论误差，提高了地层孔隙压力预测效率；
(2)本发明提供的地层压力预测方法有效地解决了Fillippone公式方法区域适应性较差的问题，拓宽了地层孔隙压力预测的应用范围。
(3)本发明提供的地层压力预测方法兼具Fillippone公式方法和Eaton模型方法的优点，克服了Eaton模型方法和Fillippone公式假设前提不一致的问题，具有较强的区域适用性，提高了地层孔隙压力预测结果的准确性。
附图说明
图1是本发明实施例提供的基于指数约束的地层孔隙压力预测方法流程图。
图2是本发明实施例提供的基于指数约束的地层孔隙压力预测系统示意图。
图中：1、C测井拟合系数获取模块；2、地震层速度、密度参数获取模块；3、时深关系获取模块；4、地层孔隙压力预测模块。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
本发明涉及的专业术语：
(1)上覆岩层压力：是指计算点所在深度以上地层中的岩石基质和孔隙流体的重量所产生的压力，在计算海底地层上覆岩层压力时还需考虑海水产生的压力。
(2)地层孔隙压力：指地层孔隙和缝洞中的流体(水、油、气)所具有的压力。
(3)静水压力：由地层水的重量所引起的压力，代表着正常的地层孔隙压力。
(4)正常压实趋势线法：利用测井资料检测异常地层孔隙压力的传统方法，主要一句“泥质沉积物不平衡压实造成地层欠压实并产生异常高压”这一最普遍的机制。
(5)岩石骨架速度：岩层有效孔隙度近于零时的声速。
(6)孔隙流体速度：岩层刚性近于零时的声速。
本发明实施例提供的基于指数约束的地层孔隙压力预测方法为：
基于地震层速度、密度、时深关系和C系数，利用指数约束公式计算获得地层孔隙压力数据体，进行地层孔隙压力的预测；其中，指数约束公式如下所示：
式中，vmax为岩石骨架速度，vmin为孔隙流体速度，vi为地层纵波速度，Pf为地层孔隙压力，Pov为上覆地层压力，Pw为静水压力，C为测井拟合参数。
图1，本发明实施例提供的基于指数约束的地层孔隙压力预测方法具体包括：
S101：根据测井所得数据，选取其中的纵波速度vi，岩石骨架速度vmax，孔隙流体速度vmin，密度ρ、深度h和地层孔隙压力Pf，根据上覆地层压力的计算公式计算得到上覆地层压力，最后再通过指数约束公式计算获得C测井拟合系数；其中，h表示上覆地层深度，ρ(h)为深度为h时所对应的密度，g为重力加速度；
S102：采用基于三维叠后地震资料和叠加速度数据反演得到的地震层速度、密度；
S103：利用三维叠后地震资料与测井资料得到对应的时深关系；
S104：将步骤S101中得到的C系数代入到公式
中，结合步骤S102中三维地震数据与测井资料反演得到的速度、密度和步骤S103得到的时深关系利用指数约束公式计算出地层压力数据体，进行地层孔隙压力预测。
如图2所示，本发明提供一种基于指数约束的地层孔隙压力预测系统，包括：
C测井拟合系数获取模块1，用于根据测井所得数据，选取其中的纵波速度vi，岩石骨架速度vmax，孔隙流体速度vmin，密度ρ、深度h和地层孔隙压力Pf，根据上覆地层压力的计算公式计算得到上覆地层压力，最后再通过指数约束公式计算获得C测井拟合系数；其中，h表示上覆地层深度，ρ(h)为深度为h时所对应的密度，g为重力加速度；
地震层速度、密度参数获取模块2，用于通过基于三维叠后地震资料和叠加速度数据反演得到的地震层速度、密度参数；
时深关系获取模块3，用于通过三维叠后地震资料与测井资料得到对应的时深关系；
地层孔隙压力预测模块4，用于通过C测井拟合系数获取模块得到的C系数代入到公式
中，结合地震层速度、密度参数获取模块获得的的速度、密度参数和时深关系获取模块得到的时深关系，利用指数约束公式计算出地层压力数据体，进行地层孔隙压力预测。
下面结合具体分析对本发明作进一步描述。
1、指数约束公式的推导
(1)常规的地层孔隙压力预测公式如下：
Pf＝Pov-σ (3)
式中，Pf为地层孔隙压力，Pov为上覆地层压力，σ为垂直有效应力。
(2)首先对公式(4)的Fillippone公式进行分析：
式中，Pf为地层孔隙压力，单位：MPa；h表示上覆地层深度，单位：m；ρ为上覆地层平均密度，单位：g/cm3；g为重力加速度；vmax为最大层速度，单位：m/s；vmin为最小层速度，单位：m/s；vi为第i层的层速度，单位：m/s。结合上式(3)和(4)可知，该方法的有效应力与速度成线性关系，其关系如下所示：
σ＝a+bV (5)
式中V为速度，σ为垂直有效应力，a、b为参数。该线性关系使得Fillippone公式的应用范围变得局限，不能较好的应用于所有区域。
(3)再对公式(6)的Eaton模型方法进行分析：
式中Pf为地层孔隙压力，Pcv为上覆地层压力、Pw为静水压力、Vn为正常压实速度，Vi为地层速度，C为测井拟合参数，结合上式(3)和(6)可知，该模型方法的假设条件是有效应力与速度成指数关系即公式(7)：
σ＝bVc (7)
式中V为速度，σ为垂直有效应力，b为参数，C为测井拟合参数。
(4)综合上述公式(5)和公式(7)可见，Fillippone公式的假设条件与Eaton模型方法的假设条件不一致，且Fillippone公式的应用范围较窄，针对以上问题本发明将二者的假设条件进行统一，将Fillippone公式的假设条件修改成指数型关系，得到如下关系式：
σ＝a+bVc (8)
其中V为速度，σ为垂直有效应力，a、b为参数，C为测井拟合参数。
结合公式(4)和(5)的关系，根据公式(8)对公式(5)进行修改，得到Fillippone公式的改进公式：
式中，Pf为地层孔隙压力，单位：MPa；h表示上覆地层深度，单位：m；ρ为上覆地层平均密度，单位：g/cm3；g为重力加速度；vmax为最大层速度，单位：m/s；vmin为最小层速度，单位：m/s；vi为第i层的层速度，单位：m/s，C为测井拟合参数。
(5)根据地层正常压实情况下两者满足关系PEaton＝PFillippone且静水压力相等，综合Eaton模型方法和Fillippone公式的改进公式以及正常压实下的速度关系Vi＝Vn，由Fillippone公式的改进公式(9)得到其正常压实速度Vn，如下所示：
(6)再将公式(10)的Vn代入公式(6)中进行组合，得到新的地层压力预测公式(指数约束公式)：
其中vmax为岩石骨架速度，vmin为孔隙流体速度，vi为地层纵波速度，Pf为地层孔隙压力，Pov为上覆地层压力，Pw为静水压力，C为测井拟合参数。
2、基于指数约束公式的地层孔隙压力预测流程
(1)根据测井所得数据，选取其中的纵波速度vi，岩石骨架速度vmax，孔隙流体速度vmin，密度ρ、深度h和地层孔隙压力Pf，根据上覆地层压力的计算公式计算得到上覆地层压力，最后再通过指数约束公式计算获得C测井拟合系数；其中，h表示上覆地层深度，ρ(h)为深度为h时所对应的密度，g为重力加速度。
(2)采用基于三维叠后地震资料和叠加速度数据反演得到的地震层速度、密度；
(3)利用三维叠后地震资料与测井资料得到对应的时深关系；
(4)将流程(1)中得到的C系数代入到公式(11)中，结合流程(2)中三维地震数据与测井资料反演得到的速度、密度和流程(3)得到的时深关系利用指数约束公式计算出地层压力数据体，实现地层孔隙压力预测。
本发明方法省去了建立正常压实趋势线的步骤，解决了Fillippone公式方法中区域适应性较差的问题，兼具有Fillippone公式方法和Eaton模型方法的优点，拓宽了应用范围，提高了地层孔隙压力预测结果的准确性。
在上述实施例中，可以全部或部分地通过软件、硬件、固件或者其任意组合来实现。当使用全部或部分地以计算机程序产品的形式实现，所述计算机程序产品包括一个或多个计算机指令。在计算机上加载或执行所述计算机程序指令时，全部或部分地产生按照本发明实施例所述的流程或功能。所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。所述计算机指令可以存储在计算机可读存储介质中，或者从一个计算机可读存储介质向另一个计算机可读存储介质传输，例如，所述计算机指令可以从一个网站站点、计算机、服务器或数据中心通过有线(例如同轴电缆、光纤、数字用户线(DSL)或无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器或数据中心进行传输)。所述计算机可读取存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质，(例如，软盘、硬盘、磁带)、光介质(例如，DVD)、或者半导体介质(例如固态硬盘Solid StateDisk(SSD))等。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。