2. 中国石化石油工程技术研究院, 北京 100101;
3. 中国石化江汉油田分公司石油工程技术研究院, 湖北武汉 430035;
4. 中国石油大学(华东)非常规油气与新能源研究院, 山东青岛 266580
2. Sinopec Research Institute of Petroleum Engineering, Beijing, 100101, China;
3. Petroleum Engineering Technology Research Institute, Sinopec Jianghan Oilfield Company, Wuhan, Hubei, 430035, China;
4. Research Institute of Unconventional Petroleum and Renewable Energy, China University of Petroleum(Huadong), Qingdao, Shandong, 266580, China
涪陵页岩气田所处地区属山地丘陵地貌, 为了最大限度地减少井场数量、减小单井占地面积及降低地面工程造价, 提高页岩气整体开发效益, 多采用"井工厂"技术进行开发[1-6]。为了满足集中压裂的需求, 水平段井眼方位一般设计为与最大主应力方向垂直或近似垂直, 导致实际井口与入窗点的闭合方位及水平段井眼方位不在同一条直线上, 因而该地区水平井大多为三维水平井[7-10]。相对于常规水平井, 三维水平井的井眼轨道设计更为复杂, 施工难度更大[11-15]。目前, 国内外学者建立了多种类型的三维井眼轨道设计模型, 白家祉等人[16]最早开始进行三维井眼轨道设计方法的研究, 并给出了三维井眼轨道随钻修正设计与计算方法。刘修善等人[17]提出了"空间圆弧段—直线段—空间圆弧段"的三维井眼轨道设计方法, 满足了给定井眼方向修正轨道的设计要求。唐雪平等人[18]建立了"直线段—圆弧段—直线段—圆弧段—直线段"的三维五段式井眼轨道设计模型, 得到了限定目标点井眼方向的各种轨道设计组合形式的统一表达式。黄根炉等人[19]提出了逐点寻优的设计方法, 完善了限定目标点井眼方向的待钻井眼轨道设计方法。鲁港[20]通过数值算法将三维井眼轨道的多元非线性方程组转化为一元高次多项式并求得拟解析解, 进一步完善了三维井眼轨道设计方法的求解问题。
三维井眼轨道设计理论经过不断发展已较为完善, 能满足常规三维井眼轨道设计的要求。目前, 涪陵页岩气田的三维水平井多采用常规三维五段式轨道设计方法, 即以"直线段—圆弧段—直线段—圆弧段—直线段"进行设计, 但在现场施工过程中为了降低中靶难度, 大多将既增斜又扭方位的第二个圆弧段分为2段进行施工, 即先扭方位对准靶点方位, 再增斜中靶, 但在设计阶段却并未形成对应的三维轨道设计方法, 导致实钻过程中并无切实的设计轨道进行参照, 使钻井风险增大。因此, 为了保证涪陵页岩气田三维水平井安全高效钻进, 提出了符合该地区实钻特点的轨道设计方法, 即三维轨道采用"直井段—增斜段—稳斜段—稳斜扭方位段—稳斜段—增斜段—水平段"七段式进行设计, 并以轨道总长度最短为目标对轨道进行优化设计。
1 三维水平井轨道设计模型的建立 1.1 设计思路根据涪陵页岩气田三维水平井先扭方位对准靶点方位再增斜中靶的特点, 三维水平井轨道采用"直井段—增斜段—稳斜段—稳斜扭方位段—稳斜段—增斜段—水平段"七段式进行设计。该设计方法可以降低中靶难度, 并且预留出稳斜段, 更方便调整井眼轨迹。
具体设计思路为:将常规三维五段式井眼轨道中既增斜又扭方位的第二个圆弧井段替换为稳斜扭方位井段与铅垂面上的增斜井段, 将三维井眼轨道分为水平投影与垂直剖面进行设计, 按照各段坐标增量与靶点坐标相等的原则, 建立完整的三维轨道设计模型[21]。须知条件为造斜点坐标、靶点坐标、造斜率及入靶井斜角与方位角。三维水平井七段式井眼轨道如图 1所示。图 1中:OA为直井段, A点为造斜点; AB为第一增斜段, 采用铅垂面圆弧曲线法进行设计; BC为稳斜段; CD为稳斜扭方位段, 采用圆柱螺线法进行设计, 此段保持井斜角不变全力进行扭方位作业; DE为第二稳斜段; EF为第二增斜段, 采用铅垂面圆弧曲线法进行设计, 此时方位角已调整好, 只需调整井斜角即可顺利中靶; FT为水平段, T点为靶点。
|
| 图 1 三维水平井七段式轨道示意 Fig.1 Schematic for seven-section trajectory of 3D horizontal wells |
三维水平井七段式轨道水平投影如图 2所示。图 2中, N表示方向朝北; 增斜段AB, EF与直线段BC, DE和FT在水平面上的投影均为直线, 而稳斜扭方位段CD采用圆柱螺线法进行设计, 因此该段在水平面上的投影是一条圆弧曲线。
|
| 图 2 三维水平井七段式轨道水平投影 Fig.2 Horizontal projection for seven-section trajectory of 3D horizontal wells |
稳斜扭方位段CD为圆柱螺线且井斜角保持不变, 因此CD段的曲率半径恒定。该井段水平投影的曲率半径r与方位角变化量γ2为:
|
(1) |
|
(2) |
式中:r为稳斜扭方位段CD水平投影的曲率半径, m; γ2为稳斜扭方位段CD的方位角变化量, (°); R2为稳斜扭方位段CD的曲率半径, m; αC为C点处的井斜角, (°); φC为C点处的方位角, (°); φD为D点处的方位角, (°)。
稳斜扭方位段CD水平投影的圆弧长度(水平投影长度)为:
|
(3) |
式中:SCD为稳斜扭方位段CD的水平投影长度, m。
根据各段坐标增量与靶点坐标相等的原则, 可得:
|
(4) |
式中:NA为A点的南北坐标, m; EA为A点的东西坐标, m; NT为T点的南北坐标, m; ET为T点的东西坐标, m; SAG为水平投影直线AG的长度, m; SGT为水平投影直线GT的长度, m; φB为B点处的方位角, (°); φT为T点处的方位角, (°)。
求解式(4)可得SAG。
水平投影AC段的长度为:
|
(5) |
式中:SAC为水平投影AC段的长度, m。
1.2.2 垂直剖面设计三维水平井七段式轨道的垂直剖面如图 3所示。图 3中, S表示方向朝南, H表示垂深; 稳斜扭方位段CD在垂直剖面上为一条直线, 而增斜段AB与EF在垂直剖面上为圆弧曲线。
|
| 图 3 三维水平井七段式轨道的垂直剖面 Fig.3 Vertical profile for seven-section trajectory of 3D horizontal wells |
由垂直剖面可知:
|
(6) |
|
(7) |
式中:γ1为增斜段AB的井斜角变化量, (°); γ3为增斜段EF的井斜角变化量, (°); αA为A点处的井斜角, (°); αB为B点处的井斜角, (°); αE为E点处的井斜角, (°); αF为F点处的井斜角, (°)。
已知造斜点A的坐标为(HA, NA, EA), 可知:
|
(8) |
式中:LOA为直井段OA的轨道长度, m; HA为造斜点A的垂深, m。
B点的坐标为:
|
(9) |
式中:HB为B点的垂深, m; NB为B点的南北坐标, m; EB为B点的东西坐标, m; φA为A点处的方位角, (°); R1为增斜段AB的曲率半径, m。
增斜段AB的轨道长度为:
|
(10) |
式中:LAB为增斜段AB的轨道长度, m。
增斜段AB的水平投影长度为:
|
(11) |
式中:SAB为增斜段AB的水平投影长度, m。
稳斜段BC的水平投影长度为:
|
(12) |
式中:SBC为稳斜段BC的水平投影长度, m。
稳斜段BC的轨道长度为:
|
(13) |
式中:LBC为稳斜段BC的轨道长度, m。
C点的坐标为:
|
(14) |
式中:HC为C点的垂深, m; NC为C点的南北坐标, m; EC为C点的东西坐标, m。
稳斜扭方位段CD的轨道长度为:
|
(15) |
式中:LCD为稳斜扭方位段CD的轨道长度, m。
D点的坐标为:
|
(16) |
式中:HD为D点的垂深, m; ND为D点的南北坐标, m; ED为D点的东西坐标, m; αD为D点处的井斜角, (°)。
根据各段坐标增量与靶点坐标相等的原则, 可得:
|
(17) |
式中:HT为靶点T的垂深, m; LDE为稳斜段DE的轨道长度, m; LFT为水平段FT的轨道长度, m; αT为T点处的井斜角, (°); R3为增斜段EF的曲率半径, m。
求解式(17)可得LDE与LFT。
E点的坐标为:
|
(18) |
式中:HE为E点的垂浑, m; NE为E点的南北坐标, m; EE为E点的东西坐标, m。
增斜段EF的轨道长度为:
|
(19) |
式中:LEF为增斜段EF的轨道长度, m。
F点的坐标为:
|
(20) |
式中:HF为F点的垂深, m; NF为F点的南北坐标, m; EF为F点的东西坐标, m; φE为E点处的方位角, (°); φF为F点处的方位角, (°)。
以上各式组成了涪陵页岩气田三维水平井井眼轨道设计模型, 利用该模型可以完成所有井段轨道长度及各节点参数(井斜角、方位角、井深、垂深、南北坐标和东西坐标)的设计, 将所有分段长度相加即得到轨道总长度:
|
(21) |
式中:L为三维水平井轨道的总长度, m。
2 轨道优选方法 2.1 优选目标求解上述三维水平井轨道设计模型时需要预设井斜角αB与造斜方位角φB, 当αB与φB取不同值, 会得到不同形状的设计轨道。因此, 可将αB与φB所有可能的组合均代入上述轨道设计模型中依次进行设计, 并以轨道总长度最短为目标对设计轨道进行优选。
井斜角αB的取值范围为0~90°, 方位角φB的取值范围为0~360°。为简化计算, 井斜角αB的初始值取αB1, 取值间隔为kα; 方位角φB的初始值取φB1, 取值间隔为kφ, 则αB和φB的取值集合为:
|
(22) |
|
(23) |
其中, 
1) 扫描αB与φB所有可能的组合, 将每个组合分别代入模型进行计算;
2) 代入αB与φB后井眼轨道即被确定, 即可求得CD段和AC段的水平投影长度SCD与SAC;
3) 在垂直剖面上确定各井段轨道的长度LOA, LAB, LBC, LCD, LDE, LEF和LFT及各节点参数(井深、井斜角、方位角、垂深、南北坐标和东西坐标);
4) 保存满足LBC≥0、LDE≥0和LFT≥0的设计轨道;
5) 计算满足要求轨道的总长度L, 选择轨道总长度最短的轨道作为最终的设计轨道。
3 实例计算分析涪陵页岩气田已完钻井的井史资料发现, 该地区有大量三维水平井按照三维七段式轨道进行实际钻进, 但依然按"直线段—圆弧段—直线段—圆弧段—直线段"的常规三维五段式轨道进行设计, 如焦页14-3HF井。现以该井的设计数据为基础, 采用上文的设计方法重新进行设计。该井设计造斜点位于垂深1 300.00 m处, 第一造斜率和第二造斜率均为4.8°/30m, 入靶方位角为0°, 入靶井斜角为87.42°。该井的基础数据见表 1。
| 关键点 | 纵坐标/m | 横坐标/m | 垂深/m |
| 井口O | 3 294 622.50 | 18 749 467.50 | |
| 第一靶点T1 | 3 294 973.00 | 18 750 150.00 | 2 449.00 |
| 第二靶点T2 | 3 295 683.00 | 18 750 150.00 | 2 481.00 |
| 第三靶点T3 | 3 296 473.00 | 18 750 150.00 | 2 546.00 |
采用上文设计方法设计的轨道各节点参数见表 2, 设计轨道的水平投影及垂直剖面见图 4和图 5。其中, OA段为直井段, A点为造斜点; AB段为第一增斜段, 造斜方位角保持112.00°不变, 增斜至44.00°; BC段为第一稳斜段; CD段为稳斜扭方位段, 井斜角保持44.00°不变, 方位角从112.00°扭到0°, 此时方位角恰好为3个靶点的连线方位; DE段为第二稳斜段; EF段为第二增斜段, 方位角保持不变, 井斜角增至87.42°, 此井斜角恰好为靶点T1和T2连线的倾角, 即按照此井斜角钻进可以贯穿两靶点; FT1为靶前水平段; 最后连接水平段T1T2与T2T3。
| 关键点 | 井深/ m |
井斜角/ (°) |
方位角/ (°) |
垂深/ m |
南北坐 标/m |
东西坐 标/m |
| O点 | 0 | 0 | 0 | 0 | ||
| A点 | 1 300.00 | 0 | 112.00 | 1 300.00 | 0 | 0 |
| B点 | 1 575.00 | 44.00 | 112.00 | 1 548.76 | -37.65 | 93.19 |
| C点 | 1 959.00 | 44.00 | 112.00 | 1 825.03 | -137.59 | 340.56 |
| D点 | 2 659.07 | 44.00 | 0 | 2 328.57 | 93.05 | 682.50 |
| E点 | 2 674.66 | 44.00 | 0 | 2 339.79 | 103.88 | 682.50 |
| F点 | 2 946.03 | 87.42 | 0 | 2 448.76 | 345.35 | 682.50 |
| T1点 | 2 951.26 | 87.42 | 0 | 2 449.00 | 350.57 | 682.50 |
| T2点 | 3 662.14 | 87.42 | 0 | 2 481.00 | 1 060.74 | 682.50 |
| T3点 | 4 454.81 | 85.30 | 0 | 2 546.00 | 1 850.74 | 682.50 |
|
| 图 4 焦页14-3HF井新设计轨道的水平投影 Fig.4 Horizontal projection for the newly designed trajectory of Well Jiaoye14-3HF |
|
| 图 5 焦页14-3HF井新设计轨道垂直剖面 Fig.5 Vertical profile for the newly designed trajectory of Well Jiaoye 14-3HF |
焦页14-3HF井的原始设计轨道分段参数与实钻轨迹分段参数见表 3和表 4。
| 关键点 | 井深/ m |
井斜角/ (°) |
方位角/ (°) |
垂深/ m |
南北坐 标/m |
东西坐 标/m |
| 井口 | 0 | 0 | 0 | 0 | ||
| 造斜点 | 1 300.00 | 0 | 91.64 | 1 650.00 | 0 | 0 |
| 1 510.61 | 33.70 | 91.64 | 1 498.68 | -1.73 | 60.14 | |
| 2 282.35 | 33.70 | 91.64 | 2 140.75 | -14.01 | 488.14 | |
| T1点 | 2 857.13 | 87.42 | 0 | 2 449.00 | 350.50 | 682.50 |
| T2点 | 3 567.85 | 87.42 | 0 | 2 481.00 | 1 060.50 | 682.50 |
| T3点 | 4 360.53 | 85.26 | 0 | 2 546.00 | 1 850.50 | 682.50 |
| 关键点 | 井深/ m |
井斜角/ (°) |
方位角/ (°) |
垂深/ m |
南北坐 标/m |
东西坐 标/m |
| 井口 | 0 | 0 | 0 | 0 | ||
| 造斜点 | 1 004.97 | 3.60 | 111.30 | 1 004.89 | -0.45 | 7.79 |
| 1 157.52 | 25.00 | 111.00 | 1 153.58 | -15.22 | 34.85 | |
| 1 892.59 | 40.10 | 111.80 | 1 778.22 | -168.16 | 383.47 | |
| 2 516.04 | 44.20 | 0.60 | 2 251.47 | -2.92 | 678.73 | |
| T1点 | 3 020.00 | 83.20 | 0.00 | 2 462.61 | 439.83 | 677.72 |
| T3点 | 4 520.00 | 85.00 | 0.10 | 2 598.00 | 1 932.85 | 687.79 |
由表 3和表 4可知, 原始设计轨道为常规的"直井段—增斜段—稳斜段—增斜扭方位段—水平段"三维五段式轨道, 而实钻轨迹可以看作"直井段—增斜段—长稳斜段—稳斜扭方位段—稳斜段—增斜段—水平段"三维七段式轨道。对比表 2和表 4可知, 实钻轨迹与采用上文设计方法重新设计的轨道类型一致, 且实钻轨迹中的稳斜角(40.10°)和造斜方位角(111.30°)均与新设计轨道的稳斜角(44.00°)和造斜方位角(112.00°)比较贴近, 说明了上文的三维水平井井眼轨道优化设计方法适用于涪陵页岩气田三维水平井的设计与现场施工。
4 结论1) 基于涪陵页岩气田三维水平井实钻的特点, 提出了"直井段—增斜段—稳斜段—稳斜扭方位段—稳斜段—增斜段—水平段"七段式轨道的设计方法。该设计方法的特点是将常规三维五段式轨道中既增斜又扭方位的第二个圆弧井段替换为稳斜扭方位井段与铅垂面上的增斜井段, 以避免在中靶阶段同时进行增斜与扭方位作业, 降低中靶难度。
2) 依据井眼轨道设计的基础理论, 增斜段与稳斜扭方位段分别采用圆弧线法与圆柱螺线法进行设计, 并按照各段坐标增量与靶点坐标相等的原则建立了三维水平井七段式井眼轨道设计模型。
3) 实例计算表明, 预设稳斜角和造斜方位角求解三维水平井七段式井眼轨道设计模型进行轨道设计, 以轨道总长度最短为目标优选轨道的设计方法更适用于涪陵页岩气田三维水平井的设计与现场施工。
| [1] |
牛新明.
涪陵页岩气田钻井技术难点及对策[J]. 石油钻探技术, 2014, 42(4): 1–6.
NIU Xinming. Drilling technology challenges and resolutions in Fuling Shale Gas Field[J]. Petroleum Drilling Techniques, 2014, 42(4): 1–6. |
| [2] |
葛洪魁, 王小琼, 张义.
大幅度降低页岩气开发成本的技术途径[J]. 石油钻探技术, 2013, 41(6): 1–5.
GE Hongkui, WANG Xiaoqiong, ZHANG Yi. A technical approach to reduce shale gas development cost[J]. Petroleum Drilling Techniques, 2013, 41(6): 1–5. |
| [3] |
张金成, 孙连忠, 王甲昌, 等.
"井工厂"技术在我国非常规油气开发中的应用[J]. 石油钻探技术, 2014, 42(1): 20–25.
ZHANG Jincheng, SUN Lianzhong, WANG Jiachang, et al. Application of multi-well pad in unconventional oil and gas development in China[J]. Petroleum Drilling Techniques, 2014, 42(1): 20–25. |
| [4] |
陈平, 刘阳, 马天寿.
页岩气"井工厂"钻井技术现状及展望[J]. 石油钻探技术, 2014, 42(3): 1–7.
CHEN Ping, LIU Yang, MA Tianshou. Status and prospect of multi-well pad drilling technology in shale gas[J]. Petroleum Drilling Techniques, 2014, 42(3): 1–7. |
| [5] | HUMMES O, BOND P, SYMONS W, et al. Using advanced drilling technology to enable well factory concept in the Marcellus shale[R]. SPE 151466, 2012. |
| [6] |
陈新龙, 徐军, 高迅, 等.
涪陵页岩气开发"井工厂"钻机现状及发展[J]. 石油机械, 2015, 43(9): 32–36.
CHEN Xinlong, XU Jun, GAO Xun, et al. Status and prospect of multi-well pad drilling rig for shale gas development in Fuling[J]. China Petroleum Machinery, 2015, 43(9): 32–36. |
| [7] |
周贤海.
涪陵焦石坝区块页岩气水平井钻井完井技术[J]. 石油钻探技术, 2013, 41(5): 26–30.
ZHOU Xianhai. Drilling & completion techniques used in shale gas horizontal wells in Jiaoshiba Block of Fuling Area[J]. Petroleum Drilling Techniques, 2013, 41(5): 26–30. |
| [8] |
刘匡晓, 王庆军, 兰凯, 等.
涪陵页岩气田三维水平井大井眼导向钻井技术[J]. 石油钻探技术, 2016, 44(5): 16–21.
LIU Kuangxiao, WANG Qingjun, LAN Kai, et al. Large diameter hole steering drilling technology for three-dimensional horizontal well in the Fuling Shale Gas Field[J]. Petroleum Drilling Techniques, 2016, 44(5): 16–21. |
| [9] |
艾军, 张金成, 臧艳彬, 等.
涪陵页岩气田钻井关键技术[J]. 石油钻探技术, 2014, 42(5): 9–15.
AI Jun, ZHANG Jincheng, ZANG Yanbin, et al. The key drilling technologies in Fuling Shale Gas Field[J]. Petroleum Drilling Techniques, 2014, 42(5): 9–15. |
| [10] |
宋争.
涪陵江东与平桥区块页岩气水平井井眼轨迹控制技术[J]. 石油钻探技术, 2017, 45(6): 14–18.
SONG Zheng. Wellbore trajectory control techniques for horizontal well in the Jiangdong and Pingqiao Blocks of the Fuling Shale Gas Field[J]. Petroleum Drilling Techniques, 2017, 45(6): 14–18. |
| [11] |
沈国兵, 刘明国, 晁文学, 等.
涪陵页岩气田三维水平井井眼轨迹控制技术[J]. 石油钻探技术, 2016, 44(2): 10–15.
SHEN Guobing, LIU Mingguo, CHAO Wenxue, et al. 3D trajectory control technology for horizontal wells in the Fuling Shale Gas Field[J]. Petroleum Drilling Techniques, 2016, 44(2): 10–15. |
| [12] |
陶丽杰.
添加虚拟靶区优化三维水平井井眼轨道[J]. 断块油气田, 2015, 22(2): 267–269.
TAO Lijie. Optimization on 3D trajectory of horizontal well by adding theoretical target[J]. Fault-Block Oil & Gas Field, 2015, 22(2): 267–269. |
| [13] |
张凯.
增设虚拟靶点控制水平井井眼轨道设计技术[J]. 石油钻采工艺, 2015, 37(2): 5–7.
ZHANG Kai. Add virtual target to control horizontal well borehole trajectory design techniques[J]. Oil Drilling & Production Technology, 2015, 37(2): 5–7. |
| [14] |
乐守群, 王进杰, 苏前荣, 等.
涪陵页岩气田水平井井身结构优化设计[J]. 石油钻探技术, 2017, 45(1): 17–20.
YUE Shouqun, WANG Jinjie, SU Qianrong, et al. The optimization of casing programs for horizontal wells in the Fuling Shale Gas Field[J]. Petroleum Drilling Techniques, 2017, 45(1): 17–20. |
| [15] |
孙腾飞, 高德利, 杜刚, 等.
目标垂深不确定条件下的水平井轨道设计[J]. 断块油气田, 2012, 19(4): 526–528.
SUN Tengfei, GAO Deli, DU Gang, et al. Design of horizontal well trajectory under the uncertainty condition of vertical depth of target[J]. Fault-Block Oil & Gas Field, 2012, 19(4): 526–528. |
| [16] |
白家祉, 苏义脑.
定向钻井过程中的三维井身随钻修正设计与计算[J]. 石油钻采工艺, 1991, 13(6): 1–4.
BAI Jiazhi, SU Yinao. Correction design and calculation of 3D wells while directional drilling[J]. Oil Drilling & Production Technology, 1991, 13(6): 1–4. |
| [17] |
刘修善, 石在虹.
给定井眼方向的修正轨道设计方法[J]. 石油学报, 2002, 23(2): 72–76.
LIU Xiushan, SHI Zaihong. A new method of path-correction planning with the desired direction[J]. Acta Petrolei Sinica, 2002, 23(2): 72–76. DOI:10.7623/syxb200202016 |
| [18] |
唐雪平, 苏义脑, 陈祖锡.
三维井眼轨道设计模型及其精确解[J]. 石油学报, 2003, 24(4): 90–93, 98.
TANG Xueping, SU Yinao, CHEN Zuxi. Three-dimensional well-path planning model and its exact solution[J]. Acta Petrolei Sinica, 2003, 24(4): 90–93, 98. DOI:10.7623/syxb200304021 |
| [19] |
黄根炉, 赵金海, 赵金洲.
复杂多目标井靶区轨道设计方法研究[J]. 天然气工业, 2006, 26(10): 69–71.
HUANG Genlu, ZHAO Jinhai, ZHAO Jinzhou. Study on trajectory design in target area for complex multi-target directional well[J]. Natural Gas Industry, 2006, 26(10): 69–71. DOI:10.3321/j.issn:1000-0976.2006.10.021 |
| [20] |
鲁港.
圆弧型井眼轨道设计问题的拟解析解理论[J]. 石油钻探技术, 2014, 42(1): 26–32.
LU Gang. Quasi-analytic solution theory for arc type well trajectory design[J]. Petroleum Drilling Techniques, 2014, 42(1): 26–32. |
| [21] |
韩志勇.
定向钻井设计与计算[M]. 东营: 中国石油大学出版社, 2007: 234-245.
HAN Zhiyong. Design and calculation of directional drilling[M]. Dongying: China University of Petroleum Press, 2007: 234-245. |