水利工程论文-黄河三门峡库区泥沙模型的设计.doc水利工程论文-黄河三门峡库区泥沙模型的设计.doc

收藏 分享

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

水利工程论文黄河三门峡库区泥沙模型的设计摘要本文针对多沙水库特点,首次采用黄河泥沙模型相似律和异重流运动相似条件的最新研究成果,完成了三门峡水库泥沙动床模型的模型沙选择和比尺设计。利用水库自然滞洪淤积及降水冲刷资料进行的验证试验结果表明,模型较好地复演了原型的水沙运动规律及河床变形。关键词三门峡水库模型设计异重流1引言开展小浪底水库运用方式研究,需要通过物理模型这一重要手段进行试验,为选择水库最优运用方式提供科学依据。由于小浪底水库正在兴建中,其模型不具备验证试验的条件,因此,我们选用位于小浪底水库上游,水沙条件及河床边界条件较为相近,且具有丰富观测资料的三门峡水库进行验证。亦即通过三门峡库区模型的验证试验,确定小浪底水库模型的主要比尺。本文运用黄河泥沙模型相似条件及异重流运动相似条件,开展了三门峡泥沙模型的设计。2模型相似条件开展黄河水库模型设计,其相似条件除包括水流重力相似、水流阻力相似、泥沙悬移相似、水流挟沙相似、河床冲淤变形相似、泥沙起动及扬动相似外,还要着重考虑异重流运动相似条件。对于泥沙沉速比尺ΛΩ、含沙量比尺ΛS等基本比尺,按照黄河泥沙模型相似律1开展设计。而在考虑异重流运动相似方面,则需要根据笔者的最新研究成果2进行计算,即应满足异重流发生或潜入相似条件ΛSEΓΛK11/ΓSMΓ/ΓSMSPΛK1ΛΓSΓ/ΛΓS1异重流挟沙相似条件ΛSEΛSE2水流连续相似条件ΛTEΛL/ΛV3式1~3中足标“M”、“P”、“E”分别代表模型、原型及异重流有关值;ΓS、Γ为泥沙及水的容重;ΛS为含沙量S的比尺;ΛL为水平比尺;ΛV为流速比尺;ΛT为时间比尺。式1中ΛK1为考虑浑水容重沿垂线分布不均匀性而引入的修正系数的比尺,修正系数K1的定义式为4在运用式4时,尚需引入异重流含沙量分布公式。由于紊动扩散作用及重力作用仍是决定异重流挟沙运动的一对主要矛盾,其浓度沿水深的分布及挟沙能力规律与一般挟沙水流应当类似。因此,作为模型设计可引用张红武的含沙量沿垂线分布公式1计算异重流含沙量沿垂线分布。为保证模型与原型水流流态相似,还许满足浑水有效雷诺数REM8000含沙量较小时,REM为一般雷诺数的4倍及模型水深HM15CM的限制条件。3几何比尺及模型沙选择31几何比尺确定从满足试验精度要求出发,根据原型河床条件、HM15CM的要求及对模型几何变率问题的前期研究结果,确定水平比尺ΛL300,垂直比尺ΛH45,几何变率DTΛL/ΛH667。分别采用张红武提出的变态模型相对保证率、窦国仁提出的模型变率限制条件及张瑞瑾等提出的模型变态指标等表达式进行计算,表明本模型采用DT667,在各家公式所限制的变率范围之内,几何变态的影响也有限,可以满足工程实际需要。32模型沙选择本次验证试验需要模拟原型库区冲淤调整幅度较大的情况,即要保证淤积相似,又要保证冲刷相似,因此对模型沙的基本特性有更高的要求。经验表明,有些种类的模型沙在潮湿的环境中固结严重,将使模型河床冲淤相似性明显偏差特别是影响冲刷过程的相似性。清华大学水利水电工程系曾于1990年开展了D50≤0038MM的电木粉起动流速试验,结果为H10CM时,初始条件下VC108CM/S;水下沉积两天后,VC增加到12CM/S;在水下沉积两个月后,VC21CM/S;而脱水固结二周后,即使流速增至28CM/S,电木粉也不能起动。由我们开展的郑州热电厂粉煤灰ΓS2058KN/M3,D500035MM及山西煤屑ΓS147KN/M3,D50005MM两种模型沙的起动流速试验结果图1看出,在相近水深条件下,山西煤屑的起动流速随着沉积时间增加有大幅度的增加。例如在水深同为4CM条件下,水下固结96小时后,起动流速从初始的595CM/S达到84CM/S,脱水固结96小时后可达131CM/S。而郑州热电厂粉煤灰的起动流速虽然随固结时间增加有所增大,但增大的幅度明显较小。图1不同模型沙起动流速试验结果图RESULTSOFTHRESHOLDVELOCITYFROMDIFFERENTMODELSANDS1、1′初始;2、2′水下固结48小时;3、3′水下固结96小时;4煤屑脱水固结48小时;4′煤灰脱水固结120小时;5煤屑脱水固结96小时;大量研究表明1,郑州热电厂粉煤灰的物理化学性能较为稳定,同时还具备造价低、宜选配加工等优点。因此选用郑州热电厂粉煤灰作为本动床模型的模型沙。该模型沙土力学特性试验成果见表1。表1郑州热电厂粉煤灰作为悬沙用土力学特性试验成果表TESTINGRESULTSOFSOILMECHANICSFORFLYASHESUSEDASSUSPENDEDSEDIMENTFROMTHEZHENGZHOUHEATANDPOWERPLANT容重ΓSKN/M3干容重Γ0T/M3内摩擦角水下休止角凝聚力CKG/CM22058066312295~3050064模型比尺计算41流速及糙率比尺由水流重力相似条件求得流速比尺ΛV671,由此求得流量比尺ΛQΛVΛHL90585;取水力半径比尺约等于水深比尺,即ΛR≈ΛH,由阻力相似条件求得糙率比尺ΛN073。对于黄河水库库区的模型,在回水变动区河床糙率模拟的正确与否会直接影响到回水长度及淤积分布。根据三门峡水库北村断面实测资料,其糙率值一般为0013~002,由此求得模型糙率应为NM00178~00274。为分析模型糙率是否满足该设计值,利用文献1中的公式及预备试验结果对模型糙率进行分析,即5式中Κ为卡门常数,为简便计取Κ04;若取原型水深为5M,则HM5M/450111M;Χ为校正参数,对于床面较为粗糙的模型小河,取Χ1;HS为模型的沙波高度,根据预备试验HS002~0028M。由式5求得模型糙率值NM0017~0019,与设计值接近,初步说明所选模型沙在模型上段可以满足河床阻力相似条件。至于库区近坝段,其水面线主要受水库运用的影响,而河床糙率的影响相对不大。42悬沙沉速及粒径比尺由三门峡库区北村站、茅津站水文泥沙实测资料,可求得悬浮指标Ω/ΚU015,对于Ω/ΚU≤015的细沙,其悬移相似条件可表示为1(6)将三门峡库区测验资料及有关比尺代入式6,得出ΛΩ的变化幅度为120~144,平均约为134。由于原型及模型沙都很细,可采用滞流区公式计算沉速,由此可得到悬沙粒径比尺关系式ΛDΛΩΛΝ/ΛΓSΓ1/27式中ΛΝ为水流运动粘滞系数比尺,该比尺与原型及模型水流温度及含沙量大小等因素有关,若原型及模型两者水温的差异较大,可使ΛΝ有很大的变化幅度,进而使ΛD有较大的取值范围。显然,在模型设计时给ΛD一定值是不合适的,合理的方法是在试验过程中根据原型与模型温差等条件适当调整ΛD。43模型床沙粒径黄科院的研究表明1,不同种类的模型沙,由于其容重、颗粒形状等方面存在较大差异,尚不能直接由现有的泥沙起动流速公式计算模型沙的起动流速,而且这些公式用于天然河流特别是黄河,其计算结果也会偏小不少。正因如此,对于黄河沙质河床的模型设计,不能直接采用泥沙起动流速公式推求模型床沙的粒径比尺,而不得不分别确定原型泥沙的起动流速和模型沙的起动流速,然后判断两者的比值即ΛVC是否满足起动和扬动相似条件。张红武在开展黄河河道模型设计时,根据罗国芳等收集的资料,点绘与三门峡库区河床组成相近的泥沙不冲流速与床沙质含沙量的关系曲线,并视该曲线含沙量等于零的流速为起动流速,由此曲线得出H1~2M时,VC≈090M/S。在水库淤积或冲刷过程中,床沙粒径变幅较大。据实测资料统计,床沙中径变化幅度一般为0018~008MM。由土力学知识,泥沙中径为006~008MM,可划分为中壤土或轻壤土;中径为0025~006MM,可划归为重壤土或中壤土一类。由文献3查得当水深为1M时,两者起动流速VC分别约为07M/S及09M/S。在水深为22CM时的起动流速为010~013M/S。通过模型沙起动流速试验,发现中值粒径D500018~0035MM的郑州热电厂粉煤灰作为模型沙,相应的起动流速比尺与流速比尺相等。附带指出,由我们初步点绘的郑州热电厂粉煤灰在水深为5CM时,起动流速VC与中径D50的点群关系来看图2,在D500018~0035MM的范围内,即使横坐标变化了近2倍,VC的变化并没有超出目前水槽起动试验的观测误差。由此说明模型沙粗度即使与理论值有一些偏差,也不致于对泥沙起动相似有大的影响。当水深增加时,原型沙起动流速将有所增加,由文献3可知,一般情况下,不冲流速VBVC1H1/4,式中VC1为H1M时的不冲流速。根据我们及文献1给出的郑州热电厂粉煤灰起动试验资料,可得知在原型水深为1~20M的范围内,上述初选的模型沙可以满足起动相似条件。例如当原型水深为12M时,由此求起动流速为130M/S。由模型沙的起动流速试验得出VCM175CM/S。则起动流速比尺ΛVC743,与上述ΛΝ接近。根据窦国仁及张红武水槽试验结果1,与原型情况接近的天然沙的扬动流速一般为起动流速的154~175倍。若取原型扬动流速VF165VC,可求得原型水深为3~6M的床沙扬动流速VFP165092~110152~182M/S。参阅文献1资料,模型相应的床沙扬动流速VFM为023~027M/S,则相应求出ΛVF661~674,与ΛΝ接近,表明模型所选床沙可以近似满足扬动相似条件。在多沙河流上修建水库后,水库上段及回水段必然出现再造床过程,并通过多因素的综合调整,力求实现新的均衡形态。由三门峡库区河床横断面变化过程可以看出,尽管三门峡建库后的河床平面形态受两岸的制约,但河床的调整变化仍具有冲积河流的特性。因此,对于多沙水库模型,应尽量兼顾河型相似条件。采用三门峡库区北村及茅津站相当于造床流量下的有关实测资料进行分析计算,求得的北村河段河床综合稳定指标ZW值为图2郑州热电厂粉煤灰起动流速VC与中径D50关系H5CMRLATIONSHIPBETWEENTHRESHOLDVELOCITYVCOFFLYASHESANDITSMEDIANDIAMETERD5077~102,茅津河段ZW值为105~131,表明本河段处于游荡及弯曲两种河型之间的过渡型据张红武的研究,ZW≤5为游荡型;ZW>15为弯曲型;介于两者之间为过渡型。将上述所选模型沙中径及其它相应因子代入,所得模型ZWM值与原型值相近,因此本模型可以满足河型相似条件。44含沙量比尺及时间比尺含沙量比尺可通过计算水流挟沙力比尺确定。采用文献1提出的同时适用于原型沙及轻质沙的水流挟沙力公式,即S25Ξ00022SVV3/ΚΓSΓM/ΓMGHΩLNH/6D500628式中Κ为卡门常数;ΓM为浑水容重;ΩS为泥沙在浑水中的沉速;V为流速;H为水深;D50为床沙中径;SV为以体积百分比表示的含沙量;Ξ为容重影响系数,可表示为Ξ17/ΓSΓ2259对于本次选用的模型沙ΓS约为21T/M3,则Ξ25。对原型沙,ΓS27T/M3,则Ξ1。将北村、茅津、小浪底水文站测验资料及相应的比尺值代入,可分别得到原型、模型水流挟沙力SP及SM。大量数据表明,两者之比SP/SM变化幅度在152~194之间,一般为160~180,取其平均值,ΛS约为170。另一方面,为在模型中较好地复演异重流的运动,含沙量比尺应兼顾式1,将三门峡水库异重流观测资料代入,并把由此得到的ΛΚ1表达式与式1联解,即可求出异重流含沙量比尺ΛSE145~192。在模型试验中,为保证异重流沿程淤积分布及异重流排沙特性与原型相似,还应满足异重流挟沙相似条件式2。与上述挟沙机理同理,可将异重流观测资料代入张红武水流挟沙力公式,计算原型及模型的异重流挟沙力,进而确定ΛS16~19,与式1得出的结果基本一致,并且与上述水流挟沙相似条件确定的ΛS也较为接近,因此,选用ΛSΛSE17可同时满足明渠水流及异重流挟沙相似条件,又能满足异重流发生相似条件。我们在白沙水库模型中发现2,不遵循水流连续相似条件,将导致模型水库蓄水过程严重失真,根本无法开展异重流运动和水库泄水排沙的模拟观测。对于本模型水流运动时间比尺ΛT1ΛL/ΛV447。而河床冲淤变形时间比尺ΛT2ΛΓ0/ΛSΛT1,还与泥沙干容重比尺ΛR0及含沙量比尺ΛS有关。根据郑州热电厂粉煤灰进行的沉积过程试验,测得模型沙初期干容量为066T/M3D500016~0017MM。至于原型淤积物干容重,通过三门峡库区实测资料分析认为,水库下段初始淤积物干容重一般为10~122T/M3,可取115T/M3。由原型及模型沙干容重求得ΛΓ0174,进而可以根据河床冲淤变形相似条件计算出ΛT2458,可见,与水流运动时间比尺接近,对于所要开展的非恒定流库区动床模型试验,可以避免常遇到的两个时间比尺相差甚远所带来的时间变态问题,也不致于对水库蓄水、排沙及异重流运动的模拟带来不利的影响。45模型高含沙洪水适应性预估及比尺汇总在小浪底水库的调水调沙运用中,可能出现高含沙洪水输沙状况。因此水库模型设计应考虑对高含沙洪水模拟的适应性。上述模型设计在确定含沙量比尺的过程中,已经考虑了高含沙洪水泥沙及水力因子的变化。为进一步预估模型中有关比尺在高含沙洪水期是否适应,下面以ΛS170为条件开展初步分析。
编号:201312122238086824    类型:共享资源    大小:60.26KB    格式:DOC    上传时间:2013-12-12
  
2
关 键 词:
专业文献 建筑水利 精品文档 水利工程
  人人文库网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
关于本文
本文标题:水利工程论文-黄河三门峡库区泥沙模型的设计.doc
链接地址:http://www.renrendoc.com/p-206824.html

当前资源信息

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

官方联系方式

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

精品推荐

相关阅读

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

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

copyright@ 2016-2018  renrendoc.com 网站版权所有   南天在线技术支持

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