Mechanical characteristics of high-pressure water jets impinging on coal
-
摘要: 以质量守恒与动量守恒定律为基础, 建立了高压水射流冲击煤体的力学模型。运用此模型分析了高压水射流在冲击煤体的过程中, 未破水体、破碎水体、煤体的破碎区与扩孔区的力学特征, 利用严格的力学守恒关系得出高压水射流冲击煤体的简化常微分方程组。将理论计算结果与现场实验和数值模拟结果进行对比, 结果表明:理论计算结果与数值模拟结果和实验结果基本一致。此模型具有明确的力学意义, 且能够反映真实的冲击过程。Abstract: On the basis of conservation of mass and momentum, a mechanical model was established for high-pressure water jets impinging on coal.By using this model, the mechanical characteristics were analyzed for the intact water jet, cracked water jet, and the crushing zone and crater expansion zone of the coal.The system of ordinary differential equations was obtained to describe the process of high-pressure water jets impinging on coal.And the theoretical values were compared with the numerical simulation and the experimental results.The theoretical values are consistent with the numerical simulation and the experimental results.So the established mechanical model can reflect the actual process of high-pressure water jets impinging on coal.
-
Key words:
- fluid mechanics /
- water jet /
- high pressure /
- coal /
- impinging
-
高压水射流破碎煤体是高压水射流的一项具体应用, 提高水射流冲击破碎煤体效率在水力采煤和水力冲孔治理瓦斯等领域具有重要的研究价值。若要提高水射流破碎煤岩的效率, 必须首先在理论上掌握煤体在高压水射流作用下的力学特征。因此, 在高压水射流破煤过程中, 水体与煤体的力学特征一直是高压水射流破煤研究的重点。倪红坚等[1-2]、王瑞和等[3]、廖华林等[4-5]、卢义玉等[6]、田方宝等[7]对高压水射流冲击岩石的基本力学特性进行了深入的研究, 在高压水射流冲击岩石的破孔过程、岩石表面与内部的应力分布、岩石破碎的门槛压力等方面得出一系列有价值的研究成果, 为高压水射流破岩提供了理论与实验方面有益的参考。穆朝民等对于煤体在高压水射流作用下的动态损伤机理[8]、临界破煤强度[9]及高压磨料射流破煤体的过程和数值计算方法[10]进行了探讨, 得出了在高压水射流作用下煤体的损伤形式和临界破煤压力。
目前高压水射流破岩机理大多为实验和数值分析结果, 没有涉及高压水射流破煤的力学特征, 即尚未建立高压水射流在冲击煤体的过程中, 未破水体、破碎水体、煤体的破碎区与扩孔区完整的力学方程, 因此对于高压水射流破煤的力学机理很难形成有效的指导。本文中, 拟在李永池等[11]关于长杆弹高速侵彻混凝土相关研究(主要是弹体蘑菇头系数和混凝土刚性破碎流体介质假设)的基础上得出高压水射流冲击煤体的基本力学特征。
1. 基本力学分析与假设
高压水射流冲击煤体的力学分析如图 1所示[11], A0-A0以左为高压水射流未变形破裂部分(未变形区), 面积为SA0, 长度为l; A0A0A1A1为高压水射流的变形蘑菇头区, A1-A1为高压水射流破碎前阵面, 面积为SA1, 未变形区与变形蘑菇头区合称为未破碎高压水射流区; A1A1B1A2A2B1A1为高压水射流的破碎反射水射流区, 其中B1-A2为破碎高压水扩孔终止界面, A1-B1为高压水反射界面, 高压水形成的总体环形面积为SB1, A2-A2为反射高压水前沿和煤渣后沿的交界面, 即冲击交界面, 其面积为SA2; A2B2A3A3B2A2为煤的破碎和扩孔区, A3-A3为煤的破碎前阵面, 面积为SA3, B2A3为煤渣与实体煤的交界面, 即扩孔界面其总体环形面积为S, 图中KB2A2B1称为煤渣的反向运动区, 边界A2B2上的煤渣是和冲击界面A2-A2一起以冲击速度u运动的, 并且A2-A2上的压力为冲击压力p。KB2为煤渣最终成孔的截面。
高压水射流冲击煤体很复杂, 为简化问题, 假定[11]:(1)煤体为刚性破碎流体介质, 当p(冲击压力) < Rt(煤体破碎强度)时, 煤体为刚体; 当p≥Rt时, 煤破碎区与实体煤交界面
上的压力为Rt。(2)高压水射流未变形部分与变形蘑菇头区以共同速度v运动。(3)高压水射流蘑菇头区、高压水射流破裂反射水射流区、煤破碎和扩孔区都极薄, 因此蘑菇头区质量和动量很小。以上假设适用于形状完整较坚硬的煤体, 对松散煤体不适用, 当水射流冲击压力小于煤体破碎强度时不发生破煤现象。
未破碎高压水射流(未变形区和蘑菇头区)具有共同的质点速度V=v(瞬时高压水射流速度), 高压水射流破碎前阵面A1-A1上的速度v下降到冲击界面
上的速度
, 压应力则由A1-A1上的0提高到冲击界面A2-A2的p, 即
煤破碎区中质点速度由A2-A2上的u降到煤破碎界面A3-A3处的0, V∈ (u, 0), 而压应力则由A2-A2上的p降到A3-A3处的Rt, σ∈ (p, Rt)。高压水射流未破碎区、破碎区、煤破碎区中质点速度变化分别为
。水体从A1-B1相对于与冲击界面一起运动的坐标系向左的反射水射流的速度为uf, 则在绝对坐标系中水体向左喷出的速度为uf-u。
2. 基本方程组
2.1 未破碎高压水射流的质量与动量守恒
在忽略蘑菇头区质量时, 未破碎高压水射流的质量守恒条件, 即高压水射流的消蚀方程为:
dιdt=−(v−u),dh dt=u (1) 式中:l(t)、v(t)、u(t)和h(t)分别为高压水射流的长度(忽略蘑菇头区的质量, 此时水射流近似等于未破碎水射流的长度)、速度、冲击界面速度和冲击深度。
将坐标系建立在冲击界面A2-A2上, 建立动量守恒方程为:
ρpSA0dl dt(v−u)+ρpSA0lddt(v−u)+Md(v−u)dt+(v−u)dM dt=−ρpSA1(v−u)2−(ρpSA0l+M)du dt (2) 式中:M为蘑菇头区的质量,
为高压水射流反射流出开口体系引起的动量损失, ρp为水射流的密度,
为冲击坐标系引起的惯性力。忽略蘑菇头区的质量和动量, 则式(2)可以转化为:
dv dt=−(v−u)2l(γ−1) (3) 式中:
为蘑菇头面积因数; 当γ =1时, 式(3)与Tate[11]给出的结果一致。
2.2 破碎高压水射流的质量与动量守恒
dM1 dt=ρpSA1(v−u)−ρpSB1uf (4) 式中:M1为高压水射流的破碎和反射水射流区的质量, SB1为反射水射流界面A1B1的面积, uf为反射压水射流相对于冲击坐标系的平均反射水射流速度; ρpSA1(v -u)为单位时间内通过A1A1面流入体系高压水的质量, ρpSB1uf为单位时间内通过A1B1面流出体系反射水射流的质量。
因为 dM1 dt≈0, 式(4)可以转化为:
uf=(v−u)/β (5) 式中
为水射流反射水射流区面积与蘑菇头面积之比, 称为反射水射流面积因数。
高压水射流破碎和反射水射流区动量守恒方程为:
M1d(u1−u)dt+(u1−u)dM1 dt=−(SA1+SB1)p+ρpSA1(v−u)2+ρpSB1uf(−uf)−M1du dt (6) 方程(6)最后4项分别A2-A2面上的外力、水射流通过A1-A1面流入的动量、水射流通过A1-B1面流出的动量、惯性力。根据式(4)~(5), 式(6)可以转化为:
p=1βρp(v−u)2 (7) 2.3 煤体破碎和扩孔区的质量与动量守恒
以M2为煤的破碎和扩孔区的质量, B2-A3面上破碎粒相对于冲击坐标系的速度为u, 煤破碎区质量守恒方程为:
dM2 dt=ρtSA3u−∫Su⋅nds (8) 式中:
为破碎煤介质通过A3-A3截面流入的质量, ρt为破碎煤介质密度
为通过截面B2-A3流出的质量。因为
式(8)可以转化为:
∫Su⋅nds=ρtSA3u (9) 煤破碎区质量守恒方程为:
M2ddt(u2−u)+(u2−u)dM2 dt=SB2p−Rt(SA3+Scosθ)−ρtSA3u2−(u2−u)∫Su⋅nds−M2du dt (10) 式中:
为截面B2-B2的面积,
为通过截面A3-A3流入体系的动量, (u2−u)∫Su⋅nds为通过截面B2-A3流出体系的动量
为惯性力, θ为B2-A3面与冲击轴线的夹角。显然有:
−Rt(SA3+Scosθ)=−RtSB2 (11) 将方程(9)和(11)代入(10), 可将方程(10)化为:
p=Rt+αρtu2u (12) 式中
为煤体的扩孔因数。煤体破碎和扩孔区质点平均速度u2近似为冲击速度的一半[11]:
u2=u/2 (13) 将式(13)代入式(12), 则式(12)可以转化为:
p=Rt+α2ρtu2 (14) 结合式(7)、(14), 可得:
1βρp(v−u)2=p=Rt+α2ρtu2 (15) 可见发生高压水射流破煤介质的最小临界冲击压力是p=Rt。当撞击速度足够高, 既产生水射流的破坏又产生对靶板的冲击时, 由式(15)可得冲击速度u和水射流速度v间的关系:
u=u(v)≡v−μ(v2+λ)121−μ2 (16) μ≡(αρt2ρp/β)12,λ≡2Rt(1−μ2)αρt (17) 根据穆朝民等[9]关于破碎强度参数Rt的确定, 可得:
Rt=Y[E/(3Y)1−√σtY1−μ√2]2η3 (18) 式中:E为弹性模量, σt为抗拉强度, Y为抗压强度, μ为泊松比
为压剪因数。
2.4 高压水射流冲击煤体的力学分析
在同时发生水的破碎和煤的冲击破坏的一般情况下, 联立式(1)、(3)、(16)、(18), 即高压水射流冲击煤体问题的常微分方程组如下:
{v(0)=v0,l(0)=l0,h(0)=0t=0dv dt=−(γ−1)l[v−u(v)]2=−f(v)ldl dt=−(v−u)dh dt=u (19) 由式(15)知, 当冲击速度u > 0, 则p(冲击界面压力) > Rt(煤的破坏强度), 高压水射流可破煤; 当冲击速度u=0, 水射流终止破煤, 此时, p≤Rt。
当u=0, 得:
v=√βRtρp≡vc (20) 高压水射流冲击煤体过程如下:
(1) 当v0≤vc, 则不会产生高压水射流破煤效应, 其间l~v之间的关系可令式(19)中u=0而得到:
l=l0exp{∫vv0vf(v)dv} (21) (2) 如果v0 > vc, 则产生高压水射流破煤效应, 总冲击深度可令式(19)中v=vc得出:
h0=−∫vcv0u(v)l(v)f(v)dv=h(vc) (22) 此时高压水射流的剩余长度lc可由式(21)给出:
lc=l(vc)=l0exp[∫vcv0v−u(v)f(v)dv] (23) 以后便只发生水射流的撞击破碎, 其间水射流的瞬时剩余长度l和其瞬时速度v间的关系可将方程(23)中的l0和v0分别代之以lc和vc而得出, 而水射流的最终残余长度lf则为水射流速度下降为0时水射流长度:
l=lcexp[∫vvcvf(v)dv],lf=lcexp[∫0vcvf(v)dv] (24) 3. 高压水射流冲击煤体的数值模拟
3.1 模型和材料参数
运用固流耦合的方法对煤体在高压水射流作用下的力学特征进行数值分析, 水射流冲击煤体的具体尺寸如图 2所示[8]。煤体采用含损伤J-H-C本构模型, 煤体力学参数见表 1, 其中ρ为密度, G为剪切模量, E为杨氏模量, ν为泊松比, A、B、C、N、D1、D2、K1、K2为材料常数, fc为单轴压缩强调, T为极限拉伸静水压力, εf, min为损伤常数, Smax为量刚一强度, pcrush破碎静水压力, μloc为最大体应变, μcrush为破碎体应变; plock为最大静水压力。对水射流采用状态方程:
表 1 煤体材料参数Table 1. Material parameters for coalρ0/(g·cm-3) E/GPa v G/GPa A B C N fc/GPa T/GPa εfmin 1.2 1.2 0.3 1.40 0.78 1.61 0.007 0.62 0.049 0.003 0.01 Smax pcrush/GPa μcrush plock/GPa μlock D1 D2 K1/GPa K2/GPa K3/GPa 7.1 0.001 0.001 0.002 0.11 0.041 1.0 85.0 -171.5 208.5 p=(p0+B)(ρρ∗)k1−B (25) 式中:p为压力, ρ为密度, B、p0、ρ*、k1为常数, B=30.5 GPa, p0=103.3 kPa, k1=7.147, ρ*=1 t/m3。
3.2 数值模拟结果
在30 MPa高压水射流作用下, 7.5、24.3、29.2和36.7μs时的破煤深度分别为0.7、2.2、3.5和4.0 mm, 水射流的剩余长度分别为7.48、4.34、3.26、1.60 mm, 如图 3所示[8]。将水射流出口水压30 MPa换算成水射流冲击煤体的初速度, 并和表 1所列的煤体参数一起代入式(22)~(23), 得出7.5、24.3、29.2、36.7μs时高压水射流的破煤深度分别为0.74、2.80、4.20和5.10 mm, 水射流的剩余长度分别为7.41、4.17、2.93、2.15 mm。这与数值计算结果基本符合, 反映了本次理论计算的正确性。
20、30、35 MPa高压水射流作用下的破煤深度如图 4所示[8]。高压水射流出口压力分别为20、30、35 MPa时的破煤深度分别为2.1、4.0和4.2 mm, 水射流剩余长度分别为2.1、1.6和1.4 mm。
将水射流出口水压换算成水射流冲击煤体的初速度, 并和表 1所列的煤体参数一起代入式(22)~ (23), 得出高压水射流在出口压力为20、30、35 MPa时的破煤深度分别为3.30、5.10和5.17 mm, 剩余长度分别为2.96、2.15、1.93 mm。这与数值计算结果基本符合。
比较理论分析结果和数值模拟结果可以看出:当水射流出口压力越大(水射流初速度越大), 理论计算结果与数值模拟结果的误差越小。这主要是由于水射流出口压力越大, 高压水射流蘑菇头区、高压水射流破裂反射水射流区、煤破碎和扩孔区越薄, 具有的质量和动量越小, 越接近本次理论计算的假设。由图 3可知:高压水射流在冲击煤体的过程中可以分成高压水射流未变形破裂部分(未变形区)和高压水射流的破碎反射水射流区, 这与理论假设基本符合。
4. 高压水射流冲击煤体的现场实验
实验地点选择在淮南矿业集团潘三煤矿1792(3)底板巷联巷。实验区域平均煤厚3.94 m, 煤层倾角5°~9°。2010年4月进行了4次高压水射流冲击煤体的实验, 高压水射流出口压力为30 MPa, 潘三矿13-1煤体力学参数见表 1, 现场水射流冲击煤体实验情况见表 2, 高压水射流对煤体冲击入射角为θ。由于现场对煤体进行高压水射流水力扩孔实验, 高压水射流可以反复冲击煤体, 因此理论计算时将水射流重复冲击煤体的深度进行累积计算, 对比现场和理论计算数值可以看出, 理论计算结果与现场实验基本一致。
表 2 不同冲击角度下水射流的破煤深度Table 2. Penetration depths of coal by water jet at different impact angles爆破孔
编号θ/(°) h/mm 理论 实验 1 45 280 297.3 2 45 260 297.3 3 20 290 312.5 4 45 285 297.3 5. 结论
建立了高压水射流冲击煤体的力学模型, 分析了高压水射流的变形蘑菇头区和破碎反射水射流区、煤的破碎和扩孔区对高压水射流冲击煤体的影响。引入蘑菇头面积因数、反射水射流面积因数、煤体的扩孔因数, 利用力学守恒关系导出了高压水射流冲击煤体的简化常微分方程组。理论计算结果与数值模拟和现场实验结果吻合较好, 可见建立的力学模型能够较好地反映高压水射流冲击煤体的相关规律。
-
表 1 煤体材料参数
Table 1. Material parameters for coal
ρ0/(g·cm-3) E/GPa v G/GPa A B C N fc/GPa T/GPa εfmin 1.2 1.2 0.3 1.40 0.78 1.61 0.007 0.62 0.049 0.003 0.01 Smax pcrush/GPa μcrush plock/GPa μlock D1 D2 K1/GPa K2/GPa K3/GPa 7.1 0.001 0.001 0.002 0.11 0.041 1.0 85.0 -171.5 208.5 表 2 不同冲击角度下水射流的破煤深度
Table 2. Penetration depths of coal by water jet at different impact angles
爆破孔
编号θ/(°) h/mm 理论 实验 1 45 280 297.3 2 45 260 297.3 3 20 290 312.5 4 45 285 297.3 -
[1] 倪红坚, 王瑞和, 张延庆.高压水射流作用下岩石破碎机理及过程的数值模拟研究[J].应用数学和力学, 2005, 26(12): 1445-1452. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=yysxhlx200512008Ni Hong-jian, Wang Rui-he, Zhang Yan-qing. Numerical simulation study on rock breaking mechanism and process under high pressure water jet[J]. Applied Mathematics and Mechanics, 2005, 26(12): 1445-1452. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=yysxhlx200512008 [2] 倪红坚, 王瑞和.高压水射流射孔过程及机理研究[J].岩土力学, 2004, 25(增): 29-32. http://d.wanfangdata.com.cn/Periodical/ytlx2004z1007Ni Hong-jian, Wang Rui-he. Study on progress and mechanism of high pressure water jet perforation[J]. Rock and Soil Mechanics, 2004, 25(Suppl): 29-32. http://d.wanfangdata.com.cn/Periodical/ytlx2004z1007 [3] 王瑞和, 倪红坚.高压水射流破岩机理研究[J].石油大学学报:自然科学版, 2002, 26(4): 118-121. http://www.cnki.com.cn/Article/CJFDTotal-SYDX200204034.htmWang Rui-he, Ni Hong-jian. Study on mechanism breaking rock with high pressure water jet[J]. Journal of Uni-versity of Petroleum, China: Natural Science Edition, 2002, 26(4): 118-121. http://www.cnki.com.cn/Article/CJFDTotal-SYDX200204034.htm [4] 廖华林, 李根生, 牛继磊.淹没条件下超高压水射流破岩影响因素与机制分析[J].岩石力学与工程学报, 2008, 27(6): 523-528. http://www.cqvip.com/Main/Detail.aspx?id=27489390Liao Hua-lin, Li Gen-sheng, Niu Ji-lei. Influential factors and mechanism analysis of rock breakage by ultra-high pressure water jet under submerged condition[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(6): 523-528. http://www.cqvip.com/Main/Detail.aspx?id=27489390 [5] 廖华林, 李根生.超高压水射流冲击岩石的流固耦合分析[J].水动力学研究与进展, 2004, 19(4): 452-457. http://d.wanfangdata.com.cn/Periodical/sdlxyjyjz200404007Liao Hua-lin, Li Gen-sheng. Fluid-structure interaction of high pressure water jets impinging on rock[J]. Journal of Hydrodynamics, 2004, 19(4): 452-457. http://d.wanfangdata.com.cn/Periodical/sdlxyjyjz200404007 [6] 卢义玉, 冯欣艳, 李晓红, 等.高压空化水射流破碎岩石的试验分析[J].重庆大学学报:自然科学版, 2006, 29(5): 87-91. http://d.wanfangdata.com.cn/Periodical_cqdxxb200605023.aspxLu Yi-yu, Feng Xin-yan, Li Xiao-hong, et al. Experiments on breaking rock with high-pressure cavitation water jets[J]. Journal of Chongqing University: Natural Science Edition, 2006, 29(5): 87-91. http://d.wanfangdata.com.cn/Periodical_cqdxxb200605023.aspx [7] 田方宝, 林缅.水射流辅助破岩机理研究(1):气泡空蚀[J].力学与实践, 2007, 29(1): 29-33. http://www.cnki.com.cn/Article/CJFDTotal-LXYS200701006.htmTian Fang-bao, Lin Mian. Study on mechanism of rock breaking with ressure water jet(1): Cavitation bubble[J]. Mechanics and Practice, 2007, 29(1): 29-33. http://www.cnki.com.cn/Article/CJFDTotal-LXYS200701006.htm [8] 穆朝民, 王海露.煤体在高压水射流作用下的损伤机制[J].岩土力学, 2013, 34(5): 1515-1520. http://www.cqvip.com/QK/94551X/20135/45676029.htmlMu Chao-min, Wang Hai-lu. Damage mechanism of coal under high pressure water jetting[J]. Rock and Soils Mechanics, 2013, 34(5): 1515-1520. http://www.cqvip.com/QK/94551X/20135/45676029.html [9] 穆朝民, 吴阳阳.高压水射流冲击下煤体破碎强度的确定[J].应用力学学报, 2013, 30(3): 451-456. http://d.wanfangdata.com.cn/Periodical/yylxxb201303028Mu Chao-min, Wu Yang-yang. Crushing strength of the coal against high pressure water penetration[J]. Chinese Journal of Applied Mechanics, 2013, 30(3): 451-456. http://d.wanfangdata.com.cn/Periodical/yylxxb201303028 [10] 穆朝民, 戎立帆.磨料射流冲击岩石损伤机制的数值分析[J].岩土力学, 2014, 35(5): 1475-1481. http://www.cnki.com.cn/Article/CJFDTotal-YTLX201405038.htmMu Chao-min, Rong Li-fan. Numerical simulation of damage mechanism of abrasive water jet impaction on rock[J]. Rock and Soils Mechanics, 2014, 35(5): 1475-1481. http://www.cnki.com.cn/Article/CJFDTotal-YTLX201405038.htm [11] 李永池, 于少娟, 罗春涛, 等.高速侵彻力学中改进的Tate工程分析方法[C]∥材料和结构的动态响应.合肥: 中国科学技术大学出版社, 2005: 132-139. 期刊类型引用(7)
1. 毕延森,高德利,鲜保安. 煤层气水平井筛管密集分段径向射流增透方法. 石油钻采工艺. 2024(01): 119-129 . 百度学术
2. 张强,张赫哲,田莹,王聪. 截齿辅助冲击作用下坚硬煤体的破碎特性. 煤炭学报. 2022(02): 1002-1016 . 百度学术
3. 张强,王聪,刘玉果,田莹,王亚军,吕祥锋,陈洪月,马英. 难采煤岩的高效破碎方法研究. 煤炭科学技术. 2021(02): 163-176 . 百度学术
4. 梁春苗,姚宁平,姚亚峰,张金宝,沙翠翠,李坤. 自进式动平衡围压水射流钻孔试验装置设计与试验. 煤炭学报. 2020(10): 3533-3540 . 百度学术
5. 陶云奇,周正涛,徐冬冬. 冲孔卸压抽采技术在中马村矿突出煤层中的应用. 煤矿安全. 2018(04): 65-67+71 . 百度学术
6. 许彦鹏,和林朋,王梦乐. 登封矿区低透气性松软煤层水力冲孔措施研究. 煤炭技术. 2016(07): 147-150 . 百度学术
7. 叶林征,祝锡晶,王建青,郭策. 基于SPH-FEM的超声珩磨空化微射流冲击研究. 振动与冲击. 2016(13): 72-77 . 百度学术
其他类型引用(11)
-