1. 引言
改革开放40年来,中国的经济实现飞速的增长,随着工业化进程不断推进,经济和产业的集聚,我国的城镇化水平不断提高,人口持续性往城市地区集聚,目前我国已成为全球城镇化的主动力,处于城市快速发展阶段。作为人类集中居住和生产的载体,城镇随经济、社会的发展,受到包括人口集聚、经济、交通发展等各种要素的影响 [1],而城镇数量增加,城镇规模逐渐扩大,城镇用地的空间扩展日益成为现在乃至将来土地利用变化的主要特征 [2],因此,定量化度量城镇扩张动态变化对于了解预测城镇发展是有必要的。
目前国内外对于城镇用地扩张动态分析的已经有了大量的相关研究,早期国外许多学者利用TM、SPOT等遥感数据,结合GIS技术,对城镇扩张进行了研究 [3] [4],国内利用RS和GIS进行城镇扩展研究相对较晚,但研究进展很快 [5],目前研究成果主要集中于城镇空间结构变化驱动力研究、城镇空间扩展规律及时空机制的研究、城镇空间扩展预测研究等 [6] [7],综合城镇空间扩展的国内外研究现状发现,研究尺度从单个城市(主要为省会城市)扩展到区域城市群,以宏观层面对大城市、城市群的扩展规律、特征、驱动力、时空机制等进行了研究,而对于微观尺度的城镇扩张研究仍存在很多不足。
综上所述,本文以武汉市汤逊湖流域为研究区,利用优于空间分辨率为1 m的遥感影像提取研究区2002、2006、2010、2015、2019年5个年份的城镇用地信息,采用扩张强度指数、扩张动态度、景观特征指数、路网密度等对城镇用扩张进行动态分析,并综合经济、人口等数据探索城镇用地扩张的驱动因子。
2. 研究区域概况与遥感数据源
汤逊湖流域位于武汉市东南部,跨江夏区、洪山区、东湖高新区,地理坐标介于114˚13'0.386''E~114˚27'40.977''E、30˚19'13.838''N~30˚29'51.839''N之间,总面积240平方公里,地理位置见图1。汤逊湖周边地貌属于丘陵地段,经江汉平原向大别山低山丘陵过度地带,岗垄环抱,北部多山。地形属于残丘性河湖冲积平原,地势平坦低洼,气候属热带季风气候,常年雨量充沛,雨热同季,降水多集中在6~8月。
研究中使用的数据源为空间分辨率优于1 m的遥感影像数据,包括QuickBird、WordView-3以及GF2数据,为了更清晰地表征城镇建设用地,影像尽量选择云量较少、质量较高的数据,若出现影像缺失或云量较高的情况,则选择前后两年内同区域的高质量影像,数据源情况如表1。影像的预处理过程主要包括:大气校正、几何校正、图像配准、影像融合、镶嵌处理以及图像裁切等,预处理过程决定了遥感影像产品的质量以及解译的精度。
Figure 1. Tangxun Lake basin in the Wuhan metropolitans
图1. 武汉市汤逊湖流域
3. 遥感影像信息提取
遥感图像的解译分类主要是根据图像像元所反映地物的光谱特征和空间特征进行分析,将影像中的每一个像元按照一定的算法或特定的规则划分为不同的类别,从而实现遥感图像的分类 [8]。目前比较常用的遥感影像信息提取方法包括传统的目视解译、监督分类与非监督分类以及新型的人工神经网络、面向对象分类与遥感专家解译系统等。由于研究区范围较小,采用的影像分辨率较高,为提高影像信息提取精度,采用的是人机交互解译的方法。遥感影像信息提取流程如图2,对研究区五个时相的城镇用地解译分类结果进行处理,使用ArcGIS软件出图后的结果如图3。
Figure 2. Flow chart of remote sensing image information extraction
图2. 遥感影像信息提取流程图
Figure 3. Urban land use in Tangxun Lake
图3. 汤逊湖流域城镇用地图
4. 城镇用地扩张特征分析
对研究区五个时期城镇用地面积进行统计分析,其城镇用地面积由2002年的29.5 km2迅速增长到2019年的106.73 km2,计算出汤逊湖流域从2002~2019年4个时段城镇用地面积的增长量,其中2010~2015年面积增长量最大,结果如下(图4)。
Figure 4. Statistical graph of urban land expansion from 2002 to 2019 (left: area, right: growth in various stages)
图4. 2002~2019年各阶段城镇用地统计图(左:面积,右:增长量)
4.1. 城镇用地扩张强度指数
城镇扩张强度指数(Urbanization intensive index, UII)指某一地区在城镇扩张过程中城镇用地的扩张面积与该地区土地总面积的百分比,即某一地区在面积标准化下城镇用地面积的年平均增长速率 [9],城镇扩张强度指数能够反映城镇用地面积扩张的速率,通过城镇面积的变化可直观地体现出城镇扩张演变的强弱趋势。UII指数越高,城镇扩张越快。UII的计算公式如下:
(1)
式中:UII为扩张强度指数,Sb、Sa分别表示在时间点b、a时城镇用地的面积(km2),A表示研究区总面积,T表示间隔时长。本研究所参考的UII划分标准为:0~0.28为缓慢扩张,0.28~0.59 为低速扩张,0.59~1.05为中速扩张,1.05~1.92为快速扩张,大于1.92为高速扩张 [10],根据UII计算公式以及参考的划分标准计算得到结果如表2。
Table 2. Expansion intensity index from 2002 to 2019
表2. 2002~2019年各时间段扩张强度指数
从表2可以看出,汤逊湖流域从2002~2019年城镇用地各时间段一直处于高、快速扩张类型,综合分析2002~2019年城镇用地扩张,其为快速扩张型,而扩张强度在2010~2015年达到顶峰后呈下降趋势。近年来,随着工业化进程的不断推进以及国家经济的快速发展,武汉市的社会经济发展也加快了脚步,汤逊湖流域作为武汉市快速发展的典型例子,在近二十年城镇用地实现了飞速扩张,其中工矿企业用地的增加也带动了地区经济的发展以及人口的增加,而地区经济的发展以及人口的增加促使了城镇用地的扩张,以此形成了正反馈机制。
4.2. 城镇用地扩张动态分析
扩张动态度是指在某一特定时间段内研究区的某一土地利用类型的在数量上的动态变化状况 [11],可用来描述土地利用变化的速度 [12],计算公式如下:
(2)
式中:K表示城镇用地扩张动态度,Sb、Sa表示研究区研究末期与研究初期城镇用地面积,T表示时间间隔。本研究所参考的扩张动态度划分标准为:高速扩张型(>20%);快速扩张型(14%~20%);中速扩张型(8%~14%);缓慢扩张型(<8%) [13],根据参考划分标准以及计算公式计算得到各时间段城镇用地扩张动态度,如表3。
Table 3. Expansion dynamics in various time periods from 2002 to 2019
表3. 2002~2019各时间段扩张动态度
根据表3,2002~2019年城镇用地年均变化率为15.73%,各时间段年均变化率在2002~2006年达到峰值后降低,综合来看2002~2019年城镇用地扩张属于快速扩张型。建设初期,由于的经济快速发展以及人口增加,相应的基础配套设施、房地产业、工矿企业等其他建设得到一定程度的提高,对城镇用地的需求增加,这一时期呈快速扩张趋势,而后扩张年均变化率降低。
4.3. 景观格局指数
景观格局(Landscape Pattern)是指景观组成单元类型、异质性、多样性及其空间分布规律,是大小和形状不一的景观嵌块体在景观空间上按照一定的规律组成的景观格局 [14] [15],影响景观格局的变化因素多种,包括人为及自然等。通过利用景观格局空间特征指数,可以分析该区域空间变异特征以及各种生态过程在该区域作用的结果。为了探究城镇用地扩张过程中出现的空间现象,本文采用的景观格局指数有:①景观斑块个数(NP) ②景观斑块面积(CA) ③斑块密度(PD) ④平均斑块面积(AREA_MN) ⑤聚合度指数(AI) ⑥连接度指数(COHESION),如表4。
Table 4. The meaning and formula of the landscape pattern index feature to be adopted
表4. 拟采用的景观格局指数特征含义及公式
计算各时间景观格局指数,如表5。
Table 5. Landscape pattern index from 2002 to 2019
表5. 2002~2019年景观格局指数
由表5可知,斑块数量、斑块密度逐渐减小,斑块面积、平均斑块面积、聚合度指数、连接度指数逐渐增大,在研究内,斑块的集聚程度增大,城镇用地的小斑块数量减少,城镇用地斑块不断融合吞并,城镇用地景观连通性增加,城镇用地空间格局逐步紧凑型,从而形成聚集性斑块。其中2002~2015年斑块数量、斑块密度急剧下降,城镇用地面积、平均斑块面积急剧上升,斑块的聚合度指数和连接度指数较大幅度增加,在这一期间,城镇用地扩张迅速,城镇用地紧凑度增加,城镇属于成长阶段,2015~2019年景观格局指数变化幅度较小,城镇用地逐渐趋向完整、集聚的趋势,空间格局也越来越紧凑,形成成熟的城镇建成区。
4.4. 道路网密度分析
道路网密度是指在一定区域内道路网总里程与该区域面积的比值,是评价某一地区交通状况、衡量城市道路建设水平的重要指标 [16],道路网密度可反映规划用地单位面积内所拥有的道路长度的指标 [17],道路网密度越大,可达性越好,交通更便捷。计算公式如下:
(3)
式中:D表示研究区域道路网密度(Km/Km2),L表示研究区道路长度(Km),A表示研究区面积。计算得到各年份的道路网密度,如表6。
Table 6. Road network density from 2002 to 2019
表6. 2002~2019年道路网密度
从表6看出,从2002~2019年,随着城镇用地的扩张,道路网密度不断增加,道路网更密集,道路网络更复杂,具有良好的空间渗透性,可达性较高,便捷程度更高,集聚人流更多;在城市环线和绕城高速附近,各道路与对外交通拥有较多的道路节点,被到访机会较高,对人流产生更强的吸引力,促使商业集聚以及区域经济发展。
5. 研究区域城镇用地扩张因素分析
城镇扩张的本质是为了满足人类发展的需求,不断进行土地资源类型和数量上的调整过程 [18],城镇扩张的同时能够促进经济发展以及农业人口向非农业人口的转化以及人口的增长,进而促使城市化进程。影响城镇扩张驱动机制的形成因素有多种,包括自然地理环境、人口、经济以及交通等。
决定城镇扩张的因子主要是经济、人口等。随着社会主义市场经济体系的建立,全国经济表现出向上的趋势,城市化进程不断加快。从21世纪起,武汉市产业结构不断优化,改革力度增强,经济实现了快速发展。武汉市地区生产总值从2002年的1492.72亿元上升到2019年的16,223.21亿元,比2002年增加了11倍,如图5所示。
Figure 5. Urban land area and Wuhan’s GDP in Tangxun Lake Basin from 2002 to 2019
图5. 2002~2019年汤逊湖流域城镇用地面积与武汉市生产总值图
从图可以看出,武汉市GDP增长与汤逊湖流域城镇用地扩张呈正相关,伴随着经济的发展,第二、三产业占比增加,大量的工矿企业在汤逊湖流域集聚,产业集聚促使城镇用地扩张。经济增长的同时人口不断增加,而人口增长对土地的需求是土地的本质需求,城镇用地是人类生活、工作、学习所依赖的场所,从图6可以看出自2002年起武汉市人口数量呈连续不断增长的趋势,常住人口数量从2002年的823.70万人增长至2019年的1108.10万人。
Figure 6. The number of permanent residents in Wuhan from 2002 to 2019
图6. 2002~2019年武汉市常住人口数量图
人类的生存活动需要占据一定的空间,人口的增长与城镇用地的变化存在一定的显著关系,随着城市化进程的加快,农业人口向非农业人口转化的速度加快,从而促使房地产业的快速发展,进一步加快城镇用地面积的扩张。
6. 结论与展望
本文以2002、2006、2010、2015、2019年五个时相的遥感影像为基础,结合遥感与地理信息技术,通过对遥感影像数据的处理以及城镇用地的信息提取,研究了汤逊湖流域近17年来城镇用地扩张的动态性。研究表明,RS与GIS技术是研究城镇用地动态变化的有效手段;2002~2019年,汤逊湖流域城镇用地增加了77.68 km2,道路网密度由0.67增加到4.11,扩张强度指数与城镇扩张动态度分别为1.90与15.73%,属于快速扩张型,并计算得到各时间点的城镇用地景观格局指数,结合城镇扩张指数、城镇扩张动态度以及道路网密度,分析得出2002~2015年汤逊湖流域城镇用地快速扩张,城镇用地紧凑度增加,城镇属于成长阶段,2015~2019年景观格局指数变化幅度较小,城镇用地逐渐趋向完整、集聚的趋势,空间格局也越来越紧凑,城镇逐渐走向成熟阶段;对汤逊湖流域城镇用地扩张的驱动因素进行综合分析,2019年武汉市的生产总值是2002年11倍,常住人口数量增加至1108.1万人,经济的发展与人口增长促进了城镇用地扩张。
本文虽然从多个角度分析了汤逊湖城镇用地扩张动态性,但是本文对于景观格局指数的研究仅限于城镇用地单一单元,并没有分析汤逊湖流域的总体景观格局以及景观格局演变,总体景观格局以及景观格局的演变分析有利于研究汤逊湖流域生态环境变化。由于缺乏汤逊湖流域的人口、经济、产业变化数据,只能选取整个武汉市的经济、人口数据研究城镇用地扩张的驱动因素,在一定程度上,仍存在局限性。后期将结合景观格局演变、景观扩张指数以及多尺度的梯度分析城镇用地扩张特征的时空差异性。