基于灰色预测模型的昆仑山区土地利用时空变化研究

姜云璐 曹月娥 贺忠发 宗亚飞 何浩天 赵敏 许仲林

摘要 運用2000、2005、2010和2015年LUCC土地利用数据，采用土地利用动态度模型及土地利用转移矩阵，对处于昆仑山区的新疆和田地区16年来的土地利用变化情况进行分析，并运用灰色预测模型对未来2025、2030年和田地区土地利用类型进行预测。结果表明，近16年来和田地区为代表的昆仑山区土地利用处于快速变化时期。土地利用类型变化速度与幅度反映2000—2015年和田地区各类土地利用变化中耕地和建设用地面积增长较明显，分别增加了985、21 km2，其他各类用地面积都有不同程度的减少，减少最多的是草地，减少面积555 km2，其次是未利用地，减少面积398 km2。从土地利用动态度来看，2000—2015年耕地和建设用地单一动态度较大，其他地类单一动态度较小，表明在研究期间耕地和建设用地面积增长较快，其他地类面积都有不同程度减少。从土地利用转移方面来看，土地利用类型转移主要发生在2000—2005和2010—2015年，转移地类主要是发生在草地—耕地、未利用地—耕地和耕地—建设用地。土地利用转移空间变化主要表现在沿河流纵向分布，发生在人类聚居的绿洲区域及城市周边。利用灰色预测模型对2025和2030年和田地区土地利用类型进行预测，预测精度良好，表明和田地区2025和2030年耕地和建设用地面积呈不断增加，其他地类都有不同程度减少。

关键词 昆仑山区;和田地区;灰色预测模型;土地利用;时空变化

中图分类号 F 301.2  文献标识码 A  文章编号 0517-6611（2022）01-0070-06

doi：10.3969/j.issn.0517-6611.2022.01.019

Research on Spatio-temporal Changes of Land Use in Kunlun Mountains Based on GM Model—Taking Hotan Area as an Example

JIANG Yun-lu1，2，3，CAO Yue’e1，2，3，4，HE Zhong-fa1，2，3 et al

（1.College of Resources and Environmental Science， Xinjiang University， Urumqi，Xinjiang 830046;2.Institute of Arid Ecological Environment， Xinjiang University， Urumqi，Xinjiang 830046;3.Key Laboratory of Oasis Ecology， Ministry of Education， Urumqi，Xinjiang 830046;4.College of Environmental and Geographical Sciences， Shanghai Normal University， Shanghai 200234）

Abstract Based on the LUCC land use data of 2000， 2005， 2010 and 2015，the land use dynamics model and the land use transfer matrix were used to analyze the land use changes in the Hotan area of Xinjiang in the Kunlun Mountains for 16 years，and the gray prediction model was used to predict the land use types in Hotan area in the future 2025 and 2030.The results showed that in the past 16 years， land use in the Kunlun Mountains represented by the Hotan area had been in a period of rapid change.The change rate and extent of land use types reflected that the area of arable land and construction land had increased significantly during the various land use changes in Hotan from 2000 to 2015， increasing by 985 and 21 km2， respectively. The area of other types of land had decreased to varying degrees. The largest decrease was grassland， with a reduction of 555 km2， followed by unused land， with a reduction of 398 km2.From the perspective of land use dynamics， the single dynamics of cultivated land and construction land from 2000 to 2015 were relatively large， while the single dynamics of other land types were relatively small， which indicated that the area of cultivated land and construction land increased rapidly during the study period， and the area of other land types decreased to varying degrees.From the perspective of land use transfer， the transfer of land use types mainly occurred in 2000-2005 and 2010-2015， and the transferred land types mainly occurred in grassland-cultivated land， unused land-cultivated land and cultivated land-construction land.The spatial changes of land use transfer were mainly manifested in longitudinal distribution along rivers， occurring in oasis areas where humans live and around cities.The gray prediction model is used to predict the land use types of Hotan in 2025 and 2030， and the prediction accuracy is good. The results show that the area of arable land and construction land in Hotan will continue to increase in 2025 and 2030， and other land types will decrease to varying degrees.

Key words Kunlun Mountains;Hotan area;Gray prediction model;Land use;Spatial-temporal change

基金项目 国家青年基金项目（41601572）;中国科学院战略性先导科技专项（XAD20040400）。

作者简介 姜云璐（1997—），男，山东潍坊人，硕士研究生，研究方向：人文地理与干旱区资源环境协调。通信作者，副教授，博士，硕士生导师，从事国土资源评价与国土空间规划研究。

收稿日期 2021-05-21;修回日期 2021-05-31

土地利用变化是指在各种驱动力作用下土地利用的目的和方式发生改变的体现，与人类活动密切相关[1-2]。土地利用/覆被变化（LUCC）是人类活动与自然环境相互作用最直接的表现形式，与人类活动紧密相关。20世纪90年代以来LUCC已成为全球研究的核心课题之一[3]。近年来对于土地利用动态变化的研究，主要集中在时空分布及变化[4-8]、驱动力机制[8-15]及运用马尔科夫和灰色预测模型（GM模型）对土地利用变化进行预测等研究[16-20]。土地利用时空变化研究既能反映出研究区土地利用在时间序列上发生的变化，又可以直观反映空间变化的状况，对区域土地利用变化作出更详细的描述。新疆和田地区由于其特殊地理条件和历史原因，其生态环境对气候变化的响应及承受力具有较高的敏感性和脆弱性，因此，及时掌握该地区的LUCC变化，有效预测未来土地利用类型变化，对优化该区域土地利用模式，促进当地社会、经济和生态环境协调发展具有重要参考价值。

该研究对昆仑山区的典型干旱区绿洲即和田地区土地利用进行时空变化研究，同时运用灰色预测模型对和田地区2025、2030年土地利用类型进行预测，为未来土地资源优化配置及合理规划土地资源提供借鉴。

1 资料与方法

1.1 研究区概况

和田地区位于新疆维吾尔自治区南隅，介于78°17′～82°38′E、36°52′～37°19′N。南部为喀喇昆仑山及昆仑山，北部为塔里木盆地，地势自南向北倾斜，呈南高北低、西高东低的格局。仅有小部分的绿洲分布在昆仑山和塔克拉玛干沙漠之间的内流河流域附近，海拔1 300～1 450 m。属暖带内陆荒漠型氣候区，主要气候特点：四季分明、夏长秋短、冬不严寒，全年降水稀少蒸发量极大。春季、夏季多沙暴、浮尘天气，冬季少积雪，昼夜温差大，盛行西北风和西南风。年平均降水量为33.7 mm，且分布不均，年际变化大，一年中降水以5—7月份最多，约占总量的59%。年日照时数为2 602.2 h，可照时数为4 433.5 h，年平均日照百分率为59%。

1.2 数据来源

该研究所选2000、2005、2010、2015年4期数据来自中国科学院资源环境科学数据中心（http：//www.resdc.cn/），根据研究区实际情况参考LUCC二级分类体系对数据进行分类，通过地统计分析方法获得土地利用转移矩阵及空间变化图。

1.3 研究方法

1.3.1 单一土地利用类型动态度。

土地利用动态变化的研究方法主要是土地利用动态变化模型的建立，土地利用变化主要体现在土地利用类型变化、土地利用程度变化及土地利用变化的区域差异方面[21]。

单一土地利用动态度可描述区域一定时期内某种土地利用类型的变化速度，它对比较土地利用变化的区域差异和预测未来土地利用变化趋势都有积极的作用[22]。

其表达式如下：

K=Ub-UaUa×1T×100%（1）

式中，K为研究期内研究区某一类土地类型动态度;Ua、Ub分别为研究初期和研究末期土地类型的面积;T为研究时段。当T的时段为年时，K值是该土地利用类型的年变化率。

1.3.2 综合土地利用动态度。

综合土地利用动态度反映某区域一段时期内土地利用变化的速度[23]。

其表达式如下：

LC=ni=1ΔLUi-j2ni=1LUi×1T×100%（2）

式中，LC表示该区域土地利用年变化率;ΔLUi-j表示研究时段内第i类土地利用类型转为非i类土地利用类型面积的绝对值;LUi表示研究初期某一土地利用类型的面积;T的时段为年。

1.3.3 土地利用转移矩阵。

土地利用转移矩阵能够描述各种土地利用类型之间的转换情况，不仅可以反映研究初期和研究末期的土地利用类型结构，同时还可以反映研究时段内各类用地之间的相互转移变化情况，用来刻画区域土地利用变化方向以及研究期末各土地利用类型的来源及构成。转移矩阵中的变量可以为土地利用类型面积，也可以生成区域土地利用变化的转移概率矩阵，从而可以推测一些特定情境下区域土地利用的变化趋势[24]。其数学表达如下：

式中，Sij为研究期内土地利用状态;n表示研究区内土地利用类型数;i、j分别为转移前和转移后的土地利用类型。

1.3.4 灰色预测模型。

GM模型是对原始数据系列进行一阶累加生成，建立一阶线性微分方程模型，求得拟合曲线，从而对系统进行预测，过程如下[25]：

设x（0）为原始数据系列，对x（0）进行一次累加生成，得到一个生成序列x（1），由此可以建立白化形式的微分方程：

dX（1）dt+aX（1）=u（4）

用最小二乘法拟合得到参数a和u：

微分方程式（4）所对应的时间响应函数即为此次数列预测的基础公式，如下：

x（1）（k+1）=x（0）（1）-uae-ak+ua（8）

为提高GM（1，1）模型预测值的可靠度，采用后验差检验方法进行模型精度检验，后验差比值C和小误差概率P定义为：

C=S2S1（9）

P=P（|ε（k）-|<0.674 5S1）（10）

式中，S1为原始数据标准差;S2为预测数据标准差;ε为预测数据误差;为预测误差的均值。C的大小表明模型預测值和实际值之差的离散程度，C越小越好，表示原始数据离散大，而预测误差离散性小，则预测精度高;P越大越好，即小误差的概率越大，表示拟合精度较高。

2 结果与分析

2.1 土地利用变化幅度分析

区域土地利用面积和结构变化是研究土地利用变化的一个重要方面。对和田地区2000、2005、2010、2015年4个时期的土地利用面积变化幅度进行对比，结果发现（表1），

2000—2005年和田地区耕地和建设用地面积呈不断增加的趋势，其中耕地变化幅度最大，为13.85%，耕地面积增加了387 km2，建设用地变化幅度为1.43%，面积增加了4 km2;林地、草地、水域和未利用地面积均呈减小趋势，其变化幅度分别为-0.47%、-0.38%、-0.37%、-0.08%，其中草地面积减少最多，减少了209 km2，其次是未利用地，面积减少了135 km2，林地和水域面积分别减少了6、41 km2。

2005—2010年和田地区各地类面积变化较小，耕地、水域面积均有所增加，林地和建设用地面积保持不变，未利用地和草地面积不断减少;耕地变化幅度为1.51%，面积增加了48 km2，草地和未利用地变化幅度分别为-0.04%、-0.02%，面积分别减少20、29 km2。

2010—2015年和田地区土地利用变化较剧烈，其中耕地和建设用地变化较大，面积分别增加了550、17 km2，增长幅度分别为17.03%、6.01%;林地、草地和未利用地均有不同程度的减少，减少面积分别为15、326、234 km2，变化幅度分别为-1.19%、-0.59%、-0.13%。

整体而言，和田地区2000—2015年土地利用变化中，耕地和建设用地面积持续增加，且增加幅度较大，土地利用变化主要发生在2000—2005和2010—2015年，尤其在2010—2015年和田地区各类土地面积变化较剧烈。

2.2 土地利用动态度

从单一土地利用动态度（表1）来看，和田地区2000—2015年耕地和建设用地单一动态度呈正值，呈增长趋势;尤其是耕地，2000—2005年单一动态度为13.85%，2005—2010年为1.51%，2010—2015年为17.03%，表明耕地面积在2000—2005和2010—2015年增长显著，其中2010—2015年增长最快;建设用地在2000—2005年单一动态度为1.43%，2005—2010年单一动态度基本保持不变，2010—2015年单一动态度为6.01%，表明2010—2015年建设用地面积增长速度最快。林地、草地和未利用地单一土地利用动态度都呈负值，其中2000—2005年单一动态度分别为-0.47%、-0.38%、-0.08%，2005—2010年单一动态度分别为0、-0.04%、-0.02%，2010—2015年单一动态度分别为-1.19%、-0.59%、-0.13%;2000—2005和2010—2015年各类土地利用类型数量变化速度较2005—2010年快。

从综合土地利用动态度来看，2000—2005、2005—2010、2010—2015年综合土地利用动态度分别为8.32%、8.33%和8.31%，表明2005—2010年和田地区土地变化幅度较小，各类用地变化较其他2个阶段平缓;2010—2015年和田地区土地利用变化速度大于2000—2005年，2000—2005年土地利用变化速度大于2005—2010年。

整体而言，2000—2015年和田地区耕地单一土地利用动态度最大，土地面积变化较剧烈且增加幅度较大;建设用地比重较小，但动态度在2000—2005和2010—2015年增长幅度较大，尤其在2010—2015年，建设用地扩张较大;林地、草地动态度呈负，呈减小趋势，但草地面积减小数量较大;未利用地比重最大，虽然动态度较小，但实际减少面积较大。

2.3 土地利用转移矩阵

由表2可看出，2000—2005年和田地区土地利用转移主要表现在耕地—草地、林地—草地、水域—草地和未利用地—草地，草地面积的增加中耕地占21.15%，林地占3.85%，水域占13.46%，未利用地占59.62%。建设用地的增加全部来自对耕地的侵占。水域—未利用地、草地—未利用地转移分别占未利用地面积增加的61.11%、38.89%。

从表3可以看出，2005—2010年整体来看各类土地之间的转化并不频繁，主要表现在草地—耕地和未利用地—耕地，其中增加的耕地中草地占22.92%，未利用地占77.08%。其次为草地—未利用地，未利用地的增加全部来自草地面积的退化。

从表4可以看出，2010—2015年和田地区土地利用类型转化主要表现在林地—耕地、草地—耕地、水域—耕地、建设用地—耕地和未利用地—耕地。其中增加的耕地中草地和未利用地面积较大，分别占增加面积的58.39%和37.32%，其次为林地和水域，分别占增加面积的2.68%、1.25%，建设用地占增加面积的0.36%;草地面积的增加主要表现在耕地—草地、未利用地—草地的转移，其中增加的草地中耕地面积占60.00%，未利用地面积占30.00%;建设用地的增加主要来源于耕地—建设用地、未利用地—建设用地的转移，其中耕地占10.53%，未利用地占84.21%。

整体而言，2000—2015年和田地区土地利用转移主要集中在草地—耕地、未利用地—耕地之间的相互转化。耕地主要表现为持续增加的趋势，且耕地的增加主要来自对草地和未利用地的大量开垦，造成了草地、未利用地的大面积减少;一部分耕地的减少主要发生在城区周围，为了满足城市的发展和扩张，城市周边一部分地理位置较好的耕地被开发为建设用地。这是由于一方面因为人口的持续增加，有限的耕地不能养活更多的人口，因此把矛头指向肥力和水资源条件较好的草地和未利用地，对草地和未利用地进行大面积的开荒;另一方面由于“耕地保护红线”和响应“退耕还林还草”等政策因素，部分耕地因生产力下降被撂荒转为草地，而“耕地保护红线”的政策和为促进经济的发展等需求促使耕地不断增加。未利用地的减少主导因素是垦荒，其次是转变为草地、水域和建设用地。

2.4 土地利用空间变化分析

2000—2015年和田地区土地利用类型转移及空间变化主要表现出沿河流纵向分布和绿洲人类聚居区周围分布的特点，且主要集中在河流中游人类活動的绿洲区域，其主要呈现出如下特点：

（1）草地—耕地。

草地—耕地的转移主要在克里雅河、喀拉喀什河和玉龙喀什河中游，沿河流纵向呈片状及点状分布;另外尼雅河和皮山县境内发源于喀喇昆仑山的中小山溪性河流中游沿河流纵向呈小片及点状分布，其次在墨玉县城市周边及皮山县城市周围呈大片分布（图1）。这是由于和田地区近16年来人口的不断增加和外来人口的大量涌入，对耕地需求加大，一些水肥条件较好的草地被不断开发，在增加耕地的同时，生态环境也遭到一定程度的破坏，更加剧了和田地区本不乐观的生态状况。

（2）未利用地—耕地。

耕地的增加一方面来源于对草地的大量开垦，其次是未利用地的大面积开垦。和田地区未利用地转化为耕地沿河流纵向呈点状分布，主要集中在民丰县城市周边，和田县北部，洛浦县南部、墨玉县南部与昆仑山交接的平原地带，在皮山县主要在沿河流及城市、乡镇周边呈大面积分布（图1）。和田地区未利用土地资源广阔，在大量开发未利用地的同时，满足了农业和经济发展的要求。

（3）其他各类用地空间变化。

草地—未利用地的转化主要在田县城市东部及民丰县境内的牙通古孜河下游，呈片状及点状分布;林地—耕地的转化主要集中在策勒县城市周边;耕地—建设用地的转化主要在和田市及于田县，呈点状分布;未利用地—草地的转化主要发生在于田县和皮山县。

总之，和田地区土地利用类型空间变化特征与和田地区特殊的地理环境关系巨大，主要分布在自然条件相对优越的地理单元内。干旱区绿洲人类生产生活与水资源条件密不可分，因此土地利用空间变化主要沿河流呈片状及点状分布。耕地—建设用地的转化主要分布在于田县、和田市，呈点状分布。

2.5 土地利用灰色预测

基于和田地区实际情况，该区未利用地以较难利用的沙漠、戈壁和石质沙地为主，利用面积较稳定;水域包括河流湖泊等，季节性变化较大，在此不对这2种地类进行预测。选择2000、2005、2010、2015年4期数据，运用灰色GM模型，对未来2025和2030年土地利用类型面积进行预测（表5），参照灰色预测精度检验等级标准（表6），所有土地利用类型后验差比值C均小于0.65;林地、建设用地预测结果分别为0.57、0.62，勉强合格;草地C值为0.49，合格;耕地C值为0.28，精度等级为好。小误差概率P=1>0.95，表明模型预测精度为合格，预测值有较高的可信度。

3 结论

该研究选取干旱区典型区域和田地区为研究区，利用近16年土地利用数据，分析土地利用转移矩阵，并借助GM模型预测该研究区2025和2030年土地利用类型。主要结论如下：

（1）近16年和田地区各类土地利用变化较剧烈。从土地利用变化幅度来看主要表现为耕地和建设用地的持续增加及林地、草地和未利用地的不断减少，其中2000—2005和2010—2015年增减幅度较大。

（2）从单一土地利用动态度来看，2000—2005和2010—2015年耕地和建设用地单一动态度较高，而林地、草地和未利用地单一动态度较低，呈负值，但面积变化相对较大;2005—2010年单一土地利用动态度相对较小。从综合土地利用动态度看，和田地区土地利用变化速度较大，2005—2010年土地利用综合动态度较前后2个时段小，表明2005—2010年土地利用变化较小。但各研究时段相差不大，表明研究期间土地利用相对较稳定。

（3）土地利用转移发生在草地—耕地、未利用地—耕地、耕地—草地和耕地—建设用地之间。频繁转移时间在2000—2005和2010—2015年。耕地面积增加的主要来源是草地和未利用地的大量开垦，建设用地增加主要是对城区周边耕地的侵占和对未利用地的开发。

（4）和田地区土地利用空间变化主要沿河流以及沿交通线呈扁十字型分布的特点，与干旱区绿洲生产活动受水资源条件约束相契合，人类活动主要集中在面积较小的绿洲区域，耕地—建设用地、未利用地—建设用地空间变化呈点状分布的特点。

（5）利用GM灰色预测模型对和田地区2025和2030年土地利用类型进行了预测，预测精度良好，预测结果为2025、2030年耕地和建设用地面积不断增加，林地、草地面积都有不同程度减少。

参考文献

[1] BROGAARD S，ZHAO X Y.Rural reforms and changes in land management and attitudes：A case study from Inner Mongolia，China[J].Ambio，2002，31（3）：219-225.

