Next Article in Journal
Variable Amplitude Gate Voltage Synchronous Drive Technique for Improving Dynamic Current Balancing in Paralleled IGBTs
Previous Article in Journal
Looking for an Optimal Composition of Nickel-Based Catalysts for CO2 Methanation
Previous Article in Special Issue
Study on Top Hierarchy Control Strategy of AEBS over Regenerative Brake and Hydraulic Brake for Hub Motor Drive BEVs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Power Measurement Using Adaptive Chirp Mode Decomposition for Electrical Vehicle Charging Load

1
Marketing Service Center (Measurement Center), Ningxia Electric Power Co., Ltd., National Grid of China, Yinchuan 750001, China
2
The College of Electrical Engineering, Sichuan University, Chengdu 610065, China
*
Author to whom correspondence should be addressed.
Energies 2023, 16(14), 5305; https://doi.org/10.3390/en16145305
Submission received: 23 May 2023 / Revised: 2 July 2023 / Accepted: 6 July 2023 / Published: 11 July 2023
(This article belongs to the Special Issue Energy Efficiency of Electric Vehicle Charging Process)

Abstract

:
Due to nonlinear components in the charging piles of electric vehicles, harmonics and nonstationary signals in the electric vehicle charging load bring voltage and current distortion, seriously affecting the accuracy of the power-related calculation in nonsinusoidal environments. This paper proposed a new approach to calculate the active power and root mean square values from decomposed components using the adaptive chirp mode decomposition (ACMD) method on voltage and current. The advantage of the ACMD-based method is that it correctly provides the power-related quantities of harmonics or nonstationary components for the electric vehicle charging load. The performance of the proposed method was verified using synthetic signals and simulation tests. The experimental results presented better estimations for each quantity defined in IEEE Standard 1459-2010, compared with the discrete wavelet transform approach.

1. Introduction

With the enhancement of environmental protection awareness, electric vehicles have been widely used in many countries to alleviate global warming and energy crisis [1,2]. Typically, electric vehicles obtain electrical energy from the power grid through charging piles, especially direct current (DC) charging piles, whereas due to nonlinear components contained in the batteries of the electric vehicles and charging piles, the charging behavior of the electric vehicles, especially fast charging mode, brings a large number of harmonic and nonstationary signals, resulting in voltage and current distortion on the grid side of the charging piles, i.e., the voltage or current waveform is not standard sinusoidal [3]. The nonsinusoidal waveform would seriously affect the accuracy of electricity metering; thus, the research on the power calculation for electric vehicles makes sense for fair trade between users and power companies, as well as the generation, transmission, and operation of electricity.
For the power calculation under nonsinusoidal conditions, Budeanu [4] and Fryze [5] first proposed the definition of power quantities, commonly referred to as power theories, based on frequency-domain or time-domain analysis, respectively. After that, the definitions of various power components have been developed by many researchers [6,7,8,9,10,11]. Among these, IEEE Standard 1459-2010 [12], as one of the most prominent power theories, provides definitions mainly based on Fourier transform analysis. Although Fourier transform can provide the frequency information of analytical signals, it is only suitable to analyze the stationary and periodic signals due to the missing time-domain information. Commonly, harmonics, sudden changes, and nonstationary signals are generated during the charging processing of electric vehicles [13]. For Fourier transform, it is not possible to extract the time–frequency information of nonstationary signals; thus, it would bring large errors in power calculation.
To address the limitation of Fourier transform, wavelet-related transforms, as one of the most widely used time–frequency analyses, have been applied to power calculation in nonsinusoidal environments. Discrete wavelet transform (DWT) is first applied to calculate the root mean square (RMS) and power values by decomposing the given signals into constituent frequency sub-bands [14,15]. According to the definition of IEEE Standard 1459-2010 [12], the DWT was utilized to estimate the power components in a single system [14] and a three-phase system [16]. Further, for maintaining the temporal resolution of each decomposition level, stationary discrete wavelet transform (SDWT) without downsampling was used for the power calculation [17]. To achieve uniform frequency bandwidth for each sub-band, discrete wavelet packet transform (DWPT) [18] or stationary DWPT [19] was proposed for power calculation.
Although wavelet-related transforms can decompose the signals into frequency sub-bands through a filtering process, there is the issue of energy leakage [20], i.e., the content of decomposed frequency sub-bands does not correspond to individual harmonics or special nonstationary signal components. If the sub-bands are used for the power calculation of individual harmonic components or special nonstationary signal components, the results cause inconsistencies from the real power. This limitation of wavelet-related transforms is not present in adaptive chirp mode decomposition (ACMD). The ACMD extracts signal modes one by one with good adaptability and high time–frequency resolution, widely used in many fields [21,22,23,24]. However, the ACMD has not been applied for power calculation in the literature.
Thus, according to the definition of the IEEE Standard 1459-2010 [12], this paper presents an ACMD-based method to calculate power-related quantities of the charging load of the electric vehicles through the charging piles, since the electric vehicles connect to the power grid through the charging piles. Based on a time–frequency representation of the voltage and current using the ACMD method, the power-related quantities are calculated for the fundamental, harmonic, and nonstationary components, respectively. The proposed method is suitable for estimating harmonics or nonstationary components of voltage and current signals on the grid side of the charging piles, improving the representation of power-related quantities in nonstationary conditions.

