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

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

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

kg1l\.lF0作者:徐聪 冯民权水利新闻[ _7Atnj)sE

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

水利新闻%oo6@ `7rxM5X

关键词: 年径流量; 时序变化; 丰枯演化; 中部引黄工程

vDP(Fx]c0 水利新闻8Jt@u.X qJ

中图分类号: P333.1 文献标识码: A 文章编号: 1672-643X(2014) 06-0105-04

0[To2p oL:_M0

/X)x_'oqFMUY,N _0Study on long-term evolution law of annual runoff in water source of Yellow river diversion project水利新闻|lZ F9Q(e_'G fj

$]%lg/['@,W&ka7B Bk0XU Cong,FENG Minquan

(r%km&? rc8I0 水利新闻o7O v7~*C/^CVL-E

Abstract: 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.水利新闻&p0U(]i8`

水利新闻j(p"x4W J6h-t1i

Key words: annual runoff; temporal variation; evolution of high and low state; central Yellow River diversion project水利新闻'[0~%b%nK.R;l7A

^ L P$yNX(w8z0黄河是中国的第二大河,干流长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]的一种或几种对黄河径流长期变化趋势进行了分析研究。但是,这些研究所采用的研究方法相对较为单一,其研究结果具有一定的局限性,不能全面的揭示黄河年径流量长期演变的时序变化特征及丰枯演化规律。

(DN6A&Jd[0 水利新闻dWBv~,d L

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

Q|T{D+Bz%}?/B0 水利新闻G0y7s#QtbYP.C

1 基本理论及方法水利新闻]:J.B'j6o2dT

水利新闻 zM(] ?(w6g?3c

1.1 Kendall 秩次相关检验法[8]

&R;F[%? x&J0

EB ^:P%i}_D0对于序列x1,x2,x3,…xn,先确定所有对偶值(xi,xj,j> i) 中的(xi < xj) 出现的次数,记为P,根据计算的P,代入下列公式:水利新闻2p"Fr$^D(`5t!xh

水利新闻4}WWy;Br%S

水利新闻mL;O[lIm

8v2g)v;WlD:Y0式中: τ为Kendall 秩次统计量; U为Kendall 秩次相关检验系数,随n 的增加很快收敛于标准状态分布;n 为系列长度。水利新闻^ t'K$[*i l(}{D

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

^D YJ$WH0 水利新闻6| GdNP%n;d#W2N c

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

"sC$tc2B ^0

#v\9Q2i0L^ xA01.2 R/S 分析法[3]

+N ^f~ r Z|0 水利新闻7h$Z0j_ ]&a6Rg

将R/S 分析法应用于径流时序特性分析时,判断径流时序状态的持续性,可通过计算Hurst 指数。

XP,v#dm)}%M&x0 水利新闻Yb^0Y6`0F @

R/S 分析法基本原理如下:水利新闻6j.y)B ?3}DV y!E

9h cf0W)HE A0对于时间序列{ x(t) } (t=1,2…,n),下列量可定义为:

Ld1K:v-b:C)oN0 水利新闻NPz:uY

均值序列:水利新闻 ? w`'w3z

7`gp ]*B0水利新闻 S7t4s:Y/w`

水利新闻-V'~+MbC%J A J

累计离差:水利新闻 v2shVnBi$m

水利新闻4A:?8{M |E&S sN-A

水利新闻 u;g]5A^A8Bu

水利新闻 YJ1X0tF.nh4J9O

极差序列:

2r3RH(Y0r9ToE;P H0 水利新闻7O-A-g+B8rS2R {{ D

vPaB4UM0

RV2`9z7tR!G0标准差序列:水利新闻 U j`;Q'F}:dz1BtI

V.m h/K@3GMi0水利新闻2^&u1f3I.c'H(G

水利新闻{?"] Bc6q)a)j Q_

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

/x;h;u+M,^0

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

?RET#b0

t"jQ Uq!uI u.| {01.3 游程分析法[9]

1Q8[3Cdb$QTe0 水利新闻+MPR@7?x DH

游程分析法的计算可用以下公式:

7X:GL2`3^rd v5nSq0

yzH^c;o!b0P=ρ(k-1) (1-ρ) (8)

