1. 引言

近年来湖泊水库的富营养化问题日益突出，影响人们生产生活及供水安全。学者们针对富营养化机理和预测模式，开展了大量研究工作。富营养化水质模型是研究水质动态变化的有力工具，为富营养化研究提供便捷的方法，为水华的防治提供有力的工具。自20世纪70年代Vollenweider提出了简单总磷模型以来，经过众多科学家的努力，富营养化模型已经得到了很大程度的发展。近年来，富营养化机理性模型中发展较为完善的是生态动力学模型。该类模型涵盖了生态系统中影响藻类生物量的大部分生态学过程，较为完善。国外已有相关模型[1] [2] ，如PACGAP 模型,即藻类种群生长和生产力的预测模型；PROTECH-2模型，即浮游植物与环境因子关系模型[3] 。也有部分商业软件，如美国EPA开发的WASP中富营养模块EUTRO和CE-QUAL2。自上世界80年代以来，我国学者在湖泊富营养化模型构建方面也做了大量工作，特别是在一些富营养化较为严重的湖泊，如太湖、滇池、巢湖等，构建了较好的富营养化生态动力学模型。

这些生态动力学模型虽然尝试从机理上反映藻类生物量的变化，但涉及的生态学过程复杂，参数众多，对数据量要求大，并且取样频次与取样时间对模型参数测定有较大影响，给模型的参数率定和验证造成了很大困难，且多数主要预测藻类总生物量或叶绿素a，专门针对蓝藻水华的模拟预测模型尚不多见。

而采用逐步多元回归方法，分析环境因子与藻类相关关系，方法简单易行，在实际监测资料的基础上，能较好的明确在特定水体中，环境因子与富营养化藻类的关系，为水华藻类机理模型提供较好的参数分析基础。国内学者们分别在太湖、洪泽湖、千岛湖、淀山湖、巢湖等水体[3] -[12] 进行了相关研究。但由于不同水体有各自的特性，需结合各自的情况具体分析。

深圳市茜坑水库是深圳市重要的饮用水水源地之一，富营养化问题一直是该水库重点防治问题。为了给该水库管理部门防治富营养化提供理论支持，本研究采用该水库连续5年的监测资料，采用逐步回归的统计方法，分析藻类生物量与相关环境理化因子的相关关系，研究对藻类生物量和水华蓝藻生长影响较显著的因子，为水库水华防治工作提供技术参考。同时筛选关键参数，为建立该水库富营养化蓝藻预测模型提供一定的基础。

2. 研究方法

2.1. 数据来源

深圳市茜坑水库是深圳市重要的水源地之一，为全面监控水库的水质状况，该水库管理部门自2000年开始，在水库布设采样点，每月采样1~2次，分析了氨氮、总氮、总磷、溶解氧、透明度等化学指标，以及蓝绿藻及藻类总生物量等生物指标，取得了丰富的监测资料。水库共布设6个监测点，各点位置布设如图1所示。

由于监测点VI为后来设立，数据不连续，所以选用监测点I~V连续5年的监测资料进行分析。选择的环境理化因子包括：引水量、水深、水温、透明度、pH值、化学需氧量(COD)、溶解氧(DO)、氨氮(NH₃-N)、总氮(TN)、总磷(TP)等，其中溶解氧包括表层、4m水深处和8m水深处数据。生物因子分析包括叶绿素a、藻类生物量和蓝藻的生物量。

2.2. 统计分析

借鉴陈宇炜等[3] [5] [13] [14] 的方法，采用SPSS软件，构建回归方程，具体如下：

(1) 采用SPSS软件，对各采样点，以各理化因子为自变量，叶绿素a为因变量，选择临界值95%的显著水平(FC = FJ = 3.0)筛选自变量[3] [4] ，得到回归方程和复相关系数，分析得到各采样点的与叶绿素a关系最接近的环境因子。

(2) 在上述分析结果基础上，从理化因子中，剔除pH，DO，COD，SD等5项，以剩余6项理化因子为自变量，以藻类总生物量和蓝藻生物量为因变量，采用SPSS分析，分别得到两项生物因子(藻类总生物量和蓝藻生物量)的逐步回归预测方程和复相关系数等。

(3) 对11项理化因子进行标准化处理，以其标准化值为自变量，以藻类总生物量和蓝藻生物量的标准化值为因变量，采用SPSS分析，得到两项生物因子(藻类总生物量和蓝藻生物量)标准化值的逐步回归预测方程和复相关系数等。

3. 结果和讨论

3.1. 环境因子对藻类生物量的影响

本研究中涉及的环境因子均对藻类的生长有着间接或直接的影响作用。例如，水温对藻类生长有限制作用，氮、磷营养盐通常是藻类生长的限制性因子，透明度和悬浮质的变化会影响水下日光辐射的变化进而影响藻类的生长 [3] ，同时，透明度又会因为藻类的生长而受到影响。通过对各个因子的简单相关分析，可以看出水库理化因子和藻类生物量的线性相关关系(见表1)。

从上表中可知，蓝藻为水库中的主要的藻类，其与藻类生物量的相关性最大，达0.83。而在理化因子中，藻类生物量与COD的相关性较大，说明在水库中COD对藻类的生长有较大影响。

在本研究中，对DO的分析分别设定了表层DO、4 mDO、4 mDO和DO均值这四个指标。其中，4 mDO与DO均值相关性最大，表明了在水库中，4 mDO的含量最大。而4 mDO与藻类生物量的相关因子达0.57，说明4 mDO是又一影响藻类生长的环境因子。

营养盐的影响：通过大量的研究，人们一致认为，湖泊中的氮、磷等营养物是藻类生长的主要限制性因素。在本研究中，TN和TP和藻类生物量的影响因子分别为0.41和0.38，影响较为密切。

3.2. 叶绿素a 与理化因子的逐步回归分析

在不同的水体中，环境理化因子对藻类生长的影响程度各不相同。通过多元逐步回归分析，可以筛选出相对重要的影响因子。藻类叶绿素a与理化因子的逐步回归统计结果如表2所示。

由表2可以看出，各个监测点筛选出的重要因子各不相同。

表1. 水库一些理化因子和藻类生物量的线性相关关系

表2. 叶绿素a与理化因子的逐步回归统计

P1点共有4个因子入选，其中COD和TN为正相关，引水量和NH3-N为负相关。

P2点有COD和DO均值这两个因子入选，且均为正相关。

P3最为特殊，只有表温这个因子影响较为显著。其复相关系数和综合F值均最低，分别为0.54和5.751。

P4有COD和透明度这两个因子入选。其复相关系数和综合F值均最高，分别为0.966和63.641。

而P5有COD、表层DO、水深和透明度这四个因子入选，其中透明度为负相关，其他因子为正相关。

P3点为引水入口附近监测点，受引水影响较大，P1、P2、P4和P5均有COD这个因子入选，而且均为正相关。表明，各个理化因子中，COD对叶绿素a的影响最为显著。

3.3. 藻类生物量和蓝藻生物量的逐步回归分析

上述分析结果是假设所有理化因子与叶绿素a之间是单向的线性相关，忽略了藻类生长或者大规模爆发，也会影响水体理化因子，如：藻类等光合作用能产氧气，影响水体溶解氧的变化，藻类生长过程中产生的化学物质，也会影响水体中PH的变化，这些环境因子可以说是以藻类为“因”而形成的“果”因素，难以作为自变量用于预测藻类生长 [3] [4] 。

因此，要预测藻类生物量尤其是蓝藻生物量的变化，需要剔除溶解氧、pH、SD等因素，用藻类总生物量和蓝藻生物量分别作为因变量，计算逐步回归方程(见表3)。

由表3可见，监测点4无法建立回归方程，监测点1，2，3，5的回归方程分别与水深、TP，表温有关。各个方程的复相关系数相差不大，表明，水深、温度和营养物质对藻类生物量的影响相当。

对蓝藻生物量与环境因子进行回归分析，发现，只有监测点P2能建立回归方程，入选变量为TP，回归方程为：BLUE = −0.016 + 6.187 × TP，复相关系数为0.482，综合F值为6.353。

3.4. 标准化回归分析

以监测点P2为例，对藻类生物量和各环境因子进行标准化处理，标准化后的逐步回归方程以TP为显著因子，回归方程为：

表3. 藻类生物量与相关环境因子的逐步回归统计结果

Ln(PHY + 1) = 0.659 + 2.518 × Ln (TP + 1)

复相关系数为0.566，与标准化前相近。

同样以监测点P2为例，对蓝藻生物量和各环境因子进行标准化处理，标准化后的逐步回归方程以水深为显著因子，回归方程为：

Ln(PHY + 1) = −1.078 + 0.616 × Ln (Depth + 1)

其中，复相关系数为0.469，显著性较差。与未标准化回归分析相比，未标准化回归方程显著因子为TP，更符合实际情况。

4. 结论

对藻类生物量与环境理化因子的相关性回归分析，结果表明，COD的影响最为显著，但化学需氧量受藻类生长的影响而变化，目前无证据表明COD影响藻类生长，因而难以用于预测藻类生长情况。剔除COD，DO等次要环境因子外，藻类生物量主要受水深，水温及TP的影响。但不同监测点的显著因子不同。表明，藻类生长因环境不同，限制因子也不同。为了建立蓝藻预测模型，对各个因子及藻类生物量进行标准化处理，标准化后的回归方程相对未标准化方程，相关性并无显著改善。

综上所述，逐步回归统计模型，能够寻找影响藻类生物量的关键性因子，为准确预测提供参考，但尚不能对藻类爆发性生长做出精确的预测，要对该水库蓝藻水华进行较准确的模拟预测，还需在分析显著因子的基础上，结合藻类生长机理，建立预测模型。

基金项目

国家基础科学人才培养基金《武汉大学地理科学理科基地》科研能力训练项目(J1103409)；中央高校基本科研业务费专项资金武汉大学自主科研项目(121066)。

参考文献