CN117154852A - Power grid economic dispatching method based on alternate iteration under consideration of N-1 tide constraint - Google Patents
Power grid economic dispatching method based on alternate iteration under consideration of N-1 tide constraint Download PDFInfo
- Publication number
- CN117154852A CN117154852A CN202311403718.0A CN202311403718A CN117154852A CN 117154852 A CN117154852 A CN 117154852A CN 202311403718 A CN202311403718 A CN 202311403718A CN 117154852 A CN117154852 A CN 117154852A
- Authority
- CN
- China
- Prior art keywords
- constraint
- power
- line
- solution
- key
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000005457 optimization Methods 0.000 claims abstract description 56
- 230000008569 process Effects 0.000 claims abstract description 6
- 239000000243 solution Substances 0.000 claims description 67
- 230000005283 ground state Effects 0.000 claims description 13
- 238000010248 power generation Methods 0.000 claims description 11
- 238000012546 transfer Methods 0.000 claims description 9
- 230000002068 genetic effect Effects 0.000 claims description 7
- 238000002347 injection Methods 0.000 claims description 6
- 239000007924 injection Substances 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 6
- 108010074506 Transfer Factor Proteins 0.000 claims description 3
- 230000035772 mutation Effects 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 3
- 239000013589 supplement Substances 0.000 claims description 3
- 230000009977 dual effect Effects 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 239000004820 Pressure-sensitive adhesive Substances 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000009194 climbing Effects 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009467 reduction 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- 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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas 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/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
-
- 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/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- 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]
-
- 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
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/22—The renewable source being solar energy
- H02J2300/24—The renewable source being solar energy of photovoltaic origin
-
- 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
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/28—The renewable source being wind energy
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Power Engineering (AREA)
- Life Sciences & Earth Sciences (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Biophysics (AREA)
- Evolutionary Biology (AREA)
- General Health & Medical Sciences (AREA)
- Economics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Genetics & Genomics (AREA)
- Data Mining & Analysis (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computational Linguistics (AREA)
- Physiology (AREA)
- Software Systems (AREA)
- Biomedical Technology (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The application discloses a power grid economic dispatching method based on alternate iteration under consideration of N-1 tide constraint. According to the method, all base state power flow constraints and N-1 power flow constraints in a power grid safety economic dispatching problem are ignored, an optimization target for minimizing the maximum line load rate is introduced into a model, and an economic better solution meeting safety requirements is obtained through solving the model in a hopeful and efficient manner and is used as an initial point for a follow-up further optimization dispatching plan; the follow-up optimization process adopts the idea of alternate iterative optimization, firstly, based on a safe economic dispatch model with simplified constraint, a large number of constraint conditions covered in the power grid safe economic dispatch model are processed by alternate iterative execution of a model solving-constraint checking-constraint supplementing-model perfecting mode, the problem of solving difficulty caused by considering a large number of constraints in the model at one time is avoided, the overall solving speed is accelerated, and the requirement of efficiently making a safe economic dispatch plan in a power system containing wind power is met.
Description
Technical Field
The application belongs to the technical field of economic dispatch of power systems, and particularly relates to a power grid economic dispatch method based on alternate iteration under consideration of N-1 tide constraint.
Background
In order to adjust the energy structure, the permeability of new energy, represented by wind power generation, in the electric power system is continuously increasing. Unlike conventional thermal power generation, wind power output has strong randomness and intermittence, and wind power output can be accurately predicted only in a short-term time scale, for example, wind power output after one hour in the future is predicted, and when the prediction period is increased, the prediction error is obviously increased. In order to reduce the influence of wind power output uncertainty on safe and economic operation of a power grid, safe and economic scheduling decisions can be made based on relatively accurate wind power output prediction data obtained in a short-term time scale. However, since a certain time is required for collecting the wind power output prediction data, and a certain time margin is required for scheduling decision, namely, a formulated scheduling plan is required to be issued after a certain time is extracted, the time left for scheduling decision after the wind power prediction data under a short-term time scale is acquired is very short, and the efficiency of scheduling decision is very important.
The power grid safety economic dispatching decision-making problem belongs to the large-scale complex nonlinear programming problem. The existing solving method is concentrated on a direct optimization method, namely, the running cost of a power grid is taken as an objective function, a dispatching optimization model is constructed by taking constraints such as the base state power flow of the power grid into consideration, and then, the solution is carried out through nonlinear programming or heuristic algorithm and the like, so that an economic dispatching strategy meeting the constraints is obtained.
In the prior art, an economic scheduling method of an active power distribution network based on opportunity constraint
CN 110336308B), comprising: an active power distribution network economic dispatching model considering wind power output uncertainty is established, the active power distribution network economic dispatching model takes the minimized power distribution network running cost as an objective function, and active power distribution network tide constraint, unit output constraint, unit climbing constraint, system maximum rotation standby constraint, energy storage constraint and wind power output opportunity constraint are taken as constraint conditions; and solving an economic dispatch model of the active power distribution network.
In the prior art, a regional power grid dispatching optimization method and system containing source load uncertainty utilizes an improved Gaussian fitting model to generate scenes, builds scenes of wind power, photovoltaic and load double-side uncertainty, combines a quantile regression method, and performs scene reduction according to dispatching requirements of different regions; and taking the system operation cost into consideration, constructing an interval objective function and constraint conditions, optimizing the objective function to be the minimum system operation cost, taking constraints such as unit constraint, start-stop constraint, wind-light output constraint, power flow constraint and the like as constraint conditions, and solving the model by adopting a random programming method to obtain a final scheduling optimization strategy (CN 113852069A).
However, when considering the N-1 power flow constraint, the constraint condition scale in the model is greatly increased, and at this time, all the constraints in the disposable processing model become difficult to process in terms of calculation, so that the algorithm convergence speed is slow, the solving time is long, and the requirement of efficient decision cannot be met.
In view of the above, there is a need to solve the technical problems in the related art.
Disclosure of Invention
In order to solve the problem that the safety economic dispatch decision efficiency is low after N-1 power flow constraint is considered, the application provides a power grid economic dispatch plan optimization method based on alternate iteration under N-1 power flow constraint.
The object of the application is achieved by at least one of the following technical solutions.
The power grid economic dispatching method based on alternate iteration under the consideration of N-1 tide constraint comprises the following steps:
s1, constructing a double-target optimization model M1 by taking minimized conventional unit power generation cost and minimized maximum line load rate as objective functions and considering conventional unit output constraint and power balance constraint;
s2, establishing a ground state power flow constraint and an N-1 power flow constraint which need to be considered in power grid economic dispatch, and forming a constraint set C1;
s3, solving the model M1 by adopting a multi-objective genetic algorithm to obtain a group of optimal solution setsPS;
S4, correspondingly generating the power at the cost from small to largePSIs a solution to the ordering of (2);
s5, sequentially taking out the optimal solution setPSSolution in (a)S;
S6, judging solutionSWhether the constraints in constraint set C1 are violated,if yes, recording the violated constraint, forming a key constraint set C2, and entering a step S7, otherwise, entering a step S8;
s7, judging whether the traversal is finishedPSIf yes, the step S8 is carried out, if not, the step S5 is returned;
s8, outputting solutionS;
S9, taking the minimized conventional unit power generation cost as an objective function, and constructing an optimization model M2 taking the key constraint set into consideration by taking the conventional unit output constraint, the power balance constraint and the N-1 constraint in the key constraint set C2 into consideration;
s10, bySFor initial solution,Obtaining an optimal solution by solving an optimization model M2 considering a key constraint set through nonlinear programmingS’;
S11, judging the optimal solutionS’If yes, executing step S13, and if no, executing step S12;
s12, judging the optimal solutionS’Whether the constraint in the constraint set C1 is violated, if yes, adding the violated constraint into the key constraint set C2 and returning to the step S9, and if not, outputting an optimal solutionS’Obtaining an economic optimal scheduling plan meeting all constraints, and performing power grid scheduling according to the economic optimal scheduling plan;
s13, finishing searching, and outputting that no feasible solution exists.
Further, in step S1, the objective function in the dual-objective optimization model M1 is expressed as follows:
wherein,f 1 the output cost of a conventional unit is represented;gnumbering the conventional units;NGthe number of the conventional units;P g for conventional units in dispatch plansgIs an active force of (a);a g 、b g 、c g are all conventional unitsgIs a coefficient of force;f 2 indicating maximum line load factorkThe sum of the load rates of the lines;F l representation linelThe specific expression of the line flow of the (C) is as follows:
Wherein,G l-g is a conventional unitgLine-to-linelA power injection transfer profile of (2);wnumbering the wind farm;NWthe number of wind farms in the grid;P w is a power generation fieldwIs a predicted force of (1);G l-w for wind farmswLine-to-linelA power injection transfer profile of (2);nnumbering the load nodes;G l-n is a load nodenLine-to-linelNode load transfer distribution factors of (a);d n is a nodenIs a load of (2);
the double-objective optimization model M1 includes the following constraint conditions:
1) Power balance constraint:
2) Conventional unit output constraint:
wherein,P g,min andP g,max respectively conventional unitsgAnd a maximum output.
Further, in step S2, the specific expression of the ground state power flow constraint is as follows:
wherein,F l,max is a circuitlIs a trend limit value of (2);
the specific expression of the N-1 tide constraint is as follows:
wherein,cpredicting a faulty line for the N-1 under consideration;F l,c representation linecLine after N-1 faultlIs a trend of (2);is a circuitcLine-to-linelIs a power flow transfer factor of (1).
Further, the specific implementation flow of step S3 is as follows:
s3.1, setting parameters of a multi-target genetic algorithm: population numberHCross probability factorp m Variation probability factorp c Maximum number of iterationsN max ;
S3.2, randomly generating a parent population;
s3.3, performing cross mutation operation on the parent population to obtain a child population, and combining the child population and the parent population;
s3.4, calculating the fitness of all individuals of the population, comparing the merits of the individuals according to the Pareto dominant relationship, and reserving the superior individuals to form a next generation parent population;
s3.5, removing individuals with smaller crowding degree distance from the parent population to control the parent population scale to be kept asHThe crowding degree distance formula is:
wherein,D i representing individuals in a populationiIs used for the distance of the degree of congestion,D i smaller, indicate individualiThe more crowded the distribution of (2);hrepresenting the number of optimization targets;representing individualsiAlong the first edgemThe horizontal distance between two sides of each optimization target and adjacent individuals;
and S3.6, judging whether the iteration times meet the termination condition, if so, outputting the solved set, otherwise, returning to the step S3.3.
Further, in step S2, the N-1 contemplated faulty linecIs preassigned and satisfies the disconnection linecThe lower network has no island node condition.
Further, in step S6, by checking the optimal solution setPSIdentifying and extracting a key constraint set C2 in the power grid ground state power flow constraint and the N-1 power flow constraint;
checking optimal solution setPSIn the process of (1), violation conditions of a plurality of constraints are analyzed simultaneously in a parallel computing mode.
Further, in step S9, the constraint conditions of the optimization model M2 constructed by considering the key constraint set only consider the conventional unit output constraint, the power balance constraint and the key constraint set C2, and the ground state power flow constraint and the N-1 power flow constraint outside the key constraint set C2 are ignored.
Further, in step S10, the solution obtained in step S8 is always usedSAs a starting point for the optimization model M2 taking into account the set of key constraints by non-linear programming solution.
Further, in step S10, if the solution isSAnd if the constraint in the constraint set C1 is not violated, solving the optimization model M2 considering the key constraint set by an interior point method, otherwise solving the optimization model M2 considering the key constraint set by an exterior point method.
Further, in step S11, if the test finds the optimal solution of the optimization model M2 considering the key constraint setS’If all constraints are not satisfied, recording the optimal solutionS’Violated constraints and supplement these into the optimization model M2, which considers the set of key constraints.
Compared with the prior art, the application has the following beneficial effects:
aiming at the problem that the safety economic dispatching decision efficiency is low after N-1 power flow constraint is considered, the application provides a power grid economic dispatching method based on alternate iteration under N-1 power flow constraint. According to the application, all ground state power flow constraints and N-1 power flow constraints in the power grid safety economic dispatching problem are ignored, an optimization target for minimizing the maximum line load rate is introduced into the model, the complexity of the model is reduced, meanwhile, the economical efficiency and the safety requirements are considered, and the model is solved to obtain an economical efficiency preferred solution meeting the safety requirements in a high-efficiency manner, and the economic efficiency preferred solution is used as an initial point for further optimizing a dispatching plan. In the subsequent optimization process, the method adopts the idea of alternate iterative optimization, firstly, the safe economic dispatch model after constraint simplification is based, and then, a large number of constraint conditions covered in the power grid safe economic dispatch model are processed by means of alternate iterative execution of model solving-constraint checking-constraint supplementing-model perfecting, so that the problem of solving difficulty caused by considering a large number of constraints in the model at one time is avoided, the overall solving speed is increased, and the requirement of efficiently making the safe economic dispatch plan in the wind power system is met.
Drawings
FIG. 1 is a schematic flow diagram of a power grid economic dispatch method based on alternate iteration under consideration of N-1 tide constraint provided by an embodiment of the application;
FIG. 2 is a diagram of an improved IEEE30 node system test grid structure in a grid economic dispatching method based on alternate iteration under consideration of N-1 tide constraint provided by an embodiment of the application;
fig. 3 is a schematic flow chart of solving a double-objective optimization model M1 by adopting a multi-objective genetic algorithm in the power grid economic dispatching method based on alternate iteration under consideration of N-1 tide constraint provided by the embodiment of the application.
Detailed Description
The present application will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present application more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the application.
Examples:
considering the power grid economic dispatching method based on alternate iteration under the constraint of N-1 tide, as shown in fig. 1, the method comprises the following steps:
s1, constructing a double-target optimization model M1 by taking minimized conventional unit power generation cost and minimized maximum line load rate as objective functions and considering conventional unit output constraint and power balance constraint;
the objective function in the dual-objective optimization model M1 is expressed as follows:
wherein,f 1 the output cost of a conventional unit is represented;gnumbering the conventional units;NGthe number of the conventional units;P g for conventional units in dispatch plansgIs an active force of (a);a g 、b g 、c g are all conventional unitsgIs a coefficient of force;f 2 indicating maximum line load factorkThe sum of the load rates of the lines;F l representation linelThe specific expression of the line flow is as follows:
wherein,G l-g is a conventional unitgLine-to-linelA power injection transfer profile of (2);wnumbering the wind farm;NWthe number of wind farms in the grid;P w is a power generation fieldwIs a predicted force of (1);G l-w for wind farmswLine-to-linelA power injection transfer profile of (2);nnumbering the load nodes;G l-n is a load nodenLine-to-linelNode load transfer distribution factors of (a);d n is a nodenIs a load of (2);
the double-objective optimization model M1 includes the following constraint conditions:
1) Power balance constraint:
2) Conventional unit output constraint:
wherein,P g,min andP g,max respectively conventional unitsgIs the minimum of (2)Output and maximum output.
In one embodiment, in the test system, the conventional crew at node 22 in the original IEEE30 node system is replaced with a wind farm and assumed to be predicted to have an output of 25MW, as shown in FIG. 2. The load of each node in the system is shown in table 1.
Table 1 load table for each node in test system
S2, establishing a ground state power flow constraint and an N-1 power flow constraint which need to be considered in power grid economic dispatch, and forming a constraint set C1;
the specific expression of the ground state tide constraint is as follows:
wherein,F l,max is a circuitlIs a power flow limit value of (2).
The specific expression of the N-1 tide constraint is as follows:
wherein,cpredicting a faulty line for the N-1 under consideration;F l,c representation linecLine after N-1 faultlIs a trend of (2);is a circuitcLine-to-linelIs a power flow transfer factor of (1).
N-1 contemplated fault linecIs preassigned and satisfies the disconnection linecThe lower network has no island node condition.
In one embodiment, N-1 expected failure lines for improved IEEE39 node system considerations are shown in Table 2.
Table 2N-1 expected failure line set table for improved IEEE39 node system considerations
S3, solving the model M1 by adopting a multi-objective genetic algorithm to obtain a group of optimal solution setsPSAs shown in fig. 3, the specific implementation flow is as follows:
s3.1, setting parameters of a multi-target genetic algorithm: population numberHCross probability factorp m Variation probability factorp c Maximum number of iterationsN max The method comprises the steps of carrying out a first treatment on the surface of the In one embodiment, a settingHAt the point of 50 a,p m the total number of the components is 0.9,p c the content of the acid in the solution is 0.1,N max 100;
s3.2, randomly generating a parent population;
s3.3, performing cross mutation operation on the parent population to obtain a child population, and combining the child population and the parent population;
s3.4, calculating the fitness of all individuals of the population, comparing the merits of the individuals according to the Pareto dominant relationship, and reserving the superior individuals to form a next generation parent population;
s3.5, removing individuals with smaller crowding degree distance from the parent population to control the parent population scale to be kept asHThe crowding degree distance formula is:
wherein,D i representing individuals in a populationiIs used for the distance of the degree of congestion,D i smaller, indicate individualiThe more crowded the distribution of (2);hrepresenting the number of optimization targets;representing individualsiAlong the first edgemThe horizontal distance between two sides of each optimization target and adjacent individuals;
and S3.6, judging whether the iteration times meet the termination condition, if so, outputting the solved set, otherwise, returning to the step S3.3.
S4, correspondingly generating the power at the cost from small to largePSIs a solution to the ordering of (2);
s5, sequentially taking out the optimal solution setPSSolution in (a)S;
S6, judging solutionSIf yes, recording the violated constraint in the constraint set C1, forming a key constraint set C2, and entering a step S7, otherwise, entering a step S8;
by checking the optimal solution setPSIdentifying and extracting a key constraint set C2 in the power grid ground state power flow constraint and the N-1 power flow constraint;
checking optimal solution setPSIn the process of (1), violation conditions of a plurality of constraints are analyzed simultaneously in a parallel computing mode.
S7, judging whether the traversal is finishedPSIf yes, the step S8 is carried out, if not, the step S5 is returned;
s8, outputting solutionS;
S9, taking the minimized conventional unit power generation cost as an objective function, and constructing an optimization model M2 taking the key constraint set into consideration by taking the conventional unit output constraint, the power balance constraint and the N-1 constraint in the key constraint set C2 into consideration;
the constraint conditions of the constructed optimization model M2 considering the key constraint set only consider the conventional unit output constraint, the power balance constraint and the key constraint set C2, and the ground state power flow constraint and the N-1 power flow constraint outside the key constraint set C2 are ignored.
S10, bySFor initial solution,Obtaining an optimal solution by solving an optimization model M2 considering a key constraint set through nonlinear programmingS’;
Always with the solution obtained in step S8SAs a starting point for the optimization model M2 taking into account the set of key constraints by non-linear programming solution.
If solveSAnd if the constraint in the constraint set C1 is not violated, solving the optimization model M2 considering the key constraint set by an interior point method, otherwise solving the optimization model M2 considering the key constraint set by an exterior point method.
S11, judging the optimal solutionS’If yes, executing step S13, and if no, executing step S12;
if the test finds that the optimization model M2 considers the key constraint setOptimal solutionS’If all constraints are not satisfied, recording the optimal solutionS’Violated constraints and supplement these into the optimization model M2, which considers the set of key constraints.
S12, judging the optimal solutionS’Whether the constraint in the constraint set C1 is violated, if yes, adding the violated constraint into the key constraint set C2 and returning to the step S9, and if not, outputting an optimal solutionS’Obtaining an economical optimal scheduling plan meeting all constraints;
s13, finishing searching, and outputting that no feasible solution exists.
In one embodiment, the resulting solution is shown in Table 3.
Table 3 table of the solution results using the method provided by the present application in the example
To better verify the effectiveness of the present application, the alternate optimization method provided by the present application is compared with a direct optimization method that considers all constraints in the model once. The solution results of the direct optimization method are shown in table 4.
Table 4 Table of the solution results using the direct optimization method in the example
Comparing tables 3 and 4, the application is verified to be capable of efficiently obtaining the economic dispatch plan meeting the safety constraint, and meeting the requirement of high-efficiency decision of the safety economic dispatch in the wind power system.
The above description is only of the preferred embodiments of the present application and is not intended to limit the present application, but various modifications and variations can be made to the present application by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present application should be included in the protection scope of the present application.
Claims (10)
1. The power grid economic dispatching method based on alternate iteration under the consideration of N-1 tide constraint is characterized by comprising the following steps:
s1, constructing a double-target optimization model M1 by taking minimized conventional unit power generation cost and minimized maximum line load rate as objective functions and considering conventional unit output constraint and power balance constraint;
s2, establishing a ground state power flow constraint and an N-1 power flow constraint which need to be considered in power grid economic dispatch, and forming a constraint set C1;
s3, solving the model M1 by adopting a multi-objective genetic algorithm to obtain a group of optimal solution setsPS;
S4, correspondingly generating the power at the cost from small to largePSIs a solution to the ordering of (2);
s5, sequentially taking out the optimal solution setPSSolution in (a)S;
S6, judging solutionSIf yes, recording the violated constraint in the constraint set C1, forming a key constraint set C2, and entering a step S7, otherwise, entering a step S8;
s7, judging whether the traversal is finishedPSIf yes, the step S8 is carried out, if not, the step S5 is returned;
s8, outputting solutionS;
S9, taking the minimized conventional unit power generation cost as an objective function, and constructing an optimization model M2 taking the key constraint set into consideration by taking the conventional unit output constraint, the power balance constraint and the N-1 constraint in the key constraint set C2 into consideration;
s10, bySFor initial solution,Obtaining an optimal solution by solving an optimization model M2 considering a key constraint set through nonlinear programmingS’;
S11, judging the optimal solutionS’If yes, executing step S13, and if no, executing step S12;
s12, judging the optimal solutionS’Whether the constraint in the constraint set C1 is violated, if yes, adding the violated constraint into the key constraint set C2 and returning to the step S9, and if not, outputting an optimal solutionS’Obtaining an economic optimal scheduling plan meeting all constraints, and performing power grid scheduling according to the economic optimal scheduling plan;
s13, finishing searching, and outputting that no feasible solution exists.
2. The grid economic dispatch method based on alternate iteration under consideration of N-1 power flow constraint of claim 1, wherein in step S1, the objective function in the dual objective optimization model M1 is expressed as follows:
wherein,f 1 the output cost of a conventional unit is represented;gnumbering the conventional units;NGthe number of the conventional units;P g for conventional units in dispatch plansgIs an active force of (a);a g 、b g 、c g are all conventional unitsgIs a coefficient of force;f 2 indicating maximum line load factorkThe sum of the load rates of the lines;F l representation linelThe specific expression of the line flow is as follows:
wherein,G l-g is a conventional unitgLine-to-linelA power injection transfer profile of (2);wnumbering the wind farm;NWthe number of wind farms in the grid;P w is a power generation fieldwIs a predicted force of (1);G l-w for wind farmswLine-to-linelA power injection transfer profile of (2);nnumbering the load nodes;G l-n is a load nodenLine-to-linelNode load transfer distribution factors of (a);d n is a nodenIs a load of (2);
the double-objective optimization model M1 includes the following constraint conditions:
1) Power balance constraint:
2) Conventional unit output constraint:
wherein,P g,min andP g,max respectively conventional unitsgAnd a maximum output.
3. The grid economic dispatch method based on alternate iteration under consideration of N-1 power flow constraint according to claim 1, wherein in step S2, the specific expression of the ground state power flow constraint is as follows:
wherein,F l,max is a circuitlIs a trend limit value of (2);
the specific expression of the N-1 tide constraint is as follows:
wherein,cpredicting a faulty line for the N-1 under consideration;F l,c representation linecLine after N-1 faultlIs a trend of (2);is a circuitcLine-to-linelIs a power flow transfer factor of (1).
4. The grid economic dispatch method based on alternate iteration under consideration of N-1 tide constraint according to claim 1, wherein the specific implementation flow of step S3 is as follows:
s3.1, setting parameters of a multi-target genetic algorithm: population numberHCross probability factorp m Variation probability factorp c Maximum number of iterationsN max ;
S3.2, randomly generating a parent population;
s3.3, performing cross mutation operation on the parent population to obtain a child population, and combining the child population and the parent population;
s3.4, calculating the fitness of all individuals of the population, comparing the merits of the individuals according to the Pareto dominant relationship, and reserving the superior individuals to form a next generation parent population;
s3.5, removing individuals with smaller crowding degree distance from the parent population to control the parent population scale to be kept asHThe crowding degree distance formula is:
wherein,D i representing individuals in a populationiIs used for the distance of the degree of congestion,D i smaller, indicate individualiThe more crowded the distribution of (2);hrepresenting the number of optimization targets;representing individualsiAlong the first edgemThe horizontal distance between two sides of each optimization target and adjacent individuals;
and S3.6, judging whether the iteration times meet the termination condition, if so, outputting the solved set, otherwise, returning to the step S3.3.
5. The grid economic dispatch method based on alternate iteration under consideration of N-1 tide constraint of claim 1, wherein in step S2, N-1 considered envisions a faulty linecIs preassigned and satisfies the disconnection linecThe lower network has no island node condition.
6. The grid economic dispatch method based on alternate iteration under consideration of N-1 tide constraint of claim 1, wherein in step S6, the optimal solution set is verifiedPSIs to identify the solution in (a)Extracting a key constraint set C2 in the power grid ground state power flow constraint and the N-1 power flow constraint;
checking optimal solution setPSIn the process of (1), violation conditions of a plurality of constraints are analyzed simultaneously in a parallel computing mode.
7. The grid economic dispatch method based on alternate iteration under consideration of N-1 power flow constraint according to claim 1, wherein in step S9, the constraint condition of the constructed optimization model M2 considering the key constraint set only considers the conventional unit output constraint, the power balance constraint and the key constraint set C2, and the ground state power flow constraint and the N-1 power flow constraint outside the key constraint set C2 are ignored.
8. The grid economic dispatch method based on alternate iteration under consideration of N-1 tide constraint according to claim 1, wherein in step S10, the solution obtained in step S8 is always usedSAs a starting point for the optimization model M2 taking into account the set of key constraints by non-linear programming solution.
9. The grid economic dispatch method based on alternate iteration under consideration of N-1 tidal current constraint of claim 1, wherein in step S10, if the solution isSAnd if the constraint in the constraint set C1 is not violated, solving the optimization model M2 considering the key constraint set by an interior point method, otherwise solving the optimization model M2 considering the key constraint set by an exterior point method.
10. The grid economic dispatch method based on alternate iteration under consideration of N-1 tidal current constraint according to claim 1, wherein in step S11, if the test finds the optimal solution of the optimization model M2 considering the key constraint setS’If all constraints are not satisfied, recording the optimal solutionS’Violated constraints and supplement these into the optimization model M2, which considers the set of key constraints.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311403718.0A CN117154852B (en) | 2023-10-27 | 2023-10-27 | Power grid economic dispatching method based on alternate iteration under consideration of N-1 tide constraint |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311403718.0A CN117154852B (en) | 2023-10-27 | 2023-10-27 | Power grid economic dispatching method based on alternate iteration under consideration of N-1 tide constraint |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117154852A true CN117154852A (en) | 2023-12-01 |
CN117154852B CN117154852B (en) | 2024-01-02 |
Family
ID=88910393
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311403718.0A Active CN117154852B (en) | 2023-10-27 | 2023-10-27 | Power grid economic dispatching method based on alternate iteration under consideration of N-1 tide constraint |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117154852B (en) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103473604A (en) * | 2013-08-22 | 2013-12-25 | 国家电网公司 | Method for evaluating comprehensive property of day-ahead power generation plan |
US20150378387A1 (en) * | 2013-05-14 | 2015-12-31 | Rensselaer Polytechnic Institute | Methods of computing steady-state voltage stability margins of power systems |
US20150378384A1 (en) * | 2014-06-29 | 2015-12-31 | International Business Machines Corporation | Reducing corrective actions in security-constrained optimal power flow via sparsity regularization |
CN106786529A (en) * | 2016-12-23 | 2017-05-31 | 国网山西省电力公司电力科学研究院 | A kind of distributed static security analysis method |
CN110880072A (en) * | 2019-11-15 | 2020-03-13 | 国网湖南省电力有限公司 | Real-time power grid static security risk disposal optimization method and device and storage medium |
CN114142530A (en) * | 2021-09-23 | 2022-03-04 | 重庆大学 | Risk scheduling method considering N-1 security constraint based on near-end policy optimization algorithm |
CN116934021A (en) * | 2023-07-17 | 2023-10-24 | 西安理工大学 | Electric-thermal comprehensive energy system optimal scheduling method considering N-1 safety constraint |
-
2023
- 2023-10-27 CN CN202311403718.0A patent/CN117154852B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150378387A1 (en) * | 2013-05-14 | 2015-12-31 | Rensselaer Polytechnic Institute | Methods of computing steady-state voltage stability margins of power systems |
CN103473604A (en) * | 2013-08-22 | 2013-12-25 | 国家电网公司 | Method for evaluating comprehensive property of day-ahead power generation plan |
US20150378384A1 (en) * | 2014-06-29 | 2015-12-31 | International Business Machines Corporation | Reducing corrective actions in security-constrained optimal power flow via sparsity regularization |
CN106786529A (en) * | 2016-12-23 | 2017-05-31 | 国网山西省电力公司电力科学研究院 | A kind of distributed static security analysis method |
CN110880072A (en) * | 2019-11-15 | 2020-03-13 | 国网湖南省电力有限公司 | Real-time power grid static security risk disposal optimization method and device and storage medium |
CN114142530A (en) * | 2021-09-23 | 2022-03-04 | 重庆大学 | Risk scheduling method considering N-1 security constraint based on near-end policy optimization algorithm |
CN116934021A (en) * | 2023-07-17 | 2023-10-24 | 西安理工大学 | Electric-thermal comprehensive energy system optimal scheduling method considering N-1 safety constraint |
Non-Patent Citations (1)
Title |
---|
肖逸 等: "基于微服务的大电网在线安全稳定分析系统数据中台与架构设计", 中国电机工程学报, vol. 42, pages 126 - 135 * |
Also Published As
Publication number | Publication date |
---|---|
CN117154852B (en) | 2024-01-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112217202B (en) | Distributed new energy, energy storage and power distribution network planning method considering flexibility investment | |
CN112186765B (en) | Modeling method of day-ahead scheduling model of unit combination decision | |
Ghadimi et al. | PSO based fuzzy stochastic long-term model for deployment of distributed energy resources in distribution systems with several objectives | |
CN107316113B (en) | Power transmission network planning method and system | |
CN103683326A (en) | Method for calculating optimal admitting ability for wind power multipoint access of regional power grid | |
CN105790265A (en) | AC power flow constraint-based uncertainty unit commitment model and solving method | |
CN111738519A (en) | Power distribution network planning method, system and equipment | |
CN107977744A (en) | A kind of electric system based on traditional Benders decomposition methods Robust Scheduling method a few days ago | |
CN110854891B (en) | Power distribution network pre-disaster resource allocation method and system | |
CN113393339A (en) | Power distribution network short-term load prediction method based on deep neural network combination model | |
CN115439000A (en) | Power distribution network block division method considering wind-solar-load power uncertainty and correlation | |
CN113344283B (en) | Energy internet new energy consumption capability assessment method based on edge intelligence | |
CN112926762B (en) | Clearing method and device based on safety constraint unit combination model | |
CN104915788B (en) | A method of considering the Electrical Power System Dynamic economic load dispatching of windy field correlation | |
CN117154852B (en) | Power grid economic dispatching method based on alternate iteration under consideration of N-1 tide constraint | |
CN113705964B (en) | Method and device for formulating pre-disaster plan for toughness recovery of power distribution network | |
CN112436514B (en) | Multi-microgrid interconnection optimization method considering photovoltaic uncertainty | |
CN115146870A (en) | Method, system, equipment and medium for generating day-to-day power supply plan of isolated island | |
CN114239976A (en) | Distribution network distributed energy storage planning method considering typhoon weather influence | |
CN112152267B (en) | Power grid random reactive power optimization scheduling method considering source load uncertainty | |
CN110571791B (en) | Optimal configuration method for power transmission network planning under new energy access | |
CN112446521B (en) | Multi-objective planning method for wind power plant access system considering economy and safety | |
Xie et al. | Optimal Sizing of Renewables and Energy Storage System in Microgrids: The Quest for Economic Uncertainty Management | |
Zhou et al. | Multi-Time scale unit commitment optimization under hybrid uncertainties | |
Chi et al. | Research on Distribution Network Dispatching Control Strategy based on Power Spot Market |
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 |