《目前世界上主要的水模型介绍》由会员分享，可在线阅读，更多相关《目前世界上主要的水模型介绍（14页珍藏版）》请在金锄头文库上搜索。
1、5.3CEQUALICM 模型及水质富营养化模式5.3.1CEQUALICM 模型简介 CEQUAL ICM由美国陆军工程兵团水体试验基地的Carl F.Cerco和Thomas Cole等人开发， ICM代表集成网格模型，该模型的建立最初是为了应用于美国弗吉尼亚的切萨皮克湾（Chesapeake Bay） ，它能模拟一维、二维、三维水体结构，它能够模拟大量的水质变量，如：不同种类藻、不同形态碳、不同形态氮、不同形态磷、不同形态硅、化学需氧量、溶解氧、盐度、温度、金属等，对于这些状态变量可以根据自己的需求进行开关设置。但它本身没有水动力模块，所以必须从别的模型中获得流量、扩散系数和蓄水量等信息
2、。在指定观测资料和子程序的基础上，能够模拟计算底质水界面的氧和营养盐的转化通量。如果在计算过程中计算机突然中断或发生其它类似的情况，由于程序中设置了热启文件重新启动计算机后可以继续计算，有效避免了重新计算的发生。模型对于输入输出文件没有固定时间步长的限制，可以根据自己的实际情况任意设定时间步长，这是该模型的又一大优势。模型具体结构分布情况见图27 187，模型由主程序、输入输出文件和子程序组成，在处理大量输入输出文件的时候，主程序和子程序根据各自功能都能够执行读入写出的任务，模型的主程序包括3个基本的功能：对于模型运行的输入输出文件能够制定详细的说明；3维质量平衡方程的解法；处理指定的期望输入
3、输出文件。在大部分应用中它与美国陆军工程兵团的另一个水动力模型CH3DWES(曲线网格水动力三维模型)合用。它是目前世界上发展程度最高的三维模型之一。CEQUAL ICM模型以浮游植物和水生生物的生长动力学为核心，以C:N:P这三个主要元素的比例反映浮游植物和水生生物与水体环境中营养盐之间的竞图 27CE-QUAL-ICM 模型结构图Fig.27 Model subroutines and files争转化关系。模型不仅考虑了浮游植物的三种藻类（蓝藻、绿藻和硅藻）以及用不同的动力学参数、半饱和常数、新陈代谢速率等影响因子加以区别，还考虑了有机营养盐在矿化过程中根据降解速率的不同分为难分解（RE
4、FRACTORY）的营养盐、易分解（LABILE）的营养盐和溶解（DISSOLVED）的营养盐。浮游植物和水生生物生长动力过程在富营养化作用中起着非常重要的作用，影响着其它所有系统。图28反映了浮游植物和水生生物动力学变量作用关系 188。Phytop NO3SiNH3NH+4Phytop./PeriphytonLightPO4C:N:P）植 物 ） （： 沉 降 速 率 （ 仅 对 浮 游）： 死 亡 速 率 （ ）： 生 长 速 率 （ ）的 浓 度 （： 浮 游 植 物 或 水 生 生 物1s 1dg 3ppsdg mt TRTLCCRP和Periphyton 分别代表浮游植物和水生生物
5、，在光照的条件下，吸收氮磷等营养盐后，进行光合作用而释放出氧气，浮游植物和水生生物衰老死亡后，一部分在水体中水解，一部分则逐渐沉到水底底泥内被矿化。在不同单元水体体积 中，浮游植物、水生生物的生长速率和死亡速率、jV图 28浮游植物与水生生物的生长动力学变量关系Fig.28 Phytoplankton and Periphyton kinetics沉降速率之间的关系是不同的，可用式（73）表示：jsjpjjk PkDGS414（73）式中：浮游植物动力变量的反应关系，mg carbon/L-day。jkS4： 浮游植物群落，mg carbon/L。jP：生长速率常数，day -1。j1pG：死
6、亡速率和呼吸速率常数，day -1。jD：沉降速率常数，day -1。j4sk：不同的水体单元号。5.3.2水质平衡方程控制质量平衡方程对于每一个水质状态变量都可以用下式表达：cyxzyxyxxy yxxy HSmCAzCmHACAmt x（74）式中： 各水质状态变量的浓度。： 分别为在曲线 坐标下 方向的速度分量。, zyx,：分别为 方向上的紊流扩散系数。zyxAzyx,： 每个单位体积的内外源汇项。cS: 水体深度。H：水平曲线坐标 方向上的比例因子。yxm, yx,式（74）中 这三项表示平流传输过程，CmzHCyxxy 和、这三项表示扩散传输过程， AymAxHAx zyxxy 和
7、，上述六项的物理传输过程十分相似，因此，数值解法几乎也是一样的。式（74）中 表示每个水质变量的水动力过程和外部负荷。目前这个模型cyxS求解公式（74）时使用了从物理传输项中减少动力项的分步程序来完成。 pxpm cyxzyxyxxy yxxy HSmCAmzCHACAt （75）cktS（76）kyxpx tmt CHHCt yxy（77）从公式（75）和公式（76）看，源汇项已经被分为关联流入流出的物理源汇项和动力源汇项，如果这些状态变量在一定的水体深度处与水体传输域的散度相关，对于物理传输步骤来说，动力传输步骤在同一个对应的水深处已经完成了。能够进一步将反应过程和内部源汇项分离开来的动
8、力方程式（76）可以消除水深和比例因子的影响。(78)RKCtk式中：为动力速率， 代表内部源汇项，在 和 已知的条件下，物理传输过KR程和动力方程式的解法都是比较精确的。5.3.3富营养化动力方程 187,189富营养化模型的主要问题是依据藻类和水体中溶解氧浓度对碳初级生产力的影响，初级生产力为水生态系统提供了能量来源，但是初级生产能力旺盛会对水体带来有害的作用，在水体中分解、耗氧、沉降。因此，溶解氧成为衡量水生态系统健康状态的一个重要指标。5.3.3.1藻类生物量循环过程藻类在碳和营养盐循环过程中起着非常重要的作用，藻类生物量用含碳量来统一计算，然后通过比例关系，计入到叶绿素中。为了量化藻
9、类生物量对氮磷浓度的影响，需要给定藻类生物体中碳、氮、磷的比值。由于不同藻类的特性不同，控制方程中将藻类分为三大类：蓝藻（Cyanobacteria） 、硅藻（Algae diatoms）和绿藻（Algae greens） ，这三种藻类主要是依据每个种类自身的不同特征和各自在生态系统中的功能来划分的。蓝藻，通常称为蓝绿藻，其主要特征是在含盐水中生长受到限制并在淡水水体中异常繁殖形成水华。蓝藻有时还会上浮，被捕食的压力较小。硅藻在光照条件下，吸收氮磷等营养元素进行光合作用释放氧气的时候，需要硅元素来形成细胞壁，硅藻具有较大的沉降速率，春季，硅藻沉降到底泥中，可以成为底泥耗氧的一个重要碳源。绿藻在
10、分类时是指不同于蓝藻和硅藻特征的另一类藻，绿藻的沉降速率介于蓝藻和硅藻之间，其捕食压力大于蓝藻。藻类生物量循环过程中的源汇项为：藻类生长、基础新陈代谢、捕食压力和沉降四个过程。具体公式如下：xzxxt BWSPRBM（79）式中：藻类生物量，用碳来表示（gm C m-3） 。xB：藻类生长速率（day -1） 。P：基础新陈代谢速率（day -1） 。M：藻类被捕食速率（day -1） 。R：藻类死亡沉降速率（m day-1） 。WS：分别代表三种藻类。g,dcx：垂直坐标（m）z模型中，浮游植物生长速率主要受有效营养盐量、光照强度和周围环境的温度三个因素的影响，其具体关系可用乘法公式表示如下
11、：(80)TINPMffx式中：最佳条件下的生长速率（day -1） 。x：营养盐浓度的影响系数（ ） 。Nf 1f0：光照强度的影响系数（ ） 。I：环境温度的影响系数（ ） 。Tf f蓝藻在一定的盐度水体中将会遭受快速的死亡，这个影响已被引起关注，盐分对蓝藻生长的影响可用下式表示：(81)STINPMffc式中：盐分对蓝藻生长速率的影响系数（ ） 。Sf 1f0碳、氮、磷是藻类生长最为基本的营养元素，硅藻生长还需要硅元素，无机碳在水体中十分充足，所以在模型中不予考虑它的影响限制，营养盐对藻类的生长可用莫洛方程（Monod Kinetic）表示，在莫洛方程中，藻类生长速率由可利用的低浓度营养
12、盐决定，与高浓度营养盐浓度无关，在公式中半饱和浓度是一个关键的参数，半饱和浓度为生长速率是最大生长速率一半时所对应的可利用的营养盐浓度，但并不是线性关系。 “Liebig”的最小定律也表明生长速率由最小的营养盐浓度决定。对于蓝藻和绿藻的营养盐浓度影响系数可用下式表示： dpx,NOHnxmiuf 434PKKN（82）式中：氨氮的浓度（gm N m-3） 。4H：硝态氮浓度（gm N m-3） 。3NO：藻类生长吸收氮的半饱和常数（gm N m-3） 。nxK：溶解态磷酸盐浓度（gm P m-3） 。d4P：藻类生长吸收磷的半饱和常数（gm P m-3） 。pH硅藻的生长需要氮、磷、硅营养元素
13、，对硅藻生长的营养盐限制因子采用公式（82 ）和公式（83）的最小值来确定：dsfSAKHN（83）式中：溶解的可利用硅浓度（gm Si m-3） 。dSA：硅藻吸收硅元素的半饱和常数（gm Si m-3） 。sKH光照是浮游植物和水生生物进行光合作用的必要条件，藻类的生长是随着光强的增加而增加直到达到最佳光强，超过最佳光强后，会对初级生产力产生反作用，藻类的生长会随着光强的增加而下降，这一现象可用斯蒂尔方程（Steele s equation）描述如下：s1efII（84）式中：光照强度（Langleys day-1） 。I：最佳光照强度（Langleys day-1） 。s斯蒂尔方程描述了
14、一个空间点的瞬时限光过程。假设光强随着水深的增加而呈现指数衰减，斯蒂尔方程可变为：tbezes72.fKFDI（85）ZIes0b（86）ZDKFIes0t（87） 式中：白昼时间比例（ ） 。FD10FD：全辐射衰减系数（m -1） 。esK：模型分段计算厚度（m） 。z：水面日照强度（Langleys day-1） 。0I：从水面到水下模型计算部分顶部的距离（m）Z温度也是浮游植物和水生生物进行生长的关键参数之一，藻类的生长是随着温度的增加而增加直到达到最佳温度为止，超过最佳温度后，藻类的生长会随着温度的增加而下降，温度对藻类生长的影响类似于高斯概率曲线，如下式：TmxwhenT2T-mx
15、Kg-1ef（88）式中：:温度（） 。T：为藻类生长的最佳温度（） 。mx：当 时温度对藻类生长的影响系数。1gKxT：当 时温度对藻类生长的影响系数。2碳循环过程如图29所示，三种类型的藻类均具有固碳能力，藻类通过呼吸作用产生二氧化碳而失去了一部分碳，还有一部分碳通过藻类的死亡和掠夺过程转化为溶解态有机碳和颗粒有机碳。根据有机物质分解速率的不同分为易溶有机碳和难溶有机碳颗粒。有机碳的一部分通过水解作用转化为溶解态有机碳，从系统中通过异氧呼吸作用直接将产出的溶解态有机碳通过藻类或者水解作用而消耗，最后，颗粒碳随藻类一起沉积在底泥中。蓝绿藻 硅藻 鞭毛藻及其它藻溶解态有机碳易溶有机碳难溶有机碳颗粒呼吸作用沉降作用沉降作用5.3.3.2氮循环过程藻类和浮游植物在生长期间吸收氨氮和硝态氮，在呼吸和死亡分解过程中释放出有机氮和氨氮，部分有机氮颗粒水解成溶解态有机氮，其余的沉降到底泥中，溶解态有机氮矿化为氨氮。在有氧条件下，部分氨氮经过硝化作用被氧化为硝态氮，在缺氧条件下，硝态氮经过反硝化作用被还原为氮气而溢出，底泥颗粒吸
