茜坑水库水质生态修复的数值模拟.docVIP

浅水水体水质生态修复的数值模拟 彭 虹1 王艳2 张万顺2 陈文祥3 冯桃辉2 （1武汉大学水利水电学院，武汉 430070；2武汉大学资源与环境科学学院，武汉 430079；3水利部中国科学院水工程生态研究所，武汉 430079） 摘 要： 本文基于二维水流水动力学模型、具有源汇项的对流扩散方程及水生态系统动力学模型，考虑浅水水体生物群落中影响水体水质及营养物质循环的主要种群的演化规律，建立浅水水体生态修复模型。该模型包括水动力学模型、守恒物质对流扩散模型和富营养化动力学模型。其中富营养化动力学模型考虑了浮游植物、浮游动物、蚌类、溶解氧和各种营养物质。将该模型应用于深圳市某实体水库，模拟人工养殖蚌类修复水库水质的效果。数值模拟结果表明：与蚌类养殖前相比，浮游植物的生物量以及氮、磷的负荷有明显的减少，说明将蚌类用于富营养化水体水质的修复效果是比较理想的，应当加以推广。 关键词：水环境生态学；生态修复；河蚌；富营养化 一 引言 近年来,由于国民经济迅猛发展,人口剧增,大量工业废水及生活污水排入江河湖泊中,使浅水水体中的营养负荷不断增加，富营养化已经成为影响湖泊水质的主要问题，目前对于东湖、滇池、太湖等我国的大型浅水湖泊，虽然采取各种措施来控制其富营养化，但无论是外源控制消减营养负荷，还是采用各种除藻措施（如养鱼等），均未能取得预期的效果或只在短期内有效，甚至有时因采用不适当的物理、化学方法反而加剧了原本脆弱的水生态系统的破坏。这说明富营养化产生机制是很复杂的，要对其进行控制，除了采用物理化学以及生态的方法外，更重要的是对其进行机理性的研究，数学模型为该问题的研究提供了强有力的工具。 早期的模型是单一的N、P负荷模型[1]。70年代世界各地开展了大规模的湖泊富营养化调查，在分析资料的基础上建立了一系列藻类生物量与营养物质负荷量之间的相关经验模型[2]，但这类模型不能全面考虑浮游植物的动态变化。研究生态系统内子系统间相互作用过程,综合考虑系统外部环境驱动变量的生态动力学模型模拟水体的富营养化过程是采用数学模型方法模拟系统变化，具有较强的物理概念，对系统的发展趋势有一定的预见性，因而是目前众多学者采用的方法之一。生态动力学模型的研究始于Chen, Di toro开发的简单水质动力模型[3]4].主要有GlusmΦ生态模型、Clean系列湖泊综合模式、WASP模型等。 近年来，随着对水生生态系统构成的认识不断深入，通过改变食物网结构达到控制富营养化水体藻类数量成为新的研究热点。并提出了生态系统的修复方法－生物操纵。目前国内对此生态修复方法的模拟尚不多见。本研究在实验数据的基础上，模拟了人工养殖蚌类对水库环境中藻类生物量的影响。 本文基于二维水流动力学模型，采用具有源汇项的对流扩散方程及水生态系统动力学模型，考虑浅水水体中生物群落中影响水体水质及营养物质循环的主要种群的演化规律，建立了浅水水体生态动力学模型，并将该模型应用于深圳市某实体水库，在人工养殖蚌类实验的基础上模拟了人工养殖蚌类对水库水生态环境的影响。 二 数值模型理论 2.1水动力学模型 采用单宽流量、及水位作为因变量的圣维南方程组，以及营养物质的对流扩散方程。 连续方程： (1) 动量方程 (2) (3) 2.2守恒物质对流扩散模型 考虑稀释扩散过程、迁移传递过程，以及物质转化及生物的新陈代谢过程，污染物的直接负荷、边界负荷等因素，水质成分满足如下的守恒方程： (4) 2.3 富营养化的动力学模型 本模型主要考虑11种物质成分,包括:浮游植物(Ｃ4)、浮游动物(Ｃ9)、贝类（C11）、五日生化需氧量BOD5 (Ｃ5)、有机氮(Ｃ7)、氨氮(Ｃ1)、硝氮(Ｃ2)、有机磷(Ｃ8)、无机磷(Ｃ3)、溶解氧(Ｃ6)。水体中各成分之间相互作用有如图1所示[5]。 图1 主要水质成分及其相互关系 （1）水生生物动力学系统 根据图1所示各水质成分的相互关系，结合各种生物的生长规律，给出浮游植物、浮游动物以及蚌类的生态动力学模型： 浮游植物： （5） 浮游动物: （6） 双壳蚌类： （7） （2）磷循环 包括有机磷和无机磷。 无机磷： （8） 有机磷： （9） （3）氮循环 包括有机氮、氨氮和硝氮。 有机氮： (10) 氨氮： (11) 硝氮：