WO2007140627A1

WO2007140627A1 - A method for determining the linear electrical response of a transformer, generator or electrical motor

Info

Publication number: WO2007140627A1
Application number: PCT/CH2006/000303
Authority: WO
Inventors: Martin Tiberg; Christoph Heitz; Olaf Hoenecker; Bjorn Gustavsen
Original assignee: Abb Technology Ag
Priority date: 2006-06-07
Filing date: 2006-06-07
Publication date: 2007-12-13
Also published as: CN101460856B; CN101460856A; BRPI0621703B1; ATE493668T1; DE602006019326D1; EP2024755A1; BRPI0621703A2; EP2024755B1; NO20090078L; US20090088995A1

Abstract

To characterise an electrical component (1), namely a medium or high voltage transformer, electrical motor or generator, a two-step procedure is carried out. In a first step, a set of terminal configurations are applied to the terminals (p1,..., pn) of the component (1) in order to obtain data describing the linear electrical response of the component 1 to any pattern of voltages uk or currents ik applied to the terminals (p1,..., pn). Typically, such data is e.g. expressed in terms of an admittance matrix Y or impedance matrix Z or, advantageously, by the set of current and voltage vectors (ik, uk). Using this data, the linear electrical response of the component 1 under a test terminal configuration can now be calculated in a second step. This procedure allows to determine the response under any desired test terminal configuration without the need to carry out the measurement under the test terminal configuration.

Description

A method for determining the linear electrical response of a transformer, generator or electrical motor

^■ DESCRIPTION

Technical field

The invention relates to a method for determining the linear electrical response of a transformer, generator or electrical motor as well as to a use of such a method.

Background

There are various methods to characterise the linear electrical response of transformers, generators or electrical motors. Typically, they include a measurement step under a certain terminal configuration, where the terminals of the component are connected to defined voltage or current sources or impedances or interconnected with each other. Some terminals may also be left uncon- nected (i.e. connected to an impedance of infinite value) . Then, a measurement is carried out, e.g. by determining the voltage at one terminal as a function of the voltage at another terminal .

A typical such technique is the Frequency Re- sponse Analysis (FRA) , which has evolved as a transformer diagnostic technique for quality control and detection of internal faults, which may not otherwise be detected without opening the transformer and inspecting it by hand and eye. Opening a transformer and emptying the oil therein is very costly and takes long time, whereas FRA is relatively cheap and done within less than one day.

BESTATIGUNGSKOPIE A FRA measurement consists of applying a voltage on one terminal of the transformer (the source terminal) and measuring the voltage output on one of the others (the sink terminal) for a large range of frequen- cies . The other terminals, which are not source or sink, can either be grounded or left open.

The diagnosis is carried out by studying how the voltage ratio between sink and source varies over the frequency range, and comparing these variations between the phases to check for asymmetries and/or comparing these variations to older records of the same or similar transformers to check for changes over time.

There is no standard specifying exactly how FRA measurements shall be done, even though groups both within IEEE and CIGRE work on such standardisation. Typically, the voltage ratio and sometimes the phase between a source and sink terminal is measured. The number of possible terminal configurations is very large because of the variations in the following parameters : - Type of transformer,

- Number of terminals on transformer,

-' Having the unused terminals open or grounded.

A typical delta-star or star-delta 2-winding transformer has seven terminals: 3 on the HV side, 3 on the LV side and one neutral point. Each of these 7 terminals can be the source terminal and any of the other 6 can be the sink terminal resulting in 42 different combinations and possible measurements. Having the unused terminals both open and grounded doubles the number of. meas- urements. Extra terminals or special configurations, such as having some terminals open and some closed, multiplies the number of measurements further. :

This large number of possible terminal configurations make a full characterisation of the component cumbersome if not impossible by standard FRA.. Therefore, typically, only a subset of the possible terminal configurations are measured, e.g.: - 3 on the delta side; connect any two terminals to measure over each winding,

- 3 on the star side; measure between neutral and each terminal, - 3 between HV and LV side for each phase.

These nine measurements often give a good picture of the winding conditions _> such that an assessment can be made of the transformer condition. However, there are. several important shortcomings of this kind of measurements:

1. It still takes a significant amount of time to reconnect nine times.

