MIKE21水质材料

1、MIKE 21水质培训教程 王颖 DHI China Office number 021 6417 7291 1 ECO Lab 简介 ECOLab是DHI在传统的水质模型概念发展起来的全新的水质和生态模拟工具。 ECOLab软件开发的理念和方法非常先进，用户不仅可以修改模型参数，更重要的是可 以修改模型核心程序、甚至编写新程序，然后ECOLab将其与MIKE 11/21/3的HD、 AD集成计算。 DHI已经将大部分传统的水质模块转换成 ECOLab通用模板，供用户调用或修改使用， 包括： ?水质模块 ?富营养化模块 ?重金属模块 1.1 应用领域 ?河流、湿地、湖泊、水库、河口、海岸和海洋

2、 ?各生态系统反应的空间预测 ?简单和复杂的水质研究 ?环境影响和优化研究 ?规划和可行性研究 ?水质预报 1.2 内置模板和使用手册 ? DHI预定义的ECO Lab模板在以下目录中： Templates ECOLab ?使用手册和说明在以下目录中: ManualsMIKE ZEROECOLab 2 WQ -水质模块 2.1 MIKE 21水质模块的目标 MIKE 21水质模型主要针对湖泊、海洋区域的污水排放引起的水质问题，比如 BOD/DO,富营养化和细菌污染。 2.2目前水质模块可进行以下模拟： ?大肠杆菌，粪大肠杆菌/总大肠杆菌的传输和死亡（用一级降解来表示），降解 速率取决于当地的光

3、强，温度和盐度条件等。 ? BOD-DO关系，即排放的有机物所引起的耗氧。考虑以下几个过程： BOD 一级降解 BOD降解引起的耗氧 底泥需氧量 水体中的呼吸作用 光合作用产氧 水气相互作用下的氧交换（大气复氧） ? BOD-DO模块包括不同营养物（氨氮，硝酸盐和磷）以及三种BOD形式：溶解 性，悬浮性和沉积性 BOD。使用该模块需要设置三种 BOD组分的一级降解速 率。悬浮和沉积的BOD将考虑沉降和再悬浮。该模块中氧平衡过程主要包括： BOD降解需氧量，底泥需氧量，硝化反应需氧量，光合作用产氧，呼吸作用耗 氧以及大气复氧。营养物转化的基本过程包括：BOD降解释放有机氮和磷，产 生的氨氮经硝化

4、反应变成硝酸盐氮，最终通过反硝化作用生成氮气，释放在大 气中。同时，BOD降解所释放的部分氨氮和磷可以被浮游生物，植物和细菌所 吸收。 ?用户可以按实际需求自定义多种污染物质，并定义相应的降解速率进行模拟 典型污染问题 与典型污染问题相关的污染物质有: ?近海水域中与健康相关的微生物 ?耗氧物质 ?营养物质 ?异生化合物，例如有危害性或毒性的化合物 与健康相关的微生物 对于近海水域微生物调查的主要目的在于指出其用水安全性，或是作为对该处鱼 类，贝类等生长环境的调查。一个全面的微生物风险评估包括： ?环境健康评估 包括关于排水管道或污水排放口，雨水排放口的季节性变化，水温，流量，潮汐变 化等信息

5、，以及一个报告和行动系统以确保水质恶化引起的问题能及时通知到健康 权威机构并作出相应处理。 ?指示剂生物体的出现和这些生物体的行为，包括其与物理一化学因素及相关 病原体关联的死亡速率（基于光强、盐度、水温、沉降速率和污染程度 等）。 ?病原体的呈现 耗氧物质 耗氧物质分为溶解性和悬浮性物质，与氧进行生物或生物化学作用，消耗水中的溶 解氧。这些耗氧物质主要是一些不同类型的有机物，具有不同类的降解速率。生化 需氧量（BOD）是间接反映水中能为微生物分解的有机物总量的一个综合指标。有 机物在有氧条件下为微生物分解产生 H2O、CO2和NH3。一般BOD以被检验的水样 在标准条件下5天内的耗氧量为代表

6、，称为 BOD5。 营养物质 许多营养物质都是生物生长的必要元素。适量的营养物对于水中微生物的生长及活 动是必需的，然而，一旦营养物质过量就会引起富营养化，将引起一系列的问题， 如水体污浊，河床底部缺氧，生物沉积量的增加等。富营养化模块可用来模拟这种 情况，因为该模块考虑到藻类对其它物质的直接影响。 在营养物质中氮和磷是最重要的，它们是水生植物生长的控制因子。氮以氨氮和硝 酸盐这两种无机氮的形式存在。许多国家对近海水域中的这些营养物质都设定了浓 度标准。MIKE21水质模型（WQ）就是设计用于评估和这些标准浓度相关的水质问 题。MIKE 21富营养化模块（EU）更为复杂，一般水质问题无需使用。

7、 水质模块状态变量涉及到的主要过程描述： DO: reaera （大气复氧）+ phtsyn （光合作用）-respT （呼吸作用）-BodDecay （BOD 降 解）-SOD （底泥需氧量）-OxygenConsumptionFromNitrification （硝化耗氧） TEMP : Radn （太阳辐射）-Rad_out长波辐射） AMMONIA : AmmoniaReleaseFromBOD （BOD 降解释放氨氮）-Nitrification （硝化）- Pla ntuptake 植物摄取）-bacteriaUptake 细菌摄取） NITRATE : Nitrification

8、 （硝化）-Denitrification （反硝化） BOD : - BodDecay （BOD 降解）-Sedimentation 沉降）+ Resuspension 再悬浮） OP: PPdecay （颗粒型磷降解）-PPformation （颗粒型磷合成）+ OPreleaseFromBOD （BOD降解释放溶解型磷）-OPplantUptake植物摄取溶解型磷） PP: - PPdecay颗粒型磷降解）+ PPformation （颗粒型磷合成）-PPsedimentation颗粒 型磷沉降）+ PPresuspension颗粒型磷再悬浮） FaecalColi: -FaecalCo

9、liDecay （粪大肠杆菌降解） TotColi : -TotalColiDecay （总大肠杆菌降解） 主要常数： ? 降解系数 ? 温度系数 ? 沉降和再悬浮速率 ? 沉降临界速率 ? 产氧速率 ? 呼吸速率 ? 底泥需氧量 ? 耗氧速率(如硝化过程) ? N/P产率和被吸收速率 ? 反应级数 ? 硝化/反硝化速率 主要参数经验值： 1. BOD 级降解速率:0.1 - 0.2 /day 温度系数 1.07 (1.02-1.09) 2. BOD降解过程释放氨氮的典型产出率 原污水：0.065 gNH4/gBOD (0.01-0.1) 生化处理后污水：0.3 gNH4/gBOD (0.1-

10、 0.6) 3. 硝化速率 0.05 /day (0.01 - 0.3) 温度系数1.088 硝化需氧量：4.57gO2/gNH4 4. 反硝化速率 0.1/day (0.05 - 0.3) 温度系数1.16 5. 植物吸收N(光合作用):植物0.066 gNH4/gBOD 6. BOD降解过程中细菌吸收N: 0.109 gNH4/gBOD 7. BOD降解释放无机磷的典型产出率 原污水：0.014 gP/gBOD (0.003-0.03) 生化处理后污水：0.06 gP/gBOD (0.01- 0.09) 8. 植物吸收 P: 0.0091 gP/gBOD 9. 颗粒态磷(PP)降解速率：0

