Computer Modeling in Engineering & Sciences

DOI: 10.32604/cmes.2022.017958

ARTICLE

A Novel Bidirectional Interaction Model and Electric Energy Measuring Scheme of EVs for V2G with Distorted Power Loads

Jiarui Cui1,2,*, Qing Li1,*, Bin Cao2,3, Xiangquan Li1 and Qun Yan1

1School of Automation and Electrical Engineering, University of Science and Technology Beijing, Beijing, 100083, China
2Guiyang Aluminum Magnesium Design and Research Institute Co., Ltd., Guiyang, 550081, China
3Chinalco Intelligent Technology Development Co., Ltd., Hangzhou, 311199, China
*Corresponding Authors: Jiarui Cui. Email: cuijiarui@ustb.edu.cn; Qing Li. Email: liqing@ies.ustb.edu.cn
Received: 19 June 2021; Accepted: 02 September 2021

1  Introduction

Electric vehicles (EVs) play a critical role in meeting the environmental goals of the Sustainable Development Scenario to reduce local air pollution and to address climate change [1]. With the increasing demand for petroleum resources and environmental issues, new energy electric vehicles are increasingly being put into use in the last decade. The number of electric cars on the world’s roads had swelled to 7.2 million by 2019 from only about 17000 in 2010 [2]. EVs could actively provide about 460 GWh of energy to the grid at suitable times via vehicle-to-grid services (V2G). But every coin has two sides, and there are a lot of challenges to unlock the full flexibility potential of electric vehicles, such as technologies, regulatory, market frameworks, and so on.

V2G technology has been widely concerned by experts and scholars all over the world, and has achieved a series of results which mainly focus on scheduled charging and discharging [3–6], bidirectional converter [7–9], power energy management system [10,11], battery lifecycle [12], modeling approaches [13], energy storage systems [14] and so on. It is pointed out in [15] that a large amount of energy storage needs to consider the fast-dynamic response of power demand and a modular method is given. A green V2G network is proposed in [16] for energy management which consists of three planes: management plane, control plane, and data plane. Literature [17,18] mainly studies aspects of data security and communication security. V2G technology enables electric vehicles to deliver power to the power system to provide complementary capabilities. Literature [17] established a novel and secure communication framework to meet the real-time and security of power transportation. Literature [18] proposed a member-based secure data sharing scheme in V2G networks to protect data privacy and save system resources. In [19], a battery state based scheduling algorithm is proposed to optimize the V2G charging and discharging strategy. In [20], a construction scheme based on transfer learning is proposed to update the safe charging and discharging strategy dynamically. These documents have played a pivotal role in the development of V2G technology. In [21], based on grid load time-sharing pricing, an electric vehicle scheduling model was established with the goal of reducing charging cost and network loss cost. The charging and discharging strategies of electric vehicles were obtained by real-time calculation of power flow calculation and convex optimization algorithm. Literature [22] proposes a V2G management system consisting of an ordered charging dispatch system and a V2G converter control system. The system can fully utilize the flexibility of electric vehicle load while meeting the charging needs of electric vehicles. In [23], an ordered charge and discharge control model in V2G mode is proposed, in order to reduce the charge and discharge cost of electric vehicles and the variance of grid load.

However, electric energy measuring methods for V2G have not been concerned widely, although there have been a large number of mature modeling methods [24–26], filtering algorithms [27–29], measurement techniques [30–32] and they have been successfully applied to actual systems in the field of industrial processes. Based on the idea of [33], we established the charging model and discharging model of EVs for V2G, and given the electric energy measuring formula in [34]. But they did not consider the effect of the distorted power loads that is always stayed in power grid when EVs are charging or discharging. Furthermore, the unified bidirectional model and the corresponded electric energy measuring methods did not discuss. With the widespread use of V2G technology, it will be a very effective means to increase the participation of EV owners in V2G. So, this paper proposes a bidirectional electric energy measuring scheme for V2G with distorted power loads which considered the great randomness of the charging and discharging time of EVs and the distorted power loads on the power grid.

The rest of this paper is presented as follows. Section 2 establishes the bidirectional interaction model of EVs for V2G based on electric network theory. In Section 3, a bidirectional electric energy measuring scheme is proposed for V2G with distorted power loads based on power flow analysis and energy metering theory. Section 4 gives the simulation results and analyzes the data. Finally, the conclusions are shown in Section 5.

2  A Bidirectional Interaction Model of EVs for V2G

The signals’ type is very complex because of the large number of non-linear loads in the grid. In addition to harmonic signals and inter-harmonic signals, there are many other forms of distorted signals such as impact signals and fluctuation signals which are typical in grid. They cannot be described by harmonic models. We can establish their model as follows using modulation theory:

ik(t)=(1+αk(t))Imsin⁡(ω0t+ϕk) (1)

where k is the k-th impact or fluctuation, αk(t) is a random function, ϕk are random variables, Im is the amplitude of current. In order to facilitate the measurement, the zero point of the measurement time is based on the time at which αk(t) occurs.

If the signal is an impact signal, αk(t) can be approximated as follows:

αk(t)=Mke−tτk,(Mk>0) (2)

If the signal is a fluctuating signal, αk(t) can be approximated as follows:

αk(t)={Mk,0<t<τk,0, other ,(Mk>0) (3)

where the duration τk is random since occur time of the impact or fluctuation signals is uncertain. Mk is the amplitude of the impact or fluctuating signal.

Under the condition of distortion signal, the voltage at a certain measuring point in the power grid and the current flowing into the measuring point can be expressed as follows:

u(t)=uI(t)+uS(t) (4)

i(t)=iI(t)+iS(t) (5)

where uI(t) and iI(t) are the fundamental voltage at metering point and the current flowing into the metering point, respectively, uS(t) and iS(t) are the distortion voltage and the current, respectively.

On the basis of [33] and [34], a bidirectional interaction model of EVs for V2G is established as shown in Fig. 1.

Figure 1: Bidirectional interaction model of electric vehicles in the power grid

In Fig. 1, the left side of the circuit is a simplified model of the grid, where u(t) is the grid voltage and i(t) is the grid current. Zl1,Zl2,Zl3 are line loads in the grid. Z is the other load in the grid. The right side of the circuit is an electric car, where E is the electric vehicle battery electromotive force, Zev is the internal load of electric vehicle. ua(t) is the voltage at the metering point a, and ia(t) is the current flowing into the metering point a. uev(t) is the output voltage when the electric vehicle is discharged to the grid. If there are no other loads Z, then the model is simplified to the model in [34].

3  A Bidirectional Electric Energy Measuring Scheme

Assuming that the grid voltage is assumed to be as follows:

u(t)=Umsin⁡(ω0t+ψk) (6)

where voltage phase ψk is a random variable and Um is the amplitude of voltage.

Assuming the current direction shown in Fig. 1 is the positive direction, that is, the charging direction of the electric vehicle is the positive direction. According to formula (1), the current at the metering point a when charging the electric vehicle is as follows:

ia(t)=(1+α(t))Imsin⁡(ω0t+ϕk) (7)

When the electric vehicle is discharged, the current at the measuring point a is as follows:

ia(t)=−(1+α(t))Imsin⁡(ω0t+ϕk) (8)

When the electric vehicle is charged in the power grid, the voltage at the measuring point a is as follows:

ua(t)=u(t)−Zl2ia(t) (9)

When the electric vehicle is discharged into the power grid, the voltage at the measuring point a is as follows:

ua(t)=uev(t)+Zl3ia(t) (10)

In fact, the output voltage uev(t) is the same as the grid voltage when the EVs are discharged, and is a standard sinusoidal voltage.

uev(t)=u(t) (11)

Assuming that Zl2 and Zl3 are pure resistive loads Rl. In the case of V2G charging or discharging, the voltage at the measuring point a can be expressed as follows:

ua(t)=u(t)−Rl⋅(1+αk(t))Imksin(ω0t+ϕk) (12)

Take the k-th impact and fluctuation signal as an example for analysis. For the k-th impact or fluctuation, the voltage at the measuring point a is as follows:

uak(t)=Umsin⁡(ω0t+ψk)−Rl(1+αk(t))Imsin⁡(ω0t+φk)=Umsin⁡(ω0t+ψk)−RlImsin⁡(ω0t+φk)−Rlαk(t)Imsin⁡(ω0t+φk)=Umsin⁡ω0tcos⁡ψk+Umsin⁡ψkcos⁡ω0t−RlImsin⁡ω0tcos⁡φk−RlImsin⁡φkcos⁡ω0t−RlImαk(t)sin⁡(ω0t+φk)=uakI(t)+uaks(t) (13)

According to formula (13), the voltage at the metering point a can be divided into a fundamental voltage and a distortion voltage. The fundamental voltage at the metering point a is as follows:

uakI(t)=Umsin⁡ω0tcos⁡ψk+Umsin⁡ψkcos⁡ω0t−RlImsin⁡ω0tcos⁡φk−RlImsin⁡φkcos⁡ω0t=(Umcos⁡ψk−RlImcos⁡φk)sin⁡ω0t+(Umsin⁡ψk−RlImsin⁡φk)cos⁡ω0t=UakImsin⁡(ω0t+θk) (14)

The fundamental voltage amplitude is as follows:

UakIm=(Umcos⁡ψk−RlImcos⁡φk)2+(Umsin⁡ψk−RlImsin⁡φk)2=Um2+Rl2Im2−2UmRlImcos⁡(φk−ψk) (15)

The fundamental voltage phase is as follows:

θk=arctan⁡Umsin⁡ψk−RlImsin⁡ϕkUmcos⁡ψk−RlImcos⁡ϕk (16)

The distortion voltage at the metering point a is as follows:

uakS(t)=−RlImαk(t)sin⁡(ω0t+ϕk) (17)

According to Eq. (5), the current flowing into the metering point a can be divided into a fundamental current and a distortion current. The fundamental current is as follows:

iakI(t)=−Imsin(ω0t+ϕk) (18)

The distortion current is as follows:

iakS(t)=−Imαk(t)sin(ω0t+ϕk) (19)

Assume that within the measurement time (T=N⋅T0, where T0 is the fundamental period.), the charging time of the electric vehicle is T1, and the discharging time is T2. Theoretically, there is P1>P2, where P1 is the total power charged and P2 is the total power discharged. In order to facilitate the measurement, it is assumed that in the measurement time T, the N1 times impact occurs in the charging process, and the charging measurement time is T1=N1T0, the N2 times impact occurs in the discharging process, and the discharging measurement time is T2=N2T0. The average power at point a is as follows:

Pa=1T∑k=1N1+N2∫0T1+T2pak(t)dt=1T(∑k=1N1∫0T1pak−G2V(t)dt+∑k=1N2∫T2T1+T2pak−V2G(t)dt)=1T(Pa−G2VT1+Pa−V2GT2) (20)

Therefore, the electrical energy at point a is as follows:

Wa=TPa=Pa−G2VT1+Pa−V2GT2=(PI_G2V+PIS−G2V+PSI−G2V+PS−G2V)T1+(PI−V2G+PIS−V2G+PSI−V2G+PS−V2G)T2=(PI_G2VT1+PI−V2GT2)+(PIS−G2VT1+PIS−V2GT2)+(PSI−G2VT1+PSI−G2VT2)+(PS−G2VT1+PS−V2GT2)=WI+WIS+WSI+WS (21)

According to the formula (21), the electrical energy at the point a can be divided into four parts:

WI=WI_G2V+WI_V2G=PI_G2VT1+PI_V2GT2 (22)

WIS=WIS_G2V+WIS_V2G=PIS_G2VT1+PIS_V2GT2 (23)

WSI=WSI_G2V+WSI_V2G=PSI_G2VT1+PSI_V2GT2 (24)

WS=WS_G2V+WS_V2G=PS_G2VT1+PS_V2GT2 (25)

Taking the discharge situation of electric vehicles as an example, this paper will analyze the four parts of electrical energy in turn.

3.1 Fundamental Power

Combining Eqs. (14) and (18), the fundamental components of the average power of point a is as follows:

PI_V2G=∑k=1N21T2∫T1T1+T2uakI(t)iakI(t)dt=−∑k=1N21T2∫T1T1+T2{UakImsin(ω0t+θk)⋅Imsin⁡(ω0t+φk)}dt=12∑k=1N21T2∫T1T1+T2UakImImcos⁡(2ω0t+θk+φk)dt−12∑k=1N21T2∫T1T1+T2UakImImcos⁡(θk−φk)dt (26)

In Eq. (26), the first part is the second harmonic component of the power, and the integral mean is zero. Additionally, there is UakIm>0, Im>0 and cos⁡(θk−ϕk), so the first part must be positive. In summary, the following conclusion can be drawn:

WI_V2G=T2PI_V2G=−12∑k=1N2∫T1T1+T2UkIm_V2GImcos⁡(θk_V2G−φk)dt<0 (27)

3.2 Power Generated by Fundamental Voltage and Distortion Current

Combining Eqs. (14) and (19), the instantaneous power generated by the fundamental voltage and distortion current at point a is as follows:

uakI(t)iaks(t)=−UakImsin⁡(ω0t+θk)⋅Imαk(t)sin⁡(ω0t+ϕk) (28)

For the impact or fluctuating signals, the cumulative time of energy can be considered to be infinite. Therefore, the electric energy generated by the fundamental voltage and the distorted current at point a is as follows:

WIS_V2G=−∑k=1N2∫0∞UakImsin(ω0t+θk)⋅Imαk(t)sin(ω0t+φk)dt=12∑k=1N2∫0∞UakIm⋅Imαk(t)[cos⁡(2ω0t+θk+φk)−cos⁡(θk−φk)]dt=12∑k=1N2∫0∞UakIm⋅Imαk(t)cos⁡(2ω0t+θk+φk)dt−12∑k=1N2∫0∞UakIm⋅Imαk(t)cos⁡(θk−φk)dt (29)

The first term in Eq. (29) is denoted by A, and the second term is denoted by B. When αk(t) is an impact signal, combining Eqs. (2) and (29), A and B can be simplified as follows:

A=12∑k=1N2UakImImMk∫0∞e−1τkcos⁡(2ω0t+θk+φk)dt=14∑k=1N2UakImImMk∫0∞e−1τk(ej(2ω0t+θk+φk)+e−j(2ω0t+θk+φk))dt=14∑k=1N2UakImImMk[ej(θk+φk)∫0∞e(2jω0−1τk)tdt+e−j(θk+φk)∫0∞e−(2jω0+1τk)tdt]=14∑k=1N2UakImImMk[e−j(θk+φk)2jω0+1τk−ej(θk+φk)2jω0−1τk]=14∑k=1N2UakImImMkτk(2jω0τk−1)e−j(θk+φk)−(2jω0τk+1)ej(θk+φk)(2jω0τk+1)(2jω0τk−1)=14∑k=1N2UakImImMkτk2jω0τk(−2jsin⁡(θk+φk))−2cos⁡(θk+φk)−(4ω02τk2+1)=12∑k=1N2UakImImMkτkcos⁡(θk+φk)−2ω0τksin⁡(θk+φk)4ω02τk2+1=12∑k=1N2UakImImMkτksin⁡ηk4ω02τk2+1 (30)

B=12∑k=1N2UakImImMkcos⁡(θk−φk)∫0∞e−1τkdt=12∑k=1N2UakImImMkcos⁡(θk−φk)⋅(−τk)⋅e−1τk|0∞=12∑k=1N2UakImImMkτkcos⁡(θk−φk) (31)

In Eq. (30), ηk is as follows:

ηk=arctan⁡(12ω0τk)−ϕk−θk (32)

In Eq. (31), all parameters are positive, so there is B>0. A is a small value relative to B. Moreover, A is a random variable whose cumulative of multiple impacts tends to zero. Therefore, the following conclusion can be drawn:

WIS_V2G=A−B<0 (33)

When αk(t) is a fluctuating signal, combining Eqs. (3) and (29), WIS_V2G can be simplified as follows:

WIS−V2G=12∑k=1N2Ua/HmImMk∫0τk(cos⁡(2ω0t+θk+φk)−cos⁡(θk−φk))dt=∑k=1N2UaifmImMk4ω0sin⁡(2ω0t+θk+φk)|0τk−12∑k=1N2UadimImMkτkcos⁡(θk−φk)=∑k=1N2Uai^mImMk4ω0sin⁡(2ω0τk+θk+φk)−∑k=1N2Ua/imImMk4ω0sin⁡(θk+φk)−12∑k=1N2Ua/imImMkτkcos⁡(θk−φk) (34)

In Eq. (34), the first two terms are random variables, and the sum of multiple fluctuations will tend to zero. Similar to Eq. (31), the last term must be positive. Therefore, combining Eqs. (33) and (34), the following conclusion can be drawn:

WIS_V2G<0 (35)

3.3 Power Generated by Distortion Voltage and Fundamental Current

Combining Eqs. (17) and (18), the average power generated by the distortion voltage and fundamental current of point a is as follows:

PSI−V2G=∑k=1N21T2∫T1T1+T2(−RlImαk(t)sin⁡(ω0t+φk)⋅(−Imsin⁡(ω0t+φk)))dt=12∑k=1N21T2∫T1T1+T2RlIm2αk(t)(1−cos⁡(2ω0t+2φk))dt>12∑k=1N21T2∫T1T1+T2RlIm2αk(t)dt−12∑k=1N21T2∫T1T1+T2RlIm2αk(t)dt=0 (36)

Therefore, the following conclusion can be drawn:

WSI_V2G=T2PSI_V2G>0 (37)

3.4 Distortion Power

Combining (17) and (19), the average distortion power at point a is as follows:

PS−V2G=∑k=1N21T2∫T1T1+T2{−RlImαk(t)sin⁡(ω0t+φk)⋅(−Imαk(t)sin⁡(ω0t+φk))}dt=12∑k=1N21T2∫T1T1+T2RlIm2αk(t)2(1−cos⁡(2ω0t+2φk))dt>−12∑k=1N21T2∫T1T1+T2RlIm2αk(t)2dt+12∑k=1N21T2∫T1T1+T2RlIm2αk(t)2dt=0 (38)

Therefore, the following conclusion can be drawn:

WS_V2G>0 (39)

Based on the above analysis, the following conclusions can be drawn when the electric vehicle is discharged:

WI_V2G<0, the fundamental energy is negative. Electric vehicles supply power to the grid, so this part of the electricity should be measured.

WIS_V2G<0, the energy generated by the fundamental voltage and the distorted current is negative. This part of the electrical energy is also provided by the electric vehicle, so it needs to be measured.

WSI_V2G>0, the energy generated by the distortion voltage and the fundamental current is positive. The grid side returns electric energy to the electric vehicle in the form of a fundamental current, which should be used to offset the total absorbed energy and should be metered.

WS_V2G>0, the distortion energy is positive. Although this part of the electric energy is fed back to the electric vehicle by the grid side, it is fed back in the form of distorted current, which is harmful to electric vehicles. This part of the electrical energy should not offset the total absorbed energy and should not be measured.

Energy analysis when electric vehicle charging is similar to that of discharging, and can be directly obtained by combining literature [34]. Therefore, this article will not be deduced in detail.

According to Table 1, the total electrical energy at the metering point is as follows:

W=WG2V+WV2G=WI−G2V+WIS−G2V+WSI−G2V+WI−V2G+WIS−V2G+WSI−V2G (40)

However, current metering devices simply add up four types of electrical energy without taking into account the effects of distorted electrical energy on the grid and electric vehicles. Currently, the energy at metering point a is as follows:

Wa=WI−G2V+WIS−G2V+WSI−G2V+WS−G2V+WI−v2G+WIS−V2G+WSI−v2G+WS−V2G (41)

Combined with Eqs. (40) and (41), a reasonable measurement formula can be described as follows:

W=Wa−WS−G2V−WS−V2G=Wa−WS (42)

4  Simulation and Result Verification

A reasonable measurement method for electric vehicles during charging has been verified in literature [34]. Next, this paper will verify the reasonable measurement method during discharging. The simulation was performed in the MATLAB 2015b environment. Assume the line load Zl=Rl=1 Ω. Grid voltage is standard sinusoidal voltage. f0=50 Hz, Um=220 V as follows:

u(t)=220sin⁡(100πt) (43)

During the measurement time, the number of occurrences of the distortion signal, the time of occurrence, the duration, and the amplitude are random. The amplitude of the impact signal can be several times the amplitude of the current. The duration can be several seconds or even tens of seconds. Assume that the fundamental period is T0=0.02 s. Suppose the number of occurrences of an impact signal or the fluctuation signal is N=64. There is T=NT0.

4.1 Simulation Case 1

It is assumed that only the impact signal is included in the distortion current. When an impact signal is present, according to Eqs. (2) and (8), the current at the metering point a can be expressed as follows:

ia(t)=−(1+α1(t))Imsin⁡(ω0t+ϕ) (44)

α1(t)={2e−200t,t∈[t01,t01+0.005]0, other (45)

where t01 is the time at which α1(t) occurs. Let the current amplitude Im=10 A flowing into the metering point a, and the currents’ phase ϕ=−π/50. Therefore, the current expression flowing into the metering point a can be described as follows:

ia(t)={−10sin⁡(100πt−π50)−20e−200tsin⁡(100πt−π50),t∈[t01,t01+0.005]−10sin⁡(100πt−π50), other  (46)

The voltage at point a is as follows:

ua(t)={{220sin⁡(100πt)−10sin⁡(100πt−π50)−20e−200tsin⁡(100πt−π50)},t∈[t01,t01+0.005]220sin⁡(100πt)−10sin⁡(100πt−π50), other  (47)

Taking several random impact signals as an example, the current and voltage at the metering point a when the electric vehicle is discharged are drawn, as shown in Fig. 2.

As can be seen from Fig. 2, the current distortion caused by impact signals can cause the voltage to drop. The simulation results of the electric vehicle discharge energy under the impact signal condition can be obtained, as shown in Table 2.

Figure 2: Discharge current and voltage under the condition of impact signals. (a) Current flowing into the metering point a; (b) Voltage at the metering point a

As can be seen from the data in Table 2, the simulation results of the four-part electric energy during discharge are consistent with the analysis results in Table 1, that is, WI_V2G<0, WIS_V2G<0, WSI_V2G>0 and WS_V2G>0. Therefore, under the impact signal condition, the rationality of the proposed measurement method can be verified.

4.2 Simulation Case 2

It is assumed that only the fluctuating signal is included in the distortion current. When a fluctuating signal is present, according to Eqs. (3) and (8), the current at the metering point a can be expressed as follows:

ia(t)=−(1+α2(t))Imsin⁡(ω0t+ϕ) (48)

α2(t)={0.5,t∈[t02,t02+0.06]0, other (49)

where t02 is the time at which α2(t) occurs. Let the current amplitude Im=10 A flowing into the metering point a, and the currents’ phase ϕ=−π/50. Therefore, the current expression flowing into the metering point a can be described as follows:

ia(t)={−15sin⁡(100πt−π50),t∈[t0,t0+0.06]−10sin⁡(100πt−π50), other (50)

The voltage at point a is as follows:

ua(t)={220sin⁡(100πt)−15sin⁡(100πt−π50),t∈[t0,t0+0.06]220sin⁡(100πt)−10sin⁡(100πt−π50),other  (51)

Taking several random fluctuating signals as an example, the current and voltage at the metering point a when the electric vehicle is discharged are drawn, as shown in Fig. 3.

Figure 3: Discharge current and voltage under the condition of fluctuating signals. (a) Current flowing into the metering point a; (b) Voltage at the metering point a

As can be seen from Fig. 3, the fluctuating signals can cause the metering point current to rise and the voltage to drop. Assuming that the number of occurrences of the fluctuating signal is N=64, the simulation results of the electric vehicle discharge energy under the fluctuating signal condition can be obtained, as shown in Table 3.

As can be seen from Table 3, under the fluctuating signal condition, there is WI_V2G<0, WIS_V2G<0, WSI_V2G>0 and WS_V2G>0, which is consistent with the theoretical analysis in Table 1 when electric vehicles discharging.

4.3 Simulation Case 3

Assuming that the distortion signal does not contain the impact signal or the fluctuating signal, the distortion signal degenerates into a harmonic signal. Under the condition of harmonic signals, the current expression flowing into the metering point a can be described as follows:

ia(t)=−10sin⁡(100πt−π50)−5.5sin⁡(300πt+π6)−3.5sin⁡(100πt+2π5)−2.5sin⁡(700πt+π4) (52)

The voltage at point a is as follows:

ua(t)=220sin⁡(100πt)−10sin⁡(100πt−π50)−5.5sin⁡(300πt+π6)−3.5sin⁡(100πt+2π5)−2.5sin⁡(700πt+π4) (53)

According to the above, the current and voltage at the metering point a when the electric vehicle is discharged can be drawn, as shown in Fig. 4. And the simulation results are shown in Table 4.

As can be seen from Table 4, under the harmonic signal condition, there is WI_V2G<0, WIS_V2G=0, WSI_V2G=0 and WS_V2G>0. Therefore, the proposed energy metering method is applicable not only to impact signals and fluctuating signals, but also to harmonic signals.

Figure 4: Discharge current and voltage under the condition of harmonic signals. (a) Current flowing into the metering point a; (b) Voltage at the metering point a

In the V2G system, the reasonable measurement method of electric vehicle discharging energy at the metering point is as follows:

WV2G=WI_V2G+WIS_V2G+WSI_V2G (54)

The reasonable measurement of electric vehicle charging energy has been verified in literature [34], as follows:

WG2V=WI_G2V+WIS_G2V+WSI_G2V (55)

When the circuit contains only harmonic signals, there is:

WIS_V2G=WSI_V2G=WIS_G2V=WSI_G2V=0 (56)

In summary, a reasonable method of measuring electric energy is as follows:

W=WI−G2V+WIS−G2V+WSI−G2V+WI−V2G+WIS−V2G+WSI−V2G=Wa−WS−G2V−WS−V2G=Wa−WS (57)

where Wa is the electric energy metered by the electric energy meter, and WS is the total distorted electric energy.

5  Conclusions

First, the article points out that the traditional energy metering method in the grid under the distortion signal cannot measure correctly for V2G. Secondly, this paper establishes a load model for the bidirectional interaction of electric vehicles that can be used for energy metering. The power at the measuring point is analyzed, and a reasonable electric energy metering formula is summarized. The V2G bidirectional interaction model can express the charging and discharging process of the electric vehicle more practical and calculate the electric energy more reasonably at the metering point in a certain period.

Although a reasonable metering method is given, meeting the electric energies more accurately based on the V2G bidirectional interaction model is an important issue for engineering applications. The following issues are our future work:

1.    Accurate metering methods based on modern signal processing theories, such as wavelet, time-frequency analysis, and so on.

2.    Metering instruments development with real-time demands based on embedded technologies.

3.    The bidirectional control scheme of EVs for V2G and G2V applications.

Funding Statement: This work is supported by China Postdoctoral Science Foundation (2021M690798), Guizhou Province Science and Technology Plan Project (No. [2021] General 085), National Natural Science Foundation of China (No. 61603034), the Fundamental Research Funds for the Central Universities (Nos. FRF-BD-19-002A, FRF-DF-20-14).

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.

This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.