二维泥沙数学模型及工程应用问题探讨(赵明登,李义天)

热度125票 浏览73次 【共0条评论】【我要评论 时间:2002年1月01日 12:37

水利论文 x rJ5Xh Ic
二维泥沙数学模型工程应用问题探讨
水利论文/aS3]?kEA.f

赵明登李义天
S$d*?1e0J P*I}8b0(武汉大学河流工程系)

Wd8{uLQhq0

摘要:从泥沙连续性方程和悬移质河床变形方程出发,分析恢复饱和系数的起因和悬移质含沙量垂线分布特点,给出了新的恢复饱和系数;对于动边界问题,通过阻力项的变换,给出稳定性好、精度高的迭代格式;对长江南京上元门河段河床冲淤变化进行了分析计算,结果表明本文所给泥沙数学模型和计算方法合理,能够反映河床冲淤变化规律。水利论文L"Myt\0c9PcpQ i

关键词:恢复饱和系数;动边界;河床变形;泥沙数学模型

#@+jreA+W;[0x8WL,c0

作者简介:赵明登(1963-),男,武汉大学副教授。水利论文)B0yD#^+\O6s-`

1 引言水利论文/](`Ca"[$K

随着计算机存储容量和计算速度的飞速发展,河流泥沙数值模拟也得到长足的发展。一维泥沙数学模型虽然在理论上做了许多概化,但一维数学模型实践经验丰富,比较成熟,目前已能够解决许多长时期、长河段的河床冲淤问题。二维泥沙数学模型虽然在理论上较一维模型完善,又能够给出河床上任何部位的流速和冲淤变化情况,但由于实测资料和实践经验不足,目前还处于发展完善阶段。例如,二维糙率、挟沙力沿河宽分布问题、非均匀沙分组挟沙力及级配问题、恢复饱和系数和动边界问题等都有待进一步完善。本文从泥沙连续性方程和悬移质河床变形方程出发,分析恢复饱和系数的起因和悬移质含沙量沿垂线分布特点,给出了新的恢复饱和系数,并用长江南京上元门河段河床冲淤变化实测资料分析验证。水利论文5r7_k1P6d&W2ZQ

2 泥沙数学模型

$Tqd4A/Y-T;_ EY_0

2.1 水沙运动控制方程

:i1|@%|/n$?w0

平面二维水沙运动的连续性微分方程和运动微分方程分别为水利论文7A)u8W [Qj5v#j3CK

,j=1,2(1)
 i,j=1,2(2)
 j=1,2(3)

式中H为水深,ζ为水位,Uj为垂线平均流速,S为垂线平均含沙量,Ds为泥沙紊动扩散系数。水利论文^5w+DNr.z

2.2 悬移质河床变形方程水利论文-l3vN2q3t"z&~

  在河床表面上z=zb附近取一个微小棱柱体作为控制体研究,假设在床面上z=zb处悬移质时均含沙量为sb,在dt时间内,由于悬移质变化穿过床面z=zb的泥沙质量由两部分组成,一部分由于泥沙沉降引起,一部分由于泥沙紊动扩散()引起,其大小为

水利论文K8h-}6}3fk-q @ o

(4)

]e U n8b8Cn j0
  同时,在dt时间内,由于悬移质运动引起的床面高程变化即冲淤厚度为,在控制体内其质量为

W4gtdtdxdy水利论文!mi"},k_O/Pzv:d ?

(5)

\)M0a&Ew0

根据质量守恒定律,W3=W4,即水利论文,Eb8x&A0f%RKV

水利论文F!d)@P u

(6)

)TCE$?E0

上式为悬移质河床变形方程。当悬移质含沙量处于饱和状态时,悬移质含沙量等于挟沙力,悬移质引起的河床冲淤平衡,河床变形为0,即有水利论文?eA5z]

Y zh Z&o2nF2d0

(7)

9["Sh AEv0
式中为河床底部挟沙力。悬移质河床变形方程可以变为

k#yGh+N Fdhw-L0

(8)

_H nB9y0

在平面二维计算中,假设床面含沙量和挟沙力与垂线平均含沙量S和垂线平均挟沙力S*之间的关系为水利论文~,rW0Y%]9@ A p

b1S,*b2S*

n a_3c)p X-]:o0

(9)水利论文e,s#beO4[$Z%h*j

通常认为α1≈α2≈α,并称α为恢复饱和系数[1,2],最后得平面二维悬移质河床变形方程可以变为

#i(`9s]b0

水利论文~8XR?bZV;k

(10)水利论文{f1t)`+Z)A|,?

3 恢复饱和系数水利论文F As-_/[w

3.1 原恢复饱和系数存在的问题

O?8J$j!zi ow0
  式(9、10)中恢复饱和系数α1,α2,α是由于采用垂线平均含沙量代替河底含沙量而引起的,由于河底含沙量一般是大于垂线平均含沙量的,因此从理论上讲,恢复饱和系数总是大于1的,这一点在有关教材[1]和专著[2,3]中都提到。但是,在实际计算时,根据实测资料推求的恢复饱和系数α往往都是小于1的,有的河流(河段)中,由实测资料率定的恢复饱和系数甚至为远小于1的正数[4]。有人将这一矛盾归咎于床面附近悬沙边界条件的提法和河底挟沙力的引入,但是,上述恢复饱和系数是在河底挟沙力引入之后,由于采用垂线平均含沙量代替河底含沙量而引起的,因此恢复饱和系数的大小应该与河底挟沙力的引入关系不大。笔者通过实测资料分析,认为产生这一矛盾的关键在于假设的α1≈α2≈α,即b/S=*b/S*与实际情况不符。比较式(10)和(8)的差别可以看出,α实际上是由于用垂线平均含沙量与平均挟沙力之差(S-S*)代替河底含沙量与河底挟沙力之差(b-*b)而引入的修正系数。

悬移质含沙量垂线分布规律和公式一般都是在冲淤基本平衡的饱和输沙条件下得出的,因此由此得出的恢复饱和系数只能适用于冲淤平衡时含沙量等于挟沙力的情况。例如,采用莱恩及卡林斯基公式水利论文h M+_)P6MF

**/*b=e-KZ/H,K=6ω/κu*水利论文-vj]7O,`$qT\h

(11)水利论文!qG)t(lM9X0D.U

沿垂线积分平均得恢复饱和系数为水利论文6Q9F2m'?|6b

α2=*b/S*=K(1-e-K)-1

-^^(zK;hQ&E v.NH0

(12)

\i\K X.Y_ bu0
显然,恢复饱和系数α2是大于1的。例如当K=1时,α2=**b/S*=1.58。但是应该注意到,这仅仅是在悬移质泥沙处于饱和状态、河床处于冲淤平衡条件下才成立;当悬移质泥沙处于非饱和状态、河床冲淤不平衡时,由于含沙量垂线分布与饱和状态的分布不同,因而α1≠α2≠α。

3.2 新恢复饱和系数提出

'a Mo1p3\&tdiG)Q0
  由床面附近悬移质泥沙扩散理论和河床边界条件(7)式可知,河床底部挟沙力由脉动项()引出,它与水流的强度有关,与底部含沙量梯度有关。在水流条件一定的情况下,含沙量处于次饱和、饱和、超饱和状态时的垂线平均含沙量差别较大,而底部含沙量差别较小,因此相应的系数不同,即α1=*b/S≠α2=**b/S*。显然,非饱和时的系数α1与饱和时的系数α2差别与悬移质的非饱和程度有关,非饱和程度越大,则系数差别也越大。以S*/S反映其非饱和程度,令α′=α12=(S*/S)β,β的大小一般在0和1之间,可以由实测资料推求。平面二维悬移质河床变形方程可以变为

DZ } z$P S;R3E r0

(13)

c0B3o9Y6{dX0

从以上推导过程可以看出,系数α2只是在饱和情况下由于采用垂线平均含沙量代替河底含沙量而引起的,并非调整恢复非饱和到饱和状态的系数,而真正反映非饱和程度的系数应为α′12。当含沙量处于饱和状态时,α′等于1;当含沙量处于次饱和状态时,α′大于1;当含沙量处于超饱和状态时,α′小于1。因此,严格地讲,应称α′为恢复饱和系数,而α2为饱和条件下用垂线平均含沙量代替河底含沙量而引起大于1的修正系数。

C9na6N5vW [:Z0

4 模型应用水利论文R`%j6xOE

为了检验上述泥沙数学模型和恢复饱和系数合理性,对长江南京上元门河段河床冲淤变化进行了分析计算[5]。动量方程和泥沙连续性方程采用有限分析法离散,水流连续性方程和河床变形方程采用有限差分法离散[5]水利论文 l0}~v_FW ]v)n

4.1 动边界处理方法水利论文S`8U{9{%v l

河道两岸的水边线位置是随着水位升降而变动的,这给数值计算带来了很大困难。河道非恒定流计算网格一般按高水位划分,当水位降低时,边滩露出水面,边滩上的网格处在水边线以外,从物理本质来说,此时边滩上水深、流速应为0,但在一般的迭代计算过程中,如果水深为0或很小时,常会出现计算机数值溢出或迭代发散。目前国内外许多学者常假设一个5~10cm的虚拟水深作近似处理,但由于虚拟水深较大,会引起水量、沙量和边滩冲淤量的计算误差。本文通过阻力项的变换,给出收敛性好、精度高的迭代格式,虚拟水深可以小到10-20cm,能够反映水沙运动的物理本质,新的迭代格式如下水利论文M4ZCg F1[X

水利论文&gZ wC0oV*\ YC(}

(14)水利论文rB8gtA.g bT

式中Cn、CP为有限分析系数。

.A$bP D4}/A4`3L7A0

4.2 非均匀沙处理方法

Ce N+c-GCFC)r n0

对于非均匀沙一般的处理方法是:将非均匀沙按粒径大小划分成若干个粒径组,将每个粒径组都看作是均匀沙,分别满足泥沙连续性方程和悬移质河床变形方程。对于非均匀沙模拟计算的关键问题是分组挟沙力计算问题。对于均匀沙来说,在泥沙粒径一定时,由于河床和悬移质中该粒径的含量都是百分之百,因此挟沙力只与水流强度有关。而对于非均匀沙来说,挟沙力不仅与泥沙粒径和水流强度有关,而且与河床中该粒径泥沙的含量多少有关,因为河床中含量越多,床面附近泥沙的紊动传递和悬富换机率越大。以某一粒径组泥沙在床沙中含量所占的百分比Pbn作为泥沙悬富换的机率,乘以该粒径组泥沙含量为百分之百时的可能挟沙力S′n*,即为该粒径组泥沙的实际挟沙力Sn*水利论文 F h~ FBv9O n

Sn*=PbnS′n*=PbnKn(U3/gωnH)m水利论文 `MM1|G?t)M|

(15)水利论文5@DkQ2D5VFZ

4.3 模型应用

~HP6B$Bb`1\w0

模型验证资料为长江南京河段1997年10月至1998年10月实测水沙过程和1997年10月和1998年10月实测地形[5]。模型率定参数选取为:水流挟沙力公式中系数[6]Kn=0.03~0.07,m=1.0~2.0,恢复饱和系数中β=0.4~0.6。

r-r"qp!` S:a)?0

1分别为上元门附近河段1998年10月计算地形和实测地形比较图,图2为上元门断面计算地形和实测地形比较图。计算初始地形为1997年10月实测地形,1997年10月至1998年7月整个断面均为冲刷,1998年7月至1998年10月深槽有点回淤,深槽与边滩之间河床仍在冲刷。比较结果表明,经过一个大水中沙年的冲淤后,计算地形和实测地形符合良好,说明本文所给泥沙数学模型合理,能够反映河床冲淤变化规律。水利论文:N2Af;?{TnS:y7I\3h

利用本数学模型预测了京沪高速铁路南京长江沉管隧道工程对防洪航运及河床演变的影响[5]。计算结果表明,沉管工程阻水所产生的水位最大壅高为1cm,垂线平均流速最大变化为5cm/s,沉管附近河段左边河床受沉管控制不能冲刷,右边河床最大冲深0.5m,冲刷面积50m2左右。数学模型计算结果与河工模型试验结果相互印证,可以满足工程设计需要。水利论文o%m6FSb+BU$?z

1(a)上元门附近河段1998年10月计算地形图水利论文7|xs \;TZ/E0w
Fig.1(a) River bed contour plan in Shangyuanmen reach (1998.10, computed)
1(b) 上元门附近河段1998年10月实测地形图水利论文Qk!e/a8PDx
Fig.1(b) River bed contour in Shangyuanmen reach (1998.10, measured)

5 结论

T^ N:V"["A/l0

从泥沙连续性方程和悬移质河床变形方程的推导过程可以看出,由于采用垂线平均含沙量代替河床底部含沙量而引入的系数α1、α2或α,从理论上讲,应该是大于1的,而实测资料率定的α小于1。产生这一矛盾的主要原因是,假定了非饱和悬移质垂线分布形式与饱和悬移质垂线分布形式相同,且α12。本文通过分析系数的起因和悬移质含沙量垂线分布特点,给出了新的恢复饱和系数,可以避免上述矛盾且具有真正意义的调整恢复饱和程度。系数β取值可以由实测资料率定。

2P:o{(h|h}wj\0
2 上元门断面地形比较图水利论文L/p:b)UT
Fig.2 Comparison of riverbed at Shangyuanmen cross-section

对于动边界问题,通过阻力项的变换,给出收敛性好、精度高的迭代格式,虚拟水深可以由通常的几厘米减小到10-20厘米;对长江南京上元门河段河床冲淤变化进行了分析计算,说明本文所给泥沙数学模型和计算方法合理,能够反映河床冲淤变化规律。水利论文c$v%p1u%`/h(e;C

参 考 文 献水利论文c2H ar{

[1] 谢鉴衡。河流模拟,水利水电出版社,1988.6.

%`{:ofa,X0

[2] Zhang Ruijin, and Xie Jianheng, Sedimentation Research in China, China Water & Power Press,1993.pp217-218.

.]/C``P0

[3] 杨国录。河流数学模型,海洋出版社,1993.10.

]i9wl;io$j0

[4] 韦直林.付小平。论泥沙连续性方程的建立及相关的几个问题.武汉水利电力大学学报,1997.10,pp14-17.水利论文"N*l5d{IG3j ~,E$M

[5] 赵明登。京沪高速铁路南京长江沉管隧道数值模拟研究报告,武汉水利电力大学,1999.6.

7Y [P2AK D` r.Y0

[6] 长江水利委员会长江下游水文水资源勘测局。南京河段泥沙基本特性分析。1999.1.

1p-W TC1aL lV0
TAG: 工程 李义天 泥沙 数学模型 赵明登
顶:18 踩:25
【已经有82人表态】
25票
极差
7票
很差
5票
较差
6票
稍差
13票
稍好
8票
较好
9票
很好
9票
极好
下一篇:沙漠沙的筛分粒径与粒径计粒径的比较实验(吴世亮,倪晋仁,李振山)
上一篇:泾河东庄水库泥沙模型设计(张红武, 张俊华,王国栋,柴秋峰,张隆荣,姚文艺)
查看全部回复【已有0位网友发表了看法】

广告投放

广告投放