湖泊水库水质模型ppt

第五章

湖泊水库水质模型

1.2本章内容

?湖泊水库的污染特性

?湖泊水库温度模型

?湖泊水库水质模型

?湖泊水库富营养化模型

?湖泊水库生态系统模型

?湖泊水库水质模拟通用软件介绍

?案例分析

第五章

湖泊水库水质模型

湖泊水库的污染特性

3.污染来源与途径

4?污染源

污染类型

污染物来源

工业废水

点源

城镇生活污水

固体废物处置场

矿区地表径流

外源污染物

城镇地表径流

农牧区地表径流

面源

大气降尘

大气降水

水体投饵养殖

水面娱乐活动废弃物

水土流失及土壤侵蚀

内源污染物

底泥及沉积物

?5污染物入湖途径

点源污染

降尘

河渠

流域非点源污染

湖入泊降水

地水表

库

面地下水

养殖投饵

营养源与营养负荷

6湖库营养源

（1）点源

（2）非点源

（3）内源

营养负荷

（1）地面径流的营养负荷Ijl（2）降水的营养负荷Ijp

（3）人为因素营养负荷（生活污水Ijs和工业污水Ijk）（4）湖泊水库的总营养平衡

Ij?Ijl?IjP?Ijs?Ijk

?

?7富营养化的控制因子

?湖泊富营养化定义

指湖泊水体接纳过量的氮、磷等营养性物质，使水体中藻类以及其他水生生物异常繁殖，水体透明度和溶解氧变化，造成湖泊水质恶化，加速湖泊老化，从而使湖泊生态和水功能受到阻碍和破坏。

?富营养化的控制因子

地理位置

8

人类的影响地质地形纬度地下水、农业、湖盆的形态气候矿业基底的组成深度湖底的降水风日照湖的流入的营养盐流域生态系统凹凸面积热的传氧气季湖水湖底堆积物湖外物质的流播和水的供节流动的性质入（营养盐）透明度光的透入温分层给和变和混消耗化合湖泊营养度植物的生产活动湖岸带的发育（1）限制性营养物质

9

（2）

（3）

（4）溶解氧和温度和照度

GT?Gmax?T?20T?201GT?1.8(1.066)If?KHR2.718282?eIe?If1?sIsK(e?eeHT)湖库形态

pH10湖泊水库富营养化判别标准

?判别指标

（1）物理指标

即水温、透明度、光强等，透明度最为常用。

（2）化学指标

与藻类增殖直接相关的DO，CO2，COD以及营养盐等指标。

（3）生物学指标

大致可分为藻类现存量（叶绿素）、生物指标（调查特定生物出现的状况）、多样性指数（调查群集生物的多样性）和藻类增殖的潜在能力。

11?判别标准

（1）卡尔森营养状态指数（TSI）

48TSI(SD)?10(6?lnSDln

ln2)TSI(TP)?10(6?TPln2)

TSI(chla)?10(6?2.04?0.68lnchla（2）修正的营养状态指数（ln2)TSIM）

TSIM(chla)?10(2.46?lnchla3.69?1.53lnSDln2.5)TSIM(SD)=10(2.46+ln2.5)

TSI(TP)=10(2.46+6.71+1.15lnTP

ln2.5)（3）综合营养状态指数

m2

TLI(?)??WWrijj?TLI(j)j?j?1m

?r2ijj?1

第五章

湖泊水库水质模型

湖泊水库温度模型

12.13湖泊水库水温特征

?湖泊水温受湖面以上气象条件（主要是气温与风）、湖泊容积和水深以及湖盆形态等因素的影响，呈现出具有时间与空间的变化规律，比较明显的季节性变化与垂直变化。

?一般容积大、水深大的湖泊，水温常呈垂向分层型。通常水温的垂向分布有三个层次，上层温度较高，下层温度较低，中间为过渡带，称为温跃层。

?冬季因表面水温不高，可能没有显著的温跃层。夏季的温跃层较为明显。

14?水温分层

TZ湖库中温度的竖向分布湖库中的热分层

均匀混合温度模型

15对于均匀混合型湖泊，假定水温在各个方向是均匀的，仅考虑它随时间的变化，可利用总体热量平衡模型，计算湖泊温度随时间的变化过程。

?T?t?(SR?SRb?AR?ARb?BR?E?C?HOI)tHOISRSRb

ARARb

CEBR

16分层温度模型

?垂向一维温度模型

A(z)?

??T??t?u?T????T?z?t????Z??A(z)(D?Ev)?z???S当忽略水流垂向速度及分子扩散项时

?T?t?1??A(z)?z??A(z)E?T?Sv?z???A(z)?能量平衡模型

能量平衡模型将湖库沿垂向分为

n+m个小薄层其中，

17上部

m层为混合均匀层，下部

n层为变温层，每层厚度为ΔZ。

混合层

混合层

Δz

变温层

变温层

Ek/Ep?1Ek/Ep≥?

1*18Ek??0WAsdtmEg????V(i,k)?P(m?1,k)??(i,k)??p????(m?1?i)?zi?1?

当

Ek>Ep时，水面输入的风能将转化为势能，分层不稳定，混合层水体向下扩展。

?

当

Ek>Ep时，水体形成稳定的分层，水面风能用于克服水体黏性而消失，混合层水体的厚度不增加。

第五章

湖泊水库水质模型

湖泊水库水质模型

19.湖库完全混合箱式模型——沃伦威德尔模型

20概述

创始：沃伦威德尔（R.A.Vollenweider）在

20世纪

70年代初期研究北美大湖时提出。

适用：停留时间很长，水质基本处于稳定状态的湖泊水库。

假定：湖泊中某种营养物的浓度随时间的变化率，是输入、输出和在湖泊内沉积的该种营养物量的函数。

不足：不能描述发生在湖泊内的物理、化学和生物过程，同时也不考虑湖泊和水库的热分层，是只考虑其输入—产出关系的模型。

?21?模型

VdCdt?Ic?sCV?QC

引入冲刷速度常数

r（令

r＝

Q/V），则得到

dC?Ic?sC?rC

dtV在给定初始条件，当

t＝

0，C＝

C0

时，求得上式的解析解为C?IcV(s?r)C0?IcV(s?r)?V(s?r)exp[?(s?r)t]

22

在湖泊、水库的出流、入流流量及营养物质输入稳定的情况下，当

t→∞时，可以得到营养物质的平衡浓度

Cp：

CIcp?(r?s)V

如进一步令

t1Vw??且

V?Ah

rQs，则湖泊、水库中的

营养物质平衡浓度变为：

Ccp?Lsh?htLIwcc?As23?案例分析

已知湖泊的容积

V=1.0×107m3/a，湖泊内

CODCr的本底浓度

C0=1.5mg/L，河流中

COD浓度为C1=3mg/L，COD在湖泊中的沉积速度常数

s=0.08a-1。试求湖泊中的

COD平衡浓度，及达到平衡浓度的99%所需的时间。

解答：根据题目，得到

CC?1???V?s?r?C0?1???s?r?t?I?epc?C?1??C???t??1Cp1?C1?Icp??s?rlnV?s?r?C??s?rln0V?s?r?C0?II?1cc24

根据题意已知：V=1.0×107m3，s=0.08a-1，r=Q/V=5a-1

，

C0=1.5g/m3，Ic=0.5×108×3=1.5×108g/a。

当

C/Cp=0.99时：

8

t??1?0.99?1??1.5?100.08?5ln1.0?107??0.08?5??1.5?1.5?108

??15.08ln0.02033?0.77?a?

8

C1.5?103p?

?0.08?5??107?2.95?g/m?

根据公式，计算结果为：达到

COD平衡浓度的

99%约需

0.77a；平衡浓度值为

2.95g/m3。

湖库完全混合箱式模型——吉柯奈尔—狄龙模型

25模型

引入滞留系数

Rc：

dCIc(1?Rc)dt?V?rC如给定初始条件

t＝

0，C＝

C0，得到上式的解析解：

C?Ic(1?Rc)Ic(1?rV?[CRc)0?rV]exp(?rt)若湖库得出流、入流及污染物的输入都比较稳定，当

t→∞时，可以达到营养物质的平衡浓度

Cp

：

CIc(1?Rc)Lc(1?Rc)p?rV?rh?

26可根据湖库的入流、出流近似计算出滞留系数。

?mq0jC0jRc?1?j?1?nqikCikk?1湖库完全混合箱式模型——分层箱式模型

27概述

1975年，斯诺得格拉斯（Snodgrass）等提出了一个分层的箱式模型，用以近似描述水质分层状况。分层箱式模型把上层和下层各视为完全混合模型，在上、下层之间紊流扩散的传递作用。分层箱式模型分为夏季模型和冬季模型，夏季模型考虑上、下分层现象，冬季模型则考虑上、下层之间的循环作用。模拟包含的水质组分为正磷酸盐（Po）和偏磷酸盐（Pp）的变化规律。

?

28?概化图

?Q(poj?ppj)Q(poe?ppe)?Qj(poj?ppj)Q(poe?ppe)KthpoeppeKVpoethZAthPphthZAthPpepethpoKKththArZAthPohZthPoeththpohpShAsPpephVShAsPphSAsPp?模型

（1）夏季分层模型

29

对表层正磷酸盐PoeVdPoedt??Q?PVkthejPoj?QPoeeePoe?ZAth(Poh?Poe)th对表层偏磷酸盐PpeVdPpekthedt??QjPpj?QPpe?SeAthPpe?PVeePoe?ZAth(Pph?Ppe)th对下层正磷酸盐PohVdPohkthhdt?rhVhPph?ZAth(Poe?Poh)th对下层偏磷酸盐PphVdPphkdt?SVthh?SeAthPpehAsPph?rhhPph?ZAth(Ppe?Pph)th30（2）冬季分层模型

对全湖的正磷酸盐Po：

VdPodt?QjPoj?QPo?PeuVeuPo?rVPp

对下层偏磷酸盐Pp：

VdPpdt?QjPpj?QPp?PeuVeuPo?rVPp?SAsPp

夏季的分层模型和冬季的循环模型可以用秋季或春季的“翻池”过程形成的完全混合状态作为初始条件，此时

31PPoeVe?PohVho?VPPpeVe?PphVhp?V湖库分层水质模型——卡拉乌舍夫湖泊水库水质扩散模型

32模型

?Cr?t????MQp?1?Cr?2Cr

r??H???r??r?Mr?r2当为稳定排放，且边界条件取距排放口充分远的某点

r0处的现状值

Cr0，上式求解得：

QpC?r??HMrr?Cp?(Cp?Cr0)??r?0?1/3M?H2/3d?2r?f?uh?0g?πH?2??u??湖库分层水质模型——易降解物质简化的水质模型

33

模型

将卡拉乌舍夫湖泊水库水质扩散模型中扩散项忽略掉，

并考虑污染物的降解作用，这样既可得到稳态条件下污染物在湖库中推流和生化降解共同作用下的基本方程：QdCrpdr??K1CrH?r?

当边界条件取

r=0时，Cr=Cr0（为排出口浓度），则其解析解为：

?K1?Hr2Cr?Cr0exp???172800Qp?????当考察湖库的水质指标是溶解氧时，并只考虑

BOD的耗氧因素与大气复氧因素，可推导出湖库的氧亏方程：

解析解：

dDQp?(K1L?K2D)H?rdr2?2??2?????K1L0K1?HrK2?HrK2?HrD???exp???D0exp???exp????????????K2?K1?2Q2Q2Qppp????????34

第五章

湖泊水库水质模型

湖泊水库富营养化模型

35.经验模型

36

分类

（1）单一营养物质负荷模型

（2）藻类生物量与营养物质负荷量之间的相关模型。?

37?

一般模型

根据物质平衡原理，

对于完全混合型湖泊、水库，

以磷负荷为研究对象，

磷负荷的收支平衡方程为：

dWp

dt?WQpi?kpWp?VWp解析解：

WW?pi?W?pi?1?1?W??1??pt??p0??exp???kT?kp??T?kp????Tp??t??经验模型——沃伦威德（Vollenweider）模型

38

模型

P?Wp?L

V??QHk?H

?k??PV??PT

湖泊磷的临界负荷量计算公式：

L10q?H?cr?s??1???qs??L?WPiA?

经验模型——狄龙（Dillon）模型

39

模型

RpL?kkp??wP?L(1?RL)H?wR?1?QPL??QiPi?

生态模型

40φrgensen（1976）根据丹麦Glumsoe湖的研究建立了湖泊

该模型主要包括以下几个主要成分：

（1）三个营养级，即浮游植物、浮游动物和鱼类

（2）浮游植物生物量动态

（3）水—沉积物之间营养物质的交换

（4）氮的反硝化作用

（5）温度、光照对藻类生长率的影响

（6）河流氮的输入情况等

J富营养化生态模型。

生态动力学模型

41

模型

?C?t?u?C?x?v?C?y?w?C?z???x(E?C??Cx?x)??y(Ey?y)???z(E?Cz?z)?F(C)?S?

42?功能

对营养盐循环过程的模拟，对藻类生长机理的模拟，对影响藻类生长相关因子的定量描述。

?模拟对象

浮游植物的生长与死亡，营养盐的循环过程，DO的动力过程

NO?310

P

10

PO3?43

N

8

13

11

C

11

12

NO?23

Chla

ORG-P

PHTY

7

NH?413

NH9

38

12

5

ORG-N

7

4

O1

DO

2

CBOD

2

2

13

6

CO7

水体

底部沉积物

湖库生态系统中生态变量的循环

第五章

湖泊水库水质模型

湖泊水库生态系统模型

43.44概念模型

以磷为核心的湖泊水库生态系统模型包括下述水质项目：

藻类、浮游动物、有机氮、无机磷、有机磷、氨氮、硝酸盐、亚硝酸盐、含碳有机物的生化需氧量、溶解氧、总溶解固体和悬浮物等

12个水质项目。

降解BODBOD耗氧DO水面复氧磷氨氮耗氧氮光藻1合类作呼PO3-底泥4-P有机磷用吸需氧沉--产耗厌氧释放氨氮NO-N淀2NO3-N有机氮氧氧沉淀物吸氧化氧化藻类释放藻类物吸吸收藻类收收水解藻类吸收藻类释放藻类吸收厌氧

沉淀磷浮游动光沉淀反硝浮游物捕食衰1氮释化动物减放沉淀氮释放浮游动物1悬浮泥渣一般数学表达式

45

12个水质项目都可以用下述偏微分方程表示

?C?t?(V?V?C11??C?Sint1s)?z?A??z??ADz?z???A?A(qinCin?qoutCout)上述

系统模拟

46

藻类（浮游生物）生物量

浮游动物

Z

磷

溶解态的无机磷

P1游离态的有机磷

P2沉淀态的磷

P3

CAdCAdt??CA?(??CgZ)CAdZdt??zZ?(?z?Cz)ZdP1dt???CA(App)?(I3P3?I1P1)?I2P2dP2dt??CAApp??zZApz?(I4P2?I2P2)dP3dt?I4P2?I3P3???

47