CN104933481B - Reference electric network model and method for solving for electric system assessment and progressive planning - Google Patents
Reference electric network model and method for solving for electric system assessment and progressive planning Download PDFInfo
- Publication number
- CN104933481B CN104933481B CN201510333979.9A CN201510333979A CN104933481B CN 104933481 B CN104933481 B CN 104933481B CN 201510333979 A CN201510333979 A CN 201510333979A CN 104933481 B CN104933481 B CN 104933481B
- Authority
- CN
- China
- Prior art keywords
- under
- load
- model
- node
- vector
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 22
- 230000000750 progressive effect Effects 0.000 title claims abstract description 17
- 238000010248 power generation Methods 0.000 claims abstract description 36
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 13
- 238000010276 construction Methods 0.000 claims abstract description 10
- 239000013598 vector Substances 0.000 claims description 60
- 239000011159 matrix material Substances 0.000 claims description 32
- 230000005611 electricity Effects 0.000 claims description 19
- 230000035945 sensitivity Effects 0.000 claims description 13
- 239000000243 solution Substances 0.000 claims description 11
- 230000009467 reduction Effects 0.000 claims description 8
- 238000002347 injection Methods 0.000 claims description 5
- 239000007924 injection Substances 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 241000039077 Copula Species 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 4
- 238000005520 cutting process Methods 0.000 claims description 4
- 230000009977 dual effect Effects 0.000 claims description 4
- 230000008602 contraction Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims 1
- 240000002853 Nelumbo nucifera Species 0.000 claims 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims 1
- 230000005540 biological transmission Effects 0.000 description 11
- 238000005457 optimization Methods 0.000 description 8
- 238000012545 processing Methods 0.000 description 7
- 238000011161 development Methods 0.000 description 5
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- WCUXLLCKKVVCTQ-UHFFFAOYSA-M Potassium chloride Chemical compound [Cl-].[K+] WCUXLLCKKVVCTQ-UHFFFAOYSA-M 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000004870 electrical engineering Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000002803 fossil fuel Substances 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
Classifications
-
- 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
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses the reference electric network models and method for solving for electric system assessment and progressive planning, including:It establishes and refers to power grid archetype;It will be changed with reference to power grid archetype:The operating cost part of system carries out model solution to two class variables respectively later to which the variable of entire model is divided into two classes after addition line disconnection in object function and constraints;Build primal problem model, the power generation dispatching under the model treatment normal operating condition and the optimal construction capacity problem of circuit;Build subproblem model, the system operatio cost minimization problem under the model treatment line disconnection accident condition;Contacting between structure primal problem and subproblem.To moral, this decomposition method is applied in the solution with reference to electric network model the application, model overall size is set to be greatly lowered, the reference electric network model that large-scale electrical power system can effectively be solved, improves solution efficiency so that model can be applied to the progressive programme of system.
Description
Technical field
The present invention relates to electrical engineering fields, more particularly, to the reference power grid mould of electric system assessment and progressive planning
Type and method for solving.
Background technology
Electric power enterprise is as foundation of national economy industry, and the development of itself must be ahead of the development of macroeconomic, in this way
Guarantee could be provided for economic development.Power system network planning is the premise for ensureing electric power and reliably supplying.With renewable energy
The progress of the market-oriented reform of the fast development of source power generation, the interconnection of large-scale electric system, and China, electric system
Network planning issue is faced with more stern challenge.
It is increasingly paid close attention to the influence of combustion of fossil fuel environmental pollution, various renewable energy power generations are vigorously sent out
Exhibition.Due to renewable energy power generation have the characteristics that unique distributivity, power generation it is intermittent and uncertain, lead to system
Line Flow is distributed more irregular, Line Flow fluctuation reinforcement, predictable variation.The main purpose of the application is to carry
Go out the mechanism that a kind of performance to system entirety provides rational evaluation, planning construction is carried out on this manufacturing basis, can either be carried
High energy source utilization rate, while it being avoided that the system excessively wasting of resources caused by construction again.
There is prodigious advantage for assessment system overall performance with reference to the concept of power grid.It is initially from economy with reference to power grid
What angle proposed, the qualitative calculating for System Margin cost determines optimal transmission price of network etc..Power network development is referred to later
Go out many different meanings, therefore occurs many different types of with reference to electric network model.It is an ideal system with reference to power grid,
Actual electric power networks can be compared therewith.It is used for the progressive planning of electric system with reference to the concept of power grid, planning can be made
Preferably cope with impact of the various new energy accesses to electric system.
In order to realize the use above for referring to power grid, especially implement progressive programme, it is desirable that can locate with reference to power grid
Large scale electric network model is managed, while to have higher solving speed.Therefore, the present invention is proposed to changing with reference to electric network model
Into wherein focusing in the modeling with reference to power grid and the method to model treatment.
Invention content:
In order to solve the above-mentioned technical problem, the reference power grid that the present invention is built not only allows for the load wave in time scale
Emotionally condition can be changed with the system architecture in processing space.System need to add new load bus, power generation node or
When needing to add new circuit, very easily these variations can be taken into account with reference to power grid;The present invention proposes the moral that will run quickly
This decomposition algorithm is applied in the solution procedure with reference to electric network model, to improve ability and the processing of processing big data quantity problem
Speed is realized and quickly calculates analysis.
For the reference electric network model of electric system assessment and progressive planning, specifically include:
It establishes and refers to power grid archetype;
It will be changed with reference to power grid archetype:The system after addition line disconnection in object function and constraints
Operating cost part carries out model solution to two class variables respectively later to which the variable of entire model is divided into two classes;
Primal problem model is built, the power generation dispatching and the optimal construction capacity of circuit under the model treatment normal operating condition are asked
Topic;
Subproblem model is built, the system operatio cost minimization problem under the model treatment line disconnection accident condition,
It is required that system loss load and loss generated energy are in setting range;
Contacting between structure primal problem and subproblem, that is, formed and extreme value is determined by subproblem and according to the principle of duality
Moral this cutting plane of running quickly of structure, the form that this cuts to moral can feed back to the result information that subproblem solves in primal problem, and
Effect of contraction is played to the solution of primal problem variable.
The explanation of concrete model:
The object function of the reference power grid with reference to power grid archetype is:
Establishing primal problem object function corresponding constraints in normal event is:
AkFpL(k)-PpG(k)+PpD=0 (2)
FpL(k)=H(k)(PpG(k)-PpD) (3)
|F|≤T (4)
Wherein, ng is the quantity of generating set;Nd is the quantity of system loading node;Nl is the quantity of circuit in power grid;
Np is the quantity of load Time segments division, and p is p-th of period in the load period;Nc is total anticipation in the event sets considered
Event number, k are k-th of event in event sets, and k=0 indicates that normal event, k ≠ 0 indicate anticipation event;τpFor load
The duration of period p;CgFor the operating cost of generating set g, CbFor the year cost of investment of circuit b;PgFor generating set g's
Generated energy, PpgFor the generated output of generating set g in load period p, PGGo out force vector, P for generating setpGFor load period p
Generating set go out force vector;T is the capacity vector of system line, TlFor the capacity of circuit l;F is Line Flow vector, Fplk
For trend vectors of the circuit l under the p periods under accident and non-accident condition;A is the node-of whole system under normal operating condition
Branch incidence matrix, AkFor the node-branch incident matrix of whole system under k-th of anticipation event;H is the branch of whole system
Road-node sensitivity matrix, HkFor branch-node sensitivity matrix of whole system under k-th of anticipation event;PGmaxFor power generation
Unit generation upper limit vector, PGminFor generating set power generation lower limit vector;PDFor node load vector, PpDFor the node of period p
Load vector.
Formula (1) includes that system minimizes year cost of electricity-generating and year Transmission Cost;Formula (2) is each event lower node power-balance
Equation;Formula (3) is the equation that Branch Power Flow is injected about node power under each event, and (4) indicate power transmission line under each event
Road effective power flow, which is less than, builds capacity limit;Formula (5) indicates that generating set active power output bound constrains under each event.
Note that the formula (2) write out here is the node power equation obtained according to Kirchhoff's current law (KCL), also
It is that requirement injects the power of a certain node equal to all power for flowing out the node;Formula (3) is according to Kirchhoff's second law
Relationship between obtained system branch trend and node injection rate, since sensitivity matrix has been arrived in use in the formula, so it
The information such as framework condition and the line reactance of whole system are reflected.
The capacity vector value range of circuit is:
0≤T≤∞ (6)
Sensitivity matrix H reflects the incidence relation between node injection and Branch Power Flow, and definition is as follows:
Wherein, YdIt is the diagonal matrix of branch admittance matrix, A is that node-branch of whole system under normal operating condition closes
Join matrix, YrbusIt is the matrix that the row and column of system admittance battle array removal buffer joint obtains, which is nonsingular, is being handled
In can carry out inversion operation.
Object function with reference to power grid archetype is formal and is unsatisfactory for, to the requirement of this decomposition algorithm of moral, needing to mesh
Scalar functions are changed, and are specifically changed into:
Wherein, nc is the event sets considered, dpkg:G generates electricity node generator output less than lower limit under k accident conditions
The generated energy for simulating reduction afterwards is equivalent to power generation node input load;
rpkj:Jth load bus lacks load under k accident conditions, is equivalent to load bus input generator;
ρ:The operating cost coefficient of system under different accidents.
The object function of archetype is changed to formula (8) form, it is notable that here, former problem object function is
Meet the form to this universal model of moral, the influence of the stoppage in transit accident that may occur with power generating equipment that should also count and transmit electricity,
And system needs to make corresponding emergency processing after such phenomenon occurs, emergency processing is often the adjustment of the method for operation, very
To the excision and stoppage in transit for being load or generating set, all these operating process should all have the cost of its operation, because
This is specifically contemplated that these accidental operation expenses in the structure of optimization object function.So here, the target letter of optimization problem
Number should include finally three parts:Minimize year cost of electricity-generating, year Transmission Cost, emergency episode operating cost.
The structure of primal problem model, concrete model are as follows:
Primal problem object function is:
Normal event is corresponding to be constrained to:
∑PG=∑ PD (10)
Wherein, ZpcFor the equivalent variable of operating cost under system p periods, k accidents;PpDjFor the load of p period nodes j
Amount, remaining variables, which are shown in archetype, to be illustrated.
Wherein:Formula (9) is object function, including two large divisions, when operating cost and capacity of trunk Installed capital cost,
Second is that accident condition operating cost forms the part such as z in formula (9) by transformation, target is to take this two parts expense summation minimum
Value, due to the optimization problem be finally reached the result is that any forecast accident occurs after, the operating cost of system is zero, that is,
Say that second part optimum results are zero, therefore coefficient ρ here should take bigger numerical.
Formula (10) is the equilibrium equation of the power generation and load that do not consider network loss;Formula (11) indicates the transmission of electricity active transmission energy of element
Force constraint;Formula (12) is generating set active power output restriction.
Subproblem should be the model that guarantee system occurs that various accident situations remain to normal operation, accident condition here
It is reduced to circuit and single failure, i.e. disconnection fault occurs.After line disconnection, other working line trends necessarily change, and have
The trend of a little circuits may increase, or even increase above the capacity of trunk upper limit that primal problem obtains and at this moment just will appear two kinds
Situation cannot obtain enough electricity supplies one is the load of certain nodes, and excision load phenomenon occur;Another situation
It is that generator runs connected circuit with minimum power generation lower limit and is still not able to the electricity that generator is sent out being transported to
Load at this time just will appear excision generator phenomenon.During Power System Planning, both of which should be kept away
Exempt from.Therefore, in the building of subproblem model, both of which is taken into account, that is, be equivalent to and add in load bus
Enter generator, load is added in generator node.
The structure of subproblem model:Subproblem concrete model is as follows:
Constraints has:
F=Hkg*(PGk+dk)-Hkd*(PD-rk) (14)
0≤rk≤PD (17)
|Fkl|≤T0l (20)
Wherein:
rkj:Jth load bus lacks load under k accident conditions, is equivalent to load bus input generator;
dkg:G generators node generator output is less than the generated energy for simulating reduction after lower limit under k accident conditions, quite
Load is put into power generation node;
rk:Load bus lacks load vector under k accident conditions;
dk:Generator node generator output is less than the generated energy vector for simulating reduction after lower limit under k accident conditions;
PD:Node load amount vector;
PGk:Generator node generated energy vector under k accident conditions;
Pgk:G generator nodes power generation amount lost under k accident conditions;
PG0:The generator output vector that primal problem obtains;
Fk:Line Flow vector under k accident conditions;
Fkl:Trend under k accident conditions on circuit l
T0l:The optimal capacity of circuit l that primal problem determines;
Hk:Branch Power Flow injects the sensitivity matrix of electricity about node under k accident conditions, here, is opened for certain circuit
Matrix after disconnected.It is the incidence relation of generator node and load bus and branch that subscript g, d, which indicates corresponding,.
The constraints of subproblem is specifically contemplated that Line Flow cannot surmount in the capacity of trunk acquired under normal condition
Limit, such as formula (20);The output bound of generating set notices that output here should not solve to obtain more than normal operating condition
Generating set contribute, such as formula (15), the output upper limit of generator should be also must not exceed, such as formula (16);Simulate load bus
Losing load should not overload the load of node, such as formula (17);Simulating generator node loss generated energy should not
More than the node generated energy upper limit, such as formula (18);And the power-balance of system entirety pays attention to energy balance here such as formula (19)
Contain generated energy, load, load amount lost, four part of power generation amount lost.
After subproblem model foundation, need to build contacting between primal problem and subproblem, that is, formed and asked by son
Topic determines that extreme value, moral this cutting plane of running quickly according to principle of duality structure, specific cut set expression formula are:
πgk1×PG+(πfk1-πfk2)×T+E≤Zpk (21)
pg0:Primal problem it needs to be determined that generator output;πf1、πf2:Branch Power Flow constrains corresponding Lagrange multiplier,
Corresponding 20);
E:Rest part in subproblem dual problem object function other than copula related to primal problem aleatory variable.
For the method for solving of electric system assessment and the reference electric network model of progressive planning, specially:
Step 1:It calculates and obtains initial power generation vector sum circuit appearance with reference to the primal problem under power grid normal operating condition
Measure the upper limit;It is UB=+ ∞ to give the initial upper bound, and initial lower bound is LB=-∞, the convergence mark using the difference of bound as algorithm
It is accurate;
Step 2:The generating set generated energy and line construction capacity that primal problem determines are brought into as known quantity to moral
This is decomposed in subproblem, solves the operating cost for obtaining system to this subproblem of moral and specific generating set loses generated energy
Load is lost with load bus, the system operation and circuit that the system operatio expense after subproblem is determined is determined with primal problem
Construction cost total is assigned to LB;
Step 3:According to subproblem acquire as a result, formed run quickly moral this cut such as following formula, bring into primal problem as constraint item
Part:
πg1×pg0+(πf1-πf2)×T+E≤Zpc
To moral, this cuts the primal problem of constraint for the addition looked for novelty, and obtains new power generation capacity, capacity of trunk and newly optimal
Target function value assigns new primal problem target function value to LB, if UB=LB, algorithm terminates, and otherwise, goes to step (2-2);
Wherein, pg0:Primal problem it needs to be determined that generator output;T is the capacity vector of circuit;The two is herein as waiting for
Determine variable;
πgk1:The Lagrange multiplier vector that the upper limit that generates electricity under k accident conditions constrains;
πfk1、πfk2:Branch Power Flow constrains corresponding Lagrange multiplier vector under k accident conditions;
Zpk:Second part in primal problem object function is corresponded under k accident conditions;
E:Under k accident conditions in subproblem dual problem object function other than copula related to primal problem aleatory variable
Remaining real part, be the form of vector here.
Beneficial effects of the present invention:
1. the reference power grid that the present invention is built not only allows for the load fluctuation situation in time scale, sky can also be handled
Between on system architecture change.When system needs to add new load bus, power generation node or need to add new circuit,
Very easily these variations can be taken into account with reference to power grid, realize and quickly calculate analysis;
2. the present invention proposes this decomposition method to moral being applied in the solution with reference to electric network model, process problem is improved
Ability and calculating speed, propose that progressive programme provides practical tool for after.
Description of the drawings
Fig. 1 is the flow chart of the present invention.
Specific implementation mode:
The present invention is described in detail below in conjunction with the accompanying drawings:
As shown in Figure 1, the application proposes the concept with reference to power grid, and construct model, to the model propose with run quickly moral this
The method of decomposition is solved, and solution procedure builds the model of primal problem and subproblem respectively, between primal problem and subproblem
By run quickly moral this cut contact, improve the efficiency of model solution.To which the implementation for follow-up progressive programme provides effectively
Handling implement.Including following components:
(1) power grid archetype is referred to;
(2) this decomposition primal problem model of moral of running quickly;
(3) this decomposition subproblem model of moral of running quickly;
(4) to moral, this cuts;
(5) power grid is referred to this decomposition method step of moral.
Refer to power grid archetype in (1), model foundation system operation cost is united with system line construction cost
The one coordination optimization mechanism considered together ensure that system has maximum social benefit;
To the primal problem model of this decomposition model of moral in (2), the object function of model includes two large divisions:A part
It is operating cost and capacity of trunk Installed capital cost, another part is accident condition operating cost.Constraint includes not considering net
The power generation of damage and the active transmittability constraint of the equilibrium equation of load, transmission line of electricity, generating set active power output restriction.
To the subproblem model of this decomposition model of moral in (3), the object function of model includes two parts:A part is
Load bus loses load, and a part is power generation node loss generated energy, and two parts take the minimum of sum.Constraints considers
The output bound that the capacity of trunk upper limit, generating set that are acquired under normal condition cannot be surmounted to Line Flow, pays attention to here
Contribute should not be more than the generating set that solves of normal operating condition contribute, must not exceed generator the output upper limit,
Simulation load bus lose load should not overload node load, simulation generator node loss generated energy do not answer
This is more than the power-balance of the node generated energy upper limit and system entirety, notices that energy balance contains generated energy, load here
Amount, load amount lost, four part of power generation amount lost.
The structure that this cuts to moral in (4), the result information that the form that this cuts to moral should be able to solve subproblem
It feeds back in primal problem, and effect of contraction is played to the solution of primal problem variable.
More detailed scheme is as follows:
The concrete model and method for solving of the present invention is described as follows:
1. with reference to the mathematical description of electric network model
Indicate as follows with reference to the object function of power grid:
Normal event and anticipation event is corresponding is constrained to
AkFpL(k)-PpG(k)+PpD=0 (2)
FpL(k)=H(k)(PpG(k)-PpD) (3)
|F|≤T (4)
Wherein, ng is the quantity of generating set;Nd is the quantity of system loading node;Nl is the quantity of circuit in power grid;
Np is the quantity of load Time segments division, and p is p-th of period in the load period;Nc is total anticipation in the event sets considered
Event number, k are k-th of event in event sets, and k=0 indicates that normal event, k ≠ 0 indicate anticipation event;τpFor load
The duration of period p;CgFor the operating cost of generating set g, CbFor the year cost of investment of circuit b;PgFor generating set g's
Generated energy, PpgFor the generated output of generating set g in load period p, PGGo out force vector, P for generating setpGFor load period p
Generating set go out force vector;T is the capacity vector of system line, TlFor the capacity of circuit l;F is Line Flow vector, Fplk
For trend vectors of the circuit l under the p periods under accident and non-accident condition;A is the node-of whole system under normal operating condition
Branch incidence matrix, AkFor the node-branch incident matrix of whole system under k-th of anticipation event;H is the branch of whole system
Road-node sensitivity matrix, HkFor branch-node sensitivity matrix of whole system under k-th of anticipation event;PGmaxFor power generation
Unit generation upper limit vector, PGminFor generating set power generation lower limit vector;PDFor node load vector, PpDFor the node of period p
Load vector.
Formula (1) includes that system minimizes year cost of electricity-generating and year Transmission Cost;Formula (2) is each event lower node power-balance
Equation;Formula (3) is the equation that Branch Power Flow is injected about node power under each event, and (4) indicate power transmission line under each event
Road effective power flow, which is less than, builds capacity limit;Formula (5) indicates the constraint of generating set active power output bound.
Since optimization problem is to solve the optimal capacity of circuit, the value range of the variable is:
0≤T≤∞ (6)
Note that the formula (2) write out here is the node power equation obtained according to Kirchhoff's current law (KCL), also
It is that requirement injects the power of a certain node equal to all power for flowing out the node;Formula (3) is according to Kirchhoff's second law
Relationship between obtained system branch trend and node injection rate, since sensitivity matrix has been arrived in use in the formula, so it
The information such as framework condition and the line reactance of whole system are reflected.
Sensitivity matrix H reflects the incidence relation between node injection and Branch Power Flow, and definition is as follows:
Wherein, YdIt is the diagonal matrix of branch admittance matrix.YrbusBe system admittance battle array removal buffer joint row and column obtain
Matrix, which is nonsingular, can carry out inversion operation in processes.
2. to this decomposition method of moral
Object function with reference to power grid archetype is formal and is unsatisfactory for, to the requirement of this decomposition algorithm of moral, needing to mesh
Scalar functions carry out some changes:
Wherein,
dkg:G generates electricity node generator output less than the generated energy for simulating reduction after lower limit under k accident conditions, is equivalent to
The node that generates electricity puts into load;
rkj:Jth load bus lacks load under k accident conditions, is equivalent to load bus input generator;
ρ:The operating cost coefficient of system under different accidents.
The object function of archetype is changed to formula (8) form, it is notable that here, former problem object function is
Meet the form to this universal model of moral, the influence of the stoppage in transit accident that may occur with power generating equipment that should also count and transmit electricity,
And system needs to make corresponding emergency processing after such phenomenon occurs, emergency processing is often the adjustment of the method for operation, very
To the excision and stoppage in transit for being load or generating set, all these operating process should all have the cost of its operation, because
This is specifically contemplated that these accidental operation expenses in the structure of optimization object function.So here, the target letter of optimization problem
Number should include finally three parts:Minimize year cost of electricity-generating, year Transmission Cost, emergency episode operating cost.
The structure of 2.1 primal problem models
Concrete model is as follows:
Primal problem object function is:
Normal event is corresponding to be constrained to:
∑PG=∑ PD(10)
Wherein:
Formula (9) is object function, including two large divisions, when operating cost and capacity of trunk Installed capital cost, second is that
Accident condition operating cost forms the part such as z in formula (9) by transformation, and target is to take this two parts expense summation minimum value,
Due to the optimization problem be finally reached the result is that any forecast accident occurs after, the operating cost of system is zero, that is to say, that
Second part optimum results are zero, therefore coefficient ρ here should take bigger numerical.Formula (10) be do not consider the power generation of network loss with
The equilibrium equation of load;Formula (11) indicates the active transmittability constraint of transmission of electricity element;Formula (12) limits for generating set active power output
Restrict beam.
The structure of 2.2 subproblem models
Subproblem should be the model that guarantee system occurs that various accident situations remain to normal operation, accident condition here
It is reduced to circuit and single failure, i.e. disconnection fault occurs.After line disconnection, other working line trends necessarily change, and have
The trend of a little circuits may increase, or even increase above the capacity of trunk upper limit that primal problem obtains and at this moment just will appear two kinds
Situation cannot obtain enough electricity supplies one is the load of certain nodes, and excision load phenomenon occur;Another situation
It is that generator runs connected circuit with minimum power generation lower limit and is still not able to the electricity that generator is sent out being transported to
Load at this time just will appear excision generator phenomenon.During Power System Planning, both of which should be kept away
Exempt from.Therefore, in the building of subproblem model, both of which is taken into account, that is, be equivalent to and add in load bus
Enter generator, load is added in generator node.
Subproblem concrete model is as follows:
Constraints has:
F=Hkg*(PGk+dk)-Hkd*(PD-rk) (14)
0≤rk≤PD (17)
|Fkl|≤T0l (20)
Wherein:
rkj:Jth load bus lacks load under k accident conditions, is equivalent to load bus input generator;
dkg:G generators node generator output is less than the generated energy for simulating reduction after lower limit under k accident conditions, quite
Load is put into power generation node;
rk:Load bus lacks load vector under k accident conditions;
dk:Generator node generator output is less than the generated energy vector for simulating reduction after lower limit under k accident conditions;
PD:Node load amount vector;
PGk:Generator node generated energy vector under k accident conditions;
Pgk:G generator nodes power generation amount lost under k accident conditions;
PG0:The generator output vector that primal problem obtains;
Fk:Line Flow vector under k accident conditions;
Fkl:Trend under k accident conditions on circuit l
T0l:The optimal capacity of circuit l that primal problem determines;
Hk:Branch Power Flow injects the sensitivity matrix of electricity about node under k accident conditions, here, is opened for certain circuit
Matrix after disconnected.It is the incidence relation of generator node and load bus and branch that subscript g, d, which indicates corresponding,.
After subproblem model foundation, need to build contacting between primal problem and subproblem, that is, formed and asked by son
Topic determines that extreme value, moral this cutting plane of running quickly according to principle of duality structure, specific cut set expression formula are:
πgk1×PG+(πfk1-πfk2)×T+E≤Zpk (21)
Wherein,
πgk1:The Lagrange multiplier vector that the upper limit that generates electricity under k accident conditions constrains, corresponding (3-10);
πfk1、πfk2:Branch Power Flow constrains corresponding Lagrange multiplier vector under k accident conditions;
Zpk:Second part in primal problem object function is corresponded under k accident conditions, sees formula (3-4b);
E:Under k accident conditions in subproblem dual problem object function other than copula related to primal problem aleatory variable
Remaining real part, be the form of vector here, do not write out expression.
2.3 algorithm basic steps
Step 1:The primal problem under normal operating condition is calculated, pays attention to primal problem here only and include the in formula (3-4b)
A part.Obtain the initial power generation vector sum capacity of trunk upper limit.It is UB=+ ∞ to give the initial upper bound, and initial lower bound is LB=-
∞, upper bound lower bound are used for determining whether to restrain, and cycle is completed to need to update every time.
Step 2:By current fixed variate-value, solve to this subproblem of moral.It enables:
Obtained variable value will be solved by, which paying attention to, brings into formula, if UB=LB, algorithm terminates.Otherwise 3 are gone to step.
Step 3:Constraint formula is added in the relaxation problem of primal problem:
πg1×pg0+(πf1-πf2)×T+E≤Zpc
The relaxation problem looked for novelty, and obtain new power generation, capacity of trunk and new optimal objective function value.By new examination in chief
It inscribes target function value and assigns LB.If UB=LB, algorithm terminates, and otherwise enables t=t+1, goes to step 2.
Above-mentioned, although the foregoing specific embodiments of the present invention is described with reference to the accompanying drawings, not protects model to the present invention
The limitation enclosed, those skilled in the art should understand that, based on the technical solutions of the present invention, those skilled in the art are not
Need to make the creative labor the various modifications or changes that can be made still within protection scope of the present invention.
Claims (6)
1. the reference electric network model for electric system assessment and progressive planning, characterized in that specifically include:
It establishes and refers to power grid archetype;
It will be changed with reference to power grid archetype:In object function and constraints be added line disconnection after system operation
Cost components carry out model solution to two class variables respectively later to which the variable of entire model is divided into two classes;
Build primal problem model, the power generation dispatching under the model treatment normal operating condition and the optimal construction capacity problem of circuit;
Build subproblem model, the system operatio cost minimization problem under the model treatment line disconnection accident condition, it is desirable that
System loss load and loss generated energy are in setting range;
Contacting between structure primal problem and subproblem, that is, formed and determined extreme value by subproblem and built according to the principle of duality
Moral this cutting plane of running quickly, the form that this cuts to moral can feed back to the result information that subproblem solves in primal problem, and to master
Effect of contraction is played in the solution of problem variable;
The object function of the reference power grid with reference to power grid archetype is:
Establishing primal problem object function corresponding constraints in normal event is:
AkFpL(k)-PpG(k)+PpD=0 (2)
FpL(k)=H(k)(PpG(k)-PpD) (3)
|F|≤T (4)
Wherein, ng is the quantity of generating set;Nd is the quantity of system loading node;Nl is the quantity of circuit in power grid;Np is
The quantity of load Time segments division, p are p-th of period in the load period;Nc is total anticipation event in the event sets considered
Quantity, k are k-th of event in event sets, and k=0 indicates that normal event, k ≠ 0 indicate anticipation event;τpFor load period p
Duration;CgFor the operating cost of generating set g, CbFor the year cost of investment of circuit b;PgFor the power generation of generating set g
Amount, PpgFor the generated output of generating set g in load period p, PGGo out force vector, P for generating setpGFor the hair of load period p
Electric unit output vector;T is the capacity vector of system line, TlFor the capacity of circuit l;F is Line Flow vector, FplkFor line
Trend vectors of the road l under the p periods under accident and non-accident condition;A is node-branch of whole system under normal operating condition
Incidence matrix, AkFor the node-branch incident matrix of whole system under k-th of anticipation event;H is branch-section of whole system
Point sensitivity matrix, HkFor branch-node sensitivity matrix of whole system under k-th of anticipation event;PGmaxIt is sent out for generating set
Electric upper limit vector, PGminFor generating set power generation lower limit vector;PDFor node load vector, PpDIt is sweared for the node load of period p
Amount;
Object function with reference to power grid archetype is formal and is unsatisfactory for, to the requirement of this decomposition algorithm of moral, needing to target letter
Number is changed, and is specifically changed into:
Wherein, nc is the event sets considered, dpkg:G generates electricity node generator output less than lower limit rear mold under k accident conditions
The quasi- generated energy cut down is equivalent to power generation node input load;
rpkj:Jth load bus lacks load under k accident conditions, is equivalent to load bus input generator;
ρ:The operating cost coefficient of system under different accidents.
2. the reference electric network model for electric system assessment and progressive planning as described in claim 1, characterized in that circuit
Capacity vector value range be:
0≤T≤∞ (6)。
3. the reference electric network model for electric system assessment and progressive planning as described in claim 1, characterized in that sensitive
Degree matrix H reflects the incidence relation between node injection and Branch Power Flow, and definition is as follows:
Wherein, YdIt is the diagonal matrix of branch admittance matrix, A is that node-branch of whole system under normal operating condition is associated with square
Battle array, YrbusIt is the matrix that the row and column of system admittance battle array removal buffer joint obtains, which is nonsingular, in processes may be used
Carry out inversion operation.
4. the reference electric network model for electric system assessment and progressive planning as described in claim 1, characterized in that examination in chief
The structure of model is inscribed, concrete model is as follows:
Primal problem object function is:
Normal event is corresponding to be constrained to:
∑PG=∑ PD (10)
Wherein, ZpcFor the equivalent variable of operating cost under system p periods, k accidents;PpDjFor the load of p period nodes j,
Remaining variable, which is shown in archetype, to be illustrated.
5. the reference electric network model for electric system assessment and progressive planning as described in claim 1, characterized in that son is asked
Inscribe the structure of model:Subproblem concrete model is as follows:
Constraints has:
F=Hkg*(PGk+dk)-Hkd*(PD-rk) (14)
0≤rk≤PD (17)
|Fkl|≤T0l (20)
Wherein:
rkj:Jth load bus lacks load under k accident conditions, is equivalent to load bus input generator;
dkg:G generators node generator output is equivalent to hair less than the generated energy for simulating reduction after lower limit under k accident conditions
Electrical nodes put into load;
rk:Load bus lacks load vector under k accident conditions;
dk:Generator node generator output is less than the generated energy vector for simulating reduction after lower limit under k accident conditions;
PD:Node load amount vector;
PGk:Generator node generated energy vector under k accident conditions;
Pgk:G generator nodes power generation amount lost under k accident conditions;
PG0:The generator output vector that primal problem obtains;
Fk:Line Flow vector under k accident conditions;
Fkl:Trend under k accident conditions on circuit l
T0l:The optimal capacity of circuit l that primal problem determines;
Hk:Branch Power Flow injects the sensitivity matrix of electricity about node under k accident conditions, is the square after certain line disconnection
Battle array, it is the incidence relation of generator node and load bus and branch that subscript g, d, which indicates corresponding,.
6. based on the method for solving described in claim 1 for electric system assessment and the reference electric network model of progressive planning,
It is characterized in that specially:
Step 1:It calculates with reference to the primal problem under power grid normal operating condition, obtains in initial power generation vector sum capacity of trunk
Limit;It is UB=+ ∞ to give the initial upper bound, and initial lower bound is LB=-∞, using the difference of bound as the convergence of algorithm;
Step 2:The generating set generated energy and line construction capacity that primal problem determines are brought into as known quantity to this point of moral
In solution subproblem, solves the operating cost for obtaining system to this subproblem of moral and specific generating set is lost generated energy and born
Lotus node loss load, the system operation and line construction that the system operatio expense after subproblem is determined is determined with primal problem
Expense total is assigned to LB;
Step 3:According to subproblem acquire as a result, formed run quickly moral this cut such as following formula, bring into and be used as constraints in primal problem:
πg1×pg0+(πf1-πf2)×T+E≤Zpc
To moral, this cuts the primal problem of constraint for the addition looked for novelty, and obtains new power generation capacity, capacity of trunk and new optimal objective
Functional value assigns new primal problem target function value to LB, if UB=LB, algorithm terminates, and otherwise, goes to step (2-2);
Wherein, pg0:Primal problem it needs to be determined that generator output;T is the capacity vector of circuit;The two is herein as change undetermined
Amount;
πgk1:The Lagrange multiplier vector that the upper limit that generates electricity under k accident conditions constrains;
πfk1、πfk2:Branch Power Flow constrains corresponding Lagrange multiplier vector under k accident conditions;
Zpk:Second part in primal problem object function is corresponded under k accident conditions;
E:Its under k accident conditions in subproblem dual problem object function other than copula related to primal problem aleatory variable
Remaining real part is the form of vector here.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510333979.9A CN104933481B (en) | 2015-06-16 | 2015-06-16 | Reference electric network model and method for solving for electric system assessment and progressive planning |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510333979.9A CN104933481B (en) | 2015-06-16 | 2015-06-16 | Reference electric network model and method for solving for electric system assessment and progressive planning |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104933481A CN104933481A (en) | 2015-09-23 |
CN104933481B true CN104933481B (en) | 2018-07-17 |
Family
ID=54120641
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510333979.9A Active CN104933481B (en) | 2015-06-16 | 2015-06-16 | Reference electric network model and method for solving for electric system assessment and progressive planning |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104933481B (en) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105846456B (en) * | 2016-05-13 | 2018-08-10 | 清华大学 | A kind of AC-DC interconnecting power network wind, fire coordinate dynamic economic dispatch optimization method |
CN106920012B (en) * | 2017-03-10 | 2021-02-19 | 国网山东省电力公司经济技术研究院 | Reference power grid evaluation method containing power flow control equipment |
CN106971239B (en) * | 2017-03-10 | 2020-09-11 | 国网山东省电力公司经济技术研究院 | Improved reference power grid evaluation method |
CN107862160B (en) * | 2017-12-06 | 2020-08-11 | 清华大学 | Method and device for generating future power grid evolution model of compressed air energy storage system |
CN108879675A (en) * | 2018-08-09 | 2018-11-23 | 国网山东省电力公司经济技术研究院 | A kind of electric network composition adaptability teaching algorithm |
CN109784605A (en) * | 2018-11-20 | 2019-05-21 | 国网浙江省电力有限公司 | Electric system spare capacity joint optimal operation method based on Transient Stability Constraints |
CN109710882B (en) * | 2018-12-21 | 2023-06-09 | 重庆大学 | Off-grid micro-grid electric vehicle ordered charge-discharge load modeling and solving method based on optimized operation |
CN113872184B (en) * | 2021-08-26 | 2023-09-01 | 国网山东省电力公司济南供电公司 | Power grid operation optimization method and system for solving overcurrent problem of power transmission line after fault |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101222141A (en) * | 2007-12-29 | 2008-07-16 | 湖南大学 | Method for selecting capacity and installation position of active electric power wave filter |
CN102195362A (en) * | 2011-05-26 | 2011-09-21 | 中国电力科学研究院 | Dynamic probability dispatching method for calculating system operation reliability |
CN102545258A (en) * | 2012-01-12 | 2012-07-04 | 西安交通大学 | Power grid optimal planning method of large-scale grid-connected wind farm |
CN103544372A (en) * | 2013-08-05 | 2014-01-29 | 中国科学院电工研究所 | Unit type selection method of micro gas turbine combined heat and power system |
-
2015
- 2015-06-16 CN CN201510333979.9A patent/CN104933481B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101222141A (en) * | 2007-12-29 | 2008-07-16 | 湖南大学 | Method for selecting capacity and installation position of active electric power wave filter |
CN102195362A (en) * | 2011-05-26 | 2011-09-21 | 中国电力科学研究院 | Dynamic probability dispatching method for calculating system operation reliability |
CN102545258A (en) * | 2012-01-12 | 2012-07-04 | 西安交通大学 | Power grid optimal planning method of large-scale grid-connected wind farm |
CN103544372A (en) * | 2013-08-05 | 2014-01-29 | 中国科学院电工研究所 | Unit type selection method of micro gas turbine combined heat and power system |
Also Published As
Publication number | Publication date |
---|---|
CN104933481A (en) | 2015-09-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104933481B (en) | Reference electric network model and method for solving for electric system assessment and progressive planning | |
CN108961018A (en) | Micro-capacitance sensor electricity market decentralization mechanism of exchange based on intelligent contract | |
CN108683192A (en) | A kind of power spot market goes out clearing method, system, equipment and storage medium | |
Xu et al. | Trajectory sensitivity analysis on the equivalent one‐machine‐infinite‐bus of multi‐machine systems for preventive transient stability control | |
CN103199521B (en) | Power network planning construction method based on network reconstruction and optimized load-flow simulating calculation | |
CN102436534B (en) | Simulating method for electric power system | |
Ghaemi et al. | Risk aversion energy management in the networked microgrids with presence of renewable generation using decentralised optimisation approach | |
Saberi et al. | A heuristic benders-decomposition-based algorithm for transient stability constrained optimal power flow | |
CN111079972A (en) | Method, device and medium for planning reliability of active power distribution network | |
Hu et al. | Collaborative optimization of distributed scheduling based on blockchain consensus mechanism considering battery-swap stations of electric vehicles | |
CN105162143B (en) | A kind of determination method of optimal trade capacity between large regional grid | |
CN110502838A (en) | Spare parts management strategy optimization model based on emulation | |
Xu et al. | Blockchain-based trading and settlement framework for electricity markets | |
CN114741956A (en) | Carbon emission flow tracking method based on graph neural network | |
CN113139271B (en) | Aggregation resource modeling method and device, electronic device and storage medium | |
CN107423259A (en) | A kind of GPU of domino optimization accelerates trigonometric equation group back substitution method on electric power | |
Li et al. | Hybrid calculation architecture of cyber physical power system based on correlative characteristic matrix model | |
Murali et al. | Estimation of locational marginal price in a restructured electricity market with different loss cases using seed genetic algorithm | |
Ding | Power System Operation with Large Scale Stochastic Wind Power Integration: Interval Arithmetic Based Analysis and Optimization Methods | |
Enshaee et al. | Approach to evaluate active loss contributions for transmission systems | |
Matcha et al. | Regular paper LMP Calculation with Distributed Loss using GA based DCOPF | |
CN109274092A (en) | Energy management method and device based on energy internet | |
CN111897237B (en) | Design method of simulation model of photovoltaic cluster based on FPGA | |
CN114357743B (en) | Edge cloud collaborative optimization method and device for regional energy Internet | |
CN115986828A (en) | Active power distribution network source storage coordination optimization operation method and system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |