水利工程论文-关于黄河三门峡年径流量频率线型分析.doc

水利工程论文-关于黄河三门峡年径流量频率线型分析论文关键词：年径流量统计分析频率线型三门峡黄河论文摘要：通过三门峡1919至1989年71a年径流量的统计分析，对P-、正态分布线型和对数分布线型与经验频率点据拟合结果进行了比较。作为一种尝试，采用四参数对数分布线型拟合三门峡年径流量资料系列进行频率特性分析，取得了较为合理的分析研究成果。研究结果表明四参数对数分布线型优于P-和正态分布线型，验证了三门峡年径流量频率分布服从四参数对数分布线型。年径流量统计分析是水资源科学利用和水利电力工程规划设计的重要依据。在水文频率计算中,规范1推荐采用P-线型,但同时规定,经分析论证,也可采用其他线型。笔者在文献广义分布的特性和应用(金光炎、董秀颖,2003)2研究的基础上，验证了沙坪站年最大日雨量、白蕉站年最低潮位以及沙坪水闸年最高水位观测资料系列服从四参数对数分布345，本文尝试应用四参数对数分布线型拟合三门峡年径流量资料系列进行频率特性的统计分析。1.本文采用的三门峡年径流量实测资料概况三门峡1919至1989年71a年径流量实测资料系列和推求的14701918年年径流量系列出自文献6，实测资料系列见表1。实测资料系列可通过统计学游程检验满足样本随机性要求。实测资料系列的最小值、最大值分别为241.4亿m3和802.6亿m3,极差561.2亿m3。为此将年径流量实测资料系列分为12个组限,其中组距宽度hi为61.5亿m3的10组,其他组距2组。计算在每个组限内的年径流量出现的频数fi，图1。相应的年径流量各组距区间上的经验频率为fi/n(n=71,i=1,2,12)。三门峡1919至1989年71a年径流量实测资料系列经验频率分组计算详见表3。在各组距区间上作以平均频率密度fi/(nhi)为高的直方图,见图3。从图1、3可看出实测资料系列直方图呈近似的正态分布形态,采用单样本K-S检验可验证其近似服从正态分布，图1中的曲线为按正态分布拟合。实测资料系列直方图有一个中间偏左主峰,两头低。中位数为491.2亿m3，在主峰区间445.95x507.45亿m3上；年径流量平均值为510.93亿m3,在主峰右侧区间507.45x568.95亿m3上。用矩法计算三门峡年径流量实测资料系列统计量分别为:标准差=114.92亿m3，变差系数Cv=0.225，偏态系数Cs=0.258，峰态系数Ce=-0.045。(a)实测系列(b)推求系列图1三门峡年径流量资料系列直方图2.P-、正态分布和对数分布线型比较分析水文分析与计算(刘光文,1963)中提出了径流频率分析适线的线型选择原则7:在计算简便的同时,具有尽量高的精度和弹性；曲线与经验频率点据得到最好的拟合；曲线的形状大致符合水文现象的一定物理性质,如曲线应该有一定的极限,不该出现负特征值。其中第点的实质是对资料系列进行统计判别，第点的实质是合理性要求。水文频率分析述评(金光炎,1999)也认为按照水文物理概念,曲线应该有上限,并对Slade(1936)、谢家泽(1958)等人的研究观点作了介绍8。文献9、10也认为现行的频率线型在合理性方面需要完善。因此本文在上述原则的基础上对备选的3种频率线型进行比较。2.1P-和对数分布线型计算方法在采用P-、四参数对数分布线型频率累积曲线函数拟合适线过程时，由于参数都为非线性关系形式，不能通过某种转换变为线性形式，因此只能采用非线性迭代回归的办法求解。进行非线性迭代回归时，首先确定频率累积曲线函数的表达式，确定参数的初始值，然后根据某种方法进行搜索迭代，反复调整初始值，按规范应用最小二乘法原理使得观测值与拟合值的离差平方和最小时(或者结合其他一些条件)结束迭代过程，得到各参数的最后计算结果。四参数对数分布线型比P-线型多了1个参数，因此适线弹性要比P-线型大。在当今计算机技术被广泛应用于各领域的条件下,就计算方法难易程度与精度比较而言与P-分布线型相当。本例适线结果如图2所示，点绘在概率格纸上的经验频率点据接近直线分布形态。图中点据旁标注数字是径流量对应的年份。图2三门峡年径流量-频率曲线线型比较图2.2P-、正态分布和对数分布线型拟合优劣的统计判别本例适线计算结果见表1。表1最后一行离差的SUMSQ值是观测值点据与序号按年份排序年份径流量按径流量排序年份径流量频率bpm(%)P-线型拟合值a离差c正态线型拟合值a离差c对数线型拟合值a离差c11919415.71967802.61.389817.34-14.74773.5329.07794.278.3321920510.91964799.82.778770.4929.31739.4160.39755.6044.2031921523.819377064.167741.46-35.46717.59-11.59730.82-24.8241922390.219836965.556719.95-23.95701.07-5.07712.07-16.0751923421.31949679.26.944702.66-23.46687.55-8.35696.76-17.5661924306.119586678.333688.07-21.07675.98-8.98683.68-16.6871925430.31968661.29.722675.38-14.18665.78-4.58672.18-10.9881926315.41989654.311.111664.09-9.79656.60-2.30661.85-7.5591927403.81940645.912.500653.88-7.98648.21-2.31652.44-6.54101928241.41975645.613.889644.531.07640.455.15643.761.84111929347.21976644.615.278635.878.73633.2111.39635.688.92121930364.41981627.916.667627.790.11626.381.52628.09-0.19131931336.51946626.218.056620.206.00619.926.28620.935.27141932321.41961621.319.444613.028.28613.767.54614.12781935619.420.833606.2013.20607.8611.54607.6211.78161934479.81938618.322.222599.6918.61602.1916.11601.3916.9111935619.41961623.59323.59620.59521.743.9611.4545.7218.4050181936455.71984616.325.000587.4528.85591.4224.88589.6226.681919377061955607.126.389581.6625.44586.2820.82584.0323.07201938618.31954575.927.778576.07-0.17581.28-5.38578.60-2.70211939428.81963568.929.167570.64-1.74576.39-7.49573.32-4.42221940645.91985568.830.556565.363.44571.61-2.81568.180.62231941368.51966564.931.944560.234.67566.93-2.03563.151.75241942411.71959551.233.333555.21-4.01562.33-11.13558.23-7.03251943616.91933541.834.722550.31-8.51557.81-16.01553.41-11.61261944477.21945541.636.111545.50-3.90553.35-11.75548.68-7.08271945541.61951535.437.500540.79-5.39548.96-13.56544.02-8.62281946626.21982529.238.889536.16-6.96544.61-15.41539.43-10.23291947492.31978524.240.278531.60-7.40540.31-16.11534.91-10.71301948480.91921523.841.667527.10-3.30536.04-12.24530.45-6.65311949679.21988519.543.056522.67-3.17531.81-12.31526.03-6.53321950511.41979514.644.444518.29-3.69527.60-13.00521.66-7.06331951535.41950511.445.833513.95-2.55523.42-12.02517.33-5.93341952491.21920510.947.222509.651.25519.24-8.34513.03-271947492.348.611505.39-13.09515