Next Article in Journal
Modeling Magnetostrictive Transducers for Structural Health Monitoring: Ultrasonic Guided Wave Generation and Reception
Previous Article in Journal
Vehicle Trajectory Estimation Based on Fusion of Visual Motion Features and Deep Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Calibration of On-Board Energy Measurement Systems Installed in Locomotives for AC Distorted Current and High Voltage Waveforms and Determination of Its Uncertainty Budget

1
FFII-LCOE, Eric Kandel Street 1, Getafe, 28906 Madrid, Spain
2
Department of Electrical Engineering, Universidad Politécnica de Madrid, Rda. de Valencia 3, 28012 Madrid, Spain
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(23), 7967; https://doi.org/10.3390/s21237967
Submission received: 18 September 2021 / Revised: 22 November 2021 / Accepted: 23 November 2021 / Published: 29 November 2021
(This article belongs to the Section Electronic Sensors)

Abstract

:
Periodic calibrations of Energy Measurement Systems (EMS) installed in locomotives must be carried out to demonstrate the required accuracy established in the EN 50463-2 standard according to European Parliament and Council Directive 2008/57/EC on the interoperability of rail systems within the Community. As a result of the work performed in the “MyRailS” EURAMET project an AC calibration facility was developed consisting of a fictive power source was developed. This fictive power source can generate distorted sinusoidal voltages up to 25 kV-50 Hz and 15 kV-16.7 Hz as well as distorted sinusoidal currents up to 500 A with harmonic content up to 5 kHz or phase-fired current waveform stated in EN50463-2 standard. These waveforms are representative of those that appear during periods of acceleration and breaking of the train. Reference measuring systems have been designed and built consisting of high voltage and high current transducers adapted to multimeters, which function as digital recorders to acquire synchronized voltage and current signals. An approved procedure has been developed and an in-depth uncertainty analysis has been performed to achieve a set of uncertainty formulas considering the influence parameters. Different influence parameters have been analyzed to evaluate uncertainty contributions for each quantity to be measured: rms voltage, rms current, active power, apparent power and non-active power of distorted voltage and current waveforms. The resulting calculated global expanded uncertainty for the developed Energy Measuring Function calibration set up has been better than 0.5% for distorted waveforms. This paper is focused on presenting the complete set of expressions and formulas developed for the different influence parameters, necessary for uncertainty budget calculation of an Energy Measuring Function calibration.

1. Introduction

The opening of the railway infrastructures of the different countries belonging to the European Economic Community for the free movement of trains of different companies requires the on-board measurement of electrical energy in trains to achieve the European objectives stated in [1,2]. Therefore, regular calibration and verification of on-board Energy Measuring Function (EMF) for energy billing will be a necessity for years to come. These EMF must operate within the maximum permissible error limits under actual operating conditions according to the applicable standards [3]. Strong distortions in current and voltage waveforms are caused by speed control systems (thyristors, IGBTs, etc.) of locomotives when they increase or decrease their speed by accelerating or breaking their motors. Therefore, the calculation of the uncertainty of the reference EMS used for the calibration of EMF installed in the trains for energy billing is required. This calibration must be carried out in accordance with the requirements of EN 50463-2 [3] applying voltage and current waveforms with harmonics and distortions expected in railway electrical networks [4].
The MyRailS European Project [5] has not only made it possible to develop fictitious power generation facilities in the high voltage ranges or grid frequency and harmonic content of current and voltage waveforms that are representative of the service conditions, but also it has allowed to establish a robust calibration procedure with an uncertainty of active, apparent and non-active power measurements better than 0.5%. These special calibration facilities are now available and their calibration metrological capabilities (CMC) are evaluated and presented in this paper to demonstrate traceable measurements beyond the current EURAMET recognized CMCs for any metrological institute.
The reference EMS, as part of a European project [5], was introduced and presented in [6] and a calibration setup with its capabilities was fully described in [7]. In [8], the calibration setup described in [7] was compared with another alternative, and the corresponding estimation uncertainty results were also compared. However, detailed information related to the determination of uncertainty contributions, especially those related to signal processing, were not provided. The current paper defines different installation configurations for 16.7 Hz and 50 Hz grid frequencies from the calibration configuration described in [7] and emphasizes each formula developed ad hoc for the uncertainty contribution of electrical quantities involved in the calibration procedure. In addition, a practical case of “uncertainty budget” is also developed and presented. The procedure involves transferring the voltage and current waveforms from the high voltage circuit to measurable low voltage signals, as well as using two high-precision digital multimeters as recorders that synchronously sample distorted voltage and current signals. The high voltage and high current generators used for the calibration set up generate synchronously distorted voltage and current waveforms up to 25 kV at 50 Hz or 15 kV at 16.7 Hz and up to 500 A at 50 Hz or 16.7 Hz with different frequency content up to 5 kHz to reproduce actual or normative operation conditions.
Although many authors are concerned with the measurement of active, apparent and non-active power including reactive power (see [9,10,11,12]) even of distorted alternating currents and voltages (see [13,14,15,16,17,18,19,20,21,22,23,24]), there is no literature on the evaluation of uncertainty in power measurement with harmonics or with a phase-fired current waveform as required by standard [3]. The evaluation of uncertainty due to synchronous sampling of digital multimeters has been discussed in depth in [13,14], but no study has been found in the literature on the sensitivity of power measurements to certain processing parameters, such as the integration of the signals at a time other than the full period, nor the influence of the synchronization errors between the two DMMs on the power measurement, nor the method of signal integration.
This paper also presents novel model functions for active, apparent and non-active power measurements that represent the applied measuring procedure with its influence parameters. These model functions are fundamental to establish the “uncertainty budget” through the sensitivity coefficients of each one of the influence parameters. The formulas developed in this paper for the model functions and for the determining the uncertainty together with the calculated curves that relate the influence parameters with the measurement uncertainty are the main contributions of this paper. The uncertainty analysis according to the BIPM guide [22] is presented for the developed calibration set up to determine its best Calibration Metrological Capabilities. This work is intended to contribute to future reviews of on-board EMF calibration standards to establish internal guidelines for good calibration practices that can be used by accredited calibration laboratories.

2. Calibration Set-Up

A new calibration facility for EMS calibrations traced to national standards has been developed according to EN 50463-2 standard [3] through the European Project [5,6,7,8]. The voltage, current, and frequency values for the developed calibration set up were chosen to meet the standard requirements for laboratory calibrations or for on-site calibrations (see Figure 1). It consists of a phantom power generator (see Figure 2) composed of two independent synchronized sources: (1) a sinusoidal current source (50 Hz or 16.7 Hz) or a phase-fired current waveform up to 500 A rms (Figure 3) and (2) a high voltage source up to 25 kV, 50 Hz or 15 kV, 16.7 Hz. The current source is connected to the primary winding of a current injection transformer (4) and fed by a programmable calibration source (1′). The secondary of the current injection transformer consists of a current loop (3), composed by a high voltage cable (5) with its two cable terminations (6) short-circuited by means of a bus bar (7), where the EMS under calibration (8) is installed (see Figure 4 and Figure 5). The high voltage source (2) consists of two voltage transformers connected to the current loop, fed by a second programmable calibration source (2′). Depending on the operating conditions, the transformers are connected in parallel for 25 kV, 50 Hz to avoid excessive temperature rise in their windings, or in series for 15 kV-16.7 Hz, to avoid their magnetic saturation (see Figure 6).
This phantom generator works in the two different ways:
  • Mode C “Current with harmonics”: fictive electrical power is generated by injecting a sinusoidal current (50 Hz or 16.7 Hz) or a phase-fired current waveform up to 500 A rms with harmonic content up to 5 kHz in the current loop (see Figure 5). Simultaneously, a 50 Hz or 16.7 Hz high voltage is applied.
  • Mode V “Voltage with harmonics”: fictive electrical power is generated by applying a sinusoidal high voltage up to 25 kV, 50 Hz or 15 kV 16.7 Hz with harmonic content up to 5 kHz (see Figure 4) to the current loop. Simultaneously, a 50 Hz or 16.7 Hz current or a phase-fired current waveform up to 500 A is induced on the loop.
The two programmable voltage sources (1′) and (2′) are California models: type CSW5550 of 312 V, 16 A and type 5001-ix of 300 V, 16.6 A respectively. One source is connected to the low voltage winding of the high voltage transformer (9) (Arteche VEG-24 22,000/110 V, 50 Hz), while the other source is connected to primary winding of the current transformer (4) (Mimaven 160/4000 A, 50 Hz). Both current and voltage transformers were characterized in the frequency domain up to 5 kHz. This characterization is used to compensate the attenuation ratio at higher frequencies than the fundamental frequency: 50 Hz or 16.7 Hz. The frequency of the low voltage and current sources are synchronized to allow a controlled phase shift between both signals. This synchronization is performed by a device (SYNC) whose output is a trigger signal with the same frequency as the input signal. The input signal, taken as reference, can be chosen from mains network (see Figure 6a) or from another programmable or non- programmable source (see Figure 6b). The shift delay between both sources is an independent parameter that is set by a dedicated control software.

2.1. Reference Energy Measuring System

The reference energy measuring system performs synchronized measurements using independent voltage and current measuring systems. The current measuring system consists of a Fluxgate sensor type LEM ITN-900-S with a shunt resistor as the current transducer (10) and the high voltage measuring system consists of an improved R/C high voltage divider [25] used as a high voltage transducer (11). Two identical Keithley DMM7510 multimeters, working as recorders of 1 MS/s sampling rate, were set up to acquire both voltage and current waveforms. One multimeter works as the master, while the other is controlled by the first, working as a slave. The sampling rate and the starting time of slave DMM is the same as the master DMM and is configured by the control software.
Moreover, special software was developed to measure both voltage and current signals and to calculate the active, apparent, non-active power for sinusoidal waveforms and distorted waveforms (including harmonics and 90° phase-fired current waveform).
Table 1 shows the components of both voltage and current reference measuring systems used for the calibration set up. The frequency responses of the Scale Factor relative error for the high voltage divider developed by FFII-LCOE and for the Fluxgate current transducer, are shown in Figure 7 and Figure 8, respectively.

2.2. Characterization and Implementation

The frequency response of current and high voltage transformers is the ratio between the output of the transformer and the input signal supplied by its California source. These transformers used in the calibration set up causes a significant attenuation and phase displacement on the signal supplied by California sources. A frequency characterization of these current and voltage transformers makes it possible to compensate the attenuation and the phase shift influence. The transformer attenuation of each harmonic component was determined during these characterizations. With this purpose, the sinusoidal signal generated by each California source was varied between the main frequency (50 Hz or 16.7 Hz) and 5 kHz by steps of the main frequency value to determine both frequency responses. The limitation for the harmonic amplitude of the voltage transformer is 5% of the main component up to 50th harmonic for 50 Hz and up to 2% for 10th harmonic, while the harmonic limitation for current transformer is not less 10% up to 50th harmonic for both 50 Hz and 16.7 Hz and more than 10% up to 100th harmonic for 16.7 Hz. These data ensure the allowable tolerance on the rise time of leading edge 0.2 ms ± 0.1 ms of the phase-fire current waveform stated in the standard EN 50,463 because not more than 4% amplitude is required up to 21st harmonic, (see Figure 3). Therefore, different correction factors were applied for each harmonic component, resulting in higher values as a function of the harmonic frequency. To achieve the required phase-fired current signal in the current transformer secondary winding, the California source generated a distorted signal considering the characterization results. More detailed information on attenuation of the injector transformers and the form of compensation is provided in [8].
The required harmonic components are measured by the current and voltage measuring systems and regulated by means of California sources considering the deviation from the target harmonic amplitudes and phase shifts. The frequency response of the relative error of the Scale factor of the High Voltage divider and of the Fluxgate current transducer are shown in Figure 7 and Figure 8 respectively.

3. Uncertainty Analysis

An important aim of this calibration is to verify the error of the on-board trains energy measurement systems. The calibration uncertainty should provide lower figures than the maximum permissible errors for an accepted device. The acceptable limits for the energy measurement systems, EMS, of on-board trains limits are given in [3]. Therefore, next sections are addressed to determine the calibration uncertainty of EMS calibration set up according to [8]. The definitions for electric power quantities under non-sinusoidal or unbalanced conditions are used according to [23] and studied in [24]. Appropriate signal processing treatment of the acquired AC voltage and current signals is followed according to [15,16,17,18,19,20] by synchronous digital synthesis and sampling.

3.1. Model Functions

3.1.1. Model Functions for Rms Voltage and Current, Active, Apparent and Non-Active Power Obtained by Digital Sampling at the Low Voltage Side

A low voltage signal v(t) is sampled by means of a digital multimeter (DMM) configured for direct voltage measurements (DCV). As the voltage is sampled n times per period, the trapezoidal rule can be used to calculate the root mean square of the voltage Vrms(tj) at the low voltage side for each specific period, starting at the instant tj:
V r m s ( t j ) = 1 T t j t j + T v ( t ) 2 d τ 1 n k = 1 n 1 v k 2 + v n 2 + v 0 2 2
t j = t 0 , t 0 + T , t 0 + 2 T , , t 0 + ( N c 1 ) T
where:
  • T = 1/f: period of the power system.
  • f: fundamental frequency of the power system.
  • t0: time origin of the starting measuring times.
  • v0: acquired voltage in the first sampling interval of the period under consideration.
  • vk: acquired voltage in the kth sampling interval of the period under consideration.
  • vn: acquired voltage in the last sampling interval of the period under consideration.
  • n: Number of sample intervals per period. n = T / h s
  • hs = 1/fs: Sampling interval.
  • fs: Sampling rate.
  • Nc: Number of periods included in the measuring time interval used to determine the energy.
Considering that the current i(t) measured at the low voltage side causes a voltage drop ν’(t) through a shunt resistance Rs a similar approach can be applied to calculate the root mean square of the current Irms(tj):
I r m s ( t j ) = 1 T t j t j + T i ( t ) 2 d τ 1 n k = 1 n 1 i k 2 + i n 2 + i 0 2 2 = 1 R s 1 n k = 1 n 1 v k 2 + v n 2 + v 0 2 2
where:
  • ik: acquired current in the kth sampling interval of the period under consideration.
  • i0: acquired current in the first sampling interval of the period under consideration.
  • in: acquired current in the last sampling interval of the period under consideration.
