颗粒相对隐蔽度对泥沙输移的影响.docxVIP

颗粒相对隐蔽度对泥沙输移的影响 1 维泥沙颗粒相对隐蔽度的概念 在研究沉积物的流动速度时，我们需要考虑是否暴露于河床的位置或河床的表面。如何定义标准床的表面（基面）。为回避这个难题,1965年,长江科学院在四川五通桥茫溪河野外水槽试验中,韩其为等研究人员提出用(绝对)暴露度Δ来表示所研究颗粒在床面的位置(图1(a))。它直接反映了颗粒滚动时力臂的大小,而这个力臂恰好反映了所研究颗粒滚动起动的临界条件和滚动速度,还反映了颗粒起跳时的碰撞效果。Δ愈小表示暴露愈充分,在实际应用中,常用Δ′=Δ/D来表示,即为(二维)相对暴露度。自Einstein(1950)提出用隐蔽系数概念之后,隐蔽度、暴露度和隐暴度等名称相继出现,其实它们都是表达同一含义。为从物理概念上理解方便,本文尝试对以上名称做一个统一:①暴露度为床面泥沙颗粒在水中的裸露程度;②隐蔽度为由于受其他颗粒的影响,所研究颗粒不能完全暴露在水中的程度。不难推断泥沙颗粒的暴露度大,则其隐蔽度就小,更易起动;反之,若其暴露度较小,则其隐蔽度就大,更难以起动。 近50年来,无论是在河流泥沙起动还是在风沙起动研究中,泥沙颗粒相对隐蔽度的应用得到很大发展。但是为研究方便,我国大多研究者都将泥沙颗粒概化为二维“圆饼”模型研究颗粒隐蔽度,即二维泥沙颗粒相对隐蔽度。然而泥沙颗粒实际上多可概化为球体,为研究泥沙颗粒的二维“圆饼”状与三维球状概化对泥沙颗粒隐蔽度的计算有何影响,2011年长江科学院提出了三维泥沙颗粒相对隐蔽度概念。即床面上所研究颗粒的最低点到与下游颗粒接触点或者下游2个颗粒交点之间的一个“虚拟点”之间沿床面的垂向距离(隐蔽长度)与所研究颗粒粒径的比值,记为η,同时可定义ξ=1-η为相对暴露度。不仅如此,2011年长江科学院又提出并解决了受床面坡度影响下的泥沙相对隐蔽度最小值的问题,即非均匀沙隐蔽度的最小值不是从0算起,而是至少从(1-cosθ)/2算起,其中θ为床面坡度,只当θ=0时,相对隐蔽度的最小值才可能取0。下面就基于不同维度的相对隐蔽度、隐蔽度的分布规律等相关研究进行简要阐述。 2 相对隐藏度研究 2.1 维泥沙颗粒隐蔽度的计算 如图1(a)所示,D为所研究颗粒直径,O为该颗粒的球心,E为该颗粒沿床面垂向的最低点,J2为与下游颗粒的交点,B为该颗粒与下游颗粒的交点J2在OE线上的垂足,O1和O2为下游颗粒的球心,θ为线OO2与线OE之间的夹角。根据二维泥沙颗粒隐蔽度的定义,可以得出均匀沙隐蔽度η的计算式为 η=BE2ΟE=D-Dcosθ2D=(1-cosθ)2。(1) 不难论证θ?[π/6,π/2],由此推得表层泥沙η的取值区间是[1/2-√3/4,1/2]。 如图1(b)所示,O′是O沿河床垂向在3颗“次表层”泥沙球心O1,O2,O3所构成平面上的投影,BJ2与线O′O2平行,其余符号含义与图1(a)相同。根据三维泥沙颗粒隐蔽度概念,同样可由式(1)得出均匀沙隐蔽度η,并不难论证θ?[arcsin(√3/3),π/2],由此推得该定义下η的取值区间为[1/2-√6/6,1/2]。 根据上述分析,可以发现三维泥沙颗粒隐蔽度在数值上较之经典的二维泥沙颗粒隐蔽度为大。显然,由于在三维模型中,“次表层”颗粒会为所研究的“表层”颗粒提供一个更深的“凹槽”来稳定该颗粒。这里也有必要说明:当泥沙颗粒隐蔽度大于 0.5时,则所研究颗粒处于“次表层”,该颗粒不会瞬时起动,在统计意义上,无须研究泥沙颗粒隐蔽度大于0.5的情况。 2.2 非均匀沙隐蔽度 由于缺乏丰富的实测资料检验,泥沙颗粒相对隐蔽度分布形式至今没有得到很好的解决。即便均匀沙相对隐蔽度分布形式也没有一致的结论,非均匀沙隐蔽度分布问题更是未知。对于均匀沙来讲,Paintal、韩其为等一些开创者都认为颗粒隐蔽度在理论上应该符合均匀分布,并把这一假设推广到非均匀沙隐蔽度研究领域中。均匀沙隐蔽度是非均匀沙的特殊情况,目前本文仅能对比较简单的均匀沙颗粒相对隐蔽度的分布问题进行一些讨论。为简单起见,在此以韩其为定义的二维相对隐蔽度为研究对象。 2.2.1 基于隐蔽度密度函数的密度函数 由于泥沙颗粒隐蔽度的取值区间为[1/2-√3/4,1/2]。韩其为认为泥沙颗粒等可能排列,即认为泥沙颗粒隐蔽度符合均匀分布,进而可以求出对应的泥沙颗粒隐蔽度密度函数为 f(η)={112-(12-√34)12η12-√34?0其它。(2) 2.2.2 密度函数的性质[2] 分析公式(1),如果理论上认为泥沙颗粒隐蔽度符合均匀分布成立,那么也没有理由否定“角度θ在[π/6,π/2]也有可能符合均匀分布”这一论点。由于颗粒排列是随机的,不同的排列对应不同的角度,若假定角度θ对应的区间上符合均匀分布,即θ的密度函数为 f(θ)={1π2-π6=3ππ6θπ2,0其它。