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 PDF

Info

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
Application number
CN201510341645.6A
Other languages
Chinese (zh)
Other versions
CN106329537A (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.)
State Grid Corp of China SGCC
State Grid Zhejiang Electric Power Co Ltd
China Electric Power Research Institute Co Ltd CEPRI
Original Assignee
State Grid Corp of China SGCC
State Grid Zhejiang Electric Power Co Ltd
China Electric Power Research Institute Co Ltd CEPRI
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 State Grid Corp of China SGCC, State Grid Zhejiang Electric Power Co Ltd, China Electric Power Research Institute Co Ltd CEPRI filed Critical State Grid Corp of China SGCC
Priority to CN201510341645.6A priority Critical patent/CN106329537B/en
Publication of CN106329537A publication Critical patent/CN106329537A/en
Application granted granted Critical
Publication of CN106329537B publication Critical patent/CN106329537B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/30Reactive 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

A kind of idle work optimization method adapting to bulk power grid automatism voltage control
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 δijij, 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,…,UN1,…,θ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 δijij, 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,…,UN1,…,θ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 δiji- θ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,…,UN1,…,θ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.
CN201510341645.6A 2015-06-18 2015-06-18 A kind of idle work optimization method adapting to bulk power grid automatism voltage control Active CN106329537B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (3)

* Cited by examiner, † Cited by third party
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
CN104300552B (en) A kind of take into account correction voltage out-of-limit with reduce active loss idle work optimization method
CN105244901B (en) A kind of nonlinear decentralized control method of HVDC transmission system
CN105226665A (en) The method of busbar voltage daily fluctuation rate is reduced by automatic plant stand cooperation control
CN107732917B (en) A kind of closed loop network turn power supply Load flow calculation optimization method
CN107465198B (en) A kind of wind power field automatic voltage control and system
CN106329537B (en) A kind of idle work optimization method adapting to bulk power grid automatism voltage control
CN104158199B (en) Power system real-time status is carried out the system and method for reactive power and voltage control
CN106773644B (en) A kind of AGC control systems and its method changed based on the heat supply amount of drawing gas
CN104300546A (en) Voltage stability constraint reactive power optimization method based on wide-area measurement information
CN105262112B (en) Wind power plant cluster formula static var compensator control method
CN106159957B (en) A kind of optimization method exchanging weak section steady-state overvoltage Pre-control measures
CN108321848A (en) A kind of distributed electrical operated control method being adapted with OLTC voltage regulating modes
CN107069757B (en) A kind of mixing multi-infeed HVDC system secondary voltage control method
Sharma et al. Optimal location for shunt connected FACTS devices in a series compensated long transmission line
CN108173277A (en) A kind of change Target Control Method of flexible direct current power transmission system
CN105071397A (en) Coordinated reactive voltage control method of different reactive compensation devices of wind power delivery
CN104917170A (en) Method for adjusting voltage and frequency through micro power grid self-adaption droop control based on PI control
CN107465213A (en) A kind of power station AVC systems and its idle regulated quantity computational methods in real time
CN102904265B (en) Transformer station voltage reactive control method and system based on region tidal current
CN106712044A (en) Multi-objective optimization method of reactive voltage in wind power plant
CN107171354A (en) The sagging coefficient calculation method in flexible direct current line commutation station
CN106169760A (en) Main website AVC system and substation SVG system coordination control method
CN106253305A (en) A kind of full-separate isolated island powerless control method for distributed power source

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
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 Science Research Institute Co., Ltd.

Applicant after: State Grid Corporation of China

Applicant after: Zhejiang Electric Power Co., Ltd.

Address before: 100192 Beijing city Haidian District Qinghe small Camp Road No. 15

Applicant before: China Electric Science Research Institute Co., Ltd.

Applicant before: State Grid Corporation of China

TA01 Transfer of patent application right
CB02 Change of applicant information

Address after: 100192 Beijing city Haidian District Qinghe small Camp Road No. 15

Applicant after: China Electric Science 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

GR01 Patent grant
GR01 Patent grant