Next Article in Journal
Statistics Using Neural Networks in the Context of Sustainable Development Goal 9.5
Previous Article in Journal
Sustainable Mining and Processing of Mineral Resources
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on the Location Selection Problem of Electric Bicycle Battery Exchange Cabinets Based on an Improved Immune Algorithm

1
School of Management, Shanghai University, Shanghai 200444, China
2
School of Business and Public Administration, Barry University, Miami Shores, FL 33161, USA
3
School of Cultural Heritage and Information Management, Shanghai University, Shanghai 200444, China
*
Author to whom correspondence should be addressed.
Sustainability 2024, 16(19), 8394; https://doi.org/10.3390/su16198394
Submission received: 19 August 2024 / Revised: 14 September 2024 / Accepted: 24 September 2024 / Published: 26 September 2024

Abstract

:
The rise of new energy technologies has accelerated progress towards sustainable development, and many companies are beginning to invest in renewable resource-related facilities. Electric bicycles have always been an important mode of green transportation; however, they also have problems such as slow charging, difficult charging, and that burning and short circuiting may occur during charging. Electric bicycle battery exchange cabinets effectively solve these problems by exchanging low batteries with full batteries instead of charging. However, current battery exchange cabinets face the problems of insufficient construction and unreasonable site selection. Therefore, this paper proposes a location selection model for electric bicycle battery exchange cabinets based on point demand theory, aiming to maximize rider satisfaction and the service capacity of exchange cabinets. The immune algorithm is introduced to solve the location model; however, the traditional immune algorithm has some problems such as poor stability and slow convergence. In this paper, the mutation process of the traditional immune algorithm is improved by introducing multi-point mutation, guided mutation, and local search. Finally, based on the data of electric bicycle riders in Shanghai, we verify that the location model based on point demand theory performs well on two objective functions of rider satisfaction and battery exchange cabinet service capability. We also expand the application of point demand theory to location models. Then, by conducting experiments with different parameter groups, through sensitivity analysis and convergence analysis, we verified that the improved immune algorithm performs better than the traditional immune algorithm in its accuracy, search accuracy, stability, and convergence.

1. Introduction

In the wake of the COVID-19 outbreak, economic development in many countries around the world has been hit hard. For example, China implemented lockdowns, bans, and quarantines to control the spread of the pandemic. Although these measures limited the spread of the virus, they also restricted the movement of people and goods, significantly hindering the country’s economic and social development [1]. Similarly, COVID-19 has had a significant impact on the energy sector. It has not only disrupted the supply and demand balance of traditional energy commodities, but it has also affected the supply chain of new energy commodities, such as batteries and solar photovoltaic panels, especially on the intercontinental trade routes between China and other countries. Restrictions on the movement of goods and people also led to the suspension of many projects and reduced demand for energy products, which has reduced the ability of companies and governments to invest in energy projects [2]. Therefore, businesses need to weigh the pros and cons of investing in energy projects more carefully. Now, in the post-pandemic era, the global economy is gradually recovering. This also brings some opportunities for businesses. On the one hand, the demand for energy commodities has increased sharply due to the movement of people and commodities after the pandemic. On the other hand, the reshuffling of the energy market will create a higher-quality market. Thus, the energy industry will face the most significant return shocks [3]. In addition, while the development of traditional energy sources will be under greater pressure after the pandemic, in contrast, investing in renewable energy infrastructure as a forward-looking strategy has proven to have great potential for development [4].
The post-pandemic era has accelerated the pace of sustainable development, and the economic recovery of countries around the world will inevitably have an impact on sustainable development in the future. Therefore, we should consider green economic recovery, promote economic recovery, and improve the environment together [5]. Green transportation has always been regarded as an effective way to improve the environment, as it can reduce the large amount of oil consumption, relieve pressure on resources, and promote energy conservation and emission reductions. In recent years, electric vehicles have been a research hotspot in achieving environmental sustainability. Similarly, the choice of transportation for green travel has become an emerging topic. Electric bicycles have always been one of the options for green transportation [6]. Compared to other traditional four-wheeled vehicles, electric bicycles alleviate traffic congestion and reduce exhaust emissions. Additionally, for shorter distances within 5 km, electric bicycles are a more convenient means of transportation. Compared to non-electric bicycles, electric bicycles have advantages in terms of environmental protection and travel speed. Even shared bicycles can potentially harm the environment, as their delivery and distribution are carried out by trucks, which produce harmful gases that pollute the environment. Clearly, researching issues related to electric bicycles will be a key focus for countries looking to enhance green transportation and energy consumption [7].
Electric bicycles serve as a low-carbon, environmentally friendly, convenient, and affordable means of transportation, making them an excellent advocate for energy conservation, emission reduction, and green travel [8]. In China, the development of electric bicycles has clearly become mature. The government actively promotes electric bicycles as a mode of travel while considering the economic, social, and environmental benefits brought by electric bicycles, and promotes relevant policy research to encourage the exploration and development of shared electric bicycles and Internet-based electric bicycles while ensuring urban traffic safety. As of 2023, the number of electric bicycles in China has exceeded 450 million, with bicycle delivery personnel being one of the main user groups [9]. Although electric bicycles are favored for their ease of operation, economic and environmental benefits, and maneuverability, they also face issues such as slow charging, difficulties in charging, and risks of short circuits and burning during charging. These problems significantly hinder the development of electric bicycles [10].
Electric bicycle battery exchange cabinets combine Internet of Things technology with the sharing economy mindset. They are placed on street corners and in alleys, resembling rectangular cabinets over two meters tall with typically 12 compartments. These compartments can store batteries and charge them. Riders of electric bicycles can connect to the cabinet via a mobile app, allowing them to open a compartment and exchange their low battery with a full battery within three minutes. For delivery and logistics workers who ride long distances daily, ranging from 100 to 150 km, the typical lithium battery electric bike has a range of about 50 km per charge. Therefore, riders usually need to charge their batteries two to three times a day on average [11]. Slow charging batteries takes 4 to 6 h, while fast charging times range from 30 min to 2 h. This significantly impacts the riders’ work. Therefore, riders increasingly need public battery exchange cabinets that allow them to quickly replace low batteries with full batteries throughout their daily travels, thereby replacing the need for traditional charging [12]. At the same time, this battery exchange approach effectively avoids various safety hazards that may occur during the charging process of electric bicycles. While the cabinet is charging the batteries, it can also instantly report any abnormal battery statuses to backend supervisors for immediate handling. This capability helps effectively eliminate safety risks and prevents fire incidents [13].
After 2012, China’s food delivery industry entered a period of rapid growth with frequent capital investments and intense competition among enterprises. By 2016, the market was predominantly dominated by two major players, Elema and Meituan [14]. In the same year, a group of Chinese enterprises began independently developing battery exchange cabinets specifically for charging electric bicycle batteries. Companies providing battery exchange services also emerged within five years. As of now, products like those from these manufacturers like Tower Energy, Hello Bike, Didi Chuxing, and Qingju Bike in China, which represent battery exchange cabinets, have been launched into the market. They have effectively alleviated the battery exchange demand for food delivery riders [15]. Especially in the post-pandemic era, the demand for food delivery and outdoor activities continues to rise, leading to a gradual increase in demand for battery exchange cabinets. However, within the battery exchange cabinet industry, there is intense competition, with many enterprises currently experiencing a surge in capital influx. Many cities still face issues such as insufficient construction of battery exchange cabinets and unreasonable site selection, which greatly hinder the development of the electric bicycle industry. Therefore, it is crucial to address these issues by constructing an adequate number of proper battery exchange cabinets to meet the battery exchange demand of electric bicycles. This will further mitigate the problems of slow and difficult charging, enhance the purchasing willingness of electric bicycle riders, and promote the development of the electric bicycle industry.
The significance of this study lies in providing a solution for the site selection of battery exchange cabinets for electric bicycles, manifested in two main aspects. Firstly, this paper uses kernel density analysis to identify the hot spot distribution of rider demand. Secondly, we build a site selection model based on point demand theory, which aims to maximize rider satisfaction and maximize battery exchange cabinet service capacity. Compared to models based on a single objective, the objective function in this study is more comprehensive and reasonable, offering a method for determining demand volume based on point demand theory. Finally, we innovatively propose an improved immune algorithm that combines two different types of mutation operators and local search. It is compared with the traditional immune algorithm in terms of convergence speed, stability, and search accuracy.

