CN103296649B - The lossy electric transmission line current traveling-wave differential protection method of the saturated impact of anti-current instrument transformer - Google Patents

Abstract

The invention discloses the lossy electric transmission line current traveling-wave differential protection method of the saturated impact of a kind of anti-current instrument transformer.First the method gathers voltage, the current traveling wave component of each sampling instant of transforming plant protecting installation place of transmission line two ends; take into account the impact of line loss; by the current traveling wave component of prediction algorithm computing electric power line one end, the current traveling wave component calculated is utilized to form current traveling wave differential protection criterion with this end line current traveling wave component obtained of sampling.The inventive method, without the need to being asked for the current traveling wave component of each sampling instant by interpolation arithmetic, decreases algorithm operation quantity, improves protection act speed.The inventive method can complete line fault traveling-wave protection function in the short-data windows in 1/4 cycle, and protection act speed is fast, and takes into account line loss impact, is therefore applicable to the traveling-wave protection of the whole transient fault process of various electric pressure transmission line.<!--1-->

Description

The lossy electric transmission line current traveling-wave differential protection method of the saturated impact of anti-current instrument transformer

Technical field

The present invention relates to Relay Protection Technology in Power System field, specifically relate to the lossy electric transmission line current traveling-wave differential protection method of the saturated impact of a kind of anti-current instrument transformer.

Background technology

Owing to not affecting by system operation mode and electric network composition and having natural phase-selecting function, current differential protection is the main protection of various electric pressure transmission line always.In 220kV and following electric pressure transmission line, because transmission line capacitance current along the line is very little, distributed capacitance is very little on the impact of current differential protection performance.The voltage of ultra-high/extra-high voltage transmission line of alternation current, current delivery have obvious wave process; capacitance current along the line is very large; the amplitude of the vector of two ends fundamental frequency steady-state current component is utilized to be faced with current differential protection starting current as the conventional current differential protection of actuating quantity large; and in order to prevent false protection; improving startup set point can cause again protection sensitivity not enough, governs the application of conventional current differential protection on ultra-high/extra-high voltage transmission line of alternation current.

When failure point of power transmission line is near current transformer installation place; because fault current is very large; current transformer is easily saturated; current differential protection misoperation when CT saturation can cause line protection external area error, current differential protection refused action when CT saturation can cause line protection troubles inside the sample space.

Traveling-wave differential protection considers the impact of distributed capacitance in protection algorism Mathematical Modeling, and the impact not by transmission line distributed capacitance in traveling-wave differential protection principle, has very high performance.Application number 200910034669.1 patent of invention " is applicable to the traveling-wave differential protection method of series capacitor compensated line " and solves the impact of distributed capacitance on differential protection performance; but for the situation that the time delay of row wave traveling is the non-integral multiple sampling interval; need the electric parameters obtained by interpolation arithmetic on each time point; very high to the requirement of protective device sample frequency; therefore very high to protective device hardware requirement; and each sampling time to carry out interpolation arithmetic; the required operand of protection algorism itself is large, cannot meet the requirement of protection quick-action." traveling-wave differential protection of UHV Transmission Line with Shunt Reactor " of Su Bin, Dong Xinzhou and Sun Yuan Zhang Fabiao and " traveling-wave differential protection based on wavelet transformation " of Su Bin, Dong Xinzhou and Sun Yuan Zhang Fabiao and application number 200410079501.X patent of invention " detection method of voltage zero cross near fault in travelling wave protection " are the electric parameters that the situation in non-integral multiple sampling interval also needs to be obtained by interpolation arithmetic on each time point for the time delay of row wave traveling, there is the problem that operand is large equally; Need to carry out wavelet transformation, desired data window is large, and protection detects that fault generation required time is long.

At present; the transmission line travelling wave differential protecting method that many scholars have proposed is that the situation in non-integral multiple sampling interval all needs to carry out interpolation arithmetic and asks electric parameters on its each time point to the time delay of row wave traveling; the operand of protection algorism own is large, requires high to protective device sampling hardware.Part transmission line travelling wave differential protecting method even needs to carry out wavelet transformation, and desired data window is large, extends protection and the time that fault occurs detected, cannot meet the requirement of relaying protection to quick-action.

Summary of the invention

The object of the invention is to the deficiency overcoming prior art existence, a kind of lossy electric transmission line current traveling-wave differential protection method without the need to carrying out interpolation arithmetic, the saturated impact of anti-current instrument transformer is provided.

The lossy electric transmission line current traveling-wave differential protection method of the saturated impact of anti-current instrument transformer, is characterized in that, comprises following sequential steps:

(1) protective device utilize the traveling wave electric amount being positioned at the m transforming plant protecting installation place at transmission line two ends and each sampling instant of n transforming plant protecting installation place calculate the m transforming plant protecting installation place of t sampling instant 0, α, β mould current traveling wave component i ' _m0(t), i ' _{m α}(t), i ' _{m β}(t):

$\begin{array}{c}{i}_{m0}^{\&prime;}\left(t\right)=\frac{u_{m0}\left(t\right)(1-\cos \left({\mathrm{\&omega;\&tau;}}_0\right))-u_{m0}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_0\right)}{(Z_{c0}+\frac{R_0}{4})(1+\cos \left({\mathrm{\&omega;\&tau;}}_0\right))}-\frac{i_{m0}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_0\right)}{1+\cos \left({\mathrm{\&omega;\&tau;}}_0\right)}\\ -i_{n0}\left(t\right)-\frac{i_{n0}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_0\right)}{1+\cos \left({\mathrm{\&omega;\&tau;}}_0\right)}+\frac{u_{n0}\left(t\right)(1-\cos \left({\mathrm{\&omega;\&tau;}}_0\right))-u_{n0}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_0\right)}{(Z_{c0}+\frac{R_0}{4})(1+\cos \left({\mathrm{\&omega;\&tau;}}_0\right))}\end{array}$

$\begin{array}{c}{i}_{\mathrm{m\&alpha;}}^{\&prime;}\left(t\right)=\frac{u_{\mathrm{m\&alpha;}}\left(t\right)(1-\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right))-u_{\mathrm{m\&alpha;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right)}{(Z_{\mathrm{c\&alpha;}}+\frac{R_{\mathrm{\&alpha;}}}{4})(1+\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right))}-\frac{i_{\mathrm{m\&alpha;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right)}{1+\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right)}\\ -i_{\mathrm{n\&alpha;}}\left(t\right)-\frac{i_{\mathrm{n\&alpha;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right)}{1+\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right)}+\frac{u_{\mathrm{n\&alpha;}}\left(t\right)(1-\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right))-u_{\mathrm{n\&alpha;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right)}{(Z_{\mathrm{c\&alpha;}}+\frac{R_{\mathrm{\&alpha;}}}{4})(1+\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right))}\end{array}$

$\begin{array}{c}{i}_{\mathrm{m\&beta;}}^{\&prime;}\left(t\right)=\frac{u_{\mathrm{m\&beta;}}\left(t\right)(1-\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right))-u_{\mathrm{m\&beta;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right)}{(Z_{\mathrm{c\&beta;}}+\frac{R_{\mathrm{\&beta;}}}{4})(1+\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right))}-\frac{i_{\mathrm{m\&beta;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right)}{1+\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right)}\\ -i_{\mathrm{n\&beta;}}\left(t\right)-\frac{i_{\mathrm{n\&beta;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right)}{1+\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right)}+\frac{u_{\mathrm{n\&beta;}}\left(t\right)(1-\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right))-u_{\mathrm{n\&beta;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right)}{(Z_{\mathrm{c\&beta;}}+\frac{R_{\mathrm{\&beta;}}}{4})(1+\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right))}\end{array}$

Wherein, t is the sampling time; , , ; L is the transmission line connecting m transformer station and n transformer station

Length; T is the cycle time of fundamental component; Z _c0, Z _{c α}, Z _{c β}be respectively the characteristic impedance of transmission line 0, α, β line wave component; v _c0, v _{c α}, v _{c β}be respectively the propagation velocity of transmission line 0, α, β line wave component; ω is electric power system angular frequency; R ₀, R _α, R _βbe respectively the resistance of transmission line 0, α, β line wave component; u _m0(t), u _{m α}(t), u _{m β}t () is respectively 0 of the t sampling instant of m transforming plant protecting installation place, the voltage traveling wave component of α, β mould; i _m0(t), i _{m α}(t), i _{m β}t () is respectively 0 of the t sampling instant of m transforming plant protecting installation place, the current traveling wave component of α, β mould; u _n0(t), u _{n α}(t), u _{n β}t () is respectively 0 of the t sampling instant of n transforming plant protecting installation place, the voltage traveling wave component of α, β mould; i _n0(t), i _{n α}(t), i _{n β}t () is respectively 0 of the t sampling instant of n transforming plant protecting installation place, the current traveling wave component of α, β mould; , , be respectively m transforming plant protecting installation place 0 of sampling instant, the voltage traveling wave component of α, β mould; , , be respectively n transforming plant protecting installation place 0 of sampling instant, the voltage traveling wave component of α, β mould; , , be respectively m transforming plant protecting installation place 0 of sampling instant, the current traveling wave component of α, β mould; , , be respectively n transforming plant protecting installation place 0 of sampling instant, the current traveling wave component of α, β mould;

(2) protective device is by i ' _m0(t), i ' _{m α}(t), i ' _{m β}t () carries out the three-phase current traveling-wave component i ' that phase mould inverse transformation obtains the m transforming plant protecting installation place of t sampling instant _mA(t), i ' _mB(t), i ' _mC(t); To i ' _mA(t), i ' _mB(t), i ' _mCt () adopts Fourier algorithm to calculate the three-phase Fundamental-frequency Current component of the t sampling instant of m transforming plant protecting installation place , , ; To the current traveling wave component i that the three-phase of the m transforming plant protecting installation place of t sampling instant is surveyed _mA(t), i _mB(t), i _mCt () adopts Fourier algorithm to calculate the Fundamental-frequency Current component of the three-phase actual measurement of the t sampling instant of m transforming plant protecting installation place , , ;

(3) protective device judges $\cos {\mathrm{\&theta;}}_A\sqrt{{\left[\mathrm{Re}({\stackrel{.}{I}}_{\mathrm{mA}}^{\&prime;}+{\stackrel{.}{I}}_{\mathrm{mA}})\right]}^2+{\left[\mathrm{Im}({\stackrel{.}{I}}_{\mathrm{mA}}^{\&prime;}+{\stackrel{.}{I}}_{\mathrm{mA}})\right]}^2}>k\sqrt{{\left[\mathrm{Re}({\stackrel{.}{I}}_{\mathrm{mA}}^{\&prime;}-{\stackrel{.}{I}}_{\mathrm{mA}})\right]}^2+{\left[\mathrm{Im}({\stackrel{.}{I}}_{\mathrm{mA}}^{\&prime;}-{\stackrel{.}{I}}_{\mathrm{mA}})\right]}^2}$ Whether set up, if set up, then protective device sends trip signal, the circuit breaker at tripping A phase transmission line two ends to the circuit breaker at A phase transmission line two ends; Wherein, θ _afor leading angle; K is tuning coefficient; for imaginary part; for real part; for real part; for imaginary part;

(4) protective device judges $\cos {\mathrm{\&theta;}}_B\sqrt{{\left[\mathrm{Re}({\stackrel{.}{I}}_{\mathrm{mB}}^{\&prime;}+{\stackrel{.}{I}}_{\mathrm{mB}})\right]}^2+{\left[\mathrm{Im}({\stackrel{.}{I}}_{\mathrm{mB}}^{\&prime;}+{\stackrel{.}{I}}_{\mathrm{mB}})\right]}^2}>k\sqrt{{\left[\mathrm{Re}({\stackrel{.}{I}}_{\mathrm{mB}}^{\&prime;}-{\stackrel{.}{I}}_{\mathrm{mB}})\right]}^2+{\left[\mathrm{Im}({\stackrel{.}{I}}_{\mathrm{mB}}^{\&prime;}-{\stackrel{.}{I}}_{\mathrm{mB}})\right]}^2}$ Whether set up, if set up, then protective device sends trip signal, the circuit breaker at tripping B phase transmission line two ends to the circuit breaker at B phase transmission line two ends; Wherein, θ _bfor leading angle; K is tuning coefficient; for imaginary part; for real part; for real part; for imaginary part;

(5) protective device judges $\cos {\mathrm{\&theta;}}_C\sqrt{{\left[\mathrm{Re}({\stackrel{.}{I}}_{\mathrm{mC}}^{\&prime;}+{\stackrel{.}{I}}_{\mathrm{mC}})\right]}^2+{\left[\mathrm{Im}({\stackrel{.}{I}}_{\mathrm{mC}}^{\&prime;}+{\stackrel{.}{I}}_{\mathrm{mC}})\right]}^2}>k\sqrt{{\left[\mathrm{Re}({\stackrel{.}{I}}_{\mathrm{mC}}^{\&prime;}-{\stackrel{.}{I}}_{\mathrm{mC}})\right]}^2+{\left[\mathrm{Im}({\stackrel{.}{I}}_{\mathrm{mC}}^{\&prime;}-{\stackrel{.}{I}}_{\mathrm{mC}})\right]}^2}$ Whether set up, if set up, then protective device sends trip signal, the circuit breaker at tripping C phase transmission line two ends to the circuit breaker at C phase transmission line two ends; Wherein, θ _cfor leading angle; K is tuning coefficient; for imaginary part; for real part; for real part; for imaginary part.

The present invention compared with prior art, has following positive achievement:

The inventive method, without the need to being asked for the current traveling wave component of each sampling instant by interpolation arithmetic, decreases algorithm operation quantity, improves protection act speed.The inventive method can complete transmission line malfunction traveling-wave protection function in the short-data windows in 1/4 cycle, and protection act speed is fast, and takes into account line loss impact, is therefore applicable to the traveling-wave protection of the whole transient fault process of various electric pressure transmission line.The inventive method eliminates the impact of CT saturation on travelling wave current differential protection performance; no matter whether current transformer is saturated; line protection external area error the inventive method is reliably failure to actuate, the inventive method correct action message during line protection troubles inside the sample space.

Accompanying drawing explanation

Fig. 1 is application multi-line power transmission system schematic of the present invention.

Embodiment

According to Figure of description, technical scheme of the present invention is expressed in further detail below.

Fig. 1 is application multi-line power transmission system schematic of the present invention.In the present embodiment, first protective device gathers transmission line at the three-phase voltage traveling-wave component of each sampling instant of m transforming plant protecting installation place and three-phase current traveling-wave component; Gather transmission line at the three-phase voltage traveling-wave component of each sampling instant of n transforming plant protecting installation place and three-phase current traveling-wave component.

Protective device adopt phase-model transformation the three-phase voltage traveling-wave component of each sampling instant of m transforming plant protecting installation place, three-phase current traveling-wave component are converted to m transforming plant protecting installation place 0, α, β mode voltage traveling-wave component and 0, α, β mould current traveling wave component.

Protective device adopt phase-model transformation the three-phase voltage traveling-wave component of each sampling instant of n transforming plant protecting installation place, three-phase current traveling-wave component are converted to n transforming plant protecting installation place 0, α, β mode voltage traveling-wave component and 0, α, β mould current traveling wave component.

Protective device utilize the traveling wave electric amount being positioned at the m transforming plant protecting installation place at transmission line two ends and each sampling instant of n transforming plant protecting installation place calculate the m transforming plant protecting installation place of t sampling instant 0, α, β mould current traveling wave component i ' _m0(t), i ' _{m α}(t), i ' _{m β}(t):

$\begin{array}{c}{i}_{m0}^{\&prime;}\left(t\right)=\frac{u_{m0}\left(t\right)(1-\cos \left({\mathrm{\&omega;\&tau;}}_0\right))-u_{m0}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_0\right)}{(Z_{c0}+\frac{R_0}{4})(1+\cos \left({\mathrm{\&omega;\&tau;}}_0\right))}-\frac{i_{m0}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_0\right)}{1+\cos \left({\mathrm{\&omega;\&tau;}}_0\right)}\\ -i_{n0}\left(t\right)-\frac{i_{n0}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_0\right)}{1+\cos \left({\mathrm{\&omega;\&tau;}}_0\right)}+\frac{u_{n0}\left(t\right)(1-\cos \left({\mathrm{\&omega;\&tau;}}_0\right))-u_{n0}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_0\right)}{(Z_{c0}+\frac{R_0}{4})(1+\cos \left({\mathrm{\&omega;\&tau;}}_0\right))}\end{array}$

$\begin{array}{c}{i}_{\mathrm{m\&alpha;}}^{\&prime;}\left(t\right)=\frac{u_{\mathrm{m\&alpha;}}\left(t\right)(1-\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right))-u_{\mathrm{m\&alpha;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right)}{(Z_{\mathrm{c\&alpha;}}+\frac{R_{\mathrm{\&alpha;}}}{4})(1+\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right))}-\frac{i_{\mathrm{m\&alpha;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right)}{1+\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right)}\\ -i_{\mathrm{n\&alpha;}}\left(t\right)-\frac{i_{\mathrm{n\&alpha;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right)}{1+\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right)}+\frac{u_{\mathrm{n\&alpha;}}\left(t\right)(1-\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right))-u_{\mathrm{n\&alpha;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right)}{(Z_{\mathrm{c\&alpha;}}+\frac{R_{\mathrm{\&alpha;}}}{4})(1+\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right))}\end{array}$

$\begin{array}{c}{i}_{\mathrm{m\&beta;}}^{\&prime;}\left(t\right)=\frac{u_{\mathrm{m\&beta;}}\left(t\right)(1-\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right))-u_{\mathrm{m\&beta;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right)}{(Z_{\mathrm{c\&beta;}}+\frac{R_{\mathrm{\&beta;}}}{4})(1+\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right))}-\frac{i_{\mathrm{m\&beta;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right)}{1+\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right)}\\ -i_{\mathrm{n\&beta;}}\left(t\right)-\frac{i_{\mathrm{n\&beta;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right)}{1+\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right)}+\frac{u_{\mathrm{n\&beta;}}\left(t\right)(1-\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right))-u_{\mathrm{n\&beta;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right)}{(Z_{\mathrm{c\&beta;}}+\frac{R_{\mathrm{\&beta;}}}{4})(1+\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right))}\end{array}$

Wherein, t is the sampling time; , , ; L is the transmission line connecting m transformer station and n transformer station

Length; T is the cycle time of fundamental component; Z _c0, Z _{c α}, Z _{c β}be respectively the characteristic impedance of transmission line 0, α, β line wave component; v _c0, v _{c α}, v _{c β}be respectively the propagation velocity of transmission line 0, α, β line wave component; ω is electric power system angular frequency; R ₀, R _α, R _βbe respectively the resistance of transmission line 0, α, β line wave component; u _m0(t), u _{m α}(t), u _{m β}t () is respectively 0 of the t sampling instant of m transforming plant protecting installation place, the voltage traveling wave component of α, β mould; i _m0(t), i _{m α}(t), i _{m β}t () is respectively 0 of the t sampling instant of m transforming plant protecting installation place, the current traveling wave component of α, β mould; u _n0(t), u _{n α}(t), u _{n β}t () is respectively 0 of the t sampling instant of n transforming plant protecting installation place, the voltage traveling wave component of α, β mould; i _n0(t), i _{n α}(t), i _{n β}t () is respectively 0 of the t sampling instant of n transforming plant protecting installation place, the current traveling wave component of α, β mould; , , be respectively m transforming plant protecting installation place 0 of sampling instant, the voltage traveling wave component of α, β mould; , , be respectively n transforming plant protecting installation place 0 of sampling instant, the voltage traveling wave component of α, β mould; , , be respectively m transforming plant protecting installation place 0 of sampling instant, the current traveling wave component of α, β mould; , , be respectively n transforming plant protecting installation place 0 of sampling instant, the current traveling wave component of α, β mould.

Protective device is by i ' _m0(t), i ' _{m α}(t), i ' _{m β}t () carries out the three-phase current traveling-wave component i ' that phase mould inverse transformation obtains the m transforming plant protecting installation place of t sampling instant _mA(t), i ' _mB(t), i ' _mC(t).

Protective device is to i ' _mA(t), i ' _mB(t), i ' _mCt () adopts Fourier algorithm to calculate the three-phase Fundamental-frequency Current component of the t sampling instant of m transforming plant protecting installation place , , .

The current traveling wave component i that protective device is surveyed the three-phase of the m transforming plant protecting installation place of t sampling instant _mA(t), i _mB(t), i _mCt () adopts Fourier algorithm to calculate the Fundamental-frequency Current component of the three-phase actual measurement of the t sampling instant of m transforming plant protecting installation place , , .

Protective device judges $\cos {\mathrm{\&theta;}}_A\sqrt{{\left[\mathrm{Re}({\stackrel{.}{I}}_{\mathrm{mA}}^{\&prime;}+{\stackrel{.}{I}}_{\mathrm{mA}})\right]}^2+{\left[\mathrm{Im}({\stackrel{.}{I}}_{\mathrm{mA}}^{\&prime;}+{\stackrel{.}{I}}_{\mathrm{mA}})\right]}^2}>k\sqrt{{\left[\mathrm{Re}({\stackrel{.}{I}}_{\mathrm{mA}}^{\&prime;}-{\stackrel{.}{I}}_{\mathrm{mA}})\right]}^2+{\left[\mathrm{Im}({\stackrel{.}{I}}_{\mathrm{mA}}^{\&prime;}-{\stackrel{.}{I}}_{\mathrm{mA}})\right]}^2}$ Whether set up, if set up, then protective device sends trip signal, the circuit breaker at tripping A phase transmission line two ends to the circuit breaker at A phase transmission line two ends; Wherein, θ _afor leading angle; K is tuning coefficient; for imaginary part; for real part; for real part; for imaginary part.

Protective device judges $\cos {\mathrm{\&theta;}}_B\sqrt{{\left[\mathrm{Re}({\stackrel{.}{I}}_{\mathrm{mB}}^{\&prime;}+{\stackrel{.}{I}}_{\mathrm{mB}})\right]}^2+{\left[\mathrm{Im}({\stackrel{.}{I}}_{\mathrm{mB}}^{\&prime;}+{\stackrel{.}{I}}_{\mathrm{mB}})\right]}^2}>k\sqrt{{\left[\mathrm{Re}({\stackrel{.}{I}}_{\mathrm{mB}}^{\&prime;}-{\stackrel{.}{I}}_{\mathrm{mB}})\right]}^2+{\left[\mathrm{Im}({\stackrel{.}{I}}_{\mathrm{mB}}^{\&prime;}-{\stackrel{.}{I}}_{\mathrm{mB}})\right]}^2}$ Whether set up, if set up, then protective device sends trip signal, the circuit breaker at tripping B phase transmission line two ends to the circuit breaker at B phase transmission line two ends; Wherein, θ _bfor leading angle; K is tuning coefficient; for imaginary part; for real part; for real part; for imaginary part.

Protective device judges $\cos {\mathrm{\&theta;}}_C\sqrt{{\left[\mathrm{Re}({\stackrel{.}{I}}_{\mathrm{mC}}^{\&prime;}+{\stackrel{.}{I}}_{\mathrm{mC}})\right]}^2+{\left[\mathrm{Im}({\stackrel{.}{I}}_{\mathrm{mC}}^{\&prime;}+{\stackrel{.}{I}}_{\mathrm{mC}})\right]}^2}>k\sqrt{{\left[\mathrm{Re}({\stackrel{.}{I}}_{\mathrm{mC}}^{\&prime;}-{\stackrel{.}{I}}_{\mathrm{mC}})\right]}^2+{\left[\mathrm{Im}({\stackrel{.}{I}}_{\mathrm{mC}}^{\&prime;}-{\stackrel{.}{I}}_{\mathrm{mC}})\right]}^2}$ Whether set up, if set up, then protective device sends trip signal, the circuit breaker at tripping C phase transmission line two ends to the circuit breaker at C phase transmission line two ends; Wherein, θ _cfor leading angle; for tuning coefficient; for imaginary part; for real part; for real part; for imaginary part.

The inventive method, without the need to being asked for the current traveling wave component of each sampling instant by interpolation arithmetic, decreases algorithm operation quantity, improves protection act speed.The inventive method can complete line fault traveling-wave protection function in the short-data windows in 1/4 cycle, and protection act speed is fast, and takes into account line loss impact, is therefore applicable to the traveling-wave protection of the whole transient fault process of various electric pressure transmission line.The inventive method eliminates the impact of CT saturation on travelling wave current differential protection performance; no matter whether current transformer is saturated; line protection external area error the inventive method is reliably failure to actuate, the inventive method correct action message during line protection troubles inside the sample space.

The foregoing is only preferred embodiment of the present invention; but protection scope of the present invention is not limited thereto; anyly be familiar with those skilled in the art in the technical scope that the present invention discloses, the change that can expect easily or replacement, all should be encompassed within protection scope of the present invention.

Claims (1)

1. the lossy electric transmission line current traveling-wave differential protection method of the saturated impact of anti-current instrument transformer, is characterized in that, comprise the following steps:

(1) protective device utilize the traveling wave electric amount being positioned at the m transforming plant protecting installation place at transmission line two ends and each sampling instant of n transforming plant protecting installation place calculate the m transforming plant protecting installation place of t sampling instant 0, α, β mould current traveling wave component i ' _m0(t), i ' _{m α}(t), i ' _{m β}(t):

$\begin{array}{c}{i}_{m0}^{\&prime;}\left(t\right)=\frac{u_{m0}\left(t\right)\left(1-\cos \left({\mathrm{\&omega;\&tau;}}_0\right)\right)-u_{m0}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_0\right)}{\left(Z_{c0}+\frac{R_0}{4}\right)(1+\cos ({\mathrm{\&omega;\&tau;}}_0\left)\right)}-\frac{i_{m0}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_0\right)}{1+\cos \left({\mathrm{\&omega;\&tau;}}_0\right)}\\ -i_{n0}\left(t\right)-\frac{i_{n0}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_0\right)}{1+\cos \left({\mathrm{\&omega;\&tau;}}_0\right)}+\frac{u_{n0}\left(t\right)\left(1-\cos \left({\mathrm{\&omega;\&tau;}}_0\right)\right)-u_{n0}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_0\right)}{\left(Z_{c0}+\frac{R_0}{4}\right)(1+\cos ({\mathrm{\&omega;\&tau;}}_0\left)\right)}\end{array}$

