CN118157104A - Comprehensive performance evaluation method for power grid rolling scheduling strategy - Google Patents
Comprehensive performance evaluation method for power grid rolling scheduling strategy Download PDFInfo
- Publication number
- CN118157104A CN118157104A CN202410193952.3A CN202410193952A CN118157104A CN 118157104 A CN118157104 A CN 118157104A CN 202410193952 A CN202410193952 A CN 202410193952A CN 118157104 A CN118157104 A CN 118157104A
- Authority
- CN
- China
- Prior art keywords
- evaluation
- weight
- calculating
- index
- indexes
- 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
- 238000011156 evaluation Methods 0.000 title claims abstract description 67
- 238000005096 rolling process Methods 0.000 title claims abstract description 17
- 238000000034 method Methods 0.000 claims abstract description 38
- 239000011159 matrix material Substances 0.000 claims description 34
- 238000004364 calculation method Methods 0.000 claims description 14
- 230000000295 complement effect Effects 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 11
- 238000010606 normalization Methods 0.000 claims description 10
- 238000005265 energy consumption Methods 0.000 claims description 6
- 238000013210 evaluation model Methods 0.000 claims description 6
- 230000003749 cleanliness Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000011158 quantitative evaluation Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/007—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources
- H02J3/0075—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources for providing alternative feeding paths between load and source according to economic or energy efficiency considerations, e.g. economic dispatch
-
- 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
- G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/003—Load forecast, e.g. methods or systems for forecasting future load demand
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/004—Generation forecast, e.g. methods or systems for forecasting future energy generation
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Power Engineering (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Public Health (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to the technical field of generalization evaluation of power grid active cruising scheduling, in particular to a comprehensive performance evaluation method for a power grid rolling scheduling strategy. The method carries out active cruising scheduling on a scheduling strategy under different source load scenes by a parallel computing and evaluating method, and computes each sub-index of a scheduling result; and finally, calculating the weight of each index by using an AHP-inverse entropy weight model, and comprehensively grading the scheduling results in different scenes by using a fuzzy evaluation method.
Description
Technical Field
The invention relates to the technical field of generalization evaluation of active cruising scheduling of a power grid, in particular to a comprehensive performance evaluation method for a rolling scheduling strategy of the power grid, which is used for evaluating decision performance of the scheduling strategy in various complex source load scenes and partial fault scenes.
Background
Along with the gradual promotion of the construction of novel electric power systems in China, the number of main bodies involved in the optimized dispatching of the electric power systems is increased, the fluctuation of power supplies is enhanced, and the difficulty of the prospective dispatching of the power grid is greatly improved. To address this challenge, various data-driven based or partially data-driven intelligent algorithms are applied to grid look-ahead scheduling aid decisions to promote the automation and intelligence level of grid look-ahead scheduling in multiple operating scenarios. Under the background, aiming at the characteristic of strong adaptability of the prospective scheduling intelligent agent to multiple scenes, the decision effect of the intelligent agent is comprehensively evaluated, so that the construction of the prospective scheduling intelligent algorithm is guided, and the method has practical significance.
In the field of power grid economic dispatch evaluation, the existing research mainly uses the economic dispatch agent design as a link of the economic dispatch agent design, and is used for comparing with the traditional dispatch method to verify the effectiveness of the agent design, the whole flow is simpler, and the research specially aiming at the economic dispatch evaluation is less. The system is developed for daily power generation dispatching planning, a power grid economic dispatching pre-evaluation system is provided, a basic data acquisition subsystem, a standardized reference building subsystem and a deviation rate calculation subsystem are constructed, and the automation level of power grid economic dispatching pre-evaluation is improved. The method is researched and oriented to various data driving methods, and the intelligent algorithm for power grid optimization scheduling is evaluated from the dimensions of decision optimality, computational complexity and the like, so that the computational efficiency advantage of the data driving method is fully embodied. The influence of bad data in the power grid look-ahead scheduling is researched, and a sensitivity matrix of the look-ahead scheduling on the bad data is provided so as to rapidly evaluate the influence of data loss or errors on the power grid look-ahead scheduling decision when a line is blocked. However, the research is not performed on the quality of the scheduling policy decision under a large number of different source load and fault scenes, so that the scheduling policy is evaluated according to the performance of a plurality of source load scenes with different characteristics.
Disclosure of Invention
Therefore, it is necessary to provide a parallel computing and evaluating method for performing active cruising scheduling on a scheduling strategy under different source load scenes, and computing each sub-index of a scheduling result; and finally, calculating the weight of each index by using an AHP-inverse entropy weight model, and comprehensively grading the scheduling results in different scenes by using a fuzzy evaluation method.
In order to achieve the above purpose, the present invention provides the following technical solutions:
a comprehensive performance evaluation method for a power grid rolling scheduling strategy comprises the following steps:
S1, selecting a plurality of new energy day-ahead output curves and load prediction curves, multiplying the curves by Gaussian distribution sampling curves with different mean values and variances, and generating a plurality of groups of source load curve scenes with different deviation degrees; the scheduling strategy is subjected to parallel decision calculation under the scenes, and a daily rolling scheduling result obtained by scheduling calculation is obtained;
S2, constructing generalized sub-indexes comprising safety, economy and cleanliness, calculating sub-indexes under different scenes, and finally carrying out normalization treatment and extremely-large treatment on the indexes;
s3, establishing a fuzzy judgment matrix F by adopting an improved fuzzy AHP analytic hierarchy process, and calculating the weight of each index by using the AHP method to obtain subjective weight
S4, calculating the weight of each index by adopting an inverse entropy weight model to obtain objective weight
S5, subjective weight calculated according to the steps S3 and S4And objective weight/>Calculating a variable weight coefficient omega i of the equalization function, and forming an index weight matrix omega;
S6, calculating membership functions of the generalization indexes by using the normalized result in the step S2;
and S7, carrying out quantitative evaluation on the scheduling results of different source load scenes according to membership functions of the generalization indexes, and obtaining scheduling performance results under different source load scenes.
According to further optimization of the technical scheme, the step S1 specifically comprises the following steps:
S1.1, taking a group of new energy output curves PG renewable, wherein PG i is the output value of the new energy at the ith moment, and taking one-day output curves at 15 minutes intervals, so i takes 96 as shown in the following formula:
PGrenewable=[pg1,pg2,...,pgi],i=96
Taking a set of load prediction curves P L, wherein pl i is a predicted value of the load at the ith moment, and the formula is as follows:
PL=[pl1,pl2,...,pli],i=96
And generating N sets of bias curves a= { a 1,a2,...,ai,...aN}T by using a gaussian distribution model, wherein a i is a set of matrices which contain 96 bias values and conform to gaussian distribution, N j is a bias value, and the formula is as follows:
ai=[n1,n2,...nj],nj~Normal(μj,σj 2),j∈[1,96]
Finally, generating N groups of source load scene curves PG renewable,PL with different deviation degrees, wherein:
PL′=AN×96·PL96×1;PGrenewable′=AN×96·PGrenewable,96×1
S1.2, respectively loading the source load curves into different threads, and then loading scheduling strategies to be evaluated into the threads in parallel; and finally, carrying out 96-point parallel calculation to obtain an active cruising path result.
The technical scheme is further optimized, and the step S2 specifically refers to: the calculation method for giving the index by comprehensively considering various factors such as system cost, section power flow out-of-limit condition, load balance, new energy consumption rate, balance machine out-of-limit power and the like comprises the following steps:
S2.1, system cost:
wherein P i is the decision output of the corresponding unit i at a certain moment, C i is the unit running cost of the corresponding unit, N units are in total in the system, C max is a normalization constant, the lower the system cost is, the higher the score is, and the upper limit is 100 score;
S2.2, section tide out-of-limit:
Wherein S is the total section number, r s is the rewarding value of the S-th section, and the calculation method of r s is as follows:
Wherein P s is the tide of section s, P s max is the upper limit of section s, and P s min is the lower limit of section s;
s2.3, load balancing:
Wherein L represents the total system real load at the moment t, P i is the output of the unit i, N is the total number of the units, and the reward is 0 when the load is balanced;
s2.4, new energy consumption rate:
Wherein, P i represents the output of the current new energy, P i max represents the predicted maximum output of the current new energy, and M new energy units are all arranged.
The technical scheme is further optimized, and the step S3 specifically comprises the following steps:
S3.1, constructing a fuzzy complementary judgment matrix F, namely constructing F m×n by a 1-9 scale method, wherein the definition of the 1-9 scale method is shown in the following table:
s3.2, calculating a weight formula by using the fuzzy complementary matrix, wherein the weight formula is as follows:
Wherein F ij represents an element in the fuzzy complementary determining matrix F m×n, and W i represents a weight of the i-th index;
s3.3, carrying out consistency test through fuzzy complementary judgment matrix F:
the matrix F is first summed row by row as follows:
Then transformed by the result of the above equation: Finally obtaining a fuzzy consistency matrix R
Wherein, the relation of each element in R needs to satisfy:
0≤rij≤1,rii=0.5,rij+rji=1
And is also provided with There is r ij=rik-rjk +0.5.
The technical scheme is further optimized, and the step S4 specifically comprises the following steps:
S4.1, forward processing of evaluation indexes: converting the minimum index, the intermediate index and the interval index into the maximum index; the formula is as follows:
Very small- > very large: x i=max{x1,x2,...,xi}-xi, if all elements are positive, the reciprocal can be directly taken
Intermediate- > very large: m=max { |x i-xbest | },
M=max{a-min{xi},max{xi}-b}
Regional- > very large:
S4.2, a normalization process, in which n evaluation objects, m indexes, i.e., x=x n×m, are normalized as follows:
s4.3, calculating the proportion of each element, wherein the sum of the proportion is 1, namely normalizing the standardized matrix, and calculating a probability matrix p, wherein the formula is as follows:
S4.4, calculating the information entropy e of each index
S4.5, calculating the weight, wherein 1-e is the useful value of the information, normalizing to obtain the weight,
The technical scheme is further optimized, and the step S5 specifically comprises the following steps:
S5.1, calculating the fixed weight
S5.2, calculating an evaluation value v i:
Wherein, The constant weight coefficient of the jth key index in the ith item is given, and x ij is the evaluation value of the jth key index in the ith item;
S5.3, variable weight coefficient
S5.4 equalization function
S5.5, the variable weight coefficient of the equalization function is as follows:
Wherein a is an equalization factor.
The technical scheme is further optimized, and the step S6 specifically comprises the following steps:
Calculating membership functions of all indexes by using the result x i obtained by normalizing and summing the step S2 to be the maximum, wherein the fuzzy relation of all indexes and state evaluation indexes is reflected by adopting a normal membership function A (x i) according to the characteristic decision of the active cruise generalization evaluation indexes of the scheduling strategy
Where a represents the boundary of the membership function and σ represents the standard deviation of the normal curve.
The technical scheme is further optimized, and the step S7 specifically comprises the following steps:
S7.1, calculating a single factor evaluation matrix R through the membership function in the step S6:
Wherein m is the number of generalization indexes;
S7.2, establishing a factor set of comprehensive evaluation: u= (U 1,u2,...,um) where element U 1 represents the ith factor affecting the evaluation object, an evaluation set of comprehensive evaluation is established: v= (V 1,v2,...,vm) where element V 1 represents the j-th evaluation result;
s7.3, building a comprehensive evaluation model:
B=ΩR=(ω1,ω2,…ωm)R=(b1,b2,b3,b4)
Wherein the method comprises the steps of
And S7.4, according to the maximum membership degree principle and the evaluation set V, obtaining the evaluation results of the scheduling strategy under different source load scenes through the evaluation model B.
Compared with the prior art, the technical scheme has the beneficial effects that: firstly, constructing a multi-operation source load scene set based on actual data, and carrying out multi-scene active scheduling on a scheduling strategy in parallel; secondly, the invention adopts a subjective and objective comprehensive weighting method, which is more scientific compared with other weighting design methods; third, the improved fuzzy comprehensive evaluation method of the invention uses a normal membership function to reflect the fuzzy relation between each sub-index and the state evaluation.
Drawings
FIG. 1 is a flow chart of parallel decision making of a scheduling strategy in multiple scenarios;
FIG. 2 is a graph of a scheduling policy generalization evaluation index system;
FIG. 3 is a flow chart of a scheduling policy state evaluation method.
Detailed Description
In order to describe the technical content, constructional features, achieved objects and effects of the technical solution in detail, the following description is made in connection with the specific embodiments in conjunction with the accompanying drawings.
As shown in fig. 1, a comprehensive performance evaluation method for a rolling scheduling strategy of a power grid includes the following sequential steps:
S1, selecting a plurality of new energy day-ahead output curves and load prediction curves, multiplying the curves by Gaussian distribution sampling curves with different mean values and variances, and generating a plurality of groups of source load curve scenes with different deviation degrees; and the scheduling strategies to be evaluated are calculated in parallel under the scenes, and a calculated daily rolling scheduling result is obtained.
The step S1 specifically comprises the following steps:
S1.1, a group of new energy output curves PG renewable are taken, wherein PG i is the output value of the new energy at the ith moment, and the output curves of one day are taken at intervals of 15 minutes. So i takes 96, as shown in the following formula:
PGrenewable=[pg1,pg2,…,pgi],i=96
Taking a set of load prediction curves P L, wherein pl i is a predicted value of the load at the ith moment, and the formula is as follows:
PL=[pl1,pl2,...,pli],i=96
And generating N sets of bias curves a= { a 1,a2,...,ai,...aN}T by using a gaussian distribution model, wherein a i is a set of matrices which contain 96 bias values and conform to gaussian distribution, and N j is a bias value. The formula is as follows:
ai=[n1,n2,…nj],nj~Normal(μj,σj 2),j∈[1,96]
And finally generating N groups of source load scene curves PG renewable,PL with different deviation degrees. Wherein:
S1.2, respectively loading the source load curves into different threads, and then loading scheduling strategies to be evaluated into the threads in parallel; and finally, carrying out 96-point parallel calculation to obtain an active cruising path result.
S2, constructing generalized sub-indexes including safety, economy and cleanliness, and calculating sub-indexes under different scenes. And finally, carrying out normalization treatment and huge treatment on the indexes.
The step S2 specifically refers to: the invention comprehensively considers the system cost, the cross-section power flow out-of-limit condition, the load balance, the new energy consumption rate, the out-of-limit power of the balancing machine and other factors to give out the calculation method of the index.
S2.1, system Cost th (positive rewards):
where P i is the decision output of the corresponding unit i at a certain moment, C i is the unit running cost of the corresponding unit, and the system has N units in total. C max is the normalization constant. The lower the system cost, the higher the score, with an upper limit of 100 points.
S2.2, section power flow out-of-limit Cost pf (positive rewards):
wherein S is the total section number, and r s is the prize value of the S-th section. The calculation method of r s is as follows:
Wherein P s is the power flow of section s, Is the upper limit of section s,/>Is the lower limit of the section s.
S2.3, load balancing Cost bal (positive rewards):
Wherein L represents the total system real load at time t, P i is the output of the unit i, and N is the total number of units. The prize is 0 when the load is balanced.
S2.4, new energy consumption rate (positive rewards):
Wherein, P i represents the output of the current new energy, P i max represents the predicted maximum output of the current new energy, and M new energy units are all arranged.
S3, establishing a fuzzy judgment matrix F by adopting an improved fuzzy AHP analytic hierarchy process, and calculating the weight of each index by using the AHP method to obtain subjective weight
The step S3 specifically comprises the following steps:
S3.1, constructing a fuzzy complementary judgment matrix F, namely constructing F m×n by a 1-9 scale method, wherein the definition of the 1-9 scale method is shown in the following table:
s3.2, calculating a weight formula by using the fuzzy complementary matrix, wherein the weight formula is as follows:
where F ij represents an element in the fuzzy complementary determining matrix F m×n, and W i represents a weight of the i-th index.
S3.3, carrying out consistency test through fuzzy complementary judgment matrix F:
the matrix F is first summed row by row as follows:
Then transformed by the result of the above equation: Finally obtaining a fuzzy consistency matrix R
Wherein, the relation of each element in R needs to satisfy:
0≤rij≤1,rii=0.5,rij+rji=1
And is also provided with With r ij=rik-rjk +0.5
S4, calculating the weight of each index by adopting an inverse entropy weight model to obtain objective weight
The step S4 specifically includes the following steps:
S4.1, forward processing of evaluation indexes: converting the minimum index, the intermediate index and the interval index into the maximum index; the formula is as follows:
very small- > very large: x i=max{x1,x2,...,xi}-xi. If all elements are positive numbers, the reciprocal can be directly taken
Intermediate- > very large: m=max { |x i-xbest | },
M=max{a-min{xi},max{xi}-b}
Regional- > very large:
S4.2, normalization (balancing because of differences between indices or differences in dimensions), where n evaluation objects, m indices, i.e. x=x n×m. The normalization formula is as follows:
S4.3, calculating the proportion of each element (the sum of the proportion is 1, namely, the normalized matrix is normalized), and calculating the probability matrix p, wherein the formula is as follows:
s4.4, calculating the information entropy e (uncertainty) of each index
S4.5, calculating the weight (1-e is the useful value of the information, and normalizing to obtain the weight)
S5 subjective weight calculated according to the steps S3 and S4And objective weight/>And calculating a variable weight coefficient omega i of the equalization function, and forming an index weight matrix omega.
The step S5 specifically includes the following steps:
S5.1, calculating the fixed weight
S5.2, calculating an evaluation value v i:
Wherein, And x ij is the evaluation value of the jth key index in the ith item.
S5.3, variable weight coefficient
S5.4 equalization function
S5.5, the variable weight coefficient of the equalization function is as follows:
Wherein a is an equalization factor.
S6, calculating membership functions of the generalization indexes by using the normalized result in the step S2.
The step S6 specifically includes the following steps:
S6.1, calculating membership functions of all indexes by utilizing the result x i obtained by normalizing and summing to be maximum in the step S2, wherein the fuzzy relation of all indexes and state evaluation indexes is reflected by adopting a normal membership function A (x i) according to the characteristic decision of the scheduling strategy active cruise generalization evaluation index
Where a represents the boundary of the membership function and σ represents the standard deviation of the normal curve.
S7, carrying out quantitative evaluation on the scheduling results of different source load scenes according to membership functions of various generalization indexes, and obtaining scheduling performance results under different source load scenes;
The step S7 specifically includes the following steps:
S7.1, calculating a single factor evaluation matrix R through the membership function in the step S6:
Wherein m is the number of generalization indexes.
S7.2, establishing a factor set of comprehensive evaluation: u= (U 1,u2,...,um), where element U 1 represents the ith factor affecting the evaluation object. Establishing an evaluation set of comprehensive evaluation: v= (V 1,v2,...,vm) where element V 1 represents the j-th evaluation result.
S7.3, building a comprehensive evaluation model:
B=Ω·R=(ω1,ω2,…ωm)·R=(b1,b2,b3,b4)
Wherein the method comprises the steps of
And S7.4, according to the maximum membership degree principle and the evaluation set V, obtaining the evaluation results of the scheduling strategy under different source load scenes through the evaluation model B.
According to the invention, a plurality of operation scenes of new energy and load prediction fluctuation are constructed according to an actual data curve, and a scheduling strategy is operated in a plurality of scenes in parallel by using a computer multithreading technology. Active scheduling results under a plurality of scenes are obtained; meanwhile, the weight of each sub-target of the scheduling strategy in the active scheduling cruising of the power grid is calculated by combining the AHP-inverse entropy weight model and an improved fuzzy analytic hierarchy process, and compared with the traditional analytic hierarchy process, the method is more accurate and scientific.
It is noted that relational terms such as first and second, and the like are used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Moreover, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or terminal that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or terminal. Without further limitation, an element defined by the statement "comprising … …" or "comprising … …" does not exclude the presence of additional elements in a process, method, article, or terminal device that includes the element. Further, herein, "greater than," "less than," "exceeding," and the like are understood to not include the present number; "above", "below", "within" and the like are understood to include this number.
While the embodiments have been described above, other variations and modifications will occur to those skilled in the art once the basic inventive concepts are known, and it is therefore intended that the foregoing description and drawings illustrate only embodiments of the invention and not limit the scope of the invention, and it is therefore intended that the invention not be limited to the specific embodiments described, but that the invention may be practiced with their equivalent structures or with their equivalent processes or with their use directly or indirectly in other related fields.
Claims (8)
1. The comprehensive performance evaluation method for the power grid rolling scheduling strategy is characterized by comprising the following steps of:
S1, selecting a plurality of new energy day-ahead output curves and load prediction curves, multiplying the curves by Gaussian distribution sampling curves with different mean values and variances, and generating a plurality of groups of source load curve scenes with different deviation degrees; the scheduling strategy is subjected to parallel decision calculation under the scenes, and a daily rolling scheduling result obtained by scheduling calculation is obtained;
S2, constructing generalized sub-indexes comprising safety, economy and cleanliness, calculating sub-indexes under different scenes, and finally carrying out normalization treatment and extremely-large treatment on the indexes;
s3, establishing a fuzzy judgment matrix F by adopting an improved fuzzy AHP analytic hierarchy process, and calculating the weight of each index by using the AHP method to obtain subjective weight
S4, calculating the weight of each index by adopting an inverse entropy weight model to obtain objective weight
S5, subjective weight calculated according to the steps S3 and S4And objective weight/>Calculating a variable weight coefficient omega i of the equalization function, and forming an index weight matrix omega;
S6, calculating membership functions of the generalization indexes by using the normalized result in the step S2;
and S7, carrying out quantitative evaluation on the scheduling results of different source load scenes according to membership functions of the generalization indexes, and obtaining scheduling performance results under different source load scenes.
2. The comprehensive performance evaluation method for the power grid rolling scheduling strategy according to claim 1, wherein the step S1 specifically includes the following steps:
S1.1, taking a group of new energy output curves PG renewable, wherein PG i is the output value of the new energy at the ith moment, and taking one-day output curves at 15 minutes intervals, so i takes 96 as shown in the following formula:
PGrenewable=[pg1,pg2,...,pgi],i=96
Taking a set of load prediction curves P L, wherein pl i is a predicted value of the load at the ith moment, and the formula is as follows:
PL=[pl1,pl2,...,pli],i=96
And generating N sets of bias curves a= { a 1,a2,...,ai,...aN}T by using a gaussian distribution model, wherein a i is a set of matrices which contain 96 bias values and conform to gaussian distribution, N j is a bias value, and the formula is as follows:
ai=[n1,n2,...nj],nj~Normal(μj,σj 2),j∈[1,96]
Finally, generating N groups of source load scene curves PG renewable,PL with different deviation degrees, wherein:
PL′=AN×96·PL96×1;PGrenewable′=AN×96·PGrenewable,96×1
S1.2, respectively loading the source load curves into different threads, and then loading scheduling strategies to be evaluated into the threads in parallel; and finally, carrying out 96-point parallel calculation to obtain an active cruising path result.
3. The comprehensive performance evaluation method for the power grid rolling scheduling strategy according to claim 1, wherein the step S2 specifically refers to: the calculation method for giving the index by comprehensively considering various factors such as system cost, section power flow out-of-limit condition, load balance, new energy consumption rate, balance machine out-of-limit power and the like comprises the following steps:
S2.1, system cost:
wherein P i is the decision output of the corresponding unit i at a certain moment, C i is the unit running cost of the corresponding unit, N units are in total in the system, C max is a normalization constant, the lower the system cost is, the higher the score is, and the upper limit is 100 score;
S2.2, section tide out-of-limit:
Wherein S is the total section number, r s is the rewarding value of the S-th section, and the calculation method of r s is as follows:
Wherein P s is the power flow of section s, Is the upper limit of section s,/>Is the lower limit of the section s;
s2.3, load balancing:
Wherein L represents the total system real load at the moment t, P i is the output of the unit i, N is the total number of the units, and the reward is 0 when the load is balanced;
s2.4, new energy consumption rate:
Wherein, P i represents the output of the current new energy, P i max represents the predicted maximum output of the current new energy, and M new energy units are all arranged.
4. The comprehensive performance evaluation method for the power grid rolling scheduling strategy according to claim 1, wherein the step S3 specifically includes the following steps:
S3.1, constructing a fuzzy complementary judgment matrix F, namely constructing F m×n by a 1-9 scale method, wherein the definition of the 1-9 scale method is shown in the following table:
s3.2, calculating a weight formula by using the fuzzy complementary matrix, wherein the weight formula is as follows:
Wherein F ij represents an element in the fuzzy complementary determining matrix F m×n, and W i represents a weight of the i-th index;
s3.3, carrying out consistency test through fuzzy complementary judgment matrix F:
the matrix F is first summed row by row as follows:
Then transformed by the result of the above equation: Finally obtaining a fuzzy consistency matrix R
Wherein, the relation of each element in R needs to satisfy:
0≤rij≤1,rii=0.5,rij+rji=1
And is also provided with There is r ij=rik-rjk +0.5.
5. The comprehensive performance evaluation method for the power grid rolling scheduling strategy according to claim 1, wherein the step S4 specifically comprises the following steps:
S4.1, forward processing of evaluation indexes: converting the minimum index, the intermediate index and the interval index into the maximum index; the formula is as follows:
Very small- > very large: x i=max{x1,x2,...,xi}-xi, if all elements are positive, the reciprocal can be directly taken
Intermediate- > very large:
M=max{a-min{xi},max{xi}-b}
regional- > very large:
S4.2, a normalization process, in which n evaluation objects, m indexes, i.e., x=x n×m, are normalized as follows:
s4.3, calculating the proportion of each element, wherein the sum of the proportion is 1, namely normalizing the standardized matrix, and calculating a probability matrix p, wherein the formula is as follows:
S4.4, calculating the information entropy e of each index
S4.5, calculating the weight, wherein 1-e is the useful value of the information, normalizing to obtain the weight,
6. The comprehensive performance evaluation method for the power grid rolling scheduling strategy according to claim 1, wherein the step S5 specifically comprises the following steps:
S5.1, calculating the fixed weight
S5.2, calculating an evaluation value v i:
Wherein, The constant weight coefficient of the jth key index in the ith item is given, and x ij is the evaluation value of the jth key index in the ith item;
S5.3, variable weight coefficient
S5.4 equalization function
S5.5, the variable weight coefficient of the equalization function is as follows:
Wherein a is an equalization factor.
7. The comprehensive performance evaluation method for the power grid rolling scheduling strategy according to claim 1, wherein the step S6 specifically includes the following steps:
Calculating membership functions of all indexes by using the result x i obtained by normalizing and summing the step S2 to be the maximum, wherein the fuzzy relation of all indexes and state evaluation indexes is reflected by adopting a normal membership function A (x i) according to the characteristic decision of the active cruise generalization evaluation indexes of the scheduling strategy
Where a represents the boundary of the membership function and σ represents the standard deviation of the normal curve.
8. The method for evaluating comprehensive performance of a rolling scheduling strategy for a power grid according to claim 7, wherein the step S7 specifically includes the following steps:
S7.1, calculating a single factor evaluation matrix R through the membership function in the step S6:
Wherein m is the number of generalization indexes;
S7.2, establishing a factor set of comprehensive evaluation: u= (U 1,u2,...,um) where element U 1 represents the ith factor affecting the evaluation object, an evaluation set of comprehensive evaluation is established: v= (V 1,v2,...,vm) where element V 1 represents the j-th evaluation result;
s7.3, building a comprehensive evaluation model:
B=ΩR=(ω1,ω2,…ωm)R=(b1,b2,b3,b4)
Wherein the method comprises the steps of
And S7.4, according to the maximum membership degree principle and the evaluation set V, obtaining the evaluation results of the scheduling strategy under different source load scenes through the evaluation model B.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410193952.3A CN118157104A (en) | 2024-02-21 | 2024-02-21 | Comprehensive performance evaluation method for power grid rolling scheduling strategy |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410193952.3A CN118157104A (en) | 2024-02-21 | 2024-02-21 | Comprehensive performance evaluation method for power grid rolling scheduling strategy |
Publications (1)
Publication Number | Publication Date |
---|---|
CN118157104A true CN118157104A (en) | 2024-06-07 |
Family
ID=91300848
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410193952.3A Pending CN118157104A (en) | 2024-02-21 | 2024-02-21 | Comprehensive performance evaluation method for power grid rolling scheduling strategy |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN118157104A (en) |
-
2024
- 2024-02-21 CN CN202410193952.3A patent/CN118157104A/en active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104408562B (en) | A kind of photovoltaic system generating efficiency comprehensive estimation method based on BP neural network | |
CN111178616B (en) | Wind speed prediction method based on negative correlation learning and regularization extreme learning machine integration | |
CN107274067B (en) | Distribution transformer overload risk assessment method | |
CN114648176B (en) | Wind-solar power supply consumption optimization method based on data driving | |
CN113361946B (en) | Power quality assessment method and device based on distributed photovoltaic grid-connected system | |
CN113344589B (en) | Intelligent identification method for collusion behavior of power generation enterprise based on VAEGMM model | |
CN105631532A (en) | Power system load prediction method using fuzzy decision-based neural network model | |
CN112308427A (en) | New energy consumption restriction factor evaluation method and system based on combined empowerment-grey correlation | |
CN115271253B (en) | Method and device for constructing water-wind-solar power generation power prediction model and storage medium | |
CN112149976B (en) | Power grid accurate investment project decision method based on artificial intelligence | |
CN117081063A (en) | Distributed charging load prediction method and system based on GCN-Crossformer model | |
CN115358437A (en) | Power supply load prediction method based on convolutional neural network | |
CN115879602A (en) | Ultra-short-term photovoltaic output prediction method based on transient weather | |
CN113065715A (en) | Multi-load ultra-short-term prediction method for comprehensive energy system | |
CN113762591B (en) | Short-term electric quantity prediction method and system based on GRU and multi-core SVM countermeasure learning | |
CN117078047B (en) | LSTM-based heat load prediction and distribution optimization method and system | |
CN114037209A (en) | Comprehensive benefit analysis method and device for distributed photovoltaic access direct-current power distribution system | |
CN108694475A (en) | Short-term time scale photovoltaic cell capable of generating power amount prediction technique based on mixed model | |
CN117113086A (en) | Energy storage unit load prediction method, system, electronic equipment and medium | |
CN118157104A (en) | Comprehensive performance evaluation method for power grid rolling scheduling strategy | |
CN115860552A (en) | Rural energy Internet comprehensive benefit evaluation method | |
CN111815085A (en) | Method for realizing comprehensive energy efficiency evaluation of rail transit energy management system | |
CN113537575A (en) | Method for predicting trend load of grid connection of distributed photovoltaic and electric automobile | |
CN113283638A (en) | Load extreme curve prediction method and system based on fusion model | |
CN112949936A (en) | Short-term photovoltaic power prediction method based on similar-day wavelet transform and multilayer perceptron |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |