基于GIS的县域尺度土地利用动态变化——以福建省建瓯市为例

基于GIS 的县域尺度土地利用动态变化
——
—以福建省建瓯市为例赵莉莉范勋承*
（邵阳学院城乡建设学院，
湖南邵阳422000）20世纪90年代，International Geosphere-Bisphere Progra-mme （IGBP ）和International Human Dimensions Programme
（IHDP ）提出的土地利用与覆被变化
（Land Use and Land Cover Change ，LUCC ）已成为目前全球环境变化研究的前沿和热门课题之一[1]。

土地利用变化是LUCC 研究的三大核心课题之一[2]，它深刻影响着全球可持续发展状况，因此土地利用变化及其对周围环境的影响已成为可持续发展研究中的重要研究领域[3-4]。

近
年来，国内学者就土地利用变化方面开展了诸多研究，
研究区域主要集中在生态敏感区或经济发达区域，
如马帅等[5]利用3S 技术对江淮生态经济区2000～2015年的土地利用景观格局时空
变化情况进行研究；
沈振涛等[6]以大汶河流域作为研究区域，利用土地利用相关模型对1985～2017年研究区的土地利用情况
进行了分析，
并对其2030年的土地利用变化进行了预测；常小燕等[7]以微山县的11个乡镇作为研究区，
对2000～2017年矿区的土地利用变化进行了研究；刘雨等[8]以山东省为研究区，探究
了土地开发和效益的耦合协同关系；甄江红等[9]
利用定量分析的方法，对呼和浩特市的城市用地扩展及驱动机制进行了深入
探讨。

以上研究大多将研究区域选定在大、
中尺度上，对于以县为研究对象的小尺度研究较少[10-11]。

然而近年来，
随着经济的不断发展，以县域为代表的小尺度区域城镇化速度明显加快，
对县域的土地利用情况展开深入分析，对于地方长期稳定发展具有重要意义。

建瓯市简称“芝城”，是福建省北部的农业大市，
省级园林城市、省级生态城市。

研究期间，建瓯市地区生产总值从1995年的24.70×108元增至2015年的198.67×108元，20年间翻了8倍
多；建瓯市历来人口基数大，
是南平市（闽北）所辖的10个地区中人口最多的城市，2015年常住人口45.1×104人，城镇化水平
48.1%。

目前，县域已经被作为推进新型城镇化的重点区域，
这势必引起土地利用结构的改变，如不予以重视，合理规划及利
用，最终会影响到该市的人居环境和生态安全。

因此，
本文利用1995和2015年2期遥感影像数据，
研究20年来建瓯市的土地利用状况，以期为该市土地资源的可持续利用提供科学参考。

1研究区概况
建瓯市属福建省南平市（闽北）下辖地区（图1），土地总面积4233km 2，是福建省陆地面积最大的县级市。

属东南沿海低山丘陵区，山地占土地总面积的57.3%，丘陵占29.9%。

境内河流发
基金项目：邵阳学院博士专项基金资助项目（17ZX03）；湖南省教育厅一般项目（19C1670）。

作者简介:赵莉莉（1985-），女，籍贯：黑龙江青冈，学历：博士，职称：助理研究员，从事土地利用与景观生态研究。

通讯作者:范勋承，学历：博士，职称：讲师，从事土地利用与3S 技术应用。

摘要：土地利用变化深刻影响着全球可持续发展状况，目前已成为全球环境变化研究的热门课题之一。

近年来，县域成为了
推进新型城镇化的重点区域，土地利用十分活跃，掌握县域土地利用状况对于区域的可持续发展具有重要意义。

因此，
本文以县级市-福建省建瓯市作为研究案例，对该市1995耀2015年的土地利用情况进行定量分析，以期为小尺度区域土地的合理利用提供
参考。

研究结果表明：（1）总体上看，建瓯市的土地利用变化速度较为平稳。

从各土地利用类型上看，20年间建瓯市各土地利用变
化活跃，其中耕地是面积增加最多的用地类型，动态变化最为剧烈，动态度达42.4%。

（2）研究期间，建瓯市的土地资源开发强度和
集约化程度均在全国水平之上，土地利用的综合水平较高，土地利用开发指数和垦殖指数在20年间均呈持续增长态势。

（3）建瓯
市耕地面积的大量增加，主要得益于林地、
未利用地和草地的大量转入。

关键词：土地利用与覆被变化；小尺度；可持续发展；闽北
Abstract:At present,land use change has become one of the hot topics in the study of global environmental change,because it profoundly affects the global sustainable development.Since the county has become a key area for promoting new urbanization in recent years,land use has been very active.Therefore,mastering the land use status of the county is of great significance to the sustainable development of the region.This paper takes Jian'ou City (county-level City),Fujian Province as a research case to quantitatively analyze the land use of this city from 1995to 2015,in order to provide a reference for the rational use of land in small-scale areas.The results showed that (1)overall,the rate of land use change is relatively stable in Jian'ou City.From the perspective of land use types,the land use changes have been active in the past 20years.Among them,cultivated land is the land type with the largest increase in area,with the most dramatic dynamic changes,with a value of 42.4%.(2)During the study period,the intensity and intensification of land resources in Jian'ou City were above the national level,and the comprehensive level of land use was high.Both the land use development index and the reclamation index have shown continuous growth in the past 20years.(3)The large increase in the area of cultivated land in Jian'ou City is mainly due to the massive transfer of forest land,unused land and grassland.
Key words:Land Use and Land Cover Change;Small scale;Sustainable development;Northern Fujian Province 中图分类号:TP79，K921文献标识码:A 文章编号:2096-4390(2021)07-0004-04
育完全，属闽江建溪流域，总长1079km。

建瓯市属亚热带海洋性季风气候，四季分明，年均气温14-20℃。

截至2019年末，户籍人口54.91×104人，是闽北人口最多的县级市，实现地区生产总值（地区GDP）291.72×108元。

图1研究区示意图
2研究方法
2.1数据来源与处理
根据研究目的和研究意义，同时考虑到数据的可获取性等因素，建瓯市1995年Landsat-5TM和2015年Landsat-7ETM遥感影像来源于地理空间数据云平台（/），空间分辨率均为30m×30m。

参照我国土地利用分类系统并结合南平市的地理特征，本文将研究区划分为林地、草地、耕地、建筑用地、水域和未利用地共6种土地利用类型。

经前人的经验及研究分类的需要发现，RGB组合为波段5、4、3能够更加清晰的区分本研究中需要的6种土地利用类型。

因此，利用ENVI Classic中的Layer stacking功能将每期遥感影像图按照波段5、4、3进行合并。

而后，采用直方图匹配法将图像进行镶嵌，导入研究区shp矢量图，执行裁剪操作，最终得到所需图像及相关数据。

2.2土地利用变化的研究方法
2.2.1土地利用动态度
土地利用动态度模型是表征土地利用变化速度的模型。

表示的是某种地类在一定时间段内面积变化的快慢，正值表示面积增加，负值反之。

该模型包括两种类型，分别是单一土地利用动态度和综合土地利用动态度。

2.2.1.1单一土地利用动态度
该指标反映的是研究期间某种地类的动态变化程度[12]，即某种地类增长或减少的速率，其数学表达式为
（1）
式中，K为某一地类的动态度；U a和U b则分别表示开始和结束时该地类的面积；T则为研究的时间长度，单位为a，下同。

2.2.1.2综合土地利用动态度
综合土地利用动态度可反映研究区内各地类综合的变化速率，是研究区总体土地利用变化指标，其值越低，表明土地利用变化越平稳[13]。

其数学表达式为：
（2）
式中，IR j为研究区综合动态度；n为地类的数量；Z ija、Z ijb是研究区某地类的初始面积和研究末期的面积。

2.2.2土地利用程度
土地利用程度可以表征人类开发土地的强度，是多种土地利用类型变化的结果[14]，其模型如下：
（3）
式中，L j代表土地利用程度；n代表土地利用分类数；A i为第i级土地利用的分级数；C i为第i级土地利用分类占总面积的百分比。

本研究利用刘纪远先生[15]提出的土地利用程度的分级数，并结合建瓯市土地利用类型的特点，把未利用地的分级数定为1，林地、草地和水域定为2，耕地为3，建筑用地为4。

土地利用程度变化量（率）的数学表达式分别为：
（4）
（5）
式中，为土地利用程度的变化量；R是土地利用程度的变化率；L a和L b为研究开始和研究结束时的土地利用变化程度指数；A i为第i级土地利用的分级数；C i,a和C i,b分别为研究开始和研究结束时第i级土地利用分类占总面积的百分比。

如果，表明研究区土地利用处于发展阶段，反之，则处于调整期或衰退期。

2.2.3土地利用结构
土地利用结构可用结构状态特征指数来表示，即可用一种或几种地类占研究区总面积的比值来表示[15]。

根据研究目的需要，本章对土地利用结构特征指数分别定义为土地利用开发指数和土地利用垦殖指数。

具体情况如下：
土地利用开发指数的数学表达式如下：
（6）
式中，b i为建筑用地的面积；B则是开发指数。

该指数是表示土地被利用为建筑用地的比例。

土地利用垦殖指数计算公式如下：
（7）
式中，a i为耕地类型的面积，A为研究区面积，K是垦殖指数，该公式表示耕地占研究区面积的百分比。

2.2.4土地利用转移
矩阵
土地利用转移矩阵是目前研究土地利用变化特征、方向、结构的常用方法，它不仅可以表达土地利用的结构特征，还可以表示研究时段内的各地类的来源及组成[15]。

除此之外，转移矩阵还可以生成转移概率矩阵，它可以按照百分比进行排序，从该排序中可以看出引起该地类变化的主要土地利用类型和次要类型，从而分析出引起该种地类变化的原因。

土地利用转移矩阵的数学公式为：
（8）
式中，P代表面积；n代表土地利用类型的个数；i和j分别代表研究时段初期与末期的土地利用类型。

3结果与分析
3.1土地利用类型的面积变化
利用ArcGIS10.2软件的统计分析工具对1995年和2015年的遥感解译数据进行统计，得到两期各地类面积，整理计算得到各地类面积占整个研究区面积的比例变化图（图2）。

从图2可知，20年间建瓯市土地利用类型间变化频繁，总体表现为“三增三减”具体为耕地、建筑用地和水域面积有所增加，其中耕地增加比例最大，达7.70%。

建瓯市的人口数量在闽北地区排名第二位，仅次于延平区，2015年常住人口为45.1万人，同时，建瓯市也是闽北基本农田主要集中区，因此，在保护基本农田的基础上，不断进行土地整理和复垦开发，以满足人口增长的需要。

这是建瓯市在20年中耕地面积增加的主要原因。

研究期间，草地、未利用地和林地呈现减少趋势，面积比例依次减少了5.30%、3.99%和0.72%。

主要因为这些面积减少的用地类型部分转入了耕地所致。

图21995～2015年建瓯市不同地类面积比例变化
3.2土地利用类型的动态度
由图3可知，研究期间建瓯市各土地利用类型单一动态度值有正有负，说明其土地利用类型处于不断调整中。

与其他用地类型相比，耕地和水域的动态度为正值，且耕地动态变化最大，动态度大42.4%，说明此期间建瓯市的耕地面积增加显著。

除耕地和水域外，其余用地类型的动态度均为负值，说明建瓯市在此期间其余用地类型面积有所减少，具体减少顺序由多到少依次为：未利用地（-4.1）＞草地（-3.37）＞建筑用地（-1.04）＞林地（-0.61）。

这是由于20年间，随着经济的不断发展，建瓯市对未利用地的开发强度持续加大，未利用地不断向其他用地类型转换的结果。

根据公式（2）的综合土地利用动态度数学模型，计算出建瓯市1995～2015年综合动态度为0.56，说明建瓯市总体土地利用变化速度较为平稳。

图3建瓯市各土地利用类型单一动态度
3.3土地利用程度
按照公式（3）的土地利用程度模型对建瓯市1995年和2015年的土地利用程度进行了计算，其值分别为202.70和218.19。

从计算结果可以看出，20年间建瓯市土地利用的深度和广度不断增大。

另外，与全国平均水平201.86[16]相比，无论是1995年还是2015年都高于全国平均水平，说明建瓯市土地资源的开发强度和集约化程度在全国平均水平之上。

利用公式（4）和（5）计算可知，建瓯市的土地利用程度变化量和变化率分别为15.49和0.08，其值均大于0，说明该市土地利用较为活跃，土地利用的综合水平较高，处于蓬勃发展期。

3.4土地利用结构
根据公式（6）和（7）计算得到研究区土地利用开发指数和垦殖指数（图5）。

从图4可以看出，1995-2015年间，建瓯市土地开发指数从1995年的2.78增至2015年的4.47，说明20年间建瓯市土地利用开发强度不断增大，新增建筑用地持续增加。

土地垦殖指数从2.65增至9.83，20年间翻了3倍多。

建瓯市是福建省陆地面积最大、闽北人口最多的县级市，同时，建瓯市还是闽北基本农田集中区，在国家保护基本农田政策的影响下，20年间该市耕地面积不断增加，因此垦殖指数持续增大。

图41995～2015年建瓯市土地利用结构
特征
3.5土地利用类型的转移矩阵
从1995～2015年建瓯市土地利用变化转移概率矩阵上看（表1），建瓯市林地保持率是66.96%，转为林地最高的用地类型是水域，转入率为61.05%，建筑用地次之，为58.45%。

建瓯市耕地保持率为34.62%，加之其他地类的转入也较高，其中，有21.64%的林地转为了耕地，而林地是研究区土地面积占比最大的用地类型，因此研究期间建瓯市耕地面积增加最多。

各土地利用类型具体转出概率排序为：未利用地（99.76%）＞草地（92.63%）＞建筑用地（88.29%）＞水域（71.06%）＞耕地（65.38%）＞林地（33.04%）；除了第一转出方向是林地外，未利用地和草地的主要转出方向还有耕地，其向耕地的转入概率分别为33.87%和35.52%。

表11995～2015年建瓯市土地利用变化转移概率矩阵（%）
4主要结论
土地利用变化可以反映某地区的人类活动对地表覆盖物的影响程度，尤其是最近几十年，县域城镇化进程加快，人类对土地资源开发和利用的规模猛增。

因此，以县级市-福建省建瓯市作为研究案例，利用土地利用相关计算模型，对该市20年间的土地利用情况进行定量分析，以期为该市及其它县域的土地可持续利用提供科学参考。

其主要研究结论如下：
4.120年间，建瓯市各土地利用变化总体表现为“三增三减”，其中耕地和草地分别是面积增加和减少最大的用地类型。

与其他用地类型相比，建瓯市虽然耕地动态变化最为剧烈，其值达42.4%，但总体而言，该市的土地利用变化速度较为平稳。

4.2研究期间，建瓯市的土地资源开发强度和集约化程度均在全国水平之上，土地利用的综合水平较高。

由于耕地和建筑用地的面积均呈增加趋势，因此，其土地利用开发指数和垦殖指数在20年间均呈持续增长态势。

4.3建瓯市耕地保持率虽不高，但林地、未利用地和草地的大量转入，使得20年间耕地成为了面积增加最多的用地类型。

参考文献
[1]Canadell J G,Steffen W L,White P S.IGBP/GCTE terrestrial transects:Dynamics of terrestrial ecosystems under environmental change[J].Journal of vegetation science，2002,13 (3):298-300.
[2]Ii B L T.Toward integrated land-change science:advances in1.5decades of sustained international research on land-use and land-cover change[J].Gaia:okologische perspektiven in natur-,geistes-und wirtschaftswissenschaften，2002,10(4):21-26.
[3]李秀彬.全球环境变化研究的核心领域冥冥土地利用/土地覆被变化的国际研究动向[J].地理学报，1996(6):553-558. [4]Ii B L T,Robbins nd-change science and political ecology:similarities,differences,and implications for sustainability science[J].Annual review of environment&resources，2008,33(1): 295-316.
[5]马帅,程浩,林晨,等.江淮生态经济区土地利用景观格局时空变化[J].水土保持研究,2021,28(1):292-299.
[6]沈振涛,赵强,高倩,等.大汶河流域土地利用变化特征及预测[EB/OL].济南大学学报(自然科学版),2021,(4):1-9[2021-02-04]. https:///10.13349/ki.jdxbn.20210201.004.
[7]常小燕,李新举,刁海亭.矿区土地利用时空变化及重心转移分析[EB/OL].内蒙古农业大学学报(自然科学版):1-11[2021-02-04]./kcms/detail/15.1209.S.20210126.1550.018. html.
[8]刘雨,朱晓伟,刘凯,等.山东省土地开发度与利用效益的时空耦合协调研究[EB/OL].西安理工大学学报:1-12[2021-02-04]. /kcms/detail/61.1294.N.20210201.1612.002.html.
[9]甄江红,冯琰玮,峰一.呼和浩特城市用地扩展及其驱动因素分析[J].内蒙古农业大学学报(自然科学版),2020,41(6):31-39. [10]刘晶晶,王静,戴建旺,等.黄河流域县域尺度生态系统服务供给和需求核算及时空变异[J].自然资源学报,2021,36(1):148-161.
[11]胡毅,乔伟峰,万懿,等.江苏省县域土地利用效益综合评价及其分异特征[J].经济地理,2020,40(11):186-195.
[12]李毅,胡文敏,易敏,等.基于土地利用变化的湖南省生态系统服务价值分析[J].林业经济,2017,39(07):91-94+102. [13]赵莉莉,何东进,覃德华,等.2000~2014年闽东海岸带土地利用动态变化[J].热带作物学报,2017,38(5):803-810. [14]宋开山,刘殿伟,王宗明,等.1954年以来三江平原土地利用变化及驱动力[J].地理学报,2008(1):93-104.
[15]刘纪远.中国资源环境遥感宏观调查与动态研究[M].北京:中国科学技术出版社,1996.
[16]侯西勇,徐新良.21世纪初中国海岸带土地利用空间格局特征[J].地理研究,2011,30(8):1370-1379.
1995 2015 建筑
用地
林地水域耕地未利
用地
草地
建筑用地11.71 6.29 4.55 8.70 5.72 4.90 林地58.45 66.96 61.05 43.95 51.01 50.20 水域 6.04 1.66 28.94 4.95 0.97 1.76 耕地21.09 21.64 4.96 34.62 33.87 35.52 未利用地0.14 0.18 0.02 0.47 0.24 0.25 草地 2.57 3.27 0.48 7.31 8.19 7.37。