随着社会经济发展和城镇化进程加快，人类活动对生态系统的质量和数量产生严重破坏，水土流失、物种减少、森林锐减等生态问题突出，国土空间生态安全遭到显著威胁^[1]。生态空间作为人类生产生活的基础场所，保障其安全成为缓解人类社会发展与生态环境保护之间矛盾的重要解决途径^[2]。构建生态安全格局作为城市生态安全保障的关键环节，同时也是生态系统管理的关键步骤。国家以城市化、农业发展和生态安全三大战略格局为目标的国土开发格局作为生态文明建设的首要任务，其重要性不言而喻^[3]。

生态安全格局源于景观安全格局，将景观空间格局与生态过程进行耦合，既保留了格局与过程相互作用原理，又强调了生态问题的作用机制^[4]。自20世纪末以来，生态安全格局这一概念逐渐得到国内外学者的重视，许多学者从理论方法、指标制定、模型构建等方面开展了深入研究，其成果在生物多样性保护^[5]、生态红线划定^[6]、矿区生境治理等^[7]方面得到广泛应用。目前，采用生态源地确定—构建阻力面—生态廊道识别模式对生态安全格局进行构建已成为该研究的基本范式^[8]。生态源地是承担重要生态系统服务关键生态用地^[9]。多数学者将大型水域和林地斑块^[2]、生态红线^[10]、风景区^[11]、重要生态区^[12]作为生态源地。生态阻力面构建是构建生态安全格局的重要一步，反映物种迁移过程及信息能量流动等生态功能受到土地利用方式和人为干扰的程度，已有研究多从土地利用类型以及坡度、高程等角度进行表征，缺少反映人类活动的指标因子。生态廊道是具有保持生态系统稳定的通道，它涉及生态、社会、文化多功能并连接陆域和水域生态源地^[13]。最小累计阻力模型(minimum cumulative resistance，MCR)是定量识别生态廊道最有效的方法，相较其他模型能很好地模拟景观对空间运动过程的阻碍作用，从而更好的表达景观格局和生态过程之间的相互作用关系^[8]，因而得到广泛应用。近年来已有学者对干旱区^[14]、矿区^[7]、火山熔岩湿地^[15]等生态敏感区以及北京市^[16]、雄安新区^[17]、黔中城市群等^[18]重点经济发展城市或城市群进行生态安全格局的构建。例如，冯琰玮等^[19]对内蒙古生态承载力进行评价，将识别的生态优势斑块和自然保护区作为源地，通过MCR模型和“三生”用地类型构建、优化生态安全格局；李涛等^[20]以湖南省衡阳市为研究区，基于InVEST模型HQ模块和电路理论构建了城市生态安全格局；王浩等^[21]从生态系统服务重要性和敏感性两个角度识别广东省的生态源地，构建广东省生态安全格局。

济南市由于其独特地质地貌条件形成“清泉石上流”的绮丽风光，成为中国岩溶泉的典型代表^[22]。近年来，济南市在快速推进城市化过程中对生态环境产生一定破坏，环境污染加剧，泉水停涌的危机日益显现，“家家泉水，户户垂杨”的泉城景色已不复存在^[23-24]。资料显示，由于过量开采地下水、地面硬化面积增加、砍伐树木等原因，近50 a来济南的泉水补给量已减少约30%，面临枯水年泉水断流或停喷的威胁^[25]。当前，国内外专家多从泉水来源^[26]、泉流量影响因素^[27]、保泉措施等^[28]方面对济南泉域进行研究，尚未将泉水生态安全与济南生态安全格局构建进行联系。本研究以济南市为研究区域，以“山水林田湖草泉”系统为研究对象，识别生态源地并通过MCR模型构建生态廊道，从而构建具有济南市城市特色的生态安全格局，以期为济南市实现生态环境的保护与修复、区域协调发展提供思路。

1 研究区概况

济南市位于山东省的中西部、黄河下游、华北平原东南部，地理位置为北纬34°46′21″—37°32′13″, 东经116°13′45″—117°58′10″。东与淄博市接壤，西与聊城市、德州市相连，南与泰安市毗邻，北与滨州市相接，是山东省省会、济南都市圈核心城市、中国的历史文化名城，也是重要的交通枢纽。全市下辖12个行政区，总面积10 244 km²，2019年全市常住人口8.91×10⁶人。济南南依泰山北跨黄河，地势南高北低，最高海拔1 108.4 m，最低海拔5 m。济南年平均气温13.6 ℃，年平均降水量614.0 mm，降水集中于夏季，属温带季风气候。境内河流较多，有黄河、小清河、徒骇河等，湖泊有大明湖、白云湖、雪野湖等。区域内泉水及泉群众多，以清冽甘美闻名于世，主要风景名胜有趵突泉、五龙潭等。

2 数据及研究方法 2.1 数据来源

本研究涉及的数据主要包括土地利用数据、DEM数据、人口空间分布公里网格数据、经济社会数据以及基础地理数据。具体数据来源及使用见表 1。

2.2 研究方法 2.2.1 生态源地综合识别

参考国际城市与规划学会^[29]在南京市江北新区城市规划设计中提出的“蓝绿”系统的概念，结合济南市“山水林田湖草泉”系统的特点，通过模型法和数据提取等方法，从“蓝、绿、泉”3个角度选择生态源地。其中，“绿”指以保护绿色生态资源为目标，选取能够对维持陆地生态系统稳定、促进区域可持续发展具有重要作用的生态斑块，通过InVEST模型计算得出；“蓝”指河流、湖泊、水库等对城市水资源生态保护起到关键作用的生态要素，“泉”是对维持济南泉水持续喷涌起到重要作用的生态斑块，“蓝”和“泉”的生态源地通过基础地理数据提取得出。

(1) “绿源”识别。近年来济南市快速城镇化过程中将大量的生态用地用于城市建设，特别是林地和耕地被硬化成不透水面，导致生态环境恶化，生物多样性和物质、能量流动受到干扰，生态系统为城市提供生态资源的能力降低。生境质量是衡量区域能够提供生物生存繁衍条件能力高低的重要指标。生态源地是生态系统中具有重要作用的斑块，其本身应具有较高的生境质量，能够为物种生长、发育繁衍提供良好的条件，在协调区域生态安全与可持续发展过程中具有重要作用^[30]。本研究选择生境质量指标对济南市陆地生态环境进行评价进而识别生态源地。InVEST模型是由美国斯坦福大学开发的、能够权衡区域生态系统服务的开源代码软件，其中的生境质量模块在假设生境质量是连续的条件下，以生境质量的高低代表空间上生物多样性的分布情况^[31-32]。该模型基于土地利用数据、生境质量威胁因子的影响距离、不同生境对威胁因子的敏感性4个变量信息生成区域的生境质量地图，并且考虑了人为活动对生态系统的干扰。计算方法如下：

$ D_{i j}=\sum\limits_{r=1}^{R} \sum\limits_{y=1}^{Y_{r}}\left(w_{i} / \sum\limits_{r=1}^{R} w_{r}\right) r_{y} i_{r i y} \beta_{x} S_{j r} $

(1)

$ Q_{x j}=H_{j}\left(1-\frac{D_{x j}^{2}}{D_{x j}^{2}+k^{2}}\right) $

(2)

式中：D_ij为LUCC中j类型土地中栅格i的生境退化度；R为选取的威胁因子个数；Y_r为r威胁因子栅格图上的栅格的数量；r_y为栅格图层y上的威胁因子值；w_r为威胁因子的权重，通常对生态系统影响越高的威胁因子的权重越接近于1；i_riy表示生态系统生境栅格i与威胁因子r栅格的距离；β_x为威胁因子的可达性；S_jr为生境类型j对威胁因子r的敏感度；Q_xj为生境类型为j的栅格x的生境质量指数；H_j为生境类型j的生境适宜性；D_xj为生境类型j的栅格x的生境威迫水平；k为半饱和常数。

如表 2—3所示，研究选取建设用地、裸地、其他用地作为生态系统的威胁因子，并且考虑在城市建设过程中交通用地会对生境质量产生负面影响，将高速道路、国道、省道共同作为威胁因子。威胁因子的最大影响距离及权重、济南市各生境类型对威胁因子的敏感性等模型参数的确定参考InVEST模型用户指导手册以及多位专家学者^[32]的研究。生态源地的破碎化程度对生态系统的物质和能量传递会产生一定影响^[33]。为降低济南市生态源地破碎化程度，将图斑进行聚合，参考康洁铭等^[34]的研究，将高等级生境质量累计面积占总面积的80%的区域作为“绿源”。

(2) “蓝源”识别。水生态系统是城市生态安全保护的一个重要方面，济南市水生态系统可分为河流、湖泊、水库、泉4种类型^[35]。其中，黄河、小清河两大流域共包含120条大、中、小型河流。济南市水生态文明建设提出“河湖联通惠民生、五水统筹润泉城”的发展格局，将河流、林地连接构成生态网络，维持城市水生态廊道畅通。参考济南市水生态文明建设成果以及《济南市城市发展战略规划(2018—2050)》^[36]，将市域范围内玉符河、大辛河在内的六条河流廊道以及白云湖、玫瑰湖在内的5个湖泊作为“蓝源”，同时将水库累计面积占总面积75%的卧龙山水库、雪野水库、鹊山水库等纳入蓝色生态空间，构建水源生态网络。泉水出露也会在周围形成小型湿地，将市域范围内形成公园绿地的趵突泉泉群、五龙潭泉群、珍珠泉泉群、百脉泉泉群等四大泉群范围内的公园绿地作为湿地系统，形成泉水蓝色源地。

(3) “泉源”识别。济南市天然泉眼共计800余处，泉水保护是彰显城市特色、建设城市名片、弘扬泉水文化的重要举措。“保水源、保路径”成为维持泉水生态循环过程的关键方向。重点渗漏带的土地利用、生境质量等对维持泉域生态系统平衡、济南特色资源保护以及防洪排涝等具有重要作用，是“保水源”的核心区域。重点渗漏带是指由于地形地质条件、含水层特性等因素形成的地表水向地下水补给能力强、速度快的区域^[24]。根据管清花等^[37]的研究，重点渗漏带可根据受到人为干扰的程度划分为12处近自然型重点渗漏带和12处人为干扰型重点渗漏带。近自然型渗漏带人为干扰较小，基本保持了较好的植被覆盖度、发育完善的构造和裂隙，渗漏性强，地下水补给条件较好。人工型重点渗漏带中7处处于正在开发建设过程，地面硬化、植被破坏，部分丧失地下水补给能力，5处已经被建设用地覆盖，地下水补给功能已经完全丧失。生态源地的选择考虑具有良好地下水补给条件的12处重点渗漏带，包括玉符河渗漏带、店子—二仙渗漏带、龙洞渗漏带、涝坡—郑家窝坡渗漏带、大涧沟渗漏带、车脚山—太平庄渗漏带、河东—河西渗漏带、泉泸—钱家庄渗漏带等，不再考虑城市建筑部分或完全覆盖并且丧失渗漏补给地下水功能的羊头峪、双龙等12处重点渗漏带。

2.2.2 构建阻力面

区域内物质和能量的流动和相互联系需要克服景观的阻碍作用，土地覆被类型、地势地貌和人为活动都会对源地扩张产生负面影响^[8]。参考已有的研究^{[12, 18]}并结合济南的自然和社会经济条件，从自然地理、人为干扰两个方面选取指标构建源地扩张的阻力面(表 4)。各指标因子的权重通过AHP法求得，采用加权求和方法构建阻力面并按照从低到高值划分为6个等级。

2.2.3 生态廊道与生态安全缓冲区识别

1992年Knaapen提出最小累积阻力模型，可以计算从源到目的地克服阻力面障碍所做的功^[38]。本研究采用MCR模型定量识别济南市的生态廊道、构建生态安全格局，基本公式(3)如下：

$ \mathrm{MCR}=f_{\min } \sum\limits_{j=n}^{i=m} D_{i j} \times R_{i} $

(3)

式中：MCR为最小累积阻力值; D_ij为物种从源地j到目的地i的空间距离; R_i为目的地i的阻力系数; f为最小累积阻力与生态过程的正相关关系。

生态安全缓冲区作为生态源地与城镇建设开发区的过渡带，对保护生态源地起着至关重要的作用。根据相关研究^{[18, 32]}以及济南市实际情况，将济南市最小累积阻力面按照自然断点法划分为其他区域、高水平缓冲区、中水平缓冲区、低水平缓冲区。

3 结果与分析 3.1 生态源地提取 3.1.1 “绿源”识别结果与分析

基于InVEST模型得到济南市生境质量评价结果(图 1a)，栅格值范围在0~1之间，值越高代表生境质量越好，维持生物多样性的功能越强。参考郝月等^[32]的分类标准，将生境质量等间距分为低、较低、中等、较高和高5个等级，0~2为低等级，该范围的区域生态系统功能遭到严重破坏，受到人为干扰较大，难以满足生物多样性和生态系统可持续发展的要求，0.8~1范围内生态系统能够为区域提供良好的生态资源，生态功能稳定。

通过ArcGIS 10.5的面积制表功能，得到济南市不同等级生境质量的面积和比例，结果如表 4所示。济南市平均生境质量为0.433，总体生境质量在中等级以下，面积为8 401.58 km²，全市较低等级生境质量面积最大，占总面积的62.79%，较高和高等级分别占5.55%，12.44%。分别统计各区各乡镇的生境质量，生境质量最高的区县为历城区，平均值达到0.53，生态系统能够满足一般的城市服务需求；生境质量最低的是槐荫区，处于较低等级，由于该区域内建设用地面积较大，人类活动对生境质量产生的负面影响较高。从各区县不同等级生境所占比例来看，历城区的高等级生境质量最高，占区域的29.04%，其次为历下区，占区域面积的20.12%。从空间上看(图 1b)，济南市生境质量呈“盆地型”分布，南部各区县生境质量较高，其中西营街道平均生境质量为0.844，最低的是工人新村南村街道，历下区、天桥区、市中区、槐荫区和历城区五区交界处各乡镇、街道的生境质量均小于0.3，生态系统服务能力较低。济南市生境质量南高中低的原因主要是南部为泰山余脉，森林覆盖广泛，这些区域的生态用地受人为活动威胁较小，生态系统服务功能稳定；中间为济南市主城区，大部分为不透水面，绿化面积较小且呈零星分布，城市道路密度较大，成为生态系统功能的重要威胁因子。

“绿源”识别结果如图 2a所示。“绿源”总面积为977.81 km²，各地类面积及比例如表 5所示。其中90.07%的土地利用类型为林地，多数分布在南部山区，呈现出连片分布的特点，并且生态源地之间的连通性较好，部分分布在平阴县南部、莱芜区北部和钢城区北部。山前平原区以耕地为主，人为活动对生境质量产生较大影响，未能识别出生态源地。草地、灌木林地和耕地分别占总面积的4.20%，2.92%和2.61%。“绿源”包括了济南市柳埠自然保护区、药乡国家森林公园、子房洞景区等一级生态控制线内的区域，表明识别结果较为合理。

3.1.2 “蓝源”及“泉源”识别结果与分析

本研究将黄河、小清河、嬴汶河、大汶河、徒骇河、秀江河、西八漏河等主要河流作为河流生态源地，主要分布在济南市中部、东北部；识别出面积累计占总面积80%以上的30个大中型水库，水库源地以南部山区、莱芜区、钢城区为主要分布地；玉清湖、大明湖、鹊山龙湖、华山湖、白云湖五大湖泊为湖泊源地。将趵突泉泉群、五龙潭泉群、珍珠泉泉群、百脉泉泉群等四大泉群范围内的公园绿地作为泉水蓝色源地(图 2b)。重点渗漏带作为保泉的核心区域对保障泉水持续喷涌具有重要作用，12个“泉源”主要分布在南部山区北侧(图 3)。

3.1.3 生态源地综合识别结果与分析

将“蓝、绿、泉”源地叠加，获得济南市总生态源地(图 4)，总面积为1 249.32 km²，占济南市总面积的12.20%。其中，历城区分布面积为397.50 km²，占总面积31.82%，其次为莱芜区，占总面积的231.09%。商河县仅有7.54 km²的生态源地，占0.60%。

3.2 阻力面构建结果与分析

从自然地理、人为干扰两个方面构建的源地扩张的阻力面阻力值位于10~100之间，按照自然断点法将数据分为低、较低、一般、较高、高5个等级，结果如图 5所示。济南市阻力面呈现出“中间高，南北低，多团块”的特点，即从中心城区向南北降低的趋势。济南市主城区以及莱芜区、章丘区、钢城区等区县中心城区建设用地较为集中，并且人口密集，物质流和能量流受到的阻力较高。低阻力区域集中分布在南部山区，以林地、草地、耕地等为主，人类活动频率较低，虽然坡度较陡，但综合阻力小，对周边的物种扩散较为有利。

3.3 生态廊道识别结果与分析

研究共识别出生态源点111个。通过Cost-Distance工具提取的生态廊道共有119条，全长1 395.37 km，最长生态廊道187.92 km，最短为0.97 km。空间上西部地区生态廊道呈西南向东北走向分布，东部为东南—西北走向，南部呈东西分布，连通了东南西北四方的生态源地。分布广泛的河流构成天然廊道，其中黄河面积宽分布长，在研究区起到了重要、稳定的物质能量信息流动带的作用，小清河沟通了南北走向的韩仓河、巨野河、大辛河等，形成环状的廊道网络；南部山区由于生态源地广布成为“绿源”、“泉源”廊道的聚集区，各生态源地之间的廊道通过水库、河流连通，且由于源地分布面积较广，各生态源地连通性好，廊道较为稀疏；章丘区生态源地较为分散，廊道以陆地廊道为主；长清区、平阴县、莱芜区和钢城区廊道由河流与陆地共同构成，河流协调了物质能量流动功能，为生物迁移起到较好的桥梁作用。

整体上济南市形成了联通南泰山北黄河的“山河通廊”，稠密的河网有效提高各生态用地的连通性。生态廊道的宽度受到物种种类、生态环境等有关因素的影响，济南市适合物种迁移的陆地生态廊道宽度设置为12~30 m^[39]。

3.4 生态安全缓冲区识别结果与分析

济南市生态安全缓冲区结果如图 6所示。高、中、低水平安全缓冲区面积分别为1 495.5，1 382.19，1 262.87 km²，占济南市总面积的40.42%。其中高水平安全缓冲区为生态源地重要保护区，具有重要的生态保护价值，应禁止进行任何开发建设活动，保障区域内泉水渗漏、生物多样性保护等生态系统服务功能不受人类活动的影响。

城市发展建设、人们生活与济南的泉水特色息息相关，生态泉城建设应重点关注泉水生态环境，为避免城市建设蚕食重点渗漏带，政府应评估该区域的建设活动，保障泉水与城市可持续发展。中水平安全缓冲区为一般保护区，限制该区域城镇开发活动，尽量避免天然岩溶地质环境破坏，影响大气降水渗入，同时应注意一般保护区内生境质量较高的斑块，鼓励进行植被修复、山体绿化等生态保护活动。低水平安全缓冲区为建设缓冲区，在该区域可以进行建设活动，但需要注意建设的规模和强度，科学合理进行规划，保证生态系统功能不受影响。

3.5 济南市生态安全格局

通过InVEST模型、基础地理数据提取等方法，提取济南市生态源地、生态廊道的空间分布形态，根据《济南市城市发展战略规划(2018—2050)》、水生态文明建设成果，将研究区湿地生态区、泉水景观区和生态保护区作为重要分区，构建济南市生态安全格局。济南市整体上呈现“一带、三区、两翼、九点、多廊”的格局(图 7)。“一带”是指黄河生态保护带，贯穿济南市，保障区域的水资源环境，与沿岸湿地共同构成生态景观带，成为济南的“生态脊”。“三区”指济南南部生态保护区、城市泉水景观区、沿河湿地生态区。其中，南部生态保护区包含城市大部分绿源和泉源，保护城市生物多样性和泉水生态安全；泉水景观区以天然河流廊道为主，通过水库湖泊连接小清河、黄河与南部“绿源”，提供泉水观景功能；沿河湿地生态区主要以黄河北岸形成的湿地为主，提供滞蓄洪水、景观游览功能。“两翼”指长清区、平阴县和莱芜区、钢城区，生态廊道连接“绿源”和“蓝源”。“九点”指环绕中心城区的卧虎山、锦绣川、狼猫山、杜张4座水库和玉清湖、大明湖、鹊山龙湖、华山湖、白云湖5个湖泊及毗邻的相关湿地等。“多廊”指城市的天然河流廊道和潜在生态廊道，主要河流廊道包括德惠新河、徒骇河、小清河、玉符河、北大沙河、大辛河、韩仓河、巨野河、绣川河、大汶河，潜在生态廊道以识别出南部山区以及各区县的“绿源”、“蓝源”连接廊道为主。

4 讨论与结论 4.1 结论

本文选取济南市作为研究对象，基于生态学理论以及泉水相关知识，从“山水林田湖草泉”系统出发选取生态源地，基于MCR模型识别研究区生态廊道，构建济南市“一带、三区、两翼、九点、多廊”的生态安全格局。

(1) 通过生境质量评价以及生态保护需求，识别出济南生态源地共1 249.32 km²，“绿源”主要分布在南部泰山余脉，“蓝源”以区域内主要河流、湖泊、大中型水库和四大泉群的绿地为主，“泉源”分布在南部山区北部。

(2) 研究区阻力面以主城区为中心向四面降低，各区县城区呈现小范围的阻力高值区，研究表明人类活动对物种迁移影响最大。

(3) 济南市构建了“一带、三区、两翼、九点、多廊”的点—线—面生态安全格局。生态源点共111个，为生态源地的斑块几何中心，生态廊道共有119条，全长1 395.37 km，黄河天然廊道从研究区西南部延伸到东北部，距离较长，联通了济南市北部与南部生态斑块，成为区域重要廊道。

4.2 讨论

(1) 本研究基于“蓝—绿—泉”系统识别济南市生态源地，避免了主观干扰，并且考虑了济南泉城特色，有利于缓解城市发展与生态安全、泉水保护之间的矛盾。本文采用了InVEST模型提取生态源地，受到模型和数据精度的影响，最终结果可能会存在一定的局限性；其次，生态系统不完全行政边界的制约，周边区域的物质、能量流动也会对生态安全产生影响，因此在今后研究中应重点关注行政单元周边各要素的相互作用。

(2) 根据《济南市发展战略规划(2018—2050)》内容，未来将形成轨道支撑的多中心组团式结构，老城、长清以及雪野等发展片区按照规划内容将会与生态源地部分重叠。其中老城发展片区主要影响保泉核心区域—重点渗漏带，包括人工型渗漏带和近自然型渗漏带，泉水的喷涌受到一定威胁；雪野和长清等发展片区对部分林地产生一定影响，生态环境可能会遭到一定破坏；城市轨道交通的发展也会占用部分生态用地，影响正常的生态系统服务功能。因此，在城市发展过程中要注意对生态源地的保护，特别是重点渗漏带区域，建设用地尽量避开生态源地以及生态廊道，为物种生存及迁徙提供良好的条件。以城市生态生态系统安全为底线，基于生态安全格局进行城市规划，构建城市与自然和谐发展格局，以实现“蓝绿交织，清新明亮，泉城共融”的济南城市布局。

(3) 生态安全格局构建不仅为济南2019年行政区划调整后生态规划工作提供参考借鉴，也可为政策实施等工作提供指导建议。根据本文研究结论，南部山区作为城市的生态源地集聚区对城市安全至关重要，因此济南市政府相关部门在进行生态环境保护时，要结合城市特色，落实《济南南部山区保护与发展规划》，贯彻“南控”方针，以“蓝天工程”、大环境绿化工程为重点提高森林覆盖率、涵养城市水源，构建“山水林田湖草泉”一体的现代化生态济南。