第一作者简介: 王海庆(1980-),男,博士,高级工程师,主要从事遥感地质应用方面的研究。Email:whq0705@126.com。
为了研究采矿沉陷危害,选择山东省济宁市东部煤矿矿集区为研究区,采用多期光学遥感数据与历史资料相结合的方法,辅以必要的现场调查和走访,进行采矿沉陷调查与危害性研究。结果表明: ①研究区内采矿沉陷危害严重,截止到2013年,至少破坏了25.095 km2的土地,并造成24个村庄被迫搬迁; ②总体来说,研究区内采矿沉陷灾害发展迅速,从2009—2013年,采矿沉陷积水面积共增加4.747 km2,2013年相对于2009年的增长率为23.33%; ③就单个采矿沉陷来说,开始阶段来势猛、发展快、危害大; 但快速发展之后,扩展的速度将逐步变缓,并逐渐过渡到自然平衡; ④光学遥感技术完全可以用于采矿沉陷调查和危害性研究,并且其宏观、廉价等特点较好地弥补了地面调查的不足。
In this paper , aimed at the study of the hazards caused by the mining subsidence using remote sensing images,according to the research objective,the coal mine ore concentration located in the east of Jining City, Shandong Province of China was chosen as the study area. Multi stage optical remote sensing images and historical data were used, and necessary field survey were supplemented also. The results show that: ①The mining subsidence in the study area is very serious. By 2013, at least 25.095 km2 of land was destroyed, 24 villages were forced to move; ②In generally, mining subsidence developed rapidly in the study area. From 2009 to 2013, 4.747 km2 of mining subsidence water area was increased in total, and the growth rate was 23.33%; ③ For a single mining subsidence, it is fierce, fast development, great hazard in the beginning. But after the rapid development, the expansion rate will gradually slow down, and gradually transit to the natural balance; ④Optical remote sensing technique can be used to research and investigation of mining subsidence hazard, and the characteristics such as macroscopic and cheap could make up for the lack of ground survey.
煤矿资源是我国重要的矿产资源。煤矿资源的开采利用已为区域经济发展做出了重要贡献, 但同时也造成了严重的地面沉陷, 破坏了大量的农田、建筑和道路, 严重影响了当地人民群众的生活。
国内对采矿沉陷的研究已取得不少的成果, 如阎跃观等[1]研究了急倾斜多煤层开采条件下围岩破坏机理并将地表沉陷分为4种类型区域; Wu等[2]研究了山东省龙口地区煤矿开采引起的地面沉陷及其对环境的影响; 李成尊等[3]对采矿塌陷的遥感特征进行了研究; 王晓红等[4]对采矿塌陷的遥感可识别性进行了研究; 王钦军等[5]探讨了矿山塌陷的遥感识别与边界提取方法; 姚丹丹等[6]和高小六[7]探讨了使用D-InSAR技术监测煤矿区沉陷的方法和效果; 王勤发等[8]对万年县老白马煤矿采空塌陷的现状进行了调查和预测分析; 王创业等[9]在神东矿区锦界煤矿对煤矿开采的松散层厚度、采厚、覆岩综合硬度等影响沉陷的因素进行了研究; 迟占国等[10]在平朔井工一矿研究了煤矿区地表沉陷的规律; 杨成奎[11]研究了地面沉陷的预测方法; 王永辉等[12]和龙四春等[13]探讨了煤矿沉陷灾害的预测方法; 赵晓霞等[14]将采煤塌陷损毁耕地分为轻度、中度、重度3个级别; 刘哲荣等[15]研究了采煤沉陷区土壤性质的演变; 黄晓娜等[16]研究了煤矿塌陷区不同复垦年限土壤颗粒组成的分形特征; 柴华彬等[17]对煤矿塌陷区地基稳定性与承载力的研究现状进行了分析。上述研究虽然广泛涉及围岩破坏机理、危害影响、遥感识别与监测方法、预测分析、土壤性质演变和地基稳定性等方面, 但对应用高分辨率遥感资料研究地面沉陷危害的报道较少。笔者[18]曾利用2009年获取的GeoEye-1高分辨率遥感数据对济宁市煤矿矿集区地面沉陷的现状进行过遥感调查, 并利用多期遥感数据分别对兴隆庄煤矿和龙堌煤矿周边的地面沉陷面进行过调查与监测[19, 20]; 汪宝存等[21]利用遥感数据对开滦煤矿地面塌陷积水进行了动态监测; 许冬等[22]运用遥感与地理信息技术探讨了济宁煤矿区地表塌陷积水的时空演变。但上述研究只注重对采矿沉陷积水面积的调查, 对采矿沉陷的危害性研究不足。本文采用多期遥感数据, 结合地形图和地质图等历史资料, 辅以实地调查验证, 利用ArcGIS平台开展研究区内的采矿沉陷调查和危害性研究。调查和研究过程中充分发挥了遥感技术宏观、真实、高效的特点, 结果真实可靠。
选择山东省济宁市东部煤矿矿集区为研究区(图1)。
研究区具备下述特点和研究条件: ①地下开采煤矿集中, 在较小的范围内有较大的煤矿开采和生产规模; ②采矿沉陷严重, 已经形成了一定的采矿沉陷灾害, 并且采矿沉陷持续扩展, 其危害不断增加; ③有多期可供选择利用的遥感数据。
研究区位于济宁市东部, 行政区域上涉及曲阜市、兖州市和邹城市3个县级市, 面积396 km2。区内地势平坦, 为典型的平原地貌。区内有多个国有大型煤矿, 主要涉及兴隆庄、杨庄、杨村、东滩、鲍店、横河和南屯共7个煤矿(其中南屯煤矿仅涉及其北部); 这些煤矿都采用地下开采, 生产规模较大。以往研究结果表明, 该区采矿沉陷严重, 而且扩展快, 危害大[12, 14, 18, 19, 20, 22], 是典型的煤矿采矿沉陷区。
采矿沉陷是指在矿产资源开采利用时造成的地面沉陷, 不同于采矿塌陷。采矿塌陷往往形成地表高程急剧变化的塌陷坑, 且面积较小(往往以m2为单位计算), 在遥感图像中受阳光照射角度和卫星遥感图像拍摄角度的影响, 可以见到较为明显的阴影, 从而易于识别。采矿沉陷的地面高程变化相对较缓, 且面积较大, 在光学遥感图像中难以直接识别, 但可根据采矿沉陷形成的积水坑间接识别。采矿沉陷作为一种灾害, 如果没有积水, 仅使沉陷区的土地丧失部分使用功能(如不再具有建房居住的功能), 但仍可生长草木。而在济宁地区, 由于人口密度大、耕地资源严重紧张, 当地居民不顾危险, 仍在采矿沉陷的危险区内耕种, 直至积水淹没耕地才不得已而放弃耕种。本文着重研究采矿沉陷的危害, 侧重于研究有积水的采矿沉陷区, 故对采矿沉陷的遥感识别也将从采矿沉陷积水坑入手。
采矿沉陷积水坑具有独特的影像特征, 在高空间分辨率遥感图像中易于识别, 完全可以与露天的采沙坑、采土坑区别开[18]。其主要特征有: ①边界比较圆滑, 不会出现明显的拐角; 但采沙坑、采土坑的边界往往较平直, 有明显的拐角。②边坡比较平缓; 但采沙坑、采土坑的边坡往往较陡(甚至直立)。③形成和发育受地下采空区的控制, 往往可间接反映地下采空区的空间形态(比如采煤巷道的走向与规模)。另外, 通过对比历史资料与现场调查, 还可以将历史上形成的湖泊等地表水体与地面沉陷积水坑区别开。但需要注意的是, 某些采矿沉陷可能已被治理和利用, 其积水坑也会出现较平直的边界和明显的拐角。因此, 在利用光学遥感手段进行采矿沉陷积水坑调查时, 需要对多个识别标志进行综合判断分析, 必要时应进行相应的现场调查和走访。
根据多个识别标志进行综合判断分析, 从研究区的模拟真彩色遥感图像中圈定出采矿沉陷积水坑, 并通过野外调查验证进一步修改完善其识别结果。图3为研究区内采矿沉陷积水坑的分布情况。
经研究发现, 研究区内采矿沉陷积水坑分布广、面积大、危害严重。采矿沉陷积水坑的面积在逐年扩大, 在2009年获取的遥感图像中, 采矿沉陷积水面积为20.349 km2; 在2010年的遥感图像中面积为22.385 km2, 比2009年增加了2.036 km2, 增长率超过了10%; 在2011年的遥感图像中变为23.090 km2, 又比2010年增加了0.705 km2, 相对于2010年的增长率为3.15%; 在2012年的遥感图像中, 采矿沉陷积水面积为24.097 km2, 比2011年增长了1.007 km2, 相对于2011年的增长率为4.36%; 在2013年的遥感图像中, 采矿沉陷积水面积为25.095 km2, 比2012年增长了0.998 km2, 相对于2012年的增长率为4.14%(图4及表2)。
从2009— 2013年, 根据获取的遥感图像分析, 采矿沉陷积水面积共增加了4.747 km2, 2013年相对于2009年的增长率为23.33%。在增加的采矿沉陷积水面积方面, 兴隆庄煤矿(矿山范围内, 下同)增加最多, 为1.758 km2; 其次为鲍店煤矿, 为1.401 km2。在2013年相对于2009年的增长率方面, 杨村煤矿增长最快, 达56.91%; 其次是鲍店煤矿, 为31.02%。
根据调查得到的采矿沉陷积水面积, 结合矿山面积, 可计算出采矿沉陷积水比率。该比率可在一定程度上反映采矿沉陷的严重程度。研究区内采矿沉陷积水比率最大的是兴隆庄煤矿, 达17.20%; 其次是鲍店煤矿, 为15.99%(表3)。
采矿沉陷是一个动态变化过程, 采矿沉陷积水坑是采矿沉陷的结果之一。即有采矿沉陷未必有采矿沉陷积水坑, 只有当沉陷达到一定的程度, 使地表高程低于地下水位时, 才能形成有长期积水的采矿沉陷积水坑。即使如此, 在积水区域与非沉陷区之间也有一定宽度的过渡地带。对于过渡地带的宽度, 笔者也开展了一些调查和研究, 然而从目前的研究进展来看, 还没有统一的规律性数据或公式可以推荐, 但可明确的是其影响因素主要有: 采空区埋深、采空区高度、采空区宽度、采空区至地表的物质结构及其力学性质等。另外, 采矿沉陷积水面积也受降水量的影响, 连续性降雨和特大暴雨会使采矿沉陷积水面积扩大。研究区内地势非常平坦, 没有自然因素造成的地表高程变化, 故除非发生洪涝灾害使区内地面普遍被淹, 采矿沉陷积水坑的范围不会超出采矿沉陷的范围。在本文使用的遥感图像中, 未发现洪涝灾害的迹象; 气象资料也未显示在获取遥感数据之前短期内有长期连续性降雨和特大暴雨的记录, 故本次研究可以排除气象因素的影响。因此, 采矿沉陷的范围往往大于采矿沉陷积水坑的范围。本文遵循保守的原则, 利用采矿沉陷积水坑的面积估算采矿沉陷破坏的土地面积, 虽然不足以全面反映采矿沉陷破坏土地的实际面积, 但至少可作为一种参考, 并且也符合当地居民珍惜每一寸土地、辛勤耕种的实际情况。
根据获取的遥感图像, 本文认为研究区内的采矿沉陷至少破坏了25.095 km2的土地。其中, 兴隆庄煤矿(矿山范围内, 下同)破坏的土地最多, 为9.675 km2; 其次是东滩煤矿, 为6.106 km2; 鲍店煤矿为5.917 km2; 横河煤矿为1.282 km2; 南屯煤矿北部为1.212 km2; 杨村煤矿为0.863 km2; 杨庄煤矿为0.040 km2。
就单个采矿沉陷(对应单独的地下采空区)而言, 开始阶段来势猛、发展快、危害大, 往往在短时间内破坏大量土地及土地上的道路、房屋、庄稼和树木等; 快速发展后, 扩展速度将逐步变缓, 逐渐过渡到自然平衡状态。以兴隆庄镇东侧(大施村东北)的采矿沉陷(图5)为例, 2009年, 有一条NNW走向的采矿沉陷积水坑, 长1 070 m(图5(a)); 2010年, 该采矿沉陷积水坑急剧向北扩展, 长度扩展为1 953 m, 积水淹没了兴隆庄煤矿向外运煤的铁路(图5(b)中呈圆弧状的线状地物); 2011年相对2010年没有明显变化(图5(c)); 2012年, 该采矿沉陷积水坑仍没有明显的变化, 但其西侧出现了一个新的采矿沉陷积水坑, 长1 589 m(图5(d)); 2013年相对2012年没有明显变化(图5(e))。
采矿沉陷在破坏大量土地的同时, 也迫使多个村庄搬迁。根据多期遥感图像, 结合历史资料, 可以很容易地识别出已消失的村庄。但要识别出因采矿沉陷消失的村庄, 还需要综合利用煤矿开采和采矿沉陷等信息。为了区别因土地耕作集约化、并村合集、新农村建设等非避灾性搬迁的村庄, 笔者对室内解译结果进行了实地走访调查, 以进一步确认村庄搬迁的原因。如图6所示, 位于左上角的四新村所在地, 在2009年的遥感图像中已经沦为采矿沉陷积水坑; 位于右下角的后屯村所在地, 在2009年的遥感图像中已位于采矿沉陷积水坑的边沿, 村庄已基本搬迁完毕, 但仍依稀可见村庄的痕迹, 而在2013年的遥感影像中已经完全沦为采矿沉陷积水坑; 位于右上角的程家庄、东程家庄和西程家庄3个村庄所在地, 在2009年的遥感图像中还是完好的村庄(但在其南部、西部和西北部都有严重的采矿沉陷), 在2013年的遥感图像中已完全消失。
对照历史资料, 结合采矿沉陷积水坑分布情况, 圈定了24个因采矿沉陷消失的村庄, 并在现场调查和走访中得到了证实。24个村庄中, 位于兴隆庄煤矿矿山范围内的有9个, 东滩煤矿矿山范围内有6个, 鲍店煤矿矿山范围内有5个(表4)。
采矿沉陷灾害严重破坏了当地的生态环境, 大量土地被损毁, 多个村庄被迫搬迁, 不少建筑倒塌或被淹没(图7), 多条道路遭毁坏, 给采矿沉陷区原居民的生产和生活造成了灾难性影响, 也给当地的社会治安和安全稳定埋下了隐患。
根据采矿沉陷遥感调查结果, 笔者分析了研究区内采矿沉陷扩展趋势, 就单个采矿沉陷区域研究了其扩展方向。根据历年发展变化情况, 尤其是在近2 a快速发展的区域, 圈定了7处危险区(图3), 预测这些区域可能会发生新的采矿沉陷。在这7处危险区内, 涉及3个村庄, 预测这些村庄将需要搬迁。以图5所示的兴隆庄镇东侧采矿沉陷区为例, 从2009―2013年的发展变化来看, 采矿沉陷正在向SW方向快速发展。故笔者预测该采矿沉陷区域的西南部为危险区, 有可能发生新的采矿沉陷; 且其西南部的大施村处在危险区域内, 预测该村庄需要搬迁。
总体来说, 采矿沉陷受地下采空区控制, 采矿沉陷积水坑的形态和规模也在一定程度上反映了地下采空区的形态和规模。若要根治采矿沉陷, 仍需从地下采空区入手。建议在矿山开采之前制订切实有效的采空区治理方案, 并在矿山开采过程中严格执行。
1)所选择的研究区内采矿沉陷危害严重, 截止到2013年, 至少破坏了25.095 km2的土地, 并迫使24个村庄搬迁, 毁坏多条道路, 给沉陷区原居民的生产和生活造成了灾难性影响, 也给当地的社会治安和安全稳定埋下了隐患。
2)总体来说, 研究区内采矿沉陷灾害发展迅速。2009―2013年间, 根据获取的遥感图像分析, 采矿沉陷积水面积共增加4.747 km2, 2013年相对于2009年的增长率达到23.33%。
3)就单个采矿沉陷(对应单独的地下采空区)来说, 刚开始来势猛、发展快、危害大, 但快速发展后, 扩展的速度将逐步变缓, 并逐渐过渡到自然平衡。
4)本文采用的光学遥感技术具有宏观、廉价、快速、高效的特点, 可以用于采矿沉陷的危害调查和研究, 弥补常规地面调查的不足。
然而, 关于采矿沉陷危害的研究仍有一些重要问题需要进一步开展。比如: 本文中提到的在积水区域与非沉陷区之间也有一定宽度的过渡地带, 如何能够更加科学合理、符合实际地计算该宽度, 还需要深入探讨和研究。另外, 如何权衡矿产资源开采的利与弊?如何评价矿产资源开采利用带来的经济效益与采矿沉陷等矿山地质灾害的危害对社会整体发展的影响?是需要长期探讨的难题。限于笔者收集到的资料有限, 未能在这方面开展工作, 期待有志学者深入研究。
致谢: 本文研究过程中得到了杨金中、王晓红、周英杰、汪洁、王斌、姚维岭、李飞等同仁的帮助, 在此表示衷心的感谢!
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|
[21] |
|
[22] |
|