水利工程论文-黄河口拦门沙疏浚对尾闾河道泄洪排沙影响的研究.doc

水利工程论文-黄河口拦门沙疏浚对尾闾河道泄洪排沙影响的研究摘要：本文研究和探讨了黄河口拦门沙疏浚后，尾闾河道水位、流速和水面比降的变化以及疏浚部位的回淤过程。试验结果表明，河口拦门沙疏浚后，在一定时期和一定范围内，可以使尾闾河道发生溯源冲刷。在施放不同典型年的水沙条件下，尾闾河道可在疏浚点以上2135km范围内发生溯源冲刷，同流量水位最大下降0.240.47m。同时，疏浚河槽始终处于回淤状态，经过一个汛期回淤量可达到70%80%。关键词：黄河口拦门沙疏浚泄洪排沙黄河口地区河床主要处于一个堆积的环境，当河流来沙不能全部输送至深海，则在河口地区发生沉积。从纵剖面上看，一般都存在着突出于上下游河段河底连线之上的成型堆积体，其中淤积部位处于河口段与口外滨海段的交接地区，亦即口门附近，称之为拦门沙。拦门沙形成之后，侵蚀基面抬高，对河口泄水排沙不利，导致水位壅高，泥沙沉积，产生溯源淤积，对尾闾河道具有负面的反馈影响。近年来，随着经济社会的持续发展，黄河两岸工农业用水量的持续增长，进入河口的水沙条件发生变异，尾闾河道主槽萎缩严重，防洪压力增大。因此，为减轻黄河下游河道及河口地区的淤积，一些专家提出采取疏浚、拖淤等措施治理河口拦门沙，尽可能保持口门畅通，以利泄洪排沙入海。然而，河口拦门沙演变剧烈，影响因素极为复杂，疏浚效果难于预测。目前黄河口尚未建立实体模型，同时在黄河口拦门沙地区还缺乏挖沙疏浚的实践，尚有很多关键技术问题和理论问题亟待解决。为此，本文采用概化物理模型1对河口拦门沙疏浚效果进行试验研究，为黄河口综合治理提供参考。1试验的基本情况1.1模型比尺模型设计为水平比尺L=1000，垂直比尺H=50，模型下边界设在海底高程达-14m的海域，上边界设在西河口以上1.5km处，模拟尾闾河道及河口长66km(如图1所示)，采用电木粉作为模型沙。由于模型的变率较大，必然造成一定的相似性损失，为尽可能保证模型与天然相似，模型试验按照模型相似率理论进行设计，并通过水位及冲淤量等模型率定试验对有关比尺进行适当的调整，确保主要物理量相似，使试验结果满足研究内容的要求。图1模型平面布置天然河道的水流一般处于紊流状态，保证水流相似的条件为重力相似和阻力相似，则流速比尺V=1/2H=7.07，糙率n=1/6H(H/L)1/2=0.43。在天然情况下，原型试验段的糙率系数一般在0.0090.014之间，则要求模型糙率一般在0.02090.0325之间。对于动床模型，电木粉的糙率系数在0.02左右，小于设计要求的糙率。但由于概化模型的变率较大，河道的形态阻力和边壁阻力较大，模型的综合糙率基本能够满足设计要求，这一点在验证试验中得到了印证。悬移质运动相似的条件为悬移、起动和冲淤相似。起动相似条件要求起动流速比尺与流速比尺相等，即Vc=V。沉降相似比尺=V(H/L)0.75=0.75，粒径比尺d=(V/s-)1/2=0.45。为保证模型中河床冲淤与原型相似，含沙量比尺经过多次验证试验调整后，模型试验最终采用s=0.9，即s=s*=0.9。冲淤时间比尺：t2=L/VS=341。1.2验证试验自1996年黄河口汊河流路行河以来，进入黄河口的水沙连年偏枯，只有1996年发生了较为典型的洪水，同时又恰逢汊河流路改道行水当年，尾闾河道发生了较大的冲淤变化。因此，选择1996年汛前实测河道和河口地形作为验证试验的初始地形，1996年610月实测水沙资料作为验证试验的水沙条件。表1给出了各级典型流量模型试验水位与天然实测水位，可以看出，误差在0.12m以内，表明模型基本达到了阻力相似；表2给出了分段冲淤量的验证结果，基本上能够满足河床冲淤相似性的要求。1.3方案试验条件1.3.1地形条件方案试验的初始地形，河道部分根据1996年汛后河道大断面测量资料及2000年实测的11万的河道地形图确定。河口及海域部分根据1996年实测的12.5万河口及海域地形确定。在此基础上，根据清水沟流路顶点高程最大的拦门沙实测地形资料，塑造了拦门沙I+1(见图2)。拦门沙I+1顶点约在S7处，顶坡段长约2km，顶点高程0m，自顶坡起始点向上游为一倒坡，坡度约为-1/2200，倒坡段长约4km，最低高程为-1.8m，坡顶与坡脚高差为1.8m。1.3.2水沙条件根据试验内容，河道上游来水来沙条件采用恒定水沙过程和非恒定水沙过程。恒定流量分为大、中、小三级，根据该河段的水沙特点确定为1000m3/s、2000m3/s和3000m3/s，相应的含沙量分别为15kg/m3、30kg/m3和40kg/m3，来流持续时间分别为26d、13d和8d，下边界采用固定潮位(分别为平均中潮位0.5m和平均高潮位1.0m)。非恒定水沙过程采用黄委会设计院提供的小浪底水库正常运用期丰、平、枯不同典型年(利津站)资料(见表3)，下边界采用典型潮位过程。表1不同流量级水位验证(单位：m)流量/(m3/s)100020003000西河口原型6.757.508.18模型6.727.608.07误差-0.030.10-0.11丁字路口原型3.384.244.62模型3.504.184.56误差0.12-0.06-0.06表2冲淤量验证结果(单位：万m3)河段CS7清3清3清7清7汊3CST汊3模型-975-665-361-2001原型-842-581-490-1913表3典型年水沙特征值典型年丰水丰沙(1988)中水平沙(1995)枯水枯沙(1991)径流量/亿m3295.4148.790.1输沙量/亿t9.764.531.241.3.3疏浚概况在拦门沙2.0km长的顶坡段顺水流方向开挖一宽500m、深约2.0m的河槽，同时在拦门沙4km长的倒坡段做相应疏通开挖，使得在拦门沙顶坡段和倒坡段形成一个顺畅的河槽，开挖量相当于原型约为380万m3。图2河口拦门沙I+1地形图3疏浚前后的流速分布2试验成果分析2.1典型流量的疏浚试验成果分析2.1.1疏浚前后的流速变化在试验中观测到，河口拦门沙疏浚前，进入河口的径流自寻低洼捷径入海，在径流和潮流的作用下形成口门众多的支汊，每股水流虽有一定的河槽，但比较散乱，出汊摆动比较频繁。河口拦门沙疏浚后，在口门地区人为塑造了一个比较明显的河槽，水流归顺。在试验中，量测了不同径流和潮位的组合下，河口拦门沙疏浚前后的流速变化。由图3可见，拦门沙疏浚后与疏浚前相比，由于水流集中入海，疏浚部位的流速都有所增大，疏浚后行水之初流速可以增加10%20%，但随着疏浚部位的回淤，流速逐渐减小。2.1.2疏浚前后尾闾河道水位的变化疏浚河口拦门沙的主要目的是降低侵蚀基准面，使尾闾河段主槽河床降低，增加主河槽的过水断面，降低水位，提高尾闾河道的泄洪排沙能力。因此，水位的分析包括以下两方面的内容：一是疏浚前后同流量级下水位的升降，反映了同一断面槽蓄能力的变化；二是疏浚前后水面比降的调整，反映了尾闾河段泄洪排沙能力的变化。图4给出了疏浚前后沿程水位的变化情况。由于疏浚部位处于河海的交汇处，受河流动力和海洋动力的双重作用。可以看出，河口拦门沙的疏浚，降低了侵蚀基准面，在疏浚部位以上的尾闾河道产生了溯源冲刷；由于受平均潮位的控制，疏浚区同流量水位落差向海域方向快速减小至零，即达到平均潮位。因此，疏浚前后同流量的水位落差沿程分布呈V字形。由图可见，拦门沙疏浚后，尾闾河道溯源冲刷的大小和范围与黄河的水沙条件和潮位密切相关。在造床流量范围内，随着流量的增大，溯源冲刷范围增大；相同流量下，潮位增高，溯源冲刷的大小和范围有所减小。2.1.3疏浚前后尾闾河道水面比降的变化疏浚河口拦门沙的目的是降低侵蚀基准面，在尾闾河道的下端形成一定的水位落差，并由此产生自下而上的溯源冲刷，下面就这一物理图景进行分析。对于恒定均匀流，一维水流和泥沙方程可以表述如下。水流连续方程：Q=BhU(1)水流运动方程：(2)泥沙连续方程：(3)式中：z为河床高程；C为谢才系数；为干容重。输沙率QS的计算公式较多，如采用QS与河床比降J的一次方成正比的公式24，代入式(3)，并联解上述三式可得河床高程变化的偏微分方程(4)式中：K为综合系数。河口拦门沙疏浚前后，在尾闾河道末端形成了水位落差，其边界条件可以考虑为,0x，0tz(x,0)=0z(0,t)=h(t)(5)假设h(t)=-hH(t)，其中H(t)=0，t0；H(t)=1,t0也就是h(0)=-h和解得(6)因此，发生溯源冲刷后的水面线为一下凹型曲线。由图5可见，河口拦门沙疏浚后，在上游的影响段内，发生了溯源冲刷，并呈下凹曲线型式；在疏浚河段及以下，由于受平均潮位的控制，水面比较平缓。图4疏浚前后沿程水位差值的变化由表4可以看出，河口拦门沙疏浚后，在上游溯源冲刷影响的河段内，水面比降变陡。在平均潮位为0.5m的