清华大学数学建模竞赛.pdfVIP

福田红树林自然保护区湿地生态系统 模型框架的构建及应用实例研究 摘要： 结合具体生态保护区的实例进行生态模拟分析，通过对种群、能量、生态因子三个方面 的建模进行数学模拟，讨论需要收集的数据指标，并通过这些指标进行生态区健康程度 的定量检测。此外，也针对外部因素对模型的影响和稳态建立条件进行了简单分析。 第一部分 能量流动模型 深圳福田红树林自然保护区位于深圳湾北岸，区域范围在北纬 22°30’- 22°32’， 东经 113°56’ - 114°3’之间，总面积 3.68 平方公里，是我国面积最小的红树林保 护区，也是我国唯一处于城市腹地的国家级自然保护区。由于这个特殊的地理位置，福 田红树林生态系统对于深圳一带的气候、湿度、生态调节起到了至关重要的作用，红树 林湿地生态系统由红树植物、其它陆地植被、鸟类、两栖爬行动物、昆虫、底栖动物、 浮游生物等生物以及其赖以生存的土壤、大气、海水等环境要素共同组成，通过物质循 环和能量转换发挥生态功能。但是今年来随着深圳经济的告诉发展，这种以环境为代价 的发展模式渐渐出现了问题，整个生态系统在空间结构与生态过程、自我调节与更新能 力以及对外部胁迫的恢复能力等方面，均呈现出不稳定性和不可持续性征兆。 下面就能量方面对于福田红树林的生态系统进行评价，以研究该生态系统的健康状 况。 模型建立： 1 能量的流动与物质传递时相辅相成的，物质是能量的载体，由于能量的输入输出是 守恒的，遵循能量守恒定律，为了研究红树林生态系统的能量利用情况，这里采用用 C 元素的积累来间接能量的积累，Q = F(C),(Q 表示能量，C 代表 c 元素)，初步可以这样 认为生产者固定的碳元素（A）等于呼吸排出的二氧化碳（S）、沉积在土壤的有机物（W）、 消费者捕食量（E），即可认为该生态系统处于稳定状态。等式： A=S+W+E. 现在一项项分析，首先是生产者固定的 C ，可以采取样方法，在一定面积大小的区 域，测定在一定时间内的光合作用强度，而光合作用和太阳辐射 j 及叶绿素 t 含量等有 关，故 A=f （j，t）在计算全年的光合作用产量时刻取太阳辐射w 欸全年平均水平，而 叶绿素含量，由于地处北纬 22°30’ - 22°32’，福田红树林基本属于常绿阔叶林， 全年的叶绿素含量变化幅度很小，只需要样方法测出单位面积叶绿素含量即可。 S= ∗ （面积） 呼吸作用，叶片呼吸量满足该公式： R1 = 0.083 2e0.082 （R 2= 0.82） 式中， R1 为叶片单位面积的呼吸速率（μmol ·m-2 ·s-1）； T 为气温（℃）。进而通过叶面积指数求出单位群落面积的叶片呼吸量，公式如下： = R ⋅ B / LMA 1 1 式中， RL 为单位群落面积叶片呼吸量（μmol ·m-2 ·s-1）；R1 为叶片单位面积的 呼吸速率（μmol ·m-2 ·s-1）；LMA 为叶片单位面积质量（m2 ·t-1），通过对秋茄群 落的测定，LMA 为 1.094 9 m2 ·t-1；B1 为单位群落面积的叶生物量.。 沉淀在土壤的有机物主要有三部分组成：枯立木，枯倒木，凋落物， 枯立木采用软尺和树高测定仪分别测量枯立木的胸径（1.3 m 处）和树高，利用修 改的活立木生物量异速生长方程计算生物量。 ( 2 ) = a (1 − ) 式中，Ba 表示枯立木的生物量（t/hm2）；D、H、a、b 等系数含义见公式； ω 表示枯立木缺枝少叶的折算系数，对于上述Ⅰ类枯立木，ω=0.02~0.03；对于Ⅱ 类，ω=0.2；对于Ⅲ类和Ⅳ类枯立木，ω 值分别设定为 0.3 和 0.5。 枯倒木：