The active power measured at the low voltage side, P l v ( t j ) , is calculated by the following formula:
P l v ( t j ) = 1 T t j t j + T v ( t ) i ( t ) d t 1 n k = 1 n 1 v k i k + v n i n + v 0 i 0 2
The apparent power measured at the low voltage side, S l v ( t j ) , is given by the formula:
S l v ( t j ) = V r m s ( t j ) I r m s ( t j )
and the non-active power, N l v ( t j ) , measured at the low voltage side by the following expression:
N l v ( t j ) = S l v 2 ( t j ) P l v 2 ( t j )

3.1.2. Model Functions of Voltage and Current Sampled Values at the High Voltage Side

The model function corresponding to the high voltage sample acquired at the kth sampling interval, V k , depends on the high voltage divider and DMM.
V k = v k 1 + δ v 1 + δ v 2 V F S v k S F V D 1 + j δ j , V D 1 + j c j δ j , t
where:
  • v k : sampled voltage value acquired at the low voltage side.
  • δ v 1 : constant term of error of the DMM, for DCV measurements, typically a percentage of the voltage reading.
  • δ v 2 : additional term of error of the DMM, affected by the ratio between the full-scale voltage of the DMM, V F S , and the absolute value of the sampled voltage v k .
  • S F V D : calibrated scale factor of the high voltage divider.
  • δ j , V D : corrections of the voltage divider scale factor, such us drift ( δ 1 , V D ), temperature coefficient ( δ 2 , V D ) , short term stability ( δ 3 , V D ) and non-linearity with voltage ( δ 4 , V D ).
  • c j : sensitivity coefficient of the voltage phase displacement.
  • δ j , t : corrections of the sampled voltage value due to phase displacements. The phase shift of the DMM (negligible), the high voltage divider phase calibration ( δ 1 , t ), the high voltage divider phase drift (negligible), the high voltage divider residual phase correction ( δ 2 , t ).
The model function corresponding to the high current measurement acquired at the kth sampling interval depends on the current transducer (Fluxgate sensor), shunt and DMM.
I k = v k ´ 1 + δ v 1 + δ v 2 U F S v k ´ S F C T 1 + j δ j , C T 1 R s 1 + j δ j , R s 1 + j c j δ j , t
where:
  • v k ´ : sampled voltage value measured at the shunt.
  • δ v 1 : constant term of error of the DMM connected to the shunt, for DCV measurements.
  • δ v 2 : additional term of error of the DMM connected to the shunt, affected by the ratio between the full-scale voltage of the DMM, V F S , and the absolute value of the sampled voltage v k ´ .
  • S F C T : calibrated scale factor of the current transducer.
  • δ j , C T : corrections of the current transducer scale factor, such us drift ( δ 1 , C T ), temperature coefficient ( δ 2 , C T ) and non-linearity with current ( δ 3 , C T ) (see note below).
  • R s : calibrated value of the shunt resistor.
  • δ j , R s : corrections of the calibrated shunt value, such us: drift ( δ 1 , R s ), temperature coefficient ( δ 2 , R s ), and variation versus frequency ( δ 3 , R s ).
  • c j : sensitivity coefficient of the voltage phase displacement in the current circuit.
  • δ j , t : corrections of the sampled voltage value due to phase displacements. The phase shift of the DMM (negligible), the current sensor phase calibration ( δ 3 , t ), the current sensor phase drift ( δ 4 , t ), the current sensor residual phase correction ( δ 5 , t ) and the shunt phase error ( δ 6 , t ).
Note: The effect of leakage flux from the current transformer on the scale factor of the current sensor is negligible (lower than 0.0025%).
Amplitude corrections in the scale factor for the harmonic components (Figure 7 and Figure 8) are not considered in Formulas (6) and (7) because amplitude of harmonic components bigger than 5th is not more than 10% of main component and its scale factor correction due to frequency dependence is less than 0.1% (see Figure 7 and Figure 8). Consequently, these uncertainty contributions can be neglected. For 3rd harmonic component the scale factor is very close to the rated scale factor for the main frequency and no correction is needed.

3.1.3. Model Functions for the Quantities at the High Voltage Side

The function model for the square of the rms voltage value calculated at the high voltage side is defined by the following formula:
V R M S ( t j ) = 1 n k = 1 n 1 V k 2 + V n 2 + V 0 2 2 1 + l S p l
where:
  • 1 n k = 1 n 1 V k 2 + V n 2 + V 0 2 2 is the rms voltage value calculated by trapezoidal rule from a finite number of samples Vk (k = 1: n) taken from the recorded voltage signal.
Spl: are the corrections of the rms voltage value due to signal processing restrictions, such as finite sampling rate ( S p 1 ) and integration time different of the period ( S p 2 ). The uncertainty contribution of these terms is described in Sections Uncertainty Contribution Due to Trapezoidal Integration Rule ( S p 1 ) and Uncertainty Contribution Due to an Integration Time Different to the Complete Period ( S p 2 ), respectively.
Similarly, the function model for the current rms value calculated at the high voltage side will be:
I R M S ( t j ) = 1 n k = 1 n 1 I k 2 + I n 2 + I 0 2 2 1 + l S p l
where:
  • 1 n k = 1 n 1 I k 2 + I n 2 + I 0 2 2 is the current rms value calculated by trapezoidal rule from a finite number of samples Ik (k = 1:n) taken from the recorded voltage signal.
S p l : are the corrections of the rms current value due to signal processing restrictions, such as finite sampling rate ( S p 1 ) and integration time different of the period ( S p 2 ). The uncertainty contribution of these terms is described in Sections Uncertainty Contribution Due to Trapezoidal Integration Rule ( S p 1 ) and Uncertainty Contribution Due to an Integration Time Different to the Complete Period ( S p 2 ), respectively.
The apparent power at the high voltage side, SHV(tj), is directly determined from the rms voltage and current values determined by (8) and (9):
S H V ( t j ) = V R M S ( t j ) I R M S ( t j )
The function model for the active power measured at the high voltage side by acquisition of voltage and current samples is given by the following expression:
P H V t j = 1 n k = 1 n 1 V k I k + V n I n + V 0 I 0 2 · 1 + l S p l + l S p l + S p 3
where:
  • Sp3: correction of the active power value due to synchronization error between both DMM.
And the non-active power, NHV(tj), will be given by:
N H V ( t j ) = S H V 2 ( t j ) P H V 2 ( t j )

3.2. Uncertainty Analysis

3.2.1. Uncertainty of the Calculated Rms Voltage and Current Values at the High Voltage Side

Considering Formula (8) as the function model of rms voltage value, considering the auxiliary function y(tj) and assuming as the best estimations of S p 1 and S p 2 equal to zero:
y ( t j ) = 1 n k = 1 n 1 V k 2 + V n 2 + V 0 2 2 V R M S 2 ( t j )
The V R M S ( t j ) uncertainty can be expressed by:
u [ V R M S ( t j ) ] = V R M S ( t j ) y ( t j ) 2 u 2 [ y ( t j ) ] + V R M S 2 ( t j ) l u 2 ( S p l )
where:
V R M S ( t j ) y ( t j ) = 1 2 V R M S ( t j )
Taking into account the expression given in (6) as the function model of each voltage sample, V k , measured at the high voltage side in the sampling interval kth and assuming a fully correlation between voltage measurements (correlation coefficient = 1), each calculated rms voltage uncertainty at the high voltage side, VRMS(tj), can be expressed in per unit value (p.u.):
u p u [ V R M S ( t j ) ] = u 2 ( δ v 1 ) + V A R V ( t j ) V R M S ( t j ) V F S S F V D V R M S ( t j ) 2 u 2 ( δ v 2 ) + u 2 ( S F V D ) + l u 2 ( δ l , V D ) + l c l 2 u 2 ( δ l , t ) + l u 2 ( S p l )
where the average rectified high voltage value is:
V A R V ( t j ) = S F V D 1 n k = 1 n 1 v k + v n + v 0 2
  • u ( δ v 1 ) : represents the uncertainty in p.u. of the correction term δv1.
  • u ( δ v 2 ) : represents the uncertainty in p.u. of the correction term δv2.
  • u S F V D : represents the calibration uncertainty in p.u. of the high voltage divider scale factor,   S F V D .
  • u ( δ l , V D ) : represents the uncertainty in p.u. factor δj,VD.
  • u ( δ l , t ) : represents the uncertainty in p.u. of the correction terms related to the phase displacements due to the digitizer and the divider.
  • c l : represents the sensitivity coefficient of each lth correction term.
  • u ( S p l ) : represents the uncertainty in p.u. of the correction terms related to the signal processing restrictions related to a rms of the voltage in a period T starting at tj.
Considering a sinusoidal voltage waveform with a low harmonic content (THD < 2%):
V A R V ( t j ) V R M S ( t j ) = 1 F F v ( t j ) 2 2 π 0.9
The final formula of the uncertainty expressed in p.u. using the form factor of the voltage signal in the period T starting at the tj instant, FFv(tj), is the following, where FFv(tj) ≈ 1.11.
u p u [ V R M S ( t j ) ] = u 2 ( δ v 1 ) + V F S 2 S F V D 2 1.11 2 V R M S 2 ( t j ) u 2 ( δ v 2 ) + u 2 ( S F V D ) + l u 2 ( δ l , V D ) + l c l 2 u 2 ( δ l , t ) + l u 2 ( S p l )
Applying the same considerations and assumptions as the ones followed for the uncertainty of a calculated rms voltage value, the uncertainty of a calculated rms current value can be expressed by the formula:
u p u [ I R M S ( t j ) ] = u 2 ( δ v 1 ) + I F S 2 S F C T 2 F F i 2 ( t j ) I R M S 2 ( t j ) u 2 ( δ v 2 ) + u 2 ( S F C T ) + l u 2 ( δ l , C T ) + u 2 ( R s ) + l u 2 ( δ l , R s ) + l c l , 2 u 2 ( δ l , t ) + l u 2 ( S p l )
where:
I A R V t j = S F C T I a r v t j ;   I R M S t j = S F C T I r m s t j ;   I F S = V F S R s
  • u S F C T : represents the calibration uncertainty in p.u. of the current transducer scale factor,   S F C T .
  • u ( δ l , C T ) : represents the uncertainty of the correction term related to the current transducer δl,CT.
  • u R s : represents the calibration uncertainty in p.u. of the current shunt,   R s .
  • u ( δ l , t ) : represents the uncertainty of the correction terms related to the phase displacements due to current transducer, shunt resistance, digitizer.
  • c l : represents the sensitivity coefficient of each lth correction term.
  • u ( S p l ) : represents the uncertainty in p.u. of the correction terms related to the signal processing restrictions related to a rms of the current in a period T starting at tj.
When injecting a sinusoidal current waveform with a low harmonic content (THDr < 2%) the form factor to be applied in the expression (20) will is F F i t j 1.11 .
For a phase-fired current waveform of angle α, as shown in Figure 3, the quotient between the rectified average current and its rms value in a period follows the expression:
I a r v ( t j ) I r m s ( t j ) = 2 π 1 + cos ( α ) ( 1 α π ) + 1 2 π sin ( 2 α ) w h e r e : I a r v ( t j ) = I p k π 1 + cos ( α ) I r m s ( t j ) = I p k 2 ( 1 α π ) + 1 2 π sin ( 2 α )
For example, for the case of a α = 90° phase-fired waveform the form factor is:
F F i t j = I r m s t j I a r v t j 1.56

3.2.2. Uncertainty of the Calculated Power Quantities at the High Voltage Side

Active Power

Considering Formula (11) as the function model of the active power and assuming as the best estimations of S p 3 equal to zero:
u [ P H V ( t j ) ] = k P H V ( t j ) V k 2 u 2 ( V k ) + k P H V ( t j ) I k 2 u 2 ( I k ) + P H V 2 ( t j ) l u 2 ( S p l ) + P H V 2 ( t j ) l u 2 ( S p l ) + P H V 2 ( t j ) u 2 ( S p 3 )
To simplify editing, this expression is divided in four uncertainty terms: the first term is related to independent samples at low voltage side of both signals: voltage and current. The second and third terms are related to the high voltage and high current scale factors, respectively, and the fourth term to the signal processing:
u p u [ P H V ( t j ) ] = u 1 2 + u 2 2 + u 3 2 + u 4 2
Assuming a correlation coefficient among sampled measurements of +1 the first uncertainty term can be expressed in per unit value (p.u.) as:
u 1 = u 2 ( δ v 1 ) + u 2 ( δ v 1 ) + V F S 2 I a r v 2 P l v 2 ( t j ) u 2 ( δ v 2 ) + I F S 2 V a r v 2 P l v 2 ( t j ) u 2 ( δ v 2 )
This expression can be particularized for sinusoidal current with low harmonic content and for 90° phase-fired current waveform.
The active power at the low voltage side for a sinusoidal voltage or current with low harmonic content can be calculated by the formula:
P l v ( t j ) I r m s _ 1 ( t j ) V r m s _ 1 ( t j ) c o s φ
and the u1 uncertainty term can be expressed by the formula:
u 1 = u 2 ( δ v 1 ) + u 2 ( δ v 1 ) + 1 F F i 2 ( t j ) c o s 2 φ V F S 2 V r m s 2 ( t j ) u 2 ( δ v 2 ) + 1 F F v 2 ( t j ) c o s 2 φ I F S 2 I r m s 2 ( t j ) u 2 ( δ v 2 )
For the case of α phase-fired current waveform and an undistorted sinusoidal voltage with angle phase shift of ϕ = 0° between voltage and current, the active power can be calculated by:
P l v ( t j ) = ( 1 α π ) + 1 2 π s i n ( 2 α ) I r m s ( t j ) V r m s ( t j )
For α = 90°:
P l v ( t j ) 0.707 I r m s ( t j ) V r m s ( t j )
and using the form factor for voltage and current 1.11 and 1.56 respectively, the u1 uncertainty term will be:
u 1 = u 2 ( δ v 1 ) + u 2 ( δ v 1 ) + 1 1.5 6 2 × 0.70 7 2 V F S 2 V r m s 2 ( t j ) u 2 ( δ v 2 ) + 1 1.1 1 2 × 0.70 7 2 I F S 2 I r m s 2 ( t j ) u 2 ( δ v 2 )
resulting:
u 1 = u 2 ( δ v 1 ) + u 2 ( δ v 1 ) + 0.82 V F S 2 V r m s 2 ( t j ) u 2 ( δ v 2 ) + 1.62 I F S 2 I r m s 2 ( t j ) u 2 ( δ v 2 )
To determine the active power uncertainty calculated at the high voltage side, PHV(tj), the rest of the uncertainty terms u2, u3 and u4 should be considered:
u 2 = u 2 ( S F V D ) + l u 2 ( δ l , V D ) + l c l 2 u 2 ( δ l , t )
u 3 = u 2 ( S F C T ) + l u 2 ( δ l , C T ) + u 2 ( R s ) + l u 2 ( δ l , R s ) + l c l , 2 u 2 ( δ l , t )
u 4 = l u 2 S p l + l u 2 S p l + u 2 S p 3 2

