基于EOS_MODIS资料的太湖藻类动态模拟

基于ＥＯＳ／ＭＯＤＩＳ资料的太湖藻类动态模拟颜润润，等

基于ＥＯＳ／ＭＯＤＩＳ资料的太湖藻类动态模拟

颜润润，

逄勇

２１００９８）

（河海大学环境科学及工程学院，南京

摘

要：建立了太湖藻类生长的动态模型，并将其与水动力模型和水质模型相耦合，利用２００１年７～８月太湖的实测资料对模型

进行了率定。借助该模型，对２００４年８月太湖水体中ＴＮ、ＴＰ的变化以及藻类生长过程进行了模拟。利用中分辨率成像光谱仪

ＥＯＳ／ＭＯＤＩＳ的数据对太湖叶绿素ａ浓度进行遥感定量，将遥感监测数据和模型计算结果进行了比较。结果表明：该模型可进行风生湖流、ＴＮ、ＴＰ的模拟，以叶绿素ａ浓度描述的藻类浓度的模拟值能较好地拟合遥感监测值，且遥感监测图和模型模拟图所反映的全太湖叶绿素浓度分布基本一致。最后根据遥感和模拟对太湖全区的藻类分布作了具体的分析。

关键词：太湖；

叶绿素ａ浓度；

ＥＯＳ／ＭＯＤＩＳ数据；遥感监测；动态模型

中图分类号：Ｘ８３５文献标识码：Ａ文章编号：１００３－６５０４（２００７）０７－００２９－０４

ＥＯＳ卫星是美国对地面进行连续和综合观测的系列低轨卫星，目前已发射了上午星Ｔｅｒｒａ和下午

星Ａｑｕａ。这两颗卫星都装有中分辨率成像光谱仪－

供有效的技术支持。此外，通过对遥感数据及模拟结果的分析，更全面的掌握全太湖的藻类分布状况。

ＭＯＤＩＳ，它是ＥＯＳ系列卫星上的主要仪器，同时也是ＥＯＳＴｅｒｒａ平台上唯一进行直接广播的对地观测仪

器。ＭＯＤＩＳ是当前世界上新一代“图谱合一”的光学遥感仪器，共有３６个光谱通道，光谱范围从０．４Ｌｍ（可见光）到１４．４Ｌｍ（热红外），覆盖了ＡＶＨＲＲ和ＳｅａＷ－

１１．１

水量水质藻类生长模型建立水量模型

假设湖水为均匀不可压的流体，垂直方向上服从

静水压力分布。采用笛卡尔左手直角坐标系，见图１。ｘ轴和ｙ轴位于湖水的平均水平面上，ｘ轴向东为正，ｙ轴向北为正，ｚ轴向下为正，得其流体动力学方程为：

ｉＦＳ的光谱通道。ＭＯＤＩＳ第１和第２通道的空间分辨率为２５０ｍ，第３到第７通道的空间分辨率为５００ｍ，其余通道的空间分辨率为１０００ｍ，扫描幅宽达２３３０ｋｍ［１－２］。ＥＯＳ／ＭＯＤＩＳ卫星数据以其较高的空间分

辨率和时间分辨率，以及星载多个专业海洋水色传感器的特点成为水色遥感研究的首选数据源。

本文在前人研究基础上建立太湖藻类生长的动态模型，并将该模型与水动力模型和水质模型相耦合，综合考虑了风场、水动力特征、光照、水温等环境因子的影响。以２００１年７～８月太湖２３个采样点同步监测的水文水质资料对模型进行了率定。借助该模型，对２００４年８月太湖水体中ＴＮ、ＴＰ的变化以及藻类生长过程进行了模拟。从模拟结果得知，各采样点ＴＮ，ＴＰ以及藻类叶绿素ａ的模拟值能较好地拟合实测值，模型对各状态变量的描述是理想的。利用中分辨率成像光谱仪ＥＯＳ／ＭＯＤＩＳ的数据对太湖叶绿素ａ浓度进行定量遥感，将该遥感结果和模型模拟结果比较得知：遥感监测图和模型模拟图所反映的全太湖叶绿素ａ浓度分布基本一致。说明该模型对太湖富营养化的产生过程具有较好的描述能力，能为太湖富营养化的研究提

基金项目：国家自然科学基金（５０５７９０１５）；国家“８６３”计划资助（２００３ＡＡ

!ｕ＋ｕ!ｕ＋ｖ!ｕ＋ｗ!ｕ＝ｆｖ－１!ｐ＋

!ｔ!ｘ!ｙ!ｚ

２

!２ｕ］＋Ａ!２ｕ＋（１）Ａｈ［!ｕｚ

!ｖ＋ｕ!ｖ＋ｖ!ｖ＋ｗ!ｖ＝－ｆｕ－１!ｐ＋２

!２ｖ］＋Ａ!２ｖ＋（２）Ａｈ［!ｖｚ

!ｕ＋!ｖ＋!ｗ＝０（３）（４）ｐ＝!ｗｇ（

式中：ｕ、ｖ、ｗ分别为ｘ、ｙ、ｚ轴方向上的流速分量；η为垂直方向上湖面相对与平均水面的高度；ρＷ为水体密度；Ａｚ、Ａｈ分别为垂直和水平涡动粘滞系数；ｆ为柯氏参数，ｆ＝７．２３×１０－５；ｇ为重力加速度；ｐ为水的压强。模型的边界条件、初始条件、差分格式等见文献［３－５］。１．２水质模型

水质方程是以质量平衡方程为基础的。由于三维水质输移方程包含很多不可确定的参数，在现有条件下，模型的率定验证存在困难，故采用二维水质模型。二维水质迁移扩散方程为［６－７］：

６０１１００－２）

作者简介：颜润润（１９８２－），女，博士研究生，从事水环境评价及规划的研究工作，（手机）１３９５１６２４９２１（电子信箱）ｙｒｒ１１３７＠ｈｈｕ．ｅｄｕ．ｃｎ。

!Ｃ＋Ｕ!Ｃ＋Ｖ!Ｃ＝!（Ｅ!Ｃ）＋

ｘ

!ｔ!ｘ!ｙ!ｘ!ｘ!（Ｅ!Ｃｉ）－ＫＣ＋Ｓ

ｘｉｉｉ

（５）

环境科学与技术

第３０卷第７期２００７年７月

式中：Ｃｉ为ＴＮ、ＴＰ的浓度；ｕ、ｖ分别为ｘ、ｙ方向上的流速分量；Ｅｘ、Ｅｙ分别为ｘ、ｙ向上的扩散系数；Ｋｉ为ＴＮ、ＴＰ的降解系数；Ｓｉ为ＴＮ、ＴＰ的底泥释放项。

１．３藻类生长动态模型

由于浮游植物的大量繁殖是富营养化水体的重

要标志，也是造成危害的主要因素，而叶绿素ａ是藻类现存量的重要指标，水体叶绿素ａ含量的高低能够反映水体的营养状况。因此，叶绿素ａ常作为湖泊富营养化调查的一个重要参数，并且在水体富营养状况评价中起关键作用。故本文以叶绿素ａ的含量来说明水体藻类的含量及评价水体的富营养水平。根据质量守恒原理可得富营养化变量的基本方程为［６，８］：

!Ｃ!ｔ＋Ｕ!Ｃ!ｘ＋Ｖ!Ｃ!ｙ＝!!ｘ（Ｅｘ

!!ｘ）＋

!（Ｅｙ!Ｃｃｈｌ－ａ）＋Ｓｃｈｌ－ａ（６）式中：Ｃｃｈｌ－ａ为叶绿素ａ的浓度；Ｕ、Ｖ分别为ｘ、ｙ轴方向上的流速分量，由水量模型计算得出，Ｅｘ，Ｅｙ分别为藻类的纵向扩散系数和横向扩散系数；Ｓｃｈｌ－ａ为叶绿素ａ的转化项。其中：

Ｓｃｈｌ－ａ＝ＧＰ１（ｔ）－ＤＰ１（ｔ）－Ｖｓ（７）

ＧＰ１（ｔ）＝ｕｍａｘｆ（Ｔ（ｔ））ｆ（Ｌ（ｔ））ｆ（ＴＮ（ｔ））ｆ（ＴＰ（ｔ））（８）

ｆ（Ｔ（ｔ））＝ｅｘｐ（－２．３×｜Ｔ（ｔ）｜－Ｔｏｐｔ）（９）ｆ（Ｌ（ｔ））＝Ｌ（ｔ）（１０）

Ｌ

ｆ（ＴＮ（ｔ））＝ＴＮ（ｔ）

（１１）Ｎ

ｆ（ＴＰ（ｔ））＝ＴＰ（ｔ）

（１２）Ｐ

Ｄｐ１（ｔ）＝Ｋ１Ｒ!１ＲＴ－２０＋Ｋ１Ｄ（１３）式中：ＧＰ１为藻类生长项；ＤＰ１为藻类死亡项；Ｖｓ为藻类的沉降项；Ｔ为温度；Ｔｏｐｔ为藻类生长的最佳温度；Ｄ为水深；ＫＮ、ＫＰ为氮和磷的半饱和常数；Ｌ为水面上的光照强度；Ｋ１Ｒ为藻类的呼吸率；θ１Ｒ为系数；Ｋ１Ｄ为藻

类死亡率。

１．４模型的率定

本文采用２００１年７月、８月全太湖２３个采样点

每个月１次实测的气温、

水温、日照、风速、风向资料和水质监测资料对模型进行了率定。这２３个采样点分别位于污染严重的梅梁湖区域、流速较大的太湖西岸区域、湖心区和污染较轻的东太湖区域等，遍布了整个太湖，具有一定的代表性，具体的分布位置见图１。

水动力模型和水质模型的时间步长均取为１２０ｓ，取２００１年７月的入太湖河道的实测流量、水位、流速资料分别作为水动力模型的开边界条件，各入湖河道的排污量和主要排污口的排污量作为水质模型的开边界条件，以２００１年８月实测的气象、水文、水质资料作为判别基准，判别要素ＴＮ、ＴＰ和叶绿素ａ三项，率定结果见图２。模型参数共１２项，符号、单位、取值如表１。从图２可以看出：模型的率定结果与实测结果都比较接近，

说明率定所得的参数是合理的

，

可以用来进行预测。

表１模型主要参数率定结果

参数符号

单位

数值

方法总磷的降解系数

ＫＴＰ（ｌ／ｄ）

０．０３０

文中率定总氮的降解系数ＫＴＮ（ｌ／ｄ）０．０３４文中率定磷的半饱和常数ＫＰ（ｍｇ／Ｌ）０．０２５文中率定氮的半饱和常数ＫＮ（ｍｇ／Ｌ）

０．３４文中率定温度系数!１Ｒ１．０６８文中率定藻类的呼吸率Ｋ１Ｒ（ｌ／ｄ）０．４文献［６］藻类最大生长率

