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

沙坪水闸年最高水位的频率线型分 析 论文关键词：年最高水位统计分析频 率线型沙坪水闸 论文摘要：以沙坪水闸年最高水位的统 计分析为例，本文对 P-m、指数r分布线 型和对数r分布线型与经验频率点据拟合 结果进行了比较。作为一种尝试，采用指数 r分布线型和对数 r分布线型拟合沙坪水 闸年最高水位资料系列进行频率特性分析， 取得了较为合理的分析研究成果。研究结果 表明对数「分布线型优于P-m和指数「分 布线型，验证了沙坪水闸年最局水位频率分 布线型服从对数r分布线型。 洪水频率统计分析是防洪排涝科学决 策和规划设计的重要依据。在水文频率计算 中，规范［1］推荐采用P-m线型，但同时规 定，经分析论证，也可采用其他线型。本文 在文献《广义r分布的特性和应用〉〉（金光 炎、董秀颖,2003）研究的基础上，尝试应 用指数r分布线型和对数 r分布线型拟合 沙坪水闸年最高水位资料系列进行频率特 性的统计分析。 1.沙坪水闸年最高水位的描述性统计 分析 沙坪水闸位于西江下游右岸沙坪河入 江河口处，见图1。该水闸是广东省鹤山市 的一座以防洪、排涝为主，结合蓄水灌溉和 改善航运的综合运用的中型水闸。该闸有从 1964年自建闸以来至2006年43年的闸外年 最局水位观测资料系列，见表 1。资料系列 的最小值、最大值分别为和，为此将年最高 水位系列分为11个组限，其中组距宽度hi 为的9组,其他组距2组。计算在每个组限 内的年最高水位出现的频数 fi ,相应的年最 高水位各组距区间上的经验频率为 fi/n(n=43,i=1,2, …,11)。沙坪水闸年最高 水位经验频率分组计算详见表 3。在各组距 区间上作以平均频率密度 fi/(nhi) 为高的 直方图，见图4。从图4可看出直方图呈近似 的正态分布形态，它有一个中间主峰，两头 低，主峰在区间v xw上；年最高水位平均 值、中位数都落在区间＜ xw上，经计算平均 值为，中位数为。用矩法计算年最高水位经 验系列统计量分别为：标准差b 12n\变差系数Cv=，偏态系数Cs=-, 峰态系数Ce=-。 鹤山市 鹤山市 图1沙坪水闸地理位置示意图 2. P- :rn、指数r分布和对数 r分布 线型比较分析 《水文分析与计算》已有的模比系数Kp 值表通过变量坐标的平移转换得到适线结 果。四参数指数r分布和对数r分布线型 都比式⑴的P-m线型多了 1个参数，因此适 线弹性要比P-m线型大。在当今计算机技术 被广泛应用于各领域的条件下，就计算方法 难易程度与精度比较而言与 P-m分布线型 相当。本例适线结果如图 2、3所示，点绘 在概率格纸上的经验频率点据上部呈现出 向下凹的分布 形态。图中点据旁标注数字是水位发生 年份。 图3沙坪水闸年最高水位-频率曲线指数r 分布、对数r分布线型比较图 P- 口、指数r分布和对数 r分布线型 拟合优劣分析及统计判别 频率线型拟合优劣分析可分为：一是从 适用范围上对频率线型是否与给定区域内 的所有水文数据系列拟合得好进行分析， 这 需要对多站数据系列进行分析；二是对单站 数据样本拟合优劣进行分析。本例仅就沙坪 水闸年最高水位单站观测数据样本拟合优 劣进行分析，因此可以应用统计学的距离分 析和统计判别方法。统计判别分析是根据事 物特点的变量值和它们所属的类求出判别 函数，依据判别函数对未知所属类别的事物 进行分类的一种分析方法。常用的判别方法 有距离判别、贝叶斯判别、费歇尔判别以及 逐步判别。而距离判别法较为直观，适用面 广，对各类(或总体)的分布无特定要求，因 此本文采用距离判别法。 多个总体的距离判别是计算样本 x到每个总 体的距离d2(x C Gi,i=1,2, ???*),然后比较 这些距离，如x到总体Gi的距离最短，贝U 判x属于总体Gi。明氏距离特别是其中的欧 氏距离是人们较为熟悉的，也是使用最多的 距离。在本例计算 序号 按时间序列排序 年份水位a 按水位大小排序 年份水位a 频率b p m (%) p-m线型 拟合值a离差c 指数r线型 拟合值a离差c 对数r线型 拟合值a离差 C119641994---219652005--319661968--41 9671998519681976--619691997-719701978 8197119889197219741019731966111974197 525-1219751973--1319761981--141977198 2-15197819921619791979171980197018198 119831919822001---2019831996---211984 1977---221985196450---2319861993---24--2519882006---261989197