长江中游散乱水位流量关系单值化处理方法研究(戴清 韩其为 毛继新)

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

长江中游散乱水位流量关系单值化处理方法研究

^cX:V(rM Pb0

戴 清 韩其为 毛继新

'q#{ v5Lf0

(中国水利水电科学研究院)水利论文 N? d"i@v#o+a j#y

摘 要本研究在重点分析荆江洞庭湖防洪重要水文站监利、城陵矶(七里山)、螺山、武汉关水位流量长系列实测资料的基础上,初步建立了考虑河段水面比降因子模拟水位流量的单一关系模式:即J=J0+JQ。该式待定系数意义明确,简便易求,并可借此对河道水流特性进行深入分析,利用该方法所求结果与实测值完全吻合。

FKF*r2x L]`0

关键词长江 水位 单值化

Q\-k gl/M!ot#g0

1 问题的提出水利论文 w%tZXjF(I:dc@

  长江荆江河段水位流量关系变化既受荆江裁弯河床演变、洞庭湖分流入汇变化、上游葛洲坝兴建,又将受未来三峡水库建库的影响,变得十分复杂。荆江包括洞庭湖区行洪水道,是长江中游防洪最严峻的河段,水位流量的变化,尤其是城陵矶(洞庭湖出湖七里山站)、监利水位流量关系是荆江与洞庭湖防洪研究的重要特征参数。

$O.| t{K R8G0

  近来研究水位流量关系特性采用的方法主要有:

8J,s8Qt"Sz0

  (1)综合落差指数法:对城陵矶七里山站流量进行单值化处理,公式形式为水利论文"VF5kw)i&IWs

Q=Qm/ΔHβ=f(H)

R |!H"s5Vc3Y,p0

(1)水利论文;R2pmVN/Uj

ΔH=α1ΔH12ΔH23ΔH3

8c:\3@(T@0

(2)水利论文jUg[.^9i#u

式中Qm为实测流量,H为七里山站水位,ΔH1、ΔH2、ΔH3分别为湖区鹿角至七里山、七里山至上游监利、七里山至下游螺山的落差,α1、α2、α3为相应的距离(km),β为待定指数。选取典型年作各大时段单一水位流量关系。

e#w$ZOl-`A0

  (2)幂函数拟合法:依据裁弯前后的一部分水文资料,分别建立在某一顶托流量下七里山站、监利站水位流量关系,采用幂函数拟合,公式形式为

&Bix5}V0

Z=αQβ水利论文0DBB-Z$H

(3)

6o}9R:])L0

式中Z为水位,Q为流量,α、β为待定系数,不同流量,α、β值也不同。

&JUO2h(]|^0

  研究者认为,上述方法一般相关关系较好,除个别点外,相关系数可达到0.9以上。水利论文/cAc{FSc2Zn4Y

  本文研究对相应河段乃至下游螺山站、武汉关站水位流量关系进行了深入分析,提出关系模式成果,将克服上述两个方法不足。并(1)为仅需一站水位资料。(2)变分级流量关系式为统一的关系式。

u/Aefl#F1z;^$k0

2 分析方法及基本关系式确定

ww E ^%l"r5g%Tb![0

2.1 水位随流量变化特性水利论文aUbL-H

  点绘监利站1980年至1987年逐日水文关系,如图1所示。由图可知,监利同流量下水位变幅最大可达4~5m,这种散乱关系,是受下游洞庭湖出流影响的结果,同流量条件下,下游水位高,该站水位就高,反之亦然。既使在关系点相对集中的流量范围,如15000m3/s流量以下,流量由4000m3/s增加到13000m3/s时,水位相应由21.2m抬高至27.8m,变化幅度达6.6m。这种变化的主要因素也是下游水位抬高引起的。下游水位变化小时,流量的变化引起上游水位抬高的变化要小的多,以下游螺山水位20~21m为例,流量变化在3960m3/s至12900m3/s,监利水位由23.4m增至25.8m,变化仅2.4m。

Fp wP0S [2X$e#H0

如将下游螺山站的水位按大小划分若干组,则可知每一组的监利Q、H均存在着图2所示特性。水利论文$f WrH.V,R(U

990310t1.gif (2402 bytes)水利论文 v&T;KR2U4E'`gK

图1 1980至1987年监利站逐日水位流量关系

a,Pt#}.z.F4V6i0

Relation between water level and dischargeat Jianli station from 1980 to 1987水利论文;Z+^3A;|G#T3b`+WH

990310t2.gif (2227 bytes)

4dJ%}0R@J,H_0

图2 监利站水位流量关系(H螺山=20~21m)

g9} YA@:f0

Relation between water level and dischargeat Jianli station (Hluoshan=20~21m)水利论文C&A:q,u#n*i

分析认为,监利的水位由以下两部分组成,一是图2中的点延长线至Q=0的交点,记为H0;二是因流量的变化引起的水位变化,记为ΔHQ,即水利论文/e-|$nBf~

H=ΔHQ+H0水利论文m&aV5Pc]9pu`

(4)

l1Vl~1a0si P0

2.2 ΔHQ、H0的确定水利论文O&Pc%q-~B$~U

  将监利站统计年份逐日流量与监利至螺山水面比降关系点绘在图3中,点群是很散乱的,同样采取下游水位分组方法,如H=20~21m组点,见图4“+”点所示。H=27~28m组的点也一并点绘于图4中,用“*”点表示。水利论文&pfq.Kc ~

990310t3.gif (2734 bytes)

{Gd o `5i;RO0

图3 监利流量与监利至螺山水面比降关系水利论文 O5tQYt Gh m9O

Relation between discharge at Jianli station and water surface slope in reach from Jianli to Luoshan水利论文en+Wl NeFS

990310t4.gif (2329 bytes)水利论文|F _ ^"E Xl

图4 不同螺山水位监利流量与水面比降关系水利论文Nn2R;w2s+Z(n1A

Relation between discharge and water surface slope at Jianli station under different water levels at Luoshan station水利论文6b1\8W9C \4p ? ?G5D

 由图中点可知,图3中很散乱的点是由下游水位不同组数据汇合而成的,实际每组的Q~J关系是很有规律的。同组条件下,比降随流量的增大而单调增加。进一步分析可知,J的增加与Q可由线性关系表示。如将Q=0交点的比降值用J0表示,流量引起的比降变化用JQ表示,即可得到

9T1z uo(O:x0

J=JQ+J0水利论文L,Rx4C+p6s} I

(5)水利论文i+M7X$m(n;~

对比(4)式可知,监利水位ΔHQ、H0分别可表示为JQL、J0L+H(L为两站间距,H为下站水位)。即有

TvM7Kst0

H=JQL+J0L+H水利论文+|:tct5t"I6\*@i

(6)

WdU*y%N*C'?2|q I0

2.3 经验拟合

#n!G`z h"kW[ f0

  由上述可知,在同组数据J与Q呈线性关系,即可表示为

9HY&| P}k4O0

JQ=φ(H)Q水利论文?B*dN%@-bp

(7)水利论文(K7h^-hY ?-m5Y*b

分析数站实测水文资料,研究后发现(H)、J表达式分别为水利论文 xH+GsKf!d^4~9[

(H)=10aHb+c水利论文 j6|{%cS

(8)

N)P1Kf ]!zh0

J0=d+eH水利论文z!I6S6v)^

(9)

D @U#Y cL0

将(8)、(9)式代入(5)、(6)式中,即为本文建立的水位流量关系基本表达式水利论文O6q+opJ4^ jk

H=(10aHb+c)QL+(d+eH)L+H水利论文-iy&?Gr"S

(10)水利论文A:GV9o:KDZ9x)WeHFv1[

水利论文U1iCpn

J=(10aHb+c)Q+(d+eH)

j"a0WWJN$Y0

(11)水利论文F*atGB7}

其中 a、b、c、d、e五个为待定系数,且a、b、c与d、e互不相关,每个系数均有明确的数学含义和其简单的计算方法,后文中将详细介绍。其他符号意义同前。

/K1Xhv_bD0

3 实例计算水利论文x]h&^0L

3.1 系数确定水利论文 G CuQPc%q)U^F:I

  水文资料按监利1980~1987年下游螺山站水位大小分组,计算用两次最小二乘法,分别对每组数据点拟合求J、φ(H),结果点绘于图5及图6中。

Z p t4l"l7zrW2s0

990310t5.gif (2333 bytes)

s(B j8x?/j0

图5 螺山水位与监利J关系水利论文;c RZ2lOL EGZ

Relation between water level at Luoshan and J0at Jianli station水利论文"o#h"p*n5~w1]

990310t6.gif (2199 bytes)

w e4A+_8`H0

图6 螺山水位与监利φ(H)关系

%GU5C#C.d;G0

Relation between water level at Luoshan and φ(H) at Jianli station

F6V`%xJ,p0

  由上述图点关系可分别求得a、b、c、d、e五个待定系数(比降J用1/1000值,流量Q用103值,距离L用公里数)。系数值列于表1。以同样方法,对七里山站(下站螺山)、螺山站(下站龙口)、武汉关站(下站黄石)1980年~1987年逐日水文资料进行分析,确定各站系数并一一列入表1中。

S Fk JO0

表1 系数取值(1980~1987)

Eq;N9e i/Z0

Values of coefficient(1980~1987)

IzxJx]0

站名

9C:Q/X$E)ul0

监利(螺山)水利论文J L`"?Lq?

七里山(螺山)水利论文 mQ!I,a/D&S"G6L

螺山(龙口)

i5Chc;M*b1S0

武汉关(黄石)

`@]I(F u0

a水利论文-E8W @/SlrI)U

4.22

'Hk;\Z(uK0

3.7634

&Pxa)S5c,q$g0

4.6267水利论文9_"y8N*uU-Q5J

-0.85460

:^~-j'F[1X9C0

a2

|H_'Uc-d0or0

-5.333

u#G/x1iS];p*q6i!y#S0

-5.2259水利论文5w,xz[M a _~g)n

-6.1998水利论文+e-P9T[L:w/q6T

-1.9122

t!q1tm _F*I0

a3水利论文+BI|0dNv1A

0.000308水利论文-kU.}U%h

0.000043304

n4ZZP$Rt"r;N0~M0

0.00021897水利论文-r W2Y]&e:W%r

-0.000026073水利论文7Q%pLp$J;V

a4水利论文b'Z/W;x*d'H(d

-0.00061278

5Z Jd9M)eo4t?P1Y0

-0.00071889

k(O Rr I \0

-0.00062114水利论文I;g `1wu

0.00055556水利论文^!] SI7GDa

a5

$V/X1Nme7x0

0.031678

YtRv(yo qlr0

0.049650

:p&GTY| `$|/`}`0

0.028696水利论文 O3t/[h'G

0.0006667

7a:nf9?@3VbMv*v"o0

3.2 验算结果水利论文 ?c8|y+x&JY)N

将统计年份的点计算对比结果点绘于图7(a、b、c、d)中,可见,计算与实测点完全吻合。水利论文]YY4rKn"aBk,A

综上分析,得出以下认识,(1)无论水文站流量水位点分散或单一,该站水位高低主要受下游站水位变化的影响;(2)利用上述(10)式关系,当已知本站流量、水位、下站水位三个参数中的任意两个,均可计

"ho)F-e*u `F0

990310t7.gif (6106 bytes)水利论文2d&RA$lA|:CQ [

图7 实测水位与计算值对比

i4W@WE,d5A0

Comparison of observed and computed stages

9|p+Dm!@ C"N(d_ {0

990310t8.gif (2345 bytes)水利论文zkFo!vmw

图8 监利站实测流量与计算值对比

T"_2XU t-qkXc)r0F0

Comparison of observed and computeddischarge at Jianli station水利论文(D1OG2W5K%g Uz

算出第三个参数值。仍以监利站为例,由监利站水位及下游水位计算得到监利流量,并与实测流量对比,点绘于图8中,点子落在45°线上;(3)天然河道均匀流概念很难定义,而由均匀流基本公式J=Q2n2/A2R4/3;仍看不出公式中各参数之间的因果关系。实际上,Q、n、J均是独立量。如果边壁条件不一样,水深或水位变化将引起n值的变化,反映河道沿程阻力特性,可认为流量变化仅引起水深(水位)变化,但对n值不应有所影响。本文推荐的水位流量公式中包括与流量无关的比降项,即J0(=d+eh)。天然河道中均匀流时,Q、n值无法确定,仍有待今后研究,但该项能够具体计算,确实对该问题的解决大大前进一步,这方面研究有待进一步深入。

rh(e5BW5[4USX2t0

4 监利水位年代变化分析水利论文6G}!fe1P Tq

  建立了上述水位流量关系式,即可很方便地对比分析不同年代各种水文组合条件下的变化特点。

qwII8J*j&]s0

4.1 水文频率分析

#B1r6|G{)OV!n0

  受实测资料所限,与80年代相对比时段选为1954~1959年系列,对监利站水位流量及下游螺山站水位进行频率分析,结果列于表2。水利论文 Y!q*a(COecV

  由表中数据可看出,50年代与80年代相比,同频率监利流量增大,水位反而有所下降。其原因之一是螺山水位同频率水位也是在下降的。

ic,xp{wV$\-J$J \0

4.2 监利流量同频率条件水位变化

!z3DgP\ @F2?:O0

对于表2中监利各频率流量在不同下游水位条件下用水位流量关系计算出本站水位变化,如图9(a,b),图10(a,b)所示,由图可知:水利论文4XHHJ*kgy2y

(1)不同年代出现最大流量时,下游同水位条件下80年代水位具有高于50年代水位的趋势,但实际上因下游螺山水位在80年代低于50年代水位,因此监利水位在80年代低于50年代。水利论文.gp;mi9P5i'Bq

表2 监利流量水位、螺山水位频率分析对比

1?&s"jS/aK.U#['Y ? G0

Analyses of frequency of discharge and water level at Jianli and the water level at Luoshan水利论文:L^.dA$R M


水利论文 P-Yj&u"IT

年代

/[(M8n9QV.nh r(\0

频率

+Z+m9CD0tg0

Max

?*QY o nK T0

1%

"e0a.W@$~ g8d1Wx0

5%水利论文$]_,p$};v1JS0a

10%水利论文{yo&mm8f?g([P

50%

B s'V1A-v`2@1B1o Z0

min水利论文V5zFgS!mJ


监利水利论文 om;Me.U

1980~1987水利论文8T Zo-~A

Q水利论文MU;CBV

35200

'wT~;N9W0

29400

o&Ho(M.Z ~0

23200

!t+?,Q!G)@?0

18800

;?;s(v CJ%w0

7810

D2is WDe0

3140

7ME:b/AnU/^*c0

H

+Q8m-f1J0W*Z.D?0

34.7水利论文J tS)aIG

33.89水利论文5Gn2wX$j#WF_3W_

32.31

Y:|N NLu0

31.26水利论文.I Xc/hR A3M k)i

25.87

'D V\2s-R"GUK0

22.08

K"C4d I/jUt6q0

1954~1959

f?Q)a,zb9hO+{0

Q水利论文+q:{)y9N$vB*I6X@

45700水利论文Vj+a1fL@Z(s

37900

g$pfMxyD `3}2t3m0

29200

a%K A#NfHp0

26000水利论文G4U%e8r&e(m(_8Oc/?

9080水利论文O4W ]H&LXs

3150

v:Q?7KM @4T0

H

^"r'K$V],Df0

34.51水利论文8d\!Pms ?

33.36水利论文 H.m;u^:_#N4z,_ Dy

31.95

2X0{;E-TOm WA0

31.20

`k(p A_'@0

26.00水利论文7fRMg,pY/H[%\

21.16水利论文e3W(`:m8W F s


螺山水利论文?7DiA#y4D0[(Z-f[(m

1980~1987水利论文 zyQy3]v kx

31.32

AZM$x%O9`g q6dS$y0

31.27

,H'W*TsD l$ZD0

28.84水利论文+p~9~ulJ tR F

27.90

h0l&_/Y+N0

21.28水利论文(o/Rm^KI&w(A

13.99

C!v%_-c+s0

1954~1959

*q6b$cfx7q0

30.98

8maU|&w9v#`|&@~m0

29.81水利论文4@(H nJYA1ZH-a

28.10水利论文lT%?b MK

27.31水利论文K4bj#i&_

21.89

@5zVx2d,m/Q$s0

15.32水利论文.X1IM K7U|u/a


990310t9.gif (2979 bytes)水利论文6fVm2c^7K3g*o^

图9 监利站80年代与50年代水位变化趋势水利论文[s[Z X/C;`

Tendency of stages at Jianli station in 80s and 50s

1bZPB-l.@0
990310t10.gif (3213 bytes)

图10 0.10频率流量下监利站水位变化趋势与实测值对比

I,jN]POe0

Comparison of computed and observed stages at Jianli station at a flow frequency of 0.10

Y5Zr eulS0

(2)其它频率时,监利水位变化参见图9(b)及图10(a,b)。

wg6j$A7`(bN0

 对于最大流量,实测点为下游水位在80年代为27.0~28.0m,相应监利水位33.5m,50年代下游水位31m左右,监利水位达到34.6m,高出1.1m。

6x]6]2mAO0

  由图10(a,b)可知,将实测点绘在图上,说明计算趋势与实测点相符,鉴于此可以说明计算关系式(10)是可靠的,图中所示的变化趋势也同样可靠。水利论文3w(R[{:`7D#G0w

5 流量水位变化趋势特性分析

'`:t;?,F @Zu0

  在下游水位不变时,流量变化将引起本站水位(或比降)变化,流量变幅与本站水位(河段比降)变幅关系怎样,即流量增加10%,本站水位相应抬高多少,河段比降变化多少,如能知道,将不仅可以对河道特性进行深入认识,而且可以直接而简便地知道流量变化引起本站的水位变化。这个问题解决同样可以利用前述关系式求得。  水利论文8DEQ%Tc b%W

 设下游水位在某数量级条件下,对应的实测流量平均值为990310e1.gif (873 bytes),河段平均比降为990310e2.gif (857 bytes),则

990310e3.gif (1282 bytes)水利论文-_,g%F@6I8Yvh9}

(12)

Y$r!UxcSM0

990310e4.gif (1190 bytes)水利论文2W-n*l}'] Xpd

(13)水利论文$kE(Vm)h;}R

代入关系式(11)中有

Zm{3p a.h&v0

990310e5.gif (1376 bytes)水利论文vV!K,g5j$n'l

(14)水利论文6AQ a I6_G%n

水利论文3Q&c7?yI1TSPU

ηJQ[1-1/J(d+eh)]

:?S [1L M(Qon;Zu/J0

(15)水利论文VkjT\v

其中990310e1.gif (873 bytes)990310e2.gif (857 bytes)分别为平均流量及河段平均比降,且两者同样具有(11)式关系;ΔQ、ΔJ分别为流量及河段比降变幅;ηQ、ηJ分别为流量及河段比降变幅百分数。

S(Y"N^ v g J-m$G0

  该式反映流量变幅百分数与比降变幅百分数关系。将80年代监利站H=23~24m,实测点子与用上式计算结果对比,发现变化关系基本吻合。即可知在此下游水位条件下,流量增加,比降相应增大,当流量增加80%时(平均流量9710m3/s,上游相应平均水位H=27.49m),比降增加26%,流量增加20%,比降增加7%左右。即流量增加80%,从9710m3/s增至17500m3/s时,水位由27.49m增至28.54m(+1.05m);流量增加20%,从9710m3/s增至11700m3/s,水位将由27.49m增至27.78m(+0.28m)。同样可对50年代进行分析,或将其它站进行相应的分析。水利论文 t+KW!wHpFWv G3z5P`

6 几点推荐建议

bA{Cm9V4VR0

  1.本研究所建立的关系式(10)所应用的资料,均建立在百公里左右的上下游站河段范围,当建立更长河段时,如螺山至武汉关,关系一样成立,且计算误差仅为±0.3m,显示了关系式内含规律性的普遍适用性。

a9i/L_,V xg(@0

  2.本研究成文于1998年6月。当1998年7月、8月长江中游出现特大洪水时,由文中80年代所展延的特性关系(见图9(a)),在防总所报98水情数据,均在特性线上,误差均为±0.1m左右。

l_%H6Ml8@ f0水利论文1o%rc)w3[7y/r-C

  3.本研究成果对长江中游水文特性分析,尤其对江湖关系、荆江裁弯河床演变、三峡工程兴建其下游发展预测等重大问题的深入研究,均具有重要意义。

[B+{4e ^o!_{0

:~+\R6H~ g[ AN0 水利论文-HI|7}%ZZl+\

TAG: 长江 水位 戴清 韩其为 毛继新
顶:22 踩:32
【已经有156人表态】
24票
极差
18票
很差
14票
较差
19票
稍差
22票
稍好
21票
较好
18票
很好
20票
极好
下一篇:黑龙江省北部引嫩工程江道演变与治理措施(周晓东 王国志 李亚清)
上一篇:渭河下游近期淤积发展情况的分析研究(唐先海)
查看全部回复【已有0位网友发表了看法】

广告投放

广告投放