CN106099907A - Meter and transient state and the online emergent control decision method of static security scleronomic constraint - Google Patents

Abstract

The invention discloses a kind of meter and transient state and the online emergent control decision method of static security scleronomic constraint, belong to power system safety and stability emergent control technical field.The present invention is directed to the feature that bulk power grid multiclass safety and stability problem interweaves and the on-line decision time is short, propose to select different types of control measure according to different safety and stability problems, and respectively all kinds of control measure are screened by Control performance standard, combination is enumerated in measure to controllable devices different after screening again, and overload check based on power sensitivity and steady frequency Security Checking based on power-frequency characteristic screen enumerating combination measure, finally use PC cluster pattern that the combination measure of enumerating after screening is carried out transient state and static security stablizes quantitative evaluation, in conjunction with the calculative strategy that safety and stability Question Classification is preferential, achieve online emergent control quick of meter and transient state and static security scleronomic constraint, Optimal Decision-making.

Description

Meter and transient state and the online emergent control decision method of static security scleronomic constraint

Technical field

The invention belongs to power system safety and stability emergent control technical field, more precisely, the present invention relates to one The method being applicable to the scale grid line emergent control Optimal Decision-making of meter and transient state and static security scleronomic constraint.

Background technology

Bulk power grid running status is changeable, and safety and stability characteristic is complicated, and multiclass safety and stability problem interweaves, safety and stability problem Affect wide, be increasingly difficult to weigh safe operation of electric network stability and economy by calculated off line design control strategy Requirement, even is all difficult to design being able to ensure that the control plan of safe operation of electric network stability in the condition not considering to control cost Slightly, the specific aim using the pattern of On-line Control to improve control strategy it is badly in need of, it is ensured that safe operation of electric network stability, and effectively Reduce and control cost.

Can the calculating speed of control decision be one of the key index being implemented in line traffic control.For guaranteeing the safety of bulk power grid Stable and high effective operation, it is sometimes desirable at alternating current-direct current, sending receiving end to take control measure, control decision space is big simultaneously.Pacify in multiclass Full stable problem interweaves, control device is various, span of control is wide, control accuracy is high and control cost little under conditions of, existing The calculating speed of emergent control decision making algorithm is not met by the basic demand of On-line Control.

Patent " electric power system transient stability and thermally-stabilised coordination emergency control method " (ZL 200710135088.8) proposes Preferentially carry out the emergent control decision-making of apparatus overload safety, then carry out the calculating plan of the emergent control decision-making of transient safe and stable Slightly, in apparatus overload security control performance indications, transient safe and stable Control performance standard is considered.Steady in the safety considered Determining problem aspect, this patent is not counted and the out-of-limit safety of quiescent voltage and the constraint of frequency out-of-limit safety, it is impossible to adapt to cause big The failure condition of power shortage, such as, the sending end Power System Steady-state frequency out-of-limit safety problem of DC bipolar block initiation and receiving end Power System Steady-state voltage out-of-limit safety problem；In terms of the kind of control measure, do not account for straight-flow system power and promptly modulate and arrange Execute to throw with capacitor/reactor and move back measure；In terms of computation schema, do not take PC cluster pattern.

Patent " the PC cluster method of Transient Security for Power Systems stably online emergency control policy " (patent No.: ZL201310145394.5) calculating of online for transient safe and stable emergency control policy is divided into two stages, first carries out transient state The calculating of angle stability emergency control policy, then on the basis of the Optimal Control Strategy searched, carry out transient voltage and Frequency security stablizes the calculating of emergency control policy；Use the steady frequency Security Checking of electrically-based system power frequency characteristic Screen combination measure of enumerating, reduce the combination measure scale of transient safe and stable assessment, and use PC cluster pattern to add Calculate speed soon.But, this patent does not accounts for apparatus overload and the out-of-limit problem of steady state voltage, the most accurately processes stable state frequency The out-of-limit problem of rate.

Additionally, above-mentioned patent does not the most provide the processing method of the bulk power grid for Asynchronous Interconnection, at Control performance standard The aspect such as screening of design, the design of control measure combination and combination measure also need to strengthen, to further speed up control certainly The calculating speed of plan.

Summary of the invention

The present invention seeks to:, peace changeable for alternating current-direct current series-parallel connection bulk power grid (bulk power grid containing Asynchronous Interconnection) running status Full stability characteristic (quality) is complicated, multiclass safety and stability problem interweaves, the impact of safety and stability problem is wide, is badly in need of taking On-line Control pattern To support the actual demand of safety and stability economical operation of bulk power grid, it is provided that a kind of meter and transient state and static security scleronomic constraint Online emergent control decision method, control device is various, span of control is wide, control accuracy is high and controls little the wanting of cost meeting Ask down, reach On-line Control and control decision is calculated the requirement of speed.

The present invention rationale here is that: relative to static security stable problem, transient safe and stable problem is more complicated, shadow The face of sound more greatly, more urgent, be more difficult to decision-making.Therefore, the present invention proposes the calculative strategy using safety and stability Question Classification preferential, excellent First calculate the control strategy of transient safe and stable, then calculate the control strategy that static security is stable.Identical control measure are to difference The impact of safety and stability problem is different, the controlled quentity controlled variable identical different types of control measure shadow to same class safety and stability problem Ring also different.Therefore, the present invention propose according to the type of safety and stability problem be controlled measure Control performance standard point Class evaluation and the type according to control measure carry out control measure category filter based on Control performance standard.

Specifically, the present invention realizes by the following technical solutions, comprises the following steps:

1) according to the control strategy model of safety stabilization control system (referred to as Safety system), and the pressure of Safety system Board status, definite value and real measured data, generate the controlled measure collection CC corresponding to forecast failure defendd with Safety system, and will Under forecast failure, the online emergent control measure collection OC of meter and transient state and static security scleronomic constraint is set to empty set, by be checked Control measure combination increment collection CA is set to empty set, is set to need not take control measure by control measure combination of sets CB to be checked Set, enter step 2)；

Described controlled measure includes controllable device, control direction and controlled quentity controlled variable；

2) for up-to-date operation of power networks state S obtained by operation of power networks section data integration₀, carry out Safety system Transient safe and stable TSS quantitative evaluation and static security under the forecast failure defendd stablize SSS quantitative evaluation, by TSS and SSS quantitative evaluation result is designated as TSR, if TSS and SSS of this forecast failure meets safety and stability requirement, then using CB as OC, Method ends, otherwise, is designated as S by the operation of power networks state obtained based on SSS quantitative evaluation accordingly₁, enter step 3)；

Described TSS includes transient rotor angle stability TAS, Transient Voltage Stability TVS, the safe TVDS of transient voltage dip, transient state Frequency falls safe TFDS and transient frequency rises safe TFRS；

The quantitative evaluation of described TAS refers to the dominant pattern W by being calculated TAS_tsaAnd nargin η_a, for by directly The asynchronous electrical network that streaming system interconnection is constituted, also includes the dominant pattern W of the TAS of each synchronised grids_tsa.iAnd nargin η_a.i, its In, i=1,2 ..., n, n are the number of synchronised grids, by the nargin minimum of any one TAS dominant pattern in n synchronised grids Synchronised grids is defined as TAS key electrical network, and W_tsa、η_aIt is respectively TAS dominant pattern and the nargin of TAS key electrical network；Described master Waveguide mode includes that oscillation center, electromotor hive off, load hives off and the participation factors of electromotor and the participation factors of load, its In, in neck pre-group, in the participation factors of electromotor and delayed group, the participation factors of load is just, the ginseng of electromotor in delayed group It is negative with the participation factors of load in the factor and neck pre-group；

The quantitative evaluation of described TVS refers to the critical load collection W by being calculated TVS_tvsAnd the nargin of critical load, The critical load of described TVS refers to that its TVS nargin and TVS dominate critical load TVS nargin η_tvsNegative less than setting value of difference Lotus, described TVS dominates the load that critical load refers to that TVS nargin is minimum in all loads；

The quantitative evaluation of described TVDS refers to the key node collection W by being calculated TVDS_tvdAnd key node is abundant Degree, the key node of described TVDS refers to that its TVDS nargin and TVDS dominate key node TVDS nargin η_tvdDifference less than set The node of value, described TVDS dominates the node that key node refers to that TVDS nargin is minimum in all nodes；

The quantitative evaluation of described TFDS refers to the key node by being calculated TFDS and key power generator collection W_tfdAnd close Key node and the nargin of key power generator, the key node of described TFDS or key power generator refer to its TFDS nargin and TFDS master Lead key point TFDS nargin η_tfdDifference less than the node of setting value or electromotor, described TFDS dominates key point and refers to all The node of minimum TFDS nargin or electromotor in node and electromotor；

The quantitative evaluation of described TFRS refers to the key node by being calculated TFRS and key power generator collection W_tfrAnd close Key node and the nargin of key power generator, the key node of described TFRS or key power generator refer to its TFRS nargin and TFRS master Lead key point TFRS nargin η_tfrDifference less than the node of setting value or electromotor, described TFRS dominates key point and refers to all TFRS nargin minimum node or electromotor in node and electromotor；

Described SSS includes the safe OLS of apparatus overload, variation safety and frequency shift (FS) safety, wherein variation peace Entirely being divided into again voltage to get over the safe VLLS of lower limit and the safe VULS of Over High-Limit Voltage, frequency shift (FS) is also classified into safely frequency and gets over lower limit peace Full FLLS and frequency get over the safe FULS of the upper limit；

The quantitative evaluation of described OLS refers to the key equipment collection W by being calculated OLS_olAnd the nargin of key equipment, institute The key equipment stating OLS refers to that its OLS nargin and OLS dominate key equipment OLS nargin η_olDifference less than the equipment of setting value, institute State OLS and dominate the equipment that key equipment refers to that OLS nargin is minimum in all devices；

The quantitative evaluation of described VLLS refers to the key node collection W by being calculated VLLS_vlAnd the nargin of key node, The key node of described VLLS refers to that its VLLS nargin and VLLS dominate key node VLLS nargin η_vlDifference less than setting value Node, described VLLS dominates the node that key node refers to that VLLS nargin is minimum in all nodes；

The quantitative evaluation of described VULS refers to the key node collection W by being calculated VULS_vuAnd the nargin of key node, The key node of described VULS refers to that its VULS nargin and VULS dominate key node VULS nargin η_vuDifference less than setting value Node, described VULS dominates key node and refers to node minimum for VULS in all nodes；

The quantitative evaluation of described FLLS refers to nargin η by being calculated FLLS_fl, for being made up of straight-flow system interconnection Asynchronous electrical network, also include FLLS nargin η of each synchronised grids_fl.i, wherein, i=1,2 ..., n, η_flFor each synchronised grids FLLS nargin in minima；

The quantitative evaluation of described FULS refers to nargin η by being calculated FULS_fu, for being made up of straight-flow system interconnection Asynchronous electrical network, also include FULS nargin η of each synchronised grids_fu.i, wherein, i=1,2 ..., n, η_fuFor each synchronised grids FULS nargin in minima；

Described TSS and SSS all meet safety and stability require refer to TAS, TVS, TVDS, TFDS, TFRS, OLS, VLLS, The nargin of VULS, FLLS and FULS totally 10 class safety and stability is all respectively greater than equal to meeting the nargin that safety and stability requires accordingly Threshold value a_m, wherein, m=1,2,3 ..., 10, m equal to 1 time, a₁Represent the nargin threshold value of TAS, when m is equal to 2, a₂Represent TVS Nargin threshold value, m equal to 3 time, a₃Represent the nargin threshold value of TVDS, when m is equal to 4, a₄Represent the nargin threshold value of TFDS, m During equal to 5, a₅Represent the nargin threshold value of TFRS, when m is equal to 6, a₆Represent the nargin threshold value of OLS, when m is equal to 7, a₇Represent The nargin threshold value of VLLS, when m is equal to 8, a₈Represent the nargin threshold value of VULS, when m is equal to 9, a₉Represent the nargin threshold of FLLS Value, when m is equal to 10, a₁₀Represent the nargin threshold value of FULS；

3) straight-flow system power to be combined is promptly modulated measure collection TD and be set to empty set, by capacitor/electricity to be combined The throwing of anti-device is moved back measure collection TX and is set to empty set, and electromotor excision measure collection TG to be combined is set to empty set, by load to be combined Excision measure collection TL is set to empty set, enters step 4)；

4) if in CB all control measure combination η_aBoth less than a₁, the most therefrom choose η_aMaximum control measure combination, makees For OC, and update TSR, S according to OC₁And CC, enter step 5)；

If the η of all control measure combination in CB_aMore than or equal to a₁Control measure combination in have η_tvsLess than a₂Or η_tvdLittle In a₃Control measure combination, then from these control measure combine choose η_tvsWith η_tvdThe control measure combination that sum is maximum, As OC, and update TSR, S according to OC₁And CC, enter step 6)；

If the η of all control measure combination in CB_aMore than or equal to a₁And η_tvsMore than or equal to a₂And η_tvdMore than or equal to a₃Control Combined measure processed there is η_tfdLess than a₄Or η_tfrLess than a₅Control measure combination, then from these control measure combine choose η_tfd With η_tfrThe control measure combination that sum is maximum, as OC, and updates TSR, S according to OC₁And CC, enter step 6)；

If the η of all control measure combination in CB_aMore than or equal to a₁And η_tvsMore than or equal to a₂And η_tvdMore than or equal to a₃And η_tfdMore than or equal to a₄And η_tfrMore than or equal to a₅Control measure combination in have η_olsLess than a₆Control measure combination, then from these η is chosen in control measure combination_olsMaximum control measure combination, as OC, and updates TSR, S according to OC₁And CC, enter step Rapid 10)；

If the η of all control measure combination in CB_aMore than or equal to a₁And η_tvsMore than or equal to a₂And η_tvdMore than or equal to a₃And η_tfdMore than or equal to a₄And η_tfrMore than or equal to a₅And η_olsMore than or equal to a₆Control measure combination in have η_vlLess than a₇Or η_vuIt is less than a₈Control measure combination, then from these control measure combine choose η_vlWith η_vuThe control measure combination that sum is maximum, as OC, and update TSR, S according to OC₁And CC, enter step 10)；

Otherwise, the η of all control measure combination from CB_aMore than or equal to a₁And η_tvsMore than or equal to a₂And η_tvdIt is more than or equal to a₃And η_tfdMore than or equal to a₄And η_tfrMore than or equal to a₅And η_olsMore than or equal to a₆And η_vlMore than or equal to a₇And η_vuMore than or equal to a₈'s η is chosen in control measure combination_flWith η_fuThe control measure combination that sum is maximum, as OC, and updates TSR, S according to OC₁And CC, Enter step 14)；

