土地资源是人类赖以生存并为人类生产生活提供所需基础物质的不可再生资源^[1-2], 其可持续利用是社会可持续发展的关键^[3-4]。土地适宜性是土地利用情况和土地利用规划的基本依据^[5], 目前对土地适宜性的研究多集中于损毁矿山的土地复垦适宜性评价^[6-7], 而对建设用地开发适宜性研究少, 且研究区多为山地等地势较高的地区, 如赵筱青等^[8]对云南喀斯特山区进行了建设用地开发适宜性评价及分区研究。因此, 对于经济社会快速发展、城镇化进程加快和城乡土地利用结构失衡问题突出的城市地区, 亟需开展以城乡建设用地开发适宜性评价为内容的相关研究。

土地承载力是指在一定的物质生活水准下, 土地资源所能承载的最大人口和社会活动规模^[9-10]。20世纪80年代末, 我国开始展开围绕人、地、粮关系的土地资源承载力研究^[11-13], 随着土地开发程度的加深^[14-15], 土地利用空间不均衡等问题突出, 相关研究逐渐拓展到土地综合承载力^[16-18]。近年来, 随着土地承载力评价受到愈来愈多的关注, 其评价对象、评价方法、评价体系等都愈加丰富。评价对象方面, 主要为城市土地承载力^[19]、区域土地承载力^[20-22]以及特殊地区土地承载力^[23-24]。研究方法上, 系统动力学^[25]、生态足迹模型^[26]、多指标综合评价^[27-28]等方法被广泛应用到不同地区的土地资源承载力评价。从评价体系来看, 土地承载力评价主要从环境容量^[29]、生态承载力^[30]、自然基础要素^[31]等方面选取评价指标, 较少考虑到区位条件及人口密度等社会经济因素, 无法反映在人类活动对土地承载力的影响。

不同学者分别从不同方面开展了土地适宜性和土地承载力评价研究, 但以土地适宜性评价为基础进行城乡建设用地承载力研究还较少, 更缺少城乡建设用地承载力和现状开发状态之间的空间匹配关系研究。因此, 理清城乡建设用地承载力和现状开发状态之间的空间匹配关系, 对于合理利用土地和优化国土空间布局, 以及加快生态文明建设和实现区域可持续利用具有重要的意义^[32-34]。以江西省南昌市为例, 在城乡建设用地开发适宜性评价基础上, 测算城乡建设用地承载力和可承载城乡建设用地丰度, 并探讨城乡建设用地承载力和现状开发状态之间的空间匹配关系, 以期为南昌市城乡建设用地合理开发利用和国土空间优化提供参考。

1 研究区概况与数据来源 1.1 研究区概况

南昌市是江西省的省会, 地处江西中部偏北, 位于115°27′—116°35′E、28°09′—29°11′N之间, 总面积达7402.36 km², 城乡建设用地面积为1237.57 km², 占总面积的16.72%, 下辖五区(东湖区、西湖区、青云谱区、湾里区、青山湖区)和四县(南昌县、新建县、安义县、进贤县), 地理位置如图 1所示。境内地形东南平坦, 西北丘陵起伏, 气候属于亚热带湿润季风气候, 光照充足, 雨水充沛, 自然资源丰富。近年来, 南昌市经济社会发展迅猛, 2018年全市生产总值达到5274.67亿元, 常住人口达554.55万人, 城镇化率74.2%。

1.2 数据来源

本研究数据主要包括：南昌市2017年土地利用变更调查数据库、永久基本农田划定成果、土地利用总体规划(2006—2020年)数据库、生态红线划定成果等, 还有江西省地面塌陷分布图、地质灾害调查评价图、南昌市统计年鉴、道路矢量数据等, 以及通过地理空间数据云(http://www.gscloud.cn/)获取的数字高程模型(DEM)(空间分辨率为30 m)。通过ArcGIS 10.2软件, 将数据投影到统一坐标系下, 并划分为30 m×30 m的格网作评价单元, 以便对数据进行空间分析。

2 研究方法 2.1 研究思路

随着经济快速发展和城镇化不断加快, 大量耕地开发为建设用地, 国土空间开发利用格局发生明显重构, 自然资源的过度开发以及土地利用不合理等使区域空间开发失衡问题日益凸显。为解决区域土地利用空间失衡问题, 协调人口、资源和社会发展的关系, 需要开展城乡建设用地承载力和现状开发状态之间的空间匹配关系研究。本文研究思路如图 2所示, 具体如下：

(1) 指标选取及量化处理：适宜性评价从生态安全、自然条件、社会经济等方面选取11个评价指标构建评价体系, 生态红线、永久基本农田区、地面塌陷为强限制性指标, 高程、蓄滞洪区、土地利用类型、区位条件为较强限制性指标, 坡度大于25°为强限制性指标, 小于25°为较强限制性指标, 并对指标赋值, 进行量化处理。

(2) 划分适宜性等级：以ArcGIS 10.2软件为处理平台, 进行建设用地开发适宜性评价, 并通过SPSS 22.0软件对适宜性分值进行聚类分析, 划分适宜性等级。

(3) 计算城乡建设用地承载力：利用ArcGIS 10.2软件进行空间叠加分析, 计算城乡建设用地可承载临界值和可承载城乡建设用地丰度。

(4) 空间匹配度分析和相关分析：计算城乡建设用地空间匹配度, 并将城乡建设用地可承载临界值、可承载城乡建设用地丰度和空间匹配度进行相关分析, 研究现状利用状态与其合理开发状态之间的关系。

2.2 南昌市城乡建设用地开发适宜性评价

城乡建设用地开发适宜性评价是指评定一定条件下土地开发为建设用地的适宜性程度^[35]。根据数据的可获取性和南昌市实际, 从生态安全、自然条件、社会经济等3方面选取11个评价指标, 构建城乡建设用地开发适宜性评价指标体系, 并将指标分为强限制性指标和较强限制性指标两类, 如表 1所示。

在评价时, 考虑到经济发展与生态保护的矛盾, 将生态红线作为强限制性指标处理。生态红线是指为维持生态安全和生态系统完整性, 通过科学合理的方法划定的必须实施特殊保护的国土空间界线^[36]。同时为贯彻最严格的耕地保护政策, 选取永久基本农田区为强限制性指标。地面塌陷区开发难度大且多为自然灾害易发区, 所以将地面塌陷区也列为强限制性指标。此外, 坡度直接影响建设用地开发成本和开发难度, 坡度过大会影响建筑物的稳定性及增加潜在危险性。相关研究表明, 在进行建设用地开发时, 坡度大于25°的陡坡地一般被视为“禁止建设区”^[37-38]。地质灾害直接威胁到人类的生命财产安全, 区位条件和人口密度影响社会经济发展以及信息的可达性和居民生活的便利程度。因此, 在适宜性评价时, 都将其作为评价指标考虑在内。

对于强限制性指标, 根据“一票否决”原则进行0和1赋值。对于较强限制性指标, 在参考相关研究^[39-41]的基础上, 根据适宜性程度进行0—100赋值, 采用层次分析法和专家打分法确定各较强限制性指标的权重。由于坡度对建设用地开发影响特殊, 对于坡度大于25°的陡坡地进行0和1赋值, 而坡度小于25°的地区则按限制性等级进行0—100赋值(见表 1)。

采用限制系数评价模型, 计算各评价单元的城乡建设用地开发适宜性分值, 计算公式如下：

(1)

式中E为综合适宜性分值；i为强限制性指标序号；j为较强限制性指标序号；F_i为第i个强限制性指标适宜性分值；f_j为第j个较强限制性指标适宜性分值；w_j为第j个较强限制性指标权重；m为强限制性指标个数；n为较强限制性指标个数。

2.3 南昌市城乡建设用地承载力计算

城乡建设用地可承载临界值是指区域土地资源所能承载的人类活动的最大规模和强度。基于城乡建设用地开发适宜性评价结果, 将最适宜开发建设区视为最适宜承载区, 较适宜开发建设区视为较适宜承载区, 结合现状建设用地情况, 通过ArcGIS 10.2软件, 测算城乡建设用地可承载临界值。计算公式如下：

(2)

式中, LDI为城乡建设用地可承载临界值, E₁为最适宜承载区面积, E₂为较适宜承载区面积, C为现状城乡建设用地面积, U为研究区总面积。

城乡建设用地已开发强度反映区域现状建设用地利用状况, 以各县(区)现状城乡建设用地总量(C)占各县(区)区域总面积(U_T)的比例表示。

(3)

以城乡建设用地可承载临界值为城乡建设用地极限开发强度, 将其与已开发强度相减, 差值为剩余开发强度, 即城乡建设用地的开发潜力。

2.4 南昌市可承载城乡建设用地丰度计算

可承载城乡建设用地丰度(A)反映的是研究区内可承载建设用地的丰富程度^[42], 主要从数量和空间分布情况考虑。数量方面(A₁), 以各县(区)可承载城乡建设用地面积(E₁+E₂)占县域总面积(U_T)的比例表示；空间方面(A₂), 以各县(区)可承载城乡建设用地面积(E₁+E₂)占研究区可承载城乡建设用地面积(U_S)的比例表示, 数据标准化处理后求和得到各县(区)可承载城乡建设用地丰度, 计算公式如下：

(4)

采用自然断点法将各县(区)的丰度差异划分为Ⅰ级、Ⅱ级和Ⅲ级^[43-44]。其中, Ⅰ级城乡建设用地丰度为“丰富”, A>1.0；Ⅱ级城乡建设用地丰度为“中等”, 0.6 < A < 1.0；Ⅲ级城乡建设用地丰度为“缺乏”, A < 0.6。

2.5 南昌市城乡建设用地空间匹配度分析

城乡建设用地空间匹配度是指现状城乡建设用地与可承载城乡建设用地的空间匹配情况。将现状城乡建设用地与最适宜承载区和较适宜承载区进行空间叠加, 分析现状城乡建设用地的可承载适宜性情况, 计算公式如下：

(5)

式中, M为现状城乡建设用地与可承载城乡建设用地的空间匹配度, C₁为现状城乡建设用地处在最适宜承载区面积, C₂为现状城乡建设用地处在较适宜承载区面积, C为现状城乡建设用地总面积。

3 结果与分析 3.1 城乡建设用地适宜性评价分析

通过ArcGIS 10.2软件的处理和SPSS 22.0软件的聚类分析, 南昌市城乡建设用地开发适宜性评价结果见表 2和图 3, 最适宜区、较适宜区、较不适宜区和不适宜区面积分别为1466.19、1346.01、353.69 km²和3770.70 km²。其中, 不适宜区面积最大, 占研究区面积的54.36%；较不适宜区面积最小, 占研究区面积的5.10%；最适宜区和较适宜区分别占比为21.14%和19.40%。

由图 3可知, 南昌市城乡建设用地开发适宜性空间分异显著, 不同等级的空间集散特征差异较大。最适宜区分布较为集中, 主要分布在中心城区的西湖区、青云谱区、青山湖区、东湖区等, 由于这些地区地势平坦、交通便利、经济发达且较少处于生态保护区域, 适宜进行建设用地开发。较适宜区分散分布在适宜区周边地区, 如新建县、南昌县等, 在对较适宜区进行建设用地开发时, 应根据实际情况和生态保护要求, 避免过度开发。较不适宜区面积小且分布较为分散。由于受到永久基本农田区和生态保护的限制, 研究区大部分地区不适宜开发为建设用地, 作为面积占比最大的适宜性等级, 不适宜区广泛分布在各县(区)。

3.2 城乡建设用地承载力分析 3.2.1 城乡建设用地开发强度分析

以城乡建设用地可承载临界值为极限开发强度, 以现状建设用地占区域总面积比例为已开发强度, 以极限开发强度与已开发强度差值为剩余开发强度, 测算研究区城乡建设用地开发强度, 结果如表 3所示。

南昌市城乡建设用地极限开发强度为41.25%, 各县(区)的极限开发强度介于19.43%—98.71%之间, 差异较大。其中, 东湖区、青山湖区和青云谱区的极限开发强度大于90%, 这些区域地理条件优越, 交通便利, 为南昌市的经济、政治、文化中心, 大部分适宜开发为建设用地；西湖区的极限开发强度为86.90%, 也比较适宜进行建设用地开发；南昌县的极限开发强度为50.50%；新建县、进贤县和安义县等极限开发强度在30%—40%之间。湾里区依托梅岭国家森林公园, 坚持生态立区、旅游兴区, 实施严格的生态保护制度, 因此湾里区大部分地区不宜开发为建设用地, 其极限开发强度为全市最小值19.43%。

南昌市城乡建设用地已开发强度为17.20%, 由于区域地理条件和经济状况差别较大, 各县(区)的已开发强度差异明显, 范围在10.41%—86.70%。由于地理优势明显, 位于中心城区的市区, 其已开发强度大于其他县(区)。其中, 青云谱区作为南昌市“银三心”繁荣板块之一, 其已开发强度最大, 为86.70%；西湖区和青山湖区分别作为南昌市主要历史文化区和工业区, 已开发强度分别为76.19%和72.37%；东湖区、南昌县已开发强度适中, 分别为43.51%和20.19%。处于生态保护限制区和农业发展区的县(区)已开发强度一般较小, 安义县、新建县、进贤县和湾里区的已开发强度均小于20%, 其中进贤县的已开发强度为最小。

将城乡建设用地极限开发强度与已开发强度相减, 得到南昌市的剩余开发强度为24.05%, 略大于现状城乡建设用地规模, 表明南昌市还有较大的开发空间。各县(区)的剩余开发强度介于5.65%—49.83%之间, 差异较大。其中, 湾里区的剩余开发强度最小, 为5.65%, 由于受到永久基本农田区和生态保护的限制, 不宜进行建设用地开发。西湖区和青云谱区的剩余开发强度介于10%—20%之间, 由于已开发强度大导致其开发潜力较小。大部分县(区)的剩余开发强度范围在20%—30%, 包括青山湖区、新建县、安义县和进贤县。安义县、进贤县和新建县虽然受到永久基本农田区和生态保护的限制, 但由于已开发强度较小, 仍有部分地区可进行建设用地开发。青山湖区基于发展经济的需求, 大部分面积已经开发为建设用地, 剩余开发强度相对较小。剩余开发强度大于30%的县(区)有南昌县和东湖区, 其中东湖区的剩余开发强度最大, 为49.83%, 在进行建设用地规划时, 在不破坏区域生态环境的情况下可优先考虑指标划分。

3.2.2 可承载城乡建设用地丰度分析

南昌市可承载城乡建设用地丰度情况整体良好, 但各县(区)差异较大, 丰度范围介于0.19—1.50之间, 丰度最大的是南昌县, 最小的是湾里区, 具体如图 4所示。可承载城乡建设用地丰度为“丰富”(即Ⅰ级)的县(区)有新建县、青山湖区、青云谱区、南昌县和东湖区；丰度为“中等”(即Ⅱ级)的县(区)有进贤县和西湖区；丰度为“缺乏”(即Ⅲ级)的县(区)有安义县、湾里区。

3.3 城乡建设用地空间匹配度分析

将现状城乡建设用地与适宜性评价结果进行空间叠加, 得到现状城乡建设用地的可承载适宜性情况, 结果如图 5所示。经统计, 南昌市现状城乡建设用地位于最适宜、较适宜、较不适宜和不适宜承载区的面积分别为963.31、120.14、24.67 km²和85.46 km², 分别占现状城乡建设用地面积的80.70%、10.67%、2.07%和7.16%。通过计算各县(区)现状城乡建设用地位于最适宜和较适宜承载区的面积占区域总面积的比例, 得到各县(区)的空间匹配度在67.41%—99.25%之间, 反映各县(区)的空间匹配度较好, 且差异不大。其中, 空间匹配度大于90%的县(区)有东湖区、青云谱区、西湖区和青山湖区, 均为中心城区, 区位条件好且地势平坦, 大部分地区适宜进行建设用地开发, 因此现状开发适宜性情况相对较好；安义县、进贤县、南昌县和新建县的空间匹配度在70%—90%, 空间匹配度较好；受地理位置及生态环境的限制, 湾里区的空间匹配度最小。南昌市的空间匹配度为91.37%, 说明研究区现状城乡建设用地的承载适宜性整体情况良好。

将南昌市各县(区)的城乡建设用地可承载临界值、可承载城乡建设用地开发丰度与城乡建设用地空间匹配度进行相关性分析, 发现空间匹配度与城乡建设用地可承载临界值和可承载城乡建设用地丰度均呈正相关关系(表 4), 反映可承载临界值和丰度越高, 其空间匹配度越高, 这说明城乡建设用地的现状利用状态与其合理开发状态存在一定的联系, 揭示两者之间的空间匹配程度。

4 结论与讨论 4.1 讨论

随着经济快速发展和城镇化进程不断加快, 土地开发强度日益增大, 导致生态环境不断恶化, 分析城乡建设用地承载能力及空间格局匹配情况, 对解决土地利用不合理、生态破坏等问题具有重要意义。本研究以江西省南昌市为例, 探讨了城乡建设用地承载力和现状开发状态之间的空间匹配关系, 为南昌市国土空间优化和可持续发展提供参考。

在适宜性评价中, 增加了生态红线、地质灾害、区位条件、人口密度等评价指标, 考虑比较全面, 但通过层次分析法和专家打分法确定指标权重, 存在一定的主观性, 在今后的研究中有待进一步改进。从整体来看, 南昌市可承载城乡建设用地丰度和空间匹配情况良好, 各县(区)城乡建设用地开发强度差异显著, 安义县、进贤县的极限开发强度和已开发强度小于区位条件好的青山湖区、青云谱区, 剩余开发强度与区位条件无明显关联, 这与易丹等^[39]针对环鄱阳湖城市群的研究结果较为吻合。

通过对城乡建设用地承载力和现状开发状态之间的空间匹配关系分析发现, 南昌市及各县(区)的现状开发适宜性情况较为合理, 但空间分异显著, 应根据国土空间规划相关要求, 对土地资源进行合理利用和优化配置。湾里区由于大部分位于不适宜开发区且开发潜力较小, 应控制建设用地开发, 加强生态环境保护；安义县、新建县和进贤县等虽然也有较大部分位于不适宜开发区, 但其已开发强度小, 开发潜力大, 应结合地区实际情况, 在不破坏生态环境和保障农业生产的基础上, 对其适宜建设区进行适量开发, 承接中心城区的部分功能, 带动经济发展；对于开发潜力较大、空间匹配度较高的东湖区和南昌县, 由于位于都市发展核心区, 可充分利用区位优势, 在不破坏生态环境的基础上可优先安排城乡建设用地开发；对于已开发强度大、开发潜力小的青山湖区、青云谱区和西湖区, 为缓解土地承载压力, 不宜进行城乡建设用地开发, 同时应加强城市改造和生态建设, 改善人居生态环境。

4.2 结论

研究通过对南昌市进行城乡建设用地开发适宜性评价, 根据评价结果确定可承载区域, 计算城乡建设用地可承载临界值和可承载城乡建设用地丰度, 并与现状城乡建设用地进行空间匹配度分析, 探讨城乡建设用地现状开发状态与其合理开发状态之间的联系。主要结论如下：

(1) 根据城乡建设用地适宜性评价结果, 南昌市不适宜区面积最大且广泛分布在各县(区), 较不适宜区面积最小且分布分散, 最适宜区和较适宜区面积共为2812.20 km², 占总面积的40.54%, 集中分布在地理条件良好的中部市区。

(2) 南昌市城乡建设用地可承载临界值为41.25%, 各县(区)可承载临界值在19.43%—98.71%之间, 各县(区)可承载城乡建设用地丰度范围为0.19—1.50, 各县(区)值均有较大差异。其中, 东湖区、青云谱区、青山湖区由于经济发达、地理条件优越, 可承载临界值和可承载丰度均大；而湾里区由于生态保护的限制, 在进行建设用地开发时要注意与生态环境的保护相结合。

(3) 南昌市剩余开发强度为24.05%, 略大于现状建设用地规模, 表明南昌市还有较大的开发潜力。各县(区)现状城乡建设用地与可承载城乡建设用地的空间匹配度介于67.41%—99.25%之间, 说明城乡建设用地现状开发状态较好。对现状适宜性与开发状态相关性研究表明, 在进行后续城乡建设用地开发时, 应注意与可承载临界值和可承载丰度相结合, 合理规划国土空间布局, 避免过度开发。