赵贵章，徐远志,王莉莉，王艳艳，张世英

(华北水利水电大学地球科学与工程学院，河南 郑州 450045)

基流是黄河水资源的重要组成部分，它不仅能够反映流域地下水对河流的补给情况，同时在提高水资源利用率、保护生态环境以及计算生态需水量、河流输沙量等方面具有重要意义。 20世纪60年代以来，由于修建水利枢纽等一系列人类活动因素的影响[1]，黄河径流量大幅度减少[2]，改变了原有的地下水的时空分布规律，黄河基流量呈现明显减少的趋势，导致湖泊萎缩、湿地减少、草场退化、黄河断流等一系列环境问题[3]，严重影响流域经济发展和生态环境安全[4]。因此，研究水利工程的修建对黄河基流的变化量以及变化趋势的影响，有助于进一步研究黄河地下水循环规律以及保护流域生态环境。

基于不同的理论和方法，国内外学者提出了多种基流的分割方法，但由于缺乏试验资料以及涉及气候、水文、地质等多种科学，无法验证其结果的合理性，因此国内外学者尚未找出一种普遍认同的合理的基流分割方法[5-6]。图解法是依据流量过程线的几何特征，按照水文地质学者的经验以及流域的水文地质的特征来分割基流的传统基流分割方法。该方法操作简单，但计算结果粗略且只能对完整的波形进行分割。刘昌明等[7]进行基流划分时将最小月份的月流量作为基流。水量平衡法通过计算流域各水均衡项的水量得出基流量[8-10]，其优点是参数具有明确的物理意义，缺点是操作难度高且参数较多。熊立华等[11]通过对牧马河流域进行基流分割发现：非线性水库法与数字滤波法得到的结果相差不大，认为非线性水库法可以客观合理地反映实际退水过程。物理化学法主要应用于圈闭性和单一性的山区流域，具有可操作性和可信度高的优点，是未来基流分割方法研究的方向之一，但成本较高，目前只应用于小流域。陈宗宇等[12]利用稳定同位素及化学方法估算出张掖—正义峡河道段地下水向河道平均排泄量占该段河流径流量的69%。数字滤波法[13-14]通过计算机自动实现，具有操作性强、可重复、普遍适用，且人为因素干扰等优点，被广泛应用于黄河流域等广大地区。宋小圆等[5]运用4种数字滤波方法探究了不同滤波方法滤波次数、滤波参数变化对基流分割结果的影响，结果表明Lyne-Holick数字滤波法模拟效果较好，且省时高效。李瑞等[15]等通过研究黄河上游药水河流域发现Lyne-Holick数字滤波法具有较好的应用效果。党素珍等[16]通过数字滤波法研究发现：黑河上游基流量年内变化呈先增后减的趋势，年际变化呈增加的趋势。王雁林等[17]在探讨黄河河川基流量演化规律及驱动因子时指出，水库修建对河川基流量衰减影响不大。白乐等[18]在研究降水和人类活动对秃尾河流域基流的影响时发现，水利工程的建设在一定程度上削减了地下水对径流的补给。王蕊等[19]分析了近40 a黄河中游径流情势变化，指出1970年以来黄河中游岔巴沟流域淤地坝的大量修建造成河流年径流减少30%以上。林学钰等[20]认为在修建有大型水利工程的黄河干、支流的径流资料，需要经过人工还原后，方可进行基流分割。窟野河流域内现存坝库约850处，总容量近55万m3，这些大规模的水利工程改变了该流域原有的下垫面水循环规律，其水文效应在20世纪80年代初期开始凸显，对河川基流量的衰减有一定程度的影响[3]。前人对黄河基流量的研究中，研究水利工程的建设对河川基流量影响的较少。黄河是人类活动最为频繁的流域之一，而青铜峡水利枢纽位于黄河上游东部季风区与西部干旱区域的交汇地带，此地区气候干旱、降雨量少，该水利工程的修建改变了地表水与地下水的补排给关系[21]，因此，研究青铜峡水利枢纽建坝前后流域基流量的变化规律对维持该流域的生态平衡以及实现水资源优化配置具有重要作用。

