CN106407543B - 一种设计条件下双峰型洪水过程线形状推求方法 - Google Patents
一种设计条件下双峰型洪水过程线形状推求方法 Download PDFInfo
- Publication number
- CN106407543B CN106407543B CN201610813782.XA CN201610813782A CN106407543B CN 106407543 B CN106407543 B CN 106407543B CN 201610813782 A CN201610813782 A CN 201610813782A CN 106407543 B CN106407543 B CN 106407543B
- Authority
- CN
- China
- Prior art keywords
- flood
- design condition
- parameter
- peb process
- under
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- E—FIXED CONSTRUCTIONS
- E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
- E02B—HYDRAULIC ENGINEERING
- E02B1/00—Equipment or apparatus for, or methods of, general hydraulic engineering, e.g. protection of constructions against ice-strains
-
- 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
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Mechanical Engineering (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种设计条件下双峰型洪水过程线形状推求方法,通过统计水库历年最大洪水过程线资料并计算得到相应的时段洪量,归一化得到对应的单位洪水过程线,在采用混合概率密度函数对归一化的单位洪水过程线进行优化拟合得到洪水形状参数系列的基础上,对洪水形状参数系列进行频率分析,计算获得设计条件下的洪水形状参数值,进而推求设计条件下双峰型洪水过程线形状。本发明提供了基于混合概率密度函数的设计条件下双峰型洪水过程线形状的推求方法,统计理论基础较强,更加符合设计条件下的暴雨洪水产汇流特性,更好地满足大中型水库工程设计中对双峰型洪水过程的实际需求。
Description
技术领域
本发明属于水库防洪规划设计领域,特别涉及一种设计条件下双峰型洪水过程线形状推求方法。
背景技术
我国大中型水库的防洪安全设计一般以设计洪水过程线为依据。现行拟定设计洪水过程线方法的程序是先通过洪水频率分析得到某一指定重现期的设计洪峰流量和设计时段洪量,然后从实测资料中选择一种典型洪水过程线形状,采用同倍比或同频率方法放大推求设计洪水过程线[1]。这就使得典型洪水过程线形状的选择结果将对水利工程建设规模及运行管理产生重要影响。
在工程实践中,通常基于“峰高量大、主峰靠后、对水库防洪偏于不利”等定性原则选取典型洪水过程线形状,理论依据不足,且选择的典型往往因人而异,任意性较大[2-3]。事实上,洪水过程线的形状千变万化,在稀遇频率(如千年一遇)的设计条件下由于设计暴雨往往超过实测的常遇暴雨很多,其时空分布和产汇流规律与实测的典型暴雨不同,从而导致产生的设计洪水形状同典型洪水之间可能存在较大的差异[4]。因此,推求符合设计条件下产汇流特性的洪水过程线形状并作为工程设计依据更加科学合理。
考虑到随机变量概率密度函数的形状与洪水过程线具有很大的相似性,有水文学者提出采用概率密度函数来表示洪水过程线的形状[5]。该法采用概率密度函数对归一化的单位洪水过程线进行拟合,并可以通过调整概率分布的参数组来控制概率密度函数的形状,进而模拟多种可能的洪水过程形状类型。此外,概率密度函数法的另一个优点是将洪水过程线形状转化为易于处理的形状参数随机变量来表征,通过对形状参数变量进行频率分析就可以计算得到设计条件下的形状参数值,在此基础上推求设计洪水过程线形状,统计理论基础较强[6]。
然而,现有技术只能提供设计条件下单峰型洪水过程线形状,存在较大的局限性,无法满足大中型水库工程设计中对双峰型洪水过程的实际需求。目前,没有文献研究了设计条件下双峰型洪水过程线形状的推求方法。
本发明涉及的参考文献如下:
[1]张波,张敏,范瑞琪,等.推求设计洪水过程线的方法,CN104182650A.
[2]李敏,薛志春,张艳平.水库典型洪水选择的可变模糊模式识别方法[J].水利科技与经济,2013,19(6):3-5.
[3]钟平安,张阳,张宇,等.一种保持典型洪水形态的设计洪水过程解析推求方法及系统,CN104598676A.
[4]葛慧,黄振平,王银堂,等.基于模糊识别理论的典型洪水选择[J].水电能源科学,2011,29(3):54-56.
[5]Yue S,Ouarda T B M J,Bobée B,et al.Approach for describingstatistical properties of flood hydrograph[J].Journal of HydrologicEngineering,2002,7(2):147-153.
[6]肖义,郭生练,方彬,等.设计洪水过程线方法研究进展与评价[J].水力发电,2006,32(7):61-63.
发明内容
针对现有技术存在的不足,本发明提供了一种设计条件下双峰型洪水过程线形状推求方法。
为解决上述技术问题,本发明采用如下的技术方案:
一种设计条件下双峰型洪水过程线形状推求方法,其特征在于包括以下步骤:
步骤1,收集水库历年最大洪水过程线资料;
步骤2,根据步骤1中控制时段水库流量数据,选定历年最大洪水过程线,进而计算出历年最大洪水过程的时段洪量,对历年最大洪水过程线进行归一化得到对应的单位洪水过程线;
步骤3,根据步骤2中的单位洪水过程线,选取混合Beta分布作为混合概率密度函数线型,并估计混合概率密度函数的参数,获得洪水形状参数系列;
步骤4,对步骤3中的洪水形状参数系列进行频率分析,选取Gumbel分布作为洪水形状参数的概率分布函数线型,并估计概率分布函数的参数,计算设计条件下的洪水形状参数值;
步骤5,将步骤4中设计条件下的洪水形状参数值代入步骤3中的混合概率密度函数线型,推求设计条件下双峰型洪水过程线形状。
作为优选,所述步骤1中收集水库历年最大洪水过程线资料的时间尺度为△t,所述△t范围为0-24小时。
作为优选,所述步骤3中,估计混合概率密度函数的参数采用非线性优化算法、模拟退火算法、极大似然法、以及极大似然的EM算法和最大熵准则法的任意一种或几种组合。
作为优选,所述步骤4中,估计概率分布函数的参数采用的方法为线性矩法、权函数法或概率权重矩法。
本发明在采用混合概率密度函数对归一化的单位洪水过程线进行优化拟合得到洪水形状参数系列的基础上,对洪水形状参数系列进行频率分析,计算获得设计条件下的洪水形状参数值,进而推求设计条件下双峰型洪水过程线形状。
与现有技术相比,本发明的有益效果在于:提供了基于混合概率密度函数的设计条件下双峰型洪水过程线形状的推求方法,统计理论基础较强,更加符合设计条件下的暴雨洪水产汇流特性,更好地满足大中型水库工程设计中对双峰型洪水过程的实际需求。
附图说明
图1是本发明方法的流程图。
图2是模拟退火算法的流程图。
图3是设计条件下双峰型洪水过程线形状示意图。
具体实施方式
下面通过实施例,并结合附图对本发明作进一步说明。
如图1-图3所示,一种设计条件下双峰型洪水过程线形状推求方法,收集水库历年最大洪水过程线资料,统计得到相应的时段洪量,并归一化得到对应的单位洪水过程线,在采用混合概率密度函数对归一化的单位洪水过程线进行优化拟合得到洪水形状参数系列的基础上,对洪水形状参数系列进行频率分析,计算获得设计条件下的洪水形状参数值,进而推求设计条件下双峰型洪水过程线形状。图1是本实施例的计算流程图,按照以下步骤进行:
1.收集水库历年最大洪水过程线资料。
本具体实施中收集的水库流量数据的时间尺度为△t(一般小于日,如1小时,3小时等)。历年最大洪水过程线的控制时段根据流域暴雨洪水特性和水库调洪能力确定,是一个固定历时。
2.归一化得到对应的单位洪水过程线。
根据步骤1中控制时段水库流量数据,选定历年最大洪水过程线,进而计算出历年最大洪水过程线的时段洪量,本具体实施中控制时段为D=m·△t,有n年最大洪水过程线系列,采用复化梯形公式计算第j年的时段洪量Vj(j=1,2,…,n):
其中,qij为第j年,第i时刻(i=0,1,…,m)的洪水流量,i、j、n和m均为自然数。
由于混合Beta概率密度函数的定义域为[0,1],且纵坐标概率密度函数值与横轴包围的面积为1个单位,为了使历年最大洪水过程线满足混合Beta概率密度函数这些统计特性,需要将历年最大洪水过程线进行归一化得到对应的单位洪水过程线。
本具体实施中将横坐标除以控制时段总历时,纵坐标乘以控制时段总历时与时段洪量之比,计算公式分别如下:
其中,xi为归一化的横坐标,fij为归一化的第j年,第i时刻(i=0,1,…,m)的纵坐标洪水流量。
3.确定混合概率密度函数线型和参数,得到洪水形状参数系列。
本步骤包括两个子步骤:
3.1选择混合概率密度函数线型
由于单位洪水过程线的概率密度线型是未知的,通常选用能较好拟合实测单位洪水过程线的线型。经过水文学者的比较与研究,发现Beta概率密度函数适应性强,拟合使用灵活,能够给出多种可能的形状类型。
本具体实施中采用由2个Beta概率密度函数加权组成的混合Beta概率密度函数作为混合概率密度函数线型,表达式为
f(x)=αf1(x)+(1-α)f2(x)(0≤x≤1) (4)
其中,f(x)为混合Beta概率密度函数;α为权重系数;f1(x)、f2(x)分别为Beta概率密度函数,表达式分别为
其中,B(a,b)为Beta函数,计算公式如下
其中,a1和b1、a2和b2分别为Beta函数的参数。
3.2估计混合概率密度函数线型的参数
在子步骤3.1中选定混合Beta概率密度函数线型后,接着需要进行估计混合Beta概率密度函数线型的参数。每年的单位洪水过程线采用一个混合Beta概率密度函数进行表示和模拟,每个混合Beta概率密度函数共有α、a1、b1、a2和b2等5个参数需要进行估计。
混合分布的参数估计方法主要有非线性优化算法、极大似然法、极大似然的EM算法和最大熵准则法等。模拟退火算法是一种有效的全局优化算法,在组合优化问题的求解中,与以往的近似算法相比,其具有描述简单、使用灵活、运用广泛、运行效率高及受初始条件限制较少等特点。模拟退火算法具有渐近收敛性,已在理论上被证明是一种以概率l收敛于全局最优解的全局优化算法。
本具体实施中采用模拟退火算法估计混合Beta概率密度函数的参数。寻求最优化参数的目标函数为
其中,Ej表示第j年的目标函数值;fij^、fij分别为第j年,第i时刻模拟和实测的单位洪水过程线流量值。
如图2所示,给出了模拟退火算法的流程图。基本步骤如下:
(1)初始化:初始温度T(充分大),初始解状态S(算法迭代的起点),每个T值的迭代次数L;
(2)对k=1,…,L做第3至第6步:
(3)产生新解S′;
(4)计算增量ΔT=C(S′)-C(S),其中C(S)为目标函数;
(5)若ΔT<0则接受S′作为新的当前解,否则按Metropolis准则接受新解,即以概率exp(-ΔT/T)接受S′作为新的当前解;
(6)如果满足终止条件则输出当前解作为最优解,结束程序;终止条件通常取为连续若干个新解都没有被接受时终止算法;
(7)T逐渐减少,且T=T-1且T>0,然后转第2步。
对历年的单位洪水过程线,分别估计其混合Beta概率密度函数的参数,就可以获得样本容量为n的洪水形状参数系列:{αj},{a1j},{b1j},{a2j}和{b2j}(j=1,2,…,n)。
4.对洪水形状参数系列进行频率分析,计算设计条件下的洪水形状参数值。
根据对步骤3中的洪水形状参数系列进行频率分析,选取适当的概率分布函数线型,并估计概率分布函数的参数,计算设计条件下的洪水形状参数值,本步骤包括三个子步骤:
4.1选择概率分布函数线型
由于洪水形状参数的总体分布频率线型是未知的,通常选用能较好拟合洪水形状参数资料系列的线型。经过研究发现,Gumbel分布对于洪水形状参数资料拟合较好。
本具体实施中采用Gumbel分布作为洪水形状参数的概率分布函数线型,其概率分布函数表达式为
其中μ和β为Gumbel分布函数的参数。
4.2估计概率分布函数线型的参数
在子步骤4.1中采用Gumbel分布作为洪水形状参数线型后,接着需要进行估计Gumbel分布的参数。线性矩法是目前国内外公认的有效参数估计方法,最大特点是对序列的极大值和极小值没有常规矩那么敏感,求得的参数估计值比较稳健。
本具体实施中采用线性矩法估计概率分布函数线型的参数。
4.3计算设计条件下的洪水形状参数值
经过子步骤4.1和4.2后就可以得到对应的洪水形状参数的概率分布函数曲线。对于指定的设计重现期T,设计条件下的洪水形状参数值由下式计算:
yT=μ-βln(lnT) (10)
5.推求设计条件下双峰型洪水过程线形状。
将步骤4中设计条件下的洪水形状参数值αT,a1T,b1T,a2T和b2T代入步骤3中的混合Beta概率密度函数式(4)、(5)和(6),推求设计条件下双峰型洪水过程线形状。
设计条件下混合Beta概率密度函数表达式为
f(x)=αTf1(x)+(1-αT)f2(x)(0≤x≤1) (11)
如图3所示,给出了设计条件下双峰型洪水过程线形状示意图。
综上,本发明通过收集水库历年最大洪水过程线资料,统计得到相应的时段洪量,并归一化得到对应的单位洪水过程线,在采用混合概率密度函数对归一化的单位洪水过程线进行优化拟合得到洪水形状参数系列的基础上,对洪水形状参数系列进行频率分析,计算获得设计条件下的洪水形状参数值,进而推求设计条件下双峰型洪水过程线形状。本发明提供了基于混合概率密度函数的设计条件下双峰型洪水过程线形状的推求方法,统计理论基础较强,更加符合设计条件下的暴雨洪水产汇流特性,更好地满足大中型水库工程设计中对双峰型洪水过程的实际需求。
Claims (4)
1.一种设计条件下双峰型洪水过程线形状推求方法,其特征在于包括以下步骤:
步骤1,收集水库历年最大洪水过程线资料;
步骤2,根据步骤1中控制时段水库流量数据,选定历年最大洪水过程线,进而计算出历年最大洪水过程的时段洪量,对历年最大洪水过程线进行归一化得到对应的单位洪水过程线;
步骤3,根据步骤2中的单位洪水过程线,选取混合Beta分布作为混合概率密度函数线型,并估计混合概率密度函数的参数,获得洪水形状参数系列;
步骤4,对步骤3中的洪水形状参数系列进行频率分析,选取Gumbel分布作为洪水形状参数的概率分布函数线型,并估计概率分布函数的参数,计算设计条件下的洪水形状参数值;
步骤5,将步骤4中设计条件下的洪水形状参数值代入步骤3中的混合概率密度函数线型,推求设计条件下双峰型洪水过程线形状。
2.如权利要求1所述一种设计条件下双峰型洪水过程线形状推求方法,其特征在于:所述步骤1中收集水库历年最大洪水过程线资料的时间尺度为△t,所述△t范围为0-24小时。
3.如权利要求2所述一种设计条件下双峰型洪水过程线形状推求方法,其特征在于:所述步骤3中,估计混合概率密度函数的参数采用非线性优化算法、模拟退火算法、极大似然法、以及极大似然的EM算法和最大熵准则法的任意一种或几种组合。
4.如权利要求3所述一种设计条件下双峰型洪水过程线形状推求方法,其特征在于:所述步骤4中,估计概率分布函数的参数采用的方法为线性矩法、权函数法或概率权重矩法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610813782.XA CN106407543B (zh) | 2016-09-09 | 2016-09-09 | 一种设计条件下双峰型洪水过程线形状推求方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610813782.XA CN106407543B (zh) | 2016-09-09 | 2016-09-09 | 一种设计条件下双峰型洪水过程线形状推求方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106407543A CN106407543A (zh) | 2017-02-15 |
CN106407543B true CN106407543B (zh) | 2019-08-20 |
Family
ID=57999500
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610813782.XA Active CN106407543B (zh) | 2016-09-09 | 2016-09-09 | 一种设计条件下双峰型洪水过程线形状推求方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106407543B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109685334B (zh) * | 2018-12-10 | 2020-07-10 | 浙江大学 | 一种新的基于多尺度理论的水文模型模拟评估方法 |
CN111708980A (zh) * | 2020-06-22 | 2020-09-25 | 江西省水利科学研究院 | 一种考虑历史洪水信息的分期设计洪水计算方法 |
CN115510631B (zh) * | 2022-09-15 | 2024-03-19 | 长江水利委员会水文局 | 考虑多场次洪水形态的洪水过程线设计方法及系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102930357A (zh) * | 2012-11-20 | 2013-02-13 | 中铁第四勘察设计院集团有限公司 | 岩溶隧道暗河暴雨涌水洪峰峰值及峰值时间的预测方法 |
CN103276686A (zh) * | 2013-05-24 | 2013-09-04 | 武汉大学 | 一种梯级水库下游设计洪水的推求方法 |
CN104021312A (zh) * | 2014-06-23 | 2014-09-03 | 武汉大学 | 基于过程线形心的洪水过程相似性分析方法及系统 |
CN104182650A (zh) * | 2014-09-11 | 2014-12-03 | 中国电建集团成都勘测设计研究院有限公司 | 推求设计洪水过程线的方法 |
CN104615907A (zh) * | 2015-03-11 | 2015-05-13 | 武汉大学 | 基于多变量条件最可能组合推求设计洪水过程线的方法 |
CN105887751A (zh) * | 2016-04-28 | 2016-08-24 | 武汉大学 | 一种考虑水库调洪演算的多变量洪水重现期计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160055593A1 (en) * | 2014-08-21 | 2016-02-25 | David P. Groeneveld | System and Method to Predict Field Access and the Potential for Prevented Planting Claims for Use by Crop Insurers |
-
2016
- 2016-09-09 CN CN201610813782.XA patent/CN106407543B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102930357A (zh) * | 2012-11-20 | 2013-02-13 | 中铁第四勘察设计院集团有限公司 | 岩溶隧道暗河暴雨涌水洪峰峰值及峰值时间的预测方法 |
CN103276686A (zh) * | 2013-05-24 | 2013-09-04 | 武汉大学 | 一种梯级水库下游设计洪水的推求方法 |
CN104021312A (zh) * | 2014-06-23 | 2014-09-03 | 武汉大学 | 基于过程线形心的洪水过程相似性分析方法及系统 |
CN104182650A (zh) * | 2014-09-11 | 2014-12-03 | 中国电建集团成都勘测设计研究院有限公司 | 推求设计洪水过程线的方法 |
CN104615907A (zh) * | 2015-03-11 | 2015-05-13 | 武汉大学 | 基于多变量条件最可能组合推求设计洪水过程线的方法 |
CN105887751A (zh) * | 2016-04-28 | 2016-08-24 | 武汉大学 | 一种考虑水库调洪演算的多变量洪水重现期计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106407543A (zh) | 2017-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105512767B (zh) | 一种多预见期的洪水预报方法 | |
CN105887751B (zh) | 一种考虑水库调洪演算的多变量洪水重现期计算方法 | |
CN106407543B (zh) | 一种设计条件下双峰型洪水过程线形状推求方法 | |
CN107358060B (zh) | 一种对风电功率预测误差区间进行估计的方法 | |
CN105260803B (zh) | 一种系统用电量预测方法 | |
CN108711847B (zh) | 一种基于编码解码长短期记忆网络的短期风电功率预测方法 | |
CN108428017A (zh) | 基于核极限学习机分位数回归的风电功率区间预测方法 | |
CN103020434A (zh) | 一种基于粒子群优化的最小二乘支持向量机组合预测方法 | |
Zhang et al. | Multivariate and multimodal wind distribution model based on kernel density estimation | |
CN109615120B (zh) | 一种分布式光伏发电出力预测方法 | |
CN107066425A (zh) | 一种变化环境下超定量洪水非一致性分析方法 | |
CN101899820B (zh) | 面向河流生态系统保护的流域地表水可利用量确定方法 | |
CN110378504A (zh) | 一种基于高维Copula技术的光伏发电爬坡事件概率预测方法 | |
CN105631520A (zh) | 一种提取风速随机模糊不确定特征的建模新方法 | |
CN114912671A (zh) | 基于气象预报的气象型骤旱预警方法 | |
CN109299853A (zh) | 一种基于联合概率分布的水库调度函数提取方法 | |
CN107330538A (zh) | 一种气候变化条件下水库适应性调度规则编制的方法 | |
CN112482488B (zh) | 一种城市错峰供水方法及系统 | |
CN112966221A (zh) | 一种南方丰水地区考核年份用水总量折算方法 | |
CN109614642A (zh) | 一种水文频率分析的模型选择方法 | |
CN115510631A (zh) | 考虑多场次洪水形态的洪水过程线设计方法及系统 | |
CN105243446A (zh) | 一种基于粒子群算法的电量组合预测方法 | |
CN115239029A (zh) | 考虑功率时序及气象相依特性的风电功率预测方法及系统 | |
CN113705091B (zh) | 考虑变化环境和调度影响的非一致性设计洪水计算方法 | |
CN109902266A (zh) | 一种基于Copula函数的河道流量演算方法 |
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 |