CN117411082A - Combined optimization solving method and device for safety constraint unit - Google Patents
Combined optimization solving method and device for safety constraint unit Download PDFInfo
- Publication number
- CN117411082A CN117411082A CN202311406837.1A CN202311406837A CN117411082A CN 117411082 A CN117411082 A CN 117411082A CN 202311406837 A CN202311406837 A CN 202311406837A CN 117411082 A CN117411082 A CN 117411082A
- Authority
- CN
- China
- Prior art keywords
- line
- constraint
- aggregation
- safety constraint
- safety
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000005457 optimization Methods 0.000 title claims abstract description 26
- 230000002776 aggregation Effects 0.000 claims abstract description 158
- 238000004220 aggregation Methods 0.000 claims abstract description 158
- 239000011159 matrix material Substances 0.000 claims abstract description 97
- 238000009826 distribution Methods 0.000 claims abstract description 44
- 238000012546 transfer Methods 0.000 claims abstract description 31
- 239000013598 vector Substances 0.000 claims description 66
- 230000006870 function Effects 0.000 claims description 58
- 230000009194 climbing Effects 0.000 claims description 21
- 108010074506 Transfer Factor Proteins 0.000 claims description 13
- 238000004590 computer program Methods 0.000 claims description 11
- 238000003860 storage Methods 0.000 claims description 10
- 238000005259 measurement Methods 0.000 claims description 9
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 6
- 230000003247 decreasing effect Effects 0.000 claims description 5
- 230000002349 favourable effect Effects 0.000 abstract 1
- 238000004891 communication Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000004146 energy storage Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000013473 artificial intelligence Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000001953 sensory effect Effects 0.000 description 1
- 238000010977 unit operation Methods 0.000 description 1
- 230000000007 visual effect 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- 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
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/24—Arrangements for preventing or reducing oscillations of power in 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
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Power Engineering (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Operations Research (AREA)
- Computing Systems (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The application provides a method and a device for combined optimization solving of a safety constraint unit, wherein the method comprises the following steps: acquiring a line safety constraint space based on a source power transfer distribution factor matrix, and measuring the similarity of the line safety constraint in the line safety constraint space by calculating an included angle cosine value matrix; acquiring a line safety constraint similarity matrix according to the similarity and a line similarity threshold, and acting on the similarity matrix through an integer programming model to determine a specific line aggregation mode so as to obtain an aggregation line safety constraint set; acquiring economic operation basic data of an electric power system, and constructing a safety constraint unit combination model applicable to line safety constraint aggregation according to an aggregation line safety constraint set; and solving the safety constraint unit combination model by adopting a mixed integer linear programming solver to obtain a line safety constraint aggregation unit combination result. The method and the device are favorable for improving the optimal solving efficiency of the power system and guaranteeing the stable operation of the power system.
Description
Technical Field
The application relates to the technical field of power systems, in particular to a method and a device for combined optimization solving of a safety restraint unit.
Background
The safety constraint unit combination is a mathematical foundation for related problems such as safe scheduling operation of a power system, efficient market clearing, reasonable planning and the like. The unit combination is essentially a combination optimization problem and is an important tool for the spot market clearance of the electric power market and the dispatching of the electric power system. With the continuous improvement of the duty ratio of new energy and the expansion of the scale of the interconnected large power grid, the power system expands the decision space in the coming time and space, and the constraint quantity and the constraint coupling degree of the unit combination problem are greatly increased, so that great challenges are brought to the solution of the unit combination problem, and a research model dimension reduction and constraint reduction method is needed urgently.
The research of the traditional unit combined model simplification method still has the defects, part of the model simplification methods are only suitable for the power grid mainly of the thermal power unit, the redundant constraint identification method tends to be homogeneous, and the new research of reducing the line safety constraint from the mathematical and physical perspectives is lacking.
Disclosure of Invention
The present application aims to solve, at least to some extent, one of the technical problems in the related art.
To achieve the above objective, an embodiment of a first aspect of the present application provides a method for optimizing and solving a safety constraint unit combination, including the following steps:
Acquiring a line safety constraint space based on a source power transfer distribution factor matrix, and measuring the similarity of line safety constraints in the line safety constraint space by calculating an included angle cosine value matrix;
acquiring a line safety constraint similarity matrix according to the similarity degree of the line safety constraint and a line similarity threshold, and acting on the similarity matrix through an integer programming model to determine a specific line aggregation mode so as to obtain an aggregation line safety constraint set;
acquiring economic operation basic data of an electric power system, and constructing a safety constraint unit combination model applicable to line safety constraint aggregation according to the aggregation line safety constraint set;
and solving the safety constraint unit combination model by adopting a mixed integer linear programming solver to obtain a line safety constraint aggregation unit combination result.
Optionally, the obtaining the line safety constraint space based on the source power transfer distribution factor matrix includes:
calculating to obtain a power transfer distribution factor matrix according to a relation equation between the node active power and the voltage phase angle;
according to the power transfer distribution factor matrix, combining nodes where power sources in an actual power grid are located to obtain the source power distribution transfer factor matrix;
Taking each row vector in the source power distribution transfer factor matrix as a point in Euclidean space, and establishing a line safety constraint space;
wherein the line safety constraint space shares N G Dimension, N G For all power supply amounts in the network.
Optionally, the measuring the similarity degree of the line safety constraint in the line safety constraint space by calculating an included angle cosine value matrix includes:
calculating an included angle cosine value matrix of all the line safety constraint vectors based on Euclidean distance and two norms in the line safety constraint space;
and judging whether different line safety constraint vectors are in the same direction according to a preset symbol indication matrix.
Optionally, the obtaining the line safety constraint similarity matrix according to the similarity degree of the line safety constraint and the line similarity threshold value includes:
setting the line safety constraint similarity matrix according to the difference value of the absolute value of the cosine value of the included angle between the line similarity threshold and the two line safety constraint vectors, and formulating as:
wherein, l andfor two different line safety constraint vectors, < ->For two line safety constraint vectors l and +.>Cosine value of included angle,/">For two line safety constraint vectors l and +. >Is a line safety constraint similarity matrix of (a) T For the set line similarity threshold value, N L Is the number of lines.
Optionally, the acting on the similarity matrix through an integer programming model to determine a specific line aggregation mode to obtain an aggregate line security constraint set includes:
determining the number of line safety constraints which maximize aggregation through an objective function, wherein the expression of the objective function is as follows:
in the formula, v l For line safety constraintsVector l indicates whether or not to participate in line aggregation, v l A line security constraint vector l is represented by 1 and is taken as a main aggregation line to participate in line aggregation;
and constraining the similar lines, if one line is used as a main aggregation line, all the similar lines participate in the aggregation and are aggregated into one line, wherein the constraint equation is as follows:
wherein,constraint vector +.>An indicator variable of whether to participate in line aggregation;
solving the objective function and the constraint equation to obtain a main aggregation line set participating in line aggregation;
acquiring a corresponding matched aggregation line set for the main aggregation lines in the main aggregation line set;
and generating an aggregation line safety constraint set according to the main aggregation line set and the aggregation line set.
Optionally, the power system economic operation basic data includes:
the power output upper limit and the power output lower limit of the unit, the climbing rate upper limit of the unit increase and decrease, the minimum continuous start-stop time of the unit, the running cost function of the unit, the start-stop cost function of the unit, the tide transfer distribution factor, the line transmission capacity and the positive and negative reserve rate of the system.
Optionally, the constructing a security constraint unit combination model applicable to line security constraint aggregation according to the aggregated line security constraint set includes:
determining a safety constraint unit combination objective function applicable to line safety constraint aggregation, wherein the expression is as follows:
wherein, the combined objective function of the safety restraint unit is used for representing the total running cost and start-stop cost of each variable-scale period to be minimized, and p i，t 、u i，t The output variable and the start-stop 0-1 variable of the generator set i in the period t are respectively, u i，t When the value is 0, the machine is stopped, 1 is started,for the start-stop cost function of the unit i, u i，0 Defining an initial start-stop state of the generator set i;
based on the safety constraint unit combination objective function, constructing a safety constraint unit combination constraint condition based on line aggregation constraint, wherein the constraint condition comprises a system power balance constraint, a system positive and negative standby constraint, a system positive and negative climbing standby constraint, a generator unit positive and negative climbing constraint, a generator unit minimum continuous start-stop time constraint, an aggregated line safety constraint and an un-aggregated line safety constraint;
And integrating the safety constraint unit combination constraint conditions, and establishing the safety constraint unit combination model.
To achieve the above object, an embodiment of a second aspect of the present application provides a safety restraint unit combination optimization solving device, including:
the similarity measurement module is used for acquiring a line safety constraint space based on the source power transfer distribution factor matrix and measuring the similarity of line safety constraints in the line safety constraint space by calculating an included angle cosine value matrix;
the safety constraint set acquisition module is used for acquiring a circuit safety constraint similarity matrix by combining the similarity degree of the circuit safety constraint with a circuit similarity threshold value, and acting on the similarity matrix through an integer programming model to determine a specific circuit aggregation mode so as to obtain an aggregated circuit safety constraint set;
the combined model construction module is used for acquiring the basic data of the economic operation of the power system and constructing a safety constraint unit combined model suitable for the line safety constraint aggregation according to the aggregation line safety constraint set;
and the solving module is used for solving the safety constraint unit combination model by adopting a mixed integer linear programming solver to obtain a line safety constraint aggregation unit combination result.
Optionally, the similarity measurement module is further configured to:
calculating to obtain a power transfer distribution factor matrix according to a relation equation between the node active power and the voltage phase angle;
according to the power transfer distribution factor matrix, combining nodes where power sources in an actual power grid are located to obtain the source power distribution transfer factor matrix;
taking each row vector in the source power distribution transfer factor matrix as a point in Euclidean space, and establishing a line safety constraint space;
wherein the line safety constraint space shares N G Dimension, N G For all power supply amounts in the network.
Optionally, the similarity measurement module is further configured to:
calculating an included angle cosine value matrix of all the line safety constraint vectors based on Euclidean distance and two norms in the line safety constraint space;
and judging whether different line safety constraint vectors are in the same direction according to a preset symbol indication matrix.
Optionally, the security constraint set acquisition module is further configured to:
setting the line safety constraint similarity matrix according to the difference value of the absolute value of the cosine value of the included angle between the line similarity threshold and the two line safety constraint vectors, and formulating as:
Wherein, l andfor two different line safety constraint vectors, < ->For two line safety constraint vectors l and +.>Cosine value of included angle,/">For two line safety constraint vectors l and +.>Is a line safety constraint similarity matrix of (a) T For the set line similarity threshold value, N L Is the number of lines.
Optionally, the security constraint set acquisition module is further configured to:
determining the number of line safety constraints which maximize aggregation through an objective function, wherein the expression of the objective function is as follows:
in the formula, v l Indicating variable, v, for whether line security constraint vector l participates in line aggregation l A line security constraint vector l is represented by 1 and is taken as a main aggregation line to participate in line aggregation;
and constraining the similar lines, if one line is used as a main aggregation line, all the similar lines participate in the aggregation and are aggregated into one line, wherein the constraint equation is as follows:
wherein,constraint vector +.>An indicator variable of whether to participate in line aggregation;
solving the objective function and the constraint equation to obtain a main aggregation line set participating in line aggregation;
acquiring a corresponding matched aggregation line set for the main aggregation lines in the main aggregation line set;
And generating an aggregation line safety constraint set according to the main aggregation line set and the aggregation line set.
Optionally, the power system economic operation basic data includes:
at least one of the upper limit and the lower limit of the output of the unit, the upper limit of the increasing and decreasing climbing rate of the unit, the minimum continuous start-stop time of the unit, the running cost function of the unit, the start-stop cost function of the unit, the power flow transfer distribution factor, the transmission capacity of a line and the positive and negative reserve rate of a system.
Optionally, the combined model building module is further configured to:
determining a safety constraint unit combination objective function applicable to line safety constraint aggregation, wherein the expression is as follows:
wherein, the combined objective function of the safety restraint unit is used for representing the total running cost and start-stop cost of each variable-scale period to be minimized, and p i，t 、u i，t The output variable and the start-stop 0-1 variable of the generator set i in the period t are respectively, u i，t When the value is 0, the machine is stopped, 1 is started,for the start-stop cost function of the unit i, u i，0 Defined as generatorAn initial start-stop state for group i;
based on the safety constraint unit combination objective function, constructing a safety constraint unit combination constraint condition based on line aggregation constraint, wherein the constraint condition comprises a system power balance constraint, a system positive and negative standby constraint, a system positive and negative climbing standby constraint, a generator unit positive and negative climbing constraint, a generator unit minimum continuous start-stop time constraint, an aggregated line safety constraint and an un-aggregated line safety constraint;
And integrating the safety constraint unit combination constraint conditions, and establishing the safety constraint unit combination model.
To achieve the above object, an embodiment of a third aspect of the present application provides an electronic device, including: a processor, and a memory communicatively coupled to the processor;
the memory stores computer-executable instructions;
the processor executes computer-executable instructions stored by the memory to implement the method of any one of the first aspects above.
To achieve the above object, an embodiment of a fourth aspect of the present application proposes a computer-readable storage medium having stored therein computer-executable instructions for implementing the method according to any of the above first aspects when being executed by a processor.
To achieve the above object, an embodiment of a fifth aspect of the present application proposes a computer program product comprising a computer program which, when executed by a processor, implements a method as described in any of the above first aspects.
The technical scheme provided by the embodiment of the application at least brings the following beneficial effects:
the method has the advantages that the line safety constraint in the unit combination optimization problem is reduced based on line aggregation, the number of the unit combination line safety constraints is greatly reduced, the complexity of a model is reduced, meanwhile, the accuracy and the economy of a unit combination result are basically guaranteed, the optimization solving efficiency of the power system based on the safety constraint unit combination is improved, the optimized operation of the power system after a large number of introduced new energy sources, energy storage, power electronic equipment and the like are supported, and the stable operation of the power system is guaranteed.
Additional aspects and advantages of the application will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the application.
Drawings
The foregoing and/or additional aspects and advantages of the present application will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings, in which:
FIG. 1 is a flow chart illustrating a method of combined optimization solution for a safety restraint unit according to an embodiment of the present application;
FIG. 2 is a flow chart illustrating a method of combined optimization solution for a safety restraint set according to another embodiment of the present application;
FIG. 3 is a flow chart illustrating a method of combined optimization solution for a safety restraint set according to another embodiment of the present application;
FIG. 4 is a flow chart illustrating a method of combined optimization solution for a safety restraint set according to another embodiment of the present application;
FIG. 5 is a block diagram of a security constraint unit combination optimization solver according to an embodiment of the present application;
fig. 6 is a block diagram of an electronic device.
Detailed Description
Embodiments of the present application are described in detail below, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to like or similar elements or elements having like or similar functions throughout. The embodiments described below by referring to the drawings are exemplary and intended for the purpose of explaining the present application and are not to be construed as limiting the present application.
The following describes a method and a device for solving the combination optimization of a safety constraint unit in an embodiment of the application with reference to the accompanying drawings.
A method for optimizing and solving a safety constraint unit set according to an embodiment of the application is described with reference to the accompanying drawings.
FIG. 1 is a flow chart illustrating a method of combined optimization solution for a security restraint set according to an embodiment of the present application.
As shown in fig. 1, the method comprises the steps of:
step 101, obtaining a line safety constraint space based on a source power transfer distribution factor matrix, and measuring the similarity of the line safety constraint in the line safety constraint space by calculating an included angle cosine value matrix.
In the embodiment of the application, the similarity degree of the circuit safety constraint in the circuit safety constraint space needs to be calculated for the reduction of the subsequent circuit safety constraint.
Specifically, as shown in fig. 2, step 101 further includes:
step 201, calculating to obtain a power transfer distribution factor matrix according to a relation equation between the node active power and the voltage phase angle.
In the embodiment of the application, a relation equation between node active power and voltage phase angle is given first:
P n ＝B x θ
wherein B is x For the node admittance matrix calculated based on kirchhoff's law, the relation of node current and voltage, X, is described i，j Is the line reactance between node i and node j.
And B is x Can be represented by the following expression:
thus, the relationship between line flow and voltage phase angle can be described by the following equation:
P l ＝B l θ
wherein B is l For the line-node admittance matrix, the relationship between line current and node voltage difference is described.
In one possible embodiment, assuming that a line connection relationship exists between node 1 and node 2, and between node 1 and node 3, then B line Can be represented as (rows represent lines,columns represent line two end nodes):
thus, by combining the formulas set forth above, one can obtain:
P l ＝B l (B x ) -1 P n
G＝B l (B x ) -1
wherein G is a power transfer distribution factor matrix.
Note that B x Is a singular matrix, can not directly invert B x Setting the corresponding item of the middle reference node to zero, thus obtainingTherefore, the power transfer distribution factor matrix G is obtained through the following equation in actual calculation:
step 202, according to the power transfer distribution factor matrix, combining the nodes where the power supplies in the actual power grid are located to obtain a source power distribution transfer factor matrix.
In the embodiment of the present application, as can be seen from the formula shown in step 201, the power transfer distribution factor matrix G is a matrix of size N L ×N N Wherein N is a real number matrix of L N is the number of network lines N Is the number of network nodes.
Therefore, based on the power transfer distribution factor matrix and the node where the power supply is located in the actual power grid, the following source power distribution transfer factor matrix T can be obtained:
wherein N is c For networksAll power supply numbers, k n Numbering the nodes where the nth unit is located.
In step 203, each row vector in the source power distribution transfer factor matrix is used as a point in the euclidean space, and a line safety constraint space is established.
In the embodiment of the application, N can be established G The line of dimensions safely constrains the space.
And 204, calculating an included angle cosine value matrix of all the line safety constraint vectors based on Euclidean distance and two norms in the line safety constraint space.
It will be appreciated that the closer the absolute value of the cosine value of the angle between the two line safety constraint vectors is to 1, the more similar the two line safety constraint vectors are.
In the embodiment of the present application, in the line safety constraint space established in the foregoing, the cosine value matrix a of the included angle of all the line safety constraint vectors is directly calculated based on the euclidean distance and the two norms:
in the method, in the process of the invention,for two line safety constraint vectors l and +.>Cosine of the included angle.
It is understood that, in step 101, since the directions of different line safety constraint vectors may be different, the present application considers the following symbol indication matrix S to indicate whether a certain line safety constraint vector is in the same direction as all other line safety constraint vectors:
wherein,let 1 denote the line safety constraint vectors l and +.>In the same direction, otherwise, different directions are indicated.
Step 102, obtaining a line safety constraint similarity matrix according to the similarity degree of the line safety constraint and a line similarity threshold, and acting on the similarity matrix through an integer programming model to determine a specific line aggregation mode so as to obtain an aggregation line safety constraint set.
Specifically, as shown in fig. 3, step 102 further includes:
step 301, setting a line safety constraint similarity matrix according to the difference between the line similarity threshold and the absolute value of the cosine value of the included angle of the two line safety constraint vectors.
In this embodiment of the present application, the cosine value matrix a of the included angle obtained above measures the similarity of the security constraints of each line, and in this step, a line similarity threshold α is set T The line safety constraint similarity matrix can thus be obtained as follows:
wherein B is l，l` =1 represents that line safety constraint l is similar to line safety constraint l', and line aggregation can be performed.
Step 302, determining, by an objective function, a number of line security constraints that maximize participation in the aggregation.
In the embodiment of the present application, in the case that repeated similarity may occur in these lines, in order to maximally aggregate the lines and avoid the occurrence of repeated aggregation of the lines as much as possible, it may be determined which lines specifically participate in aggregation based on the following integer programming model, where the objective function expression is as follows:
wherein v is l Indicating variable, v, for whether line security constraint vector l participates in line aggregation l And 1 represents a line safety constraint vector l as a main aggregation line to participate in line aggregation.
In step 303, constraint is performed on similar lines, and if one line is used as the main aggregation line, all similar lines participate in the current aggregation and are aggregated into one line.
In the embodiment of the application, two similar line safety constraints are limited by the constraint formula, and the two similar line safety constraints cannot be simultaneously used as a main aggregation line to participate in line safety constraint aggregation, wherein the constraint formula is as follows:
in the method, in the process of the invention,constraint vector +.>An indicator variable of whether to participate in line aggregation.
And step 304, solving the objective function and the constraint equation to obtain a main aggregation line set participating in line aggregation.
In the embodiment of the application, solving the above model to obtain a main aggregation line set participating in line aggregation:
wherein N is mc Is the number of main aggregate lines.
Step 305, for the main aggregation lines in the main aggregation line set, obtaining a corresponding matched aggregation line set.
In the embodiment of the present application, for each main aggregation line, the main aggregation lines are concentratedAll have an aggregate line set matching it +.>The properties are defined as follows
Taking all the aggregation line sets as a union set to obtain the following total aggregation line set:
likewise, the set of bus lines { l } pertains to the set of total aggregate lines { l } c Complement to get the set of lines not participating in aggregation { l } e }：
{l e }＝C {l} {l c }
Step 306, generating an aggregate line security constraint set according to the main aggregate line set and the aggregate line set.
In the embodiment of the application, a general unit combined line safety constraint based on direct current is considered:
wherein P is l，max For maximum power flow of line l, p k，t For the output of the kth unit at the time t, n k Is the node where it is located. d, d j，t For the j-th load at time t, n j Is the node where it is located.
It can be appreciated that the unit combination line safety constraint based on the direct current power flow can be converted into:
since for each main aggregation line Has an aggregate line set matched with it>Which in effect constitutes a series of line safety constraints
All the line safety constraint vectors are converted into vectors which are in the same direction as the main aggregation line safety constraint vector based on the symbol indication matrix S, and the obtained states are:
and replaced as follows:
subsequently, in the aggregate line setRespectively find l ci min And->So that one of the following two sets of relations is established, and for all the aggregation lines, only one mode is selected to find l ci min And->
The first search method is as follows:
the second search mode:
thus, embodiments of the present application replace all line safety constraints with an aggregate line safety constraint set shown below.
The aggregate line security constraint set is:
all line safety constraints are:
and 103, acquiring economic operation basic data of the electric power system, and constructing a safety constraint unit combination model suitable for line safety constraint aggregation according to the aggregation line safety constraint set.
Specifically, as shown in fig. 4, step 103 further includes:
step 401, obtaining basic data of economic operation of the power system.
In the embodiment of the application, the economic operation basic data of the power system comprise an upper limit and a lower limit of unit output, an upper limit of increasing and decreasing climbing rate of the unit, minimum continuous start-stop time of the unit, a unit operation cost function, a unit start-stop cost function, a tide transfer distribution factor, line transmission capacity and positive and negative reserve rate of the system.
Step 402, determining a security constraint set combination objective function applicable to line security constraint aggregation.
In the embodiment of the application, the expression of the safety restraint unit combined objective function is as follows:
wherein, the combined objective function of the safety restraint unit is used for representing the total running cost and start-stop cost of each variable-scale period to be minimized, and p i，t 、u i，t The output variable and the start-stop 0-1 variable of the generator set i in the period t are respectively, u i，t When the value is 0, the machine is stopped, 1 is started,for the start-stop cost function of the unit i, u i，0 Defined as the initial start-stop state of genset i.
And step 403, constructing a safety constraint unit combination constraint condition based on the line aggregation constraint based on the safety constraint unit combination objective function.
In the embodiment of the application, constraint conditions comprise system power balance constraint, system positive and negative standby constraint, system positive and negative climbing standby constraint, generator set positive and negative climbing constraint, minimum continuous start-stop time constraint of the generator set, aggregated line safety constraint and line safety constraint which does not participate in aggregation.
Wherein, the system power balance constraint is:
wherein D is t Is the system payload during period t.
The positive and negative standby constraint of the system is as follows:
wherein r is + 、r - For the positive and negative standby rate of the system, t Dmaximum and minimum, respectively, of the system payload during period t> i PThe upper and lower limits of the output of the generator set i are respectively set.
The standby constraint of positive and negative climbing of the system is as follows:
and the following represents the generator set output range constraint
RU in i 、RD i The upper limits of the output increase and decrease are respectively given to the generator set i.
The positive and negative climbing constraint of the generator set is as follows:
in SU i The upper limit of the output of the generator set i during starting is set; SD (secure digital memory card) i The upper limit of the output before the generator set i is turned off.
The minimum continuous start-stop time constraint of the generator set is as follows:
in the YU i 、TD i The minimum continuous on-off time (in hours) of the generator set i is given, and T is the number of time periods of the set combination.
It will be appreciated that the positive and negative hill climbing constraints of the generator set and the minimum continuous on-off time constraints of the generator set are constraints specific to the thermal power plant, and that for other possible power sources, such as hydroelectric power plants, wind power plants, etc., it is also necessary to meet their respective physical characteristic constraints.
The aggregated line security constraints are:
in { l } c i The main aggregate line set, l ci max Andif the first searching mode or the second searching mode is met, a looser aggregation line safety constraint is obtained, the corresponding unit combination model solving acceleration is more remarkable, but the error of the obtained feasible solution is larger; if the second searching mode is met, a more compact aggregation line safety constraint is obtained, the obtained error of the feasible solution is smaller, an infeasible model is possibly obtained, and the solving speed is not as high as that of the former model.
The line security constraints not participating in aggregation are:
in the formula, { l e And is a set of lines that do not participate in the aggregation.
And step 404, integrating the safety restraint unit combination constraint conditions, and establishing a safety restraint unit combination model.
In the embodiment of the application, a safety restraint unit combination model is constructed according to various restraint conditions mentioned in the application.
And 104, solving the safety constraint unit combination model by adopting a mixed integer linear programming solver to obtain a line safety constraint aggregation unit combination result.
In the embodiment of the present application, the combined result of the line safety constraint aggregate unit is obtained by solving the safety constraint unit combination model by using a combined integer linear programming solver related to the integer programming model decision in step 102.
According to the safety constraint unit combination optimization solving method, the reduction of the line safety constraint in the unit combination optimization problem is realized based on line aggregation, the number of the unit combination line safety constraints is greatly reduced, the complexity of a model is reduced, meanwhile, the accuracy and economy of a unit combination result are basically guaranteed, the optimization solving efficiency of an electric power system based on the safety constraint unit combination is improved, the optimized operation of a large number of introduced electric power systems such as new energy, energy storage and power electronic equipment is supported, and the stable operation of the electric power system is guaranteed.
Next, a safety constraint unit combination optimization solving device according to an embodiment of the application is described with reference to the accompanying drawings.
FIG. 5 is a block diagram illustrating a security constraint set combined optimization solver according to an embodiment of the present application.
As shown in fig. 5, the block diagram of the security constraint set-up optimization solver 10 includes a similarity measurement module 100, a security constraint set acquisition module 200, a combined model construction module 300, and a solution module 400.
The similarity measurement module 100 is configured to obtain a line safety constraint space based on the source power transfer distribution factor matrix, and measure a similarity of line safety constraints in the line safety constraint space by calculating an angle cosine value matrix
The safety constraint set acquisition module 200 is configured to acquire a line safety constraint similarity matrix by combining a similarity degree of line safety constraints and a line similarity threshold, and act on the similarity matrix through an integer programming model to determine a specific line aggregation mode, so as to obtain an aggregate line safety constraint set;
the combined model construction module 300 is used for acquiring the basic data of the economic operation of the power system and constructing a safety constraint unit combined model suitable for the line safety constraint aggregation according to the aggregation line safety constraint set;
And the solving module 400 is used for solving the safety constraint unit combination model by adopting a mixed integer linear programming solver to obtain a line safety constraint aggregation unit combination result.
Further, in one embodiment of the present application, the similarity measurement module 100 is configured to:
calculating to obtain a power transfer distribution factor matrix according to a relation equation between the node active power and the voltage phase angle;
according to the power transfer distribution factor matrix, combining nodes where power sources in an actual power grid are located to obtain a source power distribution transfer factor matrix;
taking each row vector in the source power distribution transfer factor matrix as a point in Euclidean space, and establishing a line safety constraint space;
wherein the line safety constraint space shares N G Dimension, N c For all power supply quantities in the network;
calculating an included angle cosine value matrix of all the line safety constraint vectors based on Euclidean distance and two norms in the line safety constraint space;
and judging whether different line safety constraint vectors are in the same direction according to a preset symbol indication matrix.
Further, in one embodiment of the present application, the security constraint set acquisition module 200 is configured to:
setting a line safety constraint similarity matrix according to the difference value between the line similarity threshold and the absolute value of the cosine value of the included angle of the two line safety constraint vectors, and formulating as:
Wherein, l andfor two different line safety constraint vectors, < ->For two line safety constraint vectors l and +.>Cosine value of included angle,/">For two line safety constraint vectors l and +.>Is a line safety constraint similarity matrix of (a) T For the set line similarity threshold value, N L Is the number of lines.
In one embodiment of the present application, the security constraint set acquisition module 200 is further configured to:
determining the number of line safety constraints which maximize aggregation through an objective function, wherein the expression of the objective function is as follows:
/>
in the formula, v l Indicating variable, v, for whether line security constraint vector l participates in line aggregation l A line security constraint vector l is represented by 1 and is taken as a main aggregation line to participate in line aggregation;
and constraining the similar lines, if one line is used as a main aggregation line, all the similar lines participate in the aggregation and are aggregated into one line, wherein the constraint equation is as follows:
wherein,constraint vector +.>An indicator variable of whether to participate in line aggregation;
solving the objective function and the constraint equation to obtain a main aggregation line set participating in line aggregation;
acquiring a corresponding matched aggregation line set for a main aggregation line in the main aggregation line set;
And generating an aggregation line security constraint set according to the main aggregation line set and the aggregation line set.
Optionally, the power system economic operation basic data includes:
at least one of the upper limit and the lower limit of the output of the unit, the upper limit of the increasing and decreasing climbing rate of the unit, the minimum continuous start-stop time of the unit, the running cost function of the unit, the start-stop cost function of the unit, the power flow transfer distribution factor, the transmission capacity of a line and the positive and negative reserve rate of a system.
Further, in one embodiment of the present application, the combined model construction module 300 is configured to:
determining a safety constraint unit combination objective function applicable to line safety constraint aggregation, wherein the expression is as follows:
wherein, the combined objective function of the safety restraint unit is used for representing the total running cost and start-stop cost of each variable-scale period to be minimized, and p i，t 、u i，t The output variable and the start-stop 0-1 variable of the generator set i in the period t are respectively, u i，t When the value is 0, the machine is stopped, 1 is started,for the start-stop cost function of the unit i, u i，0 Defining an initial start-stop state of the generator set i;
based on a safety constraint unit combination objective function, constructing a safety constraint unit combination constraint condition based on line aggregation constraint, wherein the constraint condition comprises a system power balance constraint, a system positive and negative standby constraint, a system positive and negative climbing standby constraint, a generator unit positive and negative climbing constraint, a generator unit minimum continuous start-stop time constraint, an aggregated line safety constraint and an un-aggregated line safety constraint;
And integrating the safety restraint unit combination constraint conditions, and establishing a safety restraint unit combination model.
The specific manner in which the various modules perform the operations in the apparatus of the above embodiments have been described in detail in connection with the embodiments of the method, and will not be described in detail herein.
According to the safety constraint unit combination optimization solving device, line safety constraint is reduced in the unit combination optimization problem based on line aggregation, the number of the unit combination line safety constraints is greatly reduced, the complexity of a model is reduced, meanwhile, accuracy and economy of unit combination results are basically guaranteed, the optimization solving efficiency of an electric power system based on the safety constraint unit combination is improved, the optimized operation of a large number of introduced electric power systems such as new energy, energy storage and power electronic equipment is supported, and the stable operation of the electric power system is guaranteed.
Fig. 6 shows a schematic block diagram of an example electronic device 700 that may be used to implement embodiments of the present application. Electronic devices are intended to represent various forms of digital computers, such as laptops, desktops, workstations, personal digital assistants, servers, blade servers, mainframes, and other appropriate computers. The electronic device may also represent various forms of mobile devices, such as personal digital processing, cellular telephones, smartphones, wearable devices, and other similar computing devices. The components shown herein, their connections and relationships, and their functions, are meant to be exemplary only, and are not meant to limit implementations of the application described and/or claimed herein.
As shown in fig. 6, the apparatus 700 includes a computing unit 701 that can perform various appropriate actions and processes according to a computer program stored in a Read Only Memory (ROM) 702 or a computer program loaded from a storage unit 708 into a Random Access Memory (RAM) 703. In the RAM 703, various programs and data required for the operation of the device 700 may also be stored. The computing unit 701, the ROM 702, and the RAM 703 are connected to each other through a bus 704. An input/output (I/O) interface 705 is also connected to bus 704.
Various components in device 700 are connected to I/O interface 705, including: an input unit 706 such as a keyboard, a mouse, etc.; an output unit 707 such as various types of displays, speakers, and the like; a storage unit 708 such as a magnetic disk, an optical disk, or the like; and a communication unit 709 such as a network card, modem, wireless communication transceiver, etc. The communication unit 709 allows the device 700 to exchange information/data with other devices via a computer network, such as the internet, and/or various telecommunication networks.
The computing unit 701 may be a variety of general and/or special purpose processing components having processing and computing capabilities. Some examples of computing unit 701 include, but are not limited to, a Central Processing Unit (CPU), a Graphics Processing Unit (GPU), various specialized artificial intelligence (AT) computing chips, various computing units running machine learning model algorithms, a Digital Signal Processor (DSP), and any suitable processor, controller, microcontroller, etc. The computing unit 701 performs the respective methods and processes described above, such as a voice instruction response method. For example, in some embodiments, the voice instruction response method may be implemented as a computer software program tangibly embodied on a machine-readable medium, such as storage unit 708. In some embodiments, part or all of the computer program may be loaded and/or installed onto device 700 via ROM 702 and/or communication unit 709. When the computer program is loaded into RAM 703 and executed by computing unit 701, one or more steps of the voice instruction response method described above may be performed. Alternatively, in other embodiments, the computing unit 701 may be configured to perform the voice instruction response method by any other suitable means (e.g., by means of firmware).
Various implementations of the systems and techniques described here above may be implemented in digital electronic circuitry, integrated circuit systems, field Programmable Gate Arrays (FPGAs), application Specific Integrated Circuits (ASICs), application Specific Standard Products (ASSPs), systems On Chip (SOCs), load programmable logic devices (CPLDs), computer hardware, firmware, software, and/or combinations thereof. These various embodiments may include: implemented in one or more computer programs, the one or more computer programs may be executed and/or interpreted on a programmable system including at least one programmable processor, which may be a special purpose or general-purpose programmable processor, that may receive data and instructions from, and transmit data and instructions to, a storage system, at least one input device, and at least one output device.
Program code for carrying out methods of the present application may be written in any combination of one or more programming languages. These program code may be provided to a processor or controller of a general purpose computer, special purpose computer, or other programmable data processing apparatus such that the program code, when executed by the processor or controller, causes the functions/operations specified in the flowchart and/or block diagram to be implemented. The program code may execute entirely on the machine, partly on the machine, as a stand-alone software package, partly on the machine and partly on a remote machine or entirely on the remote machine or server.
In the context of this application, a machine-readable medium may be a tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device. The machine-readable medium may be a machine-readable signal medium or a machine-readable storage medium. The machine-readable medium may include, but is not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any suitable combination of the foregoing. More specific examples of a machine-readable storage medium would include an electrical connection based on one or more wires, a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.
To provide for interaction with a user, the systems and techniques described here can be implemented on a computer having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to a user; and a keyboard and pointing device (e.g., a mouse or trackball) by which a user can provide input to the computer. Other kinds of devices may also be used to provide for interaction with a user; for example, feedback provided to the user may be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user may be received in any form, including acoustic input, speech input, or tactile input.
The systems and techniques described here can be implemented in a computing system that includes a background component (e.g., as a data server), or that includes a middleware component (e.g., an application server), or that includes a front-end component (e.g., a user computer having a graphical user interface or a web browser through which a user can interact with an implementation of the systems and techniques described here), or any combination of such background, middleware, or front-end components. The components of the system can be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: local Area Networks (LANs), wide Area Networks (WANs), the internet, and blockchain networks.
The computer system may include a client and a server. The client and server are typically remote from each other and typically interact through a communication network. The relationship of client and server arises by virtue of computer programs running on the respective computers and having a client-server relationship to each other. The server can be a cloud server, also called a cloud computing server or a cloud host, and is a host product in a cloud computing service system, so that the defects of high management difficulty and weak service expansibility in the traditional physical hosts and VPS service ("Virtual Private Server" or simply "VPS") are overcome. The server may also be a server of a distributed system or a server that incorporates a blockchain.
It should be appreciated that various forms of the flows shown above may be used to reorder, add, or delete steps. For example, the steps described in the present application may be performed in parallel, sequentially, or in a different order, so long as the desired results of the technical solutions of the present application are achieved, and the present application is not limited herein.
The above embodiments do not limit the scope of the application. It will be apparent to those skilled in the art that various modifications, combinations, sub-combinations and alternatives are possible, depending on design requirements and other factors. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present application are intended to be included within the scope of the present application
It should be appreciated that various forms of the flows shown above may be used to reorder, add, or delete steps. For example, the steps described in the present application may be performed in parallel, sequentially, or in a different order, so long as the desired results of the technical solutions of the present application are achieved, and the present application is not limited herein.
The above embodiments do not limit the scope of the application. It will be apparent to those skilled in the art that various modifications, combinations, sub-combinations and alternatives are possible, depending on design requirements and other factors. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present application are intended to be included within the scope of the present application.
Claims (17)
1. The combined optimization solving method of the safety constraint unit is characterized by comprising the following steps of:
acquiring a line safety constraint space based on a source power transfer distribution factor matrix, and measuring the similarity of line safety constraints in the line safety constraint space by calculating an included angle cosine value matrix;
acquiring a line safety constraint similarity matrix according to the similarity degree of the line safety constraint and a line similarity threshold, and acting on the similarity matrix through an integer programming model to determine a specific line aggregation mode so as to obtain an aggregation line safety constraint set;
acquiring economic operation basic data of an electric power system, and constructing a safety constraint unit combination model applicable to line safety constraint aggregation according to the aggregation line safety constraint set;
and solving the safety constraint unit combination model by adopting a mixed integer linear programming solver to obtain a line safety constraint aggregation unit combination result.
2. The method of claim 1, wherein the obtaining the line safety constraint space based on the source power transfer distribution factor matrix comprises:
calculating to obtain a power transfer distribution factor matrix according to a relation equation between the node active power and the voltage phase angle;
According to the power transfer distribution factor matrix, combining nodes where power sources in an actual power grid are located to obtain the source power distribution transfer factor matrix;
taking each row vector in the source power distribution transfer factor matrix as a point in Euclidean space, and establishing a line safety constraint space;
wherein the line safety constraint space shares N G Dimension, N G For all power supply amounts in the network.
3. The method of claim 1, wherein said measuring the similarity of the line safety constraints in the line safety constraint space by calculating an angle cosine value matrix comprises:
calculating an included angle cosine value matrix of all the line safety constraint vectors based on Euclidean distance and two norms in the line safety constraint space;
and judging whether different line safety constraint vectors are in the same direction according to a preset symbol indication matrix.
4. The method of claim 1, wherein the obtaining a line security constraint similarity matrix from the line security constraint similarity to a line similarity threshold comprises:
setting the line safety constraint similarity matrix according to the difference value of the absolute value of the cosine value of the included angle between the line similarity threshold and the two line safety constraint vectors, and formulating as:
Wherein, l andfor two different line safety constraint vectors, < ->For two line safety constraint vectors l and +.>Cosine value of included angle,/">For two line safety constraint vectors l and +.>Is a line safety constraint similarity matrix of (a) T For the set line similarity threshold value, N L Is the number of lines.
5. The method of claim 4, wherein the acting on the similarity matrix by an integer programming model to determine a specific line aggregation manner, to obtain an aggregate line security constraint set, comprises:
determining the number of line safety constraints which maximize aggregation through an objective function, wherein the expression of the objective function is as follows:
in the formula, v l Indicating variable, v, for whether line security constraint vector l participates in line aggregation l A line security constraint vector l is represented by 1 and is taken as a main aggregation line to participate in line aggregation;
and constraining the similar lines, if one line is used as a main aggregation line, all the similar lines participate in the aggregation and are aggregated into one line, wherein the constraint equation is as follows:
wherein,constraint vector +.>An indicator variable of whether to participate in line aggregation;
solving the objective function and the constraint equation to obtain a main aggregation line set participating in line aggregation;
Acquiring a corresponding matched aggregation line set for the main aggregation lines in the main aggregation line set;
and generating an aggregation line safety constraint set according to the main aggregation line set and the aggregation line set.
6. The method of any one of claims 1-5, wherein the power system economic operation base data comprises:
at least one of the upper limit and the lower limit of the output of the unit, the upper limit of the increasing and decreasing climbing rate of the unit, the minimum continuous start-stop time of the unit, the running cost function of the unit, the start-stop cost function of the unit, the power flow transfer distribution factor, the transmission capacity of a line and the positive and negative reserve rate of a system.
7. The method of claim 1, wherein constructing a security constraint set combination model suitable for line security constraint aggregation from the aggregated line security constraint set comprises:
determining a safety constraint unit combination objective function applicable to line safety constraint aggregation, wherein the expression is as follows:
wherein, the combined objective function of the safety restraint unit is used for representing the total running cost and start-stop cost of each variable-scale period to be minimized, and p i，t 、u i，t The output variable and the start-stop 0-1 variable of the generator set i in the period t are respectively, u i，t When the value is 0, the machine is stopped, 1 is started,for the start-stop cost function of the unit i, u i，0 Defining an initial start-stop state of the generator set i;
based on the safety constraint unit combination objective function, constructing a safety constraint unit combination constraint condition based on line aggregation constraint, wherein the constraint condition comprises a system power balance constraint, a system positive and negative standby constraint, a system positive and negative climbing standby constraint, a generator unit positive and negative climbing constraint, a generator unit minimum continuous start-stop time constraint, an aggregated line safety constraint and an un-aggregated line safety constraint;
and integrating the safety constraint unit combination constraint conditions, and establishing the safety constraint unit combination model.
8. The utility model provides a safe restraint unit combination optimization solution device which characterized in that includes:
the similarity measurement module is used for acquiring a line safety constraint space based on the source power transfer distribution factor matrix and measuring the similarity of line safety constraints in the line safety constraint space by calculating an included angle cosine value matrix;
the safety constraint set acquisition module is used for acquiring a circuit safety constraint similarity matrix by combining the similarity degree of the circuit safety constraint with a circuit similarity threshold value, and acting on the similarity matrix through an integer programming model to determine a specific circuit aggregation mode so as to obtain an aggregated circuit safety constraint set;
The combined model construction module is used for acquiring the basic data of the economic operation of the power system and constructing a safety constraint unit combined model suitable for the line safety constraint aggregation according to the aggregation line safety constraint set;
and the solving module is used for solving the safety constraint unit combination model by adopting a mixed integer linear programming solver to obtain a line safety constraint aggregation unit combination result.
9. The apparatus of claim 8, wherein the similarity measurement module is further configured to:
calculating to obtain a power transfer distribution factor matrix according to a relation equation between the node active power and the voltage phase angle;
according to the power transfer distribution factor matrix, combining nodes where power sources in an actual power grid are located to obtain the source power distribution transfer factor matrix;
taking each row vector in the source power distribution transfer factor matrix as a point in Euclidean space, and establishing a line safety constraint space;
wherein the line safety constraint space shares N G Dimension, N G For all power supply amounts in the network.
10. The apparatus of claim 8, wherein the similarity measurement module is further configured to:
Calculating an included angle cosine value matrix of all the line safety constraint vectors based on Euclidean distance and two norms in the line safety constraint space;
and judging whether different line safety constraint vectors are in the same direction according to a preset symbol indication matrix.
11. The apparatus of claim 8, wherein the security constraint set acquisition module is further to:
setting the line safety constraint similarity matrix according to the difference value of the absolute value of the cosine value of the included angle between the line similarity threshold and the two line safety constraint vectors, and formulating as:
wherein, l andfor two different line safety constraint vectors, < ->For two line safety constraint vectors l and +.>Cosine value of included angle,/">For two line safety constraint vectors l and +.>Is a line safety constraint similarity matrix of (a) T For the set line similarity threshold value, N L Is the number of lines.
12. The apparatus of claim 11, wherein the security constraint set acquisition module is further configured to:
determining the number of line safety constraints which maximize aggregation through an objective function, wherein the expression of the objective function is as follows:
in the formula, v l Indicating variable, v, for whether line security constraint vector l participates in line aggregation l A line security constraint vector l is represented by 1 and is taken as a main aggregation line to participate in line aggregation;
and constraining the similar lines, if one line is used as a main aggregation line, all the similar lines participate in the aggregation and are aggregated into one line, wherein the constraint equation is as follows:
wherein,constraint vector +.>An indicator variable of whether to participate in line aggregation;
solving the objective function and the constraint equation to obtain a main aggregation line set participating in line aggregation;
acquiring a corresponding matched aggregation line set for the main aggregation lines in the main aggregation line set;
and generating an aggregation line safety constraint set according to the main aggregation line set and the aggregation line set.
13. The apparatus of any one of claims 8-12, wherein the power system economic operation base data comprises:
at least one of the upper limit and the lower limit of the output of the unit, the upper limit of the increasing and decreasing climbing rate of the unit, the minimum continuous start-stop time of the unit, the running cost function of the unit, the start-stop cost function of the unit, the power flow transfer distribution factor, the transmission capacity of a line and the positive and negative reserve rate of a system.
14. The apparatus of claim 8, wherein the combined model construction module is further configured to:
Determining a safety constraint unit combination objective function applicable to line safety constraint aggregation, wherein the expression is as follows:
wherein, the combined objective function of the safety restraint unit is used for representing the total running cost and start-stop cost of each variable-scale period to be minimized, and p i，t 、u i，t The output variable and the start-stop 0-1 variable of the generator set i in the period t are respectively, u i，t When the value is 0, the machine is stopped, 1 is started,for the start-stop cost function of the unit i, u i，0 Defining an initial start-stop state of the generator set i;
based on the safety constraint unit combination objective function, constructing a safety constraint unit combination constraint condition based on line aggregation constraint, wherein the constraint condition comprises a system power balance constraint, a system positive and negative standby constraint, a system positive and negative climbing standby constraint, a generator unit positive and negative climbing constraint, a generator unit minimum continuous start-stop time constraint, an aggregated line safety constraint and an un-aggregated line safety constraint;
and integrating the safety constraint unit combination constraint conditions, and establishing the safety constraint unit combination model.
15. An electronic device, comprising: a processor, and a memory communicatively coupled to the processor;
The memory stores computer-executable instructions;
the processor executes computer-executable instructions stored in the memory to implement the method of any one of claims 1-7.
16. A computer readable storage medium having stored therein computer executable instructions which when executed by a processor are adapted to carry out the method of any one of claims 1-7.
17. A computer program product comprising a computer program which, when executed by a processor, implements the method of any of claims 1-7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311406837.1A CN117411082A (en) | 2023-10-26 | 2023-10-26 | Combined optimization solving method and device for safety constraint unit |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311406837.1A CN117411082A (en) | 2023-10-26 | 2023-10-26 | Combined optimization solving method and device for safety constraint unit |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117411082A true CN117411082A (en) | 2024-01-16 |
Family
ID=89486750
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311406837.1A Pending CN117411082A (en) | 2023-10-26 | 2023-10-26 | Combined optimization solving method and device for safety constraint unit |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117411082A (en) |
-
2023
- 2023-10-26 CN CN202311406837.1A patent/CN117411082A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN117332897A (en) | Artificial intelligent driven new energy small time scale power interpolation integrated prediction method | |
| CN117411082A (en) | Combined optimization solving method and device for safety constraint unit | |
| CN116245253A (en) | Substation planning scheme determining method and device, electronic equipment and storage medium | |
| CN117131315B (en) | Out-of-tolerance electric energy meter determining method and medium based on solving multi-element quadratic function extremum | |
| CN117131353B (en) | Method and device for determining out-of-tolerance electric energy meter, electronic equipment and storage medium | |
| CN117200266B (en) | Distribution line photovoltaic collection point energy storage capacity configuration method, device and equipment | |
| Shi et al. | Wind Power Prediction Based On Rough Set Conditional Entropy Weighted Markov Chain | |
| CN117094537B (en) | Power grid planning method, device, electronic equipment and storage medium | |
| CN116014759A (en) | New energy-oriented multistage self-adaptive edge intelligent control method and device | |
| CN116094059A (en) | Layout optimization method, device and equipment for distributed power supply and storage medium | |
| CN115511047B (en) | Quantification method, device, equipment and medium of Softmax model | |
| CN116128204A (en) | Power distribution network scheduling method and device, electronic equipment and storage medium | |
| CN117767353A (en) | Power support cooperative control method and device for virtual synchronous power distribution network | |
| CN116404650A (en) | Power grid power supply area dividing method, device, equipment and storage medium | |
| CN116260392A (en) | Method, device, equipment, storage medium and program product for obtaining direct current line loss of photovoltaic power station | |
| CN116365523A (en) | Tidal current calculation method, device, equipment and medium applied to power distribution network | |
| CN116054135A (en) | Reactive power prediction method and device for micro-grid, electronic equipment and storage medium | |
| CN117422238A (en) | Method and device for determining carbon reduction demand response strategy, electronic equipment and storage medium | |
| CN117559532A (en) | Modeling method and device for virtual synchronous power distribution network and control platform of modeling device | |
| CN117252444A (en) | Power distribution network planning method and device, electronic equipment and storage medium | |
| CN117422260A (en) | State estimation method, device, equipment and medium based on graph neural network | |
| CN117035734A (en) | Importance degree calculation method, device, equipment and medium for substation node | |
| CN115549213A (en) | Distributed power supply access scale determination method, device, electronic device, storage medium, and program product | |
| CN117422259A (en) | Distributed resource prediction aggregation method, device, equipment and storage medium | |
| CN117477543A (en) | Power distribution system stability evaluation system, method and electronic equipment |
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 |