CN114662390A - 一种基于svr算法的钻井漏失压力预测方法 - Google Patents

Abstract

本发明属于钻井漏失预测领域和人工智能机器学习领域，具体涉及一种基于SVR算法的钻井漏失压力预测方法。所述方法为：从钻井井史数据库中采集数据进行预处理，创建样本集与测试集，训练SVR模型，引入松弛因子，引入拉格朗日函数进行无约束化，对偶化求导，使用SMO算法选择最优化参数，构建回归超平面得到预测函数。将预处理后得到的数据集输入训练好的SVR模型，能够得到钻井漏失压力的预测值。该方法解决了目前漏失压力数据对于堵漏防漏有重大意义但是无法直接获取的困难，而且具有结果可靠，成本较低，便于操作，可推广性强的优势，在防漏堵漏领域具有一定的指导意义和参考价值。

Description

一种基于SVR算法的钻井漏失压力预测方法

技术领域

本发明属于钻井漏失预测领域和人工智能机器学习领域，具体涉及一种基于SVR算法的钻井漏失压力预测方法。

背景技术

井漏在国内外各大油气产区普遍存在的难题，由此引起的井下复杂事故，对钻完井工程造成极大危害，产生严重的经济损失，一直是石油工程界备受关注的问题。目前，国内外防漏堵漏技术以随钻、桥堵、水泥固化堵漏为主，科学判断漏失压力、漏失通道、漏速、漏失层位等漏失性质是提升防漏堵漏一次成功率的必要前提之一。目前，井漏判断及防漏堵漏决策对个人经验依赖性强，技术可复制性差，尚无法有效满足现场井漏治理需求。漏失压力作为是井漏的三大主要特征参数之一，准确预测判断漏失压力，将有力支撑防漏堵漏技术科学筛选及施工。但是，由于漏失压力的井下数据获取难度非常大，还无法实现基于人工智能技术直接预测漏失压力，值得注意的是，与漏失压力有强相关性的地层压力、破裂压力及当量循环密度的相关预测研究已经逐渐开展，对降低漏失风险有明显的进步意义，为间接预测漏失压力提供了可能。

因此，结合循环当量密度和地层压力影响因子岩性密度、孔隙度、声波时差、钻压、转速、钻速、钻井液密度等7个参数以及破裂压力影响因子孔隙压力、钻压、钻速、转速、扭矩、钻井液密度等6个参数，最终整合形成利用循环当量密度、岩性密度、孔隙度、声波时差、钻压、转速、钻速、钻井液密度、孔隙压力、扭矩10个参数，通过及时掌握钻井施工时的这些数据变化情况，采取合理的预测方法对漏失压力进行预测。

基于人工智能机器学习领域对井漏进行预测的一种方法是采用SVR算法，传统的回归方法对于数据噪音要求高，泛化性能不佳，并且很难拟合复杂的概率分布模型。因此，提升回归模型的鲁棒性和泛化性能就是当前研究中的挑战。

SVR模型是将SVM算法与传统的回归算法相结合形成的新的更加鲁棒的回归模型，它具备以下四点优势，一是在于解决高维特征的分类问题和回归问题很有效,在特征维度大于样本数时依然有很好的效果；二是仅仅使用一部分支持向量来做超平面的决策，无需依赖全部数据；三是有大量的核函数可以使用，从而可以很灵活的来解决各种非线性的分类回归问题；四是样本量不是海量数据的时候，分类准确率高，泛化能力强。通过该方法对处理后的测井数据进行钻井漏失压力预测具有良好的应用前景。

发明内容

本发明的目的在于提供一种使用SVR算法预测漏失压力参数的方法，通过采集相同区块下不同井队的钻井井史数据，进行预处理后应用SVR模型得到漏失压力参数，解决了漏失压力的井下数据难以获取的难题，有力支撑防漏堵漏技术科学筛选及施工。为了实现上述目的，基于SVR算法的钻井漏失压力预测方法具体实施步骤如下：

1)采集钻井井史数据

收集样本数据，通过对目标区块钻井井史数据的收集、整理和验证，确定了10个输入参数的参数集：循环当量密度、岩性密度、孔隙度、声波时差、钻压、转速、钻速、钻井液密度、孔隙压力、扭矩。

2)将采集到的井史数据进行数据预处理，数据预处理内容包括数据清洗、数据集成、数据标准化、数据归约，最终得到经过预处理后的井史数据资料

上述数据预处理的进一步特征在于：

2-1)所述数据清洗的内容包括缺失值处理、噪声数据处理、数据值替换、数据转换及异常值检测；

2-2)所述数据集成的内容包括实体识别、数据冗余处理、元组重复处理及数据值冲突检测与处理；

2-3)所述数据标准化的内容包括离差标准化数据、标准差标准化数据；

2-4)所述数据归约的内容包括维归约、数量归约和数据压缩；

3)将预处理后的钻井井史数据资料作为输入，钻井漏失压力作为输出，以静液柱压力与循环压耗之和作为有效压力标准值，构建SVR漏失压力预测模型；

上述技术方案的进一步特征在于，所述的SVR漏失压力预测模型构建方案如下：

3-1)定义原始数据集为S＝A₁,A₂,A₃,…,A_n，其中A_i表示区块内一个井的井史数据：

A_i＝{(x_i,1,x_i,2,x_i,3,x_i,4,x_i,5,x_i,6,x_i,7,x_i,8,x_i,9,x_i,10),i≥10}；

式中，x₁表示循环当量密度，x₂表示岩性密度，x₃表示孔隙度，x₄表示声波时差，x₅表示钻压，x₆表示转速，x₇表示钻速，x₈表示钻井液密度，x₉表示孔隙压力，x₁₀表示扭矩。

3-2)确定SVR模型；

y＝<w，φ(x)>+b

式中，yER是模型输出，x∈Rn是模型输入，w∈R^N是表示超平面的法向量，φ(·):Rⁿ→R^N是将x从n维映射到N维的函数，b∈R是相对于过原点超平面的偏移。

3-3)确定SVR模型的原始优化问题引入松弛因子，利用拉格朗日函数将目标函数变为无约束形式，对偶化后再求导，使用SMO算法得到最优参数，构建回归超平面得到预测函数：

y＝sgn(w*^Tx+b*)

4)采用验证集对训练结果进行验证

5)正钻井漏失压力预测

根据获取的钻井井史数据，将相关参数输入该SVR漏失压力预测模型，获取相应的钻井漏失压力。

附图说明

为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

附图1是井史数据预处理框图。

附图2是SVR结构示意图。

附图3是钻井漏失压力预测方法实施流程图。

具体实施方式

为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。

图1为本发明技术方案中的井史数据预处理框图，数据预处理内容包括数据清洗、数据集成、数据标准化、数据归约，其中数据清洗包括缺失值处理、噪声数据处理、数据值替换、利用函数或映射进行数据转换及异常值检测，从而提高数据质量；数据集成包括实体识别、数据冗余处理、相关分析、元组重复处理及数据值冲突检测与处理，从而将多个数据源中的数据合并，存放于一个一致的数据存储中；数据标准化包括离差标准化数据、标准差标准化数据，从而消除特征之间量纲和取值范围的差异可能造成的影响；数据归约包括维归约、数量归约和数据压缩，从而使得尽可能保持数据完整的基础上得到数据的归约表示。

图2是本发明技术方案中的SVR结构示意图，本实例为了全面考虑主要参数对钻井漏失压力预测的影响，将预处理后的井史数据资料作为输入，漏失压力作为输出，以有效压力为标准值，构建钻井漏失压力预测SVR模型。

图3是钻井漏失压力预测方法实施流程图，具体如下：

1)采集钻井井史数据资料；

通过对钻井漏失压力所需井史数据资料整理和验证，确定了预测钻井漏失压力所需井史数据资料应包含的10个参数集：循环当量密度、岩性密度、孔隙度、声波时差、钻压、转速、钻速、钻井液密度、孔隙压力、扭矩，将所获取到的井史数据资料作为训练模型的标准样本。

2)将采集到的井史数据进行数据预处理，数据预处理内容包括数据清洗、数据集成、数据标准化、数据归约，最终得到经过预处理后的井史数据资料；

上述数据预处理的进一步特征在于：

2-1)所述数据清洗的内容包括缺失值处理、噪声数据处理、数据值替换、数据转换及异常值检测；

针对上述数据清洗的进一步特征在于：

2-1-1)使用回归方法基于推理工具填补最有可能的缺失值；

2-1-2)通过聚类方法检测离群点，将类似的值组织成簇，直观地将落在簇集合之外的值视为离群点；

2-1-3)在Pandas利用replace()方法将查询到的数据进行值的替换；

2-1-4)利用box-cox方法进行数据变换，在一定程度上减小不可观测的误差；

2-1-5)利用散点图方法绘制图像，可以较为直观地观测到是否存在异常值；

2-2)所述数据集成的内容包括实体识别、数据冗余处理、元组重复处理及数据值冲突检测与处理；

针对上述数据集成的进一步特征在于：

2-2-1)注意多来源数据的结构，确保源系统中函数依赖和参数约束与目标系统匹配；

2-2-2)使用卡方检验处理标称属性，使用相关系数处理数值属性；

2-2-3)检测重复的元组，如给定唯一的数据实体，检测是否存在两个或多个相同元组，如果存在则只保留一个；

2-2-4)对于不同数据来源的数据需要进行合并，以此确保信息传递不存在障碍；

2-3)所述数据标准化的内容包括离差标准化数据、标准差标准化数据；

针对上述数据标准化的进一步特征在于：

2-3-1)离差标准化是对原始数据所做的一种线性变换，将原始数据的数值映射到区间，统匹配，如下式：

上式进一步特征在于，其中max为样本最大值，min为样本最小值，x为样本原始值，y为样本映射值；

2-3-2)标准差标准化是对原始数据处理均值为0，标准差为1，如下式：

上式进一步特征在于，其中mean为平均值，std为标准值，x为样本原始值，y为标准差化后的样本映射值；

2-4)所述数据归约的内容包括维归约、数量归约和数据压缩；

针对上述数据归约的进一步特征在于：

2-4-1)对于上述的维归约，采用属性子集选择的方法，将不相关、弱相关或冗余的属性(或维)检测出来以后进行删除；

2-4-2)对于上述的数量归约，采用数据立方体聚类的方式，将低粒度的数据分析聚合成汇总粒度的数据分析；

2-4-3)对于上述的数据压缩，基于傅里叶变换，采用小波变换的方式，对于数据向量进行小波变换将会得到具有相同长度的小波系数，由于数据的主要能量集中在低频区域，因此仅存储一小部分最强的小波系数，就能保留近似的压缩数据。

3)将预处理后的钻井井史数据资料作为输入，钻井漏失压力作为输出，以静液柱压力与循环压耗之和作为有效压力标准值，构建SVR漏失压力预测模型；

上述技术方案的进一步特征在于，所述的SVR漏失压力预测模型构建方案如下：

3-1)定义原始数据集为S＝A₁,A₂,A₃,…,A_n，其中A_i表示区块内一个井的井史数据：

A_i＝{(x_i,1,x_i,2,x_i,3,x_i,4,x_i,5,x_i,6,x_i,7,x_i,8,x_i,9,x_i,10),i≥10}；

式中，x₁表示循环当量密度，x₂表示岩性密度，x₃表示孔隙度，x₄表示声波时差，x₅表示钻压，x₆表示转速，x₇表示钻速，x₈表示钻井液密度，x₉表示孔隙压力，x₁₀表示扭矩。

3-2)确定SVR模型；

y＝<w，φ(x)>+b

式中，y∈R是模型输出，x∈Rⁿ是模型输入，w∈R^N是表示超平面的法向量，φ(·)：Rⁿ→R^N是将x从n维映射到N维的函数，b∈R是相对于过原点超平面的偏移。

3-3)确定SVR模型的原始优化问题引入松弛因子，利用拉格朗日函数将目标函数变为无约束形式，对偶化后再求导，使用SMO算法得到最优参数，构建回归超平面得到预测函数：

y＝sgn(w*^Tx+b*)

针对上述原始优化问题的进一步特征在于：

3-3-1)SVR模型的原始优化问题：

s.t.|y_i-(w^Tx_i+b)|≤ε，i＝1，2，…，N

式中，ε为大于0的任意常数。

3-3-2)引入松弛因子的形式是：

式中，

是上界松弛因子，

是下界松弛因子，C是惩罚参数。

3-3-3)引入拉格朗日函数的形式如下：

式中，

都是拉格朗日系数。

3-3-4)优化问题满足KKT条件可以等价于对偶问题，形式如下：

3-3-5)求导后通过SMO算法求得最优参数，计算方式如下：

3-3-6)寻找一个满足

的样本点(x_k，y_k)，计算方式如下：

b^*＝y_k+∈-w*^Tx_k

b*＝y_k-∈-w*^Tx_k

为了保证求出的b*的更加具备鲁棒性，选取上述两式求取多个下界和上界的支持向量解b后再取平均。

3-3-7)构建回归超平面得到预测函数：

y＝sgn(w*^Tx+b*)

4)采用验证集对训练结果进行验证

5)正钻井漏失压力预测

根据获取的钻井井史数据，将相关参数输入该SVR漏失压力预测模型，获取相应的钻井漏失压力。

以上所述具体技术方案仅用以说明本发明，并不用于限制本发明；尽管参照上述具体技术方案对本发明进行了详细的说明，本领域的普通技术人员应当理解，凡是基于本发明的精神和原则之内，所做的任何修改和等同替换均应包含在本发明的保护范围之内。

Claims (5)

1.一种基于SVR算法的钻井漏失压力预测方法，其特征在于，包括如下步骤：

S1：收集样本数据，从钻井井史数据库中采集相关数据并进行预处理，创建样本集与测试集；

S2：从样本集数据中确定特征属性，确定特征属性的取值为SVR模型的输入，预测结果为输出；

S3：使用S2训练集求解SVR模型的原始优化问题，加入松弛因子，引入拉格朗日函数，对偶化后求导，使用SMO算法得到最优参数，最后构建回归超平面得到预测函数；

S4：利用验证集对训练结果进行验证；

S5：SVR模型的实际运用，获取钻井区块相关参数的实际钻井数据输入所述SVR模型，获取相应的钻井漏失压力。

2.根据权利要求1所述的一种基于SVR算法的钻井漏失压力预测方法，其特征在于，所述S1中预处理包括：数据清洗，数据集成，数据标准化，数据归约。

3.根据权利要求1所述的一种基于SVR算法的钻井漏失压力预测方法，其特征在于，所述S1中创建样本集与测试集包括，将预处理后的钻井参数数据进行随机分层抽样，按照预设比例10：1形成样本集与测试集。

4.根据权利要求1所述的一种基于SVR算法的钻井漏失压力预测方法，其特征在于，所述S2具体为：

确定钻井漏失参数数据为特征属性，包括循环当量密度，岩性密度，孔隙度，声波时差，钻压，转速，钻速，钻井液密度，孔隙压力，扭矩，确定以特征属性的取值为特征输入，确定以相应特征属性取值计算的钻井漏失压力为特征输出，根据确定的特征属性，依次将特征取值记为x_i，循环当量密度记为x₁,岩性密度记为x₂，孔隙度记为x₃，声波时差记为x₄,钻压记为x₅,转速记为x₆,钻速记为x₇,钻井液密度记为x₈，孔隙压力记为x₉,扭矩记为x₁₀,确定的特征属性取值记为向量X：X＝(x₁,x₂,x₃,…x₉,x₁₀)，用Y表征输出值。

5.根据权利要求1所述的一种基于SVR算法的钻井漏失压力预测方法，其特征在于，所述S3具体为：

S31、SVR模型的形式如下：

y＝<w，φ(x)>+b

式中，y∈R是模型输出，x∈Rⁿ是模型输入，w∈R^N是表示超平面的法向量，φ(·)：Rⁿ→R^N是将x从n维映射到N维的函数，b∈R是相对于过原点超平面的偏移；

S32、SVR模型的原始优化问题形式如下：

s.t.|y_i-(w^Tx_i+b)|≤ε，i＝1，2，…，N

式中，在SVR模型基础上新增的ε为大于0的任意常数；

S33、对于SVR模型选择一个惩罚参数，引入松弛因子的形式是：

式中，

是上界松弛因子，

是下界松弛因子，C是惩罚参数；

S34、对于SVR模型引入拉格朗日函数的形式如下：

式中，

都是拉格朗日系数；

S35、对于SVR模型优化问题满足KKT条件可以等价于对偶问题，形式如下：

S36、对于SVR模型对偶问题求导后通过SMO算法求得最优参数，计算方式如下：

S37、寻找一个满足

的样本点(x_k，y_k)，计算方式如下：

为了保证求出的b*的更加具备鲁棒性，选取上述两式求取多个下界和上界的支持向量解b后再取平均。

S38、构建回归超平面将前述计算所得代入即可得到预测函数：

得到的最优化后的预测函数后即可输入相关参数得到输出的漏失压力值。

Images