CN106446547B - 一种基于广义第二类贝塔分布的洪水频率分析方法 - Google Patents

一种基于广义第二类贝塔分布的洪水频率分析方法 Download PDF

Info

Publication number
CN106446547B
CN106446547B CN201610856618.7A CN201610856618A CN106446547B CN 106446547 B CN106446547 B CN 106446547B CN 201610856618 A CN201610856618 A CN 201610856618A CN 106446547 B CN106446547 B CN 106446547B
Authority
CN
China
Prior art keywords
flood
distribution
parameter
frequency analysis
function
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201610856618.7A
Other languages
English (en)
Other versions
CN106446547A (zh
Inventor
陈璐
何典灿
黄康笛
熊丰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN201610856618.7A priority Critical patent/CN106446547B/zh
Publication of CN106446547A publication Critical patent/CN106446547A/zh
Application granted granted Critical
Publication of CN106446547B publication Critical patent/CN106446547B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00—Subject matter not provided for in other main groups of this subclass
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Complex Calculations (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于广义第二类贝塔分布的洪水频率分析方法，属于水文学中的洪水预报领域。该发明方法率先将广义第二类贝塔分布引入洪水频率分析中，建立洪水概率密度预测模型，并依据分布函数特点采用最大熵原理实现洪水概率密度预测模型的参数推导，将河流历年洪峰序列代入参数推导公式中求出参数，再将参数代入洪水概率密度预测模型，最后利用完整的洪水概率密度预测模型预测大规模洪水的概率预测。本发明的洪水概率分析方法非常适合水文频率分析，其拟合效果基本优于水文中的其它传统分布，本发明为水文频率分析提供了一种更加有效的途径。

Description

一种基于广义第二类贝塔分布的洪水频率分析方法
技术领域
本发明属于洪水预报领域，更具体地，涉及一种基于广义第二类贝塔分布(Generalized beta distribution of the second kind，GB2)的洪水频率分析方法。
背景技术
洪水一般用洪水发生的概率密度定义洪水大小，如20年一遇、50年一遇和百年一遇等。洪水频率分析的目就是通过频率曲线的外延，推求T年一遇洪水的设计洪水值。如计算的设计洪水值较大，则规模过大，会增加投资，造成浪费；如计算的设计洪水值较低，则规模过小，又可能在不利水文条件下导致工程失事造成损失。
因此，开展高精度的洪水频率分析是水利水电工程设计与规划的首要问题，选择合适的频率分布线型与参数估计方法是其重要内容。
目前多数研究仅选取单一的频率分布线型，如指数分布(Exponential，EXP)、Weibull分布、Gamma分布、Gumbel分布、广义极值分布(GEV)、皮尔逊III型分布(P-III)、对数皮尔逊III型分布(LP-III)及对数正态分布(LN)等分布进行洪水频率分析，其设计洪水结果具有较大的不确定性，设计洪水值的高估与低估会导致投资过多或安全风险增大等后果。
发明内容
针对现有技术的以上缺陷或改进需求，本发明提供了一种基于广义第二类贝塔分布(GB2)的洪水频率分析方法，其目的在于将广义第二类贝塔分布引入洪水频率分析中，并依据分布函数特点采用最大熵原理实现分布函数的参数估计，在此基础上开展基于广义第二类贝塔分布的设计洪水计算，由此解决现有分析技术的洪水频率分析精度不高的技术问题。
为实现上述目的，按照本发明的一个方面，提供了一种基于广义第二类贝塔分布的洪水频率分析方法，该方法包括以下步骤：
(1)采用年最大取样方法,采集年最大流量样本序列；
(2)应用广义第二类贝塔分布建立洪水概率密度模型：
其中，变量x表示年最大洪水流量，f(x)表示流量为x的概率密度；B()为贝塔函数；r1、r2、r3是形状参数，且r1>0,r2>0,r3＞0，r3决定全局形状，r1控制左尾形状，r2控制右尾形状，参数r1和r2共同决定分布的偏态特性；b为位置参数，b＞0；
(3)推求基于最大熵原理的广义第二类贝塔分布的概率密度函数参数值，包括以下子步骤：
(31)引入最大熵原理，最大熵原理指出，某一随机变量x的概率密度函数f(x)可通过最大化熵值得到：
其中，gi(x)为x的函数；Ci为gi(x)的期望；
依据拉格朗日乘子法，f(x)可表示为：
f(x)＝exp(-λ0-λ1g1(x)-λ2g2(x)…λmgm(x)) (3)
其中：m为约束个数；λi,i＝0,1,2,...,m为拉格朗日乘子；
(32)基于最大熵原理，广义第二类贝塔分布约束条件表达式为：
其中，q为构造的约束条件参数；E为期望；
则广义第二类贝塔分布的概率密度函数可构造为：
其中，λ2'为构造的拉格朗日乘子；
(33)推导拉格朗日乘子与广义第二类贝塔分布函数的约束条件之间关系，将式(5)带入到(4a)中可得：
令p＝b-q、和t＝pxq，则上式可化简为：
令则因为y(0)＝0且y(∞)＝1，y∈[0,1]，则：
即有下式成立：
同时，计算λ0的另一种方法为：
令随后对式(9)、(10)中的λ1和λ2分别求导，可得：
式中：为digamma函数；
则拉格朗日乘子与广义第二类贝塔分布函数约束条件之间满足：
对式(11)求二阶导可得下式：
(34)推导拉格朗日乘子与广义第二类贝塔分布函数参数之间关系，将式(8)带入到式(6)中，并与式(1)作比较，则下式成立：
(35)推求广义第二类贝塔分布函数参数与约束条件之间关系，将式(16)代入式(13)和式(14)中，可得广义第二类贝塔分布函数参数与约束条件之间满足下式：
式中，var为方差；将步骤(1)中年最大流量样本序列,即x序列代入公式(17)，求出参数r1、r2、r3和b；
(4)将参数r1、r2、r3和b代入公式(1)中，即可得到洪水概率密度函数，可利用该函数得到洪水为某一量级x时的概率值f(x)；还可利用洪水概率密度函数，推求概率f(x)为T年一遇的设计洪水值x。
总体而言，通过本发明所构思的以上技术方案与现有技术相比，具有以下技术特征及有益效果：
四参数的广义第二类贝塔分布包含了指数分布、威布尔(Weibull)分布及伽马(Gamma)分布，具有足够的灵活性模拟复杂多变的数据集，因此非常适于水文频率分析。本发明通过引入广义第二类beta分布，并利用最大熵原理推求估计了分布函数参数，实现了精度较高的洪水频率分析。其拟合效果优于水文中的其它常用分布；最大熵原理能够有效地估计分布函数的参数，得到精度高的洪水频率分布函数。
附图说明
图1为本发明方法流程图；
图2为本发明方法拟合的洪水序列边缘分布图及概率密度曲线效果图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
如图1所示，本发明包括以下流程：
(1)采用年最大取样方法,采集年最大流量样本序列；
(2)应用广义第二类贝塔分布建立洪水概率密度模型：
其中，变量x表示年最大洪水值，f(x)表示流量为x的概率密度；B()为贝塔函数；r1、r2、r3是形状参数，且r1>0,r2>0,r3＞0，r3决定全局形状，r1控制左尾形状，r2控制右尾形状，参数r1和r2共同决定分布的偏态特性；b为位置参数，b＞0；
(3)基于最大熵原理，推导出广义第二类贝塔分布函数的参数计算公式：
将步骤(1)中采集年最大流量样本序列数据作为序列x，代入以上公式，求出参数r1、r2、r3和b；
(4)将参数r1、r2、r3和b代入公式(1)中，得到洪水概率密度函数，可利用该函数输入需预测洪水流量x，计算得到洪水流量达到x的概率密度f(x)；还可利用洪水概率密度函数，输入洪水发生的概率密度f(x)，计算发生该概率下洪水流量x。
实施例：现利用某流域典型水文站点A、B和C的年最大日流量数据，检验广义第二类贝塔分布的模型效果。基于最大熵原理的方法，由A、B和C三个水文站点的年最大日流量序列求出模型的参数，表1给出了A、B和C三个水文站点洪水概率密度函数的参数，图2给出了基于GB2分布的3站年最大洪水序列边缘分布和概率密度曲线的拟合结果，结果显示GB2分布拟合效果良好，可用于设计洪水的计算。
表1
采用GB2分布、正态分布、指数分布(EXP)、Gamma分布、Gumbel分布、广义正态分布(Generalized normal，GN)、P-III分布、广义Pareto分布(Generalized Pareto，GP)及Weibull分布等拟合典型水文站C的年最大日流量序列，并基于K-S(Kolmogorov-Smirnov，K-S)检验法、均方根误差以及AIC准则(Akaike information criterion，AIC)对各分部的拟合结果进行比较分析，确定拟合最优的分布函数，进行洪水频率分析。
取K-S检验显著性水平为α＝0.05，P值大于0.05时通过检验，RMSE和AIC值越小说明拟合效果越好，表2给出了各分布K-S检验统计量P值、RMSE和AIC值。结果表明，所有分布均通过了K-S检验。比较各分部的RMSE和AIC值可知，GB2分布的拟合效果最优。
表2
表3给出了基于上述分布函数计算的水文站点C的设计洪水值。结果表明，当重现期较大时，各分布计算结果差异明显，考虑到工程设计中，低估的设计值将大大增加大坝及下游的洪水风险，由表中数据比较分析可知GB2的频率分析结果更优。
表3
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (2)

1.一种基于广义第二类贝塔分布的洪水频率分析方法，其特征在于，该方法包括以下步骤：
(1)采用年最大取样方法,采集年最大流量样本序列；
(2)应用广义第二类贝塔分布建立洪水概率密度模型：
其中，变量x表示年最大洪水流量，f(x)表示流量为x的概率密度；B()为贝塔函数；r1、r2、r3是形状参数，且r1>0,r2>0,r3＞0，r3决定全局形状，r1控制左尾形状，r2控制右尾形状，参数r1和r2共同决定分布的偏态特性；b为位置参数，b＞0；
(3)基于最大熵原理，推导出广义第二类贝塔分布函数的参数计算公式：
其中，E为期望；表示digamma函数；var()为方差；将步骤(1)中采集河流历年最大洪峰流量序列数据即x序列代入以上公式，求出参数r1、r2、r3和b；
(4)将参数r1、r2、r3和b代入公式(1)中，得到洪水概率密度函数，利用该函数输入洪水发生的概率密度f(x)，推求发生该概率下洪水流量x。
2.根据权利要求1所述的一种基于广义第二类贝塔分布的洪水频率分析方法，其特征在于，所述步骤(4)中还可利用洪水概率密度函数，输入需预测洪水流量x，计算得到洪水流量达到x的概率密度f(x)。
CN201610856618.7A 2016-09-27 2016-09-27 一种基于广义第二类贝塔分布的洪水频率分析方法 Active CN106446547B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610856618.7A CN106446547B (zh) 2016-09-27 2016-09-27 一种基于广义第二类贝塔分布的洪水频率分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610856618.7A CN106446547B (zh) 2016-09-27 2016-09-27 一种基于广义第二类贝塔分布的洪水频率分析方法

Publications (2)

Publication Number Publication Date
CN106446547A CN106446547A (zh) 2017-02-22
CN106446547B true CN106446547B (zh) 2018-12-28

Family

ID=58170618

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610856618.7A Active CN106446547B (zh) 2016-09-27 2016-09-27 一种基于广义第二类贝塔分布的洪水频率分析方法

Country Status (1)

Country Link
CN (1) CN106446547B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110069750A (zh) * 2019-04-12 2019-07-30 电子科技大学 基于四参数威布尔分布的非精确分析方法
CN110598181B (zh) * 2019-08-30 2021-07-02 华中科技大学 一种基于最大熵的极端水文事件风险分析方法和系统
CN112149893A (zh) * 2020-09-17 2020-12-29 中国水利水电科学研究院 一种设计洪水值预测方法和装置
CN112214734A (zh) * 2020-10-10 2021-01-12 东北大学 一种基于统计物理学与人工智能的电力负荷预测方法
CN116011687B (zh) * 2023-03-30 2023-08-11 山东锋士信息技术有限公司 一种基于Copula函数的洪水预报方法、系统及介质

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4362042B2 (ja) * 2003-05-15 2009-11-11 東京電力株式会社 洪水ハザードマップ生成方法及びシステム
KR100707101B1 (ko) * 2006-07-19 2007-08-10 이화여자대학교 산학협력단 센서 및 이동 네트워크에서 유사 엔트로피를 이용하여플러딩 공격을 탐지하는 방법
CN104615845B (zh) * 2014-12-19 2017-12-26 兰州大学 一种超临界洪水风险分析方法
CN105740969B (zh) * 2016-01-21 2019-08-06 水利部交通运输部国家能源局南京水利科学研究院 一种基于数据驱动的小流域洪水实时预报方法

Also Published As

Publication number Publication date
CN106446547A (zh) 2017-02-22

Similar Documents

Publication Publication Date Title
CN106446547B (zh) 一种基于广义第二类贝塔分布的洪水频率分析方法
CN111898831B (zh) 一种实时洪水概率预报实用化方法
CN110276104B (zh) 一种集合气候模式下的分期设计洪水推求方法
CN108830419B (zh) 一种基于ecc后处理的梯级水库群入库流量联合预报方法
Liu et al. A spatial panel statistical analysis on cultivated land conversion and Chinese economic growth
CN101697172B (zh) 一种水文频率线型参数估计方法
CN110598181B (zh) 一种基于最大熵的极端水文事件风险分析方法和系统
CN107886160B (zh) 一种bp神经网络区间需水预测方法
CN109635372B (zh) 基于改进先验概率的贝叶斯模型的设计洪水方法
CN109902267B (zh) 一种受下游湖泊顶托影响的河道安全泄量计算方法
CN111709147A (zh) 一种基于水文失事机理的设计洪水地区组成计算方法
CN106557614B (zh) 一种基于Halphen IB分布的洪水频率分析方法
CN106529157B (zh) 一种基于Halphen B分布的洪水频率分析方法和系统
CN114969655A (zh) 流域输沙量的模拟估算方法
Ma et al. Bayesian statistic forecasting model for middle-term and long-term runoff of a hydropower station
CN112116200A (zh) 基于动态比例代换和层次贝叶斯的缺灾情资料城市洪灾损失函数的构建方法
CN117077420A (zh) 基于Copula函数的荒漠河岸林生态保护阈值确定方法
CN106547981B (zh) 一种基于Halphen A分布的洪水频率分析方法和系统
CN109614642A (zh) 一种水文频率分析的模型选择方法
CN106599566B (zh) 一种基于广义伽马分布的洪水频率分析方法和系统
CN115130736B (zh) 基于贝叶斯层次模型的水文生态响应关系、生态流量确定方法及系统
CN109390946B (zh) 一种基于多参数规划理论的最优概率潮流快速计算方法
Ali et al. Statistical downscaling of precipitation and temperature using gene expression programming
Xing et al. Nash model parameter uncertainty analysis by AM-MCMC based on BFS and probabilistic flood forecasting
CN117390125B (zh) 一种产流模式智能适配的流域水文预报方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant