文章信息
- 石丹, 关婧文, 刘吉平
- SHI Dan, GUAN Jingwen, LIU Jiping
- 基于DPSIR-EES模型的旅游型城镇生态安全评价研究
- Ecological security evaluation of tourism towns based on DPSIR-EES-matter element
- 生态学报. 2021, 41(11): 4330-4341
- Acta Ecologica Sinica. 2021, 41(11): 4330-4341
- http://dx.doi.org/10.5846/stxb202012253275
-
文章历史
- 收稿日期: 2020-12-25
- 修订日期: 2021-05-11
生态安全是指区域内生态系统的完整情况和健康水平, 其本质在于强调生态系统的自身结构是否受到破坏, 同时强调生态系统对人类生活所提供的服务是否得到满足[1], 即区域内自然生态环境状况能够维系其经济、社会可持续发展时, 这种状态下生态处于安全水平, 反之则处于不安全水平[2]。随经济、社会发展进入纵深阶段, 旅游业在国民经济中的比重越来越大, 旅游型城镇建设日趋成熟, 但在加快发展过程中粗放式、无节制地开发, 使得旅游型城镇生态环境保护形势十分严峻, 已经逼近旅游型城镇生态安全承载力阈值, 资源环境持续恶化与生态系统可持续发展的目标不相符[3], 严重影响旅游型城镇生态系统的自我完整性及可持续发展, 还威胁其生态安全[4], 如何实现旅游型城镇发展与资源环境的全面脱钩是现阶段亟待解决的重要问题, 生态安全问题引发业界对旅游型城镇生态安全和环境保护的理性思考。对于生态安全的研究, 大多数学者主要集中在生态安全的概念和内涵[5]、评价指标模型[6]和格局构建与优化[7]等方面的研究, 评价方法涉及模糊分析法[8]、层次分析法[9]、TOPSIS法[10]、灰色聚类分析法[11]等, 指标构建主要选取PSR模型[12-13]、DPSIR模型[14-16]、EES模型[17]等。研究尺度大多为省域[18]、市域[19]或区域[20], 研究分析方法主要集中在生态足迹法[21]、正态云模型[22]、ARIMA模型[23]、物元模型[24-27]、GIS技术[28-29]等。
对于旅游型城镇生态安全评价分析而言, 仍存在以下三点问题亟待解决:(1)现有研究针对生态安全研究领域, 大多数采用中、大尺度, 虽从不同角度丰富和完善了生态安全的研究体系, 但缺乏通过从单一旅游型城镇小尺度视角进行生态安全评价分析; (2)大部分生态安全评价体系多以单一模型构建为主, 其只能反映各子系统间的相互关系, 不能更好地实现各子系统内部相互作用关系; (3)多数生态安全评价模型仅考虑各评价指标之间的相互影响, 未考虑各评价指标因素间具有不相容性且存在交叉作用, 忽视评价过程中信息屏蔽和主观性等问题。针对上述问题, 基于区域可持续发展理论, 运用DPSIR-EES复合模型进一步完善生态安全评价指标体系, 以二道白河镇为例, 通过熵值法和均方差决策法对数据进行处理并赋权, 引入物元分析模型, 评价2010-2019年生态安全水平, 分析各指标之间的生态安全水平变化趋势, 以期更好地实现二道白河镇生态系统可持续发展, 对其自身及同类型旅游型城镇改善生态安全环境质量、促进生态系统良性循环、实现经济社会环境可持续发展具有指导意义。
1 研究区域概况与数据来源 1.1 研究区域概况二道白河镇位于吉林省东南部, 地理位置处于东经127°53′-128°34′、北纬42°01′-42°48′之间, 地处长白山脚下, 使得二道白河镇旅游业发展独具优势[30]:(1)境内林地面积272883hm2, 森林覆盖率达94%, 拥有丰富的森林资源; (2)境内除长白山外, 地势开阔、起伏较小, 草原面积12342hm2, 适合人类开展各类活动; (3)区域内有120余种植物种类, 其中30多种树种属于经济价值高的种类, 具有丰富的自然生态资源。独特的地理位置和丰富自然生态资源为二道白河镇的生态环境发展奠定了良好的基础, 同时为该区域生态环境建设发挥了极其重要的作用。
1.2 数据来源生态安全评价指标可分为原始指标数据以及构建指标数据, 选择2010-2019年作为研究期。原始指标数据主要通过查阅《中国区域经济统计年鉴(2011-2020年)》、《吉林统计年鉴(2011-2020年)》、《长白山保护开发区年鉴(2014-2019年)》, 其中个别数据来源于以下三个方面:(1)《延边年鉴(2011-2015年)》、《延吉年鉴(2011-2019年)》; (2)吉林省延边州、安图县市(州)旅游局、统计局等官方网站; (3)吉林省国民经济和社会发展统计公报和安图县、长白山保护区其他相关政府网站等有关旅游业发展的政府工作报告。构建指标数据通过原始指标数据复合计算得出, 例如人均生产总值=生产总值/人口数。对于个别缺失数据, 根据已有数据通过算法计算进行合理性赋值。
2 旅游型城镇生态安全物元评价模型 2.1 DPSIR-EES生态安全评价模型框架由“驱动力(D)、压力(P)、状态(S)、影响(I)、响应(R)”五个方面所构成的DPSIR模型[14-16], 用于表达影响生态安全整个系统各因素之间的信息耦合关系; 由“环境(E)、经济(E)、社会(S)”三个层面所构成的EES模型[17], 用于描述整个人类活动系统中各属性之间的协同关系。DPSIR模型能够较好地反映生态各系统之间的因果影响, 但无法实现各系统内部的分析; EES模型能够较好地反应生态系统内部复合属性, 进行系统内部结构分析, 但无法反映系统内因果关系。利用两个单项模型的优点实现有机结合, 形成DPSIR-EES模型(图 1), 将各要素有机结合, 能更清晰地反映之间相互作用关系。
![]() |
图 1 二道白河镇生态安全评价DPSIR-EES模型框架 Fig. 1 DPSIR-EES model framework for ecological security assessment in Erdaobaihe town |
初步选取45项具体评价指标, 为保证指标体系的科学性, 参考SDGs(可持续发展目标指标)[31]、EPI(环境绩效评估指标)[32]以及综合相关研究成果[2-3, 9, 16, 25-27], 运用SPSS22.0和Delphi法进行主成分分析和因子分析, 最终确定36个指标(表 1)。明确各指标属性, 分为效益型(+)或成本型(-), 结合熵权法和均方差决策法综合计算得出更为精确的综合权重(Ci)。
目标层 Target layer |
准则层 Criteria layer |
要素层 Element layer |
指标层 Index layer |
指标符号 Indicator symbol |
属性 Attribute |
熵值权重 Entropy weight |
均方差权重 Mean square error weight |
综合权重 Comprehensive weight |
旅游型城镇生态安全 | 驱动力 | 环境ED1 | 森林覆盖率/% | ED11 | + | 0.0291 | 0.0242 | 0.0266 |
Ecological security of | Driving force | 建成区绿化率/% | ED12 | + | 0.0285 | 0.0254 | 0.0269 | |
tourist towns(ES) | (D) | 空气质量优良天数水平/% | ED13 | + | 0.0275 | 0.0252 | 0.0263 | |
经济ED2 | 经济发展水平/亿元 | ED21 | + | 0.0289 | 0.0238 | 0.0264 | ||
人均生产总值/元 | ED22 | + | 0.0287 | 0.0245 | 0.0266 | |||
旅游业发展水平/% | ED23 | + | 0.0258 | 0.0339 | 0.0299 | |||
社会SD3 | 人口密度/(人/km2) | SD31 | - | 0.0286 | 0.0266 | 0.0276 | ||
人口自然增长率/% | SD32 | - | 0.0263 | 0.0297 | 0.0280 | |||
城镇化水平/% | SD33 | + | 0.0287 | 0.0286 | 0.0287 | |||
压力 | 环境EP1 | 农药化肥使用量/万吨 | EP11 | - | 0.0281 | 0.0271 | 0.0276 | |
Pressure(P) | 工业固体废物综合水平/% | EP12 | + | 0.0278 | 0.0295 | 0.0287 | ||
经济EP2 | GDP增长率/% | EP21 | + | 0.0284 | 0.0255 | 0.0269 | ||
人均可支配收入/元 | EP22 | + | 0.0285 | 0.0263 | 0.0274 | |||
社会SP3 | 社会保障建设水平/% | SP31 | + | 0.0280 | 0.0298 | 0.0289 | ||
城镇登记失业率/% | SP32 | - | 0.0273 | 0.0317 | 0.0295 | |||
状态State(S) | 环境ES1 | 人均水资源占有量/(m3/人) | ES11 | + | 0.0274 | 0.0287 | 0.0281 | |
旅游资源禀赋情况/个 | ES12 | + | 0.0290 | 0.0247 | 0.0269 | |||
经济ES2 | 经济密度/(亿元/hm2) | ES21 | + | 0.0289 | 0.0238 | 0.0264 | ||
人均固定资产投资额/元 | ES22 | + | 0.0279 | 0.0291 | 0.0285 | |||
社会SS3 | 人均消费能力/元 | SS31 | + | 0.0285 | 0.0282 | 0.0283 | ||
人均建设用地面积/hm2 | SS32 | + | 0.0285 | 0.0258 | 0.0271 | |||
影响Impact(I) | 环境EI1 | 水土流失率/% | EI11 | - | 0.0273 | 0.0297 | 0.0285 | |
清洁能源使用率/% | EI12 | + | 0.0279 | 0.0303 | 0.0291 | |||
经济EI2 | 第三产业比重/% | EI21 | + | 0.0259 | 0.0332 | 0.0295 | ||
旅游接待能力/万人次 | EI22 | + | 0.0275 | 0.0265 | 0.0270 | |||
社会SI3 | 就业与经济增长关联度/% | SI31 | + | 0.0284 | 0.0304 | 0.0294 | ||
科技进步对经济贡献率/% | SI32 | + | 0.0275 | 0.0271 | 0.0273 | |||
响应 | 环境ER1 | 生活污水处理能力/% | ER11 | + | 0.0290 | 0.0274 | 0.0282 | |
Response(R) | 生活垃圾处理能力/% | ER12 | + | 0.0289 | 0.0279 | 0.0284 | ||
土地退化治理率/% | ER13 | + | 0.0284 | 0.0265 | 0.0274 | |||
经济ER2 | 社会消费品零售总额/亿元 | ER21 | + | 0.0289 | 0.0277 | 0.0283 | ||
节能环保力度/% | ER22 | + | 0.0268 | 0.0270 | 0.0269 | |||
财政支出总额/亿元 | ER23 | + | 0.0252 | 0.0294 | 0.0273 | |||
社会SR3 | 公共设施建设水平/% | SR31 | + | 0.0280 | 0.0307 | 0.0293 | ||
居民人均居住用地/(m2/人) | SR32 | + | 0.0255 | 0.0272 | 0.0264 | |||
区域开发比例/% | SR33 | + | 0.0245 | 0.0271 | 0.0258 |
旅游型城镇生态安全N, 生态安全特征w和特征量值v共同构成生态安全物元。假设生态安全N中有多个特征值则表示为:
![]() |
式中, R为生态安全的n维物元, 可表示为R=(N, w, v)。
对于二道白河镇生态安全综合性评价物元表示为:
![]() |
对于二道白河镇生态安全准则层评价, 如驱动力环境层面ED1生态安全水平评价的物元表示为:
![]() |
旅游型城镇生态安全的经典域物元Roj=(Noj, Wi, Vo), 其矩阵表示为:
![]() |
式中, Roj为生态安全经典域物元, Noj为第j个生态安全评价等级, Wi为第i个生态安全评价指标, (aoji, boji)为生态安全对应评价等级j的量值范围, 即旅游型城镇生态安全评价经典域。
旅游型城镇生态安全的节域物元Rp=(Np, Wn, Vp), 其矩阵表示为:
![]() |
式中, Rp为生态安全节域物元, Np为生态安全评价等级, (apn, bpn)为生态安全节域物元对应Wn的取值范围, 即旅游型城镇生态安全评价节域。
旅游型城镇生态安全待评价对象Nx的物元可表示为RX=(N, W, V), 其矩阵表示为:
![]() |
旅游型城镇生态安全关联函数K(x)数值表示评价对象与某一标准范围的符合程度, 物元模型的逻辑值从模糊数学的[0, 1]取值范围扩展到(-∞, +∞)实数轴, 在实数轴上的数值表示评价对象与某标准范围的隶属程度:K(x)<-1, 则表示被评价对象不符合标准对象的要求且不具备转化对标准对象的条件; -1≤K(x)<0, 则表示被评价对象不符合标准对象的要求, 但具备转化为标准对象的条件; 0≤K(x)<1, 则表示被评价对象符合标准对象的要求, 数值越大则表示越接近标准上限; K(x)≥1, 则表示被评价对象超过标准对象的要求上限, 数值越大则表示具有越大的开发潜力[26-27, 33-35]。
旅游型城镇生态安全评价指标关联度函数K(x)表示为:
![]() |
其中实变函数距离表示为:
![]() |
式中, ρ(X, Xo)为点X与有界区间Xo=[ao, bo]的距离, ρ(X, Xp)为点X与对应节域的有界区间Xp=[ap, bp]的距离, X、Xo、Xp为旅游型城镇生态安全评价物元的量值、经典域及节域物元量值范围。
2.3.4 计算综合关联度并确定评价等级旅游型城镇生态安全待评价对象Nx关于等级j的综合关联度Kj(Nx)表示为:
![]() |
式中, Kj(xi)为各评价等级水平, Ci为第i个指标的权重。确定评价对象归属等级标准采用最大值原则Kj=maxKj(Nx), 则可评定评价对象Nx属于对应评价等级标准j级。
2.4 生态安全评价经典域、节域的确定确定经典域是物元评价的基础, 依据生态安全综合属性, 参照国家环保总局《生态县、生态市、生态省建设指标(试行)》标准值、国际通行标准、全国及区域平均水平及同类相关指标等标准值, 参考相关研究[2-6, 26-27, 30-35], 将二道白河生态安全等级分为五个等级:Ⅴ级(安全)、Ⅳ级(较安全)、Ⅲ级(临界安全)、Ⅱ级(较不安全)、Ⅰ级(不安全)。评价主要标准从以下方面选取:(1)国际、国家、地方所规定的标准, 如ED1、EP1、EI1、ER1等与生态环境质量标准为主要统计指标, 采取对应标准值划分等级; (2)类比同类标准, 如ED2、SD3、EP2、SP3、ES1、ES2、SS3、EI2、SI3、ER2、SR3等与经济发展、社会生活为主要统计指标, 参考国内外具有良好特色的同类区域现状值作为标准值或趋势值外推设定“Ⅴ级(安全)”, 并以国内外同类区域现状较劣值或国际公认的警戒线指标设定“Ⅰ级(不安全)”, 最后综合二道白河镇实际情况, 得出该区域生态安全评价经典域复合物元量值范围以及节域物元量值范围评价标准(见表 2), 通过关联度函数计算判断待评价对象所属等级。
目标层 Target layer |
准则层 Criteria Layer |
要素层 Element layer |
指标层 Index layer |
评价标准(经典域)Evaluation criteria (classical domain) | 节域 Joint domains |
||||
Ⅴ-安全 Safety |
Ⅳ-较安全 Relatively safety |
Ⅲ-临界安全 Critical safety |
Ⅱ-较不安全 Less safety |
Ⅰ-不安全 Unsafety |
|||||
旅游型城镇生态安全 | 驱动力(D) | 环境ED1 | 森林覆盖率ED11 | 0.8-1 | 0.6-0.8 | 0.5-0.6 | 0.3-0.5 | 0-0.3 | 0-1 |
Ecological security of | 建成区绿化率ED12 | 0.75-1 | 0.6-0.75 | 0.45-0.6 | 0.25-0.45 | 0-0.25 | 0-1 | ||
tourist towns(ES) | 空气质量优良天数水平ED13 | 0.8-1 | 0.6-0.8 | 0.4-0.6 | 0.3-0.4 | 0-0.3 | 0-1 | ||
经济ED2 | 经济发展水平ED21 | 0.7-1 | 0.5-0.7 | 0.4-0.5 | 0.2-0.4 | 0-0.2 | 0-1 | ||
人均生产总值ED22 | 0.8-1 | 0.7-0.8 | 0.5-0.7 | 0.25-0.5 | 0-0.25 | 0-1 | |||
旅游业发展水平ED23 | 0.7-1 | 0.5-0.7 | 0.4-0.5 | 0.2-0.4 | 0-0.2 | 0-1 | |||
社会SD3 | 人口密度SD31 | 0.8-1 | 0.65-0.8 | 0.5-0.65 | 0.35-0.5 | 0-0.35 | 0-1 | ||
人口自然增长率SD32 | 0.75-1 | 0.6-0.75 | 0.45-0.6 | 0.25-0.45 | 0-0.25 | 0-1 | |||
城镇化水平SD33 | 0.8-1 | 0.65-0.8 | 0.45-0.65 | 0.3-0.45 | 0-0.3 | 0-1 | |||
压力(P) | 环境EP1 | 农药化肥使用量EP11 | 0.8-1 | 0.6-0.8 | 0.4-0.6 | 0.2-0.4 | 0-0.2 | 0-1 | |
工业固体废物综合水平EP12 | 0.85-1 | 0.7-0.85 | 0.5-0.7 | 0.35-0.5 | 0-0.35 | 0-1 | |||
经济EP2 | GDP增长率EP21 | 0.7-1 | 0.6-0.7 | 0.45-0.6 | 0.2-0.45 | 0-0.2 | 0-1 | ||
人均可支配收入EP22 | 0.8-1 | 0.6-0.8 | 0.4-0.6 | 0.2-0.4 | 0-0.2 | 0-1 | |||
社会SP3 | 社会保障建设水平SP31 | 0.85-1 | 0.7-0.85 | 0.5-0.7 | 0.25-0.5 | 0-0.25 | 0-1 | ||
城镇登记失业率SP32 | 0.8-1 | 0.6-0.8 | 0.4-0.6 | 0.2-0.4 | 0-0.2 | 0-1 | |||
状态(S) | 环境ES1 | 人均水资源占有量ES11 | 0.8-1 | 0.7-0.8 | 0.5-0.7 | 0.25-0.5 | 0-0.25 | 0-1 | |
旅游资源禀赋情况ES12 | 0.85-1 | 0.65-0.85 | 0.45-0.65 | 0.25-0.45 | 0-0.25 | 0-1 | |||
经济ES2 | 经济密度ES21 | 0.75-1 | 0.6-0.75 | 0.45-0.6 | 0.2-0.45 | 0-0.2 | 0-1 | ||
人均固定资产投资额ES22 | 0.8-1 | 0.6-0.8 | 0.4-0.6 | 0.2-0.4 | 0-0.2 | 0-1 | |||
社会SS3 | 人均消费能力SS31 | 0.8-1 | 0.6-0.8 | 0.4-0.6 | 0.15-0.4 | 0-0.15 | 0-1 | ||
人均建设用地面积SS32 | 0.75-1 | 0.65-0.75 | 0.45-0.65 | 0.3-0.45 | 0-0.3 | 0-1 | |||
影响(I) | 环境EI1 | 水土流失率EI11 | 0.8-1 | 0.6-0.8 | 0.3-0.6 | 0.1-0.3 | 0-0.1 | 0-1 | |
清洁能源使用率EI12 | 0.75-1 | 0.6-0.75 | 0.4-0.6 | 0.25-0.4 | 0-0.25 | 0-1 | |||
经济EI2 | 第三产业比重EI21 | 0.9-1 | 0.8-0.9 | 0.5-0.8 | 0.1-0.5 | 0-0.1 | 0-1 | ||
旅游接待能力EI22 | 0.8-1 | 0.65-0.8 | 0.4-0.65 | 0.25-0.4 | 0-0.25 | 0-1 | |||
社会SI3 | 就业与经济增长关联度SI31 | 0.85-1 | 0.7-0.85 | 0.6-0.7 | 0.4-0.6 | 0-0.4 | 0-1 | ||
科技进步对经济贡献率SI32 | 0.85-1 | 0.65-0.85 | 0.45-0.65 | 0.25-0.45 | 0-0.25 | 0-1 | |||
响应(R) | 环境ER1 | 生活污水处理能力ER11 | 0.9-1 | 0.8-0.9 | 0.6-0.8 | 0.3-0.6 | 0-0.3 | 0-1 | |
生活垃圾处理能力ER12 | 0.85-1 | 0.7-0.85 | 0.6-0.7 | 0.4-0.6 | 0-0.4 | 0-1 | |||
土地退化治理率ER13 | 0.8-1 | 0.6-0.8 | 0.4-0.6 | 0.2-0.4 | 0-0.2 | 0-1 | |||
经济ER2 | 社会消费品零售总额ER21 | 0.85-1 | 0.65-0.85 | 0.5-0.65 | 0.3-0.5 | 0-0.3 | 0-1 | ||
节能环保力度ER22 | 0.8-1 | 0.6-0.8 | 0.4-0.6 | 0.25-0.4 | 0-0.25 | 0-1 | |||
财政支出总额ER23 | 0.8-1 | 0.6-0.8 | 0.3-0.6 | 0.1-0.3 | 0-0.1 | 0-1 | |||
社会SR3 | 公共设施建设水平SR31 | 0.9-1 | 0.75-0.9 | 0.45-0.75 | 0.2-0.45 | 0-0.2 | 0-1 | ||
居民人均居住用地SR32 | 0.65-1 | 0.3-0.65 | 0.25-0.3 | 0.1-0.25 | 0-0.1 | 0-1 | |||
区域开发比例SR33 | 0.8-1 | 0.6-0.8 | 0.2-0.6 | 0.1-0.2 | 0-0.1 | 0-1 |
计算得2010-2019年二道白河镇生态安全评价指标关联度(表 3)及综合关联度(表 4)结果以及对应安全等级。
指标及年份 Index and yea |
2010 | 2010 | 2011 | 2012 | 2013 | 2014 | 2015 | 2016 | 2017 | 2018 | 2019 | |||||
Ⅴ-安全 Safety |
Ⅳ-较安全 Relatively safety |
Ⅲ-临界安全 Critical safety |
Ⅱ-较不安全 Less safety |
Ⅰ-不安全 Unsafety |
等级 Grade | |||||||||||
kj(ED11) | 0.0333 | -0.9667 | -0.9800 | -0.9833 | -0.9875 | Ⅰ | Ⅱ | Ⅱ | Ⅱ | Ⅱ | Ⅲ | Ⅳ | Ⅳ | Ⅴ | Ⅴ | |
kj(ED12) | -0.1699 | 0.3217 | -0.3015 | -0.4761 | -0.5809 | Ⅱ | Ⅰ | Ⅱ | Ⅲ | Ⅲ | Ⅳ | Ⅳ | Ⅰ | Ⅴ | Ⅴ | |
kj(ED13) | 0.0333 | -0.9667 | -0.9750 | -0.9833 | -0.9875 | Ⅰ | Ⅰ | Ⅱ | Ⅰ | Ⅱ | Ⅱ | Ⅴ | Ⅴ | Ⅰ | Ⅱ | |
kj(ED21) | 0.0500 | -0.9667 | -0.9750 | -0.9800 | -0.9857 | Ⅰ | Ⅱ | Ⅱ | Ⅲ | Ⅳ | Ⅳ | Ⅳ | Ⅴ | Ⅴ | Ⅴ | |
kj(ED22) | 0.0400 | -0.9600 | -0.9800 | -0.9857 | -0.9875 | Ⅰ | Ⅰ | Ⅱ | Ⅱ | Ⅲ | Ⅲ | Ⅲ | Ⅲ | Ⅴ | Ⅴ | |
kj(ED23) | 0.4348 | -0.5652 | -0.7826 | -0.8260 | -0.8758 | Ⅰ | Ⅰ | Ⅰ | Ⅰ | Ⅱ | Ⅱ | Ⅳ | Ⅴ | Ⅴ | Ⅴ | |
kj(SD31) | -0.1185 | 0.3624 | -0.1913 | -0.3779 | -0.4946 | Ⅱ | Ⅰ | Ⅰ | Ⅱ | Ⅱ | Ⅲ | Ⅳ | Ⅴ | Ⅴ | Ⅴ | |
kj(SD32) | 0.0400 | -0.9600 | -0.9778 | -0.9833 | -0.9867 | Ⅰ | Ⅰ | Ⅰ | Ⅰ | Ⅰ | Ⅱ | Ⅱ | Ⅴ | Ⅴ | Ⅴ | |
kj(SD33) | 0.0333 | -0.9667 | -0.9778 | -0.9846 | -0.9875 | Ⅰ | Ⅲ | Ⅱ | Ⅱ | Ⅰ | Ⅳ | Ⅴ | Ⅴ | Ⅴ | Ⅴ | |
kj(EP11) | 0.0500 | -0.9500 | -0.9750 | -0.9833 | -0.9875 | Ⅰ | Ⅰ | Ⅱ | Ⅱ | Ⅲ | Ⅲ | Ⅳ | Ⅳ | Ⅴ | Ⅴ | |
kj(EP12) | 0.0286 | -0.9714 | -0.9800 | -0.9857 | -0.9882 | Ⅰ | Ⅰ | Ⅰ | Ⅰ | Ⅲ | Ⅳ | Ⅳ | Ⅳ | Ⅴ | Ⅴ | |
kj(EP21) | -0.8167 | -0.7334 | -0.6335 | -0.5113 | 0.4887 | Ⅴ | Ⅴ | Ⅴ | Ⅱ | Ⅱ | Ⅰ | Ⅰ | Ⅲ | Ⅱ | Ⅲ | |
kj(EP22) | 0.0500 | -0.9500 | -0.9750 | -0.9833 | -0.9875 | Ⅰ | Ⅰ | Ⅱ | Ⅲ | Ⅲ | Ⅲ | Ⅳ | Ⅴ | Ⅴ | Ⅴ | |
kj(SP31) | -0.7588 | -0.6381 | -0.3969 | 0.2062 | -0.1460 | Ⅳ | Ⅴ | Ⅴ | Ⅲ | Ⅴ | Ⅳ | Ⅰ | Ⅲ | Ⅰ | Ⅰ | |
kj(SP32) | 0.0500 | -0.9500 | -0.9750 | -0.9833 | -0.9875 | Ⅰ | Ⅰ | Ⅰ | Ⅱ | Ⅲ | Ⅳ | Ⅳ | Ⅴ | Ⅴ | Ⅴ | |
kj(ES11) | 0.0400 | -0.9600 | -0.9800 | -0.9857 | -0.9875 | Ⅰ | Ⅰ | Ⅰ | Ⅱ | Ⅱ | Ⅱ | Ⅲ | Ⅳ | Ⅴ | Ⅴ | |
kj(ES12) | 0.0400 | -0.9600 | -0.9778 | -0.9846 | -0.9882 | Ⅰ | Ⅰ | Ⅱ | Ⅱ | Ⅲ | Ⅲ | Ⅳ | Ⅳ | Ⅴ | Ⅴ | |
kj(ES21) | 0.0500 | -0.9500 | -0.9778 | -0.9833 | -0.9867 | Ⅰ | Ⅱ | Ⅱ | Ⅲ | Ⅲ | Ⅲ | Ⅲ | Ⅳ | Ⅴ | Ⅴ | |
kj(ES22) | 0.4055 | -0.4055 | -0.7028 | -0.8018 | -0.8514 | Ⅰ | Ⅰ | Ⅰ | Ⅲ | Ⅱ | Ⅳ | Ⅳ | Ⅴ | Ⅴ | Ⅴ | |
kj(SS31) | 0.0667 | -0.9333 | -0.9778 | -0.9833 | -0.9867 | Ⅰ | Ⅰ | Ⅱ | Ⅲ | Ⅳ | Ⅳ | Ⅳ | Ⅴ | Ⅴ | Ⅴ | |
kj(SS32) | 0.0333 | -0.9667 | -0.9778 | -0.9846 | -0.9867 | Ⅰ | Ⅱ | Ⅰ | Ⅲ | Ⅲ | Ⅲ | Ⅲ | Ⅲ | Ⅴ | Ⅴ | |
kj(EI11) | 0.1000 | -0.9000 | -0.9667 | -0.9833 | -0.9875 | Ⅰ | Ⅰ | Ⅱ | Ⅱ | Ⅲ | Ⅳ | Ⅲ | Ⅳ | Ⅴ | Ⅴ | |
kj(EI12) | 0.0400 | -0.9667 | -0.9750 | -0.9833 | -0.9867 | Ⅰ | Ⅰ | Ⅰ | Ⅲ | Ⅴ | Ⅴ | Ⅴ | Ⅴ | Ⅰ | Ⅳ | |
kj(EI21) | 0.3616 | -0.6385 | -0.9277 | -0.9548 | -0.9598 | Ⅰ | Ⅰ | Ⅰ | Ⅰ | Ⅱ | Ⅳ | Ⅴ | Ⅳ | Ⅳ | Ⅲ | |
kj(EI22) | 0.0400 | -0.9600 | -0.9750 | -0.9846 | -0.9875 | Ⅰ | Ⅰ | Ⅰ | Ⅱ | Ⅱ | Ⅱ | Ⅲ | Ⅳ | Ⅴ | Ⅴ | |
kj(SI31) | 0.2055 | -0.7945 | -0.8630 | -0.8826 | -0.9033 | Ⅰ | Ⅰ | Ⅰ | Ⅲ | Ⅳ | Ⅴ | Ⅳ | Ⅳ | Ⅴ | Ⅴ | |
kj(SI32) | 0.0400 | -0.9600 | -0.9778 | -0.9846 | -0.9882 | Ⅰ | Ⅰ | Ⅰ | Ⅱ | Ⅱ | Ⅲ | Ⅳ | Ⅲ | Ⅳ | Ⅴ | |
kj(ER11) | -0.2291 | 0.4228 | -0.2886 | -0.4665 | -0.5257 | Ⅱ | Ⅳ | Ⅲ | Ⅲ | Ⅰ | Ⅳ | Ⅴ | Ⅰ | Ⅴ | Ⅲ | |
kj(ER12) | 0.0250 | -0.9750 | -0.9833 | -0.9857 | -0.9882 | Ⅰ | Ⅰ | Ⅴ | Ⅴ | Ⅳ | Ⅳ | Ⅳ | Ⅲ | Ⅲ | Ⅴ | |
kj(ER13) | 0.0500 | -0.9500 | -0.9750 | -0.9833 | -0.9875 | Ⅰ | Ⅰ | Ⅱ | Ⅱ | Ⅲ | Ⅲ | Ⅳ | Ⅳ | Ⅴ | Ⅴ | |
kj(ER21) | 0.0333 | -0.9667 | -0.9800 | -0.9846 | -0.9882 | Ⅰ | Ⅰ | Ⅱ | Ⅲ | Ⅳ | Ⅳ | Ⅳ | Ⅴ | Ⅴ | Ⅴ | |
kj(ER22) | -0.6897 | -0.6122 | -0.4183 | 0.1635 | -0.1232 | Ⅳ | Ⅴ | Ⅳ | Ⅲ | Ⅱ | Ⅰ | Ⅰ | Ⅱ | Ⅰ | Ⅰ | |
kj(ER23) | 0.1000 | -0.9000 | -0.9667 | -0.9833 | -0.9875 | Ⅰ | Ⅰ | Ⅰ | Ⅱ | Ⅱ | Ⅱ | Ⅲ | Ⅲ | Ⅴ | Ⅴ | |
kj(SR31) | 0.0500 | -0.9500 | -0.9778 | -0.9867 | -0.9889 | Ⅰ | Ⅰ | Ⅱ | Ⅲ | Ⅲ | Ⅲ | Ⅳ | Ⅳ | Ⅴ | Ⅴ | |
kj(SR32) | 0.0561 | -0.0561 | -0.6224 | -0.6854 | -0.8548 | Ⅰ | Ⅱ | Ⅳ | Ⅰ | Ⅰ | Ⅱ | Ⅱ | Ⅲ | Ⅴ | Ⅴ | |
kj(SR33) | 0.1000 | -0.9000 | -0.9500 | -0.9833 | -0.9875 | Ⅰ | Ⅰ | Ⅰ | Ⅱ | Ⅱ | Ⅱ | Ⅱ | Ⅱ | Ⅳ | Ⅴ |
Ⅰ-不安全 | Ⅱ-较不安全 | Ⅲ-临界安全 | Ⅳ-较安全 | Ⅴ-安全 | 等级 | |
kj(w2010) | -0.0005 | -0.7542 | -0.8479 | -0.8394 | -0.8486 | 向Ⅰ级转化 |
kj(w2011) | 0.0455 | -0.5376 | -0.7038 | -0.7753 | -0.7607 | Ⅰ级 |
kj(w2012) | -0.0948 | -0.1025 | -0.4256 | -0.5478 | -0.6178 | 向Ⅰ级转化 |
kj(w2013) | -0.2255 | -0.0496 | -0.1157 | -0.3911 | -0.5188 | 向Ⅱ级转化 |
kj(w2014) | -0.3213 | -0.1241 | -0.1240 | -0.2744 | -0.4529 | 向Ⅲ级转化 |
kj(w2015) | -0.4788 | -0.2539 | -0.1428 | -0.1676 | -0.3685 | 向Ⅲ级转化 |
kj(w2016) | -0.5124 | -0.4068 | -0.2533 | -0.1075 | -0.2811 | 向Ⅳ级转化 |
kj(w2017) | -0.6087 | -0.5297 | -0.3635 | -0.1746 | -0.2052 | 向Ⅳ级转化 |
kj(w2018) | -0.7243 | -0.7395 | -0.6578 | -0.4826 | 0.0484 | Ⅴ级 |
kj(w2019) | -0.8212 | -0.8243 | -0.8005 | -0.8523 | -0.0476 | 向Ⅴ级转化 |
从表 3可知, 在2010-2019年ED11、ED21、EP11、EP22、SP32、ES11、ES12、ES21、ES22、SS31、EI11、EI22、ER13、ER21、SR31均呈现出不同程度的从Ⅰ级(不安全)向Ⅴ级(安全)逐年转化的上升趋势, 说明以上指标对于二道白河镇生态安全水平的提升具有积极作用。而ED13、EP21、SP31、SI31、ER11、ER22在2019年仍处于Ⅰ级(不安全)、Ⅱ级(较不安全)或Ⅲ级(临界安全)水平, 且呈现明显波动变化, 其中EP21、SP31、ER22呈现逐年下降的趋势, 已成为制约二道白河镇生态安全进一步提升的主要因素。
由表 4看出, 2010-2019年生态安全等级呈现逐年递增转化且接近向Ⅴ级(安全)的趋势, 由2010年向Ⅰ级(不安全)转化发展至2019年向Ⅴ级(安全)转化, 其中2018年已达到Ⅴ级(安全)。划分为四个阶段:(1)2010-2012年:关联度量值kj(w2012)<kj(w2010)<kj(w2011), 说明2011年生态安全水平较2010年有所提升且属于Ⅰ级(不安全), 2012年生态安全水平较2010-2011年有所下滑; (2)2013-2015年:2013年生态安全关联度最大值-0.0496<0, 有向Ⅱ级(较不安全)转化的趋势, 2014-2015年关联度量值kj(w2015)<kj(w2014)<0, 说明2015年生态安全水平较2014年有所下降且均属于向Ⅲ级(临界安全)转化; (3)2016-2017年:关联度量值kj(w2017)<kj(w2016)<0, 说明2017年生态安全较2016年有所下降且均属于向Ⅳ级(较安全)转化; (4)2018-2019年:2018年最大值0.0484属于Ⅴ级(安全), 但2019年关联度量值-0.0476<0, 整体而言已趋近于0, 则表示具备转化为Ⅴ级(安全)的趋势。
综上, 自2016年二道白河镇的生态安全水平向Ⅳ级(较安全)、Ⅴ级(安全)转化, 整体生态安全等级不断上升且趋安全等级发展, 为了使二道白河镇生态安全水平日后发展过程中处于稳定的Ⅴ级(安全), 需针对制约因素及时合理调整, 制定相应政策措施。
3.2 生态安全各层面影响要素变化分析同理可计算出二道白河镇2010-2019年按DPSIR-EES层面的各年间综合关联度(图 2)。由图 2可知, 除EP2在2010年为Ⅴ级(安全)外, 其余层面关联度量值等级均为Ⅰ级(不安全), 具体结果:(1)EP1、ES1、ES2、EI1关联度量值, 按照Ⅰ级(不安全)-Ⅱ级(较不安全)-Ⅲ级(临界安全)-Ⅳ级(较安全)-的趋势逐年发展; (2)EP2在2010年为Ⅴ级(安全), 2011年降至Ⅰ级(不安全)-2012年至2014年提升至Ⅱ级(较不安全)-2015年、2016年回降至Ⅰ级(不安全)-2017年提升至Ⅲ级(临界安全)-2018、2019年升回Ⅴ级(安全), 呈现明显的波动; (3)SP3、ER1在2015年第一次达到Ⅴ级(安全), 此外SP3在2017、2019年达到Ⅴ级(安全)、2018年降至Ⅰ级(不安全), 而ER1则在2018、2019年再次达到Ⅴ级(安全); (4)除上述层面外, 其余层面各年间变化均呈现上下一至两级波动变化趋势。
![]() |
图 2 2010-2019年二道白河镇生态安全DPSIR-EES评价结果 Fig. 2 DPSIR-EES evaluation results of ecological security in Erdaobaihe town from 2010 to 2019 |
无论从DPSIR层面或EES层面分析, 2010-2019年二道白河镇生态安全水平均大致呈现出从Ⅰ级(不安全)向Ⅴ级(安全)逐年转化的上升趋势。从DPSIR层面看, 驱动力(D)、状态(S)层面生态安全水平在2019年处于Ⅴ级(安全), 其余层面在2019年均处于向Ⅴ级(安全)转化水平; 影响(I)层面的关联度量值明显低于其余四个层面, 除2010-2012年外均<0, 说明一直处于转化水平, 是影响二道白河镇生态安全水平基于DPSIR层面不断提升的制约因素。从EES层面看, 环境(E)在2019年处于Ⅴ级(安全), 其余层面在2018年处于Ⅴ级(安全), 但在2019年处于向Ⅴ级转化水平; 社会(S)层面的关联度量值明显低于其余两个层面, 除2012、2018年外均<0, 说明一直处于转化水平, 是影响二道白河镇生态安全水平基于EES层面不断提升的制约因素。
年份 Year DPSIR各维度 DPSIR dimensions |
驱动力 Driving force | 压力 Pressure | 状态 State | 影响 Impact | 响应 Response | |||||||||
关联度 Correlation degree |
等级 Grade |
关联度 Correlation degree |
等级 Grade |
关联度 Correlation degree |
等级 Grade |
关联度 Correlation degree |
等级 Grade |
关联度 Correlation degree |
等级 Grade |
|||||
2010 | 0.0462 | Ⅰ级 | -0.2299 | 向Ⅰ级转化 | 0.1082 | Ⅰ级 | 0.1339 | Ⅰ级 | -0.0555 | 向Ⅰ级转化 | ||||
2011 | 0.0769 | Ⅰ级 | 0.0769 | Ⅰ级 | 0.1013 | Ⅰ级 | 0.2658 | Ⅰ级 | -0.0348 | 向Ⅰ级转化 | ||||
2012 | 0.0164 | Ⅰ级 | -0.1568 | 向Ⅱ级转化 | 0.0114 | Ⅱ级 | 0.1646 | Ⅰ级 | -0.1389 | 向Ⅱ级转化 | ||||
2013 | 0.0104 | Ⅱ级 | 0.0663 | Ⅱ级 | 0.0507 | Ⅲ级 | -0.1170 | 向Ⅱ级转化 | -0.0886 | 向Ⅲ级转化 | ||||
2014 | 0.0346 | Ⅱ级 | 0.0686 | Ⅲ级 | -0.0237 | 向Ⅲ级转化 | -0.1322 | 向Ⅱ级转化 | -0.2434 | 向Ⅱ级转化 | ||||
2015 | 0.0034 | Ⅲ级 | -0.0918 | 向Ⅳ级转化 | 0.0511 | Ⅲ级 | -0.2206 | 向Ⅳ级转化 | -0.2181 | 向Ⅱ级转化 | ||||
2016 | -0.0576 | 向Ⅳ级转化 | -0.0961 | 向Ⅳ级转化 | 0.0319 | Ⅳ级 | -0.1640 | 向Ⅴ级转化 | -0.2013 | 向Ⅳ级转化 | ||||
2017 | -0.0981 | 向Ⅴ级转化 | -0.1145 | 向Ⅳ级转化 | 0.0665 | Ⅳ级 | -0.0150 | 向Ⅳ级转化 | -0.2358 | 向Ⅳ级转化 | ||||
2018 | 0.0798 | Ⅴ级 | -0.0868 | 向Ⅳ级转化 | 0.3927 | Ⅴ级 | -0.2096 | 向Ⅴ级转化 | 0.0577 | Ⅴ级 | ||||
2019 | 0.0094 | Ⅴ级 | -0.1613 | 向Ⅴ级转化 | <0.0001 | Ⅴ级 | -0.0412 | 向Ⅴ级转化 | -0.0630 | 向Ⅴ级转化 |
年份 Year |
环境 Environment | 经济 Economics | 社会 Sociology | |||||
关联度 Correlation degree |
等级 Grade |
关联度 Correlation degree |
等级 Grade |
关联度 Correlation degree |
等级 Grade |
|||
2010 | 0.0037 | Ⅰ级 | 0.0136 | Ⅰ级 | -0.0184 | 向Ⅰ级转化 | ||
2011 | 0.1752 | Ⅰ级 | 0.0043 | Ⅰ级 | -0.0420 | 向Ⅰ级转化 | ||
2012 | -0.0354 | 向Ⅱ级转化 | -0.0817 | 向Ⅱ级转化 | 0.0073 | Ⅰ级 | ||
2013 | -0.0716 | 向Ⅱ级转化 | -0.0491 | 向Ⅱ级转化 | -0.0284 | 向Ⅱ级转化 | ||
2014 | -0.0672 | 向Ⅲ级转化 | 0.0468 | Ⅱ级 | -0.2001 | 向Ⅱ级转化 | ||
2015 | -0.1099 | 向Ⅳ级转化 | -0.1858 | 向Ⅲ级转化 | -0.1134 | 向Ⅲ级转化 | ||
2016 | -0.0157 | 向Ⅳ级转化 | -0.1951 | 向Ⅳ级转化 | -0.1122 | 向Ⅳ级转化 | ||
2017 | -0.1396 | 向Ⅳ级转化 | -0.1124 | 向Ⅳ级转化 | -0.1772 | 向Ⅴ级转化 | ||
2018 | 0.0188 | Ⅴ级 | 0.0490 | Ⅴ级 | 0.0772 | Ⅴ级 | ||
2019 | 0.0108 | Ⅴ级 | -0.1402 | 向Ⅴ级转化 | -0.0142 | 向Ⅴ级转化 |
(1) 物元分析法将研究问题概括为关联或不关联并进行转化处理, 能更为清晰地揭示当中分异信息。将旅游型城镇生态安全特征及特征量值作为物元, 对评价对象进行量化处理, 结合二道白河镇实际状况及相关指标标准得到经典域、节域, 建立物元评价模型, 借助熵权法和均方差法得出综合权重, 计算出生态安全关联度量值及对应等级, 克服评价过程中主观人为影响, 提高生态安全评价的精度。
(2) 基于DPSIR-EES层面, 结合物元分析模型对2010-2019年二道白河镇生态安全进行评价, 结果表明:2010-2013年二道白河镇生态安全水平较差, 其中GDP增长率、社会保障建设水平、节能环保力度等三项指标影响最为突出, 自2010年处于安全级后呈现逐年下降趋势, 但在2016年后二道白河镇的各项政策与措施发挥显著作用, 生态安全水平逐年向安全水平转化。基于DPSIR-EES评价指标体系进行各层面要素变化分析, 经济压力(EP2)层面的波动变化趋势明显, 影响(I)层面、社会(S)层面已成为二道白河镇生态安全进一步提升的制约因素。
(3) 综合DPSIR-EES-物元分析, 二道白河镇下一阶段的工作重点:第一, 应不断提升经济发展水平, 同时重视GDP增长率的变化情况以及人均可支配收入水平; 第二, 注重人民生活水平质量的提升以及社会保障体系的建设, 重视环境开发与保护, 实现环境、经济、社会和谐发展; 第三, 政府积极采取响应措施弥补由于生产、生活行为的消极影响给生态环境系统带来的压力, 目前应加强节约资源和能源、加大节能环保力度、推进清洁生产和循环经济等, 从根本上降低区域生态环境压力, 才能保证二道白河镇生态安全迈向更高的水平。
4.2 讨论引入物元分析模型对旅游型城镇生态安全进行评价, 基本实现分异及综合指标评价信息的目的。但物元模型要求必须确定评价的经典域及节域, 而目前较多标准和准则中仅有标准等级值, 而并非对应标准等级的量值范围, 因此对于旅游型城镇生态安全评价经典域和节域的界定等问题还有待进一步研究和探讨。同时, 生态系统的动态平衡往往由内力和外力共同确定, 对于旅游型城镇而言, 空间范围尺度明显差别于区域层面(大尺度或中等尺度), 致使其难以构成一个完整的生态系统, 在探究过程中更多的考虑内力影响, 对于外力的影响分析仅做表面探究, 还有待进一步深挖。另外, 由于旅游产业边界模糊, 综合交叉所涉及的内容复杂, 以旅游业为主导产业的旅游型城镇研究过程中很难将因单一原因造成生态环境破坏、环境污染等数据一一剥离, 致使指标体系设定和数据统计工作难度加大, 这将是日后研究与旅游相关的生态安全需面临的挑战。
[1] |
秦晓楠, 卢小丽, 武春友. 国内生态安全研究知识图谱——基于Citespace的计量分析. 生态学报, 2014, 34(13): 3693-3703. |
[2] |
马克明, 傅伯杰, 黎晓亚, 关文彬. 区域生态安全格局: 概念与理论基础. 生态学报, 2004, 24(4): 761-768. DOI:10.3321/j.issn:1000-0933.2004.04.017 |
[3] |
郭付友, 吕晓, 于伟, 任嘉敏, 初楠臣. 山东省绿色发展水平绩效评价与驱动机制——基于17地市面板数据. 地理科学, 2020, 40(2): 200-210. |
[4] |
Ruan W Q, Li Y Q, Zhang S N, Liu C H. Evaluation and drive mechanism of tourism ecological security based on the DPSIR-DEA model. Tourism Management, 2019, 75: 609-625. DOI:10.1016/j.tourman.2019.06.021 |
[5] |
庞雅颂, 王琳. 区域生态安全评价方法综述. 中国人口·资源与环境, 2014, 24(S1): 340-344. |
[6] |
陈昕, 彭建, 刘焱序, 杨旸, 李贵才. 基于"重要性-敏感性-连通性"框架的云浮市生态安全格局构建. 地理研究, 2017, 36(3): 471-484. |
[7] |
彭建, 赵会娟, 刘焱序, 吴健生. 区域生态安全格局构建研究进展与展望. 地理研究, 2017, 36(3): 407-419. |
[8] |
翁钢民, 潘越, 李凌雁. 基于改进DPSIR-DS模型的旅游生态安全等级测度及时空演变分析——以"丝绸之路"沿线五省区为例. 旅游科学, 2018, 32(6): 17-32. |
[9] |
张建清, 张岚, 王嵩, 范斐. 基于DPSIR-DEA模型的区域可持续发展效率测度及分析. 中国人口·资源与环境, 2017, 27(11): 1-9. |
[10] |
王磊, 郭灿, 李慧明. 基于PSR-TOPSIS模型的宁夏回族自治区土地生态安全评价. 水土保持研究, 2016, 23(6): 154-159. |
[11] |
柯小玲, 向梦, 冯敏. 基于灰色聚类法的长江经济带中心城市生态安全评价研究. 长江流域资源与环境, 2017, 26(11): 1734-1742. DOI:10.11870/cjlyzyyhj201711002 |
[12] |
黄烈佳, 杨鹏. 长江经济带土地生态安全时空演化特征及影响因素. 长江流域资源与环境, 2019, 28(8): 1780-1790. |
[13] |
呙亚玲, 李巧云. 基于改进PSR模型的洞庭湖区生态安全评价及主要影响因素分析. 农业现代化研究, 2021, 42(1): 132-141. |
[14] |
Gari S R, Newton A, Icely J D. A review of the application and evolution of the DPSIR framework with an emphasis on coastal social-ecological systems. Ocean & Coastal Management, 2015, 103: 63-77. |
[15] |
Ehara M, Hyakumura K, Sato R, Kurosawa K, Araya K, Sokh H, Kohsaka R. Addressing maladaptive coping strategies of local communities to changes in ecosystem service provisions using the DPSIR framework. Ecological Economics, 2018, 149: 226-238. DOI:10.1016/j.ecolecon.2018.03.008 |
[16] |
崔馨月, 方雷, 王祥荣, 康俊锋. 基于DPSIR模型的长三角城市群生态安全评价研究. 生态学报, 2021, 41(1): 302-319. |
[17] |
王大海, 张荣群, 艾东, 孙玮健. 基于EES-PSR的土地生态安全物元模型评价方法实证研究. 农业机械学报, 2017, 48(S1): 228-237. DOI:10.6041/j.issn.1000-1298.2017.S0.036 |
[18] |
马金涛, 苗红, 温胜强, 卜新春. 基于P-S-R模型的甘肃省旅游生态安全时空格局演变. 湖南科技大学学报: 自然科学版, 2020, 35(3): 118-124. |
[19] |
赵文祯, 韩增林, 闫晓露, 钟敬秋. 基于生态系统服务多情景权衡的生态安全格局构建——以大连市瓦房店为例. 自然资源学报, 2020, 35(3): 546-562. |
[20] |
刘晓阳, 曾坚, 贾梦圆, 张森. 闽三角城市群生态安全格局构建及城镇扩展模拟. 生态学报, 2020, 40(21): 7873-7885. |
[21] |
魏黎灵, 李岚彬, 林月, 吴仪, 祁新华. 基于生态足迹法的闽三角城市群生态安全评价. 生态学报, 2018, 38(12): 4317-4326. |
[22] |
何刚, 王雯雯, 夏业领, 朱艳娜. 基于正态云模型的区域生态安全评价. 安全与环境学报, 2018, 18(2): 807-813. |
[23] |
王毅, 魏江超, 孙启元, 王菲凤. 基于ARIMA-ANN模型的生态安全评价及预测——以河西走廊城市群为例. 生态学杂志, 2020, 39(1): 326-336. |
[24] |
何刚, 金兰, 朱艳娜. 基于云物元的合肥市生态安全评价及生态阻力分析. 地域研究与开发, 2020, 39(6): 171-176. DOI:10.3969/j.issn.1003-2363.2020.06.030 |
[25] |
黄和平, 杨宗之. 基于PSR-熵权模糊物元模型的森林生态安全动态评价——以中部6省为例. 中国农业资源与区划, 2018, 39(11): 42-51. |
[26] |
曹秉帅, 邹长新, 高吉喜, 何萍. 生态安全评价方法及其应用. 生态与农村环境学报, 2019, 35(8): 953-963. |
[27] |
张凤太, 王腊春, 苏维词. 基于物元分析-DPSIR概念模型的重庆土地生态安全评价. 中国环境科学, 2016, 36(10): 3126-3134. DOI:10.3969/j.issn.1000-6923.2016.10.042 |
[28] |
雷金睿, 陈宗铸, 陈毅青, 陈小花, 李苑菱, 吴庭天. 1990-2018年海南岛湿地景观生态安全格局演变. 生态环境学报, 2020, 29(2): 293-302. |
[29] |
苏凯, 王茵然, 孙小婷, 岳德鹏. 基于GIS与RS的东北森林带景观格局演变与模拟预测. 农业机械学报, 2019, 50(12): 195-204. DOI:10.6041/j.issn.1000-1298.2019.12.022 |
[30] |
赵宏波, 马延吉. 基于变权-物元分析模型的老工业基地区域生态安全动态预警研究——以吉林省为例. 生态学报, 2014, 34(16): 4720-4733. |
[31] |
朱婧, 孙新章, 何正. SDGs框架下中国可持续发展评价指标研究. 中国人口·资源与环境, 2018, 28(12): 9-18. |
[32] |
文博, 魏微, 王杰才, 尚英男, 江沂璟, 王琴, 李琳. 环境绩效评估中多种赋权方法的对比分析. 环境保护科学, 2020, 46(1): 41-46. |
[33] |
吕广斌, 廖铁军, 姚秋昇, 邓薇. 基于DPSIR-EES-TOPSIS模型的重庆市土地生态安全评价及其时空分异. 水土保持研究, 2019, 26(6): 249-258, 266-266. |
[34] |
杨建宇, 张欣, 李鹏山, 欧聪, 马瑞明, 朱德海. 基于物元分析的区域土地生态安全评价方法研究. 农业机械学报, 2017, 48(S1): 238-246. DOI:10.6041/j.issn.1000-1298.2017.S0.037 |
[35] |
陈伊多, 杨庆媛, 杨人豪, 曾黎. 基于熵权物元模型的土地生态安全评价——重庆市江津区实证. 干旱区地理, 2018, 41(1): 185-194. |