Next Article in Journal
Optimal Scheduling of Combined Heat and Power Generation Units Using the Thermal Inertia of the Connected District Heating Grid as Energy Storage
Next Article in Special Issue
Optimal Decision-Making to Charge Electric Vehicles in Heterogeneous Networks: Stackelberg Game Approach
Previous Article in Journal
Flow and Fast Fourier Transform Analyses for Tip Clearance Effect in an Operating Kaplan Turbine
Previous Article in Special Issue
Applications of a Strong Track Filter and LDA for On-Line Identification of a Switched Reluctance Machine Stator Inter-Turn Shorted-Circuit Fault
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Intelligent Energy Management Algorithms for EV-charging Scheduling with Consideration of Multiple EV Charging Modes

1
State Key Laboratory of HVDC, Electric Power Research Institute, China Southern Power Grid, Guangzhou 510663, China
2
School of Electrical and Electronic Engineering, Nanyang Technological University, 50 Nanyang Avenue, Singapore 639798, Singapore
*
Author to whom correspondence should be addressed.
Energies 2019, 12(2), 265; https://doi.org/10.3390/en12020265
Submission received: 30 November 2018 / Revised: 3 January 2019 / Accepted: 10 January 2019 / Published: 16 January 2019
(This article belongs to the Special Issue Optimization Methods Applied to Power Systems)

Abstract

:
Electric vehicles (EVs) are now attracting increasing interest from both industries and countries as an environmentally friendly and energy efficient mode of travel. Therefore, the EV charging and/or discharging issue has become an important challenge and research topic in power systems in recent years. An advanced and economic EV charging process, however, should employ smart scheduling, which depends on effective and robust algorithms. To that end, a comprehensive intelligent scatter search (ISS) algorithm within the frame of a basic scatter search has been designed with both unidirectional and bidirectional charging considered. The ISS structure also supports both a flexible and constant charging power rate by respectively employing filter-SQP (sequential quadratic programming) and mixed-integer SQP as local solvers with module control. The detailed design of ISS is presented and the objectives of smoothing the daily load profile and minimizing the charging cost have been tested. Compared with methods based on GS (global search), GA (genetic algorithm), and PSO (particle swarm optimization), the outcome-verified ISS can produce attractive results with a significantly short computational time. Moreover, to handle a large scale EV charging scenario, a hybrid method comprised of a GA and ISS approach has been further developed. Simulation results also verified its prominent performance, plus superbly low computational time.

1. Introduction

The prominent advantages of electric vehicles (EVs) lie in their high energy transfer efficiency and lack of carbon dioxide (CO2) and other air pollutant emissions. The transition to replace petrol vehicles by EVs, however, is never an easy task. From the aspect of power systems, conventional power distribution networks are typically designed for non-mobile loads, meaning that current power system infrastructures may not be resilient enough to accommodate the integration of EVs [1]. For instance, the daily load peak could be amplified by the aggregated load of unregulated EVs [2,3], which is certainly a threat for power grids. Other potential adverse impacts including voltage drops, frequency oscillations, network congestions, power imbalances, etc., as well concern for power system operators [4,5]. In addition, the implementation of V2G (Vehicle-to-grid) also affects the dynamics and performance of the entire power grid [6,7].
Reinforcing power infrastructures is a typical approach for addressing the above challenges, yet such a method appears to be inefficient and environmentally unfriendly. Meanwhile, it is also impossible to upgrade current power facilities for the accommodation of EV penetration, due to the expensive investment. As a more economical and attractive alternative, smart scheduling strategies can be employed to postpone the construction of power infrastructures to support EV charging demand, through which the goal of wise load regulation as well as economic benefits can be reached [8,9,10]. Nevertheless, owing to the differences in operational characteristics between EV and the power system, smart EV scheduling becomes a rather complex task that must rely on an efficient and robust charging algorithm [11]. In general, the operation of EVs should take into account, but not be limited to, a charging and/or discharging power rate limit, initial state of charge (SOC), customer travelling habits, final energy demand, and battery capacity. On the other hand, the operation of a power system is constrained from the limitations of generation units, network structure, transformer capacities, voltage and frequency requirements, etc. In addition, power grids themselves are dynamic networks with unpredictable internal complexities and volatilities, indicating that real-time charging schemes are more practical and superior than merely day-ahead negotiation. As of now, several algorithms using combinations of different solving techniques have also been developed to solve real-time EV scheduling problems [12,13]. With all the above considerations, the demand for powerful algorithms to perform intelligent EV charging scheduling is apparent [14].
As of now, several approaches or algorithms attempting to regulate EV charging have been designed and reported. Typically, the constraint functions related to EV scheduling are usually, or particularly, linear. For example, the EV charging/discharging power should be bounded within the lower and upper limits, whereas the aggregated load must be constrained by the power supply. In the meantime, the objective functions (e.g., cost minimization, load smoothing) can be also formed or reformulated as convex functions. Hence, the EV scheduling issue can be modeled as a convex optimization problem and solved through CVX [15], GAMS [16] or CPLEX [17]. However, these tools come with license concerns, i.e., one must pay for their license to legally use these tools. On the other hand, many conventional methods, such as the alternating direction method of multipliers [18] and tree-based dynamic programming [19], have also been employed to solve the intractable issues of EVs. With multi-constraints and multi-objectives, this kind of approach often divides the problem into several sub-problems, thereafter solving them over each stage. When the number of EVs is increased, however, these methods can get stuck, due to “curse of dimensionality” [20]. This can limit their usage in large-scale EV charging. Furthermore, some heuristic algorithms inspired from nature process such as a particle swarm optimization (PSO) [21], CRO (chemical reaction optimization) [22], and genetic algorithm (GA) [23,24] have also been utilized. These artificial intelligent algorithms, as verified, can outperform conventional methods in terms of computation overhead, with promising results. For example, the PSO approach in [21] for implementing EVs demand response is almost 2700 times faster than the MINLP (mixed integer nonlinear programming) technique under different scenarios. The deficiency of such methods is that they are intrinsically stochastic and can stagnate prematurely into local optima.
Although many algorithms have been utilized for scheduling EV charging, the emerging approaches are limited in a sense that the methods either only consider EV charging by assigning its states at an arbitrary value from 0 to a specified charging power limit (flexible power rate charging) [10,21], can just manage the EV charging pattern through adjusting the charging status (constant power rate charging) [25], or focus merely on unidirectional charging. Nevertheless, in the long run, both V2G and a flexible charging rate can occur simultaneously with unidirectional charging and a constant power charging rate, due to the fact that the future smart grid will accommodate various applications and services [26]. All the above issues indicate that a generic algorithm is desirable to meet the current EV charging demand and also support future bidirectional charging.
Hence, the aim of this work is to design a comprehensive and universal algorithm that can fit multiple charging modes and diverse charging rate scenarios for distribution-side management. With such complex considerations, an intelligent scatter search (ISS) algorithm framework has been designed. In principal, this novel method is essentially a hybrid SS (scatter search) framework, integrated with sequential quadratic programming (SQP), based local solvers. SS in nature is a population-based evolutionary approach [27]. Instead of depending on massive random components and traditional evolution techniques, like mutation and crossover, SS only concentrates on a small solution set. Besides, SS is attractive for its well-designed flexibility, enabling various local search methods to complete different optimization tasks. In the design, SQP techniques are chosen as the local solvers to cope with different charging power rates, since they are efficient in solving nonlinear constrained problems. More specifically, a trust-region SQP method based on a filter technique (filter-SQP) [28] is adopted for solving the flexible power rate charging, while the MISQP (Mixed-integer SQP) technique evolved in [29] is in charge of constant power rate charging.
The main contributions of this work lie in: the general SS framework is redefined and adapted to be a new algorithm framework; for the scheduling of EV charging, two algorithms are proposed and applied, i.e., an ISS framework to deal with single EV charging, and a GA-ISS method comprised of GA theory and the proposed ISS approach for massive EV charging; and the proposed algorithms can support both V2G and G2V (grid-to-vehicle) for different power rate configurations.
For comprehensive comparison, various algorithms are compared in the simulation part, specifically:
(1)
For the ISS framework, GA, PSO, and global search (GS) are compared;
(2)
For the GA-ISS framework, three methods are compared, including: a global control that calculates all the variables and constraints through CVX; DCM (dumbing control method) that charges all EVs as soon as they are plugged in; and a GA-PSO hybrid method.
The following paper is organized as follows: the modeling of the constraints and objective functions for the optimization of EV charging is described in Section 2, the detail designs of the proposed algorithms are presented in Section 3, simulation results to demonstrate and verify the effectiveness of the proposed ISS and GA-ISS are described in Section 4, and finally, a conclusion is given in Section 5.

2. Modelling for the Optimization of EV Charging

2.1. Assumptions

Normally, an EV completes its charging within a certain time window, rather than a whole day. Within this time frame, the charging status (either charging, discharging, or idling) together with the charging power rate, are required to be determined for each timeslot. The following assumptions have been made for the algorithm development:
(1)
All participated EVs are with the ideal charging/discharging efficiency of 100%.
(2)
The charging price and discharging tariffs follow the spot prices of basic load at each timeslot.
(3)
The time period of study starts from 8:00AM to 8:00AM of following day with 15 min per interval. Therefore, the time window is divided into 96 timeslots.

2.2. Definition of Charging Modes

Since both unidirectional and bidirectional charging, combined with either flexible or constant charging power rate, are considered, to better elaborate the design, four charging modes are herein defined:
(1)
CD-F mode: charging/discharging with a flexible charging rate,
(2)
CD-C mode: charging/discharging with a constant charging rate,
(3)
C-F mode: charging only with a flexible charging rate,
(4)
C-C mode: charging only with a constant charging rate.
The detailed comparisons for the different charging modes are listed in Table 1. Based on the type of variables, the charging problems can be transformed into NLP (nonlinear programming) or MINLP sub-problems, which will be introduced later.

2.3. Constraints

With practical considerations, EV scheduling is subjected to its charging power limit, battery capacity, and customer demand, while the total load is constrained by the supply capacity of power grid.
(1) Charging power limit
Generally, the value of P E V , j max is equal to the battery nominal charging power rate, P E V , j min is set to be zero for charging only cases while set to be same as P E V , j max when V2G is permitted, i.e.,
P E V , j min P E V , j i P E V , j max i ST , j SE .
(2) Restriction of the battery capacity
For the sake of prolonging the lifetime of EV batteries, S o c j , i is often limited between S o c j , min and S o c j , max . Take commonly used Li-ion batteries as an example, the range S o c j , min and S o c j , max are normally set to be 20% and 90% of their nominal state-of-charge (SOC) E R C , j , respectively [30]. Therefore, the following constraints are considered:
S o c j , min S o c j , i S o c j , max E R C , j ,
S o c j , i = S o c j , s t a r t + 1 C E V , j i = s t a r t c u r r e n t P E V , j i .
(3) SOC expectation level for journey requirement
To fulfill the journey requirement, the final energy state S o c j , end of the battery must reach E exp , j , i.e.,
S o c j , end = S o c j , s t a r t + 1 C E V , j i = s t a r t e n d P E V , j i = E exp , j .
(4) Capacity constraint of power grid
To ensure normal operation of the power grid, the total load cannot exceed the maximum capacity of generation units, i.e.,
P t o t a l i P g i ,
P t o t a l i = P l d i + j SE P E V , j i .

2.4. Objective Functions

In terms of optimization objectives, the emerging applications for EV regulation normally consider: (i) load curve flattening, (ii) energy cost minimization, (iii) running cost reduction for the power system, (iv) benefit maximization for aggregators, (v) reduction of CO2 emissions, and (vi) combinations of the above. The goal of multiple purpose optimizations can be set by utilizing weight factors, based on the emphases of different entities. For exemplification’s sake, only the first two objectives are considered to verify the proposed scheduling algorithm.
(1) Load curve flattening
This optimization objective aims to minimize the power load fluctuation, and the standard deviation of load profile can be utilized for this purpose:
Minimize F o b j , 1 = [ 1 n 1 i ST ( P t o t a l i P t o t a l A v e ) 2 ] 1 / 2 ,
where P t o t a l A v e denotes the daily load mean value and calculated as:
P t o t a l A v e = 1 n i ST P t o t a l i .
(2) Energy cost minimization
The objective function to minimize power energy cost is given as:
Minimize F o b j , 2 = i ST P t o t a l i C e i
The two mostly studied pricing strategies, namely time-of-use (TOU) price and real-time price, will be used for the simulation to determine the electricity price. The former usually sets the price C T O U i for different timeslots, while the latter determines an electricity tariff strictly according to time and load demand. For convenient purposes, the linear-price (LP) given in Equation (10) is adopted to model the dynamic real-time price. By incorporating Equation (10) into Equation (9), the final form of the function for cost minimization with the LP price is obtained and given in Equation (11), which is in a quadratic polynomial form.
C L P i = ψ P t o t a l i + γ
Minimize F o b j , 2 = i ST P t o t a l i C L P i = i ST P t o t a l i ( ψ P t o t a l i + γ ) = i ST ( ψ P t o t a l i 2 + γ P t o t a l i )
As can be seen, the optimization objectives in Equations (7), (9), and (11) are in three different forms, which will be separately simulated.

2.5. Transform the EV Charging into NLP/MINLP Problem

