Evaluation of Accommodation Capability for Electric Vehicles of a Distribution System Considering Coordinated Charging Strategies

: In recent years, the ever-increasing charging demand of electric vehicles (EVs) imposes challenges on both power supply security and reliability in the distribution system. In this paper, an EV accommodation capability evaluation model of a distribution system, with high penetrations of ﬂexible resources, is established. Firstly, according to the actual classiﬁcations of EVs and transportation rules, a Monte Carlo simulation is used to simulate the charging behaviors of EVs so as to obtain the relevant parameters of EV charging. Then, a coordinated charging optimization model for various types of EVs is proposed based on the charging characteristics of EVs. The presented model comprises a mixed-integer linear programming problem and a constrained optimization problem which are respectively solved by CPLEX (the Simplex method implemented in the ◦ C programming language) and the particle swarm optimization (PSO) algorithm. Last of all, a real-life distribution system in the coastal areas of China is served for demonstrating the feasibility and e ﬃ ciency of the proposed approach. Moreover, the impacts of ﬂexible resources, distribution network zoning rules, and EV growth on the EV accommodation capability of a distribution system are also discussed.


Introduction
The rapid development of electric vehicles (EVs) has received extensive attention [1][2][3], due to their great potential for saving energy, and reducing emissions. The ever-increasing penetration of EV charging demand has been observed in many actual power distribution systems. As a new type of load (with a high degree of uncertainty), the integration of numerous EVs could bring great challenges to the secure and economic operations of the distribution system [4][5][6]. This may result in increasing energy losses, voltage deviations, and even overloads in feeders and/or substations. The relationship among feeder losses, load factor, and load variance is explored in the context of coordinated charging of plug-in hybrid vehicles (PHEV) in Reference [1]. In Reference [7], the impacts of EV charging on power loss and the voltage threshold-crossing of the distribution network at different scales are analyzed. The impact of EV integration on the life cycles of transformers is investigated in Reference [8]. The mathematical model for harmonic evaluation of multiple EV charging machines is established in Reference [9]. Therefore, it is necessary to conduct a quantitative evaluation of the accommodation 1.
The ACE and the OECL of a distribution system are comprehensively presented so as to consider the EV accommodation and charging load distribution holistically. Furthermore, EVs are divided into four major categories to discuss the corresponding OECL in detail.

2.
The proposed co-optimization model considers both the optimal dispatch of flexible resources and the optimal charging scheme for EVs with respect to security constraints, such as power balance, voltage limits and N − 1 security criterion.

3.
A data-driven method is proposed to model EV charging behaviors, which is important in coordinating charging schedule. The distribution system is divided into multiple subareas according to its load characteristics. Coordinated charging optimization is performed in each subarea to meet its reliability standard while addressing its load characteristics.
The remainder of this paper is organized as follows. The framework of the model is illustrated in Section 2. The ACE evaluation model and the coordinated charging model are proposed in Section 3. Case studies are given to demonstrate the proposed method and to examine the impact of different factors in Section 4. Finally, conclusions and future research are presented in Section 5.

Framework of the Model
Different types of EV propulsion technology like pure EVs and PHEVs/EREVs are available, and the design of various DC-DC converter topologies for Battery Electric Vehicles (BEVs) and Plug-in Hybrid Electric Vehicles (PHEVs) are reviewed in Reference [21]. An overall picture of the current EV technology and prospects of future development are provided in Reference [22]. According to the different technical standards of the charging piles employed at present, EV charging modes are broadly divided into a quick mode (e.g., DC charging mode) and slow mode (e.g., constant AC power charging mode). Specifically, rapid DC charging applies to electric buses, cars, taxies and trucks, while slow constant power charging is targeted for small-scale EVs which obtain charging service mostly in public parking lots, large shopping centers, community garages, and household charging piles.
On this basis, an ACE evaluation model, with respect to power supply security constraints of the distribution system penetrated with flexible resources, is established at first. A coordinated EV charging optimization method is then proposed, taking into account the uncertainties in EV charging behaviors. EV charging load patterns are optimized according to real-life EV charging rules and EV ownership, and thereby a comparison between the OECL and the ACE evaluation is conducted. A schematic diagram of the model framework is shown in Figure 1.

