Next Article in Journal
Compact UWB MIMO Antenna for 5G Millimeter-Wave Applications
Next Article in Special Issue
Configuration of the Geometric State of Railway Tracks in the Sustainability Development of Electrified Traction Systems
Previous Article in Journal
Design of Hardware Accelerators for Optimized and Quantized Neural Networks to Detect Atrial Fibrillation in Patch ECG Device with RISC-V
Previous Article in Special Issue
CNN-LSTM vs. LSTM-CNN to Predict Power Flow Direction: A Case Study of the High-Voltage Subnet of Northeast Germany
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pearson Correlation in Determination of Quality of Current Transformers

1
Elektrogen d.o.o., 10000 Zagreb, Croatia
2
Power Engineering Department, Faculty of Electrical Engineering, Computer Science and Information Technology, University of Osijek, 31000 Osijek, Croatia
3
Faculty of Dental Medicine and Health, University of Osijek, 31000 Osijek, Croatia
4
SFB Education Center for Technology and Management, Higher Technical School, 8953 Dietikon-Zürich, Switzerland
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(5), 2704; https://doi.org/10.3390/s23052704
Submission received: 17 January 2023 / Revised: 20 February 2023 / Accepted: 24 February 2023 / Published: 1 March 2023
(This article belongs to the Special Issue Sensors for Measurements and Diagnostic in Electrical Power Systems)

Abstract

:
The article elaborates on the accuracy of current transformers (CT) in interaction with temperature and frequency using Pearson’s correlation. The first part of the analysis compares the accuracy of the mathematical model of the current transformer and the result of the measurement on the real CT using the Pearson correlation calculation. The mathematical model of CT is determined by deriving the formula of the functional error with the display of the accuracy of the measured value. The accuracy of the mathematical model is affected by the accuracy of current transformer model parameters and the calibration characteristic of the ammeter used to measure the CT current. Variables that cause deviation in the accuracy of CT are temperature and frequency. The calculation shows the effects on accuracy in both cases. The second part of the analysis refers to the calculation of the partial correlation of three quantities: (1) CT accuracy, (2) temperature, and (3) frequency on a set of 160 measurements. First, the influence of temperature on the correlation of CT accuracy and frequency is proven, following the proof of the influence of frequency on the correlation of CT accuracy and temperature. In the end, the analysis is combined by comparing the measured results of the first and second part of the analysis.

1. Introduction

This paper considers the process system identification methodology applied to the current transformer [1]. The theory of identification of process systems uses two basic methods: non-parametric and parametric identification procedures [2].
  • Non-parametric methods are performed in the time domain. Fourier analysis and correlation analysis are used in that case.
  • Identification using the parameter estimation procedure uses static variables. This method is characterized by the need to evaluate the initial model, which must be performed in order to obtain convergence towards the final solution. The process is, finally, described by a transfer function.
The specificity of the experiment is the use of mathematical procedures of Pearson correlation of measured and calculated values to obtain the most accurate initial model possible. The analysis highlights anomalies detected during statistical analysis of input and output quantities. The partial correlation calculation includes the third quantities as functions of disturbance, namely temperature and frequency.
The accuracy of measuring transformers is defined by accuracy class, transmission ratio, and rated power. These are numerical quantities. The error of the measuring transformer is a variable quantity that depends on the secondary load. The class index represents the maximum allowable error of the instrument when used in reference conditions. The error limits can be expressed in percentages of the scale length, the true value, or, most often, in percentages of the maximum value of the measuring range. The accuracy classes of current transformers are 0.1; 0.2; 0.5; 1; 3; and 5. Some omissions and approximations were made in the derived model. Physical phenomena such as hysteresis losses and eddy current losses are neglected in the model due to the small estimated amount. Frequency and temperature affect the amounts of core and winding losses. This is the reason for using the Pearson correlation as a statistical method [3,4]. As a result of the analysis, partial correlation factors are presented in the article as evidence of the interaction of three variables: CT accuracy, frequency [5,6], and temperature. The correction amounts obtained by partial correlation can be used as a justification for neglecting these influences, depending on the accuracy requirements [7].

2. Methodology of CT Accuracy Analysis Using Pearson Correlation Factors and Regression Factors