3j3q@,x0N^q4DdB{0 水利新闻"\\$h_.~P5b"e(O0uO7]

ρ=(S1-S2)/S (9)

)U fd+n9}&{*prM([0

4B)FA4? aC6o{ Su0式中: P 为连续k 年枯水(丰水) 的概率; k 为连枯(连丰) 水年的年数; ρ为分布参数,其值在0 ~ 1 之间; S1为时间序列中连枯年(连丰年) 总年数; S2为时间序列中各种统计长度的连枯(连丰) 年发生的累积频次; S 为时间序列枯水年(丰水年) 的累积频次。

NDOI WZ[#PC0

t2g`j&Q]01.4 马尔可夫过程基本理论[10-11]

[ u'R&r f8b'F^0

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

\!MC)T G8W+[0 水利新闻&Z2BG[#u

水利新闻:U@nsAz-OZ

水利新闻\Eie)q5y WZ

(au$`6qEO0

vZ4@l-`wb0一步转移概率矩阵与m 步转移概率矩阵的关系为:

W6Y1o-g%A^5U~0 水利新闻(U\'~0|inT6CU'C

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

#y p Vp6iR6n|xv ?6i0 水利新闻R,{)w OFiC

马尔可夫过程(有限状态的) 在经过长时间的转移后,会达到平稳状态,此后过程的状态将不再随时间的变化而变化。此时的概率被称为稳定概率,又称为极限概率,可定义为:水利新闻@F'd6C&~;epf#E;f

p#aLs s7\;q/@E0水利新闻%k MsD"H1l0R

9WbW~~{+s7{0系统处于任意特定状态的概率可以由极限概率来代表,该特定状态重复再现的平均时间可用1/P表示。水利新闻 j ?3EPET4u

mI(sxv ~*^02 结果及分析水利新闻0I+|tM| t

水利新闻5q vBe3B h#B

由于黄河河曲水文站是国家基本水文站、中部引黄工程水源地天桥水库入库站,本研究采用黄河河曲站1978-2010 年的实测时序年径流量作为基本数据,具体如图1 所示。水利新闻]:Lfm!Y#[W q

水利新闻&ZQ'x9p_8e{

水利新闻5S2O;dX-MjQr ld&N

水利新闻/[@!J#HX3U[

2.1 径流时序的年际变化趋势性水利新闻{n1E[ Ih

水利新闻"[I/K8M~az

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

;k*n A\c%s(Fm0 水利新闻 Q!Q~|M@*j7uX

,W]x!Q,v D Ua0 水利新闻*}8V'F a'^

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

] mZ,["Q(w(Z0 水利新闻 ?S.CQ2RT;_

2.2 径流时序的年际变化持续性

3_bq(e+L0

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

r*Ys0eJ+x!L0

;e~1X0W h3Y0水利新闻-n4m C5h,t@0aj,n,E

水利新闻7k)UA#\$u

2.3 径流的连丰、连平与连枯概率水利新闻 m6nOka+GX`

水利新闻N {`4PQa

采用常用的频率分析法来划分径流丰枯,表2为划分标准[11],划分结果如表3 所示。

l%T hS/za0

(I%Cylae0应用条件概率分析法对黄河河曲水文站的连丰、连枯概率进行分析,结果如表4 所示。水利新闻h lD&`"aifdb

7v BQ|Y0水利新闻/g? O|1u)e9Hm:H

水利新闻w!O,c"O^2B;?

水利新闻&}{BxB0Y8C

水利新闻{ k3k:k1cN(U

从表4 中可知:水利新闻&Oq\c_W r

水利新闻obWWRh)o!]

(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%。

$zX0{]Hv0 水利新闻`g!O`*o[sKk

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

*sAIu2? WO `YL0

wx,| L0J|*n^02.4 径流丰枯状态的马尔可夫过程分析水利新闻b#rMOv

水利新闻%S$t(fKVr

计算黄河河曲水文站年径流丰枯状态的结果如表5、表6 所示。水利新闻%~)i O'} W`

IK/c5Tb-f H-Os y2eg0水利新闻/hJ2Vd~-?\%k

2Z2uHj8~3o U,W!\u l0从表5 可以看出:水利新闻2g%w'L*G(J

水利新闻z t"~y%Z/h$o

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

9}cg8v~!?xtc J0 水利新闻j1x e@+u1b%v {&o

(2) 从各态互转移概率的角度来分析,丰水年向偏丰年转移的概率较大,为28.57%,偏丰年向各态转移的概率相同,均为25.00%,平水年向偏枯年转移的概率较大,为33.32%,偏枯年向平水年的转移概率较大,为33.33%,枯水年向偏枯年转移的概率最大,为42.86%。水利新闻 n5N-~Xn%~)x:@

水利新闻-x[Ru+c'b%B

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

XX.V|:H0

ZbJ"D%K%b]0从表6 可以看出: 黄河河曲水文站出现丰水年、偏丰年、平水年、偏枯年、枯水年的概率分别为19.05%、8.32%、17.29%、28.81%、26.52%; 其平均重现时间分别为5.25、12.02、5.78、3.47、3.77 a。这说明在黄河河曲站年径流长期丰枯演化过程中,偏枯年的重现时间最短,偏丰年的重现时间最长。

"fZ#E3r aB P~5a0 水利新闻l[,`j\3{

3 结语水利新闻dZ&E7l5BC&]/W

水利新闻 K3hVP } ^S

(1) 黄河河曲站年径流量时序在1978-2010年整个时段内呈现出显著的减少趋势,并且可以推测出在未来黄河河曲站的年径流量将继续呈较显著减少趋势。

6{oC&Z+P-r _.{V0 水利新闻.u5A`DT!w

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

Ogcs&m3E K0

'LZ7gu{'kx0(3) 黄河河曲站年径流丰水状态与枯水状态的自保守性最强,偏丰水状态的自保守性最弱。该站年径流无论处于何种初始状态,其向偏枯水状态转移的概率均大于偏丰水状态转移的概率。在该站年径流长期丰枯演化过程中,偏枯年的重现时间最短,偏丰年的重现时间最长。

B x:CH8J{P j:I0 水利新闻~cJ$O2h3Hg

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

I+y2iGGd0 水利新闻+G6z `4A!Sy"X

参考文献:水利新闻jIS x6q'i&d

水利新闻WcBUq-~UH_E)u

[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.水利新闻x.o;{yzX

水利新闻#|mkqTv

[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.

^!q(C3AX*Pb0 水利新闻 QJA LL]A7{lt

[3]赵雪花,安莉莉,袁旭琦.基于HHT 和R/S 分析的黄河上游年径流序列演变模式分析[J].水电能源科学,2013,31(7) : 9-12.水利新闻A,iZ'I/zT r"C\n

水利新闻%k2x~9TP)U2Po

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

水利新闻 D^mJ'`l.W#XA5~C

[5]安莉莉,赵雪花.黄河上游月径流序列多尺度特性分析[J].人民黄河,2012,34(11) : 26-28 + 31.

]7V.gw j~0

-f"^O/}'C'f0[6]赵国辉,张耀南,蓝永超.黄河源区径流长期演变特征与趋势预测模型研究[J].冰川冻土,2010,32 (1) : 189-195.

zpt8o-MIm:Z+@0

)SgJR HF0[7]李彦彬,尤凤,徐建新,等.黄河径流变化规律的小波分析[J].水利水电技术,2012,43(1) : 17-20 + 26.水利新闻$w8NE"Z!v kO;^k/V@ S2o

水利新闻9Ob$??5eiq]9o

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

w"b(SAD0[9]丁志宏,冯平,牛军宜.黑河莺落峡年径流量时序变化的趋势特性及丰枯演化规律研究[J].干旱区资源与环境,2009,23(10) : 59-63.水利新闻 YE5vL%[:yhh*t xt

9Lp(r!^'U(G`0[10]刘招,吴新,燕爱玲,等.关中地区主要河流丰枯变化特征分析[J].干旱区资源与环境,2013,27 (3) : 99-103.

\nFedv0 水利新闻)Q%C#h8D+d'n:J#F

[11]张少文,张学成,王玲,等.黄河天然年径流长期丰枯状态变化特性研究[J].人民黄河,2005,27(5) : 9-10 +46.水利新闻"\wT\@

3\qEBQ2n!c0作者简介: 徐聪(1990-),男,河南长葛人,硕士研究生,研究方向: 水环境模拟与预测。

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

广告投放区


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


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

广告投放区

广告投放区