三峡水库整体一维水质数学模拟研究.docVIP

三峡水库整体一维水质数学模拟研究 李锦秀1，廖文根1，黄真理2(1.中国水利水电科学研究院 水环境研究所;2.国务院三峡工程建设委员会办公室) 摘要：以三峡整个库区作为研究对象体，根据三峡水库河道型特点，开发了一维非恒定水流水质数学模型，利用稳定性比较好的四点隐格式有限差分法进行模型数值求解。以主要污染指标BOD5、DO为例，介绍了三峡水库一维水流水质数学模型的开发、研制过程。模型中充分考虑了三峡水库建成前后水流条件巨大变化对库区水流水质的影响，对主要水质模型参数，如污染物纵向离散系数、有机污染物衰减系数和大气复氧系数等建立了与水流条件相关的经验关系式，建立的数学模型经过三峡库区江段长河段水文水质和污染源同步观测资料验证，具有较高的模拟精度，可为预测三峡水库水流水质总体变化趋势提供技术支持。关键词：三峡水库；水质模型；水质参数；经验公式中图分类号：X11? ??????????? 文献标识码：A 三峡大坝建成以后，库区水质将如何变化一直是社会各界广泛关注的环境问题之一。三峡库区从库首三斗坪至库尾重庆全长约650km左右，水库建成后，尽管水位抬高、河面加宽，但水库蓄水后仍属于典型的河道型水库(如图1所示)，库区沿程总体水流、水质运动特性仍符合一维水流水质运动规律。为此，作者进行了三峡库区一维水流水质数学模型的开发研究，为模拟预测三峡水库建成以后库区不同江段平均水质变化趋势提供技术支持。　　一维水流水质数学模型开发研制的关键技术是模型参数的精确取值问题，尤其是三峡水库建成前后水流条件将发生巨大变化，开发研制的数学模型能否预测水库建成以后的水流水质状况，是数学模型开发研究的难点。本文将以三峡库区主要污染指标BOD5、DO为例，介绍三峡库区一维水流水质数学模型的开发研究成果。 图1 三峡库区河道示意 1 三峡库区一维水质数学模型　　描述物质在水体中输移离散的一维不恒定流基本运动方程为： (1) 式中：C为污染物浓度；A为过水断面面积；E为纵向离散系数；Q为断面流量；dC/dt为生化反应项；S为源、汇项。　　在对方程(1)进行数值求解时，采用稳定性比较好的四点加权隐格式进行数值离散，用双追赶法求方程收敛解。　　根据三峡水库实际污染情况，水质模型主要水质模拟指标确定为：DO、BOD5、CODmn、NH3-N、NO2-N、NO3-N、TP、TN、悬浮物(SS)、酚、汞、铬、油、水温等14项1。本文将以三峡库区主要水质指标BOD5、DO为例，介绍三峡水质数学模型的建立。根据三峡水库实际情况，BOD5、DO生化反应项表示如下： (2) =K2(Cs-CD)-K1CB-α5KN1CN1-α6KN2CN2? (3) 式中：CB为BOD5的浓度；K1为BOD5降解速率常数；K2为大气复氧系数；CD为溶解氧浓度；CS为饱和溶解氧浓度；CN1、CN2分别为氨氮和亚硝酸氮浓度；α5为单位氨氮氧化时的耗氧率；KN1为氨氮氧化速率，α6为单位亚硝酸盐氮氧化时的耗氧率；KN2为亚硝酸盐氮氧化速率。　　由于三峡库区天然河道中藻类含量比较低，简单分析计算表明，天然库区江段光合作用产氧量和呼吸耗氧量与大气复氧量和有机污染物生化耗氧量相比非常小，可以忽略不计。三峡水库建成以后，随着库区流速减缓，泥沙将发生淤积，库区水体透明度会比建库前天然河道状况增大，浮游植物光合作用产氧量也有可能增大。由于藻类光合作用产氧量与光照强度、水体透明度、温度和藻类含量等因素有关，十分复杂。受资料等多种条件限制，本文从偏于保守考虑，在方程(3)的溶解氧生化反应项中忽略了浮游植物光合作用产氧和呼吸耗氧过程。1三峡水库一维水流水质数学模拟研究，中国水利水电科学研究院水环境研究所，2001年8月2 一维水质模型参数研究　　水质模型中各类模型参数通常代表着污染物质在水体中发生的物理、生物和化学反应过程。受诸多因素的影响，水质模型参数精确取值一直成为水质模型的开发难点。尤其是三峡水库建成前后水流条件将发生巨大变化，必将影响水质的物理和生化过程。因此，水质模型参数精确取值是三峡水库水质数学模型开发研制所需解决的关键技术。以水质指标BOD5-DO为例，一维BOD5、DO水质模型中主要待定水质参数包括纵向离散系数E、BOD5衰减系数K1和大气复氧系数K2，关于氨氮和亚硝酸氮氧化过程模拟研究成果的介绍在本文中从略。2.1 纵向离散系数 纵向离散系数是代表河流纵向混合特性的重要参数。三峡库区江段属山区型河道，江段横向和纵向断面形态十分复杂，水流的纵向离散特征应该比较显著。作者通过理论分析，结合三峡库区实际河道特征建立的三峡库区纵向离散系数计算公式如下，有关的研究成果见参考文献[1]。 (4) 2.2 BOD5衰减系数K1 BOD5衰减系数反映了有机污染物质复杂的耗