第三章水环境化学-第四节水质模型

1、n 水质模型，是一个用于描述物质在水环境中的混合、水质模型，是一个用于描述物质在水环境中的混合、迁移、扩散和转化过程（包括物理、化学、生物作用过迁移、扩散和转化过程（包括物理、化学、生物作用过程）的数学方程（或方程组）程）的数学方程（或方程组） .第四节第四节 水质模型水质模型n 水质模型的基本原理是质量守恒原理；建立水质模水质模型的基本原理是质量守恒原理；建立水质模型的目的是用来描述污染物数量与水环境影响因素之间型的目的是用来描述污染物数量与水环境影响因素之间的定量关系，从而为水质分析、预测和水环境管理提供的定量关系，从而为水质分析、预测和水环境管理提供基础的量化依据。基础的量化依据。 n

2、本节讨论的水质模型主要是：氧平衡模型、湖泊富本节讨论的水质模型主要是：氧平衡模型、湖泊富营养化模型和有毒有机污染物归趋模型。营养化模型和有毒有机污染物归趋模型。1. 1. Streeter-PhelpsStreeter-Phelps（S-PS-P）模型（河流水质自净模型）模型（河流水质自净模型）nS-PS-P模型的建立基于两项假设：模型的建立基于两项假设：（1 1）只考虑好氧微生物参加的有机物降解反应，并）只考虑好氧微生物参加的有机物降解反应，并认为该反应为一级反应。认为该反应为一级反应。（2 2）河流中的耗氧只是有机物降解反应引起的。有）河流中的耗氧只是有机物降解反应引起的。有机物的降解反应

3、速率与河水中溶解氧机物的降解反应速率与河水中溶解氧( (DO)DO)的减少速的减少速率相同，大气中的氧进入水体的复氧速率与河水中率相同，大气中的氧进入水体的复氧速率与河水中的亏氧量的亏氧量 D D 成正比。成正比。一、氧平衡模型一、氧平衡模型式中式中:L河水中的河水中的BOD值，值，mg/L；D河水亏氧值，河水亏氧值，mg/L,是饱和溶解氧浓度是饱和溶解氧浓度Cs （mg/L）与实际溶解氧浓度与实际溶解氧浓度C（mg/L）的差值的差值D=CS-C；k1河水河水耗氧速度常数，耗氧速度常数，1d；k2河水复氧速度常数，河水复氧速度常数，1d；uu河水平均流速河水平均流速km/d;x-顺河水流动方向