2. It happens easily that measurement errors are made when reconnecting. The contact resistance can vary greatly depending on how well the measurement cables are connected to the transformer terminals. Changing the position of the cables between two reconnections can significantly impact the measurement response in the high frequency range . 3. The measurement is incomplete. From just the LV-HV measurements per each phase, you do not see the coupling between the phases, which may reveal important information. 4. Most commercial FRA measurement equipment has 50 Ohm impedance both on the source and sink channel, which is added in series to the transformer measurement and damps the response for low impedances. Thus, impedance variations between say 1-3 Ohm may hardly be seen due to the 50 Ohm impedance in series. 5. Each of the 9. measurements is performed on a different state of the transformer. Consequently, the terminals are- alternatively open or terminated with 50 Ohm, which completely changes the transformer state. Even if all non-measured terminals are grounded, the termi- nation is still changed from ground to 50 Ohm on some terminals between each measurement, which invalidates advanced analysis and modeling that assumes that the transformer is a constant well-defined system.

Summary of the invention

Hence, it is an object to provide a method, which allows to determine the linear response of a transformer, generator or electrical motor having several ter- minals, i.e. at least two terminals, in particular at least three terminals, under a given test terminal configuration. The method should have improved ease of use and/or reliability.

This object is achieved by the method accord- ing to claim 1.

Accordingly, in a first step a) , a set of terminal configurations are applied to the terminals of the component in order to obtain data describing the lin- ear electrical response of the component to any pattern of voltages or currents applied to the terminals. Such data is e.g. expressed in terms of an admittance or impedance matrix or, advantageously, by a set of current and voltage vector pairs as described below. Using this data, the linear electrical response of the component under the test terminal configuration can now be calculated in a second step b) .

This procedure has the advantage that no actual measurement under the test terminal configuration is required. Rather, the measurement can take place under any suitable set of terminal configurations, which allows to choose the most suitable measurement process at hand. Advantageously, the same data can be used to calculate the response of the component to a plurality of different test configurations.

In a further advantageous embodiment, a measuring device is connected simultaneously to all terminals of the component. The measuring device is adapted to generate the set of terminal configurations and to measure the response of the component to each of these terminal configurations. For example, the measuring device may be equipped to apply different values of voltages, currents and/or impedances to each terminal. This allows to generate the set of terminal configurations without the need to change the cables attached to the component, which increases the accuracy of the measurement. In particular, such a measuring device can be operated automatically, which allows to increase the measurement speed and reliability.

Advantageously, the linear electrical response at the test terminal configuration is calculated as the ratio and/or phase shift between two voltages at different terminals, in particular as a function of frequency. This type of information is used in the so-called frequency response analysis (FRA) , which is applied when assessing the status or ageing of a transformer. The present method al- lows to carry out FRA even if no direct measurement of the ratio and/or phase shift between two voltages at different terminals was made.

A particularly advantageous implementation of step a) comprises an "estimation procedure" in which an estimated admittance matrix Y' is determined by applying voltages to the terminals of the component and measuring the response of the component. The estimation procedure can e.g. consist of a conventional measurement of the admittance matrix Y' by applying a voltage to one terminal, grounding all other terminals, measuring the current at each terminal, and repeating this procedure for all terminals. The estimation procedure is followed by a "measurement procedure" , in which several voltage patterns Uj₅- are applied to the terminals. The voltage patterns corre- spond to the eigenvectors v^ of the estimated admittance matrix Y⁷ , wherein "correspond" is to express that the pattern u^- is substantially (but not necessarily exactly) parallel to the (normalised) eigenvector v^- corresponding to each eigenvalue λ_]ς. For each applied voltage pattern U£, the response of the component is measured.

As it has been found, applying voltage patterns U]^ corresponding to the eigenvectors vj_j- of the admittance matrix Y' allows to obtain a more accurate description of the component, especially when the eigenvalues λjς of the admittance matrix Y' differ substantially from each other. The invention is particularly useful for high-voltage or medium-voltage components, i.e. for components suited for operation at voltages exceeding 1 kV. The method can e.g. be used for characterising the electrical component. In order to do so, a reference can be provided, e.g. measured at an earlier time (prior to step a) or measured on a reference component, which reference describes the response (the "first response" ) of the component under a given test terminal configuration. A measurement according to the present method is then carried out to determine the actual state of the component, and the data from this measurement is used to calculate a "second response" of the component under the test terminal configuration. The first and the second response are then compared for checking the actual status of the component.

