BioInspired Learning and Adaptation for Optimization and Control of Complex Systems
View this Special IssueResearch Article  Open Access
Yimin Zhou, Zhifei Li, Xinyu Wu, "The Multiobjective Based LargeScale Electric Vehicle Charging Behaviours Analysis", Complexity, vol. 2018, Article ID 1968435, 16 pages, 2018. https://doi.org/10.1155/2018/1968435
The Multiobjective Based LargeScale Electric Vehicle Charging Behaviours Analysis
Abstract
In the paper, the effect of the charging behaviours of electric vehicles (EVs) on the grid load is discussed. The residential traveling historical data of EVs are analyzed and fitted to predict their probability distribution, so that the models of the traveling patterns can be established. A nonlinear stochastic programming model with the maximized comprehensive index is developed to analyze the charging schemes, and a heuristic searching algorithm is used for the optimal parameters configuration. With the comparison of the evaluation criteria, the multiobjective strategy is more appropriate than the singleobjective strategy for the charging, i.e., electricity price. Furthermore, considering the characteristics of the normal batteries and charging piles, user behaviour and EV scale, a Monte Carlo simulation process is designed to simulate the largescale EVs traveling behaviours in longterm periods. The obtained simulation results can provide prediction for the analysis of the energy demand growth tendency of the future EVs regulation. As a precedent of opensource simulation system, this paper provides a standalone strategy and architecture to regulate the EV charging behaviours without the unified monitoring or management of the grid.
1. Introduction
With the gradually deteriorated air quality, environments and energy crisis caused by the fuelpowered vehicles around the world, renewable vehicles, i.e., electric vehicles (EVs) are greatly promoted by all the governments, and many policies have been issued related to their development, where EVs will become the main transportation tools in the future along with the increasing improved technologies and infrastructure construction. However, the continuous energy demand growth of EVs puts a heavy pressure on the power grid with serious uncertainties in grid regulation [1]. Generally, one EV connection to the grid is roughly equal to the load requirement from one small household, where the construction of the grid in California, the most popular states of EVs in America, is facing the grid update problem currently [2]. The effect of EVs on the distributed power grid has been investigated in [3–5], and the involved impact factors can be summarized as traveling patterns, battery characteristics, charging schedule, and EV penetration, as shown in Figure 1.
From the analysis of the Danish national transportation survey data, the EV traveling model can be established with the driving distance and driving periods as the statistical data, so that the power demand and expected charging time from EVs can be determined [6]. The data demonstrate that the average driving distance of the EVs is 29.48 km, while the daily driving distance of 75% of the EVs is less than 40 km. Therefore, 40 km driving distance can be used to determine whether the EV battery capacity can meet the daily traveling requirement. Moreover, EVs can be considered as the mobile load connected between the power system and the transportation system, where the random mobility and charging of the EVs are dependent on the stochastic traveling patterns of the EV owners [7]. The graph theory was used to analyze the traffic and power network based on the driving patterns, and the randomly mobilized characteristics of the EV charging load was modelled accordingly.
Besides the requirement of the power quality (i.e., voltage and frequency) during EV charging, the required charging power is the most concerned for the power supply utilities [8]. According to the standard J1772 of the Society of Automotive Engineers [9], charging power can be divided into three levels, where grade 1 (1.5–3 kW) and grade 2 (10–20 kW) are suitable for family use, and grade 3 (40 kW and above) is more suitable for EV fast charging, by which EV charging can be completed within one hour; however, it is mostly used in centralized charging scenarios due to overly high power.
The charging plan is normally determined by two factors, charging availability (charging situation) and charging strategy [10]. Since the energy consumption of the daily EV journey is less than the batteryrated capacity, it is unnecessary to charge the EV on daily basis. The EV can start to be charged only when its battery stateofcharge (SOC) is less than the threshold. The EV charging strategies play an important role in the evaluation of the EV influence on the power system [11], which can be divided into (1) simple charging (dumb charging), i.e., the unplanned “plug and play” pattern, generally charged at the end of the trip in one day or when the charging facility is available; (2) tariffdriven charging, i.e., charging during offpeak periods with cheaper cost; and (3) intelligent charging, i.e., the extra battery energy can be used as the energy resource to provide assistant service during highpeak hours, which is beneficial to the stable grid operation.
Another important aspect to analyze the EVs influence on the grid is to design the evaluation index quantitatively. In [12], based on the insular power system of Crete in Greece, it studied the impact of plugin EVs with smart and direct charging patterns on the grid energy scheduling and cost. An EV aggregator that participates in the market is proposed with the pricetaking approach for solving optimally selfscheduling problem. Compared with different scenarios of EV penetration, an irregular performance would be presented on the unit commitment schedule while the number of EVs reaches 21,000. Paper [13] introduces an architecture for the intelligent energy management system, which is composed of an admission control, a pricing module, and a power scheduling module that determines the charging sequence of EVs. The queuing problem of power dispatch is formulated as a stochastic dynamic programming problem, and a threshold admission and greedy scheduling method is adopted to minimize the electricity cost.
Based on the review of EVs impacts on the grid, it can be shown that EVs popularization will bring nonneglectful new electric load. At present, however, most research conclusions about largescale EVs charging are drawn from simulation experiments. There are obvious discrepancies of gas station facilities and EV supporting facilities among different cities. How to improve the credibility and universality of the prediction results is still a challenge. Therefore, many involved issues should be addressed to improve the load prediction accuracy, i.e., driving and charging models with appropriately practical condition, traveling probability distribution model, and multiobjective EV charging patterns.
Furthermore, most works merely consider the price incentive in the charging modes, and the relation between EVs traveling behaviours and charging schemes has not been accurately described by mathematical models [14, 15]. This paper investigates largescale EV charging behaviours as a multiobjective optimization problem, and the main contributions are summarized as follows: (1)From the benefits of the EV users and power suppliers perspectives, the comprehensive evaluation index system has been developed with three key factors: load peak value, charging bills, and traveling rate, which are used as the objective to describe the prediction and optimization problems of the EVs power usage requirement(2)Based on the actual traffic data, the probability distribution model of the traveling patterns is established; Monte Carlo (MC) is adopted to establish the largescale EVs traveling and charging simulation system to improve the load prediction accuracy(3)A multiobjective charging strategy has been developed under multiple constraints, and the involved parameters are determined via a novel constraint mechanism which transfers the equalityboundary constraints into binary feature selection(4)Besides, the nonlinear model of the charging behaviour and the MC simulation algorithm have been fully opensource [16]. Based on the proposed framework shown in Figure 2, more gridimpact assessments can be focused on the situation of EVs in different districts and their charging facilities
The remainder of the paper is organized as follows. Section 2 analyzes the energy usage status of the EVs, and the optimized indexes are determined. The statistical model of the traveling pattern and the charging strategies of the EVs are introduced in Section 3. Monte Carlo simulation algorithm is used to calculate the solution of the nonlinear programming model in the same section. In Section 4, a multiobjective charging strategy is developed with optimal parameters, and the proposed strategies are compared via simulation experiments as well. Conclusions and future work are given in Section 5.
2. Problem Formulation
2.1. System Description
According to the traffic survey in Beijing [17], the traveling periods of private EVs during working days are concentrated among 6:00–9:00 (early peak) and 16:00–19:00 (evening peak), and 40.3% of the parking time (slack time >5 hours) of EVs are distributed among 18:00–21:00, which are depicted in Figure 3. If EVs charging is performed without any guidance, peak load would be generated due to the largescale charging concentrated at parking periods. According to the 2020 Shenzhen EVs growth plan [18], the current power grid will be difficult to satisfy the energy demand of EVs, which would result in power supply tension, or even endanger the safety of the grid operation [19]. Furthermore, more investment should be put to satisfy higher load peak requirement, resulting in largescaled equipment upgrade and capital budget pressure. On the other hand, due to the constrained battery capacity and limited charging piles, the EVs traveling plan would be hindered if without appropriate charging schedule, which is currently one of the most important reasons to limit the popularity of EVs.
Therefore, an objective to study the effect of largescaled EV charging on the grid load (hereinafter referred to as “EV charging problem”) is summarized as: based on the routine traveling habits of EV users, the optimal charging strategy is proposed to reduce the grid peak load, lower the cost of EV users, and guarantee the success of traveling plan, which can be ascribed as one of the NPhard problems.
First, the relationship between the traveling patterns and charging behaviour are analyzed so as to extract the involved key factors and determine the evaluation indices. The total required power of the EV charging from the grid system at the specific moment is determined by the four main impact factors, i.e., battery characteristics , charging piles , users behaviour (), and the current EV number . As it is shown, the characteristics of the battery pack, charging piles, and current EV number can be assumed as known variables. The variables describing user behaviours are unknown variables which will be determined by mathematical models.
For this type of NPhard problems which cannot be solved directly, the MC simulation method is adopted to simulate the EVs operation system under different charging strategies for the purpose of optimal EV charging mode obtainment. Based on the statistics, different probability distributions are selected to describe the EVs driving behaviours to improve the credibility of the MC simulation results. Finally, several commonly used optimization algorithms are introduced to adjust the parameters configuration of the charging modes so as to realize the optimal comprehensive index.
According to the technical specifications of the small EVs and common characteristics of urban residents’ vehicle traveling habit, several assumptions have been made to describe the complex EVs charging problem without loss of generality. (1)The EV battery and parameters related to the charging piles are set similar to those of the commercial products(2)Each day (24 hours) will be divided equally into 96 intervals, i.e.,15 minutes as the sampling period(3)The traveling frequency, departureparking time, and driving mileage are the main factors directly affecting the EV charging behaviour, i.e., the charging period (), initial SOC (), and the required charging quantity (). So the traveling patterns and the charging strategies are the two key factors to be focused(4)In order to guarantee the success of the expected traveling plan, it is assumed that the battery packs will be fully charged at each charging(5)The traveling patterns of the personal EVs on the working days are considered; other types of EVs or driving patterns under holidays/weekends will not be addressed unless otherwise specified(6)The urban residents in Shenzhen and in Beijing share the same traveling pattern
Replacing the input module (shown in Figure 2) with specific traffic data, the problem formulation and models can be extended to other urban areas or holiday/weekend driving patterns. Thus, the designed framework has the advantage of solving a universal problem and benefits to obtain reasonable results in different situations.
2.2. The Evaluation Index System
The different evaluation indices are listed in Nomenclature after Conclusions, which are developed from the perspectives of power suppliers and EV owners . Since the averagepeak ratio (), saving rate (), and traveling rate () are the most concerned indices for the energy suppliers and EV owners [10, 20], they are used to form the comprehensive index , which is defined as where , , and are the positive weight coefficients (), and they can be set flexibly. Here, , , and , denoting load balance index taken as priority. The higher value of means higher system performance. is used to normalize the index value by equation (2), and are bound of . The normalization can improve the sensitivity of the optimization algorithm to the index variation so as to avoid premature convergence.
3. The Ev Load Estimation Based on Monte Carlo Simulation
After the optimal objective and key impact factors are determined, the EV charging problem can be configured as a nonlinear programming model. Under the constraints of the EV battery characteristics and user traveling patterns, the objective functions can be established based on different demands from the power suppliers and EV users so as to obtain the optimal charging schemes. The undetermined variables in the model, i.e., , are quantified via MC simulation tool to simulate the daily EV traveling, charging, and slack state. Furthermore, the statistic models of the traveling habit and EV charging behaviour have been established to improve the MC simulation credibility level.
3.1. The Probability Distributed Fitting Model for the Traveling Variables
According to the data collected from the GPS installed on 112 private EVs in Beijing [17], totally, 4892 data from June 2012 to March 2013 are recorded and used as the relevant information are outlined in Table 1. Here, from the analysis of the EVs driving data in working days, the probabilistic density function (PDF) and its parameters of the traveling variables can be determined as follows (see Figure 4): (1)To select the PDF from the generally used probabilistic functions, i.e., exponential, gamma, normal, and Poisson distribution functions, which are called “PDFx”(2)To estimate the scale and shape parameters of the “PDFx” through the maximum likelihood estimation (MLE) [21](3)To verify whether the generated 1000 groups of data accept the original hypothesis at a confidence level of 95% (significance level ) by the KS, Ftest, and Ttest such statistical methods [22](4)According to the value results, the fitting degree between the “PDFx” model is obtained(5)To evaluate the best fitting of the “PDF” as “PDFbest” after iterative comparison

Then the obtained “PDFbest” of each variable is applied to the simulated system so the input values satisfied specific probabilistic distributions are in accordance with the actual vehicles traveling patterns. Although there are certain biases existence at and (shown as Figures 5(a) and 5(d)), different hypotheses testing methods (including chisquare test, KolmogorovSmirnov test, and ttest) are adopted and verified that these deviations are acceptable with sufficient statistical samples, the detailed results are achievable in our opensource project [16]. Therefore, the designed MC simulation model is quite reliable with high credibility. (1)The daily traveling frequency, , follows distribution, where is the gamma function, is the random variable, is the PDF of , and and are the shape parameter and scale parameter, respectively. Through the MLE, , the expectation , and the variance can be obtained. The comparison between the actual data and the fitted curve of the traveling frequency is shown in Figure 5(a)(2)The driving mileage of each traveling, , follows BirnbaumSaunders (BS) distribution, where is the PDF of the BS distribution; and are the scale parameter and shape parameter, respectively. Through the MLE, , , , and can be obtained. The probability density distribution curves of the actual data and the BS fitting curve are shown in Figure 5(b)(3)The duration of each traveling, , follows distribution. Suppose each traveling duration also obeys gamma distribution as expressed in equation (3), the parameters can be acquired via MLE, i.e., , , , and . The fitted curve with distribution and the actual data curve are depicted in Figure 5(c). (4a)The departure time of each traveling (AM), , follows a locationscale distribution, where is the PDF of , is the gamma function, and are the location parameter, scale parameter, and shape parameter, respectively. Through the MLE, , , , and can be obtained.(4b)The departure time of each traveling (PM), , follows normal distribution, Through the MLE, the expectation and the standard variance . The fitted curve with probability density distribution and the actual data curve are shown in Figure 5(d).
(a) The distribution of the daily traveling frequency Fd
(b) The distribution of the driving mileage of each traveling Md
(c) The distribution of the duration of each traveling Td
(d) The distribution of the departure time of each traveling (AM/PM) and
3.2. Modelling the Charging Strategies
EV users would select different charging strategies, i.e., to set a specific target or to charge at optimal periods. For instance, under the incentive of the timeofuse (ToU) electricity prices [23], the users would charge EVs during parking periods with lower electricity prices, which is called the tariffguiding strategy. The usually used three charging strategies are described in Table 2.

Equation (7) describes the unified quantization formula of the users’ charging motivation. By changing the weight coefficient of the three charging strategies, the charging strategy models can be derived quantitatively as follows: where is the charging priority to be arranged at certain period, i.e., EV starts charging when ; is the random priority factor following the uniform distribution with , in order to simulate the random charging behaviour; and are the predicted time intervals required for the fully charged and the slack time intervals, respectively; denotes that if the current SOC of the EV, , the EV would be charged immediately. In order to ensure the reliability of the traveling plan, . and are the average charging cost at certain periods and the average charging cost of one day, which can be calculated as follows: where is the continuous charging duration and is the ToU price at the moment, as listed in Table 3; here, .

Through the comparison of the priority level at each interval during EV slack periods, the charging moment can be regarded as the highest priority (e.g., the maximum random value, the least cost, or the longest parking time) to achieve the maximum benefit. The pseudocode of the charging procedure is shown in Algorithm 1.

3.3. Monte Carlo Simulation
Based on the EV traveling and charging statistics, the MC is used for numerical simulation to solve the output performance index of the different charging schemes.
3.3.1. The Basic Concept of the MC Method
MC method is known as a stochastic simulation method or statistical testing method to use random sampling for mathematical function estimation. MC has statistical convergence where the fitting deviation converges to a certain threshold.
The predicted model of the EV charging capacity calculation is developed based on MC random sampling tests. The capacity of the battery charging is calculated per day evenly divided into 96 intervals, and the total charging capacity in the time interval is described as where represents the charging capacity of the EV in the time slot on the workday, represents the number of the EVs acquiring power supply from the grid in the time period, and is the total counted days. Several constraints are set in the simulation. Consider the condition that the battery should be in fully charged status before driving, the starting time of the battery charging should be limited as where and are the starting and ending slack status of the EV, is the maximum continuous charging duration, is the charging capacity, is the initial state of the battery capacity that is related to the driving distance per trip, and delegates the starting charging moment under fully charged conditions. The condition of convergence adopted in MC sampling model is expressed as where is the variance coefficient of the system indicator at the moment; , , and are the variance, expectation, and standard deviation. The repeated times in MC simulation is at least 100, and the variance coefficient is set to less than 0.5%.
3.3.2. The Comprehensive Index Calculation via MC Simulation
Firstly, according to the users traveling model, the daily driving duration and mileage are determined, while the EV charging time period is dependent on the traveling status, characteristics of the battery, and the adopted charging strategy. After the determination of the EV driving/charging period, the total driving mileage, charging power, and charging bills can be calculated in the 96 time intervals per day. Finally, the comprehensive index and related variance coefficients can be studied based on equations (1) and (12) to complete the whole day simulation.
According to the Shenzhen EVs demonstration and promotion plan, the number of private EVs will reach 240,000 in 2020, which will be used for the charging load calculation. Based on the probability distribution model of the traveling pattern and three types of charging strategies, the simulated largescale EVs operation procedure via MC simulation is shown in Figure 6. The related indices of the predicted grid load per day in 2020 are shown in Table 4.

It can be seen from Table 4 that there are large difference for the peak loads or with different charging strategies. Compared to the traveling habits of EV users, the EVs charging behaviours are more controllable via certain guidance. Hence, the charging strategy is one of the most feasible optimization objectives.
4. Multiobjective Charging Strategies for the LargeScale Evs
According to the government report, the grid load of the Shenzhen exceeded 15598.6 MW at 11:10 AM of 27th June 2016, and the EV total charging load will exceed 5000 MW per day in 2020. Moreover, with random and parkingcharging patterns, EV owners select to perform charging during peak periods, i.e., 9 AM–14 PM and 20–23 PM, which is 1.57% and 1.97% of the historical peak value, respectively. From the grid safety perspective, it could easily cause peak superposition of the electricity usages, thus putting heavy pressure on the grid accommodation capacity obviously.
On the other hand, EVs charging with tariff guidance can avoid the charging during highpeak periods but could generate another charging peak value 824.64 MW in shorter periods between 23:00–2:00, which is shown in Figure 7. Thus, with only the adoption of ToU pricing mechanism, it would result in other concentrated charging periods and abruptly steep highpeak electricity usage would endanger the stability of the grid operation seriously.
From the comprehensive index perspective, the three charging strategies cannot provide satisfied performance under single factor consideration. For instance, the charging at arrival pattern can guarantee the success of traveling but also can produce higher bills; moreover, charging too frequently would be harmful to the battery life. If only the pricingbased charging mode is adopted, the tariff could be least during midnight; however, in addition to a sharp peak would be caused in short period of time, the required electricity could not be guaranteed if the emergency traveling is required.
Here, considering the impact factors involved in the system indicator in equation (1), the three single indices and the compressive index are configured as the objective set, which is denoted by . Then a combined multiobjective charging strategy can be proposed with the developed nonlinear programming model, which is defined as where are the coefficients as described in equation (7); directly impacts the required energy supply and the traveling distance of next trip. Given that EV starts charging when , so is added as the model bias to make the charging decision threshold more reasonable. These parameters are set to be optimized.
For the multiobjective problems, previous studies suggest the gradientbased solution [24, 25]. As described in [25], when the objective function is nondifferentiable, gradientbased optimization techniques are not suitable for dynamic programming models with constraints. Moreover, there are many potential local extrema in the charging model, and the gradient descent algorithm can easily lead to pseudoconvergence. Considering that our problem has constraints of equality and inequality, the objective function is not required to be continuous and differentiable in heuristic algorithm; the multiobjective optimization algorithm based on feature selection is proposed to find the global optimal combination of system parameters.
First, the objective function is established as ; three most popular heuristic algorithms, i.e., particle swarm optimization (PSO) [26], genetic algorithm (GA) [27], and simulated annealing (SA) [28] are selected to search the nondominated set. Then the maximized comprehensive index, , can be extracted from the Paretofront space, so that a unique optimal solution can be obtained via the proposed algorithm.
4.1. Optimal Searching Algorithm under Constraints
According to the constraints expressed in equation (13), a feature selectionbased algorithm is proposed to solve the problem from continuous optimization transformed into discrete optimization. First, the involved variables are coded into binary form. Then the feature selection is applied to optimize the feature code of each individual so as to improve the quality of the variables combination under constraints. The related opensource algorithms can be referred in [16], and the detailed steps are described as follows.
Step 1. Preprocess each input variables, , , and with 8bit binary coded data, as Algorithm 2
Step 2. Initialize the particles swarm; each is expressed by an 8bit 5 (variables)
Step 3. With the adoption of feature selection algorithms, e.g., BPSO, optimize the feature code of each individual so as to improve the quality of the variable combination
Step 4. Perform the inverse procedure of Step 1 and transform the binary code of each individual into variables set form in order to calculate the individual fitness. The pseudocode (Algorithm 3) are displayed to describe the equation constraint
Step 5. Implement the MC procedure to acquire the fitness, and the converged negative index () is regarded as the individual fitness. The related parameters are set to speed up the calculation and guarantee convergence as , , and
Step 6. If the optimality condition is satisfied or the maximum number of iterations has been reached, the procedure stops and the optimal variable combination solution is obtained; otherwise, go to Step 3 to continue the next iteration
Compared to the general constrainthandling methods, such as penalty function [29], the proposed method possesses simplicity and high efficiency by transforming the continuous combination problem into a binary feature optimal problem. It can avoid invalid searching due to a variable out of bound in each iteration. Therefore, binary feature optimization can be applied in multiple intelligent optimal algorithms.
The PSO, GA, and SA are the typical algorithms for the solution of multiobjective problems with constraints. Combined with the basic feature selection method, they are used to solve the objective in equation (13) with multiconstraint. The involved parameter settings are listed in Table 5.


(a)  
 
(b)  
 
(c)  

4.2. Experimental Results and Analysis
In order to investigate the access impact of the largescale EVs to the grid load features, the simulations are kept running for 3 days with two cluster servers, whose hardware are listed in Table 6. The EV charging behaviours are analyzed with multiple charging strategies, and the performance index difference with multiobjective schemes are also discussed.

4.2.1. The Convergence Analysis of the MC Process
The results of the EVs charging in 100 days under different charging strategies via MC simulation are shown in Figure 8. The variance coefficient of the comprehensive value (calculated by equation (13)) can be converged lower than 0.5%, where the reliability is satisfied with the MC simulation requirement.
Based on the traveling model with certain credibility and different charging modes, MC is used to simulate the 240,000 EVs driving in 100 working days. After the convergence of the simulated system is achieved, the daily average charging load curves can be acquired with four different strategies, as shown in Figures 9 and 10(a).
(a) The hourly charging loads from the EVs
(b) The hourly charging bills (left Yaxis) and the saving rates (right Yaxis)
(c) The hourly aborted trips (left Yaxis) and the traveling rates (right Yaxis)
4.2.2. The Performance Comparison of the Charging Strategies
The index comparison of different charging schemes is illustrated in Figure 10. First, under the random strategy, there will be new loads added as a superposition during the peak periods of daytime, and the scheduled trips of EV could be cancelled with higher probability due to insufficient capacity (see Figure 10(c)). Moreover, a sharp peak load could be generated in the midnight with the tariffguiding strategy. As for the parkingcharging strategy, a peak load as well as the highest electricity payment during the evening period could be generated (see Figure 10(b)).
However, the multiobjective charging strategy can achieve a compromise among different indicators (see Figure 9). The essence of the multiobjective is to balance the resource allocation among indicators, and the selection of weights in the formula is rather subjective. Therefore, heuristic algorithm is used to determine the Pareto optimum of variables.
Table 7 summarizes the indexes of different charging strategies during daily driving. The is the highest with random charging mode; the charging tariff is the least with tariffguiding mode; the highest rate of success traveling can be achieved via parkingcharging model; the best comprehensive index is achieved via the multiobjective charging strategy. The results demonstrate that the combination model described in equation (13) and the parameters configuration obtained via SA can guide the charging behaviour effectively and improve the performance indices greatly.

4.2.3. Observations and Discussion
The experimental results shown in Figure 10(a) indicate that the power load will reach 245.01 MW during peak periods in 2020 with 240,000 EVs charging without any guidance in Shenzhen. The power suppliers have to take measures to update the current grid infrastructure and develop the smart grid to reduce the obviously peakvalley difference. Thus, based on the feature selection, three optimal searching algorithms considering constraints are proposed and applied for solving multiobjective charging problem. As listed in Table 7, under the multiobjective strategy and its optimized parameters configuration, the is increased by 11%, charging cost is saved by 66.2%, and the successful traveling rate is 99.5%. Therefore, the multiobjective combination of charging strategy can provide a maximized comprehensive benefit for the EV users and power suppliers. From the users perspective, it is suggested that the EV should be charged timely if its is less than .
Differentiating from other literatures that concern about the accuracy verification of the load prediction [30, 31], this paper designs a universal model to determine the optimal periods for standalone charging behaviours, which has no direct communication between the grid and EVs [11, 13] or the online adjustment on the electricity price [12]. Referring to the multiobjective strategy, the onboard programmable control system can estimate the optimal periods to charge the EV automatically during EV slack periods set by the EV owners. Not only the traveling plan can be guaranteed but also the electricity bills saving and offpeak energy usage can be achieved simultaneously.
5. Conclusions
This paper investigates the EVs largescale charging behaviours and discusses intelligent charging strategies. First, the probability distribution model of the traveling pattern based on the actual traffic data is established to improve the credibility of the prediction. Considering the benefits for the EV users and power suppliers, the evaluation indices are selected, including load peak value, charging bills, and traveling rate. Under the constraints of the EV battery features and users traveling statistics characteristics, the multiobjective charging strategy has been developed, where the Pareto optimum of the parameters are determined via a novel constrainthandling method with feature selection. Monte Carlo tool is adopted to simulate the EVs activities, and an opensource system with a universal framework is released to promote the research on access impact of the largescale EVs to the local grid load.
In summary, the adopted probability distribution statistics, MC simulation, evaluation indices, and EV charging model can provide an alternative for the discussion of largescale EVs development and grid resources assignment. Further research will investigate the MCbased modelling for different EVs (i.e., public transportation) in holidays and the EV connection mode in the realworld grid network. Besides, more work should be focused on the statistical analysis of EVs access to distributed grid nodes, which can give insights on how to realize the proposed methods in the local distribution grids.
Nomenclature
Battery:  Battery capacity, 100 Ah 
:  Voltage, 230 V 
:  Charging efficiency, 0.9 
:  Full charged duration, 5 h 
:  Energy consumption, 0.175 kW/h/km 
:  Stateofcharge = remainder capacity/rated capacity. 
:  Demand power from grid, constant power, 15 kW. 
:  Traveling frequency (times) per day 
:  The driving mileage of each traveling, unit: km 
:  The duration of each traveling 
:  The departure time of each traveling 
:  Charging periods 
:  The SOC at the initial charging moment 
:  The required charging power, unit: kW/h . 
:  The amount of EVs. 
:  Daily charging peak value 
:  Total load/time interval 
:  
:  Load rate = 
:  Charging cost = electricity price × power consumption 
:  The highest price charging quantity 
:  Saving rate = 
:  Actually driving mileage/planned mileage 
:  The comprehensive index; the higher value means higher system performance. 
:  Three single indices and the comprehensive index make up the objective set 
:  The coefficients of the factors (, i.e., randomness, price guidance, and slack time) 
:  The lowest limit of the battery SOC 
:  The model bias. 
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The work was supported by the Shenzhen Science and Technology Innovation Commission Project Grant Ref. JCYJ20160510154736343 and Ref. JCYJ20170307165442023 and Science and Technology Planning Project of Guangdong Province with Ref. 2017B010117009.
References
 D. B. Richardson, “Electric vehicles and the electric grid: a review of modeling approaches, impacts, and renewable energy integration,” Renewable and Sustainable Energy Reviews, vol. 19, pp. 247–254, 2013. View at: Publisher Site  Google Scholar
 S. Eric, “The impact of growing electric vehicle adoption on electric utility grids,” Electric Utilities, 2017. View at: Google Scholar
 S. Habib, M. Kamran, and U. Rashid, “Impact analysis of vehicletogrid technology and charging strategies of electric vehicles on distribution networksa review,” Journal of Power Sources, vol. 277, pp. 205–214, 2015. View at: Publisher Site  Google Scholar
 Y. Zhou and X. Li, “The stateofart of the EV charging control strategies,” in 2015 34th Chinese Control Conference (CCC), pp. 7916–7921, Hangzhou, China, July 2015. View at: Publisher Site  Google Scholar
 R. GarciaValle and J. P. Lopes, Eds., Electric Vehicle Integration into Modern Power Networks, Springer Science & Business Media, New York, NY, USA, 2012. View at: Publisher Site
 Q. Wu, A. H. Nielsen, J. Ostergaard et al., “Driving pattern analysis for electric vehicle (EV) grid integration study,” in 2010 IEEE PES Innovative Smart Grid Technologies Conference Europe (ISGT Europe), pp. 1–6, Gothenberg, Sweden, October 2010. View at: Publisher Site  Google Scholar
 D. Tang and P. Wang, “Dynamic modeling of electric vehicle movable loads based on driving pattern analysis,” in 2014 IEEE PES General Meeting  Conference & Exposition, pp. 1–5, National Harbor, MD, USA, July 2014. View at: Publisher Site  Google Scholar
 I. K. A. Aswantara, K. S. Ko, and D. K. Sung, “A centralized EV charging scheme based on user satisfaction fairness and cost,” in 2013 IEEE Innovative Smart Grid TechnologiesAsia (ISGT Asia), pp. 1–4, Bangalore, India, November 2013. View at: Publisher Site  Google Scholar
 Society of Automotive Engineers and C. B. Toepfer, SAE Electric Vehicle Conductive Charge Coupler, SAE J 1772, REV. MONTH01 (DOC), California Air Resources Board, 2009.
 E. Veldman and R. A. Verzijlbergh, “Distribution grid impacts of smart electric vehicle charging from different perspectives,” IEEE Transactions on Smart Grid, vol. 6, no. 1, pp. 333–342, 2015. View at: Publisher Site  Google Scholar
 V. Aravinthan and W. Jewell, “Controlled electric vehicle charging for mitigating impacts on distribution assets,” IEEE Transactions on Smart Grid, vol. 6, no. 2, pp. 999–1009, 2015. View at: Publisher Site  Google Scholar
 S. I. Vagropoulos, G. A. Balaskas, and A. G. Bakirtzis, “An investigation of plugin electric vehicle charging impact on power systems scheduling and energy costs,” IEEE Transactions on Power Systems, vol. 32, no. 3, pp. 1902–1912, 2017. View at: Publisher Site  Google Scholar
 Z. Yu, S. Chen, and L. Tong, “An intelligent energy management system for largescale charging of electric vehicles,” CSEE Journal of Power and Energy Systems, vol. 2, no. 1, pp. 47–53, 2016. View at: Publisher Site  Google Scholar
 X. Fernandes, J. Rebelo, J. Gouveia, R. Maia, and N. Bustorff Silva, “Onoff scheduling schemes for powerconstrained electric vehicle charging,” 4OR, vol. 15, no. 2, pp. 163–181, 2017. View at: Publisher Site  Google Scholar
 R. Rao, H. Cai, and M. Xu, “Modeling electric taxis’ charging behavior using realworld data,” International Journal of Sustainable Transportation, vol. 12, no. 6, pp. 452–460, 2018. View at: Publisher Site  Google Scholar
 “The open source simulation of the largescale electric vehicle charging behaviours,” https://cn.mathworks.com/matlabcentral/fileexchange/62816. View at: Google Scholar
 H. Wang, X. Zhang, L. Wu et al., “Beijing passenger car travel survey: implications for alternative fuel vehicle deployment,” Mitigation and Adaptation Strategies for Global Change, vol. 20, no. 5, pp. 817–835, 2015. View at: Publisher Site  Google Scholar
 Implementation Plan of Demonstration and Extension of Energy Conservation and New Energy Vehicles in Shenzhen (2009–2012), Report from Shenzhen Development and Reform Commission, Shenzhen, China, 2010, in Chinese.
 J. Lei, W. Duan, X. Dong, and Z. Hu, “Plugin electric vehicles charging load calculating based on Monte Carlo simulation in Shenzhen,” in The 14th China Association for Science and Technology Annual Meeting, pp. 1–6, Hebei, China, 2012. View at: Google Scholar
 E. B. Iversen, J. M. Morales, and H. Madsen, “Optimal charging of an electric vehicle using a Markov decision process,” Applied Energy, vol. 123, pp. 1–12, 2014. View at: Publisher Site  Google Scholar
 S. Guindon and O. Gascuel, “A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood,” Systematic Biology, vol. 52, no. 5, pp. 696–704, 2003. View at: Publisher Site  Google Scholar
 E. L. Lehmann and J. P. Romano, Testing Statistical Hypotheses, Springer Science & Business Media, 2006.
 Notice on Publication of Tariff and Peak Tariff Tariffs for Residents, Price Bureau of Guangdong, Government Report, 2015, in Chinese.
 C. Liu, Z. Gong, K. L. Teo, J. Sun, and L. Caccetta, “Robust multiobjective optimal switching control arising in 1,3propanediol microbial fedbatch process,” Nonlinear Analysis: Hybrid Systems, vol. 25, pp. 1–20, 2017. View at: Publisher Site  Google Scholar
 C. Liu, Z. Gong, K. L. Teo, and E. Feng, “Multiobjective optimization of nonlinear switched timedelay systems in fedbatch process,” Applied Mathematical Modelling, vol. 40, no. 2324, pp. 10533–10548, 2016. View at: Publisher Site  Google Scholar
 Y. Zhou, Y. Chen, G. Xu, Q. Zhang, and L. Krundel, “Home energy management with PSO in smart grid,” in 2014 IEEE 23rd International Symposium on Industrial Electronics (ISIE), pp. 1666–1670, Istanbul, Turkey, June 2014. View at: Publisher Site  Google Scholar
 O. Babatunde, L. Armstrong, J. Leng, and D. Diepeveen, “A genetic algorithmbased feature selection,” British Journal of Mathematics & Computer Science, vol. 4, no. 21, pp. 889–905, 2014. View at: Google Scholar
 K. A. Dowsland and J. M. Thompson, “Simulated annealing,” in Handbook of Natural Computing, G. Rozenberg, T. Bäck, and J. N. Kok, Eds., pp. 1623–1655, Springer, Berlin Heidelberg, 2012. View at: Publisher Site  Google Scholar
 A. K. Yadav, R. Ranjan, U. Mahbub, and M. C. Rotkowitz, “New methods for handling binary constraints,” in 2016 54th Annual Allerton Conference on Communication, Control, and Computing (Allerton), pp. 1074–1080, Monticello, IL, USA, September 2016. View at: Publisher Site  Google Scholar
 R. C. Leou, J. H. Teng, and C. L. Su, “Modelling and verifying the load behaviour of electric vehicle charging stations based on field measurements,” IET Generation, Transmission & Distribution, vol. 9, no. 11, pp. 1112–1119, 2015. View at: Publisher Site  Google Scholar
 S. Shafiee, M. FotuhiFiruzabad, and M. Rastegar, “Investigating the impacts of plugin hybrid electric vehicles on power distribution systems,” IEEE Transactions on Smart Grid, vol. 4, no. 3, pp. 1351–1360, 2013. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Yimin Zhou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.