Performance of Heuristic Optimization in Coordination of Plug-In Electric Vehicles Charging

A heuristic load management (H-LMA) algorithm is presented for coordination of Plug-in Electric Vehicles (PEVs) in distribution networks to minimize system losses and regulate bus voltages. The impacts of optimization period T (varied from 15 minutes to 24 hours) and optimization time interval (varied 15 minutes to one hour) on the performance, accuracy and speed of the H-LMA is investigated through detailed simulations considering enormous scenarios. PEV coordination is performed by considering substation transformer loading while taking PEV owner priorities into consideration. Starting with the highest priority consumers, HLMA will use time intervals to distribute PEV charging within three designated high, medium and low priority time zones to minimize total system losses over period T while maintaining network operation criteria such as power generation and bus voltages within their permissible limits. Simulation results generated in MATLAB are presented for a 449 node distribution network populated with PEVs in residential feeders.


Introduction
Preliminary studies by Amin et al. (2005), Amin (2008) and Lightner et al. (2010) indicate that Plug-In Electric Vehicles (PEVs) will dominate the market in the near future as pollution-free alternatives to the conventional petroleum-based transportation.However, according to Moses et al. (2010), Masoum et al. (2011) and Moses et al. (2012), uncoordinated PEV charging specially at high penetration levels during the peak load hours may cause undesirable impacts on the power grid such as unpredictable system peaks, unaccepted voltage deviations, significant increases in losses and poor power quality, as well as overloading of the distribution and substation transformers.This has motivated researchers to propose and implement different PEV coordination algorithms.
In general, PEV chargers can be controlled to operate in charge or discharge modes with the energy being transferred from grid to vehicle (V2G) or from vehicle to grid (G2V), respectively.One of the first approaches for PEV coordination based deterministic and stochastic dynamic International Journal of Renewable Energy and Biofuels 2 _______________ Sara Deilami, Amir S. Masoum, Mohammad A.S. Masoum and A. Abu-Siada (2013), International Journal of Renewable Energy and Biofuels, DOI: 10.5171/2013.898203programing was presented by Clement-Nyns et al. (2010).Masoum et al. (2011) performed peak load saving with PEV coordination without considering the random nature of PEV arrivals and departures.A relatively fast PEV coordination algorithm suitable for online applications is proposed by Deilami et al (2011).Ashtari et al. (2012) predicted PEV charging profiles and electrical range reliability based on recorded vehicle usage data.Wu et al (2012) designed a minimumcost load scheduling algorithm based on the forecasted electricity price and PEV power demands.In a research study by Khodayar et al. (2012) and Zhao et al. (2012), the PEV coordination problem is performed considering the impact of wind power generation.Wen et al. (2012) and Ma et al. (2013) have presented decentralized charging control algorithms considering large populations of PEVs.There are also many documents investigating the operation of PEVs in V2G mode to support the grid by performing frequency regulation and/or energy storage including the research performed by Bashash et al. (2012), Han (2010) and Sortomme et al. (2011).This paper will first present a heuristic load management algorithm (H-LMA) to coordinate PEV charging activities while reducing system losses and regulating bus voltages over a 24 hour period.Then, simulation results will be presented for a 449 node distribution network populated with PEVs in residential feeders.Finally, the impacts of heuristic optimization parameters including optimization period T and optimization time interval t ∆ on the accuracy and speed of H-LMA will be investigated.

Problem Formulation
PEV charge coordination is a constrained optimization problem that could be solved by using online algorithms (i.e., PEV coordination is performed as soon as vehicles are randomly plugged-in) or offline schemes (i.e., all vehicles are assumed to be plugged-in according to their pre-known/forecasted charging patterns).This paper assumes the charging patterns of all PEVs are known or forecasted and utilizes a heuristic approach to solve the optimization problem.
The optimization problem objective function is formulated based on the minimization of total system power losses: Where t ∆ and T are the optimization time interval and period used for loss minimization.loss t P is the system power loss at time t (computed using the Newtonbased power flow), k V is voltage of node k at time t, and n is total number of nodes while y + are resistance and admittance of line section between nodes k and k+1.
PEV coordination constraints are node voltage limits and system demand level at time t:

