CN105178240A - Optimization method for drawing P-III type distribution frequency curve - Google Patents
Optimization method for drawing P-III type distribution frequency curve Download PDFInfo
- Publication number
- CN105178240A CN105178240A CN201510310171.9A CN201510310171A CN105178240A CN 105178240 A CN105178240 A CN 105178240A CN 201510310171 A CN201510310171 A CN 201510310171A CN 105178240 A CN105178240 A CN 105178240A
- Authority
- CN
- China
- Prior art keywords
- frequency
- iii type
- flood
- type distribution
- frequency curve
- 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.)
- Granted
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses an optimization method for drawing a P-III type distribution frequency curve. By establishing a P-III type distribution frequency curve optimization model, theoretical frequency and experimental frequency are endowed with different weights, and the minimum square sum of a relative error weight serves as an objective function; according to the drainage basin general situation and the actual condition of a reservoir, constraint conditions of a parameter are added; and the parameter estimated through a moment method serves as an initial value, an improved very fast simulated annealing algorithm is adopted for solving the model, the optimal solution of the parameter is obtained, and therefore the P-III type distribution frequency curve is drawn. By means of the method, the parameter estimation result is unique, randomness caused by human factors is avoided, historical flood and catastrophic flood are also taken into account, the applicability of the P-III type distribution frequency curve is widened, and an important reference basis is provided for drawing the P-III type distribution frequency curve.
Description
Technical field
The present invention relates to a kind of optimization method of frequency curve, particularly a kind of optimization method drawing P-III type distribution frequency curve.
Background technology
Flood Frequency Analysis estimates the main contents of Design of Water Resources and Hydroelectric Projects flood, calculates the important evidence of design standard when being also hydraulic engineering construction.Carry out Design Flood Calculation by building Flood Frequency Analysis model, this model is made up of flood frequency distribution function and method for parameter estimation.The sixties in 20th century, China is according to the practical experience of a large amount of long-term flood serial analysis result and design work for many years, the distribution function that the distribution of regulation P-III type calculates as China's hydrological analysis, adopt the initial value of Moment method estimators parameter, adopt mathematic expectaion formulae discovery empirical Frequency, the numerical value of every flood and corresponding empirical Frequency point are painted on ruled paper, frequency curve is drawn with the first guess of statistical parameter, if curve and empirical Frequency point good not according to matching, adjustable statistical parameter, until curve and point better according to matching, till comparatively satisfied, determine the statistical parameter adopted thus, this method for parameter estimation is called as estimating fitting line method (CF), the Flood Frequency Analysis model of China is called as P-III/CF model.
According to the suitable line method of P-III/CF model, when frequency curve alignment, generally will look after historical flood and extraodinary flood point certificate as much as possible, on the frequency curve namely drawn, afterbody will be tried one's best near it.Because the experience of different project planner, parallax are variant, there is larger arbitrariness in the flood frequency curve adopting the method to obtain, in engineering reality, often occur that different design departments can calculate the situation of different designs result, is difficult to formulate unified Engineering Disign Standard.And, adopt the method manually repeatedly adjusting parameter, also consume a large amount of time and efforts.
(the Liu Li such as Liu Li, Zhou Jianzhong, Yang person of outstanding talent etc., particle group optimizing fits the application of collimation method in hydrologic(al) frequency analysis. the hydrology, 2009, 29 (2), 21-23.) curve-fitting method based on particle swarm optimization algorithm is proposed, and the method is applied in annual flood frequency analysis, although the method avoids the uncertainty of human factor to a certain extent, but the deviation absolute value of the estimate under same frequency that it estimates with the actual value of sample data and theoretical curve with minimum for object function, do not look after historical flood and extraodinary flood, do not meet the actual requirement of Chinese engineering design.(the Li Hongwei such as Li Hongwei, Song Songbai. the application of ant group algorithm in frequency curve parameter calculates. people the Yellow River, 2009,31 (4), 38-40.) a kind of ant group algorithm based on Gridding Method that can be applied to hydrological frequency parameter and calculate is proposed, the method not only parameter is numerous, complex operation, does not look after historical flood and extraodinary flood equally.
Summary of the invention
For the deficiency that prior art exists, the invention provides a kind of optimization method drawing P-III type distribution frequency curve.
The present invention is by setting up P-III type distribution frequency optimization of profile model, different weights is given to the relative error of theoretic frequency and empirical Frequency, minimum for object function with the quadratic sum of relative error weight, and the constraints of parameter is added according to the actual conditions of basin overview and reservoir, with the parameter of Moment method estimators for initial value, adopt very fast simulated annealing algorithm (MVFSA) solving model improved, obtain the optimal solution of parameter, thus draw out P-III type distribution frequency curve.
For solving the problems of the technologies described above, the present invention adopts following technical scheme:
Draw an optimization method for P-III type distribution frequency curve, comprise the steps:
Step 1, collects historical flood and observed flood data;
For the historical flood of collecting, if its maximum investigation phase is N, the historical flood number of investigation is m, and the year number of actual measurement series is n, and containing extraodinary flood number in actual measurement series is l; X
1, X
2..., X
n-l, X
n-l+1..., X
n-l+mthe series of samples of to be sample capacity be n-l+m, and by order arrangement from small to large, use X
jrepresent the jth value in series;
Step 2, calculates the empirical Frequency of Flood evolution, and by the parameter of moments method distribution function according to a preliminary estimate;
For the series of samples X obtained in step 1
j, its empirical Frequency is denoted as P
0(X
j);
Step 3, sets up P-III type distribution frequency optimization of profile model;
Flood data in step 1 are substituted into the frequency distribution function that step 2 obtains, in the hope of different floods point according to the theoretic frequency of correspondence, P can be denoted as
1(X
j), j=1,2 ..., n-l+m, thus the relative error RE (X that can obtain theoretic frequency and empirical Frequency
j), to different floods point certificate, give corresponding weight α
j; Adopt the quadratic sum of relative error weight
minimum is object function, and adds the constraints of parameter according to the actual conditions of basin overview and reservoir, sets up P-III type distribution frequency optimization of profile model;
Step 4, adopts the model in very fast simulated annealing algorithm (MVFSA) solution procedure 3 improved, thus obtains the parameters of P-III type distribution function, thus draws frequency curve.
The very fast simulated annealing algorithm (MVFSA) of described improvement is prior art, specifically see Tongji University's journal (natural science edition) the 34th volume the 8th phase in 2006, and 1121-1125 pages, author is Chen Huagen, Li Lihua, Xu Huiping, Chen Bing.
Described weight α
j, determine by the different interval importance of frequency curve, the sort order of flood sample data is larger, the weight just larger α of this flood point certificate
j> 0, and
Compared with prior art, the present invention has the following advantages and beneficial effect:
The present invention is by setting up P-III type distribution frequency optimization of profile model, different weights is given to the relative error of theoretic frequency and empirical Frequency, minimum for object function with the quadratic sum of relative error weight, and the constraints of parameter is added according to the actual conditions of basin overview and reservoir, with the parameter of Moment method estimators for initial value, adopt the very fast simulated annealing algorithm solving model improved, obtain the optimal solution of parameter, thus draw out P-III type distribution frequency curve.The method parameter estimation result is unique, avoids the arbitrariness that human factor produces, can also look after historical flood and extraodinary flood, expand the applicability of P-III type distribution curve, provides important reference frame for drawing P-III type distribution frequency curve.
Accompanying drawing explanation
Fig. 1 is the flow chart of the inventive method;
Fig. 2 is the frequency curve adopting the inventive method to obtain.
Detailed description of the invention
The present invention is by setting up P-III type distribution frequency optimization of profile model, different weights is given to the relative error of theoretic frequency and empirical Frequency, minimum for object function with the quadratic sum of relative error weight, and the constraints of parameter is added according to the actual conditions of basin overview and reservoir, with the parameter of Moment method estimators for initial value, adopt very fast simulated annealing algorithm (MVFSA) solving model improved, obtain the optimal solution of parameter, thus draw out P-III type distribution frequency curve.Below in conjunction with Fig. 1, concrete steps of the present invention are described:
Step 1, collects historical flood and observed flood data;
For the historical flood of collecting, if its maximum investigation phase is N, the historical flood number of investigation is m, and the year number of actual measurement series is n, and containing extraodinary flood number in actual measurement series is l.X
1, X
2..., X
m+l, X
m+l+1..., X
n-l+mthe series of samples of to be sample capacity be n-l+m, and by order arrangement from big to small, use X
jrepresent the jth value in series.
Step 2, calculates the empirical Frequency of Flood evolution, and by the parameter of moments method distribution function according to a preliminary estimate;
The empirical frequency formula that the present invention recommends is the empirical frequency formula not connecting sequence sample:
In formula, P
rrepresent the empirical Frequency of r field historical flood, P
jrepresent the empirical Frequency of jth field observed flood.
For the series of samples obtained in step 1, for X
j, its empirical Frequency is denoted as P
0(X
j).
Step 3, sets up P-III type distribution frequency optimization of profile model;
Flood data in step 1 are substituted into the frequency distribution function that step 2 obtains, in the hope of different floods point according to the theoretic frequency of correspondence, P can be denoted as
1(X
j), j=1,2 ..., n-l+m, thus the relative error RE (X that can obtain theoretic frequency and empirical Frequency
j), to different floods point certificate, give corresponding weight α
j.Adopt the quadratic sum of relative error weight
minimum is object function, and adds the constraints of parameter according to the actual conditions of basin overview and reservoir, so set up P-III type distribution frequency optimization of profile model.
Relative error
about weight α in Optimized model
jassignment method can adopt experience assignment method, the Bayesian model method of average, arithmetic mean method or fuzzy reasoning method etc.Generally, for weight α
j, determine by the different interval importance of frequency curve, the sort order of flood sample data is larger, the weight just larger α of this flood point certificate
j> 0, and
For constraints, the regulation of different basins to P-III distribution function parameters is distinguished to some extent, but generally requires C
v/ C
s>=2.This detailed description of the invention only adds this point constraints, and the constraints continuing to add other in practical application does not exceed protection scope of the present invention.
Step 4, with the parameter estimated in step 2 for initial value, adopts the model in very fast simulated annealing algorithm (MVFSA) solution procedure 3 improved, thus obtains the parameters of P-III type distribution function, thus draw frequency curve.
Adopting that very fast simulated annealing algorithm (MVFSA) solves local optimum problem is the ordinary skill in the art, and it will not go into details.
Fig. 2 gives the schematic diagram of the frequency curve adopting the present invention to obtain, and the upper afterbody matching of frequency curve is good as we can see from the figure, meets the actual demand of Chinese engineering.
To sum up, the present invention is by setting up P-III type distribution frequency optimization of profile model, different weights is given to the relative error of theoretic frequency and empirical Frequency, minimum for object function with the quadratic sum of relative error weight, and the constraints of parameter is added according to the actual conditions of basin overview and reservoir, with the parameter of Moment method estimators for initial value, adopt very fast simulated annealing algorithm (MVFSA) solving model improved, obtain the optimal solution of parameter, thus draw out P-III type distribution frequency curve.The method parameter estimation result is unique, avoids the arbitrariness that human factor produces, can also look after historical flood and extraodinary flood, expand the applicability of P-III type distribution curve, provides important reference frame for drawing P-III type distribution frequency curve.
Claims (5)
1. draw an optimization method for P-III type distribution frequency curve, it is characterized in that comprising the steps:
Step (1), collects historical flood and observed flood data;
Step (2), calculates the empirical Frequency of Flood evolution, and by the parameter of moments method distribution function according to a preliminary estimate;
Step (3), sets up P-III type distribution frequency optimization of profile model;
Step (4), with the parameter estimated in step (2) for initial value, adopt the model in the very fast simulated annealing algorithm solution procedure (3) improved, obtain the parameters of P-III type distribution function, thus draw frequency curve.
2. a kind of optimization method drawing P-III type distribution frequency curve as claimed in claim 1, it is characterized in that: in step (1), for the historical flood of collecting, if its maximum investigation phase is N, the historical flood number of investigation is m, the year number of actual measurement series is n, and containing extraodinary flood number in actual measurement series is l; X
1, X
2..., X
n-l, X
n-l+1..., X
n-l+mthe series of samples of to be sample capacity be n-l+m, and by order arrangement from small to large, use X
jrepresent the jth value in series.
3. a kind of optimization method drawing P-III type distribution frequency curve as claimed in claim 1, it is characterized in that: in step (2), described empirical frequency formula is the empirical frequency formula not connecting sequence sample; For the series of samples X obtained in step (1)
j, its empirical Frequency is denoted as P
0(X
j).
4. a kind of optimization method drawing P-III type distribution frequency curve as claimed in claim 1, it is characterized in that: in step (3), flood data in step (1) are substituted into the frequency distribution function that step (2) obtains, try to achieve different flood points according to the theoretic frequency of correspondence, be denoted as P
1(X
j), j=1,2 ..., n-l+m, thus the relative error RE (X obtaining theoretic frequency and empirical Frequency
j), to different floods point certificate, give corresponding weight α
j; Adopt the quadratic sum of relative error weight
minimum is object function, and adds the constraints of parameter according to the actual conditions of basin overview and reservoir, sets up P-III type distribution frequency optimization of profile model.
5. a kind of optimization method drawing P-III type distribution frequency curve as claimed in claim 4, is characterized in that: described weight α
j, determine by the different interval importance of frequency curve, the sort order of flood sample data is larger, the weight just larger α of this flood point certificate
j> 0, and
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510310171.9A CN105178240B (en) | 2015-06-08 | 2015-06-08 | A kind of optimization method of drafting P type III distribution frequency curves |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510310171.9A CN105178240B (en) | 2015-06-08 | 2015-06-08 | A kind of optimization method of drafting P type III distribution frequency curves |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105178240A true CN105178240A (en) | 2015-12-23 |
CN105178240B CN105178240B (en) | 2017-11-28 |
Family
ID=54900630
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510310171.9A Expired - Fee Related CN105178240B (en) | 2015-06-08 | 2015-06-08 | A kind of optimization method of drafting P type III distribution frequency curves |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105178240B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105887751A (en) * | 2016-04-28 | 2016-08-24 | 武汉大学 | Multi-variable flood recurrence period calculating method considering reservoir flood routing calculation |
CN107066425A (en) * | 2017-03-17 | 2017-08-18 | 中山大学 | Overdetermination amount flood nonuniformity analysis method under a kind of changing environment |
CN111210141A (en) * | 2020-01-03 | 2020-05-29 | 黄河勘测规划设计研究院有限公司 | Reservoir capacity curve correction method based on constraint mechanism particle swarm algorithm |
CN112085298A (en) * | 2020-09-23 | 2020-12-15 | 中国电建集团成都勘测设计研究院有限公司 | Non-continuous sequence flood frequency analysis method considering historical flood |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5753818A (en) * | 1995-05-15 | 1998-05-19 | North American Geotechnical Company | Method and apparatus for measuring scour around bridge foundations |
CN103882827A (en) * | 2014-04-14 | 2014-06-25 | 武汉大学 | Reservoir flood control risk rate prediction method based on runoff ensemble forecasting |
CN104182650A (en) * | 2014-09-11 | 2014-12-03 | 中国电建集团成都勘测设计研究院有限公司 | Method for inquiring into design flood hygrograph |
CN104615907A (en) * | 2015-03-11 | 2015-05-13 | 武汉大学 | Method for deriving and designing flood process line based on multi-variable most possible condition combination |
CN104674748A (en) * | 2014-12-23 | 2015-06-03 | 南京南瑞集团公司 | Quantitative method for water abandon probability of cascade hydroelectric station groups |
-
2015
- 2015-06-08 CN CN201510310171.9A patent/CN105178240B/en not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5753818A (en) * | 1995-05-15 | 1998-05-19 | North American Geotechnical Company | Method and apparatus for measuring scour around bridge foundations |
CN103882827A (en) * | 2014-04-14 | 2014-06-25 | 武汉大学 | Reservoir flood control risk rate prediction method based on runoff ensemble forecasting |
CN104182650A (en) * | 2014-09-11 | 2014-12-03 | 中国电建集团成都勘测设计研究院有限公司 | Method for inquiring into design flood hygrograph |
CN104674748A (en) * | 2014-12-23 | 2015-06-03 | 南京南瑞集团公司 | Quantitative method for water abandon probability of cascade hydroelectric station groups |
CN104615907A (en) * | 2015-03-11 | 2015-05-13 | 武汉大学 | Method for deriving and designing flood process line based on multi-variable most possible condition combination |
Non-Patent Citations (3)
Title |
---|
张静怡等: "极值分布和P-III型分布线性矩法在区域洪水频率分析中的检验", 《水文》 * |
梁忠民等: "《水文水利计算》", 31 August 2006, 中国水利水电出版社 * |
陈子全等: "基于模拟退火算法的皮尔逊III型分布参数估计", 《人民黄河》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105887751A (en) * | 2016-04-28 | 2016-08-24 | 武汉大学 | Multi-variable flood recurrence period calculating method considering reservoir flood routing calculation |
CN107066425A (en) * | 2017-03-17 | 2017-08-18 | 中山大学 | Overdetermination amount flood nonuniformity analysis method under a kind of changing environment |
CN107066425B (en) * | 2017-03-17 | 2020-08-14 | 中山大学 | Method for analyzing non-uniformity of ultra-quantitative flood in changing environment |
CN111210141A (en) * | 2020-01-03 | 2020-05-29 | 黄河勘测规划设计研究院有限公司 | Reservoir capacity curve correction method based on constraint mechanism particle swarm algorithm |
CN111210141B (en) * | 2020-01-03 | 2023-07-14 | 黄河勘测规划设计研究院有限公司 | Reservoir capacity curve correction method based on constraint machine granulation subgroup algorithm |
CN112085298A (en) * | 2020-09-23 | 2020-12-15 | 中国电建集团成都勘测设计研究院有限公司 | Non-continuous sequence flood frequency analysis method considering historical flood |
CN112085298B (en) * | 2020-09-23 | 2023-06-16 | 中国电建集团成都勘测设计研究院有限公司 | Non-sequential flood frequency analysis method considering historical flood |
Also Published As
Publication number | Publication date |
---|---|
CN105178240B (en) | 2017-11-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105887751B (en) | A kind of multivariable return period of flood computational methods for considering reservoir routing | |
CN105808868B (en) | A kind of hydrological models combination Uncertainty Analysis Method based on Copula function | |
CN105178240A (en) | Optimization method for drawing P-III type distribution frequency curve | |
CN107423546A (en) | Multivariable hydrological uncertainty processing method based on Copula functions | |
CN106875060B (en) | Flood real-time correction method based on global automatic optimization algorithm | |
CN105512767A (en) | Flood forecasting method of multiple forecast periods | |
CN104615907A (en) | Method for deriving and designing flood process line based on multi-variable most possible condition combination | |
CN107622162A (en) | A kind of rating curve calculation method based on Copula functions | |
CN107526880B (en) | Method for calculating upstream river-entering pollution flux based on river water quality observation | |
CN104615845B (en) | A kind of overcritical risk analysis on flood method | |
CN104727264A (en) | Damsite flood restoration method based on Copula function | |
CN101899820B (en) | Method for determining amount of available surface water of river basin facing to protection of river ecosystem | |
CN110276150A (en) | A kind of Mountain Area river basal flow capacity system interpolation extension method based on Copula function | |
CN108897940B (en) | Rectangular grid-based one-way coupling method for distributed hydrological model and two-dimensional hydrodynamic model | |
CN106192863A (en) | A kind of power station installed capacity and annual electricity generating capacity evaluation method | |
CN104021307B (en) | Measurement method for backstepping of reservoir water inlet flow | |
CN108875130A (en) | A kind of design flood peak amount condition most probable combined method based on Copula function | |
CN109902267B (en) | River channel safety discharge amount calculation method influenced by downstream lake jacking | |
CN104766245B (en) | A kind of cable load curve evaluation method theoretical based on weather information and entropy weight | |
CN106320255A (en) | Overbank flow floodplain and channel excess flow computing method | |
CN104091227A (en) | Method for analyzing load density of power consumers in city | |
CN106969756A (en) | River suspended sediment sedimentary loading correction computational methods | |
Hansen et al. | Flow forecasting in urban drainage systems using deterministic updating of water levels in distributed hydraulic models | |
CN104751008A (en) | L-moment parameter estimating method for fuzzy information containing historical flood | |
CN106557614A (en) | A kind of Flood Frequency Analysis method being distributed based on Halphen IB |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171128 Termination date: 20180608 |