11、.1-0.2/day 10. 光合作用最大值 1.75 -7.0 g O2/m2/d ,相当于 0.5 -2.0 g C/m2/d 11. 底泥需氧量典型值0.2 -1.0 g O2/m2/d (可通过实验获取) 无机底泥取0;河口或有多年排污历史的区域可取1 -2 g O2/m2/d;自然河道 取 0.5 g O2/m2/d 温度系数为1.07 12. 大肠杆菌一级降解速率：0.8/day 温度系数1.07 13. 沉降速率与颗粒物粒径有关：一般对于1 -10卩m取0.07 -0.7 m/d ,典型对于5 卩 m 取 0.2 m/d 14. 再悬浮速率根据率定或相关经验获得. 3 MIKE

12、21水质模型 MIKE 21水质模型所需数据资料： ?基本模型参数： 地形网格（结构网格图或非结构网格图） 时间步长和模拟时间 输出项类型和频率 ?地形和HD条件 ?耦合的AD模型：扩散系数的率定 ?初始值：各参数的浓度值 ?边界条件：各参数的浓度值 ?污染源：坐标位置、水动力条件及各参数的浓度值 ?各生物过程速率值：参考率定值、经验值或监测值等。 在MIKE 21模型中添加 ECO Lab模块 第一步：引入水质模块： 1. 在MIKE 21结构网格模型中引入 ECO Lab模板：打开 MIKE 21 Flow Model （.m21） Basic Parameters Module Sele

13、ction Hydrodynamic and ECO Lab （图 3.1）,出现 ECO Lab Parameters 在 Model Definition里选择适当的内置水质模块或自定义模块，参见图 3.2 图3.1 图3.2 2. 在MIKE 21非结构网格模型中引入 ECO Lab模块：打开 MIKE 21 Flow Model FMModule SelectionECO Lab （参见图 3.3），出现 ECO Lab Module，在 Model Definition里选择适当的内置水质模块或 自定义模块，参见图3.4。 整 New File Product Tpes; Dccun

14、ents: EE議| A FluMudel FF MIKE Zbid MIFE I I NIKE 21 NIKE 3 MIKE 21/3lrteited Models LIT PACK MIKE -LJOOD MIKE SHF FluMudsl n21) 5p=Ctidl Wwes. Bujssintsq BlktlcMH tW) 沖釈tralM, Slop# g,” Flow Model FM Cbical | 图3.3 MIC.E 211 X Dak* dr d Th# 扌 Mn 116： 阳禅 E Ej y二防cdJe 0 x ECO Ls Modd- 為 MuJi-l L efmiij

15、r! d前顾gri士顷 / SflkiiaH rnrrtriEr Csrxlats Fatiw + V DI泊sr討杵 审 IM* Crdlpcni V Dufuif em;late E-Elect cm frc File MIFE 21/3 Wq Shpb- ncLdr Edtorni i MIKE 21Z3 WQpls inckxlng I er phuIldbJiii* 洲KE 21Z3 WQ F.nph nrJudr T fripHMjr巧囲 m悼为刊 GgllbF rctclict- H 口 al 丄如I址 5 hlKF Zera JmllleA Mocfinedl -卜|反 =le

16、 Edit 甘的 E wribw 申|p _X Ld黄匪邑寻冒皆 空 ModDfinilion I Eufcor#*MRn MiitH 1 lEutophcafn仃衬去日mcL Sedimenl and B ehthic=ftcn ICutoplhR6aiifonModVl 2 MIKE H WD L臥暇 I Mll? 11 诃Q Ld 1 arid Eui MIH11 WDL#sd2 MIS-E 11 Urq LrmH 3 Mlffi 11 W3 MIFF 1 MQ 4 nnd Eoi MIKE 11 WiQ Le晶 4 dird Cofi * P Mill H WaLweiS MlkE 1

17、1 WQ 血咬 e H W2 L-wd t and Cai MllbE 11 WC| d U ivd 匚uli * P Mire 21/3WQ 刃卅他冠ert* 材IBJE 21 禺5| 胖ih 111?.讹.昌卅 nhbinphyhei NdlfdGvii im旳Mo T苗袖rgi 图3.4 第二步：在Model Definition里选择计算方法并设置水质模拟时间步长。 水质模拟时间步长通常先以0.5小时进行计算。 第三步：在ECO Lab模块中分别对需模拟的状态变量、边界水质条件、扩散系数、 污染源浓度、ECO Lab模块各参数、作用力和输出项进行设定。 4水质模型应用要点归纳 1. 二

