Operation Strategy for Electric Vehicle Battery Swap Station Cluster Participating in Frequency Regulation Service

Idle batteries in the battery swap stations (BSSs) of electric vehicles (EVs) can be used as regulated power sources. Considering the battery swap service and the frequency regulation (FR) service, this paper establishes a model of BSS cluster participating in the FR service and formulates a two-stage operation strategy. The day-ahead strategy arranges the battery charging plan and FR plan with the goal of the optimal operating economy on the next day. The intra-day strategy aims at maximizing the satisfaction degree of battery swap, minimizing the loss of planned revenue and ensuring the coordination of battery swap service and FR service by regulating the charging and discharging status of each battery in real-time. The simulation case shows that, under the prerequisite of gratifying the battery swap demand, the strategy improves the operating economy by making full use of idle batteries which bear a part in the FR service.


Introduction
Compared with charging mode, battery swapping mode of electric vehicles (EVs) has many conveniences [1], but the fundamental causes of unpopularity are its high investment and operating costs and low economic benefits. After the initial investment of the BSS is determined, how to improve the operating economy is a key issue that needs to be considered. Enabling battery swap stations (BSSs) clusters to participate in the frequency regulation (FR) service can make full use of idle batteries to gain revenue, thereby improving the operating economy of BSSs and promoting the popularization of battery swapping mode [2,3].
At present, improving the efficiency of the battery swap business is the main solution to the economic problems of BSSs, which mainly include battery charging optimization and battery swap scheduling optimization. In [4], charging period and charging power were taken as decision variables, which were researched to reduce network loss, peak-to-valley difference and power purchase costs. The authors of [5,6] used the battery in BSSs to solve the problem of photovoltaic power generation surplus, which could reduce the battery charging costs, while solving the consumption of new energy. The authors of [7] built an operating model of a centralized distribution model based on the closed-loop supply chain, established a battery logistics model using time-space network technology and optimized the management of battery charging and distribution.
To further enhance the economics of operation, some scholars have broadened the scope of research to the construction of the entire battery swap service and the relationship between the BSS operator and other entities. For example, [8] proposed a nano-grids-based charging and swapping system for EVs to improve the cleanliness of energy supply and the electrification of transportation, while [9] studied the joint control of BSSs and batteries in a dynamic energy pricing environment and proposed an optimization framework to reduce costs and energy consumption. The above-mentioned research papers only focused on the internal optimization of the service network of BSSs and did not fully consider and utilize the external environment of BSSs. The authors of [10] used reasonable control and scheduling to provide battery charging and swapping stations as active power stations to provide peak-shaving service for the power grid, whereas [11] aggregated the battery fast-charging station, BSS and energy storage system in the micro-grid into a whole and proposed a multi-time scale optimization operation strategy.
The FR market brings great profitability to industries containing energy storage technologies such as EVs [12]. The key issues are the aggregation of battery resources, the estimation of FR capabilities and the allocation of PR power. The authors of [13,14] separately studied the forecasting of disposable capacity and the development of aggregation modules for EVs to participate in the FR service, while [15] considered the uncertainty of battery charging demand and proposed a vehicle-to-grid control strategy that coordinated the development of the FR service and battery charging service. The authors of [16] proposed an aggregation model and a calculation method of disposable capacity for EVs to participate in the FR service in the electricity market and proposed a bidding strategy that considers risks. In [17], the authors proposed a two-stage self-dispatching model that took into account the day-ahead electricity market and the FR market to realize the day-ahead formulation of battery swap plans and bidding strategies. To improve the mileage score of BSS, [18] considered the uncertainty of electricity prices and FR signals and proposed a robust control model for EVs based on FR performance evaluation. The authors of [19] used Lyapunov optimization for the FR control strategy considering various factors, such as the travel demand of EVs, the cost of batteries and the access to energy to realize the optimization of decision-making without historical data. Considering the applicability of battery energy storage to participate in FR, some scholars have recently carried out research based on battery degradation. A refined model of the battery degradation principle is given in [20,21]. Based on this, the FR control strategy to minimize the degree of battery degradation is proposed and the accurate calculation method of the battery energy storage FR cost is given, both of which can promote the application of battery energy storage in the FR market.
Based on the above research inspiration, a BSS operation strategy for the collaborative development of the battery swap business and the FR service business is proposed. Under the premise of meeting the battery swap demand, it obtains benefits by making full use of idle batteries to participate in the FR service and provides a more cost-effective and efficient operation mode for BSSs.
The contributions of this article are listed as follows: 1.
The unified model of BSS cluster participating in the FR service is designed and established, which clearly describes the multi-time scale behavior of BSSs. The key items that determine the economic effect are systematically given, including battery degradation costs, power purchase costs, FR service income and battery swap service income.

2.
The optimization problem constructed in this paper is based on fewer assumptions and is more in line with engineering reality than existing research. Processes such as the battery swap service and FR service are accurately described by only integer variables and linear functions, which can avoid solving complex optimization problems and can obtain more detailed results, including charge and discharge power profiles, SOC, charger plug-in status and battery swap state.

3.
The index that characterizes the busyness of the battery swap service of BSSs is defined and, in the FR service, the FR power is optimally allocated based on it, which realizes the power support between the BSSs and minimizes the impact of the FR service on the battery swap service.
This paper is organized as follows: Section 2 introduces the model of BSS cluster participating in FR services and analyzes its operating economics. Section 3 proposes the day-ahead and intra-day two-stage strategy. Section 4 presents the case studies and analysis to verify the correctness and effectiveness of the proposed model and strategy. Section 5 concludes this paper.

Model for Battery Swap Station Cluster Participating in Frequency Regulation Service
The model for a BSS cluster participating in the FR service is designed in this section and the corresponding model is established.

Model for Battery Swap Station Cluster Participating in Frequency Regulation Service
The EV battery swapping modes include centralized charging mode and batteryswapping mode [22]. The charging and battery-swapping modes are the research object in this paper. BSS operators provide battery leasing service and battery swap service for EVs and use idle batteries to participate in the FR service to increase revenue [23]. The business model of the BSS cluster participating in the FR service is shown in Figure 1. The BSS operators are located in the middle of the industrial chain, providing a battery swap service for EVs and the FR service for the power system. EV users pay battery rental fees and battery swap service fees through BSSs to obtain batteries swap service; the power sector organizes and supervises the FR service market [24].
Processes 2021, 9, x FOR PEER REVIEW 3 of 16 3. The index that characterizes the busyness of the battery swap service of BSSs is defined and, in the FR service, the FR power is optimally allocated based on it, which realizes the power support between the BSSs and minimizes the impact of the FR service on the battery swap service.
This paper is organized as follows: Section 2 introduces the model of BSS cluster participating in FR services and analyzes its operating economics. Section 3 proposes the dayahead and intra-day two-stage strategy. Section 4 presents the case studies and analysis to verify the correctness and effectiveness of the proposed model and strategy. Section 5 concludes this paper.

Model for Battery Swap Station Cluster Participating in Frequency Regulation Service
The model for a BSS cluster participating in the FR service is designed in this section and the corresponding model is established.

Model for Battery Swap Station Cluster Participating in Frequency Regulation Service
The EV battery swapping modes include centralized charging mode and batteryswapping mode [22]. The charging and battery-swapping modes are the research object in this paper. BSS operators provide battery leasing service and battery swap service for EVs and use idle batteries to participate in the FR service to increase revenue [23]. The business model of the BSS cluster participating in the FR service is shown in Figure 1. The BSS operators are located in the middle of the industrial chain, providing a battery swap service for EVs and the FR service for the power system. EV users pay battery rental fees and battery swap service fees through BSSs to obtain batteries swap service; the power sector organizes and supervises the FR service market [24].

Model for Battery Swap Station Cluster Participating in Frequency Regulation Service
Under the model of a BSS cluster participating in the FR service, the costs borne and benefits obtained by BSS operators are analyzed in this section and the model of a BSS cluster participating in the FR service is established.

Model for Battery Swap Station Cluster Participating in Frequency Regulation Service
Under the model of a BSS cluster participating in the FR service, the costs borne and benefits obtained by BSS operators are analyzed in this section and the model of a BSS cluster participating in the FR service is established.

Operating Costs
The investment costs [25] and annual financial costs [26] of BSSs can be determined in the initial construction stage. The operating costs mainly consist of operation and maintenance costs, power purchase costs and degradation costs.
where Z day total is the daily operating costs (USD), Z 1 is the daily maintenance costs (USD), α i is the industry operation and maintenance rate of the ith equipment, P i is the investment costs of the ith equipment (USD), ε is the average daily wage of the employee (USD/person), n is the number of employees (person), Z 2 is the daily power purchase costs (USD), E i represents the charging demand of the ith BSS on that day (kWh), η i,ch is the charging efficiency of the ith BSS, η i,j is the power utilization efficiency of the ith BSS, P grid is the power purchase price of the BSS (USD/kWh), Z 3 is the daily degradation cost (USD), π ess is the investment cost per unit capacity of the battery (USD/kWh), C cap is the battery capacity (kWh), T f loat is the floating charge life of the battery (day) [27] and C cycle ess is the average daily cycle costs of the battery (USD) [28]. (2) where N eq,t is the equivalent number of full cycles of energy storage in the tth period, P ch,i (t)/P dis,i (t) is the actual value of charge/discharge power of the ith battery in the tth period [29], η ch,i /η dis,i is the charge/discharge efficiency of the ith battery and k p is the fitted energy storage characteristic parameter, generally 0.8-2.1.

Operating Income
Under the mode of BSS cluster participating in the FR service, the operating income of the BSS operator mainly comes from the battery swap business and the FR service business.
where S day total is the daily operating income (USD), S 1 is the income from daily battery swap service (USD), Q max,i is the maximum battery capacity of the ith EV (kWh) [30], SOC m.i is the state of charge (SOC) requirement of the ith EV, SOC i is the SOC of the ith EV before swapping the battery [31], P sell is the price of the battery swap service (USD/kWh), µ is the fixed charge for the battery swap service (USD), S 2 is the revenue from the daily FR service (USD), A j is the signal mileage of the BSS operator in the trading cycle of the jth FR service (kW), k j is the average of the composite FR score [32] for the jth trading cycle, P grid,j is the mileage settlement price of the jth trading cycle (USD/kW), Q j is the FR capacity awarded in the jth trading cycle (kWh), P j is the FR capacity compensation pricing (USD/kWh), S 3 is the daily residual value return (USD), C i is the charging and discharging times of the ith battery (times), C i,M is the maximum charge-discharge times of the ith battery (times) and σ is the residual value yield, that is, the proportion of the residual value income when the fixed asset is scrapped to its original value.

Model for Battery Swap Station Cluster Participating in Frequency Regulation Service
Therefore, a day-ahead and intra-day dual time scale model is established to formulate an FR plan and power allocation plan meeting the demand for power exchange. The details are as follows.
where A da n×m is the FR capacity distribution coefficient matrix of n BSSs in m periods of the next day, which is the output matrix of the model in the day-ahead stage, Q grid 1×m is the tendered FR capacity (kWh), which is the input matrix of the model in the day-ahead stage, Q da 1×n represents the maximum available capacity of each BSS (kWh),A rt n×1 is the power distribution coefficient matrix of a single FR signal for n BSSs in the intra-day stage, which is the output matrix of the model in the intra-day stage, P grid is the single FR signal in the intra-day stage (kW), which is the input matrix of the intra-day stage model, and P rt 1×n is the maximum adjustable power of BSSs (kW).
Given that uncertainty [33] affects operations, some studies adopt the robust optimization framework for system scheduling and control [34], which can lead to stable and reliable solutions, but the solutions are often conservative, that is, the system capabilities cannot be fully utilized. So, in this paper, some parameters in the model are treated as random variables, which mainly include the battery swap demand that obeys the normal distribution and the time when the battery swap demand arises that obeys the uniform distribution [35]. Further, the impact of the uncertainty on the system is reduced through the day-ahead plan and intra-day adjustments.

Two-Stage Strategy for Battery Swap Station Cluster Participating in Frequency Regulation Service
In order to ensure that the battery swap service and FR service are carried out costeffectively, this chapter proposes an optimal operation model which includes two parts, day-ahead and intra-day operation.
In the day-ahead stage, an operation plan is made with the goal of maximizing the net income of the next day; it includes a battery charging plan and an FR plan. In the intra-day stage, when the BSS operator receives the FR signal, the FR power is allocated with the goal of maximizing the satisfaction degree of battery swap service and minimizing the loss of planned income, so the charging and discharging status of each BSS is regulated. The detailed optimal model described in this paper is as follows.

Day-Ahead Operation
A BSS cluster participating in the FR service needs to coordinate the operation plan of each station in the day-ahead stage. It specifically includes the following:

1.
Determining the available FR capacity of each station on the next day; 2.
Arranging a battery charging plan. In order to achieve the above goal, it is necessary to optimize the solution based on the predicted value of battery swap demand and the FR demand issued by the power sector.
Therefore, for each battery in each period, its charging power (real), discharging power (real), charge-discharge state variable (integer), charger plug-in state variables (integer), state variable indicating whether the battery is fully charged (integer) and state variable of battery swap (integer) are used as decision variables to construct an optimization model.

Objective Function
In the model for a BSS cluster participating in the FR service, batteries are used for the FR service, so the average number of daily cycles increases, resulting in more degradation costs and potential power purchase costs, but obtains the FR service income and some Processes 2021, 9, 1513 6 of 15 residual value income. The optimization goal of the day-ahead stage is to maximize the net income of the next day, as described in the following formula: where RS is the daily net income of the BSS operator (USD).

Constraints
When optimizing the operation plan for the next day, it is necessary to comprehensively consider the constraints of battery energy constraints, FR service rules and battery swap demand.

1.
Battery energy constraints The charge and discharge power of the battery at any time is less than its maximum limit.
where P ch,i,max /P dis,i,max is the maximum limit of the charge/discharge power of the ith battery and I i (t) is the charge-discharge state variable, which takes 1 to indicate that the ith battery is charged, otherwise it is discharging.
The SOC of the battery remains within the allowable range at any time to avoid a sharp drop in life and the battery SOC has a quantitative relationship with its charges and discharges power.
where SOC i,min /SOC i,max is the lower/upper limit of the i th battery SOC, SOC i (t) is the SOC of the ith battery in the tth period and E i,0 is the ith battery rated capacity (kWh).
To ensure that the optimization strategy is continuously executed, the electricity stored in each station needs to be equal at the beginning and end of each day. The battery swap process is treated as a quick discharge that rapidly reduces the SOC of the battery to the level of the inbound battery in this paper.
where SOC i, f irst is the SOC at the beginning of the next day of the ith battery, SOC i, f inal is the SOC at the end of the next day of the ith battery and N is the collection of batteries in the Nth power station.

2.
FR service rules In this paper, we assume that, in the FR service, the BSS operator declares the equivalent upward and downward FR capacity and that the FR signals issued by the power dispatch agency in each transaction cycle are electrically neutral.
where Q t is the tendered FR capacity for the tth period (kWh), Q i,max is the total battery power for the ith power change in the tth period (kWh), Q i,t,ev is the power needed for the battery swap service of the ith BSS in the tth period (kWh), Q i,t,ch is the remaining charge of the battery that needs to be continuously charged during the tth period (kWh), Q i,rate is the total capacity of all batteries for the ith BSS (kWh), Q i,t,ch is the amount of charge that can be obtained by the battery that needs to be continuously charged in the ith BSS in the tth period (kWh) and Q grid,t is the maximum amount of FR service that an FR service provider can tender during the tth period (kWh).

Battery swap demand
Carrying out the battery swap business requires that the number of available batteries in each station at each period is greater than the battery swap demand.
where Q j is the rated capacity of the jth battery (kWh), C j,t is the state variable indicating whether the battery is fully charged, which takes 1 to indicate that the jth battery is full at the end of the tth period, otherwise, it is not fully charged, D j,t is the state variable of battery swap, taking 1 means that the jth battery was swapped out of the station in the tth period, otherwise, it was not swapped out, Q h,i is the difference in battery power before and after the ith battery swap (kWh) and dem n,k+1 is the predicted value of the battery swap demand of the nth BSS in the k + 1th period (pc).

Intra-Day Operation
The uncertainty of the battery swap demand and the FR demand make the operation plan unable to be fully implemented in the intra-day stage. To ensure the stable operation of BSSs in the intra-day stage, it is necessary to regulate the charging and discharging status of each power station in real-time to ensure the reasonable distribution of FR power.
Therefore, for each battery, when the FR signal is received, its charging power (real), discharging power (real) and response degree of the FR signal (real) are used as decision variables to construct an optimization model.

Objective Function
The goal of power allocation decision-making on the FR signal based on real-time data is to ensure the coordination of the battery swap service and FR service to maximize battery swap satisfaction and minimize planned revenue loss. Its objective function is described in the following formula: where Z 3k is the degradation cost corresponding to the kth FR (USD), S 2k is the mileage revenue for the kth FR (USD), S 3k is the residual value income for the kth FR (USD), B n is the battery swap saturation indicator for the nth BSS, P n,k is the FR power provided by the nth station in k times of FR (kW), positive means downward FR, and ∆t is the FR signal duration (s).
where x n,t,k is the number of battery swap demands that occurred in the tth period when the kth FR signal was issued.

Constraints
When the intra-day operating strategy is formulated, it is necessary to comprehensively consider the constraint conditions of the loss of revenue and the power balance of the grid and battery.
(1) Loss of revenue where v k is the response degree of FR signal, the value range is [0,1], A k is the FR signal value (kW), C st,i,k is the total number of charge and discharge times of the batteries in the ith BSS in the kth FR (times), C st,i,M is the sum of the maximum charge and discharge times of each battery in the ith BSS (times) and P st,i is the sum of the prices of each battery in the ith BSS (USD).
(2) Power constraints of the grid and battery ∑ ( where P D,max is the maximum transmission power limit for the point of common coupling (PCC) D (kW) [36], A st,n,k is the FR power of the nth BSS in kth FR (kW) and Q t,n is the available FR capacity of the nth BSS in the tth period (kWh). In addition, the optimization problem should meet Equations (6)-(8).

Two-Stage Strategy for Battery Swap Station Cluster Participating in Frequency Regulation Service
The solution sequence of the above optimization model is to solve the day-ahead stage first, then solve the intra-day stage. The specific solution process is shown in Figure 2. The scale of the optimization problem at each stage is shown in Table 1  In addition, the optimization problem should meet Equations (6)-(8).

Two-Stage Strategy for Battery Swap Station Cluster Participating in Frequency Regulation Service
The solution sequence of the above optimization model is to solve the day-ahead stage first, then solve the intra-day stage. The specific solution process is shown in Figure  2. The scale of the optimization problem at each stage is shown in Table 1 below.

Begin
Initialization data: configuration parameters of the power swap stations, frequency regulation market information, forecast value of battery swap demand.
Form the objective and constraints of the day-ahead stage.
Iteratively reduce the feasible solution space, find the optimal solution, and obtain the optimal charge-discharge curve and frequency regulation capacity of the BBS.
Are all groups solved? N Output data: available frequency regulation capacity of each BSS in each period, batteries charging plan.

Y
Enter the intra-day stage. Input data: batteries charging plan, frequency regulation capacity for each period, real-time frequency regulation signal. Form the objective and constraints of the intra-day stage.
Iteratively reduce the feasible solution space, optimizing the real-time frequency regulation power allocation scheme.
Output data: optimal frequency regulation power allocation scheme   Figure 2. Two-stage strategy solution process for BSS cluster participating in FR service. Among them, the intra-day stage is a linear programming model, which is easy to solve by using classic methods. The day-ahead stage involves multiple integer variables, such as charge and discharge state variables, full charge state variables and battery swap state variables. They are coupled in the battery swap service constraint. At the same time, the total number of batteries in the service area is usually large, so it is easy to form a more complex large-scale mixed-integer linear programming (MILP) problem [37]. In this paper, the basic method to solve the problem is the branch and bound method according to the grouping of BSSs. Its specific steps are as follows:

1.
Divide the original problem into n MILP problems according to the number of BSSs; 2.
Perform linear relaxation on the MILP problem and determine the relaxed solution space Ω 0 and the corresponding upper and lower bounds of the objective function; 3.
The substitution problem, after linear relaxation, is divided into several sub-problems W i , whose solution set is X i and requires Ω 0 ∈ X 1 ∪ X 2 ∪ · · · · · · X i . For each subproblem, if the optimal solution of the sub-problem is a feasible solution of the original problem, it is the optimal solution of the MILP problem and the calculation is completed; otherwise, the value of the objective function is regarded as the new upper bound U 1 of the MILP problem. The optimal solution of the sub-problem that is the feasible solution of the MILP problem is selected and its objective function is regarded as the lower bound L 1 of the MILP problem; 4.
Abandon the sub-problems where the objective function value of the optimal solution is less than L 1 and keep the sub-problems where the objective function value of the optimal solution is greater than L 1 ; 5.
Select the sub-problem with the largest objective function of the optimal solution and repeat 1 and 2. If the optimal feasible solution of the sub-problem is found, the maximum value of the objective function of the feasible solution and all the previously retained sub-problems is regarded as the new lower bound L 2 and d is repeated until the optimal solution is found; 6.
Judge whether to traverse all the BSSs; if yes, end, otherwise, select the next BSS and go back to 2.

Basic Data
In order to verify the effectiveness of the strategy, six BSSs established by a BSS operator are simulated and analyzed. The relevant parameter settings are shown in Table 2. Each BSS is located in different locations in the service area, so the battery swap demand of each station is different. The battery swap demand in 24 periods of a typical day is selected for simulation and the relevant data are shown in Appendix A. It is assumed that the SOC of the battery obeys the probability distribution of N (0.2, 0.042) when the user of the EV generates the battery swap demand and the generation of the battery swap demand in each period obeys the uniform distribution. There are two peaks in the battery swap demand curve, which are 9-11 h and 15-20 h. The battery swap demand at night is little, showing a low state.  Table 3 and Figure 3 are the simulation results of the day-ahead operation strategy for a BSS cluster participating in FR service [38]. It can be seen, from Table 3, that, under the strategy proposed in this paper, each BSS obtains some revenue from the FR service while obtaining the revenue from the battery swap service. Among them, BSS 3 has the highest daily battery swap income, because the battery swap demand of this station is the highest level on that day; BSS 5 has the least demand for batteries swap and the lowest daily battery swap income, but this station has more idle batteries for the FR service and its daily FR yield is the highest, which makes up for the overall profit of this station on that day.   It can be seen, from Figure 3, that the FR capacity provided by the BSS operator varies in 24 periods. Figure 4 shows the comparison between the FR capacity and the total battery swap demand. It illustrates that the available FR capacity curve of each BSS and the battery swap demand curve show an opposite trend. In general, the available FR capacity It can be seen, from Figure 3, that the FR capacity provided by the BSS operator varies in 24 periods. Figure 4 shows the comparison between the FR capacity and the total battery swap demand. It illustrates that the available FR capacity curve of each BSS and the battery swap demand curve show an opposite trend. In general, the available FR capacity is low when the demand for the battery swap service is high; the available FR capacity is high when the battery swap demand is low. Moreover, the change of available FR capacity is ahead of the total battery swap demand; for example, the sharp decrease in the available FR capacity in the 5-6 periods corresponds to the surge in the total battery swap demand in the 8-9 periods, while the sharp decrease in the available FR capacity during the 13-16 periods corresponds to the surge in the total power exchange demand during the 15-18 periods. It can be seen that, under the influence of the strategy in this paper, each BSS arranges battery charging before the period of high battery swap demand to ensure the sustainability of the battery swap business, so the available FR capacity is reduced at this time. It can be seen, from Figure 3, that the FR capacity provided by the BSS operator varies in 24 periods. Figure 4 shows the comparison between the FR capacity and the total battery swap demand. It illustrates that the available FR capacity curve of each BSS and the battery swap demand curve show an opposite trend. In general, the available FR capacity is low when the demand for the battery swap service is high; the available FR capacity is high when the battery swap demand is low. Moreover, the change of available FR capacity is ahead of the total battery swap demand; for example, the sharp decrease in the available FR capacity in the 5-6 periods corresponds to the surge in the total battery swap demand in the 8-9 periods, while the sharp decrease in the available FR capacity during the 13-16 periods corresponds to the surge in the total power exchange demand during the 15-18 periods. It can be seen that, under the influence of the strategy in this paper, each BSS arranges battery charging before the period of high battery swap demand to ensure the sustainability of the battery swap business, so the available FR capacity is reduced at this time.

Results and Comparisons of Strategy in the Intra-Day Stage
Based on the simulation design in Section 4.1, the intra-day operation strategy simulation is carried out and the strategy proposed in this paper is compared with the traditional distribution based on capacity strategy.
After the BSS operator receives the FR signal, it allocates the FR power intending to take into account the coordination of the FR service and battery swap service and realizes the response of the FR signal. The total output power of 6 stations in 24 periods is shown in Figure 5.

Results and Comparisons of Strategy in the Intra-Day Stage
Based on the simulation design in Section 4.1, the intra-day operation strategy simulation is carried out and the strategy proposed in this paper is compared with the traditional distribution based on capacity strategy.
After the BSS operator receives the FR signal, it allocates the FR power intending to take into account the coordination of the FR service and battery swap service and realizes the response of the FR signal. The total output power of 6 stations in 24 periods is shown in Figure 5. Among them, the FR power in the 8th period is shown in Figure 6 and the distribution of the battery swap demand is shown in Table 4. The red circle representing the battery swap demand generated at this time makes the actual battery swap demand in this Among them, the FR power in the 8th period is shown in Figure 6 and the distribution of the battery swap demand is shown in Table 4. The red circle representing the battery swap demand generated at this time makes the actual battery swap demand in this period greater than the predicted value. From Table 4, the battery swap saturation level of each station can be calculated when the FR signal is issued, which affects the FR power distribution. For example, the BSS cluster responds to the upward FR signal at 7:15, that is, the overall system needs to reduce the charging power. BSS 3 and BSS 6 have the highest battery swap saturation; they provide the largest FR power and reduce the minimum charging power, so that they have more sufficient power to meet possible battery swap demand. The BSS cluster responds to the downward FR signal at 7:50, that is, the overall system needs to increase the charging power. At this time, the battery swap saturation is the lowest in BSS 4 and 6, so they provide the smallest FR power, that is, their increased charging power is the smallest, which, in turn, makes the BSS with a higher degree of battery swap saturation have a higher charging power to cope with the possible battery swap demand.    Figure 7 shows the comparison of the FR power of BSS 3 in the 8th period under the strategy of this paper and the strategy of capacity allocation. It can be seen that the FR power of BSS 3 under the strategy of this paper is basically positive and the FR power under the capacity allocation strategy is determined by the positive or negative value of the FR power of the system. According to Table 4, it can be seen that the battery swap demand predicted value for this period is reached at 7:16 and the battery swap demand predicted value is exceeded at 7:50. Compared with the capacity allocation strategy, the strategy in this paper can make BSS 3 have more sufficient power at 7:50 to meet the exceeded battery swap demand.     Figure 7 shows the comparison of the FR power of BSS 3 in the 8th period under the strategy of this paper and the strategy of capacity allocation. It can be seen that the FR power of BSS 3 under the strategy of this paper is basically positive and the FR power under the capacity allocation strategy is determined by the positive or negative value of the FR power of the system. According to Table 4, it can be seen that the battery swap demand predicted value for this period is reached at 7:16 and the battery swap demand predicted value is exceeded at 7:50. Compared with the capacity allocation strategy, the strategy in this paper can make BSS 3 have more sufficient power at 7:50 to meet the exceeded battery swap demand.   Figure 7 shows the comparison of the FR power of BSS 3 in the 8th period under the strategy of this paper and the strategy of capacity allocation. It can be seen that the FR power of BSS 3 under the strategy of this paper is basically positive and the FR power under the capacity allocation strategy is determined by the positive or negative value of the FR power of the system. According to Table 4, it can be seen that the battery swap demand predicted value for this period is reached at 7:16 and the battery swap demand predicted value is exceeded at 7:50. Compared with the capacity allocation strategy, the strategy in this paper can make BSS 3 have more sufficient power at 7:50 to meet the exceeded battery swap demand.   Figure 7 shows the comparison of the FR power of BSS 3 in the 8th period under the strategy of this paper and the strategy of capacity allocation. It can be seen that the FR power of BSS 3 under the strategy of this paper is basically positive and the FR power under the capacity allocation strategy is determined by the positive or negative value of the FR power of the system. According to Table 4, it can be seen that the battery swap demand predicted value for this period is reached at 7:16 and the battery swap demand predicted value is exceeded at 7:50. Compared with the capacity allocation strategy, the strategy in this paper can make BSS 3 have more sufficient power at 7:50 to meet the exceeded battery swap demand.   Figure 7 shows the comparison of the FR power of BSS 3 in the 8th period under the strategy of this paper and the strategy of capacity allocation. It can be seen that the FR power of BSS 3 under the strategy of this paper is basically positive and the FR power under the capacity allocation strategy is determined by the positive or negative value of the FR power of the system. According to Table 4, it can be seen that the battery swap demand predicted value for this period is reached at 7:16 and the battery swap demand predicted value is exceeded at 7:50. Compared with the capacity allocation strategy, the strategy in this paper can make BSS 3 have more sufficient power at 7:50 to meet the exceeded battery swap demand.    Figure 7 shows the comparison of the FR power of BSS 3 in the 8th period under the strategy of this paper and the strategy of capacity allocation. It can be seen that the FR power of BSS 3 under the strategy of this paper is basically positive and the FR power under the capacity allocation strategy is determined by the positive or negative value of the FR power of the system. According to Table 4, it can be seen that the battery swap demand predicted value for this period is reached at 7:16 and the battery swap demand predicted value is exceeded at 7:50. Compared with the capacity allocation strategy, the strategy in this paper can make BSS 3 have more sufficient power at 7:50 to meet the exceeded battery swap demand.    Figure 7 shows the comparison of the FR power of BSS 3 in the 8th period under the strategy of this paper and the strategy of capacity allocation. It can be seen that the FR power of BSS 3 under the strategy of this paper is basically positive and the FR power under the capacity allocation strategy is determined by the positive or negative value of the FR power of the system. According to Table 4, it can be seen that the battery swap demand predicted value for this period is reached at 7:16 and the battery swap demand predicted value is exceeded at 7:50. Compared with the capacity allocation strategy, the strategy in this paper can make BSS 3 have more sufficient power at 7:50 to meet the exceeded battery swap demand. 500 1000 1500 tion power Figure 7 shows the comparison of the FR power of BSS 3 in the 8th period under the strategy of this paper and the strategy of capacity allocation. It can be seen that the FR power of BSS 3 under the strategy of this paper is basically positive and the FR power under the capacity allocation strategy is determined by the positive or negative value of the FR power of the system. According to Table 4, it can be seen that the battery swap demand predicted value for this period is reached at 7:16 and the battery swap demand predicted value is exceeded at 7:50. Compared with the capacity allocation strategy, the strategy in this paper can make BSS 3 have more sufficient power at 7:50 to meet the exceeded battery swap demand. under the capacity allocation strategy is determined by the positive or negative value of the FR power of the system. According to Table 4, it can be seen that the battery swap demand predicted value for this period is reached at 7:16 and the battery swap demand predicted value is exceeded at 7:50. Compared with the capacity allocation strategy, the strategy in this paper can make BSS 3 have more sufficient power at 7:50 to meet the exceeded battery swap demand.  The case of the intra-day operation strategy was performed by multiple optimization calculations using CPLEX/MATLAB 2016b. The longest time for a single optimization is 1.0431 s (problem-construction time + problem-solving time); the time to solve the problem is within 60 ms, which can give enough time for other links, such as data communication.

Conclusions
To increase the operating income of BSSs, a day-ahead and intra-day two-stage strategy for a BSS cluster participating in the FR service is proposed. The following conclusions are drawn through simulation analysis:

1.
This work has a two-fold contribution. In theory, it provides a systematic and achievable method for a BSS cluster to participate in the FR service. In practice, it makes full use of idle batteries to participate in the FR service, which can improve the operating economy of BSSs. In addition, different from the conservativeness of other methods, such as robust optimization [39], the method in this paper is more conducive to improving the economy, while still maintaining a certain level of robustness through the power support and the planned use of limited resources.

2.
From the results of optimized operations, it can be seen that the battery swap income and FR income are complementary. When the battery exchange income is high/low, FR income is low/high. By using limited resources in a planned way, the FR service can bring great economic benefits to BSS operators, especially when the battery swap demand is low. 3.
FR power is allocated based on the busyness of the battery swap service, which realizes the power support between the BSSs. Numerical experiments using a realistic battery swapping project date and including comparison with the traditional method show that this method is more economical and synergistic, while still maintaining fast calculation speed, which is conducive to real-time regulation. Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy restrictions.

Conflicts of Interest:
The authors declare no conflict of interest.
Appendix A Table A1. Battery swap demand on a typical day.