MIKE水质培训教程

MIKE3水质培训教程DHIChinaECOLab简介ECOLab是DHI在老式的水质模型概念发展起来的全新的水质和生态模拟工具。ECOLab软件开发的理念和措施非常先进，顾客不仅可以修改模型参数，更重要的是可以修改模型关键程序、甚至编写新程序，然后ECOLab将其与MIKE11/21/3的HD、AD集成计算。DHI已经将大部分老式的水质模块转换成ECOLab通用模板，供顾客调用或修改使用，包括：水质模块富营养化模块重金属模块应用领域河流、湿地、湖泊、水库、河口、海岸和海洋各生态系统反应的空间预测简朴和复杂的水质研究环境影响和优化研究规划和可行性研究水质预报内置模板和使用手册DHI预定义的ECOLab模板在如下目录中：C:\ProgramFiles\DHI\\MIKEZero\Templates\ECOLab使用手册和阐明在如下目录中：C:\ProgramFiles\DHI\\MIKEZero\Manuals\MIKE_ZERO\ECOLab

WQ－水质模块2.1MIKE3WQ水质模块的目的MIKE3水质模型重要针对湖泊、海洋区域的污水排放引起的水责问题，例如BOD/DO,富营养化和细菌污染。2.2目前水质模块可进行如下模拟：大肠杆菌，粪大肠杆菌/总大肠杆菌的传播和死亡（用一级降解来表达），降解速率取决于当地的光强，温度和盐度条件等。BOD-DO关系，即排放的有机物所引起的耗氧。考虑如下几种过程：BOD一级降解BOD降解引起的耗氧底泥需氧量水体中的呼吸作用光合作用产氧水气互相作用下的氧互换（大气复氧）BOD-DO模块包括不一样营养物（氨氮，硝酸盐和磷）以及三种BOD形式：溶解性，悬浮性和沉积性BOD。使用该模块需要设置三种BOD组分的一级降解速率。悬浮和沉积的BOD将考虑沉降和再悬浮。该模块中氧平衡过程重要包括：BOD降解需氧量，底泥需氧量，硝化反应需氧量，光合作用产氧，呼吸作用耗氧以及大气复氧。营养物转化的基本过程包括：BOD降解释放有机氮和磷，产生的氨氮经硝化反应变成硝酸盐氮，最终通过反硝化作用生成氮气，释放在大气中。同步，BOD降解所释放的部分氨氮和磷可以被浮游生物，植物和细菌所吸取。顾客可以按实际需求自定义多种污染物质，并定义对应的降解速率进行模拟。

经典污染问题与经典污染问题有关的污染物质有：近海水域中与健康有关的微生物耗氧物质营养物质异生化合物，例如有危害性或毒性的化合物与健康有关的微生物对于近海水域微生物调查的重要目的在于指出其用水安全性，或是作为对该处鱼类，贝类等生长环境的调查。一种全面的微生物风险评估包括：环境健康评估包括有关排水管道或污水排放口，雨水排放口的季节性变化，水温，流量，潮汐变化等信息，以及一种汇报和行动系统以保证水质恶化引起的问题能及时告知到健康权威机构并作出对应处理。指示剂生物体的出现和这些生物体的行为，包括其与物理－化学原因及有关病原体关联的死亡速率（基于光强、盐度、水温、沉降速率和污染程度等）。病原体的展现耗氧物质耗氧物质分为溶解性和悬浮性物质，与氧进行生物或生物化学作用，消耗水中的溶解氧。这些耗氧物质重要是某些不一样类型的有机物，具有不一样类的降解速率。生化需氧量（BOD）是间接反应水中能为微生物分解的有机物总量的一种综合指标。有机物在有氧条件下为微生物分解产生H2O、CO2和NH3。一般BOD以被检查的水样在原则条件下5天内的耗氧量为代表，称为BOD5。营养物质许多营养物质都是生物生长的必要元素。适量的营养物对于水中微生物的生长及活动是必需的，然而，一旦营养物质过量就会引起富营养化，将引起一系列的问题，如水体污浊，河床底部缺氧，生物沉积量的增长等。富营养化模块可用来模拟这种状况，由于该模块考虑到藻类对其他物质的直接影响。在营养物质中氮和磷是最重要的，它们是水生植物生长的控制因子。氮以氨氮和硝酸盐这两种无机氮的形式存在。许多国家对近海水域中的这些营养物质都设定了浓度原则。MIKE3水质模型(WQ)就是设计用于评估和这些原则浓度有关的水责问题。MIKE3富营养化模块（EU）更为复杂，一般水责问题无需使用。水质模块状态变量波及到的重要过程描述：DO：reaera(大气复氧)+phtsyn(光合作用)–respT(呼吸作用)–BodDecay(BOD降解)-SOD(底泥需氧量)–OxygenConsumptionFromNitrification(硝化耗氧)TEMP：Rad_in(太阳辐射)-Rad_out(长波辐射)AMMONIA：AmmoniaReleaseFromBOD(BOD降解释放氨氮)–Nitrification(硝化)–Plantuptake(植物摄取)–bacteriaUptake(细菌摄取)NITRATE：Nitrification(硝化)–Denitrification(反硝化)BOD：-BodDecay(BOD降解)–Sedimentation(沉降)+Resuspension(再悬浮)OP：PPdecay(颗粒型磷降解)–PPformation(颗粒型磷合成)+OPreleaseFromBOD(BOD降解释放溶解型磷)–OPplantUptake(植物摄取溶解型磷)PP：-PPdecay(颗粒型磷降解)+PPformation(颗粒型磷合成)-PPsedimentation(颗粒型磷沉降)+PPresuspension(颗粒型磷再悬浮)FaecalColi：-FaecalColiDecay(粪大肠杆菌降解)TotColi：-TotalColiDecay(总大肠杆菌降解)重要常数：降解系数温度系数沉降和再悬浮速率沉来临界速率产氧速率呼吸速率底泥需氧量耗氧速率（如硝化过程）N/P产率和被吸取速率反应级数硝化/反硝化速率重要参数经验值：BOD一级降解速率:0.1-0.2/day温度系数1.07(1.02-1.09)BOD降解过程释放氨氮的经典产出率原污水:0.065gNH4/gBOD(0.01-0.1)生化处理后污水:0.3gNH4/gBOD(0.1-0.6)硝化速率0.05/day(0.01-0.3)温度系数1.088硝化需氧量:4.57gO2/gNH4反硝化速率0.1/day(0.05-0.3)温度系数1.16植物吸取N(光合作用):植物0.066gNH4/gBODBOD降解过程中细菌吸取N:0.109gNH4/gBODBOD降解释放无机磷的经典产出率原污水:0.014gP/gBOD(0.003-0.03)生化处理后污水:0.06gP/gBOD(0.01-0.09)植物吸取P:0.0091gP/gBOD颗粒态磷(PP)降解速率:0.1-0.2/dayMIKE3水质模型MIKE3水质模型所需数据资料：基本模型参数：地形网格（构造网格图或非构造网格图）时间步长和模拟时间输出项类型和频率地形和HD条件耦合的AD模型：扩散系数的率定初始值：各参数的浓度值边界条件：各参数的浓度值污染源：坐标位置、水动力条件及各参数的浓度值各生物过程速率值：参照率定值、经验值或监测值等。在MIKE3模型中添加ECOLab模块第一步：引入水质模块：在MIKE3构造网格模型中引入ECOLab模块：打开MIKE3FlowModel(.m3)BasicParametersModuleSelectionEnvironmentalmodulesECOLab(图3.1)，出现ECOLabParameters，在ModelDefinition里选择合适的内置水质模块或自定义模块，参见图3.2。图3.1图3.2在MIKE3非构造网格模型中引入ECOLab模块：打开MIKE3FlowModelFMModuleSelectionECOLab（参见图3.3），出现ECOLabModule，在ModelDefinition里选择合适的内置水质模块或自定义模块，参见图3.4。图3.3图3.4第二步：在ModelDefinition里选择计算措施并设置水质模拟时间步长。水质模拟时间步长一般先以0.5小时进行计算。第三步：在ECOLab模块中分别对需模拟的状态变量、边界水质条件、扩散系数、污染源浓度、ECOLab模块各参数、作用力和输出项进行设定。水质模型应用要点归纳三维水质模型重要合用于水库、湖泊、河口和海洋等的水质模拟。模拟成果重要取决于边界条件、污染源负荷以及外部作用力（如温度、太阳辐射、盐度等）。若要减少初始值对模型计算成果的影响或继续下一模拟时间段的计算，可采用hotstart进行再次计算。也就是在模拟时间段内运行2次，第一次是将初始值设为常数进行模拟（一般采用起始模拟时间实测数据的平均值），第二次是将第一次运行后的如三维成果图*.dfs3或*.dfsu文献作为初始值进行模拟，以减少成果对于初始值的依赖性和考虑不一样研究区域内污染物浓度的梯度变化。WQ模块重要用于：水体中人为原因导致的污染负荷（点源和面源污染）占主体的状况，污染物在水体中停留的时间相对较短。若对于污染物停留时间较长的湖泊水库，富营养化程度较大，则应当考虑用EU模型，考虑浮游生物的影响。增大水质计算时间步长，可以缩短计算时间，但也会增长模型发散的也许性，因此通过方案比较选择合适的时间步长进行计算。初始选择时间间隔0.5hr进行模拟计算。进行水质模拟前需进行对流扩散模拟，确定对流扩散系数以及检查模型的发散性。AD模块可简朴模拟盐水入侵、温度变化、污染物传播和一级降解。完整的物理、化学和生物过程在ECOLab模块中实现，并与AD耦合使用进行计算。降解系数的单位在MIKE3AD模块中为/sec，在ECOLab中为/天。对ECOLab中各常数进行率定期，一般采用同一参数值进行率定。若研究区域内组分浓度梯度变化较大，则使用*.dfs3或*.dfsu文献对不一样区域的参数值进行分别设定。假如实测数据明显有问题，通过对各生态过程中的参数调整仍然不能到达很好的率定效果，则应当考虑根据实际状况对污染源的流量和负荷进行重新评估，作出合适修改后输入模型。EU－富营养化模块（附加）富营养化模块（EU）用来描述水中溶解氧状态，营养物循环，浮游植物和浮游动物的生长过程以及根系植被和大型藻类生长和分布等。重要应用在内陆水体(如湖泊、水库等)以及有关的海洋水域，污染问题重要与营养物质负荷有关。氮和磷一般是基本的营养物质,它们控制着水中浮游植物和大型藻类（如石莼）的生长，而这些植物会引起富营养化问题。水生生态系统的富营养化状态不仅取决于营养物质负荷，水温、光照、盐度、水动力条件也相称重要。温度和光辐射是藻类进行光合作用的必要条件，前者决定细胞内酶促反应的速率，后者是代谢的能源。多种植物必须在一定的温度下才能进行生理活动，温度太高或太低时都会影响浮游植物的生长。当光强增长时，浮游植物的生长率也增大，并呈线性关系，直到光合速率达饱和。海水的温度是赤潮发生的重要因子，20℃～30℃是赤潮发生的最适温度。盐度的变化是增进生物大量繁殖的原因之一，盐度变化在26～37‰范围内均有发生赤潮的也许。海水盐度在15～21.6‰时轻易形成温跃层和盐跃层，温、盐跃层的存在为赤潮生物的聚焦提供条件，易诱发赤潮。在水流缓慢的水体如湖泊、水库、河口、港湾和内海等处易发生富营养化，这是由于水体分层使得底层得不到氧气供应，而大量繁殖起来的藻类死亡后又沉入湖底需要耗氧分解，从而使底层缺氧更为严重。藻类在表层的大量繁殖，使光线不能射入底层，底层藻类不能很好的进行光合作用，也加剧了底层缺氧。海洋由于径流、涌升流、水团或海流的交汇作用，使海底层营养盐上升到水上层导致沿海水域高度富营养化。MIKE软件预设了3个富营养化模板：富营养化模拟的数据搜集需搜集的数据重要包括：浮游植物生产量，即可用（gO2/m2/d）也可用（gC/m2/d）叶绿素-a浓度（g/m3） 总氮和总磷（g/m3）无机氮和无机磷（g/m3）溶解氧（g/m3）底栖植被生物量（gC/m2）（假如包括在建模中）这些数据是必要的，此外有关数据有浮游动物生物量和腐质碳浓度。这些数据必须覆盖研究区域，采集位置必须靠近模型边界位置。为得到可靠模拟成果，边界数据尤为重要。模拟对象：浮游植物、浮游动物、有机物（腐质）、有机和无机营养物N和P，DO，底栖植被生物量等。状态变量有12个（以EU1为例）：PC、PN、PP-浮游植物碳,浮游植物氮,浮游植物磷CH-叶绿素aZC-浮游动物DC、DN、DP-腐质碳，腐质氮，腐质磷IN、IP-无机氮，无机磷DO-溶解氧BC-底栖植物碳ECOLab中的富营养化模型也属于生态动力学模型，其描述了营养盐的循环过程，浮游植物和浮游动物的生长，根生植物以及大型藻类的生长和分布。(重要转化过程如下图所示）此外还模拟水体中的氧环境。1生产，浮游植物；2沉降，浮游植物；3牧食；4死亡，浮游植物；5排泄，浮游动物；6死亡，浮游动物；7呼吸，浮游动物；8悬浮性腐质的矿化；9腐质沉积；10腐质矿化；11沉积物中的累积；12生产，底栖植被；13死亡，底栖植被；14与上覆水体的互换。以EU1为例，富营养化模块状态变量波及到的过程描述：AlgalCarbon:PRPC-GRPC-DEPC-SEPCPRPC-Productionphytoplanktoncarbon浮游植物生产的碳量GRPC-Grazingofphytoplanktoncarbon浮游动物牧食消耗的浮游植物碳量DEPC-Deathofphytoplanktoncarbon浮游植物死亡损失的碳量SEPC-Settlingofphytoplanktoncarbon浮游植物沉积损失的碳量AlgalNitrogen:UNPN-GRPN-DEPN-SEPNUNPN-Uptakeofphytoplanktonnitrogen浮游植物吸取的氮量GRPN-Grazingofphytoplanktonnitrogen浮游动物牧食消耗的浮游植物氮量DEPN-Deathofphytoplanktonnitrogen浮游植物死亡损失的氮量SEPN-Settlingofphytoplanktonnitrogen浮游植物沉积损失的氮量AlgalPhosphorous:UPPP-GRPP-DEPP-SEPPUPPP-Uptakeofphytoplanktonphosphorous浮游植物吸取的磷量GRPP-Grazingofphytoplanktonphosphorous浮游动物牧食消耗的浮游植物磷量DEPP-Deathofphytoplanktonphosphorus浮游植物死亡损失的磷量SEPP-SettlingofPhytoplanktonphosphorous浮游植物沉积损失的磷量Chlorophyll-a:PRCH-DECH-SECHPRCH-Productionphytoplanktonchlorophyll浮游植物光合作用生产的叶绿素DECH-Deathofphytoplanktonchlorophyll浮游植物死亡损失的叶绿素SECH-Settlingofphytoplanktonchlorophyll浮游植物沉积过程损失的叶绿素Zooplankton:PRZC-DEZCPRZC-Productionofzooplanktoncarbon浮游动物生长过程中生产的碳量DEZC-Deathofzooplanktoncarbon浮游动物死亡损耗的碳量DetritusCarbon:DEPC2DC+DEBC2M3+EKZC-REDC-SEDC+DEZCDEPC2DC-Deathphytoplanktontodetrituscarbon浮游植物死亡后转化为腐质的碳量DEBC2M3-Deathbenthicvegetationcarbon,perm3每立方米死亡的底栖植物中的碳量EKZC-Excretionbyzooplanktoncarbon浮游动物排泄的碳量REDC-Respirationdetrituscarbon腐质呼吸作用中消耗的碳量SEDC-Settlingofdetrituscarbon腐质沉积进入沉积物的碳量DEZC-Deathofzooplanktoncarbon死亡浮游动物的碳量DetritusNitrogen:DEPN2DN+EKZN-REDN-SEDN+DEZN+DEBNDEPN2DN-Deathphytoplanktontodetritusnitrogen浮游植物死亡后转化为腐质的氮量EKZN-Excretionbyzooplanktonnitrogen浮游动物排泄的氮量REDN-Respirationdetritusnitrogen腐质呼吸作用中消耗的氮量SEDN-Settlingofdetritusnitrogen腐质沉积进入沉积物的氮量DEZN-Deathofzooplanktonnitrogen死亡浮游动物所含的氮量DEBN-Deathbenthicvegetationnitrogen死亡底栖植物所含的氮量DetritusPhosphorous:DEPP2DP+EKZP-REDP-SEDP+DEZP+DEBPDEPP2DP-Deathphytoplanktontodetritusphosphorous浮游植物死亡后转化为腐质的磷量EKZP-Excretionbyzooplanktonphosphorus浮游动物排泄的磷量REDP-Respirationdetritusphosphorus腐质呼吸作用中消耗的磷量SEDP-Settlingofdetritusphosphorus腐质沉积进入沉积物的磷量DEZP-Deathofzooplanktonphosphorus死亡浮游动物所含的磷量DEBP-Deathbenthicvegetationphosphorous死亡底栖植物所含的磷量InorganicNitrogen:REDN+REZN+RESN-UNPN+DEPN2IN-UNBN+REBNREDN-Respirationdetritusnitrogen腐质有机氮矿化作用释放的无机氮量REZN-Respirationofzooplanktonnitrogen浮游动物呼吸释放的无机氮量RESN-Respirationofsedimentnitrogen沉积物有机氮矿化作用释放的无机氮量UNPN-Uptakeofphytoplanktonnitrogen浮游植物吸取的无机氮量DEPN2IN-Deathphytoplanktontoinorganicnitrogen死亡浮游植物在矿化作用下释放的无机氮量UNBN-Uptakebenthicvegetationnitrogen底栖植物吸取的无机氮量REBN-Respirationofbenthicvegetationnitrogen底栖植物呼吸作用释放的无机氮量InorganicPhosphorous:REDP+REZP+RESP-UPPP+DEPP2IP-UPBP+REBPREDP-Respirationdetritusphosphorus腐质有机磷矿化作用释放的无机磷量REZP-Respirationofzooplanktonphosphorous浮游动物呼吸作用释放的无机磷量RESP-Respirationofsedimentphosphorous沉积物中的有机磷矿化作用释放的无机磷量UPPP-Uptakeofphytoplanktonphosphorous浮游植物吸取的无机磷量DEPP2IP-Deathphytoplanktontoinorganicphosphorous死亡浮游植物在矿化作用下释放的无机磷量UPBP-Uptakebenthicvegetationphosphorous底