闽江流域位于东经116˚23'~119˚43'，北纬25˚23'~28˚19'之间，属于亚热带海洋季风气候区，流域内的植被较好，气候温和，雨量充沛，是福建省最大的河流，发源于武夷山脉，流经36个县(市)，全长541 km，流域面积60,992 km²。其中，福建省境内面积59,922 km²，占流域面积的98.2%，约占福建省面积的一半；属浙江省境内面积1070 km²，占流域面积的1.8%。

闽江流域状如扇形，地势西北高东南低，南平以上有建溪、富屯溪和沙溪三大支流发源于武夷山脉，平均海拔较高，地形较为复杂，险滩多，水流湍急，流态紊乱，具有水位暴涨暴落的山区河流的特点。南平以下为闽江干流，自南平到闽江入海口，沿途有吉溪、尤溪、古田溪、梅溪、大樟溪等支流汇入。大樟溪发源于德化县戴云山主峰西面，是闽江下游最大的支流，地形同样复杂，水流落差较大，具有山区性河流特点。闽江干流下游逐渐从丘陵区过渡到福州平原，海拔较三大支流明显较低，水面落差较小。闽江下游在淮安分南北两港，北港绕南台岛北面，穿过福州市区至马尾与绕过南台岛南面的南港相汇，再流经闽安镇被琅岐岛分为南北两支，南支经长乐县梅花镇，北支经连江县琯头镇流入北海(图1)。

据调查，20世纪50年代后期至今，闽江流域中上游的干流和支流修建了大型水库7座，中型水库60多处和一些小型水库 [ 15 ]，其中主要的大型水库的水利工程建成情况见表1。至今这些大中型水库相继建成蓄水，主要用于发电、防洪与灌溉 [ 16 ]，重点用于调节年内径流，对径流的年际变化影响可能较小 [ 15 ]。大型水库对洪水有一定的自然调蓄功能，此外为了下游的防洪的需要，还需进行防洪错峰调度，以降低下游的洪峰流量。大型水库对径流的年内分配也产生了一定的影响，它们在一定程度上加大了枯水季节的径流量。使下游河道枯水期日平均流量和月平均流量相对增加，径流的年内分配较均匀，从而降低了径流的年内的波动性，FI值也相应的降低。由于闽江流域范围广，年径流量大，支流分布多，水电开发的年限较长，因此水电的开发对径流的年内分布作用并不明显，但同样对降低径流的年内波动性起到了一定的影响。另外，闽江干支流绝大多数梯级电站建成于1993年，特别是2000年以后，目前流域梯级电站开发率达到85%，水口水库以上流域主要一级支流干流上已建梯级水电站24座，二、三级支流上的梯级小水电站有数百座 [ 17 ]。水电站建设过程中产生了大量的泥沙倒入河流中，建成后对河流的泥沙起到了拦截作用。

表1. 闽江流域主要水利工程统计 [ 15 ]

注：*古田溪水库分为4个阶段建设，第一阶段建成是1959年，最终完成是在1971年。

图1. 闽江流域研究区域概况图

1) Flashiness指数。水文资料选取闽江流域的干流(竹岐)、沙溪(沙县石桥)、富屯溪(洋口)、大樟溪(永泰清水)、建溪(七里街)水文站近60年来的日径流量监测数据。基于Baker等人 [ 7 ] 提出的Flashiness指数分析闽江流域径流量的时空变化，其计算公式如下：

FI = ∑ i = 1 n | q i − q i − 1 | ∑ i = 1 n q i (1)

式中：FI为Flashiness指数；q为日均径流量；i为天数，n = 365 (366)。FI是一个无量纲的指数值，介于0~2之间 [ 12 ]。FI趋于0表示流域的径流波动越小，FI趋于2表示流域的日径流量波动就越大。

2) 基流分离及基流指数(BFI)的计算。本研究采用数字滤波法进行基流的分割，滤波法较传统的图解法具有更加客观，操作简便等优点 [ 18 ]。本研究基于Lyne-Hollick算法进行基流的分离及基流指数的计算。

3) Mann-Kendall趋势检验。近年来Mann-Kendall检验被广泛应用于FI值、BFI值的变化趋势检验 [ 19 ]。Kendall τ相关系数用非参数检验的方法来度量定序变量间的线性相关关系。

Kendall τ相关系数的检验统计量计算公式如下 [ 18 ] ：

τ = 2 ( U − V ) n ( n − 1 ) (2)

式中：τ为Mann-Kendall相关系数；U是一致对数目；V为非一致对数目；n为样本数。

当U值较大、V值较小时，两变量呈较强的正相关；当U值较小，V值较大时，两变量呈较强的负相关；如果U值和V值相近时，两变量呈较弱的相关关系。

闽江干流、四条支流的Flashiness指数发生了明显的变化。近60年Flashiness指数经Mann-Kendall分析可见(表2)，其τ值均小于0，大樟溪流域和建溪流域的P值小于0.01，富屯溪的p值小于0.05。这说明近60a来闽江流域的年内径流波动性呈现出下降的趋势，其中富屯溪显著下降，大樟溪和建溪呈现出极显著的下降趋势，沙溪和闽江干流下降趋势不明显(图2)。总体而言，大樟溪和建溪的径流波动性更大，表现在FI值比其他流域更大(见表3)。

表2. 闽江干流、四条支流近60年FI值趋势检验结果

表3. 闽江干流及四条支流近60年FI值统计特征

图2. 1950~2012年闽江干流及支流Flashiness指数

近60a来，闽江干流、沙溪、大樟溪、富屯溪、建溪的基流指数(BFI)也发生了明显的变化(图3)。经Mann-Kendall分析(表4)，其τ值分别为0.2135、0.3319、0.5234、0.2367、0.3141且绝对值都在0.2以上，闽江干流p值小于0.05，其他支流的p值均小于0.01，说明近60a来闽江干流基流指数呈显著的上升趋势，其他支流的基流指数呈极显著的上升趋势。由图3和表4可见，闽江流域的基流指数呈明显的时空变异特征，即闽江干流、沙溪、富屯溪的基流指数明显高于建溪和大樟溪，且都呈现出显著的上升趋势。说明影响闽江干流和四条支流河流径流量的水文地理条件发生了明显的年际变化。

基于闽江流域1950~2012年近60年的日径流量数据由2.2计算得出近60a闽江流域Flashiness指数，并对FI指数进行变异点分析，结果见图4。

图3. 1950~2012年闽江干流与四大支流BFI指数

表4. 1950~2012年闽江干流、四大支流BFI值的M-K检验结果

图4. 闽江干流及支流FI指数变异点分析

由图4可见，闽江干流及四条支流的Flashiness指数仅在建溪有出现变异点，时间是1968年和1970年。据调查，闽江流域大坝建设、水电开发开始于50年代至今，其干支流绝大多数梯级电站建成于1993年，特别是2000年以后，目前流域梯级电站开发率达到85%，但是根据各流域变异点分析结果显示，只有建溪流域的FI值分别在1968年和1970年出现拐点，其他流域FI指数在近60年甚至没有出现拐点。这说明闽江流域径流的波动性受水电开发建设的影响并不敏感。闽江是福建省最大的河流，流域面积约占福建省陆域总面积的一半，流域尺度较大，水库分布在流域的各干支流，这可能是造成水电开发建设对径流波动性变化影响不明显的原因。另一方面，闽江流域水电开发的持续较长，也可能是闽江径流波动性变化不明显的原因。从50年代中后期起大型水库开始兴修水利，还有60座中型水库也是这段时间开始兴修。

本文基于竹岐(闽江干流)、沙县石桥(沙溪)、洋口(富屯溪)、永泰清水(大樟溪)、七里街(建溪)五个水文站点收集的1978~2012年近35年的日输沙率数据进行变异点分析，结果见图5。

图5. 闽江干流及四条支流近35年输沙率变异点

由图5可见，闽江干流、沙溪、富屯溪、大樟溪和建溪的输沙率分别在1984、1993 (2000)、1996 (1998)、1991 (2006)和1997年出现变异点。除闽江干流竹岐水文站点输沙率拐点时间为80年代，其余四条支流所对应的水文站点的输沙率拐点均分布在90年代至2000以后，即处于闽江流域梯级电站开发的高峰期。因此对1978~1992年和1993~2012年两个时间段的年平均输沙率进行对比(见表5)。由表5可知，1993年前后各水文站点所测得的年平均输沙率均发生了一定的变化，这反应了水电站的梯级开发使得闽江流域输沙率发生了明显变化，从而出现了变异点。闽江干支流绝大多数梯级电站建成于1993年，特别是2000年以后，目前流域梯级电站开发率达到85%，大规模的梯级电站建设过程中，大量的泥沙排入到流域中，使得河流的输沙量增加。而这些水库的建成后，如同一个个巨型沉沙池，沉积了上游的大量的泥沙，使得下泄的河水的含量明显减少。

表5. 闽江干流及支流水文站点1993年前后年均输沙率变化/kg/s