5) for CC gives, receiving end exchange node all each straight-flow systems in TAS key electrical network, with W_tsaBefore middle neck Mass-sending motor participation factors is weights, calculates S respectively₁In lower neck pre-group, each electromotor node exchanges between node with sending end Weighted sum Y that electrical distance is reciprocal₁, and weighted sum Y that between node, electrical distance is reciprocal is exchanged with receiving end₂；If | Y₁|/|Y₂ | more than parameter b more than 1 set₁, then power emergency lifting measure relevant to this straight-flow system in CC is joined TD In, and will (| Y₁|-|Y₂|) as the TAS Control performance standard of these measures；If | Y₂|/|Y₁| more than b₁, then by CC with should The power that straight-flow system is relevant promptly returns fall measure and joins in TD, and will (| Y₂|-|Y₁|) as these measures TAS control Performance indications；

For CC gives, receiving end exchange node and wherein have a synchronised grids to be TAS not in same synchronised grids Each straight-flow system of crucial electrical network, synchronizes electricity for boundary by sending end with the TAS dominant pattern oscillation center of its sending end synchronised grids The exchange node division of net is two set, for boundary, receiving end is same with the TAS dominant pattern oscillation center of its receiving end synchronised grids The exchange node of step electrical network is also divided into two set, if the sending end exchange node of straight-flow system and the TAS of sending end synchronised grids Dominant pattern is led pre-group electromotor node belong to same set, and the receiving end exchange node electricity Tong Bu with receiving end of straight-flow system In the TAS dominant pattern of net, delayed mass-sending motor node belongs to same set, then by merit relevant to this straight-flow system in CC Rate emergency lifting measure joins in TD, and with W_tsaMiddle neck pre-group electromotor participation factors is weights, calculates S₁Lower neck pre-group In each electromotor node exchange with the sending end of the straight-flow system at TAS key electrical network node or receiving end exchange node between electricity Weighted sum Y of gas inverse distance₃, will | Y₃| as the TAS Control performance standard of these measures；If the sending end exchange of straight-flow system In the TAS dominant pattern of node and sending end synchronised grids, delayed mass-sending motor node belongs to same set, and straight-flow system Receiving end exchange node and the TAS dominant pattern of receiving end synchronised grids lead pre-group electromotor node to belong to same set, then will Power relevant to this straight-flow system in CC promptly returns fall measure and joins in TD, and with W_tsaMiddle neck pre-group electromotor participate in because of Son is weights, calculates S₁In lower neck pre-group, each electromotor node is handed over the sending end of the straight-flow system being positioned at TAS key electrical network Weighted sum Y that between stream node or receiving end exchange node, electrical distance is reciprocal₄, and will | Y₄| as the TAS controlling of these measures Can index；

By Control performance standard in TD, the ratio of the maximum in all controlling measurement performance indications is less than setting ginseng with TD Number c₁Measure reject；

For CC belongs to W_tsaEach electromotor excision measure of neck pre-group electromotor, by it at W_tsaIn participation factors Exert oneself with unit generated power and control the ratio of cost, as the TAS Control performance standard of each electromotor excision measure, and Corresponding electromotor is added in TG；

By Control performance standard in TG, the ratio of the maximum in all controlling measurement performance indications is less than setting ginseng with TG Number c₂Measure reject；

Described unit generated power controls cost and refers to that the control cost of electromotor excision measure goes out with generated power The ratio of power；

For CC belongs to W_tsaEach load excision measure of delayed group's load, by it at W_tsaIn participation factors and list The ratio of position load real power control cost, as the TAS Control performance standard of each load excision measure, and by corresponding load Excision measure adds in TL；

By Control performance standard in TL, the ratio of the maximum in all controlling measurement performance indications is less than setting ginseng with TL Number c₃Measure reject；

Described specific load real power control cost refers to the ratio that the control cost of load excision measure is meritorious with load；

Enter step 15)；

6) the load excision measure collection TLV to be combined only considering TVS and TVDS is set to empty set, will only consider TFDS and The load excision measure collection TLF to be combined of TFRS is set to empty set, enters step 7)；

7) if the η of OC_tvsLess than a₂Or η_tvdLess than a₃, then step 8 is entered after carrying out following process), otherwise, enter step 8)；

Measure, meter and W is exited for each capacitor input/reactor in CC_tvsThe nargin of middle critical load, according to S₁Lower capacitor input/reactor exits measure to W_tvsThe sensitivity of middle critical load access node voltage, calculates each respectively Individual capacitor input/reactor exits the measure Control performance standard X to TVS_tvs, meter and W_tvdThe nargin of middle key node, root According to S₁Lower capacitor input/reactor exits measure to W_tvdThe sensitivity of middle key node voltage, calculates each electric capacity respectively Device input/reactor exits the measure Control performance standard X to TVDS_tvd；

Measure is exited, respectively by its X for each capacitor input/reactor in CC_tvsWith X_tvdSum is right as it The control performance aggregative indicator of TVS and TVDS, if this control performance aggregative indicator is more than 0, then joins this measure in TX；

By control performance aggregative indicator in TX, the ratio of the maximum in all controlling measurement composites of performance index is little with TX In setup parameter c₄Measure reject；

For each load excision measure, meter and W in CC_tvsThe nargin of middle critical load and the unit of load excision measure Load real power control cost, according to S₁The excision measure of lower load is to W_tvsThe sensitivity of middle critical load access node voltage, respectively Calculate each load excision measure Control performance standard L to TVS_tvs；Meter and W_tvdNargin and the load of middle key node are cut Except the specific load real power control cost of measure, according to S₁The excision measure of lower load is to W_tvdThe sensitivity of middle key node voltage, Calculate each load excision measure Control performance standard L to TVDS respectively_tvd；

For each load excision measure in CC, respectively by its L_tvsWith L_tvdSum, as its control to TVS and TVDS Composite of performance index, if this control performance aggregative indicator is more than 0, then joins this measure in TLV；

By control performance aggregative indicator in TLV with TLV the ratio of the maximum in all controlling measurement composites of performance index Less than setup parameter c₅Measure reject；

8) if the η of OC_tfdLess than a₄Or η_tfrLess than a₅, then step 9 is entered after carrying out following process), otherwise, enter step 9)；

For each straight-flow system in CC, meter and W_tfdThe nargin of middle key node/electromotor, calculates S respectively₁Under W_tfdMiddle key node/electromotor node exchanges weighted sum Y that between node, electrical distance is reciprocal with sending end_dfd1, and W_tfdIn Key node/electromotor node exchanges weighted sum Y that between node, electrical distance is reciprocal with receiving end_dfd2, meter and W_tfrMiddle crucial joint The nargin of point/electromotor, calculates S respectively₁Lower W_tfrMiddle key node/electromotor node exchanges between node electric with sending end Weighted sum Y of inverse distance_dfr1, and W_tfrMiddle key node/electromotor node exchanges electrical distance between node and falls with receiving end Weighted sum Y of number_dfr2；

η as OC_tfdLess than or equal to η_tfrTime, if | Y_dfd1| more than b₁|Y_dfr1|、|Y_dfd1| more than b |₁Y_dfd2| and | Y_dfr2| More than b₁|Y_dfd2|, then power relevant to this straight-flow system in CC is promptly returned fall measure and join in TD, will (| Y_dfd1|-| Y_dfd2|) as the Control performance standard of this measure；If | Y_dfd1|/|Y_dfr1|、|Y_dfd1|/|Y_dfd2| and | Y_dfr2|/|Y_dfd2| the least In 1/b₁, then power emergency lifting measure relevant to this straight-flow system in CC is joined in TD, and will (| Y_dfd2|-|Y_dfd1 |) as the Control performance standard of this measure；

η as OC_tfdMore than η_tfrTime, if | Y_dfr1| more than b₁|Y_dfd1|、|Y_dfr1| more than b₁|Y_dfr2| and | Y_dfd2| more than b₁ |Y_dfr2|, then power emergency lifting measure relevant to this straight-flow system in CC is joined in TD, will | Y_dfr1|-Y_dfr2| as The Control performance standard of this measure；If | Y_dfr1|/Y_dfd1|、|Y_dfr1|/|Y_dfr2| and | Y_dfd2|/|Y_dfr2| both less than 1/b₁, then will Power relevant to this straight-flow system in CC promptly returns fall measure and joins in TD, and will (| Y_dfr2|-|Y_dfr1|) as this measure Control performance standard；

By Control performance standard in TD, the ratio of the maximum in all controlling measurement performance indications is less than setting ginseng with TD Number c₁Measure reject；

η as OC_tfrLess than a₅Time, excise measure, meter and W first against each electromotor in CC_tfrMiddle key node/ The nargin of electromotor, calculates S respectively₁Lower W_tfrThe joint that middle key node/electromotor node is connected with electromotor excision measure Weighted sum Y that between point, electrical distance is reciprocal_gfr, meter and W_tfdThe nargin of middle key node/electromotor, calculates S respectively₁Under W_tfdThe weighted sum that between the node that middle key node/electromotor node is connected with electromotor excision measure, electrical distance is reciprocal Y_gfd；If | Y_gfr| more than b₁|Y_gfd|, then the excision measure of this electromotor is joined in TG, and will (| Y_gfr|-|Y_gfd|) arrange with this Execute unit generated power exert oneself control cost ratio as its Control performance standard；Then, by Control performance standard in TG The ratio of the maximum in all controlling measurement performance indications is less than setup parameter c with TG₁Measure reject；

η as OC_tfdLess than a₄Time, excise measure, meter and W first against each load in CC_tfdMiddle key node/send out The nargin of motor, calculates S respectively₁Lower W_tfdThe node that the excision measure of middle key node/electromotor node and load is connected it Between reciprocal weighted sum Y of electrical distance_lfd, meter and W_tfrThe nargin of middle key node/electromotor, calculates S respectively₁Lower W_tfrIn Weighted sum Y that between the node that key node/electromotor node is connected with load excision measure, electrical distance is reciprocal_lfr；If | Y_lfd| more than b₁|Y_lfr|, then the excision measure of this load is joined in TLF, and will (| Y_lfd|-|Y_lfr|) bear with this measure unit Lotus is gained merit to exert oneself and controls the ratio of cost as its Control performance standard；Then, by Control performance standard in TLF and institute in TLF The ratio having the maximum in controlling measurement performance indications is less than setup parameter c₁Measure reject；

9) using the union of TLV Yu TLF as TL, step 15 is entered)；

10) the load excision measure collection TLO to be combined only considering OLS is set to empty set, will only consider VLLS and VULS Load excision measure collection TLSV to be combined be set to empty set, enter step 11)；

11) if the η of OC_olsLess than a₆, then step 12 is entered after carrying out following process), otherwise, enter step 12)；

Measure is promptly modulated, based on S for each straight-flow system power in CC₁Lower straight-flow system power promptly modulates measure To W_olThe active power sensitivity of middle key equipment, meter and W_olThe nargin of middle key equipment, calculates each straight-flow system respectively Power promptly modulates the measure Control performance standard D to OLS_ols, by urgent for the Control performance standard straight-flow system power more than 0 Modulation measure joins in TD；

By Control performance standard in TD, the ratio of the maximum in all controlling measurement performance indications is less than setting ginseng with TD Number c₆Measure reject；

For each electromotor excision measure in CC, based on S₁The excision measure of lower electromotor is to W_olGaining merit of middle key equipment Power sensitivity, meter and W_olThe unit generated power of the nargin of middle key equipment and electromotor excision measure is exerted oneself control generation Valency, calculates each electromotor excision measure Control performance standard G to OLS respectively_ols, by Control performance standard sending out more than 0 Motor excision measure joins in TG；

For each load excision measure in CC, based on S₁The excision measure of lower load is to W_olThe active power of middle key equipment Sensitivity, meter and W_olThe nargin of middle key equipment and the specific load real power control cost of load excision measure, calculate respectively Each load excision measure Control performance standard L to OLS_ols, the Control performance standard load excision measure more than 0 is added In TLO；

If TG non-NULL or TLO non-NULL, then by Control performance standard in TG and all controlling measurement performances in both TG and TLO The ratio of the maximum in index is less than setup parameter c₇Measure reject, by Control performance standard in TLO and both TG and TLO The ratio of the maximum in all controlling measurement performance indications is less than setup parameter c₇Measure reject；

12) if the η of OC_vlLess than a₇Or η_vuLess than a₈, then step 13 is entered after carrying out following process), otherwise, enter step 13)；

Measure of moving back, meter and W is thrown for each capacitor in CC/reactor_vlThe nargin of middle key node, according to S₁Lower electricity Container/reactor is thrown and is moved back measure to W_vlThe sensitivity of middle key node voltage, calculates each capacitor/reactor throwing respectively and moves back The measure Control performance standard X to VLLS_vl, meter and W_vuThe nargin of middle key node, according to S₁Lower capacitor/reactor is thrown to move back and is arranged Execute W_vuThe sensitivity of middle key node voltage, calculates each capacitor/reactor throwing respectively and moves back the measure control to VULS Performance indications X_vu；

η as OC_vlLess than or equal to η_vuTime, if X_vlMore than b₂|X_vu|, then this measure is joined in TX, and by (X_vl-| X_vu|) as the Control performance standard of this measure, wherein b₂For set more than 1 parameter；

η as OC_vlMore than η_vuTime, if X_vuMore than b₂|X_vl|, then this measure is joined in TX, and by (X_vu-|X_vl|) Control performance standard as this measure；

By Control performance standard in TX, the ratio of the maximum in all controlling measurement performance indications is less than setting ginseng with TX Number c₈Measure reject；

η as OC_vlLess than a₇Time, first against each load excision measure, meter and W in CC_vlThe nargin of middle key node With load excision measure specific load real power control cost, according to S₁The excision measure of lower load is to W_vlMiddle key node voltage Sensitivity, calculates each load excision measure Control performance standard L to VLLS_vl, meter and W_vuThe nargin of middle key node and The specific load real power control cost of load excision measure, according to S₁The excision measure of lower load is to W_vuThe spirit of middle key node voltage Sensitivity, calculates each load excision measure Control performance standard L to VULS_vu；If L_vlMore than b₂|L_vu|, then by this measure Join in TLSV, and by (L_vl-|L_vu|) as the Control performance standard of this measure；Then, by Control performance standard in TLSV The ratio of the maximum in all controlling measurement performance indications is less than setup parameter c with TLSV₉Measure reject；

13) using the union of TLO Yu TLSV as TL, step 15 is entered)；

14) by sending in CC, receiving end exchange node do not modulates in each straight-flow system related emergency of same synchronised grids Measure joins in TD, joins in TG by all electromotors excision measure in CC, all loads excision measure in CC is added In TL, enter step 15)；

