1. 引言
最近几十年里,我国众多水库正处于富营养化过程中,水库富营养化问题已对人类生存环境及饮用水安全构成了严重威胁,并制约了社会、经济及环境的可持续发展。随着现代经济的高速发展和人口数量剧增导致了水库水体中的营养盐不断增加,超过了水库本身的自净能力,造成水库藻类大量繁殖,引起藻类水华频繁爆发,形成藻类水华现象。因此,随着世界经济的快速发展,水体富营养化引起的藻类水华问题已成为一个日趋严重的全球性问题之一 [1] [2] [3] 。
藻类水华形成过程是在一定的营养、气候、水文条件和生态环境下藻类种群过度繁殖和聚集的过程,是水体生态环境因子如TN、TP、温度、光照、pH、流速、风力、溶解氧等综合动态过程 [4] 。在最近几十年里,国内外学者分析和明确了影响藻类动态生长的关键因子,例如光照、温度、营养盐、pH值和流速等。论文 [5] 选择光照强度、温度和营养物作为限制因素研究了藻类生长规律,建立了镜泊湖藻类生长动态新模型,探析了此类因子与藻类种群动态生长的定量关系。论文 [6] 研究了光照强度、温度和营养盐对滇池铜绿微囊藻生长的影响规律,研究结果表明滇池中磷是限制藻类动态生长的主要因子,控制滇池富营养化应以控制磷为主。论文 [7] 主要研究了温度、光照强度和硝酸盐对藻类种群生长的影响规律,研究发现拟柱孢藻对温度、光照强度和氮源均有较宽的生态位,有利于在较大空间尺度上进行扩散。论文 [8] 研究了温度和光照对单细胞绿藻1102的影响特性,结果显示温度和光照对单细胞绿藻生长具有交互作用。论文 [9] 研究了温度、光照和磷酸盐脉冲输入对藻类种群生长的交互影响机制,证实了温度、光照和脉冲对藻类种群生长存在交互影响机制。论文 [10] 揭示了营养盐磷含量及不同温度对铜绿微囊藻生长的动态影响机制,研究结果较好地说明洋河水库铜绿微囊藻种群从早春到晚秋持续发生的原因。论文 [11] 研究了光照和磷的交互作用对两种淡水藻类生长的影响,研究结果显示光照和磷的交互作用对铜绿微囊藻的生长影响明显。论文 [12] 研究了光照和温度对藻类生长的影响规律,研究结果从藻类种群繁殖能力角度探索了光照强度和光照周期及其温度对藻类生长的影响特性及规律。论文 [13] 研究了光和营养盐如何影响浮游植物的大面积生长繁殖,研究结果揭示了藻类种群生长比率的大小严重依赖于光和营养盐的供给情况。论文 [14] 研究了氮,磷,铁和硅等因素对5种海洋底栖硅藻生长动态的影响情况,研究揭示了氮,磷,铁和硅等因素对硅藻生长比率的影响特性,5种硅藻可以再2~6天内到达最大生长比率。论文 [15] 主要研究了光照和营养盐在不同情况下对亚热带水域浮游植物生长动态的影响机制,结果显示不同程度的光照强度和不同比率的营养盐对浮游植物的生长动态演化有着不同程度的影响特性,且它们之间存在一定的协同交互影响机制。Lv等人 [16] 研究了氮磷、温度等因素对亚热带浅水湖泊浮游植物空间组成和生物量的影响规律,结果揭示温度和总氮是影响浮游植物生长的关键因素。综上所述,藻类种群生长动态过程是多种因素动态协同交互作用过程,影响其过程的主要因素是:1) 水体中丰富的N、P营养物质,BOD及其他微量元素;2) 适宜的温度和光照;3) 适宜的水文地理条件如缓慢的水流等。
2. 动态建模
吴家园水库位于温州苍南县藻溪镇吴家园村,是苍南县第二大饮用水水库,每年平均为60万人民群众提供清洁水源3300万方,但是 2014 年 9 月 9 日 爆发大规模藻类水华,一眼望去是绿油油的一片,被迫停止供水。此外吴家园水库在2012年10月和2014年8月也发生过不同程度的藻类水华现象。显而易见,温州地区饮用水水库富营养化问题还是比较严峻得,严重制约了当地经济的发展与周边城镇人口饮水安全。因此,水体富营养化和藻类水华预警预测及控制已成为当前迫切解决的问题之一。
众所周知,底泥磷的释放在维持藻类种群长期生长及其藻类水华爆发过程中具有举足轻重的作用。一般来说,当总磷浓度高于一个临界值时,底泥将吸收磷;当总磷浓度低于这个临界值时,底泥将释放磷。为了数学刻画这类动态机制,可以用数学函数来表示底泥的释放与吸收机制,其中都是常数,,,且表示吴家园水库常年监测到的藻类种群最大生物量。同时,函数具有一个重要特征:如果,则有具有正的递减性,如果,则有具有负的递减性。这个性质导致具有一个特性就是当一类藻类种群数量达到一定范围内,这类藻类种群将持续增加并最终成为主导藻类种群。
同时,时滞动力系统能够更好地描述水体富营养化过程中的延迟现象,特别是水体富营养化过程中的生态系统内部结构、驱动因素和交互过程中的时间延迟性。根据水域生态学所知,藻类种群吸收营养盐并转化为自身生长的营养物质需要一定时间段,进而会有一定的时间延迟现象。为了在动态建模过程中显示时间延迟的影响作用,在动态建模过程中将引入时滞效应,因此所建总磷限制下的藻类种群生长动态模型如下:
(1)
其中分别代表在时间刻总磷的浓度(mg/L)和藻类种群的生物量(ug/L),代表稀释率,代表外界输入水库的磷的浓度,代表磷的沉积速率,代表藻类种群来自其它营养源的平均增长率,代表藻类种群对磷的吸收率,表示半饱和常数,表示磷营养转化为藻类种群生长能量的转化率,代表吴家园水库水域营养提供的藻类种群最大环境容纳量,代表吴家园水库资源的最大环境容纳量,
代表藻类种群平均死亡率。同时表示区间上的有界函数,可以用来表示在时间
之间的藻类种群生长动态扩散率。现在假设时滞效应发生在种群动态交互作用过程中,也就是说藻类种群不可能直接吸收总磷营养盐后,就可以转化为维持自身生长的能量,而是需要一段时间。此外,模型(1)的初始条件满足:
,
其中。
在生态建模过程中,我们引入了分布时滞项来刻画水域生态系统过程中的一类时间延迟现象,并假设藻类种群生长严重依赖营养盐-总磷的限制,进而研究如何通过控制总磷的浓度来预防与控制藻类种群快速繁殖。
3. 理论分析
为了合理地推导出总磷营养盐与藻类种群生长动态的互动作用机制,将对所建动态模型(1)进行相关稳定性分析。藻类水华控制的实质问题就是尽可能地保持水库水域藻类种群的浓度偏低,尽可能等于零更好,但是这是不可能实现。然而借助生态数学模型可以理论上获得一些参数限制条件,保证藻类种群是灭绝的,也就是生态模型(1)的藻类种群不存在的边界平衡点是稳定的。如果生态模型(1)的所有解在边界平衡点处是趋于稳定,那么藻类种群即将面临灭绝。若藻类种群最终走向灭绝,则,得到:
1):若,则生态模型(1)存在一个边界平衡点,其中
.
2):若,则生态模型(1)存在两个平衡点和,其中
和
在范围内,下面将推导出这些边界平衡点稳定的关键阈值条件,这些理论条件是进一步推导出藻类种群走向灭绝的关键判断准则。
定理2.1. 平衡点是一个局部渐近稳定的当且仅当和。
证明:为了研究平衡点的稳定性,构建一类Lyapunov-Krasovskii函数,表示如下:
沿着生态模型(1)的轨迹对求导得到
因此,可得
如果,则
其中,
因此,当满足条件和时,平衡点是渐近稳定的。
基于上述相似的分析过程,可以得到下面两个定理,因为证明过程,所以不再给出证明过程。
定理3.2 平衡点是局部渐近稳定的当且仅当和。
定理3.3. 平衡点是局部渐近稳定的当且仅当和。
基础上述理论分析,能够得到维持边界平衡点稳定的一些关键参数的阈值表达式,而这些表达式能为下一步数值模拟提供一定的理论基础。
4. 结论
基于亚热带水库水域富营养化情况与吴家园水库水质监测情况,在动态建模过程中引入分布时滞项,建立了一类总磷限制下的藻类种群生长动态模型,对所建生态模型进行了一定的理论分析,研究了该生态模型边界平衡点的渐近稳定性,明确了该生态模型具有这些特定动力学性态的阈值条件,这些研究工作为进一步预测吴家园水库营养盐动态演化趋势和揭示水域优势藻类种群生长动态规律提供了一定的理论基础。
致谢
感谢浙江省水环境与海洋资源重点实验室的同仁们,感谢你们一直以来在实验室对我们问题的细致引导、鼓励和解答,使我们对专业有了更深更细致的理解。
感谢国家自然科学基金面上项目(31570364),浙江省公益技术研究项目基金项目(2015C33227),浙江省自然科学基金项目(LY16BO70008),温州科技计划项目(S20140028和S20140024)的资助。
参考文献