As previously introduced, the constraints considered for EV charging are nonlinear. Therefore, the EV charging cases are transformed into NLP or MINLP optimization problems [16,31], which consist of an objective function, some box bound constraints, several equal, and unequal constrained functions. Generally, the NLP problem has the form:
minimize F o b j ( x ) , subject to : H e q ( x ) = 0 G i e q ( x ) 0 l b x u b
where x represents the controllable vectors, F o b j ( x ) denotes the objective function, H e q ( x ) and G i e q ( x ) represent the equality and inequality constraints, and l b and u b are the lower and upper box limits respectively, all with compatible dimensions.
In this work, the issue of flexible charging rate (CD-F and C-F mode) is therefore transformed into an NLP problem. For the condition of constant charging rate (CD-C and C-C mode), only the charging status of EVs needs to be settled, and the Equation (1) thus evolves into Equation (13), and the optimization model turns into MINLP.
P E V , j i = P E V , j i S E V , j i i ST , j SE S E V , j i = { 1 , charging 1 , discharging , if permitted 0 , otherwise

3. Proposed ISS and GA-ISS Algorithm Frameworks

There are two algorithms proposed in this work: ISS for single EV charging and GA-ISS for massive EV charging. The proposed ISS algorithm is based on a basic scatter search (SS) framework with advanced local solvers. Therefore, to better understand the proposed algorithms, the principle of SS and the utilized local search solvers will first be described, followed by the ISS framework, and the hybrid method GA-ISS handling massive EVs.

3.1. Basics of Scatter Search (SS)

SS is a novel evolutionary-based method, since it normally avoids too many random components (e.g., mutation or crossover operators), which other methods, such as GA, rely on [27]. It has been reported that SS has been successfully applied in a variety of continuous and discrete optimization problems [27,32,33,34]. One prominent feature of a scatter search lies in its flexible structure, where a variety of ways and degrees of sophistication can be deployed for its elements [27]. Moreover, instead of operating on massive populations, SS concentrates only on a small set denoted herein as R e f S e t [32,33,34]. It obtains an optimum through balancing diversification (robustness) and intensification (efficiency) for R e f S e t systematically. The whole process of basic SS is shown in Figure 1.
In general, SS contains five strategies in total [27]: (i) a diversification generation method is used for generating an original population P , (ii) an improvement method wherein local search procedures can be embedded to examine the trial individuals and acquire high quality solutions, (iii) a reference set update method to initialize and maintain R e f S e t , normally resulting in a high quality subset R e f S e t 1 with best fitness value and a diversified subset R e f S e t 2 with optimal diverse value ( R e f S e t = R e f S e t 1 + R e f S e t 2 ), (iv) a subset generation method operates on R e f S e t , via which subsets are created and then utilized for producing combined vectors, and (v) a solution combination method for computing new testing vectors (one or more) from the results obtained from the subset generation method.
Optionally, a restart phase can be scheduled to rebuild the solutions. However, the SS framework described above is just a generic structure and all the five methods will be clearly defined in the latter designed ISS method.

3.2. The Utilized Local Search Solvers

As previously introduced, different SQP techniques, including filter-SQP and MISQP, are chosen as local solvers to complete different tasks. The process of the filter-SQP, as shown in Figure 2, essentially follows the procedure described in [28]. At each iteration, the sub-problem is replaced and solved as the trust-region sub-problem T R Q P ( x k , Δ k ) , with Δ k denoting the trust-region radius. For trust-region methods, however, a common problem lies in the difficulty of tuning penalty coefficients. A filter technique is hence incorporated to determine whether the solution outperforms the previous one and also guarantee the acceptability. The acceptance of the basic filter is defined by comparing the violation function value h ( x k ) and fitness value f ( x k ) :
h ( x k ) β h F o r f ( x k ) f F λ h F for all ( h F , f F ) F ,
where 0 < λ < β < 1 , F represents the filter set, h F and f F are the respective violation and fitness values of the element in the set. Additionally, with the aim of improving feasibility and optimality, a novel non-monotone technique for setting the filter criteria is also implemented [28]. In essence, it sets the filter by dividing the trail step d k into a quasi-normal part d k n and a tangential part d k t .
For CD-C and C-C charging modes, the problem turns out to be far more complex, since no explicit available criterion can be utilized to approach the optimal solution, nor can the optimal positions of the integers be effectively captured from relevant corresponding continuous solutions. The MISQP technique in [29] is therefore employed for solving the MINLP model in the proposed framework. Three key techniques are employed to achieve the tasks and they are: (i) a trust region technique with second-order amendments to stabilize the algorithm, (ii) a quasi-Newton formula to update the Hessian matrix, and (iii) a branch-and-bound skill to handle the sub-problems.

3.3. The Proposed ISS Algorithm for Single EV Charging

The implementation of the proposed ISS algorithm follows the generic framework of SS as shown in Figure 1. Several techniques have been introduced to enhance the overall performance.
(1) Initialization and restart mechanism
A diversification generation method is implemented with Latin hypercube uniform sampling [33]. The initial results are then filtered as the original population P , of which half ( R e f S e t 1 ) are top-ranking solutions with better fitness values, while another half ( R e f S e t 2 ) are arbitrarily selected from the rest of the diverse vectors. On top of this, a restart mechanism that aims to escape from local minima is also designed when improvement cannot be obtained after the subset updating step. For this process, while leaving the subset R e f S e t 1 unchanged, additional solutions are first created with the diversification generation method, and the subset R e f S e t 2 is reconstructed in the manner of subsequent updating method.
(2) Improvement procedure
The improvement procedure embedded with local search solvers is utilized to improve the results obtained from the diversification or the combination methods. Since different objectives and price types are considered for both NLP and MINLP problems, an intelligent module is designed to conduct the optimizations, as shown in Figure 3. The data will be categorized into different types, including: power rate type, charging type, price type selection, and objective selection, with which the charging power rate, charging type, price category and objective function would be determined respectively. All this information is then collected by the joint units where different local solvers are allocated.
(3) Population updating
For the update of the solutions, a new solution can only enter into R e f S e t 1 upon the following criteria: (i) the individual occupies a fitness value superior to the worst one in R e f S e t 1 , or (ii) the individual has a diverse value larger than the worst one in R e f S e t 2 . The diverse value is assessed according to the Euclidean distance d ( x , z ) between the candidate and the individuals in R e f S e t 1 . The new solution that maximize d m i n ( x ) will be selected to update R e f S e t 2 :
d m i n ( x ) = m i n z R e f S e t 1 { d ( x , z ) } .
(4) Subset generation and combination
A common strategy considering all pairs of the individual solutions is adopted to generate the basis for creating combined solutions. Instead of using simple linear combination, an improved hyper-rectangle based skill proposed in [34] has been chosen for this task. For each pairwise member, a bias of their relative position is utilized to define the hyper-rectangles for leading the solutions to approach “better” ones and move away from “bad” ones. This method can help to extend the search region and also expand the exploring directions without an additional memory term. Meanwhile, the “go-beyond” strategy from the same paper is also employed to enhance the intensification. The principle of this technique is to search the potential directions of the best generated individual and its parents continuously, in order to produce new individuals within two generations. Since the search area in this process is extended, the strategy promotes a diversity of solutions and accelerates convergence.
To further confine each element x k r of the solution vector within the correct search region, the repair strategy given below is employed:
if x k r > u b r , then x k r = u b r if x k r < l b r , then x k r = l b r
where l b r and u b r are respective lower and upper limits.

3.4. The Proposed GA-ISS Method for Massive EV Charging

When massive EV charging is simultaneously considered, the computation burden is enormous. Neither traditional deterministic techniques nor intelligent meta-heuristics can easily produce the optimized result, since: (i) the CPU may run out of memory due to large data demand, and (ii) long computation time is expected because of its complexity. For instance, a whole day is required to solve a single EV charging case that comes with 96 variables. Accordingly, there are 96 box constraints for Equation (1); 1 equality constraint for Equation (4), and 96 × (2 + 1) inequality constraints for Equations (2) and (5). Therefore, it is desirable to have an efficient algorithm to deliver accurate solutions.
To handle large-scale EV charging, a hybrid method integrating GA and ISS is proposed, as shown in Figure 4.
In Figure 4, N 0 , R 1 and R 2 denote the population scale, crossover, and mutation rate, respectively. The proposed method (GA-ISS) expects to benefit from the gist of both algorithms, for which GA processes present a strong global search ability and robustness, while the ISS structure is prominent in searching capability and quick convergence. During the procedure, a single EV state is determined by ISS, whilst the entire fleet co-evolves to the optimum and relies on the mutation and crossover procedures from GA.
To speed up the process, the initialization of population is first processed by ISS and copied into the N 0 populations, which will then be optimized in the main loop. The mutation operation takes the speed and optimization advantage of the ISS, through which individuals from the previous population are refined. As previous individuals are utilized as the starting points of the ISS, the time spent for refining the individuals will be substantially reduced, since the old ones have already been optimized to approach their optimal positions.

4. Simulation and Results

4.1. Parameter Setting

A case of overnight charging in a typical residential area where EV charging occurs during the night is considered in the simulations. The TOU price modeled from [35] is set on a rolling hourly basis, and the price curve as well as the tested basic daily load profile is shown in Figure 5. The LP rate coefficients are chosen to be ψ = 2 × 10 4 and γ = 0.22 . In addition, the maximum generation capacity is limited to be 20% above the daily peak load.
A single EV is used to illustrate the effectiveness of the ISS. The battery parameters and charging behavior are given in Table 2. It is also assumed that the minimum allowed SOC of the battery is 20%.
To test the GA-ISS method, an EV fleet with 100 vehicles is simulated. As EVs arrive and depart at the charging points in a random manner, the arriving time, departure time, rated capacity, rated plugging power, and the initial SOC of the test fleet are defined based on truncated normal distribution [36]. In Table 3, the standard deviation (Std.), minimum, maximum, and mean values of the EV fleet are listed. It should be noted that the parameters of EVs in this work are just adopted for demonstration purposes, and there is no restriction on the power and capacity of EV batteries for the proposed algorithms.

4.2. Simulation Results

The algorithm is coded on the MATLAB Platform (R2012a) and executed on a PC with 2.66-GHz Quad CPU (Intel Core 2), 4GB RAM and Windows 7 64-bit system. The performance of the ISS is first tested with the single EV, followed by the testing for a GA-ISS with an EV fleet.
(1) Single EV charging
For a single EV case, the algorithm can help the user determine the charging schedule to reduce the charging cost, or to better cooperate with the system operator to smooth the load profile. In order to evaluate the performance of ISS:
(a)
Two other commonly used heuristic algorithms, GA [37] and PSO [38], have been tested for flexible charging rate modes. And the GS method is utilized to produce the best reference result, since it explores all possible solutions [39].
(b)
Since a GS is unable to solve MINLP problems, only a GA and PSO [40] have been tested for constant charging rate.
(c)
The population number of the ISS is set to 30, while its iteration number of local solvers is limited to 10. The population member of the GA and PSO are set to 100. In addition, the total iteration numbers for all the above methods are set to 100.
Three different objective functions are used for the evaluation: (i) OF1—load curve flattening, refer to Equation (7), (ii) OF2—cost minimization with TOU price, refer to Equation (9), and (iii) OF3—cost minimization with LP price, refer to Equation (11).
For each condition, 50 simulations have been carried out, and Table 4 gives the mean values of the simulation results for all methods for comparison. It can be seen that the mean values for all the methods under different charging scenarios appear to be similar, however, an ISS can give slightly better performances compared with the heuristic GA and PSO algorithms for all three objective functions.
On the other hand, it is also observed that (i) a flexible charging rate enables more benefits than a constant charging rate, and (ii) the discharging ability exhibits more profits than charging-only cases in reducing charging costs and flattening daily load.
For more comparison to show the properties of the ISS, a box-plot for CD_C and CD_F mode with the OF1 objective is illustrated in Figure 6. From the distribution of the optimized outcome, it demonstrates that an ISS can achieve more consistent results than that of the other two methods. The robustness and good convergence of this novel algorithm are thus demonstrated.
Table 5 gives the average computational time of the respective algorithms for performing 50 simulations. It is clearly shown that the computational time of ISS is significantly shorter than that of the GS, GA, and PSO, which is the main advantage of the proposed algorithm. Furthermore, it can also be observed that the computational time is generally longer for (i) bidirectional charging cases as compared to the charging-only cases, and (ii) a flexible charging rate as compared to constant charging rate.
To examine the computational efficiency of the ISS scheduling algorithm, the iteration behaviors of GA, PSO and ISS, with trial simulations for the three objective functions under CD-F mode and CD-C mode are shown in Figure 7 and Figure 8. The convergence behaviors of C-F mode and C-C mode present similar properties respective to that of the CD-F and CD-C modes, hence their iteration figures are omitted here. As demonstrated in Figure 7 and Figure 8, despite the strong global search ability, the convergence process of both GA and PSO methods are very slow and may easily get trapped into local optima, whereas ISS always converges with far fewer iterations.
In terms of robustness, convergence, and efficiency, the proposed ISS method has adapted very well to single EV charging cases.
(2) Group EV charging
For massive amounts of EVs connecting to a charging station or regulated by an aggregator, the proposed GA-ISS hybrid method can help smooth the load profile or minimize the overall charging cost. For illustration purpose, CD-F mode with the aforementioned three objectives is utilized as the testing scenario for the group EV charging case. The desired SOC for all EVs is set to 90%. The maximum running time of the ISS at the initial phase is set as 0.8 s for an individual EV. For comparison and highlighting the benefit of scheduling for massive EVs:
  • The dumbing control method (DCM), which charges all EVs as soon as they are plugged in, is tested;
  • A GA-PSO hybrid method from [41] is chosen for comparison, since it can obtain better solutions along with less variation and processing time in comparison to other common heuristic methods. The population size and iteration number for GA-PSO and GA-ISS are set to be 20 and 50, respectively. Moreover, during the inner optimization for EV scheduling of GA-PSO and GA-ISS, the determination of a single EV state via PSO/ISS is set to be stopped when the minimum criteria of the solution quality is satisfied. In Table 6, the detailed parameters of the algorithms utilized for group EV charging are displayed.
  • A commercial CVX [42] toolbox is also simulated to perform the EV charging scheduling method, since it can simultaneously calculate all the variables and constraints and obtain the global optimized result. This approach is herein named global control.
Figure 9 demonstrates the simulation load profiles of these approaches with different optimization goals, chosen from one simulation result. As can be seen, the dumbing method can amplify the peak, indicating the necessity of EV charging scheduling. On the other hand, it is also shown that the proposed GS-ISS method obtains far better results than that of the GA-PSO method, and it can deliver a performance comparable to global control solved through the CVX solver, which requires a commercial licensing payment.
In addition, the final load profile of TOU price (OF2) is more fluctuated when compared to the LP price case (OF3), as shown in Figure 9b,c. As a matter of fact, the LP price load curve is almost the same as the OF1 one, which is obtained with the aim of load curve flattening. This suggests that an appropriate pricing strategy can help smooth the daily load curve with EV charging/discharging.
The average computational time as well as the optimization results for 20 simulations are listed in Table 7. It can be seen that the time spent by the hybrid GA-ISS approach is comparable to that of the global control, but outperforms the GA-PSO method. This verifies the effectiveness of the GA-ISS method for large scale EV scheduling. However, it should also be noted that the CVX solver may become stalled when more EVs are involved (e.g., when the EV number of the fleet changes to be 200, the computer may run out of memory with the tested PC), and hence solution may not be guaranteed for the global control.

5. Conclusion and Future Work

This paper presents the detailed derivation of two algorithms: an ISS framework to deal with single EV charging, and a hybrid GA-ISS method comprised of GA theory and the proposed ISS approach for handling a large-scale EV fleet. The main contribution is that the proposed algorithms can support both V2G and G2V, with either a flexible or constant charging power rate. It has been demonstrated that the ISS algorithm is the most computationally efficient way to obtain attractive performances among other methods including GS, GA, and PSO for a single EV charging scenario. On the other hand, for group EV charging, the GA-ISS approach has also been shown to be an extremely effective approach.
The proposed ISS algorithm and GA-ISS hybrid method are shown to be promising in terms of both efficiency and accuracy, thus providing potential techniques for the implementation of EV charging controls. The proposed algorithms are suitable to be implemented to schedule EV charging in many occasions, such as home charging, EV aggregations (stations), at power companies, etc. It should also be noted that this paper is focused on the issue of EV charging scheduling. Currently, energy storage is increasingly applied in power grids. The charging powers of energy storages are normally much higher than that of EVs, thus can exhibit a greater impact than EVs. Many characteristics of energy storage are similar to EVs and it is valuable to upgrade the algorithm for a wider usage, including the issue of energy storage.
Future work can be performed to incorporate more considerations for EVs, such as the various charging rates of a single EV battery, advanced modeling for SOC impacts on the EV charging rate, EV battery degradations, etc. Studies on the essential coordination between EVs and other power system components (renewable energy resources, home appliances, etc.) can be also addressed to accomplish more comprehensive and practical EV charging algorithms.

Author Contributions

T.M. and X.Z. were responsible to the conceptualization of the algorithms; T.M. conducted the investigation, designed the methodology and completed the original draft; X.Z. and B.Z. gave suggestions on the proposed idea and simulation cases; T.M. and X.Z. revised the paper based on comments from the reviewers.

Funding

This paper was supported by the project “Study on intelligent simulation technologies and panoramic simulation platform design for power market” of China Southern Power Grid (Project number: ZBKJXM20170082).

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

ST Time interval set
SE EV set
i Time index
j EV index
n Total number of timeslots
P E V , j i Charging power value for EV j in time i
P E V , j min Lower charging power limit for EV j
P E V , j max Upper charging power limit for EV j
S o c j , i SOC for EV j in time i
S o c j , s t a r t SOC value at the beginning for EV j
S o c j , end SOC value at the end for EV j
S o c j , min Minimum allowed SOC value for EV j
S o c j , max Maximum allowed SOC value for EV j
E exp , j Expectation level for the SOC of EV j
E R C , j Nominal battery SOC of EV j
P t o t a l i Total load demand in time i
P l d i Active base load in time i
P g i Maximum supply power in time i
C e i Electricity tariff in time i
C T O U i TOU price at time interval i
C L P i Linear price rate at time interval i
ψ Linear term of linear price
γ Constant term of linear price
C E V , j Nominal capacity of EV j

References

  1. Haidar, A.M.A.; Muttaqi, K.M.; Sutanto, D. Technical challenges for electric power industries due to grid-integrated electric vehicles in low voltage distributions: A review. Energy Convers. Manag. 2014, 86, 689–700. [Google Scholar] [CrossRef]
  2. De Quevedo, P.M.; Muñoz-Delgado, G.; Contreras, J. Impact of electric vehicles on the expansion planning of distribution systems considering renewable energy, storage and charging stations. IEEE Trans. Smart Grid 2017. [Google Scholar] [CrossRef]
  3. Morais, H.; Sousa, T.; Vale, Z.; Faria, P. Evaluation of the electric vehicle impact in the power demand curve in a smart grid environment. Energy Convers. Manag. 2014, 82, 268–282. [Google Scholar] [CrossRef] [Green Version]
  4. Tang, D.; Wang, P. Nodal impact assessment and alleviation of moving electric vehicle loads: From traffic flow to power flow. IEEE Trans. Power Syst. 2016, 31, 4231–4242. [Google Scholar] [CrossRef]
  5. Beaude, O.; Lasaulce, S.; Hennebel, M.; Mohand-Kaci, I. Reducing the impact of EV charging operations on the distribution network. IEEE Trans. Smart Grid 2016, 7, 2666–2679. [Google Scholar] [CrossRef]
  6. Zhang, X.; Zhong, Q.; Kadirkamanathan, V.; He, J.; Huang, J. Source-side Series-virtual-impedance Control to Improve the Cascaded System Stability and the Dynamic Performance of Its Source Converter. IEEE Trans. Power Electron. 2018. [Google Scholar] [CrossRef]
  7. Jian, L.; Zhu, X.; Shao, Z.; Niu, S.; Chan, C.C. A scenario of vehicle-to-grid implementation and its double-layer optimal charging strategy for minimizing load variance within regional smart grids. Energy Convers. Manag. 2014, 78, 508–517. [Google Scholar] [CrossRef]
  8. Xie, D.; Chu, H.; Gu, C.; Li, F.; Zhang, Y. A novel dispatching control strategy for EVs intelligent integrated stations. IEEE Trans. Smart Grid 2017, 8, 802–811. [Google Scholar] [CrossRef]
  9. Panwar, L.K.; Reddy, K.S.; Kumar, R.; Panigrahi, B.K.; Vyas, S. Strategic Energy Management (SEM) in a micro grid with modern grid interactive electric vehicle. Energy Convers. Manag. 2015, 106, 41–52. [Google Scholar] [CrossRef]
  10. Sun, B.; Huang, Z.; Tan, X.; Tsang, D.H.K. Optimal scheduling for electric vehicle charging with discrete charging levels in distribution grid. IEEE Trans. Smart Grid 2018, 9, 624–634. [Google Scholar] [CrossRef]
  11. García-Álvarez, J.; González, M.A.; Vela, C.R. Metaheuristics for solving a real-world electric vehicle charging scheduling problem. Appl. Soft Comput. J. 2018, 65, 292–306. [Google Scholar] [CrossRef]
  12. Mohamed, A.; Salehi, V.; Ma, T.; Mohammed, O. Real-time energy management algorithm for plug-in hybrid electric vehicle charging parks involving sustainable energy. IEEE Trans. Sustain. Energy 2014, 5, 577–586. [Google Scholar] [CrossRef]
  13. Deilami, S.; Masoum, S.A.; Moses, S.P.; Masoum, A.S.M. Real-time coordination of plug-in electric vehicle charging in smart grids to minimize power losses and improve voltage profile. IEEE Trans. Smart Grid 2011, 2, 436–467. [Google Scholar] [CrossRef]
  14. Wang, B.; Yang, J. Optimal electric vehicle charging scheduling with time varying profits. In Proceedings of the 2018 52nd Annual Conference on Information Sciences and Systems (CISS), Princeton, NJ, USA, 21–23 March 2018; pp. 1–6. [Google Scholar]
  15. Mao, T.; Lau, W.H.; Shum, C.; Chung, H.S.H.; Tsang, K.F.; Tse, N.C.F. A regulation policy of EV discharging price for demand scheduling. IEEE Trans. Power Syst. 2018, 33, 1275–1288. [Google Scholar] [CrossRef]
  16. Zhang, X.; Zhong, Q. Improved Adaptive-Series-Virtual-Impedance Control Incorporating Minimum Ripple Point Tracking for Load Converters in DC Systems. IEEE Trans. Power Electron. 2016, 31, 8088–8095. [Google Scholar] [CrossRef]
  17. Bourass, A.; Cherkaoui, S.; Khoukhi, L. Secure optimal itinerary planning for electric vehicles in the smart grid. IEEE Trans. Ind. Inform. 2017, 13, 3236–3245. [Google Scholar] [CrossRef]
  18. Zhao, T.; Peng, Y.; Nehorai, A. An optimal and distributed demand response strategy with electric vehicles in the smart grid. IEEE Trans. Smart Grid 2014, 5, 861–869. [Google Scholar]
  19. Leterme, W.; Ruelens, F.; Claessens, B.; Belmans, R. A flexible stochastic optimization method for wind power balancing with PHEVs. IEEE Trans. Smart Grid 2014, 5, 1238–1245. [Google Scholar] [CrossRef]
  20. Malikopoulos, A.A. Supervisory power management control algorithms for hybrid electric vehicles: A survey. IEEE Trans. Intell. Transp. Syst. 2014, 15, 1869–1885. [Google Scholar] [CrossRef]
  21. Soares, J.; Morais, H.; Sousa, T.; Vale, Z.; Faria, P. Day-ahead resource scheduling including demand response for electric vehicles. IEEE Trans. Smart Grid 2013, 4, 596–605. [Google Scholar] [CrossRef]
  22. Yu, J.J.Q.; Li, V.O.K.; Lam, A.Y.S. Optimal V2G scheduling of electric vehicles and unit commitment using chemical reaction optimization. In Proceedings of the IEEE Cong. Evolutionary Computation (CEC), Cancun, Mexico, 20–23 June 2013; pp. 392–399. [Google Scholar]
  23. Kang, Q.; Wang, J.; Zhou, M.; Ammari, A.C. Centralized charging strategy and scheduling algorithm for electric vehicles under a battery swapping scenario. IEEE Trans. Intell. Transp. Syst. 2016, 17, 659–669. [Google Scholar] [CrossRef]
  24. Alonso, M.; Amaris, H.; Germain, J.G.; Galan, J.M. Optimal charging scheduling of electric vehicles in smart grids by heuristic algorithms. Energies 2014, 7, 2449–2475. [Google Scholar] [CrossRef]
  25. Yilmaz, M.; Krein, P.T. Review of the impact of vehicle-to-grid technologies on distribution systems and utility interfaces. IEEE Trans. Power Electron. 2013, 28, 5673–5689. [Google Scholar] [CrossRef]
  26. Fang, X.; Misra, S.; Xue, G.; Yang, D. Smart grid—The new and improved power grid: A survey. IEEE Commun. Surv. Tutor. 2012, 14, 944–980. [Google Scholar] [CrossRef]
  27. Glover, F.; Laguna, M.; Martí, R. Fundamentals of scatter search and path relinking. Control Cybern. 2000, 29, 653–684. [Google Scholar]
  28. Yu, K.; Pu, D. A nonmonotone filter trust region method for nonlinear constrained optimization. J. Comput. Appl. Math. 2009, 233, 230–239. [Google Scholar]
  29. Exler, O.; Schittkowski, K. A trust region SQP algorithm for mixed-integer nonlinear programming. Optim. Lett. 2007, 1, 269–280. [Google Scholar] [CrossRef]
  30. Charging Lithium-Ion. Available online: http://batteryuniversity.com (accessed on 25 December 2018).
  31. Tsui, K.M.; Chan, S.C. Demand response optimization for smart home scheduling under real-time pricing. IEEE Trans. Smart Grid 2012, 3, 1812–1821. [Google Scholar] [CrossRef]
  32. Nebro, A.J.; Luna, F.; Alba, E.; Dorronsoro, B.; Durillo, J.J.; Beham, A. AbYSS: Adapting scatter search to multiobjective optimization. IEEE Trans. Evol. Comput. 2008, 12, 439–457. [Google Scholar] [CrossRef]
  33. McKay, M.D.; Beckman, R.J.; Conover, W.J. A comparison of three methods for selecting values of input variables in the analysis of output from a computer code. Technometrics 2000, 42, 55–61. [Google Scholar] [CrossRef]
  34. Egea, J.A.; Balsa-Canto, E.; Garcia, M.-S.G.; Banga, J.R. Dynamic optimization of nonlinear processes with an enhanced scatter search method. Ind. Eng. Chem. Res. 2009, 48, 4388–4401. [Google Scholar] [CrossRef]
  35. Qian, K.; Zhou, C.; Allan, M.; Yuan, Y. Modeling of load demand due to EV battery charging in distribution Systems. IEEE Trans. Power Syst. 2011, 26, 802–810. [Google Scholar] [CrossRef]
  36. Vagropoulos, S.I.; Bakirtzis, A.G. Optimal bidding strategy for electric vehicle aggregators in electricity markets. IEEE Trans. Power Syst. 2013, 28, 4031–4041. [Google Scholar] [CrossRef]
  37. Goldberg, D.E. Genetic Algorithms in Search, Optimization, and Machine Learning; Addison-Wesley: Reading, MA, USA, 1989. [Google Scholar]
  38. Kennedy, J.; Eberhart, R. Particle swarm optimization. In Proceedings of the IEEE International Conference on Neual Networks, IV, Perth, Australia, 27 November–1 December 1995; pp. 1942–1948. [Google Scholar]
  39. Kearfott, R.B. Rigorous Global Search: Continuous Problems; Kluwer: Dordrecht, The Netherlands, 1996. [Google Scholar]
  40. Kennedy, J.; Eberhart, R. A discrete binary version of the particle swarm algorithm. In Proceedings of the IEEE International Conference on Systems, Man, and Cybernetics, Orlando, FL, USA, 12–15 October 1997; pp. 4104–4108. [Google Scholar]
  41. Sheikhalishahi, M.; Ebrahimipour, V.; Shiri, H.; Zaman, H.; Jeihoonian, M. A hybrid GA-PSO approach for reliability optimization in redundancy allocation problem. Int. J. Adv. Manuf. Technol. 2013, 68, 317–338. [Google Scholar] [CrossRef]
  42. Michael, G.; Stephen, B. CVX: Matlab Software for Disciplined Convex Programming. Available online: http://cvxr.com/cvx (accessed on 25 December 2018).
Figure 1. Schematic of a standard scatter search (SS) template.
Figure 1. Schematic of a standard scatter search (SS) template.
Energies 12 00265 g001
Figure 2. Pseudo code of filter- sequential quadratic programming (SQP).
Figure 2. Pseudo code of filter- sequential quadratic programming (SQP).
Energies 12 00265 g002
Figure 3. The intelligent module for executing the local solvers of the intelligent scatter search (ISS).
Figure 3. The intelligent module for executing the local solvers of the intelligent scatter search (ISS).
Energies 12 00265 g003
Figure 4. A genetic algorithm (GA)-ISS based hybrid method.
Figure 4. A genetic algorithm (GA)-ISS based hybrid method.
Energies 12 00265 g004
Figure 5. Data profile for simulation.
Figure 5. Data profile for simulation.
Energies 12 00265 g005
Figure 6. Box-plot for CD_C and CD_F mode with objective OF1.
Figure 6. Box-plot for CD_C and CD_F mode with objective OF1.
Energies 12 00265 g006
Figure 7. Iteration map for single EV with different objectives under CD-F mode. (Left) Load standard deviation with OF1; (Center) Energy cost with OF2; (Right) Energy cost with OF3.
Figure 7. Iteration map for single EV with different objectives under CD-F mode. (Left) Load standard deviation with OF1; (Center) Energy cost with OF2; (Right) Energy cost with OF3.
Energies 12 00265 g007
Figure 8. Iteration map for single EV with different objectives under CD-C mode. (Left) Load standard deviation with OF1; (Center) Energy cost with OF2; (Right) Energy cost with OF3.
Figure 8. Iteration map for single EV with different objectives under CD-C mode. (Left) Load standard deviation with OF1; (Center) Energy cost with OF2; (Right) Energy cost with OF3.
Energies 12 00265 g008
Figure 9. Load profiles for group EV charging of DCM, GA-PSO, Global Control, and GA-ISS under CD-F mode.
Figure 9. Load profiles for group EV charging of DCM, GA-PSO, Global Control, and GA-ISS under CD-F mode.
Energies 12 00265 g009
Table 1. Comparisons for the four charging modes.
Table 1. Comparisons for the four charging modes.
Charging ModesCharging Rate Range (Normalized Value)Variables TypeSub-Problem
CD-F[−1,1]ContinuousNLP
CD-C[−1,1]DiscreteMINLP
C-F[0,1]ContinuousNLP
C-C[0,1]DiscreteMINLP
Table 2. Parameters for the single electric vehicle (EV).
Table 2. Parameters for the single electric vehicle (EV).
CapacityInitial SOCRated Charging PowerStart TimeEnd TimeDesired SOC
24 (kWh)58.59%3.3 kW21:156:3090%
Table 3. Parameters for the fleet with 100 EVs.
Table 3. Parameters for the fleet with 100 EVs.
ParameterMinimumMaximumMeanSt. Dev
Arriving Time (h)18:0022:0020:001:30
Departure Time (h)5:457:457:000:45
SOC (%)20905020
Capacity (kWh)1030186.93
Plugging power (kW)2103.541.48
Table 4. Mean values of the methods for single EV charging.
Table 4. Mean values of the methods for single EV charging.
ObjectivesTypeGSGAPSOProposed ISS
OF1
(standard deviation, kW)
CD-F368.19368.91368.30368.21
CD-CN/A368.40368.40368.36
C-F368.71368.96368.77368.74
C-CN/A369.00368.80368.80
OF2
(normalized value, p. u.)
CD-F44,277.1944,278.5544,278.0444,277.22
CD-CN/A44,278.6344,278.1944,277.69
C-F44,280.2344,281.5944,281.3444,280.45
C-CN/A44,280.6544,280.4244,280.39
OF3
(normalized value, p. u.)
CD-F40,457.3540,459.6640,459.1540,457.39
CD-CN/A40,459.7440,458.9840,458.40
C-F40,463.2140,464.2940,463.7240,463.24
C-CN/A40,463.5540,463.5440,463.43
Table 5. Computational time for a single EV for the methods under four charging modes (In seconds).
Table 5. Computational time for a single EV for the methods under four charging modes (In seconds).
Objective FunctionTypeGSGAPSOProposed ISS
OF1CD-F93.1611.678.511.18
CD-CN/A9.386.240.98
C-F49.3912.077.431.22
C-CN/A9.975.980.86
OF2CD-F80.3911.136.881.14
CD-CN/A9.776.150.95
C-F21.4112.317.321.03
C-CN/A9.626.470.81
OF3CD-F43.6112.656.421.25
CD-CN/A9.385.970.94
C-F31.4313.816.261.12
C-CN/A8.885.830.83
Table 6. Parameters of the algorithms utilized for group EV charging.
Table 6. Parameters of the algorithms utilized for group EV charging.
ParameterDCMGA-PSOGlobal ControlGA-ISS
Mutation rateN/A0.2N/A0.2
Crossover rateN/A0.4N/A0.4
Population sizeN/A20N/A20
Iteration numberN/A50N/A50
Table 7. Outcomes of different solving methods.
Table 7. Outcomes of different solving methods.
MethodOF1OF2OF3
ResultTimeResultTimeResultTime
GA-PSO298.02337.644,615.2291.241,274.5359.7
Global Control289.3992.644,538.187.741,222.294.2
GA-ISS291.75112.944,558.5105.941,229.2117.1

Share and Cite

MDPI and ACS Style

Mao, T.; Zhang, X.; Zhou, B. Intelligent Energy Management Algorithms for EV-charging Scheduling with Consideration of Multiple EV Charging Modes. Energies 2019, 12, 265. https://doi.org/10.3390/en12020265

AMA Style

Mao T, Zhang X, Zhou B. Intelligent Energy Management Algorithms for EV-charging Scheduling with Consideration of Multiple EV Charging Modes. Energies. 2019; 12(2):265. https://doi.org/10.3390/en12020265

Chicago/Turabian Style

Mao, Tian, Xin Zhang, and Baorong Zhou. 2019. "Intelligent Energy Management Algorithms for EV-charging Scheduling with Consideration of Multiple EV Charging Modes" Energies 12, no. 2: 265. https://doi.org/10.3390/en12020265

APA Style

Mao, T., Zhang, X., & Zhou, B. (2019). Intelligent Energy Management Algorithms for EV-charging Scheduling with Consideration of Multiple EV Charging Modes. Energies, 12(2), 265. https://doi.org/10.3390/en12020265

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop