沙坪水闸年最高水位的频率线型分析.docx

沙坪水闸年最高水位的频率线型分析 论文 关键词：年最高水位 统计分析 频率线型 沙坪水闸  论文摘要：以沙坪水闸年最高水位的统计分析为例，本文对P-Ⅲ、指数Γ分布线型和对数Γ分布线型与经验频率点据拟合结果进行了比较。作为一种尝试，采用指数Γ分布线型和对数Γ分布线型拟合沙坪水闸年最高水位资料系列进行频率特性分析，取得了较为合理的分析研究成果。研究结果表明对数Γ分布线型优于P-Ⅲ和指数Γ分布线型，验证了沙坪水闸年最高水位频率分布线型服从对数Γ分布线型。中国编辑。  洪水频率统计分析是防洪排涝 科学 决策和规划设计的重要依据。在水文频率 计算 中，规范[1]推荐采用P-Ⅲ线型，但同时规定，经分析论证,也可采用其他线型。本文在 文献 《广义Γ分布的特性和应用》(金光炎、董秀颖,XX) [2]研究的基础上，尝试应用指数Γ分布线型和对数Γ分布线型拟合沙坪水闸年最高水位资料系列进行频率特性的统计分析。  1.沙坪水闸年最高水位的描述性统计分析  沙坪水闸位于西江下游右岸沙坪河入江河口处，见图1。该水闸是广东省鹤山市的一座以防洪、排涝为主，结合蓄水灌溉和改善航运的综合运用的中型水闸。该闸有从1964年自建闸以来至XX年43年的闸外年最高水位观测资料系列，见表1。资料系列的最小值、最大值分别为和,为此将年最高水位系列分为11个组限,其中组距宽度hi为的9组,其他组距2组。计算在每个组限内的年最高水位出现的频数fi，相应的年最高水位各组距区间上的经验频率为fi/n(n=43,i=1,2,…,11)。沙坪水闸年最高水位经验频率分组计算详见表3。在各组距区间上作以平均频率密度fi/(nhi)为高的直方图,见图4。从图4可看出直方图呈近似的正态分布形态,它有一个中间主峰，两头低，主峰在区间＜x≤上；年最高水位平均值、中位数都落在区间＜x≤上,经计算平均值为,中位数为。用矩法计算年最高水位经验系列统计量分别为:标准 差σ   12m，变差系数Cv=，偏态系数Cs=-，峰态系数Ce=-。  鹤山市   鹤山市   图1 沙坪水闸地理位置示意图 . P-Ⅲ、指数Γ分布和对数Γ分布线型比较分析  《水文分析与计算》(刘光文,1963)中提出了径流频率分析适线的线型选择原则[3]:⑴在计算简便的同时,具有尽量高的精度和弹性;⑵曲线与经验频率点据得到最好的拟合;⑶曲线的形状大致符合水文现象的一定物理性质,如曲线应该有一定的极限,不该出现负特征值。这一原则同样适用在年最高水位统计分析中,其中第⑶点的实质是合理性要求,《水文频率分析述评》(金光炎,1999)也认为按照水文物理概念,曲线应该有上限,并对Slade(1936)、谢家泽(1958) 等人的研究观点作了介绍[4]。因此本文在上述原则的基础上对备选的3种频率线型进行比较。  .1 P-Ⅲ、指数Γ分布和对数Γ分布线型的密度函数表达式及计算方法  P-Ⅲ密度函数[2]为式⑴：   式⑴中有3个参数。密度函数的定义域为当β0时a0≤x 四参数指数Γ分布线型[2]的密度函数如式⑵：   四参数指数Γ分布定义域为λ≤y 三参数对数Γ分布线型[2]的密度函数如式⑶：    三参数对数Γ分布定义域为,当β0时,A0≤y0时，-∞    在采用P-Ⅲ、四参数指数Γ分布和对数Γ分布线型频率累积曲线函数拟合适线过程时，由于参数都为非线性关系形式，不能通过某种转换变为线性形式，因此只能采用非线性迭代回归的办法求解。进行非线性迭代回归时，首先确定频率累积曲线函数的表达式，确定参数的初始值，然后根据某种方法进行搜索迭代，反复调整初始值，按规范应用最小二乘法原理使得观测值与拟合值的离差平方和最小时(或者结合其他一些条件)结束迭代过程，得到各参数的最后计算结果。另外，P-Ⅲ线型可利用已有的Φ值表、四参数指数Γ分布线型可利用三参数指数Γ分布线型已有的模比系数Kp值表通过变量坐标的平移转换得到适线结果。四参数指数Γ分布和对数Γ分布线型都比式⑴的P-Ⅲ线型多了1个参数，因此适线弹性要比P-Ⅲ线型大。在当今计算机技术被广泛应用于各领域的条件下,就计算方法难易程度与精度比较而言与P-Ⅲ分布线型相当。本例适线结果如图2、3所示，点绘在概率格纸上的经验频率点据上部呈现出向下凹的分布