1
European Association for the Development of Renewable Energies, Environment and Power Quality (EA4EPQ)
International Conference on Renewable Energies and Power Quality (ICREPQ’10) Granada (Spain), 23rd to 25th March, 2010
Preparation of full paper for the International Conference on Renewableof Energies Deterministic and Probabilistic Assessment the and Power Quality Impact of the Electrical Vehicles on the Power Grid A. Author1, B. Author2 and C. Author1 1 E. Valsera-Naranjo1 , A. Sumper1,2 afila-Robles1 , A. Sudri`a-Andreu1,2 1 , P. Lloret-Gallego , R. Villaf´ Department of Electrical Engineering 1 Centre d’Innovaci´ o Tecnol`ogica en E.T.S.I.I., Vigo University Convertidors Est`atics Campus i Accionaments(CITCEA-UPC), Departament d’Enginyeria El`ectrica, of Lagoas – Marcosende, 36310 Vigo (Spain) Universitat Polit` e cnica de Catalunya. Phone/Fax number:+0034 986 812685, e-mail:
[email protected],
[email protected] EU d’Enginyeria T`ecnica Industrial de Barcelona, 2 European Association for the Development of Renewable Environment and Power Quality (EA4EPQ) Comte d’Urgell, 187; 08036 Barcelona,Energies, Spain. http://www.citcea.upc.edu 2 IREC CataloniaMailing address Institute for Energy Research Phone, fax, e-mail Josep Pla, B2, Pl. Baixa. 08019 Barcelona, Spain http://www.irec.cat
[email protected],
[email protected], A. Language
[email protected],
[email protected],
[email protected] Abstract. The abstract in placed in the left-hand column and starts at about 11 cm from de top of the page. It should give concise information on the aims, the methods, the results and the conclusions of the paper. Its length should not exceed 200 words. Abstract—This paper analyzes the impact of the charge of EVs (Electrical Vehicles) on a power grid. In order to simulate the behavior of the EV charging on the grid a model of its battery Key words has been found (EV’s charging curve). As the way of the EV is charged from the grid affects critically to the voltage levels Please the main keylines words your text to 5)are and to thewrite saturation of the two of modalities of (3 charge presented: separatednot-controlled by commas. charge and controlled charge. Finally, the simulations have been performed from two points of view: deterministic analysis and probabilistic analysis.
1. Introduction
Authors of acceptedI. abstracts will be requested to submit I NTRODUCTION a full paper for inclusion in the International Conference TheRenewable recent awareness about fuels and the environment on Energies andfossil Power Quality CDRom and hasalso arisen more sustainable alternatives regarding means an extended abstract of one sheet. The full papers are of transport. Consequently, hybrid vehicles and pure electrical due January 20, 2010. vehicle have become the main alternatives for green transportation. new trend hassent caused market or activation and it The full This papers have to be as a pdf-file a word file is (.doc) expected thatofhybrid and electrical will constitute to one the following e-mail vehicles addresses:
[email protected], majority in private
[email protected] transport. or: Moreover, one kind of hybrid vehicles, the plug-in electric
[email protected] vehicles (PHEV), will not only charge their batteries, but Pleasewill carefully to PHEV also beread ableandto follow inject these power instructions to the network ensure a uniform quality vehicles and appearance when required, as the electrical (EV) do. of Thisall fact contributions. suggests that EV penetration will affect current power system performance. Then, it is necessary to study some scenarios explain the nature of the previous work, in of Clearly penetration of such vehicles intoproblem, the electrical network purpose and contribution of the paper. order to maintain security and quality of power supply within standard limits. 2. Page Layout On one hand, in order to analyze the impact of EVs integration into the power grid, aspects related to storage Papers should notbatteries) exceed and six the double column A4 of technologies (namely charging process (21x29.7cm) pages including figures, tables and such storage devices have to be studied. In order to analyze diagrams. On the last page try to balance both columns.
English (check carefully for correct grammar and spelling)
B. Margins the impact of EV on the power grid, load flow calculations in a standard network performed following. doand that, Set the urban top and bottomare margins to 2 cm, and theToleft different penetrations of EV network are considered. right margins to 1.7 cm. into The the spacing between the two should be 1.25 are cm shown and commented lately. Thecolumns results of simulations On the other hand, an European standard which defines the Justification: Fullby EVs is needed. Until a new type of 1) connector to be used 2) Headers and Footers: None standard gets develop, it is proposed to use a SCHUKO (CEE 3) Page Numbering: None 7/4) connector type for currents up to 16A. In addition, to slow charging, the output values of the charging station should be C. 16Type Sizes up to A per plug, 230 V ± 10%, and 50 Hz ± 1% [1]. Manuscripts must be typed with a proportional serif E LECTRICAL V EHICLE typeface such II. as Time New Roman. The type sizes to be used for the different parts in the paper are shown in a In order to analyze the impact of the EVs into the grid, Table I. model of its charging curve is needed. From all the existing Table I. - Type Sizes
models of EV, the charging curve chosen is from the Mitsubishi I-MIEV [2]. This vehicle has an autonomy of 160 km SIZE APPEARANCE and Title its batteries are made of Li-ion (50Ah 16kWh 330V) 14 pt Bold centred [2]; Fig.1Author shows the charging curve of 11 theptEV when its capacity at Centred 20%.Abstract The battery is charging at full9 power during 2 hours and pt thenMain decreases its charging power exponentially since the hour text 10 pt 4; atPrimary this time it is considered that the battery is completely heading 12 pt Bold, left charged. Secondary headings 10 pt Italic, flush left The way headings of how EVs are charged theindented grid has a Tertiary 10 pt from Italic, critical influence on its impact on the voltage levels Tables 9 pt Centred and on the saturation of the lines. Considering this fact, two modalities Table and figure captions 9 pt Centred of charge are presented: not-controlled charge and controlled References 9 pt charge. In the not-controlled charge EVs start the charge as theyD.park. In Spacing the controlled charge, EVs only can charge during Line a determined period of time of the day (low-load period). In both modalities charge, the number EVs between charging Singled spaced, ofwith a double line of space paragraphs. Allow an extra half space above a line
60
40
20
0
40000 35000
2,00000
30000 25000
1,50000 20000 1,00000
15000 10000
0,50000 5000 0,00000
0
0,5
1
1,5
2
2,5
3
3,5
4
0 0:00 1:00 2:00 3:00 4:00 5:00 6:00 7:00 8:00 9:00 10:0011:0012:0013:0014:0015:0016:0017:0018:0019:0020:0021:0022:0023:00
Time [h]
(a) Not-controlled charge
Charging curve of the EV EEV's chargge power [MW]
Fig. 1.
Hour
45000
14
40000 35000
12
30000
10
25000 8 20000 6
15000
4
10000
EV's charge power (R= 0.1) EV's charge power (R= 0.2) Vehicle movement
2
5000
0
0 0:00 0:30 1:00 1:30 2:00 2:30 3:00 3:30 4:00 4:30 5:00 5:30 6:00 6:30 7:00 7:30 8:00 8:30 9:00 9:30 10:00 10:30 11:00 11:30 12:00 12:30 13:00 13:30 14:00 14:30 15:00 15:30 16:00 16:30 17:00 17:30 18:00 18:30 19:00 19:30 20:00 20:30 21:00 21:30 22:00 22:30 23:00 23:30
is defined establishing different scenarios of EV penetration over the vehicles acceding to the city (Fig.2). The curve of the movement of the vehicles has been obtained from an estimation of the curve of Barcelona [3] and has been shifted two hours earlier in order have a better representation of the Denmark’s vehicle movement (in this estimation, population has also been considered).
16
Vehicle e movement
Chargin ng power [%]
80
45000
EV's charge power (R= 0.1) EV's charge power (R= 0.2) Vehicle movement
2,50000
Vehicle e moveme ent
EEV's chargge power [MW]
3,00000 100
Hour
(b) Controlled charge Fig. 3. EV’s power demand for the not-controlled and controlled charge modality
800 700 600
EV
500
B006
B005
400
B012
B013 WT
CHP
300 200
External Grid
100
B001
B002 B003 B004
WT
0:00 0:30 1:00 1:30 2:00 2:30 3:00 3:30 4:00 4:30 5:00 5:30 6:00 6:30 7:00 7:30 8:00 8:30 9:00 9:30 10:00 10:30 11:00 11:30 12:00 12:30 13:00 13:30 14:00 14:30 15:00 15:30 16:00 16:30 17:00 17:30 18:00 18:30 19:00 19:30 20:00 20:30 21:00 21:30 22:00 22:30 23:00 23:30
L02
WT
L03 EV03
Fig. 4.
During the charging process of the EV, will be charging simultaneously EV with different states of charge. Therefore, a model is needed to take this fact into account in the steadystate simulations. Eq.1 models the superposition of the demand caused by the EV in different charging states, where R is penetration of EV, Pev is the maximum power which a car can charge, i is the counter associated to the actual hour, j is the counter associated to the previous hour, ∆vak,k−1 is the increase of EV between the hour i and the hour j and Cevk is the charging state of the EV (Fig.1).
B017
L01 EV01
L00 EV00
Movement of vehicles
i X
EV02
B016
Hour
P Hi = R · P ev ·
B015
B007 B008 B009
0
Fig. 2.
B014
B011
B010
Danish power grid
In order to study the most critical situation the Winter(week) case has been chosen from all the possible situations of electrical demand because it has the highest demand in every hour (Fig.5). Therefore, the wind power generation curve is from this period. Fig.6 depicts a representative daily generation for the three wind power generation units. 14 Winter (week) Summer (week) Winter (weekend) ( ) 12
Summer (weekend)
10
(∆vak,k−1 · Cevk )
(1) Power [MW W]
k=1
Applying this model to the vehicle movement curve, the demand curves for the both charging modalities can be obtained (Fig.3)
8
6
4
2
III. P OWER G RID 0
The scenario of the simulations is a part of the Danish sample grid (Fig.4) [4]. This grid, of 400 MW of short-circuit power, has 3 wind turbine generation units of 630 kW at busbars B013, B015 and B017. In addition, there are three combined cycle units of 3 MW each one at the busbar B005, but for purposes of the study has remained disconnected. Loads are at the busbars B005, B010, B011 and B012. Due to simplify the interpretation of the simulations and the results, EVs have been separated from the other consumptions.
0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
Hour [h]
Fig. 5.
Electrical demand
In the deterministic and in the probabilistic analysis the voltage results presented are from two busbar. The first busbar is B005 and it has been chosen because L00 is always the load with the highest power. The second busbar chosen is B014 in order to have the behavior of a closer busbar to a wind power 2
0,73 0,71
1,005
0,69 0,67
1
0,65 0,63 0,995
0,61 0,59
0,99
0,57 0,55 0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
Voltage [pu u]
Wiind Power Generattion [MW W]
0,75
23
Hour
0,985
0,98
Fig. 6.
Wind power generation 0,975
0,97
generation unit.
EV 0% EV 10% EV 20%
0,965 1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
Hour [h]
IV. D ETERMINISTIC A NALYSIS
Fig. 8. Voltage levels for the busbar B005 for the different penetrations of EV in the case of not-controlled charge for the deterministic analysis
In this section, the result from executing a deterministic analysis are shown. In this kind of analysis all loads are previously determined. The results from the two typologies of charge for the EV are presented for the different penetrations of EV.
1,05
A. Not-controlled charge
1
Fig.7 and Fig.8 show respectively the voltage levels for busbar B005 and B014 for the case of not-controlled charge. Voltage [pu u]
0,95
1,005
0,9
1 0,85
0,995 0,8
Voltage [pu u]
0,99
EV 0% EV 10% EV 20% 0,75
0,985
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
Hour [h] 0,98
Fig. 9. Voltage levels for the busbar B005 for the different penetrations of EV in the case of controlled charge for the deterministic analysis
0,975
0,97
EV 0% EV 10% EV 20%
0,965 1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
Hour [h]
1,05
Fig. 7. Voltage levels for the busbar B005 for the different penetrations of EV in the case of not-controlled charge for the deterministic analysis
1
0,95
Voltage [pu u]
B. Controlled charge Fig.9 and Fig.10 show respectively the voltage levels for busbar B005 and B014 for the case of controlled charge.
0,9
0,85
0,8
V. P ROBABILISTIC A NALYSIS
EV 0% EV 10%
From the compiled data during the week days of the winter season a mean and an standard deviation have been found for each hour of the day. In order to perform a probabilistic analysis 50 load flows for each hour have been executed generating random powers with the mean and the standard deviation founded previously.
EV 20% 0,75 1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
Hour [h]
Fig. 10. Voltage levels for the busbar B005 for the different penetrations of EV in the case of controlled charge for the deterministic analysis
3
A. Not-controlled charge 1,005
Fig.11 and Fig.12 show mean voltage levels for the different penetrations of EV in the case of not-controlled charge. In both figures are the calculated mean voltage levels from the probabilistic analysis and the mean from the deterministic analysis. In order to have a most detailed view of the critical busbar, in this case B014, Fig.13 show the the range of voltage level that can be reached at this busbar.
1
0,995
Voltage [pu u]
0,99
1,005
0,985
0,98
0,975
0 97 0,97
1
MIN
0,965
AVG
0 995 0,995
MAX 0,96
Voltage [pu u]
0
1
2
3
4
5
6
7
8
9
10
0,99
11
12
13
14
15
16
17
18
19
20
21
22
23
Hour [h]
Fig. 13. Voltage levels for the busbar B005 for the different penetrations of EV in the case of not-controlled charge
0,985
0,98 EV 0% (mean) EV 10% (mean) EV 20% (mean)
0,975
EV 0% (measured mean) EV 10% (measured mean) 1,05
EV 20% (measured mean) 0,97 1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
Hora [h] 1
Fig. 11. Voltage levels for the busbar B005 for the different penetrations of EV in the case of not-controlled charge Voltage [pu u]
0,95
1,005
0,9
0,85
1
EV 0% (mean) EV 10% (mean)
0 995 0,995
0,8
EV 20% (mean) EV 0% (measured mean)
Voltage [pu u]
EV 10% (measured mean) 0,99
EV 20% (measured mean) 0,75 1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
Hora [h]
0,985
Fig. 14. Voltage levels for the busbar B005 for the different penetrations of EV in the case of controlled charge
0,98 EV 0% (mean) EV 10% (mean) EV 20% (mean)
0,975
EV 0% (measured mean) EV 10% (measured mean) EV 20% (measured mean) 0,97 1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24 1,05
Hora [h]
Fig. 12. Voltage levels for the busbar B014 for the different penetrations of EV in the case of not-controlled charge
1
Voltage [pu u]
0,95
B. Controlled charge Fig.14 and Fig.15 show mean voltage levels for the different penetrations of EV in the case of controlled charge. As happens in the not-controlled charge, the B014 is the busbar with the lowest voltage levels, Fig.16 show the the range of voltage level that can be reached at this busbar.
0,9
0,85
EV 0% (mean) EV 10% (mean) 0,8
EV 20% (mean) EV 0% (measured mean) EV 10% (measured mean) EV 20% (measured mean)
0,75 1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
Hora [h]
VI. C ONCLUSION
Fig. 15. Voltage levels for the busbar B005 for the different penetrations of EV in the case of not-controlled charge
The objective of the paper is to analyze the impact of charging EVs from a Danish grid. Therefore, different penetrations of EV and charge modalities are proposed. In order to cover 4
1,05
deterministic analysis and probabilistic analysis. On one hand, results from the simulations show that notcontrolled charge amplifies the demand of the line at the hours of the higher electrical demands. On the other hand, simulations show that with the controlled charge a high penetration of EVs can be charged without the need of investments in the power grid infrastructure.
1
0,95
Voltage [pu u]
0,9
0,85
0,8
R EFERENCES
0,75
0,7
[1] Valsera-Naranjo Sumper Lloret-Gallego Villaf´afila-Robles, Sudri`aAndreu. Electrical vehicles: state of art and issues for their connection to the network. Electrical Power Quality and Utilisation, 2009. 1 [2] http://www.mitsubishi-motors.com/special/ev/index.html. 1 [3] Dades B`asiques 2008 v.02. Ajuntament de Barcelona. Direcci´o de Serveis de Mobilitat., 2008. 2 [4] Roberto Villaf´afila-Robles. Probabilistic modeling in normal operation and the control of distribution systems with renewable source based DG units. PhD thesis, Universitat Polit`ecnica de Catalunya, 2009. 2
MIN AVG MAX
0,65 0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
Hour [h]
Fig. 16. Voltage levels for the busbar B005 for the different penetrations of EV in the case of not-controlled charge
the maximum situations two analysis have been performed:
5