CN106451414A - Dynamic partition-based large power grid emergency state control auxiliary decision method - Google Patents

Abstract

The present invention discloses a dynamic partition-based large power grid emergency state control auxiliary decision method, and belongs to the electric power system operation control technology field. The method of the present invention is characterized by automatically obtaining the feasible dynamic partition schemes, and aiming at the plurality of feasible schemes, quantitatively evaluating an index system according to the dynamic partition, thereby obtaining an optimal control scheme. The method cooperatively considers the economy, the reliability and the safety, gives out the optimal dynamic partition control scheme rapidly and accurately aiming at a current state and the single and various emergency states of a power grid under the preconceived faults, can give out the correct auxiliary decision suggestions, and guarantees the safe and stable operation of the power grid.

Description

A kind of bulk power grid emergency state control aid decision-making method based on dynamic partition

Technical field

The invention belongs to technical field of power system operation control, it particularly relates to a kind of big based on dynamic partition Electrical network emergency state control aid decision-making method.

Background technology

Put into operation with building up of extra-high voltage alternating current-direct current circuit, the general layout of the transregional transmission of electricity of Large Copacity gradually forms.Transregional interconnection Electrical network makes the dynamical couple characteristic between each regional power grid further enhance, and partial electric grid fault is to large regional grid security and stability Impact is bigger.When electrical network enters the state of emergency, if not implementing in time to take control measure it would be possible to cause cascading failure, make Become the generation of large-scale blackout.

The operation control device of electrical network classification two class from implementation at present, a class is to be gone out by changing generator Power and workload demand dynamically adjust power system operating mode, and a class is to change trend distribution using power electronic equipment.Cut with the machine of cutting The two traditional three lines of defence measures such as load, off-the-line are compared, and load transfer only changes the topological structure of network in a small range, just Can largely reduce to subtract and supply load proportion, effect is significant, reliability is high, has caused the pass of domestic and international researcher Note.

In this context, can be using the dynamic partition skill including load transfer when electrical network is in a state of emergency Art, carries out Reasonable adjustment, and then controls the trend of related elements, voltage etc., make electrical network exit urgent shape to partition network structure State is to keep safety and stability to have certain meaning in holding electrical network.Currently for apparatus overload, voltage out-of-limit problem, engineering should Use and how to be solved by measures such as generator output adjustment, switching capacitive reactances.It is contemplated that from network topology structure, Introduce Dynamic partitioning techniques, provide the dynamic partition Adjusted Option coordinating the targets such as security, economy, to solve equipment mistake The unitary class state of emergency such as load, voltage out-of-limit, safeguards the stable operation of electrical network.

Content of the invention

The present invention seeks to：A kind of bulk power grid emergency state control aid decision-making method based on dynamic partition is provided, with Realize when the problems such as voltage out-of-limit, apparatus overload in electrical network under current state or forecast failure, start dynamic partition in time Aid decision calculates, and provides correct control strategy suggestion for dispatcher.

Specifically, the present invention employs the following technical solutions realization, comprises the following steps：

1) according to power system operation mode and corresponding model and parameter, calculated online, after obtaining forecast failure System load flow data, enter step 2)；

2) Voltage security margin is less than and sets threshold value η_v,thr1It is less than with Voltage security margin after forecast failure and set door Threshold value η_v,thr2Bus be added separately to out-of-limit bus collection and bus collection to be controlled, by overload nargin be less than set threshold value η_p,thr1It is less than with overload nargin after forecast failure and set threshold value η_p,thr2Equipment be added separately to overload prevention device collection and Apparatus to be controlled collection, wherein out-of-limit bus collection and bus collection to be controlled are referred to as bus and investigate collection, overload prevention device collection and apparatus to be controlled collection The equipment that is referred to as investigates collection；

If out-of-limit bus collection and overload prevention device collection are sky, there is not apparatus overload under electrical network current state or voltage is got over Limit, then method ends；Otherwise, enter step 3) the emergency state control aid decision calculating based on dynamic partition for the startup；

Described equipment refers to circuit and transformer；

Described dynamic partition refers to, when bulk power grid enters the state of emergency, by carrying out switching to part circuit, change Topological structure of electric, so that the state of emergency is mitigated or eliminated, including the solution cyclization of load transfer and electromagnetic looped network；

3) set the dynamic partition scheme switching number of lines, screening can switching circuit, formed alternative space, enter step Rapid 4)；

