你的位置:水利新闻 >> 资讯 >> 水利科技资讯 >> 详细内容 在线投稿

引黄工程水源地年径流量长期演变规律研究

排行榜 收藏 打印 发给朋友 举报 来源: 水资源与水工程学报 发布者:cjk3d
热度221票 浏览97次 时间:2015年5月14日 12:56

hJ7c$^ ],b0作者:徐聪 冯民权水利新闻dEf:K,`@ x

^'C Y;W g0摘要: 利用Kendall 秩次相关检验法、R/S 分析法、游程分析法和马尔可夫过程分析法对中部引黄工程水源地年径流量长期演变规律进行分析。黄河河曲站年径流量时序在1978-2010 年整个时段内呈现出显著的减少趋势,并且可以预测出未来河曲站的年径流量将继续呈较显著减少趋势; 河曲站连丰、连平和连枯状态的出现概率都随连续年数的增加而减小; 在长期丰枯演化过程中,偏枯年的重现时间最短,偏丰年的重现时间最长。中部引黄工程水源地年径流量未来呈较显著的减少趋势; 在年径流长期丰枯演化过程中,出现偏枯年的可能性较大。本研究旨在为黄河水资源的可持续开发利用和中部引黄工程供需调节提供科学依据。水利新闻f\0_#@M'b

_W4hQ^2c7W+E6|0关键词: 年径流量; 时序变化; 丰枯演化; 中部引黄工程水利新闻u J{8D/Q@Q

水利新闻@pB%y-d

中图分类号: P333.1 文献标识码: A 文章编号: 1672-643X(2014) 06-0105-04水利新闻E!fBvZZ(V-K

:`(h!W8r%B!{:T+L0Study on long-term evolution law of annual runoff in water source of Yellow river diversion project水利新闻F| }N4rW O!~/w

*WbHn,J3Pe0XU Cong,FENG Minquan

o F'i#PA O L0

*|(s1D5XEb#_0Abstract: By use of Kendall trend test method,R/S analysis method,run-length analysis method and Markov process analysis method,the paper analyzed the long-term evolution law of annual runoff at the water source of central Yellow River diversion project.The annual runoff timing at Hequ station of the Yellow River showed a significant decreasing trend during the whole period 1978 to 2010.According to the result,the runoff will continue to show a more significant decreasing trend in the future.The occurrence probability of continuous multi-annual high flow state,normal flow state and low flow state decreased with the increase of consecutive years.In the long-term evolution process of high and low state,the reproduce time of partial low flow year is the shortest and the reproduce time of partial high flow year is the longest.The annual runoff at water source of the central Yellow River diversion project will show a more significant decreasing trend in the future; In the process of high and low evolution of annual runoff,it is more likely to appear low flow year.The study can provide a scientific basis for the sustainable development and utilization of water resources on the Yellow River and the regulation of supply and demand of the central Yellow River diversion project.水利新闻e*G[){;F uUT

水利新闻Sy(OD4?`Q

Key words: annual runoff; temporal variation; evolution of high and low state; central Yellow River diversion project

EJNj~5~v0

^uc(|6b c(YI4Y+Y0V0黄河是中国的第二大河,干流长5 464 km,流域总面积达79.5 万km2。随着经济发展和人们生活水平的提高,黄河水资源供求形势日趋紧张,大型引黄调水工程也越来越多。为此,近年来许多专家学者作了大量的工作,从不同的角度来分析其长期演变规律,并且取得了大量有意义的成果。如ShiChangxing 等[1]基于长期水文气象系列、数字高程模型、土壤数据库和黄河中游水土保持措施的文件等多个方面分析了黄河中游年径流量的变化,并研究了其与水土保持和气候变化的关系。Wang Suiji等[2]利用累积异常识别法、线性回归法和改进的累积量斜率变化率分析方法(SCRCQ) 分析了气候和人类活动对黄河径流变化的贡献。还有学者分别利用希尔伯特-黄变换(HHT) 法[3,5]、重标度极差(R/S) 法[3]、Mann-Kendall 秩次相关检验[4]、流量历时曲线法[4]、双累积曲线法[4]、小波分析方法[6-7]的一种或几种对黄河径流长期变化趋势进行了分析研究。但是,这些研究所采用的研究方法相对较为单一,其研究结果具有一定的局限性,不能全面的揭示黄河年径流量长期演变的时序变化特征及丰枯演化规律。

%[9M!Z3F/d0

#tKP8]C4H _0为了能够全面深入分析中部引黄工程水源地年径流量的长期演化规律,本研究结合Kendall 秩次相关检验法、R/S 分析法、游程分析法和马尔可夫过程分析法这4 种方法的优点,并将这些方法进行了有机的结合。首先应用Kendall 秩次相关检验法对径流量的时序变化趋势进行分析,并在其基础上,应用R/S 分析法分析径流序列的持续性; 然后,应用游程分析法和马尔可夫过程分析法对径流的丰枯演化规律进行分析。目的在于为黄河水资源的可持续开发利用和中部引黄工程供需调节提供科学依据。水利新闻rz Q#e WL

水利新闻C\2k^sk |b

1 基本理论及方法

!]PgN2^g`0

z4qK Y G8E}#{01.1 Kendall 秩次相关检验法[8]水利新闻N lu2o8uE

hl.P6uCa NR0对于序列x1,x2,x3,…xn,先确定所有对偶值(xi,xj,j> i) 中的(xi < xj) 出现的次数,记为P,根据计算的P,代入下列公式:

1\swy$t^0 水利新闻F%gDb%i)D0i

,u%qX5KUS;V8k0 水利新闻,OJq+]J

式中: τ为Kendall 秩次统计量; U为Kendall 秩次相关检验系数,随n 的增加很快收敛于标准状态分布;n 为系列长度。

o:}D1DENK em[6q0 水利新闻OmC)oo0QAF&wK|

若给定显著性水平为α,则其临界值为Uα/2。若| U | > Uα/2,则说明序列具有明显趋势性; 当U > 0时,则具有上升趋势; 当U < 0 时,则具有下降趋势。

t5B6K6l+V"{0 水利新闻7NE {&h,XG

若| U | < Uα/2,则序列趋势性不明显。

Wh4tS9~+|0UA(?Q+ZwG0 水利新闻zK?;no7V6qr

1.2 R/S 分析法[3]

3c&X(ec\4PA5q0

A4`&X+d}CT \`H0将R/S 分析法应用于径流时序特性分析时,判断径流时序状态的持续性,可通过计算Hurst 指数。

(v8fb$@8Krg(~+z0

"M+~8B2h$k0Gc qg'e3j~0R/S 分析法基本原理如下:水利新闻W8x z(WFrCN\

Mm7mJ8u0对于时间序列{ x(t) } (t=1,2…,n),下列量可定义为:水利新闻C:shs6?0eG

-\ [,n$GA0均值序列:

K Y1|%j8dV0

^ua_ q!_0

vBL-_%yY0 水利新闻:IwG s.TK3z

累计离差:水利新闻SP(sG,P3b zU W

@!O(f3v2k0

"@ J KLJEm!b h,U0 水利新闻Dag VFP\%wb }I:U

极差序列:

)l(])^tA/S}n-N0

XR c} e-j0

'EYR2OJ*J+Nk.aw0 水利新闻Q/v|(v%l

标准差序列:水利新闻M-kqv5U*p;Yr

|'N I\!qL j:c/S0

5_e[*U7VX0

(x G Ao)gj6EkG0如果比值R(τ)/S(τ) ≡ R/S 存在关系: R/S ∝τH,则说明时间序列{ x(t) } (t=1,2…,n) 具有Hurst 现象,H 称为Hurst 指数。在双对数坐标系(lg(τ),lg(R/S) ) 中用最小二乘法拟合,拟合直线的斜率即为H。水利新闻yYWR%T%MlC

水利新闻J|6?'m1^l4AA U;{

根据H 的大小,可以对时间序列趋势成分所表现的持续性进行判断。Hurst 指数取值范围为[0,1],其存在以下3 种情况: H=0.5 时,表明该时间序列变化是随机的。0.5 < H < 1 时,表明该时间序列具有长程相关性,即过程具有正持续性。0 < H< 0.5 时,表明该时间序列具有反持续性。

4o s R p C4mC[0

7L0o.\0k#e01.3 游程分析法[9]水利新闻's-{]Q!u*r

)I[2E-kQJ#Vf0游程分析法的计算可用以下公式:

$}3s:]H8AE1fM @0

.W]k:l8n(k"n0P=ρ(k-1) (1-ρ) (8)水利新闻]_vm t

-Eo1JX ]0ρ=(S1-S2)/S (9)

2k MS6B2FN/a0

d*qx`YFL+n'F0S&~L0式中: P 为连续k 年枯水(丰水) 的概率; k 为连枯(连丰) 水年的年数; ρ为分布参数,其值在0 ~ 1 之间; S1为时间序列中连枯年(连丰年) 总年数; S2为时间序列中各种统计长度的连枯(连丰) 年发生的累积频次; S 为时间序列枯水年(丰水年) 的累积频次。

*W7pz H;l![ _h?'E9Y0 水利新闻cCq.dD|e@@3|

1.4 马尔可夫过程基本理论[10-11]水利新闻*A$U}ye-`0]/`C~*j

水利新闻0lw$s?3N`"GD

马尔可夫分析法是指在马尔可夫过程的假设前提下,通过分析随机变量的现时变化情况来预测这些变量未来变化情况的一种预测方法,又称为马尔可夫转移矩阵法。马尔可夫一步转移概率矩阵,用P(1) 表示; 多步转移概率矩阵,用P(m) 表示,它是由初始状态经多步(m 步) 转移从而达到相应状态的概率矩阵,则P(1) 、P(m) 可分别表示为:

&k6h'f.MQ Iazx0 水利新闻Eb5?4Lw%q/g/Sw

水利新闻Wq XN@ RoH

/i[$@,b[0水利新闻6Fees ~ys*g

!|2_sV;Yomv }0一步转移概率矩阵与m 步转移概率矩阵的关系为:水利新闻"OU*T gR)U:D*ON0S

水利新闻5LA:nCn5cu Yd

P(m)=[P(1) ]m  (12)

v1E$` I,O7A/?0

,G1[5S4[3N*c0马尔可夫过程(有限状态的) 在经过长时间的转移后,会达到平稳状态,此后过程的状态将不再随时间的变化而变化。此时的概率被称为稳定概率,又称为极限概率,可定义为:水利新闻3ku5K D3P%R

DtK&oy&x;]*xFN0水利新闻8To0U `0a\.{ i

TI],a7^~6~ t0系统处于任意特定状态的概率可以由极限概率来代表,该特定状态重复再现的平均时间可用1/P表示。水利新闻.UI#W0p8gc:mS

水利新闻 [!at-Ws

2 结果及分析

%|/[AZ E p.cu0 水利新闻~8@2M-J X9HK

由于黄河河曲水文站是国家基本水文站、中部引黄工程水源地天桥水库入库站,本研究采用黄河河曲站1978-2010 年的实测时序年径流量作为基本数据,具体如图1 所示。水利新闻4N D+S4@sn2|U

-qW+e0W A'fq%K:H0I0

_F'mF7g0

0@wjm}'T5ws#y02.1 径流时序的年际变化趋势性

)DFv,X A4i0 水利新闻pjB2`*K

河曲水文站年径流量时序的累积距平曲线如图2 所示,由图可以看出: 1978-1986 年,年径流量呈现增加趋势; 1986-1994 年,年径流量呈现震荡趋势; 1994-2010 年,年径流量呈现下降趋势。根据Kendall(坎德尔) 秩次相关检验,当n=33时,计算的结果如下表1 所示。

O m.aO@S Ahy(X0

h Y!];N q$h0水利新闻^"}(csb(gckL

水利新闻 r J%AT"_2\#zLh

当n=33 时,取置信度水平α=0.05,Uα/2 =1.96。而年径流量序列的检验统计量U =-2.884,可见黄河河曲水文站年径流量时序在1978-2010年整个时段内呈现出显著的减少趋势。

QL|{&}5X&x0 水利新闻\&GqA7T$o

2.2 径流时序的年际变化持续性水利新闻txyB"{Y

水利新闻j` GW-C| ``U

应用R/S 分析法对河曲水文站年径流量时序变化进行分析,结果如图3 所示。所求得的Hurst指数为0.9176,0.5 < 0.9176 < 1,可以看出黄河河曲水文站的年径流量存在正持续性。由上节的分析可知,黄河河曲站年径流量时序在整个观测时段内呈现出显著的减少趋势。因而,可以推测出在未来黄河河曲站的年径流量将继续呈较显著减少趋势。水利新闻3`2C,^c7E8l B

)h)iGVzz0水利新闻sl_a'Ic%`z

3kTP N,K&|%w02.3 径流的连丰、连平与连枯概率水利新闻(vO)rA~9f*\

水利新闻"z k'E-t'E-@

采用常用的频率分析法来划分径流丰枯,表2为划分标准[11],划分结果如表3 所示。水利新闻 ?5mX?Xa uZ o#j6F-r

A#iM9[!e[^)Y0应用条件概率分析法对黄河河曲水文站的连丰、连枯概率进行分析,结果如表4 所示。

Co-G5O6D8MssEP0 水利新闻9VkN*l1_}5?

$dTu:r$dt hS$U0

o1h2oA8e0\ ^"]7i0

h d&RQ)cY0 水利新闻UB4@b+@rU0O

从表4 中可知:水利新闻6xY-Tw(tQ)W

水利新闻^,DZr$Yh;m}

(1) 当连续年数为2a时,连丰概率最大,为23.14%; 连枯次之,为18.75%; 连平最小,为13.89%。当连续年数为3a时,连丰概率最大,为14.73%; 连枯次之,为14.06%; 连平最小,为2.32%。当连续年数为4a时,连枯概率最大,为10.55%; 连丰次之,为9.37%; 连平最小,为0.39%。当连续年数为5 a时,连枯概率最大,为7.91%; 连丰次之,为5.96%; 连平最小,为0.06%。

!NJ Rf$T!P"|w0 水利新闻h j|7k7b^

(2) 连丰、连平和连枯状态的出现概率都随连续年数的增加而减小。

7uZtg^h\0 水利新闻 ]6_G4Z R Ssh]w

2.4 径流丰枯状态的马尔可夫过程分析水利新闻C yh\|,eA1UO

水利新闻1K*F]"i9ev

计算黄河河曲水文站年径流丰枯状态的结果如表5、表6 所示。

(R.V3DOtJfi T0 水利新闻9W9u;R_r

'~na-NM*?V\!R0

o;u m-u3W;K] M0从表5 可以看出:水利新闻lx gt0D!|oXF

