三门峡水库运用方案模型计算.docxVIP

《三门峡水库给陕西带来的灾害及治理对策建议》附件ID 三门峡水库泥沙冲淤过程方案计算 西安理工大学 二OO—年十二月 专题名称：三门峡水库泥沙冲淤过程方案计算 课题 名称：三门峡水库给陕西带来的灾害及治理对策建议 专题编号：附件III 专题承担单位：西安理工大学水利水电学院 专题负责人：王新宏 冯普林 审 定：张润民 工作人员：王新宏 唐先海 张根广 杨武学 石长伟 马雪妍 冯普林 张强 报告撰写：王新宏 TOC \o 1-5 \h \z 摘要 1 1数学模型简介 1 1.1数学模型的控制方程式 1 1.2控制方程的离散及差分求解 2 2模型计算中有关问题的处理 4 1断面概化 4 2.2水面线计算的若干问题 5 3非均匀沙的沉速计算 7 2.4力匚流挟沙力级酉己 8 5床沙级配的调整及计算方法 8 3泥沙数学模型验证计算及分析 10 1验证依据的基本资料及有关问题的处理 11 3.2有关参数的率定 12 3 1960年7月?1965年10月验证结果与分析 13 3.4 1965年11月?1974年6月验证结果与分析 14 5 1974年7月?1990年10月验证结果与分析 16 3.6三门峡水库1960?1990年总体冲淤量验证结果 19 4方案计算及成果分析 20 1方案计算的目的 20 4.2计算条件及方案 21 3方案计算成果及分析 23 5结语 27 参考文献 28 三门峡水库泥沙冲淤过程方案计算 摘要：本文采用西安理工大学的准二维悬移质不平衡输沙数学模型，对三门峡水库4种 不同运用方案下潼关以下库区的泥沙冲淤演变过程进行了预测计算。预测计算结果表明, 三门峡水库采用4种推荐的运用方案情况下，水库运行25年后，潼关1000m3/s流量时 的水位均大幅下降，可比现状潼关高程328.47m降低3.31 ~ 4.30m,潼关断面达到冲淤平 衡约需8-12年时间。文中建议，为了使潼关高程在尽可能短的时间内降下来，三门峡水 库今后应考虑采用全年敞泄的运用方式。 1数学模型简介 西安理工大学泥沙数学模型属于准二维悬移质不平衡输沙数学模型，适 用于进行长历时长河段的河床变形研究。模型中将所研究的长河段划分成若 干小河段，可预测计算各断面的水力、泥沙要素以及相邻断面Z间的平均冲 淤厚度的沿程变化及因时变化情况。 1.1数歹儀型的滋刮方無式 本模型属于准二维恒定流非均匀悬移质非饱和输沙模型，可用来模拟水 深、流速、水面比降、悬移质含沙量及其级配、床沙级配等水力泥沙要素随 时间和沿流程的变化。模型所依据的基木控制方程有水流连续方程、水流运 动方程、泥沙连续方程、悬移质扩散方程、水流阻力计算公式和水流挟沙力 计算公式，其形式为： 水流连续方程式 ~c1l =°水流运动方程式dZ ( ~c1l =° 水流运动方程式 dZ ( d(V2} Q2 =0 (1-1) (1-2) 泥沙连续方程式（分粒径组） TOC \o 1-5 \h \z + 7 学+ %=0 伙=1,2???化） （1-3） OX Gt 悬移质扩散方程（分粒径组） 鱼二丝佝-SJ 伙=1,2 …,NJ （1-4） dx q 式中，0是断面平均流量，血是侧向入流量，Z是水位，乙是动能修正系数, 孑是局部水头损失系数，U是断面平均流速，K=ACR!/2是流量模数（其中，A 是过水断面面积，R是水力半径,C=Ri/6/n是谢才系数，“是糙率）,S、S*分 别是断而平均含沙量和水流挟沙力，以kg/n?计，厂是泥沙干容重，人是泥 沙冲淤面积，血是单位流程上第k粒径组泥沙的侧向输沙率（入为正，出为负）, P是粒径组角标，M是粒径分级总数，是恢复饱和系数，。是平均沉速，q 是单宽流量，g是重力加速度，x、f分别是流程和时间。 方程（1-1）?（1-4）的定解条件包括初始条件和边界条件。初始条件包括河 道初始地形和床沙级配。边界条件包括上游入口断而潼关站的水沙过程，河 道沿程水沙汇入（或流出）过程，以及下游出口断面史家滩站 的水位过程。上 述4个基本方程式本身并不封闭，还需增加一些补充关系式，诸如水流挟沙 力计算公式，水流挟沙力级配计算公式，动床阻力计算公式等。 1.2控制方程的离散及差分求解 模型屮求解上述基本方程时，首先将研究河段沿流程划分为若干河段，使 每一河段内的水流接近均匀流，如图1-1所示。图中，i是断面编号，是 断面总数，断面序号由下游向上游递增。为简便起见，称第i断面至笫/+/ 断面之间的河段为第i河段。其次，把长历时的来水、来沙过程概化为梯级 恒定流，使每一计算吋段内的水流接近恒定流，按恒定非均匀渐变流进行求 解。 图1-1河段划分及横断面编号示意图 计算方法采用非耦合解法。具体作法是先解水流方程一一水力计算，即 利用明渠恒定渐变流运动方程，推求库区水面线，从而得到各断而