水利工程论文-小浪底水库运用后三门峡水库运用方式中一些问题的初步研究.doc

水利工程论文-小浪底水库运用后三门峡水库运用方式中一些问题的初步研究摘要：小浪底水库投入运用后三门峡水库的运用方式将作调整，本文运用黄河水库一维水动力学泥沙数学模型通过典型条件下的方案对比计算与分析，从减少库区淤积及增加电站出力两方面综合考虑，认为水库非汛期运用的第4方案较好；从排沙效果来看,汛期水库降低水位开始排沙的入库流量应选择大于1500m3/s为宜。关键词：泥沙数学模型水库运用方式方案北村水位1引言小浪底水库投入运用后，三门峡水库所承担的防洪、防凌和调沙减淤任务将会有所变更，在新情况下，三门峡水库如何调度运用需要进行研究。在发生特殊的洪水或凌汛年份，三门峡与小浪底水库联合运用共同承担，这在小浪底水库设计中已有具体安排，故不再进行研究，本文主要就一般年份里三门峡水库运用的有关问题作初步研究。采用一维水动力学泥沙数学模型进行方案计算，并结合实测资料，对各方案库区的冲淤变化作综合分析。2水库一维泥沙数学模型2.1基本方程水流连续方程Qi+1=Qi+Ql(1)水流能量方程Zi+iV2i/2g=Zi+1+i+1V2i+1/2g+hw(2)非均匀悬沙扩散方程(3)泥沙连续方程(4)式中i、j分别为过水断面序号及悬沙粒径组序号；Ql、hw分别为相邻断面之间的入汇流量及水头损失；为动能校正系数；为泥沙恢复饱和系数；Sj、S*j分别为第j组悬沙含沙量、挟沙力；qsj为单位流程汇入或流出的第j组悬沙输沙率。2.2异重流输沙计算当异重流在底坡较陡的库底运行时，异重流的潜入条件可以采用下式判别0/m-0V20/gh0=0.6(5)式中V0、h0分别为异重流潜入断面的平均流速和水深；m、0分别为异重流和清水的容重。当库底坡度很缓时，根据韩其为等人的研究，异重流的潜入条件应为h=maxh0,hn(6)式中h0为按式(5)计算的潜入点水深；hn为异重流均匀运动时的水深。产生异重流后，从潜入断面开始要逐断面计算异重流的水力参数。在异重流运动过程中，因沿程流量减少不多，故在一个计算时段内可以作为恒定流来处理。异重流的阻力系数取平均值t=0.025。异重流水深计算选用异重流均匀运动水深计算公式进行，即hn=（t/8gJ00/Qe/Be）X)1/3(7)式中J0为分河段底坡；Qe为异重流流量；Be为异重流宽度。异重流挟沙力、含沙量以及冲淤量等的计算均与明流相同，只是此时的水力参数应采用异重流的。2.3求解过程数学模型中一些基本问题的处理，如冲积河流的阻力变化关系、悬移质水流挟沙力、床沙级配计算以及冲淤量沿横断面的分布等已在文献1中论述过，一维水流、泥沙控制方程组的离散和数值求解方法也详见文献1。值得指出的是，对于每一计算时段，当断面发生冲刷时，为了避免床沙级配计算中可能出现负值，特限制任一粒径组泥沙的冲刷量应小于活动层内床沙所能提供的量；当断面淤积时，规定每一时段的淤积总量应小于活动层面积，否则，通过缩短时段来满足。若相邻时段水流泥沙条件发生较大变化，如坝前水位骤降等，模型可以通过自动缩短计算时段来实现水沙条件的较平稳过渡，从而使溯源冲刷的模拟水平有了较大的提高。3模型验证为了使所建模型适用于预测黄河中游大型水库库区泥沙冲淤过程和趋势，本文选用三门峡水库潼关以下库区19601988年共28年的实测水沙系列进行验证计算。三门峡水库三个不同运用期库区及分河段冲淤总量的验证结果见表1，其中实测值为断面法测量结果。可以看出，各时段各河段计算的冲淤总量与实测值相差最大为0.9亿m3，一般为0.5亿m3，二者符合良好。累积冲淤过程以及每年汛期、非汛期的冲淤量也与实测结果符合较好，纵向冲淤分配在定性和定量方面与实际也较为接近。表1三门峡水库三个不同运用期库区冲淤总量验证(单位：亿m3)ComparisonofcalculatederosionanddepositionwiththemeasuredinSanmenxiareservoir河段潼关史家滩潼关太安太安史家滩时段实测计算实测计算实测计算1961.61965.1028.0827.817.918.4420.1819.361965.101974.6-3.48-3.52-0.50-0.21-2.98-3.311974.61988.10-0.83-0.54-0.87-0.480.03-0.064三门峡水库非汛期运用方式研究三门峡水库在非汛期一般不再承担或很少承担防凌任务，为了保证三门峡水库在5、6月份来水较小时能维持12台机组发电，需要配合小浪底水库进行春灌蓄水。为此，经考虑拟定如下运用原则：春灌蓄水位尽量降低，蓄水时间尽量缩短，使非汛期淤积的泥沙主要分布在大禹渡以下，以利于汛期近坝段溯源冲刷时排出库外，淤积末端限于以下，不影响潼关至河段的冲淤，不影响潼关高程的升降变化。总之，要改善三门峡水库的冲淤状况，尽量多发电。据此原则，初步设计了五种非汛期水库运用方案：(1)方案1，又称“四省会议”方案，即每年的11月次年的6月20日控制水位310m，6月21日6月30日水库水位逐渐降至305m；(2)方案2规定，每年的11月次年的6月20日控制水位312m，6月21日6月30日库水位逐渐降至305m；(3)方案3为，每年11月次年的6月20日控制水位315m，6月21日6月30日库水位逐渐降至305m；(4)方案4为，每年的11月次年的3月25日控制水位312m，3月264月20日最小出库流量700m3/s，允许最高蓄水位319m，4月21日6月20日库水位逐渐降至312m，6月21日6月30日水位逐渐降至305m；(5)方案5的蓄、降水过程和方案4相同，但水位不同，即每年的11月次年的3月25日控制水位315m，3月26日4月20日最小出库流量700m3/s，允许最高水位319m，4月21日6月20日库水位逐渐降至315m，6月21日6月30日水位逐渐降至305m。19901995年度潼关站非汛期平均来水量161.4亿m3，来沙量1.93亿t，与1993年度非汛期潼关水量155.0亿m3、沙量1.94亿t较接近，因而选择该年度非汛期潼关站的水沙条件(包括逐日流量、含沙量、悬沙级配等)作为上述五个方案计算的依据，起始地形采用1995年10月份统测的潼关以下库区大断面资料，计算的起止日期为1992年10月20日至1993年6月30日，统计该时段的实测水沙资料，潼关入库水量162.1亿m3，沙量2.0亿t。水电站发电量按下式估算E=AQHT(8)式中E为发电量，千瓦小时(度)；Q为过机组流量(m3/s)，非汛期过机组流量按出库流量考虑；A为出力系数，采用8.5，实际上不同水位该值有所不同；H为发电水头,m；T为发电时间，小时。各方案的计算结果见表2所示。由表可见，各方案潼关至NFDA2NFDA3河段计算的冲淤量基本一致，表明水库按照这些方案运用对潼NFDA2段的冲淤不产生影响。方案1至方案3的运用过程一样，只是控制水位不断抬高，从库区总冲淤量的计算结果看，方案3比方案2、方案2比方案1淤的多；五个方案中第5方案累积淤积量最大，含有春灌蓄水的方案4比方案3少淤0.121亿m3。从计算的非汛期发电量来看，基本上与累积冲淤量的大小相对应，即非汛期水库淤积量大，则发电量也大。综合起来看，方案4计算的累积冲淤量适中，发电量也偏多一些。表2三门峡水库非汛期运用各方案计算成果比较Calculatedresultsofdifferentmodesofoperationinnonfloodseasons方案特征水位(m)各河段冲淤量(亿m3)发电量桃汛前水位最高水位41-3636-3131-2222-141-1(亿度)P13103100.1658-0.1254-0.0120.58220.610711.7P23123120.164-0.11180.10790.60110.761212.45P33153150.1621-0.05210.36580.50190.977713.57P43123190.1653-0.07130.24940.51310.856513.07P53153190.165-0.01970.43940.43111.015913.895三门峡水库汛期各级流量排沙效果研究考虑到三门峡水库汛期的调沙减淤任务将主要由小浪底水库承担，这样就可以进一步研究水库在不同流量级当库水位降低时的排沙效果，因为从黄河下游河道减淤的需要来看，三门峡水库目前的泄洪排沙流量应大于2500m3/s，但近十几年的实测资料表明，这种流量出现的机率在不断减小，这就需要研究水库在较小流量时降水排沙的可能性。这次研究共选择了五级入库流量，分别为1000、1500、2000、2500及3000m3/s，持续时间均为20天，方案依流量从小到大的顺序编号，各方案均采用控制最低水位298m敞泄运用。各级入库流量的含沙量及悬沙级配均按近几年潼关站汛期的实测资料统计得出，方案1至方案5的入库含沙量分别为25、35、50、85及96kg/m3。起始库区地形采用1996年汛前统测的大断面资料。水库泄流规模按1990年以后打开9#、10#底孔以及扩装2台机组的情况考虑，采用天津设计院1993年提供的泄流曲线进行计算。模型输出结果主要有各级流量下水库持续20天敞泄库区冲刷量和北村流量为1000m3/s的水位变化过程。计算结果显示，在20天的计算时段内，库区冲刷主要集中在黄淤31断面以下，潼关以下库区累积冲刷量随着入库流量的增大而增大；当入库流量为1000m3/s时，即使库水位降到298m，库区的冲刷强度也不大，整个计算时段内北村水位降不到309m(北村站稳定水位)，相对而言，入库流量为2500m3/s时(方案4)北村水位降至稳定的时间最短，因为当入库流量达到3000m3/s时，按照目前的泄流规模敞泄，水库产生壅水，库水位超过298m，因而降低了对北村水位的冲刷作用。表3为三门峡水库19901992年汛期低水位排沙期实测的排沙情况统计。由表可见，流量在25003000m3/s和20002500m3/s这两级流量范围内水库冲刷强度最大,15002000m3/s这级流量的冲刷强度居中，流量小于1000m3/s以后库区冲刷明显减弱，这些与数模计算结果在定性上基本一致。表3三门峡水库19901992年汛期低水位排沙期排沙情况统计Sedimentdischargingatlowwaterstagesinfloodseasonsin1990-1992流量级(m3/s)天数tg(m3/s)smx(m3/s)sjt(m3/s)冲刷强度(108t/d)25003000628082740300.14-0.042200025001622552220300.51-0.038150020002917931744300.97-0.026100015003013581240299.71-0.013500100036780736300.61-0.0036结语1.小浪底水库建成投入运用后，三门峡水库非汛期运用的五个方案计算结果表明，各方案潼关至河段计算的冲淤量基本一致，表明水库按照这些方案运用对潼段的冲淤不产生影响。从计算的非汛期发电量来看，基本上与同期累积淤积量成正比。从减少库区淤积及增加电站出力两方面综合考虑，水库非汛期运用的第4方案较好。2.小浪底水库运用后，三门峡水库在汛期可以按照入库水沙条件和发电需要选择排沙流量，根据对汛期不同流量级的冲刷计算与各流量级排沙的实测结果来看，库区累积冲刷量随着入库流