Multi-wave joint detection of unfavorable geological bodies in coal mining face
-
摘要: 影响岳南煤矿工作面回采速度的主要因素为煤层上方灰岩富水区域和内部地质隐伏构造。通过采用瞬变电磁法、无线电波透视、槽波地震多种物探方法可以对回采工作面进行透明化探测。探测过程中发现：瞬变电磁法二次波场接收到的感应电压对富水低阻区域较为敏感；观测无线电波透视发射、接收的能量衰减可识别煤层变化区域横向情况；槽波地震的包络振幅、频谱分析与频散曲线综合分析成像，对判别断层走向有良好效果。通过多波场中感应因子对异常区域综合响应特征的识别，并依据多源波场圈定的异常区域进行平面图叠加和综合显现，可有效识别工作面顶板富水区域与煤层内部构造区域的关联情况，清晰判断导水通道情况。试验结果表明，多种物探方法可探明回采工作面内部异常区域并提前预测预报，对今后类似地质条件的探测提供了借鉴依据。Abstract: The main factors restricting the mining speed of the working face in Yuenan Coal Mine are the water-rich area of limestone above the coal seam and the internal geological concealed structure. The transparent detection of the working face is carried out by using transient electromagnetic method, radio wave perspective and slot wave seismic geophysical methods. In the process of detection, it is found that the induced voltage received by the secondary wave field of transient electromagnetic method is more sensitive to the area of rich water and low resistance, and the lateral situation of coal seam can be identified by observing the energy attenuation emitted and received by radio wave perspective. The integrated analysis of envelope amplitude, spectrum and dispersion curve of in-seam wave seismic imaging has good effect on identifying fault strike. Through the identification of the comprehensive response characteristics of the abnormal region by the induction factor in the multi-wave field, the plane map is superimposed according to the abnormal area delineated by the multi-source wave field, and a variety of abnormal characteristics are shown synthetically, which can effectively identify the correlation between the water-rich area of the roof of the working face and the internal structural area of the coal seam, and clearly judge the existence of the water channel. The test results show that the geophysical method can find out the internal abnormal area of the mining face and predict in advance, and the combined detection of various geophysical methods can provide a reference basis for similar geological conditions in the future.
-
0. 引言
2023年12月18日23时59分，甘肃省临夏州积石山县柳沟乡发生6.2级地震，震源深度为 10 km[1]。地震造成甘肃117人遇难，781人受伤，青海34人遇难，198人受伤。此次地震属于逆冲型，最大烈度为Ⅷ度[2]，地震发生后造成大量山体滑坡失稳，震后多次余震更增加了滑坡发生的风险，给当地居民的影响深远。
地震发生往往会诱发大量的崩塌、滑坡、泥石流等地质灾害，严重威胁当地人民的生命以及财产安全[3]。由于我国位于欧亚板块东南部，受印度板块和太平洋板块的挤压碰撞与俯冲作用的影响，地壳变形强烈，是全球陆内地震最为活跃的地区[4]。特别是青藏高原东缘，地质构造复杂，地形地貌多样，断裂带密集分布，地质灾害频发，造成当地众多人员伤亡和大量财产损失[5 − 6]，例如：1933年四川茂县叠溪地震[7]、2008年5月12日汶川地震[8 − 9]、2010年4月14日玉树地震[10 − 11]；2014年8月13日云南鲁甸地震[12]、2017年8月8日四川九寨沟地震[13]以及2022年泸定地震[14]等。受活动构造带强震影响，上述地震导致的地面运动强度超过国家设计规范最高标准（9度罕遇）的水平，导致了严重的基础设施破坏、滑坡灾害和人员伤亡。甘肃积石山地震也造成较大的损失[15]，同时该区域历史上地质灾害频发[16 − 17]，因此，进行积石山震后地质灾害危险性评价，对防灾减灾和保证人民生命财产安全具有现实意义。近年来，国内外学者对地质灾害危险性评价进行了大量的研究，其中，利用GIS技术和统计学方法相融合已成为最有效、最主要的评价方法。这些方法常见的有：层次分析法[18]、信息量法[19]、逻辑回归模型[20]、灰色模型[21]等。由于单一方法的评价结果并不高以及缺乏足够的说服力[18]，因此，一些学者也采用多种方法结合应用于地质灾害风险评价[22]，取得较好的评价效果。可见，多种评价方法交叉融合可综合多方法优势、科学准确地表达灾害风险评价结果。
然而，上述研究多致力于影响因子与滑坡等灾害的特征分析和风险评价，在影响因子对滑坡的响应分析，特别是影响因子对灾害发生风险的贡献率和重要性分析等方面鲜有研究。为此，本文以甘肃积石山Ms6.2级地震区为研究区，基于多源遥感数据，在深入探讨滑坡灾害与影响因子间的特征分布，以及影响因子对滑坡的响应分析基础上，快速进行震区地质灾害危险性评价和分析，以期为灾后救援和重建工作提供参考。
1. 研究区概况
甘肃积石山Ms6.2级地震位于甘肃和青海两省交界地带，地震中心位置为北纬35.70°、东经102.79°，震源深度为 10 km，最高震中烈度达Ⅷ度。此次地震涉及甘肃省3个市（州）9个县（市、区），涉及青海省2个市（州）4个县（市），强烈的地震造成大量人员伤亡、房屋、道路等基础设施损毁（图1）。震区地势上属青藏高原一级阶梯向黄土高原二级阶梯的过渡带，高差起伏较大，最大高差为
3000 余米，区内断裂带活动强烈，有拉脊山南、北断裂带，倒淌河-临夏断裂带，西秦岭北缘断裂带等，为地震等自然灾害的频发区[23]。其中拉脊山断裂带是穿过震区的主要断裂，受印度板块俯冲挤压的动力驱动，青藏高原东北缘和东缘持续扩张，导致该区域新构造活跃、地震频发，为我国同震地质灾害发育最频繁的地区[24]。本次同震地质灾害主要分布于拉脊山断裂两侧，就地形地貌来看，海拔西高东低，西侧主要为基岩高山区，该区域主要分布有古生界砂岩、粉砂岩和泥岩等；右侧区域主要为黄土低山丘陵，由于黄河的下切和强烈的构造抬升，形成广泛分布的黄河阶地和大量高陡边坡存在，为地震同震灾害发育创造了条件。震区的气候类型属于典型的大陆季风性气候，年均降水量650 mm/a，年均蒸发量
1110 mm/a[24]。2. 数据与方法
2.1 数据
（1）卫星影像数据
本次所用光学遥感影像为高分一号和吉林一号，其空间分辨率分别为2 m和0.75 m，震前影像为12月18日，震后影像分别为12月20日和12月19日。影像预处理主要包括辐射定标、大气校正、几何纠正和图像融合等处理，最后得到真彩色高分辨率影像，主要用于滑坡灾害点的目视解译[25]。
（2）地震点数据
本次所用地震和余震点发生的地理位置数据来自国家地震科学数据中心（https://search.asf.alaska.edu/#/）。该数据为截至2023年12月20日、震级3.0以上的资料（表1）。
表 1 甘肃积石山Ms6.2级地震及余震数据Table 1. Ms6.2 Jishishan earthquake and its aftershocks in Gansu Province序号 发震时刻 纬度 经度 深度/km 震级/Ms 地名 备注 1 2023-12-18T23:59:30.0 35°42′ 102°47′ 10 6.2 积石山县 震中 2 2023-12-19T00:24:49.9 35°44′ 102°47′ 10 3.9 积石山县 余震 3 2023-12-19T00:36:18.3 35°47′ 102°47′ 10 4.0 积石山县 余震 4 2023-12-19T00:43:12.9 35°47′ 102°46′ 10 3.4 积石山县 余震 5 2023-12-19T00:56:51.3 35°42′ 102°47′ 10 3.4 积石山县 余震 6 2023-12-19T00:59:11.3 35°44′ 102°46′ 10 3.1 积石山县 余震 7 2023-12-19T00:59:39.0 35°50′ 102°47′ 10 4.1 积石山县 余震 8 2023-12-19T01:10:31.4 35°48′ 102°47′ 10 3.2 积石山县 余震 9 2023-12-19T01:20:12.6 35°48′ 102°46′ 10 3.2 积石山县 余震 10 2023-12-19T02:10:06.4 35°50′ 102°46′ 10 3.2 积石山县 余震 11 2023-12-19T00:32:52.9 35°46′ 102°47′ 9 3.4 积石山县 余震 （3）环境变量数据
本文所使用的环境变量数据，主要包括：地形、断裂带、道路、土地利用、植被、人口、土壤质地，并通过数据处理得到14个影响因子（表2）。其中：地形数据为 ALOS 的数字高程模型（DEM），空间分辨率为 12.5 m，主要用于获取坡度、坡向、剖面曲率、平面曲率、曲率、距河流的距离和地形湿度指数（TWI）；断裂带，下载自国家地震科学数据中心，经处理得到距断裂带距离；道路数据，来自OSG（Open Street Map）官网，用于计算距道路的距离；土地利用为武汉大学的 CLCD（China land cover dataset）数据集，其空间分辨率为 30m，主要包括：耕地、林地、灌木、草地、水域、雪/冰、裸地、不透水面和湿地；归一化植被指数来自国家青藏高原科学数据中心，其数据源为 MODIS，空间分辨率为 250 m；6）人口分布，来自WorldPop 全球人口数据，空间分辨率为100 m，下载自GEE（Google Earth Engine）平台；土壤质地数据，来自世界土壤数据库，其空间分辨率为1 km。
表 2 甘肃积石山Ms6.2级地震诱发滑坡灾害影响因子数据Table 2. Disaster-causing factors of landslides induced by the Ms6.2 Jishishan earthquake in Gansu Province环境变量 影响因子 数据来源 地形因子 高程、坡度、坡向、剖面曲率、平面曲率、
曲率、距河流距离、TWI高程数据为数字高程模型（DEM），下载自ASF官网（https://search.asf.alaska.edu/#/），
其它因子数据为DEM的派生数据断裂带 距断裂带的距离 断裂带数据下载自国家地震科学数据中心（https://search.asf.alaska.edu/#/） 土壤质地 土壤质地数据 下载自世界土壤数据库（https://www.fao.org/） 道路 距道路的距离 OSM官网（https://www.openstreetmap.org/） 人口 人口分布密度 OSM官网（https://www.openstreetmap.org/） 土地利用 地表覆盖 武汉大学CLCD数据集（https://zenodo.org/） 植被 归一化植被指数 国家青藏高原科学数据中心（https://data.tpda.ac.cn/home ） 2.2 研究方法
2.2.1 MaxEnt模型原理
最大熵（MaxEnt）模型是基于最大熵原理来预测随机事件概率分布的方法，该模型在灾害危险性评价、预测物种分布等领域都得到了广泛应用[26 − 27]。利用MaxEnt模型进行地质灾害预测的原理，主要以地质灾害发生点和环境变量之间的关系为基础，建立地质灾害危险性的概率模型，从而对整个区域的地质灾害危险性进行预测。本文是以14个影响因子为约束条件，以灾害点为事件，寻求地质灾害危险性在14个约束条件下的最大熵，进而估算地质灾害的危险性。
2.2.2 数据预处理
对目视解译的
1205 个滑坡灾害点数据进行预处理，包括：1）剔除滑坡目视解译得到的重复和错误点；2）通过建立渔网去除自相关点；3）数据转换，主要是将最终得到的980个滑坡灾害点数据转成csv格式用于模型构建。对影响因子的预处理，包括将致灾因子重采样为统一像元尺度、裁剪使所有影响因子有相同的行列号，然后将其转换为ASCII格式用于构建MaxEnt模型[28]。
2.2.3 模型参数设置
模型参数调整是构建MaxEnt模型的重要环节。具体流程包括，输出格式设置为logistic，将滑坡灾害点分为测试集和训练集，其中，随机测试集设置比例为25%，75%用于模型训练[29]。精度评定采用受试者工作特征曲线（receiver operating characteristic，ROC），将正则化参数（regularization multiplier）设置为1，重复建模次数设置为10，以防止欠拟合和过拟合现象的发生，并且使模型构建达到最佳效果[30]。
2.2.4 危险评价
对滑坡危险性评价结果等级的划分，参考前人研究成果[31 − 32]，并结合高分影像目视解译得到的滑坡灾害点空间分布，将危险区等级划分为五个等级：极低危险区，低危险区，中危险区，高危险区和极高危险区。
2.2.5 精度评价
模型精度评价采用ROC曲线与坐标轴围成的曲线下面积（area under curve，AUC）值来表示，AUC值大小表示预测成功率，其值越大，则准确率越高，说明模型的预测效果就越好[32]。AUC值的具体评价标准如表3 。
表 3 AUC评价标准Table 3. AUC evaluation criteriaAUC值 精度评价 [0, 0.6） 很差 [0.6, 0.7） 较差 [0.7, 0.8） 一般 [0.8, 0.9） 好 [0.9, 1] 极好 3. 结果
3.1 滑坡点的影响因子分布特征
本文通过对比积石山Ms6.2级地震前后的高分一号和吉林一号卫星影像，对滑坡进行目视解译。由于地震前，该区域发生大面积降雪，在地震发生后，发生滑坡的区域裸露出新的土壤，因此有利于滑坡灾害点的卫星影像目视解译。此次目视解译原则为：首选区域内空间分辨率高的影像，若存在云层遮挡或地形阴影，将按照空间分辨率从高到低的顺序，选择时相相近的影像，最终达到覆盖整个震区。在滑坡点目视解译中，主要采取震前震后影像对比的方式，本次解译共编目滑坡和潜在灾害点
1205 处（图2），大多以小型崩塌和滑坡为主，主要集中于震区黄河两岸黄土丘陵区、道路与沟谷附近，多发育于建房和公路内边坡陡坎部位[30]，其威胁的主要承灾体为公路和农田等。![]()
图 2 积石山地震诱发滑坡灾害点空间分布Figure 2. Spatial distribution of landslides induced by the Jishishan earthquake为进一步分析灾害点在各影响因子上的分布特征，研究中将灾害点与各因子进行叠加，分别作直方图进行统计分析，其中横坐标为各因子分级，纵坐标为密度。
对于各地形因子（图3），分析中将高程因子按100 m间隔分为九级，叠加统计分析发现，地震诱发滑坡灾害点在高程因子上基本呈抛物线分布（R2 =
0.7394 ），主要分布在1700 ～2250 m高程带；就灾害点在坡度上的分布，将坡度按5°间隔等间距分级，统计发现灾害点主要发生在小于25°区间内，大于30°坡度范围零星分布；在地形湿度指数（TWI）上，灾害点发生基本呈现指数分布（R2 =0.6816 ）；就坡向分布来看，地震诱发地质灾害点大多数发生于东、东南和南向范围内。对于各类距离因子（图4），研究中将距道路距离因子按间隔0.5 km等间距作缓冲区分级，分析发现地震诱发的绝大多数灾害点发生于距道路1.5 km范围内；对于距断裂带距离，第一级按1 km、第二级以后按2 km进行分级统计，分析发现地震诱发地质灾害点在距断裂带距离上基本呈现指数型分布（R2 = 0.201），主要分布于小于1 km、3～4 km和8～12 km范围内，其中3～4 km和8～12 km区间为黄河两岸；就距河流距离来看，地震诱发地质灾害点分布呈指数型分布（R2 =
0.7727 ），主要分布在河流两岸及附近；对于距地震中心距离因子，研究中第一级按1 km统计，第二级按5 km统计，后面各级按10 km统计，分析发现灾害点呈指数型分布（R2 =0.8706 ）。同时，也对土地利用、土壤质地、归一化植被指数和人口分布因子进行地质灾害点分布统计（图5）。就土地利用来看，地震诱发地质灾害点主要发生于耕地、部分草地和水域附近也有一些地质灾害发生；就土壤质地来看，地震诱发地质灾害点主要分布在壤土、有少量分布于黏土和壤砂土层；就NDVI来看，地震诱发地质灾害主要分布在0.08～0.16；就震区人口分布来看，人口密度大的地方人类活动也较强烈，因而地震诱发地质灾害点的分布也相应多一些。
![]()
图 5 滑坡灾害点的其他因子分布Figure 5. Distribution of landslides based on LULC, soil texture, population density, and NDVI3.2 影响因子对滑坡的响应分析
1）模型评价精度
将滑坡灾害点数据和选取的14个影响因子数据输入MaxEnt模型，通过10次的迭代计算，最后得到AUC值为0.854（图6），模型可靠性达到“好”的水平。因此，本次研究利用解译的地质灾害点和各影响因子，通过10次迭代计算构建MaxEnt模型进行积石山Ms6.2级地震诱发滑坡危险性评价，其结果具有较好的可靠性。
2）影响因子重要性分析
置换重要性是反映模型对该变量的依赖程度的指标[33]。表4为各影响因子对滑坡灾害影响程度的贡献率和置换重要性。可见，排名前五的影响因子为距断裂带的距离、高程、人口分布、土壤质地和距河流的距离，其贡献率分别为39.0%、38.1%、17.8%、1.3%和1.2%，其累计贡献率占比高达97.4%。同时可以看出，置换重要性排名前五的影响因子为距断裂带的距离、高程、距河流的距离、人口分布以及土壤质地，其置换重要性分别为48.3%、45.1%、2.4%、1.4%和1.3%，累计值达98.5%。
表 4 滑坡灾害主要影响因子贡献率和置换重要性Table 4. Contribution rates and permutation importance of main disaster-causing factors of landslides序号 因子 贡献率/% 置换重要性/% 1 距断裂带距离 39 48.3 2 高程 38.1 45.1 3 人口分布 17.8 1.4 4 土壤质地 1.3 1.3 5 距河流距离 1.2 2.4 6 归一化植被指数 0.8 0.6 7 坡度 0.8 0.1 8 距道路距离 0.6 0.5 9 坡向 0.3 0.2 10 地形湿度指数 0.1 0 11 土地利用 0 0 12 平面曲率 0 0 13 剖面曲率 0 0 14 曲率 0 0 图7为通过刀切法检验对各影响因子重要性的检验结果。由测试增益值[34]，见图7（a），重要性位列前五的影响因子为距断裂带的距离、高程、人口分布、土壤质地和归一化植被指数，其值分别为：0.35、0.32、0.22、0.12和0.1。由AUC值，见图7（b），排位前五的影响因子为高程、距断裂带的距离、人口分布、归一化植被指数和距河流的距离，其值分别为0.74、0.73、0.67、0.65和0.59。由正则化训练增益，见图7（c），位列前五的影响因子为高程、距断裂带距离、人口分布、距道路距离和距河流距离，其值分别为0.3、0.29、0.22、0.1 和 0.08。
3）影响因子对滑坡危险性响应分析
图8和图9分别为各影响因子对滑坡发生的响应曲线，其中，纵轴代表滑坡发生的概率，横轴代表各因子的取值范围。设定参考概率阈值为0.5[35]，当大于0.5时，认为该因子取值范围有利于灾害的发生。图8可知，坡向因子对滑坡发生的响应最高。其他因子的某一取值范围对滑坡发生的响应也较敏感。例如，当地形湿度指数大于4 m时，其概率均大于0.5，极易引起滑坡的发生；同理，当高程带在
1700 ～2250 m、剖面曲率为−4.2～3、平面曲率为−3.9～4.1、综合曲率为−6～11时，该范围对滑坡发生强响应；当距断裂带距离小于1.7 km、距河流距离小于3.8 km、距道路距离小于2 km、坡度小于30°、人口分布密度大于20人/km2时，其概率均大于0.5，极易引起滑坡的发生。同时也可以看出，当归一化植被指数小于−0.04和0.06～0.15时，其概率均大于0.5，该区段滑坡灾害发生的响应较好；就土地利用因子来看（图9），耕地、草地和水域的概率均大于0.5，极易发生滑坡灾害；对于土壤质地，砂质壤土和壤土的概率大于0.5，极易发生滑坡灾害。![]()
图 9 土地利用与土壤质地的滑坡响应直方图Figure 9. Histogram of landslide responses to LULC and soil texture3.3 滑坡危险性评价
本文采用影响因子重要性和相关系数法，剔除贡献率较低（即土地利用、平面曲率、泡面曲率和曲率）和相关性高（即高程）的因子，然后将其余因子构建模型，计算最大熵结果，并按自然断点法分五级[36]，图10为本研究获取的积石山Ms6.2级地震诱发滑坡危险性评价结果。经统计得出，极高危险区面积为 49.38 km2，占研究区总面积的0.84%；高危险区面积为157.79 km2，占研究区总面积的 2.69%；中危险区面积为430.03 km2，占研究区总面积的7.33%；低危险区面积为526.07 km2，占研究区总面积的8.96%；极低危险区面积为
4699.02 km2，占研究区总面积的80.18%。可见，由于本次地震发生在冬季，大多数地方为季节性冻土，因而本次地震诱发滑坡大多为小型，其极高和高危险性分布相对较少，主要位于黄河两岸的部分地区，与文献[25]结果一致。![]()
图 10 积石山Ms6.2级地震诱发的滑坡危险性评价结果Figure 10. Landslide hazard assessment results induced by the Ms6.2 Jishishan earthquake为进一步分析危险区与地震烈度之间的关系，将危险性评价结果与地震烈度图[37]进行叠加，统计得出（表5），极高和高危险区密度主要位于地震烈度为Ⅷ区，其面积为21.2 km2，占Ⅷ区的百分比为26.38%；中危险区密度主要分布于Ⅶ和Ⅷ区，面积为341.22 km2，占Ⅶ和Ⅷ区面积百分比分别为16.92%和28.82%；低和极低危险区主要分布于Ⅶ和Ⅵ区，面积占比分别为75.33%和97.55%，该区域远离发震区，地震诱发地质灾害危险性也较低。
表 5 不同地震烈度区的危险性等级面积百分比统计Table 5. Area percentage of different risk grades in various seismic intensity zones地震烈度 极高危险/% 高危险/% 中危险/% 低危险/% 极低危险/% Ⅷ区 6.91 19.47 28.82 11.76 33.05 Ⅶ区 1.80 5.95 16.92 16.12 59.21 Ⅵ区 0.03 0.23 2.19 6.15 91.40 4. 讨论
影响因子的选择是地震诱发滑坡危险性评价可靠性的主要环节[38]。与前人研究不同，本研究从地形因子、断裂带、土壤质地、道路、人口、土地利用以及植被等环境变量出发，选取影响因子，在分析影响因子对滑坡危险性响应的基础上，通过获取贡献率和置换重要性、采用刀切法对影响因子进行重要性评估，来优选主要影响因子。研究发现，距断裂带距离、高程、人口分布、土壤质地和距河流距离的贡献率和置换重要性均大于1，其中，距断裂带距离最大，分别为39.0%和48.3%。
为进一步分析本研究的可靠性，采用刀切法对各影响因子重要性进行检验，由测试增益值、AUC 值以及正则化训练增益三个指标进行计算，得出对本次地震诱发滑坡危险性的主要响应因子为距断裂带距离、高程和人口分布。可见，除属于孕灾环境条件的高程因子外，距断裂带距离和人口分布是本次地震诱发滑坡危险性的主要驱动因子。由于震区相对高差约
3000 m，地形起伏较大，高海拔地区因地形坡度陡峭，土壤和岩石的稳定性相对较差，特别是在1700 ～2300 m区间，人类活动（例如：修路、建房、耕作等）也较强烈，遥感解译发现，大多发生在距离道路较近的地区，道路建设往往会破坏原有地形和土壤结构，为滑坡崩塌发育创造了条件。断裂带区因地壳应力集中以及地质结构脆弱，在地震等外部因素的驱动下，极易引发滑坡灾害；同时，由于河流的侵蚀作用往往会破坏斜坡的稳定性，导致滑坡灾害的发生。实地调查发现，由于本次地震主要发生在冬季，其诱发地质灾害主要以小型崩塌和滑坡为主。从灾害分布区域、规模和密集程度来看，与本研究计算确定的主要影响因素吻合性较好。例如：图11为实地地质灾害调查中，获得的典型性滑坡，其中，图11a和图11b为砂/泥岩滑坡类，大多分布于公路等的内边坡区；图11c和图11d为黄土滑坡类，大多发生于建房和公路等的切坡区；图11e为裂缝，发育于震区黄河两岸阶地居民区。
![]()
图 11 积石山震区典型滑坡注：a—b为砂/泥岩滑坡；c—d为黄土滑坡；e为裂缝。Figure 11. Typical landslides in the Ms6.2 Jishishan earthquake5. 结论
（1）从地形因子来看，本次地震诱发滑坡灾害点主要分布在
1700 ～2250 m高程带、20°～25°的坡度区间，且大多数群发于坡向为东、东南和南向的阳坡范围内，在TWI上基本呈现指数分布，同时在距道路距离1.5 km、距地震中心距离5 km范围内广布。从LULC来看，地震诱发地质灾害点主要发生于耕地，且土壤质地多为壤土区域；在植被覆盖上，主要集中发生于NDVI为0.2～0.4区间、震区人口分布密度大的地方。（2）由影响因子的贡献率和置换重要性、刀切法计算的测试效益值、AUC 值和正则化训练增益值，得到本次地震诱发地灾危险性的主要影响因子为：距断裂带的距离、高程和人口分布。当距断裂带距离小于1.7 km、人口分布密度达20人/km2时，其概率均大于0.5，对滑坡危险性发生的响应较明显。
（3）基于构建的MaxEnt模型得出，震区滑坡高危险区主要分布于黄河两岸及附近区域。其中，极高和高危险区密度主要位于地震烈度为Ⅷ区，其面积为21.2 km2，占Ⅷ区的面积百分比为26.38%；中危险区密度主要分布于Ⅶ和Ⅷ区，面积百分比分别为16.92%和28.82%；低和极低危险区主要分布于Ⅶ和Ⅵ区，面积百分比分别为75.33%和97.55%。
-
![]()
图 3 顺槽巷顶板30°多测道剖面电压曲线图
Figure 3. Voltage curve of multi-channel section at 30° in the roof of the roadway
![]()
图 11 槽波透射能量衰减分布图
Figure 11. Attenuation distribution of in-seam wave transmitted energy
![]()
图 12 瞬变电磁、无线电波叠加示意图
Figure 12. Schematic diagram of superposition of transient electromagnetic and radio waves
表 1 物探异常推断分析表
Table 1 Geophysical anomaly inference analysis
物探方法 异常区域/m 异常推断 导水程度 回采影响 瞬变电磁 0～100 富水区域 弱 一般 500～650 富水区域 强 重大 650～800 富水区域 弱 一般 无线电波 0～80 煤层变化 弱 较大 350～650 地质构造 强 重大 槽波地震 300 地质构造 弱 较大 420～650 地质构造 强 重大 700 地质构造 弱 较大 
下载: 导出CSV
-
[1] 李江华,廉玉广,焦阳,等. 综合物探技术在工作面导水构造探测中的应用[J]. 煤矿安全,2018,49(3):129 − 132. [LI Jianghua,LIAN Yuguang,JIAO Yang,et al. Application of comprehensive geophysical exploration technique in exploration of water conductive structures in working face[J]. Safety in Coal Mines,2018,49(3):129 − 132. (in Chinese with English abstract) LI Jianghua, LIAN Yuguang, JIAO Yang, et al. Application of comprehensive geophysical exploration technique in exploration of water conductive structures in working face[J]. Safety in Coal Mines, 2018, 49（3）: 129-132. （in Chinese with English abstract）
[2] 张平松,胡泽安,吴荣新,等. 煤层工作面地质构造及异常透射CT综合成像方法与应用[J]. 中国煤炭地质,2017,29(9):49 − 52. [ZHANG Pingsong,HU Zean,WU Rongxin,et al. Coal face geological structure and application of anomalous body transmitting CT integrated imaging[J]. Coal Geology of China,2017,29(9):49 − 52. (in Chinese with English abstract) DOI: 10.3969/j.issn.1674-1803.2017.09.10 ZHANG Pingsong, HU Zean, WU Rongxin, et al. Coal face geological structure and application of anomalous body transmitting CT integrated imaging[J]. Coal Geology of China, 2017, 29（9）: 49-52. （in Chinese with English abstract） DOI: 10.3969/j.issn.1674-1803.2017.09.10
[3] 马志超,杨高峰,王克南. 利用透射槽波衰减系数探查煤层内部的断层响应特征研究[J]. 煤炭技术,2021,40(3):49 − 51. [MA Zhichao,YANG Gaofeng,WANG Kenan. Fault response characteristics of coal seam are investigated by using attenuation coefficient of transmitted in-seam wave[J]. Coal Technology,2021,40(3):49 − 51. (in Chinese with English abstract) MA Zhichao, YANG Gaofeng, WANG Kenan. Fault response characteristics of coal seam are investigated by using attenuation coefficient of transmitted In-seam wave[J]. Coal Technology, 2021, 40（3）: 49-51. （in Chinese with English abstract）
[4] 辛光明,邢文彬,武凯,等. 阳城煤矿断层导水灾害“挡-堵”多体系防治技术[J]. 煤矿安全,2019,50(5):91 − 94. [XIN Guangming,XING Wenbin,WU Kai,et al. Multi-system prevention and control technology of “blocking-plugging” in fault water guiding disaster of Yangcheng coal mine[J]. Safety in Coal Mines,2019,50(5):91 − 94. (in Chinese with English abstract) XIN Guangming, XING Wenbin, WU Kai, et al. Multi-system prevention and control technology of “blocking-plugging” in fault water guiding disaster of Yangcheng coal mine[J]. Safety in Coal Mines, 2019, 50（5）: 91-94. （in Chinese with English abstract）
[5] 牟义,李江华,徐慧,等. 矿井瞬变电磁法参数优化试验及超前探测应用[J]. 煤炭科学技术,2020,48(6):184 − 190. [MU Yi,LI Jianghua,XU Hui,et al. Parameters optimization test of mine transient electromagnetic method and application of advanced detection[J]. Coal Science and Technology,2020,48(6):184 − 190. (in Chinese with English abstract) MU Yi, LI Jianghua, XU Hui, et al. Parameters optimization test of mine transient electromagnetic method and application of advanced detection[J]. Coal Science and Technology, 2020, 48（6）: 184-190. （in Chinese with English abstract）
[6] 李宏杰,黎灵,李健,等. 采动覆岩导水断裂带发育高度研究方法探讨[J]. 金属矿山,2015(4):1 − 6. [LI Hongjie,LI Ling,LI Jian,et al. Discussion on the methods for determining the height of fractured water-conducting zone[J]. Metal Mine,2015(4):1 − 6. (in Chinese with English abstract) LI Hongjie, LI Ling, LI Jian, et al. Discussion on the methods for determining the height of fractured water-conducting zone[J]. Metal Mine, 2015（4）: 1-6. （in Chinese with English abstract）
[7] 李宏杰,马君,姜鹏,等. 准格尔煤田岩溶水文地质特征及水害防治技术[J]. 煤矿安全,2018,49(9):246 − 251. [LI Hongjie,MA Jun,JIANG Peng,et al. Karst hydrogeological characteristics and prevention and control technology for water damage in Zhungar coalfield[J]. Safety in Coal Mines,2018,49(9):246 − 251. (in Chinese with English abstract) LI Hongjie, MA Jun, JIANG Peng, et al. Karst hydrogeological characteristics and prevention and control technology for water damage in zhungar coalfield[J]. Safety in Coal Mines, 2018, 49（9）: 246-251. （in Chinese with English abstract）
[8] 牟义,邱浩,牛超,等. 多源干扰条件下瞬变电磁法电性响应规律研究[J]. 地球物理学进展,2019,34(6):2493 − 2502. [MU Yi,QIU Hao,NIU Chao,et al. Study on the electrical response law of transient electromagnetic method under multi-source interference conditions[J]. Progress in Geophysics,2019,34(6):2493 − 2502. (in Chinese with English abstract) DOI: 10.6038/pg2019CC0434 MU Yi, QIU Hao, NIU Chao, et al. Study on the electrical response law of transient electromagnetic method under multi-source interference conditions[J]. Progress in Geophysics, 2019, 34（6）: 2493-2502. （in Chinese with English abstract） DOI: 10.6038/pg2019CC0434
[9] 牟义,徐慧,马志超,等. 带压综采工作面矿井综合物探及安全性评价研究[J]. 重庆大学学报,2019,42(12):41 − 49. [MU Yi,XU Hui,MA Zhichao,et al. Comprehensive geophysical exploration and safety evaluation of the mine with a fully mechanized pressurized mining face[J]. Journal of Chongqing University,2019,42(12):41 − 49. (in Chinese with English abstract) DOI: 10.11835/j.issn.1000-582X.2019.12.005 MU Yi, XU Hui, MA Zhichao, et al. Comprehensive geophysical exploration and safety evaluation of the mine with a fully mechanized pressurized mining face[J]. Journal of Chongqing University, 2019, 42（12）: 41-49. （in Chinese with English abstract） DOI: 10.11835/j.issn.1000-582X.2019.12.005
[10] 李文,牟义,邱浩. 煤矿含水异常体矿井综合物探方法及应用[J]. 煤矿安全,2017,48(7):208 − 211. [LI Wen,MU Yi,QIU Hao. Application of mine comprehensive geophysical detection methods on water bearing abnormal bodies[J]. Safety in Coal Mines,2017,48(7):208 − 211. (in Chinese with English abstract) LI Wen, MU Yi, QIU Hao. Application of mine comprehensive geophysical detection methods on water bearing abnormal bodies[J]. Safety in Coal Mines, 2017, 48（7）: 208-211. （in Chinese with English abstract）
[11] 牟义. 切片技术在煤矿回采工作面瞬变电磁法探测中的研究与应用[J]. 中国煤炭,2014,40(3):35 − 39. [MU Yi. Research and application of slicing technology in transient electromagnetic detection in coal mining face[J]. China Coal,2014,40(3):35 − 39. (in Chinese with English abstract) DOI: 10.19880/j.cnki.ccm.2014.03.008 MU Yi. Research and application of slicing technology in transient electromagnetic detection in coal mining face[J]. China Coal, 2014, 40（3）: 35-39. （in Chinese with English abstract） DOI: 10.19880/j.cnki.ccm.2014.03.008
[12] 马志超,廉玉广,牟义,等. 煤矿井下探测中影响不含水断层视电阻率变化的单一主控因素分析[J]. 能源与环保,2020,42(12):106 − 111. [MA Zhichao,LIAN Yuguang,MU Yi,et al. Analysis of single main controlling factor affecting change of apparent resistivity of water-free faults in underground coal mine exploration[J]. China Energy and Environmental Protection,2020,42(12):106 − 111. (in Chinese with English abstract) MA Zhichao, LIAN Yuguang, MU Yi, et al. Analysis of single main controlling factor affecting change of apparent resistivity of water-free faults in underground coal mine exploration[J]. China Energy and Environmental Protection, 2020, 42（12）: 106-111. （in Chinese with English abstract）
[13] 廉玉广,马志超,李江华,等. 岩石单轴加载破坏全过程波速变化特征研究[J]. 煤炭科学技术,2019,47(8):64 − 69. [LIAN Yuguang,MA Zhichao,LI Jianghua,et al. Study on variation characteristics of wave velocity in whole process of rock uniaxial loading failure[J]. Coal Science and Technology,2019,47(8):64 − 69. (in Chinese with English abstract) LIAN Yuguang, MA Zhichao, LI Jianghua, et al. Study on variation characteristics of wave velocity in whole process of rock uniaxial loading failure[J]. Coal Science and Technology, 2019, 47（8）: 64-69. （in Chinese with English abstract）
[14] 李江华,廉玉广,马志超. 受载岩体破坏全过程声波响应特征及工程意义[J]. 煤炭科学技术,2022,50(2):106 − 112. [LI Jianghua,LIAN Yuguang,MA Zhichao. Acoustic response characteristics and engineering significance of loaded rock during failure process[J]. Coal Science and Technology,2022,50(2):106 − 112. (in Chinese with English abstract) LI Jianghua, LIAN Yuguang, MA Zhichao. Acoustic response characteristics and engineering significance of loaded rock during failure process[J]. Coal Science and Technology, 2022, 50（2）: 106-112. （in Chinese with English abstract）
[15] 焦阳,廉玉广,李梓毓,等. 综合矿井物探技术在陷落柱探测中的应用[J]. 煤矿开采,2018,23(6):16 − 18. [JIAO Yang,LIAN Yuguang,LI Ziyu,et al. Application of integrated geophysical techniques in collapse column exploring[J]. Coal Mining Technology,2018,23(6):16 − 18. (in Chinese with English abstract) JIAO Yang, LIAN Yuguang, LI Ziyu, et al. Application of integrated geophysical techniques in collapse column exploring[J]. Coal Mining Technology, 2018, 23（6）: 16-18. （in Chinese with English abstract）
[16] 李忠华,梁影,包思远,等. 断层冲击地压的影响因素分析[J]. 中国地质灾害与防治学报,2020,31(3):126 − 131. [LI Zhonghua,LIANG Ying,BAO Siyuan,et al. Analysis on influence factors of the fault rock burst[J]. The Chinese Journal of Geological Hazard and Control,2020,31(3):126 − 131. (in Chinese with English abstract) LI Zhonghua, LIANG Ying, BAO Siyuan, et al. Analysis on influence factors of the fault rock burst[J]. The Chinese Journal of Geological Hazard and Control, 2020, 31（3）: 126-131. （in Chinese with English abstract）
[17] 张平松,欧元超,李圣林. 我国矿井物探技术及装备的发展现状与思考[J]. 煤炭科学技术,2021,49(7):1 − 15. [ZHANG Pingsong,OU Yuanchao,LI Shenglin. Development quo-status and thinking of mine geophysical prospecting technology and equipment in China[J]. Coal Science and Technology,2021,49(7):1 − 15. (in Chinese with English abstract) DOI: 10.13199/j.cnki.cst.2021.07.001 ZHANG Pingsong, OU Yuanchao, LI Shenglin. Development quo-status and thinking of mine geophysical prospecting technology and equipment in China[J]. Coal Science and Technology, 2021, 49（7）: 1-15. （in Chinese with English abstract） DOI: 10.13199/j.cnki.cst.2021.07.001
[18] 赵虎,张泉,谭建秋,等. 基于综合物探方法的公路不稳定边坡潜在滑面探测[J]. 中国地质灾害与防治学报,2022,33(1):117 − 122. [ZHAO Hu,ZHANG Quan,TAN Jianqiu,et al. Detection of the potential failure zones of the unstable slopes along highway using comprehensive geophysical methods[J]. The Chinese Journal of Geological Hazard and Control,2022,33(1):117 − 122. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2022.01-14 ZHAO Hu, ZHANG Quan, TAN Jianqiu, et al. Detection of the potential failure zones of the unstable slopes along highway using comprehensive geophysical methods[J]. The Chinese Journal of Geological Hazard and Control, 2022, 33（1）: 117-122. （in Chinese with English abstract） DOI: 10.16031/j.cnki.issn.1003-8035.2022.01-14
[19] 王瑞丰,温来福,程久龙,等. 高密度电法与瞬变电磁法联合勘查河北承德地区基岩裂隙水[J]. 地球科学与环境学报,2020,42(6):784 − 790. [WANG Ruifeng,WEN Laifu,CHENG Jiulong,et al. Joint detection of bedrock fissure water using high-density electrical method and transient electromagnetic method in Chengde area of Hebei,China[J]. Journal of Earth Sciences and Environment,2020,42(6):784 − 790. (in Chinese with English abstract) WANG Ruifeng, WEN Laifu, CHENG Jiulong, et al. Joint detection of bedrock fissure water using high-density electrical method and transient electromagnetic method in Chengde area of Hebei, China[J]. Journal of Earth Sciences and Environment, 2020, 42（6）: 784-790. （in Chinese with English abstract）
[20] 邓凤茹,赵艳. 基于无线电波坑道透视算法数学模型的研究[J]. 北华航天工业学院学报,2012,22(3):12 − 13. [DENG Fengru,ZHAO Yan. Study of mathematical model based on the algorithm of radio waves tunnels perspective[J]. Journal of North China Institute of Aerospace Engineering,2012,22(3):12 − 13. (in Chinese with English abstract) DENG Fengru, ZHAO Yan. Study of mathematical model based on the algorithm of radio waves tunnels perspective[J]. Journal of North China Institute of Aerospace Engineering, 2012, 22（3）: 12-13. （in Chinese with English abstract）
-
期刊类型引用(1)
1. 刘星宇，朱立峰，孙建伟，贾煦，刘向东，黄虹霖，程贤达，孙亚柯，胡超进，张晓龙. 沟谷型泥石流特征参数的等代面积递归精细求解. 西北地质. 2024(03): 272-284 . 
百度学术
其他类型引用(0)
计量
- 文章访问数: 1983
- HTML全文浏览量: 1123
- PDF下载量: 196
- 被引次数: 1
邮件订阅
RSS