会员注册 | 登录 | 微信快捷登录 支付宝快捷登录 QQ登录 微博登录 | 帮助中心 人人文库renrendoc.com美如初恋!
站内搜索 百度文库

热门搜索: 直缝焊接机 矿井提升机 循环球式转向器图纸 机器人手爪发展史 管道机器人dwg 动平衡试验台设计

   首页 人人文库网 > 资源分类 > DOC文档下载

水利工程论文-钱塘江河口围垦回淤过程预测探讨.doc

  • 资源星级:
  • 资源大小:59.47KB   全文页数:7页
  • 资源格式: DOC        下载权限:注册会员/VIP会员
您还没有登陆,请先登录。登陆后即可下载此文档。
  合作网站登录: 微信快捷登录 支付宝快捷登录   QQ登录   微博登录
友情提示
2:本站资源不支持迅雷下载,请使用浏览器直接下载(不支持QQ浏览器)
3:本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰   

水利工程论文-钱塘江河口围垦回淤过程预测探讨.doc

水利工程论文钱塘江河口围垦回淤过程预测探讨摘要在河口进行围垦,由于边界条件的改变和纳潮量的变化,一般都会引起河床的冲淤变化,对淤积过程的预测是一个非常重要且困难的问题。本文根据定床潮流计算的结果,利用河床变形方程求得围垦后的初始淤积速率,利用最小能耗原理求得极限淤积平衡状态,而由灰色理论求得整个淤积过程,经动床物理模型的印证和实测资料的检验以及工程实施后的发展过程的验证,均表明本文提出的预测方法是可行的,具有一定的可靠性和预测精度。关键词河床变形回淤过程预测灰色模型1前言钱塘江河口潮强流急,宽浅多变,根据若干年来河口治理实践,对钱塘江河口主要采用趁潮淤围涂、固定江道、以围代坝等治江与围涂相结合的措施来进行河口整治。围垦之后,由于边界条件的改变,导致了水沙运行环境的改变,涨落潮流特性随之变化,净输沙格局也相应变化。一般地,围垦后,槽蓄量减小,涨、落潮动力减小,潮流所挟带的沙量也随之减小,与之相适应,河床过水断面将减小,产生回淤现象,回淤的幅度和时间过程对围垦工程而言是一个非常重要且困难的问题。由于泥沙问题的复杂性和不成熟性,采用泥沙数学模型进行长历时的计算受到诸多影响因素的制约,计算精度也就受到影响采用物理模型代价又很大。因此,如何估算泥沙回淤过程是一个很有价值和值得研究的问题。本文从半理论半经验的途径出发,在定床潮流计算的基础上,利用河床变形方程、最小能耗原理、灰色模型对此问题进行了探讨,与物理模型结论和实测资料分析结论相一致,并且与工程实施后的回淤过程相吻合。2回淤机理概述河流系统是一个开放的非线性系统,它与周围环境存在着物质和能量的交换,系统内部的自组织、自我调整作用是其演变发展的内因,外部环境的改变则是其变化的外因,外因要通过内因而起作用。围垦工程的实施,改变了来水来沙条件和水沙运行环境,必将引起河床的自我调整。决定河床形态变化趋势和程度的自我调节机制是受河流体系能量分配和耗散原理的制约,能量分配的调整使体系内能量分配具有最大或然率。这要求体系在调整过程中力求使熵达到最大,因而使各个可能出现状态的或然率相等。根据最小能耗原理,当流域的来水来沙条件发生变化时,河流将主要通过调整能量耗散率VJ来影响这种变化,而水流能量耗散率的沿程调整将由河宽B和水深H的沿程调整均等地承担,从而有VJVV2n2/H4/3q3n2/H13/3min1由式1可知,为适应边界条件的变化,河流系统中反应最为敏感、最易调整的因子是水深H,其次是单宽流量q,再次是糙率n,而变化相对缓慢、滞后的是min,它与施加于系统的约束有关。在围垦后,边界条件发生改变,在河流系统的自我调整作用下,各水力要素均会相应地发生变化以适应变化了的环境条件,由1式易得H2/H1mq2/q19/132mn2/n16/13min1/min23/133式中1代表工程前,2代表工程后,H为水深,q为单宽流量。该结构形式已在钱塘江河口和其它航道整治工程中得到了广泛的应用,本文也将据此来推求河床调整平衡后的可能冲淤厚度。钱塘江河口以悬移质造床为主,在围垦前后,其糙率变化不大,即n1≈n2并且局部的围垦尚不足以引起河口系统的能耗最小值的大幅度相对变化,因此也可假设min1/min23/13≈1,为此近似地取m≈1。根据2式,围垦后纳潮量减小,即q2减小,则H2应减少,河床将发生淤积,其回淤关系将满足式2,利用定床水流计算所得的水力要素和2式可求得围垦后可能的平衡淤积厚度。另一方面,对于潮汐河口,由于潮波变形,涨、落潮不对称,对输水输沙的作用是不一样的。在某时刻的单宽输沙率可写成4式中,,分别为潮周期的垂线平均流速、含沙量和水深,V′,S′,H′分别为某一时刻平均流速、含沙量和水深与垂线平均流速、含沙量和水深的偏差。以潮周期为计算单元时,有∑V′∑S′∑H′05故潮周期平均单宽输沙率为T1T2T3T4T56式中T1项中,,恒为正,正负两可,因此,T1的正负方向仅取决于涨、落潮流的强度对比。T1表征了优势流对净输沙的贡献。T2项反映了stokes漂流效应对净输沙的贡献,在河口地区该项为负值。T3的符号依V′S′的正负而定,后者既取决于流速过程线的含沙量过程线和相位关系,还与涨、落潮掀沙效应的差异有关。涨潮含沙量大于落潮含沙量时,T3为负,反之为正。故T3反映了涨、落潮流挟沙强度对净输沙的贡献。T4为含沙量与水深变化的相关项,因研究区段处在动力平衡带,AKV的绝对值通常较小,故T4的绝对值也较小,T5项依赖于V,S,H的相关性,其绝对值与前三项相比小得多。由方程4可知,潮流对净输沙的贡献取决于涨、落潮流特性的相对变化。围垦之后,由于流速、潮量的减小,涨、落潮动力减小,潮流所挟带的沙量也随之减小,涨、落潮优势流和挟沙强度的相对变化决定了围垦之后河口的净输沙格局的变化和演变发展趋势。工程前后,河道的冲淤主要取决于优势流与挟沙强度的相对变化,为此可用输水输沙的相对变化来判断可能发生的冲淤情况。比如,对涨潮冲刷槽而言,当涨潮输沙量大于落潮输沙量时,将发生冲刷,可用优势沙表示GVHStf/VHSte∝Vf3tf/Ve3teS∝V2/H7式中下标f表示涨潮,下标e表示落潮,S为含沙量。同样,定义优势流为涨潮潮量与落潮潮量之比,即优势流为涨、落潮相对输水量之比FQf/Qe8G2G1,F2F1时,河床将发生冲刷,当G2G1,F2F1时,河床将发生淤积。由优势流和优势沙的概念,可以利用定床计算的水力要素来判断围垦后河床发生冲、淤的部位及可能发生的演变趋势。根据上述概念所得到的冲、淤平面分布与工程实施后出现的分布情况相一致。3回淤预测公式一般形式根据河床变形方程1可建立如下冲淤模式9式中k为沉降机率在钱塘江河口可取1~2ω为泥沙沉降速度4~5104m/ss为含沙量s为挟沙力γ′为泥沙干容重。在围垦之前天然情况下,可认为河床处于相对冲淤平衡状态,即,有s1s1,在潮汐河口,有sαv2/gh10围垦后初期,假定来沙不变,即s2s1,由于围垦引起水流条件的改变,必然引起挟沙能力的相应变化,造成河床的冲淤变化,其冲淤强度可由9式的差分而求得,即pkws2s2/γ′11代入挟沙关系10可得pkws1/γ′1v2/v12h1/h212上式中,当p→0时,可认为得到冲淤极限平衡状态,即h2→H2,从而可求得平衡水深H2H2q2/q1aH113式中α2/3,q1v1h1,q2v2h2,水力要素均由定床潮流计算求得。上式在钱塘江河口中已达到了广泛的应用,它与前文根据最小能耗原理所得的平衡关系式2相接近,但该式隐含了工程前后挟沙力不变的假设,从而也说明了最小能耗原理更具有一般性,可用于一般的河床演变分析。当围垦工程实施后,由于河床变形相对滞后,可假设围垦后初期含沙量不变,从而可利用定床计算所得的水力要素,由式12可求得围垦后初期的淤积速率p。对于潮汐河口,正如前述,由于涨、落潮对泥沙的输移作用是不一样的,为了利用上述由冲积河流所得的公式12,可以将公式中对应的水力要素取全潮平均进行计算。当取半潮平均时,可以分别得到半潮淤积速率p1涨潮和p2落潮,则淤积速率取二者相比较大者,pmaxp1,p2,或根据当地的冲淤特性来选取,对涨潮占优者取p1,对落潮占优取p2,对涨、落持平者,取pp1p2/2作为初期淤积速率。4回淤过程的灰色预测上文已求得围垦后初期的可能淤积速率和平衡淤积厚度。由于河床的自我调整作用,随着淤积的发展,水流要素也会受到影响而发生变化,从而进一步地影响河床的回淤速率。从原则上讲,可由12式进行迭代求解整个回淤过程,但由于在迭代过程中,水沙要素的变化不易确定,为进一步求解带来了一定的困难。在此,采用灰色预测方法来预估回淤过程。由于泥沙系统是一个本征性的灰色系统,可以采用GM1,12模型来进行预测,利用上文已求得的初始淤积速率和平衡淤积厚度,可以求得回淤过程和淤积速率的变化过程。对GM1.1模型有14其中,为待定系数,z1为累积淤积厚度。对灰色模型14而言,有解析解15对上式求导,可得淤积速率满足16为方便起见,以初始地形为基准,即取z00。显然,当淤积得到平衡时,有在淤积初期有,根据上述条件和前文已求得的平衡淤积厚度和初始淤积速率,可以求得待定系数,从而得到回淤过程。5应用实例根据上述思想,在钱塘江河口尖山河湾南股槽整治研究中,建立了平面二维潮流模型进行数值模拟,利用该定床模型所求得的水流参数和实测含沙量资料,对围垦后的河床变形和演变趋势进行了分析预测。宋立松等.钱塘江尖山河湾南股槽整治数学模型整治研究之三.浙江省河口海岸研究所内部报告,1996.12。5.1定床潮流模型针对钱塘江尖山河湾特点,建立了平面二维潮流模型,其控制方程为连续方程17运动方程1819式中z为潮位,即水面到达一基准面的距离U,V分别为x,y方向上的垂线平均流速分量h为水深g为重力加速度f为柯氏力参数f2ωesinφ,φ为纬度,ωe为地球自转速度Cz为谢才系数εx,εy分别为水流在x,y方向的涡动扩散系数τx,τy分别为x,y方向的风应力分量。上述方程组的初始条件zx,y|t0z0x,yUx,y|t0U0x,yVx,y|t0V0x,y图1钱塘江尖山河湾南股槽治理方案示意图SketchofsouthbranchcontrolschemeofJianshanbendinQiantangestuary边界条件水边界zx,y,tzx,y,t表示已知值,陆边界法线方向流速为零在南股槽整治研究中,选用了无结构网格的差分方法对控制方程进行离散,即对式18,19应用不规则网格差分离散,对离散方程进行了数值求解,求得围垦前后的有关水力要素U1,V1,h1,U2,V2,h2,从而可得到,q1v1h1,,q2v2h2,有了以上参数,即可应用本文方法进行围垦后回淤过程的预测。5.2回淤预测从优势流F与优势沙G的相对比较来看,二者所反映的趋势是基本一致的,淤积基本上以东进闸为中心图1,处于淤积环境中,个别地方反而略有冲刷,下游较远处临海浦以下基本上可维持平衡。利用定床水流计算的水力结果v2,v1,h2,h1,q1v1h1,q2v2h2和实测资料分析所得的泥沙参数s1,ω,由式2和式9求得平衡淤积厚度和初始淤积速率如表1。表1河床变形预测Forecastingresultsofriverbeddeformation方案位置东进闸临海浦近期平衡淤积厚度cm44493方案初始淤积速率cm/d4.100.77中期平衡淤积厚度cm839164方案初始淤积速率cm/d5.091.11

注意事项

本文(水利工程论文-钱塘江河口围垦回淤过程预测探讨.doc)为本站会员(奋斗不息)主动上传,人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知人人文库网([email protected]),我们立即给予删除!

温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。

copyright@ 2015-2017 人人文库网网站版权所有
苏ICP备12009002号-5