Construction and implementation of emission reduction and treatment technology system in deep mining of high salt mine water
-
摘要:
近年来，随着我国煤矿矿井水排放标准的不断提高，针对高盐矿井水处理及资源化利用方面的技术创新需求进一步提升。减少矿井水尤其是高盐矿井水的排放、提高矿井水的可利用性是提高煤炭生产安全、加快绿色矿山建设、促进煤炭企业可持续发展的重要保障。针对张双楼煤矿深部开采过程中面临的矿井水涌水量大、矿化度高、处理利用成本高等问题，提出了该矿高盐矿井水以“减量、储存、净化”为核心的减排治理技术体系，即：通过底板涌水封堵减量技术对已有残余涌水点进行减量封堵，从源头上减少高盐矿井水的排放；通过矿井水深层回灌储存与保水修复技术将高盐矿井涌水进行深部转移储存，助力实现矿井水“零排放”与深层地下水的保水修复；通过采空区煤岩自净与调蓄技术进行高盐矿井水的煤岩自净预处理，并进一步通过地面深度脱盐处理技术，实现高盐矿井水净化处理以及资源化利用或达标排放。同时，结合该矿的水文地质背景、矿井涌水规律以及矿井水处理利用现状，对以上技术在该矿应用实现的可行性和途径进行了综合分析与评价：通过底板涌水封堵减量技术的实施，实现了188 m3/h的底板四灰涌水减量治理；通过深层回灌储存技术，西翼采区能够实现对奥灰含水层200 m3/h的回灌储存与保水修复；通过采空区煤岩自净预处理与调蓄技术，东翼采区可蓄存约165万m3矿井水并进行特征组分的井下预处理；通过地面深度脱盐处理技术，矿井水处理能力达到700 m3/h，每年可减少盐排量14 271.5 t。最后，提出并形成了张双楼煤矿高盐矿井水减排治理的综合技术思路。张双楼煤矿高盐矿井水减排治理技术体系的提出，可为该矿及其他类似深部开采煤矿高盐矿井水的减量治理以及资源化利用提供理论基础与技术借鉴依据。
Abstract:With the improvement of mine water discharge standards, the demand for technological innovation in mine water treatment and resource utilization has been further enhanced. Reducing the discharge of mine water, especially high salt mine water, and improving the availability of mine water are important guarantees for improving coal production safety, accelerating green mine construction, and promoting sustainable development of coal enterprises. Focused on the problems of high water inflow, high salinity and high cost of treatment and utilization during deep mining of Zhangshuanglou Coalmine, a treatment technology system of mine water emission reduction in high salt mine with “reduction, storage and purification” as the core was put forward: Blocking the water inflow points through the floor water inflow blocking reduction technology to reduce the discharge of high salt mine water from the source; Through the deep reinjection technology of mine water, some mine water will be transferred and stored in depth to realize the “zero discharge” and water conservation of mine water; Through the technology of coal and rock self-purification in goaf, the high-salt mine water is carried out by regulation and storage; Through surface deep desalting treatment technology, high salt mine water purification treatment, resource utilization or standard discharge can be realized. Further, combined with the hydrogeological background of the mine, the law of mine water gushing and the status quo of mine water treatment and utilization, the feasibility of the application of the above technology in the mine was comprehensively analyzed and evaluated: Through the implementation of floor water gushing plugging reduction technology, the reduction treatment of 188 m3/h floor four ash gushing water was realized; Through deep reinjection technology, 200 m3/h recharge storage and water conservation restoration of Ordovician limestone aquifer can be realized in the west wing mining area; About 1.65 million m3 of mine water can be stored in the east wing mining area through self-cleaning pretreatment of coal and rock in goaf and regulation storage technology, and the underground pretreatment of characteristic components can be carried out; Through the surface deep desalination treatment technology, the mine water treatment capacity reaches 700 m3/h, and the annual salt discharge can be reduced by 14 271.5 t. Finally, the comprehensive technical idea of water emission reduction treatment of high salt mine in Zhangshuanglou Coalmine is put forward. The proposal of the water emission reduction treatment technology system of the high-salt mine in Zhangshuanglou Coalmine can provide a theoretical basis and technical reference for the water reduction treatment and resource utilization of high-salt mine in Zhangshuanglou Coalmine and other similar deep mining coal mines.
-
0. 引 言
煤岩体中裂隙的赋存状态、空间结构形态、起裂形式及断裂韧性等属性是影响煤岩体的物理与力学性质的重要因素。对于裂隙多角度形态特征与力学属性的研究一直是工程应用的重要内容。
煤岩体裂隙形态的定量化描述可以表征裂隙演化特征及煤岩受扰动的损伤情况，而可视化技术在工程方面的应用[1]为煤岩体裂隙形态结构的研究开启了一个新方向。许多学者以裂隙形态可视化及属性计算为出发点开展了相关研究。王登科等[2]利用VG Studio软件计算了重构煤岩体裂隙的体积、表面积及裂隙率等特征参数，研究了受载条件下煤岩裂隙演化规律及形态变化特征。刘玉龙等[3]基于X射线层析扫描技术及三维重构技术，实现了煤岩孔裂隙分布的三维可视化精细化表征。刘国方等[4]利用Soble算子和编程将扫描数据进行处理并重构出三维裂隙结构，建立了精确的含原生裂隙的煤岩三维细观数值模型。张平等[5]基于工业CT扫描并重构含裂隙煤岩体，结合VG Stutio MAX软件实现了裂隙宽度、表面积和体积等的定量计算。CAI等[6]采用CT扫描技术研究了三轴压缩条件下裂隙几何特征、连通率、密度及应力变化对裂隙煤岩渗透率的影响规律。
对于存在初始缺陷的煤岩体，内部裂隙参数是影响其力学性能的重要参考依据。因此，学者们从裂隙倾角、长度及分布位置等角度开展了相应研究。赵延林等[7]在三轴压缩条件下，研究了灰岩倾角的不同对岩石破坏过程及峰值的影响规律。秦楠等[8]制作了单裂隙类岩石试件来研究低蠕变作用对试件强度及裂隙扩展规律的影响。王莉等[9]研究了初始裂隙倾角及长度参数对岩石抗拉特性的影响。针对不同张开型裂隙，NIU[10]利用RFPA数值模拟软件模拟发现煤岩体强度随着张开型裂隙长度增大而下降。由此得出，裂隙的力学性质、形态结构及分布特征参数对岩石力学研究具有重要意义，随着研究的不断深入，裂隙的属性数据及相关研究成果对工程研究具有很好的指导作用。
现有研究成果表明，可利用数据库将裂隙结构形态和力学属性进行整合、管理与存储，不断深化裂隙研究。在力学与岩土工程等领域，国内外专家学者已建立了多种数据库。康红普等[11]通过整合多源煤矿地应力数据，建立了中国煤矿井下地应力数据库，分析了我国煤矿地应力分布特征。王宏伟等[12]构筑了弹性力学经典解答数据库，实现了弹性力学问题的可视化描述。李震等[13]建立了露天矿道路网数据库，实现了道路网数据的可视性并提高了露天矿管理效率。LANG等[14]建立了岩土工程勘察数据库系统，实现了勘察工程数据的可视化。上述研究工作表明，数据库的出现使得工程及试验数据得到有效利用，可在实现数据分析及可视化的同时极大地提高工作效率。但是，煤岩体裂隙属性的智能计算分类、裂隙结构形态的可视化、裂隙信息实时传输，以及构筑裂隙数据库的相关算法研究，当前还处于未成熟阶段。
因此，笔者从数据库概念和建立数据库基本方法出发，以实地取样的煤岩体裂隙重构模型为研究对象，构建了裂隙结构形态和力学属性的可视化数据库，为研究裂隙形态多样性识别提供一个新的方向和思路，不仅可以减少研究人员实地取样的繁重工作，提高科研成员的工作效率，对于需要重点考察煤岩裂隙结构形态的研究还能提供可实时获取结果的服务，满足资源共享的需求。
1. 数据库基本构架
1.1 数据库构筑的基本思路
数据库是以数据结构来进行数据读取、信息存储等操作的媒介，它是将现实数据先抽象为信息模型，再经过信息模型转换为数据模型的过程[15]。
建立数据库中各个实体之间的联系是建立数据库基本框架的必要条件。在煤岩体裂隙结构形态和力学属性的可视化数据库中，煤岩样品和裂隙是数据库中的两大实体。一个煤岩样品实体通常会存在多个裂隙实体，其关系是一对多的映射关系。所有煤岩样品实体的属性包括煤岩样品力学性质、赋存条件、地理位置、地质构造类型及加载条件等；裂隙实体的属性包含裂隙编号、张开位移、裂隙产状、应力强度因子、连续性、密度、体积、长宽和断裂类型等。两大实体之间的关系如图1所示。
![]()
图 1 煤岩样品与裂隙实体关系映射图Figure 1. Mapping diagram of relationship between coal and rock mass and cracks entities构筑裂隙数据库可提高研究含裂隙煤岩体损伤与断裂演化机制的效率。图2为数据库构筑流程，主体部分的搭建大致可以分为以下3个步骤：
![]()
图 2 裂隙结构形态和力学属性可视化数据库的构筑流程Figure 2. Logical graph of visual database construction of cracks structure morphology第1步，多源数据的智能获取。通过不同类型数据传递接口自动获取煤岩赋存勘测数据、模型参数与各类开采监测传感器等多源时空数据，统一存放入数据库用于可视化分析。其中，重构模型及裂隙的工程可用于实验及数值模拟研究；勘测信息中钻孔方位角、倾角和深度、煤岩开采水平及地质构造的产状等用于建立岩层面坐标系以确定岩层面与裂隙相对位置，判断局部地区裂隙赋存状况。
第2步，属性的智能分类与存储。数据库以建立裂隙属性存储码为出发点，整合裂隙、岩层及监测数据，完成裂隙智能计算与存储。通过融合多元线性回归、插值拟合与统计分析等多种算法，设计目标子程序担任裂隙形态分类属性计算任务；结合裂隙赋存条件与损伤断裂力学基本理论，计算各类断裂力学指标，确定煤岩体裂隙力学性能，完成裂隙力学属性与结构形态属性的高效自动分类与入库。计算完毕将裂隙各类信息汇总存储于具有动态更新的条形码与二维码中，以方便用户即时查阅。
第3步，形态可视化分析。基于查询索引的功能，数据库可用于属性图表分析及数据导出、自定义模型、裂隙面的点线面操作、力学属性的计算与结构形态数据的可视化。结合工程地质研究数据，基于大数据分析和协同管理等技术[16-18]，数据库建立专家知识库和裂隙动态演绎平台，通过实时监测裂隙演化数据，可实现煤样动态演化信息的协同和高效管理。通过自动判断、匹配裂隙发育类型，分析不同地质构造裂隙发育规律，以评估煤岩安全系数并提供最优支护开采方案。
通过以上流程，结合数据库计算裂隙特征参数，按照裂隙受力特征确定裂隙断裂类型、应力模型函数、应力强度因子及J积分等，计算断裂能、裂隙扩展方向与规律，可实现能量场、应力−应变场及裂隙形态演化过程的可视化分析。最终形成一个集模型构建与更新、智能计算与传输、工程应用于一体的综合动态可视化分析数据库。在实现各类时空信息可视化展示的同时，还可以根据工程研究进展不断对数据库进行扩充与完善，逐步提高数据库决策的准确性。
对于试验或数值模拟，数据库可提供多种多样的裂隙类型，解决了试验制备裂隙形态及裂隙组合的单一化问题，节省了取样、试样加工的时间。在数据库不断充实的情况下，可根据裂隙形态结构与发育特征推测出地质构造类型。
1.2 裂隙结构形态数据的获取
裂隙结构形态数据的获取是数据库可视化及属性智能识别计算的基础。
1.2.1 CT扫描
目前，CT扫描是广泛用于检测煤岩体内部结构及其组分的无损检测技术。付裕等[19]采用CT扫描重构煤样内部矿物结构模型，分析了煤样细观非均质性对其力学性能的影响。毛灵涛[20]等将工业CT扫描与DVC法相结合，分析了受载过程中红砂岩内部三维变形特征。宋红华等[21]基于CT扫描技术重构含原生裂隙的非均质煤岩模型，研究了不同方向加载时煤样非均质性对其破坏特征的影响。
CT扫描的本质是X射线断层扫描，利用煤岩体中不同成分对射线的吸收程度不同来工作。探测线经过扫描断层形成带有大量灰度信号的二维图像，采用阈值分割及形态学处理等算法将灰度图像不同组分进行划分。结合三维重构技术可提取出含有裂隙、煤基质及杂质的三维可视化模型。
1.2.2 核磁共振
核磁共振技术广泛应用于实验室各种岩样的孔裂隙分析。ZHAO等[22]采用核磁共振低温计量法测量了煤样孔径并对结果进行校准，实现了不同等级煤样孔隙结构形态的准确表征。李治豪等[23]采用核磁共振仪扫描了煤岩裂隙图像并重构出三维裂隙渗流模型用于数值模拟分析。周宏伟等[24]利用核磁共振成像技术实现了煤样注水过程中孔隙结构的可视化。
其工作原理是原子核在外加磁场与某特定射频磁场的作用下改变原有运动方式，产生共振吸收造成磁场发生变化[25-26]，接收器与扫描仪感知到此种变化并接收后，经过特定的空间编码形成带有灰度信号的核磁共振二维图。结合核磁共振弛豫机制与试验获得的图谱，可分析出煤岩体裂隙分布及结构属性特征参数。与CT扫描三维重建原理类似，采用核磁共振图像三维重建相关算法，可实现煤样内部的三维可视化重建。
1.2.3 聚焦离子束扫描电子显微镜（FIB-SEM）
基于扫描获得的岩样连续切片图，不仅可以进行矿物学成分分析，还可以通过重构三维数字岩样分析微纳米级尺度裂隙。杨永飞[27]等用三维聚焦离子束扫描图像作为训练样本完成了页岩数字岩芯的构建。JIANG等[28]采用聚焦离子束扫描电子显微镜和三维重构技术定性定量研究了页岩内部细观孔隙形态结构及其多相模型。LIU等[29]基于聚焦离子束扫描电子显微镜技术扫描图提取了煤样特征参数，构建了孔裂隙三维数字模型。
此技术的工作原理是聚焦离子束与待测样品相互作用激发出新的电子，探测器接收激发出来的二次电子进行二维成像。扫描精度更高，可识别微纳米尺度裂隙。
研究表明，数据库采用以上多种手段均可实现裂隙结构形态数据的获取及其三维可视化重现。同时，数据库可以基于三维裂隙展布形态与力学属性，结合有限元方法研究加载过程中含裂隙煤岩体的动态力学演化特征，并基于数据监测处理与理论计算，实现煤岩裂隙演化过程中应力、应变与能量场等数据的可视化。
2. 数据库特色功能
2.1 裂隙结构形态的智能识别
裂隙结构形态的智能识别是计算裂隙属性、判断裂隙类型及计算断裂指标的前提。
图3给出了数据库智能识别流程图，主要分为重构过程中裂隙结构形态的智能识别与裂隙分类属性计算时的智能识别2部分。基于边缘检测及形态学思想，采用编程算法将二维图像与背景区分。针对图像较暗难以区分裂隙和煤岩样品基质的情况，可通过提高图片的亮度以提升裂隙与基质对比度，获得具有明显裂隙形态的二维图像，结合三维重构软件实现裂隙的智能识别与三维重构。此方法操作简单，编程识别图片速度快，可较为准确地提取出裂隙的结构形态。
基于扫描与重构裂隙模型，可完成裂隙分类属性的智能识别。首先，通过数值分析软件提取三维数字模型的单元及坐标，将裂隙模型及二维断层图导入数据库。采用图像处理手段及分段自适应拟合曲线算法计算二维断层图像裂隙长宽、裂隙条数、分形维数及张开度等属性并确定裂隙所属类别。
其次，根据裂隙单元的位置关系判断其属于单裂隙还是多裂隙。基于最小二乘法将裂隙单元坐标进行拟合，计算出对应的线性回归面并算出裂隙面的法向量以确定局部坐标系下裂隙的倾角。多裂隙以分裂隙法向量及坐标关系确定多裂隙类型。根据线性拟合面对应的坐标与实际坐标的残差值确定裂隙平整度，通过三维裂隙所占单元等计算体积、长度、宽度及高度等属性参数。
最后，基于数据库中岩层及钻孔取样数据建立以岩层面为基准的整体坐标系。经空间直角坐标系间的转换，求出统一坐标系下裂隙面与岩层面法向量分别为n1和n2，根据二者间的夹角$ \alpha $确定两者的位置关系。同时，数据库也可根据煤岩加载情况与裂隙形态结构分析裂隙的受力特征，以确定裂隙类型并采用相应的断裂力学理论公式计算其应力强度因子。计算完毕，输出裂隙属性及分类结果。
裂隙智能识别分类算法流程如图4所示，其是以算法程序的方式实现二维和三维裂隙结构形态的智能识别，同时完成其属性的智能计算与分类。其中，m为循环读取过程中二维图像的号数；S、L、W分别为二维图像中实际裂隙的总面积、总周长和平均宽度；N、N0、N1分别为二维裂隙像素总数，二维图像像素总数，二维裂隙骨架像素总数； S0为二维图像实际总面积；（xu，yu）、（xe'，ye'，ze'）分别为二维图像中裂隙骨架化后的第u个像素坐标和三维裂隙第e个单元的中心坐标，（u和e取1，2，3，······）；Fj为根据裂隙坐标拟合的第j个曲线拟合方程式；p1，p2， p3分别为拟合曲线方程的二次项系数，一次项系数和常数项；B1、B2分别三维裂隙总单元数，三维煤样总单元数；d、h、V1分别为三维煤样实际底面直径，高度，总体积；V为三维裂隙实际总体积；z''与q1、q2、q3为三维裂隙面的拟合方程式及其各项系数；Qmin为所有裂隙三维坐标与拟合平面差的平方和的最小值。$ \left| {{\delta _e}} \right| $max为三维裂隙单元中心坐标与拟合平面距离的最大值的绝对值；n为以煤样底面为XY面，高为Z轴时裂隙法向量。
![]()
图 4 裂隙智能分类算法流程Figure 4. Logical graph of the intelligent classification algorithm of cracks2.2 分类属性的多角度判别
裂隙分类属性的多角度判别是数据库的第二项特色功能，基于此可判断裂隙在不同演化阶段的变化特征并预测裂隙扩展方向。
针对现有煤岩样品中提取的裂隙，结合已有裂隙分类标准，结合裂隙的扩展与受力情况，从裂隙几何属性、力学属性、稳定性及能量角度出发，将裂隙分为以下6类。图5为裂隙部分属性及分类标准。
1）根据裂隙面与岩层面的关系可分为水平裂隙、垂直裂隙、斜交裂隙三大类[30]。当裂隙面法向量所在直线与岩层面形成角度$ \alpha $=90°时为水平裂隙；当$ \alpha $=0°时为垂直裂隙；0°<$ \alpha $<90°时为斜交裂隙，在煤岩中分布最为广泛。
2）根据裂隙间的相对位置，可将裂隙分为单裂隙和多裂隙2种。两条裂隙距离最近的两点间的垂直距离D>d (d为煤样的直径)时划分为单裂隙。2条及以上裂隙间的距离0≤D≤d时分为多裂隙，多裂隙可分为交叉式、平行式和复合式。
3）根据裂隙面的平整度可分为平面型裂隙、曲面型裂隙和复合型裂隙。对每条单裂隙进行计算，裂隙各点坐标与其线性拟合面残差值的绝对值的最大值$ \left| {{\delta _e}} \right| $max≤2 mm时为平面型裂隙；$ \left| {{\delta _e}} \right| $max>2 mm分为曲面型裂隙。根据裂隙面与岩层面的夹角，平面型与曲面型裂隙又可分为斜交和水平裂隙两种。
4）结合形态学图像识别算法，计算二维扫描断层面裂隙长度和宽度的大小，按照长度及宽度所在范围分类，将其划分为A型、B型、C型和D型四种裂隙类型[31]。
5）根据裂隙受力情况分为张开型裂隙（I型）、滑移型裂隙（II型）、撕开型裂隙（III型）及复合型裂隙[32]。根据煤岩体的加载方式分析确定裂隙所受的应力特征，进而判断裂隙类型。结合煤岩体几何结构、裂隙形态尺寸形状与位置可计算其应力强度因子及应变能密度因子，得出临界应变能释放率的断裂判据，研究煤岩体的失稳破坏特征。
通过受力分析计算裂隙前缘各点的应力，采用裂隙尖端应力场的渐近解表达式，可推导出二三维裂隙尖端的应力强度因子为[32]：
$$ {K_{\text{I}}} = \mathop {\lim }\limits_{r \to 0} \sqrt {2{\text{π}} r} {\sigma _z}{|_{\theta = 0}} $$ (1) $$ {K_{{\text{II}}}} = \mathop {\lim }\limits_{r \to 0} \sqrt {2{\text{π}} r} {\tau _{zx}}{|_{\theta = 0}} $$ (2) $$ {K_{{\text{III}}}} = \mathop {\lim }\limits_{r \to 0} \sqrt {2{\text{π}} r} {\tau _{zy}}{|_{\theta = 0}} $$ (3) 式中：KI，KII，KIII分别为I型，II型，III型裂隙尖端应力强度因子，常用单位为MPa·m1/2；$ {\sigma _{\text{z}}} $、$ {\tau _{zx}} $、$ {\tau _{zy}} $分别为垂直于裂隙面的正应力，沿裂隙端部圆弧的法向和切向剪应力，MPa；x，y，z为直角坐标； $ （r,\theta ) $为裂隙尖端附近点的极坐标。
对于平面应力与平面应变状态，裂隙尖端应力强度因子的一般表达式为：
$$ {K_{\text{i}}} = Y\sigma \sqrt {\pi a} $$ (4) 式中：Ki为i型裂隙应力强度因子，其中i=I，II，III；Y为应力强度因子修正系数，与裂隙形状尺寸有关；a为裂隙尺寸，mm；$ \sigma $为名义应力，根据Ki确定其为正应力或剪应力，MPa。
主要研究单轴压缩条件下的煤岩体裂隙力学属性。假设式（4）中的Y取1，对于可以简化为平面应力与平面应变状态的裂隙，其应力强度因子K可表示为：
$$ {K_{\text{I}}} = - {\sigma _1}\sqrt {{\text{π}} a} {\cos ^2}\beta $$ (5) $$ {K_{{\text{II}}}} = {\tau _{\rm{e}}}\sqrt {{\text{π}} a} $$ (6) 式中：$ {\sigma _1} $为轴压，MPa；$ {\tau _{\rm{e}}} $为作用在裂隙面的有效剪应力，MPa；$ \beta $为裂隙与轴压加载方向的夹角，（°）。
在三轴压缩条件下，平面状态的裂隙应力强度因子的计算，应先求出作用在裂隙面上的正应力$ {\sigma _{\rm{n}}} $与有效剪应力$ {\tau _{\rm{e}}} $为
$$ {\sigma _{\rm{n}}} = \frac{1}{2}\left[ {\left( {{\sigma _1} + {\sigma _3}} \right) + \left( {{\sigma _1} - {\sigma _3}} \right)\cos\; 2\beta } \right] $$ (7) $$ {\tau _{\rm{e}}} = \frac{1}{2}\left( {{\sigma _1} - {\sigma _3}} \right)\sin \;2\beta - f{\sigma _{\rm{n}}} - {c_{\rm{w}}} $$ (8) 式中：$ {\sigma _3} $为围压，MPa；f为裂隙面的摩擦因数；cw为煤岩黏聚力，MPa。将式（7）和式（8）分别代入式（4）可得应力强度因子KI和KII。
对于三维裂隙，在单轴或三轴压缩条件下，其应力强度因子要通过式（1）—式（3）进行求解。
6）根据裂隙扩展的稳定性将其分为稳定扩展型裂隙和失稳扩展型裂隙，以预测裂隙扩展稳定程度。以基本断裂理论为依据，结合不同时段扫描裂隙的形态特征参数及受载条件计算断裂指标。当应力强度因子Ki≥KiC或能量释放率Gi≥G iC，（i=Ⅰ，Ⅱ，Ⅲ）时，为失稳型裂隙，反之则为稳定扩展型裂隙。其中，KiC代表i型裂隙尖端断裂韧度；Gi和GiC分别为i型裂隙尖端能量释放率和临界能量释放率，MN/m。
2.3 裂隙几何展布形态的可视化
裂隙几何结构展布形态的可视化是煤岩体裂隙结构形态和力学属性可视化数据库的第三项特色功能。
文中通过CT扫描和三维重构技术来实现裂隙几何结构展布形态的可视化功能。利用 CT扫描煤岩样品形成二维断层图像，基于编程算法对二维图像进行预处理，获得裂隙与背景具有明显对比度的图像。在三维重构软件中，根据裂隙和煤岩基质灰度范围的不同进行阈值选择，并以不同的颜色分别标记。对于难以区分的部分可以采用擦除颜色的方式来确认裂隙的边缘形态。经三维计算得出裂隙及煤岩样品其他成分的三维结构数值可视化模型。当钻孔垂直于岩层面且煤岩样品横截面与钻孔横截面平行时，煤样的XY面可看成岩层面，结合属性分类判别标准可确定裂隙类别。部分重构裂隙的几何展布形态可视化模型如图6所示。
![]()
图 6 裂隙几何结构形态可视化Figure 6. Visualization of cracks geometric structure and distribution三维可视化模型是裂隙形态特征参数计算与分类的依据。将重构的可视化模型导入数据库进行计算存储，用户在数据库中可基于裂隙特征参数查询筛选符合要求的裂隙属性及三维可视化模型。另外，用户可以通过交互功能查看空间裂隙的点坐标、位置分布及几何参数等并进行模型导出，利用相应的建模软件进行裂隙的切割、不同裂隙结构形态的重组、力学参数的赋值等。
2.4 裂隙属性的智能传输
裂隙属性信息的智能传输可对考察的裂隙区域提供实时获取结果的服务。用户扫描裂隙上传至数据库，调用数据库特色功能，完成裂隙属性的智能计算与自动判别，经条形码及二维码存储信息后上传云端并发送至用户移动端，完成裂隙信息的智能传输功能。数据库裂隙存储码生成流程如图7所示。
![]()
图 7 裂隙条形码和二维码生成流程Figure 7. Logical graph for the generation of crack barcodes and QR codes为初步区分大量裂隙，数据库计算裂隙特征参数后自动生成含有裂隙扫描时间、赋存条件及类型属性信息的特有的条形码序列号，形成具备存储及溯源裂隙研究信息的条形码，并为属性快速存储至二维码提供一个准确的属性索引标志。裂隙条形码序列号的编码规则如下：第1位代表煤或岩石，分别以“M”与“Y”表示，后2位代表扫描裂隙的煤矿代码，将不同地区的煤矿按照名称顺序编码，第5—6位是煤样所在地质构造类型，将不同地质构造类型编号。第8—15位是扫描时间，第17—19位是同一时间扫描的裂隙顺序码，最后一位是裂隙类型码，其中0代表多裂隙，1代表单裂隙。
基于裂隙的条形码序列号，经编码可生成黑白相间的条形码。条形码是对裂隙基本信息的存储，通过将裂隙信息与用户身份绑定，可在数据库追溯裂隙对应的研究进展。对应的条形码序列号相当于裂隙的ID，与下文中裂隙名称一一对应，通过裂隙ID可以查询裂隙的基本属性。将裂隙ID与数据库进行连接，经扫描解码即可查询出条形码所对应的裂隙信息。
二维码是对裂隙所有属性的综合描述，包含裂隙基本属性特征参数及动静态可视化研究。数据库连接二维码生成平台，以条形码序列号为标志进行裂隙各属性数据索引，存储在二维码平台创建的模板中。待信息编码整合后，将含有裂隙属性及研究信息的动态码上传至云端并发送至用户移动端供用户查询使用。其中，条形码与二维码可以根据数据库信息进行批量修改及实时更新。由此实现裂隙属性的智能传输。
3. 数据库应用实例
3.1 数据库构筑基础
数据库构筑过程主要包括裂隙重构及命名和功能实现2个方面。
3.1.1 数据库裂隙重构及命名
将取得的煤岩试样经标准加工后采用CT扫描裂隙断层，经Image J软件处理输出二维图像。采用Matlab编程智能识别处理二维灰度图，结合Mimics三维重构软件计算提取出煤岩试样内部结构的三维数字模型。
裂隙的命名采用“煤矿名称−煤岩性质及编号−裂隙编号”格式。前两个字母以取样煤矿首字母大写代表，第一个连字符后用于表示样品的性质，煤用“M”表示，岩石用“Y”表示，后面以两位数字格式表示样品的编号，第二个连字符后表示煤岩中提取出来的分裂隙，从上到下依次编号，以三位数字格式表示。如DA-M13-001代表大安山煤矿中编号为13的煤中提取的第一个裂隙。计算的部分裂隙属性及其分类见表1。
表 1 不同裂隙的参数及分类Table 1. Parameters and classification of different cracks裂隙编号 体积/mm3 角度/（°） 按平整度分 单/多裂隙 埋深/m
构造倾角/（°）构造类型 裂隙二维码 DA-M13-001 196.70 0 平面型 单裂隙 400
10~90倒转背斜 
DT-M20-001 81.62 60 曲面型 单裂隙 800
70~90向斜 
DT-M37-001 118.29 68 平面型 单裂隙 800
70~90向斜 
CJ-M01-001 162.4 30 曲面型 单裂隙 680
8~15背斜 
CJ-M01-002 40.05 28 平面型 单裂隙 680
8~15背斜 
CJ-M01-003 47.65 50 平面型 单裂隙 680
8~15背斜 
DA-M01-004 506.63 — 平面型 多裂隙 400
50~90倒转向斜 
CJ-M02-001 1704.93 — 曲面型 多裂隙 680
8~15背斜 
CJ-M03-001 687.79 — 曲面型 多裂隙 680
8~15背斜 
3.1.2 数据库功能实现基础
煤岩体裂隙结构形态和力学属性的可视化数据库应具有编辑、录入和筛选等数据访问与智能计算功能，而实现这些功能的基础是接口的转化。
接口的类型众多，数据库通过系统操作接口来实现数据的编辑、录入和查询等[33-34]。以数据库查询模块[35]为例，接口的作用是将服务器与用户端进行连接供用户使用操作数据库。裂隙可视化数据库是基于Matlab与Visual C++进行的数据库构筑，基于Matlab GUI界面编程进行数据库各个模块的设计与优化。利用C++设计数据库采集接口实现与外部软件的连接进行数据与模型的更新导入；以Matlab编程实现数据库中裂隙面的智能识别计算；数据导入、筛选与可视化分析等功能通过编辑设计控件按钮实现；结合C#编译接口实现不同库表及模块的交互、属性的自动分类和数据存储等功能；连接条形码与二维码生成平台实现含有动态信息的裂隙专属二维码与条形码的自动生成等，以便用户随时随地进行相关信息查询。
3.2 数据库案例分析
为了进一步探讨数据库的操作与应用，以现有的北京昊华能源股份有限公司大安山煤矿取得的煤岩样品为例，阐述数据索引入库与管理模块的特点，进行裂隙可视化数据库的初探。
3.2.1 地理位置及数据来源
北京昊华能源股份有限公司大安山煤矿位于北京城正西距50 km。大安山煤矿的矿井范围内赋存大量逆冲断层和倒转褶皱等复杂地质构造，如图8所示[36]。井田构造主要是以西南−东北向复式褶曲构造形态为主，背斜多呈紧闭状，向斜较为宽缓，伴有一定量的倾向、走向、断层，次级构造十分发育。其中煤岩内部裂隙结构形态具有代表性，如裂隙扩展方向及角度具有多样性、裂隙体积变化多样、简单与复杂形态结构的裂隙有时可共存且可能为层状与交叉等多种位置形式、主裂隙与次生裂隙分布明显、具有平面曲面型裂隙等，总体来说裂隙形态丰富多样且有一定的规律可循。因此，对裂隙进行重构提取，进行数据库各模块的构筑。
3.2.2 裂隙检索及数据管理
数据库具有资源共享性的特点。如图9所示，在数据库前端界面，用户通过账号密码授权登录进入数据分析与管理模块。图10为初步设计的数据库各模块窗口。
数据分析模块主要是对数据库裂隙的查询分析。具有筛选裂隙属性类型与数据可视化分析的作用，可查看裂隙的研究进展，为后续研究提供基础。数据库提供多个裂隙分类词条供用户选择，如选择表面平整度为曲面、裂隙面与岩层面所成夹角在60°到80°之间的单裂隙进行可视化分析。符合条件的裂隙会以列表形式体现，选择相应的裂隙可查看其类型、研究次数及可视化模型等详细属性，对于满足要求的裂隙可批量导出分析。
数据管理模块主要是对数据的导入、属性计算与数据更新，是保证数据库质量好坏的核心。对于大量裂隙，可批量导入计算与存储。数据库根据导入的图像、勘测数据和裂隙模型计算裂隙的体积、张开度、角度及断裂类型等属性并进行自动分类，同时结合煤岩体的荷载及边界条件进行受力分析，计算应力强度因子并分别存放于相应的库表中。对数据库中原有数据的更新，可通过裂隙编号的索引对裂隙的属性进行修改、编辑、录入与保存等。
4. 结 论
1）从裂隙结构形态的获取及智能识别、裂隙智能分类与存储、形态结构的可视化分析与实现三个方面构筑了裂隙可视化数据库的基本框架，实现了裂隙结构形态智能识别、分类属性多角度判别、几何结构展布形态可视化、基于条形码与二维码的裂隙属性智能传输等特色功能。
2）数据库可基于裂隙几何结构形态的坐标，采用图像拟合等算法实现裂隙多角度分类属性的智能计算与分类；数据库与条形码及二维码生成平台相连接，可实时获取裂隙计算结果。调用数据库特色功能，将计算的裂隙力学属性及形态结构等信息编码整合后存储至实时更新的条形码与二维码中并发送至用户端，实现了裂隙属性的智能传输。
3）以北京昊华能源股份有限公司大安山煤矿为地质背景，进行了煤岩样品的可视化重构与统一命名，结合模型计算方法对裂隙的角度、体积和形态类型等属性参数进行了计算与分类。基于数据库构筑流程设计了数据库登录模块、数据筛选与分析模块和数据管理模块，实现了数据库内部模块的基本功能。
-
![]()
图 2 矿井综合水文地质结构示意
Figure 2. Comprehensive hydrogeological schematic histogram of the coal mine
![]()
图 3 高盐矿井水的减排治理技术思路与框架
Figure 3. Technical system for emission reduction and treatment of high salt mine water
![]()
图 4 底板涌水封堵减量技术示意
Figure 4. Schematic of water inrush plugging and reduction technology of floor
![]()
图 5 矿井水深层回灌顶板含水层和底板含水层示意
Figure 5. Schematic of the roof aquifer and floor aquifer of the deep reinjection of mine water
![]()
图 7 不同减量治理方案下水位恢复预测趋势
Figure 7. Prediction trend of water level recovery under different reduction control schemes
![]()
图 9 采空区调蓄系统水化学成分对比
Figure 9. Comparison of chemical composition of water in goaf regulation and storage system
表 1 近5年矿井涌水量情况
Table 1 Water inflow of mine in recent 5 years
年份 2018 2019 2020 2021 2022 月均涌水量/(m3·h−1) 769 748 671 695 540 
下载: 导出CSV
表 2 13采区西翼四灰含水层疏降情况预计
Table 2 The dewatering of four limestone aquifer in the west wing of 13 mining area
时间t/d 总疏降强度/(m3·h−1) 总降深/m 1 200 132.8 2 141.6 5 146.8 7 153.3 14 157.6 20 170.9 40 179.8 60 184.9 $ \vdots $ $ \vdots $ 120 193.7 180 198.9
下载: 导出CSV
表 3 奥灰含水层各指标情况
Table 3 Indicators of Ordovician limestone aquifer
含水层 组 与9煤
距离/m层厚/m 渗透系数/
(m·d−1)奥灰含水层 阁庄组 207 (105.5～121.16 m)/113.33 0.32～4.83 马家沟组 198.53 0.55～1.70
下载: 导出CSV
表 4 回灌后奥灰水位抬升预测
Table 4 Ordovician limestone water level rise after reinjection
时间/d 回灌量/
(m3·h−1)水位抬升/m r=200 m处 r=500 m处 1 200 18.78 15.86 2 19.88 16.97 5 21.34 18.42 7 21.87 18.96 10 22.44 19.52 20 23.54 20.63 40 24.64 21.73 60 25.29 22.37 120 26.39 23.47 180 27.03 24.12 300 27.84 24.93 480 28.59 25.68 600 28.95 26.03 1000 29.76 26.84 1500 30.40 27.49
下载: 导出CSV
表 5 不同含水层水质指标对比情况
Table 5 Comparison of water quality indicators of different aquifers
序号 水样来源 阳离子浓度/(mg·L−1) 阴离子浓度/(mg·L−1) TDS/(mg·L−1) K++Na+ Ca2+ Mg2+ Cl− SO4 2− HCO3 − 1 四灰水 361.65～441.05 199.80～465.53 49.68～110.98 184.04～324.25 1 047.93～1 893.36 169.94～203.44 1 926～3384 2 矿井水 468.33～783.32 210.76～473.86 78.64～188.95 252.7～336.6 1 638.68～2 296.21 155.29～858.69 3 112～4418 3 奥灰水 579.41 731.38 143.91 815.84 2 066.29 219.90 4447
下载: 导出CSV
-
[1] 武 强,申建军,王 洋. “煤–水”双资源型矿井开采技术方法与工程应用[J]. 煤炭学报,2017,42(1):8−16. WU Qiang,SHEN Jianjun,WANG Yang. Mining techniques and engineering application for “Coal-Water”dual-resources mine[J]. Journal of China Coal Society,2017,42(1):8−16.
[2] 马莲净,赵宝峰,吕玉广,等. 煤矿矿井水水化学形成作用与灌溉适宜性评价[J]. 环境化学,2023:1–10. MA Lianjing,ZHAO Baofeng,LYU Yuguang. Hydrogeochemical processes and irrigation applicability of mine water in coal mine[J]. Environmental Chemistry,2023:1–10.
[3] YANG L,TIANTIAN W,JIAN Y. Evaluating the quality of mine water using hierarchical fuzzy theory and fluorescence regional integration[J]. Mine Water and the Environment,2019,38(2):243−251. doi: 10.1007/s10230-018-0567-4
[4] 王 皓,董书宁,尚宏波,等. 国内外矿井水处理及资源化利用研究进展[J]. 煤田地质与勘探,2023,51(1):222−236. WANG Hao,DONG Shuning,SHANG Hongbo, et al. Domestic and foreign progress of mine water treatment and resource utilization[J]. Coal Geology & Exploration,2023,51(1):222−236.
[5] 孙亚军,陈 歌,徐智敏,等. 我国煤矿区水环境现状及矿井水处理利用研究进展[J]. 煤炭学报,2020,45(1):304−316. SUN Yajun,CHEN Ge,XU Zhimin, et al. Research progress of water environment,treatment and utilization in coal mining areas of China[J]. Journal of China Coal Society,2020,45(1):304−316.
[6] TONG L,FAN R,YANG S,et al. Development and Status of the Treatment Technology for Acid Mine Drainage[J]. Mining,Metallurgy & Exploration,2021,38(1):315–327.
[7] 曾繁富,左明星,宋洪柱,等. 乌审旗一带刘家沟组作为高矿化度矿井水回灌目的层的可行性分析[J]. 煤炭与化工,2020,43(11):59−62. ZENG Fanfu,ZUO Mingxing,SONG Hongzhu, et al. Feasibility analysis of the Liujiagou Group in the Wushen Banner area as a target layer for water recharge in highly mineralized mine[J]. Coal and Chemical Industry,2020,43(11):59−62.
[8] 顾大钊,李 庭,李井峰,等. 我国煤矿矿井水处理技术现状与展望[J]. 煤炭科学技术,2021,49(1):11−18. GU Dazhao,LI Ting,LI Jingfeng, et al. Current status and prospects of coal mine water treatment technology in China[J]. Coal Science and Technology,2021,49(1):11−18.
[9] 海 晶,孙玉芳,马永祥. 宁夏石炭井矿区矿井涌水量与充水因素关系研究[J]. 地下水,2019,41(6):12−14. HAI Jing,SUN Yufang,MA Yongxiang. Study on the relationship betwenen well inflow and water filling factors in carbonferous mine area in ningxia[J]. Ground Water,2019,41(6):12−14.
[10] 曾一凡,武 强,赵苏启,等. 我国煤矿水害事故特征、致因与防治对策[J]. 煤炭科学技术,2023,51(7):1−14. doi: 10.13199/j.cnki.cst.2023-0500 ZENG Yifan,WU Qiang,ZHAO Suqi, et al. Characteristics, causes, and prevention measures of coal mine water hazard accidents in China[J]. Coal Science and Technology,2023,51(7):1−14. doi: 10.13199/j.cnki.cst.2023-0500
[11] 张 保,郝秀强,李会强. 神东矿区采区式煤矿地下水库系统设计研究[J]. 煤炭工程,2022,54(12):1−6. ZHANG Bao,HAO Xiuqiang,LI Huiqiang. System design of mining area-coal mine underground reservoir in Shendong mining[J]. Coal Engineering,2022,54(12):1−6.
[12] 智国军,鞠金峰,刘 润,等. 水岩相互作用对煤矿地下水库水质影响机理研究[J]. 采矿与安全工程学报,2022,39(4):779−785. ZHI Guojun,JU Jinfeng,LIU Run, et al. Water-rock interaction and its influence on water quality in the underground reservoir[J]. Journal of Mining & Safety Engineering,2022,39(4):779−785.
[13] 董书宁,王 海,黄选明,等. 基于保障生态地下水位的露天煤矿主动保水技术研究[J]. 煤炭科学技术,2021,49(4):49−57. DONG Shuning,WANG Hai,HUANG Xuanming, et al. Research on active water conservation technology in open-pit coal mine based on ecological protection groundwater level[J]. Coal Science and Technology,2021,49(4):49−57.
[14] 武 文. 煤矿矿井水处理站排放达地表水Ⅲ类标准改造工艺选择[J]. 能源与节能,2023(2):156−158. WU Wen. Selection of reforming process for discharging surface water up to class Ⅲ standard in mine water treatment stations[J]. Energy and Energy Conservation,2023(2):156−158.
[15] 李 鑫,孙亚军,陈 歌,等. 高矿化度矿井水深部转移存储介质条件及影响机制[J]. 煤田地质与勘探,2021,49(5):17−28. LI Xin,SUN Yajun,CHEN Ge, et al. Medium conditons and influence mechanism of high salinity mine water transfer and storage by deep well recharge[J]. Coal Geology & Exploration,2021,49(5):17−28.
[16] 赵春虎,杨 建,王世东,等. 矿井水深层回灌过程量质耦合模拟分析[J]. 煤田地质与勘探,2021,49(5):36−44. ZHAO Chunhu,YANG Jian,WANG Shidong, et al. Coupling simulation of groundwater dynamics and solute transfer in the process of deep reinjection of mine water[J]. Coal Geology & Exploration,2021,49(5):36−44.
[17] 顾大钊,李井峰,曹志国,等. 我国煤矿矿井水保护利用发展战略与工程科技[J]. 煤炭学报,2021,46(10):3079−3089. GU Dazhao,LI Jingfeng,CAO Zhiguo, et al. Technology and engineering development strategy of water protection and utilization of coal mine in China[J]. Journal of China Coal Society,2021,46(10):3079−3089.
[18] 汪 佩,李 鑫,张金陵,等. 城郊煤矿底板太灰疏降流场演化的数值模拟与应用[J]. 煤炭科技,2021,42(4):38−43. WANG Pei,LI Xin,ZHANG Jinling, et al. Numerical simulation and application of the floor evolution of flow field in Chengjiao coal mine[J]. Coal Science & Technology Magazine,2021,42(4):38−43.
[19] 戴 磊,段李宏,杨 斌,等. 顺和煤矿工作面底板太灰水动态引流方案设计应用[J]. 矿业安全与环保,2022,49(6):106−111. DAI Lei,DUAN Lihong,YANG Bin, et al. Design and application of dynamic drainage scheme for Taiyuan formation limestone aquifer in working face floor of Shunhe coal mine[J]. Mining Safety & Environmental Protection,2022,49(6):106−111.
[20] 孙亚军,徐智敏,李 鑫,等. 我国煤矿区矿井水污染问题及防控技术体系构建[J]. 煤田地质与勘探,2021,49(5):1−16. SUN Yajun,XU Zhimin,LI Xin, et al. Mine water drainage pollution in China’s coal mining areas and the construction of prevention and control technical system[J]. Coal Geology & Exploration,2021,49(5):1−16.
[21] 刘毅涛,王国文,龙陆军,等. 矿井水深部回灌运移机理及扩散规律研究[J]. 中国煤炭地质,2021,33(12):36−41. LIU Yitao,WANG Guowen,LONG Lujun, et al. Study on migration mechanism and diffusion pattern in mine water deep influsion[J]. Coal Geology of China,2021,33(12):36−41.
[22] 郭再峰. 煤矿矿井水净化处理与回灌治理技术[J]. 探矿工程(岩土钻掘工程),2011,38(12):79−82. GUO Zaifeng. Water decontamination and recharge treatment in coal mine[J]. Drilling Engineering,2011,38(12):79−82.
[23] 孙亚军,李 鑫,冯 琳,等. 鄂尔多斯盆地煤–水协调开采下矿区水资源异位回灌–存储技术思路[J]. 煤炭学报,2022,47(10):3547−3560. SUN Yajun,LI Xin,FENG Lin, et al. Technical thinking on ectopic injection and storage of mine area water resources under the coordinated exploitation of coal and water background in Ordos basin[J]. Journal of China Coal Society,2022,47(10):3547−3560.
[24] 赵彩凤,杨 建. 高矿化度矿井水地下回灌对毛乌素沙漠地下水水质的影响[J]. 煤矿安全,2018,49(3):29−32. ZHAO Caifeng,YANG Jian. Influence of groundwater recharge with high salinity mine water on groundwater quality in Mu Us desert[J]. Safety in Coal Mines,2018,49(3):29−32.
[25] 姜 素,黄敬军,徐智敏,等. 徐州煤矿采空区地下水库建库可行性研究[J]. 水文地质工程地质,2018,45(5):17−23. JIANG Su,HUANG Jingjun,XU Zhimin, et al. Feasibility study of the construction of groundwater reservoir in the goaf of the Xuzhou Coal Mines[J]. Hydrology & Engineering Geology,2018,45(5):17−23.
[26] 袁慧卿,徐智敏,罗武贤,等. 张双楼煤矿井下四灰涌水减量治理与节能减排技术[J]. 煤炭科技,2022,43(4):157−163. YUAN Huiqing,XU Zhimin,LUO Wuxian, et al. Reduction treatment and energy saving and emission reduction technology of underground four limestone water inflow in Zhangshuanglou coal mine[J]. Coal Science & Technology Magazine,2022,43(4):157−163.
-
期刊类型引用(1)
1. 李小辉，邢楷，彭群. 含水纵向裂隙煤岩破裂特征与演化规律研究. 矿业研究与开发. 2024(07): 191-196 . 
百度学术
其他类型引用(2)
计量
- 文章访问数: 110
- HTML全文浏览量: 19
- PDF下载量: 38
- 被引次数: 3
