受人类活动与气候波动的共同干扰,土地退化正成为全世界面临的重大环境问题[1],遏制土地退化对保障生态系统服务及人类福祉至关重要[2]。煤炭资源开发过程会对土地造成直接损毁,如挖损、压占、塌陷[3-4]等,显著改变了生态环境要素的表现、组成和结构,严重制约矿区生态系统的完整性[5],进而诱发矿区土地退化[6]。另外,当前中国能源重心已转移至地处干旱半干旱区、生态环境脆弱的西部省份[7],高强度的煤炭开采可能引发比中部、东部地区更为严重的土地退化[8]。
国外对土地退化主要围绕碳源/碳汇[9-10]、植被恢复演替[11]、退化土地重建[12]等方面开展研究,在方法上包括遥感识别与解译[13-15]、土壤理化性质测定实验[9]等。但土地退化过程复杂,实验测定多采用小尺度采样,且监测点密度不足,难以实现实时动态监测。GLASOD,ASSOD,RUSSIA等[16]方法依据区域的异质性,结合土地退化成因、特点、主要表现等划分不同类型。国内对矿区土地退化的研究多集中于因素、分类、调查[8,17],或围绕土壤质量、生态健康、景观尺度等方面[18-19],以及利用3S技术开展土地退化的较大尺度监测[20-21]。
土地退化具有复杂性、区域性、动态性等特点,且导致土地退化的因素具有多样性。因此人为干扰生态系统内各环境因子及其生态学过程都可能引起土地退化[22]。同时,矿区土地退化受气候波动与开采扰动[23-24]等自然及人为因素的共同影响,表现为土壤侵蚀、沙漠化/石漠化、植被退化等。结合科技基础性工作专项项目“西部重点矿区土地退化因素调查”的部分成果,统筹考虑开采类型、地域分布、矿区的典型性与代表性等条件,从黄土高原区、西南山地丘陵区、北方草原区和西北干旱区选择12个重点煤矿区(图1)开展煤矿区土地退化研究,并将煤矿区土地退化定义为:以受煤炭资源开发活动的扰动为主,兼受自然因素与其他人为活动干扰等多种营力作用的影响,矿区及其周边土地结构遭受损毁,质量下降及生产力衰退的演替过程。其中,黄土高原区包括宁夏石嘴山矿区、灵武矿区,陕西榆神矿区、旬耀矿区,甘肃华亭矿区,内蒙古神东矿区和准格尔矿区;西南山地丘陵区包括贵州盘县矿区和云南小龙潭矿区;北方草原区包括内蒙古平庄矿区和胜利矿区;西北干旱区包括新疆准东矿区。
采用现场调查和遥感解译方法,提取矿区2013,2015,2017年的挖损、压占和塌陷3类土地损毁信息,分析土壤侵蚀、沙漠化/石漠化、植被退化等土
图1 西部重点煤矿分区
Fig.1 Location of 12 coal mine areas in Western China
地退化的主要表现,判断各煤矿区土地退化的关键因素,建立西部重点煤矿区土地退化综合评价模型并按照退化程度分等定级,为西部煤炭资源的保护性开发提供决策依据。
(1)土地损毁信息提取
露天矿山为露天采场挖损、排土场压占;井工矿山为煤矸石压占、采空区塌陷[3]。结合实地调查并辅以高分辨率与高光谱遥感数据解译[25],提取挖损、压占、塌陷等信息,如图2所示。
图2 土地损毁遥感影像特征解译
Fig.2 Interpretation of land damage characteristics using remote sensing image
(2)井工矿万吨塌陷率与露天矿万吨损毁率
万吨塌陷率是矿区塌陷土地面积与累计煤炭采出量的比值,可采用式(1)[26]计算:
(1)
式中,S为塌陷区面积,hm2;Q为煤炭开采量,万t;a和b分别为采区倾向和走向在平面上的投影长度,m;δ,β,τ分别为走向、倾向下山、倾向上山地表边界角;dδ,dβ,dτ则对应为主断面开采边界至塌陷边界的水平距离,m;α为煤层的倾角;M为开采厚度,m;γ为煤的密度,t/m3;c为采区采出率,%;P为万吨塌陷率,hm2/万t。
露天矿土地损毁程度既可用剥采比来描述,也可用万吨损毁率表示。由于剥采比不够直观,且与井工开采万吨塌陷率无法比较,已知剥采比后可采取换算的方法得到露天矿万吨损毁率,见式(2):
(2)
式中,η为剥采比,m3/t;Ve为煤炭开采剥离下来的表土体积,m3;Qm为开采的煤炭质量,t;ωe和ωo分别为露天煤矿的万吨挖损率(hm2/万t)与万吨压占率(hm2/万t);ω为露天煤矿的土地损毁率,hm2/万t;hm为煤层厚度,m;he为表土平均厚度即煤层埋深,m;h为煤层厚度与表土厚度之和,m;Sm为煤田底部投影面积;Se为开采土体投影面积,m2;So为排土场投影面积,m2;ho为排土场平均高度,m;ε为土体膨胀系数,无量纲,取ε=1.1;γ=1.5 t/m3。
(1)植被退化信息提取
植被退化信息可通过植被覆盖度变化反映。笔者运用像元二分法模型估算植被覆盖度(FVC):
(3)
其中,NDVIveg,NDVIsoil分别为植被和裸土的NDVI值。利用3期Landsat 8 OLI影像测算矿区范围内植被覆盖度。将植被覆盖度划分为裸土地(0~10%)、低覆盖度(10%~30%),中覆盖度(30%~45%)、和高覆盖度(45%~100%)4个等级。
(2)土壤侵蚀信息提取方法
考虑模型成熟性及资料的可获取性,笔者采用RUSLE模型[27]测算土壤侵蚀:
A=RKLSCP
(4)
其中,A为预测土壤侵蚀量,t/(hm2·a);R为降雨侵蚀力因子,采用国家气象信息中心的基于日降雨资料的降雨侵蚀力估算方法;K为土壤可蚀性因子,基于世界土壤数据库(HWSD)土壤数据集(v1.1)运算得出;L与S分别为坡长与坡度因子(无量纲),依据DEM高程数据提取;C为植被覆盖因子(无量纲);P为水土保持措施因子,将裸地、草地、建设用地、水域分别赋值为0.6,0.1,0.5,0。结合SL190—2007《土壤侵蚀分类分级标准》,将强度、极强度和剧烈侵蚀合并为强度侵蚀,测算结果划分为微度、轻度、中度、强度侵蚀4个等级。
(3)沙漠化/石漠化信息提取方法
沙漠化是土地退化的重要表现。构建反照率(Albedo)-NDVI特征空间[28],见式(5):
I=aNDVI-Albedo
(5)
根据Albedo-NDVI回归方程求解a。利用自然断裂法将I值分为4级,依次为无沙漠化、轻度沙漠化、中度沙漠化、重度沙漠化。
石漠化是由多种营力导致地表基岩裸露的土地退化过程,多发生在湿润半湿润气候区。基于Landsat 8 OLI影像反演NDVI、归一化退化指数(NDDI)、归一化湿度指数(NDMI)和地表温度(Ts)4项因子,运用主成分分析构建石漠化评价模型并测算[29],并将结果划分为无石漠化、轻度石漠化、中度石漠化、重度石漠化共4个等级。
采用灰色关联分析法,分析不同煤矿区土地退化的主要表现[4]。因神东、准格尔2个煤矿区有较长时间、较高分辨率的跟踪监测数据,笔者还运用随机森林法分析了土地退化因素及其主要表现。以NDVI值为参考指标,运用ArcGIS 10.6的PCA计算重分类后的植被退化、沙漠化、土壤侵蚀权重,并对3者加权求和,可得矿区土地退化结果,见式(6):
LD=a1B1+a2B2+a3B3
(6)
其中,LD为土地退化的测算结果;a1,a2,a3分别为植被退化、沙漠化/石漠化、土壤侵蚀的权重;B1,B2,B3对应其分级测算值。LD取值为1~4。将土地退化测算结果划分为无退化(1.00~1.50)、轻度退化(1.50~2.00)、中度退化(2.00~2.75)、重度退化(2.75~4.00)4个等级。
结合调研及测算结果,得到土地损毁面积、万吨损毁率与万吨塌陷率等相关参数,见表1,2。
表1为井工煤矿万吨塌陷率相关参数结果。数据显示:旬耀矿区万吨塌陷率较高,盘县矿区万吨塌陷率最低。采深较大、倾角较小、煤层厚度较小是旬耀矿区万吨塌陷率较高的原因;而盘县矿区是由于煤层厚度大、采深相对较小以及倾角较大所致,故万吨塌陷率较低。结合式(1),影响万吨塌陷率的主要因素有采厚(多煤层开采时,累计采厚)、埋深、倾角、上覆岩性与基岩和表土厚度等,岩性、基岩和表土厚度决定了边界角的大小,从而决定了塌陷影响的范围。
表2为露天矿万吨损毁率测算结果。小龙潭矿区万吨损毁率最低,其特点是煤层埋藏浅、煤层较厚、排土场堆垫高度较高致使剥采比小、损毁率低;石嘴山矿区万吨损毁率最高,这是由于汝箕沟煤矿煤层埋藏深、煤层相对较薄、排土场堆垫高度小,导致剥采比高、损毁率高。因此,结合式(2)得出:露天煤矿区的万吨土地损毁率主要与煤层厚度、剥离岩土厚度及岩土堆放方式有关。
依据信息识别与测算结果,按土地损毁面积加权计算,得出西部重点煤矿区露天煤矿万吨损毁率与井工煤矿万吨塌陷率,结果见表3。
表1 井工煤矿万吨塌陷率及相关参数
Table 1 Parameters and results of the collapse rate of per 10 000 tons for subterranean coal mine
矿区及矿山损毁面积/hm2上覆基岩厚度/m平均倾角/(°)平均埋深/m煤层(累计)开采厚度/m可采煤层数量万吨塌陷率/(hm2·万t-1)榆神矿区常乐堡14349.50.51176.6190.300 0 灵武矿区羊场湾51291.810.01007.7040.257 3 旬耀矿区青岗坪234174.03.03358.2040.324 0 华亭矿区大柳986360.015.075016.4450.246 0 神东矿区补连塔2 738123.02.02054.4270.244 7 盘县矿区红果154230.040.029021.5890.238 0
表2 露天煤矿万吨损毁率及相关测算结果
Table 2 Parameters and results of the damage rate of per 10 000 tons for open-pit coal mine
矿区及矿山损毁面积/hm2剥采比/(m3·t-1)煤层厚度/m煤层埋深/m煤层厚度与埋深之和/m排土场平均高度/m万吨挖损率/(hm2·万t-1)万吨压占率/(hm2·万t-1)万吨损毁率/(hm2·万t-1)石嘴山矿区汝箕沟9976.9026.53137.54164.07600.138 7 0.379 3 0.518 0 准格尔矿区黑岱沟4243.4628.8056.1484.94800.146 0 0.142 7 0.288 7 小龙潭矿区小龙潭2261.2972.0039.72111.722000.052 7 0.021 3 0.074 0 平庄矿区元宝山1283.0984.29113.21197.502140.060 0 0.047 3 0.107 3 胜利矿区胜利1号1 4572.6633.4262.7896.20600.104 0 0.146 0 0.250 0 准东矿区五彩湾2 9233.1360.1574.28134.43800.084 7 0.129 3 0.214 0
表3 西部重点煤矿区土地万吨损毁率和万吨塌陷率
Table 3 Situation of land destruction of 12 coal mine areas in western China hm2/万t
参数黄土高原区西南山地丘陵区北方草原区西北干旱区西部煤矿区平均露天煤矿万吨损毁率0.487 90.074 00.235 20.214 00.304 9井工煤矿万吨塌陷率0.257 60.238 0——0.257 5
(1)植被退化分析
植被退化的解译结果如图3(a)与图4所示。分析发现:① 西北干旱区裸土地与低覆盖度占比超过90%,且无高覆盖度区域,表明植被生长状况较差;西南山地丘陵区高植被覆盖面积约占2/3,植被生长状况良好;黄土高原区与北方草原区则以中、低植被覆盖为主。② 从变化状况看,研究期内:西南山地丘陵区高覆盖度减少3.32%,其他等级稍有增加,表明该地区植被退化较为不明显;黄土高原区中覆盖度面积增加11.74%,其余等级面积有波动,主要由高覆盖度向中覆盖度转移,存在植被退化现象;北方草原区高、中覆盖度分别减少1.60%和3.01%,其他等级面积有所增加,植被退化较为明显;西北干旱区裸土地呈大幅度增加趋势,中覆盖度减少7.87%,表明该区域植被退化状况极为严峻。
(2)土壤侵蚀分析
土壤侵蚀解译结果如图3(b)与图5所示,分析发现:① 黄土高原区、北方草原区、西北干旱区的微度和轻度侵蚀面积占比较高,西南山地丘陵区则以强度侵蚀为主。② 从变化状况看,黄土高原区研究期间微度侵蚀减少10.20%、轻度及中度侵蚀共增加10.34%、强度侵蚀小范围下降,表明侵蚀状况有加重的趋势;西南山地丘陵区3个年度的侵蚀状况有波动,但是2017年较2013年强度侵蚀范围增加4.81%,其他等级侵蚀范围有不同程度的下降,说明土壤侵蚀呈加重的趋势;北方草原区3个年度的侵蚀状况存在波动,2013~2017年间微度侵蚀占比降低2.74%,其他侵蚀等级均有所增加,强度侵蚀主要集中于采矿扰动区周围,侵蚀程度呈加重趋势;西北干旱区微度侵蚀大幅减少9.22%,强度、中度侵蚀面积增加且均分布于采矿扰动区周围,表明区域受采煤扰动导致土壤侵蚀加剧。
图3 西部煤矿区不同植被覆盖度、土壤侵蚀、沙漠化/石漠化程度比较
Fig.3 Comparison of the proportion of different vegetation coverage,soil erosion,and desertification /stony desertification in western CMAs
图4 西部煤矿区植被覆盖分布((a)~(l)同图1)
Fig.4 Distribution of vegetation coverage of 12 coal mine areas((a)~(l) are the same as Fig.1)
图5 西部煤矿区土壤侵蚀分布((a)~(l)同图1)
Fig.5 Distribution of soil erosion of 12 coal mine areas((a)~(l) are the same as Fig.1)
图6 西部煤矿区沙漠化、石漠化分布((a)~(l)同图1)
Fig.6 Distribution of desertification and stony desertification of 12 coal mine areas((a)~(l) are the same as Fig.1)
(3)沙漠化/石漠化分析
沙漠化/石漠化的解译结果如图3(c)与图6所示。分析表明:① 西南山地丘陵区绝大多数土地的石漠化程度较低,西北干旱区主要是中度沙漠化,其他2个区域则以无沙漠化及轻度沙漠化为主。② 从变化状况看,黄土高原区沙漠化存在波动,2013~2017年轻度沙漠化面积增加5.79%,其他沙漠化程度面积小幅减少,中度和重度沙漠化区域向轻度沙漠化转移,表明沙漠化稍有缓解,这可能与黄土高原区近年来的气候变化相关;西南山地丘陵区石漠化程度有微弱的波动,2013,2015,2017年无石漠化面积占比分别为74.41%,83.42%,71.71%,其他程度的石漠化面积稍有增加,石漠化程度基本稳定;北方草原区在研究期内无沙漠化面积减少16.26%,其他程度沙漠化面积有所增加,沙漠化呈加重的趋势;西北干旱区3个年度的中度沙漠化面积变化存在波动,但研究期内中度沙漠化面积减少10.22%,重度沙漠化面积增幅达5.47%,沙漠化呈加剧的态势。
研究区土地的主要损毁形式及退化表现的结果见表4。
(1)黄土高原自然环境较为特殊,易受人为活动影响而加剧土壤侵蚀。石嘴山和准格尔土地损毁的关键因素仍是土地挖损,其次是排土场压占土地;其他井工煤矿的关键因素则为煤矸石堆积压占与采空区塌陷,因此要重视煤矸石的处理和利用,如合理堆放、填埋或充填采空区以及资源化利用,这样既可减轻矸石无序堆放造成的土地退化,又可延长产业链带来一定的经济效益。
(2)西南山地丘陵区的小龙潭煤矿煤层较厚、采坑较深,但剥离岩土堆积压占土地较多,因此小龙潭露天开采导致土地损毁的主要因素为排土场压占;盘县矿区煤矸石压占对土地退化影响较大。该地区山地较多,强降雨可能导致水土流失加重,要关注土壤侵蚀与石漠化演变,预防地质灾害。
(3)北方草原区是典型的草原生态系统,胜利与平庄矿区土地损毁主要为挖损、排土场压占,采矿剥离活动对植被的清除加速了草地逆向演替,区域内水土保持能力下降,有轻度沙漠化及植被退化的倾向。
(4)准东矿区地处戈壁地区,是典型的荒漠生态系统,植被覆盖率极低,自我修复能力较弱,极易受外界干扰发生逆向演替,挖损、压占均对土地退化影响较大,土地退化主要表现为沙漠化和植被退化。
矿区包括多个煤矿,是由煤炭资源正在开采或规划开采的边界圈定的范围,涵盖了煤炭开采、加工等生产及人类生活的聚集区;而矿山则是由单一煤矿开采及其受到影响的范围构成的区域。露天开采与井工开采方式对地表扰动的过程存在差异,故在矿山尺度分析需分别讨论。露天矿山选择矿区范围内的全部露天矿分析,包括露天采坑、排土场、建设用地等,井工矿采动选取1个典型煤矿分析(表4)。土地退化综合评价结果见表5~7及图7。
表4 西部重点煤矿区土地的主要损毁形式及主要退化表现
Table 4 Land damage and external manifestations of land degradation in 12 coal mine areas
区域矿区典型矿山开采方式主要损毁形式主要退化表现石嘴山矿区汝箕沟露天矿露天挖损、压占植被退化、土壤侵蚀灵武矿区羊场湾矿井井工压占、塌陷土壤侵蚀榆神矿区常乐堡矿井井工压占、塌陷沙漠化黄土高原区旬耀矿区青岗坪矿井井工塌陷、压占土壤侵蚀华亭矿区大柳矿井井工压占土壤侵蚀神东矿区补连塔矿井井工压占、塌陷土壤侵蚀准格尔矿区黑岱沟露天矿露天挖损、压占植被退化、土壤侵蚀西南山地丘陵区盘县矿区红果矿井井工压占土壤侵蚀小龙潭矿区小龙潭露天矿露天压占土壤侵蚀北方草原区平庄矿区元宝山西露天矿露天挖损土壤侵蚀、沙漠化胜利矿区胜利露天矿露天挖损、压占植被退化、沙漠化西北干旱区准东矿区五彩湾露天矿露天挖损、压占植被退化、沙漠化、土壤侵蚀
表5 矿区尺度的土地退化数据
Table 5 Statistical of land degradation on the CMA scale
强度分级黄土高原2013年面积比例2015年面积比例2017年面积比例西南山地丘陵区2013年面积比例2015年面积比例2017年面积比例无退化2 349.1523.722 838.5428.662 085.4721.05506.2412.28480.3211.65468.0311.35轻度退化4 874.3749.213 951.5939.904 453.5844.962 785.6867.552 923.5670.892 671.9464.79中度退化2 211.2722.322 460.6824.842 700.8627.27797.7019.34646.9015.69893.2621.66重度退化470.154.75654.136.60665.036.7134.410.8373.251.7890.802.20强度分级北方草原区2013年面积比例2015年面积比例2017年面积比例西北干旱区2013年面积比例2015年面积比例2017年面积比例无退化66.544.8761.384.4955.234.04000000轻度退化288.2121.08266.7919.51270.1119.7642.2611.9537.5410.6217.454.93中度退化560.5941.01481.0235.18634.6046.42242.4468.56237.4567.15240.2567.94重度退化451.7833.05557.9340.81407.1829.7868.9219.4978.6322.2495.9227.13
注:面积单位是km2,比例单位是%,下同。
表6 露天矿山尺度的土地退化数据
Table 6 Statistical of land degradation on the open-pit coal mine scale
强度分级黄土高原2013年面积比例2015年面积比例2017年面积比例西南山地丘陵区2013年面积比例2015年面积比例2017年面积比例无退化12.4823.8211.1021.189.7218.552.8520.852.2716.612.5718.80轻度退化9.1217.409.7018.514.468.513.0922.604.0129.333.0021.95中度退化21.6841.3721.3240.6921.9141.815.2238.194.0229.415.4239.65重度退化9.1217.4010.2819.6216.3131.132.5118.363.3724.652.6819.60强度分级北方草原区2013年面积比例2015年面积比例2017年面积比例西北干旱区2013年面积比例2015年面积比例2017年面积比例无退化3.072.402.952.311.661.30000000轻度退化7.285.705.824.558.506.653.478.022.245.182.164.99中度退化43.0533.6842.3733.1532.4825.4123.6954.7522.3451.6321.3449.32重度退化74.4158.2276.6759.9985.1766.6416.1137.2318.6943.1919.7745.69
表7 井工矿山尺度的土地退化数据
Table 7 Statistical of land degradation on the underground coal mine scale
强度分级黄土高原2013年面积比例2015年面积比例2017年面积比例西南山地丘陵区2013年面积比例2015年面积比例2017年面积比例无退化66.4929.4862.7727.8352.1723.130.2522.940.2220.180.2220.18轻度退化75.0333.2775.2833.3875.7333.580.4339.450.3935.780.4137.61中度退化68.9430.5770.2431.1472.7132.240.2623.850.2926.610.2623.85重度退化15.086.6917.257.6524.9311.050.1513.760.1917.430.2018.35
(1)黄土高原区。① 矿区尺度以轻度退化为主,露天矿山尺度则以中度、重度退化为主,露天矿山重度退化比例约17.40%~31.13%,井工矿山尺度以轻度、中度退化为主,重度退化比例远小于露天矿山,约6.69%~11.05%,表明露天采煤扰动比井工采煤更为剧烈。② 研究期内,矿区无退化和轻度退化面积占比减少2.67%,中度及重度退化面积共增加6.91%;露天矿山重度退化面积增加13.73%,无退化面积减少5.27%;井工矿山重度退化面积增加4.36%,无退化面积减少6.35%。由此说明矿区和矿山尺度土地退化均呈加剧趋势。
(2)西南山地丘陵区。① 矿区尺度以轻度退化为主,露天矿山以中度退化为主,井工矿山以轻度退化为主。② 研究期内,2015年在矿区和矿山尺度都表现出波动,可能与当年的气候有关。2013—2017年矿区及井工矿山的各个退化强度上变化不大,土地退化趋势不明显。小龙潭露天矿所处特殊的地理位置及煤炭赋存特点,采坑面积基本稳定、排土场面积小幅增加,土地退化稍有增加。
(3)北方草原区。① 矿区尺度以中度退化为主,矿山尺度则以重度退化为主,表明草原矿区煤炭开采扰动剧烈,露天矿山土地退化比矿区尺度更为明显。大量剥离物堆积成的排土场土壤瘠薄,加速了土地沙化,使得矿山土地退化态势严峻。② 研究期内矿区与矿山尺度的土地退化变化趋势并不一致,矿区重度退化面积减少3.27%,中度退化面积增加5.41%,退化程度稍有减缓;矿山重度退化面积增加8.42%,中度退化面积减少8.27%,土地退化加重。
(4)西北干旱区。① 矿区中度退化约占矿区总面积2/3,矿山以中度、重度退化为主,2者之和超过矿山总面积的90%,表明该区域土地重度退化主要集中在矿山及其周边。② 研究期内,矿区重度退化面积增加7.64%,矿山重度退化面积增加4.59%,矿山与矿区尺度土地退化均表现为加重的趋势。值得关注的是,相比其他3个区域,该区域不存在无退化区域,土地退化形势更为严峻。这是由于该区域位于典型大陆性温带极干旱气候区,多年平均降水量仅35 mm,而蒸发量却达2 800 mm,水资源极度匮乏,戈壁地貌有机质积累量少,植被覆盖度较低,生态环境极为脆弱,易损毁、难修复,因此新疆等地区煤炭资源的开发强度必须适应其环境承载力。
如果用煤矿扰动范围与矿区总面积之比表示煤炭资源开发强度,那么露天开采引起矿区土地退化的程度就与煤炭资源开发强度有关,见表8。准格尔露天矿区包括黑岱沟和哈尔乌素2个矿山群,煤矿扰动范围与矿区总面积之比高达35.54%,矿山尺度与矿区尺度的土地退化变化基本处于同等水平,以重度退化面积增加为主。平庄矿区煤矿扰动范围与矿区总面积之比相对较小,仅为7.40%,虽矿山尺度重度退化面积有所增加,但矿区尺度的重度退化面积减少。准东矿区处于戈壁的脆弱生态环境,煤矿扰动范围与矿区总面积之比为12.24%,虽然不高,但矿区尺度与矿山尺度重度退化面积均呈增加趋势。因此,为防止矿区土地过度退化,根据不同区域的生态承载力,露天煤矿存在合理的开发强度,以约束采煤影响的空间扩散。
图7 西部12个重点煤矿区土地退化分级((a)~(l)同图1)
Fig.7 Distribution of land degradation of 12 coal mine areas((a)~(l) are the same as Fig.1)
表8 2013—2017年露天采矿在矿山与矿区尺度的土地退化数据
Table 8 Changes of land degradation on the open-pit coal mine scale & CMA scale from 2013 to 2017
矿区名称开发强度/%尺度面积变化/km2无退化轻度退化中度退化重度退化矿区名称开发强度/%尺度面积变化/km2无退化轻度退化中度退化重度退化石嘴山14.54矿山-0.14-0.89-0.461.49平庄7.40矿山0-1.82-6.248.06矿区-12.84-28.3738.103.11矿区0.28-12.1861.28-49.38准格尔35.54矿山-2.62-3.770.695.70胜利13.76矿山-1.413.04-4.332.70矿区-9.12-0.821.128.82矿区-11.59-5.9212.734.78小龙潭20.21矿山-0.28-0.090.200.17准东12.24矿山0-1.31-2.353.66矿区-16.977.3810.11-0.52矿区0-11.357.673.68
(1)揭示了西部煤矿区万吨塌陷率与万吨损毁率的主要影响因素。万吨塌陷率:除旬耀矿区较高外,其余矿区的万吨塌陷率差异较小,其主要影响因素有累计采厚、平均采深、煤层倾角、上覆岩性及基岩和表土厚度。万吨损毁率:西南山地丘陵区的小龙潭矿区万吨损毁率最低(0.518 0 hm2/万t),黄土高原区的石嘴山矿区万吨损毁率最高(0.074 0 hm2/万t),其主要影响因素包括煤层厚度、剥离岩土厚度及岩土堆放方式。
(2)辨识了西部煤矿区土地损毁及其退化的主要外在表现。导致矿区土地损毁的主要因素包括:露天矿山的挖损与排土场压占,井工矿山的煤矸石压占与采空区塌陷。不同区域的退化表现也存在差异,黄土高原区采矿扰动加剧土地退化,主要表现为土壤侵蚀和植被退化;西南山地丘陵区土地退化主要表现为土壤侵蚀;北方草原区土与西北干旱区环境本底条件较差,土地退化的主要表现兼有土壤侵蚀、植被退化与沙漠化。
(3)分析了矿区与矿山两个尺度土地退化程度、趋势及其差异性。矿区尺度上,黄土高原区与西南山地丘陵区轻度退化面积占比较高,北方草原区与西北干旱区则以中度退化为主。露天矿山尺度:西南山地丘陵区以中度退化为主,西北干旱区、黄土高原区以中度、重度退化为主,北方草原区以重度退化为主。井工矿山尺度:黄土高原区主要为轻度及中度退化,西南山地丘陵区主要为轻度退化。整体看,矿山尺度的土地退化比矿区尺度更为剧烈,表明煤炭开采是导致矿山土地退化的关键因素,且露天开采对土地退化的扰动强度远强于井工开采。因此,针对不同区域、不同采煤方法造成的土地退化需要采取不同的生态修复策略。
(4)发现矿区尺度的土地退化程度与煤炭资源开发强度有关,开发强度越大,矿区土地退化程度越严重,而矿区土地退化程度与开发强度并不呈正比,它还与矿区的生态脆弱性或承载力有关。为防止过度土地退化,需要根据矿区生态承载力确定合适的开发强度。因此,需要根据矿区生态的脆弱性和承载力的差异性,将煤炭开采规模与开发强度控制在合理区间,以约束采煤影响的空间扩散。
[1] 江泽慧.中国干旱地区土地退化综合监测与评价指标体系建设研究[M].北京:中国林业出版社,2013.
[2] 郭晓娜,陈睿山,李强,等.土地退化过程、机制与影响——以土地退化与恢复专题评估报告为基础[J].生态学报,2019,39(17):6567-6575.
GUO Xiaona,CHEN Ruishan,LI Qiang,et al.Processes,mechanisms,and impacts of land degradation in the IPBES Thematic Assessment[J].Acta Ecologica Sinica,2019,39(17):6567-6575.
[3] 卞正富,雷少刚,金丹,等.矿区土地修复的几个基本问题[J].煤炭学报,2018,43(1):190-197.
BIAN Zhengfu,LEI Shaogang,JIN Dan,et al.Several basic scientific issues related to mined land remediation[J].Journal of China Society,2018,43(1):190-197.
[4] 于昊辰,牟守国,卞正富,等.北方草原露天煤矿区植被退化因素分析[J].生态与农村环境学报,2019,35(1):1-8.
YU Haochen,MU Shouguo,BIAN Zhengfu,et al.Analysis for vegetation degradation factors in opencast coal mines in Northern Grass-Land Area,China[J].Journal of Ecology and Rural Environment,2019,35(1):1-8.
[5] 白中科,周伟,王金满,等.再论矿区生态系统恢复重建[J].中国土地科学,2018,32(11):1-9.
BAI Zhongke,ZHOU Wei,WANG Jinman,et al.Rethink on ecosystem restoration and rehabilitation of mining areas[J].China Land Science,2018,32(11):1-9.
[6] 卞正富,沈渭寿.西部重点矿区土地退化因素调查[J].生态与农村环境学报,2016,32(2):173-177.
BIAN Zhengfu,SHEN Weishou.Driving factors of land degradation in major mining areas in West China[J].Journal of Ecology and Rural Environment,2016,32(2):173-177.
[7] 王丽,雷少刚,卞正富.系统视角下中国西部煤炭开采生态损伤与自然修复研究综述[J].资源开发与市场,2017,33(10):1188-1192.
WANG Li,LEI Shaogang,BIAN Zhengfu.Review on study of ecological damage and natural recovery in the coal mining subsidence area in Western China[J].Resource Development & Market,2017,33(10):1188-1192.
[8] 李海东,沈渭寿,司万童,等.中国矿区土地退化因素调查:概念、类型与方法[J].生态与农村环境学报,2015,31(4):445-451.
LI Haidong,SHEN Weishou,SI Wantong,et al.Investigation of driving factors of land degradation in mine areas in China:Concept,types and approaches[J].Journal of Ecology and Rural Environment,2015,31(4):445-451.
[9] AHIRWAL Jitendra,MAITI Subodh Kumar,SINGH Ashok Kumar.Changes in ecosystem carbon pool and soil CO flux following post-mine reclamation in dry tropical environment,India[J].Science of the Total Environment,2017(583):153-162.
[10] CHUAI Xiaowei,QI Xinxian,ZHANG Xiuying,et al.Land degradation monitoring using terrestrial ecosystem carbon sinks/sources and their response to climate change in China[J].Land Degradation & Development,2018,29(10):3489-3502.
[11] OUEDRAOGO Issa,BARRON Jennie,TUMBO Siza,et al.Land co-ver transition in northern Tanzania[J].Land Degradation and Development,2015,27(3):682-692.
[12] NOVIANTI Vivi,MARRS Rob H,CHOESIN Devi N,et al.Natural regeneration on land degraded by coal mining in a tropical climate:Lessons for ecological restoration from Indonesia[J].Land Degradation & Development,2018,29(11):4050-4060.
[13] MAO Dehua,WANG Zongming,WU Bingfang,et al.Land degradation and restoration in the arid and semiarid zones of China:Quantified evidence and implications from satellites[J].Land Degradation & Development,2018,29(11):3841-3851.
[14] IBEZ J,CONTADOR J F,LAVADO Schnabel S,et al.A model-based integrated assessment of land degradation by water erosion in a valuable Spanish rangeland[J].Environmental Modelling & Software,2014,55:201-213.
[15] ABDEL-KADER Fawzy Hassan.Assessment and monitoring of land degradation in the northwest coast region,Egypt using Earth observations data[J].The Egyptian Journal of Remote Sensing and Space Science,2019,22(2):165-173.
[16] 孙武,李森.土地退化评价与监测技术路线的研究[J].地理科学,2000(1):92-96.
SUN Wu,LI Sen.Technical framework on monitoring and assessing of land degradation[J].Scientia Geographica Sinica,2000(1):92-96.
[17] 沈渭寿,曹学章,沈发云.中国土地退化的分类与分级[J].生态与农村环境学报,2006(4):88-93.
SHEN Weishou,CAO Xuezhang,SHEN Fayun.Classification and grading of land degradation in China[J].Journal of Ecology and Rural Environment,2006(4):88-93.
[18] 杨超,王金亮.流域土地退化研究综述[J].云南师范大学学报(自然科学版),2015,35(1):73-78.
YANG Chao,WANG Jinliang.Progress in the research on basin land degradation[J].Journal of Yunnan Normal Unviersity,2015,35(1):73-78.
[19] WU Zhenhua,LEI Shaogang,HE BaoJie,et al.Assessment of landscape ecological health:A case study of a mining city in a semi-arid steppe[J].International Journal of Environmental Research and Public Health,2019,16(5):752.
[20] 郭山川,汤傲,李效顺,等.融合主被动遥感的乌海矿区土地损伤测度[J].生态与农村环境学报,2018(8):678-685.
GUO Shanchuan,TANG Ao,LI Xiaoshun,et al.Measurement of land damage based on active and passive remote sensing in Wuhai Mining Area[J].Journal of Ecology and Rural Environment,2018(8):678-685.
[21] 李心慧,雷少刚,田雨,等.基于时间序列的准格尔典型矿区土地退化特征分析[J].河南农业大学学报,2019,53(2):273-281.
LI Xinhui,LEI Shaogang,TIAN Yu,et al.Analysis of land degradation characteristics in typical mining areas of Jungar Banner based on time sequence[J].Journal of Henan Agricultural University,2019,53(2):273-281.
[22] 陈彬,俞炜炜,陈光程,等.滨海湿地生态修复若干问题探讨[J].应用海洋学学报,2019,38(4):464-473.
CHEN Bin,YU Weiwei,CHENG Guangcheng,et al.Coastal wetland restoration:An overview[J].Journal of Applied Oceanography,2019,38(4):464-473.
[23] 杨永均,张绍良,侯湖平,等.基于非线性动力学模型的矿山土地生态系统恢复力机制[J].煤炭学报,2019,44(10):3174-3184.
YANG Yongjun,ZHANG Shaoliang,HOU Huping,et al.Resilience mechanism of land ecosystem in mining area based on nonlinear dynamic model[J].Journal of China Coal Society,2019,44(10):3174-3184.
[24] 卞正富.矿山生态学导论[M].北京:煤炭工业出版社,2015.
[25] 姚维岭,余江宽,路云阁.基于ZY-3卫星数据的神东煤矿区土地退化人为影响因素调查与评价[J].生态与农村环境学报,2016,32(3):355-360.
YAO Weiling,YU Jiangkuan,LU Yunge.Investigation and assessment of artificial influencing factors of land degradation in Shendong Coal Mining Area Based on ZY-3 satellite data[J].Journal of Ecology and Rural Environment,2016,32(3):355-360.
[26] 尹国勋,邓寅生,李栋臣,等.煤矿环境地质灾害与防治[M].北京:煤炭工业出版社,1997.
[27] 李天宏,郑丽娜.基于RUSLE模型的延河流域2001—2010年土壤侵蚀动态变化[J].自然资源学报,2012,27(7):1164-1175.
LI Tianhong,ZHENG Lina.Soil erosion changes in the Yanhe Watershed from 2001 to 2010 based on RUSLE model[J].Journal of Natural Resource,2012,27(7):1164-1175.
[28] MA Zongyi,XIE Yaowen,JIAO Jizong,et al.The construction and application of an ALEDO-NDVI based desertification monitoring model[J].Procedia Environmental Sciences,2011,10:2029-2035.
[29] 李辉霞,李森,周红艺.粤北典型岩溶山区土地石漠化程度遥感评价[J].水土保持通报,2011,31(6):155-159.
LI Huixia,LI Sen,ZHOU Hongyi.Assessment of rocky desertification in typical Karst Mountain Area of North Guangdong Province based on RS[J].Bulletin of Soil and Water Conservation,2011,31(6):155-159.