盐间地层压力预测方法、装置、电子设备及介质与流程

1.本发明涉及油气地球物理技术领域，更具体地，涉及一种盐间地层压力预测方法、装置、电子设备及介质。


背景技术：

2.地层压力又称孔隙压力，是重要的工程甜点参数。近年来，随着国内对页岩油气藏的勘探研究的逐步深入，人们对国内海相页岩气资源的开采有了更深的认识。特别是随着涪陵页岩气藏商业化开采的逐渐深入，人们发现页岩气藏的产量通常与地层压力正相关。准确的地层压力预测结果会为钻井过程中泥浆密度的确定提供重要的信息，合适的泥浆密度配比可以降低泥浆对原状地层的污染，同时降低井喷等事故的发生概率。
3.现今认可度最高的地层孔隙压力计算理论是terzaghi提出的欠压实理论。该理论认为地层孔隙流体压力(pp：pore pressure)等于上覆地层压力(obp：over burdern pressure)及垂直有效应力(ves：vertical effective stress)之差。obp可通过对上覆岩层的密度积分计算得到。因此，孔隙压力计算的核心在于如何获取颗粒间的垂向应力值(ves)。准确的ves值求取是计算pp的核心。
4.ves计算的核心需要构建一条合理的正常压实趋势(nct:normal compaction trend)。砂泥岩地层在沉积过程中，随着沉积物的逐渐压实，泥岩孔隙中的地层流体因沉积压实作用而逐渐被排出，相应地在测井曲线上显示为一条正常压实趋势线，即密度曲线逐渐线性增加，声波曲线逐渐减少，电阻率曲线因矿化度提高而逐渐增加。如果沉积物在压实过程中受到额外的压力引起地层流体不能正常排出，而处于欠压实状态，相应地测井曲线将偏离开正常趋势线。
5.基于欠压实理论，目前对地层孔隙压力的求取主要包括经验系数法、等效深度法、eaton法和bowers法。
6.经验系数法适用于已经有过一定数量地层孔隙压力实测数据的地区。利用该地区已钻井的中途测试、完井试油、rft测试等数据，建立声波时差正常趋势线方程，然后回归出经验系数公式来计算地层孔隙压力。但该方法的局限性在于必须对工区内的地层压力展布有相当的了解，现阶段非常规勘探探区大多面临资料少，开发程度低等问题，因此，该方法由较强的局限性。等效深度法又称平衡深度法，在地层压力预测和检测方面是世界上多数盆地最有效的方法之一。由岩石力学的原理可知，无论在正常压实带还是欠压实带，同样的孔隙度(或其它可反映孔隙度的物性参数)值对应相同的有效应力。平衡深度法认为欠压实带内观察点的孔隙度值之所以得以保存，是因为该点地层在埋深到正常段内的等效深度时地层完全封闭，以后增加的上覆负荷都施加到孔隙流体上。eaton法是国内外油田公司普遍采用的地层孔隙压力计算方法，它具有计算精度高，使用范围广等特点。eaton法主要利用地震速度与垂直有效应力的关系，该方法主要基于正常底层压实条件下的速度变化趋势线来求取地层压力。bowers方法是在确定了ves的基础上，在考虑了因流体膨胀所引起的卸载状态，利用obp来求取pp，实践证实该方法对于不均衡压实产生的高压预测较为准确。
7.以上提及的常规的地层压力方法均基于欠压实理论，而欠压实理论最大的弊端在于需要人为构建nct曲线，为了避免nct曲线所带入的人为影响，一些学者尝试跳出欠压实理论进行压力预测。比较具有代表性的是fillippone针对墨西哥湾地区提出的fillippone法，该方法不需要构建nct，通过计算地层最大和最小的压实速度，进而通过dix公式转化的地震层速度，最终计算地层压力。该方法后续也得到了一定的发展，比较著名的是刘震提出的刘震公式法(1993)和云美厚提出的云美厚公式法(1996)。近年来，孙炜(2015)提出api法，利用地层压力与其他弹性参数之间的对应关系，进行参数筛选，构建适合特定地区的地层压力预测模型。
8.现有方法大多基于欠压实理论，但该方法的弊端在于其预测结果很大程度上依赖正常压实趋势线(nct)的构建，不同的nct曲线的构建会直接影响预测结果。常规的nct趋势线构建方法大多适用于连续沉积的海相地层，而针对陆相盐间地层的nct构建需要进一步探究。虽然有些学者尝试跳出欠压实理论进行压力预测，但是由于压实是异常压力形成的重要原因之一，如果完全不考虑压实的影响，会使得大多方法具有很大的区域局限性。特别是针对压力成因较为复杂的陆相盐间地层，由于其岩性纵向变化异常剧烈，盐岩地层对异常压力有很好的封堵效应，如何准确预测盐间地层的异常压力则面临更多的挑战。因此，有必要开发一种盐间地层压力预测方法、装置、电子设备及介质。
9.公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。


技术实现要素：

10.本发明提出了一种盐间地层压力预测方法、装置、电子设备及介质，其能够通过构建盐间岩石物理模型，获取常压条件的弹性信息，提高了压实趋势线的精度，并结合eaton地层压力预测方法，能够有效提高地层压力预测精度。
11.第一方面，本公开实施例提供了一种盐间地层压力预测方法，包括：
12.构建盐间地层的岩石物理模型；
13.根据所述岩石物理模型，计算盐间地层的弹性张量；
14.根据所述弹性张量，计算每一个深度的正常压实纵波速度；
15.根据每一个深度的正常压实纵波速度，构建正常压实趋势线，计算目标深度的地层压力。
16.优选地，构建盐间地层的岩石物理模型包括：根据vrh平均理论，将盐岩、钙芒硝与石英矿物均匀混合，获得盐间混合物；根据dem理论，将干酪根加入所述盐间混合物中，获得富有机质盐间混合物；计算湿孔隙度，进而计算湿黏土骨架；根据dem理论，将所述湿黏土骨架加入所述富有机质盐岩混合物，获得湿矿物骨架；根据gassmann理论，将干孔隙加入所述湿矿物骨架，获得岩石物理模型。
17.优选地，通过公式(1)计算湿孔隙度：
18.19.其中，为黏土的湿孔隙度，f
c
为考虑孔隙时黏土所占体积百分比，f
cm
为不考虑孔隙时黏土所占矿物的体积百分比，为总孔隙度。
20.优选地，通过公式(2)计算湿黏土骨架：
[0021][0022]
其中，m
wet
为湿黏土骨架的弹性张量，m
clay
为黏土的弹性张量。
[0023]
优选地，通过公式(3)计算岩石物理模型：
[0024][0025]
其中，为孔隙流体的弹性参数，参考表1，m
mix2
为岩石物理模型的弹性张量，m
mix1
为湿矿物骨架的弹性张量。
[0026]
优选地，通过公式(4)计算正常压实纵波速度：
[0027][0028]
其中，v
p_depth
为正常压实纵波速度，k和μ为弹性张量，ρ为对应深度的总密度。
[0029]
优选地，通过公式(5)计算目标深度的地层压力：
[0030][0031]
其中，p
p
为目标深度的待求压力值，p
ov
为上覆地层压力值，p
w
为待求深度的静水压力，v为目标深度的实测速度值。
[0032]
作为本公开实施例的一种具体实现方式，
[0033]
第二方面，本公开实施例还提供了一种盐间地层压力预测装置，包括：
[0034]
构建模块，构建盐间地层的岩石物理模型；
[0035]
弹性张量计算模块，根据所述岩石物理模型，计算盐间地层的弹性张量；
[0036]
纵波速度计算模块，根据所述弹性张量，计算每一个深度的正常压实纵波速度；
[0037]
地层压力计算模块，根据每一个深度的正常压实纵波速度，构建正常压实趋势线，计算目标深度的地层压力。
[0038]
优选地，构建盐间地层的岩石物理模型包括：根据vrh平均理论，将盐岩、钙芒硝与石英矿物均匀混合，获得盐间混合物；根据dem理论，将干酪根加入所述盐间混合物中，获得富有机质盐间混合物；计算湿孔隙度，进而计算湿黏土骨架；根据dem理论，将所述湿黏土骨架加入所述富有机质盐岩混合物，获得湿矿物骨架；根据gassmann理论，将干孔隙加入所述湿矿物骨架，获得岩石物理模型。
[0039]
优选地，通过公式(1)计算湿孔隙度：
[0040]
[0041]
其中，为黏土的湿孔隙度，f
c
为考虑孔隙时黏土所占体积百分比，f
cm
为不考虑孔隙时黏土所占矿物的体积百分比，为总孔隙度。
[0042]
优选地，通过公式(2)计算湿黏土骨架：
[0043][0044]
其中，m
wet
为湿黏土骨架的弹性张量，m
clay
为黏土的弹性张量。
[0045]
优选地，通过公式(3)计算岩石物理模型：
[0046][0047]
其中，为孔隙流体的弹性参数，参考表1，m
mix2
为岩石物理模型的弹性张量，m
mix1
为湿矿物骨架的弹性张量。
[0048]
优选地，通过公式(4)计算正常压实纵波速度：
[0049][0050]
其中，v
p_depth
为正常压实纵波速度，k和μ为弹性张量，ρ为对应深度的总密度。
[0051]
优选地，通过公式(5)计算目标深度的地层压力：
[0052][0053]
其中，p
p
为目标深度的待求压力值，p
ov
为上覆地层压力值，p
w
为待求深度的静水压力，v为目标深度的实测速度值。
[0054]
第三方面，本公开实施例还提供了一种电子设备，该电子设备包括：
[0055]
存储器，存储有可执行指令；
[0056]
处理器，所述处理器运行所述存储器中的所述可执行指令，以实现所述的盐间地层压力预测方法。
[0057]
第四方面，本公开实施例还提供了一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的盐间地层压力预测方法。
[0058]
其有益效果在于：
[0059]
(1)压实趋势线的构建是压力预测的核心。由于盐间地层岩性纵向变化十分剧烈，基于岩石物理建模的压力预测方法考虑了岩性的速度的影响，可以更好地刻画速度随岩性变化的趋势，得到更为准确的压实趋势线，为后续压力预测奠定基础。
[0060]
(2)地下不同深度速度值的变化，可能是由于压力的不同导致，也可能是由岩性的剧烈变化导致，存在多解性，针对盐间地层的特征，这种多解性尤甚。本方法可以使压实趋势线排除岩性对其干扰，使得压力预测更为准确。
[0061]
(3)基于岩石物理的压力预测流程，所有输入参数均基于测井数据，由数据驱动，避免了常规方法人为确定压实趋势线所带来的人为因素干扰。
[0062]
(4)针对盐间地层构建的盐间岩石物理模型可以较好的反映常压条件下盐间地层的弹性信息。
[0063]
本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
[0064]
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
[0065]
图1示出了欠压实理论的示意图。
[0066]
图2示出了根据本发明的一个实施例的盐间地层岩石物理建模流程图。
[0067]
图3示出了根据本发明的一个实施例的盐间地层压力预测方法的步骤的流程图。
[0068]
图4示出了根据本发明的一个实施例的盐间地层的典型地层的示意图。
[0069]
图5a-图5d示出了根据本发明的一个实施例的测井曲线的示意图。
[0070]
图6示出了根据本发明的一个实施例的关于纵波速度的正常压实趋势线的示意图。
[0071]
图7示出了根据本发明的一个实施例的压力预测结果的示意图。
[0072]
图8示出了根据本发明的一个实施例的一种盐间地层压力预测装置的框图。
[0073]
附图标记说明：
[0074]
201、构建模块；202、弹性张量计算模块；203、纵波速度计算模块；204、地层压力计算模块。
具体实施方式
[0075]
下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。
[0076]
图1示出了欠压实理论的示意图。该理论认为上覆地层压力等于地层孔隙流体压力和垂直有效应力之和。其物理意义可表示为，孔隙内流体的压力值等于上覆地层压力之和减去骨架颗粒之间的垂向应力。obp相对比较好计算，利用密度测井曲线，通过对上覆岩层的密度积分即可计算得到。因此，孔隙压力计算的核心在于如何获取颗粒间的垂向应力值(ves)。准确的ves值求取是计算pp的核心。ves可以利用地层在常压下的速度和异常压力下的速度之差来计算得到。异常压力下的地层速度即地下实测得到的速度，可以通过声波时差测井获取，常压下的地层速度则需要基于nct曲线获取。总结来说，地层压力计算的关键在于如何获得ves值，而ves值计算的关键在于nct曲线的确定。
[0077]
图2示出了根据本发明的一个实施例的盐间地层岩石物理建模流程图。
[0078]
本发明提供一种盐间地层压力预测方法，包括：
[0079]
构建盐间地层的岩石物理模型；在一个示例中，构建盐间地层的岩石物理模型包括：根据vrh平均理论，将盐岩、钙芒硝与石英矿物均匀混合，获得盐间混合物；根据dem理论，将干酪根加入盐间混合物中，获得富有机质盐间混合物；将湿孔隙加入黏土，获得湿黏
土；根据dem理论，将湿黏土加入富有机质盐岩混合物，获得湿矿物骨架；根据gassmann理论，将干孔隙加入湿矿物骨架，获得岩石物理模型，如图2所示。
[0080]
vrh平均理论公式为：
[0081][0082][0083][0084]
其中，i表示第i种矿物，f
i
为该矿物的体积含量百分比，各矿物体积百分比随深度变化而变化，m
vrh
为盐间混合物，m
i
为第i种矿物的弹性模量，包括体积模量和剪切模量，具体数值如表1所示，不同地区矿物弹性模量可能会略有区别，具体以本地实验室测试样品为准。
[0085]
表1
[0086][0087][0088]
dem理论为公式(9)，通过公式(9)向盐间混合物中加入干酪根，最终得到富有机质岩盐混合物的弹性张量：
[0089][0090]
其中，k1，μ1为盐间混合物的体积模量与剪切模量，k2，μ2为干酪根的体积模量与剪切模量，k
mix1
，μ
mix1
为盐间混合物的体积模量与剪切模量，y为干酪根的体积百分比，随深度变化，p和q为控制包裹体形状的形状因子。
[0091]
通过公式(1)计算湿孔隙度：
[0092][0093]
其中，为黏土的湿孔隙度，f
c
为考虑孔隙时黏土所占体积百分比，f
cm
为不考虑孔隙时黏土所占矿物的体积百分比，为总孔隙度，可通过测
井曲线直接获取。
[0094]
则通过公式(2)计算湿黏土骨架：
[0095][0096]
其中，m
wet
为湿黏土骨架的弹性张量，m
clay
为黏土的弹性张量。
[0097]
利用dem理论，将湿黏土骨架加入富有机质岩盐混合物，得到湿矿物骨架。
[0098]
通过公式(10)计算干孔隙度：
[0099][0100]
利用公式(3)将干孔隙加入湿矿物骨架，获得岩石物理模型：
[0101][0102]
其中，为孔隙流体的弹性参数，参考表1，m
mix2
为岩石物理模型的弹性张量，m
mix1
为湿矿物骨架的弹性张量。
[0103]
各深度的矿物组分、孔隙度和密度测井曲线输入至岩石物理模型，计算该深度盐间地层的弹性张量。
[0104]
根据盐间地层的弹性张量，计算每一个深度的正常压实纵波速度；在一个示例中，通过公式(4)计算正常压实纵波速度：
[0105][0106]
其中，v
p_depth
为正常压实纵波速度，k和μ为弹性张量，ρ为对应深度的总密度。
[0107]
根据每一个深度的正常压实纵波速度，构建关于纵波速度的正常压实趋势线，计算目标深度的地层压力；在一个示例中，通过公式(5)计算目标深度的地层压力：
[0108][0109]
其中，p
p
为目标深度的待求压力值，p
ov
为上覆地层压力值，p
w
为待求深度的静水压力，v为目标深度的实测速度值。
[0110]
现有技术的正常压实趋势线(nct)，通常是一条随深度缓慢增大的速度曲线，将趋势线上的速度带入公式(5)可计算得到常规方法的压力值。但该方法忽略了岩性变化对于nct趋势的影响，本方法将不同深度岩性变化的影响，带入压力预测，使得新构建的正常压实趋势线，不再是一条近似直线，而是随深度不断变化的曲线，将本方法趋势线上的速度带入公式(5)，即可计算得到目标深度的地层压力值。
[0111]
本发明还提供一种盐间地层压力预测装置，包括：
[0112]
构建模块，构建盐间地层的岩石物理模型；在一个示例中，构建盐间地层的岩石物理模型包括：根据vrh平均理论，将盐岩、钙芒硝与石英矿物均匀混合，获得盐间混合物；根据dem理论，将干酪根加入盐间混合物中，获得富有机质盐间混合物；将湿孔隙加入黏土，获得湿黏土；根据dem理论，将湿黏土加入富有机质盐岩混合物，获得湿矿物骨架；根据
gassmann理论，将干孔隙加入湿矿物骨架，获得岩石物理模型。
[0113]
vrh平均理论公式为：
[0114][0115][0116][0117]
其中，i表示第i种矿物，f
i
为该矿物的体积含量百分比，各矿物体积百分比随深度变化而变化，m
vrh
为盐间混合物，m
i
为第i种矿物的弹性模量，包括体积模量和剪切模量，具体数值如表1所示，不同地区矿物弹性模量可能会略有区别，具体以本地实验室测试样品为准。
[0118]
dem理论为公式(9)，通过公式(9)向盐间混合物中加入干酪根，最终得到富有机质岩盐混合物的弹性张量：
[0119][0120]
其中，k1，μ1为盐间混合物的体积模量与剪切模量，k2，μ2为干酪根的体积模量与剪切模量，k
mix1
，μ
mix1
为盐间混合物的体积模量与剪切模量，y为干酪根的体积百分比，随深度变化，p和q为控制包裹体形状的形状因子。
[0121]
通过公式(1)计算湿孔隙度：
[0122][0123]
其中，为黏土的湿孔隙度，f
c
为考虑孔隙时黏土所占体积百分比，f
cm
为不考虑孔隙时黏土所占矿物的体积百分比，为总孔隙度，可通过测井曲线直接获取。
[0124]
则通过公式(2)计算湿黏土骨架：
[0125][0126]
其中，m
wet
为湿黏土骨架的弹性张量，m
clay
为黏土的弹性张量。
[0127]
利用dem理论，将湿黏土骨架加入富有机质岩盐混合物，得到湿矿物骨架。
[0128]
通过公式(10)计算干孔隙度：
[0129][0130]
利用公式(3)将干孔隙加入湿矿物骨架，获得岩石物理模型：
[0131][0132]
其中，为孔隙流体的弹性参数，参考表1，m
mix2
为岩石物理模型的弹性张量，m
mix1
为湿矿物骨架的弹性张量。
[0133]
弹性张量计算模块，各深度的矿物组分、孔隙度和密度测井曲线输入至岩石物理模型，计算该深度盐间地层的弹性张量。
[0134]
纵波速度计算模块，根据盐间地层的弹性张量，计算每一个深度的正常压实纵波速度；在一个示例中，通过公式(4)计算正常压实纵波速度：
[0135][0136]
其中，v
p_depth
为正常压实纵波速度，k和μ为弹性张量，ρ为对应深度的总密度。
[0137]
地层压力计算模块，根据每一个深度的正常压实纵波速度，构建关于纵波速度的正常压实趋势线，计算目标深度的地层压力；在一个示例中，通过公式(5)计算目标深度的地层压力：
[0138][0139]
其中，p
p
为目标深度的待求压力值，p
ov
为上覆地层压力值，p
w
为待求深度的静水压力，v为目标深度的实测速度值。
[0140]
本发明还提供一种电子设备，电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述的盐间地层压力预测方法。
[0141]
本发明还提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现上述的盐间地层压力预测方法。
[0142]
为便于理解本发明实施例的方案及其效果，以下给出四个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
[0143]
实施例1
[0144]
图3示出了根据本发明的一个实施例的盐间地层压力预测方法的步骤的流程图。
[0145]
如图3所示，该盐间地层压力预测方法包括：步骤101，构建盐间地层的岩石物理模型；步骤102，根据岩石物理模型，计算盐间地层的弹性张量；步骤103，根据弹性张量，计算每一个深度的正常压实纵波速度；步骤104，根据每一个深度的正常压实纵波速度，构建正常压实趋势线，计算目标深度的地层压力。
[0146]
图4示出了根据本发明的一个实施例的盐间地层的典型地层的示意图，由两套盐岩，一套盐间地层所组成，盐岩厚度通常在10-25米左右，盐岩地层孔隙度低，渗透率差，可以起到较好的压力封堵作用。盐间地层通常厚度在5-38米之间，由于有机质热演化等，使得盐间地层孔隙度相对较高，是自生自储的非常规储层。
[0147]
图5a-图5d示出了根据本发明的一个实施例的测井曲线的示意图。由左向右依次
为伽马、纵波速度、密度、孔隙度。盐岩地层的测井相应表现为低gr，高速度，低密度，低孔隙度，而盐间地层则表现为高gr，低速度，中密度，高孔隙度。
[0148]
图6示出了根据本发明的一个实施例的关于纵波速度的正常压实趋势线的示意图。由图可知，虚线为常规方法构建的nct线，灰色实线为本方法计算得到的nct线，而黑色实线为实测的纵波速度。三条曲线的特征十分明显，虚线反映了速度随深度变化的趋势，大体表现为一条直线；而灰色实线与实测的黑色实线则呈现出剧烈的变化趋势，与地质分层结果比对可知，该深度段为岩性变化剧烈的盐间韵律层。因此，实测曲线的周期性抖动很大程度上是由岩性纵向的剧烈变化而导致的。速度是孔隙压力和岩性共同影响的参数，使用单一的nct线显然无法描述岩性变化剧烈的盐韵律层，其只考虑孔隙压力对速度的影响，而忽略了岩性的影响。最终使得孔隙压力预测结果含有岩性的信息，导致最终压力预测结果忽高忽低。然而，利用本方法的盐间岩石物理模型所构建的灰色nct线，随着岩性的变化而变化，其结果包含了纵向岩性的变化信息，从而使得孔隙压力的预测结果更为准确。
[0149]
图7示出了根据本发明的一个实施例的压力预测结果的示意图。空心点为泥浆密度，灰色曲线为常规孔隙压力预测结果，黑色粗实线为新方法的预测结果。对比可知，灰色所示的常规预测结果由于受到了盐间岩性纵向变化剧烈的影响，孔隙压力预测结果呈现周期性抖动，且在许多深度段均超过泥浆密度。泥浆密度侧面反映了地层压力的上限。当孔隙压力大于泥浆密度时，会发生井涌甚至井喷等事故，因此，在压力预测中，通常利用其来衡量预测结果的合理性。基于本方法的预测结果剔除了岩性对压力结果预测的影响，其压力预测值大部分低于泥浆密度值，相较于常规方法，预测结果更为准确合理。
[0150]
实施例2
[0151]
图8示出了根据本发明的一个实施例的一种盐间地层压力预测装置的框图。
[0152]
如图8所示，该盐间地层压力预测装置，包括：
[0153]
构建模块201，构建盐间地层的岩石物理模型；
[0154]
弹性张量计算模块202，根据岩石物理模型，计算盐间地层的弹性张量；
[0155]
纵波速度计算模块203，根据弹性张量，计算每一个深度的正常压实纵波速度；
[0156]
地层压力计算模块204，根据每一个深度的正常压实纵波速度，构建正常压实趋势线，计算目标深度的地层压力。
[0157]
作为可选方案，构建盐间地层的岩石物理模型包括：根据vrh平均理论，将盐岩、钙芒硝与石英矿物均匀混合，获得盐间混合物；根据dem理论，将干酪根加入盐间混合物中，获得富有机质盐间混合物；计算湿孔隙度，进而计算湿黏土骨架；根据dem理论，将湿黏土骨架加入富有机质盐岩混合物，获得湿矿物骨架；根据gassmann理论，将干孔隙加入湿矿物骨架，获得岩石物理模型。
[0158]
作为可选方案，通过公式(1)计算湿孔隙度：
[0159][0160]
其中，为黏土的湿孔隙度，f
c
为考虑孔隙时黏土所占体积百分比，f
cm
为不考虑孔隙时黏土所占矿物的体积百分比，为总孔隙度。
[0161]
作为可选方案，通过公式(2)计算湿黏土骨架：
[0162][0163]
其中，m
wet
为湿黏土骨架的弹性张量，m
clay
为黏土的弹性张量。
[0164]
作为可选方案，通过公式(3)计算岩石物理模型：
[0165][0166]
其中，为孔隙流体的弹性参数，参考表1，m
mix2
为岩石物理模型的弹性张量，m
mix1
为湿矿物骨架的弹性张量。
[0167]
作为可选方案，通过公式(4)计算正常压实纵波速度：
[0168][0169]
其中，v
p_depth
为正常压实纵波速度，k和μ为弹性张量，ρ为对应深度的总密度。
[0170]
作为可选方案，通过公式(5)计算目标深度的地层压力：
[0171][0172]
其中，p
p
为目标深度的待求压力值，p
ov
为上覆地层压力值，p
w
为待求深度的静水压力，v为目标深度的实测速度值。
[0173]
实施例3
[0174]
本公开提供一种电子设备包括，该电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述盐间地层压力预测方法。
[0175]
根据本公开实施例的电子设备包括存储器和处理器。
[0176]
该存储器用于存储非暂时性计算机可读指令。具体地，存储器可以包括一个或多个计算机程序产品，该计算机程序产品可以包括各种形式的计算机可读存储介质，例如易失性存储器和/或非易失性存储器。该易失性存储器例如可以包括随机存取存储器(ram)和/或高速缓冲存储器(cache)等。该非易失性存储器例如可以包括只读存储器(rom)、硬盘、闪存等。
[0177]
该处理器可以是中央处理单元(cpu)或者具有数据处理能力和/或指令执行能力的其它形式的处理单元，并且可以控制电子设备中的其它组件以执行期望的功能。在本公开的一个实施例中，该处理器用于运行该存储器中存储的该计算机可读指令。
[0178]
本领域技术人员应能理解，为了解决如何获得良好用户体验效果的技术问题，本实施例中也可以包括诸如通信总线、接口等公知的结构，这些公知的结构也应包含在本公开的保护范围之内。
[0179]
有关本实施例的详细说明可以参考前述各实施例中的相应说明，在此不再赘述。
[0180]
实施例4
[0181]
本公开实施例提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的盐间地层压力预测方法。
[0182]
根据本公开实施例的计算机可读存储介质，其上存储有非暂时性计算机可读指
令。当该非暂时性计算机可读指令由处理器运行时，执行前述的本公开各实施例方法的全部或部分步骤。
[0183]
上述计算机可读存储介质包括但不限于：光存储介质(例如：cd－rom和dvd)、磁光存储介质(例如：mo)、磁存储介质(例如：磁带或移动硬盘)、具有内置的可重写非易失性存储器的媒体(例如：存储卡)和具有内置rom的媒体(例如：rom盒)。
[0184]
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
[0185]
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。