18、维水质模型主要适用于水库、湖泊、河口和海洋等的水质模拟。模拟结果主 要取决于边界条件、污染源负荷以及外部作用力（如温度、太阳辐射、盐度 等）。 2. 若要降低初始值对模型计算结果的影响或继续下一模拟时间段的计算，可采用 hotstart进行再次计算。也就是在模拟时间段内运行2次，第一次是将初始值设 为常数进行模拟（一般采用起始模拟时间实测数据的平均值），第二次是将第 一次运行后的如二维结果图*dfs2或*.dfsu文件作为初始值进行模拟，以减少结 果对于初始值的依赖性和考虑不同研究区域内污染物浓度的梯度变化。 3. WQ模块主要用于：水体中人为因素造成的污染负荷（点源和面源污染）占主体 的情况

19、，污染物在水体中停留的时间相对较短。若对于污染物停留时间较长的 湖泊水库，富营养化程度较大，则应该考虑用EU模型，考虑浮游生物的影响。 4增大水质计算时间步长，可以缩短计算时间，但也会增加模型发散的可能性， 所以通过方案比较选择合适的时间步长进行计算。初始选择时间间隔0.5hr进行 模拟计算。 5. 进行水质模拟前需进行对流扩散模拟，确定对流扩散系数以及检验模型的发散 性。AD模块可简单模拟盐水入侵、温度变化、污染物传输和一级降解。完整的 物理、化学和生物过程在 ECO Lab模块中实现，并与AD耦合使用进行计算。 6. 降解系数的单位在 MIKE 21 AD模块中为/sec,在ECO Lab

20、中为/天。 7. 对ECO Lab中各常数进行率定时，一般采用同一参数值进行率定。若研究区域 内组分浓度梯度变化较大，则使用*dfs2或*.dfsu文件对不同区域的参数值进行 分别设定。 8. 如果实测数据明显有问题，通过对各生态过程中的参数调整仍然不能达到较好 的率定效果，则应该考虑根据实际情况对污染源的流量和负荷进行重新评估， 作出适当修改后输入模型。 5 EU 富营养化模块（附加） 富营养化模块（EU）用来描述水中溶解氧状态，营养物循环，浮游植物和浮游动 物的生长过程以及根系植被和大型藻类生长和分布等。主要应用在内陆水体（如湖 泊、水库等）以及相关的海洋水域，污染问题主要与营养物质负荷有

21、关。 氮和磷通常是基本的营养物质，它们控制着水中浮游植物和大型藻类（如石莼）的 生长，而这些植物会引起富营养化问题。水生生态系统的富营养化状态不仅取决于 营养物质负荷，水温、光照、盐度、水动力条件也相当重要。温度和光辐射是藻类 进行光合作用的必要条件，前者决定细胞内酶促反应的速率，后者是代谢的能源。 各种植物必须在一定的温度下才能进行生理活动，温度太高或太低时都会影响浮游 植物的生长。当光强增加时，浮游植物的生长率也增大，并呈线性关系，直到光合 速率达饱和。海水的温度是赤潮发生的重要因子，20C30C是赤潮发生的最适温 度。盐度的变化是促进生物大量繁殖的原因之一，盐度变化在2637%o范围内均

