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

太湖流域降水不均匀性研究

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

Ma4q-[ h3Z0作者:邓鹏鑫 胡庆芳 王银堂 王磊之 崔婷婷水利新闻]&z0b6fQH EG

水利新闻cI!~}!G1dz2R

摘要: 基于1951—2011年逐日降水数据,采用信息熵分配紊乱指数、基尼系数和变差系数3 种不均匀性指标,分析旬、月时间尺度上太湖流域及各分区降水分布的不均匀性。比较3 种不均匀性指数的差异,并探讨指数与流域洪涝灾害的相关关系;同时采用基于预置白处理的Mann-Kendall 方法和Hurst 系数,诊断降水不均匀性指数的变化趋势及持续性,分析对流域旱涝灾害情势的可能影响。结果表明:3 种不均匀性指数间的差异较小;当年或汛期降水总量和不均匀性指数值均较大时,太湖流域较容易遭遇洪涝灾害,而在降水总量较小且不均匀性指数值较大时,则容易遭遇干旱灾害。在月时间尺度上,太湖流域汛期降水的不均匀性将持续增加,降水分布具有更为集中的态势,发生洪涝灾害的风险将有所增加。因此,有必要继续分析降水在空间上的不均匀性,以揭示降水在空间尺度上的分布差异,继而为流域防灾减灾提供技术支持。

Uvs"e:~"N*b0

Rmh/yI-}/Lkd0关键词: 太湖流域; 降水; 不均匀性; 基尼系数; 信息熵分配紊乱指数;旱涝灾害

)E.`|&s p(A0

xH D.Q/B0中图分类号: P338 文献标志码: A 文章编号: 1009-640X(2014)05-0034-07水利新闻 W x^ce^q.`

水利新闻!O0a Cp(S2G#hF#k

在气候变化和城市化背景下,我国区域降水要素的变化及对旱涝情势的可能影响引起了广泛关注[1]。然而已有的研究多关注全年或年内特定时段降水总量的变化特征,对年内降水时程分布特征变化的研究相对较少[2-6]。而某一区域的水旱灾害,既与降水总量有关,也与降水年内分布密切相关。因此,在分析降水总量变化的同时,还应分析降水要素在年内分布的时间不均匀性及可能变化。水利新闻(B{Y3T4TCC~o7Y

水利新闻zH0~9ZRX

太湖流域是我国典型的平原河网地区,位于东亚季风气候区,降水时空分布不均,水旱灾害比较频繁。研究流域降水的时间不均匀性,可以更全面深入地理解太湖流域降水的时程分布特性,为流域水资源利用、水旱灾害分析提供科学依据。因此,本文同时采用信息熵分配紊乱指数(Apportionment entropy disorder index,ADI)[7-8]、基尼系数(Gini coefficient,GC)[9]和变差系数(Coefficient of variation,Cv)3 种不均匀性指数,在旬、月时间尺度上分析全流域及分区域降水要素的不均匀性,并探讨不均匀指数变化对流域遭遇旱涝灾害情势的可能影响。

BYQ2xvbWUj0

HW j(b7B01 研究区域与数据水利新闻Z l4Li2[

-fnABI!AI01.1 研究区域水利新闻6]6t(Y'k5~7XHaOQ

水利新闻 A,D'[f T4L,p.Wm-N

太湖流域位于东经119°11'~121°53',北纬30°28'~32°15'之间,北抵长江,东临东海,南濒杭州湾,西以天目山、茅山等山区为界,流域总面积约36 895 km2。流域地形呈西高东低、四周高中间低之势,地貌特征以平原为主,山地丘陵仅分布在流域西部(图1)。太湖流域属典型亚热带季风气候,冬季干冷、夏季湿热,年均降水量为1 189 mm(1951—2011)。降水年内分布不均,汛期多年平均降水量为710 mm,约占全年的60%[10-11]。其中,5—7月为梅汛期,8—9月为台风期。

U0@)Asr7\%g0

#A-f|0xx6o j o E Xw~F0

%J3DG:]'B3W0

ALYpo(J$^q$h-{f0根据流域地形及水系特点,将全流域分为湖西、浙西、太湖湖区、杭嘉湖、武澄锡虞、阳澄淀泖和浦东浦西7个分区。水利新闻.Q/Ge$y/n{_S

pHN'he01.2 研究数据

9i3Rnb6N0

;o#z7YG*x \\z0收集全流域和7个分区1951—2011年逐日面降水量数据(太湖流域管理局水文局提供)。降水资料经过了质量控制,包括一致性、可靠性和代表性的分析处理。在逐日降水数据基础上,得到全流域及各分区年、汛期、月、旬降水数据。

W"IzN a@0

0]X\ U)L ]02 研究方法水利新闻!Av{5F u

水利新闻.xlTJ-nu g BXC

通常多采用变差系数(Cv)对降水的不均匀性进行分析,这一指标是描述样本相对均值的离散程度的一般性指标。近年来,一些学者相继探讨采用信息熵紊乱指数(ADI)[7-8]、基尼系数(GC)[9]等指标分析降水的不均匀性。其中,ADI 能揭示序列变化的深层随机性,起源于热力学、统计力学领域;GC 最早由意大利经济学家基尼(Gini)于1922年提出,用于社会财富分配平等性评价[12-15]。两项指标在降水不均匀性分析上已有成功案例[7-8]。

&Z/iz8qX R"~~0

nPqxo!s0Cv采用序列均方差δ与均值x(—)的比值衡量序列相对离差程度表示:水利新闻[Y4wWK&FZ2mg/u9BE

bw,zk5_D m0(1)

_Y6c'I Ai#V'lXG+g0

+d+USf `2U+u0式中:xi(i = 1,2,3,…,n)为样本序列;x(—)为序列的均值。Cv为实数。对于降水要素而言,Cv值越大,降水分布的不均匀性程度越高;反之,不均匀性程度越低。

-B&q u)~"tvW0

1V N/E@Z|0GC 主要通过构建洛伦兹曲线进行不均匀性的定量评价[12-13]。流域降水时间分布不均匀性的洛伦兹曲线的构建方法可参考文献[9]。在对比多种GC 计算方法[14]的基础上,本文采用张建华提出的简易公式[15],该公式利用定积分的定义将洛伦兹曲线的积分分成n个等高梯形:

o0Ele2Z(l0

7Sx:Ai8ob0(2)水利新闻p;Y@)xGW

g3|Uy.G{b0式中:Wt(t = 1,2,3,…,n-1)为累积降水量的百分比,即将降水序列从小到大排列,依次计算累计降水量占总降水量的比值。水利新闻Hs0qV{*[3m7U ~uu

水利新闻/]6M%H7Q7]s

ADI 是建立在信息熵概念的基础上,揭示要素深层随机性的指标。根据文献[7],对于一组有限的随机序列x,信息熵H 可表示为:

r S\4iX#o0

z Ib~!p }6]%C8Cn0(3)

7K.@"M,_.gynI"b t0 水利新闻;D y(hG1S Gi_(?C0A P;m

式中:xk(k = 1,2,3,…,n)为样本序列;P(xk)为样本序列xk对应的经验频率。水利新闻w$i!t$zq v

水利新闻.|Mq`*p&Z

在信息熵定义的基础上,进一步计算分配熵(Apportionment entropy,AE):

8vfu^gheP/r0 水利新闻z+SxBwU

(4)

_%u I!P o0

(GP#Ob+z6qzVFlh0式中:rk为时段k 对应的降水量(k = 1,2,…,6);时段降水总量R=。对于月或旬降水序列,n=12或36,AE 的最大值为log2(12)或log2(36),最小值为0。在AE 的基础上,文献[8]又引入了信息熵分配紊乱指数,即:水利新闻Z T$L2Fg9N3R GK(w m

/pi:RV6m*uXM+t(yw ]0ADI = Max(AE)-AE (5)水利新闻O E7A-Y!p*D aS-T

水利新闻h3]G/P X P:?,OAz v

本文采用带预置白处理的Mann-Kendal 趋势检验方法TFPW-MK ( Trend-free pre-whitening Mann-Kendall)分析降水不均匀性指数的变化趋势。TFPW-MK 与通常的MK(Mann-Kendal)方法相似,在给定某一置信水平时,当统计量的绝对值大于标准正态分布临界值时,则认为原序列趋势显著;否则,原序列趋势不显著。相比MK,TFPW-MK 可消除序列正向自相关性,避免趋势显著性被放大,因而在序列数小于50,且变化趋势不是非常明显的情况下对趋势性的检测更为有效[17]。具体方法步骤可参考文献[16]。此外,引入Hurst 系数来反映趋势性的持续性,当Hurst 系数大于0.5,则认为未来趋势与当前趋势保持一致[7]。水利新闻'?y-Mcj

水利新闻f]CrE

3 结果分析与讨论水利新闻V%t#qu4WC`

水利新闻Qd_t Z'r

3.1 3 种降水不均匀性指数的比较

y aM3u/x+G b0

-L7Y2Txx ~-M3@ }0本文基于全流域及7个分区1951—2011年逐旬、逐月面降水量时间序列,得到相应的ADI,GC 和Cv,并计算各指数间的相关系数。表1 中r1表示GC 与Cv的相关系数,r2表示ADI 与GC 的相关系数,r3表示ADI与Cv的相关系数。

2oQQ%F#yN{?D0 水利新闻].| `$_2]-Sn;C

表1 3 种不均匀性指数的相关系数

0vmI{AQU't [P-Wg0 水利新闻A2v \1c,U)w

+CsA"Ex;u1y8x0 水利新闻9DHf0t)yG@1Z U

由表1 可知,3 种不均匀指数间的相关程度较高,仅武澄锡虞区和湖西区在全年旬尺度下,GC 与Cv的相关系数值略低外,其他各分区指数间相关系数均大于0.93,其中,ADI 与GC 和Cv的相关系数大于0.95,最高可达0.98,整体反应出各指数间的差异很小。这与文献[8]的结果一致。但文献[8]认为,ADI 基于概率密度函数的离散,避免了受到均值的限制,更能揭示要素的深层随机性。

(W ?P3d(^{b1I0

ZUm9j I+A^5y0本文以ADI 为例,对降水的时程不均匀性作进一步阐述。其他指数的结果类似,本文不再赘述。

tCgCZ8UZ.q0 水利新闻(j`!G)s'P

3.2 流域降水不均匀性指数与旱涝相关关系

wJNg1eG&}0

4ABC[ yL0旱涝灾害不仅与年或汛期降水总量有关,而且与降水的年内分布不均匀性密切相关。本文针对太湖流域年或汛期降水总量和不均匀指数与流域旱涝灾害的相关关系进行了初步分析。图2 给出了1951—2011年全流域ADI 值、标准化降水(即年或汛期降水量与多年最大值之比) 的变化曲线。由图2( a) 和( c) 可知,在月尺度上,ADI 值和年或汛期降水量在1954,1957,1962,1983,1991,1999 和2010年均出现了较大值。由于降水量和ADI 均较大,月降水分布不均,某些月份集中过量降水,较容易遭遇较大的洪涝灾害。如1954年、1991年和1999年均发生典型的梅雨型洪水。同时,ADI 在1967,1971,1978,1994,2003 和2011年出现较大值,而年或汛期降水量低于多年平均值。在降水总量较小,而ADI 值较大的情况下,月降水量偏少,且分布不匀,较易产生干旱灾害。如1967 和1978年为流域特枯干旱年。此外,若遭遇旱涝急转年时,年降水量偏低且过量集中于汛期,易造成月降水在汛期的分配较年分配均匀,呈现出在年或汛期分配下ADI 值均较大的情况下,汛期降水量大于多年均值而相应年降水量却低于多年平均值。如太湖流域在2011年发生了60年来最严重的干旱以及典型的旱涝急转,其相应ADI 值和汛期降水量高于多年平均值,而年降水量较多年平均值低。水利新闻v3b@/H+k8gq,b

+@a b\EUB!b|E4X0水利新闻"aZ`/b `

Z+U2?6btns0旬尺度下,流域降水不均匀性指数与旱涝的指示关系除与月尺度具有一致性外,还易反映出短历时强降水带来的灾害影响。由图2(b)和(d)可知,相比月尺度,旬尺度的降水分配差异更加显著。汛期旬尺度下,1963,1970,1991 和2001年的ADI 较汛期月尺度大,表明汛期旬降水分布较集中,遭遇短历时强降水,特别是台风型降水洪涝灾害的可能性较大,如1963 和1970年为台风型洪涝年。水利新闻t6T5vQ(UUq+J}

-R ]Y&nh8J3\ ]U0可见,结合降水总量,可初步反映出不均匀指数与旱涝具有一定的相关关系;但其与旱涝灾害对应的精确性还需作更深一步的探讨。水利新闻Q~/O+p%V|$h Cg6L1G

r{2j*Vj Z/Z@03.3 降水不均匀性指数的变化特征水利新闻.U9Yf9v G;VTp

R?c/q.m {6a4RM0表2 给出了不同时间尺度下ADI 的均值和极值。在旬时间尺度上,太湖流域全年和汛期的ADI 最大值分别为1.12 和1.38,最小值仅0.38 和0.17,极值比为2.95 和8.12;而月尺度下,全年和汛期的ADI 极值比高达5.93 和24.0。在各分区中,北部的武澄锡虞区和湖西区的特征值较大,全年旬尺度下的均值达到0.84和0.80,极大值为1.32 和1.65,而南部浙西区和杭嘉湖区较小,相同尺度下均值仅0.67 和0.71,极大值为1.31 和1.32。这反映出太湖流域旬、月降水量在年内和汛期分配的不均匀性在不同年份的差异较大,且降水具有一定的空间分布差异,北部武澄锡虞区和湖西区降水的年内分布的不均匀性较其他各分区大。水利新闻\f3F n0X j O-W

u W i,lF$u0表2 各分区及全流域不同时间尺度下信息熵分配紊乱指数的特征值

)oJ'^)d3hy0 水利新闻P)@p5p9KP

Im.V v6\m1K0 水利新闻U+vy9X{K

采用TFPW-MK 和MK 两种非参数检验方法,结合Hurst 系数[8],进一步检验了不均匀性指数的变化趋势性以及趋势的持续性。由表3 可知,TFPW-MK和MK统计结果比较一致。两种方法都说明在汛期月尺度下,除湖西区和浦东浦西区外,其他各分区ADI 指数均具有显著上升趋势,其中浙西区、太湖湖区、武澄锡虞区的TFPW-MK 值甚至超过95%的置信度,Hurst 系数分别为0.71,0.64 和0.65,均大于0.5,上升趋势的持续性较强。

0R%AX`x;Q;y U;{gn0

\2kV5y,m0z|0本文进一步对汛期各月降水量的变化趋势进行了检测,结果发现太湖流域5 和9月降水呈现显著下降趋势,汛期降水愈趋集中于6—8月。这是汛期降水不均匀性指数增强的主要原因。而这也说明汛期流域降水分布趋于更集中的态势,较容易导致汛期洪涝灾害的发生。水利新闻0x3Tp7hZ*R&P8E

水利新闻HofXmge j

表3 信息熵分配紊乱指数的Hurst,TFPW-MK 和MK 计算结果水利新闻 i&L j[f,i A GU

M/^y)w9x%k0水利新闻py!n CA;IkV

水利新闻 P:g)W K\E$o P6b

旬尺度下,流域降水不均匀性指数与旱涝的指示关系除与月尺度具有一致性外,还易反映出短历时强降水带来的灾害影响。由图2(b)和(d)可知,相比月尺度,旬尺度的降水分配差异更加显著。汛期旬尺度下,1963,1970,1991 和2001年的ADI 较汛期月尺度大,表明汛期旬降水分布较集中,遭遇短历时强降水,特别是台风型降水洪涝灾害的可能性较大,如1963 和1970年为台风型洪涝年。

O tSY:Y7p/f@0 水利新闻 E:G9IR'b-TR:X

可见,结合降水总量,可初步反映出不均匀指数与旱涝具有一定的相关关系;但其与旱涝灾害对应的精确性还需作更深一步的探讨。

/s@^\aUE7wp ?0 水利新闻C\3e4{-S(L5Tz

3.3 降水不均匀性指数的变化特征水利新闻/Qg1K q.kLAh C

水利新闻_J:ipo"f

表2 给出了不同时间尺度下ADI 的均值和极值。在旬时间尺度上,太湖流域全年和汛期的ADI 最大值分别为1.12 和1.38,最小值仅0.38 和0.17,极值比为2.95 和8.12;而月尺度下,全年和汛期的ADI 极值比高达5.93 和24.0。在各分区中,北部的武澄锡虞区和湖西区的特征值较大,全年旬尺度下的均值达到0.84和0.80,极大值为1.32 和1.65,而南部浙西区和杭嘉湖区较小,相同尺度下均值仅0.67 和0.71,极大值为1.31 和1.32。这反映出太湖流域旬、月降水量在年内和汛期分配的不均匀性在不同年份的差异较大,且降水具有一定的空间分布差异,北部武澄锡虞区和湖西区降水的年内分布的不均匀性较其他各分区大。

9@ FFZ,v R#G0

IG8YT:F0采用TFPW-MK 和MK 两种非参数检验方法,结合Hurst 系数[8],进一步检验了不均匀性指数的变化趋势性以及趋势的持续性。由表3 可知,TFPW-MK 和MK 统计结果比较一致。两种方法都说明在汛期月尺度下,除湖西区和浦东浦西区外,其他各分区ADI 指数均具有显著上升趋势,其中浙西区、太湖湖区、武澄锡虞区的TFPW-MK 值甚至超过95%的置信度,Hurst 系数分别为0.71,0.64 和0.65,均大于0.5,上升趋势的持续性较强。水利新闻0GBa{ V

水利新闻]} G4Ws c|jH

本文进一步对汛期各月降水量的变化趋势进行了检测,结果发现太湖流域5 和9月降水呈现显著下降趋势,汛期降水愈趋集中于6—8月。这是汛期降水不均匀性指数增强的主要原因。而这也说明汛期流域降水分布趋于更集中的态势,较容易导致汛期洪涝灾害的发生。水利新闻&{ Pqg,Bz-\

CCr[%u04 结语水利新闻:B-q ` XC'\!B

水利新闻y"X7{!T_!nrp

基于1951—2011 年太湖流域及各分区长系列降水资料,采用基尼系数、信息熵分配紊乱指数和变差系数分析了旬、月时间尺度上全年以及汛期降水时程分布的不均匀性。结果表明:ADI,GC 和Cv这3 种不均匀性指数的差异很小,相关系数可达0. 93 以上;太湖流域及各分区年内或汛期降水的不均匀性在年际间的差异较大,且北部武澄锡虞区和湖西区降水的年内分布不均匀性大于其他各分区;当降水量和指数值均较大时,流域较容易遭遇洪涝灾害,而在降水总量较小,不均匀性指数值较大时,则容易引发干旱灾害;在月尺度上,汛期降水不均匀性呈持续上升趋势,更加向6—8 月集中,太湖流域遭遇洪涝的风险将加大。

Ur*MU e~0

%~ U \ OsM5xe0此外,随着当前气候变化、人类活动的加剧以及流域城镇化进程的演进,区域性洪涝灾害愈加频繁。经分析表明武澄锡虞区和湖西区降水的年内分布不均匀性较大,凸显出降水空间分布的差异。因此,笔者认为有必要继续分析降水在空间上的不均匀性,以揭示降水在空间尺度上的分布差异,进一步为流域防灾减灾提供技术支撑。

I)Inh`:t.m0 水利新闻g![1] s` xx

作者简介:邓鹏鑫(1988-),男,福建建瓯人,硕士研究生,主要从事水文水资源方面的研究。水利新闻S\5}oP*B

顶:27 踩:41
对本篇资讯内容的质量打分:
当前平均分:-0.5 (152次打分)
【已经有98人表态】
15票
感动
16票
路过
11票
高兴
12票
难过
14票
搞笑
11票
愤怒
12票
无聊
7票
同情
下一篇:考虑土地利用类型的黑河中游水量平衡模型上一篇:淮河流域海水利用发展与对策

广告投放区


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


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

广告投放区

广告投放区