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

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

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

戴 清 韩其为 毛继新水利论文V*bf&^k8Q'~

(中国水利水电科学研究院)水利论文 H]+}Z1M#{ YyZ

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

]9y+E%Yyq8]N0

关键词长江 水位 单值化

5@ X;Z9f\e$D!Z$~0{;|0

1 问题的提出

MKX8@B `$aK5w0

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

  近来研究水位流量关系特性采用的方法主要有:水利论文2aQi&V_#i

  (1)综合落差指数法:对城陵矶七里山站流量进行单值化处理,公式形式为

Au(EK6B!O0

Q=Qm/ΔHβ=f(H)水利论文]T}3y t2b6A0[ {

(1)

4Hmi&`QD w)Z0

ΔH=α1ΔH12ΔH23ΔH3水利论文J&I@&O^z

(2)水利论文,F;ep]){X

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

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

Z=αQβ水利论文*Xkq-o lK

(3)

J#O0OL&l#a0

式中Z为水位,Q为流量,α、β为待定系数,不同流量,α、β值也不同。水利论文L$L#_6w'_

  研究者认为,上述方法一般相关关系较好,除个别点外,相关系数可达到0.9以上。

SP%]-Nm qC&[Q7AL0

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

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

e \,_cff3ovsT)v0

2.1 水位随流量变化特性水利论文 h"@N)E|`

  点绘监利站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。

&}%Bb+N,R4uWt1B0x0

如将下游螺山站的水位按大小划分若干组,则可知每一组的监利Q、H均存在着图2所示特性。水利论文N)jZYDo0t

990310t1.gif (2402 bytes)

.s#sF4I$B0

图1 1980至1987年监利站逐日水位流量关系水利论文&X$h%w)Kw

Relation between water level and dischargeat Jianli station from 1980 to 1987水利论文&V X2rQsu

990310t2.gif (2227 bytes)水利论文1jCJ;r2eU

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

e QK#U)bwA0

Relation between water level and dischargeat Jianli station (Hluoshan=20~21m)

g/p'P/lRc f x^ d0

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

kR-b[.O|8]U0

H=ΔHQ+H0水利论文Lq_9MIp

(4)水利论文G+A$\\'@A o(\

2.2 ΔHQ、H0的确定

6} [R$mMb0O0

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

990310t3.gif (2734 bytes)水利论文 v?F2cr

图3 监利流量与监利至螺山水面比降关系水利论文-Q,K&u|-GR

Relation between discharge at Jianli station and water surface slope in reach from Jianli to Luoshan水利论文!j$^j ?*?,_ Kx

990310t4.gif (2329 bytes)

M+d&J'fI'h0

图4 不同螺山水位监利流量与水面比降关系水利论文*~Bf Iiu8b!U

Relation between discharge and water surface slope at Jianli station under different water levels at Luoshan station水利论文lv&tx#rhd#D:EI

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

J=JQ+J0

(n5xB5R.fMO0

(5)

!|3rc!m q;cp0

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

Z ?~;g/_Z?`+_^.}0

H=JQL+J0L+H

5Q:fU6I ]0

(6)

)fy0{(P5z0

2.3 经验拟合水利论文t@P^|

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

)^\9l1x%@R0

JQ=φ(H)Q

#mesx7ff'S z0

(7)水利论文C_-?:R(SlU

分析数站实测水文资料,研究后发现(H)、J表达式分别为水利论文&t3tlg~C

(H)=10aHb+c

:E d7pI-mU)yI6A,v0

(8)

yg}g KC\NK0

J0=d+eH水利论文KL+G/bB_-i5L6Q

(9)水利论文0Tjb R'b `\ ~p

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

3K!Q:i`d0f8f0

H=(10aHb+c)QL+(d+eH)L+H

B5F&H7Ba@g0

(10)水利论文 AR4DY V"J

M0R~ m'P9C7S%Oq0

J=(10aHb+c)Q+(d+eH)水利论文 ?2W:j.R!M$}g

(11)水利论文5V:p6_!aRbo

其中 a、b、c、d、e五个为待定系数,且a、b、c与d、e互不相关,每个系数均有明确的数学含义和其简单的计算方法,后文中将详细介绍。其他符号意义同前。水利论文1[8s"{4g#\s{I)G3H#R Q

3 实例计算

#h:[ GP8RYll0

3.1 系数确定

S9DO+f[ XB0

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

~7F3P E2w#f9\+@0

990310t5.gif (2333 bytes)

/ruRy{:gpN0

图5 螺山水位与监利J关系水利论文P,P(_EDp*D-u

Relation between water level at Luoshan and J0at Jianli station水利论文 F#|q` \^ U2R

990310t6.gif (2199 bytes)

{9s|D#Ja(P"O0

图6 螺山水位与监利φ(H)关系水利论文|pU]7Wx W}yhS

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

T/~:k7r%u@Wb0

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

表1 系数取值(1980~1987)

+Q8N:?$D?)V0

Values of coefficient(1980~1987)

/d+D7B}Z+A A'G4D6q(d0

站名水利论文 S$c f,E+w9J;}-I+?

监利(螺山)

;`4UwXUw$f0

七里山(螺山)水利论文C$eRy&v#GA/? J

螺山(龙口)

.ScX:mb&Fh v y6a0

武汉关(黄石)水利论文^)uJop*y A


a水利论文z'v&|h:u"X

4.22水利论文(tN7~K p]Z_W

3.7634

Hus!q&~&e;R,zdZ HQ0

4.6267水利论文;f Zxor'F0rOl

-0.85460水利论文0G4i_J-RrXmo

a2水利论文WX&r%P5p5z'G)A2r

-5.333

d6?*d(@)m1ks0

-5.2259

"Y['g0X8FS_0

-6.1998水利论文%A!b m-DCm)h1B

-1.9122

&lV[.UT0

a3

^Td#yy(|#Q0

0.000308水利论文,VK*}9U2T

0.000043304水利论文$Q7U@0H!s C$Q

0.00021897

#U*JO9E2W{l6e'W0

-0.000026073

(}*^!h2`:l0

a4

]J GwQ He5bz0

-0.00061278水利论文 J7Zmu#W

-0.00071889

5C|$vt1D| w0Pu Q0

-0.00062114

8U@Q*O(|A0

0.00055556

'sr6GMu&I J0

a5水利论文8ZA2E ])S-QL,D

0.031678水利论文gYa(R4h

0.049650

eK)XGL7KXo4Q:|9P0

0.028696水利论文Wk k&kW

0.0006667

)Fw'C.iL0

3.2 验算结果

)nD-c1jb"}cs0

将统计年份的点计算对比结果点绘于图7(a、b、c、d)中,可见,计算与实测点完全吻合。水利论文y(`q.A&W ]:g xWe

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

*M ^k%CT Sd0

990310t7.gif (6106 bytes)

4ICw/}nb FH0

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

6P0t!j w&GQ0

Comparison of observed and computed stages

$W!F+tfnAc%f0

990310t8.gif (2345 bytes)

Tdq/f3g#z"FI0

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

Hw"L FU.TM2f0

Comparison of observed and computeddischarge at Jianli station水利论文$`WM%U$t!Y

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

4 监利水位年代变化分析水利论文B-q3EhdN)Bd

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

4.1 水文频率分析水利论文?3l0j6p!?B;D/e

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

h+A tZ+e5XA4M+m X0

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

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

.p `f{2H0

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

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

9TMd5N+IVb0

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

7Ki/@ oRg!d2M J0

Analyses of frequency of discharge and water level at Jianli and the water level at Luoshan

'j3z+T&k4M[!bp0

5{b^5RT/KC0

年代水利论文)Y*iMU*]x(z`W\d

频率

x!n1`Ih0

Max

4m5M nq2y/I GU0

1%水利论文 K%|4\(`7gK0]y;TO

5%水利论文$U+M7_q*U&G

10%水利论文5XOcs%\hX4x

50%

6E3?2TG0N0

min

.jk+l]/L*DV lA/r0

监利水利论文@/aSyhPec!y,G+?:C

1980~1987

0H6Y#`:` ` ?/x0

Q水利论文XtG?m3i

35200水利论文8yKA#fj2kbw

29400水利论文$I)h:gy7Avo:Z

23200

0upG*yGJ0

18800水利论文u1a)~*[j

7810

|,Zy8C_ W0

3140

:WWAsBI3hs6I0

H水利论文 `9R,qJ!xGR

34.7水利论文L&I!O%S#zf3r8LK

33.89水利论文4Ra OZ;Z5\

32.31水利论文-~.h:PGR

31.26水利论文3A:i$])e(TW

25.87

I8XO i7Dm'q4ok.v0

22.08水利论文/PX x t Vg8R@Z1tJC

1954~1959

:C|$K_u H4Y8K0

Q

6lAB0a ~%R5d h7B0

45700水利论文6C Q)@%Q'M[5N%t7o

37900水利论文"F,TFM8R'V c

29200

x6Cx[H;p0

26000水利论文HHE4^uLAj3C|

9080

!G9F1x u1@5X9k0

3150

`[/d1e(as3j/T0

H水利论文QM Ho mYb5x

34.51水利论文&J"_ y3Wq

33.36

PJ1P9w:@g0

31.95

ew'AG4L0

31.20

uqqL2NDb0

26.00

hpe}ap X0

21.16

q#O/?{#v3];R:D2co w;xI0

螺山水利论文iq f O g)O-kcO^0^F E

1980~1987水利论文9SphE3?tg

31.32水利论文 C#J;HPix

31.27水利论文 rE-@fuBDJZt

28.84

X8D'A;w K*{R.c0

27.90水利论文{-m$u"?[\

21.28水利论文3i3].K#A?6JZ?4S

13.99

/U.^,lUl0

1954~1959水利论文%U*S/H7|6S@(FC Q\b

30.98水利论文2t,D#ag/nx

29.81

*KK#~(~#vdT0

28.10

5KgBL0W&L0

27.31水利论文g:\8uMk%_6gV

21.89

L9lz1G4D!~0

15.32水利论文2R7} B:W2Y\R.A


990310t9.gif (2979 bytes)水利论文&bt OuL:F5z[

图9 监利站80年代与50年代水位变化趋势

k)ZF_$z,w C0

Tendency of stages at Jianli station in 80s and 50s水利论文g e.l T4_

990310t10.gif (3213 bytes)

图10 0.10频率流量下监利站水位变化趋势与实测值对比水利论文2^M z%^(n

Comparison of computed and observed stages at Jianli station at a flow frequency of 0.10水利论文+VhCg"Uj,v8C4I3` WF

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

a/{O n#|[-q$D0

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

  由图10(a,b)可知,将实测点绘在图上,说明计算趋势与实测点相符,鉴于此可以说明计算关系式(10)是可靠的,图中所示的变化趋势也同样可靠。

J'iL3hk|g~9x_J0

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

7WI$O"BG:t1c}+r+E0

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

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

990310e3.gif (1282 bytes)水利论文ak,^u*C%N

(12)水利论文;g9]\#Q8z6ZK

990310e4.gif (1190 bytes)

&Q C$ylIU2u"Up0

(13)水利论文/Bf@Umu7^

代入关系式(11)中有

S7r'{o D @q;Us!?(|0

990310e5.gif (1376 bytes)

/};e&`*d _~uG0

(14)水利论文'V*E z%oXTA,o

水利论文1l2j:?o,\,[Lq$[O

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

M|;Ks4rBV}e^`0

(15)

n h'p ca'B|h$?0

其中990310e1.gif (873 bytes)990310e2.gif (857 bytes)分别为平均流量及河段平均比降,且两者同样具有(11)式关系;ΔQ、ΔJ分别为流量及河段比降变幅;ηQ、ηJ分别为流量及河段比降变幅百分数。水利论文/U!_;v'iI b:Mt

  该式反映流量变幅百分数与比降变幅百分数关系。将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年代进行分析,或将其它站进行相应的分析。水利论文UL'G'a7m M9R4L {:\:l

6 几点推荐建议

z;^;k1D2{ v;v&G At0

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

r yg"z!R^%g~0

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

bTcC;p!j r0水利论文+q0X#O%dfy

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

3C p2e1T8u-q {4[n0 

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

广告投放

广告投放