符 莎,成静清,刘章君
(1.曹娥江大闸运行管理中心,浙江 绍兴,312000;2.江西省水利科学院,江西 南昌,330029)
降水作为水文循环里的一个要素,影响着各个流域的径流、水位变化,在水文循环中起着至关重要的作用。近年来众多专家学者分析了降水量演变特征,其中王兆礼[1]等为阐明珠江流域降水量的时空演变特征,用Mann-Kendall法分析降水序列的变化趋势及其突变分量。结果表明:40年来,流域总降水量呈微弱的增加趋势,且具有明显的11年主周期振荡;东半部以增加为主,西南部以减少为主;1、2、3、7四个月份是增加明显的月份,其余各月都呈不显著减少趋势。刘贵花[2]等采用Mann-Kendall突变检验等统计分析方法分析流域降水、径流等水文气象要素及水库库容量变化趋势及突变特征,结果表明:流域年均降水逐渐上升,径流与降水年际变化一致,降水是引起径流年际变化的主要因素,而人类活动是引起径流季节变化的主要原因。曾小凡[3]等对长江流域年降水量进行了空间分区和时间演变分析,利用奇异谱分析方法来检测区域降水量发生突变的情况,结果表明,长江流域年降水量有比较明显的空间区域特征,各个区域年降水量的时间演变规律也不一致。衡彤[4]等利用小波变换对降水量时间序列的多时间尺度变化及突变特征进行了探讨,计算表明,其年际及年代际时间尺度在时域中分布不均匀,具有明显的局部化特征;主汛期降水与年降水的时间尺度变化比较接近。
鄱阳湖是我国最大的淡水湖,是一个过水性、季节性和吞吐型的浅淡水湖,位于江西省北部、长江中下游南岸。鄱阳湖纳赣江、抚河、信江、饶河和修水五水系以及博阳河、东河和西河诸河来水,经鄱阳湖调蓄后由湖口汇入长江。随着经济的快速增长,鄱阳湖的水文气象要素正逐年发生着巨大变化。由于气候变化和人类活动的影响不断加剧,鄱阳湖区水资源虽然较为丰富,但由于降雨时空分布不均,鄱阳湖水位变幅大,水资源开发利用程度低、难度大,自2003年三峡水库建设运行以来,长江对鄱阳湖水文水资源的影响也日益显现。因此,本文以鄱阳湖环湖区为研究对象,探究环湖区的降水演变特征。
本文根据全国水资源综合规划水资源分区划分,以鄱阳湖环湖平原区为研究对象,范围包括:赣江从外洲水文站以下、抚河李家渡水文站以下、信江梅港水文站以下、乐安河石镇街水文站以下、昌江古县渡水文站以下、修水永修水位站以下至湖口县的湖口水文站。
由于20世纪80年代之前进行围湖造田,80年代后的退耕还湖带来的变化以及2003年之后建造的三峡水库对环湖区影响变化,本文将对鄱阳湖环湖区9个雨量站以1980年和2003年为时间点分三个时间段对年内、年际变化进行分析,各雨量站分布如图1所示,表1列出详细的水文数据。
表1 雨量站资料汇总表
图1 鄱阳湖环湖区雨量站分布图
1.2.1 年内变化分析
采用年内分配不均匀系数(Cn)和完全调节系数(Cr)等指标作为衡量年内分配的标度,从不同角度定量分析年内分配特征的时空变化规律[5,6]。不均匀系数和完全调节系数是水文统计中用来说明降水分配不均匀性的重要参数,能够体现水文要素调控的难易程度,其值越大年内分配越不均匀。
式中:Ri为第i月平均降水量;为年内月平均降水量。
式中:R(i)为第i月平均降水量,R为月平均降水量,n为序列样本长度,下同。
1.2.2 年际变化分析
1.2.2.1 趋势性
衡量年际变化的特征值常用变差系数(Cv)和年际极值比(K)表示[5,6]。它们能够反映相对变化程度,其值越大,表示年际丰枯变化越剧烈。常用线性回归法和M-K趋势检验的趋势变化。
(1)线性回归法
如果水文时间序列中存在线性趋势,可采用线性相关方程进行描述,即线性回归检验,其数学模型为
采用最小二乘法对参数a、b进行估计,其值分别为
式中,xt为水文时间序列;n为序列样本长度;a、b为线性回归检验参数;t为时间;x为序列均值[7]。
(2)Mann-Kendall趋势检验
Mann-Kendall(M-K)非参数检验法常用于分析降水、径流、气温等要素时间序列的趋势变化,其优点在于受少数异常值干扰较少,计算简便[7-13]。
假设有 n个样本量(X1,…,Xn)的时间序列,对于所有 k,j≤n,且 k≠j,Xk和 Xj的分布是不同的,计算检验统计量S,公式如下:
S服从均值为0,方差为var(S)的正态分布。当n>10时,标准的正态统计变量Z通过下式计算:
当统计值Z大于0时,时间序列存在增加趋势;小于0时,则存在减少趋势。Z的绝对值大于1.28、1.96和2.32,分别表示通过了置信度90%、95%和99%的显著性检验。
1.2.2.2 突变性
Pettitt法是由Pettitt于1979年提出的一种非参数变点检验方法,该方法计算简便且受少数异常值干扰较小,在水文、气象等领域应用十分广泛。对于一个水文时间序列x1,x2,…,xn,假设序列中突变点最有可能发生在t时刻,因此,以t为分割点,将样本序列分为x1,x2,…,xt和 xt+1,xt+2,…,xn两个部分,分别服从分布 F1(x)和F2(x),Pettitt检验即在显著性水平α下,检验F1(x)和F2(x)是否为同一分布[14]。计算统计量Ut,n:
于是得到可能突变点位置t对应的统计量Kt为:
突变点的显著性水平为:
原假设为实测样本序列无突变,当pt>α时,则接受原假设,认为序列在位置t处不存在显著的突变点;当pt<α时,则拒绝原假设,认为序列在位置t处存在显著的突变点。
1.2.2.3 周期性
称ψa,b(t)为子小波;a为尺度因子,反映小波的周期长度;b为时间因子,反映时间上的平移。
小波变换若 ψa,b是(14)式给出的子小波,对于时间序列f(t)∈L2(R),其连续小波变换定义为:
式(14)中,ψ(t)为ψ(t)的复共轭函数;Wf(a,b)称小波(变换)系数。实际工作中,时间序列常常是离散的。如f(k△t)(k=1,2,…,N;△t为取样时间间隔),则式(15)的离散形式为:
Wf(a,b)随参数a和b变化,可作出以b为横坐标,a为纵坐标的关于Wf(a,b)的二维等值线图。将时间域上的关于a的所有小波系数的平方进行积分,即为小波方差,小波方差随尺度a变化过程称小波方差图[4,15]。
本文以9个雨量站的降水量的平均值作为环湖区的降水量,首先计算各雨量站1956~2016年月平均降水量占全年降水量的百分比,结果见表2,并绘制各雨量站的月平均分配的年内过程图(图2)。由表2可知各站连续4个月最大降雨量分布在4~7月。
图2 鄱阳湖环湖区雨量站1~12月降水量月平均分配年内过程图
表2 鄱阳湖环湖区雨量站1~12月月平均降水量占全年降水量百分比 %
从月平均降水量时间分布来看,各雨量站的最大月平均降水量均分布在6月,各站6月月平均降水量分别占全年降水量17.29%、18.87%、16.53%、16.25%、16.89%、14.73%、17.52%、17.23%、16.89%,庐山站6月平均降水量占年降水量的比重最小,进贤站6月降水量占年降水量的比重最大。各站月平均降水量最小的月份分布在12月,各站12月月平均降水量分别占全年降水量 3.00%、3.01%、2.93%、2.85%、3.02%、2.76%、2.89%、2.58%、2.71%,因此,鄱阳湖环湖区降水季节差异较大,年内分配极不均匀,这种变化加剧了降水年内分布的不平衡,增加了流域洪涝和干旱灾害发生的几率。
从月平均降水量的空间分布来看,各站连续最大4个月总降水量占年降水量分别是55.58%、56.45%、55.94% 、53.85% 、54.42% 、50.71% 、57.43% 、57.03% 、56.65%,其中石门街站最大,其次是星子站,环湖区的东北部较西南部降雨季节性差异更加明显;由图2可以看到,庐山站的各月平均降水量较其他站的降水量偏多,由于庐山站位于海拔1 164m的地方,受地形影响,不同高程降水量不同,庐山地区的海拔高程与降水量有较大关系,随高程增加而降水量也随之增大。
下面采用不均匀系数、完全调节系数的计算方法进一步对降水年内分配作定量分析,以1980年和2003年为界对比分析,计算成果见表3。从计算结果来看,昌邑站的不均匀系数从1956~1979的0.679上升到1980~2002的 0.734,再到 2003~2016年的 0.832,湖口站是逐步下降,其他站都是在1980年前后是下降的,2003年前后是上升的;而从完全调节系数来看,昌邑站和庐山站逐步上升,进贤站、湖口站、石门街以及星子站逐步下降,其余各站先下降后上升,各站降水量的不均匀系数和完全调节系数变化规律具有较好的同步性。
表3 鄱阳湖环湖区雨量站降水量年内分配特征值
从时间上来看,对比1980年以前,1980年以后除昌邑站降水量分配更加不均匀,其他站降水量的不均匀性有所缓解的趋势,但整体的变化不明显;对比2003年以后,各站降水量都较2003年以前不均匀性加剧,但不显著。例如都昌站在1956~1979年的不均匀系数和完全调节系数是 0.813、0.329,1980~2002年下降到0.763、0.316,2003~2016 年又上升到 0.847、0.330,而整个环湖区的降水量年内分配逐年趋向均匀。
从空间上来看,石门街站、星子站以及都昌站的不均匀系数和完全调节系数较其他雨量站更大,说明鄱阳湖环湖区西南部较东北部的降雨年内分配更均匀,同时环湖区东西部或者南北部的降水量相互消涨,使整个湖区降水量年内分配趋向均匀。
运用变差系数和年际极值比定量分析降水的相对变化程度,其值越大,说明降水的丰枯变化越剧烈。由表4可以看出庐山站的年最大降水量为3 034.8mm,发生在1975年,大多数雨量站的最大年降雨量发生在1998年,都昌站的最大年降水量最小,为699.1mm,发生在1978年,同时大多数雨量站的最小年降水量发生在1978年。各站的变差系数在0.19~0.45之间,其中庐山站的变差系数最大,各站的极值比在2.40~2.89之间,其中都昌站极值比最大。各站的年极值比和变差系数都较小,说明鄱阳湖环湖区的年降水量总体变化不显著。
表4 1956~2016年鄱阳湖环湖区雨量站降水量年际变化特征值
运用线性回归法和M-K趋势检验法对鄱阳湖环湖区的雨量站进行趋势分析,并将结果列于表5、表6,绘制了年降水量变化明显的代表站年际变化趋势图。表5中可以看到线性回归法求得的各雨量站的相关系数 R2,其中,超过 95%的置信水平用“*”表示。从图3中可以看,其中石门街站7月呈明显的上升趋势。而利用M-K趋势检验各站的各月平均降水量和年降水量,发现星子站和湖口站1月、石门街站的7月以及昌邑站的8月通过95%的显著性检验,其他站未发生明显趋势变化,说明鄱阳湖环湖区东北部降水量年际变化明显,并且丰枯变化较西南部更加显著。
图3 代表站降雨趋势变化图
表5 1956~2016年鄱阳湖环湖区雨量站1~12月年降水量线性回归分析结果
表6 鄱阳湖雨量站降水量M-K趋势检验结果
Pettitt检验在近几年被广泛用于水文序列的突变检验,对于变异点在两端的序列比较敏感,且适合检验大容量的样本序列,下面利用这一方法对雨量站进行突变检验,结果列于表7,并绘制出突变明显的雨量代表站(星子站、石门街站)前后均值的突变过程图,见图4,其中梓坊、湖口、庐山、石门街、星子等雨量站在1月降水量均发生显著突变,突变点都在1988年;石门街站7月降水量超过了95%的置信水平,突变发生在1995年。
图4 星子站及石门街站降水突变过程图
表7 鄱阳湖环湖区雨量站1~12月月平均降水量pettitt突变检验结果
根据鄱阳湖环湖区年降雨量进行周期分析,图5是鄱阳湖环湖区年降水量的实部(a)和模方(b)的时频分布图,模的大小表示特征时间尺度信号的强弱,实部表示不同特征时间尺度信号在不同时间上的分布和位相两方面的信息。从模方的时频分布图来看各时段尺度的强弱分布,其中3~9年时间尺度变化较强,主要发生在1960~1986年,振荡中心在1981年;22~28年的时间尺度变化也比较明显,主要出现在20世纪60年代到20世纪90年代,振荡中心在1981年左右。
图5 年降水量小波变换的实部(a)和模方(b)时频分布图
由实部等值线图可以看出年降水量时间尺度变化、突变点分布及其位相结构。其中3~7年尺度表现十分明显,其中心时间尺度为5年左右,正负位相交替出现;10~14年时间尺度表现突出,其尺度中心在13年左右,其次25~30年时间尺度表现十分明显,其尺度中心在28年左右。
根据图6的小波方差图可以明显地看出年降水量序列主要存在28年左右的主周期,图7绘制的28年时间尺度的小波变换系数变化过程,有四个多水期(1956~1964年、1974~1982年、1993~2002年以及 2011年以后)和三个少水期(1964~1974年、1982~1993年以及2002~2011年),旱涝突变点有1964年、1974年、1982年、1993年、2002年以及2011年,且鄱阳湖环湖区的年降水量的平均周期在19年左右,大约经历了3.5个变化期。
图6 年降水量小波方差图
图7 年降水量28年时间尺度的小波变换实部变化过程
通过对环湖区月尺度和年尺度的降水量分析结果,得出以下结论:
(1)环湖区降水量年内季节差异性大,主汛期占全年降雨的50%~60%,降水量年内分配极不均匀,12月降水量在全年最少,6月降水量全年最多,东北部较西南部降雨季节性差异更加明显,庐山站的降水量受地形影响较大。相较1956~1980年,1980~2003年区域降雨量年内分配趋向均匀,而2003~2016年,降雨量年内分配的不均匀性有所增长,从1956~2016年来看降水量的季节差异性有所下降。
(2)环湖区最大的年降水量在1998年,达2107.1mm,最小年降水量在1978年,为911.8mm,极值比为2.31,降水量的年际变化波动不大。1956~2016年星子站和湖口站1月降水量趋势变化明显,渡峰坑站8月以及石门街站7月的趋势变化明显;梓坊、湖口、庐山、石门街、星子等雨量站1月降水量均发生显著突变,突变点都在1988年,石门街站7月降水量超过了95%的置信水平,突变发生在1995年。而环湖区1~12月总体趋势以及突变并不明显,说明湖区内不同区域的降水量此消彼长,同时非汛期和汛期相互消涨,造成总体突变以及趋势变化不显著的现象。
(3)环湖区的年降水量序列存在28年左右的主周期,根据28年时间尺度的小波变换系数变化过程,在28年尺度下有4个多水期、3个少水期。
本文分析了环湖区近60年的降水变化,降水是影响径流变化的主要因素,当地降水量的变化会直接或间接影响洪水的变化,本文研究成果可为鄱阳湖洪旱灾害的预防治理和未来水资源调节防控提供借鉴和参考。