３０文献［６］藻生长的最佳光强的半饱和常数ＫＬ（

３００文献［４］藻类最大死亡率Ｋ１Ｄ（ｍｇ／Ｌ）０．１７

文献［４］

２藻类动态模拟及验证２．１

太湖水质遥感结果

根据太湖水体的光谱特性，选用ＥＯＳ／ＭＯＤＩＳ第

９波段（４３８～４４８ｎｍ）数据开展研究，利用该波段数据湖面离水反射率的不同来区分太湖流域叶绿素浓度的分布。通常湖水叶绿素浓度越大，第９波段离水反射率的值越高，将该波段的离水反射率值划分出不同的区间，并根据区间内值的范围来监测太湖流域水体

叶绿素的浓度分布。将２００４年８月１０日上午１０时４０分ＥＯＳ／ＭＯＤＩＳ资料太湖流域６、２、９三通道伪彩

色合成图中第９波段光谱数据离水辐射率值作分段增强，得到ＥＯＳ／ＭＯＤＩＳ卫星数据监测太湖流域叶绿素浓度分布图（见图３），将其与中国科学院南京地理与湖泊研究所的科研人员对太湖多个站点叶绿素浓度实测资料进行对比（见表２）。从表２中可以看出反演值和实测值吻合较好，可见ＥＯＳ／ＭＯＤＩＳ遥感数据的反演结果能够较为准确地反映太湖流域叶绿素的浓度分布情况。

表２２００４年８月１０日遥感数据的反演结果和实测值比较

监测点实测值（ｍｇ／ｍ３）

反演值（ｍｇ／ｍ３）

相对误差（％）

五里湖２１２２００６．０梅园１９６１７５１２．０小湾里１４８１３０１３．８胥口１３．６１３４．６漫山

２４．２

２５

３．２

２．２模型模拟结果

利用上述将藻类生长与水动力、水质相耦合的模

型对太湖２００４年８月太湖的叶绿素ａ浓度值进行了模拟。模型参数利用表１率定结果，风速和风向均取８月实测的平均值：２．８ｍ／ｓ，东南风；边界处的流速和流量取８月的月平均值。各采样点的实测值与计算值的对比见图４。从图４可以看出：无论叶绿素ａ浓度还是

ＴＮ、

ＴＰ

浓度

，模拟值和实测值基本一致。

２．３遥感与模型模拟结果比较及分析

基于ＥＯＳ／ＭＯＤＩＳ资料的太湖藻类动态模拟颜润润，等

利用上述模型模拟出的２００４年８月太湖流域叶绿素浓度分布状况见图５，将遥感监测图（图３）和该模型模拟图比较发现，两者反映的全太湖叶绿素浓度分布基本一致。通过模型模拟和遥感手段得到的太湖中不同浓度叶绿素区域所占的面积见表３，从表３中可以看出：不同叶绿素浓度对应的区域面积也基本一致。夏季太湖流域的叶绿素分布以中、高浓度为主，占很大面积，而低浓度区主要集中在东太湖一带，总体看来所占面积远小于中、高浓度区。

表３２００４年８月１０日遥感和模拟得到不同浓度叶绿素区域

所占面积比较

叶绿素浓度（ｍｇ／ｍ３

）

遥感结果（ｋｍ２）模拟结果（ｋｍ２）相对误差（％）

＞２００５５０．６５００１０．１１５０～２００１８２．２１５０２１．５１００～１５０２００．４１８０１１．３６０～１００４０９．７３５０１７．１４０～６０４２２．９４５０６．０２５～４０２７０．５２８０３．４１５～２５２１５．６１８０１９．８０～１５

７０．３１００２９．７

从遥感监测图和模型模拟结果可以看出：北部湖区均要高于南部湖区，高值出现在梅梁湾入湖河口及太湖沿岸地带，太湖西北面大浦附近也出现高浓度区及次高浓度区，低值出现在湖心区和东太湖，并大致呈现从入湖河口往梅梁湾内、湾口及大太湖、东太湖逐渐递减的趋势，至太湖东南角的东太湖湖区出现最低值。出现这样的分布状况主要原因在于：

（１）在太湖北面，竺山湖、

梅梁湾渔业村等湖区都是养殖带，造成富营养化程度较高，水体叶绿素浓度也相应较高；太湖沿岸的工业生产用水及居民生活用水很大部分直接排入太湖，水体污染与富营养化严重，造成了沿岸水体叶绿素浓度高。

（２）根据调查，由于太湖北面与西北面的无锡市和武进市排放大量工农业和生活污水进入太湖，腐殖质浓度非常高所致，

大量污水进入西太湖后在此处沉淀过滤，因而太湖西北面悬浮物和叶绿素浓度较高。

（３）东北角的贡湖湾大片区域属于太湖的一个草型湖区，水体相对清澈，悬浮物浓度较低，富营养化程度也相应较低，因而叶绿素浓度值低于北太湖和沿岸

（下转第５１页）

工业节水综合效益的分析郝桂珍，等

式中：Ｖ－相应节水水平下单位节水量的投资（万元／ｍ３・（ｍ３／万ｄ）；Ｘ－城市工业平均万元产值取水量元）。

由图１、图２可知节约单位水量投资随万元产值取水量的降低而增大。在节水的初期，由于节水技术相对容易，在节水投资相对来说较少的情况下，可获得较显著的经济效果，这时由于节水的直接经济效益大于节水投资，一般用户都能积极配合当地政府，主动采取节水措施，节水效果显著。

