Loading [MathJax]/jax/output/SVG/jax.js
  • 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊
  • Scopus数据库收录期刊

露天潜孔钻机随钻测量系统研发与应用研究

侯仕军, 丁伟捷, 田帅康, 梁书锋, 刘殿书, 郭大超

侯仕军, 丁伟捷, 田帅康, 梁书锋, 刘殿书, 郭大超. 露天潜孔钻机随钻测量系统研发与应用研究[J]. 岩土工程学报, 2024, 46(7): 1509-1515. DOI: 10.11779/CJGE20230477
引用本文: 侯仕军, 丁伟捷, 田帅康, 梁书锋, 刘殿书, 郭大超. 露天潜孔钻机随钻测量系统研发与应用研究[J]. 岩土工程学报, 2024, 46(7): 1509-1515. DOI: 10.11779/CJGE20230477
HOU Shijun, DING Weijie, TIAN Shuaikang, LIANG Shufeng, LIU Dianshu, GUO Dachao. Development and application of MWD system for DTH drilling rig[J]. Chinese Journal of Geotechnical Engineering, 2024, 46(7): 1509-1515. DOI: 10.11779/CJGE20230477
Citation: HOU Shijun, DING Weijie, TIAN Shuaikang, LIANG Shufeng, LIU Dianshu, GUO Dachao. Development and application of MWD system for DTH drilling rig[J]. Chinese Journal of Geotechnical Engineering, 2024, 46(7): 1509-1515. DOI: 10.11779/CJGE20230477

露天潜孔钻机随钻测量系统研发与应用研究  English Version

详细信息
    作者简介:

    侯仕军(1995—),男,博士研究生,主要从事岩石破碎和工程爆破方面的研究。E-mail: hsj131723@163.com

    通讯作者:

    梁书锋, E-mail: 201638@cumtb.edu.cn

  • 中图分类号: TU43;TD164

Development and application of MWD system for DTH drilling rig

  • 摘要: 露天采矿随钻测量技术研究目前尚处于初期阶段,研发随钻参数测量系统是实现随钻测量技术服务于采矿工程行业的首要任务。基于TAIYE-390-Ⅱ潜孔钻机工作原理,确定了随钻参数,提出了一套通用的测量方法,并研发了测量装置。装置主要由数据采集系统、数据转换系统和数据处理系统组成,可实现对钻进时间、钻进深度、钻进速度、钻杆轴压、回转力矩、回转速度、冲击风压和冲击风量等参数的实时准确测量、采集、计算、存储、显示和远程传输。结合某石灰岩矿钻孔数据,提出了一种反映岩性变化的数据整合方法,并借助钻孔成像试验分析了随钻测量技术在岩体结构面识别中的可行性。结果表明钻进速度和回转力矩对岩体结构响应最为灵敏,冲击风量和冲击风压在钻凿含孔洞的破碎岩体时表现突出,而钻杆轴压和回转速度在岩体结构识别中作用较小。研究成果在潜孔钻机随钻测量系统研发和随钻岩性识别等方面具有重要的理论意义和应用价值。
    Abstract: The researches on the measurement while drilling (MWD) technology in the field of open-pit mining are still in the early stage, and the development of MWD systems is the primary task to realize the application of MWD technology in the mining engineering industry. Based on the working principles of the TAIYE-390-Ⅱ hydraulic DTH drilling rig, the MWD parameters and a universal measurement method are respectively determined, and the corresponding measurement devices are developed. The devices mainly include a data acquisition system, a data conversion system, and a data processing system, which can realize the real-time and accurate measurement, acquisition, calculation, storage, display and remote transmission of parameters such as drilling time, depth, penetration rate, axial force, rotation torque, rotation speed, percussive air pressure and air quantity. With the MWD data from a limestone mine, a data integration method that can reflect lithological variations is proposed, and the feasibility of MWD in the structural identification of rock mass is analyzed through the borehole imaging experiments. The results indicate that the penetration rate and rotation torque exhibit the highest sensitivity to rock mass structures. The air quantity and air pressure are more sensitive only when drilling fragmented rock masses with cavities, while the axial force and rotation speed have limited influences on the structural identification of rock mass. The research findings are of important theoretical significance and practical value in the development of MWD systems for DTH drilling rigs and lithological identification.
  • 岩体中结构面相互交错形成了复杂的块体系统,地下工程的开挖将会破坏原有系统的平衡,出现块体掉落或滑移现象,导致工程事故发生。这些不稳定的块体被称为关键块体。块体理论自正式提出以来,一直受到广泛关注。随着块体理论的发展,许多学者将其应用于岩体工程的各个领域,并取得了满意的效果。对于块体理论的研究,最早由Goodman等[1]提出,之后由Priest[2]、Delport等[3]和Lin等[4]进行了发展和补充,形成了不同的研究方向,提出了多种研究方法。不同的方法适用于不同的情况,但其本质是结构面的空间位置是否固定,由此分为定位结构面间确定性块体的识别和不定位结构面间潜在最大关键块体求解。

    定位结构面间的块体识别指的是已知结构面和开挖面的参数识别出区域内的所有块体,从而确定开挖面上的关键块体。对此类研究目前存在方法:①基于交线回路的拓扑分析方法[1-4];②基于平面对凸多面体切割算法而提出的组合方法[5-8]。然而在工程设计阶段,由于勘察技术及经济条件的限制,所获取的岩体内部结构面信息有限,仅能根据出露面信息获取其分组和产状信息。因此有必要考虑结构面和开挖面间所有的组合,并找到最不利情况,为工程开挖提供预测及评价,这在工程设计中具有重要意义。最大关键块体的概念最早由石根华提出,并通过向量法确定了块体的有限性和可动性,但由于当时计算机及三维可视化技术的限制,此方法并不便于应用和理解。为此,石根华提出了赤平投影方法,通过图形清晰展示了块体的有限性和可动性。然而,做图方法工作量大,精度有限。随着计算机技术的发展,张子新等[9]在此基础上提出了赤平解析法,将投影图用数学解析式表达,通过构建位置参量矩阵来判断块体的有限性和可移动性。赤平解析法是赤平投影法的发展和延伸,吸收了向量方法和做图方法的优点,但在计算位置参量时,需获取块体对应区域内的一点,这是解析法无法实现的,因此需要在赤平投影内采点试算,当代表块体的投影区较小时可能出现误差,需提高采点数量来弥补。本文在向量方法的基础上,提出了无效节理锥的删除方法,便于程序开发,且无精度问题。在应用方面,张奇华等[10]将最大块体用于地下厂房稳定性评估中;杨勇等[11]通过投影相切法对单一曲线开挖面进行了最大块体求解,并将其应用于隧道工程中。然而,对于多曲线开挖面上最大块体的求解方法却鲜有报道。另外在隧道工程中,关键块体的临空面均为曲面,计算机图形学中对于曲面的处理方式较多。在块体识别的文献中,石根华采用多平面替代开挖面;Elmouttie等[12]提出了一种由面单元组成的多面体开挖模型,可以描述曲面开挖;Li等[13]和Liu等[14]使用了解析表达式来描述曲面。然而在计算机图形学中,对元素的表达应简单化、单元化,这有利于程序开发和算法统一。

    鉴于此,本文针对柱状开挖面的隧道工程,提出了一种投影平移方法,将复杂的空间问题转化为简单的平面问题,实现了最大可移动块体的快速求解。在此基础上,采用离散-切割-组合方法完成了曲面块体的精细化构建。数据格式统一且简洁,易于存储和程序开发。基于此开发的三维可视化程序显示出算法具有良好的精度与鲁棒性。

    本文将勘察中易于获取的结构面产状视为已知,考虑结构面定位的随意性求解出潜在最大块体。在实际工程中存在3组和4组结构面的情况较为多见,实际上多于4组结构面所形成的最大块体体积均不会超过3组或者4组所形成的最大块体体积。因此本文以四组结构面岩体为研究对象,进行了潜在最大可移动块体的求解与构建。最大可移动块体是由节理锥与开挖面围限而成,因此本文针对节理锥的构建、定位以及曲面块体的重构3个方面进行详细阐述。

    本文采用Goodman等[1]提出的向量法构建节理锥。在含四组结构面的岩体工程中,以四面体和五面体为主。四面体是由其中3组结构面与开挖面形成的,而五面体是由四组结构面与开挖面形成的。4组结构面分别为F1,F2,F3和F4,那么可能存在5种组合情况,分别为F1-F2-F3-F4,F1-F2-F3,F1-F2-F4,F1-F3-F4和F2-F3-F4。因此在确定最大可移动块体时,应综合考虑3组及4组结构面参与形成的节理锥。

    节理锥(JP)是在结构面通过同一点条件下之间的交线组合而成。以交点为起点,在空间中每条交线存在两个方向,因此在不同结构面之间可以形成多种节理锥。本文采用向量法实现不同节理锥的构建。以每组结构面参数的均值为代表,则结构面法向量为

    ni=(sinαisinβi,cosαisinβi,cosβi) (1)

    式中:αiβi为第i组结构面平均倾向和倾角;i为结构面组号。

    (1)节理锥构建方法

    构造C2m×m阶矩阵M,其中m为参与构建节理锥的结构面组数量。当m=4时,M

    M=(M121M122M123M124M131M132M133M134M141M142M143M144M231M232M233M234M241M242M243M244M341M342M343M344) (2)

    m=3时,M

    M=(M121M122M123M131M132M133M231M232M233) (3)

    式中:

    Mijk=sgn(ni×njnk) (4)

    将节理锥JPr表示为JP(a1 a2  am),其中r为节理锥的编号,r=1,2,,2mai为节理锥的位置参数;ai=0表示节理锥位于第i条结构面的上方,ai=1表示节理锥位于第i条结构面的下方。

    按下式构造m×m阶对角矩阵:

    GJP(a1 a2... am)=diag(2a1+1,2a2+1,,2am+1) (5)

    令矩阵

    Br=MGJP(a1 a2 ... am)=(brst) (6)

    式中:BrC2m×m阶矩阵,其中s=1,2,…,C2mt=1,2,…,m;则JPr的棱线方向参数矩阵可表示为

    Er=diag(er1,er2,,erN) (7)

    式中:N=C2mersBr矩阵中第s行元素的最大值与最小值之和。4条结构面构建节理锥时,结构面之间存在六条交线,其方向向量可以用N表示为

    N=(n1×n2,n1×n3,n1×n4,n2×n3,n2×n4,n3×n4)T (8)

    3条结构面构建节理锥时,结构面之间存在3条交线,其方向向量可以用N表示为

    N=(n1×n2,n1×n3,n2×n3)T (9)

    则得到节理锥JPr棱线向量集的表达式为

    Vr=ErN (10)

    (2)删除无效节理锥

    Vr中的0向量代表此条交线不参与节理锥的构建,若非零向量的数量小于3则表示此节理锥为有限节理锥,无法与开挖面形成块体。将节理锥的所有棱线向隧道断面所在平面投影,若存在这样一条棱线,使得其余所有棱线的投影均位于其所在投影直线的一侧,则此节理锥可以与隧道开挖面形成块体,否则无法形成块体,为无效节理锥。如图 1所示,图 1(a)中所有棱线均位于红色虚线一侧,此为有效节理锥。而图 1(b)棱线投影均位于各虚线两侧,则此节理锥无法与开挖面形成块体,为无效节理锥。

    图  1  节理锥有效性判断
    Figure  1.  Validity judgment of joint pyramids

    在三维坐标系中,Y轴为隧道走向,Z轴竖直向上。形成最大块体时,开挖面与JP的边部棱线相切。最大块体的顶点是4条棱线的交点,即4条结构面的交点。当开挖面由单一曲线表示时,那么很容易通过相切条件确定棱线定位,但开挖面由多条曲线组合而成时,首先应判断边缘棱线与哪条曲线相切,之后才能进行定位。本文提出了一种边缘棱线与任意开挖曲面相切的定位方法,适用于任何凸开挖面。

    将节理锥的棱线向XOZ平面投影,通过最大夹角可以确定节理锥的边缘棱线。开挖面在XOZ平面上的投影由n条曲线组成,将曲线顺时针编号,令曲线Πj的表达式为

    Πjfj(x,z)=0 (11)

    图 2所示,曲线ΠjΠj+1相交于点Hj,其坐标为(xj, zj),其中j=1, 2, …, n。当j=1时, 下标j-1被n取代。当j=n时,下标j+1被1取代,若坐标系原点在开挖曲线内,则编号为j的曲线上(x, z)点处的单位外法向量为

    mj(x,z)=sgn(F)(fj(x,z)x fj(x,z)z)(fj(x,z)x)2+(fj(x,z)z)2 (12)
    图  2  切点的确定
    Figure  2.  Determination of tangential points

    式中:F为向量的方向判别式,

    F=fj(x,z)xx+fj(x,z)zz (13)

    那么曲线Πj两端点处的单位外法线向量分别为

    mj,1=mj(xj1,zj1)  (14)
    mj,2=mj(xj,zj) (15)

    由式(14),(15)可知,曲线在交点Hj处的法向量会发生突变,从mj,2变为mj+1,1n段曲线间存在n个交点,那么2n个向量将XOZ平面划分为2n个区域。假设u是节理锥边缘棱线投影的单位外法向向量,它指向节理锥投影的外部。若节理锥内侧任一棱线的投影为e,应满足:

    u=mj(x,z) ue<0} (16)

    如果向量u位于mj,2mj+1,1之间,那么切点为点Hj。当隧道开挖曲线为凸时,式(16)有唯一解,可计算出边缘棱线与开挖面的切点,从而将边缘棱线定位。之后可以求解出两个边缘棱线投影的交点。由于Y轴坐标不影响最大块体的形态,那么可将交点的Y轴坐标设定为0。将节理锥的所有棱线的起点设置为此交点,即可构建出最大可移动区域。

    根据被定位的棱线可求解出各棱线与开挖面的切点及交点,由切点、交点及节理锥的顶点可以构建出锥形块体。如图 3所示,六面体ABCDE为所构建的锥形块体,其6个面均为平面多边形。但实际开挖面为曲面,因此BCECDE面实际上属于同一个曲面。针对所构建的六面体ABCDE,本文采用离散-切割-聚合的方法实现曲面块体的重构,所得的曲面块体体积和滑动面面积更为精准。

    图  3  锥体的构建
    Figure  3.  Construction of joint pyramid

    (1)虚拟结构面

    本研究通过薄片单元的组合来构建曲面块体,而虚拟结构面则是一组用来将锥形块体切割离散为薄片单元的假想平面,并不真实存在,它是为了获取更多块体内部信息而设定的辅助平面。单元的离散是以平面对多面体的切割算法为基础,此算法针对凸多面体而设计,因此虚拟结构面仅适用于凸多面体。

    (2)曲面块体重构过程

    上一节所构建的锥体为凹多面体,需将其转化为两个四面体。图 3中六面体ABCDE为凹体,将其分解为四面体ABCEAECD,并分别重构。首先,设定一组通过同一点的放射型虚拟结构面,采用平面对凸多面体切割算法,将两个四面体分别进行离散。其次,采用曲面对凸多面体切割算法实现隧道开挖面对单元块体的切割,以开挖面为界,切割完成后的单元将被分为内外两部分。最后,判断被切割后的单元与开挖面的位置关系,开挖面外部的单元组合成为曲面块体,此时完成了曲面块体的重构。以图 3中的子块ACED为例进行说明,具体过程详见图 4

    图  4  曲面块体的重构
    Figure  4.  Reconstruction of surface blocks

    a)在开挖面曲线内部选定一点O,设置3个平面PEPDPC,使它们分别通过OEODOC且平行于隧道轴线。

    b)在平面PEPD之间设置一组平面,使他们均通过点O且与隧道轴线平行,相邻平面间夹角为φ,如图 4中的P1~P4。在PDPC之间设置另外一组平面,使他们均通过点O且与隧道轴线平行,相邻平面间夹角为η,如图 4中的P5,P6,令φη的值接近

    c)棱锥体AECD被所设置的平面分割为薄片单元,如图 4(a)中的B1~B8。

    d)开挖面将每个单元分割成内、外两部分。图 4(a)中,B4被分为B4a和B4b,开挖面外侧的所有单元被组合成为曲面块体,如图 4(b)中的Block。

    在平面对多面体切割的基础上,本文提出了柱状曲面对单元块体的切割算法,可实现开挖面的精细化构建。

    (1)单元中面的分类

    单元块体的数据按照点-线-面-体的格式进行存储,故切割的实质是面、线与曲面关系的判别。以简单的六面体单位为例,对单元面的类型进行说明。如图 5所示,y轴为隧道走向,在单元内部,平行于y轴的面为B型面,此类型的面是由虚拟结构面切割而产生的,其余的面为A型面。为切割计算准备,按照与切割曲面位置关系将单元的面划分为五类,其判断方法及图示见表 1图 5中的实心圆点代表边与开挖面的交点。

    图  5  面的类型示意图
    Figure  5.  Diagram of face type
    表  1  面的类型
    Table  1.  Types of face
    类别 判断方法 图示
    A型面 B型面
    Ⅰ类面 除落在曲面上的点外,其余顶点均在开挖面外侧
    Ⅱ类面 除落在曲面上的点外,其余顶点均在开挖面内侧
    Ⅲ类面 顶点分布在曲面内外两侧
    Ⅳ类面 位于开挖面外,以边的方式与开挖面接触
    Ⅴ类面 位于开挖面内,以边的方式与开挖面接触
    下载: 导出CSV 
    | 显示表格

    (2)曲面切割单元块体流程

    由于B型面平行于隧道轴线,不会影响单元与开挖曲面间的位置关系,在判断位置关系时可用A型面代替单元块体,这将判断点的数量从八减少到四,从而减少运算量。当A型面为Ⅰ类时,单元块体位于开挖面以内,无需切割;当A型面为Ⅱ类时,单元块体位于开挖面以外,无需切割;当A型面为Ⅲ类时,单元块体与开挖面相交,需要对其按如下方式进行切割。

    曲面将与其相交的单元块体切割为内外两部分,因此需定义U1和U2两个几何,分别存储切割后位于曲面外侧和内侧的两部分子面,定义面单元变量P用于存储切割形成的新面,即切口,对不同的面进行如下操作:Ⅰ类面:完整的放置于U1中。Ⅱ类面:完整的放置于U2中。Ⅲ类面:开挖曲面与此类面存在两个交点,两交点的连线将此面切割为两部分,按照与开挖面的位置关系被分别放入U1和U2中,除此之外将两交点连线放入P中。如图 6所示,Ⅲ类面EFGH与曲面相交于AB两点,直线AB将其分割为两部分,AFGHB面放入U1中,ABE面放入U2中,形成的新边AB放入P中,用以组建单元的切口面。Ⅳ类面:完整的放入U1中,并且将与曲面接触的边放入P中。Ⅴ类面:完整的放入U2中。将Ⅳ,Ⅴ类面与曲面接触的边放入P中。

    图  6  曲面切割单元块体
    Figure  6.  Curve surface cutting element block

    开挖曲面对单元内所有面切割完成后,P内存储的边形成闭合回路ABCD,为一平面四边形。将其分别放入U1和U2中,此时U1和U2将形成封闭体系,即切割后形成的两个块体。

    采用VC++实现了本文所提算法,通过OpenGL完成可视化,开发了隧道最大可移动块体程序MKB V1.0。

    (1)以4组结构面为例,在马蹄形开挖面的隧道工程中进行测试。图 7展示了MKB V1.0构建的五面体锥体。由2.1节可知,所构建的五面体是由两个四面体(Ⅰ和Ⅱ)组合而成。

    图  7  锥形块体的构建
    Figure  7.  Construction of pyramid blocks

    (2)选取虚拟结构面间夹角φ=3°,采用虚拟结构面对块体Ⅰ和Ⅱ进行切割,将两个块体切割为薄片单元,每个单元被开挖面分割为内、外两部分,将所有外部单元组合为曲面块体,图 8展示了曲面块体的重构过程。可以看出本文提出的方法可以较好地描述曲面,图 8中显示块体曲面与开挖面吻合度较高,块体的边与开挖面相切。计算结果表明,重构后与重构前的体积比为0.513,块体失稳形式为掉落,无滑动面。重构后与重构前锥体侧面面积的比值分别为0.792,0.813,0.745,0.832。重构前后块体体积和侧面积比均有显著差异,这将影响沿面滑动块体的安全系数。

    图  8  曲面块体重构
    Figure  8.  Reconstruction of curve surface blocks

    (3)以1°为间隔,将隧道走向从0°变换到180°,可以得到两条关于隧道走向的曲线(图 9)。①重构后与重构前块体的体积比曲线;②重构后与重构前滑动面面积比曲线。图 9中显示,不同走向时块体体积比值均在0.4左右,滑动面的面积比值主要分布在0.5~0.9之间。在面积比曲线中,比值为0表示无滑动面,块体失稳模式为掉落。重构前后块体的参数差异较大,因此重构曲面块体是必要的。

    图  9  重构后和重构前参数的比值
    Figure  9.  Ratios of parameters before and after reconstruction

    辽宁省内某隧道,开挖面为马蹄状,半径为8.8 m。该区岩体主要成分为太古代混合岩片麻岩,节理裂缝发育,节理略开放,开放节理多充填碎屑和泥浆。该区域有4组结构面,参数取每组结构面的平均值,见表 2

    表  2  结构面参数
    Table  2.  Parameters of the discontinuities
    节理组 产状代表值 走向
    F1 NE15°∠52° NW105
    F2 SW242°∠43° NW152°
    F3 SW203°∠39° NW113°
    F4 NW347°∠45° NE77°
    下载: 导出CSV 
    | 显示表格

    此工程存在四组节理,在构建块体时存在五种结构面的组合情况,即A组F1-F2-F3-F4;B组F1-F2-F3,C组F1-F2-F4;D组F1-F3-F4和E组F2-F3-F4。以隧道走向为0度为例,对最大可移动块体进行求解。求得节理锥集合的交线向量矩阵N和方向参数矩阵H分别为

    N=(N1N2N3N4N5N6)=(+0.75650.5217+0.3944+0.94820.3099+0.0690+0.30470.6470+0.6990+0.4048+0.6671+0.62540.7939+0.33640.50640.9625+0.05120.2664)  (17)
    H=(e1e2e3e4e5e6)=(011001100110011010100101101001011100001111000011101010100101010111001100001100110000111111110000) (18)

    方向参数矩阵H中,每列代表一个节理锥内不同棱线的方向参数,由4条节理构成的锥体列号分别为1,6,7,10,11和16,其中,第7和第10列所代表的节理锥无法与开挖面形成块体,为无效节理锥。在矩阵H中通过式(7)构建矩阵Er,结合式(8),(10)可求得节理锥所含棱线集的向量Vr。以JP-Xr表示节理锥,其中X为组号,r的含义见1.1节。在A组中,JP-A1和JP-A16顶角较大,与开挖面围成的块体体积较小,对工程稳定无影响。将A组中有效节理锥的基本信息展示在表 3中。以同样方法可求得3条结构面参与构建的节理锥,其中节理锥JP-B2,JP-B3,JP-B5,JP-B6,JP-C2,JP-C3,JP-C5,JP-C6,JP-D2,JP-D3,JP-D5,JP-D6,JP-E1,JP-E3,JP-E4和JP-E6的顶角较大,与开挖面围成的块体体积较小,对工程稳定无影响。将最大体块重构后在图 10中进行三维展示。在组成锥体的棱线中,若存在方向向下(其Z轴分量为负数)的棱线,那么在重力作用下存在滑落的可能性。若一个锥体内,所有棱线的Z向量均为正,那么此锥体仅在自重作用下是无法滑落的。因此,JP-A11,JP-B1,JP-C1,JP-D1和JP-E2仅在重力作用下无法滑动。4条结构面构建的最大块体为JP-A6,体积为120.7 m3,任意3条结构面构建的最大块体为JP-E5,体积为85.5 m3

    表  3  F1-F2-F3-F4构建的有效节理锥棱线向量
    Table  3.  Vectors of effective edges of joint pyramid constructed by F1-F2-F3-F4
    交线 棱线向量(x, y, z
    JP-A1 JP-A6 JP-A11 JP-A16
    F1F2 -N1 +N1
    F1F3 +N2 -N2 +N2 -N2
    F1F4 -N3 +N3
    F2F3 -N4 +N4
    F2F4 +N5 +N5 -N5 -N5
    F3F4 +N6 -N6
    体积小 体积小
    下载: 导出CSV 
    | 显示表格
    图  10  最大块体形态
    Figure  10.  Shape of maximum block

    单元块体的短边尺寸由虚拟结构面的夹角φ决定,夹角越小则块体的临空面越接近真实情况,但同时也将占用更大的内存和花费更多的计算时间,因此需要选择合适的夹角。在本工程中,将隧道走向θ分别设置为0°,45°,90°,135°,计算得到的块体体积随虚拟结构面的夹角φ值而变化。以φ为横坐标,以重构之后与重构之前的体积比为纵坐标得到曲线图 11图 11中显示随φ值的减小,体积比趋于一个定值。拟合4条曲线,R2均在0.997以上。因此,拟合曲线可以准确的预测不同隧道走向下φ=0°时的体积比。表 4展示了计算结果,并给出了φ为5°,0°时纵坐标的差值百分比,均小于5%。因此本工程将φ取值为5°,足以满足工程计算的精度。

    图  11  夹角与计算精度间的关系
    Figure  11.  Relationship between intersecting angle and calculation
    表  4  φ为5°,0°的体积比
    Table  4.  Volumetric ratios under φ of 5°and 0°
    隧道走向θ 45° 90° 135°
    φ=5° 0.4438 0.3395 0.4134 0.3784
    φ=0° 0.4404, 0.3357 0.4094 0.3748
    差值百分比/% 0.77 1.13 0.978 0.96
    下载: 导出CSV 
    | 显示表格

    由于结构面的隐蔽性,其力学参数很难准确获得,这给块体稳定性计算带来了困难。采用最大可移动块体的体积作为指标,对隧道走向进行评价。以1°为间隔,依次将隧道走向从0°改变到179°得到了两条角度与最大块体体积的曲线,如图 12。蓝色曲线代表 4条结构面构建的五面体的体积,红色曲线表示仅有3条结构面参与构建的最大四面体的体积。当隧道走向变化到特定值时,块体与开挖面相切的边将改变,导致最大块体的形成方式改变,从而引起体积变化率的突变。两条曲线中均有4次曲率的突变,θ分别为77°,105°,113°,152°。其中θ=113°时四面体体积的曲率变化并不明显。图 1314展示了θ为77°,105°,152°等3个极值点处的最大四面体和最大五面体形态和体积。当θ∈[0, 38]或θ∈[174,180)时,由4条结构面构建的五面体体积大于由任意3条结构面构建的四面体体积。而目前的文献均采用四面体进行最大块体的评估,从工程设计角度考虑,这种做法并不能涵盖所有情况。因此在求解最大可移动块体时,应该同时考虑四面体和五面体,从而进行包络设计。曲线图清晰地显示了隧道不同走向时块体的最大体积,在θ∈[77°, 105°]及θ=152°附近将形成体积较大的可移动四面体,因此在选择走向时应尽量避免。

    图  12  最大块体体积与隧道走向间的关系
    Figure  12.  Relationship between volume of maximum block and tunnel direction
    图  13  曲线极值处四面体形态
    Figure  13.  Tetrahedral morphologies at extreme of curve
    图  14  曲线极值处五面体形态
    Figure  14.  Pentahedral morphologies at extreme of curve

    MKB V1.0仅考虑形成最大块体的边界结构面。若结构面密度较大,则其余结构面被认为是块体内部裂隙。边界结构面的产状取本组结构面产状的均值。当4组以上结构面相互切割时,除单个锥体外,还可能存在两个及以上锥体相交的情况,此时它们联合可能会形成复杂的可移动区域,如图 15中的锥体Ⅰ和Ⅱ,对此进行讨论。

    图  15  存在交集的两锥体
    Figure  15.  Two pyramids with intersection

    (1)当锥体顶点距离较近甚至重合时,可将它们组合为一个大的锥体,这个大锥体是可以通过MKB V1.0求得的。如图 16所示,当5个结构面均通过A点(或接近)时,结构面P5(蓝色阴影面)两侧的锥体A-EDCGFA-GFB将组合为大椎体A-BCDE,这个锥体可利用结构面P1,P2,P3和P4求解得到。

    图  16  锥体的合并
    Figure  16.  Conflation of pyramids

    (2)无法忽略锥体顶点间距离时,锥体间的联合可能会形成复杂的可移动区域。此时图 15中的组合区域与锥体顶点的距离有关,这是定位结构面所面临的问题。需要采用概率模型模拟生成所有结构面。当确定了每条结构面的参数后,再根据确定信息进行块体搜索,通过大量模拟获取最大区域。MKB V1.0考虑不同结构面组合下开挖面上潜在的最大可移动块体,以最大块体体积为指标从块体理论角度来评价开挖面方向,因此它不是一种确定性的块体识别方法,是一种定量的风险评价指标。

    MKB V1.0将勘察中易于获取的结构面产状信息作为已知参数,考虑结构面定位的任意性,求解不同组之间可能形成的最大关键锥形块体。以最大块体体积为指标从块体理论角度来评价开挖方向的风险,它不是一种确定性的块体识别方法,而是一种风险评价指标。

    (1)提出的多曲线开挖面下的结构面定位方法可以有效的确定最大可滑动区域,通过算例对此方法进行了验证,显示了此方法对于多曲线组合而成的柱状开挖面的适用性。

    (2)在平面与曲面切割凸多面体算法的基础上,采用了一种离散-切割-组合算法来重构曲面块体,三维可视化算例显示重构的曲面与开挖面匹配度较高。与简化方法相比,采用此方法获得的滑动面面积以及块体体积更为准确。算例显示,重构后与重构前块体的体积比为0.4左右,滑动面面积比为0.5~0.9,重构前后差异较大。因此,有必要对块体进行重构。此外,文中是以多面体为单元进行运算,多面体数据存储格式为“点-线-多边形-多面体”,统一了数据格式,提高了数据处理的效率,也符合计算机图形学的基本思想。

    (3)通过工程实例进行了最大块体求解,其中包括四组结构面参与构建的五面体以及任意3组结构面参与构建的四面体,结果显示在隧道走向相同的情况下,存在五面体体积大于四面体体积的情况,因此在对隧道进行最大可移动块体分析过程中不能忽视五面体的影响,应综合考虑两种块体。

    (4)本文所提方法尚无法考虑以下两种情况:①由于结构面定位的不确定性,两个以上锥体存在重合区域时,求解潜在最大联合区域是一项较为复杂的工作,目前鲜有文献关注,是一项值得研究的问题,但应在结构面统计分布的基础上对其进行研究才具有工程意义。②非平面结构面之间形成的块体。非平面的凹凸结构面具有非规则性,多采用离散的三角单元或样条曲线曲面进行描述,它对围岩的切割应建立在围岩离散单元基础之上。本文仅利用这种思想做了初步探索,采用开挖曲面对离散后的薄片单元进行切割,尚未实现非平面结构面对围岩的切割,但离散-切割思想应用于此方面是一项值得研究且有意义的工作。

  • 图  1   潜孔钻机随钻测量系统组成

    Figure  1.   Installation position of sensors of TAIYE-390-Ⅱ MWD system

    图  2   潜孔钻机随钻参数测量系统原理

    Figure  2.   Principles of MWD system based on TAIYE-390-Ⅱ DTH drilling rig

    图  3   随钻测量系统数据采集逻辑框图

    Figure  3.   Logical flowchart of data acquisition for MWD system

    图  4   原始的深度序列随钻测量数据

    Figure  4.   Raw MWD data of depth series

    图  5   整合后的深度序列随钻测量数据及炮孔成像图

    Figure  5.   Integrated MWD data of depth series and borehole imaging

    图  6   钻孔成像现场试验

    Figure  6.   Borehole imaging tests

    图  7   整合后的深度序列随钻测量数据及孔洞岩体炮孔成像图

    Figure  7.   Integrated MWD data of depth series and borehole imaging of cavity rock mass

    表  1   TAIYE-390-Ⅱ潜孔钻机参数表

    Table  1   Working parameters of TAIYE-390-Ⅱ DTH drilling rig

    参数 参考值 参数 参考值
    钻孔直径/mm 115~203 液压泵工作压力/MPa 35
    钻进深度/m 21 推进/提升力/N 18000
    钻杆数量/根 7 回转力矩/(N·m) 3300
    单次推进/m 3 回转速度/RPM 120
    滑架行程/mm 4115 工作风压/MPa 1.05~2.46
    下载: 导出CSV

    表  2   随钻参数定义

    Table  2   Definition of MWD parameters

    名称 定义
    钻进深度 沿钻孔方向钻头底部离地表的轴向长度
    钻进速度 钻头钻进岩体的速度
    钻杆轴压 施加于钻头上的轴压力
    回转力矩 钻头转动所需要的的扭矩
    回转速度 钻头每分钟转动的圈数
    冲击风压 作用于冲击器上的风的压力
    冲击风量 作用于冲击器上的风的流量
    下载: 导出CSV
  • [1] 刘红岐, 张元中. 随钻测井原理与应用[M]. 北京: 石油工业出版社, 2018.

    LIU Hongqi, ZHANG Yuanzhong. Principle and Application of Logging While Drilling[M]. Beijing: Petroleum Industry Press, 2018. (in Chinese)

    [2] 侯仕军, 丁伟捷, 田帅康, 等. 随钻测量技术在非油气工程领域的应用现状与展望[J]. 矿业研究与开发, 2022, 42(12): 41-49. https://www.cnki.com.cn/Article/CJFDTOTAL-KYYK202212009.htm

    HOU Shijun, DING Weijie, TIAN Shuaikang, et al. Application status and prospects of MWD technology in non-oil and gas engineering field[J]. Mining Research and Development, 2022, 42(12): 41-49. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-KYYK202212009.htm

    [3]

    KHORZOUGHI B M, HALL R, APEL D. Rock fracture density characterization using measurement while drilling (MWD) techniques[J]. International Journal of Mining Science and Technology, 2018, 28(6): 859-864. doi: 10.1016/j.ijmst.2018.01.001

    [4]

    MANZOOR S, LIAGHAT S, GUSTAFSON A, et al. Establishing relationships between structural data from close-range terrestrial digital photogrammetry and measurement while drilling data[J]. Engineering Geology, 2020, 267: 105480. doi: 10.1016/j.enggeo.2020.105480

    [5]

    KOSOLAPOV A I. Modern methods and tools for determining drillability and blastability of rocks[J]. IOP Conference Series: Earth and Environmental Science, 2020, 459(2): 022097. doi: 10.1088/1755-1315/459/2/022097

    [6] 韩新平, 邹伟, 侯成恒. 基于BP神经网络的回转切削钻机凿岩速度预测[J]. 应用泛函分析学报, 2015, 17(1): 86-90. https://www.cnki.com.cn/Article/CJFDTOTAL-YYFH201501013.htm

    HAN Xinping, ZOU Wei, HOU Chengheng. Prediction on drilling speed of rotary drilling rigs based on BP neural network[J]. Acta Analysis Functionalis Applicata, 2015, 17(1): 86-90. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YYFH201501013.htm

    [7] 邹伟. 基于潜孔钻机钻进参数的岩石硬度预测方法研究[D]. 阜新: 辽宁工程技术大学, 2015.

    ZOU Wei. Study on Drediction Method of Rock Hardness Based on Drilling Parameters of DTH Drilling Rig[D]. Fuxin: Liaoning Technical University, 2015. (in Chinese)

    [8] 岳中琦. 钻孔过程监测(DPM)对工程岩体质量评价方法的完善与提升[J]. 岩石力学与工程学报, 2014, 33(10): 1977-1996. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201410005.htm

    YUE Zhongqi. Drilling process monitoring for refining and upgrading rock mass quality classification methods[J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(10): 1977-1996. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201410005.htm

    [9] 李宁, 李骞, 宋玲. 基于回转切削的岩石力学参数获取新思路[J]. 岩石力学与工程学报, 2015, 34(2): 323-329. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201502013.htm

    LI Ning, LI Qian, SONG Ling. Acquiring mechanical parameters of rock based on rotational cutting[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(2): 323-329. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201502013.htm

    [10] 王琦, 高红科, 蒋振华, 等. 地下工程围岩数字钻探测试系统研发与应用[J]. 岩石力学与工程学报, 2020, 39(2): 301-310. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202002010.htm

    WANG Qi, GAO Hongke, JIANG Zhenhua, et al. Development and application of a surrounding rock digital drilling test system of underground engineering[J]. Chinese Journal of Rock Mechanics and Engineering, 2020, 39(2): 301-310. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202002010.htm

    [11]

    LAKSHMINARAYANA C R, TRIPATHI A K, PAL S K. Experimental investigation on potential use of drilling parameters to quantify rock strength[J]. International Journal of Geo-Engineering, 2021, 12(1): 23. doi: 10.1186/s40703-021-00152-5

    [12]

    RAJESH KUMAR B, VARDHAN H, GOVINDARAJ M. Prediction of uniaxial compressive strength, tensile strength and porosity of sedimentary rocks using sound level produced during rotary drilling[J]. Rock Mechanics and Rock Engineering, 2011, 44(5): 613-620.

    [13]

    KALANTARI S, BAGHBANAN A, HASHEMALHOSSEINI H. An analytical model for estimating rock strength parameters from small-scale drilling data[J]. Journal of Rock Mechanics and Geotechnical Engineering, 2019, 11(1): 135-145.

    [14] 郭勇, 周振华. 潜孔钻机的应用现状与发展趋势[J]. 矿业快报, 2008(4): 13-15. https://www.cnki.com.cn/Article/CJFDTOTAL-KYKB200804006.htm

    GUO Yong, ZHOU Zhenhua. Application actualities and development trend of down-the-hole drill[J]. Express Information of Mining Industry, 2008(4): 13-15. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-KYKB200804006.htm

    [15] 徐小荷, 余静. 岩石破碎学[M]. 北京: 煤炭工业出版社, 1984.

    XU Xiaohe, YU Jing. Rock Fragmentation[M]. Beijing: China Coal Industry Publishing House, 1984. (in Chinese)

    [16] 赵宏强, 李美香, 高斌, 等. 潜孔钻机凿岩过程自动防卡钻理论与方案研究[J]. 机械科学与技术, 2008, 27(6): 739-743. https://www.cnki.com.cn/Article/CJFDTOTAL-JXKX200806011.htm

    ZHAO Hongqiang, LI Meixiang, GAO Bin, et al. On automatic anti-jamming control in the rock drilling process for a down the hole drill[J]. Mechanical Science and Technology for Aerospace Engineering, 2008, 27(6): 739-743. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JXKX200806011.htm

图(7)  /  表(2)
计量
  • 文章访问数:  0
  • HTML全文浏览量:  0
  • PDF下载量:  0
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-05-26
  • 网络出版日期:  2024-07-11
  • 刊出日期:  2024-06-30

目录

/

返回文章
返回