[发明专利]基于分位数回归的非一致性设计洪水计算方法

摘要

本发明公开了一种基于分位数回归的非一致性设计洪水计算方法，通过确定与流域站点年最大洪峰流量系列存在相关关系的物理因子作为自变量，收集年最大洪峰流量和物理因子的同步资料系列，构建年最大洪峰流量与物理因子之间的分位数回归模型，并估计分位数回归模型系数，在进行拟合优度检验和回归系数的显著性检验的基础上，推求非一致性条件下指定重现期的设计洪水。本发明不需要假设特定的洪水概率分布类型，且以与因变量最大洪峰流量相关的物理因子作为自变量，可以解释导致洪水序列非一致性的原因。另外，可以同时描述年最大洪峰流量序列的趋势性和跳跃式变化，避免了计算的设计洪水无限增大或者无限减小的不合理现象。

主权项

1.一种基于分位数回归的非一致性设计洪水计算方法，其特征在于包括以下步骤：步骤1，确定与流域站点年最大洪峰流量系列存在相关关系的物理因子作为分位数回归模型的自变量，所述物理因子包括气候变量和人类活动变量，其中，气候变量包括降水和气温，人类活动变量包括水土保持面积和水库因子；水库因子是一种用来反映大坝和水库对洪水变化影响的水库指标，计算公式如下：其中，RI为水库因子；i为水文站点上游水库编号，i＝1,2,…,N；N为水文站点上游水库数量；Ai为水库集水面积，km2；AT为水文站点集水面积，km2；Ci为水库调洪库容，亿m3；CT为水文站点汛期多年平均径流量，亿m3；步骤2，收集因变量年最大洪峰流量和步骤1中确定的自变量物理因子的同步资料系列；步骤3，根据步骤2中的资料系列，构建年最大洪峰流量与物理因子之间的分位数回归模型，并估计分位数回归模型的系数；其中，采用多元线性分位数回归模型作为分位数回归模型，其中因变量个数为1，自变量个数为4；步骤4，对步骤3中建立的分位数回归模型进行拟合优度检验和回归系数的显著性检验；步骤5，基于步骤4中统计检验通过的分位数回归模型，推求非一致性条件下指定重现期的设计洪水。