洱海水质预测模型的研究.pdfVIP

洱海水质预测模型研究 彭文启 (中国水利水电科学研究院水环境研究所，北京。100038) 摘要本文在二维水动力学基础上提出了洱海水质预测模型．论文系统地模拟洱海 湖流运动，水动力学模型能较好的反映洱海湖流特征，与已有的测验资料基本吻合．论文 构建三类水质项目COD、TN和TP的水质模型，并对2000年的水质监测资料进行验证， 吻合较好．洱海水质预测模型可以较好地揭示洱海水质时空变化规律，为洱海水污染控制 规划提供了一个较好的平台． 关键词瑚泊，水质，横墨} 1洱海流域自然环境概况 洱海位于大理市境内，是云南省第二大高原淡水湖泊。洱海是大理市主要饮用水源地。 又是苍山洱海国家级自然保护区和风景名胜区的核心，具有调节气候，提供工农业生产用 水、水生生物多样性等多种功能，是整个流域乃至大理州经济社会可持续发展的重要基础， 堪称大理人民的“母亲湖”． 流域面积2 年输送水量8．13×108m3，湖水停留时间75年。 洱海上游有茈碧湖、海西海、西湖等高原断陷湖泊，区内河流河网较多，入湖河流大 小拭117条。洱海北承弥苴河注入，西纳苍山十八溪。东收波罗江、风尾阱等溪水，西出 洱河汇入黑惠江，最后汇入澜沧江。 2洱海水动力学模型控制方程 湖泊湖流运动是湖泊中污染物质扩散、沉积物再悬浮与悬移质输送等过程的主要水动 力因素，是湖泊水体水质时空分异特征形成的主要参与者，因此也是湖泊水质管理的重要 ·48· 科学依据。 现有研究成果表明，洱海的湖流研究仍然停留在定性描述阶段，主要原因是缺少湖流 监测资料。湖流监测工作量十分庞大，且在监测断面设置、湖区风情同步监测、监测数据 的统计分析等方面存在诸多技术难题，因此国内进行过湖流监测的湖泊屈指可数，这也是 中国湖泊研究的诸多不足之一。 深水湖的一个重要物理现象是温度分层，温度分层的一个显著标志是温跃层的出现。 研究成果表明，中等维度的湖泊的季节分层深度间于5-15m，形成于春季或初夏，并维持 至秋季。根据温跃层深度也可以界定“深水湖”和“浅水湖”，当湖泊水深低于正常季节温 跃层深度时，可以称为“水力浅湖”，或简称“浅水湖”。如果属于“浅水湖”，湖泊水体垂 向温度分布接近均匀，其水动力模拟研究，可以不考虑分层，如果出现温度分层，则最好 采用分层模型。 最大水深20．9m，已经超过区分深水湖和浅水湖的4m水深界限，因此从水深角度分析， 洱海应该属于深水湖。但根据几次湖泊水温监测研究成果[1，2]，未能发现明显的温度分 层现象，因此本次洱海水质预测采用平面二维模型。 2．1二维水动力学控制方程 水动力学方程采用水深平均的二维非恒定模型，紊流模型采用混长模型。水流与水质 模型方程采用守恒形式，控制方程连续方程和水平向动量方程[3]： 誓+∥挚+∥孥：一∥善．晏(。冬)+昙(。誓)一 讲 ox 印 ox出出 却 印 r—。—-‘ 。 。 f11 g迪高攀+乞肌厢 、。 挈ot+∥孥ox+卢鲁=一妒筹+导ox卵 卯 cs挈ox，+熹@等卜泖 钟 g巫高攀+去肌厢 (2) 等+等专咆 (3) 式中，U、v为水深平均流速分量，q，=uh，qr--vH为Ⅳ、Y向单宽流量，f为水位，H 为总水深，口为水平向流速垂宣分布非均匀分布修正系数，g为重力加速度，C为谢才系 ·49· 数，声风摩阻系数，¨一，Vw分别为X、y向风速度。Q。为源汇项，主要用来模拟湖泊出 入流。f为水深平均涡黏系数，由混长紊流模型计算。 2．2二维水质模型控制方程 湖泊二维非恒定水质模型控制方程为： 警+警+等=去c破日》叫C。(DyHa吵S)+S。+HD。S， 式中，s为水质项目水深平均浓度， 百dS为生化反应项，