Joint Virtual Energy Storage Modeling with Electric Vehicle Participation in Energy Local Area Smart Grid

In this research, the joint virtual energy storage modeling with electric vehicle participation in energy local area Smart Grid is considered.-is article first constructs a virtual energy storage model and a joint virtual energy storage model for air conditioning and electric vehicles. -erefore, for the optimization problem of virtual energy storage power, a continuous rolling optimization algorithm to determine the feasible solution of the high-dimensional complex constraint optimization problem is proposed to solve the optimization problem. Finally, the analysis, for example, illustrates the economics of joint virtual energy storage in the Smart Grid. -e results prove that air conditioning and electric vehicles have the ability to jointly participate in virtual energy storage, and the comparison proves that joint virtual energy storage can effectively improve the economics of electricity consumption.


Introduction
Due to the access of large-scale distributed power equipment, especially the strong randomness of wind power and solar energy, energy storage capacity is of vital importance in the Smart Grid [1][2][3][4][5]. However, batteries are expensive and easily cause pollution [6]. In recent years, a large number of scholars have put forward the concept of virtual energy storage and thus have studied the virtual energy storage with air conditioning participation. e study in [7] analyzes the way of virtual energy storage in the Smart Grid from two aspects: electricity price control and direct control. e study in [8] used the building's heat storage capacity to establish a building-based virtual energy storage system model, by managing the charging (or discharging) power of the virtual energy storage system according to the user's indoor temperature limit and applying it to the microgrid; this literature successfully improved running economy. e study in [9] established a thermodynamic model of air conditioning load with the air conditioning set value as the dependent variable and used it to absorb the strong random distributed new energy power generation equipment output, which improved the energy utilization as well as the stability of the system operation. e study in [10] conducted an in-depth analysis of the direct controllability of some demand-side loads and thus dug into the potential of air conditioning load virtual energy storage. Based on the premise of satisfying the user's comfort requirements, relying on its virtual energy storage to suppress the output of distributed wind turbines with strong randomness, in order to reduce the configuration capacity of energy storage equipment and reduce the cost of energy Internet management regulation, for the controllable load equipment air conditioning, with consideration of its load-adjustable space, the study in [11] established an air conditioning load double-layer optimal scheduling and control model and maximized the interests for both parties through the coordination and optimization between the macrolevel power company and the microlevel direct load control agents. e study in [12][13][14] improves the economy of energy consumption through the demand response of air conditioning and refrigerator, respectively, on the premise that user comfort can be guaranteed. Although the virtual energy storage is not explicitly proposed, describing control objects with controllable load and demand response instead, the control method is consistent with that in the research of virtual energy storage.
Furthermore, there are many optimization algorithms applied to multienergy optimization [15,16] in Smart Grid. In literature [17], a hybrid constraint handling strategy (HCHS) based on nondominated sorting genetic algorithm II (NSGAII) is proposed to deal with the typical constraints, by which the constraint violations can be removed in several steps during the evolutionary process. e study in [1] proposed an artificial shark optimization (ASO) method to remove the limitation of existing algorithms for solving the economical operation problem of microgrid. e study in [18] presents a new metaheuristic optimization algorithm, the firefly algorithm, and an enhanced version of it, called chaos mutation firefly algorithm (CMFA), for solving power economic dispatch problems while considering various power constraints such as valve-point effects, ramp rate limits, prohibited operating zones, and multiple generator fuel options. In literature [19], an improved genetic algorithm is proposed to overcome the obstacle of infeasibility and exhibit better economic dispatch. e study in [20] proposes a new distributed method, namely, the neuralnetwork-based Lagrange multiplier selection (NN-LMS), to prominently reduce the iterations and avoid an oscillation. A calibrated building simulation model was developed in literature [21] and used to assess the effectiveness of demand response strategies under different time-of-use electricity tariffs in conjunction with zone thermal control.
is article first constructs a virtual energy storage model and a joint virtual energy storage model for air conditioning and electric vehicles. erefore, for the optimization problem of virtual energy storage power, a continuous rolling optimization algorithm to determine the feasible solution of the high-dimensional complex constraint optimization problem is proposed to solve the optimization problem. Finally, the analysis, for example, illustrates the economics of joint virtual energy storage in the Smart Grid.  Figure 1. e virtual energy storage of air conditioning and building can be abstracted as a closed-loop system with feedback as shown in Figure 2.

ermodynamic Model of ermoelectric Conversion of Air Conditioning.
e commonly used thermodynamic equivalent equations [22,23] are expressed as follows: In equation (1), T in (t) and T out (t) are the indoor and outdoor temperatures, respectively, in units of°C; Q ac (t) is the air conditioning cooling capacity, in units of kW; R is the equivalent thermal resistance, in units of°C/kW; C is the building equivalent heat capacity, in units of kJ/°C; S is the building area. e above environmental parameters all have an impact on the virtual energy storage capacity of the air conditioning.
e relationship between the electric power of the air conditioning equipment and the cooling capacity is as follows: P ac (t) is the air conditioning power in kW; η is the air conditioning thermoelectric conversion coefficient.
Besides, the expression of air conditioner virtual energy storage rated capacity S ac− virtual− N is From the above derivation, the relationship among the remaining power of virtual energy storage, the rated capacity, and the equivalent heat capacity of the building as well as the comfortable temperature range of the human body can be clearly shown.

Natural Power Consumption in Virtual Energy
Storage with Air Conditioning Participation. For the virtual energy storage with air conditioner participation, when the indoor and outdoor temperatures are different, there is always heat exchange in the air conditioning and building virtual energy storage system. At this time, the air conditioning needs a certain amount of power to maintain the indoor temperature equal to the set temperature, that is, natural power consumption for virtual energy storage of air conditioning: P ac− o (t) will change according to the change of indoor and outdoor temperatures. When the air conditioner power is P ac− o (t), the balance of indoor and outdoor temperature can be maintained.

Virtual Charge and Discharge Power of Virtual Energy
Storage with Air Conditioning Participation. When the air conditioning is running, in the time period which contains time t, it can be known from equation (2) that its average operating power is P ac (t). Considering the natural loss of electric energy of the air conditioner virtual energy storage P ac− o (t), the virtual power of the air conditioning virtual energy storage is shown in the following equation: In equation (6), P ac− virtual (t) is the virtual power of the virtual energy storage that the air conditioning participates in. For equation (1), within a period of time, the temperature in the room where the air conditioner works will change to the set temperature, where τ is the length of the time period. en, equation (1) has Δt � τ.
After organization, After calculation with respect to equation (2), the relationship between the air conditioning power and the indoor and outdoor temperature as well as the set temperature and the length of the time period in the t th time period is obtained as follows: Equations (2), (5)- (7), and (10) are combined to obtain the relationship between the virtual charging power and temperature of the virtual energy storage as follows:

e Overall Structure of the Virtual Energy Storage System with the Participation of Electric Vehicles.
e schematic diagram of virtual energy storage with electric vehicle participation is shown in Figure 3. e total capacity of the electric vehicle battery is fixed. e battery capacity of each electric vehicle can be used for energy storage in addition to the daily required power. erefore, the amount of electricity that can be used for virtual energy storage of a single electric vehicle can be determined, and then the total virtual energy storage power of the electric vehicle cluster can be determined.

Important Parameters of Virtual Energy
en in the area i, the number of electric vehicles participating in the virtual energy storage in the time period t is shown in en for the remaining energy in the virtual energy storage at different times in the area i, equation (14) is as follows: In equation (14), S EV− virtual,i (t) is the virtual energy storage capacity of electric vehicles in the period of t in the i th region; W j,i (t) is the remaining power of the j th electric vehicle in the i th region in the time period t; W 0 j,i is the minimum electric power required for the j th electric car in the i th area, which does not change with time.
en, the rated energy storage of electric vehicles in the i th area is where S EV− virtual,i,N (t) is the rated energy storage capacity of the virtual energy storage in the time period t and W j,i,N (t) is the rated battery capacity of the j th electric vehicle in the area i in the time period t.

Virtual Charging and Discharging Power for Virtual
Energy Storage of Electric Vehicles. e charge and discharge power of the virtual energy storage of electric vehicles can be expressed by the following equation: In equation (16), P EV is the charging and discharging power of the electric vehicle charging and discharging piles, P EV− virtual,i (t) is the virtual charging power of the virtual energy storage of the electric vehicle in the area i with time period which includes time t, and NE i (t) is the number of electric vehicles participating in the virtual energy storage in the area i with time period which includes time t.
From the perspective of the electric vehicle's virtual energy storage capacity, its charging power meets From the perspective of the virtual energy storage capacity of electric vehicles, the discharge power meets In equations (17) and (18), τ is the length of the discharge time period.

Construction of the Mathematical Model of Joint Virtual Energy Storage Optimization Based on the Rolling Optimization Method
Aiming at the problem of determination of high-dimensional complex constraint decision variables in global optimization of electric vehicle virtual energy storage model, this paper proposes a feasible solution continuous rolling correction algorithm. Aiming at the global optimization problem of electric vehicle virtual energy storage output power, during the initialization process, due to the complexity of mutual coupling of local constraints and global constraints, the probability of obtaining a feasible solution is low, and initialization cannot be completed. In order to solve the above problems, this paper proposes constructing a standard function to measure the global constraint error that satisfies the local constraint decision variables during initialization and thus continuously revising the decision variables through continuous rolling optimization. Experiments show that this method can effectively improve the initialization efficiency of complex constrained optimization problems with high-dimensional decision variables and then obtain the global optimized output power of electric vehicle virtual energy storage, which verifies the economics of electric vehicle with the combination of virtual energy storage.

Rolling Optimization Process.
For the value of the virtual energy storage power in different time periods, the optimization result of the previous period will affect the virtual energy storage power range in the next period. Drawing on the idea of rolling optimization, the value of virtual energy storage power can be adjusted continuously with the advancement of the sampling time to determine the power of virtual energy storage in each period. e process is shown in Figure 4. e state in the time period t mainly includes the outdoor temperature and the number of arriving and leaving For the virtual energy storage in which air conditioning participate, For electric vehicles, for the virtual energy storage state at time t + 1 during the rolling optimization process in the region i, when the virtual energy storage state S EV− virtual,i (t), the input power P EV− virtual− in (t), and the output power P EV− virtual− out (t) at the previous time are known, the status of the virtual energy storage capacity of electric vehicles at the beginning of the next time period can be determined:

Objective Function.
In the case where other conditions are certain and uncontrollable, for the virtual energy storage optimization scheduling in which air conditioning and electric vehicles jointly participate, the optimization objective function has the lowest total cost of electricity consumption throughout the day: In equation (21), M is the total number of residential area and units in the area; in the area i at the time t, C WT,i (t) and C PV,i (t) are the total cost of wind power and photovoltaic power, respectively; C grid,i (t) is the cost of purchasing electricity from the power grid; C EV,i is the cost of virtual energy storage with electric vehicles participation in the area i throughout the day; C AC,i is the cost of air conditioning virtual energy storage throughout the day in the area i.
In addition, In equation (22), p WT,i,repare and p PV,i,repare are the power generation maintenance costs of wind turbines and photovoltaic power generation, respectively. And their units are CNY/kWh; E WT,i (t) and E PV,i (t) are the power generation capacity of wind turbines and photovoltaic power generation during the time period which contains time t: In equation (23), p grid (t) is the price of electric energy in the time period which contains time t, which can be step price, peak-valley price, time-sharing price, or real-time price, and E grid,i (t) represents time t: In equation (24), p EV is the cost of virtual energy storage with electric vehicles participation; the unit is CNY/kWh; E EV,i is the total amount of virtual energy storage with electric vehicles participation in the i th area in the whole day: In equation (25), p AC is the cost of in virtual energy storage with air conditioning participation; the unit is CNY/ kWh; E AC,i is the total amount of virtual energy storage with air conditioning participation in the i th area throughout the day. Complexity 5 equal to the electricity consumption at that time, as shown in equation (26): P ac− o is the natural power consumption load curve of the virtual energy storage of the air conditioning, P EV− o is the charging power load curve of the electric vehicle, demand is the daily load demand, P WT is the wind power, P PV is the photovoltaic power, P ac− virtual is the input (or output) power of the virtual energy storage with air conditioning participation, P grid is the power of the power grid, and P EV− virtual is the input (or output) power of the virtual energy storage involved in electric vehicles.

Output Power Constraints.
Any distributed power supply and virtual energy storage device have a limitation for power rating, as shown in equation (27): In this equation, P i, min (t) and P i, max (t) are the minimum and maximum values of the average output power of the i th device in the period of time t, respectively.

Restrictions on the Basic Driving Power Demand of the Next Day after the Electric Vehicle Participates in the Virtual Energy Storage.
e remaining power of any electric vehicle can at least meet the needs of traveling the next day, as shown in equation (28): W j,i is the remaining power of the j th electric car in the i th area, and W 0 j,i is the minimum power required for the j th electric vehicle in the i th area. e remaining power constraint for electric vehicles is participating in virtual energy storage. e remaining power of virtual energy storage cannot be negative: S EV− virtual,i is the remaining power of the i th electric vehicle.

Continuous Rolling Optimization of Feasible Solution Initialization for High-Dimensional Complex Constrained Optimization Problem
By analyzing the rolling optimization model proposed above, we can find that the model built in this chapter has the following characteristics: (1) e dimension of decision variables is high (2) ere are local constraints described by the rolling optimization model among decision variables (3) ere are global constraints on decision variables is section proposes a continuous rolling optimization algorithm for initializing feasible solutions for high-dimensional complex constrained optimization problems for rolling optimization models and complex constraints, so as to obtain decision variables that simultaneously satisfy multiple types of constraints.

Treatment Strategy for Multiconstraint Problems.
A large number of inequality constraints will result in multiple feasible regions with irregular shapes in the feasible solution space, and the decision variables must satisfy all the constraints, which requires the determination of the range of feasible regions. Take the feasible region of a two-dimensional optimization problem as an example, as shown in Figure 5.
In the figure shown above, A 1 and A 2 are, respectively, feasible regions of decision variables for two different constraints. Since the area of A 2 is small, processing the constraints in parallel will greatly reduce the probability of the decision variable falling in A 2 . is will make it difficult for a large number of decision variables to meet the A 2 constraint conditions, resulting in inefficient decision variable selection. However, when the problem to be solved is a high-dimensional optimization problem, the shape of the feasible region corresponding to the constraint conditions is more difficult to determine, and the intersection between the feasible regions is also more difficult to obtain. For the optimization problems shown in this chapter, the dimension of its solution space is shown in equation (30) as follows: In equation (30), N 1 and N 2 are the number of time periods divided per hour, which are both 4 in the questions in this chapter. Besides, T 1 and T 2 are the total time of the optimization cycle, which are both 24 in this paper, so for the optimization problem given in this chapter, its dimension is 192. In practical problems, as the influence factors considered increase, its dimensions may become larger.

Solution Feasible Probability.
For the determination of the feasible solution with multiple constraints, this section uses the solution feasible probability to evaluate its efficiency, as shown in equation (31): e following is an analysis of its value, first of all for the one-dimensional optimization problem. e feasible solution selection process is shown in Figure 6.
In the two one-dimensional constraints given in the above figure, if N feasible solutions are determined in parallel, the process is to first generate N feasible solutions for A 1 constraint in the [x 1 , x 2 ] interval and then judge whether it meets A 2 . According to the range covered by A 1 and A 2 , we can know the solution feasible probability when the feasible solutions are evenly distributed: 6 Complexity Claim: after choosing a feasible solution with A 1 constraint, the probability of satisfying A 2 constraint is P (A 1 ∩ A 2 )/(A 1 ) . e following dimension of the optimization problem is increased. In the two-dimensional optimization problem, the solution feasible probability is discussed. e simpler inequality constraint is selected; thus, the selection process of a feasible solution is shown in Figure 7. e feasible probability is shown in According to the solution feasible probability of the lowdimensional complex constrained optimization problem, if its P (A 1 ∩ A 2 )/(A 1 ) in each dimension is P, then the solution feasible probability of the N-dimensional multiconstrained optimization problem should be P N . When N is too large, the probability of its feasible solution selection probability will be very small, and it is easy to make the optimization process fall into an infinite loop. e following experiment determines the feasible probability of the solution of the optimization problem in different dimensions and sets P to 0.5. e results are shown in Table 1.
In view of the dimensional disaster in this paper, the undeterminable problem of feasible solutions in this chapter could be analyzed and solved. erefore, this section proposes a continuous rolling optimization algorithm to solve the inefficiency problem of multiconstrained feasible solutions.

Algorithm Flow.
For the high-dimensional complex constrained optimization problem, the analysis is first based on the two-dimensional optimization problem. When the optimization problem is faced with global constraints and local constraints, the optimization process is shown in Figure 8, where A 1 represents the constraint condition 1 and A 2 represents the constraint condition 2.
In the figure shown above, the red and blue paths are the optimization paths of the results with two consecutive A 1   optimizations. Since A 1 represents local constraints and x 1 has a mutual influence with x 2 , combined with the rolling optimization model proposed in the previous section, the evolutionary process of the continuous rolling optimization algorithm determined by the local constraint and the global constraint common feasible region of the optimization problem can be determined as shown in Figure 9. e specific steps are as follows: (1) Determine the value of the N-dimensional decision variable x 1 , . . ., x i .
Determination of X′ i+1 value based on rolling optimization Determination of X′ i+1 value based on rolling optimization

Optimization Effect Analysis.
In the problem solved in this chapter, one of the decision variables is the input and output power of the virtual energy storage with electric vehicle participation in each time period, and the remaining power of the electric vehicle virtual energy storage is used as the standard to detect whether the global optimization is completed. e distribution diagram of the arrival and departure time of electric vehicles in this simulation is shown in Figure 10. e results of the comparison between the unoptimized and optimized global constraint detection standard function values are shown in Figure 11.
As can be seen from Figure 10, the last two vehicles in the area left within the time period of 10:00-10:15 and the total power 47.78 kWh can be determined by calculation.
As can be seen from Figure 11, at 10:00 before rolling optimization, when the last two electric vehicles have not left, the remaining power of virtual energy storage is 154.4 kWh. According to simulation calculations, the remaining power of the two electric vehicles is 47.78 kWh.
at is, after the last two vehicles have left, the remaining energy of the virtual energy storage is still not 0, which obviously does not meet the global constraints, indicating that the decision variables fall outside the feasible domain of the global constraints, and the initialization fails. If the continuous rolling optimization algorithm is not adopted, the decision variables need to be initialized again.
After the continuous rolling optimization determined for the high-dimensional complex constraint optimization feasible solution, the results are consistent with the electric vehicle arrival and departure states, indicating that the continuous rolling optimization algorithm [5,24] for the high-dimensional optimization of global constraints and local constraints that jointly govern the feasible domain can effectively limit the decision variables within the global constraints and local constraints that jointly govern the feasible domain for the problems in this chapter, which helps improve the optimization efficiency.

Analysis
is section first simulates the virtual energy storage capacity of air conditioning and electric vehicles and illustrates the feasibility of virtual energy storage. en virtual energy storage is optimized for scheduling. In the example calculation process, the virtual energy storage power range of air conditioning and electric vehicles is first determined, and then the day is divided into 96 time periods to optimize the virtual energy storage power within the power range.

Analysis of Air Conditioning and Electric Vehicle Virtual
Energy Storage Capacity

Output Power of the Virtual Energy Storage Model with EV Participation.
In the virtual energy storage in which electric vehicles participate, the input and output power of the virtual energy storage and the arrival and departure of the electric vehicle will affect the remaining power of the virtual energy storage. Under rolling optimization, the relationship can be obtained by Figure 12.
As can be seen from Figure 12, the continuous energy optimization of the electric vehicle's virtual energy storage output meets the local constraints; the virtual energy storage remaining power becomes 0 after all the electric vehicles leave, indicating that the virtual energy storage output power meets the global constraint conditions.

Output Power of Virtual Energy Storage Model with Air
Conditioning Participation. In the virtual energy storage with air conditioning participation, the input and output power of the virtual energy storage and the natural power loss of the virtual energy storage based on the indoor and outdoor temperature difference will affect the remaining power of the virtual energy storage and the change of room temperature. Under rolling optimization, the changes of the above four parameters are shown in Figures 13 and 14.
It can be seen from the comparison that, within the temperature constraint range, the virtual energy storage of the air conditioner can realize the virtual power change of the virtual energy storage through the change of the air conditioner power.

Determination of the Range of Virtual Energy Storage
Power for Air Conditioning and Electric Vehicles

Virtual Energy Storage Charge and Discharge Power
Range with EV Participation. For the virtual energy storage of electric vehicles, according to the spatial and temporal distribution of electric vehicles given in Figure 12 and the power demand for driving on the day according to equations (17) and (18), the input and output power range of virtual energy storage with electric vehicles participation can be determined. e results are shown in Figure 15.

Determination of the Range of Charge and Discharge Power at the Starting Time in the Virtual Energy Storage with
Air Conditioning Participation. e daily temperature change curve in this area is shown in the red line in Figure 13. Considering the natural consumption rate of the virtual energy storage of the air conditioning, the output power of the virtual energy storage reaches the maximum when the air conditioner is not running, which is P ac− O (t); when the air conditioner is running, its input and output power are shown in equation (11). e maximum value of the virtual energy storage discharge power changes according to the change of the set temperature. e temperatures are set to 19.5°C, 20.5°C, and 21.5°C, respectively; Complexity in the following figure, the upper limit of the virtual energy storage output power changes as shown in Figure 16.
It can be seen from the above figure that the output power of the virtual energy storage decreases as temperature increases.

Economic Analysis of Optimization for Joint Virtual
Energy Storage. Based on air conditioning, electric vehicles have the ability to adjust the operating power within a certain range to convert electrical energy into thermal energy storage or perform bidirectional energy exchange with the power grid to achieve the virtual energy storage capacity of energy transfer. e following is the simulation of the effect of virtual energy storage of air conditioning and electric vehicles in this example. Regarding the load of the resident users in the constant temperature community, consider a total of about 650 small high-rise (five-story) households with a total construction area of 65,000 square meters, a total air conditioning power of 1000 kW, and an average of about 1.5 kW per household; on the top floor, there is about 6500 square meters of photovoltaic panels installed, and 10 wind turbines equipped; the power of its distributed power generation equipment is shown in Figure 17.
Considering that virtual energy storage with electric vehicles participation only charges 0.4883 CNY/kWh, so when actually using electric vehicles for virtual energy storage, the discharge cost of electric vehicles is 0.6883 CNY/ kWh, and the discharge efficiency is 80%. In addition, the maintenance cost for wind turbines and photovoltaic power generation is shown in Table 2 [8].     Complexity e total load demand curve of the residents is shown in Figure 18.
is environment corresponds to four scenes, respectively. Power (kW)

Scenario
The upper limit of virtual energy storage output power at 19.5 °C The upper limit of virtual energy storage output power at 20.5 °C The upper limit of virtual energy storage output power at 21.5 °C     e economics of these four scenarios are compared as follows.
In Scenario1, the results obtained by air conditioning and electric vehicles participating in virtual energy storage optimization scheduling and the results of using electric vehicles or air conditioning alone to participate in virtual energy storage are shown in Figure 19.
In Scenario2, the results obtained by air conditioners independently participating in virtual energy storage optimization scheduling are shown in Figure 20.
In Scenario3, the results of electric vehicles independently participating in virtual energy storage optimization scheduling are shown in Figure 21.
In the above four scenarios, the comparison results of the total daily electricity charges in the residential area are shown in Table 3.
From the above comparative analysis, the following conclusions can be drawn: (1) When the joint virtual energy storage of air conditioners and electric vehicles is adopted, the total cost of electricity charges in residential areas can be reduced by 10%. (2) e photovoltaic power generation in residential areas is highly compatible with the operating time of air conditioners. erefore, when participating in virtual energy storage alone, the effect of air conditioning participating in virtual energy storage is better than electric vehicles participating in virtual energy storage.
(3) Analyze the reasons for the different total electricity costs in residential areas under different scenarios on a certain day, mainly because air conditioning and electric vehicles can store the electrical energy generated by new energy power generation equipment, which reduces the phenomenon of wind and heat rejection, realizes the efficient utilization of new energy power generation equipment, and also reduces the cost of purchasing electricity from the power grid. erefore, it has also proved that the joint

Summary
is paper proposes the use of air conditioning and electric vehicles to jointly participate in virtual energy storage to realize the economic dispatch of energy local area Smart Grid in view of the current status of the controllable load of air conditioning and the load growth of electric vehicles.
Firstly, the virtual energy storage model of thermoelectric conversion with the participation of air conditioning load and the electric energy transfer virtual energy storage model with electric vehicle participation are established. Besides, the rolling optimization idea is introduced to restrict the input and output power of the virtual energy storage with air conditioning and electric vehicle participation.
Secondly, a joint virtual energy storage power optimization model is constructed, and for the solution of the model, in the initialization process of decision variables, due to the condition that dimensional disaster occurred and the initial feasible solution cannot be obtained, a continuous rolling optimization method is proposed for the feasible solution of the high-dimensional complex constraint optimization problem, so that the virtual energy storage power can simultaneously meet the local and global constraints under rolling optimization during the initialization process and thus improve the initialization efficiency.
Finally, the capacity and economy of joint virtual energy storage in residential areas are calculated. e results prove that air conditioning and electric vehicles have the ability to jointly participate in virtual energy storage, and the comparison proves that joint virtual energy storage can effectively improve the economics of electricity consumption. erefore, the conclusion is that the joint virtual energy storage with the participation of electric vehicles and air conditioning helps improve the economics of consumer electricity consumption.
Data Availability e data used in this paper comes from the statistics of urban travel data, which has not been published publicly. is paper only takes this data as an example to prove the effectiveness of the proposed method.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this article.

Authors' Contributions
Rui-Cheng Dai and Xiao-di Zhang participated in the algorithm simulation and draft writing. Bi Zhao and Xiao-di Zhang participated in the concept, design, interpretation, and commenting on the manuscript, and the critical revision   Figure 20: Output power of virtual energy storage with independent participation of air conditioning in Scenario2.
14 Complexity of this paper. Jun-Wei Yu, Bo Fan, and Biao Liu participated in the data collection and analysis of the paper.