, 徐曾和1,2, 杨鑫1,21. 东北大学 深部金属矿山安全开采教育部重点实验室, 辽宁 沈阳 110819;
2. 东北大学 资源与土木工程学院, 辽宁 沈阳 110819
收稿日期：2017-08-11
基金项目：国家自然科学基金资助项目(U1710253，51574059)。
作者简介：杨斌(1989-)，男，吉林白城人，东北大学博士研究生;
杨天鸿(1968-)，男，辽宁沈阳人，东北大学教授, 博士生导师;
徐曾和(1953-)，男，上海人，东北大学教授, 博士生导师。
摘要：以陕北榆横矿区为例, 研制全程可视化水沙两相高速流动试验装置, 开展颗粒起动试验、溃沙试验和水沙两相流动试验.试验结果表明:细颗粒的流失是溃沙的先决条件, 为粗颗粒的运移创造膨胀空间.0.3~0.6 mm风积沙颗粒之间主要以碰撞的形式传递能量, 而0.15~0.3 mm的风积沙颗粒之间的相互摩擦和挤压作用占优, 相比于碰撞作用能量损失较少, 相同时间内出沙量随着颗粒粒径的增大呈递减趋势.随着轴向应力的不断增大, 风积沙所需要的溃沙临界流速也越来越大.当风积沙粒径d＜0.6 mm时, 溃沙时渗流状态为线性层流, 临界流速为0.03~0.4 cm/s.单位时间内溃沙量与水力梯度成正比, 表明含水层水压力是决定溃沙灾害程度的关键因素.
关键词：厚松散层突水溃沙临界流速轴向应力拖曳力体积膨胀
Critical Velocity of Sand Inrush and Flow Characteristics of Water-Sand in Thick Unconsolidated Formations of Mine in Western China
YANG Bin1,2, YANG Tian-hong1,2
, XU Zeng-he1,2, YANG Xin1,21. Key Laboratory of Ministry of Education on Safe Mining of Deep Metal Mines, Northeastern University, Shenyang 110819, China;
2. School of Resources & Civil Engineering, Northeastern University, Shenyang 110819, China
Corresponding author: YANG Tian-hong, E-mail: yangtianhong@mail.neu.edu.cn
Abstract: Taking the Yuheng mine in Northern Shaanxi as an example, aeolian-sand incipience experiment and water-sand inrush experiment were expansion by high-velocity seepage and entire visual equipment of water-sand two-phase flow. The experiment results indicated that, the loss of fine particles is the precondition of water-sand inrush and creates expansive space for the coarse particle migration. The process of water-sand inrush is a physical process, when the aeolian-sand size between 0.3~0.6 mm, the energy transmits mainly in the form of particles collision.But it is different for the aeolian-sand size between 0.15~0.3 mm, because the friction and squeeze are in the dominant role, that will lose less energy compared with collision. So the variation mass of aeolian-sand inrush is inversely proportional to the particle size. The critical velocity of aeolian-sand inrush is becoming higher with the increasing of axial stress. When the aeolian-sand begins to move, the seepage state is linear laminar flow and the critical velocity is 0.03~0.4 cm/s for the aeolian-sand size d < 0.6 mm.The mass of sand inrush is proportional to the hydraulic gradient in unit time. It can be deduced from this result that water pressure of aquifer is the key factor, which greatly decisive the disaster degree of water-sand inrush.
Key words: thick unconsolidated formationwater-sand inrushcritical velocityaxial stressdrag forcevolumetric expansion
近年来西部矿区成为我国煤炭资源开发的重心.中国西部的浅埋煤层位于含水量大的第四纪冲积层下方, 松散含沙层较厚、沙粒之间胶结程度较差, 采动引起的突水溃沙是对安全生产的主要威胁.2012年11月17日15时, 陕西省神木县隆德煤矿由于打钻诱发突水溃沙事故, 导致大山消失, 地面下塌约25m, 造成上亿元的经济损失[1].由于地下采煤工程的隐蔽性和不可视性, 现场观测突水溃沙机理和影响因素比较困难, 所以通过室内试验模拟水沙两相通过不动多孔介质(粗)骨架的渗流试验, 研究水沙两相渗流机理是积极措施之一.
国内对于突水溃沙问题的研究主要是基于地下水动力学理论以及颗粒流理论, 采用室内试验、数值模拟和理论分析等方法进行研究, 并以临界水力梯度作为判据.水沙通过多孔介质的渗流从更广泛的意义上是液、固两相混合流体通过多孔介质的渗流, 介于液、固两相流动和多相流体流动之间.以往的研究均是将水沙混合物当作单相流体进行研究, 忽略了水和沙在流动过程中的相互作用.
Shirvanian等[2]和Ookawara等[3]利用FLUENT有限元软件研究了管道内水沙的运动规律.Wang和Ni等基于固液两相流理论研究了水平方管内水沙两相流动特征以及颗粒浓度的分布形式[4-6].上述对于水沙两相流的理论研究主要是对河流以及管道中的水沙两相运动规律的研究.由于试验设备的限制, 不能够观察突水溃沙的全程, 不能实现水沙精确的分相测量, 缺少对于突水溃沙临界流速、高水力梯度条件下水沙两相流动特征的研究.
本文以榆横矿区为例, 研制全程可视化水沙两相高速流动试验装置, 研究溃沙临界流速与粒径、应力、渗流状态的关系, 得出该地区水沙两相流动规律、水沙相互作用机理及能量传递机理, 提出突水溃沙灾害的临界流速判据.为煤矿开采过程中突水溃沙灾害的预测与防治提供科学依据.
1 试验装置及测试方法1.1 试验装置本试验采用东北大学自主研发的水沙两相高速流动试验装置, 试验设备及原理图见文献[7].
对于水沙两相混合流动试验, 要研究其运动机理, 必须得到水、沙分相流速, 分相测量装置如图 1所示, 容器1中预先盛满水, 当水沙混合流体进入容器1中, 沙颗粒沉降堆积在容器底部并替换出等体积的水, 水通过容器1的溢流口进入容器2, 由电子天平进行实时测量, 可以得到单位时间内水的质量增量, 被水密度与过流断面的面积除, 就能得到水的流动速度(比流量).电子天平读数的变化反映了进入容器1中风积沙的质量与排出的等体积水质量之差, 由此可以得到沙的运动速度.
图 1(Fig. 1)
 | 图 1 固液两相分相测量装置Fig.1 Solid-liquid two-phase phase separation measuring device |
1.2 试验材料陕北榆横矿区特有的地表风积沙地貌与发达的地表水系是造成高强度、大规模开采条件下发生突水溃沙灾害的主要原因.风积沙的密度为2.62 g/cm3, 通过高频振筛机筛分后可知该地区风积沙主要颗粒粒径为0.15~0.6 mm, 其中0.15~0.3 mm粒径的颗粒质量分数为58.84%, 0.3~0.6 mm粒径的颗粒质量分数为39.84%, 粒径小于0.15 mm颗粒质量分数仅为1.6%.由于试验样本比较少, 所以增加了颗粒形状、表面粗糙度与风积沙相似, 密度为2.65 g/cm3, 粒径为0.6~1.0和1.0~2.0 mm的海沙(图 2).
图 2(Fig. 2)
 | 图 2 试验材料Fig.2 Experimental materials (a)—风积沙；(b)—海沙；(c)—石英砂. |
1.3 试验方案与方法1) 颗粒起动试验(图 3).采用下部填装沙样, 上覆钢珠, 在两端施加轴压的条件下, 测定沙粒起动流速, 即溃沙临界流速.试验装样时为了减小“粮仓效应”[8]对轴压效果的影响, 沙柱的高度控制在100 mm, 小于直径的两倍, 并在有机玻璃装样筒内壁涂抹一层薄薄的凡士林, 减小壁面效应.在沙样与钢珠间添加多孔板, 防止钢珠沉入沙中, 上部钢珠强度大且填装密实, 认为是不可压缩多孔介质, 所以轴压柱塞向下移动的距离完全是沙柱的压缩量.沙粒粒径范围:0.075~0.15 mm, 0.15~0.3 mm，0.3~0.6 mm, 0.6~1 mm, 1~2 mm, 2~2.36 mm共6组.轴向应力设置为:0, 0.1, 0.5, 1.0, 1.5, 1.8 MPa共6级.
图 3(Fig. 3)
 | 图 3 风积沙起动试验Fig.3 Experiment of aeolian-sand incipience |
2) 溃沙试验(图 4).预先将粒径为2.36~4.75 mm的石英砂填入长为200 mm的有机玻璃管中作为骨架, 其孔隙率为0.46.然后再将粒径为0.075~0.6 mm的风积沙填入长为322 mm的有机玻璃管中, 沙样的孔隙率为0.34.用石英砂骨架模拟破碎岩体突水通道, 风积沙模拟含沙含水层, 石英砂和风积沙之间用多孔板分隔.研究其在不同水力梯度下的出沙量, 实测水力梯度为1.02, 2.04, 3.06, 4.08, 5.10.
图 4(Fig. 4)
 | 图 4 溃沙试验Fig.4 Experiment of sand inrush |
3) 水沙两相流动试验(图 5).在长为322 mm的有机玻璃管中, 分别装填粒径为0.15~0.3和0.3~0.6 mm的风积沙作两次试验.采用由下进水、由上出水的方式研究水沙两相流动机理.具体试验方案见表 1.
图 5(Fig. 5)
 | 图 5 水沙两相流动试验Fig.5 Experiment of two-phase flow constituted by water and sand |
表 1(Table 1)
 表 1 试验方案设计Table 1 Experimental scheme
| 表 1 试验方案设计 Table 1 Experimental scheme |
每次试样装填完毕后, 抽真空20 min, 然后打开进水截止阀, 采用由底部向上逐层饱和的方法, 令水缓慢地进入试样, 在进水过程中保持真空泵运转以维持真空度不变, 以保证颗粒间的空气和颗粒表面的气泡完全被水置换, 颗粒表面被水完全浸润.试验时水温基本恒定在20 ℃左右, 流体动力黏滞系数μ取1.01×10-3 Pa·s.
2 试验结果及分析2.1 溃沙临界流速与粒径、应力、渗流状态的关系在没有发生突水时, 承压含水层风积沙在围压与自重作用下紧密压实, 此时风积沙颗粒堆具有静态性质, 风积沙颗粒在重力、上层沙颗粒的压力、下层沙颗粒的支持力、静摩擦力等力的作用下保持现有的形态, 能够抵抗一定的应力而不发生屈服和变形.突水发生时, 还要受到水的浮力、拖曳力、渗透压力以及动摩擦力等力的作用.拖曳力是溃沙的主要动力, 随着突水速度逐渐增大, 拖曳力也越来越大, 当拖曳力增加到一定程度后, 破坏了颗粒之间通过力链(或拱)形成的力的网络, 颗粒开始运动.风积沙的粒径非常小, 所以颗粒起动速度即为溃沙临界速度.试验中不断增加试样上游水压力, 水的流速随之增加, 通过有机玻璃管可以观察到沙粒群首先发生体积膨胀, 其次是颗粒在水流拖曳力的作用下产生翻滚运动, 逐渐进入上覆钢珠堆积的多孔介质中, 此时的流速视为溃沙的临界流速.
由图 6可知, 对于同一种粒径, 随着轴向应力不断增大, 风积沙颗粒之间镶嵌咬合越紧密, 阻止风积沙颗粒运动的关键力即摩擦阻力逐步增大, 风积沙所需要的起动流速也越来越大.值得注意的是0~0.075mm风积沙颗粒的起动速度随着轴压的增大先不变再减小最后趋于稳定, 起动流速在0.33~0.029 mm/s之间, 这是由于该粒径范围的风积沙非常细小, 饱水堆积状态下非常密实, 沙柱整体压缩性差, 很小的扰动就会使其悬浮在水中.
图 6(Fig. 6)
 | 图 6 不同粒径范围沙的轴向应力与溃沙临界流速关系Fig.6 Relationship between the axial stress and velocity of aeolian-sand inrush for different particles size range (a)—0.075~0.6 mm; (b)—0.6~2.36 mm. |
对于相同轴压条件下, 从小颗粒风积沙到大颗粒海沙, 随着颗粒自重应力的增加, 起动所需要的拖曳力也越来越大, 起动流速逐步增大.从雷诺数可以看出, 当沙粒径d＜0.6 mm时, 颗粒运动时渗流状态为线性层流, 临界流速为0.03~0.4 cm/s, 临界雷诺数为0.29~1.8;当2.36 mm>d>0.6 mm时, 颗粒运动时渗流状态符合紊流, 临界流速为3~11 cm/s, 此时临界雷诺数变化范围比较大, 为20~240.
2.2 含水层水压力对溃沙的影响溃沙试验中通过有机玻璃管可以观察到, 当渗流速度很小时(大约0.3 mm/s), 水流就开始携带着粒径为0~0.075 mm的风积沙流经石英砂骨架突水通道到达出口处, 导致出口处水流变浑浊, 随着水流中含沙比越来越大, 风积沙的孔隙率逐渐增大.陆坤权等认为局部颗粒堆积密度的降低即体积膨胀是颗粒流动的前提[9], 0~0.075 mm粒径范围风积沙的流出, 为大颗粒风积沙运动提供了膨胀空间.这是一个由量变到质变的物理过程.可见, 细颗粒的流失是溃沙的先决条件, 且细颗粒的含量影响着溃沙量和溃沙范围.现场可以通过观察巷道渗水的浑浊程度, 提前采取相应的防治措施.
图 7为风积沙在不同压力梯度下出沙量随时间变化曲线, 由图可知在同一压力梯度条件下, 出沙量与时间呈线性关系, 可见溃沙过程中很快就达到了水沙两相稳定流动.水力梯度为1.02和2.04时90 s的出沙量占总量的20%~40%, 而水力梯度为4.08和5.10时出沙量占总量的70%~100%, 这是一个阶跃性变化, 且单位时间内溃沙量与水力梯度成正比, 水力梯度越大, 溃沙量越多, 表明含水层水压力是决定溃沙灾害程度的关键因素, 含沙量大、速度快是高水压溃沙的特点.初始水力梯度越大, 溃沙发生越迅速, 带来的灾难也越严重, 因此实际工程中含水层水压力可以作为实时监测和预报突水溃沙灾害的前兆信息.
图 7(Fig. 7)
 | 图 7 不同水力梯度条件下溃沙量变化规律Fig.7 Variation mass of sand inrush with different hydraulic gradients |
2.3 水沙两相流动机理图 8中两组试验均是在水力梯度为3.57条件下进行的.粒径范围为0.15~0.3 mm的风积沙溃沙试验中, 沙的流速基本大于水的流速, 表明沙对水的流动具有促进作用, 而0.3~0.6 mm的试验结果正好与0.15~0.3 mm相反, 水的流动速度比沙的速度要大很多, 沙对水的流动起到阻碍作用.从颗粒流的角度解释这种现象可以认为0.3~0.6 mm风积沙颗粒之间主要以碰撞的形式传递能量, 正是这种非弹性碰撞导致颗粒速度的突变, 其本质是稀疏流和密集流交替出现, 而沙粒间隙中水流的存在对沙粒碰撞起着阻尼作用.Wang等认为颗粒群之间的相互作用除了碰撞以外, 还有颗粒之间的相对滑动及相互挤压作用[5].0.15~0.3 mm风积沙溃沙试验中, 沙柱在上游水压作用下整体被托起来了(图 5), 颗粒之间形成了一种类似于“液桥力”[10], 颗粒之间的相互摩擦和挤压作用占优, 相比于碰撞作用能量损失较少, 其本质是持续的密集流, 所以相同时间内出沙量与颗粒粒径呈反比的趋势(图 9).
图 8(Fig. 8)
 | 图 8 风积沙和水两相流动速度随时间变化曲线Fig.8 Curves of aeolian-sand and water flow rate changing with time (a)—0.15~0.3 mm; (b)—0.3~0.6 mm. |
图 9(Fig. 9)
 | 图 9 不同粒径风积沙出沙量随时间变化曲线Fig.9 Curves of mass of aeolian-sand inrush with time for different particle size |
3 结论1) 当粒径大于0.075 mm时, 随着轴向应力的不断增大, 风积沙溃沙临界流速也越来越大.对于粒径小于0.075 mm的风积沙, 轴压对其溃沙临界流速的影响可以忽略.
2) 相同轴压条件下, 从小颗粒风积沙到大颗粒海沙, 溃沙临界流速逐步增大.当沙粒径d＜0.6 mm时, 发生突水溃沙瞬时渗流状态为线性层流, 临界流速为0.03~0.4 cm/s, 临界雷诺数为0.29~1.8;当2.36>d>0.6 mm时, 颗粒运移时渗流状态符合紊流, 临界流速为3~11 cm/s, 此时临界雷诺数变化范围比较大, 为20~240.
3) 突水溃沙灾害发生时, 粒径为0.3~0.6 mm风积沙颗粒之间主要以碰撞的形式传递能量, 而0.15~0.3 mm的风积沙颗粒之间形成了一种类似于“液桥力”, 颗粒之间的相互摩擦和挤压作用占优, 相比于碰撞作用能量损失较少, 相同时间内出沙量随着颗粒粒径的增大而减少.
参考文献
| [1] | 张蓓, 张桂民, 张凯, 等. 钻孔导致突水溃沙事故机理及防治对策研究[J].采矿与安全工程学报, 2015, 32(2): 219–225. ( Zhang Bei, Zhang Gui-min, Zhang Kai, et al. Water and sands bursting mechanism induced by geological borehole and control measures[J].Journal of Mine & Safety Engineering, 2015, 32(2): 219–225.) |
| [2] | Shirvanian P A, Calo J M, Hradil G. Numerical simulation of fluid-particle hydrodynamics in a rectangular spouted vessel[J].International Journal of Multiphase Flow, 2006, 32(6): 739–753.DOI:10.1016/j.ijmultiphaseflow.2006.02.009 |
| [3] | Ookawara S, Street D, Ogawa K. Numerical study on development of particle concentration profiles in a curved microchannel[J].Chemical Engineering Science, 2006, 61(11): 3714–3724.DOI:10.1016/j.ces.2006.01.016 |
| [4] | Wang G Q, Ni J R. The kinetic theory for dilute solid-liquid two-phase flow[J].International Journal of Multiphase Flow, 1991, 17(2): 273–281.DOI:10.1016/0301-9322(91)90020-4 |
| [5] | Wang G Q, Ni J R. Kinetic theory for particle concentration distribution in two-phase flows[J].Journal of Engineering Mechanics, 1990, 116(12): 2738–2748.DOI:10.1061/(ASCE)0733-9399(1990)116:12(2738) |
| [6] | Ni J R, Wang G Q, Borthwick A G L. The kinetic theory for vertical profile of paticle concentration in dilute and dense solid-liquid flows[J].Journal of Hydraulic Engineering, 2000, 126(12): 893–903.DOI:10.1061/(ASCE)0733-9429(2000)126:12(893) |
| [7] | Yang B, Yang T H, Xu Z H, et al. Numerical simulation of the free surface and water inflow of a slope, considering the nonlinear flow properties of gravel layers:a case study[J].Royal Society Open Science, 2018, 5(2): 172109.DOI:10.1098/rsos.172109 |
| [8] | Hidalgo R C, Lozano C, Zuriguel I, et al. Force analysis of clogging arches in a silo[J].Granular Matter, 2013, 15(6): 1–8. |
| [9] | 陆坤权, 刘寄星. 颗粒物质(上)[J].物理, 2004, 33(9): 629–634. ( Lu Kun-quan, Liu Ji-xing. Static and dynamic properties of granular matter(Ⅰ)[J].Physics, 2004, 33(9): 629–634.DOI:10.3321/j.issn:0379-4148.2004.09.002) |
| [10] | Washino K, Chan E L, Matsumoto T, et al. Normal viscous force of pendular liquid bridge between two relatively moving particles[J].Journal of Colloid and Interface Science, 2017, 494: 255–265.DOI:10.1016/j.jcis.2017.01.088 |
