水质keepass教程

1MIKE21水质培训教程目录MIKE21水质培训教程.
11ECOLab简介11.
1应用领域11.
2内置模块和使用手册12WQ-水质模块22.
1MIKE21WQ水质模块的目标.
22.
2目前水质模块可进行以下模拟:23EU-富营养化模块.
64MIKE21水质模型114.
1基本水质模块设置114.
2仅用AD进行水质模拟.
134.
3ECOLab与AD耦合进行水质模拟.
135MIKE21水质模型练习155.
1No.
1-对已有内置模块的简单修改.
155.
2No.
2-创建简单的大肠杆菌模块.
175.
3No.
3-BOD/DO和大肠杆菌.
185.
4No.
4-使用ECOLab模拟盐度和CODcr.
225.
5No.
5-暴雨污水溢流和氮过程模拟.
235.
6No.
6-水质率定.
255.
7No.
7-渤海湾水质练习.
276水质模型应用要点归纳3811ECOLab简介ECOLab软件开发的理念和方法非常先进,用户不仅可以修改模型参数,更重要的是可以修改模块核心程序,甚至编写新程序,然后ECOLab将其与MIKE11/21/3的HD、AD集成计算,实现将各水生生态系统转化为可靠的数值模拟用于精确的预报.
1.
1应用领域河流、湿地、湖泊、水库、河口、海岸和海洋各生态系统反应的空间预测简单和复杂的水质研究环境影响和优化研究规划和可行性研究水质预报1.
2内置模块和使用手册DHI预定义的ECOLab模块在以下目录中:C:\ProgramFiles\DHI\MIKEZero\Templates\ECOLab使用手册和说明在以下目录中:C:\ProgramFiles\DHI\MIKEZero\Manuals\MIKE_ZERO\ECOLab\ECOLAB_UserGuide.
pdf(中文版本参见ECOLAB使用手册.
doc)C:\ProgramFiles\DHI\MIKEZero\Manuals\MIKE_21\FlowModel\EL\M21EL.
pdf22WQ-水质模块2.
1MIKE21WQ水质模块的目标MIKE21水质模型主要针对海洋区域的污水排放引起的水质问题,比如BOD/DO,富营养化和细菌污染.
该模型还可以简单模拟富营养化对光合作用的短期效应(以叶绿素a来表示).
如需模拟营养物质和浮游生态系统的基本转化,则需使用MIKE21EU富营养化模块.
2.
2目前水质模块可进行以下模拟:大肠杆菌,粪大肠杆菌/总大肠杆菌的传输和死亡(用一级降解来表示),降解速率取决于当地的光强,温度和盐度条件等.
BOD-DO关系,即排放的有机物所引起的耗氧.
考虑以下几个过程:BOD一级降解BOD降解引起的耗氧底泥需氧量水体中的呼吸作用光合作用产氧水气相互作用下的氧交换(大气复氧)BOD-DO模块包括不同营养物(氨氮,硝酸盐和磷)以及三种BOD形式:溶解性,悬浮性和沉积BOD.
使用该模块需要设置三种BOD组分的一级降解速率.
悬浮和沉积的BOD将考虑沉降和再悬浮.
该模块中氧平衡过程主要包括:BOD降解需氧量,底泥需氧量,硝化反应需氧量,光合作用产氧,水体呼吸作用耗氧以及大气复氧.
营养物转化的基本过程包括:BOD降解释放有机氮和磷,产生的氨氮经硝化反应变成硝酸盐氮,最终通过反硝化作用生成氮气,释放在大气中.
同时,BOD降解所释放的部分氨氮和磷可以被浮游生物,植物和细菌所吸收.
用户可以实际需求自定义多种污染物质,并定义相应的降解速率进行模拟.
3典型污染问题与典型污染问题相关的污染物质有:近海水域中与健康相关的微生物耗氧物质营养物质异生化合物,例如有危害性或毒性的化合物与健康相关的微生物对于近海水域微生物调查的主要目的在于指出其用水安全性,或是作为对该处鱼类,贝类等生长环境的调查.
一个全面的微生物风险评估包括:环境健康评估包括关于排水管道或污水排放口,雨水排放口的季节性变化,水温,流量,潮汐变化等信息,以及一个报告和行动系统以确保水质恶化引起的问题能及时通知到健康权威机构并作出相应处理.
指示剂生物体的出现和这些生物体的行为,包括其与物理-化学因素及相关病原体关联的死亡速率(基于光强、盐度、水温、沉降速率和污染程度等).
病原体的呈现耗氧物质耗氧物质分为溶解性和悬浮性物质,与氧进行生物或生物化学作用,消耗水中的溶解氧.
这些耗氧物质主要是一些不同类型的有机物,具有不同类的降解速率.
生化需氧量(BOD)是间接反应水中能为微生物分解的有机物总量的一个综合指标.
有机物在有氧条件下为微生物分解产生H2O、CO2和NH3.
一般BOD以被检验的水样在标准条件下5天内的耗氧量为代表,称为BOD5.
营养物质许多营养物质都是生物生长的必要元素.
适量的营养物对于水中微生物的生长及活动是必需的,然而,一旦营养物质过量就会引起富营养化,将引起一系列的问题,4如水体污浊,河床底部缺氧,生物沉积量的增加等.
MIKE21富营养化模块可用来模拟这种情况,因为该模块考虑到藻类对其它物质的直接影响.
在营养物质中氮和磷是最重要的,它们是水生植物生长的控制因子.
氮以氨氮和硝酸盐这两种无机氮的形式存在.
许多国家对近海水域中的这些营养物质都设定了浓度标准.
MIKE21水质模型(WQ)就是设计用于评估和这些标准浓度相关的水质问题.
MIKE21富营养化模块(EU)更为复杂,一般水质问题无需使用.
在目前国内的实际项目中,进入河流的污染负荷主要来源于生活、工业和农业,因此多用四级WQ模块进行基本水质指标的模拟和率定.
水质模块状态变量涉及到的过程描述(以WQLevel4为例):DO:reaera+phtsyn-respT-BodDecay-SOD–OxygenConsumptionFromNitrificationTEMP:Rad_in-Rad_outAMMONIA:AmmoniaReleaseFromBOD-Nitrification-Plantuptake–bacteriaUptakeNITRATE:Nitrification–DenitrificationBOD:-BodDecay-Sedimentation+ResuspensionOP:PPdecay-PPformation+OPreleaseFromBOD–OpplantUptakePP:-PPdecay+PPformation-PPsedimentation+PpresuspensionFaecalColi:-FaecalColiDecayTotColi:-TotalColiDecay主要常数:降解系数温度系数沉降和再悬浮速率沉降临界速率产氧速率呼吸速率底泥需氧量5耗氧速率(如硝化过程)N/P产率和被吸收速率反应级数硝化/反硝化速率主要参数经验值:1.
BOD一级降解速率:0.
1-0.
5/day温度系数1.
02-1.
092.
BOD降解过程释放氨氮的典型产出率原污水:0.
065gNH4/gBOD(0.
01-0.
1)生化处理后污水:0.
3gNH4/gBOD(0.
1-0.
6)3.
硝化速率0.
05/day(0.
01-0.
3)温度系数1.
088硝化需氧量:4.
57gO2/gNH44.
反硝化速率0.
1/day(0.
05-0.
3)温度系数1.
165.
植物吸收N(光合作用):植物0.
066gNH4/gBOD6.
BOD降解过程中细菌吸收N:0.
109gNH4/gBOD.
7.
BOD降解释放无机磷的典型产出率原污水:0.
014gP/gBOD(0.
003-0.
03)生化处理后污水:0.
06gP/gBOD(0.
01-0.
09)8.
植物吸收P:0.
0091gP/gBOD9.
颗粒态磷(PP)降解速率:0.
1-0.
2/day63EU-富营养化模块本章节主要介绍两部分内容.
一是描述目前MIKE21富营养化模块,二是大致描述富营养化和污染问题.
MIKE21富营养化模块的目标富营养化模块(EU)用来描述水中溶解氧状态,营养物循环,浮游植物和浮游动物的生长过程以及根系植被和大型藻类生长和分布等.
主要应用在相关的海洋水域,污染问题主要与营养物质的负荷有关.
氮和磷通常是基本的营养物质,它们控制着水中浮游植物和大型藻类(如石莼)的生长.
这些植物会引起富营养化问题.
目前富营养化模块可进行以下模拟:浮游植物的生长,主要与氮和磷的营养物负荷、被研究系统中与水动力过程相结合的内部生物和化学过程有关.
由于水中浮游植物的数量以及死亡有机物(腐质)所造成的水体透明度的下降.
腐质主要来自生物过程中产生的有机物质或污染源排放产生的腐质.
底栖植被(如大型藻类)的生长,主要取决于与浮游植物竞争所获得的营养物质,这也将影响光的穿透性以及能够到达底部供大型藻类生长所需的光照.
由于冷却水排放造成的海洋温度变化给生物过程带来的影响.
与各生物活动(产氧、藻类呼吸作用,死亡有机物即腐质的降解和沉积物呼吸作用)相关的氧平衡.
典型污染问题与正常、自然的水体相比,若水生态系统中营养物过剩则出现富营养化.
通常氮和磷营养物质的过剩是由人类活动造成的.
自然界中营养物质的来源主要包括:未经处理或经过不同程度处理的生活污水7粪便播撒以及化肥使用等农业活动工业废水大气沉降,来源于矿物燃料的燃烧以及粪便/化肥的蒸发作用(如氨氮)海洋系统的富营养化状态不仅仅取决于营养物质的负荷,该区域的水动力条件也相当重要.
将等量的营养物质分别以强流形式排入远海,或以低水量交换方式排入浅海礁湖,两者对当地造成的影响不尽相同.
以下简要描述一些富营养化模块的典型运用.
废水管理全球许多接受城市污水和工业废水排放的沿海地区正在饱受着由于富营养化引起的藻华、大型藻类生长失控以及海藻腐烂时大量耗氧等灾害之苦.
这最终导致了生态平衡被打破,海洋生态系统被恶化.
对这些区域的生活和工业废水进行必需的治理通常会考虑以下问题:若无污染负荷进入的自然海洋生态系统是怎样的,其藻类生长以及其他生物过程的叁数是如何设置的藻类生长的控制因子是磷、氮还是两者相互作用的结果仅减少某一种营养物(如磷)还是同时减少氮磷两种营养物对生态系统更为有效为获得合适的生态条件,营养物质排放浓度降低到何种程度是必需的与污染源排放的污水相比,诸如大气沉降等产生的营养物质会起到什么作用在污水处理的规划阶段,会有以上和其他许多问题产生,但最终都会选择一个合适的方案来改善海洋环境效应.
使用富营养化模块,并与可行性技术方案、相关成本和环境法律法规相结合,进行结果模拟.
海岸生态系统管理8许多海洋生态系统,比如河口,峡湾,咸水湖和海湾,已经由于农业污水排放和地表径流而导致严重的富营养化.
尤其在一些浅水区域中,浮游植物和底栖植被之间的相互作用和竞争已成为生物状态的主导因子.
自然状态下这些系统的初级生产通常受根系植被的影响,如鳗草(绣眼儿属),导致浮游植物生长所需的营养物质很少,这就意味着水中微藻类的浓度较低,水体透明度较高.
当营养物输入增加时,水体透明度下降,导致浮游植物增长.
随着光照条件的下降,在较深海域内根系植被逐渐消失.
而在较浅的海域中可以看见大型藻类的生长.
在富营养化阶段,浮游植物在较深海域的生态系统中占主导地位.
此时透光度低,水体呈"绿色".
根系植被仅在浅水区域里稀疏地生长,导致初级生产力的是主要因子是大型藻类.
尤其在沉积物中,所产生的有机物在降解过程中会增加耗氧量,导致在沉积物中产生缺氧期.
缺氧状态以及由于在沉积物中厌氧条件产生的有毒气体(如硫化氢),会对底部的动物群和鱼类繁衍有害,而这些海域通常是鱼类重要的栖息地.
海洋生态系统的修复项目通常会考虑以下问题:哪种营养物排放该对富营养化贡献最大与季节性变化明显的农业负荷相比,连续排放的污水负荷对水体会产生何种效果推荐何种补救措施来恢复生态系统,比如污水处理、减少农业化肥使用量、沉积物中营养物的疏浚等.
提高水量交换会产生何种效果例如:在沿海泻湖区是否可以改善生态条件营养物在进入沿海生态系统前有多少部分将被滞留,有多少部分可到达公海,系统的同化能力是多少在实际例子中可能出现其他疑问,可能必须考虑许多其他因素,在此不一一进行介绍.
9冷却水研究来自发电厂的冷却水排放对沿海地区的海洋生态系统有各种不同的影响.
主要有:由于出口区域温度升高对初级生产的影响,经过发电厂冷却系统会导致浮游动物的大量死亡.
较小的影响可能表现为水流变化,尤其在许多周边区域.
值得注意的是,除冷却水效应之外,发电厂使用的有毒物质可能会影响生物体.
例如:管道防垢剂的使用.
岸上建筑的环境效应港口、桥梁、堤坝和垃圾填埋场等建筑物在建设期内或在竣工后都会给海洋生态带来一定影响.
在工程期间,尤其用开挖操作时,开挖沙石和倾倒材料可能对植物和动物群产生有害影响.
这种操作不仅会在物理上破坏和减少天然栖息地,而且悬挂材料会降低光的穿透力,限制底部植物的生长.
富营养化模拟的数据收集模型率定所收集的必要数据包括:浮游植物生产量,既可用(gO2/m2/d)也可用(gC/m2/d)叶绿素-a浓度(g/m3)无机氮和磷(g/m3)总氮和磷(g/m3)氧(g/m3)底栖植被生物量(gC/m2)(如果包含在建模中)这些数据是必要的,此外其他相关数据有浮游动物生物量和腐质碳浓度.
这些数据必须覆盖研究区域,采集位置必须靠近模型边界位置.
为得到可靠模拟结果,边界数据尤为重要.
模拟对象:10浮游植物、浮游动物、有机物(腐质)、有机和无机的营养物N和P,DO,底栖植被生物量等.
状态变量有12个(以EU1为例):PC、PN、PPCHZCDC、DN、DPIN、IPDOBC富营养化模块状态变量涉及到的过程描述:1)AlgalCarbon:PRPC-GRPC-DEPC-SEPC2)AlgalNitrogen:PRPN-GRPN-DEPN-SEPN3)AlgalPhosphorous:PRPP-GRPP-DEPP-SEPP4)Chlorophyll-a:PRCH-DECH-SECH5)Zooplankton:PRZC-EKZC-DEZC6)DetritusCarbon:DEPC+SLMC+SLAC+EKZC+DEZC-REDC-SEDC7)DetritusNitrogen:DEPN+SLMN+SLAN+EKZN+DEZN-REDN-SEDN8)DetritusPhosphorous:DEPP+SLMP+SLAP+EKZP+DEZP-REDP-SEDP9)InorganicNitrogen:REDN+REZN+RESN-UNMN+REBN10)InorganicPhosphorous:REDP+REZP+RESP-UNMP+REBP11)DissolvedOxygen:ODPC+ODMC+ODAC-ODDC-ODZC-ODSC+REAR12)Macroalgae:PRMC-SLMC-REMC114MIKE21水质模型基本水质模块设置MIKE21水质模型所需数据资料:基本模型参数:模型网格大小和范围;时间步长和模拟时间;输出项类型和频率地形和HD条件耦合的AD模型:扩散系数的率定初始值:各参数的浓度值边界条件:各参数的浓度值污染源:坐标位置、水动力条件及各参数的浓度值各生物过程速率值:参考率定值、经验值或监测值等.
在MIKE21模型中添加ECOLab模块第一步:引入水质模块:MIKE21FlowModelBasicParametersModuleSelectionHydrodynamicandECOLab,出现ECOLabParameters,参见图4.
1.
第二步:在ModelDefinition里选择适当的内置水质模块或自定义模块.
水质模拟时间步长通常先以0.
5小时进行计算.
第三步:在ECOLab模块中分别对需模拟的状态变量、边界水质条件、扩散系数、污染源浓度、ECOLab模块各参数、作用力和输出项进行设定.
12图4.
113仅用AD进行水质模拟AD模拟物质在水体中的对流和扩散过程,可以设定一个恒定的衰减常数模拟非保守物质,满足一级反应方程式KCdtdC.
可以把MIKE21AD作为简单的水质模型,其衰减系数仅是简单的一级衰减反应系数.
真正的水质模型和生态模型是ECOLab.
使用AD进行水质模拟包括以下几个方面:1.
组分定义:Decaying,conservative,Heatdissipating和HeatExchange2.
初始浓度值设定:常数或dfs2文件初始条件对模型的结果影响较小,一般能在较短的模拟时间段内趋于稳定.
一般初始值取模拟起始时间水质指标的监测平均值.
若使用hotstart文件进行AD计算,则初始值不再参与计算.
3.
边界浓度值的设定:常数、dfs0文件或dfs1文件4.
衰减系数的设定:常数或dfs0文件(考虑随时间步长的改变)5.
扩散系数的设定:扩散系数为率定参数,是进行水质模拟相当重要的一个基础参数,一般可以利用保守物质,如盐度或示踪剂进行率定.
在MIKE21中,网格间距大的地形其扩散系数要比网格间距小的地形的扩散系数大,而且X和Y方向的取值可以不同.
ECOLab与AD耦合进行水质模拟1.
进行MIKE21的深度水质模拟,必须使用ECOLab模块,此时AD模块会耦合在ECOLab模块中.
对水质模拟起关键作用的是水质组分和扩散系数的设定.
2.
在ECOLab模块中详细描述了各状态变量的物理、生化过程.
水质模块包括对状态变量初始浓度值、边界条件、扩散系数、源汇、ECOLab模板参数、作用力以及输出结果类型等的设定.
3.
ECOLab根据模拟水质过程从简单到复杂可分为许多级.
每级模块都分别定义了需模拟的组分及其排列顺序.
4.
DHI自定义的水质模块主要有以下一些模块:14若对标准DHI模块以外的水质组分进行模拟,必须在ECOLab里预先定义该状态变量以及相关的物理、生化过程,再将用户自定义的模块导入MIKE21模型.
155MIKE21水质模型练习No.
1-对已有内置模块的简单修改练习目的:熟悉ECOLab模块的6大组成部分(STATE_VARIABLES,CONSTANTS,FORCINGS,AUXILIARY_VARIABLES,PROCESSES和DERIVED_OUTPUTS)以及如何进行相互转化和修改.
打开C:\ProgramFiles\DHI\MIKEZero\Templates\ECOLabWQlevel4ColiPhos.
ecolab,另存在Mike21Exercises\BasicTemplate\WQlevel4ColiPhos_update.
ecolab练习内容:1、添加衍生变量DERIVED_OUTPUTSTIN=AMMONIA+NITRATETIP=OP+PPCODcr=K*BOD,常数K=32、将温度TEMP从STATE_VARIABLES变成FORCING进行模拟计算.
3、创建盐度和CODcr的ECOLab模块MIKEZeroECOLab,如图5.
1OK,双击打开ECOLab参数文件编辑器,进入空白的水质模块编辑对话框,如图5.
2所示.
假设盐度为保守物质,不考虑衰减过程.
CODcr考虑简单一级衰减过程,满足一级反应方程式CODcrKcoddtdCODcr*,分别定义状态变量、衰减系数和过程,创建一个简单的盐度和CODcr降解过程的ECOLab模块,文件命名为S_CODcr.
ecolab,保存在Mike21Exercises\BasicTemplate目录中.
16图5.
1图5.
217No.
2-创建简单的大肠杆菌模块练习目的:创建一个简单的SimpleColi.
ecolab模板,描述水生态环境中的大肠杆菌变化.
进一步掌握ECOLab的模块结构.
假定大肠杆菌在水中的降解是采用一级衰减公式,衰减与温度有关.
为了简化,在本例中我们不考虑盐度与光辐射对降解的影响,虽然这两项指标对大肠杆菌的降解有着重要影响.
模块保存在Mike21Exercises\BasicTemplate目录中.
参数信息见下表5.
1:表5.
1大肠杆菌模块参数信息一览表状态变量过程常数常数作用力辅助变量变量名FCFCdegrKcoliFtetaTtempfTemp描述FaecalColif-ormBacteriaDegradationofFaecalColiform1storderdecay20°C,(faecal)TemperaturecoefficientfordecayrateTemperatureFactorfortemperatureinfluence范围WCWC传输ADVECTION_DISPERSION空间变化HORIZONTAL_AND_VERTICALNONENONEHORIZONTAL_AND_VERTICALHORIZONTAL_AND_VERTICAL过程类型TransformationEUM类型Undefined1storderrateWQmodelUndefinedUndefinedUndefined单位1/100mlperdayDimensionlessdegreesCDimensionless用户自定义YesYesNo内置IDTEMPERATURE默认值00.
81.
0910最小值0010最大值1000000000101.
2100表达式-FCdegrKcoliF*FC*fTempPOW(tetaT,temp-20)输出YesYes18No.
3-BOD/DO和大肠杆菌练习目的:研究单一排污口对水中BOD、DO和大肠杆菌浓度的影响.
第一部分打开Mike21Exercises\DO-BODwithcoliforms\Masn\masn.
21,另存为masn_WQ1.
21.
文件内容:地形:bathy.
dfs2边界:流量west.
dfs0;水位east.
dfs0外力:wind.
dfs0水质模块:Wqsimple.
ecolab,仅简单定义了BOD/DO过程.
注意:当加载一个新的ECOLab模块时,常数、初始浓度值、外力和边界条件均设定为系统默认值,需根据当地实际情况和率定值进行进一步修改.
模型简介:1)该模型描述了两个岛屿之间的海峡区域,覆盖16.
5*7.
2km的范围.
2)有一个西边界和两个东边界,西边界为随潮汐变化的流量边界,东边界为水位边界.
3)在点(22,13)处有个污染源排放口,连续进行排放.
排放量为1.
0m3/s.
参见图5.
3.
19图5.
3参数设定:1)网格间距:300m2)时间步长:240s,更新频率:1个(ADTimeStep)3)时间范围:0~18004)初始条件:DO=10mg/l,BOD=0mg/l5)开边界:DO=10mg/l,BOD=1mg/l6)污染源(22,13):Q=1.
0m3/s,V=0.
1m/s,排放角度:0DO=0mg/l,BOD=50mg/l7)扩散系数:选择"independentofthecurrent",X方向5.
4m2/s,Y方向5.
4m2/s8)"Deposition"选项在Mike21中无需设定9)ECOLab模块参数设定参见下表5.
210)作用力设置:温度T=10℃,盐度S=20psu,其余均采用Built-in值.
饱和溶解氧=F(温度,盐度)11)输出时间步长:2hr20表5.
2练习No.
3中ECOLab模块的参数设定值编号常数描述数值1LatitudeBuilt-in2BODProcesses:1storderdecayrateat20deg.
celcius1.
5/day3BODprocesses:Temperaturecoefficientfordecayrate1.
0474BODProcesses:Half-saturationoxygenconcentration2.
0mg/l5Oxygenprocesses:Maximumoxygenproductionatnoon,m22.
0gO2/m2/day6Oxygenprocesses:Secchidiskdepths2.
65m7Oxygenprocesses:Timecorrectionforatnoon0hours8Oxygenprocesses:Respirationrateofplants,m21.
0gO2/m2/day9Oxygenprocesses:Temperaturecoefficient,respiration1.
010Oxygenprocesses:Half-saturationconc.
forrespiration2.
0mg/l11Oxygenprocesses:SedimentOxygenDemandperm21.
5gO2/m2/day12Oxygenprocesses:TemperaturecoefficientforSOD1.
0713Oxygenprocesses:Half-saturationconc.
forSOD2.
0mg/l结果查看和分析:提取A(16,13)和B(31,15)两点进行分析.
这两点的位置参见图5.
3.
可以直接在结果文件中输出这两点的结果或者根据.
dfs2结果文件提取A点和B点BOD和DO浓度的时间序列结果文件.
文件命名为x16_y13_10_20.
dfs0和x31_y15_10_20.
dfs0.
使用PlotComposer对结果进行绘制.
小结:要想获得合理的水质结果,边界条件对模型设置起到控制作用.
经过很短时间的计算,污染物浓度值即可趋于稳定,因此初始值对模型结果的影响较小.
BOD和DO浓度的变化遵循潮汐流的变化规律.
21第二部分从ECOLab模板中选择MIKE21/3WQSimpleincludingColiforms,按照第一部分的练习更新BOD/DO各参数值.
创建masn_WQ2.
21模型文件.
大肠杆菌组分参数设定:初始值:FC=0;TC=0;边界:FC=104个/100ml;TC=105个/100ml;在点源(22,13)处,FC=105个/100ml;TC=106个/100ml;使用系统默认值模拟粪大肠杆菌和总大肠杆菌的浓度分布.
此时,A(16,13)和B(31,15)两点的浓度值为多少第三部分基于第二部分练习的模型设置,进行灵敏度分析:1.
水温设为20℃,模拟在水温较高情况下BOD和大肠杆菌的变化.
结果与第二部分练习有何不同2.
模拟水温为20℃,盐度为30PSU的情况下,模拟BOD和大肠杆菌的变化.
此时,结果又有何不同22No.
4-使用ECOLab模拟盐度和CODcr练习目的:利用ECOLab-AD耦合的水质模块,对研究区域内的盐度和CODcr进行模拟,观测其模拟结果.
盐度和CODcr的水质模块使用练习No.
1中创建的新模块S_CODcr.
ecolab.
实际项目中往往利用盐度的率定来确定扩散系数的取值,这对后续水质模型的率定相当重要.
打开Mike21Exercises\MIKE21Bohai\BigModel\Bigmodel.
m21,另存为Bigmodel_Salinity_CODcr.
m21,进行模型设置.
Step1:导入新创建的盐度和CODcr的ECOLab模块-S_CODcr.
ecolab.
Step2:在MIKE21中进行参数设定,参数设定如下所述:1.
ECOLab更新频率:1个ADTimeStep2.
初始值:Salinity=15psu,CODcr=10mg/l,然后再提取结果文件以.
dfs2文件形式作为初始值,进行第二次模拟;3.
边界:Salinity=15psu;CODcr=10mg/l.
4.
扩散系数:X方向=200m2/s,Y方向=200m2/s5.
污染源Source:坐标(390,275),Salinity=0psu;CODcr=100mg/l6.
Deposition:不考虑7.
衰减系数:Salinity=0;CODcr=0.
3perday8.
结果输出:间隔一小时输出一个水质模拟值.
Step3:运行模型和结果分析.
比较不同类型初始值(常数和dfs2文件)对模型结果的影响.
23No.
5-暴雨污水溢流和氮过程模拟练习目的:该练习主要是研究水质模型中的氮循环过程.
详细介绍氨化作用、硝化和反硝化过程对N的影响.
打开Mike21Exercises\Nitrogenprocesses\Masn\masn.
21,另存为masn_WQnutrients.
21.
在ECOLab模板中选择MIKE21/3WQwithnutrients,根据练习no.
3中的参数值设定BOD/DO各参数.
注意:BODsediment项与传输扩散无关,其Transport项为NONE.
无需在边界条件中设定.
模拟时间:1987年9月1日00:00:00-9月6日00:00:00,时间步长240s.
初始值采用系统默认值.
边界水质浓度:BOD=0mg/lDO=10mg/lNH4-N=0.
5mg/lNO3-N=1.
0mg/lPO4-P=0.
1mg/l暴雨污水溢流发生在模拟时段的24hr后,即从1987年9月2日00:00:00开始.
发生地点在点(22,13)处.
整个污水溢流的持续时间为225分钟.
溢流详细信息如下:污水排放量=80m3/s,流速4m/sBODdiss=50mg/lBODsusp=50mg/lDO=2.
0mg/lNH4-N=5mg/lNO3-N=0mg/lPO4-P=1mg/l24根据以上信息创建污水溢流事件的污染源HD的时间序列文件:HD-spill.
dfs0,包括流量和流速两项,时间步长为15分钟.
在MIKE21中时间要选择timeequidistantcalendaraxis.
第一部分在水质参数对话框中将硝化反应和反硝化反应一级衰减速率均设成0.
硝化反应的温度系数设为1.
02,运行模型.
观察污水溢流前后DO、氨氮和硝酸盐氮的变化情况.
选择点(14,13)进行分析.
第二部分将硝化反应一级衰减速率设成2.
0/day,观察污水溢流前后DO、氨氮和硝酸盐氮的变化情况.
选择点(14,13)进行分析.
第三部分将硝化反应一级衰减速率设成2.
0/day,反硝化反应一级衰减速率设成1.
5/day,观察污水溢流前后DO、氨氮和硝酸盐氮的变化情况.
选择点(14,13)进行分析.
25No.
6-水质率定学习目的:以污水溢流和氮过程为例,进行水质模型率定,使模拟值与实测值达到较好的吻合.
该练习主要对水质模型中的BOD和氮过程进行率定,主要包括氨化、硝化和反硝化过程.
打开Mike21Exercises\WQCalibration\Lake.
m21HD模块地形:Lake-bathy.
dfs2时间步长:30s模拟时间:10天污染源溢流点(1,1):是由于化肥厂操作失误造成外泄,参见"spill.
dfs0".
WQ模块选择MIKE21/3WQwithnutrients,污染物组分浓度参见表5.
3.
表5.
3各污染物组分浓度一览表污染物组分[mg/l]BOD-dissBODsuspBODsedDONH4NO3P初始值/边界条件000100.
010.
10.
01溢流浓度300-0100.
13参数设置如下所述:AD扩散系数:1m/s(x和y方向)水温T=20°C盐度S=0其他参数选择缺省值,运行模型.
26水质监测数据在点(4,4)处有水质的实测值.
因此,通过比较该点的实测值与模型计算值,进行水质模型率定.
实测值的水质数据见表5.
4.
可以直接观察.
dfs0结果文件或者使用plotcomposer将实测值与模拟值放在一起进行比较.
提示:先观察BOD及其降解系数,然后再观察氮的数据变化.
表5.
4水质监测站(4,4)的水质指标一览表时间BOD-disBODsuspBODsedDONH4NO3PSt(4,4)02/01:004.
7004.
10.
681.
00.
49St(4,4)02/01:125.
3001.
70.
451.
40.
56St(4,4)03/01:125.
0000.
880.
151.
50.
54St(4,4)07/01:003.
7001.
660.
0151.
10.
46St(4,4)11/01:122.
4002.
50.
0100.
680.
38降解系数的首次选择:BOD-0.
05/day氨氮-3/day硝酸盐氮-0.
05/day27No.
7-渤海湾水质练习1、大、小模型地形900m网格大模型和300m网格小模型的研究区域见图5.
4.
图5.
42、大模型水质模型的建立水质模型设置信息等如下所述:模拟时间:2005-01-0100:00至2005-01-0200:00时间步长:100s,时间步数:864Source:坐标(390,275)Discharge:100m3/s;Velocity:0.
5m/s;OutletDir.
:45DO=0mg/l;Ammonia=25mg/l;Nitrate=10mg/l;BOD=30mg/l初始条件:DO=10mg/l;Ammonia=1mg/l;Nitrate=1mg/l;BOD=2mg/l边界条件:DO=10mg/l;Ammonia=0.
05mg/l;Nitrate=0.
25mg/l;BOD=0.
5mg/l28扩散系数:X方向=200m2/s;Y方向=200m2/sECOLab:引用WQlevel4.
ecolab,并将温度TEMP从StateVariable变成Forcing.
主要参数设定见表5.
5,其他参数设定均采用默认值.
表5.
5渤海湾水质练习参数设定一览表编号参数数值单位5OxygenProcesses:No.
ofreaerationexpression2-11Degradation:1.
orderdecayrateat20deg.
C0.
1/day31NitrogenContent:RatioofammoniareleasedatBODdecay0.
1gNH4/gBOD35Nitrification:Ammoniadecayrateat20degCelcius0.
05/day40Denitrification:Denitrificationrate,conversionofnitrateintofreenitrogenN20.
05/dayForcing:温度T=6℃,盐度S=5psu结果输出:间隔一小时输出一个水质模拟值,文件命名为result900m-WQ.
dfs2.
RunStartSimulation,运行模型.
3、小模型水质模型的建立Step1:提取小模型水质初始值从900m网格大模型结果文件中提取300m网格小模型的初始值.
提取大模型dfs2结果中最后一个时间步长的各水质浓度值作为小模型初始值.
(1)打开大模型水质结果文件result900m-WQ.
dfs2ToolsCrop…,提取最后一个时间步长的水质浓度值,生成WQini900m.
dfs2文件.
29(2)将WQini900m.
dfs2转换成小模型区域的初始值文件WQini300m.
dfs2NewMikeZeroToolboxTransformationRotateGrid,生成WQini900mTo300m.
MZT工具文件.
30选择结果输入文件WQini900m.
dfs2:SelectItem包括:溶解氧、氨氮、硝酸盐氮和BOD.
31参数设定如下:输出文件:WQini300m.
dfs232Step2:提取小模型水质边界小模型研究范围和坐标见图5.
4.
小模型位于大模型内部,小模型的原点在大模型(265,152)处.
从大模型的计算结果提取小模型计算所需的四个二维水质开边界.
以西边界为例:NewMike21ToolboxEnvironmentalTransferBoundary,弹出下面对话框,生成Extractboundaryfor300m.
21t工具文件.
33定义WEST边界条件:34引入大模型的水质结果文件:边界参数设置35Boudnary:WestLine:265,152(小模型在大模型中的原点)Angle:90(边界线的角度,顺时针到X正向的夹角)NumberofGridPoints:450(小模型边界的网格数)Gridspacingalongline:300(小模型的X空间步长)Gridspacingperpendicular:300(小模型的Y空间步长)36指定提取的时段(小模型模拟时段1月1号-1月2号)指定输出文件名West_WQ.
dfs1及其存储路径,并执行:37其他三个边界,做法同上.
四个边界的提取参数列表如下:边界名称边界起点旋转角度网格数X空间步长Y空间步长输出文件West265,15290450300m300mWest_WQ.
dfs1South265,1520550300m300mSouth_WQ.
dfs1North265,3020322300m300mNorth_WQ.
dfs1East448,15290274300m300mEast_WQ.
dfs1Step3:创建小模型水质模型水质模型设置信息等如下所述:模拟时间:2005-01-0100:00至2005-01-0200:00时间步长:100s,时间步数:864Source:坐标(375,380)Discharge:100m3/s;Velocity:0.
5m/s;OutletDir.
:45DO=0mg/l;Ammonia=25mg/l;Nitrate=10mg/l;BOD=30mg/l初始条件:从大模型提取的dfs2文件边界条件:从大模型提取的dfs1文件扩散系数:X方向=80m2/s;Y方向=80m2/sECOLab:参数设定同大模型.
主要参数设定见表5.
5,其他参数设定均采用默认值.
Forcing:温度T=6℃,盐度S=5psu结果输出:间隔一小时输出一个水质模拟结果,文件命名为result300m-WQ.
dfs2.
RunStartSimulation,运行模型.
386水质模型应用要点归纳1.
二维水质模型主要适用于局部河流、湖泊、河口和海洋等的水质模拟.
模拟结果主要取决于边界条件、污染源负荷以及外部作用力(如温度、太阳辐射、盐度等).
2.
若要降低初始值对模型计算结果的影响或继续下一模拟时间段的计算,可采用hotstart进行再次计算.
3.
WQ模块主要用于:水体中人为因素造成的污染负荷(点源和面源污染)占主体的情况,污染物在水体中停留的时间相对较短.
若对于污染物停留时间较长的湖泊水库,富营养化程度较大,则应该考虑用EU模型,考虑浮游生物的影响.
4.
增大水质计算时间步长,可以缩短计算时间,但也会增加模型发散的可能性,所以通过方案比较选择合适的时间步长进行计算.
初始选择时间间隔0.
5hr进行模拟计算.
5.
进行水质模拟前需进行对流扩散模拟,确定对流扩散系数以及检验模型的发散性.
AD模块可简单模拟盐水入侵、温度变化、污染物传输和一级降解.
完整的物理、化学和生物过程在ECOLab模块中实现,并与AD耦合使用进行计算.
6.
降解系数的单位在MIKE21AD模块中为/妙,在ECOLab中为/天.
7.
在MIKE21中,模型率定通常在模拟时间段内运行2次,第一次是将初始值设为常数进行模拟(一般采用起始模拟时间实测数据的平均值),第二次是将第一次运行后的最终二维结果平面图*.
dfs2文件作为初始值进行模拟,以减少结果对于初始值的依赖性和考虑不同研究区域内污染物浓度的梯度变化.
8.
对ECOLab中各常数进行率定时,一般采用同一参数值进行率定.
若研究区域内组分浓度梯度变化较大,则使用dfs2文件对不同区域的参数值进行分别设定.
9.
如果实测数据明显有问题,通过对各生态过程中的参数调整仍然不能达到较好的率定效果,则应该考虑根据实际情况对污染源的流量和负荷进行重新评估,作出适当修改后输入模型.
10.
造成模拟结果偏差的主要原因:39模型边界浓度定义不准确;输入生态系统的负荷定义不准确;初始浓度不准确;水动力模型不准确;选择的ECOLab过程描述不恰当;ECOLab过程表达式中的率定参数值不准确;实测值不准确;数值计算不稳定.
11.
率定耦合的常微分方程一般非常复杂,建议先从较独立的状态变量开始率定.
WQ水质模块建议的率定顺序:BOD-DO关系降解速率大气复氧产生和呼吸作用硝化营养物质与BOD的关系硝化/反硝化过程吸附作用EU富营养化模块建议的率定顺序:浮游植物潜在生长,硅藻和蓝绿藻的优势种群转化沉降消光作用温度系数调整营养物质腐质的沉降/降解40沉积(反硝化/固氮作用)氧大型植物12.
嵌套网格水质模型和基本MIKE21水质模型的区别和注意点:1)初始值应在每个area分别设定;2)扩散系数若取常数,则在整个areas适用,若为二维dfs2文件,则应在每个area分别设定;3)ECOLab模块的参数设定在全部areas保持一致;4)定义source时给出污染源在area1或area2中的坐标位置即可;5)输出结果和输出格式可以分别设置;6)使用嵌套网格计算速度一般比基本MIKE21模型低10%左右.