1. 引言
洪峰为一次洪水过程中测流断面上的最大流量。洪峰流量是洪水的重要特征之一,它不仅反映洪水的强度,也预示防洪的级别。洪峰流量也往往对应最高洪水位,是筑堤防洪参照量之一,洪峰流量的分析和预测对防洪建设至关重要 [1]。洪水灾害是影响经济发展最严重的自然灾害之一。博格达山北麓水系主要以冰川融雪和降雨补给为主,该区域水系肩负着阜康、吉木萨尔、奇台、木垒和巴里坤等县的工农业用水及人畜饮水的重任。且该区域常有春旱夏洪之灾,洪水灾害严重影响、制约着该区域农牧业及工业生产。因此,对博格达山北麓河区洪水变化趋势进行分析,对于预防洪水、减少旱灾有着十分重要的意义。
2. 研究区概况及数据
2.1. 研究区域
博格达山北麓河区是指发源于天山支脉博格达山北麓的诸多河流及湖泊,流经昌吉回族自治州的阜康市、吉木萨尔县、奇台县、木垒县及哈密地区的巴里坤县,其地理位置为87˚55'~91˚43'E,43˚30'~44˚35'N [2]。
Figure 1. Distribution diagram of rivers and lakes in the north of Bogda Mountain
图1. 博格达山北麓河区河湖分布示意图
整个流域地势南高北低,最高点为海拔5445 m的博格达峰,峰顶常年积雪,冰川总数213条,冰川面积91.5 km2。流域内自西向东依次有水磨沟河、三工河、天山天池、四工河、干河子、白杨河、黄山河、二工河、西大龙口河、东大龙口河、开垦河等诸多河流,且绝大多数河流水量,流程短,消失在北部的古尔班通古特沙漠南缘,博格达山北麓河区河湖分布示意图见图1。
2.2. 数据来源
本文选取博格达山北麓河区3处历史洪水资料较长的水文站进行分析。
3. 研究方法
3.1. Man-Kendall趋势检验
在Man-Kendall检验中 [1] [3] [4] [5],原假设H0为时间序列数据
,是n个独立的,随机变量同分布的样本;备择假设H1是双边检验。对于所有的i,j ≤ n,且i ≠ j,xi = xj的分布是不相同的。定义检验统计量S:
(1)
其中:sign()为符号函数。当xi − xj小于、等于或大于零时,sign(xi − xj)分别为−1、0或1。S为正态分布,其均值为0,方差
。Mann-Kendal统计了公式S大于、等于、小于零时分别为:
(2)
在双边趋势检验中,对于给定的置信水平α,若
,则原假设H0是不可接受的,即在置信水平α,时间序列数据存在明显的上升或下降的趋势。Z为正值表示增加趋势,负值表示减少趋势。Z的绝对值在大于等于1.28、1.64、2.32时表示分别通过了信度90%、95%、99%显著性检验。
3.2. Spearman秩相关法
许景璇,吴树诚等采用Spearman秩序相关检验法对年最大流量序列的趋势性特征进行分析 [6]。宋毅,施瑜等也采用该方法对年最大流量的趋势进行检验 [7]。该方法是一种非参数检验法,该方法计算简单,精确性较高。根据计算出的统计量T,选择显著性水平α,当
,表明序列趋势显著,反之
,则序列趋势不显著。
4. 研究结果与分析
4.1. 研究区洪峰流量量级及出现时段频次分析
经统计该区域3个水文站建站至今的实测洪峰流量:白杨河站59a系列中,超过警戒流量(70 m3/s)的洪峰流量出现8次,其中大于保证流量(160 m3/s)的1次;五圣宫站42a系列,超过警戒流量(70 m3/s)的洪峰流量出现2次,其中大于保证流量(100 m3/s)的2次;开垦河站65a系列中,超过警戒流量(130 m3/s)的洪峰流量出现7次,其中大于保证流量(250 m3/s)的3次。
年最大洪峰流量主要出现在6~7月(见表1)。从年最大洪峰流量出现的月份来看,阜康白杨河、奇台白杨河出现在7月的频次最高,开垦河出现在6月的频次最高。从年最大洪峰流量出现的时段来看,年最大洪峰流量主要出现在主汛期(6~7月)的60d时间内,阜康白杨河、奇台白杨河和开垦河最大洪峰流量出现频次分别占统计年数的90%、88%、68%。由此可见,以7月份为中心的前后60d,是该区出现大洪水的主汛期。
Table 1. The statistical characteristics of the annual maximum peak flow series in the northern Mountain Bogda Rivers
表1. 博格达山北麓河区各河流年最大洪峰流量序列统计特征值
4.2. 趋势分析
运用线性趋势法和5年滑动平均法对白杨河站、五圣宫站和开垦站建站以来年最大流量序列的趋势性进行初步分析,结果见图2。经过计算分析可知,1963~2021年,白杨河站年最大洪峰流量序列的均值为44.7 m3/s,最大值出现在2016年(320 m3/s),最小值出现在1963年(14.5 m3/s);1980~2021年,五圣宫站年最大洪峰流量序列的均值为41 m3/s,最大值出现在2007年(136 m3/s),最小值出现在2021年(16.2 m3/s);1957~2021年,开垦站年最大洪峰流量序列的均值为86.8 m3/s,最大值出现在1987年(357 m3/s),最小值出现在1968年(24.5 m3/s)。3站的年最大洪峰流量序列的线性系数分别为0.7891、0.1076、0.5596,均大于0,表明系列均呈上升趋势。
Figure 2. Trend line of maximum annual flow about the hydrological stations in the northern Mountain Bogda Rivers
图2. 博格达山北麓河区各水文站年最大洪峰流量趋势过程线
采用Mann-Kendall趋势检验法、Spearman秩相关法对博格达北麓河区白杨河站、五圣宫站和开垦河站年最大洪峰流量序列的趋势性进行趋势检验,取显著水平
,计算结果见表2。白杨河站年最大流量序列Mann-Kendall趋势检验法统计量Z的绝对值为2.55,大于
显著性水平临界值2.32;开垦河站年最大流量序列Mann-Kendall趋势检验法统计量Z的绝对值为2.12,大于
显著性水平临界值1.64;从而表明2站的年最大洪峰流量序列变化趋势显著;白杨河站年最大流量序列通过了Spearman秩相关法95%的显著性检验,其余2站均未通过90%的显著性检验,说明开垦河站、五圣宫站的年最大洪峰流量序列变化趋势不显著。
Table 2. Trend test results of annual maximum flow series about the hydrological stations in the district of northern Mountain Bogda Rivers
表2. 博格达山北麓河区各水文站年最大洪峰流量序列趋势检验结果
注:*,**和***分别表示通过了90%,95%和99%的显著性检验。
4.3. 突变分析
采用Mann-Kendall突变检验方法绘制UF和UB曲线图对该河区年最大洪峰流量序列进行图标分析(见图3)。由图2可知,白杨河的年最大洪峰流量在1996年之前呈缓慢上升的趋势,自1996年后呈明显上升的趋势且通过(α = 0.001)的显著性检验。五圣宫站1980年至2000年呈缓慢上升的趋势,2000年之后呈下降趋势,但趋势不明显,但在2004年存在突变。开垦河站1977年之前呈缓慢下降趋势,1975年至2007年呈现缓慢上升趋势,且趋势不明显。根据UF和UB曲线间交点的位置,开垦河年最大洪峰流量在1964、1983、1986、2008、2018年均存在突变。
Figure 3. The Mann-Kendall mutation test for maximum annual flow about the hydrological stations in the district of northern Mountain Bogda Rivers
图3. 博格达山北麓河区各水文站年最大洪峰流量Mann-Kendall突变检验
5. 结论
1) 从年最大洪峰流量出现的月份及频次来看,主汛期(6~7月)的60d时间内,阜康白杨河、奇台白杨河和开垦河最大洪峰流量出现频次很高,分别占统计年数的90%、88%、68%。由此可见,以7月份为中心的前后60d,是该区出现大洪水的主汛期;
2) 运用线性趋势法和5年滑动平均法对白杨河站、五圣宫站和开垦站年最大洪峰流量序列的趋势性进行初步分析,3站的年最大流洪峰量序列的线性系数分别为0.7891、0.1076、0.5596,均大于0,表明系列均呈上升趋势;
3) 采用Mann-Kendall趋势检验法,白杨河站、开垦河站的年最大洪峰流量序列通过了Mann-Kendal趋势检验法统计量Z的绝对值大于α = 0.05显著性水平临界值1.64、2.32,表明2站的年最大洪峰流量序列变化趋势显著;
4) 采用Spearman秩相关法,白杨河站年最大洪峰流量序列通过了Spearman秩相关法95%的显著性检验,表明该序列变化趋势显著;其余2站均未通过90%的显著性检验,说明这2站序列变化趋势不显著;
5) 采用Mann-Kendall突变检验方法,白杨河的年最大洪峰流量在1996年之前呈缓慢上升的趋势,自1996年后呈明显上升的趋势且通过(α = 0.001)显著性检验。五圣宫站1980年至2000年呈缓慢上升的趋势,自2000年之后呈下降趋势,但趋势不明显,但在2004年存在突变。开垦河站1977年之前呈缓慢下降趋势,1975年至2007年呈现缓慢上升趋势,且趋势不明显,但在1964、1983、1986、2008、2018年均存在突变。