水利新闻/nP4Gtau

(1) 从各态自转移概率的角度来分析,丰水年、偏丰年、平水年、偏枯年、枯水年的状态值依次为57.14%、0、16.67%、33.33%、57.14%。这说明丰水状态与枯水状态的自保守性最强,偏丰水状态的自保守性最弱。

^"O.|M*q LGF"w0 水利新闻7]vE+| DN

(2) 从各态互转移概率的角度来分析,丰水年向偏丰年转移的概率较大,为28.57%,偏丰年向各态转移的概率相同,均为25.00%,平水年向偏枯年转移的概率较大,为33.32%,偏枯年向平水年的转移概率较大,为33.33%,枯水年向偏枯年转移的概率最大,为42.86%。水利新闻b'ZI F e,m:P*u-kkn

水利新闻5rd k&zd;I'PA

(3) 年径流从各种状态向偏枯水状态转入的平均概率最大,为26.9%,向偏丰水状态转移的平均概率最小,为9.05%。这说明年径流量无论处于何种初始状态,其向偏枯水状态转移的概率均大于偏丰水状态转移的概率。水利新闻IJ7Y l`)w9Ml

DeW tH1e \%|0从表6 可以看出: 黄河河曲水文站出现丰水年、偏丰年、平水年、偏枯年、枯水年的概率分别为19.05%、8.32%、17.29%、28.81%、26.52%; 其平均重现时间分别为5.25、12.02、5.78、3.47、3.77 a。这说明在黄河河曲站年径流长期丰枯演化过程中,偏枯年的重现时间最短,偏丰年的重现时间最长。水利新闻 r:R!Ab D:DC.K5Oii

*U.c#|U#??%yTS03 结语

v\'Wn6eI ?0 水利新闻(kl7N#@%qH| t&E

(1) 黄河河曲站年径流量时序在1978-2010年整个时段内呈现出显著的减少趋势,并且可以推测出在未来黄河河曲站的年径流量将继续呈较显著减少趋势。水利新闻8@'p!p#hz!]J7qH

水利新闻Sk0b{#HJ(S

(2) 黄河河曲站连丰、连平和连枯状态的出现概率都随连续年数的增加而减小。水利新闻oNi9B[_ [-M @9t q

pAA3Ubz3[yP3tW0(3) 黄河河曲站年径流丰水状态与枯水状态的自保守性最强,偏丰水状态的自保守性最弱。该站年径流无论处于何种初始状态,其向偏枯水状态转移的概率均大于偏丰水状态转移的概率。在该站年径流长期丰枯演化过程中,偏枯年的重现时间最短,偏丰年的重现时间最长。水利新闻7ev3e-hC"cIm

+]]G5p#GBy%x,A0由此可见: 中部引黄工程水源地年径流量未来呈较显著的减少趋势; 在年径流长期丰枯演化过程中,出现偏枯年的可能性较大。

khG6qd[.p)p"|m0 水利新闻gd b8z%isWu K}n

参考文献:水利新闻*^Q/U zn;y*?

6\G;vzk6\{k {0[1]Shi Changxing,Zhou Yuanyuan,Fan Xiaoli,et al.A study on the annual runoff change and its relationship with water and soil conservation practices and climate change in the middle Yellow River basin[J].Catena,2013,100: 31-41.水利新闻?)K0v BVQ

水利新闻uFR"u/KU

[2]Wang Suiji,Yan Ming,Yan Yunxia,et al.Contributions of climate change and human activities to the changes in runoff increment in different sections of the Yellow River [J].Quaternary International,2012,282: 66-77.水利新闻%L$Z?HvF:Qt^

水利新闻6lY.mPK,[

[3]赵雪花,安莉莉,袁旭琦.基于HHT 和R/S 分析的黄河上游年径流序列演变模式分析[J].水电能源科学,2013,31(7) : 9-12.水利新闻m ~ c&n5p\}

水利新闻\"bcI,C!AS ? k`

