本发明提供了一种碳酸盐岩地层超压预测方法，其方法包括：针对碳酸盐岩地层选择油气水测井综合解释模型获取基础物性参数；然后利用Voigt‑Reuss‑Hill平均模量模型计算岩石基质体积模量，利用Wood模型或者Patchy模型计算孔隙流体体积模量，利用BISQ模型计算岩石骨架体积模量；最后通过多孔弹性力学量化模型预测超压。本发明的有益效果是：本发明所提供的技术方案针对碳酸盐岩地层提出一种新的超压预测方法，一定程度上解决了碳酸盐岩岩性致密且极不均匀带来的超压预测困难，具有工业应用价值。
1.一种碳酸盐岩地层超压预测方法，其特征在于：包括以下步骤：
S101：根据碳酸盐岩储层物性解释方法，利用测井资料对待预测碳酸盐岩地层的基础物性参数进行解释，获得待预测碳酸盐岩地层的基础物性参数；所述基础物性参数包括：待预测碳酸盐岩地层的岩石泥质矿物含量Vsh、各矿物组成体积含量Vi、岩石基质孔隙度φs和岩石裂缝孔隙度φc；
S102：根据待预测碳酸盐岩地层的岩石泥质矿物含量Vsh和各矿物组成体积含量Vi，采用Voigt-Reuss-Hill平均模量模型，计算得到待预测碳酸盐岩地层的岩石基质体积模量Ks；
S103：采用Wood模型或者Patchy模型，计算得到待预测碳酸盐岩地层的孔隙流体体积模量Kf；
S104：采用BISQ模型，计算得到待预测碳酸盐岩地层的岩石骨架体积模量Kd；
S105：根据岩石孔隙度φ、岩石基质体积模量Ks、孔隙流体体积模量Kf和岩石骨架体积模量Kd，采用多孔弹性力学量化模型，计算获得待预测碳酸盐岩地层的超压预测值p；其中，岩石孔隙度φ为岩石基质孔隙度φs和岩石裂缝孔隙度φc的和；
所述步骤S105中，根据岩石孔隙度φ、岩石基质体积模量Ks、孔隙流体体积模量Kf和岩石骨架体积模量Kd，计算获得待预测碳酸盐岩地层的超压预测值p的计算公式如公式(1)所示：上式中，为平均主应力，为先验值；β为无量纲经验性参数系数，为先验值；A为复合弹性模量组合项；φ＝φs+φc。
2.如权利要求1所述的一种碳酸盐岩地层超压预测方法，其特征在于：步骤S101中，根据碳酸盐岩储层物性解释方法，利用测井资料对待预测碳酸盐岩地层的基础物性参数进行解释，获得待预测碳酸盐岩地层的基础物性参数的步骤，包括：S201：获取碳酸盐岩储层的录井数据和测井数据，所述测井数据包括声波测井数据、密度测井数据、中子测井数据和电阻率测井数据；
S202：利用公式(2)计算碳酸盐岩储层的岩石泥质矿物含量Vsh：
上式中，Vsh为岩石泥质矿物含量，GCUR为地层常数，SH为目的层自然伽马射线强度指数；
S203：根据录井数据中的钻井液滤液电导率，采用阿尔奇公式计算得到碳酸盐岩储层的岩石裂缝孔隙度φc，计算公式如公式(3)所示：上式中，φc为岩石裂缝孔隙度，σLLD和σLLS分别为深侧向电导率和浅侧向电导率，是电阻率测井数据的倒数，σmf和σw分别为钻井液滤液电导率和地层水电导率，为录井数据，mf为裂缝的孔隙度指数，为先验值；
S204：根据岩石泥质矿物含量Vsh和岩石裂缝孔隙度φc，采用宏观岩石体积平衡模型，计算得到碳酸盐岩储层的岩石基质孔隙度φs和碳酸盐岩储层中各矿物组成体积含量Vi，计算公式如公式(4)所示：上式中，φs为岩石基质孔隙度，Vi为第i种矿物体积含量，Δtf、Δtsh、Δtima和Δt分别为流体时差值、泥质时差值、第i种矿物时差值和声波测井数据中的声波时差值；ρf、ρsh、ρima和ρ分别为流体密度值、泥质密度值、第i种矿物密度值和密度测井数据中的密度值；CNLf、CNLsh、CNLima和CNL分别为流体中子值、泥质中子值、第i种矿物中子值和中子测井数据中的中子值；其中Δt、ρ和CNL分别为声波测井、密度测井、中子测井的测井数据，Δtf、Δtsh、Δtima、ρf、ρsh、ρima、CNLf、CNLsh和CNLima为先验值；i＝1,2,3,…,N，N为碳酸盐岩储层中的矿物种类数量；
S205：将岩石基质孔隙度φs、岩石裂缝孔隙度φc和各矿物组成体积含量Vi,作为最终的碳酸盐岩储层物性解释结果。
3.如权利要求2所述的一种碳酸盐岩地层超压预测方法，其特征在于：步骤S202中，SH的计算公式如公式(5)所示：上式中，GRmax和GRmin分别为自然伽马曲线极大值和极小值，GR为含泥质目的层自然伽马读数，为自然伽马测井数据；GRmax、GRmin为先验值。
4.如权利要求2所述的一种碳酸盐岩地层超压预测方法，其特征在于：步骤S202中，地层常数GCUR的值为2。
5.如权利要求2所述的一种碳酸盐岩地层超压预测方法，其特征在于：步骤S204中，碳酸盐岩中矿物含量包括：泥质含量、方解石含量、白云石含量和膏盐含量；在泥质含量已知的情况下，Vi分别表示三种矿物含量：方解石含量V1、白云石含量V2和膏盐含量V3、；宏观岩石体积平衡模型的计算步骤如下：S301：由公式(4)整理得到未知数个数与方程数相等的非奇异线性恰定方程组；所述未知数包括：φs、V1、V2和V3；
S302：通过粒子群算法对所述非奇异线性恰定方程组进行搜索，得到一个接近于解的初值；
S303：根据所述初值，采用非线性约束优化方法对所述非奇异线性恰定方程组进一步求解，得到岩石基质孔隙度φs和V1、V2、V3的准确解。
6.如权利要求1所述的一种碳酸盐岩地层超压预测方法，其特征在于：步骤S102中，根据待预测碳酸盐岩地层的岩石泥质矿物含量Vsh和各矿物组成体积含量Vi，采用Voigt-Reuss-Hill平均模量模型，计算待预测碳酸盐岩地层的岩石基质体积模量Ks的计算公式如公式(6)所示：上式中，Mi为各矿物成分等效模量，为先验值； N为待预测碳
酸盐岩地层中含有的矿物种类数量，其中，f2、f3和f4分别代表灰岩、白云岩和膏盐岩的体积含量。
7.如权利要求1所述的一种碳酸盐岩地层超压预测方法，其特征在于：步骤S103中，计算待预测碳酸盐岩地层的孔隙流体体积模量Kf的计算公式如公式(7)所示：上式中，Ki为待预测碳酸盐岩地层中各混合物流体成分体积模量，为先验值；xi为组成混合物流体的n种组分的各部分体积百分数，计算公式为： m、n、a、b为岩电参数，为先验值；Rt为地层真电阻率曲线，为测井数据；Ri为地层水电阻率，为录井数据；i＝1,
一种碳酸盐岩地层超压预测方法
技术领域
[0001] 本发明涉及地质勘探领域，尤其涉及一种碳酸盐岩地层超压预测方法。
背景技术
[0002] 含油气盆地中的超压成因已有大量的研究和共识，超压可以由各种物理和化学过程产生，欠压实作用与构造侧向挤压等直接改变有效应力；生烃、黏土矿物脱水、烃类裂解等新生孔隙流体体积增大；浮力与流体流动增压作用等。大规模异常高压的成因机理大致可分为三类：欠压实增压(地层压实不平衡)、生烃增压和构造挤压增压等。超压可能造成井喷、井壁失稳、井眼冲蚀和钻井液循环漏失等多种工程安全问题，超压预测研究对钻井工程安全及其成本效益至关重要。地层超压作为自然产能对常规和非常规油气开发有重要影响，也是储层建模和油气藏数值模拟的重要参数。
[0003] 碎屑岩地层超压预测方法主要是在分析测井、地震响应参数的基础上，结合Terzaghi有效应力理论提出的经验性的方法。Taizaghi有效应力理论描述了土壤在上覆应力的作用下的压实过程，对于具有一定压实规律的碎屑岩来说较为适用，但不适用于成岩过程中化学压实作用显著的碳酸盐岩地层。
[0004] 碳酸盐岩地层沉积过程中后生改造强烈(重结晶作用、白云岩化作用、溶蚀作用)，造成原始孔隙结构特征改变显著，岩石成分和岩性也是多重化学反应叠加的结果，这种岩性和物性的多重非均质性使得碳酸盐岩超压预测仍然是国内外尚未解决的研究难题。
发明内容
[0005] 为了解决上述问题，本发明提供了一种碳酸盐岩地层超压预测方法、设备及存储设备，一种碳酸盐岩地层超压预测方法，主要包括以下步骤：
[0006] S101：根据碳酸盐岩储层物性解释方法，利用测井资料对待预测碳酸盐岩地层的基础物性参数进行解释，获得待预测碳酸盐岩地层的基础物性参数；所述基础物性参数包括：待预测碳酸盐岩地层的岩石泥质矿物含量Vsh、各矿物组成体积含量Vi、岩石基质孔隙度φs和岩石裂缝孔隙度φc；
[0007] S102：根据待预测碳酸盐岩地层的岩石泥质矿物含量Vsh和各矿物组成体积含量Vi，采用Voigt-Reuss-Hill平均模量模型，计算得到待预测碳酸盐岩地层的岩石基质体积模量Ks；
[0008] S103：采用Wood模型或者Patchy模型，计算得到待预测碳酸盐岩地层的孔隙流体体积模量Kf；
[0009] S104：采用BISQ模型，计算得到待预测碳酸盐岩地层的岩石骨架体积模量Kd；
[0010] S105：根据岩石孔隙度φ、岩石基质体积模量Ks、孔隙流体体积模量Kf和岩石骨架体积模量Kd，采用多孔弹性力学量化模型，计算获得待预测碳酸盐岩地层的超压预测值p；其中，岩石孔隙度φ为岩石基质孔隙度φs和岩石裂缝孔隙度φc的和。
[0011] 进一步地，步骤S101中，根据碳酸盐岩储层物性解释方法，利用测井资料对待预测碳酸盐岩地层的基础物性参数进行解释，获得待预测碳酸盐岩地层的基础物性参数的步骤，包括：
[0012] S201：获取碳酸盐岩储层的录井数据和测井数据，所述测井数据包括声波测井数据、密度测井数据、中子测井数据和电阻率测井数据；
[0013] S202：利用公式(1)计算碳酸盐岩储层的岩石泥质矿物含量Vsh：
[0014]
[0015] 上式中，Vsh为岩石泥质矿物含量，GCUR为地层常数，SH为目的层自然伽马射线强度指数；
[0016] S203：根据录井数据中的钻井液滤液电导率，采用阿尔奇公式计算得到碳酸盐岩储层的岩石裂缝孔隙度φc，计算公式如公式(2)所示：
[0017]
[0018] 上式中，φc为岩石裂缝孔隙度，σLLD和σLLS分别为深侧向电导率和浅侧向电导率，是电阻率测井数据的倒数，σmf和σw分别为钻井液滤液电导率和地层水电导率，为录井数据，mf为裂缝的孔隙度指数，为先验值；
[0019] S204：根据岩石泥质矿物含量Vsh和岩石裂缝孔隙度φc，采用宏观岩石体积平衡模型，计算得到碳酸盐岩储层的岩石基质孔隙度φs和碳酸盐岩储层中各矿物组成体积含量Vi，计算公式如公式(3)所示：
[0020]
[0021] 上式中，φs为岩石基质孔隙度，Vi为第i种矿物体积含量，Δtf、Δtsh、Δtima和Δt分别为流体时差值、泥质时差值、第i种矿物时差值和声波测井数据中的声波时差值；ρf、ρsh、ρima和ρ分别为流体密度值、泥质密度值、第i种矿物密度值和密度测井数据中的密度值；CNLf、CNLsh、CNLima和CNL分别为流体中子值、泥质中子值、第i种矿物中子值和中子测井数据中的中子值；其中Δt、ρ和CNL分别为声波测井、密度测井、中子测井的测井数据，Δtf、Δtsh、Δtima、ρf、ρsh、ρima、CNLf、CNLsh和CNLima为先验值；i＝1,2,3,…,N，N为碳酸盐岩储层中的矿物种类数量；
[0022] S205：将岩石基质孔隙度φs、岩石裂缝孔隙度φc和各矿物组成体积含量Vi,作为最终的碳酸盐岩储层物性解释结果。
[0023] 进一步地，步骤S202中，SH的计算公式如公式(4)所示：
[0024]
[0025] 上式中，GRmax和GRmin分别为自然伽马曲线极大值和极小值，GR为含泥质目的层自然伽马读数，为自然伽马测井数据；GRmax、GRmin为先验值。
[0026] 进一步地，步骤S202中，地层常数GCUR的值为2。
[0027] 进一步地，步骤S204中，碳酸盐岩中矿物含量包括：泥质含量、方解石含量、白云石含量和膏盐含量；在泥质含量已知的情况下，Vi分别表示三种矿物含量：方解石含量V1、白云石含量V2和膏盐含量V3、；宏观岩石体积平衡模型的计算步骤如下：
[0028] S301：由公式(3)整理得到未知数个数与方程数相等的非奇异线性恰定方程组；所述未知数包括：φs、V1、V2和V3；
[0029] S302：通过粒子群算法搜索一个接近于解的初值；
[0030] S303：根据初值，采用非线性约束优化方法得到岩石基质孔隙度φs和V1、V2、V3的准确解。
[0031] 进一步地，步骤S102中，根据待预测碳酸盐岩地层的岩石泥质矿物含量Vsh和各矿物组成体积含量Vi，采用Voigt-Reuss-Hill平均模量模型，计算待预测碳酸盐岩地层的岩石基质体积模量Ks的计算公式如公式(5)所示：
[0032]
[0033] 上式中，Mi为各矿物成分等效模量，为先验值； i＝1,2,…,N，N为待预测碳酸盐岩地层中含有的矿物种类数量，其中，f2、f3和f4分别代表灰岩、白云岩和膏盐岩的体积含量。
[0034] 进一步地，步骤S103中，计算待预测碳酸盐岩地层的孔隙流体体积模量Kf的计算公式如公式(6)所示：
[0035]
[0036] 上式中，Ki为待预测碳酸盐岩地层中各混合物流体成分体积模量，为先验值；xi为组成混合物流体的n种组分的各部分体积百分数，计算公式为： m、n、a、b为岩电参数，为先验值；Rt为地层真电阻率曲线，为测井数据；Ri为地层水电阻率，为录井数据；i＝1,2,…,n。
[0037] 进一步地，步骤S105中，根据岩石孔隙度φ、岩石基质体积模量Ks、孔隙流体体积模量Kf和岩石骨架体积模量Kd，计算获得待预测碳酸盐岩地层的超压预测值p的计算公式如公式(7)所示：
[0038]
[0039] 上式中，为平均主应力，为先验值；β为无量纲经验性参数系数，为先验值；A为复合弹性模量组合项；φ＝φs+φc。
[0040] 本发明提供的技术方案带来的有益效果是：本发明所提供的技术方案针对碳酸盐岩地层提出一种新的超压预测方法，一定程度上解决了碳酸盐岩岩性致密且极不均匀带来的超压预测困难，具有一定的工业应用价值。
附图说明
[0041] 下面将结合附图及实施例对本发明作进一步说明，附图中：
[0042] 图1是本发明实施例中一种碳酸盐岩地层超压预测方法的流程图；
[0043] 图2是本发明实施例中双庙1井碳酸盐岩层段测井资料超压预测结果的示意图。
具体实施方式
[0044] 为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图详细说明本发明的具体实施方式。
[0045] 本发明的实施例提供了一种碳酸盐岩地层超压预测方法。
[0046] 请参考图1，图1是本发明实施例中一种碳酸盐岩地层超压预测方法的流程图，具体包括如下步骤：
[0047] S101：根据碳酸盐岩储层物性解释方法，利用测井资料对待预测碳酸盐岩地层的基础物性参数进行解释，获得待预测碳酸盐岩地层的基础物性参数；所述基础物性参数包括：待预测碳酸盐岩地层的岩石泥质矿物含量Vsh、各矿物组成体积含量Vi、岩石基质孔隙度φs和岩石裂缝孔隙度φc；
[0048] S102：根据待预测碳酸盐岩地层的岩石泥质矿物含量Vsh和各矿物组成体积含量Vi，采用Voigt-Reuss-Hill平均模量模型，计算得到待预测碳酸盐岩地层的岩石基质体积模量Ks；
[0049] S103：采用Wood模型或者Patchy模型，计算得到待预测碳酸盐岩地层的孔隙流体体积模量Kf；
[0050] S104：采用BISQ模型，计算得到待预测碳酸盐岩地层的岩石骨架体积模量Kd；(具体方法见文献：Dvorkin J.,Nur A.Dynamic poroelasticity:A unified model with the squirt and the Biot mechanisms[J].Geophysics,1993,58(4):524-533.和Mavko G.,Mukerji T.,Dvorkin J.The rock physics handbook:tools for seismic analysis of porous media[M].New York:Cambridge University Press,2009.)；
[0051] S105：根据岩石孔隙度φ、岩石基质体积模量Ks、孔隙流体体积模量Kf和岩石骨架体积模量Kd，采用多孔弹性力学量化模型，计算获得待预测碳酸盐岩地层的超压预测值p；其中，岩石孔隙度φ为岩石基质孔隙度φs和岩石裂缝孔隙度φc的和。
[0052] 步骤S101中，根据碳酸盐岩储层物性解释方法，利用测井资料对待预测碳酸盐岩地层的基础物性参数进行解释，获得待预测碳酸盐岩地层的基础物性参数的步骤，包括：
[0053] S201：获取碳酸盐岩储层的录井数据和测井数据，所述测井数据包括声波测井数据、密度测井数据、中子测井数据和电阻率测井数据；
[0054] S202：利用公式(1)计算碳酸盐岩储层的岩石泥质矿物含量Vsh：
[0055]
[0056] 上式中，Vsh为岩石泥质矿物含量，GCUR为地层常数，SH为目的层自然伽马射线强度指数；
[0057] S203：根据录井数据中的钻井液滤液电导率，采用阿尔奇公式计算得到碳酸盐岩储层的岩石裂缝孔隙度φc，计算公式如公式(2)所示：
[0058]
[0059] 上式中，φc为岩石裂缝孔隙度，σLLD和σLLS分别为深侧向电导率和浅侧向电导率，是电阻率测井数据的倒数，σmf和σw分别为钻井液滤液电导率和地层水电导率，为录井数据，mf为裂缝的孔隙度指数，为先验值；
[0060] S204：根据岩石泥质矿物含量Vsh和岩石裂缝孔隙度φc，采用宏观岩石体积平衡模型，计算得到碳酸盐岩储层的岩石基质孔隙度φs和碳酸盐岩储层中各矿物组成体积含量Vi，计算公式如公式(3)所示：
[0061]
[0062] 上式中，φs为岩石基质孔隙度，Vi为第i种矿物体积含量，Δtf、Δtsh、Δtima和Δt分别为流体时差值、泥质时差值、第i种矿物时差值和声波测井数据中的声波时差值；ρf、ρsh、ρima和ρ分别为流体密度值、泥质密度值、第i种矿物密度值和密度测井数据中的密度值；CNLf、CNLsh、CNLima和CNL分别为流体中子值、泥质中子值、第i种矿物中子值和中子测井数据中的中子值；其中Δt、ρ和CNL分别为声波测井、密度测井、中子测井的测井数据，Δtf、Δtsh、Δtima、ρf、ρsh、ρima、CNLf、CNLsh和CNLima为先验值；i＝1,2,3,…,N，N为碳酸盐岩储层中的矿物种类数量；
[0063] S205：将岩石基质孔隙度φs、岩石裂缝孔隙度φc和各矿物组成体积含量Vi,作为最终的碳酸盐岩储层物性解释结果。
[0064] 步骤S202中，SH的计算公式如公式(4)所示：
[0065]
[0066] 上式中，GRmax和GRmin分别为自然伽马曲线极大值和极小值，GR为含泥质目的层自然伽马读数，为自然伽马测井数据；GRmax、GRmin为先验值。
[0067] 步骤S202中，地层常数GCUR的值为2。
[0068] 步骤S204中，碳酸盐岩中矿物含量包括：泥质含量、方解石含量、白云石含量和膏盐含量；在泥质含量已知的情况下，Vi分别表示三种矿物含量：方解石含量V1、白云石含量V2和膏盐含量V3、；宏观岩石体积平衡模型的计算步骤如下：
[0069] S301：由公式(3)整理得到未知数个数与方程数相等的非奇异线性恰定方程组；所述未知数包括：φs、V1、V2和V3；
[0070] S302：通过粒子群算法搜索一个接近于解的初值；
[0071] S303：根据初值，采用非线性约束优化方法得到岩石基质孔隙度φs和V1、V2、V3的准确解。
[0072] 步骤S102中，根据待预测碳酸盐岩地层的岩石泥质矿物含量Vsh和各矿物组成体积含量Vi，采用Voigt-Reuss-Hill平均模量模型，计算待预测碳酸盐岩地层的岩石基质体积模量Ks的计算公式如公式(5)所示：
[0073]
[0074] 上式中，Mi为各矿物成分等效模量，为先验值； i＝1,2,…,N，N为待预测碳酸盐岩地层中含有的矿物种类数量，其中，f2、f3和f4分别代表灰岩、白云岩和膏盐岩的体积含量。
[0075] 步骤S103中，计算待预测碳酸盐岩地层的孔隙流体体积模量Kf的计算公式如公式(6)所示：
[0076]
[0077] 上式中，Ki为待预测碳酸盐岩地层中各混合物流体成分体积模量，为先验值；xi为组成混合物流体的n种组分的各部分体积百分数，计算公式为： m、n、a、b为岩电参数，为先验值；Rt为地层真电阻率曲线，为测井数据；Ri为地层水电阻率，为录井数据；i＝1,2,…,n。
[0078] 步骤S105中，根据岩石孔隙度φ、岩石基质体积模量Ks、孔隙流体体积模量Kf和岩石骨架体积模量Kd，计算获得待预测碳酸盐岩地层的超压预测值p的计算公式如公式(7)所示：
[0079]
[0080] 上式中， 为平均主应力，为先验值；β为无量纲经验性参数系数，为先验值；A为复合弹性模量组合项；φ＝φs+φc。
[0081] 由于上述模型和参数都是用来表征多孔介质的线弹性行为的，实际岩石往往具有黏弹性(非线弹性性质)和能量耗散。而黏弹性介质中的应力-应变本构关系比经典线弹性介质中的应力-应变本构关系复杂得多，且多个参数通常没有明确的物理意义，很难在实验和工程中测量和应用。因此，由于黏弹性理论唯象性，一般地解决方法是将黏弹性简化为经典弹性力学问题，基于这种思路，碳酸盐岩地层超压预测理论模型需要考虑黏弹性，但可通过引入无量纲经验性参数系数β可简化处理黏弹性复杂问题，进而得到实际用来预测碳酸盐岩地层超压的方程(7)。无量纲经验性参数系数β可理解为与黏弹性有关的校正系数，该无量纲经验性参数系数β可利用岩石物理模拟实验数据、地层测试数据、碳酸盐岩地层沉积相、岩性和物性等资料经验性的获取。
[0082] 本发明实施例的应用研究选取双庙1井，该井位于川东北地区NNE向毛坝-双庙场构造上，为雷口坡组、嘉陵江组、飞仙关组海相碳酸盐岩地层发育超压或强超压的典型钻井，利用录井资料、常规测井(孔隙度，密度，中子，自然伽马，电阻率等)和全阵列声波测井等资料预测和模拟计算碳酸盐岩层段的超压情况，并与钻杆实测压力进行对比。
[0083] 通过提取碳酸盐岩岩性段测井密度、纵横波速度和测井解释岩石物性等测井参数，结合录井数据在不同小层取各参数的平均值，利用Voigt-Reuss-Hill平均模量模型方程和Wood模型方程分别计算岩石基质体积模量和流体体积模量，利用BISQ模型计算岩石骨架体积模量，通过孔隙压力量化模型方程预测孔隙压力(图2)。由图2、表1可看出，预测孔隙压力值与6个钻杆实测深度处的压力值接近，相对误差范围在2％～10％；泥浆密度换算压力和随钻Sigma监测压力变化趋势基本与预测孔隙压力高值变化相一致，可看作预测孔隙压力高值的包络线即泥浆密度压力可平衡地层超压的高范围，说明预测孔隙压力的变化幅度与实际情况较接近。由图2还可看出，预测孔隙压力随深度的变化幅度跳跃明显，蒸发岩发育段超压显著降低，灰岩和白云岩段超压发育，该预测结果可能较好的反映了与双庙1井碳酸盐岩层段多重非均质性相一致的超压频繁变化的特点。
[0084] 表1 双庙1井测井资料预测超压与DST实测压力的误差分析表
[0085]
[0086] 本发明的有益效果是：本发明所提供的技术方案针对碳酸盐岩地层提出一种新的超压预测方法，一定程度上解决了碳酸盐岩岩性致密且极不均匀带来的超压预测困难，具有一定的工业应用价值。
[0087] 以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
