CN101807798B - Section data integration method for power system safety and stability online analysis - Google Patents

Abstract

The invention belongs to the technical field of the power system and automation, which comprises the following steps: 1) calculating impact factor of a generator and a bus on safety and stability of the grid according to the margin and model of safety and stability of the previous operating state power grid; 2) for the power grid (internal grid) under the dispatching jurisdiction, adopting the weight number impact factor of the generator and the bus on safety and stability of the grid to correct weight number measured, performing state estimation to obtain the internal grid section data for security and stability online analysis; and 3) for the power grid (external grid) out of the dispatching jurisdiction, by modifying output power and load of an internal generator of the external gird, ensuring that power flow of a contact line between the internal grid and the external grid is basically constant with the power flow of the internal gird, and adjusted quantity of output power and load of the generator is calculated to be minimum after impact factor of the generator and the bus, and obtaining adjusted external gird section data; and 4) combining the data of the internal grid section and the external grid section data, arranging a balancing machine in the external grid so as to obtain the section data meeting the internal safety and stability online analysis.

Description

The section data integration method that is used for power system safety and stability online analysis

Technical field

The invention belongs to electric power system and technical field of automation.

Background technology

The operation of power networks profile data is the data basis of power system safety and stability online analysis, and its accuracy directly has influence on power system safety and stability online analysis result's practicality.For control centre in the interconnected network transfer pipe subnet (referred to as: Intranet) carry out safety and stability online analysis, usually directly adopt the state estimation result of this control centre's EMS (EMS) for the profile data of Intranet, for this control centre transfer the outer electrical network of scope tube (referred to as: profile data outer net) also is employing state estimation result usually.The data acquisition time that is used for Power Network Status Estimation is inconsistent, larger, the different weights of the weights degree of correlation of measurement amount in the result of calculation of state estimation and the residual equation, the result of state estimation is different, the weights of measurement amount are larger, and state estimation result is just near corresponding measuring value.The degree of correlation of the running status of each generator and bus is different in the security and stability of electrical network and the electrical network, if the estimated accuracy on the running status of the large generator of electricity net safety stable impact and bus is high, then the analysis result of electricity net safety stable will be more accurate., need to the state estimation algorithm of routine be improved the precision of analyzing to improve electricity net safety stable for this reason.Equally also need to be optimized adjustment to the influence degree of safety and stability according to generator and bus for the outer net profile data.

Therefore, the present invention proposes the section data integration method for power system safety and stability online analysis, optimize and revise at the state estimation algorithm of routine and outer net profile data and to introduce generator and bus in the algorithm to the factor of influence of safety and stability, realization can improve the precision of electricity net safety stable on-line analysis towards the section data integration of electricity net safety stable on-line analysis.

Summary of the invention

The present invention seeks to: realize the section data integration towards the electricity net safety stable on-line analysis, improve the precision of electricity net safety stable on-line analysis.

The present invention adopts following technical scheme to realize, comprises the steps:

1) generator, bus calculate the factor of influence of electricity net safety stable

Stable for microvariations, according to the microvariations stable mode analysis result of previous operation of power networks state, adopt formula (1) calculating generator to the stable factor of influence of electrical network microvariations.

$G_{s,i}=1+\underset{j=1}{\overset{m}{\mathrm{\Σ}}}\frac{p_{a.\mathrm{ij}}(\mathrm{\σ}-{\mathrm{\σ}}_j)}{\mathrm{\σ}}---\left(1\right)$

Wherein, σ is for being used for filtering out from all low frequency oscillation modes the damping ratio threshold value of relatively serious pattern, m be damping ratio less than the pattern count of σ, σ _jBe the damping ratio of j pattern, p _A.ijIt is the participation factors of generator i in j the pattern.

For static voltage stability, according to the static voltage stability pattern analysis result of previous operation of power networks state, adopt formula (2) to calculate bus to the factor of influence of static voltage stability.

$B_{s.i}=1+\underset{j=1}{\overset{n}{\mathrm{\Σ}}}\frac{p_{v.\mathrm{ij}}(\mathrm{\λ}-{\mathrm{\λ}}_j)}{\mathrm{\λ}}---\left(2\right)$

Wherein, λ is for being used for filtering out from all voltage stabilization patterns the characteristic value threshold value of relatively serious pattern, n be characteristic value less than the pattern count of λ, λ _jBe the characteristic value of j pattern, p _V.ijIt is the participation factors of j pattern median generatrix i.

For transient rotor angle stability and transient voltage safe and stable, transient rotor angle stability nargin, pattern and transient voltage safe and stable nargin, the pattern of each forecast failure of the previous operation of power networks state of foundation adopt respectively formula (3) and (4) calculating generator and bus to the factor of influence of electrical network transient safe and stable.

$G_{t.i}=1+\underset{j=1}{\overset{k}{\mathrm{\Σ}}}\frac{\left|q_{a.\mathrm{ij}}\right|({\mathrm{\η}}_a-{\mathrm{\η}}_{a.j})}{{\mathrm{\η}}_a}---\left(3\right)$

$B_{t.i}=1+\underset{j=1}{\overset{l}{\mathrm{\Σ}}}\frac{q_{v,\mathrm{ij}}({\mathrm{\η}}_v-{\mathrm{\η}}_{v,j})}{{\mathrm{\η}}_v}---\left(4\right)$

Wherein, η _aFor being used for filtering out from all forecast failures the stability margin threshold value of the relative catastrophe failure of transient rotor angle stability, k is that transient rotor angle stability nargin is less than η _aThe forecast failure number, η _A.jBe the transient rotor angle stability nargin (1≤η after the fault j generation _A.j≤ 1, η _A.j＜0 expression merit angle unstability, η _A.j=0 expression merit angle neutrality, η _A.j＞0 expression angle stability), q _A.ijBe the participation factors of generator i in the transient rotor angle stability pattern of fault j (participation factors of Critical Group generator for just, the participation factors of remaining mass-sending motor is for negative); η _vFor being used for filtering out from all forecast failures the security margin threshold value of the relative catastrophe failure of transient voltage safe and stable, l is that transient voltage safe and stable nargin is less than η _vThe forecast failure number, η _V.jTransient voltage safe and stable nargin (1≤η of rear bus i occurs for fault j _V.j≤ 1, η _V.jDangerous or the unstability of＜0 expression voltage, η _V.j=0 expression voltage critical safety and stability, η _V.j＞0 expression voltage security is stable), q _V.ijParticipation factors for the transient voltage safe and stable pattern median generatrix i of fault j.

2) be stabilized in the Intranet section state estimation of line analysis towards intranet security

State estimation algorithm to Intranet is made amendment, and adopts generator and bus the factor of influence of electricity net safety stable to be come the weights of correction measurement, namely in the residual equation of state estimation, for generator measurement amount weights, adopts (k _sG _S.i+ k _tG _T.i) and reflect that the product of the weights of error in measurement replaces, for bus measurement amount weights, adopt (l _sB _S.i+ l _tB _T.i) replace k wherein with the product of weights of reflection error in measurement _s, k _t, l _sAnd l _tIt is set point.Again Intranet is carried out state estimation and calculate, can obtain the Intranet profile data for safety and stability online analysis;

3) the outer net profile data that is stabilized in line analysis towards intranet security is optimized and revised

Take into account generator, bus to the factor of influence of electricity net safety stable in the profile data to outer net is optimized the target function of adjustment, formula (5) is the meritorious model of optimizing and revising of outer net profile data, and the idle work optimization adjustment model is similar with it.

Wherein, M is the adjustable generator number of outer net, Δ P _iBe the active power adjustment amount of generator, N is the deferrable load bus number of outer net, Δ P _jThe active power adjustment amount of load, C _KiThat i platform generator is to the active po wer sensitivity of k bar interconnection, C _KjThat j load is to the active po wer sensitivity of k bar interconnection, Δ P _L.kMeritorious poor of corresponding interconnection before to be that k bar Intranet interconnection is meritorious adjust with the outer net profile data, L is the number of internal, external network interconnection.Adopt Novel Algorithm to calculate gaining merit and the reactive power adjustment amount of generator and load, and outer net section primary data is made amendment, namely obtain being stabilized in towards intranet security the outer net profile data of line analysis.

4) the internal, external network profile data that is stabilized in line analysis towards intranet security merges

Step 2) and step 3) the internal, external network profile data that is stabilized in line analysis towards intranet security that obtains is combined, and will be arranged on outer net for the balancing machine that the whole network trend is calculated, then be met the whole network section tidal current data that intranet security is stabilized in line analysis.

Additional disclosure be: when the electricity net safety stable on-line analysis system starts for the first time, need to initial value be set to the factor of influence of generator and bus.

Effect and advantage:

The present invention is according to the different characteristic of the degree of correlation of the running status of each generator in the security and stability of electrical network and the electrical network and bus, by online calculating generator and the bus factor of influence to electricity net safety stable, optimize and revise at the state estimation algorithm of routine and outer net profile data and to introduce generator and bus in the algorithm to the factor of influence of safety and stability, realized towards the section data integration of electricity net safety stable on-line analysis, can improve the precision of electricity net safety stable on-line analysis, solve towards the difficult problem of the section data integration of large electricity net safety stable on-line analysis.

Description of drawings

Fig. 1 is the flow chart of the inventive method.

Embodiment

Below in conjunction with accompanying drawing 1, the inventive method is described in detail.

Step 1 is described among Fig. 1 is electricity net safety stable on-line analysis system when starting, owing to there not being previous electricity net safety stable analysis result in the continuous running, so generator and bus need to artificially be set to the factor of influence of electricity net safety stable;

Step 2 is described among Fig. 1 be calculating generator, bus to the factor of influence of electricity net safety stable, comprise following 3 parts:

1) stable for microvariations, according to the microvariations stable mode analysis result of previous operation of power networks state, adopt formula (1) calculating generator to the stable factor of influence of electrical network microvariations.

$G_{s.i}=1+\underset{j=1}{\overset{m}{\mathrm{\Σ}}}\frac{p_{a,\mathrm{ij}}(\mathrm{\σ}-{\mathrm{\σ}}_j)}{\mathrm{\σ}}---\left(1\right)$

Wherein, σ is for being used for filtering out from all low frequency oscillation modes the damping ratio threshold value of relatively serious pattern, m be damping ratio less than the pattern count of σ, σ _jBe the damping ratio of j pattern, p _A.ijIt is the participation factors of generator i in j the pattern.

2) for static voltage stability, according to the static voltage stability pattern analysis result of previous operation of power networks state, adopt formula (2) to calculate bus to the factor of influence of static voltage stability.

$B_{s.i}=1+\underset{j=1}{\overset{n}{\mathrm{\Σ}}}\frac{p_{v.\mathrm{ij}}(\mathrm{\λ}-{\mathrm{\λ}}_j)}{\mathrm{\λ}}---\left(2\right)$

Wherein, λ is for being used for filtering out from all voltage stabilization patterns the characteristic value threshold value of relatively serious pattern, n be characteristic value less than the pattern count of λ, λ _jBe the characteristic value of j pattern, p _V.ijIt is the participation factors of j pattern median generatrix i.

3) for transient rotor angle stability and transient voltage safe and stable, transient rotor angle stability nargin, pattern and transient voltage safe and stable nargin, the pattern of each forecast failure of the previous operation of power networks state of foundation adopt respectively formula (3) and (4) calculating generator and bus to the factor of influence of electrical network transient safe and stable.

$G_{t.i}=1+\underset{j=1}{\overset{k}{\mathrm{\Σ}}}\frac{\left|q_{a.\mathrm{ij}}\right|({\mathrm{\η}}_a-{\mathrm{\η}}_{a.j})}{{\mathrm{\η}}_a}---\left(3\right)$

$B_{t.i}=1+\underset{j=1}{\overset{l}{\mathrm{\Σ}}}\frac{q_{v,\mathrm{ij}}({\mathrm{\η}}_v-{\mathrm{\η}}_{v.j})}{{\mathrm{\η}}_v}---\left(4\right)$

Wherein, η _aFor being used for filtering out from all forecast failures the stability margin threshold value of the relative catastrophe failure of transient rotor angle stability, k is that transient rotor angle stability nargin is less than η _aThe forecast failure number, η _A.jBe the transient rotor angle stability nargin (1≤η after the fault j generation _A.j≤ 1, η _A.j＜0 expression merit angle unstability, η _A.j=0 expression merit angle neutrality, η _A.j＞0 expression angle stability), q _A.ijBe the participation factors of generator i in the transient rotor angle stability pattern of fault j (participation factors of Critical Group generator for just, the participation factors of remaining mass-sending motor is for negative); η _vFor being used for filtering out from all forecast failures the security margin threshold value of the relative catastrophe failure of transient voltage safe and stable, l is that transient voltage safe and stable nargin is less than η _vThe forecast failure number, η _V.jTransient voltage safe and stable nargin (1≤η of rear bus i occurs for fault j _V.j≤ 1, η _V.jDangerous or the unstability of＜0 expression voltage, η _V.j=0 expression voltage critical safety and stability, η _V.j＞0 expression voltage security is stable), q _V.ijParticipation factors for the transient voltage safe and stable pattern median generatrix i of fault j.

What step 3 was described among Fig. 1 is the state estimation of the Intranet section being carried out being stabilized in towards intranet security line analysis, namely the state estimation algorithm of Intranet is made amendment, adopt generator and bus the factor of influence of electricity net safety stable to be come the weights of correction measurement, namely in the residual equation of state estimation, for generator measurement amount weights, adopt (k _sG _S.i+ k _tG _T.i) and reflect that the product of the weights of error in measurement replaces, for bus measurement amount weights, adopt (l _sB _S.i+ l _tB _T.i) replace k wherein with the product of weights of reflection error in measurement _s, k _t, l _sAnd l _tIt is set point.Again Intranet is carried out state estimation and calculate, can obtain the Intranet profile data for safety and stability online analysis;

What step 4 was described among Fig. 1 is that the outer net section is carried out being stabilized in optimizing and revising of line analysis towards intranet security, namely in being optimized the target function of adjustment, the profile data to outer net takes into account generator, bus to the factor of influence of electricity net safety stable, formula (5) is the meritorious model of optimizing and revising of outer net profile data, and the idle work optimization adjustment model is similar with it.

Wherein, M is the adjustable generator number of outer net, Δ P _iBe the active power adjustment amount of generator, N is the deferrable load bus number of outer net, Δ P _jThe active power adjustment amount of load, C _KiThat i platform generator is to the active po wer sensitivity of k bar interconnection, C _KjThat j load is to the active po wer sensitivity of k bar interconnection, Δ P _L.kMeritorious poor of corresponding interconnection before to be that k bar Intranet interconnection is meritorious adjust with the outer net profile data, L is the number of internal, external network interconnection.Adopt Novel Algorithm to calculate gaining merit and the reactive power adjustment amount of generator and load, and outer net section primary data is made amendment, namely obtain being stabilized in towards intranet security the outer net profile data of line analysis.

Step 5 is described among Fig. 1 be internally, the net profile data carries out being stabilized in towards intranet security the merging of line analysis, namely step 3) and step 4) the internal, external network profile data that is stabilized in line analysis towards intranet security that obtains is combined, and will be arranged on outer net for the balancing machine that the whole network trend is calculated, then be met the whole network section tidal current data that intranet security is stabilized in line analysis.

Claims (5)

1. be used for the section data integration method of power system safety and stability online analysis, may further comprise the steps:

1) according to security margin and the pattern of previous operation of power networks state, calculates respectively generator and bus to the factor of influence of electricity net safety stable;

The electrical network of 2) administering for a certain control centre in the interconnected network, it is Intranet, adopt generator and bus the factor of influence of electricity net safety stable to be come the weights of correction measurement, carry out state estimation, obtain the Intranet profile data for safety and stability online analysis;

3) for the outer electrical network of this control centre's administration, it is outer net, by revising generator output and load in the outer net, the trend of guaranteeing the interconnection between the internal, external network is consistent with the Intranet trend, and it is minimum to take into account generator and bus generator output and the adjustment amount of load after to the factor of influence of electricity net safety stable, the outer net profile data after being adjusted;

4) the internal, external network profile data is combined, and will be arranged on outer net for the balancing machine that trend is calculated, then be met the section tidal current data that intranet security is stabilized in line analysis.

2. the section data integration method for power system safety and stability online analysis according to claim 1 is characterized in that step 1) in generator and bus the calculation procedure of the factor of influence of electricity net safety stable is comprised;

Stable for microvariations, according to the microvariations stable mode analysis result of previous operation of power networks state, adopt formula (1) calculating generator to the stable factor of influence of electrical network microvariations;

$G_{s.i}=1+\underset{j=1}{\overset{m}{\mathrm{\Σ}}}\frac{p_{a.\mathrm{ij}}(\mathrm{\σ}-{\mathrm{\σ}}_j)}{\mathrm{\σ}}---\left(1\right)$

Wherein, σ is for being used for filtering out from all low frequency oscillation modes the damping ratio threshold value of relatively serious pattern, m be damping ratio less than the pattern count of σ, σ _jBe the damping ratio of j pattern, p _A.ijIt is the participation factors of generator i in j the pattern;

For static voltage stability, according to the static voltage stability pattern analysis result of previous operation of power networks state, adopt formula (2) to calculate bus to the factor of influence of static voltage stability;

$B_{s.i}=1+\underset{j=1}{\overset{n}{\mathrm{\Σ}}}\frac{p_{v.\mathrm{ij}}(\mathrm{\λ}-{\mathrm{\λ}}_j)}{\mathrm{\λ}}---\left(2\right)$

Wherein, λ is for being used for filtering out from all voltage stabilization patterns the characteristic value threshold value of relatively serious pattern, n be characteristic value less than the pattern count of λ, λ _jBe the characteristic value of j pattern, p _V.ijIt is the participation factors of j pattern median generatrix i;

For transient rotor angle stability and transient voltage safe and stable, transient rotor angle stability nargin, pattern and transient voltage safe and stable nargin, the pattern of each forecast failure of the previous operation of power networks state of foundation adopt respectively formula (3) and (4) calculating generator and bus to the factor of influence of electrical network transient safe and stable;

$G_{t.i}=1+\underset{j=1}{\overset{k}{\mathrm{\Σ}}}\frac{\left|q_{a.\mathrm{ij}}\right|({\mathrm{\η}}_a-{\mathrm{\η}}_{a.j})}{{\mathrm{\η}}_a}---\left(3\right)$

$B_{t,i}=1+\underset{j=1}{\overset{l}{\mathrm{\Σ}}}\frac{q_{v,\mathrm{ij}}({\mathrm{\η}}_v-{\mathrm{\η}}_{v.j})}{{\mathrm{\η}}_v}---\left(4\right)$

Wherein, η _aFor being used for filtering out from all forecast failures the stability margin threshold value of the relative catastrophe failure of transient rotor angle stability, k is that transient rotor angle stability nargin is less than η _aThe forecast failure number, η _A.jBe the transient rotor angle stability nargin after the fault j generation, wherein ,-1≤η _A.j≤ 1, η _A.j＜0 expression merit angle unstability, η _A.j=0 expression merit angle neutrality, η _A.j＞0 expression angle stability, q _A.ijBe the participation factors of generator i in the transient rotor angle stability pattern of fault j, wherein the participation factors of Critical Group generator is for just, and the participation factors of remaining mass-sending motor is for negative; η _vFor being used for filtering out from all forecast failures the security margin threshold value of the relative catastrophe failure of transient voltage safe and stable, l is that transient voltage safe and stable nargin is less than η _vThe forecast failure number, η _V.jThe transient voltage safe and stable nargin of rear bus i occurs for fault j, wherein ,-1≤η _V.j≤ 1, η _V.jDangerous or the unstability of＜0 expression voltage, η _V.j=0 expression voltage critical safety and stability, η _V.j＞0 expression voltage security is stable, q _V.ijParticipation factors for the transient voltage safe and stable pattern median generatrix i of fault j.

3. the section data integration method for power system safety and stability online analysis according to claim 1, it is characterized in that, step 2) state estimation algorithm of Intranet is made amendment, adopt generator and bus the factor of influence of electricity net safety stable to be come the weights of correction measurement, namely in the residual equation of state estimation, for generator measurement amount weights, adopt (k _sG _S.i+ k _tG _T.i) and reflect that the product of the weights of error in measurement replaces, for bus measurement amount weights, adopt (l _sB _S.i+ l _tB _T.i) replace k wherein with the product of weights of reflection error in measurement _s, k _t, l _sAnd l _tIt is set point.

4. the section data integration method for power system safety and stability online analysis according to claim 1, it is characterized in that, step 3) take into account generator and bus to the factor of influence of electricity net safety stable in the profile data to outer net is optimized the target function of adjustment, formula (5) is the meritorious model of optimizing and revising of outer net profile data;

Wherein, M is the adjustable generator number of outer net, Δ P _iBe the active power adjustment amount of generator, N is the deferrable load bus number of outer net, Δ P _jThe active power adjustment amount of load, C _KiThat i platform generator is to the active po wer sensitivity of k bar interconnection, C _KjThat j load is to the active po wer sensitivity of k bar interconnection, Δ P _L.kMeritorious poor of corresponding interconnection before to be that k bar Intranet interconnection is meritorious adjust with the outer net profile data, L is the number of internal, external network interconnection.

5. the section data integration method for power system safety and stability online analysis according to claim 1, it is characterized in that, step 4) Intranet, outer net profile data are combined, and will be arranged on outer net for the balancing machine that the whole network trend is calculated, then be met the whole network section tidal current data that intranet security is stabilized in line analysis.

Images