CN103279638B - Large power grid overall situation on-line integrated quantitative evaluation method based on response - Google Patents

Abstract

The invention provides a large power grid overall situation on-line integrated quantitative evaluation method based on a response. The method comprises a first step of acquiring power grid topology structural information from an SCADA system and an EMS system and establishing a corresponding relation between the power grid topology structural information and a WAMS system power grid component, a second step of acquiring present power grid operation method trend data from the SCADA system or the EMS system or a WAMS system or acquiring various preconceived trends or transient state fault time domain data from a DSA system, and a third step of dividing responding data into two operation scenes including a steady state (or a quasi-stable state) and a transient state (or a dynamic state) in a macroscopic mode, wherein static state stabilization situation assessment is carried out on a power grid through an on-line node-facing method, and transient state stabilization situation assessment is carried out on the power grid through an on-line unit-facing method. Static state and transient state comprehensive assessment indicators are constructed based on component class thermostabilization, electric parameter acceptable range and system class stabilization, the comprehensiveness and reasonability of the overall situation assessment indicators are improved, and the efficiency of overall situation integrated assessment is improved through the method that tasks are carried out at the same time.

Description

A kind of full situation of bulk power grid based on response integrated quantitative estimation method online

Technical field

The present invention relates to scale grid line security monitoring with early warning field and in particular to a kind of based on response bulk power grid complete Situation integrated quantitative estimation method online.

Background technology

Guarantee that bulk power grid safety and economical operation are the targets that work about electric power person lays siege to for many years.With Power System Interconnection The continuous expansion of scope, the power industry marketization the deepening continuously of system reform, and UHV transmission, renewable energy power generation Lasting access, increased the uncertainty of operation of power networks so that operation of power networks environment and dynamic behaviour are more complicated, it entered Row stability analysis and the difficulty controlling are also bigger.Modern bulk power grid once has an accident, if processed not in time, consequence will very Seriously.Since only 2000, a lot of large-scale blackout are just there occurs in the world, wherein " 8.14 " beautiful ash moss is in history A big power outage, up to 30,000,000,000 dollars of the daily economic loss of outage, the large-scale blackout frequently occurring is to electricity Net security and stability analysis propose more urgent requirement with control, operation of power networks realtime monitoring.

Currently, it is primarily present following 4 problems in scale grid line security evaluation and early warning field：1. powernet safety Assessment mainly adopts legacy network modeling and simulation form although the power grid security assessment based on modeling and simulating is Electric Power Network Planning With run in requisite important tool, but this kind of method is by factor systems such as electric network model, parameter and numerical computations About, application scale, speed and with the aspects such as true electrical network Conditions Matching degree be difficult to adapt to electrical network monitor in real time requirement it is necessary to Seek new theory and more efficiently solution；2. large-scale blackout is often and disturbs by unpredictalbe cascading failure or at random Move and cause, existing " offline decision-premaking, On-line matching " electrical network peace control pattern based on contingency set is it is impossible to mate electrical network Real working condition, and affected by model and parameter, gained scheme is sometimes overly conservative（Or it is optimistic）, especially fault set multiple shot array Problem seriously limits the quantity that may consider operating mode.Bulk power grid monitor in real time and when preventing large-scale blackout it is possible to aobvious Obtain unable to do what one wishes；3. existing powernet safety estimation system has very strong independence when for particular problem, often Electrical network angle stability and synchronism stability are isolated to come and is processed, or be specifically divided into static problem and transient problems.However, Actually the stabilization of power grids is a kind of unified nonlinear dynamic system stable problem, from the point of view of grid responsive angle, generator rotor angle and electricity Pressure is a kind of external expressive form, and static problem and transient problems have debating of a kind of " disorder of internal organs having quiet, motion in quiescence " again Card thought.So, existing a lot of appraisal procedures and index have very big plyability, need more in run-time environment Evaluation index system that is effective and can truly reflecting grid operation situation；4. the one of the main reasons having a power failure on a large scale puts into operation Powernet monitoring system does not provide grid operation situation quantitative evaluation directly perceived, effective in the slow advanced stage of electric network state Early warning information, also not according to electrical network real working condition change provide corresponding control measure and missed make vigorous efforts to turn the tide optimal when Machine.Therefore, electrical network needs to set up brand-new safety on line appraisal procedure system, realizes bulk power grid increasingly automated, intelligentized Safety on line assessment and early warning.

With the development of computer, communication and network technology, and WAMS (Wide Area Measurement System, extensively Domain measurement system) extensively applied in power network schedule automation, it is that the full situation of operation of power networks based on WAMS quantifies online Assessment and real-time adaptive control research bring new opportunity.In the face of above-mentioned opportunities and challenges, merge advanced measurement and passed The intelligent grid (Smart Grid) of multiple cutting edge technology such as sense technology, control technology, communication technology, computer technology becomes existing Only way for electric power industry development.One of important function of intelligent grid is it is simply that improve visualization and the early warning energy of electrical network Power, finally realizes intelligentized closed loop control, makes operation of power networks more safe and reliable, economical.And intelligent scheduling is intelligent grid The core content built, intelligent scheduling technology supports that system will have " with dispatcher's thinking model as framework, with visualization interface For functional module, with interaction as the system core " function.

It is a grand system engineering that intelligent grid wants to realize scale grid line security evaluation with early warning, in conjunction with current Ripe WAMS system and be based on SCADA（Supervisory Control And Data Acquisition, data acquisition with Monitoring）\EMS（Energy Management System, EMS）DSA with forecast failure collection emulation（Dynamic Security Assessment, dynamic secure estimation）System, for different operation of power networks scenes, needs to further investigate comprehensively Bulk power grid based on response message full situation quantitative estimation method data fusion treatment mechanism, sets up and is based on multiple response information source The current state of bulk power grid and the anticipation safety on line evaluation system that combines of state.

Content of the invention

The present invention relates to a kind of full situation of bulk power grid based on response integrated quantitative estimation method online, including：

Step S1, obtains topological structure of electric information from SCADA system, EMS system, and sets up and WAMS system electrical network unit The corresponding relation of part；

Step S2, obtains electrical network current method of operation trend number from described SCADA system, EMS system or WAMS system According to or obtain trend or the transient fault time domain data of various anticipations from DSA system；

Step S3, carries out steady stability Situation Assessment to electrical network in the way of online nodes oriented, with online towards unit To mode transient stability Situation Assessment is carried out to electrical network；

Wherein, the content carrying out steady stability Situation Assessment to electrical network includes：Generator stable margin index, circuit are stable Margin index, node stability margin index, generator thermal stablize qualification rate, the thermally-stabilised qualification rate of circuit, node voltage qualification rate, Load bus power factor qualification rate, steady stability situation synthesis index and the horizontal index of node reactive-load compensation；

The content carrying out transient stability Situation Assessment to electrical network includes：Transient stability margin index, transient stability estimate finger Mark, node voltage keep qualification rate, node stability margin index, the thermally-stabilised qualification rate of circuit and transient stability situation synthesis to refer to Mark.

In the first preferred embodiment that the present invention provides：Various electrical network steady-error coefficient data include envision the offline trend of state, N-1, N-2 trend convergence after result of calculation, SCADA EMS system state estimation result, PMU measurement information；

Various electrical network transient state transient process response datas include various anticipation transient fault collection time-domain simulation results, WAMS in fact When the grid disturbance procedural information that measures；

According to its previous moment or the change of periphery P MU metric data, reasonability anticipation or filter are carried out for PMU measurement information Ripple is processed；

For N-1, N-2 forecast failure, it is allocated according to parallel machine quantity.

In the second preferred embodiment that the present invention provides：In described step S3, electrical network is carried out steady stability Situation Assessment and During transient stability Situation Assessment, all for electrical network nodes are allocated by parallel machine quantity；

According to operation of power networks trend section information and combine topological structure of electric, try to achieve each node equivalent by active flow direction defeated Electrical power P_E+jQ_E；

Using the tracking parameter discrimination method being measured based on local, equivalent defeated under the current state of each node of on-line identification electrical network Electrical model parameters, described equivalent transmission of electricity model parameter comprises equivalent source potential E_E, equivalent branch impedance mould Z_E, impedance angle α, real Now the node equivalent transmission of electricity model reflection based on electrical network current operating conditions is penetrated.

In the third preferred embodiment that the present invention provides：The described generator stable nargin of electrical network is obtained in described step S3 The method of index, circuit stability margin index and node stability margin index is：

According to maximum transmitted power thought, it is respectively directed to electromotor branch road, interconnection branch road and node equivalent transmission of electricity model, Obtain described electromotor branch road, interconnection branch road or the current method of operation of node corresponding steady stability margin index S_Li, S_Li Computing formula is： $S_{\mathrm{Li}}=\frac{Z_L-Z_E}{Z_L};$

Wherein Z_LRepresent described electromotor branch road, interconnector or node equivalent load impedance, according to terminal voltage and work( Rate is tried to achieve.

In the 4th preferred embodiment that the present invention provides：The described generator thermal obtaining electrical network in described step S3 is surely qualified The method of rate and the steady qualification rate of interconnection heat includes：

According to the thermally-stabilised operation constraint of each electromotor, interconnection itself, statistics generator thermal stablizes qualification rate G_PRAnd line Thermally-stabilised qualification rate L in road_PRIt is respectively： $G_{\mathrm{PR}}=\frac{G_{\mathrm{QN}}}{G_N};$ $L_{\mathrm{PR}}=\frac{L_{\mathrm{QN}}}{L_N};$

Wherein, G_N、L_NIt is respectively electrical network total electromotor number, always get in touch with line number, G_QN、L_QNRespectively meet each self-operating heat steady Conclude a contract or treaty electromotor, the interconnection number of bundle；

The described node voltage qualification rate of electrical network and the side of load bus power factor qualification rate is obtained in described step S3 Method includes：

The node voltage being given according to electrical network normal operating mode and load bus power factor acceptability limit bound, system Meter grid nodes rate of qualified voltageLoad bus power factor qualification rate

Wherein N_N、P_FNIt is respectively the total nodes of electrical network and load bus number, V_QN、P_QFNRespectively meet the section of voltage range Points and the load bus number meeting power factor scope.

Obtain the quiet surely out-of-limit rate method of electrical network to include：

Set S_LiA certain early warning threshold value S_LC, count electrical network all interconnections branch road and node equivalent branch road be more than S_LC? Number, node includes electromotor node, middle contact node and load bus, calculates the quiet steady qualification rate of circuit respectively The quiet steady qualification rate of node $N_{\mathrm{SR}}=\frac{N_{\mathrm{CN}}}{N_N};$

Wherein L_N、N_NThe system of being respectively always gets in touch with line number and nodes, L_CN、N_CNIt is respectively and be more than S_LCContact line number and section Points.

In the 5th preferred embodiment that the present invention provides：Electrical network entirety steady stability situation synthesis are obtained in described step S3 Refer to calibration method to include：

Electrical network entirety steady stability situation synthesis index S_SI=α₁H_QR+β₁R_QR+χ₁S_SR；

Wherein, α₁、β₁、χ₁It is respectively electrical network to press thermally-stabilised constraint, run the weight system of electric parameters constraint and quiet steady constraint Number, rule of thumb or analysis need flexible configuration, α₁、β₁、χ₁Meet：

$\left\{\begin{array}{c}0\≤{\mathrm{\α}}_1\≤1\\ 0\≤{\mathrm{\β}}_1\≤1\\ 0\≤{\mathrm{\χ}}_1\≤1\\ 0\≤{\mathrm{\α}}_1+{\mathrm{\β}}_1+{\mathrm{\χ}}_1\≤1\end{array}\right.$

H_QRFor the steady qualification rate of electrical network overall thermal, H_QRComputing formula is：

R_QRFor electrical network overall operation electric parameters qualification rate, R_QRComputing formula is：

S_SRFor electrical network integrally quiet steady qualification rate, S_SRComputing formula is：

In the 6th preferred embodiment that the present invention provides：Being tried to achieve using the method for virtual idle change in described step S3 is made Equivalent transmission of electricity model corresponds to required idle size during active loss minimum, obtains idle benefit under the current method of operation of each node The level of repaying refers to calibration method：

In conjunction with described each node equivalent transmission of electricity model and its equivalent transmission of electricity active-power P_E, obtain described node equivalent model Lower circulation P_ENo work value Q that corresponding active loss needs when minimum_min, then contrast described Q_minWith the described node tried to achieve etc. Effect transmission of electricity reactive power Q_E, try to achieve described reactive-load compensation vacancy index Q under the current method of operation of described each node_IFor：

$Q_I=\frac{Q_E-Q_{\min }}{Q_E}$

Wherein, Q_I>Idle deficient benefit, Q is represented when 0_IRepresent idle optimal when=0, Q_I<Indicate when 0 that no merits and demerits are mended.

In the 7th preferred embodiment that the present invention provides：Temporary from described WAMS system or various anticipation in described step S3 Obtain electrical network transient process track time domain data in state fault set, directly adopt the method envelope electrical network of disturbed serious unit pair temporary State behavior analysiss, obtain the described transient stability margin index of electrical network, including step S101- step S103：

Step S101, inertia time constant M according to each electromotor and its angular velocity omega, generator rotor angle δ, generator terminal busbar voltage phase Angle θ and generator terminal electromagnetism active-power P_GEThe situation of change in moment in front and back, quickly identifies disturbed relatively advanced x platform and the most delayed x Platform electromotor, constitutes disturbed most serious unit set omega relatively；

Described numeral x can be manually set as needed, and the scope of x is 5≤x≤10；

Step S102, for any one the lead and lag unit pair in described set omega, calculates generating set i and sends out Transient stability margin index between group of motors j

Wherein, δ_ijRepresent generator rotor angle δ between described generating set i and generating set j_iAnd δ_jDifference, δ_ij=δ_i-δ_j；δ_seq For the mechanical output function P between described generating set i and generating set j_MeqWith electromagnetism active-power P_EeqIntersection point be virtual Stable equilibrium operating point；

Step S103, calculates between any lead and lag generating set i and generating set j in described set omega successively Transient stability margin index T_SIij, with the T of minimum unit pair_SIijTransient stability margin index as this time domain data S_TSI, i.e. S_TSI=min { T_SIij,i,j∈Ω}.

In the 8th preferred embodiment that the present invention provides：Calculate described δ_seqMethod include step S1021- step S1023：

Step S1021, the equivalent unit rotor inertia time constant between described generating set i and described generating set j ForCalculate the mechanical output function between described generating set i and described generating set j $P_{\mathrm{Meq}}=\frac{M_j}{M_i+M_j}{P}_{\mathrm{Mi}}-\frac{M_i}{M_i+M_j}{P}_{\mathrm{Mj}};$

Wherein M_iAnd M_jIt is respectively the inertia time constant of described generating set i and generating set j, P_MiAnd P_MjIt is respectively institute State the power that is mechanically injected of generating set i and generating set j, described generating set be mechanically injected power

Step S1022, calculates the electromagnetism active power between described generating set i and described generating set j $P_{\mathrm{Eeq}}=\frac{M_j}{M_i+M_j}{P}_{\mathrm{Ei}}-\frac{M_i}{M_i+M_j}{P}_{\mathrm{Ej}}=P_{\mathrm{EM}}\sin ({\mathrm{\&delta;}}_{\mathrm{ij}}+C);$

Wherein,A=Dcos (δ_j-θ_i)+Ecos(δ_i-θ_j), B=Dsin (δ_j-θ_i)-Esin(δ_i- θ_j), C=tan^-1(B/A), U_i、U_j、θ_i、θ_j、X_i、X_jRepresent electromotor respectively The generator terminal busbar voltage amplitude of group i and generating set j, phase angle and equivalent internal reactance；

Step S1023, calculates ${\mathrm{\&delta;}}_{\mathrm{seq}}=a\sin \left(\frac{P_{\mathrm{Meq}}}{P_{\mathrm{EM}}}\right)-C.$

In the 9th preferred embodiment that the present invention provides：For any lead and lag unit pair in set omega, foundation Unit after calculating is to equivalent mechanical output P_MeqWith electromagnetism active-power P_EeqChange path, using curve fitting technique and from energy The angle of amount builds transient stability and estimates index, specifically includes step S201- step S203：

Step S201 is it is assumed that described generating set is to equivalent mechanical output P_MeqConstant in transient process, described electromotor Group is to equivalent electromagnetism active-power P_EeqChange path is fitted with following SIN function：Y=x₀sin(δ_ij+x₁)+x₂；

Wherein, x₀, x₁, x₂For matching SIN function coefficient to be asked for, y represents generating set to equivalent P_Eeq；

Step S202, the transient stability calculating between described generating set i and generating set j estimates index T_SEIij；

Step S203, the T of any lead and lag unit pair in set of computations Ω successively_SEIij, with minimum unit pair T_SEIijTransient stability as this time domain data estimates index S_TEI, S_TEI=min { T_SEIij,i,j∈Ω}；

Calculate the transient stability between described generating set i and generating set j in described step S202 and estimate index and include Step S2021- step S2023：

Step S2021, according to P_EeqChange path fitting function y, tries to achieve the unstable equilibrium point δ after failure removal_iju, institute Stating unstable equilibrium point is described generating set to equivalent mechanical output P_MeqWith electromagnetism active-power P_EeqThe intersection point of change path, δ_ijuComputing formula is： ${\mathrm{\&delta;}}_{\mathrm{iju}}=\mathrm{\&pi;}-\arctan \frac{P_{\mathrm{Meq}}-x_2}{x_0}-x_1;$

Step S2022, the transient state kinetic energy V of failure removal moment unit pair_TAFor：

$V_{\mathrm{TA}}={\&Integral;}_{{\mathrm{\&delta;}}_{\mathrm{ij}0}}^{{\mathrm{\&delta;}}_{\mathrm{ijc}}}[P_{\mathrm{Meq}}-x_0\sin ({\mathrm{\&delta;}}_{\mathrm{ij}}+x_1)-x_2]d{\mathrm{\&delta;}}_{\mathrm{ij}}$

$=(P_{\mathrm{Meq}}-x_2)({\mathrm{\&delta;}}_{\mathrm{ijc}}-{\mathrm{\&delta;}}_{\mathrm{ij}0})+x_0[\cos ({\mathrm{\&delta;}}_{\mathrm{ijc}}+x_1)-\cos ({\mathrm{\&delta;}}_{\mathrm{ij}0}+x_1)]$

The critical potential that after failure removal, system is absorbed（Retardation area）V_TBFor：

$V_{\mathrm{TB}}={\&Integral;}_{{\mathrm{\&delta;}}_{\mathrm{ijc}}}^{{\mathrm{\&delta;}}_{\mathrm{iju}}}[x_0\sin ({\mathrm{\&delta;}}_{\mathrm{ij}}+x_1)+x_2-P_{\mathrm{Meq}}]d{\mathrm{\&delta;}}_{\mathrm{ij}}$

$=(x_2-P_{\mathrm{Meq}})({\mathrm{\&delta;}}_{\mathrm{iju}}-{\mathrm{\&delta;}}_{\mathrm{ijc}})-x_0[\cos ({\mathrm{\&delta;}}_{\mathrm{iju}}+x_1)-\cos ({\mathrm{\&delta;}}_{\mathrm{ijc}}+x_1)]$

Wherein, δ_ij0And δ_ijcThe unit in expression fault generation moment and failure removal moment is poor to generator rotor angle respectively.

Step S2023, the described transient stability calculating between described generating set i and described generating set j estimates index T_SEIij： $T_{\mathrm{SEIij}}=1-\frac{V_{\mathrm{TA}}}{V_{\mathrm{TB}}};$

Wherein, T_SEIijIndex transient stability meaning is：T_SEIij>0 transient stability, T_SEIij=0 critical transient stability, T_SEIij< 0, Transient Instability.

In the tenth preferred embodiment that the present invention provides：Obtain grid nodes voltage in described step S3 and keep qualification rate Method includes：

Given transient process interior joint voltage holding capability investigates scope：Voltage Drop lower limit V_LWith persistent period V_C, statistics The all nodes of electrical network are in the interstitial content V of the horizontal tolerance interval of transient process fullness in the epigastrium and abdomen afc voltage_QN, calculate in described transient process Node voltage meets the qualification rate of holding capacity

The method obtaining the steady qualification rate of the heat of circuit described in electrical network in described step S3 includes：Given transient process interior joint Steady stability margin index and circuit heat steady investigation scope：The quiet steady nargin limit value S of node_IC, the steady limit value L of circuit heat_PRWith lasting when Between T_C, count the interstitial content S that electrical network meets quiet surely sustainable scope_QN, meet number of lines L of the surely sustainable scope of heat_QN, meter Calculate qualification rate S that described transient process interior joint is quiet steady and circuit heat is steady_TR, S_TRComputing formula is：

In the 11st preferred embodiment that the present invention provides：Described electrical network entirety transient stability state is obtained in described step S3 Gesture aggregative indicator $S_{\mathrm{TI}}={\mathrm{\&alpha;}}_2{S_T}_{\mathrm{SEI}}+{\mathrm{\&beta;}}_2{V}_{\mathrm{TR}}+{\mathrm{\&chi;}}_2{S}_{\mathrm{TR}};$

Wherein, α₂、β₂、χ₂It is respectively electrical network to press thermally-stabilised constraint, run the weight system of electric parameters constraint and quiet steady constraint Number, rule of thumb or analysis need flexible configuration, α₂、β₂、χ₂Meet：

$\left\{\begin{array}{c}0\&le;{\mathrm{\&alpha;}}_2\&le;1\\ 0\&le;{\mathrm{\&beta;}}_2\&le;1\\ 0\&le;{\mathrm{\&chi;}}_2\&le;1\\ 0\&le;{\mathrm{\&alpha;}}_2+{\mathrm{\&beta;}}_2+{\mathrm{\&chi;}}_2\&le;1\end{array}\right.$

A kind of beneficial effect of full situation of bulk power grid based on response that the present invention provides integrated quantitative estimation method online Fruit includes：

1st, a kind of full situation of bulk power grid based on response that the present invention provides integrated quantitative estimation method online, can be directed to Various emulation data under the real-time measurement of the current state of electrical network or anticipation state are it is achieved that bulk power grid is static and two kinds of typical cases of transient state The full situation of Run-time scenario integrated quantitative evaluation online.Achieve the Situation Assessment under the current state of electrical network, also fully combine Existing DSA envisions the Situation Assessment that emulation advantage achieves anticipation state, is conducive to management and running personnel to understand that electrical network is worked as in time Front operation situation and the potential system risk of electrical network.Also various offline anticipation modes or fault can be applied indirectly to Sample Method The later stage intelligent assessment of collection time-domain-simulation data, greatly reduces planning personnel or mode formulates the workload of personnel.This Bright core methed used is all based on grid responsive data, has stronger independence, and full method for situation assessment and algorithm directly, Simply, weak node or the weakness zone under any method of operation of electrical network can quickly be identified, to run or planning personnel's reference, Prevent trouble before it happens, voltage stabilization accident is reduced or avoided and occurs, be suitable for online engineer applied.The present invention can make conventional " modeling Defense mode is progressively transitioned into " track excavation " response control pattern, is the extension of traditional electrical network prevention and control thought and method for emulation ", Scale grid line intelligent assessment and early warning level can be effectively improved.

2nd, obtain topological structure of electric relation from EMS system, be so conducive to the later stage to be directed to various response datas and carry out Electrical network full situation quantitative evaluation and statistical analysiss.

3rd, from the point of view of response angle, electric network element and voltage levvl are directly directly carried out simple statistics from acceptability limit and divides Analysis.And in traditional angle stability and voltage stabilization situation qualitative assessment, static stability focuses on grid power transmission Ability, transient stability focuses on preservation of energy, sets up the full situation of the bulk power grid online unified quantization assessment thought based on response.

4th, arbitrarily measure in real time for electrical network or the trend section of various forecast failure form carries out steady stability form Situation Assessment, when concrete node is stablized with situation or fault set analysis, the method using tasks in parallel improves comprehensive assessment efficiency.

5th, all for electrical network nodes are pressed current trend distribution and topological relation, count this node actual circulation（Or transfer） " single supply list load " that node equivalent is simplification is simplified transmission of electricity model by power based on this, and using local measurement letter The online tracking of breath recognizes its equivalent transmission of electricity parameter, realizes the equivalent dummy model mapping of all nodes of electrical network, is that next step is static Stablize Situation Assessment to lay the foundation.

6th, the impedance modular character being used uniformly across expressing total transfer capability thought carries out steady stability Situation Assessment, and will Modulus of impedance method is applied in electromotor branch road, power transmission line and node equivalent transmission of electricity model.Modulus of impedance method is suitable for quick online Calculate, and can effectively identify weak unit under the current method of operation, critical circuits and weak node.

7th, count the current method of operation of electrical network from the thermally-stabilised angle running acceptability limit of electromotor and circuit The steady qualification rate of heat, can fast and effeciently screen overload element, also comply with the custom of conventional Load Flow mode visual assessment.

8th, from node voltage levvl and load bus power factor angle are set out the electricity under the statistics current method of operation of electrical network Tolerance qualification rate, can fast and effeciently evaluate current method of operation voltage and the qualified level of power factor, also comply with conventional Load Flow Mode quality of power supply check requirements.

9th, nodes oriented equivalent transmission of electricity model and circuit, realizes steady stability qualification rate statistics using impedance modular character.

10th, fully combine statistical result it is achieved that from component-level heat is steady, transportation load acceptability limit is to the quiet steady synthesis of system Property quiet stablize very much Situation Assessment, and can need flexibly to select three's weight according to actual motion state and assessment emphasis.

11st, the node equivalent transmission of electricity model by construction and its parameter, can try to achieve when node currently circulates active so that having Work(is lost the virtual idle size needed for minimum, and then electrical network all nodes reactive-load compensation level under the current method of operation of realizing is commented Estimate

12nd, carry out transient stability margin for transient state time domain procedures information and when transient stability estimates, all directly adopt The dynamic behaviour envelope reflection electrical network overall dynamics behavior of disturbed most serious unit pair, need not be processed using point group's equivalence etc., thing Reason meaning is clearly directly perceived, and assessment result degree of accuracy is high.

13rd, in order to reduce disturbed most serious unit to the error selecting, disturbed most serious relatively are selected using multiple angles Unit set omega, then select the unit pair of most serious from set omega, because of the weakest unit, electrical network is caused to selection error to reduce The erroneous judgement of actual situation.

14th, for the lead and lag unit pair in arbitrary collection Ω, according to each unit electric quantity information of itself, pass through Simple algebraic manipulation tries to achieve unit to transient stability margin index T_SI, and the T with unit pair minimum in set omega_SIAs this The transient stability margin index S of time domain data_TSI, and then realize the transient stability margin quantitative evaluation of current transient track.

15th, adopt trigonometric function form matching unit reflection reasonable to equivalent generator rotor angle function curve electrical network transient overshoot Track characteristic, can easily calculate energy type transient stability by this matched curve and estimate index T_SEI, physical significance is clear Intuitively, the T and with unit pair minimum in set omega_SEITransient stability as this time domain data estimates index S_TEI, Jin Ershi The transient stability of existing current transient track quantifies to estimate.

16th, use for reference Transient Voltage Stability habituation determination methods, weighed using a certain restriction of Voltage Drop and persistent period Node voltage holding capacity in the given transient process of amount and qualification rate level.

17th, the same method using " limit value+persistent period ", and combine quiet steady in certain methods, try to achieve transient process Qualification rate S that interior joint is quiet steady and circuit heat is steady_TR, be easy to assess node static stability margin in given transient state transient process and Line load level condition.

18th, fully combine statistical result it is achieved that quiet from dynamical system preservation of energy, transient voltage holding capacity, node Comprehensive transient stability Situation Assessment index S steady and that circuit heat is steady_TI, spirit is needed according to actual dynamic trajectory and assessment emphasis Live and select weight, calculate the corresponding S of different faults_TIAnd realize their order of severity sequence, improve forecast failure collection and divide online Analysis efficiency, S_TIEffectively the stable situation of multiple transient faults can be carried out with effective quantitative evaluation and be easy to serious transient fault sieve Choosing.

19th, a large amount of forecast failures are concentrated to carry out transient stability Study on Trend transient fault using the method for tasks in parallel.

Brief description

A kind of full situation of bulk power grid based on response integrated quantitative evaluation number online of present invention offer is provided According to flow graph；

The comprehensive index system elemental map of present invention offer is provided；

The electrical network local nodes schematic diagram of present invention offer is provided；

The node equivalent transmission of electricity model schematic of present invention offer is provided；

The maximum transmitted power border modulus of impedance schematic diagram of present invention offer is provided；

The node reactive-load compensation proficiency assessment schematic diagram of present invention offer is provided；

The electrical network unit being illustrated in figure 7 present invention offer is to schematic diagram；

The unit being illustrated in figure 8 present invention offer is to equivalent power-angle curve matching schematic diagram；

The energy type transient stability being illustrated in figure 9 present invention offer estimates index schematic diagram；

The integrated quantitative evaluation flow process online of the full situation of bulk power grid based on response providing for the present invention as shown in Figure 10 Figure.

Specific embodiment

Below in conjunction with the accompanying drawings the specific embodiment of the present invention is described in further detail.

The present invention provides a kind of full situation of bulk power grid based on response integrated quantitative estimation method online, as shown in Figure 1 A kind of full situation of bulk power grid based on response that the present invention provides integrated quantitative evaluation data flow diagram online.As shown in Figure 1, should Appraisal procedure includes：

Step S1, obtains topological structure of electric information from SCADA system, EMS system, and sets up and WAMS system electrical network unit The corresponding relation of part.

Step S2, from SCADA system, EMS system or WAMS system obtain electrical network current method of operation flow data or Person obtains trend or the transient fault time domain data of various anticipations from DSA system.

Step S3, carries out steady stability Situation Assessment to this electrical network in the way of online nodes oriented, with online towards machine Group to mode transient stability Situation Assessment is carried out to this electrical network.

Build situation synthesis evaluation index system when considered electric network element level thermostability, run electric parameters qualified Each factor with system-level stability three, is illustrated in figure 2 the comprehensive index system elemental map of present invention offer.

Wherein, the content carrying out steady stability Situation Assessment to electrical network includes：Generator stable margin index, circuit are stable Margin index, node stability margin index, generator thermal stablize qualification rate, the thermally-stabilised qualification rate of circuit, node voltage qualification rate, Load bus power factor qualification rate, steady stability situation synthesis index and the horizontal index of node reactive-load compensation.

The content carrying out transient stability Situation Assessment to electrical network includes：Transient stability margin index, transient stability estimate finger Mark, node voltage keep qualification rate, node stability margin index, the thermally-stabilised qualification rate of circuit and transient stability situation synthesis to refer to Mark.

Various electrical network steady-error coefficient data include envision the offline trend of state, N-1, N-2 trend convergence after result of calculation, SCADA EMS system state estimation result, PMU measurement information；Various electrical network transient state transient process response datas include various pre- Think the grid disturbance procedural information of transient fault collection time-domain simulation results, WAMS measurement in real time；For PMU measurement information according to it Previous moment or the change of periphery P MU metric data carry out reasonability anticipation or Filtering Processing；For N-1, N-2 forecast failure, root It is allocated according to parallel machine quantity.

When steady stability Situation Assessment and transient stability Situation Assessment being carried out to this electrical network in step S3, by all for electrical network sections Press parallel machine quantity to be allocated, be illustrated in figure 3 the present invention offer electrical network local nodes schematic diagram, node may with send out The Various Components such as motor, load, capacitor, circuit are connected.

According to operation of power networks trend section information and combine topological structure of electric, try to achieve each node equivalent by active flow direction defeated Electrical power P_E+jQ_E, it is illustrated in figure 4 node equivalent transmission of electricity model schematic.Using the tracking parameter identification being measured based on local Method, the equivalent transmission of electricity model parameter under the current state of each node of on-line identification electrical network, mainly comprise equivalent source potential E_E, equivalent Branch impedance mould Z_E, impedance angle α, realize based on electrical network current operating conditions node equivalent transmission of electricity model reflection penetrates.

Generator stable margin index, circuit stability margin index and the node stability margin of this electrical network is obtained in step S3 Refer to calibration method in, using maximum transmitted power thought as shown in Figure 5 impedance modular character as steady stability margin index, Including：

According to maximum transmitted power thought, it is respectively directed to electromotor branch road, interconnection branch road and node equivalent transmission of electricity model, Obtain this electromotor branch road, interconnection branch road or the current method of operation of node corresponding steady stability margin index S_Li, S_LiMeter Calculating formula is：

$S_{\mathrm{Li}}=\frac{Z_L-Z_E}{Z_L}$

Z_LEach can represent electromotor branch road, interconnector or node equivalent load impedance, can be according to terminal voltage and work( Rate is tried to achieve.S_LiWeigh the size cases of each electromotor branch road, transmission of electricity linear load and node static stability margin, and can be used to position Weak unit, critical circuits and weak node.

Specifically, by the steady stability margin index S of all nodes_LiCarry out average, electricity under the current method of operation can be tried to achieve Net averaged static stability margin index S_AI, computing formula is：

$S_{\mathrm{AI}}=\frac{\underset{i=1}{\overset{N}{\mathrm{\&Sigma;}}}{S}_{\mathrm{Li}}}{N_N}$

The method of the steady qualification rate of generator thermal and the steady qualification rate of interconnection heat that obtain this electrical network in step S3 includes：

According to the thermally-stabilised operation constraint of each electromotor, interconnection itself, statistics generator thermal stablizes qualification rate G_PRAnd line Thermally-stabilised qualification rate L in road_PRIt is respectively： $G_{\mathrm{PR}}=\frac{G_{\mathrm{QN}}}{G_N};$ $L_{\mathrm{PR}}=\frac{L_{\mathrm{QN}}}{L_N}.$

Wherein, G_N、L_NIt is respectively electrical network total electromotor number, always get in touch with line number, G_QN、L_QNRespectively meet each self-operating heat steady Conclude a contract or treaty electromotor, the interconnection number of bundle.

The method of the node voltage qualification rate and load bus power factor qualification rate that obtain this electrical network in step S3 includes：

The node voltage being given according to electrical network normal operating mode and load bus power factor acceptability limit bound, system Meter grid nodes rate of qualified voltageLoad bus power factor qualification rate

Wherein N_N、P_FNIt is respectively the total nodes of electrical network and load bus number, V_QN、P_QFNRespectively meet the section of voltage range Points and the load bus number meeting power factor scope.

Obtain the quiet surely out-of-limit rate method of electrical network to include：

Set S_LiA certain early warning threshold value S_LC, count electrical network all interconnections branch road and node（Electromotor node, middle connection Network node, load bus）Equivalent branch road is more than S_LCNumber, respectively calculate the quiet steady qualification rate of circuitNode is quiet steady Qualification rate $N_{\mathrm{SR}}=\frac{N_{\mathrm{CN}}}{N_N};$

Wherein L_N、N_NThe system of being respectively always gets in touch with line number and nodes, L_CN、N_CNIt is respectively and be more than S_LCContact line number and section Points.

Obtain this electrical network entirety steady stability situation synthesis in step S3 and refer to calibration method including：Electrical network entirety steady stability Situation synthesis index S_I=α₁H_QR+β₁R_QR+χ₁S_SR.

Wherein, α₁、β₁、χ₁It is respectively electrical network to press thermally-stabilised constraint, run the weight system of electric parameters constraint and quiet steady constraint Number, can rule of thumb or analysis need flexible configuration.α₁、β₁、χ₁Meet following condition：

$\left\{\begin{array}{c}0\&le;{\mathrm{\&alpha;}}_1\&le;1\\ 0\&le;{\mathrm{\&beta;}}_1\&le;1\\ 0\&le;{\mathrm{\&chi;}}_1\&le;1\\ 0\&le;{\mathrm{\&alpha;}}_1+{\mathrm{\&beta;}}_1+{\mathrm{\&chi;}}_1\&le;1\end{array}\right.$

H_QRFor the steady qualification rate of electrical network overall thermal, H_QRComputing formula is:

R_QRFor electrical network overall operation electric parameters qualification rate, R_QRComputing formula is：

S_SRFor electrical network integrally quiet steady qualification rate, S_SRComputing formula is：

S_SIThe comprehensive electrical network steady stability situation level weighed under current flow state, can be effectively to multiple trend modes Steady stability situation carries out effective quantitative evaluation.

Tried to achieve using the method for virtual idle change in step S3 and make equivalent transmission of electricity model corresponding active loss minimum when institute The idle size needing, node reactive-load compensation proficiency assessment schematic diagram is as shown in fig. 6, obtain under the current method of operation of each node no Work(compensation level refers to calibration method and includes：

In conjunction with each node equivalent transmission of electricity model and its equivalent transmission of electricity active-power P_E, circulation under node equivalent model can be obtained P_ENo work value Q that corresponding active loss needs when minimum_min, then contrast Q_minWith the node equivalent transmission of electricity reactive power tried to achieve Q_E, reactive-load compensation vacancy index Q under the current method of operation of each node can be tried to achieve_IFor：

$Q_I=\frac{Q_E-Q_{\min }}{Q_E}$

Wherein, Q_I>Idle deficient benefit, Q is represented when 0_IRepresent idle optimal when=0, Q_I<Indicate when 0 that no merits and demerits are mended.

Concentrate from WAMS system or various anticipation transient fault in step S3 and obtain electrical network transient process track time domain number According to, directly adopt disturbed serious unit pair method envelope electrical network transient behavior analysis, be illustrated in figure 7 electrical network unit to signal Figure, obtains the transient stability margin index of this electrical network, including step S101- step S103：

Step S101, inertia time constant M according to each electromotor and its angular velocity omega, generator rotor angle δ, generator terminal busbar voltage phase Angle θ and generator terminal electromagnetism active-power P_GEThe situation of change in moment in front and back, quickly identifies disturbed relatively advanced x platform and the most delayed x Platform electromotor, constitutes disturbed most serious unit set omega relatively.

This digital x can be manually set as needed, and the scope of x is 5≤x≤10.

Step S102, for any one the lead and lag unit pair in set omega, calculates generating set i and electromotor Transient stability margin index between group j

Wherein, δ_ijRepresent generator rotor angle δ between generating set i and generating set j_iAnd δ_jDifference, δ_ij=δ_i-δ_j.

δ_seqFor the mechanical output function P between generating set i and generating set j_MeqWith electromagnetism active-power P_EeqIntersection point Virtual stable equilibrium, operating point, calculated δ_seqMethod include step S1021- step S1023：

Step S1021, the equivalent unit rotor inertia time constant between generating set i and generating set j isCalculate the mechanical output function between generating set i and generating set j

Wherein M_iAnd M_jIt is respectively the inertia time constant of generating set i and generating set j, P_MiAnd P_MjIt is respectively electromotor Group i and generating set j is mechanically injected power, generating set be mechanically injected power

Step S1022, calculates the electromagnetism active power between generating set i and generating set j $P_{\mathrm{Eeq}}=\frac{M_j}{M_i+M_j}{P}_{\mathrm{Ei}}-\frac{M_i}{M_i+M_j}{P}_{\mathrm{Ej}}=P_{\mathrm{EM}}\sin ({\mathrm{\&delta;}}_{\mathrm{ij}}+C).$

Wherein,A=Dcos (δ_j-θ_i)+Ecos(δ_i-θ_j), B=Dsin (δ_j-θ_i)-Esin(δ_i- θ_j), C=tan^-1(B/A), θ_iAnd θ_jRepresent generating set i and generating respectively The generator terminal busbar voltage phase angle U of unit j_i、U_j、θ_i、θ_j、X_i、X_jRepresent the generator terminal bus of generating set i and generating set j respectively Voltage magnitude, phase angle and equivalent internal reactance.

Step S1023, calculates ${\mathrm{\&delta;}}_{\mathrm{seq}}=a\sin \left(\frac{P_{\mathrm{Meq}}}{P_{\mathrm{EM}}}\right)-C.$

Step S103, calculates between any lead and lag generating set i and generating set j in described set omega successively Transient stability margin index T_SIij, with the T of minimum unit pair_SIijTransient stability margin index as this time domain data S_TSI, i.e. S_TSI=min { T_SIij,i,j∈Ω}.S by minimum unit pair_TSICan achieve the positioning of weakness unit.

Obtain in step S3 this electrical network transient stability estimate finger calibration method include：For arbitrarily advanced in set omega With delayed unit pair, the generating set after foundation calculating is to equivalent mechanical output P_MeqWith electromagnetism active-power P_EeqChange path, Using curve fitting technique and build transient stability from the angle of energy and estimate index, be illustrated in figure 8 unit to equivalent work( Angular curve matching schematic diagram.Obtain this electrical network transient stability estimate finger calibration method specifically include step S201- step S203：

Step S201, generating set is to equivalent mechanical output P_MeqTransient process can be assumed that constant, generating set equity Effect electromagnetism active-power P_EeqChange path can use following SIN function to be fitted：Y=x₀sin(δ_ij+x₁)+x₂.

Wherein, x₀, x₁, x₂For matching SIN function coefficient to be asked for, y represents generating set to equivalent P_Eeq.

Step S202, the transient stability calculating between generating set i and generating set j estimates index T_SEIij, including step S2021- step S2023：

Step S2021, according to P_EeqChange path fitting function y, tries to achieve the unstable equilibrium point δ after failure removal_iju, δ_ijuComputing formula is： ${\mathrm{\&delta;}}_{\mathrm{iju}}=\mathrm{\&pi;}-\arctan \frac{P_{\mathrm{Meq}}-x_2}{x_0}-x_1.$

Step S2022, the transient state kinetic energy of failure removal moment unit pair（Accelerate area）V_TAFor：

$V_{\mathrm{TA}}={\&Integral;}_{{\mathrm{\&delta;}}_{\mathrm{ij}0}}^{{\mathrm{\&delta;}}_{\mathrm{ijc}}}[P_{\mathrm{Meq}}-x_0\sin ({\mathrm{\&delta;}}_{\mathrm{ij}}+x_1)-x_2]d{\mathrm{\&delta;}}_{\mathrm{ij}}$

$=(P_{\mathrm{Meq}}-x_2)({\mathrm{\&delta;}}_{\mathrm{ijc}}-{\mathrm{\&delta;}}_{\mathrm{ij}0})+x_0[\cos ({\mathrm{\&delta;}}_{\mathrm{ijc}}+x_1)-\cos ({\mathrm{\&delta;}}_{\mathrm{ij}0}+x_1)]$

The critical potential that after failure removal, system is absorbed（Retardation area）V_TBFor：

$V_{\mathrm{TB}}={\&Integral;}_{{\mathrm{\&delta;}}_{\mathrm{ijc}}}^{\mathrm{\&delta;iju}}[x_0\sin ({\mathrm{\&delta;}}_{\mathrm{ij}}+x_1)+x_2-P_{\mathrm{Meq}}]d{\mathrm{\&delta;}}_{\mathrm{ij}}$

$=(x_2-P_{\mathrm{Meq}})({\mathrm{\&delta;}}_{\mathrm{iju}}-{\mathrm{\&delta;}}_{\mathrm{ijc}})-x_0[\cos ({\mathrm{\&delta;}}_{\mathrm{iju}}+x_1)-\cos ({\mathrm{\&delta;}}_{\mathrm{ijc}}+x_1)]$

Wherein, δ_ij0And δ_ijcThe unit in expression fault generation moment and failure removal moment is poor to generator rotor angle respectively.

Step S2023, energy type transient stability as shown in Figure 9 estimates index schematic diagram, calculates generating set i and generating Transient stability between unit j estimates index T_SEIij：

Wherein, T_SEIijIndex transient stability meaning is as follows, T_SEIij>0 transient stability, T_SEIij=0 critical transient stability, T_SEIij <0, Transient Instability.

Step S203, the T of any lead and lag unit pair in set of computations Ω successively_SEIij, with minimum unit pair T_SEIijTransient stability as this time domain data estimates index S_TEI, S_TEIComputing formula is：S_TEI=min { T_SEIij,i,j∈ Ω}.

Obtaining this grid nodes voltage in step S3 keeps the method for qualification rate to include：

Given transient process interior joint voltage holding capability investigates scope（Voltage Drop lower limit V_LWith persistent period V_C）, system The meter all nodes of electrical network are in the interstitial content V of the horizontal tolerance interval of transient process fullness in the epigastrium and abdomen afc voltage_QN, calculate in this transient process Node voltage meets qualification rate V of holding capacity_TR, V_TRComputing formula is：

The method obtaining the steady qualification rate of circuit heat in this electrical network in step S3 includes：Node static in given transient process Stability margin index and circuit heat steady investigation scope：The quiet steady nargin limit value S of node_IC, the steady limit value L of circuit heat_PRAnd duration T_C, Statistics electrical network meets the interstitial content S of quiet surely sustainable scope_QN, meet number of lines L of the surely sustainable scope of heat_QN, calculating should Qualification rate S that transient process interior joint is quiet steady and circuit heat is steady_TR, S_TRComputing formula is：

This electrical network entirety transient stability situation synthesis index S is obtained in step S3_TI=α₂S_TEi+β₂V_TR+χ₂S_TR；

Wherein, α₂、β₂、χ₂It is respectively electrical network to press thermally-stabilised constraint, run the weight system of electric parameters constraint and quiet steady constraint Number, can rule of thumb or analysis need flexible configuration.α₂、β₂、χ₂Meet following condition：

$\left\{\begin{array}{c}0\&le;{\mathrm{\&alpha;}}_2\&le;1\\ 0\&le;{\mathrm{\&beta;}}_2\&le;1\\ 0\&le;{\mathrm{\&chi;}}_2\&le;1\\ 0\&le;{\mathrm{\&alpha;}}_2+{\mathrm{\&beta;}}_2+{\mathrm{\&chi;}}_2\&le;1\end{array}\right.$

S_TIThe Power Network Transient Stability situation level of current transient procedural information can comprehensively be weighed, can be effectively to the event of multiple transient state The stable situation of barrier carries out effective quantitative evaluation.

Finally it should be noted that：Above example is only not intended to limit in order to technical scheme to be described, to the greatest extent Pipe has been described in detail to the present invention with reference to above-described embodiment, and those of ordinary skill in the art should be understood：Still The specific embodiment of the present invention can be modified or equivalent, and any without departing from spirit and scope of the invention Modification or equivalent, it all should be covered in the middle of scope of the presently claimed invention.

Claims (11)

1. a kind of full situation of bulk power grid based on response online integrated quantitative estimation method it is characterised in that described assessment side Method includes：

Step S1, obtains topological structure of electric information, and sets up and WAMS system electric network element from SCADA system, EMS system Corresponding relation；

Step S2, from described SCADA system, EMS system or WAMS system obtain electrical network current method of operation flow data or Person obtains trend or the transient fault time domain data of various anticipations from DSA system；

Step S3, carries out steady stability Situation Assessment to electrical network in the way of online nodes oriented, with online towards unit pair Mode carries out transient stability Situation Assessment to electrical network；

When steady stability Situation Assessment and transient stability Situation Assessment being carried out to electrical network in described step S3, by all for electrical network nodes It is allocated by parallel machine quantity；

According to operation of power networks trend section information and combine topological structure of electric, try to achieve each node equivalent transmission of electricity work(by active flow direction Rate P_E+jQ_E；Wherein, P_EFor node equivalent transmission of electricity active power, Q_EFor node equivalent transmission of electricity reactive power；

Using the tracking parameter discrimination method being measured based on local, the equivalent transmission of electricity mould under the current state of each node of on-line identification electrical network Shape parameter, described equivalent transmission of electricity model parameter comprises equivalent source potential E_E, equivalent branch impedance mould Z_E, impedance angle α, realize base Node equivalent transmission of electricity model reflection in electrical network current operating conditions is penetrated；

Wherein, the content carrying out steady stability Situation Assessment to electrical network includes：Generator stable margin index, circuit stability margin Index, node stability margin index, generator thermal stablize qualification rate, the thermally-stabilised qualification rate of circuit, node voltage qualification rate, load Node power factor qualification rate, steady stability situation synthesis index and the horizontal index of node reactive-load compensation；

The content carrying out transient stability Situation Assessment to electrical network includes：Transient stability margin index, transient stability estimate index, section Point voltage keeps qualification rate, node stability margin index, the thermally-stabilised qualification rate of circuit and transient stability situation synthesis index.

2. the method for claim 1 it is characterised in that various electrical network steady-error coefficient data include envision state damp offline Stream, N-1, N-2 trend convergence after result of calculation, SCADA EMS system state estimation result, PMU measurement information；

Various electrical network transient state transient process response datas include various anticipation transient fault collection time-domain simulation results, WAMS surveys in real time The grid disturbance procedural information of amount；

PMU measurement information is carried out at reasonability anticipation or filtering according to its previous moment or the change of periphery P MU metric data Reason；

For N-1, N-2 forecast failure, it is allocated according to parallel machine quantity.

3. the method for claim 1 is it is characterised in that the described generator stable obtaining electrical network in described step S3 is abundant Degree index, the method for circuit stability margin index and node stability margin index are：

According to maximum transmitted power thought, it is respectively directed to electromotor branch road, interconnection branch road and node equivalent transmission of electricity model, obtains Described electromotor branch road, interconnection branch road or the current method of operation of node corresponding steady stability margin index S_LiS_LiCalculate Formula is：

Wherein Z_LRepresent described electromotor branch road, interconnector or node equivalent load impedance, asked according to terminal voltage and power ?.

4. the method for claim 1 is it is characterised in that the described generator thermal obtaining electrical network in described step S3 is stable The method of qualification rate and the thermally-stabilised qualification rate of circuit includes：

According to the thermally-stabilised operation constraint of each electromotor, interconnection itself, statistics generator thermal stablizes qualification rate G_PRWith circuit heat Stablize qualification rate L_PRIt is respectively：

Wherein, G_N、L_NIt is respectively electrical network total electromotor number, always get in touch with line number, G_QN、L_QNRespectively meet each self-operating thermally-stabilised about The electromotor of bundle, interconnection number；

The described node voltage qualification rate of electrical network and the method bag of load bus power factor qualification rate is obtained in described step S3 Include：

The node voltage being given according to electrical network normal operating mode and load bus power factor acceptability limit bound, statistics electricity Net node voltage qualification rateLoad bus power factor qualification rate

Wherein N_N、P_FNIt is respectively the total nodes of electrical network and load bus number, V_QN、P_QFNRespectively meet the nodes of voltage range With the load bus number meeting power factor scope；

Obtain the quiet steady qualification rate method of electrical network to include：

Set S_LiA certain early warning threshold value S_LC, count electrical network all interconnections branch road and node equivalent branch road be more than S_LCNumber, section Point includes electromotor node, middle contact node and load bus, calculates the quiet steady qualification rate of circuit respectivelyNode is quiet Steady qualification rate

Wherein L_N、N_NThe system of being respectively always gets in touch with line number and nodes, L_CN、N_CNIt is respectively and be more than S_LCContact line number and node Number.

5. method as claimed in claim 4 it is characterised in that obtain in described step S3 electrical network entirety steady stability situation comprehensive Conjunction refers to calibration method and includes：

Electrical network entirety steady stability situation synthesis index S_SI=α₁H_QR+β₁R_QR+χ₁S_SR；

Wherein, α₁、β₁、χ₁It is respectively electrical network to press thermally-stabilised constraint, run the weight coefficient of electric parameters constraint and quiet steady constraint, root According to empirically or analytically needing flexible configuration, α₁、β₁、χ₁Meet：

$\left\{\begin{array}{c}0\&le;{\mathrm{\&alpha;}}_1\&le;1\\ 0\&le;{\mathrm{\&beta;}}_1\&le;1\\ 0\&le;{\mathrm{\&chi;}}_1\&le;1\\ 0\&le;{\mathrm{\&alpha;}}_1+{\mathrm{\&beta;}}_1+{\mathrm{\&chi;}}_1\&le;1\end{array}\right.$

H_QRFor the steady qualification rate of electrical network overall thermal, H_QRComputing formula is：

R_QRFor electrical network overall operation electric parameters qualification rate, R_QRComputing formula is：

S_SRFor electrical network integrally quiet steady qualification rate, S_SRComputing formula is：

6. the method for claim 1 is it is characterised in that tried to achieve using the method for virtual idle change in described step S3 Make required idle size when minimum of the corresponding active loss of equivalent transmission of electricity model, obtain idle under the current method of operation of each node Compensation level refers to calibration method：

In conjunction with described each node equivalent transmission of electricity model and its equivalent transmission of electricity active-power P_E, obtain described node equivalent model dirty Logical P_ENo work value Q that corresponding active loss needs when minimum_min, then contrast described Q_minDefeated with the described node equivalent tried to achieve Electric reactive power Q_E, try to achieve described reactive-load compensation vacancy index Q under the current method of operation of described each node_IFor：

$Q_I=\frac{Q_E-Q_{min}}{Q_E}$

Wherein, Q_I>Idle deficient benefit, Q is represented when 0_IRepresent idle optimal when=0, Q_I<Indicate when 0 that no merits and demerits are mended.

7. the method for claim 1 is it is characterised in that from described WAMS system or various anticipation in described step S3 Transient fault is concentrated and is obtained electrical network transient process track time domain data, directly adopts the method envelope electrical network of disturbed serious unit pair Transient behavior is analyzed, and obtains the described transient stability margin index of electrical network, including step S101- step S103：

Step S101, inertia time constant M according to each electromotor and its angular velocity omega, generator rotor angle δ, generator terminal busbar voltage phase angle theta With generator terminal electromagnetism active-power P_GEThe situation of change in moment in front and back, quickly identifies that disturbed relatively advanced x platform and the most delayed x platform are sent out Motor, constitutes disturbed most serious unit set omega relatively；

Described numeral x can be manually set as needed, and the scope of x is 5≤x≤10；

Step S102, for any one the lead and lag unit pair in described set omega, calculates generating set i and electromotor Transient stability margin index between group j

Wherein, δ_ijRepresent generator rotor angle δ between described generating set i and generating set j_iAnd δ_jDifference, δ_ij=δ_i-δ_j；δ_seqFor institute State the mechanical output function P between generating set i and generating set j_MeqWith electromagnetism active-power P_EeqIntersection point be virtual stable Balance movement point；

Step S103, calculates temporary between any lead and lag generating set i and generating set j in described set omega successively State stability margin index T_SIij, with the T of minimum unit pair_SIijTransient stability margin index S as this time domain data_TSI, that is, S_TSI=min { T_SIij,i,j∈Ω}.

8. method as claimed in claim 7 is it is characterised in that calculate described δ_seqMethod include step S1021- step S1023：

Step S1021, the equivalent unit rotor inertia time constant between described generating set i and described generating set j isCalculate the mechanical output function between described generating set i and described generating set j

Wherein M_iAnd M_jIt is respectively the inertia time constant of described generating set i and generating set j, P_MiAnd P_MjIt is respectively described sending out Group of motors i and generating set j are mechanically injected power, described generating set be mechanically injected power

Step S1022, calculates the electromagnetism active power between described generating set i and described generating set j

Wherein,A=Dcos (δ_j-θ_i)+Ecos(δ_i-θ_j), B=Dsin (δ_j-θ_i)-Esin(δ_i-θ_j), C =tan^-1(B/A),U_i、U_j、θ_i、θ_j、X_i、X_jRespectively represent generating set i and The generator terminal busbar voltage amplitude of generating set j, phase angle and equivalent internal reactance；E_i、E_jRepresent that generating set i's and j is virtual interior respectively Electromotive force；

Step S1023, calculates

9. method as claimed in claim 7 is it is characterised in that for any lead and lag unit pair in set omega, according to According to the unit after calculating to equivalent mechanical output P_MeqWith electromagnetism active-power P_EeqChange path, using curve fitting technique and from The angle of energy builds transient stability and estimates index, specifically includes step S201- step S203：

Step S201 is it is assumed that described generating set is to equivalent mechanical output P_MeqConstant in transient process, described generating set pair Equivalent electromagnetism active-power P_EeqChange path is fitted with following SIN function：Y=x₀sin(δ_ij+x₁)+x₂；

Wherein, x₀, x₁, x₂For matching SIN function coefficient to be asked for, y represents generating set to equivalent P_Eeq；

Step S202, the transient stability calculating between described generating set i and generating set j estimates index T_SEIij；

Step S203, the T of any lead and lag unit pair in set of computations Ω successively_SEIij, with the T of minimum unit pair_SEIij Transient stability as this time domain data estimates index S_TEI, S_TEI=min { T_SEIij,i,j∈Ω}；

Calculate the transient stability between described generating set i and generating set j in described step S202 and estimate index and include step S2021- step S2023：

Step S2021, according to P_EeqChange path fitting function y, tries to achieve the unstable equilibrium point δ after failure removal_iju, described not Stable equilibrium point is described generating set to equivalent mechanical output P_MeqWith electromagnetism active-power P_EeqThe intersection point of change path, δ_iju Computing formula is：

Step S2022, the transient state kinetic energy V of failure removal moment unit pair_TAFor：

$\begin{array}{c}{V}_{TA}={\&Integral;}_{{\mathrm{\&delta;}}_{ij0}}^{{\mathrm{\&delta;}}_{ijc}}\&lsqb;P_{Meq}-x_0\sin ({\mathrm{\&delta;}}_{ij}+x_1)-x_2\&rsqb;{\mathrm{d\&delta;}}_{ij}\\ =(P_{Meq}-x_2)({\mathrm{\&delta;}}_{ijc}-{\mathrm{\&delta;}}_{ij0})+x_0\&lsqb;\cos ({\mathrm{\&delta;}}_{ijc}+x_1)-cos({\mathrm{\&delta;}}_{ij0}+x_1)\&rsqb;\end{array}$

The critical potential V that after failure removal, system is absorbed_TBFor：

$\begin{array}{c}{V}_{TB}={\&Integral;}_{{\mathrm{\&delta;}}_{ijc}}^{{\mathrm{\&delta;}}_{iju}}\&lsqb;x_0\sin ({\mathrm{\&delta;}}_{ij}+x_1)+x_2-P_{Meq}\&rsqb;{\mathrm{d\&delta;}}_{ij}\\ =(x_2-P_{Meq})({\mathrm{\&delta;}}_{iju}-{\mathrm{\&delta;}}_{ijc})-x_0\&lsqb;\cos ({\mathrm{\&delta;}}_{iju}+x_1)-cos({\mathrm{\&delta;}}_{ijc}+x_1)\&rsqb;\end{array}$

Wherein, δ_ij0And δ_ijcThe unit in expression fault generation moment and failure removal moment is poor to generator rotor angle respectively；

Step S2023, the described transient stability calculating between described generating set i and described generating set j estimates index T_SEIij：

Wherein, T_SEIijIndex transient stability meaning is：T_SEIij>0 transient stability, T_SEIij=0 critical transient stability, T_SEIij<0, temporarily State unstability.

10. method as claimed in claim 9 keeps qualified it is characterised in that obtaining grid nodes voltage in described step S3 The method of rate includes：

Given transient process interior joint voltage holding capability investigates scope：Voltage Drop lower limit V_LWith persistent period V_C, count electrical network All nodes are in the interstitial content V of the horizontal tolerance interval of transient process fullness in the epigastrium and abdomen afc voltage_QN, calculate described transient process interior joint Voltage meets the qualification rate of holding capacity

The method obtaining the steady qualification rate of the heat of circuit described in electrical network in described step S3 includes：Given transient process interior joint is static Stability margin index and circuit heat steady investigation scope：The quiet steady nargin limit value S of node_IC, the steady limit value L of circuit heat_PRAnd duration T_C, Statistics electrical network meets the interstitial content S of quiet surely sustainable scope_QN, meet number of lines L of the surely sustainable scope of heat_QN, calculate institute State qualification rate S that transient process interior joint is quiet steady and circuit heat is steady_TR, S_TRComputing formula is：L_N、N_NRespectively Always get in touch with line number and nodes for system.

11. methods as claimed in claim 10 are it is characterised in that obtain electrical network entirety transient stability situation in described step S3 Aggregative indicator S_TI=α₂S_TEI+β₂V_TR+χ₂S_TR；

Wherein, α₂、β₂、χ₂It is respectively electrical network to press thermally-stabilised constraint, run the weight coefficient of electric parameters constraint and quiet steady constraint, root According to empirically or analytically needing flexible configuration, α₂、β₂、χ₂Meet：

$\left\{\begin{array}{c}0\&le;{\mathrm{\&alpha;}}_2\&le;1\\ 0\&le;{\mathrm{\&beta;}}_2\&le;1\\ 0\&le;{\mathrm{\&chi;}}_2\&le;1\\ 0\&le;{\mathrm{\&alpha;}}_2+{\mathrm{\&beta;}}_2+{\mathrm{\&chi;}}_2\&le;1\end{array}\right..$