22、 有发生赤潮的可能。海水盐度在1521.6%。时容易形成温跃层和盐跃层，温、盐跃 层的存在为赤潮生物的聚焦提供条件，易诱发赤潮。在水流缓慢的水体如湖泊、水 库、河口、港湾和内海等处易发生富营养化，这是因为水体分层使得底层得不到氧 气供应，而大量繁殖起来的藻类死亡后又沉入湖底需要耗氧分解，从而使底层缺氧 更为严重。藻类在表层的大量繁殖，使光线不能射入底层，底层藻类不能很好的进 行光合作用，也加剧了底层缺氧。海洋由于径流、涌升流、水团或海流的交汇作 用，使海底层营养盐上升到水上层造成沿海水域高度富营养化。 MIKE软件预设了 3个富营养化模板： EutrciphicationModell .eoo

23、la b Edr叩 hicationModell5edinneritBeriitiicV&getatior.ecolab Eutr ap hi cat io nTdodelz r ecoh b 富营养化模拟的数据收集 需收集的数据主要包括： ?浮游植物生产量，即可用（gO2/m2/d）也可用（gC/ m2/d） ?叶绿素-a浓度（g/m3） ?总氮和总磷（g/m3） ?无机氮和无机磷（g/m3） ?溶解氧（g/m3） ?底栖植被生物量（gC/m2）（如果包含在建模中） 这些数据是必要的，此外相关数据有浮游动物生物量和腐质碳浓度。 这些数据必须覆盖研究区域，采集位置必须靠近模型边界位置。为得到可

24、靠模拟结 果，边界数据尤为重要。 模拟对象： 浮游植物、浮游动物、有机物（腐质）、有机和无机营养物N和P，DO,底栖植 被生物量等。 状态变量有12个（以EU1为例）： PC、PN、PP -浮游植物碳,浮游植物氮,浮游植物磷 CH -叶绿素a ZC -浮游动物 DC、DN、DP -腐质碳，腐质氮，腐质磷 IN、IP -无机氮，无机磷 DO -溶解氧 BC -底栖植物碳 ECO Lab中的富营养化模型也属于生态动力学模型，其描述了营养盐的循环过程, 浮游植物和浮游动物的生长，根生植物以及大型藻类的生长和分布。（主要转化过 程如下图所示）此外还模拟水体中的氧环境。 14 11 ft 浮游动构 1生

25、产，浮游植物；2沉降，浮游植物；3牧食；4死亡，浮游植物；5排泄，浮游动 物；6死亡，浮游动物；7呼吸，浮游动物；8悬浮性腐质的矿化；9腐质沉积；10 腐质矿化；11沉积物中的累积；12生产，底栖植被；13死亡，底栖植被；14与上覆 水体的交换。 以EU1为例，富营养化模块状态变量涉及到的过程描述： 1）Algal Carbo n: PRPC-GRPC-DEPC-SEPC PRPC-Production phytoplankton carbon 浮游植物生产的碳量 GRPC-Grazing of phytoplankton carbon 浮游动物牧食消耗的浮游植物碳量 DEPC-Death o

26、f phytoplankton carbon 浮游植物死亡损失的碳量 SEPC-Settling of phytoplankton carbon 浮游植物沉积损失的碳量 2) Algal Nitroge n: UNPN-GRPN-DEPN-SEPN UNPN-Uptake of phytopla nkton nitroge n 浮游植物吸收的氮量 GRPN-Grazi ng of phytopla nkton nitroge n 浮游动物牧食消耗的浮游植物氮量 DEPN-Death of phytopla nkton nitroge n 浮游植物死亡损失的氮量 SEPN-Settling of

27、 phytoplankton nitrogen 浮游植物沉积损失的氮量 3) Algal Phosphorous: UPPP-GRPP-DEPP-SEPP UPPP-Uptake of phytoplankton phosphorous 浮游植物吸收的磷量 GRPP-Grazing of phytoplankton phosphorous 浮游动物牧食消耗的浮游植物磷量 DEPP-Death of phytoplankton phosphorus 浮游植物死亡损失的磷量 SEPP-Settling of Phytoplankton phosphorous 浮游植物沉积损失的磷量 4) Chlo

28、rophyll-a: PRCH-DECH-SECH PRCH-Production phytopla nkton chlorophyll 浮游植物光合作用生产的叶绿素 DECH-Death of phytopla nkton chlorophyll浮游植物死亡损失的叶绿素 SECH-Settli ng of phytopla nkton chlorophyll 浮游植物沉积过程损失的叶绿素 5) Zoopla nkto n: PRZC-DEZC PRZC-Production of zooplankton carbon 浮游动物生长过程中生产的碳量 DEZC-Death of zooplank

29、ton carbon 浮游动物死亡损耗的碳量 6) Detritus Carbon: DEPC2DC+DEBC2M3+EKZC-REDC-SEDC+DEZC DEPC2DC-Death phytopla nkton to detritus carbon 浮游植物死亡后转化为腐质的碳量 DEBC2M3-Death be nthic vegetation carbo n, per m3 每立方米死亡的底栖植物中的碳量 EKZC-Excretio n by zoopla nkton carbon 浮游动物排泄的碳量 REDC-Respiration detritus carbon腐质呼吸作用中消耗的

30、碳量 SEDC-Settling of detritus carbon腐质沉积进入沉积物的碳量 DEZC-Death of zooplankton carbon 死亡浮游动物的碳量 7) Detritus Nitroge n: DEPN2DN+EKZN-REDN-SEDN+DEZN+DEBN DEPN2DN- Death phytopla nkton to detritus nitroge n 浮游植物死亡后转化为腐质的氮量 EKZN- Excretio n by zoopla nkton nitroge n 浮游动物排泄的氮量 REDN- Respiration detritus nitro

31、gen 腐质呼吸作用中消耗的氮量 SEDN- Settli ng of detritus nitroge n 腐质沉积进入沉积物的氮量 DEZN- Death of zoopla nkton nitroge n 死亡浮游动物所含的氮量 DEBN- Death be nthic vegetation nitroge n 死亡底栖植物所含的氮量 8) Detritus Phosphorous: DEPP2DP+EKZP-REDP-SEDP+DEZP+DEBP DEPP2DP- Death phytoplankton to detritus phosphorous 浮游植物死亡后转化为腐质的磷量 E

32、KZP- Excretio n by zoopla nkton phosphorus 浮游动物排泄的磷量 REDP-Respiration detritus phosphorus腐质呼吸作用中消耗的磷量 SEDP-Settling of detritus phosphorus腐质沉积进入沉积物的磷量 DEZP-Death of zooplankton phosphorus 死亡浮游动物所含的磷量 DEBP-Death benthic vegetation phosphorous 死亡底栖植物所含的磷量 9) Ino rga nic Nitroge n: REDN+REZN+RESN-UNPN+

33、DEPN2IN-UNBN+REBN REDN- Respiration detritus nitroge n 腐质有机氮矿化作用释放的无机氮量 REZN- Respiration of zoopla nkton nitroge n 浮游动物呼吸释放的无机氮量 RESN- Respiration of sediment nitrogen沉积物有机氮矿化作用释放的无机氮量 UNPN- Uptake of phytopla nkton nitroge n 浮游植物吸收的无机氮量 DEPN2IN- Death phytopla nkton to in orga nic nitroge n 死亡浮游植物

34、在矿化作用下释放的无机 氮量 UNBN- Uptake be nthic vegetation nitroge n 底栖植物吸收的无机氮量 REBN- Respiration of be nthic vegetation nitroge n 底栖植物呼吸作用释放的无机氮量 10) Inorganic Phosphorous: REDP+REZP+RESP-UPPP+DEPP2IP-UPBP+REBP REDP- Respiration detritus phosphorus腐质有机磷矿化作用释放的无机磷量 REZP-Respiration of zooplankton phosphorous 浮游动物呼吸作用释放的无机磷量 RESP-Respiration of sediment phosphorous沉积物中的有机磷矿化作用释放的无机磷量 UPPP-Uptake of phytoplankton phosphorou