Summary of the invention
The embodiment of the present invention provides a kind of region power spot market distributing and goes out clearing method, device, equipment and medium, energy
It effectively solves subdispatch center in the prior art and needs to receive and process a large amount of data, Optimization Solution calculation amount is excessive, behaviour
Make inconvenient problem.
One embodiment of the invention provides a kind of region power spot market distributing and goes out clearing method, comprising:
Overall area DC flow model is inputted, by replicating the method for boundary node variable so that the overall area direct current is damp
Flow model is divided into several subregion DC flow models;
Each subregion by the inside equation of DC flow model and inequality constraints condition and each subregion it
Between coupling constraint conditional combination obtain energy and assisted hatching goes out clear model;
The energy and assisted hatching are gone out into clear model by iterative calculation, so that each subregional boundary node phase
Angle value is equal to the angle values of overall area, exports the power output plan of each subregional generating set.
Further, further includes:
By augmentation Lagrangian Relaxation and Benders decomposition method so that the energy and assisted hatching are clear out
Model decomposition coordinates submodel, subregion scene prediction submodel and subregion error scene submodel at overall area;
The overall area coordinates submodel and exchanges boundary phase angle information with subregion scene prediction submodel, to the prediction
Scene submodel and boundary phase angle optimize;
Boundary phase angle after optimization is transferred to error submodel, by error scene submodule to prediction scene submodule
It is modified, until meeting convergence criterion, exports each subdispatch plan.
Further, the energy and assisted hatching go out clear model and includeWherein, faRepresent the first total expense of subregional prediction scene
With fa,sThe first subregional error scene abandonment expense is represented, first constraint represents subregion intra-prediction context restrictions,
Second constraint condition represents subregion internal error context restrictions formula, and third constraint condition and the 4th constraint condition represent
Coupling constraint between each subregion, PaGo out force vector, θ for the first subregional conventional generatoraFor each inside in the first subregion
Node voltage phase angle vector,WithFor the phase angle vector of the first subregional boundary node.
Further, the subregion scene prediction submodel is
Wherein,It is the phase angle correction value of kth time iterative boundary node m,It is corresponding Lagrange
Multiplier,It is corresponding quadratic penalty function multiplier,To represent the first subregional intermediate variable, total D, e is element
It is all 1 column vector, FaFor optimal cutling coefficient vector, MaAnd NaFor optimal cutling coefficient matrix.
Further,
The subregion error scene submodel includesWherein,
Pa,l, θa,lThe unit output value and node angle values that scene submodel is found out are predicted for the l times iteration.
Further, the overall area primal coordination submodel includesWherein,It is that kth time iteration subregion uploads to overall area
The boundary node angle values of Consultation Center.
Further, the convergence criterion isWherein, indicate kth time repeatedly
Generation the calculated node m in the first subregion period t angle values,Indicate the calculated section in kth time the second subregion of iteration
Point m period t angle values,Indicate the kth time calculated node n in the first subregion of iteration in the angle values of period t,
Indicate the kth time calculated node n in the second subregion of iteration in the angle values of period t.
Another embodiment of the present invention is corresponding to be provided a kind of region power spot market distributing and goes out clear device, comprising:
Divide module, for the overall area DC flow model to be divided into several subregion DC flow models;
Establish model module, for each subregion by the inside equation of DC flow model and inequality constraints condition with
And the coupling constraint conditional combination between each subregion obtains energy and assisted hatching goes out clear model;
Output module is iterated to calculate for the energy and assisted hatching to be gone out clear model, so that each subregional
Boundary node angle values are equal to the angle values of overall area, export the power output plan of each subregional generating set.
Another embodiment of the present invention is corresponding to be provided a kind of region power spot market distributing and goes out clear terminal device, including
Processor, memory and storage in the memory and are configured as the computer program executed by the processor, institute
It states and realizes that a kind of region electric power is existing as claimed in any of claims 1 to 7 in one of claims when processor executes the computer program
Goods market distributing goes out clearing method.
Another embodiment of the present invention correspondence provides a kind of computer readable storage medium, the computer-readable storage medium
Matter includes the computer program of storage, wherein controls the computer readable storage medium in computer program operation
Place equipment executes a kind of region power spot market distributing as claimed in any of claims 1 to 7 in one of claims and goes out clearing method.
Compared with prior art, a kind of region power spot market distributing disclosed by the embodiments of the present invention go out clearing method,
Device, equipment and medium, by the way that the trend constraint model of overall area is divided into several subregional tributary tide models, and root
According to the coupling constraint item between the inside equation and inequality constraints condition and each subregion DC flow model of tide model
Part combines to obtain energy and assisted hatching goes out clear model, is iterated by going out clear model to energy and assisted hatching
It calculates so that each subregional boundary node angle values are equal to the angle values of overall area, each subregional generating set of output goes out
The technical solution of power plan, so that the calculation amount of overall area reduces, it is easy to operate;Each subregion independent operating and it is subregional go out
Power plan is still issued by overall area.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
It is that a kind of region power spot market distributing that one embodiment of the invention provides goes out the stream of clearing method referring to Fig. 1
Journey schematic diagram, comprising:
S1, input overall area DC flow model, by replicating the method for boundary node variable so that the overall area is straight
Stream tide model is divided into several subregion DC flow models;
S2, each subregion pass through the inside equation and inequality constraints condition of DC flow model and each subregion
Between coupling constraint conditional combination obtain energy and assisted hatching goes out clear model;
S3, the energy and assisted hatching are gone out to clear model by iterative calculation, so that each subregional boundary section
Point angle values are equal to the angle values of overall area, export the power output plan of each subregional generating set.
Further, further includes:
By augmentation Lagrangian Relaxation and Benders decomposition method so that the energy and assisted hatching are clear out
Model decomposition coordinates submodel, subregion scene prediction submodel and subregion error scene submodel at overall area;
The overall area coordinates submodel and exchanges boundary phase angle information with subregion scene prediction submodel, to the prediction
Scene submodel and boundary phase angle optimize;
Boundary phase angle after optimization is transferred to error submodel, by error scene submodule to prediction scene submodule
It is modified, until meeting convergence criterion, exports each subdispatch plan.
S1, input overall area DC flow model, by replicating the method for boundary node variable so that the overall area is straight
Stream tide model is divided into several subregion DC flow models.
It should be noted that in the present embodiment, overall area is subdispatch center, and subregion is each province, the first subregion
Domain is to save a, and the second subregion is to save b, and provincial power output plan is reflected by the power generating value under prediction scene and error scene;By region
Control centre's decoupling is two independent region province a and saves b, and the boundary node for defining two regions is m and n, first by side
The phase angle variable duplication of boundary's node is primary, respectivelyWithAndBelong to province a,Belong to
Save b.Boundary node phase angle variable after duplication needs to meet coupling constraint:To duplication variable and establish coupling
Two province's systems after constraint establish DC power flow equation i.e. subregion DC power flow equation, mathematically and before duplication variable
DC power flow equation is mathematically equivalent.
Dispersion optimization is realized by subdispatch center.By the boundary node phase angle variable of all areas in upper layer coordinator
Second of middle duplication, is usedIt indicates, simultaneouslyTherefore, to the same boundary node m,
Saving a, province b and upper layer coordinator has one group of variable to indicate its phase angle.
S2, each subregion pass through the inside equation and inequality constraints condition of DC flow model and each subregion
Between coupling constraint conditional combination obtain energy and assisted hatching goes out clear model.
It should be noted that minimize all conventional power units in the whole network each province (such as thermoelectricity and Hydropower Unit) in dispatching cycle
Interior total power generation expense and to dispatch the sum of expense desired value again be target, it may be assumed that
Wherein, N represents the number of all provincial power networks;NTTo dispatch total period;Represent the conventional generator number for saving a
Mesh;For active power output of the conventional power unit i in period t for saving a under prediction scene;For the routine for saving a under prediction scene
Unit i period t power generation expense, WithIt is the conventional power unit i of province a respectively
Generate electricity cost coefficient quadratic term, first order and constant term;For spinning reserve capacity of the conventional power unit i in period t for saving a;
For the spinning reserve cost coefficient of the conventional power unit i of province a;SaFor the error scene number for saving a;psFor the probability of s-th of scene,
Take 1/Sa;Cost coefficient is dispatched again to save the conventional power unit i of a;For under error scene s save a conventional power unit i when
Section t's dispatches power generating value again upwards;Power output is dispatched again downwards in period t to save the conventional power unit i of a under error scene s
Value.
It should be noted that the present embodiment for saving a, saves the equation of inside and inequality constraints includes:
Save internal node DC power flow equation:
In formula,It is to predict that conventional power unit of a in period t is saved under scene goes out force vector;It is to save a under prediction scene to exist
The node load vector of period t;It is to predict to save a under scene in the output of wind electric field vector of period t;BaIt is that province a ignores branch
Road resistance and the node admittance matrix that ground leg is set up;Be predict scene under save a period t node phase angle to
Amount;It is to save conventional power unit of a in period t under error scene s to go out force vector;It is that a is saved under error scene s in period t
Node load vector;It is that a is saved under error scene s in the wind power output vector of period t;It is to save a under error scene s to exist
The node phase angle vector of period t.
Predict that unit output bound constrains under scene:
In formula,WithFor the active power output lower and upper limit of the unit i of province a.
Wind turbines power output bound constraint:
Unit climbing is constrained with landslide:
In formula,Respectively save the active power output climbing rate of the unit i of a.
The constraint of unit spinning reserve regulating power:
Units limits are dispatched again under machine grouping error scene:
Unit up/down recalls power limit constraint:
Line transmission power constraint:
In formula,It is the maximum transmission power value of province a route kl, connection node k and node l;It is route j
Reactance value;It is phase angle of the node k in moment t.Line transmission power constraint includes saving internal wiring and getting in touch between saving
Line;It is phase angle of the scene s lower node k in moment t.
Predict that the power output regulations speed under scene and error scene constrains:
-Δi≤Pi,t-Pi,t,s≤Δi (15)
In formula, ΔiIndicate the power output increment that fired power generating unit i can be adjusted rapidly in 10 minutes.
Predict the boundary node phase angle consistency constraint under scene and error scene:
In formula, TLaaIndicate the node that province a is only belonged in boundary node, TLabIndicate its being connected in boundary node with province a
The boundary node that he saves.Formula (16) ensures the boundary node that the boundary node angle values under error scene are equal under prediction scene
Angle values.
Coupling constraint between each province
If saved between a and province b, there is interconnections between province, and their boundary node is m and n, then saving a and province
The coupling constraint of b is as shown in formula (17).
Going out clear model by the region that formula (1)-(17) are constituted is a convex programming problem, can be written as follow compact form:
In formula, faRepresent the prediction scene total cost for saving a;fa,sRepresent the error scene abandonment expense for saving a;1st constraint
It represents and saves intra-prediction context restrictions, the 2nd constraint condition represents the internal error field scape of province and constrains formula;3rd and the 4th constraint article
Part represents interdepartmental coupling constraint.PaFor the force vector out of the conventional generator of province a;θaTo save each internal node voltages phase angle of a
Vector;WithFor the phase angle vector of the boundary node of province a.
S3, the clear model out is passed through into iterative calculation, so that each subregional boundary node angle values are equal to overall area
Angle values, export the power output plan of each subregional generating set;
Further include: by augmentation Lagrangian Relaxation and Benders decomposition method so that it is described go out clear model decomposition at
Overall area primal coordination submodel, subregion scene prediction submodel and subregion error scene submodel;
It should be noted that (1), which saves power grid, predicts scene submodule
The problem of prediction scene submodule mainly reflects this province, the objective function of power saving host's problem is in this province total cost
On the basis of, coupling constraint (17) is relaxed in objective function using Augmented Lagrangian Functions, while increasing representative and missing
The intermediate variable of poor scene submodule objective function, constraint condition are generated comprising the constraint under prediction scene and by error scene
Optimal cutling, compact form are as follows:
In formula,It is the phase angle correction value of kth time iterative boundary node m, is found out using k-1 iterationWithIt calculates;It is the corresponding Lagrange multiplier of coupling constraint (17),It is corresponding quadratic penalty function multiplier;For the intermediate variable for representing region a, total D;E is the column vector that element is all 1;FaFor optimal cutling coefficient vector;MaAnd Na
For optimal cutling coefficient matrix.
Error scene is not only corrected the unit output of prediction scene, while also carrying out to the phase angle of boundary node
Amendment, after this leads to each province's random optimization, the phase angle of boundary node not necessarily be can satisfy it can also happen that change
Coupling constraint needs to carry out thus dispersion optimization again, and boundary node phase angle information is passed to phase neighbouring province by each province, willDomain power grid predicts scene submodel again after amendment.
Each zone boundary node phase angle is coordinated, it is ensured that meets coupling constraint, i.e., alternately solves the whole network dispersion optimization and ask
Topic and region stochastic optimization problems.
(2) power grid error field scape submodel is saved
The error scene submodel of the iteration of the l times random optimization, region a can be written as follow compact form:
In formula, Pa,l, θa,lThe unit output value and node angle values that scene submodel is found out are predicted for the l times iteration.
After s-th of error scene submodel solves, Benders cutting constraint item will be generated using the information solved
Part is added in prediction scene submodel.If the dual variable value of the constraint condition of the l times iteration function (17) is πa,s,l, then
The Benders that s-th of error scene generates cuts mathematical expression are as follows:
Each random optimization iteration, each error scene can generate a Benders cutting and be added to prediction scene
In model, when there are many scene number, predict that the scale of scene submodule will expand rapidly, and solve the time with iterations going on
It can increase sharply.Scene submodule is predicted in order to Efficient Solution, it is partially polymerized to Benders cutting progress, first by Sa
A error scene is divided into D class, correspondingly, by every one kind SaThe optimal cutling that/D error scene generates aggregates into 1, every in this way
Secondary iteration is added to the optimal cutling number of prediction scene submodule from SaIt is a to be reduced to D, delay prediction scene submodule
The speed that block size increases, while remaining more optimal cutling information.Optimal cutling mathematics after polymerization is expressed as:
(3) overall area coordinates submodule
Regional coordination primal problem is interior and generator is not present, and in objective function and power generation expense and abandonment fee item is not present
fa;The Augmented Lagrangian Functions item of all boundary node phase angles is contained in objective function simultaneously, regional coordination submodule can
It is write as:
In formula,It is the boundary node angle values that kth time iteration provincial power network uploads to regional coordination center.
Overall area coordinates submodule and plays coordinative role, it is therefore an objective to so that the calculated boundary node phase angle of each regional power grid
The boundary node angle values that value is calculated close to upper layer coordinator as far as possible, thus the boundary section for calculating each provincial power network
Point angle values are closer to each other, last equal.
Further, the overall area coordination submodel exchanges boundary phase angle information with subregion scene prediction submodel,
The prediction scene submodel and boundary phase angle are optimized;
Boundary phase angle after optimization is transferred to error submodel, by error scene submodule to prediction scene submodule
It is modified, until meeting convergence criterion, exports each subdispatch plan.
It should be noted that solution procedure is as follows after considering wind-powered electricity generation error scene to the amendment for respectively saving clear problem
It is shown:
(0) it initializes.Setting the number of iterations is k=1, and parameter alpha=1.05 are arranged, Convergence precision ε=10 are set-3。
(1) to not considering that the deterministic models of wind-powered electricity generation randomness carry out dispersion optimization, until meeting convergence criterion (24).By
This provides a preferable Lagrange multiplier initial value to be subsequent, is convenient for subsequent calculating, reduces and calculates the time.By length institute
Limit, repeats no more.
Wherein, indicate the kth time calculated node m in the first subregion of iteration in the angle values of period t,Indicate kth time
The calculated node m in the second subregion of iteration period t angle values,Indicate that kth time the first subregion of iteration is calculated
Node n period t angle values,Indicate the kth time calculated node n in the second subregion of iteration in the angle values of period t.
(2) random optimization is carried out to a in each regional model, a=1, small the number of iterations l=1 is set;
A) it solves and saves a prediction scene submodel (19), and the target function value for predicting scene submodel is assigned to lower limit La;
B) it solves and saves a error scene submodel (20), and obtain upper limit Ua;
C) judge small Iterations of Multi:
If formula (26) meets, then it is assumed that the random optimization for saving a has been restrained.A=a+1, into the d) step.
Otherwise, it is generated and optimal cutling and is added in primal problem, l=l+1 using the information of subproblem, return to a) step.
D) judge whether that the random optimization of all areas is restrained:
If a > N, into the e) step;
Otherwise a) step is returned;
(3) decentralized coordinating:
E) angle values of each region to adjacent area transmitting one's respective area boundary nodeWith
F) each region updates one's respective area boundary node angle values using formula (20)With
G) judge convergence:
If formula (24) meets, algorithmic statement exports each region operation plan a few days ago;
Otherwise, formula (23) are solved and update penalty function coefficient, k=k+1 returns to (2) step.
In formula, α is the parameter for adjusting step-length, generally takes 1 to 3 some value.
A kind of region power spot market distributing disclosed by the embodiments of the present invention goes out clearing method, device, equipment and medium,
By the way that the DC flow model of overall area is divided into several subregional DC flow models, and according to the constraint of tide model
Condition obtains out clear model, by the iterative calculation to clear model out so that each subregional boundary node angle values are equal to total area
The angle values in domain export the technical solution of the power output plan of each subregional generating set, so that the calculation amount of overall area reduces,
It is easy to operate;Each subregion independent operating and subregional power output plan is still issued by overall area.
It referring to fig. 2, is that a kind of region power spot market distributing that one embodiment of the invention provides goes out clear device/terminal
The structural schematic diagram of equipment, comprising:
Divide module, for the overall area DC flow model to be divided into several subregion DC flow models;
Model module is established, for the inside equation and inequality constraints item according to each subregion DC flow model
Coupling constraint conditional combination between part and each subregion DC flow model obtains energy and assisted hatching goes out clear mould
Type;
Output module is iterated to calculate for the energy and assisted hatching to be gone out clear model, so that each subregional
Boundary node angle values are equal to the angle values of overall area, export the power output plan of each subregional generating set.
Another embodiment of the present invention is corresponding to provide terminal device, including processor, memory and is stored in described deposit
In reservoir and it is configured as the computer program executed by the processor, the processor executes real when the computer program
Now a kind of region power spot market distributing goes out clearing method as claimed in any of claims 1 to 7 in one of claims.
Another embodiment of the present invention correspondence provides a kind of computer readable storage medium, the computer-readable storage medium
Matter includes the computer program of storage, wherein controls the computer readable storage medium in computer program operation
Place equipment executes a kind of region power spot market distributing as claimed in any of claims 1 to 7 in one of claims and goes out clearing method.
It referring to fig. 2, is that a kind of region power spot market distributing that one embodiment of the invention provides goes out clear device/terminal
The schematic diagram of equipment.A kind of Regional Dispersion formula of the embodiment go out clear device/terminal device include: processor, memory and
Store the computer program that can be run in the memory and on the processor.The processor executes the computer
The step in above-mentioned each embodiment of the method is realized when program, alternatively, the processor is realized when executing the computer program
The function of each module/unit in above-mentioned each Installation practice, such as establish model module.
Illustratively, the computer program can be divided into one or more module/units, one or more
A module/unit is stored in the memory, and is executed by the processor, to complete the present invention.It is one or more
A module/unit can be the series of computation machine program instruction section that can complete specific function, and the instruction segment is for describing institute
State the implementation procedure that computer program goes out in a kind of Regional Dispersion formula in clear device/terminal device.For example, the computer
Program can be divided into segmentation module, establish model module, output module, and each module concrete function is as follows: segmentation module is used
In the overall area DC flow model is divided into several subregion DC flow models;Model module is established, for according to institute
It states between the inside equation of each subregion DC flow model and inequality constraints condition and each subregion DC flow model
Coupling constraint conditional combination obtain energy and assisted hatching goes out clear model;Output module, for by the energy and auxiliary
Service market is helped to go out clear model iterative calculation, so that each subregional boundary node angle values are equal to the angle values of overall area, it is defeated
The power output plan of each subregional generating set out.
A kind of region power spot market distributing, which goes out clear device/terminal device, can be desktop PC, pen
Note sheet, palm PC and cloud server etc. calculate equipment.A kind of region power spot market distributing goes out clear device/end
End equipment may include, but be not limited only to, processor, memory.It will be understood by those skilled in the art that the schematic diagram is only
A kind of region power spot market distributing goes out the example of clear device/terminal device, does not constitute to a kind of region electric power stock
Market distributing goes out the restriction of clear device/terminal device, may include components more more or fewer than diagram, or combine certain
Component or different components, for example, a kind of region power spot market distributing go out clear device/terminal device can be with
Including input-output equipment, network access equipment, bus etc..
Alleged processor can be central processing unit (Central Processing Unit, CPU), can also be it
His general processor, digital signal processor (Digital Signal Processor, DSP), specific integrated circuit
(Application Specific Integrated Circuit, ASIC), ready-made programmable gate array (Field-
Programmable Gate Array, FPGA) either other programmable logic device, discrete gate or transistor logic,
Discrete hardware components etc..General processor can be microprocessor or the processor is also possible to any conventional processor
Deng, the processor is the control centre that a kind of region power spot market distributing goes out clear device/terminal device, utilize
Entirely a kind of region power spot market distributing goes out the various pieces of clear device/terminal device for various interfaces and connection.
The memory can be used for storing the computer program and/or module, and the processor is by operation or executes
Computer program in the memory and/or module are stored, and calls the data being stored in memory, described in realization
A kind of region power spot market distributing goes out the various functions of clear device/terminal device.The memory can mainly include depositing
Store up program area and storage data area, wherein storing program area can application program needed for storage program area, at least one function
(such as sound-playing function, image player function etc.) etc.;Storage data area, which can be stored, uses created number according to mobile phone
According to (such as audio data, phone directory etc.) etc..In addition, memory may include high-speed random access memory, can also include
Nonvolatile memory, such as hard disk, memory, plug-in type hard disk, intelligent memory card (Smart Media Card, SMC), safety
Digital (Secure Digital, SD) card, flash card (Flash Card), at least one disk memory, flush memory device or
Other volatile solid-state parts.
Wherein, a kind of region power spot market distributing goes out the integrated module/unit of clear device/terminal device such as
Fruit is realized in the form of SFU software functional unit and when sold or used as an independent product, and can store can in a computer
It reads in storage medium.Based on this understanding, the present invention realizes all or part of the process in above-described embodiment method, can also
It is completed with instructing relevant hardware by computer program, the computer program can be stored in one and computer-readable deposit
In storage media, the computer program is when being executed by processor, it can be achieved that the step of above-mentioned each embodiment of the method.Wherein, institute
Stating computer program includes computer program code, and the computer program code can be source code form, object identification code shape
Formula, executable file or certain intermediate forms etc..The computer-readable medium may include: that can carry the computer journey
Any entity or device of sequence code, recording medium, USB flash disk, mobile hard disk, magnetic disk, CD, computer storage, read-only storage
Device (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), electric carrier signal,
Telecommunication signal and software distribution medium etc..
It should be noted that the apparatus embodiments described above are merely exemplary, wherein described be used as separation unit
The unit of explanation may or may not be physically separated, and component shown as a unit can be or can also be with
It is not physical unit, it can it is in one place, or may be distributed over multiple network units.It can be according to actual
It needs that some or all of the modules therein is selected to achieve the purpose of the solution of this embodiment.In addition, device provided by the invention
In embodiment attached drawing, the connection relationship between module indicate between them have communication connection, specifically can be implemented as one or
A plurality of communication bus or signal wire.Those of ordinary skill in the art are without creative efforts, it can understand
And implement.
The above is a preferred embodiment of the present invention, it is noted that for those skilled in the art
For, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also considered as
Protection scope of the present invention.