MIKE21在水库氨氮运移过程中应用.pdfVIP

MIKE2 1在水库氨氮运移过程中的应用 张神铭徐征和李福杰杜敏徐秀凤 (济南大学资源与环境学院，山东济南250022) 调蓄水库进行了水流一水质的耦合模拟。以2008年3月为例分别模拟并分析了水库在引水排水、 引水不排水、不引水排水、不引水不排水四种不同的水库调度情况下水库流速场变化，选用氨 氮为主要控制指标进行水质模拟并分析水库氨氮的浓度场变化和氨氮的运移过程，可以直观观 测到氨氮在水库中的运移降解情况，为水库管理者对水库引排水调度和污染物的控制提供了决 策依据。 关键词：MIKE21；氨氮运移；济南市引黄供水鹊山调蓄水库 0前言 氨氮是引黄水库的主要污染物，也是水质模型研究的主要对象，利用计算机软件模拟的方 法，模拟水库氨氮运移过程，计算结果较为精确，提高了计算效率，经过实测值与模拟值的对 比，表明mike21软件自带的水质模型可以正确模拟水库的氨氮运移过程。影响氨氮运移降解的 因素有很多，利用mike21软件模拟水库在不同引排水情况下氨氮的运移降解过程，直观反映出 流速和引排水口氨氮浓度对降解过程的影响，可以直观模拟氨氮在水库中的运移情况，为决策 者对水库的调控提供了科学依据。 1软件介绍 用于模拟河流、 湖泊、 河口、 水库等地表水体流动和海湾、 海岸、 海洋的水流、 波浪、 泥沙运移及水环境变化n1。软件中的ECOLab是一个完备的、用于生态模拟的数值实验室，它 提供了从简单到复杂的解决方案。而且，还提供了一系列的模板，用户可根据自己的具体应用 选择使用模板。模型可以描述物理沉降过程，也可以描述化学、生物、生态过程以及状态变量 之间的相互作用。一个ECOLab模板包括生态系统的数学描述或者一系列通用的微分方程。 ECO Labqb这个数学描述分为6个部分，分别是状态变量、常量、作用力、辅助变量、过程、 衍生结果。 2模型构建 2．1水库概况 济南市引黄供水鹊山调蓄水库位于黄河北展区下首，属于黄河下游冲积平原。东l晦黄河 大堤，南靠原津浦铁路，北靠黄河北展宽大堤，平面上近三角形，库区内地势平缓，地面高程 一般为23．5—24．5米，平均库底高程22．4米。由于人工取土筑堤、挖塘养鱼以及地表河流的冲 刷，致使四周及中部河沟处地势低洼，芦苇塘、沼泽地随处可见，区内最低位于大吴排水闸出 口，高程为22．4米，在库区周围，北展堤堤顶高程32．6—4．0米，高出地面9米。黄河大堤堤顶 高程约36米，临河滩面29米，高出背河地面3—5米，津浦铁路路堤堤顶高程约26米，高出地面 2-3米。 水库经l#泵站提水进入沉砂条渠，通过轴线长度为3．2km的沉砂条渠、人库输水涵洞送 至2#泵站前池，再由2#泵站提升到调蓄水库，最后由3#出库泵站增压经输水钢管从河底穿 越黄河送至黄河地面水厂，规模为44万m3／dt21。 2．2地形构建 地形前处理是水动力学数值模拟的重要内容之一，也是水动力学数值模拟快速发展的一个 重要环节[3】。济南市鹊山水库地形的确定，是根据现有资料，确定水库的经纬度坐标，经转化 以后，导入软件中，经过反复调试使网格大小尽量均匀，考虑到网格大小、网格数量与计算的 拟区域最大网格边长不超过200米，最小网格边长不小于80米，在生成模拟区域的进口和出口 添加线段来定义进口边界和出口边界。 插值是地行前处理的中不可少的步骤，从模型离散格式上看，地形前处理的插值方法主要 有结构网格法、非结构网格法和非结构、有结构网格相结合发。MIKE21模型所涉及的属于结 构网格【4】。根据资料，选取不同位置的13个点的水深值进行了内插，得到水库地形图。 2．3有关问题处理 2：3．1糙率的选择 水库的糙率与水库形态，库底的粗糙情况、植被生长状况、水位的高低，库区的冲积以及 整治的人工建筑物等诸多因素有关。糙率值n的大小或变化规律，会影响到方程中阻力项的大 小，从而对计算结果带来直接影响[5】。由于缺少实测资料，这里采用软件中默认的一般的经值 32／s，软件中所提到的曼宁系数和其他文献提到的曼宁数为倒数关系。 2．3．2时间空间步长的选取 离散格式的选定以离散后的方程与原方程相容为前提，且要求其收敛性和稳定性均能够得 到满足，而具体如何合理选取空间时间步长，对计算能否顺利进行将有直接影响[5】。此处经过 调试， 我们选取的步长为10800s，模拟步数为80。 与2008年3月