4) if out-of-limit bus collection and the equal non-NULL of overload prevention device collection, show under current state electrical network occur the multiclass state of emergency and There is apparatus overload and voltage out-of-limit simultaneously, be then directed to each the dynamic partition scheme in alternative space, calculate it comprehensive Performance indications R_e,i, integrated performance index is more than and sets threshold value R_e,i,threScheme add to treat check scheme concentrate, enter step Rapid 5)；

If out-of-limit bus collection is empty and overload prevention device collection non-NULL, show that under current state, electrical network occurs the urgent shape of unitary class , there is apparatus overload in state, then be directed to each the dynamic partition scheme in alternative space, calculate its integrated performance index ω_p,i, integrated performance index is more than and sets threshold value ω_p,i,threScheme add to treat check scheme concentrate, enter step 5)；

If out-of-limit bus collection non-NULL and overload prevention device collection are sky, that is, under current state there is the unitary class state of emergency in electrical network, Have that busbar voltage is out-of-limit, be then directed to each the dynamic partition scheme in alternative space, calculate its integrated performance index ω_u,i, integrated performance index is more than and sets threshold value ω_u,i,threScheme add to treat check scheme concentrate, enter step 5)；

R_e,i、ω_p,iAnd ω_u,i, for characterizing i-th kind of dynamic partition scheme to institute under electrical network current state or forecast failure The improvement degree of the specific state of emergency occurring, its value is bigger, and effect is better；

5) if treating check scheme collection non-NULL, for each the dynamic partition scheme treating that check scheme is concentrated, using exchange tide After stream calculation dynamic partition scheme is implemented, bus investigates the Voltage security margin concentrating each bus and the equipment investigation each equipment of concentration Overload nargin, if the overload nargin of the Voltage security margin of each bus and equipment be all higher than set threshold value, should Scheme adds feasible program to concentrate；Enter step 6)；

If treating, check scheme collection is sky, judges whether switching circuit number has reached the upper limit, if the not up to upper limit, increases Circuit switching bar number, redefines alternative space, enters step 3)；If switching circuit number has reached the upper limit, terminate This method；

6) if feasible program collection non-NULL, for multiple dynamic partition schemes of feasible program concentration, based on dynamic partition Quantitatively evaluating index system, carries out overall merit, selects optimum control scheme, method ends；If feasible program collection is sky, Each scheme treating check scheme concentration is ranked up by the size of integrated performance index, takes maximum the moving of integrated performance index State partition scheme is as last control program, method ends.

Technique scheme is further characterized by, described step 2) in, calculate bus according to formula (1)～(2) Voltage security margin：

In formula：η_v,iVoltage security margin for bus i；U_iVoltage magnitude for bus i；WithIt is respectively its upper limit And lower limit；

Described equipment includes circuit and transformer, calculates the overload of circuit and transformer respectively according to formula (3) and (4) Margin of safety：

In formula：η_pl,iOverload nargin for circuit i；I_liCurrent value for circuit i；I_li.limThermally-stabilised for circuit i Limit value；

In formula：η_pt,iOverload nargin for transformer i；S_tiApparent energy for transformer i；S_ti.limFor transformer i Rated capacity.

Technique scheme is further characterized by, described step 3) in, according to the alternative of switching circuit formation Space, need to follow following principle：Switching for 220kV circuit；The circuit putting into is to be in the line stopped transport and overhaul done state Road；Folding scheme need to ensure power supply reliability it is ensured that isolated network will not be formed；Need to arrange switching circuit according to electrical network actual conditions The upper limit.

Technique scheme is further characterized by, described step 4) in, R is calculated according to formula (5)_e,i：

In formula：R_e,iIt is the integrated performance index of i-th kind of dynamic partition scheme for voltage out-of-limit and apparatus overload；I= 1,2 ..., I, I investigate the number of devices concentrated, j=1 for equipment, and 2 ..., J, J investigate the bus number concentrated, η ' for bus_p,iWith η′_v,jIt is respectively the apparatus overload margin of safety after dynamic partition scheme is implemented and busbar voltage margin of safety, k_pvFor overload Conversion coefficient between nargin and Voltage security margin, shows the attention degrees different with voltage out-of-limit to apparatus overload；In phase Equipment is answered to be η during circuit_p,iEqual to η_pl,i, the η when relevant device is for transformer_p,iEqual to η_pt,i；

ω is calculated according to formula (6)_u,i：

In formula：ω_u,iIt is the integrated performance index of i-th kind of dynamic partition scheme for voltage out-of-limit；V investigates for bus Collection, including out-of-limit bus collection and bus collection to be controlled；γ_vImportance modifying factor for different electric pressure buses；N is current shape The bus number of voltage out-of-limit, λ under state_i,j,vVoltage security margin for wherein j-th bus；S_i,j,vDynamic for implementing i-th kind The Voltage security margin variable quantity of j-th bus after partition scheme；W is the sum of forecast failure, N_kFor under wherein k-th fault There is the bus number of voltage out-of-limit, λ_i,l,k,vVoltage security margin for l-th bus under k-th fault；S_i,l,k,vFor i-th Plant the variable quantity of the Voltage security margin of l-th bus under k-th fault after dynamic partition scheme is implemented；

ω is calculated according to formula (7)_p,i：

In formula：ω_p,iIt is the integrated performance index of i-th kind of dynamic partition scheme for apparatus overload；V investigates for equipment Collection, that is, the equipment intersection transshipped；γ_vImportance modifying factor for different electric pressure equipment；N is overload under current state Equipment number, θ_i,j,vOverload nargin for wherein j-th equipment；μ_i,j,vFor implementing jth after i-th kind of dynamic partition scheme The variable quantity of the overload nargin of individual equipment；W is the sum of forecast failure, N_kFor there is overload under wherein k-th fault Equipment number, θ_i,l,k,vOverload nargin for l-th equipment under k-th fault；μ_i,l,k,vFor i-th kind of dynamic partition scheme The variable quantity of the overload nargin of l-th equipment under k-th fault after enforcement.

Technique scheme is further characterized by, described step 6) in, based on dynamic partition quantitatively evaluating index body The dynamic partition scheme that system concentrates to feasible program carries out overall merit, carries out in the following manner：

For each scheme, it is first determined the weight of its index at different levels, and then obtained by certain mathematical algorithms The comprehensive evaluation index of upper level index, successively calculates upwards, finally obtains the comprehensive evaluation result of dynamic partition scheme；Dynamically Zonal quantization assessment indicator system includes security, reliability and economy three class index, as follows：

The first kind, safety indexes：

Line load index x_L

In formula：I=1,2 ..., n_l；η_pl,iOverload nargin for circuit i；n_lFor the number of lines；α_liImportant for circuit Property modifying factor；

Transformer load index x_T

In formula：I=1,2 ..., n_t；η_pt,iOverload nargin for transformer i；n_tFor transformer number；α_tiFor transformation Device importance modifying factor；

Voltage levvl index λ_U

In formula：I=1,2 ..., n_u；n_uFor bus quantity；η_v,iVoltage security margin for bus i；α_uiFor busbar voltage Importance modifying factor；

Short circuit current level index x_I

In formula：I=1,2 ..., n_sc；n_scFor short trouble number；I_i,scFor the short-circuit current value under i-th short trouble； I_i,sc.limThe minimum of a value of the rupturing capacity of the breaker being connected for short-circuit busbar；

Equations of The Second Kind, reliability index：

Current supply circuit number S_M

In formula：S_iFor being not more than the capacity of the 220kV end transformer station of two supplying channels；

Spinning reserve index x_re

In formula：P_SRFor active spinning reserve capacity；P_SR.limMinimum limit value for active spinning reserve capacity；

3rd class, economic index：

System losses index S_loss

S_loss=∑ P_G-∑P_D(14)

In formula：P_GFor generating capacity；P_DFor load capacity；

Operation cost index F_c

In formula：N is the number of times of circuit folding；f_iThe operating cost producing for each folding, cost and the use with switch Life-span is relevant.

Beneficial effects of the present invention are as follows：The present invention has taken into full account that under electrical network current state and forecast failure, electrical network may The voltage out-of-limit occurring and apparatus overload problem, relate to Corrective control and prevention and control.In addition, the present invention be based on dynamic Partitioning technique (includes falling load, the solution cyclization of electromagnetic looped network), and with circuit switching as core, dynamic adjustment runs the topology of rack Structure, speed is fast and controls cost little.Finally, the present invention dynamically divides to feasible according to dynamic partition quantitatively evaluating index system Area's control program carries out preferentially, considering the impact to electric network security, economy and reliability for each scheme, selecting optimal Control program.

Brief description

Fig. 1 is the FB(flow block) of the present invention.

Fig. 2 is the quantitatively evaluating index system of dynamic partition scheme.

Specific embodiment

Simultaneously referring to the drawings the present invention is described in further detail with reference to embodiment.

Embodiment 1：

One embodiment of the present of invention, its step is as shown in Figure 1：

In Fig. 1, step 1 describes, and according to power system operation mode and corresponding model and parameter, is counted online Calculate, obtain the system load flow data after forecast failure, enter step 2)；

In Fig. 1, step 2 describes, and 2) Voltage security margin is less than setting threshold value η_v,thr1With electricity after forecast failure Pressure margin of safety is less than and sets threshold value η_v,thr2Bus be added separately to out-of-limit bus collection, bus collection to be controlled, by overload Nargin is less than and sets threshold value η_p,thr1It is less than with overload nargin after forecast failure and set threshold value η_p,thr2Equipment respectively It is added to overload prevention device collection, apparatus to be controlled collection, wherein out-of-limit bus collection and bus collection to be controlled are referred to as bus and investigate collection, and overload sets Standby collection and apparatus to be controlled collection are referred to as equipment and investigate collection；

If out-of-limit bus collection and overload prevention device collection are sky, there is not apparatus overload under electrical network current state or voltage is got over Limit, then terminate to settle accounts；Otherwise, start the emergency state control aid decision based on dynamic partition to calculate, enter step 3)；

Described equipment refers to circuit and transformer.

Described dynamic partition refers to, when bulk power grid enters the state of emergency, by carrying out switching to part circuit, change Topological structure of electric, to be mitigated or eliminated the technology of the state of emergency.Dynamic partitioning techniques specifically include two kinds, and one is that load turns For two is the solution cyclization of electromagnetic looped network.

The busbar voltage margin of safety of described electrical network, apparatus overload margin of safety computational methods as follows：

The state of emergency that the present invention considers includes two aspects：One is voltage out-of-limit under current state, apparatus overload；Two It is voltage out-of-limit and apparatus overload to occur after forecast failure, relate separately to Corrective control and the prevention and control of electrical network.According to formula (1)～(2) calculate the Voltage security margin of bus.

In formula：η_v,iVoltage security margin for bus i；U_iVoltage magnitude for bus i；WithIt is respectively its upper limit And lower limit.

Calculate the overload nargin of circuit and transformer according to formula (3) and (4) respectively.

In formula：η_pl,iThe overload nargin of circuit i；I_liCurrent value for circuit i；I_li.limThermally-stabilised limit for circuit i Value.

In formula：η_pt,iOverload nargin for transformer i；S_tiApparent energy for transformer i；S_ti.limFor transformer i Rated capacity.

In Fig. 1, step 3 describes, set the dynamic partition scheme switching number of lines, screening can switching circuit, formation can Select solution space, enter step 4)；

The formulation in the described alternative space being formed by switching circuit need to follow following principle：

(1) circuit of switching is 220kV circuit；

(2) circuit putting into is to be in the circuit stopped transport and overhaul done state；

(3) folding scheme need to ensure power supply reliability it is ensured that isolated network will not be formed；

(4) consider that only relying on dynamic partition scheme possibly cannot be fully solved the electrical network state of emergency, and as switching simultaneously Circuit number excessive, not only can increase amount of calculation, and actual schedule is difficult to operate it is therefore desirable to actual according to electrical network in running Situation arranges the upper limit of switching circuit；

Typically initial switching circuit number is set to 1, incremented by successively on this basis.

In Fig. 1, step 4 describes, if out-of-limit bus collection and the equal non-NULL of overload prevention device collection, that is, under current state, electrical network is same When there is apparatus overload and voltage out-of-limit, then be directed to alternative space in each dynamic partition scheme, calculate it comprehensive Can index R_e,i, integrated performance index is more than and sets threshold value R_e,i,threScheme add to treat check scheme concentrate, enter step 5)；If out-of-limit bus collection is sky, overload prevention device collection non-NULL, that is, under current state, electrical network there occurs apparatus overload, then be directed to optional Each dynamic partition scheme in solution space, calculates its integrated performance index ω_p,i, integrated performance index is more than and sets door Threshold value ω_p,i,threScheme add to treat check scheme concentrate, enter step 5)；If out-of-limit bus collection non-NULL, overload prevention device collection is Sky, that is, under current state, electrical network there occurs that busbar voltage is out-of-limit, then be directed to each the dynamic partition scheme in alternative space, Calculate its integrated performance index ω_u,i, integrated performance index is more than and sets threshold value ω_u,i,threScheme add to and treat check side Case is concentrated, and enters step 5)；

Described integrated performance index (R_e、ω_u,iAnd ω_p,i) characterize i-th kind of dynamic partition scheme to current state or anticipation The improvement degree of the specific state of emergency being occurred under fault, its value is bigger, and effect is better.

There is the multiclass state of emergency in electrical network, when that is, voltage out-of-limit and apparatus overload simultaneously in electrical network, counted according to formula (5) Calculate a set of dynamic partition scheme to the effect improving the multiclass state of emergency.

In formula：R_e,iIntegrated performance index for voltage out-of-limit and i-th kind of dynamic partition scheme of apparatus overload；I=1, 2 ..., I, I investigate the number of devices concentrated, j=1 for equipment, and 2 ..., J, J investigate the bus number concentrated, η ' for bus_p,iWith η '_v,j It is respectively the apparatus overload margin of safety after dynamic partition scheme is implemented and busbar voltage margin of safety, k_pvFor overload nargin Conversion coefficient and Voltage security margin between, shows that the attention degrees different with voltage out-of-limit to apparatus overload (typically take 0.5 Between～1.5).The η when relevant device is for circuit_p,iEqual to η_pl,i, the η when relevant device is for transformer_p,iEqual to η_pt,i.

There is the unitary class state of emergency in electrical network, when that is, voltage out-of-limit or apparatus overload in electrical network, respectively according to formula (6), formula (7) calculates the improvement to the state of emergency for the dynamic partition scheme.

In formula：ω_u,iIt is the integrated performance index of i-th kind of dynamic partition scheme for voltage out-of-limit；V investigates for bus Collection, including out-of-limit bus collection and bus collection to be controlled；γ_vImportance modifying factor for different electric pressure buses；N is current shape The bus number of voltage out-of-limit, λ under state_i,j,vVoltage security margin for wherein j-th bus；S_i,j,vDynamic for implementing i-th kind The Voltage security margin variable quantity of j-th bus after partition scheme；W is the sum of forecast failure, N_kFor under wherein k-th fault There is the bus number of voltage out-of-limit, λ_i,l,k,vVoltage security margin for l-th bus under k-th fault；S_i,l,k,vFor i-th Plant the variable quantity of the Voltage security margin of l-th bus under k-th fault after dynamic partition scheme is implemented.

In formula：ω_p,iIt is the integrated performance index of i-th kind of dynamic partition scheme for apparatus overload；V investigates for equipment Collection, that is, equipment (circuit/transformer) intersection transshipped；γ_vImportance modifying factor for different electric pressure equipment；N is to work as The equipment number of overload, θ under front state_i,j,vOverload nargin for wherein j-th equipment；μ_i,j,vDynamic for implementing i-th kind The variable quantity of the overload nargin of j-th equipment after partition scheme；W is the sum of forecast failure, N_kFor wherein k-th event The lower equipment number that overload occurs of barrier, θ_i,l,k,vOverload nargin for l-th equipment under k-th fault；μ_i,l,k,vFor i-th Plant the variable quantity of the overload nargin of l-th equipment under k-th fault after dynamic partition scheme is implemented.

In Fig. 1, step 5 describes, if treating check scheme collection non-NULL, for treat check scheme concentrate each dynamically divide Area scheme, after being implemented using AC power flow calculating dynamic partition scheme, bus is investigated and is concentrated the Voltage security margin of each bus and set Appendix examines the overload nargin concentrating each equipment, if the overload nargin of the Voltage security margin of each bus and equipment is all higher than setting Set threshold value, then add feasible program to concentrate the program.Enter step 6).

If treating, check scheme collection is sky, judges whether switching circuit number has reached the upper limit, if the not up to upper limit, increases Circuit switching bar number, redefines alternative space, enters step 3)；If switching circuit number has reached the upper limit, terminate Calculate；

In Fig. 1, step 6 describes, if feasible program collection non-NULL, for multiple dynamic point that feasible program is concentrated Area scheme, based on dynamic partition quantitatively evaluating index system, carries out overall merit, selects optimum control scheme, terminates to calculate；If Feasible program collection is sky, then be ranked up each scheme treating check scheme concentration by the size of integrated performance index, take comprehensive Close the maximum dynamic partition scheme of performance indications as last control program, terminate to calculate.

The method that the described dynamic partition scheme to feasible program concentration carries out overall merit is as follows：

Based on dynamic partition quantitatively evaluating index system, for each scheme, it is first determined the weight of its index at different levels, enter And the comprehensive evaluation index of upper level index is obtained by certain mathematical algorithms, successively calculate upwards, finally obtain dynamic The comprehensive evaluation result of state partition scheme.

Dynamic partition quantitatively evaluating index system includes security, reliability and economy three class index, as shown in Fig. 2 Specific targets are as follows：

The first kind, safety indexes：

(1) line load index x_L

In formula：I=1,2 ..., n_l；η_pl,iOverload nargin for circuit i；n_lFor the number of lines；α_liImportant for circuit Property modifying factor.

(2) transformer load index x_T

In formula：I=1,2 ..., n_t；η_pl,iOverload nargin for transformer i；n_tFor transformer number；α_tiFor transformation Device importance modifying factor.

(3) voltage levvl index λ_U

In formula：I=1,2 ..., n_u；n_uFor bus quantity；η_v,iVoltage security margin for bus i；α_uiFor busbar voltage Importance modifying factor.

(4) short circuit current level index x_I

In formula：I=1,2 ..., n_sc；n_scFor short trouble number；I_i,scFor the short-circuit current value under i-th short trouble. I_i,sc.limThe minimum of a value of the rupturing capacity of the breaker being connected for short-circuit busbar.

Equations of The Second Kind, reliability index：

(1) current supply circuit number S_M

In formula：S_iFor being not more than the capacity of the 220kV end transformer station of two supplying channels.

(2) spinning reserve index x_re

In formula：P_SRFor active spinning reserve capacity；P_SR.limMinimum limit value for active spinning reserve capacity.

《National grid safety and stability computing technique specification》Point out that the general actual load by no more than 2%～5% determines Active spinning reserve.Therefore, P_SR.limThe 2% of desirable system peak load.

3rd class, economic index：

(1) system losses index S_loss

S_loss=∑ P_G-ΣP_DFormula (14)

In formula：P_GFor generating capacity；P_DFor load capacity.

(2) operate cost index F_c

In formula：N is the number of times of circuit folding；f_iThe operating cost producing for each folding, cost and the use with switch Life-span is relevant.

The determination of the weight of index is the core of comprehensive evaluation algorithm, typically has subjective weighting method, objective weighted model and comprehensive Close enabling legislation.The present embodiment determines the weight of indexs at different levels using the analytic hierarchy process (AHP) being under the jurisdiction of subjective weighting method.

The index being comprised in view of the present invention carried state zonal quantization assessment indicator system is type index less than normal, that is, Index evaluation value is lower, and Criterion Attribute is better, single evaluation index is being carried out assigning temporary, should be assigned to even more important index Lower weight.

After weight determines, the comprehensive evaluation index of every grade of index can be determined by operation methods such as fuzzy mathematicses.This reality Apply the method acquisition comprehensive evaluation result that example is considered as linear weighted function, can be designated as：

In formula：Weight for each index set；β_jWeight for single index in each index set；x_jFor the list in index set Individual index；N is each index set middle finger target number.

In sum, the big I according to comprehensive evaluation result carries out trap queuing to multiple dynamic partition schemes, comprehensive Evaluation result its feasibility less is better.

Although the present invention is open as above with preferred embodiment, embodiment is not for limiting the present invention's.Not Depart from the spirit and scope of the present invention, any equivalence changes done or retouching, also belong to the protection domain of the present invention.Cause The content that this protection scope of the present invention should be defined with claims hereof is as standard.

Claims (5)

1. a kind of bulk power grid emergency state control aid decision-making method based on dynamic partition is it is characterised in that include following walking Suddenly：

1) according to power system operation mode and corresponding model and parameter, calculated online, obtain forecast failure after be System flow data, enters step 2)；

2) Voltage security margin is less than and sets threshold value η_v,thr1It is less than with Voltage security margin after forecast failure and set threshold value η_v,thr2Bus be added separately to out-of-limit bus collection and bus collection to be controlled, by overload nargin be less than set threshold value η_p,thr1 It is less than with overload nargin after forecast failure and set threshold value η_p,thr2Equipment be added separately to overload prevention device collection and to be controlled set Standby collection, wherein out-of-limit bus collection and bus collection to be controlled are referred to as bus and investigate collection, and overload prevention device collection and apparatus to be controlled collection are referred to as Equipment investigates collection；

If out-of-limit bus collection and overload prevention device collection are sky, there is not apparatus overload or voltage out-of-limit under electrical network current state, Then method ends；Otherwise, enter step 3) the emergency state control aid decision calculating based on dynamic partition for the startup；

Described equipment refers to circuit and transformer；

Described dynamic partition refers to, when bulk power grid enters the state of emergency, by part circuit is carried out with switching, change electrical network Topological structure, so that the state of emergency is mitigated or eliminated, including the solution cyclization of load transfer and electromagnetic looped network；

3) set the dynamic partition scheme switching number of lines, screening can switching circuit, formed alternative space, enter step 4)；

4) if out-of-limit bus collection and the equal non-NULL of overload prevention device collection, show under current state that electrical network occurs the multiclass state of emergency and simultaneously There is apparatus overload and voltage out-of-limit, be then directed to each the dynamic partition scheme in alternative space, calculate its combination property Index R_e,i, integrated performance index is more than and sets threshold value R_e,i,threScheme add to treat check scheme concentrate, enter step 5)；

If out-of-limit bus collection is empty and overload prevention device collection non-NULL, shows that under current state, electrical network occurs the unitary class state of emergency, deposit In apparatus overload, then it is directed to each the dynamic partition scheme in alternative space, calculates its integrated performance index ω_p,i, will be comprehensive Close performance indications and be more than setting threshold value ω_p,i,threScheme add to treat check scheme concentrate, enter step 5)；

If out-of-limit bus collection non-NULL and overload prevention device collection are sky, that is, under current state there is the unitary class state of emergency in electrical network, exists Busbar voltage is out-of-limit, then be directed to each the dynamic partition scheme in alternative space, calculate its integrated performance index ω_u,i, will Integrated performance index is more than setting threshold value ω_u,i,threScheme add to treat check scheme concentrate, enter step 5)；

R_e,i、ω_p,iAnd ω_u,i, occurred for i-th kind of dynamic partition scheme of sign is lower to electrical network current state or forecast failure The specific state of emergency improvement degree, its value is bigger, and effect is better；

5) if treating check scheme collection non-NULL, for each the dynamic partition scheme treating that check scheme is concentrated, using AC power flow meter Calculating bus investigation after dynamic partition scheme is implemented concentrates the Voltage security margin of each bus and equipment to investigate the mistake concentrating each equipment Carry margin of safety, if the overload nargin of the Voltage security margin of each bus and equipment is all higher than setting threshold value, by the program Feasible program is added to concentrate；Enter step 6)；

If treating, check scheme collection is sky, judges whether switching circuit number has reached the upper limit, if the not up to upper limit, increases circuit Switching bar number, redefines alternative space, enters step 3)；If switching circuit number has reached the upper limit, terminate we Method；

6) if feasible program collection non-NULL, for multiple dynamic partition schemes of feasible program concentration, quantified based on dynamic partition Assessment indicator system, carries out overall merit, selects optimum control scheme, method ends；If feasible program collection is sky, will treat Each scheme that check scheme is concentrated is ranked up by the size of integrated performance index, takes dynamic point that integrated performance index is maximum Area scheme is as last control program, method ends.

2. the bulk power grid emergency state control aid decision-making method based on dynamic partition according to claim 1, its feature It is, described step 2) in, the Voltage security margin of bus is calculated according to formula (1)～(2)：

In formula：η_v,iVoltage security margin for bus i；U_iVoltage magnitude for bus i；WithBe respectively its upper limit and under Limit；

Described equipment includes circuit and transformer, calculates the overload of circuit and transformer respectively according to formula (3) and (4) Nargin：

In formula：η_pl,iOverload nargin for circuit i；I_liCurrent value for circuit i；I_li.limThermally-stabilised limit for circuit i Value；

In formula：η_pt,iOverload nargin for transformer i；S_tiApparent energy for transformer i；S_ti.limVolume for transformer i Constant volume.

3. the bulk power grid emergency state control aid decision-making method based on dynamic partition according to claim 1, its feature It is, described step 3) in, according to the alternative space of switching circuit formation, following principle need to be followed：Switching for 220kV Circuit；The circuit putting into is to be in the circuit stopped transport and overhaul done state；Folding scheme need to ensure power supply reliability it is ensured that not Isolated network can be formed；Need to arrange the upper limit of switching circuit according to electrical network actual conditions.

4. the bulk power grid emergency state control aid decision-making method based on dynamic partition according to claim 1, its feature It is, described step 4) in, R is calculated according to formula (5)_e,i：

In formula：R_e,iIt is the integrated performance index of i-th kind of dynamic partition scheme for voltage out-of-limit and apparatus overload；I=1, 2 ..., I, I investigate the number of devices concentrated, j=1 for equipment, and 2 ..., J, J investigate the bus number concentrated, η ' for bus_p,iWith η '_v,j It is respectively the apparatus overload margin of safety after dynamic partition scheme is implemented and busbar voltage margin of safety, k_pvFor overload nargin Conversion coefficient and Voltage security margin between, shows the attention degrees different with voltage out-of-limit to apparatus overload；Accordingly setting Standby for η during circuit_p,iEqual to η_pl,i, the η when relevant device is for transformer_p,iEqual to η_pt,i；

ω is calculated according to formula (6)_u,i：

In formula：ω_u,iIt is the integrated performance index of i-th kind of dynamic partition scheme for voltage out-of-limit；V investigates collection, bag for bus Include out-of-limit bus collection and bus collection to be controlled；γ_vImportance modifying factor for different electric pressure buses；N is under current state The bus number of voltage out-of-limit, λ_i,j,vVoltage security margin for wherein j-th bus；S_i,j,vFor implementing i-th kind of dynamic partition The Voltage security margin variable quantity of j-th bus after scheme；W is the sum of forecast failure, N_kFor occurring under wherein k-th fault The bus number of voltage out-of-limit, λ_i,l,k,vVoltage security margin for l-th bus under k-th fault；S_i,l,k,vDynamic for i-th kind The variable quantity of the Voltage security margin of l-th bus under k-th fault after the enforcement of state partition scheme；

ω is calculated according to formula (7)_p,i：

In formula：ω_p,iIt is the integrated performance index of i-th kind of dynamic partition scheme for apparatus overload；V investigates collection for equipment, that is, The equipment intersection of overload；γ_vImportance modifying factor for different electric pressure equipment；N is the equipment of overload under current state Number, θ_i,j,vOverload nargin for wherein j-th equipment；μ_i,j,vFor implementing to set for j-th after i-th kind of dynamic partition scheme The variable quantity of standby overload nargin；W is the sum of forecast failure, N_kFor there is the equipment transshipping under wherein k-th fault Number, θ_i,l,k,vOverload nargin for l-th equipment under k-th fault；μ_i,l,k,vImplement for i-th kind of dynamic partition scheme Afterwards under k-th fault the overload nargin of l-th equipment variable quantity.

5. the bulk power grid emergency state control aid decision-making method based on dynamic partition according to claim 1, its feature It is, described step 6) in, based on dynamic partition quantitatively evaluating index system, the dynamic partition scheme that feasible program is concentrated is entered Row overall merit, is carried out in the following manner：

For each scheme, it is first determined the weight of its index at different levels, and then obtain upper one by certain mathematical algorithms The comprehensive evaluation index of level index, successively calculates upwards, finally obtains the comprehensive evaluation result of dynamic partition scheme；Dynamic partition Quantitatively evaluating index system includes security, reliability and economy three class index, as follows：

The first kind, safety indexes：

Line load index x_L

In formula：I=1,2 ..., n_l；η_pl,iOverload nargin for circuit i；n_lFor the number of lines；α_liRepair for circuit importance Positive divisor；

Transformer load index x_T

In formula：I=1,2 ..., n_t；η_pt,iOverload nargin for transformer i；n_tFor transformer number；α_tiThink highly of for transformation The property wanted modifying factor；

Voltage levvl index λ_U

${\mathrm{\λ}}_U=\frac{1}{n_u}\underset{i=1}{\overset{n_u}{\Σ}}{\mathrm{\α}}_{ui}(1-{\mathrm{\η}}_{v,i})---\left(10\right)$

In formula：I=1,2 ..., n_u；n_uFor bus quantity；η_v,iVoltage security margin for bus i；α_uiFor busbar voltage importance Modifying factor；

Short circuit current level index x_I

$x_I=\frac{1}{n_{sc}}\underset{i=1}{\overset{n_{sc}}{\Σ}}\frac{I_{i,sc}}{I_{i,sc,\lim }}---\left(11\right)$

In formula：I=1,2 ..., n_sc；n_scFor short trouble number；I_i,scFor the short-circuit current value under i-th short trouble； I_i,sc.limThe minimum of a value of the rupturing capacity of the breaker being connected for short-circuit busbar；

Equations of The Second Kind, reliability index：

Current supply circuit number S_M

$S_M=\underset{i\∈R}{\Σ}{S}_i---\left(12\right)$

In formula：S_iFor being not more than the capacity of the 220kV end transformer station of two supplying channels；

Spinning reserve index x_re

$x_{re}=\frac{P_{SR.\lim }}{P_{SR}}---\left(13\right)$

In formula：P_SRFor active spinning reserve capacity；P_SR.limMinimum limit value for active spinning reserve capacity；

3rd class, economic index：

System losses index S_loss

S_loss=∑ P_G-∑P_D(14)

In formula：P_GFor generating capacity；P_DFor load capacity；

Operation cost index F_c

$F_c=\underset{i=1}{\overset{n}{\Σ}}{f}_i---\left(15\right)$

In formula：N is the number of times of circuit folding；f_iThe operating cost producing for each folding, is had with the cost of switch and service life Close.