摘要
随着定向钻进技术应用越来越广泛,钻孔轨迹的控制技术变得尤为重要。为了研究钻孔自然弯曲规律,找到顶角和方位角的变化趋势,进而预测待钻钻孔轨迹的变化,使用SPSS软件对安徽省一号矿区的钻孔测斜数据进行多元方程回归统计分析。在考虑孔深、开孔点空间位置信息(x、y、z坐标)对顶角和方位角影响的基础上,建立了四元八次多项式回归方程,并采用均角全距法,将回归拟合的钻孔轴线空间位置与实际钻孔空间位置进行对比。结果表明:将钻孔原始数据进行线性处理后,用四元八次多项式进行回归,能得到较好的效果,若按照钻孔所处地层产状的不同进行分类后进行回归,拟合程度将进一步提高。利用拟合度高的回归方程,能用来预测未开孔地区的钻孔自然弯曲轨迹,对工程施工控制钻孔轨迹也起到了良好的作用,能更好指导工程的施工。
随着各种定向钻进技术的应用越来越广泛,对于定向钻孔而言,保证钻孔轨迹满足设计要求,并准确进入目标区域(靶区)是定向钻探工程中的关键,因此钻孔轨迹的掌握与控制显得更加重
对于钻孔自然弯曲规律分析,常用的方法是统计作图法或相关分析法。其中相关分析法是通过回归方程,以孔深为自变量,以顶角、方位角为因变量,分别得出顶角、方位角与孔深之间的关系,来反映钻孔的自然弯曲趋势。这种方法的好处在于能以较高的准确度反映钻孔空间位置的变化规律,且能随时掌握顶角、方位角两个参数的变化,当变化超过限度时,能及时对钻孔的偏斜进行纠正与预防,避免轨迹偏离严重后难以进行补救措
SPSS软件可用于统计学分析运算、数据挖掘、预测分析和决策支持任务等方
除了上述领域运用到SPSS,在钻探领域SPSS软件也有所应用。葛子明
SPSS软件能很好的对多元方程进行回归分析,可以探究多个变量对顶角、方位角的影响。目前,在钻孔轨迹自然弯曲规律分析中还没有广泛运用SPSS进行回归拟合分析。
一号矿区位于潘集背斜及东部倾伏转折端南翼。潘集背斜为一轴部开阔南缓北陡的不对称背斜。轴迹NWW向,枢纽总体向SEE倾伏,倾伏角5°,沿走向稍有起伏。地层走向自东向西为N30°E至N60°W,倾向SE~SW。背斜北翼倾角较陡,为20°~30°,南翼倾角较缓,为5°~23°,井田内地层倾角由浅入深逐渐变缓(23°~5°)。一号矿区内主要发育两组张扭性断层,倾角均为50°~75°。其次一号矿区内还发育有压扭性断层,该断层属于逆断层,断层的走向与背斜的轴向基本一致或者形成20°~30°的夹角。断层的落差较大,是确定一号矿区边界和采矿区域边界的地质依据。
据钻探揭露及邻区资料,一号矿区内自老至新的地层有奥陶系中下统、石炭系、二叠系、新近系和第四系。含煤的地层主要分布在二叠系山西组、上石盒子组和下石盒子组。
由于泥岩、泥灰岩、灰岩、煤层的岩性差别很大,在钻进过程中,这些地层交替出现,形成地层的软硬互层现象,会使得钻孔发生弯曲。而且矿区断层、破碎地层发育较多。在这些地层处钻进时,会因地层掉块、脱落造成钻具受力不平衡,从而造成钻孔超径,增大了孔壁间隙。这些现象都会导致钻孔轨迹的弯曲。
本文选取了一号矿区内已有的56个钻孔的共849组测斜数据进行回归分析,每组数据包括钻孔钻进的孔深L,钻孔初始点经转化后的横、纵、竖轴相对坐标x、y、z,顶角和方位角。
顶角和方位角随孔深L的关系选择常用的多项式函数拟合,次数越高,其拟合程度越好,本文矿区钻孔的顶角和方位角随孔深的关系选择六次多项式函数来做回归方程,可写为
(1) |
式中:——因变量顶角(方位角),(°);——多项式系数;L——孔深,m。
考虑到钻孔开孔坐标的不同,会决定钻孔所处的地层岩性、地质情况也有所不同,进而影响到顶角和方位角随孔深的变化规律,当钻孔初始点坐标发生改变,该多项式也会随之改变,即不同钻孔的顶角、方位角随孔深的变化趋势不同,这符合实际规律。因此钻孔初始点坐标x、y、z影响着该多项式的系数,
(2) |
式中:x、y、z——分别为钻孔初始点的横、纵、竖坐标,m;——系数。
将
对一号矿区钻孔轨迹数据进行回归,因深度L,孔口坐标x、y、z数据的数量级各不一致,过大数量级的数据如孔口坐标x、y,在高次方程的拟合中会使得高次项的系数很小,SPSS软件会显示不出准确的系数值,进而使得误差偏大,故需要缩小4个变量数据的数量级使其一致。吴鹏
由于钻孔深度L存在0数据,无法对数处理,且对数变换后,各变量的相关性变强,显著变量变少,回归方程的拟合度不高。因此本文采用线性变换,将4个变量数据的数量级调整一致,在不改变原多项式的情况下,将过小的系数也回归出来。
在使用SPSS软件对矿区所有钻孔数据的顶角、方位角进行多元回归分析时,先对深度L、孔口坐标x、y、z的数据进行线性变换,统一缩小数量级,再使用四元八次多项式进行多元回归分析,变量筛选策略选择后退筛选。发现顶角的调整后
地层倾角/(°) | 调整后 | 调整后 |
---|---|---|
0~10 | 0.848 | 0.511 |
10~15 | 0.833 | 0.787 |
15~25 | 0.785 | 0.573 |
25~75 | 0.705 | 0.644 |
以地层倾角<10°时的钻孔数据,因变量为顶角的回归过程为例,自变量采用后退筛选策略,通过SPSS软件进行多元回归分析后能得到以下输出结果。
首先给出了采用后退回归法变量的引入和剔除情况,并且显示引入与剔除的判别标准。引入的标准是将变量全输入方程,直到容差为0,剔除的标准是使用F的概率P>0.1。一共建立了15个回归模型,第一个模型引入了四元八次多项式中的30个显著变量后,使得容差为0,不再继续引入变量,然后在这30个变量中逐步剔除显著性不大的项,最终剔除了14个变量,保留下16个与因变量顶角最显著的相关变量,以第15个回归模型作为最终采用的回归模型。
由第15个模型的模型摘要可知模型15的复相关系数R为0.924,判定系数
模型15的回归方程显著性检验结果见
模型15 | 平方和 | 自由度 | 均方 | F | 显著性 |
---|---|---|---|---|---|
回归 | 2088.938 | 16 | 130.559 | 127.380 | 0.000 |
残差 | 355.660 | 347 | 1.025 | ||
总计 | 2444.597 | 363 |
根据回归方程的系数表,可知回归方程中16个变量的系数,这16个回归系数t统计量相应的伴随概率P均值小于0.05,即有统计学意义,其回归方程可写为
(3) |
在回归分析中,残差需要服从正态分布。带正态曲线的残差分布直方图见