Apparent Power

Considering Formula (10), the apparent power uncertainty can be calculated as:
u S H V t j = S H V t j V R M S 2 u 2 V R M S t j + S H V t j I R M S 2 u 2 I R M S t j
and using per unit values this formula is transformed to:
u p u [ S H V ( t j ) ] = u p u 2 [ V R M S ( t j ) ] + u p u 2 [ I R M S ( t j ) ]

Non-active power

Considering Formula (12), the non-active power uncertainty is:
u [ N H V ( t j ) ] = 1 N H V ( t j ) S H V 2 ( t j ) u 2 [ S H V ( t j ) ] + P H V 2 ( t j ) u 2 [ P H V ( t j ) ]
and using per unit values this formula is transformed to:
u p u [ N H V ( t j ) ] = S H V ( t j ) N H V ( t j ) 4 u p u 2 [ S H V ( t j ) ] + P H V ( t j ) N H V ( t j ) 4 u p u 2 [ P H V ( t j ) ]

3.2.3. Uncertainty of the Signal Processing Contributions

Uncertainty Contribution Due to Trapezoidal Integration Rule ( S p 1 )

To achieve a null error due to numerical integration applying the trapezoidal integration rule when sinusoidal waveforms are used the Nyquist criterion must be satisfied up to and including the maximum harmonic hmax, n ≥ 2 · hmax + 1. For example, if hmax = 100, then n ≥ 201 samples are required. This assumes that, for the fundamental period of 20 ms, 201 sampling intervals (202 samples) are required. Therefore, the sampling rate must be at least fs,min slightly higher than the minimum sampling rate according to the Nyquist criterion:
f N y q u i s t = 2 h m a x 20   ms = 2 100 20   = 10   kHz f s , m i n = 202 20   ms = 10.1 kHz
However, the relative error due to the application of the trapezoidal rule when the rms value of a 90° phase-fired waveform is much bigger because of their large frequency content. This relative error can be determined as a function of the number of samples per period, n. For example, for a sampling rate of 500 kHz that means n = 10,000 when the period is 20 ms (50 Hz), the relative error is less than 0.02% (see Figure 9a).
The relative error curve for the active power measurement is twice (0.04%) when the current signal is a 90° phase-fired waveform using a sinusoidal voltage (see Figure 9b).

Uncertainty Contribution Due to an Integration Time Different to the Complete Period ( S p 2 )

If the numeric integration applying the trapezoidal rule is performed up to an integer n1 of sampling intervals different from the integer n of sampling intervals whose duration is the same as the entire period, a relative error is caused when rms voltage, rms current, active power, apparent power or non-active power are measured.
For sinusoidal voltage and current signals assuming n >> 1, the rms value relative error is given by the following analytical formula:
ε ( λ 1 ) c o s 2 2 π f t j + λ π s i n ( 2 π λ ) 2 π
where λ = n1/n.
The curve family of the relative error, ε, expressed in%. for different starting times tj is shown in Figure 10a. The starting time tj should be close to zero, while λ should be as close as possible to 1 as much as possible to achieve a negligible error. For instance, for tj = 0.01·T and n = 10,000 (fs = 500 kHz for f = 50 Hz) the relative error is less than 5 × 10−5 % (ε < 5 × 10−5 %). When the difference between n1 and n is known and tj ≤ 0.01·T, the relative error can be determined using the curve shown in Figure 10b. For example, if n1n = 10 samples (λ = 0.999 for Figure 10a) then ε < 5 × 10−5 %.
The following analytical formula for the relative error of the measured active power of sinusoidal voltage and current signals with a ϕ phase shift can be determined, assuming n >> 1:
ε = λ 1 c o s 2 ( 2 π f t j + λ π ) φ s i n ( 2 π λ ) 2 π
For ϕ = 0 this formula is transformed to (39) and the same conclusions as the ones achieved for rms voltage and current measurements can be extended for the active power.
In the case of a 90° phase-fired waveform current signal using a sinusoidal voltage, same relative error ε < 5 × 10−5 % is obtained by numerical methods, for n1n = 10 samples.

Uncertainty Due to the Synchronization Error between Multimeters ( S p 3 )

Considering that the maximum synchronization error between the digital multimeters is ±1 μs, the relative error of the active power expressed in (%) due to this synchronization error of ±1 μs for the cosφ range between from 0.85 and 1 is calculated. Figure 11a shows the error trend for 90° phase-fired current waveform and Figure 11b for sinusoidal current waveform with a low harmonic content (THD < 2%) and in both cases for a sinusoidal voltage signal with harmonic content of 10% for the 5th harmonic and 3% for the 11th harmonic.

4. Application of the Uncertainty Estimation of the EMS Calibration Set Up

Applying Formulas (19) and (20), the rms voltage and current uncertainties can be determined. Based on these results, the apparent power uncertainty can be determined by applying Formula (36). The active power uncertainty is given by Formula (24) and the non-active power uncertainty is calculated by applying Formula (38). Table 2 shows a summary of all calculated expanded uncertainties for the EMS calibration set up developed by LCOE, for 90° phase-fired current waveform. For this waveform, a sensitivity coefficient of −6.26 × 10−5%/μrad was obtained from the relative error of −0.01966% obtained for 1 μs shift between voltage and current signals of a 90° phase-fired current waveform, corresponding to cosφ = 1, as shown in Figure 11a. Therefore, taking into account that 1 μs corresponds to 314.16 μrad for 50 Hz signal, the sensitivity coefficient is equal to −0.01966%/314.16 = −6.26 × 10−5%/μrad.
The uncertainty estimation of the EMS calibration set up developed by LCOE is analyzed in detail for active power measurements by applying Formulas (31)–(34) (see Table 3). The uncertainty contributions due to independent low voltage sampling of both voltage and current signals are given by Formula (31), in which the sensitivity coefficients corresponding to the constant errors for voltage and current measurements, δv1 and δ’v1, are 1 and the ones for the digitizer variable errors, δv2 and δ’v2, depend on the voltage or current rms value to be measured by the DMMs, which are assumed to be higher than the 10% of the full scale in both cases. The uncertainty contributions due to the relation between the high voltage and current scale factors are given by Formulas (32) and (33), respectively. The standard uncertainties of the different influence parameters such as drift, temperature coefficient, short term stability, non-linearity and frequency dependence, for voltage divider, current transducer and shunt, are collected from the manufacturer data sheet. A rectangular probability distribution is assumed in this case. Influence parameter such as the voltage divider, current transducer and shunt calibration uncertainties are collected from the calibration certificates, assuming a normal probably distribution. The standard uncertainties related to signal processing method are obtained by applying Formula (34) from the analysis carried out in Section 3.2.2. The standard uncertainty corresponding to the trapezoidal integration method, Sp1 and Sp’1, depends on the sampling rate. For example, assuming a rectangular probability distribution, for fs = 500 kHz it is 0.04% for active power (Figure 9). This total error contribution is split in two equal contributions for both the voltage and the current, equal to √2 × 0.04%/2 = 0.028% (Sp1 = 0.028% and Sp2 = 0.028%). The standard uncertainty due to integrating in an integer n1 of sampling intervals different from the integer n of sampling intervals whose duration is the same as the entire period, Sp2 and Sp’2, is determined by applying Formula (40). Assuming a rectangular distribution, if the difference between the acquired number of samples and the number of samples corresponding to a complete period is 10 samples, when n = 10,000 samples and the sampling rate, fs, is 500 kSamples/s, the relative error is less than 5 × 10−5. Finally, the standard uncertainty due to the synchronization between both digital multimeters, Sp3, is derived of the manufacturer data sheet. Considering the maximum synchronization time given by the manufacturer of 1 μs, and assuming a rectangular probability distribution, this standard uncertainty results in a value of 314μrad/√3. From the uncertainty contributions referred in Table 3 and applying Formula (24), the active power expanded uncertainty is lower than 0.3%.

5. Conclusions

A new EMS calibration set up is available for AC voltages up to 25 kV addressed to on board energy metering under distorted conditions. It consists of a phantom power generator composed of two independent synchronized sources capable of generating high voltages with superimposed harmonics components up to 5 kHz, and a phase-fired current waveform or a sinusoidal current up to 500 A with harmonics up to 5 kHz, according to the EN 50463-2 requirements.
Based on two synchronized digital multimeters, to sample the voltage and current waveforms to be measured, two precise traceable measuring systems, one for high voltage up to 15 kV (16.7 Hz) or 25 kV (50 Hz) and the other for high current up to 500 A, with a bandwidth up to 5 kHz, were developed.
Function models were introduced for rms voltage, rms current and power quantities (active, apparent and non-active) and analytical formulas have been derived to analyze the uncertainty contribution due to different influence parameters. Particular attention was paid to uncertainty contributions of signal processing, such as integration by applying the trapezoidal rule, the integral at a set of number of sampling intervals different to the full period duration, and synchronization error between digital multimeters. This uncertainty analysis provides an expanded uncertainty for the EMS calibrations for active power measurement and apparent power better than 0.3%, and for non-active power better than 0.5% according to the target of the “MyRailS” European project [5,6,7,8], solving the traceability lack for this type of measurements.

Author Contributions

Conceptualization and methodology of fictive power sources P.S. and F.G.; software and validation A.K. and J.R.; implementations A.K. and J.R.; analysis and research A.K. and F.G.; writing F.G., P.S. and A.R.; supervision F.G., P.S. and A.K.; uncertainty formulation A.K. All authors have read and agreed to the published version of the manuscript.

Funding

The latter received funding from the EMPIR program co-financed by the Participating States and from the European Union’s Horizon 2020 research and innovation program. The EMPIR program belongs to EURAMET association.

Acknowledgments

The research here described has been developed in the framework of the 16ENG04 MyRailS EMPIR Project.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Commission Regulation (EU) No 1301/2014 of 18 November 2014 on the Technical Specifications for Interoperability Relating to the ‘Energy’ Subsystem of the Rail System in the Union; EU: Brussels, Belgium, 2014.
  2. Commission Regulation (EU) No 1302/2014 of 18 November 2014 Concerning a Technical Specification for Interoperability Relating to the ‘Rolling Stock—Locomotives and Passenger Rolling Stock’ Subsystem of the Rail System in the European Union; EU: Brussels, Belgium, 2014.
  3. EN 50463-2. Railway Applications—Energy Measurement on Board Trains—Part 2: Energy Measuring. Available online: https://standards.globalspec.com/std/10259948/EN%2050463-2 (accessed on 1 October 2017).
  4. EN 50163. Railway Applications—Supply Voltages of Traction Systems. Available online: https://standards.iteh.ai/catalog/standards/clc/de298b48-41ab-46b9-9328-99740496e4a2/en-50163-2004-a2-2020 (accessed on 7 February 2020).
  5. 16ENG04 MyRailS. Metrology for Smart Energy Management in Electrical Railway Systems; EURAMET H2020 Project; EMPIR Program, 2017–2021. Available online: https://www.euramet.org/research-innovation/search-research-projects/details/?tx_eurametctcp_project[project]=1484&tx_eurametctcp_project[controller]=Project&tx_eurametctcp_project[action]=show (accessed on 25 November 2021).
  6. Garnacho, F.; Rovira, J.; Khamlichi, A.; Simón, P.; Garrido, I. Reference Energy Measuring System for On-Board Calibration of EMS Installed in Locomotives; SEIA: Washington, DC, USA, 2020. [Google Scholar]
  7. Garnacho, F.; Rovira, J.; Khamlichi, A.; Simón, P.; García, T.; Istrate, D. Calibration set-up for energy measuring systems installed in AC railway systems. In Proceedings of the IEEE Vehicle Power and Propulsion Conference (VPPC), Gijon, Spain, 18 November–16 December 2020. [Google Scholar]
  8. Istrate, D.; Khamlichi, A.; Soccalingame, S.; Rovira, J.; Fortune, D.; Sira, M.; Simon, P.; Garnacho, F. Laboratory calibration of energy measurement systems (EMS) under AC distorted waveforms. Sensors 2020, 20, 6301. [Google Scholar] [CrossRef] [PubMed]
  9. Kaczmarek, M.; Stano, E. Measuring system for testing the transformation accuracy of harmonics of distorted voltage by medium voltage instrument transformers. Measurement 2021, 181, 109628. [Google Scholar] [CrossRef]
  10. Serov, A.N.; Novitskiy, A.; Shatokhin, A.A.; Schlegel, S.; Dolgacheva, E.A.; Westermann, D. The Influence of Power Frequency Deviation on the Active and Reactive Power Measurement Error with the Application of DFT. In Proceedings of the IEEE 2019 20th International Symposium on Power Electronics (Ee), Novi Sad, Serbia, 23–26 October 2019. [Google Scholar]
  11. Filipović-Grčić, D.; Filipović-Grčić, B.; Krajtnerc, D. Frequency response and harmonic distortion testing of inductive voltage transformer used for power quality measurements. Procedia Eng. 2017, 202, 159–167. [Google Scholar] [CrossRef]
  12. Vlahinić, S.; Brnobić, D.; Vučetić, D. Measurement and analysis of harmonic distortion in power distribution systems. Electr. Power Syst. Res. 2009, 79, 1121–1126. [Google Scholar] [CrossRef]
  13. Kürten Ihlenfeld, W.G. Maintenance and Traceability of AC Voltages by Synchronous Digital Synthesis and Sampling; PTB Report E-75; Physikalisch-Technische Bundesanstalt: Braunschweig, Germany, 2001. [Google Scholar]
  14. Kürten Ihlenfeld, W.G.; Mohns, E.; Bachmair, H.; Ramm, G.; Moser, H. Evaluation of the Synchronous Generation and Sampling Technique. IEEE Trans. Instrum. Meas. 2003, 52, 371–374. [Google Scholar] [CrossRef]
  15. Krajewski, M. Constructing an uncertainty budget for voltage RMS measurement with a sampling voltmeter. Metrologia 2018, 55, 95–105. [Google Scholar] [CrossRef]
  16. L’vova, E.V.; Glukhova, R.M.; Vagarina, N.S.; Umnova, E.G.; Melnikova, N.I.; Sytnik, A.A. Comparison of high-accuracy measurement methods for power quality characterization. In Proceedings of the 2017 IEEE Conference of Russian Young Researchers in Electrical and Electronic Engineering, St. Petersburg and Moscow, Russia, 1–3 February 2017; pp. 1558–1562. [Google Scholar]
  17. Leniček, I. Realization and traceability of AC power standard at frequency of 50 Hz. In Proceedings of the 2018 First International Colloquium on Smart Grid Metrology, Split, Croatia, 24–27 April 2018; pp. 1–5. [Google Scholar]
  18. Guilherme, W.; Ihlenfeld, K.; Ramm, G.; Bachmair, H.; Moser, H. Evaluation of the synchronous generation and synchronous sampling technique for the determination of low frequency AC quantities. In Proceedings of the Conference Digest Conference on Precision Electromagnetic Measurements, Ottawa, ON, Canada, 16–21 June 2002; pp. 398–399. [Google Scholar]
  19. Petrovic, P.; Stevanovic, M. Measuring active power of synchronously sampled AC signals in presence of interharmonics and subharmonics. IEE Proc. Electr. Power Appl. 2006, 153, 227–235. [Google Scholar] [CrossRef]
  20. Chen, K.-L.; Chen, S.-F.; Ho, T.-H. New AC Power Measurement Standard Based on Digital Sampling at CMS. In Proceedings of the 2020 Conference on Precision Electromagnetic Measurements (CPEM), Denver, CO, USA, 24–28 August 2020; pp. 1–2. [Google Scholar]
  21. Iuzzolino, R.; Ihlenfeld, W.G.K. High-Accuracy Methods and Measurement Procedures for Power Quality Parameters Using the Digital Synchronous Sampling Technique. IEEE Trans. Instrum. Meas. 2007, 56, 426–430. [Google Scholar] [CrossRef]
  22. JCGM 100, Evaluation of Measurement Data—Guide to the Expression of Uncertainty in Measurement; Bureau International des Poids et Measures: Paris, France, 2008.
  23. 1459-2010—IEEE Standard Definitions for the Measurement of Electric Power Quantities under Sinusoidal, Nonsinusoidal, Balanced, or Unbalanced Conditions; IEEE: Piscataway, NJ, USA, 2010.
  24. Digital Methods in Waveform Metrology; NBS Special Publications 707; National Bureau of Standards: Gaithersburg, MD, USA, 1985.
  25. Garnacho, F.; Khamlichi, A.; Rovira, J. The Design and Characterization of a Prototype Wideband Voltage Sensor Based on a Resistive Divider. Sensors 2017, 17, 2657. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. EMS calibrations according to EN 50463-2 standard: (a) Calibration AC facility for laboratory; (b) On site Calibration AC facility for on board.
