山西省浊漳河北源流域洪水峰量演变规律及联合分布研究.docxVIP

?

?

山西省浊漳河北源流域洪水峰量演变规律及联合分布研究

?

?

宋佳佳张芳齐

摘要：为研究山西省浊漳河北源流域洪水特征，采用浊漳河北源石栈道水文站1957～2016年共60a逐日流量资料，对洪峰流量、1日洪量进行研究。通过Mann-Kendall非参数检验法，揭示了该水文站洪峰流量和1日洪量演变特征及规律，并通过Copula函数研究了浊漳河北源流域洪水峰量联合分布情况。结果表明：洪峰流量和1日洪量均呈现明显的下降趋势;时间序列上洪峰流量在1971年后有显著的下降趋势，1日洪量在1974年后呈现显著的下降趋势;Copula函数可较好拟合浊漳河北源洪水峰量的相关关系。研究成果可为该流域水利工程规划设计和风险评估提供科学依据。

关键词：洪峰流量;1日洪量;Mann-Kendall;Copula函数;浊漳河北源流域

中图法分类号：TV122.5文献标志码：ADOI：10.15974/j.cnki.slsdkb.2021.10.001

：1006-0081（2021）10-0006-06

0引言

在全球气候变化和人类活动的影响下，流域内洪水情势发生变化，需对其变化和发展趋势进行分析。洪水峰量是研究洪水的重要特征，它不仅反映洪水的强度，也预示防洪的级别，因此对洪水峰量的演变规律研究尤为重要[1]。洪水需要多个特征变量才能完整描述，如洪峰、峰现时间、时段洪量和洪水历时等存在相关性的特征变量。多个特征变量的联合分析相比于单变量分析能更好地描述洪水特性并将有利于建设安全的防洪工程系统。Copula函数是多变量联合分布构建理论与方法的重大突破，DeMichele[2]采用Copula函数建立了降雨强度和降雨历时的联合概率分布。熊立华等[3]采用Copula函数构建了长江流域某站点洪峰和洪量的联合分布。近十几年来，利用Copula函数研究洪水的成果众多，主要集中在洪水峰量、洪水历时等变量的联合分布和重现期研究上[4]。多变量重现期以联合重现期、同现重现期和Kendall重现期应用最为广泛[5]，结构荷载重现期[6]也引起了关注。

为较全面研究浊漳河北源流域洪水极端事件，需研究多变量的联合分布，Copula函数是实现这种相关性分析的有效方法。本文以浊漳河北源流域为例，采用石栈道水文站实测水文资料，利用Mann-Kendall非参数检验法，揭示了浊漳河北源洪水峰量的演变特征及规律。洪水峰量均单独服从P-Ⅲ型分布，通过Copula函数将浊漳北源流域洪水峰量结合起来建立联合分布模型，为水利工程规划设计和风险评估提供科学依据。

1研究区概况

浊漳河属漳卫南运河水系，上游有北、南、西三大支流，称浊漳北源、浊漳南源、浊漳西源。浊漳河北源是浊漳河三源之一，发源于山西省晋中市榆社县社城镇焦红寺大牛村，在山西省长治市襄垣县小峧村南汇入浊漳河干流（图1）。浊漳北源流域总面积3797km2，主河道长135km。

2数据与方法

2.1数据来源

本文数据来源为浊漳河北源河道榆社县石栈道水文站，利用石栈道水文站1957～2016年共60a逐日流量资料，采取年内日最大选样法选取洪峰流量、1日洪量为研究对象，运用Mann-Kendall检验法对洪水峰量进行演变规律分析，并采用Copula函数对石栈道水文站站址处洪水峰量的联合分布进行研究。

2.2数据资料“三性分析”

2.2.1可靠性

石栈道水文站位于晋中市榆社县石栈道村浊漳河北源河道上，地理坐标为东经112°58′、北纬37°04′，于1957年6月设立。该站属于山西省水文水资源勘测局，为基本水文站。该站观测编制水文资料均符合有关规程、规范要求，资料质量可靠，精度满足水文分析要求。

2.2.2一致性

石栈道水文站上游约15km处建有双峰水库，控制流域面积190.2km2。双峰水库自1975年兴建以来，共经过两次续建，于2008年完成全部工程，达到现状规模。由于库区移民问题，2017年才开始蓄水。本文选取石栈道水文站1957～2016年实测流量系列资料具有一致性。

2.2.3代表性

石栈道水文站实测流量系列资料为1957～2016年，流量系列相对较长，包括了丰、平、枯年。绘制石栈道水文站模比差积曲线及洪峰流量和1日洪量变化过程线（图2～3）。石栈道水文站年最大洪峰流量和1日洪量变化规律基本一致，即大洪峰大洪量，小洪峰小洪量，但年最大洪峰流量和最大1日洪量不一定出现在同一日。石栈道水文站实测流量系列资料具有代表性。

2.3研究方法

2.3.1Mann-Kendall非參数检验法

本文采用Mann-Kendall非参数检验法研究水资源演变特征[7]，选择95%的置信水平检验，取值为1.96。同时，利用Mann-Kendall非参数检验法进行突变识别。

2.3.2