20th Iranian Conference on Electrical Engineering, (ICEE2012), May 15-17, Tehran, Iran
Charging of plug-in electric vehicles: Stochastic modelling of load demand within domestic grids Ehsan Pashajavid
Masoud Aliakbar Golkar
K. N. Toosi University of Technology
K. N. Toosi University of Technology
Tehran, Iran
Tehran, Iran
[email protected]
[email protected]
Abstract-This paper proposes a stochastic approach based on
•
D eriving the hourly load demand of PEVs in order to determine the penetration level of the vehicles by evaluating the impacts on the grid.
•
Quantifying the load delivered through a transformer which can be useful for various distribution system applications such as network planning, load
Monte Carlo simulation to derive the load demand of a fleet of domestic commuter plug-in electric vehicles. At first, appropriate non-Gaussian probability density functions are fitted to the employed datasets to generate random samples required in the Monte Carlo simulation. The datasets include home arrival time, daily travelled distance and home departure time of randomly
management and probabilistic load flow as well as sitting and sizing issues [6].
selected private ICE vehicles. In each iteration, extraction of the charging profile is carried out for the individual PEVs in order to derive the hourly aggregated load profile of the fleet. Then,
The most effective factors on the load demand of PEVs are home arrival time, daily travelled distance, home departure time, driving habits and road traffic conditions [7]. Besides, battery capacity of PEVs and efficiency of the battery chargers can be noted as well. Some of the mentioned factors, due to the related uncertainties, are not deterministic. Thus, it is essential to employ stochastic and probabilistic modeling approaches [8].
probability density function of the aggregated load of the PEVs within each hour is estimated. Eventually, the expected value of the hourly load demand can be calculated regarding the achieved power distributions.
The PEVs are assumed to be charged
through a distribution transformer. Thus, profile of the power delivered through the transformer to the PEVs is attained which can be useful for various distribution system applications such as network planning, load management and probabilistic load flow as well as sitting and sizing issues.
This paper proposes a stochastic method based on Monte Carlo simulation to achieve two abovementioned goals, namely extracting the hourly aggregated load demand of PEVs and quantifying the power delivered to a fleet of PEVs through a domestic transformer. Home arrival time, daily travelled distance and home departure time of a set of commuter ICE
Keywords-Electric vehicles; stochastic modeling; smart grid; distribution network; load modeling.
I.
INTRODUCTION
Recently, a great attention has been paid to plug-in electric vehicles (PEVs) as one of the hopeful and effective means in order to tackle the energy crisis and to alleviate the environmental problems [1,2]. A considerably extensive amount of regulatory and technical efforts have been making to provide affordable PEVs. Therefore, a numerous number of PEVs will be undoubtedly on the road all around the world [3].
vehicles in Tehran is used as the input datasets. [n order to generate the required random samples, the aforesaid random variables (RVs) are modeled by probability density functions (pdfs) fitted to them. Unlike most of the earlier researches that used the Normal pdf for the sake of simplicity [9, I 0], this paper suggests appropriate non-Gaussian pdfs to fit to the RVs. In addition, the effectiveness of the modeled pdfs in comparison
PEVs can be generally divided into two categories, namely
with the Normal pdf is illustrated. Afterwards, the devised stochastic modeling algorithm with two scenarios to extract the initial state-of-charge (SOCinit) of PEV batteries is thoroughly explained. By applying the approach, the hourly load distribution functions of the PEVs is derived. Eventually, the
battery electric vehicles (BEV) which are pure electric vehicles and plug-in hybrid electric vehicles (PREV) that include both electric motors and internal combustion engines (ICEs) [4]. The possibility of bidirectional power transfer with the grid is one of the salient characteristics of PEVs. Actually, PEVs can absorb power from the grid when the state-of-charge (SOC) of the battery is low. On the other hand, it is possible for PEVs to inject power into the grid (V2G), especially, during the peak load hours [5]. However, V2G requires a longer time to be practical [1].
distributions are employed to extract the demand profile of the fleet. The remainder of this paper is organized as follows. The employed datasets are reviewed and modeled by fitting appropriate pdfs to them in section [I. Section III elaborates the proposed stochastic methodology. Afterwards, the simulation results are discussed in section [v. Eventually, the paper is summarized in section V.
Estimating the power consumption of the PEVs may be considered as one of the critical challenges. This issue can be addressed in two ways:
978-1-4673-1148-9112/$3l.00 ©2012 IEEE
535
II.
generated applying the extracted pdfs in section II. The available charging time (t.vi) for each of the PEVs is calculated by subtracting the departure time of the next day (dn+l) from
DATASETS OF THE ICE VEHICLES
The employed datasets in the modeling algorithm have been gathered using questionnaires filled-out by the randomly selected owners of the commuter light duty ICE vehicles in Tehran. The owners were asked to give merely the data. The datasets include home arrival time travelled distance (trd) and home departure time vehicles during weekdays. The datasets can be seen
the arrival time of today (a,J The battery capacity (Capbat), SOCinib power rating (Prat) and efficiency (Cchr) of the battery chargers determine the necessary time (tful1) to fully charge the battery. In case tfull would be less than t.vi, the complete charging of the battery can be accomplished. In any other way, it would be impossible to fully charge the battery. The bigger the battery capacity and the lesser the power ating of the
commuting (at), daily (dt) of the in Fig.l.
To generate random samples required in the Monte Carlo simulation, it is essential to fit pdfs to the datasets. Most of the earlier researches have used the Gaussian (Normal) pdf as a straightforward distribution that did not match with the datasets
charger, the longer time is necessary to fully charge the battery. Accordingly, the hourly power consumption of the PEVs can
and this led to inaccurate results. It is suggested in this paper to
of the load profile of the individual PEVs (DPIPn) is fulfilled in
fit appropriate non-Gaussian pdfs to the three mentioned RVs. In fact, a number of pdfs are tested on them and then, the best is selected. As may be seen in Fig. I (a), the Weibull pdf Ctdt(t)) is suggested as the most appropriate function to be fitted to the departure time RV as below:
order to estimate the demand profile of the fleet (DPFm).
be estimated after arriving at home. In each iteration, extraction
�
.
�
o
1
� :0
(I)
'"
-g 0.5
Ii. �
To model the daily travelled distance a type III Generalized expected value (Gev) pdf (2) is derived. The result is illustrated in Fig.l(b).
8.5 7.5 ) pa=rt u re'-' ti=me D e= '-i(=ho=u r'= O.06i---r___,---___r-' - -----,-��:====;_r �
.;;;
� � �
1
r.tid (d)
.
1 = -
(J,,;,
(1 + k tr,/
(d
(d-fitrd)� 1 -j1 ) -(1+-) -(l+k"d--) d 0"" tod) k", e ,
'"
(2)
c::::::J Data - - - Normal pdf
0.05 0.04
--
Gevpdf
0.03 0.02
Ii. e: 0.01
(JII;,
o
As it is shown in Fig. I (c), a type III Generalized expected value (Gev) pdf (3) is fitted to the home arrival time RV as well.
-
10
20
40
30
50
Trawlled distance (km)
0.6i-----,-----,-------,----'---�--,------;:;��===il
(3) Arrival time (hour)
Figure I.
Moreover, the extracted distribution of each RV is compared
The data and pdfs of (a) home departure time d" (b) daily travelled distance trd and (c) home arrival time a,.
with the Normal pdf fitted to the same dataset in order to verify the mentioned claim (Fig. I). It is seen that fitted non-Gaussian pdfs provide a better approximation of the original datasets. The parameters of the fitted pdfs are given in TABLE I.
TABLE 1.
In the next section, the extracted pdfs are utilized to generate the random samples.
III. A.
Datasets
MODELING METHODOLOGY
THE PARAMETERS OF THE FITTED PDFS The Normal pdf
a=7.67454
(T.vd, = 0.43178
jJ=21.3812
fly", = 21.4150 (TN'" = 8.58711
The proposed algorithm
The main objective of this paper is to extract the hourly load demand of the PEVs with uncoordinated charging. It is assumed that only the home charging is available and the PEVs are plugging into the grid as soon as their arrival. Figure 2 illustrates the overall procedure of the proposed approach. First of all, the random samples of at (�), trd (tr,,) and dt (d,,) are
at
536
The suggested pdf
flVd, = 7.48436
fly" = 17.7170 (TN" = 1.01385
k" = -0.052368 , fl", = 17.6568 (T", = 7.1222
k" = -0.060798 fl" = 17.2700 (T", = 0.84832
[t is notable that the PEVs are supposed to be charged through a distribution transformer. Therefore, profile of the power delivered to the PEVs through the transformer is attained as well. Considering [N as the iteration number of the Monte Carlo
DPF
simulation, the explained procedure is carried out for IN times in order to derive the distribution of the aggregated power consumption of the fleet within each hour (D APFh). Afterwards, the expected values of the hourly load demand of the PEVs can be calculated regarding the D APF". Eventually,
EN = Total number of PEVs IN = Iteration number m=1
demand profile of the fleet (DPF) is estimated by employing the extracted hourly expected values. B.
The battery initial SOC
Extraction of the battery SOC of the PEV s at the charging start time (SOCinit) is done with two following scenario cases:
Random generation (an, trdn, dn)
•
tavi = dn - an tfulln=(1 OO-SOCinitn)xCapbat/(1 OOxPratXCchr)
Case 1: As the worst case, the battery initial SOC of the PEVs (SOCinit) can be assumed to be a constant. This constant is determined according to the depth-of discharge (D OD ) of the PHEVs as follows:
SOC;n;,n =100-DOD
(4)
However, the mentioned assumption is not feasible in reality. Therefore, it is suggested to alter (4) by using randomly generated positive numbers through employing an exponential pdf as below:
{
\
SOCinim = 00
��OD
fexp (x ) = -- e f.1exp •
+ Rand (fexp (x
)) (5)
POXI'
Case 2: The SOCinit is determined based on the daily travelled distances. Regarding the fact that tavi in home charging is usually bigger than tfulb it is rational to assume that the battery SOC are 100% at the departure time of the PEVs. Hence, the SOCinit of the PEV can be derived as follows:
SOCmlln , =100
lrdn Ceif X Caphm
xlOO
(6)
where Ceff is the efficiency coefficient of the PEVs during driving which is dependent on the driving patterns and traffic conditions as well as driver efficiency of the electric motors.
IV.
Figure 2,
SIMULATIONS AND DISCUSSIONS
D ue to the assumption that the fleet is delivered power through a distribution transformer, it seems rational to perform the simulation for 20 PEVs. The above-explained procedure is repeated 10000 times to deal with the related uncertainties. The simulation parameters can be found in TABLE 11.
Flowchart of the proposed approach,
537
TABLE II.
THE SIMULATION PARAMETERS
10000 random samples. Then, the SOCinit is determined accordingly which can be considered as the worst case. The synthetic results can be seen in Fig. 4(a). As it is obvious, a number of the PEVs have fully charged batteries (SOCinitn=100%) correspondent to the untraveled vehicles.
10000
IN
EN
20
DOD
70% 3kW 20kWh 2km/kWh
Cchr
0.9
j.1cv..p
8
A.
Synthetic datasets
Figure 3 presents distribution of the randomly generated samples (the synthetic data) for PEVs by using the three suggested pdfs in section n. Each generated set includes 10000 samples. As it is obvious from Fig. 3(a), the generated samples of dt is distributed between 4:30 a.m. and 8:30 a.m. with the
C.
transformer. According to Fig. 2, the procedure is repeated for 10000 times in order to derive the distribution of the aggregated power consumption of the fleet within each hour (D APFh). The estimated D APFh at h=16 olclock, h=19 olclock and h=22 o'clock regarding two scenarios of the battery initial SOC explained in section IV.B is shown in Fig. 5(a), Fig. 5( b) and
The battery initial SOC of the PEVs
To derive the battery initial SOC of the PEVs the mentioned scenarios are employed: Case 1: Regarding to (5),
•
fexp(x)
Fig. 5(c) respectively. Moreover, the demand profile of the fleet (DPF) consisted by the expected values of the D APFh is presented in TABLE 111.
used to generate
IS
As may be expected, the DPF values calculated by adopting the first scenario (case1) are generally greater in comparison with the estimated results employing the second scenario (case 2).
PEVS
20 4000
g �
i3
o
:§:
3000
18
2000
4.5
4000
� i3 25
PEVs
16
1000
O L...L
g
._- , --,--,_ , 5
5.5
.oL_· �[_ ']____
•
6.5
oooo ��---,---�--�-�--�----,---__,
..' . _---_'--___ _ 8.5
7.5
" me (hour) __,__----,---,------ --' i- _ '-__,__-___,_--,---- --,
.---_
3000
5000
14
g
12
10
1000
- 1000
--
---
'--------
-L--.- ---·--� --.-'""50
Travelled distance (km)
4000
�
25
"'-
�
60
__:jo
4000
� 6 3000 (5 :§: 2000
2000
:0
g
The aggregated power consumption
The demand profile of the fleet (DPFm) is achieved by estimating the demand profile of the individual PEVs (DPIPn) in each iteration of the devised procedure. As it is mentioned in section lll, the estimated profile can be adopted as the profile of the power delivered to the PEVs through the distribution
highest num ber of occurrence between 7:30 a.m. and 8 a.m.. Fig. 3( b) presents the generated samples of trd distributed between 0 km and 76 km. 0 km means the PEV departure is not happened. Fig. 3(c) shows the generated samples of at distributed between 15 p.m. and 24 p.m.. The highest num ber of occurrence appears between 17 p.m. and 17:30 p.m..
B.
Case 2: The SOCinit is determined based on the daily travelled distances. According to (6), the SOCinitn is dependent on the travelled distance (trdn), the battery capacity (Capbat) and the efficiency coefficient of PEVs during driving (Ceff). The distribution of the synthetic initial SOCs of the PEVs is shown in Fig. 4( b).
•
14
o -- · ---- --40 30
�----.- ·---- ..�----'-50 80 70 80
12
-- ' --
90
100
SOCinil
3000��---,---�--�-�- -___,_--�-__.
80
2500
3000
2000
1000
500 � ' ----
0 15
16
----- ---- ' ---17
·--------"--...... 22 21 lime (hour)
�
o 30
24
40
50
60
70
80
90
100
sOeinit
Figure 4. Battery initial SOC generated by using (a) scenario No. 1, (b) scenario No. 2.
Figure 3. The randomly generated samples of (a) home departure time db (b) daily travelled distance trd and (c) home arrival time al.
538
10
2000 � 1500
domestic commuter plug-in electric vehicles has been devised. A series of datasets, including home arrival time, daily travelled distance and home departure time of randomly selected private ICE vehicles have been employed. To avoid any mismatching between the original dataset and their probabilistic models, appropriate non-Gaussian probability density functions has been fitted to them. Extraction of the charging profile has been performed in each iteration for the individual PEVs in order to estimate the hourly aggregated load profile of the fleet. Afterwards, probability density function of the aggregated load of the PEVs within each hour has been derived by applying the Monte Carlo simulation. The expected value of the hourly load demand has finally been
_Case 1 c::::::JCase 2
c
i"
8
8
1000
� 500 0 0
15
10
20
25 (kW)
Load demand
30
35
40
45
" 4000 g � 3000 "
82000 :£ 1000 0 40
calculated regarding the achieved power distributions. 42
44
46
48
Load demand
60
(kW)
REFERENCES
Load demand
(kW)
Figure 5. The distribution of the aggregated power consumption of the PEVs (DAPFh) regarding two scenarios of the battery initial SOC at (a) h=16 o'clock (b) h=19 o'clock, (c) h=22 o'clock.
TABLE III.
DEMAND PROFILE OF THE FLEET (DPF) DELIVERED THROUGH A DISTRIBUTION TRANSFORMER (W)
Hours
Case 1
Case 2
Hours
Case 1
Case 2
2 3 4 5 6 7 8 9 10 II 12
4806.55 ll22.78 212.01 32.38 3.66 0 0 0 0 0 0 0
3106.79 709.22 121.69 15.92 2.70 0 0 0 0 0 0 0
13 14 15 16 17 18 19 20 21 22 23 24
0 0 205.88 22623.76 52557.08 57000.38 56678.69 54044.72 48384.19 37408.21 26518.46 14151.84
0 0 199.01 22022.61 51048.34 54728.31 52702.32 48386.52 42582.69 3384l.04 22466.78 10270.33
V.
[1]
W. Su, H. Rahimi Eichi, W. Zeng, M. Chow, "A Survey on the Electrification of Transportation in a Smart Grid Environment", IEEE Transaction on Industrial Informatics, VOL. 8, NO. 1, pp. 1-10, FEB. 2012.
[2]
M. A Golkar, E. Pashajavid, "Analytical assessment of mutual impacts between PHEVs and power grid", 21st International Conference on Electricity Distribution (CIRED),Frankfurt,6-9 June 20ll.
[3]
AY. Saber, G.K. Venayagamoorthy, "One million plug-in electric vehicles on the road by 2015", Proceeding of the 12th International Conference on Intelligent Transportation Systems,2009,pp. 141- 147.
[4]
A Emadi, Young Joo Lee, K. Rajashekara, "Power Electronics and Motor Drives in Electric, Hybrid Electric, and Plug-In Hybrid Electric Vehicles", IEEE Transaction on Industrial Electronics, Vol. 55, No. 6, pp. 2237-2245,2008.
[5]
J.R. Pillai, B. Bak-Jensen, "Integration of Vehicle-to-Grid in the Western Danish Power System", IEEE Transaction on Sustainable Energy, Vol. 2,No. 1,pp. 1949-3029,Sep. 2010.
[6]
AH. Hajimiragha, C.A. Canizares, M.W. Fowler, S. Moazeni, A Elkamel, "A Robust Optimization Approach for Planning the Transition to Plug-in Hybrid Electric Vehicles", IEEE Transaction on Power Systems, pp. 0885-8950,Feb. 2011.
[7]
M. Ferdowsi, "Plug-in Hybrid Vehicles - A Vision for the Future",IEEE Vehicle Power and Propulsion Conference, VPPC,pp. 457-462,2007.
[8]
D. Wu, D.C. Aliprantis, K. Gkritza, "Electric Energy and Power Consumption by Light-Duty Plug-In Electric Vehicles", IEEE Transaction on Power Systems, Vol. 26,No. 2,pp. 738 - 746 ,2011.
[9]
L. Trovao, H.M. Jorge, "Power Demand Impacts of the Charging of Electric Vehicles on the Power Distribution Network in a Residential Area", Proceedings of the 2011 3rd International Youth Conference on Energetics (IYCE), 7-9 July 2011.
[10] K. Qian, C. Zhou, M. Allan, Y. Yuan, "Modeling of Load Demand Due to EV Battery Charging in Distribution Systems", IEEE Transaction on Power Systems, Vol. 26, No. 2, pp. 802 - 810, Aug. 2010.
CONCLUSION
Tn this contribution, a stochastic approach based on Monte Carlo simulation to derive the load demand of a fleet of
539