4、的纵向距离顺河水流动方向的纵向距离km。S-PS-P模型的基本方程为：模型的基本方程为：0,0,00 xCCxLL)()(/2101/0/02121uxkuxkuxkssuxkeekkLkeCCCCeLL当边界条件当边界条件S-PS-P模型基本方程模型基本方程的解析解为的解析解为:05101520250246810012345678X kmL mg/LL mg/LDOmg/LDOmg/LS-P 模型的临界点和临界点氧浓度模型的临界点和临界点氧浓度n 以以S-P解析解中溶解氧解析解中溶解氧C对距离作图可得一条下垂曲线，对距离作图可得一条下垂曲线，称为氧垂曲线。称为氧垂曲线。n 氧垂曲线最低点所对

5、应的溶解氧，称为极限溶解氧氧垂曲线最低点所对应的溶解氧，称为极限溶解氧Cc，出出现现Cc的距离称为极限距离的距离称为极限距离xc；n 在极限距离处在极限距离处xc，溶解氧变化率为零，由，溶解氧变化率为零，由S-P方程可得：方程可得：极限距离：极限距离：极限溶解氧：极限溶解氧：（DC为极限氧亏）为极限氧亏）2 2托马斯托马斯(Thomas)模型模型n对于一维静态河流，在对于一维静态河流，在SPSP模型的基础上考虑沉淀、絮模型的基础上考虑沉淀、絮凝、冲刷和再悬浮过程对凝、冲刷和再悬浮过程对BODBOD变化的影响，引入了变化的影响，引入了BODBOD沉沉浮系数浮系数k k3 3DkLkdxdDuLk

6、kdxdLu2131)(n 当边界条件当边界条件0,0,00 xCCxLL求取解析解可得：求取解析解可得：二、水体富营养化预测模型二、水体富营养化预测模型1 1水体的富营养化问题水体的富营养化问题n 水的停留时间较长（可达数月至数年），属于缓流水域，水的停留时间较长（可达数月至数年），属于缓流水域，其中的化学和生物学过程保持一个比较稳定的状态。其中的化学和生物学过程保持一个比较稳定的状态。n 进入湖泊和水库中的营养物质在其中容易积累，致使水进入湖泊和水库中的营养物质在其中容易积累，致使水质发生富营养化。质发生富营养化。n 在水深较大的湖、库中，水温和水质是竖向分层的。在水深较大的湖、库中，水温

7、和水质是竖向分层的。湖泊湖泊(水库水库)的水质特征：的水质特征：n 指在人类活动的影响下，氮、磷等生物所需营养物质大指在人类活动的影响下，氮、磷等生物所需营养物质大量进入湖泊、水库、海湾等缓流水体，引起藻类及其他浮量进入湖泊、水库、海湾等缓流水体，引起藻类及其他浮游生物迅速繁殖，水体溶解氧量下降，水质恶化，引起鱼游生物迅速繁殖，水体溶解氧量下降，水质恶化，引起鱼类及其他生物大量死亡的现象。类及其他生物大量死亡的现象。 雨水：雨水： 雨水中硝酸盐的含量在雨水中硝酸盐的含量在 0.16 -1.06 0.16 -1.06 mg/Lmg/L间，氨氮间，氨氮含量在含量在 0.04 - 1.70 0.04

8、 - 1.70 mg/Lmg/L间，大面积湖泊和水库从雨水间，大面积湖泊和水库从雨水受纳氮数量相当可观。受纳氮数量相当可观。 农业排水：农业排水： 土壤中的氮、磷通过地表水径流被引入到湖土壤中的氮、磷通过地表水径流被引入到湖泊、水库。泊、水库。 水中营养物质的来源水中营养物质的来源城市污水：城市污水：包括：排泄物、食品污物，含磷洗涤剂。污水厂包括：排泄物、食品污物，含磷洗涤剂。污水厂通过厌氧处理污泥的方法可去除通过厌氧处理污泥的方法可去除20205050的氮和大部分的磷的氮和大部分的磷，但在污水处理中也使用到多种含氮、磷的化学试剂，如氯，但在污水处理中也使用到多种含氮、磷的化学试剂，如氯胺、有

9、机絮凝剂、无机助絮凝剂、多聚磷酸钠等。胺、有机絮凝剂、无机助絮凝剂、多聚磷酸钠等。其他来源：其他来源：包括工业废水等。包括工业废水等。湖泊水质模型的类型：湖泊水质模型的类型：n 湖泊水质模型可划分为：多元相关模型；输入输出湖泊水质模型可划分为：多元相关模型；输入输出模型；富营养化预测模型和扩散模型，这里仅讨论富模型；富营养化预测模型和扩散模型，这里仅讨论富营养化预测模型。营养化预测模型。2. 富营养化预测模型富营养化预测模型n 对于停留时间很长、水质基本处于稳定状态的中小对于停留时间很长、水质基本处于稳定状态的中小型湖泊和水库，可视为一个均匀混合的水体。型湖泊和水库，可视为一个均匀混合的水体。

10、n 沃兰伟德假定，湖泊中某种营养物的浓度随时间的沃兰伟德假定，湖泊中某种营养物的浓度随时间的变化率，是输入、输出和在湖泊内沉积的该种营养物量变化率，是输入、输出和在湖泊内沉积的该种营养物量的函数，用质量平衡方程表示就是：的函数，用质量平衡方程表示就是：cVqcIdtdcVPP)(cPVIdtdcPWP)(简化可得：简化可得：湖泊富营湖泊富营养物质的养物质的变化变化=单位时间输单位时间输入湖泊营养入湖泊营养物质的量物质的量单位时间输单位时间输出湖泊营养出湖泊营养物质的量物质的量单位时间单位时间营养物质营养物质沉积的量沉积的量即：即：式中：式中：c - c - 湖水平均总磷浓度湖水平均总磷浓度 m

11、g/L, mg/L, I IP P 输入湖泊磷的浓度输入湖泊磷的浓度 g/dg/d P PW W - - 水力冲刷系数水力冲刷系数 P PW W = q / V (d= q / V (d-1-1) ) q - q - 出湖河道流量出湖河道流量 m m3 3/d, V/d, V 湖泊容积湖泊容积 m m3 3 P P - - 磷的沉降速率常数磷的沉降速率常数 d d-1-1 t - t - 河水入湖时间河水入湖时间 d d令边界条件令边界条件: t = 0, C = C0 ； t=t, C = C对上式积分可得对上式积分可得:n 当时间足够长，利用上面的积分解可求取湖泊、当时间足够长，利用上面的

12、积分解可求取湖泊、水库中污染物水库中污染物(营养物营养物)的平衡浓度的平衡浓度.n 在水质分析和水质预测中在水质分析和水质预测中,利用该式可求出湖泊、利用该式可求出湖泊、水库中污染物达到一指定浓度水库中污染物达到一指定浓度CT所需时间所需时间t0。表征环境特征：取决于实际水环境。表征环境特征：取决于实际水环境。( (水流量水流量, ,流速流速, ,pH,pH,水温水温, ,风速风速, ,细菌数细菌数, ,光强等）光强等）三、有毒有机物的归趋模型三、有毒有机物的归趋模型n 通过研究水环境中各种有机毒物的动力学模型，可通过研究水环境中各种有机毒物的动力学模型，可以预测污染物在环境中浓度的时空分布，

13、以及通过各种以预测污染物在环境中浓度的时空分布，以及通过各种迁移转化过程后的最终归趋，这对合理使用有机有毒物迁移转化过程后的最终归趋，这对合理使用有机有毒物质意义重大。质意义重大。有机污染物迁移转化的动力学机理有机污染物迁移转化的动力学机理表征化合物固有性质：可由实验室测得。表征化合物固有性质：可由实验室测得。（溶解度（溶解度, ,蒸汽蒸汽 压，辛醇水分配系数等）压，辛醇水分配系数等）模型中的水模型中的水质参数：质参数：化合物迁移转化过程：化合物迁移转化过程：n负载过程（输入过程）负载过程（输入过程）来源：污水人为排放，来源：污水人为排放，大气沉降，陆地径流大气沉降，陆地径流等将有机毒物引入水

14、等将有机毒物引入水体。体。n迁移过程迁移过程 沉淀沉淀- -溶解作用：有机或无机物溶解度对其迁移转化可溶解作用：有机或无机物溶解度对其迁移转化可利用性的影响，将改变迁移的速率。利用性的影响，将改变迁移的速率。 对流作用：水流作用使污染物进入（排出）特定水生对流作用：水流作用使污染物进入（排出）特定水生生态系统。生态系统。 挥发作用：大气、水中浓度分布。挥发作用：大气、水中浓度分布。 沉积作用：底部沉积物的吸附沉积作用：底部沉积物的吸附- -解吸。解吸。n 形态过程形态过程 酸碱平衡：酸碱平衡：PHPH值、有机酸碱分数值、有机酸碱分数( (离子、分子离子、分子) )，挥发，挥发作用等。作用等。

15、吸着吸着- -解吸平衡：有机物滞留在悬浮颗粒物或沉积物上解吸平衡：有机物滞留在悬浮颗粒物或沉积物上的过程最终将影响其归趋。的过程最终将影响其归趋。n转化过程转化过程 生物降解：微生物代谢将改变污染物和它们的毒性。生物降解：微生物代谢将改变污染物和它们的毒性。 光解作用：破坏有毒有机物分子的结构。光解作用：破坏有毒有机物分子的结构。 水解作用：使污染物分子变成简单分子，低毒或无毒化水解作用：使污染物分子变成简单分子，低毒或无毒化合物。合物。 氧化还原：微生物催化氧化氧化还原：微生物催化氧化, ,光催化氧化光催化氧化, ,均将改变有机均将改变有机分子的结构。分子的结构。n 生物积累过程生物积累过程

16、 生物浓缩：通过可能的生物浓缩手段（如鱼腮吸附），生物浓缩：通过可能的生物浓缩手段（如鱼腮吸附）， 摄取有机物进入生物体。摄取有机物进入生物体。 生物放大：高营养级生物以低营养级生物为食物，使生生物放大：高营养级生物以低营养级生物为食物，使生 物体中有机毒物的浓度随营养级的提高而逐步增大。物体中有机毒物的浓度随营养级的提高而逐步增大。污染物归趋模式研究思路污染物归趋模式研究思路n从研究单个的主要迁移转化过程着手，建立单个过程从研究单个的主要迁移转化过程着手，建立单个过程的模式，作为整体模式的基础。的模式，作为整体模式的基础。n假设某化合物从水环境消失的总速率，是各个单过程假设某化合物从水环境消

17、失的总速率，是各个单过程消失速率之和，每个单过程的速率均是一级反应过程，消失速率之和，每个单过程的速率均是一级反应过程，总速率亦同。总速率亦同。 n模式中环境具有固定特征，环境参数不变，速率方程为模式中环境具有固定特征，环境参数不变，速率方程为准一级方程，并假定有机污染物不会改变环境参数。准一级方程，并假定有机污染物不会改变环境参数。n假定吸着过程不是瞬间能完成的，但其速率远快于挥发假定吸着过程不是瞬间能完成的，但其速率远快于挥发和其他过程转化速率。因此，模式不适用于污染源近地和其他过程转化速率。因此，模式不适用于污染源近地距离，只适用于长时间大范围。距离，只适用于长时间大范围。模型建立步骤：

18、模型建立步骤：1.1.计算有机物因转化和挥发从水环境中消失速率计算有机物因转化和挥发从水环境中消失速率CEKRRtiiiKi表示表示i过程速率常数，过程速率常数，Ei表示表示i过程环境参数过程环境参数KT=Kvm (挥发挥发)+Kb(生物降解生物降解)+KP (光光 降解降解)+Kh (水解水解)CKCKRTiT当当Ei固定时，固定时，n有机物在颗粒物上的吸着会降低有机物在水中的浓度，有机物在颗粒物上的吸着会降低有机物在水中的浓度，吸着也会发生转化（如微生物转化代谢），但在这里吸着也会发生转化（如微生物转化代谢），但在这里不考虑转化过程或转化很慢（比溶液中慢），并且吸不考虑转化过程或转化很慢（比溶液中慢），并且吸着过程具有可逆性。着过程具有可逆性。2.2.吸着过程对有机物消失的影响吸着过程对有机物消失的影响WSPCCKWPSTCCCCn 当有机物含量很低时，它在水和颗粒物之间的分配往当有机物含量很低时，它在水和颗粒物之间的分配往往可以用分配系数往可以用分配系数(