Brief description of the drawings

Further embodiments, advantages and applica- tions of the invention are disclosed in the dependent claims as well as in the following description, which makes reference to the figures, wherein:

Fig. 1 is an example of a transformer to be characterised by the present method, Fig. 2 shows the component of Fig. 1 connected to a measuring device, Fig. 3 is a schematic illustration of a component to be characterised,

Fig. 4 is a block circuit diagram of a first embodiment of a device for a measuring device for charac- terising the component,

Fig. 5 is a second embodiment of a measuring device, and

Fig. 6 is a third embodiment of a measuring device .

Embodiments of the invention

1. Definitions:

The term "terminal configuration" refers to a defined state of all terminals of the component. The state of a terminal k can be defined by

- the current ij_j- that is flowing through it (ij_j- = 0 corresponds to an open terminal) or

- the voltage uj_j- that is applied to it (Uj₅- = 0 corre- sponds to a grounded terminal) or

- the impedance Zjς and a voltage φjς of a voltage source that are connected in series to it, or

- the index m of another terminal that the given terminal k is connected to (or a series of indices m.χ, m₂, ..., if terminal k is connected to several other terminals) .

2. A first example:

To illustrate the present invention, an exam- pie thereof is described in the following. This example relates to the characterisation and, in particular, to the quality control of a transformer 1 as it is e.g. shown in Fig . 1.

The transformer 1 of Fig. 1 is of the star- delta configuration and has n = 7 terminals Ll, L2, L3 , N, Hl, H2, H3. To characterise transformer 1, a two-step procedure is carried out, namely a measurement step a) and a calculation step b) .

In measurement step a) all seven terminals Ll, L2, L3, N, Hl, H2 , H3 of transformer 1 are connected to a measuring device 2 as shown in Fig. 2. Measuring device 2 comprises, for each terminal Ll, L2 , L3 , N, Hl, H2, H3, an adjustable voltage source and/or an adjustable current source and/or an adjustable impedance. Measuring device 2 automatically applies a set of terminal configurations to transformer 1 by repetitively adjusting the voltage sources, current sources and/or impedances. For each terminal configuration, the response of transformer 1 is measured, e.g. by measuring the voltages and cur- rents at all terminals Ll, L2 , L3 , N, Hl, H2 , H3. In a simplest embodiment, n different terminal configurations are applied, e.g. by applying n different, linearly independent voltage vectors u^ (i = 1 ... n) and measuring the corresponding current vectors i^ (i = 1 ... n) . Each such voltage vector u^ has n elements {u^±, ... v._n±) , and each current vector i^ has n elements {i-i±i — i-ni^ indicating the voltage u^ and current i_]ζj_ at terminal k in terminal configuration i. The knowledge of the vectors u^ and i^ for the n terminal configurations allows to estimate the admittance matrix Y from a set of n vector equations

±± = Y-U₁ (1)

The measurements are carried out as a func- tion of frequency.

It must be noted though that this type of measurement yields inaccurate results only, and a more refined measurement method is described in the section "improved measurement method" below. In general, the result of measurement step a) is data describing the linear electrical response of the component 1 (such as transformer 1) to any pattern of voltages or currents applied to the terminals Ll, L2 , L3, N, Hl, H2, H3. This data can e.g. be expressed by the admittance matrix Y, the corresponding impedance matrix Z, or by n linearly independent current and voltage vector pairs, each pair describing the voltages and corresponding currents at all terminals Ll, L2 , L3 , N, Hl, H2 , H3. For reasons explained under the section ^λλimproved measurement method" below, the n linearly independent current and voltage vector pairs is considered to be the most ad- vantageous representation of the data.

In the calculation step b) following measuring step a) the linear electrical response of transformer 1 is calculated for a given test terminal configuration. In general, this test terminal configuration will not be among the set of terminal configurations used in step a) for measuring the component .

For example, when FRA is to be carried out on transformer 1, the test configuration will be a configuration where all except two terminals are grounded (or open) . The two terminals are assumed to be terminated by known impedances, such as 50 Ohm. A voltage is applied to one of the two terminals, while the voltage at the other terminal is measured. In other words, the test terminal configuration corresponds to a typical FRA measurement terminal configuration. However, instead of carrying out an actual measurement at the test terminal configuration, the result of such a measurement is simulated by calculating the component's behavior from the data obtained in step a) . ^• The calculated result can then e.g. be used to determine the ratio and/or voltage phase difference between two terminals as a function of frequency to obtain a graph as used in FRA.

For example, to calculate the response of an FRA measurement, where the first terminal (e.g. Ll in Fig. 1) is the source, the fourth terminal (e.g. Hl in Fig. 1) is the sink, and all the other terminals L2,- L3 , N, H2, H3 are open, and the first and fourth terminal are terminated with a R_]_ = R4 = 50 Ohm impedance, the current vector i is given by

The FRA response U4/U1 is straight-forward to calculate from the relationship

with y±j being components of the admittance matrix Y. The easiest way to proceed is to multiply both sides of the equation (3) with the impedance matrix Z = Y"¹ and then to reduce the system of equations to just the first and fourth row. Then there are three unknowns, ui, U4 and i-^, and two equations. Since we are interested in the complex ratio U4/U1 (i.e. ratio of the voltage amplitudes as well as their mutual phase shift) , the unknowns are further reduced to two, and we get the following result:

H_{l =} . R₄ ^■ Z₄₁

(4)

U₁ ^₁₁(Z₄₄ + R₄) — Z₄₁Z₁₄

with z£j being the components of the impedance matrix Z. If the same FRA with the unused terminals L2 , L3, N, H2, H3 being grounded instead of being open were to be calculated, their terminal configurations would just change from i_x = 0 to U_x = 0, yielding the expres- sion

^U ₁ _ &. ^• y*i (5)

U₁ 1 + R₄₁ - y₄₄

Once the data is obtained in step a) , the re- sponse of the component 1 under any test configuration can be calculated. In particular, the response of the component 1 under several test configurations (e.g. under the nine configurations typically used for FRA) can be calculated easily and quickly.

3. General applications :

In the first embodiment above, the application of the present method for FRA measurements on a transformer 1 having seven terminals Ll, L2 , L3 , N, Hl, H2, H3 has been described. It must be noted, though, that the method can be used for numerous other applications. In particular:

Instead of using the method on a transformer 1, the component 1 under test may also be an electrical generator and/or an electrical motor.

In general, the number n of terminals pi, ..., P_n (Fig.3, 4) of the component 1 may vary, and e.g. be three (e.g. for a three-phase motor in delta configuration) or four (e.g. for a three phase generator in star configuration or for a one phase transformer) . The present method is in particular useful for components 1 with more than two terminals pi, ..., p_n, where there is a potentially large number of different terminal configurations.

The present method can be used for various purposes. A typical application is quality control, e.g. by using FRA as described above, or by simulating other measurements using the data obtained in measurement step a) . Another application is network modeling: Some network models require the measurement of the linear response of a component 1 under certain given terminal configurations - the present method can obviate the need to actually carry out these measurements by using the data from measurement step a) in order to calculate the response under the given terminal configurations .

4. Improved measurement method:

This section describes an improved method for obtaining the data in measurement step a) .

Fig. 3 shows a multi-terminal component 1 having n > 1 terminals p^ through p_π, which may be a transformer, electrical motor or generator. When linear voltages u^ through u_n are applied to the terminals pi through p_π, currents i_]_ through i_n will flow. The linear electrical response of component 1 is characterised by its admittance matrix Y or, equivalently, by its imped- ance matrix Z. In admittance notation, applying the voltage vector u = (U_]_ ... u_n) of voltages -at the terminals Pl through p_π generates a current vector i = (i^ ... I_n) as follows :

i = Y-u. (11)

The general principle of the improved measurement method is based on an estimation procedure and a measurement procedure. In the estimation procedure, an estimated admittance matrix Y' is determined, in the measurement procedure a more accurate measurement is carried out.

