全 文 ：
摇 摇 摇 摇 摇 生 态 学 报
摇 摇 摇 摇 摇 摇 摇 渊杂匀耘晕郧栽粤陨 载哉耘月粤韵冤
摇 摇 第 猿源卷 第 员圆期摇 摇 圆园员源年 远月摇 渊半月刊冤
目摇 摇 次
中国景观生态学发展历程与未来研究重点 陈利顶袁李秀珍袁傅伯杰袁等 渊猿员圆怨冤噎噎噎噎噎噎噎噎噎噎噎噎噎
城市景观格局演变的水环境效应研究综述 黄摇 硕袁郭青海 渊猿员源圆冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
多功能景观研究进展 汤摇 茜袁丁圣彦 渊猿员缘员冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
空间形态受限型城市紧凑发展研究要要要以厦门岛为例 黄摇 硕袁郭青海袁等 渊猿员缘愿冤噎噎噎噎噎噎噎噎噎噎噎
紫金山森林公园降温效应影响因素 闫伟姣袁孔繁花袁尹海伟袁等 渊猿员远怨冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
城市公园景观空间结构对其热环境效应的影响 冯悦怡袁胡潭高袁张力小 渊猿员苑怨冤噎噎噎噎噎噎噎噎噎噎噎噎
基于 韵宰粤的低丘缓坡建设开发适宜性评价要要要以云南大理白族自治州为例
刘焱序袁彭摇 建袁韩忆楠袁等 渊猿员愿愿冤
噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
生态安全条件下土地利用格局优化要要要以皇甫川流域为例 喻摇 锋袁李晓兵袁王摇 宏 渊猿员怨愿冤噎噎噎噎噎噎噎
新疆玛纳斯河流域 圆园园园要圆园员园年土地利用 辕覆盖变化及影响因素 刘金巍袁靳甜甜袁刘国华袁等 渊猿圆员员冤噎噎噎
基于 郧陨杂和 砸杂的赣江上游流域土地利用动态趋势分析 鲁燕飞袁彭摇 芳袁万摇 韵袁等 渊猿圆圆源冤噎噎噎噎噎噎噎
员怨缘源要圆园员园年三江平原土地利用景观格局动态变化及驱动力 刘吉平袁赵丹丹袁田学智袁等 渊猿圆猿源冤噎噎噎噎
基于斑块评价的三峡库区腹地坡耕地优化调控方法与案例研究 王永艳袁李阳兵袁邵景安袁等 渊猿圆源缘冤噎噎噎噎
贵州省山地鄄坝地系统土地利用与景观格局时空演变 李阳兵袁姚原温袁谢摇 静袁等 渊猿圆缘苑冤噎噎噎噎噎噎噎噎
中国西南地区土地覆盖情景的时空模拟 李摇 婧袁范泽孟袁岳天祥 渊猿圆远远冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
基于移动窗口法的岷江干旱河谷景观格局梯度分析 张玲玲袁赵永华袁殷摇 莎袁等 渊猿圆苑远冤噎噎噎噎噎噎噎噎噎
基于植被覆盖度的藏羚羊栖息地时空变化研究 赵海迪袁刘世梁袁董世魁袁等 渊猿圆愿缘冤噎噎噎噎噎噎噎噎噎噎噎
西南峡谷型喀斯特坡地土壤微生物量 悦尧晕尧孕 空间变异特征 范夫静袁黄国勤袁宋同清袁等 渊猿圆怨猿冤噎噎噎噎噎
峡谷型喀斯特不同生态系统的土壤微生物数量及生物量特征 谭秋锦袁宋同清袁彭晚霞袁等 渊猿猿园圆冤噎噎噎噎噎
长三角地区土地利用时空变化对生态系统服务价值的影响 刘桂林袁张落成袁张摇 倩 渊猿猿员员冤噎噎噎噎噎噎噎
基于视觉廊道的青藏铁路沿线旅游动态景观评价 张瑞英袁席建超袁姚予龙袁等 渊猿猿圆园冤噎噎噎噎噎噎噎噎噎噎
基于 砸杂与 郧陨杂的农村居民点空间变化特征与景观格局影响研究 任摇 平袁洪步庭袁刘摇 寅袁等 渊猿猿猿员冤噎噎噎
生态系统保护现状及保护等级评估要要要以江西省为例 樊乃卿袁张育新袁吕一河袁等 渊猿猿源员冤噎噎噎噎噎噎噎噎
崇明东滩盐沼植被变化对滩涂湿地促淤消浪功能的影响 任璘婧袁李秀珍袁杨世伦袁等 渊猿猿缘园冤噎噎噎噎噎噎噎
基于气候尧地貌尧生态系统的景观分类体系要要要以新疆地区为例 师庆东袁王摇 智袁贺龙梅袁等 渊猿猿缘怨冤噎噎噎噎
黄土丘陵沟壑区景观格局演变特征要要要以陕西省延安市为例 钟莉娜袁赵文武袁吕一河袁等 渊猿猿远愿冤噎噎噎噎噎
不同干扰背景下农业景观异质性要要要以巩义市为例 张晓阳袁梁国付袁丁圣彦 渊猿猿苑愿冤噎噎噎噎噎噎噎噎噎噎
山西高原草地景观的数量分类与排序 张先平袁李志琴袁王孟本袁等 渊猿猿愿远冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
山区夏季地表温度的影响因素要要要以泰山为例 孙常峰袁孔繁花袁尹海伟袁等 渊猿猿怨远冤噎噎噎噎噎噎噎噎噎噎噎
典型岩溶洼地土壤水分的空间分布及影响因素 张继光袁苏以荣袁陈洪松袁等 渊猿源园缘冤噎噎噎噎噎噎噎噎噎噎噎
基于移动窗口法的豫西山地丘陵地区景观异质性分析 李栋科袁丁圣彦袁梁国付袁等 渊猿源员源冤噎噎噎噎噎噎噎噎
桂西北喀斯特区域植被变化趋势及其对气候和地形的响应 童晓伟袁王克林袁岳跃民袁等 渊猿源圆缘冤噎噎噎噎噎噎
喀斯特与非喀斯特区域植被覆盖变化景观分析要要要以广西壮族自治区河池市为例
汪明冲袁王兮之袁梁钊雄袁等 渊猿源猿缘冤
噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
不同干扰背景下景观指数与物种多样性的多尺度效应要要要以巩义市为例
董翠芳袁梁国付袁丁圣彦袁等 渊猿源源源冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
石栎鄄青冈常绿阔叶林土壤有机碳和全氮空间变异特征 杨摇 丹袁项文化袁方摇 晰袁等 渊猿源缘圆冤噎噎噎噎噎噎噎
湘中丘陵区南酸枣阔叶林群落特征及群落更新 易摇 好袁邓湘雯袁项文化袁等 渊猿源远猿冤噎噎噎噎噎噎噎噎噎噎噎
基于 砸月云晕的桂西北喀斯特区植被碳密度空间分布影响因素分析 张明阳袁王克林袁邓振华袁等 渊猿源苑圆冤噎噎噎
期刊基本参数院悦晕 员员鄄圆园猿员 辕 匝鄢员怨愿员鄢皂鄢员远鄢猿缘圆鄢扎澡鄢孕鄢 预 怨园郾 园园鄢员缘员园鄢猿远鄢圆园员源鄄园远
室室室室室室室室室室室室室室
封面图说院 空间发展受限城市的厦门要要要在我国城市化进程中袁中小城市在城镇体系建设中处于中间环节袁起到了联系大城市
和小城镇的作用遥 但是袁每个城市由于发展历史尧社会经济结构尧自然地理形态等因素的不同袁都有其发展的特性袁
这些问题都必须要因地制宜地去把握遥 例如袁厦门岛相对隔离袁没有多余的发展空间袁该城市以居住功能为主袁城市
功能较为单一袁公共服务功能和商业服务功能比例较小遥 研究这样紧凑型的城市发展必须要考虑该城市结构转换
的承受力袁周边社会经济环境以及居民的生活习惯等遥
彩图及图说提供院 陈建伟教授摇 北京林业大学摇 耘鄄皂葬蚤造院 糟蚤贼藻泽援糟澡藻灶躁憎岳 员远猿援糟燥皂
第 34 卷第 12 期
2014年 6月
生 态 学 报
ACTA ECOLOGICA SINICA
Vol.34,No.12
Jun.,2014
http: / / www.ecologica.cn
基金项目:国土资源部公益性行业科研专项经费课题(201311001鄄2)
收稿日期:2013鄄12鄄23; 摇 摇 修订日期:2014鄄04鄄01
*通讯作者 Corresponding author.E鄄mail: jianpeng@ urban.pku.edu.cn
DOI: 10.5846 / stxb201312233004
刘焱序.基于 OWA的低丘缓坡建设开发适宜性评价———以云南大理白族自治州为例.生态学报,2014,34(12):3188鄄3197.
Liu Y X, Peng J, Han Y N, Wei H, Du Y Y.Suitability assessment for building land consolidation on gentle hillside based on OWA operator:a case in Dali
Bai Nationality Borough in Yunnan, China.Acta Ecologica Sinica,2014,34(12):3188鄄3197.
基于 OWA的低丘缓坡建设开发适宜性评价
———以云南大理白族自治州为例
刘焱序,彭摇 建*,韩忆楠,魏摇 海,杜悦悦
(北京大学城市与环境学院,地表过程分析与模拟教育部重点实验室, 北京摇 100871)
摘要:山区城市在城镇化、粮食安全、生态保护的多重要求下,有必要探索将基本农田以外的低丘缓坡山地开发为建设用地的可
行性和限制性。 研究选取涵盖景观风险和灾害风险的共 12 个空间化指标,基于有序加权平均( ordered weighted averaging,
OWA)方法,以大理白族自治州为例进行低丘缓坡区域建设开发适宜性评价。 通过设置不同位序权重表征评价者偏好,结果表
明:在指标风险程度的可信度较高的情况下,农用地、城镇、断裂带、河流等高风险指标权重被加大,空间特征明显;在指标风险
程度的可信度较低的情况下,评价结果的空间均质性增强。 通过模拟城镇建设导向、维持现有政策、风险控制导向 3种适宜性
分区情景,得出在城镇建设导向下应优先开发中小城镇,在风险控制导向下应优先保证大城市合理扩展。 评价结果可以描述在
城市的不同发展阶段如何构建政策权衡后的城镇化空间格局,满足了不同决策思路下的建设用地开发布局需求。
关键词:低丘缓坡; 建设开发适宜性; 有序加权平均; 情景分析
Suitability assessment for building land consolidation on gentle hillside based on
OWA operator:a case in Dali Bai Nationality Borough in Yunnan, China
LIU Yanxu, PENG Jian*, HAN Yinan, WEI Hai, DU Yueyue
Laboratory for Earth Surface Processes, Ministry of Education, College of Urban and Environmental Sciences,Peking University,Beijing 100871,China
Abstract: Under the multiple requirements of urbanization, food security and ecosystem management, it is necessary to
explore the feasibilities and limitations on how to develop building land on gentle hillside in mountain cities. In the process
of mountain development, the ecological sensitive area should be prevented from human interferences and human health from
disaster threats. In this study, 12 spatial indices were selected which can represent the landscape and disaster risks in order
to find the suitable degree of building land development on gentle hillside. Ordered weighted averaging (OWA) operators
are usually used to model a family of parameterized decision strategies, which can be applied to route planning to provide
more realistic results. In this method the optimistic degree can be quantified by the ordered weight. In an optimistic
preference, the function of the relatively important indices is magnified; while in a pessimistic preference, the importance is
relatively ignored. OWA algorithm can be used to conduct the scenario simulations of the feasibility of various construction
land under different preferences,which is apparently able to reduce the impact of a single result of a subjective perception
by decision鄄makers, and thus can well reflect the changes of evaluation result led by the slight adjustment of regional
policies. In this study, Dali Bai Nationality Autonomous Prefecture was chosen as the study area, and OWA operators were
adopted to measure the suitability of building land construction on gentle hillside. By setting up different ordered weights,
scenarios can be set with the quantified optimistic or pessimistic preference of the decision makers. The result showed that in
http: / / www.ecologica.cn
the optimistic scenario, the characteristics of high risk indices like farm land, urban, fault zone, rivers can be obviously
detected spatially; in the pessimistic scenario where the degree of reliability for risk indices was decreased, the spatial
homogeneity was enhanced on the result layer. Then, three scenarios were built up by the OWA operators, which support
three different strategy orientations—urban construction, keeping the present policies and risk control. After defining the
suitability sub region, it was found that priority development of small towns should be encouraged in urban construction
scenario, while ensuring the land supply for big cities in the risk control scenario. With the assessments, we can describe
how to build up the spatial pattern of urbanization after a policy trade鄄off among different stages of city development, and the
results can satisfy the different decisions from the government and various requirements for the layout of construction land
from different stakeholders. This study tended to draw a comprehensive conclusion through a combination of inter鄄scene.
Although in a mathematical level, the three policy鄄perspectives provided by this paper are described as the preference of
evaluator忆s, they are actually a reflection of decision makers忆 different strategic tendencies and thus are no clear good or
bad. This study didn忆 t attempt to select an optimal solution, but tends to describe how to build a spatial pattern of
urbanization with policy trade鄄offs in different stages of urban development, so as to meet the needs of construction land忆s
development under different decisions or layout ideas, with relatively high practical significance.
Key Words: gentle hillside; suitability of building land consolidation; ordered weighted averaging; scenario analysis
摇 摇 土地适宜性评价是评定土地对于某种用途是否
适宜以及适宜的程度,中国传统的土地评价研究多
以农林牧业综合适宜性评价为主[1]。 近年来,土地
适宜性评价的领域正在不断扩展,针对特色功能区
域进行的单一用途土地评价备受重视[2]。 在土地类
型层面上,对耕地质量[3]、都市开发[4]、城区重建[5]、
植被恢复[6]、绿色基础设施[7]等单一目标地类适宜
等级评价层出不穷,推动了传统评价理论框架的不
断完善。 土地生态适宜性评价以生物多样性和生境
保护为导向设置评价指标[8],可以有效指导生态脆
弱区域的土地开发相关规划。 尤其对于丘陵地区或
山区而言,城市往往分布在地势相对较低的河谷地
带,在建设用地空间布局上受到较多的地理本底条
件限制。 由于山区本身耕地面积较小,耕地保护压
力大,而为数不多的优质耕地又距离城市近,容易被
开发为建设用地。 为了保障优质耕地,维护粮食安
全,国家和地方对此类耕地的开发有严格的指标约
束,大幅提升了在平整河谷布局建设用地的难度。
因此,为了优化土地资源配置,缓解供需矛盾,有必
要探索山区城市将基本农田以外的低丘缓坡山地开
发为建设用地的生态可行性和限制性,确定区域建
设用地开发的适宜范围。
随着 GIS技术与多准则评价方法相结合在土地
适宜性评价中被广泛应用,大量研究采用层次分析
法(AHP)对栅格图层进行赋权并叠加的方法得出评
价结果[9]。 由于评价选取的地理或生态因子往往是
存在相互关联的,简单的对加权后的图层进行叠加
往往会造成一些有明显限制作用的因子被另一些指
标所补偿,作为一个折中后的结果往往失去了评价
者所要表述的限制性。 因此,美国学者 Yager提出了
有序加权平均(OWA)的概念[10],对因子属性的重要
程度进行排序,通过排序权重的变换模拟不同程度
的评价者偏好下得出的决策结果。 目前,OWA 方法
已逐步应用到国内对地理鄄生态过程研究中,在循环
经济综合评价[11]、购房地理空间决策[12]等研究视角
上取得了一些进展,但该方法目前在土地适宜性评
价中的应用并不多见。 基于 OWA 算法可以对多种
不同偏好下的建设用地开发可行性作以情景模拟,
显然能够降低某一种单一结果受决策者的主观认知
的影响,良好地反映区域政策微调导致的评价结果
变化。 基于此,本研究拟应用 OWA算法对低丘缓坡
山区的建设用地开发适宜性作以情景模拟,探索
OWA算法在土地适宜性评价中的适用性,并为山区
城市的城镇化布局提供决策支持。
1摇 研究区概况
大理白族自治州地处云南省西部,地处云贵高
原与横断山脉结合部位,全州山区面积占全州国土
面积的 90%以上,是典型的西南山区城市,如图 1 所
示。 大理州当地居民将面积较大的山间小盆地、河
9813摇 12期 摇 摇 摇 刘焱序摇 等:基于 OWA的低丘缓坡建设开发适宜性评价———以云南大理白族自治州为例 摇
http: / / www.ecologica.cn
谷沿岸和山麓地带称为“坝子冶,坝区内土壤肥沃、灌
溉便利,是优质耕地的集中区域。 然而近年来,随着
地区工业与旅游业的快速发展,大理州在城镇化过
程中有大量的建设用地需求。 由于坝区人口稠密交
通便利,往往是建设用地开发的首选区域,致使研究
区内为数不多的优质良田迅速减少。 “积极培育中
心城市,择优发展小城镇冶是大理州城镇体系规划的
明确要求。 然而,中心城市大理市毗邻苍山洱海,生
态保护压力大,致使城市建设用地的布局受到明显
制约。 而另一方面中心城市的极化作用吸引了大量
人口的流入和企业的进驻,对建设用地的刚性需求
不断上升。 因此,有必要通过开发市区周边的低丘
缓坡山地满足城市化的需要。 同时,周边各县的中
心镇一般位于坝区中央,在发展过程中容易侵占坝
区的优质基本农田,对全州未来的粮食安全构成威
胁。 而非坝区部分一般生境质量较好且有一定坡
度,如何控制生态风险和灾害风险是建设开发选址
的核心环节。 因此,有必要在空间上对研究区低丘
缓坡建设开发的适宜程度作以定量化,判别出建设
用地开发的优先区域。
图 1摇 研究区位置
Fig.1摇 location of the study area
在对大理州低丘缓坡区域的界定过程中,必须
考虑到我国东西部的地理条件差异较大,对于西部
地区低丘缓坡范围的界定以及低丘缓坡建设用地开
发的限制指标均不能直接照搬东部沿海地区的经
验[13]。 由于滇西地区地势相对较高,地表起伏剧
烈,过多的在地形角度进行限制可能会导致建设用
地布局的极度破碎化,不符合现实需要。 因此,本研
究根据大理州的海拔高度与地形起伏,暂时不将绝
对高程作为限制性因素,只将坡度 25毅作为低丘缓坡
范围的上限。 同时,山区中城镇周边河谷部分相对
平坦,有大面积基本农田,不满足研究目标,因此,本
研究将坝区范围作为低丘缓坡界定的下限。
2摇 研究方法
2.1摇 数据预处理
在数据准备过程中,按照城镇化、粮食安全、生
态保护和灾害防治的多重导向,选择土地利用、地质
灾害、地下水、地震、归一化植被指数(NDVI)、地形、
降水等多个数据集。 其中农地、城镇、公路、河流数
据解译自 2010年春季 ETM+影像, 并转换为矢量格
式对照 Google Earth进行人工修正;林地类型源自云
南大学生态学与地植物学研究所编制的云南省 1颐25
万植被图和林业科学数据中心相关森林资源数据;
地下水保护区位置由大理州地下水功能区划的二级
功能区图矢量化得到;地质灾害类型和断裂带位置
矢量化自大理州地质灾害防治相关图件;地震位置
取自历史震害资料中 1500 年以来 5 级以上地震;
NDVI数据源自地理空间数据云平台的 MOD13Q1数
据集,预处理过程中对 2010 全年的 NDVI 产品求年
平均;DEM数据源自地理空间数据云平台的 GDEM
数据集,在 ArcGIS 中生成坡度图层;降水数据源自
中国气象科学数据共享网的多年平均站点数据集,
按薄板样条函数插值方法以经纬度作为自变量、高
程作为协变量在 ANUSPLIN 软件中插值所得[14]。
所有数据在归一化处理后均统一生成 30 m 分辨率
栅格图层。
2.2摇 适宜性指标构建
在山地开发的过程中,既要保证生态敏感区域
不受干扰,又要防止人类健康受到灾害威胁。 本研
究将这两方面指标集归纳为景观风险和灾害风险,
风险越大则越不适宜建设用地开发。 按照景观风险
和灾害风险的分类以及指标的可获得性,本研究共
选取 12个指标参与评价。 暂定景观风险和灾害风
险作为两个准则层对山地建设开发的贡献为同等重
要,分别赋权为 0.5。 每个准则层包含 6 个指标,指
标相对于准则层的权重由 AHP 法确定。
在景观风险的选取和赋值中主要考虑三个层
面。 首先,采用农地类型与林地类型表征土地利用
方式影响的建设开发限制性。 农地类型按农地景观
0913 摇 生摇 态摇 学摇 报摇 摇 摇 34卷摇
http: / / www.ecologica.cn
被开发所造成的粮食安全风险赋值;林地类型主要
按照林地保护区等级赋值。 其次,采用公路通达性
和城镇距离表征新开发建设用地与已有城乡聚落的
空间联系。 显然距离公路景观越远则越不适宜布局
建设用地,建设用地在公路两侧开发的通达性最高;
而距离大型聚落越近则建设用地开发的需求也会相
应增大。 这里大型聚落最终选取大于 5 hm2的城镇
与农村居民点图斑作为阈值。 最后,选择植被绿度
和地下水保护作为建设开发的辅助性指标。 通过
NDVI反映植被生物量,可以作为建设开发成本的侧
面衡量;建设用地开发增大了地下水的供给压力,因
此地下水保护区也应作为辅助性限制条件。
在灾害风险的选取和赋值中也主要考虑三个方
面。 首先,在灾害易发性角度考虑地质灾害密度和
地震密度。 在地质灾害易发性的表征中,以核密度
函数[15]表征灾害点空间集聚,给不同类型地质灾害
赋予不同权重,本研究按照灾害的危害差异赋权为
泥石流 100,滑坡 80,崩塌 70,塌陷与地裂缝 60,不稳
定斜坡 50,以默认的 6871 m 为搜索半径,生成栅格
图层;地震易发性的表征也采用核密度方法,其中由
于地震点的两两距离相对较远,参照默认半径取 20
km为搜索范围。 其次,在灾害的危险性角度考虑坡
度和降雨侵蚀力的作用,显然坡度越高、降雨侵蚀力
越大,地质灾害的发生强度可能会随之增大。 在处
理过程中虽然坡度 25毅以上不应属于低丘缓坡范围,
但默认为零星的散点可以开发,因此赋值暂不考虑
范围限制;降雨侵蚀力的计算是基于多年月平均数
据,按照常用的 Wischmeier 经验公式得出年降雨侵
蚀力[16]。 最后,立足于规划视角考虑了断裂带距离
和河流距离两种限制性因素,这两种线状指标并不
存在空间衰减关系,而是作为限制条件。 本研究仅
考虑研究区内较明显的几条大断裂带,暂定断裂带
左右的缓冲区内不宜大规模建设;在河流引发的灾
害角度主要考虑洪水的影响,考虑到用水需求,仅设
置了相对较小的缓冲半径作为限制。
每个指标的归一化方法如表 1 所示,其中权重
栏显示了乘以 0.5 后每个指标对于总目标的权重。
两种风险的 AHP 判断矩阵一致性比例 CR 分别为
0郾 058和 0.087,均小于 0.1,通过一致性检验。 权重
较高的指标包括农地类型、城镇距离、断裂带距离、
河流距离等,这是由于:在景观风险指标集的赋权过
程中,本研究认为农用地的保护和与城镇的空间联
系是最为重要的两个因子,前者在供给上决定了低
丘缓坡开发的限制,后者在需求上决定了低丘缓坡
开发的目标;在灾害风险指标集的赋权过程中,本研
究认为断裂带距离和河流距离是以二值化形式出现
的,本身就是明确的限制性条件,因而相对其他因子
更为重要。 权重较低的指标包括植被绿度、地震密
表 1摇 指标归一化方法
Table 1摇 The method of indices normalization
类型 Type 指标 Index 归一化方法 The method of normalization 权重 Weight
景观风险 农地类型 赋值水田为 1,水浇地 0.9,旱地 0.75,园地 0.5,茶园 0.4,属于基本农田则直接赋值为 1 0.164
Landscape risk 林地类型 特殊保护为 1,重点保护 0.75,一般保护 0.5 0.073
公路通达性 10 km以外直接赋值 1,0—10 km按欧氏距离归一化 0.083
城镇距离 赋值方法同上 0.120
植被绿度 按结果分布赋 NDVI均值在 0.45以下为 0,在 0.8以上为 1,0.45—0.8区间内归一化 0.029
地下水保护 不宜开发区和脆弱区为 1,集中供水区 0.9,分散供水区 0.8,应急供水区 0.7,储备区 0.5 0.032
灾害风险 地质灾害密度 观察结果分布将核密度值大于 8赋值为 1,0—8区间内归一化 0.058
Disaster risk 地震密度 按结果分布将核密度值大于 0.8赋值为 1后,0—0.8区间内归一化 0.018
坡度 坡度 25毅以上赋为 1,8毅以下为 0,然后 8—25毅间归一化 0.080
降雨侵蚀力 按结果分布在 3000—4500 MJ·mm·hm
-2·h-1·a-1区间进行正向归一化,区
间以上为 1,区间以下为 0
0.068
断裂带距离 距离断裂带 1 km以内为 1,其他全赋 0 0.130
河流距离 距离河流 500 m以内为 1,其他全赋 0 0.146
1913摇 12期 摇 摇 摇 刘焱序摇 等:基于 OWA的低丘缓坡建设开发适宜性评价———以云南大理白族自治州为例 摇
http: / / www.ecologica.cn
度等,这是由于:在景观风险指标集的赋权过程中
NDVI只是植被长势的一个指标,但植物长势并不能
完全表明生物多样性的重要性,仅是辅助性指标,相
比而言林区的保护等级更为重要;在灾害风险指标
集的赋权过程中地震发生的不确定性因素强,其危
害更多取决于房屋构造和坡度,在市域范围内地震
密度不一定是城市建设的重要限制指标。
2.3摇 有序加权平均法
OWA算法的核心在于对指标按照属性重要性
重新排序,对不同的排序位次赋予不同的次序权重。
传统的图层叠加法实际上是 OWA 算法中默认次序
权重相等的特殊情况。 在 OWA评价中,在评价者较
为乐观的情况下,会给重要性较高的次序图层赋予
较高的权重,评价结果能直接显示指标中最重要的
属性;在评价者较为悲观的情况下,会给重要性较低
的次序图层赋予较高的权重,评价结果不能反映指
标的重要属性。
关于 OWA 权重的确定有 10 余种算法,其中单
调规则递增(RIM)的定量方式出现较早且易于理
解[17]。 按照 Yager的定义,位序权重可表述为[17鄄19]。
w j = QRIM
jæ
è
ç
ö
ø
÷
n
- QRIM
j - 1æ
è
ç
ö
ø
÷
n
,
j = 1,2,...,n , QRIM ( )r = ra (1)
式中,j为位序,w j为位序权重,n为指标数量。 r 为自
变量,a为表征乐观程度的幂指数。
可知,在 a = 1 的情况下位序权重相等,计算转
化为普通的准则权重叠加。 若 a<1,则最重要的属
性位序权重越大,评价者对指标的属性持乐观态度,
只需要重要的前几层算子即可判定;若 a>1,则越重
要的属性位序权重越小,评价者对指标的属性持悲
观态度,后面几层不重要的算子更被重视。 在本研
究中,“乐观冶表明评价者认为生态风险的限制条件
能够制约建设用地的开发,而“悲观冶表明评价者不
认为生态风险指标可以限制建设用地的开发。 根据
上述公式可以将次序权重的取值如表 2所示。
表 2摇 次序权重计算结果
Table 2摇 Result of the ordered weight
位序权重
Ordered
weight
a= 0.0001
最乐观
Most optimistic
a= 0.1
乐观
Better optimistic
a= 0.5
较乐观
Optimistic
a= 1
无偏好
No preference
a= 2
较悲观
Pessimistic
a= 10
悲观
More essimistic
a= 1000
最悲观
Most pessimistic
W1 1.000 0.781 0.286 0.0833 0.006 0.000 0.000
W2 0.000 0.056 0.120 0.0833 0.021 0.000 0.000
W3 0.000 0.035 0.092 0.0833 0.035 0.000 0.000
W4 0.000 0.025 0.077 0.0833 0.049 0.000 0.000
W5 0.000 0.020 0.068 0.0833 0.063 0.000 0.000
W6 0.000 0.017 0.062 0.0833 0.076 0.001 0.000
W7 0.000 0.014 0.057 0.0833 0.090 0.004 0.000
W8 0.000 0.013 0.053 0.0833 0.104 0.013 0.000
W9 0.000 0.011 0.050 0.0833 0.118 0.039 0.000
W10 0.000 0.010 0.047 0.0833 0.132 0.105 0.000
W11 0.000 0.009 0.045 0.0833 0.146 0.257 0.000
W12 0.000 0.009 0.043 0.0833 0.160 0.581 1.000
3摇 结果分析
3.1摇 不同偏好适宜性对比
图 2为赋值并归一化的 12 种适宜性指标,分别
将其带入表 2 所设计的 7 种偏好情况,结果如图 3
所示。 在评价者无偏好的状态下,计算过程中不区
分重要图层,评价结果直接依赖于数据预处理中的
准则权重。 结果表征了海拔较低、坝区面积较大的
大理市区、宾川县和巍山县在坝区周边低丘缓坡地
带的建设用地开发适宜性较低,而同样海拔相对不
高、距离大理市区较远的弥渡县和永平县坝区周边
则显示出了相对较高的开发适宜性。 在评价者最为
乐观的状态下,除了半径 10 km 之内的城镇周边有
微弱的开发适宜性,其他大面积区域呈现出适宜性
为 0 的基质状态,除了城镇指标外其他指标难以体
现。这证明评价者在对风险极端重视的情况下,城
2913 摇 生摇 态摇 学摇 报摇 摇 摇 34卷摇
http: / / www.ecologica.cn
图 2摇 适宜性指标空间分布
Fig.2摇 Spatial pattern of the suitability indices
图 3摇 不同偏好的适宜性评价
Fig.3摇 suitable assessments for different preferences
3913摇 12期 摇 摇 摇 刘焱序摇 等:基于 OWA的低丘缓坡建设开发适宜性评价———以云南大理白族自治州为例 摇
http: / / www.ecologica.cn
郊以外的广大区域是不适宜进行开发的,仅有城市
周边鉴于居民需求的压力不得不进行微量的开发。
在评价者乐观的状态下,风险被高度重视,河流、断
裂带、坡度与降雨侵蚀力高值部分依然适宜性极低,
而其他大多数栅格显示出具有较低适宜性的格局。
在评价者较乐观的状态下,适宜性低值部分和较低
部分出现分化,作为限制性条件的河流和断裂带可
以良好地在空间上被分离出来,成为评价者认知的
生态风险重要来源。 在评价者较悲观的状态下,风
险评价的指标体系可信度低,大多数像元出现适宜
性大于 0.9的结果。 至于后两种评价者悲观和最悲
观的状态下,数据预处理中的风险赋值已不可信,坝
区外全部是适宜开发区域,评价失去意义。
3.2摇 基于多情景的适宜性分区
在选取 a = 0.5 和 a = 2 作为偏好适宜性的次序
权重时,结果显示其数据分离程度是相对较差的,适
宜度比较集中。 这种结果并不符合根据适宜性差异
进行分区管制的决策目标,同时站在决策的角度上,
评价人员的偏好也不可能出现如此大的波动。 因
此,基于情景模拟可用性的导向,本研究暂时默认为
a= 1是在粮食安全、生态安全、社会发展、居民健康
等视角权衡下的正常赋权结果,是现有城镇化、耕地
保护、生态保护等政策要求的权衡值,命名为维持现
有政策视角。 则 a= 1.2可以作为一种较为忽视风险
而重视城镇扩张的政策倾向,命名为城镇建设导向
视角;a= 0.8 可以作为一种较为重视风险而严格限
制城镇扩张的政策倾向,命名为风险控制导向视角。
为排除栅格值向某一方向偏移造成的分区比例严重
失调,采用自然断点法自动将栅格分为 4 类,保证了
a较大的情况下断点阈值也同步增大,将 4类分区按
适宜程度从高到低依次命名为优先开发区、适度开
发区、限制开发区和不开发区。 同时,对坡度图进行
重分类,25毅以上为 1,其他为 0,转为矢量后将面积 1
km2以下的图斑融合到邻接的面积最大图斑中,再提
取出属性值为 1 的部分作为低丘缓坡的上限,从而
排除了一些坡度大于 25毅但不连续的破碎图斑干扰。
将提取出的 25毅以上“陡坡冶范围与坝区范围融合,
得出了非低丘缓坡区域作为不开发区。 最终情景假
设结果如图 4所示。
图 4摇 不同情景的建设开发适宜性分区
Fig.4摇 Suitability of building land consolidation zoning for different scenarios
摇 摇 在城镇建设导向下,永平县虽然目前经济实力
相对较低,但由于地势相对平缓、基本农田总量不
大、坝区面积较小,因而有较大的扩展空间,优先开
发区面积最大;祥云县经济实力相对较强,其同样有
较为平整的地形条件,由于路网密布、地质灾害密度
较低,因而也有较为广阔的优先开发区范围;同时,
并不在城镇附近的巍山县西部由于坡度较缓、耕地
保护压力小,也形成了一定面积的优先开发区;反
之,研究区的人口与经济中心大理市由于邻接苍山
洱海,只有南部和东部有大面积低丘缓坡空间,致使
城镇建设受到制约,相比其他各县优先开发区面积
反而最小。
维持现有政策的模式则体现出略有不同的分区
格局。 地势较低的永平县与巍山县西部的优先开发
区面积缩小一半以上,而北部地势较高的剑川县和
洱源县优先开发区面积减小幅度甚微,这体现出了
远离城镇的风险正在加强;相比之下,同样地势较低
但城镇化水平相对高的祥云县优先开发区范围并未
明显缩小,证实了城镇距离指标的重要作用;另外,
相比城镇建设导向,不宜开发区的面积并未有太多
4913 摇 生摇 态摇 学摇 报摇 摇 摇 34卷摇
http: / / www.ecologica.cn
变化,分区类型的转换主要发生在从优先开发区向
适度开发区过渡。
风险控制导向则体现出了与城镇建设导向较大
的不一致性。 巍山县西部、云龙县北部较大规模的
优先开发区全部消失;永平县的优先开发区缩减至
中心镇周围 10 km 半径范围内,祥云县成为优先开
发区面积最大的地区;城镇建设导向下优先开发区
面积最小的大理市其优先开发区范围变化反而也较
小,且并没有像云龙县、永平县、南涧县一样形成较
大面积的限制开发区。 显然,在风险控制控制视角
下不仅远离城镇的风险继续加强,农地林地保护、坡
度与降雨侵蚀力等因素也在空间上得到体现。
对比 3种情景视角可以发现,由于客观地理条
件的限制,在倾向于城镇建设的导向下,研究区大城
市依然不适于大面积的进行建设用地开发,反而是
相对开发风险较小的中小城镇可以作为产业园区的
优先选址地点。 而在风险控制导向下则不然,自然
环境相对较好的中小城镇出于耕地和生境保护的目
的不能进行大规模扩张;而大城市周边由于存在刚
性需求,即便有生态保护的压力依然会预留一部分
建设用地开发的空间。 目前,在山地城市生态保护
与城镇化的权衡中,这种结论往往可以反映在土地
规划中。 例如建设用地指标多则会把多出的那一部
分规划到较偏远的乡镇或大型村落中,若建设用地
指标少则优先保证城区的发展,侧面证明了本研究
适宜性分区的情景合理性。
4摇 讨论
在多指标综合评价中,指标赋值和权重判定是
结果不确定性的最主要来源。 首先需要明确,评价
本身就是一个基于研究者认知的主观过程,即使是
完全客观化赋权的熵值法也是评价者为了刻意拉开
属性值域差异而做出的一种主观选择,不一定是描
述指标权值的最佳方式。 在无法确认我们所设置的
指标属性值和权值是否足够贴近真实情况的时候,
采用 OWA算法的偏好设置可以有效规避一部分误
差。 例如,采用乐观偏好可以将已有的属性较大的
指标进一步扩大化,采用悲观偏好可以将属性较大
的指标平均化。 因此,可以通过偏好的设置对已设
置的权值做调整,使评价的逻辑性更加严密。
然而,OWA 方法的应用并不意味着准则权重不
再重要。 在本研究中,准则权重较大的几个指标图
层直接决定了乐观偏好下的适宜性评价结果。 按照
对低丘缓坡建设用地开发背景的理解,本研究认为
政策作用具有较大的影响力,那么有必要将城镇化
的需求和耕地保护的压力赋予较高的权重。 此外,
对于一些明显的限制性强但并未通过法律进行约束
的要素例如断裂带、河流,理应赋予较大的权重,否
则其限制作用有可能被其他指标所抵消。 由此可
见,确定准则权重是 OWA评价的关键环节,OWA 方
法的分析结果依然高度依赖于 AHP 或者其他方法
对准则权重设置的合理性。 OWA 方法只能规避部
分误差而不能纠正错误的准则权重设置,这也使其
应用会受到一定程度的制约。
情景分析作为地理-生态过程的常用手段,有两
种不同的结果归纳方式。 一种是在所假设的几种情
景里选取最优值,作为最终结果[20鄄21];而另一种是通
过几种情景的相互印证,归纳出一个综合全部利弊
的结论[22鄄23]。 在土地适宜性评价中,纯粹的生态建
设、耕地保护、城市扩展对于地区而言都是不符合统
筹发展需要的,现实中案例极少。 作为决策支持不
仅要考虑极端情景,还应当重视不同情景间的组合
与阈值变化的过渡,通过互竞权衡发挥土地利用情
景分析的现实意义[24]。 在维持城市可持续发展的
导向下,不仅要关注城镇化的生态与环境效应[25],
还要认识到生态空间格局的动态演变[26]。 随着时
间演变,进行政策微调在现实中是易于发生的,因此
本研究倾向于通过情景间的组合得出一个综合性的
结论。 显然,应用 OWA方法进行土地生态适宜性评
价恰恰满足情景组合分析的理论要求。 将评价者偏
好作为政策波动因素带入权重体系中,一方面避免
了传统情景分析中权重不确定性导致的误差;另一
方面避免了情景的极化,对地区未来发展模式有良
好的预见作用。 因此,使用 OWA方法评价土地生态
适宜性不仅是可行的,而且是可以有效跟实际政策
对接的。 在 OWA的应用体系下,本研究并不试图选
出一种最优方案,而是倾向于描述在城市的不同发
展阶段如何构建政策权衡后的城镇化空间格局,以
满足不同决策思路下的建设用地开发布局需求,具
有相对较高的现实意义。
5摇 结论
本研究基于所选择的 12 个低丘缓坡建设用地
5913摇 12期 摇 摇 摇 刘焱序摇 等:基于 OWA的低丘缓坡建设开发适宜性评价———以云南大理白族自治州为例 摇
http: / / www.ecologica.cn
开发适宜性指标,按照 OWA算法对图层按照属性重
新排序,对不同的排序位次赋予不同的次序权重,表
征了从乐观到悲观 7种评价者偏好下的建设用地开
发适宜性评价。 在乐观情况下,农用地、城镇、断裂
带、河流等高风险指标权重被加大,空间特征明显;
在悲观情况下,指标风险程度的可信度下降,评价结
果的空间均质性增强。
按照 OWA算法模拟城镇建设导向、维持现有政
策、风险控制导向 3 种适宜性分区情景。 在倾向于
城镇建设的导向下,相对开发风险较小的中小城镇
可以作为产业园区的优先选址地点;而在风险控制
导向下,自然环境相对较好的中小城镇出于耕地和
生境保护的目的不能进行大规模扩张,而大城市周
边由于存在刚性需求依然会有一部分建设用地开发
的空间。 通过描述在城市的不同发展阶段如何构建
政策权衡后的城镇化空间格局,可以满足不同决策
思路下的建设用地开发布局需求。
致谢:感谢北京大学城市与环境学院李双成教授对
写作的帮助。
References:
[ 1 ]摇 Ni S X. Progress in the research on land evaluation in China
during the latest ten years. Journal of Natural Resources, 2003,
18(6): 672鄄683.
[ 2 ] 摇 Shi T G, Zheng G Q, Wang Z Y, Wang L L. Progress in
Research on Land Suitability Evaluation in China. Progress in
Geography, 2007, 26(2): 106鄄115.
[ 3 ] 摇 Liu Y L, Jiao L M, Liu Y F, He J H. A self鄄adapting fuzzy
inference system for the evaluation of agricultural land.
Environmental Modelling & Software, 2013, 40: 226鄄234.
[ 4 ] 摇 Gong J Z, Liu Y S, Chen W L. Land suitability evaluation for
development using a matter鄄element model: A case study in
Zengcheng, Guangzhou, China. Land Use Policy, 2012, 29(2):
464鄄472.
[ 5 ] 摇 Wang H, Shen Q P, Tang B S, Skitmore M. An integrated
approach to supporting land鄄use decisions in site redevelopment for
urban renewal in Hong Kong. Habitat International, 2013, 38:
70鄄80.
[ 6 ] 摇 McVicar T R, Van Niel T G, Li L T, Wen Z M, Yang Q K, Li
R, Jiao F. Parsimoniously modelling perennial vegetation
suitability and identifying priority areas to support China忆 s re鄄
vegetation program in the Loess Plateau: Matching model
complexity to data availability. Forest Ecology and Management,
2010, 259(7): 1277鄄1290.
[ 7 ] 摇 Greca P L, Rosa D L, Martinico F, Privitera R. Agricultural and
green infrastructures: The role of non鄄urbanised areas for eco鄄
sustainable planning in a metropolitan region. Environmental
Pollution, 2011, 159(8 / 9): 2193鄄2202.
[ 8 ] 摇 Ferretti V, Pomarico S. Ecological land suitability analysis through
spatial indicators: An application of the Analytic Network Process
technique and Ordered Weighted Average approach. Ecological
Indicators, 2013, 34: 507鄄519.
[ 9 ] 摇 Jeong J S, Garc侏a鄄Moruno L, Hern佗ndez鄄Blanco J. A site planning
approach for rural buildings into a landscape using a spatial multi鄄
criteria decision analysis methodology. Land Use Policy, 2013,
32: 108鄄118.
[10] 摇 Yager R R. On ordered weighted averaging aggregation operators in
multicriteria decisionmaking. IEEE Transactions on Systems,
Man, and Cybernetics, 1988, 18(1): 183鄄190.
[11] 摇 Fan Y, Liu G P. Comprehensive evaluation of circular economy
development in Gansu Province. Areal Research and
Development, 2012, 31(1): 20鄄23.
[12] 摇 Hu Z W, Liu X X, Peng C, Wei L. Spatial multi鄄criteria decision
making for house location selection by ordered weighted averaging
approach. Geographical Research, 32(3): 476鄄486.
[13] 摇 Liu W D, Yan W. Mild slope land exploitation of low mountains
and hills in developed regions: with Yongkang City as an example.
Scientific and Technological Management of Land and Resources,
2007, 24(3): 1鄄5.
[14] 摇 Hutchinson M F. Anusplin Version 4. 2 User Guide. Cauberra:
Australian National University, 2001.
[15] 摇 Shi X. Selection of bandwidth type and adjustment side in kernel
density estimation over inhomogeneous backgrounds. International
Journal of Geographical Information Science, 2010, 24 ( 5 ):
643鄄660.
[16] 摇 Wischmeier W H. A soil erodibility nomograph for farmland and
construction sites. Journal of Soil and Water Conservation, 1971,
26(5): 189鄄193.
[17] 摇 Zarghami M, Szidarovszky F. Fuzzy quantifiers in sensitivity
analysis of OWA operator. Computers and Industrial Engineering.
2008, 54(4): 1006鄄1018.
[18] 摇 Yager R R. Quantifier guided aggregation using OWA operators.
International Journal of Intelligent Systems, 1996, 11(1): 49鄄73.
[19] 摇 Yager R R. Nonmonotonic OWA operators. Soft Computing, 1999,
3(3): 187鄄196.
[20] 摇 Tong K K, Ma K M. Significant impact of job鄄housing distance on
carbon emissions from transport: a scenario analysis. Acta
Ecologica Sinica, 2012, 32(10): 2975鄄2984.
[21] 摇 Bai Y, Zheng H, Zhuang C W, Ouyang Z Y. Ecosystem services
valuation and its regulation in Baiyangdian baisn: Based on
InVEST model. Acta Ecologica Sinica, 2013, 33(3): 711鄄717.
[22] 摇 Yin H W, Kong F H, Qi Y, Wang H Y, Zhou Y N, Qin Z M.
Developing and optimizing ecological networks in urban
agglomeration of Hunan Province, China. Acta Ecologica Sinica,
6913 摇 生摇 态摇 学摇 报摇 摇 摇 34卷摇
http: / / www.ecologica.cn
2011, 31(10): 2863鄄2874.
[23] 摇 Xiao X, Zhang J, Lu J Y, Zhong S E, Yin L J. Analysis on
spatial structure and scenarios of carbon dioxide emissions from
tourism transportation. Acta Ecologica Sinica, 2012, 32 ( 23):
7540鄄7548.
[24] 摇 Li P, Jiang L G, Feng Z M, Yu X B. Research progress on trade鄄
offs and synergies of ecosystem services: an overview. Acta
Ecologica Sinica, 2012, 32(16): 5219鄄5229.
[25] 摇 Xie M M, Wang Y L, Fu M C. An overview and perspective about
causative factors of surface urban heat island effects. Progress in
Geography, 2011, 30(1): 35鄄41.
[26] 摇 Gong J, Xie Y C, Sun P, Jie Y W. Spatio鄄temporal change on
township scale of Jinta oasis in the latest 20 years. Acta Ecologica
Sinica, 2013, 33(11): 3470鄄3479.
参考文献:
[ 1 ]摇 倪绍祥. 近 10年来中国土地评价研究的进展. 自然资源学报,
2003, 18(6): 672鄄683.
[ 2 ] 摇 史同广, 郑国强, 王智勇, 王林林. 中国土地适宜性评价研究
进展. 地理科学进展, 2007, 26(2): 106鄄115.
[11] 摇 樊元, 刘国平. 甘肃省循环经济发展水平综合评价. 地域研究
与开发, 2012, 31(1): 20鄄23.
[12] 摇 胡卓玮, 刘晓旭, 彭程, 魏铼. 基于次序权重平均法的购房选
择地理空间多准则决策. 地理研究, 2013, 32(3): 476鄄486.
[13] 摇 刘卫东, 严伟. 经济发达地区低丘缓坡土地资源合理开发利
用: 以浙江省永康市为例. 国土资源科技管理, 2007, 24(3):
1鄄5.
[20] 摇 童抗抗, 马克明. 居住鄄就业距离对交通碳排放的影响. 生态
学报, 2012, 32(10): 2975鄄2984.
[21] 摇 白杨, 郑华, 庄长伟, 欧阳志云, 徐卫华. 白洋淀流域生态系
统服务评估及其调控. 生态学报, 2013, 33(3): 711鄄717.
[22] 摇 尹海伟, 孔繁花, 祈毅, 王红扬, 周艳妮, 秦正茂. 湖南省城
市群生态网络构建与优化. 生态学报, 2011, 31 ( 10 ):
2863鄄2874.
[23] 摇 肖潇, 张捷, 卢俊宇, 钟士恩, 尹立杰. 旅游交通碳排放的空
间结构与情景分析. 生态学报, 2012, 32(23): 7540鄄7548.
[24] 摇 李鹏, 姜鲁光, 封志明, 于秀波. 生态系统服务竞争与协同研
究进展. 生态学报, 2012, 32(16): 5219鄄5229.
[25] 摇 谢苗苗,王仰麟, 付梅臣. 城市地表温度热岛影响因素研究进
展. 地理科学进展, 2011, 30(1): 35鄄41.
[26] 摇 巩杰, 谢余初, 孙朋, 颉耀文. 乡镇尺度金塔绿洲时空格局变
化. 生态学报, 2013, 33(11): 3470鄄3479.
7913摇 12期 摇 摇 摇 刘焱序摇 等:基于 OWA的低丘缓坡建设开发适宜性评价———以云南大理白族自治州为例 摇
粤悦栽粤 耘悦韵蕴韵郧陨悦粤 杂陨晕陨悦粤 灾燥造援猿源袁晕燥援员圆 允怎灶藻袁圆园员源渊杂藻皂蚤皂燥灶贼澡造赠冤
悦韵晕栽耘晕栽杂
阅藻增藻造燥责皂藻灶贼 澡蚤泽贼燥则赠 葬灶凿 枣怎贼怎则藻 则藻泽藻葬则糟澡 责则蚤燥则蚤贼蚤藻泽 燥枣 造葬灶凿泽糟葬责藻 藻糟燥造燥早赠 蚤灶 悦澡蚤灶葬
悦匀耘晕 蕴蚤凿蚤灶早袁 蕴陨 载蚤怎扎澡藻灶袁 云哉 月燥躁蚤藻袁 藻贼 葬造 渊猿员圆怨冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
砸藻泽藻葬则糟澡 则藻增蚤藻憎 燥灶 藻枣枣藻糟贼泽 燥枣 怎则遭葬灶 造葬灶凿泽糟葬责藻 责葬贼贼藻则灶 糟澡葬灶早藻泽 燥灶 憎葬贼藻则 藻灶增蚤则燥灶皂藻灶贼 匀哉粤晕郧 杂澡怎燥袁郧哉韵 匝蚤灶早澡葬蚤 渊猿员源圆冤噎噎噎噎
粤 则藻增蚤藻憎 燥枣 皂怎造贼蚤枣怎灶糟贼蚤燥灶葬造 造葬灶凿泽糟葬责藻 栽粤晕郧 匝蚤葬灶袁 阅陨晕郧 杂澡藻灶早赠葬灶 渊猿员缘员冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
悦燥皂责葬糟贼 凿藻增藻造燥责皂藻灶贼 燥枣 泽责葬糟藻鄄造蚤皂蚤贼藻凿 糟蚤贼赠院 葬 糟葬泽藻 泽贼怎凿赠 燥枣 载蚤葬皂藻灶 陨泽造葬灶凿 匀哉粤晕郧 杂澡怎燥袁 郧哉韵 匝蚤灶早澡葬蚤袁 栽粤晕郧 蕴蚤灶葬 渊猿员缘愿冤噎噎噎
粤灶葬造赠泽蚤泽 燥枣 枣葬糟贼燥则泽 糟燥灶贼则蚤遭怎贼蚤灶早 贼燥 贼澡藻 糟燥燥造蚤灶早 藻枣枣藻糟贼泽 燥枣 孕怎则责造藻 酝燥怎灶贼葬蚤灶 云燥则藻泽贼 孕葬则噪
再粤晕 宰藻蚤躁蚤葬燥袁 运韵晕郧 云葬灶澡怎葬袁 再陨晕 匀葬蚤憎藻蚤袁 藻贼 葬造 渊猿员远怨冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
陨皂责葬糟贼泽 燥枣 泽贼则怎糟贼怎则藻 糟澡葬则葬糟贼藻则蚤泽贼蚤糟泽 燥灶 贼澡藻 贼澡藻则皂葬造 藻灶增蚤则燥灶皂藻灶贼 藻枣枣藻糟贼 燥枣 糟蚤贼赠 责葬则噪泽
云耘晕郧 再怎藻赠蚤袁 匀哉 栽葬灶早葬燥袁 在匀粤晕郧 蕴蚤曾蚤葬燥 渊猿员苑怨冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
杂怎蚤贼葬遭蚤造蚤贼赠 葬泽泽藻泽泽皂藻灶贼 枣燥则 遭怎蚤造凿蚤灶早 造葬灶凿 糟燥灶泽燥造蚤凿葬贼蚤燥灶 燥灶 早藻灶贼造藻 澡蚤造造泽蚤凿藻 遭葬泽藻凿 燥灶 韵宰粤 燥责藻则葬贼燥则院葬 糟葬泽藻 蚤灶 阅葬造蚤 月葬蚤 晕葬贼蚤燥灶葬造蚤贼赠
月燥则燥怎早澡 蚤灶 再怎灶灶葬灶袁 悦澡蚤灶葬 蕴陨哉 再葬灶曾怎袁 孕耘晕郧 允蚤葬灶袁 匀粤晕 再蚤灶葬灶袁 藻贼 葬造 渊猿员愿愿冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
韵责贼蚤皂蚤扎葬贼蚤燥灶 燥枣 造葬灶凿 怎泽藻 责葬贼贼藻则灶 遭葬泽藻凿 燥灶 藻糟燥鄄泽藻糟怎则蚤贼赠院 葬 糟葬泽藻 泽贼怎凿赠 蚤灶 贼澡藻 澡怎葬灶早枣怎糟澡怎葬灶 憎葬贼藻则泽澡藻凿
再哉 云藻灶早袁 蕴陨 载蚤葬燥遭蚤灶早袁 宰粤晕郧 匀燥灶早 渊猿员怨愿冤
噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
粤灶葬造赠泽蚤泽 燥枣 造葬灶凿 怎泽藻 辕 糟燥增藻则 糟澡葬灶早藻 枣则燥皂 圆园园园 贼燥 圆园员园 葬灶凿 蚤贼泽 凿则蚤增蚤灶早 枣燥则糟藻泽 蚤灶 酝葬灶葬泽 砸蚤增藻则 月葬泽蚤灶袁 载蚤灶躁蚤葬灶早
蕴陨哉 允蚤灶憎藻蚤袁 允陨晕 栽蚤葬灶贼蚤葬灶袁 蕴陨哉 郧怎燥澡怎葬袁 藻贼 葬造 渊猿圆员员冤
噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
阅赠灶葬皂蚤糟 贼则藻灶凿 葬灶葬造赠泽蚤泽 燥枣 造葬灶凿 怎泽藻 糟澡葬灶早藻 蚤灶 贼澡藻 早葬灶躁蚤葬灶早 怎责泽贼则藻葬皂 憎葬贼藻则泽澡藻凿 遭赠 怎泽蚤灶早 砸杂 葬灶凿 郧陨杂 贼藻糟澡灶蚤择怎藻泽
蕴哉 再葬灶枣藻蚤袁 孕耘晕郧 云葬灶早袁 宰粤晕 再怎灶袁 藻贼 葬造 渊猿圆圆源冤
噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
蕴葬灶凿泽糟葬责藻 责葬贼贼藻则灶 凿赠灶葬皂蚤糟泽 葬灶凿 凿则蚤增蚤灶早 枣燥则糟藻泽 葬灶葬造赠泽蚤泽 蚤灶 贼澡藻 杂葬灶躁蚤葬灶早 孕造葬蚤灶 枣则燥皂 员怨缘源 贼燥 圆园员园
蕴陨哉 允蚤责蚤灶早袁 在匀粤韵 阅葬灶凿葬灶袁 栽陨粤晕 载怎藻扎澡蚤袁 藻贼 葬造 渊猿圆猿源冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
韵责贼蚤皂蚤扎蚤灶早 贼澡藻燥则赠 葬灶凿 糟葬泽藻 泽贼怎凿蚤藻泽 燥枣 糟怎造贼蚤增葬贼藻凿 泽造燥责藻 造葬灶凿 蚤灶 贼澡藻 糟藻灶贼藻则 燥枣 贼澡则藻藻 早燥则早藻泽 则藻泽藻则增燥蚤则 葬则藻葬 遭葬泽藻凿 燥灶 责葬贼糟澡鄄泽糟葬造藻 造葬灶凿
藻增葬造怎葬贼蚤燥灶 宰粤晕郧 再燥灶早赠葬灶袁蕴陨 再葬灶早遭蚤灶早袁 杂匀粤韵 允蚤灶早葬灶袁藻贼 葬造 渊猿圆源缘冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
杂责葬贼蚤葬造鄄贼藻皂责燥则葬造 藻增燥造怎贼蚤燥灶 燥枣 造葬灶凿 怎泽藻 葬灶凿 造葬灶凿泽糟葬责藻 责葬贼贼藻则灶 燥枣 贼澡藻 皂燥怎灶贼葬蚤灶鄄遭葬泽蚤灶 泽赠泽贼藻皂 蚤灶 郧怎蚤扎澡燥怎 孕则燥增蚤灶糟藻
蕴陨 再葬灶早遭蚤灶早袁再粤韵 再怎葬灶憎藻灶袁载陨耘 允蚤灶早袁 藻贼 葬造 渊猿圆缘苑冤
噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
杂责葬贼蚤燥鄄贼藻皂责燥则葬造 泽蚤皂怎造葬贼蚤燥灶 燥枣 造葬灶凿 糟燥增藻则 泽糟藻灶葬则蚤燥泽 蚤灶 泽燥怎贼澡憎藻泽贼藻则灶 燥枣 悦澡蚤灶葬 蕴陨 允蚤灶早袁 云粤晕 在藻皂藻灶早袁 再哉耘 栽蚤葬灶曾蚤葬灶早 渊猿圆远远冤噎噎噎噎
郧则葬凿蚤藻灶贼 葬灶葬造赠泽蚤泽 燥枣 凿则赠 增葬造造藻赠 燥枣 酝蚤灶躁蚤葬灶早 砸蚤增藻则 造葬灶凿泽糟葬责藻 责葬贼贼藻则灶袁 遭葬泽藻凿 燥灶 皂燥增蚤灶早 憎蚤灶凿燥憎 皂藻贼澡燥凿
在匀粤晕郧 蕴蚤灶早造蚤灶早袁 在匀粤韵 再燥灶早澡怎葬袁 再陨晕 杂澡葬袁 藻贼 葬造 渊猿圆苑远冤
噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
杂贼怎凿赠 燥灶 泽责葬贼蚤燥鄄贼藻皂责藻则葬造 糟澡葬灶早藻 燥枣 栽蚤遭藻贼葬灶 粤灶贼藻造燥责藻忆泽 澡葬遭蚤贼葬贼 遭葬泽藻凿 燥灶 增藻早藻贼葬贼蚤燥灶 糟燥增藻则葬早藻
在匀粤韵 匀葬蚤凿蚤袁 蕴陨哉 杂澡蚤造蚤葬灶早袁 阅韵晕郧 杂澡蚤噪怎蚤袁 藻贼 葬造 渊猿圆愿缘冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
杂责葬贼蚤葬造 澡藻贼藻则燥早藻灶藻蚤贼赠 燥枣 泽燥蚤造 皂蚤糟则燥遭蚤葬造 遭蚤燥皂葬泽泽 糟葬则遭燥灶袁 灶蚤贼则燥早藻灶袁 葬灶凿 责澡燥泽责澡燥则怎泽 蚤灶 泽造燥责蚤灶早 枣蚤藻造凿 蚤灶 葬 早则燥早藻 运葬则泽贼 则藻早蚤燥灶袁
杂燥怎贼澡憎藻泽贼 悦澡蚤灶葬 云粤晕 云怎躁蚤灶早袁 匀哉粤晕郧 郧怎燥择蚤灶袁 杂韵晕郧 栽燥灶早择蚤灶早袁 藻贼 葬造 渊猿圆怨猿冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
悦澡葬则葬糟贼藻则蚤泽贼蚤糟泽 燥枣 泽燥蚤造 皂蚤糟则燥遭蚤葬造 责燥责怎造葬贼蚤燥灶泽 葬灶凿 遭蚤燥皂葬泽泽 怎灶凿藻则 凿蚤枣枣藻则藻灶贼 藻糟燥泽赠泽贼藻皂泽 蚤灶 葬 糟葬灶赠燥灶 噪葬则泽贼 则藻早蚤燥灶
栽粤晕 匝蚤怎躁蚤灶袁 杂韵晕郧 栽燥灶早择蚤灶早袁 孕耘晕郧 宰葬灶曾蚤葬袁 藻贼 葬造 渊猿猿园圆冤
噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
杂责葬贼蚤葬造 葬灶凿 贼藻皂责燥则葬造 凿赠灶葬皂蚤糟泽 燥枣 造葬灶凿 怎泽藻 葬灶凿 蚤贼泽 蚤灶枣造怎藻灶糟藻 燥灶 藻糟燥泽赠泽贼藻皂 泽藻则增蚤糟藻 增葬造怎藻 蚤灶 再葬灶早贼扎藻 砸蚤增藻则 阅藻造贼葬
蕴陨哉 郧怎蚤造蚤灶袁在匀粤晕郧 蕴怎燥糟澡藻灶早袁在匀粤晕郧 匝蚤葬灶 渊猿猿员员冤
噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
耘增葬造怎葬贼蚤燥灶 燥枣 贼燥怎则蚤泽皂 凿赠灶葬皂蚤糟 造葬灶凿泽糟葬责藻 葬造燥灶早 匝蚤灶早澡葬蚤鄄栽蚤遭藻贼 则葬蚤造憎葬赠 遭葬泽藻凿 燥灶 贼澡藻 增蚤泽怎葬造 糟燥则则蚤凿燥则
在匀粤晕郧 砸怎蚤赠蚤灶早袁载陨 允蚤葬灶糟澡葬燥袁再粤韵 再怎造燥灶早袁藻贼 葬造 渊猿猿圆园冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
粤 泽贼怎凿赠 燥枣 泽责葬贼蚤葬造 藻增燥造怎贼蚤燥灶 糟澡葬则葬糟贼藻则蚤泽贼蚤糟泽 燥枣 则怎则葬造 泽藻贼贼造藻皂藻灶贼泽 葬灶凿 蚤灶枣造怎藻灶糟藻泽 燥枣 造葬灶凿泽糟葬责藻 责葬贼贼藻则灶泽 燥灶 贼澡藻蚤则 凿蚤泽贼则蚤遭怎贼蚤燥灶 怎泽蚤灶早 郧陨杂
葬灶凿 砸杂 砸耘晕 孕蚤灶早袁 匀韵晕郧 月怎贼蚤灶早袁 蕴陨哉 再蚤灶袁藻贼 葬造 渊猿猿猿员冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
粤泽泽藻泽泽蚤灶早 贼澡藻 藻糟燥泽赠泽贼藻皂 糟燥灶泽藻则增葬贼蚤燥灶 泽贼葬贼怎泽 葬灶凿 责则蚤燥则蚤贼赠院 葬 糟葬泽藻 泽贼怎凿赠 枣则燥皂 允蚤葬灶早曾蚤 孕则燥增蚤灶糟藻袁 悦澡蚤灶葬
云粤晕 晕葬蚤择蚤灶早袁 在匀粤晕郧 再怎曾蚤灶袁 蕴譈 再蚤澡藻袁 藻贼 葬造 渊猿猿源员冤
噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
栽澡藻 蚤皂责葬糟贼 燥枣 泽葬造贼 皂葬则泽澡 糟澡葬灶早藻 燥灶 泽藻凿蚤皂藻灶贼 葬糟糟怎皂怎造葬贼蚤燥灶 葬灶凿 憎葬增藻 葬贼贼藻灶怎葬贼蚤燥灶 葬贼 贼澡藻 耘葬泽贼 悦澡燥灶早皂蚤灶早 陨泽造葬灶凿
砸耘晕 蕴蚤灶躁蚤灶早袁 蕴陨 载蚤怎扎澡藻灶袁 再粤晕郧 杂澡蚤造怎灶袁 藻贼 葬造 渊猿猿缘园冤
噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
蕴葬灶凿泽糟葬责藻 糟造葬泽泽蚤枣蚤糟葬贼蚤燥灶 泽赠泽贼藻皂 遭葬泽藻凿 燥灶 糟造蚤皂葬贼藻袁 造葬灶凿枣燥则皂袁 藻糟燥泽赠泽贼藻皂院 葬 糟葬泽藻 泽贼怎凿赠 燥枣 载蚤灶躁蚤葬灶早 葬则藻葬
杂匀陨 匝蚤灶早凿燥灶早袁 宰粤晕郧 在澡蚤袁 匀耘 蕴燥灶早皂藻蚤袁 藻贼 葬造 渊猿猿缘怨冤
噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
粤灶葬造赠泽蚤泽 燥枣 造葬灶凿泽糟葬责藻 责葬贼贼藻则灶 藻增燥造怎贼蚤燥灶 糟澡葬则葬糟贼藻则蚤泽贼蚤糟 蚤灶 贼澡藻 澡蚤造造赠 葬灶凿 早怎造造赠 葬则藻葬 燥枣 造燥藻泽泽 责造葬贼藻葬怎院 葬 糟葬泽藻 泽贼怎凿赠 蚤灶 再葬灶忆葬灶 悦蚤贼赠袁
杂澡葬葬灶曾蚤 孕则燥增蚤灶糟藻 在匀韵晕郧 蕴蚤灶葬袁 在匀粤韵 宰藻灶憎怎袁 蕴灾 再蚤澡藻袁 藻贼 葬造 渊猿猿远愿冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
粤灶葬造赠泽蚤泽 燥枣 贼澡藻 糟澡葬则葬糟贼藻则蚤泽贼蚤糟泽 燥枣 葬早则燥鄄造葬灶凿泽糟葬责藻 澡藻贼藻则燥早藻灶藻蚤贼赠 怎灶凿藻则 贼澡藻 凿蚤枣枣藻则藻灶贼 凿蚤泽贼怎则遭葬灶糟藻泽院葬 糟葬泽藻 泽贼怎凿赠 燥枣 郧燥灶早赠蚤 悦蚤贼赠
在匀粤晕郧 载蚤葬燥赠葬灶早袁 蕴陨粤晕郧 郧怎燥枣怎袁 阅陨晕郧 杂澡藻灶早赠葬灶 渊猿猿苑愿冤
噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
悦造葬泽泽蚤枣蚤糟葬贼蚤燥灶 葬灶凿 燥则凿蚤灶葬贼蚤燥灶 燥枣 早则葬泽泽造葬灶凿 造葬灶凿泽糟葬责藻 蚤灶 贼澡藻 杂澡葬灶曾蚤 孕造葬贼藻葬怎
在匀粤晕郧 载蚤葬灶责蚤灶早袁 蕴陨 在澡蚤择蚤灶袁 宰粤晕郧 酝藻灶早遭藻灶袁 藻贼 葬造 渊猿猿愿远冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
粤灶葬造赠泽蚤泽 燥枣 枣葬糟贼燥则泽 葬枣枣藻糟贼蚤灶早 皂燥怎灶贼葬蚤灶燥怎泽 造葬灶凿 泽怎则枣葬糟藻 贼藻皂责藻则葬贼怎则藻 蚤灶 贼澡藻 泽怎皂皂藻则院 葬 糟葬泽藻 泽贼怎凿赠 燥增藻则 酝燥怎灶贼 栽葬蚤
杂哉晕 悦澡葬灶早枣藻灶早袁 运韵晕郧 云葬灶澡怎葬袁 再陨晕 匀葬蚤憎藻蚤袁 藻贼 葬造 渊猿猿怨远冤
噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
砸藻泽藻葬则糟澡 燥灶 泽责葬贼蚤葬造 凿蚤泽贼则蚤遭怎贼蚤燥灶 葬灶凿 蚤灶枣造怎藻灶糟蚤灶早 枣葬糟贼燥则 燥枣 泽燥蚤造 皂燥蚤泽贼怎则藻 蚤灶 贼赠责蚤糟葬造 凿藻责则藻泽泽蚤燥灶 葬则藻葬 燥枣 噪葬则泽贼 则藻早蚤燥灶
在匀粤晕郧 允蚤早怎葬灶早袁 杂哉 再蚤则燥灶早袁 悦匀耘晕 匀燥灶早泽燥灶早袁 藻贼 葬造 渊猿源园缘冤
噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
蕴葬灶凿泽糟葬责藻 澡藻贼藻则燥早藻灶藻蚤贼赠 燥枣 皂燥怎灶贼葬蚤灶燥怎泽 葬灶凿 澡蚤造造赠 葬则藻葬 蚤灶 贼澡藻 憎藻泽贼藻则灶 匀藻灶葬灶 孕则燥增蚤灶糟藻 遭葬泽藻凿 燥灶 皂燥增蚤灶早 憎蚤灶凿燥憎 皂藻贼澡燥凿
蕴陨 阅燥灶早噪藻袁阅陨晕郧 杂澡藻灶早赠葬灶袁蕴陨粤晕郧 郧怎燥枣怎袁藻贼 葬造 渊猿源员源冤
噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
栽则藻灶凿泽 蚤灶 增藻早藻贼葬贼蚤燥灶 葬灶凿 贼澡藻蚤则 则藻泽责燥灶泽藻泽 贼燥 糟造蚤皂葬贼藻 葬灶凿 贼燥责燥早则葬责澡赠 蚤灶 灶燥则贼澡憎藻泽贼 郧怎葬灶早曾蚤
栽韵晕郧 载蚤葬燥憎藻蚤袁 宰粤晕郧 运藻造蚤灶袁再哉耘 再怎藻皂蚤灶袁 藻贼 葬造 渊猿源圆缘冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
蕴葬灶凿泽糟葬责藻 责葬贼贼藻则灶 葬灶葬造赠泽蚤泽 燥灶 糟澡葬灶早藻 燥枣 枣则葬糟贼蚤燥灶葬造 增藻早藻贼葬贼蚤燥灶 糟燥增藻则 遭藻贼憎藻藻灶 噪葬则泽贼 葬灶凿 灶燥鄄噪葬则泽贼 葬则藻葬泽院 葬 糟葬泽藻 泽贼怎凿赠 蚤灶 匀藻糟澡蚤
阅蚤泽贼则蚤糟贼袁 郧怎葬灶早曾蚤 在澡怎葬灶早 粤怎贼燥灶燥皂燥怎泽 砸藻早蚤燥灶 宰粤晕郧 酝蚤灶早糟澡燥灶早袁 宰粤晕郧 载蚤扎澡蚤袁 蕴陨粤晕郧 在澡葬燥曾蚤燥灶早袁 藻贼 葬造 渊猿源猿缘冤噎噎噎噎噎
酝怎造贼蚤鄄泽糟葬造藻 藻枣枣藻糟贼泽 枣燥则 造葬灶凿泽糟葬责藻 皂藻贼则蚤糟泽 葬灶凿 泽责藻糟蚤藻泽 凿蚤增藻则泽蚤贼赠 怎灶凿藻则 贼澡藻 凿蚤枣枣藻则藻灶贼 凿蚤泽贼怎则遭葬灶糟藻院 葬 糟葬泽藻 泽贼怎凿赠 燥枣 郧燥灶早赠蚤 悦蚤贼赠
阅韵晕郧 悦怎蚤枣葬灶早袁 蕴陨粤晕郧 郧怎燥枣怎袁 阅陨晕郧 杂澡藻灶早赠葬灶袁 藻贼 葬造 渊猿源源源冤
噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
杂责葬贼蚤葬造 澡藻贼藻则燥早藻灶藻蚤贼赠 燥枣 泽燥蚤造 燥则早葬灶蚤糟 糟葬则遭燥灶 葬灶凿 贼燥贼葬造 灶蚤贼则燥早藻灶 糟燥灶糟藻灶贼则葬贼蚤燥灶泽 蚤灶 葬 蕴蚤贼澡燥糟葬则责怎泽 早造葬遭藻则鄄悦赠糟造燥遭葬造葬灶燥责泽蚤泽 早造葬怎糟葬
藻增藻则早则藻藻灶 遭则燥葬凿造藻葬增藻凿 枣燥则藻泽贼 再粤晕郧 阅葬灶袁 载陨粤晕郧 宰藻灶澡怎葬袁 云粤晕郧 载蚤袁 藻贼 葬造 渊猿源缘圆冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
栽澡藻 糟澡葬则葬糟贼藻则蚤泽贼蚤糟泽 葬灶凿 则藻早藻灶藻则葬贼蚤燥灶 燥枣 贼澡藻 悦澡燥藻则燥泽责燥灶凿蚤葬泽 葬曾蚤造造葬则蚤藻泽 遭则燥葬凿鄄造藻葬增藻凿 糟燥皂皂怎灶蚤贼赠 蚤灶 贼澡藻 澡蚤造造赠 则藻早蚤燥灶 燥枣 糟藻灶贼则葬造 匀怎灶葬灶
孕则燥增蚤灶糟藻袁 悦澡蚤灶葬 再陨 匀葬燥袁 阅耘晕郧 载蚤葬灶早憎藻灶袁 载陨粤晕郧 宰藻灶澡怎葬袁 藻贼 葬造 渊猿源远猿冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
云葬糟贼燥则泽 蚤灶枣造怎藻灶糟蚤灶早 贼澡藻 泽责葬贼蚤葬造 凿蚤泽贼则蚤遭怎贼蚤燥灶 燥枣 增藻早藻贼葬贼蚤燥灶 糟葬则遭燥灶 凿藻灶泽蚤贼赠 蚤灶 噪葬则泽贼 造葬灶凿泽糟葬责藻泽 燥枣 晕燥则贼澡憎藻泽贼 郧怎葬灶早曾蚤院 葬 糟葬泽藻 泽贼怎凿赠
遭葬泽藻凿 燥灶 则葬凿蚤葬造 遭葬泽蚤泽 枣怎灶糟贼蚤燥灶 灶藻贼憎燥则噪 皂燥凿藻造 在匀粤晕郧 酝蚤灶早赠葬灶早袁 宰粤晕郧 运藻造蚤灶袁阅耘晕郧 在澡藻灶澡怎葬袁 藻贼 葬造 渊猿源苑圆冤噎噎噎噎噎噎噎
园愿源猿 摇 生摇 态摇 学摇 报摇 摇 摇 猿源卷摇
叶生态学报曳圆园员源年征订启事
叶生态学报曳是由中国科学技术协会主管袁中国生态学学会尧中国科学院生态环境研究中心主办的生态学
高级专业学术期刊袁创刊于 员怨愿员年袁报道生态学领域前沿理论和原始创新性研究成果遥 坚持野百花齐放袁百家
争鸣冶的方针袁依靠和团结广大生态学科研工作者袁探索生态学奥秘袁为生态学基础理论研究搭建交流平台袁
促进生态学研究深入发展袁为我国培养和造就生态学科研人才和知识创新服务尧为国民经济建设和发展服务遥
叶生态学报曳主要报道生态学及各分支学科的重要基础理论和应用研究的原始创新性科研成果遥 特别欢
迎能反映现代生态学发展方向的优秀综述性文章曰研究简报曰生态学新理论尧新方法尧新技术介绍曰新书评价和
学术尧科研动态及开放实验室介绍等遥
叶生态学报曳为半月刊袁大 员远开本袁圆愿园页袁国内定价 怨园元 辕册袁全年定价 圆员远园元遥
国内邮发代号院愿圆鄄苑袁国外邮发代号院酝远苑园
标准刊号院陨杂杂晕 员园园园鄄园怨猿猿摇 摇 悦晕 员员鄄圆园猿员 辕 匝
全国各地邮局均可订阅袁也可直接与编辑部联系购买遥 欢迎广大科技工作者尧科研单位尧高等院校尧图书
馆等订阅遥
通讯地址院 员园园园愿缘 北京海淀区双清路 员愿号摇 电摇 摇 话院 渊园员园冤远圆怨源员园怨怨曰 远圆愿源猿猿远圆
耘鄄皂葬蚤造院 泽澡藻灶早贼葬蚤曾怎藻遭葬燥岳 则糟藻藻泽援葬糟援糟灶摇 网摇 摇 址院 憎憎憎援藻糟燥造燥早蚤糟葬援糟灶
编辑部主任摇 孔红梅摇 摇 摇 执行编辑摇 刘天星摇 段摇 靖
生摇 态摇 学摇 报渊杂匀耘晕郧栽粤陨摇 载哉耘月粤韵冤渊半月刊摇 员怨愿员年 猿月创刊冤
第 猿源卷摇 第 员圆期摇 渊圆园员源年 远月冤
粤悦栽粤 耘悦韵蕴韵郧陨悦粤 杂陨晕陨悦粤摇渊杂藻皂蚤皂燥灶贼澡造赠袁杂贼葬则贼藻凿 蚤灶 员怨愿员冤摇灾燥造郾 猿源摇 晕燥郾 员圆 渊允怎灶藻袁 圆园员源冤
编摇 摇 辑摇 叶生态学报曳编辑部
地址院北京海淀区双清路 员愿号
邮政编码院员园园园愿缘
电话院渊园员园冤远圆怨源员园怨怨憎憎憎援藻糟燥造燥早蚤糟葬援糟灶泽澡藻灶早贼葬蚤曾怎藻遭葬燥岳 则糟藻藻泽援葬糟援糟灶
主摇 摇 编摇 王如松
主摇 摇 管摇 中国科学技术协会
主摇 摇 办摇 中国生态学学会
中国科学院生态环境研究中心
地址院北京海淀区双清路 员愿号
邮政编码院员园园园愿缘
出摇 摇 版摇
摇 摇 摇 摇 摇 地址院北京东黄城根北街 员远号
邮政编码院员园园苑员苑
印摇 摇 刷摇 北京北林印刷厂
发 行摇
地址院东黄城根北街 员远号
邮政编码院员园园苑员苑
电话院渊园员园冤远源园猿源缘远猿耘鄄皂葬蚤造院躁燥怎则灶葬造岳 糟泽责早援灶藻贼
订摇 摇 购摇 全国各地邮局
国外发行摇 中国国际图书贸易总公司
地址院北京 猿怨怨信箱
邮政编码院员园园园源源
广告经营
许 可 证摇 京海工商广字第 愿园员猿号
耘凿蚤贼藻凿 遭赠摇 耘凿蚤贼燥则蚤葬造 遭燥葬则凿 燥枣
粤悦栽粤 耘悦韵蕴韵郧陨悦粤 杂陨晕陨悦粤
粤凿凿院员愿袁杂澡怎葬灶早择蚤灶早 杂贼则藻藻贼袁匀葬蚤凿蚤葬灶袁月藻蚤躁蚤灶早 员园园园愿缘袁悦澡蚤灶葬
栽藻造院渊园员园冤远圆怨源员园怨怨
憎憎憎援藻糟燥造燥早蚤糟葬援糟灶
泽澡藻灶早贼葬蚤曾怎藻遭葬燥岳 则糟藻藻泽援葬糟援糟灶
耘凿蚤贼燥则鄄蚤灶鄄糟澡蚤藻枣摇 宰粤晕郧 砸怎泽燥灶早
杂怎责藻则增蚤泽藻凿 遭赠摇 悦澡蚤灶葬 粤泽泽燥糟蚤葬贼蚤燥灶 枣燥则 杂糟蚤藻灶糟藻 葬灶凿 栽藻糟澡灶燥造燥早赠
杂责燥灶泽燥则藻凿 遭赠摇 耘糟燥造燥早蚤糟葬造 杂燥糟蚤藻贼赠 燥枣 悦澡蚤灶葬
砸藻泽藻葬则糟澡 悦藻灶贼藻则 枣燥则 耘糟燥鄄藻灶增蚤则燥灶皂藻灶贼葬造 杂糟蚤藻灶糟藻泽袁 悦粤杂
粤凿凿院员愿袁杂澡怎葬灶早择蚤灶早 杂贼则藻藻贼袁匀葬蚤凿蚤葬灶袁月藻蚤躁蚤灶早 员园园园愿缘袁悦澡蚤灶葬
孕怎遭造蚤泽澡藻凿 遭赠摇 杂糟蚤藻灶糟藻 孕则藻泽泽
粤凿凿院员远 阅燥灶早澡怎葬灶早糟澡藻灶早早藻灶 晕燥则贼澡 杂贼则藻藻贼袁
月藻蚤躁蚤灶早摇 员园园苑员苑袁悦澡蚤灶葬
孕则蚤灶贼藻凿 遭赠摇 月藻蚤躁蚤灶早 月藻蚤 蕴蚤灶 孕则蚤灶贼蚤灶早 匀燥怎泽藻袁
月藻蚤躁蚤灶早 员园园园愿猿袁悦澡蚤灶葬
阅蚤泽贼则蚤遭怎贼藻凿 遭赠摇 杂糟蚤藻灶糟藻 孕则藻泽泽
粤凿凿院员远 阅燥灶早澡怎葬灶早糟澡藻灶早早藻灶 晕燥则贼澡
杂贼则藻藻贼袁月藻蚤躁蚤灶早 员园园苑员苑袁悦澡蚤灶葬
栽藻造院渊园员园冤远源园猿源缘远猿
耘鄄皂葬蚤造院躁燥怎则灶葬造岳 糟泽责早援灶藻贼
阅燥皂藻泽贼蚤糟 摇 摇 粤造造 蕴燥糟葬造 孕燥泽贼 韵枣枣蚤糟藻泽 蚤灶 悦澡蚤灶葬
云燥则藻蚤早灶 摇 摇 悦澡蚤灶葬 陨灶贼藻则灶葬贼蚤燥灶葬造 月燥燥噪 栽则葬凿蚤灶早
悦燥则责燥则葬贼蚤燥灶
粤凿凿院孕援韵援月燥曾 猿怨怨 月藻蚤躁蚤灶早 员园园园源源袁悦澡蚤灶葬
摇 陨杂杂晕 员园园园鄄园怨猿猿悦晕 员员鄄圆园猿员 辕 匝 国内外公开发行 国内邮发代号 愿圆鄄苑 国外发行代号 酝远苑园 定价 怨园郾 园园元摇