A study of the development regularity and formation mechanism of ground subsidence in shallow coal seam mining of Yangchangwan coal mine, Ningxia
-
摘要:
以宁夏羊场湾煤矿Y110207工作面为研究对象，采用无人机遥感技术、野外调查与有限差分软件模拟方法研究浅埋煤层开采的地面塌陷类型、发育规律及其形成机理。①浅埋煤层开采地面塌陷以地表裂缝发育为主，地表破坏严重。②平行切眼裂缝间隔性出现，展布于整个工作面内，间隔距离为10~120m，局部裂缝形成错台高度约为15cm。平行顺槽裂缝为拉张型裂缝，发育在顺槽至外围一定范围。③采煤活动导致地下形成采空区，上覆岩层发生移动破坏，破坏区分为剪切破坏区、拉张破坏区及剪-拉破坏区，分别对压应力区、拉应力区和压-拉转化区。④当应力扰动传递至地表，应力值超过覆盖层抗拉强度时地表产生裂缝。随着工作面推进，覆岩内部裂缝带上行裂缝与地表下行裂缝贯通，形成错台。研究成果丰富了该区浅埋煤层的地面塌陷理论知识，为地面塌陷防治提供了理论依据。
Abstract:The mechanism of ground subsidence induced by coal mining is an important theoretical basis for the prevention and is also one of the hotspots in the study of geological environment. Subsidence type, features and the mechanism were studied by methods of the UAV remote sensing technology, field observation and numerical simulation software on the Y110207 working face in Yangchangwan coal mine, Ningxia. The results show that, during shallow coal seam mining, cracks and surface collapses are developed, which cause ground surface destruction. In the working face, the parallel cutting cracks are distributed with an interval distance between 10~120m. The parallel channel fracture is a tensile fracture, which is developed in a certain range from the trough to the periphery. The coal mining activities lead to the move and destruction of the overlying strata. The damage can be divides into three types, i.e., shear failure zone, tensile failure zone and shear-tensile failure zone, corresponding to compressive stress region, tensile stress zone and compression tensile transition zone. When the stress is more than the tensile strength of the covering layer transferring to the surface, the cracks will occur. With the advancing of the working face, the water flowing fracture zone will develop to the surface. After that cracks will perforate and the staggered platform will emerge. It enriches the theoretical knowledge of ground subsidence in shallow coal seam and provides a theoretical for the prevention.
-
地下水氮污染是世界性难题[1-2]，而中国90%的城市地下水不同程度遭受N等元素污染，地下水污染已呈现由点向面、由城市向农村扩展的趋势[3]。近年来，随着煤炭资源开发促进能源化工基地经济的发展，煤炭开采过程中工业污染物的排放极有可能对地下水水质安全造成严重的威胁[4-12]。90%的煤炭含氮量为0.52%~1.41%[13]。前人探讨了煤矸石回填对地下水氮的影响[14]、采煤沉陷积水区地表水与浅层地下水氮的影响[6]，对能源化工基地煤炭资源开发区地下水三氮污染影响因素尚未开展研究。
煤炭资源开发过程中，改变原有包气带中所含的潜在污染物改变了原有的循环模式，而与地表水联系最紧密的是浅层地下水，开采过程在地表、浅层地下水中形成了特有的转化模式[15]，尤其是在采动过程中可能形成的一系列波状地貌，如大量发育的裂隙、裂缝、阶地等。研究区为沙地区，会形成大量的波状起伏的丘陵地貌，对地表、包气带的结构及水、污染物的运移产生了潜在的影响，区域物源补给直接影响氮的含量。故本文立足于地下水径流系统和物源条件探讨转化机理。
宁东煤炭基地鸳鸯湖矿区是国家煤炭规划建设的重要矿区，梅花井井田为鸳鸯湖矿区的5个大型煤矿之一。宁东煤炭基地属2个含煤时代，石炭纪—二叠纪和侏罗纪煤田8个矿区共含煤38层，梅花井矿区主要产侏罗纪煤炭。煤炭开采过程中工业污染物可能含有大量的氮污染物，导致含水层氮污染，生态环境脆弱，水资源相对匮乏，煤炭开采区地下水资源的保护及矿区可持续发展具有重要意义。前人研究了梅花井矿区水文地质条件[16]、采矿对含水层水量及结构的破坏类型和程度[17-19]，而对煤矿区乃至周边的影响区水质尚未探讨。梅花井井田位于宁东夹山梁地貌单元，本文基于地下水补径排条件选择宁东夹山梁地貌为基本研究单元（包括了梅花井井田和周边环境），研究风积沙滩地区煤矿井开采地下水氮污染的程度和影响因素。
1. 研究区水文地质条件
夹山梁位于宁夏回族自治区灵武市以东33km处。梅花井井田位于夹山梁的中东部，鸳鸯湖背斜东翼中部，矿区面积78.96km2，开采深度1300~ 200m。属半沙漠低山丘陵地形。地表为沙丘掩盖，多系风成新月形和垄状流动沙丘。区内无常年地表径流，雨季降水及矿井排水在井田北部低洼地区形成盐碱湖，冬季干枯为碱滩，基岩被第四系风积沙覆盖。
夹山梁地区含水层按岩性组合特征及地下水水力性质、埋藏条件等，主要含水层由上而下划分为第四系孔隙潜水含水层和侏罗系直罗组、延安组砂岩裂隙-孔隙承压水含水层（图 1）。
![]()
图 1 水文地质图及样点分布Q—第四系；K—白垩系；J—侏罗系；T—三叠系Figure 1. The hydrological geological map and sample distribution（1）本区第四系厚2.75~17m，平均厚5.15m①，地下水主要赋存于风积沙、小型洼地及沟谷冲洪积层中。按地下水赋存条件，可分为风积沙潜水层、风积-冲洪积潜水层。地下水主要赋存于风积沙、小型洼地及沟谷冲洪积层中。含沙漠凝结水，地下水位埋深1.00~3.48m，多随地形起伏而异，水位、水量随季节变化，矿化度为0.8~8.03g/L，水化学类型主要为Na- HCO3、SO4 · HCO3- Na · Mg、Cl·SO4-Na·Mg。
（2）侏罗系含水层组是影响矿区煤矿开采的主要含水层，包括上侏罗统直罗组含水层、中统延安组含水层，厚度为13.01~618.62m。煤层厚度主要分布在6.85~337.45m之间，垂向上上部直罗组含水层富水性较强，对梅花井煤层开采影响较大，下部延安组含水层结构较致密，裂隙不发育，富水性较差，对煤层开采影响较小，水化学类型为SO4· Cl-Na、HCO3-Mg·Na·Ca或SO4·HCO3-Na·Mg。
（3）隔水层岩性多为中细砂岩与粉砂岩、泥岩互层，岩性致密，与煤层共同形成良好的隔水层。
2. 研究方法
采样点位于鸳鸯湖矿区夹山梁，采集第四系含水层水样10个和侏罗系含水层样品3个（图 1），利用GPS精确定位，潜水水位埋深一般为6~8m，侏罗系含水层水位埋深一般为65~100m。采样时间为2016年9月12日—9月17日，统测1次。水样采集、固定及保存均按照《水质采样样品的保存和管理技术规定》（HJ 493—2009）②进行，并及时运往实验室检测。依据水和废水监测分析方法[20]对浅层地下水中的氮分布及限制性特征进行分析；现场测试指标为水温，室内测试分析TDS、NO3-、NO2–、NH4+共4项。水化学测试由中国地质调查局西安地质调查中心实验测试中心完成。评价方法参照《矿山地质环境调查评价规范DD2014— 05》 8.3.2.2条款，地下水污染评价采用单项超标倍数法。
3. 结果
3.1 三氮空间污染特征
对夹山梁地区浅层地下水中主要离子及部分水化学指标的数据统计显示（表 1），本区浅层地下水pH值范围为7.13~9.63，浅层地下水普遍处于偏碱性环境，且酸碱度空间差异不大。矿化度变化范围为185.2~3278.22mg/L，表明矿化度指标空间分布差异较大。水质为淡水-咸水之间。其中一处水样中NO3–毫克当量百分数超过25%（表 1），对水化学类型产生较大影响。
表 1 地下水测试数据Table 1. The list of groundwater test data样号 地貌类型 水温/℃ pH 井水埋深/m 含水层类型 总矿化度mg/L NH4+mg/L NO3–mg/L NO2–mg/L 水化学类型 J1 风积沙滩 16 7.73 65 承压水 491.3 0.06 15.7 < 0.01 HCO3-SO4-Na-Mg-Ca J2 风积沙滩 15 8.53 70 承压水 569.32 0.07 20.2 < 0.01 HCO3-SO4-Na-Ca-Mg J3 风积沙滩 14 7.96 100 承压水 630.02 0.06 6.71 < 0.01 HCO3-Mg-Na-Ca J4 风积沙滩 18 7.71 0 泉水 826.62 0.08 80.8 0.74 HCO3-SO4-Na-Mg-Ca J5 风积沙滩 17 7.89 5 潜水 1368.28 0.12 200 < 0.01 SO4-HCO3-Cl-Na-Mg J6 风积沙滩 17 7.95 2 潜水 1038.91 0.11 234 0.92 HCO3-NO3-Na-Mg J7 风积沙滩 16 8.01 4 潜水 1538.07 0.11 67.1 0.49 HCO3-Cl-SO4-Na-Mg J8 风积沙滩 12 7.94 4.3 潜水 1668.22 0.1 171 2.01 Na-HCO3 J9 冲积平原 7.13 1 潜水 185.2 0.06 7.97 < 0.01 HCO3-SO4-Ca-Mg J10 冲积平原 17 8.51 潜水 1622.24 0.08 43.2 0.84 Cl-SO4-HCO3-Na-Mg J11 冲洪积平原 16 9.63 2 潜水 1549.54 0.1 6.3 0.69 Cl-SO4-Na-Mg J12 冲洪积平原 8.19 3.5 潜水 3278.22 0.1 79.6 < 0.01 Cl-SO4-Na-Mg 由表 1、图 2可以看出，研究区13组水样中三氮NH4+、NO3–、NO2–水化学特征，含量范围分别为0.06~0.12mg/L、4.67~234mg/L、 < 0.01~2.01mg/L，与国家地下水质量标准Ⅲ类水质限值对比，氨氮均检出，但未污染；NO2–均检出，其中6组含量为重度或极严重污染；NO3–的7组样品未检出，其余5组NO3–污染级别为中度、轻度污染。空间上无论矿权范围还是矿权外，污染样点均存在，不同的是NO3–污染主要在第四系潜水含水层中，而NO3–污染已经迁移至承压水中。超标样点占调查样点的75%。
![]()
图 2 宁东梅花井地下水三氮污染评价Figure 2. Evaluation graph of groundwater three nitrogen pollutants in Meihuajing field, Ningdong3.2 矿区地下水氮物源
本次采集了煤矸石和区域土壤，pH为碱性。煤矸石产生量约占煤炭开采量的5%~25%[20]。从表 2可知，区域土壤氮含量平均值为346mg/kg，氮含量为5级，氮养分缺乏；煤矸石氮含量平均值为1213mg/kg，氮含量为3级，氮养分足量。矸石中含有丰富的N元素，即矿业开发为地下水氮污染提供了丰富的氮源。刘钦甫等[21]在煤矸石中氮溶出的动态淋滤实验中，得出煤矿区地表矸石中含有吸附状态的硝态氮，而这部分氮较容易被水所溶出，特别是中性水比酸性水更容易使其溶出。在3种形态的氮中，硝态氮明显占优势，且其变化趋势与总氮一致。
表 2 煤矸石、表土中N元素含量Table 2. The content of nitrogen in gangue and soil物源类型 样品编号 N pH 附近布设井孔 GS1 301 8.7 GS2 1556 7.58 研究区内 GS3 1406 8.02 GS4 257 6.78 GS5 519 8.31 GS6 424 8.70 GS7 899 8.15 GS8 859 8.23 GS9 1523 8.12 煤矸石 GS10 173 9.22 GS11 2382 8.15 研究区外 GS12 1668 7.78 GS13 703 8.39 GS14 2331 7.89 GS15 3123 7.52 GS16 2131 7.13 GS17 1595 7.14 GS18 770 10.97 GS19 418 8.36 煤矸石平均值 1213 8.17 土壤 T1 472 9.21 J3 T2 290 9.34 J5 T3 440 9.26 J2 T4 386 9.41 J1 T5 170 8.97 J7 T6 320 8.95 J11 土壤平均值 346 9.19 中国土壤普查技术含量分级[22] 6级 < 50 5级 50~750 4级 750~1000 3级 1000~1500 2级 1500~2000 1级 >2000 粉煤灰场旁水样取自电厂灰渣贮灰场旁边的水塘，废水排放标准暂无对硝酸盐的规定，本文参照国家地表水环境质量标准（GB3838—2002）和国家地下水环境质量标准（GB/T 14848—93）Ⅲ类水限值。从表 3可知，淋滤液、河水、湖水、降雨中含有硝酸盐和氨氮成分，但相对国家标准其水质是安全的。
表 3 粉煤灰淋滤液、地表水、雨水中N元素含量Table 3. The content of nitrogen in fly ash filter, surface water and rain water水样类型 样品编号 NH4+ NO3– NO2– 粉煤灰场旁水样 F1 0.06 10.8 < 0.01 F2 0.07 8.09 < 0.01 F3 0.07 11.4 < 0.01 河水 H1 < 0.02 1.36 < 0.01 H2 < 0.02 4.09 < 0.01 湖水 H3 0.07 1.9 < 0.01 雨水 Y1 0.04 2 0 国家地表水环境质量标准（GB38382002）Ⅲ类水限值 1 10 研究区为牧区，大量的牧羊粪便随着地表或雨水进入地下水，施用动物粪便等的地区，地表水回补后的地下水NO3–值较高；地下水中N元素的浓度逐年增加[22-24]。
4. 讨论
地下水中三种主要的无机氮形态NH4+、NO3–、NO2–，以NO3–、NO2–为主，其中地下水中的NO3-占无机氮的比例最高，氨氮未出现超标，超标成分为NO3–、NO2–。那么地下水中NH4+不超标而NO3–、NO2–超标的影响因素是什么？
首先，从地下水补给-径流-排泄条件考虑：①垂向补给条件，由前述可知，目前淋滤液、河水、湖水中三氮相对国家标准水质是安全的，引起地下水硝酸盐和亚硝酸盐超标的可能性较小。但是由于有丰富的物源，在地表水或降雨长期淋滤作用下，沿垂向可能引起地下水中硝酸盐和亚硝酸盐的超标。②受地形地貌影响明显，地下水补给区位于夹山梁和布朗山，研究区只有J9井孔位于丘陵高地，水质良好，而处于经过长期的蒸发作用和溶滤作用的排泄区的样点，如J8、J10、J11井点中，NO3和NO2–含量增加，出现超标现象，J8井点中的NO3–、NO2–超标倍数分别为99.5倍、7.55倍，J10井点中的NO3–、NO2–超标倍数分别为41、1.16倍，J11井点中的NO3–超标倍数为33.5倍。③煤矿采掘扰动改变含水层围岩和径流条件，也会促进N的释放。如J1与J2井点具相似的地形条件，不同的是J2位于矿权范围，井孔中NO3–含量出现超标现象，而在矿权范围外的J1井水质良好，煤炭资源在开采过程中，煤层、围岩中的氮化物矿物与氧气和水接触，在微生物的催化作用下，经过一系列复杂的地球化学反应，可促进N元素的释放。由于矿井水在地下与围岩裂隙水存在着一定的水力联系，这些含氮物质的释放将会对水环境造成严重的威胁[25-28]。
其次，地下水水文地球化学特性的众多研究表明，地下水氮污染形式主要是NO3–的污染，其是国内外最普遍、污染面积最大的地下水污染问题[29-31]，在NH4+、NO3–、NO2–系统中三氮变化关系通常总是向NO3–转化，而NO2–浓度增高的罕见情况也只是具有暂时性意义[6]。因此长期以来，NO2–在潜水中出现增高甚至超标，取决于复杂的地下水中氮的转化过程，主要包括有机氮的矿化作用、NH4+吸附作用、硝化作用、反硝化作用、异化还原等[31]。同时温度、pH、Eh、土层介质厚度、土壤透气性等也是重要的影响因素[25, 32-35]。Sitaula等[35]认为，硝化反应的最佳pH值范围为8~8.4，pH值低于7时，硝化速率明显降低，低于6或高于9.6时，硝化反应几乎停止，研究区pH值范围为7.73~9.63，是硝化反应的最佳区域。王晓娟等[36]认为，在氮的转化过程中，细菌起着重要的作用，几乎所有微生物在好氧环境转变为厌氧环境后均可参与硝化过程；有研究表明，某些反硝化细菌在好氧条件下也可以进行反硝化[37]，一定温度范围内，氮矿化（氮由有机态转化为无机态NH4+或NH3的过程）随温度的升高而升高，随土层深度增加而降低，随土壤通透性的降低而降低，这也可能是潜水中NO3–、NO2–含量高而承压水中只有NO3–含量高的原因。
5. 结论
（1）NO3–、NO2–与国家地下水质量标准Ⅲ类水质限值对比污染严重，NO2–污染主要在第四系潜水含水层中，而NO3–在潜水和承压水中均有污染。研究区超标样点占调查样点的76.92%。根据物源和氮污染空间分布情况推测，煤矿区氮的污染很可能是研究区高地球化学背景引起的。
（2）NO3–在研究区浅层地下水中的相对含量较高，NO3–毫克当量百分数超过25%的水样，对水化学类型产生影响。
（3）NO3–、NO2–的影响因素很大程度受限于煤矿开采、地形地貌条件、垂向补给及水文地球化学条件。低山丘陵高地及煤矿开采未影响到的地方，水质较好。温度、pH、Eh、土层介质厚度、土壤透气性等也是重要的影响因素。
致谢: 地表裂缝的野外调查工作由羊场湾煤矿地测科张仲杰等工作人员协助完成, 在此表示衷心的感谢。 -
![]()
图 4 Y110207工作面地表裂缝展布
Figure 4. Surface cracks distribution map of Y110207 working face
![]()
图 6 覆岩竖向位移及竖向应力云图
Figure 6. Vertical displacement and vertical stress nephogram of overlying strata
![]()
图 7 覆岩及地表破坏分区
A—走向；B—倾向
Figure 7. Destructive zoning map of overlying strata and ground
![]()
图 10 地表走向水平应力云图及位移曲线
Figure 10. Horizontal movement map of parallel mining direction on the surface
![]()
图 11 地表倾向水平应力云图及位移曲线
Figure 11. Horizontal movement map of vertical mining direction on the surface
表 1 MD4-1000四旋翼低空无人机指数
Table 1 MD4-1000 four indexes of the UAV
性能指标 参数 爬升速率 7.5m/s 最大功率 1000W 飞行时间 ＜50min/电池 飞行半径 5000m 遥控距离 5000m 巡航速度 15.0m/s 机身自重 2650g 任务载荷 2000g(最大) 飞行高度 1000m 抗风能力 ＜12m/s 
下载: 导出CSV
表 2 岩石力学参数
Table 2 Rock mechanical parameter table
岩性 自重密度/(g·cm3) 体积模量/GPa 切变模量/GPa 粘结力/MPa 内摩檫角/° 抗拉强度/MPa 松散层 1780 2.1 0.4 0.9 28 0.6 粉砂岩 2600 3.0 1.1 1.2 33 0.8 中砂岩 2650 4.6 1.2 1.0 35 1.1 粉砂岩 2600 3.3 1.8 1.3 33 0.8 细砂岩 2750 6.4 2.2 2.0 40 1.0 粉砂岩 2600 3.3 2.5 1.3 33 0.8 砂岩 2900 4.6 1.7 2.0 40 1.3 煤层 2030 2.2 1.0 1.1 31 0.6 底板砂岩 2860 3.9 3.3 3.5 35 1.0
下载: 导出CSV
-
何芳, 徐友宁, 袁汉春, 等.煤矿地面塌陷区的防治对策[J].煤炭工程, 2003, 7:10-13. http://d.old.wanfangdata.com.cn/Periodical/mtgc200307003 初影.采煤诱发地表裂缝数值模拟研究[D].辽宁工程技术大学硕士学位论文, 2009. http://cdmd.cnki.com.cn/Article/CDMD-10147-1011025019.htm 蔡怀恩, 侯恩科, 张强骅, 等.黄土丘陵区房柱式开采地表塌陷特征及机理分析——以陕北府谷县新民镇小煤矿为例[J].地质灾害与环境保护, 2010, 21(2):101-104. doi: 10.3969/j.issn.1006-4362.2010.02.023 蔡怀恩.彬长矿区地面塌陷特征及形成机理研究[D].西安科技大学硕士学位论文, 2008. http://cdmd.cnki.com.cn/Article/CDMD-10704-2008145586.htm 贺卫中, 向茂西, 刘海南, 等.榆神府矿区地面塌陷特征及环境问题[J].煤田地质与勘探, 2016, 44(5):131-135. doi: 10.3969/j.issn.1001-1986.2016.05.025 赵兵朝, 刘宾, 王建文, 等.柠条塔煤矿叠置开采地表岩移参数分析[J].煤矿安全, 2016, 47(9):213-216. http://d.old.wanfangdata.com.cn/Periodical/mkaq201609060 赵兵朝, 刘飞, 凡奋元, 等.黄土沟壑区下斜交叠置开采地表下沉系数研究[J].矿业安全与环保, 2016, 43(5):54-57. doi: 10.3969/j.issn.1008-4495.2016.05.013 杨帆, 余海锋, 郭俊廷.采动地表裂缝形成机理的数值模拟[J].辽宁工程技术大学学报, 2016, 6:566-570. doi: 10.11956/j.issn.1008-0562.2016.06.002 吴侃, 胡振琪, 常江, 等.开采引起的地表裂缝分布规律[J].中国矿业大学学报, 1997, 2:56-59. http://cdmd.cnki.com.cn/Article/CDMD-10290-2010280128.htm 胡振琪, 王新静, 贺安民.风积沙区采煤沉陷地裂缝分布特征与发生发育规律[J].煤炭学报, 2014, 39(1):11-18. http://d.old.wanfangdata.com.cn/Periodical/mtxb201401002 胡振琪, 龙精华, 王新静.论煤矿区生态环境的自修复、自然修复和人工修复[J].煤炭学报, 2014, 39(8):1751-1757. http://d.old.wanfangdata.com.cn/Periodical/mtxb201408048 杜善周.神东矿区大规模开采的地表移动及环境修复技术研究[D].中国矿业大学(北京)博士学位论文, 2010. http://cdmd.cnki.com.cn/Article/CDMD-11413-2010240648.htm 邓喀中, 王刘宇, 范洪冬.基于InSAR技术的老采空区地表沉降监测与分析[J].采矿与安全工程学报, 2015, 32(6):918-922. http://d.old.wanfangdata.com.cn/Periodical/ksylydbgl201506008 王瑞国.基于WorldView-2数据的乌东煤矿地质灾害遥感调查及成因分析[J].国土资源遥感, 2016, 28(2):132-138. http://d.old.wanfangdata.com.cn/Periodical/gtzyyg201602022 刘辉, 何春桂, 邓喀中, 等.开采引起地表塌陷型裂缝的形成机理分析[J].采矿与安全工程学报, 2013, 30(3):380-384. http://d.old.wanfangdata.com.cn/Periodical/ksylydbgl201303011 余学义, 邱有鑫.沟壑切割浅埋区塌陷灾害形成机理分析[J].西安科技大学学报, 2012, 32(3):269-274. doi: 10.3969/j.issn.1672-9315.2012.03.001 刘辉, 刘小阳, 邓喀中, 等.基于UDEC数值模拟的滑动型地裂缝发育规律[J].煤炭学报, 2016, 41(3):625-632. http://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201603015.htm 黄庆享, 张文忠.浅埋煤层条带充填保水开采岩层控制[M].北京:科学出版社, 2014. 周文生, 吴振宇, 刘海燕.无人机遥感在矿山地质环境调查中的应用[J].地下水, 2014, 2:128-129. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=DXSU201402053&dbname=CJFD&dbcode=CJFQ 张启元.无人机航测技术在青藏高原地质灾害调查中的应用[J].青海大学学报, 2015, 2:67-72. http://d.old.wanfangdata.com.cn/Periodical/qhdxxb-zr201502012 侯恩科, 首召贵, 徐友宁, 等.无人机遥感技术在采煤地面塌陷监测中的应用[J].煤田地质与勘探, 2017, 6:102-110. http://d.old.wanfangdata.com.cn/Periodical/mtdzykt201706017 赵坤阳.煤矿地下开采诱发地表裂缝与导水裂缝分布规律预测[D].中国地质大学(北京)硕士学位论文, 2015. http://cdmd.cnki.com.cn/Article/CDMD-11415-1015391331.htm 余学义, 张恩强.开采损害学[M].北京:煤炭工业出版社, 2010.
计量
- 文章访问数: 3131
- HTML全文浏览量: 579
- PDF下载量: 2596
