CN109271664B - Design rainstorm safety analysis method based on assumed conditions of consistency and non-consistency - Google Patents

Design rainstorm safety analysis method based on assumed conditions of consistency and non-consistency Download PDF

Info

Publication number
CN109271664B
CN109271664B CN201810877331.1A CN201810877331A CN109271664B CN 109271664 B CN109271664 B CN 109271664B CN 201810877331 A CN201810877331 A CN 201810877331A CN 109271664 B CN109271664 B CN 109271664B
Authority
CN
China
Prior art keywords
rainstorm
design
value
consistency
sequence
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
CN201810877331.1A
Other languages
Chinese (zh)
Other versions
CN109271664A (en
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.)
South China University of Technology SCUT
Original Assignee
South China University of Technology SCUT
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 South China University of Technology SCUT filed Critical South China University of Technology SCUT
Priority to CN201810877331.1A priority Critical patent/CN109271664B/en
Publication of CN109271664A publication Critical patent/CN109271664A/en
Application granted granted Critical
Publication of CN109271664B publication Critical patent/CN109271664B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Abstract

The invention discloses a design rainstorm safety analysis method based on assumed conditions of consistency and non-consistency, which comprises the following steps: extracting an annual maximum value sequence from a perennial daily precipitation sequence; extracting a series of subsequences from the annual maximum value sequence by adopting a 30-year window sliding method and a 30-year window increasing method; fitting each subsequence through generalized extremum distribution and calculating a design rainstorm value sequence; calculating a trend value of a designed rainstorm sequence by adopting a linear regression method; and (4) making a difference value of the trend values of the designed rainstorm sequence under the two extraction methods, and carrying out safety analysis according to the difference value. The invention provides a new method for researching the difference of the non-uniformity and the uniformity hypothesis to the design rainstorm, compares the annual difference of the design rainstorm calculated under the two hypotheses, is favorable for evaluating the safety of the design standard, and can provide reference significance for hydraulic engineering construction, water resource and safety management, flood control and disaster reduction.

Description

Design rainstorm safety analysis method based on assumed conditions of consistency and inconsistency
Technical Field
The invention relates to the field of hydrological meteorology and hydrological statistics, in particular to a designed rainstorm safety analysis method based on assumed conditions of consistency and non-consistency.
Background
Due to the influence of global climate change, many hydrological and meteorological elements show obvious variation with time, and precipitation is one of the elements with variation. In actual hydraulic engineering construction, design rainstorm is generally calculated by a maximum annual method and a hydrological frequency curve based on consistency assumption, and typical hydrological frequency curves include generalized extreme value distribution (GEV), pearson type iii, generalized normal distribution, generalized logistic distribution and the like. Therefore, the existence of the inconsistency can cause the difference between the actually adopted design value and the expected value, and the discussion of the potential difference has important significance on the safety of the hydraulic engineering. In the related research for designing rainstorm differences under consistent and inconsistent assumptions, most of the research focuses on discussing the calculation method and the variation trend of the inconsistent rainstorm design value. To investigate the effect of non-uniformity, there are probably two approaches: (1) considering the statistical parameters of the distribution function as functions of time or other factors (such as precipitation, air temperature and other factors), although the method can discuss the change trend of historical and future design rainstorm, the type of the distribution of the function met by the parameters needs to be assumed in advance; (2) the influence of non-uniformity is weakened by adopting a fixed year window sliding method, such as a 30-year sliding window method, namely, a hydrological frequency curve is fitted to a precipitation sequence which is continuous for 30 years each time, which is different from the traditional method of directly fitting all precipitation sequences, and the method is suitable for discussing the change trend of historical design rainstorm. However, for the differences in designed rainstorms under the consistent and inconsistent assumptions, most studies consider the designed rainstorm values, which are derived based on the consistent assumptions, as fixed values, i.e. fitting the hydrological frequency curves to all available precipitation sequences, ignoring their effect as the observed years increase (precipitation sequences grow).
Disclosure of Invention
The invention aims to overcome the defects and shortcomings in the prior art, provides a design rainstorm safety analysis method based on the assumed conditions of consistency and inconsistency, provides reference for the safety of large-scale hydraulic engineering design, and further provides scientific and practical values for hydraulic engineering construction, water resource and safety management, flood control and disaster reduction.
In order to realize the purpose, the invention adopts the following technical scheme:
the design rainstorm safety analysis method based on the assumed conditions of consistency and inconsistency comprises the following steps:
s1, extracting an annual maximum value sequence X ═ X from the annual daily precipitation sequence 1 ,X 2 ,...,X i ,...,X n Is given, wherein n is the sequence length, X i Represents the maximum rainfall in the ith year in the sequence, 0<i<N, and i ∈ N +
S2, extracting n-29 subsequences from the annual maximum value sequence by respectively adopting a 30-year window sliding method and a 30-year window increasing method, wherein the general formula of the set of the subsequences is S ═ S { (S) } 1 ,S 2 ,...,S j ,...,S n-29 In which S is j Denotes the jth group of subsequences, 0<j<N-29, and j ∈ N +
S3, adopting generalized extremum distribution to fit each subsequence of the set S, and calculating the subsequence S j Designing a rainstorm value in a given reproduction period, and further solving a designed rainstorm sequence;
s4, performing trend analysis on the designed rainstorm sequence obtained in the step S3 by using a linear regression method, and calculating the annual change trend of the designed rainstorm value;
s5, extracting the subsequence S based on the 30-year window sliding method and the 30-year window increasing method respectively j The annual variation trends of the design rainstorm values under the two extraction methods are obtained according to the annual variation trend of the design rainstorm values obtained in step S4, the two trends are subtracted, and the safety of the design rainstorm values calculated under the conditions of consistency assumption and non-consistency assumption is analyzed according to the trend difference.
Preferably, in step S2, for n-29 subsequences extracted by 30-year window sliding method, each subsequence has 30 consecutive years of data, and the subsequence S extracted by 30-year window sliding method j Can be expressed as: s j ={X j ,X j+1 ,...,X j+29 },0<j<N-29 and j ∈ N +
For n-29 sub-sequences extracted by a 30-year window-growing method, the initial sequence is the same as the initial sequence extracted by a 30-year window-sliding method, and then each sub-sequence is gradually increased by one value, so that the sub-sequence S extracted by the 30-year window-growing method j Can be expressed as: s j ={X 1 ,X 2 ,...,X j+29 },0<j<N-29 and j ∈ N +
As a preferable embodiment, in step S2, the assumed non-conformity and conformity conditions are simulated by using a 30-year window sliding method and a 30-year window increasing method.
As a preferred technical solution, in step S3, each subsequence of the generalized extremum distribution fitting set S is adopted, a maximum likelihood method is adopted to estimate parameters of the generalized extremum distribution, and then a design rainstorm value is calculated according to a given rainstorm recurrence period, wherein each subsequence obtains the design rainstorm value in the given recurrence period, and since the subsequences are extracted respectively based on a 30-year window sliding method and a 30-year window increasing method, two groups of design rainstorm sequences are obtained in each recurrence period, and each group of sequences contains n-29 design rainstorm values.
As a preferred technical solution, in step S3, the specific process is as follows:
the generalized extreme value distribution is based on an extreme value theory, has 3 parameters in total, and is used for quantitatively designing the intensity and the frequency of the rainstorm; the generalized extremum distribution function is as follows:
Figure GDA0003814793790000031
wherein the content of the first and second substances,
Figure GDA0003814793790000032
x represents time, and μ, σ and k represent a position parameter, a scale parameter and a shape parameter, respectively;
calculating by adopting a maximum likelihood estimation method to obtain 3 parameters of the generalized extreme value distribution function, and further finishing the quantification of designing the rainstorm according to the generalized extreme value distribution function; for each subsequence S j Given a design rainstorm in T years, the quantified design rainstorm value is:
Figure GDA0003814793790000033
for a given recurrence period, i.e., T years of design rainstorm, set S ═ S by subsequence 1 ,S 2 ,...,S j ,...,S n-29 And subsequence S j Calculating the design rainstorm value to obtain a design rainstorm value sequence:
Figure GDA0003814793790000037
wherein 0<j<N-29 and j ∈ N +
Preferably, the step S3 is based on the 30-year window sliding method and the 30-year windowSubsequence S extracted by mouth growth method j In contrast, the generalized extremum distribution is adopted for the subsequence S j When fitting is carried out, fitting parameters obtained by calculation through a maximum likelihood estimation method are different, and the quantized design rainstorm value y j T In contrast, two different design rainstorm values can be obtained for each given recurrence period T, and two different sets of design rainstorm sequences can be obtained.
As a preferable technical solution, in step S4, a linear regression equation is used to calculate the variation trend of the design rainstorm value, where:
Figure GDA0003814793790000034
Figure GDA0003814793790000035
Figure GDA0003814793790000036
wherein
Figure GDA0003814793790000041
And
Figure GDA0003814793790000042
respectively representing intercept and slope, x is time, y is precipitation, x i Which represents the time of the i-th year,
Figure GDA0003814793790000043
arithmetic mean, y, representing the year sequence i Indicates the precipitation amount of the i-th year,
Figure GDA0003814793790000044
denotes y i The arithmetic mean of (a); the annual variation trend of the design rainstorm values can be determined by designing the rainstorm sequence y T The following formula is used to obtain:
Figure GDA0003814793790000045
wherein
Figure GDA0003814793790000046
Is a sequence y T The jth group of (a) design storm values,
Figure GDA0003814793790000047
is a sequence y T Is calculated as the arithmetic mean of the average of the values,
Figure GDA0003814793790000048
and designing an annual trend value of the rainstorm for the first time of T years.
As a preferable technical solution, in step S4, when trend analysis is performed on the designed rainstorm sequence, F distribution is used to check the significance of the trend, and the confidence level is 0.05.
Preferably, in step S5, the subsequence S is extracted based on the 30-year window sliding method and the 30-year window growing method, respectively j By quantified design rainstorm values
Figure GDA0003814793790000049
And the annual change trend of the design storm value
Figure GDA00038147937900000410
Two trend values can be obtained and are respectively recorded as
Figure GDA00038147937900000411
The difference between the two is:
Figure GDA00038147937900000412
wherein the content of the first and second substances,
Figure GDA00038147937900000413
indicating a design rainstorm value trend value based on a 30-year window sliding method,
Figure GDA00038147937900000414
represents a design rainstorm annual trend value based on a 30-year window growth law,
Figure GDA00038147937900000415
representing the difference of the two trend values;
if it is
Figure GDA00038147937900000416
That is, the design rainstorm value calculated based on the non-consistency assumption is smaller than the design rainstorm value calculated based on the consistency assumption, which indicates that the design rainstorm value calculated by adopting the consistency assumption is safer for engineering construction in the long run; if it is
Figure GDA00038147937900000417
It indicates that in the long run, it is not advantageous to design a rainstorm method using the consensus hypothesis calculation.
Compared with the prior art, the invention has the following advantages and effects:
the invention provides a new analysis method for discussing the design rainstorm safety under the assumption of non-uniformity and consistency, and the 30-year window sliding method and the 30-year window increasing method are used for respectively simulating the non-uniformity and the consistency assumption to plan the rainstorm, thereby not only representing the annual change of the non-uniformity to the design rainstorm, but also considering the influence of the extension of a precipitation sequence to the consistency assumption to plan the rainstorm, overcoming the defect that the traditional engineering design ignores the non-uniformity of the precipitation sequence, and perfecting the defect that the influence of the increase of the precipitation sequence to the design rainstorm value in the existing research. By comparing the difference of the annual change trend of the design rainstorm under two assumptions, the method is beneficial to evaluating the change of the actual standard of related infrastructure along with the increase of the service life, and evaluating the influence (namely the safety) of the change on the normal operation of the infrastructure on the basis, so that the method can provide reference significance for hydraulic engineering construction, water resource and safety management, flood control and disaster reduction.
Drawings
FIG. 1 is a flow chart of a design storm safety analysis method based on assumed conditions of consistency and inconsistency according to the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not limiting.
Examples
As shown in fig. 1, the method for analyzing the safety of design rainstorm based on assumed conditions of consistency and inconsistency comprises the following steps:
s1, extracting an annual maximum value sequence X ═ X from the annual daily precipitation sequence 1 ,X 2 ,...,X i ,...,X n Is given, wherein n is the sequence length, X i Represents the maximum rainfall in the ith year in the sequence, 0<i<N, and i ∈ N +
S2, adopting a 30-year window sliding method and a 30-year window increasing method to simulate assumed conditions of non-consistency and consistency, respectively adopting the 30-year window sliding method and the 30-year window increasing method to extract n-29 subsequences from the annual maximum value sequence, wherein the general set formula of the subsequences is S ═ { S { (S) } 1 ,S 2 ,...,S j ,...,S n-29 In which S is j Denotes the jth group of subsequences, 0<j<N-29, and j ∈ N +
For n-29 groups of subsequences extracted by a 30-year window sliding method, each subsequence has 30 continuous years of data, and the subsequence S extracted by the 30-year window sliding method j Can be expressed as: s j ={X j ,X j+1 ,...,X j+29 },0<j<N-29 and j ∈ N +
For n-29 groups of subsequences extracted by the 30-year window increasing method, the initial sequence is the same as the initial sequence extracted by the 30-year window sliding method, and then each subsequence is increased by one value, so that the subsequence S extracted by the 30-year window increasing method j Can be expressed as: s j ={X 1 ,X 2 ,...,X j+29 },0<j<N-29 and j ∈ N +
S3, adopting generalized extremum distribution to fit each subsequence of the set S, and calculating the subsequence S j Designing a rainstorm value in a given reappearance period, and further solving a designed rainstorm sequence; the specific process is as follows:
generalized extreme value distribution (GEV) is based on an extreme value theory, has 3 parameters in total, and is used for quantitatively designing the intensity and frequency of the rainstorm; the generalized extremum distribution function is as follows:
Figure GDA0003814793790000051
wherein the content of the first and second substances,
Figure GDA0003814793790000061
x represents time, and μ, σ and k represent a position parameter, a scale parameter and a shape parameter, respectively;
calculating by adopting a maximum likelihood estimation method to obtain 3 parameters of the generalized extreme value distribution function, and further finishing the quantification of the designed rainstorm according to the generalized extreme value distribution function; for each subsequence S j Given a design rainstorm in T years, the quantified design rainstorm value is:
Figure GDA0003814793790000062
for a given recurrence period, i.e., T years of design rainstorm, set S ═ S by subsequence 1 ,S 2 ,...,S j ,...,S n-29 And subsequence S j Calculating the design rainstorm value to obtain a design rainstorm value sequence:
Figure GDA0003814793790000063
wherein 0<j<N-29 and j ∈ N +
S4, performing trend calculation on the designed rainstorm sequence obtained in the step S3 by using a linear regression method to obtain the annual variation trend of the designed rainstorm value; the specific process is as follows:
calculating the change trend of the designed rainstorm value by adopting linear regression, wherein the linear regression equation is as follows:
Figure GDA0003814793790000064
Figure GDA0003814793790000065
Figure GDA0003814793790000066
wherein
Figure GDA0003814793790000067
And
Figure GDA0003814793790000068
respectively representing intercept and slope, x is time, y is precipitation, x i Which represents the time of the i-th year,
Figure GDA0003814793790000069
arithmetic mean, y, representing the year sequence i Indicates the precipitation amount of the i-th year,
Figure GDA00038147937900000610
denotes y i The arithmetic mean of (a); the annual variation trend of the design storm values can be determined by designing the storm sequence y T The following formula is used to obtain:
Figure GDA00038147937900000611
wherein
Figure GDA00038147937900000612
Is a sequence y T The jth group of (a) design storm values,
Figure GDA00038147937900000613
is a sequence y T Is calculated as the arithmetic mean of the average of the values,
Figure GDA00038147937900000614
designing an annual trend value of rainstorm for T year; trend significance was tested using an F-distribution (i.e., F-distribution, presented in 1924 by the statistician r.a. fisher), with a confidence level of 0.05.
S5, and subsequence S extracted respectively based on 30-year window sliding method and 30-year window growing method j Obtaining the annual variation trend of the designed rainstorm value under two extraction methods according to the annual variation trend of the designed rainstorm value obtained in S4, making a difference between the two trends, and analyzing the safety of the designed rainstorm value calculated under the conditions of consistency assumption and non-consistency assumption according to the trend difference; the specific process is as follows:
subsequence S respectively extracted based on 30-year window sliding method and 30-year window growing method j Design of rainstorm values by quantification
Figure GDA0003814793790000071
And annual variation trend of design storm values
Figure GDA0003814793790000072
Two trend values can be obtained and are respectively recorded as
Figure GDA0003814793790000073
Figure GDA0003814793790000074
The difference between the two is:
Figure GDA0003814793790000075
wherein the content of the first and second substances,
Figure GDA0003814793790000076
represents a design rainstorm annual trend value based on a 30-year window sliding method,
Figure GDA0003814793790000077
represents a design rainstorm annual trend value based on a 30-year window growth law,
Figure GDA0003814793790000078
representing the difference of the two trend values;
if it is
Figure GDA0003814793790000079
That is, the design rainstorm value calculated based on the non-consistency assumption is smaller than the design rainstorm value calculated based on the consistency assumption, which indicates that the design rainstorm value calculated by adopting the consistency assumption is safer for engineering construction in the long run; if it is
Figure GDA00038147937900000710
It indicates that in the long run, it may be disadvantageous to design a rainstorm method using the consensus hypothesis calculation.
The following is a specific application example of the method of the present invention:
based on the 1961-2013-year Chinese grid daily rainfall data set, a 30-year window sliding method and a 30-year window increasing method are adopted to simulate the assumption of non-consistency and consistency, the historical evolution trend of designing rainstorm of 3825 grid points is discussed, and the difference of calculation results of the two methods in China is explored.
The specific calculation steps are as follows:
(1) extracting a annual maximum daily precipitation sequence (total 53 years from 1961 to 2013) from daily precipitation data of 3825 grid points;
(2) respectively extracting 24 subsequences from the annual maximum value sequence by adopting a 30-year window sliding method and a 30-year window growing method;
(3) and fitting each subsequence by adopting GEV distribution, estimating GEV parameters by adopting a maximum likelihood method, calculating design rainstorm values by taking 20-year-one-chance, 50-year-one-chance and 100-year-one-chance in the recurrence period as examples, and obtaining 3 design rainstorms in the recurrence period for each subsequence. Thus, based on the 30-year window sliding method and the 30-year window growing method, two sets of designed rainstorm sequences were obtained for each recurrence period, each set of sequences containing 24 values.
(4) Based on design violenceAnd 4, calculating a designed rainstorm change trend by adopting a linear regression equation in the rain value sequence. The spatial distribution of 3825 grid point trends in the country is the annual change trend values of 20-year-first, 50-year-first and 100-year-first obtained by adopting a 30-year window sliding method
Figure GDA0003814793790000081
A value; adopting 30-year window increasing method to obtain annual change trend values of 20-year-one, 50-year-one and 100-year-one
Figure GDA0003814793790000082
A value; grid points that failed the significance test (significance of 0.05) are represented by black scatter dots, and trend values (units: mm/yr, i.e., mm/year) are represented by slope.
(5) And (5) calculating the difference of the two change trends according to the formula (5), and analyzing the safety of the two methods for hydraulic engineering design.
Utilizing the difference of annual variation trend of 20-year-one, 50-year-one and 100-year-one
Figure GDA0003814793790000083
And the calculation result shows that the difference of the rainstorm trend has obvious spatial difference. In general, most grid point trend differences in northeast, south and northwest China
Figure GDA0003814793790000084
The trend calculated by the 30-year window sliding method is greater than that calculated by the 30-year window increasing method, which indicates that the designed rainstorm value calculated by the non-uniformity assumption is higher than that calculated by the uniformity scenario, and the traditional calculation method is not beneficial to large-scale hydraulic engineering design; trend difference of most grid points in northern and southwest China
Figure GDA0003814793790000085
Namely, the trend calculated by the 30-year window sliding method is smaller than that calculated by the 30-year window increasing method, which shows that the designed rainstorm value calculated by the non-uniformity assumption is lower than that calculated by the uniformity situation, and the traditional calculation method is safer for hydraulic engineering design.
The above-mentioned embodiments only express several embodiments of the present invention, and the description thereof is specific and detailed, but not to be understood as limiting the scope of the present invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the inventive concept, which falls within the scope of the present invention. Therefore, the protection scope of the present patent shall be subject to the claims.

Claims (7)

1. The method for analyzing the safety of the design rainstorm based on the assumed conditions of consistency and non-consistency is characterized by comprising the following steps of:
s1, extracting an annual maximum value sequence X ═ X from the annual daily precipitation sequence 1 ,X 2 ,...,X i ,...,X n Where n is the sequence length, X i Represents the maximum rainfall of the ith year in the sequence, 0 < i < N, and i belongs to N +
S2, extracting n-29 subsequences from the annual maximum value sequence by respectively adopting a 30-year window sliding method and a 30-year window increasing method, wherein the general formula of the set of the subsequences is S ═ S 1 ,S 2 ,...,S j ,...,S n-29 In which S is j Represents the jth group of subsequences, 0 < j < N-29, and j ∈ N +
S3, adopting generalized extremum distribution to fit each subsequence of the set S, and calculating the subsequence S j Designing a rainstorm value in a given reappearance period, and further solving a designed rainstorm sequence;
s4, performing trend analysis on the designed rainstorm sequence obtained in the step S3 by using a linear regression method, and calculating the annual change trend of the designed rainstorm value;
s5, extracting the subsequence S based on the 30-year window sliding method and the 30-year window growing method respectively j Obtaining the annual variation trend of the designed rainstorm value under two extraction methods according to the annual variation trend of the designed rainstorm value obtained in S4, subtracting the two trends, and analyzing the safety of the designed rainstorm value calculated under the conditions of consistency assumption and non-consistency assumption according to the trend difference;
in step S2, for n-29 sub-sequences extracted by 30-year window sliding method, each sub-sequence has 30 continuous years of data, and the 30-year window sliding method extracts the sub-sequence S j Can be expressed as: s j ={X j ,X j+1 ,...,X j+29 0 < j < N-29 and j is equal to N +
For n-29 groups of subsequences extracted by the 30-year window increasing method, the initial sequence is the same as the initial sequence extracted by the 30-year window sliding method, and then each subsequence is increased by one value, so that the subsequence S extracted by the 30-year window increasing method j Can be expressed as: s j ={X 1 ,X 2 ,...,X j+29 0 < j < N-29 and j is equal to N +
In step S5, subsequences S extracted based on 30-year window sliding method and 30-year window growing method respectively j Design of rainstorm values by quantification
Figure FDA0003730645320000011
And annual variation trend of design storm values
Figure FDA0003730645320000012
Two trend values can be obtained and are respectively recorded as
Figure FDA0003730645320000013
The difference between the two is:
Figure FDA0003730645320000014
wherein the content of the first and second substances,
Figure FDA0003730645320000015
represents a design rainstorm value trend value based on a 30-year window sliding method,
Figure FDA0003730645320000016
represents a design rainstorm annual trend value based on a 30-year window growth law,
Figure FDA0003730645320000017
representing the difference of the two trend values;
if it is
Figure FDA0003730645320000018
That is, the design rainstorm value calculated based on the non-consistency assumption is smaller than the design rainstorm value calculated based on the consistency assumption, which indicates that the design rainstorm value calculated by adopting the consistency assumption is safer for engineering construction in the long run; if it is
Figure FDA0003730645320000019
It is shown that designing a rainstorm method using the consistency assumption calculation is disadvantageous in the long run.
2. The method for analyzing the safety of design rainstorm based on the assumed conditions of consistency and non-consistency as claimed in claim 1, wherein in step S2, the assumed conditions of non-consistency and consistency are simulated by using a 30-year window sliding method and a 30-year window increasing method.
3. The method for analyzing the safety of design rainstorm based on assumed conditions of consistency and non-consistency according to claim 1, wherein in step S3, each subsequence of the set S is fitted with the generalized extremum distribution, and parameters of the generalized extremum distribution are estimated by using a maximum likelihood method, and then the design rainstorm value is calculated for a given rainstorm recurrence period, wherein each subsequence obtains the design rainstorm value for the given recurrence period, and two sets of design rainstorm sequences are obtained for each recurrence period, each set containing n to 29 design rainstorm values, because the subsequences are extracted based on a 30-year window sliding method and a 30-year window growing method, respectively.
4. The method for analyzing the safety of the design rainstorm based on the assumed conditions of consistency and inconsistency according to claim 3, wherein the specific process in step S3 is as follows:
the generalized extreme value distribution is based on an extreme value theory, has 3 parameters in total, and is used for quantitatively designing the intensity and the frequency of the rainstorm; the generalized extremum distribution function is as follows:
Figure FDA0003730645320000021
wherein the content of the first and second substances,
Figure FDA0003730645320000022
x represents time, and mu, sigma and k represent a position parameter, a scale parameter and a shape parameter respectively;
calculating by adopting a maximum likelihood estimation method to obtain 3 parameters of the generalized extreme value distribution function, and further finishing the quantification of the designed rainstorm according to the generalized extreme value distribution function; for each subsequence S j Given a design rainstorm in T years, the quantified design rainstorm value is:
Figure FDA0003730645320000023
for a given recurrence period, i.e., T years of design rainstorm, set S ═ S by subsequence 1 ,S 2 ,...,S j ,...,S n-29 And subsequence S j Calculating the design rainstorm value to obtain a design rainstorm value sequence:
Figure FDA0003730645320000024
wherein j is more than 0 and less than N-29 and j is equal to N +
5. The method for analyzing the safety of design rainstorm based on the assumption of consistency and non-consistency as claimed in claim 3, wherein in step S3, the subsequence S extracted based on the 30-year window sliding method and the 30-year window growing method is used as a result of the extraction j In contrast, the subsequence S is given a generalized extremum distribution j When fitting is carried out, fitting parameters obtained by calculation by adopting a maximum likelihood estimation method are notIn the same way, the design rainstorm value of the quantification
Figure FDA0003730645320000031
Differently, two different design rainstorm values can be obtained for each given recurrence period T, and thus two different sets of design rainstorm sequences can be obtained.
6. The method for analyzing the safety of the design storm according to the claim 1, wherein in the step S4, the change trend of the design storm value is calculated by using the linear regression equation:
Figure FDA0003730645320000032
wherein
Figure FDA0003730645320000033
And
Figure FDA0003730645320000034
respectively representing intercept and slope, x is time, y is precipitation, x i Which represents the time of the i-th year,
Figure FDA0003730645320000035
arithmetic mean, y, representing the year sequence i Indicates the precipitation amount of the i-th year,
Figure FDA0003730645320000036
denotes y i The arithmetic mean of (a); the annual variation trend of the design storm values can be determined by designing the storm sequence y T The following formula is used to obtain:
Figure FDA0003730645320000037
wherein
Figure FDA0003730645320000038
Is a sequence y T The jth group of (a) design storm values,
Figure FDA0003730645320000039
is a sequence y T Is calculated as the arithmetic mean of the average of the values,
Figure FDA00037306453200000310
and designing an annual trend value of the rainstorm for the first time of T years.
7. The method for analyzing the safety of the design rainstorm based on the assumed conditions of consistency and non-consistency of the design rainstorm according to claim 6, wherein in the step S4, when the trend analysis is carried out on the design rainstorm sequence, the F distribution is adopted to check the significance of the trend, and the confidence level is 0.05.
CN201810877331.1A 2018-08-03 2018-08-03 Design rainstorm safety analysis method based on assumed conditions of consistency and non-consistency Active CN109271664B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810877331.1A CN109271664B (en) 2018-08-03 2018-08-03 Design rainstorm safety analysis method based on assumed conditions of consistency and non-consistency

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810877331.1A CN109271664B (en) 2018-08-03 2018-08-03 Design rainstorm safety analysis method based on assumed conditions of consistency and non-consistency

Publications (2)

Publication Number Publication Date
CN109271664A CN109271664A (en) 2019-01-25
CN109271664B true CN109271664B (en) 2022-09-27

Family

ID=65153347

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810877331.1A Active CN109271664B (en) 2018-08-03 2018-08-03 Design rainstorm safety analysis method based on assumed conditions of consistency and non-consistency

Country Status (1)

Country Link
CN (1) CN109271664B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103033856A (en) * 2012-12-06 2013-04-10 南京信息工程大学 Rainfall frequency estimation method based on hydrometeorology consistency geographical analysis
CN106777985A (en) * 2016-12-19 2017-05-31 中国水利水电科学研究院 A kind of Hydrologic Series nonuniformity diagnostic method based on classification
CN107797969A (en) * 2017-10-30 2018-03-13 中国水利水电科学研究院 A kind of determination method and device of flood retrieving parameter

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103033856A (en) * 2012-12-06 2013-04-10 南京信息工程大学 Rainfall frequency estimation method based on hydrometeorology consistency geographical analysis
CN106777985A (en) * 2016-12-19 2017-05-31 中国水利水电科学研究院 A kind of Hydrologic Series nonuniformity diagnostic method based on classification
CN107797969A (en) * 2017-10-30 2018-03-13 中国水利水电科学研究院 A kind of determination method and device of flood retrieving parameter

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"Lack of uniform trends but increasing spatial variability in observed Indian rainfall extremes";Subimal Ghosh等;《Nature Climate Change》;20111218;全文 *

Also Published As

Publication number Publication date
CN109271664A (en) 2019-01-25

Similar Documents

Publication Publication Date Title
CN107423546B (en) Multivariable hydrological uncertainty processing method based on Copula function
CN109272146B (en) Flood prediction method based on deep learning model and BP neural network correction
Wang et al. An improved ARIMA model for precipitation simulations
CN107066425B (en) Method for analyzing non-uniformity of ultra-quantitative flood in changing environment
CN107092744B (en) Surface subsidence prediction method based on EMD-SVR
Caillouet et al. Ensemble reconstruction of spatio-temporal extreme low-flow events in France since 1871
CN106529814B (en) Distributed photovoltaic ultra-short term prediction method based on Adaboost clustering and Markov chain
Li et al. A comparison of methods for estimating climate change impact on design rainfall using a high-resolution RCM
CN107423496B (en) Novel random rainfall event generation method
CN110378504B (en) Photovoltaic power generation climbing event probability prediction method based on high-dimensional Copula technology
CN112330065A (en) Runoff forecasting method based on basic flow segmentation and artificial neural network model
KR101463492B1 (en) Error correction method for global climate model using non-stationary quantile mapping
CN113780657A (en) Flood event identification method and device, computer equipment and storage medium
Wang et al. Uncertainty analysis of hydrological processes based on ARMA-GARCH model
CN105653728A (en) Wireless sensor network data collection method based on double-layer model system
CN108694472B (en) Prediction error extreme value analysis method, device, computer equipment and storage medium
CN109271664B (en) Design rainstorm safety analysis method based on assumed conditions of consistency and non-consistency
Jahanbaksh Asl et al. Frequency analysis of climate extreme events in Zanjan, Iran
Gragne et al. Recursively updating the error forecasting scheme of a complementary modelling framework for improved reservoir inflow forecasts
CN109919362B (en) Medium-and-long-term runoff forecasting method considering hydraulic engineering scheduling influence
CN113435630B (en) Basin hydrological forecasting method and system with self-adaptive runoff yield mode
CN110807475B (en) Flood classification, identification and forecast method based on certainty coefficient
CN116432370A (en) Drought prediction method suitable for climate change background
CN110322063B (en) Power consumption simulation prediction method and storage medium
Shahreza et al. Estimating Stationary Inflow Time Series for Planning Water Resources at the Basin Level (Case Study: Karkheh Watershed)

Legal Events

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