Next Article in Journal
A Study on Urban Ethnic Segmentation in Kabul City, Afghanistan
Previous Article in Journal
The Dynamic Correlation and Volatility Spillover among Green Bonds, Clean Energy Stock, and Fossil Fuel Market
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Planning of Electric Vehicle Fast-Charging Stations Considering Uncertain Charging Demands via Dantzig–Wolfe Decomposition

1
College of Economics and Management, Southwest University, Chongqing 400715, China
2
School of Economics, Chongqing Financial and Economic College, Chongqing 401320, China
3
College of Mathematics and Statistics, Chongqing Jiaotong University, Chongqing 400074, China
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(8), 6588; https://0-doi-org.brum.beds.ac.uk/10.3390/su15086588
Submission received: 19 March 2023 / Revised: 7 April 2023 / Accepted: 10 April 2023 / Published: 13 April 2023

Abstract

:
This study investigates the planning problem of fast-charging stations for electric vehicles with the consideration of uncertain charging demands. This research aims to determine where to build fast-charging stations and how many charging piles to be installed in each fast-charging station. Based on the multicommodity flow model, a chance-constrained programming model is established to address this planning problem. A scenario-based approach as well as a big-M coefficients generation algorithm are applied to reformulate the programming model into tractable one, then the Dantzig–Wolfe decomposition method is leveraged to find its optimal solution. Finally, a numerical experiment is conducted in a 25-node network to assess the efficiency of the proposed model and solution approach.

1. Introduction

Electric vehicles (EVs) are widely acknowledged as a promising solution for reducing greenhouse gas emissions from road transportation and reducing reliance on finite fossil energy sources. Consequently, numerous governments have implemented incentive policies to encourage widespread adoption of EVs, resulting in significant progress towards this goal. However, there are several major obstacles that hinder the adoption of EVs on a larger scale. Among them, the “range anxiety” may be the most momentous challenge, which reflects the fear of running out of energy before reaching the destination or charging station. In practice, an efficient way to alleviate such negative psychology for popularizing EVs is to reasonably plan fast-charging stations in transportation networks. This situation motivates us to investigate the optimal planning problem of fast-charging stations, that is, answer the questions of where to build fast-charging stations and how many charging piles are required in one fast-charging station [1,2,3].
The planning of fast-charging stations has been studied in various contexts, with the charging demands of EVs being a critical component. In this paper, we address the challenge of uncertain charging demands, for which probability distributions are unknown. Instead, we rely on large amounts of collected data on charging demands and adopt a data-driven approach to formulate a siting and sizing model for fast-charging stations. Our approach involves formulating data-driven chance constraints for planning decisions, by incorporating sample data on uncertain charging demands. Specifically, we aim to determine the appropriate number of charging piles to be installed at fast-charging stations, such that they can accommodate the allocated charging demands while meeting a prescribed service level.
Chance constraints are known to be intractable due to their general nonconvexity [4]. To overcome this computational difficulty, we adopt a scenario-based reformulation approach based on a finite number of scenarios extracted from samples of the charging demands. By considering these scenarios, we reformulate the chance-constrained program into a mixed-integer linear program (MILP) that can be efficiently solved using Dantzig–Wolfe decomposition. However, in practice, the MILP is often weak due to the introduction of big-M coefficients through the scenario-based approach. To address this, we propose a big-M coefficient strengthening method to identify appropriate big-M coefficients that result in better computational performance for the subproblems.
The remainder of this paper is organized as follows. Section 2 reviews the literature on the chance-constrained siting and sizing problem of fast-charging stations and highlights the main contributions of this paper. Section 2.5 formulates the chance-constrained programming model for such problem. Section 3 presents a mixed-integer linear reformulation using scenario-based approach and Dantzig–Wolfe decomposition, and provides an efficient way to strengthen the big-M coefficients that can be exploited to accelerate the computation. Section 4 conducts numerical experiments and illustrates managerial insights in the context of both virtual and real world transportation networks. Section 5 concludes this paper with a summary of the most important findings and opportunities for future work.

2. Literature Review

This section reviews the literature on the fast-charging station planning problem and on Dantzig–Wolfe decomposition. We close this section with a statement of the main contributions of this paper.

2.1. Fast-Charging Station Planning Problem

The optimization models that tackle the refueling station planning problem are generally extended from the facility location problem. There are two dominant modeling approaches: the node-based approach and the flow-based approach [5]. The node-based approach considers the node-based demands and utilizes the generalization of p-median model, whose goal is to maximally cover the potential charging demands from a given area by planning a minimum number of fast-charging stations [6,7,8,9,10]. However, the limited driving ranges of EVs make the complete coverage impossible, since EVs may require to be charged multiple times to sustain their long trips, and there always exist unexpected charging demands from other areas. In this regard, the flow-based approach provides an alternative modeling way which considers the charging demands on the path between O-D pairs. The flow-based modeling idea originates from [11], who proposes a flow-capturing location model (FCLM) which aims to plan a given number of stations to capture as much flow as possible. To account for the limited driving ranges, [12] proposes the flow refueling location model (FRLM), where a sequence of refueling stations along a path between an O-D pair that enables a successful trip is defined with respect to a feasible combination of refueling facilities. The FRLM has been extended to be applicable for various theoretical and engineering aspects ([13] and the references therein).
The aforementioned works merely consider the deterministic charging demands; however, the charging demands usually encounter the uncertain fluctuation in reality, which brings about difficulties in obtaining accurate charging demands. The pioneering work that handles the stochastic charging demands is reported by [14], in which the charging demands are assumed to follow Poisson distribution. They develop a cost-concern model as well as a goal-driven model by the distributionally robust optimization approach, where the constraints that relate to the worst-case expected cost and the worst-case occurring probability are approximated by the standard normal distribution. Following this work, Refs. [15,16] consecutively consider the fast-charging station planning problem with the consideration of Poisson-distribution-based charging demands. Such two works simultaneously formulate a second-order cone programming planning model, where the differences appear in the method of solving: Ref. [15] directly uses the CPLEX solver, while Ref. [16] leverages Benders decomposition as well as the CPLEX solver.
In practice, the distribution of uncertain charging demands is often unknown, which has posed challenges to the stochastic modeling approach [17,18,19]. Instead, massive amounts of data are routinely available in transportation networks, for which the distributionally robust optimization and chance-constrained programming approach is more suitable. Ref. [20] formulates a two-stage distributionally robust optimization planning model that incorporates the uncertain charging demands, which is a generalization of FCLM. They employ a sample average approximation method and an averaged two-replication procedure to solve the planning model. Ref. [21] formulates a distributionally robust optimization planning model based on ϕ -divergence, where the potential distributions of uncertain charging demands are included in a ϕ -divergence-based ambiguous set that is constructed using the available data of charging demands. In order to be efficiently solved via off-the-shelf commercial solvers, they reformulate the planning model into its tractable counterpart using standard reformulation techniques of distributionally robust programs. In addition, Ref. [22] formulates a bi-level chance-constrained programming (CCP) model to tackle the uncertain charging demands, in which the chance constraints guarantee that the probability the EVs successfully get charged in a fast-charging station at least reaches a specified service level. The chance-constrained planning model is finally reformulated into the second-order cone counterparts that can be efficiently solved by off-the-shelf commercial solvers.

2.2. Dantzig–Wolfe Decomposition

Although the planning models have concise formulations, the massive number of decision variables related to large-scale transportation networks makes the solving process computationally demanding, and considering an arbitrary subset of them may bring about a suboptimal solution [23,24,25]. Thereby, the column generation method is introduced to screen out the variables that do not have the potential to consist the optimal solution. Ref. [26] considers the refueling station location problem with routing. The authors apply a branch-and-price algorithm, in which the column generation procedure is employed to manage the size of the location model so as to be solved efficiently. Ref. [27] investigates the optimal deployment problem of charging stations, in which a branch-and-price algorithm is devised for the nonlinear deployment model whose linear relaxations are solved by column generation. Ref. [28] develops a Benders-and-Price algorithm to solve the developed charging station location model, where the column generation procedure is proposed to solve the Benders subproblem. It should be pointed out that Benders decomposition is not applicable for the scenario-based reformulation of chance constraints, because this approach introduces the auxiliary binary variables and big-M coefficient, which are obstacles for obtaining Benders cuts through the duality method.
In this aspect, Dantzig–Wolfe decomposition is applicable, which also employs a column generation procedure [29,30]. Different from Benders decomposition that adds new constraints to the Benders master problem according to the solution of the Benders subproblem, Dantzig–Wolfe decomposition adds new variables to Dantzig–Wolfe master problem according to Dantzig–Wolfe subproblems, which are column generation (pricing) subproblems that aim to find the variables with positive/negative reduced costs [31]. Briefly speaking, Benders decomposition adds rows to the master problem, while Dantzig–Wolfe decomposition adds columns to the master problem. Ref. [32] adapts Dantzig–Wolfe decomposition to the multiple depot vehicle scheduling problem. They apply two types of Dantzig–Wolfe decomposition to the inventory formulation: one for each depot, the other for each trip. Ref. [33] considers the facility location and production planning problem, in which Dantzig–Wolfe decomposition is proposed to detect and compare the lower bounds between the two formulated mathematical models. Ref. [34] describes two versions of the chance-constrained bin packing problem, the first one is a scenario-based reformulation, and the second one is a distributionally robust reformulation. Both versions are solved by a Dantzig–Wolfe formulation suited to a branch-and-price algorithm. Ref. [35] utilizes Dantzig–Wolfe decomposition as well as a branch-and-price procedure to solve an MILP model for the vehicle allocation problem.

2.3. Summary of Main Contributions

When implementing the standard Dantzig–Wolfe decomposition, one master problem and multiple subproblems should be formulated. The master problem should gather the optimal solutions (extreme points) of all the subproblems to represent its feasible solution as the convex combination of such extreme points, which may give rise to unexpected demanding computation in a large-scale network (so-called “dimensional curse”), because massive decision variables have to be introduced corresponding to the extreme points of the subproblems. To overcome this “dimensional curse”, we intend to separate the original planning problem into disjoint sub-components corresponding to the O-D pairs, then leverage a Dantzig–Wolfe decomposition to find their optimal solutions, in which the master problem and subproblem are one-to-one with respect to each O-D pair. Thus, a moderate scale of decision variables are involved in each pair of master problem and subproblem, which is conducive to improving the computational performance of Dantzig–Wolfe decomposition. As long as all the sub-components are solved, we are able to obtain the optimal solution of the original problem by synthesizing the optimal solution of all the sub-components. In addition, the disjoint property of the sub-components makes their facilitation to be implemented in a decentralized manner using parallel machines, which helps improve the computational efficiency as well. To the best of our knowledge, this is the first attempt to leverage Dantzig–Wolfe decomposition for the optimal planning of electric vehicle fast-charging stations considering uncertain charging demands. We contribute to existing works as follows:
(i)
The data-driven chance constraints are introduced to address the uncertain charging demands, since it is difficult to extract the exact values and accurate probability distributions of the uncertain charging demands in practice. In this circumstance, date-driven chance constraints are viable alternatives to measure the service abilities of fast-charging stations.
(ii)
The original planning model is separated into multiple disjoint sub-components corresponding to the O-D pairs in the transportation network, which are facilitated to be implemented in a decentralized manner using parallel machines. Then, the Dantzig–Wolfe decomposition is utilized to solve the planning model, where the master problems and subproblems are one-to-one. Each possesses a moderate scale of decision variables, which makes it conducive to improve the computational performances.
(iii)
The computational performance is tested using the 25-node benchmark network, in which the proposed Dantzig–Wolfe decomposition is shown to be efficient for solving the planning problem over a large transportation network.

2.4. Expanded Network Generation

Prior to establishing the planning model, a finite set of candidate siting nodes should be identified according to the driving ranges of EVs, a part of which is then chosen to build fast-charging stations and install charging piles so as to satisfy the potential charging demands of EVs. Here, we use the expanded network approach. The expanded network approach is originally proposed in [36], where the charging constraints that express the charging demands on the paths are presented by the flow balance constraints, that is, an available flow on the expanded network ensures a feasible charging pattern. Although the expanded network expands the scale of the original network, its numerical efficiency is demonstrated to outperform some existing models, and simultaneously, it removes the big-M constants from the formulations that yield poor linear relaxation bounds [28]. Due to the above-mentioned advantages, the expanded network model is widely employed and extended in the literature, such as in [28,37,38,39]. On the other hand, another modeling idea called sub-path is also proposed in [14] and employed and extended in [15,16,26,27]. As the name suggests, the sub-path model segments a complete path into several sub-paths with the consideration of driving ranges, then the optimization model for siting and sizing the fast-charging stations is formulated based on such sub-paths. Owing to the fact that the paths in the original transportation network are segmented into small pieces, the flow balance constraints that are considered in the expanded network model cannot be incorporated. Therefore, the sub-path model often considers the charging demands on each sub-path or individual nodes instead. In practice, however, the transportation networks are usually composed by a large amount of nodes that may result in countless sub-paths, which leads to computational difficulties. To avoid this situation, we employ the expanded network model, which is generated based on the following assumptions:
(A1)
More than half of the battery is full when travel starts at the original node.
(A2)
More than half of the battery is full when travel ends at the destination node.
(A3)
The driving ranges of electric vehicles are deterministic and fixed.
When finishing their previous trips, the drivers of electric vehicles always have the willingness to reserve enough remaining mileage for their next trips. Owing to the driving range anxiety of electric vehicle drivers, there exists a threshold of remaining mileage for the drivers: If the remaining mileage stage is above the threshold, the drivers would continue their next trip; otherwise, they would choose to charge in fast-charging stations. Assumption (A1) describes the situation in which the threshold of remaining mileage is set to be half of the full driving range. Assumption (A2) describes another situation that, when a driver needs to complete a round trip, they should reserve adequate remaining mileage after finishing a single-loop trip. Assumption (A3) indicates that the driving ranges are always a constant R that describes the driving distance of a fully charged electric vehicle.
Let G V , A , Q be a given transportation network, where V and A are the set of nodes and arcs, respectively. We consider an O-D pair q on G , which starts at the origin node O and terminates at the destination node D. We denote the set of nodes and arcs on O-D pair q as V q and A q , respectively, and the driving distance on arc ( i , j ) A q as d i j q . Then, the generation of the expanded network G ¯ q V ¯ q , A ¯ q based on O-D pair q involves the following steps:
Step 1 
Before the original node O, add a source node s q . After the destination node D, add a sink node t q . Note that V ¯ q is V ¯ q V q { s q , t q } .
Step 2 
Connect s q and O, D and t q . Note that A ¯ q A q { s q , O , D , t q } .
Step 3 
Ordering the nodes in V ¯ q sequentially. Let o r d ( · ) denote the ordering index of the nodes.
Step 4 
Add arc ( s q , j ) to A ¯ q if d O j q is strictly less than R 2 ; add arc ( i , t q ) to A ¯ q if d i D q is strictly less than R 2 . Add arc ( i , j ) to A ¯ q if o r d ( i ) < o r d ( j ) and d i j q is strictly less than R.
We take a toy network in Figure 1, Figure 2 and Figure 3 as an example to clarify how to generate the expanded network according to Step 1Step 4, where the driving range R is set to be 120 km.

2.5. Mathematical Model

Based on the expanded network, an optimization model for the fast charging station planning problem is proposed, which is based upon a generalized multicommodity flow model. The conventional multicommodity flow model considered in [28] utilizes binary variables for the arcs to represent whether the arcs are used by the EVs, while we utilize continuous variables for the arcs to represent the proportion of charging demands that allocated on the arcs. In practice, due to the rapid growth of electric vehicle ownership, their charging demands have emerged a salient feature of uncertainty; even their probability distributions are difficult to be accurately obtained. Thus, we consider the uncertain charging demands F ˜ q on path q Q . For any path q Q , we generate the corresponding expanded network G ¯ q V ¯ q , A ¯ q . Let c 1 , i and c 2 , i be the building cost of a fast-charging station and installing cost of charging piles at node i V . Then, the mathematical model for the fast-charging station planning problem is formulated as follows:
(1a) ( P ) min i V c 1 , i u i + c 2 , i y i , (1b) s . t . j N i x i j q j N i + x j i q = 1 , if i = s q , 1 , if i = t q , 0 , if i s q , i t q , q Q , (1c) P q Q F ˜ q j N i + x j i q y i α , i V , (1d) y i Υ u i , i V , (1e) u i { 0 , 1 } , i V , (1f) y i Z + , i V , (1g) x i j q 0 , ( i , j ) A ¯ q , q Q .
where N i + and N i are the set of in-neighbors and out-neighbors of node i on the expanded network of O-D pair q, that is, N i + = j | ( j , i ) A ¯ q , q Q and N i = j | ( i , j ) A ¯ q , q Q . Objective function (1a) aims to minimize the total setup cost of fast-charging stations, which includes building the fast-charging stations and installing the charging piles. Constraints (1b) are the flow balance constraints that allocate the charging demands to every node that the path traverses, which describe the process that the charging demands flow losslessly from the source node to the sink node through the intermediate nodes. Constraints (1c) are the chance constraints that ensure that the allocated charging demands are accommodated by building fast-charging stations and installing charging piles that at least reach a specified service level α . Constraints (1d) enforce the upper bound of installed charging piles at a fast-charging station. Constraints (1e)–(1g) restrict the domains of the siting, sizing, and allocation decision variables, respectively. Figure 4 illustrates a toy network to clarify the generalized multicommodity flow model. The flow conservation constraints (1b) ensures that the out-flow and in-flow at each node are balanced, that is, x O 1 = x 13 + x 41 = 1 2 , x O 2 = x 2 D = 1 2 , x 13 = x 3 D = 1 4 , x 14 = x 4 D = 1 4 , which indicates that 1 2 , 1 2 , 1 4 , and 1 4 of total charging demands are allocated at Node 1–Node 4, respectively. Thus, fast-charging stations should be planned properly at these nodes to satisfy such allocated charging demands.
Even though the number of charging piles installed in a fast-charging station has to be an integer, it can be relaxed as a continuous variable without affecting the optimality of P. The following theorem supports this claim.
Theorem 1. 
Denote RP the relaxed model of P, where variables y i in model P are relaxed to be continuous variables v i in RP for all i V . Let u i * , v i * i V , x j i q * q Q , ( j , i ) A ¯ q be the optimal solution of  RP. Then, u i * , y i * i V , x j i q * q Q , ( j , i ) A ¯ q with y i * = v i * is an optimal solution of  P, where · stands for the ceiling function of real numbers.
Proof. 
We consider the following optimization model
(2a) ( P ˜ ) min i V c 1 , i u i , s . t . ( 1 b ) , ( 1 e ) , ( 1 g ) , (2b) P q Q F ˜ q j N i + x j i q Υ u i α , i V .
Suppose the optimal solution of P ˜ to be u i * i V , x j i q * q Q , ( j , i ) A ¯ q . Notice that we are capable of equivalently transforming Constraints (2b) into Constraints (1c) and (1d) by introducing the sizing variables y i for all i V . Therefore, the optimal values of variables u i and x i j q in RP and P must be identical to that in model ( P ˜ ). Thus, we suppose the optimal solution of RP and P to be u i * , v i * i V , x j i q * q Q , ( j , i ) A ¯ q and u i * , y i * i V , x j i q * q Q , ( j , i ) A ¯ q , respectively. Since the sizing variables are restricted to be integers according to Constraints (1f), in order to minimize the objective (1a), y i * has to be the lowest integer not less than v i * , that is, y i * = v i * for all i V .   □
Hereafter, we consider the relaxed model RP instead of P
(3a) ( RP ) min i V c 1 , i u i + c 2 , i v i , s . t . ( 1 b ) , ( 1 e ) , ( 1 g ) , (3b) P q Q F ˜ q j N i + x j i q v i α , i V , (3c) v i Υ u i , i V , (3d) v i 0 , i V .
Although RP provides a concise formulation of planning the fast-charging stations, the nonconvexity of chance constraints as well as the tremendous candidate nodes give rise to its intractability. In the next section, we initially reformulate the chance constraints (3b) into their tractable counterparts, then propose a Dantzig–Wolfe-decomposition-based approach to solve the reformulations to optimality.

3. Scenario-Based Approach

Solving optimization model RP requires tackling chance constraints (3b) prior to this. Here, we leverage the scenario-based technique to do such work. The Abbreviations section lists the nomenclatures of the used symbols and notations of this section.

3.1. Scenario-Based Reformulation

Denote Ω as the finite set of all possible scenarios based on the sample average approximation (SAA), where ω indexes the realizations of charging demands F ˜ q for all q Q as F ^ 1 q , F ^ 2 q , , F ^ | Ω | q . With the consideration of scenarios of charging demands, we have to further introduce auxiliary binary variables z i ω to reformulate Constraints (3b) into their linear counterparts, where z i ω = 1 if fast-charging station at node i provides adequate service abilities for the allocated charging demands in scenario ω , that is, the ith chance constraint is satisfied, and z i ω = 0 otherwise. Then, a tractable reformulation of RP can be expressed as follows
(4a) ( SRP ) min i V c 1 , i u i + c 2 , i v i , s . t . ( 1 b ) , ( 1 e ) , ( 1 g ) , ( 3 c ) , ( 3 d ) , (4b) q Q F ω q j N i + x j i q v i + M i ω 1 z i ω , i V , ω Ω , (4c) ω Ω z i ω α | Ω | , i V , (4d) z i ω { 0 , 1 } , i V , ω Ω ,
where M i ω are big-M coefficients for all i V , ω Ω . Constraints (4b) and (4c) linearize Constraints (3b), which jointly enforce the service abilities of fast-charging stations.
The big-M coefficients M i ω in Constraints (4b) have to be not less than q Q F ω q j N i + x j i q Υ for all ω Ω and i V , so as to ensure the feasibility of SRP-SP for all i V , since Constraints (4b) can be satisfied when z i ω equals 0. However, the big-M coefficients must not be unnecessarily large, which degrades the performance of SRP-SP. Accordingly, if the columns are fixed to satisfy chance constraints (Constraints (3c)), M i ω could be strengthened more tightly.
We recall P ˜ as a simplified version of RP, where the sizing decisions are not involved. The scenario-based reformulation SP ˜ of P ˜ is stated as follows:
( SRP ˜ ) min i V c 1 , i u i , s . t . ( 1 b ) , ( 1 e ) , ( 1 g ) , q Q F ω q j N i + x j i q Υ u i + M ˜ i ω 1 z i ω , i V , ω Ω , ω Ω z i ω α | Ω | , i V , z i ω { 0 , 1 } , i V , ω Ω .
Indeed, although SRP ˜ does not include the sizing variables v i , the sizing decisions can also be made according to the binary variable z i ω . For ω Ω and i V , if z i ω = 1 , then the sizing decision at node i could be q Q F ω q j N i + x ^ j i q , where x ^ j i q is the optimal value of variable x j i q ; otherwise, the sizing decision at node i could be Υ . Consequently, SRP ˜ is actually equivalent to SRP.
In the following, we have to find a valid M ˜ i ω in SRP ˜ for all ω Ω and simultaneously make their values as small as possible. We are able to achieve such a goal by solving the following optimization problem for a particular ω Ω ,
M ¯ i ω = max q Q F ω q j N i + x j i q Υ , s . t . ( 1 b ) , ( 1 g ) , q Q F ω q j N i + x j i q Υ + M ˜ i ω 1 z i ω , i V , ω Ω , ω Ω z i ω α | Ω | , i V , z i ω { 0 , 1 } , i V , ω Ω .
Then, M ˜ i ω can be valued greater or equal to M ¯ i ω for all ω Ω . In practice, however, it is computationally intensive to determine the exact value of M ¯ i ω . There exist several approaches reported that are efficiently compute moderate upper bounds of M ¯ i ω . We adapt the approach established in [40], which proposes a simple yet efficient algorithm that performs how to generate M ˜ i ω . This approach is a generalization of the quantile approximation approach utilized in [34,41]. We consider the following continuous knapsack problem for particular ω Ω
(5a) max q Q F ω q j N i + x j i q , (5b) s . t . q Q F ω q j N i + x j i q Υ , ω Ω , (5c) x j i q 0 , j N i + , q Q .
Because x j i q are continuous variables, we leverage the greedy-based algorithm to obtain the optimal solution of problem (5), in which we denote L = ( 1 α ) | Ω | as the upper bound of the total number of violated scenarios for i V , where · is the flooring function of real numbers. Then, we present the big-M coefficients generation algorithm as follows:
In Algorithm 1, we follow a decreasing order of F ω q / F ω q for all q Q (O-D pair q should traverse node i) to fill up the value of allocation variables x j i q to their upper bounds 1, which correspond to Steps 3–15. Since the chance constraints can be violated in at most L scenarios, we merely require the ( L + 1 )-th smallest value among all m ˜ i ω ω ’s, rather than max ω Ω m ˜ i ω ω . This corresponds to Step 18 in Algorithm 1.
Algorithm 1 Big-M Coefficients Generation algorithm.
for  i V do
for ω Ω do
  for ω Ω do
    for q Q do
     if q traverses node i then
      Compute ζ i q = F ω q / F ω q .
     end if
    end for
    Sort all ζ i q ’s as a descending order, ζ i q 1 ζ i q 2 ζ i q | Q i | .
    1 , a 0 , b Υ .
    while b > 0 do
     c min 1 , b / F ω q .
     a a + c F ω q , b b c F ω q .
     + 1 .
    end while
    m ˜ i ω ω a .
  end for
   M ˜ i ω the ( L + 1 ) - th smallest number of m ˜ i ω ω for all ω Ω .
end for
end for
It is worth mentioning that [28] considers the fast-charging station location problem using Benders decomposition. However, Benders decomposition is inapplicable for SRP, because the big-M coefficients M i ω and integer variables z i ω introduced in Constraints (4b) and (4c) impede obtaining Benders cuts through the duality method. To tackle this obstacle, we exploit Dantzig–Wolfe decomposition [30], which decomposes SRP into a master problem and several subproblems, and simultaneously, the column generation procedures are employed to obtain the optimal solutions.

3.2. Dantzig–Wolfe Decomposition for SRP

In this section, we apply Dantzig–Wolfe decomposition for solving SRP. Different from the existing works that employ Dantzig–Wolfe decomposition, which usually generate one master problem and several subproblems [29,30,31,32,33,34,42], we separate SRP into Q sub-components corresponding to the Q O-D pairs, then generate Q pairs of master problem and subproblem using Dantzig–Wolfe decomposition to find their optimal solutions, where the sub-components are decomposed by the available arcs between O-D pairs. Lastly, we obtain the optimal solution of SRP by synthesizing the optimal solutions of the sub-components. The sub-components of SRP are as follows
(6a) ( SRP ( q ) ) min i V q c 1 , i u i q + c 2 , i v i q , (6b) s . t . j N i x i j q j N i + x j i q = 1 , if i = s q , 1 , if i = t q , 0 , if i s q , i t q , (6c) F ω q j N i + x j i q v i + M i ω 1 z i ω q , i V , ω Ω , (6d) v i q Υ u i q , i V q , (6e) ω Ω z i ω q α | Ω | , i V q , (6f) z i ω q { 0 , 1 } , i V q , ω Ω , (6g) u i q { 0 , 1 } , i V q , (6h) v i q 0 , i V q . (6i) x i j q 0 , ( i , j ) A ¯ q .
In the master problem for O-D pair q, the allocation variable x i j q ’s are represented by the convex combination of the extreme points of the feasible region in subproblems. Let H q be the index set of extreme points of allocation variables on available arc ( i , j ) between O-D pair q. We then express the allocation variables in the master problem for O-D pair q as
(7a) h H q j N i λ i j h q x ^ i j h q = j N i x i j q , ( i , j ) A ¯ q , (7b) h H q j N i + λ j i h q x ^ j i h q = j N i + x j i q , ( j , i ) A ¯ q , (7c) h H q λ i j h q = 1 , ( i , j ) A ¯ q , (7d) λ i j h q 0 , h H q , ( i , j ) A ¯ q .
where x ^ j i h q represents the extreme points of allocation variables of subproblems. We then apply Dantzig–Wolfe decomposition to SRP using (6), and denote the multiscenario, column-oriented master problem as SRP-MP(q), where H q is the index set of all corner points of O-D pair q.
(8a) ( SRP - MP ( q ) ) min h H q Δ ^ h q ( i , j ) A ¯ q λ i j h q , (8b) s . t . h H q j N i λ i j h q x ^ i j h q h H q j N i + λ j i h q x ^ j i h q = 1 , if i = s q , 1 , if i = t q , 0 , if i s q , i t q , (8c) h H q λ i j h q = 1 , ( i , j ) A ¯ q , (8d) λ i j h q 0 , h H q , ( i , j ) A ¯ q ,
where Δ ^ h q = i V c 1 , i u ^ i h q + c 2 , i v ^ i h q denotes the planning cost at node i for corner point h H q . Objective function (7a) minimize the combination of the planning costs of all paths. Constraints (7b) ensure that the charging demands of each path allocated at each node are balanced. Constraints (7c) and (7d) are convex combination constraints.
The huge number of extreme points of x j i q result in countless variables λ j i h q in SRP-MP(q), which makes it difficult to explicitly investigate all of them. In addition, considering any arbitrary subset of index set H q may cause a suboptimal solution. Thus, it is necessary to employ the column generation method to screen out the indexes that are useful to obtain the optimal solution. Accordingly, we consider a restricted master problem SRP-RMP(q) that is identical to SRP-MP(q) except the index set H q is replaced by a subset of it. Then, we solve the SRP-RMP(q) to extract the optimal solutions and dual variables. Lastly, we generate new columns corresponding to allocation variables with negative reduced cost according to the solutions of subproblem. We generate | Q | subproblems with respect to the reduced cost of λ j i h q for all i V and j N i + , which are referred to as SRP-SP(q). Moreover, the independency among the subproblems for different paths make them ideal to be parallelized. For simplicity, we suppose that SRP-SP(q) is always feasible, that is, H q and q Q .
( SRP - SP ( q ) ) ζ q S R P S P = min Δ q π ^ i s q q j N i s q x i s q j q π ^ i t q q j N i t q + x j i t q q (9a) i V π ^ i q j N i x i j q j N i + x i j q ( i , j ) A ¯ q μ ^ i j q , (9b) s . t . F ω q j N i + x j i q v i q + M i ω 1 z i ω q , i V , ω Ω , (9c) ω Ω z i ω q α | Ω | , i V , (9d) v i q Υ u i q , i V , (9e) j N i x i j q 1 , i V , (9f) j N i + x j i q 1 , i V , (9g) j N i s q x i s q j q = 1 , (9h) j N i t q + x j i t q q = 1 , (9i) z i ω q { 0 , 1 } , i V , ω Ω , (9j) u i q { 0 , 1 } , v i q 0 , i V , (9k) x i j q 0 , ( i , j ) A ¯ q ,
where Δ q = i V c 1 , i u i q + c 2 , i v i q is defined by the aggregated planning cost of fast-charging stations on O-D pair q, and π ^ i q and μ ^ i j q are the value of dual variables of Constraints (7b) and (7c), respectively. Objective (8a) minimizes the total planning cost minus the summation of dual prices associated with the column of O-D pair q, which is the reduced cost of variable λ i j h q in SRP-MP(q). Constraints (9e)–(9h) originate from the flow conservation Constraints (1b), which ensures that the unitary allocation on the arcs be respected.
After subproblem SRP-SP(q) is solved with optimal values u i q * , v i q * , and x i j q * , and it generates an improving column h + 1 to SRP-MP(q) for negative ξ q S R P S P , we then solve SRP-MP(q) with new column that u i , h + 1 q = u i q * , v i , h + 1 q = v i q * , and x ^ i j , h + 1 q = x i j q * . The repeated procedures between solving SRP-MP(q) and SRP-SP(q) operate until no column is favorable, that is, ξ q S R P S P is positive for all q Q . If none of the columns is improving toward SRP-MP(q), it indicates SRP is already solved to optimality with u i * = max q Q i u i q * and v i * = q Q i v i q * . In fact, we have multiple improving columns for SRP-MP(q), such as adding the first improving column, the most improving column, or all of these columns [43]. In practice, however, u i * and v i * may not satisfy Constraints (3c), because the sizing bound Υ can not be separated with respect to the paths. A feasible way to overcome this obstacle is to choose a proper Υ , such that the fast-charging station is able to accommodate the allocated charging demands from all possible paths and simultaneously retains an acceptable solution time of SRP-SP(q).
Remark 1. 
When conducting the numerical experiment, we are aware that  SRP  is computationally demanding in a large-scale network using off-the-shelf solves. Even when we use the conventional Dantzig–Wolfe decomposition, it still consumes massive amount of CPU time to find the optimal solution of the master problem because  SRP  incorporates a multitude of decision variables λ i j h q all over the expanded network. Compared with the existing works, our proposed approach enables us to improve the computational efficiency twofold: one is that  SRP  is separated into Q disjoint sub-components, which are facilitated to be solved using parallel machines; the other is that the sub-components merely incorporate the decision variables corresponding to available arcs between specific O-D pairs, which contain far fewer decision variables in the Dantzig–Wolfe master problem, and thus simplifies the computation to a extremely large extent.

4. Numerical Experiments

In this section, a network topology has been used to evaluate the performance of the proposed model and Dantzig–Wolfe decomposition method. The solving algorithms are coded in Matlab R2017b called YALMIP (version: R20210331) [44] using CPLEX 12.10 on a Windows machine with Intel(R) Core(TM) i9-11950H and 64 GB RAM. The building cost of the fast-charging stations and the installing cost of the charging piles are set to be 100 and 10, respectively.
We test the proposed model in a hypothetical 25-node network shown in Figure 5, which is originally proposed by [45] and has been extensively used as a benchmark network for the charging station location problem [26,27,36]. The circled numbers denote the labels of the nodes, the line between two nodes denotes the link that connects them, and the number in each link denotes the length of the link. This network possesses 25 nodes and 86 arcs (43 undirected links). Every node is considered as an origin, destination, or candidate site for the fast-charging stations, yielding in total 300 O-D pairs and 25 candidate sites. We suppose that the EVs start their trips with fully charged batteries, since they can be charged through slow chargers during their leisure time at home or work place [46]. For the scenario-based reformulation, we investigate 10 scenarios of charging demands between O-D pairs. In each scenario ω , the charging demands between each O-D pair are equal, F ω q = 5 + w for all q Q .
To evaluate the performance of Dantzig–Wolfe decomposition reformulations, we generate 120 instances by considering 4 types of driving ranges, 5 levels of service, abilities and 6 values of sizing limits:
  • Driving ranges: R = 8 , 9 , 10 , 11 .
  • Service levels: α = 10 % , 20 % , 30 % , 40 % , 50 % .
  • Sizing limits: Υ = 800 , 900 , 1000 , 1100 , 1200 , 1300 .
Table A1, Table A2, Table A3, Table A4 and Table A5 in Appendix A illustrate the results of numerical experiments in the 25-node network for the five different service levels α = 10 % , 20 % , 30 % , 40 % , 50 % , respectively. In each table, the rows correspond to the results of a particular instance indicated by the service level ( α ), the driving ranges (R), and the sizing limit ( Υ ). These tables compare the numerical performances between the original model SRP(q) and Dantzig–Wolfe decomposition SRP-MP(q) and SRP-SP(q), in which the total planning costs of fast-charging stations (Obj.) and the runtime (in seconds) of such two models as well as the iteration rounds (in Num.) of Dantzig–Wolfe decomposition are provided. Especially, the maximal (Max.), median (Med.), average (Ave.), and minimal (Min.) runtime and iteration rounds are reported in detail. In general, the CPLEX solver terminates with a relatively optimal tolerance in 10 4 . The upper bound of runtime is set to be 14,400 s (4 h); the results marked by superscript 14 , 400   * indicate that optimal solutions cannot be found within 14,400 s, while the results marked by—means that we cannot obtain the exact output within 14,400 s. Finally, the initial values of the allocation variable are selected randomly with respect to Constraints (1b) and (1g).
According to the results in Table A1, Table A2, Table A3, Table A4 and Table A5 in Appendix A, we observe that Dantzig–Wolfe decomposition SRP-MP(q) and SRP-SP(q) can significantly reduce the runtime of solving the original model SRP(q) directly. For instance, in Table A3, the row under service level α = 10 % , driving range R = 8 , and sizing bound Υ = 800 manifests that, after leveraging Dantzig–Wolfe decomposition, the maximal runtime is reduced from 7364.90 s to 0.22 s, the median runtime is reduced from 14.27 s to 0.16 s, and the average runtime is reduced from 100.28 s to 0.16 s. The similar runtime reduction effect can also be observed through the numerical results shown in other tables, which demonstrates that Dantzig–Wolfe decomposition gives rise to higher numerical efficiency than solving the original model directly. In addition, for some instances, such as service level α = 10 % and driving range R = 11 , we cannot obtain the optimal planning strategies directly, whereas Dantzig–Wolfe decomposition can still effectively obtain the optimal planning results. Table A1, Table A2, Table A3, Table A4 and Table A5 also show the iteration rounds of SRP-MP(q) and SRP-SP(q), which depend on the selected initial values of the allocation variables. We observe that Dantzig–Wolfe decomposition can find the optimal planning strategy within 10 rounds of iteration and thus is able to guarantee its numerical efficiency.
In addition to the computational efficiency of Dantzig–Wolfe decomposition, Table A1, Table A2, Table A3, Table A4 and Table A5 also report the effects of service level α and the driving range R on the optimal planning costs. Although the optimal planning costs of the original model and Dantzig–Wolfe decomposition are nonidentical, they share a common changing tendency with respect to the change of α and R. On the whole, the optimal planning costs grow with the increase of service level α . For instance, under driving range R = 8 , the optimal planning costs grow with respect to the increasing α . The reason is that larger α usually implies more charging piles should be provided to the EVs, which stimulates the increase of planning costs. Moreover, the optimal planning costs decline with the increase of driving range R. For instance, under driving range α = 20 % , the optimal planning costs decline with respect to the increasing R. The reason is that larger R means EVs are able to drive a longer distance without refueling, which is conducive to reducing the charging demands and thus decreasing the planning costs. On the other hand, the sizing bound Υ seems to not affect the planning costs of fast-charging stations according to Table A1, Table A2, Table A3, Table A4 and Table A5. In practice, however, Υ determines the feasibility of the planning models, that is, some fast-charging stations may be unable to accommodate the allocated charging demands due to the limitation of installing charging piles. Thus, it is important to set a proper Υ before establishing the planning model. Finally, the number of scenarios | Ω | profoundly affects the runtime of our models. It is obvious that, if Ω contains more scenarios, the models will encounter more constraints, and solving the model thus requires more runtime. Nevertheless, the models with more scenarios can reflect the real world more closely. Hence, a moderate size of scenarios that makes a good trade-off between the computational efficiency and the practical significance plays an indispensable role in establishing the planning models.

5. Conclusions

We addressed the modeling and solution approaches of the fast-charging station planning problem. We formulated a CCP model that exploits the multicommodity flow model to handle the uncertain charging demands with unknown probability distributions. We first applied a scenario-based reformulation approach that reformulates the CCP model into a mixed-integer linear program, and further enhanced the performance of the scenario-based formulation using a coefficients-strengthening algorithm, which aimed to find the appropriate big-M coefficients. We then proposed Dantzig–Wolfe decomposition to find the optimal solutions of the planning model. Finally, we showed that the proposed Dantzig–Wolfe decomposition method considerably improved the computational efficiency for solving a large-scale problem instance in reality. It should be mentioned that Dantzig–Wolfe decomposition can also be applied to the real-world transportation network, such as the Texas highway network and the California Road network. Such applications will be considered as a future work.

Author Contributions

Conceptualization, B.Z.; methodology, L.W.; software, L.W.; validation, B.Z.; formal analysis, L.W.; investigation, B.Z.; resources, B.Z.; data curation, B.Z.; writing—original draft preparation, L.W.; writing—review and editing, B.Z.; visualization, L.W.; supervision, B.Z.; project administration, B.Z.; funding acquisition, B.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported in part by the National Natural Science Foundation of China [61803056].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data will be made available on request.

Conflicts of Interest

There is no conflict of interest for this study.

Abbreviations

The following abbreviations are used in this manuscript:
SetsDescriptions
G V , A Original network G , where V and A are sets of nodes and arcs.
G ¯ q V ¯ q , A ¯ q Expanded network on O-D pair q, where V ¯ q and A ¯ q are sets of nodes and arcs.
QSet of O-D pairs in original network, indexed by q.
N i + , N i Set of in-neighbors and out-neighbors of node i on expanded network.
Ω Set of scenarios, indexed by ω .
V q Set of nodes that traversed by the path between O-D pair q.
H q Index set of extreme points of allocation variables on available arc ( i , j )
between O-D pair q.
ParametersDescriptions
s q Added source node of O-D pair q on expanded network.
t q Added sink node of O-D pair q on expanded network.
F ˜ q Uncertain charging demands on O-D pair q.
c 1 , i Cost of building a fast-charging station at node i V .
c 2 , i Cost of intalling a charging pile at node i V .
d i j q Driving distance of arc ( i , j ) A ¯ q .
Υ Sizing bound in fast-charging stations.
α Service level of the fast-charging stations.
RDriving ranges of electric vehicles.
F ^ ω q ω Ω Realizations of F ˜ q in scenario ω .
M i ω Big-M coefficient at node i in scenario ω .
π ^ i q Value of dual variables.
μ ^ q Value of dual variables.
x ^ i j h Extreme points of allocation variables.
Δ ^ h q Extreme points of planning cost on O-D pair q.
| Ω | The cardinality of set Ω .
FunctionsDescriptions
· Ceiling function, the lowest integer not less than a real number.
VariablesDescriptions
x i j q Continuous variables, proportion of charging demands allocated on arc ( i , j ) A ¯ q .
y i Integer variables, the sizing of fast-charging stations at node i V .
v i Continuous variables, relaxation of y i , the size of fast-charging stations at node i V .
u i Binary variables, whether fast-charging station sites at node i V .
u i q Binary variable, whether fast-charging station sites at node i to accommodate
the charging demands from O-D pair q.
v i q Continuous variable, the size of fast-charging stations at node i to accommodate
the charging demands from O-D pair q.
λ i j h q Continuous variable, combination weight of the extreme points for all h H q ,
( i , j ) A ¯ q , q Q .
z i ω Binary variable, whether chance constraints is satisfied for all i V , ω Ω .

Appendix A. Experimental Results of the 25-Node Network

Table A1. Results of scenario-based reformulation in the 25-node network with α = 10 % , where superscript * indicates that optimal solutions cannot be found within 14,400 s.
Table A1. Results of scenario-based reformulation in the 25-node network with α = 10 % , where superscript * indicates that optimal solutions cannot be found within 14,400 s.
α = 10 % Performance of SRP(q)Performance of Dantzig–Wolfe Decomposition
R Υ Obj.Runtime (in Sec.)Obj.Runtime of MP (in Sec.)Runtime of SP (in Sec.)Iter. (in Num.)
Max.Med.Ave.Min.Max.Med.Ave.Min.Max.Med.Ave.Min.Max.Med.Ave.Min.
880054,1807364.9014.27100.280.1638,7600.220.160.160.130.280.210.210.167441
90054,1807383.0322.6697.670.2138,7600.290.130.130.080.650.240.240.149561
100054,1807901.7714.91106.170.1538,7600.300.220.220.170.520.340.340.229431
110054,1807880.4217.23117.960.1438,7600.320.160.190.130.450.220.250.167321
120054,1807533.7016.5098.950.1738,7600.310.140.160.120.400.250.250.155341
130054,1807750.6923.94123.610.2438,7600.360.170.180.140.420.210.240.138661
980046,42010,551.686.72150.300.1436,9200.320.190.200.160.390.250.250.207321
90046,42010,266.4929.46174.910.1136,9200.340.180.140.110.340.240.240.206341
100046,42010,365.1910.93150.000.1536,9200.320.190.190.130.310.250.270.215341
110046,42010,206.2734.75186.170.1736,9200.330.160.190.110.360.260.270.229441
120046,42010,451.6041.91164.010.2336,9200.310.180.160.120.360.260.260.217561
130046,42010,553.1816.10198.160.1836,9200.350.120.120.100.390.220.230.199651
1080037,38013,421.004.45214.820.1635,2200.360.240.250.210.460.300.300.247231
90037,38013,137.684.54222.840.0635,2200.350.260.230.200.420.310.320.235221
100037,38013,108.0813.84273.300.1035,2200.340.290.290.220.470.330.340.279231
110037,38013,633.322.03220.920.1435,2200.370.250.250.210.440.340.360.269341
120037,38013,863.3625.54275.850.0835,2200.380.260.260.230.470.330.350.258231
130037,38013,654.7316.97220.270.1735,2200.320.280.250.200.460.320.340.275331
1180014,400 *0.1133,5400.460.280.290.260.490.340.340.289211
90014,400 *0.2333,5400.420.260.250.220.430.370.370.277551
100014,400 *0.1233,5400.450.290.280.240.470.380.370.265341
110014,400 *0.2433,5400.430.270.270.230.470.360.350.299671
120014,400 *0.2033,5400.420.260.260.220.490.340.340.267541
130014,400 *0.1333,5400.470.280.260.210.480.330.320.245211
Table A2. Results of scenario-based reformulation in the 25-node network with α = 20 % .
Table A2. Results of scenario-based reformulation in the 25-node network with α = 20 % .
α = 20 % Performance of SRP(q)Performance of Dantzig–Wolfe Decomposition
R Υ Obj.Runtime (in Sec.)Obj.Runtime of MP (in Sec.)Runtime of SP (in Sec.)Iter. (in Num.)
Max.Med.Ave.Min.Max.Med.Ave.Min.Max.Med.Ave.Min.Max.Med.Ave.Min.
880062,8605308.9010.3432.170.1544,7200.470.350.350.300.740.410.410.345221
90062,1605283.0315.6036.770.1244,7200.440.370.370.330.710.420.420.396341
100062,1605101.7710.6234.790.1744,7200.410.350.340.300.750.430.430.348211
110062,1605280.4210.5932.390.1244,7200.420.380.370.320.760.480.460.358451
120062,1605233.7011.5932.390.1244,7200.490.320.340.310.780.450.450.325231
130062,1605250.699.7936.090.1644,7200.470.320.320.280.700.470.470.399781
980052,9505520.189.7941.940.1342,6900.540.390.390.350.790.450.450.375331
90052,9505533.7513.3842.050.1542,6900.570.330.330.300.740.440.440.379661
100052,9505512.587.8244.120.1442,6900.540.360.360.320.780.410.410.365231
110052,9505490.6812.2938.280.1742,6900.520.390.380.310.720.400.400.386331
120052,9505512.499.8544.090.0942,6900.510.380.360.330.710.440.430.399781
130052,9505510.958.7746.170.0942,6900.560.350.350.300.700.450.440.346231
1080041,8605639.599.7456.850.1040,6900.600.420.420.380.740.470.480.426231
90041,8605519.5910.2654.110.1440,6900.660.400.400.340.790.480.490.425321
100041,8605561.047.4457.170.1540,6900.680.420.420.320.750.460.470.427341
110041,8605656.068.5957.120.1340,6900.630.480.480.330.700.470.470.416451
120041,8605513.6410.4153.250.1140,6900.640.470.470.320.720.410.410.368231
130041,8605659.9210.0952.610.1040,6900.670.410.400.360.710.430.440.407341
1180036,6805759.5010.3866.790.1438,7300.680.470.480.430.810.520.530.467331
90036,6805644.646.4665.520.1238,7300.690.480.490.430.890.570.570.447231
100036,6805618.2810.6459.030.0938,7300.630.470.470.430.840.540.540.465211
110036,6805638.4816.0966.790.0938,7300.660.490.490.420.880.500.510.445331
120036,6805500.0212.5861.660.1538,7300.670.480.480.420.840.570.560.427541
130036,6805632.6113.3367.850.1138,7300.650.480.480.430.850.560.560.458561
Table A3. Results of scenario-based reformulation in the 25-node network with α = 30 % .
Table A3. Results of scenario-based reformulation in the 25-node network with α = 30 % .
α = 30 % Performance of SRP(q)Performance of Dantzig–Wolfe Decomposition
R Υ Obj.Runtime (in Sec.)Obj.Runtime of MP (in Sec.)Runtime of SP (in Sec.)Iter. (in Num.)
Max.Med.Ave.Min.Max.Med.Ave.Min.Max.Med.Ave.Min.Max.Med.Ave.Min.
880071,640137.778.8515.120.0850,8800.760.510.520.460.750.560.570.509341
90071,640129.827.1513.840.0850,8800.700.590.590.490.780.520.520.485231
100071,640139.598.9312.730.1350,8800.710.510.580.450.770.540.540.519761
110071,640103.665.2214.230.1250,8800.790.520.500.420.790.570.570.538321
120071,640119.826.5315.350.1350,8800.750.580.570.410.700.520.520.518231
130071,640127.835.9012.550.0950,8800.770.500.500.440.750.540.540.507441
980060,660193.9910.1919.380.1248,5601.340.570.570.500.890.630.630.559541
90060,660195.6113.5516.260.0948,5601.380.570.570.510.880.650.650.565231
100060,660200.7012.2917.890.1348,5601.360.580.580.520.870.640.640.549541
110060,660196.6914.8816.270.1748,5601.300.540.540.500.880.620.620.568541
120060,660203.7912.2118.310.1648,5601.310.550.550.480.820.640.650.596211
130060,660190.0710.2516.690.1248,5601.310.570.570.520.890.650.650.527231
1080048,950279.1914.9820.680.1346,1600.550.060.060.040.440.120.130.075221
90048,950276.7914.1822.730.1346,1600.590.110.120.070.460.140.140.108441
100048,950268.5613.8621.070.1346,1600.520.090.090.060.400.180.170.128331
110048,950263.1313.2221.560.1146,1600.500.120.110.050.460.120.120.076211
120048,950262.9019.4222.640.1446,1600.510.080.070.030.470.130.140.067221
130048,950267.4214.8321.020.1046,1600.580.070.070.040.420.170.160.096321
1180041,620336.8519.6725.840.1643,9200.160.100.100.080.300.160.170.106231
90041,620323.5020.0726.280.1043,9200.170.120.120.060.310.150.140.056431
100041,620313.9724.1525.110.1043,9200.190.130.120.050.300.160.160.099551
110041,620314.8118.3126.150.1043,9200.230.100.120.040.320.170.170.078441
120041,620331.3818.6326.550.1343,9200.190.110.120.040.390.130.140.047211
130041,620324.5024.6924.180.1343,9200.170.100.120.080.350.120.110.078211
Table A4. Results of scenario-based reformulation in the 25-node network with α = 40 % .
Table A4. Results of scenario-based reformulation in the 25-node network with α = 40 % .
α = 40 % Performance of SRP(q)Performance of Dantzig–Wolfe Decomposition
R Υ Obj.Runtime (in Sec.)Obj.Runtime of MP (in Sec.)Runtime of SP (in Sec.)Iter. (in Num.)
Max.Med.Ave.Min.Max.Med.Ave.Min.Max.Med.Ave.Min.Max.Med.Ave.Min.
880080,32042.1412.9015.400.0956,9400.210.140.150.120.320.200.210.146341
90079,52042.4112.8317.820.1356,9400.250.110.110.090.360.220.220.105211
100079,52039.1410.7413.510.1556,9400.230.140.140.110.310.240.240.117541
110079,52041.009.5413.090.1356,9400.280.150.160.120.390.210.210.186231
120079,52040.699.4316.070.1656,9400.220.100.100.070.380.280.270.175321
130079,52043.4011.4214.330.1756,9400.210.120.120.080.370.250.250.188651
980067,08048.8110.7214.360.1154,2300.290.190.190.160.330.250.250.187221
90066,98047.6810.7417.490.1854,2300.240.140.140.110.390.250.240.176431
100066,98048.3912.3813.190.0954,2300.270.180.180.140.300.220.210.136231
110066,98049.419.8814.660.1554,2300.280.170.170.120.350.210.220.148651
120066,98047.399.5717.780.1354,2300.230.190.180.130.430.260.250.158331
130066,98048.069.6917.960.1554,2300.280.150.150.100.310.250.250.116451
1080053,72052.8410.4513.770.1051,6300.730.230.230.200.620.290.290.226221
90053,72052.9011.9213.480.1151,6300.740.260.260.210.620.260.260.206341
100053,62053.3912.4115.470.1351,6300.700.240.240.200.610.280.280.226451
110053,62054.219.6316.780.1651,6300.730.270.260.220.620.250.250.206331
120053,62055.8710.3017.660.1151,6300.750.260.260.220.610.250.250.219341
130053,62055.2610.6517.840.1651,6300.700.220.220.170.690.240.240.218341
1180046,76057.139.0710.920.1449,1100.370.270.270.240.470.330.340.276231
90046,76056.0410.489.440.0949,1100.340.240.240.200.400.300.310.228321
100046,76058.5010.2415.010.0849,1100.350.210.210.170.450.380.380.319771
110046,76055.5210.6910.100.1649,1100.350.220.220.160.420.320.320.259661
120046,76058.3910.1317.790.1549,1100.390.230.230.140.480.350.350.279321
130046,76059.849.7213.270.1049,1100.320.220.220.180.410.320.320.215211
Table A5. Results of scenario-based reformulation in the 25-node network with α = 50 % .
Table A5. Results of scenario-based reformulation in the 25-node network with α = 50 % .
α = 50 % Performance of SRP(q)Performance of Dantzig–Wolfe Decomposition
R Υ Obj.Runtime (in Sec.)Obj.Runtime of MP (in Sec.)Runtime of SP (in Sec.)Iter. (in Num.)
Max.Med.Ave.Min.Max.Med.Ave.Min.Max.Med.Ave.Min.Max.Med.Ave.Min.
880089,00034.4812.1314.280.1562,8000.680.470.470.420.760.520.510.447441
90088,60033.6312.6314.750.1762,8000.630.490.480.460.720.560.560.465331
100088,20032.779.9314.980.1362,8000.600.460.460.410.730.540.540.427211
110088,20032.8912.2213.130.1362,8000.670.450.450.400.780.590.590.439771
120088,20032.9310.9113.730.1262,8000.630.490.490.450.740.500.500.456331
130088,20032.0913.1412.850.0862,8000.660.410.420.380.730.560.560.446231
980075,43035.6513.3815.160.1160,0000.510.070.070.040.590.130.130.075211
90075,43034.8312.0514.400.160,0000.550.080.080.040.590.140.140.075221
100075,43032.9811.7613.460.1260,0000.530.090.090.030.510.100.110.058441
110075,43032.3210.499.000.1760,0000.550.080.080.050.550.150.140.075331
120075,43032.207.5211.270.1660,0000.530.070.070.030.520.130.130.055211
130075,43032.997.7312.170.1760,0000.500.090.090.060.500.190.190.086211
1080060,90033.7712.0014.240.1557,1000.590.110.110.080.880.170.170.117341
90060,90032.8310.6714.480.1157,1000.580.140.140.080.810.180.170.147341
100060,90032.5110.8212.810.1857,1000.560.120.120.040.860.150.150.077331
110060,90032.488.1413.320.1057,1000.590.120.120.070.810.240.240.138561
120060,90033.6910.9315.370.1757,1000.560.100.100.060.860.220.220.165221
130060,90033.767.4114.150.1757,1000.530.110.110.090.880.190.190.136211
1180051,50037.427.367.790.0954,3000.530.060.070.045.260.140.170.096341
90051,50037.877.1913.460.1454,3000.530.150.150.114.390.240.240.145341
100051,50036.119.2714.330.1054,3000.590.190.190.154.590.210.210.146231
110051,50038.2811.2814.570.0854,3000.540.130.130.104.450.240.240.117551
120051,50035.5413.1615.160.1254,3000.530.120.130.094.370.230.240.196321
130051,50038.6811.9112.360.1854,3000.260.180.180.154.410.250.240.195221

References

  1. Alshareef, S.M. Analyzing and mitigating the impacts of integrating fast-charging stations on the power quality in electric power distribution systems. Sustainability 2022, 14, 5595. [Google Scholar] [CrossRef]
  2. Yang, Y.; Tan, Z.F.; Ren, Y.L. Research on factors that influence the fast charging behavior of private battery electric vehicles. Sustainability 2020, 12, 3439. [Google Scholar] [CrossRef] [Green Version]
  3. Chang, M.; Bae, S.; Cha, G.; Yoo, J. Aggregated electric vehicle fast-charging power demand analysis and forecast based on LSTM neural network. Sustainability 2021, 13, 13783. [Google Scholar] [CrossRef]
  4. Luedtke, J.; Ahmed, S. A sample approximation approach for optimization with probabilistic constraints. SIAM J. Optim. 2008, 19, 674–699. [Google Scholar] [CrossRef] [Green Version]
  5. Upchurch, C.; Kuby, M. Comparing the p-median and flow-refueling models for locating alternative-fuel stations. J. Transp. Geogr. 2010, 18, 750–758. [Google Scholar] [CrossRef]
  6. Abareshi, M.; Zaferanieh, M. A bi-level capacitated p-median facility location problem with the most likely allocation solution. Transp. Res. Part B Methodol. 2019, 123, 1–20. [Google Scholar] [CrossRef]
  7. An, Y.; Zeng, B.; Zhang, Y.; Zhao, L. Reliable p-median facility location problem: Two stage robust models and algorithms. Transp. Res. Part B Methodol. 2014, 64, 54–72. [Google Scholar] [CrossRef]
  8. Jackson, L.E.; Rouskas, G.N.; Stallmann, M.F.M. The directional p-median problem: Definition, complexity, and algorithms. Eur. J. Oper. Res. 2007, 179, 1097–1108. [Google Scholar] [CrossRef] [Green Version]
  9. Scutaru, M.L.; Vlase, S.; Marin, M.; Modrea, A. New analytical method based on dynamic response of planar mechanical elastic systems. Bound. Value Probl. 2020, 1, 104. [Google Scholar] [CrossRef]
  10. Alzahrani, F.; Hobiny, A.; Abbas, I.; Marin, M. An eigenvalues approach for a two-dimensional porous medium based upon weak, normal and strong thermal conductivities. Symmetry 2020, 12, 848. [Google Scholar] [CrossRef]
  11. Hodgson, M.J. A flow-capturing location-allocation model. Geogr. Anal. 1990, 22, 270–279. [Google Scholar] [CrossRef]
  12. Kuby, M.; Lim, S. The flow-refueling location problem for alternative-fuel vehicles. Socio Econ. Plan. Sci. 2005, 39, 125–145. [Google Scholar] [CrossRef]
  13. Shen, Z.J.M.; Feng, B.; Mao, C.; Ran, L. Optimization models for electric vehicle service operations: A literature review. Transp. Res. Part B Methodol. 2019, 128, 462–477. [Google Scholar] [CrossRef]
  14. Mak, H.Y.; Rong, Y.; Shen, Z.J.M. Infrastructure planning for electric vehicles with battery swapping. Manag. Sci. 2013, 59, 423–443. [Google Scholar] [CrossRef] [Green Version]
  15. Zhang, H.C.; Moura, S.J.; Hu, Z.C.; Qi, W.; Song, Y.H. A second-order cone programming model for planning PEV fast-charging stations. IEEE Trans. Power Syst. 2019, 33, 2763–2777. [Google Scholar] [CrossRef] [Green Version]
  16. Zhang, H.C.; Moura, S.J.; Hu, Z.C.; Qi, W.; Song, Y.H. Joint PEV charging network and distributed PV generation planning based on accelerated generalized Benders decomposition. IEEE Trans. Transp. Electrif. 2018, 4, 789–803. [Google Scholar] [CrossRef]
  17. Ghelichi, Z.; Gentili, M.; Mirchandani, P.B. Drone logistics for uncertain demand of disaster-impacted populations. Transp. Res. Part C Emerg. Technol. 2022, 141, 103735. [Google Scholar] [CrossRef]
  18. Yang, D.T.; Sarma, N.J.S.; Hyland, M.F.; Jayakrishnan, R. Dynamic modeling and real-time management of a system of EV fast-charging stations. Transp. Res. Part C Emerg. Technol. 2021, 128, 103186. [Google Scholar] [CrossRef]
  19. Ye, B.J.; Ni, C.; Tian, Y.; Ochieng, W.Y. Data-driven distributionally robust generation of time-varying flow corridor networks under demand uncertainty. Transp. Res. Part C Emerg. Technol. 2022, 136, 103546. [Google Scholar] [CrossRef]
  20. Wu, F.; Soishansi, R. A stochastic flow-capturing model to optimize the location of fast-charging stations with uncertain electric vehicle flows. Transp. Res. Part D Transp. Environ. 2017, 53, 354–376. [Google Scholar] [CrossRef]
  21. Zhou, B.; Chen, G.; Huang, T.W.; Song, Q.K.; Yuan, Y.F. Planning PEV fast-charging stations using data-driven distributionally robust optimization approach based on Φ-divergence. IEEE Trans. Transp. Electrif. 2020, 6, 170–180. [Google Scholar] [CrossRef]
  22. Zhou, B.; Chen, G.; Song, Q.K.; Dong, Z.Y. Robust chance-constrained programming approach for the planning of fast-charging stations in electrified transportation networks. Appl. Energy 2020, 262, 114480. [Google Scholar] [CrossRef]
  23. Li, Y.; Kim, A.M. Allocating and scheduling resources for a mobile photo enforcement program. Transp. Res. Part C Emerg. Technol. 2019, 125, 103000. [Google Scholar] [CrossRef]
  24. Zhou, S.; Jamshidi, A.; Núñez, A.; Baldi, S.; Schutter, B.D. Integrated condition-based track maintenance planning and crew scheduling of railway networks. Transp. Res. Part C Emerg. Technol. 2019, 105, 359–384. [Google Scholar]
  25. Liu, J.T.; Mirchandani, P.; Zhou, X.S. Integrated vehicle assignment and routing for system-optimal shared mobility planning with endogenous road congestion. Transp. Res. Part C Emerg. Technol. 2020, 117, 102675. [Google Scholar] [CrossRef]
  26. Yıldız, B.; Arslan, O.; Karaşan, Q.E. A branch and price apporach for routing and refueling station location model. Eur. J. Oper. Res. 2016, 248, 815–826. [Google Scholar] [CrossRef] [Green Version]
  27. Xu, M.; Meng, Q. Optimal deployment of charging stations considering path deviation and nonlinear elastic demand. Transp. Res. Part B Methodol. 2020, 135, 120–142. [Google Scholar] [CrossRef]
  28. Lee, C.; Han, J. Benders-and-price approach for electric vehicle charging station location problem under probabilistic travel range. Transp. Res. Part B Methodol. 2017, 106, 130–152. [Google Scholar] [CrossRef]
  29. Appelgren, L.H. A column generation algorithm for a ship scheduling problem. Transp. Sci. 1969, 3, 53–68. [Google Scholar] [CrossRef]
  30. Dantzig, G.B.; Wolfe, P. Decomposition principle for linear programs. Oper. Res. 1960, 8, 101–111. [Google Scholar] [CrossRef]
  31. Singh, K.J.; Philpott, A.B.; Wood, R.K. Dantzig-Wolfe decomposition for solving multistage stochastic capacity-planning problems. Oper. Res. 2009, 57, 1271–1286. [Google Scholar] [CrossRef] [Green Version]
  32. Kulkarnia, S.; Krishnamoorthyd, M.; Ranadef, A.; Ernstc, A.T.; Patilb, R. A new formulation and a column generation-based heuristic for the multiple depot vehicle scheduling problem. Transp. Res. Part B Methodol. 2018, 118, 457–487. [Google Scholar] [CrossRef]
  33. Wu, T.; Shi, Z.S.; Liang, Z.; Zhang, Z.N.; Zhang, C.R. Dantzig-Wolfe decomposition for the facility location and production planning problem. Comput. Oper. Res. 2020, 124, 105068. [Google Scholar] [CrossRef]
  34. Zhang, Z.; Dention, B.T.; Xie, X.L. Branch and price for chance-constrained bin packing. INFORM J. Comput. 2020, 32, 547–564. [Google Scholar] [CrossRef]
  35. Cruz, C.A.; Munari, P.; Morabito, R. A branch-and-price method for the vehicle allocation problem. Comput. Ind. Eng. 2020, 149, 106745. [Google Scholar] [CrossRef]
  36. MirHassani, S.A.; Ebrazi, R. A flexible reformulation of the refueling station location problem. Transp. Sci. 2013, 47, 617–628. [Google Scholar] [CrossRef]
  37. Wang, C.Z.; He, F.; Lin, X.; Shen, Z.-J.M.; Li, M. Designing locations and capacities for charging stations to support intercity travel of electric vehicles: An expanded network approach. Transp. Res. Part C Emerg. Technol. 2019, 102, 210–232. [Google Scholar] [CrossRef]
  38. Zhang, H.C.; Moura, S.J.; Hu, Z.C.; Song, Y.H. PEV fast-charging station siting and sizing on coupled transportation and power networks. IEEE Trans. Smart Grid 2018, 9, 2595–2605. [Google Scholar] [CrossRef]
  39. Xu, M.; Yang, H.; Wang, S.A. Mitigate the range anxiety: Siting battery charging stations for electric vehicle drivers. Transp. Res. Part C Emerg. Technol. 2020, 114, 164–188. [Google Scholar] [CrossRef]
  40. Deng, Y.; Shen, S.Q.; Denton, B. Chance-constrained surgery planning under conditions of limited and ambiguous data. INFORMS J. Comput. 2019, 31, 559–575. [Google Scholar] [CrossRef]
  41. Song, Y.; Luedtke, J.R.; Küçükyavuz, S. Chance-constrained binary packing problems. INFORMS J. Comput. 2014, 26, 735–747. [Google Scholar] [CrossRef] [Green Version]
  42. Liu, J.T.; Zhou, X.S. Observability quantification of public transportation systems with heterogeneous data sources: An information-space projection approach based on discretized space-time network flow models. Transp. Res. Part B Methodol. 2019, 128, 302–323. [Google Scholar] [CrossRef]
  43. Lamiri, M.; Xie, X.; Zhang, S. Column generation approach to operating theater planning with elective and emergency patients. IIE Trans. 2008, 40, 838–852. [Google Scholar] [CrossRef]
  44. Lofberg, J. YALMIP: A toolbox for modeling and optimization in MATLAB. In Proceedings of the IEEE International Symposium on Computer Aided Control Systems Design, Taipei, Taiwan, 2–4 September 2004; pp. 284–289. [Google Scholar]
  45. Simchi-Levi, D.; Berman, O. A heuristic algorithm for the traveling salesman location problem on networks. Oper. Res. 1988, 36, 478–484. [Google Scholar] [CrossRef]
  46. Trigg, T.; Telleen, P.; Boyd, R.; Cuenot, F.; D’Ambrosio, D.; Gaghen, R.; Gagné, J.; Hardcastle, A.; Houssin, D.; Jones, A. Global EV outlook: Understanding the electric vehicle landscape to 2020. Int. Energy Agency 2013, 1–40. Available online: https://iea.blob.core.windows.net/assets/af46e012-18c2-44d6-becd-bad21fa844fd/Global_EV_Outlook_2020.pdf (accessed on 1 March 2023).
Figure 1. The original transportation network G ( V , A , Q ) , where the set of nodes is V = { O , A , B , C , E , D } and the set of arcs is A = { ( O , A ) , ( A , B ) , ( B , C ) , ( C , E ) , ( E , D ) } . The distances between nodes are labeled above the arcs.
Figure 1. The original transportation network G ( V , A , Q ) , where the set of nodes is V = { O , A , B , C , E , D } and the set of arcs is A = { ( O , A ) , ( A , B ) , ( B , C ) , ( C , E ) , ( E , D ) } . The distances between nodes are labeled above the arcs.
Sustainability 15 06588 g001
Figure 2. Step 1: add a source node s before original node O and a sink node t after destination node D. The expanded node set is V ¯ = { s , O , A , B , C , E , D , t } . Step 2: connect s and O, D and t, and add ( s , O ) and ( D , t ) to A. Then, A ¯ = { ( s , O ) , ( O , A ) , ( A , B ) , ( B , C ) , ( C , E ) , ( E , D ) , ( D , t ) } .
Figure 2. Step 1: add a source node s before original node O and a sink node t after destination node D. The expanded node set is V ¯ = { s , O , A , B , C , E , D , t } . Step 2: connect s and O, D and t, and add ( s , O ) and ( D , t ) to A. Then, A ¯ = { ( s , O ) , ( O , A ) , ( A , B ) , ( B , C ) , ( C , E ) , ( E , D ) , ( D , t ) } .
Sustainability 15 06588 g002
Figure 3. Step 3: order the nodes from s to t sequentially. The number in the bracket denotes the ordering index of node. Step 4: add ( s , A ) and ( s , B ) to A ¯ since the distances between s, A and s, B are strictly less than 60 km; add ( E , t ) and ( D , t ) to A ¯ since the distances between E, t and D, t are strictly less than 60 km; add ( O , B ) , ( O , C ) , ( A , C ) , ( B , E ) , ( C , D ) to A ¯ since the distances between these node pairs are strictly less than 120 km. Then, we formulate the expanded arc set A ¯ = { ( s , O ) , ( s , A ) , ( s , B ) , ( O , A ) , ( O , B ) , ( O , C ) , ( A , B ) , ( A , C ) , ( B , C ) , ( B , E ) , ( C , E ) , ( C , D ) , ( E , D ) , ( E , t ) , ( D , t ) } . Combining V ¯ and A ¯ , we obtain the expanded network G ¯ V ¯ , A ¯ .
Figure 3. Step 3: order the nodes from s to t sequentially. The number in the bracket denotes the ordering index of node. Step 4: add ( s , A ) and ( s , B ) to A ¯ since the distances between s, A and s, B are strictly less than 60 km; add ( E , t ) and ( D , t ) to A ¯ since the distances between E, t and D, t are strictly less than 60 km; add ( O , B ) , ( O , C ) , ( A , C ) , ( B , E ) , ( C , D ) to A ¯ since the distances between these node pairs are strictly less than 120 km. Then, we formulate the expanded arc set A ¯ = { ( s , O ) , ( s , A ) , ( s , B ) , ( O , A ) , ( O , B ) , ( O , C ) , ( A , B ) , ( A , C ) , ( B , C ) , ( B , E ) , ( C , E ) , ( C , D ) , ( E , D ) , ( E , t ) , ( D , t ) } . Combining V ¯ and A ¯ , we obtain the expanded network G ¯ V ¯ , A ¯ .
Sustainability 15 06588 g003
Figure 4. Example of a network graph in a toy instance composed of an O-D pair { O , D } , four intermediate nodes { 1 , 2 , 3 , 4 } and seven arcs { ( O , 1 ) , ( O , 2 ) , ( 1 , 3 ) , ( 1 , 4 ) , ( 2 , D ) , ( 3 , D ) , ( 4 , D ) } . The number on each arc represents the value of allocation variables.
Figure 4. Example of a network graph in a toy instance composed of an O-D pair { O , D } , four intermediate nodes { 1 , 2 , 3 , 4 } and seven arcs { ( O , 1 ) , ( O , 2 ) , ( 1 , 3 ) , ( 1 , 4 ) , ( 2 , D ) , ( 3 , D ) , ( 4 , D ) } . The number on each arc represents the value of allocation variables.
Sustainability 15 06588 g004
Figure 5. A hypothetical 25-node network.
Figure 5. A hypothetical 25-node network.
Sustainability 15 06588 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, L.; Zhou, B. Optimal Planning of Electric Vehicle Fast-Charging Stations Considering Uncertain Charging Demands via Dantzig–Wolfe Decomposition. Sustainability 2023, 15, 6588. https://0-doi-org.brum.beds.ac.uk/10.3390/su15086588

AMA Style

Wang L, Zhou B. Optimal Planning of Electric Vehicle Fast-Charging Stations Considering Uncertain Charging Demands via Dantzig–Wolfe Decomposition. Sustainability. 2023; 15(8):6588. https://0-doi-org.brum.beds.ac.uk/10.3390/su15086588

Chicago/Turabian Style

Wang, Luyun, and Bo Zhou. 2023. "Optimal Planning of Electric Vehicle Fast-Charging Stations Considering Uncertain Charging Demands via Dantzig–Wolfe Decomposition" Sustainability 15, no. 8: 6588. https://0-doi-org.brum.beds.ac.uk/10.3390/su15086588

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