CN113047775A - 一种井壁稳定性约束下的钻进轨迹多目标优化方法 - Google Patents
一种井壁稳定性约束下的钻进轨迹多目标优化方法 Download PDFInfo
- Publication number
- CN113047775A CN113047775A CN202110327665.3A CN202110327665A CN113047775A CN 113047775 A CN113047775 A CN 113047775A CN 202110327665 A CN202110327665 A CN 202110327665A CN 113047775 A CN113047775 A CN 113047775A
- Authority
- CN
- China
- Prior art keywords
- constraint
- track
- drilling
- wall stability
- individuals
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000005553 drilling Methods 0.000 title claims abstract description 46
- 238000005457 optimization Methods 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000013461 design Methods 0.000 claims abstract description 23
- 239000002002 slurry Substances 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000002068 genetic effect Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 2
- 230000008859 change Effects 0.000 abstract description 2
- 230000002349 favourable effect Effects 0.000 abstract 1
- 230000009286 beneficial effect Effects 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B7/00—Special methods or apparatus for drilling
- E21B7/04—Directional drilling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Mining & Mineral Resources (AREA)
- Pure & Applied Mathematics (AREA)
- Geology (AREA)
- Computational Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Earth Drilling (AREA)
Abstract
本发明提供一种井壁稳定性约束下的钻进轨迹多目标优化方法,以轨迹长度、钻柱扭矩作为优化目标,以中靶误差,安全泥浆密度上限和下限为约束,建立井壁稳定性约束下的钻进轨迹优化模型;利用结合拐点和反世代距离的多目标进化算法求解井壁稳定性约束下的轨迹优化问题,得到一组互不支配的轨迹设计方案集;选择在各点处安全泥浆密度值变化最小的方案作为最终的轨迹设计方案,本发明建立的井壁稳定性约束下的钻井轨迹优化问题模型,能够在轨迹设计中给出安全钻进方向;有效解决井壁稳定性约束下的钻井轨迹优化问题,有利于在实际钻井工程的应用。
Description
技术领域
本发明涉及复杂地质钻进过程智能控制领域,尤其涉及一种井壁稳定性约束下的钻进轨迹多目标优化方法。
背景技术
定向水平井是石油、天然气、页岩气等资源勘探开发的主要形式,如何在钻前进行轨迹设计来提高钻进效率、保证钻进安全,是钻进轨迹优化要研究的问题。
许多已有的研究通过最小化钻进轨迹长度来降低工程成本。例如用序列无约束最小化方法、数论序贯优化算法等优化造斜点和造斜率。这些基于梯度的优化方法对模型的要求较高,或者需要简化条件。对于复杂的轨迹模型,遗传算法、粒子群算法等智能优化算法能够有效解决其优化问题。然而更短的轨迹长度往往意味着更复杂的井眼结构。这将会产生快速变化的井斜角和方位角,以及较大的狗腿角,从而加大钻进安全风险。因此,有必要将钻进轨迹优化问题归为多目标优化问题进行求解。
多目标进化算法在单次计算中得到一系列非劣解集,且不必对优化问题模型进行处理,有利于保留问题信息。已有研究以最小化轨迹长度和钻柱扭矩为目标,来得到同时考虑了效率和安全性的轨迹设计方案,但是尚未在轨迹设计方案中考虑到井壁稳定性的约束。而有关井壁稳定性与钻进方向的研究,仅给出了各类地应力对应的安全钻进方向。
因此,本发明将井壁稳定性约束引入到钻进轨迹多目标优化模型中,从而使轨迹设计方案给出的钻进方向都能保证井壁稳定。
发明内容
有鉴于此,本发明提供一种井壁稳定性约束下的钻进轨迹多目标优化方法,包括以下步骤:
S1、以轨迹长度L、钻柱扭矩T作为优化目标,以中靶误差e、井深s处的安全泥浆密度上限MWhigh和下限MWlow为约束,建立井壁稳定性约束下的钻进轨迹多目标优化问题模型;
S2、利用结合拐点和反世代距离的多目标进化算法求解步骤S1中井壁稳定性约束下的轨迹优化问题,得到一组互不支配的轨迹设计方案集;
S3、计算轨迹设计方案集在各关键点处的安全泥浆密度上限和下限,选择上限和下限标准差最小的方案作为最终的轨迹设计方案。
本发明提供的技术方案带来的有益效果是:
1)本发明建立的井壁稳定性约束下的钻井轨迹优化问题模型,能够在轨迹设计中给出安全钻进方向;
2)本发明提供的优化方法能够有效解决井壁稳定性约束下的钻井轨迹优化问题,有利于在实际钻井工程的应用。
附图说明
图1是本发明一种井壁稳定性约束下的钻进轨迹多目标优化方法的流程;
图2是本发明实施例中待优化的钻进轨迹的示意图;
图3是本发明实施例中得到的帕累托解集中的解在目标空间内的分布图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地描述。
请参考图1,本发明的实施例提供了一种井壁稳定性约束下的钻进轨迹多目标优化方法,包括以下步骤:
S1、以轨迹长度L、钻柱扭矩T作为优化目标,以中靶误差e、井深s处的安全泥浆密度上限MWhigh和下限MWlow为约束,建立井壁稳定性约束下的钻进轨迹优化问题模型;
上式用于计算第i段轨迹上井深s处的井斜角和方位角,基于上式,可以得到井深s处某点的方位坐标为:
待优化的钻进轨迹如图2所示,建立的钻进轨迹多目标优化问题模型为:
min[f1(x)=L,f2(x)=T]
其中,优化目标L为轨迹长度,T为钻柱扭矩,轨迹长度为
L=LOA+LAB+LBC+LCF
直井段的钻柱扭矩T和所受拉力F分别为
Ti=μbwpLirp sinαi
Fi=Fi-1+bwpLicosαi
上式中,井壁与钻柱间的摩擦系数μ=0.66,浮力系数b=0.7,钻柱的重力系数wp=0.98kN/m,钻柱的半径rp=0.33m,αi是井斜角,弧线井段的钻柱扭矩Ti和所受拉力Fi分别为:
Ti=μrpFi/ωi
上式中,ωi是弧线段起始点的工具面角,所以轨迹的总钻柱扭矩为:
T=TOA+TAB+TBC+TCF
约束函数e是中靶误差:
其中,N(L)、E(L)是轨迹终点的北、东方向坐标,Nta=16m、Eta=175m是目标靶区中心的北、东坐标。
约束函数安全泥浆密度上限安全泥浆密度下限重力常数g=9.8m/s2,plow(s)和phigh(s)分别是最低和最高的极限压力,D(s)是井深为s时对应的井的垂直深度;MWhb和MWlb分别是安全泥浆密度上限和下限的边界值。经分析,设置MWlb=1.4,而所有深度和钻进方向对应的phigh值均大于5×103kg/m3,也就说一般情况下设置的泥浆密度都不会超过安全值,所以本发明在此实施例中不限制MWhb;
S2、结合拐点和反世代距离的多目标进化算法求解井壁稳定性约束条件下轨迹优化问题,具体步骤如下:
S21、输入参数:井壁稳定性约束下的钻进轨迹优化问题、种群大小N=200、迭代终止条件为迭代200次(根据具体情况设定迭代次数)、阈值参数e=0.5、P=0.5;
S22、初始化:随机产生N个体(以各井段长度、曲率、起始点工具面角为决策参数),初始化计算反世代距离的参考点、存档集合、拐点,对随机个体进行非支配排序;
f1 max和f1 min是第k非支配等级中的解对应的轨迹长度最大值和最小值,f2 max和f2 min是解所对应钻柱扭矩的最大值和最小值,qgen是一个调节邻域大小的参数:
P是一个控制拐点数量的阈值参数,q0=1,p0=0。
为了找到对反转世代距离贡献相对更小的解,定义
R是参考点的集合,r是参考点,X是一些同等级的非劣解集,x为同等级的非劣解,X*是X中对反转世代距离无贡献的解集,x*为对反转世代距离无贡献的解,一个较小的IGD-NS要求X中含有尽可能少的对反转世代距离无贡献的解。
S23、产生子代:根据约束e≤16、MWlow(sA)≤1.4、MWlow(sB)≤1.4、MWlow(sC)≤1.4,计算个体的违约值;用锦标赛策略选择亲代,优先选择违约值小的个体,若违约值一致则选择支配等级更高的个体,若支配等级一致则选择是拐点的个体。用基因算法中的交叉变异产生子代;
S24、非支配排序:对于可行解,根据轨迹长度和钻柱扭矩进行非支配排序,非支配等级高的个体具有更小的轨迹长度和钻柱扭矩,互不支配的个体至少有一个目标函数值不小于其它个体的目标函数值;对于不可行解,根据违约值排序,违约值小的个体具有更高的非支配等级;可行解支配不可行解;
S25、环境选择:选择前k-1个支配等级的个体,对于第k级的个体,进行精英选择策略来删减,使种群数量达到N;
若当前可行解(违约值≤0的解)的数量大于e*N,(0<e<1),则保留第k级中可能成为拐点的点,即将第k级中的点按与邻域内直线M的距离排序,按顺序删除距离更小的点,直到种群数量达到N;
若当前可行解的数量小于等于e*N,在X中随机选择解xp,将其排除后计算得到的IGD-NS值记为Ip,则最小的I所对应的解将会被认为对逆世代距离贡献最小,不断删除对逆世代距离贡献最小的解,直到种群数量达到N;
S26、返回到S23进行循环,直到迭代次数达到200次,输出当前种群。
得到的解集所对应的帕累托前沿如图3所示,从图中可以看出本发明提出的方法在井壁稳定性约束下的钻进轨迹优化问题上,得到的解集具有较好的收敛性和分布性。
S3、计算轨迹设计方案在各关键点处的MWlow,选择标准差最小的方案作为最终的轨迹设计方案,所得最终轨迹设计方案如表2所示,是轨迹设计方案在各点处MWlow的平均值,括号内是标准差;其中,所述关键点是根据需要人为选择的轨迹点。
表2安全泥浆密度变化最小的轨迹设计方案
在不冲突的情况下,本文中上述实施例及实施例中的特征可以相互结合。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种井壁稳定性约束下的钻进轨迹多目标优化方法,其特征在于,包括以下步骤:
S1、以轨迹长度L、钻柱扭矩T作为优化目标,以中靶误差e、井深s处的安全泥浆密度上限MWhigh和下限MWlow为约束,建立井壁稳定性约束下的钻进轨迹多目标优化问题模型;
S2、利用结合拐点和反世代距离的多目标进化算法求解步骤S1中井壁稳定性约束下的轨迹优化问题,得到一组互不支配的轨迹设计方案集;
S3、计算轨迹设计方案集在各关键点处的安全泥浆密度上限和下限,选择上限和下限标准差最小的方案作为最终的轨迹设计方案。
6.根据权利要求1所述的一种井壁稳定性约束下的钻进轨迹多目标优化方法,其特征在于,步骤S2中结合拐点和反世代距离的多目标进化算法为:
S21、输入参数:井壁稳定性约束下的钻进轨迹优化问题、种群大小N、迭代终止条件;
S22、初始化:随机产生N个体,初始化计算反世代距离的参考点、存档集合、拐点,对随机个体进行非支配排序;
S23、产生子代:根据中靶误差、安全泥浆密度约束计算个体的违约值,用锦标赛策略选择亲代,优先选择违约值小的个体,若违约值一致则选择支配等级更高的个体,若支配等级一致则选择是拐点的个体,用基因算法中的交叉变异产生子代;
S24、非支配排序:对于可行解,根据轨迹长度和钻柱扭矩进行非支配排序,非支配等级高的个体具有更小的轨迹长度和钻柱扭矩,互不支配的个体至少有一个目标函数值不小于其它个体的目标函数值;对于不可行解,根据违约值排序,违约值小的个体具有更高的非支配等级;可行解支配不可行解;
S25、环境选择:选择前k-1个支配等级的个体,对于第k级的个体,若当前可行解的数量大于e*N,则选择第k级中的拐点;否则,删除对逆世代距离没有贡献的解,直到种群数量达到N;
S26、若达到终止条件,则输出当前种群,否则返回到步骤S23进行循环。
7.根据权利要求1所述的一种井壁稳定性约束下的钻进轨迹多目标优化方法,其特征在于,步骤S3中给出最终的轨迹设计方案的步骤为:
S31、计算所有设计方案所对应轨迹的关键点处安全泥浆密度上限MW1 high,MW2 high,…MWn high,下限MW1 low,MW2 low,…MWn low;n为安全泥浆密度上限和下限的个数;
S32、求取各关键点处安全泥浆密度上限MWstd high和下限的标准差MWstd low;
S33、选择标准差最小的设计方案作为最终的轨迹设计方案。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110327665.3A CN113047775B (zh) | 2021-03-26 | 2021-03-26 | 一种井壁稳定性约束下的钻进轨迹多目标优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110327665.3A CN113047775B (zh) | 2021-03-26 | 2021-03-26 | 一种井壁稳定性约束下的钻进轨迹多目标优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113047775A true CN113047775A (zh) | 2021-06-29 |
CN113047775B CN113047775B (zh) | 2022-01-18 |
Family
ID=76515879
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110327665.3A Active CN113047775B (zh) | 2021-03-26 | 2021-03-26 | 一种井壁稳定性约束下的钻进轨迹多目标优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113047775B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2910759B1 (ja) * | 1998-05-19 | 1999-06-23 | 日本電気株式会社 | プラズマの粒子解析方法及びプラズマの粒子解析プログラムを記録するコンピュータ読取可能な記録媒体 |
CN101719331A (zh) * | 2009-12-08 | 2010-06-02 | 西南石油大学 | 钻井模拟器钻进工艺模拟方法 |
CN105740583A (zh) * | 2016-03-17 | 2016-07-06 | 成都创源油气技术开发有限公司 | 一种适用于碳酸盐岩地层井壁稳定性测井评价方法 |
CN110397431A (zh) * | 2019-05-15 | 2019-11-01 | 上海大学 | 一种钻具接头防断预警方法 |
-
2021
- 2021-03-26 CN CN202110327665.3A patent/CN113047775B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2910759B1 (ja) * | 1998-05-19 | 1999-06-23 | 日本電気株式会社 | プラズマの粒子解析方法及びプラズマの粒子解析プログラムを記録するコンピュータ読取可能な記録媒体 |
CN101719331A (zh) * | 2009-12-08 | 2010-06-02 | 西南石油大学 | 钻井模拟器钻进工艺模拟方法 |
CN105740583A (zh) * | 2016-03-17 | 2016-07-06 | 成都创源油气技术开发有限公司 | 一种适用于碳酸盐岩地层井壁稳定性测井评价方法 |
CN110397431A (zh) * | 2019-05-15 | 2019-11-01 | 上海大学 | 一种钻具接头防断预警方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113047775B (zh) | 2022-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110348046B (zh) | 一种基于模糊综合评价的钻进轨迹多目标优化方法 | |
CN104615862B (zh) | 基于进化算法的高含水油田确定井位的方法 | |
US8005658B2 (en) | Automated field development planning of well and drainage locations | |
AU2017208259B2 (en) | Controlling wellbore operations | |
Chen et al. | Well placement optimization using an analytical formula-based objective function and cat swarm optimization algorithm | |
US20100243328A1 (en) | Continuous geomechanically stable wellbore trajectories | |
CN103177185B (zh) | 一种pdc钻头钻进参数多目标优化方法及装置 | |
CN112462803B (zh) | 一种基于改进nsga-ii的无人机路径规划方法 | |
CN108894768B (zh) | 一种基于蝙蝠算法和井壁稳定的钻进轨迹设计方法与系统 | |
CN110134006A (zh) | 基于改进的多目标粒子群算法的复杂井眼轨迹优化方法 | |
CN105117809B (zh) | 海上丛式井组钻井顺序优化方法 | |
CN109885061B (zh) | 一种基于改进nsga-ⅱ的动力定位多目标优化方法 | |
CN106940742A (zh) | 基于快速自适应量子遗传算法的复杂井眼轨迹优化方法 | |
CN111046593A (zh) | 一种基于钻进时间最短的三维水平井轨道优化设计方法 | |
CN111219175A (zh) | 考虑应力敏感的裂缝性碳酸盐岩酸压裂缝匹配性优化方法 | |
CN113047775B (zh) | 一种井壁稳定性约束下的钻进轨迹多目标优化方法 | |
CN111144053A (zh) | 一种基于遗传算法优化的钻进过程造斜率预测方法及系统 | |
CN105257338B (zh) | 一种储气库盐穴垂直分段溶腔方法 | |
CN113312688B (zh) | 一种页岩气水平井的轨迹设计方法 | |
CN111550184A (zh) | 一种超深水平井井眼轨道设计方法 | |
CN111259506A (zh) | 基于改进遗传算法的车辆发动机机体加工工艺路径优化方法 | |
CN103993829B (zh) | 一种二维小靶前距水平井的中长半径轨道 | |
CN116245236A (zh) | 一种海洋油田水下生产系统多目标生产方案优化方法 | |
CN111625925B (zh) | 一种基于色谱分离的三元复合驱注采优化方法 | |
CN110306968A (zh) | 不规则井网优化方法及其计算机可读存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |