后5-徐老师--三峡水库调峰运行下的香溪河水质数值模拟.ppt

 三峡水库调峰运行下的香溪河水质数值模拟 主要内容 1 前言 香溪河是长江一级支流，也是三峡水库库区最大的支流，河长106 km，总流域面积3099 km2，其河口距三峡水库大坝约30km。 随着三峡水库蓄水，其下游近河口约20km的河段水位随之升高，水环境由典型的河流水体转变为类似湖泊的缓流水体，水体理化指标变化较大，水质问题相对较为严重，多次出现大面积水华现象，造成沿岸部分居民饮水困难。 2 基本控制方程 2.1 平面二维水流方程 水流连续方程： 2 基本控制方程 2．2 一维水质方程 式中：h为水深，m；u和v分别为x和y方向流速，m/s；g为重力加速度，m/s2；zb为床面高程，m；ρ为水密度，kg/m3； τsx和τsy分别为表面切应力的x和y方向分量，N/m2； τbx和τby分别为底部切应力的x和y方向分量，N/m2；εij为涡粘系数；f为科氏力系数；c为污染物浓度，mg/m3 ；A为河道过水断面面积，m2 ；Q为河道流量，m3/s；Ex为纵向离散系数，m2/s；dc/dt为生化反应项，mg/（m3s）；Sext为由于外源在单位时间单位体积内增加或减少的污染物质量，mg/（m3s）。 3 数值离散求解方法 平面二维水流方程在空间上采用有限元加权余量法离散，此方法网格划分灵活，适于处理复杂多变的边界地形条件。时间上采用隐式非线性有限差分方法离散，以保证计算稳定性和精度。对控制方程进行离散后得到非线性代数方程组，采用Newton-Raphson迭代方法求其数值解。 一维水质方程采用有限差分法中的六点交替隐式格式离散。此离散格式为时间向前插分，空间中心差分，具有二阶精度。离散后采用追赶法（双扫描法）求解。 4 模拟区域 5 计算结果与分析 应用上述水质模型，计算了水库不同谷荷流量、调峰流量、水位及污染物浓度对香溪河水质的影响。通过计算得到香溪河平均污染物浓度及水质速率变化规律。 所谓水质速率定义为初始平均浓度与一天后平均浓度的差值。 水质速率为负时，表示污染物浓度增加，水质恶化；水质速率为正时，表示污染物浓度减小，水质得到改善。 5.1 不同坝前水位对浓度场的影响 5.1 不同坝前水位对浓度场的影响 5.2 不同谷荷流量对浓度场的影响 5.2 不同谷荷流量对浓度场的影响 5.3 不同调峰流量对浓度场的影响 5.3 不同调峰流量对浓度场的影响 6 结语 Tianjin University 天津大学建筑工程学院 徐国宾 教授 2 基本控制方程 1 前言 3 数值离散求解方法 4 模拟区域 5 计算结果与分析 6 结语 水流运动方程： 图1 模拟区域网格观测点布置 模拟区域为香溪河汇流口至上游25km，三斗坪坝址至长江干流上游40km 范围内。计算区域概化以1：5000的精细河道地形图为基础，根据河道边界及 地形的特点，划分二次四边形单元共8340个网格。网格及观测点布置如图1。 175 165 155 145 0 0 0 0 10 0.03930 0.03945 0.03971 0.03914 8 0.09824 0.09862 0.09927 0.09786 5 0.13754 0.13806 0.13898 0.13700 3 0.19649 0.19724 0.19854 0.19571 0 10000 5000 水位（m） 浓度 （mg/L) 调峰流量 （m3/s） 谷荷流量 （m3/s） 表1 不同坝前水位下的香溪河水质速率（单位：mg/L） 计算结果表明：在调峰运行方式相同情况下，不同的坝前水位对香溪河水质速率影响 不大。不同水位下香溪河的总体污染物浓度约有0.002mg/L的差别，对香溪河水质改善 不明显。同时可以看到，干流污染物浓度边界对污染物浓度变化的影响比较大，当干 流污染物浓度边界相差3mg/L时，香溪河水质速率相差0.06mg/L左右，见图2。 图2 不同干流浓度下香溪河水质速率变化 0.03648 0.03630 0.03582 0.03522 8 0.09121 0.09075 0.08954 0.08806 5 0.12769 0.12705 0.12535 0.12328 3 0.18242 0.18149 0.17908 0.17611 0 175 0.03662 0.03641 0.03597 0.03548 8 0.09155 0.09102 0.08992 0.08870 5 0.12817 0.12743 0.12589 0.12418 3 0.18310 0.18205 0.17985 0.17740 0 165 0.03661