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

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

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

2Zj#t7[ U/^2\:m"P0

戴 清 韩其为 毛继新水利论文_DyQ]K Qn

(中国水利水电科学研究院)水利论文S(F"[_+W6V S%j#Jk LA'{

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

关键词长江 水位 单值化

1Q}1CUC&R \&H0

1 问题的提出

Ddf%g"vW0

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

  近来研究水位流量关系特性采用的方法主要有:水利论文|oBgqd5k)Q

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

Zrs{C/wU$f1Tw0

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

6cyneqWd0

(1)水利论文o3I:RL3}X+xR

ΔH=α1ΔH12ΔH23ΔH3水利论文h-d,`.f7U4x F$J\.cK

(2)

pIX9Ng Z,o0

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

zS3`%uO2o;\R0

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

Z=αQβ水利论文6a%t1\|&w+G-Td @*iNr

(3)

!U9fw:AW*B'{b0

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

D)Qi D3p?3B k0J-z0

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

cE/]/J] @j0

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

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

} L T wS*F)}?0

2.1 水位随流量变化特性

1j-i#@'E F0

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

z4y}0IS}J s8RM0

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

G$K?!q.@/G!SM8oP0

990310t1.gif (2402 bytes)

:o[7a}(\.}S"Z"L3c6i0

图1 1980至1987年监利站逐日水位流量关系水利论文 J8j!yJu}_x+p6o3[

Relation between water level and dischargeat Jianli station from 1980 to 1987

(n#Yj/^ T [e0

990310t2.gif (2227 bytes)

!} q;p"e'Y1DCP.{0

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

'~6rp1`LA0

Relation between water level and dischargeat Jianli station (Hluoshan=20~21m)水利论文QW1dijv g^

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

H=ΔHQ+H0

)k6_a3i;p4q"Q;a'Np0

(4)水利论文w7O"A(@{+lW8G,Z!C

2.2 ΔHQ、H0的确定

u1]W f%Xe0

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

990310t3.gif (2734 bytes)水利论文0^V%B&][

图3 监利流量与监利至螺山水面比降关系

?$~,s(\w0

Relation between discharge at Jianli station and water surface slope in reach from Jianli to Luoshan水利论文fu.za h-w

990310t4.gif (2329 bytes)

9O3h4LN6U[0

图4 不同螺山水位监利流量与水面比降关系水利论文/TBO~:B a

Relation between discharge and water surface slope at Jianli station under different water levels at Luoshan station

8Spp;Yv%`;e0

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

J=JQ+J0

/n#qg R9qY0

(5)

]{ m6A%N0

对比(4)式可知,监利水位ΔHQ、H0分别可表示为JQL、J0L+H(L为两站间距,H为下站水位)。即有水利论文Ed M^0[1z"e;Du]e

H=JQL+J0L+H

6M@5~t6C}"RI0

(6)水利论文rI&Y0]c7E&If6S

2.3 经验拟合水利论文em!Ur!] u-X

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

JQ=φ(H)Q

%HI` f8Lz1K0

(7)

|3YMol$n hx3t0

分析数站实测水文资料,研究后发现(H)、J表达式分别为水利论文i\t3kJ9h'i

(H)=10aHb+c水利论文cS%o!q:P7_ ?6CI

(8)

R*us Cl9Gr ^0

J0=d+eH水利论文${u:m-r(jfrk)R.NI

(9)

` G AS\J9}+J0

将(8)、(9)式代入(5)、(6)式中,即为本文建立的水位流量关系基本表达式水利论文-S OtZ;@f4g

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

cJn$hCb2PW^"To0

(10)水利论文-u3wV P+X7SI

水利论文8hd?%Yu hj7l

J=(10aHb+c)Q+(d+eH)水利论文.H@N4g/]f0w%qo VM6o

(11)

vpnA8SgV0VZ+cA0

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

c+GG$r/s*VMX0

3 实例计算水利论文:~z ^6n(M y`V

3.1 系数确定

8t ^nZ h:` V0

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

7v v|5k`C2d0

990310t5.gif (2333 bytes)水利论文;l,Am5CTF

图5 螺山水位与监利J关系

Lk T B&T0

Relation between water level at Luoshan and J0at Jianli station

"y;[GTT'AY1v3t5v0

990310t6.gif (2199 bytes)

3ea*ZlNSW0

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

G/s+x(~`;LRb0

Relation between water level at Luoshan and φ(H) at Jianli station水利论文J[!bP"Y:@|,ErEt

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

表1 系数取值(1980~1987)水利论文@o-M0JCHG:`*d:aK

Values of coefficient(1980~1987)水利论文H5C Xq/cis!H


站名

l*xe7fhO p8uk"?0

监利(螺山)

fz1RUP(Y(qEX s&Y Z0

七里山(螺山)

7Cuoe E] V Fo0

螺山(龙口)

}s L:a^ s0

武汉关(黄石)

+d E5F O t'Z)r4x0

a

"L"`4H%R~S#Z&M0

4.22

@,FizbCl!n0

3.7634水利论文@,p@6[+HS|:H

4.6267水利论文+L,V vQ"m9nBf

-0.85460水利论文8[$M u*T;K7Z(V

a2水利论文3BSa#TgD

-5.333水利论文2u"FQJ1J/daX%R

-5.2259

A|8U5ot5a(f0

-6.1998水利论文Up:w5{|6QgmQ2\

-1.9122水利论文-iG'w;k8u6SzZ

a3水利论文o;P qs:U!m wOm

0.000308

q5Rgl@5dx+C!o0

0.000043304

B-~b:Mm[r(G*q(T0

0.00021897水利论文DYx0D q9`0B#aN*W

-0.000026073水利论文$w'DD,N+n4?*@

a4

N4B-ON9P+A m/I5u0

-0.00061278

4S4vGJkC-U0

-0.00071889

-Z)UF U/W vn0

-0.00062114水利论文0yn} yR

0.00055556水利论文)cM"F}A c:H$Y

a5

tj5O~ _0

0.031678

1IN$H:hTj thW0

0.049650水利论文@5T6X'N(Qp

0.028696水利论文#]"u/uf$h)|f,{jv

0.0006667水利论文2m/vMm(V*fo!PB


3.2 验算结果水利论文DxvK2B*f7Q-F

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

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

7B6V#QMekg0

990310t7.gif (6106 bytes)水利论文E)|,`fkKz nx!X

图7 实测水位与计算值对比水利论文S4G(DY']scC.r

Comparison of observed and computed stages水利论文:s ST _;@:p

990310t8.gif (2345 bytes)

k*]*A6C8Ct0

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

F` t;`#l0

Comparison of observed and computeddischarge at Jianli station

tH5|}Lb'Qa0

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

4 监利水位年代变化分析水利论文XEYH5_@c

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

4.1 水文频率分析

&J%J&i.h[6O&z0

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

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

e*?.r \Z-d+r4R0

4.2 监利流量同频率条件水位变化水利论文 M? AO!Vw2Y6H4Z

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

'm h9dFA4RuY0

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

pg/\8i:M"~0

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

7wJ(o.lv0

Analyses of frequency of discharge and water level at Jianli and the water level at Luoshan水利论文 yJU V b0O


@TRAW0

年代水利论文1{_ |!Y eX

频率

W St*DDb0

Max

tr'KYl!@0

1%

;m_N~G0

5%

A)H KdJ6C8F&W cD`0

10%

y)arzQ&c$C0

50%水利论文?q8K HG gb

min水利论文o`#ks;d1s0u'q)N!m


监利水利论文h&xA`1o"C J

1980~1987

8x7N0F \#w }'[0

Q水利论文(nU_&P$@YA[

35200水利论文n y3xUb3ic

29400

9Z};L IY0

23200水利论文.W:g6j V%A2Y-q

18800水利论文Mj^6H)[F^d L.SG

7810水利论文'X'b.c/{0B6g+c%X

3140水利论文(IY KR `u

H

SE/QV4@0

34.7水利论文F5L/Os&{$R`CE

33.89

$Q1u8b)hl{Lyn0

32.31

*Tr#w[\l0

31.26

}"iB!IZ)_}@0

25.87

Fys.rf }"Z GK2}/P0

22.08

W9H$F3SM]"U`0

1954~1959水利论文5E1J w+J{ N%f1u p

Q

9_C? N3a}:H I5]*D0

45700水利论文x/|qzkQK

37900水利论文Q~V&A&S&W9p&k H

29200水利论文W5N6^.Qh*y H

26000水利论文 WER]*}.?[

9080水利论文.}I? k;t9z1]6G)ju

3150

*YMc}Dg6\"xv+z0

H水利论文#w Ca2L?s9If

34.51水利论文.c#Mnn-U2x

33.36

NAa kGH&` T0

31.95

p/qk N+QyE0

31.20

IQ6r CA~c6O'Q.\M0

26.00

IC$\fLi0

21.16

X/pX%DM/?&QOC0

螺山

.@)Fs"c\H1fNE0

1980~1987

X5E0vd){w2ae0

31.32

:Q&|5r:X(_o$PF8Y0

31.27

,e0E+C%m3XZmz0

28.84水利论文$K m%Z Pt B8xz

27.90

?&X.zvjV p3m*Ro}0

21.28

Xy+d?t6i#fS u?l0

13.99水利论文2v4^-C/V,j(H.l:d7V

1954~1959水利论文+_X2?)hL~?

30.98

3O6}i+^:|9cBBZX0

29.81

W'`.NVV{#Js0

28.10

2J1z[J4s3c0

27.31水利论文-e-?4Yy}A3z4K

21.89

5OU"y$T;d yp0

15.32

[)W$L#~ o9g[w,x0

990310t9.gif (2979 bytes)水利论文0^ue4Z&@$c0A

图9 监利站80年代与50年代水位变化趋势水利论文8W6B3@ob w3`O.?

Tendency of stages at Jianli station in 80s and 50s水利论文;wC5^!ij[,D

990310t10.gif (3213 bytes)

图10 0.10频率流量下监利站水位变化趋势与实测值对比水利论文2u.u8R;c/o#|

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

e1tzr[X zab%[0

(2)其它频率时,监利水位变化参见图9(b)及图10(a,b)。水利论文\4NSw'O0@k9T

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

2Oa W qIFfY0

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

5 流量水位变化趋势特性分析水利论文 c&s9]2qDEv6A

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

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

990310e3.gif (1282 bytes)水利论文 s8Y9n}S],XH#QJa

(12)

M"FGh%V-\0

990310e4.gif (1190 bytes)

(gL4? O [Pbf6t a0

(13)水利论文 ]!m$N9Z0v_)U5V

代入关系式(11)中有水利论文E"[o_qzt

990310e5.gif (1376 bytes)水利论文5cjvj[)B

(14)

1\6RuhP hq#Om0

.yT'v1CN ]S0

ηJQ[1-1/J(d+eh)]水利论文Gm0H3R!?-F

(15)水利论文|4w$Xr,Ek

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

u7US Zd7KQ/R*B0

  该式反映流量变幅百分数与比降变幅百分数关系。将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年代进行分析,或将其它站进行相应的分析。水利论文X @.GZI@SzOy

6 几点推荐建议

nyI"n/Q0

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

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

水利论文 N;R,qD%[3F)w

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

4S2_c5?-\F;f|4\0水利论文 b$S-wvlp

 水利论文0M i+Ic?

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

广告投放

广告投放