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

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

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

-D ]Vy-]`m.deEK:k0作者:徐聪 冯民权

9wQ,@q9]Si0

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

m/Z,S ~DO2A7Qu*X0 水利新闻 P(C r3V,j L8R0Cm

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

}9^s6}h0 水利新闻Ik&bU/mX

中图分类号: P333.1 文献标识码: A 文章编号: 1672-643X(2014) 06-0105-04水利新闻#Y-TSo z0ssN@y

水利新闻9M%m p1N;{8H

Study on long-term evolution law of annual runoff in water source of Yellow river diversion project

&buLz-rm zh6C0 水利新闻4~H8L7G1D v%t

XU Cong,FENG Minquan水利新闻.rAA9`n#E-@3WLV

q'E:ePf `5DxzpH0Abstract: 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.

0FXfBBID&k j4| G0

6i+pNcRuL0Key words: annual runoff; temporal variation; evolution of high and low state; central Yellow River diversion project

xn@b&Y W/m7S e0

X F%fV+F[ U0黄河是中国的第二大河,干流长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]的一种或几种对黄河径流长期变化趋势进行了分析研究。但是,这些研究所采用的研究方法相对较为单一,其研究结果具有一定的局限性,不能全面的揭示黄河年径流量长期演变的时序变化特征及丰枯演化规律。水利新闻8C2IhRy Yfb

水利新闻y-g4RK,Q

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

j-B ud(x'g vi0

\q5@f)V7[01 基本理论及方法水利新闻9j9oDg a7k|3P\

水利新闻,c%GR jK-{"w{V

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

@ Q`Z+Jl0 水利新闻SsJ*w}J2`

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

`+aEJn.{K0

ms;FT8FH'm,J iZ0水利新闻&v'u B6~ n

6i;_#t:`Kr4ig0式中: τ为Kendall 秩次统计量; U为Kendall 秩次相关检验系数,随n 的增加很快收敛于标准状态分布;n 为系列长度。水利新闻9@8S3hG k

RYCe oa6L0若给定显著性水平为α,则其临界值为Uα/2。若| U | > Uα/2,则说明序列具有明显趋势性; 当U > 0时,则具有上升趋势; 当U < 0 时,则具有下降趋势。水利新闻rV T]'F

水利新闻 Rd n?2KoRl5Tj

若| U | < Uα/2,则序列趋势性不明显。水利新闻`VaFbq

|;aN'|[ a Q01.2 R/S 分析法[3]水利新闻 k%I~)Y,p(xVH

水利新闻)jRY%AXl S#X.bO

将R/S 分析法应用于径流时序特性分析时,判断径流时序状态的持续性,可通过计算Hurst 指数。水利新闻`\ |)o/kg2k5a

水利新闻3_{Of3^ U7F(b

R/S 分析法基本原理如下:

&x6K.N8j?^0

FQz [xk0对于时间序列{ x(t) } (t=1,2…,n),下列量可定义为:

n-n q3yBR0

HI-v$[ c6r`0均值序列:水利新闻;NQ-F&u)\Y'Z/T9b$ka~

6}^$gk eV0水利新闻 Q)S6`4W7R

;Y};lwD2`5_K0累计离差:

g:lJ.cL-y1p1];u0

O7{!c%b d.I9B0水利新闻 F-|F.f.tm't

水利新闻4Wn-j0B@M D

极差序列:水利新闻SW-u9QPn.[g/R

水利新闻Y qLz5[xhX,Bo

4A8MjX%a*q)s0 水利新闻*h N'd5b-~{2b!x.U%j9J

标准差序列:

IJ wry)J,r's0p0

N0w2u.C%Cv(^6}e0水利新闻 Hd Op1K6{ K

N`Lo m;Mpate:s[0Q0如果比值R(τ)/S(τ) ≡ R/S 存在关系: R/S ∝τH,则说明时间序列{ x(t) } (t=1,2…,n) 具有Hurst 现象,H 称为Hurst 指数。在双对数坐标系(lg(τ),lg(R/S) ) 中用最小二乘法拟合,拟合直线的斜率即为H。水利新闻P`4z~M6j9d

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

3M,k)L"xY(}dHyfd0

htx9u6j2ds)Y1EU01.3 游程分析法[9]水利新闻+y iy qf1YF

水利新闻"m Za~ m;}3u

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

nj0Ig:h l0 水利新闻_g,fQ4E5c

P=ρ(k-1) (1-ρ) (8)

,hK-zDw6_a0

w5o j*qD*e i0ρ=(S1-S2)/S (9)水利新闻 L'[p;kPH4a

水利新闻| k#]3U8CG6pyV

式中: P 为连续k 年枯水(丰水) 的概率; k 为连枯(连丰) 水年的年数; ρ为分布参数,其值在0 ~ 1 之间; S1为时间序列中连枯年(连丰年) 总年数; S2为时间序列中各种统计长度的连枯(连丰) 年发生的累积频次; S 为时间序列枯水年(丰水年) 的累积频次。水利新闻 w,W }P \}OA

CbkV3Qv+J+j01.4 马尔可夫过程基本理论[10-11]

*rG1e z2kEq C;G0V0

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

y*h_2V5};f]k0

)hlf&]0TH0

R3g3[ Qoz7t!w N0水利新闻1Y;qO!IQ

水利新闻2Ro'p Rfn(i

一步转移概率矩阵与m 步转移概率矩阵的关系为:

/{&@&Jx\N&XxX0

m)| Ti G$C l/y0P(m)=[P(1) ]m  (12)水利新闻"Ni$q(I*p}/_*Q1hP&I

水利新闻]1V!mKJ9IC

马尔可夫过程(有限状态的) 在经过长时间的转移后,会达到平稳状态,此后过程的状态将不再随时间的变化而变化。此时的概率被称为稳定概率,又称为极限概率,可定义为:

8E5N:}#XWs8o;i-G0

K[-eCG|0水利新闻j%a7g+nMe Sb"\

水利新闻#p5\t-i2AS#i j6`}"a

系统处于任意特定状态的概率可以由极限概率来代表,该特定状态重复再现的平均时间可用1/P表示。水利新闻:^$F#~0H @6E3T-lf[1R#R]

.] |Rb%^"|1@02 结果及分析

BA7X(Ep-K|Y.vM%g0

{%\:Dviy(PN |0由于黄河河曲水文站是国家基本水文站、中部引黄工程水源地天桥水库入库站,本研究采用黄河河曲站1978-2010 年的实测时序年径流量作为基本数据,具体如图1 所示。水利新闻7I_e6O mR.G0^S

水利新闻,^c%JHQW

0tUo S2`0

[2y~6hJ,Ux?}02.1 径流时序的年际变化趋势性水利新闻:Jt%`-K!o8@

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

P}Rk sA2p4LT0水利新闻BK)J I;u

水利新闻HQ^)u$W L5C(j

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

[O2W6^3y#ZO%d02.2 径流时序的年际变化持续性

p4o7q@Y U5jW(Q#~0 水利新闻y Bz+l4C3^

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

水利新闻:q3p(QP }7|"q

rr(@7v4~0 水利新闻0\"]QY-?L,{ wf*}!h

2.3 径流的连丰、连平与连枯概率水利新闻/H9bT)LMZ#d3yy

水利新闻/K@{}S;Hb

采用常用的频率分析法来划分径流丰枯,表2为划分标准[11],划分结果如表3 所示。水利新闻 nr+b%f v%H Ae

水利新闻E!t0mwB

应用条件概率分析法对黄河河曲水文站的连丰、连枯概率进行分析,结果如表4 所示。水利新闻*MV0^q#TQ)t

+| U` @0e,I(Ra f-n0

8OmD:AxS'h0

t3_c y-{w%jl [0水利新闻#bw9}Y0F%y

I!MWv6h8@G0从表4 中可知:水利新闻 ?X,mUv

#V3Oh)?1YZA R0(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%。水利新闻3iZd4z'O(F&os#NTS

水利新闻h5uc$DG1K!`

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

%uWS)Za`6Yy"@9@V0

Z!Os P)c nYU5G02.4 径流丰枯状态的马尔可夫过程分析水利新闻:R!jB1|Q-d3v&Y

水利新闻2E6n{@`&P UP

计算黄河河曲水文站年径流丰枯状态的结果如表5、表6 所示。水利新闻5rOZ ^.\k

水利新闻)|M"hZi%B(Z

:\7M7^3GS D0

u5g:cx S2b c2a0从表5 可以看出:水利新闻5y(~C;l'x,ft'qy v

水利新闻Z _9Z/q[No(d3?

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

9`9k x ^2Dx1|{f(Z0 水利新闻7v+BI]Sg5q N.{

(2) 从各态互转移概率的角度来分析,丰水年向偏丰年转移的概率较大,为28.57%,偏丰年向各态转移的概率相同,均为25.00%,平水年向偏枯年转移的概率较大,为33.32%,偏枯年向平水年的转移概率较大,为33.33%,枯水年向偏枯年转移的概率最大,为42.86%。

5L$`Z1Fp/o3q0 水利新闻ZG-fR1DVf

(3) 年径流从各种状态向偏枯水状态转入的平均概率最大,为26.9%,向偏丰水状态转移的平均概率最小,为9.05%。这说明年径流量无论处于何种初始状态,其向偏枯水状态转移的概率均大于偏丰水状态转移的概率。水利新闻o~ k$O:\*F9Qs

水利新闻\i[ k,X}s%x] ^

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

-]]v!~$c.N?%I Q0

e.K_6e6z\T,Qfl^03 结语

6g:|9LS.F;~C0 水利新闻vWrPD

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

lSjA4W YM+b {0 水利新闻6z}:C1Ut"KZ

(2) 黄河河曲站连丰、连平和连枯状态的出现概率都随连续年数的增加而减小。水利新闻 f`3x$d.C!|M.Jm5\*V

X/KR%^ Bo p?'w0(3) 黄河河曲站年径流丰水状态与枯水状态的自保守性最强,偏丰水状态的自保守性最弱。该站年径流无论处于何种初始状态,其向偏枯水状态转移的概率均大于偏丰水状态转移的概率。在该站年径流长期丰枯演化过程中,偏枯年的重现时间最短,偏丰年的重现时间最长。水利新闻 R R6_8M?+oh

水利新闻ylNZ+fW

由此可见: 中部引黄工程水源地年径流量未来呈较显著的减少趋势; 在年径流长期丰枯演化过程中,出现偏枯年的可能性较大。水利新闻(QW z\2bq}bw d;m

+k3A X,F$K jo0参考文献:水利新闻iBXiUp#wL

z5@J&VHy0[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.

*W?/NO/a"zu/k0 水利新闻pmW:H3`

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

7N MS+M"_y0 水利新闻!P[.CQk| l

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

B"}[ RV%u+y \7l0 水利新闻pq\3E%o/l

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

水利新闻2cS+s&Wvwe(E

[5]安莉莉,赵雪花.黄河上游月径流序列多尺度特性分析[J].人民黄河,2012,34(11) : 26-28 + 31.水利新闻:d(v*xG9Q`` d(q;[4k:h

水利新闻c)| @J+l}-S S?f~

[6]赵国辉,张耀南,蓝永超.黄河源区径流长期演变特征与趋势预测模型研究[J].冰川冻土,2010,32 (1) : 189-195.水利新闻oW9eJ QVo(I:~a

水利新闻s R L3vZ&_

[7]李彦彬,尤凤,徐建新,等.黄河径流变化规律的小波分析[J].水利水电技术,2012,43(1) : 17-20 + 26.水利新闻6\@"oS)M%u*gH

水利新闻q9r/}]0L:SLYk}

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

2y9bD\x1y#_k%NJ+s_0 水利新闻m BMChd(d)VqF

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

luY&Xn ^0 水利新闻l-[E)T/NM"C

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

N6Opd fL,W0

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

_o"}({4O&s)e$i&s0

4f E^.i/^tmI8z:n0作者简介: 徐聪(1990-),男,河南长葛人,硕士研究生,研究方向: 水环境模拟与预测。水利新闻'X2Ke\:y%QVD!N4y&i

顶:17 踩:13
对本篇资讯内容的质量打分:
当前平均分:-1.48 (73次打分)
【已经有72人表态】
31票
感动
4票
路过
4票
高兴
6票
难过
12票
搞笑
3票
愤怒
4票
无聊
8票
同情
下一篇:淮河流域海水利用发展与对策上一篇:非洲水图(一二一)

广告投放区


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


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

广告投放区

广告投放区