CN104835080B - The modeling method of micro-capacitance sensor in a few days operation plan mixed-integer programming model - Google Patents

The modeling method of micro-capacitance sensor in a few days operation plan mixed-integer programming model Download PDF

Info

Publication number
CN104835080B
CN104835080B CN201510220659.2A CN201510220659A CN104835080B CN 104835080 B CN104835080 B CN 104835080B CN 201510220659 A CN201510220659 A CN 201510220659A CN 104835080 B CN104835080 B CN 104835080B
Authority
CN
China
Prior art keywords
active
power
constraint
maximum
restriction
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.)
Expired - Fee Related
Application number
CN201510220659.2A
Other languages
Chinese (zh)
Other versions
CN104835080A (en
Inventor
汪隆君
许海林
王钢
丁茂生
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
South China University of Technology SCUT
State Grid Ningxia Electric Power Co Ltd
Original Assignee
South China University of Technology SCUT
State Grid Ningxia Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by South China University of Technology SCUT, State Grid Ningxia Electric Power Co Ltd filed Critical South China University of Technology SCUT
Priority to CN201510220659.2A priority Critical patent/CN104835080B/en
Publication of CN104835080A publication Critical patent/CN104835080A/en
Application granted granted Critical
Publication of CN104835080B publication Critical patent/CN104835080B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Control Of Eletrric Generators (AREA)

Abstract

The invention discloses the modeling methods of micro-capacitance sensor in a few days operation plan mixed-integer programming model, include the following steps:1st, structure micro-capacitance sensor in a few days the object function 2 of operation plan mixed-integer programming model, build diesel oil synchronous generator active and reactive output linearisation constraint equation;3rd, the linearisation constraint equation of the active and reactive output through grid-connected inverters distributed generation resource is built;4th, the linearisation constraint equation of the active and reactive output of structure double-fed induction wind turbine;5th, structure considers the linearisation power flow equation of voltage magnitude;6th, structure tributary capacity linearisation constraint equation;7th, with reference to the existing linearisation constraint of micro-capacitance sensor, micro-capacitance sensor in a few days operation plan mixed-integer programming model is formed, solving software with mixed integer programming is solved.The modeling method of the present invention can effectively reduce the solving complexity of micro-capacitance sensor in a few days operation plan problem, have the characteristics that calculating speed is fast, solving precision is high and engineering has strong operability in fact.

Description

The modeling method of micro-capacitance sensor in a few days operation plan mixed-integer programming model
Technical field
The present invention relates to the modeling method of micro-capacitance sensor in a few days operation plan model, more particularly to a kind of micro-capacitance sensor is in a few days dispatched Plan the modeling method of mixed-integer programming model.
Background technology
As the effective means for playing distributed power generation efficiency, micro-capacitance sensor technology becomes current intelligent grid theory and technology The hot spot of research.Micro-grid operation mode is flexible, has grid-connected and two kinds of methods of operation of isolated island, and can realize two kinds as needed The seamless switching of the method for operation.However, microgrid energy management has its complexity, not only need meeting workload demand and electricity Can ensure micro-capacitance sensor economy, safe and stable operation under the premise of quality, also to coordinate different response speed power supplys output, The uncertainty of wind-powered electricity generation and photoelectricity in smooth net, particularly distributed power generation high permeability exacerbate this complexity significantly.
At present, domestic and international researcher has made related exploration to microgrid energy management strategy.A few days ago-rolling-real-time 3 layers of frame Frame and a few days ago-in a few days 2 layer frame are successively suggested.3 layers of frame and 2 layers of frame are based on predictive control theory, according to current fortune Row operating mode and newest prediction result adjust last layer controlled quentity controlled variable, the interference for the factor that removes uncertainty.Wherein, in a few days scheduling is counted Draw be in advance some time plan a few days ago in a certain period is corrected according to newest operation information, substantially consider operation pacify Full optimal load flow model.
Distributed generation resource as power generation source main in micro-capacitance sensor, mainly have diesel-driven generator (diesel engine, DE), miniature gas turbine (micro-turbine, MT), fuel cell (fuel cell, FC), wind-driven generator (wind Turbine, WT), photovoltaic cell (photo-voltaic, PV), accumulator (storage battery, ST) etc., these distributions Formula power supply can also convey the reactive power of certain capacity while active power is provided for micro-capacitance sensor.Since micro-capacitance sensor carries Certain load or burden without work, if circuit actual transmission power will likely be caused by ignoring reactive power in micro-capacitance sensor in a few days operation plan The problem of out-of-limit and electrical equipment actual capacity is out-of-limit.And the single reactive-load compensation for considering capacitor switching is then held in micro-capacitance sensor Easily cause variation.When the capacitor switching response time is slower, the idle output or absorption of distributed generation resource can also be effective The idle surplus and deficiency of smooth micro-capacitance sensor improves the voltage level of micro-capacitance sensor and reduces network loss, so cannot ignore distributed generation resource The idle effect played in micro-capacitance sensor Optimized Operation of contributing.And if idle units limits only simply establish pole to distributed generation resource Limit value limits equation, then has ignored the idle output characteristics of distributed generation resource complexity, and it is practical to be detached from engineering.
Invention content
In order to overcome the disadvantages mentioned above of the prior art, the purpose of the present invention is to provide a kind of micro-capacitance sensors in a few days to adjust with insufficient The modeling method of degree plan mixed-integer programming model.
The purpose of the present invention is achieved through the following technical solutions:
The modeling method of micro-capacitance sensor in a few days operation plan mixed-integer programming model, includes the following steps:
The object function of step 1, structure micro-capacitance sensor in a few days operation plan mixed-integer programming model is meeting normal fortune Under conditions of row constraint, by the output of each controllable of reasonable arrangement, make system operation the lowest cost;
The linearisation constraint equation of step 2, the active and reactive output of structure diesel oil synchronous generator;
Step 3, the linearisation constraint equation for building the active and reactive output through grid-connected inverters distributed generation resource;It is described Include single-shaft variant miniature gas turbine, fuel cell, accumulator, photovoltaic cell through grid-connected inverters distributed generation resource;
The linearisation constraint equation of step 4, the active and reactive output of structure double-fed induction wind turbine;
Step 5, structure consider the linearisation power flow equation of voltage magnitude;
Step 6, structure tributary capacity linearisation constraint equation;
Step 7 is constrained with reference to the existing linearisation of micro-capacitance sensor, forms micro-capacitance sensor in a few days operation plan mixed integer programming mould Type solves software with mixed integer programming and is solved.
The object function of the in a few days operation plan mixed-integer programming model of micro-capacitance sensor described in step 1, it is specific as follows:
In formula, Cdeii,t、Cmtjj,tAnd Cfck,tRespectively diesel-driven generator ii, miniature gas turbine jj and fuel cell k exist The operation totle drilling cost of period t, including fuel cost, startup with shutting down cost, maintenance cost;T ∈ T represent period t in scheduling week In phase T, ii ∈ I represent that diesel-driven generator ii belongs to diesel-driven generator set I, jj ∈ J and represents that miniature gas turbine jj belongs to micro- Type gas turbine set J, k ∈ K represent that fuel cell k belongs to fuel cell set K, s ∈ S and represents that accumulator s belongs to accumulator Set S;Css,tIt is accumulator s in period t charge and discharge voltage ageing cost;Cl,tFor load l curtailing expenditure in period t,zltFor the operating status of period t interruptible load l, pL,ltAnd PltIt is negative to represent that period t can interrupt respectively The unit power economic loss and active power of lotus l excisions;Cm,tPurchase sale of electricity expense during to be incorporated into the power networks between micro-capacitance sensor and major network With Cm,t=pgPPCC,t, pgFor period t electricity transaction price, PPCC,tIt is positive number for the active power that period t is exchanged with major network It represents power purchase, is negative number representation sale of electricity.
The linearisation constraint equation of the active and reactive output of diesel oil synchronous generator described in step 2 is synchronized by diesel oil and sent out Motor maximum stator current restriction, diesel oil synchronous generator maximum rotor voltage restriction, diesel oil synchronous generator are most Big load angle restriction, diesel oil synchronous generator minimax prime mover active power restriction composition;
Diesel oil synchronous generator maximum stator current restriction is as follows:
In formula, VsgThe voltage of node is connect by synchronous generator stator,For maximum stator current, PsgAnd QsgRespectively For active and reactive output power,For the maximum active output power of synchronous generator,For maximum rotor voltage,For Longitudinal axis reactance, PbAnd QbIt is maximum stator current constraint and the intersection point power of maximum rotor voltage constraint, PcAnd QcIt is maximum stator The intersection point power of restriction of current and prime mover maximum active power restriction, the inequality symbol of inequality constraints needs root in formula Change according to concrete numerical value and change;
Diesel oil synchronous generator maximum rotor voltage restriction is as follows:
In formula,For the minimum active output power of synchronous generator, PaAnd QaIt is the constraint of maximum rotor voltage and prime mover The intersection point power of minimum active power restriction, in formula the inequality symbol of inequality constraints need to be changed according to concrete numerical value and Variation;
Diesel oil synchronous generator peak load angle restriction is as follows:
In formula,For the maximum power angle of synchronous generator,For quadrature-axis reactance;
Diesel oil synchronous generator minimax prime mover active power restriction is as follows:
The linearisation constraint equation of the active and reactive output through grid-connected inverters distributed generation resource is by inverse described in step 3 Become device and flow through current limit constraint, the maximum active export-restriction constraint composition of distributed generation resource;
It is as follows that inverter flows through current limit constraint:
In formula, VdgThe voltage of node, I are met by inverterirTo flow through the electric current of inverter, the tune of inverter can be passed through It is whole to keep constant, Pdg、QdgThe reactive power that the active power and inverter that respectively distributed generation resource is sent out are sent out, α are inversion Device flows through the vectorial angle of the intersection point of current limit constraint and the maximum active export-restriction constraint of distributed generation resource, and θ, β are respectively nothing 2 points of vectorial angle on work(constraint curve, m linearize segments for idle circular arc, and k is integer;
The maximum active export-restriction constraint of distributed generation resource is as follows:
In formula,For the maximum output active power of distributed generation resource,Minimum output for distributed generation resource is active Power.
The linearisation constraint equation of the active and reactive output of double-fed induction wind turbine described in step 4 is by wind turbine maximum stator Current limit constraint, the constraint of wind turbine maximum rotor current limit, the active export-restriction constraint composition of wind turbine minimax;
Wind turbine maximum stator current restriction is as follows:
In formula, VwtThe voltage of node is connect by wind turbine,For the maximum stator current of wind turbine, PwtAnd QwtRespectively have Work(, idle output power,For the maximum active output power of wind turbine, swtFor the rotor speed of wind turbine, swt∈ (- 1,1), with the wind The fluctuation of speed and change, but it can be made for definite value by adjusting, Ph、QhIt is maximum for maximum stator current restriction and wind turbine The intersection point power of active export-restriction constraint, Pi、QiFor maximum stator current restriction and the minimum active export-restriction of wind turbine about The intersection point power of beam, the inequality symbol of inequality constraints need to be determined according to concrete numerical value in formula;
The constraint of wind turbine maximum rotor current limit is as follows:
In formula, the inequality symbol of inequality constraints need to be determined according to concrete numerical value, Pf、QfIt is limited for maximum rotor electric current Restrict beam and the intersection point power of the minimum active export-restriction constraint of wind turbine, Pg、QgFor the constraint of maximum rotor current limit and wind turbine most The intersection point power of big active export-restriction constraint,For excitation impedance,For stator impedance, a, b, c, d are the multiplier introduced The factor, value are as follows:
In formula,For excitation resistance,For excitation reactance,For stator resistance,For stator reactance.
The active export-restriction constraint of wind turbine minimax is as follows:
Consider that the linearisation power flow equation of voltage magnitude is as follows described in step 5:
In formula, PGiAnd QGiThe respectively generated power and reactive power of node i, PDiAnd QDiThe respectively load of node i Active and reactive power;ViAnd VjThe respectively voltage magnitude of node i and j;J ∈ i represent that node j is connected directly with node i, θij、GijAnd BijRespectively node i and phase angle difference, conductance and the susceptance of node j, GiiAnd BiiThe respectively reality of node i self-admittance Portion and imaginary part, PijAnd QijThe respectively active and reactive power of branch ij;Gi0And Bi0Respectively branch ij interior joints i is over the ground Conductance and susceptance;
Using special sort set pairMake linearization process.
The linearisation constraint equation of tributary capacity described in step 6 is as follows:
In formula, SijRepresent the transmission capacity limits value of branch ij, α ', β ' are tributary capacity restrained circleOn 2 points of vectorial angle, m' are the equilateral number of edges of inscribed, and k' is integer.
The operation principle of the present invention:By micro-capacitance sensor, in a few days operation plan nonlinear model is converted into the modeling method of the present invention Micro-capacitance sensor based on analytic geometry in a few days operation plan mixed-integer programming model.In the base of traditional micro-capacitance sensor economic load dispatching model On plinth, establish and consider the linearisation constraint equation of the active and reactive output relation of each distributed generation resource of micro-capacitance sensor, consider voltage amplitude Linearisation power flow equation, the tributary capacity linear equation based on analytic geometry of value.
Compared with prior art, the present invention has the following advantages and beneficial effect:
Firstth, the present invention is according to diesel oil synchronous generator, the distributed generation resource through grid-connected inverters and double-fed induction wind Machine establishes the difference of idle output control and linearizes idle output characteristics equation accordingly respectively, meets distributed generation resource reality Characteristic, making micro-capacitance sensor, in a few days operation plan result is more accurate.
Secondth, the present invention is based on the constraint of the linearisation tributary capacity of analytic geometry, conventional linear tide model is overcome Only consider the limitation of the active transmission capacity constraint of branch.
Third, the present invention propose the mixed-integer programming model of microgrid in a few days planned dispatching a kind of, and not linearized Model is compared, and is put forward the result of model acquisition under the premise of deviation is little, can be quickly converged on globally optimal solution, and do not depend on In initial value.In a few days operation plan is substantially mixed-integer nonlinear programming model to micro-capacitance sensor.The original of solving model continuous variable Beginning dual interior point, constringency performance depend on the distance of given initial value and optimized operation point.Since linear programming problem belongs to Convex programming problem, method for solving is ripe, and locally optimal solution is globally optimal solution;Branch and bound method is a kind of in solution sky Between gradually segment the integer programming derivation algorithm of feasible zone on tree, have global optimizing ability, with linear programming for solution algorithm phase Software, such as CPLEX are solved with reference to mixed integer programming is widely used in.Therefore, by micro-capacitance sensor in a few days operation plan model line Property after solve again, can not only efficient stable obtain globally optimal solution, also reconciled the contradiction of solving speed and precision.
Description of the drawings
Fig. 1 is the in a few days operation plan mixed-integer programming model of the micro-capacitance sensor based on analytic geometry of the embodiment of the present invention Flow chart.
Fig. 2 (a) is the active and reactive output curve diagram of not linearized diesel-driven generator of the embodiment of the present invention.
Fig. 2 (b) is the active and reactive output curve diagram of linearized diesel-driven generator of the embodiment of the present invention.
Fig. 3 (a) is the active and reactive curve of output schematic diagram with inverter distributed generation resource of the embodiment of the present invention.
Fig. 3 (b) is the active and reactive curve of output schematic diagram of the double-fed induction wind turbine of the embodiment of the present invention.
Fig. 4 (a) is the idle circle constraint schematic diagram of the linearisation with inverter distributed generation resource of the embodiment of the present invention.
Fig. 4 (b) is the tributary capacity linearisation constraint schematic diagram of the embodiment of the present invention.
Fig. 5 is the micro-capacitance sensor topological diagram of the embodiment of the present invention.
Fig. 6 is the power prediction value schematic diagram of photovoltaic, wind-powered electricity generation, load in the embodiment of the present invention.
The micro-capacitance sensor of Fig. 7 the embodiment of the present invention in a few days operation plan mixed-integer programming model solving result schematic diagram.
Specific embodiment
With reference to embodiment, the present invention is described in further detail, but the implementation of the present invention is not limited to this.
Embodiment
The in a few days operation plan mixed-integer programming model of the micro-capacitance sensor based on analytic geometry of the present invention, as shown in Figure 1, packet Include following steps:
The object function of step 1, structure micro-capacitance sensor in a few days operation plan mixed-integer programming model is meeting normal fortune Under conditions of row constraint, by the output of each controllable of reasonable arrangement, make system operation the lowest cost, specific steps are such as Under:
It is specific as follows:
In formula, Cdeii,t、Cmtjj,tAnd Cfck,tRespectively diesel-driven generator ii, miniature gas turbine jj and fuel cell k exist The operation totle drilling cost of period t, including fuel cost, startup with shutting down cost, maintenance cost;T ∈ T represent period t in scheduling week In phase T, ii ∈ I represent that diesel-driven generator ii belongs to diesel-driven generator set I, jj ∈ J and represents that miniature gas turbine jj belongs to micro- Type gas turbine set J, k ∈ K represent that fuel cell k belongs to fuel cell set K, s ∈ S and represents that accumulator s belongs to accumulator Set S;Css,tIt is accumulator s in period t charge and discharge voltage ageing cost;Cl,tFor load l curtailing expenditure in period t,zltFor the operating status (1 represents excision, and 0 represents operation) of period t interruptible load l, pL,ltAnd Plt The unit power economic loss and active power of period t interruptible load l excision are represented respectively;Cm,tMicro- electricity during to be incorporated into the power networks Purchase sale of electricity expense between net and major network, Cm,t=pgPPCC,t, pgFor period t electricity transaction price, PPCC,tIt is handed over for period t and major network The active power changed represents power purchase for positive number, is negative number representation sale of electricity.
The linearisation constraint equation of step 2, the active and reactive output of structure diesel oil synchronous generator:
The linearisation constraint equation of the active and reactive output of the diesel oil synchronous generator is by diesel oil synchronous generator maximum Stator current restriction, diesel oil synchronous generator maximum rotor voltage restriction, diesel oil synchronous generator peak load angle Restriction, diesel oil synchronous generator minimax prime mover active power restriction composition;
Diesel oil synchronous generator maximum stator current restriction is as follows:
In formula, VsgThe voltage of node is connect by synchronous generator stator,For maximum stator current, PsgAnd QsgRespectively Active and reactive output power,For the maximum active output power of synchronous generator,For maximum rotor voltage,For Longitudinal axis reactance, PbAnd QbIt is maximum stator current constraint and the intersection point power of maximum rotor voltage constraint, PcAnd QcIt is maximum stator The intersection point power of restriction of current and prime mover maximum active power restriction, the inequality symbol of inequality constraints needs root in formula Change according to concrete numerical value and change;
Diesel oil synchronous generator maximum rotor voltage restriction is as follows:
In formula,For the minimum active output power of synchronous generator, PaAnd QaIt is the constraint of maximum rotor voltage and prime mover The intersection point power of minimum active power restriction, in formula the inequality symbol of inequality constraints need to be changed according to concrete numerical value and Variation;
Diesel oil synchronous generator peak load angle restriction is as follows:
In formula,For the maximum power angle of synchronous generator,For quadrature-axis reactance;
Diesel oil synchronous generator minimax prime mover active power restriction is as follows:
Specific derivation process is as follows:
Fig. 2 (a)~Fig. 2 (b) is the active and reactive output power curve figure of diesel oil synchronous generator, by maximum stator electricity Flow restriction C1, maximum rotor voltage restriction C2, peak load angle restriction C3, minimax prime mover wattful power Rate restriction C4Composition;
The constraint C1It is represented by:
In formula, VsgThe voltage of node is connect by synchronous generator stator,For maximum stator current, PsgAnd QsgRespectively For active and reactive output power.
The constraint C2It is represented by:
In formula,For maximum rotor voltage,For longitudinal axis reactance.
The constraint C3It is represented by:
In formula,For the maximum power angle of synchronous generator,For quadrature-axis reactance.
The constraint C4It is represented by:
In formula,Respectively synchronous generator minimum, maximum active output power.
Dash area then constrains surrounded irregular figure to meet formula (2)-(5) in Fig. 2 (a).Notice formula (2)- (3) it is constrained for quadratic inequality, is linearized, by the way that the intersection point between each function is obtained, two are replaced with the straight line of approximately equivalent Secondary inequation, active and reactive output power curve figure have been reformed into such as the polygon in Fig. 2 (b).A in Fig. 2 (b) can be obtained (Pa,Qa)、B(Pb,Qb)、C(Pc,Qc) point coordinate it is as follows:
Then 2. the linear equation in Fig. 2 (b) is 1. written as respectively:
Note that the inequality symbol in formula (9)-(10) remakes judgement after need to substituting into specific data, it is only for reference herein.
This method is substantially that will constrain (2)-(3) using the method for edge fit in circular arc to carry out piece-wise linearization, due to same It is smaller to walk generator constraint equation radian, the mode that can be directly connected under the premise of certain required precision is met with endpoint, when Computational accuracy need to consider to carry out multistage linear when being unsatisfactory for.
Step 3, the linearisation constraint equation for building the active and reactive output through grid-connected inverters distributed generation resource;It is described Include single-shaft variant miniature gas turbine, fuel cell, accumulator, photovoltaic cell through grid-connected inverters distributed generation resource;The warp The linearisation constraint equation of the active and reactive output of grid-connected inverters distributed generation resource by inverter flow through current limit constraint, The maximum active export-restriction constraint composition of distributed generation resource.
It is as follows that inverter flows through current limit constraint:
In formula, VdgThe voltage of node, I are met by inverterirTo flow through the electric current of inverter, the tune of inverter can be passed through It is whole to keep constant, Pdg、QdgThe reactive power that the active power and inverter that respectively distributed generation resource is sent out are sent out, α are inversion Device flows through the vectorial angle of the intersection point of current limit constraint and the maximum active export-restriction constraint of distributed generation resource, and θ, β are respectively nothing 2 points of vectorial angle on work(constraint curve, m linearize segments for idle circular arc, and k is integer;
The maximum active export-restriction constraint of distributed generation resource is as follows:
In formula,For the maximum output active power of distributed generation resource,Minimum output for distributed generation resource is active Power.
Specific derivation process is as follows:
Single-shaft variant miniature gas turbine, fuel cell, accumulator, photovoltaic cell after inverter commutation inversion by accessing Micro-capacitance sensor, Gu its active and reactive output is constrained by inverter.Micro- electricity is accessed as Fig. 3 (a) show source inventer with voltage The active and reactive output power curve figure of the distributed generation resource of net, influencing constraint has inverter to flow through current limit constraint C5, band The maximum active export-restriction constraint C of inverter distributed generation resource6
The constraint C5It is represented by:
In formula, VdgThe voltage of node, I are met by inverterirTo flow through the electric current of inverter, the tune of inverter can be passed through It is whole to keep constant, Pdg、QdgThe idle work(that the active power and inverter respectively sent out with inverter distributed generation resource is sent out Rate.
The constraint C6It is represented by:
In formula,For the maximum output active power of distributed generation resource,Minimum output for distributed generation resource is active Power, generally 0.
Formula (11) mathematically represents that a radius R is | Vdg|IirRound inside, therefore using analytic geometry Method linearizes the circular arc portion of Fig. 3 (a), as shown in Fig. 4 (a).Power diagram is by line segment A'B', C'D' and circular arc A' D', B'C' are surrounded, and only need to carry out piece-wise linearization to circular arc.By circular arc A'D'm deciles, each branch is linked successively, obtains circular arc Inscribed isometric broken line.As m → ∞, broken line equivalence circular arc.
If X, Y is adjacent on circular arc A'D' 2 points, the vectorial angle that D' can be obtained by the coordinate of D' points is α, then Fig. 4 (a) In β=α+(pi/2-α) k/m, θ=α+(pi/2-α) (k+1)/m, k ∈ [0, m-1], that is, put X-coordinate be (Rcos β, Rsin β), Point Y coordinate is (Rcos θ, Rsin θ), then the equation of arbitrary 1 broken line XY is represented by:
Based on analytic geometry theory, the part equation of below circular arc A'D' can be equivalent to:
(sinθ-sinβ)Pdg-(cosθ-cosβ)Qdg-sin(θ-β)R≤0 (14)
Similarly, same processing can be done to circular arc B'C'.
The linearisation constraint equation of step 4, the active and reactive output of structure double-fed induction wind turbine;The double-fed induction wind The linearisation constraint equation of the active and reactive output of machine is by wind turbine maximum stator current restriction, wind turbine maximum rotor electric current The active export-restriction constraint composition of restriction, wind turbine minimax;
Wind turbine maximum stator current restriction is as follows:
In formula, VwtThe voltage of node is connect by wind turbine,For the maximum stator current of wind turbine, PwtAnd QwtRespectively have Work(, idle output power,For the maximum active output power of wind turbine, swtFor the rotor speed of wind turbine, swt∈ (- 1,1), with the wind The fluctuation of speed and change, but it can be made for definite value by adjusting, Ph、QhIt is maximum for maximum stator current restriction and wind turbine The intersection point power of active export-restriction constraint, Pi、QiFor maximum stator current restriction and the minimum active export-restriction of wind turbine about The intersection point power of beam, the inequality symbol of inequality constraints need to be determined according to concrete numerical value in formula;
The constraint of wind turbine maximum rotor current limit is as follows:
In formula, the inequality symbol of inequality constraints need to be determined according to concrete numerical value, Pf、QfIt is limited for maximum rotor electric current Restrict beam and the intersection point power of the minimum active export-restriction constraint of wind turbine, Pg、QgFor the constraint of maximum rotor current limit and wind turbine most The intersection point power of big active export-restriction constraint,For excitation impedance,For stator impedance, a, b, c, d are multiplying for introducing The number factor, value are as follows:
In formula,For excitation resistance,For excitation reactance,For stator resistance,For stator reactance.
The active export-restriction constraint of wind turbine minimax is as follows:
Specific derivation process is as follows:
Fig. 3 (b) show the active and reactive output power curve figure of double-fed induction wind turbine, and influencing constraint has maximum stator Current limit constrains C7, maximum rotor current limit constraint C8With the maximum active export-restriction constraint C of wind turbine9
The constraint C7It is represented by:
In formula, VwtThe voltage of node is connect by wind turbine,For the maximum stator current of wind turbine, PwtAnd QwtRespectively have Work(, idle output power, swtFor the rotor speed of wind turbine, swt∈ (- 1,1), changes, but can pass through with the fluctuation of wind speed Adjusting makes it be definite value.
The constraint C8It is represented by:
In formula,For excitation impedance,For excitation resistance,For excitation reactance,For stator impedance,For stator resistance,For stator reactance;Turn for wind turbine maximum Electron current;γwtFor power of fan factor angle, it is represented by:
Due to γwtFor positive number, then sin γwtAlso it is positive number, introducing factor a, b then has:
Factor c, d are introduced, then formula (16) can be equivalent to:
The maximum active output of wind turbine is expressed as:
In formula,For the maximum active output power of wind turbine.
Each circular arc is directly connected into straight line equivalence using intersection point with as synchronous generator power curve linear method, can be asked Go out Fig. 3 (b) midpoints F (Pf,Qf)、G(Pg,Qg)、H(Ph,Qh)、I(Pi,Qi) coordinate it is as follows:
Similarly constrain C7With constraint C8It approximately linear can turn to respectively:
Note that the inequality symbol in formula (25)-(26) remakes judgement after need to substituting into specific data, it is only for reference herein.
Step 5, structure consider that the linearisation power flow equation of voltage magnitude is as follows:
In formula, PGiAnd QGiThe respectively generated power and reactive power of node i, PDiAnd QDiThe respectively load of node i Active and reactive power;ViAnd VjThe respectively voltage magnitude of node i and j;J ∈ i represent that node j is connected directly with node i, θij、GijAnd BijRespectively node i and phase angle difference, conductance and the susceptance of node j, GiiAnd BiiThe respectively reality of node i self-admittance Portion and imaginary part, PijAnd QijThe respectively active and reactive power of branch ij;Gi0And Bi0Respectively branch ij interior joints i is over the ground Conductance and susceptance;Using special sort set pairMake linearization process.
Specific derivation process is as follows:
Micro-capacitance sensor in a few days operation plan substantially consider operational safety optimal load flow model.Trend constraint is in a few days adjusts The equality constraint of planning model is spent, using based on polar tide model:
In formula, PGiAnd QGiThe respectively generated power and reactive power of node i, PDiAnd QDiThe respectively load of node i Active and reactive power;ViAnd VjThe respectively voltage magnitude of node i and j;J ∈ i represent that node j is connected directly with node i, θij、GijAnd BijRespectively node i and phase angle difference, conductance and the susceptance of node j, GiiAnd BiiThe respectively reality of node i self-admittance Portion and imaginary part, PijAnd QijThe respectively active and reactive power of branch ij;Gi0And Bi0Respectively branch ij interior joints i is over the ground Conductance and susceptance.
It improves power quality, improve the important content that voltage level is not only micro-capacitance sensor economic load dispatching, also relationship micro-capacitance sensor The operation characteristic of middle micro battery and load.Therefore it is presumed that the DC flow model that voltage magnitude is 1.0pu is not directly applicable Micro-capacitance sensor economic load dispatching.For this purpose, the monomial linear approximation processing one by one formed to power flow equation variable.Micro-capacitance sensor radius of electricity supply It is small, meet 0.95p.u.≤Vi≤ 1.05p.u., | θij|≤40o conditions, state variable form monomial original value and approximation Maximum deviation, as shown in table 1.The maximum deviation of power flow equation monomial and approximate expression is within 3% as seen from Table 1, and passes The cos θ of systemij1.0 worst errors of ≈ are up to 23.4%.With reference to trigonometric function formulaThen formula (27) it can be turned to formula (28):
The error of 1 power flow equation monomial of table and approximate expression
It noticesBe not linear equation there are formula (31) and formula (32), using special sort set (special Order set) it is rightMake linearization process.While formula (31) and formula (32) are converted into linear equation, also retain voltage magnitude, Be conducive to consider the influence that voltage runs micro-capacitance sensor.
Step 6, structure tributary capacity linearisation constraint equation are as follows:
In formula, SijRepresent the transmission capacity limits value of branch ij, α ', β ' are tributary capacity restrained circleOn 2 points of vectorial angle, m' are the equilateral number of edges of inscribed, and k' is integer.
Specific derivation process is as follows:
The transmission capacity of branch ij must meet:
In formula, SijRepresent the transmission capacity limits value of branch ij.Analytically geometric angle is seen, formula (33) represents point (Pij, Qij) in radius be SijCircle in, as shown in Fig. 4 (b).If by circumference m' deciles, each branch is linked successively, obtains inscribed Wait m' sides shape.As m' → ∞, point (Pij,Qij) point (P is equivalent in circleij,Qij) in the shape of the m' such as inscribed sides.
If A, B is adjacent 2 points of the m' sides shape such as inscribed, the vectorial angle of A, B is respectively α ', β ', then the π k'/m' of α '=2, The π (k'+1) of β '=2/m', k' ∈ [0, m'-1], i.e. point A coordinates are (SijCos α ', SijSin α '), point B coordinates are (SijCos β ', SijSin β '), then as m', when shape is 1 arbitrary, AB linear equations are represented by:
Based on analytic geometry theory, formula (33) can turn to:
(sinβ'-sinα')Pij-(cosβ'-cosα')Qij-sin(β'-α')Sij≤0 (35)
Step 7 is constrained with reference to the existing linearisation of micro-capacitance sensor, forms micro-capacitance sensor in a few days operation plan mixed integer programming mould Type solves software with mixed integer programming and is solved.
The present embodiment in a few days dispatches the micro-capacitance sensor based on analytic geometry of the present invention by taking microgrid figure shown in fig. 5 as an example The idiographic flow of plan mixed-integer programming model illustrates:
Microgrid figure shown in fig. 5 by power supply, load, circuit, bus nodes, reactive power compensator, major network and its with public affairs Transformer composition between tie point (point of common couple, PCC) altogether.Wherein, power supply includes 1 miniature diesel Machine, 1 asynchronous wind turbine, 1 photovoltaic generator and 2 accumulators (including the 1st accumulator and the 2nd accumulator);Load includes 1 Important load, 5 interruptible loads (including the 1st~5 interruptible load) and 1 elastic load;Circuit has 8 bars of (the 1st~8 lines Road), model is LJ_95, wherein the length of the 1st circuit is 500m, the length of the 2nd circuit is 100m, the length of the 3rd circuit For 50m, the length of the 4th circuit is 150m, and the length of the 5th circuit is 120m, and the length of the 6th circuit is 100m, the 7th circuit Length is 200m, the length of the 8th circuit is 150m,;Bus nodes have 10, are respectively designated as A~J;Reactive power compensator connects On bus nodes B;The voltage class of major network is 110kV, and the voltage class of microgrid is 10kV.
It is 0.4kV, power reference value 100kVA to take micro-capacitance sensor voltage a reference value.In a few days the time span of operation plan is 4h is divided into 16 periods, per the period a length of 15min.Specifically choose certain day 14:00-18:00, photovoltaic, the wind-powered electricity generation of the time range And load power short-term forecast value is as shown in Figure 6.
The idiographic flow of model construction is as follows:
(1) object function of structure micro-capacitance sensor in a few days operation plan mixed-integer programming model, the present embodiment fuel unit Using diesel engine as representative, Gu involved operating cost is taken including diesel engine operating cost, accumulator aging cost, load reduction With, be incorporated into the power networks when micro-capacitance sensor and major network between purchase sale of electricity expense.
The object function of the micro-capacitance sensor in a few days operation plan is:
(2) the linearisation constraint equation of the active and reactive output of structure diesel oil synchronous generator, diesel oil synchronous generator Related operating parameter (perunit value) such as table 2:
The related operating parameter of 2 diesel oil synchronous generator of table
The linearisation constraint equation of structure is as follows:
(3) build accumulator, photovoltaic cell active and reactive output linearisation constraint equation, accumulator and photovoltaic electric Related operating parameter (perunit value) such as table 3 in pond:
3 accumulator of table operating parameter related to photovoltaic cell
It is 1 to take the hop count m of segmentation, the storage battery active power of structure, it is idle go out force linearizing constraint equation it is as follows:
The photovoltaic cell of structure it is active and reactive go out force linearizing constraint equation it is as follows:
(4) the linearisation constraint equation of the active and reactive output of structure double-fed induction wind turbine, the correlation of double-fed induction wind turbine Operating parameter (perunit value) such as table 4:
The related operating parameter of 4 double-fed induction wind turbine of table
The double-fed induction wind turbine of structure it is active and reactive go out force linearizing constraint equation it is as follows:
(5) structure considers the linearisation power flow equation of voltage magnitude, and the resistance of circuit is 0.069 Ω/km, and reactance is 0.099 Ω/km, ignores branch admittance over the ground, and G is calculated with reference to Fig. 5 dataii、Bii, it is as follows to list linearisation power flow equation:
It is right that special sort set (special order set) can be usedMake linearization process.
(6) structure tributary capacity linearisation constraint equation, branch maximum capacity SijFor 100kW, it is as follows to establish equation:
(sinβ'-sinα')Pij-(cosβ'-cosα')Qij-100sin(β'-α')≤0 (42)
(7) with reference to other existing linearisation constraints of micro-capacitance sensor, micro-capacitance sensor in a few days operation plan mixed integer programming is formed Model solves software with mixed integer programming and is solved, and the in a few days active plan of each distributed generation resource is as shown in Figure 7.
Above-described embodiment is the preferable embodiment of the present invention, but embodiments of the present invention are not by the embodiment Limitation, other any Spirit Essences without departing from the present invention with made under principle change, modification, replacement, combine, simplification, Equivalent substitute mode is should be, is included within protection scope of the present invention.

Claims (1)

1. the modeling method of micro-capacitance sensor in a few days operation plan mixed-integer programming model, which is characterized in that include the following steps:
The object function of step 1, structure micro-capacitance sensor in a few days operation plan mixed-integer programming model;
The object function of the micro-capacitance sensor in a few days operation plan mixed-integer programming model, it is specific as follows:
In formula, Cdeii,t、Cmtjj,tAnd Cfck,tRespectively diesel-driven generator ii, miniature gas turbine jj and fuel cell k are in period t Operation totle drilling cost, including fuel cost, startup with shut down cost, maintenance cost;T ∈ T represent period t in dispatching cycle T, Ii ∈ I represent that diesel-driven generator ii belongs to diesel-driven generator set I, jj ∈ J and represents that miniature gas turbine jj belongs to miniature gas Turbine set J, k ∈ K represent that fuel cell k belongs to fuel cell set K, s ∈ S and represents that accumulator s belongs to accumulator set S; CSs, tIt is accumulator s in period t charge and discharge voltage ageing cost;Cl,tFor load l curtailing expenditure in period t,zltFor the operating status of period t interruptible load l, pL,ltAnd PltIt is negative to represent that period t can interrupt respectively The unit power economic loss and active power of lotus l excisions;Cm,tPurchase sale of electricity expense during to be incorporated into the power networks between micro-capacitance sensor and major network With Cm,t=pgPPCC,t, pgFor period t electricity transaction price, PPCC,tIt is positive number for the active power that period t is exchanged with major network It represents power purchase, is negative number representation sale of electricity;
The linearisation constraint equation of step 2, the active and reactive output of structure diesel oil synchronous generator;
The linearisation constraint equation of the active and reactive output of the diesel oil synchronous generator is by diesel oil synchronous generator maximum stator Current limit constraint, diesel oil synchronous generator maximum rotor voltage restriction, the limitation of diesel oil synchronous generator peak load angle Constraint, diesel oil synchronous generator minimax prime mover active power restriction composition;
Diesel oil synchronous generator maximum stator current restriction is as follows:
In formula, VsgThe voltage of node is connect by synchronous generator stator,For maximum stator current, PsgAnd QsgRespectively have Work(, idle output power,For the maximum active output power of synchronous generator,For maximum rotor voltage,For the longitudinal axis Reactance, PbAnd QbIt is maximum stator current constraint and the intersection point power of maximum rotor voltage constraint, PcAnd QcIt is maximum stator current Constraint and the intersection point power of prime mover maximum active power restriction, the inequality symbol of inequality constraints need to be according to tool in formula Body numerical value change and change;
Diesel oil synchronous generator maximum rotor voltage restriction is as follows:
In formula,For the minimum active output power of synchronous generator, PaAnd QaIt is that the constraint of maximum rotor voltage and prime mover are minimum The intersection point power of active power restriction, the inequality symbol of inequality constraints need to change according to concrete numerical value and be become in formula Change;
Diesel oil synchronous generator peak load angle restriction is as follows:
In formula,For the maximum power angle of synchronous generator,For quadrature-axis reactance;
Diesel oil synchronous generator minimax prime mover active power restriction is as follows:
Step 3, the linearisation constraint equation for building the active and reactive output through grid-connected inverters distributed generation resource;It is described through inverse Become the grid-connected distributed generation resource of device and include single-shaft variant miniature gas turbine, fuel cell, accumulator, photovoltaic cell;
The linearisation constraint equation of the active and reactive output through grid-connected inverters distributed generation resource flows through electric current by inverter The maximum active export-restriction constraint composition of restriction, distributed generation resource;
It is as follows that inverter flows through current limit constraint:
In formula, VdgThe voltage of node, I are met by inverterirTo flow through the electric current of inverter, can be kept by the adjustment of inverter It is constant, Pdg、QdgThe reactive power that the active power and inverter that respectively distributed generation resource is sent out are sent out, α are flowed through for inverter Current limit constrains and the vectorial angle of the intersection point of the maximum active export-restriction constraint of distributed generation resource, and θ, β are respectively idle constraint 2 points of vectorial angle on curve, m linearize segments for idle circular arc, and k is integer;
The maximum active export-restriction constraint of distributed generation resource is as follows:
In formula,For the maximum output active power of distributed generation resource,Minimum active power of output for distributed generation resource;
The linearisation constraint equation of step 4, the active and reactive output of structure double-fed induction wind turbine;
The linearisation constraint equation of the active and reactive output of the double-fed induction wind turbine is limited about by wind turbine maximum stator current Beam, the constraint of wind turbine maximum rotor current limit, the active export-restriction constraint composition of wind turbine minimax;
Wind turbine maximum stator current restriction is as follows:
In formula, VwtThe voltage of node is connect by wind turbine,For the maximum stator current of wind turbine, PwtAnd QwtRespectively active, nothing Work(output power,For the maximum active output power of wind turbine, swtFor the rotor speed of wind turbine, swt∈ (- 1,1), with wind speed It fluctuates and changes, but it can be made for definite value by adjusting, Ph、QhIt is active for maximum stator current restriction and wind turbine maximum The intersection point power of export-restriction constraint, Pi、QiIt is constrained for maximum stator current restriction and the minimum active export-restriction of wind turbine Intersection point power, the inequality symbol of inequality constraints need to be determined according to concrete numerical value in formula;
The constraint of wind turbine maximum rotor current limit is as follows:
In formula, the inequality symbol of inequality constraints need to be determined according to concrete numerical value, Pf、QfFor maximum rotor current limit about Beam and the intersection point power of the minimum active export-restriction constraint of wind turbine, Pg、QgHave for the constraint of maximum rotor current limit and wind turbine maximum The intersection point power of work(export-restriction constraint,For excitation impedance,For stator impedance, a, b, c, d are the multiplier factor introduced, Its value is as follows:
In formula,For excitation resistance,For excitation reactance,For stator resistance,For stator reactance,
The active export-restriction constraint of wind turbine minimax is as follows:
Step 5, structure consider the linearisation power flow equation of voltage magnitude;
The linearisation power flow equation for considering voltage magnitude is as follows:
In formula, PGiAnd QGiThe respectively generated power and reactive power of node i, PDiAnd QDiThe load of respectively node i is active And reactive power;ViAnd VjThe respectively voltage magnitude of node i and j;J ∈ i represent that node j is connected directly with node i, θij、Gij And BijRespectively node i and phase angle difference, conductance and the susceptance of node j, GiiAnd BiiThe respectively real part and void of node i self-admittance Portion, PijAnd QijThe respectively active and reactive power of branch ij;Gi0And Bi0Respectively the conductance over the ground of branch ij interior joints i and Susceptance;
Using special sort set pairMake linearization process;
Step 6, structure tributary capacity linearisation constraint equation;
The tributary capacity linearisation constraint equation is as follows:
In formula, SijRepresent the transmission capacity limits value of branch ij, α ', β ' are tributary capacity restrained circleUpper 2 points Vectorial angle, PijAnd QijThe respectively active and reactive power of branch ij, m' are the equilateral number of edges of inscribed, and k' is integer;
Step 7 is constrained with reference to the existing linearisation of micro-capacitance sensor, forms micro-capacitance sensor in a few days operation plan mixed-integer programming model, Software is solved with mixed integer programming to be solved.
CN201510220659.2A 2015-04-30 2015-04-30 The modeling method of micro-capacitance sensor in a few days operation plan mixed-integer programming model Expired - Fee Related CN104835080B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510220659.2A CN104835080B (en) 2015-04-30 2015-04-30 The modeling method of micro-capacitance sensor in a few days operation plan mixed-integer programming model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510220659.2A CN104835080B (en) 2015-04-30 2015-04-30 The modeling method of micro-capacitance sensor in a few days operation plan mixed-integer programming model

Publications (2)

Publication Number Publication Date
CN104835080A CN104835080A (en) 2015-08-12
CN104835080B true CN104835080B (en) 2018-06-22

Family

ID=53812950

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510220659.2A Expired - Fee Related CN104835080B (en) 2015-04-30 2015-04-30 The modeling method of micro-capacitance sensor in a few days operation plan mixed-integer programming model

Country Status (1)

Country Link
CN (1) CN104835080B (en)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105870912A (en) * 2016-03-03 2016-08-17 广州华工弈高科技有限公司 Modeling method for island-division mixed integer programming model for active power distribution network
CN105955931B (en) * 2016-05-10 2019-02-19 东南大学 Regional Energy network optimization dispatching method towards the consumption of high density distributed photovoltaic
CN106026074A (en) * 2016-05-25 2016-10-12 安徽大学 Low-voltage DC micro-grid experiment platform
CN106712032B (en) * 2017-02-23 2018-12-25 国家电网有限公司 A kind of optimal load flow model building method considering Wind turbines active voltage regulating power
CN108173265B (en) * 2018-01-15 2020-09-25 国网吉林省电力有限公司电力科学研究院 Power distribution network reconstruction method based on linearized power flow
CN108509378B (en) * 2018-02-11 2021-07-06 哈尔滨理工大学 Parameter matching method for composite energy storage system in electric vehicle test platform
CN108448636B (en) * 2018-05-10 2021-02-05 合肥工业大学 Alternating current-direct current hybrid micro-grid optimization planning method considering line factors
CN112818498A (en) * 2019-10-29 2021-05-18 新奥数能科技有限公司 Energy system modeling optimization method and device
CN111130118B (en) * 2020-01-09 2021-02-02 清华大学 Optimal power flow calculation method of power system based on piecewise linearization
CN111641226B (en) * 2020-05-19 2021-12-17 浙江工业大学 Building type micro-grid photovoltaic utilization rate improvement method considering automatic demand response
CN111884266A (en) * 2020-07-02 2020-11-03 河海大学 Gas turbine intraday rolling unit combination optimization method
CN114167131B (en) * 2021-12-06 2023-07-18 国网湖南省电力有限公司 Energy storage system power envelope diagram test method and system based on dynamic characteristic adjustment

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7853911B1 (en) * 2005-11-04 2010-12-14 Altera Corporation Method and apparatus for performing path-level skew optimization and analysis for a logic design
CN102104251A (en) * 2011-02-24 2011-06-22 浙江大学 Microgrid real-time energy optimizing and scheduling method in parallel running mode
CN102184475A (en) * 2011-05-11 2011-09-14 浙江大学 Optimizing and dispatching method for microgrid economical operation on basis of multiple time scale coordination
CN102751728A (en) * 2012-07-26 2012-10-24 浙江大学 Energy management method for isolated network running mode in micro network based on load interruption model

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9014868B2 (en) * 2012-03-29 2015-04-21 International Business Machines Corporation Power factor
US10296988B2 (en) * 2013-08-19 2019-05-21 Board Of Trustees Of Michigan State University Linear optimal power flow system and method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7853911B1 (en) * 2005-11-04 2010-12-14 Altera Corporation Method and apparatus for performing path-level skew optimization and analysis for a logic design
CN102104251A (en) * 2011-02-24 2011-06-22 浙江大学 Microgrid real-time energy optimizing and scheduling method in parallel running mode
CN102184475A (en) * 2011-05-11 2011-09-14 浙江大学 Optimizing and dispatching method for microgrid economical operation on basis of multiple time scale coordination
CN102751728A (en) * 2012-07-26 2012-10-24 浙江大学 Energy management method for isolated network running mode in micro network based on load interruption model

Also Published As

Publication number Publication date
CN104835080A (en) 2015-08-12

Similar Documents

Publication Publication Date Title
CN104835080B (en) The modeling method of micro-capacitance sensor in a few days operation plan mixed-integer programming model
Khosravi et al. A novel control approach to improve the stability of hybrid AC/DC microgrids
CN102751728B (en) Energy management method for isolated network running mode in micro network based on load interruption model
Wang et al. Power-flow control and stability enhancement of four parallel-operated offshore wind farms using a line-commutated HVDC link
CN106487042B (en) A kind of Multiple Time Scales micro-capacitance sensor voltage power-less optimized controlling method
CN104092250B (en) The distributed economic load dispatching of micro-grid system and control method for coordinating
CN102856899B (en) Method of reducing network loss of micro power grid
CN110206686A (en) A kind of adaptive maximum power tracking and controlling method for wind power generating set
CN106026113A (en) Micro-grid system monitoring method having reactive automatic compensation function
CN103593711B (en) A kind of distributed power source Optimal Configuration Method
CN103219732A (en) reactive voltage controlling method of power distribution network with variable speed constant frequency wind farm
Rekik et al. A micro-grid ensuring multi-objective control strategy of a power electrical system for quality improvement
Ishraque et al. Optimized design of a hybrid microgrid using renewable resources considering different dispatch strategies
El-Bahay et al. Computational methods to mitigate the effect of high penetration of renewable energy sources on power system frequency regulation: a comprehensive review
CN112994051A (en) Capacity optimization configuration method for seawater pumped storage power station with variable-speed unit
CN106780119A (en) Power distribution network regenerative resource based on many active managements is dissolved method
Tadjeddine et al. VRE integrating in Piat grid with AFRR using PSS, MPPT, and PSO-based techniques: A case study kabertene
CN106786734A (en) Power grid wind accesses site selecting method
Yao et al. Active power dispatch strategy of the wind farm based on improved multi‐agent consistency algorithm
CN116316666A (en) Reactive power compensation device and energy storage coordination optimization configuration method and system
Gupta et al. Narendra kumar,“
CN105262107A (en) Wind power grid connected reactive optimization method
Azis et al. Fuzzy logic based active power generation dispatching considering intermittent wind power plants output
Zhang et al. Hierarchical coordinated control method for multiload DC microgrid units
Blaabjerg et al. Power electronics in wind generation systems

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180622

Termination date: 20190430