密集井网直井段井眼轨道交碰风险计算新方法

李琪 刘毅 王六鹏 高云文 张燕娜 张明

李琪,刘毅,王六鹏,高云文,张燕娜,张明. 密集井网直井段井眼轨道交碰风险计算新方法[J]. 石油钻采工艺,2021,43(1):29-33 doi:  10.13639/j.odpt.2021.01.005
引用本文: 李琪,刘毅,王六鹏,高云文,张燕娜,张明. 密集井网直井段井眼轨道交碰风险计算新方法[J]. 石油钻采工艺,2021,43(1):29-33 doi:  10.13639/j.odpt.2021.01.005
LI Qi, LIU Yi, WANG Liupeng, GAO Yunwen, ZHANG Yanna, ZHANG Ming. A new method for calculating the wellbore collision risk of vertical sections of dense cluster wells[J]. Oil Drilling & Production Technology, 2021, 43(1): 29-33 doi:  10.13639/j.odpt.2021.01.005
Citation: LI Qi, LIU Yi, WANG Liupeng, GAO Yunwen, ZHANG Yanna, ZHANG Ming. A new method for calculating the wellbore collision risk of vertical sections of dense cluster wells[J]. Oil Drilling & Production Technology, 2021, 43(1): 29-33 doi:  10.13639/j.odpt.2021.01.005

密集井网直井段井眼轨道交碰风险计算新方法

doi: 10.13639/j.odpt.2021.01.005
基金项目: 国家自然科学基金“基于多源信息和智能计算的钻井异常自适应预警方法研究”(编号:51574194)
详细信息
    作者简介:

    李琪(1963-),2002年获西南石油大学石油与天然气工程专业博士学位,主要从事石油钻井信息技术、旋转导向钻井等研究工作,教授,博士生导师。通讯地址:(710065)陕西省西安市电子二路东段18号西安石油大学。Email:liqi@xsyu.edu.cn

    通讯作者:

    刘毅(1995-),主要从事井眼轨道优化设计与井眼防碰风险分析研究工作。通讯地址:(719000)陕西省榆林市榆阳区高新技术产业园长庆油田分公司第二采气厂。E-mail:2724581053@qq.com

  • 中图分类号: TE24

A new method for calculating the wellbore collision risk of vertical sections of dense cluster wells

  • 摘要: 井眼防碰是密集丛式井网定向井施工面临的一个主要难题,为了有效减小密集井网加密井直井段井眼交碰风险,综合考虑名义间距、基准井和比较井的井径之和以及邻井数等因素,在分析Brooks运用概率评价井眼交碰方法的基础上,提出了比较全面的密集丛式井网井眼交碰风险计算新方法。该方法以基准井井眼中心轴线为圆心、基准井与比较井间名义间距为半径作圆,通过基准井、比较井井径之和与此圆周长的比值计算偏移概率,采用井眼测点位置坐标的正态分布概率密度函数沿井深方向的定积分计算交碰概率,基准井与比较井的实际交碰概率为偏移概率与交碰概率的乘积,通过基准井与各比较井实际交碰概率求和确定综合交碰风险,并将综合交碰风险按照行业标准划分为3个风险等级,得到以综合交碰风险为评价指标的风险评价方法。该方法可定量评价一口基准井与多口比较井之间的综合交碰风险,对加密井钻井防碰理论研究及指导生产实践具有重要意义。
  • 图  1  多口比较井工况下的交碰风险分析

    Figure  1.  Collision risk analysis in the working condition of multiple comparison wells

    图  2  基准井与比较井相对位置示意图

    Figure  2.  Sketch of the relative position of comparison well to reference well

    表  1  密集井网井眼防碰风险分析量化表

    Table  1.   Quantitative wellbore collision risk analysis of dense cluster wells

    RWCCP交碰等级交碰风险预防措施
    10−3RWCCP<3一级重大风险停钻直到风险消除
    10−4RWCCP<10−3二级较小风险关闭邻井制定防碰措施
    0<RWCCP<10−4三级警戒时刻监测邻井情况
    下载: 导出CSV

    表  2  调整加密井部分测斜数据

    Table  2.   Partial directional data of adjustment infill wells

    测点垂深/m井斜角/(°)方位角/(°)N/mE/m
    10.000.000.000.000.00
    239.931.00230.58−0.22−0.27
    368.490.91223.84−0.54−0.62
    497.140.67212.47−0.85−0.87
    5125.730.61240.17−1.07−1.09
    6154.320.50259.84−1.16−1.34
    7182.791.0046.19−1.01−1.29
    8211.340.8828.97−0.65−1.00
    9239.921.2544.76−0.24−0.67
    10268.511.2549.320.19−0.22
    11285.970.83108.460.270.05
    12313.970.98102.700.160.47
    13326.331.0598.130.120.69
    14354.911.76145.44−0.281.20
    15383.482.35181.67−1.231.43
    16411.913.96181.14−2.791.39
    17440.356.25175.56−5.331.49
    18468.698.65173.30−9.021.86
    19496.7812.72175.49−14.302.35
    20524.5714.20180.27−20.942.58
    下载: 导出CSV

    表  3  井深470 m处邻井井身轨迹基本参数

    Table  3.   Basic wellbore trajectory parameters of the neighboring wells at the well depth of 470 m

    井号垂深/m井斜角/(°)方位角/(°)N/mE/m
    #147210.38216.53−7.222.12
    #246911.73157.26−1.213.14
    下载: 导出CSV
  • [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.
  • [1] 黄忠桥, 罗雄, 魏红燕, 吕瑞典.  三维井眼抽油杆扶正器间距配置计算方法 . 石油钻采工艺, 2010, 32(3): 60-63.
    [2] 宋广健, 严建奇, 王丽珍, 王春耘, 卢印生.  油气上窜速度计算方法的改进与应用 . 石油钻采工艺, 2010, 32(5): 17-19.
    [3] 张建兵, 吕祥鸿, 杨勇, 刘歆.  油井管强度计算方法新ISO标准讨论 . 石油钻采工艺, 2009, 31(2): 68-73.
    [4] 程木林, 杨万有, 赵骊川, 谢建华, 王继东.  加密井区油井压裂增产技术 . 石油钻采工艺, 2006, 28(2): 69-71. doi: 10.3969/j.issn.1000-7393.2006.02.021
    [5] 张振华, 鄢捷年.  储层损害对大斜度定向井产能影响的计算方法研究 . 石油钻采工艺, 2001, 23(6): 40-43,83. doi: 10.3969/j.issn.1000-7393.2001.06.012
    [6] 曲占庆, 董长银, 张琪.  影响高压注气管柱变形的主要因素及计算方法 . 石油钻采工艺, 2000, 22(1): 53-55,72. doi: 10.3969/j.issn.1000-7393.2000.01.014
    [7] 周生田, 张琪.  水平井筒压降计算方法 . 石油钻采工艺, 1997, 19(1): 53-59. doi: 10.3969/j.issn.1000-7393.1997.01.013
    [8] 徐华义, 余志清.  适应小井眼的反循环井控方法及其计算 . 石油钻采工艺, 1997, 19(6): 10-15. doi: 10.3969/j.issn.1000-7393.1997.06.003
    [9] 张义堂, 闫亚茹.  垂直裂缝井三线性流模型与Cinco-Ley模型 . 石油钻采工艺, 1996, 18(6): 64-68,79. doi: 10.3969/j.issn.1000-7393.1996.06.012
    [10] 窦宏恩.  预测水平井产能的一种新方法 . 石油钻采工艺, 1996, 18(1): 76-81. doi: 10.3969/j.issn.1000-7393.1996.01.018
    [11] 徐玉兵, 崔孝秉.  注蒸汽井温度场分析计算设计原理 . 石油钻采工艺, 1996, 18(4): 53-61. doi: 10.3969/j.issn.1000-7393.1996.04.014
    [12] 龚伟安.  斜井与水平井套管扶正器计算的三弯矩方程解法 . 石油钻采工艺, 1996, 18(4): 15-24. doi: 10.3969/j.issn.1000-7393.1996.04.004
    [13] 赵洪激.  水平井钻柱摩擦力的计算 . 石油钻采工艺, 1995, 17(3): 6-1. doi: 10.3969/j.issn.1000-7393.1995.03.002
    [14] 张先普, 练章华, 王钢, 何明烛, 卫玉祥.  荆丘油田套管变形机理及防治方法 . 石油钻采工艺, 1994, 16(2): 1-9. doi: 10.3969/j.issn.1000-7393.1994.02.001
    [15] 刘从彪, 刘翔鹗.  利用油田开发资料计算就地三向主应力 . 石油钻采工艺, 1994, 16(5): 45-49. doi: 10.3969/j.issn.1000-7393.1994.05.012
    [16] 胡明君.  水平井造斜轨迹的设计通式 . 石油钻采工艺, 1994, 16(5): 29-32,37. doi: 10.3969/j.issn.1000-7393.1994.05.006
    [17] 周效全.  物质溶解度参数的计算方法 . 石油钻采工艺, 1991, 13(3): 63-70. doi: 10.3969/j.issn.1000-7393.1991.03.015
    [18] 苏义脑.  求定向井井底真实钻压值的理论分析和初步计算 . 石油钻采工艺, 1991, 13(2): 31-36,42. doi: 10.3969/j.issn.1000-7393.1991.02.006
    [19] 周开吉, 杨金华, 钟兵, 唐林.  井内波动压力 . 石油钻采工艺, 1990, 12(2): 1-10. doi: 10.3969/j.issn.1000-7393.1990.02.001
    [20] 李福军, 蒋汉青, 齐振林.  抽油机井系统效率计算方法的改进 . 石油钻采工艺, 1989, 11(6): 75-78,90. doi: 10.3969/j.issn.1000-7393.1989.06.012
  • 加载中
图(2) / 表 (3)
计量
  • 文章访问数:  7
  • HTML全文浏览量:  132
  • PDF下载量:  43
  • 被引次数: 0
出版历程
  • 修回日期:  2020-12-19
  • 网络出版日期:  2021-05-24
  • 刊出日期:  2021-01-20

密集井网直井段井眼轨道交碰风险计算新方法

doi: 10.13639/j.odpt.2021.01.005
    基金项目:  国家自然科学基金“基于多源信息和智能计算的钻井异常自适应预警方法研究”(编号:51574194)
    作者简介:

    李琪(1963-),2002年获西南石油大学石油与天然气工程专业博士学位,主要从事石油钻井信息技术、旋转导向钻井等研究工作,教授,博士生导师。通讯地址:(710065)陕西省西安市电子二路东段18号西安石油大学。Email:liqi@xsyu.edu.cn

    通讯作者: 刘毅(1995-),主要从事井眼轨道优化设计与井眼防碰风险分析研究工作。通讯地址:(719000)陕西省榆林市榆阳区高新技术产业园长庆油田分公司第二采气厂。E-mail:2724581053@qq.com
  • 中图分类号: TE24

摘要: 井眼防碰是密集丛式井网定向井施工面临的一个主要难题,为了有效减小密集井网加密井直井段井眼交碰风险,综合考虑名义间距、基准井和比较井的井径之和以及邻井数等因素,在分析Brooks运用概率评价井眼交碰方法的基础上,提出了比较全面的密集丛式井网井眼交碰风险计算新方法。该方法以基准井井眼中心轴线为圆心、基准井与比较井间名义间距为半径作圆,通过基准井、比较井井径之和与此圆周长的比值计算偏移概率,采用井眼测点位置坐标的正态分布概率密度函数沿井深方向的定积分计算交碰概率,基准井与比较井的实际交碰概率为偏移概率与交碰概率的乘积,通过基准井与各比较井实际交碰概率求和确定综合交碰风险,并将综合交碰风险按照行业标准划分为3个风险等级,得到以综合交碰风险为评价指标的风险评价方法。该方法可定量评价一口基准井与多口比较井之间的综合交碰风险,对加密井钻井防碰理论研究及指导生产实践具有重要意义。

English Abstract

李琪,刘毅,王六鹏,高云文,张燕娜,张明. 密集井网直井段井眼轨道交碰风险计算新方法[J]. 石油钻采工艺,2021,43(1):29-33 doi:  10.13639/j.odpt.2021.01.005
引用本文: 李琪,刘毅,王六鹏,高云文,张燕娜,张明. 密集井网直井段井眼轨道交碰风险计算新方法[J]. 石油钻采工艺,2021,43(1):29-33 doi:  10.13639/j.odpt.2021.01.005
LI Qi, LIU Yi, WANG Liupeng, GAO Yunwen, ZHANG Yanna, ZHANG Ming. A new method for calculating the wellbore collision risk of vertical sections of dense cluster wells[J]. Oil Drilling & Production Technology, 2021, 43(1): 29-33 doi:  10.13639/j.odpt.2021.01.005
Citation: LI Qi, LIU Yi, WANG Liupeng, GAO Yunwen, ZHANG Yanna, ZHANG Ming. A new method for calculating the wellbore collision risk of vertical sections of dense cluster wells[J]. Oil Drilling & Production Technology, 2021, 43(1): 29-33 doi:  10.13639/j.odpt.2021.01.005
    • 近年来海洋油气勘探进入高速发展阶段[1],人工岛是实现“海油陆采”的一种高效开发模式。以冀东油田为例,已先后开发高尚堡、柳赞、老爷庙、唐海、南堡5座人工岛。其中南堡1-3D人工岛占地0.13 km2,截止2018年底已完钻井200口以上,但开发中后期采油速度大幅降低,为完成生产任务需在人工岛上进行丛式井组加密调整。目前地面可利用空间和井口十分有限,地下已完钻密集井段安全穿越也愈发困难,新建加密井数量越多相邻井眼之间的距离越小,在钻进时出现井眼交碰的概率就越大,任何非预期的井眼交碰都会给油田造成严重损失,因此防碰风险评价是密集井网加密井设计一个必不可少的环节[2-4]。常用的井眼交碰风险评价方法有井眼分离系数法、井眼交碰概率法以及井眼最近距离法[5-7],这些方法的共同点是评价体系中只有一口比较井,而密集井网条件下在防碰分析范围内的井眼往往存在多口,要规避其中任意井,势必造成与其他井的交碰风险增加。可见,在密集井网条件下采用现有方法进行井眼交碰风险评价不够全面准确。为此,根据密集井网加密井设计的特点,提出了一种可定量评价一口基准井与多口比较井的密集井网加密井井眼交碰综合风险评价新方法。

    • 1996年Brooks等首次提出井眼交碰概率的概念并给出了数学计算方法[8],该方法通过目标点和参考点的协方差矩阵计算井眼交碰概率,可定量评价2个误差椭球的分离程度,但是该方法只考虑到目标点和参考点的位置不确定性矩阵,忽略了两井眼名义间距对交碰概率的影响,评价结果难以满足防碰钻井要求。

      笔者改进了Brooks的交碰概率模型,综合考虑名义间距、基准井和比较井的井径之和以及邻井数等因素,建立了一个比较全面的密集丛式井网井眼交碰风险计算方法。通过几何概率模型计算井眼偏向某邻井方向的概率θi,运用Brooks的井眼交碰概率模型计算相邻井眼的交碰概率Pi,井眼实际交碰风险值为θiPi乘积,则密集井网井眼交碰综合风险值为基准井与各比较井井眼实际交碰风险值的和,和值越小井眼综合交碰风险程度越低,否则越高。图1为密集井网某一深度基准井、比较井相对位置示意图,图中O为基准井,O1O2O3为比较井,d1为基准井井径,d2为比较井井径,S1S2S3分别为基准井到各比较井的名义间距,则基准井向O1方向倾斜的概率等于长度为d1+d2的弧长与半径为S1的圆周周长的比值,记为θ1,运用Brooks交碰概率计算模型可计算得到基准井O与比较井O1的交碰概率P1,同理可以得到θ2θ3P2P3,则图示密集井网的综合交碰风险值为$\displaystyle\sum\nolimits_{i = 1}^3 {{\theta _i}} {P_i}$

      图  1  多口比较井工况下的交碰风险分析

      Figure 1.  Collision risk analysis in the working condition of multiple comparison wells

    • 受地层非均质性、钻具组合、钻头类型等因素影响,钻井过程中会出现方位偏移,偏移方向具有随机性但服从几何概率分布模型,基准井与某一比较井的相对位置如图2所示,图中O为基准井、O1为比较井,基准井、比较井井眼直径分别为d1d2,名义间距为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  基准井与比较井相对位置示意图

      Figure 2.  Sketch of the relative position of comparison well to reference well

      图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]。假设通过井眼轨迹位置不确定性分析方法,得到目标点和参考点的位置不确定矩阵C1C2,其中不确定矩阵C1C2相互独立,则两个矩阵之和表示目标点和参考点位置的相对不确定性,记为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)
    • 通过上述计算方法得到目标井与各比较井的θiPi,则密集井网井眼交碰综合风险值(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−3RWCCP<3一级重大风险停钻直到风险消除
      10−4RWCCP<10−3二级较小风险关闭邻井制定防碰措施
      0<RWCCP<10−4三级警戒时刻监测邻井情况
    • 冀东南堡3号人工岛某设计调整加密井部分测斜数据如表2所示,将该轨迹数据导入Campass采用最近距离扫描,扫描结果显示与邻井的分离系数均大于1,但按照设计进行施工在井深470 m处却发生了非预期的井眼交碰,给油田带来较大的经济损失。

      表 2  调整加密井部分测斜数据

      Table 2.  Partial directional data of adjustment infill wells

      测点垂深/m井斜角/(°)方位角/(°)N/mE/m
      10.000.000.000.000.00
      239.931.00230.58−0.22−0.27
      368.490.91223.84−0.54−0.62
      497.140.67212.47−0.85−0.87
      5125.730.61240.17−1.07−1.09
      6154.320.50259.84−1.16−1.34
      7182.791.0046.19−1.01−1.29
      8211.340.8828.97−0.65−1.00
      9239.921.2544.76−0.24−0.67
      10268.511.2549.320.19−0.22
      11285.970.83108.460.270.05
      12313.970.98102.700.160.47
      13326.331.0598.130.120.69
      14354.911.76145.44−0.281.20
      15383.482.35181.67−1.231.43
      16411.913.96181.14−2.791.39
      17440.356.25175.56−5.331.49
      18468.698.65173.30−9.021.86
      19496.7812.72175.49−14.302.35
      20524.5714.20180.27−20.942.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/mE/m
      #147210.38216.53−7.222.12
      #246911.73157.26−1.213.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)该方法可定量评价一口基准井与多口比较井之间的综合交碰风险,对加密井钻井防碰理论研究及指导生产实践具有重要意义。

参考文献 (14)

目录

    /

    返回文章
    返回