马兰煤矿矿井水水质变化特征及成因

李文生, 孙红福

李文生, 孙红福. 马兰煤矿矿井水水质变化特征及成因[J]. 煤田地质与勘探, 2013, 41(4): 46-49. DOI: 10.3969/j.issn.1001-1986.2013.04.012
引用本文: 李文生, 孙红福. 马兰煤矿矿井水水质变化特征及成因[J]. 煤田地质与勘探, 2013, 41(4): 46-49. DOI: 10.3969/j.issn.1001-1986.2013.04.012
LI Wensheng, SUN Hongfu. Variation characteristics and cause of mine water quality in Malan coal mine[J]. COAL GEOLOGY & EXPLORATION, 2013, 41(4): 46-49. DOI: 10.3969/j.issn.1001-1986.2013.04.012
Citation: LI Wensheng, SUN Hongfu. Variation characteristics and cause of mine water quality in Malan coal mine[J]. COAL GEOLOGY & EXPLORATION, 2013, 41(4): 46-49. DOI: 10.3969/j.issn.1001-1986.2013.04.012

 

马兰煤矿矿井水水质变化特征及成因

基金项目: 

国家自然科学基金项目(40972110;41102096);国家国际科技合作项目(2012DFG71060)

详细信息
    作者简介:

    李文生(1968-),男,山西运城人,硕士,高级工程师,从事煤矿防治水技术和管理工作.

  • 中图分类号: P641;X142

Variation characteristics and cause of mine water quality in Malan coal mine

  • 摘要: 通过对山西省马兰煤矿2号煤层采掘面在开采和封闭时期的矿井水和沉积物的研究,揭示采掘面封闭前后对矿井水水质和沉积物的影响机理。研究结果表明:马兰煤矿矿井水均为SO4-Ca型水质,矿井水均富含SO42-和Fe离子;随着上部煤层的不断开采,3处矿井水呈现相同的变化规律,矿井水的pH值升高,Eh值降低,SO42-、Fe、Mn和Zn离子浓度随之下降,其中北一暗斜井处的矿井水水质变化最显著;矿井水水质指标和流速变化能够控制其沉积物的矿物组成和结晶程度,北一暗斜井处的沉积物在两次采样中由斯沃特曼铁矿变为针铁矿,而其他两处的矿井水沉积物矿物组分没有发生变化,主要由针铁矿组成。研究结果能够提高对老空区积水水质的预测精度,并对煤矿突水水源判识具有重要意义。
    Abstract: Quality characteristics of mine water and mineral compositions of sediments at different time from the Malan coal mines, Shanxi Province, China were investigated. The results show that water quality types of all mine water samples were SO4-Ca, and mine water samples had high concentrations of sulfate and iron. With mining of the upper coal seam, pH increased and Eh value, SO42-, Fe, Mn and Zn ions concentrations of mine water samples decreased. Mine water samples in three sampling sites showed the same change patterns and the most significant changes of mine water quality happened in No.1 northern inclined shaft. Mineral composition and crystalline degree of sediments are controlled by water quality index and velocity of mine water. The minerals of sediment in No.1 northern inclined shaft have changed form schwertmannite to goethite in twice sampling. There were no changes of mineral constituents in other two sediment sampling sites, which were composed of goethite.
  • 地下矿产采出后,开采区周围岩土体的原始应力平衡状态被破坏,岩土体上出现位移与变形[1],这是一个十分复杂的过程,是采矿、地质等许多因素及多场耦合综合影响的结果,在地表上表现为塌陷坑与地裂缝[2],属于矿山地质灾害的范畴,具有渐进性特点,本文称为沉陷灾害。沉陷灾害是地下采煤扰动导致地表快速或缓慢沉降的一种地质灾害,在矿区分布广泛,既造成大量土地资源损毁,又严重破坏当地的地形地貌景观,还会严重影响周边植物生长[2],也是诱发斜坡地质灾害和影响区域生态地质环境的重要因素。因此,开展地下采煤区沉陷灾害发育敏感性分区研究对国土综合整治与生态修复、下伏老采空区重大工程项目设计与稳定性监测、地质灾害防治与预警等至关重要,同时对地下采煤区生态环境修复战略[3]具有助力作用。

    灾害敏感性是指在特殊地形或某些因素作用下发生灾害的可能性[4],灾害敏感性评价则是针对某一地区现存或潜在灾害空间分布概率的定性或定量分析。目前,地质灾害敏感性评价工作已经取得众多研究成果[5],针对地表沉陷灾害敏感性评价,主要集中在开采地下水引起的地表沉降[6-8]和废弃矿区地表沉降的灾害敏感性评价研究[9-11],提出了众多评价模型并被应用和开展对比分析[12-14],为沉陷灾害发育敏感性分区研究提供了评价模型选择参考。敏感性评价模型一般可分为定性模型和定量模型,随着研究的不断深入,更多评价模型被引入该领域[15-16],为灾害风险评估与定量评价提供了良好的模型基础。目前,常用的评价模型主要包括:频率比(Frequency Ratio, FR) [17]、确定性系数(Certainty Factor, CF) [18]、支持向量机(Support Vector Machine, SVM) [19]、随机森林(Random Forest, RF) [7-8]和自适应神经模糊推理系统(Adaptive-Network-Based Fuzzy Inference System, ANFIS) [20]等。在模型应用中,确定地表沉陷影响因子是沉陷灾害发育敏感性分区评价的前提,如有研究者认为覆岩岩性、地表坡度、采空区范围、岩体力学性质、地下水等是诱发地表沉陷的主要影响因素[21-22],沉陷灾害敏感性评价则是在上述基础上所做的进一步研究工作[23]。评价因子的可靠性选择,评价模型的适用性对比,为开展地下采煤区沉陷灾害发育敏感性分区研究提供了丰富的理论与实践参考[21-23]

    山西省太原市西山地区地下采煤强度高,沉陷灾害及次生地质灾害分布广泛,给人民生活安全及经济发展带来较大影响。文献[24]开展了崩塌、滑坡、不稳定斜坡和地面塌陷4种地质灾害与发育特征因子之间关联性的空间统计分析,为太原西山地区沉陷灾害发育敏感性评价提供了特征因子选择参考。文献[25]针对斜坡地质灾害敏感性评价中地势起伏度提取最佳尺度开展了研究,为太原西山地区地质灾害敏感性评价时地势起伏度因子可靠提取提供了依据。

    针对地下采煤区沉陷灾害发育敏感性分区,笔者选择在小样本情况下应用广泛的SVM模型为敏感性评价模型,分别以不同年度核查编录的沉陷灾害作为模型建立与预测能力验证数据,开展太原市西山前山区沉陷灾害发育敏感性分区研究,以期为地下采煤区沉陷灾害及次生地质灾害防治提供科学的方法与数据支持。

    太原市西山前山区位于山西省太原市西部山区(图1),北西以尖草坪区、万柏林区和晋源区行政界线为边界,南东以山前坡脚线为界,总面积441.063 km2。区内采煤活动历史悠久,已近百年,煤矿开采主要为地下开采方式,既有国有大矿,包括西铭矿、杜儿坪矿、官地矿、白家庄矿和西峪煤矿,还有分布广泛且数量庞大的乡办、村办、联营小型煤矿,逐渐在空间上形成了大矿与小矿既有邻近、又有重叠的地下采煤区,地表沉陷情况极其复杂,沉陷灾害及次生地质灾害分布广泛,部分地裂缝发育于山顶基岩或山体出露岩层(图2)。

    图  1  太原市西山前山区地理位置
    Figure  1.  Geographical map of Xishan area of Taiyuan City
    图  2  太原市西山前山区沉陷灾害
    Figure  2.  Subsidence disaster in Xishan area of Taiyuan City

    1) 沉陷灾害数据

    收集到研究区2012年、2014年核查编录的地面塌陷与地裂缝灾害共113处,提取地裂缝的中心位置,以点形式与地面塌陷构建沉陷灾害点图层,其空间分布如图1所示。其中,2012年79处,2014年34处,数据来源于山西省太原市规划和自然资源局。

    2) DEM数据

    在地理空间数据云网站下载了研究区的ASTER GDEM V2数据,并重采样为30 m分辨率的DEM数据,作为研究区的数字地貌因子提取基础数据,高程中误差为±22.723 m[25]

    3) 地质数据

    研究区1∶200000数字地质图来源于山西省太原市规划和自然资源局。

    沉陷灾害发育敏感性分区步骤如下:

    (1) 收集研究区实际核查编录的2012年和2014年的沉陷灾害空间分布数据,以2012年发育的沉陷灾害为训练样本建立模型,以2014年发育的沉陷灾害为验证样本评价模型预测精度。

    (2) 从地形地貌和地质因素两个方面出发建立评价因子组合。

    (3) 以评价因子提取结果输入基于机器学习思想的SVM模型开展沉陷灾害发育敏感性评价,以受试者特征曲线(Receiver Operating Characteristic, ROC)下面积(Area Under the ROC Curve, AUC)作为定量评价指标对分区结果进行精度验证。

    (4) 以自然间断点法进行等级划分并制作研究区沉陷灾害发育敏感性分区图。

    总体技术路线如图3所示。

    图  3  总体技术路线
    Figure  3.  Technology roadmap for the study

    针对地下采煤活动对矿区微地貌、土地利用等带来的扰动影响,文献[26-27]选择研究区中西部的杜儿坪煤矿南三盘区为对象,其中,文献[26]采用1∶10000地形图所生产的3期DEM数据分别代表矿区采前、采中和采后3个时期的地形表面,开展了数字地貌时空演变特征分析,结果表明坡度与地势起伏度统计结果随时间推进而变化明显,坡向也有一定的影响。因此,DEM数据的获取时间对数字地貌因子分析结果是否可靠具有重要约束,灾害敏感性预测研究中应尽可能采用与沉陷灾害发育时间相匹配的DEM数据进行数字地貌因子提取。文献[27]则以采前、采中和采后的5期遥感影像为数据源,采用最小距离分类法进行了矿区土地利用信息提取与演变分析,结果表明,土地利用类型转移过程和方向与煤炭开采扰动过程相吻合。因此,遥感影像时相不同,所提取的土地利用信息则不同,以特定年份的遥感影像所提取的土地利用数据和植被覆盖数据与历史编录的沉陷灾害数据进行相关分析,其结果不确定性极大。

    沉陷灾害影响因子众多,在开采沉陷预计中多采用地形地貌、地质因素、采煤方法、煤层厚度、埋深、工作面尺寸等因子[1,28],专家学者基于不同目的建立了不同的影响因子序列。受地下采煤相关数据获取困难、定量表述不足制约,并考虑方法的通用性,本文在现有研究的基础上[24-25],综合考虑地下采煤区数字地貌特征和土地利用时空演变的影响[26-27],从地形地貌与地质因素出发,选取高程、坡度、坡向、地势起伏度、地面曲率、地层岩组、地质构造共7个因子作为沉陷灾害发育敏感性评价因子组合(表1)。

    表  1  研究区沉陷灾害敏感性评价因子
    Table  1.  Assessment factors of subsidence disaster sensitivity in the study area
    影响因素评价因子数据分类孕灾特征备注
    地形地貌高程/m784~1 829影响植被类型、植被覆盖度连续值
    坡度/(°)0~66.20坡度大易引起滑移破坏连续值
    坡向平面;北;东北;东;东南;南;
    西南;西;西北
    影响光照、风的干湿、降雨、植被覆盖和
    土壤潮湿程度等
    分类值
    地势起伏度/m9~343地形起伏大,易引起滑移破坏连续值,9×9网格
    提取[25]
    地面曲率−17.529~17.072定量表征地形表面某点的扭曲变化程度连续值
    地质因素距地质构造距离/m<200;≥200~400;≥400~600;
    ≥600~800;≥800~1 000;
    ≥1 000~1 200;≥1 200
    地质构造附近易形成地裂缝分类值,距构造距离
    地层岩组洪积;坡洪积;保德组;石千峰组;上石盒子组三段;上石盒子组二段;上石盒子组一段;下石盒子组;山西组;太原组;本溪组;峰峰组;上马家沟组覆岩性质对地表破坏形式具有重要影响分类值,组
    下载: 导出CSV 
    | 显示表格

    以ASTER GDEM V2数据在ArcGIS软件中提取高程、坡度、坡向、地势起伏度和地面曲率,基于1∶200000地质图提取地层岩组与地质构造信息,利用ArcGIS软件制作研究区2012年沉陷灾害发育敏感性评价因子图(图4)。

    图4统计,沉陷灾害在各行政区之间的空间分布规模差异性较大,2012年的平均发育密度为0.179处/km2,万柏林区的发育密度最高,达到0.336处/km2,高于研究区平均水平。其次是晋源区,为0.034处/km2,发育密度相对较低。

      4  太原市西山前山区沉陷灾害评价因子
      4.  Assessment factors of the subsidence disaster in Xishan area of Taiyuan City

    支持向量机是一种基于最大间隔的线性判别分类方法,其功能强大,可以最大化模型的预测准确度,被认为是目前针对小样本统计估计和预测学习的最佳理论,最早被用于研究小样本情况下的机器学习[29],在地质灾害敏感性评价研究中应用较多[13,30]

    SVM通过引入核函数有效地解决了非线性分类问题,使敏感性评价中的非线性分类计算复杂度不再取决于空间维数,而是取决于用于建模的沉陷灾害样本数量,尤其是其中支持向量的灾害点数量。显然,核函数是SVM克服维数灾难实现沉陷灾害敏感性评价有效分类的关键。SVM常用核函数(表2)主要有线性核函数(Linear Kernel, LN)、多项式核函数(Polynomial Kernel, PL)、径向基核函数(Radial Basis Function, RBF)和Sigmoid核函数(Sigmoid Kernel, SIG)[13]

    表  2  SVM常用核函数[31]
    Table  2.  Common kernel functions of SVM[31]
    核函数类型核函数解决问题
    线性核函数${k}\left({x}_{i},{x}_{j}\right)={x}_{i}^{{\rm{T}}}{x}_{j}$线性
    多项式核函数${k}\left({x}_{i},{x}_{j}\right)={ {(\gamma x}_{i}^{{\rm{T}}}{x}_{j}+C)}^{d}$偏线性
    径向基核函数${k}\left({x}_{i},{x}_{j}\right)=\mathrm{e}\mathrm{x}\mathrm{p}(-{\gamma }{x}_{i}-{x}_{j}^{2})$偏非线性
    Sigmoid核函数${k}\left({x}_{i},{x}_{j}\right)=\mathrm{t}\mathrm{a}\mathrm{n}\mathrm{h}({(\gamma x}_{i}^{{\rm{T}}}{x}_{j}+C)$非线性
    下载: 导出CSV 
    | 显示表格

    目前,地质灾害敏感性评价中径向基核函数应用较多[30],但研究区内沉陷灾害样本数量较少,应对比4种核函数SVM模型的预测精度,选择最优核函数建立SVM预测模型。

    沉陷灾害敏感性分区SVM模型预测步骤如下:

    (1) 建立研究区沉陷灾害数据集,包括训练样本集和验证样本集。

    (2) 提取沉陷灾害点在评价因子图中的值,作为属性值记录下来,包括:高程、坡度、坡向、地势起伏度、地面曲率、地层岩组和距地质构造的距离。

    (3) 建立4种SVM模型,包括LN-SVM模型、PL -SVM模型、RBF-SVM模型和SIG-SVM模型。

    (4) 采用训练样本对4种SVM模型进行训练,统计评价因子在每一种模型中的权重排序,分析模型预测正确率与错误率,避免出现过拟合和欠拟合。

    (5) 选择预测正确率最高的SVM模型作为沉陷灾害敏感性分区预测模型。

    (6) 采用验证样本评价最优SVM模型的预测精度,对比训练样本的预测精度,分析预测结果的合理性。

    敏感性评价是一个二值分类问题,一般将样本分为正类和负类,现有研究中通常采用真正类率和真负类率构建受试者特征曲线(ROC),利用ROC曲线下的面积(AUC值)作为敏感性分区预测精度的度量[13,31],AUC值的范围为[0,1],值越大越好,表明模型预测能力越强,灾害敏感性分区预测结果精度越高。

    根据AUC值的大小可以划分敏感性分区预测模型的精度等级(表3)[13]

    表  3  敏感性预测模型精度等级[13]
    Table  3.  Accuracy levels of sensitivity prediction model[13]
    AUC0~0.6>0.6~0.7>0.7~0.8>0.8~0.9>0.9~1.0
    精度等级一般非常好极好
    下载: 导出CSV 
    | 显示表格

    对于二值分类问题,假设C1为正类,C2为负类,对应的2×2混淆矩阵见表4。其中,TP表示分类方法准确预测为C1的点的数目;FP表示分类方法预测为C1,但实际属于C2的点的数目;FN表示分类方法预测为C2,但实际属于C1的点的数目;TN表示分类方法准确预测为C2的点的数目。

    表  4  二值分类混淆矩阵
    Table  4.  Confusion matrix of binary classification
    预测类真实类
    正类(C1)负类(C2)
    正类(C1)真正类 (TP)假正类 (FP)
    负类(C2)假负类 (FN)真负类 (TN)
    下载: 导出CSV 
    | 显示表格

    敏感性Sp,也称为真正类率,是C1中所有点被正确预测的比例,计算公式[32]如下:

    $$ {{S}}_{{p}}\text=\frac{\text{TP}}{\text{TP+FN}} $$ (1)

    特异性(Xp),也称为真负类率,计算公式[33]如下:

    $$ {{X}}_{{p}}\text=\frac{\text{TN}}{\text{FP+TN}} $$ (2)

    AUC值计算公式[33]如下:

    $$ \mathrm{A}\mathrm{U}\mathrm{C}={\sum }_{p=1}^{n}({X}_{p+1}-{X}_{p})\times ({S}_{p+1}-{S}_{p}-{S}_{p}/2) $$ (3)

    本文按照图3的技术流程开展研究区沉陷灾害发育敏感性分区预测。

    训练样本和验证样本分别随机生成同等数量的非沉陷灾害样本,建立由79(沉陷灾害)+79(非沉陷灾害)个样本点组成的训练样本数据集,共有训练样本点158个。建立由34(沉陷灾害)+34(非沉陷灾害)个样本点组成的验证样本数据集,共有验证样本点68个。其中,沉陷灾害样本点与非沉陷灾害样本点之间的距离为500 m。

    利用SPSS软件统计4种常用核函数SVM模型评价因子的重要性排序,结果对比如图5所示。

    图  5  研究区不同核函数SVM模型沉陷灾害敏感性评价因子权重对比
    Figure  5.  Comparison of the weight of subsidence disaster sensitivity evaluation factors in different kernel function SVM models in the study area

    图5显示,4种模型中评价因子的权重顺序并不相同,其权重从高到低的排列顺序分别为:

    (1) LN-SVM模型:地层岩组,地质构造,高程,坡向,地面曲率,坡度,地势起伏度。

    (2) PL-SVM模型:地层岩组,坡向,地质构造,坡度,地面曲率,地势起伏度,高程。

    (3) RBF-SVM模型:地层岩组,地质构造,坡向,坡度,地面曲率,高程,地势起伏度。

    (4) SIG-SVM模型:地层岩组,地面曲率,坡度,地势起伏度,坡向,高程,地质构造。

    4种核函数SVM模型预测正确率在SPSS软件中的统计结果见表5

    表  5  研究区不同核函数SVM模型预测正确率
    Table  5.  Prediction accuracy of SVM models with different kernel functions in the study area
    模型预测正确率/%错误率/%
    LN-SVM75.3224.68
    PL-SVM92.417.59
    RBF-SVM72.1527.85
    SIG-SVM62.6637.34
      注:正确率=(样本正确数/样本数)×100%;错误率=1−正确率。
    下载: 导出CSV 
    | 显示表格

    表5可以看出,4种核函数SVM模型的预测正确率由高到低排列为:PL-SVM模型、LN-SVM模型、RBF-SVM模型、SIG-SVM模型,其中,PL-SVM模型表现最优。

    通过因子权重分析及预测正确率对比,PL-SVM模型的因子重要性排序更合理,预测正确率最高,明显优于其他3种模型。因此,本文采用PL-SVM模型进行研究区沉陷灾害发育敏感性分区预测,并在ArcGIS中采用自然间断点法将研究区划分为4个敏感性等级,分别为:极高敏感区、高敏感区、中等敏感区和低敏感区,得到研究区沉陷灾害敏感性分区结果(图6a),验证结果如图6b所示,对训练样本和验证样本的分区统计结果见表6

    图  6  研究区沉陷灾害敏感性分区与验证
    Figure  6.  Sensitivity zoning and verification of subsidence disaster in the study area
    表  6  研究区沉陷灾害敏感性分区统计
    Table  6.  Statistics on sensitivity zoning of subsidence disaster in the study area
    敏感性
    等级
    分区
    栅格数
    栅格比例/
    %
    训练
    样本数
    灾害点比例/
    %
    验证
    样本数
    灾害点比例/
    %
    训练样本
    频率比
    验证样本
    频率比
    极高8 89920.193949.371338.242.451.89
    7 68017.432329.11926.471.671.52
    9 33421.181417.72720.590.840.97
    18 16041.2033.80514.710.090.36
    下载: 导出CSV 
    | 显示表格

    表6可以看出, PL-SVM模型通过训练样本和验证样本分别统计得到的频率比值(频率比=灾害点比例/分区栅格比例)均随敏感性等级升高而递增,经对比其趋势均符合线性函数,拟合方程分别为:

    (1) 训练样本频率比趋势:y = −0.791x + 3.24,R² = 0.999 7。

    (2) 验证样本频率比趋势:y = −0.514x + 2.47,R² = 0.988 7。

    其中,y为频率比值;x为敏感性等级,极高为1,高为2,中为3,低为4。

    可以看出,2个拟合方程的决定系数R²的值均接近于1,表明PL-SVM模型通过训练样本和验证样本所计算的频率比值与敏感性等级之间呈良好的正相关,均以更小的极高与高敏感性分区面积分布了更多的沉陷灾害。训练样本频率比的趋势较验证样本频率比的趋势相关性更好,表明2012年与2014年发育的沉陷灾害在研究区内的空间分布略有不同。

    利用式(3)分别计算训练样本与验证样本在PL-SVM模型敏感性分区预测结果中的AUC值,其中,以2014年沉陷灾害作为验证样本的PL-SVM预测模型的AUC值为0.755,低于以2012年沉陷灾害作为训练样本的PL-SVM预测模型的AUC值(0.854),表明PL-SVM模型在验证样本集的预测精度略低于训练样本集,符合预测模型在测试数据集上的精度一般低于训练数据集精度的规律。但2种模型预测精度均较高,训练模型精度等级为非常好,预测模型精度等级为好。预测精度表明,模型在训练样本与验证样本上既没有过拟合,也没有欠拟合,结果较为合理。

    图5可以看出,地层岩组在4种模型中均是最重要的评价因子,表明沉陷灾害发育与地层岩组分布关系密切,是地面塌陷与地裂缝发育的最重要影响因素。尽管在SIG-SVM模型中地质构造表现为影响最弱的因子,但从4种模型中前3位因子出现的频率可以看出,地质构造对沉陷灾害发育的影响仅次于地层岩组,表明在构造断裂附近,岩层破碎,地下采煤扰动再次打破断裂附近的岩体力学平衡,诱发地表塌陷或开裂。

    图5可知,地势起伏度因子在LN-SVM和RBF-SVM模型中均是沉陷灾害发育影响最弱的因子,在PL-SVM模型中为次最弱因子,表明区域范围内高程最大值与最小值之间的差值大小对沉陷灾害发育影响较小。一般情况下,当沉陷灾害处于沟谷边缘时极易诱发崩塌与滑坡灾害,而沉陷灾害仅包括地面塌陷与地裂缝,因此,地势起伏度因子在敏感性评价中贡献较小。

    图5同时可以看出,LN-SVM和RBF-SVM模型中地层岩组的权重远超其他因子, SIG -SVM模型中各因子的权重相对平衡,地层岩组略高于其他因子,而PL-SVM模型中地层岩组与坡向因子的权重明显高于其他因子,在评价因子组合分析中,既没有过分依赖某一个因子,也没有分散因子权重,其合理性高于其他3种核函数模型。

    SVM预测模型的性能优劣主要取决于核函数的选择,本文分别对比了LN-SVM模型、PL-SVM模型、RBF-SVM模型和SIG-SVM模型在研究区沉陷灾害敏感性分区中的预测精度,预测正确率分别为75.32%、92.41%、72.15%、62.66%,其中,PL-SVM模型预测正确率明显高于其他3种核函数SVM模型,分别比LN-SVM模型、RBF-SVM模型和SIG-SVM模型的预测正确率高22.69%、28.08%和47.48%。SIG-SVM模型的预测正确率是4种模型中最低的,分别比LN-SVM模型、PL-SVM模型和RBF -SVM模型的预测正确率低16.81%、32.19%和13.15%。

    图5可知, PL-SVM模型的评价因子权重排序合理性最好,而SIG-SVM模型则最差,与预测正确率完全一致,选择PL-SVM模型作为研究区沉陷灾害发育敏感性分区预测模型是可靠的。

    敏感性评价研究中通常采用总样本的70%或80%作为训练点建立预测模型,剩余30%或20%作为验证点验证预测模型精度[30-31]。本文以研究区2012年和2014年实际发育的沉陷灾害分别作为训练样本和验证样本,基于PL-SVM模型进行沉陷灾害发育敏感性分区预测。结果表明,2012年沉陷灾害发育在极高、高、中、低4个等级的数量百分比分别为49.37%、29.11%、17.72%、3.80%,极高和高敏感区分布了共78.48%的沉陷灾害。2014年沉陷灾害发育在极高、高、中、低4个等级的数量百分比分别为38.24%、26.47%、20.59%、14.71%,极高和高敏感区分布了共64.71%的沉陷灾害。

    图6可知,2012年沉陷灾害在研究区内有明显的聚集特征,2014年则演变为随机分布,部分沉陷灾害出现在低敏感区,导致敏感性分区预测结果对验证样本(AUC=0.755)的精度略低于训练样本(AUC=0.854),但图6b可以看出,模型的预测结果良好,2014年新发育的沉陷灾害大部分位于极高和高敏感区,为沉陷灾害核查提供了重点区预测,可极大减少野外调查的盲目性与野外工作量。

    目前,地下采煤区地质灾害敏感性评价研究中一般将多年灾害点混合作为样本进行模型训练,但忽视了受采煤扰动影响的矿区地貌特征和土地利用时空演变[26-27],采用特定年份的灾害数据进行建模又需要面对样本数量过少对敏感性评价结果所带来的不确定性影响。今后应就此问题继续展开研究,利用遥感解译手段增加特定年份灾害样本数量,提升区域灾害点小样本情况下沉陷灾害发育敏感性分区预测精度。

    研究区含煤地层出露,极高和高敏感区主要分布于王封煤业、西铭煤矿、杜儿坪煤矿、西峪煤矿、白家庄煤矿等多个大型煤矿区域,相比较地形地貌与地质因素,大、小矿地下采煤扰动是研究区沉陷灾害发育的主控因素。本文基于方法适用性考虑,未将地下采矿扰动影响列为敏感性评价指标,今后应继续研究如何定量化引入地下采矿相关指标参与沉陷灾害敏感性评价。

    太原市西山地区近年所开展的采煤沉陷区治理效果明显,玉泉山城郊森林公园、长风城郊森林公园、四达沟生态恢复景区等分布于研究区东缘,生态地质环境恢复较好,但沉陷灾害极高和高敏感区与植被覆盖区重叠度较高,对区内地下水资源保护形成潜在威胁,应防止水资源破坏导致植被衰退而引起矿区生态环境再次恶化[28],也是矿山生态环境治理需要重点关注的区域。同时,分布于林地内的沉陷灾害点给遥感解译和野外调查带来困难,灾害数据集完备性受限,给敏感性分区评价结果的精准性与可靠性带来一定的影响。

    a. PL-SVM模型是研究区沉陷灾害敏感性分区预测的最优模型,分别以2012年和2014年实际发育沉陷灾害点为训练模型和验证模型,其训练精度与验证精度均较高,且训练精度较验证精度高,对地下采煤区沉陷灾害发育敏感性分区具有良好的适用性和预测性。

    b. PL-SVM模型划分了研究区沉陷灾害极高、高、中和低敏感区4个等级的面积比例,极高与高敏感区主要分布于万柏林区和晋源区边山一带,与研究区内大矿与小矿的重合带高度吻合。灾害点比例随敏感性等级升高而递增,与敏感性等级之间呈现良好的线性正相关。分区结果合理,对研究区沉陷灾害普查与防治重点区确定具有指导价值,对开展同类研究具有参考意义。

    c. 应进一步改进、完善评价因子组合,并选择更多评价模型开展对比,组建更可靠的沉陷灾害发育敏感性预测体系,同时收集研究区最新的沉陷灾害数据进行模型验证,为地下采煤区沉陷灾害发育敏感性分区预测提供最适用的结果。

  • 期刊类型引用(4)

    1. 夏玉成,孙学阳,苗霖田,郭晨,杜荣军. 智能时代的矿井地质工作展望——矿井开采智能地质保障技术体系架构. 煤田地质与勘探. 2025(01): 64-76 . 本站查看
    2. 杨智强,魏中凯,李锐鹏. 渑池县地质灾害敏感因子分析及易发性评价. 粉煤灰综合利用. 2024(03): 86-90+96 . 百度学术
    3. 陈刚,郝社锋,蒋波,喻永祥,车增光,刘汉湖,杨容浩. 基于机载LiDAR技术植被茂密区小型滑坡识别与评价. 自然资源遥感. 2024(03): 196-205 . 百度学术
    4. 牛晨昊,焦润成,韩建锋,王晟宇,郭学飞,刘畅,韩玉成,宋国峰. 基于信息量-机器学习耦合模型的采煤沉陷区滑坡敏感性评价. 煤田地质与勘探. 2024(09): 140-153 . 本站查看

    其他类型引用(3)

计量
  • 文章访问数:  56
  • HTML全文浏览量:  6
  • PDF下载量:  4
  • 被引次数: 7
出版历程
  • 收稿日期:  2012-12-06
  • 网络出版日期:  2021-10-22

目录

/

返回文章
返回