In the estimation procedure, the elements of the estimated admittance matrix Y' can e.g. be measured directly using conventional methods. The diagonal elements Y'±± can e.g. by measured by applying a voltage u± to terminal p-£ and measuring the current i± at the same terminal p^ while all other terminals are short-circuited to zero volt, i.e. Y'±± = i±/u-± while UJ - 0 for i ≠ j. The other elements Y'±ή of the matrix can be measured by applying a voltage u-[ at terminal p^ while setting all other terminals to zero volt and measuring the current ±ή at terminal p-y, Y'±ή = ij/^u± while Uj = 0 for ± ≠ j.

Other conventional methods for measuring the estimated impedance matrix Y' in the estimation procedure can be used as well. In general, the estimated admittance matrix

Y' has n eigenvalues λ-± ... X_n and n corresponding to (normalised) eigenvectors vj_ ... V_n for which

^γ/ 'Vk = Λfc-vfc. (12)

Once the estimated admittance matrix Y' is known, its eigenvectors v/_j- can be calculated.

In a measurement procedure following the estimation procedure, several (in general n) voltage pat- terns u^- = (u^ ... u_n]ς) are applied to terminals p_]_ ... P_n of component 1. Each voltage pattern u^ corresponds to one of the eigenvectors V£. For each applied voltage pattern uji_j-, a response of the component 1 is measured, in particular by measuring the induced current pattern ij_ς. As mentioned above, voltage pattern Uj₅- corresponds to the (normalised) eigenvector vj^ (which is one of the n normalised eigenvectors of the admittance matrix) , namely in the sense that the voltage pattern uj_j- is substantially parallel to the eigenvector vj^- correspond- ing to eigenvalue A^. Theoretically, using Uj₅- °= Vj₅- would be the best solution, but a measuring device 2, 3 generating the voltage patterns Uj₅- will, in general, not be able to generate voltage patterns matching the eigenvectors vj_j- exactly due to discretisation errors. Methods for handling devices 2, 3 with limited resolution for generating the voltage patterns will be addressed below. Once the measurement procedure is complete, the voltage patterns u^ and the corresponding current patterns ij_j- , i.e. a set of n voltage and current vector pairs u^-, i^-, fully characterises the linear response of the component 1.

In general, the admittance matrix Y is frequency dependent. Hence, in many applications, the linear electrical response of component 1 should be known for an extended frequency range, e.g. from 50 Hz to more than 500 kHz. For this reason, the estimation procedure is carried out at a plurality of frequencies ύ)± in the given frequency range.

Advantageously, for each estimation procedure, the eigenvalues X]^[COj) at the given frequency co± are calculated. Then, the most critical frequencies are determined, which are those frequencies where the eigenvalues λfciCD-i) reach a local maximum or minimum or, in particular, where the absolute ratio between the largest and smallest eigenvalue has a maximum or exceeds a given threshold. These critical frequencies are of particular interest, either because they are indicative of a resonance of component 1 or because they show that some of the estimated eigenvalues may be of poor accuracy and the described measurement procedure is required to increase the accuracy.

It is principally possible to divide the desired frequency range into a number of frequency windows and to calculate the most critical frequencies in each frequency window. For each or at least some of the critical frequencies, the measurement procedure described above is carried out to refine the measurement. In addition or alternatively thereto, the measurement procedure can be carried out at other points within the frequency range of interest.

The frequencies 0)± where measurements are carried out can be distributed linearly or logarithmi- cally over the range of frequencies of interest. In. an advantageous embodiment, though, the density of measurement frequencies co± close to the critical frequencies as mentioned above is larger than the density of measurement frequencies ύ)± in spectral regions far away from the critical frequencies. This allows to obtain a more reliable characterisation of the component 1.

4.1 The measuring device :

A general measuring device 2, 3 for carrying out the improved measurement method is disclosed in Fig. 4. In a most general case, measuring device 2, 3 comprises n adjustable voltage sources generating volt- ages φi to φ_n, which are fed to the terminals pi to p_n through impedances Z^ to Z_n. The voltages ψi to φ_n all have equal frequency and known phase relationship. The impedances Z-^ through Z_n may be practically zero or, as described below, they may be adjustable and potentially non-zero. A control unit 3 is provided for automatically adjusting the voltage sources and, where applicable, the impedances Z-± to Z_n.

For the device of Fig. 4 we have

φ = u + D-i, (13)

where φ = (φ^ ... φ_n) axe. the voltages of the voltage sources, u = (ιzχ ... U_n) the input voltages at the terminals P₁, ..., P_n, and D is a diagonal matrix with the diago- nal elements Z^ to 2_π.

Combining equations (11) and (13) gives the following relationship between the input voltages u and the applied voltages φ:

^■ u = (I + D-Y)^"1-^). (14)

where I is the π X Ji identity matrix. As mentioned above, the applied voltages u should correspond to the eigenvalues v^ of the estimated admittance matrix Y'. In general, however, it will not be possible to match this condition exactly, because the voltage sources will not be able to generate any arbitrary voltage values, but only a discrete set of values. If the number of voltage values that can be generated is small, the impedances Z-^ to Z_n can be designed to be adjustable as well in order to obtain a larger number of different input voltages u.

The input voltage vector u^ can be expressed as a linear combination of the eigenvalues v^, i.e.

i=l

with coefficients cc± . Combining equations (15), (11) and (12) yields

n i = ]jT λ±a_jy_jL . (5) i=l

Hence, in order to maximise the influence of the k-th eigenvalue λfc on the input current vector i in proportion to the other eigenvalues, the following error function must be minimised

In other words, for each eigenvalue X]^₁ the coefficients a^ ... O_n must be found (among the set of possible coefficients, which is a finite set due to the discretisation inherent to measuring device 2) for which the term of equation (16) is smallest. If measuring device 2 has adjustable voltage sources and impedances as shown in Fig. 5, we have

α = [ T"i V_n]-¹ • (I + D • Y-) ,-1 φ (17)

A measuring device 2, 3 for carrying out the above method should, in general, comprise π voltage generators 10 that are programmable to apply the voltage pattern u to the n terminals of the component 1 undre test. Further, it should comprise n current sensors 11 to measure the currents i. It should be adapted to apply at least n suitable voltage patterns u to the terminals pi, ..., p_n consecutively for measuring the linear response of the component 1 automatically. This is especially advan- tageous for components 1 having more than two terminals Pi/ •••/ Pn/ because using this kind of automatic measurement on components 1 with n > 2 terminals pi, ..., p_n provides substantial gains in speed and accuracy while reducing the costs . Advantageously, the measuring device 2, 3 should comprise a control unit 3 for carrying out the measurement using the estimation and measurement procedures outlined above.

One possible embodiment of a measuring device 2, 3 is shown in Fig. 5. In this device 2, 3, a voltage generator 10 for generating an individual voltage φ± of adjustable amplitude and phase is provided for each input terminal pi, ..., p_n. It also comprises n current sensors 11, one for measuring the current to/from each terminal P₁, ..., p_n. Control unit 3 is able to set the applied input voltage u directly by controlling the voltage generators 10. If the number of voltage values that can be generated by each voltage generator 10 is small, an optimum voltage for a given eigenvector Vj₅- can be calculated by minimis- ing the term of equation (16) . For each applied voltage pattern U£, control unit 3 measures the currents i through the terminals pi, ..., p_n by means of the current sensors 11.

Another possible embodiment of a measuring device is shown in Fig. 6. This device comprises a single voltage source 4 only. The voltage φ from the voltage source 4 is fed to n voltage converters 5 controlled by control unit 3, the voltage source 4 and voltage converters 5 being used instead of the voltage generators 10 of the previous embodiment. Each voltage converter 5 selec- tively connects one terminal P₁, ..., p_n to either the voltage φ directly, to the voltage φ through a damping circuitry 6, to ground via an impedance 7, to ground directly, or leaves the terminal pi, ..., p_n open (infinite impedance) . This measuring circuit has the advantage that it requires a single voltage source 4 only. Suitable settings of the voltage converters for each value can be calculated form equations (16) and (17) .

Further processing of the results:

As mentioned above, the described improved measurement procedure yields, for a given frequency, a set of voltage patterns u^ and the corresponding current patterns i^, which fully characterise the linear response of component 1 at the given frequency.

The values u^- and i^ for k = 1 ... n can, in principle, be converted into a more accurate estimate of the admittance matrix Y or the corresponding impedance matrix Z. However, if the smallest and largest eigenvalues of admittance matrix Y differ by several orders of magnitude, such a matrix Y is difficult to process numerically with floating point calculations due to rounding errors and limited accuracy of the numerical algo- rithms. Hence, in an advantageous embodiment of the present invention, the vector pairs u^- and i^ are used di- rectly for further processing, without prior conversion to an admittance matrix Y or impedance matrix Z.

For example, the results of the measurement procedure can e.g. be used for FRA as described above or for modeling the electrical properties the component 1 or of a network that component 1 is part of. Such a model can e.g. be used to analyse the stability of the network in general or its response to given events in particular.

List of reference numbers

1 component under test

2 measuring device 3 control unit, part of measuring device

4 single voltage source

5 voltage converter

6 damping circuit

7 impedance 10 voltage generator

11 current sensor

Claims

1. A method for determining a linear electrical response of a component (1) under at least one test terminal configuration, wherein said component (1) is a transformer, generator or electrical motor comprising several terminals (pi, ..., p_n) , said method comprising the steps of step a) applying a set of terminal configura- tions to said terminals (pi, ..., p_n) to obtain data descriptive of the linear electrical response of said component (1) to any pattern of voltages or currents applied to the terminals (pi, ..., p_n) , wherein the set of terminal configurations does not comprise the test terminal con- figuration, and step b) calculating the response under the test terminal configuration from said data.

2. The method of claim 1, wherein step b) comprises the calculation of the response of the compo- nent (1) under a plurality of different test terminal configurations .

3. The method of any of the preceding claims, wherein said step a) comprises the step of simultaneously connecting a multitude of the terminals (pi, ..., p_n) , in particular all of the terminals (pi, ..., p_n) , of said component (1) to a measuring device (2, 3), which is adapted to generate said set of terminal configurations, and to measure the response of said component (1) to said terminal configurations .

4. The method of claim 3, wherein said measuring device (2, 3) comprises, for each terminal Cp₁, ...,

Pn) , an adjustable voltage source (10; 4, 5) , and/or an adjustable current source, and/or an adjustable impedance (Z₁, .„, Z_n) , wherein said set of terminal configurations is generated by adjusting said voltage sources (10; 4, 5), current sources and/or impedances (Zi, ..., Z_n) , respectively. 5. The method of any of the claims 3 or 4 wherein said voltage sources (10; 4,

5), current sources and/or impedances (Z₁, ..., Z_n) are adjusted automatically under control of said measuring device (2, 3) .

6. The method of any of the preceding claims, wherein said step b) comprises the step of calculating a voltage ratio and/or voltage phase difference between two different terminals (p_1# ..., p_n) .

7. The method of claim 6 , wherein said step b) comprises the step of calculating a voltage ratio and/or voltage phase difference between two different terminals (p_x, ..., p_n) as a function of frequency.

8. The method of any of the preceding claims, wherein said component (1) has more than two terminals

9. The method of any of the preceding claims, wherein said component (1) is a high-voltage or medium- voltage device (1) .

10. The method of any of the preceding claims, wherein said data comprises a set of N linearly independ- ent current and voltage vector pairs i.]_ς, U£, each pair describing the voltages and corresponding currents at said terminals (pi, ..., p_n) , wherein N is the number of terminals (pi, ..., p_n) being measured of said component (1) .

11. The method of any of the preceding claims, wherein said data is descriptive of the linear electrical response of said component (1) over a frequency range between less than 100 Hz and more than 500 kHz.

12. The method of any of the preceding claims, wherein said step a) comprises an estimation procedure comprising the step of determining an estimated admittance matrix Y' of said component (1) by applying voltages to said terminals (pi, ..., p_n) and measuring a response of said component (1) and a measurement procedure comprising the step of applying several voltage patterns U^- to the terminals (pi, ..., P_n) of said component (1) , each voltage pattern u^ corresponding to an eigenvector V£ of said estimated admittance matrix Y' , and determining, for each applied voltage pattern u^, a response of said component (1) .

13. The method of claim 12 wherein said voltage patterns Uj₅- are generated by means of a measuring de- vice (2, 3) capable of applying a discrete set of different voltage patterns u^- to said terminal (pi, ..., p_n) , wherein each voltage pattern u/_j- corresponds to that member of said set that has the property that the term

n i=l

[λ_ka_kf

is minimal, wherein λ^ __ X_n are n eigenvalues of the estimated admittance matrix Y' and

n U^ = ∑ Of₁V₁ i=l

with coefficients Of^.

14. A use of the method of any of the preceding claims for characterising said component (1) by the steps of providing a reference describing a first response of said component (1) at said test terminal configuration, using said step a) to measure an actual state of said component (1) , and using said step b) to calculate a second response at said test terminal configuration, and comparing said first and said second response.

15. The use of claim 14, wherein said reference was derived from a measurement carried out prior to said step a) .