干旱区城市生态安全格局优化研究——以内蒙古呼和浩特市为例

冯琰玮, 甄江红, 马晨阳

(内蒙古师范大学 地理科学学院, 呼和浩特 010022)

快速城镇化背景下的土地利用变化对生态环境的负面效应不断加剧，这正深刻改变着区域自然环境基底，具体表现在生态用地规模的减少[1]、耕地面源污染加剧[2]、环境承载力水平下降[3-4]、生境适宜性降低[5]等方面。与此同时，区域生态安全水平对城市用地扩展的约束作用也进一步加强，科学合理的生态安全格局可以有效缓解经济发展和环境保护之间的矛盾。构建科学合理的生态安全格局是保障区域生态安全的重要举措。规划先行，让绿水青山变成金山银山的顶层设计也成为城市地理学所关注的热点话题[6]。

国外对生态安全格局的研究主要集中在城市增长边界的划定[7]、生态网络的构建[8]等角度进行。国内学者对武汉市[9]、天水市[10]、关中城市群[11]等地的生态安全格局构建进行了诸多有益探索。从研究内容来看，多数研究单从生态用地扩张的角度构建区域生态安全格局，对于城市这类经济、人口密度高度集中的区域，如何加大对生态用地的保护成为关键。对于地处西北内陆的呼和浩特市来说，涉及其区域生态安全的研究主要集中在城市用地扩展的生态环境效应及城市增长边界研究[12-13]，有关研究区生态安全格局的报道还鲜见纸端。鉴于此，本文试构建建设用地与生态用地扩展阻力评价指标体系，通过MCR模型建立区域生态安全格局，识别潜在生态廊道与生态战略点，构建区域生态网络框架，对城市区域空间管制进行分区，以期为生态脆弱区的生态安全格局优化研究进行有益补充。

1 研究区概况与数据来源

1.1 研究区概况

呼和浩特市地处内蒙古自治区的中西部，其气候属典型中温带大陆性季风气候，年均气温6.4℃，年均降水量335.2～534.6 mm。2016年市区人均水资源占有量仅为465 m3/年，是全国平均水平的1/6，水资源供需矛盾突出，属生态严重失衡的地区之一。本文以呼和浩特市4个市辖区(赛罕区、新城区、回民区、玉泉区)及下属6个乡镇(榆林镇、金河镇、黄合少镇、保合少镇、攸攸板镇、小黑河镇)为研究区，总面积为2 083.68 km2。2016年末，呼和浩特市建成区面积达到260 km2，人均国民生产总值为10.32万元，全市常住人口城镇化率达69.08%，均高于全国与自治区平均水平。高强度的经济开发活动使得研究区生态问题愈发凸显，以呼和浩特市为研究区，对探索生态脆弱区生态安全格局的优化也具有一定的代表性和典型性。

1.2 数据来源

采用2016年8月29日Landsat-8 OLI遥感影像作为土地利用类型信息提取的基础数据；DEM数字高程数据选择GDEMV2，分辨率为30 m的数据产品。遥感数据均来源地理空间数据云(http:∥www.gscloud.cn/)。净初级生产力NPP数据来源于北纬18°以北中国陆地生态系统逐月净初级生产力1 km栅格数据集[14]。采用人机交互目视解译法对研究区土地利用类型进行分类，参照《土地利用现状分类》国家标准(GB/T21010～2017)，分为林地、耕地、草地、建设用地、水域和未利用地6大类型。其他空间数据提取主要运用ArcGIS，Fragstats，GeoDa等软件进行处理。

2 研究方法

2.1 城市用地与生态用地扩张阻力评价指标体系的建立

建立简明科学的用地扩张阻力评价指标体系对区域发展规划具有重要意义。鉴于研究区生态环境脆弱的实际，在城市发展过程中强调对生态、资源环境的保护也有利于推动区域人地关系和谐发展。考虑到城市环境的自然特征与人工特性，从研究区的自然环境、资源条件、人类干扰实际出发，选取地形坡度、土地覆盖、生物敏感性、生态风险性、区域交通5类影响因子构建城市用地“源”与生态用地“源”扩张阻力评价指标体系，包含二级阻力因子15个。选择面积大于5 km2的建设用地斑块作为城市用地扩张“源”，总面积为390.29 km2；选择面积大于2 km2的林地和草地斑块[15]作为生态用地扩张“源”，总面积为479.35 km2。将阻力因子分为5级进行赋值，具体赋值从大到小依次为9，7，5，3，1。各因子权重的选择应用Delphi法进行赋值，具体评价指标及赋值情况见表1。对各阻力因子的解释如下：

(1) 地形坡度。a.高程：海拔越高，生态环境越脆弱，就越不利于城市用地扩张。采用自然断裂法以1 141，1 300，1 511，1 757 m为界对城市用地扩张进行阻力赋值，生态用地扩张阻力赋值与之相反；b.坡度：坡度越高，越不利于建设用地扩张。根据《城市用地竖向规划规范》，采用自然断裂法以3°，8°，15°，25°为界对城市用地扩张进行阻力赋值，生态用地扩张阻力赋值与之相反；c.地表崎岖度：地表越粗糙，越不适宜开发建设活动。采用自然断裂法以1.03°，1.09°，1.19°，1.37°为界对城市用地扩张进行阻力赋值，生态用地扩张阻力赋值与之相反。(2) 土地覆盖。d.土地利用类型：源地土地利用类型与当前土地利用类型越接近，就越利于源地用地扩张。根据土地利用类型对建设用地与生态用地扩张进行阻力赋值；e.土壤稳定性：土壤越稳定，就越有利于城市用地与生态用地扩张。建设用地与生态用地扩张阻力赋值一致。(3) 生物敏感性。鉴于研究区位于干旱区与半干旱区的过渡地带，降水少而变率大，植被覆盖偏低。采用f.植被覆盖率(NDVI)、g.环境检测指数(GEMI)及h.净初级生产力(NPP)表征研究区植被的固碳释氧能力。植被分布越宽泛的地区越适宜生态用地扩张。采用自然断裂法对NDVI及GEMI以0.15，0.35，0.55，0.75为界，NPP以12.8，26.9，36.5，47.7为界向生态用地扩张进行阻力赋值，建设用地扩张阻力赋值与之相反。其中，NDVI=(NIR-R)/(NIR+R)；GEMI=η(1～0.25η)-(R-0.125)/(1-R)。其中，NIR，R分别为Landsat遥感影像中的近红、红波段，η=[2(NIR2-R2)+1.5NIR+0.5R]/(NIR+R+0.5)。(4) 生态风险性。i.土地退化率：土地退化越严重，就越不适宜各类用地扩张。建设用地扩张阻力赋值与生态用地阻力赋值一致。j.景观生态风险：采用景观干扰度表征研究区的景观生态风险程度。将景观生态风险指数赋予1.5 km×1.5 km格网单元后进行克里金插值。结果表明城市景观生态风险的高值主要集中在城市建设用地边缘地带，景观生态风险值越高，越说明该地带是城市用地扩张的潜力区域。采用自然断裂法以0.35，0.27，0.2，0.13为界对城市用地扩张进行阻力赋值，生态用地扩张阻力赋值与之相反。(5) 区域交通。城市用地扩展多沿交通线进行定向推进，具有重要交通意义的地带往往具有很高的经济性。城市交通的发展对城市结构的内部调整也具有战略性意义。本文对研究区的k.水体、l.风景名胜区、m.城市中心、n.高速路口、o.铁路站点等矢量数据进行缓冲区分析。其中，水体以200，400，600，800 m做缓冲区对其赋值；风景名胜区以2 500，2 000，1 500，1 000 m做缓冲区对其赋值；城市中心以1 000，1 500，2 000，2 500 m做缓冲区对其赋值；高速路口以300，600，900，1 200 m做缓冲区对其赋值；铁路站点以500，1 000，1 500，2 000 m做缓冲区对其赋值。

表1 城市用地“源”与生态用地“源”扩张阻力指标评价体系

2.2 最小累积阻力面模型(MCR)的构建

最小累积阻力面模型(minimal Cumulative Resistance，MCR)由荷兰学者Knappen提出，最初始于对物种扩散的研究。近年来国内学者将其广泛应用于城市地理学、景观生态学等领域，其实质是各“源”地经过不同阻力景观所耗费的费用或克服阻力所做的功[16]。采用俞孔坚修正的MCR模型计算得到城市用地扩张阻力面和生态用地扩张阻力面[17]。其中，MCR模型的计算方式如下：

(1)

式中：MCR表示源扩展到某景观处的最小累积阻力值；f是个未知的负函数；Dij表示源i到景观j的空间距离；Ri表示源i在扩展过程中所受到的阻力值大小。本文将生态用地扩张阻力面与建设用地扩张阻力面进行空间叠加，从而构建研究区生态安全格局。具体计算方式为：MCR生态安全=MCR生态用地扩张-MCR建设用地扩张，将二者的差值作为生态安全格局的建立依据。当MCR=0时，是适宜生态用地与建设用地扩张的分界线；当MCR>0时，适宜城市用地扩张；当MCR<0时，适宜生态用地扩张。

2.3 生态廊道的识别

生态廊道是区域生态安全格局的重要组成部分，作为区域生态系统能量流动与物质循环的重要载体，对于改善区域生境质量、完善生态系统功能、加强生物多样性保护具有重要意义。识别重要生态廊道是维护区域生态要素流动的有效手段，构建区域生态安全网络利于增强生态系统功能的完整性。以生态源地的几何中心点作为生态源点，从每个源点出发，运用成本路径模块构建迭代模型提取源点(n)到剩余n-1个中心点的最小阻力路径，即潜在生态廊道。大型斑块和较宽的廊道生境适宜性比较髙，物种迁移的阻力较小，能大大增加物种在迁移过程中的存活率和成功率[18]。而潜在廊道的重要性可以用源地斑块间的相互作用强度加以判别。本文基于重力模型(Gravity Method)相互作用原理，定量评价斑块之间相互作用力的强弱，引入相互作用力指数判别重要廊道，对潜在生态廊道进行分级[6,19]。具体计算方式如下：

(2)

式中：Gab是斑块a与b的相互作用指数；φa是斑块a的权重值；Ra是斑块a的阻力值；Sa是斑块a的面积；Lab是斑块a与b之间廊道的累积阻力值；Lmax是所有廊道累积阻力最大值。

3 结果与分析

3.1 生态安全格局的建立

对各阻力因子图层进行加权叠加，在成本距离模块中生成城市用地和生态用地的最小累积阻力面，将生态用地最小累积阻力面与城市用地最小累积阻力面做差值处理，分级构建区域生态安全格局(图1)。采用自然断裂法将差值分为高生态安全水平(-116 562.34～-67 733.19)、较高生态安全水平(-67 733.19～-35 831.48)、较低生态安全水平(-35 831.48～0)、低生态安全水平(0～14 299.77)、建设用地扩展区(14 299.77～49 456.75)5种生态安全类型。结果表明各生态安全类型呈“圈层式”分布，建设用地扩展区主要分布在城市建成区及南部地带；低生态安全水平区作为建设用地的“储备区”，主要分布在建设用地扩展边缘地带；较高以上安全水平区主要分布在北部及东部边缘带；较低生态安全水平区则作为高水平与低水平生态安全水平区域的过渡带。

3.2 潜在生态廊道及生态战略点的识别

基于生态用地扩张阻力面，运用成本路径模块生成各生态源地间的阻力路径，对源汇相同的阻力路径两两比较得到最小阻力路径。为识别重要生态廊道，计算各生态源地间的相互作用指数(表2)，对潜在生态廊道进行分级。若相互作用力指数大于1，则为一级生态廊道；若该值小于1，则为二级生态廊道。最后运用网络分析模块建立拓扑关系提取最小阻力路径的交汇点，作为研究区的生态战略点，从而形成“生态战略点—生态廊道—生态斑块”的“点、线、面”生态网络架构。

图1 生态安全格局分级

表2 各生态源地间的相互作用力指数

注：1号生态“源”位于玉泉区南湖湿地公园；2号生态“源”为横跨榆林镇和黄合少镇的大面积草地斑块区；3号生态“源”为黄合少镇赛音不浪村、苏计村的大面积林地管护区；4号生态“源”位于乌素图国家森林公园；5号生态“源”为保合少镇与榆林镇交界的大面积林地斑块区；6号生态“源”位于大青山森林保护区。

根据统计分析，共提取15条潜在生态廊道，总长度为315.63 km，生态战略点共9个。相互作用力指数大于1的一级生态廊道有1→3，2→3，2→5，2→6，3→5，4→6，5→6共7条廊道，长度为129.31 km，占总长度的44.14%，大致分为北线、东线和南线(图2)。北线主要位于大青山山前地带，是城市生态环境的天然屏障。其中，5→6号廊道的相互作用力指数最高，为70.669，表明5号到6号生态源地的累积阻力值最小，生境适应性最强，对于物种的迁移、扩散至关重要。南线则以大黑河为媒介，将城乡、河流、林地斑块连为一体，对南线廊道的建设也是提高城市人居环境的重要工程。东线是连接南北廊道的关键纽带。其中，2→3，2→5，3→5号廊道共同构成“生态三角”，利于生态流的交换。此外，对于一级生态廊道沿线地带，应加强生态绿化建设，强化生态保育。对于二级生态廊道沿线地带，在生态规划中也应加强公共绿地与湿地建设，不断提高城市人居环境。

图2 潜在生态廊道及生态战略点分布

3.3 生态网络框架的建立

根据现有建设用地分布、生态安全格局及现阶段城市空间结构规划图，构建生态网络框架，即“一环、两带、两轴、多中心”的生态网络发展模式(图3)。“一环”即基于一级生态廊道，采用邻域分析工具构建宽度为2 km的生态缓冲区，此环通过连接两条生态涵养带构成环城“生态屏障”；“两带、两轴”与城市空间结构规划相一致，即“大青山生态涵养带、大黑河生态涵养带、城市中心十字发展轴”；“多中心”即生态源点及生态战略点共同组成城市发展的生态绿心。“一环、两带、两轴、多中心”的生态网络框架是在原有“两带、两轴、三心”的城市空间结构规划基础上，以生态源地为核心，多生态战略点为线索，一级生态廊道为骨架构建的城乡一体化绿色生态网络，对于提升区域生境质量、提高生态服务价值、保护生物多样性及优化城市人居环境大有裨益。

图3 生态网络框架的建立

3.4 城市区域空间管制分区

根据所建立的生态安全格局及生态网络框架(图1，2，3)，基于主体功能分区思想，以乡镇尺度为单位识别禁止开发、限制开发、重点开发以优化开发4类主体功能分区(图4)。

图4 城市区域管制分区

其中，禁止开发的乡镇为保合少镇，占研究区总面积的23.2%。禁止开发区作为城市重点生态功能区，原则上禁止任何形式的土地开发活动，需引导人口、经济活动向中心城市集聚，提升大青山生态源地的生境适宜性，加快推进大青山南坡生态修复工作，提升大青山生态功能区的抗风险能力；限制开发的乡镇有攸攸板镇、榆林镇及黄合少镇，占研究区总面积的36.4%。与中心城区相较，限制开发区人口、经济功能较弱，生态建设易被忽视。故在城市建设过程中要严格保护耕地，大力发展绿色生态农业，严控耕地面源污染。此外，榆林镇与黄合少镇作为连接南北生态廊道的关键纽带，起到提升生态源地连通性的重要功能，对稳定城市内部绿色低碳的功能不容忽视，需全力维护该区域内的生态战略点，强化生态绿心建设；重点开发的乡镇有赛罕区、金河镇、小黑河镇，占研究区总面积的28.8%。重点开发区是未来城市用地扩展的主要发展方向，要引导城市重心向南部、东南部迁移，发展城市新区，加快基础设施投资建设，引导人口、产业合理集聚，权衡城市用地扩张、耕地保护、生态保育的土地利用格局优化。同时，应加强对大黑河的污染防治，特别是途经中心城市的河段，提升水质并加强对水资源的循环利用至关重要；优化开发的乡镇有回民区、玉泉区及新城区，占研究区总面积的11.6%。作为开发密度较高、人口过度集聚的地区，要加快转变产业发展结构，并注重城市内部挖潜，推行“空间职能置换”战略。

4 结论与讨论

本文以呼和浩特市区为研究区，从研究区自然环境、资源条件、人类干扰等角度出发，构建区域用地扩展阻力评价指标体系，在此基础上通过MCR模型构建区域生态安全格局，识别潜在生态廊道与生态战略点，构建区域生态网络框架，对城市区域空间管制进行分区，得到的主要结论有：

(1) 对区域生态安全格局进行分级，各生态安全格局类型区呈“圈层式”分布。较高以上生态安全水平区合占研究区总面积的22.47%，生态安全水平总体偏低。

(2) 根据生态用地扩张阻力面共提取15条潜在生态廊道，总长度为315.63 km，生态战略点共9个。相互作用力指数大于1的一级生态廊道有1→3，2→3，2→5，2→6，3→5，4→6，5→6共7条廊道，长度为129.31 km，占总长度的44.14%。

(3) “一环、两带、两轴、多中心”的生态网络框架是在原有“两带、两轴、三心”的城市空间结构规划基础上，以生态源地为核心，多生态战略点为线索，一级生态廊道为骨架构建的城乡一体化绿色生态网络架构。对于提升区域生境质量、提高生态服务价值、保护生物多样性及优化城市人居环境具有重要意义。

(4) 重点开发的乡镇有赛罕区、金河镇、小黑河镇，占研究区总面积的28.8%。重点开发区是未来城市用地扩展的主要发展方向，要引导城市重心向南部、东南部迁移，发展城市新区，加快基础设施投资建设，引导人口、产业合理集聚。

构建科学合理的生态安全格局对于改善干旱地区生态环境脆弱，推动区域可持续发展具有重要的现实意义[12]。本文所建立的生态网络框架，是鉴于现行规划中生态战略点、生态廊道研究不足的实际。诚然，在构建区域生态安全格局的过程中，对于用地扩张阻力评价指标体系还未形成统一标准。鉴于研究区生态环境基底脆弱，因而在指标选取过程中更加侧重于对植被、水土等生态因子的考量。结果证明所识别的潜在生态廊道、生态战略点符合研究区实际，验证了对用地指标选取的合理性。此外，文本所优化的生态安全格局是基于研究区发展现状，若能与研究区近景年城市用地扩展预测模拟相结合，以乡镇为单位分区进行土地利用合理规划，会使研究结果更具参考价值，特别是对于控制城市蔓延并遏制生态破坏问题[20-21]。