第25卷第6期2009年11月源护水资保
wATER矗EsouRcEsPRoTEcnoN
V01.25
No.6
Nov.2009
邯郸地区年降水量统计特性
杜富慧1,郝振纯1,陈新美2,王晓丽3
(1.河海大学水文水资源与水利工程科学国家重点实验室,江苏南京210098;2.邯郸市水利局,河北邯郸056000;3.河北工程大学图书馆,河北邯郸056000)
摘要:运用线性回归分析与M帅n.Kendall趋势检验法对邯郸地区及各水资源分区的降雨资料进行分析,揭示了该地区及各水资源分区年降水量不断减少的变化趋势及显著性水平,并得出如下结论:各水资源分区年降水量变化规律不一致;随着海拔高度的降低,区域年降水量变化趋势越明显;低海拔区域年平均降水量减少是海河流域年径流量减少的重要原因。
关键词:年降水量;统计特性;线性回归分析;M柚n.Kendall检验;邯郸地区中图分类号:P333.1
文献标识码:A
文章编号:1004—6933(2009】06一0016—05
StatiSticalcharacteristicsof
DU
a皿ualprecipitation
inHandanregion
(1.&眦J时儿6D埘D叮矿跏切一耽ter胍DⅡ恻口以上砸h以幻啦i凡凹咖,舶胁i踟洫邝如,肫彬略210098,
c^in0;2.黜盯cD心脱,吖舰reⅡ“旷觑,渤n胁‘n幻咖口m,,,讹,渤056000,c^讹口;3.肌掘洗池邢蚵旷
precipitation
Fu-hIlil,HAO
Zhen.ch吼1。C腿N】【in.mei2,WANG
xiao-u3
啦i他e^,lg,讹n也056000。饥i胍)
Abstract:The肋nual
in山eHaIld柚mgionandeverywater阳sourc∞subregionw∞arIalyzed奶thlinear
te8t,
regres6ion肌alysis柚dtheMaIln-Kendaudecre够ing肌nual
卸dme
vaIiation
tendency蛐dBi印i6canceleveI“thecontinuously
was
p坨cipiⅢoninthatre舀on肌devery
water
water
re∞urcessubregionofwe陀陀vealed.Theconclusion
thatthevaJiouBpattemBofdifkr;entI.e∞urces8ubre舀ondifkred;thevariation0fre画onalprecipitationbecamemore
m08timporLant
at
apparent试thdecre鹪ingelevationaboves朗1evel;如dthe
HaihemVerw驰t11edecre曲eKe_ywords:aIlnual
rea∞nforthedecre舾ingdischargeofthe
0f咖ual
aVemgeprecipitation
loweleva“on8.
precipjtation;
statisticalcharacteri8tics;ljne缸regression蛐alysis;Mann-KendaU
test;
Handan
re西on
邯郸市位于华北平原中部,属于干旱半干旱地区,暖温带半湿润半干旱大陆性季风气候,全市多年平均降水量548.9mm(1956—2000年系列)。降水量时空分布不均,年际变化悬殊是其主要特征。该地区地形地貌特征差异明显,对华北地区而言具有代表性。地区年降水量的70%一80%集中在6—9月,其中又主要集中在7月下旬和8月上旬。该地区是我国水资源匮乏地区之一,人均水资源占有量只是全国平均水平的l/6,耕地单位面积水资源量只是全国的1/lO。而天然降水是水资源最重要的来源,降水是产生地表径流和补给地下水的主要来源,降
基金项目:国家自然科学基金(40830639;4080l叭2;5∞79016)
水量的大小及时空变化间接反映了一个地区的天然水资源状况。水资源与人们的生产和生活密切相关,受到广泛关注。因此,分析地区降水量的变化趋势,对研究地区水资源量的变化趋势具有参考价值,对区域水资源评价与再分配具有指导意义。笔者运用线性回归分析与Mann.Kendau趋势检验法对邯郸地区及各水资源分区的降雨资料进行了分析,得出了各水资源分区年降水量变化规律不一致的结论,并对同一地区不同水资源分区年降水量变化过程及趋势不一致的原因进行了分析。
作者简介:牡富慧(1974一)。男.山东济宁人.博士研究生.研究方向为水生态系统演化机理。B皿矗il:dufuhlli@hhu.edu.∞
・
16・
万方数据
1资料及方法
1.1资料
水文观测站点选择以观测降水站资料质量较好、系列较长、面上分布均匀且能反映地形变化影响等条件为原则。邯郸市境内共选出实测年数大于或等于34的分析雨量站37处,市境外选用雨量站10处。共选用雨量站47处。据统计,境内雨量站平均站网密度为325.6km2/站。其中:山区雨量站17处,
站网密度为262.4km2/站;平原雨量站20处。站网密度为379.4km2/站。
邯郸地区属太行山中南部中低山向河北平原西南部过渡地带,地形地貌复杂多变,形式多样,中低山、丘陵、盆地、平原和洼地均有分布,地形变化大。影响不同地形地貌特征的气象因素也存在较大差异,因此,各区域问的降水及水资源分布也有较大差异。
根据邯郸市的流域水系情况及地形的影响,将该地区的2个Ⅱ级流域分区重新划分为6个Ⅲ级流域分区,详细分区情况如图l所示。
1.2分析方法
笔者主要对邯郸地区降水量进行趋势分析。趋势分析研究的是时间序列顺序递增或递减的变化规律。这里仅介绍常用的线性回归方法、肯德尔秩相关检验法。其中线性回归方法属于参数方法,而肯德尔秩相关检验法属于非参数方法。
1.2.1线性回归分析
线性回归是一种参数估计方法[1引,设时间序列存在线性的趋势,即
毛=口+6t
(I=1,2,…,n)(1)
式中:毛为I时段的原始时间序列;n为时段数;口、6为使毛与口+缸误差最小的最小二乘估计量。可采用回归分析方法估计模型参数,即
“型可亍打
6=————■———■再广
宝(i一掣)(”牙)蚤(1.孚)。
Q)
Lz,
万
方数据a=i一占(孚)
(3)
式中二为随机变量的均值。为了检验线性趋势,可检验统计量ZL是否服从自由度为n一2的t分布。
z。=砉
(4)
若原假设为无趋势,选择置信水平口,计算统计量zI,。若lzLl>£口/2,拒绝原假设,从而认为线性趋势存在,否则接受原假设,线性趋势不存在。其中;6
指估计6的方差。可采用下式估计:
(5)
1.2.2肯德尔检验【3-5】
先计算统计量S:
^・l
^
s=∑∑sgn(巧一茗;)
(6)
‘-l』-‘+l
其中:
r
1(口>O)sgn(日)={o
(口=o)
(7)
L—l
(口<O)
不存在趋势的序列的统计量S满足正态分布,均值和方差分别为
E(s)=0
(8)
匕(s)=——————景L——一
n(n一1)(2n+5)一∑fJ(i—1)(2i+5)
(9)
最终可定义肯德尔检验统计量zK:
S
>
OZK=
赫。尚
S=
O,LlO
S
<
0zK满足标准正态分布【6]。利用zK的值进行趋势统计的显著性检验。zK值为正,表明上升趋势;zK值为负,表明下降趋势。在双尾趋势检验中,如果
zK
I≤zl一。/2,则接受零假设(无变化趋势);如果
zKI>z1.。/2,则拒绝零假设。zl一。/2从标准正态分布函数获得。口为检验的显著水平。
采用肯德尔检验时,时间序列的趋势幅度越大,或者序列长度越长,检验能力越强;而方差越小,则检验能力越弱。总的说来,对于没有趋势或者中等趋势的序列,如果序列长度较短,正相关使得误认为存在趋势的概率变大,而负相关使得误认为存在趋
势的概率变小;样本长度足够长,序列的自相关性基本不影响肯德尔趋势检验结果。
・
17
・
表l邯郸地区及各水资源分区年降雨系列线性检验殛Ma日昏Kendall检验计算结果
序列
区域名称序列年份长度/
.
线性趋势检验
zL
Mnn.Kend且ll检验
z‘
序列
区域d漳卫河平原1956。20∞马爝河平原1956.20∞
4S4545
线性趋势检验
置信水平
0.950.950.950.95
显著水平
称序列1日70、●。-……置信水平
一2.墙—2.15
M一.KendalI检验
zt
显著水平
O.os0,05
漳河山区
1蛳。2∞0
1956.20∞
45454545
一1.打
.1.52.0.94一1.58
一1.僻—0.∞—0.挣一1.∞
一2.潞
—2.40
0.95
滏阳河山区19耶。20∞滏西平原
邯郸地区1956。20∞
-1.弱一O.%
黑龙港平原1956—20∞
原
平贩
j粤
■
N
年份
围2邯郸地区及各水资源分区Z丘值变化过程
2.2滏阳河山区
2分析结果
邯郸地区及各水资源分区一元线性回归检验及M衄mKendau检验计算结果如表l所示,计算期内M蚰r卜Kencklu值(称为ZK值)的变化过程如图2所示。
2.1
滏阳河山区多年年平均降水量566.0栅,降水
量平均每lOa减少O.7rm,其趋势线方程见图4,线性趋势检验统计量ZL为一1.52,在95%的置信区间之内,线性变化趋势不明显。该区域内的M蚰n-Kendall检验统计量zK为一O.90,大于临界值一1.96,可以接受零假设。由z。值的变化过程可以看出,该
漳河山区
漳河山区多年平均年降水量576.9mm,降水量
地区肷值的变化过程与漳河山区非常接近,说明这
2个地区降水量分布具有较高的一致性。
平均每10a减少O.5咖,其趋势线方程见图3,线性
趋势检验统计量ZL为一1.27,在95%的置信区间之内(区间临界值为±2.O),因此,可以认为线性变化趋势不明显。该区域内的M柚n.Kendall检验统计量zK为一1.08,大于临界值一1.96,可以接受零假设。但由zx值的变化过程可以看出,该地区从1977年开始,年降水量减少明显,1980年超出了临界值,1980年后年降水量逐步增加。
暑E
删*进廿
枷咖
E
年份
图4滏阳河山区年降水量变化过程
{
蚓苌世廿
咖啪枷瑚
。
19561960
2.3滏西平原
滏西平原多年年平均降水量524.3姗,降水量平均每10a减少0.4舢,其趋势线方程见图5,线性
趋势检验统计量zL为一0.94,为各水资源分区中的最小值,远小于置信区间的临界值,线性变化趋势不
19“1968
19721976198019841988
1992
1996
20∞
年份
囝3漳河山区年降水量变化过程
明显。该区域内的M卸n.Kendall检验统计量zK为一0.29。也是各分区中的最小值,接受零假设,即无变化的假设。由ZK值的变化过程可以看出,该地区除1980年降水量减少明显外,其全过程较平稳。
・
18
・
万方数据
EE
由{
*溢廿
年份
图5滏西平原年降水量变化
2.4黑龙港平原
黑龙港平原多年年平均降水量535.4
m,降水
量平均每10a减少0.6衄,其趋势线方程见图6,线
性趋势检验统计量zL为一1.58,在置信区间之内,
线性变化趋势不明显。该区域内的M咖.Kendall检
验统计量zK为一1.08,在临界值以内,可以接受零假设。由zK值的变化过程可以看出,该地区1968年和1972年降水量减少明显,其他年份较稳定,无明显变化。
0
98
l
7面
6*5
进4
廿
3
2
1956
1960
19“196819721976198019841988199219962000
年份
图6黑龙港平原年降水量变化
2.5漳卫河平原
漳卫河平原多年年平均降水量548.5mm,降水
量平均每10a减少0.8帆,其趋势线方程见图7,线
性趋势检验统计量ZL为一2.18,为各分区中的最大
枷咖
EE
咖倒*鞋伽廿
枷瑚
。
19561960
19“196819721976198019841988199219962000
年份
围7漳卫河平原年降水量变化过程
万
方数据值,超出了95%的置信区问,因此,可以认为线性变
化趋势明显。该区域内的M蛐.Kendall检验统计量
ZK为一2.08,小于lI缶界值一1.96,因此,拒绝零假设,减少趋势明显。由z。值的变化过程可以看出,该地区从1977年开始,年降水量减少明显,1982年前后超出了临界值,1986年后年降水量进一步减少,并在较低的降水量值范围内波动。
2.6马颊河平原
马颊河平原多年年平均降水量547.7咖,降水
量平均每10a减少0.7mm,其趋势线方程见图8,线性趋势检验统计量ZL为一2.15,超出了置信区间,
线性变化趋势明显。该区域内的M咖.Kendau检验
统计量zK为一2.40,为各分区的最大值,小于I临界值一1.96,因此,拒绝零假设,呈趋势变化明显。由ZK值的变化过程可以看出,该地区从1977年开始,年降水量减少明显,1981年超出了临界值,1981年后年降水量均明显减少。
987
Eg
6豳{
*5
澄4
廿
3
2
l
∞∞∞∞∞∞∞∞∞0
19561960
19“19681972197619∞19¨19881992199620∞
年份
图8马颊河平原年降水量变化过程
3结论
a.邯郸地区各水资源分区年降水量均呈减少
趋势,整个地区平均每108减少0.6衄,线性趋势检
验统计量zL为一1.58,在置信区间之内,线性变化
趋势不明显。该区域内的M咖.Kendall检验统计量
ZK为一0.96,在临界值以内,即无明显变化,但各水资源分区的变化过程存在差异性,年降水量变化的空间分布不一致。
b.1980年前后年降水量的均值有明显减少。这与文献[7—8]中海河流域径流变化的趋势的结论相一致。其原因是黑龙港平原处于海河流域的上游,黑龙港平原的年降水量减少。这也从侧面揭示了海河流域径流减少的原因。
c.漳河山区、滏阳河山区、滏西平原、漳卫河平原等区域平均海拔都在100m以上,而黑龙港平原及马颊河平原的平均海拔在100m以下。这几个水
・
19
・
资源分区中,黑龙港平原及马颊河平原变化趋势明271.
显,超出了临界值,其他区域变化趋势不明显,可以[4]吴师.高嵩蝾.安徽省地表水水质变化趋势分析[J]。水
认为区域年降水量的变化趋势与区域的平均海拔高资源保护,2006,22(5):74—76.
度有关。[5]李俊峰,盛东,程晓如,等.玛纳斯河流域水质变化趋势
分析[J】.水资源保护,2008。24(4):10.13。
参考文献:
[6】刘攀,郭生练,肖义.等.水文时间序列趋势和跳跃分析
[1]丁晶,邓育仁.随机水文学[M].成都:成都科技大学出版
的再抽样方法研究[J].水文,200r7。27(2):49—53.
社.1988.
[7]刘春骧,刘志雨,谢正辉,等.近50年海河流域径流的变
[2]覃爱基,陈雪英,郑艳霞.宜昌径流时问序列的统计分析
化趋势研究[J].应用气象学报,2004.15(4):385。393.[J].水文,1993,13(5):15-21.[8]张建云,章四龙.王金星.等.近50年来中国六大流域年
[3】Y1厄S。Pr【DN
of际径流变化趋势研究[J].水科学进展,2007,18(2):230‘
P,CAVADIASG.POMertIIeMann-KendalI衄d
Spe咖’8
234.
rIlo
te马协fh
dctectillg哪∞toIIic№rlds
iⅡ
hydrologjcal鲥∞[J].Jo岫1且l
(收稿日期:20∞-10一15编辑:陈吉平)
of
Hyd∞lo盯,2002.259:254。
-●一+-+-+-・-・-・●—-—●--・●--+-+・・卜-・卜-・●・・・--・・●呻-+—+-+・—●-・-●・-・●--+-・●-・-●--+--●・-・●・-・●・-—●--+斗-—卜-—●-・・●・-—●—-●-・・。・.--—●-・+斗--_--・●-——●。・・■’-+-—._-‘●。‘。卜-+
(上接第ll页)
[9]贾玉霞.规模化畜禽养殖环境影响及主要防治问题[J].
[3]尚银生.吴有志,宋尚孝.水库水源保护区戈!f分方法初探
环境保护科学.2002,28(6):42_47.
[J].山西水利科技,1999(1):“146.
[10]赵虹,韩永升.沈阳市畜禽养殖环境污染分析及防治对
[4]王晓燕,王晓峰.汪清平,等.北京密云水库小流域非点
策[J].环境保护科学,2007,33(4):118一120.
源污染负荷估算[J].地理科学.2004,24(2):227—23I.【11]刘凌岩.王红军.大伙房水库上游禽畜业污染现状及其
【5]wlscHMEIER
wH,sMITH
DD.Rainfau
ene‘gy蚰di协
防治对策[J].环境保护科学,2002,28(3):21.22.他h山onBhipto∞iI
lo明[J].1hn距ctiom
of
Americ蛐
[12]李学灵。刘新媛.马岳雄,等.珠江流域片水资源保护规
G∞Ph8y8ic且I
Uni∞,1988,39(3):285—291.
划报告[R].广州:水利部珠江水利委员会,200l:65.[6]张力.为畜禽排污现状十分惊人,正成为中国新的污染
[13]洪晓瑜、逢勇.黎曼近似解模型在太湖藻类浓度场求解
大户破题:养殖业与种植业要“破镜重圆”[J].中国动物
中的应用[J].河海大学学报:自然科学版,2005,33(1):
保健.20cr7(9):5蛳1.
4l叫H.
[7]吴淑杭,姜震芳,俞清英.禽畜粪便污染现状与发展趋势
[14]韩菲,陈永灿,刘昭伟.湖泊及水库富营养化模型研究
[J].上海农业科技,2002(1):9-10.
综述[J].水科学进展,2003。14(6):785—791.
[8]苏杨.我国集约化畜禽养殖场污染月题研究[J].中国生
[15]张征,沈珍瑶,邵景力,等.环境评价学[M].北京:高等
态农业学报。2006,14(2):15一18.
教育出版社,2004:198.
(收稿日期:2008—1l-∞编辑:徐娟)
—+_..・・..——+——●—一.——■——+——卜・..-—■——■——’——+——+——卜——卜—■——■—・●——■—-■-・■--■_-+・■・-■—・■——●——●呻—■・・■・—●—・’——+——■—_.-・■-—一.——+・-■—・■——●——+——卜・_卜—..-—■—-●
(上接第15页)人流的处理。与其他河网地区河道[c]//黄真理.中国环境水力学2006.北京:中国水利水
相比,太浦河旁侧入流量很大。超过了东太湖来水水电出版社,2006:46-52。
量,如何处理旁侧人流对干流水环境容量计算结果[2]张红举.感潮河流水环境容量计算方法探讨[c]//房玲
的影响是较关键的问题。笔者提出的“控制污染负娣.水资源管理创新理论与实践.北京:中国水利水电出
荷量”的概念为解决这个问题提供了新的思路。③版社,2006:495-501.
河道不均匀混合系数的确定。由于采用分段零维公[3]程文辉,王船海。朱琰.太湖流域模型[M].南京:河海大
学出版社,2006:117一176.
式计算,不可避免地会带来由于污染物不均匀混合[4]孙卫红,姚国金,逢勇.基于不均匀系数的水环境容量计
造成的误差,因此,必须对结果进行修正。算方法探讨[J].水资源保护。200l,17(2):25—26.
参考文献:
(收稿日期:2008一12—30编辑:陈吉平)
[1]翟淑华。张红举.人工调控河道水环境容量计算研究
・20・
万
方数据
邯郸地区年降水量统计特性
作者:作者单位:
杜富慧, 郝振纯, 陈新美, 王晓丽, DU Fu-hui, HAO Zhen-chun, CHEN Xin-mei,WANG Xiao-li
杜富慧,郝振纯,DU Fu-hui,HAO Zhen-chun(河海大学水文水资源与水利工程科学国家重点实验室,江苏,南京,210098), 陈新美,CHEN Xin-mei(邯郸市水利局,河北,邯郸,056000), 王晓丽,WANG Xiao-li(河北工程大学图书馆,河北,邯郸,056000)水资源保护
WATER RESOURCES PROTECTION2009,25(6)
刊名:英文刊名:年,卷(期):
参考文献(8条)
1.丁晶;邓育仁 随机水文学 1988
2.张建云;章四龙;王金星 近50年来中国六大流域年际径流变化趋势研究[期刊论文]-水科学进展 2007(02)3.刘春臻;刘志雨;谢正辉 近50年海河流域径流的变化趋势研究[期刊论文]-应用气象学报 2004(04)4.刘攀;郭生练;肖义 水文时间序列趋势和跳跃分析的再抽样方法研究[期刊论文]-水文 2007(02)5.李俊峰;盛东;程晓如 玛纳斯河流域水质变化趋势分析[期刊论文]-水资源保护 2008(04)6.吴师;高嵩嵘 安徽省地表水水质变化趋势分析[期刊论文]-水资源保护 2006(05)
7.YUE S;PILON P;CAVADIAS G Power of the Mann-Kendall and Spearman' s rho tests for detectingmonotonic trends in hydrological series 2002
8.覃爱基;陈雪英;郑艳霞 宜昌径流时间序列的统计分析[期刊论文]-水文 1993(05)
本文读者也读过(10条)
1. 张国建.尹红美.李坷凌.ZHANG Guo-jian.YIN Hong-mei.LI Ke-ling 大气降水时空分布的数理统计与分析[期刊论文]-水资源保护2005,21(5)
2. 刘晓群.郝振纯 湖泊围垦对城陵矶水位抬升作用模拟分析[期刊论文]-安徽农业科学2009,37(26)3. 《水资源保护》稿约[期刊论文]-水资源保护2010,26(3)
4. 黄夏坤.郝振纯.梁才贵.朱长军.HUANG Xia-kun.HAO Zhen-chun.LIANG Cai-gui.ZHU Chang-jun 浔江防洪堤对西江中下游洪水的归槽影响[期刊论文]-桂林工学院学报2009,29(1)
5. 郝振纯.张晓鹏.张磊磊.王璐.史学丽.HAO Zhen-chun.ZHANG Xiao-peng.ZHANG Lei-lei.WANG Lu.SHI Xue-li 气候变暖下黄河源区冻土变化的数值模拟[期刊论文]-黑龙江水专学报2009,36(3)
6. 王建群.王尧.郝振纯 粒子群算法在小浪底水库补偿调度中的应用[期刊论文]-人民黄河2009,31(12)
7. 任政.郝振纯.REN Zheng.HAO Zhen-chun 滑动窗口二次自回归模型在径流预报中的应用[期刊论文]-河海大学学报(自然科学版)2009,37(3)
8. 郝振纯.谷军方.杜富慧.HAO Zhenchun.GU Junfang.DU Fuhui 随机GP法在概念性水文模型参数优选中的应用[期刊论文]-中山大学学报(自然科学版)2009,48(6)
9. 黄夏坤.郝振纯.罗惠先.贺小珍.HUANG Xia-kun.HAO Zhen-chun.LUO Hui-xian.HE Xiao-zheng 悬移质取样垂线流量权重法测定断面平均含沙量的应用研究[期刊论文]-泥沙研究2009(4)
10. 鞠琴.余钟波.郝振纯.欧耿鑫.佘超.刘德东.JU Qin.YU Zhong-bo.HAO Zhen-chun.SHE Chao.OU Geng-xin.LIUDe-dong 小波网络模型及其在日流量预测中的应用[期刊论文]-西安建筑科技大学学报(自然科学版)2009,41(1)
引证文献(1条)
1.杨肖丽.任立良.江善虎.袁山水.石琳 西辽河源头流域径流变化趋势及影响因素分析[期刊论文]-河海大学学报(自然科学版) 2012(1)
本文链接:http://d.g.wanfangdata.com.cn/Periodical_szybh200906005.aspx
第25卷第6期2009年11月源护水资保
wATER矗EsouRcEsPRoTEcnoN
V01.25
No.6
Nov.2009
邯郸地区年降水量统计特性
杜富慧1,郝振纯1,陈新美2,王晓丽3
(1.河海大学水文水资源与水利工程科学国家重点实验室,江苏南京210098;2.邯郸市水利局,河北邯郸056000;3.河北工程大学图书馆,河北邯郸056000)
摘要:运用线性回归分析与M帅n.Kendall趋势检验法对邯郸地区及各水资源分区的降雨资料进行分析,揭示了该地区及各水资源分区年降水量不断减少的变化趋势及显著性水平,并得出如下结论:各水资源分区年降水量变化规律不一致;随着海拔高度的降低,区域年降水量变化趋势越明显;低海拔区域年平均降水量减少是海河流域年径流量减少的重要原因。
关键词:年降水量;统计特性;线性回归分析;M柚n.Kendall检验;邯郸地区中图分类号:P333.1
文献标识码:A
文章编号:1004—6933(2009】06一0016—05
StatiSticalcharacteristicsof
DU
a皿ualprecipitation
inHandanregion
(1.&眦J时儿6D埘D叮矿跏切一耽ter胍DⅡ恻口以上砸h以幻啦i凡凹咖,舶胁i踟洫邝如,肫彬略210098,
c^in0;2.黜盯cD心脱,吖舰reⅡ“旷觑,渤n胁‘n幻咖口m,,,讹,渤056000,c^讹口;3.肌掘洗池邢蚵旷
precipitation
Fu-hIlil,HAO
Zhen.ch吼1。C腿N】【in.mei2,WANG
xiao-u3
啦i他e^,lg,讹n也056000。饥i胍)
Abstract:The肋nual
in山eHaIld柚mgionandeverywater阳sourc∞subregionw∞arIalyzed奶thlinear
te8t,
regres6ion肌alysis柚dtheMaIln-Kendaudecre够ing肌nual
卸dme
vaIiation
tendency蛐dBi印i6canceleveI“thecontinuously
was
p坨cipiⅢoninthatre舀on肌devery
water
water
re∞urcessubregionofwe陀陀vealed.Theconclusion
thatthevaJiouBpattemBofdifkr;entI.e∞urces8ubre舀ondifkred;thevariation0fre画onalprecipitationbecamemore
m08timporLant
at
apparent试thdecre鹪ingelevationaboves朗1evel;如dthe
HaihemVerw驰t11edecre曲eKe_ywords:aIlnual
rea∞nforthedecre舾ingdischargeofthe
0f咖ual
aVemgeprecipitation
loweleva“on8.
precipjtation;
statisticalcharacteri8tics;ljne缸regression蛐alysis;Mann-KendaU
test;
Handan
re西on
邯郸市位于华北平原中部,属于干旱半干旱地区,暖温带半湿润半干旱大陆性季风气候,全市多年平均降水量548.9mm(1956—2000年系列)。降水量时空分布不均,年际变化悬殊是其主要特征。该地区地形地貌特征差异明显,对华北地区而言具有代表性。地区年降水量的70%一80%集中在6—9月,其中又主要集中在7月下旬和8月上旬。该地区是我国水资源匮乏地区之一,人均水资源占有量只是全国平均水平的l/6,耕地单位面积水资源量只是全国的1/lO。而天然降水是水资源最重要的来源,降水是产生地表径流和补给地下水的主要来源,降
基金项目:国家自然科学基金(40830639;4080l叭2;5∞79016)
水量的大小及时空变化间接反映了一个地区的天然水资源状况。水资源与人们的生产和生活密切相关,受到广泛关注。因此,分析地区降水量的变化趋势,对研究地区水资源量的变化趋势具有参考价值,对区域水资源评价与再分配具有指导意义。笔者运用线性回归分析与Mann.Kendau趋势检验法对邯郸地区及各水资源分区的降雨资料进行了分析,得出了各水资源分区年降水量变化规律不一致的结论,并对同一地区不同水资源分区年降水量变化过程及趋势不一致的原因进行了分析。
作者简介:牡富慧(1974一)。男.山东济宁人.博士研究生.研究方向为水生态系统演化机理。B皿矗il:dufuhlli@hhu.edu.∞
・
16・
万方数据
1资料及方法
1.1资料
水文观测站点选择以观测降水站资料质量较好、系列较长、面上分布均匀且能反映地形变化影响等条件为原则。邯郸市境内共选出实测年数大于或等于34的分析雨量站37处,市境外选用雨量站10处。共选用雨量站47处。据统计,境内雨量站平均站网密度为325.6km2/站。其中:山区雨量站17处,
站网密度为262.4km2/站;平原雨量站20处。站网密度为379.4km2/站。
邯郸地区属太行山中南部中低山向河北平原西南部过渡地带,地形地貌复杂多变,形式多样,中低山、丘陵、盆地、平原和洼地均有分布,地形变化大。影响不同地形地貌特征的气象因素也存在较大差异,因此,各区域问的降水及水资源分布也有较大差异。
根据邯郸市的流域水系情况及地形的影响,将该地区的2个Ⅱ级流域分区重新划分为6个Ⅲ级流域分区,详细分区情况如图l所示。
1.2分析方法
笔者主要对邯郸地区降水量进行趋势分析。趋势分析研究的是时间序列顺序递增或递减的变化规律。这里仅介绍常用的线性回归方法、肯德尔秩相关检验法。其中线性回归方法属于参数方法,而肯德尔秩相关检验法属于非参数方法。
1.2.1线性回归分析
线性回归是一种参数估计方法[1引,设时间序列存在线性的趋势,即
毛=口+6t
(I=1,2,…,n)(1)
式中:毛为I时段的原始时间序列;n为时段数;口、6为使毛与口+缸误差最小的最小二乘估计量。可采用回归分析方法估计模型参数,即
“型可亍打
6=————■———■再广
宝(i一掣)(”牙)蚤(1.孚)。
Q)
Lz,
万
方数据a=i一占(孚)
(3)
式中二为随机变量的均值。为了检验线性趋势,可检验统计量ZL是否服从自由度为n一2的t分布。
z。=砉
(4)
若原假设为无趋势,选择置信水平口,计算统计量zI,。若lzLl>£口/2,拒绝原假设,从而认为线性趋势存在,否则接受原假设,线性趋势不存在。其中;6
指估计6的方差。可采用下式估计:
(5)
1.2.2肯德尔检验【3-5】
先计算统计量S:
^・l
^
s=∑∑sgn(巧一茗;)
(6)
‘-l』-‘+l
其中:
r
1(口>O)sgn(日)={o
(口=o)
(7)
L—l
(口<O)
不存在趋势的序列的统计量S满足正态分布,均值和方差分别为
E(s)=0
(8)
匕(s)=——————景L——一
n(n一1)(2n+5)一∑fJ(i—1)(2i+5)
(9)
最终可定义肯德尔检验统计量zK:
S
>
OZK=
赫。尚
S=
O,LlO
S
<
0zK满足标准正态分布【6]。利用zK的值进行趋势统计的显著性检验。zK值为正,表明上升趋势;zK值为负,表明下降趋势。在双尾趋势检验中,如果
zK
I≤zl一。/2,则接受零假设(无变化趋势);如果
zKI>z1.。/2,则拒绝零假设。zl一。/2从标准正态分布函数获得。口为检验的显著水平。
采用肯德尔检验时,时间序列的趋势幅度越大,或者序列长度越长,检验能力越强;而方差越小,则检验能力越弱。总的说来,对于没有趋势或者中等趋势的序列,如果序列长度较短,正相关使得误认为存在趋势的概率变大,而负相关使得误认为存在趋
势的概率变小;样本长度足够长,序列的自相关性基本不影响肯德尔趋势检验结果。
・
17
・
表l邯郸地区及各水资源分区年降雨系列线性检验殛Ma日昏Kendall检验计算结果
序列
区域名称序列年份长度/
.
线性趋势检验
zL
Mnn.Kend且ll检验
z‘
序列
区域d漳卫河平原1956。20∞马爝河平原1956.20∞
4S4545
线性趋势检验
置信水平
0.950.950.950.95
显著水平
称序列1日70、●。-……置信水平
一2.墙—2.15
M一.KendalI检验
zt
显著水平
O.os0,05
漳河山区
1蛳。2∞0
1956.20∞
45454545
一1.打
.1.52.0.94一1.58
一1.僻—0.∞—0.挣一1.∞
一2.潞
—2.40
0.95
滏阳河山区19耶。20∞滏西平原
邯郸地区1956。20∞
-1.弱一O.%
黑龙港平原1956—20∞
原
平贩
j粤
■
N
年份
围2邯郸地区及各水资源分区Z丘值变化过程
2.2滏阳河山区
2分析结果
邯郸地区及各水资源分区一元线性回归检验及M衄mKendau检验计算结果如表l所示,计算期内M蚰r卜Kencklu值(称为ZK值)的变化过程如图2所示。
2.1
滏阳河山区多年年平均降水量566.0栅,降水
量平均每lOa减少O.7rm,其趋势线方程见图4,线性趋势检验统计量ZL为一1.52,在95%的置信区间之内,线性变化趋势不明显。该区域内的M蚰n-Kendall检验统计量zK为一O.90,大于临界值一1.96,可以接受零假设。由z。值的变化过程可以看出,该
漳河山区
漳河山区多年平均年降水量576.9mm,降水量
地区肷值的变化过程与漳河山区非常接近,说明这
2个地区降水量分布具有较高的一致性。
平均每10a减少O.5咖,其趋势线方程见图3,线性
趋势检验统计量ZL为一1.27,在95%的置信区间之内(区间临界值为±2.O),因此,可以认为线性变化趋势不明显。该区域内的M柚n.Kendall检验统计量zK为一1.08,大于临界值一1.96,可以接受零假设。但由zx值的变化过程可以看出,该地区从1977年开始,年降水量减少明显,1980年超出了临界值,1980年后年降水量逐步增加。
暑E
删*进廿
枷咖
E
年份
图4滏阳河山区年降水量变化过程
{
蚓苌世廿
咖啪枷瑚
。
19561960
2.3滏西平原
滏西平原多年年平均降水量524.3姗,降水量平均每10a减少0.4舢,其趋势线方程见图5,线性
趋势检验统计量zL为一0.94,为各水资源分区中的最小值,远小于置信区间的临界值,线性变化趋势不
19“1968
19721976198019841988
1992
1996
20∞
年份
囝3漳河山区年降水量变化过程
明显。该区域内的M卸n.Kendall检验统计量zK为一0.29。也是各分区中的最小值,接受零假设,即无变化的假设。由ZK值的变化过程可以看出,该地区除1980年降水量减少明显外,其全过程较平稳。
・
18
・
万方数据
EE
由{
*溢廿
年份
图5滏西平原年降水量变化
2.4黑龙港平原
黑龙港平原多年年平均降水量535.4
m,降水
量平均每10a减少0.6衄,其趋势线方程见图6,线
性趋势检验统计量zL为一1.58,在置信区间之内,
线性变化趋势不明显。该区域内的M咖.Kendall检
验统计量zK为一1.08,在临界值以内,可以接受零假设。由zK值的变化过程可以看出,该地区1968年和1972年降水量减少明显,其他年份较稳定,无明显变化。
0
98
l
7面
6*5
进4
廿
3
2
1956
1960
19“196819721976198019841988199219962000
年份
图6黑龙港平原年降水量变化
2.5漳卫河平原
漳卫河平原多年年平均降水量548.5mm,降水
量平均每10a减少0.8帆,其趋势线方程见图7,线
性趋势检验统计量ZL为一2.18,为各分区中的最大
枷咖
EE
咖倒*鞋伽廿
枷瑚
。
19561960
19“196819721976198019841988199219962000
年份
围7漳卫河平原年降水量变化过程
万
方数据值,超出了95%的置信区问,因此,可以认为线性变
化趋势明显。该区域内的M蛐.Kendall检验统计量
ZK为一2.08,小于lI缶界值一1.96,因此,拒绝零假设,减少趋势明显。由z。值的变化过程可以看出,该地区从1977年开始,年降水量减少明显,1982年前后超出了临界值,1986年后年降水量进一步减少,并在较低的降水量值范围内波动。
2.6马颊河平原
马颊河平原多年年平均降水量547.7咖,降水
量平均每10a减少0.7mm,其趋势线方程见图8,线性趋势检验统计量ZL为一2.15,超出了置信区间,
线性变化趋势明显。该区域内的M咖.Kendau检验
统计量zK为一2.40,为各分区的最大值,小于I临界值一1.96,因此,拒绝零假设,呈趋势变化明显。由ZK值的变化过程可以看出,该地区从1977年开始,年降水量减少明显,1981年超出了临界值,1981年后年降水量均明显减少。
987
Eg
6豳{
*5
澄4
廿
3
2
l
∞∞∞∞∞∞∞∞∞0
19561960
19“19681972197619∞19¨19881992199620∞
年份
图8马颊河平原年降水量变化过程
3结论
a.邯郸地区各水资源分区年降水量均呈减少
趋势,整个地区平均每108减少0.6衄,线性趋势检
验统计量zL为一1.58,在置信区间之内,线性变化
趋势不明显。该区域内的M咖.Kendall检验统计量
ZK为一0.96,在临界值以内,即无明显变化,但各水资源分区的变化过程存在差异性,年降水量变化的空间分布不一致。
b.1980年前后年降水量的均值有明显减少。这与文献[7—8]中海河流域径流变化的趋势的结论相一致。其原因是黑龙港平原处于海河流域的上游,黑龙港平原的年降水量减少。这也从侧面揭示了海河流域径流减少的原因。
c.漳河山区、滏阳河山区、滏西平原、漳卫河平原等区域平均海拔都在100m以上,而黑龙港平原及马颊河平原的平均海拔在100m以下。这几个水
・
19
・
资源分区中,黑龙港平原及马颊河平原变化趋势明271.
显,超出了临界值,其他区域变化趋势不明显,可以[4]吴师.高嵩蝾.安徽省地表水水质变化趋势分析[J]。水
认为区域年降水量的变化趋势与区域的平均海拔高资源保护,2006,22(5):74—76.
度有关。[5]李俊峰,盛东,程晓如,等.玛纳斯河流域水质变化趋势
分析[J】.水资源保护,2008。24(4):10.13。
参考文献:
[6】刘攀,郭生练,肖义.等.水文时间序列趋势和跳跃分析
[1]丁晶,邓育仁.随机水文学[M].成都:成都科技大学出版
的再抽样方法研究[J].水文,200r7。27(2):49—53.
社.1988.
[7]刘春骧,刘志雨,谢正辉,等.近50年海河流域径流的变
[2]覃爱基,陈雪英,郑艳霞.宜昌径流时问序列的统计分析
化趋势研究[J].应用气象学报,2004.15(4):385。393.[J].水文,1993,13(5):15-21.[8]张建云,章四龙.王金星.等.近50年来中国六大流域年
[3】Y1厄S。Pr【DN
of际径流变化趋势研究[J].水科学进展,2007,18(2):230‘
P,CAVADIASG.POMertIIeMann-KendalI衄d
Spe咖’8
234.
rIlo
te马协fh
dctectillg哪∞toIIic№rlds
iⅡ
hydrologjcal鲥∞[J].Jo岫1且l
(收稿日期:20∞-10一15编辑:陈吉平)
of
Hyd∞lo盯,2002.259:254。
-●一+-+-+-・-・-・●—-—●--・●--+-+・・卜-・卜-・●・・・--・・●呻-+—+-+・—●-・-●・-・●--+-・●-・-●--+--●・-・●・-・●・-—●--+斗-—卜-—●-・・●・-—●—-●-・・。・.--—●-・+斗--_--・●-——●。・・■’-+-—._-‘●。‘。卜-+
(上接第ll页)
[9]贾玉霞.规模化畜禽养殖环境影响及主要防治问题[J].
[3]尚银生.吴有志,宋尚孝.水库水源保护区戈!f分方法初探
环境保护科学.2002,28(6):42_47.
[J].山西水利科技,1999(1):“146.
[10]赵虹,韩永升.沈阳市畜禽养殖环境污染分析及防治对
[4]王晓燕,王晓峰.汪清平,等.北京密云水库小流域非点
策[J].环境保护科学,2007,33(4):118一120.
源污染负荷估算[J].地理科学.2004,24(2):227—23I.【11]刘凌岩.王红军.大伙房水库上游禽畜业污染现状及其
【5]wlscHMEIER
wH,sMITH
DD.Rainfau
ene‘gy蚰di协
防治对策[J].环境保护科学,2002,28(3):21.22.他h山onBhipto∞iI
lo明[J].1hn距ctiom
of
Americ蛐
[12]李学灵。刘新媛.马岳雄,等.珠江流域片水资源保护规
G∞Ph8y8ic且I
Uni∞,1988,39(3):285—291.
划报告[R].广州:水利部珠江水利委员会,200l:65.[6]张力.为畜禽排污现状十分惊人,正成为中国新的污染
[13]洪晓瑜、逢勇.黎曼近似解模型在太湖藻类浓度场求解
大户破题:养殖业与种植业要“破镜重圆”[J].中国动物
中的应用[J].河海大学学报:自然科学版,2005,33(1):
保健.20cr7(9):5蛳1.
4l叫H.
[7]吴淑杭,姜震芳,俞清英.禽畜粪便污染现状与发展趋势
[14]韩菲,陈永灿,刘昭伟.湖泊及水库富营养化模型研究
[J].上海农业科技,2002(1):9-10.
综述[J].水科学进展,2003。14(6):785—791.
[8]苏杨.我国集约化畜禽养殖场污染月题研究[J].中国生
[15]张征,沈珍瑶,邵景力,等.环境评价学[M].北京:高等
态农业学报。2006,14(2):15一18.
教育出版社,2004:198.
(收稿日期:2008—1l-∞编辑:徐娟)
—+_..・・..——+——●—一.——■——+——卜・..-—■——■——’——+——+——卜——卜—■——■—・●——■—-■-・■--■_-+・■・-■—・■——●——●呻—■・・■・—●—・’——+——■—_.-・■-—一.——+・-■—・■——●——+——卜・_卜—..-—■—-●
(上接第15页)人流的处理。与其他河网地区河道[c]//黄真理.中国环境水力学2006.北京:中国水利水
相比,太浦河旁侧入流量很大。超过了东太湖来水水电出版社,2006:46-52。
量,如何处理旁侧人流对干流水环境容量计算结果[2]张红举.感潮河流水环境容量计算方法探讨[c]//房玲
的影响是较关键的问题。笔者提出的“控制污染负娣.水资源管理创新理论与实践.北京:中国水利水电出
荷量”的概念为解决这个问题提供了新的思路。③版社,2006:495-501.
河道不均匀混合系数的确定。由于采用分段零维公[3]程文辉,王船海。朱琰.太湖流域模型[M].南京:河海大
学出版社,2006:117一176.
式计算,不可避免地会带来由于污染物不均匀混合[4]孙卫红,姚国金,逢勇.基于不均匀系数的水环境容量计
造成的误差,因此,必须对结果进行修正。算方法探讨[J].水资源保护。200l,17(2):25—26.
参考文献:
(收稿日期:2008一12—30编辑:陈吉平)
[1]翟淑华。张红举.人工调控河道水环境容量计算研究
・20・
万
方数据
邯郸地区年降水量统计特性
作者:作者单位:
杜富慧, 郝振纯, 陈新美, 王晓丽, DU Fu-hui, HAO Zhen-chun, CHEN Xin-mei,WANG Xiao-li
杜富慧,郝振纯,DU Fu-hui,HAO Zhen-chun(河海大学水文水资源与水利工程科学国家重点实验室,江苏,南京,210098), 陈新美,CHEN Xin-mei(邯郸市水利局,河北,邯郸,056000), 王晓丽,WANG Xiao-li(河北工程大学图书馆,河北,邯郸,056000)水资源保护
WATER RESOURCES PROTECTION2009,25(6)
刊名:英文刊名:年,卷(期):
参考文献(8条)
1.丁晶;邓育仁 随机水文学 1988
2.张建云;章四龙;王金星 近50年来中国六大流域年际径流变化趋势研究[期刊论文]-水科学进展 2007(02)3.刘春臻;刘志雨;谢正辉 近50年海河流域径流的变化趋势研究[期刊论文]-应用气象学报 2004(04)4.刘攀;郭生练;肖义 水文时间序列趋势和跳跃分析的再抽样方法研究[期刊论文]-水文 2007(02)5.李俊峰;盛东;程晓如 玛纳斯河流域水质变化趋势分析[期刊论文]-水资源保护 2008(04)6.吴师;高嵩嵘 安徽省地表水水质变化趋势分析[期刊论文]-水资源保护 2006(05)
7.YUE S;PILON P;CAVADIAS G Power of the Mann-Kendall and Spearman' s rho tests for detectingmonotonic trends in hydrological series 2002
8.覃爱基;陈雪英;郑艳霞 宜昌径流时间序列的统计分析[期刊论文]-水文 1993(05)
本文读者也读过(10条)
1. 张国建.尹红美.李坷凌.ZHANG Guo-jian.YIN Hong-mei.LI Ke-ling 大气降水时空分布的数理统计与分析[期刊论文]-水资源保护2005,21(5)
2. 刘晓群.郝振纯 湖泊围垦对城陵矶水位抬升作用模拟分析[期刊论文]-安徽农业科学2009,37(26)3. 《水资源保护》稿约[期刊论文]-水资源保护2010,26(3)
4. 黄夏坤.郝振纯.梁才贵.朱长军.HUANG Xia-kun.HAO Zhen-chun.LIANG Cai-gui.ZHU Chang-jun 浔江防洪堤对西江中下游洪水的归槽影响[期刊论文]-桂林工学院学报2009,29(1)
5. 郝振纯.张晓鹏.张磊磊.王璐.史学丽.HAO Zhen-chun.ZHANG Xiao-peng.ZHANG Lei-lei.WANG Lu.SHI Xue-li 气候变暖下黄河源区冻土变化的数值模拟[期刊论文]-黑龙江水专学报2009,36(3)
6. 王建群.王尧.郝振纯 粒子群算法在小浪底水库补偿调度中的应用[期刊论文]-人民黄河2009,31(12)
7. 任政.郝振纯.REN Zheng.HAO Zhen-chun 滑动窗口二次自回归模型在径流预报中的应用[期刊论文]-河海大学学报(自然科学版)2009,37(3)
8. 郝振纯.谷军方.杜富慧.HAO Zhenchun.GU Junfang.DU Fuhui 随机GP法在概念性水文模型参数优选中的应用[期刊论文]-中山大学学报(自然科学版)2009,48(6)
9. 黄夏坤.郝振纯.罗惠先.贺小珍.HUANG Xia-kun.HAO Zhen-chun.LUO Hui-xian.HE Xiao-zheng 悬移质取样垂线流量权重法测定断面平均含沙量的应用研究[期刊论文]-泥沙研究2009(4)
10. 鞠琴.余钟波.郝振纯.欧耿鑫.佘超.刘德东.JU Qin.YU Zhong-bo.HAO Zhen-chun.SHE Chao.OU Geng-xin.LIUDe-dong 小波网络模型及其在日流量预测中的应用[期刊论文]-西安建筑科技大学学报(自然科学版)2009,41(1)
引证文献(1条)
1.杨肖丽.任立良.江善虎.袁山水.石琳 西辽河源头流域径流变化趋势及影响因素分析[期刊论文]-河海大学学报(自然科学版) 2012(1)
本文链接:http://d.g.wanfangdata.com.cn/Periodical_szybh200906005.aspx