经验频率曲线计算参数
4第四章 水文统计基础知识
P-Ⅲ型概率密度曲线的 特点：
（1）单峰型； （2）与x轴有一交点，对应水文 变量的最小值； （3）后端与x轴不相交
P-Ⅲ 型曲线的应用
将P-Ⅲ型曲线的方程式进行一定的积分演算，就可以得到频率曲 线纵坐标值 的计算公式，即频率曲线的方程式（分布函数）为： P
x
xP (Cv 1) x K P x
3、偏态系数
偏态系数是反映随机变量系列中各随机变量对其均 值对称性的参数。 对于总体
Cs
Cs
( xi x) 3
i 1
n
n x C v3
3
对于样本
Cs
( xi x) 3
i 1
n
(n 3) x Cv3
3
频率曲线的三个参数，其中均值（ x ）一般直接采用矩 法计算值；变差系数（Cv）可先用矩法估算，并根据适线拟 合最优的准则进行调整；偏态系数（Cs）一般不进行计算， 而直接采用倍比，我国绝大多数河流可采用 Cs=(2~3)Cv。
Ki
2）.中值 x
xi x
中值的大小能反映系列中间项和密度曲线的位置。
3）.众值
x
众值的大小能反映系列中最大几率项和密度曲线的位置。
4）.均值、中值、众值的位置关系
y y
y
o
xxx a)
x 0
xxx b)
x 0
xxx
c)
x
a）正偏态；
b）正态； 密度曲线图
c）负偏态
水文现象为不对称分布，年洪峰流量频率分布多为正偏。
频率曲线（ P-Ⅲ 型曲线） 设计洪水水位 设计洪水流量
工程设计标准
设计洪水频率 （洪水重现期） 经验频率曲线
公式4-18~公式4-20
随机变量的统计参数资料
解：尽管甲 乙, 但由于x甲 x乙， 所以不能用来说明离散程度的大小，
而必须用相对值Cv来衡量，即：
Cv甲
甲
x甲
360 1200
0.3
Cv乙
乙
x乙
320 800
0.4
这说明甲地区的年雨量离散程度比乙地区的小。
3、偏态系数CS 衡量随机变量分布在均值两边是否对
称以及不对称（偏态）程度的参数。
n
二、矩法
§4—4 频率计算 一、经验频率曲线 （一）经验频率的计算公式
p m 100% n
经验频率的修正公式
p m 100% n 1
（二）经验频率曲线的绘制及其应用
1、绘制经验频率曲线
（1）排序（由大到小）；
（2）计算经验频率（
p
m n 1
100%
）；
（3）在机率格纸上点绘各相应点据（pi , xi）；
三、随机变量的统计参数
1、算术平均数 x
x
x1 x2
x3 xn n
1 n
n i 1
xi
均值对密度曲线的影响
f(x)
1
2
o
x1
x
x2
均值对频率曲线的影响
x
x1 x2
2 1
频率p（%）
2、均方差σ与变差系数CV 均值相等的不同系列，它们的离散程
度可用均方差σ来衡量。
n
(xi x)2
i1
四、抽样误差
抽样误差：由随机抽样而引起的误差。 均方误计算公式：
x
n
cv
Cv 2n
1
2Cv 2
3 4
Cs 2
2CvCs
Cs
6 n
水文分析计算课程设计
《水文分析与计算》课程设计指导书———设计年径流及设计洪水的计算一、课程设计的目的1．掌握PIII型频率曲线的制作方法2. 掌握设计年径流及其年内分配的计算方法3．掌握考虑历史特大洪水的设计洪水及其过程的计算方法二、课程设计任务1．根据所给资料推求设计年径流与设计年内分配过程表1是某站1958～1976年各月径流量资料，根据所给资料推求P=10%的设计丰水年、P=50%的设计平水年、P=90%的设计枯水年的设计年径流量；并计算P=90%的设计枯水年径流年内分配过程。
要求：理论频率曲线采用PIII型分布，由矩法作参数无偏估计，并以估计值为初值，用目估适线法选配理想的理论频率曲线，注意比较验证均值X a、变差系数C V、偏态系数C S对频率曲线的影响效果。
检查所选最终的理论频率曲线的合理性，并计算所求设计频率的相应设计年径流，年径流分配过程采用典型年同倍比放大法。
3三、课程设计成果要求要求提交设计成果：一份电子文档，一份打印文档。
设计中的计算可采用采用excel 或编程计算，编程语言可采用FORTRAN 语言、C 语言、Basic 语言或同等功能的语言编程。
要求程序正确、可靠、可运行，符合结构化程序设计思想，具有易读性、可修改性、可验证性、通用性，关键变量应作注释说明。
计算结果要表格化，便于检查、保存和打印。
设计设计报告，其重点是对计算成果的说明和合理性分析及其有关问题的讨论。
要求文字流畅，简明扼要；图表整齐清楚，名称、编号齐全；封面统一，最后装订成册。
四、课程设计的考核平日考勤、设计报告，加上抽查提问及上机操作，对成绩进行综合评定。
五、课程设计时间与地点时间： 2013年5月9日星期四 地点： 学院六、实验原理1.经验频率计算经验频率：P=m/(n+1)*100%，模比系数：Q Q Ki i = 2．线型选择频率曲线一般应采用皮尔逊Ⅲ型。
3.频率曲线参数估计平均值：n1∑==ni iQQ变差系数：()1n 112--=∑=ni iv K C4.偏态系数：Cs=2-3Cv七、实验步骤1、将测站所得数据年份及年平均流量数据复制与Excel 表格中，并列出序号，同时计算出年平均流量的均值。
公路桥涵水文计算基本方法
第二节 水文勘测分析计算基本途径
•
桥涵水文计算、分析基本途径如下：
•
1， 有水文观测资料—— —— 水文统计法
•
2， 无水文观测资料—— --- 形态断面法
•
3， 无水文观测资料（无居民)—经验公式法
• 一， 有水文系列观测资料时水文统计法：
• (一)，资料搜集和准备：
• 1，外业勘测前的准备工作
桥涵水文与本行业其它专业有所不同，桥涵水文调查、分析、计算本 身并无精度指标要求，特别强调的是将通过各种途径和方法得到的计算 结果进行比较、论证后确定最终设计流量，使其更接近实际，更趋于合 理。本次交流着重于桥涵水文分析、计算的基本方法和途径。
以上内容是桥涵水文工作主要工作内容，但重中之重是设计流量的推 算。至于桥长、冲刷、调治构造物、桥面标高计算相比之下要简单得多， 因此，本次交流的重点放在设计流量的推算、外业调查、勘测的主要过 程以及内业工作的主要内容和步骤。
系列流量。上表括号内（流量y）为插补后分析站流量y的系
列流量,插补延长所得资料不宜用于第三站，可能引起较大误差。
4400 4000 3600
分 3200 析 2800 站 2400 流 2000 量 1600
1200 800 400
相关分析方程图像
方程：Y=0.63x+98.04
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
* *
*
500 1000 1500 2000 2500 3000 3500 4000 4500 5000 5500 6000
参证站流量
示例 2,水位、流量关系曲线法示例： 水位、流量关系曲线法就是利用两系列的水位、流量对应关系曲线,直接对分
工程1_水文频率分析计算(P-Ⅲ型曲线)
计算说明书━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━工程名称：工程1计算类型：水文频率分析计算(P-Ⅲ型曲线)一、计算原理1.适用范围本程序可一次完成一个水文系列频率计算的全部工作，对连续系列和不连续系列均为适用。
本程序完成的工作内容包括：系列排队、计算经验频率及统计参数值、通过优选P-Ⅲ型曲线的参数Cv、Cs值进行适线或用目估法适线、绘制频率曲线图、计算所采用的频率曲线的各设计频率下的设计值等。
为满足工程的实际需要，本程序除可用优选统计参数的方法适线外，还可用目估适线法进行适线。
因为本程序在用优选法适线时，对各经验点据是给以等权重的处理。
而当需要对各点据给以非等权重的处理时（如：设计洪水中要求多照顾首几项洪水；在年径流计算时要多照顾末端；或由于基本资料精度差等），单用优选法就不合适，此时可改用目估适线法。
为了减少目估适线时的盲目性，实际使用时，一般采用优选与目估适线相结合的方法，即先用优选法选出一条通过点群中心的频率曲线。
在此基础上再用目估的方法对优选出的参数Cv、Cs做少许调整，重新适线，以达到对各点据给以不同权重的目的，获得满意的结果。
2.计算方法和公式3.规范规程(1)《水利水电工程设计洪水计算规范》（SL44-2006）(2)《水利水电工程水文计算规范》（SL278-2002）4.参考文献(1)《水利水电工程设计洪水计算手册》水利部长江水利委员会水文局水利部南京水文水资源研究所主编,1995年10月(2)《工程水文学》(第三版)武汉大学叶守泽,河海大学詹道江合编,中国水利水电出版社,2000年10月(3)《最优化理论与算法》(第二版)陈宝林编著,清华大学出版社,2005年10月(4)《水利水电工程设计计算程序集 A-3 水文频率计算程序》作者马明（新疆水利水电勘测设计院）(5)《Visual Basic常用数值算法集》何光渝编著,科学出版社,2002年(6)《科学与工程数值算法〔Visual Basic版〕》周长发编著,清华大学出版社,2002年二、基本数据连续系列项数n=33序位系列值编号1 1952 114.002 1953 118.003 1954 116.004 1955 105.005 1956 122.006 1957 88.807 1958 141.008 1959 132.009 1960 107.0010 1961 94.8011 1962 94.0012 1963 113.0013 1964 114.0014 1965 101.0015 1966 104.0016 1967 92.8017 1968 97.1018 1969 116.0019 1970 122.0020 1971 145.0021 1972 119.0022 1973 111.0023 1974 83.1024 1975 93.5025 1976 104.0026 1977 88.5027 1978 95.3028 1979 92.5029 1980 115.0030 1981 94.5031 1982 107.0032 1983 90.9033 1984 89.10三、计算结果1.统计参数值：均值 Xa=106.694均方差 S=15.308变差系数Cv=0.143偏态系数Cs=0.708Cs/Cv=4.934经验频率值表序位系列值频率(%)编号1 1971 145.000 2.9412 1958 141.000 5.8823 1959 132.000 8.8244 1970 122.000 11.7655 1956 122.000 14.7066 1972 119.000 17.6477 1953 118.000 20.5888 1969 116.000 23.5299 1954 116.000 26.47110 1980 115.000 29.41211 1964 114.000 32.35312 1952 114.000 35.29413 1963 113.000 38.23514 1973 111.000 41.17615 1960 107.000 44.11816 1982 107.000 47.05917 1955 105.000 50.00018 1966 104.000 52.94119 1976 104.000 55.88220 1965 101.000 58.82421 1968 97.100 61.76522 1978 95.300 64.70623 1961 94.800 67.64724 1981 94.500 70.58825 1962 94.000 73.52926 1975 93.500 76.47127 1967 92.800 79.41228 1979 92.500 82.35329 1983 90.900 85.29430 1984 89.100 88.23531 1957 88.800 91.17632 1977 88.500 94.11833 1974 83.100 97.0592.优选P-Ⅲ型曲线的参数Cv、Cs值(离差平方和准则)：离差平方和S=148.063Xa=107.218Cv=0.157Cs=1.053Cs/Cv=6.7113.理论频率曲线设计值(目估适线)：Xa=107.218Cv=0.157Cs=1.053Cs/Cv=6.711理论频率曲线设计值表频率P(%) 模比系数Kp 设计值Xp B 绝对误差δXp 相对误差δ'Xp(%)0.01 1.954 209.467 11.46 33.551 16.017 0.1 1.723 184.710 8.09 23.694 12.828 0.2 1.651 177.039 7.06 20.665 11.6730.5 1.555 166.681 5.71 16.734 10.0401 1.480 158.636 4.72 13.822 8.7132 1.402 150.357 3.77 11.053 7.3513 1.356 145.388 3.13 9.166 6.3045 1.296 138.942 2.74 8.038 5.78510 1.210 129.76620 1.118 119.84630 1.058 113.48140 1.012 108.54250 0.973 104.32160 0.937 100.48370 0.903 96.77280 0.867 92.92090 0.825 88.42995 0.796 85.39297 0.781 83.71999 0.757 81.12299.9 0.728 78.104四、图形结果────────────────────────────────────────────────────────计算软件：SGGH-Tools 2011 计算者：校核者：计算日期：2020/8/24。
3.5 理论频率曲线
环境与市政工程学院水文学第3章 水文学统计基本原理与方法本章提纲3.1 水文统计概述3.2 频率与概率3.3 经验频率曲线3.4 随机变量的统计参数3.5 理论频率曲线3.6 频率曲线的分析与抽样误差3.7 相关分析3.5.1 皮尔逊III型曲线（1）曲线方程 英国生物学家皮尔逊在统计分析了大量随机现象后，发现概率密度曲线均为类似于铃形的曲线，l 只有一个众数，在众数处曲线的斜率等于零，若把纵坐标移到均值处，即当：l 曲线的两端或一端以横轴为渐近线 ，即当：0d d ,=-=xy d x 0d d ,0==xy y 这种曲线有两个特点：3.5.1 皮尔逊III型曲线（1）曲线方程式中：a 0，系列起点到坐标原点的距离；α，代换参数，α－1＝a /d ；β，代换参数，β＝1/d ；e ，自然对数的底。
())(100e )(a x a x y ----Γ=βαααβØ 根据上述两点，皮尔逊建立了概率密度曲线微分方程，经过参数代换、分离变量积分等得出皮尔逊III 型曲线方程：Ø也可将这些待定参数用统计参数表示，代入皮尔逊III 型曲线的方程式中，则方程式可以写成：),,,(x C C x f y s v =3.5.1 皮尔逊III型曲线（2）皮尔逊III型曲线的绘制 对密度函数进行积分得：()x a x x x P p x a x p d e )()()(100⎰∞----Γ=≥βαααβ 为简化计算，引入离均系数，进行参数代换，制成数表便于查用vC x x x -=Φ在频率计算时，先由已知C s 值查Φ值表，得到不同频率P 的离均系数Φp 值，然后将Φp 值及已知的 、C v 代入下式，即可得到相应的水文特征值：()xC Φx C ΦK v P P v P P 11+=+=x3.5 理论频率曲线3.5.1 皮尔逊III型曲线（2）皮尔逊III型曲线的绘制绘制频率曲线步骤：•由实测资料，统计计算 、C v ；•确定C s ；•由C s 查附录3，得不同频率P 的离均系数Φp 值；•由Φp C v ＋1＝K p ，求K p ；•求不同频率P 对应的x p ，在海森概率置上绘制理论点据(P , x p )；•根据理论点据分布趋势，目估并绘制一条光滑曲线。
(完整版)第二节设计洪峰流量及设计洪量的推求
第二节设计洪峰流量及设计洪量的推求由流量资料推求设计洪峰及不同时段的设计洪量，可以使用数理统计方法，计算符合设计标准的数值，一般称为洪水频率计算。
一、资料审查在应用资料之前，首先要对原始水文资料进行审查，洪水资料必须可靠，具有必要的精度，而且，具备频率分析所必须的某些统计特性，例如洪水系列中各项洪水相互独立，且服从同一分布等。
除在第三章谈到审查资料的可靠性之外，还要审查资料的一致性和代表性。
为使洪水资料具有一致性，要在调查观测期中，洪水形成条件相同，当使用的洪水资料受人类活动如修建水工建筑物、整治河道等的影响有明显变化时，应进行还原计算，使洪水资料换算到天然状态的基础上。
洪水资料的代表性，反映在样本系列能否代表总体的统计特性，而洪水的总体又难获得。
一般认为，资料年限较长，并能包括大、中、小等各种洪水年份，则代表性较好。
此可见，通过古洪水研究，历史洪水调查，考证历史文献和系列插补延长等增加洪水列的信息量方法，是提高洪水系列代表性的基本途径。
根据我国现有水文观测资料情况，SL44—93规定坝址或其上下游具有较长期的实测水资料(一般需要30年以上)，并有历史洪水调查和考证资料时，可用频率分析法计算计洪水。
二、样本选取河流上一年内要发生多次洪水，每次洪水具有不同历时的流量变化过程，如何从历洪水系列资料中选取表征洪水特征值的样本，是洪水频率计算的首要问题。
根据SL44—93规定，应采用年最大值原则选取洪水系列，即从资料中逐年选取一个大流量和固定时段的最大洪水总量，组成洪峰流量和洪量系列。
固定时段一般采用l、3、5、7、15、30天。
大流域、调洪能力大的工程，设计时段可以取得长一些；小流域、调洪能力小的工程，可以取得短一些。
在设计时段以内，还必须确定一些控制时段，即洪水过程对工程调洪后果起控制作用的时段，这些控制时段洪量应具有相同的设计频率。
同一年内所选取的控制时段洪量，可发生在同一次洪水中，也可不发生在同一次洪水中，关键是选取其最大值。
水文频率计算方法
x1x3
13
（3-25）
Sf(Cs,P)
S表
附录C P-III曲线三点法 S 与 Cs 关系表
P368
P = 5 ~ 50 ~ 95%
求参步骤
√ 三点法：
① 据(3-25)式左端计算得S值，且已知P1,P2,P3,查 附录C,求得参数Cs；
② 据Cs查附录B
(3-24)式,求得
。
x,Cv ,Cs
2. 优化适线法
最小二乘估计法：
n
S(){ [xi f(Pi,)]2} i1
S(ˆ)minS()
S 0
(3-26) (3-27) (3-28)
End
水文频率计算方法
频率计算推求 xp
参数初估
适线法
一. 统计参数初估方法 1 矩法 2 经验关系法 3 三点法 4 权函数法 5 概率权重矩法
矩法
1.矩法
依据实测系列计算三个统计参数：x,Cv ,Cs
经验法
2. 经验关系法
x
1 n
n
xi
i1
Cv
(K i 1) 2
n 1
Cs初值：
设C 计s 暴 雨量( n
x
2
x (1 C v 2 )
x 3 x (1 C v 3 )
（3-22）
三点的取法： 1～59% 或 3～50～97% 或 5～50～95% 或 10～50～90%
参数计算公式
x 1x3 3x1 1 3
Cv
x1 1x3
x3 3x1
（3-23） （3-24）
x1x32x2 1322 =S 偏度系数
计算表 题 P66
p m 100% n 1
年最大洪峰流量 Q（m3/s)
水文频率适线法
三、适线法
适线法旳原理：根据经验频率点据，找 出配合最佳之频率曲线，相应旳分布参 数为总体分布参数旳估计值。
1、计算环节：
（1）点绘经验点据： 纵坐标为变量值，横坐标为经验 频率，采用期望值公式估计。
思索题
1、用配线法进行频率计算时，判断配线是否 良好所遵照旳原则是[_c___]。
a、抽样误差最小旳原则； b、统计参数误差最小旳原则； c、理论频率曲线与经验频率点据配合最佳旳
原则； d、设计值偏于安全旳原则。
是非题
1水文频率计算中配线时，增大Cv能够使频率 曲线变陡。（ √ ）
2 给经验频率点据选配一条理论频率曲线，目 旳之一是便于频率曲线旳外延。（ √ ）
第六节 水文频率计算适线法
x
1200
1000
W（X≥xi）＝i / n
800
0
20
40
60
80
100 W（%）
某地年降雨量经验分布曲线
二、经验频率
假如用P（X≥xi）＝m/n 旳经验分布曲
线估计总体分布曲线，存在不合理现象。
当m＝n时，最末项旳频率为100%，样
本末项值为总体中旳最小值，不符合事实。
5.变化参数，选
定Cv＝0.30， Cs＝0.75，查 表计算出各xP
值。
绘制频率曲线， 该线与经验点据 配合很好，取为 最终采用旳频率 曲线。
配线法得到旳成果仍具有抽样误差，而这种 误差目前还难以精确估算，所以对于工程上 最终采用旳频率曲线及相应旳统计参数，不 但要从水文统计方面分析，而且还要亲密结 合水文现象旳物理成因及地域规律进行综合 分析。
利用Excel绘制P_型频率曲线
良好接地。
四、西潮河闸防雷方案西潮河闸位于射阳县黄沙港镇境内，年雷暴日为46天，属多雷区，加之地处水陆交界处，建筑物及其内部设备极易遭受雷击。
所以本着“安全、可靠、先进、经济、适用”的原则，根据国家有关防雷技术规范的要求，结合工程实际情况，西潮河闸采取了如下防雷方案：1.直接雷防护利用办公楼、启闭机房屋顶的避雷带和四周墙面内的柱钢筋作为引下线，与梁钢筋相互焊接，组成一个完整的接地系统，把进入建筑物的水管、金属管道等金属构件作良好电气连接。
这样，整座水闸、办公楼就形成了一个理想的“法拉第笼”屏蔽网，不但能使雷电流有良好的散流途径，均压分流，而且整座建筑物形成统一的等电位系统，保持均压作用。
为改善接地效果，在办公楼背面及西潮河闸南北端各增设一组接地体，并与现有接地连成一体，以形成共用接地装置，其接地电阻不大于1.0Ω。
2.感应雷与雷电侵入波防护在直击雷防护措施完善合格的前提下，对雷电入侵的可能通道进行必要的防护。
（1）在综合布线时避免将线缆靠近避雷带或引下线布设，因为避雷带及其引下线有雷电流通过时会在周围的邻近导体上产生较强的感应电动势，从而对线路上的设备造成损害。
进出户线路采用埋地缆进入，并用金属导管屏蔽，屏蔽金属管在进入建筑物或机房前重复接地，铠装电缆两端金属外皮接地，其冲击接地电阻均小于30Ω。
（2）电源系统的防护①在总配电房配电柜安装一只电源防雷箱，作为电源一级防护，使大部分雷电能量在该级得到初步释放。
②在办公楼配电箱及闸控室两配电盘分别安装一只电源防雷箱，作为电源二级防护，以降低电源线的雷电波残压，防止其进一步危害设备。
③在各视频监控线路分别安装一只防雷器，各电脑采用电源防雷防雷专用插座，以进一步降低雷电能量，保护设备。
五、结语雷电防护不是简单的避雷设施的安装，而是一项要求高、难度大的系统性工程。
除了要依据系统防雷的科学理论采取相应的防护措施，进行有针对性的防护外，还应委托有资质的专业防雷检测机构定期检测防雷设施，设立防范雷电灾害责任人，对新增建筑和新增设备的防雷系统应进行重新设计和建设。
工程水文水力学思考题和计算题(25题思考问答题,20题计算题)答案
工程水文水力学思考题和计算题一、思考问答1、水文现象是一种自然现象，它具有什么特性，各用什么方法研究？1）成因分析法:据水文变化成因规律，由其影响因素预报、预测水文情势的方法。
如降雨径流预报法、河流洪水演算法等。
2）数理统计法：据水文现象的统计规律，对水文观测资料统计分析，进行水文情势预测、预报的方法。
如设计年径流计算、设计洪水计算、地区经验公式等。
水文计算常常是二种方法综合使用，相辅相成，例如由暴雨资料推求设计洪水，就是先由数理统计法求设计暴雨，再按成因分析法将设计暴雨转化为设计洪水。
此外，当没有水文资料时，可以根据水文现象的变化在地区分布上呈现的一定规律（水文现象在各流域、各地区的分布规律）来研究短缺和无资料地区的水文特征值。
2、何谓水量平衡？试叙闭合流域水量平衡方程在实际工作中的应用和意义。
答：对任一地区、任一时段进入的水量与输出的水量之差，必等于其蓄水量的变化量，这就是水量平衡原理，是水文计算中始终要遵循的一项基本原理。
依此，可得任一地区、任一时段的水量平衡方程。
对一闭合流域：设 P 为某一特定时段的降雨量，E 为该时段内的蒸发量，R 为该时段该流域的径流量，则有：P=R+EC+△U△U为该时段流域内的蓄水量，△U=U1+U2。
对于多年平均情况，△U ＝0，则闭合流域多年平均水量平衡方程变为：⎺P=⎺R+⎺E影响水资源的因素十分复杂，水资源的许多有关问题，难于由有关的成因因素直接计算求解，而运用水量平衡关系，往往可以使问题得到解决。
因此，水量平衡原理在水文分析计算和水资源规划的分析计算中有广泛的应用。
如利用水量平衡式可以用已知的水文要素推求另外的未知要素。
例如：某闭合流域的多年平均降雨量⎺P=1020mm ，多年平均径流深⎺R=420mm，试求多年平均蒸发量⎺E 。
⎺E=⎺P-⎺R＝600mm。
3、何谓年径流？它的表示方法和度量单位是什么？径流深度、径流总量、平均流量、径流模数的概念及相互关系。
2019-6水文统计-PIII型分布-精选文档
0 l
n
使 L 达到极值时，有u i的极大似然估计值和极大似然估计
量如下：
0 0 ˆ ˆ u u x , x , , x ) i i( 1 2 n
0 0 ˆ ˆ u u X , X , , X ) i i( 1 2 n
i 1
p.178 例6
*矩法和极大似然法的异同
4/41
权函数法 —— 加权矩法” p.187
进一步推导，还有
X M 0; M 1 1 Cv H ( ) ; M 2 2 M 0 M 2 3 R M 2 M 1 2
x
这里的H和R都与Cs有关，可通过查表求出，再计算得到x、 Cv和Cs。
7/41
水文适线方法：
目估适线法（配线法、图解法）—— 以经验频率点据为基础，在一定的适线准则下，求解与经验 点据拟合最优的频率曲线参数，是集线型选配和参数估计“合二 为一”的方法(p.179～p.185) 。 基本步骤： 1）将实测水文系列从大到小排列，用经验频率公式(期望公式)
2 2 1 x
1
x
式中，f(t)和f(x)均为P-III型分布的概率密度函数。
6/41
其中，
1 1 1 t x S ( ) x [ t e dt ] x e dx 1 ( ) ( ) 0 0
x
1 1 t 2 1 x S ( ) x [ t e dt ] x e dx 2 ( ) ( ) 0 0
n 1 2 ˆ2 s (x x ) i n1
ˆ C s
3 ( k 1 ) i 1
n
ˆ 3 nC v
3/41
极大似然法 ——“发生概率最大的事件可能发生”
水文统计PIII型分布
§6 P-III型分布参数估计和经验频率曲线（p.173-196）
6.1 估计理论
参数估计 —— X、Cv、Cs 非参数估计 —— 当没有确切的数学分布时
点估计的思想：假设样本X的分布函数为F(x; u )，u 为未知参数，
构造估计量U＝U(X1,X2,…,Xn)，若将一个具体样本(x1,x2, …,xn)带入U， 则可以得到u 的估计值为u＝u(x1,x2, …,xn)。 区间估计的思想：即估计参数所在的区间，或永一个区间区区估计未 知参数。
Qp＝1％＝1246×(1＋3.33×0.6)＝3730m3/s
14/41
三点法 ——
在经验频率曲线上任取3个点(xp1,P)、(xp2,P)、(xp3,P)，代入P-III
i 1 n 2 2 ( K 1 ) i i 1 n
C v
n n 1
n
n 1
12/41
3）点绘经验频率曲线（图4-9）
13/41
4）选配理论频率曲线（表4-4）
5）推求百年一遇的设计洪峰流量
由图4-9，查 p＝1％对应的流量 Qp＝3730m3/s
或用式 x 直接计算，得到 x ( 1 Cv ) p p
0
0 l
n
使 L 达到极值时，有u i的极大似然估计值和极大似然估计
量如下：
0 0 ˆ ˆ u u x , x , , x ) i i( 1 2 n
0 0 ˆ ˆ u u X , X , , X ) i i( 1 2 n
i 1
p.178 例6
*矩法和极大似然法的异同
4/41
权函数法 —— 加权矩法” p.187
n 1 2 ˆ2 s (x x ) i n1
水文频率计算适线法
5.改变参数，选
定Cv＝0.30， Cs＝0.75，查 表计算出各xP
值。
绘制频率曲线， 该线与经验点据 配合较好，取为 最后采用的频率 曲线。
配线法得到的成果仍具有抽样误差，而这种 误差目前还难以精确估算，因此对于工程上 最终采用的频率曲线及相应的统计参数，不 仅要从水文统计方面分析，而且还要密切结 合水文现象的物理成因及地区规律进行综合 分析。
全相关和零相关之间
直线相关
将对应的 xi 与 yi（ i = 1,2,…n ）对应
点绘在方格纸上，如果点据的平均趋势
线为直线，说明变量x与y为线性相关, 满足方程： y ＝ a ＋ bx
相关系数
y 倚 x 回归方程
X 倚y 回归方程
y 倚 x 回归方程
X 倚 y 回归方程
水文频率计算适线法
x
1200
W（X≥xi）＝i / n
1000
800 0
20
ห้องสมุดไป่ตู้
40
60
80
100 W（%）
某地年降雨量经验分布曲线
二、经验频率
如果用P（X≥xi）＝m/n 的经验分布曲
线估计总体分布曲线，存在不合理现象。
当m＝n时，最末项的频率为100%，样
本末项值为总体中的最小值，不符合事实。
水文上用期望值公式估计频率
作业4
某水库坝址处的年平均流量资料如表1 所列，现要求用适线法（Cs＝2Cv）推 求设计标准P＝90%的设计年径流量。
表1 某水库坝址处年平均流量表 单位：m3/s
年份 1959 1960 1961 1962 1963 1964 1965 1966 1967 1968 流量 11.9 7.78 10.0 9.64 14.4 4.73 7.83 10.4 10.2 10.9 年份 1969 1970 1971 1972 1973 1974 1975 1976 1977 1978 流量 12.6 10.3 15.1 7.24 11.3 11.7 8.42 16.9 6.82 5.74 年份 1979 1980 1981 1982 1983 1984 1985 1986 1987 1988 流量 15.6 10.7 5.15 7.27 13.1 7.72 6.42 12.9 13.2 9.49
巧用excel实现P-Ⅲ型频率曲线绘制
巧用excel实现P-Ⅲ型频率曲线绘制摘要：本文提出了用Excel软件解决水文计算中P－Ⅲ型频率曲线的特殊函数计算方法。
实践证明，该方法计算准确、快捷和规范，绘图标准、精美，有较高的实用价值。
关键词: Excel；频率计算；图形绘制Abstract: This paper puts forward using Excel software to solve the calculation of hydrology P - Ⅲtype frequency curve calculation method of special function. Proved by practice, the method to calculate accurately, fast and standardized, drawing standards, exquisite, have high practical value.Key words: Excel; frequency calculation; graphics rendering1概述在工程设计的水文计算以及实际运用的水文复核工作中，经常推求不同频率洪水或各种保证率的水位及流量等，这就要对一些样本系列进行排频并绘制PⅢ频率曲线。
一般都采用手工计算的方法绘制经验频率曲线，然后进行适线，推求理论频率曲线，计算查表工作很繁琐。
Excel是全世界最广泛使用的办公软件之一。
Excel本身功能强大，包括打印，文件处理，格式化和文本编辑。
常用的Excel内置大量函数，界面熟悉，可连接到多种数据库。
经过一些简单设置，完全可以胜任此项工作，既简便、又经济。
2 机率格纸的绘制普通坐标纸上绘制的频率曲线存在着两端曲率大，曲线不便外延的问题，1913年海森提出使用正态机率格纸，即以x为纵坐标，以为横坐标，使正态曲线在这种图纸上呈直线形式。
尽管水文变量大多不呈正态分布，绘在海森机率格纸上的频率曲线仍是曲线，但曲率已大为减缓，所以生产上一直使用已印制好的海森机率格纸进行频率计算。
水文学与水文地质学3经验曲线-理论曲线
曲线绘制步骤
曲线绘制步骤：
① 据实测水文系列资料，不论年序，数值从大到小 排列，统计m，确定n；
② 计算（查表） P = m/(n+1)（%）； ③ 以实测值x为纵坐标，P(%)为横坐标，点绘
(P1, x1 ), (P2 , x2 ), (Pn , xn )，连线；
④ 若资料充分，可据设计频率标准[P]，在该
2
2390
9.0909091 13
1050 59.0909091
3
1860
13.6363636 14
1050 63.6363636
4
1740
18.1818182 15
980
68.1818182
5
1540
22.7272727 16
980
72.7272727
6
1520
27.2727273 17
883
77.2727273
1952
2750
1963
1953
762
1964
1954
2390
1965
1955
1210
洪峰流量Q（m3/s）
1270 1200 1740 883 1260 408 1050 1520 483 794
表整理
m
序号
1 2 3 4 5 6 7 8 9 10 11
洪峰流量
2750 2390 1860 1740 1540 1520 1270 1260 1210 1200 1140
理论频率曲线
纵坐标：
xp x(1 Cv p )
横坐标为理论点P(%)或
P( X
xp)
( )