钱塘江河口围垦回淤过程预测探讨

1、钱塘江河口围垦回淤过程预测讨论摘要：在河口进展围垦，由于边界条件的改变和纳潮量的变化，一般都会引起河床的冲淤变化，对淤积过程的预测是一个非常重要且困难的问题。本文根据定床潮流计算的结果，利用河床变形方程求得围垦后的初始淤积速率，利用最小能耗原理求得极限淤积平衡状态，而由灰色理论求得整个淤积过程，经动床物理模型的印证和实测资料的检验以及工程施行后的开展过程的验证，均说明本文提出的预测方法是可行的，具有一定的可靠性和预测精度。关键词：河床变形回淤过程预测灰色模型1前言钱塘江河口潮强流急，宽浅多变，根据假设干年来河口治理理论，对钱塘江河口主要采用趁潮淤围涂、固定江道、以围代坝等治江与围涂相结合的措施

2、来进展河口整治。围垦之后，由于边界条件的改变，导致了水沙运行环境的改变，涨落潮流特性随之变化，净输沙格局也相应变化。一般地，围垦后，槽蓄量减小，涨、落潮动力减小，潮流所挟带的沙量也随之减小，与之相适应，河床过水断面将减小，产生回淤现象，回淤的幅度和时间过程对围垦工程而言是一个非常重要且困难的问题。由于泥沙问题的复杂性和不成熟性，采用泥沙数学模型进展长历时的计算受到诸多影响因素的制约，计算精度也就受到影响；采用物理模型代价又很大。因此，如何估算泥沙回淤过程是一个很有价值和值得研究的问题。本文从半理论半经历的途径出发，在定床潮流计算的根底上，利用河床变形方程、最小能耗原理、灰色模型对此问题进展了讨

3、论，与物理模型结论和实测资料分析结论相一致，并且与工程施行后的回淤过程相吻合。2回淤机理概述河流系统是一个开放的非线性系统，它与周围环境存在着物质和能量的交换，系统内部的自组织、自我调整作用是其演变开展的内因，外部环境的改变那么是其变化的外因，外因要通过内因此起作用。围垦工程的施行，改变了来水来沙条件和水沙运行环境，必将引起河床的自我调整。决定河床形态变化趋势和程度的自我调节机制是受河流体系能量分配和耗散原理的制约，能量分配的调整使体系内能量分配具有最大或然率。这要求体系在调整过程中力求使熵到达最大，因此使各个可能出现状态的或然率相等。根据最小能耗原理，当流域的来水来沙条件发生变化时，河流将主

4、要通过调整能量耗散率VJ来影响这种变化，而水流能量耗散率的沿程调整将由河宽B和水深H的沿程调整均等地承当，从而有VJ=VV2n2/H4/3=q3n2/H13/3=in(1)由式(1)可知，为适应边界条件的变化，河流系统中反响最为敏感、最易调整的因子是水深H，其次是单宽流量q，再次是糙率n，而变化相对缓慢、滞后的是in，它与施加于系统的约束有关。在围垦后，边界条件发生改变，在河流系统的自我调整作用下，各水力要素均会相应地发生变化以适应变化了的环境条件，由(1)式易得H2/H1=(q2/q1)9/13(2)=(n2/n1)6/13(in1/in2)3/13(3)式中1代表工程前，2代表工程后，H为

5、水深，q为单宽流量。该构造形式已在钱塘江河口和其它航道整治工程中得到了广泛的应用，本文也将据此来推求河床调整平衡后的可能冲淤厚度。钱塘江河口以悬移质造床为主，在围垦前后，其糙率变化不大，即n1n2;并且部分的围垦尚缺乏以引起河口系统的能耗最小值的大幅度相对变化，因此也可假设(in1/in2)3/131，为此近似地戎1。根据(2)式，围垦后纳潮量减小，即q2减小，那么H2应减少，河床将发生淤积，其回淤关系将满足式(2)，利用定床水流计算所得的水力要素和(2)式可求得围垦后可能的平衡淤积厚度。另一方面，对于潮汐河口，由于潮波变形，涨、落潮不对称，对输水输沙的作用是不一样的。在某时刻的单宽输沙率可写

6、成(4)式中,分别为潮周期的垂线平均流速、含沙量和水深，V，S，H分别为某一时刻平均流速、含沙量和水深与垂线平均流速、含沙量和水深的偏向。以潮周期为计算单元时，有V=S=H=0(5)故潮周期平均单宽输沙率为(T1)(T2)(T3)(T4)(T5)(6)式中T1项中，恒为正，正负两可，因此，T1的正负(方向)仅取决于涨、落潮流的强度比照。T1表征了优势流对净输沙的奉献。T2项反映了stkes漂流效应对净输沙的奉献，在河口地区该项为负值。T3的符号依VS的正负而定，后者既取决于流速过程线的含沙量过程线和相位关系，还与涨、落潮掀沙效应的差异有关。涨潮含沙量大于落潮含沙量时，T3为负，反之为正。故T3

7、反映了涨、落潮流挟沙强度对净输沙的奉献。T4为含沙量与水深变化的相关项，因研究区段处在动力平衡带，AKV-的绝对值通常较小，故T4的绝对值也较小，T5项依赖于V，S，H的相关性，其绝对值与前三项相比小得多。由方程(4)可知，潮流对净输沙的奉献取决于涨、落潮流特性的相对变化。围垦之后，由于流速、潮量的减小，涨、落潮动力减小，潮流所挟带的沙量也随之减小，涨、落潮优势流和挟沙强度的相对变化决定了围垦之后河口的净输沙格局的变化和演变开展趋势。工程前后，河道的冲淤主要取决于优势流与挟沙强度的相对变化，为此可用输水输沙的相对变化来判断可能发生的冲淤情况。比方，对涨潮冲刷槽而言，当涨潮输沙量大于落潮输沙量时

8、，将发生冲刷，可用优势沙表示G=(VHSt)f/(VHSt)e(Vf3tf)/(Ve3te)SV2/H(7)式中下标f表示涨潮，下标e表示落潮，S为含沙量。同样，定义优势流为涨潮潮量与落潮潮量之比，即优势流为涨、落潮相对输水量之比F=(Qf/Qe)(8)3回淤预测公式一般形式根据河床变形方程1可建立如下冲淤形式(9)式中k为沉降机率(在钱塘江河口可取12)；为泥沙沉降速度(4510-4/s);s为含沙量；s*为挟沙力；为泥沙干容重。在围垦之前(天然情况下)，可认为河床处于相对冲淤平衡状态，即，有s1=s1*，在潮汐河口，有s*=(v2/gh)(10)围垦后初期，假定来沙不变，即s2=s1，由于

9、围垦引起水流条件的改变，必然引起挟沙才能的相应变化，造成河床的冲淤变化，其冲淤强度可由(9)式的差分而求得，即p=k(s2-s2*)/(11)代入挟沙关系(10)可得p=ks1/(1-(v2/v1)2(h1/h2)(12)上式中，当p0时，可认为得到冲淤极限平衡状态，即h2H2，从而可求得平衡水深H2H2=(q2/q1)aH1(13)式中=2/3,q1=v1h1,q2=v2h2，水力要素均由定床潮流计算求得。上式在钱塘江河口中已到达了广泛的应用，它与前文根据最小能耗原理所得的平衡关系式(2)相接近，但该式隐含了工程前后挟沙力不变的假设，从而也说明了最小能耗原理更具有一般性，可用于一般的河床演变

10、分析。当围垦工程施行后，由于河床变形相对滞后，可假设围垦后初期含沙量不变，从而可利用定床计算所得的水力要素，由式(12)可求得围垦后初期的淤积速率p。对于潮汐河口，正如前述，由于涨、落潮对泥沙的输移作用是不一样的，为了利用上述由冲积河流所得的公式(12)，可以将公式中对应的水力要素取全潮平均进展计算。当取半潮平均时，可以分别得到半潮淤积速率p1(涨潮)和p2(落潮)，那么淤积速率取二者相比拟大者，p=ax(p1,p2)，或根据当地的冲淤特性来选取，对涨潮占优者取p1，对落潮占优取p2，对涨、落持平者，取p=(p1+p2)/2作为初期淤积速率。4回淤过程的灰色预测上文已求得围垦后初期的可能淤积速

11、率和平衡淤积厚度。由于河床的自我调整作用，随着淤积的开展，水流要素也会受到影响而发生变化，从而进一步地影响河床的回淤速率。从原那么上讲，可由(12)式进展迭代求解整个回淤过程，但由于在迭代过程中，水沙要素的变化不易确定，为进一步求解带来了一定的困难。在此，采用灰色预测方法来预估回淤过程。由于泥沙系统是一个本征性的灰色系统，可以采用G(1，1)2模型来进展预测，利用上文已求得的初始淤积速率和平衡淤积厚度，可以求得回淤过程和淤积速率的变化过程。对G(1.1)模型有(14)其中,为待定系数，z(1)为累积淤积厚度。对灰色模型(14)而言，有解析解(15)对上式求导，可得淤积速率满足(16)为方便起见

12、，以初始地形为基准，即取z0=0。显然，当淤积得到平衡时，有；在淤积初期有，根据上述条件和前文已求得的平衡淤积厚度和初始淤积速率，可以求得待定系数，从而得到回淤过程。转贴于论文联盟.ll.5应用实例*论文联盟.LL.编辑。根据上述思想，在钱塘江河口尖山河湾南股槽整治研究中，建立了平面二维潮流模型进展数值模拟，利用该定床模型所求得的水流参数和实测含沙量资料，对围垦后的河床变形和演变趋势进展了分析预测。*宋立松等.钱塘江尖山河湾南股槽整治数学模型整治研究之三.浙江省河口海岸研究所(内部报告)，1996.12。5.1定床潮流模型针对钱塘江尖山河湾特点，建立了平面二维潮流模型，其控制方程为连续方程：(

13、17)运动方程：(18)(19)式中z为潮位，即水面到达一基准面的间隔 ；U，V分别为x,y方向上的垂线平均流速分量；h为水深；g为重力加速度；f为柯氏力参数(f=2esin,为纬度，e为地球自转速度)；z为谢才系数；x,y分别为水流在x,y方向的涡动扩散系数；x,y分别为x,y方向的风应力分量。上述方程组的初始条件z(x,y)|t=0=z0(x,y)U(x,y)|t=0=U0(x,y)V(x,y)|t=0=V0(x,y)图1钱塘江尖山河湾南股槽治理方案示意图SkethfsuthbranhntrlsheefJianshanbendinQiantangestuary边界条件水边界:z(x,y,t

14、)=z*(x,y,t)“*表示值，陆边界:法线方向流速为零；在南股槽整治研究中，选用了无构造网格的差分方法对控制方程进展离散，即对式(18)，(19)应用不规那么网格差别离散，对离散方程进展了数值求解，求得围垦前后的有关水力要素(U1,V1,h1,U2,V2,h2),从而可得到：,q1=v1h1,q2=v2h2,有了以上参数，即可应用本文方法进展围垦后回淤过程的预测。5.2回淤预测从优势流(F)与优势沙(G)的相比照拟来看，二者所反映的趋势是根本一致的，淤积根本上以东进闸为中心(图1)，处于淤积环境中，个别地方反而略有冲刷，下游较远处(临海浦以下)根本上可维持平衡。利用定床水流计算的水力结果(

15、v2,v1,h2,h1,q1=v1h1,q2=v2h2)和实测资料分析所得的泥沙参数(s1,)，由式(2)和式(9)求得平衡淤积厚度和初始淤积速率如表1。表1河床变形预测Freastingresultsfriverbeddefratin方案位置东进闸临海浦近期平衡淤积厚度()44493方案初始淤积速率(/d)4.100.77中期平衡淤积厚度()839164方案初始淤积速率(/d)5.091.11根据表1，由(15)和(16)求得东近闸附近的回淤过程满足如下方程对近期方案(围垦6万亩)有(20)(21)对中期方案有(围垦8.6万亩：即近期6万亩及中期2.6万亩)(22)(23)由上述方程可知，近

16、期方案施行后的淤积平衡厚度约为4.44，初始淤积速率为4.1/d，得到平衡所需时间约为300500天左右；中期方案施行后的淤积平衡厚度为8.39，初始淤积速率约为5.09/d，得到平衡的时间约为500700天左右(如图2所示)，可见其淤积速度是较快的。图2近期、中期方案施行后东进闸附近的回淤过程PreditedaggradatinpressnearDngjingbarrage6结果评价与讨论根据70年代的实测资料分析*，尖山河湾的时段冲淤量(划分冬半年、夏半年)约0.30.61083。从尖山河湾屡次汊道发生淤积的实测资料说明(1969年7.5洪水北沙串通，南股槽淤积和1973年，1985年尖山

17、一带主槽由南靠北，南岸汊道的淤积)，当其中的一股汊道出现淤积环境时，其淤积速度是很快的，一般约0.51年。当床面高程为-10时，淤积率23/d。南股槽内水深较大，因此估计围垦后其淤积速率前期会略大于上述实测淤积率，因此，按上述形式所求得的淤积强度略大是可以承受的，并且由上述方法所得到的冲淤部位，淤积平衡厚度、回淤时间等都与动床物理模型*的结论相一致，也与实测资料分析结果相符合，19961997年近期方案施行后出现的淤积相貌及其开展过程均与本文预测吻合良好(图2)，从而也说明上述预测的河床变形有一定的参考价值和可信度。而且利用所得到的灰色模型，还可以从数学的角度3，4对河床演变进展敏感性、稳定性

18、等各种动态行为的分析，为进一步提醒河口泥沙运动规律提供了一种途径。*余祈文等，钱塘江尖山河湾南股槽的形成及河床演变特性研究，浙江省河口海岸研究所(内部报告)，1992.12。*熊绍隆等，钱塘江尖山河湾南股槽整治动床模型整治研究之二，浙江省河口海岸研究所(内部报告)，1996.12。当然，由于河口泥沙问题的复杂性和特殊性，要准确地预测河床长历时的变形是很困难的，目前尚无有效手段。采用泥沙数学模型进展长历时的计算受到诸多影响因素的制约，计算可靠性及精度也就受到影响；采用物理模型费时又费钱，代价很大；定床潮流模型只能反映水力变化，不能反映河床变形；最小能耗原理或河相关系只能反映河床的平衡状态，难以提醒其开展变化过程；灰色模型可以描绘过程，但其动态行为对初值、终值有一定的依赖性，也不能完全确定河床演