水利工程论文-黄河下游引黄灌溉对策研究.doc

水利工程论文-黄河下游引黄灌溉对策研究摘要：在新的形势下，对黄河下游引黄灌溉中泥沙处理和利用应提高到一个新的认识高度，变泥沙“包袱”为“资源”。泥沙长距离输送，使绝大部分泥沙进入田间，就可实现引黄灌溉的良性循环。除了以往有效的措施外，本文在总结灌区渠道不同部位设置提水泵站设施作用与效果的同时，通过建立灌区渠道悬移质冲刷计算模型，对其减淤清淤的作用与效果进行了理论分析研究。长距离输送；提水灌溉关键词：长距离输送提水灌溉1泥沙是制约引黄灌溉发展的重要因素黄河下游河南、山东两省的引黄灌区，自1985年至1990年底，累计引水2333亿m3，引沙38.65亿t。其中，约有33.22%的泥沙沉积在沉沙池内，35.32%淤积在灌溉渠系，22.9%被输入田间，8.56%进入排水河道。这就是说，除了被输入田间的约23%的泥沙成为资源，其余近80%的泥沙需要进行处理。这种长期采取的泥沙处理方式(沉沙池集中沉沙，渠系两岸堆积清淤泥沙)，造成了引黄泥沙在平面上分布不合理的状况，给沿黄两岸的经济、社会和环境带来了严重的负面影响。引黄泥沙已成为制约引黄灌溉发展的主要因素。随着黄河整个流域经济的迅速发展，用水量不断增加，而黄河下游水资源却日趋贫乏。引黄灌区面临需水期无水可引的局面，不得不引用洪水期高含沙水流。但目前，大多数灌区渠首无可供高含沙水流沉沙的空间，以进行沉沙输水。面对这一形势，对泥沙的处理和利用应有新的认识，这一认识的基本点是变泥沙“包袱”为“资源”。从技术上讲就是要实现泥沙的长距离输送。实现泥沙长距离输送，关键在于提高渠道输送泥沙的能力。要减少输沙渠道的泥沙淤积，除了以往行之有效的措施外，可采用改输沙渠水流自流入沉沙池变为扬水入沉沙池的方式，这种措施必将增加渠道的水面比降，加大流速，增加渠道内水流输送泥沙的能力。如若在渠首设置提水泵站扬水入渠，同样可起到加大渠道水面比降的作用。利用提水泵站实现泥沙长距离输送的可行性已在部分灌区得到证实。例如，渠首设置提水灌溉设施的夹马口灌区和尊村灌区及渠尾(中)设置提水泵站设施的曹店灌区和麻湾灌区。这些灌区在不同部位设置了扬水设施来加大水面比降，提高水流挟沙能力，达到了少淤或不淤的目的。本文在参照典型灌区冲淤变化状况的基础上，建立了适用于灌溉渠道输沙计算的数学模型，就提水灌溉的效果与作用进行了定性的研究与分析。2渠道悬移质冲刷模型的建立2.1基本方程的建立及数值解法基于灌溉渠道水流泥沙运动规律与明渠水流泥沙运动规律一致，模型采用非耦合饱和输沙一维模型的计算方法，计算中假定渠道发生冲刷淤积过程，计算时段内河床变形对水流条件影响不大，也不考虑水体中含沙量随时间的变化，则该方程组为水流连续方程(1)水流运动方程(2)河床变形方程(3)水流挟沙力公式S*=S*(J,u,d)(4)式中Q为流量；A为过水面积；B为水面宽；Y为水深；u为流速；Y0为断面渠底平均高程；S*为水流挟沙力；G为断面总输沙率；为泥沙干密度；g为重力加速度；n为河床糙率；x,t为距离及时间。上述简化后的基本方程，即为建立模型理论依据。采用Pressman的四点偏心差格式，其中时间偏导采用前差，空间偏导采用加权中心差，可得函数和导数离散关系式式中为加权系数01，f表示任意函数；上标n为时间坐标序列号，下标j为断面序号。按上述差分格式推演水流连续方程(1)和水流运动方程(2)中各项，对计算点的未知量按泰勒级数展开，忽略增量fj的高次项，可使多元高次方程化为线性方程组。为简化计算方法，考虑渠道断面形态为梯形，经整理，对式(1)和式(2)分别有(5)(6)以上系数A1，A2，A3，A4，A5及B1，B2，B3，B4，B5均为时段初相邻断面变量值Zjn,Qjn,及进库水流条件的函数。上述两个方程中，共有四个未知数，对M个计算点共有2M个未知数，由差分方程可列出2M-2个类似于(5)、(6)的方程组。因此需要补充两个边界条件可以方程组闭合。使用交错网格追赶法求解差分方程。为将Q与Y分别布置在两个不同断面上，在断面间再增加一个节点(以K表示)，并假定增加节点断面流量与上一个断面相等。2.2初始条件及边界条件初始条件：Q(x)t=0=Q0(x),h(x)t=0=h0(x)上边界条件：Q1=Q(t)下边界条件：Ym=Y(t)2.3河床变形方程及水流挟沙能力河床变形方程式(3)，采用有限差分法进行离散，得到差分方程，如下(G2-G1)T=xBY0(7)式中G1、G2分别为进出口断面输沙率；为渠床泥沙干密度；Y0为计算河段平均河床淤积或冲刷厚度，正值为冲，负值为淤。对于悬移质变化引进的冲淤计算，取G=QS。水流挟沙能力采用模型相似类比法计算渠道断面挟沙力，由下式确定S*=S*0(8)其中S*=J7/4u3/2/d2式中J、u、d分别为水面比降、流速、悬移质中值粒径；S*为挟沙能力比尺；S*0为原始渠道平衡挟沙力。3提水灌溉作用与效果的分析研究3.1提水灌溉的作用以曹店灌区干渠为模拟计算对象，取代主要参数为：设计引水流量30m3/s,梯形断面，渠底宽度为14m，边坡系数为12，渠底比降为1/7000，渠道水深1.373m，断面平均流速1.342m/s,过水面积22.36m2。计算了下游泵站水位由7.373m下降至6.873m，即下降水位50cm，相应分流比约50%条件下渠道沿程含沙量的变化，见图1、图2。图1降水过程中含沙量的变化图2降水过程结束后含沙量的变化SedimentconcentrationvariationalongthecanalduringpumpingwaterattheendofitThesedimentconcentrationvariationafterpumpingended由图可看出，在降水过程中，沿程含沙量随降水幅度的增加而增长；降水过程结束后，沿程含沙量又趋于平缓。经计算分析，沿程水面比降和平均流速也存在着这样的规律。这是由于泵站提水，造成局部水位下降，泵站的上渠段水面比降加大，导致水流流速增加，该渠段发生冲刷，出口含沙量迅速增加；当降水结束时，断面最大含沙量达到150kg/m3；降水结束后，随着冲刷的不断发展，下游水深增大，水面比降和流速逐渐减小，加上冲刷过程中悬移质的粗化，水流挟沙能力减小，制约了冲刷的进一步发展。因此，建立提水泵站设施，可以对渠道起到一定的拉沙作用，但其拉沙效果则与降水幅度、降水历时及入流条件、渠道床面形态等有关。3.2提水灌溉效果的单因素分析利用建立的模型，进行了包括降水模式(包括降水方式、幅度、历时)、入流量、渠道断面形态、糙率等诸因素对提水灌溉效果的计算，见表1、表2、表3、表4。表1降水模式对冲刷效果的影响Theeffectoferosionindifferentdurationofpumping引水模式降水历时(h)降低出口水位值(cm)延伸距离(km)平衡历时(h)冲刷深度(cm)冲刷量(万m3)分段降水2.0305.53.615.700.770集中降水0.2307.55.011.800.760分段降水2.0508.05.237.462.072集中降水0.2509.05.035.402.297表1说明，长时段的分段降水和短时段的集中降水对冲刷效果的影响差别不大，而降水幅度对冲刷效果的影响则很大。表2渠道入流量对冲刷效果的影响Theeffectoferosionfordifferentflowdischarge河宽(m)流量(m3/s)挟沙力(kg/m3)延伸距离(km)平衡历时(h)冲刷深度(cm)冲刷量(万m3)分段降水50cm入流量对冲刷的影响142023.9368.838.771.87114303085.237.462.072144035.159.54.437.162.442145039.6911438.982.775集中降水50cm入流量对冲刷的影响142023.9386.834.092.0281430309534.952.297144035.15104.435.42.422145039.6910437.762.585显然，渠道入流量对降水冲刷效果的影响比较大。流量越大冲刷越深、拉沙距离越远、所需平衡历时越少、冲刷量也越多；即越易冲刷。表3不同宽深比B/h条件下的冲刷效果Theeffectoferosionfordifferentvalueof/hofcanal渠底宽(m)水深(m)水面宽(m)宽深比/h过水面积(m2)挟沙力S*(kg/m3)延伸距离(km)平衡历时(h)冲刷深度(cm)冲刷量(万m3)分段降水101.59615.322.4520.2136.8585.638.482.762121.47416.912.7921.3130.0085.237.462.07250cm不同宽深比条件下的冲刷效果141.37318.583.1422.3626.4884.836.781.735181.21222.043.8724.2822.8774.835.941.379集中降水50cm不同宽深比条件下的冲刷效果101.59615.322.4520.2136.859.55.434.13.115121.47416.912.7921.3130.009.5534.92.297141.37318.583.1422.3626.4894.835.41.836181.21222.043.8724.2822.878.54.837.11.454从表3中可看出，宽深比越小冲刷效果越好，且有随宽深比的减小而递增的趋势。这说明了当下游条件改变时，渠道越宽浅冲刷效果越差。相反，窄深断面形态的渠道相对而言更易于冲刷。表4降低水位50cm糙率的影响Theeffectoferosionfordifferentroughnesscoeff.降水模式糙率延伸距离(km)平衡历时(h)平衡冲深(cm)平衡冲刷(万m3)分段降水0.0118.56.631.91.818分段降水0.0206.53.219.880.999分段降水0.0305.02.612.040.572集中降水0.0119.07.831.221.842集中降水0.0209.02.417.961.025集中降水0.0308.01.410.240.579从表4可看出，糙率越小冲刷效果越显著。反之，糙率越大冲刷效果越差，而且回淤也越严重。各表说明水位下降幅度是影响冲