粘性泥沙冲刷源项公式的改进.pdf

《粘性泥沙冲刷源项公式的改进.pdf》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、泥沙研究2013年10月JournalofSedimentResearch第5期粘性泥沙冲刷源项公式的改进雷文韬，夏军强，谈广鸣(1．长江科学院河流研究所，湖北武汉430010；2．武汉大学水资源与水电工程科学国家重点实验室，湖北武汉430072)摘要：为研究粘性泥沙的冲刷规律，本文引入床沙相对干密度以反映冲刷源项与床沙干密度的关系，提出分别适用于弱固结和强固结河床的冲刷源项公式。采用室内粘性泥沙的冲刷试验资料确定公式中各系数，然后采用天然沙及细颗粒石英沙的冲刷试验资料验证强固结河床的冲刷源项公式。计算得到的冲刷源项与实

2、测值符合较好，结果表明这两个公式能够用来计算水沙数学模型中粘性泥沙的输移过程。关键词：粘性泥沙；冲刷源项；相对干密度；弱固结；强固结中图分类号：TV149．1文献标识码：A文章编号：0468—155X(2013)05-0035-06在水沙数学模型中，泥沙输移方程中源项的研究对于含沙量过程及河床变形的数值计算具有重要意义。其通常表示为泥沙冲刷源项与淤积源项的差值，其中淤积源项可根据泥沙沉速及含沙量进行计算，而冲刷源项用于描述水流冲刷作用与床沙抗冲能力之间的关系，其计算方法较为复杂。粘性泥沙的物理和化学性质决定了其在水流作用

3、下难以冲刷，冲刷性质较非粘性泥沙存在很大不同。现有数学模型多采用恢复饱和系数形式的公式计算均匀沙或非均匀非粘性沙的冲刷源项¨，，这些公式能够较好地描述非粘性沙的冲刷过程。已有关于粘性沙数学模型的研究表明对于粘性和非粘性沙的冲刷源项分别采用不同处理方法能够更好地描述非均匀沙的冲刷过程。Mehta等认为河床的冲刷类型可分为有限冲刷和无限冲刷两类。据此Sanford等提出同时适用于这两类冲刷的粘性沙源项公式，但该公式很难应用于数学模型。用于描述河床无限冲刷的源项公式在粘性沙数学模型中应用广泛，但公式系数的确定具有较大的经验性，

4、不同研究区域的取值往往相差较大。为了描述粘性沙冲刷源项与床沙干密度的关系，本文基于Mehta等对于粘性沙冲刷源项计算公式的研究，在其公式的冲刷系数中引入床沙相对干密度改进得到分别适用于弱固结和强固结河床的冲刷源项公式。然后采用室内粘性沙的冲刷试验资料率定公式中的各系数；最后采用已有的基于天然沙及细颗粒石英沙的冲刷试验资料对改进后的强固结河床冲刷源项公式进行验证。1已有冲刷源项公式1．1冲刷源项的影响因素从物理意义及量纲来看，粘性沙冲刷源项等同于冲刷率，表示在水流作用下单位时间单位面积床面冲起的泥沙质量，其主要影响因素可归

5、结为反映水流对泥沙冲刷作用及泥沙抗冲能力的各种因素。对于河床有限冲刷，现有的现场和室内冲刷试验研究’均观测到河床冲刷过程逐渐趋于停滞的现象。Sanford等指出影响粘性沙冲刷源项的床沙干密度及临界冲刷切应力均与水深及冲刷历时相关，当经过一定历时临界冲刷切应力与床面水流切应力相等时，河床冲刷随之停止。Tan等。。指出影收稿日期：2013-01-21基金项目：国家自然科学基金资助项目(51079103)；高等学校博士学科点专项科研基金资助项目(20120141110011)；中央级公益性科研院所基本科研业务费项目(CKSF2

6、013006／HL)．作者简介：雷文韬(1986一)，男，湖北武汉人，博士研究生，主要从事河口泥沙运动等方面研究。E-mail：leiwt@whu．edu．el35响粘性沙冲刷源项的床沙干密度与其固结历时相关，床沙干密度在固结过程中逐渐增大并接近其稳定值，河床随之变得越来越难以冲刷。现有关于河床无限冲刷的研究¨。。”同样表明粘性沙冲刷源项与水流切应力及床沙干密度相关。其中Owen最早研究了粘性沙冲刷源项与水流切应力的关系，并提出采用剩余水流切应力描述冲刷源项的线性公式。Robe．s等¨提出关于水流切应力及床沙干密度的冲刷

7、源项公式，并认为当床沙中值粒径d。0．222mm时冲刷源项与干密度无关。Lick等¨在对天然沙冲刷规律的研究中指出粘性沙冲刷源项与床沙干密度显著相关。1．2冲刷源项公式有限冲刷通常发生于河床表面弱固结层，而大多数情况下河床表现为无限冲刷占优，因此粘性沙数学模型更多的是采用适用于无限冲刷的源项公式。二维悬移质泥沙不平衡输移方程可写为(s)+(^Js)+芸(s)：(sOS)+(s_OS)+F(1)aaoyaOX,oyoy式中F为泥沙输移源项；为时间；h为水深；、分别为、Y方向上的流速；S为含沙量；为泥沙紊动扩散系数。由泥沙不

8、平衡输移引起的河床变形可由下式表示AZ：一F(2)P式中△z为河床冲淤厚度；At为计算时间步长；P为床沙干密度。泥沙输移源项F可由冲刷源项E及淤积源项D表示为F=E—D，其中D可由泥沙沉速及含沙量s表示为D=S。对于粘性沙的冲刷源项，Mehta等提出分别适用于弱固结及强固结河床的公式为E=E／exp[(f—r)“]f