快速城市化背景下的城市生态安全问题已成为近些年规划领域关注的热点。识别城市中重要的生态用地, 构建生态安全格局已成为维护城市生态系统平衡, 促进城市可持续发展的重要途径^[1-3]。就定量化的构建方法而言, 主要是利用生态学理论的景观格局优化模型, 例如可通过最小累积阻力模型, 形成生态用地空间运动趋势的阻力面, 以此来判别研究区域生态安全格局^[4-5]。而随着人地矛盾的加剧, 相关学者也在研究作为自然与人工耦合而成的复合生态系统的城市, 应如何协调生态保护与建设用地扩张之间的关系。例如, 以不同水平的生态安全格局用地范围作为城市空间管制分区的用地划分依据^[6-7]; 或以生态安全格局为基础, 在空间上确定城市扩展的合理范围与发展方向, 提出建设用地最优扩展方案^[8]。

本文的研究对象为山水型城市。该类城市具有得天独厚丰富的山水资源, 且城市形态结构较大程度地受自然地形的制约。近些年, 由于盲目开发, 该类城市“摊大饼”式扩张导致其山体被蚕食, 水域面积缩小, 使其失去了原有“山”“水”相连的固有格局, 而城市内部空间使用效率低下。因此, 构建山水型城市生态安全格局, 探讨生态安全目标导向下建设用地开发适宜性就更具现实意义。安庆市规划区三面环水, 一面临山, 属典型的山水城市, 但同时也是长江中游城市群主要成员之一, 由于近些年工业用地扩张, 使城市滨江滨湖岸线70%被占据, 因石化而使城市原有生态资源遭到严重破坏(图 1)。在此背景下, 本文结合山水型城市的主要特征, 通过重要生态用地识别, 来构建符合城市特色风貌的生态安全格局, 在此基础上结合城市空间增长需求, 以乡镇为基本评价单元, 进行可建设用地开发适宜性分区, 并提出相应的开发策略。

1 研究区概况与研究方法 1.1 研究区概况

安庆市位于长江下游上段北岸, 安徽省西南部, 支流甚为发育。本次研究的范围为安庆市规划区, 是城市发展核心区域, 也是城市建设与生态保护矛盾突出区域, 总面积886平方公里, 包括迎江、大观、宜秀3个行政区及皖河农场(图 2)。研究区西北部为低山丘陵、岗地和阶地, 东南为长江冲积平原, 其间湖泊广布。中山、低山、丘陵、岗地及平原地貌俱全, 为典型的山水型城市。城市防洪、排涝、自然保护区、风景名胜区等多种限制性因素对城镇空间利用方式影响很大, 城市建设用地发展受到限制。因此如何利用山水型城市特有的生态资源, 构建符合城市特色风貌的生态安全格局, 并以此来提出城市可建设用地开发策略就显得尤为重要。

1.2 数据来源

本研究数据源主要为安庆市2015年航空影像图(空间分辨率15m)、安庆市规划区土地利用总体规划(2006-2020)、规划区地形图(1:10000)、以及安庆市水源保护区、自然保护区及工程地质分布图等。基于上述数据资料并结合实地调研, 将研究区划分为城市建设用地、道路用地、村庄建设用地、林地、耕地、园地、湿地、河流湖泊、水库、坑塘沟渠、草地、其他建设用地等类型(图 3)。社会经济评价指标主要数据来源为《安庆市统计年鉴》。

1.3 研究方法

首先为重要生态用地即“源”的识别, 利用GIS空间分析技术, 结合研究区“山、湖、城、江”现状自然生态特征, 选择识别生物多样性保护、水资源安全、地质灾害规避3类生态用地, 并综合叠加上述分析结果, 将研究区生态用地划分为极重要、较重要、一般重要3个级别; 最后, 将极重要生态用地与相关法规、标准及政策所规定的禁建区作为“源”, 利用最小阻力模型构建不同水平生态用地安全格局^[9]。在此基础上, 结合城市空间扩展动力因素, 以各乡镇为评价单元通过生态约束和发展需求双重评价进行建设用地开发分区, 并提出相应开发策略。

1.3.1 单因子的生态用地识别

(1)生物多样性保护

该类生态用地的划分依据为生物多样性, 通过谢高地等人制定的生物多样性服务当量即单位面积内每种土地利用类型所能提供生物多样性保护相对贡献大小的潜在能力^[10-11]予以间接表征。依据研究区不同土地利用类型计算出各评价单元生物多样性服务价值, 并可通过保护级别予以修正。计算式如下：

式中, S为生物多样性服务价值, l为土地利用, n为生物多样性服务当量, m为保护级别修正值, 其中大龙山森林公园、高程大于120m林地、重要湿地保护区等区域赋值为1.5, 高程50-120m林地赋值为1.25, 其他区域为1。通过GIS空间分析计算得到研究区各土地利用类型的S值, 并运用自然断裂法对该类生态用地进行敏感性分区。

(2)水资源安全

结合研究区现状水资源特征, 选取水源保护、洪水淹没区范围、地表重要水源安全、洪水调蓄区作为水资源安全重要评价因子^[12]。其划分标准(表 1)主要为参考相关研究成果得到。水源保护区及洪水淹没区范围主要依据为安庆市各类专项规划、《安庆市城市总体规划(2010-2030)》、《安庆市城市给水工程规划(2012-2030)》。地表水源安全与洪水调蓄区类型分别是利用ArcGIS 10.0缓冲区分析与水文分析工具获得。最终将4个指标叠加, 基于“木桶”原理进行水资源安全敏感性分区。

(3)地质灾害规避

崩塌、泥石流、采空塌陷、滑坡为研究区地质灾害的主要类型。而其主要受到坡度、地形起伏度、高程及人类活动强度等的影响。基于前人研究成果^[8], 得到地质灾害敏感性划分标准(表 2), 并将敏感值以150、120、100作为其划分的断裂点^[9]。在此基础上, 叠加已有的地质灾害易发区, 最终获得研究区地质灾害规避敏感性等级划分范围。

1.3.2 综合生态用地的识别

本研究运用GIS加权总和工具将生物多样性保护、水资源安全、地质灾害规避三类生态用地分析结果进行等权重叠加, 即认为上述3个生态过程对于维护研究区生态安全具有同等作用。将叠加的分析结果以自然断裂法进行分类, 最终将研究区生态用地分为极重要、较重要、一般重要3个级别。

1.3.3 生态用地安全格局构建

(1) “源”的确定

“源”从物种保护的角度来看, 应为保护物种生存的栖息地^[13]。本研究认为“源”是生态保护核心区域, 因此应将上述分析所确定的极重要生态用地与相关法规、标准及政策所规定的禁建区用地叠加得出的最大用地范围作为“源”。

(2)阻力面的建立

研究利用最小累积阻力模型(Minimum Cumulative Resistance, 简称MCR)来建立反映生态用地空间运动趋势阻力面^[14-16]。该模型需要3方面因素, 即源、距离和景观介面特征, 计算式如下：

该公式是根据Knaapen等人的模型及地理信息系统中常用的费用距离模块修改而来^[17]。MCR为物种由出发点到任意一目标点之间所耗费的最小累积阻力。式中, f是未知函数, 反映空间中任一点的MCR值与(D_ij×R_i)正相关关系。D_ij是物种从源j到空间某目标点所穿越的景观基面i的运动距离; R_i为景观i对某物种运动的阻力。(D_ij×R_i)可被看作是物种从源到空间目标点的某一路径的相对易达性的衡量。因此, 阻力面反映了生态用地水平运动的潜在可能性及趋势^{[4, 18]}。

(3)依据阻力面判别研究区生态安全格局

依据生成的阻力面, 可判别除源之外的生态安全格局的其他组成部分, 主要包括^[19-20]：

缓冲区其为MCR阻力面上源周围的低阻力区, 相关研究表明依据阻力面得到的MCR值与面积的关系曲线存在着某些阶段性的“门槛值”^{[4, 9, 19]}, 研究以其门槛值作为不同安全水平生态缓冲区划分依据。

“源”间廊道即任意“源”之间低累积阻力谷线, 为连通相邻“源”最便捷、最高效生态廊道。

战略点是联系相邻“源”之间有关键意义的节点, 在控制生态过程方面起到重要的作用。通过增加战略点, 可提高区域生态用地之间的连通性。

1.3.4 基于生态安全格局的可建设用地开发分区

在得到生态用地基础上, 基于可建设用地生态影响及开发潜力评价, 进行建设用地开发适宜性分区(图 4)。首先基于生态约束的角度进行生态影响分区, 即选取生态底线用地与易损性用地占比来衡量生态环境影响指数水平。其比值越高则生态底线用地规模越大, 相应生态环境影响力就越大。与此同时进行经济发展需求分析, 通过多元统计因子分析法, 得到各评价单元的经济发展需求能力排名。最终通过综合发展分区分级矩阵将建设用地分为优先发展区、适度发展区、控制发展区以及禁止发展区。

2 结果与分析 2.1 综合生态用地识别分析 2.1.1 基于单因子的生态用地识别

生物多样性保护生态用地识别中(表 3和图 5), 高敏感区域面积为201.7km², 占规划区总面积22.8%, 主要分布于大龙山森林公园等森林茂密, 野生动植物种类繁多区域; 水资源安全生态用地识别中(表 3和图 5), 高敏感区域面积268.5km², 所占比例为30.3%, 分布于长江流域、石塘湖等重要的地表水源涵养区; 地质灾害规避生态用地识别方面(表 3和图 5), 高敏感区面积为152.8km², 所占比例为17.2%, 此区域是崩塌与山体滑坡等地质灾害易发区。

2.1.2 多因子综合评价的生态用地分级

基于自然断裂法将GIS叠加分析的结果进行分类(表 3和图 5), 得出极重要即底线型生态用地范围, 其为安庆市规划区生态保护核心区域, 面积394.1km², 所占比例为44.5%, 需加以严格保护, 禁止开发建设。

2.2 生态用地安全格局的构建 2.2.1 景观生态流阻力面构建

