CN117522012A - Runoff scene generation method based on seasonal period characteristics - Google Patents
Runoff scene generation method based on seasonal period characteristics Download PDFInfo
- Publication number
- CN117522012A CN117522012A CN202311447184.1A CN202311447184A CN117522012A CN 117522012 A CN117522012 A CN 117522012A CN 202311447184 A CN202311447184 A CN 202311447184A CN 117522012 A CN117522012 A CN 117522012A
- Authority
- CN
- China
- Prior art keywords
- runoff
- season
- scene
- data
- time 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 73
- 230000001932 seasonal effect Effects 0.000 title claims abstract description 27
- 238000005070 sampling Methods 0.000 claims abstract description 44
- 238000005315 distribution function Methods 0.000 claims abstract description 27
- 230000001186 cumulative effect Effects 0.000 claims abstract description 23
- 238000004088 simulation Methods 0.000 claims abstract description 20
- 238000013507 mapping Methods 0.000 claims abstract description 18
- 230000009466 transformation Effects 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims description 28
- 238000005215 recombination Methods 0.000 claims description 6
- 230000006798 recombination Effects 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 7
- 238000013486 operation strategy Methods 0.000 abstract description 5
- 238000011160 research Methods 0.000 abstract description 5
- 238000000605 extraction Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 3
- 238000000546 chi-square test Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000001276 Kolmogorov–Smirnov test Methods 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011158 quantitative evaluation Methods 0.000 description 2
- 238000010998 test method Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000011426 transformation method Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
-
- 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)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Water Supply & Treatment (AREA)
- Game Theory and Decision Science (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Public Health (AREA)
- Quality & Reliability (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Primary Health Care (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Complex Calculations (AREA)
Abstract
A runoff scene generation method based on seasonal characteristics comprises the following steps: step 1, obtaining N k Runoff time sequence data of each season of the year divided into T time periodsStep 2, obtaining the cumulative distribution function cdf of each season j The method comprises the steps of carrying out a first treatment on the surface of the Step 3, expanding the original runoff time sequence data interval to be an extreme interval; step 4, the original runoff time sequence data is processedMapping to extreme intervals to obtainStep 5, obtaining the cumulative distribution function cdf of each period of each season j,t The method comprises the steps of carrying out a first treatment on the surface of the Step 6, selecting LHS sampling times N; step 7, LHS sampling inverse transformation; step 8, quantifying the fluctuation amplitude sigma of each season of the original data j The method comprises the steps of carrying out a first treatment on the surface of the And 9, recombining the sampled data to obtain a random simulation result. The time sequence scene generated by the method is more similar to the original scene in terms of upper and lower limit intervals and scene fluctuation amplitude, can improve the simulation effect of scene simulation and the rationality of further scene extraction, and has important popularization and use values for the method research of formulating the optimal operation strategy of the reservoir or new energy.
Description
Technical Field
The invention relates to the technical field of scene generation, in particular to a runoff scene generation method based on seasonal characteristics.
Background
Scene analysis is needed in the fields of reservoir dispatching optimization operation, grid dispatching system optimization operation and the like, and the optimization typical scene generation technology has important popularization and use values for method research for making the system optimal operation strategy. At present, methods such as Monte Carlo (Monte Carlo simulation) and Latin hypercube (Latin hypercube sampling, LHS) are widely applied in scene generation technology, and compared with random sampling in a conventional Monte Carlo method, LHS sampling is more robust, sampling space covered by the LHS is wider and calculation efficiency is higher under the same sampling scale. However, the LHS sampling method has higher correlation among random variables, and the current methods for data recombination include genetic algorithm, column pair algorithm, gram-Schmidt orthogonalization method, cholesky decomposition method and the like. The method only ranks the data internally aiming at the correlation among the reduced sampling matrixes, and does not consider the time sequence process and fluctuation condition of the original scene.
Therefore, it is necessary to design a runoff scene generating method based on seasonal characteristics to overcome the above-mentioned problems.
Disclosure of Invention
In order to avoid the problems, the runoff scene generation method based on seasonal characteristics is provided, the generated time sequence scene is more similar to an original scene from the upper limit interval, the lower limit interval and the scene fluctuation range, the simulation effect of scene simulation and the rationality of further scene extraction can be improved, and the method has important popularization and use values for the method research of formulating the optimal operation strategy of the reservoir or new energy.
The invention provides a runoff scene generation method based on seasonal characteristics, which comprises the following steps:
step 1, obtaining N k Runoff time sequence data of each season of the year divided into T time periodsJ=1, 2,3,4, which is a season identification, and corresponds to spring, summer, autumn and winter; k=1, 2, [ N ] N k ,N k Is the total number of scenes; t=1, 2, the terms T and T, T is the total time period number;
step 2, for the runoff time sequence dataPerforming distribution fitting to obtain a cumulative distribution function cdf of each season j ;
Step 3, setting a percentile point alpha in each season, and calculating the runoff quantity corresponding to the upper alpha component point and the lower alpha component point by adopting a reverse operation (inverse cumulative distribution function, ICDF) methodAnd will->As a runoff extreme section;
step 4, the original runoff time sequence data is processedMapping to extreme interval, the formula isWherein (1)>Is the time sequence data of the original runoff,the maximum value and the minimum value of the original runoff data; obtaining mapping result, which is new runoff time sequence data
Step 5, for new runoff time sequence dataPerforming distribution fitting to obtain cumulative distribution function cdf of each period of each season j,t ;
Step 6, selecting proper LHS sampling times N, and accumulating a distribution function cdf of each period of each season j,t At [0,1]Equally dividing the probability interval into N sub-intervals, and recording the midpoint of each probability sub-interval as the probability value for inverse transformation sampling(n=1,2,···,N);
Step 7, cumulative distribution function cdf of each period of each season j,t At the position ofPerforming LHS sampling inverse operation to obtain sampling result, which is marked as +.> Obtaining LHS sampling result matrix x of NxT j ;
Step 8, quantifying the fluctuation amplitude sigma of each season of the original data j ;
Step 9, according to the fluctuation range sigma of each season j Sampling the LHS result matrix x j And carrying out random transformation data recombination to obtain a random simulation result reflecting the fluctuation characteristics of the original data.
Preferably, the method for performing the distribution fitting in the step 2 and the step 5 is a kernel density estimation method, a Weibull distribution or a Normal distribution, and when the actual distribution fitting is performed, a hypothesis test method (chi-square test, KS test) is used to compare the fitting effect of the kernel density estimation method (kernel density estimation) and a parameter estimation method (for example, weibull distribution and Normal distribution) on the raw data of runoffs in each season, and a better fitting method is selected.
Preferably, the value interval of the percentile alpha in the step 3 is [95%, 100%).
Preferably, mapping is performed in a normalized manner in step 4.
Preferably, the amplitude sigma of the seasonal fluctuations of the raw data is quantized in step 8 j The method of (1) is as follows: firstly, the original runoff time sequence data is mappedMapping to LHS sample times [1, N ]]Within the interval, get->For mapping result, the formula is +.>Then, calculate the fluctuation of time period before and after each season +.>The formula isAnd calculates the fluctuation of the time period before and after each season +.>Standard deviation sigma of j The fluctuation range of each season is obtained.
Preferably, step 9 comprises the following sub-steps:
9.1 generating compliance mean value 0 and standard deviation sigma for different seasons j Normal random number matrix m of (2) j Matrix m j Is N× (T-1), where m j All the data in (a) are rounded up,
9.2sampling LHS result matrix x j And carrying out random exchange on each row in the row to obtain a random simulation result reflecting the fluctuation characteristic of the original data, wherein the formula is expressed as follows:
wherein,sampling the result matrix x for LHS j Inner nth row and nth column sample data; />For m under j seasons j Random numbers of nth row and nth column in matrix, ">The size change of the number of lines n is limited to the interval [1, N ]]An inner part; />And (5) simulating the extreme runoff scene in the j seasons.
Compared with the prior art, the invention has the following beneficial effects: the method is mainly applied to rationally describing possible uncertainty by using time sequence runoff data or other data with seasonal characteristics, and a scene generation method based on Latin hypercube (Latin hypercube sampling, LHS) probability distribution inverse transformation of seasonal division is adopted for time sequence data to generate a time sequence scene under extreme conditions. And obtaining a cumulative distribution function by adopting two different methods of parameter estimation, hypothesis test and non-parameter estimation (nuclear density estimation), and generating a time sequence scene by an inverse transformation method. The innovation provides that the fluctuation quantitative evaluation is carried out on different seasonal time sequence scenes and the method is used for simulating seasonal fluctuation characteristics during data recombination. The scene generation method for simulating the seasonal fluctuation characteristics can enable the generated time sequence scene to be more similar to the original scene from the upper limit interval, the lower limit interval and the scene fluctuation range, and can improve the simulation effect of scene simulation and the rationality of further scene extraction. Has important popularization and use values for the research of the method for making the optimal operation strategy of the reservoir or the new energy.
Drawings
FIG. 1 is a flow chart of a method for generating a runoff scenario based on seasonal characteristics according to a preferred embodiment of the present invention;
FIG. 2 is a diagram of historical runoff timing data according to a preferred embodiment of the present invention;
FIG. 3 is a graph showing the fitting result of the historical runoff amount distribution in each season according to a preferred embodiment of the present invention;
FIG. 4 is a diagram illustrating the result of the historical traffic map according to a preferred embodiment of the present invention;
FIG. 5 is a schematic diagram of the inverse LHS sample transform result in accordance with a preferred embodiment of the present invention;
FIG. 6 is a schematic diagram of a wave random simulation process according to a preferred embodiment of the present invention;
FIG. 7 is a schematic diagram of the result of random simulation of a runoff extreme scene according to a preferred embodiment of the present invention.
Detailed Description
As shown in fig. 1, the runoff scene generating method based on seasonal characteristics provided in this embodiment includes the following steps:
step 1, obtaining N k Runoff time sequence data of each season of the year divided into T time periodsJ=1, 2,3,4, which is a season identification, and corresponds to spring, summer, autumn and winter; k=1, 2, [ N ] N k ,N k Is the total number of scenes; t=1, 2, the terms T and T, T is the total time period number;
step 2, for the runoff time sequence dataPerforming distribution fitting to obtain a cumulative distribution function cdf of each season j ;
Step 3, setting a percentile point alpha in each season, and calculating the runoff quantity corresponding to the upper alpha component point and the lower alpha component point by adopting a reverse operation (inverse cumulative distribution function, ICDF) methodAnd will->As a runoff extreme section;
step 4, the original runoff time sequence data is processedNormalized mapping to extreme intervals, the formula isWherein (1)>Is the time sequence data of the original runoff,the maximum value and the minimum value of the original runoff data; obtaining mapping result, which is new runoff time sequence data
Step 5, for new runoff time sequence dataPerforming distribution fitting to obtain cumulative distribution function cdf of each period of each season j,t ;
Step 6, selecting proper LHS sampling times N, and accumulating a distribution function cdf of each period of each season j,t At [0,1]Equally dividing the probability interval into N sub-intervals, and recording the midpoint of each probability sub-interval as the probability value for inverse transformation sampling(n=1,2,···,N);
Step 7, cumulative distribution function cdf of each period of each season j,t At the position ofPerforming LHS sampling inverse operation to obtain sampling result, which is marked as +.> Obtaining LHS sampling result matrix x of NxT j ;
Step 8, quantifying the fluctuation amplitude sigma of each season of the original data j The method comprises the steps of carrying out a first treatment on the surface of the Firstly, the original runoff time sequence data is mappedMapping to LHS sample times [1, N ]]Within the interval, get->For the mapping result, the formula isThen, calculate the fluctuation of time period before and after each season +.>The formula isAnd calculates the fluctuation of the time period before and after each season +.>Standard deviation sigma of j The fluctuation range of each season is obtained;
step 9, according to the fluctuation range sigma of each season j Sampling the LHS result matrix x j Carrying out random transformation data recombination to obtain a random simulation result reflecting the fluctuation characteristics of the original data; comprises the following stepsThe sub-steps are as follows:
9.1 generating compliance mean value 0 and standard deviation sigma for different seasons j Normal random number matrix m of (2) j Matrix m j Is N× (T-1), where m j All data in (2) are rounded;
9.2 sampling result matrix x for LHS j And carrying out random exchange on each row in the row to obtain a random simulation result reflecting the fluctuation characteristic of the original data, wherein the formula is expressed as follows:
wherein,sampling the result matrix x for LHS j Inner nth row and nth column sample data; />For m under j seasons j Random numbers of nth row and nth column in matrix, ">The size change of the number of lines n is limited to the interval [1, N ]]An inner part; />And (5) simulating the extreme runoff scene in the j seasons.
Preferably, the method for performing the distribution fitting in the step 2 and the step 5 is a kernel density estimation method, a Weibull distribution or a Normal distribution, and when the actual distribution fitting is performed, a hypothesis test method (chi-square test, KS test) is used to compare the fitting effect of the kernel density estimation method (kernel density estimation) and a parameter estimation method (for example, weibull distribution and Normal distribution) on the raw data of runoffs in each season, and a better fitting method is selected.
The runoff scene based on seasonal period characteristics is generated by using the time sequence data of a certain hydrologic station and a certain 20-year runoff amount ten-day scale, and the method comprises the following steps:
step 1, season division. In the present exampleUsing time sequence data of a certain hydrologic station in a certain 20-year runoff amount in ten-day scale, dividing the original data according to time stamps corresponding to seasons to obtain a value of a kth scene in a jth season at a t moment, and marking the value as(j=1, 2,3,4, season identification, spring, summer, autumn, winter; k=1, 2, & N k The method comprises the steps of carrying out a first treatment on the surface of the t=1, 2, T; (wherein T and N k Respectively representing the number of time periods and the number of scenes, in this example Tfetch 9,N k 20, i.e. 20 pieces of runoff time sequence scene of 9 ten days in each season). The historical runoff data are shown in the following table 1, the flow units (m/s), and the historical runoff time sequence diagram is shown in fig. 1.
Table 1 20 year runoff ten-day scale timing data
And 2, comparing the fitting method. Performing distribution fitting on historical runoff data by using Kernel density estimation (Kernel) method, weibull distribution and Normal distribution, wherein the function expression of Kernel density estimation is as follows Representing an estimate of the cumulative distribution function F (x)And h represents window width, n represents sample capacity, and K (·) represents kernel function. There are many kinds of kernel functions, which in this example select Gaussian kernel with the expression +.>The fitting results are shown in fig. 3. Chi-square statistics are core indicators of chi-square tests, and the smaller the chi-square statistics are, the smaller the deviation between sample data and theoretical distribution is. The chi-square statistic result pair is shown in table 2.
TABLE 2 chi-square statistics comparison
As can be seen from the chi-square statistic comparison, the fitting effect of using the kernel density estimation method in this example is better than the Weibull distribution and the Normal distribution, so the kernel density estimation method is selected as the preferred fitting method.
Step 3, acquiring a cumulative distribution function cdf of each season j Carrying out distribution fitting on the raw runoff data of each season by using a Kernel fitting method to obtain a cumulative distribution function cdf of each season j The results are shown in Table 3.
TABLE 3 estimation of core Density to obtain the cumulative distribution function estimate cdf for each season j
And 4, expanding an original data interval. Based on the obtained cumulative distribution function cdf j In this embodiment, the dividing point α=99.9% is selected, and the new upper and lower limit data of each season is calculated by performing the inverse operation (ICDF), and recorded asWherein (1)>Corresponding to the lower alpha positionThe point diameter flow; />Corresponds to the alpha split site diameter flow rate and willThe extreme interval is marked; extreme runoff interval +.>Interval with original runoffIs shown in table 4.
TABLE 4 comparison of extreme runoff intervals with original runoff intervals
By comparing extreme runoff intervalsInterval with original runoffIt is understood that, in the present embodiment, when the quantile α=99.9% is selected, the runoff data interval is enlarged by about 18% on average in each season after the original data interval enlarging method is used.
And 5, mapping the original runoff data to an extreme interval. Expanding the original runoff data interval to an extreme interval in a mapping mode by a normalization mapping mode, wherein the formula is that Wherein (1)>For the original runoff time sequence data, < > a->The maximum value and the minimum value of the original runoff data; obtaining mapping result, which is new radial flow time sequence data-> The schematic diagram is shown in fig. 4.
Step 6, obtaining the cumulative distribution function cdf of each period j,t Radial flow data mapped to extreme intervals using a nuclear density estimation methodCarrying out distribution fitting on the data of each period of each season to obtain a cumulative distribution function cdf of each period of each season j,t 。
And 7, determining the LHS sampling times N. In this embodiment, the sampling number n=100 is selected, and then each period of each season will accumulate [0,1 ] of the distribution function according to the LHS sampling number N]The probability interval is divided into 100 sub-intervals, and the midpoint of each probability sub-interval is recorded as the probability value for inverse transformation sampling(n=1,2,···,N)。
Step 8, inverse LHS sampling transformation. The cumulative distribution function of each period of each season is respectively in the probability valuePerforming inverse operation to obtain sampling result, which is marked as +.> Obtaining LHS sampling result matrix x of NxT j The method comprises the steps of carrying out a first treatment on the surface of the The sampling result is shown in fig. 5.
Step 9, quantifying the fluctuation amplitude sigma of each season of the original data j . First, original data is mapped to LHS sampling times [1, N ] by means of normalized mapping]Within the interval, the formula isWherein->Is normalized result. Calculating fluctuation of time period before and after each season>As the calculation basis of the amplitude of the analog data fluctuation, the formula is that Then calculate every season +.>Standard deviation sigma of j And as the magnitude of its fluctuation; every season in this embodiment +.>Standard deviation sigma of j The calculation results are shown in Table 5.
TABLE 5 seasons of the yearStandard deviation sigma of j
Step 10, reorganizing the sampled data. According to the quantized fluctuation amplitude result sigma of each season j Generating a compliance mean value of 0 and a standard deviation of sigma for different seasons j Normal random number matrix m of (2) j Matrix m j Is N× (T-1), m j All data in (1) are rounded and then the result matrix x is sampled by LHS j (matrix x j Is of a rank size of N x T) to obtain a random simulation result reflecting the fluctuation characteristics of the original data by performing random exchange of each row in the rank, the formula can be expressed as,wherein (1)>Sampling the result matrix x for LHS j Inner nth row and nth column sample data; />For m under the season j Random numbers of nth row and nth column in matrix, ">The size change of the number of lines n is limited to the interval [1, N ]]An inner part; />And (5) simulating the extreme runoff scene in the j seasons. The wave random simulation process is shown in fig. 6, the extreme runoff scene simulation result is shown in fig. 7, and the whole flow of the embodiment of the invention is shown in fig. 1.
The method is mainly applied to rationally describing possible uncertainty by using time sequence runoff data or other data with seasonal characteristics, and a scene generation method based on Latin hypercube (Latin hypercube sampling, LHS) probability distribution inverse transformation of seasonal division is adopted for time sequence data to generate a time sequence scene under extreme or non-extreme conditions. And obtaining a cumulative distribution function by adopting two different methods of parameter estimation, hypothesis test and non-parameter estimation (nuclear density estimation), and generating a time sequence scene by an inverse transformation method. The innovation provides that the fluctuation quantitative evaluation is carried out on different seasonal time sequence scenes and the method is used for simulating seasonal fluctuation characteristics during data recombination. The scene generation method for simulating the seasonal fluctuation characteristics can enable the generated time sequence scene to be more similar to the original scene from the upper limit interval, the lower limit interval and the scene fluctuation range, and can improve the simulation effect of scene simulation and the rationality of further scene extraction. Has important popularization and use values for the research of the method for making the optimal operation strategy of the reservoir or the new energy.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solution of the present invention, and not for limiting the same; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some or all of the technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit of the invention.
Claims (6)
1. The runoff scene generation method based on the seasonal period characteristics is characterized by comprising the following steps of:
step 1, obtaining N k Runoff time sequence data of each season of the year divided into T time periodsJ=1, 2,3,4, which is a season identification, and corresponds to spring, summer, autumn and winter; k=1, 2, [ N ] N k ,N k Is the total number of scenes; t=1, 2, the terms T and T, T is the total time period number;
step 2, for the runoff time sequence dataPerforming distribution fitting to obtain a cumulative distribution function cdf of each season j ;
Step 3, setting a percentile alpha in each season,calculating the runoff corresponding to the upper alpha split point and the lower alpha split point by adopting an inverse operation methodAnd will->As a runoff extreme section;
step 4, the original runoff time sequence data is processedMapping to the extreme interval to obtain new runoff time sequence data
Step 5, for new runoff time sequence dataPerforming distribution fitting to obtain cumulative distribution function cdf of each period of each season j,t ;
Step 6, selecting proper LHS sampling times N, and accumulating a distribution function cdf of each period of each season j,t At [0,1]Equally dividing the probability interval into N sub-intervals, and recording the midpoint of each probability sub-interval as the probability value for inverse transformation sampling
Step 7, cumulative distribution function cdf of each period of each season j,t At the position ofPerforming LHS sampling inverse operation to obtain sampling result, which is marked as +.>Obtaining LHS sampling result matrix x of NxT j ;
Step 8, quantifying the fluctuation amplitude sigma of each season of the original data j ;
Step 9, according to the fluctuation range sigma of each season j Sampling the LHS result matrix x j And carrying out random transformation data recombination to obtain a random simulation result reflecting the fluctuation characteristics of the original data.
2. A runoff scene generating method based on seasonal characteristics as defined in claim 1, wherein: the method for performing distribution fitting in the step 2 and the step 5 is a nuclear density estimation method, a Weibull distribution or a Normal distribution.
3. A runoff scene generating method based on seasonal characteristics as defined in claim 1, wherein: and 3, the value interval of the percentile alpha in the step 3 is [95%, 100%).
4. A runoff scene generating method based on seasonal characteristics as defined in claim 1, wherein: and 4, mapping in a normalization mode.
5. A runoff scene generating method based on seasonal characteristics as defined in claim 1, wherein: step 8, quantifying the fluctuation amplitude sigma of each season of the original data j The method of (1) is as follows: firstly, the original runoff time sequence data is mappedMapping to LHS sample times [1, N ]]Within the interval, get->Then, calculate the fluctuation of time period before and after each season +.>The formula is->And calculates the fluctuation of the time period before and after each season +.>Standard deviation sigma of j The fluctuation range of each season is obtained.
6. A runoff scene generating method based on seasonal characteristics as defined in claim 1, wherein: step 9 comprises the following sub-steps:
9.1 generating compliance mean value 0 and standard deviation sigma for different seasons j Normal random number matrix m of (2) j Matrix m j Is N× (T-1), where m j All data in (2) are rounded;
9.2 sampling result matrix x for LHS j And carrying out random exchange on each row in the row to obtain a random simulation result reflecting the fluctuation characteristic of the original data, wherein the formula is expressed as follows:
wherein,sampling the result matrix x for LHS j Inner nth row and nth column sample data; />For m under j seasons j Random numbers of nth row and nth column in matrix, ">The size change of the number of lines n is limited to the interval [1, N ]]An inner part; />And (5) simulating the extreme runoff scene in the j seasons.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311447184.1A CN117522012A (en) | 2023-11-02 | 2023-11-02 | Runoff scene generation method based on seasonal period characteristics |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311447184.1A CN117522012A (en) | 2023-11-02 | 2023-11-02 | Runoff scene generation method based on seasonal period characteristics |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117522012A true CN117522012A (en) | 2024-02-06 |
Family
ID=89763691
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311447184.1A Pending CN117522012A (en) | 2023-11-02 | 2023-11-02 | Runoff scene generation method based on seasonal period characteristics |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117522012A (en) |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100299126A1 (en) * | 2009-04-27 | 2010-11-25 | Schlumberger Technology Corporation | Method for uncertainty quantifiation in the performance and risk assessment of a carbon dioxide storage site |
KR101486798B1 (en) * | 2013-12-13 | 2015-02-04 | 세종대학교산학협력단 | Step-wise scaling method for correcting bias of climate information |
CN107294087A (en) * | 2017-06-23 | 2017-10-24 | 清华大学 | A kind of integrated energy system typical scene set creation method containing meteorological energy sources |
CN110276104A (en) * | 2019-05-23 | 2019-09-24 | 武汉大学 | A kind of seasonal design flood calculation method under set climatic model |
CN110311420A (en) * | 2019-06-29 | 2019-10-08 | 南京理工大学 | A kind of generation method of scene joint power output timing scene |
CN112116149A (en) * | 2020-09-17 | 2020-12-22 | 河海大学 | Multi-station medium and long term runoff rolling probability prediction method considering forecast uncertainty associated evolution characteristics |
CN114357865A (en) * | 2021-12-17 | 2022-04-15 | 重庆大唐国际彭水水电开发有限公司 | Hydropower station runoff and associated source load power year scene simulation and prediction method thereof |
CN114757548A (en) * | 2022-04-22 | 2022-07-15 | 国网福建省电力有限公司电力科学研究院 | Wind power energy storage equipment adjusting performance evaluation method adopting scene construction |
CN115081693A (en) * | 2022-06-08 | 2022-09-20 | 天津大学 | Two-stage capacity configuration method of water-wind-light complementary system with uncertain energy |
CN116579120A (en) * | 2023-03-09 | 2023-08-11 | 中国南方电网有限责任公司 | Synthetic method-based solar-wind power output sequence scene generation method |
CN116720448A (en) * | 2023-08-09 | 2023-09-08 | 长江三峡集团实业发展(北京)有限公司 | Wind power generation random simulation method, device, equipment and medium |
CN116821585A (en) * | 2023-06-28 | 2023-09-29 | 中国地质大学(北京) | Non-parameter time downscaling method and system for long-sequence reconstruction runoff |
-
2023
- 2023-11-02 CN CN202311447184.1A patent/CN117522012A/en active Pending
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100299126A1 (en) * | 2009-04-27 | 2010-11-25 | Schlumberger Technology Corporation | Method for uncertainty quantifiation in the performance and risk assessment of a carbon dioxide storage site |
KR101486798B1 (en) * | 2013-12-13 | 2015-02-04 | 세종대학교산학협력단 | Step-wise scaling method for correcting bias of climate information |
CN107294087A (en) * | 2017-06-23 | 2017-10-24 | 清华大学 | A kind of integrated energy system typical scene set creation method containing meteorological energy sources |
CN110276104A (en) * | 2019-05-23 | 2019-09-24 | 武汉大学 | A kind of seasonal design flood calculation method under set climatic model |
CN110311420A (en) * | 2019-06-29 | 2019-10-08 | 南京理工大学 | A kind of generation method of scene joint power output timing scene |
CN112116149A (en) * | 2020-09-17 | 2020-12-22 | 河海大学 | Multi-station medium and long term runoff rolling probability prediction method considering forecast uncertainty associated evolution characteristics |
CN114357865A (en) * | 2021-12-17 | 2022-04-15 | 重庆大唐国际彭水水电开发有限公司 | Hydropower station runoff and associated source load power year scene simulation and prediction method thereof |
CN114757548A (en) * | 2022-04-22 | 2022-07-15 | 国网福建省电力有限公司电力科学研究院 | Wind power energy storage equipment adjusting performance evaluation method adopting scene construction |
CN115081693A (en) * | 2022-06-08 | 2022-09-20 | 天津大学 | Two-stage capacity configuration method of water-wind-light complementary system with uncertain energy |
CN116579120A (en) * | 2023-03-09 | 2023-08-11 | 中国南方电网有限责任公司 | Synthetic method-based solar-wind power output sequence scene generation method |
CN116821585A (en) * | 2023-06-28 | 2023-09-29 | 中国地质大学(北京) | Non-parameter time downscaling method and system for long-sequence reconstruction runoff |
CN116720448A (en) * | 2023-08-09 | 2023-09-08 | 长江三峡集团实业发展(北京)有限公司 | Wind power generation random simulation method, device, equipment and medium |
Non-Patent Citations (6)
Title |
---|
IL-WON JUNG等: "Uncertainty assessment of climate change impacts for hydrologically distinct river basins", JOURNAL OF HYDROLOGY, vol. 466, 15 August 2012 (2012-08-15), pages 73 - 87, XP028941751, DOI: 10.1016/j.jhydrol.2012.08.002 * |
代倩;曾平良;周勤勇;赵峰;李柏青;: "多风电场与梯级水电站协调运行对电力系统可靠性的影响", 电网技术, no. 06, 5 June 2015 (2015-06-05), pages 1679 - 1684 * |
周研来;梅亚东;张代青;: "一种新的径流过程随机模拟方法", 水利水电科技进展, no. 03, 20 June 2011 (2011-06-20), pages 9 - 12 * |
李炳锋等: "金沙江流域实际蒸散发GRACE重力卫星遥感重构不确定性分析", 水资源保护, vol. 39, no. 4, 20 July 2023 (2023-07-20), pages 159 - 166 * |
李驰;刘纯;黄越辉;王伟胜;: "基于波动特性的风电出力时间序列建模方法研究", 电网技术, no. 01, 5 January 2015 (2015-01-05), pages 208 - 214 * |
王玲玲;王昕;郑益慧;孙洪波;李佳春;赵明阳;: "计及多个风电机组出力相关性的配电网无功优化", 电网技术, no. 11, 6 July 2017 (2017-07-06), pages 3463 - 3469 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111091139B (en) | Photovoltaic prediction method, device and equipment for similar day clustering and readable storage medium | |
CN110674864B (en) | Wind power abnormal data identification method comprising synchronous phasor measurement device | |
CN114792156A (en) | Photovoltaic output power prediction method and system based on curve characteristic index clustering | |
CN116128141A (en) | Storm surge prediction method and device, storage medium and electronic equipment | |
CN112651543A (en) | Daily electric quantity prediction method based on VMD decomposition and LSTM network | |
Li et al. | GMM-HMM-based medium-and long-term multi-wind farm correlated power output time series generation method | |
CN114819374A (en) | Regional new energy ultra-short term power prediction method and system | |
CN114357737B (en) | Agent optimization calibration method for time-varying parameters of large-scale hydrologic model | |
CN106909798A (en) | A kind of Daily rainfall multi-mode collection approach based on cumulative probability curve | |
CN116739172B (en) | Method and device for ultra-short-term prediction of offshore wind power based on climbing identification | |
CN117522012A (en) | Runoff scene generation method based on seasonal period characteristics | |
CN110414776B (en) | Quick response analysis system for power utilization characteristics of different industries | |
CN116845875A (en) | WOA-BP-based short-term photovoltaic output prediction method and device | |
CN116307139A (en) | Wind power ultra-short-term prediction method for optimizing and improving extreme learning machine | |
CN114676931B (en) | Electric quantity prediction system based on data center technology | |
CN115456286A (en) | Short-term photovoltaic power prediction method | |
CN112581311B (en) | Method and system for predicting long-term output fluctuation characteristics of aggregated multiple wind power plants | |
CN111861259A (en) | Load modeling method, system and storage medium considering time sequence | |
CN116627953B (en) | Method for repairing loss of groundwater level monitoring data | |
CN117909888B (en) | Intelligent artificial intelligence climate prediction method | |
CN110570322B (en) | Agricultural meteorological index insurance rate calibrating method based on time sequence simulation | |
CN114282431B (en) | Runoff interval prediction method and system based on improved SCA and QRGRU | |
CN113656953B (en) | Wind power sequence modeling method based on state number optimal decision model | |
Kumar | Solar Radiation Analysis for Predicting Climate Change Using Deep Learning Techniques | |
CN117609791A (en) | Method for generating medium-long time scale sequence by considering wind-light-load correlation |
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 |