节约用水相当于水资源生产，符合传统工业经济的边际效益递减规律，也就是说，随着节水达到一定水平后，其投入就不再与节水量成正比。因此，就要根据节水曲线进行核算，当节水的投入高于调水或节水的投入高于水价时，这时应考虑调水或调整水价了［７］。

节水投入的限度也靠水价，当节水投入边际效益递减规律起作用的时侯，如果水的供需仍不能平衡，就应该提高水价，使得节水有更大的投入，促进“谁节约，谁受益；谁浪更大量的节水［８］。节水应本着

费，谁受罚”的原则，鼓励各行各业和各个部门更大量的节水。

节水是企业节能降耗的重要途径，除了节省用水量、减少水资源开采外，还减少了污水的排放量，减轻了环境的污染。建议政府部门积极支持工业节水项目，实行优惠政策加以鼓励，以进一步缓解我国水危机，解决工业用水的供需矛盾。

［参考文献］

［１］伍振毅，宋宗文．我国工业节水创新技术的研究进展［Ｊ］．工

业水处理，２００６，１２．

［２］林佩凤，王丽琼，张江山．水资源价值损失模糊估算模型的

改进与应用［Ｊ］．环境科学与技术，２００７，１．

［３］王习元，李文青，刘文英，等．农村河流环境保护市场化探析

［Ｊ］．环境科学与技术，２００６，９．

［４］周彤．污水回用是解决城市缺水的有效途径［Ｊ］．给水排水，

２００１，２７（１１）．

［５］叶维钧．城市基础设施效益分析［Ｊ］．科学技术文献出版社重

庆分社，１９８９．

［６］陈冬辰．宾馆污水处理与回用工程简述［Ｊ］．给水排水，１９９９，

２５（４）．

［７］王文革．循环经济激励制度存在问题与对策［Ｊ］．环境科学与

技术，２００６，１２．

［８］谭绍仲，王运莲．城市供水企业授水价格调价要素分析［Ｊ］．

城镇供水，２００６，１．

３结论

（收稿２００７－０１－２８；修回２００７－０３－２０）

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!（上接第３１页）

地带，而位于太湖东南角的东太湖狭长湖湾有着茂盛的水生植物，能够对水体有过滤、净化、加速悬浮物沉降等作用，使得营养化程度最低，因而东太湖一直是全太湖水质最好的湖区，叶绿素浓度值也相应为最低。

动态变化，通过水质遥感图，可以更全面直观地了解整个区域的藻类分布状况。

［参考文献］

［１］王其茂，马超飞，唐军武，等．ＥＯＳ／ＭＯＤＩＳ遥感资料探测海洋

赤潮信息方法［Ｊ］．遥感技术与应用，２００６，２１（１）：６－１０．

３结论

建立了三维风生湖流模型、水质模型和富营养化

［２］殷青军，扬英莲．基于ＥＯＳ／ＭＯＤＩＳ数据的青海湖遥感监测

［Ｊ］．湖泊科学，２００５，１７（４）：３５６－３６０．

［３］逄勇，濮培民．太湖风生流三维数值模拟试验［Ｊ］．地理学报，

１９９６，５１（４）：３２８－３３２．

［４］许秋瑾，秦伯强，陈伟民，等．太湖藻类生长模型研究［Ｊ］．湖泊

科学，２００１，１３（２）：１４９－１５７．

模型耦合的数学模型，该模型不仅可以对太湖的风生湖流，ＴＮ、ＴＰ、ＣＯＤ等水质要素进行模拟，还可以模拟藻类在太湖中的生长和消亡情况以及随风生湖流迁移的规律。利用中分辨率成像光谱仪ＥＯＳ／ＭＯＤＩＳ的数据对太湖叶绿素浓度进行定量遥感，将该遥感结果用于太湖藻类生长模拟结果的验证，可以从空间分布上验证藻类生长的模拟结果是否切实符合实际情况，这一点用常规实地测量手段难以实现。通过比较得知：模拟值和实测值在分布趋势上和分布规律上都基本一致。说明模型能够较好地描述太湖藻类的时空分布特征，对太湖富营养化的产生过程具有较好的描述能力，从而能为制定太湖富营养化的有效控制措施提供技术支持。此外，水质遥感可以随时掌握藻类的

［５］吴旭光．分裂特征线法及收敛稳定性分析［Ｊ］．河海大学学报，

１９９２，２０（６）：７１－７９．

［６］逄勇，姚琪，濮培民，等．太湖地区大气－水环境的综合数值研

究［Ｍ］．北京：气象出版社，１９９８．４８－７０．

［７］朱维斌．河流混合带污染物质二维浓度分布的计算［Ｊ］．河海

大学学报，１９８９，１７（４）：１０８－１１３．

［８］胡维平，濮培民，秦伯强．太湖水动力学三维数值试验研究

［Ｊ］．湖泊科学，１９９８，１０（４）：２６－３３．

［９］李一平，逄勇，崔广柏．太湖藻类动态模拟［Ｊ］．河海大学学报，

２００４，３２：１８５－１８９．

（收稿２００６－０９－３０；修回２００６－１１－２０）

ＥｎｖｉｒｏｎｍｅｎｔａｌＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙＶｏｌ．３０，Ｎｏ．７，Ｊｕｌｙ２００７

ｏｆｙｅｌｌｏｗｐｈｏｓｐｈｏｒｕｓ．Ｐａｒａｍｅｔｅｒｓｏｆａｄｓｏｒｂｉｎｇｐｕｒｉｆｉｃａｔｉｏｎｗｅｒｅｏｐｔｉｍｉｚｅｄｗｉｔｈｏｒｔｈｏｇｏｎａｌｅｘｐｅｒｉｍｅｎｔ，ａｎｄｍｏｄｉｆｉｅｄａｃｔｉｖａｔｅｄｃａｒｂｏｎｓａｍｐｌｅｓｗｅｒｅｃｈａｒａｃｔｅｒｉｚｅｄｂｙＸＰＳ．Ｒｅｓｕｌｔｓｓｈｏｗｅｄｔｈａｔｍｏｄｉｆｉｅｄａｃｔｉｖａｔｅｄｃａｒｂｏｎｗｉｔｈｃｏｐｐｅｒｃａｔｉｏｎｓｈａｓｔｈｅｈｉｇｈｅｓｔｐｕｒｉｆｉｃａｔｉｏｎｅｆｆｉｃｉｅｎｃｙ．Ｏｐｔｉｍｕｍｐａｒａｍｅｔｅｒｓｏｆｒｅａｃｔｉｏｎｗｅｒｅｔｅｍｐｅｒａｔｕｒｅ９５℃，ｏｘｙｇｅｎｃｏｎｔｅｎｔ３％，ｇａｓｆｌｏｗｒａｔｅ０．２１／ｍｉｎ，ｃｏｎｃｅｎｔｒａｔｉｏｎｏｆＣｕ０．３ｍｏｌ／Ｌ，ｒａｄｉｕｓｏｆａｃｔｉｖａｔｅｄｃａｒｂｏｎ４ｍｍａｎｄｃｒｅｍａｔｉｎｇｔｅｍｐｅｒａｔｕｒｅ３５０℃．Ｋｅｙｗｏｒｄｓ：ＰＨ３；ｍｏｄｉｆｉｅｄａｃｔｉｖａｔｅｄｃａｒｂｏｎ；ｏｒｔｈｏｇｏｎａｌｅｘｐｅｒｉｍｅｎｔ；ｃｈａｒａｃｔｅｒｉｚａｔｉｏｎ

ＤｙｎａｍｉｃＳｉｍｕｌａｔｉｏｎｏｆＡｌｇａｅｉｎＴａｉｈｕＬａｋｅＢａｓｅｄｏｎ

ＥＯＳ／ＭＯＤＩＳＤａｔａＹＡＮＲｕｎ－ｒｕｎ，ＰＡＮＧＹｏｎｇ

（ＳｃｈｏｏｌｏｆＥｎｖｉｒｏｎｍｅｎｔＳｃｉｅｎｃｅａｎｄＥｎｇｉｎｅｅｒｉｎｇ，ＨｏｈａｉＵｎｉｖｅｒｓｉｔｙ，

Ｎａｎｊｉｎｇ２１００９８）

ＡｂｓｔｒａｃｔＩＤ：１００３－６５０４（２００７）０７－００２９－０４－ＥＡ

Ａｂｓｔｒａｃｔ：ＡｄｙｎａｍｉｃｍｏｄｅｌｏｆａｌｇａｅｇｒｏｗｔｈｗａｓｅｓｔａｂｌｉｓｈｅｄｆｏｒＴａｉｈｕＬａｋｅ，ｗｈｉｃｈｃｏｕｐｌｅｄｗｉｔｈｈｙｄｒｏｄｙｎａｍｉｃｍｏｄｅｌａｎｄｗａｔｅｒｑｕａｌｉｔｙｍｏｄｅｌ．ＴｈｅｎｕｍｅｒｉｃａｌｍｏｄｅｌｗａｓｃａｌｉｂｒａｔｅｄｂａｓｅｄｏｎｍｏｎｉｔｏｒｉｎｇｄａｔａｏｆｔｈｅＬａｋｅｆｒｏｍＪｕｌｙｔｏＡｕｇｕｓｔｉｎ２００１．ＷａｔｅｒｑｕａｌｉｔｙｉｎｃｌｕｄｉｎｇＴＮａｎｄＴＰａｓｗｅｌｌａｓｇｒｏｗｔｈｏｆａｌｇａｅｉｎＡｕｇｕｓｔ２００４ｗａｓｓｉｍｕｌａｔｅｄｗｉｔｈｔｈｅｍｏｄｅｌ．Ｃｈｌｏｒｏｐｈｙｌｌ－ａｃｏｎｃｅｎｔｒａｔｉｏｎｗａｓｄｅｔｅｒｍｉｎｅｄｂｙｕｓｅｏｆｍｏｄｅｒａｔｅｒｅｓｏｌｕｔｉｏｎｉｍａｇｉｎｇｓｐｅｃｔｒｏｒａｄｉｏｍｅｔｅｒ－ＥＯＳ／ＭＯＤＩＳｄａｔａ，ａｎｄｗａｓｃｏｍｐａｒｅｄｗｉｔｈｒｅｍｏｔｅｓｅｎｓｉｎｇｍｏｎｉｔｏｒｉｎｇｄａｔａ．Ｒｅｓｕｌｔｓｉｎｄｉｃａｔｅｄｔｈａｔｔｈｅｃａｌｃｕｌａｔｅｄｗｉｎｄｃｕｒｒｅｎｔ，ＴＮａｎｄＴＰａｎｄｍｏｄｅｌｅｄｃｏｎｃｅｎｔｒａｔｉｏｎｏｆａｌｇａｅｄｅｓｃｒｉｂｅｄｂｙｃｈｌｏｒｏｐｈｙｌｌ－ａｃｏｎｃｅｎｔｒａｔｉｏｎｗａｓｉｎｗｅｌｌａｃｃｏｒｄａｎｃｅｗｉｔｈｔｈｅｍｏｎｉｔｏｒｉｎｇｄａｔａ，ｃｈｌｏｒｏｐｈｙｌｌ－ａｃｏｎｃｅｎｔｒａｔｉｏｎｄｉｓｔｒｉｂｕｔｉｏｎｒｅｆｌｅｃｔｅｄｆｒｏｍｒｅｍｏｔｅｓｅｎｓｉｎｇｍｏｎｉｔｏｒｉｎｇｃｈａｒｔａｎｄｍｏｄｅｌｓｉｍｕｌａｔｉｎｇｃｈａｒｔｗａｓａｌｓｏｉｎａｃｃｏｒｄａｎｃｅ．

Ｋｅｙｗｏｒｄｓ：ＴａｉｈｕＬａｋｅ；ｃｈｌｏｒｏｐｈｙｌｌ－ａｃｏｎｃｅｎｔｒａｔｉｏｎ；ＥＯＳ／ＭＯＤＩＳｄａｔａ；ｒｅｍｏｔｅｓｅｎｓｉｎｇｍｏｎｉｔｏｒｉｎｇ；ｄｙｎａｍｉｃｍｏｄｅｌｉｎｇ

ＥｕｔｒｏｐｈｉｃａｔｉｏｎＡｎａｌｙｓｉｓｏｆＳｅｗａｇｅＴｒｅａｔｅｄｂｙＳｕｓｐｅｎｄｅｄＣａｒｒｉｅｒｓＢｉｏ－ｆｉｌｍＢｅｄＲｅａｃｔｏｒＫＡＮＧＹａｎ－ｈｏｎｇ，ＺＥＮＧＪｉｎｇ，ＹＵＬｏｎｇ，ＺＨＡＯＬｉａｎｇ

（ＳｃｈｏｏｌｏｆＣｈｅｍｉｓｔｒｙａｎｄＬｉｆｅＳｃｉｅｎｃｅ，ＳｈｅｎｙａｎｇＮｏｒｍａｌＵｎｉｖｅｒｓｉｔｙ，

Ｓｈｅｎｙａｎｇ１１００３４）

ＡｂｓｔｒａｃｔＩＤ：１００３－６５０４（２００７）０７－００３２－０３－ＥＡ

Ａｂｓｔｒａｃｔ：Ｗａｔｅｒｅｕｔｒｏｐｈｉｃａｔｉｏｎｗｉｔｈａｄｄｉｔｉｏｎｏｆｓｅｗａｇｅｔｒｅａｔｅｄｂｙｓｕｓｐｅｎｄｅｄｃａｒｒｉｅｒｂｉｏ－ｆｉｌｍｂｅｄｒｅａｃｔｏｒｗａｓａｎａｌｙｚｅｄ．Ｂｙｓｔｕｄｙｉｎｇｔｈｅｒｅｐｒｅｓｅｎｔａｔｉｖｅｐａｒａｍｅｔｅｒｓｉｎｃｌｕｄｉｎｇｔｏｔａｌｎｉｔｒｏｇｅｎ，ｔｏｔａｌｐｈｏｓｐｈｏｒｕｓａｎｄｃｈｌｏｒｏｐｈｙｌｌ－ａ，ｔｈｅｃｏｒｒｅｌａｔｉｏｎｂｅｔｗｅｅｎｔｅｍｐｅｒａｔｕｒｅ，ｐＨ，ｔｕｒｂｉｄｉｔｙａｎｄｃｈｌｏｒｏｐｈｙｌｌ－ａｗａｓｄｉｓｃｕｓｓｅｄ．Ｔｈｅｅｕｔｒｏｐｈｉｃａｔｉｏｎｐｒｏｃｅｓｓａｎｄｅｆｆｉｃｉｅｎｃｙｏｆｓｕｓｐｅｎｄｅｄｃａｒｒｉｅｒｂｉｏ－ｆｉｌｍｂｅｄｒｅａｃｔｏｒｗａｓｅｖａｌｕａｔｅｄａｃｃｏｒｄｉｎｇｔｏｗａｔｅｒ

ｑｕａｌｉｔｙ

ａｓｓｅｓｓｍｅｎｔ

ｓｔａｎｄａｒｄ．

Ｒｅｓｕｌｔｓ

ｉｎｄｉｃａｔｅｄ

ｔｈａｔ

ｔｏｔａｌ

ｐｈｏｓｐｈｏｒｕｓｗａｓｔｈｅｒｅｓｔｒｉｃｔｉｖｅｆａｃｔｏｒ，ａｎｄｔｈｅｗａｔｅｒｗａｓｃｌｏｓｅｔｏｍｉｄｄｌｅｅｕｔｒｏｐｈｉｃａｔｉｏｎｌｅｖｅｌ．

Ｋｅｙｗｏｒｄｓ：ｓｕｓｐｅｎｄｅｄｃａｒｒｉｅｒｂｉｏ－ｆｉｌｍｂｅｄｒｅａｃｔｏｒ；ｅｕｔｒｏｐｈｉｃａｔｉｏｎ；ＴＮ；ＴＰ；ｃｈｌｏｒｏｐｈｙｌｌ－ａ

ＰｒｅｌｉｍｉｎａｒｙＡｎａｌｙｓｉｓａｎｄＣｏｎｔｒｏｌｏｆＴＮｉｎＤａｎｊｉａｎｇｋｏｕ

Ｒｅｓｅｒｖｏｉｒ

Ａｂｓｔｒａｃｔ

ＹＩＮＭｉｎｇ１，ＳＨＩＭｉｎ－ｆａｎｇ２，ＬＩＵＣｈｅｎｇ－ｆｕ３

（１．ＳｈｉｙａｎＭｕｎｉｃｉｐａｌＥｎｖｉｒｏｎｍｅｎｔＭｏｎｉｔｏｒｉｎｇＣｅｎｔｒｅ，Ｓｈｉｙａｎ４４２０００；２．

ＨｕｂｅｉＥｎｖｉｒｏｎｍｅｎｔＭｏｎｉｔｏｒｉｎｇＣｅｎｔｒｅ，Ｗｕｈａｎ４３００７２；３．Ｈｕｂｅｉ

ＥｎｖｉｒｏｎｍｅｎｔＳｃｉｅｎｃｅＡｃａｄｅｎｕｇ，Ｗｕｈａｎ４３００７２）

ＡｂｓｔｒａｃｔＩＤ：１００３－６５０４（２００７）０７－００３５－０３－ＥＡ

Ａｂｓｔｒａｃｔ：ＲｅａｓｏｎｓｏｆｔｏｔａｌｎｉｔｒｏｇｅｎｏｖｅｒｔｈｅｃｒｉｔｅｒｉａｉｎＤａｎｊｉａｎｇｋｏｕＲｅｓｅｒｖｏｉｒｗｅｒｅｐｒｅｌｉｍｉｎａｒｙａｎａｌｙｚｅｄｂａｓｅｄｏｎａｓｓｅｓｓｍｅｎｔｏｆｙｅａｒｌｙｃｈａｎｇｅｔｅｎｄｅｎｃｙｏｆｗａｔｅｒｑｕａｌｉｔｙｉｎｒｅｃｅｎｔ２０ｙｅａｒｓａｓｗｅｌｌａｓｕｐｓｔｒｅａｍｉｎｃｏｍｉｎｇｗａｔｅｒａｎｄｓｅｗａｇｅ．Ｍｅａｓｕｒｅｓｗｅｒｅａｌｓｏｐｒｏｐｏｓｅｄｔｏｃｏｎｔｒｏｌｔｈｅｐｏｌｌｕｔｉｏｎ．

Ｋｅｙｗｏｒｄｓ：ｔｏｔａｌｎｉｔｒｏｇｅｎ；ｒｅｓｅｒｖｏｉｒ

ＤｅｖｅｌｏｐｍｅｎｔｏｆＴｅｃｈｎｉｃａｌＭａｎａｇｅｍｅｎｔＳｙｓｔｅｍｏｆＥｍｅｒｇｅｎｃｙＭｏｎｉｔｏｒｉｎｇｏｎＥｎｖｉｒｏｎｍｅｎｔａｌＰｏｌｌｕｔｉｏｎ

Ａｃｃｉｄｅｎｔ

ＷＥＮＪｉａｎ－ｈｕｉ，ＨＵＡＮＧＨｕｉ－ｌａｉ，ＳＨＩＨｕｉ，ＫＡＮＧＸｉｎｇ－ｋｕａｉ

（ＧｕｉｌｉｎＥｎｖｉｒｏｎｍｅｎｔＭｏｎｉｔｏｒｉｎｇＣｅｎｔｒｅ，Ｇｕｉｌｉｎ５４１００２）

ＡｂｓｔｒａｃｔＩＤ：１００３－６５０４（２００７）０７－００３７－０３－ＥＡ

Ａｂｓｔｒａｃｔ：Ａｔｅｃｈｎｉｃａｌｍａｎａｇｅｍｅｎｔｓｙｓｔｅｍｏｆｅｍｅｒｇｅｎｃｙｍｏｎｉｔｏｒｉｎｇｏｎｅｎｖｉｒｏｎｍｅｎｔａｌｐｏｌｌｕｔｉｏｎａｃｃｉｄｅｎｔｗａｓｄｅｖｅｌｏｐｅｄ，ｗｉｔｈｄｉｓｃｕｓｓｉｏｎｏｎｔｈｅｄｅｖｅｌｏｐｍｅｎｔｐｒｏｃｅｓｓａｎｄｂｒｉｅｆｌｙｉｎｔｒｏｄｕｃｔｉｏｎｔｏｓｙｓｔｅｍｉｃｆｕｎｃｔｉｏｎ，ｃｈａｒａｃｔｅｒｉｓｔｉｃｓａｎｄｄａｔａｂａｓｅｓｔｒｕｃｔｕｒｅａｓｗｅｌｌａｓｓｕｍｍａｒｙｏｆｅｘｐｅｒｉｅｎｃｅａｎｄｅｆｆｅｃｔ．Ｆｕｒｔｈｅｒｄｅｖｅｌｏｐｍｅｎｔａｎｄｉｍｐｒｏｖｅｍｅｎｔｗａｓａｌｓｏｐｒｏｓｐｅｃｔｅｄ．Ｋｅｙｗｏｒｄｓ：ｅｎｖｉｒｏｎｍｅｎｔａｌｐｒｏｔｅｃｔｉｏｎ；ｐｏｌｌｕｔｉｏｎａｃｃｉｄｅｎｔ；ｅｍｅｒｇｅｎｃｙｍｏｎｉｔｏｒｉｎｇ；ｃｏｍｐｕｔｅｒｓｏｆｔｗａｒｅ；ｄａｔａｂａｓｅ

ＡｎａｌｙｓｉｓｏｆＣａｌｏｒｉｃＶａｌｕｅｏｆＦｌａｍｍａｂｌｅＲｅｆｕｓｅＷＥＮＫｅ－ｊｕｎ，ＷＵＬｉ－ｐｉｎｇ，ＹＡＮＧＬｉ，ＷＡＮＧＴｉｎｇ

（ＤｅｐａｒｔｍｅｎｔｏｆＥｎｖｉｒｏｎｍｅｎｔａｌａｎｄＭｕｎｉｃｉｐａｌＥｎｇｉｎｅｅｒｉｎｇ，Ｔｉａｎｊｉｎ

ＩｎｓｔｉｔｕｔｅｏｆＵｒｂａｎＣｏｎｓｔｒｕｃｔｉｏｎ，Ｔｉａｎｊｉｎ３００３８４）

ＡｂｓｔｒａｃｔＩＤ：１００３－６５０４（２００７）０７－００４０－０３－ＥＡＡｂｓｔｒａｃｔ：

Ｆｌａｍｍａｂｌｅｒｅｆｕｓｅｉｓａｎｏｒｇａｎｉｃｗａｓｔｅｔｈａｔｕｎｅａｓｉｌｙ

ｂｉｏｄｅｇｒａｄａｂｌｅａｎｄｆｌａｍｍａｂｌｅ．Ｉｍｐａｃｔｆａｃｔｏｒｓｏｆｃａｌｏｒｉｃｖａｌｕｅｏｆｆｌａｍｍａｂｌｅｒｅｆｕｓｅｉｎｃｌｕｄｉｎｇｉｎｃｉｎｅｒａｔｉｏｎｏｆｒａｗｍａｔｅｒｉａｌｓ，ｂｒｏｋｅｎｐａｒｔｉｃｌｅｓｉｚｅａｎｄｍｏｉｓｔｕｒｅｃｏｎｔｅｎｔｗｅｒｅａｎａｌｙｚｅｄｂｙｏｒｔｈｏｇｏｎａｌｔｅｓｔ．Ｔｈｅｃｏｎｔｒｉｂｕｔｏｒｙｅｆｆｉｃｉｅｎｃｙｏｆｔｈｒｅｅｆａｃｔｏｒｓｗａｓｇａｉｎｅｄｂｙｏｒｔｈｏｇｏｎａｌａｎａｌｙｓｉｓ，ａｎｄｔｈｅｌａｒｇｅｓｔｃａｌｏｒｉｃｖａｌｕｅｃｏｍｂｉｎａｔｉｏｎｏｆｅｘｐｅｒｉｍｅｎｔａｌｍａｔｅｒｉａｌｗａｓｏｂｔａｉｎｅｄ．Ｋｅｙｗｏｒｄｓ：

ｆｌａｍｍａｂｌｅｒｅｆｕｓｅ；ｃａｌｏｒｉｃｖａｌｕｅ；ｏｒｔｈｏｇｏｎａｌｔｅｓｔ；

ｃｏｎｔｒｉｂｕｔｏｒｙｅｆｆｉｃｉｅｎｃｙ

ＨｅａｖｙＭｅｔａｌＰｏｌｌｕｔｉｏｎｉｎＶｅｇｅｔａｂｌｅＰｌａｎｔａｔｉｏｎｉｎ

Ｇｕｉｌｉｎ

ＳＨＩＨｕｉ，ＨＵＡＮＧＨｕｉ－ｌａｉ，ＷＥＮＪｉａｎ－ｈｕｉ

（ＧｕｉｌｉｎＭｕｎｉｃｉｐａｌＥｎｖｉｒｏｎｍｅｎｔＭｏｎｉｔｏｒｉｎｇＣｅｎｔｒｅ，Ｇｕｉｌｉｎ５４１００２）ＡｂｓｔｒａｃｔＩＤ：１００３－６５０４（２００７）０７－００４３－０３－ＥＡ

Ａｂｓｔｒａｃｔ：ＣｏｎｃｅｎｔｒａｔｉｏｎｏｆＣｕ，Ｐｂ，ＺｎａｎｄＣｄｉｎｉｒｒｉｇａｔｉｏｎｗａｔｅｒ，ｓｏｉｌａｎｄｖｅｇｅｔａｂｌｅｗａｓｄｅｔｅｒｍｉｎｅｄｔｏｅｖａｌｕａｔｅｈｅａｖｙｍｅｔａｌｐｏｌｌｕｔｉｏｎｉｎｖｅｇｅｔａｂｌｅｐｌａｎｔａｔｉｏｎｏｆＧｕｉｌｉｎＣｉｔｙ，ＧｕａｎｇｘｉＰｒｏｖｉｎｃｅ．ＩｔｗａｓｆｏｕｎｄｔｈａｔｉｒｒｉｇａｔｉｏｎｗａｔｅｒｑｕａｌｉｔｙｍｅｔｃｌａｓｓＩＩＩｏｆｉｒｒｉｇａｔｉｏｎｗａｔｅｒｑｕａｌｉｔｙｓｔａｎｄａｒｄ，ａｎｄｈｅａｖｙｍｅｔａｌｃｏｎｃｅｎｔｒａｔｉｏｎｉｎｓｏｉｌｗａｓｈｉｇｈｅｒｔｈａｎｂａｃｋｇｒｏｕｎｄｖａｌｕｅ，