非恒定流泥沙数学模型原理及其应用(梁国亭,高懿堂,梁跃平,龚坚)

热度179票 浏览38次 【共0条评论】【我要评论 时间:1999年7月01日 14:27

非恒定流泥沙数学模型原理及其应用

'?aaA.G3y.C0

梁国亭高懿堂梁跃平,龚 坚

C/C#g"m1R:eVJK2N0

(黄委会黄河水利科学研究院; 三门峡水利枢纽管理局)

3}H%c%T8Q6J0

摘要:本文介绍了一维非恒定泥沙数学模型的原理和计算方法,并且对于泥沙数学模型中的一些关键技术等问题进行了较详细的讨论。利用已建立的一维非恒定泥沙数学模型,对黄河小北干流1981年洪水资料进行了模拟计算,计算结果表明:黄河干流、渭河、北洛河计算水位、流量过程线与实测值符合良好,可用于黄河的洪水演进计算。水利论文[4Wn"S lfCl

关键词:非恒定流;泥沙数学模型;河床变形

!_d4Pmi#d p_6}0

基金项目国家自然科学基金和水利部联合资助重大项目(59890200)

K+h ZT]2z0

作者简介:梁国亭(1963-),男,河南省淅川县人,黄河水利科学研究院高级工程师。水利论文/{ _VN8~O

1 泥沙数学模型基本方程

'JSYB4az{)p9TfjZ0

  明渠或天然河流常被考虑作为一维流动,根据洪水波运动的圣维南方程、泥沙连续方程和泥沙扩散方程,可以简化推导出一维非恒定流泥沙数学模型的基本方程,其形式为水流连续方程

8k a+k/y7g,On&tw1d0

990408e1.gif (1111 bytes)

2oYH#|q6x\8G;Z0

(1)

$@d&LWqm#v6Y0
水流动量方程

990408e2.gif (1615 bytes)水利论文 NQ2YcNmx

(2)

7L]/v"P!Y(e V0
泥沙连续方程

990408e3.gif (1317 bytes)水利论文`og,~BwDT

(3)水利论文x8De wl

不平衡输沙方程

a1tH"uf|h$~0

(4)

_xbj-K7q_H8L.R0
挟沙力方程

S*=f(Q,A,B,ω,S)

%@n@z_)S5b8H.J}0

(5)

m2T![3_W3B0

其中 Q为流量,A为断面面积,B为断面宽度,Z为水位,K为流量模数,S为断面平均含沙量,S*为断面平均挟沙力,g为重力加速度,α为泥沙非平衡恢复饱和系数,ω为泥沙颗粒沉速,Ad为断面冲淤面积,γ′s为泥沙干容重。水利论文uN"k1\0v|,{S8o

2 计算方法简介水利论文!\%fP'm.\ W+b}C

  一维非恒定泥沙数学模型的计算采用非耦合方法,首先求解水流连续方程和动量方程,然后求解水流挟沙力、泥沙不平衡输沙方程和泥沙连续方程,具体求解过程如下。水利论文(F%[ qb*~4vi'w

2.1 水流方程的求解水利论文K7`5\r;z jS

  首先利用Preissmann隐式差分格式将水流连续方程和动量方程转化为差分方程,然后对差分方程进行线性化,在线性化过程中,略去增量的乘积项,最后得到以下线性方程组水利论文cV+@G N7]N$k IQ

A1jΔQj+B1jΔZj+C1jΔQj+1+D1jΔZj+1=E1j水利论文We(P3P [%UQ

(6)

6zc~ VD2c#K0G,l:u0

A2jΔQj+B2jΔZj+C2jΔQj+1+D2jΔZj+1=E2j水利论文5P#]W3AHx)z.hM

(7)水利论文.{ m ^\$TU

其中 Aij、Bij、Cij、Dij、Eij(i=1,2)为第j单元河段差分方程的系数(j=1,2,......,N-1,其中N为断面个数)。水利论文DsEvG#],v3n

  给定边界条件水利论文/Zy9jZ2}%P4kxI

ΔQ1=Q1n+1-Qn1=Q1(tn+1)-Q1n

2MA%F/b:mn,}/l q!d-p0

(8)水利论文b4GwCdp8U`

ΔZN=ZNn+1-ZNn=ZN(tn+1)-ZNn

v+z UZ$B0

(9)水利论文'Q4{L7h\0y

方程(6)、(7)式及边界条件,共有2N个未知数,2N个方程,可以求解。由于差分方程中的系数包含有未知数,方程求解不能直接求出未知变量,因此方程求解时必须进行迭代处理。下面给出用追赶法求解的步骤,追赶方程为

9eg$p8T%B#hn g _0

ΔQj=FjΔZj+Gj

&[.],[ U{ dz$W&H([0

(10)水利论文 o g `\S e"S!?D~

ΔZj=HjΔQj+1+IjΔZj+1+Jj水利论文\ {&A8z ~L%T.\ v

(11)

-GF\,r7N g0

其中Hj、Ij、Jj、,Fj、Gj为追赶系数。

e`+x;N[G&`7e0

2.2 水流挟沙力、动床阻力、河宽变化、床沙级配调整采用文献[2]的方法水利论文 t1v2Y s7[Qe/[

2.3 不平衡输沙方程求解水利论文%i Kl `'iyWz

  利用迎风格式,将(3)式离散为差分方程,整理后得水利论文:C0b_ B$jE

990408e12.gif (3380 bytes)

5B.U i,j/m3i0

(12)水利论文 IyBB cu'db

当Q≥0时,利用上边界条件,自上而下计算各断面含沙量;当Q<0时,利用下边界条件由下至上计算各断面含沙量。

e'K.Q3R Z:vh0

2.4 河床变形及淤积量计算

"g3Fd4Ma2d9q#M0

  水利论文8g}^a.u;J uY

    由式(3)与式(4)相减可以得到河床变形方程为

2YmI`#`P0

990408e13.gif (1174 bytes)水利论文 x1?"s)vP{M

(13)

T]#V6w!v0
    将上式写成差分方程,各断面淤积面积为

990408e14.gif (1352 bytes)

gJy(hlDMt'z0

(14)

jK6P"]`7n8_HB0
    第j河段的淤积量为

ΔWj=(ΔAj+ΔAj+1)Δxj/2水利论文6l Ge*Mq4]6T"tO H

(15)水利论文0^xzvl L0J

3 水流内边界的处理水利论文2U)HVi Vs)LeK

  水流内边界是指河道的几何形状的不连续或水力特性的不连续点。例如,河流的汇合点、河流分流、局部河段内生产堤决口等等。在这些内部边界处,圣维南方程组和单一河道泥沙不平衡输沙方程等都不再适用,必须根据其水力特性作特殊处理。内边界条件通常包含两个相容条件:即流量的连续条件和能量守恒条件(或动量守恒条件)。本模型主要考虑了以下几个类型的内边界处理。

^(A A^"j.gG0

3.1 水沙的汇入或汇出

8g7a f~ {5s.\8I0

  如图1所示,假设汇入或汇出点上下断面满足以下条件水利论文&ig M1XKJieC

990408e15a.gif (1001 bytes)水利论文`%m%E w i"K'w)L^@

990408e15b.gif (1065 bytes)

_JLt1n0

990408e15c.gif (1214 bytes)

,e2u6rv$r%g5N k0

3.2 支流从干流分流

;qK @+O4M*b1LZ0

如图2所示,干流和支流上断面之间应满足连续方程和能量方程水利论文5~yO{~|9}}1v7P

990408e15d.gif (1078 bytes)

"cbN{p3sS0

990408e15e.gif (1688 bytes)

1wHc$zM4dM.]L0

990408e15e.gif (1688 bytes)

H:|3}5xID~0

990408t1.gif (1665 bytes)水利论文oG s!ec

990408t2.gif (1796 bytes)

-KF;aU)RG E$Q/I3v0

990408t3.gif (1804 bytes)

Zax;Kbasg%J6n0

图1 水沙的汇入或汇出示意图
b]9nLm0
Sketch of water/sediment inflow or outflow

2H&m{O n$d#w ]0

图2 支流从干流分流示意图水利论文*Wa'\9a+Zm
Outflow from the main stem水利论文|9BPj6Ue

图3 支流汇入干流示意图
&v5M.O4h@4~(Q0Inflow into the main stem水利论文)q1tAu)[#X.M

  支流从干流分流时,将干流分流断面按干支流流量比分为两部分,忽略时变项,将方程(4)式直接写成差分形式,求得干流和支流下游断面的含沙量为

2D(t0QYM4](i0

990408e15g.gif (1950 bytes)水利论文Sh/@ln]r

990408e15h.gif (2027 bytes)水利论文t2S AH:Z$kpU

其中φ为分流系数,由实测资料确定。水利论文FGHs@1S9Q(M-j w

3.3 支流汇入干流

,t!V)\"e*f0

  当支流汇入干流时,与支流从干流分流类似,干流和支流上断面之间也满足连续方程和能量方程水利论文2w(`vP#a*L.S,T

990408e15i.gif (1074 bytes)

+_7F"im$\EL0

990408e15j.gif (1679 bytes)

K`+?W's*y k0

990408e15k.gif (1688 bytes)

&@Yn(X.D:e'e3E:d]0

  将汇流断面按干支流流量比分为两部分,忽略时变项,将方程(4)式直接写成差分形式,求得汇流断面的含沙量为

^2Y!dz9aP:PML0

990408e15l.gif (1932 bytes)

`#G3y4lsp0`+?0

990408e15m.gif (2041 bytes)水利论文1BH B3s:BQ9C$hl;b

990408e15n.gif (1327 bytes)

+^'rJ+f-n"A0

其中 φ为分流系数,由实测资料确定。水利论文9~IA_turF.X

4 泥沙数学模型的应用

w z Q2G-I[\$V0

4.1 计算区域和时段水利论文rz&bW)M}X

  计算区域为龙门、华县、河津、状头至潼关河段,黄河干流上有渭河和汾河汇入,在渭河上有北洛河汇入。本模型同时模拟黄河干流、渭河、北洛河三条河流的洪水演进过程,各汇入点作为内边界处理,汾河仅作为已知水、沙过程线汇入黄河干流。上边界条件为龙门、华县、河津、NFDA1头各站,水沙为进口控制的已知条件,下边界条件为潼关站出口控制水位。

r4v*xf}w0

  计算时段为1981年汛初第一场洪水,洪水时间为7月3号至7月14号,洪水持续时间为12天,龙门最大流量为6400m3/s、最大含沙量为298.0kg/m3;华县最大流量为970m3/s、最大含沙量为117.0kg/m3

)mjcp8j? sP9i [+IN0

  原始大断面资料采用1981年汛前实测大断面资料。水利论文U*u!Ok v8D

4.2 计算结果与实测值的比较

jDV1O`L&] R.M0

  图4为潼关站计算流量与实测值的比较、计算含沙量与实测值的比较。由图4可以看出,计算的潼关出口流量过程线与实测过程线比较符合,计算洪峰最大值和相应洪峰传播时间与实测值比较接近;计算潼关出口含沙量过程线与实测值也是比较一致的。由图4分析得出,本模型能够比较好地模拟出三条河流的洪水传播过程和泥沙冲淤调整过程。水利论文C,]7eGz1]c-k C

990408t4.gif (6843 bytes)水利论文*c3B/T PK N:a

图4 潼关站流量、含沙量计算值与实测值的比较水利论文2ZA6VR5OPy"Tl&v3X^

Comparison of calculated and measured discharge and concentrations at Tongguan station水利论文 e;H7Q(s(m,_a

990408t5.gif (5773 bytes)水利论文q:b9m{ OCSo]

图5 华阴站流量、水位计算值与实测值的比较

4Q%W N}S#`v$u u+`0

Comparison of calculated discharge and water levels with measured at Huayin station

!s$gT4_+e4]6Sb7SXv0

  图5为华阴站计算水位与实测值的比较、计算流量与实测值的比较。由图中可以看出,计算水位、流量过程线与相应实测过程线比较符合。同时也看出,在洪水初期,华阴站出现了倒灌现象,流量出现了负值,本模型也比较好地模拟出了黄河干流倒灌渭河的现象。

x#o+p~3YW)t H0

  图6为朝邑站计算水位与实测值的比较、计算流量与实测值的比较。可以看出,计算水位、流量过程线与相应实测过程线比较一致。由于计算北洛河各断面的水位和流量,是通过一系列内边界条件计算得到的,因此也可以说明本模型所采用的内边界处理技术是比较合理的。水利论文!Q#? y~`rr6F-TF

5 结语

6@pwQ-i1sB(S/U?f0  本文较详细地介绍了一维非恒定泥沙数学模型的原理和计算方法,并对泥沙数学模型中的一些关键技术等问题进行了讨论。利用一维非恒定泥沙数学模型,对黄河小北干流1981年洪水资料进行了模拟计算,结果表明:模型比较好地模拟了黄河干流、渭河、北洛河洪水传播特性,洪峰传播时间和洪峰最大值与实测值比较一致;模型对于汇流区水流和泥沙计算所采用的处理技术也是比较合理的,对于干流倒灌支流或支流倒灌干流的现象,模型也是比较能够适应的,因此,模型可用于多条河流的洪水演进计算。水利论文5H+j#h F&w1c7M*M

990408t6.gif (5810 bytes)水利论文,T,PO i c\ lmp)k

图6 朝邑站流量、水位计算值与实测值的比较

1}Y\(D9p7~/D4m"o&O$p#U0

Comparison of measured and calculated discharge and water levels at Chaoyi station水利论文,O:U+yd'Eo X

参 考 文 献水利论文*]0\-`9W"v(Pmbp,k,?

[1] Lyn,D.A. Goodwin,P. Stabling of a eneral Preissmann scheme. J. of Hydra. Eng. Vol. 113,No.1,1987.水利论文C7JOa}7r6_!SF

[2] 梁国亭,张仁。黄河小北干流一维泥沙冲淤数学模型,人民黄河,1996年9月。

&q2P2Xvv ~0

 

O0zX&]'_ l-ux0
TAG: 数学模型 高懿堂 龚坚 梁国亭 梁跃平
顶:27 踩:21
【已经有131人表态】
15票
极差
18票
很差
13票
较差
19票
稍差
15票
稍好
10票
较好
20票
很好
21票
极好
下一篇:高含沙洪水输水输沙特性及对河道的破坏作用与机理研究(张林忠,江恩惠,赵连军,王…
上一篇:冲积河流悬移质泥沙与床沙交换机理及计算方法研究(赵连军,张红武,江恩惠)
查看全部回复【已有0位网友发表了看法】