小浪底水库运用后三门峡水库运用方式中一些问题的初步研究.docVIP

小浪底水库运用后三门峡水库运用方式中一些问题的初步研究 目录 TOC \o 1-9 \h \z \u 目录 1 正文 1 文1：小浪底水库运用后三门峡水库运用方式中一些问题的初步研究 2 1 引言 2 310 9 310 9 312 9 312 9 315 9 315 9 312 9 319 9 315 10 319 10 5 三门峡水库汛期各级流量排沙效果研究 10 文2：小浪底水库运用后三门峡水库运用方式中一些问题的初步研究 13 1 引言 13 2 水库一维泥沙数学模型 13 3 模型验证 15 4 三门峡水库非汛期运用方式研究 16 5 三门峡水库汛期各级流量排沙效果研究 18 6 结语 20 参考文摘引言： 20 原创性声明（模板） 21 文章致谢（模板） 21 正文 小浪底水库运用后三门峡水库运用方式中一些问题的初步研究 文1：小浪底水库运用后三门峡水库运用方式中一些问题的初步研究 1 引言 小浪底水库投入运用后，三门峡水库所承担的 防洪 、防凌和调沙减淤任务将会有所变更，在新情况下，三门峡水库如何调度运用需要进行研究。在发生特殊的洪水或凌汛年份，三门峡与小浪底水库联合运用共同承担，这在小浪底水库设计中已有具体安排，故不再进行研究，本文主要就一般年份里三门峡水库运用的有关问题作初步研究。采用一维水动力学泥沙数学模型进行方案计算，并结合实测资料，对各方案库区的冲淤变化作综合分析。 2 水库一维泥沙数学模型 基本方程 水流连续方程 Q i+1 =Q i +Q l (1) 水流能量方程 Z i +α i V 2 i /2g=Z i+1 +α i+1 V 2 i+1 /2g+h w (2) 非均匀悬沙扩散方程 (3) 泥沙连续方程 (4) 式中 i 、 j 分别为过水断面序号及悬沙粒径组序号； Q l 、 h w 分别为相邻断面之间的入汇流量及水头损失；α为动能校正系数；β为泥沙恢复饱和系数； S j 、 S *j 分别为第 j 组悬沙含沙量、挟沙力； q sj 为单位流程汇入或流出的第 j 组悬沙输沙率。 异重流输沙计算 当异重流在底坡较陡的库底运行时，异重流的潜入条件可以采用下式判别 γ 0 /γ m -γ 0 ·V 2 0 /gh 0 = (5) 式中 V 0 、 h 0 分别为异重流潜入断面的平均流速和水深；γ m 、γ 0 分别为异重流和清水的容重。 当库底坡度很缓时，根据韩其为等人的研究 [ ２ ] ，异重流的潜入条件应为 h=max[h 0 ,h n ] (6) 式中 h 0 为按式 (5) 计算的潜入点水深； h n 为异重流均匀运动时的水深。 产生异重流后，从潜入断面开始要逐断面计算异重流的水力参数。在异重流运动过程中，因沿程流量减少不多，故在一个计算时段内可以作为恒定流来处理。异重流的阻力系数取平均值λ t= 。异重流水深计算选用异重流均匀运动水深计算公式进行，即 h n = （ λt/8gJ 0 ·γ 0 /Δγ·Qe/Be ） X)] 1/3 (7) 式中 J 0 为分河段底坡； Q e 为异重流流量； B e 为异重流宽度。 异重流挟沙力、含沙量以及冲淤量等的计算均与明流相同，只是此时的水力参数应采用异重流的。 求解过程 数学模型中一些基本问题的处理，如冲积河流的阻力变化关系、悬移质水流挟沙力、床沙级配计算以及冲淤量沿横断面的分布等已在 文献 [1] 中论述过，一维水流、泥沙控制方程组的离散和数值求解方法也详见文献 [1] 。值得指出的是，对于每一计算时段，当断面发生冲刷时，为了避免床沙级配计算中可能出现负值，特限制任一粒径组泥沙的冲刷量应小于活动层内床沙所能提供的量；当断面淤积时，规定每一时段的淤积总量应小于活动层面积，否则，通过缩短时段来满足。若相邻时段水流泥沙条件发生较大变化，如坝前水位骤降等，模型可以通过自动缩短计算时段来实现水沙条件的较平稳过渡，从而使溯源冲刷的模拟水平有了较大的提高。 3 模型验证 为了使所建模型适用于预测黄河中游大型水库库区泥沙冲淤过程和趋势，本文选用三门峡水库潼关以下库区 1960 ～ 1988 年共 28 年的实测水沙系列进行验证 计算 。三门峡水库三个不同运用期库区及分河段冲淤总量的验证结果见表 1 ，其中实测值为断面法测量结果。可以看出，各时段各河段计算的冲淤总量与实测值相差最大为 亿 m 3 ，一般为 亿 m 3 ，二者符合良好。累积冲淤过程以及每年汛期、非汛期的冲淤量也与实测结果符合较好，纵向冲淤分配在定性和定量方面与实际也较为接近。