Page 22 - 2024年第55卷第6期
P. 22

图 1 长江上游干支流 29座重点大型梯级水库概化图


















                       图 2 长江上游大型水库的水库系数                           图 3 三峡坝址年最大洪峰和 15d洪量变化趋势

              3.2  非一致性洪水频率分析  三峡水库考虑历史洪水的非一致性洪水频率分析资料系列从 1153—
              2022年,其中历史洪水调查期内(1153—1876年)确认的历史洪水、水文站实测大洪水见表 2,最终纳
              入历史洪水、作为特大值处理的洪水峰量均在括号中标明相应排位。需要强调的是在时变矩模型中,
                                                          、    的取值会影响时变 P - Ⅲ型分布的函数值,因此
              虽然洪水峰量和 RI系列已经给定,但是参数 ω 11 ω 21
                                                 的取值而变化。为了合理限制协变量 RI对洪水频率的影响,本
              同一个洪水对应的排位会随着 ω 11             、ω 21
              研究对时变 P - Ⅲ型分布进行参数估计时,维持特大洪水数量始终为 a且对应的排位不变(即式(1)所
              示历史洪水与其经验频率不变)。
                  采用年最大值取样法对 1877—2022年三峡水库坝址(宜昌站)日流量资料系列进行取样,并与表 2
              所示历史洪水共同组成一个具有 870a调查考证期的不连序年最大系列,基于 GAMLSS模型构建时变
              P - Ⅲ型分布参数与水库系数 RI之间的广义可加模型,由于三峡水库设计洪水初设成果的 C值普遍小
                                                                                                  s
              于 1.0,可以采用极大似然法估计时变 P - Ⅲ型分布参数,除此之外不对参数优化边界设限,使最终优
              化参数不在边界处取得。表 3对比分析了极大似然法和时变 P - Ⅲ型适线法的参数估计值与拟合优度评
              价指标,可以看出:①K - S检验的 p值均大于 0.05,接受原假设,即认为实测洪量数据的经验频率分
              布与时变 P - Ⅲ型分布均服从(0,1)区间上的均匀分布,因此可以采用时变 P - Ⅲ型适线法进行参数估
                                                                  的估计值均为负,而根据位置、尺度参数与
              计。②极大似然法和时变 P - Ⅲ型适线法对参数 ω 11                  和 ω 21
                                        - 1
              样本总体均值的关系 EY = αβ + a 可知,水库系数 RI对洪水峰量的均值具有削减作用。③Q 、W 、
                                                                                                    m
                                            0
                                                                                                         3d
                     8
                —  6 4  —
   17   18   19   20   21   22   23   24   25   26   27