CN103001218B - Reactive voltage coordination control method based on comprehensive benefits of transmission and distribution power grid - Google Patents

Abstract

The invention provides a reactive voltage coordination control method based on comprehensive benefits of a transmission and distribution power grid. The reactive voltage coordination control method comprises enabling a subordinate power grid to upload gate characteristics to a superior control center; enabling the superior control center to establish a superior and subordinate reactive voltage coordinated optimization global model according to the gate characteristics; performing virtual control partitioning on a superior administered power grid, and obtaining a coordinated optimization calculation model of a plurality of virtual control partitions; enabling the superior control center to obtain optimized values of gate coordination variables of the subordinate power grid by calculation according to the coordinated optimization calculation model and send the optimized values to the subordinate power grid; enabling the subordinate power grid to add coordination variable restraint in an automatic voltage control (AVC) algorithm according to the optimized values of the gate coordination variables, and obtaining a virtual control partition optimization strategy based on the comprehensive benefits of the transmission and distribution power grid; and enabling the superior power grid and the subordinate power grid to perform reactive voltage coordination control according to the virtual control partition optimization strategy. The method improves the running state of the whole network, improves system benefits, reduces model calculation scales and complexity, and is good in reliability and high in compatibility.

Description

Based on the reactive voltage control method for coordinating of transmission and distribution network comprehensive benefit

Technical field

The present invention relates to learning algorithms technical field, particularly a kind of reactive voltage control method for coordinating based on transmission and distribution network comprehensive benefit.

Background technology

In power system reactive power voltage control field, dispatch system different in transmission and distribution network determines the executor of interregional cooperation control.In the dispatch system of layering and zoning, cooperation control is performed by superior control centre, and subordinate's electrical network is by receiving the critical point set point that issues of superior control centre or interval, uneoupled control between feasible region.Superior control centre, not only as Consultation Center, also can have the controllable resources of electric power networks and the correspondence of oneself administering simultaneously, and such as, in China's electrical network system, economizing tune is exactly the higher level that each lower possession is adjusted, the equipment of simultaneously self also administering to some extent.In order to ensure network-wide security, high-quality, economical operation, superior control centre, when being optimized administered regional power grid and controlling, also needs the regulating power, the demand for control and benefit etc. that take into full account subordinate's electrical network.

Such as, economize ground AVC (Automatic Voltage Control, automatism voltage control) coordinate to belong to the skimble-scamble situation of control centre of the superior and the subordinate information modeling, control centre is adjusted to introduce subordinate's electrical network full detail if economize, then real time data amount is too huge, therefore general is equivalent to load bus subordinate's electrical network, does not carry out complete modeling to subordinate's electrical network.If only have the imperfect information of subordinate's electrical network will be difficult to from the whole network angle calculation optimal load flow, thus each critical point specification of variables value cannot be obtained.Therefore traditional way only makes province's ground two-stage operation of power networks in reasonable level.

In order to ensure network-wide security, high-quality, economical operation, power transmission network side control centre is while being optimized institute's compass of competency electrical network and controlling, need to consider and the ability of its information modeling skimble-scamble subordinate electrical network, demand and benefit, and coordinate, and prior art can not well solve the problem.

Summary of the invention

The present invention is intended at least to solve one of technical problem existed in prior art.For this reason, the object of the invention is to propose a kind of reactive voltage control method for coordinating based on transmission and distribution network comprehensive benefit.This method improves the whole network running status, improves system benefit, reduces model and calculates scale and complexity, and good reliability, compatible high.

For achieving the above object, embodiments of the invention propose a kind of reactive voltage control method for coordinating based on transmission and distribution network comprehensive benefit, comprise the following steps: electrical network superior control centre of subordinate uploads critical point characteristic, wherein, described critical point characteristic comprises regulating power, demand for control and optimization benefit; Described critical point characteristic is converted into the coordination variable constraint of critical point variable sign and target correction term by described superior control centre, sets up reactive voltage coordination optimization world model of the superior and the subordinate; According to reactive voltage coordination optimization world model of described the superior and the subordinate, electrical network is administered to described higher level and carries out virtual controlling subregion, obtain the coordination optimization computation model of multiple virtual controlling subregion; Described superior control centre, according to described coordination optimization computation model, calculates the optimal value of the critical point coordination variable of described subordinate electrical network and is sent to described subordinate electrical network; Described higher level's electrical network and described subordinate electrical network carry out reactive voltage cooperation control according to described virtual controlling partition zone optimizing strategy.

The critical point characteristic of subordinate's electrical network is considered according to the reactive voltage control method for coordinating based on transmission and distribution network comprehensive benefit of the embodiment of the present invention, control partition is carried out to higher level's power transmission network, and computation model and the Practical resolving Algorithm of control partition are proposed, improve the running status of the whole network, improve the system benefit of the whole network.By by bulk power grid reactive voltage Harmonic Control decomposition and inversion be dimensionality reduction control partition in the superior and the subordinate's electric network reactive-load voltage coordination problem, reduce computation model scale and computation complexity, and improve the compatibility of traditional algorithm to the benefit expression-form that subordinate uploads, reduce the calculation cost that subordinate's electrical network is extra.

In one embodiment of the invention, described regulating power is for representing that described subordinate electrical network is to the idle controllable scope in described critical point.

In one embodiment of the invention, described demand for control is qualified for controlling described subordinate line voltage, and the restriction range that shown critical point voltage needs higher level to control.

In one embodiment of the invention, described optimization benefit is when the described virtual controlling partition zone optimizing strategy of control does not clash in benefit, each subordinate grid benefit considered in coordination optimization computation model.

In one embodiment of the invention, described optimization benefit, according to subordinate electrical network AVC, with described critical point variable for independent variable, is stated by curve form or range format.Use different method for expressing as the case may be, what consider that subordinate's network optimization benefit uploads is multi-form, improves and uploads to subordinate the compatibility optimizing benefit expression-form

In one embodiment of the invention, also comprise innovatory algorithm, for virtual controlling subregion is optimized calculating, integrates and generate new virtual partition, wherein, the number of described new virtual partition is less than the virtual partition number before optimization.By bulk power grid reactive voltage Harmonic Control decomposition and inversion be dimensionality reduction control partition in the superior and the subordinate's electric network reactive-load voltage coordination problem, and further virtual partition is optimized, reduce computation model scale and computation complexity, subregion uneoupled control improves the reliability of AVC algorithm simultaneously.If certain subregion occurs abnormal, the normal coordination optimization that generally also can not have influence on other subregions controls.

In one embodiment of the invention, described superior control centre is optimized calculating, obtains each control partition VACE set point.

In one embodiment of the invention, described virtual controlling subregion is optimized calculating, integrates and generate new virtual partition and comprise further:

Step 1: select constraint subspace S _in' start to calculate;

Step 2: the described constraint subspace making current calculating is S _ij', solve following computation model:

min f _H(x',u',Q _Li,U _Hi)

s.t.g _i(x',u',U _Hi,Q _Li)=0

h _j(x',u',U _Hi,Q _Li)≤0

x' _min<x'<x' _max

u' _min<u'<u' _max

Q _Lij∈S' _ij=[l' _ij,u' _ij]i=1,2,…,N,1≤j≤n

U _Hi，min<U _Hi<U _Hi，maxi＝1,2,…,N'

-δ≤I _i,VACE+C _i,VACEΔu' _i≤δ

Wherein, x ', u ' be the subset of higher level's electric network state amount x and controlled quentity controlled variable u, U is described critical point voltage, and Q is that described critical point is idle, min and max represents corresponding minimum and maximum threshold value respectively, I _{i, VACE}for the VACE index of described constraint subspace i, C _{i, VACE}for described constraint subspace i is to the control sensitivity of VACE, Δ u _i' be control deviation, δ is the VACE maximum deviation that region voltage control system allows.Calculate locally optimal solution f _{ij, opt}', can obtain coordinating subarea i and retrain subspace S _in' locally optimal solution drop point x _ij' and locally optimal solution f _{hL, ij, opt}, wherein,

f _HL,ij,opt=f' _ij,opt+f' _ij

Step 3: if drop point x _ijbelong to other subspaces S _ik', then basis relation, obtain subspace S _ik' locally optimal solution f _{hL, ik, opt}, wherein, f _{hL, ik, opt}=f' _{ij, opt}+ f ' _ik.

If described drop point x _ij' belong to multiple subspace, then interval rank is numbered minimum interval numbering and be designated as N _tEMP.Due to f ' _ik<f' _ij, described interval rank numbering N _tEMPlocally optimal solution f _{hL, ik, opt}for the solution that the network loss in described multiple space is relatively minimum.

Step 4: if first time calculates, remember described interval rank numbering N _tEMPdescribed locally optimal solution f _{hL, ik, opt}for f _{hL, i, temp}, and remember current drop point x _ij' be x _{i, temp}.Once calculated if carried out, more described f _{hL, ik, opt}with described f _{hL, i, temp}if meet f _{hL, ik, opt}<f _{hL, i, temp}, then f is made _{hL, i, temp}=f _{hL, ik, opt}, x _{i, temp}=x _ij'; Otherwise, if meet f _{hL, i, temp}<f' _{ij, opt}+ f _i1', then innovatory algorithm iteration terminates.

Step 5: if described interval numbering N _tEMP+ 1 is non-vanishing, gets corresponding described constraint subspace perform step 2.

If described interval numbering N _tEMP+ 1 is zero, then innovatory algorithm iteration terminates.

Step 6: single described control partition i calculates and terminates, and obtains the described optimal solution f of described control partition i _{hL, i, opt}with optimal solution drop point x _i, wherein, f _{hL, i, opt}=f _{hL, i, temp}, x _i=x _{i, temp}.

According to the embodiment of the present invention.。。。。

Additional aspect of the present invention and advantage will part provide in the following description, and part will become obvious from the following description, or be recognized by practice of the present invention.

Accompanying drawing explanation

Above-mentioned and/or additional aspect of the present invention and advantage will become obvious and easy understand from accompanying drawing below combining to the description of embodiment, wherein:

Fig. 1 is the flow chart of the reactive voltage control method for coordinating based on transmission and distribution network comprehensive benefit according to the embodiment of the present invention;

Fig. 2 is that the corresponding control partition of IEEE39 node system adopted in simulation example divides schematic diagram

Total active power loss contrast schematic diagram of the system of algorithms of different that Fig. 3 is simulation example total system when being node load or burden without work sustainable growth situation;

The power transmission network active power loss contrast schematic diagram of algorithms of different when Fig. 4 is node load or burden without work sustainable growth situation;

The power distribution network of algorithms of different total active power loss contrast schematic diagram when Fig. 5 is node load or burden without work sustainable growth situation;

The subsystem A power transmission network side active power loss contrast schematic diagram of algorithms of different when Fig. 6 is node load or burden without work sustainable growth situation;

The subsystem B power transmission network side active power loss contrast schematic diagram of algorithms of different during Fig. 7 node load or burden without work sustainable growth situation;

No. 1 power distribution network active power loss contrast schematic diagram of algorithms of different when Fig. 8 is node load or burden without work sustainable growth situation;

No. 2 power distribution network active power loss contrast schematic diagrames of algorithms of different when Fig. 9 is node load or burden without work sustainable growth situation;

No. 1 power distribution network critical point voltage-contrast schematic diagram of algorithms of different when Figure 10 is node load or burden without work sustainable growth situation;

No. 2 power distribution network critical point voltage-contrast schematic diagrames of algorithms of different when Figure 11 is node load or burden without work sustainable growth situation;

The idle contrast schematic diagram in No. 1 power distribution network critical point of algorithms of different when Figure 12 is node load or burden without work sustainable growth situation; With

The idle contrast schematic diagram in No. 2 power distribution network critical points of algorithms of different when Figure 13 is node load or burden without work sustainable growth situation.

Embodiment

Be described below in detail embodiments of the invention, the example of described embodiment is shown in the drawings, and wherein same or similar label represents same or similar element or has element that is identical or similar functions from start to finish.Being exemplary below by the embodiment be described with reference to the drawings, only for explaining the present invention, and can not limitation of the present invention being interpreted as.

Describe the reactive voltage control method for coordinating based on transmission and distribution network comprehensive benefit according to the embodiment of the present invention below with reference to Fig. 1, comprise the following steps:

Step S110: electrical network superior control centre of subordinate uploads critical point characteristic, wherein, critical point characteristic comprises regulating power, demand for control and optimization benefit.

This method introduces three critical point characteristics of subordinate's electrical network, comprising: regulating power, demand for control and optimization benefit.

(1), regulating power.Regulating power is for representing that subordinate's electrical network is to the idle controllable scope in critical point.

Critical point is idle Q _limain by subordinate's electric network influencing, therefore need subordinate's electrical network to upload each critical point regulating power.Critical point is idle Q _liin, power distribution network is system loading place, and because inner general not containing generator, idle control device is mainly capacity reactance device group, the first-class voltage adjusting device of on-load voltage regulation tap, and therefore power distribution network is mainly had an impact to power distribution network critical point is idle by the switching of capacity reactance device group.In addition, higher level's power transmission network also can affect subordinate's power distribution network and inject idle, mainly contains following three aspects:

1, change in voltage causes capacitive reactance device state to change, and then it is idle to affect the injection of power distribution network critical point.

2, the acting in conjunction of change in voltage and static load characteristics causes the idle change in distribution critical point.

3, change in voltage causes the reactive loss of distribution network to change.

When change in voltage is comparatively large, during as voltage collapse, 1,2 liang of aspect just has considerable influence to critical point is idle.In reality under normal circumstances, power transmission network change in voltage amplitude is less, and therefore power transmission network belongs to the adjustment that critical point is idle and regulates among a small circle, and critical point is idle generally to be regulated by power distribution network.

(2), demand for control.Demand for control is qualified for controlling subordinate's line voltage, and the restriction range that shown critical point voltage needs higher level to control.Qualified in order to ensure subordinate's line voltage, the main critical point voltage U by power transmission network regulation and control _hneed to control in certain restriction range.Therefore, subordinate's electrical network is needed to upload U according to current loads level _hrange of needs.In the voltage of critical point, power transmission network adjusts main transformer impedance much smaller in the equiva lent impedance of power distribution network with comparing, so the control of power distribution network relevant device is very little to the effect of critical point voltage, and critical point voltage mainly affects by power transmission network.

(3), benefit is optimized.Optimize benefit for when control virtual controlling partition zone optimizing strategy does not clash in benefit, each subordinate grid benefit considered in coordination optimization computation model.Wherein, virtual controlling partition zone optimizing, can think a kind of method of the superior and the subordinate's reactive voltage coordination optimization.Each subordinate grid net loss summation accounts for the whole network network loss ratio and sometimes cannot ignore.In order to ensure that the superior and the subordinate's reactive power and voltage control does not clash in benefit, need the optimization benefit introducing each subordinate electrical network in coordination optimization computation model.Wherein, coordination optimization computation model and virtual controlling partition zone optimizing recited above, and a virtual controlling subregion, include a part corresponding with this subordinate's electrical network in subordinate's electrical network and higher level's electrical network.

In one embodiment of the invention, optimize benefit according to subordinate electrical network AVC, with critical point variable for independent variable, by curve form or range format statement.

In upper level model, subordinate's electrical network occurs with the form of critical point variable.When electrical network superior control centre of subordinate uploads critical point characteristic, regulating power, demand for control are changed into two critical points idle Q in variable critical point by subordinate's electrical network respectively _liwith critical point voltage U _hiconstraint.The expression-form optimizing benefit can have following three kinds:

(1), perfect information model tormulation.Optimize benefit by subordinate electrical network perfect information model f _li(x) accurate expression.This form is suitable for the unified situation of the superior and the subordinate's electric network information modeling.Such as the situation that in coordination, power transmission network and power transmission network are coordinated economized by net.The application benefit of subordinate's electrical network perfect information model is also comparatively remarkable, but and be not suitable for the situation that the skimble-scamble power transmission network of the superior and the subordinate's information model and power transmission network coordinate because in this case, subordinate's electrical network many non-critical points information is at higher level's electrical network not observable.

(2), parameter of curve is expressed.Based on current power distribution net state, according to real-time virtual mass data, optimization benefit is carried out parameter fitting, the algebraic expression that to be then expressed as with critical point variable be independent variable, as: f _li(U _hi, Q _li).

Wherein, due to critical point voltage U _himain to control by power transmission network, power distribution network is smaller to its effect, and in different coordination optimization strategy, critical point voltage value generally changes not quite.Generally can critical point voltage U _hiregard constant as.Therefore f _li(U _hi, Q _li) can be similar to and think the idle Q in critical point _livariable expression, shown in formula specific as follows:

f _Li(U _Hi,Q _Li)=f _Li(Q _Li|U _Hi≈1.0)=a _inQ _Li ⁿ+a _in-1Q _Li ^n-1+...+a _i1Q _Li+a _i0

Wherein, a _im, a _im-1a _i0for coefficient.The key of this expression way is to need real-time virtual mass data.Virtual data provides different critical point variablees by supposition superior control centre, under current distribution state, repeatedly calls subordinate AVC algorithm and calculate by regulating available optimum benefit.Distribution state can not be enumerated by total number at off-line, and needs the coefficient above real-time fitting in multinomial, adds calculation cost extra in a large number.

(3), the form of piecewise function is expressed.Part do not built or built the power distribution network of AVC, adopts the expression-form of the matched curve described in (two) to realize difficulty comparatively large, therefore, being easier in actual conditions provide and carrying out subsection fuzzy expression benefit with range format.Such as, for the power distribution network being still in the reactive power compensator construction period, need superior control centre to support, idle power grid control can not regard merely continuous control process as, and therefore its benefit is expressed and also should be adopted piecewise function form.For built AVC system but the more weak situation of its ability, subordinate's electrical network is difficult to increase real-time virtual data and carry out the functions such as parameter fitting to current benefit, piecewise function form also can be adopted to state and optimize benefit.Piecewise function form is as follows:

$f_{\mathrm{Li}}\left(Q_{\mathrm{Li}}\right)=\left\{\begin{array}{cc}{f}_{i1}& Q_{\mathrm{Li}}\&Element;S_{i1}=[l_{i1},u_{i1}]\\ f_{i2}& Q_{\mathrm{Li}}\&Element;S_{i2}=[l_{i2},u_{i2}]\\ & ...\\ f_{\mathrm{in}}& Q_{\mathrm{Li}}\&Element;S_{i2}=[l_{\mathrm{in}},u_{\mathrm{in}}]\end{array}\right.,$

Wherein, S ₁∪ S ₂∪ ...=S _i, S _ifor the idle Q in critical point _lifeas ible space.Identical with distribution demand for control, subinterval non-overlapping copies.In above formula, the benefit in each interval gets its average level usually.The given of average level not only can based on real time data, also can based on Expert Rules or operating experience, and applicability is stronger.

Step S120: critical point characteristic is converted into the coordination variable constraint of critical point variable sign and target correction term by superior control centre, sets up reactive voltage coordination optimization world model of the superior and the subordinate.

In the superior and the subordinate's reactive voltage cooperation control, coordination variable carries out mutual key between the superior and the subordinate, containing the information realizing global coordination optimization in cooperation control.The selection of coordination variable, the formulation of coordination variable constraint, coordination variable, on the impact coordinating target, determine the implementation effect of the superior and the subordinate's cooperation control jointly.

The coordination variable of this method characterizes with critical point variable, according to the critical point characteristic received, adopts the idle Q in the on high-tension side critical point of system interconnection transformer conventional at present _liwith the critical point voltage U of power transmission network root node _hrepresent.

Such as, suppose existence 1 superior control centre, under connect N number of net, wherein N ₁individual sheet network optimization benefit adopts parameter of curve to express, N-N ₁individual distribution is optimized benefit and is adopted segmentation letter segment table to reach.In conjunction with critical point recited above characteristic, consider that the reactive voltage coordination optimization world model of the superior and the subordinate of transmission and distribution network comprehensive benefit can be written as,

min f _H(x,u,Q _Li,U _Hi)+∑f _Li(Q _Li)

s.t.g _i(x,u,U _Hi,Q _Li)=0

h _j(x,u,U _Hi,Q _Li)≤0

x _min<x<x _max

u _min<u<u _max

Q _Li，min<Q _Li<Q _Li，maxi＝1,2,…,N

U _Hi，min<U _Hi<U _Hi，maxi＝1,2,…,N

Wherein, x, u are higher level's electric network state amount and controlled quentity controlled variable respectively, equality constraint mainly trend constraint, and inequality constraints is the critical point adjustable range uploaded of the constraint of higher level's reactive-load adjusting device and subordinate and critical point demand for control mainly, the cumulative item ∑ f in target function _li(Q _li) be the optimization benefit of each power distribution network of subordinate, be specifically expressed as:

$f_{\mathrm{Li}}\left(Q_{\mathrm{Li}}\right)=\left\{\begin{array}{c}{y}_i=a_{\mathrm{in}}{Q_{\mathrm{Li}}}^n+a_{\mathrm{in}-1}{Q_{\mathrm{Li}}}^{n-1}+...+a_{i0},i\&Element;C_{N_1}\\ \\ ...\\ y_j=\left\{\begin{array}{cc}{f}_{j1}& Q_{{\mathrm{Li}}^{\&prime;}}\&Element;S_{j0}=[l_{j0},u_{j0}]\\ f_{j2}& Q_{{\mathrm{Lj}}^{\&prime;}}\&Element;S_{j1}=[l_{j1},u_{j1}]\\ & .\\ & .\\ & .\\ & \\ & \\ f_{\mathrm{jn}}& Q_{{\mathrm{Lj}}^{\&prime;}}\&Element;S_{\mathrm{jn}}=[l_{\mathrm{jn}},u_{\mathrm{jn}}]\end{array}\right.,j\&Element;C_{N-N_1}\\ ...\end{array}\right.$

Wherein, S ₁∪ S ₂∪ ...=S _i, S _ifor the idle Q in critical point _lifeas ible space, a _im, a _im-1... a _i0for coefficient.

Step S130: according to reactive voltage coordination optimization world model of the superior and the subordinate, administers electrical network to higher level and carries out virtual controlling subregion, obtains the coordination optimization computation model of multiple virtual controlling subregion.

In one embodiment of the invention, superior control centre is optimized calculating, obtains each control partition VACE set point.

Reactive voltage coordination optimization world model of the superior and the subordinate is containing piecewise function, and model is complicated.Introduce control partition below, and provide Decoupled Model and Practical resolving Algorithm thereof, to reduce calculating scale and algorithm complex further, improve computational efficiency and reliability.

VACE(Voltage/Var Area Control Error, region REACTIVE POWER/VOLTAGE CONTROL Deviation Indices), be expressed as:

I _VACE=(Q _B-Q _Bref)+E _S(U _B-U _Bref)。

Wherein, E _s=-K _qSbe portion of district network paramter matrix, have nothing to do with its external network, Q _bfor critical point reactive power vector, U _bfor critical point voltage vector.When one's respective area does not have reactive disturbance, no matter how outside changes, and affects one's respective area I by critical point variable _vACEvalue is all zero.If control deviation all approximate zero of each element in the rear idle deviation VACE in each region, then each intra-zone reactive load change can be ignored to other regional effects outside.Therefore, by controlling the idle VACE vector of each sub regions near null value, can uneoupled control between feasible region.Retrain below all subregion AVC control system adds:

-δ≤I _i,VACE+C _i,VACEΔu' _i≤δ

Wherein, C _{i, VACE}the control sensitivity of one's respective area Controlling vertex to VACE, Δ u _i' be control deviation, δ is the VACE maximum deviation that region voltage control system allows, and must get suitable minimum, ensure existing feasible solution, also reach decoupling effect according to actual conditions.Size comparison definition in formula be vector element in order correspondence carry out numeric ratio comparatively.

In one embodiment of the invention, according to the Decoupling Characteristics of VACE, to considering that the reactive voltage coordination optimization world model of the superior and the subordinate of transmission and distribution network comprehensive benefit decomposes, obtain each control partition coordination optimization computation model.

Wherein, each control partition coordination optimization computation model is obtained:

min f _H(x',u',Q _Li,U _Hi)+∑f _Li(Q _Li)

s.t.g _i(x',u',U _Hi,Q _Li)=0

h _j(x',u',U _Hi,Q _Li)≤0

x' _min<x'<x' _max

u' _min<u'<u' _max

Q _Li，min<Q _Li<Q _Li，maxi＝1,2,…,N'

U _Hi，min<U _Hi<U _Hi，maxi＝1,2,…,N'

-δ≤I _i，VACE+C _i,VACEΔu' _i≤δ

Wherein, x ', u ' be the subset of higher level's electric network state amount x and controlled quentity controlled variable u, represent quantity of state and the controlled quentity controlled variable of i inside, control partition, U _hi, Q _lithe critical point coordination variable of subordinate's electrical network corresponding in the i of control partition, Δ u _i' be control deviation in the i of control partition.

Step S140: superior control centre, according to the coordination optimization computation model of virtual controlling subregion, calculates the optimal value of the critical point coordination variable of subordinate's electrical network and is sent to subordinate's electrical network.

Calculate the optimal value of the critical point coordination variable of subordinate's electrical network.Be optimized calculating at superior control centre, obtain the voltage U on the interconnection of each control partition _{iB, opt}with idle wherein, computation model is as follows:

min f _H(x,u,Q _Li,U _Hi)+∑f _Li(Q _Li)

s.t.g _i(x,u,U _Hi,Q _Li)=0

h _j(x,u,U _Hi,Q _Li)≤0

x _min<x<x _max

u _min<u<u _max

For control partition i, be U with the current critical point voltage vector on outside link line _iB, current critical point reactive power vector is Q _iB, region VACE coefficient matrix is E _iS, obtain the VACE vector of control partition i.Each control partition set up this control divide about corresponding coordination optimization computation model.Wherein, if if in the coordination optimization computation model of control partition there is not piecewise function in target function, namely not with piecewise function statement subordinate grid benefit characteristic, then directly obtain higher level's power grid control amount by control partition coordination optimization computation model and instruct the critical point coordination variable of subordinate's optimisation strategy.

If target function exists piecewise function in the coordination optimization computation model of control partition, then just needing to upload benefit before subregion is that the sheet net of segmented version processes, and marks separately a subregion and its coordination.For the latter, can by enumerating the optimum benefit in each feasible subspace, obtain corresponding computation model as shown in the formula:

min f _H(x',u',Q _Li,U _Hi)+f _ij

s.t.g _i(x',u',U _Hi,Q _Li)=0

h _j(x',u',U _Hi,Q _Li)≤0

x' _min<x'<x' _max

u' _min<u'<u' _max

Q _Lij∈S _ij=[l _ij,u _ij] i=1,2,…,N,1≤j≤n

U _Hi，min<U _Hi<U _Hi,maxi＝1,2,…,N'

-δ≤I _i,VACE+C _i,VACEΔu' _i≤δ

Wherein, Q _{li, min}<Q _li<Q _{li, max}be set to the subspace Q of critical point variable _lij∈ S _ij=[l _ij, u _ij], target function f _li(Q _li) become constant f from piecewise function _ij, can in the hope of the optimal solution x of this control partition i at constraint subspace j _{ij, opt}, and the optimum network loss f of correspondence _{hL, ij, opt}.Finally more all retrain the optimal solution f of subspace _{hL, i, opt}, the optimal solution in whole coordination subarea can be obtained, and then obtain the coordination subarea optimisation strategy considering the superior and the subordinate electric network synthetic benefits, and each critical point coordination variable optimal value of subordinate's power distribution network corresponding to control partition.

Superior control centre is by the optimal value calculating the critical point coordination variable of subordinate's electrical network and be sent to subordinate's electrical network.

In one embodiment of the invention, also comprise innovatory algorithm, for virtual controlling subregion is optimized calculating, integrates and generate new virtual partition, wherein, the number of new virtual partition is less than the virtual partition number before optimization.

When the interval number optimizing benefit is more, need corresponding Optimal calculation Load flow calculation number of times more, such as, containing 10 intervals, then need 10 suboptimums to calculate Load flow calculation.Therefore, according to the feature of feasible constraint space and piecewise function, propose to improve derivation algorithm, effectively can reduce control partition Optimal calculation Load flow calculation number of times.The interval non-overlapping copies of general benefit expression form, can be summarized as follows formula:

$f_{\mathrm{Li}}\left(Q_{\mathrm{Li}}\right)=\left\{\begin{array}{cc}{f}_{i1}& Q_{\mathrm{Li}}\&Element;S_{i1}=[l_{i1},u_{i1}]\\ f_{i2}& Q_{\mathrm{Li}}\&Element;S_{i2}=[l_{i2},u_{i2}]\\ & ...\\ f_{\mathrm{in}}& Q_{\mathrm{Li}}\&Element;S_{i2}=[l_{\mathrm{in}},u_{\mathrm{in}}]\end{array}\right.$

Wherein, f _ijfor namely i-th is coordinated the network loss level that subarea jth retrains subspace, S ₁∪ S ₂∪ ...=S _i, S _ifor the idle Q in critical point _lifeas ible space, and 1<=j, k<=n.General conventional network loss, as performance indicator, is carried out mass data to the relation of power distribution network when given critical point variable between the idle and minimum network loss in critical point virtual, can be obtained the relation that both are single paddy function in constraint space.

In constraint space, according to the horizontal f of the network loss of each subspace _ij, resequenced in interval, and adjacent interval merged, obtain the following modification optimizing interval expression benefit:

Changed in each interval of constraint space further and optimize the interval form expressing the amendment of benefit, the meaning of interval representative becomes: the more excellent benefit that can reach in interval.Along with the increase that interval rank is numbered, interval is progressively relaxed, and network loss level declines.

State according to revised benefit, computation model below adopting in this control partition:

min f _H(x',u',Q _Li,U _Hi)

s.t.g _i(x',u',U _Hi,Q _Li)=0

h _j(x',u',U _Hi,Q _Li)≤0

x' _min<x'<x' _max

u' _min<u'<u' _max

Q _Lij∈S' _ij=[l' _ij,u' _ij] i=1,2,…,N,1≤j≤n

U _Hi，min<U _Hi<U _Hi，maxi＝1,2,…,N'

-δ≤I _i,VACE+C _i,VACEΔu' _i≤δ

Innovatory algorithm for interval statement comprises the steps:

Step S141: select constraint subspace (i.e. staff cultivation space) S _in' start to calculate.

Step S142: make the constraint subspace of current calculating be S _ij', ask following computation model:

min f _H(x',u',Q _Li,U _Hi)

s.t.g _i(x',u',U _Hi,Q _Li)=0

h _j(x',u',U _Hi,Q _Li)≤0

x' _min<x'<x' _max

u' _min<u'<u' _max

Q _Lij∈S' _ij=[l' _ij,u' _ij] i=1,2,…,N,1≤j≤n

U _Hi，min<U _Hi<U _Hi，maxi＝1,2,…,N'

-δ≤I _i,VACE+C _i,VACEΔu' _i≤δ

Wherein, x ', u ' be the subset of higher level's electric network state amount x and controlled quentity controlled variable u, U is critical point voltage, and Q is that critical point is idle, min and max represents corresponding minimum and maximum threshold value respectively, I _{i, VACE}for retraining the VACE index of subspace i, C _{i, VACE}retrain subspace i to the control sensitivity of VACE, Δ u _i' be control deviation, δ is the VACE maximum deviation that region voltage control system allows.Calculate locally optimal solution f _{ij, opt}', can obtain coordinating subarea i and retrain subspace S _in' locally optimal solution drop point x _ij' and locally optimal solution f _{hL, ij, opt}, wherein,

f _HL,ij,opt=f' _ij,opt+f' _ij

Step S143: if drop point belongs to other subspaces S _ik', then basis relation, obtain subspace S _ik' locally optimal solution f _{hL, ik, opt}

f _H,ik，opt=f' _ij,opt+f' _ik

If drop point x _ij' belong to multiple space, interval rank is numbered minimum interval numbering and be designated as N _tEMP.Due to f ' _ik<f' _ij, interval rank numbering N _tEMPlocally optimal solution f _{hL, ik, opt}for the solution that the network loss in above-mentioned multiple space is relatively minimum.

Step S144: if first time calculates, remember interval rank numbering N _tEMPlocally optimal solution f _{hL, ik, opt}for f _{hL, i, temp}, and remember current drop point x _ij' be x _{i, temp}.Once calculated if carried out, more current f _{hL, ik, opt}with the f of record _{hL, i, temp}if meet f _{hL, ik, opt}<f _{hL, i, temp}, then f is made _{hL, i, temp}=f _{hL, ik, opt}, x _{i, temp}=x _ij'; Otherwise, if meet f _{hL, i, temp}<f' _{ij, opt}+ f _i1', then innovatory algorithm iteration terminates.

Step S145: if interval numbering N _tEMP+ 1 is non-vanishing, gets corresponding constraint subspace perform step S2.

If interval numbering N _tEMP+ 1 is zero, then innovatory algorithm iteration terminates.

Step S146: single control partition i calculates end, obtains the optimal solution f of control partition district i _{hL, i, opt}with optimal solution drop point x _i, wherein, f _{hL, i, opt}=f _{hL, i, temp}, x _i=x _{i, temp}.

Above-mentioned innovatory algorithm, according to the feature of feasible constraint space and piecewise function, devises Practical resolving Algorithm, Optimal calculation Load flow calculation number of times in subregion is reduced.

It should be noted that the interval expression formula of benefit above adjusts benefit to be the hypothesis of single paddy trend at feas ible space based on general distribution.When not being single paddy, by judging the new drop point that point paddy last rank interval computation obtains, the iteration of other paddy can be carried out after completing the iteration of this drop point place paddy again, then result compared.

Step S150: subordinate's electrical network, according to the optimal value of critical point coordination variable, adds coordination variable constraint, obtains the virtual controlling partition zone optimizing strategy based on transmission and distribution network comprehensive benefit in AVC algorithm.By each critical point coordination variable set point U _{hi, ref}and Q _{li, ref}be handed down to control centre of each corresponding subordinate, All Control Center application critical point voltage U _{hi, ref}, and in respective AVC algorithm, add the idle constraint in tracking critical point, thus obtain the virtual controlling partition zone optimizing strategy based on transmission and distribution network comprehensive benefit further.

Step S160: higher level's electrical network and subordinate's electrical network carry out reactive voltage cooperation control according to virtual controlling partition zone optimizing strategy.

Be further described feasibility of the present invention and advantage below by a simulation example, be understandable that, only for illustrative purposes, embodiments of the invention are not limited thereto described simulation example.

This example adopts IEEE39 node system to carry out simulation calculation.First demonstrate and introduce the feasibility that VACE index realizes uneoupled control between control partition, and the result of global optimization and coordination optimization is contrasted, then demonstrate respectively when single power distribution network accesses and when multiple power distribution network accesses, the feasibility of the main website coordination optimization control strategy that the present invention proposes, traditional is not considered power distribution network critical point ability, algorithm and the algorithm of the present invention of demand and optimization benefit contrast, finally consider that multiple power distribution network accesses, when load or burden without work increases continuously, the implementation effect of the algorithm that traditional algorithm and the present invention propose is contrasted.

The decoupling zero of higher level's power grid control subregion is realized based on VACE index.IEEE39 node system is adopted to carry out simulation calculation.Be preset with 17 interconnections, 1-2,3-4,16-17 tri-interconnections cut, obtain two subregions of two IEEE39 systems, containing node { 2,3,17,18,25,26,27,28,29,30,37,38} is designated as subsystem A, containing node { Isosorbide-5-Nitrae, 5,6,7,8,9,10,11,12,13,14,15,16,19,20,21,22,23,24,31,32,33,34,35,36,39} is subsystem B, as shown in Figure 2.

Be optimized calculating at total system, can be optimized result, P _loss=0.27288; Two systems carry out independent optimization respectively, obtain respective optimum results: P _aLoss=0.11937, P _bLoss=0.15438.Now, subsystem network loss summation is, P _loss'=P _aLoss+ P _bLoss=0.27375

Visible, the optimum results obtained before and after subregion is close, can think that global optimization is consistent with the result of coordination optimization.According to other cutting mode, such as, by { 1,3,9} is included into subsystem A, also can obtain similar results.Prove that the subregion realized based on VACE does not affect the control effects of optimisation strategy.

Coordination optimization control strategy when single power distribution network access.The feasibility of main website coordination optimization control strategy when verifying that single power distribution network accesses below.Choose five kinds of algorithms to compare: the algorithm not considering power distribution network critical point ability, demand and optimization benefit; Only consider the algorithm of power distribution network critical point ability and demand; Consider power distribution network critical point ability and demand, the algorithm being with curve statement benefit; Consider power distribution network critical point ability and demand, the algorithm being with interval statement benefit; With the innovatory algorithm that interval statement benefit is.Wherein, the algorithm that the present invention proposes is latter two algorithm.

For IEEE39 node system.No. 17 nodes are coordinated critical point as transmission and distribution network, and namely power transmission network root node carries out load equivalence at No. 17 nodes to power distribution network.Coordination variable selects power transmission network root node voltage U _hidle Q is injected with power distribution network critical point _l, connect being assumed to by beginning regulating power and demand for control of power distribution network:

1.00<U _H<1.07

-8<Q _L<-1

Algorithm effect for convenience of description, makes the optimum benefit of power distribution network be about the half of the optimum benefit of power transmission network.Grid benefit is assessed by network loss, is set to

f _Distr=0.1475-0.01Q _L+0.01Q _L ²

In order to mutual verification computation result, be the section process to above formula with the function of range format statement benefit, i.e. following formula:

$f_{\mathrm{Distr}}=\left\{\begin{array}{cc}0.785,& Q_L\&Element;(-8,-7]\\ 0.635,& Q_L\&Element;(-7,-6]\\ 0.505,& Q_L\&Element;(-6,-5]\\ 0.395,& Q_L\&Element;(-5,-4]\\ 0.305,& Q_L\&Element;(-4,-3]\\ 0.235,& Q_L\&Element;(-3,-2]\\ 0.185,& Q_L\&Element;(-2,-1]\end{array}\right.$

That note does not adjust critical point ability, demand and optimization benefit with not considering is algorithm A ₀; That only considers power distribution network critical point ability and demand is designated as algorithm A ₁; Consider power distribution network critical point ability and demand, be designated as algorithm A with curve statement benefit form ₂; Algorithm A is designated as with interval statement benefit form ₃; A is designated as with the innovatory algorithm of interval statement benefit form ₄.Result of calculation is shown in that IEEE39 node system transmission and distribution network coordinates example table 1.

Table 1

As seen from Table 1, algorithm A ₀although the result obtained is best, because not considering demand and the constraint of subordinate's electrical network, critical point voltage gets over the receivable upper limit of subordinate's electrical network, if implement this optimisation strategy will cause subordinate's power distribution network key node voltage out-of-limit.Algorithm A ₃, A ₄in distribution network loss only represent the level of the attainable optimization network loss in each interval, only for calculating power transmission network optimisation strategy, obtain critical point optimization solution, actual distribution network loss still needs to issue value depending on critical point is idle.Algorithm A ₁, A ₂, A ₃, A ₄result is close, and table 1 fails to show its nuance because retaining result length.Because the expression function of subordinate's grid net loss in its constraint space without the minimum point Q of the idle coordination constraint in its critical point _l=0.5, below demand for control and regulating power are modified, make constraint space toward the direction continuation of subordinate's grid net loss minimum point.

$f_{\mathrm{Distr}}=\left\{\begin{array}{cc}0.785,& Q_L\&Element;(-8,-7]\\ 0.635,& Q_L\&Element;(-7,-6]\\ 0.505,& Q_L\&Element;(-6,-5]\\ 0.395,& Q_L\&Element;(-5,-4]\\ 0.305,& Q_L\&Element;(-4,-3]\\ 0.235,& Q_L\&Element;(-3,-2]\\ 0.185,& Q_L\&Element;(-2,-1]\\ 0.155,& Q_L\&Element;(-\mathrm{1,0}]\\ 0.145,& Q_L\&Element;(\mathrm{0,1}]\\ 0.155,& Q_L\&Element;\left(\mathrm{1,2}\right]\end{array}\right.$

1.00<U _H<1.07

Wherein,

-8<Q _L<2

Recalculate according to above-mentioned amendment below, obtain the table 2 that IEEE39 node system transmission and distribution network coordinates example:

Table 2

After constraint space comprises subordinate's grid net loss minimum point, contrast table 2 can find, only considers the algorithm A of the algorithm of subordinate's ability and demand ₁result there occurs change.This algorithm further considers the algorithm A of subordinate's benefit at that point with other ₂, A ₃, A ₄compare, although Losses is less, subordinate's distribution network loss is higher, and result is that the total network loss of total system increases to some extent, reflects the superior and the subordinate's optimisation strategy conflict that such algorithm is implicit.Algorithm A ₂, A ₃, A ₄consider subordinate's power distribution network benefit, although therefore Losses increases to some extent, the overall network loss of total system is optimized.

Whether comprise in twice calculating in the minimum front and back of subordinate's grid net loss in adjustment, contrast algorithm A ₃with algorithm A ₄, the optimum results before and after visible discrete model amendment is basically identical.Such as, A ₃losses be 0.281296, A ₄be 0.281298.But variant in computational speed, such as, in twice example, A ₃needing to calculate all interval, as having 7 in the 1st example, in the 2nd example, having 10, and A ₄only need calculate 2 and 3.

Coordination optimization control strategy when multiple power distribution network access.When verifying that multiple power distribution network accesses below, the feasibility of main website coordination optimization control strategy proposed by the invention.Have chosen three kinds of algorithms in sample calculation analysis to compare: the algorithm A not considering power distribution network critical point ability, demand and optimization benefit ₀', wherein, A ₀' be algorithms most in use; Only consider the algorithm A of power distribution network critical point ability and demand ₁'; Consider the algorithm A of power distribution network critical point ability, demand and benefit ₂', method namely proposed by the invention.Wherein, algorithm A ₂' be carry out calculating each control partition VACE set point at total system, then carry out control partition and optimize calculating.Each control partition forms corresponding control constraints based on this VACE set point, forms coordination constraint and target correction term based on critical point characteristic, thus sets up control partition coordination optimization computation model, calculates and obtain each subordinate power distribution network coordination variable to issue value.

Below for IEEE39 node system, be split as two subsystems A and B as shown in Figure 2.Consider that simple two power distribution networks are linked into the situation of power transmission network, wherein No. 1 power distribution network is connected to No. 17 critical points, and power distribution network AVC ability is higher, can provide demand for control and regulating power and can state its benefit with curve.No. 2 power distribution networks are connected to No. 13 critical points, and No. 2 power distribution network AVC abilities are slightly not enough, can provide demand for control and regulating power, state its benefit with interval.Subsystem A and No. 1 power distribution network composition control subregion A1, subsystem B and No. 2 power distribution network composition control subregion B2.And the total system that note does not comprise subordinate's power distribution network is S, the total system comprising subordinate's power distribution network is S'.By calculating IEEE39 node system transmission and distribution network tuning algorithm contrast table 3.

Table 3

In order to make three kinds of algorithm network loss possess comparativity, make the computational methods of each distribution network loss consistent, optimize respectively hanging oneself calculating try to achieve critical point inject idle after, substitute into f further _distr=0.1475-0.01Q _l+ 0.01Q _l ²calculate.

As can be seen from Table 3, algorithm A ₀' owing to not considering subordinate's power distribution network situation, although Losses is optimum, overall total system network loss is comparatively large, and subordinate's electrical network exists the out-of-limit situation of node voltage; Algorithm A ₁' consider control ability and the demand for control of subordinate's electrical network, total system voltage security can be ensured, but owing to reckoning without subordinate's grid net loss, so overall network loss is larger; Algorithm A ₂' introduce consideration to subordinate's grid net loss, although Losses is a little more than algorithm A ₁', but have a clear superiority in overall total system S ' network loss.

Each algorithm implementation effect when load or burden without work increases continuously.Consider that multiple power distribution network accesses below, when load or burden without work increases continuously, compare the implementation effect of traditional algorithm and algorithm proposed by the invention.Sample calculation analysis comprises the algorithm not considering power distribution network critical point ability, demand and optimization benefit; Only consider the algorithm of power distribution network critical point ability and demand; Consider the algorithm of power distribution network critical point ability, demand and benefit.

For IEEE39 node system, consider that simple two power distribution networks are linked into power transmission network situation, wherein No. 1 power distribution network is connected to No. 18 critical points, and No. 2 power distribution networks are connected to No. 13 critical points.Algorithm A ₀', A ₁', A ₂' implication is the same; No. 1, No. 2 power distribution network critical point characteristics are also same as above.Make system from the 2nd control cycle, namely during t=2, load bus load or burden without work by 1% sustainable growth, below to the transmission and distribution network network loss of each algorithm in 14 control cycles, critical point is idle, critical point voltage is analyzed.

Can be seen by Fig. 3, fall damage angle from total system, algorithm A ₂' good compared with other algorithms.As shown in Figure 4, fall in power transmission network side and damage in effect angle, algorithm A ₂' slightly poor compared with other algorithms.Reason can be contrasted by the total network loss of Fig. 5 power distribution network and learn, algorithm A ₂' consider power distribution network benefit, after coordination optimization controls, the network loss connecing power distribution network under No. 13 critical points and No. 18 critical points is all optimized.Figure Fig. 6 to Fig. 9 gets access to grid under illustrating power transmission network side sub-system A, B and No. 13, No. 18 network loss situation.Figure 10 and Figure 11 shows the contrast of two power distribution network critical point voltages.As we can see from the figure, algorithm A ₀' under not having to consider, connecing the voltage requirements of power distribution network, critical point voltage, by general node processing, only controls 0.9 to 1.1, connects power distribution network builtin voltage out-of-limit under may causing.Algorithm A ₁', A ₂' then according under connect distribution network voltage demand, well voltage control in 1.0 to 1.07 scopes.Figure 12 and Figure 13 shows two power distribution network critical points and injects idle contrast.As we can see from the figure, algorithm A ₀' under not having to consider, connect the Reactive-power control ability of power distribution network, only think that it is adjustable, under the 13rd, 14 control cycle Optimal Setting values exceed, connect power distribution network Reactive-power control lower limit; Algorithm A ₁', A ₂' consider down and connect power distribution network critical point and inject Reactive-power control ability, the optimization issued is calculated set value calculation in the reasonable scope.

To sum up, the algorithm that the present invention proposes effectively can improve system losses, and consider subordinate to the demand for control of critical point voltage and Reactive-power control ability, thus ensure transmission and distribution network both sides reactive balances, avoid because critical point voltage is inclined or the power distribution network internal node that causes on the low side is out-of-limit.Owing to have employed zoning concept, when a certain coordination subarea optimal load flow result calculate disperse or optimum results unavailable, be can not have influence on other to coordinate the calculating in subareas and even control; When adjusting the AVC system or algorithm that do not have or only have insufficiency of function in multiple coordination subarea, when needing to adopt piecewise function to express network loss level, total algorithm performance also can not be made to decline to a great extent, and namely algorithm has certain compatibility.

Consider subordinate's operation of power networks demand, ability and benefit according to the reactive voltage control method for coordinating based on transmission and distribution network comprehensive benefit of the embodiment of the present invention, show the running status improving the whole network, improve the system benefit of the whole network;

The algorithm that the present invention proposes, owing to improving the running status of the whole network, improves the system benefit of the whole network; In algorithm, control partition is carried out to higher level's power transmission network, and computation model and the Practical resolving Algorithm of control partition are proposed, by bulk power grid reactive voltage Harmonic Control decomposition and inversion be dimensionality reduction control partition in the superior and the subordinate's electric network reactive-load voltage coordination problem, reduce computation model scale and computation complexity; Simultaneously subregion uneoupled control improves the reliability of AVC algorithm, if certain subregion occurs abnormal, such as optimal load flow disperse or optimum results unavailable, the normal coordination optimization that also can not have influence on other subregions controls.Take into full account that subordinate's benefit is uploaded in addition in algorithm multi-form, taken into full account that subordinate uploads and optimize the multi-form of benefit, improve the compatibility of traditional algorithm to the benefit expression-form that subordinate uploads, reduce the calculation cost that subordinate's electrical network is extra.

In the description of this specification, specific features, structure, material or feature that the description of reference term " embodiment ", " some embodiments ", " example ", " concrete example " or " some examples " etc. means to describe in conjunction with this embodiment or example are contained at least one embodiment of the present invention or example.In this manual, identical embodiment or example are not necessarily referred to the schematic representation of above-mentioned term.And the specific features of description, structure, material or feature can combine in an appropriate manner in any one or more embodiment or example.

Although illustrate and describe embodiments of the invention, for the ordinary skill in the art, be appreciated that and can carry out multiple change, amendment, replacement and modification to these embodiments without departing from the principles and spirit of the present invention, scope of the present invention is by claims and equivalency thereof.

Claims (8)

1., based on a reactive voltage control method for coordinating for transmission and distribution network comprehensive benefit, it is characterized in that, comprise the following steps:

Electrical network superior control centre of subordinate uploads critical point characteristic, and wherein, described critical point characteristic comprises regulating power, demand for control and optimization benefit;

Described critical point characteristic is converted into the coordination variable constraint of critical point variable sign and target correction term by described superior control centre, sets up reactive voltage coordination optimization world model of the superior and the subordinate;

According to reactive voltage coordination optimization world model of described the superior and the subordinate, electrical network is administered to higher level and carries out virtual controlling subregion, obtain the coordination optimization computation model of multiple virtual controlling subregion;

Described superior control centre, according to described coordination optimization computation model, calculates the optimal value of the critical point coordination variable of described subordinate electrical network and is sent to described subordinate electrical network;

Described subordinate electrical network, according to the optimal value of described critical point coordination variable, adds coordination variable constraint, obtains the virtual controlling partition zone optimizing strategy based on transmission and distribution network comprehensive benefit in automatism voltage control AVC algorithm; And

Higher level's electrical network and described subordinate electrical network carry out reactive voltage cooperation control according to described virtual controlling partition zone optimizing strategy.

2. as claimed in claim 1 based on the reactive voltage control method for coordinating of transmission and distribution network comprehensive benefit, it is characterized in that, described regulating power is for representing that described subordinate electrical network is to the idle controllable scope in described critical point.

3., as claimed in claim 1 based on the reactive voltage control method for coordinating of transmission and distribution network comprehensive benefit, it is characterized in that, described demand for control is qualified for controlling described subordinate line voltage, and the restriction range that shown critical point voltage needs higher level to control.

4. as claimed in claim 1 based on the reactive voltage control method for coordinating of transmission and distribution network comprehensive benefit, it is characterized in that, described optimization benefit is when the described virtual controlling partition zone optimizing strategy of control does not clash in benefit, each subordinate grid benefit considered in coordination optimization computation model.

5. as claimed in claim 1 based on the reactive voltage control method for coordinating of transmission and distribution network comprehensive benefit, it is characterized in that, described optimization benefit, according to subordinate electrical network AVC, with described critical point variable for independent variable, is stated by curve form or range format.

6. as claimed in claim 1 based on the reactive voltage control method for coordinating of transmission and distribution network comprehensive benefit, it is characterized in that, also comprise innovatory algorithm, for virtual controlling subregion is optimized calculating, integrate and generate new virtual partition, wherein, the number of described new virtual partition is less than the virtual partition number before optimization.

7., as claimed in claim 1 based on the reactive voltage control method for coordinating of transmission and distribution network comprehensive benefit, it is characterized in that, described superior control centre is optimized calculating, obtains each control partition VACE set point.

8. as claimed in claim 6 based on the reactive voltage control method for coordinating of transmission and distribution network comprehensive benefit, it is characterized in that, described virtual controlling subregion be optimized calculating, integrate and generate new virtual partition and comprise further:

Step 1: select constraint subspace S _in' start to calculate;

Step 2: the described constraint subspace making current calculating is S _ij', solve following computation model:

minf _H(x',u',Q _Li,U _Hi)

s.t.g _i(x',u',U _Hi,Q _Li)＝0

h _j(x',u',U _Hi,Q _Li)≤0

x' _min<x'<x' _max

u' _min<u'<u' _max

Q _Lij∈S' _ij＝[l' _ij,u' _ij]i＝1,2,…,N,1≤j≤n

U _Hi,min<U _Hi<U _Hi,maxi＝1,2,…,N'

-δ≤I _i,VACE+C _i,VACEΔu' _i≤δ

Wherein, x ', u ' be the subset of higher level's electric network state amount x and controlled quentity controlled variable u, U is described critical point voltage, and Q is that described critical point is idle, min and max represents corresponding minimum and maximum threshold value respectively, f _hfor target function, g _i=0 is equality constraint, h _j≤ 0 is inequality constraints, u' _ij, l' _ijfor control partition S _icritical point is idle Q _lifeas ible space bound, I _{i, VACE}for the VACE index of described constraint subspace i, C _{i, VACE}for described constraint subspace i is to the control sensitivity of VACE, Δ u _i' be control deviation, δ is the VACE maximum deviation that region voltage control system allows, and calculates locally optimal solution f _{ij, opt}', can obtain coordinating subarea i and retrain subspace S _in' locally optimal solution drop point x _ij' and locally optimal solution f _{hL, ij, opt}, wherein,

f _HL,ij,opt＝f' _ij,opt+f' _ij

Step 3: if drop point x _ij' belong to other subspaces S _ik', then basis relation, obtain subspace S _ik' locally optimal solution f _{hL, ik, opt}, wherein, f _{hL, ik, opt}=f' _{ij, opt}+ f' _ik;

If described drop point x _ij' belong to multiple subspace, then interval rank is numbered minimum interval numbering and be designated as N _tEMP; Due to f' _ik<f' _ij, described interval rank numbering N _tEMPlocally optimal solution f _{hL, ik, opt}for the solution that the network loss in described multiple space is relatively minimum;

Step 4: if first time calculates, remember described interval rank numbering N _tEMPdescribed locally optimal solution f _{hL, ik, opt}for f _{hL, i, temp}, and remember current drop point x _ij' be x _{i, temp}if carried out once calculating, more described f _{hL, ik, opt}with described f _{hL, i, temp}if meet f _{hL, ik, opt}<f _{hL, i, temp}, then f is made _{hL, i, temp}=f _{hL, ik, opt}, x _{i, temp}=x _ij'; Otherwise, if meet f _{hL, i, temp}<f' _{ij, opt}+ f _i1', then innovatory algorithm iteration terminates;

Step 5: if described interval numbering N _tEMP+ 1 is non-vanishing, gets corresponding described constraint subspace perform step 2;

If described interval numbering N _tEMP+1be zero, then innovatory algorithm iteration terminates;

Step 6: single described control partition i calculates and terminates, and obtains the described optimal solution f of described control partition i _{hL, i, opt}with optimal solution drop point x _i, wherein, f _{hL, i, opt}=f _{hL, i, temp}, x _i=x _{i, temp}.