MIKE3-水质材料.doc

Four short words sum up what has lifted most successful individuals above the crowd: a little bit more.------------------------------------------author------------------------------------------dateMIKE3-水质材料MIKE 11 HD、Ecolab完全实例教程MIKE 3水质培训教程王 颖DHI ChinaInland Water DepartmentOffice number 021 6417 7291way@1 ECO Lab简介ECOLab是DHI在传统的水质模型概念发展起来的全新的水质和生态模拟工具ECOLab软件开发的理念和方法非常先进，用户不仅可以修改模型参数，更重要的是可以修改模型核心程序、甚至编写新程序，然后ECOLab将其与MIKE 11/21/3的HD、AD集成计算DHI已经将大部分传统的水质模块转换成ECOLab通用模板，供用户调用或修改使用，包括：v 水质模块v 富营养化模块v 重金属模块1.1 应用领域v 河流、湿地、湖泊、水库、河口、海岸和海洋v 各生态系统反应的空间预测v 简单和复杂的水质研究v 环境影响和优化研究v 规划和可行性研究v 水质预报1.2 内置模板和使用手册· DHI预定义的ECO Lab模板在以下目录中：C:\Program Files\DHI\2009\MIKE Zero\Templates\ECOLab· 使用手册和说明在以下目录中：C:\Program Files\DHI\2009\MIKE Zero\Manuals\MIKE_ZERO\ECOLab2 WQ－水质模块2.1 MIKE 3 WQ 水质模块的目标MIKE 3水质模型主要针对湖泊、海洋区域的污水排放引起的水质问题，比如BOD/DO, 富营养化和细菌污染。

2.2目前水质模块可进行以下模拟：· 大肠杆菌，粪大肠杆菌/总大肠杆菌的传输和死亡（用一级降解来表示），降解速率取决于当地的光强，温度和盐度条件等· BOD-DO关系，即排放的有机物所引起的耗氧考虑以下几个过程： BOD一级降解 BOD降解引起的耗氧 底泥需氧量 水体中的呼吸作用 光合作用产氧 水气相互作用下的氧交换 （大气复氧）· BOD-DO模块包括不同营养物（氨氮，硝酸盐和磷）以及三种BOD形式：溶解性，悬浮性和沉积性BOD使用该模块需要设置三种BOD组分的一级降解速率悬浮和沉积的BOD将考虑沉降和再悬浮该模块中氧平衡过程主要包括：BOD降解需氧量，底泥需氧量，硝化反应需氧量，光合作用产氧，呼吸作用耗氧以及大气复氧营养物转化的基本过程包括：BOD降解释放有机氮和磷，产生的氨氮经硝化反应变成硝酸盐氮，最终通过反硝化作用生成氮气，释放在大气中同时，BOD降解所释放的部分氨氮和磷可以被浮游生物，植物和细菌所吸收· 用户可以按实际需求自定义多种污染物质，并定义相应的降解速率进行模拟典型污染问题与典型污染问题相关的污染物质有：· 近海水域中与健康相关的微生物· 耗氧物质· 营养物质· 异生化合物，例如有危害性或毒性的化合物与健康相关的微生物对于近海水域微生物调查的主要目的在于指出其用水安全性，或是作为对该处鱼类，贝类等生长环境的调查。

一个全面的微生物风险评估包括：· 环境健康评估包括关于排水管道或污水排放口，雨水排放口的季节性变化，水温，流量，潮汐变化等信息，以及一个报告和行动系统以确保水质恶化引起的问题能及时通知到健康权威机构并作出相应处理· 指示剂生物体的出现和这些生物体的行为，包括其与物理－化学因素及相关病原体关联的死亡速率（基于光强、盐度、水温、沉降速率和污染程度等）· 病原体的呈现耗氧物质耗氧物质分为溶解性和悬浮性物质，与氧进行生物或生物化学作用，消耗水中的溶解氧这些耗氧物质主要是一些不同类型的有机物，具有不同类的降解速率生化需氧量（BOD）是间接反映水中能为微生物分解的有机物总量的一个综合指标有机物在有氧条件下为微生物分解产生H2O、CO2和NH3一般BOD以被检验的水样在标准条件下5天内的耗氧量为代表，称为BOD5营养物质许多营养物质都是生物生长的必要元素适量的营养物对于水中微生物的生长及活动是必需的，然而，一旦营养物质过量就会引起富营养化，将引起一系列的问题，如水体污浊，河床底部缺氧，生物沉积量的增加等富营养化模块可用来模拟这种情况，因为该模块考虑到藻类对其它物质的直接影响在营养物质中氮和磷是最重要的，它们是水生植物生长的控制因子。

氮以氨氮和硝酸盐这两种无机氮的形式存在许多国家对近海水域中的这些营养物质都设定了浓度标准MIKE 3 水质模型 (WQ) 就是设计用于评估和这些标准浓度相关的水质问题MIKE 3富营养化模块（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产率和被吸收速率· 反应级数· 硝化/反硝化速率主要参数经验值：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- 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.066 gNH4/gBOD6. 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/gBOD9. 颗粒态磷(PP)降解速率: 0.1-0.2/day 3 MIKE 3水质模型MIKE 3水质模型所需数据资料：· 基本模型参数：地形网格（结构网格图或非结构网格图） 时间步长和模拟时间 输出项类型和频率· 地形和HD条件· 耦合的AD模型：扩散系数的率定· 初始值：各参数的浓度值· 边界条件：各参数的浓度值· 污染源：坐标位置、水动力条件及各参数的浓度值· 各生物过程速率值：参考率定值、经验值或监测值等。

在MIKE 3模型中添加ECO Lab模块第一步：引入水质模块：1． 在MIKE 3结构网格模型中引入ECO Lab模块：打开 MIKE 3 à Flow Model (.m3) à Basic Parameters à Module Selection à Environmental modules à ECO Lab (图3.1)，出现ECO Lab Parameters，在Model Definition里选择适当的内置水质模块或自定义模块，参见图3.2 图3.1 图3.22． 在MIKE 3非结构网格模型中引入ECO Lab模块：打开 MIKE 3 à Flow Model FM à Module Selection à ECO Lab（参见图3.3），出现ECO Lab Module，在Model Definition里选择适当的内置水质模块或自定义模块，参见图3.4 图3.3 图3.4第二步：在Model Definition里选择计算方法并设置水质模拟时间步长。

水质模拟时间步长通常先以0.5小时进行计算第三步：在ECO Lab模块中分别对需模拟的状态变量、边界水质条件、扩散系数、污染源浓度、ECO Lab模块各参数、作用力和输出项进行设定4 水质模型应用要点归纳1. 三维水质模型主要适用于水库、湖泊、河口和海洋等的水质模拟模拟结果主要取决于边界条件、污染源负荷以及外部作用力（如温度、太阳辐射、盐度等） 2. 若要降低初始值对模型计算结果的影响或继续下一模拟时间段的计算，可采用hotstart进行再次计算也就是在模拟时间段内运行2次，第一次是将初始值设为常数进行模拟（一般采用起始模拟时间实测数据的平均值），第二次是将第一次运行后的如三维结果图*.dfs3或 *.dfsu文件作为初始值进行模拟，以减少结果对于初始值的依赖性和考虑不同研究区域内污染物浓度的梯度变化3. WQ模块主要用于：水体中人为因素造成的污染负荷（点源和面源污染）占主体的情况，污染物在水体中停留的时间相对较短若对于污染物停留时间较长的湖泊水库，富营养化程度较大，则应该考虑用EU模型，考虑浮游生物的影响4. 增大水质计算时间步长，可以缩短计算时间，但也会增加模型发散的可能性，所以通过方案比较选择合适的时间步长进行计算。

初始选择时间间隔0.5hr进行模拟计算5. 进行水质模。