洪水频率计算(规范方法).docxVIP

PAGE PAGE # 附录A 洪水频率计算 A1洪水频率曲线统计参数的估计和确定 A1.1 参数估计法 A1.1.1矩法。对于n年连序系列，可采用下列公式计算各统计参数 均值均万差或变差系数偏态系数或式中lUi-X)2 均值 均万差 或 变差系数 偏态系数 或 式中 lUi-X)2 n-1 二 Xi2-n([Xi)2 n7 (Xi - X)3 i ￡ (n —1)( n —2)X3C； n n n n n2 v X； _3 n^ Xi X2 2(^ X J3 i # i￡ i住 i仝 ——3 n(n -1)( n - 2) X C： Xi 系列变量（i=1，…，n）； (A1) (A2) (A3) (A4) n 系列项数。 对于不连序系列，其统计参数的计算与连序系列的计算公式有所不同。 如果 在迄今的N年中已查明有a个特大洪水（其中有I个发生在n年实测或插补系列 中），假定（n-l）年系列的均值和均方差与除去特大洪水后的（N-a）年系列的 相等，即Xn』= Xn4,Sn』=Sn4，可推导出统计参数的计算公式如下： — 1 a N — a n X 二丄C Xj Xi） （A5） nJ j n — I 4 Cs =N_1 一)2 Cs = N_1 一)2 N JXj—X)3 活二Xi -对 (A6) (A7) (N -1)( N _2)X Cv 式中 Xj 特大洪水变量（j=1，…，a）; Xi――实测洪水变量（i=l+1,…，n）o A1.1.2概率权重矩法。概率权重矩定义为 1 .M 1 . M . = oxFJ(x)dF j=0,1,2,… (A8) 皮尔逊川型频率曲线的三个统计参数不能用概率权重矩的显式表达。但经 推导有： CvH(1M2 -M0/3M^ M0/2式中，H和R都和Cs有关，并已有近似的经验关系如下: Cv H( 1 M2 -M0/3 M^ M0/2 式中，H和R都和Cs有关，并已有近似的经验关系如下: 广 2 3 4 Cs =16.41u-13.51U 州0.72u +94.54U R—1 (4/3-R)012 2 3 4 H =3.545+29.857 —29.15V +363.8V +6093V (R -1)2 (4/3 -R)0.14 (1 R ： 4) 3 (A9) (A10) (A11) (A12) (A13) 为保证Cv和Cs有二位小数准确，要求在用式（A11）计算R时，M。、M1 和M2的计算值至少达到5位有效数字。 1根据连序系列计算概率权重矩。将洪水系列按从大到小顺序排列，样本 概率权重矩按下式计算: M o 彳 n 七 Xi n住 M 1 彳 n 七 X i n -i n住 n —1 彳 n 一1 y Y (n — i)( n -i -1) Ivl 2 — 厶 X i n住 (n —1)( n — 2) (A14) 2根据含历史洪水特大值的不连序样本计算概率权重矩。 才C1(N _ j)(n _ j -1)N -a n _ IX j C2N _a n —IXi(A15)(N _1)(N _2)1 -N-n I(n —I —1) (n —I —2)式中，C1，C2都是对不连序系列中实测洪水概率权重的修正系数。N _a +1 C1N 1「N _a +1 YN 1(A16)C2 才C1 (N _ j)(n _ j -1) N -a n _ I X j C2 N _a n —I Xi (A15) (N _1)(N _2) 1 -N - n I (n —I —1) (n —I —2) 式中，C1，C2都是对不连序系列中实测洪水概率权重的修正系数。 N _a +1 C1 N 1 「N _a +1 Y N 1 (A16) C2 A1.1.3双权函数法。均值仍用矩法，如式（A1）计算。而Cv和Cs的计算公式 丄 _-E1_ C2 _ hX k2H1 v = △.旦 -D1 H 1 (A17) Cs —(Xcv △-丄) Cv D1 h (A18) 式中，k、h是待优选的系数，可采用未加权的、数值积分计算的 Cv,按下式选 疋：h=Cv， K=1/Cv。 qQ E1 二(X - X)住 1(x) f (x)dx (A19) 第一权函数 第二权函数 积分式 二 (X -乂)2「1(x) f (x)dx A1 = . 1 (x) f (x)dx D! =「(X —X)- !(x) f (x)dx G(x) ■ 1(x) (A19) (A20) (A21) (A22) -k exD—k2(X -X)2 _ X 2 二 e x P = exp-吨垃 X (A23) ——2 2X (A24) ~式（A22 ）可用数值积分公式计算。例如，当 n为奇数时, 采用权积分系数：8，-4, 8，1, 4