基于理论引导神经网络的破裂压力预测方法技术

本发明专利技术公开了基于理论引导神经网络的破裂压力预测方法，包括：收集已钻井的测井数据，并进行测井解释得到地质力学参数和破裂压力数据；将数据集划分为训练集和测试集；构建基于理论引导的神经网络模型；利用训练集的数据对基于理论引导的神经网络模型进行训练，并对神经网络模型进行优化，获得优化后的神经网络模型；根据测试集的数据对优化后的神经网络模型预测效果进行测试，并通过评价指标进行模型预测结果评估。本发明专利技术采用理论引导的神经网络模型进行直井破裂压力的预测；该模型在数据驱动的基础上增加了理论知识的约束，进一步提高了模型的预测精度并增强了模型的可解释性，对于钻井工程、水力压裂的设计与施工等具有十分重要的作用。重要的作用。重要的作用。

【技术实现步骤摘要】
基于理论引导神经网络的破裂压力预测方法


[0001]本专利技术涉及基于理论引导神经网络的破裂压力预测方法，属于石油勘探开发


技术介绍

[0002]地层破裂压力是当井内钻井液柱所产生的压力升高至足以压裂地层，使其原有的裂缝张开延伸或形成新的裂隙时的井内流体压力，通常用破裂压力p
f
或当量钻井液密度ρ
f
来表示。在石油工程领域，破裂压力是钻井、固井、水力压裂设计与施工的重要基础参数。随着我国将油气勘探开发的重点转移到非常规油气，所面临的地层环境越来越复杂，破裂压力的研究愈发受到关注。
[0003]目前，获取破裂压力的方法主要有通过矿场试验直接测试破裂压力和通过地球物理测井解释预测破裂压力两大类。矿场试验主要有地层完整性测试和地破试验，这种方法主要用于测试每一开次裸眼井段顶部的破裂压力，得到的破裂压力精度较高但费用高、实测数据少；测井解释方法先通过测井资料预测地层岩石力学、孔隙压力和地应力等地质力学参数，然后将地质力学参数带入到破裂压力计算模型，从而可以预测出破裂压力剖面，这种方法费用较低，而且可以获得相对连续的单井剖面，但精度较低、普适性不强。
[0004]近年来随着机器学习方法在各个领域的广泛应用，许多学者利用机器学习方法解决石油工程领域中的诸多问题。其中，神经网络因其能够有效构建输入参数与输出参数之间的非线性关系，实现准确预测而被广泛使用。而传统的人工神经网络是一种纯粹的数据驱动模型，其基本规则是一系列线性运算加非线性激活函数的组合，如果有足够多的样本数据便可以充分学习到数据的多样性和复杂性，从而推理出未知数据的预测值。但是，传统的神经网络模型需要大量的数据进行训练才能学习到样本数据的多样性，有限的样本数量往往无法使模型得到理想的结果；另外，模型的可解释性较差，神经网络通常由许多层神经元组成，这些层之间的复杂交互很难解释，因此被称为“黑盒模型”。
[0005]目前，许多学者将理论知识、机理模型或函数关系嵌入神经网络模型中，利用这种理论与数据相融合的新范式打破模型效果提升的瓶颈，并取得了良好的效果。但是将理论知识结合神经网络预测破裂压力的研究尚未见报道。因此，专利技术了一种基于理论引导神经网络的破裂压力预测方法，该方法相较于数据驱动神经网络具有更高的精度和更强的可解释性。

技术实现思路

[0006]为了克服现有技术中存在的缺陷，本专利技术旨在提供基于理论引导神经网络的破裂压力预测方法，该方法可根据测井数据直接、准确地预测破裂压力。
[0007]本专利技术解决上述技术问题所提供的技术方案是：基于理论引导神经网络的破裂压力预测方法，包括：
[0008]步骤S1、收集已钻井的测井数据，并进行测井解释得到地质力学参数和破裂压力
数据，构建神经网络模型训练和测试所需的数据集；
[0009]步骤S2、对数据集的数据进行预处理，并将数据集划分为训练集和测试集；
[0010]步骤S3、构建基于理论引导的神经网络模型；
[0011]所述基于理论引导的神经网络模型包含数据驱动模块和理论引导模块，所述数据驱动模块包括输入层、隐藏层、输出层；
[0012]步骤S4、利用训练集的数据对基于理论引导的神经网络模型进行训练，并对神经网络模型进行优化，获得优化后的神经网络模型；
[0013]步骤S5、根据测试集的数据对优化后的神经网络模型预测效果进行测试，并通过评价指标进行模型预测结果评估。
[0014]进一步的技术方案是，所述测井数据包括井深、横波时差、纵波时差、补偿中子、岩性密度、井径、自然伽马；所述地质力学参数包括最小水平地应力、最大水平地应力、孔隙压力、Biot系数和抗拉强度。
[0015]进一步的技术方案是，所述步骤S1中通过一维地质力学模型计算得到地质力学参数，通过直井破裂压力计算模型计算得到破裂压力。
[0016]进一步的技术方案是，所述一维地质力学模型包括：
[0017]Biot系数计算公式：
[0018][0019]式中：α为Biot系数，无因次；Δt
cma
为岩石骨架的纵波时差，μs/ft；Δt
sma
为岩石骨架的横波时差，μs/ft；ρ
ma
为岩石骨架密度，g/cm3；ρ
b
为岩石密度，g/cm3；Δt
s
为横波时差，μs/ft；Δt
c
为纵波时差，μs/ft；
[0020]抗拉强度计算公式：
[0021][0022]式中：S
t
为岩石抗拉强度，MPa；S
c
为岩石单轴强度，MPa；K
tc
为比例系数；
[0023]孔隙压力计算公式：
[0024][0025]式中：p
p
为地层孔隙压力，MPa；p
w
为静水液柱压力，MPa；e为伊顿指数，无因次；σ
v
为垂向地应力，MPa；Δt
c
为纵波时差，μs/ft；Δt
nc
为正常趋势线上对应的纵波时差，μs/ft；
[0026]水平地应力计算公式：
[0027][0028]式中：σ
H
为最大水平地应力，MPa；σ
h
为最小水平地应力，MPa；ω1为最大水平地应力方向构造应力系数，无因次；ω2为最小水平地应力方向构造应力系数，无因次。
[0029]进一步的技术方案是，所述直井破裂压力计算模型为：
[0030]p
f
＝3σ
h
‑
σ
H
‑
αp
p
+S
t
[0031]式中：p
f
为破裂压力，MPa；σ
H
为最大水平地应力，MPa；σ
h
为最小水平地应力，MPa；α为Biot系数，无因次；S
t
为岩石抗拉强度，MPa；p
p
为地层孔隙压力，MPa。
[0032]进一步的技术方案是，所述步骤S2中对测井数据进行归一化处理，其归一化公式如下：
[0033][0034]式中：X
*
为归一化的参数值；X
max
和X
min
分别为参数的最大和最小值。
[0035]进一步的技术方案是，所述基于理论引导的神经网络模型的具体预测过程包括：
[0036]步骤S31、将测井参数输入到神经网络模型当中，通过数据驱动模块预测得到地质力学参数和破裂压力数据；
[0037]步骤S32、将基于数据驱动模块预测得到的破裂压力数据与实际的破裂压力数据做均方误差，构成数据约束项损失函数；
[0038]步骤S33、将基于数据驱动模块预测得到的地质力学参数代入理论引导模块中，根据直井破裂压力计算公式得到破裂压力的计算值，并与破裂压力实际值做均方误差构成理论本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.基于理论引导神经网络的破裂压力预测方法，其特征在于，包括：步骤S1、收集已钻井的测井数据，并进行测井解释得到地质力学参数和破裂压力数据，构建神经网络模型训练和测试所需的数据集；步骤S2、对数据集的数据进行预处理，并将数据集划分为训练集和测试集；步骤S3、构建基于理论引导的神经网络模型；所述基于理论引导的神经网络模型包含数据驱动模块和理论引导模块，所述数据驱动模块包括输入层、隐藏层、输出层；步骤S4、利用训练集的数据对基于理论引导的神经网络模型进行训练，并对神经网络模型进行优化，获得优化后的神经网络模型；步骤S5、根据测试集的数据对优化后的神经网络模型预测效果进行测试，并通过评价指标进行模型预测结果评估。2.根据权利要求1所述的基于理论引导神经网络的破裂压力预测方法，其特征在于，所述测井数据包括井深、横波时差、纵波时差、补偿中子、岩性密度、井径、自然伽马；所述地质力学参数包括最小水平地应力、最大水平地应力、孔隙压力、Biot系数和抗拉强度。3.根据权利要求2所述的基于理论引导神经网络的破裂压力预测方法，其特征在于，所述步骤S1中通过一维地质力学模型计算得到地质力学参数，通过直井破裂压力计算模型计算得到破裂压力。4.根据权利要求3所述的基于理论引导神经网络的破裂压力预测方法，其特征在于，所述一维地质力学模型包括：Biot系数计算公式：式中：α为Biot系数，无因次；Δt
cma
为岩石骨架的纵波时差，μs/ft；Δt
sma
为岩石骨架的横波时差，μs/ft；ρ
ma
为岩石骨架密度，g/cm3；ρ
b
为岩石密度，g/cm3；Δt
s
为横波时差，μs/ft；Δt
c
为纵波时差，μs/ft；抗拉强度计算公式：式中：S
t
为岩石抗拉强度，MPa；S
c
为岩石单轴强度，MPa；K
tc
为比例系数；孔隙压力计算公式：式中：p
p
为地层孔隙压力，MPa；p
w
为静水液柱压力，MPa；e为伊顿指数，无因次；σ
v
为垂向地应力，MPa；Δt
c
为纵波时差，μs/ft；Δt
nc
为正常趋势线上对应的纵波时差，μs/ft；水平地应力计算公式：
式中：σ
H
为最大水平地应力，MPa；σ
h
为最小水平地应力，MPa；ω1为最大水平地应力方向构造应力系数，无因次；ω2为最小水平地应力方向构造应力系数，无因次。5.根据权利要求3所述的基于理论引导神经网络的破裂压力预测方法，其特征在于，所述直井破裂压力计算模型为：p
f
＝3σ
h
‑
σ
H
‑
αp
p
+S
t
式中：p
f
为破裂压力，MPa；σ
H
为最大水平地应力，MPa；σ
h
为最小水平地应力，MPa；α为Biot系数，无因次；S
t
为岩石抗拉强度，MPa；p
p
为地层孔隙压力，MPa。6.根据权利要求1所述的基于理论引导神经网络的破裂压力预测方法，其特征在于，所述步骤S2中对测井数据进行归一化处理，其归一化公式如下：式中：X
*
为归...

【专利技术属性】
技术研发人员：马天寿，张东洋，刘阳，陈春宇，万夫磊，徐文，谢祥锋，
申请(专利权)人：西南石油大学，
类型：发明
国别省市：