根据各地多年的监测资料,近年来我国多数城市湖泊水体呈严重富营养化状态⋯。2007年,28个国控重点湖(库)中,V类的5个,占17.9%,劣V类的11个,占39.3%旧J。水体富营养化的重要特征就是发生水华。水华是当水体中出现富营养状况并具备适宜的温度、光照、气候及合适的水文条件等有利于藻类生长或聚集的环境条件时,水体藻类大量生长繁殖或聚集并达到*定浓度的现象。水华*旦发生,就会使水体透明度下降,溶解氧降低,水体出现黑臭等现象,而有些类型的水华还会产生藻毒素,给人类居住环境和人体健康造成损害。水华的发生是突发性的,而水华*旦发生,控制难度就会加大,治理成本成倍提高,因此如果能够预见到水华的发生并及时采取相应措施会取得事半功倍的效果。水华预警是水质预警中的*种突发型预警类别,是指在*定范围内,对藻类生长状况进行分析、评价,对其未来发展状况进行预测。水华预警系统具有超前性预报的功能,能够提前预测出水质演化趋势、方向、速度和后果,在发生水华之前及早发出警报,为水华控制提供科学依据。国内外对水华预警的研究主要围绕三个方面展开:(1)利用单变量或多变量营养指标对水体营养程度进行预测;(2)利用水质模型对水体富营养化程度进行模拟和预测;(3)利用地理信息系统或遥感系统对水华的发生进行预测。水华预警的方法有模糊评价法、人工神经网络、遗传算法、支持向量机(support vector machines,SVM)等,而神经网络方法由于具有较强的适应能力、学习能力和真正的多输人多输出系统的特点而受到人们的重视。
1 预警因子
对于*个确定的水体环境,藻类数量的某*个特定值,是对应了水体环境的某*个状态,这些状态值与藻类的生长需求密切相关,包括营养因子、环境
因子和生态因子,例如TN、,I’P、光照强度、pH值、溶解氧、氧化还原电位等。因此藻类数量的变化可以在这些状态值的变化中表现出来。同时藻类数量的变化也是这些状态值综合作用的结果HJ。水华的发生是由于许多因素如营养盐、水温、光照、pH值、生物因素等共同作用的结果,发生时又有多种水质指针如pH值、溶解氧、氧化还原电位、氮磷浓度等同时发生变化,因此在水质预警模型中参数较多,在模型设计时需要在这些影响因素中筛选出合适的因素作为预警因子。*般选择那些受周围环境影响小的、适合于所选择的模型的、监测方便并且与藻类生长密切相关的这些因素作为预警因子。在众多的影响因子中,氮磷及氮磷比、pH值、DO和ORP等都与藻类生长密切相关,且与藻类的生物量之间具有很好的相关性,可以作为预警因子。
1.1氮磷及氮磷比
氮磷是藻类生长重要的营养因子。氮磷及氮磷比与藻类的生长之间有很好的相关性。可以利用它们之间的这种相关关系建立预测模型,通过监测氮磷浓度来预测水华是否发生。V.H.Smith通过研究指出:对藻类生长来说,总氮质量浓度P(tn)与总磷质量浓度P(tp)之比在20:1以上时,表现为氮过量,磷为限制因子,藻种群密度高峰值主要受磷含量的影响;当小于13:1时,表现为氮不足,氮成为限制因子,藻密度高峰值主要与氮含量有关。本湖沼学者坂本研究发现,当湖水的总氮质量浓度P(tn)与总磷质量浓度P(tp)之比在10:1~25:1范围时,藻类生长与氮和磷的浓度存在直线相关关系。王志红等通过研究初始总氮、总磷、氮磷比等营养因子对“水华”生物量的影响,发现当初始总氮质量浓度小于2.0mg/L、初始总磷质量浓度小于0.1mg/L时,藻生长高峰值与总氮总磷质量浓度比之间具有良好的相关性,并提出了不同氮磷比值与对应藻类“水华”生物量回归模型,可以对藻类“水华”的生物量进行预测。钟卫鸿等"1研究了N、P等对铜山源水库*势藻类(绿球藻和舟形硅藻)生长的影响,发现当N/P为25时,绿球藻生长量*高。杨广杏等对里湖水体进行实地调查,分析水中浮游藻类叶绿素a(chla)与氮磷营养盐含量,进行回归统计,发现它们之间是正相关关系,N、P营养盐变化趋势与浮游藻类叶绿素a变化趋势相吻合。
1.2 pH值
水体pH值与藻类生长关系密切。在碳源丰富的水体中,藻类光合作用影响缓冲体系,从而影响水体pH值。*常见的是藻类大量吸收CO:引起水体pH上升,同时部分藻类对水体中有机酸的吸收和重碳酸盐的利用,也会引起pH的升高归1;而藻类的呼吸作用产生的CO2溶于水中促进H+的生成,会引起pH下降。水体酸碱度也会影响藻类的生长,例如碱性环境有利于藻类光合作用,因为碱性系统易于捕获大气中的CO2,因而较高的生产力往往出现在碱性水体中。每种藻类都有其适合的pH值范围,因而pH不仅会影响藻类的生长繁殖速度,还会影响种类的演替。刘春光等引研究了淡水藻类在不同pH下的生长情况和种类变化。研究结果表明,在pH 8.0—9.5的范围内,pH=8.5藻类生长状况*好,pH=9.5生长*差,表明藻类有适合其生长的pH值,且人为改变pH值会影响藻类的生长。王志红等研究了水库水藻类生长过程中pH的变化,发现pH值随藻类数量的增长呈现出有规律的变化,并建立了藻类数量与pH之间的数学模型,可以通过监测pH来预测藻类的水华现象。游亮等以北京什刹海原水作为培养对象,进行了*系列实验,研究水体中pH的变化,并对实验过程中的pH与藻类生长情况进行相关性分析,发现它们之间的相关系数为0.9312。这也说明了pH值与藻类生长关系非常密切,是*个合理的水华预警参数。
1.3溶解氧(DO)和氧化还原电位(ORP)
水体中DO含量受多因素影响,例如水温、溶解离子、微生物等,而在富营养水体中,DO则主要受生物过程的控制。当藻类数量上升到*定数量*时,其数量的多少、生命活动的旺盛程度对水体的DO变化起主导作用。氧化还原电位是反映介质(土壤、天然水、培养基等)氧化还原状况的*个指标,在湖泊形成水华期间表面水华会造成水体中溶解氧含量降低,氧化还原电位也会随之降低,从而改变藻类的生长环境。张民等钊研究了铜绿微囊藻和栅藻在单培养
和混合培养条件下降低氧化还原电位对两种藻的影响,结果表明单培养下,降低氧化还原电位对两种藻的生长速率没有影响;在混合培养条件下,降低氧化还原电位提高了铜绿微囊藻的生长速率,而降低了栅藻的生长速率。同时试验也发现,ORP降低使得铜绿微囊藻体积变大,生理参数发生改变,有可能是铜绿微囊藻迅速增值的原因。S.Marsili.Libelli对Orbetello礁湖的为期*年的监测数据进行分析,比较了春季和冬季的两个变量的每日变化情况,发
现在春季的时候DO和ORP在每天的不同时段变化很大,呈现周期的变化趋势;而在冬季DO和ORP在每天的不同时段变化不大。由此得出这样的结论:在春季的时候水体温度上升,藻类开始复苏,水体中藻类数量不断增加,藻类的生长影响了DO和ORP,使之呈现周期性的变化趋势。冬季水温低,水体中藻类非常少,使得DO和ORP的变化不是很明显。同时说明,DO、ORP和藻密度是相互影响的,因此可以通过检测DO、ORP数据来对藻类的数量进行预测。
2 国外水华预警模型研究进展
上世纪70年代湖泊学家们通过建立简单的磷负荷模型,用于评价、预测湖泊水体的营养状态。这类模型的典型代表是加拿大湖泊专家Vollenweider提出的Vollenweider模型。Vollenweider模型假定,湖泊中随时间而变化的总磷浓度值等于单位容积内输入的磷减去输出的磷及其在湖内沉积的磷。至80年代,随着对水华和富营养化研究的不断深入,不少专家建立了*系列藻类生物量与营养物质负荷量之间的相关经验模型,其中比较经典的有Rast和Lee的经验模型。这类经验模型简单直观,使用方便,但都假定水体混合均匀、稳态,且限制性营养物质是唯*的,与实际情况往往有较大差别,更不能反映藻类生长的机理。进入90年代后,国外出现了更多对湖泊藻类的预测模型,较有*的有PACGAP类型即藻类种群生长和生产力的预测模型和PROTECH-2类型即浮游植物与环境因子关系模型[20。22I,以及由美国国家环境保护局提出*种多参数综合水质生态模型(water quality analysis simulation program,WQASP)。此后,越来越多的有针对性的水华预警模型被建立并得到成功应用。Scardi和HardingⅢ1研究美国东部的切萨皮克海湾(Chesapeake Bay)的富营养化问题时,采用多层传感器,运用概括方法构造了两个人工神经网络模型,对处于富营养化的初*生产力进行了成功预测。Bin Wei心纠等建立了Kasumi—gaura湖的多因子水质关系模型,利用人工神经网络成功预测到了几种主要*势藻微囊藻(Microcystis)、席藻(Phormidium)和针杆藻(Synedra)的爆发。Friedrieh Recknagel等心刮根据12年以上的环境监测资料,利用神经网络建立了四个系统的淡水水华预测模型,对藻类发生的时间、数量*等的成功预测显示该类模型对复杂的非线性的生态现象进行预测的准确率达到了很高的程度。此外,Nitin Muttil,Joseph H.W.Lee悼刊运用遗传算法对香港铜锣湾3年的叶绿素a、溶解氧和气象水文资料建立实时预测模型进行水华超前预测,也得到令人满意的效果。S.Marsili.Libell利用15个月的水质监测值(每日溶解氧(DO)、pH值、氧化还原电位(ORP)和温度等)的变化来预测水华发生的可能性,运用模糊评
价的方法建立了Orbetello湖的水华预警模型,进行了成功的水华预测。
3 国内水华预警模型应用研究
国内对水华预警的研究起步较晚。近年来随着水体水质的恶化及不断发生的严重水污染事件,人们对环境问题越来越重视。国内对水质预报的研究工作已经全面展开,但更多着眼于大流域的水质预警和湖泊的综合水质预报。朱继业等口引在研究物元分析理论的基础上,运用综合评判模型对南京市秦淮外河进行综合水质评定,并建立回归预警模型进行综合水质预报,在实际应用中取得了较满意的结果。董志颖等悼引采用模糊综合评判法对吉林地区的潜层地下水水质进行预警评价后,结合GIS系统得出了该地区的水质预警结果图。王东云等运用多层前馈神经网络模型和B.P算法,对我国某海域的水质富营养化水平进行了评价,只要将观测结果提供给网络,模型可自动将评价结果输出。刘载文等旧¨利用算法改进型的BP(back propagation)神经网络,选择叶绿素含量、磷、氮磷比、电导率和水温五个参数作为模型输人,以预测ld:3d和5d后的叶绿素含量为目标,构建了北京市长河水系水华短期预报系统,对该水系三个周期的预测精度分别达到了97.2%、94%、88.3%。王洪礼等∞21利用支持向量机理论对海水水质富营养化的程度进行评价,并与BP人工神经网络方法所得结果进行比较,发现SVM理论能更好地解决小样本的分类评价问题,评价效果良好,在水质评价领域有较好的应用前景。韩涛等以MATLAB为工具,建立了评价湖泊水体富营养化状态的BP神经网络模型,应用此模型对我国9个湖泊富营养化程度进行评价。通过对比用分*评分法、模糊数学法、Fuzzy—Grey决策法(F.G决策法)的评价结果,BP神经网络的评价结果更为准确。由于采用了足够多的学习样本对网络进行了训练,*大限度地避免了人为主观因素的影响,并经过样本的检验证明了网络具有很强的泛化能力,所以其评价结果更客观、可靠。曾勇等旧4o采用决策树方法和非线性回归方法建立湖泊水华预警模型。应用决策树方法预测水华爆发时机,非线性回归方法预测水华爆发强度。以北京“六海”为例,利用分段线性多元统计回归预测公式,建立了三个由叶绿素a、水量Q、水温r以及总磷TP组成的回归方程。通过这几个回归方程来计算叶绿素a的含量,从而达到预测水华的目的。近年来,随着科技的发展,更多的高新技术应用于水质预警中。GIS,RS应用比较广泛并取得了良好的效果。窦明等旧纠综合运用GIS、RS、网络、多媒体及计算机仿真等现代高新科技手段,对汉江流域的地形地貌、水质状况、生态环境、水资源分布等各种信息进行动态监测与处理,建立全流域水质基础信息平台、不同功能的水质模型及其相应的管理系统。汉江水质预警系统具备对汉江水质实时监控、水污染事故应急响应、水资源*化调度和水环境综合管理等功能。朱灿等以GIS和数据库管理系统(DBMS)为开发平台,建立了数字西江水质预警预报系统,在发生水污染突发事故后,能够快速预报污染物向下游的扩散时间、扩散面积、确定污染范围、污染程度及对下游取水口等所造成的影响,为决策部门提供决策支持。丰江帆等13刊针对太湖的蓝藻爆发引起的太湖水质不断恶化,结合预警模型和GIS技术以太湖历年来的连续监测资料为基础,运用多元逐步回归统计方法,选择水温等多项环境理
化因素与叶绿素a、藻类生物量、蓝藻生物量等生物因素进行逐步回归分析,建立起多元逐步回归方程,对太湖藻类生物量的变化情况进行预测预报。
4存在问题及研究方向
近20年来,富营养化模型和水华预警模型得到了很大发展:状态变量由*初的几个发展到几十个;水体维数由*维稳态发展到多维动态;研究角度由简单的营养盐吸收发展到对生态系统分析模拟;研究对象由单*的藻类生长模拟发展到综合考虑水体的动力学、热力学及生物动态过程等。但是,在建模过程中仍存在许多问题:(1)建模所依据的数据量不足,缺乏统*详细的水体水化学方面的数据,这给模型的校正、验证造成很大困难,降低了水华预警模型的可靠度和适用性。(2)模型缺乏真正生态系统所具有的灵活性,不能实时模拟环境的突变,因而预测结果不能反映水体生态系统的真实性。(3)模型的模拟对象主要是营养盐的循环、浮游植物的生长和死亡的动态过程,水华预警模型在整个生态系统中非常独立,没有形成*整套水体管理决策支持体系。为了克服上述问题提出了*些新的方法。例如:用模糊数据方法克服数据量不足的问题,用人工智能方法进行参数估计,用混沌与分形理论增强参数估计的能力,用灾变理论模拟系统结构变化,建立生态参数数据库,使用目标函数等。富营养化模拟的发展趋势以学科相互渗透与交错为主,如水体物理环境与藻类生态行为相结合,藻类生态学与分子生物学相结合。富营养化模型也将从单*的预测和评价发展成为多目标管理*化模型。随着新技术的发展,*些新的研究思路和技术也开始逐渐应用到水体富营养化模型中,比如在模型中综合考虑社会学、经济学和心理学因素,结合人工智能方法或GIS技术,从而使模型的适用性和可靠性得到进*步加强。综合运用GIS、RS、计算机仿真技术、多媒体技术等,对水体预警因子进行实时监测与处理,建立更符合水体实际的多维动态水华预警模型,进*步提高预警模型的可靠性,是水华预警模型研究的方向。