Formulation of the Proposed Model
In this section, we first give some reasonable assumptions: 1 All the charging loads of EVs at the charging stations or charging facilities that are connected to one bus are grouped into one load, and the ACE for each bus is evaluated in our model. 2 EVs are divided into the following four categories in this paper: Electric buses, taxies, private cars, and social vehicles (mainly referring to rental cars). Each type of vehicle is affiliated with a corresponding charging rate. 3 It is assumed that slow charging piles can interact with distribution system operators in real-time and the operators are capable of controlling the charging power and the duration for each EV once it is plugged into the system.

Formulation of the Proposed Model
In this section, we first give some reasonable assumptions: 1.
All the charging loads of EVs at the charging stations or charging facilities that are connected to one bus are grouped into one load, and the ACE for each bus is evaluated in our model.

2.
EVs are divided into the following four categories in this paper: Electric buses, taxies, private cars, and social vehicles (mainly referring to rental cars). Each type of vehicle is affiliated with a corresponding charging rate.

3.
It is assumed that slow charging piles can interact with distribution system operators in real-time and the operators are capable of controlling the charging power and the duration for each EV once it is plugged into the system.

Evaluation Model of the ACE in a Distribution System
In order to maintain a secure and economic operation of the distribution system and accommodate as many EVs as possible, the objective function is to simultaneously maximize the EV accommodation capability and minimize operation cost of the distribution system with respect to the economic operational constraints.
(1) The objective function The first objective of the model is to maximize the ACE, while the second objective is to minimize the operational costs of the distribution system which include the power purchase cost, the generation cost of DG units, the compensation for loads involved in the demand side response, and the ESS compensation cost. Then these two objective functions are normalized and integrated into one objective function based on the linear combination.
and O bj all are respectively the ACE, the operational costs and the aggregative objectives of the model; T is the number of time periods in a day-ahead dispatch scheme; P EV i,t is the EV charging power at node i in time period t; N EV is the number of nodes with EV charging facilities in the system; N node , N DG , N IL , N TL , and N ESS are the numbers of nodes, DG units, interruptible loads, transferable loads, and the energy storage devices in the distribution system, respectively; C S t , C DG g,t , C IL t , C TL t , and C ESS t are the electricity purchasing price from the upstream transmission system, the unit generation cost of DG unit g, the unit compensation price for interruptible loads, the unit compensation price for transferable loads, and the unit compensation cost of the energy storage device in time period t, respectively; P G i,t is the active power injected at node i in time period t, and if node i is the slack bus then the pertinent active power is injected by the upstream supply transformer, otherwise the injected active power at node i is set as 0; P DG g,t is the actual output power of DG unit g (a wind turbine unit or photovoltaic unit) in time period t; P IL d,t is the active power curtailment of interruptible load d in time period t; P TL r,t is the transferred active power of transferable load r in time period t; P ESSC m,t and P ESSD m,t represents the active charging and discharging power of energy storage device m in time period t, respectively; α c m,t and α d m,t respectively represent the charging and discharging states of energy storage device m in time period t, and both are binary variables subject to α c m,t + α d m,t ≤ 1; and w 1 , w 2 > 0 are the weights of the linear combination of the two normalized objective functions.
(2) Constraints (a) Power balance constraints The piecewise linearized DistFlow equations [19,23] are adopted in this paper to describe nodal power balances. where Ω E and Ω L represent the sets of nodes and branches in the distribution system, respectively; Ω i is the set of nodes directly connected to node i; P D i,t and Q D i,t are the active and reactive load power at node i in time period t; Q G i,t is the reactive power injected at node i in time period t; P ij,t and Q ij,t are respectively the active and reactive power on branch ij in time period t; P QU ij,t and Q QU ij,t are used to estimate the quadratic terms of P ij,t and Q ij,t , and both can be approximated using the piecewise linearization approximation (PLA) method [19]; r ij and x ij are respectively the resistance and reactance of branch ij; and V 0 is the voltage amplitude at the slack bus. The network topology will change if an N − 1 contingency occurs-and, consequently, the power balance constraints will change accordingly. The state variable λ ij,t is then introduced to represent the working state of branch ij in time period t; λ ij,t = 0 or 1 stands for that branch ij is out of operation or in operation.
(b) Branch and node voltage constraints where V i,t is the voltage amplitude at node i in time period t; Ω S represents the sets of the supply substation nodes; V G is defined as the voltage of the supply substation; V max i and V min i are respectively the upper and lower limits of the voltage amplitude at node i; M is a very large positive number.
When λ ij,t = 0, i.e., branch ij is out of operation, Equation (5) will never be binding, since M is very large; in other words, this constraint will always be met. When λ ij,t = 1, i.e., branch ij is in operation, the left-hand side of the inequality in Equation (5) must be zero, indicating that the voltage drop constraint for branch ij is linearized.

(c) Power supply margin constraints
The N − 1 security inspection is widely used in power system planning for reliability testing, which requires that the system not only effectively supplies the load in the normal operation, but also can guarantee the load supply in case of failure of any single component connected in the system. Besides, feeder or unit contingencies will affect the power supply of the power system, which in turn affects the ACE of the distribution system.
According to the differences of the power backup condition and the power supply reliability level in the wiring mode, N − 1 security inspection is carried out based on the corresponding wiring mode. The changes of feeder connection in case of N − 1 contingencies are illustrated with two specific wiring examples.
For a line which is connected to two or more power sources at the same time and has the ability to transfer power, as shown in Figure 2a, when N − 1 contingencies occur in region B, the connection at the interconnection switch will be reconnected and the residual loads of region B, originally brought by Substation 2, will be transferred to Substation 1 [24,25]. Then, the load rate of Substation 1 is increased, which limits the ACE of nodes in this region. In this case, the power supply margin constraints, when N − 1 contingencies occur, can be summarized as constraints (8): where Ω J and Ω F respectively represent the sets of nodes and feeders in wiring group l; and Ω A and Ω B are respectively the sets of nodes in regions A and B. However, for a radial connection or multi-segment wiring type, like the single-end power supply line shown in Figure 2b, when N − 1 contingencies occur to a certain section of region B, the line in this area will lose power supply and all the loads connected will be cut off. In this situation, constraints (8) can be transformed into Equation (9).
Energies 2019, 12, x FOR PEER REVIEW 6 of 21 (d) Nodal charging power constraints EV SME , where P SME i is the sum of the rated charging power of all charging piles connected to node i.

(e) Transformer capacity constraints
The active and reactive power constraints of a power supply transformer are as follows: where P where Pij,t is the active power of feeder ij in time period t; P LCM ij is the active power capacity limit of feeder ij.
(g) DG output constraints pv max , where P max w and P wind w,t are, respectively, the capacity and actual active power output of wind turbine w in time period t; P max p and P pv p,t are, respectively, the capacity and the actual active power output of PV unit p in time period t.
(h) DR constraints (d) Nodal charging power constraints where P SME i is the sum of the rated charging power of all charging piles connected to node i.

(e) Transformer capacity constraints
The active and reactive power constraints of a power supply transformer are as follows: where P Gmax k and Q Gmax k are respectively the upper limits of the active and reactive power of the transformer k.
(f) Feeder capacity constraints where P ij,t is the active power of feeder ij in time period t; P LCM ij is the active power capacity limit of feeder ij. where P max w and P wind w,t are, respectively, the capacity and actual active power output of wind turbine w in time period t; P max p and P pv p,t are, respectively, the capacity and the actual active power output of PV unit p in time period t.

(h) DR constraints
In this paper, the interruptible loads (ILs) and the transferable loads (TLs) are employed to implement DR. The parameters associated with the ILs include the duration and quantity of the interruptible loads and the summation of the TLs are kept constant within one scheduling period.
where δ IL t and δ TL t are respectively the maximum allowed proportion of the ILs and the TLs over the total loads in time period t; T d , T max   (20) and (21), which impose the upper and lower limits of the charged state of the energy storage devices.
where S m is the charged state of the energy storage device m; P max m is the capacity of the energy storage device m; S max m and S min m are the upper and lower limits of the charged state of device m, respectively.

The Coordinated Charging Optimization Model
Based on available data of public charging piles and special charging stations for EVs, in a Chinese coastal province, electric buses, taxies, and social vehicles usually prefer quick charging, and terminate charging after an expected level of the state of charge (SOC) is attained. Therefore, these kinds of EVs are not suitable for participating in charging scheduling optimization. Private EV cars usually adopt slow charging at office locations or at home. Parking time of private EV cars at charging piles is normally long, and is then suitable for participation in the coordinated charging optimization. Parameters of different types of EVs discussed in this paper are listed in Table 1. Specific implementation methods of the coordinated charging optimization model include the following steps: (1) EV ownership and charging rules are statistically obtained, and a probability model of charging behaviors is attained; (2) the well-established Monte Carlo simulation is used to generate the charging behavior parameters; (3) the coordinated charging optimization is finally obtained based on sampling results.

