# Research on ISFLA-Based Optimal Control Strategy for the Coordinated Charging of EV Battery Swap Station

Authors: Xueliang Huang; Hao Qiang; Qidong Zhang; Haijuan Li
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

As an important component of the smart grid, electric vehicles (EVs) could be a good measure against energy shortages and environmental pollution. A main way of energy supply to EVs is to swap battery from the swap station. Based on the characteristics of EV battery swap station, the coordinated charging optimal control strategy is investigated to smooth the load fluctuation. Shuffled frog leaping algorithm (SFLA) is an optimization method inspired by the memetic evolution of a group of frogs when seeking food. An improved shuffled frog leaping algorithm (ISFLA) with the reflecting method to deal with the boundary constraint is proposed to obtain the solution of the optimal control strategy for coordinated charging. Based on the daily load of a certain area, the numerical simulations including the comparison of PSO and ISFLA are carried out and the results show that the presented ISFLA can effectively lower the peak-valley difference and smooth the load profile with the faster convergence rate and higher convergence precision.


---

## Body

## 1. Introduction

Over the past decades, many issues, such as energy shortages, serious environmental pollution, and global warming, have increasingly become worldwide concerns. EVs have emerged as a new traffic tool. Compared with internal combustion engine vehicles (ICEVs), which burn fossil fuels, EVs are driven by electricity. They demonstrate considerable advantages in solving the energy crisis and reducing the emissions of carbon dioxide, as well as in providing a means to drastically reduce the man-made pollution. More and more governments, car manufacturers, and energy companies are getting active in development and production of EVs [1].With the large-scale introduction of EVs, the power grid will face a significant challenge. Many domestic and foreign scholars have carried out researches on the impact of EVs on power distribution system [2–7], which mainly focus on the coordinated charging of EVs but rarely involve the optimal control strategy for the coordinated charging of EV battery swap station. Here, the coordinated charging means that the batteries are controlled to orderly charge for achieving an optimal objective such as minimizing the power losses or maximizing the grid load factor. Tian et al. established a mathematical model of dispatching strategy based on different objective functions which were solved by particle swarm optimization (PSO). The results show that the coordinated charging can lower the peak-valley difference and smooth the load profile [8].The SFLA is a global optimization algorithm proposed by Eusuff et al. [9, 10]. It is a memetic metaheuristic that is based on the evolution of memes carried by interactive individuals and a global exchange of information among the frog population. The SFLA draws its formulation from two other search techniques: the local search of the “particle swarm optimization” technique and the competitiveness mixing of information of the “shuffled complex evolution” technique. It locates a global optimum by combining global information exchange and local search, simulating the process of a group of frogs’ population-based cooperative seeking food. The strategies of local search and global information exchange make frog leap out of local optimum towards the global optimum solution (food place) with high probability. The algorithm has been developed to arrive at near-optimum solutions to complex and large-scale optimization problems which cannot be solved by gradient based mathematical programming techniques [11]. The SFLA has also been applied successfully to solving engineering problems such as traveling salesman problem (TSP) [12, 13], unit commitment (UC) problem [14–16], embedding virtual networks [17], speech emotion recognition [18], cognitive radio system [19], and other types of issues.In this paper, based on the characteristics of the EV battery swap station, the ISFLA-based optimal control strategy for coordinated charging has been investigated. Through the discretization of the solution vector, an ISFLA based on the reflecting method to deal with the boundary constraint is proposed. The comparison of PSO and ISFLA shows that the presented ISFLA can lower the peak-valley difference and smooth the load profile with the faster convergence rate and higher convergence precision.


## 2. Shuffled Frog Leaping Algorithm

As a bioinspired optimization technique, the SFLA is a metaheuristic optimization method, which imitates and models the behavior of frogs searching for food laid on discrete stones randomly located in a pond. In SFLA there is a population, which consists of a set of frogs (solutions). The set of frogs is partitioned into subsets referred to as memeplexes. Different memeplexes are considered as different cultures of frogs, and each memeplex performs a local search. Within each memeplex, the individual frogs hold ideas, which can be influenced by those of other frogs, and evolve through a process of memetic evolution. Frog leaping improves an individual’s meme and enhances its performance towards the goal. After a predefined number of memetic evolution steps, ideas are passed among memeplexes in a shuffling process. The local search and the shuffling processes continue until the defined convergence criteria are satisfied.The SFLA involves a population of possible solutions defined by a set ofprandomly generated frogs denoted asQ=[X1,X2,…,Xp]T, whereXi=(x1,x2,…,xiS)is the position of theith frog in theSdimension search space. All the frogs are sorted in a descending order according to their fitness, and the population is divided intommemeplexes, each containingnfrogs (i.e.,p=m×n), in such a way that the first frog goes to the first memeplex, the second frog goes to the second memeplex, themth frog goes to themth memeplex, and the(m+1)th frog goes back to the first memeplex, and the process continues in this manner.In each memeplex the frogs with the best and the worst fitness are represented byXbandXw, respectively. The best frog in the whole population is denoted byXg. During memeplex evolution, the worst frogXwleaps toward the best frogXb, which is formulated as the following updating rule:
(1)D=rand(0,1)×(Xb-Xw),Xw_new=Xw+D,-Dmax≤D≤Dmax,
whereDdenotes the updated step size for frog leaping, rand(0,1) generates a random number between[0,1], andDmaxis the maximum distance in one leaping.IfXw_newhas better fitness, the worst frogXwwill be replaced. Otherwise, the calculations in (1) are repeated with replacement ofXbbyXg. If no improvement occurs in this case, a new solution is randomly generated within the feasible space to replace the worst frogXw. Then, the calculations continue for a specific number of iterations. After a prespecified number of memetic evolutionary steps within each memeplex, to ensure global exploration, ideas passed within memeplexes are combined in the shuffling process. All the frogs are resorted, and the population is redivided intommemeplexes. The concurrently implemented local search and global shuffling continue alternatively until predefined convergence criteria are satisfied.


## 3. Charging Optimization Model for EV Battery Swap Station

Large-scale charging behavior of the EV will have a serious impact on the grid. Coordinated charging of the EV battery swap station can reduce the difference between the growing load peak and off-peak and save the costs of grid operation.


### 3.1. Charging Power Model of the EV Battery Charger

In the EV battery swap station, the typical strategy for battery charging is a two-stage method. The first stage has a constant current and limited pressure and the second stage has a constant pressure and limited current. Charging load power can be expressed as(2)P(t)={ImaxU(t),0≤t≤TC,UmaxI(t),TC<t≤TF,
whereImaxandUmaxdenote the maximum charging current and voltage, respectively.TCis the duration of charging with a constant current andTFis the total charging duration. These parameters are constants determined by the battery type and charging characteristics.In the first stage, to facilitate the modeling and analysis, the charging voltage is treated as the linear representation ofUmax. And in the second stage, the charging current is exponentially declined with time [20]. Then (2) can be rewritten as
(3)P(t)={ImaxUmax[(1-k)tTC+k],0≤t≤TC,UmaxImaxe-α(t-TC),TC<t≤TF,
whereαandkare also constants determined by the battery type and charging characteristics.In order to investigate the optimal control strategy for coordinated charging, the total charging durationTFis divided equally intoMsections, and the durationTMof each section can be derived asT/M. According to the equal area rule, the charging power of each section can be expressed with the average power:
(4)P-j=∫(j-1)TMjTMP(t)dtTM,1≤j≤M.Then an integerCis defined to satisfy the following expression:
(5)CTM≤TC≤(C+1)TM.Substituting (3) and (5) into (4) yields the multisection charging power model of battery charger:
(6)P-j={UmaxImax[(1-k)(2j-1)TM2TC+k]0<j≤CUmaxImax[(1-k)(TC2-C2TM2)2TC+k(TC-CTM)-eα(TC-CTM-TM)-1α(1-k)(TC2-C2TM2)2TC]C<j≤C+1UmaxImaxeα[TC-(j-1)TM]-eα(TC-jTM)αC+1<j≤M.Obviously, the accuracy of (6) is relative toM. A largerMimplies the higher precision of the model with the greater amount of computation, and vice versa.


### 3.2. Optimization Objective Function of the EV Battery Charging Power

According to the “technical guide for electric vehicle battery-swap station” presented by the State Grid Corporation of China (SGCC), a single EV battery charger is used to charge a single battery box. Inside the battery box, the battery pack consists of a plurality of battery cells. The battery mentioned in this paper refers to the battery pack. There are two types of batteries in the EV battery swap station, namely, the regulated battery connected to the EV battery charger and participating in the grid charging optimization by controlling its charging start time and the full charged reserve battery disconnected to the charger and used when not meeting the demand for swapping battery. Only the regulated-battery is considered in this paper.In order to investigate the optimal control strategy for coordinated charging of the EV battery swap station, some assumptions are made as follows.(
1
)All batteries are of the same type. The initial-charging SOC of each battery is the same and set to 20%, and the end-charging SOC is also the same and set to 100%.(
2
)Corresponding to the above multisection charging power model of the battery charger, one optimization periodTOis also divided into multisections(MO)with each section having the same durationTM. AndMO=TO·M/T. Then in one optimal periodTO, the numberNbof batteries needed for charging can be expressed as
(7)Nb=∑j=1MObj1≤j≤MO,
wherebjis the predicted number of batteries swapped into the station during thejth section. For optimal control, it is necessary to number these batteries needed for optimal charging. From theb0batteries, uncharged and reserved in the swap station at the beginning of theTO, to the batteries swapped into the station during theTO, the batteries are numbered consecutively untilNb. The other batteries swapped into the station of theTOare not indexed again and charged during the next optimization period.(
3
) For the batteryi,tin(i)andtout(i)denote its swap-in time and swap-out time, respectively:tin(i)=0
: batteryiis one of the batteriesb0;tin(i)=tj: batteryiis swapped in during thejth optimization section;tout(i)=tk: batteryiis swapped out during thekth optimization section;tout(i)=MO+1: batteryiis remained in the station.The battery can be swapped into and swapped out of the station at any time in one section. To facilitate the optimal control, we suppose that the battery could be charged from the section next to that of swapping in and be swapped out of the station from the section next to that of full charged. Then the start-charging timeTS(i)of the batteryiis betweentin(i)+1andtout(i)-M:
(8)tin(i)+1≤TS(i)≤tout(i)-M.To ensure that all batteries are full charged at the end of theTO, the swapped-in batteries during the lastMsections of theTOare no longer charged. Then the start-charging timeTS(Nb)of the batteryNbwould meet
(9)b0+∑j=1TS(Nb)-1bj≤Nb≤b0+∑j=1TS(Nb)bj,TS(Nb)≤MO-M.If there areif0full charged batteries at the beginning of theTOand all swapped out in thetfsection, the new full charged batteries would be swapped out from thetfsection. Then
(10)∑j=1tf-1bj≤if0≤∑j=1tfbj,tf≥M+1.(
4
)The battery is the core of EV and quite expensive. Many times of start and stop charging would greatly reduce its life, in which case the battery charging continuity constraint is introduced.Based on assumptions(1)and(4), the charging duration of all batteries is the same. In this paper, the minimum variance of the load profile including the charging loads generated by the EV battery swap station is the optimization objective and is expressed as
(11)minF=min[∑j=1MO(Pd,j+∑i=1NCi,jPi,j-Pavr)2]MOS.T.{b0+∑j=1TS(Nb)-1bj≤Nb≤b0+∑j=1TS(Nb)bjTS(Nb)≤MO-M∑j=1tf-1bj≤if0≤∑j=1tfbjtf≥M+1tin(i)+1≤TS(i)≤tout(i)-M,
where
(12)Pavr=∑j=1MO(Pd,j+∑i=1NCi,jPi,j)MO,Pd,j: the predicted daily load power of thejth section excluding the charging load generated by EV swap station;N: the total number of EV battery chargers;Ci,j: the charging decision variable of theith battery during thejth section: “1” means charging and “0” is discharging. Based on assumption (4),
(13)Ci,j={1TS(i)≤j≤TS(i)+M-1,0else;Pi,j: the charging load of theith battery during thejth section.


## 3.1. Charging Power Model of the EV Battery Charger

In the EV battery swap station, the typical strategy for battery charging is a two-stage method. The first stage has a constant current and limited pressure and the second stage has a constant pressure and limited current. Charging load power can be expressed as(2)P(t)={ImaxU(t),0≤t≤TC,UmaxI(t),TC<t≤TF,
whereImaxandUmaxdenote the maximum charging current and voltage, respectively.TCis the duration of charging with a constant current andTFis the total charging duration. These parameters are constants determined by the battery type and charging characteristics.In the first stage, to facilitate the modeling and analysis, the charging voltage is treated as the linear representation ofUmax. And in the second stage, the charging current is exponentially declined with time [20]. Then (2) can be rewritten as
(3)P(t)={ImaxUmax[(1-k)tTC+k],0≤t≤TC,UmaxImaxe-α(t-TC),TC<t≤TF,
whereαandkare also constants determined by the battery type and charging characteristics.In order to investigate the optimal control strategy for coordinated charging, the total charging durationTFis divided equally intoMsections, and the durationTMof each section can be derived asT/M. According to the equal area rule, the charging power of each section can be expressed with the average power:
(4)P-j=∫(j-1)TMjTMP(t)dtTM,1≤j≤M.Then an integerCis defined to satisfy the following expression:
(5)CTM≤TC≤(C+1)TM.Substituting (3) and (5) into (4) yields the multisection charging power model of battery charger:
(6)P-j={UmaxImax[(1-k)(2j-1)TM2TC+k]0<j≤CUmaxImax[(1-k)(TC2-C2TM2)2TC+k(TC-CTM)-eα(TC-CTM-TM)-1α(1-k)(TC2-C2TM2)2TC]C<j≤C+1UmaxImaxeα[TC-(j-1)TM]-eα(TC-jTM)αC+1<j≤M.Obviously, the accuracy of (6) is relative toM. A largerMimplies the higher precision of the model with the greater amount of computation, and vice versa.


## 3.2. Optimization Objective Function of the EV Battery Charging Power

According to the “technical guide for electric vehicle battery-swap station” presented by the State Grid Corporation of China (SGCC), a single EV battery charger is used to charge a single battery box. Inside the battery box, the battery pack consists of a plurality of battery cells. The battery mentioned in this paper refers to the battery pack. There are two types of batteries in the EV battery swap station, namely, the regulated battery connected to the EV battery charger and participating in the grid charging optimization by controlling its charging start time and the full charged reserve battery disconnected to the charger and used when not meeting the demand for swapping battery. Only the regulated-battery is considered in this paper.In order to investigate the optimal control strategy for coordinated charging of the EV battery swap station, some assumptions are made as follows.(
1
)All batteries are of the same type. The initial-charging SOC of each battery is the same and set to 20%, and the end-charging SOC is also the same and set to 100%.(
2
)Corresponding to the above multisection charging power model of the battery charger, one optimization periodTOis also divided into multisections(MO)with each section having the same durationTM. AndMO=TO·M/T. Then in one optimal periodTO, the numberNbof batteries needed for charging can be expressed as
(7)Nb=∑j=1MObj1≤j≤MO,
wherebjis the predicted number of batteries swapped into the station during thejth section. For optimal control, it is necessary to number these batteries needed for optimal charging. From theb0batteries, uncharged and reserved in the swap station at the beginning of theTO, to the batteries swapped into the station during theTO, the batteries are numbered consecutively untilNb. The other batteries swapped into the station of theTOare not indexed again and charged during the next optimization period.(
3
) For the batteryi,tin(i)andtout(i)denote its swap-in time and swap-out time, respectively:tin(i)=0
: batteryiis one of the batteriesb0;tin(i)=tj: batteryiis swapped in during thejth optimization section;tout(i)=tk: batteryiis swapped out during thekth optimization section;tout(i)=MO+1: batteryiis remained in the station.The battery can be swapped into and swapped out of the station at any time in one section. To facilitate the optimal control, we suppose that the battery could be charged from the section next to that of swapping in and be swapped out of the station from the section next to that of full charged. Then the start-charging timeTS(i)of the batteryiis betweentin(i)+1andtout(i)-M:
(8)tin(i)+1≤TS(i)≤tout(i)-M.To ensure that all batteries are full charged at the end of theTO, the swapped-in batteries during the lastMsections of theTOare no longer charged. Then the start-charging timeTS(Nb)of the batteryNbwould meet
(9)b0+∑j=1TS(Nb)-1bj≤Nb≤b0+∑j=1TS(Nb)bj,TS(Nb)≤MO-M.If there areif0full charged batteries at the beginning of theTOand all swapped out in thetfsection, the new full charged batteries would be swapped out from thetfsection. Then
(10)∑j=1tf-1bj≤if0≤∑j=1tfbj,tf≥M+1.(
4
)The battery is the core of EV and quite expensive. Many times of start and stop charging would greatly reduce its life, in which case the battery charging continuity constraint is introduced.Based on assumptions(1)and(4), the charging duration of all batteries is the same. In this paper, the minimum variance of the load profile including the charging loads generated by the EV battery swap station is the optimization objective and is expressed as
(11)minF=min[∑j=1MO(Pd,j+∑i=1NCi,jPi,j-Pavr)2]MOS.T.{b0+∑j=1TS(Nb)-1bj≤Nb≤b0+∑j=1TS(Nb)bjTS(Nb)≤MO-M∑j=1tf-1bj≤if0≤∑j=1tfbjtf≥M+1tin(i)+1≤TS(i)≤tout(i)-M,
where
(12)Pavr=∑j=1MO(Pd,j+∑i=1NCi,jPi,j)MO,Pd,j: the predicted daily load power of thejth section excluding the charging load generated by EV swap station;N: the total number of EV battery chargers;Ci,j: the charging decision variable of theith battery during thejth section: “1” means charging and “0” is discharging. Based on assumption (4),
(13)Ci,j={1TS(i)≤j≤TS(i)+M-1,0else;Pi,j: the charging load of theith battery during thejth section.


## 4. ISFLA-Based Optimal Control Strategy for Coordinated Charging

There areNchargers controlled by the control center of the EV battery swap station. From (11), the optimal control of the battery coordinated charging is a multivariable, nonlinear, and discrete integer optimization problem. Using the reflecting method to deal with the boundary constraint, an improved shuffled frog leaping algorithm (ISFLA) is proposed to achieve the optimal control strategy for coordinated charging of EV battery swap station.


### 4.1. Definition of Frog Position

In this paper, the position of the virtual frog consists of the start-charging time of every battery needed for coordinated charging, so the dimension of solution space is equal to the numberNbof the regulated battery. Theith frog can be expressed as
(14)Xi=(Tis(1),Tis(2),…,Tis(Nb)).Based on the assumptions,(15)Tis(h)∈{[1,MO-M]1≤h≤b0,[K+2,MO-M]∑j=0Kbj+1≤h≤∑j=0K+1bj,
where0≤K≤MO-M-2.


### 4.2. Inner-Memeplex Evolution

T
isis an integer, and the position of the frog with the worst fitness is adjusted as the following modified rule:
(16)dj=rand(0,1)×(Tbs(j)-Tws(j)),0≤j≤Nb,Tw_new(j)=Tws(j)+[dj],-dj,max≤dj≤dj,max,
where [dj] means an integer rounded todj.


### 4.3. Boundary Constraint

Dealing with the boundary constraint may affect the performance of the algorithm. In general, there are three basic approaches, namely, absorbing, reflecting, and damping. Using the reflecting method, when the updating position of the frog aftertiterations is beyond a certain boundary, the frog would return within the boundary with the original speed but in the opposite direction. In this paper, the reflecting method is chosen to deal with the boundary constraint and the position evolution of the frog obeys the following rule.Whendj>0,
(17)Twsnew(j)={Tws(j)-[dj]Tws(j)=MO-M,min(Tws(j)+[dj],MO-M)Tws(j)<MO-M;whendj<0,
(18)Twsnew(j)={Tws(j)-[dj]Tws(j)=tin(j)+1,max(Tws(j)+[dj],tin(j)+1)Tws(j)>tin(j)+1.


### 4.4. Steps of the ISFLA

The overall process of the ISFLA can be described in the following main steps.Step 0 (initialize the parameters).
Set the parameters ofp,n,ite_m, andite_p.pis number of frogs in the whole population.nis the number of frogs in each memeplex.ite_mis the maximum number of iterations for each memeplex.ite_pis the maximum number of shuffling iterations for whole population. Therefore,m=p/nis the number of memeplexes.Step 1 (initialize the population).
The initial population is formed byprandomly generated frogs:X1,X2,…,Xp. Compute the fitness valueF(Xi)according to (11) for each frogXi.Step 2 (shuffle the frogs).
Firstly sort thepfrogs in order of decreasing fitness, and then partition all thepfrogs intommemeplexes. The first frog goes to the first memplex, the second frog goes to the second memeplex, themth frog goes to themth memeplex, the (m+1)th frog goes back to the first memeplex, and at last thepth frog goes to themth memeplex.Step 3 (memetic evolution).
Repeat the following operations forite_mtimes: firstly find the worst frogXwand the best frogXbin each memeplex, then updateXwaccording to (16)–(18), and yieldXw_new. IfXw_newhas a better fitness, replaceXwbyXw_new. Otherwise, repeat the updated strategy withXgreplacingXb. If this still cannot produce a better solution, replace theXwby a randomly generated frog.Step 4 (check convergence).
If the convergence criteria are met, stop and output the best frogXgin the whole population. Otherwise, return to Step 2.In this ISFLA, the convergence criteria are defined as the maximum number of shuffling iterations for whole population (ite_p).


## 4.1. Definition of Frog Position

In this paper, the position of the virtual frog consists of the start-charging time of every battery needed for coordinated charging, so the dimension of solution space is equal to the numberNbof the regulated battery. Theith frog can be expressed as
(14)Xi=(Tis(1),Tis(2),…,Tis(Nb)).Based on the assumptions,(15)Tis(h)∈{[1,MO-M]1≤h≤b0,[K+2,MO-M]∑j=0Kbj+1≤h≤∑j=0K+1bj,
where0≤K≤MO-M-2.


## 4.2. Inner-Memeplex Evolution

T
isis an integer, and the position of the frog with the worst fitness is adjusted as the following modified rule:
(16)dj=rand(0,1)×(Tbs(j)-Tws(j)),0≤j≤Nb,Tw_new(j)=Tws(j)+[dj],-dj,max≤dj≤dj,max,
where [dj] means an integer rounded todj.


## 4.3. Boundary Constraint

Dealing with the boundary constraint may affect the performance of the algorithm. In general, there are three basic approaches, namely, absorbing, reflecting, and damping. Using the reflecting method, when the updating position of the frog aftertiterations is beyond a certain boundary, the frog would return within the boundary with the original speed but in the opposite direction. In this paper, the reflecting method is chosen to deal with the boundary constraint and the position evolution of the frog obeys the following rule.Whendj>0,
(17)Twsnew(j)={Tws(j)-[dj]Tws(j)=MO-M,min(Tws(j)+[dj],MO-M)Tws(j)<MO-M;whendj<0,
(18)Twsnew(j)={Tws(j)-[dj]Tws(j)=tin(j)+1,max(Tws(j)+[dj],tin(j)+1)Tws(j)>tin(j)+1.


## 4.4. Steps of the ISFLA

The overall process of the ISFLA can be described in the following main steps.Step 0 (initialize the parameters).
Set the parameters ofp,n,ite_m, andite_p.pis number of frogs in the whole population.nis the number of frogs in each memeplex.ite_mis the maximum number of iterations for each memeplex.ite_pis the maximum number of shuffling iterations for whole population. Therefore,m=p/nis the number of memeplexes.Step 1 (initialize the population).
The initial population is formed byprandomly generated frogs:X1,X2,…,Xp. Compute the fitness valueF(Xi)according to (11) for each frogXi.Step 2 (shuffle the frogs).
Firstly sort thepfrogs in order of decreasing fitness, and then partition all thepfrogs intommemeplexes. The first frog goes to the first memplex, the second frog goes to the second memeplex, themth frog goes to themth memeplex, the (m+1)th frog goes back to the first memeplex, and at last thepth frog goes to themth memeplex.Step 3 (memetic evolution).
Repeat the following operations forite_mtimes: firstly find the worst frogXwand the best frogXbin each memeplex, then updateXwaccording to (16)–(18), and yieldXw_new. IfXw_newhas a better fitness, replaceXwbyXw_new. Otherwise, repeat the updated strategy withXgreplacingXb. If this still cannot produce a better solution, replace theXwby a randomly generated frog.Step 4 (check convergence).
If the convergence criteria are met, stop and output the best frogXgin the whole population. Otherwise, return to Step 2.In this ISFLA, the convergence criteria are defined as the maximum number of shuffling iterations for whole population (ite_p).


## 5. Numerical Simulation

In this paper, the optimization objective is one area’s daily load profile which includes the charging loads generated by the EV battery swap station. The optimization periodTOis set to a day and divided into 24 sections with the durationTMof one hour.


### 5.1. Parameters of EV Battery

In the EV market, there are different battery types such as NiMH, Lead Acid, and Li-Ion. And the market share of Li-Ion with its technical advantages has increased annually. The “E6 pioneer” EV developed by BYD Co., Ltd., has been configured a Li-Ion battery with the energy density of 100 W·h per kilogram (i.e., a battery of 600 kg can store power of 60 kW·h for each charging). Its parameters are shown in Table1.Table 1
Basic parameters of battery.
Type
U
max  (V)
I
max  (A)
Constant current charging duration (min)
Constant voltage charging duration (min)
Li-Ion
58.8
120
90
150In consideration of the charging continuity constraint and based on the battery parameters, the average charging power of each section in the multisection charging model can be calculated and shown in Table2.Table 2
Average charging power of each section.
Section
1
2
3
4
Charging power (kW)
6.45
6.18
2.52
0.85


### 5.2. Area Daily Load

A certain area daily load power before regulation can be predicted and shown in Table3.Table 3
Certain area daily load before regulation.
Time (h)
Power (MW)
Time (h)
Power (MW)
Time (h)
Power (MW)
Time (h)
Power (MW)
1
5.25
7
9
13
8
19
17.75
2
4.25
8
10.5
14
7.5
20
16
3
4
9
9.5
15
7.75
21
14.25
4
4
10
9.25
16
8.25
22
9.75
5
4.25
11
9.75
17
13
23
7.25
6
5.75
12
12
18
18
24
5.75


### 5.3. Parameters of the EV Battery Swap Station

At the beginning of one optimization period, there are 600 batteries in the swap station, including 350 full charged batteries, 200 unfull changed batteries, and 50 reserved batteries. All batteries except the reserved batteries participate in the optimal control of coordinated charging. By predicting the demand of EV owners for swapping batteries during each optimization section, the numberbjcan be derived and shown in Table 4.Table 4
Predicted number of batteries swapped in during each optimization section.
Time (h)
1
2
3
4
5
6
7
8
9
10
11
12
Number
3
3
1
2
10
18
53
68
46
31
17
24
Time (h)
13
14
15
16
17
18
19
20
21
22
23
24
Number
30
30
37
50
67
42
28
28
25
18
13
6


### 5.4. Simulation of ISFLA

The main parameters of ISFLA are set and shown in Table5. Based on ISFLA and PSO, the area load profiles are compared and shown in Figure 1.Table 5
Parameters of ISFLA.
Parameter
p
m
n
ite_m
ite_p
Value
50
10
5
5
100Figure 1
Area load profile under different conditions.In Figure1, ISFLA and PSO can lower the load peak-valley difference and smooth the load profile. But the convergence characteristics of two algorithms in the numerical simulation differ greatly and are shown in Figure 2.Figure 2
Convergence characteristics of ISFA and PSO.Clearly, the convergence rate of ISFLA is much faster than that of PSO and the minimum variance of the load profile generated by ISFLA is less than that generated by PSO. Under different conditions, the variances of the load profile are shown in Table6.Table 6
Variances under different conditions.
Condition
Uncoordinated charging
Coordinated charging
ISFLA
PSO
Variance (kW2)
1.836
1.615
1.737
Result
—
Declining 12.04%
Declining 5.40%With ISFLA, the distribution of batteries according to the start-charging time in the optimization periodTOis shown in Figure 3.Figure 3
Distribution of batteries starting charging.In Figure3, during some optimization sections, there is a very small or no difference between the number of starting-charging batteries and that of those swapped into the station. This can explain that the optimized load profile with ISFLA and the load profile with uncoordinated charging almost coincide in some sections in Figure 1.For a fixed populationp, the different numbermof partitioned memeplexes will affect the convergence rate and the global optimal solution. Some simulations are made of ISFLA with the parameters of a fixedp(200) and differentmand the results are shown in Figure 4 and Table 7.Table 7
Global convergence value.
Memeplex
5
10
20
40
F (kW2) ×107
1.616
1.615
1.6139
1.6137Figure 4
Iterations of ISFLA with different memeplexes.In Figure4 and Table 7, the greater the number of memeplex is, the faster the algorithm converges with a higher convergence performance. However, the relationship between them is not linear, and how to choose an appropriate memeplex considering the convergence rate and global convergence value is a main goal of our further work.


## 5.1. Parameters of EV Battery

In the EV market, there are different battery types such as NiMH, Lead Acid, and Li-Ion. And the market share of Li-Ion with its technical advantages has increased annually. The “E6 pioneer” EV developed by BYD Co., Ltd., has been configured a Li-Ion battery with the energy density of 100 W·h per kilogram (i.e., a battery of 600 kg can store power of 60 kW·h for each charging). Its parameters are shown in Table1.Table 1
Basic parameters of battery.
Type
U
max  (V)
I
max  (A)
Constant current charging duration (min)
Constant voltage charging duration (min)
Li-Ion
58.8
120
90
150In consideration of the charging continuity constraint and based on the battery parameters, the average charging power of each section in the multisection charging model can be calculated and shown in Table2.Table 2
Average charging power of each section.
Section
1
2
3
4
Charging power (kW)
6.45
6.18
2.52
0.85


## 5.2. Area Daily Load

A certain area daily load power before regulation can be predicted and shown in Table3.Table 3
Certain area daily load before regulation.
Time (h)
Power (MW)
Time (h)
Power (MW)
Time (h)
Power (MW)
Time (h)
Power (MW)
1
5.25
7
9
13
8
19
17.75
2
4.25
8
10.5
14
7.5
20
16
3
4
9
9.5
15
7.75
21
14.25
4
4
10
9.25
16
8.25
22
9.75
5
4.25
11
9.75
17
13
23
7.25
6
5.75
12
12
18
18
24
5.75


## 5.3. Parameters of the EV Battery Swap Station

At the beginning of one optimization period, there are 600 batteries in the swap station, including 350 full charged batteries, 200 unfull changed batteries, and 50 reserved batteries. All batteries except the reserved batteries participate in the optimal control of coordinated charging. By predicting the demand of EV owners for swapping batteries during each optimization section, the numberbjcan be derived and shown in Table 4.Table 4
Predicted number of batteries swapped in during each optimization section.
Time (h)
1
2
3
4
5
6
7
8
9
10
11
12
Number
3
3
1
2
10
18
53
68
46
31
17
24
Time (h)
13
14
15
16
17
18
19
20
21
22
23
24
Number
30
30
37
50
67
42
28
28
25
18
13
6


## 5.4. Simulation of ISFLA

The main parameters of ISFLA are set and shown in Table5. Based on ISFLA and PSO, the area load profiles are compared and shown in Figure 1.Table 5
Parameters of ISFLA.
Parameter
p
m
n
ite_m
ite_p
Value
50
10
5
5
100Figure 1
Area load profile under different conditions.In Figure1, ISFLA and PSO can lower the load peak-valley difference and smooth the load profile. But the convergence characteristics of two algorithms in the numerical simulation differ greatly and are shown in Figure 2.Figure 2
Convergence characteristics of ISFA and PSO.Clearly, the convergence rate of ISFLA is much faster than that of PSO and the minimum variance of the load profile generated by ISFLA is less than that generated by PSO. Under different conditions, the variances of the load profile are shown in Table6.Table 6
Variances under different conditions.
Condition
Uncoordinated charging
Coordinated charging
ISFLA
PSO
Variance (kW2)
1.836
1.615
1.737
Result
—
Declining 12.04%
Declining 5.40%With ISFLA, the distribution of batteries according to the start-charging time in the optimization periodTOis shown in Figure 3.Figure 3
Distribution of batteries starting charging.In Figure3, during some optimization sections, there is a very small or no difference between the number of starting-charging batteries and that of those swapped into the station. This can explain that the optimized load profile with ISFLA and the load profile with uncoordinated charging almost coincide in some sections in Figure 1.For a fixed populationp, the different numbermof partitioned memeplexes will affect the convergence rate and the global optimal solution. Some simulations are made of ISFLA with the parameters of a fixedp(200) and differentmand the results are shown in Figure 4 and Table 7.Table 7
Global convergence value.
Memeplex
5
10
20
40
F (kW2) ×107
1.616
1.615
1.6139
1.6137Figure 4
Iterations of ISFLA with different memeplexes.In Figure4 and Table 7, the greater the number of memeplex is, the faster the algorithm converges with a higher convergence performance. However, the relationship between them is not linear, and how to choose an appropriate memeplex considering the convergence rate and global convergence value is a main goal of our further work.


## 6. Conclusion

As a main way of energy supply to EVs, the optimal control strategy for coordinated charging of the swap station is very important in smoothing the load profile. Based on the characteristics of the EV battery swap station, a multisection charging power model of battery is presented and an ISFLA in reflecting method to deal with the boundary constraint is proposed to achieve coordinated charging of batteries. In numerical simulations, the comparison of PSO and ISFLA is made, and the results show that the presented ISFLA can effectively lower the peak-valley difference and smooth the load profile with the faster convergence rate and higher convergence precision.



---

*Source: 613404-2013-11-26.xml*