Heuristic Load Management Algorithm (H-LMA)
A MATLAB based algorithm has been developed to perform PEV scheduling based on H-LMA (Fig. 1).The algorithm will perform loss minimization over the optimization period T using time interval t ∆ based on Eqs. 1-2while considering the system constraints (Eqs.3-4).Three charging time zones are defined: • Red charging zone (18:00h-22:00h) coinciding with most of the on-peak period and is designated for highpriority PEV owners willing to pay higher tariff rates in order to charge their vehicles as soon as possible.
• Blue charging zone (18:00h-01:00h) intended for medium-priority consumers that prefer to charge their vehicles at partially off-peak periods and pay lower tariff rates.
• Green charging zone (18:00h-08:00h) when most PEV charging will probably take place due to the cheapest tariff rates as most low-priority consumers will require their vehicles fully charged for the following day.
The algorithm assumes all PEVs are plugged in at 18:00 (6pm).It begins by first reading the input parameters (e.g., bus and branch impedance data, nodes with PEVs, optimization period T, optimization time interval t ∆ , designated priority time zones, load profiles for PEV chargers and residential loads as well as system constraints) and performing initialization (e.g., selecting the highest priority group, time zone and PEV).
The main program loop is progressing from high to low PEV priority groups (e.g., red zone to green zone).Within the selected priority group, individual PEVs are temporarily activated to determine system performance at all possible PEV nodes and charging time combinations within that priority charging time zone.From these combinations, the algorithm selects the PEV within the group and the charging start time resulting in the minimum system losses, taking into consideration the charging duration and the current demand level.The physical node location at which PEV charging occurs is an important factor as it impacts the load flow, power losses in the cables and system voltage profile.Therefore, the H-LMA determines the PEV node and charging time that will result in the least system losses (Eq.2).
If at any time the load flow indicates a constraint violation at any node (Eqs.3-4), the algorithm will try the next possible charging start time such that the constraints are satisfied.Therefore, it may not be possible for all PEV owners to be accommodated in their preferred charging zones and must be deferred to the next possible hour.Once it has been determined which PEV node in that priority group can begin charging and at what time resulting in minimum system losses, the selected PEV scheduling is permanently assigned and the system load curve updated ready for the next iteration.This process is repeated for all nodes in that priority group before advancing to the next prioritycharging zone (e.g., blue zone subscribers).
At the end of this process, the H-LMA arrives at individual schedules assigned to all PEV chargers.The program then exits the main loop and computes the 24 hour load flow to print new system performances (e.g., all node voltage profiles and power losses).

Smart Grid Test System
The selected test system is a modification of the IEEE 31 bus 23 kV distribution systems (Deilami et al (2011)) combined with 22 residential 19 nodes LV 415 V networks populated with PEVs.The resulting 449 node system is supplied from the HV main bus via a 23kV/415V 100 kVA distribution transformer as shown in Fig. 2. System data are listed in the Appendix.
The peak power consumption of a house is assumed to be on average 2 kW with a power factor of 0.9.Four PEV penetration levels are selected including 16% (with nodes "o", "b" and "q" randomly designated with red, blue and green priorities, respectively), 32% (with nodes "o", "b, r" and "f, h, q" randomly designated with red, blue and green priorities, respectively), 47% (with nodes "o", "b, j, r" and "f, g, h, m, q" randomly designated with red, blue and green priorities, respectively) and 63% (with nodes "o, s", "b, d, j, r" and "f, g, h, k, m, q" randomly designated with red, blue and green priorities, respectively).
For this study, a 10 kWh battery capacity per PEV with a depth of discharge (DOD) of 70% and battery charger efficiency of 88% is assumed (Deilami et al (2011)) which will require a total of 8 kWh of energy from the grid to charge a single PEV.A standard single-phase 240V outlet (Australia) can typically supply a maximum of 2.4 kW.
There are also 15A and 20A outlets (singlephase and three-phase) which can supply approximately 4 kW and 14.4 kW, respectively.In this paper, a fixed charging power of 4 kW is used.

Simulation Results and Discussion
Simulation results for uncoordinated and coordinated (using H-LMA of Fig. 1) PEV charging for the smart grid system of Fig. 2 are presented in Figs.3-5 and Tables 1 and2.

A. Case A: Random PEV Charging
Simulation results of Fig. 3

B. Case B: H-LMA Coordinated PEV Charging
Coordinated PEV charging is performed with (Fig. 4-5) and without (Table 1) PEV owner preferred time zone priorities.Compared to Case A, a significant improvement in smart grid performance is achieved.Most notably, the system demand peak has been reduced (Figs. 3(a) and 4(a)) which is more favorable from a standpoint of generation dispatch and preventing overloads.
Comparison of results also indicate the significant impacts of coordinated (H-LMA) PEV charging on voltage profile where the unacceptable voltage deviations of about 17% (Fig. 3(b)) at the worst bus for uncoordinated PEV charging is compensated to less than 10% (Fig. 4(b)) which is within the regulation limits.However, there is a trade off in that a few PEV subscribers who designated a preferred priority charging time zone were not accommodated in their requested charging zone (Fig. 4(a)) because the system reached a point where PEV loading caused voltage regulation to be violated.H-LMA handled these cases by attempting to schedule the PEV owners causing the violations to a charging time where the system is not under strain, thereby satisfying constraints.The improvements in system efficiency with H-LMA coordination strategy are also evident in Table 1.Energy losses for the high penetration (63%) with H-LMA are limited to 2.59% of system consumption versus the worst uncoordinated charging scenario with losses of 3.09%.Furthermore, peak power losses are limited to less than a third of the worst case random uncoordinated charging (Fig. 4(c)).The H-LMA charging also has positive impacts on peak transformer load currents.
For many of the uncoordinated random charging scenarios (Table 1), distribution transformers are experiencing load currents of up to 0.88 pu, while with H-LMA coordination, transformer currents are reduced to levels of approximately 0.54 pu (Table 1).

Coordination
Detailed simulations are presented and compared in Table 2 to highlight impacts of t ∆ and T (Eq. 1) on the performance of H- LMA.In general, the speed and accuracy of the PEV coordination algorithms will depend on the selection of optimization time interval ( t ∆ ) and period (T).
The accuracy can be improved by using shorter time intervals (e.g., checking the status of PEVs and network as quickly as possible based on online information and measurements available through smart meters) and performing loss minimization over a long period (e.g., 24 hours).However, the drawback is the computing time will dramatically increase, especially in realistic large smart grids with many nodes and high penetration levels of PEVs.Therefore, a compromise should be made between the solution accuracy and computation time considering system size and the anticipated PEV penetration level.
Based on the results of Table 2, the practical options may be to use moderate time intervals with large optimization periods for offline PEV coordination (e.g., t ∆ =60 min and T=24 hours for applications where all vehicles are plugged-in or their charging patterns are known/forecasted before the start of optimization) and select small values for online PEV coordination (e.g., t ∆ =T=15 min to start charging batteries as soon as vehicles are randomly plugged-in).

Conclusion
This paper investigates the impacts of optimization parameters including optimization period T and optimization time interval t ∆ on the accuracy and the speed of a heuristic load management algorithm (H-LMA) that coordinates PEV charging activities while reducing system losses and regulating bus voltages over a 24 hour period.Main conclusions are: • H-LMA will limit overall system overloads and voltage fluctuations while reducing stress on distribution circuits such as cables and transformers.
• The speed and accuracy of H-LMA depend on the selected values for T and t ∆ .
• It is showed that optimization accuracy can be improved by using shorter time intervals and performing loss minimization over long periods (e.g,.24 hours).This will however, require long computing times.Therefore, a compromise should be made between the solution accuracy and the associated computation time considering system size and the anticipated PEV penetration levels.

Sara Deilami ,
Fig. 1: Proposed H-LMA for Coordination of Pevs to Minimize Total System Losses over Period T Using Optimization Time Interval t ∆ Considering Node Voltage Profiles and Maximum Demand Level Sara Deilami, Amir S. Masoum, Mohammad A.S. Masoum and A. Abu-Siada (2013), International Journal of Renewable Energy and Biofuels, DOI: 10.5171/2013.898203

Fig. 4 :
Fig. 2: The 449 Node Smart Grid Test System Consisting of the IEEE 31 Node 23 Kv System with Several 415 V Residential Feeders.Each Low Voltage Residential Network Has 19 Nodes Representing Customer Households Populated with Pevs Randomly Arriving within 24 Hours.

Table 1 : Comparison of Simulation Results for Uncoordinated and Coordinated (H-Lma, =15 Min, T=24 Hours) Pev Charging for the Smart Grid Test System of Fig. 3: Pevs are Assumed to be Randomly Arriving at Each Time Interval. For Comparison, Consumer Priorities are not considered and the Same Gaussian Random Distributions are Used in the Simulations.
*) Ratio of system losses over 24 hours compared to total power consumption over 24 hours.**) Voltage devataion at the worst bus.***) Maximum of all distribution transformer load current.