Statistical Analysis
Based on the statistical analysis of EV charging data, transportation rules of electric buses, taxies, social vehicles, and private cars are obtained by piecewise fitting using standard fitting methods. From the analysis of EV charging data, the charging duration of buses approximately follows a lognormal distribution, and its probability density function can be expressed as in Equation (22) where f bus s is the probability distribution function of the charging duration of buses; t bin is the charging duration of a bus; A bus is the function amplitude; μ bus and σ bus are respectively the mean and the standard deviation of the logarithmic charging duration.
The probability density function of the initial SOC of electric taxies is divided into two stages based on their charging start time; each approximately follows a lognormal distribution, as shown in Equations (23) and (24). Based on the distributions of the initial SOC, the charging duration can be estimated as in Equation (25). The probability distribution function of the initial SOC for private cars is consistent with Equation (24), due to similar transportation patterns. Based on the statistical analysis of EV charging data, transportation rules of electric buses, taxies, social vehicles, and private cars are obtained by piecewise fitting using standard fitting methods. Figures 3 and 4 show the sample distributions of the initial charging time of EVs for buses and social vehicles. The probability distributions of the charging durations for four types of EVs are approximated in Equations (22)(23)(24)(25)(26). From the analysis of EV charging data, the charging duration of buses approximately follows a lognormal distribution, and its probability density function can be expressed as in Equation (22) where f bus s is the probability distribution function of the charging duration of buses; t bin is the charging duration of a bus; A bus is the function amplitude; μ bus and σ bus are respectively the mean and the standard deviation of the logarithmic charging duration.
The probability density function of the initial SOC of electric taxies is divided into two stages based on their charging start time; each approximately follows a lognormal distribution, as shown in Equations (23) and (24). Based on the distributions of the initial SOC, the charging duration can be estimated as in Equation (25). The probability distribution function of the initial SOC for private cars is consistent with Equation (24), due to similar transportation patterns.  From the analysis of EV charging data, the charging duration of buses approximately follows a lognormal distribution, and its probability density function can be expressed as in Equation (22).
where f bus s is the probability distribution function of the charging duration of buses; t bin is the charging duration of a bus; A bus is the function amplitude; µ bus and σ bus are respectively the mean and the standard deviation of the logarithmic charging duration. The probability density function of the initial SOC of electric taxies is divided into two stages based on their charging start time; each approximately follows a lognormal distribution, as shown in Equations (23) and (24). Based on the distributions of the initial SOC, the charging duration can be estimated as in Equation (25). The probability distribution function of the initial SOC for private cars is consistent with Equation (24), due to similar transportation patterns. are respectively the mean and the standard deviation of the driving distance at different times; P taxi ch and P taxi bat are respectively the charging rate and the battery capacity of taxies. Charging rules of social vehicles are analyzed using data from public charging stations, as most social vehicles get access to charging services at public charging stations. The statistics show that the charging duration approximately follows a normal distribution shown in Equation (26).
where f social s is the probability distribution function of the charging duration of social vehicles; A social is the function amplitude; t sin is the charging duration of a social vehicle; µ social and σ social are respectively the mean and the standard deviation of the charging duration.
Combining transportation and charging data for different types of EVs, MATLAB's fitting toolbox is used to fit the probability distribution parameters. The statistical analysis of other charging behavior parameters can also be adopted in the same way.

Sampling
Monte Carlo simulation is used to sample the initial SOC, the charging duration, and the arrival time of EVs, with input data of EV ownership, performance parameters, and a probability distribution of charging behavior parameters. The sampling results of the initial SOC and arrival time of taxies are shown in Figure 5 as an example.

Sampling
Monte Carlo simulation is used to sample the initial SOC, the charging duration, and the arrival time of EVs, with input data of EV ownership, performance parameters, and a probability distribution of charging behavior parameters. The sampling results of the initial SOC and arrival time of taxies are shown in Figure 5 as an example.

Coordinated Charging Optimization
The objective function of the coordinated charging optimization model is to minimize the load's peak-valley difference within the system scheduling period.

Coordinated Charging Optimization
The objective function of the coordinated charging optimization model is to minimize the load's peak-valley difference within the system scheduling period. min P cut = P peak − P valley (27) where P cut is the peak-valley difference of the system load; P peak and P valley are the load peak and valley during the system scheduling period, respectively. The charging power of electric buses, taxies, and social vehicles are respectively formulated in Equations (28)-(30) based on sampling results.
where P bus i,t , P taxi i,t , and P social i,t refer to the total charging power outputs of buses, taxies, and social vehicles connected to node i within time period t, respectively; P bus ch and P social ch are the charging rates of buses and social vehicles, respectively; b bus v,t , b taxi v,t , and b social v,t are respectively the charging states of the vth bus, taxi, and social vehicle in time period t with values of 1 or 0 standing for the charging and the idle condition, respectively; N bus i , N taxi i , and N social i are the numbers of buses, taxies, and social vehicles connected to node i for charge. Equation (31) shows the initial charging time constraint of private cars, while (32) ensures that EVs be charged to the expected SOC before departure.
where t arr v , t real v , t dep v , and t last v are respectively the times when the vth EV arrives at the charging pile, the optimized initial charging time of the vth EV, departure time of the vth EV, and the charging duration. The charging power of private cars can be calculated by Equation (33). The OECL is then computed with the charging power of other EVs, which should not exceed the corresponding ACE.
where P car i,t refers to the charging power of private cars connected to node i within time period t; P car ch is the rated charging power of a private car; b car v,t is the charging state of the vth private car in time period t; N car i is the number of private cars connected to node i for charge; P REV i,t represents the OECL connected to node i within time period t.
The optimization procedure is shown in Figure 6.  Figure 6. Flowchart of the optimization procedure.

Relationship between the Two Models
The accommodation capability for EVs, calculated by the ACE evaluation model, is transferred to the coordinated charging optimization problem through P EV i,t . Then sampling for different types of EVs is conducted, and the OECL is obtained by the coordinated charging optimization model using the real-life EV charging and transportation rules. When there is a significant mismatch between the charging demand and the optimized charging power attained in the upper level, if the charging demand is lower than that in ACE it will be fully supplied-otherwise the charging demand will be partially supplied with the amount of the upper limit.

Relationship between the Two Models
The accommodation capability for EVs, calculated by the ACE evaluation model, is transferred to the coordinated charging optimization problem through P EV i,t . Then sampling for different types of EVs is conducted, and the OECL is obtained by the coordinated charging optimization model using the real-life EV charging and transportation rules. When there is a significant mismatch between the charging demand and the optimized charging power attained in the upper level, if the charging demand is lower than that in ACE it will be fully supplied-otherwise the charging demand will be partially supplied with the amount of the upper limit.

Solving Method
The ACE evaluation model is a mixed-integer linear programming problem. The YALMIP platform and solver CPLEX (the Simplex method implemented in the • C programming language) [26] are used to solve the problem. The Particle Swarm Optimization (PSO) algorithm [27,28] is then adopted to solve the coordinated charging optimization problem.

System Settings
In this section, an actual distribution system in a coastal area of China is employed to demonstrate the proposed method. There are seven 110 kV substations, 376 nodes, 368 branches, 6 PV units, and two wind farms, with 51 nodes connected to EV charging facilities in this system. The system is divided into eight power grid units according to the division rules mentioned in Section I. The numbers of EVs in different parts of the distribution system are listed in Table 2. The basic load profile of a typical day is selected to represent the system basic load. Based on wind speed and light intensity data in a certain time period, the forecasted outputs of the DGs can be obtained by employing the Latin hypercube sampling algorithm [29].

