全 文 ：
摇 摇 摇 摇 摇 生 态 学 报
摇 摇 摇 摇 摇 摇 摇 渊杂匀耘晕郧栽粤陨 载哉耘月粤韵冤
摇 摇 第 猿源卷 第 猿期摇 摇 圆园员源年 圆月摇 渊半月刊冤
目摇 摇 次
前沿理论与学科综述
灾后生态恢复评价研究进展 刘孝富袁王文杰袁李摇 京袁等 渊缘圆苑冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
基于生态能量视角的我国小水电可持续性分析 庞明月袁张力小袁王长波 渊缘猿苑冤噎噎噎噎噎噎噎噎噎噎噎噎噎
个体与基础生态
北部湾北部海域夏季微型浮游动物对浮游植物的摄食压力 马摇 璐袁曹文清袁张文静袁等 渊缘源远冤噎噎噎噎噎噎
鲶鱼和胡子鲶的两性异形与雌性个体生育力 樊晓丽袁林植华袁丁先龙袁等 渊缘缘缘冤噎噎噎噎噎噎噎噎噎噎噎噎
环境温度对白头鹎代谢产热和蒸发失水的影响 林摇 琳袁曹梦婷袁胡益林袁等 渊缘远源冤噎噎噎噎噎噎噎噎噎噎噎
灌溉对沙拐枣幼苗生长及氮素利用的影响 黄彩变袁曾凡江袁雷加强袁等 渊缘苑圆冤噎噎噎噎噎噎噎噎噎噎噎噎噎
孕粤匀泽污染土壤植物修复对酶活性的影响 朱摇 凡袁洪湘琦袁闫文德袁等 渊缘愿员冤噎噎噎噎噎噎噎噎噎噎噎噎噎
基于修正 砸陨杂运耘模型的重庆岩溶地区地下水脆弱性评价 魏兴萍袁蒲俊兵袁赵纯勇 渊缘愿怨冤噎噎噎噎噎噎噎噎
排水沟蓄水条件下农田与排水沟水盐监测 潘延鑫袁罗摇 纨袁贾忠华袁等 渊缘怨苑冤噎噎噎噎噎噎噎噎噎噎噎噎噎
种群尧群落和生态系统
高寒退化草地不同海拔梯度狼毒种群分布格局及空间关联性 高福元袁赵成章袁卓马兰草 渊远园缘冤噎噎噎噎噎噎
捕食者对空心莲子草叶甲种群的生物胁迫 刘雨芳袁李摇 菲袁桂芳艳袁等 渊远员猿冤噎噎噎噎噎噎噎噎噎噎噎噎噎
夏尧冬季南海北部浮游植物群落特征 马摇 威袁孙摇 军 渊远圆员冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
磨盘山天然次生林凋落物数量及动态 范春楠袁郭忠玲袁郑金萍袁等 渊远猿猿冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
持续干旱对樱桃根际土壤细菌数量及结构多样性影响 刘方春袁邢尚军袁马海林袁等 渊远源圆冤噎噎噎噎噎噎噎噎
随机森林算法基本思想及其在生态学中的应用要要要以云南松分布模拟为例
张摇 雷袁王琳琳袁张旭东袁等 渊远缘园冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
基于水文平衡的湿地退化驱动因子定量研究 侯摇 鹏袁申文明袁王摇 桥袁等 渊远远园冤噎噎噎噎噎噎噎噎噎噎噎噎
华北低丘山地人工林蒸散的控制因子 黄摇 辉袁孟摇 平袁张劲松袁等 渊远远苑冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
新疆天山高寒草原不同放牧管理下的 悦韵圆袁悦匀源和 晕圆韵通量特征 贺桂香袁李凯辉袁宋摇 韦袁等 渊远苑源冤噎噎噎噎
景观尧区域和全球生态
宁夏生态足迹影响因子的偏最小二乘回归分析 马明德袁马学娟袁谢应忠袁等 渊远愿圆冤噎噎噎噎噎噎噎噎噎噎噎
引黄灌区土壤有机碳密度剖面特征及固碳速率 董林林袁杨摇 浩袁于东升袁等 渊远怨园冤噎噎噎噎噎噎噎噎噎噎噎
自养微生物同化 悦韵圆的分子生态研究及同化碳在土壤中的转化 吴小红袁简摇 燕袁陈晓娟袁等 渊苑园员冤噎噎噎噎
资源与产业生态
基于能值分析法的矿区循环经济系统生态效率分析 孙玉峰袁 郭全营 渊苑员园冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎
基于 粤阅鄄粤杂模型的海岸带生态系统综合承载力评估要要要以舟山海岸带为例
苏盼盼袁叶属峰袁过仲阳袁等 渊苑员愿冤
噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
城乡与社会生态
基于增强回归树和 蕴燥早蚤泽贼蚤糟回归的城市扩展驱动力分析 李春林袁刘摇 淼袁胡远满袁等 渊苑圆苑冤噎噎噎噎噎噎噎
研究简报
陕西省不同生态区大气氮素干湿沉降的时空变异 梁摇 婷袁同延安袁林摇 文袁等 渊苑猿愿冤噎噎噎噎噎噎噎噎噎噎
不同覆盖方式对旱地果园水热特征的影响 刘小勇袁李红旭袁李建明袁等 渊苑源远冤噎噎噎噎噎噎噎噎噎噎噎噎噎
长白山苔原带土壤动物群落结构及多样性 王振海袁殷秀琴袁蒋云峰 渊苑缘缘冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
最大可允许填海面积模拟要要要厦门西海域案例研究 林琛琛袁 饶欢欢袁 刘摇 岩袁等 渊苑远远冤噎噎噎噎噎噎噎噎
学术信息与动态
圆园员猿年水文土壤学与自然资源可持续利用国际学术研讨会述评 张摇 骁袁赵文武 渊苑苑源冤噎噎噎噎噎噎噎噎噎
期刊基本参数院悦晕 员员鄄圆园猿员 辕 匝鄢员怨愿员鄢皂鄢员远鄢圆缘圆鄢扎澡鄢孕鄢 预 怨园郾 园园鄢员缘员园鄢圆怨鄢圆园员源鄄园圆
室室室室室室室室室室室室室室
封面图说院 云南松树冠要要要云南松为松科松属裸子植物袁多生长在海拔 员园园园要猿缘园园皂的高山袁喜光尧耐干旱尧耐瘠薄袁适应酸性
的红壤尧黄壤袁在其他树种不能生长的贫瘠石砾地或冲刷严重的荒山坡分布袁易于天然更新遥 主要分布于四川西南
部尧云南尧西藏东南部尧贵州西部尧广西西部袁常形成大面积纯林袁尤以云南分布最广袁故有云南松之称遥 云南松树高
可达 猿园皂袁胸径达 员皂袁树皮呈灰褐色袁叶通常 猿针一束袁鲜有两针袁球果圆锥状卵圆形袁种子近卵圆形或倒卵形遥 树干
通直袁木质轻软细密袁是优质造纸尧人造板原料袁富含松脂是云南松的重要特点之一遥
彩图及图说提供院 陈建伟教授摇 北京林业大学摇 耘鄄皂葬蚤造院 糟蚤贼藻泽援糟澡藻灶躁憎岳 员远猿援糟燥皂
第 34 卷第 3 期
2014年 2月
生 态 学 报
ACTA ECOLOGICA SINICA
Vol.34,No.3
Feb.,2014
http: / / www.ecologica.cn
基金项目:国家海洋公益资助项目 (201105006);国家社会科学基金资助项目(12BJ064);福建省自然科学基金资助项目(No. 2010J01360)
收稿日期:2012鄄09鄄04; 摇 摇 修订日期:2013鄄03鄄04
*通讯作者 Corresponding author.E鄄mail: brpeng@ xmu.edu.cn
DOI: 10.5846 / stxb201209041246
林琛琛, 饶欢欢, 刘岩, 李京梅, 江毓武,彭本荣.最大可允许填海面积模拟———厦门西海域案例研究.生态学报,2014,34(3):766鄄773.
Lin C C, Rao H H, Liu Y, Li J M, Jiang Y W, Peng B R.Modeling the total allowable area for coastal reclamation:a case study of Xiamen, China.Acta
Ecologica Sinica,2014,34(3):766鄄773.
最大可允许填海面积模拟
———厦门西海域案例研究
林琛琛1, 饶欢欢1, 刘摇 岩2, 李京梅3, 江毓武1,彭本荣1,*
(1. 福建省海陆界面生态环境重点实验室,厦门大学环境与生态学院,厦门摇 361005;
2. 国家海洋局海洋发展战略研究所,北京摇 100860; 3. 中国海洋大学经济学院,青岛摇 266100)
摘要:填海造地在解决海岸带地区发展空间短缺问题的同时,严重损害了人类赖以生存和发展的海洋与海岸带生态系统。 如何
控制大规模填海造地是我国面临的重大课题。 构建了估算海湾最大可允许填海造地面积的分析框架和模型。 填海造地的各种
收益和成本,包括环境和生态成本被系统性地量化并嵌入到框架和模型之中。 模型的估算结果可以为实施控制填海造地的约
束机制提供科技支撑。 利用建立的框架和模型对厦门西海域模拟的结果显示:当贴现率为 4.5%时,厦门西海域最大可允许填
海面积是 89.44 hm2,其面积至少应该保持在 44.52 km2;即便是使用 8%的高贴现率,最大可允许填海面积也只占规划填海面积
的 35.83%;如果采用更加重视环境与生态效益的低贴现率 2%,则厦门西海域最大可允许填海面积是 0。
关键词:最大可允许填海面积;最优填海位置;收益和成本分析;厦门
Modeling the total allowable area for coastal reclamation:a case study of Xiamen,
China摇
LIN Chenchen1, RAO Huanhuan1, LIU Yan2, LI Jingmei3, JIANG Yuwu1, PENG Benrong1,*
1 Joint Key Laboratory of Coastal Study, College of Environment and Ecology Xiamen University, Xiamen 361005, China
2 China Institute for Marine Affairs, State Oceanic Administration, Beijing 100860, China
3 College of Economics, Ocean University of China, Qingdao 266100, China
Abstract: While creating useful space for agriculture, industry, and urban development, coastal reclamation damages
considerably the marine ecosystems which human鄄kind depends on. How to control the large鄄scale coastal reclamation is a
great challenge in China. This paper presents an analytical framework and model to estimate the Total Allowable Area for
Coastal Reclamation ( TAACR) to provide scientific support for the implementation of a coastal reclamation restriction
mechanism. Various benefits and costs, including the ecological and environmental costs of coastal reclamation, are
systematically quantified in the framework. Model simulations are developed using data from West Sea of Xiamen. The
results show that the TAACR in West Sea is 89.44 hm2 at a discount rate of 4.5%. Area of West Sea must be maintained at a
minimum of 44.52 km2 so that the social well鄄being can be maximized. Under low discount rate, 2%, the TAACR is zero.
Even with high discount rate, 8%, the optimal reclamation area is 157.97 hm2, about 35.83% of the total planned area.
Key Words: total allowable reclamation; optimal reclamation area; benefit and cost analysis; Xiamen City
http: / / www.ecologica.cn
摇 摇 因其高生物生产力和便利的交通条件,海岸带
地区一直是人类活动的中心[1]。 世界大约 40%的人
口、一半的生产和消费活动都集中在仅占陆地总面
积 10%的海岸带地区[2]。 所有高人口密度的海岸带
地区都面临着发展“空间赤字冶问题,填海造地成为
海岸带地区人们解决土地短缺问题的共同方式。 在
岛屿国家如日本、韩国、新加坡等,非岛国家如荷兰、
德国、美国等均可观察到大规模填海造地活动[3]。
填海造地意味着创造了更多的农业、工业和城市发
展用地,然而这通常是以损害环境为代价的。 海洋
转变成陆地使得海洋和海岸带地区的自然属性发生
永久性改变,严重破坏了人类赖以生存的海洋与海
岸带生态系统[4鄄10]。
从 20世纪中期开始的大规模填海造地,已经使
我国约 1.34万 km2滩涂和浅海成为农业、盐场、海洋
养殖、工业和城镇发展用地[11]。 这导致近 50%的滨
海湿地和 73%的红树林丧失,进而引起水生生物栖
息地破坏、生物多样性减少和海岸带地区生产力下
降,威胁到海岸带地区可持续发展[12鄄13]。
我国目前相关控制填海造地的法律法规的实施
效果都不是很理想。 2001 年颁布的《中华人民共和
国海域使用管理法》第 33条规定的海域有偿制度为
控制填海造地提供了经济刺激手段。 但现有的填海
造地海域使用金标准偏低影响了其作用的发挥[14];
而第 35条和 36条关于用海项目海域使用金减免条
款限制了经济刺激手段控制政府所实施的大型填海
造地项目的功能;第 18条规定的大型用海项目中央
政府审批制度要求面积超过 50 公顷填海项目必须
得到国务院批准。 在实践中,大型填海造地项目的
开发者往往把整个大项目分成若干个小项目来规避
中央政府审批制度。 1982 年颁布(1999 年修订)的
《中华人民共和海洋环境保护法》第 43 条和 47 条规
定海岸工程和海洋工程建设项目,包括填海造地项
目必须进行项目环境影响评价(EIA)。 项目环境影
响评价往往只考虑单一项目的环境影响而没有考虑
到所有相关项目的累积环境影响。 单一项目对环境
影响往往比较小因而容易通过环境影响评价审查。
更为重要的是,大多数的填海造地项目都源于地方
政府的发展规划。 所以尽管中央政府再三要求严格
限制填海造地,沿海省市仍然提出了大量的填海造
地计划。
随着新一轮的沿海发展战略的开展淤,到 2020
年我国还有 5780 km2填海造地需求,这无疑将对海
洋与海岸带生态环境造成严重影响[15]。 因此急需
探索控制填海造地的有效途径,以协调填海造地需
求和海洋与海岸生态系统保护的关系。 其中的一种
解决方案是建立类似国家耕地保护政策的填海造地
“红线制度冶 [15]。 填海造地红线制度意味着在充分
考虑填海造地的多重收益和成本的基础上,确定每
一个具体海湾必须保留的最小海域面积于。 并以法
律法规的形式将这一最小必须保留的海域面积确定
下来。 任何试图突破所设定的临界值的个人和单
位,包括政府,必须首先修改法律。 而修改法律这个
既困难而又漫长的过程可以有效地控制大规模填海
造地。
在海湾要实施填海造地红线制度,关键是在权
衡经济发展和环境保护的基础上,估算出该海湾最
大可允许的填海造地的面积。 论文第二部分构建了
估算最大可允许的填海面积的分析框架和模型;第
三部分应用建立的方法和模型对厦门西海域最大可
允许填海面积进行了研究;结论和讨论包括在论文
的第四部分。
1摇 研究方法和模型
1.1摇 基本模型
理论上,估算最大可允许的填海造地的面积的
基本思路是直观的。 填海造地所形成的土地可用于
工农业生产以及城市发展。 这些可以促进社会福利
的改善,增加社会收益。 同时,填海造地涉及到一系
列成本,包括环境成本,使得社会福利的减少。 从整
个社会角度来看,填海造地规划者的目标是通过选
取最优的填海造地位置和面积来最大化社会净收益
(总收益扣除总成本)。
社会净收益最大化问题可用以下非线性规划模
型来描述:
Max. NB =移
n
i = 1
[pi 伊 f(xi) - c(xi)] - 移
n
i = 1
移
m
k = 1
Ck(xi)
(1)
767摇 3期 摇 摇 摇 林琛琛摇 等:最大可允许填海面积模拟———厦门西海域案例研究 摇
淤
于
http: / / news. xinhuanet. com / local / 2009鄄06 / 24 / content _ 11573216.
htm
考虑到海湾环境比开放海域的环境更加脆弱,所以本论文主要
研究海湾最大可允许填海面积
http: / / www.ecologica.cn
摇 摇 摇 摇 摇 摇 摇 摇 摇 s.t.
hi(xi) 臆 himax (2)
式中,NB表示海岸围海造地的净效益; i ( = 1, 2,
3,…, n) 代表填海位置;pi代表位置 i填海形成的土
地上生产的产品的价格;f( xi)是位置 i 填海形成的
土地的生产函数;c(xi) 是位置 i 填海形成的土地生
产的成本函数;k = (1, 2, 3,…,m)代表填海造地成
本项;Ck(xi)表示位置 i 填海第 k 种成本的成本函
数;函数 hi是位置 i 填海的约束条件,如填海规模约
束、关键生境面积不减少等。 它是关于 x 的线性
函数。
以上非线性规划问题的最优解,x*i (在位置 i 的
填海面积),必须满足 Kuhn鄄Tucker 条件。 但是,
Kuhn鄄Tucker条件是潜在最优解的必要条件,而不是
全局最优解的充分条件。 根据微观经济学一般理
论,填海造地总成本函数以及所有单项成本函数为
拟凸函数,填海造地收益函数为凹函数。 这样目标
函数是拟凹的。 因此,拉格朗日量 Hessian 矩阵是半
负定的,即存在唯一最大解。
方程 ( 3 )—( 6) 列出了所需的 Kuhn鄄Tucker
条件。
pi 伊
鄣f(x*i )
鄣xi
-
鄣c(x*i )
鄣xi
- 移
m
k = 1
鄣Ck(x*i )
鄣xi
-
姿 i
鄣h(x*i )
鄣xi
= 0
(3)
hi(x*i ) 臆 himax (4)
姿 ihi(x*i ) = 0 (5)
姿 i 逸 0摇 摇 摇 摇 ( = 0若 hi(x*i ) < himax) (6)
方程(3)中是最优性条件,方程(4)是可行性条
件,方程(5)是松驰互补条件,方程(6)是非负条件。
姿 i是拉格朗日因子,在经济学理论上被认为代表影
子价格。
1.2摇 填海造地收益
在位置 i,填海造地的收益体现在方程(1)的 [pi
伊 f(xi) - c(xi)] 项中。 它是填海形成土地上生产的
产品的价值与将土地转化为产品的成本之差。 在自
然资源经济学中,这一差额是土地的租金或者土地
的价格[16]。 考虑到填海形成土地占现有土地总量
的比重较小,可以合理假定填海形成土地的供给将
不会影响现有土地的市场价格。 因此填海形成土地
的收益可以用相邻的、同类型的土地的市场价格来
代表。
如上所述,填海土地能够有不同的用途。 不同
用途的土地价格是不同的。 而且,不同区域的土地
价格变化是显著的。 如果知道填海形成土地各种用
途的权重,位置 i填海形成土地价格可用方程(7)来
进行计算:
P i =移
j = J
j = 1
P ij 伊 sij (7)
式中,j= (1, 2, 3,…,J)表示土地利用方式;i = (1,
2, 3,…,n) 代表填海位置;P ij表示在位置 i中 j使用
方式的土地价格;sij位置 i 填海形成土地 j 使用方式
的权重。
1.3摇 填海造地成本
填海造地的成本可以分为 4个部分:(1)工程成
本;(2)生态系统损害成本;(3)填海导致泥沙淤积
产生的成本;(4)填海导致周边海域海洋环境容量下
降而产生的成本。 工程成本是指将被填海域变成可
以马上使用的土地而发生的成本,可通过现场调查
或者工程分析来估算。
海岸生态系统为人类提供了各种服务,比如提
供食物和材料,营养物的存储和循环,净化内陆淡水
系统的污染物,保护海岸线免遭侵蚀和风暴潮破坏,
调节全球水文和气候,以及容纳和吸收废弃物等。
填海造地将导致被填海域提供这些服务能力的消
失,并对相邻海域及海岸生态系统提供上述服务的
能力产生负面影响。 估算填海造地对生态系统服务
损害时,必须对海岸生态系统服务价值进行评估。
填海造地将改变海湾水动力条件,降低海湾纳
潮量,从而将导致航道泥沙淤积的增加,增加清淤成
本。 清淤成本可以用来估算填海造地导致淤积而产
生的成本。 海湾纳潮量减少将降低海湾水体吸收、
稀释、溶解污染物的能力。 社会必须采取措施提高
治理这些污染的能力。 增加的污染处理成本可以用
以估算因填海造地引起的环境容量下降而产生的
成本。
1.4摇 收益和成本的年价值
在中国,土地价格反映了法律规定的土地最高
使用年限的现值。 例如,法律规定商业、住宅和工业
用地的最高使用年限分别为 40、70a 和 50a。 其他的
估算,例如生态系统损害成本,通常以年金形式表
示。 为了能够比较,论文分别估算了与填海造地有
关的收益和成本的年价值。 一次性收益(如土地价
867 摇 生摇 态摇 学摇 报摇 摇 摇 34卷摇
http: / / www.ecologica.cn
格)和成本使用以下公式转化成年价值。
A = S (1
+ r) nr
(1 + r) n - 1
(8)
式中,A是收益或成本的年价值;S 是相关的一次性
收益和成本的现值;r是贴现率;n是年限。
一次性发生的填海造地的工程成本(CENG)在未
来几乎无限期时期里将带来收益。 相当于公式(8)
中的 n趋向于无穷大。 所以填海造地工程成本的年
价值(CENGA)计算公式为:
CENGA = CENG 伊 r (9)
2摇 案例研究
厦门市位于福建省东南,台湾海峡西岸。 作为
海岸带城市,厦门的社会经济发展和人们福利的改
善依赖海洋提供的资源、产品和服务。 20世纪 50 年
代以来的大规模填海造地使厦门 125.74 km2的海域
变成陆地(截止到 2005 年)。 填海造地导致 90%天
然红树林消失,多种生物资源的自然栖息地遭到破
坏或改变。 填海造地与其他人类活动一起,使得海
洋环境质量恶化、水体富营养化和相应的赤潮增加、
渔业资源及其他海洋物种下降,以及航道淤积等。
严重威胁到厦门的可持续发展。
即使在这种情况下,为应对经济快速发展、人口
增长和城市化进程的加快而出现的空间短缺问题,
有关部门仍然提出了 20km2的填海计划,其中西海
域的 7个规划填海区(从 WS1 到 WS7)规划填海面
积 4.41 km2(见图 1,表 1第 2列)。 西海域位于厦门
城市的核心部位,长期以来是厦门社会经济发展的
重心。 大规模填海造地已经使海域面积从 108km2
缩小到 45.41km2,纳潮面积减少了 58%[17]。 厦门迫
切需要研究海湾最大可允许填海造地面积及其位
置,在此基础上,寻找有效途径确保各海域最小需要
维持的面积,以维护海洋生态系统的健康和恢复力。
图 1摇 厦门西海域、同安湾填海规划区域和面积
Fig. 1 摇 location and area of planned reclamation in West Sea
and Tongan Bay of Xiamen
图中阴影区为围海规划区,数字为规划填海面积(km2)
本部分利用建立的模型对厦门西海域最大可允
许的填海面积及其位置进行了估算。 论文以 20a 长
期政府债券利率 4.5%作为基准社会贴现率[18鄄19]。
2.1摇 填海土地的收益函数
根据规划,厦门西海域填海形成的土地主要是
用于商业、工业和基础设施建设。 为了建立填海土
地的收益函数,本研究调查了与规划填海区域相邻
的商业和工业用地的价格(表 1)。 填海形成的土地
用于商业和工业的比例见表 1。 根据调查得到的数
据和公式(7)、(8),计算得到各规划的填海区形成
土地的加权平均年价格 PAi (表 1)。 从而填海造地的
收益函数为:
B i(xi) = PAi xi (10)
表 1摇 土地价格和填海造地收益
Table 1摇 Revenue of reclaimed land
区域
Location
面积 Area淤
/ hm2
不同使用的比例
Share of different uses于
商业 Commercial 工业 Industrial
土地价格
Price of land盂 / (元 / m2)
商业 Commercial 工业 Industrial
平均加权价格
Weighted average price
(Pi) / (元 / m2年)
WS1 48 0.2 0.8 2950 250 37.95
WS2 147 0.35 0.65 2950 250 56.92
WS3 22 0.45 0.55 2950 250 69.57
WS4 44 0.4 0.6 2950 250 63.25
WS5 114 0.35 0.65 2950 250 56.92
WS6 39 0.4 0.6 2950 250 63.25
WS7 27 0.75 0.25 4150 700 155.66
摇 摇 数据来源: 淤, 于 厦门市发展与改革委员会和厦门市海洋与渔业局;盂 厦门市国土与房地产管理局和实地调查(用于基础设施的土地价格
同工业用地价格)
967摇 3期 摇 摇 摇 林琛琛摇 等:最大可允许填海面积模拟———厦门西海域案例研究 摇
http: / / www.ecologica.cn
2.2摇 填海土地的成本函数
(1)工程成本
根据关于厦门市填海项目的研究[20鄄21],填海造
地的工程成本为 660 元 / m2。 利用公式(9),计算得
到工程成本的年价值为 29.7元 / m2。 因此工程造价
的成本函数 CENGAi (xi) 是:
CENGAi (xi) = 29.7xi (11)
(2)生态系统损害成本
关于厦门海洋生态系统服务价值及海洋生态损
害评估的研究较多[22鄄24],但是具体到每个海域的生
态系统服务价值的评估尚少。 厦门大学海洋与海岸
带发展研究院[25]对厦门各海域的生态系统进行了
识别,并对每一个海区的生态系统服务价值进行了
评估(表 2)。 填海造地破坏了被填海域提供的所有
这些服务。 生态系统服务价值的损失是填海造地导
致海岸带生态系统损害的成本。 从而海岸带生态系
统损害的成本函数 CEDi (xi) 为:
表 2摇 规划填海区的生态系统服务价值
Table 2摇 Value of ecosystem services of planned reclamation area
区域
Location
生态系统服务价值(ei) / (元·m-2·a-1)
Value of ecosystem services
WS1 9.34
WS2 8.61
WS3 9.46
WS4 9.46
WS5 9.46
WS6 9.46
WS7 9.62
摇 摇 数据来源:厦门大学海洋与海岸带发展研究院[25]
CEDi (xi) = eixi (12)
其中,ei是表 2 第 2 列所列的生态系统服务
价值。
(3)环境容量下降和淤积成本
福建海洋研究所和海洋环境科学国家重点实验
室(厦门大学) [26]建立了厦门湾的水动力、泥沙沉积
物以及水质模型。 利用这一模型可以模拟各种填海
造地方案及其组合造成的泥沙淤积和纳潮量的变
化。 填海造地方案的组合是根据填海土地需求的先
后顺序以累积叠加的形式设计的(表 3)。 通过模型
估算得到的不同填海方案下每年航道泥沙淤积量和
每个潮周期淤的纳潮量结果见表 3。
厦门航道清淤的平均成本是 29.7 元 / m3于。 以
单位成本乘以泥沙淤积量,可得到不同填海方案下
的清淤成本(表 3)。 不同填海方案下的环境容量的
价值用下列公式计算,计算结果见表 3:
V = 365 伊 2 伊 驻c 伊 v 伊 c (13)
式中,吟c 是西海域内部和外部 COD 的浓度差(吟c
的值是 0. 035mg / L[27] );v 是纳潮量;c 是去除 COD
的单位成本(c = 4300元 / t[28鄄29])。
使用表 3的成本估算的结果,运用普通最小二
乘法对边际成本函数进行估计。
驻C
驻x
= a 伊 x + b (14)
dCdc
dx
= 0.288x + 2.91 摇 (R2 = 0.97) (15)
( t= 13.06 P<0.001) (1.71, P<0.1)
淤厦门海域每天有 2个潮周期
于资料来源:厦门市港务局
表 3摇 不同填海方案的清淤成本和环境容量价格价值
Table 3摇 Dredging cost and value of carrying capacity of different reclamation scenarios
情景
Scenarios
方案
Schemes
填海面积
Reclamation area
/ hm2
淤积量
Siltation / m3
清淤成本
Dredging cost
/ 104元
纳潮量
Tidal volume
/ 106m3
环境容量价值
Value of environmental
carrying capacity
/ 104元
1 无填海 0 1528597 4540 416.5 8890.28
2 WS6 39 1534061 4556 414.7 8852.00
3 WS6+WS1 87 1540786 4576 412.5 8804.90
4 WS6+WS1+WS3 109 1543439 4584 411.3 8779.98
5 WS6+WS1+WS3+WS7 136 1546695 4594 409.9 8749.40
6 WS6+WS1+WS3+WS7+WS5 250 1559404 4631 404.1 8625.60
7 WS6+WS1+WS3+WS7+WS5+WS4 294 1563642 4644 402.3 8587.17
8 WS6+WS1+WS3+WS7+WS5+WS4+WS2 441 1578295 4688 396.4 8460.38
077 摇 生摇 态摇 学摇 报摇 摇 摇 34卷摇
http: / / www.ecologica.cn
dCecc
dx
= 0.839x + 5.995 摇 (R2 = 0.99) (16)
( t= 31.02 P<0.001)摇 (2.88, P<0.05)
方程(14)是边际成本函数的一般形式。 吟C 表
示增加的清淤成本或因环境容量下降的增加的清污
成本;吟x代表填海面积的增加。 方程(15)和(16)
分别是泥沙淤积边际成本函数和环境容量减少的边
际成本函数。
2.3摇 最大可允许填海面积及位置
将前文得到的填海造地收益和成本方程的信息
带入模型(1),并设计 GAMS 程序来求解模型(1)的
填海造地净收益最大化问题。 其中约束条件方程
(4)hi max是每个区域的规划填海面积(表 1)淤。 结
果表明在贴现率为 4.5%时西海域最大可允许填海
面积是 89.44 hm2(表 4)。 此最优面积占规划的填海
面积的 20.28%。 也就是西海域面积必须至少维持
在 44.52 km2才能达到社会净收益的最大化。 这意
味着从整个社会角度来看,现有的填海规划经济效
率是低下的。
最优的填海位置和填海面积如表 4 所示。 结果
表明,在 WS1 区域填海造地是不可行的,WS2 和
WS5区最大填海面积必须控制在规划填海面积的
8%以下。
表 4摇 估算结果和敏感性分析
Table 4摇 Estimation results and sensitive analysis
区域
Location
规划面积 / hm2
Planned area
最大可允许填海面积 Optimal reclamation area / hm2
4.5% 8% 2%
WS1 48 0 0 0.00
WS2 147 8.6 23 0.00
WS3 22 19.07 22 0.00
WS4 44 13.46 31.86 0.00
WS5 114 7.85 22.25 0.00
WS6 39 13.46 31.86 0.00
WS7 27 27 27 0.00
合计 441 89.44 157.97 0.00
淤规划的填海区避开了环境敏感区,所以没有关键生境面积约束
摇 摇 为检查不同贴现率下最大可允许填海面积的变
化,论文分别使用 2%和 8%两个贴现率进行了敏感
性分析。 如表 4 的第 4 列和第 5 列所示,分析结果
对贴现率非常敏感。 在 2%的低贴现率下,最大可允
许填海面积为 0。 在 8%的高贴现率下,最大可允许
填海面积是 157.97hm2。 由于海洋生态系统提供的
服务分布的时间较长,这些环境效益是下一代以及
下好几代人了利益。 使用较低的社会贴现率更能体
现下一代人的利益。 所以在制定填海造地规划时应
该采用较低的社会贴现率。 同时需要指出的是,即
使使用较高的贴现率 8%,最优填海面积也只占规划
面积的 35. 82%。 西海域的面积必须至少维持在
43郾 83 km2。
3摇 结论和讨论
填海造地在为农业、工业和城市发展创造有用
空间的同时,永久性改变了海洋和海岸环境的自然
属性,损害了人类赖以生存和发展的海洋生态系统。
决策者必需仔细权衡短期利益(满足土地需求)与长
期利益(确保海洋和海岸带生态系统可持续利用)二
者之间的关系。 通过估算最大可允许填海面积及其
位置,可以找到既能够满足当前社会经济发展的需
又可以维持海岸生态系统的健康和恢复力的解决方
案。 论文构建的分析框架和模型涵盖了填海造地的
多重利益及成本部分。 从整个社会角度出发来考虑
社会净福利,在分析框架中同时包含填海造地多种
用途和环境生态效应、内部效应和外部效应是十分
必要的。
估算结果显示当贴现率为 4.5%,厦门西海域最
大可允许填海面积是 89.44hm2。 即便是使用 8%的
高贴现率,最大可允许填海面积也只占规划填海面
积的 35.83%。 如果采用低贴现率 2%,则最大可允
许填海面积是 0。 为了保持海岸生态系统为人类提
供各种服务的能力,厦门市应重新审视目前的填海
造地规划,使西海域面积至少 44.52 km2。 采用 2%
177摇 3期 摇 摇 摇 林琛琛摇 等:最大可允许填海面积模拟———厦门西海域案例研究 摇
http: / / www.ecologica.cn
的低贴现率意味着更加重视长期的环境效益和后代
的利益。 因此政府决策时要尽可能采用低贴现率。
需要指出的是,由于数据的缺乏,填海造地的某
些成本,如填海造地对相邻海岸生态系统服务造成
的损害、填海造地导致环境容量减少造成的其他污
染物(如氮、磷等营养盐)处理成本增加等,在经验估
计中没有被包括。 而海洋营养盐污染是我国近岸海
域面临的严重环境问题之一。 所以估算出来的最大
可允许填海面积是能够填海的最大值。
虽然本研究是以厦门为案例,但是论文提出的
分析框架和模型可适用于其他区域,也适用于不同
规模填海造地方案的评估。 最后需注意的是,模型
估算出的最大可允许填海面积,或者最小必须保持
的海域面积,必须通过法律规定下来,使之成为填海
造地红线制度的一部分。 这样才能确保填海造地得
到有效控制。
References:
[ 1 ]摇 World Resources Institute ( WRI ) . Pilot Analysis of Globe
Ecosystems: Coastal Ecosystems. Washington: WRI
Publication, 2001.
[ 2 ] 摇 Pernetta J, Elder D. Cross鄄Sectoral, Integrated Coastal Area
Planning (CICPAP): Guidelines and Principle for Coastal Area
Development. Gland: IUCN, 1993.
[ 3 ] 摇 United Nations Atlas of the Oceans ( UNAOO) ( 2006) Land
reclamation by the coast. [ 2012鄄 02鄄 04 ] . http: / / www.
oceansatlas. org / servlet / CDSServlet? status =
ND0xODAwOCY2PWVuJjMzPSomMzc9a29z.
[ 4 ] 摇 Xia Z, Jia P H, Lei Y, Chen Y Z. Dynamics of coastal land use
patterns of inner Lingdingyang bay in the Zhujiang river estuary.
Chinese Geographical Science, 2007, 17(3): 222鄄228.
[ 5 ] 摇 Hoeksema R J. Three stages in the history of land reclamation in
the Netherlands. Irrigation and Drainage, 2007, 56 ( S1 ):
S113鄄S126.
[ 6 ] 摇 Airoldi L, Beck M W. Loss, status and trends for coastal marine
habitats of Europe. Oceanography and Marine Biology, 2007, 45:
345鄄405.
[ 7 ] 摇 Halpern B S, Walbridge S, Selkoe K A, Kappel C V, Micheli F,
D忆Agrosa C, Bruno J F, Casey K S, Ebert C, Fox H E, Fujita
R, Heinemann D, Lenihan H S, Madin E M P, Perry M T, Selig
E R, Spalding M, Steneck R, Watson R. A global map of human
impact on marine ecosystems. Science, 2008, 319 ( 5865 ):
948鄄952.
[ 8 ] 摇 Editorial M, Cutts N D. Marine habitats: loss and gain, mitigation
and compensation. Marine Pollution Bulletin, 2004, 49(9 / 10):
671鄄674.
[ 9 ] 摇 Li T H, Peng H, Zhao Z J. Impact analysis of coastal engineering
projects on mangrove wetland area change with remote sensing.
China Ocean Engineering, 2008, 22(2): 347鄄358.
[10] 摇 Zhang Q, Xu C Y, Chen Y D. Wavelet鄄based characterization of
water level behaviors in the Pearl River estuary, China. Stochastic
Environmental Research and Risk Assessment, 2010, 24 ( 1):
81鄄92.
[11] 摇 Fu Y B, Cao K, Wang F, Zhang F S. On the quantitative
evaluation method of intensity and potential of reclamation. Ocean
Development and Management, 2010, 27(1): 27鄄30.
[12] 摇 State Oceanic Administration (SOA). Bulletin of China忆s Marine
Environmental Quality 2002. Beijing, 2003.
[13] 摇 GEF / UNDP / IMO Regional Programmer on Partnerships in
Environmental Management for the Seas of East Asia (PEMSEA).
The Development of National Coastal and Marine Policies in the
People忆s Republic of China: A Case Study. Printed in Quezon
City, 2003.
[14] 摇 Peng B R, Chen W Q, Hong H S. Integrating ecological damages
into the user charge for land reclamation: a case study of Xiamen,
China. Stochastic Environmental Research and Risk Assessment,
2011, 25(3): 341鄄351.
[15] 摇 Task Force of the China Council for International Cooperation on
Environment and Development (CCICED). Ecosystem Issues and
Policy Options addressing the Sustainable Development of China忆s
Ocean and Coast. Beijing: China Environmental Science Press,
2012: 30,50.
[16] 摇 Hartwick J N, Olewiler N D. The Economics of Natural Resource
Use. 2nd ed. New York: Addison鄄Wesley, 1997.
[17] 摇 Zhang L P, Chen W, Fang Q. Retrospective Assessment of
Environmental Impact of Reclamation Planning in Bays of Fujian
Province. Beijing: Science Press, 2008:40.
[18] 摇 MPP鄄EAS (Regional Program for the Prevention and Management
of Marine Pollution in the East Asian Seas ) . Socioeconomic
Assessment Framework and Guidelines for Integrated Coastal
Management. Quezon: Philippines: MPP鄄EAS / info / 99 / 199.
GEF / UNDP / IMO, 1999.
[19] 摇 Environment Protection Agency (EPA). Guidelines for Preparing
Economic Analyses. EPA 240鄄R鄄00鄄003. 2000.
[20] 摇 4th Institute of Navigational Engineering Exploration and Design,
Ministry of Communication ( 4th INEED). Feasibility Research
Report of Xiangyu Free Trade Zone ( second phase ) Project.
Xiamen, 2002.
[21] 摇 Panning and Design Research Institute of Navigational Bureau of
Yangtze River ( PDRINBYR). Feasibility Research Report of
Removing Silt and Shoreline Treatment Project in Seas Around
Xiamen Bridge. Xiamen, (2003) .
[22] 摇 Peng B R, Hong H, Chen W Q, Xue X Z, CaoA X L, Peng J P.
Ecological damage appraisal of sea reclamation: Theory, method
and application. Journal of Natural Resources, 2005, 20 ( 5):
714鄄726.
[23] 摇 Li J M, Liu T Y. The key points and empirical analysis of the
reclamation忆 s exterior ecology cost compensation. Ecological
Economy, 2010, (3): 143鄄146.
[24] 摇 Peng B R, Hong H. Valuation of Coastal Ecosystem Services:
Theory and Application. Beijing: Ocean Press, 2006: 155鄄158.
[25] 摇 Coastal and Ocean Management Institute, Xiamen University
(COMI). The Compensation Standard of Ecological Damage in
277 摇 生摇 态摇 学摇 报摇 摇 摇 34卷摇
http: / / www.ecologica.cn
Sea Area Use and Implementation Scheme in Xiamen.
Xiamen, 2012.
[26] 摇 Fujian marine institute and the State Key Laboratory of Marine
Environmental Science ( Xiamen University) . Study on the total
quantity control of sea reclamation in West Sea and Tongan Bay of
Xiamen. Xiamen, 2006.
[27] 摇 Xiamen Environmental Protection Bureau ( EPB ). Xiamen
Environmental State Report 2006. Xiamen, 2007.
[28] 摇 Third Institute of Oceanography, State Oceanic Administration
(TIO). Study on the Environmental Carrying Capacity of Xiamen
Seas. Xiamen, 1995.
[29] 摇 Chen W Q, Zhang L P, Hong H S, Xue X Z. A preliminary study
on evaluation of coastal sea water environmental capacity. Journal
of Xiamen University: Natural Science, 1999, 38(6): 896鄄901.
参考文献:
[11]摇 付元宾, 曹可, 王飞, 张丰收. 围填海强度与潜力定量评价方
法初探. 海洋开发与管理, 2010, 27(1): 27鄄30.
[12] 摇 国家海洋局. 2002 年中国海洋环境质量公报. 北京, 2003.
(请核对本条文献格式)
[15]摇 中国海洋可持续发展的生态环境问题与政策研究课题组. 中
国海洋可持续发展的生态环境问题与政策研究. 北京: 中国
环境科学出版社, 2012: 30鄄50.
[17] 摇 张骆平,陈伟琪,方秦华. 福建省海湾围填海规划环境影响回
顾性评价. 北京: 科学出版社, 2008: 40鄄40.
[20] 摇 中交第四航务工程勘察设计院. 厦门象屿保税区二期工程可
行性研究报告. 厦门, 2002. (请核对本条文献格式)
[21] 摇 长江航道规划设计研究院. 厦门大桥东侧海域清淤及岸线整
治(嘉庚广场)工程可行性研究报告. 厦门, 2003. (未链接到
本条信息, 请核对)
[22] 摇 彭本荣, 洪华生, 陈伟琪, 薛雄志, 曹秀丽, 彭晋平. 填海造
地生态损害评估: 理论、方法及应用研究. 自然资源学报,
2005, 20(5): 714鄄726.
[23] 摇 李京梅, 刘铁鹰. 填海造地外部生态成本补偿的关键点及实
证分析. 生态经济, 2010, (3): 143鄄146.
[24] 摇 彭本荣, 洪华生. 海岸带生态系统服务价值评估理论与应用
研究. 北京: 海洋出版社, 2006: 155鄄158.
[25] 摇 厦门大学海洋与海岸带发展研究院(COMI). 厦门海洋生态损
害补偿 /赔偿标准及实施方案. 厦门, 2012. (未链接到本条信
息, 请核对)
[26] 摇 福建海洋研究所, 近海海洋环境国家重点实验室(厦门大
学) . 厦门西海域及同安湾海域面积(围填海)总量控制研究
报告. 厦门, 2006. (未链接到本条信息, 请核对)
[27] 摇 厦门市环境保护局. 2006 年厦门环境状况公报. 厦门, 2007.
(请核对本条文献格式)
[28] 摇 国家海洋局第三海洋研究所. 厦门海域环境容量研究. 厦门,
1995. (未链接到本条信息, 请核对)
[29] 摇 陈伟琪, 张珞平, 洪华生, 薛雄志. 近岸海域环境容量的价值
及其价值量评估初探. 厦门大学学报: 自然科学版, 1999, 38
(6): 896鄄901.
377摇 3期 摇 摇 摇 林琛琛摇 等:最大可允许填海面积模拟———厦门西海域案例研究 摇
粤悦栽粤 耘悦韵蕴韵郧陨悦粤 杂陨晕陨悦粤 灾燥造援猿源袁晕燥援猿 云藻遭援袁圆园员源渊杂藻皂蚤皂燥灶贼澡造赠冤
悦韵晕栽耘晕栽杂
云则燥灶贼蚤藻则泽 葬灶凿 悦燥皂责则藻澡藻灶泽蚤增藻 砸藻增蚤藻憎
砸藻泽藻葬则糟澡 葬凿增葬灶糟藻泽 葬灶凿 责则燥泽责藻糟贼泽 燥枣 责燥泽贼鄄凿蚤泽葬泽贼藻则 藻糟燥造燥早蚤糟葬造 则藻泽贼燥则葬贼蚤燥灶 葬泽泽藻泽泽皂藻灶贼
蕴陨哉 载蚤葬燥枣怎袁宰粤晕郧 宰藻灶躁蚤藻袁 蕴陨 允蚤灶早袁 藻贼 葬造 渊缘圆苑冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
栽澡藻 泽怎泽贼葬蚤灶葬遭蚤造蚤贼赠 葬灶葬造赠泽蚤泽 燥枣 泽皂葬造造 澡赠凿则燥责燥憎藻则 责造葬灶贼泽 蚤灶 悦澡蚤灶葬 遭葬泽藻凿 燥灶 藻糟燥造燥早蚤糟葬造 藻灶藻则早藻贼蚤糟 葬糟糟燥怎灶贼蚤灶早
孕粤晕郧 酝蚤灶早赠怎藻袁 在匀粤晕郧 蕴蚤曾蚤葬燥袁 宰粤晕郧 悦澡葬灶早遭燥 渊缘猿苑冤
噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
粤怎贼藻糟燥造燥早赠 驭 云怎灶凿葬皂藻灶贼葬造泽
粤灶 藻糟燥造燥早蚤糟葬造 泽贼怎凿赠 燥灶 扎燥燥责造葬灶噪贼燥灶 蚤灶 贼澡藻 灶燥则贼澡藻则灶 月藻蚤遭怎 郧怎造枣 灾院 贼澡藻 藻枣枣藻糟贼泽 燥枣 皂蚤糟则燥扎燥燥责造葬灶噪贼燥灶 早则葬扎蚤灶早 燥灶 责澡赠贼燥责造葬灶噪贼燥灶 蚤灶
泽怎皂皂藻则 酝粤 蕴怎袁悦粤韵 宰藻灶择蚤灶早袁在匀粤晕郧 宰藻灶躁蚤灶早袁 藻贼 葬造 渊缘源远冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
杂藻曾怎葬造 泽蚤扎藻 凿蚤皂燥则责澡蚤泽皂 葬灶凿 枣藻皂葬造藻 蚤灶凿蚤增蚤凿怎葬造 枣藻糟怎灶凿蚤贼赠 燥枣 杂蚤造怎则怎泽 葬泽燥贼怎泽 葬灶凿 悦造葬则蚤葬泽 枣怎泽糟怎泽
云粤晕 载蚤葬燥造蚤袁 蕴陨晕 在澡蚤澡怎葬袁 阅陨晕郧 载蚤葬灶造燥灶早袁 藻贼 葬造 渊缘缘缘冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
耘枣枣藻糟贼 燥枣 藻灶增蚤则燥灶皂藻灶贼葬造 贼藻皂责藻则葬贼怎则藻 燥灶 贼澡藻则皂燥早藻灶藻泽蚤泽 葬灶凿 藻增葬责燥则葬贼蚤增藻 憎葬贼藻则 造燥泽泽 蚤灶 悦澡蚤灶藻泽藻 遭怎造遭怎造泽 渊孕赠糟灶燥灶燥贼怎泽 泽蚤灶藻灶泽蚤泽冤
蕴陨晕 蕴蚤灶袁悦粤韵 酝藻灶早贼蚤灶早袁匀哉 再蚤造蚤灶袁 藻贼 葬造 渊缘远源冤
噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
耘枣枣藻糟贼泽 燥枣 蚤则则蚤早葬贼蚤燥灶 燥灶 责造葬灶贼 早则燥憎贼澡 葬灶凿 灶蚤贼则燥早藻灶 怎泽藻 糟澡葬则葬糟贼藻则蚤泽贼蚤糟泽 燥枣 悦葬造造蚤早燥灶怎皂 糟葬责怎贼鄄皂藻凿怎泽葬藻 杂糟澡则藻灶噪 泽藻藻凿造蚤灶早泽
匀哉粤晕郧 悦葬蚤遭蚤葬灶袁 在耘晕郧 云葬灶躁蚤葬灶早袁 蕴耘陨 允蚤葬择蚤葬灶早袁 藻贼 葬造 渊缘苑圆冤
噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
耘灶扎赠皂葬贼蚤糟 葬糟贼蚤增蚤贼赠 凿怎则蚤灶早 责澡赠贼燥则藻皂藻凿蚤葬贼蚤燥灶 燥枣 责燥造赠糟赠糟造蚤糟 葬则燥皂葬贼蚤糟 澡赠凿则燥糟葬则遭燥灶 蚤皂责葬糟贼藻凿 泽燥蚤造
在匀哉 云葬灶袁 匀韵晕郧 载蚤葬灶早择蚤袁 再粤晕 宰藻灶凿藻袁 藻贼 葬造 渊缘愿员冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
粤泽泽藻泽泽皂藻灶贼 燥枣 噪葬则泽贼 早则燥怎灶凿憎葬贼藻则 增怎造灶藻则葬遭蚤造蚤贼赠 蚤灶 悦澡燥灶早择蚤灶早 遭葬泽藻凿 燥灶 则藻增蚤泽藻凿 砸陨杂运耘 皂燥凿藻造
宰耘陨 载蚤灶早责蚤灶早袁孕哉 允怎灶遭蚤灶早袁在匀粤韵 悦澡怎灶赠燥灶早 渊缘愿怨冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
酝燥灶蚤贼燥则蚤灶早 泽葬造贼 葬灶凿 憎葬贼藻则 凿赠灶葬皂蚤糟泽 蚤灶 枣葬则皂造葬灶凿 葬灶凿 凿则葬蚤灶葬早藻 凿蚤贼糟澡 蚤灶 葬 泽葬造蚤灶藻 藻灶增蚤则燥灶皂藻灶贼 怎灶凿藻则 则藻凿怎糟藻凿 凿则葬蚤灶葬早藻 蚤灶贼藻灶泽蚤贼赠
孕粤晕 再葬灶曾蚤灶袁 蕴哉韵 宰葬灶袁 允陨粤 在澡燥灶早澡怎葬袁 藻贼 葬造 渊缘怨苑冤
噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
孕燥责怎造葬贼蚤燥灶袁 悦燥皂皂怎灶蚤贼赠 葬灶凿 耘糟燥泽赠泽贼藻皂
杂责葬贼蚤葬造 凿蚤泽贼则蚤遭怎贼蚤燥灶 葬灶凿 泽责葬贼蚤葬造 葬泽泽燥糟蚤葬贼蚤燥灶 燥枣 杂贼藻造造藻则葬 糟澡葬皂葬藻躁葬泽皂藻 责燥责怎造葬贼蚤燥灶 蚤灶 贼澡藻 凿蚤枣枣藻则藻灶贼 葬造贼蚤贼怎凿藻 蚤灶 蚤灶 凿藻早则葬凿藻凿 葬造责蚤灶藻 早则葬泽泽造葬灶凿
郧粤韵 云怎赠怎葬灶袁 在匀粤韵 悦澡藻灶早扎澡葬灶早袁 在匀哉韵 酝粤 蕴葬灶糟葬燥 渊远园缘冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
月蚤燥贼蚤糟 泽贼则藻泽泽 燥枣 责则藻凿葬贼燥则泽 燥灶 责燥责怎造葬贼蚤燥灶 燥枣 葬造造蚤早葬贼燥则 憎藻藻凿 枣造藻葬 遭藻藻贼造藻袁 粤早葬泽蚤糟造藻泽 澡赠早则燥责澡蚤造葬 渊 悦燥造援院 悦澡则赠泽燥皂藻造蚤凿葬藻冤
蕴陨哉 再怎枣葬灶早袁蕴陨 云藻蚤袁郧哉陨 云葬灶早赠葬灶袁 藻贼 葬造 渊远员猿冤
噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
悦澡葬则葬糟贼藻则蚤泽贼蚤糟泽 燥枣 责澡赠贼燥责造葬灶噪贼燥灶 糟燥皂皂怎灶蚤贼赠 蚤灶 贼澡藻 灶燥则贼澡藻则灶 杂燥怎贼澡 悦澡蚤灶葬 杂藻葬 蚤灶 泽怎皂皂藻则 葬灶凿 憎蚤灶贼藻则 酝粤 宰藻蚤袁 杂哉晕 允怎灶 渊远圆员冤噎噎噎
栽澡藻 葬皂燥怎灶贼 葬灶凿 凿赠灶葬皂蚤糟泽 燥枣 造蚤贼贼藻则枣葬造造 蚤灶 贼澡藻 灶葬贼怎则葬造 泽藻糟燥灶凿葬则赠 枣燥则藻泽贼 蚤灶 酝燥责葬灶 酝燥怎灶贼葬蚤灶
云粤晕 悦澡怎灶灶葬灶袁 郧哉韵 在澡燥灶早造蚤灶早袁 在匀耘晕郧 允蚤灶责蚤灶早袁 藻贼 葬造 渊远猿猿冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
耘枣枣藻糟贼泽 燥枣 糟燥灶贼蚤灶怎燥怎泽 凿则燥怎早澡贼 燥灶 泽燥蚤造 遭葬糟贼藻则蚤葬 责燥责怎造葬贼蚤燥灶泽 葬灶凿 糟燥皂皂怎灶蚤贼赠 凿蚤增藻则泽蚤贼赠 蚤灶 泽憎藻藻贼 糟澡藻则则赠 则澡蚤扎燥泽责澡藻则藻
蕴陨哉 云葬灶早糟澡怎灶袁 载陨晕郧 杂澡葬灶早躁怎灶袁 酝粤 匀葬蚤造蚤灶袁 藻贼 葬造 渊远源圆冤
噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
栽澡藻 遭葬泽蚤糟 责则蚤灶糟蚤责造藻 燥枣 则葬灶凿燥皂 枣燥则藻泽贼 葬灶凿 蚤贼泽 葬责责造蚤糟葬贼蚤燥灶泽 蚤灶 藻糟燥造燥早赠院 葬 糟葬泽藻 泽贼怎凿赠 燥枣 孕蚤灶怎泽 赠怎灶灶葬灶藻灶泽蚤泽
在匀粤晕郧 蕴藻蚤袁 宰葬灶早 蕴蚤灶造蚤灶袁 在匀粤晕郧 载怎凿燥灶早袁 藻贼 葬造 渊远缘园冤
噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
匝怎葬灶贼蚤贼葬贼蚤增藻 葬灶葬造赠泽蚤泽 燥枣 凿则蚤增蚤灶早 枣葬糟贼燥则泽 枣燥则 憎藻贼造葬灶凿 凿藻早则葬凿葬贼蚤燥灶 遭葬泽藻凿 燥灶 澡赠凿则燥造燥早赠 遭葬造葬灶糟藻
匀韵哉 孕藻灶早袁 杂匀耘晕 宰藻灶皂蚤灶早袁宰粤晕郧 匝蚤葬燥袁 藻贼 葬造 渊远远园冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
杂贼燥皂葬贼葬造 葬灶凿 藻灶增蚤则燥灶皂藻灶贼葬造 糟燥灶贼则燥造 燥灶 藻增葬责燥贼则葬灶泽责蚤则葬贼蚤燥灶 蚤灶 葬 责造葬灶贼葬贼蚤燥灶 蚤灶 贼澡藻 造燥憎藻则 皂燥怎灶贼葬蚤灶 葬则藻葬泽 燥枣 晕燥则贼澡 悦澡蚤灶葬
匀哉粤晕郧 匀怎蚤袁 酝耘晕郧 孕蚤灶早袁 在匀粤晕郧 允蚤灶泽燥灶早袁 藻贼 葬造 渊远远苑冤
噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
栽澡藻 枣造怎曾藻泽 燥枣 糟葬则遭燥灶 凿蚤燥曾蚤凿藻袁 皂藻贼澡葬灶藻 葬灶凿 灶蚤贼则燥怎泽 燥曾蚤凿藻 蚤灶 葬造责蚤灶藻 早则葬泽泽造葬灶凿 燥枣 贼澡藻 栽蚤葬灶泽澡葬灶 酝燥怎灶贼葬蚤灶泽袁 载蚤灶躁蚤葬灶早
匀耘 郧怎蚤曾蚤葬灶早袁 蕴陨 运葬蚤澡怎蚤袁 杂韵晕郧 宰藻蚤袁 藻贼 葬造 渊远苑源冤
噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
蕴葬灶凿泽糟葬责藻袁 砸藻早蚤燥灶葬造 葬灶凿 郧造燥遭葬造 耘糟燥造燥早赠
粤灶葬造赠泽蚤泽 贼澡藻 则藻造葬贼蚤燥灶泽澡蚤责 遭藻贼憎藻藻灶 藻糟燥造燥早蚤糟葬造 枣燥燥贼责则蚤灶贼 渊耘云冤 燥枣 灶蚤灶早曾蚤葬 葬灶凿 蚤灶枣造怎藻灶糟蚤灶早 枣葬糟贼燥则泽院 孕葬则贼蚤葬造 蕴藻葬泽贼鄄杂择怎葬则藻泽 砸藻早则藻泽泽蚤燥灶 渊孕蕴杂冤
酝粤 酝蚤灶早凿藻袁 酝粤 载怎藻躁怎葬灶袁 载陨耘 再蚤灶早扎澡燥灶早袁 藻贼 葬造 渊远愿圆冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
孕则燥枣蚤造藻 凿蚤泽贼则蚤遭怎贼蚤燥灶 责葬贼贼藻则灶泽 燥枣 泽燥蚤造 燥则早葬灶蚤糟 糟葬则遭燥灶 葬灶凿 贼澡藻 则葬贼藻 燥枣 糟葬则遭燥灶 泽藻择怎藻泽贼则葬贼蚤燥灶 蚤灶 晕蚤灶早曾蚤葬 陨则则蚤早葬贼蚤燥灶 在燥灶藻
阅韵晕郧 蕴蚤灶造蚤灶袁 再粤晕郧 匀葬燥袁 再哉 阅燥灶早泽澡藻灶早袁 藻贼 葬造 渊远怨园冤
噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
酝燥造藻糟怎造葬则 皂藻糟澡葬灶蚤泽皂 燥灶 糟葬则遭燥灶 凿蚤燥曾蚤凿藻 葬泽泽蚤皂蚤造葬贼蚤燥灶 燥枣 葬怎贼燥贼则燥责澡蚤糟 皂蚤糟则燥燥则早葬灶蚤泽皂 葬灶凿 糟葬则遭燥灶 贼则葬灶泽造燥糟葬贼蚤燥灶 蚤灶 葬早则蚤糟怎造贼怎则葬造 泽燥蚤造泽
宰哉 载蚤葬燥澡燥灶早袁 允陨粤晕 再葬灶袁 悦匀耘晕 载蚤葬燥躁怎葬灶袁 藻贼 葬造 渊苑园员冤
噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
砸藻泽燥怎则糟藻 葬灶凿 陨灶凿怎泽贼则蚤葬造 耘糟燥造燥早赠
耘糟燥造燥早蚤糟葬造 藻枣枣蚤糟蚤藻灶糟赠 葬灶葬造赠泽蚤泽 燥枣 贼澡藻 糟蚤则糟怎造葬则 藻糟燥灶燥皂赠 泽赠泽贼藻皂 蚤灶 皂蚤灶蚤灶早 葬则藻葬 遭葬泽藻凿 燥灶 藻皂藻则早赠 葬灶葬造赠贼蚤糟 葬责责则燥葬糟澡
杂哉晕 再怎枣藻灶早袁郧哉韵 匝怎葬灶赠蚤灶早 渊苑员园冤
噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
粤泽泽藻泽泽蚤灶早 泽赠灶贼澡藻贼蚤糟 糟葬则则赠蚤灶早 糟葬责葬糟蚤贼赠 遭葬泽藻凿 燥灶 粤阅鄄粤杂 皂燥凿藻造院 葬 糟葬泽藻 泽贼怎凿赠 蚤灶 悦燥葬泽贼葬造 在燥灶藻袁 在澡燥怎泽澡葬灶
杂哉 孕葬灶责葬灶袁 再耘 杂澡怎枣藻灶早袁 郧哉韵 在澡燥灶早赠葬灶早袁 藻贼 葬造 渊苑员愿冤
噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
哉则遭葬灶袁 砸怎则葬造 葬灶凿 杂燥糟蚤葬造 耘糟燥造燥早赠
阅则蚤增蚤灶早 枣燥则糟藻泽 葬灶葬造赠泽蚤泽 燥枣 怎则遭葬灶 藻曾责葬灶泽蚤燥灶 遭葬泽藻凿 燥灶 遭燥燥泽贼藻凿 则藻早则藻泽泽蚤燥灶 贼则藻藻泽 葬灶凿 蕴燥早蚤泽贼蚤糟 则藻早则藻泽泽蚤燥灶
蕴陨 悦澡怎灶造蚤灶袁 蕴陨哉 酝蚤葬燥袁 匀哉 再怎葬灶皂葬灶袁 藻贼 葬造 渊苑圆苑冤
噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
砸藻泽藻葬则糟澡 晕燥贼藻泽
杂责葬贼蚤葬造鄄贼藻皂责燥则葬造 增葬则蚤葬遭蚤造蚤贼赠 燥枣 凿则赠 葬灶凿 憎藻贼 凿藻责燥泽蚤贼蚤燥灶 燥枣 葬贼皂燥泽责澡藻则蚤糟 灶蚤贼则燥早藻灶 蚤灶 凿蚤枣枣藻则藻灶贼 藻糟燥造燥早蚤糟葬造 则藻早蚤燥灶泽 燥枣 杂澡葬葬灶曾蚤
蕴陨粤晕郧 栽蚤灶早袁 栽韵晕郧 再葬灶忆葬灶袁 蕴陨晕 宰藻灶袁 藻贼 葬造 渊苑猿愿冤
噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
栽澡藻 藻枣枣藻糟贼泽 燥枣 凿蚤枣枣藻则藻灶贼 皂怎造糟澡蚤灶早 憎葬赠 燥灶 泽燥蚤造 憎葬贼藻则 贼澡藻则皂葬造 糟澡葬则葬糟贼藻则蚤泽贼蚤糟泽 蚤灶 责藻葬则 燥则糟澡葬则凿 蚤灶 贼澡藻 葬则蚤凿 葬则藻葬
蕴陨哉 载蚤葬燥赠燥灶早袁 蕴陨 匀燥灶早曾怎灶袁蕴陨 允蚤葬灶皂蚤灶早袁 藻贼 葬造 渊苑源远冤
噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
杂贼则怎糟贼怎则藻 葬灶凿 凿蚤增藻则泽蚤贼赠 燥枣 泽燥蚤造 枣葬怎灶葬 糟燥皂皂怎灶蚤贼蚤藻泽 蚤灶 贼澡藻 贼怎灶凿则葬 燥枣 贼澡藻 悦澡葬灶早遭葬蚤 酝燥怎灶贼葬蚤灶泽袁 悦澡蚤灶葬
宰粤晕郧 在澡藻灶澡葬蚤袁 再陨晕 载蚤怎择蚤灶袁 允陨粤晕郧 再怎灶枣藻灶早 渊苑缘缘冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
酝燥凿藻造蚤灶早 贼澡藻 贼燥贼葬造 葬造造燥憎葬遭造藻 葬则藻葬 枣燥则 糟燥葬泽贼葬造 则藻糟造葬皂葬贼蚤燥灶院葬 糟葬泽藻 泽贼怎凿赠 燥枣 载蚤葬皂藻灶袁 悦澡蚤灶葬
蕴陨晕 悦澡藻灶糟澡藻灶袁 砸粤韵 匀怎葬灶澡怎葬灶袁 蕴陨哉 再葬灶袁 藻贼 葬造 渊苑远远冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
愿苑苑 摇 生摇 态摇 学摇 报摇 摇 摇 猿源卷摇
叶生态学报曳圆园员源年征订启事
叶生态学报曳是由中国科学技术协会主管袁中国生态学学会尧中国科学院生态环境研究中心主办的生态学
高级专业学术期刊袁创刊于 员怨愿员年袁报道生态学领域前沿理论和原始创新性研究成果遥 坚持野百花齐放袁百家
争鸣冶的方针袁依靠和团结广大生态学科研工作者袁探索生态学奥秘袁为生态学基础理论研究搭建交流平台袁
促进生态学研究深入发展袁为我国培养和造就生态学科研人才和知识创新服务尧为国民经济建设和发展服务遥
叶生态学报曳主要报道生态学及各分支学科的重要基础理论和应用研究的原始创新性科研成果遥 特别欢
迎能反映现代生态学发展方向的优秀综述性文章曰研究简报曰生态学新理论尧新方法尧新技术介绍曰新书评价和
学术尧科研动态及开放实验室介绍等遥
叶生态学报曳为半月刊袁大 员远开本袁圆愿园页袁国内定价 怨园元 辕册袁全年定价 圆员远园元遥
国内邮发代号院愿圆鄄苑袁国外邮发代号院酝远苑园
标准刊号院陨杂杂晕 员园园园鄄园怨猿猿摇 摇 悦晕 员员鄄圆园猿员 辕 匝
全国各地邮局均可订阅袁也可直接与编辑部联系购买遥 欢迎广大科技工作者尧科研单位尧高等院校尧图书
馆等订阅遥
通讯地址院 员园园园愿缘 北京海淀区双清路 员愿号摇 电摇 摇 话院 渊园员园冤远圆怨源员园怨怨曰 远圆愿源猿猿远圆
耘鄄皂葬蚤造院 泽澡藻灶早贼葬蚤曾怎藻遭葬燥岳 则糟藻藻泽援葬糟援糟灶摇 网摇 摇 址院 憎憎憎援藻糟燥造燥早蚤糟葬援糟灶
编辑部主任摇 孔红梅摇 摇 摇 执行编辑摇 刘天星摇 段摇 靖
生摇 态摇 学摇 报渊杂匀耘晕郧栽粤陨摇 载哉耘月粤韵冤渊半月刊摇 员怨愿员年 猿月创刊冤
第 猿源卷摇 第 猿期摇 渊圆园员源年 圆月冤
粤悦栽粤 耘悦韵蕴韵郧陨悦粤 杂陨晕陨悦粤摇渊杂藻皂蚤皂燥灶贼澡造赠袁杂贼葬则贼藻凿 蚤灶 员怨愿员冤摇灾燥造郾 猿源摇 晕燥郾 猿 渊云藻遭则怎葬则赠袁 圆园员源冤
编摇 摇 辑摇 叶生态学报曳编辑部
地址院北京海淀区双清路 员愿号
邮政编码院员园园园愿缘
电话院渊园员园冤远圆怨源员园怨怨憎憎憎援藻糟燥造燥早蚤糟葬援糟灶泽澡藻灶早贼葬蚤曾怎藻遭葬燥岳 则糟藻藻泽援葬糟援糟灶
主摇 摇 编摇 王如松
主摇 摇 管摇 中国科学技术协会
主摇 摇 办摇 中国生态学学会
中国科学院生态环境研究中心
地址院北京海淀区双清路 员愿号
邮政编码院员园园园愿缘
出摇 摇 版摇
摇 摇 摇 摇 摇 地址院北京东黄城根北街 员远号
邮政编码院员园园苑员苑
印摇 摇 刷摇 北京北林印刷厂
发 行摇
地址院东黄城根北街 员远号
邮政编码院员园园苑员苑
电话院渊园员园冤远源园猿源缘远猿耘鄄皂葬蚤造院躁燥怎则灶葬造岳 糟泽责早援灶藻贼
订摇 摇 购摇 全国各地邮局
国外发行摇 中国国际图书贸易总公司
地址院北京 猿怨怨信箱
邮政编码院员园园园源源
广告经营
许 可 证摇 京海工商广字第 愿园员猿号
耘凿蚤贼藻凿 遭赠摇 耘凿蚤贼燥则蚤葬造 遭燥葬则凿 燥枣
粤悦栽粤 耘悦韵蕴韵郧陨悦粤 杂陨晕陨悦粤
粤凿凿院员愿袁杂澡怎葬灶早择蚤灶早 杂贼则藻藻贼袁匀葬蚤凿蚤葬灶袁月藻蚤躁蚤灶早 员园园园愿缘袁悦澡蚤灶葬
栽藻造院渊园员园冤远圆怨源员园怨怨
憎憎憎援藻糟燥造燥早蚤糟葬援糟灶
泽澡藻灶早贼葬蚤曾怎藻遭葬燥岳 则糟藻藻泽援葬糟援糟灶
耘凿蚤贼燥则鄄蚤灶鄄糟澡蚤藻枣摇 宰粤晕郧 砸怎泽燥灶早
杂怎责藻则增蚤泽藻凿 遭赠摇 悦澡蚤灶葬 粤泽泽燥糟蚤葬贼蚤燥灶 枣燥则 杂糟蚤藻灶糟藻 葬灶凿 栽藻糟澡灶燥造燥早赠
杂责燥灶泽燥则藻凿 遭赠摇 耘糟燥造燥早蚤糟葬造 杂燥糟蚤藻贼赠 燥枣 悦澡蚤灶葬
砸藻泽藻葬则糟澡 悦藻灶贼藻则 枣燥则 耘糟燥鄄藻灶增蚤则燥灶皂藻灶贼葬造 杂糟蚤藻灶糟藻泽袁 悦粤杂
粤凿凿院员愿袁杂澡怎葬灶早择蚤灶早 杂贼则藻藻贼袁匀葬蚤凿蚤葬灶袁月藻蚤躁蚤灶早 员园园园愿缘袁悦澡蚤灶葬
孕怎遭造蚤泽澡藻凿 遭赠摇 杂糟蚤藻灶糟藻 孕则藻泽泽
粤凿凿院员远 阅燥灶早澡怎葬灶早糟澡藻灶早早藻灶 晕燥则贼澡 杂贼则藻藻贼袁
月藻蚤躁蚤灶早摇 员园园苑员苑袁悦澡蚤灶葬
孕则蚤灶贼藻凿 遭赠摇 月藻蚤躁蚤灶早 月藻蚤 蕴蚤灶 孕则蚤灶贼蚤灶早 匀燥怎泽藻袁
月藻蚤躁蚤灶早 员园园园愿猿袁悦澡蚤灶葬
阅蚤泽贼则蚤遭怎贼藻凿 遭赠摇 杂糟蚤藻灶糟藻 孕则藻泽泽
粤凿凿院员远 阅燥灶早澡怎葬灶早糟澡藻灶早早藻灶 晕燥则贼澡
杂贼则藻藻贼袁月藻蚤躁蚤灶早 员园园苑员苑袁悦澡蚤灶葬
栽藻造院渊园员园冤远源园猿源缘远猿
耘鄄皂葬蚤造院躁燥怎则灶葬造岳 糟泽责早援灶藻贼
阅燥皂藻泽贼蚤糟 摇 摇 粤造造 蕴燥糟葬造 孕燥泽贼 韵枣枣蚤糟藻泽 蚤灶 悦澡蚤灶葬
云燥则藻蚤早灶 摇 摇 悦澡蚤灶葬 陨灶贼藻则灶葬贼蚤燥灶葬造 月燥燥噪 栽则葬凿蚤灶早
悦燥则责燥则葬贼蚤燥灶
粤凿凿院孕援韵援月燥曾 猿怨怨 月藻蚤躁蚤灶早 员园园园源源袁悦澡蚤灶葬
摇 陨杂杂晕 员园园园鄄园怨猿猿悦晕 员员鄄圆园猿员 辕 匝 国内外公开发行 国内邮发代号 愿圆鄄苑 国外发行代号 酝远苑园 定价 怨园郾 园园元摇