依据相关研究, 各生态保护源地与土地覆被类型越接近, 其运动所耗费的阻力值就越小, 反之则越大。本文将极重要生态用地与相关法规、标准及政策所规定的禁建区作为“源”(表 4), 以土地覆被类型为阻力因子, 参考前人研究成果获得不同土地覆被类型的阻力系数值^{[2, 9]}(表 5)。在此基础上运用ArcGIS 10.0距离分析模块, 计算出各生态保护源地与各土地覆盖类型之间的费用距离, 最终获得生态保护源地最小累积阻力面(图 6)。

2.2.2 生态用地安全格局构建

在获得阻力面的基础上, 参考GIS水文分析原理, 根据研究区现状景观特征, 运用GIS中的水文分析模块确定连接相邻“源”之间的最小耗费路径, 进而获得生态廊道的空间分布。通过提取“源”间生态流水平运动的最大耗费路径和最小耗费路径的交汇点, 并结合现状特征, 获得“战略点”的位置。最终获得安庆市高、中、低3种水平的综合生态安全格局(图 7)。“源”与低水平生态安全用地总面积为582.9km², 是保障城市生态系统服务的核心区域, 占规划区总面积的65.8%, 此区域应加强对原有生态资源的保护, 严格禁止进行与生态保护无关的开发建设活动。中、高水平生态安全用地总面积为163.1km², 所占比例为18.4%, 该区域可作为未来安庆市规划区生态基础设施培育区, 同时也是城市严格的限制建设区, 应严格控制城市在该区域的开发强度与用地开发性质。分析结果表明依山临水、组团式城市用地布局是安庆市基于生态安全格局的、合理的发展模式, 大龙山镇、老峰镇、海口镇为未来城市空间扩展的主要区域。除去“源”与低水平的生态安全用地, 城市可建设用地为303.1km², 满足2030年城市总体规划用地规模需求。

2.3 基于生态安全格局的建设用地开发策略

传统的建设用地选择主要依据经济发展需求, 本文先构建生态安全格局, 确定哪些生态用地是不可开发的, 再结合经济发展需求确定建设用地开发分区。选取乡镇为基本评价单元, 运用生态约束和经济发展需求2类指标评价^[21]。经济发展需求选择的主成分因子分别为城镇建设与经济发展因子、经济效率与人口集聚因子、对外交通因子(表 6)。依据生态约束和发展需求双重评价(表 7), 获得安庆市规划区综合发展分区分级矩阵(表 8)。

基于此将研究区分为优先发展区, 包括老峰镇、大龙山镇、白泽湖乡、长风乡、海口镇, 此区域经济发展需求强或较强、生态环境影响力中等或较低, 是城市化发展重点区域, 提倡高强度开发, 可进行较大规模工业开发; 适度发展区包括杨桥镇, 该区域经济发展需求较高、生态环境影响力较高, 应适度控制开发强度, 尤其是控制建设用地扩展规模和速度, 避免过度开发, 适当发展无污染工业项目; 控制发展区包括皖河农场、五横乡、山口乡、罗岭镇, 该区域经济发展需求一般或低, 生态环境影响力高或较高, 因此应控制开发强度和规模, 有节制地开发休闲、旅游、度假等同时具有自然生态保护意义与经济开发效益的绿色产业, 禁止一般工业项目建设; 禁止发展区包括新洲乡, 此区域经济发展需求低、生态环境影响力高(图 8)。

3 结论与讨论

(1)构建山水型城市高、中、低3种水平生态安全格局, 可有效协调该类城市空间扩展与生态要素限制之间的矛盾。除去“源”与低水平的生态安全用地, 城市可建设用地为303.1km², 即能维护该类城市固有的“山”“水”景观格局, 延续和发展研究区原有山水特色风貌, 也可满足2030年城市总体规划用地规模需求。

(2)通过生态约束和发展需求双重评价, 开展可建设用地适宜性分区, 弥补了传统可建设用地选择仅考虑经济发展、忽略生态保护的弊端, 评价结果具有较强的科学性。同时, 评价单元以乡镇行政管理范围为依据可提高可建设用地管控的可操作性。

(3)本研究在构建研究区生态安全格局基础上进行可建设用地开发策略分析, 而基于土地生态服务的城市土地分类体系应不仅仅包括生态用地与建设用地, 而应存在一些缓冲带即生态交错地, 包括城乡结合部土地、农业用地等, 对于该类用地应当如何进行开发、管理与保护还应进一步完善。与此同时, 出于对山水型城市生态要素重要性、研究尺度等方面的考虑, 生态用地识别仅选取生物多样性保护、水资源安全、地质灾害规避三类生态因子, 没有将气候调节、水土保持等因素一起纳入指标体系, 是否具有全面性与代表性还有待深入探讨, 且单一生态过程的影响因子体系和经济发展需求能力评价因子的选择、单一生态过程敏感性和生态安全格局分级所涉及的“门槛值”与“断点法”由于缺乏相关标准难免会有遗漏, 这都需要在今后的研究工作中不断完善。