渭河上游水沙变化特征分析.pdf

DOI:1016616/jcnki10-1326/TV20230810渭河上游水沙变化特征分析唐骄宇1,2 卢 洋2(1 长安大学,陕西 西安 710000;2 黄河水利委员会三门峡库区水文水资源局,河南 三门峡 472000)【摘 要】本文以渭河流域天水测区多年水文资料为基础,运用线性趋势法、Mann-Kendall 秩次检验法、启发式分割法等多种数理统计方法,研究了 31 年来渭河上游降水、径流、泥沙的趋势性及突变性。EMD 分解利用非平稳信号分析降水、径流的趋势变化,通过双累积曲线分析径流-泥沙关系,累积输沙量的趋势线存在变缓趋势,说明在径流增加的情况下,输沙量逐渐减少,这与流域内的生态景观工程和近些年的水沙治理都有着密不可分的关系。【关键词】渭河上游;水文要素;水沙变化;趋势性;突变性;诊断分析中图分类号:TV14 文献标志码:B 文章编号:2096-0131(2023)08-051-09Analysis of the Characteristics of Water and Sediment Changes inthe Upstream of the Weihe RiverTANG Jiaoyu1,2,LU Yang2(1.Changan University,Xian 710000,China;2.Hydrology and Water Resources Bureau of Sanmenxia Reservoir of YRCC,Sanmenxia 472000,China)收稿日期:2022-10-18作者简介:唐骄宇(1996),女,助理工程师,硕士,主要从事水文要素规律分析等工作。Abstract:Based on long-term hydrological data in the Tianshui region of the Weihe River Basin,various mathematical andstatistical methods such as linear trend analysis,Mann-Kendall rank test,and heuristic segmentation method are employedto investigate the trend and abrupt changes in precipitation,runoff,and sediment in the upstream of the Weihe River over aperiod of 31 years.The EMD decomposition is used to analyze the non-stationary signals and identify the trend changes inprecipitation and runoff.The cumulative sediment discharge trend line exhibits a decreasing trend,indicating a gradualreduction in sediment transport with increasing runoff.This phenomenon is closely related to the ecological landscapeprojects and recent sediment control measures implemented in the watershed.Key words:upstream of the Weihe River;hydrologic features;water and sediment changes;trend;abrupt change;diagnostic analysis 渭河流域为干旱半干旱区域,在气候变化和人类活动的影响下,渭河上游降水、径流和泥沙均有不同程度的变化。渭河上游流域下垫面为黄土沟壑,水土流失严重,水资源供需矛盾突出,降水、径流和泥沙的年际和季节变化规律对于研究地区水资源开发利用有重要影响。因此,研究区域内降水、径流、泥沙变化规律,为渭河上游防汛抗旱、生态环境保护、水资源规划和开发利用提供技术支撑,对于分析水循环规律、预测极端15水文事件等具有重要意义。本文以渭河上游天水水文勘测区(以下简称“天水测区”)为研究区域。天水测区上起渭河武山,下至天水市麦积区东岔镇,包括区间支流散渡河、葫芦河、藉河、牛头河流域,测区内有黄土高原、秦岭、六盘山,年降水量地域分布极不均匀。采用的水文资料系列为19902020 年,年降水量采用测区各雨量站实测数据,据此计算出流域平均面雨量;年径流量、年输沙量采用北道水文站实测值,利用线性趋势法、肯德尔秩次检验法、启发式分割、有序聚类法、季-海哈林法、双累积曲线法等方法研究了降水、径流和泥沙的变化特点、趋势性、突变性,进而分析变化的影响原因,分析成果可为流域水资源合理开发利用与配置提供参考。1 分析计算方法1 1 趋势性分析分别采用线性回归法、非参数检验法、肯德尔秩次检验法等进行趋势性分析。1 1 1 线性回归法线性回归法构建水文变量 x 与其时间 t 的一元线性回归方程:x(t)=a+bt(1)回归系数 b 的显著性可通过 t 检验进行判断,构建的 T统计量:T=rn-21-r2(2)服从自由度为(n-2)的t 分布,给定显著性水平,可查得 t/2,如果 T T/2,则认为回归效果显著,存在线性趋势;否则,线性趋势不显著。1 1 2 非参数检验法非参数检验法是对原始数据的秩进行分析,不需要样本服从一定的分布,根据样本信息对总体分布情况进行推断和处理,并针对其不同条件下的参数进行优化和整理。1 1 3 肯德尔秩次检验法利用肯德尔秩次检验法(Mann-Kendall 法,简称M-K)1对水文序列的趋势性采取显著性检验方式,是一种非参数方法(亦称无分布检验),其优点是样本不需要遵从一定的分布,也不受少数特异值的干扰,计算过程相对简单,被广泛应用于水文变量、要素等非正态分布的趋势分析中,计算出相应的检验值,对水文序列的趋势性和突变性进行有效的检验,其原理见表 1。表 1 M-K 趋势判断表(显著性水平=0 05)统计量z z1-/2z z1-/2Z0显著上升趋势不显著的变化趋势Z 00,s=0s+1var(s),s 00(xj-xi)=0-1(xj-xi)10 时,s 为近似服从正态分布;Var(s)为方差;m 为数据相同的组数;tk为与第 k 组的数据相同的个数。给定显著性水平,依据表 1 可以判断序列的上升或下降趋势显著情况。以 Kendall 倾斜度 量化序列,可表示为=Medianxi-xj()i-j j i(7)式中:l j i 0时,呈现上升趋势,反之呈下降趋势。251 2 跳跃成分分析跳跃成分分析分别采用有序聚类法、季-海哈林法、启发式分割法对降水、径流进行突变点诊断,年输沙量除采用前述三种方法外,还采用了双累积曲线法进行突变点诊断分析。1 2 1 有序聚类法2有序聚类法利用最优分割点来推求突变点,基本原则是使同要素之间的离差平方和较小。设有水文序列x1,x2,xn,假设可能的突变点为(2 n-1),则突变点前后的离差平方和分别为V=i=1(xi-x)2(8)Vn-=ni=+1(xi-xn-)2(9)式中:x和xn-为 前后两个序列的均值,即前后两个序列的聚类中心。有序聚类法常用的目标函数为S=min2n-1S0()=min2n-1(V+Vn-)(10)当式(10)中 S 取极小值时,对应的 为最优分割点,可推断为突变点。1 2 2 季-海哈林检验法2对于水文时间序列 x(t),假定总体为正态分布,且可能变异点 的先验分布为均匀分布的情况下,推得 的后验分布为f(x1,x2,xn)=k n/(n-R()-n-22(1 n-1)(11)R=t=1(xt-x)2+nt=+1(xt-x)2/nt=1(xt-xn)2(12)式中:k 为比例常数;n 为样本容量。当统计值 f()取最大值时,满足后验分布条件,此时对应的 值为最可能分割点,由此确定最可能变异年份。1 2 3 启发式分割法相关研究3-4表明启发式分割能在不同时间尺度高效地把非平稳时间序列划分成多个不同均值的子序列,并且相比多种常用突变检测方法,能有效地排除虚假的变异点。对于时间序列 XX1,X2,Xn,其各点的合并偏差可采用下式计算:SD=(N1-1)S21+(N2-1)S22N1+N2-()21/21N1+1N()21/2(13)式中:SD为合并偏差;N1、N2分别为该点左边部分和右边部分的点数;S1、S2分别为该点左边部分和右边部分的标准差。各点的统计量 T 计算公式为T=1-2SD(14)式中:1、2分别为该点左边部分和右边部分的均值。统计量 T 最大值 Tmax的置信度 P(Tmax)的计算式为P(Tmax)1-Iv/(v+Tmax2(v,)(15)当 P(Tmax)达到指定的置信度时,则该点为显著的均值突变点,对原序列进行分割,之后对分割后的子序列继续检测。若 P(Tmax)未能达到指定的置信度区间,则不再分割。2 计算结果分析2 1 降水2 1 1 降水趋势变化天水测区平均年降水量年际变化见图 1,其多年平均年降水量为 534 9mm,最大为 755 4mm(2003年),最小为 361 2mm(1997 年)。年降水量趋势线表明 31 年来降水总体呈平稳微弱上升趋势。线性回归法、非参数检验法、肯德尔秩次检验法均一致诊断为趋势不显著,但计算检验值与标准值较接近(见表 2),说明上升与非上升趋势处于临界状态。图 1 天水测区年降水量年际变化35唐骄宇等/渭河上游水沙变化特征分析 表 2 天水测区年降水量序列趋势诊断结果诊 断 方 法检 验 值标 准 值检 验 结 果线性回归法|T|=1 62T(0 05/2)=1 64趋势不显著,微弱上升非参数检验法|Z|=0 239Z(0 05)=0 248趋势不显著,微弱上升肯德尔秩次检验法|U|=1 89U(0 05/2)=1 96趋势不显著,微弱上升2 1 2 降水突变点诊断对测区平均年降水量进行 M-K 突变检验,并设置0 05 显著性水平,即得到置信区间,临界值为-1 96 和1 96,分析结果见图 2。根据分析,年降水量 UF、UB 曲线在 2010 年后出现多个交点,相交于置信期间内,且两曲线趋势一致,通过 UB、UF 交叉点可知,突变点在2017 年前后。采用有序聚类法和季-海哈林法进行同步验证,得出 相 同 检 验 值,T=2 42,T(0 05/2)=1 96,T T(0 05/2),根据检验值趋势图,有序聚类法分析出其 2017 年资料系列达到最小值,季-海哈林法分析出其 2017 年资料系列达到最大值,根据突变点判别,在 2017 年前后均值发生显著跳跃,见图 3、图 4。图 2 天水测区年降水量 M-K 秩次检验图 3 天水测区年降水量有序聚类法检验2 2 径流2 2 1 径流趋势变化对北道水文站 19902020 年年径流数据进行统计分析,得到径流变化序列,见图 5。可知渭河上游多图 4 天水测区年降水量季-海哈林法检验年平均年径流量为 7 32 亿 m3,最大年径流量为 18 97亿 m3(2020 年),最小年径流量为 1 29 亿 m3(1997年)。由图 5 可知,北道水文站年径流量总体呈现升高趋势,线性变化率为0 0952 亿 m3/a,5 年滑动平均显示45图 5 天水测区年径流量年际变化近 31 年来年径流量呈现缓慢上升趋势。线性回归法、非参数检验法、肯德尔秩次检验法均一致诊断为“趋势不显著”,但计算检验值与标准值较接近(见表 3),说明上升与非上升趋势处于临界状态。2 2 2 径流突变点诊断北道水文站年径流量 M-K 突变检验结果见图 6,设置 0 05 显著性水平,即得到置信区间临界值为-1 96和1 96,根据分析,得出2013 年年径流量 UF U(0 05/2)=1 96,根据检验值趋势图可分析出,在 2017 年资料系列达到最值,根据方法的突变点判别,均在 2017 年前后均值发生显著跳跃,见图7、图 8,与 M-K 检验结果一致。根据以上方法,设定置信度为 0 95,最小分割尺度为 25,根据年径流变异分析结果,天水测区年径流量启发式分割检验见图 9。由图 9 可知,2018 年对应的 T 值最大,并且相对应的 P(Tmax)为 0 61,小于临界值 0 95,说明 2018 年为第一突变点5-6。2 2 3 降水径流 EMD 结果分析运用 EMD 方法对天水测区多年来降水和径流数据进行多尺度分解,结果见图 10、图 11,其中包含有多个具有物理意义的平稳固有模态函数(Intrinsic ModeFunctions,IMF)和具有单一性的趋势项(Re-sidual,Res),从中可得出如下重要结论:a 根据 EMD 原理,渭河上游天水测区降水量,原始分量除外,可分解成 3 个不同波动周期的振荡分量,55唐骄宇等/渭河上游水沙变化特征分析 图 7 天水测区年径流量有序聚类法检验图 8 天水测区年径流量季-海哈林法检验图 9 天水测区年径流量启发式分割检验径流量可分解成 4 个不同波动周期的振荡分量,反映了该测区降水和径流在不同变化程度上的多时间尺度性。b 20 年代初及 2015 年前后降水和径流 IMF1 分量变动幅度大,降水是径流变化的主要影响因素,周期较为一致。IMF1 是第一个本征模态函数,波动振幅最大,频率最高,波长最短,周期 2 5 年,随后的分解过程中,本征模态函数振幅逐渐减小,频率逐渐降低,波长逐渐变大,其对原始时间序列的影响程度逐渐降低。c Res 分量显示的是天水测区年降水量和年径流图 10 年降水量序列 EMD 分解65图 11 年径流量序列 EMD 分解量的整体变化趋势,近 31 年来整体呈急剧上升趋势。EMD 通过信号逐步去噪和趋势剔除7,实现了序列的平稳化处理。2 3 泥沙变化2 3 1 输沙量趋势变化对北道水文站 19902020 年年输沙量数据进行统计分析,得到变化序列,见图 12。渭河上游多年平均年输沙量为 3490 万 t,最大年输沙量为 15600 万 t(1992 年),最小年输沙量为 298 万 t(2014 年)。北道水文站年输沙量总体呈现下降趋势,线性变化率为198 万 t/a,5 年滑动平均显示近 31 年来输沙量呈现持续下降趋势。采用线性回归法、非参数检验、肯德尔秩次检验方法分析结果一致,均为下降趋势显著,见表 4。75唐骄宇等/渭河上游水沙变化特征分析 图 12 天水测区年输沙量年际变化表 4 天水测区年输沙量序列趋势诊断结果诊断方法检 验 值标 准 值检验结果线性回归法|T|=3 93T(0 05/2)=1 64趋势显著非参数检验法|Z|=0 462Z(0 05)=0 248趋势显著肯德尔秩次检验法|U|=3 65U(0 05/2)=1 96趋势显著2 3 2 泥沙突变点诊断对北道水文站年输沙量进行 M-K 突变检验,见图13。设置 0 05 显著性水平,即得到置信区间临界值为-1 96 和 1 96。根据分析得出,在 19902020 年,UF0,说明年输沙量呈显著下降趋势。20 世纪 90 年代以来,UF 和 UB 曲线存在较多交点,整体来看,相关检验值 U=3 651 U(0 05/2)=1 96,根据检验值趋势图分析出,在 1992 年资料系列达到最值,根据方法的突变点判别,资料系列在 1992 年前后均值发生显著跳跃,见图 14 和图 15,与 M-K 检验结果一致。图 13 天水测区年输沙量 M-K 秩次检验图 14 天水测区年输沙量有序聚类法检验2 3 3 径流-泥沙关系分析通过绘制径流-泥沙双累积曲线(见图 16),可以图 15 天水测区年输沙量季-海哈林法检验看出累积曲线显著分为 3 个阶段,19901992 年为一个变化时段,该时段相关系数为 0 9994;1993200285年为一个变化时段,该时段相关系数为 0 9935;20032020 年为一个变化时段,该时段相关系数为 0 9563;每个阶段相关性均较好。由图 16 可知,1992 年、2002年前后均为输沙量的突变点,累积输沙量的趋势线存在变缓趋势,说明在径流增加的情况下,输沙量逐渐减少,这与流域内的生态景观工程和近些年的水沙治理成效都有着密不可分的关系。图 16 径流-泥沙双累积曲线3 结 语经过分析,天水测区属干旱半干旱区,径流主要由降水补给,南北两岸支流有不同水文特性,南岸面积小,地形陡峻,地表多为森林覆盖或为剥蚀山地,降水较多,植被好,支流距渭河干流距离短,水多沙少;北岸支流水少沙多,水量较贫,主要为暴雨洪水,陡涨陡落。流域降水、径流分别采用线性回归、非参数检验、肯德尔秩次检验三种方法分析,结果一致,均为趋势性不明显,但统计值与判断值较接近,说明处在临界状态,而采用滑动平均曲线均有微弱上升趋势。分别采用 M-K秩次检验、启发式分割、有序聚类法、季-海哈林突变诊断方法进行同步验证,得出降水和径流均在 2017 年前后产生突变一致性结论。对于年输沙量分别采用线性回归、非参数检验、肯德尔秩次检验三种方法分析,结果一致,均为下降趋势明显。采用启发式分割、有序聚类、季-海哈林突变诊断方法得出,在 1992 年前后产生突变的一致结论,采用径流-泥沙双累积曲线法分析,明显发现年输沙量存在 1992 年和 2002 年两个突变点。天水测区在降水、径流无明显减少的情况下,年输沙量减少趋势明显,说明流域内土地利用、生态保护、植被恢复等人类活动剧烈,特别是水土保持工作成效显著。采用径流-泥沙双累积曲线法发现年输沙量存在两个突变点,但基于数理统计分析判断的启发式分割、有序聚类、季-海哈林突变诊断方法年输沙量诊断为一个,可见数理统计分析判断方法对于多个突变点的诊断有待进一步完善。参考文献1 姜瑶,徐宗学,王静.基于年径流序列的五种趋势检测方法性能对比J.水利学报,2020,51(7):845-857.2 郭明,杨志勇,周玖辉,等.非一致性条件下嘉陵江上游水文频率分析J.华北水利水电大学学报(自然科学版),2018,39(1):42-47.3 陈广才,谢平.基于启发式分割算法的水文变异分析研究J.中山大学学报(自然科学版),2008(5):122-125.4 龚志强,封国林,万仕全,等.基于启发式分割算法检测华北和全球气候变化的特征J.物理学报,2006,55(1):477-484.5 陆鹏翔.基于启发式分割和 TFPW-MK 的绥江流域极端降雨变异诊断J.陕西水利,2021(2):44-46.6 黄生志,黄强,王义昆,等.基于启发式分割和近似熵法的径流序列变异诊断J.中山大学学报(自然科学版),2014,53(4):154-160.7 冯平,丁志宏,韩瑞光.基于 EMD 的洮河年径流量变化多时间尺度分析J.干旱区资源与环境,2008,22(12):73-76.95唐骄宇等/渭河上游水沙变化特征分析