基于弹性阻抗的未固结地层渗透率预测方法及处理终端与流程

1.本发明涉及渗透率预测技术领域，具体涉及基于弹性阻抗的未固结地层渗透率预测方法及处理终端。


背景技术：

2.利用地震数据对渗透率预测研究长期以来是国际研究的前沿和难点，现有技术中，主要有经验法、地震属性分析法和地质统计学模拟法来分析研究渗透率。这些方法都存在其局限和不足，例如，经验法过于依赖于研究人员自身的经验，效果不稳定；地质统计学模拟法通过拟合方式实现渗透率预测，但由于地层储藏影响因素众多，使得其拟合预测精度受限，预测结果不尽如人意。在油气藏的渗透率预测方面，油藏地球物理学者甘利灯等人通过采集储层段中子孔隙度、弹性参数及岩石骨架柔度因子进行岩相分析实现渗透率预测，但这类渗透率预测方法是针对固结地层和较深地层。
3.总体来说，现有水合物储层渗透率预测存在以下不足：
4.a.还未形成针对未固结水合物地层的渗透率定量预测方法。现有技术中，针对水合物储层针对性较弱，由于水合物储层大多富集在浅软地层中，导致常规油气储层的渗透率预测方法不能直接应用到水合物储层的渗透率上，且现有油气储层的渗透率预测大多集中在定型方面的研究，定量研究极少。
5.b.现有渗透率预测中，对渗透率敏感的敏感性参数主要基于avo反演这种间接方式进行计算，导致预测稳定性差。
6.c.现有采用拟合法建立起渗透率的敏感性参数和渗透率之间关系的误差较大，计算结果的准确性还需要提高。


技术实现要素：

7.针对现有技术的不足，本发明的目的之一提供一种基于弹性阻抗的未固结地层渗透率预测方法，其能够解决根据地震数据未固结地层预测的问题；
8.本发明的目的之二提供一种处理终端，其能够解决根据地震数据未固结地层预测的问题。
9.实现本发明的目的之一的技术方案为：一种基于弹性阻抗的未固结地层渗透率预测方法，包括如下步骤：
10.步骤1：按公式
①
计算出水合物储层的流体模量
[0011][0012]
式中，
[0013]
公式计算得到的k
f
作为其中，ei表示水合物储层的弹性阻抗，ei0表示入射角度θ为0时的弹性阻抗，θ表示入射角度，k
f0
表示k
f
的平均值φ表示中子孔隙度，φ表示中子孔隙度，φ0表示φ的平均值，f
m
表示固体刚度系数，f
m0
表示f
m
的平均值，ρ表示密度，ρ0表示ρ的平均值，
[0014]
基于碎屑岩饱和岩石的固解液解耦方程计算出水合物储层的流体模量上限
[0015]
步骤2：按公式
②
计算水合物储层的流体模量k
′
f
：
[0016][0017]
式中，c1和c2均为正数的权重系数，且c1 c2＝1；
[0018]
步骤3：按公式
③
计算出水合物储层的双状态模量差pci:
[0019][0020]
其中，c为常数，μ表示剪切模量，φ
c
表示临界孔隙度，双状态模量差pci表征水合物储层为具有可压缩性且处于未固结状态的岩石，
[0021]
通过中子孔隙度反演计算得到中子孔隙度；
[0022]
步骤4：将双状态模量差pci和中子孔隙度作为输入参数，输入预测神经网络进行处理，预测神经网络的输出即为水合物储层的渗透率。
[0023]
进一步地，水合物储层的弹性阻抗ei通过avo地震反射特征方程推演得到，avo地震反射特征方程如下：
[0024]
[0025]
实现本发明的目的之二的技术方案为：一种处理终端，其包括：
[0026]
存储器，用于存储程序指令；
[0027]
处理器，用于运行所述程序指令，以执行所述基于弹性阻抗的未固结地层渗透率预测方法的步骤。
[0028]
本发明的有益效果为：本发明能够很好地对作为未固结地层的水合物储层的渗透率进行预测，且能够进行定量计算，预测结果稳定可靠，实用性强。
附图说明
[0029]
图1为较佳实施例的原理流程示意图；
[0030]
图2是与其他弹性参数(中子孔隙度、纵波速度)交汇分析示意图；
[0031]
图3是计算出的流体模量k
′
f
和真实流体模量的对比示意图；
[0032]
图4是本实施例得到的流体模量上下限以及最终整体的流体模量和真实中的流体模量上下限及整体的流体模量对比示意图；
[0033]
图5是计算得到的渗透率和实际测量的渗透率对比示意图；
[0034]
图6是中国南海某测线的海底浅层渗透率预测结果示意图；
[0035]
图7为处理终端的示意图。
具体实施方式
[0036]
下面，结合附图以及具体实施方案，对本发明做进一步描述。
[0037]
如图1～图6所示，一种基于弹性阻抗的未固结地层渗透率预测方法，包括如下步骤：
[0038]
步骤1：构建对水合物储层渗透率敏感的弹性参数pci。
[0039]
在地球物理研究中，水合物储层(即水合物沉积层)可以看作是一种岩石，饱和岩石的体积模量k
sat
、干岩石的体积模量k
dry
和岩石骨架的体积模量k0具有如下关系：
[0040][0041]
式中，φ表示中子孔隙度，k
φ
表示孔隙刚度，也即表征有效干岩石孔隙可压缩性，其定义为常孔隙压力下岩石体积对内部应力的导数，
[0042]
另外，若岩石处于较软状态，根据reuss公式(即公式
②
)可计算岩石的平均体积模量k
r
：
[0043][0044]
k
r
为岩石的体积模量下限，即岩石处于颗粒完全悬浮状态。对比公式
①
和
②
，饱和
岩石与软岩石的两种状态岩石的体积模量的差别体现在k
φ
′
和上。其中，本实施例的软岩石是指具有可压缩性且处于未固结状态的岩石，水合物储层即属于软岩石，软岩石也即是未固结岩石。据此，将上述两种状态岩石的体积模量的差值作为水合物储层的双状态模量差pci，经推导，该双状态模量差等于k
φ
，也即是有效干岩石孔隙可压缩性。因此，有公式
③
的如下关系：
[0045]
pci＝k
sat
‑
k
r
∞＝k
φ
‑‑‑‑‑‑③
[0046]
如果水合物储层正常压实，根据han and batzle等人研究结果和临界中子孔隙度模型，流体模量k
f
用公式
④
表示：
[0047][0048]
式中，f表示gassmann流体因子，φ
c
表临界中子孔隙度，为常数，本实施例取0.7，c为biot弹性常数，v
p
表示纵波速度，v
s
表示横波速度，ρ表示密度，μ表示剪切模量。
[0049]
如果水合物储层处于完全松软，岩石的平均体积模量k
r
和流体模量k
f
具有如公式
⑤
的关系：
[0050][0051]
根据公式
④
和
⑤
，可以得到水合物储层的双状态模量差pci的计算公式
⑥
：
[0052][0053]
其中，k
f
只与孔隙内的流体类型和流体饱和度有关，φ
c
表示临界孔隙度。
[0054]
参考图2，图2是弹性参数pci与其他弹性参数(中子孔隙度、纵波速度)交汇分析示意图，从该图可以很好地看出，弹性参数pci对高渗透水合物储层具有很好的识别效果，中子孔隙度对高渗透水合物储层也同样具有很好的敏感性，因此，可以利用弹性参数pci和中子孔隙度来预测水合物储层的渗透率。
[0055]
步骤2：进行反演得到弹性参数pci和计算中子孔隙度。
[0056]
真实的水合物储层的流体模量介于状态一和状态二之间，状态一是指水合物储层完全松软，此时水合物储层的流体模量为下限状态二是指水合物固结成岩，此时水合物储层的流体模量为上限基于此假设，可利用公式
④
和
⑤
分别计算得到流体模量下限和流体模量上限也即是根据公式
④
反演计算得到的k
f
作为流体模量下限根据公式
⑤
反演计算得到的k
f
作为流体模量上限进而根据公式
⑦
计算得到真实的水合物储层的流体模量k
′
f
：
[0057][0058]
式中，c1和c2均为正数的权重系数，且c1 c2＝1，本实施例，分别取值为0.63和0.37。
[0059]
为了计算上述两种状态下各自的水合物储层的流体模量，可通过地震反演方法得到。
[0060]
参考图3，图3是按公式
⑦
计算出的流体模量k
′
f
和真实流体模量的对比示意图，两者几乎一致，从而说明根据本实施例是科学有效的。
[0061]
水合物储层在完全松软条件下，水合物储层的avo地震反射特征方程如公式
⑧
所示：
[0062][0063]
式中，r(θ)表示反射系数，θ表示入射角度，f
m
表示固体刚度系数。
[0064]
需要说明的是，公式
⑧
中的v
p
、v
s
、ρ也分别是指纵波速度、横波速度和密度，但公式
⑧
是指松软条件下的纵波速度、横波速度和密度。
[0065]
为实现弹性阻抗反演，需要在公式
⑧
基础上进一步推导弹性阻抗方程，为此，用波阻抗的对数值表示反射系数r(θ)，有：
[0066][0067]
ei表示水合物储层的弹性阻抗，根据数学知识可知，可将
δln(x)
替换为因此可以得到下式：
[0068][0069]
对上式取积分后指数化处理，可得下式：
[0070][0071]
上式在入射角度θ变化时，弹性阻抗ei量纲变化剧烈，因此需要进行标准化处理，引入ei0，ei0＝v
p
ρ，ei0表示入射角度θ为0时的弹性阻抗，从而得到公式
⑨
：
[0072][0073]
式中，k
f0
,φ0,f
m0
,ρ0依次为k
f
,φ,f
m
,ρ的平均值。在实际使用中，获得的测井数据往往有多个，每一个测井数据有对应的一个k
f
,φ,f
m
,ρ，因此，k
f0
,φ0,f
m0
,ρ0也即是所有这些k
f
,φ,f
m
,ρ的平均值。
[0074]
直接从弹性阻抗ei提取地震反演所需的参数，可以减小如现有技术中通过avo反演等间接计算代理的累积误差。为此，对公式
⑨
等式两边求对数，得到不同入射角度θ的弹性阻抗的多个方程，例如，四个入射角度(θ1,θ2,θ3,θ4)有公式
⑩
：
[0075][0076]
式中，
[0077]
式中，可以通过反演出的井旁弹性阻抗以及实测的测井曲线可以计算出四个系数(a(θ),b(θ),c(θ),d(θ))，从而可以计算出公式
⑩
中的16个系数，进而可以求解出任一个采样点下的k
f
,φ,f
m
,ρ，此处得到的k
f
为
[0078]
根据yin和zhang于2014年给出的碎屑岩饱和岩石的固解液解耦方程(yin x y,zhang s x.bayesian inversion for effective pore
‑
fluid bulk modulus based on fluid
‑
matrix decoupled amplitude variation with offset approximation.geophysics,2014,79(5):r221
‑
r232)，也即是正常压实状态下的地震发射特征方程可求得流体模量上限从而根据求得的和计算得到真实的水合物储层的流体模量k
′
f
，进而根据公式
⑥
计算出pci。
[0079]
参考图4，图4是按本实施例的方法进行反演得到的流体模量上下限以及最终整体的流体模量和真实中的流体模量上下限及整体的流体模量对比示意图，从图中可以看出，反演的结果与实际得到的结果几乎一致，再次说明按本实施例进行反演是可行的。
[0080]
至于水合物储层的中子孔隙度，可通过中子孔隙度反演计算得到，中子孔隙度反演为现有技术，本实施例采用基于阻抗反射系数与中子孔隙度反射系数的关系进行反演。
[0081]
步骤3：将得到的pci水合物储层的中子孔隙度作为输入参数，输入预测神经网络进行处理，预测神经网络的输出即为水合物储层的渗透率，也即得到预测的水合物储层渗透率。
[0082]
其中，预测神经网络包括分别由正弦隐藏神经元和线性输出神经元的两层前馈网络构成，预测神经网络通过利用测井数据作为样本进行训练，输入为pci水合物储层和中子孔隙度，输出为水合物储层渗透率。由于pci水合物储层、中子孔隙度与水合物储层渗透率之间存在密切的关系，预测神经网络能够快速找出输入与输出之间的映射关系，从而可以预测出水合物储层渗透率。
[0083]
参考图5和图6，图5是采用本实施例的预测神经网络计算得到的水合物储层渗透
率和实际测量得到的水合物储层渗透率对比示意图，从该图可看出，预测结果和预测结果吻合。图6是中国南海某测线的海底浅层渗透率预测结果示意图，图6的纵坐标表示时间/毫秒，其中已标注出水合物储层发育位置，水合物位于气烟囱上方，由于水合物富集导致渗透率明显降低，而与之临近的块体流沉积底部高渗透率较高，这与实际的地质规律符合，说明本实施例的方法可以应用在实际应用中，具有很好的实用性。
[0084]
其中，本实施例的预测神经网络采用包括隐藏层和输出层的两层神经网络结构。
[0085]
本实施例提供的渗透率预测方法可以作为一项高技术服务，用于为海洋、测绘、地质勘探等目标客户提供高端服务，也即为目标客户提供在水合物物储层的渗透率预测专业技术服务，为需要水合物储层渗透率预测提供技术是的预测方案。同时，也可以作为海洋工程中勘探等自用，可作为海洋工程中工程勘察船、海底资源调查船等上的技术终端，以实现未固结地层中的渗透率预测。
[0086]
如图7所示，本发明还涉及一种处理终端100，其包括：
[0087]
存储器101，用于存储程序指令；
[0088]
处理器102，用于运行所述程序指令，以执行所述基于弹性阻抗的未固结地层渗透率预测方法的步骤。
[0089]
本说明书所公开的实施例只是对本发明单方面特征的一个例证，本发明的保护范围不限于此实施例，其他任何功能等效的实施例均落入本发明的保护范围内。对于本领域的技术人员来说，可根据以上描述的技术方案以及构思，做出其它各种相应的改变以及变形，而所有的这些改变以及变形都应该属于本发明权利要求的保护范围之内。