A Sustainability Improvement Strategy of Interconnected Data Centers Based on Dispatching Potential of Electric Vehicle Charging Stations

: With the rapid development of information technology, the electricity consumption of Internet Data Centers (IDCs) increases drastically, resulting in considerable carbon emissions that need to be reduced urgently. In addition to the introduction of Renewable Energy Sources (RES), the joint use of the spatial migration capacity of IDC workload and the temporal ﬂexibility of the demand of Electric Vehicle Charging Stations (EVCSs) provides an important means to change the carbon footprint of the IDC. In this paper, a sustainability improvement strategy for the IDC carbon emission reduction was developed by coordinating the spatial-temporal dispatch ﬂexibilities of the IDC workload and the EVCS demand. Based on the Minkowski sum algorithm, a generalized ﬂexible load model of the EVCSs, considering trafﬁc ﬂow and Road Impedance (RI) was formulated. The case studies show that the proposed method can effectively increase the renewable energy consumption, reduce the overall carbon emissions of multi-IDCs, reduce the energy cost of the DCO, and utilize the EV dispatching potential. Discussions are also provided on the relationship between workload processing time delay and the renewable energy consumption rate.


Background
With the rapid development of information technologies such as big data, cloud computing, and so on, internet data centers (IDCs) have been widely deployed around the world. Rapidly growing data load causes the explosive growth of electricity load, resulting in an increasingly severe carbon pollution situation [1]. Currently, China's data center industry uses about 117.181 billion kWh of thermal power throughout the year, resulting in 98.55 million tons of carbon emissions. It is reported that by 2025, carbon emissions from the IDCs will account for 3.2% of global carbon emissions. In 2040, this figure will increase to 14% [2]. This growing trend runs counter to global energy sustainability goals. Europe's strategy is to drive data centers to be carbon neutral by 2030 [3]. Events like the COVID-19 pandemic also place a special burden on information and communication technology (ICT) [4]. The sustainability of the energy and environment of the IDCs becomes a priority for the ICT industry. Reducing carbon emissions is the main goal and social responsibility of the data center, which is one of the most energy-intensive industries [5]. However, the implementation of a scientific strategy must take into account the economic costs of its wide-scale implementation and the required infrastructure [6]. In addition to the goal of reducing carbon emissions, one of the important driving forces for data center operators to purchase renewable energy is that the price of renewable energy power generation is lower than the current grid price due to the policy support and low-cost generation (thermal power generation accounted for 73.36% of China's public power grid in March 2022 [7], and the feed-in tariff of thermal power generation is higher than that of solar power generation and wind power generation). At present, low-carbon energy sources are favored by cloud service operators in the ICT industry, such as wind and solar energy. These renewable energy sources (RES) are introduced into the IDCs to supply power, thereby alleviating their carbon pollution problems. A high proportion of RES generation has replaced traditional thermal power in a large proportion. However, due to the intermittent output and the severe lack of adjustment capacity of RES, the IDCs cannot reduce their heavy dependence on traditional energy. Relying solely on purchasing renewable energy is not the most effective solution to the IDC sustainability issues.
To solve this problem, it is necessary to make full use of the flexible dispatching capacity of data load among multiple interconnected IDCs to realize the transfer of electricity load on the spatial scale, which can improve the flexibility of the IDC power consumption and reduce the negative impact of the intermittent RES. The IDC workload has considerable dispatching potential. According to a research report on the IDC energy efficiency released by the National Resource Defense Council, the average utilization rate of servers in a typical IDC is only 12% to 18% of its computing power [8]. Therefore, making full use of the flexible dispatching capability of the IDCs for their workloads is one of the important ways to improve the sustainable development of interconnected IDCs. The IDC workload can be divided into two categories: interactive workload and batch workload. The interactive workload is not time-shiftable. At present, most of the IDC workload dispatching systems under the cloud computing framework will back up the data according to the importance of Interactive workload in the offsite servers in advance. This process is called "redundant backup" [9]. In non-emergency situations, this technique can be used to alter the spatial distribution of workloads. Therefore, it is possible to realize the remote processing of online loads in a very short time by changing the calculation dispatching instructions, which provides the feasibility for the realization of spatial dispatching of workload.
Based on these methods, it is necessary to fully exploit the potential of local adjustable resources in the IDCs to further realize the efficient use of RES and fundamentally reduce the carbon emission pollution caused by traditional energy. Under the vision of sustainable global energy development, the method of replacing gasoline vehicles with electric vehicles (EVs) for decarbonization has been emphasized. The phasing out of sales of new conventional gasoline vehicles by a specified date is one of a series of regulatory policies announced over the past few years. The UK Department for Transport's report (the "Road to Zero" strategy, or R2Z) proposes an ambition for ultra-low emission vehicles (ULEVs) sales of 50-70% by 2030, ahead of a ban on sales of diesel and petrol cars by 2040 [10]. The development of plug-in electric vehicle (PEV) technology and quantity provide the opportunity to become a flexible and adjustable resource [11]. Plenty of EVs forms a key element of global sustainability strategies.

Previous Work
At present, many scholars have studied the application of data load spatial dispatching in pieces of literature. Hu. C proposed a Dynamic Time Scale-based server Provision (DTSP) method, considering the variability of workloads when providing servers for workload demands [12]. Yu L. et al. proposed a cross-domain distributed IDC and power system collaborative operation framework to achieve the optimal cost of energy and energy storage investment [13]. Li J. et al. defined the demand responsiveness of the IDC to participate in the electricity market, minimizing electricity costs through dynamic IDC server consolidation and dispatching [14]. Li Y. et al. integrated the utility grid with the IDCs through dynamic pricing, directing workload to areas where green energy is more readily available, reducing the IDC costs [15]. Renugadevi  and carbon emissions [16]. Yang. T proposed a spatio-temporal task migration mechanism to pursue low carbon in dual-dimension without considering the regulation ability of EVs in the future [17].
Some studies utilize UPS backup battery bank response to optimize the IDC performance and improve RES utilization [18,19]. In addition to power components such as wind turbines, solar panels, etc., DATAZERO proposed by PIERSON J M et al. conceives a new IDC electric-gas structure considering supercapacitors and H2 fuel cells [20]. Yu L. et al. integrated the IDC and RES power generation into a microgrid for coordinated operation, and the transfer of data load can reduce the uncertainty of RES power generation [21].
EVs have been widely studied as flexible and adjustable resources to improve the consumption rate of RES. Dai Q. et al. established an integrated system including EV charging stations, small photovoltaic systems, and energy storage systems, effectively improving the utilization of solar energy resources and sustainable urban efficiency [22]. Longo M. et al. showed that the significant introduction of RES and EVs in Canada and Italy can provide huge opportunities for energy production capacity and demand [23]. Yue J. et al. established a microgrid aggregator model including RES and EVs, coordinating the economic desire of microgrid (MG) owners and the stability operation demand of the distribution system operator (DSO) through the Stackelberg game method [24]. In reference [25], Yu L. et al. proposed a total cost minimization problem of the IDC operators by jointly managing the IDC load and EV charging, using ADMM distributed algorithm to achieve joint energy management, reducing the demand charge of the IDC. However, the scenario considered in this study is the charging of EVs of the IDC employees where the IDC can collect the information of each EV and make a dispatching plan. Considering the issue of user privacy protection, this scenario does not apply to the EVCSs open to the public. With the increasing number of EVs and the improvement of load aggregation technology and market, a wider range of EV flexibility resources can be used to reduce the carbon emissions of the IDCs and improve their sustainability. In addition, Li M. analyzed the cyber security risks that may have existed when EVs were connected to the cyber-physical power system [26].
From the previous work, the strategies for joint dispatching of the IDC workload mostly focus on reducing costs, and the impact of their carbon emissions is rarely considered. In addition, the potential of the EVCS has not been fully stimulated, and the research on its use in the IDC energy dispatching needs to be improved. In this study, the carbon emission reduction problem of the spatial-temporal coupling of the IDCs is solved considering the workload space adjustment potential of the IDC and the dispatching potential model of the EVCS.

Contributions
This work aims to minimize the DCO carbon emissions. Energy-intensive workloads are transferred to locations with sufficient RES for processing, so that excess carbon emissions caused by high-intensity loads are shifted, and these carbon emissions are offset by RES. The data center operator (DCO) can purchase the flexibility resources of an entire EVCS from a load aggregator, so that the DCO can fully utilize the load temporal adjustment capability of the EVCS, for better matching the RES output curve. The main contributions of this paper are presented, as follows: • A spatial-temporal coupled dispatching strategy aimed at DCO carbon emission reduction is proposed. The spatial dispatching ability of workload in interconnected IDCs and temporal dispatching ability of the EVCSs as generalized flexible load are considered. • A generalized flexible load dispatching potential evaluation method for the EVCSs and their surrounding road networks based on the Minkowski sum algorithm is proposed. Combined with the Road Impedance (RI) model, the dispatching potential assessment of the EVCS is more realistic. • Based on the proposed strategy, the impact of the time delay of processing workload in the IDC on the carbon emission and renewable energy consumption capacity is analyzed.

Organization of This Article
The remainder of this paper is organized, as follows. Section 2 introduces the dispatching strategy of the system and DCO carbon emission model. The IDC power consumption model and the EVCS dispatching potential model considering traffic flow are presented in Section 3. Subsequently, the model and solution of the DCO carbon emission optimization problem are completed in Section 4. A case study is conducted and discussed in Section 5. Finally, Section 6 concludes this paper.

Problem Formulation
A DCO operating N geographically distributed IDCs is considered. These IDCs are interconnected through the data transmission network, and the workloads can be transmitted between different IDCs, as shown in Figure 1. Due to the structure of power systems, they are located in different regional power grids or distribution networks with their own local RES and EVCS. There is no power flow among different IDCs. A generalized flexible load dispatching potential evaluation method for the and their surrounding road networks based on the Minkowski sum algorithm posed. Combined with the Road Impedance (RI) model, the dispatching po assessment of the EVCS is more realistic.

•
Based on the proposed strategy, the impact of the time delay of processing wo in the IDC on the carbon emission and renewable energy consumption capa analyzed.

Organization of This Article
The remainder of this paper is organized, as follows. Section 2 introduces th patching strategy of the system and DCO carbon emission model. The IDC powe sumption model and the EVCS dispatching potential model considering traffic flo presented in Section 3. Subsequently, the model and solution of the DCO carbon em optimization problem are completed in Section 4. A case study is conducted and dis in Section 5. Finally, Section 6 concludes this paper.

Problem Formulation
A DCO operating N geographically distributed IDCs is considered. These ID interconnected through the data transmission network, and the workloads can be mitted between different IDCs, as shown in Figure 1. Due to the structure of pow tems, they are located in different regional power grids or distribution networks wit own local RES and EVCS. There is no power flow among different IDCs. The IDCs consume a lot of power due to the processing workload and they get from the utility grid, where the carbon emission comes from. To reduce the amo carbon emission, integrated RES including photovoltaic resources and wind pow sources are introduced to the local IDC as one of the main energy sources. There a ways for IDCs to purchase renewable energy. One is to build solar panels, wind tu and supporting energy storage equipment to provide energy for itself. Another wa buy one-day power generation from local RES aggregators. This paper adopts the s way. With the growing amount of EVs, the demand for EVCS can be introduced to f reduce carbon emission. The EVCS contains multiple chargers that can be charged The IDCs consume a lot of power due to the processing workload and they get power from the utility grid, where the carbon emission comes from. To reduce the amount of carbon emission, integrated RES including photovoltaic resources and wind power resources are introduced to the local IDC as one of the main energy sources. There are two ways for IDCs to purchase renewable energy. One is to build solar panels, wind turbines, and supporting energy storage equipment to provide energy for itself. Another way is to buy one-day power generation from local RES aggregators. This paper adopts the second way. With the growing amount of EVs, the demand for EVCS can be introduced to further reduce carbon emission. The EVCS contains multiple chargers that can be charged by EV users. Through the unified management of the EVCS, the whole charging station is modeled as a flexible load that can participate in the joint dispatching of the IDC. As the load aggregator, the economic driving force of the EVCSs participating in joint dispatching needs to be discussed. The EVCSs provide auxiliary services for the DCO and receive the remuneration paid by the DCO. This part of the reward is called capacity price [27]. This tariff is paid for the energy regulated by the EVCSs in response to DCO dispatching instructions. For the DCO, although the EVCSs are paid for auxiliary services, joint dispatching reduces the power purchase of the DCO from the power grid, thereby reducing carbon emissions and power purchase costs. Therefore, this strategy will bring benefits to the EVCSs and the DCO. The scenario premise of this paper is that the chargers in the EVCSs and the intelligent control equipment participating in the joint dispatching are existing, and the EVCSs do not need to pay additional equipment construction costs for this strategy.
For sustainable IDCs, their energy comes mainly from the photovoltaic and wind power resources they purchase, and the energy consumption is mainly used for the operation of the IDC itself, including the operation of servers, refrigeration systems, lighting infrastructure, etc. [28] and the charging power of EVs in the charging station. After the temporal dispatching of the EVCS and the spatial dispatching of the IDC workload, the part of the energy shortage of the proposed system is obtained from the utility grid. Therefore, the power balance constraint of the IDC i (1 ≤ i ≤ N) is established as Equation (1): where P PV t,i is the average power generated by the photovoltaic power supply in the IDC i in slot t. P opr t,i is the average operating power consumption of processing workload in slot t of the IDC i. P EVCS t,i is the average charging power of the EVCS i connected to the IDC i in slot t.
It can be assumed that the carbon emissions from wind power and PV are negligible. It is known that the large amount of energy consumption generated by the driving and charging process of EVs will produce carbon emissions that cannot be ignored. However, in this paper, the EVCS is not directly connected to the utility grid to get power. It is integrated into the power supply and consumption balance system of the IDC so that the energy required for EV charging will be included in the energy consumption of the IDC. Therefore, there is no need to consider the carbon emission of EVs.
The total carbon emission C t generated by the DCO with N IDCs in slot t can be calculated as the sum of the carbon emissions generated by the power purchased by each IDC i from the utility grid in slot t. The equation of carbon emission C t is as follows: where e is the carbon emission rate of utility grid; ∆t is the duration of time slot; P grid t,i is the average power injected by the utility grid into the IDC i in slot t.
The calculation methods and constraints of the IDC operation power consumption model and the EVCS aggregation model will be discussed in detail in the next section.

Models for the IDC and the Traffic Flow EVCS
In this section, the mathematical models of the main components in the system are listed, including the IDC operation power consumption model and dispatching potential aggregation model of the EVCS considering traffic flow.

The IDC Operation Power Consumption Model
The operation power consumption of the IDC mainly includes server power consumption, cooling power consumption, lighting, and other auxiliary infrastructure power consumption. Among them, server power consumption is the energy consumed by the processing workload. The power consumption of the IDC can be regarded as a linear function of server power consumption. For the IDC i with M servers, the equation is as follows: where P idle i is the fixed power consumption necessary to maintain the operation of the IDC i; k i is the power consumption increased by every increase in the number of active servers of the IDC i, which is set as a constant in this paper; m t,i is the number of active servers in the IDC i in slot t, and it meets: where M i is the total number of servers in the IDC i. The workload of the IDC includes batch workload and interactive workload. The batch workload is also called offline load. After it reaches the IDC, it does not need to be processed immediately, usually within a few hours. The interactive workload needs to be processed within a very short specified delay, and due to the "redundant backup" technology, they can be migrated between different IDCs, as described in Section 1. The workload mentioned in this article is interactive. Therefore, this strategy is only applicable to interactive workloads. The IDCs provide data processing services to users, and the time delay limit D of data processing will be specified in the service level agreement (SLA) signed in advance. The data load is distributed by the front-end to the active server in the IDC, and its average stay time in the IDC cannot exceed the time delay limit D. µ i is set as the service rate of a single active server in the IDC i, which is a fixed parameter depending on server performance. The constraint on time delay can be written as follows: where L t,i is the is the amount of data load reaching the IDC i in slot t. From Equation (3), the relationship between P opr t,i and m t,i can be deduced and substituted into Equation (4) to obtain Equation (6). By substituting it into Equation (5) and summing i on both sides of the formula, Equation (7) can be obtained.
where µ i , D, k i , P idle i and M i are fixed parameters for the IDC operation.
L t,i represents the total amount of workload that the DCO needs to process in slot t.L t,i represents the data load that the IDC i needs to process in slot t after dispatching. Equation (6) gives the limits of power consumption of the IDC i. Equation (7) shows the relationship between the power consumption of the IDC i and the amount of workload and the average residence time of workload in the IDC i. Equation (8) specifies that the total amount of workload remains unchanged before and after dispatch. From the above equations, it can be concluded that when the processing workload of a single IDC i remains the same, the shorter the average residence time of the data load in the IDC i (i.e., shorter processing time), the greater the power consumption generated by the IDC i. When the average residence time of all workload in the IDC i is equal to the time delay limit D specified by the user and the IDC i, the constraint on time delay is: In this case, Equation (7) can be written as: That is, the power consumption of the IDC i is linear to the number of workload.

Evaluation Model of Dispatching Potential of the EVCS Considering Traffic Flow
In this subsection, a method to form the aggregated EV model in an EVCS is proposed. The DCO purchases EV flexibility resource services from the EVCS which need to be aggregated. In this paper, the Minkowski sum method is used to model. This method can evaluate the dispatching potential of the EVCS while retaining the original variables and constraints of EVs. In addition, since the charging station is open to the public, the dispatching potential of the EVCS considers the EVs driving in the road network around the EVCS, and a generalized flexible load model considering the road network is established. The DCO can effectively use the dispatching resources of EV polymers on the premise of knowing as little user privacy information as possible.
The following assumptions are considered here: 1.
The EVCSs obtain information from the intelligent transportation system (ITS) [29], which includes the historical flow data of the road network and the real-time monitoring data of vehicles. The electric vehicle terminal will upload the real-time battery power information, vehicle location information, and the maximum charging power allowed by the battery to the ITS.

2.
When an EV is connected to the charger in the EVCS, it is considered to obey the dispatching instructions of the EVCS. The specific implementation of dispatching command and control issued by the EVCS to a single charger is not discussed in this paper.

3.
Discrete time is applied in the proposed model [30]. The ITS provides the EVCSs with the monitoring data of EVs that need to be charged. The EVCSs calculate their capacity and power changes every 15 min.
Traffic flow is closely related to road congestion. Therefore, this paper assumes that the average congestion x t of the road near the EVCS represents the average distance between every two adjacent vehicles on the road. Combined with the road type and road length in the actual road network, the traffic flow in the nearby area at a certain time, i.e., the total number of electric vehicles, can be calculated as follows: where J i,t is the total number of vehicles in the area near the EVCS i in slot t; R i is the total length of roads in the area near the EVCS i; type represents each road type in the area near the EVCS i, which can be obtained in Equation (12): That is, the main trunk road indicates that there are two-way eight lanes; the secondary trunk road indicates that there are two-way four lanes; the branch road indicates that there are two-way two lanes.
This paper takes a part of the actual road network in Beijing, China as the research object. ArcMap10.7 is used to extract and simplify the critical roads from the actual road network. The processed road network map is classified according to Equation (12). As well, there are three charging stations in the simplified road network, as shown in Figure 2b. Figure 2a shows the actual road network in ArcMap10.7 in the blue rectangle in Figure 2b. This paper takes a part of the actual road network in Beijing, China as the research object. ArcMap10.7 is used to extract and simplify the critical roads from the actual road network. The processed road network map is classified according to Equation (12). As well, there are three charging stations in the simplified road network, as shown in Figure  2b. Figure 2a shows the actual road network in ArcMap10.7 in the blue rectangle in Figure  2b. When the historical congestion of the road is known, the data-driven method can be used to predict the total number of vehicles at a certain time in the future. Time series prediction models such as SARIMA and LSTM can be used to predict road congestion. This part is not the focus of this paper and will not be described. The EV j position information j L and battery state is lower than the threshold HT , the EV j needs to be charged. As the EV will consume electricity and time during driving, the BPR road impedance (RI) model is introduced to calculate the SOC loss and arrival time of the EV to be charged from the current location to the charging station for the EVCS dispatching potential evaluation.

Road Impedance (RI) Model
When comparing a road to a circuit, the capacity of each road can be compared to the conductivity of it, the loss on the road is similar to the impedance in the circuit, which is called road impedance. There are many factors affecting RI, including traffic time, traffic cost, road speed limit, weather environment, toll stations, and urban nodes, but most of them can be indirectly transformed into time-dependent. Therefore, traffic time is used to define RI. The most used RI model is the Federal Highway Administration RI model (BPR), which is defined as: where r t is the impedance of road r ; t x is the road traffic flow; 0 t is the passage time with no impedance; v is the average speed; c is the traffic capacity of the road and recorded as the number of EVs passing per hour; λ and δ is the impedance parameter, recommended as  When the historical congestion of the road is known, the data-driven method can be used to predict the total number of vehicles at a certain time in the future. Time series prediction models such as SARIMA and LSTM can be used to predict road congestion. This part is not the focus of this paper and will not be described. The EV j position information L j and battery state SOC j,moni at a certain time are monitored. According to the charging habits and psychology of EV users, when the battery SOC j,moni is lower than the threshold HT, the EV j needs to be charged. As the EV will consume electricity and time during driving, the BPR road impedance (RI) model is introduced to calculate the SOC loss and arrival time of the EV to be charged from the current location to the charging station for the EVCS dispatching potential evaluation.

Road Impedance (RI) Model
When comparing a road to a circuit, the capacity of each road can be compared to the conductivity of it, the loss on the road is similar to the impedance in the circuit, which is called road impedance. There are many factors affecting RI, including traffic time, traffic cost, road speed limit, weather environment, toll stations, and urban nodes, but most of them can be indirectly transformed into time-dependent. Therefore, traffic time is used to define RI. The most used RI model is the Federal Highway Administration RI model (BPR), which is defined as: where t r is the impedance of road r; x t is the road traffic flow; t 0 is the passage time with no impedance; v is the average speed; c is the traffic capacity of the road and recorded as the number of EVs passing per hour; λ and δ is the impedance parameter, recommended as λ = 0.15, δ = 4. Therefore, the RI model is proportional to the current traffic flow of the road. The passage time of the EV passing through the road is calculated by the RI model, and the time to reach the charging station and the SOC are calculated by the driving mileage and power consumption.
T j,arrive = T j,moni + t j,r , where T arrive and SOC j,arrive represent the arrival time and SOC of EV j, respectively; T j,moni and SOC j,moni are the monitoring time and SOC of EV j, respectively; ∆s j is the average driving power consumption of EV j. Due to the difference in charging demand, the EVs connected to the grid can be divided into two categories. One is the EVs that need fast charging, such as taxis. They have a short charging time and the power cannot be adjusted. It is difficult to participate in the dispatching strategy. The other category is the private EVs that can be connected to the grid at night or during long working hours, and their charging power is adjustable. Therefore, this category of load can be dispatched and participate in flexible load dispatching. This paper is mainly focused on the second category of EVs.
Equation (17) is the charging power constraint; Equation (18) is the timing constraint of battery capacity; Equation (19) is the battery capacity limit constraint.
To sum up, the parameter set of an EV driving on the road and about to charge at the charging station is T moni , SOC moni , L, P cha,max , T leave , s max .
This set includes user location information, battery power information, expected departure time, and other information related to user privacy. The EVCS cannot share the information with the IDC for joint energy dispatching. Therefore, the load aggregator should process this information.

The Generalized Flexible Load Model Based on Minkowski Sum
The dispatching power of the generalized charging station needs to be superimposed on the individual decisions of EVs in the area to form an envelope space, which can hide the user's privacy information while retaining the constraints of individual EVs, and reduce the number of variables and constraints. Minkowski sum is suitable for calculating the sum of point sets in Euclidean space, and its physical essence is the expansion set of multiple spaces [31]. The Minkowski sum of two sets A and B is as follows: The Minkowski sum of the two convex hulls is the envelope of the two variable spaces, as shown by the red line in Figure 3. The area formed by G in Figure 3 is defined as a generalized charging station decisions space. It establishes a flexible load model for the charging station. As well, its parameters are the dispatching potential parameters of the charging station. When the decision variables are two-dimensional, the model forms an area. When the decision variables  The area formed by G in Figure 3 is defined as a generalized charging station decisions space. It establishes a flexible load model for the charging station. As well, its parameters are the dispatching potential parameters of the charging station. When the decision variables are two-dimensional, the model forms an area. When the decision variables have higher dimensions, the model forms a hypercube space.
Due to the different grid access time slots of each EV, the definition domain of variable space is different. Therefore, the Minkowski sum cannot be performed directly. When the EV j is connected to the charger q, the information of the EV j is given to the charger q. A state variable X q,t is introduced to characterize the activation state of the charger q: Thus, the time definition domain of each charging point is extended from the active state to the full dispatching period, so that it can be calculated by Minkowski sum. The energy variation constraint of the charger q can be written as: s t,q = s t−1,q + s q,arrive X t,q X t,q − X t−1,q − s q,leave X t−1,q X t−1,q − X t,q + η cha P cha t,q ∆t, ∀t ∈ T.
The electric quantity variation of the charger q is described as follows: ∆s t,q = s q,arrive X t,q X t,q − X t−1,q − s q,leave X t−1,q X t−1,q − X t,q .
The total charging dispatching power of the EVCS i is the sum of the maximum charging power of each charger, which can be expressed as: where P cha,max j represents the maximum charging power of the charger q. Similarly, the minimum and maximum of the total dispatching power of the EVCS i are as follows: where s min q and s max q represent the minimum and maximum electric quantity of the charger q, respectively. Therefore, the parameter of the EVCS generalized flexible load model can be defined as: And the constraints of the model can be expressed as: This method greatly reduces the dimension of the model. The hypercube space contains all feasible charging decisions of the EVCS. The above model is the envelope space composed of all charging decisions of the EVCS. In the case study, a baseline needs to be set up. In this baseline, the temporal dispatching ability of EVs is not considered. This requires some corrections to some parameters reported in the above model.
In this case, the EV drives into the EVCS to charge at the maximum power and automatically disconnects the power connection with the piles after it is fully charged. Due to the change of grid access time, the X matrix, in this case, will change, resulting in the change of other parameters related to X. Then, the activation time of the charger q can be calculated as: where T stay q is the activation time of the charger q in the case of maximum power charging, so the X matrix can be modified according to Equation (22). The new upper limit of charging power P cha,max i,t of the EVCS i is obtained by Equation (25). The constraint condition only needs to retain the following equation: It should be noted that the activation time of the charger (electric vehicle grid access time) has been considered in the above equation. Equation (33) is the constraint that the EVCS needs to comply with in Baseline1 in Section 5.

Optimization Problem and Solution of DCO Carbon Emission Reduction
Based on the IDC and the EVCS model in Section 3, a spatial-temporal coupling dispatching mechanism for carbon emission reduction is proposed, which takes into account the spatial migration capacity of workloads in the IDCs and the temporal dispatch capacity of the generalized flexible load of the EVCS. The spatial migration of workload among the IDCs is completed while meeting the constraints of energy balance, server power consumption, and user time delay requirements. Considering the temporal dispatch capacity of the generalized flexible load model of the EVCS, it can better match the output curve of RES and further reduce carbon emissions.
Therefore, according to the analysis in Section 2, to achieve the sustainability of geographically dispersed IDCs, that is, the carbon emission reduction target, the total carbon emission of each single IDC owned by the DCO needs to be calculated, and the total optimization target is formulated as follow: where e = 0.968kg · (kW · h) −1 [32]. By establishing the spatial and temporal migration mechanism, the power consumption of the IDC can track the maximum output of RES to the greatest extent in each control period to reduce the power purchase from the utility grid, and obtain the minimum carbon emission of the DCO in Equation (34).
The constraints of the spatial-temporal coupled dispatching optimization problem are as follows: 1.
The constraints of the IDC power supply and consumption balance.
For each single IDC, the balance of energy supply and consumption must be met. In Section 2, the power supply and consumption of the single IDC are stated. The actual consumption of renewable energy is also constrained by the current maximum output. 2.
The constraints of the IDC energy consumption and workload processing time delay.
It should be noted that when the average stay time of the workload in the IDC is exactly equal to the time delay limit D specified by the user and the IDC, the workload time delay constraint should be written as: 3.
The conservation constraint of the total workload.
The total workload handled by the DCO before and after load dispatching must remain unchanged, and the workload can be transferred between the IDCs: 4. The constraints of the EVCS.
According to the above constraints, it can be found that the difference between the workloads of the IDCs and the traditional flexible load is that the constraints of the workloads of the IDCs are spatial, as shown in Equation (40). While the constraint of the EVCS is temporal, as shown in Equation (42). The decision variables of the model are the power consumption P opr t,i of the IDC i, the workload L t,i allocated to the IDC i after migration, the electricity S t,i of the EVCS i and the charging power P EVCS t,i of the EVCS i. Through the modeling work in Section 3, a large number of EVs that may produce a large number of constraints and variables is aggregated into the generalized flexible load model, which greatly reduces the dimension of the model and the complexity of the joint dispatching solution. In this process, Boolean variables X t,q are generated. Through further analysis, they are encapsulated in the variables ∆S t,q , to avoid the problem becoming a mixed integer programming problem and further reduce the difficulty of solving.

Case Study
In this Section, the spatial-temporal coupling dispatching problem of carbon emission reduction in DCO is validated, which takes into account the workload spatial migration capacity of three IDCs and the temporal dispatching capacity of a generalized flexible load of three EVCSs (N = 3). Simulation is based on the real workload curves [33]. The time slot is set as ∆t = 15min, T = 96 within one dispatch day. For performance comparisons, the following strategies are adopted as the cases: • Baseline1(B1): The workload of the IDC is processed locally, that is, the spatial migration of workload is not considered. The EV does not consider the orderly charging with adjustable power, that is, after each EV arrives at the charging station, it will be charged at the maximum power, and the charging will be stopped automatically after it is full charged. Therefore, the load of the EVCS will no longer have the flexibility of temporal adjustment. • Baseline2(B2): The workload of the IDC is processed locally, that is, the spatial migration of workload is not considered. As a dispatchable flexible load, the temporal dispatching ability of the EVCS under constraints is considered. • Proposed(P1): Considering the spatial migration ability of workloads between IDCs and the temporal dispatching ability of generalized flexible load of the EVCS, the spatial-temporal coupling joint dispatching of the IDCs is carried out.

Carbon Emissions and RES Consumption of the IDCs
The curves of RES purchased locally by the three IDCs are shown in Figure 4. The experiments compare the RES consumption and carbon emissions of the above three dispatching strategies. Figure 4 shows the matching results between the power consumption of each IDC and the maximum power generation of RES in one dispatching day. It can be derived from Figure 4 that without migration dispatching, due to the fluctuation of RES output, the RES supply in some periods cannot meet the power consumption of the IDC, as shown in dispatching slots (0-33) and (66-96) of the IDC1 in Figure 4a. During this two periods, the power consumption of the IDC is higher than the output of local RES. The part of the power shortage needs to be supplemented by thermal power generation from the utility grid, resulting in large carbon emissions. In the period with high RES output, as shown in dispatching slots (34-65) of the IDC1 in Figure 4a. During this period, the oversupply of RES has led to the abandonment of wind and solar, resulting in a low consumption rate of renewable energy.  In B2, the flexibility resource of the EVCS is introduced. It can be seen in Figure 4 that the dispatching strategy of introducing flexible resources of the EVCS has a higher consumption rate of RES. Due to the impact of the current charging station infrastructure construction, the equipment capacity of the EVCS adopted in our case is limited, but as a pilot, it has produced the effect of improving the consumption rate of RES. In the future, In B2, the flexibility resource of the EVCS is introduced. It can be seen in Figure 4 that the dispatching strategy of introducing flexible resources of the EVCS has a higher consumption rate of RES. Due to the impact of the current charging station infrastructure construction, the equipment capacity of the EVCS adopted in our case is limited, but as a pilot, it has produced the effect of improving the consumption rate of RES. In the future, with the increase in the number of EVs, the increase in battery capacity and charging power of EVs, and the improvement of charging station infrastructure, the use of the EVCS flexibility resources for the IDC carbon emission reduction strategy will bring greater benefits. On this basis, the spatial-temporal coupling optimization strategy proposed in this paper has greater consumption of RES than the B2 strategy, as shown in Figure 4. Among the three strategies, the total power consumption of the IDC remains unchanged, but the B2 and P1 strategies consume more RES, which reduces the power purchased from the utility grid and thus reduces carbon emissions. The carbon emissions of the three strategies are shown in Figure 5. In order to verify the economic feasibility of this strategy, the economy of DCO under this strategy needs to be discussed. The cost of DCO in B1 and P1 strategies are compared. The cost of the DCO includes the cost of purchasing RES, the cost of purchasing electricity from the utility grid, and the capacity cost of auxiliary services for the EVCSs. The purchase price of PV is 5.229 cents per kWh, the purchase price of wind power is 4.333 cents per kWh [34], and the purchase price of the utility grid is 5.827 cents per kWh [35]. The capacity price of the EVCSs participating in the auxiliary service market is 5.976 cents per kWh [36]. It is calculated that under the B1 strategy, the costs of the three IDCs are US $6311.80, $5582.18, and $5806.18, respectively, and the total cost of the DCO is $17,700. 16. Under the P1 strategy, the costs of the three IDCs are $6162.70, $5506.40, and $5838.02, respectively, and the total cost of DCO is $17,507.11. The P1 strategy reduces the DCO's total operating cost by 1.09% while reducing its carbon emissions. Among them, the cost of the IDC 3 has increased due to the high matching degree of the IDC 3's original energy consumption and renewable energy output curve, and the effect of the P1 strategy on reducing power purchase of the utility grid is not obvious. However, we need to treat the DCO as an entire stakeholder in this strategy. The proposed strategy can reduce the cost of the DCO in this scenario. In addition, the reduction of carbon emissions enables DCO to participate in the green certificate market and obtain more benefits, which is beyond the scope of this article.

Analysis of the EVCS Dispatching Result and RI Model Influence
The analysis of the dispatching of the EVCSs is carried out under the spatial-temporal coupling dispatching strategy (P1) proposed in this paper. Figure 6 shows the maximum energy limit max , i t S and minimum energy limit min , i t S of the three EVCSs, as well as their energy and power changes in the dispatching process. Different from flexible loads such as air conditioning, the temporal charging flexible load of the EVCS has its battery capacity limit. Therefore, its capacity change must be considered in dispatching. The proposed The results show that the total carbon emissions of DCO simulated by the three strategies are 22,448.27 kg, 20,535.29 kg, and 19,580.20 kg. The carbon emission of the B2 strategy, which introduces the flexibility resources of the EVCSs, is 8.52% lower than that of the B1 strategy. Under the proposed spatial-temporal coupling dispatching strategy P1, the carbon emission is further reduced, which is 4.65% lower than B2 and 12.78% lower than B1. This proves that the established generalized flexible load model of the EVCS and the proposed spatial-temporal coupling model considering the spatial migration capacity of the workloads of the IDCs and the temporal dispatch capacity of the generalized flexible load of the EVCSs are effective in reducing carbon emissions of the IDCs.
In order to verify the economic feasibility of this strategy, the economy of DCO under this strategy needs to be discussed. The cost of DCO in B1 and P1 strategies are compared. The cost of the DCO includes the cost of purchasing RES, the cost of purchasing electricity from the utility grid, and the capacity cost of auxiliary services for the EVCSs. The purchase price of PV is 5.229 cents per kWh, the purchase price of wind power is 4.333 cents per kWh [34], and the purchase price of the utility grid is 5.827 cents per kWh [35]. The capacity price of the EVCSs participating in the auxiliary service market is 5.976 cents per kWh [36]. It is calculated that under the B1 strategy, the costs of the three IDCs are US $6311.80, $5582.18, and $5806.18, respectively, and the total cost of the DCO is $17,700. 16. Under the P1 strategy, the costs of the three IDCs are $6162.70, $5506.40, and $5838.02, respectively, and the total cost of DCO is $17,507.11. The P1 strategy reduces the DCO's total operating cost by 1.09% while reducing its carbon emissions. Among them, the cost of the IDC 3 has increased due to the high matching degree of the IDC 3's original energy consumption and renewable energy output curve, and the effect of the P1 strategy on reducing power purchase of the utility grid is not obvious.
However, we need to treat the DCO as an entire stakeholder in this strategy. The proposed strategy can reduce the cost of the DCO in this scenario. In addition, the reduction of carbon emissions enables DCO to participate in the green certificate market and obtain more benefits, which is beyond the scope of this article.

Analysis of the EVCS Dispatching Result and RI Model Influence
The analysis of the dispatching of the EVCSs is carried out under the spatial-temporal coupling dispatching strategy (P1) proposed in this paper. Figure 6 shows the maximum energy limit S max i,t and minimum energy limit S min i,t of the three EVCSs, as well as their energy and power changes in the dispatching process. Different from flexible loads such as air conditioning, the temporal charging flexible load of the EVCS has its battery capacity limit. Therefore, its capacity change must be considered in dispatching. The proposed generalized flexible load model of the EVCS based on the Minkowski sum algorithm can describe the capacity boundary of the EVCS and make it meet the capacity demand of users in dispatching. On this basis, the aggregate charging power of the charging station is dispatched in the time dimension, focusing on the time when the RES output is high at noon, to better improve the RES consumption rate and reduce the carbon emission of the IDC. The RI model proposed in this paper describes the time loss and capacity loss of each EV driving in the road network to the charging station. As well, the grid access time and capacity state when an EV connects to a charger can be described more accurately, which has an impact on describing the dispatching potential of the EVCS. Figure 7 shows the deviation of generalized flexible load dispatching potential caused by considering the RI model. If the road network area is larger, it will lead to a longer vehicle driving distance, resulting in more time loss and capacity loss. If the traffic flow increases, i.e., the number of vehicles increases, the slower the vehicle will travel, which will affect the calculation of vehicle time loss and capacity loss. The RI model proposed in this paper describes the time loss and capacity loss of each EV driving in the road network to the charging station. As well, the grid access time and capacity state when an EV connects to a charger can be described more accurately, which has an impact on describing the dispatching potential of the EVCS. Figure 7 shows the deviation of generalized flexible load dispatching potential caused by considering the RI model. If the road network area is larger, it will lead to a longer vehicle driving distance, resulting in more time loss and capacity loss. If the traffic flow increases, i.e., the number of vehicles increases, the slower the vehicle will travel, which will affect the calculation of vehicle time loss and capacity loss.

Analysis of Workload Processing Time Delay and the Impact on RES Consumption
In the comparative analysis in Section 5.1, the power changes of the three strategies are shown. Among them, the P1 case has a positive effect on the improvement of renewable energy consumption rate, which shows that considering the migration ability of the IDC workload has a positive effect on the IDC's carbon emission reduction. To reveal the workload dispatching in the IDC under this strategy, the workload processed by the IDC before and after migration dispatching is analyzed, as shown in Figure 8. When the total workload processed by DCO remains unchanged, the workload in the IDC can be transferred among the interconnected IDCs. The migration of workload leads to the migration of power consumption. Therefore, the workload can migrate to the IDC with rich renewable energy output.

Analysis of Workload Processing Time Delay and the Impact on RES Consumption
In the comparative analysis in Section 5.1, the power changes of the three strategies are shown. Among them, the P1 case has a positive effect on the improvement of renewable energy consumption rate, which shows that considering the migration ability of the IDC workload has a positive effect on the IDC's carbon emission reduction. To reveal the workload dispatching in the IDC under this strategy, the workload processed by the IDC before and after migration dispatching is analyzed, as shown in Figure 8. When the total workload processed by DCO remains unchanged, the workload in the IDC can be transferred among the interconnected IDCs. The migration of workload leads to the migration of power consumption. Therefore, the workload can migrate to the IDC with rich renewable energy output. In the above cases, the waiting time of the workload in the IDC is assumed to equal the time delay limit 100ms D = signed between the user and the DCO. That is, Equation (10) in Section 3 is satisfied. In this subsection, the impact of workload processing time delay on RES consumption and carbon emission of the IDC under the spatial-temporal coupling dispatching model is discussed. In the simulation, the Equation (37) in the original model is changed to Equation (39). That is, the IDC is allowed to process the workload In the above cases, the waiting time of the workload in the IDC is assumed to equal the time delay limit D = 100 ms signed between the user and the DCO. That is, Equation (10) in Section 3 is satisfied. In this subsection, the impact of workload processing time delay on RES consumption and carbon emission of the IDC under the spatial-temporal coupling dispatching model is discussed. In the simulation, the Equation (37) in the original model is changed to Equation (39). That is, the IDC is allowed to process the workload more efficiently with a higher power, shorten the workload processing time delay as much as possible and improve the user experience on the premise of meeting the carbon emission reduction target. Figure 9 shows the IDC server power and the corresponding workload processing time delay in the cases of D = 100ms and D < 100ms. The simulation results show that when the time delay satisfies D < 100ms, the overall power consumption of the server in the IDC 1 on a dispatching day increases by 1.66%, and the total carbon emission of DCO on the dispatching day remains unchanged. Figure 9 can well explain this phenomenon. Due to the existence of RES, there is a part of unused energy during the peak period of power generation. This part of the energy can be used to improve the calculation efficiency of the IDC, reduce the workload processing time delay, and improve the user experience while ensuring that the total carbon emission remains unchanged.

Conclusions
In this paper, the impact of the interactive workload migration mechanism of the IDCs and the temporal dispatch ability of generalized flexible load of the EVCSs on DCO carbon emission reduction is evaluated. A generalized flexible load model of the EVCSs considering traffic flow is proposed based on the Minkowski sum algorithm. A DCO with three geographically dispersed IDCs and three EVCSs in the Beijing road network are used to demonstrate the effectiveness of the proposed model.
The simulation results show that the joint dispatching of the EVCSs temporal regu-

Conclusions
In this paper, the impact of the interactive workload migration mechanism of the IDCs and the temporal dispatch ability of generalized flexible load of the EVCSs on DCO carbon emission reduction is evaluated. A generalized flexible load model of the EVCSs considering traffic flow is proposed based on the Minkowski sum algorithm. A DCO with three geographically dispersed IDCs and three EVCSs in the Beijing road network are used to demonstrate the effectiveness of the proposed model.
The simulation results show that the joint dispatching of the EVCSs temporal regulation ability and interactive workload spatial migration ability can reduce the carbon emission of the DCO. The proposed strategy can also reduce the energy cost paid by the DCO. The aggregation model of the EVCS considering traffic flow and road resistance meets the technical needs of the social trend of continuous expansion of the scale of EVs. This paper can provide technical support for the DCO carbon emission reduction with multiple geographically dispersed IDCs.
Future research of this strategy could be continued in two directions, considering the strategy of the DCO participating in the green certificate market and the strategy of the EVCS controlling EVs in real time. Participating in the green certificate market can bring more benefits to the DCO and give it a greater driving force to implement carbon emission reduction strategies. As well, the research on the strategy of the EVCS issuing dispatching instructions to each charger in real time is conducive to the application and implementation of the research strategy in this paper.