Figure 1. EMS calibrations according to EN 50463-2 standard: (a) Calibration AC facility for laboratory; (b) On site Calibration AC facility for on board.
Sensors 21 07967 g001
Figure 2. Conceptual Generation set up for 50 Hz and 16.7 Hz developed by FFII-LCOE.
Figure 2. Conceptual Generation set up for 50 Hz and 16.7 Hz developed by FFII-LCOE.
Sensors 21 07967 g002
Figure 3. 90° Phase-fired current waveform according to EN 50463: (a) Waveform in the time domain; (b) Harmonic content versus harmonic number.
Figure 3. 90° Phase-fired current waveform according to EN 50463: (a) Waveform in the time domain; (b) Harmonic content versus harmonic number.
Sensors 21 07967 g003
Figure 4. (a) Insulated current loop composed of a high voltage cable short-circuited by a bus bar where the EMS under calibration is installed; (b) High voltage source composed of two voltage transformers connected in parallel (25 kV; 50 Hz) or in series (15 kV; 16.7 Hz).
Figure 4. (a) Insulated current loop composed of a high voltage cable short-circuited by a bus bar where the EMS under calibration is installed; (b) High voltage source composed of two voltage transformers connected in parallel (25 kV; 50 Hz) or in series (15 kV; 16.7 Hz).
Sensors 21 07967 g004
Figure 5. High current (up to 500 A) with harmonics injection to the current loop: (a) using the mains network when 50 Hz is required, (b) using a non-programable source when 16.7 Hz are required. Simultaneously, a sinusoidal high voltage can be applied.
Figure 5. High current (up to 500 A) with harmonics injection to the current loop: (a) using the mains network when 50 Hz is required, (b) using a non-programable source when 16.7 Hz are required. Simultaneously, a sinusoidal high voltage can be applied.
Sensors 21 07967 g005
Figure 6. Generation of High voltage with harmonics application to the current loop, with a simultaneous injection of a sinusoidal or a phase-fired current up to 500 A: (a) Set up for 25 kV-50 Hz; (b) Set up for 15 kV-16.7 Hz.
Figure 6. Generation of High voltage with harmonics application to the current loop, with a simultaneous injection of a sinusoidal or a phase-fired current up to 500 A: (a) Set up for 25 kV-50 Hz; (b) Set up for 15 kV-16.7 Hz.
Sensors 21 07967 g006
Figure 7. Frequency response of the Scale Factor relative error for the improved high voltage divider developed by FFII-LCOE.
Figure 7. Frequency response of the Scale Factor relative error for the improved high voltage divider developed by FFII-LCOE.
Sensors 21 07967 g007
Figure 8. Frequency response of the Scale Factor relative error for the Fluxgate current transducer: LEM ITN-900-S.
Figure 8. Frequency response of the Scale Factor relative error for the Fluxgate current transducer: LEM ITN-900-S.
Sensors 21 07967 g008
Figure 9. Relative error for the 90° phase-fired current waveform (50 Hz) vs. sampling rate: (a) of the rms and (b) of the active power.
Figure 9. Relative error for the 90° phase-fired current waveform (50 Hz) vs. sampling rate: (a) of the rms and (b) of the active power.
Sensors 21 07967 g009
Figure 10. Relative error ε of the rms value of a digitized sinusoidal waveform due to integrating in an integer n1 of sampling intervals different from the integer n of sampling intervals whose duration is the same as the entire period: (a) ε(%) vs. λ; (b) ε(%) vs. n1n.
Figure 10. Relative error ε of the rms value of a digitized sinusoidal waveform due to integrating in an integer n1 of sampling intervals different from the integer n of sampling intervals whose duration is the same as the entire period: (a) ε(%) vs. λ; (b) ε(%) vs. n1n.
Sensors 21 07967 g010
Figure 11. The relative error of the calculated active power expressed in %, due to the synchronization multimeters error ±1 μs: (a) Sinusoidal (50 Hz) voltage and 90° phase-fired current waveform; (b) Sinusoidal current with low harmonic content (THD < 2%) with sinusoidal voltage with 10% for h = 5th and 3% for h = 11th.
Figure 11. The relative error of the calculated active power expressed in %, due to the synchronization multimeters error ±1 μs: (a) Sinusoidal (50 Hz) voltage and 90° phase-fired current waveform; (b) Sinusoidal current with low harmonic content (THD < 2%) with sinusoidal voltage with 10% for h = 5th and 3% for h = 11th.
Sensors 21 07967 g011
Table 1. Reference measuring systems.
Table 1. Reference measuring systems.
Measuring SystemComponentTypePerformance
VoltageResistive dividerLCOE1060/1
Measuring cableRG-5850 Ω, 4 m
Matching impedanceLCOE1.1 MΩ
Multimeter (Master)Keithley DMM7510700 V, 10 A, 1 MS/s
CurrentFluxgate current sensorLEM ITN-900-S636 A, 1500/1 A
DC sourceLCOE±15 VDC
ShuntLCOE10 Ω
Multimeter (Slave)Keithley DMM7510700 V, 10 A, 1 MS/s
Energy calculationSoftwareLCOE
Table 2. Summary of expanded uncertainties U (k = 2).
Table 2. Summary of expanded uncertainties U (k = 2).
Expanded Uncertainty%
U[P] (%) 0.23
U[VRMS] (%) 0.15
U[IRMS] (%) 0.10
U[S] (%) 0.18
U[N] (%) 0.43
Table 3. Uncertainty budget for the active power for 90° phase-fired current waveform.
Table 3. Uncertainty budget for the active power for 90° phase-fired current waveform.
Influence Parameters
Xi
SymbolEstimateStandard Uncertainty
u(xi)
Probability DistributionSensitivity Coefficient
ci
Uncertainty Contribution
ui(y)
Digitization of Voltage
If VRMS(tj) = 10% VFS
δ v 1 % 00.012/√3% (I)Rectangular10.0069%
δ v 2 % 00.0075/√3% (I)Rectangular 0.82 10 = 9.06 0.039%
Digitization of Current
If IRMS(tj) = 10% IFS
δ v 1 % 00.012/√3% (I)Rectangular10.0069%
δ v 2 % 00.0075/√3% (I)Rectangular 1.62 10 = 12.73 0.055%
VD SF calibrat. uncertainty S F V D 10590.1/2% (II)Normal10.05%
VD SF drift δ 1 VD 00.02/√3% (I)Rectangular10.011%
VD temperature coefficient δ 2 , VD 00.0075/√3% (I)Rectangular10.0043%
VD short term stability δ 3 , VD 00.06/√3% (I)Rectangular10.035%
VD non-linearity δ 4 , VD 00.03/√3% (I)Rectangular10.017%
VD phase uncertainty δ 1 , t 058/2 μrad (II)Normal−6.26 × 10−5%/μrad0.0018%
VD residual phase error δ 2 , t 01452/2 μrad (II)Normal−6.26 × 10−5%/μrad0.045%
CT SF calibration uncertainty S F C T 15000.0070/2% (II)Normal10.0035%
CT SF drift δ 1 , C T 00.005/√3% (I)Rectangular10.0029%
CT temperature coefficient δ 2 , C T 00.006/√3% (I)Rectangular10.0034%
CT SF non-linearity δ 3 , C T 00.0036/√3% (I)Rectangular10.0021%
CT S. F. phase calibration δ 3 , t 058/2 μrad (II)Normal−6.26 × 10−5%/μrad0.0018%
CT S. F. phase error drift δ 4 , t 0293/√3 μrad (I)Rectangular−6.26 × 10−5%/μrad0.0106%
CT residual phase error δ 5 t 0668/2 μrad (II)Normal−6.26 × 10−5%/μrad0.021%
Shunt calibration uncertainty R s 9.979 Ω0.05/2% (II)Normal10.025%
Shunt drift δ 1 , R 00.01/√3% (I)Rectangular10.0057%
Shunt temperature coeffic. δ 2 , R 00.05/√3% (I)Rectangular10.029%
Shunt resistance variation vs. frequency δ 3 , R 00.01/√3% (I)Rectangular10.0057%
Shunt residual phase error δ 6 , t 0NegligibleNormal10%
V integration trapezoidal rule S p 1 , t 00.028/√3% (III)Rectangular10.0161%
I integration trapezoidal rule S p 1 , t 00.028/√3% (III)Rectangular10.0161%
V Samples taken for a period S p 2 , t 05·10−5/√3% (IV)Rectangular10.00003%
I Samples taken for a period S p 2 , t 05·10−5/√3% (IV)Rectangular10.00003%
Synchroniz. between DMM S p 3 , t 0314μrad/√3 (V)Rectangular−6.26 × 10−5%/μrad0.011%
Combined variance u c 2 ( y ) = u i 2 ( y ) 0.0136%
Combined standard uncert. u c ( y ) 0.117%
Expanded uncertainty U [ P ] ( % ) = κ u c ( y ) ( κ = 2 ) 0.23%
(I) Data from manufacturer data sheet. (II) Data from calibration certificate. (III) As indicated in Uncertainty Contribution Due to Trapezoidal Integration Rule ( S p 1 ), an error of 0.04% for active power is obtained for 90° phase-fired waveform using a sinusoidal voltage. This total error contribution is split in two equal contributions for both the voltage and the current, equal to √2 × 0.04%/2 = 0.028%. (IV) The value of 5 × 10−5% comes from Uncertainty Contribution Due to an Integration Time Different to the Complete Period ( S p 2 ). (V) 314μrad corresponds to the maximum synchronization error between the digital multimeters of 1 μs for 50 Hz, as indicated in Uncertainty Due to the Synchronization Error between Multimeters ( S p 3 ).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Khamlichi, A.; Garnacho, F.; Simon, P.; Rovira, J.; Ramirez, A. Calibration of On-Board Energy Measurement Systems Installed in Locomotives for AC Distorted Current and High Voltage Waveforms and Determination of Its Uncertainty Budget. Sensors 2021, 21, 7967. https://doi.org/10.3390/s21237967

AMA Style

Khamlichi A, Garnacho F, Simon P, Rovira J, Ramirez A. Calibration of On-Board Energy Measurement Systems Installed in Locomotives for AC Distorted Current and High Voltage Waveforms and Determination of Its Uncertainty Budget. Sensors. 2021; 21(23):7967. https://doi.org/10.3390/s21237967

Chicago/Turabian Style

Khamlichi, Abderrahim, Fernando Garnacho, Pascual Simon, Jorge Rovira, and Angel Ramirez. 2021. "Calibration of On-Board Energy Measurement Systems Installed in Locomotives for AC Distorted Current and High Voltage Waveforms and Determination of Its Uncertainty Budget" Sensors 21, no. 23: 7967. https://doi.org/10.3390/s21237967

APA Style

Khamlichi, A., Garnacho, F., Simon, P., Rovira, J., & Ramirez, A. (2021). Calibration of On-Board Energy Measurement Systems Installed in Locomotives for AC Distorted Current and High Voltage Waveforms and Determination of Its Uncertainty Budget. Sensors, 21(23), 7967. https://doi.org/10.3390/s21237967

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