一种钻进状态变量的预测方法、设备及存储设备与流程

本发明涉及地质勘探领域，尤其涉及一种钻进状态变量的预测方法、设备及存储设备。

背景技术：

我国主要矿产资源与能源大量依靠国外进口，国家安全与经济发展受制于人，而保障国家资源能源安全必需立足于国内。而我国矿产资源能源储量丰富，对其进行全面开采可以缓解国内供应的相对不足。

然而在发展地质勘探开发的过程中，由于深部钻探遇到的地层层位较多，压力体系复杂，岩石类型复杂多变，具有高地应力、高地温、高渗透压及开采扰动的复杂地层力学环境，导致钻进过程中常用的传感器设备出现设备无法正常使用、信息传输困难、钻进过程干扰强烈、各系统间的耦合严重等问题。这些问题导致整个钻进系统成为了一个黑箱系统，操作人员在很大程度上只能根据经验去判断当前的钻进状态。然而由于个人经验不同，对当前的钻进状态的判断也会不同，这往往会给钻进过程带来截然不同的控制效果。根据经验性，使得对钻进状态的判断和控制具有很强的个人主观性。因此，亟需研究一种预测钻进状态的方法，对钻进状态的判断和控制具有重要的意义。

技术实现要素：

为了解决上述问题，本发明提供了一种钻进状态变量的预测方法、设备及存储设备，一种钻进状态变量的预测方法，主要包括以下步骤：

s101：根据钻进机理，确定出与钻进状态变量相关的操作参数，得到所述操作参数的数据；

s102：利用快速傅里叶变化的方法，对所述数据进行低通滤波，滤除所述数据中存在的扰动；

s103：根据滤波处理过的数据与所述操作参数对应的钻进状态变量，采用支持向量回归算法建立钻进状态变量的预测模型；

s104：利用模型更新方法对所述模型对进行更新，得到最终的钻进状态变量的预测模型；

s105：将当前时刻的操作参数数据和钻进状态变量数据输入到最终的钻进状态变量的预测模型，预测出下一时刻的钻进状态变量，得到下一时刻的钻进状态。

进一步地，在步骤s101中，操作参数包括：钻压、转速、泵量与钻井液密度。

进一步地，在步骤s102中，利用所述快速傅里叶滤波方法进行数据滤波的过程为：

(1)预设截止频率的度量参数fc，根据公式(1)和截止频率的度量参数fc，求解截止频率fc：

其中，fc表示截止频率，fs表示采样频率，fc为截止频率的度量参数；

(2)根据截止频率fc，滤除所述数据中存在的大于截止频率fc的高频扰动，保留需要的信号频率。

进一步地，在步骤s103中，所述钻进状态变量的预测模型是以当前时刻的钻进操作参数数据与钻进状态变量数据作为模型输入，以下一时刻的钻进状态变量数据作为模型输出进行训练得到；所述钻进状态变量包括当前钻速与总池体积。

一种存储设备，所述存储设备存储指令及数据用于实现一种钻进状态变量的预测方法。

一种钻进状态变量的预测设备，包括：处理器及所述存储设备；所述处理器加载并执行所述存储设备中的指令及数据用于实现一种钻进状态变量的预测方法。

本发明提供的技术方案带来的有益效果是：能够实现钻进状态变量的有效预测，为钻进操作提供了良好的指导意见。

附图说明

下面将结合附图及实施例对本发明作进一步说明，附图中：

图1是本发明实施例中一种钻进状态变量的预测方法的流程图；

图2是本发明实施例中钻进状态变量预测方法结构图；

图3是本发明实施例中模型更新方法的示意图；

图4是本发明实施例中硬件设备工作的示意图。

具体实施方式

为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图详细说明本发明的具体实施方式。

本发明的实施例提供了一种钻进状态变量的预测方法、设备及存储设备。

请参考图1，图1是本发明实施例中一种钻进状态变量的预测方法的流程图，具体包括以下步骤：

s101：根据钻进机理，确定出与钻进状态变量相关的操作参数，得到所述操作参数的数据；钻进过程中，钻速与泥浆池总池体积是对钻进效率与安全性衡量的重要指标；钻速的调整由钻压与转速这两大操作量来决定，泥浆池总池体积则由泵量与泥浆密度来决定；为了实现对钻速与泥浆池总池体积的准确预测，得到与钻进状态变量相关的操作参数，操作参数包括：钻压、转速、泵量与钻井液密度；

s102：利用快速傅里叶变化的方法，对所述数据进行低通滤波，滤除所述数据中存在的扰动；利用所述快速傅里叶滤波方法进行数据滤波的过程为：

(1)预设截止频率的度量参数fc，根据公式(1)和截止频率的度量参数fc，求解截止频率fc：

其中，fc表示截止频率，fs表示采样频率，fc为截止频率的度量参数。

(2)根据截止频率fc，滤除所述数据中存在的大于截止频率fc的高频扰动，保留需要的信号频率。

s103：根据滤波处理过的数据与所述操作参数对应的钻进状态变量数据，采用支持向量回归算法建立钻进状态变量的预测模型；所述钻进状态变量的预测模型是以当前时刻的钻进操作参数数据与钻进状态变量数据作为模型输入，以下一时刻的钻进状态变量数据作为模型输出进行训练得到；所述钻进状态变量包括当前钻速与总池体积；

s104：利用模型更新方法对所述模型对进行在线更新，得到最终的钻进状态变量的预测模型；利用模型更新方法对所述模型对进行更新的过程为：将滤波处理过的数据与所述操作参数对应的钻进状态变量分为测试集与训练集，训练集与测试集分别具有n组数据与k组数据，每组数据中包括当前时刻的钻进操作参数数据、钻进状态变量数据和下一时刻的实际钻进状态变量数据；每当预测模型预测完测试集的一组数据时，记录得到的预测结果，并将该组数据转移到训练集中；根据训练集中的所有数据对预测模型进行更新，并将更新后的预测模型用于测试集中下一组数据的处理；最终训练集中包含n+k-1组数据；通过将测试集的数据不断向训练集中转移，确保训练集中不断包含新的数据信息，实现对预测模型的不断更新，得到最终的钻进状态变量的预测模型；以便得到精度更高的预测结果；所以，对模型进行不断的更新，提高了模型的预测精度；

s105：将当前时刻的操作参数数据和钻进状态变量数据输入到最终的钻进状态变量的预测模型，预测出下一时刻的钻进状态变量数据，得到下一时刻的钻进状态。

请参考图2，图2是本发明实施例中钻进状态变量预测方法结构图，利用模型进行预测的具体步骤为：

(1)获取操作参数数据

收集钻进现场2805米至2833米井段实时获取的数据，剔除其中停钻与起下钻的数据，共1100组数据，选取钻速与泥浆池总池体积作为状态变量，选取钻压、转速、泵量与钻井液密度作为操作参数。

(2)基于快速傅里叶变换的数据滤波

所获数据中存在一些高频扰动，因此利用快速傅里叶变化的方法对数据进行滤波，滤除数据中存在的高频扰动。

(3)建立钻进状态变量预测模型

不同变量的幅值也是不同的，幅值的不同导致了样本数据之间的差异很大，为了消除幅值的影响，根据公式(2)对每一样本数据进行归一化处理：

其中，x为样本数据的值，xnorm为经过归一化处理的输入值，在本发明实施例中，归一化后数据的变化范围为[-1,1]。

将当前操作参数与钻进状态变量作为预测模型的输入，将下一时刻的钻进状态变量作为模型输出。从1100组样本中选择1080组来训练模型，用剩余20组来测试模型性能。利用支持向量回归(svr)方法构建模型，然后利用模型更新方法对模型进行不断更新。

请参考图3，图3是本发明实施例中模型更新方法的示意图，首先将数据按时间顺序分为测试集与训练集，假设测试集与训练集分别具有n组数据与k组数据；然后每当所述模型预测完测试集的一组数据时，就将所述数据转移到训练集中，并根据训练集中的数据对所述模型进行更新，记录得到的预测结果；通过预测集的数据不断向训练集中转移，最终训练集中将包含n+k-1组数据；训练集中不断包含新的数据信息，实现对所述模型进行不断的更新，得到精度更高的预测结果；所以，对模型进行不断的更新，提高了模型的预测精度。

采用均方误差的方法对模型的预测性能进行衡量，根据公式(3)得到真实值与预测值间的均方误差：

其中，mse为均方误差，yi表示真实值，表示预测值；mse越小，表示模型的预测性能越好，预测值越接近真实值，模型的预测精度越高。

在本发明实施例中，预测得到钻速与总池体积的mse分别为0.000465和0.00110，这说明本发明提供的技术方案具有较高的预测精度。

请参见图4，图4是本发明实施例的硬件设备工作示意图，所述硬件设备具体包括：一种钻进状态变量的预测设备401、处理器402及存储设备403。

一种钻进状态变量的预测设备401：所述一种钻进状态变量的预测设备401实现所述一种钻进状态变量的预测方法。

处理器402：所述处理器402加载并执行所述存储设备403中的指令及数据用于实现所述一种钻进状态变量的预测方法。

存储设备403：所述存储设备403存储指令及数据；所述存储设备403用于实现所述一种钻进状态变量的预测方法。

本发明的有益效果是：能够实现钻进状态变量的有效预测，为钻进操作提供了良好的指导意见。

以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。