EP0619922A1

EP0619922A1 - Reconstruction of saturated current transformer signals

Info

Publication number: EP0619922A1
Application number: EP93901453A
Authority: EP
Inventors: Magnus BJÖRKLUND; Bengt Carlsson; Murari Mohan Saha
Original assignee: Asea Brown Boveri AB
Current assignee: ABB AB
Priority date: 1992-01-03
Filing date: 1992-12-14
Publication date: 1994-10-19
Also published as: SE469676B; WO1993013581A1; CA2126867A1; SE9200011L; SE9200011D0

Abstract

L'invention se rapporte à un procédé numérique qui, par l'intermédiaire de mesures en temps discontinu du courant secondaire (is(k)) à partir d'un transformateur de courant, reconstitue le courant primaire (ipr(k)) et le courant secondaire (isr(k)) d'un transformateur de courant saturé qui a lieu dans un bloc B1. L'invention se rapporte également à un procédé permettant de détecter la saturation et l'arrêt de la saturation qui a lieu dans les blocs B2, B3, B4, B5, B6 et B7. Les signaux de sortie du procédé sont obtenus par un bloc B8 et sur une décision prise quant à la saturation du transformateur de courant. Les signaux de sortie se composent du courant primaire reconstitué (ipn(k)) et du courant secondaire reconstitué (isn(k)) sont obtenus et sur une décision prise quant à la non-saturation du transformateur de courant. Les signaux de sortie se composent du courant primaire reconstitué (ipn(k)) et du courant secondaire mesuré (ip(k)). L'invention se rapporte également à un dispositif mettant en oeuvre ces procédés.The invention relates to a digital method which, by means of discontinuous time measurements of the secondary current (is(k)) from a current transformer, reconstructs the primary current (ipr(k)) and the secondary current (isr(k)) of a saturated current transformer which takes place in a block B1. The invention also relates to a method making it possible to detect the saturation and the cessation of the saturation which takes place in the blocks B2, B3, B4, B5, B6 and B7. The output signals of the method are obtained by a block B8 and on a decision taken as to the saturation of the current transformer. The output signals consist of the reconstituted primary current (ipn(k)) and the reconstituted secondary current (isn(k)) are obtained and on a decision taken as to whether the current transformer is not saturated. The output signals consist of the reconstituted primary current (ipn(k)) and the measured secondary current (ip(k)). The invention also relates to a device implementing these methods.

Description

Reconstruction of saturated current transformer signals

TECHNICAL FIELD Within largely all distribution of electric power, monitoring systems for detecting short circuits and other abnormal states are needed. Instrument transformers for measuring current and voltage therefore constitute an important and integrated part of control, monitoring and protective devices in most power distribution systems. It is, of course, very important that these measuring devices correctly reproduce the quantities they are designed to measure, both with regard to static and dynamic values. Since both control and monitoring devices and particularly protective devices are nowadays based on instantaneous value measurement, very high demands are placed on good dynamic conformity during the measurements.

Due to practical and economic reasons it is not possible to dimension current transformers so as to avoid saturation at all fault current levels. Therefore, during a fault condition a current of such a high value may arise that the current transformer becomes magnetically saturated, whereby the waveform of the delivered current becomes distorted. This leads to considerable problems, especially for current measuring numeral protective relays, since these will then receive false information on which to base the decisions to take action. A typical example is a differential protection device which, when a current transformer saturation occurs, may lead to an external fault being incorrectly detected as an intersystem fault. It is therefore of great importance to be able to detect whether a current transformer has become subjected to such a current that it has become magnetically saturated. The present invention describes a method for reconstructing, during saturation, with the aid of the distorted secondary current, the primary current which causes the saturation. A reconstruction is also performed of the secondary current which would have been obtained unless the current transformer had become saturated. The invention also comprises a method for detecting saturation and detecting when the saturation ceases. The invention also comprises a device for carrying out the method.

BRIEF DESCRIPTION OF THE DRAWINGS

Figure 1 shows an equivalent diagram for a current transformer seen from the secondary side.

Figure 2 shows a flow diagram according to the invention.

BACKGROUND ART, THE PROBLEMS One way of reducing the effect of a current transformer saturation is to use some method for detecting saturation. Several such methods exist, which is clear, inter alia, from SE 90100917-5. When a saturation has been determined, the conditions for breaking the current are made more conservative in existing protective relays. This assumption, however, has some apparent shortcomings; for example, it means that a short-circuit current may be allowed to act for an unacceptably long period of time, which in turn may lead to negative effects on the network. A further problem in this connection is to detect when the saturation ceases.

A general description of current transformers, both in steady and transient states, is given by A Wright in "Current Transformers, Their Transient and Steady State Performance", Chapman and Hall Ltd., 1968. A model which is suitable for further study of the behaviour of current transformers at different magnetizing states is described, inter alia, by T Conrad, J Schlabbach and R Speh in "Verfahren zur Korrektur der verzerrten Sekundarstrόme von Stromwandlern", etzArchiv Bd., Vol. 6, 1984, in which it is referred to an equivalent diagram according to Figure 1 where i_p = stepped down primary current i_s = secondary current L_μ = magnetizing inductance

R_fe = loss resistance in the iron core R₂ = resistance of burden, winding, and conductors L₂ = inductance of burden and leakage inductance N₂/N₁ = current transformer ratio The differential equations for the circuit according to Figure 1 will then be

(1) h

(2)

It should be pointed out here that the values of both L_μ and R_fe vary in dependence on the primary current via the magnetizing current. In the following description, however, R_fe will be assumed to be constant. When the magnetizing current exceeds a certain value, L_μ will decrease drastically. The greater part of the stepped down current will then be conducted via the magnetizing inductance, which results in the severely distorted secondary current. In the following description, the term "primary current" will be used only for "the stepped down primary current".

It should be pointed out here that the magnetizing current i_μ is needed in explicit form to be able to calculate L_μ. The differential equations (1) and (2) can be solved with numerical methods, which, inter alia, is clear from an article by T Conrad, J. Schlabbach and R Speh: "Verfahren zur Korrektur der verzerrten Sekundarstrόme von Stromwandlern" published in etzArchiv Bd., Vol 6, pp 77-79, 1984. In these calculations a special algorithm has been used for L_μ which provides a saturated value which is only 60 times lower than the unsaturated value. The result of the calculations and the simulations which have been made show a fairly good reconstruction of the primary current. The value of the method described must, however, be considered to be relatively limited. This is due, among other thing, to the saturated value of L_μ in reality being rather 1000 times lower than the unsaturated value. Considering, in addition, the transfer functions corresponding to the differential equations, it will be seen that a high gain is required to reconstruct the primary current. This means that the method used is probably sensitive to measurement disturbances and since no disturbances are used during the simulations, the usefulness of this method in practice must be considered doubtful.

Differential equations of the types according to (1) and (2) can be described in state-space form according to known technique, see, for example, T Kailath, Linear Systems, Prentice-Hall. Inc., Englewood Cliffs, N. J. 1980, pp. 50- 60. A general state-space model for the current transformer is given by dx/dt = Ac(L_μ)x + B_c(L_μ)i_p (3) i_s = C_s x (4) i_μ = C_μ x ( 5 ) where it has been indicated that the state matrix A_c and the vector B_c depend on L_μ, which in turn is dependent on the magnetizing current i_μ. The selections of state-space representations which are utilized in the invention will be described under the summary of the invention.

Equations (3), (4) and (5) describe the model in continuous time. In "Computer Controlled Systems", p. 37, Prentice Hall, Inc., 1984, by K J Åstrόm and B Wittenmark, there is described a method which discretizes the equations, whereby the following sampled state model is obtained. x(k + 1) = Ad(L_μ) x(k) + B_d(L_μ) i_p(k) (6) i_s(k) = C_s x(k) (7) i_μ(k)= Cμ x(k) (8)

Here

A_d(L_μ) = _eA_c(L_μ)h (9) B_d(L_μ) = _eA_c(L_μ)s_ds B_c(L_μ) (10)

where h is the sampling interval. To calculate (9) and (10), it is assumed that L_μ, which is a function of the magnetizing current i_μ, is constant during the sampling interval.

A model of a primary circuit upon a short circuit, which is often described in the literature (see, e.g. "Digital protective relaying through recursive least-squares identification" by A Isaksson, published in IEE Proc. Pt. C, Vol. 135, No. 5, 1988, p. 441) describes the primary current as i_p(t) = d₀e-t/τ + ( c_nsin(ω_cnt + Φ_n)) (11) n=1 where τ is the time constant of the circuit, m is the number of harmonic components, ω_c is the fundamental frequency, and Φ_n the phase shift.

To describe the primary current during a short circuit, in, inter alia, "Estimation of the primary current in a saturated transformer", by K W Chen and S T Glad: Report LiTh-ISY-I 120, 1991, Department of Electrical Engineering, Linköping University, Linkόping, Sweden, the assumption is made that the exponential term is replaced by a Taylor series expansion. However, it is obvious to use a fixed, predetermined value of τ. This reduces the number of parameters to be estimated while at the same time the model becomes less sensitive to unmodelled harmonics and noise.

For a general continuous time signal i(t), the corresponding sampled value is denoted i(k), that is, i(k) =i(t) |_t=kh (12) where k is a running index.

It is part of the prior art that a sampled version of (11) via linear regression can be written in state-space form as i_p(k) = φ^Tθ(k) (13) θ(k + 1) = F θ(k) (14)

According to conventional recursive technique, T means transposition, θ is a parameter vector, and φ a regression vector. According to M S Sachdev, H C Wood and N C Johnson: "Kalman filtering applied to power system measurements for relaying", IEEE Trans, on PAS, Vol. PAS-104, No. 12, 1985 pp. 3565-3571, φ and F for the case m = 1 can be written as φ ^T = (1 1 0) (15) -^{h/τ 0 0}

F = cos(ω_ch) sin(G_ch) (16)

-sin(ω_ch) cos (ω_ch

Now, if the output signal i_p(k) from the signal model according to (13) and (14) is taken as input signal to the current transformer model according to (6), (7) and (8), the resultant system according to known technique (see, e.g., "Control System Toolbox", The Math Works, Inc, p. 2.149) can be written as z(k + 1) = A_t z(k) (17) i_p(k) = C_pt z(k) (18) i_s(k) = C_st z(k) (19) i_μ(k) = C_μt z (k) (20) where

(21)

C_pt = (φ^T 0 0) (22)

C_st = (0 0 0 C_s) (23) C_μt = (0 0 0 C_μ) (24 ) and where the used state vector is given by z(k) = (θ(k)^T x(k)^T)^T (25)

In similar model contexts it is common to use a Kalman filter to reconstruct the state vector z(k). The general method for such filtering is described, inter alia, in "Optimal Filtering", by B D O Anderson and J B Moore, Prentice Hall, 1979, p. 44.

If a fixed value of L_μ is used, A_t in equation (17) will become a constant matrix. It is then possible to utilize a stationary Kalman filter, which provides the reconstruction (k) = At (k-1) + K( i_s(k) - C_st A_t (k-1)) (26)

Here, "K" is a gain vector. One way of calculating "K" is shown in B Carlsson, "Digital Differentiating Filters and Model Based Fault Detection" , PhD thesis , Department of Technology, Uppsala University, Uppsala, Sweden, 1989 , pp . 72-76. In a report LiTh-ISY-I 120 from the Department of Electrical Engineering, Linköping University, Linköping entitled "Estimation of the primary current in a saturated transformer" , 1991, presented by K W Chen and S T Glad, a method is described based on a Kalman filter applied to a signal model which is based on Taylor series expansion and connected to a current transformer model . Otherwise , continuous time is used here with a time-varying Kalmangain, which per se requires a high numerical capacity. To detect saturation, a simple detector is used, based on the variance of the estimated fault . SUMMARY OF THE INVENTION

The invention relates to a numerical method which, via discrete-time measurements of the secondary current i_s(k) from a current transformer, reconstructs the primary current and the secondary current in case of saturation of current transformers and a method for detecting saturation as well as for detecting when saturation ceases. The reconstruction takes place with the aid of current transformer models of an unsaturated and a saturated current transformer as well as a signal model of the primary current, wherein the the current transformer models are fed from the signal model and saturation is decided by a decision logic unit when the output signal from the saturated current transformer model gives a better description of the secondary current than what the output signal from the unsaturated current transformer model gives. The output signals i_pr(k) and i_sr(k) in case of a decision about saturated current transformer consist of a reconstructed primary current i_pn(k) and a reconstructed secondary current i_sn(k) and, in case of a decision about unsaturated current transformer, of a reconstructed primary current i_pn(k) and an actually measured secondary current i_s(k).

The invention also relates to a device for carrying out the above-mentioned methods.

DESCRIPTION OF THE PREFERRED EMBODIMENTS

The invention is based on the current transformer model which is described by equations (1) and (2) and the general state-space model described by equations (3), (4) and (5).

The following state-space description is used in equations (3), (4) and (5) A_c(L_μ) = (27)

B_c(L_μ) = (28)

C_s = (0 1) (29)

C_μ = (1 0) (30)

Further, the invention is based on discretization according to equations (9) and (10), a sampled signal model according to equations (13) ... (16), and on the above-described combined signal and current transformer model which is described by equations (17) ... (25). According to the invention, two values of the inductance Lμ are used, which is denoted

L_μn = the nominal value of the magnetizing inductance L_μs = the saturated value of the magnetizing inductance This means that according to the invention two state-space descriptions or models are obtained, whereby also the matrix A_t according to equation (21) is given two values according to

0

Ant = B_d(L_μn)φ^T Ad(L_μn) (31) (32)Ast = B_d(L_μs)φ^T Ad^(Lμs⁾

Now, operating with two different values of the inductance L_μ means at the same time that two different values of A_d and B_d are obtained according to the following

(9n)

A_d(L_μn) = _eA_c(L_μn)h

(9s)

A_d(L_μs) = _eA_c (L_μs) h

(10n)

B_d(L_μn) = _eA_c(L_μn) s_ds B_c(L_μn)

(10s) eA_c(L_μs)s_ds B_c(L_μs)

Bd⁽Lμs^{) )} =

(27n)

A_c(L_μn) =

(27s)

A_c(L_μs) =

(28n)

B_c(Lμn)=

(28s)

B_c ((L_μs) =

The numerical method according to the invention completes a number of steps for each sampled measurement of the secondary current.

The different steps comprise reconstruction from an unsaturated model of a current transformer, reconstruction from a saturated model, detection of saturation, as well as decision about returning from a saturated to an unsaturated state.

When saturation is not detected, a parallel prediction takes place with both the unsaturated and the saturated current transformer model. With the unsaturated model, the reconstruction is performed from a Kalman filter with a pre-calculated gain K, with A_nt according to equation (31), and with the C-parameters according to equations (22) and (23) as follows: n(k) = A_nt (k-1) + K.g.(i_s(k) - C_st A_nt (k-1)) (33) i_pn(k) = C_pt n(k) (34) i_sn(k) = C_st n(k) (35)

Here, i_pn(k) is a reconstruction of the primary current and i_sn(k) is a reconstruction of the secondary current.

The invention is further based on the fact that the effect of i_s(k) on the state reconstruction (k) is shut off as soon as a saturation has occurred. This is done with the aid of the gain parameter "g" which, the first time a saturated state is detected, is changed from 1 to 0. This means that during saturation, the currents will be reconstructed based on measurements of the secondary current prior to saturation. To detect that saturation has occurred, a saturated model according to the invention can be used and the following relationships be utilized: (k) = (A_st)^d+1 (k - d -1) (36)

i_ss(k) = C_st (k) (37) with A_st according to equation (32). The variable "d" will be defined below after equation (40).

As is clear, the invention calculates and uses a state reconstruction (k) for the saturated model which is determined by the state reconstruction (k-d-1), that is, a state reconstruction obtained d+1 samples earlier.

According to the invention, the reconstruction error for the respective model is then calculated according to e_n = i_s(k) - i_sn(k) (38) e_s = i_s(k) - i_ss(k) (39)

In the following the invention will now be described with reference to a flow diagram according to Figure 2. The reconstruction method as described above and in accordance with equations (27) ... (39) takes place in a block B1, which is supplied with the sampled discrete-time secondary current i_s (k).

Next, a comparison is made between the absolute value of the reconstruction errors for the two models. According to the flow diagram in Figure 2, this takes place in block B2. If the reconstruction error calculated on the basis of the saturated model is smaller than the reconstruction error for the unsaturated model, that is, if

| e_s | < | e_n | (40) this indicates that a saturation has occurred. To prevent short disturbances from denoting a saturated state, it is required according to the invention that the inequality remains for a predetermined number (f) of consecutive samples for saturation to be determined. The control whether the saturation indication has existed during the above- mentioned number of samples is made by storing, in a variable "d" in the block B3, the number of consecutive times that equation (40) is true. This variable is set to zero if equation (40) is false. The first time equation (40) is true, that is, when d=1, the gain parameter g will be changed from g = 1 to g = 0, which is made in the block B4, which then restores this information to the block B1, equation (33), to prevent the distorted current signal from being allowed to influence the reconstruction. Now, if the counter for the variable d shows that equation (40) is true for a number of consecutive samples d > f according to block B5, a decision about saturation is made, which takes place in the block B6 where also the variable M is set equal to 0.

If the saturation indication determines that no saturation has occurred, block B7 will generate M = 1, which is passed to a block B8. To this block are supplied the sampled secondary current i_s (k), the primary current i_pn(k) reconstructed by block Bl as well as the reconstructed secondary current i_sn(k). By an M = 1 from block B7, block B8 will deliver a value of the primary current i_pr(k) which is equal to the reconstructed primary current i_pn(k) and a value of the reconstructed secondary current i_sr(k) which is equal to the measured secondary current i_s(k).

Now, if the saturation indication remains for the preset number of samples f, the method according to the invention will understand that a saturation, that is, M = 0 according to the flow diagram, has occurred. The output signal i_pr(k) from block B8 will then be represented by the reconstructed primary current i_pn(k) and the output signal i_sr(k) by the reconstructed secondary current i_sn(k). The iteration goes on continuously in the usual way by setting k = k + 1, shown in principle by block B9.

As will be clear from the above, according to the invention a preferred embodiment with two current transformer models is used, one representing an unsaturated model and the other a saturated model. Within the scope of the invention, several current transformer models, based on fixed magnetizing inductances with a value which lies between those corresponding to a nominal and a saturated value, may be used.

Within the scope of the invention, the signal model may also consist of harmonics according to equation (11).

Within the scope of the invention, a variant of an embodiment with regard to detection of saturation, as described by means of blocks B2, B3 and B4, may be replaced by another detection principle, for example a method described in SE 9100917-5.

It is, of course, also important to be able to detect when a saturation ceases. One embodiment is described by the inequality (40) and the block B2, respectively. Also in this case, similar detection principles, capable of replacing the contents of blocks B2, B3 and B5 , fall within the scope of the invention. As an example, the reconstructed magnetizing current i_μ, when saturation has been determined, may be stored. During saturation, the magnetizing current will be practically equal to the reconstructed primary current i_p(k). This means that the magnetizing current increases drastically during saturation. Information to the effect that saturation has ceased may be obtained when the magnetizing current decreases below the stored value. It is, of course, also possible to use i_s(k) directly to make a decision about restoration to an unsaturated state. A device for carrying out the method according to the invention preferably forms an integral part of a protective relay where current measurement takes place with the aid of current transformers. From the design point of view, the device may be arranged in the form of a number of blocks according to Figure 2 or constitute an integrated device.

In an integrated embodiment all equations and logical decisions are included in a computer-based program which has the discrete-time, sampled secondary current i_s (k) as input signal and whose output signals, upon a decision on an unsaturated state, consist of i_pr(k) = i_pn(k) and i_sr(k) = i_s (k) and whose output signals, upon a decision on a saturated state, consist of i_pr(k) = i_pn(k) and i_sr(k) = i_sn(k).

Claims

1. A method for reconstruction of saturated current transformer signals, characterized in that the reconstructed current transformer signals which comprise reconstructed primary current (i_pr(k)) and reconstructed secondary current (i_sr(k)) are obtained via sampled discrete-time measurement of the secondary current (i_s(k)) and with the aid of at least two current transformer models and a signal model of the primary current.

2. A method for reconstruction of saturated current transformer signals according to claim 1, characterized in that the reconstruction is obtained with the aid of current transformer models of a saturated and an unsaturated current transformer.

3. A method for reconstruction of saturated currrent transformer signals according to claim 1, characterized in that the current transformer models are fed from the signal model.

4. A method for reconstruction of saturated current transformer models according to claim 1, characterized in that saturation is decided by a decision logic unit when the output signal from the saturated current transformer model gives a better description of the secondary current than what the output signal from the unsaturated current transformer model gives.

5. A method for reconstruction of saturated current transformer signals according to claims 1 to 4 and wherein the reconstruction and the saturation decisions are based on a combination of the signal model and the unsaturated and the saturated current transformer model and wherein the method is characterized by the following relationships _n(k) = A_nt (k-1) + K.g.(i_s(k) - C_st A_nt -zn(k-1)) (33)

i_pn(k) = C_pt _n(k) (34) i_sn (k) = C_st (k) (35) where (k) is the state reconstruction based on an unsaturated current transformer model

F 0

A_nt = (31)

Bd(L_μn)φ^T Ad(L_μn)

e-h/τ 0 F = 0 cos(ω_ch) sin(ω_ch) (16) 0 -sin(ω_ch) cos(ω_ch)

B_d(L_μn) = _eA_c(L_μn)s_ds B_c(L_μn) (10n) ω_c is the frequency of the fundamental tone L_μn is the nominal magnetizing inductance φT = (1 1 0) (15)

A_d(L_μn) = e^Ac^(Lμn)h (9n) A_c(L_μn)=

(27n)

B_c(L_μn)= (28n)

R_fe = loss resistance in the iron core

R₂ = resistance of burden, winding, and conductors

L₂ = inductance of burden and leakage inductance

K is a predetermined gain to the reconstruction according to equation (33) g is a gain parameter which may assume the value 0 or 1 and which, upon a decision about unsaturated state, assumes the value 1 i_pn(k) is a reconstructed primary current i_sn(k) is a reconstructed secondary current

C_pt = (φ^T 0 0) (22)

C_st = (0 0 0 C_s) (23)

C_s = (0 1) (29) h = sampling interval τ = time constant of primary circuit and for detecting saturation, the following reconstruction, based on a saturated current transformer model, is used (k) = (A_st)^d+1

(k - d -1) (36) i_ss (k) = C_st (k) (37) where (k) is the state reconstruction based on a saturated current transformer model

A_st = (32)

B_d(L_μs)φ^T A_d(L_μs)

L_μs is the magnetizing inductance for a saturated current transformer

A_d (L_μs) = e A_c(L_μs)h (9s)

B_d (L_μs) = _eA_c (L_μs ) s_ds B_c (L_μs) (10s)

A_c(L_μs) = (27s)

B_c(L_μs)= (28s) after which the reconstruction error for the two models is determined according to e_n = i_s(k) - i_sn(k) (38) e_s = i_s(k) - i_SS(k) (39) and when | e_s | < | en | (40) is true, a variable "d" is adapted to store the number of consecutive true samples and when d = 1, g is changed from 1 to 0 and if d > f, where "f" is a predetermined consecutive number of samples, it is decided that saturation has occurred, whereby the measured secondary current i_s(k) is replaced by the reconstructed secondary current i_sn(k) and when equation (40) is no longer true, or if the predetermined number of consecutive samples "f" is not achieved, the current transformer is considered not to be saturated.

6. A device for carrying out the method for reconstruction of saturated current transformer signals according to claims 1-5, characterized in that the device comprises a block Bl which has the sampled discrete-time secondary current i_s (k) as input signal and which is programmed with the equations

A_d(L_μ) = _eA_c(L_μ)h (9)

B_d(Lμ) = _eA_c(L_μ)s_ds B_c(L_μ) (10) A_c(L_μ) = (27)

B_c(L_μ) = (28)

C_s = (0 1) (29)

C_μ = (1 0) (30)

A_nt = (31)

B_d (L_μn) φ^T A_d (L_μn)

A_st = (32)

B_d(L_μs)φ^T A_d(Lμ_s)

n(k) = A_nt (k-1) + K.g. (i_s(k) - C_st A_nt (k-1)) (33)

i_pn(k) = C_pt (k)

(34) i_sn (k) = C_st (k) (35)

(k) = (A_st)^d+1 (k - d -1)

(36) i_ss(k) = C_st (k)

(37) e_n = i_s(k) - i_Sn(k) (38) e_s = i_s(k) - i_ss(k) (39) and with corresponding equations (9n), (9s), (10n), (10s), (27n), (27s), (28n) and with variables included defined according to claim 5 and that the values of e_n and e_s obtained in block B1 are arranged as input values to a block B2 which is adapted to determine whether

| e_s | < | e_n | (40) is true, and if this is the case, information about this is adapted to be supplied to a block B3 which stores the number of consecutive samples that equation (40) is true in a variable "d" and the first time that d = 1, this information is adapted to be supplied to a block B4 which is adapted to change the variable "g" in equation (33) from g = 1 to g = 0, which change is adapted to be returned to block B1 and that the consecutive value stored in the variable "d" in block B3 is adapted to be supplied to a block B5 adapted to check whether d > f where " f" is a preset value, and if this is the case, information about this is adapted to be supplied to a block

B6 for decision that saturation has occurred, whereby a variable M is set equal to 0 and that if equation ( 40 ) is false, information about this together with information as to whether the inequality of block B5 is false are adapted to be supplied to a block B7 for information that no saturation has occurred, whereby the variable M is set equal to 1 and that the current value of the variable M is adapted to be supplied to a block B8 together with the sampled measured value i_s(k) of the secondary current, the primary current i_pn(k) reconstructed by block B1 and the reconstructed secondary current i_sn(k) and that when M = 0, block B8 is adapted to deliver a value i_pr (k) of the primary current which is equal to the reconstructed value i_pn(k) of the primary current and a value i_sr(k) of the secondary current which is equal to the reconstructed value i_sn(k) of the secondary current and that when M = 1, block B8 is adapted to deliver a value i_pr (k) of the primary current which is equal to the reconstructed value i_pn(k) of the primary current and a value i_sr(k) of the secondary current which is equal to the sampled measured value i_s (k) and that the iteration is adapted to continue by setting, in a block B9, k = k + 1.