In the analysis of the measured data, we will use the regression factor “r2” and the correlation factor “r”, terms used in statistics and error theory. There is an interrelation between the regression factor “r2” and the correlation “r”, and the definitions of these two quantities should be distinguished because these values are related, but their interpretations are different:
  • The amount “r2” indicates how well the regression line approximates a certain set. It is a percentage value (%) that determines how much of the observed set is within the given variation. It ranges from 0 to 1, i.e., from 0 to 100%.
  • The correlation factor “r” is used when one wants to prove the relationship between two variables and the strength of that relationship. It is used with a set of two quantities and ranges from −1 to +1.
The analysis of the two quantities em and ef, i.e., the functional error of the CT and the measured error of the CT, are shown below. Those two quantities are dependent variables that depend on the percentage (%) of CT load. If we consider them as variables that describe the real and ideal measuring system, we can write them as dependent variables on the ordinate axis and the abscissa axis. They are not equal, but they are similar in terms of amounts, and we assume their connection, which we will try to prove statistically.
The goal of the analysis is to detect deviations and anomalies in the entire measurement area of the CT and to prove it by measurement. The diagram in Figure 1 graphically interprets data listed in Table 1:
The ideal characteristic represents the direction in which real measurements match the calculation model, i.e., functional amounts, 100%. The closer the real measurements are to the ideal characteristic cause the higher the quality of the computational model. Measurements are performed at a frequency of 100 Hz.
The formula of the measured characteristic is:
  y ^ = b 0 + b 1 · x = 0.373 + 1.292 · x
b0 = −0.373; b1 = 1.292
The formula of the ideal characteristics is:
y ^ = b 0 + b 1 · x = x
b0 = 0; b1 = 1
It can be seen that the ideal characteristic, b1 = 1, which, considering Formula (3), means that the correlation must be r = 1, and the standard deviations sx and sy must have the same value. For the real characteristic, some deviation results are shown in Formula (1).

Trajectory Obtained Using “Spline” Interpolation

A specific curve in Figure 2 is obtained by connecting the points by “spline” type interpolation in the order in Table 1. The independent variable is the primary CT current, which increases linearly from 0 to 100% with a step of 10%. Two functions with dependent variables em and ef were used to graphically display the measure of deviation from the ideal characteristic, one on the ordinate axis and the other on the abscissa axis. This way of displaying results was applied in the rest of the work to all experiments.
The benefit of “spline” interpolation:
  • visualization of the measurement characteristics of the current transformer
  • display of deviations on the entire measuring area
  • determination of the place on the curve that contributes to the deterioration of the correlation factor
  • insight into the sensitivity of trajectory
  • insight into the slope of the trajectory
Figure 3 displays the instruments used for this experiment. The parts used for the test are personal computer, current injection device, current transformer, and high precision measuring ammeter.

3. Pearson Correlation—Application to the Test

Correlation represents the mutual relationship between different sets of data represented by variables. Data sets can be stochastic or determined by a function. In this case, these are the percentage amounts of two errors: complex function errors ef and measured CT errors em.
The correlation coefficient is determined by the Formula (5). The value of the correlation coefficient ranges from “+1” to “−1”, i.e., from perfect positive correlation to perfect negative correlation. The correlation coefficient is based on the comparison of the interaction of two variables in relation to the maximum possible influence of the two variables. This correlation coefficient is also called the Pearson correlation coefficient [3,4].
r = ( X i X ¯ ) ( Y i Y ¯ ) i n ( X i X ¯ ) 2 + i n ( Y i Y ¯ ) 2
The use of correlation in measurement accuracy should be carried out with some limitations. The paper considers the relationship between function points. These function points can be declared as a set of stochastic points, so the correlation of these two sets of declared stochastic points is calculated according to the Formula (5).
Correlation does not determine measurement accuracy and error limits, which is shown in the following example.
Task: Calculate the Pearson correlation coefficient between two hyperbolas defined by the equations y1(x) = a/x + b and y2(x) = c/x + d on a limited interval.
Excerpt:
Two following equations:
y 1 ( x ) = a · f ( x ) + b
y 2 ( x ) = c · f ( x ) + d
are equated by using the common function f(x). It follows:
y 1 ( x ) b a = y 2 ( x ) d c  
y 1 ( x ) = y 2 ( x )   a c   + ( b d a c )
A linear dependence of y1(x) and y2(x) was obtained, which represents the direction of regression and assumes a perfect correlation of “+1” if the coefficient a/c is positive. Otherwise, the correlation is “−1”. Accordingly, two functions with different amounts can have a correlation coefficient of +1 i.e., perfect correlation.
The reverse is also valid: if two different functions y1(x) and y2(x) have a “+1” correlation, then they can be written in the form of linear equations and have a common function f(x).
We can conclude that the correlation is perfect, i.e., equal to “+1” for the entire family of curves:
y ( x ) = a x + b

4. Functional Error Analysis of the Current Transformer Transfer Ratio

The replacement model of the current transformer is shown in Figure 4. Table 2 gives the factory testing results of current transformer.
The variable transformation procedure can be seen in Table 3. The result is the curve Lm = f (Ict) given in Figure 5 defined by Formula (11) obtained by the interpolation of points in the transformation in Table 3. This equation will be needed to determine Lm as a function of current. It should be pointed out that CT is analyzed with external load Zb = 0. Information about the protective current transformer is given in Table 2. The material used for magnetic core was “Grain Oriented Silicon Steel Strips; Grade VM 97–27”.
Figure 6 is connected with Table 3. Variables Es and Ie are manually inserted in Table 3. Figure 6 and Table 3 were necessary to obtain Formula (11), which will be used for calculation functional error of CT (27).
A curve is obtained by the inclusion Lm = f(Ict):
Following the formula obtained using the 4th degree polynomial approximation from the data processing program:
Lm = −0.0438·Ict4 + 0.2079·Ict3 − 0.3636·Ict2 + 0.3834·Ict + 0.031

5. Derivation of the Current Transformer Functional Error

This part of the analysis refers to the determination of the functional fault on the secondary current fault. The amounts of resistance Rh and Reddy in Figure 1 are many times higher than the inductive resistance XL = jωLm and can be treated as infinitely large, i.e., ignored. The reduced model of the current transformer for Rh>> and Reddy>> is shown in Figure 7:
Below is a derivation of the error function of an unloaded current transformer with a short-circuited secondary circuit:
e = I c t I s t I s t ·
It follows:
I e x I c t = Z c t Z e x    
I c t = I s t I e x = I s t I c t Z c t Z e x
Z e x = j · ω L m ;   Z c t = R c t + j · ω L s
I c t = I s t Z e x Z e x + Z c t = I s t j ω L m j ( ω L m + ω L s ) + R c t
I c t = I s t   j · ω   L m R c t + ω 2 L m 2 + ω 2 L m L s R c t 2 + ω 2 ( L m + L s ) 2  
I c t = I s t [     j ·   ω   L m R c t R c t 2 + ω 2 ( L m + L s ) 2   +   ω 2 L m 2 + ω 2 L m L s R c t 2 + ω 2 ( L m + L s ) 2   ]
I c t = I s t [   ω   L m R c t R c t 2 + ω 2 ( L m + L s ) 2   ] 2 + [   ω 2 L m 2 + ω 2 L m L s R c t 2 + ω 2 ( L m + L s ) 2   ] 2
The following is the expression for the amplitude of the functional error of the measurement of a current transformer with a short-circuited secondary, i.e., for Zb = 0:
e = ( I c t I s t I s t   ) · 100 % = ( I c t I s t 1     ) · 100 %
e = [ [   ω   L m R c t R c t 2 + ω 2 ( L m + L s ) 2   ] 2 + [   ω 2 L m 2 + ω 2 L m L s R c t 2 + ω 2 ( L m + L s ) 2   ] 2 1 ] · 100 %
With:
  • Lm = secondary inductance Lmain of the current transformer, dependence on the current is given by (1)
  • ω = circular frequency
  • Rct = ohmic resistance of the secondary
  • Ls = inductance of the secondary connection lines
Based on the previous formula, the display of the CT error for the frequency of 50 Hz follows. The largest amount in Figure 8 is −1.6%. We notice that the error function is shifted towards the negative part of the scale.
Manufacturers of current transformers raise the error diagram of CT for the desired value by correcting the number of windings due to the requirement for better accuracy.

6. Correction of the Number of Secondary Windings

By correcting the number of secondary windings, the curve in Figure 9 is raised, the maximum error on the interval (0–100%) of In is reduced, and the expression for the new corrected secondary current is obtained:
I S T = I S T n n k
e = I ct I st I st · 100 % = ( n · I ct n k · I st 1 ) · 100 %
Inserting the following expression for the CT error:
e = [ n n k · [   ω   L m R c t R c t 2 + ω 2 ( L m + L s ) 2   ] 2 + [   ω 2 L m 2 + ω 2 L m L s R c t 2 + ω 2 ( L m + L s ) 2   ] 2 1 ] · 100 %
With:
  • Lm = secondary inductance Lmain of the current transformer, and dependence on the current is given by (11)
  • ω = circular frequency
  • Rct = ohmic resistance of the secondary
  • Ls = inductance of the secondary connection lines
  • nk = corrected number of windings
  • n = number of secondary windings
Figure 9 shows the CT diagram with the number of corrected windings. The number of windings was reduced from 100 to 99. This reduces the maximum measurement error for nominal values to 0.74%, and the interval within which the error is located is symmetrical around the abscissa of axis X. The comparation can be performed with a factory testing result of CT, which is displayed in Figure 10 for three specific loads.

7. Calibration of the Measuring Instrument

The measuring instrument can cause inaccuracy of the measurement result, so it is necessary to calibrate it. The CT current measurement system consists of a certified current source, a k = 100/1 transmission ratio current transformer, a class 10P10, and a measuring instrument, Fluke 86. It follows that:
  • A; constant current source (Amp)
  • B = (A/k) · (1 − ε1); “ε1” is functional error of the current transformer, “B” is secondary current of CT, “k” is current ratio
  • C = B · (1 − ε2); “ε2” measuring instrument error, “C” is current displayed on ammeter
By including the above we obtain:
  • C = (A/k) · (1 − ε1) · (1 − ε2) = (A/k) · (1 − ε1ε2 + ε1 · ε2) (Amp)
Follows a secondary current of CT corrected by the complex measurement error:
C ≈ (A/k) · (1 − ε1 − ε2)
with the assumption ε1, ε2 <<, C is measured secondary current.
The Fluke 86 ammeter was calibrated for the specified current and frequency. The goal was to obtain the complex error of the measuring system. This error consists of the functional error of the current transformer and the calibration error of the measuring device. Figure 11 shows a set of calibration curves for four frequencies 50 Hz, 100 Hz, 150 Hz, and 200 Hz measured at room temperature 26.3 °C.

8. Functional Dependence of CT Measurement Error on Temperature

By including the temperature in the expression for the functional error (24) of the current transformer, the error values are obtained depending on the temperature. The formula for the temperature coefficient of electrical resistance is:
Rct(T) = R0 · [1 + α · (TT0)]
with:
  • α = temperature coefficient of electrical resistance, α = 0.00386 1/K
  • T0 = initial temperature, 20 °C
  • R0 = electrical resistance at temperature T0, R0 = 0.445 Ω
In the formula for the CT error (27), the resistance is predicted to change in the variable Rct, so the temperature values for currents from 0.1 A to 1 A are included. The expression for the error is defined as follows:
e = [ n n k · [   ω   L m R c t R c t 2 + ω 2 ( L m + L s ) 2   ] 2 + [   ω 2 L m 2 + ω 2 L m L s R c t 2 + ω 2 ( L m + L s ) 2   ] 2 1 ] · 100 %

9. Temperature 26.3 °C—Results

Figure 12 shows the diagrams of complex function error, measurement error, and error difference.

10. Functional Dependence of CT Measurement Error on Frequency

In the formula for CT error (28), the frequency is predicted to change in the variable w, so the frequency values are 50 Hz to 200 Hz. The current range ranges from 0.1 A to 1 A. The expression for the error is, in this case [5,6]:
e = [ n n k · [   ω   L m R c t R c t 2 + ω 2 ( L m + L s ) 2   ] 2 + [   ω 2 L m 2 + ω 2 L m L s R c t 2 + ω 2 ( L m + L s ) 2   ] 2 1 ] · 100 %
By correcting the number of windings, a symmetrical distribution of the curves around the abscissa axis X was obtained. The error ehz ranges from −0.65% to 0.74%.

Frequency 50, 100, 150, 200 Hz—Results

Follows a presentation of the obtained results for frequencies 50 Hz, 100 Hz, 150 Hz, and 200 Hz. The values obtained from the model are compared with the measured values, and the data defining the accuracy and quality of the model parameters are shown. Figure 13 shows the diagrams of complex function error, measurement error, and error difference.
Figure 14 represents four curves of function error ef and measured error em. Differences between curves can be seen depending on frequency [8,9].

11. Discussion about Results of Measurements

The curves from Figure 14 are divided into 9 segments shown in Table 4 using the result of derivation (12), the analysis of the obtained results was performed. Each segment is correlated with an ideal characteristic. The sign of the curve slope is marked with + or −. It determines the correlation of the segment of the measured curve with the segment of the ideal characteristic. Negative results are marked as “−1”. On these parts, the proposed model behaves differently from the measured model, i.e., the correlation is negative. Since the example that considers frequency dependence is taken, it is possible to refine the CT model and cancel negative correlations on all segments. The parts marked as “−1” point out the imperfection of the CT or measuring instrument, but also the shortcomings of the model. The largest correlation oscillations are at 200 Hz.
The quality of this analysis also depends of the number of segments. The higher the number cause the more precise the segment determination. Moreover, the omissions that were made at the beginning should not be ignored and can have influence on the marked points.

12. Partial Correlation of Temperature, Frequency, and CT Measurement Error

If we calculate the connection between two phenomena, it is sometimes necessary to exclude the importance of a third variable that can affect the amount of connection between the first two variables. Partial correlation shows the correlation between two variables where the influence of the third variable is excluded. It is calculated according to the formula [4]:
r 12 / 3 = r 12 r 13 · r 23 1 r 13 2   ·   1 r 23 2
  • r12–correlation coefficient between the 1st and 2nd variables; r13–correlation coefficient between the 1st i 3rd variables; r23–correlation coefficient between the 2nd and 3rd variables
The correlation coefficient is defined by the formula:
r x y = ( X i X ¯ ) ( Y i Y ¯ ) i n ( X i X ¯ ) 2 + i n ( Y i Y ¯ ) 2
In the analysis of the influence of temperature and frequency on the accuracy of the current transformer, we will use the measurement results. Ten measurements were performed for four frequencies and four temperatures, for a total of 4 × 4 × 10 = 160 measurements. In Section 8, the functional dependence of the CT error on temperature was shown, and in Section 10, the functional dependence of the CT error on frequency was shown.
A prerequisite for the calculation of the partial correlation is the calculation of the matrix of all correlations, which is given in Table 5.
  • r12–correlation coefficient between temperature and frequency, r12 = 0
  • r13–correlation coefficient between temperature and CT error, r13 = −0.3137
  • r23–correlation coefficient between frequency and error CT-a, r23 = −0.0137
When calculating partial correlation, we obtain the “correct” correlation between the quantities, which, in the specific case, means:
  • The influence of frequency on the relationship between temperature and error is r13/2 = −0.3138 (Table 6), which is very close to the correlation of these two quantities, r13 = −0.3137 (Table 5). The conclusion is that the frequency does not influence the relationship between CT error and measurement temperature.
  • The influence of temperature on the frequency and CT error is r23/1 = −0.0144 (Table 6), which is very close to the linear correlation of these two quantities, r13 = −0.0137 (Table 5). The conclusion is that temperature has no significant influence on the relationship between CT error and measurement frequency.
It should be noted that the partial correlation r12/3 was not considered due to the illogicality of the results of that measurement.
  • r12/3–partial correlation coefficient between temperature and frequency, r12/3 = 0
  • r13/2–partial correlation coefficient between temperature and CT error, r13/2 = −0.3137
  • r23/1–partial correlation coefficient between frequency and CT error, r23/1 = −0.0137
At the end, Table 7 displays the changes after the correction of correlation and eliminates influences of the frequency (the change: −0.003%) and temperature (the change: −5.1%).

13. Conclusions

As is pointed out in Section 11, the Pearson correlation can be used to improve the model of CT. Using the proposed model, we can define the point where the deviation from the model is significant. There are two ways to be near to the ideal characteristic: To correct the mathematical model or to change the performance of CT. The Pearson correlation in this work shows possible errors which are recognized as changes of slope in comparison with the ideal characteristic.
Generally, the deviation of the measured and calculated value shown in Figure 13b and Figure 14a–d was the impetus for determining the partial correlation of temperature and frequency to the measurement error. The regression factor, same as Pearson correlation, has high values. The calculation proved that there is no significant influence of each of these quantities on the interconnection of another quantity with measurement error, except in the cases pointed out in Section 11.
Moreover, the influence of the digital measuring instrument is the limiting factor. The measurement limit of three decimal places can disturb the total result and bring in a disturbance.
At the end of experiment, the calculation of partial correlation between three factors is performed: temperature, frequency, and function of error. This verification is useful to estimate the influence of external variables such as temperature and frequency, to be sure they have an influence on the mutual interference of the two monitored variables.
The whole experiment can be considered successful from the perspective of developing a method of determining the parameters of the current transformer using a statistical method [10,11].

Author Contributions

Conceptualization D.B.; methodology, D.B.; audit parameters, D.B. and D.G.; validation, D.B. and S.N.; formal analysis, D.B. and N.M.; investigation, D.B.; resources, D.G. and N.M.; writing—original draft preparation, D.B.; writing—review and editing, S.N. and N.M.; supervision, D.B., S.N. and N.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. CT-Analyzer-Test Ceritficate nr. 111/12/C; Elektrosklop Ltd.: Bestovje, Croatia, 2018.
  2. Perić, N.; Petrović, I. Identifikacija Procesa—Predavanja; Fakultet Elektrotehnike i Računarstva—Zagreb: Zagreb, Croatia, 2002. [Google Scholar]
  3. Galić, R. Vjerojatnost i Statistika; Elektrotehnički Fakultet Sveučilišta J.J. Strossmayera u Osijeku: Osijek, Croatia, 1999. [Google Scholar]
  4. Luchka, A.Y. The Method of Averaging Functional Correction; Academic Press: New York, NY, USA; London, UK, 1963. [Google Scholar]
  5. Saleh, S.M.; El-Hoshy, S.H.; Gouda, O.E. Proposed diagnostic methodology using the cross-correlation coefficient factor technique for power transformer fault identification. In IET Electric Power Applications; 2017; pp. 412–422. Available online: https://ietresearch.onlinelibrary.wiley.com/doi/pdfdirect/10.1049/iet-epa.2016.0545 (accessed on 17 January 2023).
  6. Behjat, V.; Mahvi, M.; Rahimpour, E. New statistical approach to interpret power transformer frequency response analysis: Non-parametric statistical methods. In IET Science, Measurement & Technology; 2016; pp. 364–369. Available online: https://ietresearch.onlinelibrary.wiley.com/doi/epdf/10.1049/iet-smt.2015.0204 (accessed on 17 January 2023).
  7. Herceg, D. Numeričke i Statističke Metode u Obradi Eksperimentalnih Podataka; Institut za Matematiku: Novi Sad, Serbia, 1989. [Google Scholar]
  8. Yachikov, I.M.; Nikolaev, A.A.; Zhuravlev, P.Y.; Karandaeva, O.I. Estimate of Correlation between Transformer Diagnostic Variables with a Time Lag. Int. Conf. Ind. Eng. ICIE 2017, 206, 1794–1800. [Google Scholar] [CrossRef]
  9. Wang, X.; Li, Q.; Yang, R.; Li, C.; Zhang, Y. Diagnosis of solid insulation deterioration for power transformers with dissolved gas analysis-based time series correlation. IET Sci. Meas. Technol. 2014, 9, 393–399. [Google Scholar] [CrossRef]
  10. Guillen, D.; Esponda, H.; Vasquez, E.; Idarraga-Ospina, G. Idarraga-Ospina. Algorithm for transformer differential protection based on wavelet correlation modes. In IET Generation, Transmission & Distribution; 2016; pp. 2871–2879. Available online: https://ietresearch.onlinelibrary.wiley.com/doi/epdf/10.1049/iet-gtd.2015.1147 (accessed on 17 January 2023).
  11. Yu, T.; Yang, J.; Lin, J.C.; Geng, C. Correlation Analysis of Transformer Parameters Based on Pair-Copula. IOP Conf. Ser. Earth Environ. Sci. 2020, 571, 012005. [Google Scholar] [CrossRef]
Figure 1. Ideal characteristic, measured point and regression line of CT at f = 100 Hz, T = 26.3 °C.
Figure 1. Ideal characteristic, measured point and regression line of CT at f = 100 Hz, T = 26.3 °C.
Sensors 23 02704 g001
Figure 2. Interpolation “spline” curve in statistical experiment.
Figure 2. Interpolation “spline” curve in statistical experiment.
Sensors 23 02704 g002
Figure 3. Setting of the experiment.
Figure 3. Setting of the experiment.
Sensors 23 02704 g003
Figure 4. Equivalent circuit of the current transformer.
Figure 4. Equivalent circuit of the current transformer.
Sensors 23 02704 g004
Figure 5. Dependance of the inductance Lm on the secondary current of CT.
Figure 5. Dependance of the inductance Lm on the secondary current of CT.
Sensors 23 02704 g005
Figure 6. Test result of CT magnetization.
Figure 6. Test result of CT magnetization.
Sensors 23 02704 g006
Figure 7. Reduced model of current transformer.
Figure 7. Reduced model of current transformer.
Sensors 23 02704 g007
Figure 8. CT error function diagram.
Figure 8. CT error function diagram.
Sensors 23 02704 g008
Figure 9. Figure of the error function with the correction of the number of windings nk = 99.
Figure 9. Figure of the error function with the correction of the number of windings nk = 99.
Sensors 23 02704 g009
Figure 10. Factory testing results of ratio error.
Figure 10. Factory testing results of ratio error.
Sensors 23 02704 g010
Figure 11. Display of the ammeter calibration curve at a temperature of 26.3 °C.
Figure 11. Display of the ammeter calibration curve at a temperature of 26.3 °C.
Sensors 23 02704 g011
Figure 12. CT error depending on the secondary current for a temperature of 26.3 °C, f = 50 Hz. (a) measured (red); function (blue); difference (gray). (b) “Spline” curve of function error ef and measured error em; R2 = 0.9624.
Figure 12. CT error depending on the secondary current for a temperature of 26.3 °C, f = 50 Hz. (a) measured (red); function (blue); difference (gray). (b) “Spline” curve of function error ef and measured error em; R2 = 0.9624.
Sensors 23 02704 g012
Figure 13. CT error curves depending on the secondary current for the frequency: (a) f = 50 Hz measured (red); function (blue); difference (gray). (b) f = 100 Hz measured (red); function (blue); difference (gray). (c) f = 150 Hz measured (red); function (blue); difference (gray). (d) f = 200 Hz measured (red); function (blue); difference (gray).
Figure 13. CT error curves depending on the secondary current for the frequency: (a) f = 50 Hz measured (red); function (blue); difference (gray). (b) f = 100 Hz measured (red); function (blue); difference (gray). (c) f = 150 Hz measured (red); function (blue); difference (gray). (d) f = 200 Hz measured (red); function (blue); difference (gray).
Sensors 23 02704 g013
Figure 14. Curve: 1. f = 50 Hz; R2 = 0.9657, 2. f = 100 Hz; R2 = 0.9607, 3. f = 150 Hz; R2 = 0.9727, 4. f = 200 Hz; R2 = 0.916.
Figure 14. Curve: 1. f = 50 Hz; R2 = 0.9657, 2. f = 100 Hz; R2 = 0.9607, 3. f = 150 Hz; R2 = 0.9727, 4. f = 200 Hz; R2 = 0.916.
Sensors 23 02704 g014
Table 1. Example of functional and measured error for f = 100 Hz.
Table 1. Example of functional and measured error for f = 100 Hz.
No.1.2.3.
I (%)em (%)ef (%)
1.10%−1.000−1.574
2.20%0.000−0.620
3.30%0.3330.157
4.40%0.5000.165
5.50%0.6000.332
6.60%0.6670.432
7.70%0.7140.499
8.80%0.7500.547
9.90%0.6670.694
10.100%0.7000.713
st. dev. s sx = 0.5126sy = 0.6761
Table 2. Results of factory test of current transformer.
Table 2. Results of factory test of current transformer.
Nr.Test SettingsNr.Excitation TestNr.Ratio Test
1I-pn100 A14V-kn27.663 V27Ratio100:0.99766
2I-sn1 A15V-kn 2n.a.28e−0.2341%
3Rated burden 2.5 VA/116Ls0.0005523 H29ec1.2535%
4Operating burden2.5 VA/117Kr89.07%30Df41.82 min
5Applied standardIEC 61869-218I-kn0.08860231Polarityok
6Core type (P/M)P19I-kn 2n.a.32N99
7Class10 P20Lm1.2104 H
8ALF10.0
9f50.0 Hznr.Result with rated burden
10Tsn.a.21ALF11.63
11max. Rct0.531 W22Ts0.399 s
23ALFi11.4
nr.Resistance test
12Rmeas (25 °C)0.44532 Wnr.Result with operating burden
13Rref (75 °C)0.53112 W24ALF11.63
25Ts0.399 s
26ALFi11.4
Table 3. Factory data of the saturation magnetic curve.
Table 3. Factory data of the saturation magnetic curve.
Ie (A)0.0010.0020.0030.0040.0050.0060.0070.0080.0090.01
Es (V)0.0110.040.080.150.230.330.450.60.750.95
Lm = Es/(Ie·314) (H)0.0350.0640.0850.1190.1460.1750.2050.2390.2650.303
|Rct + jwLs| (Ω)0.500.500.500.500.500.500.500.500.500.50
Ict = Es/|Rct + jwLs| (A)0.020.080.160.300.460.660.901.201.501.90
Table 4. Correlation of ideal curve segments and measured curve segments.
Table 4. Correlation of ideal curve segments and measured curve segments.
Point1.
Segment
2.
I (%)
3.
Ideal
Curve SlopeSign
4.
Curve Slope Sign at 50 Hz/
Correlation 3. and 4.
5.
Curve Slope Sign at 100 Hz/
Correlation 3. and 5.
6.
Curve
Slope Sign at 150 Hz/
Correlation 3. and 6.
7.
Curve
Slope Sign at 200 Hz/
Correlation 3. and 7.
1. 10%
2.1–220%++/1+/1+/1+/1
3.2–330%++/1+/1+/1+/1
4.3–440%++/1+/1+/1+/1
5.4–550%+−/−1+/1+/1+/1
6.5–660%++/1+/1−/−1−/−1
7.6–770%++/1+/1+/1−/−1
8.7–880%++/1+/1+/1+/1
9.8–990%+−/−1−/−1+/1−/−1
10.9–10100%++/1+/1+/1−/−1
PEARSON CORRELATION “em ∧ ef0.98270.980.98620.957
Table 5. Correlation matrix.
Table 5. Correlation matrix.
Temp.; R1Freq.; R2Error CT; R3
Temp.; R11.00000.0000−0.3137
Freq.; R20.00001.0000−0.0137
Error CT; R3−0.3137−0.01371.0000
Table 6. Matrix of partial correlation–influence on connection.
Table 6. Matrix of partial correlation–influence on connection.
Freq.; R2CT Error; R3Temp.; R1
Temp/CT Error; R13−0.3138
Temp/Freq; R12 −0.0045
Freq./CT Error; R23 −0.0144
Table 7. Relations of values.
Table 7. Relations of values.
CorrelationValue Before
Correction
Value After
Correction
D%
1. Temp./CT Error−0.3137−0.3138 −0.003%
2. Freq./CT Error−0.0137 −0.0144−5.1%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Burgund, D.; Nikolovski, S.; Galić, D.; Maravić, N. Pearson Correlation in Determination of Quality of Current Transformers. Sensors 2023, 23, 2704. https://doi.org/10.3390/s23052704

AMA Style

Burgund D, Nikolovski S, Galić D, Maravić N. Pearson Correlation in Determination of Quality of Current Transformers. Sensors. 2023; 23(5):2704. https://doi.org/10.3390/s23052704

Chicago/Turabian Style

Burgund, Davorin, Srete Nikolovski, Dario Galić, and Nedeljko Maravić. 2023. "Pearson Correlation in Determination of Quality of Current Transformers" Sensors 23, no. 5: 2704. https://doi.org/10.3390/s23052704

APA Style

Burgund, D., Nikolovski, S., Galić, D., & Maravić, N. (2023). Pearson Correlation in Determination of Quality of Current Transformers. Sensors, 23(5), 2704. https://doi.org/10.3390/s23052704

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop