Next Article in Journal
Nanoarchitectonics of Layered Metal Chalcogenides-Based Ternary Electrocatalyst for Water Splitting
Next Article in Special Issue
Energy Efficiency in Electromagnetic and Electro-Permanent Lifting Systems
Previous Article in Journal
Detection of Blockages of the Belt Conveyor Transfer Point Using an RGB Camera and CNN Autoencoder
Previous Article in Special Issue
The Use of Hydrogen for Traction in Freight Transport: Estimating the Reduction in Fuel Consumption and Emissions in a Regional Context
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigating Empirical Mode Decomposition in the Parameter Estimation of a Three-Section Winding Model †

by
Daniel Marc Banks
,
Johannes Cornelius Bekker
* and
Hendrik Johannes Vermeulen
Department of Electrical and Electronic Engineering, Stellenbosch University, Stellenbosch 7600, South Africa
*
Author to whom correspondence should be addressed.
This paper is an extended version of our paper published in 2022 IEEE International Conference on Environment and Electrical Engineering and 2022 IEEE Industrial and Commercial Power Systems Europe (EEEIC/I&CPS Europe).
Energies 2023, 16(4), 1668; https://doi.org/10.3390/en16041668
Submission received: 30 December 2022 / Revised: 23 January 2023 / Accepted: 3 February 2023 / Published: 7 February 2023

Abstract

:
Parameter estimation represents an important aspect of modeling electromagnetic systems, and a wide range of parameter estimation strategies has been explored in literature. Most parameter estimation methodologies make use of either time-domain or frequency-domain responses as measured from the terminals of the device under test. Very limited research has, however, been conducted into exploring the use of modal decomposition strategies on the time-domain waveforms for parameter estimation applications. In this paper, the use of Empirical Mode Decomposition for estimating the parameters of a three-section lumped parameter transformer model is explored. A novel approach is proposed to define the optimization cost function in terms of the intrinsic modes of simulated time-domain waveforms. The results are compared with results obtained using classical time-domain and frequency-domain approaches. It is shown through an impulse response test that weighting the modes obtained from the Inferred Empirical Mode Decomposition approach presented in this work offers advantages in terms of accurately representing the target model transfer function dynamics and can assist in interpreting the various dynamic modes associated with the target model.

1. Introduction

The harmonics introduced by inverter-based renewable generation mean that transformers are being operated at voltages with a high level of high-frequency harmonic content [1]. This creates a need for accurate wideband equivalent-circuit transformer models that represent the expected transfer function characteristics under these operating conditions, especially for simulation studies aimed at characterizing system behavior. Additionally, being able to model transformer responses during transient conditions allows for early detection of failure and electrical faults using techniques such as Frequency Response Analysis (FRA) [2]. Aguglia [3] highlights that a model with good parameter accuracy is important to a transformer designer, while for control purposes, an accurate mathematical model relating the input to the output of a system is sufficient.
High frequency modeling of transformer windings ideally requires distributed winding models. These models are, however, commonly approximated using lumped-parameter equivalent circuit models, and the topologies of these lumped parameter models have been studied extensively in literature [4,5]. Various methodologies have been proposed for estimating the parameters of wideband transformer models. Keyhani [4,6] performed maximum likelihood estimation on a six-section transformer winding model and concluded that unique solutions can only be estimated if the initial values are close to the target values, or if sectional voltage measurements were included. Brozio [5] used a constrained Sequential Quadratic Programming (SQP) optimization algorithm to estimate the parameters of a two-winding power transformer from measured frequency responses. Time-domain signals were not employed due to computational limitations, and the overall simplicity of representing the frequency responses analytically. In a more recent contribution, Chanane [7,8,9] applied metaheuristic optimization algorithms, such as Particle Swarm Optimization (PSO) [10], Crow Search Optimization (CSO) [11], and Grey Wolf Algorithm (GWA) [12], to estimate the parameters of a fully interleaved continuous-disc winding. The parameter estimation process was simplified by making use of measured values for the DC resistance, equivalent inductance, and effective capacitance of the transformer winding.
Empirical Mode Decomposition (EMD), proposed by Huang, represents a methodology for decomposing a multi-modal signal into simpler components known as Intrinsic Mode Functions (IMFs), which can be used to extract features from nonlinear and non-stationary waveforms [13,14]. Alternative noise-assisted EMD methods, such as Ensemble Empirical Mode Decomposition (EEMD) and signal masking techniques, have been subsequently proposed to address the issue of mode-mixing [13,14,15]. EMD and variations thereof have been used extensively in literature for transformer vibration analysis [16,17,18,19]. No existing research applies EMD in the parameter estimation of transformer winding or complete transformer models.
This paper investigates the use of EMD for estimating the parameters of a three-section lumped-parameter transformer winding model. The main research objective of the work is to explore how the classical parameter estimation methodologies can be amended to incorporate EMD, and whether EMD can improve the performance of these approaches. As a second objective, the work investigates the performance of the Pseudo-Random Impulse Sequence (PRIS) [20,21,22] for the parameter estimation of transformer models.
The model is perturbed using a PRIS, and the simulated time-domain voltage waveforms of the winding are decomposed into IMFs. A novel approach is proposed to derive the cost function from these IMFs, the approach is further expanded to investigate weighting of the derived intrinsic modes. The accuracy of the estimated model parameters is interpreted by comparing the impulse responses of the targeted and estimated models in addition to consideration of classical error metrics.
This paper is structured as follows. In Section 2, the transformer winding model used for the study is presented together with its analytical input impedance frequency response. Section 3 presents the PRIS source used to perturb the transformer winding. Section 4 discusses the parameter estimation methodologies being applied. The different cost function formulations used for the parameter estimation algorithm are presented in Section 5. Section 6 presents the results obtained from the different parameter estimation approaches. Additional analysis and results are also presented to provide insight into the results. The paper is concluded in Section 7 with recommendations for further research.

2. Transformer Winding Model

The investigation targets the model proposed by Chanane [8] for a fully interleaved continuous-disc winding with 5 discs, consisting of 30 turns per disc, of copper conductor with a cross-sectional area of 10 mm2. The winding features an insulated core with insulation thickness and duct spacing corresponding to a 30 kV insulation rating. Figure 1 shows the model topology, while Table 1 summarizes the target model parameters used by Chanane and in this investigation [8]. The model has a total of 16 parameters to be estimated. These can, however, be reduced by assuming that the winding sections have the same physical dimensions such that R 1 = R 2 = R 3 = R s and L 1 = L 2 = L 3 = L s . The dimensions between sections are also assumed to be similar, and thus C 12 = C 23 = C 34 = C s . It is furthermore assumed that the capacitances from the individual sections to ground are equal, aside from the sections closest to the grounded core and tank, which is assumed to be half the ground capacitance, so C g 1 = C g 4 = C g 2 and C g 2 = C g 3 = C g . These assumptions are used extensively in literature and have shown to accurately model the wideband responses of the transformer winding [7,8,9,23,24,25].

Analytical Input Impedance Frequency Response

The input impedance Z i ( s ) of the circuit model shown in Figure 1 can be derived from the impedance matrix Z of the model, using the relationship
Z i ( s ) = V i n ( s ) I i n ( s ) = Δ Z C Z ( 7 , 7 ) ,
where Δ Z denotes the determinant of Z and C Z ( 7 , 7 ) denotes the (7,7)-th cofactor of Z . The resulting analytical expression for the transfer function Z i ( s ) can be expressed in the form
Z i ( s ) = α 1 s 5 + α 2 s 4 + α 3 s 3 + α 4 s 2 + α 5 s + α 6 β 1 s 6 + β 2 s 5 + β 3 s 4 + β 4 s 3 + β 5 s 2 + β 6 s + β 7 ,
where α and β denote constant coefficients defined in terms of the circuit parameters.

3. Target Model Perturbation

The winding model is excited at the terminals in MATLAB Simulink using the PRIS perturbation arrangement presented in Figure 2. The operation of the PRIS source is extensively discussed in literature [20,21,22], it is shown that the spectral content and SNR of the PRIS source is controllable. The controllability makes it well suited for exciting electromagnetic equipment such as transformer systems where direct current or low-frequency components of the perturbation signal are undesirable as these can lead to core saturation [20,21,22]. The RLC values of the source, i.e., R p r i s = 150 Ω , L p r i s = 4 mH and C p r i s = 1   μ F are obtained by following the source design considerations presented in [20], and V D C = 100 V.
The PRIS is generated using a Pseudo-Random Binary Sequence (PRBS) clock frequency, f c l k , of 100 kHz, and a 12th order topology for frequency response measurements and a 6th order topology for time-domain measurements. The higher order PRBS for the frequency-domain measurements is required to accommodate a high frequency resolution over a wide frequency band. Figure 3 shows typical simulated time-domain waveforms of the input voltage and input current, together with the PRBS gating signals. The simulation is run for three 6th order topology cycles. In order to remove the effects of the initial transients in the simulation, only the last sequence is captured for analysis.

4. Parameter Estimation Methodology

The parameter estimation approaches target the parameter vector given in (3), where the parameters are normalized by applying the normalization constants given in (4). The optimization bounds for parameters other than the coupling coefficients are chosen to be an entire order of magnitude so as to assume little a priori information about the parameter values. The coupling coefficient bounds are selected close to unity to reflect the assumption of a well designed transformer [26]. Table 2 summarizes the bounds of the search space adopted in the investigation.
θ = [ R s , C s , C g , L s , k 12 , k 13 , k 23 ]
θ n o r m = [ 10 3 , 10 12 , 10 9 , 10 6 , 10 2 , 10 2 , 10 2 ]
Figure 4 and Figure 5 present an overview of the parameter estimation strategies implemented in the investigation. In simulation, the target model is perturbed using a PRIS source [20]. In practice, it cannot be guaranteed that the time-domain samples obtained from the device under test align perfectly with the samples obtained through simulation of the estimated model. Two strategies to align these time-domain waveforms are considered. The first strategy, as presented in Figure 4a and Figure 5a, uses the sampled input current, i i n ( t ) , produced during the target model perturbation as an input to the estimated model. The second strategy, as presented in Figure 4b and Figure 5b, simulates the target and estimated model PRIS perturbation arrangements separately, and then aligns the waveforms by aligning the bipolar PRBS voltages produced by the H-Bridge, v p r b s ( t ) , in each simulation.
The initial studies considered various optimisation algorithms, including interior-point, patternsearch, particle swarm and genetic algorithm. The particle swarm optimization algorithm, which is a global solver that generates a population of points to find a global minimum in a widely bounded search space, is used throughout the investigation [8]. The algorithm is chosen based on the assumption that limited a priori knowledge is available on the model parameters to be estimated and the superior runtimes in comparison to the genetic algorithm. The algorithm runs 250 search iterations, with a population of 1000 points generated using a random seed. After the global optimization procedure is completed, the optimization results are refined using a local solver, fmincon. The classical time-domain approach to parameter estimation, whereby the measured output voltage waveforms from the target and estimated models are compared and used in a cost function formulation, is presented in Figure 4.
Figure 5 presents the proposed EMD approach, whereby the time-domain voltage waveforms are first decomposed into IMFs, from which the cost function is derived. Cost functions defined in terms of both Root Mean Squared Error (RMSE) and correlation coefficient ( ρ ) based metrics are implemented for all approaches.
Three implementations of modal decomposition are investigated. The first performs standard EMD on the target model and estimated model voltage waveforms, and formulates the cost function in terms of the resulting IMFs and residuals. The second computes inferred IMFs to represent the estimated model, from which the cost function is derived. The third derives the cost function from weighted inferred IMFs to represent the estimated model.

5. Parameter Estimation Cost Function Formulations

5.1. Frequency-Domain Approach

To validate the Simulink model, the frequency response of the input impedance of the target model is obtained from the simulated input voltage and input current signals during PRIS perturbation using Welch’s method [27], and compared to the analytical frequency response. Figure 6 shows the frequency responses of the input impedance of the target model obtained analytically and through perturbation. The responses exhibit a strong correlation between the analytical and estimated frequency responses.
For the frequency-domain parameter estimation, the RMSE based cost function, C r m s e z , is calculated as
C r m s e z = 1 N k k = 1 N k ( ϵ ( ω k ) ) 2 ,
where
ϵ ( ω k ) = log 10 Z i ( ω k ) log 10 Z i ˜ ( ω k ) , k = 1 , . . . , N k .
The error ϵ ( ω k ) between the input impedance frequency responses of the target model Z i ( ω ) and estimated model Z i ˜ ( ω ) is calculated on a point-by-point basis for each frequency sample up to N k samples. The sample frequencies ω k are distributed logarithmically such that the contributions from ϵ ( ω k ) to the cost function are distributed equitably across the frequency range of interest.
The correlation coefficient metric is used to compute the correlation based cost function, C ρ z , as
C ρ z = ρ { Z i ( ω ) , Z ˜ i ( ω ) } + 1
by adjusting the correlation coefficient, ρ , between Z i ( ω ) and Z ˜ i ( ω ) such that a value of zero represents a perfect match in frequency responses.

5.2. Time-Domain Approach

The time-domain target voltage waveform, shown in Figure 3, is defined by v i n ( t ) . The cost function using the RMSE metric, C r m s e v , is calculated as
C r m s e v = 1 N k k = 1 N k ( ϵ ( t k ) ) 2 ,
where
ϵ ( t k ) = v i n ( t k ) v ˜ i n ( t k ) , k = 1 , 2 , . . . , N k .
The errors between the voltage waveforms of the target model v i n ( t ) , and estimated model v ˜ i n ( t ) are calculated on a point-by-point basis for each time sample t k up to N k samples.
Similar to the frequency-domain approach in (7), the correlation coefficient based cost function, C ρ v , is defined in terms of the correlation between the target and estimated voltage waveforms as
C ρ v = ρ { v i n ( t ) , v ˜ i n ( t ) } + 1 .

5.3. Empirical Mode Decomposition Approach

The simulated time-domain voltage waveform shown in Figure 3 is decomposed into IMFs using EMD [13]. This yields six IMFs and one residual, as shown in Figure 7.
The voltage waveforms associated with the target model and estimated model, denoted by v ( t ) and v ˜ ( t ) , respectively, can be represented by the relationships
v ( t ) = i = 1 N i v i m ( t ) + v r ( t )
and
v ˜ ( t ) = i = 1 N ˜ i v ˜ i m ( t ) + v ˜ r ( t ) ,
where v i m ( t ) and v ˜ i m ( t ) denote the i t h IMFs, N i and N ˜ i denotes the total number of IMFs, and v r ( t ) and v ˜ r ( t ) denote the residuals, of the target model and estimated model, respectively. The IMFs and residuals for the target model and estimated model are represented by the sets Γ and Γ ˜ , respectively. The set sizes of the target model and estimated model are denoted as N r and N ˜ r , respectively, as shown in (13) and (14).
Γ = { γ r ( t ) r = 1 , 2 . . . N r }
Γ ˜ = { γ r ˜ ( t ) r = 1 , 2 . . . N ˜ r }
The set with fewer IMFs is padded with zeros such that the sets are of equal length and can be compared.
The Normalized RMSE (NRMSE) based cost function, C n r m s e m , for the Empirical Mode Decomposition approach is calculated as
C n r m s e m = r = 1 N r 1 ψ r 1 N k k = 1 N k ( ϵ r ( t k ) ) 2 ,
where
ϵ r ( t k ) = γ r ( t k ) γ ˜ r ( t k ) .
The error ϵ r ( t k ) is thus defined as the difference between the r t h waveforms in sets Γ and Γ ˜ at sample time t k .
The respective errors are then normalized by applying the normalization constant
ψ r = max { γ r ( t ) } min { γ r ( t ) } ,
where max { γ r ( t ) } and min { γ r ( t ) } denote the maximum and minimum amplitudes of the r t h target IMF waveform, respectively. This normalization is required to counter biasing of the cost function towards waveforms with higher amplitude, as the measurement set consists of waveforms with varying amplitudes. The normalization constant is omitted in cases when the target data set is zero-padded as it results in a ψ r = 0 .
In the case of the correlation coefficient based cost function, the larger set between the target and estimated sets is truncated to the length of the shorter set, denoted by N s , as correlation coefficient metrics do not allow for zero-padded waveforms. Due to this manipulation, a variance occurs in the number of IMFs from iteration to iteration in the optimization procedure, and the result is thus normalized by N s . The correlation coefficient based cost function, C ρ m , is thus defined as
C ρ m = 1 N s r = 1 N s [ ρ { γ r ( t ) , γ ˜ r ( t ) } + 1 ] .

5.4. Inferred Empirical Mode Decomposition

In this method, the IMFs produced by the estimated waveforms are never computed. Instead, a set of inferred IMFs (iIMFs) is produced to represent the estimated model. The target model set, Γ , remains the same as in (13). The r t h iIMF, v ˜ r i ( t ) , for the estimated model voltage waveform, v ˜ ( t ) , and is determined by subtracting waveforms 1 to r 1 and waveforms r + 1 to N r in Γ from v ˜ ( t ) . This is represented by the relationship
v ˜ r i ( t ) = v ˜ ( t ) k = 1 r 1 γ k ( t ) k = r + 1 N r γ k ( t ) .
The estimated model set, while in the same form as in (14), consists of the iIMFs of the estimated voltage waveform as given in (19). The number of iIMFs produced is always the same as the number of target IMFs, and the length of Γ ˜ thus becomes N r , such that
Γ ˜ = { γ r ˜ ( t ) r = 1 , 2 . . . N r } .
This decomposition approach creates a set of iIMFs that each contain the error dynamics of the estimated signal, whilst also containing the corresponding target IMF dynamics that are present in the estimated waveform. The cost function for the NRMSE based metric is calculated as
C n r m s e i = r = 1 N r 1 ψ r 1 N k k = 1 N k ( ϵ r ( t k ) ) 2 ,
where the errors, ϵ r ( t k ) , and normalization constants, ψ r , are defined as in (16) and (17), respectively.
The correlation coefficient based cost function in the context of the Inferred EMD (IEMD) approach does not require truncation of the IMFs nor normalization. The cost function C ρ i is thus reduced to the form
C ρ i = r = 1 N r [ ρ { γ r ( t ) , γ ˜ r ( t ) } + 1 ] .

5.5. Weighted Inferred Empirical Mode Decomposition

In this approach, the set of correlation coefficients between the target model IMF waveform and the respective iIMF waveform of the estimated model are weighted using a binary weighting vector, v ( r ) . The weighted correlation coefficient cost function is defined as
C ρ w i = r = 1 N r v ( r ) [ ρ { γ r ( t ) , γ ˜ r ( t ) } + 1 ] ,
where v ( r ) is set to either 0 or 1 for the r t h waveform in both the target model and estimated model sets. This allows various combinations of correlation coefficients between the IMFs and iIMFs to be included in the cost function formulation. The residual IMF waveform in the set of target model IMFs predominantly contains low-order dynamics. It is expected that the model parameters to not have a significant effect on the low-frequency dynamics, as the resonant points are located above 10 kHz. The residual IMF is therefore excluded in the weighted cost function formulation.

6. Results

Table 3 presents the results of the frequency-domain parameter estimation investigations, while Table 4 presents the results of the parameter estimation investigations that utilise the time-domain, EMD and IEMD waveforms. The results are presented as parameter error percentages, where the error (err%) of the respective parameters is calculated using the relationship
err % ( n ) = θ ( n ) θ ˜ ( n ) θ ( n ) × 100 [ % ] ,
where θ ˜ denotes the final estimated parameter value and θ denotes the target value of the parameter, as presented in Table 1.
The frequency-domain RMSE based optimization yields errors below 10% for all parameters aside from R s . Parameter R s does not influence the cost function significantly, as it predominantly affects the damping at the resonant points of the frequency response. In considering Table 3 and Table 4, it is observed that the run times of frequency-domain parameter estimation, in most cases, are longer compared to time-domain and modal decomposition based approaches. This is attributed to the need for longer simulation times to obtain a high frequency resolution over a wide frequency band.
It is clear from Table 4 that cost functions using the EMD approach are not able to obtain parameters as accurately as the classical time-domain approach. This is attributed to the content within each IMF being dependent on the dynamics present in the original signal. EMD does not guarantee the presence of specific modes within a specific IMF. Although the presence of certain modes within the estimated waveform may be correct, it cannot be predicted which IMF these will be decomposed into. This characteristic makes the optimization challenging when only comparing the n t h IMF of the target with the n t h IMF of the estimated waveforms. The EMD approach is thus very sensitive to minor changes in the time-domain waveforms, thereby making the search space extremely stochastic, as shown in Figure 8. Figure 8, Figure 9, Figure 10 and Figure 11 present graphical representations of some of the search spaces. The figures are created by incrementally moving the parameter away from its target value by 0.2%. The search spaces are one-dimensional, meaning that each parameter is varied individually whilst all others are kept at their target values. For legibility purposes Figure 8, Figure 9, Figure 10 and Figure 11 show the cost function as R s , C s , C g and L s are varied. Figure 8 shows the search space of C n r m s e m generated using the first alignment strategy, but this stochasticity is typical of all of the search spaces generated using the EMD approach.
The RMSE based cost functions using the first alignment strategy for the classical time-domain and IEMD methodologies converge to the same parameter values. This is attributed to the search spaces being very similar with slight differences in amplitude. The IEMD approach does, however, provide a noticeable improvement in runtime, as well as lower error percentages when comparing the results of the second time-domain alignment strategy. The IEMD search space of the NRMSE based cost function, C n r m s e i , using the first alignment strategy is shown in Figure 9. The lower parameter errors obtained through the time-domain and IEMD approaches, in comparison with the EMD approach, can be attributed to the overall downward trend in the cost function as the parameter values range closer to the target model parameters in Figure 9, as opposed to the stochastic search space shown in Figure 8.
In the case of correlation coefficient based cost functions, the time-domain and IEMD approaches converge to different results, where the IEMD approach obtains parameter values with lower error percentages in alignment strategy 1. This can be attributed to the correlation based cost function formulations C ρ v and C ρ i producing different search spaces, where the IEMD search space of alignment strategy 1, shown in Figure 10, places equal emphasis on all modes of the system. The weighted IEMD approach is introduced to investigate whether using different combinations of target IMFs and estimated iIMFs can assist in improving the parameter estimation results. This investigation is limited to binary weights that either include or exclude certain IMFs and iIMFs from the cost function formulation. The parameter estimation procedure is executed for all possible weighting combinations for the vector v ( r ) . The best result is presented in Table 5.
The best results for alignment strategy 1 are obtained with v ( r ) = [ 1 , 1 , 0 , 0 , 0 , 0 , 0 ] , whilst the best results for alignment strategy 2 are obtained with v ( r ) = [ 0 , 1 , 0 , 1 , 0 , 1 , 0 ] . The weighted IEMD approach cycles through all possible weighting options and completes the parameter estimation procedure for each option. The runtimes are, therefore, significantly longer in comparison with the approaches presented in Table 3 and Table 4. In comparison with the error percentages achieved with the correlation coefficient based IEMD approach, there is improvement in the resistive, inductive, and coupling elements of the model using alignment strategy 1, by only including the higher order dynamics of the first and second IMF. Alignment strategy 2 yields no significant improvement in the individual parameter error percentages by comparing only the 2nd, 4th, and 6th IMFs and iIMFs. The search space in Figure 11 is generated using the weighted IEMD approach cost function formulation for time-domain alignment strategy 2.
Figure 11 is less stochastic compared to the search spaces generated using alignment strategy 1. The stochasticity of the previous search spaces is due to the discretized time-domain current waveform subsequently being used as an input to the estimated model simulation. Due to the discretization of the current waveform, the simulated input current to the estimated model contains less frequency content compared to the simulated current waveform of the target model. Strategy 2 allows the estimated model to simulate the input current over a wider range of frequencies.
As discussed by Keyhani [6] and Banks [22], determining unique transformer model parameters cannot be guaranteed when only terminal measurements are available. This is due to the fact that some of the individual parameters may not be fully observable from terminal measurements alone. While obtaining the correct model parameters is important for the condition monitoring of equipment, obtaining a model with a transfer function that is representative of the device under test is adequate for many simulation purposes [3]. Therefore, in addition to considering the error percentages, the performance of the estimated model is also interpreted by considering the impulse response. Impulse response tests are commonly applied to linear dynamic systems in control theory to determine their time-domain properties [29]. The ideal impulse exhibits a flat frequency response magnitude across the entire frequency spectrum, and the impulse response represents the transfer function of the system in the time domain. The impulse response waveforms of the target and estimated models are generated through the inverse Laplace transform of their analytical input impedance transfer functions. The waveforms produced by the target model and the estimated model are compared through an RMSE metric, and the results are presented in Table 6.
With respect to the results presented in Table 4, Table 5 and Table 6, it is clear that the lowest individual error percentages do not always guarantee the most accurate impulse response compared to the target model. For example, the RMSE-based cost function formulation of the time-domain and IEMD approaches using strategy 1 yields the lowest individual error percentages. However, in terms of the impulse response test, the estimated models obtained from the RMSE frequency-domain, correlation coefficient IEMD, and weighted IEMD approaches reproduce the dynamic modes of the target model more accurately. No significant improvement in parameter error percentages is achieved when applying the weighted IEMD approach compared to the IEMD approach when using alignment strategy 2. The impulse response RMSE metrics, however, show significant improvement. This is an indication that the model obtained through the weighted IEMD approach represents the transfer function dynamics of the target model more accurately. This can be due to combinations of various parameters affecting the transfer function, rather than the values of individual parameters. It is also possible that not all parameters influence the model’s response equitably, with some parameters influencing the overall response in the frequency range of interest more than others. It is clear from Table 6 that the weighted IEMD approach using time-domain alignment strategy 2 produces the lowest RMSE, highlighted in green. It can, therefore, be concluded that this approach produces transfer function dynamics most similar to the target model.
The best impulse response results for the different alignment strategies seem to differ, and it is therefore not possible to form a general assumption as to which alignment strategy is preferred.

7. Conclusions

This paper explores the use of EMD to estimate the parameters of a three-section transformer winding from time-domain waveforms obtained whilst applying a PRIS perturbation signal. A novel approach is proposed for deriving the cost function from the IMFs.
The results show that making use of standard EMD hinders the optimization procedure and that some pre-processing of the IMFs is required to arrive at a successful result. Using the target model IMFs and removing them from the estimated model voltage waveform to create iIMFs produced similar results to time-domain parameter estimation approaches. The results are cross-validated through calculating the impulse response of the target model and estimated model for each approach. It is shown that the inclusion of a weighting vector that includes only certain iIMFs in the correlation coefficient cost function improves the impulse responses, and therefore the overall transfer function accuracy, of the estimated model. It is shown that the most accurate set of model parameters does not provide the most accurate impulse response, and, therefore, should not be the only metric of the accuracy of a model when transfer function characteristics are the main consideration.
The research introduces a methodology for determining the optimal IMF weighting vector by cycling through all possible combinations. The methodology proves the application of modal decomposition in the parameter estimation of transformer winding models. Improved modal decomposition methodologies such as Ensemble Empirical Mode Decomposition (EEMD) and Iterative Filtering should be investigated as potential options for further improvements to the methodology. EEMD, however, introduces a significant additional computational burden. Determining the optimal weighting vector is, however, time consuming. There is therefore scope to investigate other methodologies for determining a set of weighting vectors that improves the parameter estimation results obtained from the EMD approach. In this work, the weighting vector is limited to binary values of either 0 or 1 in order to decrease complexity. Other weighting values should be investigated.

Author Contributions

Conceptualization, J.C.B. and H.J.V.; methodology, D.M.B.; software, D.M.B.; validation, D.M.B., H.J.V. and J.C.B.; formal analysis, D.M.B.; investigation, D.M.B.; resources, H.J.V.; data curation, D.M.B.; writing—original draft preparation, D.M.B.; writing—review and editing, H.J.V. and J.C.B.; visualization, D.M.B.; supervision, J.C.B. and H.J.V.; project administration, J.C.B.; funding acquisition, J.C.B. and D.M.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by an institutional bursary fund of Stellenbosch University, i.e., the DW-Ackermann-Trust.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available because it is confidential.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Peterson, B.; Rens, J.; Botha, M.G.; Desmet, J. On Harmonic Emission Assessment: A Discriminative Approach. Saiee Afr. Res. J. 2017, 108, 165–173. [Google Scholar] [CrossRef]
  2. Gomez-Luna, E.; Aponte, M.G.; Gonzalez-Garcia, C.; Pleite, G.J. Current Status and Future Trends in Frequency-Response Analysis With a Transformer in Service. IEEE Trans. Power Deliv. 2013, 28, 1024–1031. [Google Scholar] [CrossRef]
  3. Aguglia, D.; Viarouge, P.; de Almeida Martins, C. Frequency-Domain Maximum-Likelihood Estimation of High-Voltage Pulse Transformer Model Parameters. IEEE Trans. Ind. Appl. 2013, 49, 2552–2561. [Google Scholar] [CrossRef]
  4. Keyhani, A.; Chua, S.W.; Sebo, S.A. Maximum likelihood estimation of transformer high frequency parameters from test data. IEEE Trans. Power Deliv. 1991, 6, 858–865. [Google Scholar] [CrossRef]
  5. Brozio, C.C.; Vermeulen, H.J. Wideband equivalent circuit modelling and parameter estimation methodology for two-winding transformers. IEE Proc.-Gener. Transm. Distrib. 2003, 6, 487–492. [Google Scholar] [CrossRef]
  6. Keyhani, A.; Tsai, H.; Abur, A. Maximum likelihood estimation of high frequency machine and transformer winding parameters. IEEE Trans. Power Deliv. 1990, 5, 212–219. [Google Scholar] [CrossRef]
  7. Chanane, A.; Houassine, H.; Bouchhida, O. Enhanced modelling of the transformer winding high frequency parameters identification from measured frequency response analysis. IET Gener. Transm. Distrib. 2019, 13, 1339–1345. [Google Scholar] [CrossRef]
  8. Chanane, A.; Bouchhida, O.; Houassine, H. Investigation of the transformer winding high-frequency parameters identification using particle swarm optimisation method. IET Electr. Power Appl. 2016, 10, 923–931. [Google Scholar] [CrossRef]
  9. Chanane, A.; Belazzoug, M. On accuracy of a mutually coupled ladder network model high-frequency parameters identification for a transformer winding using gray wolf optimizer method. COMPEL-Int. J. Comput. Math. Electr. Electron. Eng. 2019, 40, 40–50. [Google Scholar] [CrossRef]
  10. Eberhart, R.; Kennedy, J. A new optimizer using particle swarm theory. In Proceedings of the Sixth International Symposium on Micro Machine and Human Science, Nagoya, Japan, 4–6 October 1995. [Google Scholar]
  11. Askarzadeh, A. A novel metaheuristic method for solving constrained engineering optimization problems: Crow search algorithm. Comp. Struct. 2016, 169, 1–12. [Google Scholar] [CrossRef]
  12. Mirjalili, S.; Mirjalili, S.M.; Lewis, A. Grey Wolf Optimizer. Adv. Eng. Softw. 2014, 69, 46–61. [Google Scholar] [CrossRef]
  13. Huang, N.E.; Shen, Z.; Long, S.R.; Wu, M.C.; Shih, H.H.; Zheng, Q.; Yen, N.-C.; Tung, C.C.; Liu, H.H. The Empirical Mode Decomposition and the Hilbert Spectrum for Nonlinear and Non-Stationary Time Series Analysis. Proc.-Math. Phys. Eng. Sci. 1998, 454, 903–995. [Google Scholar] [CrossRef]
  14. Stallone, A.; Cicone, A.; Materassi, M. New insights and best practices for the successful use of Empirical Mode Decomposition, Iterative Filtering and derived algorithms. Sci. Rep. 2020, 10, 15161. [Google Scholar] [CrossRef] [PubMed]
  15. Deering, R.; Kaiser, J.F. The use of a masking signal to improve empirical mode decomposition. In Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing, Philadelphia, PA, USA, 18–23 March 2005. [Google Scholar]
  16. Geng, C.; Wang, F.; Zhang, J.; Zhijian, J. Modal parameters identification of power transformer winding based on improved Empirical Mode Decomposition method. Electr. Power Syst. Res. 2014, 108, 331–339. [Google Scholar] [CrossRef]
  17. Yu, H.; Kong, L. Research on Modal Parameter Identification Method of Power Transformer Winding. In Proceedings of the 37th Chinese Control Conference (CCC), Wuhan, China, 25–27 July 2018. [Google Scholar]
  18. Xiong, W.; Li, J.; Pan, H. Analysis of transformer winding vibration based on modified empirical mode decomposition. In Proceedings of the 8th World Congress on Intelligent Control and Automation, Jinan, China, 7–9 July 2010. [Google Scholar]
  19. Xiong, W.; Pan, H. Transformer Winding Vibration Enveloping for Empirical Mode Decomposition Based on Non-uniform B-Spline Fitting. In Proceedings of the Second International Symposium on Knowledge Acquisition and Modeling, Wuhan, China, 30 November–1 December 2009. [Google Scholar]
  20. Mwaniki, F.M.; Vermeulen, H.J. Characterization and Application of a Pseudo-random Impulse Sequence for Parameter Estimation Applications. IEEE Trans. Instrum. Meas. 2020, 69, 3917–3927. [Google Scholar] [CrossRef]
  21. Gerber, I.P.; Mwaniki, F.M.; Vermeulen, H.J. Parameter Estimation of a Ferro-Resonance Damping Circuit using Pseudo-Random Impulse Sequence Perturbations. In Proceedings of the 56th International Universities Power Engineering Conference (UPEC), Middlesbrough, UK, 31 August–3 September 2021. [Google Scholar]
  22. Banks, D.M.; Bekker, J.C.; Vermeulen, H.J. Parameter Estimation of a Two-Section Transformer Winding Model using Pseudo-Random Impulse Sequence Perturbation. In Proceedings of the 56th International Universities Power Engineering Conference (UPEC), Middlesbrough, UK, 31 August–3 September 2021. [Google Scholar]
  23. Satish, L.; Sahoo, S.K. Locating faults in a transformer winding: An experimental study. Electr. Power Syst. Res. 2009, 79, 89–97. [Google Scholar] [CrossRef]
  24. Mukherjee, P.; Satish, L. Construction of equivalent circuit of a transformer winding from driving-point impedance function-analytical approach. IET Electr. Power Appl. 2012, 6, 172–180. [Google Scholar] [CrossRef]
  25. Ragavan, K.; Satish, L. Localization of Changes in a Model Winding Based on Terminal Measurements: Experimental Study. IEEE Trans. Power Deliv. 2007, 22, 1557–1565. [Google Scholar] [CrossRef]
  26. Wang, F.; Ren, R.; Yin, X.; Yu, N.; Yang, Y. A transformer with high coupling coefficient and small area based on TSV. Integration 2021, 81, 211–220. [Google Scholar] [CrossRef]
  27. MATLAB Signal Processing Toolbox. Available online: https://www.mathworks.com/products/dsp-system.html (accessed on 12 June 2022).
  28. Banks, D.M.; Bekker, J.C.; Vermeulen, H.J. Parameter Estimation of a Three-Section Winding Model using Intrinsic Mode Functions. In Proceedings of the IEEE International Conference on Environment and Electrical Engineering and IEEE Industrial and Commercial Power Systems Europe (EEEIC/I&CPS Europe), Prague, Czech Republic, 28 June–1 July 2022. [Google Scholar]
  29. Franklin, G.F.; Powell, D.J.; Emami-Naeini, A. Feedback Control of Dynamic Systems, 4th ed.; Prentice Hall PTR: Hoboken, NJ, USA, 2001. [Google Scholar]
Figure 1. Three-Section Lumped Parameter Winding Model [8].
Figure 1. Three-Section Lumped Parameter Winding Model [8].
Energies 16 01668 g001
Figure 2. PRIS Test Arrangement to Perturb the Three-Section Winding Model.
Figure 2. PRIS Test Arrangement to Perturb the Three-Section Winding Model.
Energies 16 01668 g002
Figure 3. Typical Simulated Input Voltage and Input Current Waveforms.
Figure 3. Typical Simulated Input Voltage and Input Current Waveforms.
Energies 16 01668 g003
Figure 4. Classical Time-domain Parameter Estimation Methodology: (a) Alignment Strategy 1 and (b) Alignment Strategy 2.
Figure 4. Classical Time-domain Parameter Estimation Methodology: (a) Alignment Strategy 1 and (b) Alignment Strategy 2.
Energies 16 01668 g004
Figure 5. Time-domain Modal Decomposition Parameter Estimation Methodology: (a) Alignment Strategy 1 and (b) Alignment Strategy 2.
Figure 5. Time-domain Modal Decomposition Parameter Estimation Methodology: (a) Alignment Strategy 1 and (b) Alignment Strategy 2.
Energies 16 01668 g005
Figure 6. Frequency Responses of the Input Impedance of the Target Model Obtained Analytically and Through Perturbation [28]. (© 2023 IEEE. Reprinted, with permission, from EEEIC/I&CPS Europe).
Figure 6. Frequency Responses of the Input Impedance of the Target Model Obtained Analytically and Through Perturbation [28]. (© 2023 IEEE. Reprinted, with permission, from EEEIC/I&CPS Europe).
Energies 16 01668 g006
Figure 7. Target Model Voltage Intrinsic Mode Functions [28]. (© 2023 IEEE. Reprinted, with permission, from EEEIC/I&CPS Europe).
Figure 7. Target Model Voltage Intrinsic Mode Functions [28]. (© 2023 IEEE. Reprinted, with permission, from EEEIC/I&CPS Europe).
Energies 16 01668 g007
Figure 8. Empirical Mode Decomposition Cost Function C n r m s e m versus Distance from the Target Model Parameter [28]. (© 2023 IEEE. Reprinted, with permission, from EEEIC/I&CPS Europe).
Figure 8. Empirical Mode Decomposition Cost Function C n r m s e m versus Distance from the Target Model Parameter [28]. (© 2023 IEEE. Reprinted, with permission, from EEEIC/I&CPS Europe).
Energies 16 01668 g008
Figure 9. Inferred EMD Cost Function C n r m s e i versus Distance from the Target Model Parameter [28]. (© 2023 IEEE. Reprinted, with permission, from EEEIC/I&CPS Europe).
Figure 9. Inferred EMD Cost Function C n r m s e i versus Distance from the Target Model Parameter [28]. (© 2023 IEEE. Reprinted, with permission, from EEEIC/I&CPS Europe).
Energies 16 01668 g009
Figure 10. Inferred EMD Cost Function C ρ i versus Distance from the Target Model Parameter [28]. (© 2023 IEEE. Reprinted, with permission, from EEEIC/I&CPS Europe).
Figure 10. Inferred EMD Cost Function C ρ i versus Distance from the Target Model Parameter [28]. (© 2023 IEEE. Reprinted, with permission, from EEEIC/I&CPS Europe).
Energies 16 01668 g010
Figure 11. Weighted Inferred EMD Cost Function C ρ w i versus Distance from Target Parameter.
Figure 11. Weighted Inferred EMD Cost Function C ρ w i versus Distance from Target Parameter.
Energies 16 01668 g011
Table 1. Parameter values for the model shown in Figure 1.
Table 1. Parameter values for the model shown in Figure 1.
R s C s C g L s K 12 K 13 K 23
83 m Ω 970 p F 1.6 n F 181 μ H 0.8792 0.8729 0.8564
Table 2. Parameter Boundary Constraints.
Table 2. Parameter Boundary Constraints.
Parameters R s C s C g L s K 12 K 13 K 23
Lower Bound1111858585
Upper Bound1000100010001000959595
Table 3. Estimated Parameter Error Percentages obtained through the Frequency-Domain Approach.
Table 3. Estimated Parameter Error Percentages obtained through the Frequency-Domain Approach.
ApproachFrequency-Domain
Cost Function C rmse z C ρ z
R s ( err % ) −242.5045.23
C s ( err % ) −3.09−0.50
C g ( err % ) 0.6737.50
L s ( err % ) 2.55−40.25
K 12 ( err % ) −8.83−3.73
K 13 ( err % ) 0.03-8.80
K 23 ( err % ) −5.29−6.70
Runtime (h)29.8811.81
Table 4. Estimated Model Parameter Error Percentages obtained through the Time Domain Approaches.
Table 4. Estimated Model Parameter Error Percentages obtained through the Time Domain Approaches.
ApproachTime-DomainEmpirical Mode DecompositionInferred Empirical Model Decomposition
AlignmentStrategy 1Strategy 2Strategy 1Strategy 2Strategy 1Strategy 2
Cost Function C rmse v C ρ v C rmse v C ρ v C nrmse m C ρ m C nrmse m C ρ m C nrmse i C ρ i C nrmse i C ρ i
R s ( err % ) 1.1483.13−392.754.6236.90−941.37−575.99−395.221.14−84.61−392.63−101.221
C s ( err % ) −1.9147.17−2.01−0.93−3.0325.25−2.30−3.03−1.9121.580.9910.62
C g ( err % ) −0.59−723.111.308.31−59.5537.50−33.1934.50−0.59−3.74−0.17−1.45
L s ( err % ) 1.7386.261.90−8.8033.47−51.2226.99−74.001.73−0.911.771.47
K 12 ( err % ) −4.78−4.85−3.602.31−5.030.42−5.03−5.66−4.780.67−1.15−2.24
K 13 ( err % ) −3.15−5.452.631.152.50−3.10−7.940.33−3.151.88−3.642.43
K 23 ( err % ) 0.53−8.03−8.63−4.93−0.50−0.41−6.460.650.530.74−3.66−6.36
Runtime (h)12.324.316.666.225.921.7614.659.078.338.876.729.49
Table 5. Estimated Model Parameter Error Percentages obtained through the Weighted IEMD approach.
Table 5. Estimated Model Parameter Error Percentages obtained through the Weighted IEMD approach.
Cost Function C ρ wi
AlignmentStrategy 1Strategy 2
R s ( err % ) –52.55–133.64
C s ( err % ) 17.2310.82
C g ( err % ) –3.28–2.11
L s ( err % ) 0.0381.62
K 12 ( err % ) –2.35–6.63
K 13 ( err % ) 0.240.54
K 23 ( err % ) 0.73–2.99
Runtime (h)297.68487.35
Table 6. RMSE Metrics to Cross-Validate the Performance of the Parameter Estimation Results presented in Table 3, Table 4 and Table 5 using an Impulse Response Test.
Table 6. RMSE Metrics to Cross-Validate the Performance of the Parameter Estimation Results presented in Table 3, Table 4 and Table 5 using an Impulse Response Test.
TestImpulse Response
AlignmentStrategy 1Strategy 2
C r m s e z 6.847 × 10 7
C ρ z 2.104 × 10 8
C r m s e v 1.002 × 10 8 1.205 × 10 8
C ρ v 3.114 × 10 8 2.030 × 10 8
C n r m s e m 1.646 × 10 8 1.867 × 10 8
C ρ m 2.535 × 10 8 6.648 × 10 8
C n r m s e i 1.002 × 10 8 9.874 × 10 7
C ρ i 9.173 × 10 7 7.723 × 10 7
C ρ w i 9.447 × 10 7 6.784 × 10 7
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

Banks, D.M.; Bekker, J.C.; Vermeulen, H.J. Investigating Empirical Mode Decomposition in the Parameter Estimation of a Three-Section Winding Model. Energies 2023, 16, 1668. https://doi.org/10.3390/en16041668

AMA Style

Banks DM, Bekker JC, Vermeulen HJ. Investigating Empirical Mode Decomposition in the Parameter Estimation of a Three-Section Winding Model. Energies. 2023; 16(4):1668. https://doi.org/10.3390/en16041668

Chicago/Turabian Style

Banks, Daniel Marc, Johannes Cornelius Bekker, and Hendrik Johannes Vermeulen. 2023. "Investigating Empirical Mode Decomposition in the Parameter Estimation of a Three-Section Winding Model" Energies 16, no. 4: 1668. https://doi.org/10.3390/en16041668

APA Style

Banks, D. M., Bekker, J. C., & Vermeulen, H. J. (2023). Investigating Empirical Mode Decomposition in the Parameter Estimation of a Three-Section Winding Model. Energies, 16(4), 1668. https://doi.org/10.3390/en16041668

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