Research on the development law of karst caves on water conducting fractures under the influence of mining in Southwest Karst Mining Areas
-
摘要:
黔西南矿区为我国典型的岩溶发育矿区，采动影响下岩溶顶板导水裂隙高度异常发育，雨季大气降水极易通过超高导水通道进入井下工作面，造成工作面涌水灾害，严重影响矿井的正常生产。为此，以新田煤矿煤层开采的导水裂隙为研究对象，在分析矿区溶洞赋存特征基础上，采用现场实测、室内模拟及理论分析等手段，研究了工作面开采过程中溶洞下顶板导水裂隙带的发育规律，揭示了超高导水裂隙的发育机理，结果表明：①研究区域顶板溶洞在地表—玉龙山段—长兴组等层位由上而下具有明显的分带现象，地表沿沟谷处落水洞发育，玉龙山段强含水层中上部溶洞发育，且呈现不同形态的串珠状，长兴组仅局部赋存直径较小的溶洞；②岩溶含水层中溶洞对导水裂隙有着重要的影响，无溶洞条件下，导水裂隙发育高度为43.1 m，裂采比为14.4，含溶洞条件下，导水裂隙发育异常，其高度可达173.1 m，裂采比高达57.7，与玉龙山段灰岩强含水层沟通；③溶洞下顶板导水裂隙由采动上行裂隙和溶洞失稳下行裂隙两部分构成，采动影响下，溶洞在集中应力和采动附加应力的共同作用下发生失稳，并容易形成下行裂隙，与采动上行裂隙沟通，最终形成黔西南矿区特殊的超高导水裂隙。
Abstract:Southwest Guizhou mining area is a typical karst development mining area in China. Under the influence of mining, the height of karst roof water conducting cracks is abnormally developed. During the rainy season, atmospheric precipitation is extremely easy to enter the underground working face through ultra-high water conducting channels, causing water inrush disasters at the working face, seriously affecting the normal production of the mine. Therefore, based on the analysis of the occurrence characteristics of karst caves in Xintian Coal Mine, the development rules of karst roof water conducting fracture zones during mining were studied by means of on-site measurement, indoor simulation, and theoretical analysis. The development mechanism of ultra-high water conducting fractures was revealed. The results show that: ①The roof karst caves in the study area have obvious zonation phenomenon from top to bottom in the layers such as the surface, the Yulongshan section, and the Changxing Formation. The surface water holes develop along the gullies, and the karst caves in the upper and middle parts of the strong aquifer in the Yulongshan section develop, presenting different forms of beads. The Changxing Formation only locally hosts karst caves with smaller diameters; ② Karst caves in karst aquifers have an important impact on water conducting fractures. Without karst caves, the development height of water conducting fractures is 43.1 m, and the fracture mining ratio is 14.4. Under karst caves, the development of water conducting fractures is abnormal, with a height of 173.1 m, and a fracture mining ratio of 57.7, which communicates with the strong limestone aquifer in the Yulong Mountain section; ③ Karst roof water conducting fissures consist of two parts: mining upward fissures and karst cave instability downward fissures. Under the influence of mining, karst caves become unstable under the combined action of concentrated stress and mining additional stress, and are prone to form downward fissures, which communicate with mining upward fissures, ultimately forming a special ultra-high water conducting fissure in southwestern Guizhou mining area.
-
0. 引 言
近些年，我国聚焦“双碳”目标和绿色发展之路的不断深入推进，选煤是煤炭清洁利用的必经之路，煤泥脱水作为选煤厂生产的核心环节之一，直接影响水资源的和尾煤资源的回收利用，对选煤厂的产品水分、生产效率和环境保护具有重要意义，随着采煤机械化、智能化的发展，细粒煤泥不断增加，原煤煤质变差，煤炭废水的固液分离已经成为行业焦点问题[1-3]。
国内外学者针对煤泥压滤效果、滤饼结构等开展了较多研究，比如任晓汾等[4]研究发现，在加压脱水中添加OP系列的表面活性剂对煤泥有一定的助滤作用，药剂浓度为400 g/t时，滤饼水分降至最低点；夏畅斌等[5]发现在细粒煤浆加入阴离子表面活性剂SDS后，滤饼水分降低至10%，加入阳离子表面活性剂DAB滤饼水分降低至3%～4%；李玉霞等[6]选用不同的表面活性剂对山西选煤厂细粒煤泥水进行了一系列脱水试验研究，研究结果表明，非离子表面活性剂硅油在用量仅为20 g/t脱水效果明显，可降低滤饼水分4.1%。闫奋飞等[7]对比了7种不同表面活性剂的脱水试验，还研究了煤油作为溶剂和表面活性剂混合的新过滤介质对煤泥脱水的影响。结果表明：非离子表面活性剂在试验过程中助滤效率最好，Span-80用量仅为60 g/t时，滤饼水分较空白试验降低3.1%。邓铃等[8-9]通过试验分析了煤泥压滤脱水阶段压力对污泥滤饼水分的影响，结果表明过高的压力会增加固液分离的消耗，加快滤板和滤布的磨损，煤泥水压滤的压力一般控制在0.6~0.8 MPa之间较为合适；石常省[10]等人采用扫描电镜法分析了细粒煤压滤滤饼的微观结构，分析不同粒度细粒煤压滤滤饼结构。赵扬等[11]采用分层切片和扫描电镜法测试了滤饼孔隙率。IRITANI等[12-14]研究了压滤过程中过滤压力和滤饼结构的关系，随着过滤压力的增加，滤饼的平均质量比阻增大，滤饼结构的渗透性能降低。整体而言，针对加压过滤过程的加压脱水效果，尤其是滤饼结构特性仍需进一步研究[15-16]。
利用自主设计的新型加压脱水效果智能监测试验装置，对比了加入2种表面活性剂（OP-10、TX-10）后煤泥水的压滤脱水效果，测量分析了过滤速度、滤饼水分、毛细吸水时间、滤饼平均质量比阻、润湿热和压缩脱水指数等，最后采用CT扫描技术结合迂曲度、渗透率等一系列指标，来微观分析滤饼结构，为煤泥高效加压脱水提供理论和技术指导。
1. 材料与方法
1.1 试验材料
1.1.1 煤泥样品
试验所用煤泥采集自山西马兰矿选煤厂浓缩机底流，经烘干、密封保存后备用。此试验根据国家标准对干煤泥进行工业分析（GB/T 212—2008）、小筛分分析（GB/T 477—2008），并对其进行了元素分析，结果见表1和表2，由表可知，煤泥灰分为63.83%，矿物成分质量分数较高，煤泥挥发分为18.25%，为中高变质程度烟煤，且亲水性元素S、O、N质量分数较高，0.045～0.074 mm产率为26.78%，而煤泥中−0.074 mm的细粒级质量分数为58.79%，说明马兰矿选煤厂的煤泥中细煤泥质量分数高，总体而言，该煤泥体现出典型的高灰细煤泥特性。
表 1 煤泥工业和元素分析Table 1. Proximate and ultimate analysis of coal slime% 工业分析 元素分析 Mad Aad Vd FCd St,d Nd Od Hd Cd 2.68 63.83 18.25 15.24 1.412 0.6 10.02 1.802 22.83 表 2 煤泥粒度组成Table 2. Coal slurry size composition粒径/mn 产率/% 灰分/% 正累计/% 产率 灰分 0.25～0.5 4.77 25.56 4.77 25.56 0.125～0.25 16.95 48.36 21.72 43.35 0.074～0.125 19.49 58.28 41.21 50.41 0.045～0.074 26.78 68.61 67.99 57.58 −0.045 32.01 77.12 100.00 63.83 合计 100.00 63.83 1.1.2 试验药剂
试验所用药剂为OP-10（辛基酚聚氧乙烯醚，C34H62O11，非离子表面活性剂）和TX-10（壬基酚聚氧乙烯醚，C32H58O10，非离子表面活性剂）表面活性剂，无锡市亚泰联合化工有限公司生产，均为分析纯AR级别。
试验用水为马兰矿选煤厂循环水，经过滤除杂后备用。
1.2 加压脱水试验装置
自制试验室智能加压过滤脱水装置如图1所示，装置由压力传动装置、多种传感器及终端数据采集软件等构成，空气压缩机压力范围为0～1 MPa，该装置可实时记录脱水过程中的滤饼厚度、滤液量等指标。
![]()
图 1 智能加压脱水机1—启动按钮；2—压力表；3—压力传感器及其显示屏；4—位移传感器及其显示屏； 5—重量传感器及其显示屏Figure 1. Intelligent pressurised dewatering machines1.3 试验方法
1.3.1 加压脱水试验
1）采用烘干后的煤泥和选煤厂循环水配置100 mL煤泥水（浓度400 g/L），利用电动搅拌器搅拌30 min，加入定量OP-10或TX-10表面活性剂，继续搅拌30 min；
2）打开加压脱水机的空气压缩泵，调节压力至0.6 MPa，将搅拌后的煤泥水倒入样品池中，将样品池安装固定至脱水机内；
3）打开数据采集软件和智能加压脱水机启动开关，开始加压过滤，实时记录滤液体积（Vi）、滤饼厚度。
1.3.2 毛细吸水时间CST
CST测试[17-18]采用台式毛细吸水时间测试仪（厂家：天津恒祥科技有限公司，型号：TYPE 304M）进行，取2.5 mL待分析煤泥水，置于装置的特定漏斗内，记录煤泥内水分从内圈探头扩散到外圈探头所用时间，即为该煤泥的毛细吸水时间。测定出煤泥在调理前后的毛细吸水时间，计算毛细吸水时间降低率。
$$ R = \frac{{{T_{{\rm{CS}},0}} - {T_{{\rm{CS}},1}}}}{{{T_{{\rm{CS}},0}}}} \times 100\text{%} $$ (1) 式中，R为煤泥毛细吸水时间降低率%；TCS, 0为煤泥调理前毛细吸水时间，s；TCS, 1为煤泥调理后毛细吸水时间，s。
1.3.3 滤饼平均质量比阻
根据Ruth过滤理论，在恒定的压力作用下，公式如下：
$$ \int_{0}^{t} {\rm{d}} t=\dfrac{\mu c \alpha}{A^{2} \Delta p}=\int_{0}^{V} V {\rm{d}} V+\dfrac{\mu R_{{\rm{m}}}}{A \Delta p} \int_{0}^{\nu} {\rm{d}} V $$ (2) 积分并整理后可得到下列的线性化抛物线方程：
$$ \frac{t}{V}=\dfrac{\mu c \alpha}{2 A^{2} \Delta p} V+\dfrac{\mu R_{{\rm{m}}}}{A \Delta p} $$ (3) 且有
$$ c=\dfrac{s \rho}{1-s m} $$ (4) 式中：V为滤液量，m3；t为过滤时间，s；μ为滤液黏度，mPa·s；A为过滤面积，m2；c为单位体积滤液的干滤饼质量，kg/m3；ΔP为过滤压力，Pa；Rm为过滤介质的阻力，m−1；m为滤饼质量湿干比；s为物料固体颗粒的质量分数，%；ρ为滤液密度，kg/m3；α为平均质量比阻，m/kg。
上式t/V可以看作非独立变量，V为独立变量。需要先做系列压滤试验测得对应公式的对应参数，对应压滤脱水时间t的滤液体积V，进而在直角坐标系中绘制t/V-V线性拟合图，计算斜率和截距，从而得到平均质量比阻α。过滤速度公式为
$$ \nu=\dfrac{V_{{\rm{i}}}}{{t_{\rm{i}}} {\text{π}} r^{2}} $$ (5) 式中：v为过滤速度，mL/(s·m2)；Vi为滤液体积保持不变，mL；ti为滤液体积保持不变的时间，s；r为样品池半径，cm。
1.3.4 红外光谱分析
采用Thermo Scientific Nicolet iS20 FTIR仪器对煤泥样品进行表面官能团分析。测试流程：通过研磨将KBr和样品以质量比10∶1的比例充分混合，置于模具中，在油压机上压成透明薄片，将样品放入红外光谱仪中测试。波数测量范围为4000～400 cm−1。
1.3.5 润湿热分析
采用法国Setaram公司的C80微热量仪测定2种表面活性剂OP-10、TX-10与煤样两相物质发生反应的润湿热，设定30 ℃恒温状态[19-20]，量热分辨率为0.10 µW。在测试过程中，将200 mg的煤样和2 mL的表面活性剂分别置于膜混合反应池的底部和上部，用锡纸将两者分隔开，反应数小时后等待温度平衡后，刺破锡纸，使二者混合反应即可测量，对记录的热流曲线进行积分，得到OP-10和TX-10的润湿热值。
1.3.6 CT分析方法
样品的CT分析由天津三英精密仪器有限公司nanoVoxel-4000高分辨X射线三维显微镜完成主要由X射线源、平板探测器、精密样品台、图像采集系统、三维图像重建和处理系统等组成。
传统的微观结构表征方法，如光学显微镜可以对滤饼表面进行二维可视化表征，但是无法对滤饼进行三维结构特征，如滤饼迂曲度和连通性。
X射线计算机断层扫描技术（CT扫描技术）具有高分辨率、精度高、扫描范围广的特点[21-23]，利于三维多平面重建。采用独特的X光光学显微成像技术，利用不同角度的X射线透视图像，结合计算机三维数字重构技术，提供样品内部复杂结构的高分辨率三维数字图像，对样品内部的微观结构进行数字化三维表征。
CT扫描技术的基本原理是X射线的衰减作为生成X射线图像的信号。这意味着对比度信号与电子密度密切相关。由于X射线衰减成像的信噪比很大程度上取决于3D图像的每个体素检测到的X射线光子数量，即X射线曝光必须随着空间分辨率的增加而增加，因此，它的敏感性最终受到其潜在的辐射损伤的限制。
2. 结果与讨论
2.1 煤泥样品化学性质
由图2可知，煤泥中含有大量的高岭石、石英、方解石等矿物质，由于高岭石是亲水性矿物，在煤泥水中发生吸附，颗粒表面有较高的负电荷，进而在煤泥脱水过程中，颗粒的保水性及空间位阻效应，导致煤泥水脱水困难。
为探究煤泥的官能团，对煤样进行了红外分析，如图3所示。在3694 cm−1与3620 cm−1处所观察到的吸收峰为游离羟基（—OH）的伸缩振动；波数为2920 cm−1的谱峰为亚甲基（—CH2）反对称伸缩振动；在1610 cm−1形成的谱峰为羰基（—C=O）的伸缩振动和羰基（—C=O）与羟基（—OH）形成的氢键的共振形成的吸收峰；在1430 cm−1出现谱峰，说明了此处出现了芳香 C=C、—CH2、—CH3的伸缩振动；1040～910 cm−1区域内和797 cm−1、696 cm−1的峰为煤样中脉石矿物所形成的吸收峰，主要是 Si—O—Si、Si—O—Al 等的伸缩振动；754 cm−1的峰为Al3+在铝氧八面体的振动；在537 cm−1处的谱峰为芳香双硫醚（—S—S—）的吸收峰。
从上面的分析可知，本试验样品表面上具有丰富的含氧官能团（—OH、C=O等）以及矿物杂质基团（Si—O—Si、Si—O—Al等），这些基团的存在说明所用样品表面亲水性较强，从而致使过滤脱水困难，滤饼水分偏高。因此，选择合适的脱水药剂来降低滤饼水分是关键因素。
2.2 煤泥加压脱水效果
2.2.1 压滤速度与滤饼水分
图4为OP-10和TX-10作用下的煤泥脱水速度曲线，由图可知，未加药剂时，煤泥脱水速度为111.36 mL/(s·m2)，进入压缩阶段时间为150 s（拐点位置）。OP-10和TX-10作用下，随着药剂量增加，煤泥脱水速度均呈现先增加后降低的趋势，OP-10在药剂量为200 g/t时，过滤速度达到最大值209.86 mL/(s·m2)，进入压缩阶段时间为85 s，在药剂量为100 g/t和500 g/t时，脱水速度分布为154.76 mL/(s·m2)和190.37 mL/(s·m2)；TX-10药剂量为400 g/t时，过滤速度达到最大值196.63 mL/(s·m2)，进入压缩阶段时间为100 s，在药剂量为100 g/t和500 g/t时，脱水速度分布为134.85 mL/(s·m2)和176.89 mL/(s·m2)。
![]()
图 4 OP-10和TX-10作用下的煤泥脱水速度曲线Figure 4. Dewatering speed curve of coal slurry under the action of OP-10 and TX-10由图5可知，OP-10和TX-10药剂量对滤饼水分的影响规律与对脱水速度的影响规律较为一致，未加药剂情况下，原煤泥滤饼水分为20.5%，加入OP-10情况下，药剂量200 g/t时，滤饼水分最低（16.7%），比未加药剂情况下降低3.8%，TX-10在药剂浓度为400 g/t时，滤饼水分达到最低值18.2%，相比未加药剂情况下降2.2%。综上所述，2种药剂均起到较为显著的助滤作用，且OP-10比TX-10更有助于提高煤泥水的压滤速度、并降低滤饼水分。
![]()
图 5 OP-10和TX-10作用下的滤饼水分Figure 5. Filter cake moisture in the presence of OP-10 and TX-102.2.2 滤饼比阻与毛细吸水时间
平均质量比阻和毛细吸水时间是定量评价物料脱水性能的重要指标，由图6可知，2种表面活性剂中，脱水容易的是OP-10。OP-10浓度在0～200 g/t区间，平均质量比阻急速下降，浓度为200 g/t时滤饼平均质量比阻最小，其比阻为5.21×106 m/kg，相比空白试验降低了2.74×107 m/kg。脱水较难得是TX-10，在浓度为400 g/t时比阻达到最低点9.31×106 m/kg，相比浓度200 g/t的OP-10比阻高了4.09×106 m/kg。滤饼平均质量比阻越小，水分就越易从滤饼中脱除，加压脱水速度越快，滤饼水分也越低。
![]()
图 6 滤饼平均质量比阻与毛细吸水时间Figure 6. Average filter cake mass specific resistance vs. capillary wicking time CST加入OP-10、TX-10 这2种表面活性剂调理煤泥30 min后，CST降低率均升高，毛细吸水时间反映了煤泥压滤脱水性能[24]。2种表面活性剂中，脱水性能较好的是OP-10，OP-10随着药剂浓度增加CST降低率先增加后降低，浓度在200 g/t CST降低率达到最大值14.25%。脱水性能较差的是TX-10，当TX-10用量为100 g/t时，随着药剂浓度的不断增加，CST降低率随之不断增加，到400 g/t时CST降低率达到最大值13.36%，再减小。
2.3 滤饼压缩性变化规律
2.3.1 滤饼厚度
为研究不同表面活性剂对滤饼厚度的影响，加入不同浓度的TX-10，300、400和500 g/t时达到滤饼厚度稳态时用时较短，同时也研究了表面活性剂不同浓度的影响（图7）。由图7可得，加入TX-10药剂压滤而成的滤饼厚度范围在6.497～6.717 mm，形成的滤饼厚度相比加入OP-10药剂的滤饼厚度较大。随着TX-10的不断加入，滤饼厚度先减小后增大，在浓度为400 g/t时，达到滤饼厚度稳态时用时最短为276 s，滤饼厚度达到同组最小为6.497 mm。李满等[25]以−0.5 mm粒级原煤为试验煤样，分析了随着絮凝剂用量的增加，煤泥水过滤时间显著降低，滤饼水分和滤饼厚度存在正相关的关系。
加入不同浓度的OP-10，药剂量为200 g/t时滤饼厚度保持不变的时间最快249.5 s，在煤泥压滤过程中，随着粗重颗粒的快速沉降，细颗粒及水分从颗粒之间的孔隙穿过，脱水速度越快，滤饼水分越小，滤饼厚度越小。OP-10在用量为200 g/t时滤饼厚度最小，仅为6.409 mm。
2.3.2 压缩脱水指数
根据滤饼的压力−形变曲线可以定义压缩脱水指数(ICD)：
$$ {I_{{\rm{CD}}}} = \frac{{\sigma (R) - \sigma (O)}}{{\varepsilon (R) - \varepsilon (O)}} $$ (6) 式中，σ(R)、ε(R)为R点对应的压力值和形变值。
图8为2种表面活性剂的压力位移曲线，O点是压滤脱水阶段结束，也是滤饼压缩阶段开始，R点对应的是滤饼压缩结束时，在滤饼压缩这一阶段计算出压缩脱水指数，来反应滤饼的可压缩性，压缩脱水指数的几何意义是压力位移曲线上直线段的斜率。在压缩脱水时具有弹性的滤饼将通过释放水来减少滤液体积，以便消散多余的应力。由图可以看出过滤速度越快的煤泥，所需的压力越小。加入不同浓度的OP-10，浓度为200 g/t时形成滤饼所需的压差最小，仅需5110.83 Pa。加入不同浓度的TX-10，浓度为400 g/t时形成滤饼所需的压差为5505.73 Pa。
煤泥水在压缩脱水形成滤饼时是一种弹性材料，是一种可压缩性材料，ICD从宏观角度来说是体积随压力变化的尺度，从微观角度来说是固体颗粒的微观迁移和团聚能力的大小。图9为表面活性剂对压缩脱水指数的影响，根据图9可得，2种表面活性剂的ICD变化趋势与其滤饼水分一致，2种表面活性剂在可压缩性方面，OP-10的效果最好，在药剂浓度在药剂浓度为200 g/t时，ICD由未加药剂空白试验的238.35降低至147.18，降低了91.17。相比可压缩性较差的是TX-10，在药剂浓度为400 g/t时ICD才达到最低点155.35。ICD越小，在相同压力下滤饼厚度降低幅度越大，滤饼水分越低，针对多孔介质压缩性强，更易过滤脱水。
2.4 药剂与煤泥作用的润湿热分析
OP-10和TX-10作用下的润湿热流线和润湿热值如图所示。润湿热反映的是2种物质之间的接触能力，润湿热值越大，说明2种物质之间的作用程度越大，结合的也就越紧密。将200 mg的煤样和2 mL的表面活性剂（OP-10浓度：200 g/t，TX-10浓度：400 g/t）分别置于膜混合反应池的底部和上部，等待温度平衡，刺破锡纸，让两者充分混合反应，用Data acquisition软件记录润湿过程的热流曲线，如图10所示，从图中可以看出TX-10的热流线峰值高于OP-10，表明在润湿过程中TX-10的放热速率大于OP-10。图10b中可以看出，其润湿热值都为负值，表明这是自发的放热过程。
![]()
图 10 OP-10和TX-10的润湿热差异分析Figure 10. Analysis of the difference in the heat of wetting between OP-10 and TX-10OP-10的润湿热为1.196 J/g，TX-10的润湿热为4.11 J/g， TX-10比OP-10的润湿热大70.9%，表明煤泥和TX-10之间的亲和力强，煤泥表面容易被TX-10润湿，为亲水性表面，相对OP-10难脱水，因此OP-10表面活性剂的压滤脱水效果好[26]。
2.5 滤饼结构特性CT分析
2.5.1 CT图像处理及三维重建
为了进一步研究过滤滤饼的特点，选取了原煤和加入浓度为200 g/t的OP-10形成的滤饼，通过nanoVoxel-4000高分辨X射线3维显微镜对所要分析的煤样进行扫描，按照CT试验步骤后，得到相关数据和图像。将数据导入Avizo软件中，进行图像处理，图11为CT的原始二维切片，图12为中值滤波处理效果图，为了更好地区分孔隙和颗粒骨架，因此对图像进行中值滤波的操作，为了保护图像的细节，并且降噪，从而提升了图像的质量，以便于后续的图像处理。
2.5.2 CT图像阈值分割
为了对煤样孔隙结构的分析，通过裁减分割选取合适的分析范围来代表煤样，随后对中值滤波后的样品进行交互式阈值分割[27-28]，分割成孔隙可以计算出总孔隙率。图13为对煤样切片使用阈值分割法得到的图像，蓝色代表孔隙，黑色代表颗粒骨架。
2.5.3 样品三维重建及滤饼孔隙尺寸分析
通过二维的煤样切片进行煤样内部结构研究有一定的难度，所以需要对煤样切片进行三维重建，图14为样品三建重建结果，图15为三维滤饼孔隙空间分割结果，图15a，15b是原样分割前后的孔隙空间，图15c，15d为OP-10样品分割前后的孔隙空间。在三维重建后的图像中添加Label Analysis命令对滤饼的每个孔隙进行单独分析，通过等效直径可以得到三维滤饼孔隙空间的孔隙数量、孔隙体积及表面积分布[28]。通过计算得到，原样和加了OP-10药剂滤饼孔的数量分别为1448和1352，滤饼的孔体积27964.98 μm3和45025.91 μm3，孔面积为52994.88 μm2和68423.45 μm2，比表面积分别为1.90 μm−1和1.52 μm−1，原样的比表面积大，根据固体表面吸附理论，比表面积大的多孔物，吸附能力强，所以原样滤饼的水分明显高于加了OP-10表面活性剂的样品。
2.5.4 滤饼孔隙连通性及迂曲度分析
由表3可知，压滤脱水成型的滤饼都存在一定的孤立孔隙，加入OP-10所形成的滤饼内部孤立孔隙所占总孔隙体积比例与原样相比很小（OP-10为2.53%，原样为8.14%），它对孔隙的整体连通性影响非常小，而未加药剂的原样中孤立孔隙占据较大的比例8.14%，连通孔体积相比OP-10也少了18198 μm3，导致原样孔隙的连通性较差，其有效孔隙率也明显下降。由上述分析可知，颗粒的粒度对滤饼的连通性也有一定的影响，不加药剂的原样相比加OP-10的过滤速度慢，原样颗粒体积（186885 μm3）大于OP-10作用后的体积（169823 μm3），以致于未加药剂的滤饼孔道较窄，压滤脱水更容易形成孤立的孔隙，有效孔隙占比较低（11.96%），滤液流动的道路减少，导致压滤脱水困难[29]。
表 3 孔隙连通性分析结果Table 3. Analysis results of pore connectivity样品 连通孔
体积/μm3孤立孔
体积/μm3颗粒
体积/μm3总孔隙
率/%有效孔隙
率/%孤立孔
比例/%原样 25691 2274 186885 13.02 11.96 8.14 OP-10 43889 1137 169824 20.96 20.43 2.53 孔隙的迂曲度是指液体流动的实际长度与压力梯度方向上线距离的比值，流体在多孔介质中的流动不是沿直线前进，迂曲度就是反映这种迂回曲折的程度。图16为孔隙中轴线提取结果，由Avizo计算得原样和OP-10的迂曲度分别为2.37和1.9。加入OP-10样品的滤饼孔隙率大，迂曲度减小，滤饼的渗透率逐渐增大，平均质量比阻就越小，过滤速度就越快。
3. 结 论
1）OP-10和TX-10两种表面活性剂对马兰浓缩机底流提取的样品均有良好的助滤效果，与空白试验结果相比，当药剂用量为200 g/t时，OP-10使滤饼水分降低至16.7％，过滤速度提高200%，CST降低率最大为14.25%。
2）在压缩脱水阶段进一步定义了压缩脱水指数CDI，CDI越低，滤饼可压缩性越高，相同压力下滤饼厚度降低幅度越大，水分越低，最佳浓度下的OP-10和TX-10作用下的压缩脱水指数分别为147.18和155.35。
3）通过CT扫描技术，利用Avizo软件可以得到三维立体的孔隙结构，微观分析原样和加入OP-10的脱水效果，加入OP-10浓度200 g/t时，可以通过相应的运算模块得到有效孔隙率为20.43%，连通性好，原样的有效孔隙率11.96%，连通性差，孤立孔比例大，过滤脱水困难。
4）由计算得原样和OP-10的迂曲度分别为2.37和1.9。加入OP-10样品的滤饼孔隙率大，迂曲度减小，滤饼的渗透率逐渐增大，平均质量比阻就越小，过滤速度就越快。
-
![]()
图 1 新田矿井落水洞与岩溶通道分布
Figure 1. Distribution of sinkholes and karst channels in Xintian Mine
![]()
图 3 工作面涌水量与大气降水关系
Figure 3. Relationship between water inflow and atmospheric precipitation at working face
![]()
图 13 工作面开挖100 m覆岩裂隙发育
Figure 13. Development of overlying rock fractures after excavation of 100 m on the working face
![]()
图 14 工作面开挖180 m覆岩裂隙发育
Figure 14. Development of overlying rock fractures after excavation of 180 m on the working face
![]()
图 15 工作面开挖290 m覆岩裂隙发育
Figure 15. Development of overlying rock fractures after excavation of 290 m on the working face
![]()
图 16 工作面开挖350 m覆岩裂隙发育
Figure 16. Development of overlying rock fractures after excavation of 350 m on the working face
![]()
图 17 工作面开挖480 m覆岩裂隙发育
Figure 17. Development of overlying rock fractures after excavation of 480 m on the working face
表 1 矿井岩溶发育特征
Table 1 Characteristics of karst development in mines
钻孔编号 孔口标高/m 岩溶发育层位与深度 岩溶发育程度 ZK1 +1 293 玉龙山段（140 m以上段） 局部有溶孔及溶蚀裂隙发育，岩溶发育较差 ZK2 +1 238 玉龙山段（132 m以上段） 溶洞、裂隙多，且多有砾石、泥砂充填, 岩溶发育强烈 长兴组(269 m) 局部发育有小溶孔,层面有溶蚀痕迹,岩溶发育差 ZK3 +1 298 玉龙山段（153 m） 局部发育有小溶孔及溶蚀裂隙，岩溶发育较差 检2孔 +1 245 玉龙山段（72～211 m） 溶孔及裂隙较发育，岩溶发育中等 长兴组(211～290 m) 局部有溶孔及裂隙发育，岩溶发育较差 检3孔 +1 290 玉龙山段（172～186 m） 为溶洞，下部有黏土、卵石充填，岩溶发育强烈 检4孔 +1 244.54 玉龙山段（0～200 m） 0～127 m为落水洞充填物，127～200 m溶孔、裂隙发育，岩溶发育强烈 301 +1 238 长兴组（209 m以上段） 蜂窝状溶孔,岩溶发育较差 
下载: 导出CSV
表 2 相似模拟试验配比
Table 2 Comparison of similar simulation tests
序号 岩性 密度/（kg·m−3） 抗拉强度/MPa 岩层厚度/m 模型厚度/cm 材料质量/kg 沙子 石膏 碳酸钙 1 钙质泥岩 2 659 2.54 18 18 153.6 9.2 1.3 2 燧石灰岩 2 612 2.3 26 26 210.6 21.1 21.1 3 泥质粉砂岩 2 520 1.7 2 2 17.3 2.4 1.0 4 含燧石灰岩 2 612 2.54 6 6 50.6 3.0 7.1 5 泥质粉砂岩 2 520 1.7 6 6 52.0 7.3 3.1 6 粉砂质泥岩 2 595 2.54 5 5 43.3 2.6 6.1 7 泥质粉砂岩 2 520 3.9 5 5 43.3 2.6 6.1 8 细砂岩 2 550 2.5 5 5 42.5 2.6 6.0 9 泥质粉砂岩 2 520 1.7 2 2 17.3 2.4 1.0 10 泥质灰岩 2 600 2.70 2 2 16.3 1.0 2.3 11 泥质粉砂岩 2 520 2 3 3 26.0 2.6 2.6 12 4号煤 1 400 0.28 3 3 14.0 2.0 0.8 13 泥质粉砂岩 2 520 1.7 25 25 216.7 30.3 13.0
下载: 导出CSV
表 3 钻孔施工参数
Table 3 Drilling parameters
孔号 方位角/(°) 倾角/(°) 孔深/m 垂高/m 1号 35 64 96.40 68 2号 45 63 96.90 68 3号 215 66 94.50 68
下载: 导出CSV
表 4 相似模拟试验配比
Table 4 Comparison of similar simulation tests
序号 岩性 密度/（kg·m−3） 抗拉强度/MPa 岩层厚度/m 模型厚度/cm 材料质量/kg 沙子 石膏 碳酸钙 1 石灰岩 2 700 4.40 180.0 120.0 2451.5 183.9 429.0 2 钙质泥岩 2 659 2.54 18.0 12.0 258.0 11.1 25.8 3 燧石灰岩 2 612 2.30 26.0 17.3 353.1 25.2 25.2 4 泥质粉砂岩 2 520 1.70 2.0 1.3 28.4 2.8 1.2 5 含燧石灰岩 2 612 2.54 6.0 4.0 85.0 3.6 8.5 6 泥质粉砂岩 2 520 1.70 6.0 4.0 87.4 8.7 3.7 7 粉砂质泥岩 2 595 2.54 5.0 3.3 71.9 3.1 7.2 8 泥质粉砂岩 2 520 3.90 5.0 3.3 70.6 3.5 8.2 9 细砂岩 2 550 2.50 5.0 3.3 70.7 3.0 7.1 10 泥质粉砂岩 2 520 1.70 2.0 1.3 28.4 2.8 1.2 11 泥质灰岩 2 600 2.70 2.0 1.3 26.8 1.1 2.7 12 泥质粉砂岩 2 520 2.0 3.0 2.0 43.7 3.1 3.1 13 4号煤 1 400 0.28 3.0 2.0 23.9 2.1 0.9 14 泥质粉砂岩 2 520 1.70 25.0 16.7 364.7 36.5 15.6
下载: 导出CSV
-
[1] 周 泽. 岩溶矿区采动裂隙发育及其地表塌陷规律研究[D]. 湘潭: 湖南科技大学, 2017. ZHOU Ze. Study on the development of mining induced fractures and the law of surface subsidence in karst mining areas [D]. Xiantan: Hunan University of Science and Technology, 2017.
[2] 余龙哲. 溶洞地貌下浅埋煤层开采覆岩破断规律研究[D]. 贵阳: 贵州大学, 2021. YU Longzhe. Research on overburden breaking law of shallow coal seam mining under karst cave landform [D]. Guiyang: Guizhou University, 2021.
[3] 朱川曲,崔栋歌,周 泽,等. 岩溶矿区采动裂隙发育及溶洞破坏特征相似模拟[J]. 地下空间与工程学报,2019,15(1):93−100,124. ZHU Chuanqu,CUI Dongge,ZHOU Ze,et al. Similar simulation of mining-induced fissure development law and cave destruction character when mining in karst area[J]. Journal of Underground Space and Engineering,2019,15(1):93−100,124.
[4] 王文强. 新田矿岩溶顶板覆岩运动规律及涌水控制技术研究[D]. 焦作: 河南理工大学, 2021. WANG Wenqiang . Research on the movement law of karst roof overburden and water burst control technology in Xintian Mine [D]. Jiaozuo: Henan Polytechnic University, 2021
[5] 范立民. 保水采煤面临的科学问题[J]. 煤炭学报,2019,44(3):667−674. FAN Limin. Scientific problems faced with water conservation and coal mining[J]. Journal of China Coal Society,2019,44(3):667−674.
[6] 侯恩科,袁 峰,王双明,等. 导水裂隙带发育特征地震识别方法[J]. 煤炭学报,2023,48(1):414−429. HOU Enke,YUAN Feng,WANG Shuangming,et al. Seismic identification method for development characteristics of water conducting fracture zones[J]. Journal of China Coal Society,2023,48(1):414−429.
[7] 鞠金峰,马 祥,赵富强,等. 东胜煤田导水裂隙发育及其分区特征研究[J]. 煤炭科学技术,2022,50(2):202−212. doi: 10.13199/j.cnki.cst.qntk21-1396 JU Jinfeng,MA Xiang,ZHAO Fuqiang,et al. Study on the development and zoning characteristics of water conducting fractures in Dongsheng coalfield[J]. Coal Science and Technology,2022,50(2):202−212. doi: 10.13199/j.cnki.cst.qntk21-1396
[8] 陈 萍. 岩溶山区矿坑顶板冒落导水裂隙带高度估算有关问题讨论: 以贵州省黔西县青龙煤矿为例[J]. 贵州地质,2015,32(1):53−58. CHEN Ping. Discussion on the estimation of the height of water conducting fractures in the roof of mining pits in karst mountain areas: taking qinglong coal mine in qianxi county, guizhou province as an example[J]. Guizhou Geology,2015,32(1):53−58.
[9] 胡炳南, 张华兴, 申宝宏. 建筑物、水体、铁路及主要井巷煤柱留设与压煤开采指南[M]. 北京: 煤炭工业出版社, 2017. HU Bingnan, ZHANG Huaxing, SHEN Baohong. Guidelines for coal pillar retention and coal mining in buildings, water bodies, railways and main shafts [M]. Beijing: Coal Industry Press, 2017.
[10] 刘天泉. 矿山岩体采动影响与控制工程学及其应用[J]. 煤炭学报, 1995, 20(1): 1−5. LIU Tianquan. Mining influence and control engineering of mine rock mass and its application [J] Journal of China Coal Society, 1995, 20 (1): 1−5.
[11] 许家林,朱卫兵,王晓振. 基于关键层位置的导水裂隙带高度预计方法[J]. 煤炭学报,2012,37(5):762−769. XU Jialin,ZHU Weibing,WANG Xiaozhen. Prediction method for the height of water conducting fracture zone based on the location of key layers[J]. Journal of China Coal Society,2012,37(5):762−769.
[12] 杨达明,郭文兵,赵高博,等. 厚松散层软弱覆岩下综放开采导水裂隙带发育高度[J]. 煤炭学报,2019,44(11):3308−3316. YANG Daming,GUO Wenbing,ZHAO Gaobo,et al. Development height of water-conducting fracture zone in fully mechanized top-coal caving mining under soft overburden of thick loose layer[J]. Journal of China Coal Society,2019,44(11):3308−3316.
[13] 李超峰,刘英锋,李抗抗. 导水裂隙带高度井下仰孔探测装置改进及应用[J]. 煤炭科学技术,2018,46(5):166−172. LI Chaofeng,LIU Yingfeng,LI Kangkang. Improvement and application of underground up-hole detection device at the height of water-conducting fracture zone[J]. Coal Science and Technology,2018,46(5):166−172.
[14] 张国奇,吕晓磊,李东发,等. 半固结砂岩含水层下综放开采“两带”高度实测研究[J]. 煤炭科学技术,2017,45(S1):183−185,195. ZHANG Guoqi,LYU Xiaolei,LI Dongfa,et al. Measured research on the height of “two zones” in fully mechanized caving mining under semi-consolidated sandstone aquifer[J]. Coal Science and Technology,2017,45(S1):183−185,195.
[15] 张玉军,李凤明. 高强度综放开采采动覆岩破坏高度及裂隙发育演化监测分析[J]. 岩石力学与工程学报,2011,30(S1):2994−3001. ZHANG Yujun,LI Fengming. Monitoring and analysis of overburden failure height and fracture development evolution caused by high intensity fully mechanized caving mining[J]. Journal of Rock Mechanics and Engineering,2011,30(S1):2994−3001.
[16] 张 军,王建鹏. 采动覆岩“三带”高度相似模拟及实证研究[J]. 采矿与安全工程学报,2014,31(2):249−254. ZHANG Jun,WANG Jianpeng. High similarity simulation and empirical study of mining overburden “three zones”[J]. Journal of Mining and Safety Engineering,2014,31(2):249−254.
[17] 刘跃俊,张宁波. 水体下急倾斜煤层开采覆岩破断规律[J]. 煤炭科学技术,2018,46(S1):72−74,111. LIU Yuejun,ZHANG Ningbo. Overburden breaking law of mining steep seam under water[J]. Coal Science and Technology,2018,46(S1):72−74,111.
[18] WANG Yongqian, WANG Xuan, ZHANG Jiasheng, et al. Roof subsidence and movement law of composite strata mining: insights from physical and numerical modeling[J]. Minerals, 2021. DOI: 10.3990/min12010003"> 10.3990/min12010003
[19] 李蕊瑞,陈陆望,欧庆华,等. 考虑覆岩原生裂隙的导水裂隙带模拟[J]. 煤田地质与勘探,2020,48(6):179−185,194. LI Ruirui,CHEN Luwang,OU Qinghua,et al. Simulation of water-conducting fracture zone considering primary fractures of overburden[J]. Coal Geology and Exploration,2020,48(6):179−185,194.
[20] 周宏伟,张 涛,薛东杰,等. 长壁工作面覆岩采动裂隙网络演化特征[J]. 煤炭学报,2011,36(12):1957−1962. ZHOU Hongwei,ZHANG Tao,XUE Dongjie,et al. Evolution of mining-induced crack network in overburden strata of longwall face[J]. Journal of China Coal Society,2011,36(12):1957−1962.
[21] 高明忠,金文城,郑长江,等. 采动裂隙网络实时演化及连通性特征[J]. 煤炭学报,2012,37(9):1535−1540. GAO Mingzhong,JIN Wencheng,ZHENG Changjiang,et al. Real-time evolution and connectivity of mined crack network[J]. Journal of China Coal Society,2012,37(9):1535−1540.
[22] 李振华,许延春,李龙飞,等. 基于BP神经网络的导水裂隙带高度预测[J]. 采矿与安全工程学报,2015,32(6):905−910. LI Zhenhua,XU Yanchun,LI Longfei,et al. Prediction of the height of hydraulic fracture zone based on BP neural network[J]. Journal of Mining and Safety Engineering,2015,32(6):905−910.
[23] 娄高中,谭 毅. 基于PSO-BP神经网络的导水裂隙带高度预测[J]. 煤田地质与勘探,2021,49(4):198−204. LOU Gaoxue,TAN Yi. Prediction of the height of water-conducting fracture zone based on PSO-BP neural network[J]. Coal Geology and Exploration,2021,49(4):198−204.
[24] 谭 毅,郭文兵,杨达明,等. 非充分采动下浅埋坚硬顶板“两带”高度分析[J]. 采矿与安全工程学报,2017,34(5):845−811. TAN Yi,GUO Wenbing,YANG Daming,et al. Height analysis of "two zones" of shallow-buried hard roof under insufficient mining[J]. Journal of Mining and Safety Engineering,2017,34(5):845−811.
[25] 王 旭, 尹尚先, 徐 斌, 等. 综采工作条件下覆岩导水裂隙带高度预测模型优化研究[J/OL]. 煤炭科学技术: 1−13. [2023-06-27]. DOI: 10.13199/j.cnki. cst.2022-1530. WANG Xu, YIN Shangxian, XU Bin, et al. Research on optimization of prediction model of overburden water conducting fracture zone height under fully mechanized mining conditions [J/OL]. Coal Science and Technology: 1−13. [2023-06-27]. DOI: 10.13199/j.cnki.cst.2022-1530.
-
期刊类型引用(5)
1. 马晓敏，孙筱淞，董宪姝，樊玉萍，常明. 压滤过程中煤泥滤饼结构演化规律及组合药剂优化机理. 煤炭学报. 2024(03): 1611-1624 . 
百度学术
2. 杨长华，董宪姝，陈茹霞，樊玉萍，马晓敏，冯泽宇. 气泡对浮选精煤压滤压力及脱水性能影响机理. 洁净煤技术. 2024(05): 172-179 . 百度学术
3. 史博文，骆建营，郑利祥，王建，郭中权，韩朝. 响应面法优化CTAB-APAM复合助滤剂强化煤泥脱水效果. 煤炭技术. 2024(12): 264-268 . 百度学术
4. 宋洋，王宏帅，李昂，王鑫，肖作明，杜春生，李志新. 富水粉细砂层盾尾同步注浆滤饼生长速率研究. 铁道科学与工程学报. 2023(08): 2974-2986 . 百度学术
5. 冯泽宇，董宪姝，陈茹霞. 基于分形特征的煤泥滤饼孔渗关系模型研究. 煤炭科学技术. 2023(10): 312-322 . 本站查看
其他类型引用(0)
计量
- 文章访问数: 187
- HTML全文浏览量: 45
- PDF下载量: 98
- 被引次数: 5
