Deterministic and Probabilistic Assessment of the Impact of ... - icrepq

8 downloads 0 Views 1MB Size Report
[2] http://www.mitsubishi-motors.com/special/ev/index.html. 1. [3] Dades B`asiques 2008 v.02. Ajuntament de Barcelona. Direcció de Serveis de Mobilitat., 2008.
1

International Conference on Renewable Energies and Power Quality (ICREPQ’10) Granada (Spain), 23rd to 25th March, 2010

European Association for the Development of Renewable Energies, Environment and Power Quality (EA4EPQ)

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

English (check carefully for correct grammar and spelling)

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. Page Layout 2. 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. https://doi.org/10.24084/repqj08.704

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 extra half space a line 1505 paragraphs. Allow an RE&PQJ, Vol.1, No.8,above April 2010

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. https://doi.org/10.24084/repqj08.704

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

1506

2

RE&PQJ, Vol.1, No.8, April 2010

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. https://doi.org/10.24084/repqj08.704

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

1507

3

RE&PQJ, Vol.1, No.8, April 2010

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

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 https://doi.org/10.24084/repqj08.704

2

1508

4

RE&PQJ, Vol.1, No.8, April 2010

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:

https://doi.org/10.24084/repqj08.704

1509

5

RE&PQJ, Vol.1, No.8, April 2010