Calculation of Optimal Distance Between Electrode and Probe in Relief Well Magnetic Ranging
-
摘要:
测距精度是影响救援井与事故井连通的关键因素。为提高基于注入电流主动磁测距系统的测量精度,分析了各介质中的电场分布情况,建立了事故井套管电流密度分布模型;通过分析事故井套管电流密度的分布规律,建立了基于测点处磁感应强度最大原则的电极与探管最优距离计算模型,并且通过与试验结果对比验证了模型的有效性。实例计算表明,电极与探管的最优距离不仅与电极到事故井的距离有关,也与相对井斜角有关,但与注入电流的强度无关;事故井套管电流密度的峰值点距事故井坐标点的距离,近似等于电极到事故井探管的距离;电极与探管的最优距离,近似等于电极到事故井距离与相对井斜角余割的乘积。研究结果表明,合理设计电极到探管的距离,有利于提高基于注入电流主动磁测距系统的测量精度。
Abstract:Ranging accuracy is a key factor affecting the successful connection between relief wells and accident wells. In order to improve the measurement accuracy of the active magnetic ranging system based on injected current, the electric field distribution in each medium was analyzed, and the current density distribution model of the accident well casing was established. By analyzing the current density distribution law of accident well casing, the calculation model of the optimal distance between the electrode and the probe based on the principle of the maximum magnetic induction intensity at the measuring point was established and the effectiveness of the model was verified by comparison test results. Calculations with examples show that the optimal distance between the electrode and the probe is related to not only the distance between the electrode and the accident well but also the relative well inclination angle. In addition, it has nothing to do with the injected current intensity. The distance between the peak point of the casing current density of the accident well and the coordinate point of the accident well is approximately equal to that between the electrode and the probe in the accident well. The optimal distance between the electrode and the probe is approximately equal to the product of the distance between the electrode and the accident well and the cosecant of the relative inclination angle. According to the calculation results of the optimal distance, a reasonable design of the distance between the electrode and the probe can contribute to the improvement measurement accuracy of the active magnetic ranging system based on injected current.
-
目前,钻救援井是解决钻井平台井喷、着火等问题的主要方法,精确测控救援井与事故井的相对位置是确保救援井与事故井成功连通的关键[1–3]。传统的测斜工具和测距方法都有累计误差,测距误差不仅随着井深增加而增大,还会因事故井套管的磁干扰而增大[4–6]。磁测距技术具有以下优势:测量的磁信号是具有一定频率的矢量磁信号,可避免其他磁信号的干扰,提高了测量精度;计算的是测点到事故井套管的相对位置关系,没有累计误差的产生,计算精度高;可通过改变电流的强度,增大测距范围;测距时,所有的仪器均在救援井中,不需要在事故井中下入任何探测仪器,特别适用于无法接近事故井井口的救援任务;但磁测距技术难以精确地计算事故井套管上聚集电流的强度,其测量精度会受到影响[7–10]。国外以Wellspot工具为代表的磁测距工具,避免了累计误差的产生,并在救援井导向钻进中得到了推广应用[11–12]。我国在这方面的研究起步较晚,加上国外的技术封锁,使得自主研发该类测距系统的难度进一步加大。
国内在磁测距技术理论和算法上取得了一定的成果[13–18]。基于三电极系救援井与事故井连通探测系统可在一定程度上提高测距精度[19–20];基于径向梯度测量的邻井距离算法,将复杂的物理模型转化为简单的几何模型,但其仅适用近距离探测[21];考虑钻井液和原始地层2种介质对事故井套管上电流分布影响的双层径向介质邻井距离算法,进一步揭示了各介质中电场和电流密度的分布规律[22]。以上均是对测距算法的研究,但其应与现场实际测距有效结合,以进一步提高测距精度。分析可知,探管对测点的磁感应强度有一定的要求,而且磁信号的测量精度直接影响测距系统的测距精度与范围。
笔者提出了电极到探管最优距离的计算方法,意在确定测点的最优位置,使测点的磁感应强度达到最大,从而提高测距精度。为实现探管位置的合理布置,通过分析介质中电场的分布及事故井套管的电流密度分布,揭示了事故井套管周围磁感应强度的分布规律,建立了基于测点处磁感应强度最大原则的电极到探管最优距离的确定方法,并且验证了该方法的正确性。
1. 磁测距原理
磁测距系统主要由地面供电设备、信号采集设备、数据处理系统、地表电极、井下电极和井下探管组成[15–18],如图1所示。测距时,调整好电极到探管的距离,用承载电缆将电极和探管放入救援井至适当位置,注入的低频交变电流在钻井液和地层中传播,并且在事故井套管上聚集。聚集的电流沿套管流动并且在套管周围产生低频的交变磁场。利用探管测量测点的磁感应强度,结合空间几何关系推导出救援井与事故井的相对位置关系,为定向井工程师调整救援井的井眼轨迹提供依据。
2. 模型与方法
2.1 介质中电场强度的计算
电场分布模型见图2。其中,r为电极到事故井套管的距离,m;rr为救援井的井眼半径,m;σm为钻井液的电导率,S/m。
注入的低频交变电流通过钻井液和地层到达事故井套管,电流在事故井套管上聚集并传播。假设钻井液和地层均为均匀介质,具有轴对称性,电位与方位角φ无关;事故井套管的平均半径rc远小于电极到事故井套管的距离r。
介质中电位U的边界条件为:
{∇2U=0U|∞=0Um|r=r+r=Ue|r=r−rσm∂Um∂r|r=r+r=σe∂Ue∂r|r=r−r (1) 式中:U为电位,V;Um为钻井液的电位,V;Ue为地层的电位,V;σe为地层的电导率,S/m。
由式(1)可得:
{∂2U∂r2+1r∂U∂r+∂2U∂z2=0Ue=∫∞0A2(λ)K0(λr)cos(λz)dλI2π 2σm∫∞0K0(λrr)cos(λz)dλ+∫∞0A1(λ)I0(λrr)cos(λz)dλ=∫∞0A2(λ)K0(λrr)cos(λz)dλσm[−I2π 2σm∫∞0K1(λrr)cos(λz)dλ+∫∞0A1(λ)I0(λrr)cos(λz)dλ]=−σe∫∞0A2(λ)K0(λrr)cos(λz)dλ (2) 式中:I0为零阶第一类变形Bessel函数;K0为零阶第二类变形Bessel函数;I1为一阶第一类变形Bessel函数;K1为一阶第二类变形Bessel函数;A1(λ)和A2(λ)为待定系数;I为注入电流的强度,A。
当ρ→0时,钻井液电位的表达式为:
Um=I4π σm1ρ+∫∞0A1(λ)I0(λr)cos(λz)dλ (3) 式中:ρ为测点与电极的距离,m。
由式(2)和式(3)可得钻井液和地层的电位:
{Um=I4π σm[1r + 2π σm−σeσe∫∞0K0(λrr)K1(λrr)I0(λr)I0(λrr)K1(λrr)+σmσeK0(λrr)I1(λrr)cos(λz)dλ]Ue=I2π 2∫∞0K0(λrr)I1(λrr)+K1(λrr)I0(λrr)I1(λrr)σmK0(λrr)I1(λrr)+σeK1(λrr)I0(λrr)K0(λr)cos(λz)dλ (4) 介质中任意测点电场强度的表达式为:
Eρ = −∇Ue(r,z) (5) 当z = 0时,r方向上的电场强度Er可以用方程表示为:
Er={I2π 2σm∫∞0[K1(λr)−σm−σeσeK0(λrr)K1(λrr)I1(λr)I0(λrr)K1(λrr)+σmσeK0(λrr)I1(λrr)]λdλ (r⩽ (6) 地层中,r方向上的电场强度E0为:
E_0=\frac{I}{2\text{π }^2}\int_0^{\infty}\frac{\text{K}_0\left(\lambda r_{\mathrm{r}}\right)\text{I}_1\left(\lambda r\mathrm{_r}\right)+\text{K}_1\left(\lambda r\mathrm{_r}\right)\text{I}_0\left(\lambda r_{\mathrm{r}}\right)\text{I}_1\left(\lambda r_{\mathrm{r}}\right)}{\sigma_{\text{m}}\text{K}_0\left(\lambda r\mathrm{_r}\right)\text{I}_1\left(\lambda r\mathrm{_r}\right)+\sigma_{\text{e}}\text{K}_1\left(\lambda r_{\mathrm{r}}\right)\text{I}_0\left(\lambda r\mathrm{_r}\right)}\text{K}_1\left(\lambda r\right)\lambda d\lambda\ \ \ \ \left(r > r\mathrm{_r}\right) (7) 2.2 事故井套管电场强度的计算
假设地层为均匀介质,事故井套管无限长,事故井套管的平均半径rc远小于电极到事故井套管的距离r,空心的事故井套管可用具有相同传导面积的均匀圆柱体代替[11],则注入电流在事故井套管的传播如图3所示。各参数间的关系式为:
2{\text{π}}{r_{\text{c}}}{h_{\text{c}}}{\sigma _{\text{c}}}{{ = {\text{π}} }}{r^2}_{{\text{cz}}}{\sigma _{\text{c}}} (8) 式中:σc为套管的电导率,S/m;hc为事故井套管壁厚度,m;rcz为等效半径,m;
以电极在事故井轴线上的镜像位置为坐标原点,事故井套管轴线为z轴建立柱状坐标系,Δr为远离圆柱体轴线的径向距离,则地层中电位满足:
{\nabla ^2}U = 0 (9) 地层中电场强度满足:
E{\text{ = }} - \nabla U\left( {\Delta r,\varphi } \right) (10) 电位和电流密度在
\Delta r = {r_{{\mathrm{cz}}}} 处连续。当Δr=0时,电场的表达式为[11]:\frac{{E\left( z \right)}}{{{E_0}}} = \frac{1}{{\text{π }}}\int_0^\infty {\frac{{u{{\text{K}}_0}\left( u \right)\sin {\dfrac{{uz}}{r}} }}{{1 - \dfrac{1}{4}\dfrac{{{\sigma _{\mathrm{c}}}{r_{{\mathrm{cz}}}}^2}}{{{\sigma _{\mathrm{e}}}{r^2}}}{u^2}\ln \dfrac{{u{r_{{\mathrm{cz}}}}}}{r}}}}{\text{d}}u (11) 根据式(7)和式(11),得到事故井套管上的电场强度为:
E\left( z \right) = {E_0}\frac{1}{{\text{π }}}\int_0^\infty {\dfrac{{u{{\text{K}}_0}\left( u \right)\sin {\dfrac{{uz}}{r}} }}{{1 - \dfrac{1}{4}\dfrac{{{\sigma _{\mathrm{c}}}{r_{{\mathrm{cz}}}}^2}}{{{\sigma _{\mathrm{e}}}{r^2}}}{u^2}\ln \dfrac{{u{r_{{\mathrm{cz}}}}}}{r}}}}{\text{d}}u (12) 式中:u=λr;E(z)为事故井套管上的电场强度,V/m。
2.3 电极与探管最优距离的确定
事故井套管上的电流密度为:
J\left( z \right) = E\left( z \right){\sigma _{\text{c}}} (13) 事故井套管上的电流强度为:
I\left( z \right) = E\left( z \right){\text{π }}{r_{{\text{cz}}}}^2{\sigma _{\text{c}}} (14) 事故井套管周围的磁感应强度可以近似为:
B\left( z \right) = {\mu _0}\frac{{I\left( z \right)}}{{2{\text{π }}R}} (15) 式中:J(z)为事故井套管上的电流密度,A/m2;I(z)为事故井套管上的电流强度,A; B(z)为事故井套管周围的磁感应强度,T;μ0为真空磁导率,H/m;R为测点与事故井套管的距离,m。
通过分析事故井套管电流密度的分布规律,可得电极到事故井套管的距离与电极到探管距离的关系,进而根据测点磁感应强度最大的原则,确定电极与探管的最优距离为:
L\cos \alpha = z\left( r \right) (16) z\left( r \right) \approx r (17) 式中:L为电极与探管的距离,m;α为相对井斜角,即事故井井斜角与救援井井斜角的差值,计算时取差值的绝对值。
2.4 模型验证
基于低频注入电流的磁测距原理,李翠等人[15]进行了邻井距离的模拟试验。模拟井组包括裸眼井和套管井,裸眼井和套管井的井口间距为5 m,注入电流的强度为13~23 A,频率为0.25 Hz。本文模型所涉及的参数与文献[15]中的试验参数基本相同。为了验证模型的有效性,笔者将本文模型的计算结果与文献[15]中的试验结果进行了比较。
利用本文模型,计算了电极与探管距离L为10,15和20 m时不同注入电流强度下裸眼井和套管井井口的间距,并与文献[15]中的试验结果进行比较(见图4)。从图4可以看出:电极与探管距离一定时,测距误差随着注入电流强度增大而减小;注入电流强度一定时,测距误差随着电极与探管距离增大而减小;模型计算结果的趋势与试验结果相同。
与试验结果相比,模型计算结果的误差更小。当L=20 m时,测距的相对误差在5.92%~7.23%,平均相对误差为6.58%。试验结果误差偏高的原因是:计算邻井距离的模型仅考虑了电流在地层中的传播规律,未考虑钻井液对电流的影响;试验现场的地层并非均质,而模型中的地层为均质地层;测量工具本身的误差不可避免。以上分析可知:理论模型中考虑了钻井液和地层双层介质对电流传播的影响,与实际情况接近,验证了模型的有效性。
3. 算例分析
根据文献[11]提供的地层和套管的参数,计算事故井套管的电流密度,确定电极到探管的最优距离,分析注入电流、电极与事故井距离和相对井斜角对电极与探管最优距离的影响。部分原始数据为:地层电导率1 S/m,套管电导率107 S/m,钻井液电导率200 S/m,套管平均半径0.125 m,套管壁厚0.013 m,救援井井眼半径0.200 m,注入电流强度20 A,真空磁导率4π×10−7 H/m。计算结果如图5所示。
从图5可以看出,当电极与事故井探管的距离为30 m、相对井斜角为10°时,磁感应强度的峰值(B=0.15 nT)大约在±31 m处,此时电极与探管的最优距离为31.48 m。可知,事故井套管电流密度的峰值点距事故井坐标点的距离近似等于电极与事故井探管的距离。
3.1 注入电流强度的影响
图6所示为事故井套管电流密度随注入电流强度变化的计算结果。由图6可知,随注入电流强度增大,事故井套管的电流密度不断增大。增大注入电流的强度从整体上增大了测点的磁感应强度,但不影响电极与探管最优距离的确定。此时,电极与探管的最优距离仍是31 m。
3.2 电极与事故井距离的影响
图7所示为事故井套管电流密度随电极与事故井距离变化的计算结果。由图7可知,随着电极与事故井距离增大,事故井套管上电流的密度不断减小。电极与事故井距离增大不仅导致测点的磁感应强度减小,还影响电极与探管的最优距离。电极与探管的最优距离随电极与事故井距离增大而增大。
3.3 相对井斜角的影响
图8所示为电极到探管的最优距离随相对井斜角变化的计算结果。由图8可知,电极到探管的最优距离随相对井斜角增大而增大。
4. 结 论
1)钻救援井过程中,采用磁测距技术测量救援井与事故井距离时,探管在救援井中的位置将影响测量精度与范围,合理设计电极与探管的距离有利于提高测量精度、增大测距范围。
2)电极与探管的最优距离,不仅与电极与事故井的距离有关,也与相对井斜角有关,但是与注入电流的强度无关。随着电极与事故井距离增大,事故井套管上电流的密度不断减小,电极与探管的最优距离不断增大;电极与探管的最优距离随相对井斜角增大而增大。
3)事故井套管电流密度的峰值点距事故井坐标点的距离,近似等于电极与事故井探管的距离;电极与探管的最优距离,近似等于电极与事故井距离与相对井斜角余割的乘积。
-
-
[1] FLORES V, DAILEY P, TODD D, et al. Relief well planning[R]. SPE 168029, 2014.
[2] 车阳,乔磊,王建国,等. X 井精准治理技术研究与应用[J]. 断块油气田,2023,30(2):347–352. CHE Yang, QIAO Lei, WANG Jianguo, et al. Research and application of precise treatment technology for well X[J]. Fault-Block Oil & Gas Field, 2023, 30(2): 347–352.
[3] 刁斌斌,高德利,唐海雄,等. 救援井与事故井邻井距离探测技术[C]//第十六届全国探矿工程(岩土钻掘工程)技术学术交流年会论文集. 北京:地质出版社,2011:201-205. DIAO Binbin, GAO Deli, TANG Haixiong, et al. The distance detection technology between rescue wells and accident wells [C]//Proceedings of the 16th National Academic Exchange Conference on Exploration Engineering (Rock and Soil Drilling Engineering) Technology. Beijing: Geological Publishing House, 2011: 201-205.
[4] WOLFF C J M, DE WARDT J P. Borehole position uncertainty-analysis of measuring methods and derivation of systematic error model[J]. Journal of Petroleum Technology, 1981, 33(12): 2338–2350. doi: 10.2118/9223-PA
[5] WILLIAMSON H S. Accuracy prediction for directional measurement while drilling[J]. SPE Drilling & Completion, 2000, 15(4): 221–233.
[6] TORKILDSEN T, HÅVARDSTEIN S T, WESTON J, et al. Prediction of wellbore position accuracy when surveyed with gyroscopic tools[J]. SPE Drilling & Completion, 2008, 23(1): 5–12.
[7] DIAO Binbin, GAO Deli. Study on a ranging system based on dual solenoid assemblies, for determining the relative position of two adjacent wells[J]. Computer Modeling in Engineering & Sciences, 2013, 90(1): 77–90.
[8] KUCKES A F, HAY R T, MCMAHON J, et al. New electromagnetic surveying/ranging method for drilling parallel horizontal twin wells[J]. SPE Drilling & Completion, 1996, 11(2): 85–90.
[9] 高德利,刁斌斌. 复杂结构井磁导向钻井技术进展[J]. 石油钻探技术,2016,44(5):1–9. GAO Deli, DIAO Binbin. Development of the magnetic guidance drilling technique in complex well engineering[J]. Petroleum Drilling Techniques, 2016, 44(5): 1–9.
[10] 李翠,高德利,刘庆龙,等. 邻井随钻电磁测距防碰计算方法研究[J]. 石油钻探技术,2016,44(5):52–59. LI Cui, GAO Deli, LIU Qinglong, et al. A method of calculating of avoiding collisions with adjacent wells using electromagnetic ranging surveying while drilling tools[J]. Petroleum Drilling Technique, 2016, 44(5): 52–59.
[11] WEST C L, KUCKES A F, RITCH H J. Successful ELREC logging for casing proximity in an offshore Louisiana blowout[R]. SPE 11996, 1983.
[12] KUCKES A F, LAUTZENHISER T, NEKUT A G, et al. An electromagnetic survey method for directionally drilling a relief well into a blown out oil or gas well[J]. SPE Journal, 1984, 24(3): 269–274.
[13] 李翠,高丽萍,李佳,等. 邻井随钻电磁测距防碰工具模拟试验研究[J]. 石油钻探技术,2017,45(6):110–115. LI Cui, GAO Liping, LI Jia, et al. Experiment research on an electromagnetic anti-collision detection tool while drilling adjacent wells[J]. Petroleum Drilling Techniques, 2017, 45(6): 110–115.
[14] DIAO Binbin, GAO Deli, LI Genkui. Development of static magnetic detection anti-collision system while drilling[C]//Proceedings of the 2016 International Conference on Artificial Intelligence and Engineering Applications. Amsterdam: Atlantis Press, 2016: 543-551.
[15] WU Zhiyong, GAO Deli, DIAO Binbin. An investigation of electromagnetic anti-collision real-time measurement for drilling cluster wells[J]. Journal of Natural Gas Science and Engineering, 2015, 23: 346–355. doi: 10.1016/j.jngse.2015.02.016
[16] 刁斌斌,高德利. 邻井定向分离系数计算方法[J]. 石油钻探技术,2012,40(1):22–27. doi: 10.3969/j.issn.1001-0890.2012.01.005 DIAO Binbin, GAO Deli. Calculation method of adjacent well oriented separation factors[J]. Petroleum Drilling Techniques, 2012, 40(1): 22–27. doi: 10.3969/j.issn.1001-0890.2012.01.005
[17] DOU Xinyu, LIANG Huaqing, LIU Yang. Anticollision method of active magnetic guidance ranging for cluster wells[J]. Mathematical Problems in Engineering, 2018, 2018: 7583425.
[18] INGRAM S R, LAHMAN M, PERSAC S. Methods improve stimulation efficiency of perforation clusters in completions[J]. Journal of Petroleum Technology, 2014, 66(4): 32–36. doi: 10.2118/0414-0032-JPT
[19] 李翠,高德利. 救援井与事故井连通探测方法初步研究[J]. 石油钻探技术,2013,41(3):56–61. LI Cui, GAO Deli. Preliminary research on detection method for connecting relief well to blowout well[J]. Petroleum Drilling Techniques, 2013, 41(3): 56–61.
[20] LI Cui, GAO Deli, WU Zhiyong, et al. A method for the detection of the distance & orientation of the relief well to a blowout well in offshore drilling[J]. Computer Modeling in Engineering and Sciences, 2012, 89(1): 39–56.
[21] ZHANG Sen, DIAO Binbin, GAO Deli. A new method of anti-collision while drilling based on radial gradient measurement[R]. ARMA-2019-0114, 2019.
[22] ZHANG Sen, DIAO Binbin, GAO Deli. Numerical simulation and sensitivity analysis of accurate ranging of adjacent wells while drilling[J]. Journal of Petroleum Science and Engineering, 2020, 195: 107536. doi: 10.1016/j.petrol.2020.107536