15) if TD non-NULL or TX non-NULL or TG non-NULL or TL non-NULL, then first against emergent control in TD, TX, TG and TL Measure carries out enumerating combination, obtains CA, and rejects and wherein have two and controllable device phase corresponding to above emergent control measure Same combination measure；Then, by apparatus overload Security Checking and steady frequency Security Checking, sieve CA combines measure Choosing, enters step 16), otherwise it is assumed that search is less than OC, method ends；

16) if CA non-NULL, then by the order that control cost is ascending, combination measure in CA is ranked up, wherein, right In controlling the identical combination measure of cost, preferential by the increment of each measure real power control amount relative to OC in combination measure it With ascending sequence, then to the equal combination measure of wherein real power control amount sum by idle relative to OC of its each measure The ascending sequence of increment sum of controlled quentity controlled variable, enters step 17), otherwise it is assumed that search is less than OC, method ends；

Wherein, straight-flow system power is promptly modulated measure and capacitor/reactor and is thrown the control cost of measure of moving back and be set to 0；

17) respectively the union of each control measure combination in CA with OC is combined as control measure to be checked, Generate CB, and the sequence wherein control measure combination sequence in CA combined as control measure to be checked Number, for S₀, under the forecast failure defendd by Safety system, consider that the control measure combination that in CB, each is to be checked is real respectively TSS quantitative evaluation after executing and SSS quantitative evaluation are as a calculating task, and by the sequence of control measure combination to be checked Corresponding calculating task is ranked up, forms scheduling queue, submit to group system and carry out parallel computation；

During parallel computation, if the control measure combination to be checked controlling cost little is able to ensure that this forecast failure TSS and SSS all meet safety and stability requirement, then terminate all control costs and arrange more than the control to be checked of this control cost Execute the calculating task of combination；

Treat that all calculating tasks complete, enter step 18)；

18) if the control measure combination needing in CB to be checked is able to ensure that TSS and SSS of this forecast failure meets safety Stability requirement, then using the control measure combination to be checked that wherein controls Least-cost as final OC, method ends, no Then, step 3 is entered).

Further, step 2) described in the nargin of 10 class safety and stability all standardize in [-1,1] scope, and nargin is 0, Representing that criticality safety is stable, nargin is more than 0, represents safety and stability, and nargin is less than 0, represents and loses safety and stability, and nargin is the biggest, Represent that safety and stability degree is the highest.

Further, step 4) described in TSR, S₁Renewal with CC specifically includes:

Using with update after corresponding for OC TSS and SSS quantitative evaluation result as TSR, the base corresponding for OC after updating The operation of power networks state obtained in SSS quantitative evaluation is as S₁, from CC, reject and control direction identical with controlling equipment in OC Identical and controlled quentity controlled variable less than or equal to the controlled measure of controlled quentity controlled variable in OC, and with OC controls equipment is identical but control direction is contrary Controlled measure.

Further, step 7) in by formula (1)-(4) calculate X_tvs、X_tvd、L_tvsAnd L_tvd；

$X_{tvs}=\frac{1}{J_1}\underset{j_1=1}{\overset{J_1}{\Σ}}\[{(1-{\mathrm{\η}}_{tvs.j_1})}^k{s}_{Qvs.j_1}\]---\left(1\right)$

$X_{tvd}=\frac{1}{J_2}\underset{j_2=1}{\overset{J_2}{\Σ}}\[{(1-{\mathrm{\η}}_{tvd.j_2})}^k{s}_{Qvd.j_2}\]---\left(2\right)$

$L_{tvs}=\frac{1}{J_1}\underset{j_1=1}{\overset{J_1}{\Σ}}\[{(1-{\mathrm{\η}}_{tvs.j_1})}^k{s}_{Lvs.j_1}\]/C_l---\left(3\right)$

$L_{tvd}=\frac{1}{J_2}\underset{j_2=1}{\overset{J_2}{\Σ}}\[{(1-{\mathrm{\η}}_{tvd.j_2})}^k{s}_{Lvd.j_2}\]/C_l---\left(4\right)$

Wherein, J₁For W_tvsMiddle critical load number,For W_tvsMiddle jth₁The TVS nargin of individual critical load, Measure is exited to W for capacitor input/reactor_tvsMiddle jth₁The voltage sensibility of individual critical load access node, J₂For W_tvdIn Key node number,For W_tvdMiddle jth₂The TVDS nargin of individual key node,Move back for capacitor input/reactor Go out measure to W_tvdMiddle jth₂The voltage sensibility of individual key node, k is the parameter more than 1 set,Excise for load Measure is to W_tvsMiddle jth₁The voltage sensibility of individual critical load access node,For load excision measure to W_tvdMiddle jth₂ The voltage sensibility of individual key node, C_lSpecific load real power control cost for load excision measure.

Further, step 8) in by formula (5)-(10) calculate Y_dfd、Y_dfr、Y_gfd、Y_gfr、Y_lfdAnd Y_lfr；

$Y_{dfd}=\frac{1}{J_3}\underset{j_3=1}{\overset{J_3}{\Σ}}\[{(1-{\mathrm{\η}}_{tfd.j_3})}^k{y}_{dfd.j_3}\]---\left(5\right)$

$Y_{dfr}=\frac{1}{J_4}\underset{j_4=1}{\overset{J_4}{\Σ}}\[{(1-{\mathrm{\η}}_{tfr.j_4})}^k{y}_{bfr.j_4}\]---\left(6\right)$

$Y_{gfd}=\frac{1}{J_3}\underset{j_3=1}{\overset{J_3}{\Σ}}\[{(1-{\mathrm{\η}}_{tfd.j_3})}^k{y}_{gfd.j_3}\]---\left(7\right)$

$Y_{gfr}=\frac{1}{J_4}\underset{j_4=1}{\overset{J_4}{\Σ}}\[{(1-{\mathrm{\η}}_{tfr.j_4})}^k{y}_{gfr.j_4}\]---\left(8\right)$

$Y_{lfd}=\frac{1}{J_3}\underset{j_3=1}{\overset{J_3}{\Σ}}\[{(1-{\mathrm{\η}}_{tfd.j_3})}^k{y}_{lfd.f_3}\]---\left(9\right)$

$Y_{lfr}=\frac{1}{J_4}\underset{j_4=1}{\overset{J_4}{\Σ}}\[{(1-{\mathrm{\η}}_{tfr.j_4})}^k{y}_{lfr.j_4}\]---\left(10\right)$

Wherein, J₃For W_tfdMiddle key node/electromotor number,For W_tfdMiddle jth₃Individual key node/electromotor TFDS nargin,Node and W is exchanged for straight-flow system sending end node or receiving end_tfdMiddle jth₃Individual key node/electromotor The inverse of electrical distance, J between node₄For W_tfrMiddle key node/electromotor number,For W_tfrMiddle jth₄Individual crucial joint The TFRS nargin of point/electromotor,Node and W is exchanged for straight-flow system sending end node or receiving end_tfrMiddle jth₄Individual key The inverse of electrical distance between node/electromotor node, k is the parameter more than 1 set,Connected by electromotor Node and W_tfdMiddle jth₃The inverse of electrical distance between individual key node/electromotor node,Connected by electromotor The node connect and W_tfrMiddle jth₄The inverse of electrical distance between individual key node/electromotor node,Connected by load The node connect and W_tfdMiddle jth₃The inverse of electrical distance between individual key node/electromotor node,Connected by load The node connect and W_tfrMiddle jth₄The inverse of electrical distance between individual key node/electromotor node.

Further, step 11) in by formula (11)-(13) calculate D_ols、G_olsAnd L_ols；

$D_{ols}=\underset{j_5=1}{\overset{J_5}{\Σ}}\[{(1-{\mathrm{\η}}_{ols.j_5})}^k{s}_{PD.j_5}\]---\left(11\right)$

$G_{ols}=\underset{j_5=1}{\overset{J_5}{\Σ}}\[{(1-{\mathrm{\η}}_{ols.j_5})}^k{s}_{PD.j_5}\]/C_g---\left(12\right)$

$L_{ols}=\underset{j_5=1}{\overset{J_5}{\Σ}}\[{(1-{\mathrm{\η}}_{ols.j_5})}^k{s}_{PL.j_5}\]/C_l---\left(13\right)$

Wherein, J₅For W_olMiddle key equipment number,For W_olMiddle jth₅The OLS nargin of individual key equipment,For Straight-flow system is to W_olMiddle jth₅The active power sensitivity of individual key equipment,For electromotor to W_olMiddle jth₅Individual key sets Standby active power sensitivity,For load to W_olMiddle jth₅The active power sensitivity of individual key equipment, k is set Parameter more than 1, C_gUnit generated power for electromotor excision measure controls cost, C_lUnit for load excision measure Load real power control cost.

Further, step 12) in by formula (14)-(17) calculate X_vl、X_vu、L_vlAnd L_vu；

$X_{vl}=\frac{1}{J_6}\underset{j_6=1}{\overset{J_6}{\Σ}}\[{(1-{\mathrm{\η}}_{vl.j_6})}^k{s}_{Qvl.j_6}\]---\left(14\right)$

$X_{vu}=\frac{1}{J_7}\underset{j_7=1}{\overset{J_7}{\Σ}}\[{(1-{\mathrm{\η}}_{vu,j_7})}^k{s}_{Qvu.j_7}\]---\left(15\right)$

$L_{vl}=\frac{1}{J_6}\underset{j_6=1}{\overset{J_6}{\Σ}}\[{(1-{\mathrm{\η}}_{vl.j_6})}^k{s}_{Lvl.j_6}\]/C_l---\left(16\right)$

$L_{vu}=\frac{1}{J_7}\underset{j_7=1}{\overset{J_7}{\Σ}}\[{(1-{\mathrm{\η}}_{vu.j_7})}^k{s}_{Lvu.j_7}\]/C_l---\left(17\right)$

Wherein, J₆For W_vlMiddle key node number,For W_vlMiddle jth₆The VLLS nargin of individual key node,For Capacitor/reactor is thrown and is moved back measure to W_vlMiddle jth₆The voltage sensibility of individual key node, J₇For W_vuMiddle key node number,For W_vuMiddle jth₇The VULS nargin of individual key node,Throw for capacitor/reactor and move back measure to W_vuMiddle jth₇ The voltage sensibility of individual key node, k is the parameter more than 1 set,For load excision measure to W_vlMiddle jth₆Individual The voltage sensibility of key node,For load excision measure to W_vuMiddle jth₇The voltage sensibility of individual key node, C_l Specific load real power control cost for load excision measure.

Further, step 15) described in by apparatus overload Security Checking, CA is carried out screening and specifically includes:

Each combination measure being respectively directed in CA, based on S₀, use Sensitivity Analysis Method, calculate group under forecast failure The union of conjunction measure and OC implement after electrical network OLS nargin η '_olIf, η '_olMore than or equal to a₆Or η '_olLess than a₆And η '_olIt is more than In η_ol, then retain this combination measure, otherwise, this combination measure rejected from CA.

Further, step 15) described in by steady frequency Security Checking, CA is carried out screening and specifically includes:

Each combination measure being respectively directed in CA, based on S₀, use the frequence estimation of electrically-based system merit frequency characteristic Method, calculate combine under forecast failure FLLS nargin η of electrical network after the union of measure and OC is implemented '_flWith FULS nargin η '_fu, And carry out following process:

If η '_flLess than a₉、η'_fuLess than a₁₀And η '_flMore than or equal to η_fl、η'_fuMore than or equal to η_fu, then retain this combination and arrange Execute, otherwise, this combination measure is rejected from CA；

If η '_flLess than a₉、η'_fuMore than or equal to a₁₀And η '_flMore than or equal to η_fl, then retain this combination measure, otherwise, should Combination measure is rejected from CA；

If η '_flMore than or equal to a₉、η'_fuLess than a₁₀And η '_fuMore than or equal to η_fu, then retain this combination measure, otherwise, should Combination measure is rejected from CA；

If η '_flMore than or equal to a₉、η'_fuMore than or equal to a₁₀, then this combination measure is retained, otherwise, by this combination measure from CA Middle rejecting.

By using technique scheme, the present invention achieves following technique effect:

The present invention by relatively more complicated, affect face more greatly, safety and stability problem more urgent, that be more difficult to decision-making preferentially locates Reason, it will usually the decision-making more optimized, overall calculation speed also can be faster；Type according to safety and stability problem is controlled The evaluation of classification of the Control performance standard of measure processed and carry out control based on Control performance standard according to the type of control measure Measure category filter processed, and use apparatus overload Security Checking based on power sensitivity and electrically-based system power frequency The control measure enumerating combination are screened by the out-of-limit Security Checking of steady frequency of characteristic, and transient state and quiet can be greatly lowered The amount of calculation of state safety and stability quantitative evaluation, is effectively improved the speed that emergent control decision-making calculates.Big in electrical network scale, the most steady Determining the condition that problem interweaves, controlled measure type measure many, controlled distribution is wide, required precision is high, the effectiveness of this invention will more For significantly.

Accompanying drawing explanation

Fig. 1 is the flow chart of the step 1-step 4 of the inventive method.

Fig. 2 is the flow chart of the step 5-step 18 of the inventive method.

Detailed description of the invention

With reference to the accompanying drawings the present invention is described in further detail.

Step 1 in Fig. 1: according to the control strategy model of Safety system, and pressing plate state, definite value and the reality of Safety system Measurement information, generates the controlled measure collection CC corresponding to forecast failure defendd with Safety system, and by meter under forecast failure and The online emergent control measure collection OC of transient state and static security scleronomic constraint is set to empty set, control measure combination to be checked is increased Quantity set CA is set to empty set, and control measure combination of sets CB to be checked is set to need not take the set of control measure, enters step Rapid 2；

Described controlled measure includes controllable device, control direction and controlled quentity controlled variable；Such as, for m power can be divided into The same straight-flow system of modulation gear, then be processed as m controlled measure (straight-flow system power promptly modulates measure)；For collection Middle cutting load measure can also require to be processed as the controlled measure (load excision measure) of multiple gear according to decision accuracy；

Step 2 in Fig. 1: for up-to-date operation of power networks state S obtained by operation of power networks section data integration₀, enter Transient safe and stable TSS quantitative evaluation and static security under the forecast failure that row Safety system is defendd are stablized SSS and are quantified to comment Estimate, TSS and SSS quantitative evaluation result is designated as TSR, if TSS and SSS of this forecast failure meets safety and stability requirement, then Using CB as OC, method ends, otherwise, the operation of power networks state obtained based on SSS quantitative evaluation accordingly is designated as S₁, enter Enter step 3；

Described TSS includes transient rotor angle stability TAS, Transient Voltage Stability TVS, the safe TVDS of transient voltage dip, transient state Frequency falls safe TFDS and transient frequency rises safe TFRS；

The quantitative evaluation of described TAS refers to the dominant pattern W by being calculated TAS_tsaAnd nargin η_a, for by directly The asynchronous electrical network that streaming system interconnection is constituted, also includes the dominant pattern W of the TAS of each synchronised grids_tsa.iAnd nargin η_a.i, its In, i=1,2 ..., n, n are the number of synchronised grids, by the nargin minimum of any one TAS dominant pattern in n synchronised grids Synchronised grids is defined as TAS key electrical network, and W_tsa、η_aIt is respectively TAS dominant pattern and the nargin of TAS key electrical network；Described master Waveguide mode includes that oscillation center, electromotor hive off, load hives off and the participation factors of electromotor and the participation factors of load, its In, in neck pre-group, in the participation factors of electromotor and delayed group, the participation factors of load is just, the ginseng of electromotor in delayed group It is negative with the participation factors of load in the factor and neck pre-group；

The quantitative evaluation of described TVS refers to the critical load collection W by being calculated TVS_tvsAnd the nargin of critical load, The critical load of described TVS refers to that its TVS nargin and TVS dominate critical load TVS nargin η_tvsDifference less than setting value (generally Be set to 0.2) load, described TVS dominates the load that critical load refers to that in all loads TVS nargin is minimum；

The quantitative evaluation of described TVDS refers to the key node collection W by being calculated TVDS_tvdAnd key node is abundant Degree, the key node of described TVDS refers to that its TVDS nargin and TVDS dominate key node TVDS nargin η_tvdDifference less than set The node of value (being usually arranged as 0.15), described TVDS dominates the joint that key node refers to that TVDS nargin is minimum in all nodes Point；

The quantitative evaluation of described TFDS refers to the key node by being calculated TFDS and key power generator collection W_tfdAnd close Key node and the nargin of electromotor, the key node of described TFDS or key power generator refer to that its TFDS nargin and TFDS are leading and close Key point TFDS nargin η_tfdDifference less than the node of setting value (being usually arranged as 0.15) or electromotor, the leading key of described TFDS Point refers to node or the electromotor that TFDS nargin is minimum in all nodes and electromotor；

The quantitative evaluation of described TFRS refers to the key node by being calculated TFRS and key power generator collection W_tfrAnd close Key node and the nargin of key power generator, the key node of described TFRS or key power generator refer to its TFRS nargin and TFRS master Lead key point TFRS nargin η_tfrDifference dominate less than the node of setting value (being usually arranged as 0.15) or electromotor, described TFRS Key point refers to node or the electromotor that TFRS nargin is minimum in all nodes and electromotor；

Described SSS includes the safe OLS of apparatus overload, variation safety and frequency shift (FS) safety, wherein variation peace Entirely being divided into again voltage to get over the safe VLLS of lower limit and the safe VULS of Over High-Limit Voltage, frequency shift (FS) is also classified into safely frequency and gets over lower limit peace Full FLLS and frequency get over the safe FULS of the upper limit；

The quantitative evaluation of described OLS refers to the key equipment collection W by being calculated OLS_olAnd the nargin of key equipment, institute The key equipment stating OLS refers to that its OLS nargin and OLS dominate key equipment OLS nargin η_olDifference (generally arrange less than setting value Be 0.2) equipment, described OLS dominates the equipment that key equipment refers to that in all devices OLS nargin is minimum；

The quantitative evaluation of described VLLS refers to the key node collection W by being calculated VLLS_vlAnd the nargin of key node, The key node of described VLLS refers to that its VLLS nargin and VLLS dominate key node VLLS nargin η_vlDifference (logical less than setting value Standing be set to 0.15) node, described VLLS dominates the node that key node refers to that VLLS nargin is minimum in all nodes；

The quantitative evaluation of described VULS refers to the key node collection W by being calculated VULS_vuAnd the nargin of key node, The key node of described VULS refers to that its VULS nargin and VULS dominate key node VULS nargin η_vuDifference (logical less than setting value Standing be set to 0.15) node, described VULS dominates the node that key node refers to that VULS nargin is minimum in all nodes；

The quantitative evaluation of described FLLS refers to nargin η by being calculated FLLS_fl, for being made up of straight-flow system interconnection Asynchronous electrical network, also include FLLS nargin η of each synchronised grids_fl.i, wherein, i=1,2 ..., n, η_flFor each synchronised grids FLLS nargin in minima；

The quantitative evaluation of described FULS refers to nargin η by being calculated FULS_fu, for being made up of straight-flow system interconnection Asynchronous electrical network, also include FULS nargin η of each synchronised grids_fu.i, wherein, i=1,2 ..., n, η_fuFor each synchronised grids FULS nargin in minima；

Described TSS and SSS all meet safety and stability require refer to TAS, TVS, TVDS, TFDS, TFRS, OLS, VLLS, The nargin of VULS, FLLS and FULS totally 10 class safety and stability is all respectively greater than equal to meeting the nargin that safety and stability requires accordingly Threshold value a_m, wherein, m=1,2,3 ..., 10, m equal to 1 time, a₁Represent the nargin threshold value of TAS, when m is equal to 2, a₂Represent TVS Nargin threshold value, m equal to 3 time, a₃Represent the nargin threshold value of TVDS, when m is equal to 4, a₄Represent the nargin threshold value of TFDS, m During equal to 5, a₅Represent the nargin threshold value of TFRS, when m is equal to 6, a₆Represent the nargin threshold value of OLS, when m is equal to 7, a₇Represent The nargin threshold value of VLLS, when m is equal to 8, a₈Represent the nargin threshold value of VULS, when m is equal to 9, a₉Represent the nargin threshold of FLLS Value, when m is equal to 10, a₁₀Represent the nargin threshold value of FULS；

The nargin of described 10 class safety and stability is all standardized in [-1,1] scope, and nargin is 0, represents that criticality safety is steady Fixed, nargin is more than 0, represents safety and stability, and nargin is less than 0, represents and loses safety and stability, and nargin is the biggest, represents safety and stability journey Spend the highest；

Step 3 in Fig. 1: straight-flow system power to be combined is promptly modulated measure collection TD and is set to empty set, by be combined Capacitor/reactor throwing is moved back measure collection TX and is set to empty set, electromotor excision measure collection TG to be combined is set to empty set, will treat group The load excision measure collection TL closed is set to empty set, enters step 4；

Step 4 in Fig. 1: if the η of all control measure combination in CB_aBoth less than a₁, the most therefrom choose η_aMaximum control is arranged Execute combination, as OC, and update TSR, S according to OC₁And CC, enter step 5 in Fig. 2；

If the η of all control measure combination in CB_aMore than or equal to a₁Control measure combination in have η_tvsLess than a₂Or η_tvdLittle In a₃Control measure combination, then from these control measure combine choose η_tvsWith η_tvdThe control measure combination that sum is maximum, As OC, and update TSR, S according to OC₁And CC, enter step 6 in Fig. 2；

If the η of all control measure combination in CB_aMore than or equal to a₁And η_tvsMore than or equal to a₂And η_tvdMore than or equal to a₃Control Combined measure processed there is η_tfdLess than a₄Or η_tfrLess than a₅Control measure combination, then from these control measure combine choose η_tfd With η_tfrThe control measure combination that sum is maximum, as OC, and updates TSR, S according to OC₁And CC, enter step 6 in Fig. 2；

If the η of all control measure combination in CB_aMore than or equal to a₁And η_tvsMore than or equal to a₂And η_tvdMore than or equal to a₃And η_tfdMore than or equal to a₄And η_tfrMore than or equal to a₅Control measure combination in have η_olsLess than a₆Control measure combination, then from these η is chosen in control measure combination_olsMaximum control measure combination, as OC, and updates TSR, S according to OC₁And CC, enter Fig. 2 Middle step 10；

If the η of all control measure combination in CB_aMore than or equal to a₁And η_tvsMore than or equal to a₂And η_tvdMore than or equal to a₃And η_tfdMore than or equal to a₄And η_tfrMore than or equal to a₅And η_olsMore than or equal to a₆Control measure combination in have η_vlLess than a₇Or η_vuIt is less than a₈Control measure combination, then from these control measure combine choose η_vlWith η_vuThe control measure combination that sum is maximum, as OC, and update TSR, S according to OC₁And CC, enter step 10 in Fig. 2；Otherwise, the η of all control measure combination from CB_aIt is more than Equal to a₁And η_tvsMore than or equal to a₂And η_tvdMore than or equal to a₃And η_tfdMore than or equal to a₄And η_tfrMore than or equal to a₅And η_olsIt is more than or equal to a₆And η_vlMore than or equal to a₇And η_vuMore than or equal to a₈Control measure combination in choose η_flWith η_fuThe control measure group that sum is maximum Close, as OC, and update TSR, S according to OC₁And CC, enter step 14 in Fig. 2；

Described TSR, S₁Renewal with CC specifically includes:

Using with update after corresponding for OC TSS and SSS quantitative evaluation result as TSR, the base corresponding for OC after updating The operation of power networks state obtained in SSS quantitative evaluation is as S₁, from CC, reject and control direction identical with controlling equipment in OC Identical and controlled quentity controlled variable less than or equal to the controlled measure of controlled quentity controlled variable in OC, and with OC controls equipment is identical but control direction is contrary Controlled measure；

Step 5 in Fig. 2: for sending in CC, receiving end exchange node all each straight-flow systems in TAS key electrical network, with W_tsaMiddle neck pre-group electromotor participation factors is weights, calculates S respectively₁In lower neck pre-group, each electromotor node is handed over sending end Weighted sum Y that between stream node, electrical distance is reciprocal₁, and weighted sum Y that between node, electrical distance is reciprocal is exchanged with receiving end₂； If | Y₁|/|Y₂| more than parameter b more than 1 set₁(being usually arranged as 1.5), then by merit relevant to this straight-flow system in CC Rate emergency lifting measure joins in TD, and will (| Y₁|-|Y₂|) as the TAS Control performance standard of these measures；If | Y₂|/| Y₁| more than b₁, then power relevant to this straight-flow system in CC is promptly returned fall measure and join in TD, and will (| Y₂|-|Y₁|) TAS Control performance standard as these measures；

For CC gives, receiving end exchange node and wherein have a synchronised grids to be TAS not in same synchronised grids Each straight-flow system of crucial electrical network, synchronizes electricity for boundary by sending end with the TAS dominant pattern oscillation center of its sending end synchronised grids The exchange node division of net is two set, for boundary, receiving end is same with the TAS dominant pattern oscillation center of its receiving end synchronised grids The exchange node of step electrical network is also divided into two set, if the sending end exchange node of straight-flow system and the TAS of sending end synchronised grids Dominant pattern is led pre-group electromotor node belong to same set, and the receiving end exchange node electricity Tong Bu with receiving end of straight-flow system In the TAS dominant pattern of net, delayed mass-sending motor node belongs to same set, then by merit relevant to this straight-flow system in CC Rate emergency lifting measure joins in TD, and with W_tsaMiddle neck pre-group electromotor participation factors is weights, calculates S₁Lower neck pre-group In each electromotor node exchange with the sending end of the straight-flow system at TAS key electrical network node or receiving end exchange node between electricity Weighted sum Y of gas inverse distance₃, will | Y₃| as the TAS Control performance standard of these measures；If the sending end exchange of straight-flow system In the TAS dominant pattern of node and sending end synchronised grids, delayed mass-sending motor node belongs to same set, and straight-flow system Receiving end exchange node and the TAS dominant pattern of receiving end synchronised grids lead pre-group electromotor node to belong to same set, then will Power relevant to this straight-flow system in CC promptly returns fall measure and joins in TD, and with W_tsaMiddle neck pre-group electromotor participate in because of Son is weights, calculates S₁In lower neck pre-group, each electromotor node is handed over the sending end of the straight-flow system being positioned at TAS key electrical network Weighted sum Y that between stream node or receiving end exchange node, electrical distance is reciprocal₄, and will | Y₄| as the TAS controlling of these measures Can index；

By Control performance standard in TD, the ratio of the maximum in all controlling measurement performance indications is less than setting ginseng with TD Number c₁The measure of (being usually arranged as 0.3) is rejected；

For CC belongs to W_tsaEach electromotor excision measure of neck pre-group electromotor, by it at W_tsaIn participation factors Exert oneself with unit generated power and control the ratio of cost, as the TAS Control performance standard of each electromotor excision measure, and Corresponding electromotor is added in TG；

By Control performance standard in TG, the ratio of the maximum in all controlling measurement performance indications is less than setting ginseng with TG Number c₂The measure of (being usually arranged as 0.4) is rejected；

Described unit generated power controls cost and refers to that the control cost of electromotor excision measure goes out with generated power The ratio of power；

For CC belongs to W_tsaEach load excision measure of delayed group's load, by it at W_tsaIn participation factors and list The ratio of position load real power control cost, as the TAS Control performance standard of each load excision measure, and by corresponding load Excision measure adds in TL；

By Control performance standard in TL, the ratio of the maximum in all controlling measurement performance indications is less than setting ginseng with TL Number c₃The measure of (being usually arranged as 0.5) is rejected；

Described specific load real power control cost refers to the ratio that the control cost of load excision measure is meritorious with load；

Enter step 15；

Step 6 in Fig. 2: the load excision measure collection TLV to be combined only considering TVS and TVDS is set to empty set, will only Consider that the load excision measure collection TLF to be combined of TFDS and TFRS is set to empty set, enter step 7；

Step 7 in Fig. 2: if the η of OC_tvsLess than a₂Or η_tvdLess than a₃, then carry out entering after following process step 8, otherwise, Enter step 8；

Measure, meter and W is exited for each capacitor input/reactor in CC_tvsThe nargin of middle critical load, according to S₁Lower capacitor input/reactor exits measure to W_tvsThe sensitivity of middle critical load access node voltage, calculates each respectively Individual capacitor input/reactor exits the measure Control performance standard X to TVS_tvs, meter and W_tvdThe nargin of middle key node, root According to S₁Lower capacitor input/reactor exits measure to W_tvdThe sensitivity of middle key node voltage, calculates each electric capacity respectively Device input/reactor exits the measure Control performance standard X to TVDS_tvd；

Measure is exited, respectively by its X for each capacitor input/reactor in CC_tvsWith X_tvdSum is right as it The control performance aggregative indicator of TVS and TVDS, if this control performance aggregative indicator is more than 0, then joins this measure in TX；

By control performance aggregative indicator in TX, the ratio of the maximum in all controlling measurement composites of performance index is little with TX In setup parameter c₄The measure of (being usually arranged as 0.4) is rejected；

For each load excision measure, meter and W in CC_tvsThe nargin of middle critical load and the unit of load excision measure Load real power control cost, according to S₁The excision measure of lower load is to W_tvsThe sensitivity of middle critical load access node voltage, respectively Calculate each load excision measure Control performance standard L to TVS_tvs；Meter and W_tvdNargin and the load of middle key node are cut Except the specific load real power control cost of measure, according to S₁The excision measure of lower load is to W_tvdThe sensitivity of middle key node voltage, Calculate each load excision measure Control performance standard L to TVDS respectively_tvd；

For each load excision measure in CC, respectively by its L_tvsWith L_tvdSum, as its control to TVS and TVDS Composite of performance index, if this control performance aggregative indicator is more than 0, then joins this measure in TLV；

By control performance aggregative indicator in TLV with TLV the ratio of the maximum in all controlling measurement composites of performance index Less than setup parameter c₅The measure of (being usually arranged as 0.6) is rejected；

X is calculated by formula (1)-(4)_tvs、X_tvd、L_tvsAnd L_tvd；

$X_{tvs}=\frac{1}{J_1}\underset{j_1=1}{\overset{J_1}{\Σ}}\[{(1-{\mathrm{\η}}_{tvs.j_1})}^k{s}_{Qvs.j_1}\]---\left(1\right)$

$X_{tvd}=\frac{1}{J_2}\underset{j_2=1}{\overset{J_2}{\Σ}}\[{(1-{\mathrm{\η}}_{tvd.j_2})}^k{s}_{Qvd.j_2}\]---\left(2\right)$

$L_{tvs}=\frac{1}{J_1}\underset{j_1=1}{\overset{J_1}{\Σ}}\[{(1-{\mathrm{\η}}_{tvs.j_1})}^k{s}_{Lvs.j_1}\]/C_l---\left(3\right)$

$L_{tvd}=\frac{1}{J_2}\underset{j_2=1}{\overset{J_2}{\Σ}}\[{(1-{\mathrm{\η}}_{tvd.j_2})}^k{s}_{Lvd.j_2}\]/C_l---\left(4\right)$

Wherein, J₁For W_tvsMiddle critical load number,For W_tvsMiddle jth₁The TVS nargin of individual critical load, Measure is exited to W for capacitor input/reactor_tvsMiddle jth₁The voltage sensibility of individual critical load access node, J₂For W_tvdIn Key node number,For W_tvdMiddle jth₂The TVDS nargin of individual key node,Move back for capacitor input/reactor Go out measure to W_tvdMiddle jth₂The voltage sensibility of individual key node, k is the parameter more than 1 set,Cut for load Except measure is to W_tvsMiddle jth₁The voltage sensibility of individual critical load access node,For load excision measure to W_tvdIn Jth₂The voltage sensibility of individual key node, C_lSpecific load real power control cost for load excision measure；

Step 8 in Fig. 2: if the η of OC_tfdLess than a₄Or η_tfrLess than a₅, then carry out entering after following process step 9, otherwise, Enter step 9；

For each straight-flow system in CC, meter and W_tfdThe nargin of middle key node/electromotor, calculates S respectively₁Under W_tfdMiddle key node/electromotor node exchanges weighted sum Y that between node, electrical distance is reciprocal with sending end_dfd1, and W_tfdIn Key node/electromotor node exchanges weighted sum Y that between node, electrical distance is reciprocal with receiving end_dfd2, meter and W_tfrMiddle crucial joint The nargin of point/electromotor, calculates S respectively₁Lower W_tfrMiddle key node/electromotor node exchanges between node electric with sending end Weighted sum Y of inverse distance_dfr1, and W_tfrMiddle key node/electromotor node exchanges electrical distance between node and falls with receiving end Weighted sum Y of number_dfr2；

η as OC_tfdLess than or equal to η_tfrTime, if | Y_dfd1| more than b₁|Y_dfr1|、|Y_dfd1| more than b₁|Y_dfd2| and | Y_dfr2| More than b₁|Y_dfd2|, then power relevant to this straight-flow system in CC is promptly returned fall measure and join in TD, will (| Y_dfd1|-| Y_dfd2|) as the Control performance standard of this measure；If | Y_dfd1|/|Y_dfr1|、|Y_dfd1|/|Y_dfd2| and | Y_dfr2|/|Y_dfd2| the least In 1/b₁, then power emergency lifting measure relevant to this straight-flow system in CC is joined in TD, and will (| Y_dfd2|-|Y_dfd1 |) as the Control performance standard of this measure；

η as OC_tfdMore than η_tfrTime, if | Y_dfr1| more than b₁|Y_dfd1|、|Y_dfr1| more than b₁|Y_dfr2| and | Y_dfd2| more than b₁ |Y_dfr2|, then power emergency lifting measure relevant to this straight-flow system in CC is joined in TD, will | Y_dfr1|-|Y_dfr2| make Control performance standard for this measure；If | Y_dfr1|/|Y_dfd1|、|Y_dfr1|/|Y_dfr2| and | Y_dfd2|/|Y_dfr2| both less than 1/b₁, Then power relevant to this straight-flow system in CC is promptly returned fall measure to join in TD, and will (| Y_dfr2|-|Y_dfr1|) as being somebody's turn to do The Control performance standard of measure；

By Control performance standard in TD, the ratio of the maximum in all controlling measurement performance indications is less than setting ginseng with TD Number c₁Measure reject；

η as OC_tfrLess than a₅Time, excise measure, meter and W first against each electromotor in CC_tfrMiddle key node/ The nargin of electromotor, calculates S respectively₁Lower W_tfrThe joint that middle key node/electromotor node is connected with electromotor excision measure Weighted sum Y that between point, electrical distance is reciprocal_gfr, meter and W_tfdThe nargin of middle key node/electromotor, calculates S respectively₁Under W_tfdThe weighted sum that between the node that middle key node/electromotor node is connected with electromotor excision measure, electrical distance is reciprocal Y_gfd；If | Y_gfr| more than b₁|Y_gfd|, then the excision measure of this electromotor is joined in TG, and will (| Y_gfr|-|Y_gfd|) arrange with this Execute unit generated power exert oneself control cost ratio as its Control performance standard；Then, by Control performance standard in TG The ratio of the maximum in all controlling measurement performance indications is less than setup parameter c with TG₁Measure reject；

η as OC_tfdLess than a₄Time, excise measure, meter and W first against each load in CC_tfdMiddle key node/send out The nargin of motor, calculates S respectively₁Lower W_tfdThe node that the excision measure of middle key node/electromotor node and load is connected it Between reciprocal weighted sum Y of electrical distance_lfd, meter and W_tfrThe nargin of middle key node/electromotor, calculates S respectively₁Lower W_tfrIn Weighted sum Y that between the node that key node/electromotor node is connected with load excision measure, electrical distance is reciprocal_lfr；If | Y_lfd| more than b₁|Y_lfr|, then the excision measure of this load is joined in TLF, and will (| Y_lfd|-|Y_lfr|) bear with this measure unit Lotus is gained merit to exert oneself and controls the ratio of cost as its Control performance standard；Then, by Control performance standard in TLF and institute in TLF The ratio having the maximum in controlling measurement performance indications is less than setup parameter c₁Measure reject；

Y is calculated by formula (5)-(10)_dfd、Y_dfr、Y_gfd、Y_gfr、Y_lfdAnd Y_lfr；

$Y_{dfd}=\frac{1}{J_3}\underset{j_3=1}{\overset{J_3}{\Σ}}\[{(1-{\mathrm{\η}}_{tfd.j_3})}^k{y}_{dfd.j_3}\]---\left(5\right)$

$Y_{dfr}=\frac{1}{J_4}\underset{j_4=1}{\overset{J_4}{\Σ}}\[{(1-{\mathrm{\η}}_{tfr.j_4})}^k{y}_{dfr.j_4}\]---\left(6\right)$

$Y_{gfd}=\frac{1}{J_3}\underset{j_3=1}{\overset{J_3}{\Σ}}\[{(1-{\mathrm{\η}}_{tfd.j_3})}^k{y}_{gfd.j_3}\]---\left(7\right)$

$Y_{gfr}=\frac{1}{J_4}\underset{j_4=1}{\overset{J_4}{\Σ}}\[{(1-{\mathrm{\η}}_{tfr.j_4})}^k{y}_{gfr.j_4}\]---\left(8\right)$

$Y_{lfd}=\frac{1}{J_3}\underset{j_3=1}{\overset{J_3}{\Σ}}\[{(1-{\mathrm{\η}}_{tfd.j_3})}^k{y}_{lfd.f_3}\]---\left(9\right)$

$Y_{lfr}=\frac{1}{J_4}\underset{j_4=1}{\overset{J_4}{\Σ}}\[{(1-{\mathrm{\η}}_{tfr.j_4})}^k{y}_{lfr.j_4}\]---\left(10\right)$

Wherein, J₃For W_tfdMiddle key node/electromotor number,For W_tfdMiddle jth₃Individual key node/electromotor TFDS nargin,Node and W is exchanged for straight-flow system sending end node or receiving end_tfdMiddle jth₃Individual key node/electromotor The inverse of electrical distance, J between node₄For W_tfrMiddle key node/electromotor number,For W_tfrMiddle jth₄Individual crucial joint The TFRS nargin of point/electromotor,Node and W is exchanged for straight-flow system sending end node or receiving end_tfrMiddle jth₄Individual key The inverse of electrical distance between node/electromotor node, k is the parameter more than 1 set,Connected by electromotor Node and W_tfdMiddle jth₃The inverse of electrical distance between individual key node/electromotor node,Connected by electromotor The node connect and W_tfrMiddle jth₄The inverse of electrical distance between individual key node/electromotor node,Connected by load The node connect and W_tfdMiddle jth₃The inverse of electrical distance between individual key node/electromotor node,Connected by load The node connect and W_tfrMiddle jth₄The inverse of electrical distance between individual key node/electromotor node.

Step 9 in Fig. 2: using the union of TLV Yu TLF as TL, enters step 15；

Step 10 in Fig. 2: the load excision measure collection TLO to be combined only considering OLS is set to empty set, will only consider The load excision measure collection TLSV to be combined of VLLS and VULS is set to empty set, enters step 11；

If the η of step 11 OC in Fig. 2_olsLess than a₆, then enter step 12 after carrying out following process, otherwise, enter step 12；

Measure is promptly modulated, based on S for each straight-flow system power in CC₁Lower straight-flow system power promptly modulates measure To W_olThe active power sensitivity of middle key equipment, meter and W_olThe nargin of middle key equipment, calculates each straight-flow system respectively Power promptly modulates the measure Control performance standard D to OLS_ols, by urgent for the Control performance standard straight-flow system power more than 0 Modulation measure joins in TD；

By Control performance standard in TD, the ratio of the maximum in all controlling measurement performance indications is less than setting ginseng with TD Number c₆The measure of (being usually arranged as 0.5) is rejected；

For each electromotor excision measure in CC, based on S₁The excision measure of lower electromotor is to W_olGaining merit of middle key equipment Power sensitivity, meter and W_olThe unit generated power of the nargin of middle key equipment and electromotor excision measure is exerted oneself control generation Valency, calculates each electromotor excision measure Control performance standard G to OLS respectively_ols, by Control performance standard sending out more than 0 Motor excision measure joins in TG；

For each load excision measure in CC, based on S₁The excision measure of lower load is to W_olThe active power of middle key equipment Sensitivity, meter and W_olThe nargin of middle key equipment and the specific load real power control cost of load excision measure, calculate respectively Each load excision measure Control performance standard L to OLS_ols, the Control performance standard load excision measure more than 0 is added In TLO；

If TG non-NULL or TLO non-NULL, then by Control performance standard in TG and all controlling measurement performances in both TG and TLO The ratio of the maximum in index is less than setup parameter c₇The measure of (being usually arranged as 0.6) is rejected, by Control performance standard in TLO The ratio of the maximum in all controlling measurement performance indications is less than setup parameter c with both TG and TLO₇Measure reject；

D is calculated by formula (11)-(13)_ols、G_olsAnd L_ols；

$D_{ols}=\underset{j_5=1}{\overset{J_5}{\Σ}}\[{(1-{\mathrm{\η}}_{ols.j_5})}^k{s}_{PD.j_5}\]---\left(11\right)$

$G_{ols}=\underset{j_5=1}{\overset{J_5}{\Σ}}\[{(1-{\mathrm{\η}}_{ols.j_5})}^k{s}_{PD.j_5}\]/C_g---\left(12\right)$

$L_{ols}=\underset{j_5=1}{\overset{J_5}{\Σ}}\[{(1-{\mathrm{\η}}_{ols.j_5})}^k{s}_{PL.j_5}\]/C_l---\left(13\right)$

Wherein, J₅For W_olMiddle key equipment number,For W_olMiddle jth₅The OLS nargin of individual key equipment,For Straight-flow system is to W_olMiddle jth₅The active power sensitivity of individual key equipment,For electromotor to W_olMiddle jth₅Individual key sets Standby active power sensitivity,For load to W_olMiddle jth₅The active power sensitivity of individual key equipment, k is set Parameter more than 1, C_gUnit generated power for electromotor excision measure controls cost, C_lUnit for load excision measure Load real power control cost；

Step 12 in Fig. 2: if the η of OC_vlLess than a₇Or η_vuLess than a₈, then carry out entering after following process step 13, otherwise, Enter step 13；

Measure of moving back, meter and W is thrown for each capacitor in CC/reactor_vlThe nargin of middle key node, according to S₁Lower electricity Container/reactor is thrown and is moved back measure to W_vlThe sensitivity of middle key node voltage, calculates each capacitor/reactor throwing respectively and moves back The measure Control performance standard X to VLLS_vl, meter and W_vuThe nargin of middle key node, according to S₁Lower capacitor/reactor is thrown to move back and is arranged Execute W_vuThe sensitivity of middle key node voltage, calculates each capacitor/reactor throwing respectively and moves back the measure control to VULS Performance indications X_vu；

η as OC_vlLess than or equal to η_vuTime, if X_vlMore than b₂|X_vu|, then this measure is joined in TX, and by (X_vl-| X_vu|) as the Control performance standard of this measure, wherein b₂For set more than 1 parameter (being usually arranged as 1.5)；

η as OC_vlMore than η_vuTime, if X_vuMore than b₂|X_vl|, then this measure is joined in TX, and by (X_vu-|X_vl|) Control performance standard as this measure；

By Control performance standard in TX, the ratio of the maximum in all controlling measurement performance indications is less than setting ginseng with TX Number c₈The measure of (being usually arranged as 0.6) is rejected；

η as OC_vlLess than a₇Time, first against each load excision measure, meter and W in CC_vlThe nargin of middle key node With load excision measure specific load real power control cost, according to S₁The excision measure of lower load is to W_vlMiddle key node voltage Sensitivity, calculates each load excision measure Control performance standard L to VLLS_vl, meter and W_vuThe nargin of middle key node and The specific load real power control cost of load excision measure, according to S₁The excision measure of lower load is to W_vuThe spirit of middle key node voltage Sensitivity, calculates each load excision measure Control performance standard L to VULS_vu；If L_vlMore than b₂|L_vu|, then by this measure Join in TLSV, and by (L_vl-|L_vu|) as the Control performance standard of this measure；Then, by Control performance standard in TLSV The ratio of the maximum in all controlling measurement performance indications is less than setup parameter c with TLSV₉The measure of (being usually arranged as 0.7) Reject；

X is calculated by formula (14)-(17)_vl、X_vu、L_vlAnd L_vu；

$X_{vl}=\frac{1}{J_6}\underset{j_6=1}{\overset{J_6}{\Σ}}\[{(1-{\mathrm{\η}}_{vl.j_6})}^k{s}_{Qvl.j_6}\]---\left(14\right)$

$X_{vu}=\frac{1}{J_7}\underset{j_7=1}{\overset{J_7}{\Σ}}\[{(1-{\mathrm{\η}}_{vu,j_7})}^k{s}_{Qvu.j_7}\]---\left(15\right)$

$L_{vl}=\frac{1}{J_6}\underset{j_6=1}{\overset{J_6}{\Σ}}\[{(1-{\mathrm{\η}}_{vl.j_6})}^k{s}_{Lvl.j_6}\]/C_l---\left(16\right)$

$L_{vu}=\frac{1}{J_7}\underset{j_7=1}{\overset{J_7}{\Σ}}\[{(1-{\mathrm{\η}}_{vu.j_7})}^k{s}_{Lvu.j_7}\]/C_l---\left(17\right)$

Wherein, J₆For W_vlMiddle key node number,For W_vlMiddle jth₆The VLLS nargin of individual key node,For Capacitor/reactor is thrown and is moved back measure to W_vlMiddle jth₆The voltage sensibility of individual key node, J₇For W_vuMiddle key node number,For W_vuMiddle jth₇The VULS nargin of individual key node,Throw for capacitor/reactor and move back measure to W_vuMiddle jth₇ The voltage sensibility of individual key node, k is the parameter more than 1 set,For load excision measure to W_vlMiddle jth₆Individual The voltage sensibility of key node,For load excision measure to W_vuMiddle jth₇The voltage sensibility of individual key node, C_l Specific load real power control cost for load excision measure；

Step 13 in Fig. 2: using the union of TLO Yu TLSV as TL, enters step 15；

Step 14 in Fig. 2: by sending in CC, receiving end exchange node is not correlated with in each straight-flow system of same synchronised grids Urgent modulation measure joins in TD, joins in TG by all electromotors excision measure in CC, is excised by all loads in CC Measure joins in TL, enters step 15；

Step 15 in Fig. 2: if TD non-NULL or TX non-NULL or TG non-NULL or TL non-NULL, then first against in TD, TX, TG and TL Emergent control measure carries out enumerating combination, obtains CA, and reject wherein have two and corresponding to above emergent control measure can The combination measure that control equipment is identical；Then, by apparatus overload Security Checking and steady frequency Security Checking, combination in CA is arranged Execute and screen, enter step 16, otherwise it is assumed that search is less than OC, method ends；

Described by apparatus overload Security Checking, CA carried out screening and specifically includes:

Each combination measure being respectively directed in CA, based on S₀, use Sensitivity Analysis Method, calculate group under forecast failure The union of conjunction measure and OC implement after electrical network OLS nargin η '_olIf, η '_olMore than or equal to a₆Or η '_olLess than a₆And η '_olIt is more than In η_ol, then retain this combination measure, otherwise, this combination measure rejected from CA；

Described by steady frequency Security Checking, CA carried out screening and specifically includes:

Each combination measure being respectively directed in CA, based on S₀, use the frequence estimation of electrically-based system merit frequency characteristic Method, calculate combine under forecast failure FLLS nargin η of electrical network after the union of measure and OC is implemented '_flWith FULS nargin η '_fu, And carry out following process:

If η '_flLess than a₉、η'_fuLess than a₁₀And η '_flMore than or equal to η_fl、η'_fuMore than or equal to η_fu, then retain this combination and arrange Execute, otherwise, this combination measure is rejected from CA；

If η '_flLess than a₉、η'_fuMore than or equal to a₁₀And η '_flMore than or equal to η_fl, then retain this combination measure, otherwise, should Combination measure is rejected from CA；

If η '_flMore than or equal to a₉、η'_fuLess than a₁₀And η '_fuMore than or equal to η_fu, then retain this combination measure, otherwise, should Combination measure is rejected from CA；

If η '_flMore than or equal to a₉、η'_fuMore than or equal to a₁₀, then this combination measure is retained, otherwise, by this combination measure from CA Middle rejecting；

Step 16 in Fig. 2: if CA non-NULL, then by the order that control cost is ascending, combination measure in CA is arranged Sequence, wherein, the combination measure identical for controlling cost, preferential by each measure real power control relative to OC in combination measure The ascending sequence of increment sum of amount, then the combination measure to wherein real power control amount sum is equal is relative by its each measure In the ascending sequence of increment sum of the idle controlled quentity controlled variable of OC, enter step 17, otherwise it is assumed that search is less than OC, terminate this Method；

Wherein, straight-flow system power is promptly modulated measure and capacitor/reactor and is thrown the control cost of measure of moving back and be set to 0；

Step 17 in Fig. 2:

Respectively the union of each control measure combination in CA with OC is combined as control measure to be checked, generate CB, and the sequence number that wherein control measure combination sequence in CA is combined as control measure to be checked, For S₀, under the forecast failure defendd by Safety system, consider that in CB, each is after the control measure combination checked is implemented respectively TSS quantitative evaluation and SSS quantitative evaluation calculate task as one, and by the sequence of control measure combination to be checked to phase The calculating task answered is ranked up, and forms scheduling queue, submits to group system and carry out parallel computation；

During parallel computation, if the control measure combination to be checked controlling cost little is able to ensure that this forecast failure TSS and SSS all meet safety and stability requirement, then terminate all control costs and arrange more than the control to be checked of this control cost Execute the calculating task of combination；

Treat that all calculating tasks complete, enter step 18；

Step 18 in Fig. 2: if the control measure combination needing to be checked in CB is able to ensure that TSS and SSS of this forecast failure All meet safety and stability requirement, then the control measure to be checked wherein controlling Least-cost are combined as final OC, knot Bundle this method, otherwise, enters step 3.

Although the present invention is open as above with preferred embodiment, but embodiment is not for limiting the present invention's.Not Depart from the spirit and scope of the present invention, any equivalence change 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 (9)

1. meter and transient state and the online emergent control decision method of static security scleronomic constraint, it is characterised in that include following step Rapid:

1) according to the control strategy model of Safety system, and pressing plate state, definite value and the real measured data of Safety system, generate with Controlled measure collection CC corresponding to the forecast failure that Safety system is defendd, and by meter and transient state and static security under forecast failure The online emergent control measure collection OC of scleronomic constraint is set to empty set, and control measure combination increment collection CA to be checked is set to sky Collection, is set to need not take the set of control measure by control measure combination of sets CB to be checked, and enters step 2)；

Described controlled measure includes controllable device, control direction and controlled quentity controlled variable；

2) for up-to-date operation of power networks state S obtained by operation of power networks section data integration₀, carry out Safety system and prevented Imperial transient safe and stable TSS quantitative evaluation under forecast failure and static security stablize SSS quantitative evaluation, are measured by TSS and SSS Change assessment result and be designated as TSR, if TSS and SSS of this forecast failure meets safety and stability requirement, then using CB as OC, terminate This method, otherwise, is designated as S by the operation of power networks state obtained based on SSS quantitative evaluation accordingly₁, enter step 3)；

Described TSS includes transient rotor angle stability TAS, Transient Voltage Stability TVS, the safe TVDS of transient voltage dip, transient frequency Fall safe TFDS and transient frequency rises safe TFRS；

The quantitative evaluation of described TAS refers to the dominant pattern W by being calculated TAS_tsaAnd nargin η_a, for by straight-flow system The asynchronous electrical network that interconnection is constituted, also includes the dominant pattern W of the TAS of each synchronised grids_tsa.iAnd nargin η_a.i, wherein, i= 1,2 ..., n, n are the number of synchronised grids, by synchronization electricity minimum for the nargin of any one TAS dominant pattern in n synchronised grids Net is defined as TAS key electrical network, and W_tsa、η_aIt is respectively TAS dominant pattern and the nargin of TAS key electrical network；Described dominant pattern Hive off including oscillation center, electromotor, load hives off and the participation factors of electromotor and the participation factors of load, wherein, and neck In pre-group electromotor participation factors and in delayed group the participation factors of load be just, the participation factors of electromotor in delayed group And the participation factors of load is negative in neck pre-group；

The quantitative evaluation of described TVS refers to the critical load collection W by being calculated TVS_tvsAnd the nargin of critical load, described The critical load of TVS refers to that its TVS nargin and TVS dominate critical load TVS nargin η_tvsDifference less than the load of setting value, institute State TVS and dominate the load that critical load refers to that TVS nargin is minimum in all loads；

The quantitative evaluation of described TVDS refers to the key node collection W by being calculated TVDS_tvdAnd the nargin of key node, institute The key node stating TVDS refers to that its TVDS nargin and TVDS dominate key node TVDS nargin η_tvdDifference less than the joint of setting value Point, described TVDS dominates the node that key node refers to that TVDS nargin is minimum in all nodes；

The quantitative evaluation of described TFDS refers to the key node by being calculated TFDS and key power generator collection W_tfdAnd key joint Point and the nargin of key power generator, the key node of described TFDS or key power generator refer to that its TFDS nargin and TFDS are leading and close Key point TFDS nargin η_tfdDifference less than the node of setting value or electromotor, described TFDS dominates key point and refers at all nodes The node of minimum TFDS nargin or electromotor with in electromotor；

The quantitative evaluation of described TFRS refers to the key node by being calculated TFRS and key power generator collection W_tfrAnd key joint Point and the nargin of key power generator, the key node of described TFRS or key power generator refer to that its TFRS nargin and TFRS are leading and close Key point TFRS nargin η_tfrDifference less than the node of setting value or electromotor, described TFRS dominates key point and refers at all nodes With TFRS nargin minimum node or electromotor in electromotor；

Described SSS includes the safe OLS of apparatus overload, variation safety and frequency shift (FS) safety, and wherein variation safety is again Being divided into voltage to get over the safe VLLS of lower limit and the safe VULS of Over High-Limit Voltage, frequency shift (FS) is also classified into safely frequency and gets over lower limit safety FLLS and frequency get over the safe FULS of the upper limit；

The quantitative evaluation of described OLS refers to the key equipment collection W by being calculated OLS_olAnd the nargin of key equipment, described The key equipment of OLS refers to that its OLS nargin and OLS dominate key equipment OLS nargin η_olDifference less than the equipment of setting value, described OLS dominates the equipment that key equipment refers to that OLS nargin is minimum in all devices；

The quantitative evaluation of described VLLS refers to the key node collection W by being calculated VLLS_vlAnd the nargin of key node, described The key node of VLLS refers to that its VLLS nargin and VLLS dominate key node VLLS nargin η_vlDifference less than the node of setting value, Described VLLS dominates the node that key node refers to that VLLS nargin is minimum in all nodes；

The quantitative evaluation of described VULS refers to the key node collection W by being calculated VULS_vuAnd the nargin of key node, described The key node of VULS refers to that its VULS nargin and VULS dominate key node VULS nargin η_vuDifference less than the node of setting value, Described VULS dominates key node and refers to node minimum for VULS in all nodes；

The quantitative evaluation of described FLLS refers to nargin η by being calculated FLLS_fl, different for be made up of straight-flow system interconnection Step electrical network, also includes FLLS nargin η of each synchronised grids_fl.i, wherein, i=1,2 ..., n, η_flFor each synchronised grids Minima in FLLS nargin；

The quantitative evaluation of described FULS refers to nargin η by being calculated FULS_fu, different for be made up of straight-flow system interconnection Step electrical network, also includes FULS nargin η of each synchronised grids_fu.i, wherein, i=1,2 ..., n, η_fuFor each synchronised grids Minima in FULS nargin；

Described TSS and SSS all meet safety and stability require refer to TAS, TVS, TVDS, TFDS, TFRS, OLS, VLLS, VULS, The nargin of FLLS and FULS totally 10 class safety and stability is all respectively greater than equal to meeting the nargin threshold that safety and stability requires accordingly Value a_m, wherein, m=1,2,3 ..., 10, m equal to 1 time, a₁Represent the nargin threshold value of TAS, when m is equal to 2, a₂Represent the abundant of TVS Degree threshold value, when m is equal to 3, a₃Represent the nargin threshold value of TVDS, when m is equal to 4, a₄Representing the nargin threshold value of TFDS, m is equal to When 5, a₅Represent the nargin threshold value of TFRS, when m is equal to 6, a₆Represent the nargin threshold value of OLS, when m is equal to 7, a₇Represent VLLS Nargin threshold value, m equal to 8 time, a₈Represent the nargin threshold value of VULS, when m is equal to 9, a₉Represent the nargin threshold value of FLLS, m During equal to 10, a₁₀Represent the nargin threshold value of FULS；

3) straight-flow system power to be combined is promptly modulated measure collection TD and be set to empty set, by capacitor/reactor to be combined Throwing is moved back measure collection TX and is set to empty set, electromotor excision measure collection TG to be combined is set to empty set, by load excision to be combined Measure collection TL is set to empty set, enters step 4)；

4) if in CB all control measure combination η_aBoth less than a₁, the most therefrom choose η_aMaximum control measure combination, as OC, and update TSR, S according to OC₁And CC, enter step 5)；

If the η of all control measure combination in CB_aMore than or equal to a₁Control measure combination in have η_tvsLess than a₂Or η_tvdLess than a₃ Control measure combination, then from these control measure combine choose η_tvsWith η_tvdThe control measure combination that sum is maximum, as OC, and update TSR, S according to OC₁And CC, enter step 6)；

If the η of all control measure combination in CB_aMore than or equal to a₁And η_tvsMore than or equal to a₂And η_tvdMore than or equal to a₃Control arrange Execute and combination has η_tfdLess than a₄Or η_tfrLess than a₅Control measure combination, then from these control measure combine choose η_tfdWith η_tfrThe control measure combination that sum is maximum, as OC, and updates TSR, S according to OC₁And CC, enter step 6)；

If the η of all control measure combination in CB_aMore than or equal to a₁And η_tvsMore than or equal to a₂And η_tvdMore than or equal to a₃And η_tfdIt is more than Equal to a₄And η_tfrMore than or equal to a₅Control measure combination in have η_olsLess than a₆Control measure combination, then arrange from these controls Execute and combination is chosen η_olsMaximum control measure combination, as OC, and updates TSR, S according to OC₁And CC, enter step 10)；

If the η of all control measure combination in CB_aMore than or equal to a₁And η_tvsMore than or equal to a₂And η_tvdMore than or equal to a₃And η_tfdIt is more than Equal to a₄And η_tfrMore than or equal to a₅And η_olsMore than or equal to a₆Control measure combination in have η_vlLess than a₇Or η_vuLess than a₈Control Combined measure, then choose η from these control measure combine_vlWith η_vuThe control measure combination that sum is maximum, as OC, and root TSR, S is updated according to OC₁And CC, enter step 10)；Otherwise, the η of all control measure combination from CB_aMore than or equal to a₁And η_tvs More than or equal to a₂And η_tvdMore than or equal to a₃And η_tfdMore than or equal to a₄And η_tfrMore than or equal to a₅And η_olsMore than or equal to a₆And η_vlIt is more than Equal to a₇And η_vuMore than or equal to a₈Control measure combination in choose η_flWith η_fuThe control measure combination that sum is maximum, as OC, And update TSR, S according to OC₁And CC, enter step 14)；

5) for CC gives, receiving end exchange node all each straight-flow systems in TAS key electrical network, with W_tsaMass-send before middle neck Motor participation factors is weights, calculates S respectively₁In lower neck pre-group, each electromotor node exchanges between node electric with sending end Weighted sum Y of inverse distance₁, and weighted sum Y that between node, electrical distance is reciprocal is exchanged with receiving end₂；If | Y₁|/|Y₂| big In parameter b more than 1 set₁, then power emergency lifting measure relevant to this straight-flow system in CC is joined in TD, and Will (| Y₁|-|Y₂|) as the TAS Control performance standard of these measures；If | Y₂|/|Y₁| more than b₁, then by CC with this direct current system The power that system is relevant promptly returns fall measure and joins in TD, and will (| Y₂|-|Y₁|) refer to as the TAS control performance of these measures Mark；

For CC gives, receiving end exchange node and wherein have a synchronised grids to be that TAS is crucial not in same synchronised grids Each straight-flow system of electrical network, with the TAS dominant pattern oscillation center of its sending end synchronised grids for boundary by sending end synchronised grids Exchange node division is two set, for boundary, receiving end is synchronized electricity with the TAS dominant pattern oscillation center of its receiving end synchronised grids The exchange node of net is also divided into two set, if the sending end exchange node of straight-flow system is leading with the TAS of sending end synchronised grids Pattern is led pre-group electromotor node belong to same set, and the receiving end of straight-flow system exchanges node and receiving end synchronised grids In TAS dominant pattern, delayed mass-sending motor node belongs to same set, then by tight for power relevant to this straight-flow system in CC Anxious measure for improvement joins in TD, and with W_tsaMiddle neck pre-group electromotor participation factors is weights, calculates S₁In lower neck pre-group respectively Individual electromotor node exchange with the sending end of the straight-flow system at TAS key electrical network between node or receiving end exchange node electrically away from From weighted sum Y reciprocal₃, will | Y₃| as the TAS Control performance standard of these measures；If the sending end exchange node of straight-flow system Delayed mass-sending motor node belongs to same set, and the receiving end of straight-flow system with the TAS dominant pattern of sending end synchronised grids Exchange node and the TAS dominant pattern of receiving end synchronised grids lead pre-group electromotor node to belong to same set, then by CC The power relevant to this straight-flow system promptly returns fall measure and joins in TD, and with W_tsaMiddle neck pre-group electromotor participation factors is Weights, calculate S₁In lower neck pre-group, each electromotor node exchanges joint with the sending end of the straight-flow system being positioned at TAS key electrical network Weighted sum Y that between point or receiving end exchange node, electrical distance is reciprocal₄, and will | Y₄| the TAS control performance as these measures refers to Mark；

By Control performance standard in TD, the ratio of the maximum in all controlling measurement performance indications is less than setup parameter c with TD₁'s Measure is rejected；

For CC belongs to W_tsaEach electromotor excision measure of neck pre-group electromotor, by it at W_tsaIn participation factors and list Position generated power is exerted oneself and is controlled the ratio of cost, as the TAS Control performance standard of each electromotor excision measure, and by phase The electromotor answered adds in TG；

By Control performance standard in TG, the ratio of the maximum in all controlling measurement performance indications is less than setup parameter c with TG₂'s Measure is rejected；

Described unit generated power controls cost and refers to what control cost and the generated power of electromotor excision measure were exerted oneself Ratio；

For CC belongs to W_tsaEach load excision measure of delayed group's load, by it at W_tsaIn participation factors and unit bear The ratio of lotus real power control cost, as the TAS Control performance standard of each load excision measure, and excises corresponding load Measure adds in TL；

By Control performance standard in TL, the ratio of the maximum in all controlling measurement performance indications is less than setup parameter c with TL₃'s Measure is rejected；

Described specific load real power control cost refers to the ratio that the control cost of load excision measure is meritorious with load；

Enter step 15)；

6) the load excision measure collection TLV to be combined only considering TVS and TVDS is set to empty set, will only consider TFDS and TFRS Load excision measure collection TLF to be combined be set to empty set, enter step 7)；

7) if the η of OC_tvsLess than a₂Or η_tvdLess than a₃, then step 8 is entered after carrying out following process), otherwise, enter step 8)；

Measure, meter and W is exited for each capacitor input/reactor in CC_tvsThe nargin of middle critical load, according to S₁Lower electricity Container input/reactor exits measure to W_tvsThe sensitivity of middle critical load access node voltage, calculates each electric capacity respectively Device input/reactor exits the measure Control performance standard X to TVS_tvs, meter and W_tvdThe nargin of middle key node, according to S₁Under Capacitor input/reactor exits measure to W_tvdThe sensitivity of middle key node voltage, calculates each capacitor respectively and throws Enter/reactor exits the measure Control performance standard X to TVDS_tvd；

Measure is exited, respectively by its X for each capacitor input/reactor in CC_tvsWith X_tvdSum, as its to TVS and The control performance aggregative indicator of TVDS, if this control performance aggregative indicator is more than 0, then joins this measure in TX；

By control performance aggregative indicator in TX with TX the ratio of the maximum in all controlling measurement composites of performance index less than setting Determine parameter c₄Measure reject；

For each load excision measure, meter and W in CC_tvsThe nargin of middle critical load and the specific load of load excision measure Real power control cost, according to S₁The excision measure of lower load is to W_tvsThe sensitivity of middle critical load access node voltage, calculates respectively Go out each load excision measure Control performance standard L to TVS_tvs；Meter and W_tvdThe nargin of middle key node and load excision are arranged The specific load real power control cost executed, according to S₁The excision measure of lower load is to W_tvdThe sensitivity of middle key node voltage, respectively Calculate each load excision measure Control performance standard L to TVDS_tvd；

For each load excision measure in CC, respectively by its L_tvsWith L_tvdSum, as its control performance to TVS and TVDS Aggregative indicator, if this control performance aggregative indicator is more than 0, then joins this measure in TLV；

By control performance aggregative indicator in TLV, the ratio of the maximum in all controlling measurement composites of performance index is less than with TLV Setup parameter c₅Measure reject；

8) if the η of OC_tfdLess than a₄Or η_tfrLess than a₅, then step 9 is entered after carrying out following process), otherwise, enter step 9)；

For each straight-flow system in CC, meter and W_tfdThe nargin of middle key node/electromotor, calculates S respectively₁Lower W_tfdIn Key node/electromotor node exchanges weighted sum Y that between node, electrical distance is reciprocal with sending end_dfd1, and W_tfdMiddle crucial joint Point/electromotor node exchanges weighted sum Y that between node, electrical distance is reciprocal with receiving end_dfd2, meter and W_tfrMiddle key node/send out The nargin of motor, calculates S respectively₁Lower W_tfrMiddle key node/electromotor node exchanges electrical distance between node and falls with sending end Weighted sum Y of number_dfr1, and W_tfrMiddle key node/electromotor node exchanges reciprocal the adding of electrical distance between node with receiving end Power and Y_dfr2；

η as OC_tfdLess than or equal to η_tfrTime, if | Y_dfd1| more than b₁|Y_dfr1|、|Y_dfd1| more than b₁|Y_dfd2| and | Y_dfr2| more than b₁ |Y_dfd2|, then power relevant to this straight-flow system in CC is promptly returned fall measure and join in TD, will (| Y_dfd1|-|Y_dfd2|) Control performance standard as this measure；If | Y_dfd1|/|Y_dfr1|、|Y_dfd1|/|Y_dfd2| and | Y_dfr2|/|Y_dfd2| both less than 1/ b₁, then power emergency lifting measure relevant to this straight-flow system in CC is joined in TD, and will (| Y_dfd2|-|Y_dfd1|) make Control performance standard for this measure；

η as OC_tfdMore than η_tfrTime, if | Y_dfr1| more than b₁|Y_dfd1|、|Y_dfr1| more than b₁|Y_dfr2| and | Y_dfd2| more than b₁| Y_dfr2|, then power emergency lifting measure relevant to this straight-flow system in CC is joined in TD, will | Y_dfr1|-|Y_dfr2| as The Control performance standard of this measure；If | Y_dfr1|/|Y_dfd1|、|Y_dfr1|/|Y_dfr2| and | Y_dfd2|/|Y_dfr2| both less than 1/b₁, then Power relevant to this straight-flow system in CC is promptly returned fall measure join in TD, and will (| Y_dfr2|-|Y_dfr1|) arrange as this The Control performance standard executed；

By Control performance standard in TD, the ratio of the maximum in all controlling measurement performance indications is less than setup parameter c with TD₁'s Measure is rejected；

η as OC_tfrLess than a₅Time, excise measure, meter and W first against each electromotor in CC_tfrMiddle key node/generating The nargin of machine, calculates S respectively₁Lower W_tfrThe node that the excision measure of middle key node/electromotor node and electromotor is connected it Between reciprocal weighted sum Y of electrical distance_gfr, meter and W_tfdThe nargin of middle key node/electromotor, calculates S respectively₁Lower W_tfdIn Weighted sum Y that between the node that key node/electromotor node is connected with electromotor excision measure, electrical distance is reciprocal_gfd；If |Y_gfr| more than b₁|Y_gfd|, then the excision measure of this electromotor is joined in TG, and will (| Y_gfr|-|Y_gfd|) and this measure unit Generated power is exerted oneself and is controlled the ratio of cost as its Control performance standard；Then, by Control performance standard in TG and TG The ratio of the maximum in all controlling measurement performance indications is less than setup parameter c₁Measure reject；

η as OC_tfdLess than a₄Time, excise measure, meter and W first against each load in CC_tfdMiddle key node/electromotor Nargin, calculate S respectively₁Lower W_tfdElectricity between the node that middle key node/electromotor node is connected with load excision measure Weighted sum Y of gas inverse distance_lfd, meter and W_tfrThe nargin of middle key node/electromotor, calculates S respectively₁Lower W_tfrMiddle key Weighted sum Y that between the node that node/electromotor node is connected with load excision measure, electrical distance is reciprocal_lfr；If | Y_lfd| More than b₁|Y_lfr|, then the excision measure of this load is joined in TLF, and will (| Y_lfd|-|Y_lfr|) have with this measure specific load Merit is exerted oneself and is controlled the ratio of cost as its Control performance standard；Then, arrange all with TLF for Control performance standard in TLF The ratio executing the maximum in Control performance standard is less than setup parameter c₁Measure reject；

9) using the union of TLV Yu TLF as TL, step 15 is entered)；

10) the load excision measure collection TLO to be combined only considering OLS is set to empty set, will only consider treating of VLLS and VULS The load excision measure collection TLSV of combination is set to empty set, enters step 11)；

11) if the η of OC_olsLess than a₆, then step 12 is entered after carrying out following process), otherwise, enter step 12)；

Measure is promptly modulated, based on S for each straight-flow system power in CC₁Lower straight-flow system power promptly modulates measure to W_ol The active power sensitivity of middle key equipment, meter and W_olThe nargin of middle key equipment, calculates each straight-flow system power respectively The urgent modulation measure Control performance standard D to OLS_ols, the Control performance standard straight-flow system power more than 0 is promptly modulated Measure joins in TD；

By Control performance standard in TD, the ratio of the maximum in all controlling measurement performance indications is less than setup parameter c with TD₆'s Measure is rejected；

For each electromotor excision measure in CC, based on S₁The excision measure of lower electromotor is to W_olThe active power of middle key equipment Sensitivity, meter and W_olThe unit generated power of the nargin of middle key equipment and electromotor excision measure is exerted oneself control cost, point Do not calculate each electromotor excision measure Control performance standard G to OLS_ols, by the Control performance standard electromotor more than 0 Excision measure joins in TG；

For each load excision measure in CC, based on S₁The excision measure of lower load is to W_olThe active power of middle key equipment is sensitive Degree, meter and W_olThe nargin of middle key equipment and the specific load real power control cost of load excision measure, calculate each respectively The load excision measure Control performance standard L to OLS_ols, the Control performance standard load excision measure more than 0 is joined TLO In；

If TG non-NULL or TLO non-NULL, then by Control performance standard in TG and all controlling measurement performance indications in both TG and TLO In the ratio of maximum less than setup parameter c₇Measure reject, will be all in Control performance standard in TLO and both TG and TLO The ratio of the maximum in controlling measurement performance indications is less than setup parameter c₇Measure reject；

12) if the η of OC_vlLess than a₇Or η_vuLess than a₈, then step 13 is entered after carrying out following process), otherwise, enter step 13)；

Measure of moving back, meter and W is thrown for each capacitor in CC/reactor_vlThe nargin of middle key node, according to S₁Lower electric capacity Device/reactor is thrown and is moved back measure to W_vlThe sensitivity of middle key node voltage, calculates each capacitor/reactor throwing respectively and moves back and arrange Execute the Control performance standard X to VLLS_vl, meter and W_vuThe nargin of middle key node, according to S₁Lower capacitor/reactor is thrown and is moved back measure To W_vuThe sensitivity of middle key node voltage, calculates each capacitor/reactor throwing respectively and moves back the measure controlling to VULS Can index X_vu；

η as OC_vlLess than or equal to η_vuTime, if X_vlMore than b₂|X_vu|, then this measure is joined in TX, and by (X_vl-|X_vu|) make For the Control performance standard of this measure, wherein b₂For set more than 1 parameter；

η as OC_vlMore than η_vuTime, if X_vuMore than b₂|X_vl|, then this measure is joined in TX, and by (X_vu-|X_vl|) as being somebody's turn to do The Control performance standard of measure；

By Control performance standard in TX, the ratio of the maximum in all controlling measurement performance indications is less than setup parameter c with TX₈'s Measure is rejected；

η as OC_vlLess than a₇Time, first against each load excision measure, meter and W in CC_vlThe nargin of middle key node is with negative The specific load real power control cost of lotus excision measure, according to S₁The excision measure of lower load is to W_vlMiddle key node voltage sensitive Degree, calculates each load excision measure Control performance standard L to VLLS_vl, meter and W_vuThe nargin of middle key node and load The specific load real power control cost of excision measure, according to S₁The excision measure of lower load is to W_vuMiddle key node voltage sensitive Degree, calculates each load excision measure Control performance standard L to VULS_vu；If L_vlMore than b₂|L_vu|, then this measure is added Enter in TLSV, and by (L_vl-|L_vu|) as the Control performance standard of this measure；Then, by Control performance standard in TLSV with In TLSV, the ratio of the maximum in all controlling measurement performance indications is less than setup parameter c₉Measure reject；

13) using the union of TLO Yu TLSV as TL, step 15 is entered)；

14) by sending in CC, receiving end exchange node promptly do not modulates at each straight-flow system related power of same synchronised grids Measure joins in TD, joins in TG by all electromotors excision measure in CC, all loads excision measure in CC is added In TL, enter step 15)；

15) if TD non-NULL or TX non-NULL or TG non-NULL or TL non-NULL, then first against emergent control measure in TD, TX, TG and TL Carry out enumerating combination, obtain CA, and reject wherein have two and controllable device corresponding to above emergent control measure identical Combination measure；Then, by apparatus overload Security Checking and steady frequency Security Checking, screen CA combines measure, Enter step 16), otherwise it is assumed that search is less than OC, method ends；

16) if CA non-NULL, then combination measure in CA is ranked up, wherein, for control by the order that control cost is ascending The combination measure that cost processed is identical, preferential by the increment sum of each measure real power control amount relative to OC in combination measure by Little to big sequence, then its each measure idle control relative to OC is pressed in the combination measure to wherein real power control amount sum is equal The ascending sequence of increment sum of amount, enters step 17), otherwise it is assumed that search is less than OC, method ends；

Wherein, straight-flow system power is promptly modulated measure and capacitor/reactor and is thrown the control cost of measure of moving back and be set to 0；

17) respectively the union of each control measure combination in CA with OC is combined as control measure to be checked, generate CB, and the sequence number that wherein control measure combination sequence in CA is combined as control measure to be checked, For S₀, under the forecast failure defendd by Safety system, consider that in CB, each is after the control measure combination checked is implemented respectively TSS quantitative evaluation and SSS quantitative evaluation calculate task as one, and by the sequence of control measure combination to be checked to phase The calculating task answered is ranked up, and forms scheduling queue, submits to group system and carry out parallel computation；

During parallel computation, if the control measure combination to be checked controlling cost little is able to ensure that this forecast failure TSS and SSS meets safety and stability requirement, then terminate all control costs control measure to be checked more than this control cost The calculating task of combination；

Treat that all calculating tasks complete, enter step 18)；

18) if the control measure combination needing in CB to be checked is able to ensure that TSS and SSS of this forecast failure meets safety and stability Requirement, then using the control measure combination to be checked that wherein controls Least-cost as final OC, method ends, otherwise, Enter step 3).

Meter the most according to claim 1 and transient state and the online emergent control decision method of static security scleronomic constraint, its Be characterised by, step 2) described in the nargin of 10 class safety and stability all standardize in [-1,1] scope, and nargin is 0, and expression is faced Boundary's safety and stability, nargin is more than 0, represents safety and stability, and nargin is less than 0, represents and loses safety and stability, and nargin is the biggest, represents peace Full degree of stability is the highest.

Meter the most according to claim 1 and transient state and the online emergent control decision method of static security scleronomic constraint, its Be characterised by, step 4) described in TSR, S₁Renewal with CC specifically includes:

Using with update after corresponding for OC TSS and SSS quantitative evaluation result as TSR, OC after updating corresponding based on SSS The operation of power networks state that quantitative evaluation obtains is as S₁, from CC reject with in OC control equipment identical and control direction is identical and Controlled quentity controlled variable is less than or equal to the controlled measure of controlled quentity controlled variable in OC and identical with controlling equipment in OC but that control direction is contrary is controlled Measure.

Meter the most according to claim 1 and transient state and the online emergent control decision method of static security scleronomic constraint, its Be characterised by, step 7) in by formula (1)-(4) calculate X_tvs、X_tvd、L_tvsAnd L_tvd；

$X_{tvs}=\frac{1}{J_1}\underset{j_1=1}{\overset{J_1}{\mathrm{\Σ}}}\[{(1-{\mathrm{\η}}_{tvs.j_1})}^k{s}_{Qvs.j_1}\]---\left(1\right)$

$X_{tvd}=\frac{1}{J_2}\underset{j_2=1}{\overset{J_2}{\mathrm{\Σ}}}\[{(1-{\mathrm{\η}}_{tvd.j_2})}^k{s}_{Qvd.j_2}\]---\left(2\right)$

$L_{tvs}=\frac{1}{J_1}\underset{j_1=1}{\overset{J_1}{\mathrm{\Σ}}}\[{(1-{\mathrm{\η}}_{tvs,j_1})}^k{s}_{Lvs.j_1}\]/C_l---\left(3\right)$

$L_{tvd}=\frac{1}{J_2}\underset{j_2=1}{\overset{J_2}{\mathrm{\Σ}}}\[{(1-{\mathrm{\η}}_{tvd,j_2})}^k{s}_{Lvd.j_2}\]/C_l---\left(4\right)$

Wherein, J₁For W_tvsMiddle critical load number,For W_tvsMiddle jth₁The TVS nargin of individual critical load,For electric capacity Device input/reactor exits measure to W_tvsMiddle jth₁The voltage sensibility of individual critical load access node, J₂For W_tvdMiddle crucial joint Count,For W_tvdMiddle jth₂The TVDS nargin of individual key node,Measure is exited for capacitor input/reactor To W_tvdMiddle jth₂The voltage sensibility of individual key node, k is the parameter more than 1 set,Measure pair is excised for load W_tvsMiddle jth₁The voltage sensibility of individual critical load access node,For load excision measure to W_tvdMiddle jth₂Individual key The voltage sensibility of node, C_lSpecific load real power control cost for load excision measure.

Meter the most according to claim 1 and transient state and the online emergent control decision method of static security scleronomic constraint, its Be characterised by, step 8) in by formula (5)-(10) calculate Y_dfd、Y_dfr、Y_gfd、Y_gfr、Y_lfdAnd Y_lfr；

$Y_{dfd}=\frac{1}{J_3}\underset{j_3=1}{\overset{J_3}{\mathrm{\Σ}}}\[{(1-{\mathrm{\η}}_{tfd.j_3})}^k{y}_{dfd.j_3}\]---\left(5\right)$

$Y_{dfr}=\frac{1}{J_4}\underset{j_4=1}{\overset{J_4}{\mathrm{\Σ}}}\[{(1-{\mathrm{\η}}_{tfr.j_4})}^k{y}_{dfr.j_4}\]---\left(6\right)$

$Y_{gfd}=\frac{1}{J_3}\underset{j_3=1}{\overset{J_3}{\mathrm{\Σ}}}\[{(1-{\mathrm{\η}}_{tfd.j_3})}^k{y}_{gfd.j_3}\]---\left(7\right)$

$Y_{gfr}=\frac{1}{J_4}\underset{j_4=1}{\overset{J_4}{\mathrm{\Σ}}}\[{(1-{\mathrm{\η}}_{tfr.j_4})}^k{y}_{gfr.j_4}\]---\left(8\right)$

$Y_{lfd}=\frac{1}{J_3}\underset{j_3=1}{\overset{J_3}{\mathrm{\Σ}}}\[{(1-{\mathrm{\η}}_{tfd.j_3})}^k{y}_{lfd.j_3}\]---\left(9\right)$

$Y_{lfr}=\frac{1}{J_4}\underset{j_4=1}{\overset{J_4}{\mathrm{\Σ}}}\[{(1-{\mathrm{\η}}_{tfr.j_4})}^k{y}_{lfr.j_4}\]---\left(10\right)$

Wherein, J₃For W_tfdMiddle key node/electromotor number,For W_tfdMiddle jth₃The TFDS of individual key node/electromotor is abundant Degree,Node and W is exchanged for straight-flow system sending end node or receiving end_tfdMiddle jth₃Individual key node/electromotor node it Between the inverse of electrical distance, J₄For W_tfrMiddle key node/electromotor number,For W_tfrMiddle jth₄Individual key node/generating The TFRS nargin of machine,Node and W is exchanged for straight-flow system sending end node or receiving end_tfrMiddle jth₄Individual key node/send out The inverse of electrical distance between motor node, k is the parameter more than 1 set,The node connected by electromotor with W_tfdMiddle jth₃The inverse of electrical distance between individual key node/electromotor node,The node connected by electromotor With W_tfrMiddle jth₄The inverse of electrical distance between individual key node/electromotor node,The node connected by load With W_tfdMiddle jth₃The inverse of electrical distance between individual key node/electromotor node,The node connected by load With W_tfrMiddle jth₄The inverse of electrical distance between individual key node/electromotor node.

Meter the most according to claim 1 and transient state and the online emergent control decision method of static security scleronomic constraint, its Be characterised by, step 11) in by formula (11)-(13) calculate D_ols、G_olsAnd L_ols；

$D_{ols}=\underset{j_5=1}{\overset{J_5}{\mathrm{\Σ}}}\[{(1-{\mathrm{\η}}_{ols.j_5})}^k{s}_{PD.j_5}\]---\left(11\right)$

$G_{ols}=\underset{j_5=1}{\overset{J_5}{\mathrm{\Σ}}}\[{(1-{\mathrm{\η}}_{ols.j_5})}^k{s}_{PG.j_5}\]/C_g---\left(12\right)$

$L_{ols}=\underset{j_5=1}{\overset{J_5}{\mathrm{\Σ}}}\[{(1-{\mathrm{\η}}_{ols.j_5})}^k{s}_{PL.j_5}\]/C_l---\left(13\right)$

Wherein, J₅For W_olMiddle key equipment number,For W_olMiddle jth₅The OLS nargin of individual key equipment,For direct current System is to W_olMiddle jth₅The active power sensitivity of individual key equipment,For electromotor to W_olMiddle jth₅Individual key equipment Active power sensitivity,For load to W_olMiddle jth₅The active power sensitivity of individual key equipment, k is being more than of setting The parameter of 1, C_gUnit generated power for electromotor excision measure controls cost, C_lSpecific load for load excision measure Real power control cost.

Meter the most according to claim 1 and transient state and the online emergent control decision method of static security scleronomic constraint, its Be characterised by, step 12) in by formula (14)-(17) calculate X_vl、X_vu、L_vlAnd L_vu；

$X_{vl}=\frac{1}{J_6}\underset{j_6=1}{\overset{J_6}{\mathrm{\Σ}}}\[{(1-{\mathrm{\η}}_{vl.j_6})}^k{s}_{Qvl.j_6}\]---\left(14\right)$

$X_{vu}=\frac{1}{J_7}\underset{j_7=1}{\overset{J_7}{\mathrm{\Σ}}}\[{(1-{\mathrm{\η}}_{vu.j_7})}^k{s}_{Qvu.j_7}\]---\left(15\right)$

$L_{vl}=\frac{1}{J_6}\underset{j_6=1}{\overset{J_6}{\mathrm{\Σ}}}\[{(1-{\mathrm{\η}}_{vl.j_6})}^k{s}_{Lvl.j_6}\]/C_l---\left(16\right)$

$L_{vu}=\frac{1}{J_7}\underset{j_7=1}{\overset{J_7}{\mathrm{\Σ}}}\[{(1-{\mathrm{\η}}_{vu.j_7})}^k{s}_{Lvu.j_7}\]/C_l---\left(17\right)$

Wherein, J₆For W_vlMiddle key node number,For W_vlMiddle jth₆The VLLS nargin of individual key node,For electric capacity Device/reactor is thrown and is moved back measure to W_vlMiddle jth₆The voltage sensibility of individual key node, J₇For W_vuMiddle key node number,For W_vuMiddle jth₇The VULS nargin of individual key node,Throw for capacitor/reactor and move back measure to W_vuMiddle jth₇Individual crucial joint The voltage sensibility of point, k is the parameter more than 1 set,For load excision measure to W_vlMiddle jth₆Individual key node Voltage sensibility,For load excision measure to W_vuMiddle jth₇The voltage sensibility of individual key node, C_lCut for load Specific load real power control cost except measure.

Meter the most according to claim 1 and transient state and the online emergent control decision method of static security scleronomic constraint, its Be characterised by, step 15) described in by apparatus overload Security Checking, CA is carried out screening and specifically includes:

Each combination measure being respectively directed in CA, based on S₀, use Sensitivity Analysis Method, calculate combination under forecast failure and arrange Grant OC union implement after electrical network OLS nargin η '_olIf, η '_olMore than or equal to a₆Or η '_olLess than a₆And η '_olIt is more than or equal to η_ol, then retain this combination measure, otherwise, this combination measure rejected from CA.

Meter the most according to claim 1 and transient state and the online emergent control decision method of static security scleronomic constraint, its Be characterised by, step 15) described in by steady frequency Security Checking, CA is carried out screening and specifically includes:

Each combination measure being respectively directed in CA, based on S₀, use the frequency estimation method of electrically-based system merit frequency characteristic, Calculate combine under forecast failure FLLS nargin η of electrical network after the union of measure and OC is implemented '_flWith FULS nargin η '_fu, go forward side by side Row is following to be processed:

If η '_flLess than a₉、η'_fuLess than a₁₀And η '_flMore than or equal to η_fl、η'_fuMore than or equal to η_fu, then retain this combination measure, no Then, this combination measure is rejected from CA；

If η '_flLess than a₉、η'_fuMore than or equal to a₁₀And η '_flMore than or equal to η_fl, then this combination measure is retained, otherwise, by this combination Measure is rejected from CA；

If η '_flMore than or equal to a₉、η'_fuLess than a₁₀And η '_fuMore than or equal to η_fu, then this combination measure is retained, otherwise, by this combination Measure is rejected from CA；

If η '_flMore than or equal to a₉、η'_fuMore than or equal to a₁₀, then retain this combination measure, otherwise, this combination measure picked from CA Remove.