$\begin{array}{c}{i}_{m\mathrm{\&alpha;}}^{\&prime;}\left(t\right)=\frac{u_{m\mathrm{\&alpha;}}\left(t\right)\left(1-\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right)\right)-u_{m\mathrm{\&alpha;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right)}{\left(Z_{c\mathrm{\&alpha;}}+\frac{R_{\mathrm{\&alpha;}}}{4}\right)(1+\cos ({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\left)\right)}-\frac{i_{m\mathrm{\&alpha;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right)}{1+\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right)}\\ -i_{n\mathrm{\&alpha;}}\left(t\right)-\frac{i_{n\mathrm{\&alpha;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right)}{1+\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right)}+\frac{u_{n\mathrm{\&alpha;}}\left(t\right)\left(1-\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right)\right)-u_{n\mathrm{\&alpha;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\right)}{\left(Z_{c\mathrm{\&alpha;}}+\frac{R_{\mathrm{\&alpha;}}}{4}\right)(1+\cos ({\mathrm{\&omega;\&tau;}}_{\mathrm{\&alpha;}}\left)\right)}\end{array}$

$\begin{array}{c}{i}_{m\mathrm{\&beta;}}^{\&prime;}\left(t\right)=\frac{u_{m\mathrm{\&beta;}}\left(t\right)\left(1-\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right)\right)-u_{m\mathrm{\&beta;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right)}{\left(Z_{c\mathrm{\&beta;}}+\frac{R_{\mathrm{\&beta;}}}{4}\right)(1+\cos ({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\left)\right)}-\frac{i_{m\mathrm{\&beta;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right)}{1+\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right)}\\ -i_{n\mathrm{\&beta;}}\left(t\right)-\frac{i_{n\mathrm{\&beta;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right)}{1+\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right)}+\frac{u_{n\mathrm{\&beta;}}\left(t\right)\left(1-\cos \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right)\right)-u_{n\mathrm{\&beta;}}(t-\frac{T}{4})\sin \left({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\right)}{\left(Z_{c\mathrm{\&beta;}}+\frac{R_{\mathrm{\&beta;}}}{4}\right)(1+\cos ({\mathrm{\&omega;\&tau;}}_{\mathrm{\&beta;}}\left)\right)}\end{array}$

Wherein, t is sampling instant; l is the transmission line length connecting m transformer station and n transformer station; T is the cycle time of fundamental component; Z _c0, Z _{c α}, Z _{c β}be respectively the characteristic impedance of transmission line 0, α, β line wave component; ν _c0, ν _{c α}, ν _{c β}be respectively the propagation velocity of transmission line 0, α, β line wave component; ω is electric power system angular frequency; R ₀, R _α, R _βbe respectively the resistance of transmission line 0, α, β line wave component; u _m0(t), u _{m α}(t), u _{m β}t () is respectively 0 of the t sampling instant of m transforming plant protecting installation place, the voltage traveling wave component of α, β mould; i _m0(t), i _{m α}(t), i _{m β}t () is respectively 0 of the t sampling instant of m transforming plant protecting installation place, the current traveling wave component of α, β mould; u _n0(t), u _{n α}(t), u _{n β}t () is respectively 0 of the t sampling instant of n transforming plant protecting installation place, the voltage traveling wave component of α, β mould; i _n0(t), i _{n α}(t), i _{n β}t () is respectively 0 of the t sampling instant of n transforming plant protecting installation place, the current traveling wave component of α, β mould; be respectively m transforming plant protecting installation place 0 of sampling instant, the voltage traveling wave component of α, β mould; be respectively n transforming plant protecting installation place 0 of sampling instant, the voltage traveling wave component of α, β mould; be respectively m transforming plant protecting installation place 0 of sampling instant, the current traveling wave component of α, β mould; be respectively n transforming plant protecting installation place 0 of sampling instant, the current traveling wave component of α, β mould;

(2) protective device is by i ' _m0(t), i ' _{m α}(t), i ' _{m β}t () carries out the three-phase current traveling-wave component i ' that phase mould inverse transformation obtains the m transforming plant protecting installation place of t sampling instant _mA(t), i ' _mB(t), i ' _mC(t); To i ' _mA(t), i ' _mB(t), i ' _mCt () adopts Fourier algorithm to calculate the three-phase Fundamental-frequency Current component of the t sampling instant of m transforming plant protecting installation place to the current traveling wave component i that the three-phase of the m transforming plant protecting installation place of t sampling instant is surveyed _mA(t), i _mB(t), i _mCt () adopts Fourier algorithm to calculate the Fundamental-frequency Current component of the three-phase actual measurement of the t sampling instant of m transforming plant protecting installation place

(3) protective device judges ${\mathrm{cos\&theta;}}_A\sqrt{{\&lsqb;\mathrm{Re}({\stackrel{\&CenterDot;}{I}}_{mA}^{\&prime;}+{\stackrel{\&CenterDot;}{I}}_{mA})\&rsqb;}^2+{\&lsqb;\mathrm{Im}({\stackrel{\&CenterDot;}{I}}_{mA}^{\&prime;}+{\stackrel{\&CenterDot;}{I}}_{mA})\&rsqb;}^2}>k\sqrt{{\&lsqb;\mathrm{Re}({\stackrel{\&CenterDot;}{I}}_{mA}^{\&prime;}-{\stackrel{\&CenterDot;}{I}}_{mA})\&rsqb;}^2+{\&lsqb;\mathrm{Im}({\stackrel{\&CenterDot;}{I}}_{mA}^{\&prime;}-{\stackrel{\&CenterDot;}{I}}_{mA})\&rsqb;}^2}$ Whether set up, if set up, then protective device sends trip signal, the circuit breaker at tripping A phase transmission line two ends to the circuit breaker at A phase transmission line two ends; Wherein, θ _afor leading angle; K is tuning coefficient; for imaginary part; for real part; for real part; for imaginary part;

(4) protective device judges ${\mathrm{cos\&theta;}}_B\sqrt{{\&lsqb;\mathrm{Re}({\stackrel{\&CenterDot;}{I}}_{mB}^{\&prime;}+{\stackrel{\&CenterDot;}{I}}_{mB})\&rsqb;}^2+{\&lsqb;\mathrm{Im}({\stackrel{\&CenterDot;}{I}}_{mB}^{\&prime;}+{\stackrel{\&CenterDot;}{I}}_{mB})\&rsqb;}^2}>k\sqrt{{\&lsqb;\mathrm{Re}({\stackrel{\&CenterDot;}{I}}_{mB}^{\&prime;}-{\stackrel{\&CenterDot;}{I}}_{mB})\&rsqb;}^2+{\&lsqb;\mathrm{Im}({\stackrel{\&CenterDot;}{I}}_{mB}^{\&prime;}-{\stackrel{\&CenterDot;}{I}}_{mB})\&rsqb;}^2}$ Whether set up, if set up, then protective device sends trip signal, the circuit breaker at tripping B phase transmission line two ends to the circuit breaker at B phase transmission line two ends; Wherein, θ _bfor leading angle; K is tuning coefficient; for imaginary part; for real part; for real part; for imaginary part;

(5) protective device judges ${\mathrm{cos\&theta;}}_C\sqrt{{\&lsqb;\mathrm{Re}({\stackrel{\&CenterDot;}{I}}_{mC}^{\&prime;}+{\stackrel{\&CenterDot;}{I}}_{mC})\&rsqb;}^2+{\&lsqb;\mathrm{Im}({\stackrel{\&CenterDot;}{I}}_{mC}^{\&prime;}+{\stackrel{\&CenterDot;}{I}}_{mC})\&rsqb;}^2}>k\sqrt{{\&lsqb;\mathrm{Re}({\stackrel{\&CenterDot;}{I}}_{mC}^{\&prime;}-{\stackrel{\&CenterDot;}{I}}_{mC})\&rsqb;}^2+{\&lsqb;\mathrm{Im}({\stackrel{\&CenterDot;}{I}}_{mC}^{\&prime;}-{\stackrel{\&CenterDot;}{I}}_{mC})\&rsqb;}^2}$ Whether set up, if set up, then protective device sends trip signal, the circuit breaker at tripping C phase transmission line two ends to the circuit breaker at C phase transmission line two ends; Wherein, θ _cfor leading angle; K is tuning coefficient; for imaginary part; for real part; for real part; for imaginary part.