图1 流域位置Fig.1 Location of catchment

本文采用Lyne-Holick单参数数字滤波法对黄河上游下河沿水文站、青铜峡水文站以及石嘴山水文站的流量过程进行基流分割，运用Mann-Kendall非参数检验法分析基流的变化趋势和突变点，分析青铜峡建坝前后基流的变化特征以及变化规律，揭示水利枢纽工程对河川基流量变化的影响程度，以期为该流域水资源优化配置和水文模拟提供理论支持。

1 研究区概况

1.1 流域概况

研究区域属于黄河上游宁夏段，从宁夏回族自治区中卫市下河沿水文站到石嘴山市石嘴山水文站(图1)。该流域位于36°0′N～39°23′N，104°17′E～107°39′E，南邻甘肃省，西部、北部与内蒙古自治区相连，东部结合陕西省。该流域地处中国地质、地貌“南北中轴”的北端，地形结构复杂，包含高原、山地、沙漠、平原等多种地形地貌。全流域从南向北表现出由流水地貌向风蚀地貌过渡的特征，黄河上游宁夏段属于大陆性半湿润半干旱气候，降水量少且年内分布不均匀，蒸发强烈[22-23]。

1.2 水利工程及水文站概况

青铜峡水电站位于黄河上游宁夏回族自治区青铜峡市青铜峡峡谷出口处，是一座以灌溉、发电为主，兼有防洪、防凌、城市供水等效益的综合性水利枢纽。坝址以上流域面积27.5万km2，多年平均年径流量324亿m3,多年平均流量1 028 m3/s，水库正常蓄水位1 156 m，相应设计库容为6.06亿m3，水库面积为113 km2，现存库容0.56亿m3。1958年8月青铜峡水电站正式开始动工，1967年12月投入使用，1968年2月第一台机组(2号机组)正式发电，1978年8台机组安装完毕。黄河上游石嘴山水文站位于宁夏回族自治区与内蒙古的交界处(106°47′E，39°15′N)，集水面积3 215.12 km2，是控制黄河水量从宁夏进入内蒙古的重要水文站。青铜峡水文站位于105°99′E，37°89′N，是青铜峡水电站的出水控制站。下河沿水文站是黄河宁夏段的入境站，位于105°05′E，37°15′N，集水面积为2 643.08 km2。

2 数据来源及研究方法

2.1 数据来源

本文选取下河沿、青铜峡和石嘴山水文站1951—2000年共50 a的日径流量时间序列作为基础数据，对实测径流数据进行水文频率分析等处理，采用Lyne-Holick单参数数字滤波法对日径流量进行基流分割得到日基流量，对基流量进行M-K检验以及趋势分析。

2.2 研究方法

2.2.1 数字滤波法

本文采用数字滤波法中的Lyne-Holick数字滤波法，该方法是由Lyne和Holick于1979年首次提出，1990年Nathan等[24]将其改进并首次应用于水文研究中[25]。

滤波方程为

(1)

基流方程为

bt=Qt-qt

(2)

式中：qt——t时刻的地面径流,m3/s；Qt——t时刻的总径流，m3/s；qt-1——t-1时刻的地面径流，m3/s;Qt-1——t-1时刻的总径流，m3/s；bt——t时刻的基流；α——滤波参数。根据Nathan和McMahon等的研究以及结合该流域气候、水文地质条件，通过实际切割效果对比比较，得到0.95作为滤波参数[16]，采用第3次滤波所得的基流量[19]更符合实际情况。

2.2.2Mann-Kendall非参数检验法

