-
近年来海洋油气勘探进入高速发展阶段[1],人工岛是实现“海油陆采”的一种高效开发模式。以冀东油田为例,已先后开发高尚堡、柳赞、老爷庙、唐海、南堡5座人工岛。其中南堡1-3D人工岛占地0.13 km2,截止2018年底已完钻井200口以上,但开发中后期采油速度大幅降低,为完成生产任务需在人工岛上进行丛式井组加密调整。目前地面可利用空间和井口十分有限,地下已完钻密集井段安全穿越也愈发困难,新建加密井数量越多相邻井眼之间的距离越小,在钻进时出现井眼交碰的概率就越大,任何非预期的井眼交碰都会给油田造成严重损失,因此防碰风险评价是密集井网加密井设计一个必不可少的环节[2-4]。常用的井眼交碰风险评价方法有井眼分离系数法、井眼交碰概率法以及井眼最近距离法[5-7],这些方法的共同点是评价体系中只有一口比较井,而密集井网条件下在防碰分析范围内的井眼往往存在多口,要规避其中任意井,势必造成与其他井的交碰风险增加。可见,在密集井网条件下采用现有方法进行井眼交碰风险评价不够全面准确。为此,根据密集井网加密井设计的特点,提出了一种可定量评价一口基准井与多口比较井的密集井网加密井井眼交碰综合风险评价新方法。
-
1996年Brooks等首次提出井眼交碰概率的概念并给出了数学计算方法[8],该方法通过目标点和参考点的协方差矩阵计算井眼交碰概率,可定量评价2个误差椭球的分离程度,但是该方法只考虑到目标点和参考点的位置不确定性矩阵,忽略了两井眼名义间距对交碰概率的影响,评价结果难以满足防碰钻井要求。
笔者改进了Brooks的交碰概率模型,综合考虑名义间距、基准井和比较井的井径之和以及邻井数等因素,建立了一个比较全面的密集丛式井网井眼交碰风险计算方法。通过几何概率模型计算井眼偏向某邻井方向的概率θi,运用Brooks的井眼交碰概率模型计算相邻井眼的交碰概率Pi,井眼实际交碰风险值为θi与Pi乘积,则密集井网井眼交碰综合风险值为基准井与各比较井井眼实际交碰风险值的和,和值越小井眼综合交碰风险程度越低,否则越高。图1为密集井网某一深度基准井、比较井相对位置示意图,图中O为基准井,O1、O2、O3为比较井,d1为基准井井径,d2为比较井井径,S1、S2、S3分别为基准井到各比较井的名义间距,则基准井向O1方向倾斜的概率等于长度为d1+d2的弧长与半径为S1的圆周周长的比值,记为θ1,运用Brooks交碰概率计算模型可计算得到基准井O与比较井O1的交碰概率P1,同理可以得到θ2、θ3、P2、P3,则图示密集井网的综合交碰风险值为
$\displaystyle\sum\nolimits_{i = 1}^3 {{\theta _i}} {P_i}$ 。 -
受地层非均质性、钻具组合、钻头类型等因素影响,钻井过程中会出现方位偏移,偏移方向具有随机性但服从几何概率分布模型,基准井与某一比较井的相对位置如图2所示,图中O为基准井、O1为比较井,基准井、比较井井眼直径分别为d1、d2,名义间距为S,∠AOB为α,则名义间距、井径之和与交碰概率的关系如式(1)~(3)所示。
$$S = \sqrt {{{\left( {{x_o} - {x_r}} \right)}^2} + {{\left( {{y_o} - {y_r}} \right)}^2} + {{\left( {{{\textit{z}}_o} - {{\textit{z}}_r}} \right)}^2}} $$ (1) 式中,
$\left( {{x_o},{y_o},{{\textit{z}}_o}} \right)$ 为目标点坐标,$\left( {{x_r},{y_r},{{\textit{z}}_r}} \right)$ 为参考点坐标。由图2可知,当基准井向AOB区域偏移且偏移距离大于S−0.5(d1+d2)时两口井发生交碰,当基准井向其他区域偏移时,则与比较井O1发生交碰事故的概率为0。根据几何代数AOB区域满足如下关系
$$\sin \alpha = \frac{{{d_1} + {d_2}}}{S}$$ (2) 由几何概率定义,某事件发生的概率等于该事件的几何尺寸与总几何尺寸的比值,则基准井向AOB区域偏移的概率为
$$\theta = \frac{\alpha }{{360{\rm{^\circ }}}} = \frac{{\arcsin \left[ {{{\left( {{d_1} + {d_2}} \right)}/S}} \right]}}{{360{\rm{^\circ }}}}$$ (3) -
误差源是井眼轨迹误差分析的基本要素,主要分为随机误差源、系统误差源,任一测点处的井眼位置不确定性矩阵可由不同误差源贡献累加获得,ISCWSA提供了不同误差传播模式下的协方差矩阵计算公式[9-11]。假设通过井眼轨迹位置不确定性分析方法,得到目标点和参考点的位置不确定矩阵C1、C2,其中不确定矩阵C1、C2相互独立,则两个矩阵之和表示目标点和参考点位置的相对不确定性,记为C。
$${{{C}}_1} = \left[ \begin{gathered} \begin{array}{*{20}{c}} {\sigma _{1xx}^2}&{{\sigma _{1xy}}}&{{\sigma _{1x{\textit{z}}}}} \end{array} \\ \begin{array}{*{20}{c}} {{\sigma _{1yx}}}&{\sigma _{1yy}^2}&{{\sigma _{1y{\textit{z}}}}} \end{array} \\ \begin{array}{*{20}{c}} {{\sigma _{1{\textit{z}}x}}}&{{\sigma _{1{\textit{z}}y}}}&{\sigma _{1{\textit{z}}{\textit{z}}}^2} \end{array} \\ \end{gathered} \right]$$ (4) $${{{C}}_2} = \left[ \begin{gathered} \begin{array}{*{20}{c}} {\sigma _{2xx}^2}&{{\sigma _{2xy}}}&{{\sigma _{2x{\textit{z}}}}} \end{array} \\ \begin{array}{*{20}{c}} {{\sigma _{2yx}}}&{\sigma _{2yy}^2}&{{\sigma _{2y{\textit{z}}}}} \end{array} \\ \begin{array}{*{20}{c}} {{\sigma _{2{\textit{z}}x}}}&{{\sigma _{2{\textit{z}}y}}}&{\sigma _{2{\textit{z}}{\textit{z}}}^2} \end{array} \\ \end{gathered} \right]$$ (5) $${{C}} = {{{C}}_1} + {{{C}}_2} = \left[ \begin{gathered} \begin{array}{*{20}{c}} {\sigma _{xx}^2}&{{\sigma _{xy}}}&{{\sigma _{x{\textit{z}}}}} \end{array} \\ \begin{array}{*{20}{c}} {{\sigma _{yx}}}&{\sigma _{yy}^2}&{{\sigma _{y{\textit{z}}}}} \end{array} \\ \begin{array}{*{20}{c}} {{\sigma _{{\textit{z}}x}}}&{{\sigma _{{\textit{z}}y}}}&{\sigma _{{\textit{z}}{\textit{z}}}^2} \end{array} \\ \end{gathered} \right]$$ (6) 式中,σxx、σyy、σzz为坐标轴上的方差,其余各项为坐标轴间的协方差。
在井眼轨迹测量领域通常假定井眼轨迹位置坐标服从正态分布,根据概率论与数理统计学,三维空间某点位置坐标的正态分布概率密度函数为[12]
$$f\left( {x,y,{\textit{z}}} \right) = \frac{{\exp \left[ {{{ - \left( {{{{r}}^{\rm{T}}}{{{C}}^{ - 1}}{{r}}} \right)}/{\rm{2}}}} \right]}}{{{{\left( {2{\rm{\pi }}} \right)}^{1.5}}\sqrt \Delta }}$$ (7) 式中,r为某点坐标位置增量,C−1为井眼不确定矩阵的逆矩阵,∆为相对位置不确定性矩阵的行列式,该不确定性矩阵是坐标位置的函数,与坐标位置增量无关。当坐标位置确定时,式(7)分母为一个常数,则该测点处三维空间正态分布概率密度值相等的点满足下式
$${{{r}}^{\rm{T}}}{{{C}}^{ - 1}}{{r}} = {k^2}$$ (8) 式(8)为一误差椭球族,其中k为放大因子,表示井眼位于椭球区域的概率。
目标点是基准井上最接近参考点的点,目标井的局部方向与两点连线正交,井眼交碰定义为目标点与参考点间的距离小于等于两口井井眼半径之和,若两口井井眼半径之和与井眼名义间距数量级相差不大,则井眼交碰概率表示为
$$P = \frac{{\rm{\pi }}}{4}{\left( {{d_{\rm{o}}} + {d_{\rm{r}}}} \right)^2}\int\limits_{ - \infty }^\infty {f\left( {x,y,{\textit{z}}} \right){\rm{d}}L} $$ (9) 式中,do为基准井井径,dr为比较井井径,f (x,y,z)为三维空间某点位置坐标的正态分布概率密度函数,dL为井眼微元段。
目标井眼局部为直线,则坐标位置增量可由一点ro和单位向量ν线性表示,即
$$r = {{{r}}_{\rm{o}}} + L{{\nu}} $$ (10) 联立式(7)、(10)可得
$$f\left( {x,y,{\textit{z}}} \right){\rm{ = }}\dfrac{{\exp \left[ { - \dfrac{1}{2}\left( {\gamma - \dfrac{{{\beta ^2}}}{\alpha }} \right)} \right]}}{{{{\left( {{\rm{2\pi }}} \right)}^{{\rm{1}}{\rm{.5}}}}\sqrt \Delta }}\exp \left[ { - \dfrac{1}{2}{{\left( {L\sqrt \alpha + \dfrac{\beta }{{\sqrt \alpha }}} \right)}^2}} \right]$$ (11) 其中
$$ \alpha = {{{\nu}} ^{\rm{T}}}{{{C}}^{ - 1}}{{\nu}} ,\;\beta = {{r}}_{\rm{o}}^{\rm{T}}{{{C}}^{ - 1}}{{\nu}},\; \gamma = {{r}}_{\rm{o}}^{\rm{T}}{{{C}}^{ - 1}}{{{r}}_{\rm{o}}} $$ 将式(11)代入式(9)化简整理得
$$P = \frac{{{{\left( {{d_0} + {d_{\rm{r}}}} \right)}^2}}}{{8\sqrt {\alpha \Delta } }}\exp \left[ { - \frac{1}{2}\left( {\gamma - \frac{{{\beta ^2}}}{\alpha }} \right)} \right]$$ (12) -
通过上述计算方法得到目标井与各比较井的θi和Pi,则密集井网井眼交碰综合风险值(wellbore collision comprehensive probability)为
$$\begin{gathered} {R_{{\rm{WCCP}}}} = {\theta _1}{P_1} + {\theta _2}{P_2} + {\theta _3}{P_3}{\rm{ + }} \cdot \cdot \cdot + {\theta _n}{P_n} = \sum\limits_{i = 1}^n {{\theta _i}} {P_i} \end{gathered} $$ (13) 式中,θi为目标井偏向第i口比较井的概率;Pi为目标井与第i口比较井交碰概率;n为密集丛式井防碰分析区域内的邻井数。
基准井井眼偏移概率、井眼交碰概率均表示事件发生的可能性,取值范围均为(0,1)。由《密集丛式井上部井段防碰设计与施工技术规范》(Q/SY 1296—2010)防碰区域划分标准,密集丛式井直井段防碰分析区域为以基准井井眼中心轴线为圆心、半径为15 m的空间,当两井眼中心轴线距离超过15 m时为安全范围,无需进行防碰分析,同时根据开发调整井网钻井井眼轨道防碰控制要求[13-14],设计井深在同一防碰井段两井眼中心轴线距离小于15 m的井不能多于3口,否则即为高交碰风险井,因此RWCCP的取值范围为(0,3)。针对密集丛式井网现场施工实际,综合中国石油天然气集团公司生产安全风险防控管理办法及相关专家意见,依据RWCCP划分密集井网井眼交碰风险评价等级,制定防碰风险量化表,如表1所示。
表 1 密集井网井眼防碰风险分析量化表
Table 1. Quantitative wellbore collision risk analysis of dense cluster wells
RWCCP 交碰等级 交碰风险 预防措施 10−3≤RWCCP<3 一级 重大风险 停钻直到风险消除 10−4≤RWCCP<10−3 二级 较小风险 关闭邻井制定防碰措施 0<RWCCP<10−4 三级 警戒 时刻监测邻井情况 -
冀东南堡3号人工岛某设计调整加密井部分测斜数据如表2所示,将该轨迹数据导入Campass采用最近距离扫描,扫描结果显示与邻井的分离系数均大于1,但按照设计进行施工在井深470 m处却发生了非预期的井眼交碰,给油田带来较大的经济损失。
表 2 调整加密井部分测斜数据
Table 2. Partial directional data of adjustment infill wells
测点 垂深/m 井斜角/(°) 方位角/(°) N/m E/m 1 0.00 0.00 0.00 0.00 0.00 2 39.93 1.00 230.58 −0.22 −0.27 3 68.49 0.91 223.84 −0.54 −0.62 4 97.14 0.67 212.47 −0.85 −0.87 5 125.73 0.61 240.17 −1.07 −1.09 6 154.32 0.50 259.84 −1.16 −1.34 7 182.79 1.00 46.19 −1.01 −1.29 8 211.34 0.88 28.97 −0.65 −1.00 9 239.92 1.25 44.76 −0.24 −0.67 10 268.51 1.25 49.32 0.19 −0.22 11 285.97 0.83 108.46 0.27 0.05 12 313.97 0.98 102.70 0.16 0.47 13 326.33 1.05 98.13 0.12 0.69 14 354.91 1.76 145.44 −0.28 1.20 15 383.48 2.35 181.67 −1.23 1.43 16 411.91 3.96 181.14 −2.79 1.39 17 440.35 6.25 175.56 −5.33 1.49 18 468.69 8.65 173.30 −9.02 1.86 19 496.78 12.72 175.49 −14.30 2.35 20 524.57 14.20 180.27 −20.94 2.58 调整加密井井深470 m处在密集丛式井网防碰分析范围的邻井有2口,分别为#1井、#2井,该深度处#1井、#2井的轨迹基本参数如表3所示。通过井眼轨迹不确定分析方法得到井眼直径为244.5 mm的加密井、#1井以及#2井协方差矩阵。
表 3 井深470 m处邻井井身轨迹基本参数
Table 3. Basic wellbore trajectory parameters of the neighboring wells at the well depth of 470 m
井号 垂深/m 井斜角/(°) 方位角/(°) N/m E/m #1 472 10.38 216.53 −7.22 2.12 #2 469 11.73 157.26 −1.21 3.14 $$ {{{C}}}_{{\text{加密井}}}\rm=\left[\begin{array}{ccc}4.18& 0.22& 0.27\\ 0.22& 6.08& \rm-0.03\\ 0.27& \rm-0.03& 0.53\end{array}\right]$$ (14) $${{{C}}_{{\rm{\# }}1}}{\rm{ = }}\left[ {\begin{array}{*{20}{c}} {6.06}&{{\rm{ - }}1.65}&{0.35} \\ {{\rm{ - }}1.65}&{7.08}&{0.35} \\ {0.35}&{{\rm{ - }}0.03}&{0.57} \end{array}} \right]$$ (15) $${{{C}}_{{\rm{\# }}2}}{\rm{ = }}\left[ {\begin{array}{*{20}{c}} {5.93}&{1.43}&{0.51} \\ {1.43}&{8.74}&{{\rm{ - }}0.23} \\ {0.51}&{{\rm{ - }}0.23}&{0.61} \end{array}} \right]$$ (16) 采用本文的密集井网加密井井眼综合交碰风险计算方法对该设计井进行风险分析,470 m处井眼交碰综合风险值为1.019×10−4 (二级风险),此时应关闭邻井制定防碰措施,便可规避该交碰事故的发生,可见该方法对于指导密集井网加密井防碰施工具有重要意义。
-
(1)井眼分离系数法、井眼交碰概率法和井眼最近距离法不适用于密集井网条件下的井眼防碰分析。
(2)改进了Brooks的交碰概率计算模型,交碰概率计算应综合考虑名义间距、井径以及邻井数等因素,在此基础上提出了一个比较全面的密集丛式井网井眼交碰风险计算方法,并给出了密集井网井眼防碰风险量化等级标准。
(3)该方法可定量评价一口基准井与多口比较井之间的综合交碰风险,对加密井钻井防碰理论研究及指导生产实践具有重要意义。
A new method for calculating the wellbore collision risk of vertical sections of dense cluster wells
-
摘要: 井眼防碰是密集丛式井网定向井施工面临的一个主要难题,为了有效减小密集井网加密井直井段井眼交碰风险,综合考虑名义间距、基准井和比较井的井径之和以及邻井数等因素,在分析Brooks运用概率评价井眼交碰方法的基础上,提出了比较全面的密集丛式井网井眼交碰风险计算新方法。该方法以基准井井眼中心轴线为圆心、基准井与比较井间名义间距为半径作圆,通过基准井、比较井井径之和与此圆周长的比值计算偏移概率,采用井眼测点位置坐标的正态分布概率密度函数沿井深方向的定积分计算交碰概率,基准井与比较井的实际交碰概率为偏移概率与交碰概率的乘积,通过基准井与各比较井实际交碰概率求和确定综合交碰风险,并将综合交碰风险按照行业标准划分为3个风险等级,得到以综合交碰风险为评价指标的风险评价方法。该方法可定量评价一口基准井与多口比较井之间的综合交碰风险,对加密井钻井防碰理论研究及指导生产实践具有重要意义。Abstract: Wellbore anti-collision is an important difficulty in the construction of directional wells of dense cluster wells. In order to effectively reduce the collision risks of the vertical sections of the infill wells of dense cluster wells, a new comprehensive method for calculating the wellbore collision risks of dense cluster wells was put forward by considering nominal spacing, diameter sum of reference well and comparison well and number of neighboring wells comprehensively after the probability based Brooks wellbore collision evaluation method was analyzed. In this new method, a circle is set by taking the central axis of the wellbore of reference well as the center and the nominal spacing between reference well and comparison well as the radius. Migration probability is calculated by using the ratio of the diameter sum of reference well and comparison well to the circumference of the circle, and the collision probability is calculated by using the definite integral of the probability density function of normal distribution at the coordinate of wellbore measurement point in the direction of well depth. The actual collision probability between reference well and comparison well is the product of the migration probability and the collision probability. The comprehensive collision probability is determined by summing the collision probability between the reference well and each comparison well, and it is classified into three risk levels according to industry standards. In this way, a collision risk evaluation method that takes the comprehensive collision risk as an evaluation index is developed. This method can quantitatively evaluate the comprehensive collision risk between a reference well and multiple comparison wells, and is of important significance to anti-collision theory research and production practice guidance of infill wells.
-
Key words:
- dense cluster well /
- infill well /
- collision probability /
- calculation method
-
表 1 密集井网井眼防碰风险分析量化表
Table 1. Quantitative wellbore collision risk analysis of dense cluster wells
RWCCP 交碰等级 交碰风险 预防措施 10−3≤RWCCP<3 一级 重大风险 停钻直到风险消除 10−4≤RWCCP<10−3 二级 较小风险 关闭邻井制定防碰措施 0<RWCCP<10−4 三级 警戒 时刻监测邻井情况 表 2 调整加密井部分测斜数据
Table 2. Partial directional data of adjustment infill wells
测点 垂深/m 井斜角/(°) 方位角/(°) N/m E/m 1 0.00 0.00 0.00 0.00 0.00 2 39.93 1.00 230.58 −0.22 −0.27 3 68.49 0.91 223.84 −0.54 −0.62 4 97.14 0.67 212.47 −0.85 −0.87 5 125.73 0.61 240.17 −1.07 −1.09 6 154.32 0.50 259.84 −1.16 −1.34 7 182.79 1.00 46.19 −1.01 −1.29 8 211.34 0.88 28.97 −0.65 −1.00 9 239.92 1.25 44.76 −0.24 −0.67 10 268.51 1.25 49.32 0.19 −0.22 11 285.97 0.83 108.46 0.27 0.05 12 313.97 0.98 102.70 0.16 0.47 13 326.33 1.05 98.13 0.12 0.69 14 354.91 1.76 145.44 −0.28 1.20 15 383.48 2.35 181.67 −1.23 1.43 16 411.91 3.96 181.14 −2.79 1.39 17 440.35 6.25 175.56 −5.33 1.49 18 468.69 8.65 173.30 −9.02 1.86 19 496.78 12.72 175.49 −14.30 2.35 20 524.57 14.20 180.27 −20.94 2.58 表 3 井深470 m处邻井井身轨迹基本参数
Table 3. Basic wellbore trajectory parameters of the neighboring wells at the well depth of 470 m
井号 垂深/m 井斜角/(°) 方位角/(°) N/m E/m #1 472 10.38 216.53 −7.22 2.12 #2 469 11.73 157.26 −1.21 3.14 -
[1] 许军富, 徐文浩, 耿应春. 渤海人工岛大型丛式井组加密防碰优化设计技术[J]. 石油钻探技术, 2018, 46(2):24-29. doi: 10.11911/syztjs.2018048 XU Junfu, XU Wenhao, GENG Yingchun. Anti-Collision optimization design technology for Large-Scale infill drilling for cluster well groups in the artificial island of the Bohai sea[J]. Petroleum Drilling Techniques, 2018, 46(2): 24-29. doi: 10.11911/syztjs.2018048 [2] 李兵, 胥豪, 牛洪波, 等. 丛式井组加密井防碰技术及应用[J]. 钻采工艺, 2018, 41(4):12-15. doi: 10.3969/J.ISSN.1006-768X.2018.04.04 LI Bing, XU Hao, NIU Hongbo, et al. Anti-collision technology for infill well drilling[J]. Drilling & Production Technology, 2018, 41(4): 12-15. doi: 10.3969/J.ISSN.1006-768X.2018.04.04 [3] 刘永旺, 管志川, 王伟, 等. 丛式井组直井段交碰风险评价及设计优化[J]. 中国安全生产科学技术, 2015, 11(10):85-89. doi: 10.11731/j.issn.1673-193x.2015.10.015 LIU Yongwang, GUAN Zhichuan, WANG Wei, et al. Collision risk evaluation and design optimization on vertical section of well in cluster wells[J]. Journal of Safety Science and Technology, 2015, 11(10): 85-89. doi: 10.11731/j.issn.1673-193x.2015.10.015 [4] 边瑞超, 周洪林, 曹华庆. 冀东油田人工岛丛式井钻井防碰技术[J]. 石油钻探技术, 2017, 45(5):19-22. BIAN Ruichao, ZHOU Honglin, CAO Huaqing. Anti-collision technology for drilling cluster wells in the artificial island of Jidong Oilfield[J]. Petroleum Drilling Techniques, 2017, 45(5): 19-22. [5] 柳贡慧, 董本京, 高德利. 误差椭球(圆)及井眼交碰概率分析[J]. 钻采工艺, 2000, 23(3):5-12. doi: 10.3969/j.issn.1006-768X.2000.03.002 LIU Gonghui, DONG Benjing, GAO Deli. Probability analysis of error ellipsoid (ellipse) and hole intersection[J]. Drilling & Production Technology, 2000, 23(3): 5-12. doi: 10.3969/j.issn.1006-768X.2000.03.002 [6] 史玉才, 张晨, 薛磊, 等. 井眼分离系数计算新方法[J]. 石油学报, 2015, 36(12):1580-1585. doi: 10.7623/syxb201512013 SHI Yucai, ZHANG Chen, XUE Lei, et al. A new calculation method of wellbore separation factor[J]. Acta Petrolei Sinica, 2015, 36(12): 1580-1585. doi: 10.7623/syxb201512013 [7] BROOKS A G. A new look at wellbore-collision probability[R]. SPE 116155, 2010. [8] BROOKS A G, WILSON H. An improved method for computing wellbore position uncertainty and application to collision and target intersection probability Analysis[R]. SPE 36863, 1996. [9] WOLFF C J M, DE WARDT J P. Borehole position uncertainty analsysis of measuring method and deviation of systematic error model[J]. Journal of Petroleum Technology, 1981, 33(12): 339-350. doi: 10.2118/9223-PA [10] 唐宁, 熊祖根, 王贵刚. 定向随钻测量误差分析及应用[J]. 钻采工艺, 2016, 39(5):22-25. doi: 10.3969/J.ISSN.1006-768X.2016.05.08 TANG Ning, XIONG Zugen, WANG Guigang. Error analysis for directional measurement while drilling and its appication[J]. Drilling & Production Technology, 2016, 39(5): 22-25. doi: 10.3969/J.ISSN.1006-768X.2016.05.08 [11] 刘修善. 井眼轨迹不确定性表征方法[J]. 石油勘探与开发, 2019, 46(2):391-396. doi: 10.11698/PED.2019.02.20 LIU Xiushan. Borehole trajectory uncertainty and its characterization[J]. Petroleum Exploration and Development, 2019, 46(2): 391-396. doi: 10.11698/PED.2019.02.20 [12] 汪荣鑫. 数理统计[M]. 西安: 西安交通大学出版社, 1986. WANG Rongxin. Mathematical statistics[M]. Xi’an: Xi’an Jiaotong University Press, 1986. [13] 国家能源局. SY/T 6396—2014 丛式井平台布置及井眼防碰技术要求[S]. 北京: 石油工业出版. National Energy Administration. SY/T 6396—2014 Slot layout and anti-collision technical requirement of platform cluster well drilling [S]. Beijing: Petroleum Industry Press, 2014. [14] POEDJONO B, LOMBARDO G J, PHILLIPS W. Anti-collision risk management standard for well placement[R]. SPE 121040, 2009. -