[2] LAMBIN E F.Modelling and monitoring land-cover change processes in tropical regions[J].Progress in psychical geography，1997，21（3）：375-393.

[3] 陈佑启，杨鹏.国际上土地利用/土地覆盖变化研究的新进展[J].经济地理，2001，21（1）：95-100.

[4] 刘纪远，匡文慧，张增祥，等.20世纪80年代末以来中国土地利用变化的基本特征与空间格局[J].地理学报，2014，69（1）：3-14.

[5] 刘纪远，张增祥，徐新良，等.21世纪初中国土地利用变化的空间格局与驱动力分析[J].地理学报，2009，64（12）：1411-1420.

[6] 罗娅，杨胜天，刘晓燕，等.黄河河口镇—潼关区间1998—2010年土地利用变化特征[J].地理学报，2014，69（1）：42-53.

[7] 马才学，孟芬，赵利利.1990—2005年武汉市土地利用时空变化及其政策驱动因素分析[J].水土保持研究，2015，22（2）：117-122.

[8] 王薇，王昕，黄乾，等.黄河三角洲土地利用时空变化及驱动力研究[J].中国农学通报，2014，30（32）：172-177.

[9] 佟光臣，林杰，陈杭，等.1986—2013年南京市土地利用/覆被景观格局时空变化及驱动力因素分析[J].水土保持研究，2017，24（2）：240-245.

[10] 高凌寒，赵鹏祥，张晓莉，等.西宁市主城区土地利用时空变化驱动力分析[J].水土保持研究，2017，24（2）：234-239.

[11] 张月，郄睿卿，张蔚，等.吉林市土地利用时空变化及驱动力分析[J].测绘科学，2015，40（7）：73-77.

[12] 张军辉，顾蕾，岳永兵.沧州市土地利用时空变化及驱动力分析[J].中国国土资源经济，2014，27（12）：69-72.

[13] 于海影，韦安胜，陈竹君.基于RS和GIS的杨凌区土地利用变化及驱动力分析[J].水土保持研究，2014，21（5）：79-83.

[14] 吕桐，吴永波.基于RS的崇明岛土地利用变化及驱动力分析[J].中国农学通报，2014，30（5）：203-207.

[15] 史利江，王圣云，姚晓军，等.1994～2006年上海市土地利用时空变化特征及驱动力分析[J].长江流域资源与环境，2012，21（12）：1468-1479.

[16] 郭年冬，陈召亚，李恒哲，等.基于土地利用变化下的县域生态敏感性及灰色预测：以河北省平山县为例[J].水土保持研究，2016，23（5）：229-234.

[17] 李正，王军，白中科，等.贵州省土地利用及其生态系统服务价值与灰色预测[J].地理科学进展，2012，31（5）：577-583.

[18] 闫沛禄.基于信息熵的酒泉市土地利用结构分析及其灰色预测[J].甘肃农业大学学报，2011，46（6）：129-134.

[19] 王晓娇，陈英，齐鹏，等.基于信息熵的张掖市土地利用结构分析及其灰色预测[J].干旱区研究，2011，28（1）：92-97.

[20] 嚴明，廖铁军.基于Excel的灰色预测模型在土地利用规划中的应用[J].安徽农业科学，2007，35（12）：3627-3628.

[21] 李宁，岳彩荣.基于GIS与RS的肇源县土地利用变化分析[J].中国农学通报，2011，27（29）：224-228.

[22] 唐家荣，张世熔，吴光碧，等.基于RS和GIS的土地利用时空动态变化分析：以重庆永川区为例[J].山东农业大学学报（自然科学版），2014，45（4）：552-558.

[23] 李忠锋，王一谋，冯毓荪，等.基于RS与GIS的榆林地区土地利用变化分析[J].水土保持学报，2003，17（2）：97-99，140.

[24] 韩会然，杨成凤，宋金平.北京市土地利用变化特征及驱动机制[J].经济地理，2015，35（5）：148-154，197.

[25] 徐建华.现代地理学中的数学方法[M].2版.北京：高等教育出版社，2002：342-344.