PⅢ型分布在设计潮位中应用.docVIP

PⅢ型分布在设计潮位中应用 　　摘要以实测潮位资料为依据,对海河口潮位进行数理统计分析,按P-Ⅲ型分布和Gumbel型分布给出了不同频率设计潮位,同时对2种分布型的计算精度进行统计研究。结果表明:P-Ⅲ型设计拟合效果较优,对设计潮位的分析计算具有重要的参考价值。 　　关键词P-Ⅲ型分布;Gumbel型分布;精度;潮位 　　 　　海河干流是向海河流域供水的重要入海通道,河口位于天津市滨海新区,为感潮河道,汛期行洪易受渤海高潮位顶托。确定合理的河口设计潮位对海河干流治理和天津市防洪规划具有重要的意义。 　　1潮位频率计算方法 　　潮位变化复杂,影响潮位变化的因素很多。目前仍采用数理统计方法进行潮位分析,国内外学者的研究成果表明,由实测高潮位组成的样本,各元素是相互独立的随机变量[1-2]。随着样本元素增多,当元素个数n30时,统计特征趋于稳定。潮位分布线型主要有P-Ⅲ型分布和Gumbel分布2种形式[3]。 　　1.1P-Ⅲ型分布 　　分布密度函数为: 　　f(x)=(x-α0)α-1e(1) 　　式中:α=;β=;α0=EX(1-)。 　　1.2Gumbel型分布 　　分布密度函数为: 　　f(x)=αexp[-α(x-)-e](2) 　　式中:α为参数;为众值;EX=+;σ=;CS=1.139 55。 　　2海河口设计潮位计算 　　2.1资料的选择整理 　　海河口曾有多处验潮站,但多数站点的观测年限较短,且资料不完整,计算选用观测年限在30年以上且资料完整的六米和海河闸2个验潮站的1959―1992年最高潮位及汛期最高潮位资料,同时对北炮台站年最高潮位和汛期最高潮位资料进行分析。根据历年地面沉降观测值对1959年以来的潮位观测值进行修正,测站基本情况详见表1。 　　2.2高潮位重现期的确定 　　在六米和海河闸站的观测资料中,1992年的年最高潮位分别达到4.83 m,较一般年份高出很多,应做特大值处理。在近百年中1885年曾出现过1次强潮,但具体潮位值难以确定。1939年北炮台站观测值为4.63 m,但北炮台站处在感潮河段(当时未建海河闸),其多年平均值较海河闸站低0.22 m,所以确定1992年高潮位重现期约为50~100年。 　　2.3参数估算 　　均值(数学期望)采用矩法计算: 　　x=XN+Xi(3) 　　CV值采用矩法估算: 　　CV= (4) 　　Cs应用配线法估计。 　　经验频率计算公式为: 　　Pm=(5) 　　分布参数详见表2。 　　2.4设计潮位成果βα 　　采用P-Ⅲ型和Gumbel型2种分布型式分别分析出六米站、海河闸和北炮台站不同设计频率相应的设计潮位,成果详见表3。 　　由拟合频率曲线可以看出,P-Ⅲ型与Gumbel型计算出的结果相近,但拟合成果P-Ⅲ型优于Gumbel型。所以推荐采用P-Ⅲ型分布推求设计潮位。 　　3成果的可靠性分析 　　3.1假设检验 　　参数估算后,对所选用的分布线型是否合理进行分析。 　　3.1.1显著性检验。由于北炮台、六米站及海河闸的距离很近,对是否可将北炮台与六米站或海河闸做为同一系列进行计算,这里应用t检验法进行分析[4]。 　　由计算成果可知,北炮台站年潮位系列参数: 　　X=3.78,Cv=0.07,S=0.26,n=53。 　　六米站年潮位系列参数:X1=3.95,Cv1=0.08,S1=0.316,n1=34。 　　海河闸年潮位系列参数:X2=4.00,Cv2=0.074,S2=0.31,n2=35。 　　则统计量: 　　T北炮台~六米站=#8226;=2.67 　　T北炮台~海河闸=#8226;=3.55 　　由自由度n1+n-2=85取显著水平α=0.05得t=1.98,则P(T北炮台~六米站1.98)=0.95;由自由度n2+n-2=86取显著水平α=0.05得t=1.98,则P(T北炮台~海河闸1.98)=0.95。所以,前后两段的潮位是取自不同的系列,经检验可知,北炮台与六米站或与海河闸的潮位是相互独立的,来自2个总体的2个样本,故不能将北炮台与六米站或与海河闸组合计算。 　　3.1.2参数检验。①数学期望(均值):按照上文2.3中的(3)式计算。②方差σ2的检验:X2检验法。已知:总体方差S=0.316,样本方差σ0=0.28,则统计量ns2/σ02=43.3。取显著水平α=0.05,由自由度n-1=33及所对应的2个临界值得: 　　X=50.73及X=19.05, 　　PX≤≤X=1-a=0.95 　　说明方差一致。即:样本参数与总体是一致的,具有较好的代表性。 　　3.1.3非参数的检验。在P-Ⅲ型分布假设下,六米站年最高潮位系列3个参数估计为: