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

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

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

8W;c JR5enk0

戴 清 韩其为 毛继新

,{_lgAQ2`Z[0

(中国水利水电科学研究院)水利论文o9DG(QM7V_]+f;E

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

c!V"NWH0

关键词长江 水位 单值化

a#? r3p Z f0

1 问题的提出

{ BoKe0

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

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

M^N I MG:P.u4\J8~0

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

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

dz)Nz"I ^t.]0

(1)水利论文"y h"f|$y&y

ΔH=α1ΔH12ΔH23ΔH3水利论文-u!W LR/Z;i E ^(Ca

(2)水利论文-d9E^q;Ed

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

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

*n+I8Q:vS^n @3_A0

Z=αQβ水利论文 A-pr9DX`5A

(3)

L9R G/U3b0zv0

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

i f4d+W*O'F0

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

'](IE~$}SO"E0?0

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

;O;\'D6n#eHg0

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

;Ug*_K/QloV0

2.1 水位随流量变化特性水利论文;uK$v/l"i

  点绘监利站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。水利论文U,m R x ~'`?j

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

990310t1.gif (2402 bytes)

?ko Le5Rl#x0

图1 1980至1987年监利站逐日水位流量关系水利论文eE&|u.{` v

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

7J5ZV DN7_v0

990310t2.gif (2227 bytes)水利论文r%D5n.F!N0Ajek,A

图2 监利站水位流量关系(H螺山=20~21m)水利论文\ Qvw ^ [Z[z

Relation between water level and dischargeat Jianli station (Hluoshan=20~21m)水利论文8`\7h9sH'f/q

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

Y6K @3@KY0

H=ΔHQ+H0水利论文&},t FZPRLEy

(4)水利论文KUB9U!qXd

2.2 ΔHQ、H0的确定水利论文7L#S:XF*D.E/Xo

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

990310t3.gif (2734 bytes)

L_,F%meJ0

图3 监利流量与监利至螺山水面比降关系水利论文Uck]9y~

Relation between discharge at Jianli station and water surface slope in reach from Jianli to Luoshan

c5y?S5`0|7h P0

990310t4.gif (2329 bytes)

"y|V:n ]0

图4 不同螺山水位监利流量与水面比降关系水利论文[2f@w] d%FL

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

+wE/K2}@}^#C0

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

R:e*cz,I qk0

J=JQ+J0

!psQ'a hn0

(5)水利论文b;s4orL9Vl

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

H=JQL+J0L+H

b0G%v0|:`'a0

(6)水利论文DtcS2o%_3t e's

2.3 经验拟合水利论文*o)f5od)ez:NMSj

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

ZX2i q$[pW0

JQ=φ(H)Q水利论文Un7F*Q1X3G0J/}`

(7)水利论文Z+l ^6`ie[c

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

mdD9j0}5L O$c0

(H)=10aHb+c水利论文Z/Stp;qqk,T

(8)水利论文MwUrTb

J0=d+eH

h s7B7p5n _4m"vl Qg0

(9)水利论文 }k8\^+q6J

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

L)O e"{NLo0

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

AJ%Y?'}An"[0

(10)水利论文XdOX)Z+i-t

水利论文"w1m`Ll t

J=(10aHb+c)Q+(d+eH)水利论文] C1i:vYdi#Ki1}

(11)

}7sLg0J-bc1z0

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

$yW{&@ [ lt/w0

3 实例计算

0W#A9c2g}%R!i N1T0

3.1 系数确定水利论文ZRk%^k(_

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

990310t5.gif (2333 bytes)水利论文5TbcF1y j}b?

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

!t}8b1GWK R Ee0

Relation between water level at Luoshan and J0at Jianli station水利论文:vjB5z1qf

990310t6.gif (2199 bytes)

joh(R2Y0

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

B]u1^&mU b%veM!SR0

Relation between water level at Luoshan and φ(H) at Jianli station水利论文 {Ql MzW0A

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

p$S Us4YK){$W-_0

表1 系数取值(1980~1987)

/w,T/i/@ITfs(o0

Values of coefficient(1980~1987)

GcO+p^eh0

站名

~ydy"JN$JEh0

监利(螺山)

dj'Oo*t5C/G0

七里山(螺山)

DXR"^]TrA0

螺山(龙口)

-K _?EqL5a"Ci0

武汉关(黄石)

L6p X [-~ h0

a

D_upo0

4.22

OIY a l I)j%zS0

3.7634

hL;X^*I LX0[0

4.6267

noo7[ _f7{+~ v$g0

-0.85460水利论文r M,g8L6p I V_:O~

a2

'cEBg6QB]0

-5.333水利论文]2ih6E Mo5[

-5.2259水利论文6SHkj r

-6.1998

0w$UBZ-\8f/Q2^'I-d0

-1.9122水利论文w Z}HNnQa0c

a3

4l!S3h4gsAaP0

0.000308

1E*R:g NE/|2Hjfn0

0.000043304水利论文0boF*P7EH

0.00021897水利论文\l\X$L Pw!z,J

-0.000026073水利论文N!P2eep`7Es}

a4水利论文:F8r4t4xPdtx

-0.00061278水利论文C A)E/JRI@

-0.00071889

!^#J+GxOKE0

-0.00062114

D"V/H/E S2bye$c0

0.00055556水利论文T2}%^9N4tK9L

a5

*]t{9Rf3S|:Y.t0

0.031678

ke7N7U'M7v0

0.049650水利论文c^3O.X:{

0.028696

.G:^3E+V$\"E;k8v0

0.0006667水利论文3N"e q6F&q'z^u


3.2 验算结果水利论文 |r jJ2[)s{

将统计年份的点计算对比结果点绘于图7(a、b、c、d)中,可见,计算与实测点完全吻合。

U3K I wwnFc:m0

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

_0I;E.m yhN0

990310t7.gif (6106 bytes)水利论文X EB-AV2L

图7 实测水位与计算值对比水利论文:?]'P yg/s

Comparison of observed and computed stages

:aK1G#Y {0

990310t8.gif (2345 bytes)

!S t bp |0

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

[M]-ZfJPa"D$J0

Comparison of observed and computeddischarge at Jianli station

$`Fd#r!t9_:rG}0

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

4 监利水位年代变化分析水利论文L#r D&qO cu

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

4.1 水文频率分析水利论文nv3Xb? Q;W3P/K

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

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

9[ Pb%l ^?0

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

(e$U eX!|{#w0

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

I.H3L!D;L_&S n0

(1)不同年代出现最大流量时,下游同水位条件下80年代水位具有高于50年代水位的趋势,但实际上因下游螺山水位在80年代低于50年代水位,因此监利水位在80年代低于50年代。水利论文}c{9N"x R;h

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

7P\6b2i4a_0

Analyses of frequency of discharge and water level at Jianli and the water level at Luoshan水利论文n:HFn@d


水利论文&Q;d(D"m!|Mq

年代

%j xglUp-d0

频率水利论文k Z;RW?~$l

Max水利论文D\2O']:h(l7O

1%

sK!~({&v8|0

5%水利论文Z8W5g}L&@6T

10%

9}!f:A(NW,Q E0

50%

k9?fRou-RY8sR0

min水利论文M!HAt}C0_*Qh


监利水利论文6c;s1E$qh P.V

1980~1987水利论文8h @0[T&pH1f

Q

8WL5H:Vx!z0

35200

Y _*LNS9bRj-PmT0

29400

b?3Qprd0

23200水利论文1a i)k+p9X%K ZWu

18800

j)uR y e5bYyj0

7810水利论文] ^)S-YTL

3140

-O@(y2FC;f3Z0

H

:JCe{XY*YD+m0

34.7

/ZcG+R(@ H0

33.89水利论文iN)v8W#Ir i

32.31

kiP'YoGH8Z;K b0

31.26水利论文(S!y8iz }0L5q5~

25.87

:iK._BN0

22.08水利论文A!}/h/jw

1954~1959

.S"OA*~ @a ^r9N0

Q水利论文#G!`qDh0z}

45700

p,bG'X @$i0

37900水利论文[ L;V"Oq0{B K

29200水利论文5|6~*MD7U-i5aP%c

26000水利论文.M)Qd!a1OM)\V3s-z zY-q/Y"`

9080

4G%Zj,u S0z0

3150

WX[)HA4mji8@0

H水利论文hb `zO&Pct

34.51水利论文Un-X6mKN7m1X4w#c

33.36

;FepB@v)b0

31.95水利论文"L,^/Y/|a M

31.20

sZf'?4V0F^0

26.00

~u9H1]8bZ PO}3s \0

21.16水利论文j C%B{5^


螺山

uP)I*WE8S V+n0

1980~1987

)Uj'Y!k*fo/A[ G9Yc0

31.32水利论文)Ka_q*k4r

31.27水利论文m%Cycbb0y

28.84水利论文O%K*l$@^

27.90

Q:yMJ3]V0

21.28

#gH&C'P4N#k0

13.99

-?ws0K*U1B{5F0

1954~1959水利论文}/Vr({iK"}

30.98

N&lNW7U7k`0D0

29.81

Q;o$E7}:GRZ9Ym0

28.10

_ nq#^Ho2N8Y0

27.31水利论文 y&h!Kgw3wF/H

21.89

5u\K)wT0

15.32水利论文7g(xp FK wV x3Is


990310t9.gif (2979 bytes)水利论文hF2c0y%t@

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

M V6e2k"qi1z0

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

G Q}4AE[G`}0
990310t10.gif (3213 bytes)

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

/_&zy4u!^K#Vc0

Comparison of computed and observed stages at Jianli station at a flow frequency of 0.10水利论文THbO$zIW(S1g

(2)其它频率时,监利水位变化参见图9(b)及图10(a,b)。水利论文$m4Y+?"pX?

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

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

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

+Y.bM,I:r-{0

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

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

990310e3.gif (1282 bytes)水利论文?4~O6z SW

(12)

1P3O N|(F0

990310e4.gif (1190 bytes)水利论文Z.p*O#i1A

(13)

5?!c^Yp3V"I3A |'b1?0

代入关系式(11)中有

B7C'YyO(P`0

990310e5.gif (1376 bytes)

Q)c QU0\.U([(f0

(14)

UIC?7SX0

水利论文 S4@(|l!v2Z&bD

ηJQ[1-1/J(d+eh)]水利论文\b vts.fv3T

(15)水利论文MaN"]o]

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

TE7dIE5c:^r0

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

6f'Q[h8h0

6 几点推荐建议水利论文q}ez7V"Y0x0i9lnI(N

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

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

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

水利论文 _:l L ~$rK |8p

 水利论文1`k?:j.iI&i'T

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

广告投放

广告投放