随着经济发展与人口、资源、环境之间的矛盾日益加剧，人类改变土地资源的程度不断增大，导致土地利用形式及地表覆被状况发生较大变化，进而改变了自然地理环境状况^[1-2]。诸多研究表明，土地利用/土地覆被变化是影响全球及区域环境变化的主要因素^[3-4]。近年来，土地利用/土地覆被变化在城市、绿洲等区域的研究成果较多，而在大型矿业基地的研究较少，国内外学者对大型矿业基地的研究集中在景观格局变化^[5-6]、土地利用变化、环境影响评价^[7-10]等方面。因此对大型煤炭基地土地利用/土地覆被变化与人文驱动因素的耦合关系进行研究，掌握土地利用/土地覆盖在不同时空尺度上的变化趋势，揭示土地利用变化的驱动机制及其对矿区经济与环境所产生的影响，提出针对性的防控政策，对矿产资源持续开发及矿区环境保护具有重要的现实意义。

矿产资源开发活动会对矿山土地利用类型及植被类型产生重要影响，资源毁损、地质灾害和环境污染直接或间接破坏土地和植被类型^[11]。作为国家规划建设的14处亿吨级煤炭资源开采基地之一的宁东煤炭基地，面积3484km²，开发始于2003年，目前处于大规模开发建设期间，煤炭资源开采工业建筑物建设占压土地、采空塌陷破坏土地、煤矸石堆积压占林草地等，土地与煤炭资源开采之间已演化为一种互耗的关系，影响宁东煤炭基地的可持续发展^[12-13]。在宁东地区土地利用变化研究方面，马斌等^[14]认为，宁东煤炭基地2000—2007年土地利用结构发生了较大变化，林地和天然草地面积增加，耕地、裸土地和沙地面积减少；闫军^[15]等利用遥感技术，研究土地利用特征、地表生态环境、矿山地质环境等问题，提出基地矿区土地利用与生态恢复治理对策。前人对宁东煤炭基地的相关研究主要集中在2010年以前，并且采用的监测底图分辨率低，信息提取精度不高，没有开展LUCC驱动力研究。

煤炭资源大规模开发建设会使基地范围内土地利用及植被类型发生何种变化，这种变化会对宁东生态脆弱区产生何种影响，已成为政府及科学家关注的重大社会环境问题。本文主要采用高分辨率的卫星影像数据，开展宁东煤炭基地区域土地利用结构的演变过程及驱动机制的研究，旨在为政府规划建设及矿区生态环境恢复治理提供依据。

1.   研究区概况

宁东煤炭基地位于宁夏中东部，地理坐标北纬37°07'15″~38°17'17″、东经106°18'49″~106°57'44″，规划面积3484km²，涉及灵武市、盐池县、红寺堡开发区和同心县（图 1）。基地属中温带干旱半荒漠气候区，干燥少雨、蒸发量大、水资源缺乏，土壤类型以淡灰钙土和风沙土为主，区域生物量低、植被恢复困难，植被以荒漠草原和草原带沙生植被为主。基地主要包括红石湾煤矿、羊场湾煤矿、梅花井煤矿、清水营煤矿、灵新煤矿、石沟驿煤矿、金凤煤矿、枣泉煤矿8个大中型煤矿。已探明煤炭资源地质储量273.14 × 10⁸t，探明储量占全宁夏探明储量的88.6%，含煤面积3500km²，远景储量达1394.3×10⁸t，是中国规划建设的14个亿吨级煤炭开采基地之一。但基地范围内，大部分地区为荒漠型草地，2000年后煤炭资源处于开发规模不断加大的过程，各矿区均不同程度地出现了裂缝、塌陷等变形现象，煤矸石、废渣侵占土地资源，河道洗砂、洗煤厂等污染水土资源，地表植被退化、水土流失等各类环境问题，持续开发对土地利用类型和植被类型产生了明显的影响。

图  1  研究区地理位置示意图

Figure  1.  Geographical location of the study area

下载: 全尺寸图片 幻灯片

2.   研究方法

2.1   数据来源

研究采用分辨率优于2.5m的法国SPOT-5卫星影像数据，融合后空间分辨率为2.5m，接收日期分别是2005年6月20日和2010年7月11日，波段合成为R4G1B3；2015年选用的是2015年7月28日国产高分二号卫星影像，空间分辨率为1m，波段合成为R4G1B3。利用ERDAS10.4对影像数据进行校正、融合、镶嵌等处理，3期影像数学基础一致、色调均匀，拼接处无明显接缝、云覆盖低于5%。收集2000年土地利用遥感解译数据，分类标准、精度及数学基础与2005年、2010年、2015年3期解译土地利用信息一致，满足研究需求。气象、人口等数据均来自宁夏统计年鉴（2000—2016年）。

2.2   土地利用信息提取

采用人工目视解译结合野外调查验证手段，提取土地利用专题信息。经过野外实地查证，遥感解译正确率87.16%。以ERDAS IMAGINE10.4为平台，提取草地覆盖度，建立NDVI模型，得到草地的NDVI灰度图，利用clump和eliminate功能获取草地覆盖信息，分析草地变化。

2.3   土地利用变化综合指数

借鉴刘纪远^[16]在西藏土地利用研究中从生态学角度提出的一种综合分析方法——土地利用程度综合指数模型，研究土地利用变化。其原理为：将土地自然综合体在社会因素影响下的自然平衡状态分为4级，即未利用级、自身再生利用级、认为再生利用级、非再生利用级，并赋予分级指数1~4，逐级表达土地的利用程度。土地利用程度越高，人类对其干扰程度就越大，可恢复性越差。通过以上含义指数，利用各地类土地利用率、土地利用程度分级指数，评价矿区土地利用类型结构总体变化特征，进而反映自然环境和人类对土地利用结构影响的耦合关系^[17-19]。

式中，L表示土地利用程度综合指数，L∈100~ 400，L越接近400，开发利用程度越高；A_i表示土地利用类型分级指数，按照土地利用程度指数的分级标准，建设用地的土地利用程度为第四级，包括居民、采矿、工业、交通等用地；第三级为耕地等农用地；第二级为林地、牧草地和水域；第一级为沙地和其他地类（表 1）。C_i表示土地利用类型面积百分比；△L_b-a表示土地利用变化综合指数变化量；L_a、L_b表示a、b两个时间段的土地利用程度综合指数；C_ia、C_ib表示a、b两个时间段第i类土地类型面积百分比；R表示土地利用程度变化率。R＞0为发展阶段；R＜0为衰退阶段；R=0则为稳定或调整阶段。

表  1  宁东煤炭基地土地利用程度综合指数及变化

Table  1.  Comprehensive index of land use degree and its change rate in Ningdong coal base

分级类型	分级指数	土地利用类型	2000年	2005年	2010年	2015年
农用地	3	耕地	51.65	49.53	41.78	43.46
林地、草地水域用地	2	林地	0.75	19.64	20.51	30.29
		草地	105.64	128.98	130.67	115.52
		水域	0.83	2.91	3.28	2.74
建设用地	4	居民地	5.93	6.08	6.79	7.74
		采矿用地	0.77	2.33	4.5	7.04
		工业用地	0.54	1.04	4.9	9.91
		交通用地	1.27	2.09	3.53	4.46
未利用地	1	沙地	6.47	4.25	3.54	2.94
		其他地类	20.85	0.59	0.37	1.01
土地利用程度综合指数L			194.7	217.45	219.88	225.12

下载: 导出CSV 

| 显示表格

2.4   土地利用动态度

土地利用类型在单位时间内面积变化情况称土地利用动态度，反映了土地利用变化速度，可定量描述区域土地利用变化的速度^[20-24]。

式中，H为土地利用动态度；T₂-T₁为研究初期和末期时间间隔值；W₂-W₁为研究时段内该类土地面积差值。

2.5   土地利用转换趋势

采用土地利用转移矩阵反映土地利用类型转变趋势的分析方法^[25-26]，以及土地利用类型变化特征和各类型之间的转移情况。运用Arcgis10.1中的空间分析模块，结合Excel数据透视表功能建立转移矩阵，利用2005年、2015年2期遥感影像提取土地利用信息后叠加，得到基地土地利用类型转移特征，确定土地利用类型互相转换的数量和方向。

3.   结果与分析

3.1   土地利用程度

宁东煤炭基地开发开始于2003年，2005年已经完成了部分建设。为了掌握基地开发前后土地利用结构的变化情况，收集2000年宁东煤炭基地的土地利用遥感解译数据作为开发建设前的数据，与2010年、2015年的数据一起进行分析。宁东煤炭基地2000—2015年土地利用结构发生了较大变化（表 1），2000年、2005年、2010年、2015年4期土地利用程度综合指数分别为194.7、217.45、219.88、225.12，在100~400范围内，说明土地利用一直处于合理开发阶段^[13]。2000—2015年土地利用程度综合指数从194.70持续增加至225.12，增加了30.42，土地利用程度变化量大于零，说明土地利用处于发展阶段。第一阶段（2000—2010年）土地利用程度变化量为25.18，第二阶段（2010—2015年）土地利用程度变化量为5.24，表明第二阶段开发利用程度较第一阶段大幅度减小，从土地利用程度综合指数增加幅度总体可知，建设逐渐进入稳定阶段，土地利用压力逐渐下降^[27]。

3.2   土地利用动态度特征

2005—2010年，基地主要土地利用动态度大小为工业用地（18.53%）＞采矿用地（4.65%）＞交通用地（3.45%）＞沙地（-0.83%）＞耕地（-0.78%），即工业、采矿及交通用地增速显著，沙地和耕地缓慢减少。分析其原因为，宁东煤炭基地开发建设初期（第一阶段）, 主要进行工业园区主体建设及基础设施建设，使工业园区（工业用地）、基础设施用地（连接各工矿企业之间的道路及连通其他主要城市的铁路、高速路等交通运输用地）等建设用地呈现增速显著的特征；由于退耕还林措施的实施及基地建设占用部分沙地及耕地，使耕地呈现缓慢减少的特征。

2010—2015年主要土地利用动态度大小为工业用地（5.10%）＞采矿用地（2.82%）＞林地（2.38%）＞交通用地（1.32%）＞沙地（-0.85%），即工业用地、采矿用地及交通用地增速放缓，林地加速扩展，沙地和草地持续低速减少（图 2）。分析其原因为，基地建设第二阶段，化工园区规划建设框架基本完成，在初期建设的基础上进行了一些大中型能源化工业类、煤炭资源开发类项目及基础设施建设，使基地工业用地、采矿用地及建设用地呈现持续增加但增速较缓的特征；加强基地生态环境建设，进一步加大退耕还林等生态措施力度，部分耕地、草地种植了大面积柠条，使林地呈现增速显著的特征；开发建设的第一阶段及第二阶段沙地呈现匀速持续减少的特征，分析其原因主要为，基地开发建设及基地生态建设使得沙地面积逐渐减少。

图  2  宁东煤炭基地主要地类土地利用年均动态图

Figure  2.  The average annual rate of the main land use type in Ningdong coal base

下载: 全尺寸图片 幻灯片

3.3   土地利用转换特征

2005—2015年, 基地土地变化面积大小依次为建设用地（采矿用地、工业用地、交通用地、建制镇）、草地、林地及耕地（图 3）。特征一：以草地、林地及耕地向采矿用地、工业用地、交通用地及建制镇为主的建设用地转换。建设用地流出面积26km²，转变为草地（19.3km²），占建设用地总流出面积的74.23%；建设用地流进面积176.4km²，主要由2005年的草地、耕地、林地及沙地向建设用地转换，其中草地转进112.4km²，占总流进面积的63.72%，其次为林地和耕地，比例分别为21.60%和10.15%。特征二：耕地、草地向林地转换使其显著增加，草地面积大幅度减少是基地土地利用变化的另一个主要特征。林地流进总面积312.2km²，由2005年的草地、耕地转进，比例分别为85.75%和10.08%；草地流出面积为229.2km²，转为林地的面积最大，其次为工业用地和沙地，其他地类较小。

图  3  2005—2015年宁东煤炭基地主要地类变化空间分布

Figure  3.  Spatial distribution of main land use types in Ningdong coal base in 2005-2015

下载: 全尺寸图片 幻灯片

2005—2015年，土地利用转换主要集中在基地北部。其中，北部宁东镇部分草地建设成了综合工业园区及煤化工园区，宁东镇第三居委会西南部草地建设成了临州石油天然气化工园区、庆华煤化工园区；北部杨家庄、赵路圈、中部郭记场和南部王家沟附近耕地退耕为草地；北部方家塘、中部周家沟附近耕地种植柠条变成了林地；北部任家庄与范家庄子草地种植柠条变成了林地；北部宁东镇第二居委会的李家圈、甜水井、中部的马家庄子、马家窑与南部的新庄、孙家山附近的沙地变成了草地。可见，该基地土地利用转换的空间分布特征表明，其主要表现为边开发占用、边补偿生态用地的特征。

3.4   驱动力分析

（1）定量分析

矿区土地利用变化是一个复杂的过程，受多种自然和人为因素影响，其中人类活动是矿区土地利用变化的主要驱动力^[28-32]。本文以矿业开发区域土地利用变化驱动因素的理论体系为基础，基于基地特有的资源依托性及化工产业发展过程中对生态的约束，从人口、经济、煤炭产业、生态约束等方面选取指标，构建宁东基地土地利用变化驱动力定量评价的指标体系，利用主成分分析法定量分析宁东煤炭基地土地利用变化驱动力。运用SPSS22软件得到主成分载荷矩阵（表 2）。

由表 2看出，第一主成分与x5，x7，x10，x12，x13有较大正相关，而这5个变量与基地资源型产业及经济水平有关性，因此第一主成分可看作煤炭资源产业化与经济水平的代表；第二主成分与x1和x2有较大的相关性，而这2个变量与基地人口有关，因此第二主成分可看作基地人口的代表。以上分析结果表明，选取的16个变量指标可由煤炭资源产业化规模与经济水平及人口规模3个主要因子代表。

表  2  主成分载荷矩阵

Table  2.  The component matrix

变量指标	主成分
	第一主成分	第二主成分
总人口(x1)	0.432	0.883*
第二产业从业人员(x2)	0.037	0.947*
第三产业从业人员(x3)	-0.026	0.055
人口增长率(x4)	0.542	0.039
人均生产总值(x5)	0.839*	0.011
工业生产总值(x6)	0.746	0.008
第三产业生产总值(x7)	0.932*	0.005
固定资产投资额(x8)	0.642	0.007
消费性支出(x9)	0.537	0.025
能源生产总量(x10)	0.978*	0.013
能源消耗总量(x11)	0.732	0.014
采矿业就业人数(x12)	0.895*	0.033
煤产值(x13)	0.917*	0.021
工业废水排放量(x14)	0.286	0.423
工业废气排放量(x15)	0.483	0.581
工业固体废弃物产生量(x16)	0.356	0.761
注：*表示与第一、第二主成分相关较大的影响因素

下载: 导出CSV 

| 显示表格

① 产业化因素：宁夏持续以煤炭、电力和煤化工为支撑，建设国家级重点能源化工基地，目前，已建成并投入生产的包括煤矿、焦炉煤气制甲醇、洗煤厂、电厂、电解铝等一大批大中型能源化工业项目。化工园区占用大面积的草地、耕地等，近10年工矿用地面积增加了118.26km²，年变化率达到了40.20%，交通运输用地增加了21.02km²，年变化率达到11.37%，新增的工矿及交通用地主要占用草地、林地、耕地等农用地。

② 人口因素：表 2中与第二主成分有较大正相关性的是第二产业从业人员及总人口，均为人口因素，基地主要涉及6个乡镇，2005—2015年人口数量从7.3万增加到11.6万人，人口增长导致城镇化、交通运输、区位条件、市场需求等变化，并通过对人们生活方式和价值观的影响改变着土地利用的格局^[33]，城镇建设用地面积从15.36km²增加到21.73km²。

（2）定性分析

① 封山禁牧因素：基地从2005年开始实施封山禁牧、退耕还林等措施，十年间耕地面积减少了70.43km²，较2005年减少了12.24%；林地面积大幅增加，较2005年增加了54.29%。政府政策对土地利用格局的变化有强制性的影响^[34]，退耕还林、还草的实施，是基地林地增加的主要原因。

② 自然因素：基地2005—2015年气温呈缓慢起伏上升趋势（图 4），平均升温率为0.029，极端最低气温-28.0℃，极端最高气温41.4℃，明显增温导致地表蒸发和植物蒸腾作用加强，对植物生长造成一定影响^[35]；同时，基地降水总体呈显著增加趋势（图 5），平均降水量为203.27mm，降水量年际间变化较大，最大降水量为2012年的341.1mm，最小降水量为2005年的80.4mm。气温及降水量的变化对基地土地利用格局的变化产生了一定影响。

图  4  宁东地区2005—2015年均温变化

Figure  4.  The average temperature curve of Ningdong in 2005-2015

下载: 全尺寸图片 幻灯片

图  5  宁东地区2005—2015年降水量变化

Figure  5.  The precipitation change curve of Ningdong in 2005-2015

下载: 全尺寸图片 幻灯片

4.   结论

（1）2005—2015年，基地草地、耕地、沙地面积有明显减少，采矿、工业及交通用地面积有较大幅度增加，水域面积总体变化不明显，自然水域面积不断减少而人工水域面积有所增加。2005—2010年，工业、采矿及交通用地增速显著，沙地和耕地缓慢减少；2010—2015年，工业、采矿及交通用地持续中速增长，林地扩展速度较快，沙地和草地减少速度较慢。

（2）宁东煤炭基地土地利用程度综合指数表明，2005—2015年，土地利用一直处于合理开发状态，2010年以后开发利用程度较之前略有加大，但总体稳定，土地利用压力趋于下降。

（3）基地土地类型转移特征为：以草地、林地及耕地向采矿用地、工业用地、交通用地及建制镇为主的建设用地转换。耕地、草地向林地转换，草地面积大幅度减少。

（4）基地土地利用变化主要由能源化工产业化规模化、人口增长和城镇化驱动及封山禁牧政策实施，自然因素对土地利用类型的变化也产生了一定影响。