垂线平均含沙量两种表述方法之比较 .docx

工程建筑学相关资料ENGINEERING ARCHITECTURE RELATED INFORMATION 工程建筑学相关资料 ENGINEERING ARCHITECTURE RELATED INFORMATION PAGE PAGE 1 垂线平均含沙量两种表述方法之比较  摘要：表达沿垂线平均含沙量的方法通常有流量法与面积法两种，但在实际运用中常发生混淆。本文讨论了两种垂线平均含沙量之间的关系，指出面积平均含沙量要大于流量平均含沙量。含沙量及流速分布越不均匀，两者差别就越大，当悬浮指标Z≤0.3时两者相差不到10%。最后还对两种垂线平均含沙量的应用作了初步探讨。 关键词：悬移质 悬浮指标 流量平均含沙量 面积平均含沙量 1 问题的提出　　悬移质垂线平均含沙量是泥沙运动力学中的一个重要内容，是研究悬移质输沙率与水流挟沙力的基础，同时也是研究河床变形问题的依据，在生产实践中具有重要意义。　　为了研究、计算上的方便，常需考虑悬移质含沙量的垂线平均值，由于含沙量及流速分布在垂线上是不均匀的，因而存在多种垂线平均含沙量的表达方式。垂线平均含沙量的表达通常有所谓的流量法与面积法两种，它们的定义分别为[1]S1=/(1)S1=1/h·(2)式中 S1、S2分别为按流量法与面积法计算的垂线平均含沙量(以下分别简称为流量平均含沙量与面积平均含沙量)，u、s为距河底y处的当地流量与含沙量，h为垂线水深。　　两种垂线平均含沙量均具有一定的物理意义。前者S1表示的是流量加权平均值，其物理意义是单位时间内通过单宽断面的悬移质输沙率与单宽流量之比，即分部输沙率平均；而后者S2表示的是面积加权平均值，其物理意义是单位时间内通过单宽断面的沙量与单宽面积之比，即分部沙量平均。在实际应用中，两种垂线平均含沙量也各具特色，流量法由于同时考虑了流速和含沙量这两个因素，且流速为向量因素，这样求得的垂线平均含沙量S1，常用于与悬移质输沙率有关的计算中；而面积法仅考虑含沙量这一非向量因素，这样的垂线平均含沙量S2一般多用于理论研究中及一维或二维垂线平均悬移质扩散计算中。　　但在一些关于悬移质泥沙运动的研究文献中，经常会出现将S1与S2混淆，甚至是相互代用的现象。例如在计算悬移质输沙率gs时，常用单宽流量乘以面积平均含沙量S2来推求，实际上此处应以流量平均含沙量S1代替S2。又如在一维悬移质扩散方程中[2]dS/dx=-[αωs/q(S-S*)](3)式中α为泥沙恢复饱和系数，ωs为泥沙沉速，q为单宽流量。式中的垂线平均含沙量S及挟沙能力S*应取S1还是S2，常出现混乱。同样在垂线平均的平面二维悬移质扩散方程中，也常有类似的混乱之处。　　无论是在数学形式上，还是在物理意义上，流量平均含沙量S1与面积平均含沙量S2是两个截然不同的概念，一般不可相互代用。本文试图对这两种垂线平均含沙量之间的差异作一分析，以求能澄清某些概念上的问题。2 流量平均与面积平均含沙量的差异　　垂线流速分布采用著名的普朗特对数流速分布公式u/U*=(1/k)·lnζ+C(4)或写成u/U=［1+1/kC0(1+lnζ)］(5)式中 k为卡门常数，u是相对水深ζ(y/h)处的点流速，是垂线平均流速，U*为摩阻流速，C0为无因次谢才系数(C0=C/，C=/U*)。　　垂线含沙量分布采用著名的罗斯公式S/sa={ζa(1-ζa)·[（1-ζ）/ζ]}Z(6)式中 s为相对水深ζ处的点含沙量，sa为ζ=ζa处参考含沙量，Z=ωskU*为悬浮指标。式(5)及(6)所表示的流速及含沙量垂线分布规律，尽管还存在一些不足之处，但由于两式的计算精度较高，且有一定的理论基础，因而得到了广泛的应用。利用式(5)及式(6)，可得流量平均含沙量S1=∫10usdy/∫10udy=sa[ζa/（1-ζa）]Z∫10{(1+1/kC0)·[(1-ζ)/ζ]Z+1/kC0·[(1-ζ)/ζ]Zlnζ}dζ(7)若令∫10[(1-ζ)/ζ]Zdζ=J1(8)∫10[(1-ζ)/ζ]Zlnζdζ=J2(9)则式(7)变为S1=sa[ζa/(1-ζa)Z]{(1+1/kC0)J1+1/kC0J2}(10)　　又利用式(6)可得面积平均含沙量 S2=∫10sdζ=[ζa/(1-ζa)]Z·sa∫10[(1-ζ)/ζ]Zdζ(11)同样将式(8)代入上式可得 S2=sa[ζa(1-ζa)]ZJ1(12)　　流量平均含沙量与面积平均含沙量之比为β=S1/S2=1+1/kC0(1+J2/J1)(13)由于J1、J2仅与Z有关，因而β也仅为Z与kC0的函数，即β=f(kC0，Z)(14)　　从数学上可以证明，函数J1与J2具有如下性质：①J1＞0，J2＜0;②｜J1｜＜｜J2｜。所以有β=S1/S2＜1(1