太湖水质富营养化特征研究(.docVIP

太湖水质富营养化特征研究? ? 基金项目: 国家自然科学基金资助项目（No. 作者简介: 许旭峰(1984－)，男，江苏省人，硕士研究生. 许旭峰，刘青泉 (中国科学院力学研究所环境力学实验室，北京 100190) 摘 要：本文以太湖富营养化问题为对象，选择叶绿素a、总磷和总氮三个水质变量，建立以藻类生长为核心的生态动力学模型，结合二维水动力模型，建立了太湖生态水动力学模型，并运用该模型对太湖藻类生长的特征和影响因子进行数值模拟研究。研究结果表明温度和总磷是太湖藻类增长的主要控制因子，给出了藻类增长及水华爆发的环境因子阈值。进一步通过对太湖不同湖区的藻类生长状况的模拟分析，发现太湖梅梁湾和西部沿岸区是太湖藻类生长和水华暴发的危险区域。 关键词：太湖，富营养化，水华，生态动力学，二维水动力模型 1 引言 太湖是我国第三大淡水湖泊，位于经济发达的长江三角洲，是当地经济发展和人民生活的重要淡水资源。太湖是典型的大型浅水湖泊，湖泊面积约2338，平均水深只有1.9，湖水滞留期约300天，各湖区水动力差异显著。由于周边工农业的迅速发展，太湖污染日趋严重，造成了湖泊富营养化，水质恶化，蓝藻水华频繁暴发。 为了治理太湖富营养化问题，需要对太湖富营养化特征进行深入的研究，其中数值模拟是重要的研究方法之一。近二十年来，在太湖富营养化模型选择上，主要有两种思路。一种是以水动力学方程和对流扩散方程为基础建立模型，重点描述营养盐的时空分布。这种模型没有考虑生态过程，参数较少。另一种是生态动力学模型，考虑了藻类生长以及其它生态过程，而不考虑水流引起的空间变化。目前，越来越多的研究将这两种模型结合起来，建立生态水动力学耦合模型，这种模型参数众多，计算量大。 虽然对太湖做了较多数值模拟研究，模型越来越完善，但是往往考虑了较多的生态过程和水质变量，模型非常复杂，真正适合太湖实际情况的模型还较少。本文以二维浅水模型和输移扩散模型为基础，加入以藻类生长为核心的生态动力学模型，只考虑较少的关键的水质变量，建立太湖富营养化模型，对太湖富营养化特征进行数值模拟研究。 2 生态水动力学模型 如图1所示，本文主要考虑藻类（叶绿素）、总磷以及总氮三个水质变量，以及藻类生长、磷循环、氮循环等生态动力学过程。通过将这三个水质变量的生化源项放入输移扩散模型，结合二维水动力模型，建立了生态水动力学模型。下面分别给出太湖二维浅水模型和生态水质模型控制方程。 图1 湖泊藻类生态动力学 2.1 水动力方程 设由平均湖面算起的水位起伏为，湖泊水流在方向的速度分量分别为和，则考虑地球自转柯氏力和表面风应力的垂线平均二维浅水水动力控制方程为： listnum a listnum a listnum a 其中为平均湖面到湖底的深度；为实际水深。为Coriolis参数，（为地球自转角速度，为计算区域的维度）；坐标x和y分别取正东和正北方向；为重力加速度；为水平涡粘性系数；为风应力，可采用下列公式计算， listnum a 式中为风应力拖拽系数， 为大气密度， 为风速（一般取水面以上10米处的风速值）， 为风向与x坐标的夹角；是底部摩擦力，可用下式计算： listnum a 式中，为Chezy系数。 2.2 物质输移扩散方程 对三维对流扩散方程进行垂向积分，得到沿垂线平均的平面二维对流扩散方程： listnum a 其中：为水质因子的浓度（mg/L）；，为有效水深；、分别是x和y方向的水平扩散系数（）；Q为源汇项（）；表示水质因子的生化项。 2.3 生态动力学方程 总结前人的研究结果，分别给出叶绿素a、总磷、总氮三种水质因子在物质输移扩散模型中的生化项。 (1) 叶绿素a 浮游植物藻类生物量一般可用叶绿素a来表征，藻类生长主要受到温度、光照和营养盐的影响，具体表达式如下： listnum a listnum a listnum a listnum a listnum a listnum a 其中：为叶绿素a的浓度；为藻类生长率，为藻类最大生长率，、、分别为光照、温度和营养盐对藻类生长的限制作用，为光照强度，为最佳光照强度，是藻类生长的温度修正系数，、为总磷和总氮的浓度；、分别为总磷和总氮对藻类生长的半饱和常数；为藻类损失率，包括呼吸率和死亡率；为藻类沉降率，。 (2) 总磷 本文主要