中部引黄工程水源地年径流长期演变规律分析.docxVIP

中部引黄工程水源地年径流长期演变规律分析 黄河是中国的第二条河流，全长564公里，流域面积79.5万公里。随着经济发展和人们生活水平的提高, 黄河水资源供求形势日趋紧张, 大型引黄调水工程也越来越多。为此, 近年来许多专家学者作了大量的工作, 从不同的角度来分析其长期演变规律, 并且取得了大量有意义的成果。如Shi Changxing等基于长期水文气象系列、数字高程模型、土壤数据库和黄河中游水土保持措施的文件等多个方面分析了黄河中游年径流量的变化, 并研究了其与水土保持和气候变化的关系。Wang Suiji等利用累积异常识别法、线性回归法和改进的累积量斜率变化率分析方法 (SCRCQ) 分析了气候和人类活动对黄河径流变化的贡献。还有学者分别利用希尔伯特-黄变换 (HHT) 法、重标度极差 (R/S) 法、Mann-Kendall秩次相关检验、流量历时曲线法、双累积曲线法、小波分析方法的一种或几种对黄河径流长期变化趋势进行了分析研究。但是, 这些研究所采用的研究方法相对较为单一, 其研究结果具有一定的局限性, 不能全面的揭示黄河年径流量长期演变的时序变化特征及丰枯演化规律。 为了能够全面深入分析中部引黄工程水源地年径流量的长期演化规律, 本研究结合Kendall秩次相关检验法、R/S分析法、游程分析法和马尔可夫过程分析法这4种方法的优点, 并将这些方法进行了有机的结合。首先应用Kendall秩次相关检验法对径流量的时序变化趋势进行分析, 并在其基础上, 应用R/S分析法分析径流序列的持续性;然后, 应用游程分析法和马尔可夫过程分析法对径流的丰枯演化规律进行分析。目的在于为黄河水资源的可持续开发利用和中部引黄工程供需调节提供科学依据。 1 基本理论和方法 1.1 u3000da-qpb的显著性检验 对于序列x1, x2, x3, ...xn, 先确定所有对偶值 (xi, xj, ji) 中的 (xixj) 出现的次数, 记为P, 根据计算的P, 代入下列公式: 式中:τ为Kendall秩次统计量;U为Kendall秩次相关检验系数, 随n的增加很快收敛于标准状态分布;n为系列长度。 若给定显著性水平为α, 则其临界值为Uα/2。若|U|Uα/2, 则说明序列具有明显趋势性;当U0时, 则具有上升趋势;当U0时, 则具有下降趋势。若|U|Uα/2, 则序列趋势性不明显。 1.2 时间序列趋势结构的正持续性分析 将R/S分析法应用于径流时序特性分析时, 判断径流时序状态的持续性, 可通过计算Hurst指数。R/S分析法基本原理如下: 对于时间序列{x (t) } (t=1, 2…, n) , 下列量可定义为: 均值序列: 累计离差: 极差序列: 标准差序列: 如果比值R (τ) /S (τ) ≡R/S存在关系:R/S∝τH, 则说明时间序列{x (t) } (t=1, 2…, n) 具有Hurst现象, H称为Hurst指数。在双对数坐标系 (lg (τ) , lg (R/S) ) 中用最小二乘法拟合, 拟合直线的斜率即为H。 根据H的大小, 可以对时间序列趋势成分所表现的持续性进行判断。Hurst指数取值范围为[0, 1], 其存在以下3种情况:H=0.5时, 表明该时间序列变化是随机的。0.5H1时, 表明该时间序列具有长程相关性, 即过程具有正持续性。0H0.5时, 表明该时间序列具有反持续性。 1.3 连枯连丰年统计长度s1 游程分析法的计算可用以下公式: 式中:P为连续k年枯水 (丰水) 的概率;k为连枯 (连丰) 水年的年数;ρ为分布参数, 其值在0~1之间;S1为时间序列中连枯年 (连丰年) 总年数;S2为时间序列中各种统计长度的连枯 (连丰) 年发生的累积频次;S为时间序列枯水年 (丰水年) 的累积频次。 1.4 马尔可夫过程与多步转移概率矩阵 马尔可夫分析法是指在马尔可夫过程的假设前提下, 通过分析随机变量的现时变化情况来预测这些变量未来变化情况的一种预测方法, 又称为马尔可夫转移矩阵法。马尔可夫一步转移概率矩阵, 用P (1) 表示;多步转移概率矩阵, 用P (m) 表示, 它是由初始状态经多步 (m步) 转移从而达到相应状态的概率矩阵, 则P (1) 、P (m) 可分别表示为: 一步转移概率矩阵与m步转移概率矩阵的关系为: 马尔可夫过程 (有限状态的) 在经过长时间的转移后, 会达到平稳状态, 此后过程的状态将不再随时间的变化而变化。此时的概率被称为稳定概率, 又称为极限概率, 可定义为: 系统处于任意特定状态的概率可以由极限概率来代表, 该特定状态重复再现的平均时间可用1/P表示。 2 入库站和工程工况 由于黄河河曲水文站是国家基本水文站、中部引黄工程水源地天桥水库入库站, 本研究采用黄河河曲站1978-2