Study on influence of mining ground fissures on soil anti-erodibility in northern Shaanxi coal mining area
-
摘要:
采动地裂缝作为黄河流域中游陕北煤矿区最突出且典型的采动损害问题，引发的水土流失效应已经不容忽视。为了研究采动地裂缝发育对周边土壤抗蚀性的影响，以陕北煤矿区内宽度分别为0～10，10～20，20～30 cm的采动地裂缝为研究对象，采集周边水平距离80 cm以内，垂直深度40 cm以浅的土壤，测定了土壤水稳性团聚体、微团聚体、无机黏粒、有机黏粒、物理性质类等14个国内外学者普遍关注的用于量化表征土壤抗蚀性的指标，采用层次分析法、敏感性分析和因子分析相结合的方法甄选出陕北采动地裂缝发育区土壤抗蚀性的重要量化指标，构建了陕北采动地裂缝发育区土壤抗蚀性综合指数模型。结果表明：①综合层次分析法、敏感性分析和因子分析，甄选确定>0.25 mm水稳性团聚体含量、平均重量直径、<0.001 mm细黏粒含量、土壤团聚度、土壤有机质含量5个指标作为量化表征陕北采动地裂缝发育区土壤抗蚀性的重要指标；②采动地裂缝会降低周围土壤的5个土壤抗蚀性重要指标，其中>0.25 mm水稳性团聚体含量变化最为明显；③根据土壤抗蚀性重要指标，基于因子分析原理构建了陕北采动地裂缝发育区土壤抗蚀性综合指数模型；④采动地裂缝会降低周围土壤的抗侵蚀综合能力，且该效应随着裂缝宽度的增大和水平距离的减小而增强，当距采动地裂缝的水平距离超过170 cm时，采动地裂缝对周围土壤抗侵蚀综合能力的负效应基本消失，其可作为陕北采动地裂缝发育区土壤侵蚀防控的关键区域。研究结果可为黄河中游陕北矿区水土流失的精准防控提供科学依据。
Abstract:As the most prominent and typical mining damage problem in the coal mining area of northern Shaanxi in the middle reaches of the Yellow River Basin, the soil erosion effect caused by mining ground fissures can not be ignored. In order to study the influence of the development of mining ground fissures on the anti-erodibility of the surrounding soil, the mining ground fissures with widths of 0−10, 10−20 and 20−30 cm in the northern Shaanxi coal mining area were taken as the research object. The soil with a horizontal distance of less than 80 cm and a vertical depth of 40 cm was collected. Fourteen indicators widely concerned by domestic and foreign scholars for quantitatively characterizing soil erosion resistance, including soil-water stable aggregates, micro-aggregates, inorganic clays, organic clays, and physical properties, were measured. The important quantitative indexes of soil anti-erodibility in the developed area of ground fissures in northern Shaanxi were selected by the methods of analytic hierarchy process, sensitivity analysis and factor analysis, and the comprehensive index model of soil anti-erodibility in the developed area of ground fissures in northern Shaanxi was constructed. The results show that:①Comprehensive analytic hierarchy process, sensitivity analysis and factor analysis, the five indexes of >0.25 mm water-stable aggregate content, mean weight diameter, <0.001 mm fine clay content, soil reunion degree and organic matter content were selected as the key indexes for quantitative characterization of soil anti-erodibility in the mining ground fissure development area of northern Shaanxi.②Mining ground fissures will reduce the five key indicators of soil anti-erodibility in the surrounding soil, among which the content of >0.25 mm water-stable aggregates changes most obviously.③According to the important indexes of soil anti-erodibility, based on the principle of factor analysis, a comprehensive index model of soil anti-erodibility in the developed area of mining ground fissures in northern Shaanxi was established.④The mining ground fissures will reduce the comprehensive anti-erosion ability of the surrounding soil, and this effect increases with the increase of the crack width and the decrease of the horizontal distance. When the horizontal distance from the mining ground fissures exceeds 170 cm, the negative effect of the mining ground fissures on the comprehensive anti-erosion ability of the surrounding soil basically disappears, which can be used as a key area for soil erosion prevention and control in the mining-induced ground fissure development area in northern Shaanxi. The research results can provide scientific basis for accurate prevention and control of soil erosion in northern Shaanxi mining area in the middle reaches of the Yellow River.
-
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）以北京昊华能源股份有限公司大安山煤矿为地质背景，进行了煤岩样品的可视化重构与统一命名，结合模型计算方法对裂隙的角度、体积和形态类型等属性参数进行了计算与分类。基于数据库构筑流程设计了数据库登录模块、数据筛选与分析模块和数据管理模块，实现了数据库内部模块的基本功能。
-
![]()
图 3 不同宽度采动地裂缝周边土壤水稳性团聚体含量变化
Figure 3. Variation of soil water-stable aggregate content around ground fissures with different widths
![]()
图 4 不同宽度采动地裂缝周围土壤平均重量直径变化
Figure 4. Variation of soil a mean weight diameter around ground fissures with different widths
![]()
图 5 不同宽度采动地裂缝周边土壤无机黏粒含量变化
Figure 5. Variation of soil inorganic clay content around ground fissures with different widths
![]()
图 6 不同宽度采动地裂缝周边土壤团聚度变化
Figure 6. Variation of soil aggregation around ground fissures with different widths
![]()
图 7 不同宽度采动地裂缝周边土壤有机质变化
Figure 7. Variation of soil organic matter around groundfissures with different widths
![]()
图 8 不同宽度采动地裂缝周边土壤抗蚀性综合指数变化
Figure 8. Variation of comprehensive index of soil anti-erodibility around ground fissures with different widths.
表 1 土壤抗蚀性指标
Table 1 Soil erodibility indexes
属性 编号 评价指标 选取依据 文献来源 水稳性团聚体类 X1 >0.25 mm水稳性团聚体含量 对维持土壤肥力有重要作用，其值越大，土壤的抗蚀性越大
反映水稳性团聚体的分布特征
反映土壤团聚体的稳定性特征，其值越大，团聚体稳定性越强[23-24] X2 >0.5 mm水稳性团聚体含量 X3 平均重量直径 无机黏粒类 X4 <0.001 mm细黏粒含量 细颗粒含量的多少与胶结程度有直接联系，胶结程度增大，土壤抗蚀性增强 [23-25] X5 <0.01 mm物理性黏粒含量 X6 <0.05 mm粉黏粒含量 X7 结构性颗粒指数 反映土壤质地的黏重情况，黏粒占比较大，抗蚀性较强 微团聚体类 X8 分散率 反映土壤在水中被破坏的程度，其值越小，土壤抗蚀性越强 [24-25] X9 分散系数 X10 团聚状况 反映了土壤团聚化程度，其值越大，则土壤抗蚀性越强 X11 团聚度 有机胶体类 X12 有机质 能够促进土壤团粒结构的形成，增强土壤抗蚀性 [24-26] 物理性质类 X13 含水率 直接影响土壤抵抗侵蚀的能力，反映土壤水土保持功能 [27-28] X14 容重 反映土壤紧实程度，通过改变土壤的孔隙率、水分入渗性，影响土壤抗蚀性 
下载: 导出CSV
表 2 层次分析法确定指标权重结果
Table 2 Analytic hierarchy process for determining indicator weights results
抗蚀性指标 权重 X1 0.2882 X2 0.0873 X3 0.1586 X4 0.0719 X5 0.0198 X6 0.0198 X7 0.0337 X8 0.0186 X9 0.0249 X10 0.0201 X11 0.0388 X12 0.1548 X13 0.0390 X14 0.0246
下载: 导出CSV
表 3 土壤抗蚀性指标敏感性分析结果
Table 3 Sensitivity analysis results of soil anti-erodibility indexes
评价指标 均值 极差 标准差 相对极差 变异系数/% X1 15.55 36.86 6.62 2.37 42.60 X2 9.81 27.71 5.40 2.82 54.98 X3 0.36 0.64 0.10 1.79 27.71 X4 2.63 2.08 0.35 0.79 13.16 X5 27.26 12.97 2.50 0.48 9.16 X6 77.59 15.60 3.38 0.20 4.35 X7 3.50 2.14 0.34 0.61 9.59 X8 76.56 17.16 4.34 0.22 5.67 X9 0.52 0.43 0.12 0.83 22.08 X10 12.98 6.89 1.93 0.53 14.89 X11 36.85 27.98 6.80 0.76 18.44 X12 3.76 5.16 1.38 1.37 36.71 X13 5.80 6.43 1.61 1.11 27.73 X14 1.45 0.18 0.04 0.13 3.15
下载: 导出CSV
表 4 土壤抗蚀性指标因子分析
Table 4 Factor analysis of soil anti-erodibility indexes
抗蚀性指标 公因子 F1 F2 X1 0.094 0.982 X2 0.162 0.968 X3 0.015 0.995 X4 0.987 −0.01 X5 0.965 −0.03 X6 0.979 −0.15 X7 0.971 0.07 X8 −0.918 −0.23 X9 −0.945 −0.27 X10 0.963 0.145 X11 0.989 0.028 X12 −0.089 0.985 X13 0.750 −0.61 X14 0.928 −0.24 注：表中数字代表公因子在各抗蚀性指标上的因子载荷系数。
下载: 导出CSV
表 5 土壤抗蚀性重要指标筛选
Table 5 Screening of important indexes of soil anti-erodibility
排序 第1 第2 第3 第4 第5 第6 第7 第8 层次分析 X1 X3 X12 X2 X4 X13 X11 X7 敏感性分析 X1 X3 X12 X13 X9 X4 X11 X10 因子分析 X3 X11 X4 X12 X1 X6 X7 —
下载: 导出CSV
表 6 不同宽度采动地裂缝周边土壤抗蚀性重要指标测定结果
Table 6 Determination results of important indexes of soil anti-erodibility around ground fissures with different widths
土层深度/cm 裂缝宽度/cm 水平距离/cm >0.25 mm水稳性
团聚体含量/%平均重量
直径/mm<0.001 mm细黏粒
含量/%团聚度/% 有机质含量/
(g·kg−1)0～20 0～10 20 20.52±1.03e 0.45±0.03d 2.33±0.12d 31.89±1.60d 5.35±0.27d 40 23.25±1.16d 0.51±0.03c 2.53±0.13cd 36.13±1.81c 5.54±0.28cd 60 27.67±1.39c 0.54±0.03bc 2.82±0.15bc 41.09±2.06b 5.93±0.30c 80 32.71±1.64b 0.58±0.03b 2.94±0.15b 43.44±1.20ab 6.52±0.33b 10～20 20 18.35±0.92f 0.41±0.02e 2.19±0.07d 26.97±1.35e 4.31±0.22e 40 19.45±0.97ef 0.42±0.02de 2.39±0.12d 31.32±1.57d 4.54±0.23e 60 20.98±1.05e 0.45±0.03d 2.67±0.14c 38.31±1.92c 5.12±0.26d 80 24.23±1.21d 0.49±0.03cd 2.81±0.14bc 39.78±1.99bc 5.46±0.28d 20～30 20 13.24±0.67h 0.33±0.02f 1.74±0.09e 22.85±1.14f 3.54±0.18f 40 14.29±0.72gh 0.34±0.02f 2.23±0.12d 27.58±1.38e 4.15±0.21e 60 15.32±0.77g 0.36±0.02f 2.40±0.14d 32.14±1.61d 4.48±0.23e 80 17.41±0.87f 0.39±0.02ef 2.65±0.14c 34.68±1.74cd 5.07±0.26d CK — 44.47±2.23a 0.88±0.05a 3.21±0.16a 45.02±1.25a 7.15±0.35a 20～40 0～10 20 7.76±0.39g 0.26±0.01de 2.81±0.14cd 39.31±1.97c 2.48±0.12cd 40 9.87±0.49ef 0.27±0.01d 2.94±0.15c 44.13±2.21b 2.66±0.13c 60 12.12±0.61d 0.30±0.01c 3.15±0.16bc 47.78±2.39ab 2.92±0.15b 80 15.88±0.79b 0.33±0.02b 3.22±0.16b 48.75±2.44ab 3.25±0.16a 10～20 20 7.68±0.38g 0.26±0.01de 2.67±0.13d 34.13±1.71d 2.41±0.12d 40 9.07±0.45f 0.26±0.01de 2.68±0.13d 37.78±1.89cd 2.48±0.12cd 60 10.19±0.51e 0.28±0.01cd 2.93±0.15c 42.11±2.11bc 2.59±0.13cd 80 14.37±0.72c 0.30±0.02c 3.02±0.15bc 45.23±2.26b 2.74±0.14bc 20～30 20 7.61±0.38g 0.24±0.01e 2.08±0.10d 27.13±1.36e 1.99±0.10e 40 8.92±0.45f 0.24±0.01e 2.52±0.13e 34.84±1.74de 2.16±0.11e 60 9.51±0.48ef 0.28±0.01cd 2.71±0.14cd 39.39±1.97c 2.26±0.11de 80 12.82±0.64d 0.29±0.01cd 2.81±0.14cd 41.11±2.06bc 2.34±0.12de CK — 17.63±0.88a 0.36±0.02a 3.82±0.19a 50.83±2.54a 3.23±0.16a 注：相同土壤深度下同一指标数据后不同字母表示p<0.05水平上差异显著。
下载: 导出CSV
表 7 不同宽度采动地裂缝周围土壤$ {{{Y}}_{\text{K}}} $与水平距离的量化关系
Table 7 Quantitative relationship between soil$ {{{Y}}_{\text{K}}} $and horizontal distance around ground fissures with different widths
裂缝宽度/cm 土壤深度/cm 拟合公式 相关系数 0～10 0～20 ${{{Y}}_{\text{K}}} = 6.78 - 6.19{{\text{e}}^{ - {{d}}/213.64}}$ 0.99 20～40 ${{{Y}}_{\text{K}}} = 6.39 - 5.49{{\text{e}}^{ - {{d}}/283.70}}$ 0.99 10～20 0～20 ${{{Y}}_{\text{k}}} = 16.69 - 16.43{{\text{e}}^{ - {{d}}/713.60}}$ 0.97 20～40 ${{{Y}}_{\text{K}}} = - 0.49 + 1.19{{\text{e}}^{ - {{d}}/713.60}}$ 0.96 20～30 0～20 ${{{Y}}_{\text{K}}} = 39.41 - 38.80{{\text{e}}^{ - {{d}}/1661.51}}$ 0.98 20～40 ${{{Y}}_{\text{K}}} = 10.88 - 11.13{{\text{e}}^{ - {{d}}/454.54}}$ 0.96 注：d为距裂缝的水平距离，cm。
下载: 导出CSV
-
[1] 王双明,申艳军,宋世杰,等. “双碳”目标下煤炭能源地位变化与绿色低碳开发[J]. 煤炭学报,2023,48(7):2599−2612. WANG Shuangming,SHEN Yanjun,SONG Shijie,et al. Change of coal energy status and green and low-carbon development under the “dual carbon” goal[J]. Journal of China Coal Society,2023,48(7):2599−2612.
[2] 宋世杰,孙 涛,郑贝贝,等. 陕北黄土沟壑区采煤沉陷对黄土坡面形态的影响及土壤侵蚀效应[J]. 煤炭科学技术,2023,51(2):422−435. SONG Shijie,SUN Tao,ZHENG Beibei,et al. Effect of coal mining subsidence on loess slope morphology and soil erosion in loess gully region of northern Shaanxi[J]. Coal Science and Technology,2023,51(2):422−435.
[3] 中国政府网. 习近平主持召开深入推动黄河流域生态保护和高质量发展座谈会并发表重要讲话[EB/OL]. (2021−10−22). https://www.gov.cn/xinwen/2021-10/22/content_5644331.htm [4] 宋世杰,杜 麟,王双明,等. 陕北采煤沉陷区不同沉陷年限黄土坡面土壤可蚀性的变化规律[J]. 煤炭科学技术,2022,50(2):289−299. SONG Shijie,DU Lin,WANG Shuangming,et al. Variation of soil erodibility on loess slope under various subsidence years in coal mining subsidence area located northern Shaanxi[J]. Coal Science and Technology,2022,50(2):289−299.
[5] 榆林市政府网. 2022年榆林经济社会综述[EB/OL]. (2023-02-16). http://www.yl.gov.cn/xwzx/ylywe/85356.htm. [6] 王双明,杜 麟,宋世杰. 黄河流域陕北煤矿区采动地裂缝对土壤可蚀性的影响[J]. 煤炭学报,2021,46(9):3027−3038. WANG Shuangming,DU Lin,SONG Shijie. Influence of mining ground fissures on soil erodibility in northern Shanxi coal mining area of Yellow River basin[J]. Journal of China Coal Society,2021,46(9):3027−3038.
[7] 宋世杰,孙 涛,杜 麟,等. 陕北煤矿区不同形态的采动地裂缝对土壤可蚀性的影响[J/OL]. 煤炭学报:1−14[2023-11-28]. http://kns.cnki.net/kcms/detail/11.2190.td.20230303.1533.003.html. SONG Shijie,SUN Tao,DU Lin,et al. Effects of mining ground fissures of different shapes on soil erodibility in northern Shaanxi coal mining area influence[J/OL]. Journal of China Coal Society:1−14[2023-11-28]. http://kns.cnki.net/kcms/detail/11.2190.td.20230303.1533.003.html.
[8] 黄 翌,汪云甲,王 猛,等. 黄土高原山地采煤沉陷对土壤侵蚀的影响[J]. 农业工程学报,2014,30(1):228−235. doi: 10.3969/j.issn.1002-6819.2014.01.029 HUANG Yi,WANG Yunjia,WANG Meng,et al. Effect of mining subsidence on soil erosion in mountainous area of the loess plateau[J]. Transactions of the Chinese Society of Agricultural Engineering,2014,30(1):228−235. doi: 10.3969/j.issn.1002-6819.2014.01.029
[9] 尘福艳,郭仲皓,张英海,等. 基于遥感与GIS技术的陕北煤矿区生态环境质量评价:以杨伙盘矿区为例[J]. 中国煤炭,2020,46(6):45−51. doi: 10.3969/j.issn.1006-530X.2020.06.010 CHEN Fuyan,GUO Zhonghao,ZHANG Yinghai,et al. Evaluation of ecological environment quality of northern Shaanxi mining area based on remote sensing and GIS technology:taking Yanghuopan mining area as an example[J]. China Coal,2020,46(6):45−51. doi: 10.3969/j.issn.1006-530X.2020.06.010
[10] KARAN Shivesh Kishore,GHOSH Somaparna,SAMADDER Sukha Ranjan. Identification of spatially distributed hotspots for soil loss and erosion potential in mining areas of Upper Damodar Basin–India[J]. CATENA,2019,182:104144. doi: 10.1016/j.catena.2019.104144
[11] 刘 英,魏嘉莉,岳 辉,等. 神东矿区土壤侵蚀时空特征及驱动力分析[J]. 测绘科学,2022,47(1):142−153. LIU Ying,WEI Jiali,YUE Hui,et al. Analysis on temporal and spatial characteristics and driving factors of soil erosion in Shendong mining area[J]. Science of Surveying and Mapping,2022,47(1):142−153.
[12] 聂小军,高 爽,陈永亮,等. 西北风积沙区采煤扰动下土壤侵蚀与养分演变特征[J]. 农业工程学报,2018,34(2):127−134. doi: 10.11975/j.issn.1002-6819.2018.02.017 NIE Xiaojun,GAO Shuang,CHEN Yongliang,et al. Characteristics of soil erosion and nutrients evolution under coal mining disturbance in aeolian sand area of northwest China[J]. Transactions of the Chinese Society of Agricultural Engineering,2018,34(2):127−134. doi: 10.11975/j.issn.1002-6819.2018.02.017
[13] 李建明,王文龙,王 贞,等. 神府煤田废弃堆积体新增水土流失研究[J]. 自然灾害学报,2014,23(2):239−249. LI Jianming,WANG Wenlong,WANG Zhen,et al. Study on newly increased soil and water loss from waste accumulation in Shenfu coal-field[J]. Journal of Natural Disasters,2014,23(2):239−249.
[14] SONG S,ZHENG B,SUN T,et al. Influence of different mining damage types on soil erodibility incoal mining areas of northern Shaanxi in the middle reaches of the Yellow River in China[J]. Sustainability,2023,15(6):5434. doi: 10.3390/su15065434
[15] 方学敏,万兆惠,徐永年. 土壤抗蚀性研究现状综述[J]. 泥沙研究,1997(2):88−92. FANG XueMing,WAN Zhaohui,XU Yongnian. Review on research status of soil anti-erodibility[J]. Journal of Sediment Research,1997(2):88−92.
[16] 朱显谟. 黄土高原水蚀的主要类型及其有关因素[J]. 水土保持通报,1982(1):25−30. ZHU Xianmo. Main types of water erosion and its related factors in the loess plateau[J]. Bulletin of Soil and Water Conservation,1982(1):25−30.
[17] 高维森,王佑民. 土壤抗蚀抗冲性研究综述[J]. 水土保持通报,1992(5):59−63. GAO Weishen,WANG Youmin. Review on soil anti-erodibility and anti-scourability[J]. Bulletin of Soil and Water Conservation,1992(5):59−63.
[18] AMEZKETA E,SINGER M J,LE BISSONNAIS Y. Testing a new procedure for measuring water-stable aggregation[J]. Soil Science Society of America Journal,1996,60(3):888−894. doi: 10.2136/sssaj1996.03615995006000030030x
[19] WOODBURN R,KOZACHYN J. A study of related erodibility of a group of Mississippi gully soils[J]. Eos Transactions American Geophysical Union,1956,37(6):749−753. doi: 10.1029/TR037i006p00749
[20] 李程程,曾全超,贾培龙,等. 黄土高原土壤团聚体稳定性及抗蚀性能力经度变化特征[J]. 生态学报,2020,40(6):2039−2048. LI Chengcheng,ZENG Quanchao,JIA Peilong,et al. Characteristics of soil aggregate stability and corrosion resistance longitude change in the loess plateau[J]. Acta Ecologica Sinica,2020,40(6):2039−2048.
[21] 毕银丽,刘 京,尚建选,等. 陕北采煤沉陷区土壤水分入渗和蒸发特征研究[J]. 中国矿业大学学报,2022,51(5):839−849. doi: 10.3969/j.issn.1000-1964.2022.5.zgkydxxb202205002 BI Yinli,LIU Jing,SHANG Jianxuan,et al. Study on the characteristics of soil moisture infiltration and evaporation in the coal mining subsidence area of coal mines in northern Shaanxi[J]. Journal of China University of Mining & Technology,2022,51(5):839−849. doi: 10.3969/j.issn.1000-1964.2022.5.zgkydxxb202205002
[22] 宋世杰,张玉玲,王双明,等. 陕北煤矿区采动地裂缝对土壤微生物和酶活性的影响[J]. 煤炭学报,2021,46(5):1630−1640. SONG Shijie,ZHANG Yuling,WANG Shuangming,et al. Influence of mining ground fissures on soil microorganism and enzyme activities in northern Shaanxi coal mining area[J]. Journal of China Coal Society,2021,46(5):1630−1640.
[23] 高维森. 土壤抗蚀性指标及其适用性初步研究[J]. 水土保持学报,1991,5(2):60−65. GAO Weisen. Primary study of classification of different indices of soil anti-erosibility and its adaptation in the hill region of the loess plateau[J]. Journal of Soil and Water Conservation,1991,5(2):60−65.
[24] 王佑民,郭培才,高维森. 黄土高原土壤抗蚀性研究[J]. 水土保持学报,1994,8(4):11−16. WANG Youmin,GUO Peicai,GAO Weishen. A study on soil anti-erodibility in loess plateau[J]. Journal of Soil and Water Conservation,1994,8(4):11−16.
[25] 郭培才,张振中,杨开宝. 黄土区土壤抗蚀性预报及评价方法研究[J]. 水土保持学报,1992(3):48−51,58. doi: 10.3321/j.issn:1009-2242.1992.03.002 GUO Peicai,ZHANG Zhenzhong,YANG Kaibao. A study on soil anti-erodibility prediction and evaluation method in loess region[J]. Journal of Soil and Water Conservation,1992(3):48−51,58. doi: 10.3321/j.issn:1009-2242.1992.03.002
[26] 刘思璇,高建恩,李文证,等. 黄土丘陵区边坡开挖对土壤抗蚀性影响[J]. 水土保持研究,2022,29(2):23−30. LIU Sixuan,GAO Jianen,LI Wenzheng,et al. Effect of slope excavation on anti-erosion of soil in loess hilly region[J]. Research of Soil and Water Conservation,2022,29(2):23−30.
[27] 魏江生,何金军,高 永,等. 黄土丘陵区土壤水分时空变化特征对采煤沉陷的响应[J]. 水土保持通报,2008,28(5):66−69,103. WEI Jiangsheng,HE Jinjun,GAO Yong,et al. Response of spatial and temporal characteristics of soil moisture to mining subsidence in loess hilly region[J]. Bulletin of Soil and Water Conservation,2008,28(5):66−69,103.
[28] 薛 萐,李占斌,李鹏等. 不同植被恢复模式对黄土丘陵区土壤抗蚀性的影响[J]. 农业工程学报,2009,25(S1):69−72. XUE Sha,LI Zhanbin,LI Peng,et al. Effects of different vegetation restoration models on soil anti-erodibility in loess hilly area[J]. Transactions of the Chinese Society of Agricultural Engineering,2009,25(S1):69−72.
[29] 许明祥,刘国彬,赵允格. 黄土丘陵区土壤质量评价指标研究[J]. 应用生态学报,2005,16(10):1843−1848. doi: 10.3321/j.issn:1001-9332.2005.10.009 XU Mingxiang,LIU Guobin,ZHO Yunge. Assessment indicators of soil quality in hilly loess plateau[J]. Chinese Journal of Applied Ecology,2005,16(10):1843−1848. doi: 10.3321/j.issn:1001-9332.2005.10.009
[30] 王 彬. 土壤可蚀性动态变化机制与土壤可蚀性估算模型[D]. 杨凌:西北农林科技大学,2014,46−48. WANG Bin. Dynamic mechanism of soil erodibility and soil erodibility calculation model[D]. Yangling: Northwest A&F University,2014,46−48.
[31] 刘 鹏,蒋忠诚,李衍青,等. 岩溶断陷盆地高原面洼地不同坡位土壤抗蚀性研究[J]. 地球学报,2021,42(3):373−381. LIU Peng,JIANG Zhongcheng,LI Yanqing,et al. A study of soil anti-erodibility of different slope positions on plateau depression in karst gabin basin[J]. Acta Geoscientica Sinica,2021,42(3):373−381.
[32] 栗 丽,王曰鑫,王卫斌. 采煤塌陷对黄土丘陵区坡耕地土壤理化性质的影响[J]. 土壤通报,2010,41(5):1237−1240. LI Li,WANG Yuexin,WANG Weibin. Effects of mining subsidence on physical and chemical properties of soil in slope land in hilly-gully region of loess plateau[J]. Chinese Journal of Soil Science,2010,41(5):1237−1240.
[33] 夏玉成,冀伟珍,孙学阳,等. 渭北煤田井工开采对土壤理化性质的影响[J]. 西安科技大学学报,2010,30(6):677−681. doi: 10.3969/j.issn.1672-9315.2010.06.009 XIA Yucheng,JI Weizhen,SUN Xueyang,et al. Effect of underground mining on soil physical and chemical properties in Weibei coalfield[J]. Journal of Xi’an University of Science and Technology,2010,30(6):677−681. doi: 10.3969/j.issn.1672-9315.2010.06.009
[34] SHI P,ZHANG Y,HU Z,et al. The response of soil bacterial communities to mining subsidence in the west China aeolian sand area[J]. Applied Soil Ecology,2017,121:1−10. doi: 10.1016/j.apsoil.2017.09.020
[35] 杜华栋,赵晓光,张 勇,等. 榆神府覆沙矿区采煤塌陷地表层土壤理化性质演变[J]. 土壤,2017,49(4):770−775. DU Huadong,ZHAO Xiaoguang,ZHANG Yong,et al. Evolution of topsoil physical-chemical properties after coal mining subsidence in Yushenfu sand covered mining area[J]. Soils,2017,49(4):770−775.
[36] 王双明,杜华栋,王生全. 神木北部采煤塌陷区土壤与植被损害过程及机理分析[J]. 煤炭学报,2017,42(1):17−26. WANG Shuangming,DU Huadong,WANG Shengquan. Analysis of damage process and mechanism for plant community and soil properties at northern Shenmu subsidence mining area[J]. Journal of China Coal Society,2017,42(1):17−26.
[37] 宋世杰,王 艺,彭芮思,等. 陕北不同地貌类型区采煤沉陷对土壤微生物和酶活性的影响[J]. 煤炭科学技术,2023,51(12):110−124. SONG Shijie,WANG Yi,PENG Ruisi,et al. Effects of coal mining subsidence on soil microorganisms and enzyme activities in different landform types of northern Shaanxi[J]. Coal Science and Technology,2023,51(12):104−124.
[38] 李全生,贺安民,曹志国. 神东矿区现代煤炭开采技术下地表生态自修复研究[J]. 煤炭工程,2012(12):120−122. LI Quansheng,HE Anmin,CAO Zhiguo. Study on surface ecological self-repair under modern coal mining technology in Shendong mining area[J]. Coal Engineering,2012(12):120−122.
[39] 胡振琪,龙精华,王新静. 论煤矿区生态环境的自修复、自然修复和人工修复[J]. 煤炭学报,2014,39(8):1751−1757. HU Zhenqi,LONG Jinghua,WANG Xinjing,et al. Self-healing natural restoation and artificial restoration ecological environment for coal mining[J]. Journal of China Coal Society,2014,39(8):1751−1757.
[40] 朱义族,李雅颖,韩继刚,等. 水分条件变化对土壤微生物的影响及其响应机制研究进展[J]. 应用生态学报,2019,30(12):4323−4332. ZHU Yizu,LI Yaying,HAN Jigang,et al. Effects of changes in water status on soil microbes and their response mechanism :a review[J]. Chinese Journal of Applied Ecology,2019,30(12):4323−4332.
[41] 杜 涛,毕银丽,邹 慧,等. 地表裂缝对沙柳根际微生物和酶活性的影响[J]. 煤炭学报,2013,38(12):2221−2226. DU Tao,BI Yinli,ZOU Hui,et al. Effects of surface cracks caused by coal mining on microorganisms and enzyme activities in rhizo-sphere of Salix psammophila[J]. Journal of China Coal Society,2013,38(12):2221−2226.
[42] 王 锐,马守臣,张合兵,等. 干旱区高强度开采地表裂缝对土壤微生物学特性和植物群落的影响[J]. 环境科学研究,2016,29(9):1249−1255. WANG Rui,MA Shouchen,ZHANG Hebing,et al. Effects of surface cracks caused by high intensity coal mining on soil microbial characteristics and plant communities in arid regions[J]. Research of Environmental Sciences,2016,29(9):1249−1255.
[43] 李 娜,韩晓增,尤孟阳,等. 土壤团聚体与微生物相互作用研究[J]. 生态环境学报,2013,22(9):1625−1632. doi: 10.3969/j.issn.1674-5906.2013.09.027 LI Na,HAN Xiaozeng,YOU Mengyang,et al. Research review on soil aggregates and microbes[J]. Ecology and Environmental Sciences,2013,22(9):1625−1632. doi: 10.3969/j.issn.1674-5906.2013.09.027
[44] BLANKINSHIP J C,FONTE S J,SIX J,et al. Plant versus microbial controls on soil aggregate stability in a seasonally dry ecosystem[J]. Geoderma,2016,272:39−50. doi: 10.1016/j.geoderma.2016.03.008
[45] ZHENG W,MORRIS E K,RILLIG M C. Ectomycorrhizal fungi in association with Pinus sylvestris seedlings promote soil aggregation and soil water repellency[J]. Soil Biology and Biochemistry,2014,78:326−331. doi: 10.1016/j.soilbio.2014.07.015
[46] 赵 红,袁培民,吕贻忠,等. 施用有机肥对土壤团聚体稳定性的影响[J]. 土壤,2011,43(2):306−311. ZHAO Hong,YUAN Peimin,LYU Yizhong,et al. Effects of organic manure application on stability of soil aggregates[J]. Soils,2011,43(2):306−311.
[47] 薛菁芳,高艳梅,汪景宽,等. 土壤微生物量碳氮作为土壤肥力指标的探讨[J]. 土壤通报,2007,38(2):247−250. doi: 10.3321/j.issn:0564-3945.2007.02.009 XUE Jingfang,GAO Yanmei,WANG Jingkuan,et al. Microbial biomass carbon and nitrogen as an indicator for evaluation of soil fertility[J]. Chinese Journal of Soil Science,2007,38(2):247−250. doi: 10.3321/j.issn:0564-3945.2007.02.009
[48] 杜 涛,毕银丽,张 姣,等. 地表裂缝对青杨根际环境的影响[J]. 科技导报,2013,31(2):45−49. DU Tao,BI YinLi,ZHANG Jiao,et al. Effects of cracks on the rhizospheric environment of populus cathayana[J]. Science & Technology Review,2013,31(2):45−49.
[49] 史沛丽,张玉秀,胡振琪,等. 采煤塌陷对中国西部风沙区土壤质量的影响机制及修复措施[J]. 中国科学院大学学报,2017,34(3):318−328. doi: 10.7523/j.issn.2095-6134.2017.03.006 SHI Peili,ZHANG Yuxiu,HU Zhenqi,et al. Infulence mechanism of coal mining subsidence on soil quality and restoration measures in west China aeolian sand area[J]. Journal of University of Chinese Academy of Sciences,2017,34(3):318−328. doi: 10.7523/j.issn.2095-6134.2017.03.006
[50] 王文龙,李占斌,张平仓. 神府东胜煤田开发中诱发的环境灾害问题研究[J]. 生态学杂志,2004(1):34−38. doi: 10.3321/j.issn:1000-4890.2004.01.008 WANG Wenlong,LI Zhanbin,ZHANG Pingcang. Environmental disaster issues induced by coal exploitation in Shenfu-Dongsheng coal field[J]. Chinese Journal of Ecology,2004(1):34−38. doi: 10.3321/j.issn:1000-4890.2004.01.008
[51] 王 库. 植物根系对土壤抗侵蚀能力的影响[J]. 土壤与环境,2001(3):250−252. WANG Ku. Effects of plant roots on soil anti-erosion[J]. Soil and Environmental Sciences,2001(3):250−252.
[52] 刘定辉,李 勇. 植物根系提高土壤抗侵蚀性机理研究[J]. 水土保持学报,2003(3):34−37,117. doi: 10.3321/j.issn:1009-2242.2003.03.010 LIU Dinghui,LI Yong. Mechanism of plant roots improving resistance of soil to concentrated flow erosion[J]. Journal of Soil and Water Conservation,2003(3):34−37,117. doi: 10.3321/j.issn:1009-2242.2003.03.010
[53] 刘均阳,周正朝,苏雪萌. 植物根系对土壤团聚体形成作用机制研究回顾[J]. 水土保持学报,2020,34(3):267−273,298. LIU Junyang,ZHOU Zhengchao,SU Xuemeng. Review of the mechanism of root system on the formation of soil aggreates[J]. Journal of Soil and Water Conservation,2020,34(3):267−273,298.
[54] 彭新华,张 斌,赵其国. 土壤有机碳库与土壤结构稳定性关系的研究进展[J]. 土壤学报,2004,41(4):618−623. doi: 10.3321/j.issn:0564-3929.2004.04.019 PENG Xinhua,ZHANG Bin,ZHAO Qiguo. A review on relationship between soil organic carbon pools and soil structure stability[J]. Acta Pedologica Sinica,2004,41(4):618−623. doi: 10.3321/j.issn:0564-3929.2004.04.019
[55] 蒙仲举,任晓萌,陈晓燕,等. 采煤塌陷对沙柳根系损伤机理研究[J]. 北方园艺,2014(1):66−68. MENG Zhongju,REN Xiaomeng,CHEN Xiaoyan. Study on damage mechanism of roots of salix psammophila by mining subsidence[J]. Northern Horticulture,2014(1):66−68.
[56] 张旭冉,张卫青. 土壤团聚体研究进展[J]. 北方园艺,2020(21):131−137. ZHANG Xuran,ZHANG Weiqing. Research progress of soil aggregates[J]. Northern Horticulture,2020(21):131−137.
[57] 刘亚龙,王 萍,汪景宽. 土壤团聚体的形成和稳定机制:研究进展与展望[J]. 土壤学报,2023,60(3):627−643. LIU Yalong,WANG Ping,WANG Jingkuan. Formation and stability mechanism of soil aggregates:progress and prospect[J]. Acta Pedologica Sinica,2023,60(3):627−643.
[58] 郑慧慧,秦佳星,桑之婷,等. 基于区域特征的煤炭开采沉陷对土壤特性影响研究进展[J]. 土壤通报,2022,53(6):1481−1491. ZHENG Huihui,QIN Jiaxing,SANG Zhiting,et al. Progress review on the influence of coal mining subsidence on soil properties based on regional characteristics[J]. Chinese Journal of Soil Science,2022,53(6):1481−1491.
-
期刊类型引用(1)
1. 李小辉，邢楷，彭群. 含水纵向裂隙煤岩破裂特征与演化规律研究. 矿业研究与开发. 2024(07): 191-196 . 
百度学术
其他类型引用(2)
计量
- 文章访问数: 94
- HTML全文浏览量: 14
- PDF下载量: 34
- 被引次数: 3