2. Materials and Methods

2.1. Related Theory of Location Selection Method

The site selection theory adopted in this paper is a point demand-based site selection theory. This theory posits that riders’ demand originates at specific nodes, typically dividing all nodes within the site selection range into demand points and candidate points. By analyzing the relationship between demand points and candidate points, the optimal site selection scheme is determined. Research based on point demand site selection theory mainly includes P-median site selection studies, coverage site selection studies, and P-center site selection studies.
The P-median site selection theory was initially proposed in 1970 and is formulated as a mixed-integer model [16]. In recent studies, scholars have used and improved the P-median theoretical model to different degrees to make the model more suitable for solving the location problem. Murad et al. studied the optimization of healthcare service locations. They propose an MC-P-Median model. The model first used the Maximum Coverage (MC) location problem model to prioritize demand locations, then calculated the optimal hospital location through the P-median function. Finally, the model was verified using the Jeddah hospital location problem as an example [17]. M H Ramadhan et al. [18] studied warehouse locations by considering factors such as distance, customers, and cost. They used the combination of set coverage and P-median to construct the location model. Then, by comparing the size of the target value obtained by the two methods, they reached the conclusion that the P-median method has a better solving effect. Janjić et al. [19] explored the location problem of electric wake-up charging piles. In the case of determining the number of charging piles, they set up a new location model by improving the P-median iteration method and combining it with a weighted analytic hierarchy process. The objective function considers the installation cost of the charging pile, the distance of the car to the charging pile, and some actual conditions of the selected location. Finally, the effectiveness of the proposed method is verified by taking Nish City, Serbia, as an example. Hanif [20] studied the application of the centroid method and the P-median method for the site selection of a single non-averse facility for a utility store company. They aimed to minimize transportation costs and determine the optimal location for a warehouse.
The principle of covering location is usually divided into two types in the study of location problems: one is set cover, and the other is maximum cover. When studying the location of electric vehicle charging facilities, Davidov et al. [21] considered the driving path of electric vehicles, the time required for charging, and the parking time after users arrive at the charging facilities, and then established a discrete location model by using the set to cover the original.
Similarly, Zhou et al. [22] also studied the location of EV charging stations. From the perspective of station construction cost, they divided the cost into comprehensive economic cost and environmental cost. The site selection model based on maximum coverage site selection theory is established. Then, taking Ireland as an example, through the genetic algorithm, the model is solved and analyzed. Fang et al. [23] studied the site selection problem for public electric vehicle charging stations, considering the limited range and charging demands of pure electric vehicles. They established a Maximal Covering Location Problem model and solved it using a dual-level mathematical programming approach and a genetic algorithm-based procedure. Pan et al. [24] considered the spatial distribution of demand points, service coverage, and construction costs in different areas to influence the layout of charging pile site selection. They aimed to minimize the number of candidate points and costs within the service coverage area. They established a Set Covering Model and used the entropy weight method and greedy algorithm to solve it. Huang [25] studied the layout optimization of public charging stations in Jin Shui District, Zhengzhou City. The study utilized charging station data, parking lot data, statistical yearbook data, population data, and urban road network data. They employed a Set Covering Model and a P-median model with capacity constraints to optimize the layout of public charging stations in Jin Shui District.
The P-center location problem was proposed by Hakimi [26], based on the relationship between feasible candidate points and grid vertices, as well as inherent constraints. The P-center problem is divided into a vertex center and an absolute center. To provide a basis for further research and improvement of the P-center model, Liu et al. [27] proposed an iterative solving framework for computing the lower bounds of the P-center model. With the development of the times and the diffusion of social networks, to better handle the impact of social factors on facility location, Ma et al. [28] proposed a partitioning algorithm based on Voronoi diagrams. This algorithm aims to establish social constraints and integrate them with spatial constraints, thereby developing an enhanced P-center location model. Yi et al. [29] considered the impact of agglomeration effects in mountainous urban areas and proposed an improved P-center model to determine the optimal locations for charging facilities. Zhou et al. [30] analyzed the relationship between agglomeration effects in mountainous urban areas and the predicted distribution of charging demand, fitting correlation parameter values. They integrated these values with an improved P-center location model, using annual construction costs, electric vehicle operating costs, annual maintenance costs, electricity purchase costs, and expanded distribution network costs as optimization indicators. This was used to establish a site layout model for charging stations in mountainous urban areas.
In addition to the point demand theories, many scholars also study site location from the perspective of business profitability or serving users, which is also part of the research on point demand theory.
Scholars consider the profitability of enterprises in the location problem, which is usually achieved by reducing costs, such as the cost of site construction or the cost in the subsequent operation stage. Fareed et al. [31] considered the impact of electric vehicle load on power loss in the distribution network. They combined the land cost index and the electric vehicle charging demand index to construct the charging station investor index. Then, based on the decision index, an optimization model aiming to achieve the minimum power loss cost was constructed, and the optimal location of the charging station was determined. Gao et al. [32] considered the impact of existing charging stations on the location of new charging stations. By introducing intermediate variables, they constructed a segmented nonlinear function to describe the relationship between existing charging stations and new charging stations, considering the total investment cost of charging stations and the distance between users and charging stations. Finally, they developed a location model that could be based on adjustments to existing charging stations in each area. Min et al. [33] aimed to minimize daily food delivery costs using pure electric trucks. They built a first-stage site selection model for electric truck charging stations and defined relevant parameters to be input into the second-stage model. The second stage model is based on the fixed capacity model of the pure electric truck charging station, considering the annual comprehensive cost of the charging station. By integrating Phase 1 and Phase 2 models, they determined the site layout and capacity of all-electric truck charging stations.
From the perspective of serving users, considerations often include their time costs and convenience. Zhang et al. proposed a charging station location problem that considers user preferences and waiting times, formulated as a multi-objective two-level programming model. Zhang et al. [34] studied the location of charging piles and considered users’ preferences for charging piles and the waiting time for charging. They built a multi-objective, two-level programming model. The first-level model aims to minimize total costs to determine the locations and number of charging piles, while the second-level model aims to minimize total travel time to match the user’s preference. Xu et al. [35] thoroughly considered the impact of transportation convenience on users driving electric vehicles and the acceptance capacity of the distribution network. They established a fixed-capacity charging station location model that incorporates both convenience and acceptance capacity considerations. They solved this model using a combined approach of particle swarm optimization and Voronoi diagrams. Tran et al. [36] are studying the location of charging piles. They give priority to the driving distance between users and charging. Secondly, they combine the location model of charging pile construction design. Zhang et al. [37] also considered the user’s driving time to the charging pile when studying the location of charging stations, and proposed the concept of social cost, including the construction cost and operating cost of the charging pile. Finally, they built a site selection model and used an improved firefly algorithm to solve the model.
In summary, this paper considers the user experience of electric bicycle riders, with the goal of maximizing rider satisfaction and battery exchange cabinet service capacity. A location model based on point demand theory is established. Compared with the single objective location model, the objective function of this paper is more comprehensive and reasonable.

2.2. Location Selection Method-Related Algorithm

Some scholars use genetic algorithms to solve the point demand model. Although it has advantages in solving problems involving mixed discrete, continuous, and integer design variables, the immune algorithm has characteristics similar to the ability of the biological immune system to produce various specific antibodies and has advantages over genetic algorithms, such as preventing degeneration, accelerating convergence speed, and simplifying operators. Therefore, immune algorithms are more suitable for solving various types of location models, vehicle scheduling models, and many other different types of models [38]. Liu et al. [39] studied the location problem of medical waste treatment facilities and solved the multi-objective optimization model of the facility location using a traditional immune algorithm. Traditional immune algorithms suffer from slow convergence and low search accuracy. To address these issues, various improvements have been made to this algorithm. Yang et al. [40] combined the immune algorithm with an ant colony algorithm, innovatively proposing a new immune ant colony optimization (IACO) algorithm. They applied this algorithm to solve multi-objective task allocation problems for warehouse robot groups. The experimental results demonstrate that integrating these two algorithms can effectively reduce computational redundancy in immune algorithms and improve convergence speed. Xu [41] improved the immune algorithm by extracting methods from artificial immune theory, such as extracting and injecting vaccines. Based on the improved immune algorithm, he solved site selection problems and proved that the algorithm can reach the optimal solution faster through iteration. Wang et al. [42] introduced a new improvement to the classical immune algorithm by incorporating a drug-assisted process. This approach reduces the size of the search space, reaching the optimal solution of the model faster. They validated the feasibility of this approach by solving test functions. Xu et al. [43] introduced several new operators into the immune algorithm to modify the process of population generation, updating, and cross-variation. To improve algorithm performance and adapt it for solving their constructed site selection model for emergency medical facilities, Zhou et al. [44] made three specific design improvements to the immune algorithm. Firstly, they optimized the generation of the immune algorithm population based on the specific characteristics of the emergency medical facility site selection model. Secondly, they designed an antibody evaluation function. Lastly, they enhanced the crossover stage by designing two neighborhood search mechanisms to reduce antibody similarity. This ensured diversity among antibodies in each iteration, facilitating the algorithm to achieve global optimal solutions and preventing premature convergence.
In this paper, the improved immune algorithm is used to solve the location model. Scholars have improved the traditional immune algorithm in two main ways: one is to combine the immune algorithm with other intelligent algorithms and the other is to refine some internal calculation processes of the immune algorithm. This paper enhances the performance of the immune algorithm through the second approach, specifically by improving the mutation phase of the algorithm. This improvement involves introducing two types of mutation operators and operators for local search.

3. Problem Description

3.1. Basic Problems

The basic problem of this paper is how to determine the number and locations of electric bicycle battery exchange cabinets. As shown in Figure 1, the steps to solve the basic problem are mainly divided into two stages and four associated problems. The first stage is to collect battery data through battery sensors to determine the battery exchange demand of e-bike riders in the area, which is the content of the first associated problem. In the second stage, based on demand points, the hot spots of demand points are identified through nuclear density analysis, which is the content of the second associated problem. Then, based on point demand theory, a location model of electric bicycle battery exchange cabinets is constructed. This model is a bi-objective model that considers two factors that affect the selection of electric bicycle battery exchange cabinets, which is the content of the third associated problem. Then, to clearly indicate the location scheme, this paper also needs to delimit the location range of battery exchange points in the area, which is the fourth associated problem. Finally, the siting of these cabinets can also influence rider travel trajectories to some extent, leading to variations in the data collected by battery sensors correspondingly. It is essential to address the following associated problems:
  • Identifying the battery exchange demand among electric bicycle riders in the area.
This paper collects riding data of electric bicycles through internal battery sensors, including battery charge levels and current location information. Batteries may require exchange at any location, but the willingness to exchange varies with the battery’s current charge level. Therefore, we define the parameters r and v. The parameter r is the demand coefficient, the parameter E is the battery charge at a particular moment, and the relationship between the two is the parameter r = 100 E . Then, we perform a weighted average for each battery based on the demand coefficient, as shown in Equation (1). The point X , Y represents the coordinates of the battery at a specific moment. The point X ¯ , Y ¯ denotes the weighted average coordinates of the battery across all moments. Each battery represents a rider and corresponds to a demand point within a specific area range.
X ¯ , Y ¯ = r X N r , r Y N r
2.
Analyzing the distribution differences and clustering of demand hotspots.
Kernel density analysis is used to measure and analyze the density of point or line features in their surrounding neighborhoods. It is one of the methods used to determine population density distribution in urban areas. Compared to traditional population density analysis, the advantage of kernel density analysis lies in its ability to intuitively reflect the distribution of discrete measurement values within continuous areas. It also overcomes the issue of uniform density values within statistical units. Therefore, this paper employs kernel density analysis to identify differences and clustering patterns in rider density hotspots within urban areas. According to the kernel density estimation method, we generate evenly spaced grid points within the coordinates of rider demand points to cover all rider locations. The density value for each grid point is computed using the kernel density estimation formula.
For the kernel function selection, the Gaussian function is commonly used as the kernel function, as shown in Equation (2):
K μ = 1 2 π e 1 2 μ 2
In this paper, the bandwidth parameter h is set to control the width of kernel function, which determines the smoothness of kernel function. Smaller bandwidths produce more detailed density estimates (which can lead to overfitting), while larger bandwidths lead to smoother density estimates (which can lead to underfitting). Silverman’s law is commonly used to calculate bandwidth, as shown in Equation (3). The parameter σ represents the variance in the sample data, and the parameter n represents the number of samples in the data. Parameter a in Silverman’s rule adjusts the bandwidth to balance between bias and variance. The term 1 / 5 acts as a scaling factor to adjust the bandwidth h with respect to the sample size n .
h = a σ n 1 5
The density estimate at each grid point (x, y) is obtained by averaging the weighted kernel function values of all data points. The equation for the density estimates at each grid point x ,   y is shown in Equation (4). The point x i , y i represents the demand points for each rider. Then, based on the density estimates, a heatmap of the density distribution is plotted to identify hotspots where rider density is high.
f x , y = 1 n h 2 K ( x x i 2 + y y i 2 h )
3.
Factors influencing the siting of electric bicycle battery exchange cabinets
Nowadays, China is in the early stages of development for electric bicycle battery exchange cabinets. The service of electric bicycle battery exchange has not yet been fully popularized. From the point of view of promoting services and increasing riders’ willingness to replace batteries, when selecting sites for battery exchange cabinets, priority should be given to considering the riders’ perspectives and factors that encourage growth in the number of battery-exchanging riders. Therefore, this paper primarily considers the following two main factors.
The first factor is riders’ satisfaction with battery exchange cabinets, which refers to whether there is a nearby cabinet when riders need to exchange batteries, and the closer the cabinet is, the higher the riders’ satisfaction. The satisfaction function of electric bicycle riders describes the relationship between objective distance and riders’ subjective evaluation. The satisfaction function of electric bicycle riders is presented in Equation (5). According to point demand theory, the parameter i is the demand point, which is the rider’s coordinate point. The parameter j denote the candidate point, which is the cabinet’s candidate coordinate point. The parameter L i indicates the farthest distance from the power change cabinet when the rider is most satisfied. The parameter U i indicates the maximum change cabinet distance the rider can tolerate, and the parameter d i j is the distance from demand point to candidate point. The function F ( d i j ) denotes the satisfaction evaluation value of electric bicycle riders at the demand point.
F d i j = 1 , d i j < L i 1 2 + 1 2 cos π U i L i d i j U i L i 2 + π 2 , L i < d i j < U i 0 , d i j > U i
The second factor is the service capacity of a battery exchange cabinet, which indicates the number of demand points that fall within the scope of the service. The higher the service capacity of a battery exchange cabinet, the more riders it can serve. Unlike power infrastructure such as power plants and substations, electric bicycle battery exchange cabinets serve as public infrastructure for the public. In addition to meeting urban layout and grid planning requirements, they also consider the maximum distance riders can travel to exchange batteries when their electric bicycle battery levels are low. In this paper, this distance is defined as the service radius R . If a coordinate point is selected as the location for building a battery exchange cabinet and it is the closest cabinet to the demand point, and if the distance between the demand point and coordinate point is less than the service radius R , then the coordinate point is considered to serve the demand point. Therefore, a higher service capacity of the battery exchange cabinet means it can serve more riders, leading to greater profitability for the enterprise.
4.
The delineation of battery exchange point siting ranges in space
According to point demand theory, we simplify the site selection area into a rectangular region of specified length and width, establishing a coordinate system with the bottom-left corner of the rectangle as the origin. Then, we map the latitude and longitude of the demand points into the system. Every coordinate point in space is considered a potential candidate point. The problem in this paper transforms into a site selection issue where p battery exchange cabinets need to be selected from J candidate points within a specific rectangular area. Site selection solutions must align with rider satisfaction and the service capacity of the battery exchange cabinets.

3.2. Basic Assumptions

To prioritize maximizing rider satisfaction and maximizing the service capacity of battery exchange cabinets for electric bicycle riders, we establish a siting model. To facilitate the solution of the model, we make the following assumptions:
  • According to the principle of proximity, riders prefer to go to the nearest battery exchange cabinet to their current location.
  • In the study area, rider demand points remain the same every day.
  • Due to the small footprint of battery exchange cabinets, they can be set up in most streets and alleys. Therefore, this paper does not consider whether the selected locations for setting up battery exchange cabinets meet the conditions for construction.
  • The number of battery exchange cabinets set up in the network is not considered, i.e., ensuring they can meet the battery exchange demand of all riders within their service range.
  • Battery data are collected from internal sensors, which accurately reflect riders’ cycling conditions.
Some parameters in the model have been simplified and replaced accordingly, as follows:
  • The distance between riders and battery exchange cabinets is their Euclidean distance, shown in Equation (6). The point x i , y i and the point x j , y j are the coordinates of the demand point and candidate point, respectively.
    d i j = x i x j 2 + y i y j 2
  • During data collection, battery levels are uniformly distributed within a specified range. We know that riders have no intention to exchange batteries when their battery level is above a certain threshold. Therefore, we assume that riders have no battery exchange demand when their battery level is greater than 60%. Data points where the battery level exceeds 60% are excluded from the original dataset. Although this exclusion may cause slight deviations in the latitude and longitude of demand points, we believe it will not significantly affect the final results of this study.
  • This paper assumes that R is the maximum tolerable service radius for all e-bike riders.

4. Model Construction and Algorithm Design

4.1. Model Construction

The decision variables in this paper are Z j and Y i j , both of which are binary variables (0 or 1). If the variable Z j is 1, it indicates that the candidate point is selected as the site for establishing a battery exchange cabinet. If the variable Z j is 0, it means the candidate point is not selected as a cabinet location. If the candidate point serves the demand point, then the variable Y i j is 1; otherwise, the variable Y i j is 0. The demand point must be within the service range of the candidate point. Therefore, the variables Y i j , Z j ϵ 0 , 1 i ϵ n i , j ϵ n j . The model’s symbols and indicators are shown in the following Table 1.
Based on the above descriptions, the expression for establishing a battery exchange cabinet siting model based on satisfaction and service capacity is as follows.
Objective function:
m a x f 1 = j = 1 J i = 1 I F d i j Y i j Z j I
m a x   f 2 = j = 1 J i = 1 I Y i j Z j I
Constraint conditions:
j = 1 J Y i j = 1   i ϵ n i , j n j
j = 1 J Z j = P   j ϵ n j
d i j R Z j = 1 , Y i j = 1
Z j , Y i j ϵ 0 , 1 i ϵ n i , j ϵ n j
In the siting model, Equation (7) represents maximizing average rider satisfaction. Equation (8) represents maximizing average battery exchange cabinet service capacity. Due to redundancy in the decision variables, certain constraints are imposed as follows:
(a)
It is ensured that each battery exchange cabinet’s candidate point serves multiple demand points and that each demand point is served by only one candidate point of a battery exchange cabinet. The constraints are formulated as in Equation (9).
(b)
In Equation (10), this paper selects p candidate points as the coordinate points set up by the battery exchange cabinet.
(c)
If the candidate point serves the demand point, then the distance between the demand point and the candidate point is less than radius R , as shown in Equation (11).

4.2. Improved Immune Algorithm

The improved immune algorithm addresses the challenging issue of premature convergence during iteration processes while maintaining population diversity. However, traditional immune algorithms suffer from slow convergence speed, low search accuracy, and poor stability, especially evident in large-scale computational problems like battery exchange cabinet siting. To tackle these issues, the following improvements were made to the traditional immune algorithm:
(a)
A multi-point mutation operator is introduced in the mutation process. Compared with the single point variation and small number of mutations in the traditional immune algorithm, more changes are introduced in the generation of new candidate solutions in the early stage of population evolution, which increases the chance of generating the optimal solution in the iterative process and speeds up algorithm convergence. In the later stage of population evolution, the algorithm is prevented from falling into the local optimal solution and prematurely converging.
(b)
Guided mutation is added to the mutation operator of the traditional immune algorithm, mainly to improve the orientation of the mutation at the initial stage of population evolution and guide the antibody towards the optimal solution region, which can greatly increase the initial convergence of the algorithm compared with the traditional immune algorithm.
(c)
In the improved immune algorithm, the solution is fine-tuned and optimized by local search and setting termination conditions. In the late stage of population evolution, traditional immune algorithms often fall into the dilemma of local optimization. Local search can greatly improve this situation and speed up the convergence of the algorithm.
Based on the characteristics of the constructed model, an immune algorithm for battery exchange cabinet siting was designed, incorporating population generation, affinity functions, two types of mutation operators, and a local search operator. The specific calculation process of the improved immune algorithm is shown in Figure 2, described in Algorithm 1.
Algorithm 1: Improved immune algorithm
Input: Population size P s , antibody characteristics (number of battery exchange cabinets) P , generations of population evolution g , clone rate r c , mutation ratio r m .
Output: Site selection scheme for electric bicycle battery exchange cabinets
Step 1. Initialize Population
    According to the point demand theory, among the candidate points with a total amount of j , we select the candidate point with a number of P as the location of the battery exchange cabinet. These candidate sites form an antibody. Therefore, antibodies are a feasible solution to the location problem, and the number of candidate points composed of antibodies is P . These candidates have zero similarity to each other. A population is formed by the clumping together of q antibodies. The immune algorithm iterates the optimal antibody through continuous population multiplication. During antibody generation, we define the distance between each candidate point and at least one demand point to be within R , which can greatly improve the convergence speed of the algorithm.
Step 2. Population evolution
    Evolve the population over multiple generations, with the specific number of iterations determined by parameter g . In each generation, perform the following operations:
    Step 2.1. Affinity calculation
        The electric bicycle battery exchange cabinet siting model is a multi-objective nonlinear model, which can be handled using linear weighting methods. The affinity value is defined as w 1 f 1 + w 2 f 2 , where f 1 and f 2 are the normalized objective values of rider satisfaction and the service capacity of battery exchange cabinets, respectively. w 1 and w 2 are the weights for each objective, satisfying w 1 , w 2 ϵ 0 , 1 , w 1 + w 2 = 1 . The normalization of f 1 and f 2 is calculated as shown in the equation below. m a x f and m i n f represent the maximum and minimum values that f can attain.
f = f m i n f m a x f m i n f
    Step 2.2. Update optimal antibody
        Identify the antibody with the highest affinity from the population. Compare the affinity of this antibody with the current optimal antibody and update the optimal site selection scheme if the affinity exceeds the current best solution.
    Step 2.3. Clonal elite antibody
        Sort antibodies based on their affinity values. Select P s × r c proportion of antibodies as elite antibodies. Clone each elite antibody, generating a few clones equal to ( 1 r c ) / r c .
    Step 2.4. Multi-point Mutation and Guided Mutation
        Perform multi-point mutation on half of the cloned antibodies, selecting new candidate points to replace the old ones. Conduct guided mutation on the other half of the cloned antibodies, mutating based on the current optimal solution to increase directionality.
    Step 2.5. Local Search
        Perform a search on elite antibodies by aggregating searches within the elite antibody domain to enhance the directionality of mutation.
    Step 2.6. Termination Condition Evaluation
        At the end of each generation, check if the termination condition is met. If satisfied, terminate the algorithm early to reduce unnecessary computation.
    Step 2.7. Update Population
        Combine the newly generated clone antibodies of selected antibodies into a new population. Add elite antibodies to the new population to ensure that the optimal solution is retained within the population.
Step 3. Population iteration
    Repeat Step 2 until the number of iterations reaches the generations of population evolution g .

4.3. Multi-Point Mutation, Guided Mutation, and Local Search

Multi-point mutation is an extension of single-point mutation in traditional immune algorithms, where mutations occur at multiple gene positions to enhance individual diversity and algorithm exploration capabilities. By mutating at multiple positions, it becomes possible to alter the genotype of individuals over a larger range, thereby improving the algorithm’s ability to escape local optima. In this study, half of the cloned antibodies are subjected to multi-point mutation. Set the multi-point mutation rate as r m , and calculate the number of mutated antibody genes as P r m . Randomly select P r m antibody genes, ensuring that the randomly generated points fall within the service range of the battery exchange cabinet, and replace the old points with the new ones. The process of multi-point mutation is illustrated in Figure 3, where we assume a randomly selected set x 2 , y 2 , x t , y t comprising P r m genes. Then, x 2 , y 2 , x t , y t is randomly generated to replace the original gene set, thereby creating new antibodies.
Guided mutation builds upon traditional single-point mutation by targeting the remaining half of clones that have not undergone multi-point mutation. It sets a mutation magnitude c with a tendency towards the center, making the mutation process more targeted and capable of converging faster to the global optimum. The mutation process is illustrated in Figure 4, where, for example, if a gene point x t , y t is randomly selected, and if x t is greater than the average value x ¯ , it reduces by c to move closer to the center; conversely, if x t is less than x ¯ , it increases by c . Similarly, this applies to y t .
The article adopts local search for elite individuals, playing a role in fine-tuning and optimizing solutions within the improved immune algorithm, accelerating the convergence speed of the algorithm and escaping local optima. The local search algorithm process is shown in Figure 5. A gene point is randomly selected, assuming x t , y t with an affinity value f t . Within its neighborhood range, with a change magnitude of a, x t , y t transforms into eight neighboring solutions: x t + a , y t , x t a , y t , x t , y t + a , x t , y t a , x t + a / 2 , y t a / 2 , x t a / 2 , y t a / 2 , x t + a / 2 , y t + a / 2 , x t a / 2 , y t + a / 2 . If the affinity of a neighboring solution exceeds that of the current solution, it accepts replacing the current solution. Repeat this process until reaching the termination condition of the search. There are two kinds of termination conditions: one is that the optimal solution is found through local search, the other is that the local search has been carried out for each position. Typically, the algorithm terminates when one of these conditions is met.

5. Computing Example

5.1. Data Pre-Processing

The example analysis in this paper takes electric bicycles in Shanghai as the research object. The data were provided by a leading Chinese company related to shared batteries. The data collected by the sensors in the electric bicycle batteries under the company’s name mainly include information such as battery power, status, and location. Second, we simplified the land area of Shanghai into a rectangle 120 km long and 100 km wide. Then, each battery is regarded as a battery exchange demand point, and after calculating the weighted average, it is found that there are 19,000 demand points. The distribution of demand points in Shanghai is shown in Figure 6. The area is divided into 16 sections corresponding to the 16 administrative districts of Shanghai. The darker the color, the more demand points there are in that area, with the deepest brown representing the area with the highest number of demand points. Conversely, lighter colors indicate fewer demand points. The demand points are concentrated in the central and eastern districts.

5.2. Nuclear Density Analysis

A kernel density analysis was conducted for Shanghai, using the natural breaks classification method to categorize different levels of rider density. The rider density distribution is illustrated in Figure 7. Darker areas represent higher rider density. In comparison with Figure 6, it can be observed that the densest areas are also concentrated in the central and eastern parts of the city. These areas coincide with the most economically developed regions of Shanghai.

5.3. Theoretical Model Analysis

Nowadays, the electric bicycle battery exchange service industry is in its early stages of development, and battery exchange services are not yet widespread. At this stage, the objective value established by the site selection model serves as a crucial criterion for evaluating the quality of electric bicycle battery exchange cabinet location schemes. We set the variable parameters. The parameter L i is 2   km , and the parameter U i is 1   km . The service radius R of the battery exchange cabinet is 2   km . Parameters ω 1 and ω 2 are both set to 0.5. The parameters for the immune algorithm are set as follows: Parameters n and g both are 1600. The parameter p is 800. The parameter r c is 0.2 and the parameter r m is 0.1. Then, to determine the number of battery exchange cabinets to construct, sensitivity analysis of the site selection schemes is conducted, as shown in Table 2. As the number of battery exchange cabinets increases, the objective value of the site selection schemes improves continuously from the users’ perspective. In addition, an average of two battery exchange cabinets are installed at each exchange cabinet site, and each exchange cabinet can store 12 batteries. Given the demand for 19,000 batteries, an average of 791 exchange cabinets are required. To balance these factors, we chose a scheme that maximizes the difference value of comprehensive objective f . Finally, we selected 800 candidate sites for setting up electric bicycle battery exchange cabinets.
Using asterisks to represent battery exchange cabinets and shading at the bottom to indicate the density range of demand points, Figure 8 illustrates the site selection plan for electric bicycle battery exchange cabinets. Combining Table 2 and Figure 8, when the parameter P is 800, the user satisfaction reaches 92.90%, and the battery exchange cabinet service capacity reaches 98.6%. The comprehensive objective function value also exceeds 95%. It shows that most riders are very satisfied with the location scheme, and most riders are within the scope of battery exchange cabinet service. Meanwhile, we can also observe from Figure 8 that the exchange cabinets effectively cover all areas of battery exchange demand, enhancing service capability. The plan meets the battery exchange demand of all riders, and the greater the rider density of the area, the more battery exchange cabinet location points. Therefore, riders’ demand for battery exchange is met and riders’ satisfaction is improved. The results show that the point demand theory method of combining the two objectives is very effective in the application of a site selection scheme. It is also an extension of the application of the point demand theory to the location problem.

5.4. Comparison of the Two Algorithms

In this paper, sensitivity analysis and convergence analysis are used to compare the accuracy, search accuracy, stability, and convergence of the two models. To simplify the analysis in Section 5.4, we use the value of the comprehensive objective f to replace the value of the two objective functions, and the comparative analysis of the algorithms in this section only refers to the value of the comprehensive objective f . Similarly, to make the verification more accurate, this paper sets 2 × 3 different parameter combinations, the parameter r c is set to 0.1 or 0.2, and the parameter r m is set to 0.1, 0.2, or 0.3. The two algorithms are each run 10 times under each parameter combination, totaling 2 × 3 × 10 × 2 = 120 .
(1)
Sensitivity analysis
We compared the accuracy, search accuracy, and stability of the algorithms through sensitivity analysis, and calculated the optimal value, average value, and standard deviation of the comprehensive value of objective f combined with the results of each run, as shown in Table 3. Accuracy is determined by comparing the optimal value of the two algorithms, where the larger the optimal value, the more accurate the algorithm. Search accuracy is measured by the difference between the average value and the optimal value. The smaller the difference, the closer the average solution is to the optimal solution, the greater the probability of the algorithm reaching the optimal solution, and the higher the search progress of the algorithm. The stability of the algorithm is measured by standard deviation, where the smaller the standard deviation, the more stable the algorithm.
First, we find that under the conditions of r c = 0.2 and r m = 0.1, the optimal value and average value of the two algorithms are the best. Under these conditions, the accuracy of the improved immune algorithm is increased by 1%, the search progress is increased by 46%, and the stability of the algorithm is increased by 34%. Secondly, we comprehensively calculated the average values under the six combinations, and the accuracy of the improved immune algorithm increased by 0.97%, the search progress increased by 32%, and the stability of the algorithm increased by 28%. We also found that under the conditions of r c = 0.1 and r m = 0.3, the stability of the improved immune algorithm is 28% lower than that of the traditional immune algorithm. Finally, the above analysis shows that the improved immune algorithm is superior to the traditional immune algorithm in the above three aspects.
(2)
Convergence analysis
In this paper, the convergence of the two algorithms is compared through convergence analysis. In Equation (13), the parameter g represents the number of generations the algorithm needs, the parameter t represents the index of the number of generations, values f b e s t t + 1 and f b e s t t are the optimal values under t + 1 and t generations, and then the rate of change of each generation needs to be summed. If the algorithm converges earlier, f b e s t t will be lower, and thus the convergence rate will be higher. Similarly, if the algorithm’s final optimal target value is larger, the difference between values f b e s t t + 1 and f b e s t t will be larger, and the convergence rate will be improved. Therefore, the convergence rate can effectively measure the convergence speed of the algorithm and the final optimal value span. In this paper, the average convergence rate of the 10 algorithm iterations run under each combination is shown in Figure 9. We find that both algorithms perform best when r c   = 0.2 and r m = 0.1, with convergence being improved by 7.29%, and perform the worst when r c = 0.1 and r m = 0.3. Moreover, the average convergence rate of the improved immune algorithm under the six parameter combinations is higher than that of the traditional immune algorithm, with an average increase of 12.4%, which proves that the convergence of the improved immune algorithm is better than that of the traditional immune algorithm under the influence of multi-point variation and guided variation.
C o n v e r g e n c e   R a t e = t = 1 g 1 f b e s t t + 1 f b e s t t f b e s t t

6. Conclusions and Future Research

This paper proposes an electric bicycle battery exchange cabinet site selection model that considers rider satisfaction and the service capability of exchange cabinets and expands on existing literature. This paper is divided into three steps. Firstly, the study uses demand point theory to establish the site selection model, considering rider satisfaction and the service capability of the exchange cabinets. Secondly, it introduces an improved immune algorithm for model solving and enhances the algorithm by adding two mutation operators and one local search operator, thereby improving convergence speed, stability, and search accuracy. Thirdly, it provides a computational instance for electric bicycle riders in Shanghai, validating the advantages of the improved immune algorithm through comparison with the traditional immune algorithm. This offers enterprises optimal solutions for exchange cabinet site selection.
Through example calculations, we verify that the location model based on point demand theory performs well on two objective functions of rider satisfaction and battery exchange cabinet service capability. The battery exchange cabinets effectively cover all areas where battery exchange demand points are distributed, achieving the maximum service capability objective for the cabinets. Moreover, the cabinets are distributed more heavily in areas with higher demand points, maximizing rider satisfaction. We expand the application practice of point demand theory in the location selection model. Then, experiments with different parameter groups, through sensitivity analysis and convergence analysis, verify that the improved immune algorithm performs better than the traditional immune algorithm in accuracy, search accuracy, stability and convergence.
In future research, various simplifying assumptions (as discussed in Section 3) could be relaxed to better represent the problem. For instance, actual route distances based on POI data could replace Euclidean distances between riders and exchange cabinets as modeled in this study. Future work will aim to more scientifically and practically design applicable site selection models or more efficient algorithms. Furthermore, considering more practical factors in exchange cabinet site selection processes, such as enterprise perspectives on exchange cabinet construction costs, operational costs, and profitability, will be crucial.

Author Contributions

Conceptualization, Z.Z. and S.Y.S.; methodology, W.Y. and X.Y.; validation, W.Y. and X.Y.; writing—original draft preparation, W.Y. and X.Y.; writing—review and editing, Z.Z., W.Y. and X.Y.; visualization, W.Y.; supervision, Z.Z.; project administration, Z.Z.; funding acquisition, Z.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are provided by the enterprise. Restrictions apply to the availability of these data, which were used under license for this study. Data are available with the permission of the enterprise.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Jackson, J.K.; Weiss, M.A.; Schwarzenberg, A.B.; Nelson, R.M.; Sutter, K.M.; Sutherland, M.D. Global Economic Effects of COVID-19. In Congressional Research Service (CRS) Reports and Issue Briefs; Congressional Research Service: College Park, MD, USA, 2021. [Google Scholar]
  2. Andrijevic, M.; Schleussner, C.; Gidden, M.; McCollum, D.; Rogelj, J. COVID-19 Recovery Funds Dwarf Clean Energy Investment Needs. Science 2020, 370, 298–300. [Google Scholar] [CrossRef] [PubMed]
  3. Farid, S.; Naeem, M.; Paltrinieri, A.; Nepal, R. Impact of COVID-19 on the Quantile Connectedness between Energy, Metals and Agriculture Commodities. Energy Econ. 2022, 109, 105962. [Google Scholar] [CrossRef] [PubMed]
  4. Irfan, M.; Zhao, Z.; Ahmad, M.; Batool, K.; Jan, A.; Mukeshimana, M.C. Competitive Assessment of Indian Wind Power Industry: A Five Forces Model. J. Renew. Sustain. Energy 2019, 11, 063301. [Google Scholar] [CrossRef]
  5. Zhao, X.; Wen, J.; Zou, X.; Wang, Q.; Chang, C. Strategies for the Sustainable Development of China in the Post-epidemic Era. Sustain. Dev. 2023, 31, 426–438. [Google Scholar] [CrossRef]
  6. Yasir, A.; Hu, X.; Ahmad, M.; Alvarado, R.; Anser, M.K.; Isik, C.; Choo, A.; Ausaf, A.; Khan, I.A. Factors Affecting Electric Bike Adoption: Seeking an Energy-Efficient Solution for the Post-COVID Era. Front. Energy Res. 2022, 9, 817107. [Google Scholar] [CrossRef]
  7. Pal, A.; Zhang, Y. Free-Floating Bike Sharing: Solving Real-Life Large-Scale Static Rebalancing Problems. Transp. Res. Part C-Emerg. Technol. 2017, 80, 92–116. [Google Scholar] [CrossRef]
  8. Liu, X. Shared Travel Is on a Low-Carbon Road. China’s Foreign Trade 2022, 2, 62–63. [Google Scholar]
  9. Zhi, X. Overall Stable and Steady Progress—Analysis of Economic Operation of China’s Bicycle Industry in 2020. Chin. Bicycl. 2021, 2, 16–21. [Google Scholar]
  10. Zhang, G. Electric Bicycle Fire Risk and Prevention Measures. Mass Stand. 2022, 3, 119–121. [Google Scholar]
  11. Zhou, Q.; Shen, M.; Wang, X.; Chen, J. Survey Report of Electric Bicycle Riders in Pudong New Area of Shanghai. Cars Saf. 2019, 6, 12. [Google Scholar]
  12. Kou, Y. Mobile Battery Replacement Path Planning Problem in Urban Shared Electric Bicycle System. Master’s Thesis, Beijing Jiaotong University, Beijing, China, 2020. [Google Scholar]
  13. Chen, Y.; Mo, G.; Shen, J. Technical Requirements and Testing Specifications for Intelligent Charging Cabinets for Electric Bicycles. Qual. Certif. 2020, 4, 3. [Google Scholar]
  14. Xu, Y.; Lin, F.; Lin, X. Analyze the Development Status and Future Trend of Take-out O2O Industry. Theory Chin. Commer. 2017, 2, 15–16. [Google Scholar]
  15. Wang, X. The Three Kingdoms of Bike-Sharing. Chin. Entrep. 2020, 10, 4. [Google Scholar]
  16. Church, R.L.; ReVelle, C.S. Theoretical and Computational Links between the P-median, Location Set-covering, and the Maximal Covering Location Problem. Geogr. Anal. 1976, 8, 406–415. [Google Scholar] [CrossRef]
  17. Murad, A.; Faruque, F.; Naji, A.; Tiwari, A.; Qurnfulah, E.; Rahman, M.; Dewan, A. Optimizing Health Service Location in a Highly Urbanized City: Multi Criteria Decision Making and P-Median Problem Models for Public Hospitals in Jeddah City, KSA. PLoS ONE 2024, 19, e0294819. [Google Scholar] [CrossRef] [PubMed]
  18. Ramadhan, M.H.; Setiawan, D.; Sulistio, J. Comparative Study on Determining Warehouse Locations Using Set-Covering and P-Median; AIP Publishing: Long Island, NY, USA, 2024; Volume 2891. [Google Scholar]
  19. Janjić, A.; Velimirović, L.; Velimirović, J.; Vranić, P. Estimating the Optimal Number and Locations of Electric Vehicle Charging Stations: The Application of Multi-Criteria p-Median Methodology. Transp. Plan. Technol. 2021, 44, 827–842. [Google Scholar] [CrossRef]
  20. Hanif, M.; Zhang, L.; Mujtaba, N.; Li, J.; Shah, A.H.; Ullah, S. A Single Non-Obnoxious Facility Location Selection for Utility Stores Corporation Using Center of Gravity and P-Median Methods. In Proceedings of the 2021 IEEE International Conference on Industrial Engineering and Engineering Management (IEEM), Singapore, 13 December 2021; pp. 693–697. [Google Scholar]
  21. Davidov, S. Optimal Charging Infrastructure Planning Based on a Charging Convenience Buffer. Energy 2020, 192, 116655. [Google Scholar] [CrossRef]
  22. Zhou, G.; Zhu, Z.; Luo, S. Location Optimization of Electric Vehicle Charging Stations: Based on Cost Model and Genetic Algorithm. Energy 2022, 247, 123437. [Google Scholar] [CrossRef]
  23. He, F.; Yin, Y.; Zhou, J. Deploying Public Charging Stations for Electric Vehicles on Urban Road Networks. Transp. Res. Part C Emerg. Technol. 2015, 60, 227–240. [Google Scholar] [CrossRef]
  24. Pan, M.; Sun, X.; Li, X.; Chen, H.; Wang, W.; Yuan, X.; Chen, Z. Charging Pile Layout Optimization Based on Set Coverage Model of Cost Control. Mod. Electr. 2021, 38, 230–236. [Google Scholar] [CrossRef]
  25. Huang, S. Spatial Layout Optimization of Urban Public Charging Stations Based on Location Model. Master’s Thesis, Henan University, Kaifeng, China, 2023. [Google Scholar]
  26. Hakimi, S. Optimum Locations of Switching Centers in a Communication Network and Medians of Graph. Oper. Res 1964, 12, 450–459. [Google Scholar] [CrossRef]
  27. Liu, Y. Two Lower-Bounding Algorithms for the p-Center Problem in an Area. Comput. Urban Sci. 2022, 2, 5. [Google Scholar] [CrossRef]
  28. Ma, Y.; Cui, N.; Jiang, Z.-Z.; Yuan, Y.; Wang, G. Group Homophily Based Facility Location Selection in Geo-Social Networks. World Wide Web 2023, 26, 33–53. [Google Scholar] [CrossRef]
  29. Yi, L.; Xianjun, Z.; Peng, K.; Ping, Z.; Chenxu, Z.; Fangxing, C. Planning of Charging Facilities Considering Terrain Characteristics; IOP Publishing: Bristol, UK, 2021; Volume 1754, p. 012057. [Google Scholar]
  30. Zhou, S. Research on Multi-Scenario Planning of EV Charging Facilities for Typical Mountain Cities. Master’s Thesis, Southwest University, Chongqing, China, 2020. [Google Scholar]
  31. Ahmad, F.; Iqbal, A.; Ashraf, I.; Marzband, M.; Khan, I. Placement of Electric Vehicle Fast Charging Stations in Distribution Network Considering Power Loss, Land Cost, and Electric Vehicle Population. Energy Sources Part A Recovery Util. Environ. Eff. 2022, 44, 1693–1709. [Google Scholar] [CrossRef]
  32. Gao, Z.; Li, X.; Jia, D.; Wang, Y. Electric Vehicle Charging Station Location Selection Based on Existing Charging Station Adjustment. Smart Power 2022, 50, 79–84, 91. [Google Scholar]
  33. Min, D.; Jiang, K.; Liu, R.; Zhou, Z. Two-Stage Siting Capacity Model of Pure Electric Truck Charging Station. J. Chongqing Univ. Technol. (Nat. Sci.) 2023, 37, 186–195. [Google Scholar]
  34. Zhang, B.; Zhao, M.; Hu, X. Location Planning of Electric Vehicle Charging Station with Users’ Preferences and Waiting Time: Multi-Objective Bi-Level Programming Model and HNSGA-II Algorithm. Int. J. Prod. Res. 2023, 61, 1394–1423. [Google Scholar] [CrossRef]
  35. Xu, J.; Ge, L.; Yang, C.; Zhang, W. The Location and Capacity of EV Charging Station Considering the Convenience of Users and the Acceptance Capacity of Distribution Network. Power Eng. Technol. 2020, 39, 112–119, 142. [Google Scholar]
  36. Tran, C.Q.; Ngoduy, D.; Keyvan-Ekbatani, M.; Watling, D. A User Equilibrium-Based Fast-Charging Location Model Considering Heterogeneous Vehicles in Urban Networks. Transp. A Transp. Sci. 2020, 17, 439–461. [Google Scholar] [CrossRef]
  37. Zhang, Q.; Mao, Q.; Liang, P. Research on Electric Vehicle Charging Station Location Based on Improved Moth Algorithm. J. Hunan Univ. Arts Sci. (Nat. Sci. Ed.) 2021, 33, 12–19. [Google Scholar]
  38. Mulay, S.R.; Desai, J.; Kumar, S.V.; Eberhard, J.N.; Thomasova, D.; Romoli, S.; Grigorescu, M.; Kulkarni, O.P.; Popper, B.; Vielhauer, V. Towards a Network Theory of the Immune System. Ann. Dimmunologie 1974, 125C, 373. [Google Scholar]
  39. Liu, Z.; Li, Z.; Chen, W.; Zhao, Y.; Yue, H.; Wu, Z. Path Optimization of Medical Waste Transport Routes in the Emergent Public Health Event of COVID-19: A Hybrid Optimization Algorithm Based on the Immune–Ant Colony Algorithm. Int. J. Environ. Res. Public Health 2020, 17, 5831. [Google Scholar] [CrossRef] [PubMed]
  40. Yang, G.; Wei, J. Warehouse Task Scheduling Based on Immune Ant Colony Optimization Algorithm. Comb. Mach. Tools Autom. Process. Technol. 2023, 1, 179–183. [Google Scholar] [CrossRef]
  41. Xu, D. Research on EV Charging Station Location Planning Based on Multi-Objective Optimization. Master’s Thesis, Shandong Jiaotong University, Jinan, China, 2023. [Google Scholar]
  42. Wang, B.; Hu, J.; Zhang, K.; Li, X.; Zhang, Y.; Sun, S. Research and Application of Improved Immune Algorithm Based on Solution Space Oriented Optimization. Fire Command Control 2020, 45, 64–67. [Google Scholar]
  43. Xu, Y.; Zhang, J. Regional Integrated Energy Site Layout Optimization Based on Improved Artificial Immune Algorithm. Energies 2020, 13, 4381. [Google Scholar] [CrossRef]
  44. Zhou, Y.; Zhang, H.; Ma, L. An Tmproved Immune Algorithm for Solving the Hierarchical and Progressive Location Problem of Emergency Medical Facilities. J. Oper. Res. 2021, 25, 15–34. [Google Scholar] [CrossRef]
Figure 1. Flowchart of electric bicycle battery exchange cabinet siting process.
Figure 1. Flowchart of electric bicycle battery exchange cabinet siting process.
Sustainability 16 08394 g001
Figure 2. Improved immune algorithm process.
Figure 2. Improved immune algorithm process.
Sustainability 16 08394 g002
Figure 3. Example of the multi-point mutation process.
Figure 3. Example of the multi-point mutation process.
Sustainability 16 08394 g003
Figure 4. Example of the guided mutation process.
Figure 4. Example of the guided mutation process.
Sustainability 16 08394 g004
Figure 5. Example of the local search process.
Figure 5. Example of the local search process.
Sustainability 16 08394 g005
Figure 6. Distribution of the number of demand point areas.
Figure 6. Distribution of the number of demand point areas.
Sustainability 16 08394 g006
Figure 7. Kernel density analysis map.
Figure 7. Kernel density analysis map.
Sustainability 16 08394 g007
Figure 8. Battery exchange cabinet site selection plan.
Figure 8. Battery exchange cabinet site selection plan.
Sustainability 16 08394 g008
Figure 9. Iteration comparison of optimal solutions between the two algorithms.
Figure 9. Iteration comparison of optimal solutions between the two algorithms.
Sustainability 16 08394 g009
Table 1. Model variables and parameters.
Table 1. Model variables and parameters.
SymbolDescription
i Represents the identification index of demand points, i ϵ n i .
n i Set representing the collection of demand points.
j Represents the identification index of candidate points,  j ϵ n j .
n j Set representing the collection of candidate points.
I Represents the number of demand points.
J Represents the number of candidate points.
P Represents the number of exchange cabinets (antibody characteristics in immune algorithm).
d i j Distance between demand points.
R The service radius of battery exchange cabinet.
L i The farthest distance from the power change cabinet when the rider is most satisfied.
U i Indicates the maximum change cabinet distance the rider can tolerate.
P s Population size in the immune algorithm.
g Generations of population evolution in the immune algorithm.
r c Clone rate in the immune algorithm.
w 1 The weight of the objective function f 1 .
w 2 The weight of the objective function f 2 .
r m Mutation ratio in the immune algorithm.
Z j Decision variables.
Y i j Decision variables.
Table 2. Quantitative sensitivity analysis of electric bicycle battery exchange cabinets.
Table 2. Quantitative sensitivity analysis of electric bicycle battery exchange cabinets.
Number of Cabinets Comprehensive   Objective   f Rider   Satisfaction   f 1 Value Battery   Exchange   Station   Service   Capability   f 2 Value Difference   Value   of   f
7900.9486885890.9204858820.976891297
8000.9574885890.9290242740.9859529040.0088
8100.9601885890.9316440090.988733170.0027
8200.9612885890.9327113080.9898658710.0011
8300.9619885890.9333904980.990586680.0007
8400.9625885890.9339726610.9912045170.0006
Table 3. Sensitivity analysis for different parameters.
Table 3. Sensitivity analysis for different parameters.
r m = 0.1 r m = 0.2 r m = 0.3
Optimal ValueAverage ValueStandard DeviationOptimal ValueAverage ValueStandard DeviationOptimal ValueAverage ValueStandard Deviation
r c = 0.1Traditional0.95370.95000.00310.93090.92240.00480.91100.90670.0039
Improved0.96020.95780.00180.93820.93390.00250.92440.92050.0050
r c = 0.2Traditional0.95400.95010.00460.93090.92050.00610.90760.90410.0031
Improved0.96350.96140.00300.94040.93490.00330.91560.91240.0022
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

Zou, Z.; Yang, W.; Sheng, S.Y.; Yan, X. Research on the Location Selection Problem of Electric Bicycle Battery Exchange Cabinets Based on an Improved Immune Algorithm. Sustainability 2024, 16, 8394. https://doi.org/10.3390/su16198394

AMA Style

Zou Z, Yang W, Sheng SY, Yan X. Research on the Location Selection Problem of Electric Bicycle Battery Exchange Cabinets Based on an Improved Immune Algorithm. Sustainability. 2024; 16(19):8394. https://doi.org/10.3390/su16198394

Chicago/Turabian Style

Zou, Zongfeng, Weihao Yang, Shirley Ye Sheng, and Xin Yan. 2024. "Research on the Location Selection Problem of Electric Bicycle Battery Exchange Cabinets Based on an Improved Immune Algorithm" Sustainability 16, no. 19: 8394. https://doi.org/10.3390/su16198394

APA Style

Zou, Z., Yang, W., Sheng, S. Y., & Yan, X. (2024). Research on the Location Selection Problem of Electric Bicycle Battery Exchange Cabinets Based on an Improved Immune Algorithm. Sustainability, 16(19), 8394. https://doi.org/10.3390/su16198394

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