2. Proposed Method

This section first presents the general mathematical mode of the voltage and current signals, according to the principle of the charging piles of electric vehicles. And then, the ACMD method is introduced to decompose the voltage and current signals into each component. Finally, decomposed components are utilized to estimate the corresponding active power and the RMS value of the voltage and current signals for the electric vehicle charging load.

2.1. General Mathematical Mode of Charging Voltage and Current

Due to the wide use of power electronic devices in the charging piles, the charging for the batteries of electric vehicles brings harmonic current and nonstationary signals to the power grid. The general structure of the charging piles is mainly composed of a rectifier circuit, a power converter, and an output filter circuit, as shown in Figure 1. Alternating current (AC) power is rectified to DC power through a three-phase bridge uncontrollable rectifier circuit and then provided for the DC/DC power converter to output the required DC of the electric vehicles. Finally, DC is filtered to charge the batteries of the electric vehicles.
More studies on the load characteristics of electric vehicles show that the voltage and current waveforms on the grid side of the charging piles, collected at charging stations [13,25], contain multiple harmonics and nonstationary signals. Typically, the current waveform mainly includes 6 h ± 1 order harmonics, and the amplitude of the harmonics gradually decreases with the increase in the harmonic order, which further can be confirmed in the simulation model in Section 3.2.
Generally, the voltage and current on the grid side of the charging piles can be respectively expressed as
u ( t ) = u b ( t ) + h = 1 u h ( t ) + u s ( t ) = 2 U b sin ( ω t + φ b ) + h = 1 2 U h sin ( ( 6 h ± 1 ) ω t + φ h ) + u s ( t )
and
i ( t ) = i b ( t ) + h = 1 i h ( t ) + i s ( t ) = 2 I b sin ( ω t + ϕ b ) + h = 1 2 I h sin ( ( 6 h ± 1 ) ω t + ϕ h ) + i s ( t ) ,
where ( u b , u h , u s ) denote the fundamental, harmonic, and nonstationary components of the voltage signal, respectively; ( i b , i h , i s ) denote the fundamental, harmonic, and nonstationary ones of the current signal, respectively; ω is the fundamental frequency in radians; and ( φ b , φ h , ϕ b , ϕ h ) are the initial phases of the fundamental voltage, the h-th harmonic voltage, the fundamental current, and the h-th harmonic current, respectively.
To accurately calculate the charging power of charging piles for the electric vehicles, the ACMD is firstly applied to decompose the fundamental, harmonic, and nonstationary components from u ( t ) and i ( t ) in our work, and then the power-related quantities are calculated according to IEEE Standard 1459-2010 [12].

2.2. ACMD-Based Power-Related Calculation

2.2.1. Adaptive Chirp Mode Decomposition (ACMD)

The ACMD, as a recursive mode decomposition technique, can adaptively extract signal components one by one using the matching pursuit method. For a nonstationary signal x ( t ) with m components, it can be modeled as
x ( t ) = i = 1 m x i ( t ) = i = 1 m A i ( t ) cos 2 π 0 t f i ( τ ) d τ + φ i ,
where A i ( t ) , f i ( t ) , and φ i denote the instantaneous amplitude, instantaneous frequency, and the initial phase of the i-th component, respectively. Based on the technique of frequency demodulation, a demodulation frequency f ˜ i ( t ) is introduced, and then Equation (3) can be reformulated as
x ( t ) = i = 1 m p i ( t ) cos 2 π 0 t f ˜ i ( τ ) d τ + q i ( t ) sin 2 π 0 t f ˜ i ( τ ) d τ ,
with
p i ( t ) = A i ( t ) cos 2 π 0 t ( f i ( τ ) f ˜ i ( τ ) ) d τ + φ i q i ( t ) = A i ( t ) sin 2 π 0 t ( f i ( τ ) f ˜ i ( τ ) ) d τ + φ i ,
where cos ( 2 π 0 t f ˜ i ( τ ) d τ ) and sin ( 2 π 0 t f ˜ i ( τ ) d τ ) denote two demodulation operators, and p i ( t ) and q i ( t ) denote two demodulated signals with a phase difference of 90 . In particular, p i ( t ) and q i ( t ) become slowly varying baseband signals and have the narrowest bandwidth when f ˜ i ( t ) = f i ( t ) .
To extract the i-th signal component x i ( t ) , the ACMD solve the minimization problem as follows
min { p i ( t ) } , { q i ( t ) } , { f i ( t ) } p i ( t ) 2 2 + q i ( t ) 2 2 + γ x ( t ) x i ( t ) 2 2 , s . t . x i ( t ) = p i ( t ) cos 2 π 0 t f ˜ i ( τ ) d τ + q i ( t ) sin 2 π 0 t f ˜ i ( τ ) d τ
where · 2 denotes L 2 norm, ( · ) denotes the second derivative, two terms p i ( t ) 2 2 and q i ( t ) 2 2 are used as smoothness constraints for measuring bandwidth of the signal, γ > 0 is the penalty factor, and the last term x ( t ) x i ( t ) 2 2 denotes the residue energy after the currently estimated i-th component is eliminated.
To implement ACMD decomposition in a computer, the signal x ( t ) is sampled at discrete time t = t 0 , . . . , t N 1 , where N is the number of samples. And then, the objective function of (6) can be written in a discrete form
L ( p i , q i , f i ) = D p i 2 2 + D q i 2 2 + γ x ( p i c i + q i s i ) 2 2 ,
where D is a second-order difference matrix, given by
D = 1 2 1 0 0 1 2 1 0 0 1 2 1 ( N 2 ) × N ,
and
c i = diag [ cos ( θ i ( t 0 ) ) , . . . , cos ( θ i ( t N 1 ) ) ] s i = diag [ sin ( θ i ( t 0 ) ) , . . . , sin ( θ i ( t N 1 ) ) ] ,
with θ i ( t ) = 2 π 0 t f ˜ i ( τ ) d τ , where d i a g [ · ] denotes the diagonal matrix.
To solve the optimization problem in (7), an iterative algorithm is used to alternately update the demodulated signals p i , q i , and the frequency f i . For the m-th iteration, with the current frequency f i ( m 1 ) , the solutions of p i and q i are updated by setting L ( p i , q i , f i ) / p i = 0 and L ( p i , q i , f i ) / q i = 0 , given by
p i ( m ) = 1 γ D T D + ( c i ( m ) ) T c i ( m ) 1 ( c i ( m ) ) T x q i ( m ) = 1 γ D T D + ( s i ( m ) ) T s i ( m ) 1 ( s i ( m ) ) T x ,
respectively, and then, the instantaneous frequency f i is updated by
f i ( m + 1 ) = f i ( m ) + 2 β D T D + I 1 Δ f ˜ i ( m ) ,
with
Δ f ˜ i ( m ) ( t ) = 1 2 π d d t tan 1 p i ( m ) ( t ) q i ( m ) ( t ) = q i ( m ) ( t ) p i ( m ) ( t ) p i ( m ) ( t ) q i ( m ) ( t ) 2 π p i ( m ) ( t ) 2 + q i ( m ) ( t ) 2 ,
where β > 0 is a weighting factor, and ( · ) denotes the derivative.
Since the iterative algorithm requires the initial frequency f i ( 0 ) in advance, synchroextracting transform (SET) [26] is adopted to estimate f i ( 0 ) in this study. The SET, as a novel postprocessing strategy of short-time Fourier transform, has the advantage of better energy concentration to provide good instantaneous initial frequency for iteration. The iterative process continues until there is little difference in the extracted signal component x i in two adjacent iterations. After the signal component x i is extracted from the raw signal, the residual signal is obtained by removing x i from the raw signal, and then the iteration can be applied to the residual signal to extract the next signal component sequentially.
The above procedure is repeated for M times, and the components of signal x ( t ) are finally obtained as
x ( t ) = k = 1 M x k ( t ) + R s ( t ) ,
where R s ( t ) denotes the residual signal after M components are separated from the original signal. More details of the iteration algorithm can be found in [27].
When the ACMD method is applied to the voltage and current signals of charging piles, the corresponding fundamental component, harmonic components, and nonstationary components are separated from each other, and then, the decomposed components are used to implement the calculation of the active power, the RMS values of the voltage and current.

2.2.2. Power Calculation

According to IEEE Standard 1459-2010 [12], the active power (P), the RMS values of the voltage (U) and current (I) are calculated as
P = 1 K T t = t 0 t 0 + K T u ( t ) × i ( t ) d t ,
U = 1 K T t = τ τ + K T u 2 ( t ) d t ,
and
I = 1 K T t = τ τ + K T i 2 ( t ) d t ,
respectively, where τ and τ + K T are the start time and end time of measurements, respectively, and k T is an integer multiple of the power system fundamental period.
Using u ( t ) in (1) and i ( t ) in (2), the active power P in (12) is expressed as
P = P b b + P b h + P b s + P h b + P h h + P h s + P s b + P s h + P s s
with
P b b = 1 K T t = t 0 t 0 + K T u b ( t ) i b ( t ) d t , P h h = 1 K T h = 1 t = t 0 t 0 + K T u h ( t ) i h ( t ) d t P s s = 1 K T t = t 0 t 0 + K T u s ( t ) i s ( t ) d t , P b h = 1 K T h = 1 t = t 0 t 0 + K T u b ( t ) i h ( t ) d t , P h b = 1 K T h = 1 t = t 0 t 0 + K T u h ( t ) i b ( t ) d t , P b s = 1 K T t = t 0 t 0 + K T u b ( t ) i s ( t ) d t , P s b = 1 K T t = t 0 t 0 + K T u s ( t ) i b ( t ) d t , P h s = 1 K T h = 1 t = t 0 t 0 + K T u h ( t ) i s ( t ) d t , P s h = 1 K T h = 1 t = t 0 t 0 + K T u s ( t ) i h ( t ) d t ,
where P b b , P h h , and P s s are the fundamental, harmonic, and nonstationary active power, respectively; and P b h , P h b , P b s , P s b , P h s , and P s h are the distortion active powers, generated by the different cross-combinations from the fundamental/harmonic/nonstationary voltage and the fundamental/harmonic/nonstationary current, respectively.
Since
0 2 π sin ( n ω t ) sin ( m ω t ) d t = 0
for n m , we have P b h = P h b = 0 . And then, P b s 0 , P s b 0 , P h s 0 , and P s h 0 can be considered in statistics, although there are minor errors. Thus, according to References [25,28], ( P b s , P s b , P h s , and P s h ) can be ignored, and then, the simplified calculation of the active power in (15) is given by
P = P b b + P h h + P s s ,
i.e., the active power is the summation of the fundamental, harmonic, and nonstationary active power.
Accordingly, we calculate the RMS values of the voltage and current as
U = U b 2 + U h 2 + U s 2 , I = I b 2 + I h 2 + I s 2 ,
where ( U b , U h , and U s ) are the RMS values of the fundamental, harmonic, and nonstationary voltage, respectively, and ( I b , I h , and I s ) are the RMS values of the fundamental, harmonic, and nonstationary current, respectively.
In practice, the analog voltage and current are first digitalized, and the continuous integration operation would be converted into the discrete ones. Thus, the active power ( P b b , P h h , and P s s ) in (17) can be calculated as
P b b = 1 N n = 0 N 1 u b ( n ) × i b ( n ) , P h h = 1 N h = 1 n = 0 N 1 u h ( n ) × i h ( n ) , P s s = 1 N n = 0 N 1 u s ( n ) × i s ( n ) ,
where ( u b , h , s ( n ) , i b , h , s ( n ) ) are the digitalized data of ( u b , h , s ( t ) , i b , h , s ( t ) ) at the n-th sampling index, respectively, and N is the integration interval. Accordingly, the RMS values of the voltage and current in (18) can be calculated as
U b = 1 N n = 0 N 1 u b 2 ( n ) , I b = 1 N n = 0 N 1 i b 2 ( n ) , U h = 1 N h = 1 n = 0 N 1 u h 2 ( n ) , I h = 1 N h = 1 n = 0 N 1 i h 2 ( n ) , U s = 1 N n = 0 N 1 u s 2 ( n ) , I s = 1 N n = 0 N 1 i s 2 ( n ) ,
for the fundamental, harmonics, and nonstationary components, respectively.

3. Experiments and Results

In this section, we evaluate the proposed ACMD-based estimation of RMS values and active power using synthetic signals and simulation tests. In the following experiments, the key parameters γ and β of the ACMD method are set as γ = 10 6 and β = 10 9 , respectively. The performance of the ACMD-based method is compared with the DWT-based method, since DWT is a time–frequency method widely used in signal analysis for power systems. The relative error of the estimated quantity is calculated as
E r r o r = | R e f E s t | R e f ,
where R e f is the reference value, and E s t is the estimated value from the ACMD or DWT method.

3.1. Nonstationary Synthetic Signal

To demonstrate the performance of the ACMD method in signal decomposition and power-related estimation, three stationary waveforms, containing the fundamental and two harmonics, and one nonstationary component, are considered in synthetic signals of the voltage and current, expressed as
u ( t ) = u 1 ( t ) + u 2 ( t ) + u 3 ( t ) + u 4 ( t ) , = 2 220 sin ( 2 π × 50 t ) + 20 sin ( 2 π × 250 t + 2 π ) + 5 sin ( 2 π × 350 t + π ) + e 0.3 t cos ( 2 π ( 450 t + 0.5 π cos ( 2 π × 30 t ) ) ) ) ,
i ( t ) = i 1 ( t ) + i 2 ( t ) + i 3 ( t ) + i 4 ( t ) , = 2 25 sin ( 2 π × 50 t ) + 10 sin ( 2 π × 250 t ) + 5 sin ( 2 π × 350 t ) + 2 e 0.3 t cos ( 2 π ( 450 t + 0.5 π cos ( 2 π × 30 t ) ) ) ) ,
where the fundamental frequency is 50 Hz, the sampling frequency is 1500 Hz for discretizing i ( t ) and u ( t ) , and the time interval of the synthetic signal is 0.5 s in this experiment.
To visually evaluate the decomposition performance of ACMD and DWT methods, we show the resultant synthetic current waveform and each component in Figure 2a as the ground-truth ones, based on (22). After the DWT with mother wavelet ‘db4’ and ACMD methods are used to decompose the current i ( t ) , the decomposed components and reconstructed signal are shown in Figure 2b,c. According to the principle of the DWT method, the fundamental component can be extracted from three-level DWT decomposition, and its frequency is within that of the highest-level approximation coefficients, i.e., ca3 in Figure 2b. From the decomposed results of the DWT method shown in Figure 2b, we observe that there is a component that does not already exist in (22) within [93.75, 187.5] Hz. In addition, the obtained components within [0, 93.75] Hz are not pure fundamental components as expected but mixed with the leakage components from other frequency ranges. Similar patterns exist in other frequency ranges, whereas the ACMD method obtains each component, more closely to the ground-truth ones, compared with the results from the DWT method.
The calculated quantities and relative errors (‰) of the fundamental, harmonics, and nonstationary components are shown in Table 1, where the reference values are directly calculated from the analytic expressions in (22), agreeing with IEEE Standard 1459-2010. The estimated quantities for the three-type components are calculated from the decomposed results from the ACMD and DWT methods. The results in Table 1 show that the DWT-based method causes small estimations for the fundamental and large ones for harmonic, especially for nonstationary components, resulting in bringing large estimated errors for each component. The main reason is its leakage issue as mentioned above, whereas compared with the DWT-based method, the ACMD-based method brings quantities closer to the reference values with fewer estimation errors.

3.2. Simulation Signal

The performance of the proposed ACMD-based calculation method is further validated by using a simulation model in MATLAB/SIMULINK, as shown in Figure 3. The current and voltage waveforms at the grid side of the charging pile are shown in Figure 4a,b, respectively. From Figure 4, we observe that the charging current is similar to a square wave, deviating from the standard sine signal, while the voltage contains fewer harmonics, closer to the sine signal. Further, considering there are nonlinear elements in the charging piles of electric vehicles, the nonstationary signals are introduced to the current and voltage for realistically simulating the charging characteristics of electric vehicles. In this experiment, sawtooth impulse and high-frequency attenuation signals are used for the nonstationary signals u s and i s , expressed as
u s ( t ) = 8 × s a w t o o t h ( 2 π × 1800 t + 0.5 π , 0.5 ) , t ( 0.4 , 0.55 ) s , 8 × e 40 t sin ( 2 π × 1500 t ) , t ( 0.75 , 0.9 ) s ,
and
i s ( t ) = 10 × s a w t o o t h ( 2 π × 1800 t + 0.5 π , 0.5 ) , t ( 0.4 , 0.55 ) s , 10 × e 40 t sin ( 2 π × 1500 t ) , t ( 0.75 , 0.9 ) s ,
respectively. Figure 4c shows the waveform of i s ( t ) . To observe the components of current, the spectrum of the current i ( t ) is shown in Figure 4d, which demonstrates that 6 h ± 1 order harmonics with successively attenuated amplitude appear on the corresponding spectral bins beside the fundamental, and the nonstationary signals appear on multiple high-frequency ranges, as shown in the enlarged drawing.
The sampling frequency is 4000 Hz in this experiment, and five-level DWT decomposition with the mother wavelet ‘db4’ is applied on current and voltage, making it that the highest-level frequency range corresponds to that of the fundamental component (i.e., ca5). Figure 5 shows the decomposed results of the current. From the decomposed fundamental and nonstationary components, it is easy to discern that both of them are mixed with other frequency components. Meanwhile, the ACMD method is utilized for current and voltage, and the decomposed results for the current are shown in Figure 6. Since the ACMD method performs decomposition recursively, two nonstationary signals are extracted separately, and close to the raw waveform in the views of waveform trend and occurrence interval. For harmonics, parts of the waveform within black windows are enlarged and displayed in small attached windows for comparison. Compared with the results from the DWT method in Figure 5, the decomposed components from the ACMD method in Figure 6 more agree with the real ones.
Using the decomposed components from the DWT and ACMD methods, the RMS values, active power, and relative error (‰) for the three components are calculated and shown in Table 2. Since the real value for each component can not be calculated directly as the case for synthetic signals, the RMS values of the fundamental and harmonics are calculated from the corresponding information of spectral magnitude, while the quantities of the nonstationary component are calculated according to u s and i s ( t ) in (24) and (25). Similar to the comparison results for synthetic signals, the ACMD-based method achieves better estimation for each component and total active power for simulation data compared with the DWT-based method.

4. Conclusions

