Next Article in Journal
Using Hydrogen Reactors to Improve the Diesel Engine Performance
Next Article in Special Issue
Study of Grid-Connected PV System for a Low Voltage Distribution System: A Case Study of Cambodia
Previous Article in Journal
Conjugate Heat Transfer Analysis of the Aero-Thermal Impact of Different Feeding Geometries for Internal Cooling in Lifetime Extension Processes for Heavy-Duty Gas Turbines
Previous Article in Special Issue
Appropriate Volt–Var Curve Settings for PV Inverters Based on Distribution Network Characteristics Using Match Rate of Operating Point
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Nature-Inspired Algorithm to Enable the E-Mobility Participation in the Ancillary Service Market

by
Davide Falabretti
* and
Francesco Gulotta
Department of Energy, Politecnico di Milano, Via Lambruschini 4, 20156 Milan, Italy
*
Author to whom correspondence should be addressed.
Energies 2022, 15(9), 3023; https://doi.org/10.3390/en15093023
Submission received: 21 March 2022 / Revised: 13 April 2022 / Accepted: 15 April 2022 / Published: 20 April 2022

Abstract

:
In the present paper, a tool is proposed to optimally schedule the charging requests of a fleet of carsharing Electric Vehicles (EVs) in an urban area, to enable their participation in the Ancillary Service Market. The centralized scheduler minimizes the imbalance of an EV fleet with respect to the power commitment declared in the Day-Ahead Market, providing also tertiary reserve and power balance control to the grid. The regulation is carried out by optimizing the initial charging time of each vehicle, according to a deadline set by the carsharing operator. To this purpose, a nature-inspired optimization is adopted, implementing innovative hybridizations of the Artificial Bee Colony algorithm. The e-mobility usage is simulated through a topology-aware stochastic model based on carsharing usage in Milan (Italy) and the Ancillary Services requests are modeled by real data from the Italian electricity market. The numerical simulations performed confirmed the effectiveness of the approach in identifying a suitable schedule for the charging requests of a large EV fleet (up to 3200 units), with acceptable computational effort. The benefits on the economic sustainability of the E-carsharing fleet given by the participation in the electricity market are also confirmed by an extensive sensitivity analysis.

1. Introduction

Environmental and technical factors are driving the transition toward a low-carbon energy scenario. Electricity and road transport, underpinned by environmental policies, are among the sectors facing the greatest transformations. As a consequence of the spreading of renewable energy sources, electric power systems are subject to an increasing variability and unpredictability of power injections and to a gradual replacement of conventional generators, historically in charge to supply Ancillary Services (ASs) to the grid, with non-dispatchable generation. To improve the power system’s reliability, recently, some countries around the world envisaged the opening of their Ancillary Service Markets (ASMs) to new subjects, such as distributed energy resources, renewable power plants, controllable loads, energy storage systems and also e-mobility [1,2,3], with the goal of increasing the number of suppliers of the ASs needed for grid operation. In this scenario, electric networks can take advantage of synergies with the transport sector, given the deep changes the latter is facing to reduce local pollutants and increase the efficiency of urban space usage [4]. In this regard, the promotion of electric carsharing could strongly favor the shift toward a sustainable transport sector [5]. Moreover, Electric Vehicles (EVs) are a very promising flexibility resource, even though their charging requests need to be properly coordinated to avoid detrimental effects on the grid [6].
In the outlined framework, the present paper proposes a novel optimization procedure based on a hybridization of the Artificial Bee Colony algorithm aimed to schedule the charging requests of an EV carsharing fleet [7]. The initial charging time of each vehicle is coordinated by adopting a Centralized Control Architecture (CCA). To preserve batteries’ lifespan, only the starting time of the charging is managed, while no preemption nor Bidirectional Vehicle to Grid (V2G) are considered. The strategy enables the aggregated participation of e-mobility in the ASM, the supply of ASs to the grid (tertiary reserve and power balance control) and the reduction of power imbalances with respect to the market commitment defined in the Day-Ahead Market (DAM).
In the following, after a literature review in Section 2, the scenario of application and its modeling are depicted in Section 3. The description of the dynamic optimization problem studied is presented in Section 4. Section 5 introduces the enhancements proposed to adapt the Artificial Bee Colony to this large-scale optimization problem. Finally, Section 6 and Section 7 present the case study and the numerical analyses carried out to assess the effectiveness of the proposed solution, while in Section 8, the techno-economic profitability of the CCA is analyzed in different case studies as well as performing a sensitivity analysis.

2. Related Works

In recent years, the exploitation of the flexibility of EVs to supply services to the power grid has shown an increasing interest in the scientific community. Regulation capabilities, scheduler architecture and optimization strategies to approach the scheduling problem are all topics widely debated today.
The capability of EVs to manage unidirectional or bidirectional power flows deeply affects the performance during the service provision to the power system. In unidirectional grid-to-vehicle (hereinafter V1G), EVs can be considered as a flexible load, since the power can flow only from the grid to the vehicle [8], whilst the bidirectional grid to vehicle (V2G) enables EVs both to absorb or inject electric power from/into the grid. As concluded in [9], the bidirectional configuration has higher performances than V1G, but battery degradation, charging infrastructure and control logics are still crucial problems [10].
From the point of view of the scheduler architecture, in the centralized control, the operator managing the regulation of EVs and offering their flexibility services on the market (in the following, also called Aggregator) collects all the information about the connected EVs, determines the optimal charging scheduling and coordinates the fleet’s power exchanges with the grid to obtain the services settled on the market. The opposite of this approach is the decentralized control logic, where each EV is capable to determine its charging parameters. Most of the literature focuses on the centralized logic [11,12], since it allows usually reaching better coordination and provides a higher amount of services to the power system. However, the centralized approach could show severe scalability limits [13], since the EV scheduling is a dynamic optimization problem [14], therefore, the solution shall be available in a few seconds to start the EV charging and reiterate the procedure for the next time steps. It is not surprising that to face this issue, in most of the works adopting a centralized architecture, the size of EVs’ fleet is very limited [15]. To address the scalability issues of centralized architectures, recently, Meta-Heuristic optimization techniques have met an increasing interest. Genetic Algorithms, Particle Swarm Optimization and Artificial Bee Colony (ABC) have been extensively applied to solve EV charging scheduling problems [16,17,18].
In this work, an innovative hybridization of the ABC algorithm is provided to dynamically schedule the charging requests of a large carsharing fleet, modeling stochastically the EVs’ behavior and the ASM. Actually, ABC is suitable for challenging optimization problems since it systematically incorporates exploration and exploitation mechanisms [19]. Nevertheless, in literature, it is a common practice to overcome the disadvantages of ABC by hybridization techniques [20]: most of the works proposing changes to the ABC algorithm focus on the procedure for setting the initial populations [16] or on the local research process [21]. Starting from these considerations, the proposed approach (h-ABC) aims to efficiently solve the scheduling problem, by improving both the population initialization and the local research process compared to the standard ABC.
In this regard, in studies [7,22], the authors presented a preliminary architecture to enable the e-mobility participation in the market for the AS provision, focusing also on the undesirable effects caused to the distribution grid. However, neither the details of the implemented optimization algorithm nor the approach adopted to simulate EV usage have been reported. These aspects are described in detail in this work, also proving the effectiveness of the proposed method by comprehensive techno-economic analyses.

3. System Model

In the present paper, Ancillary Services (ASs) are supplied to the power system by scheduling the charging requests of a fleet of EVs managed by a carsharing operator, who is assumed to act also as Balancing Service Provider (BSP) [23], aggregating the regulation capacity of the fleet and offering it on the electricity market. To this purpose, the CCA optimizes the charging scheduling and communicates with the Charging Stations (CSs) enabling the recharge of car batteries. The purpose of the envisaged architecture is to enable a carsharing operator to participate in the Ancillary Service Market (ASM), thus increasing his profits, but also taking care to avoid detrimental effects on the quality of carsharing service provided to customers.
In this scenario, schematically represented in Figure 1, at day D-1, the carsharing operator (operating as an Aggregator) submits requests on the DAM to buy the energy requested the next day (D) to charge the EVs. As a result, a binding day-ahead schedule (hereinafter Power Baseline Schedule) is defined. During the relevant ASM sessions (for instance, according to the scenario in place in Italy, 4 h before the real-time), the Aggregator submits bids for up/downward reserve, which are selected by the Transmission System Operator (TSO) through a pay-as-bid approach and, if accepted, are notified to the Aggregator shortly before the real-time (i.e., 15 min in advance). The combination of DAM and ASM schedules is defined as Power Request Schedule, which represents a commitment toward the market that the Aggregator must respect; otherwise, it is penalized by imbalance fees. In the real-time, EVs are used by carsharing customers. Each user at the end of the rent could connect the car to one of the CSs over the city, or park it in a public spot without a CS. When the EV is plugged into a CS, the CCA optimizes the initial charging time to meet the market commitments, while preserving the final battery’s State of Charge (SoC) according to a deadline defined by the carsharing operator.
To assess the performances of the proposed scheduler, in this work, suitable numerical tools have been developed to model the participation of the carsharing operator in the electricity markets (i.e., DAM and ASM) and to simulate cars’ usage by users. An in-depth description of these tools is provided in the following subsections.

3.1. Electricity Market Participation

In the DAM, the Aggregator purchases the Power Baseline Schedule, which defines for every hour of the day D the amount of power P 0 ( t ) that should be absorbed to charge the EVs’ fleet. An accurate forecast of energy requirements for the next day is essential since a mismatch between the power actually absorbed in real-time and the scheduled one implies an economic penalty (i.e., imbalance fee). The development of forecasting algorithms is out of the scope of this work; thus, a simplistic approach is adopted to estimate future power withdrawals from the grid, based on the persistence model [24]. In the session of the DAM on the day D-1, the Power Baseline Schedule is forecasted, assuming as a reference the profile of charging requests registered on the previous day (D-2). In particular, the Power Baseline Schedule ( P 0 ) is defined, assuming all the charging requests at day D-2 are satisfied as soon as the EVs were connected to the CS:
P 0 ( t ) = i = 1 N r e q D 2 ( t ) P c h g i , D 2 ( t ) t T D
where T D is the vector of all time steps in day D, while N r e q D 2 ( t ) and P c h g i , D 2 ( t ) are, respectively, the number of EVs charging at time t and the power absorbed by the i -th EV.
Then, during the ASM sessions, the Aggregator submits bids for upward and downward reserve. In case of acceptance, AS requests ( P A S ) are notified 15 min before the real-time as power variations with respect to the Power Baseline Schedule, as defined in Equation (2), where a load sign convention is adopted:
P r e q ( t ) = P 0 ( t ) + P A S ( t )
Hence, in the real-time, the CCA will adjust the power absorption of the EVs by scheduling their charging, in order to meet the corresponding Power Request Schedule (i.e., P r e q ).

3.2. Vehicle Model

The performance of the proposed CCA is evaluated by a carsharing usage model, realistically reproducing the common behavior of customers. A fleet composed by N f l e e t t o t EVs is simulated, adopting a stochastic approach aware of the city topology, as in [25]. The approach, schematically presented in Algorithm 1, simulates each vehicle of the fleet, denoted by id ( i d N f l e e t = [ 1 , N f l e e t t o t ] ), in an arbitrary number of days ( d N d a y = [ 1 , N d a y s t o t ] ) .
For each EV, firstly, the total number of daily rent requests ( n r e n t i d , d   ) is extracted from a normal integer distribution, with a mean value equal to 4. Then, the EV usage in each rent ( r [ 1 ,   n r e n t i d , d ] ) is defined by extracting the initial time ( t i n r , i d ), the duration of the rent ( Δ t r e n t r , i d ) and the distance covered ( Δ d r e n t r , i d ). All these data are obtained from statistics based on the actual carsharing usage in Italy, as described in Section 5. Incoherent rents are discarded (line 9 of Algorithm 1). At the end of the rent, the customer can choose to plug the EV into one of the CSs dispersed in the city or to park it in a public car park without a CS. To model this decision, a biased Boolean extraction is performed, assuming that one out of four customers will plug the EV into the CS, so that each EV is charged, on average, once per day. If the user opts to connect the EV to a CS, from the set of all CSs distributed in the city, denoted by C S , only the ones far from the initial position ( d c s - i n t ) less than the distance covered during the rent (i.e., d c s - i n t < Δ d r e n t r , i d ) are considered. Finally, one of these admissible CSs ( C S a d m i d ) is randomly selected:
C S a d m i d { C S   ( d   c s i n t < Δ d r e n t r , i d ) }
This way, the motion of each EV in the fleet is simulated, tracking its path in the city and updating the relevant SoC:
S o C C S , c o n r , i d = S o C i n t r , i d Δ d r e n t r , i d e ¯ c o n s i d / C b r t i d
where S o C i n t r , i d is the SoC of the EV before the rent, e ¯ c o n s i d and C b t r i d are, respectively, the average energy consumption per kilometer and the battery capacity for the id-th EV. Once the EV is plugged in, the CS collects the relevant SoC ( S o C C S , c o n r , i d ), charging power ( P c h g i d ) and charging deadline time ( D t i m e i d , r ). Coherently with the market dynamics rules, a quarter of hour resolution has been adopted ( τ = 0.25 h) for simulating the car behavior: therefore, every 15 min, data of new EVs connected to the CSs are collected and sent to the CCA, which sets the initial charging time by solving the optimization process.
Algorithm 1. Fleet Modeling
1: for d = 1:   N d a y s t o t  do
2:      for i d   = 1:   N f l e e t t o t do
3 :                       n i d , d r e n t ← number of rents
4:            for  r = 1 :   n r e n t i d , d do
5:                                     t i n r , i d ← starting rent time
6:                                                   Δ t r e n t r , i d ← rent duration
7:                                                   Δ d r e n t r , i d ← distance covered
8:                                                   v r , i d =   Δ d r e n t r , i d / Δ t r e n t r , i d ← mean speed during rent
9:                              if  v r , i d > v M A X or t i n r + 1 , i d > D t i m e i d , r
10:                              go to line 5
11:                             end if
12:                                                   r n d C S ← Boolean extraction (if TRUE, an EV is connected to a CS)
13:                             if r n d C S = TRUE
14 :                                                     select   the   arrival   CS   from   the   set   of   C S a d m i d
15:                              end if
16:                end for
17:        end for
18: end for

4. Optimization Problem

The CCA deals with a dynamic scheduling optimization problem, because the EVs’ charging requests and the relevant characteristics at the arrival are not known in advance. Therefore, at each time step ( t t t o t = [ 1 , N d a y s t o t · 24 τ ] ), the CCA shall perform an optimization based on its best knowledge of the problem. In this work, to improve the capability of the CCA to find the optimal scheduling, a forecast is implemented to predict the charging requests in the upcoming time steps ( T f ) and to plan consequently the EVs’ charging in the next future. The optimal solution found in the current time step ( t 0 ) depends also on the forecasted future events (i.e., charging requests). Hence, during each optimization process, performed at the time t 0 , a moving time window is considered, including the present time and the time steps in the next hour ( T f ):
t o p t = [ t 0 , T f ] = { t 0 , t 0 + τ , , t 0 + 4 τ }
The number of EVs expected connecting in upcoming time steps ( N e x p ( t f )   t f T f ) and the relevant characteristics (SoC, rated power, etc.) are predicted based on historical data, while the same information for EVs already connected at the time t 0 is known. Finally, in each time step t   ( t t t o t ) , the instance   P t   is defined, including the set of all EVs (predicted or already connected) whose charging time should be optimized.
The optimization problem aims at minimizing, in each time step ( t t o p t ), the costs for the Aggregator:
F o b j ( t ) = C o s t i m b ( t ) + C o s t d l n ( t )
The first term,   C o s t i m b ( t ) , models the imbalance fees applied to the Aggregator in the considered time window ( t o p t ) in case of mismatches between the actual power absorbed from the network and the power schedule resulting from the participation in DAM and ASM sessions:
C o s t i m b ( t ) = t o p t t o p t t o p t   c i m b a ( t o p t ) τ | P a b s ( t o p t ) P r e q ( t o p t ) |
Imbalance costs are assumed proportional to the amount of energy imbalance, therefore, c i m b is a unitary imbalance fee (EUR/MWh) and the function a ( t o p t ) is a hyper-parameter that weighs the foreseen imbalances with respect to the current one, having a ( t o ) = 1 and a ( t o p t ) < 1 for t o p t > t 0   t o p t t o p t . Finally, P r e q is the Power Request Schedule defined in Equation (2) and P a b s is the power absorbed by the CSs, calculated using:
P a b s ( t o p t ) = e v = 1 N c h g ( t o p t ) P c h g e v                                   t o p t t o p t
where P c h g e v is the power absorbed by the ev-th vehicle and N c h g ( t o p t ) is the total number of EVs scheduled to charge in time step t o p t . For the sake of simplicity, a constant charging power is assumed (hypothesis usually true for the slow charging strategy commonly adopted for low-cost carsharing EVs).
The second term of the objective function (Equation (6)) is addressed to preserve the quality of the carsharing service. It assumes that the carsharing operator’s loss of profit is proportional to the reduction of the EVs’ range autonomy caused by an incomplete battery charging:
C o s t d l n ( t ) = c d l n e v e v P t E r r d l n e v
with
{ E r r d l n e v = 0                                                   i f     s t e v + Δ t c h g e v D t i m e e v E r r d l n e v = s t e v t i n e v                 i f     s t e v + Δ t c h g e v > D t i m e e v
Δ t c h g e v = ( S o C F i n S o C C S , c o n e v ) C b t r e v P c h g e v
where c d l n is the cost associated with a unitary (15 min) postponement of the charging, which results in a reduction of the SoC at the CS disconnection. E r r d l n e v is a parameter different from zero if the proposed initial charging time ( s t e v ) does not allow reaching the maximum SoC expected before the disconnection in the time step D t i m e e v . The charging duration ( Δ t c h g e v ) is calculated with Equation (11), where C b t r e v is the battery capacity and S o C F i n and S o C C S , c o n e v are, respectively, the SoC to achieve during the charging (i.e., assumed equal to 1) and the battery SoC when the vehicle is plugged into the CS, obtained by means of the carsharing model described in Section 3.
It is worth noting that the adoption of the moving time windows technique in the optimization process allows obtaining a solution (optimal EVs’ charging schedule) which is a function of the future events (i.e., power variations due to AS requests in the imminent future).
The optimization problem, formulated as follows, should be solved in all time steps ( t t t o t ), updating at each step the EVs’ data and the Power Request Schedule:
M i n         F o b j ( t ) t t t o t s t e v               e v P t
subject to
t i n e v s t e v D t i m e e v   e v P t
P a b s ( t o p t ) = e v = 1 N c h g ( t o p t )   E r r d l n e v t o p t t o p t  
P r e q ( t 0 ) = P 0 ( t o ) + P A S ( t 0 )  
P r e q ( t f ) = P 0 ( t f ) t f T f  
0 < S o C F i n 1 e v P t  
The objective function in Equation (12) should be minimized by optimizing the initial charging time ( s t e v ) of each EV in the considered instance ( e v P t ) . Constraint in Equation (13) specifies the lower and upper bounds of the charging time for each EV, while constraint in Equation (14) defines the power absorbed by the Aggregator during the considered time window ( t o p t ). Equations (15) and (16) determine, respectively, the Power Request Schedule in the current time step ( t 0 ) and in the upcoming ones ( T f ). Finally, the constraint in Equation (17) models the physical limits of the EV’s battery.
To solve the described optimization problem, in this paper, a hybridization of the ABC algorithm (h-ABC), is proposed.

5. The h-ABC Approach

The Artificial Bee Colony is a population-based Meta-Heuristic algorithm inspired by the foraging behavior of bees. It is a swarm-based optimization method, so the   N F S different solutions are simultaneously evaluated and improved, thanks to the iteration of three phases: employed, onlooker and scout ones [26]. ABC is based on a random scattering of the initial solutions and then the reiteration of a searching process, which is designed without considering the nature of the problem (pseudo-random research). This can negatively affect the computational effort to find the optimal solution, a key parameter in large-scale dynamic problems. Therefore, in this work, the ABC is improved by deeply revising two phases: (1) the definition of the initial solutions; (2) the research process during the employed bee phase.

5.1. Initial Population

The standard ABC algorithm initializes each candidate solution ( p [ 1 ; N F S ]) by randomly arranging them in the allowed domain (i.e., meeting the constraints).
In the h-ABC approach, two phases are implemented to define the initial solutions: firstly, the charging order of EVs is defined in each initial solution by dispatching rules, and then, a tournament selection mechanism is applied. Indeed, Meta-Heuristic algorithms, such as ABC, have been found to be more effective if the initial population is selected by sub-optimal scheduling rules (i.e., first-in first-served approach, etc.) [27]. The dispatching rules adopted are:
  • Due to Time (DtT) in which the charge priority is computed considering the difference between D t i m e e v   and Δ t c h g e v .
  • Arrival Time (AT) in which the EV charge priority is determined considering a first-in first-served approach.
The charging order vector for each initial solution ( v p   p ) is defined by the following procedure: firstly, the EVs not yet added to the initial solution are sorted using a dispatching rule (i.e., DtT or AT), then a tournament selection is applied to select h individuals from the sorted vector, adding to the charging order v p the best of these h individuals. One-half of the initial solutions are initialized by each dispatching rule and, to further increase the population diversification, the tournament size is decreased after the definition of each charging order.

5.2. Employed Bee

After the initial population definition, both in ABC and h-ABC, an iterative process starts. The first step (employed bee) performs local research around each solution p ( p [ 1 ; N F S ]). In the ABC, the local research process initializes a new candidate solution ( e     p ) in the neighbor of the current one ( s t p ), adopting Equation (18) [28]:
e e v , p = s t e v , p + Φ ( s t e v , p s t e v , k )  
where ev and k are randomly chosen indexes with e v P t   and k [ 1 ; N F S ] and k p. Φ is a random number between −1 and +1, extracted from a uniform distribution. Then, between the new candidate solution and the old one, a greedy selection process is carried out and, in case of objective function’s improvement, the old solution is replaced with the new one.
The research process adopted in h-ABC relies on a more advanced heuristic research. The approach (Algorithm 2) can be subdivided into two phases. In the first phase, each term of Equation (7) is inspected and the one with the highest contribution is selected as the critical one ( t c r i t ). Then, in the second phase, the algorithm evaluates if, in t c r i t , a surplus or a lack of power absorbed with respect to the requested one occurs. If the power absorbed exceeds the requested one ( P a b s ( t c r i t ) > P r e q ( t c r i t ) ), the EV having the minimum tardiness, defined in Equation (19), among all the ones scheduled to charge in t c r i t , is selected. Then, its recharge is postponed in a time step t p o s t ,   with t p o s t > t c r i t , having a non-negative imbalance. In absence of a time step t p o s t satisfying both t p o s t > t c r i t and the imbalance sign constraint, the selected EV is simply postponed in t c r i t + 1. Indeed, a ( T f ) is a decreasing function, therefore, by shifting forward in the future the imbalance, the objective function in the current time step is reduced. This approach is justified considering that in upcoming time steps, the new connected EVs and possible changes in the Power Request Schedule could eliminate foreseen imbalances. In the opposite case, if P a b s ( t c r i t )   <   P r e q ( t c r i t ) , the tardiness of the EVs whose charging is scheduled in t o p t > t c r i t   t o p t t o p t is evaluated; the EV with the maximum tardiness is selected and its charging time is set to t c r i t . Obviously, during the EV selection, the constraints defined in Equation (17) shall be respected.
t a r d e v = M a x ( 0 ; s t e v + Δ t c h g e v D t i m e e v ) e v P t  
Algorithm 2. Heuristic research
1: for p = 1:   N F S    do
2:     t c r i t ← time step with the highest C o s t i m b p
3:    if  P a b s p ( t c r i t ) > P r e q   ( t c r i t )
4:         find e v   such that s t e v , p = t c r i t and select the one with M a x ( t a r d e v )
5 :             t c r i t     P a b s p ( t p o s t ) < P r e q ( t p o s t )   and   t p o s t   >   t c r i t
6:           s t e v , p = t p o s t
7:      else
8:             find e v   such that s t e v , p = t c r i t and select the one with M i n ( t a r d e v ), so     s t e v , p = t c r i t
9:    end if
10: end for

5.3. Onlooker Bee

The onlooker bee phase is designed to further explore promising regions of the domain by the sharing of information between different solutions. This phase is implemented in the same way in ABC and h-ABC. Firstly, a fitness function is calculated for each solution ( f i t p   p ) as the inverse of the corresponding objective function. Then, the probability distribution function shown in Equation (20) is used to perform a roulette wheel selection, which associates each onlooker bee ( o [ 1 ; N F S ] ) to different solutions ( p [ 1 ; N F S ]):
P r b p = f i t p p f i t p p [ 1 , N F S ]  
To initialize a new solution ( v c h g o ) around the corresponding solution p extracted in the roulette wheel process, Equation (18) is once again applied. Finally, the new solution and the old one are compared through a greedy selection process, keeping the best one. To better control the exploitation of the domain, the parameter L i m p is introduced: it is initialized to zero for all the solutions and each time an onlooker bee does not improve the current solution ( p ), L i m p is increased by one.

5.4. Scout Bee

The last phase in both ABC and h-ABC is the scout bee: it is activated when the parameter L i m p reaches a fixed upper bound ( L i m M A X ). If this occurs, the algorithm considers the area of the domain around the solution p fully exploited, discards it and the solution is replaced with a new one, randomly generated.

5.5. Convergence Criteria

In the ABC, the three phases are iteratively repeated until a maximum number of iterations ( N m a i n ); then, the solution with the lowest value of the objective function is selected as the optimal one. In h-ABC, a more advanced stopping criterion is introduced in addition to the one just described: after a fixed number of iterations ( N c y l , m i n ), the objective function trend is analyzed. If the populations do not improve their objective function for a given number of iterations (empirically set to 5 in the numerical analyses that follow), the iterative process stops, avoiding the unnecessary computation of cycles not improving the solutions.
When the convergence criterion is reached, the EVs’ charging schedule with the lowest objective function value is adopted: EVs with s t e v = t 0 will immediately start the charging, whilst the charging of the remaining ones ( s t e v > t 0 ) will be rescheduled in the optimization performed in the next time steps by considering possible updates to the connected EVs and to the Power Request Schedule.

6. Case Study

In this work, the participation of a carsharing operator, acting as an Aggregator by selling the EVs’ flexibility on the Ancillary Service Market, has been reproduced through the historical data collected on the Italian ASM in 2018. To this purpose, all the bids submitted to the ASM by three conventional Power Units for balancing and tertiary reserve have been considered in terms of offered and awarded power and prices. Data have been processed to remove periods without offers (maintenance stops, etc.) and scaled to match the EVs aggregate’s rated power. Finally, to generalize the market requests delivered to the Aggregator, the pool of AS requests has been extended by sampling the historical market data in groups of 7 consecutive days and randomly combining them together.
For the evaluation of the economic flows deriving from the market participation, the following scheme is adopted. In the case of imbalances between the Power Request Schedule and actual absorbed power ( P a b s ), an imbalance penalty is applied to the Aggregator. If P a b s > P r e q , the power surplus (not declared in the DAM) is paid by the Aggregator p D A M + c i m b , where p D A M   is the price cleared in the DAM and c i m b is the imbalance fee. Vice versa, if P a b s < P r e q , the power not absorbed by EVs is refunded to the Aggregator at price p D A M c i m b . Therefore, in each time step in which the power schedule ( P r e q ) is not fully respected, the Aggregator is subject to an imbalance cost ( C i m b ) calculated as C i m b = c i m b | P a b s P r e q | . For the sake of simplicity, p D A M and c i m b are taken constant and respectively equal to 50 EUR/MWh and 80% of p D A M ( c i m b = 40 EUR/MWh).
Concerning the AS provision, in the case of upward regulation (i.e., reduction of the power absorbed), the Aggregator is paid p D A M + p A S , to refund the energy purchased in the DAM not absorbed and to remunerate the AS supplied ( p A S ). For downward regulations (i.e., load increase), the Aggregator pays p D A M p A S , the difference between the cost related to the energy absorbed but not declared in the DAM ( p D A M ) and the AS remuneration ( p A S ). Hence, the net revenues for the Aggregator from the AS provision ( R A S ) are calculated as a product between energy awarded in the ASM and its unitary reward ( p A S ). Assuming the economic neutrality of the TSO, the AS remuneration ( p A S ) has been set equal to c i m b , since, as general principle, the imbalance fees applied to users should reflect the costs covered by the system to purchase on the market the ASs required to correct the mismatches between expected and actual power exchanges [29].
The usage of each EV in the carsharing fleet is simulated considering the stochastic model described in Section 3. The statistical parameters adopted are based on the real data of carsharing usage in Milan collected during the weekdays of 2018, shown in Figure 2 and Table 1. EVs involved in the fleet are assumed to having P c h g i d ,   C b t r i d and e ¯ c o n s i d respectively equal to 7 kW, 40 kWh and 0.2 kWh/km [25]. These hypotheses are coherent with the scenario under analysis: a carsharing fleet composed, for technical and economic reasons, by low-cost EVs having similar characteristics.
The profitability of the CCA is evaluated by a holistic approach, taking into account also the possible detrimental effects on the carsharing service’s quality. According to data collected about the real carsharing usage in Milan, the distance covered between two charging stops is, in 99% of cases, lower than 100 km ( r n g m i n ). To ensure this autonomy, considering the technical characteristics previously presented, the EVs shall leave the CSs with a SoC ≥ 0.5. Assuming that the carsharing operator revenues are proportional to the distance traveled by each EV, it is possible to conclude that each time an EV with an autonomy ( r n g v ) lower than the minimum autonomy ( r n g m i n ) is rented, the Aggregator has a loss of profit c l o s t ( r n g m i n r n g v ) , where c l o s t is the operator’s net profit per kilometer covered by the customer. This fee is estimated by subtracting from the carsharing tariff c r e n t , set equal to 0.28 EUR/min [30], the operator’s costs (i.e., electric bill, taxes, maintenance costs, etc.). As gross estimate, c l o s t is set equal to 0.3· c r e n t . A similar approach has been applied to define c d l n in Equation (9). In Section 8, a sensitivity analysis on c r e n t   and c i m b is presented to evaluate the impact of these parameters on the Aggregator’s profits.

7. ABC/h-ABC Comparison

The ABC and h-ABC algorithms require a preliminary tuning to set the value of the optimization’s parameters (i.e.,   N F S ,   N m a i n , L i m M A X ), considering the trade-off between quality of the solution reached and time required to find it [27]. To this purpose, an upper limit of 60 s has been adopted for the computational time, since in real-life when an EV is connected to a CS, its charging can start only after that the optimal schedule has been identified by the CCA. The optimal parameters for the two algorithms, obtained in output to the tuning process, are shown in Table 2.
To test the effectiveness of the hybridizations proposed in the h-ABC, a comparison with the standard ABC approach is performed. To this purpose, the participation in the ASM of a carsharing fleet, composed of 1600 EVs, is simulated over a 3-month period. During the test, the performance of the two algorithms in finding the optimum is evaluated on almost 9000 optimization runs, using the parameters shown previously. For the comparison, in each time step, the value assumed by the objective function Equation (6) has been normalized with respect to the value obtained in the same time step at the end of the optimization process by the best performing algorithm (ABC or h-ABC). Finally, the average value of the objective function evaluated on all the populations for the two algorithms in each time interval has been calculated.
According to Figure 3, despite the lower number of initial populations ( N F S ), the h-ABC algorithm is capable to initialize solutions having objective function, on average, 98% lower than ABC method, thanks to the dispatching rules and the tournament selection procedure implemented. Additionally, regarding the value reached by the objective function for the optimal solution, the h-ABC outperforms the standard ABC, since the former is able to approach the optimal value much faster than the latter and, on average, the best value found by h-ABC is almost 50 times lower than with ABC. The effectiveness of h-ABC approach is also proved through a statistical analysis of the minimum value of the objective function reached in each instance by the two algorithms. Firstly, a Shapiro–Wilk test is performed to confirm the non-normality of the data, then, a paired Wilcoxon signed rank test is applied to evaluate whether the proposed approach performs better than the standard ABC. Adopting a Wilcoxon test with a confidence level of 1%, the resulting p-value is lower than 1 × 10−30. Hence, it is possible to conclude, with high statistical confidence, that the proposed methodology outperforms the standard ABC.
In order to better appreciate the complexity of the proposed charging scheduling and the number of variables involved during each optimization process, an example of charging profile optimized by h-ABC is shown in Figure 4. Each rectangle represents an EV in charge: therefore, the power absorbed in each time step is the total height of the rectangles, while the Power Request Schedule is the black dotted line. As one can observe, the power requested is almost always overlapped with the absorbed one, therefore, the imbalances experienced are nearly negligible. The black rectangles represent EVs not reaching the maximum SoC within the deadline set by the carsharing operator. It is important to point out that, even if these EVs are not fully charged, the SoC achieved is usually enough to cover the distance required by carsharing users. A further investigation on the performances obtainable by the h-ABC is provided in the next section.

8. Numerical Results

Once having proved the better performance of h-ABC approach compared to ABC, in this section, the techno-economic feasibility of the proposed scheduler is evaluated by simulating the behavior of a carsharing fleet with 1600 EVs over a 30-day period. Simulations are repeated 18 times, randomly changing the EVs’ usage and AS requests, therefore obtaining a total of 540 simulated days and 51.840 optimization runs. The scenario with the CCA (in the following, Sch scenario) is compared to the scenario currently in place for the e-mobility (No Sch scenario), in which EVs are charged as soon as they reach the CS, without neither charging scheduling nor provision of ASs to the grid. To make the No Sch and Sch scenarios fully comparable, the same EVs’ behavior is applied in both cases. Thus, the average energy absorbed is the same and equal to 7.94 MWh. The comparison is performed according to three aspects: (i) imbalances; (ii) AS supplied; (iii) carsharing service quality.
Concerning the first quantity, Figure 5 presents the distribution of the normalized daily imbalance, I m b % ( D ) , defined in the two scenarios using Equation (21):
I m b % ( D ) = t T D | P a b s ( t ) P r e q ( t ) | E a b s ( D )  
The numerator of Equation (21) represents the total imbalance in the day D, while the denominator, E a b s ( D ) , is the total daily energy absorbed by the aggregate. As shown in Figure 5, the CCA allows drastically cutting down the imbalance, obtaining   I m b % ( D ) < 1% in the 64% of simulated days and reducing the average imbalance from 17.3 to 1.3% (see dotted vertical lines in Figure 5). Figure 6 shows the comparison between the Power Schedule and the powers absorbed in the Sch and No Sch scenarios. The data refer to 2 days in August. Despite the uncertainties that characterize the EV’s charging requests, the proposed architecture is capable of avoiding power imbalances, exchanging with the grid the power declared in the DAM (Power Request Schedule).
The statistical distributions of simulated AS requests are presented in Figure 7 and Figure 8. On average, ASs are requested to the Aggregator in 40.6% of time steps and in 79.2% of cases these are downward requests. Figure 7 reports the PDF of the ratio between the requested power variation due to the AS and the average power absorbed ( P ¯ a b s ) by the aggregate (0.33 MW), while Figure 8 shows the PDF of the duration of dispatching orders (in hours). The ASs activated by the TSO during 3 days in August are depicted in Figure 9. It is possible to notice long upward activations (up to 8 consecutive hours) and large power variations (up to 26% of P ¯ a b s ).
To prove the capability of the proposed CCA to supply ASs, the amount of imbalance that occurred in time steps when the Aggregator is selected for AS provision ( t A S ) has been also analyzed. On average, in 85.4% of cases, a zero imbalance is registered during the AS provision (median value equal to 94%) and, for the remaining time, the imbalance is higher than 5% of the corresponding AS requests ( I m b ( t A S ) > 0.05   · P A S ( t A S ) ) only in 10% of cases, with a median value equal to 6%. This result has been achieved thanks to the implementation of the approach based on Equation (7), which allows to schedule the EVs’ charging also according to the AS requests expected in the upcoming time steps.
The effects of EVs’ scheduling on the quality of the carsharing service must be also evaluated, because an excessive postponement of the recharge could cause that the vehicles are rented before a full SoC is achieved, resulting in a lower autonomy. Table 3 reports the percentage of vehicles leaving the CS fully charged or with SoC ( S o C o u t ) lower than a given threshold. It can be noticed that the CCA affects almost negligibly the quality of the carsharing service offered, since the probability to rent a vehicle with SoC < 0.5 (corresponding to the minimum autonomy threshold presented in Section 5) is only marginally increased by the scheduling.
The loss of profit on day D for the carsharing operator can be calculated by Equation (22), where S D collects the EVs scheduled in the day D:
C r e n t ( D ) = c l o s t v v S D m a x ( 0 ; r n g m i n r n g v )  
Finally, it is possible to calculate the daily electricity bill ( C B i l l ) as the product between the total energy absorbed by the EVs ( E a b s ) and the energy price ( p D A M ). Concerning the imbalance cost ( C i m b ) and the revenues from the AS provision ( R A S ), they are computed as presented in Section 6.
Table 4 reports the economic outcomes for the two scenarios, normalized with respect to the mean daily total cost for the Aggregator obtained without the CCA ( C t o t N o s c h = C B i l l + C i m b ), on average equal to 454.5 EUR/day. In the No Sch scenario, the loss of profit due to the decrease in EV autonomy ( C r e n t ) and the AS remuneration are equal to zero. In the Sch scenario, the proposed architecture allows reducing the total daily cost by about 14.3%, thanks to the reduction of imbalance cost (56.1 EUR/day in the No Sch scenario and 4.14 EUR/day in the Sch one) and the AS revenues, roughly equal to 15 EUR/day (on average, −3.3% of C t o t N o s c h ). Moreover, analyzing the impact of the CCA on the carsharing service, it is possible to conclude that it is negligibly affected, because the corresponding costs ( C r e n t ) are increased only by 0.4%.
The tests are then repeated in two other case studies, characterized by different sizes of the carsharing fleet, respectively equal to 800 and 3200 EVs. The simulations were run on a computer equipped with an Intel Core i7-8700, 16 GiB of RAM, a Windows 10 operating system and MATLAB R2020a. The techno-economic profitability of h-ABC is once again evaluated, assuming as a strict requirement reaching the optimal scheduling in less than 60 s. On average, the h-ABC process took, respectively, in the two scenarios, 22.88 and 44.83 s to find the optimal schedule.
The objective function evaluation and the heuristic research (described in Algorithm 2) are the most impacting procedures in terms of both number of operations to be performed and computational times required. In the proposed h-ABC, similarly to other swarm-based optimization methods, the number of objective function evaluations tends to be high: on average 11,500 for each time step. However, this value does not grow considerably with the increase of the fleet size (10,020 in 800 EVs case study and 14,030 in 3200 EVs case study), while the heuristic research procedure is executed, on average, 3200 times for each time step, regardless of the case study analyzed. These results motivate the limited increment in computational times observed in the case studies analyzed in this paper, where the fleet size is increased by a factor of 4 (from 800 to 3200 EVs).
In the case characterized by a fleet of 3200 EVs, the higher number of optimization variables (on average 250 in each run) caused reaching the upper limit for the computational time (60 s) in 14.6% of the runs. In these runs, even if the convergence is not achieved by the h-ABC algorithm, the solution found can be considered very near to the optimal one: as highlighted by results reported in the next paragraph, obtained by considering also these suboptimal runs, the h-ABC performance in the 3200 EVs case study is still very promising.
In the No Sch scenario, the average normalized daily imbalance ( I m b % ) amounts, respectively, to 22.6% and 11.8% of the total absorbed energy, while with the CCA, they are reduced to 1.4% in both cases. Comparing the quality of the carsharing service offered to the users, it does not significantly change; the mean percentage of vehicles that leave the CSs with SoC < 0.5 is equal, in both cases, only to 0.17%. The economic outcomes are reported in Table 5, normalized with respect to the average daily cost obtained without the scheduler (237.2 EUR/day in the 800 EVs case study and 883.3 EUR/day in the 3200 EVs case study). The provision of AS allows saving, respectively, 3.1% and 3.6% of C t o t N o s c h . Additionally, imbalance costs, which respectively account for 16.0% and 9.7% of total costs in the scenario without CCA, are reduced to 0.5% and 1% by the scheduling. Despite a slight increase of the loss of profit ( C r e n t , % ) in both the scenarios equal to +0.3%, the scheduler permits a reduction of the total costs incurred by the carsharing operator of about 17.7% and 12.1% (i.e., 42.0 and 106.8 EUR/day).
Finally, in order to evaluate the profitability of the e-mobility scheduling under different economic conditions, a sensitivity analysis of the coefficients of the objective function is performed. In particular, both the values of c r e n t   and c i m b used in the objective function and in Equations (7) and (9) are varied in a range of ±25% with respect to the reference value previously applied (i.e., c r e n t = 0.28 EUR/min and c i m b = 40 EUR/MWh). In particular, the following steps are considered: [−25%; −15%; 0%; +15%; +25%], therefore obtaining a total of 16 different combinations. For the sake of the sensitivity analysis, a fleet composed of 1600 EVs is simulated over 200 days. The simulated days have been obtained by changing randomly the AS market requests and carsharing users’ behavior. Subsequently, an economic assessment is performed according to the procedure presented in Section 6, comparing the total costs for the carsharing operator in the Sch and No Sch scenarios: to this purpose, the quantity Δ C o s t % is calculated as the difference between the costs in the scheduling scenario and the base case, given by the electricity bill ( C B i l l ), the imbalances fee ( C i m b ), the loss of profit related to the worsening of carsharing service quality ( C r e n t ), minus the revenues from the ASM participation ( R A S ).
Figure 10 shows that the proposed h-ABC scheduling, in all the case studies, is able to effectively reduce the costs compared to the No Sch scenario, allowing for savings between 8.5% and 16% of the costs without the scheduler. It is worth noting that these results marginally depend on the cost of the carsharing rent ( c r e n t ); actually, as one can observe in Figure 10, assuming a given imbalance cost, the economic profitability of the h-ABC scheduler reduces only by about 2%, even when assuming a rent’s unitary cost increase of 50% (e.g., with c i m b equal to 40 EUR/MWh, by increasing the rent cost from 0.22 EUR/min to 0.34 EUR/min, the cost reduction changes from 13% to 11.2%). Another interesting result is related to the dependence of the economic profitability (i.e., Δ C o s t % ) on the unitary imbalance cost: an increase of c i m b allows the proposed CCA to be more profitable, even if the rent cost is kept unchanged. This fact, combined with the capability of the proposed algorithm to manage a large EV fleet, highlights how the h-ABC scheduling method could provide even more interesting results in the future, when the penetration of non-controllable renewable energy sources is expected to increase and, consequently, an increase of AS remuneration and imbalance costs is also expected, pushed by greater requests of the Ancillary Services needed to guarantee the power balance on the grid [31].

9. Conclusions

In the present paper, a hybridization of the ABC algorithm (h-ABC) has been presented to optimize the scheduling of the charging requests of an electric carsharing fleet. The centralized architecture, in perspective, will enable a carsharing operator to participate in the Ancillary Service Market by exploiting the flexibility of the EVs’ charging. This way the economic sustainability and competitiveness of the electric carsharing service will improve, without detrimental impacts on the quality of service provided to users. The numerical analyses performed highlighted that the h-ABC architecture is effective in reducing the power imbalance with respect to the commitments to the market, also enabling a reliable provision of Ancillary Services to the grid. Additionally, considering binding requirements for the computational time (a time limit of 60 s has been assumed), a cost reduction (electricity bill + imbalances) between 12.1% and 17.7% has been achieved, depending on the size of the EVs’ fleet and the forecasting error committed in estimating the DAM power schedule. A pivotal result is also the effectiveness shown by the method to preserve the quality of the carsharing service: a loss of profit, related to an unwanted decrease in the EVs’ autonomy caused by the scheduling, lower than or equal to 0.6%, has been estimated in the analyzed scenarios; therefore, it has been largely compensated by the benefits of the approach in terms of reduction of the amount of imbalance fee and the revenues from the AS provision.

Author Contributions

Conceptualization, D.F. and F.G.; methodology, D.F. and F.G.; software, F.G.; validation, D.F. and F.G.; resources, D.F.; data curation, F.G.; writing—original draft preparation, D.F. and F.G.; writing—review and editing, D.F. and F.G.; visualization, D.F. and F.G.; supervision, D.F.; project administration, D.F.; funding acquisition, D.F. All authors have read and agreed to the published version of the manuscript.

Funding

F. Gulotta’s position is partially funded by RSE S.p.A. within the National Research Fund for Electric Systems in compliance with the Decree of the Italian Minister of Economic Development on 16 April 2018.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Acronyms:
ABCArtificial Bee Colony
ASAncillary Service
ASM Ancillary Services Market
AT Arrival Time
BSP Balancing Service Provider
CCA Centralized Control Architecture
CS Charging Station
DAM Day-Ahead Market
DtT Due to Time
EV Electric Vehicle
h-ABCHybrid-Artificial Bee Colony
HP Hydro Power
SoC State of Charge
TSO Transmission System Operator
V 1 G Unidirectional Grid to Vehicle
V 2 G Bidirectional Vehicle to Grid
Parameters and Variables:
C b t r   Vehicle’s battery capacity
C B i l l Daily electricity bill cost
C i m b Daily imbalance cost
C r e n t Daily cost of non-fulfillment of the charging
C S Matrix containing the topology data of Charging Stations
c i m b Imbalance fee
c l o s t Fee in case of non-fulfillment of the charging
c r e n t Carsharing tariff for the users
D t i m e i d Charging deadline for the id-th vehicle
E a b s Total daily energy absorbed
e ¯ c o n s Mean energy consumption of the EV
f i t p Fitness function of the p-th solution
I m b % Normalized daily imbalance
L i m p Exploitation value of the p-th solution
L i m M A X Maximum exploitation limit
N d a y s t o t Number of simulated days
N e x p Number of expected EV connections in the following time steps
N f l e e t t o t Number of EVs in the fleet
N F S Number of solutions simultaneously evaluated by the optimization algorithm
N m a i n Maximum number of optimization iterations
n r e n t i d , d Number of rents during the d-th day for the id-th EV
P t Matrix containing the connected EVs’ information
P 0 Power Baseline Schedule
P A S Power accepted in the ASM
P a b s Power absorbed by the carsharing operator
P i c h g Charging power of the i-th vehicle
P r e q Power Request Schedule
p A S Remuneration for the AS provision
p D A M DAM price
R A S Daily remuneration from the AS provision
r n g m i n Minimum range required by the carsharing operator
S o C i n t r , i d / S o C C S , c o n r , i d State of Charge before/after the rent
s t i d Initial charging time for the id-th EV (optimization variable)
T D   Vector containing the time steps in day D
T f Upcoming time steps considered during the optimization process performed in t 0
t Generic time step
t A S Time steps with an AS request
t i n   Rent starting time
t 0 Current time step
t t o t Vector containing all the simulated time steps
t a r d Tardiness function
Δ C o s t   Cost variation w.r.t. the benchmark scenario
Δ t r e n t r Duration of the r-th rent
Δ d r e n t r Distance covered during r-th rent
Δ t c h g i d Charging time of the id-th EV
τ Time step resolution adopted

References

  1. Italian Energy Authority (ARERA). Res. 300/2017/R/eel. Initial Opening of the Ancillary Services Market to Electricity Demand and Production Units, Also from Renewables Not yet Enabled, as Well as Storage Systems. May 2017. Available online: https://www.arera.it/it/docs/17/300-17.htm (accessed on 20 March 2022).
  2. USA Federal Energy Regulatory Commission. Docket N. RM10-17-000; Order N. 745. Demand Response Compensation in Organized Wholesale Energy Markets. March 2011. Available online: https://www.ferc.gov/major-orders-regulations (accessed on 20 March 2022).
  3. AEMC. Demand Response Mechanism and Ancillary Services Unbundling, Final Rule Determination; AEMC: Sydney, Australia, 2016. [Google Scholar]
  4. Bulkeley, H.; Betsill, M. Rethinking Sustainable Cities: Multilevel Governance and the Urban’ Politics of Climate Change. Environ. Politics 2010, 14, 42–63. [Google Scholar] [CrossRef]
  5. Martin, E.W.; Shaheen, S.A. Greenhouse Gas Emission Impacts of Carsharing in North America. IEEE Trans. Intell. Transp. Syst. 2011, 12, 1074–1086. [Google Scholar] [CrossRef] [Green Version]
  6. Godina, R.; Rodrigues, E.M.G.; Matias, J.C.O.; Catalão, J.P.S. Smart EV charging scheduler for overloading prevention of an industry client power distribution transformer. Appl. Energy 2016, 178, 29–42. [Google Scholar] [CrossRef]
  7. Gulotta, F.; Rancilio, G.; Blaco, A.; Bovera, F.; Merlo, M.; Moncecchi, M.; Falabretti, D. E-Mobility Scheduling for the Provision of Ancillary Services to the Power System. Int. J. Electr. Electron. Eng. Telecommun. 2020, 9, 349–355. [Google Scholar] [CrossRef]
  8. Mukherjee, J.C.; Gupta, A. A Review of Charge Scheduling of Electric Vehicles in Smart Grid. IEEE Syst. J. 2015, 9, 1541–1553. [Google Scholar] [CrossRef]
  9. Mozafar, M.R.; Amini, M.H.; Moradi, M.H. Innovative appraisement of smart grid operation considering large-scale integration of electric vehicles enabling V2G and G2V systems. Electr. Power Syst. Res. 2018, 154, 245–256. [Google Scholar] [CrossRef]
  10. Uddin, K.; Dubarry, M.; Glick, M.B. The viability of vehicle-to-grid operations from a battery technology and policy perspective. Energy Policy 2018, 113, 342–347. [Google Scholar] [CrossRef]
  11. Kuran, M.Ş.; Viana, A.C.; Iannone, L.; Kofman, D.; Mermoud, G.; Vasseur, J.P. A Smart Parking Lot Management System for Scheduling the Recharging of Electric Vehicles. IEEE Trans. Smart Grid 2015, 6, 2942–2953. [Google Scholar] [CrossRef]
  12. Falabretti, D.; Gulotta, F.; Spinelli, L. Participation of Aggregated DERs to the Ancillary Services Market: A Monte Carlo Simulation-Based Heuristic Greedy-Indexing Model. Energies 2022, 15, 1037. [Google Scholar] [CrossRef]
  13. Al-Ogaili, A.S.; Tengku Hashim, T.J.; Rahmat, N.A.; Ramasamy, A.K.; Marsadek, M.B.; Faisal, M.; Hannan, M.A. Review on Scheduling, Clustering, and Forecasting Strategies for Controlling Electric Vehicle Charging: Challenges and Recommendations. IEEE Access 2019, 7, 128353–128371. [Google Scholar] [CrossRef]
  14. Álvarez, J.G.; González, M.; Vela, C.R.; Varela, R. Electric Vehicle Charging Scheduling by an Enhanced Artificial Bee Colony Algorithm. Energies 2018, 11, 2752. [Google Scholar] [CrossRef] [Green Version]
  15. Hu, Z.; Zhan, K.; Zhang, H.; Song, Y. Pricing mechanisms design for guiding electric vehicle charging to fill load valley. Appl. Energy 2016, 178, 155–163. [Google Scholar] [CrossRef]
  16. Li, S.; Li, J.; Su, C.; Yang, Q. Optimization of Bi-Directional V2G Behavior With Active Battery Anti-Aging Scheduling. IEEE Access 2020, 8, 11186–11196. [Google Scholar] [CrossRef]
  17. Rasheed, M.B.; Awais, M.; Alquthami, T.; Khan, I. An Optimal Scheduling and Distributed Pricing Mechanism for Multi-Region Electric Vehicle Charging in Smart Grid. IEEE Access 2020, 8, 40298–40312. [Google Scholar] [CrossRef]
  18. Falabretti, D.; Gulotta, F. An Algorithm for the Ancillary Services Provision by E-Mobility-based Virtually Aggregated Mixed Units. In Proceedings of the 2020 IEEE International Conference on Environment and Electrical Engineering (EEEIC), Madrid, Spain, 9–12 June 2020; pp. 1–6. [Google Scholar]
  19. Arık, O.A. Comparisons of metaheuristic algorithms for unrelated parallel machine weighted earliness/tardiness scheduling problems. Evol. Intell. 2020, 13, 415–425. [Google Scholar] [CrossRef]
  20. Muthiah, A.; Rajkumar, A.; Rajkumar, R. Hybridization of Artificial Bee Colony Algorithm with Particle Swarm Optimization algorithm for Flexible Job Shop Scheduling. In Proceedings of the 2016 International Conference on Energy Efficient Technologies for Sustainability (ICEETS), Nagercoil, India, 7–8 April 2016; pp. 896–903. [Google Scholar]
  21. Shah, H.; Tairan, N.; Garg, H.; Ghazali, R. Global Gbest Guided-Artificial Bee Colony Algorithm for Numerical Function Optimization. Computers 2018, 7, 69. [Google Scholar] [CrossRef] [Green Version]
  22. Pilo, F.; Pisano, G.; Ruggeri, S.; Soma, G.G.; Falabretti, D.; Grillo, S.; Gulotta, F. Impact of e-mobility participation in the ancillary service market on the operation of high-density urban low voltage distribution networks. In Proceedings of the 2021 AEIT International Conference on Electrical and Electronic Technologies for Automotive (AEIT AUTOMOTIVE), Torino, Italy, 17–19 November 2021. [Google Scholar] [CrossRef]
  23. European Commission, Directorate-General for Energy. Commission Regulation (EU) 2017/2195 of 23 November 2017. Establishing a Guideline on Electricity Balancing. Available online: https://eur-lex.europa.eu/legal-content/EN/TXT/?uri=CELEX%3A02017R2195-20210315 (accessed on 20 March 2022).
  24. Chow, J.H.; Wu, F.F.; Momoh, J.A. Applied Mathematics for Restructured Electric Power Systems: Optimization, Control, and Computational Intelligence, 1st ed.; Springer: New York, NY, USA, 2004. [Google Scholar]
  25. Benetti, G.; Delfanti, M.; Facchinetti, T.; Falabretti, D.; Merlo, M. Real-Time Modeling and Control of Electric Vehicles Charging Processes. IEEE Trans. Smart Grid 2015, 6, 1375–1385. [Google Scholar] [CrossRef]
  26. Karaboga, D.; Basturk, B. On the performance of artificial bee colony (ABC) algorithm. Appl. Soft Comput. 2008, 8, 687–697. [Google Scholar] [CrossRef]
  27. Álvarez, J.G.; González, M.A.; Vela, C.R. A Genetic Algorithm for Scheduling Electric Vehicle Charging. In Proceedings of the 2015 Annual Conference on Genetic and Evolutionary Computation, New York, NY, USA, 11 July 2015; pp. 393–400. [Google Scholar]
  28. Akay, B.; Karaboga, D. Parameter Tuning for the Artificial Bee Colony Algorithm. In Proceedings of the 1st International Conference on Computational Collective Intelligence (ICCCI 2009), Wroclaw, Poland, 5–7 October 2009. [Google Scholar]
  29. Agency for the Cooperation of Energy Regulators (ACER). Annual Report on the Results of Monitoring the Internal Electricity Markets in 2015. September 2016. Available online: https://documents.acer.europa.eu/Official_documents/Publications# (accessed on 20 March 2022).
  30. Ampudia-Renuncio, M.; Guirao, B.; Molina-Sanchez, R.; Bragança, L. Electric Free-Floating Carsharing for Sustainable Cities: Characterization of Frequent Trip Profiles Using Acquired Rental Data. Sustainability 2020, 12, 1248. [Google Scholar] [CrossRef] [Green Version]
  31. Seel, J.; Mills, A.D.; Wilser, R.H. Impacts of High Variable Renewable Energy Futures on Wholesale Electricity Prices, and on Electric-Sector Decision Making. Lawrence Berkeley National Laboratory, 2018. Available online: https://escholarship.org/uc/item/2xq5d6c9 (accessed on 20 March 2022).
Figure 1. Scenario of application of the proposed scheduler.
Figure 1. Scenario of application of the proposed scheduler.
Energies 15 03023 g001
Figure 2. Hourly probability of carsharing rent requests.
Figure 2. Hourly probability of carsharing rent requests.
Energies 15 03023 g002
Figure 3. Mean value of the objective function in ABC and h-ABC as a function of computational time.
Figure 3. Mean value of the objective function in ABC and h-ABC as a function of computational time.
Energies 15 03023 g003
Figure 4. Charging profile obtained in the 3rd day simulated by h-ABC.
Figure 4. Charging profile obtained in the 3rd day simulated by h-ABC.
Energies 15 03023 g004
Figure 5. Comparison of the imbalance PDF in the two scenarios.
Figure 5. Comparison of the imbalance PDF in the two scenarios.
Energies 15 03023 g005
Figure 6. Comparison between the power schedule and the power absorbed, in Sch and No Sch scenarios.
Figure 6. Comparison between the power schedule and the power absorbed, in Sch and No Sch scenarios.
Energies 15 03023 g006
Figure 7. PDF of the power awarded for AS request with respect to P ¯ a b s .
Figure 7. PDF of the power awarded for AS request with respect to P ¯ a b s .
Energies 15 03023 g007
Figure 8. PDF of the AS duration requests in hours.
Figure 8. PDF of the AS duration requests in hours.
Energies 15 03023 g008
Figure 9. Trends of the AS requests during 3 days of examples.
Figure 9. Trends of the AS requests during 3 days of examples.
Energies 15 03023 g009
Figure 10. Map of cost reduction of the proposed scheduler, normalized w.r.t the No Sch total cost, for different coefficient combinations.
Figure 10. Map of cost reduction of the proposed scheduler, normalized w.r.t the No Sch total cost, for different coefficient combinations.
Energies 15 03023 g010
Table 1. Statistics adopted to simulate the carsharing usage.
Table 1. Statistics adopted to simulate the carsharing usage.
Variable ModeledDistribution and Parameters Adopted
Rent   duration   ( Δ t r e n t ) (min)Gamma distribution
k = 2.98   θ = 5.51
Rent   turnover   ( n r e n t   ) (rents/day)Integer normal distribution
μ = 4     σ = 1
Distance   covered   ( Δ d e v ) (km)Log-logistic
μ = 1.49     θ = 0.43
Table 2. Optimal hyper-parameters found for the two algorithms.
Table 2. Optimal hyper-parameters found for the two algorithms.
Algorithm ParameterOptimal Value Found
ABCh-ABC
N °   main   algorithm   iteration   ( N m a i n )200260
N °   candidate   solution   ( N F S )3022
Scout   bee   limit   ( L i m m a x )150100
Table 3. Cumulated PDF of EVs’ SoC when disconnected from a CS.
Table 3. Cumulated PDF of EVs’ SoC when disconnected from a CS.
Scenario P   ( S o C o u t =   1 ) P   ( S o C o u t <   0.9 ) P   ( S o C o u t <   0.8 ) P   ( S o C o u t <   0.7 ) P   ( S o C o u t <   0.6 ) P   ( S o C o u t <   0.5 ) P   ( S o C o u t <   0.4 )
No Sch90.1%5.4%2.4%1.1%0.5%0%0%
Sch85.0%6.5%2.9%1.3%0.6%0.04%0.02%
Table 4. Economic results of the No Sch and Sch scenarios (in % w.r.t. the total cost of the No Sch scenario; negative values are actually incomes).
Table 4. Economic results of the No Sch and Sch scenarios (in % w.r.t. the total cost of the No Sch scenario; negative values are actually incomes).
C B i l l , %   C i m b , % R A S , % C r e n t , % Δ C o s t %
No Schmean (%) (Max; std)87.7%
(98.5; 14.0)
12.3%
(55.7; 7.3)
-0%
(0; 0)
100%
SchMean (%) (Max; std)87.7%
(97.5; 10.1)
0.9%
(7.6; 1.3)
3.3%
(11.5; 2.8)
0.4%
(4.0; 0.6)
85.7%
Table 5. Economic profitability of the proposed CCA as a function of EV fleet size.
Table 5. Economic profitability of the proposed CCA as a function of EV fleet size.
800 EVs3200 EVs
No SchSch Δ No SchSch Δ
C B i l l , % (%)84.083.9−0.190.390.1−0.2
C i m b , % (%)16.00.9−15.19.71.1−8.6
R A S , % (%)-−3.1−3.1-−3.6−3.6
C r e n t , % (%)00.6+0.600.3+0.3
Δ C o s t % (%)10082.3−17.710087.9−12.1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Falabretti, D.; Gulotta, F. A Nature-Inspired Algorithm to Enable the E-Mobility Participation in the Ancillary Service Market. Energies 2022, 15, 3023. https://doi.org/10.3390/en15093023

AMA Style

Falabretti D, Gulotta F. A Nature-Inspired Algorithm to Enable the E-Mobility Participation in the Ancillary Service Market. Energies. 2022; 15(9):3023. https://doi.org/10.3390/en15093023

Chicago/Turabian Style

Falabretti, Davide, and Francesco Gulotta. 2022. "A Nature-Inspired Algorithm to Enable the E-Mobility Participation in the Ancillary Service Market" Energies 15, no. 9: 3023. https://doi.org/10.3390/en15093023

APA Style

Falabretti, D., & Gulotta, F. (2022). A Nature-Inspired Algorithm to Enable the E-Mobility Participation in the Ancillary Service Market. Energies, 15(9), 3023. https://doi.org/10.3390/en15093023

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

Article Metrics

Back to TopTop