图1 回归方程残差分布直方图
Fig.1 Regression equation residual distribution histogram

图2 标准化残差的概率图
Fig.2 Probability graph of standardized residual
最常用的几种钻孔轴线空间位置的计算方法分别是全角全距法、均角全距法、全角半距法、曲率半径法和最小曲率法等。其中全角全距法虽然计算公式简便,但误差较大,基本很少使用。李静
以一号矿区钻孔wIV-21拟合回归后预测的顶角方位角数据为例,分别用均角全距法、全角半距法、曲率半径法和最小曲率法4种方法进行计算,得到4条拟合预测的钻孔轨迹,
孔深L/ m | 拟合顶角/ (°) | 拟合方位角/(°) | 法1拟合x/m | 法1拟合y/m | 法1拟合z/m | 法2拟合x/m | 法2拟合y/m | 法2拟合z/m |
---|---|---|---|---|---|---|---|---|
0 | 0.01 | 196.63 | 667.13 | 5469.42 | 21.25 | 667.13 | 5469.42 | 21.25 |
50 | 0.31 | 194.45 | 667.00 | 5469.38 | -28.75 | 667.00 | 5469.38 | -28.75 |
100 | 0.62 | 192.42 | 666.60 | 5469.29 | -78.75 | 666.60 | 5469.29 | -78.76 |
150 | 0.93 | 190.16 | 665.94 | 5469.16 | -128.74 | 665.94 | 5469.16 | -128.75 |
200 | 1.30 | 187.42 | 664.98 | 5469.01 | -178.73 | 664.98 | 5469.01 | -178.75 |
250 | 1.69 | 184.12 | 663.68 | 5468.88 | -228.72 | 663.68 | 5468.88 | -228.74 |
300 | 2.12 | 180.25 | 662.02 | 5468.81 | -278.69 | 662.02 | 5468.81 | -278.71 |
350 | 2.56 | 175.89 | 659.98 | 5468.88 | -328.65 | 659.98 | 5468.88 | -328.67 |
400 | 3.05 | 170.79 | 657.55 | 5469.16 | -378.59 | 657.55 | 5469.16 | -378.62 |
450 | 3.55 | 166.78 | 654.73 | 5469.72 | -428.51 | 654.73 | 5469.72 | -428.54 |
500 | 4.10 | 163.22 | 651.51 | 5470.59 | -478.39 | 651.51 | 5470.59 | -478.43 |
550 | 4.64 | 159.43 | 647.90 | 5471.81 | -528.25 | 647.90 | 5471.81 | -528.29 |
600 | 5.27 | 155.46 | 643.90 | 5473.47 | -578.06 | 643.91 | 5473.47 | -578.10 |
650 | 5.96 | 150.67 | 639.54 | 5475.68 | -627.82 | 639.54 | 5475.68 | -627.87 |
700 | 6.72 | 146.65 | 634.82 | 5478.55 | -677.52 | 634.83 | 5478.55 | -677.56 |
750 | 7.55 | 141.98 | 629.78 | 5482.18 | -727.13 | 629.78 | 5482.18 | -727.18 |
800 | 8.48 | 135.87 | 624.52 | 5486.76 | -776.64 | 624.53 | 5486.76 | -776.69 |
850 | 9.45 | 127.42 | 619.34 | 5492.58 | -826.03 | 619.35 | 5492.58 | -826.09 |
900 | 10.55 | 114.99 | 614.85 | 5500.01 | -875.27 | 614.86 | 5499.99 | -875.33 |
950 | 11.79 | 99.90 | 611.94 | 5509.25 | -924.32 | 611.97 | 5509.20 | -924.39 |
975 | 12.57 | 88.55 | 611.55 | 5514.51 | -948.76 | 611.58 | 5514.46 | -948.82 |
978.54 | 12.72 | 86.20 | 611.59 | 5515.28 | -952.21 | 611.61 | 5515.23 | -952.28 |

图3 四种方法计算的预测钻孔轨迹
Fig.3 Predicted borehole trajectory calculated by four methods
在4种不同地层倾角条件下,采用均角全距法,利用钻孔实际测斜数据值和代入回归方程计算的预测值,分别计算出实际钻孔轨迹三维坐标和预测钻孔轨迹三维坐标。钻孔轴线轨迹的三维坐标的计算公式为
(4) |
式中:——分别为上测点的钻孔轴线坐标,m;——分别为下测点的钻孔轴线坐标,m;——分别为上测点的钻孔顶角和方位角,(°);——分别为下测点的钻孔顶角和方位角,(°);——上下测点之间的轴线长度,m。
根据公式计算的三维坐标分别绘制出钻孔的实际轨迹和预测轨迹的三维图。

图4 实际轨迹与预测轨迹三维对比
Fig.4 Three dimensional comparison between actual
trajectory and predicted trajectory
地层倾角为10°~15°的顶角和方位角的回归方程的调整后
地层倾角为0°~10°和地层倾角为15°~25°时,预测的钻孔轨迹形状与实际钻孔轨迹形状就有了明显的不同,误差较大。两种情况下,顶角的回归方程的调整后
地层倾角为0°~10°时,预测的钻孔轨迹和实际钻孔轨迹虽然误差较大,但轨迹形状基本预测准确,如
将4种情况下的预测曲线的起始点统一为三维坐标原点,将4个曲线画入一张三维坐标图进行比较,见

图5 不同倾角条件的预测轨迹对比
Fig.5 Comparison of predicted trajectories under different dip angles
由
造成这种情况可能与地层有关,钻孔基本为竖直孔,当地层倾角小时,地层对钻孔的造斜程度低,钻孔轨迹顶角变化小;而地层倾角增大时,地层对钻孔造斜程度增大,当地层为硬岩层时,此时遇层角较大,钻孔会朝着垂直于层面的方向弯曲,因此,钻孔轨迹顶角变化也增大;而当地层倾角过大时,钻孔的遇层角偏小,小于临界角时可能沿着岩层面下滑,反而使得钻孔轨迹的顶角变化小。
(1)在使用SPSS软件对一号矿区钻孔测斜数据进行多元回归分析时,先将自变量孔深、钻孔起始点坐标x、y、z进行线性变换,使得4个自变量的数量级缩小到同一量级后,回归模型使用四元八次多项式方程进行多元函数拟合,变量的筛选策略选择向后筛选,这样得出方程的拟合结果更好。
(2)对钻孔所处地层的倾角不同,进行数据的分类,用分类后的数据进行多元回归拟合,拟合程度明显增大。
(3)顶角、方位角的回归方程拟合度高,所计算出的预测的钻孔轨迹与实际轨迹的三维坐标误差也小。利用拟合度高的回归方程,能用来预测未开孔地区的钻孔自然弯曲轨迹,更好指导工程的施工。
(4)地层倾角为15°~25°时,造斜程度明显大于其他地层倾角,因此在15°~25°倾角地层中钻进时,要时刻注意顶角和方位角的变化,及时减斜或纠斜。
(5)本文所得出的回归方程只适用于安徽省一号矿区,钻孔数据较少,部分方程拟合度不高,可采用上述方法进一步分析不同矿区的钻孔自然弯曲规律。
参考文献(References)
黄井武,王志强,谭玉峰.煤矿工程钻孔轨迹参数计算方法及绘图[J].建井技术,2012,33(2):19-22. [百度学术]
HUANG Jingwu, WANG Zhiqiang, TAN Yufeng. Calculation method and mapping of coal mine engineering drilling path parameters[J]. Mine Construction Technology, 2012,33(2):19-22. [百度学术]
李静,张金昌,陈晓琳.地质勘探钻孔轨迹计算新模型[J].探矿工程(岩土钻掘工程),2011,38(1):22-24. [百度学术]
LI Jing, ZHANG Jinchang, CHEN Xiaolin. A new calculation model for geological exploration borehole trajectory[J]. Exploration Engineering (Rock &Soil Drilling and Tunneling), 2011,38(1):22-24 [百度学术]
陈俊松,尚亚军,杨建军,等.重庆市秀山县ZK2地热井定向井段钻进工艺[J].钻探工程,2022,49(3):37-43. [百度学术]
CHEN Junsong, SHANG Yajun, YANG Jianjun, et al. Drilling process for the directional section of ZK2 geothermal well in Xiushan county, Chongqing [J]. Drilling Engineering, 2022,49(3):37-43. [百度学术]
康治昌.利用钻孔自然弯曲的几点体会[J].探矿工程,1960(8):15. [百度学术]
KANG Zhichang. Some experiences on using natural bending of boreholes[J]. Exploration Engineering, 1960(8):15. [百度学术]
沈朝旭.地质勘探中钻孔弯曲原因及预防措施[J].现代商贸工业,2013,25(19):180-181. [百度学术]
SHEN Chaoxu. Causes of borehole bending in geological exploration and preventive measures [J]. Modern Business Trade Industry, 2013,25(19):180-181. [百度学术]
齐瑞忱.钻孔轨迹空间位置的预测方法[J].成都地质学院学报,1992(3):113-120. [百度学术]
QI Ruichen. Method for forecasting the space position of the borehole trajectory [J]. Journal of Chengdu College of Geology,1992(3):113-120. [百度学术]
汤凤林,赵荣欣,周欣,等.受控定向孔钻进用自动导航系统的分析研究[J].钻探工程,2022,49(2):67-76. [百度学术]
TANG Fenglin, ZHAO Rongxin, ZHOU xin, et al. An automatic navigation system used for controlled directional drilling [J]. Drilling Engineering, 2022,49(2):67-76. [百度学术]
王宏.利用钻孔自然弯曲样板设计定向孔[J].地质与勘探,1989(7):60-62. [百度学术]
WANG Hong. Directional hole designing by using a template constructed according to the natural bending of the hole [J]. Geology and Exploration, 1989(7):60-62. [百度学术]
王红梅.SPSS软件在课题研究中的应用——以SYB创业培训在高职院校中的应用研究为例[J].中小企业管理与科技,2016(12):124-125. [百度学术]
WANG Hongmei. Application of SPSS software in subject research—Take the application of SYB entrepreneurship training in vocational colleges as an example[J]. Management & Technology of SME, 2016(12):124-125. [百度学术]
薛薇.基于SPSS的数据分析(第3版)[M].北京:中国人民大学出版社,2014. [百度学术]
XUE Wei. Data Analysis Based on SPSS (3rd Edition)[M]. Beijing: China Renmin University Press, 2014. [百度学术]
冯岩松.SPSS22.0统计分析应用教程[M].北京:清华大学出版,2015. [百度学术]
FENG Yansong. Application Course of PSS22.0 Statistical Analysis[M]. Beijing: Tsinghua University Press, 2015. [百度学术]
张文彤,邝春伟.SPSS统计分析基础教程(第2版)[M].北京:高等教育出版社,2011:327-342. [百度学术]
ZHANG Wentong, KUANg Chunwei. Basic Course of SPSS Statistical Analysis (2nd Edition)[M]. Beijing: Higher Education Press, 2011:327-342. [百度学术]
武松,潘发明.SPSS统计分析大全[M].北京:清华大学出版社,2014:217-234. [百度学术]
WU Song, PAN Faming. Complete Statistical Analysis of SPSS[M]. Beijing: Tsinghua University Press, 2014:217-234. [百度学术]
侯海桂.关于统计分析内容分类以及相关SPSS分析方法使用的探讨[J].经济师,2014(5):72-75. [百度学术]
HOU Haigui. Discussion on the classification of statistical analysis contents and the use of relevant SPSS analysis methods [J]. China Economist, 2014(5):72-75. [百度学术]
李明华,崔少萍,罗凤明,等.统计软件SPSS在气象中的应用[J].广东气象,2007,29(1):50-52. [百度学术]
LI Minghua, CUI Shaoping, LUO Fengming, et al. Application of statistical software SPSS in meteorology[J]. Guangdong Meteorology, 2007,29(1):50-52. [百度学术]
罗凤明,邱劲飚,李明华,等.如何使用统计软件SPSS进行回归分析[J].电脑知识与技术,2008,4(2):293-294. [百度学术]
LUO Fengming, QIU Jinbiao, LI Minghua, et al. How to do regression analysis by statistical software SPSS[J]. Computer Knowledge and Technology, 2008,4(2):293-294. [百度学术]
李伟权,陈映华,郭琳晶,等.如何综合使用SPSS和Surfer绘制风向频率随时间变化图[J].电脑知识与技术,2009,5(4):954-955. [百度学术]
LI Weiquan, CHEN Yinghua, GUO Linjing, et al. How to draw chart of variation of wind direction frequency synthetically by SPSS and Surfer[J]. Computer Knowledge and Technology, 2009,5(4):954-955. [百度学术]
SANTOS V, GOMES S, NOGUEIRA M. Sustainable packaging: Does eating organic really make a difference on product-packaging interaction?[J]. Journal of Cleaner Production,2021,304. [百度学术]
陈晓丽,张庆英,郭貔,等.PBL与软件教学在医学统计学实验中的效果评价[J].中国继续医学教育,2021,13(13):20-24. [百度学术]
CHEN Xiaoli, ZHANG Qingying, GUO Pi, et al. Evaluation on the effect of PBL combined with software in experimental teaching of medical statistics[J]. China Continuing Medical Education, 2021,13(13):20-24. [百度学术]
MIRAHMAD H, NOZARI B, RAISSI V, et al. Frequency and parasite load of Toxoplasma gondii in hemodialysis patients based on RE gene by real‑time PCR[J]. Gene Reports, 2021,23. [百度学术]
JIN Long, ZHANG Jun, FU Huiqun, et al. FOXO3a inhibits the EMT and metastasis of breast cancer by regulating TWIST-1 mediated miR-10b/CADM2 axis[J]. Translational Oncology,2021,14(7). [百度学术]
LEMMA A, SALELEW E, DEMILEW D, et al. Alcohol use disorder and associated factors among University of Gondar undergraduate students: A cross‑sectional study[J]. Journal of Substance Abuse Treatment, 2021,129. [百度学术]
葛子明,杨国庆,王浩妮,等.SPSS软件及其在钻具失效的应用[J].数字化用户,2013,19(6):18. [百度学术]
GE Ziming, YANG Guoqing, WANG Haoni, et al. SPSS software and its application in drilling tool failure[J]. Digitization User, 2013,19(6):18. [百度学术]
吴鹏,王幼民,刘司诚,等.面向注塑工艺参数优化的多元高次回归分析模型建立方法[J].塑料工业,2020,48(9):71-75. [百度学术]
WU Peng, WANG Youmin, LIU Sicheng, et al. Multiple higher order regression analysis model establishment method for injection molding process parameter optimization[J]. China Plastics Industry, 2020,48(9):71-75. [百度学术]