Comparisons of Results among Multiple Cases
In order to demonstrate the impact of flexible resources on the ACE, four different cases are investigated: Case 1: The original system without flexible resources; Case 2: An updated system with DG integration into the original system; Case 3: Case 2 with ESS integrated; Case 4: An integrated system with DG, ESS, and DR.
The accommodation capability for EVs by the distribution system under different cases and outputs from flexible resources are shown in Figure 7, where the positive value of energy storage represents the released energy, and the negative value represents the stored energy. By comparing the ACEs between Case 1 and Case 2 in Figure 7, it is clear that DG integration can effectively improve the ACE of the distribution system in the daytime as the majority of DG units are photovoltaic units that mainly generate electricity during the daytime. The accommodation level for the power output from DGs can also be enhanced in the low load period as EVs can be charged at night. In Case 3, the ESS stores the unconsumed energy of DGs during the valley load period, and releases it in the peak load period, so that the ACE is enhanced. Among the results of Cases 1-4, the ACE of the distribution system is the highest in Case 4, as the integration of multiple flexible resources increases the overall duration for accommodating EV loads, and hence enhances the ACE of the distribution system. The comparison of the four cases shows the ACE promotion effects by the integration of flexible resources. However, the addition of flexible resources will not yield a particularly large increase in the ACE, as the ACE is also affected by the transmission capacity of the power distribution system. Under the condition that the power grid structure remains unchanged, the load-bearing capacity of the system is limited.  Figure 8 shows the daily system load profile and the ACE time variation curve. It can be seen that the ACE is variable with time, and the ACE is limited by the peak load of the system. Figure 9 shows the spatial distribution of the ACE and charging demand in a typical time period. Please note that only the geographical locations of the nodes connected with EV charging facilities are shown in this figure, while the remaining nodes and feeders are not shown. The color of each circular node represents the strength degree of its ACE such that higher ACE values are represented by red nodes and lower ACE values are represented by yellow nodes. The density of green dots in the figure reflects the density of EV charging demands on the road, obtained by applying the multi-agent-based microscopic traffic assignment model [30]. As shown in Figure 9, there are some nodes in a grid, where the ACEs are adequate to satisfy the charging demands. However, there are still some cases where the ACEs do not match the charging demands. For example, in Grid-8, the ACEs at some nodes are sufficient, but the EV charging demands on the road around these nodes are relatively sparse, meaning that a part of the charging facilities would be idle. On the other hand, the charging demand in the central part of the region like Grid-7 is high, but the ACE nearby is relatively low. EV users have to get EVs charged at a relatively far distance, which may increase traffic congestion.   Figure 8 shows the daily system load profile and the ACE time variation curve. It can be seen that the ACE is variable with time, and the ACE is limited by the peak load of the system.  Figure 8 shows the daily system load profile and the ACE time variation curve. It can be seen that the ACE is variable with time, and the ACE is limited by the peak load of the system. Figure 9 shows the spatial distribution of the ACE and charging demand in a typical time period. Please note that only the geographical locations of the nodes connected with EV charging facilities are shown in this figure, while the remaining nodes and feeders are not shown. The color of each circular node represents the strength degree of its ACE such that higher ACE values are represented by red nodes and lower ACE values are represented by yellow nodes. The density of green dots in the figure reflects the density of EV charging demands on the road, obtained by applying the multi-agent-based microscopic traffic assignment model [30]. As shown in Figure 9, there are some nodes in a grid, where the ACEs are adequate to satisfy the charging demands. However, there are still some cases where the ACEs do not match the charging demands. For example, in Grid-8, the ACEs at some nodes are sufficient, but the EV charging demands on the road around these nodes are relatively sparse, meaning that a part of the charging facilities would be idle. On the other hand, the charging demand in the central part of the region like Grid-7 is high, but the ACE nearby is relatively low. EV users have to get EVs charged at a relatively far distance, which may increase traffic congestion.   Figure 9 shows the spatial distribution of the ACE and charging demand in a typical time period. Please note that only the geographical locations of the nodes connected with EV charging facilities are shown in this figure, while the remaining nodes and feeders are not shown. The color of each circular node represents the strength degree of its ACE such that higher ACE values are represented by red nodes and lower ACE values are represented by yellow nodes. The density of green dots in the figure reflects the density of EV charging demands on the road, obtained by applying the multi-agent-based microscopic traffic assignment model [30]. As shown in Figure 9, there are some nodes in a grid, where the ACEs are adequate to satisfy the charging demands. However, there are still some cases where the ACEs do not match the charging demands. For example, in Grid-8, the ACEs at some nodes are sufficient, but the EV charging demands on the road around these nodes are relatively sparse, meaning that a part of the charging facilities would be idle. On the other hand, the charging demand in the central part of the region like Grid-7 is high, but the ACE nearby is relatively low. EV users have to get EVs charged at a relatively far distance, which may increase traffic congestion. It should be mentioned that access to EVs is not limited to nodes with existing charging stations or scattered charging piles. All nodes in the distribution network are supposed to be able to accept EVs for charging, since the planning of new charging facilities is taken into account in this system. The results of the ACE of the whole distribution network at a typical time point in a day are shown in Figure 10. In this figure, each part of the distribution network is colored according to the ACE, with color strength proportional to the ACE. Comparing Figures 9 and 10, the difference between the ACE of nodes with charging facilities and the ACE of the system mainly lies in Grid-7. The ACE of nodes with charging facilities in Grid-7 is weak, while the overall ACE of the whole distribution system is the strongest among all parts of the system. This demonstrates that the distribution network structure corresponding to Grid-7 is relatively strong, which can accommodate the access of charging loads of a large number of EVs with secure operation maintained. However, due to the limited number of existing charging facilities and other reasons like city construction, the current ACE of Grid-7 is lower. Then Grid-7 is a good candidate region for new EV charging stations. The evaluation of the capability of the whole It should be mentioned that access to EVs is not limited to nodes with existing charging stations or scattered charging piles. All nodes in the distribution network are supposed to be able to accept EVs for charging, since the planning of new charging facilities is taken into account in this system. The results of the ACE of the whole distribution network at a typical time point in a day are shown in Figure 10. In this figure, each part of the distribution network is colored according to the ACE, with color strength proportional to the ACE. It should be mentioned that access to EVs is not limited to nodes with existing charging stations or scattered charging piles. All nodes in the distribution network are supposed to be able to accept EVs for charging, since the planning of new charging facilities is taken into account in this system. The results of the ACE of the whole distribution network at a typical time point in a day are shown in Figure 10. In this figure, each part of the distribution network is colored according to the ACE, with color strength proportional to the ACE. Comparing Figures 9 and 10, the difference between the ACE of nodes with charging facilities and the ACE of the system mainly lies in Grid-7. The ACE of nodes with charging facilities in Grid-7 is weak, while the overall ACE of the whole distribution system is the strongest among all parts of the system. This demonstrates that the distribution network structure corresponding to Grid-7 is relatively strong, which can accommodate the access of charging loads of a large number of EVs with secure operation maintained. However, due to the limited number of existing charging facilities and other reasons like city construction, the current ACE of Grid-7 is lower. Then Grid-7 is while the overall ACE of the whole distribution system is the strongest among all parts of the system. This demonstrates that the distribution network structure corresponding to Grid-7 is relatively strong, which can accommodate the access of charging loads of a large number of EVs with secure operation maintained. However, due to the limited number of existing charging facilities and other reasons like city construction, the current ACE of Grid-7 is lower. Then Grid-7 is a good candidate region for new EV charging stations. The evaluation of the capability of the whole distribution network can explore its potential of accommodating EVs and provide some guidance for the subsequent charging facility planning.

Spatiotemporal Distribution of the ACE
4.5. The Regional Distribution of OECL EV ownership and EV types all correspond to the basic loads, and load types in the power grids, so it is necessary to explore the regional OECL distribution characteristics.
As shown in Figure 11, the load curve trough time of each subarea is identical and is around 08:00. In connection with a real-life situation, this period is generally the departure time of EVs. Besides, the OECL in Grid 1 is the largest among all the grids concerned because of the high EV ownership, and there are two load peaks around 13:00 and 22:00. This is because this subarea is mainly composed of residence and office areas, which reveals a behavior of getting out early and coming back at dusk. Therefore, there is a charging peak after EV owners go home at night; meanwhile, the charging of social vehicles and buses near the office area during the day forms another charging peak. The load structure of Grid 6 is consistent with that of Grid 1, so it presents a similar charging load profile. Grid 3 is the downtown area and charging loads mainly arrive during the day, while EV charging loads in Grid 4 mainly appear at night as this grid is dominated by residential areas. The charging load profile analyses for the rest of the grids are similar and will not be detailed. distribution network can explore its potential of accommodating EVs and provide some guidance for the subsequent charging facility planning.

The Regional Distribution of OECL
EV ownership and EV types all correspond to the basic loads, and load types in the power grids, so it is necessary to explore the regional OECL distribution characteristics.
As shown in Figure 11, the load curve trough time of each subarea is identical and is around 08:00. In connection with a real-life situation, this period is generally the departure time of EVs. Besides, the OECL in Grid 1 is the largest among all the grids concerned because of the high EV ownership, and there are two load peaks around 13:00 and 22:00. This is because this subarea is mainly composed of residence and office areas, which reveals a behavior of getting out early and coming back at dusk. Therefore, there is a charging peak after EV owners go home at night; meanwhile, the charging of social vehicles and buses near the office area during the day forms another charging peak. The load structure of Grid 6 is consistent with that of Grid 1, so it presents a similar charging load profile. Grid 3 is the downtown area and charging loads mainly arrive during the day, while EV charging loads in Grid 4 mainly appear at night as this grid is dominated by residential areas. The charging load profile analyses for the rest of the grids are similar and will not be detailed.

OECL for Different Types of EVs
The OECLs for different types of EVs are obtained and compared with the ACE to explore the impact of EV categories on the OECL.
As shown in Figure 12, electric buses have a significant impact on the EV charging load profile, due to their high charging power and quick charge mode. Taxies mainly charge at night after the evening rush hours on the basis of a shift system. The charging load of a private car is mostly distributed in the load trough stage under the coordinated charging optimization, which has a prominent effect of filling the load valley. However, the charging load distribution of social vehicles is more random, but mainly concentrated between 14:00 and 23:00.
In general, the charging load of EVs is high between 13:00-02:00 (next day) and low during 04:00-11:00, which is basically consistent with the system load profile. Therefore, it is necessary to cope with EV charging behaviors in the future.

OECL for Different Types of EVs
The OECLs for different types of EVs are obtained and compared with the ACE to explore the impact of EV categories on the OECL.
As shown in Figure 12, electric buses have a significant impact on the EV charging load profile, due to their high charging power and quick charge mode. Taxies mainly charge at night after the evening rush hours on the basis of a shift system. The charging load of a private car is mostly distributed in the load trough stage under the coordinated charging optimization, which has a prominent effect of filling the load valley. However, the charging load distribution of social vehicles is more random, but mainly concentrated between 14:00 and 23:00.  Figure 13 shows the profile of the ACE and the OECL computed by the proposed model. It is obvious that the current distribution system has a wide capacity margin to accommodate EVs.

Accommodation Potential Analysis
With rapid developments of EVs and loads, the load accommodation capability of the distribution system will change accordingly, which will affect the ACE and its spatial and temporal distributions. Therefore, in order to test and verify the reliability of the current system, it is necessary to explore changes in the ACE and OECL, with respect to EV growth and load development. In this case, the current network structure and integration of flexible resources are kept unchanged, while the load level and EV ownership are projected to grow in the next five years. The results are shown in Figure 14.  In general, the charging load of EVs is high between 13:00-02:00 (next day) and low during 04:00-11:00, which is basically consistent with the system load profile. Therefore, it is necessary to cope with EV charging behaviors in the future. Figure 13 shows the profile of the ACE and the OECL computed by the proposed model. It is obvious that the current distribution system has a wide capacity margin to accommodate EVs.  Figure 13 shows the profile of the ACE and the OECL computed by the proposed model. It is obvious that the current distribution system has a wide capacity margin to accommodate EVs.

Accommodation Potential Analysis
With rapid developments of EVs and loads, the load accommodation capability of the distribution system will change accordingly, which will affect the ACE and its spatial and temporal distributions. Therefore, in order to test and verify the reliability of the current system, it is necessary to explore changes in the ACE and OECL, with respect to EV growth and load development. In this case, the current network structure and integration of flexible resources are kept unchanged, while the load level and EV ownership are projected to grow in the next five years. The results are shown in Figure 14.

Accommodation Potential Analysis
With rapid developments of EVs and loads, the load accommodation capability of the distribution system will change accordingly, which will affect the ACE and its spatial and temporal distributions. Therefore, in order to test and verify the reliability of the current system, it is necessary to explore changes in the ACE and OECL, with respect to EV growth and load development. In this case, the current network structure and integration of flexible resources are kept unchanged, while the load level and EV ownership are projected to grow in the next five years. The results are shown in Figure 14.
Consider the average daily profile of the ACE and the OECL in Figure 14. Most of the time within a day, the distribution system maintains a certain ACE margin over the OECL. However, there are times that both the ACE and OECL are relatively close, especially around 14:00. This narrower ACE margin is a result of the projected future EV and load growth. On the supply side, this indicates that the system needs to be expanded to cope with the load growth over time. It is also important to provide guidance to manage different charging behaviors, especially the usage of quick charging that has a higher impact on the distribution system. distribution system will change accordingly, which will affect the ACE and its spatial and temporal distributions. Therefore, in order to test and verify the reliability of the current system, it is necessary to explore changes in the ACE and OECL, with respect to EV growth and load development. In this case, the current network structure and integration of flexible resources are kept unchanged, while the load level and EV ownership are projected to grow in the next five years. The results are shown in Figure 14.  Figure 14. The average profile of the ACE and the OECL with projected EV and load growth. Figure 14. The average profile of the ACE and the OECL with projected EV and load growth.

Conclusions
This paper explores the ACE and the OECL of a distribution system with high penetrations of flexible resources, such as DGs and ESSs. At the same time, the spatiotemporal distributions of the ACE are investigated. The EV accommodation potential of a distribution system is also explored from the perspective of load growth. The following conclusions are drawn: