张 林，田 波，周云轩，朱春娇

（华东师范大学 河口海岸学国家重点实验室，上海 200062）

0 引 言

随着长三角区域城市化进程的不断推进，城市人口不断增加，建成区快速扩张．以浦东新区为代表的长三角城市在促进区域经济和社会发展的同时，也面临着一系列城市生态问题．表现在不断变化的土地利用方式影响着城市生态系统的机能、自净能力和景观承载力；城市自然生态景观要素数量不断减少、类型趋于单一；区域内自然生态单元要素特别是湿地、绿地、林地不断被侵占，表现破碎化、岛屿化；城市生态空间连通性降低［1，2］，生态空间网络整体性不断破坏，生态服务功能不断下降［3，4］．

为改善城市生态现状，国内外一些学者提出了构建城市生态网络的思路［5，6］，用于加强自然环境的保护与恢复，减少城市生态用地的破碎化，提升城市生态系统健康水平和生态服务功能［7］．2011年，上海市首次提出并发布建设多层次、成网络、功能复合的城市生态网络规划，将城市生态网络建设纳入生态文明建设的重要组成部分．

生态网络是指自然的或植被稳定的以及依照自然规律而相连接的空间，主要以湿地、绿地、耕地和林园地为主，强调自然的过程和特点［8］，是由生态景观及其之间的连线所组成的结构完整、功能完善的生态系统［3］．城市生态网络建设的实质是利用植被带、河流和农田等生态单元，通过线性廊道将城市中分散孤立的各种类型生态斑块联系起来，形成一个结构成熟、功能完善的生态系统［9］．国外对于构建生态网络的研究开展较早，从20世纪80年代开始，欧美就出现了一些生态网络规划的理论和实践，如欧洲的生态基础设施（Ecological Infrastructure）、生态架构（Ecological Framework）及美国的绿道体系（Greenway System）等［10］．国内有关生态网络的研究起步较晚，近年来取得了较大发展，如通过景观指数［6－12］、网络分析［13］等方面对于生态网络结构进行评价分析，进而对生态网络优化提供参考．但国内大多研究关注单一生态要素，如绿地［8－14］或湿地［15］等，或者是仅以某一角度如景观角度对生态本底数据的自然、社会属性（如面积、保护破坏程度、植被覆盖度等）进行深层次的分析比较．本文主要是以遥感和地面调查技术为基础，在辨识浦东新区林绿地、湿地生态空间分布现状基础上，根据生态网络单元基础数据自然和社会属性，进行分等定级，构建多层次的权重指标体系，利用最小费用法构建浦东新区潜在生态廊道，采用重力模型提取重要廊道，通过缓冲区分析对生态廊道和网络结构进行分析，为浦东新区城市生态网络构建和优化提供科学依据．

1 研究区概况

上海市浦东新区位于上海市东部，长江三角洲东缘，东临东海之滨，属长江口冲击平原，为平原河网感潮区，其地理范围在（121°26′E～121°59′E，30°50′N～31°23′N）之间．它集滨江临海之利，有着较丰富的生态资源，是上海市经济发展的引擎，国家重大发展战略区之一，也是我国快速城市化的典型代表区域．随着浦东新区20世纪90年代开发开放，特别是近10年以来，浦东新区城市化进程不断加快，快速城市化给浦东新区生态环境带来了巨大影响［16］．城市建成区迅速扩张，土地利用方式发生巨大变化，许多生态用地转为建设用地［17］，生态用地连通性下降，斑块破碎化严重，城市生态网络结构和生态功能面临较大威胁．本文选取浦东新区为研究区，以现浦东新区行政区划为界（见图1），包括浦东新区所辖陆域和九段沙湿地国家级自然保护区在内的湿地水域，区域总面积27．5万hm2，其中林地总面积约1．9万hm2，绿、草地面积约1．38万hm2，湿地总面积14．71万hm2，湿地植被面积为0．68万hm2．

图1 研究区位置（底图为2011年福卫2影像）Fig．1 Study area with satellite image of FW2 acquired in 2011

2 方 法

通过遥感和地面调查相结合的方法获取浦东新区绿地、林地、湿地、道路交通和城市建设用地等地类，按照林绿地与湿地生境类型、面积、保护程度等属性信息，对生境斑块进行多层次的分等定级，选择重要的生态斑块作为核心斑块，建立多层次的生态斑块适应性评价体系，构建景观阻力图（即消费面）．按照最小费用法，空间计算各核心板块间潜在生态廊道．使用重力模型分析核心斑块间相互作用力来选择重要生态廊道．对重要生态廊道进行缓冲分析，分析生态廊道缓冲区域土地利用结构，得到研究区域生态网络布局和结构现状．采取技术路线详见图2．

图2 技术路线图Fig．2 Methodology flow chart

2．1 生态网络单元空间分布现状识别与提取

研究区域内生态网络单元空间分布的识别与提取主要采用遥感、GIS技术和地面调查相结合的计算机提取方法，采用的数据主要有2011年10月福卫2号卫星（FW2）遥感数据（空间分辨率2m）、1∶10 000地形图、上海市2011年海图等．对FW2号遥感数据进行辐射校正、几何校正和图像增强等预处理后，基于决策树的eCognition自动分类与目视解译相结合的方法进行生态网络要素的识别和土地利用解译．通过湿地数据进行百分之百野外调查验证，获得准确的湿地空间信息和属性信息．结合浦东新区林绿地详查数据，以FW2号遥感影像数据为基础，得到研究区内林绿地分布数据．参照交通图和FW2数据，利用GIS数字矢量化提取区域内道路交通数据．依照GIS空间数据库构建标准，整合生态要素数据，得到浦东新区林地、绿地、湿地等生态网络单元分布空间数据库和土地利用现状（见图3）．

图3 2011年浦东新区土地利用现状图Fig．3 The spatial distribution of land use classification in 2011 in Pudong New Area

2．2 生态网络核心斑块选择与确定

生态网络核心斑块选取主要考虑生态斑块的景观类型、斑块面积、植被覆盖度、受保护等级、生态功能和空间分布格局等多方面因素［14］．本研究将面积大于5 hm2的林、绿地、面积大于50 hm2的湿地作为备选核心斑块，筛选其中生态区位重要、受保护状况良好的斑块．最终本文选取了18个生态单元地块作为核心斑块（见图4）．核心斑块主要包括自然保护区（16、17、18号斑块）、自然保护小区（12、15号斑块等）、森林公园（1、14号斑块等）、湿地公园（3、6号斑块等）、绿地公园（2、4、5、9号斑块等）及大面积苗圃（11、12号斑块）等，总面积为25 920．42 hm2，占区域内总生态用地面积的10．71%，其中湿地面积20 113．02 hm2，林地面积2 851．61 hm2，绿地面积1 300．41 hm2，农田面积1 089．81 hm2．

图4 浦东新区核心斑块分布图Fig．4 The core ecological patches of Pudong New Area

2．3 生境斑块景观阻力消费面构建

潜在生态网络是由核心斑块间不同土地利用类型的斑块适应性和景观阻力所决定［18，19］．生境斑块适应性是指某一生境斑块对物种生存、繁衍、迁徙等活动的适应性程度．景观阻力是物种在不同景观单元之间进行迁移的难易程度，与生境适宜性成反比关系［20］．所有斑块景观阻力值的空间分布则构成景观阻力消费面．生境景观阻力值主要由景观类型决定，也受自然和社会属性影响．景观阻力值计算需要考虑景观类型、斑块面积、保护利用状况、植被覆盖度、健康度等综合因素．对此，在计算生态斑块景观阻力值时，采用计算公式如下：

其中，Ci为斑块i的景观阻力值，Mi为斑块i对应的景观类型的基准阻力值，Wj为斑块i的第j个属性分级后的景观阻力系数．在综合国内已有研究［8，20］的基础上，确定浦东新区不同景观类型基准阻力值（见表1）．对斑块不同属性及属性值对生境景观适应性影响程度进行分等定级，确定浦东新区不同景观及属性的景观阻力系数（见表2和表3）．湿地核心斑块主要选择湿地型、湿地面积、植被覆盖率和保护程度四个指标；林地绿地核心斑块，主要选择林种、面积、保护分级、郁闭度和生态功能分级等指标．根据公式（1），计算得到浦东新区区域内所有生境斑块的景观阻力消费基准面（见图5）．

图5 浦东新区景观阻力消费基准面Fig．5 The cost surface map of Pudong New Area

表1 不同景观类型基准阻力值Tab．1 Base impedance values of different land use type

表2 湿地景观阻力分级表Tab．2 Resistance coefficient of wetland

表3 林绿地景观阻力分级表Tab．3 Resistance coefficient of forest and grass land

2．4 潜在廊道分析与重要廊道提取

最小累计阻力模型是指从源到目的地运动过程中所需消耗最小代价的模型［21］，最初源于对物种扩散过程的研究，近些年已经用于适宜性评价［22］、选址［23］、生态格局规划与构建［24，25］等多方面．最小累计阻力模型基本公式如下：

其中，MCR为最小累积阻力值，Dij表示从源j到景观单元i的空间距离，Ri表示景观单元i的阻力值．以ArcGIS为空间分析软件平台，采用空间路径分析法得到18个核心斑块间最小阻力路径，即潜在生态廊道（见图6（a））．

图6 潜在生态廊道分布图（a）与利用重力模型提取的浦东新区重要廊道（b）Fig．6 The potential ecological network（a）and the important ecological network based on gravity method（b）

潜在生态廊道的有效性和连接斑块的重要性主要依据核心斑块与目标之间的相互作用强度［26］．基于重力模型［14］，本文构建了18个核心斑块间相互作用矩阵表（见表4），定量评定生态斑块间相互作用强度，以此来判断生态廊道的重要性．根据相互作用矩阵结果，提取相互作用力大于200的廊道，作为生态网络骨干廊道，得到浦东新区重要生态网络图（见图6（b））．重力模型计算公式如下：

其中，Gab是斑块a和b之间相互作用力，Na和Nb分别是两斑块的权重值，Dab是a和b两斑块间潜在廊道阻力的标准化值，pa为斑块a的阻力值，Sa是斑块a的面积，Lab是斑块a到b之间廊道的累积阻力值，Lmax是研究区内所有廊道阻力的最大值．

2．5 重要廊道缓冲区土地利用结构分析

生态廊道构建应该要有一定宽度，才能更好的发挥其生态服务价值．为研究浦东新区不同廊道宽度土地利用结构，以重要廊道为基础，通过缓冲区空间分析，分别得到10 m、50 m、100 m和500 m的缓冲区内林地、绿地、湿地、道路交通、建设用地等土地利用构成（见图7）．

表4 浦东新区生态网络核心斑块相互作用矩阵Tab．4 Node interaction in ecological network of Pudong New Area

图7 不同缓冲距离重要廊道土地利用构成Fig．7 Land use structure of different buffer distance

3 结果与讨论

（1）景观阻力分析

由浦东新区景观阻力消费面（见图5）可以看出，浦东新区景观阻力有明显的“北高南低”的特点．浦东新区北部区域因经济发展迅速，城市化扩张明显，大量建筑用地和道路使区域景观阻力偏高；在外环沿线、滨海带的防护绿地以及建成区内附属草地、林地等区域，景观阻力值较小，有明显的带状分布．

浦东新区南部区域因有大面积农田、湿地存在，整体景观阻力值较小，特别是南汇东滩禁猎区区域；景观阻力值较大区域位于各乡镇建成区；值得注意的是南汇东滩禁猎区北部的老港垃圾填埋场和临港新城区域景观阻力值较大，这将会给南汇东滩禁猎区带来一定影响．

（2）重要核心斑块分析

结合浦东新区核心斑块分布图（见图4）和景观阻力分布图（见图5）可以看出，浦东新区适宜生物生存栖息的大型生态单元斑块明显不足，且有南多北少的特点．其与区域城市化水平相关，由于经济的发展，浦东新区出现大量的非生态用地，将区域内大生境斑块割裂成小斑块所致．

