水利工程论文-倒虹吸工程河段洪水与河床变形的数值模拟.doc水利工程论文-倒虹吸工程河段洪水与河床变形的数值模拟.doc

收藏 分享

资源预览需要最新版本的Flash Player支持。
您尚未安装或版本过低,建议您

水利工程论文倒虹吸工程河段洪水与河床变形的数值模拟摘要本文利用平面二维水沙数学模型,采用水边界全区自动跟踪方法并考虑了河床局部冲刷与河道演变分析成果,对渠河交叉工程河段进行了洪水与河床变形的数值模拟。模型经过历史特征洪水验证,模拟计算了特征洪水的流场、交叉工程处的壅水及冲刷,根据流场冲刷状况提出了工程优化修改方案。关键词渠河交叉二维水沙数值模拟壅水床面冲刷1研究问题南水北调中线总干渠沿线与许多河流交叉,其中在河北穿越七里河的交叉建筑物型式为渠穿河倒虹吸工程。倒虹吸设计长度初选为700M,设计洪水标准为百年一遇,设计洪峰流量2410M3/S。渠穿河倒虹吸工程的修建对该处河道水流及河床演变产生的影响,倒虹吸工程的位置、尺寸、埋置深度是干渠设计所关心的重大问题。利用二维水沙数学模型,可以较好地模拟反映渠、河交叉工程附近水流与河床变形状况,主要研究在该河修建渠穿河倒虹吸后,交叉工程附近河段流场流速、壅水及河床冲淤变形程度。通过河流模拟,分析、评价交叉工程设计方案并提出工程修改建议。2交叉段河道特性在交叉工程附近七里河分为两汊,属宽浅型沙质河床河道北槽较大,主流靠左岸。中泓处的河床质表层为粗沙,D50在10MM左右,滩地及两岸为中沙及壤土。该河为冲积性河流,河床形态、演变规律与洪水造床作用和常年水沙特性密切相关。据河道查勘与河床地形图分析1963年洪水前期引起河床较强的冲刷,洪水后河道又普遍回淤;河段滩槽高差一般为2~3M,河槽宽850~1200M。交叉工程附近河床组成沿流程分选明显,由粗至细变化较大。1980年以来河道受人类活动影响严重,由于河道内多处挖取粗沙,原主槽回淤的泥沙被大量挖走,有的挖深达3~4M。多年小水作用及人类挖沙已经使现河道形成人为不连续窄深槽,窄槽宽度一般为100~300M。这造成在现状河道小洪水易于归槽且侧侵蚀较为明显,见图1。3平面二维水沙数值模拟31数学模型的基本理论对于宽浅型河流,水深平均的二维水沙控制方程可较好地反映河流中挟沙水流运动特征。本模型的水流基本方程由三维时均雷诺方程沿水深积分得到,并以混长紊流模型求解紊动切应力模型采用悬移质泥沙扩散方程与河床变形方程求解河床冲淤变化,由床面冲淤临界切应力判断床面泥沙冲淤状态及床面稳定条件。该数学模型已在一些复杂工程中得到成功应用3,能较准确地模拟、预测一般冲积性河流上,河流工程附近的水沙运动与河床变形。图1不同时期河床横断面的变化CHANGESOFCROSSSECTIONINDIFFERENTTIMEPERIODS311控制方程水流连续方程1水流运动方程23悬移质输运扩散方程4河床变形方程5水流挟沙力方程SKU3/GRΩM6对于散粒沙河床,床面稳定控制的辅助方程可表示成局部区域稳定控制条件ΤC≥Τ0或ΤC/Τ0≥17式中ΤCΤC/ΓSΓDFUCD/V8Τ0Τ0/ΓSΓDΓRJ/ΓSΓD9式中Ξ水位;H水深;U,VX,Y向水深平均流速;U合速度;U摩阻流速,R水力半径;E,Ν水流涡粘系数、运动粘滞系数;Β对流项修正系数,S含沙量,Z河床高程,J水力坡度;D,Α泥沙扩散系数、恢复饱和系数;Ω泥沙沉速;CCHEZY系数;S挟沙力,K、M挟沙力系数、指数,ΓS、Γ′泥沙容重与干容重;Τ0床面无量纲水流切应力,ΤC无量纲谢尔兹临界切应力;Τ0床面水流切应力,ΤC谢尔兹临界切应力;D床面分层粒径。312数值计算格式1离散网格及变量分布模型采用非均匀网格,可以在研究量变化梯度较大的重要局部区域设置细密网格,在研究量变化较平缓的非主要区域设置较稀疏的网格。在离散网格上,标量Ξ,H,S被安排在单元中央,矢量速度U,V安排在单元的四周,交错网格上物理变量的位置相互错开。分别用四个一维数组XU,YV,XH,YH来确定水位、各流速分量的坐标3。2离散格式本数学模型采用较成熟的“交替方向隐式差分逐行求解”方法,其特点是将时间步长TIMESTEP分成前后两个半步,在前半个时间步长取某一个方向为隐式;为保持对称,在后半个时间步长改变隐式方向。每一个时间步长,都这样交替的改变隐式方向计算。在前后两个时间半步对控制方程进行离散,为了物理概念上的清晰和格式的稳定有效,引进控制体同时在进行对流项离散时引起迎风格式。在前半个时间步长,将连续方程与Y向动量方程联立,对U,Ξ进行隐式求解;在得到水流条件后隐式求解关于含沙浓度的传移输运方程。在后半个时间步长,将连续方程与Y向动量方程联立,对Ν,Ξ进行隐式求解;也在获得水流条件以后隐式求解泥沙传移输运方程。313初始条件及边界条件1初始条件对于给定的计算区域,在时间T0时,令Ξ|T0Ξ0X,Y;U|T0U0X,Y;Ν|T0Ν0X,Y2开边界条件ΞX,Y,TΞOPBX,Y,T或QQOPBT以及SΞX,Y,TSOPBX,Y,T。其中ΞOPB,QOPB以及SOPB分别为开边界上已知的水位、流量以及含沙量,一般由实测水文资料确定。对于具有宽滩深槽复杂地形的河道,模型进口开边界上宜给定洪水水位过程或流量过程,便于处理全断面流量分配。上游开边界所需相应的悬移质含沙量及过程如无法实测资料,则借用附近河流相应洪水的含沙量及过程。下游开边界处的含沙量可按第二类边界条件确定。陆地边界即河流岸边,满足固壁非穿越及无滑移条件,边界上的法向、切向流速为零。3动边界处理在宽浅河流中的边滩和江心洲随着水位的变化,其水边线也不断改变,形成所谓动边界。模型相应采用水边界全区自动跟踪的处理方法,设置一个跟踪指标数组IWET。先将最大可能的淹没区域包纳在计算域区,在计算过程中根据计算单元内的水深来判断该单元是淹没或是露出,即IWET应该赋0还是赋10表示露出,为陆地;1表示淹没,为应计算的水域,凡是陆地单元均不纳入计算范围。这种动边界智能跟踪的处理方法特别适合边滩和江心洲淹没与出露频繁的水域,使流场边界条件模拟更为真实。32七里河数学模型的建立321数学模型范围及网格剖分根据数值模拟的研究任务和交叉工程附近河段的河道形态,河床组成及洪痕调查情况,充分考虑上、下游边界的水流条件,确定模型上边界取在距交叉工程上游40KM处;模型下边界取在距交叉工程下游45KM处,沿河道模型总长约85KM。侧边界在考虑洪水的最大可能淹没范围及地形条件,基本沿陡坎近岸高地选取。模型最大宽度为25KM,计算区域总面积为2125KM2。为了反映交叉工程河段河道地形,适应不同区域流场计算精度要求,交叉工程附近的重点区域,采用加密网格为5050M;对于一般非重点计算区域,网格尺寸为10050M。模型全部节点为4500个。322数学模型调试与验证1调试依据与要求模型调试主要依据是1963年8月和1996年8月该河发生的两场重要洪水的调查洪痕以下简称“638”洪水、“968”洪水;水文分析所得计算河段进出口附近断面的水位流量关系以及河床冲淤变形有关的查勘资料与河床探测资料。河道地形采用1966年1/10000地形图并参考1994、1996两年汛后河道的纵横断面测量资料。考虑“968”洪水河道侧侵蚀严重的特点,对河道地形进行了适当修正。“638”洪水属特大洪水,“968”洪水属一般大洪水,分别选用这两场洪水验证模型流场,可使模型能正确模拟不同洪水条件,不同阻力特征的流场。2调试结果二维水沙数学模型反演“638”、968”洪水的计算成果经后处理,绘制成流场流速矢量图。流场矢量图所反映的流场的主流走向、流态基本与调查情况一致。这里仅给出“968”洪水流场中主要洪痕点水位调查值与计算值的对比,见图2。河床糙率上段取0030~0035,下段取0025~0032。通过对洪水流场主流走向、流态、流势、流速分布等方面模拟与洪痕点水位校验对比,数学模型反演的“638”和“968”洪水与实际洪水调查情况基本一致,正确反映了两类洪水的洪水河势及行洪特点。数学模型在河道地形处理,糙率选择,冲淤强度控制参数方面的选择是合理的。33数学模型计算成果与分析本次计算根据交叉工程设计方案,对设计洪水进行了数值模拟,然后根据流场状况,进一步提出修改工程方案并进行数值模拟和分析对比。331自然状态计算成果在未修工程的现状条件下,对设计洪水进行了自然状态河段的数值模拟。由洪峰期计算河段平面流速分布图3可以看出,在设计洪水时,交叉工程断面附近水流收缩集中,河槽冲刷严重。受江心洲影响,河道分成两汊又渐次合一北汊较大,南岸高地下游有部分回流区。图2“968”洪水期的洪痕点对比COMPARISONBETWEENCALCULATEDANDMEASUREDWATERLEVLESINCHARACTERISTICPOINTSIN968FLOOD图3交叉河段自然状态下洪水流场设计洪水NATURALFLOODFLOWINGFIELDINCROSSREACHDESIGNEDFLOOD332设计方案计算成果交叉工程设计倒虹吸段长度为700M,按设计洪水进行了设计工程方案的计算。为了反映洪水演进过程中有关特征值的变化,特别是为确定最高洪水位,最大冲刷深度等有关成果,根据设计洪水过程的特点,设定洪水计算过程中相应的计算成果输出时刻,一般选取23~29个时刻,包含洪水演进的各特征时期。这里仅给出洪峰期流场流速矢量分布图4以及交叉工程断面处流速及冲淤分布图5。从计算结果看,总干渠交叉工程修建后,在倒虹吸工程河段的主流分为两股,河道北侧深槽略大一些。两股深槽处单宽流量和流速均较大,最大流速达38M/S。冲刷也主要集中在两股深槽,最大冲刷在北侧深槽,最大冲刷深度约为194M。从冲刷发展过程看,随洪水发展两股深槽的冲刷也不断发展变化,南槽持续冲刷而北槽则在洪峰过后又略有回淤。图4倒虹吸工程影响下的洪水流场设计洪水FLOODFLOWFIELDINTERFERINGWITHSIPHONDESIGNEDFLOOD333修改方案计算成果根据自然状态和设计工程状态下流场流态、壅水、冲刷深度及部位的计算成果分析,得知南岸滩地水流流速不高,近期调查南槽仅略有增大趋势,在此适当约束水流不会产生很大影响。因此可以考虑进行修改、优化工程布置,这里提出了分别在南、北岸倒虹吸工程进口、出处将倒虹吸段缩窄50M、和将倒虹吸工程全部安排在北槽的三种修改方案,相应倒虹吸工程段长度均为650M。对不同修
编号:201312122137044886    类型:共享资源    大小:50.85KB    格式:DOC    上传时间:2013-12-12
  
2
关 键 词:
专业文献 建筑水利 精品文档 水利工程
  人人文库网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
关于本文
本文标题:水利工程论文-倒虹吸工程河段洪水与河床变形的数值模拟.doc
链接地址:http://www.renrendoc.com/p-204886.html

当前资源信息

4.0
 
(2人评价)
浏览:22次
wenku上传于2013-12-12

官方联系方式

客服手机:17625900360   
2:不支持迅雷下载,请使用浏览器下载   
3:不支持QQ浏览器下载,请用其他浏览器   
4:下载后的文档和图纸-无水印   
5:文档经过压缩,下载后原文更清晰   

精品推荐

相关阅读

人人文库
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们

网站客服QQ:2846424093    人人文库上传用户QQ群:460291265   

[email protected] 2016-2018  renrendoc.com 网站版权所有   南天在线技术支持

经营许可证编号:苏ICP备12009002号-5