[4]李二辉,穆兴民,赵广举.1919-2010 年黄河上中游区径流量变化分析[J].水科学进展,2014,25(2) : 155-163.水利新闻?k/_1X5?l*Y

y&t+w6IZj"d\0[5]安莉莉,赵雪花.黄河上游月径流序列多尺度特性分析[J].人民黄河,2012,34(11) : 26-28 + 31.水利新闻 X(w:K/`'gD

j*JI,D _%?Q~0[6]赵国辉,张耀南,蓝永超.黄河源区径流长期演变特征与趋势预测模型研究[J].冰川冻土,2010,32 (1) : 189-195.水利新闻}9RX~1P{S!IL

pGL tV ` Z0[7]李彦彬,尤凤,徐建新,等.黄河径流变化规律的小波分析[J].水利水电技术,2012,43(1) : 17-20 + 26.

1tyl|R:t.mX:S0

/jku2zk0[8]李明新,吕孙云,徐德龙.汉江上游水资源量变化趋势分析[J].人民长江,2008,39(17) : 49-52.

L&o|q}6b!dq6G0

N\/FeN? F4_0[9]丁志宏,冯平,牛军宜.黑河莺落峡年径流量时序变化的趋势特性及丰枯演化规律研究[J].干旱区资源与环境,2009,23(10) : 59-63.水利新闻 r#fr%F hW

NkE\/O?Q&I8d0[10]刘招,吴新,燕爱玲,等.关中地区主要河流丰枯变化特征分析[J].干旱区资源与环境,2013,27 (3) : 99-103.水利新闻 H)A*Ut.]

*Ig#gi;F oWVd0[11]张少文,张学成,王玲,等.黄河天然年径流长期丰枯状态变化特性研究[J].人民黄河,2005,27(5) : 9-10 +46.水利新闻Z#xw*o(Y _

水利新闻5gm b H$V6`{1N

作者简介: 徐聪(1990-),男,河南长葛人,硕士研究生,研究方向: 水环境模拟与预测。水利新闻"et:B5w#e$x g

顶:18 踩:19
对本篇资讯内容的质量打分:
当前平均分:-0.99 (102次打分)
【已经有82人表态】
36票
感动
4票
路过
5票
高兴
7票
难过
13票
搞笑
4票
愤怒
4票
无聊
9票
同情
下一篇:淮河流域海水利用发展与对策上一篇:非洲水图(一二一)

广告投放区


水资讯网接受会员自助投稿


2015年河海大学举办EFDC_Explorer7.2培训

广告投放区

广告投放区