Mann-Kendall非参数检验法是一种不需要样本遵从一定的分布，也不受少数异常值干扰的非参数统计检验方法[26]，其能够很好地描述时间序列的趋势变化特征，具有适用范围广、定量化程度高、计算简便等优点。原假设H0为时间序列数据(x1，x2，…，xn)是n个独立的、随机变量同分布的样本，备择假设H1是双边检验[27]，对于所有的k，j≤n，且k≠j，xk和xj的分布是不相同的，检验的统计变量S为

(3)

其中

标准的正态统计变量计算公式为

(4)

其中

D(S)=n(n-1)(2n+5)/18

Z为正值表示增加趋势，Z为负值表示减少趋势，当Z的绝对值在大于等于1.28、1.64、2.32时分别表示通过90%、95%、99%置信度显著性检验[28]。

当Mann-Kendall检验进一步用于检验序列突变时，检验统计量与上述Z有所不同，通过构造一秩序列[29-30]。

(5)

(6)

式中:E(Sk)、D(Sk)——Sk的均值、方差。

按时间序列x的逆序xn,xn-1,…，x1重复上述过程，获得统计变量UBk(k=n,n-1,…，1)，同时使：UBk=-UFk。将UBk与UFk绘制成曲线，如果2条曲线出现交点,且该点处满足UBk<1.96、UFk<1.96，则认为该点为序列的突变点，检验置信水平为0.05[31]。

3 研 究 结 果

3.1 基流量和BFI的年际变化特征

在基流计算过程中，一般以年基流量以及基流量占总径流量的比重(即基流指数，baseflow index,简称BFI)来量化基流与径流之间的关系[26]。为了分析青铜峡水利枢纽工程对黄河基流的影响，绘制了下河沿、青铜峡以及石嘴山水文站的年际基流变化规律,如图2、图3所示。

图2 基流年际变化Fig.2 Annual variations of baseflow

图3 青铜峡径流基流年际变化Fig.3 Annual variations of runoff and baseflow in Qingtongxia

人为自然条件下为青铜峡大坝建设以前，人类活动较少，自然和人为条件干预下为大坝建设期，人为干预条件下为大坝建设完成投入使用以后，人类活动影响显著如图2所示(图中：Ⅰ为自然条件下，Ⅱ为自然和人为干预条件下，Ⅲ为人为干预条件下)。由图2、图3可知：3个站的河川基流量整体呈下降趋势，其中青铜峡站下降幅度最为明显。建坝前下河沿、青铜峡、石嘴山3个水文站多年平均基流量分别为214亿m3、184亿m3、192亿m3，建坝后3站多年平均基流量分别为187亿m3、111亿m3、160亿m3，相比于建坝前分别减少了12%、40%、17%，其中青铜峡站下降幅度最大。由以上数据可以看出，建坝前3个水文站的基流量相差不大，建坝施工期3个站的基流量开始出现差异，建坝后，青铜峡站基流量明显小于其他2个站，其中上游的下河沿站基流量最大。Lyne-Holick数字滤波法分割出的年平均径流与年平均基流的pearson相关系数均高达0.94，且都通过置信水平为0.01的显著性检验。

图4 BFI年际变化和建坝前后BFI对比Fig.4 Annual variations of BFI and comparison chart before and after dam construction

基流量从建坝施工期出现明显的减少，说明水电站的修建影响了地下水的时空分布的规律。青铜峡水库正常蓄水位1 156 m，水库蓄水提高了流域地下水排泄基准面，导致地下水排泄量减少，影响了地下水对河川径流的补给，从而使河川基流量减少。降水是河川基流量的主要补给来源，青铜峡水利枢纽的修建使流域气候发生了一定程度上的改变，从而使河川基流量呈现明显减少趋势。从图4中可以看出：下河沿、青铜峡、石嘴山站多年平均基流指数分别为0.65、0.54、0.60，BFI的年际变化过程基本保持一致，在20世纪60年代青铜峡水利枢纽建设时期均出现大幅度下降，其中青铜峡站BFI下降幅度最为明显，为13%，而石嘴山站总体BFI保持相对稳定。

3.2 不同水平年基流量和BFI变化特征

运用P-Ⅲ型频率曲线法对下河沿、青铜峡、石嘴山3个水文站1951—2000年的实测径流资料进行水文频率分析，划分出相应的丰水年、平水年、枯水年。根据丰水年、平水年、枯水年的划分标准(表1)，将3个站的年均基流和BFI按照建坝前和建坝后分别进行计算统计。

表1 年径流量丰水年、平水年、枯水年划分标准Table 1 Standard for division of wet, ordinary and dry annual runoff level

结合图5、图6可以看出，下河沿、青铜峡、石嘴山3个站年均基流均为丰水年>平水年>枯水年，相比于建坝前，建坝后3个水文站在丰水年、平水年、枯水年3个水平年中平均基流量均小幅度减小，其中只有下河沿水文站在枯水年基流量小幅度上升。下河沿站和石嘴山站BFI建坝前后变化不大，且在建坝前3个站BFI相差不大。青铜峡站基流指数下降幅度最明显，丰水年中下降13%、平水年下降20%，枯水年基流指数基本保持不变。丰水年降水量相对较大，径流主要由流域内的大气降雨进行补充，青铜峡水利枢纽的修建，改变了地表汇流的条件，减少了地下水对径流的补给。

图5 建坝前后基流量对比Fig.5 Baseflow comparison chart before and after dam construction

图6 建坝前后BIF对比Fig.6 BFI comparison chart before and after dam construction

3.3 趋势和突变分析

通过M-K非参数检验法进一步分析得出：下河沿、青铜峡、石嘴山3个水文站的Z值分别为-1.95、-4.23、-2.73，都小于0，表明3个站的基流量呈减少趋势，且青铜峡、石嘴山水文站Z值通过α=0.01的显著性检验，下河沿水文站通过α=0.05的显著性检验，由此可以看出青铜峡、石嘴山基流量减少趋势明显，且青铜峡减少趋势远远大于下河沿和石嘴山。由图2可以看出，1967—1968年间青铜峡径流量和基流量明显减少。

突变即序列收到某种条件的变化或受某种因素的影响，数据从一种稳定态(或稳定持续的变化趋势)跳跃式地转变到另一种稳定态(或稳定持续的变化趋势)[16]。本文采用M-K突变检验法确定青铜峡水文站基流数据的突变点。由图7可以看出，青铜峡基流量数据的突变点出现在1962—1964年，此时段处于青铜峡水利枢纽工程的建设时期，表明青铜峡水利枢纽的施工是造成该河段基流量突变的主要原因。

图7 青铜峡建坝前后基流突变分析Fig.7 Mutation analysis of baseflow of Qingtongxia dam before and after the construction

4 结 论

a. 下河沿、青铜峡、石嘴山水文站基流量总体上均呈减少趋势，建坝前3个站基流量和基流指数相差不大，基流量和基流指数分别为214亿m3、184亿m3、192亿m3和0.65、0.54、0.60，建坝后出现明显差异，下河沿站基流量最多，石嘴山站次之，青铜峡站最少，且青铜峡站减少比例最大，基流量减少40%、BFI减少13%，受青铜峡水电站影响最为显著。

b. 下河沿、青铜峡、石嘴山3个水文站的不同水平年的多年平均基流量由高到低依次为丰水年、平水年、枯水年，多年平均基流指数由高到低依次为平水年、丰水年、枯水年，但下河沿水文站多年平均基流指数为枯水年、平水年、丰水年。下河沿水文站位于青铜峡水电站的上游，受青铜峡水电站建设的影响较小。

c. 20世纪60年代以来，青铜峡水利枢纽的建设影响了该区域的水文地质条件，使原有的河流和地下水的补排给关系发生改变，导致基流量出现明显的减少趋势。通过M-K突变检验发现，突变点出现在青铜峡水电站建设施工期，表明青铜峡水电站的建设和水库蓄水，影响了当地地下水的时空分布规律，使河川基流量减少。