CN104865843A - Power system hybrid simulation fault unified processing method - Google Patents

Abstract

The invention provides a power system hybrid simulation fault unified processing method. The method includes the following steps that: A, in hybrid simulation, a power system alternating current network is established at an electromechanical transient side through utilizing an electromechanical transient super real-time simulation program TH-STBLT; B, at a fault port, the power system alternating current network is divided into a normal sub network and multi-port fault sub networks; C, in the normal sub network, a power grid part is described through adopting a fundamental wave three-sequence network, and primary equipment dynamic elements such as a generator, a dynamic load and a power electronic device are connected into a positive sequence network; and D, in the multi-port fault sub networks, fault models of power system network elements are established according to fault information, and a multi-port and multi-fault complex admittance unified module Yf is derived; E, the Yf is superimposed into an admittance module Y of the positive sequence network of the normal sub network, and positive sequence voltage and positive sequence current of a fault port of the positive sequence network are calculated; and F, the voltage of the fault ports of a negative sequence network and a zero sequence network is calculated through the current of the fault port of the positive sequence network.

Description

The method of the unified process of a kind of electric system hybrid simulation fault

Technical field

The present invention is the method for the unified process of a kind of electric system hybrid simulation fault, belongs to alternating current-direct current bulk power grid digital simulation technique field.

Background technology

For the interactional electrical network characteristic of the extensive ac and dc systems of research and analysis, ensure bulk power grid safe and stable operation, accurate and believable analogue system is absolutely necessary instrument.At present, people have developed electromagnetism-dynamo-electric hybrid real-time simulation platform, realize organically combining electromagnetism in electric system and dynamo-electric transient state process, and close-coupled emulates together simultaneously.

To the AC network containing large-capacity power electronic installation, owing to there is complicated commutation course, the analysis of this kind of electrical network faces the problem similar to above-mentioned Ac/dc Power Systems and the demand of computational analysis.Along with a large amount of power electronic equipment put into operation, a large amount of power electronics power load drops into electrical network, electrical network asymmetrical three-phase situation is serious all the more.The practical problems occurred in these electrical networks proposes requirement to operation action under characteristic and the behavior under asymmetric fault, disturbance excitation of hybrid simulation simulation system, system asymmetrical three-phase operating mode.

Traditional hybrid simulation electromechanical transient program is to the process of asymmetrical three-phase operating mode, fault disturbance, mostly only calculate for single port electrical network short circuits or open conductors, lack unified disposal route, counting yield is low, is difficult in large scale electric network, process multiport and laterally, longitudinally answers fault.

Summary of the invention

For the shortcoming of prior art, the object of this invention is to provide the method for the unified process of a kind of electric system hybrid simulation fault, provide the complex admittance uniform expression that large scale electric network multiport is horizontal, longitudinally answer fault, cover all fault types of electric system comprehensively, overcome the unicity that conventional hybrid simulated fault calculates, thus improve the efficiency of hybrid simulation calculation of fault, realizing programming standardization, is a kind of fault calculation methods for transmission of practical.

To achieve these goals, the invention provides the method for the unified process of a kind of electric system hybrid simulation fault, the method comprises the following steps:

In A, hybrid simulation, electromechanical transient side adopts electromechanical transient faster than real time simulation program TH-STBLT to set up electric power system alternating current network;

B, at non-working port place, electric power system alternating current network is divided into normal-sub network and multiport Fault sub-network;

C, in normal-sub network, grid parts adopts first-harmonic three sequence network to describe, the primary equipment dynamic element access positive sequence nets such as generator, dynamic load, power electronic equipment;

D, in multiport Fault sub-network, set up the fault model of power system network element according to failure message, derive the complex admittance unified modules Y of the multiple fault of multiport _f; Wherein, described failure message comprises short circuits information and open conductors information; Described fault model comprises short circuits model and open conductors model;

The complex admittance module of short circuits is

$Y_f={\{Z_{\mathrm{ks}}-\left[\begin{array}{cc}{Z}_{\mathrm{km}}{Y}_{p2}& Z_{\mathrm{kn}}{Y}_{p0}\end{array}\right]\×{[I+\left[\begin{array}{cc}{Z}_{\mathrm{ks}}{Y}_{p2}& Z_{\mathrm{km}}{Y}_{p0}\\ Z_{\mathrm{kn}}{Y}_{p2}& (Z_{\mathrm{ks}}+3Z_{\mathrm{kg}})Y_{p0}\end{array}\right]]}^{-1}\left[\begin{array}{c}{Z}_{\mathrm{kn}}\\ Z_{\mathrm{km}}\end{array}\right]\}}^{-1}$

The complex admittance module of open conductors is

$Y_f={\{Z_{\mathrm{ks}}-\left[\begin{array}{cc}{Z}_{\mathrm{km}}{Y}_{p2}& Z_{\mathrm{kn}}{Y}_{p0}\end{array}\right]\×{[I+\left[\begin{array}{cc}{Z}_{\mathrm{ks}}{Y}_{p2}& Z_{\mathrm{km}}{Y}_{p0}\\ Z_{\mathrm{kn}}{Y}_{p2}& Z_{\mathrm{ks}}{Y}_{p0}\end{array}\right]]}^{-1}\left[\begin{array}{c}{Z}_{\mathrm{kn}}\\ Z_{\mathrm{km}}\end{array}\right]\}}^{-1}$

E, by Y _fsuperpose in normal-sub network positive sequence net admittance module Y, calculate positive sequence voltage and the forward-order current of positive sequence net non-working port;

F, calculated the non-working port voltage of negative phase-sequence and zero-sequence network by positive sequence port-current in fault.

The present invention gives the complex admittance uniform expression that large scale electric network multiport is horizontal, longitudinally answer fault, cover all fault types of electric system comprehensively, overcome the unicity that conventional hybrid simulated fault calculates, thus improve the efficiency of hybrid simulation calculation of fault, realizing programming standardization, is a kind of fault calculation methods for transmission of practical.

According to another embodiment of the present invention, short circuits model comprises ground short circuit fault, phase fault and three-phase imbalance load, short circuits model adopts a Y to connect three phase of impedance access non-working ports, and wherein Y is connect three phase of impedance neutral points and described by impedance earth.Multiport Fault sub-network one has t, i.e. Fault sub-network 1, Fault sub-network 2 ..., Fault sub-network t.

According to three-phase voltage and three-phase current relation, the phasor module relationship formula of t short circuits sub-network can be obtained:

$\left[\begin{array}{c}{\stackrel{\·}{U}}_{\mathrm{ka}}\\ {\stackrel{\·}{U}}_{\mathrm{kb}}\\ {\stackrel{\·}{U}}_{\mathrm{kc}}\end{array}\right]=\left[\begin{array}{ccc}{Z}_{\mathrm{ka}}+Z_{\mathrm{kg}}& Z_{\mathrm{kg}}& Z_{\mathrm{kg}}\\ Z_{\mathrm{kg}}& Z_{\mathrm{kb}}+Z_{\mathrm{kg}}& Z_{\mathrm{kg}}\\ Z_{\mathrm{kg}}& Z_{\mathrm{kg}}& Z_{\mathrm{kc}}+K_{\mathrm{kg}}\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{ka}}\\ {\stackrel{\·}{I}}_{\mathrm{kb}}\\ {\stackrel{\·}{I}}_{\mathrm{kc}}\end{array}\right]---\left(1\right)$

In formula (1), Z _ka, Z _kb, Z _kcbe respectively A phase, impedance module vector that B phase, C phase access non-working port, namely z _kgfor neutral point impedance module vector, namely be respectively the module vector of A phase, B phase, C phase voltage phasor, namely ${\stackrel{\·}{U}}_{\mathrm{ka}}={\left[\begin{array}{cccc}{\stackrel{\·}{U}}_{\mathrm{ka}1}& {\stackrel{\·}{U}}_{\mathrm{ka}2}& ......& {\stackrel{\·}{U}}_{\mathrm{kat}}\end{array}\right]}^T,{\stackrel{\·}{U}}_{\mathrm{kb}}={\left[\begin{array}{cccc}{\stackrel{\·}{U}}_{\mathrm{kb}1}& {\stackrel{\·}{U}}_{\mathrm{kb}2}& ......& {\stackrel{\·}{U}}_{\mathrm{kbt}}\end{array}\right]}^T,$ ${\stackrel{\·}{U}}_{\mathrm{kc}}={\left[\begin{array}{cccc}{\stackrel{\·}{U}}_{\mathrm{kc}1}& {\stackrel{\·}{U}}_{\mathrm{kc}2}& ......& {\stackrel{\·}{U}}_{\mathrm{kct}}\end{array}\right]}^T;$ be respectively the module vector of A phase, B phase, C phase current phasor, namely ${\stackrel{\·}{I}}_{\mathrm{ka}}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{ka}1}& {\stackrel{\·}{I}}_{\mathrm{ka}2}& ......& {\stackrel{\·}{I}}_{\mathrm{kat}}\end{array}\right]}^T,{\stackrel{\·}{I}}_{\mathrm{kb}}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{kb}1}& {\stackrel{\·}{I}}_{\mathrm{kb}2}& ......& {\stackrel{\·}{I}}_{\mathrm{kbt}}\end{array}\right]}^T,$ ${\stackrel{\·}{I}}_{\mathrm{kc}}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{kc}1}& {\stackrel{\·}{I}}_{\mathrm{kc}2}& ......& {\stackrel{\·}{I}}_{\mathrm{kct}}\end{array}\right]}^T.$

Adopt phase component method, voltage and current is carried out three-phase and the conversion of three sequences, then formula (1) can be deformed into

$\left[\begin{array}{c}{\stackrel{\·}{U}}_{\mathrm{ka}1}\\ {\stackrel{\·}{U}}_{\mathrm{ka}2}\\ {\stackrel{\·}{U}}_{\mathrm{ka}0}\end{array}\right]=\left[\begin{array}{ccc}{Z}_{\mathrm{ks}}& Z_{\mathrm{km}}& Z_{\mathrm{kn}}\\ Z_{\mathrm{kn}}& Z_{\mathrm{ks}}& Z_{\mathrm{km}}\\ Z_{\mathrm{km}}& Z_{\mathrm{kn}}& Z_{\mathrm{ks}}+{3K}_{\mathrm{kg}}\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{ka}1}\\ {\stackrel{\·}{I}}_{\mathrm{ka}2}\\ {\stackrel{\·}{I}}_{\mathrm{ka}0}\end{array}\right]---\left(2\right)$

Wherein be respectively the phasor module of the positive-sequence component of A phase voltage, negative sequence component and zero-sequence component, namely ${\stackrel{\·}{U}}_{\mathrm{ka}1}={\left[\begin{array}{cccc}{\stackrel{\·}{U}}_{\mathrm{ka}1-1}& {\stackrel{\·}{U}}_{\mathrm{ka}2-1}& ......& {\stackrel{\·}{U}}_{\mathrm{kat}-1}\end{array}\right]}^T,$ ${\stackrel{\·}{U}}_{\mathrm{ka}2}={\left[\begin{array}{cccc}{\stackrel{\·}{U}}_{\mathrm{ka}1-2}& {\stackrel{\·}{U}}_{\mathrm{ka}2-2}& ......& {\stackrel{\·}{U}}_{\mathrm{kat}-2}\end{array}\right]}^T,{\stackrel{\·}{U}}_{\mathrm{ka}0}={\left[\begin{array}{cccc}{\stackrel{\·}{U}}_{\mathrm{ka}1-0}& {\stackrel{\·}{U}}_{\mathrm{ka}2-0}& ......& {\stackrel{\·}{U}}_{\mathrm{kat}-0}\end{array}\right]}^T;$ be respectively the phasor module of the positive-sequence component of A phase current, negative sequence component and zero-sequence component, namely ${\stackrel{\·}{I}}_{\mathrm{ka}1}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{ka}1-1}& {\stackrel{\·}{I}}_{\mathrm{ka}2-1}& ......& {\stackrel{\·}{I}}_{\mathrm{kat}-1}\end{array}\right]}^T,{\stackrel{\·}{I}}_{\mathrm{ka}2}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{ka}1-2}& {\stackrel{\·}{I}}_{\mathrm{ka}2-2}& ......& {\stackrel{\·}{I}}_{\mathrm{kat}-2}\end{array}\right]}^T,$ ${\stackrel{\·}{I}}_{\mathrm{ka}0}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{ka}1-0}& {\stackrel{\·}{I}}_{\mathrm{ka}2-0}& ......& {\stackrel{\·}{I}}_{\mathrm{kat}-0}\end{array}\right]}^T;$ Z _km, Z _kn, Z _ksbe diagonal element, dimension equals fault tuple (i.e. t Fault sub-network), as shown in the formula:

$\left\{\begin{array}{c}{Z}_{\mathrm{km}}=\frac{1}{3}(Z_{\mathrm{ka}}+Z_{\mathrm{kb}}+Z_{\mathrm{kc}})\\ Z_m=\frac{1}{3}(Z_{\mathrm{ka}}+a^2{Z}_{\mathrm{kb}}+a{Z}_{\mathrm{kc}})\\ Z_{\mathrm{ks}}=\frac{1}{3}(Z_{\mathrm{ka}}+a{Z}_{\mathrm{kb}}+a^2{Z}_{\mathrm{kc}})\end{array}\right.---\left(3\right)$

Wherein a is for quoting operator e is the end of natural logarithm, and j is imaginary unit;

At short circuits port, fault branch electric current and the non-working port voltage relationship of grid side Negative-sequence Net and zero sequence net are

${\stackrel{\·}{I}}_{\mathrm{ka}2}=-Y_{p2}{\stackrel{\·}{U}}_{\mathrm{ka}2}---\left(4\right)$

${\stackrel{\·}{I}}_{\mathrm{ka}0}=-Y_{p0}{\stackrel{\·}{U}}_{\mathrm{ka}0}---\left(5\right)$

Wherein Y _p2and Y _p0be respectively negative phase-sequence admittance module and the zero sequence admittance module of normal-sub network.

Formula (4) and (5) are substituted into (2), obtains

${\stackrel{\·}{U}}_{\mathrm{ka}1}=Z_{\mathrm{ks}}{\stackrel{\·}{I}}_{\mathrm{ka}1}-\left[\begin{array}{cc}{Z}_{\mathrm{km}}{Y}_{p2}& Z_{\mathrm{kn}}{Y}_{p0}\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{U}}_{\mathrm{ka}2}\\ {\stackrel{\·}{U}}_{\mathrm{ka}0}\end{array}\right]---\left(6\right)$

$\left[\begin{array}{c}{\stackrel{\·}{U}}_{\mathrm{ka}2}\\ {\stackrel{\·}{U}}_{\mathrm{ka}0}\end{array}\right]={[I+\left[\begin{array}{cc}{Z}_{\mathrm{ks}}{Y}_{p2}& Z_{\mathrm{km}}{Y}_{p0}\\ Z_{\mathrm{kn}}{Y}_{p2}& (Z_{\mathrm{ks}}+3Z_{\mathrm{kg}})Y_{p0}\end{array}\right]]}^{-1}\left[\begin{array}{c}{Z}_{\mathrm{kn}}\\ Z_{\mathrm{km}}\end{array}\right]{\stackrel{\·}{I}}_{\mathrm{ka}1}---\left(7\right)$

Formula (7) is substituted into (6), obtains

${\stackrel{\·}{I}}_{\mathrm{ka}1}={\{Z_{\mathrm{ks}}-\left[\begin{array}{cc}{Z}_{\mathrm{km}}{Y}_{p2}& Z_{\mathrm{kn}}{Y}_{p0}\end{array}\right]\×{[I+\left[\begin{array}{cc}{Z}_{\mathrm{ks}}{Y}_{p2}& Z_{\mathrm{km}}{Y}_{p0}\\ Z_{\mathrm{kn}}{Y}_{p2}& (Z_{\mathrm{ks}}+3Z_{\mathrm{kg}})Y_{p0}\end{array}\right]]}^{-1}\left[\begin{array}{c}{Z}_{\mathrm{kn}}\\ Z_{\mathrm{km}}\end{array}\right]\}}^{-1}{\stackrel{\·}{U}}_{\mathrm{ka}1}---\left(8\right)$

Wherein I is unit module,

The complex admittance module that therefore can obtain short circuits is

$Y_f={\{Z_{\mathrm{ks}}-\left[\begin{array}{cc}{Z}_{\mathrm{km}}{Y}_{p2}& Z_{\mathrm{kn}}{Y}_{p0}\end{array}\right]\×{[I+\left[\begin{array}{cc}{Z}_{\mathrm{ks}}{Y}_{p2}& Z_{\mathrm{km}}{Y}_{p0}\\ Z_{\mathrm{kn}}{Y}_{p2}& (Z_{\mathrm{ks}}+3Z_{\mathrm{kg}})Y_{p0}\end{array}\right]]}^{-1}\left[\begin{array}{c}{Z}_{\mathrm{kn}}\\ Z_{\mathrm{km}}\end{array}\right]\}}^{-1}---\left(9\right)$

According to another embodiment of the present invention, open conductors model comprises broken string, the fault branch that open conductors model increases can access faulty line by three phase of impedance and describe: under normal circumstances, fault branch resistance value is a very little value, and three phase of impedance arrange different values to be combined into different open conductors or disturbance.Multiport Fault sub-network one has r, i.e. Fault sub-network 1, Fault sub-network 2 ..., Fault sub-network r.

According to three-phase voltage and three-phase current relation, the phasor module relationship formula of r open conductors sub-network can be obtained:

$\left[\begin{array}{c}\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}}\\ \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kb}}\\ \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kc}}\end{array}\right]=\left[\begin{array}{ccc}{Z}_{\mathrm{ka}}& 0& 0\\ 0& Z_{\mathrm{kb}}& 0\\ 0& 0& Z_{\mathrm{kc}}\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{ka}}\\ {\stackrel{\·}{I}}_{\mathrm{kb}}\\ {\stackrel{\·}{I}}_{\mathrm{kc}}\end{array}\right]---\left(10\right)$

In formula (10), Z _ka, Z _kb, Z _kcbe respectively A phase, impedance module vector that B phase, C phase access non-working port, namely be respectively A phase, B phase, C phase fault port Impedance two ends voltage phasor module vector, namely $\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}}={\left[\begin{array}{cccc}\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}1}& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}2}& ......& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kar}}\end{array}\right]}^T,$ $\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kb}}={\left[\begin{array}{cccc}\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kb}1}& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kb}2}& ......& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kbr}}\end{array}\right]}^T,\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kc}}={\left[\begin{array}{cccc}{\mathrm{\Δ}\stackrel{\·}{U}}_{\mathrm{kc}1}& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kc}2}& ......& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kcr}}\end{array}\right]}^T;,$ be respectively flow through A phase, B phase, C phase of impedance electric current phasor module vector, namely ${\stackrel{\·}{I}}_{\mathrm{ka}}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{ka}1}& {\stackrel{\·}{I}}_{\mathrm{ka}2}& ......& {\stackrel{\·}{I}}_{\mathrm{kar}}\end{array}\right]}^T,{\stackrel{\·}{I}}_{\mathrm{kb}}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{kb}1}& {\stackrel{\·}{I}}_{\mathrm{kb}2}& ......& {\stackrel{\·}{I}}_{\mathrm{kbr}}\end{array}\right]}^T,$ ${\stackrel{\·}{I}}_{\mathrm{kc}}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{kc}1}& {\stackrel{\·}{I}}_{\mathrm{kc}2}& ......& {\stackrel{\·}{I}}_{\mathrm{kcr}}\end{array}\right]}^T.$

Adopt phase component method, voltage and current is carried out three-phase and the conversion of three sequences, then formula (10) can be deformed into

$\left[\begin{array}{c}\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}1}\\ \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}2}\\ \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}0}\end{array}\right]=\left[\begin{array}{ccc}{Z}_{\mathrm{ks}}& Z_{\mathrm{km}}& Z_{\mathrm{kn}}\\ Z_{\mathrm{kn}}& Z_{\mathrm{ks}}& Z_{\mathrm{km}}\\ Z_{\mathrm{km}}& Z_{\mathrm{kn}}& Z_{\mathrm{ks}}\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{ka}1}\\ {\stackrel{\·}{I}}_{\mathrm{ka}2}\\ {\stackrel{\·}{I}}_{\mathrm{ka}0}\end{array}\right]---\left(11\right)$

Wherein be respectively the phasor module of the positive-sequence component of A phase fault port Impedance both end voltage, negative sequence component and zero-sequence component, namely $\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}1}={\left[\begin{array}{cccc}\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}1-1}& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}2-1}& ......& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kat}-1}\end{array}\right]}^T,$ $\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}2}={\left[\begin{array}{cccc}\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}1-2}& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}2-2}& ......& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kat}-2}\end{array}\right]}^T,$ $\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}0}={\left[\begin{array}{cccc}\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}1-0}& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}2-0}& ......& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kat}-0}\end{array}\right]}^T;$ be respectively the phasor module of electric current positive-sequence component, negative sequence component and the zero-sequence component flowing through A phase of impedance, namely ${\stackrel{\·}{I}}_{\mathrm{ka}1}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{ka}1-1}& {\stackrel{\·}{I}}_{\mathrm{ka}2-1}& ......& {\stackrel{\·}{I}}_{\mathrm{kat}-1}\end{array}\right]}^T,{\stackrel{\·}{I}}_{\mathrm{ka}2}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{ka}1-2}& {\stackrel{\·}{I}}_{\mathrm{ka}2-2}& ......& {\stackrel{\·}{I}}_{\mathrm{kat}-2}\end{array}\right]}^T,$ ${\stackrel{\·}{I}}_{\mathrm{ka}0}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{ka}1-0}& {\stackrel{\·}{I}}_{\mathrm{ka}2-0}& ......& {\stackrel{\·}{I}}_{\mathrm{kat}-0}\end{array}\right]}^T;$ Z _km, Z _kn, Z _ksbe diagonal element, dimension equals fault tuple (i.e. r Fault sub-network), as shown in the formula:

$\left\{\begin{array}{c}{Z}_{\mathrm{km}}=\frac{1}{3}(Z_{\mathrm{ka}}+Z_{\mathrm{kb}}+Z_{\mathrm{kc}})\\ Z_{\mathrm{kn}}=\frac{1}{3}(Z_{\mathrm{ka}}+a^2{Z}_{\mathrm{kb}}+a{Z}_{\mathrm{kc}})\\ Z_{\mathrm{ks}}=\frac{1}{3}(Z_{\mathrm{ka}}+a{Z}_{\mathrm{kb}}+a^2{Z}_{\mathrm{kc}})\end{array}\right.---\left(12\right)$

Wherein a is for quoting operator e is the end of natural logarithm, _jfor imaginary unit;

At open conductors port, fault branch electric current and the fault branch voltage relationship of grid side Negative-sequence Net and zero sequence net are

${\stackrel{\·}{I}}_{\mathrm{ka}2}=-Y_{p2}\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}2}---\left(13\right)$

${\stackrel{\·}{I}}_{\mathrm{ka}0}=-Y_{p0}\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}0}---\left(14\right)$

Wherein Y _p2and Y _p0be respectively negative phase-sequence admittance module and the zero sequence admittance module of normal-sub network.Simultaneous formula (11)-(14), derivation obtains

$\left[\begin{array}{c}\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}2}\\ \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}0}\end{array}\right]={[I+\left[\begin{array}{cc}{Z}_{\mathrm{ks}}{Y}_{p2}& Z_{\mathrm{km}}{Y}_{p0}\\ Z_{\mathrm{kn}}{Y}_{p2}& Z_{\mathrm{ks}}{Y}_{p0}\end{array}\right]]}^{-1}\left[\begin{array}{c}{Z}_{\mathrm{kn}}\\ Z_{\mathrm{km}}\end{array}\right]{\stackrel{\·}{I}}_{\mathrm{ka}1}---\left(15\right)$

Then the complex admittance module of open conductors is

$Y_f={\{Z_{\mathrm{ks}}-\left[\begin{array}{cc}{Z}_{\mathrm{km}}{Y}_{p2}& Z_{\mathrm{kn}}{Y}_{p0}\end{array}\right]\×{[I+\left[\begin{array}{cc}{Z}_{\mathrm{ks}}{Y}_{p2}& Z_{\mathrm{km}}{Y}_{p0}\\ Z_{\mathrm{kn}}{Y}_{p2}& Z_{\mathrm{ks}}{Y}_{p0}\end{array}\right]]}^{-1}\left[\begin{array}{c}{Z}_{\mathrm{kn}}\\ Z_{\mathrm{km}}\end{array}\right]\}}^{-1}---\left(16\right)$

Wherein I is unit module,

According to another embodiment of the present invention, adopt phase component method, voltage and current is carried out three-phase and the conversion of three sequences.

According to another embodiment of the present invention, in the short circuits of electric system hybrid simulation calculates, by the complex admittance module Y of short circuits _fmerge and enter in the admittance module Y of electrical network positive sequence net.

According to another embodiment of the present invention, in the open conductors of electric system hybrid simulation calculates, by the complex admittance module Y of open conductors _fmerge and enter in electrical network positive sequence net admittance module Y.

Compared with prior art, the present invention possesses following beneficial effect:

Method of the present invention gives the complex admittance uniform expression of the multiple fault of bulk power grid multiport, both the calculating to multiport short circuits had been contained, also the calculating to multiport open conductors is comprised, adapt to the simulation of the various complex fault of analysis and research, thus improve hybrid simulation calculation of fault efficiency and practicality, be a kind of effective, feasible solution.Meanwhile, the method is suitable for real-time simulated animation.

Below in conjunction with accompanying drawing, the present invention is described in further detail.

Accompanying drawing explanation

Fig. 1 is that the electric power system fault of embodiment 1 unifies handling principle figure;

Fig. 2 is the electric system short circuits schematic diagram of embodiment 1;

Fig. 3 is the electric system open conductors schematic diagram of embodiment 1;

Embodiment

Embodiment 1

That the electric system hybrid simulation fault of the present embodiment unifies handling principle figure shown in Fig. 1.The method of the unified process of electric system hybrid simulation fault of the present embodiment comprises the following steps:

In A, hybrid simulation, electromechanical transient side adopts electromechanical transient faster than real time simulation program TH-STBLT to set up electric power system alternating current network;

B, at non-working port place, electric power system alternating current network is divided into normal-sub network and multiport Fault sub-network;

C, in normal-sub network, grid parts adopts first-harmonic three sequence network to describe, the primary equipment dynamic element access positive sequence nets such as generator, dynamic load, power electronic equipment;

D, in multiport Fault sub-network, set up the fault model of power system network element according to failure message, derive the complex admittance unified modules Y of the multiple fault of multiport _f; Wherein, described failure message comprises short circuits information and open conductors information; Described fault model comprises short circuits model and open conductors model;

The complex admittance module of short circuits is

$Y_f={\{Z_{\mathrm{ks}}-\left[\begin{array}{cc}{Z}_{\mathrm{km}}{Y}_{p2}& Z_{\mathrm{kn}}{Y}_{p0}\end{array}\right]\×{[I+\left[\begin{array}{cc}{Z}_{\mathrm{ks}}{Y}_{p2}& Z_{\mathrm{km}}{Y}_{p0}\\ Z_{\mathrm{kn}}{Y}_{p2}& (Z_{\mathrm{ks}}+3Z_{\mathrm{kg}})Y_{p0}\end{array}\right]]}^{-1}\left[\begin{array}{c}{Z}_{\mathrm{kn}}\\ Z_{\mathrm{km}}\end{array}\right]\}}^{-1}$

The complex admittance module of open conductors is

$Y_f={\{Z_{\mathrm{ks}}-\left[\begin{array}{cc}{Z}_{\mathrm{km}}{Y}_{p2}& Z_{\mathrm{kn}}{Y}_{p0}\end{array}\right]\×{[I+\left[\begin{array}{cc}{Z}_{\mathrm{ks}}{Y}_{p2}& Z_{\mathrm{km}}{Y}_{p0}\\ Z_{\mathrm{kn}}{Y}_{p2}& Z_{\mathrm{ks}}{Y}_{p0}\end{array}\right]]}^{-1}\left[\begin{array}{c}{Z}_{\mathrm{kn}}\\ Z_{\mathrm{km}}\end{array}\right]\}}^{-1}$

E, by Y _fsuperpose in normal-sub network positive sequence net admittance module Y, calculate positive sequence voltage and the forward-order current of positive sequence net non-working port;

F, calculated the non-working port voltage of negative phase-sequence and zero-sequence network by positive sequence port-current in fault.

Above-mentioned short circuits model, comprise ground short circuit fault, phase fault and three-phase imbalance load etc., short circuits model adopts a Y to connect three phase of impedance access non-working ports, wherein Y is connect three phase of impedance neutral points and is described by impedance earth, and multiport Fault sub-network one has t, i.e. Fault sub-network 1, Fault sub-network 2,, Fault sub-network t, as shown in Figure 2.

According to three-phase voltage and three-phase current relation, the phasor module relationship formula of t Fault sub-network can be obtained:

$\left[\begin{array}{c}{\stackrel{\·}{U}}_{\mathrm{ka}}\\ {\stackrel{\·}{U}}_{\mathrm{kb}}\\ {\stackrel{\·}{U}}_{\mathrm{kc}}\end{array}\right]=\left[\begin{array}{ccc}{Z}_{\mathrm{ka}}+Z_{\mathrm{kg}}& Z_{\mathrm{kg}}& Z_{\mathrm{kg}}\\ Z_{\mathrm{kg}}& Z_{\mathrm{kb}}+Z_{\mathrm{kg}}& Z_{\mathrm{kg}}\\ Z_{\mathrm{kg}}& Z_{\mathrm{kg}}& Z_{\mathrm{kc}}+K_{\mathrm{kg}}\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{ka}}\\ {\stackrel{\·}{I}}_{\mathrm{kb}}\\ {\stackrel{\·}{I}}_{\mathrm{kc}}\end{array}\right]---\left(1\right)$

In formula (1), Z _ka, Z _kb, Z _kcbe respectively A phase, impedance module vector that B phase, C phase access non-working port, namely z _kgfor neutral point impedance module vector, namely be respectively the module vector of A phase, B phase, C phase voltage phasor, namely ${\stackrel{\·}{U}}_{\mathrm{ka}}={\left[\begin{array}{cccc}{\stackrel{\·}{U}}_{\mathrm{ka}1}& {\stackrel{\·}{U}}_{\mathrm{ka}2}& ......& {\stackrel{\·}{U}}_{\mathrm{kat}}\end{array}\right]}^T,$ ${\stackrel{\·}{U}}_{\mathrm{kb}}={\left[\begin{array}{cccc}{\stackrel{\·}{U}}_{\mathrm{kb}1}& {\stackrel{\·}{U}}_{\mathrm{kb}2}& ......& {\stackrel{\·}{U}}_{\mathrm{kbt}}\end{array}\right]}^T,{\stackrel{\·}{U}}_{\mathrm{kc}}={\left[\begin{array}{cccc}{\stackrel{\·}{U}}_{\mathrm{kc}1}& {\stackrel{\·}{U}}_{\mathrm{kc}2}& ......& {\stackrel{\·}{U}}_{\mathrm{kct}}\end{array}\right]}^T;$ be respectively the module vector of A phase, B phase, C phase current phasor, namely ${\stackrel{\·}{I}}_{\mathrm{ka}}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{ka}1}& {\stackrel{\·}{I}}_{\mathrm{ka}2}& ......& {\stackrel{\·}{I}}_{\mathrm{kat}}\end{array}\right]}^T,$ ${\stackrel{\·}{I}}_{\mathrm{kb}}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{kb}1}& {\stackrel{\·}{I}}_{\mathrm{kb}2}& ......& {\stackrel{\·}{I}}_{\mathrm{kbt}}\end{array}\right]}^T,{\stackrel{\·}{I}}_{\mathrm{kc}}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{kc}1}& {\stackrel{\·}{I}}_{\mathrm{kc}2}& ......& {\stackrel{\·}{I}}_{\mathrm{kct}}\end{array}\right]}^T.$ Lower same.

Phase component method is adopted to obtain,

$\left[\begin{array}{c}{\stackrel{\·}{U}}_{\mathrm{ka}}\\ {\stackrel{\·}{U}}_{\mathrm{kb}}\\ {\stackrel{\·}{U}}_{\mathrm{kc}}\end{array}\right]=\left[\begin{array}{ccc}I& I& I\\ a^2& a& I\\ a& a^2& I\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{U}}_{\mathrm{ka}1}\\ {\stackrel{\·}{U}}_{\mathrm{ka}2}\\ {\stackrel{\·}{U}}_{\mathrm{ka}0}\end{array}\right]---\left(2\right)$

$\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{ka}}\\ {\stackrel{\·}{I}}_{\mathrm{kb}}\\ {\stackrel{\·}{I}}_{\mathrm{kc}}\end{array}\right]=\left[\begin{array}{ccc}I& I& I\\ a^2& a& I\\ a& a^2& I\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{ka}1}\\ {\stackrel{\·}{I}}_{\mathrm{ka}2}\\ {\stackrel{\·}{I}}_{\mathrm{ka}0}\end{array}\right]---\left(3\right)$

Wherein I is unit module, a is for quoting operator e is the end of natural logarithm, and j is imaginary unit; be respectively the phasor module of the positive-sequence component of A phase voltage, negative sequence component and zero-sequence component, namely ${\stackrel{\·}{U}}_{\mathrm{ka}1}={\left[\begin{array}{cccc}{\stackrel{\·}{U}}_{\mathrm{ka}1-1}& {\stackrel{\·}{U}}_{\mathrm{ka}2-1}& ......& {\stackrel{\·}{U}}_{\mathrm{kat}-1}\end{array}\right]}^T{,\stackrel{\·}{U}}_{\mathrm{ka}2}={\left[\begin{array}{cccc}{\stackrel{\·}{U}}_{\mathrm{ka}1-2}& {\stackrel{\·}{U}}_{\mathrm{ka}2-2}& ......& {\stackrel{\·}{U}}_{\mathrm{kat}-2}\end{array}\right]}^T,$ ${\stackrel{\·}{U}}_{\mathrm{ka}0}={\left[\begin{array}{cccc}{\stackrel{\·}{U}}_{\mathrm{ka}1-0}& {\stackrel{\·}{U}}_{\mathrm{ka}2-0}& ......& {\stackrel{\·}{U}}_{\mathrm{kat}-0}\end{array}\right]}^T;$ be respectively the phasor module of the positive-sequence component of A phase current, negative sequence component and zero-sequence component, namely ${\stackrel{\·}{I}}_{\mathrm{ka}1}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{ka}1-1}& {\stackrel{\·}{I}}_{\mathrm{ka}2-1}& ......& {\stackrel{\·}{I}}_{\mathrm{kat}-1}\end{array}\right]}^T,$ ${\stackrel{\·}{I}}_{\mathrm{ka}2}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{ka}1-2}& {\stackrel{\·}{I}}_{\mathrm{ka}2-2}& ......& {\stackrel{\·}{I}}_{\mathrm{kat}-2}\end{array}\right]}^T,{\stackrel{\·}{I}}_{\mathrm{ka}0}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{ka}1-0}& {\stackrel{\·}{I}}_{\mathrm{ka}2-0}& ......& {\stackrel{\·}{I}}_{\mathrm{kat}-0}\end{array}\right]}^T.$ Lower same.

Formula (2) and (3) are substituted into formula (1), obtains after arrangement

$\left[\begin{array}{c}{\stackrel{\·}{U}}_{\mathrm{ka}1}\\ {\stackrel{\·}{U}}_{\mathrm{ka}2}\\ {\stackrel{\·}{U}}_{\mathrm{ka}0}\end{array}\right]=\left[\begin{array}{ccc}{Z}_{\mathrm{ks}}& Z_{\mathrm{km}}& Z_{\mathrm{kn}}\\ Z_{\mathrm{kn}}& Z_{\mathrm{ks}}& Z_{\mathrm{km}}\\ Z_{\mathrm{km}}& Z_{\mathrm{kn}}& Z_{\mathrm{ks}}+{3K}_{\mathrm{kg}}\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{ka}1}\\ {\stackrel{\·}{I}}_{\mathrm{ka}2}\\ {\stackrel{\·}{I}}_{\mathrm{ka}0}\end{array}\right]---\left(4\right)$

The block elements of formula (4) as shown in the formula

$\left\{\begin{array}{c}{Z}_{\mathrm{km}}=\frac{1}{3}(Z_{\mathrm{ka}}+Z_{\mathrm{kb}}+Z_{\mathrm{kc}})\\ Z_m=\frac{1}{3}(Z_{\mathrm{ka}}+a^2{Z}_{\mathrm{kb}}+a{Z}_{\mathrm{kc}})\\ Z_{\mathrm{ks}}=\frac{1}{3}(Z_{\mathrm{ka}}+a{Z}_{\mathrm{kb}}+a^2{Z}_{\mathrm{kc}})\end{array}\right.---\left(5\right)$

Wherein Z _km, Z _kn, Z _ksbe diagonal element, dimension equals fault tuple (i.e. t Fault sub-network).

At non-working port, fault branch electric current and the non-working port voltage relationship of grid side Negative-sequence Net and zero sequence net are

${\stackrel{\·}{I}}_{\mathrm{ka}2}=-Y_{p2}{\stackrel{\·}{U}}_{\mathrm{ka}2}---\left(6\right)$

${\stackrel{\·}{I}}_{\mathrm{ka}0}=-Y_{p0}{\stackrel{\·}{U}}_{\mathrm{ka}0}---\left(7\right)$

Wherein Y _p2and Y _p0be respectively negative phase-sequence admittance module and the zero sequence admittance module of normal-sub network.

Formula (6) and (7) are substituted into (4), obtains

${\stackrel{\·}{U}}_{\mathrm{ka}1}=Z_{\mathrm{ks}}{\stackrel{\·}{I}}_{\mathrm{ka}1}-\left[\begin{array}{cc}{Z}_{\mathrm{km}}{Y}_{p2}& Z_{\mathrm{kn}}{Y}_{p0}\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{U}}_{\mathrm{ka}2}\\ {\stackrel{\·}{U}}_{\mathrm{ka}0}\end{array}\right]---\left(8\right)$

$\left[\begin{array}{c}{\stackrel{\·}{U}}_{\mathrm{ka}2}\\ {\stackrel{\·}{U}}_{\mathrm{ka}0}\end{array}\right]={[I+\left[\begin{array}{cc}{Z}_{\mathrm{ks}}{Y}_{p2}& Z_{\mathrm{km}}{Y}_{p0}\\ Z_{\mathrm{kn}}{Y}_{p2}& (Z_{\mathrm{ks}}+3Z_{\mathrm{kg}})Y_{p0}\end{array}\right]]}^{-1}\left[\begin{array}{c}{Z}_{\mathrm{kn}}\\ Z_{\mathrm{km}}\end{array}\right]{\stackrel{\·}{I}}_{\mathrm{ka}1}---\left(9\right)$

Formula (9) is substituted into (8), obtains

${\stackrel{\·}{I}}_{\mathrm{ka}1}={\{Z_{\mathrm{ks}}-\left[\begin{array}{cc}{Z}_{\mathrm{km}}{Y}_{p2}& Z_{\mathrm{kn}}{Y}_{p0}\end{array}\right]\×{[I+\left[\begin{array}{cc}{Z}_{\mathrm{ks}}{Y}_{p2}& Z_{\mathrm{km}}{Y}_{p0}\\ Z_{\mathrm{kn}}{Y}_{p2}& (Z_{\mathrm{ks}}+3Z_{\mathrm{kg}})Y_{p0}\end{array}\right]]}^{-1}\left[\begin{array}{c}{Z}_{\mathrm{kn}}\\ Z_{\mathrm{km}}\end{array}\right]\}}^{-1}{\stackrel{\·}{U}}_{\mathrm{ka}1}---\left(10\right)$

Therefore the complex admittance module of short circuits can be obtained

$Y_f={\{Z_{\mathrm{ks}}-\left[\begin{array}{cc}{Z}_{\mathrm{km}}{Y}_{p2}& Z_{\mathrm{kn}}{Y}_{p0}\end{array}\right]\×{[I+\left[\begin{array}{cc}{Z}_{\mathrm{ks}}{Y}_{p2}& Z_{\mathrm{km}}{Y}_{p0}\\ Z_{\mathrm{kn}}{Y}_{p2}& (Z_{\mathrm{ks}}+3Z_{\mathrm{kg}})Y_{p0}\end{array}\right]]}^{-1}\left[\begin{array}{c}{Z}_{\mathrm{kn}}\\ Z_{\mathrm{km}}\end{array}\right]\}}^{-1}---\left(11\right)$

In the short circuits of electric system hybrid simulation calculates, by Y _fmerge and enter in the admittance module Y of electrical network positive sequence net.

Above-mentioned open conductors model, comprises broken string etc.The fault branch that open conductors model increases can access faulty lines by three phase of impedance and describe: under normal circumstances, fault branch resistance value is a very little value, and three phase of impedance arrange different values and can be combined into different open conductors or disturbance.Multiport Fault sub-network one has r, i.e. Fault sub-network 1, Fault sub-network 2 ..., Fault sub-network r, as shown in Figure 3.

According to three-phase voltage and three-phase current relation, the phasor module relationship formula of r Fault sub-network can be obtained:

$\left[\begin{array}{c}\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}}\\ \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kb}}\\ \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kc}}\end{array}\right]=\left[\begin{array}{ccc}{Z}_{\mathrm{ka}}& 0& 0\\ 0& Z_{\mathrm{kb}}& 0\\ 0& 0& Z_{\mathrm{kc}}\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{ka}}\\ {\stackrel{\·}{I}}_{\mathrm{kb}}\\ {\stackrel{\·}{I}}_{\mathrm{kc}}\end{array}\right]---\left(12\right)$

In formula (12), Z _ka, Z _kb, Z _kcbe respectively A phase, impedance module vector that B phase, C phase access non-working port, namely be respectively A phase, B phase, C phase fault port Impedance two ends voltage phasor module vector, namely $\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}}={\left[\begin{array}{cccc}\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}1}& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}2}& ......& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kar}}\end{array}\right]}^T,\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kb}}={\left[\begin{array}{cccc}{\mathrm{\Δ}\stackrel{\·}{U}}_{\mathrm{kb}1}& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kb}2}& ......& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kbr}}\end{array}\right]}^T,$ $\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kc}}={\left[\begin{array}{cccc}{\mathrm{\Δ}\stackrel{\·}{U}}_{\mathrm{kc}1}& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kc}2}& ......& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kcr}}\end{array}\right]}^T;,$ be respectively flow through A phase, B phase, C phase of impedance electric current phasor module vector, namely ${\stackrel{\·}{I}}_{\mathrm{ka}}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{ka}1}& {\stackrel{\·}{I}}_{\mathrm{ka}2}& ......& {\stackrel{\·}{I}}_{\mathrm{kar}}\end{array}\right]}^T,$ ${\stackrel{\·}{I}}_{\mathrm{kb}}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{kb}1}& {\stackrel{\·}{I}}_{\mathrm{kb}2}& ......& {\stackrel{\·}{I}}_{\mathrm{kbr}}\end{array}\right]}^T,{\stackrel{\·}{I}}_{\mathrm{kc}}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{kc}1}& {\stackrel{\·}{I}}_{\mathrm{kc}2}& ......& {\stackrel{\·}{I}}_{\mathrm{kcr}}\end{array}\right]}^T,$ Lower same.

In like manner, phase component method is adopted to obtain,

$\left[\begin{array}{c}\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}}\\ \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kb}}\\ \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kc}}\end{array}\right]=\left[\begin{array}{ccc}I& I& I\\ a^2& a& I\\ a& a^2& I\end{array}\right]\left[\begin{array}{c}\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}1}\\ \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}2}\\ \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}0}\end{array}\right]---\left(13\right)$

$\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{ka}}\\ {\stackrel{\·}{I}}_{\mathrm{kb}}\\ {\stackrel{\·}{I}}_{\mathrm{kc}}\end{array}\right]=\left[\begin{array}{ccc}I& I& I\\ a^2& a& I\\ a& a^2& I\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{ka}1}\\ {\stackrel{\·}{I}}_{\mathrm{ka}2}\\ {\stackrel{\·}{I}}_{\mathrm{ka}0}\end{array}\right]---\left(14\right)$

Wherein be respectively the phasor module of the positive-sequence component of A phase fault port Impedance both end voltage, negative sequence component and zero-sequence component, namely $\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}1}={\left[\begin{array}{cccc}\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}1-1}& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}2-1}& ......& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kat}-1}\end{array}\right]}^T,$ $\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}2}={\left[\begin{array}{cccc}\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}1-2}& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}2-2}& ......& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kat}-2}\end{array}\right]}^T,$ $\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}0}={\left[\begin{array}{cccc}\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}1-0}& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}2-0}& ......& \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{kat}-0}\end{array}\right]}^T;$ be respectively the phasor module of electric current positive-sequence component, negative sequence component and the zero-sequence component flowing through A phase of impedance, namely ${\stackrel{\·}{I}}_{\mathrm{ka}1}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{ka}1-1}& {\stackrel{\·}{I}}_{\mathrm{ka}2-1}& ......& {\stackrel{\·}{I}}_{\mathrm{kat}-1}\end{array}\right]}^T,{\stackrel{\·}{I}}_{\mathrm{ka}2}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{ka}1-2}& {\stackrel{\·}{I}}_{\mathrm{ka}2-2}& ......& {\stackrel{\·}{I}}_{\mathrm{kat}-2}\end{array}\right]}^T,$ ${\stackrel{\·}{I}}_{\mathrm{ka}0}={\left[\begin{array}{cccc}{\stackrel{\·}{I}}_{\mathrm{ka}1-0}& {\stackrel{\·}{I}}_{\mathrm{ka}2-0}& ......& {\stackrel{\·}{I}}_{\mathrm{kat}-0}\end{array}\right]}^T.$ Lower same.

Formula (13) and (14) are substituted into formula (12), obtains after arrangement

$\left[\begin{array}{c}\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}1}\\ \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}2}\\ \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}0}\end{array}\right]=\left[\begin{array}{ccc}{Z}_{\mathrm{ks}}& Z_{\mathrm{km}}& Z_{\mathrm{kn}}\\ Z_{\mathrm{kn}}& Z_{\mathrm{ks}}& Z_{\mathrm{km}}\\ Z_{\mathrm{km}}& Z_{\mathrm{kn}}& Z_{\mathrm{ks}}\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{ka}1}\\ {\stackrel{\·}{I}}_{\mathrm{ka}2}\\ {\stackrel{\·}{I}}_{\mathrm{ka}0}\end{array}\right]---\left(15\right)$

The block elements of formula (15) as shown in the formula

$\left\{\begin{array}{c}{Z}_{\mathrm{km}}=\frac{1}{3}(Z_{\mathrm{ka}}+Z_{\mathrm{kb}}+Z_{\mathrm{kc}})\\ Z_{\mathrm{kn}}=\frac{1}{3}(Z_{\mathrm{ka}}+a^2{Z}_{\mathrm{kb}}+a{Z}_{\mathrm{kc}})\\ Z_{\mathrm{ks}}=\frac{1}{3}(Z_{\mathrm{ka}}+a{Z}_{\mathrm{kb}}+a^2{Z}_{\mathrm{kc}})\end{array}\right.---\left(16\right)$

Wherein Z _km, Z _kn, Z _ksbe diagonal element, dimension equals fault tuple (i.e. r Fault sub-network).

At non-working port, fault branch electric current and the fault branch voltage relationship of grid side Negative-sequence Net and zero sequence net are

${\stackrel{\·}{I}}_{\mathrm{ka}2}=-Y_{p2}\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}2}---\left(17\right)$

${\stackrel{\·}{I}}_{\mathrm{ka}0}=-Y_{p0}\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}0}---\left(18\right)$

Wherein Y _p2and Y _p0be respectively negative phase-sequence admittance module and the zero sequence admittance module of normal-sub network.

Simultaneous formula (15)-(18), derivation obtains

$\left[\begin{array}{c}\mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}2}\\ \mathrm{\Δ}{\stackrel{\·}{U}}_{\mathrm{ka}0}\end{array}\right]={[I+\left[\begin{array}{cc}{Z}_{\mathrm{ks}}{Y}_{p2}& Z_{\mathrm{km}}{Y}_{p0}\\ Z_{\mathrm{kn}}{Y}_{p2}& Z_{\mathrm{ks}}{Y}_{p0}\end{array}\right]]}^{-1}\left[\begin{array}{c}{Z}_{\mathrm{kn}}\\ Z_{\mathrm{km}}\end{array}\right]{\stackrel{\·}{I}}_{\mathrm{ka}1}---\left(19\right)$

$Y_f={\{Z_{\mathrm{ks}}-\left[\begin{array}{cc}{Z}_{\mathrm{km}}{Y}_{p2}& Z_{\mathrm{kn}}{Y}_{p0}\end{array}\right]\×{[I+\left[\begin{array}{cc}{Z}_{\mathrm{ks}}{Y}_{p2}& Z_{\mathrm{km}}{Y}_{p0}\\ Z_{\mathrm{kn}}{Y}_{p2}& Z_{\mathrm{ks}}{Y}_{p0}\end{array}\right]]}^{-1}\left[\begin{array}{c}{Z}_{\mathrm{kn}}\\ Z_{\mathrm{km}}\end{array}\right]\}}^{-1}---\left(20\right)$

In the open conductors of electric system hybrid simulation calculates, by Y _fmerge and enter in electrical network positive sequence net admittance module Y.

Although the present invention discloses as above with preferred embodiment, and is not used to limit scope of the invention process.Any those of ordinary skill in the art, not departing from invention scope of the present invention, when doing a little improvement, namely every equal improvement done according to the present invention, should be scope of the present invention and contained.

Claims (6)

1. a method for the unified process of electric system hybrid simulation fault, is characterized in that: the method comprises the following steps:

In A, hybrid simulation, electromechanical transient side adopts electromechanical transient faster than real time simulation program TH-STBLT to set up electric power system alternating current network;

B, at non-working port place, electric power system alternating current network is divided into normal-sub network and multiport Fault sub-network;

C, in normal-sub network, grid parts adopts first-harmonic three sequence network to describe, the primary equipment dynamic element access positive sequence nets such as generator, dynamic load, power electronic equipment;

D, in multiport Fault sub-network, set up the fault model of power system network element according to failure message, derive the complex admittance unified modules Y of the multiple fault of multiport _f; Wherein, described failure message comprises short circuits information and open conductors information; Described fault model comprises short circuits model and open conductors model;

The complex admittance module of described short circuits is

$Y_f={\{Z_{\mathrm{ks}}-\left[\begin{array}{cc}{Z}_{\mathrm{km}}{Y}_{p2}& Z_{\mathrm{kn}}{Y}_{p0}\end{array}\right]\×{[I+\left[\begin{array}{cc}{Z}_{\mathrm{ks}}{Y}_{p2}& Z_{\mathrm{km}}{Y}_{p0}\\ Z_{\mathrm{kn}}{Y}_{p2}& (Z_{\mathrm{ks}}+{3Z}_{\mathrm{kg}})Y_{p0}\end{array}\right]]}^{-1}\left[\begin{array}{c}{Z}_{\mathrm{kn}}\\ Z_{\mathrm{km}}\end{array}\right]\}}^{-1}$

The complex admittance module of described open conductors is

$Y_f={\{Z_{\mathrm{ks}}-\left[\begin{array}{cc}{Z}_{\mathrm{km}}{Y}_{p2}& Z_{\mathrm{kn}}{Y}_{p0}\end{array}\right]\×{[I+\left[\begin{array}{cc}{Z}_{\mathrm{ks}}{Y}_{p2}& Z_{\mathrm{km}}{Y}_{p0}\\ Z_{\mathrm{kn}}{Y}_{p2}& {Z_{\mathrm{ks}}Y}_{p0}\end{array}\right]]}^{-1}\left[\begin{array}{c}{Z}_{\mathrm{kn}}\\ Z_{\mathrm{km}}\end{array}\right]\}}^{-1}$

E, by Y _fsuperpose in normal-sub network positive sequence net admittance module Y, calculate positive sequence voltage and the forward-order current of positive sequence net non-working port;

F, calculated the non-working port voltage of negative phase-sequence and zero-sequence network by positive sequence port-current in fault.

2. method according to claim 1, it is characterized in that, described short circuits model comprises ground short circuit fault, phase fault and three-phase imbalance load, short circuits model adopts a Y to connect three phase of impedance access non-working ports, and wherein Y is connect three phase of impedance neutral points and described by impedance earth.

3. method according to claim 1, it is characterized in that, described open conductors model comprises broken string, the fault branch that open conductors model increases can access faulty line by three phase of impedance and describe: under normal circumstances, fault branch resistance value is a very little value, and three phase of impedance arrange different values to be combined into different open conductors or disturbance.

4. method according to claim 1, is characterized in that, adopts phase component method, voltage and current is carried out three-phase and the conversion of three sequences.

5. method according to claim 1, is characterized in that, in the short circuits of electric system hybrid simulation calculates, by the complex admittance module Y of short circuits _fmerge and enter in the admittance module Y of electrical network positive sequence net.

6. method according to claim 1, is characterized in that, in the open conductors of electric system hybrid simulation calculates, by the complex admittance module Y of open conductors _fmerge and enter in electrical network positive sequence net admittance module Y.