一种基于卷积神经网络和Eaton公式的地层孔隙压力计算方法

摘要

一种基于卷积神经网络和Eaton公式的地层孔隙压力计算方法，包括：测井曲线重叠采样；基于短时傅里叶变换的样本1D‑2D变换预处理，将一维的深度域测井曲线样本转换为二维的深频图；基于卷积神经网络的正常压实层段智能识别模型，采用数据驱动的方式，对测井曲线进行分段特征提取，识别正常压实层段；根据识别出的正常压实层段测井曲线，构建正常压实趋势线拟合方程；根据正常压实趋势线方程，采用Eaton公式进行地层孔隙压力剖面计算。本发明采用基于数据驱动和物理模型相结合的地层孔隙计算方法，避免了正常压实趋势线构建过程中存在的人为主观性，提高了地层孔隙压力的计算精度。

权利要求

1.一种基于卷积神经网络和Eaton公式的地层孔隙压力计算方法，其特征在于，包括：

1）选取已完钻井测井资料进行模型训练：对用于计算地层孔隙压力的测井曲线参数进行重叠采样，并分别进行短时傅里叶变换，进而得到各采样样本对应的深频图；

2）将步骤1）的深频图输入卷积神经网络进行训练，实现对正常压实层段的自适应特征提取，得到正常压实层段智能识别模型；

3）将目标井对应用于计算地层孔隙压力的测井曲线参数按照步骤1）进行重叠采样和短时傅里叶变换，得到各样本对应的深频图；将本步骤的样本深频图输入步骤2）中的正常压实层段智能识别模型，以识别最优正常压实层段；

4）拟合构建正常压实趋势线方程，并带入Eaton公式；

5）计算待测目标井的地层孔隙压力剖面。

2.根据权利要求1所述的一种基于卷积神经网络和Eaton公式的地层孔隙压力计算方法，其特征在于，在步骤4）中，拟合构建正常压实趋势线方程为：其中， 为正常压实层速度，km/s； 为正常压实电阻率，Ω； 为正常压实视密度，g/cm3； 为正常压实声波时差，us/ft；为井深，m； 为常数； 时地表层速度，km/s； 时地表电阻率，Ω； 时地表视密度，g/cm3；

时地表声波时差，us/ft。

3.根据权利要求2所述的一种基于卷积神经网络和Eaton公式的地层孔隙压力计算方法，其特征在于，所述步骤5）计算待测目标井的地层孔隙压力剖面包括：将得到的所述公式（II）带入Eaton公式，对全井段地层孔隙压力剖面进行预测，其中，所述Eaton公式为：其中， 为地层孔隙压力，g/cm3； 为上覆岩层压力，g/cm3； 为静液柱压力，g/cm3；

为Eaton指数，一般取1.6； 为地层层速度，km/s，地层层速度可以通过钻前地震勘探地震波反射时间求得，也可以通过声波时差测井求得。

4.根据权利要求1所述的一种基于卷积神经网络和Eaton公式的地层孔隙压力计算方法，其特征在于，所述步骤1）中所述重叠采样为：对测井曲线根据步长k进行重叠采样。

5.根据权利要求1所述的一种基于卷积神经网络和Eaton公式的地层孔隙压力计算方法，其特征在于，所述步骤1）中的短时傅里叶变换为基于短时傅里叶变换的样本1D-2D变换预处理：对采样样本根据井深利用窗函数进行截取，并对截取的测井曲线进行傅里叶变换，将每段截取样本得到的频谱根据深度轴进行排列，最终得到具有时变特性的二维深频图，其中，短时傅里叶变换公式：式中， 为测井曲线， 为中心位于井深 处的窗函数， 表示傅里叶变换的频域参数， 表示测井曲线的深度域参数。

6.根据权利要求1所述的一种基于卷积神经网络和Eaton公式的地层孔隙压力计算方法，其特征在于，步骤2）中所述训练得到正常压实层段智能识别模型的方法为：对已完钻井的各深频图对应的压实类型进行标定，将标定压实类型的时频图集合作为训练集输入卷积神经网络中进行学习。

7.根据权利要求1所述的一种基于卷积神经网络和Eaton公式的地层孔隙压力计算方法，其特征在于，对所述卷积神经网络采用反向传播算法不断更新网络参数。

说明书

一种基于卷积神经网络和Eaton公式的地层孔隙压力计算

方法

技术领域

[0001] 本发明涉及一种基于卷积神经网络和Eaton公式的地层孔隙压力计算方法，属于油气钻井的技术领域。

背景技术

[0002] 地层孔隙压力是研究井壁稳定、进行钻井工程设计和评价储层的重要参数之一，有效、准确地进行地层孔隙压力预测，对于优快钻井、保护油气储层具有重要意义。Eaton法是目前最为常用的地层孔隙压力计算方法，该方法利用地震或测井资料根据正常压实趋势线进行地层孔隙压力计算，具有连续性好、分辨率高、实用性强等特点。正常压实趋势线的构建是利用Eaton法进行地层孔隙压力计算的前提，根据沉积压实理论，泥岩欠压实地层的孔隙度比正常压实情况偏大，测井过程中表现为参数曲线偏离正常压实趋势线。现有方法通过主观设定正常压实区间并拟合得到一条正常压实趋势线，然后基于正常压实趋势线计算纵向上沿整个井深范围内的地层孔隙压力。由于钻井地质环境的复杂性、地震测井等相关解释资料的模糊性和人为判断的主观性等，使得构建的正常压实趋势线存在不确定性，最终导致地层孔隙压力预测结果存在不同程度的误差。

[0003] 近年来，人工智能和机器学习技术在智能钻井领域的应用日益完善，基于数据驱动的分析方法为技术人员在钻井设计、地质特征参数预测等工程领域的分析研究带来了新的思路。卷积神经网络是一种有监督的深度学习算法，被广泛的应用于模式识别领域。卷积神经网络不同于其他机器学习算法，它具有局部感受野、权值共享以及池化的特点，通过卷积层对图像进行特征提取，权值共享与池化不仅降低了网络模型的复杂度，也减少了过拟合的风险，从而可以构建处理海量数据的深度学习框架。现有专用专利主要针对时间域的信号通过时频变换进行卷积神经网络模式识别，然而钻井技术领域的深度域的地震、测井数据资料，无法体现卷积神经网络处理大数据方面的优势。

[0004] 因此，有必要针对钻井技术领域地层孔隙压力的计算，考虑采样规则和信号处理方式，结合卷积神经网络和Eaton公式，建立基于数据驱动和物理模型相结合的地层孔隙压力计算方法。

发明内容

[0005] 针对现有技术的不足，本发明公开一种基于卷积神经网络和Eaton公式的地层孔隙压力计算方法。

[0006] 发明概述：

[0007] 本发明针对地层孔隙压力计算中正常压实趋势线构建过程存在的不确定性，以及深度域的测井资料采样数据量不足，通过重叠采样的方式进行采样并对采样样本进行短时傅里叶变换，将测井曲线采样样本进行1D-2D变换预处理，转换为具有深度域、频域信息的二维深频图；然后将深频图采样样本输入卷积神经网络进行地层压实的特征自适应学习，建立正常压实层段智能识别模型，实现端到端的正常压实层段识别；最后将训练好的模型应用于目标井识别正常压实层段，并拟合出正常压实趋势线方程，带入Eaton公式进行纵向上沿整个井深范围内的地层孔隙压力计算。采用基于数据驱动和物理模型相结合的地层孔隙计算方法，避免了正常压实趋势线构建过程中存在的人为主观性，提高了地层孔隙压力的计算精度。

[0008] 发明具体的技术方案如下：

[0009] 一种基于卷积神经网络和Eaton公式的地层孔隙压力计算方法，其特征在于，包括：

[0010] 1）选取已完钻井测井资料进行模型训练：对用于计算地层孔隙压力的测井曲线参数进行重叠采样，并分别进行短时傅里叶变换，进而得到各采样样本对应的深频图；其中，对用于计算地层孔隙压力的测井曲线参数为已完钻井的层速度等可用于地层孔隙压力计算的测井资料，包括但不限于：地震层速度、声波时差、电阻率和视密度；

[0011] 2）将步骤1）的深频图输入卷积神经网络进行训练，实现对正常压实层段的自适应特征提取，得到正常压实层段智能识别模型；

[0012] 3）将目标井对应用于计算地层孔隙压力的测井曲线参数按照步骤1）进行重叠采样和短时傅里叶变换，得到各样本对应的深频图；将本步骤的样本深频图输入步骤2）中的正常压实层段智能识别模型，以识别最优正常压实层段；

[0013] 4）拟合构建正常压实趋势线方程，并带入Eaton公式；

[0014] 5）计算待测目标井的地层孔隙压力剖面。

[0015] 根据本发明优选的，在步骤4）中，拟合构建正常压实趋势线方程为：

[0016]

[0017] 其中， 为正常压实层速度，km/s； 为正常压实电阻率，Ω； 为正常压实视密度，g/cm3； 为正常压实声波时差，us/ft； 为井深，m； 为常数；时地表层速度，km/s； 时地表电阻率，Ω； 时地表视密度，g/
cm3； 时地表声波时差，us/ft。

[0018] 本发明优选的，对某一口地层孔隙压力待算目标井的层速度曲线进行采样，为了减少计算量可选取井段上部地层进行采样，根据现场经验，新近系等上部地层往往表现正常压实。将各样本进行1D-2D变换预处理，得到各待识别样本的深频图，然后带入建立的正常压实层段智能识别模型，得到各样本的识别结果，将正常压实吻合度最高的样本作为最终的正常压实层段。

[0019] 根据本发明优选的，所述步骤5）计算待测目标井的地层孔隙压力剖面包括：

[0020] 将得到的所述公式（II）带入Eaton公式，对全井段地层孔隙压力剖面进行预测，其中，所述Eaton公式为：

[0021]

[0022] 其中， 为地层孔隙压力，g/cm3； 为上覆岩层压力，g/cm3； 为静液柱压力，g/cm3； 为Eaton指数，一般取1.6； 为地层层速度，km/s，地层层速度可以通过钻前地震勘探地震波反射时间求得，也可以通过声波时差测井求得。

[0023] 根据本发明优选的，所述步骤1）中所述重叠采样为：对测井曲线根据步长k进行重叠采样。其中，样本长度可根据经验或计算能力选取，假设长度为L，此处的测井曲线包括层速度曲线、声波时差曲线、电阻率曲线和视密度曲线，上述测井曲线均适用于采用本发明进行正常压实趋势线的构建。

[0024] 根据本发明优选的，所述步骤1）中的短时傅里叶变换为基于短时傅里叶变换的样本1D-2D变换预处理：对采样样本根据井深利用窗函数进行截取，并对截取的测井曲线进行傅里叶变换，将每段截取样本得到的频谱根据深度轴进行排列，最终得到具有时变特性的二维深频图，其中，短时傅里叶变换公式：

[0025]

[0026] 式中， 为测井曲线， 为中心位于井深 处的窗函数， 表示傅里叶变换的频域参数，表示测井曲线的深度域参数。

[0027] 根据本发明优选的，步骤2）中所述训练得到正常压实层段智能识别模型的方法为：

[0028] 对已完钻井的各深频图对应的压实类型进行标定，将标定压实类型的深频图集合作为训练集输入卷积神经网络中进行学习。其中，进行标定的深频图属于已完钻井，地层压实类型可根据井史资料进行标定；所述卷积神经网络优选LetNet-5等经典卷积神经网络结构模型，假设最优卷积神经网络结构如图2所示。

[0029] 根据本发明优选的，对所述卷积神经网络采用反向传播算法不断更新网络参数。

[0030] 本发明的技术优势在于：

[0031] 本发明通过短时傅里叶变换和卷积深井网络，建立了数据驱动的正常压实层段智能识别模型，避免了传统地层孔隙压力计算过程中正常压实层段人为选择存在的主观性和不确定性，提高了地层孔隙压力计算精度。

附图说明

[0032] 图1为本发明技术方案流程图。

[0033] 图2为本发明设定最优卷积神经网络结构图。

[0034] 图3为A井层速度曲线图。

[0035] 图4为A井采样样本图。

[0036] 图5为A井①号样本时频图。

[0037] 图6为A井地层孔隙压力剖面图。

具体实施方式

[0038] 下面结合实施案例和说明书附图对本发明做详细的说明，但不限于此。

[0039] 实施案例：

[0040] 以某区块A井为例，以A井为地层孔隙压力待算目标井，采用A井的10口已完钻邻井作为建模训练样本，选取A井深度小于1500m的新近系上部地层作为正常压实层段采样区间进行测试样本采样，其层速度曲线如图3所示，该图数据是地震或测井工具测得的原始数据。

[0041] 一种基于卷积神经网络和Eaton公式的地层孔隙压力计算方法，包括：

[0042] 1）选取已完钻井测井资料进行模型训练：对用于计算地层孔隙压力的测井曲线参数进行重叠采样，并分别进行短时傅里叶变换，进而得到各采样样本对应的深频图；

[0043] 其中，所述重叠采样是对A井及10口已完钻邻井的层速度曲线进行重叠采样，采样滑动步长为10m，样本长度为200m，采样频率为10hz。其中，目标A井采样样本为测试集，10口邻井采样样本为训练集，对10口邻井的采样样本根据井史资料对地层压实情况进行标定。A井采样情况如图4所示，图4为A井正常压实层段采样区间采样示意图。

[0044] 所述步骤1）中的短时傅里叶变换为基于短时傅里叶变换的样本1D-2D变换预处理：对采样样本根据井深利用窗函数进行截取，并对截取的测井曲线进行傅里叶变换，将每段截取样本得到的频谱根据深度轴进行排列，最终得到具有时变特性的二维深频图，其中，短时傅里叶变换公式：

[0045]

[0046] 式中， 为测井曲线， 为中心位于井深处的窗函数， 表示傅里叶变换的频域参数，表示测井曲线的深度域参数。

[0047] 对A井和邻井的采样样本进行短时傅里叶变换，窗函数为Hamming窗，窗长设置为512，窗口每次平移步长为1，设置频率刻度为0.5:0.5:5.5。最终得到各采样样本的时频图，其中A井①号样本时频图如图5所示。此处参数的选取需要经过多次试验确定，实验发现，其中窗口长度越长深度分辨率越差，频率分辨率越好；窗长越短深度分辨率越好，频率分辨率越差。窗口平移步长越小，深频图图像深度精度越好，但计算量增大。频率刻度根据采样频率确定，一般最大刻度为采样频率的1/2。

[0048] 2）将步骤1）的深频图输入卷积神经网络进行训练，实现对正常压实层段的自适应特征提取，得到正常压实层段智能识别模型；

[0049] 其中，所述训练得到正常压实层段智能识别模型的方法为：

[0050] 对已完钻井的各深频图对应的压实类型进行标定，将标定压实类型的深频图集合作为训练集输入卷积神经网络中进行学习。其中，进行标定的深频图属于已完钻井，地层压实类型可根据井史资料进行标定；所述卷积神经网络优选LetNet-5等经典卷积神经网络结构模型，本实施案例选择最优卷积神经网络结构如图2所示，选择最优神经网络结构各层参数如表1所示，为一个7层的CNN模型：In-C1-P1-C2-P2

[0051] -F-Out。将本步骤得到的标定压实类型的邻井采样样本深频图作为训练集，选取损失函数为对数损失函数：

[0052]

[0053] 其中，y为输出变量，x为输入变量，L为损失函数，N为输入样本量，M为可能的类别数， 为一个二值指标，表示类别j是否是输入实例xi的真实类别， 为模型预测输入xi属于类别j的概率。

[0054] 将反向传播算法中的随机梯度下降算法作为优化算法对网络参数进行更新学习，最终建立基于卷积神经网络的正常压实层段智能识别模型。

[0055] 表1 最优卷积神经网络结构

[0056]

[0057] 3）将目标井对应用于计算地层孔隙压力的测井曲线参数按照步骤1）进行重叠采样和短时傅里叶变换，得到各样本对应的深频图；将本步骤的样本深频图输入步骤2）中的正常压实层段智能识别模型，以识别最优正常压实层段；

[0058] 4）拟合构建正常压实趋势线方程，并带入Eaton公式；

[0059] 在步骤4）中，拟合构建正常压实趋势线方程为：

[0060]

[0061] 其中， 为正常压实层速度，km/s； 为正常压实电阻率，Ω； 为正常压实视密度，g/cm3； 为正常压实声波时差，us/ft； 为井深，m； 为常数； 为时地表层速度，km/s； 为 时地表电阻率，Ω； 为 时地表视密度，g/cm3； 时地表声波时差，us/ft。

[0062] 具体的，将按照步骤1）采样方式得到的A井采样样本的深频图作为测试集带入步骤2）建立的正常压实层段智能识别模型，得到各样本的识别结果，如表2所示。可知①号样本对应井段正常压实吻合度最高，选取①号样本井段935 1135m作为A井正常压实井段。根~据该井段层速度曲线进行正常压实趋势线方程拟合，得到拟合方程为：

[0063]

[0064] 其中，K1、K2根据正常压实层段层速度与深度数据，采样多项式拟合方法求出。

[0065] 表2 A井正常压实层段识别结果

[0066] 样本号 ① ② ③ ④ ⑤正常压实 0.852 0.733 0.531 0.422 0.574
异常压实 0.148 0.267 0.469 0.578 0.426

[0067] 5）计算待测目标井的地层孔隙压力剖面，包括：

[0068] 将得到的所述公式（II）带入Eaton公式，对全井段地层孔隙压力剖面进行预测，其中，所述Eaton公式为：

[0069]

[0070] 其中， 为地层孔隙压力，g/cm3； 为上覆岩层压力，g/cm3； 为静液柱压力，g/3
cm ； 为Eaton指数，一般取1.6； 为地层层速度，km/s，地层层速度可以通过钻前地震勘探地震波反射时间求得，也可以通过声波时差测井求得。

[0071] 具体的，根据步骤4）得到的A井正常压实趋势线方程，计算A井层速度测点对应井深的正常压实层速度，然后带入Eaton公式计算各测点的地层孔隙压力，Eaton公式为：

[0072]

[0073] 最终得到A井地层孔隙压力剖面，如图6所示。通过对比地层孔隙压力实测点可知，采用本发明计算得到的地层孔隙压力剖面与实测点具有很好的贴合度，相对于传统经验方法主观设定正常压实区间提高了地层孔隙压力的计算精度。

[0074] 另外，对所述卷积神经网络采用反向传播算法不断更新网络参数。