一种基于岩屑运移量的PDC钻头钻进参数的优化方法

一种基于岩屑运移量的pdc钻头钻进参数的优化方法
技术领域
1.本发明属于石油开采技术领域，具体涉及一种基于岩屑运移量的pdc钻头 钻进参数的优化方法。


背景技术：

2.在石油或者天然气开采过程中，pdc钻头是常用的钻头之一，在使用过程 中，首先通过给钻头施加钻头来达到破碎岩石的目的，然后再给钻头施加转速， 用刀翼进行刮削，进而达到钻进的目的，同时通过钻井液来运移岩屑和冷却钻头。 但是由于钻头所处的工作环境非常复杂，经常处于几千米的地下，人们对于钻头 的具体工况难以预知，对于钻头所采用的钻进参数难以评估是否正确。
3.pdc钻头所采用的钻进参数对钻头的使用时长以及钻进效率有着非常大地 影响，并且各个参数之间互相影响，很难实现每个参数达到最优，因此只能收集 足够多的历史数据，从历史钻井数据中寻找规律，针对不同的岩层以及深度，确 定最优指标，对进行钻进参数选择。


技术实现要素：

4.本发明要解决的技术问题是：提供一种基于岩屑运移量的pdc钻头钻进参 数的优化方法，用于确定钻头钻进参数的最优指标。
5.本发明为解决上述技术问题所采取的技术方案为：一种基于岩屑运移量的 pdc钻头钻进参数的优化方法，包括以下步骤：
6.s0：获取历史数据；
7.s1：根据钻井实时工况，确定钻进过程的特征、标签和参数的取值，根据算 法确定重要度靠前的参数；
8.s2：根据最大互信息值法计算各个特征与标签的互信息值并排序，从而确定 各个特征与标签之间的重要度，根据重要度选择从已有特征中选择若干个重要度 大的特征以降低数据集维度、优化系统的特定指标；
9.s3：根据步骤s2选出的特征数量及取值范围设计对应的正交试验表；利用 正交试验法分析已选择的特征，通过正交试验表进行试验找出各个特征的最优搭 配；
10.s4：结合正交试验表的结果和问题采用包括观察分析、极差分析、方差分析 的分析方法分析试验结果，确定最佳参数组合。
11.按上述方案，所述的步骤s1中，具体步骤为：
12.s11：选择包括转速、岩石特性、井深、钻井液密度、钻井液喷出速度、钻 压、井底温度的参数作为特征；
13.s12：以单位时间内向井上运移的岩屑量作为标签；
14.s13：获取在不同的转速、岩石特性、井深、钻井液密度、钻井液喷出速度、 钻压、井底温度的情况下，单位时间内向井上运移的岩屑量。
15.按上述方案，所述的步骤s2中，具体步骤为：设a、b分别是在x、y方 向上的划分格子的个数，形成网格分布；设b是变量，则最大互信息值的计算 公式为：
[0016][0017]
进一步的，所述的步骤s2中，具体步骤为：b的大小设置为数据量的0.6 次方。
[0018]
按上述方案，所述的步骤s3中，具体步骤为：在正交试验表中，称所选的 特征为因素，因素是影响结果的原因；所获得特征的取值范围为水平，水平是因 素本身的不同条件和状态，是因素的取值范围或者大小。
[0019]
进一步的，所述的步骤s3中，设计正交试验表的具体步骤为：
[0020]
s31：确定优化指标为单位时间内岩屑的运移量；
[0021]
s32：确定本次正交试验的因素与水平；设l为正交表符号；x表示行数即 试验的次数；m为因素水平数，范围为1～m；n为列数，表示影响结果的因素 个数；按照如下公式设计正交表：
[0022]
l
x
(mn)；
[0023]
正交试验表的每列都由相同的数字组成；任意两列横向形成的数字对出现的 次数相等；
[0024]
s33：根据因素和水平选择正交试验表，若所选择的列数比因素的数量多时， 则让多余的列为空；
[0025]
s34：列出试验条件，正交试验表的每行表示一次试验选择的方案。
[0026]
按上述方案，所述的步骤s4中，具体步骤为：
[0027]
s41：对试验结果进行观察分析，根据每一组试验的结果，直接找出最好的 试验条件；
[0028]
s42：对实验结果进行极差分析，具体步骤为：
[0029]
s421：计算各因素的相同水平试验结果之和；
[0030]
s422：计算各因素的相同水平试验结果的平均值，比较各水平试验结果之和 的平均值的大小，值越大越好；以每个因素的实际取值顺序作为横坐标，以对应 于该水平所做试验的结果的平均值为纵坐标作趋势图；
[0031]
s423：计算各列中相同水平试验结果之和的平均值的极差值；比较极差值的 大小，极差值大的因素表示在相同水平下该因素对结果的影响大，则该因素为主 要因素；
[0032]
s43：利用正交试验表对试验结果进行方差分析，将总偏差平方和分解为各 个因素效应的偏差平方和与误差效应平方和，并固定在正交表的每一列上，计算 因素的平均变动和误差的平均变动，然后利用f检验法对因素进行显著性检验。
[0033]
进一步的，所述的步骤s43中，具体步骤为：
[0034]
s431：设yk为第k号试验条件组的实验结果，k＝1，2，3
…
，n；为n次试验 的平均值；f为自由度；t
ji
为因素j的第i个水平对应的试验结果之和， i＝1，2，3
…
m；则计算各列的偏差平方和s为：
[0035]
[0036][0037]
se＝s
t-∑sj,fe＝f
t-∑fj；
[0038]
s432：计算方差v、方差比f和影响率ρ：
[0039]vj
＝sj/fj，
[0040]ve
＝se/fe，
[0041]fj
＝vj/ve，
[0042]
ρj＝s
′j/s
t
＝(s-f
·ve
)/s
t
；
[0043]
s433：列出包括偏差平方和s、自由度f、方差v、方差比f和影响率ρ的 方差分析表；从方差分析表中找出具有显著性的fj或较大的ρj，则对应的因素 j即为对试验质量特性有最大、次大影响的因素，j＝a，b，c，
…
；根据t
ji
值找 出影响最大的因素j的最好水平，最终得出最好的因素水平组合。
[0044]
按上述方案，所述的步骤s3中，具体步骤为：
[0045]
若观察分析结果与极差分析结果一致，则判断找出了有代表性的条件；若观 察分析结果与极差分析结果不一致，则通过极差分析得出的结果和给出的因素水 平补做试验；
[0046]
若使用不同分析法得到一致的结论时，则判断找到了最宜因素水平组合；若 使用不同分析法得到不一致的结论，且方差分析或极差分析结论不处于试验计划 内，则按此结论做验证试验；若此验证试验结果显示为较佳，则判断得到了最宜 因素水平组合；
[0047]
如果验证试验结果与观察分析结果不一致，则做第二次试验，必要时重选因 素与水平。
[0048]
一种计算机存储介质，其内存储有可被计算机处理器执行的计算机程序，该 计算机程序执行一种基于岩屑运移量的pdc钻头钻进参数的优化方法。
[0049]
本发明的有益效果为：
[0050]
1.本发明的一种基于岩屑运移量的pdc钻头钻进参数的优化方法，通过收 集历史数据，从历史钻井数据中寻找规律，针对不同的岩层以及深度确定最优指 标，实现了确定钻头钻进参数的最优指标的功能。
[0051]
2.本发明通过更好地选择钻井参数，提高了钻头的使用寿命和钻进效率。
[0052]
3.本发明通过对钻头所采用的钻进参数的正确评估，实现了对钻头的具体工 况的了解。
附图说明
[0053]
图1是本发明实施例的流程图。
具体实施方式
[0054]
下面结合附图和具体实施方式对本发明作进一步详细的说明。
[0055]
参见图1，本发明实施例的一种基于岩屑运移量的pdc钻头钻进参数的优 化方法，包括以下步骤：
[0056]
s0：获取历史数据；
[0057]
s1：确定钻进过程中各种特征、标签和各个参数的取值；
[0058]
步骤s1进一步包括：根据钻井实时工况，选择尽可能多的特征，获取或者 测量各个参数的取值，将单位时间内岩屑的运移量作为优化的目标即标签。
[0059]
根据钻井过程中的实际情况，主要选择以下参数作为特征：转速、岩石特性、 井深、钻井液密度、钻井液喷出速度、钻压、井底温度，以单位时间内井上所运 移的岩屑量作为标签。
[0060]
获取在不同的转速、岩层、井深、钻压、钻井液密度、钻井液喷出速度、井 底温度的情况下，单位时间内岩屑的运移量。
[0061]
特征选择：是指从已有的m个特征(feature)中选择n个特征使得系统的特 定指标最优化，是从原始特征中选择出一些最有效特征以降低数据集维度的过 程。
[0062]
s2：确定各个特征对标签的影响程度大小，根据具体情况确定一定数量的特 征；
[0063]
在步骤s2中，根据最大互信息值法，确定各个特征与标签即岩屑运移量之 间的重要度。最大互信息值的计算公式为：
[0064][0065]
其中a,b是在x,y方向上的划分格子的个数，本质上就是网格分布，b是变 量，b的大小设置是数据量的0.6次方左右。
[0066]
根据上述公式所算出的重要度，选择一定数量的同时重要度大的特征。
[0067]
上述的特征选择采用互信息值法来衡量两个变量之间的相关度，如果两个变 量之间存在一定的相关性，那么在这两个变量的散点图上进行某种网格划分之 后，根据这两个变量在网格中的近似概率密度分布，可计算这两个变量的互信息， 正则化后，该值可用于衡量这两个变量之间的相关性。其计算公式为：
[0068][0069]
其中a,b是在x,y方向上的划分格子的个数，本质上就是网格分布，b是变 量，b的大小设置是数据量的0.6次方左右。
[0070]
通过计算各个特征与标签的互信息值，可根据互信息值的大小进行排序，从 而得出各个特征影响岩屑运移量的重要程度。
[0071]
根据所计算的互信息值大小，关注互信息值大的特征的变化对岩屑运移量的 影响。
[0072]
s3：设计正交试验表并进行试验；
[0073]
在步骤s3中，具体包括以下步骤：
[0074]
根据步骤s2中所选出的特征数量以及取值范围设计对应的正交试验表。
[0075]
在正交试验表中，所选的特征称为因素，所获得特征的取值范围为水平，确 定本次正交试验的因素数量与水平数量。
[0076]
设计正交试验表时应具有以下特性：每列都有相同的数字组成；任意两列， 其横向形成的数字对出现的次数都是相等的，当所选择的列数比因素的数量多 时，那么就让多余的列空着。
[0077]
对于已填好的实验表来说，每一行就表示一次试验，即表示一种可选择的方 案。
[0078]
利用正交实验法对已选择的特征进行分析，制作正交表来找出各个特征的最 优
搭配。
[0079]
在制作正交表时，首先确定优化指标，然后确定因素和水平。因素指的是影 响结果的原因，水平指的是因素本身地不同的条件和状态即因素的取值范围或者 大小。
[0080]
在此方法中，优化指标指的是单位时间内岩屑的运移量，因素指的是上述经 过算法所选择出的特征，水平指的是特征的取值范围。
[0081]
正交表设计时按照以下公式设计：l
x
(mn)。
[0082]
其中l为正交表符号；x表示行数，即实验的次数；m为因素水平数，范围 为1-m，n为列数，表示影响结果的因素个数。
[0083]
设计正交表。每列都由相同的数字组成；并且任选其中两列，在其横向组成 的数字对出现的次数都是相同的。
[0084]
选择正交表。根据因素多少，水平多少来选择正交表，既达到目的，又减少 工作量。当选择的正交表的列数比行数多时，那么就让多余的列空出来。
[0085]
表1为l9(34)的样表，表2为因素水平表。
[0086]
表1
[0087][0088]
表2
[0089][0090]
列出试验条件。对于已填好的实验表来说，每一行就表示一次试验，即表示 一种可选择的方案。
[0091]
s4：选择数据分析方法分析试验结果，确定最佳参数组合。
[0092]
在步骤s4中，结合正交表结果以及问题采用不同的分析方法。一般有以下 几种方法：观察分析、极差分析、方差分析。根据试验类型以及具体要求，所采 用的方法也有所不同。
[0093]
试验结果分析。根据正交试验的结果入手，结合正交表以及问题采用不同的 分析
方法。一般有以下几种方法：观察分析、极差分析、方差分析。根据试验类 型以及具体要求，所采用的方法也有所不同。
[0094]
s41：试验结果的观察分析：根据每一组试验的结果，直接找出最好的试验 条件。
[0095]
s42：实验结果的极差分析：首先计算出各因素的相同水平试验结果之和； 其次计算其平均值；最后计算各列中相同水平试验结果之和的平均值的极差值。
[0096]
比较各水平结果之和的平均值的大小，当表示不同意义时，对其值的要求叶 有所不同，在本发明中要求其值越大越好。
[0097]
以每个因素的实际取值顺序作为横坐标，以对应于该水平所做试验的结果的 平均值为纵坐标作趋势图。
[0098]
比较极差的大小，极差大的因素表示几个水平对于结果的影响大，即为主要 因素。
[0099]
判断观察分析结果与极差分析结果，当两者结果一致时，即找出了有代表性 的条件。当两者结果不一致时，必须通过极差分析得出的结果给出的因素水平补 做试验。
[0100]
s43：试验结果的方差分析：方差分析的基础是将总偏差平方和分解为各个 因素效应的偏差平方和与误差效应平方和。
[0101]
正交试验结果的方差分析是利用正交表表来进行，将这种偏差平方和分解并 固定在正交表的每一列上。进一步计算出因素的平均变动和误差的平均变动，然 后利用f检验法对因素进行显著性检验。其步骤如下所示：
[0102]
s431：计算各列的偏差平方和s；
[0103]
设yk为第k号试验条件组的实验结果，k＝1，2，3
…
，n；为n次试验的平均 值；t
ji
为因素j的第i个水平对应的试验结果之和，i＝1，2，3
…
m。则有：
[0104][0105][0106]
se＝s
t-∑sj,fe＝f
t-∑fj[0107]
s432：计算方差v、方差比f和影响率ρ：
[0108]vj
＝sj/fj[0109]ve
＝se/fe[0110]fj
＝vj/ve[0111]
ρj＝s
′j/s
t
＝(s-f
·ve
)/s
t
[0112]
s433：表3为方差分析表：
[0113]
表3
[0114][0115]
从方差分析表中找出具有显著性的fj或较大的ρj，则它们对应的因素j(j 取a，b，c，
…
)即为对试验质量特性有最大、次大影响的因素。同时根据t
ji
值找出影响最大的因素j的最好水平，最终得出最好的因素水平组合。
[0116]
当使用不同的分析法的结论一致时，就意味着找到了最宜因素水平组合，即 最宜参数组合。
[0117]
如果不同分析法的结论不一致，且方差分析或极差分析结论不处于试验计划 之内，则应按此结论做验证试验。若此验证试验结果显示为较佳，就得到了最宜 因素水平组合。
[0118]
如果验证试验结果不如直观分析的结论，则意味着：a可能存在其他因素的 作用：b本试验较为复杂；c也可能直观分析结论确实是佳结论。这时应做第二 次试验，必要时重选因素与水平。
[0119]
以上实施例仅用于说明本发明的设计思想和特点，其目的在于使本领域内的 技术人员能够了解本发明的内容并据以实施，本发明的保护范围不限于上述实施 例。所以，凡依据本发明所揭示的原理、设计思路所作的等同变化或修饰，均在 本发明的保护范围之内。