根据某水文现象的统计特性,利用现有水文资料,分析水文变量设计值与出现频率(或重现期)之间的定量关系的工作过程称为水文频率计算。

自然界的现象按发生情况可分成:必然事件,即在一定条件下必然会发生的事情,如降雨以后就要涨水是必然发生的;不可能事件,即在各条件实现之下永远不会发生的事情,如只在重力作用下的水由低处向高处流是不可能的;随机事件(也称偶然事件),即在一定条件下可能发生也可能不发生的事件,如每条河流每年出现一个流量的年最大值是必然的,但这个最大值可能是这个值也可能是那个值,它在数量上的出现是一种随机事件。频率计算中是以 1来表示必然事件出现的可能性(即百分之百出现),以0表示不可能事件出现的可能性,随机事件出现的可能性介于0与1之间。

水文要素。如降雨、流量等在量的出现方面都有随机性的特点,水文变量如年雨量、年最大洪峰流量、枯季最小流量等都属于随机事件,均可用频率分析方法来分析计算。 水文频率分析主要包括: 利用现有水文资料组成样本系列,选择合适的频率曲线线型和估计它的统计参数,根据所绘制的频率曲线推求相应于各种频率(或重现期)的水文设计值。

样本系列。无限个成因相同、相互独立的同类水文变量的集合称为该水文变量的总体。这个总体是未知的,现有水文资料只是过去发生过的和今后可能发生的整个总体中的一个样本。把现有水文资料的水文变量按大小次序排列组成一个系列,称为样本系列,其中所含水文变量的项数(系列长度)叫做样本容量。系列愈长,样本容量愈大。水文频率分析就是通过样本系列的统计特征来估计其总体的统计特征,如各种统计参数、某水文变量的频率等。因此,样本系列是水文频率分析的基础。用样本系列去推估容量很大或无限的总体的情况,会产生因抽样而引起的误差,这就是抽样误差。水文统计分析中所估计出的各种数值(如频率、分析中的各个参数、相关系数等)都有抽样误差。样本的容量越大误差越小,否则误差越大。抽样误差分析方法有两种:①解析法。用统计原理推求出抽样误差的公式,按公式求得抽样误差值。例如,均值的均方(抽样)误差值为,其中Cv为所研究变量系列的离差系数,n为系列的长度或样本容量。②统计试验法。即生成很长的资料系列,来研究样本容量一定时统计分析中各种数值的抽样误差。

经验频率。样本系列中某水文变量x大于或等于一定数值xm(即x≥xm)的可能性大小即为频率,一般用符号pm{x≥xm}来表示,其值在0与1之间。例如,某河段年最大洪峰流量系列中,出现流量Q≥1000米3/秒的可能性为百分之一,则称Q≥1000米3/秒的频率等于1%。设系列共有n项,其中第m 项xm的频率Pm常用下列公式来计算:
水文频率分析中,称上式为经验频率公式,而Pm亦称为系列中第m 项的经验频率。经验频率在绘制频率曲线的适线法中应用。

重现期。指某水文变量的取值 (x≥xm)在很长时期内平均多少年出现一次。重现期(T)与频率(P)的关系对下列两种情况有不同的表示方法:①当研究防洪治涝的暴雨、洪水时,采用设计频率P<50%,则T=1/P(年)。例如,当P=1%时,得T=100年,称为百年一遇的暴雨或洪水。②当考虑兴利或枯水问题时,采用设计频率 P>50%,则T=1/(1-P)(年)。例如,当灌溉的设计频率P=80%时,得T=5年,称为5年一遇的枯水。象暴雨或洪水那样的水文现象并无固定的周期,所谓百年一遇是指大于或等于这样的洪水在很长时期内平均每百年出现一次,而不能理解为恰好每隔百年出现一次。对于具体的 100年来说,超过这种洪水可能不止一次,也可能一次都不出现,而只是说明长时期内平均每年出现的可能性为1%。

统计参数。资料系列的数量水平和变化幅度等情况的综合特征值称为统计参数。绘制频率曲线,除了需掌握系列各项的经验频率之外,还须了解系列的统计参数。水文频率分析中,常用三个统计参数,即均值(算术平均值的简称)塣(zheng.4)、离差系数Cv (也称变差系数)和偏差系数Cs。均值是集中表示系列数量级大小或水平高低的指标,例如对降雨系列,均值大的表示雨量充沛,反之表示雨量稀少。离差系数表示系列中各项值对其均值的相对离散程度的指标,它是系列均方差与均值之比。如果离差系数Cv较大,即系列的离散程度较大,亦即系列中各项的值对均值离散较大,如果Cv较小,则系列的离散程度较小,亦即系列各项的值同均值相差较小。偏差系数是表示系列中各项的值偏于均值左右的情况的相对指标。如果大于均值的各项值占优势称为正偏(Cs>0);若小于均值的各项值占优势称为负偏(Cs<0);当大于均值和小于均值的各项值都不偏时称为对称 (Cs=0)。

频率曲线。把水文变量和频率表达成一定的数学关系式并将它画成图形,即为频率曲线。其线型常用的有:Γ分布或皮尔孙Ⅲ型分布曲线,极值Ⅰ型分布或贡贝尔分布曲线,对数正态分布曲线,对数Γ分布或对数皮尔孙Ⅲ型分布曲线等。频率曲线常画在概率格纸上。这种格纸的纵坐标为均匀分格或对数分格,表示水文变量;横坐标按某种概率分布(一般取正态分布)分格,表示频率。将水文变量同其相应的经验频率关系(称为经验频率点)点绘在概率格纸上,用一定的配线准则,可以配出一条频率曲线。从20世纪60年代以后,中国一般采用皮尔孙Ⅲ型曲线。有时(如在初估时),也有用目估方法通过经验频率点的点群中心徒手绘制频率曲线的

频率曲线的绘制方法。每种频率曲线均含有一定个数的统计参数,一般有三个,即均值、离差系数和偏差系数。在频率曲线线型选定之后,就要估计这些参数。有了这些参数,就可以绘出频率曲线。统计参数的估计方法常用下列几种:①适线法。在概率格纸上用频率曲线去配合样本系列的经验频率点据,取用配合较佳时的统计曲线。这种方法常用的有目估适线法和优化适线法。目估适线法是通过工作者的目测,以他认为曲线与点据配合较佳时的曲线为准;这种方法具有一定的任意性,不同工作者会得到不同的结果,但它能照顾精度较高或占重要位置的点据。优化适线法要用一定形式的目标函数,使其最小而得,最小二乘法和离差绝对值之和最小法属于此类;这种方法可以避免适线的任意性,在统计试验法中应用较好,而在实测资料的分析中,难以照顾精度较高或占重要地位的点据。②矩法。是用头几阶矩来估计频率曲线统计参数的方法,矩的阶数同统计参数的个数相同。③极大似然法。是使估出的统计参数代入频率曲线的密度函数中,得到系列中各水文值对应的频率密度之乘积为最大。中国主要采用适线法。

成果合理性检查。现有的水文观测资料一般较短,至多百年左右。在推求千年一遇或万年一遇水文设计值(如千年一遇或万年一遇的洪峰流量)时,必须把频率曲线外延,外延愈远,估计所得的水文设计值的误差愈大。因此,水文频率分析时,要求尽可能地调查历史上发生过的大洪水,参证审查后加入频率分析。同时,必须对频率分析成果在时间上(单站各长短时段)和空间上(情况相似的地区上)作合理性分析。例如,相同频率时短时段的水文值不能大于长时段的水文值,相邻站同类水文系列的统计参数不能相差太大等。

水文频率分析http://baike.baidu.com/view/958646.htm

水文频率计算软件http://www.civilcn.com/shuili/slrj/shuiwen/129004932064273.html,

matlab 水文频率曲线,【求助】如何计算水文频率,外行求教 - 地学 - 小木虫 - 学术 科研 互动社区...相关推荐

  1. arcgis 同名图层合并_【求助】ArcGIS中怎样合并图层 - 地学 - 小木虫 - 学术 科研 互动社区...

    这个问题 最简单的方式是用editor中的切割多边形功能  可以把一个矢量文件按照你定义的线分割  可直接在属性表中算出面积  要是想直接用一个线文件去切割面文件的话 操作有点复杂  下面是用edit ...

  2. matlab 吸波材料,如何计算多层吸波材料的反射率 - 材料 - 小木虫 - 学术 科研 互动社区...

    各位大神,小女子最近在做多层吸波材料,可以测得每层材料的电磁参数,然后请同学编了一个计算多层吸波材料反射率的MATLAB程序,但行出来总觉得好像有点问题,有没有哪位可以帮我看一下程序对不对呢? PS: ...

  3. matlab 摄氏度符号怎么打,【求助】弱弱的问下,Matlab如何在公式中输入希腊字母? - 程序语言 - 小木虫 - 学术 科研 互动社区...

    很多时候都要在matlab画图的时候添加一些公式符号之类的,有一些特殊的字符并不能直接从键盘上输入,比如希腊字母等等.但是有想用,因为这样使图看起来漂亮而且容易理解. 例如:我想输入摄氏度的符号,怎么 ...

  4. matlab 工业相机 曝光时间_Matlab CCD工业相机采集图像问题 - 信息科学 - 小木虫 - 学术 科研 互动社区...

    我通过imaqhwinfo函数来查询我的电脑里安装的图像采集设备(我的笔记本电脑自带一个摄像头,同时我外接了一个CCD相机+显微镜头),查询结果如下: 第一部分: info = InstalledAd ...

  5. matlab比Excel好在哪,matlab计算结果和excel完全不一样 - 仿真模拟 - 小木虫 - 学术 科研 互动社区...

    请教各位: 我有一段程序通过matlab计算出的结果和用Excel中VB计算的结果完全不一样,跪求指导.....(已经确定是程序中这段的问题,实在找不出2者的区别) 两种语言中中Ndata=n,  i ...

  6. matlab实现kuan滤波,用Matlab设计一个0.5-20HZ的带通滤波器! - 程序语言 - 小木虫 - 学术 科研 互动社区...

    这是我在网上看到的,希望对你有所帮助. ----------------------------------------------------------------------------- fu ...

  7. 滑动轴承matlab程序,小弟求滑动轴承非线性油膜力matlab程序,最好是差分法的,谢谢 - 机械 - 小木虫 - 学术 科研 互动社区...

    % This function is used to calculate the pressure distribution of journal bearing without taking the ...

  8. matlab decomposition filters,MATLAB小波去噪求助(附算法和显示图片)!不知自己哪个地方出了问题,求指点! - 信息科学 - 小木虫 - 学术 科研 互动社区...

    star2883 不是非常明白楼主的意思.我的理解是: 问题1:为什么中间两幅图,显示的都是空白? 答:你是指Figure2和Figure3的显示中,全是白色吗?这个是因为你用了colormap(ma ...

  9. qpsk的matlab平方根升余,【求助】求助:【QPSK功率谱密度的matlab的代码】 - 仿真模拟 - 小木虫 - 学术 科研 互动社区...

    CODE: clear; %随机产生长度大于1000的'0'.'1'信号序列,对其进行QPSK调制 %定义待仿真序列的维数 N global N N=2000; %定义产生'1'的概率为 p glob ...

最新文章

  1. linux ll以编码格式_【系列】每天一个 Linux 命令(5):rm 命令
  2. shell获取ip的值
  3. vyos User Guide
  4. 深度学习笔记(46) 深度卷积网络学习
  5. 高薪职业,英年早秃?关于程序员,多的是你不知道的事……
  6. ElasticSearch head 安装
  7. widows升级nodejs版本
  8. Latex中将图例(图1:)改为(图1-1)
  9. Mysql数据库基础知识总结,结构分明,内容详细
  10. 拼多多校招内推编程题3 六一儿童节
  11. 浮点数 C语言 IEEE754
  12. 远程桌面连接不上是什么原因?怎么解决
  13. gts250 linux驱动下载,下载:NVIDIA显卡Linux驱动190.32测试版
  14. Pt100转RS-485,热电阻温度Modbus数据采集模块 WJ125
  15. 开源正在蚕食整个世界?
  16. 使用RMAN实现异机备份恢复(WIN平台)
  17. 判断两个链表是否交叉,并求出交叉点
  18. 【Node.js】解决中文乱码问题
  19. 基于脑机接口的光感知视觉机制模型
  20. html中插入音乐播放器,网页中插入音乐播放器WMA 方法有两种:object和embed

热门文章

  1. 历时 4 年,阿里云推出金融核心系统转型实践书
  2. Docker 终失 C 位,开源少年缘何黑化
  3. c语言怎么写贪吃蛇,刚学C语言,想写一个贪吃蛇的代码
  4. iostat命令详解_对iostat输出结果的理解
  5. html伸缩布局,CSS3 伸缩布局(一)
  6. VS Code 报错Vetur can‘t find ‘tsconfig.json‘ or ‘jsconfig.json‘的解决方法
  7. Sublime Text 3 快捷键总结(简洁版本)
  8. 解决window的bat脚本执行出现中文乱码的问题
  9. IDEA解决sun.misc.BASE64Encoder找不到jar包的解决方法
  10. linux tail命令_操作手册