In this paper, a new approach based on the ACMD method was proposed to estimate the power-related quantities for electric vehicle charging load, according to the definition in the IEEE Standard 1459-2010 in the time–frequency domain. The current and voltage on the grid side of the charging piles were decomposed into the fundamental correctly, harmonic, and nonstationary components through the ACMD method. With the decomposed components, the active power and the RMS values were calculated for each component and the overall. Compared with the DWT-based method, the proposed ACMD-based method achieved better estimation of power-related quantities for the electric vehicle charging load with significantly fewer estimation errors, especially for the harmonics and nonstationary components. To further explore the advantage of the ACMD method for more applications in power systems, our future works will focus on power-related estimation of three-phase systems or other electric devices (e.g., power transformer) for abnormality detection, fault diagnosis, etc.

Author Contributions

Conceptualization, H.D. and R.T.; methodology, H.D. and X.Y.; validation, R.T.; investigation, R.T.; writing—original draft preparation, H.D. and J.W.; writing—review and editing, X.Y.; visualization, J.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Technology Research Project of Ningxia Electric Power Co., Ltd. National Grid of China, “Research on the Influence of Electric Vehicle Charging Load on the Accuracy of Electric Energy Measurement” (SGTYHT/21-JS, SGNXYX002JJS2200153).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

PActive power
URMS value of the voltage
IRMS value of the current
u ( t ) and i ( t ) ) Voltage and current on the grid side of the charging piles, respectively
u b ( t ) and i b ( t ) Fundamental voltage and current, respectively
u h ( t ) and i h ( t ) h-th harmonic voltage and current, respectively
u s ( t ) and i s ( t ) Nonstationary voltage and current, respectively
P b b , P h h , and P s s Fundamental, harmoinc, and nonstationary active power, respectively
U b , U h , and U s RMS values of the fundamental, harmoinc, and nonstationary voltage, respectively
I b , I h , and I s RMS values of the fundamental, harmoinc, and nonstationary current, respectively

References

  1. Zhao, H.; Wang, L.; Chen, Z.; He, X. Challenges of Fast Charging for Electric Vehicles and the Role of Red Phosphorous as Anode Material: Review. Energies 2019, 12, 3897. [Google Scholar] [CrossRef] [Green Version]
  2. Liu, S.; He, K.; Pan, X.; Hu, Y. Review of Development Trend of Transportation Energy System and Energy Usages in China Considering Influences of Intelligent Technologies. Energies 2023, 16, 4142. [Google Scholar] [CrossRef]
  3. Ding, X.; Shi, H.; Wang, Y.; Zhuang, Y.; Yuan, G.; Zhu, S. Research on Harmonic Management of Single-Phase AC Charging Pile Based on Active Filtering. Energies 2023, 16, 2817. [Google Scholar] [CrossRef]
  4. Budeanu, C. Puissances Ractives et Fictives; Institut National Romain pour l’Étude de l’Aménagement et de l’Utilisation d’Énergie: Bucharest, Romania, 1927. [Google Scholar]
  5. Fryze, S. Active, reactive and apparent power in circuts with nonsinusoidal voltages and currents. Przglad Elektro. 1931, 1, 193–203. [Google Scholar]
  6. Miyasaka, G.; Silvério, E.; Xavier, G.; da Silva, H.; Braz, L.; Oliveira, R.; Lima, R.; Macedo, J. Analysis of Reactive Energy Measurement Methods Under Non-Sinusoidal Conditions. IEEE Lat. Am. Trans. 2018, 16, 2521–2529. [Google Scholar] [CrossRef]
  7. Pajic, S.; Emanuel, A.E. Modern apparent power definitions: Theoretical versus practical Approach-the general case. IEEE Trans. Power Deliv. 2006, 21, 1787–1792. [Google Scholar] [CrossRef]
  8. George, S.; Agarwal, V. A novel, DSP based algorithm for optimizing the harmonics and reactive power under non-sinusoidal supply voltage conditions. IEEE Trans. Power Deliv. 2005, 20, 2526–2534. [Google Scholar] [CrossRef]
  9. De Almeida Coelho, R.; Brito, N.S.D. Power Measurement Using Stockwell Transform. IEEE Trans. Power Deliv. 2021, 36, 3091–3100. [Google Scholar] [CrossRef]
  10. Morsi, W.G.; El-Hawary, M.E. Defining Power Components in Nonsinusoidal Unbalanced Polyphase Systems: The Issues. IEEE Trans. Power Deliv. 2007, 22, 2428–2438. [Google Scholar] [CrossRef]
  11. Lev-Ari, H.; Stankovic, A.M. A decomposition of apparent power in polyphase unbalanced networks in nonsinusoidal operation. IEEE Trans. Power Syst. 2006, 21, 438–440. [Google Scholar] [CrossRef]
  12. IEEE Std 1459-2010 (Revision of IEEE Std 1459-2000); IEEE Standard Definitions for the Measurement of Electric Power Quantities under Sinusoidal, Nonsinusoidal, Balanced, or Unbalanced Conditions. IEEE: New York, NY, USA, 2010; pp. 1–50.
  13. Liu, J.; Xu, Q.; Tian, Z.; Guo, Y.; Qi, S.; Wang, L. Feature Extraction Method for DC Charging Signal of Electric Vehicle. In Proceedings of the 3rd Annual International Conference on Information System and Artificial Intelligence (ISAI2018), Hong Kong, China, 24–26 June 2016. [Google Scholar]
  14. Morsi, W.G.; El-Hawary, M.E. Reformulating Power Components Definitions Contained in the IEEE Standard 1459-2000 Using Discrete Wavelet Transform. IEEE Trans. Power Deliv. 2007, 22, 1910–1916. [Google Scholar] [CrossRef]
  15. Yoon, W.K.; Devaney, M. Reactive power measurement using the wavelet transform. IEEE Trans. Instrum. Meas. 2000, 49, 246–252. [Google Scholar] [CrossRef] [Green Version]
  16. Morsi, W.G.; El-Hawary, M.E. Reformulating Three-Phase Power Components Definitions Contained in the IEEE Standard 1459-2000 Using Discrete Wavelet Transform. IEEE Trans. Power Deliv. 2007, 22, 1917–1925. [Google Scholar] [CrossRef]
  17. Morsi, W.; El-Hawary, M. A new perspective for the IEEE standard 1459–2000 via stationary wavelet transform in the presence of nonstationary power quality disturbance. In Proceedings of the 2009 IEEE Power & Energy Society General Meeting, Calgary, AB, Canada, 26–30 July 2009; p. 1. [Google Scholar]
  18. Osipov, D.S.; Dolgikh, N.N.; Goryunov, V.N.; Kovalenko, D.V. Algorithms of packet wavelet transform for power determination under nonsinusoidal modes. In Proceedings of the 2016 Dynamics of Systems, Mechanisms and Machines (Dynamics), Omsk, Russia, 15–17 November 2016; pp. 1–5. [Google Scholar]
  19. Alves, D.K.; Costa, F.B.; de Araujo Ribeiro, R.L.; de Sousa Neto, C.M.; Rocha, T.D.O.A. Real-Time Power Measurement Using the Maximal Overlap Discrete Wavelet-Packet Transform. IEEE Trans. Instrum. Meas. 2017, 64, 3177–3187. [Google Scholar] [CrossRef]
  20. Peng, Z.K.; Jackson, M.R.; Rongong, J.A.; Chu, F.L.; Parkin, R.M. On The Energy Leakage Of Discrete Wavelet Transform. Mech. Syst. Signal Process. 2009, 23, 330–343. [Google Scholar] [CrossRef] [Green Version]
  21. Wang, X.; Tang, G.; Yan, X.; He, Y.; Zhang, X.; Zhang, C. Fault Diagnosis of Wind Turbine Bearing Based on Optimized Adaptive Chirp Mode Decomposition. IEEE Sens. J. 2021, 21, 13649–13666. [Google Scholar] [CrossRef]
  22. Zhang, X.; Liu, Z.; Kong, Y.; Li, C. Mutual Interference Suppression Using Signal Separation and Adaptive Mode Decomposition in Noncontact Vital Sign Measurements. IEEE Trans. Instrum. Meas. 2022, 71, 1–15. [Google Scholar] [CrossRef]
  23. Wang, H.; Chen, S.; Zhai, W. Data-driven adaptive chirp mode decomposition with application to machine fault diagnosis under non-stationary conditions. Mech. Syst. Signal Process. 2023, 188, 109997. [Google Scholar] [CrossRef]
  24. Ding, C.; Wang, B. Sparsity-assisted adaptive chirp mode decomposition and its application in rub-impact fault detection. Measurement 2022, 188, 110539. [Google Scholar] [CrossRef]
  25. Qi, S.R. Research on the Influence of Electrical Vehicle Charging Load on Electric Energy Measurement Accuracy. Master’s Thesis, Southeast University, Nanjing, China, 2019. [Google Scholar]
  26. Yu, G.; Yu, M.; Xu, C. Synchroextracting Transform. IEEE Trans. Ind. Electron. 2017, 64, 8042–8054. [Google Scholar] [CrossRef]
  27. Chen, S.; Yang, Y.; Peng, Z.; Dong, X.; Zhang, W.; Meng, G. Adaptive chirp mode pursuit: Algorithm and applications. Mech. Syst. Signal Process. 2019, 116, 566–584. [Google Scholar] [CrossRef]
  28. Tong, T.; Ni, W.; Gu, D.; Cao, Y.; Yang, W. Study on the method of electric energy measurement for unsteady signal. In Proceedings of the 16th IET International Conference on AC and DC Power Transmission (ACDC 2020), Online, 2–3 July 2020; Volume 2020, pp. 1634–1640. [Google Scholar]
Figure 1. Topology of the charging pile of electric vehicles.
Figure 1. Topology of the charging pile of electric vehicles.
Energies 16 05305 g001
Figure 2. Current synthetic waveform and the corresponding decomposed results by using DWT and ACMD methods: (a) synthetic current and its components; (b) decomposed and reconstructed results of the DWT method; (c) decomposed and reconstructed results of the ACMD method.
Figure 2. Current synthetic waveform and the corresponding decomposed results by using DWT and ACMD methods: (a) synthetic current and its components; (b) decomposed and reconstructed results of the DWT method; (c) decomposed and reconstructed results of the ACMD method.
Energies 16 05305 g002
Figure 3. Simulation model in MATLAB/SIMULINK.
Figure 3. Simulation model in MATLAB/SIMULINK.
Energies 16 05305 g003
Figure 4. Waveform generated from simulation model and the spectrum of current: (a) current; (b) voltage; (c) nonstationary i s ; (d) the spectrum of the current.
Figure 4. Waveform generated from simulation model and the spectrum of current: (a) current; (b) voltage; (c) nonstationary i s ; (d) the spectrum of the current.
Energies 16 05305 g004
Figure 5. Decomposed results by using DWT method for current.
Figure 5. Decomposed results by using DWT method for current.
Energies 16 05305 g005
Figure 6. Decomposed results by using ACMD method for current.
Figure 6. Decomposed results by using ACMD method for current.
Energies 16 05305 g006
Table 1. Comparisons of estimated quantities of ACMD-based and DWT-based methods for synthetic signal.
Table 1. Comparisons of estimated quantities of ACMD-based and DWT-based methods for synthetic signal.
QuantityReference ValuesACMDDWT
SymbolsParametersAnalyticalValueError (‰)ValueError (‰)
UVoltage RMS (V)220.97220.820.68220.870.45
ICurrent RMS (A)27.3927.380.3627.380.36
PActive power (W)5625.865618.001.405602.854.09
U b Fundamental voltage RMS (V)220.00219.850.68217.969.27
I b Fundamental current RMS (A)25.0024.980.8024.808.00
P b b Fundamental active power (W)5500.005492.781.315399.4618.28
U h Harmonic voltage RMS (V)20.6220.591.4535.13703.68
I h Harmonic current RMS (A)11.1811.170.8910.8529.52
P h h Harmonic active power (W)125.00124.662.72186.88495.04
U s Nonstationary voltage RMS (V)1.311.07183.216.463931.29
I s Nonstationary current RMS (A)0.660.53196.964.145272.72
P s s Nonstationary active power (W)0.860.57337.2116.5218,209.30
Table 2. Comparisons of estimated quantities of ACMD-based and DWT-based methods for simulation signals in MATLAB/SIMULINK.
Table 2. Comparisons of estimated quantities of ACMD-based and DWT-based methods for simulation signals in MATLAB/SIMULINK.
Quantity Reference ValuesACMDDWT
Symbols        ParametersAnalyticalValueError (‰)ValueError (‰)
UVoltage RMS (V)224.76224.770.0223.883.92
ICurrent RMS (A)49.5449.621.6149.491.01
PActive power (W)10,756.4610,749.330.6610,659.489.02
U b Fundamental voltage RMS (V)224.41224.460.22200.07108.46
I b Fundamental current RMS (A)47.820047.830.2142.74106.23
P b b Fundamental active power (W)10,686.1110,691.480.508512.89203.36
U h Harmonic voltage RMS (V)10.4910.7524.78100.398570.06
I h Harmonic current RMS (A)12.6012.7411.1124.69959.52
P h h Harmonic active power (W)48.4851.0753.422139.6943,135.52
U s Nonstationary voltage RMS (V)2.282.4783.333.78657.89
I s Nonstationary current RMS (A)2.862.8213.983.69290.21
P s s Nonstationary active power (W)6.545.79114.486.8953.45
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

Ding, H.; Tian, R.; Wang, J.; Yang, X. Power Measurement Using Adaptive Chirp Mode Decomposition for Electrical Vehicle Charging Load. Energies 2023, 16, 5305. https://doi.org/10.3390/en16145305

AMA Style

Ding H, Tian R, Wang J, Yang X. Power Measurement Using Adaptive Chirp Mode Decomposition for Electrical Vehicle Charging Load. Energies. 2023; 16(14):5305. https://doi.org/10.3390/en16145305

Chicago/Turabian Style

Ding, Haili, Rui Tian, Jinfei Wang, and Xiaomei Yang. 2023. "Power Measurement Using Adaptive Chirp Mode Decomposition for Electrical Vehicle Charging Load" Energies 16, no. 14: 5305. https://doi.org/10.3390/en16145305

APA Style

Ding, H., Tian, R., Wang, J., & Yang, X. (2023). Power Measurement Using Adaptive Chirp Mode Decomposition for Electrical Vehicle Charging Load. Energies, 16(14), 5305. https://doi.org/10.3390/en16145305

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