通过表4可以看到，九段沙区域（斑块16、17、18）与其他斑块的作用强度较大，主要是由于九段沙湿地国家级自然保护区面积较大，并且景观阻力值小，可以反映出九段沙核心斑块在整个浦东新区生态格局中有着重要地位．外环楔形绿地（核心斑块7）在联系城市建成区内斑块（核心斑块2～5）与郊区大型斑块起到重要节点作用．滨江森林公园（核心斑块1）及世纪公园（核心斑块4）是在建成区内重要公园绿地，但由于建成区内斑块破碎化程度高、适宜生境少，与其他核心斑块联系较弱．

（3）生态廊道土地结构和利用分析

根据生态廊道10 m缓冲分析（见图7），区域内生态网络主要构成要素为林地，占51．39%，其大部分为道路防护林，占总廊道面积的39．25%．湿地在联系重要生态斑块中也起到了很大作用，占25．08%，其中占主导地位的为内陆河道．由于生态廊道大多为人类活动影响较弱的地区，廊道内建筑用地、道路等所占比例很小，为9．73%．

通过对生态廊道不同缓冲距离的分析（见图7）可以看出，在较窄廊道宽度（如10 m）下，生态要素占据着主要部分，所占比例为90%以上，特别是重要的湿地、林绿地资源，所占比例达76．47%，在发挥生态服务功能方面起到关键作用．但随着缓冲距离的增长，缓冲区内建筑用地由5．27%上升到23．40%，所占比例上升最快；而湿地与林绿地的总和则从76．47%降到39．34%，反映出浦东新区高质量的生态廊道宽度较窄，易受到外界影响而被阻断或者造成生态功能的削弱．

（4）浦东新区生态网络结构框架

生态斑块和廊道的有机整合对于生态网络服务功能有着重要作用．根据本文得到的浦东新区重要生态廊道，结合浦东新区生态要素的空间布局，规划提出“两核、一轴、两带、多线”的生态网络框架（见图8）．

“两核”：指九段沙与南汇东滩禁猎区．这两个区域有着大片的生态用地，而且生境自然度较好，受人为影响较小，是重要的源地．同时，这两个区域位于整个浦东新区的最前沿，在维护浦东乃至整个上海的生态安全有着至关重要作用．应该得到重点保护和合理管理．

“一轴”是指浦东沿海沿江带，包括浦东新区沿海滩涂带、沿海生态防护林带，是重要的生态栖息地，同时将浦东沿海各重要生态空间在纵向上联系起来，充分发挥生态服务功能．

“两带”是指外环绿地以及沿大治河和绕城高速形成的生态走廊，是浦东新区重要的生态骨架．将城区内重要生态节点联系成一体．

“多线”是指在横向上延伸的由湿地、林绿地组成的廊道，如联系城区绿核和九段沙核心区的张家浜生态廊道，以及联通沿大治河绿带和南汇东滩核心区的南部生态廊道等，都可作为浦东新区生态网络结构的补充和加强．

图8 浦东新区生态网络框架Fig．8 The ecological network structure of Pudong New Area

4 结 语

本项研究主要侧重于生态网络各要素空间分布现状及结构分析．通过针对不同景观属性，在考虑其各级属性信息的前提下，构建了复合多层次的景观阻力评价体系，基于最小费用路径法和重力模型，提取出浦东新区重要生态廊道和生态节点，在此基础上对浦东新区生态网络结构进行了分析评价．

遥感和GIS技术在本文的研究中起到重要作用，高精度遥感影像与野外实地调查相结合的方法，提供了详尽而又准确的数据支撑，GIS强大的空间分析功能则为生态网络结构的评价分析提供了较好的实现手段．

在构建景观阻力评价体系时，虽然考虑了湿地、林绿地等不同景观在面积、植被覆盖度、保护程度等多因子不同属性值，能体现不同属性对景观阻力的影响，但由于缺少影响因子与景观阻力精确的数字关系，在系数选取时主要以专家知识为主．生态网络的构建应该是基于多目的（如物种迁移、生态保护、居民休闲等）、多尺度（乡镇级、区级、省级或者更大），对数据的详细程度和指标体系的选取有着更高的要求，未来研究将进一步深入和强化．

由于浦东新区城市化率较高，同时经济发展有着强烈诉求，用地极为紧张，因此进行大规模生态用地建设是不实际的，唯有加强结构优化调整，充分利用原有生态资源，加强生境斑块间的链接和功能构建多层次多维度的城市生态网络，以此发挥最大生态服务功能．

［1］FREEMUTH J C．Islands under Siege：National Parks and the Politics of External Threats［M］．Lawrence：Uni－versity Press of Kansas，1991．

［2］COLLINGE S K．Ecological consequences of habitat fragmentation：implications for landscape architecture and planning［J］．Landscape and UrbanPlanning，1996，36（1）：59－97．

［3］ROB J，GLORIA P．生态网络与绿道：概念，设计与实施［M］．余青，陈海沐，梁莺莺译．北京：中国建筑工业出版社，2011．

［4］陆同伟，宋珂，杨秀，等．基于生态适宜度分析的城市用地规划研究——以杭州市东南部生态带保护和控制规划为例［J］．复旦学报：自然科学版，2011，50（2）：245－251．

［5］BASCHAK L A，BROWN R D．An ecological framework for the planning，design and management of urban river greenways［J］．Landscape and Urban Planning，1995，33：211－225．

［6］王海珍，张利权．基于GIS、景观格局和网络分析法的厦门本岛生态网络规划［J］．植物生态学报，2005（1）：144－152．

［7］王祥荣．论生态城市建设的理论、途径与措施——以上海为例［J］．复旦学报：自然科学版，2001，40（4）：349－354．

［8］鲁敏，杨东兴，刘佳，等．济南绿地生态网络体系的规划布局与构建［J］．中国生态农业学报，2010（3）：600－605．

［9］陈爽，王进，詹志勇．生态景观与城市形态整合研究［J］．地理科学进展，2004，23（5）：67－77．

［10］JACK A．Greenways as a planning strategy［J］．Landscape and Urban Planning，1995，33：131－155．

［11］王云才，刘悦来．城市景观生态网络规划的空间模式应用探讨［J］．长江流域资源与环境，2009（9）：819－824．

［12］王云才．上海市城市景观生态网络连接度评价［J］．地理研究，2009（2）：284－292．

［13］郭微，俞龙生，孙延军，等．佛山市顺德中心城区城市绿地生态网络规划［J］．生态学杂志，2012（4）：1022－1027．

［14］孔繁花，尹海伟．济南城市绿地生态网络构建［J］．生态学报，2008（4）：1711－1719．

［15］傅强，宋军，毛锋，等．青岛市湿地生态网络评价与构建［J］．生态学报，2012（12）：3670－3680．

［16］韩英．经济全球化对浦东新区环境保护的压力与对策［J］．上海环境科学，2003（2）：104－109．

［17］马凤，蔡永立．上海临港新城开发对土地利用及生态系统服务价值的影响［J］．复旦学报：自然科学版，2013，52（1）：69－77．

［18］FORMANR T T．Corridors in a landscape：their ecological structure and function［J］．Ekologia，1983，2（4）：375－387．

［19］DAUBER J，HIRSCH M，SIMMERING D，et al．Landscape structure as an indicator of biodiversity：matrix effects on species richness［J］．AgricultureEcosystems and Environment，2003，98：321－329．

［20］高赛男．基于最小路径方法的南京城市绿地网络构建［J］．山东师范大学学报：自然科学版，2009，24（3）：90－93．

［21］KNAAPEN J P，SCHEFFER M，HARMS B．Estimating habitat isolation in landscape［J］．Landscape and Urban Planning，1992，23（1）：1－16．

［22］赵琳，田永中．三峡库区城镇建设用地生态适宜性评价——以涪陵新区为例［C］．山地环境与生态文明建设——中国地理学会2013年学术年会西南片区会议论文集，2013．

［23］陈元涛．基于GIS的道路智能化选线方法研究［D］．重庆：重庆交通大学，2012．

［24］蔡青．基于景观生态学的城市空间格局演变规律分析与生态安全格局构建［D］．湖南：湖南大学，2012．

［25］尹发能，王学雷．基于最小累计阻力模型的四湖流域景观生态规划研究［J］．华中农业大学学报，2010（2）：231－235．

［26］尹海伟，孔繁花，祈毅，等．湖南省城市群生态网络构建与优化［J］．生态学报，2011（10）：2863－2874．