CN106329537B - A kind of idle work optimization method adapting to bulk power grid automatism voltage control - Google Patents
A kind of idle work optimization method adapting to bulk power grid automatism voltage control Download PDFInfo
- Publication number
- CN106329537B CN106329537B CN201510341645.6A CN201510341645A CN106329537B CN 106329537 B CN106329537 B CN 106329537B CN 201510341645 A CN201510341645 A CN 201510341645A CN 106329537 B CN106329537 B CN 106329537B
- Authority
- CN
- China
- Prior art keywords
- indicate
- node
- voltage
- idle work
- work optimization
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/30—Reactive power compensation
Abstract
The present invention provides a kind of idle work optimization method adapting to bulk power grid automatism voltage control, includes the following steps：Determine the voltage stability index of branch；Establish idle work optimization model；Idle work optimization model is solved using interior point method.For there is voltage stabilization in power grid, and since the generation of Voltage-stabilizing Problems usually has locality, using analysis branch voltage steady stability index, the voltage stability margin in this region is improved by targetedly reducing the branch voltage stability index numerical value in certain region in Optimized model object function.
Description
Technical field
The invention belongs to technical field of electric power automation, and in particular to a kind of to adapt to the idle of bulk power grid automatism voltage control
Optimization method.
Background technology
The core of the automatism voltage control of bulk power grid is that the idle work optimization of electric system calculates.Idle work optimization is optimal load flow
In a typical problem.Traditional bulk power grid automatism voltage control is under given trend section, i.e. the network of system is opened up
It flutters when thinking fixed with parameter, active and idle, the generated power output of load, is meeting system operation putting equipment in service condition
In the case of system status parameters range of operation, to generator terminal voltage, reactive-load compensation equipment throwing is moved back and adjustable transformer
Gear adjust, change the reactive power distribution of system with this, to reduce the network loss of system.
The Reactive Power Optimazation Problem of extensive practical power systems at present, widely used is that nonlinear interior-point method carries out directly
Solve and pass through two kinds of technology paths of district and grade control.The method of degree and zoning stresses the validity to power grid control, in electricity
It is weaker in net analysis.Due to subregion, controlled quentity controlled variable and quantity of state and controlled quentity controlled variable and state outside region between setting area by force
Between amount without directly affect either weak influence but it is interregional in practice be interconnection or the control for adhering to two regions separately of interconnection region
Amount processed and quantity of state are also strong correlation, therefore, do not ensure that its control direction is directed towards optimal direction adjusting sometimes.It is non-thread
Property interior point method stress the direct optimization to power grid analysis, due to directly carry power flow equation constrain, generate control strategy
It is also to meet power flow equation to be.However the data that provide of appearance due to discrete device and state estimation occur deviation it is larger when
The problem of will appear solution difficulty and being difficult to solve, therefore, the prior art also needs to be improved.
Invention content
In order to overcome the above-mentioned deficiencies of the prior art, the present invention provide it is a kind of adapt to bulk power grid automatism voltage control it is idle
Optimization method specifically adopts the following technical scheme that：
The present invention provides a kind of idle work optimization method adapting to bulk power grid automatism voltage control, and the method includes following steps
Suddenly：
Step 1：Determine the voltage stability index of branch；
Step 2：Establish idle work optimization model；
Step 3：Idle work optimization model is solved using interior point method.
The step 1 specifically includes following steps：
Step 1-1：The voltage stability margin L of branch ij between calculate node i and node jij, as the electricity of branch ij
Press stability index；
Step 1-2：Judge LijWhether it is more than threshold value T, if more than respective branch is contributed to set of fingers M.
In the step 1-1, the voltage stability margin of branch ij is L between node i and node jijIt is expressed as：
Wherein, R and X indicates resistance and the reactance of branch ij, δijIndicate the phase difference of voltage of node i and node j, UiIt indicates
The voltage magnitude of node i, PjAnd QjThe outflow active power and reactive power of node j are indicated respectively；
In the step 2, idle work optimization model includes idle work optimization object function and idle work optimization constraints.
The idle work optimization constraints includes equality constraint and inequality constraints；
The equality constraint includes trend equality constraint, transformer voltage ratio equality constraint and capacitive reactance device equality constraint；
The inequality constraints includes node voltage inequality constraints, the inequality constraints of voltage slack and generator reactive
Output inequality constraints.
The idle work optimization object function is expressed as：
Wherein, f (x) indicates that idle work optimization object function, M indicate set of fingers, UiAnd UjNode i and node j are indicated respectively
Voltage magnitude, GijIndicate the conductance of branch ij, δijIndicate that the phase difference of voltage of node i and node j, N indicate node set, m
Indicate that voltage stability margin is more than the branch of threshold value；SiIndicate the voltage slack of node i,Indicate that the voltage of branch ij is steady
Determine nargin, ω indicates that variation weight, ψ indicate voltage stability index weight.
In the equality constraint, have：
(1) trend equality constraint is expressed as：
Wherein, PG,iAnd QG,iGenerator injection active power and reactive power, P at node i are indicated respectivelyD,iAnd QD,iRespectively
Indicate generated power load and load or burden without work at node i, UiAnd UjThe voltage magnitude of node i and node j, G are indicated respectivelyijWith
BijIndicate that the conductance and susceptance of branch ij, N indicate node set, δ respectivelyijFor the phase difference of voltage of node i and node j, and δij
=θi-θj, wherein θi、θjRespectively node i, j voltage phase angle；
(2) transformer voltage ratio equality constraint is expressed as：
(Tk-Tk,min)(Tk,max-TkThe k ∈ S of)=0T (4)
Wherein, TkIndicate the no-load voltage ratio of k-th of transformer, Tk,maxAnd Tk,minAbove and below the no-load voltage ratio for indicating k-th of transformer respectively
Limit, STIndication transformer set；
(3) capacitive reactance device equality constraint is expressed as：
(Bh-Bh,min)(Bh,max-BhThe h ∈ S of)=0B (5)
Wherein, BhIndicate the susceptance of h-th of capacitive reactance device, Bh,maxAnd Bh,minAbove and below the susceptance for indicating h-th of capacitive reactance device respectively
Limit, SBIndicate capacitive reactance device set.
In the inequality constraints, have：
(1) node voltage inequality constraints is expressed as：
Ui,min+Si≤Ui≤Ui,max-Si i∈N (7)
Wherein, SiIndicate the voltage slack of node i, UiIndicate the voltage magnitude of node i, Ui,maxAnd Ui,minIt indicates respectively
The voltage magnitude bound of node i, N indicate node set；
(2) inequality constraints of voltage slack is expressed as：
Si≥0 i∈N (8)
Wherein, SiIndicate that the voltage slack of node i, N indicate node set；
(3) inequality constraints of generator reactive output is expressed as：
QGi,min≤QG,i≤QGi,max i∈N (9)
Wherein, QG,iIndicate that generator injects reactive power, Q at node i respectivelyGi,maxAnd QGi,minIt indicates at node i respectively
Generator injects reactive power bound, and N indicates node set.
The step 3 specifically includes following steps：
Step 3-1：Idle work optimization model is divided into discrete device equality constraint and continuous quantity Optimized model first；
Step 3-2：Continuous quantity Optimized model is expressed as：
Wherein, f (x) indicates that idle work optimization object function, M indicate set of fingers, UiAnd UjNode i and node j are indicated respectively
Voltage magnitude, GijIndicate the conductance of branch ij, δijIndicate that the phase difference of voltage of node i and node j, N indicate node set, m
Indicate that voltage stability margin is more than the branch of threshold value；SiIndicate the voltage slack of node i,Indicate that the voltage of branch ij is steady
Determine nargin, ω indicates that variation weight, ψ indicate voltage stability index weight, Bh,iIndicate the electricity of h-th of capacitive reactance device at node i
It receives, Tk,iIndicate the no-load voltage ratio of k-th of transformer at node i, Δ Pi、ΔQiActive deviation and idle deviation at node i are indicated respectively；
Definition vector x=[QG,1,…,QG,i,…,QG,N,S1,…,Si,…,SN,U1,…,Ui,…,UN,θ1,…,θi,…,
θN]T, unconfined majorized function L is formed for continuous quantity Optimized model, is had：
Wherein, SGIndicate that generator node, μ indicate disturbance variable, and μ >=0；li、uiUnder the voltage for indicating node i respectively
Slack variable in slack variable and voltage, ypi、yqi、zi、wiIndicate the lagrange's variable of node i, definition vector l=[l1,
l2,……,lr], u=[u1,u2,……,ur], vectorial y={ yp1,…ypi,…,ypN,yq1,…,yqi,…,yqN, vectorial z=
[z1,z2,……,zr], vectorial w=[w1,w2,……,wr], r expression inequality numbers；
Then meet：
Wherein, LxIndicate majorized function L to the local derviation of vector x, LyIndicate majorized function L to the local derviation of vectorial y, LzIt indicates
Majorized function L is to the local derviation of vectorial z, LwIndicate majorized function L to the local derviation of vectorial w, LlIndicate majorized function L to the inclined of vectorial l
It leads, LuIndicate majorized function L to the local derviation of vectorial u, ▽xIndicate idle work optimization object function f (x) to the local derviation of vector x, trend
Equality constraint vector h (x)={ Δ P1,…,ΔPi,…,ΔPN, Δ Q1,…,ΔQi,…,ΔQN, g (x) indicates inequality about
Shu Xiangliang,Indicate the inequality constraints upper limit,gIndicate inequality constraints lower limit, L, Z, U, W are respectively li、zi、ui、wiThe r of composition
Diagonal matrix is tieed up, the r that E is ties up unit matrix；
Step 3-3：The discrete device equality constraint includes transformer voltage ratio equality constraint and capacitive reactance device equality constraint；
Increase duality gapThen have：
Wherein, TkIndicate the no-load voltage ratio of k-th of transformer, Tk,maxAnd Tk,minAbove and below the no-load voltage ratio for indicating k-th of transformer respectively
Limit, STIndication transformer set；BhIndicate the susceptance of h-th of capacitive reactance device, Bh,maxAnd Bh,minThe electricity of h-th of capacitive reactance device is indicated respectively
Receive bound, SBIndicate capacitive reactance device set；
Step 3-4：Idle work optimization model is solved, is specifically included：
Step 1)：Initialization vector l, z, u, w, y, each vector kind of element is respectively 0.5, -0.2,0.5,0.2, -0.01,
It is 0.99 to initialize centripetal parameter Sigma；
Step 2)：Calculate Cgap、Lx、Ly、Lz、Lw、Ll、LuIf calculated Cgap、Lx、Ly、Lz、Lw、Ll、LuRespectively less than count
Precision ε is calculated, then stops calculating；
Step 3)：According to obtained CgapCalculation perturbation variable
Step 4)：Material calculation coefficient stepPAnd stepD, have：
Step 5)：According to stepPAnd stepDX, l, u, y, z, w are corrected, is had：
Wherein, Δ x, Δ l, Δ u, Δ y, Δ z, Δ w indicate respectively vector x, l, u, y, z, w iteration step length；
Step 6)：T is corrected by solving equation (13)kAnd Bh, have：
Tk=Tk+ΔTk (18)
Bh=Bh+ΔBh (19)
Wherein, Δ TkWith Δ BhTo correct step-length；
Step 7)：By revised Tk、Bh, after x, l, u, y, z, w bring equation (12) and (13) into, return to step 2).
Compared with prior art, the beneficial effects of the present invention are：
1) according to voltage in actual motion be soft-constraint the characteristics of, idle work optimization is improved using the method for voltage relaxationization
The convergence of computational methods；
2) it is converted into equality constraint by the way that discrete device is adjusted constraint, improves the efficiency of single optimization, can adapt to
The optimization analysis of a large amount of discrete devices control in bulk power grid, and propose corresponding optimization method for equality constraint；
3) by by LijIndex converts Pj(R+X)+Qj(X-R) form has needle under conditions of to increase model nonconvex property
The low voltage stability of voltage stability margin is improved to property, to achieve the purpose that improve system voltage stability.
Description of the drawings
Fig. 1 is the idle work optimization method flow chart that bulk power grid automatism voltage control is adapted in the embodiment of the present invention.
Specific implementation mode
Invention is further described in detail below in conjunction with the accompanying drawings.
Traditional nonlinear interior-point method calculates idle work optimization and is susceptible to due to not restrained caused by data reasons, for from
It is casually arranged with for needing regular nearby or being caused to calculate inaccurately using intelligent algorithm, for this three classes of the voltage stability of consideration system
Situation.The present invention is by mathematical distortions and considers voltage operation characteristic, compares traditional nonlinear interior-point method idle work optimization, this Shen
The characteristics of being please soft-constraint according to voltage in actual motion, idle work optimization computational methods are improved using the method for voltage relaxationization
Convergence；It is converted into equality constraint by the way that discrete device is adjusted constraint, improves the efficiency of single optimization, can adapt to big
The optimization analysis of a large amount of discrete devices control in power grid, and propose corresponding optimization method for equality constraint；Passing through will
LijIndex converts Pj(R+X)+Qj(X-R) it is steady targetedly to improve voltage under conditions of to increase model nonconvex property for form
The low voltage stability of nargin is determined, to achieve the purpose that improve system voltage stability.
Such as Fig. 1, the present invention provides a kind of idle work optimization method adapting to bulk power grid automatism voltage control, the method includes
Following steps：
Step 1：Determine the voltage stability index of branch；
Step 2：Establish idle work optimization model；
Step 3：Idle work optimization model is solved using interior point method.
The step 1 specifically includes following steps：
Step 1-1：The voltage stability margin L of branch ij between calculate node i and node jij, as the electricity of branch ij
Press stability index；
Step 1-2：Judge LijWhether it is more than threshold value T, if more than respective branch is contributed to set of fingers M.
In the step 1-1, the voltage stability margin of branch ij is L between node i and node jijIt is expressed as：
Wherein, R and X indicates resistance and the reactance of branch ij, δijIndicate the phase difference of voltage of node i and node j, UiIt indicates
The voltage magnitude of node i, PjAnd QjThe outflow active power and reactive power of node j are indicated respectively；
Certain actual electric network result of calculation such as table 1：
Table 1
| Circuit name | Voltage stability index (threshold value is set as 0.6) |
| Mountain yellow line | 0.650653 |
| Other white line | 0.469824 |
| Platform tower line | 0.397849 |
| Black multi-color cord | 0.396 |
| Xu Lu lines | 0.380926 |
In the step 2, idle work optimization model includes idle work optimization object function and idle work optimization constraints.
The idle work optimization constraints includes equality constraint and inequality constraints；
The equality constraint includes trend equality constraint, transformer voltage ratio equality constraint and capacitive reactance device equality constraint；
The inequality constraints includes node voltage inequality constraints, the inequality constraints of voltage slack and generator reactive
Output inequality constraints.
The idle work optimization object function is expressed as：
Wherein, f (x) indicates that idle work optimization object function, M indicate set of fingers, UiAnd UjNode i and node j are indicated respectively
Voltage magnitude, GijIndicate the conductance of branch ij, δijIndicate that the phase difference of voltage of node i and node j, N indicate node set, m
Indicate that voltage stability margin is more than the branch of threshold value；SiIndicate the voltage slack of node i,Indicate that the voltage of branch ij is steady
Determine nargin, ω indicates that variation weight, ψ indicate voltage stability index weight.
In the equality constraint, have：
(1) trend equality constraint is expressed as：
Wherein, PG,iAnd QG,iGenerator injection active power and reactive power, P at node i are indicated respectivelyD,iAnd QD,iRespectively
Indicate generated power load and load or burden without work at node i, UiAnd UjThe voltage magnitude of node i and node j, G are indicated respectivelyijWith
BijIndicate that the conductance and susceptance of branch ij, N indicate node set, δ respectivelyijFor the phase difference of voltage of node i and node j, and δij
=θi-θj, wherein θi、θjRespectively node i, j voltage phase angle；
(2) transformer voltage ratio equality constraint is expressed as：
(Tk-Tk,min)(Tk,max-TkThe k ∈ S of)=0T (4)
Wherein, TkIndicate the no-load voltage ratio of k-th of transformer, Tk,maxAnd Tk,minAbove and below the no-load voltage ratio for indicating k-th of transformer respectively
Limit, STIndication transformer set；
(3) capacitive reactance device equality constraint is expressed as：
(Bh-Bh,min)(Bh,max-BhThe h ∈ S of)=0B (5)
Wherein, BhIndicate the susceptance of h-th of capacitive reactance device, Bh,maxAnd Bh,minAbove and below the susceptance for indicating h-th of capacitive reactance device respectively
Limit, SBIndicate capacitive reactance device set.
In the inequality constraints, have：
(1) node voltage inequality constraints is expressed as：
Ui,min+Si≤Ui≤Ui,max-Si i∈N (7)
Wherein, SiIndicate the voltage slack of node i, UiIndicate the voltage magnitude of node i, Ui,maxAnd Ui,minIt indicates respectively
The voltage magnitude bound of node i, N indicate node set；
(2) inequality constraints of voltage slack is expressed as：
Si≥0 i∈N (8)
Wherein, SiIndicate that the voltage slack of node i, N indicate node set；
(3) inequality constraints of generator reactive output is expressed as：
QGi,min≤QG,i≤QGi,max i∈N (9)
Wherein, QG,iIndicate that generator injects reactive power, Q at node i respectivelyGi,maxAnd QGi,minIt indicates at node i respectively
Generator injects reactive power bound, and N indicates node set.
The step 3 specifically includes following steps：
Step 3-1：Idle work optimization model is divided into discrete device equality constraint and continuous quantity Optimized model first；
Step 3-2：Continuous quantity Optimized model is expressed as：
Wherein, f (x) indicates that idle work optimization object function, M indicate set of fingers, UiAnd UjNode i and node j are indicated respectively
Voltage magnitude, GijIndicate the conductance of branch ij, δijIndicate that the phase difference of voltage of node i and node j, N indicate node set, m
Indicate that voltage stability margin is more than the branch of threshold value；SiIndicate the voltage slack of node i,Indicate that the voltage of branch ij is steady
Determine nargin, ω indicates that variation weight, ψ indicate voltage stability index weight, Bh,iIndicate the electricity of h-th of capacitive reactance device at node i
It receives, Tk,iIndicate the no-load voltage ratio of k-th of transformer at node i, Δ Pi、ΔQiActive deviation and idle deviation at node i are indicated respectively；
Definition vector x=[QG,1,…,QG,i,…,QG,N,S1,…,Si,…,SN,U1,…,Ui,…,UN,θ1,…,θi,…,
θN]T, unconfined majorized function L is formed for continuous quantity Optimized model, is had：
Wherein, SGIndicate that generator node, μ indicate disturbance variable, and μ >=0；li、uiUnder the voltage for indicating node i respectively
Slack variable in slack variable and voltage, ypi、yqi、zi、wiIndicate the lagrange's variable of node i, definition vector l=[l1,
l2,……,lr], u=[u1,u2,……,ur], vectorial y={ yp1,…ypi,…,ypN,yq1,…,yqi,…,yqN, vectorial z=
[z1,z2,……,zr], vectorial w=[w1,w2,……,wr], r expression inequality numbers；
Then meet：
Wherein, LxIndicate majorized function L to the local derviation of vector x, LyIndicate majorized function L to the local derviation of vectorial y, LzIt indicates
Majorized function L is to the local derviation of vectorial z, LwIndicate majorized function L to the local derviation of vectorial w, LlIndicate majorized function L to the inclined of vectorial l
It leads, LuIndicate majorized function L to the local derviation of vectorial u, ▽xIndicate idle work optimization object function f (x) to the local derviation of vector x, trend
Equality constraint vector h (x)={ Δ P1,…,ΔPi,…,ΔPN, Δ Q1,…,ΔQi,…,ΔQN, g (x) indicates inequality about
Shu Xiangliang,Indicate the inequality constraints upper limit,gIndicate inequality constraints lower limit, L, Z, U, W are respectively li、zi、ui、wiThe r of composition
Diagonal matrix is tieed up, the r that E is ties up unit matrix；
Step 3-3：The discrete device equality constraint includes transformer voltage ratio equality constraint and capacitive reactance device equality constraint；
Increase duality gapThen have：
Wherein, TkIndicate the no-load voltage ratio of k-th of transformer, Tk,maxAnd Tk,minAbove and below the no-load voltage ratio for indicating k-th of transformer respectively
Limit, STIndication transformer set；BhIndicate the susceptance of h-th of capacitive reactance device, Bh,maxAnd Bh,minThe electricity of h-th of capacitive reactance device is indicated respectively
Receive bound, SBIndicate capacitive reactance device set；
Step 3-4：Idle work optimization model is solved, is specifically included：
Step 1)：Initialization vector l, z, u, w, y, each vector kind of element is respectively 0.5, -0.2,0.5,0.2, -0.01,
It is 0.99 to initialize centripetal parameter Sigma；
Step 2)：Calculate Cgap、Lx、Ly、Lz、Lw、Ll、LuIf calculated Cgap、Lx、Ly、Lz、Lw、Ll、LuRespectively less than count
Precision ε is calculated, then stops calculating；
Step 3)：According to obtained CgapCalculation perturbation variable
Step 4)：Material calculation coefficient stepPAnd stepD, have：
Step 5)：According to stepPAnd stepDX, l, u, y, z, w are corrected, is had：
Wherein, Δ x, Δ l, Δ u, Δ y, Δ z, Δ w indicate respectively vector x, l, u, y, z, w iteration step length；
Step 6)：T is corrected by solving equation (13)kAnd Bh, have：
Tk=Tk+ΔTk (18)
Bh=Bh+ΔBh (19)
Wherein, Δ TkWith Δ BhTo correct step-length；
Step 7)：By revised Tk、Bh, after x, l, u, y, z, w bring equation (12) and (13) into, return to step 2).
The present invention, to the control feature of voltage, reactive apparatus, this class is improved using relaxationization method according to electric system
The convergence of technology.Secondly as there are many some regional discrete devices, and it is main adjustment equipment, the present invention is set discrete
Standby Controlling model introduces idle work optimization model by the way of Constraints, although increasing the convergence time of single optimization analysis
Number, but traditional two suboptimization calculating regular nearby can be converted to primary calculating.Again, electric for existing in power grid
Stable problem is pressed, and since the generation of Voltage-stabilizing Problems usually has locality, utilizes analysis branch voltage steady stability
Index is improved by targetedly reducing the branch voltage stability index numerical value in certain region in Optimized model object function
The voltage stability margin in this region.
Finally it should be noted that：The above embodiments are merely illustrative of the technical scheme of the present invention and are not intended to be limiting thereof, institute
The those of ordinary skill in category field with reference to above-described embodiment still can to the present invention specific implementation mode modify or
Equivalent replacement, these are applying for this pending hair without departing from any modification of spirit and scope of the invention or equivalent replacement
Within bright claims.
Claims (4)
1. a kind of idle work optimization method adapting to bulk power grid automatism voltage control, it is characterised in that：The method includes following steps
Suddenly：
Step 1：Determine the voltage stability index of branch；
Step 2：Establish idle work optimization model；
Step 3：Idle work optimization model is solved using interior point method；
In the step 2, idle work optimization model includes idle work optimization object function and idle work optimization constraints；
The idle work optimization constraints includes equality constraint and inequality constraints；
The equality constraint includes trend equality constraint, transformer voltage ratio equality constraint and capacitive reactance device equality constraint；
The inequality constraints includes that node voltage inequality constraints, the inequality constraints of voltage slack and generator reactive are contributed
Inequality constraints；
The idle work optimization object function is expressed as：
Wherein, f (x) indicates that idle work optimization object function, M indicate set of fingers, UiAnd UjThe electricity of node i and node j is indicated respectively
Pressure amplitude value, GijIndicate the conductance of branch ij, δijIndicate that the phase difference of voltage of node i and node j, N indicate that node set, m indicate
Voltage stability margin is more than the branch of threshold value；SiIndicate the voltage slack of node i,Indicate that the voltage stabilization of branch ij is abundant
Degree, ω indicate that variation weight, ψ indicate voltage stability index weight；
In the equality constraint, have：
(1) trend equality constraint is expressed as：
Wherein, PG,iAnd QG,iGenerator injection active power and reactive power, P at node i are indicated respectivelyD,iAnd QD,iIt indicates respectively
Generated power load and load or burden without work, U at node iiAnd UjThe voltage magnitude of node i and node j, G are indicated respectivelyijAnd BijPoint
Not Biao Shi branch ij conductance and susceptance, N indicate node set, δijFor the phase difference of voltage of node i and node j, and δij=θi-
θj, wherein θi、θjRespectively node i, j voltage phase angle；
(2) transformer voltage ratio equality constraint is expressed as：
(Tk-Tk,min)(Tk,max-TkThe k ∈ S of)=0T (4)
Wherein, TkIndicate the no-load voltage ratio of k-th of transformer, Tk,maxAnd Tk,minThe no-load voltage ratio bound of k-th of transformer, S are indicated respectivelyT
Indication transformer set；
(3) capacitive reactance device equality constraint is expressed as：
(Bh-Bh,min)(Bh,max-BhThe h ∈ S of)=0B (5)
Wherein, BhIndicate the susceptance of h-th of capacitive reactance device, Bh,maxAnd Bh,minThe susceptance bound of h-th of capacitive reactance device, S are indicated respectivelyB
Indicate capacitive reactance device set；
The step 3 specifically includes following steps：
Step 3-1：Idle work optimization model is divided into discrete device equality constraint and continuous quantity Optimized model first；
Step 3-2：Continuous quantity Optimized model is expressed as：
Wherein, f (x) indicates that idle work optimization object function, M indicate set of fingers, UiAnd UjThe electricity of node i and node j is indicated respectively
Pressure amplitude value, GijIndicate the conductance of branch ij, δijIndicate that the phase difference of voltage of node i and node j, N indicate that node set, m indicate
Voltage stability margin is more than the branch of threshold value；SiIndicate the voltage slack of node i,Indicate that the voltage stabilization of branch ij is abundant
Degree, ω indicate that variation weight, ψ indicate voltage stability index weight, Bh,iIndicate the susceptance of h-th of capacitive reactance device at node i,
Tk,iIndicate the no-load voltage ratio of k-th of transformer at node i, Δ Pi、ΔQiActive deviation and idle deviation at node i are indicated respectively；
Definition vector x=[QG,1,…,QG,i,…,QG,N,S1,…,Si,…,SN,U1,…,Ui,…,UN,θ1,…,θi,…,θN]T,
Unconfined majorized function L is formed for continuous quantity Optimized model, is had：
Wherein, SGIndicate that generator node, μ indicate disturbance variable, and μ >=0；li、uiRelaxation under the voltage of node i is indicated respectively
Slack variable on variable and voltage, ypi、yqi、zi、wiIndicate the lagrange's variable of node i, definition vector l=[l1,
l2,……,lr], u=[u1,u2,……,ur], vectorial y={ yp1,…ypi,…,ypN,yq1,…,yqi,…,yqN, vectorial z=
[z1,z2,……,zr], vectorial w=[w1,w2,……,wr], r expression inequality numbers；
Then meet：
Wherein, LxIndicate majorized function L to the local derviation of vector x, LyIndicate majorized function L to the local derviation of vectorial y, LzIndicate optimization
Function L is to the local derviation of vectorial z, LwIndicate majorized function L to the local derviation of vectorial w, LlIndicate local derviations of the majorized function L to vectorial l,
LuIndicate majorized function L to the local derviation of vectorial u, ▽xIndicate idle work optimization object function f (x) to the local derviation of vector x, trend equation
Constrained vector h (x)={ Δ P1,…,ΔPi,…,ΔPN, Δ Q1,…,ΔQi,…,ΔQN, g (x) indicate inequality constraints to
Amount,Indicate the inequality constraints upper limit,gIndicate inequality constraints lower limit, L, Z, U, W are respectively li、zi、ui、wiThe r dimensions pair of composition
Angular moment battle array, the r that E is tie up unit matrix；
Step 3-3：The discrete device equality constraint includes transformer voltage ratio equality constraint and capacitive reactance device equality constraint；
Increase duality gapThen have：
Wherein, TkIndicate the no-load voltage ratio of k-th of transformer, Tk,maxAnd Tk,minThe no-load voltage ratio bound of k-th of transformer, S are indicated respectivelyT
Indication transformer set；BhIndicate the susceptance of h-th of capacitive reactance device, Bh,maxAnd Bh,minOn the susceptance for indicating h-th of capacitive reactance device respectively
Lower limit, SBIndicate capacitive reactance device set；
Step 3-4：Idle work optimization model is solved, is specifically included：
Step 1)：Initialization vector l, z, u, w, y, each vector kind of element is respectively 0.5, -0.2,0.5,0.2, -0.01, initially
It is 0.99 to change centripetal parameter Sigma；
Step 2)：Calculate Cgap、Lx、Ly、Lz、Lw、Ll、LuIf calculated Cgap、Lx、Ly、Lz、Lw、Ll、LuRespectively less than calculate essence
ε is spent, then stops calculating；
Step 3)：According to obtained CgapCalculation perturbation variable
Step 4)：Material calculation coefficient stepPAnd stepD, have：
Step 5)：According to stepPAnd stepDX, l, u, y, z, w are corrected, is had：
Wherein, Δ x, Δ l, Δ u, Δ y, Δ z, Δ w indicate respectively vector x, l, u, y, z, w iteration step length；
Step 6)：T is corrected by solving equation (13)kAnd Bh, have：
Tk=Tk+ΔTk (18)
Bh=Bh+ΔBh (19)
Wherein, Δ TkWith Δ BhTo correct step-length；
Step 7)：By revised Tk、Bh, after x, l, u, y, z, w bring equation (12) and (13) into, return to step 2).
2. the idle work optimization method according to claim 1 for adapting to bulk power grid automatism voltage control, it is characterised in that：It is described
Step 1 specifically includes following steps：
Step 1-1：The voltage stability margin L of branch ij between calculate node i and node jij, steady as the voltage of branch ij
Determine index；
Step 1-2：Judge LijWhether it is more than threshold value T, if more than respective branch is contributed to set of fingers M.
3. the idle work optimization method according to claim 2 for adapting to bulk power grid automatism voltage control, it is characterised in that：It is described
In step 1-1, the voltage stability margin of branch ij is L between node i and node jijIt is expressed as：
Wherein, R and X indicates resistance and the reactance of branch ij, δijIndicate the phase difference of voltage of node i and node j, UiIndicate node
The voltage magnitude of i, PjAnd QjThe outflow active power and reactive power of node j are indicated respectively.
4. the idle work optimization method according to claim 1 for adapting to bulk power grid automatism voltage control, it is characterised in that：It is described
In inequality constraints, have：
(1) node voltage inequality constraints is expressed as：
Ui,min+Si≤Ui≤Ui,max-Si i∈N (7)
Wherein, SiIndicate the voltage slack of node i, UiIndicate the voltage magnitude of node i, Ui,maxAnd Ui,minNode is indicated respectively
The voltage magnitude bound of i, N indicate node set；
(2) inequality constraints of voltage slack is expressed as：
Si≥0 i∈N (8)
Wherein, SiIndicate that the voltage slack of node i, N indicate node set；
(3) inequality constraints of generator reactive output is expressed as：
QGi,min≤QG,i≤QGi,max i∈N (9)
Wherein, QG,iIndicate that generator injects reactive power, Q at node i respectivelyGi,maxAnd QGi,minIt indicates to generate electricity at node i respectively
Machine injects reactive power bound, and N indicates node set.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510341645.6A CN106329537B (en) | 2015-06-18 | 2015-06-18 | A kind of idle work optimization method adapting to bulk power grid automatism voltage control |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510341645.6A CN106329537B (en) | 2015-06-18 | 2015-06-18 | A kind of idle work optimization method adapting to bulk power grid automatism voltage control |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106329537A CN106329537A (en) | 2017-01-11 |
| CN106329537B true CN106329537B (en) | 2018-10-23 |
Family
ID=57732337
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510341645.6A Active CN106329537B (en) | 2015-06-18 | 2015-06-18 | A kind of idle work optimization method adapting to bulk power grid automatism voltage control |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106329537B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107482639B (en) * | 2017-07-19 | 2019-12-13 | 清华大学 | dynamic reactive power equipment reactive power optimization method suitable for wind power system in power grid |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103199542A (en) * | 2013-02-26 | 2013-07-10 | 中国电力科学研究院 | Method of optimal control of wind power plant reactive voltage |
| CN103795068A (en) * | 2014-03-05 | 2014-05-14 | 广东电网公司电力调度控制中心 | Optimal configuration method for high-voltage distribution network dynamic reactive power compensation equipment capacity |
| CN104300546A (en) * | 2014-10-13 | 2015-01-21 | 国家电网公司 | Voltage stability constraint reactive power optimization method based on wide-area measurement information |
-
2015
- 2015-06-18 CN CN201510341645.6A patent/CN106329537B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103199542A (en) * | 2013-02-26 | 2013-07-10 | 中国电力科学研究院 | Method of optimal control of wind power plant reactive voltage |
| CN103795068A (en) * | 2014-03-05 | 2014-05-14 | 广东电网公司电力调度控制中心 | Optimal configuration method for high-voltage distribution network dynamic reactive power compensation equipment capacity |
| CN104300546A (en) * | 2014-10-13 | 2015-01-21 | 国家电网公司 | Voltage stability constraint reactive power optimization method based on wide-area measurement information |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106329537A (en) | 2017-01-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103715700B (en) | Be applicable to powerless control system and the control method of wind farm grid-connected point voltage control | |
| CN102361324B (en) | Method for regulating terminal voltage by reactive power support of double-fed wind generator unit and system thereof | |
| CN104578084B (en) | Dynamic reactive compensating mechanism and AVC (Automatic Voltage Control) combined control system | |
| CN104300552B (en) | A kind of take into account correction voltage out-of-limit with reduce active loss idle work optimization method | |
| CN107732917B (en) | A kind of closed loop network turn power supply Load flow calculation optimization method | |
| CN104158199B (en) | Power system real-time status is carried out the system and method for reactive power and voltage control | |
| CN105226665A (en) | The method of busbar voltage daily fluctuation rate is reduced by automatic plant stand cooperation control | |
| CN107546763A (en) | Photovoltaic generation is maximum in the power distribution network of a kind of meter and different pressure regulation strategies receives capacity calculation method | |
| CN106712044A (en) | Multi-objective optimization method of reactive voltage in wind power plant | |
| CN102769299B (en) | Wind farm voltage control method based on voltage operating state | |
| CN106773644B (en) | A kind of AGC control systems and its method changed based on the heat supply amount of drawing gas | |
| CN108321848A (en) | A kind of distributed electrical operated control method being adapted with OLTC voltage regulating modes | |
| CN109617112A (en) | Modified DC voltage control strategy suitable for multiterminal flexible direct current system | |
| CN107465198B (en) | A kind of wind power field automatic voltage control and system | |
| CN108365610A (en) | A kind of reactive voltage control method coordinated based on source net | |
| CN106253305A (en) | A kind of full-separate isolated island powerless control method for distributed power source | |
| CN106329537B (en) | A kind of idle work optimization method adapting to bulk power grid automatism voltage control | |
| CN105071397B (en) | The coordination reactive voltage control method for the polymorphic type reactive-load compensation equipment that wind-powered electricity generation is sent outside | |
| CN106159957B (en) | A kind of optimization method exchanging weak section steady-state overvoltage Pre-control measures | |
| CN106169760A (en) | Main website AVC system and substation SVG system coordination control method | |
| CN107069757B (en) | A kind of mixing multi-infeed HVDC system secondary voltage control method | |
| CN108173277A (en) | A kind of change Target Control Method of flexible direct current power transmission system | |
| CN102904265A (en) | Transformer station voltage reactive control method and system based on region tidal current | |
| CN104300546A (en) | Voltage stability constraint reactive power optimization method based on wide-area measurement information | |
| CN104917170A (en) | Method for adjusting voltage and frequency through micro power grid self-adaption droop control based on PI control |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 100192 Beijing city Haidian District Qinghe small Camp Road No. 15 Applicant after: CHINA ELECTRIC POWER RESEARCH INSTITUTE Co.,Ltd. Applicant after: State Grid Corporation of China Address before: 100192 Beijing city Haidian District Qinghe small Camp Road No. 15 Applicant before: China Electric Power Research Institute Applicant before: State Grid Corporation of China |
|
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20180329 Address after: 100192 Beijing city Haidian District Qinghe small Camp Road No. 15 Applicant after: CHINA ELECTRIC POWER RESEARCH INSTITUTE Co.,Ltd. Applicant after: State Grid Corporation of China Applicant after: STATE GRID ZHEJIANG ELECTRIC POWER Co.,Ltd. Address before: 100192 Beijing city Haidian District Qinghe small Camp Road No. 15 Applicant before: CHINA ELECTRIC POWER RESEARCH INSTITUTE Co.,Ltd. Applicant before: State Grid Corporation of China |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |