Experimental study on measurement of rock strength parameters based on digital drilling technology
-
摘要: 快速准确地超前获取岩体力学参数是指导岩体工程安全施工的重要前提。利用自主研发的多功能数字钻进测试系统,进行了不同强度等级完整试样的室内数字钻进试验及单轴压缩试验研究。根据试验分析结果,基于岩石研磨钻进过程中金刚石钻头的受力特点,利用力学极限平衡与能量守恒原理,定义并推导出了岩石单位体积研磨能的理论解,获得了数字钻进参数与岩石单位体积研磨能的数学方程式,进而建立了岩石数字钻进参数与单轴抗压强度的定量关系,即RDP-Rc预测模型,并进一步验证了该模型的可行性与有效性。研究方法与结果可以为工程岩体强度参数的实时原位测定提供一种新的途径。Abstract: Quick and accurate advance acquisition of mechanical parameters of rock mass is an important prerequisite for guiding the safe construction of rock mass engineering. Using the self-developed multi-functional digital drilling test system, the indoor digital drilling tests and uniaxial compression tests on complete samples with different strength grades are carried out. Based on the test results and according to the stress characteristics of diamond bit in rock grinding drilling process, the theoretical solution to the rock grinding energy per unit volume is defined and deduced by using the force limit equilibrium and the energy conservation principle. The mathematical equation for digital drilling parameters and grinding energy per unit volume of rock is obtained, and the quantitative relationship between rock digital drilling parameters and uniaxial compressive strength is established, namely RDP-Rc prediction model, and the feasibility and validity of the model are further verified. The research methods and results can provide a new way for the real-time in-situ measurement of strength parameters of rock mass.
-
0. 引言
粉质黏土/黏质粉土都是粉粒与黏粒的混合土,是岩土工程中广泛存在的一大类土体,如中国南海神狐海域天然气水合物储层[1]、长江中下游沿岸漫滩相淤泥质粉质黏土等;其物理、力学特性与黏粒含量密切相关,本质上是非塑性粉粒和塑性黏粒间相互作用的宏观表现。以淤泥质粉质黏土为例,黏土片之间存在范德华相互作用(引力)和双电层相互作用(斥力)、黏土片与粉粒之间也存在范德华相互作用(引力),这些分别表现为引力和斥力的物理–化学作用在复杂的微观架空结构上取得平衡;土体受力变形将引起土颗粒排列结构调整和上述物理–化学作用的相应变化,在宏观上表现出结构性、高压缩性、各向异性等一系列力学性质。其中,范德华相互作用是土体真黏聚力的来源,是粉粒与黏粒间相互作用的本质特点,对于从微观机制上认识粉质黏土/黏质粉土这一大类土体的宏观力学特性具有关键作用。
以离散单元法为代表的数值模拟方法是当前研究土体宏微观力学特性的重要手段[2],而范德华作用计算方法是离散元模拟黏性土的关键。Jiang等[3]推导了考虑接触宽度的圆盘、球体间范德华力以模拟月壤的真黏聚力;Yao等[4]在三维离散元模拟黏土中采用解析方法计算矩形薄片间的范德华力。当前离散元方法主要用于模拟纯粗粒土(如净砂)或纯黏土,尚无关于同时包含黏土片和粉粒的离散元模拟报道,也无考虑黏土片与球状颗粒间范德华力的研究。为此,本文提出一种黏土片——球状颗粒间范德华作用简便计算方法,并通过典型算例说明了计算方法的精确性和高效性。
1. 宏观体间范德华作用求解方法
两分子间范德华吸引能表达式最初由London[5]给出
u(→r)=−B|→r|6, (1) 式中,B为伦敦常数,
→r 为连接两分子中心的矢量,|→r| 为两分子间距离。实际上,完整的范德华作用包括吸引和排斥两部分,分别与距离的6次方(如式(1))和12次方成反比,都是分子内部电荷分布不均匀在不同尺度的表现。对于离散元模拟黏性土,仅式(1)给出的吸引能具有实际意义,后文将进一步说明。然而,式(1)只适用稀薄介质中的范德华作用。对于凝聚态物质(如土体中各类固体矿物),范德华作用能随分子间距离的衰减比式(1)更快[6-7]。对此,Anandarajah等[8]对式(1)进行了修正以考虑范德华作用的迟滞效应(retardation effect):
u(→r)=−Bc|→r|6(|→r|+c), (2) 式中,常数c=49.363 nm。
两分子间范德华力可通过能量对距离求导得到
→f(→r)=−∂u(→r)∂→r。 (3) 对如图1所示的分子密度均匀的两宏观体,它们之间的范德华力可以通过对式(3)体积分得到[9-10]
→F=ρ1ρ2∫V1∫V2→f(→r) dΩ1dΩ2, (4) 式中,V1和V2为两个宏观体的范围,
ρ1 和ρ2 为两宏观体中的分子密度。宏观体1受到宏观体2范德华作用形成的力矩可表示为
→M=ρ1ρ2∫V1∫V2→l×→f(→r)dΩ1dΩ2, (5) 式中,
→l 为从力矩作用点O1指向微元体中心的矢量。式(4),(5)积分只有在几种简单、规则的宏观体上才能得到解析解,如点、球体、半无限体之间[9, 11]。对于薄黏土片与其他形状实体之间的解析解则难以求解。本文采用蒙特卡罗模拟近似求解式(4),(5)。式(4)可改写为
→F=ρ1ρ2∫V1∫V2→f(→r)g1g2g1g2dΩ1dΩ2, (6) 式中,
g1 ,g2 可为任意函数。若选择合适的
g1 ,g2 使得∫V1g1dΩ1=1 ,∫V2g2dΩ2= 1,则式(6)可视为对函数ρ1ρ2→f(→r)/(g1g2) 求期望,且g1 ,g2 应为矢量→r 的两端点在宏观体1和宏观体2中的概率密度函数。对分子密度均匀的宏观体应有g1=1/V1 ,g2=1/V2 。根据期望的定义,式(6)也可表达为→F=ρ1ρ2V1V2N1N2limN1→∞N2→∞N1N2∑1→f(→r), (7) 式中,N1和N2分别为两个宏观体中的取样点数,且应有
N1:N2=ρ1V1:ρ2V2 。同理,式(5)可表示为
→M=ρ1ρ2V1V2N1N2limN1→∞N2→∞N1N2∑1→l×→f(→r)。 (8) 蒙特卡罗模拟的思想是通过计算式(7),(8)近似得到
→F 和→M 。本研究将黏土片简化为薄片状长方体,粉粒简化为球体,如图2所示。然而,由于→f(→r) 对距离非常敏感,在黏土片和球体两个宏观体上采样计算式(7),(8)将导致巨大的计算量。为此,采用以下方法大幅度降低计算量:(1)将黏土片视为由规则排列的基本立方体构成,如图2所示。根据范德华作用的可叠加性,只需求得基本立方体受到的球体颗粒的范德华作用,再将构成黏土片的所有基本立方体所受作用求和,即可得到整个黏土片受到的范德华作用。因此,待求解问题转变为求解基本立方体与球体间的范德华作用。
(2)借助点与球体间范德华作用的解析解,计算式(7),(8)时只需在基本立方体上采样即可,由此避免了在球体上采样。
此处将黏土片视为立方体的集合系基于以下考虑。水化后的黏土矿物一般由多个晶层堆叠而成,层间还可能存在阳离子与水分子。当不关心晶层间相互作用细节时,可将晶层和层间可能存在的阳离子与水分子视为薄片状整体,此即黏土片;而计算范德华力所需的哈梅克常数A(见后文公式(12),(13))可视为黏土片内各成分的某种加权平均。
2. 基本立方体与球体间范德华作用
2.1 点与球体间范德华作用
参考Anandarajah等[8]给出的两球体间的范德华力,可方便推得图3中点P受到球体的范德华引力
→f1(→r) ,其方向从点P指向球心。图3中dn表示点与球表面的净距。点与球体间范德华力大小为
|→f1(→r)|=b1+b23c4d2(d2−R2)4, (9) 式中,d为点到球心距离,R为球体半径,
b1=2cR(a1+a2+a3+a4+3R8),b2=3a5(c2−d2−R2)(d2−R2)4,a1=3(c−d)d7,a2=d3(3c+2d)(4c2−6cd+3d2)R2,a3=3c(4c−d)d2R4,a4=(−2c2+3cd−6d2)R6,a5=log[(d−R)(c+d+R)(d+R)(c+d−R)]。 图4为P点与半径为1 m的球体之间的范德华引力。图4表明,点与球体间范德华力随净距增大而快速非线性减小;当净距超过4 nm后,范德华力小到可以忽略。当净距趋近0时,范德华力将趋近无穷大。但实际上,当净距小于1 nm后,范德华斥力将快速增长并成为主导力,表现为对两宏观体相互贯入的排斥,故而当净距小于1 nm后范德华引力精确值对求解土颗粒间相互作用意义不大。Yao等[4]在采用离散元模拟黏土时则假设净距小于0.9 nm后范德华引力保持恒定。
图5给出了点与不同大小球体间范德华力的关系。可以看到,在相同净距下,点与半径大于1
μm 的球体间范德华力基本与半径无关,这是由于范德华吸引能与距离的6次方成反比,图3中点P与球体内最靠近该点的局部小范围分子间作用对整个范德华力的贡献最大,而较远位置处的分子贡献几乎为零。2.2 基本立方体与球体间范德华相互作用
在式(9)的点与球体间范德华力解析解基础上,本节采用蒙特卡罗模拟计算边长L=10 nm的基本立方体受到半径为R=1
μm 球体的范德华作用。以立方体中心为原点建立图6所示坐标系,通过改变球心位置反映立方体与球体间不同的相对位置关系。为便于表达,图6中几何体尺寸未按比例绘制。原则上,本文方法适用于任何尺寸的基本立方体且精度均可保证,而此处选择立方体边长为10 nm是考虑以下物理意义和计算效率两方面因素:①黏土片厚度一般在几十纳米范围,边长10 nm的立方体可方便地构建出一系列不同厚度的黏土片,若基本立方体再大则其构建黏土片的灵活性将降低;②若基本立方体尺寸更小,对图2中众多基本立方体与球体相互作用求和将非常费时。此处选择半径为1
μm 的球体基于以下两点考虑:①土体中非塑性颗粒半径一般不会小于1μm ;②根据图5可推知,边长为10 nm的基本立方体与半径大于1μm 的球体之间的范德华作用基本与半径无关(由于基本立方体与球体相比体积非常小),这使得本文的研究结果适用于黏土片与土体中任意粒径非塑性颗粒间范德华相互作用。为遍历基本立方体与球体间可能的相对位置关系并考虑到对称性,图6中
γ 和β 的范围均为[0,0.5π] ,球心位置用3个参数(m,n,d)表征,其中m=γ /(0.5π)∈[0,1] ,n=β /(0.5π)∈[0,1] ,d为球心与基本立方体中心距离。对γ ,β 分别均匀划分为20份,在各γ ,β 组合下,以d0表示基本立方体与球体刚好接触时两者中心距离,d的取值为d0+100.065i nm,i∈[0,20] ,d取值超出此范围后范德华作用几乎为零,不必再考虑。对3个参数逐一遍历计算,总算例数为213=9261个。借助式(9),采用蒙特卡罗模拟法计算式(7),(8)时不需在球体上采样,而仅需在立方体上采样。故式(7),(8)变形为
→F=ρ1V1limN1→∞1N1N1∑1→f1(→r), (10) →M=ρ1V1limN1→∞1N1N1∑1→l×→f1(→r), (11) 式中,
→l 为从立方体中心指向立方体中采样点的矢量。为便于后续分析,定义量纲归一化力和力矩为
→F*=cA→F, (12) →M*=1A→M, (13) 式中,
A=π2Bρ1ρ2 为哈梅克常数,与黏土矿物类型、黏土孔隙介质有关[11]。用蒙特卡罗模拟法求解式(10),(11)的方法是在基本立方体内随机采样N1个点,计算这些点处
→f1(→r) 、→l×→f1(→r) 的平均值;当采样点数N1越大,计算结果越逼近真实解。试算发现,立方体与球体间距越小,所需采样点越多。本文采用以下方法确定蒙特卡罗模拟在基本立方体中的采样点数:进行10次平行模拟,计算所得范德华力的大小误差在±0.5%范围内即可认为采样点足够;检验表明,此时蒙特卡罗模拟结果可认为收敛。经试算,采样点数需满足以下条件:立方体与球体间净距离小于等于1.2 nm时,取5×106个采样点;当净距大于1.2 nm而小于1.5 nm时,取1×106个采样点;当净距大于等于1.5 nm时取2×105个采样点。以存在解析解的半径为R′ =5 nm的球体(基本立方体的内切球)与半径为1 μm 的球体间范德华力为例(此时采样仅在半径为5 nm的球体内进行),图7给出了5次平行模拟的结果与解析解的对比。结果表明,上述采样点数的选择是合理的,保证了蒙特卡罗模拟具有预期的精度。图8为采用蒙特卡罗模拟得到的基本立方体与球体间范德华作用。图8(a)中m=n=0,即球体位于立方体正上方(+z方向),立方体受到的范德华引力沿+z方向,当距离d-R超过10 nm后范德华力小到可以忽略;x,y方向作用力以及3个方向的力矩几乎为零,与预期相符。图8(b)中m=0.2,n=0.8,球体对立方体的范德华作用表现为3个方向的力与力矩,当距离d-R超过15 nm后范德华力作用小到可以忽略。对比图8(a)和(b)可知,两宏观体之间的范德华作用大小与两者相对距离(d)和位置(m,n)都有关系。
对所有算例结果统计分析表明,立方体与球体间范德华力方向几乎沿着两者中心连线方向,这是由于基本立方体与球体相比非常小,立方体本身的尺寸对范德华力方向影响极小;但范德华力的大小与m,n,d三个参数都密切相关。将基本立方体受到的力矩大小与范德华力的大小比值定义为名义力臂。图9表明,名义力臂范围在0~2 nm,且绝大多数分布在0~1 nm,相对于立方体边长(10 nm)而言比较小。后文结果将表明,在黏土片与球体间相互作用简便计算方法中,不考虑这种力矩引起的误差很小。
2.3 人工神经网络模型拟合
为获得范德华力大小与m,n,d三个参数之间的函数关系,此处基于蒙特卡罗模拟结果训练一个人工神经网络模型来拟合这种函数关系。将9261个算例结果构成的数据集划分为70%训练集,15%验证集和15%测试集,经过多次尝试后确定的网络结构包括一个输入层(3个神经元,分别为m,n,d),一个隐藏层(10个神经元)和一个输出层(1个神经元,即范德华力的大小),采用Levenberg-Marquardt算法,经过279轮训练后均方误差稳定在10-4附近。图10给出了测试集目标值(即真实值)与神经网络预测值的比较。可以看到,数据点分布在斜率为1的直线附近很小范围,说明训练得到的神经网络模型拟合精确度足够高。
为便于读者编程实现,此处给出神经网络模型的运算过程与模型参数:
x1,3×1=xgain,3×3(xin,3×1−xoffset,3×1)+xb,3×1, (14a) x2,10×1=W1,10×3x1,3×1+b1,10×1, (14b) x′2,10×1=tanh(x2,10×1), (14c) y=W2,1×10x′2,10×1+b2, (14d) yout=(y−yb)/ygain+yoffset, (14e) 其中,
xin,3×1=[m,n,d]T 为输入数据,yout为输出的预测值,x1,3×1 ,x2,10×1 ,x′2,10×1 ,y为中间变量,其余为通过神经网络训练得到的模型参数:至此,基本立方体与球体间的范德华力大小和方向全部确定,并可忽略力矩作用。
3. 黏土片与球体间范德华作用简便计算方法
以图2为例说明黏土片与球体间范德华作用简便计算方法步骤。①将黏土片划分为边长为10 nm的立方体集合。②对每个立方体建立如图6所示局部坐标系,求得m,n,d 3个参数;当球心不位于局部坐标系的第一象限(x≥0,y≥0,z≥0)时,可利用坐标变换将球心变换至第一象限。③利用神经网络模型计算各立方体受到球体作用的范德华引力,并在整体坐标系下计算整个黏土片受到的合力与合力矩。
由于黏土片与球体间范德华作用无法求得解析解,为验证本文简便方法的可靠性,图11对比了简便计算方法与蒙特卡罗模拟得到的结果,其中黏土片尺寸为长度L=200 nm,宽度W=200 nm,厚度H=10 nm,球心x,y坐标分别为80,20 nm。为保证模拟精度,此处蒙特卡罗模拟在黏土片中的采样点数为边长10 nm立方体中采样点数的400倍(即两者的体积比)。图11表明,随着黏土片与球体净距的增大,范德华力和力矩均快速减小;当净距超过10 nm后,两实体间范德华作用小到可以忽略。图11还表明,简便计算方法尽管忽略了基本立方体与球体间的力矩作用,且神经网络模型拟合也必然存在一定误差,但其计算结果与蒙特卡罗模拟结果基本一致;相对误差较大处的范德华力和力矩本身已较小。图11中简便计算方法耗时约为蒙特卡罗模拟耗时的万分之一。因此,本文提出的简便算法精度高且能极大节省计算时间,对大规模高效率离散元数值模拟的实现具有重要价值。
图12进一步讨论了黏土片尺寸与其所受范德华力的关系。图中,黏土片保持L=W=20H的三向尺寸比例,与球体净距离保持为1.82 nm,黏土片厚度从10 nm增长到70 nm,涵盖黏土片的常见尺寸范围。图12表明,当黏土片尺寸增大时,所受的范德华力呈现增长率渐增的非线性增长趋势。当参数A取1×10-19 J(高岭石—水—氯化钠电解质体系典型值)[5],黏土矿物密度取为2.63 g/cm3时,黏土片所受范德华引力与其自身重力的比值在107~108量级,可见准确计算范德华引力对离散元模拟黏性土的准确性至关重要。
4. 结语
本文针对黏土片与球形颗粒间范德华作用难以解析求解的困难,提出一种简便计算方法,其本质是计算构成黏土片的各基本立方体与球状颗粒间的范德华力,最后矢量求和获得黏土片受到的合力、并计算合力矩。各基本立方体与球状颗粒间的范德华力可采用基于蒙特卡罗模拟数据集训练得到的神经网络模型快捷计算得到。由此避开了解析解不可得的难题,也避免了重复进行蒙特卡罗模拟。采用简便方法计算的典型算例表明,对平面尺寸为200 nm×200 nm、厚度为10 nm的黏土片,与球状颗粒的最小净距在10 nm范围内应考虑两者间范德华作用。
本文给出了简便计算所需的所有参数,包括哈梅克常数A(可参考文献[11]确定)、常数c(c=49.363 nm)及公式(14)中的拟合参数(已由神经网络训练得到),读者可植入任何程序开发使用。后续研究将应用本文方法进行粉质黏土/黏质粉土的大规模(数十万颗粒)离散元数值模拟,以解决离散元模拟计算量巨大的难题。
-
表 1 不同强度等级砂浆配比参数
Table 1 Ratio parameters of mortar with different strength grades
砂浆强度等级 水泥/(kg·m-3) 标准砂/(kg·m-3) 水/(kg·m-3) M10 275 1450 300 M15 310 M20 370 注: 水泥强度等级均为42.5。表 2 初步探讨性数字钻进试验相关参数
Table 2 Related parameters of preliminary exploratory digital drilling tests
试样类型 编号 转速n/(r·min-1) 钻进压强F/MPa 扭矩M/(N·m) 钻进速度v/(mm·min-1) Rc/MPa /MPa 最优钻进速度 /(mm·min-1) M10 J1 425 2.72 0.271 1.368 3.84 4.82 2.81 J2 425 3.74 0.554 2.058 4.85 J3 425 4.76 0.812 2.706 5.43 J4 425 5.78 1.013 2.724 4.96 J5 425 6.8 1.148 2.376 5.12 J6 425 7.82 1.238 2.202 4.74 M15 S1 425 2.72 0.762 1.14 13.22 14.29 3.03 S2 425 3.74 1.024 1.644 14.67 S3 425 4.76 1.318 2.322 14.83 S4 425 5.78 1.634 2.952 15.06 S5 425 6.8 1.962 2.688 14.24 S6 425 7.82 1.953 2.538 13.71 M20 I1 425 2.72 1.094 0.882 19.68 19.71 3.35 I2 425 3.74 1.708 1.278 20.46 I3 425 4.76 2.238 1.878 18.74 I4 425 5.78 2.613 2.616 18.62 I5 425 6.8 2.961 3.306 20.91 I6 425 7.82 3.023 2.874 19.87 KII231 K1 425 2.72 0.932 0.648 86.93 98.42 1.36 K2 425 3.74 0.854 0.93 118.46 K3 425 4.76 1.195 1.116 109.14 K4 425 5.78 1.481 1.308 96.62 K5 425 6.8 1.654 1.35 92.14 K6 425 7.82 1.843 1.142 87.25 KII237 L1 425 2.72 0.096 0.192 38.97 56.78 0.42 L2 425 3.74 0.138 0.258 61.34 L3 425 4.76 0.165 0.336 56.25 L4 425 5.78 0.213 0.391 68.46 L5 425 6.8 0.432 0.348 35.75 L6 425 7.82 0.564 0.288 79.87 KII243 X1 425 2.72 0.987 0.954 134.57 158.51 2.92 X2 425 3.74 1.345 1.554 157.46 X3 425 4.76 1.867 2.184 164.75 X4 425 5.78 2.314 2.886 167.64 X5 425 6.8 2.542 2.838 171.43 X6 425 7.82 2.713 2.346 155.18 YXQ300 Y1 425 2.72 0.823 0.804 114.63 113.68 2.37 Y2 425 3.74 1.125 1.176 106.82 Y3 425 4.76 1.438 1.656 127.46 Y4 425 5.78 1.756 2.22 102.42 Y5 425 6.8 2.013 2.31 132.02 Y6 425 7.82 2.234 1.974 98.74 YXQ304 T1 425 2.72 0.712 0.522 92.14 102.63 2.14 T2 425 3.74 1.275 1.302 110.14 T3 425 4.76 1.578 1.872 122.45 T4 425 5.78 1.982 2.082 98.28 T5 425 6.8 2.176 1.803 83.85 T6 425 7.82 2.304 1.554 108.94 表 3 系统研究性数字钻进试验相关参数
Table 3 Related parameters of systematic research digital drilling tests
试样类型 编号 转速n/(r·min-1) 钻进压强F/MPa 扭矩M/(N·m) 钻进速度v/(mm·min-1) Rc/MPa M10 A1 252 2.72 0.615 1.098 4.21 A2 252 3.74 0.948 1.59 5.25 A3 252 4.76 1.184 2.136 3.10 A4 425 2.72 0.42 1.368 3.16 A5 425 3.74 0.655 1.878 3.91 M15 B1 252 2.72 0.904 0.78 15.42 B2 252 3.74 1.475 1.29 16.24 B3 252 4.76 1.813 1.728 14.05 B4 425 3.74 1.158 1.584 15.68 B5 425 4.76 1.418 2.322 13.59 M20 C1 252 3.74 1.99 1.104 19.78 C2 252 4.76 2.496 1.446 21.24 C3 425 2.72 0.81 0.882 18.92 C4 425 3.74 1.308 1.458 17.63 C5 425 4.76 1.938 1.998 20.12 灰岩 D1 252 4.76 0.78 0.186 47.89 D2 252 6.8 3.604 1.302 47.53 D3 252 7.82 4.485 1.032 58.31 D4 425 3.74 0.953 0.258 92.23 D5 425 5.78 1.643 0.39 88.35 D6 425 6.8 1.838 0.396 122.54 D7 425 7.82 2.063 0.462 100.78 D8 721 4.76 4.853 2.478 68.34 D9 721 5.78 3.248 2.406 68.73 D10 721 6.8 6.78 2.952 119.8 D11 721 7.82 3.248 2.076 85.33 花岗岩 H1 252 5.78 4.3 0.618 91.36 H2 425 4.76 3.548 0.81 148.75 H3 425 5.78 4.546 0.978 121.27 H4 425 6.8 5.18 1.068 164.14 H5 721 5.78 4.398 1.428 175.46 H6 721 6.8 4.887 1.506 194.58 H7 721 7.82 4.52 1.386 186.32 H8 1035 6.8 3.278 1.344 136.43 表 4 数字钻进过程钻进能量计算结果
Table 4 Calculated results of drilling energy of digital drilling process
试样类型 编号 轴力功 扭矩功 底面摩擦耗能量 侧面摩擦耗能量 水动能 单位体积研磨能 WF/10-5 kJ WM/kJ Ed/10-4 kJ Ec/10-2 kJ Ek/10-13 kJ ηe/(kJ·cm-3) M10 A1 0.88 0.016 0.75 0.50 0.47 3.454 A2 1.75 0.025 1.03 0.69 0.37 3.855 A4 1.10 0.019 1.26 0.84 0.40 2.520 A5 2.07 0.029 1.74 1.16 0.41 3.149 M15 A6 0.63 0.024 0.75 0.50 0.89 8.175 A7 1.42 0.039 1.03 0.69 0.87 8.414 A8 2.42 0.048 1.31 0.87 1.06 7.660 A10 3.26 0.063 2.21 1.47 0.33 7.042 M20 A11 1.22 0.053 1.03 0.69 0.29 14.011 A12 2.03 0.066 1.31 0.87 0.35 13.388 A13 0.71 0.036 1.26 0.84 1.15 10.589 A14 1.61 0.058 1.74 1.16 1.46 10.817 灰岩 D1 0.26 0.021 1.31 0.87 0.49 21.385 D3 2.38 0.118 2.15 1.44 0.84 34.152 D4 0.28 0.042 1.74 1.16 3.03 40.321 D6 0.79 0.082 3.16 2.11 2.66 51.799 D7 1.06 0.092 3.63 2.42 2.61 49.404 D8 4.22 0.366 4.56 3.04 5.92 45.978 D10 5.91 0.512 5.36 3.57 7.72 54.715 D11 4.78 0.245 6.16 1.11 4.95 33.290 花岗岩 H1 1.05 0.113 1.59 1.06 1.03 56.430 H2 1.14 0.158 2.21 1.47 2.56 59.995 H3 1.67 0.202 2.69 1.79 3.42 63.938 H5 2.43 0.332 4.56 3.04 5.63 71.631 H6 3.02 0.369 5.36 3.57 6.87 75.020 H8 2.69 0.355 7.69 5.13 8.25 76.612 表 5 RDP-Rc模型验证相关试验参数
Table 5 Parameters for RDP-Rc model verification related tests
试样类型 编号 转速n/(r·min-1) 钻进压强F/MPa 扭矩M/(N·m) 钻进速度v/(mm·min-1) Rc/MPa Rηe/MPa M10 A3 252 4.76 1.184 2.136 3.10 2.66 M15 A9 425 3.74 1.158 1.584 15.68 14.33 M20 A15 425 4.76 1.938 1.998 20.12 22.29 灰岩 D2 252 6.80 3.604 1.302 47.53 43.04 D5 425 5.78 1.643 0.390 88.35 101.46 D9 721 6.80 3.248 2.406 68.73 60.81 花岗岩 H4 425 6.80 5.180 1.068 164.14 143.49 H7 721 7.82 4.520 1.386 186.32 158.14 -
[1] 贾志欣, 汪小刚, 赵宇飞, 等. 岩石钻孔原位测试技术的应用与改进[J]. 岩石力学与工程学报, 2013, 6: 1264-1269. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201306022.htm JIA Zhi-xin, WANG Xiao-gang, ZHAO Yu-fei, et al. Application and improvement of rock borehole in-situ test technology[J]. Chinese Journal of Rock Mechanics and Engineering, 2013, 6: 1264-1269. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201306022.htm
[2] XU T, RANJITH P G, WASANTHA P L P, et al. Influence of the geometry of partially-spanning joints on mechanical properties of rock in uniaxial compression[J]. Engineering Geology, 2013, 167: 134-147. doi: 10.1016/j.enggeo.2013.10.011
[3] TANG X M, XU S, ZHUANG C X, et al. Assessing rock brittleness and fracability from radial variation of elastic wave velocities from borehole acoustic logging: rock brittleness and fracability[J]. Geophysical Prospecting, 2016, 64(4): 958-966. doi: 10.1111/1365-2478.12377
[4] WANG Q, GAO H, YU H, et al. Method for measuring rock mass characteristics and evaluating the grouting-reinforced effect based on digital drilling[J]. Rock Mechanics and Rock Engineering, 2019, 52(3): 841-851. doi: 10.1007/s00603-018-1624-9
[5] KAHRAMAN S, BILGIN N, FERIDUNOGLU C. Dominant rock properties affecting the penetration rate of percussive drills[J]. International Journal of Rock Mechanics & Mining Sciences, 2003, 40(5): 711-723.
[6] GUI M W, SOGA K, BOLTON M D, et al. Instrumented borehole drilling for subsurface investigation[J]. Journal of Geotechnical and Geoenvironmental Engineering, ASCE, 2002, 128(4): 283-291. doi: 10.1061/(ASCE)1090-0241(2002)128:4(283)
[7] SUGAWARA J, YUE Z Q, THAM L G, et al. Weathered rock characterization using drilling parameters[J]. Canadian Geotechnical Journal, 2003, 40(3): 661-668. doi: 10.1139/t03-007
[8] YUE Z Q, LEE C F, LAW K T, et al. Automatic monitoring of rotary-percussive drilling for ground characterization- illustrated by a case example in Hong Kong[J]. International Journal of Rock Mechanics & Mining Sciences, 2004, 41(4): 573-612.
[9] 谭卓英, 岳中琦, 蔡美峰. 风化花岗岩地层旋转钻进中的能量分析[J]. 岩石力学与工程学报, 2007, 26(3): 478-483. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200703005.htm TAN Zhuo-ying, YUE Zhong-qi, CAI Mei-feng. Analysis of energy for rotary drilling in weathered granite formation[J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26(3): 478-483. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200703005.htm
[10] 谭卓英, 王思敬, 蔡美峰. 岩土工程界面识别中的地层判别分类方法研究[J]. 岩石力学与工程学报, 2008, 27(2): 316-322. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200802016.htm TAN Zhuo-ying, WANG Si-jing, CAI Mei-feng. Study on discriminant classification method for ground formation in identification of geotechnical engineering interfaces[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(2): 316-322. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200802016.htm
[11] AALIZAD S A, RASHIDINEJAD F. Prediction of penetration rate of rotary-percussive drilling using artificial neural network a case study[J]. Archives of Mining Sciences, 2012, 57(3): 715-728. doi: 10.2478/v10267-012-0046-x
[12] YAşAR E, RANJITH P G, VIETE D R. An experimental investigation into the drilling and physico-mechanical properties of a rock-like brittle material[J]. Journal of Petroleum Science & Engineering, 2011, 76(3/4): 185-193.
[13] HE M M, LI N, ZHANG Z Q, et al. An empirical method for determining the mechanical properties of jointed rock mass using drilling energy[J]. International Journal of Rock Mechanics & Mining Sciences, 2019, 116(4): 64-74.
[14] HUANG S L, WANG Z W. The mechanics of diamond core drilling of rocks[J]. International Journal of Rock Mechanics and Mining Sciences, 1997, 34(3/4): 1-14.
[15] MOSTOFI M, RASOULI V, MAWULI E. An estimation of rock strength using a drilling performance model: a case study in Blacktip Field, Australia[J]. Rock Mechanics and Rock Engineering, 2011, 44(3): 305-316. doi: 10.1007/s00603-011-0142-9
[16] KALANTARI S, ALIREZA B, HASHEMALHOSSEINI H, et al. 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.
[17] 岳中琦. 钻孔过程监测(DPM)对工程岩体质量评价方法的完善与提升[J]. 岩石力学与工程学报, 2014, 33(10): 1977-1996. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201410005.htm YUE Zhong-qi. Drilling process monitoring for refining and upgranding rock mass quality classific ation 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
[18] 曾俊强, 王玉杰, 曹瑞琅, 等. 基于钻孔过程监测的花岗岩钻进比能研究[J]. 水利水电技术, 2017, 48(4): 112-117. https://www.cnki.com.cn/Article/CJFDTOTAL-SJWJ201704023.htm ZENG Jun-qiang, WANG Yu-jie, CAO Rui-lang, et al. Drilling process monitoring-based study on granite drilling specific energy[J]. Water Resources and Hydropower Engineering, 2017, 48(4): 112-117. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SJWJ201704023.htm
[19] 王玉杰, 赵宇飞, 曹瑞琅, 等. 一种合理评价岩体TBM施工适宜性的试验设备及方法[P]. 中国专利:E21B49, 106593426A, 2017-04-26. [20] ISRM. Suggested methods for determining the uniaxial compressive strength and deformability of rock materials[J]. International Journal of Rock Mechanics & Mining Sciences & Geomechanics Abstracts, 1979, 16: 135-140.
-
期刊类型引用(25)
1. 赵同彬,赵志刚,齐炎山,尹延春,牛旭,李淇凡. 基于大直径钻孔钻进多参量的煤体应力钻测方法. 煤炭科学技术. 2025(01): 122-132 . 百度学术
2. 马语航,何明明,李宁. XCY-2旋切触探系统研发及应用. 岩土力学. 2025(03): 1025-1038 . 百度学术
3. 易文豪,夏覃永,孙鸿强,何永义,王明年. 基于岩石可钻性指标的钻爆法隧道掌子面岩石坚硬程度智能判识. 中国铁道科学. 2025(02): 128-139 . 百度学术
4. 刘勇斌,张晓平,李馨芳,李德宏,陈广军,刘小波,熊雪飞,杨朗,李玉生. 基于冲击回转钻进的岩石强度随钻识别研究. 岩土力学. 2024(03): 857-866 . 百度学术
5. 何海峰. 工程勘察中钻进参数数字化技术应用探讨. 工程勘察. 2024(06): 12-16+23 . 百度学术
6. 郝建,刘河清,刘建康,吕家庆,郑义宁,刘建荣. 基于振动信号的岩石单轴抗压强度钻进预测实验研究. 岩石力学与工程学报. 2024(06): 1406-1424 . 百度学术
7. 盛益明,柴修伟,常志锋,刘建,李治全,向彬. 磷矿地下爆破随钻参数与岩体可爆性关系模型研究. 爆破. 2024(02): 120-126 . 百度学术
8. 吴俊. 基于多要素钻进信息的围岩分级方法研究. 科学技术创新. 2024(20): 157-160 . 百度学术
9. 程谞,汤华,吴振君,张勇慧,秦辉,王腾,贾泽庆,房昱纬. 矿山岩质边坡地层结构及岩体抗剪强度参数随钻智能识别方法. 工程科学与技术. 2024(05): 35-47 . 百度学术
10. 李萍丰,徐振洋,王雪帆,郭润泽,姜立春. 露天矿山爆破共享控制理论的基本方法. 爆破. 2024(03): 16-25+50 . 百度学术
11. 肖浩汉,曹瑞琅,王玉杰,赵宇飞,孙彦鹏. 随钻监测数据预处理方法研究. 水利学报. 2024(11): 1379-1390 . 百度学术
12. 王琦,高红科,冯帆. 岩体力学参数旋切钻进测试方法研究进展. 绿色矿山. 2024(04): 333-343 . 百度学术
13. 岳中文,戴诗清,李杨,岳小磊,李世辉,曹武. 煤巷液压锚杆钻机随钻参数采集系统及其应用. 矿业科学学报. 2023(01): 66-73 . 百度学术
14. 于庆磊,王宇恒,李友,曹永胜,蒲江涌. 基于随钻测量的岩体结构与力学参数表征研究进展. 金属矿山. 2023(05): 45-58 . 百度学术
15. 陈祖军,何明明,周佳佩,邓边员. 基于数字钻技术的岩石强度特性预测方法研究. 长沙理工大学学报(自然科学版). 2023(03): 91-101 . 百度学术
16. 王圆文,丘永富,杨天鸿,邓文学,李华,张紫瑞,刘洋. 基于钻孔数据库和RMR法的露天矿非均匀性岩体质量参数表征及应用. 金属矿山. 2023(06): 192-199 . 百度学术
17. 刘华吉,孙红林,张占荣,尤明龙,谭飞,李炜. 基于随钻参数的砂岩与砂质泥岩地层分界面智能识别. 隧道建设(中英文). 2023(S1): 304-312 . 百度学术
18. 程勇,王琛,刘夏临,刘继国,陈世纪,黄胜. 基于机器学习的隧道地质勘察岩性识别分析及应用研究. 隧道建设(中英文). 2023(09): 1549-1557 . 百度学术
19. 张晟斌,刘继国,刘夏临,舒恒,程勇,赵强,马保松. 一种基于水平定向钻的模型试验装置. 隧道建设(中英文). 2023(09): 1455-1462 . 百度学术
20. 乔梁,李晓昭,邓龙传,熊志勇,张弛. 基于砂浆钻进的随钻机械参数监测试验研究. 工程地质学报. 2023(06): 2020-2029 . 百度学术
21. 王少锋,吴毓萌,蔡鑫,周子龙. 基于随钻参数的岩石强度预测与可钻性识别(英文). Journal of Central South University. 2023(12): 4036-4051 . 百度学术
22. 岳中文,岳小磊,杨仁树,王煦,李为,戴诗清,李杨. 随钻岩性识别技术研究进展. 矿业科学学报. 2022(04): 389-402 . 百度学术
23. Haoteng Wang,Mingming He,Zhiqiang Zhang,Jiwei Zhu. Determination of the constant m_i in the Hoek-Brown criterion of rock based on drilling parameters. International Journal of Mining Science and Technology. 2022(04): 747-759 . 必应学术
24. 罗文超,王军伟. 基于钻进参数的围岩分类指标优化研究. 工程勘察. 2022(10): 19-23+30 . 百度学术
25. 徐振洋,吴怡璇,王雪松,刘鑫,郭连军,柴青平. 旋压钻进破岩的应力特征研究. 金属矿山. 2022(10): 24-29 . 百度学术
其他类型引用(18)