TY - JOUR AU1 - Hua, Yuanpeng AU2 - Wang, Shiqian AU3 - Wang, Yuanyuan AU4 - Zhang, Linru AU5 - Liu, Weiliang AB - Introduction The influx of a vast number of EVs accessing the grid may inevitably lead to a challenging situation of peak demand, putting a strain on grid stability. However, EVs can function as an effective energy storage device, and through the optimization of vehicle networking scheduling strategy, vehicles can help the power grid to achieve peak regulation, frequency modulation and new energy consumption. Specifically, by guiding EVs to supply electricity to the grid during peak periods and charge during low peak periods, the stability and safety of grid operation can be improved [1]. Incentive based on electricity price is a common method for electric vehicles to participate in dispatching. Wei Wu [2], studied the economic value of electric vehicles connected to the grid, established a power supply cost model, and investigated the cost reduction under three charging operation modes: random charging, controlled charging, and vehicle-to-grid (V2G) charging. Literature [3, 4], based on the previous questionnaire results, the behavior reality of electric vehicle users and empirical deduction, an electric vehicle behavior probability model was established to describe the randomness and uncertainty behaviors of EV users in travel and behavior decision-making. A multi-objective optimal scheduling model is established with the economic benefits of EV users and the load side peak load and valley filling as objective functions. Pan Xiaotian [5] established a four-objective optimization scheduling model with the aim of stabilizing load fluctuations, minimizing user costs, and extending elastic travel time and charging state to the greatest extent. With this model, user needs can be better met by providing a sufficiently flexible state of charge with higher travel time. Literature [6–8] respectively considered the problem of electric buses participating in power grid dispatching and the problem of resource allocation of electric vehicle charging piles. Compared with ordinary private cars, electric buses have the characteristics of limited operating range, long charging time, complex grid characteristics, etc. By introducing multi-objective multi depot optimal scheduling model, the total operating cost and peak load generated by concurrent charging activities can be minimized. The uncertainty of EV user preferences and decisions may affect V2G scheduling, resulting in the imbalance between the electric vehicle’s schedulable capacity and the required power. However, the charging pile resource allocation method proposed in this paper based on the two-stage classification and hierarchical scheduling framework can solve such problems in real time. Reference [9] investigated the participation of electric vehicles in the energy scheduling of virtual power plants. When the electric vehicle aggregator adopts the deterministic strategy and the virtual power plant adopts the stochastic strategy, the energy complementarity is realized and the overall operating economy is improved. Ju L [10] incorporated electric vehicles into carbon virtual power plants as a flexible resource and used the concept of electric vehicle aggregators to flexibly respond to grid operation requirements. Sheng [11],proposed a multi-time scale active distribution network (ADN) scheduling method, which includes backup coordination strategy and scheduling framework to improve the adaptive capacity of distribution network and reduce the impact of fluctuating power on the upstream transmission network. The backup coordination strategy can schedule available backup resources based on their temporal and spatial characteristics. Literature [12, 13] puts forward reactive power optimization strategies for power grid including electric vehicles, and establishes a reactive power optimization model aiming at reducing voltage deviation and network loss, so as to reduce the operating pressure of traditional reactive power compensation equipment. As the link between the electric grid and electric vehicle users, electric vehicle aggregators play a crucial role in coordinating the economic interests between the grid and users, and are indispensable participants in the interaction process of the vehicle network. In literature [14–16], aggregators focus on maximizing their own economic benefits, while also taking into account the demand response needs of users and the power grid, and they will participate in V2G as an auxiliary regulation method. In addition, the analysis of the scheduling potential of electric vehicles is the basis of the implementation of vehicle network interactive optimization scheduling technology. Literature [4, 17, 18], has established a probability model for the spatiotemporal characteristics of electric vehicles considering various travel needs of users. At present, there are limited studies considering uncertainty factors in V2G technology. Kong [19] proposed a bid-based double-layer multi-time scale scheduling method for multi-operator virtual power plants In the upper layer, a bidding equilibrium-based power allocation and internal pricing method for operators was proposed. The fluctuation cost coefficient was introduced to express the impact of the uncertainty of renewable energy generation on the bidding process. The above research results are true and effective, and provided a good reference for the research in this paper. The above paper has conducted a detailed study on the potential of electric vehicle scheduling, fully considering the impact of electric vehicle participation in grid scheduling on the reliability and economy of grid operation, and has achieved good research results. However, the research in the above papers is generally one-sided and does not fully utilize the scheduling potential of electric vehicles, and the final scheduling results do not meet the actual user needs, making practical application difficult. We have developed an EV model in this study that can effectively and reasonably describe the charging and discharging characteristics and behavior of EVs, fully considering the goals and constraints that EVs need to achieve in participating in the power grid dispatch process, introducing uncertainty description, and establishing an interval multi-objective optimization model to make the obtained optimal dispatch scheme more reasonable and practical. Methods Model for temporal and spatial distribution characteristics of electric vehicle load The model has been modified based on the probability distribution fitting of electric vehicle load time characteristics and spatial characteristics [4] using electric vehicle daily travel data from a certain district in Chongqing. The resulting model is as follows. Parking time. Based on the statistical data, the approximation of parking duration is estimated to follow a log-normal distribution [4, 17], with a mean of μp = 2.8 and a variance of σp = 0.95. The probability density function is given by the following equation: (1) Arrival/departure time. The start time of charging or discharging for electric vehicles at a charging station can be approximated as the entry time, which follows a normal distribution with a mean of μs = 17.5 and a variance of σs = 4.1. The probability density function of electric vehicle arrival time is defined by the given equation. (2) Electric vehicle departure time from the station is given by the following equation: (3) Daily driving mileage. The initial State of Charge (SOC) of an electric vehicle is determined by its daily driving mileage. By fitting data, it has been found that the daily driving mileage of electric vehicles follows a log-normal distribution with a mean value of μr = 3.4 and a variance of σr = 0.98. The probability density function for this distribution is given as follows: (4) Planned charging power. In order to ensure that the EVs have the desired SOC by the expected departure time, it is necessary to calculate the planned charging power Pevplan at time t, which can only participate in scheduling if it falls within the upper and lower limits of the EV charging power. (5) In the equation, SOCq represents the expected SOC of the user when leaving the station, represents the current SOC, Ce represents the rated capacity, tq represents the user’s expected departure time, and td represents the current time. Initial SOC. (6) In the equation, SOCinit represents the initial SOC of the electric vehicle upon arrival at the charging station, dr represents the daily travel distance of the electric vehicle, and de represents the rated travel distance of the electric vehicle. Analysis of dispatch potential for electric vehicles Charging behavior of a single electric vehicle follows the parallelogram law [4, 17]. The charging or discharging power for each time period can be determined based on the actual scheduling needs between the upper and lower limits of scheduling, and the maximum waiting time is a measure of whether the electric vehicle can accept discharging scheduling at time t (i.e., the electric vehicle must start charging when this time is reached). In the Fig 1, Cev represents the scheduling capacity. (7) In the expression, Cevq represents the desired electric quantity, and Cevs represents the current electric quantity. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 1. Schematic diagram of EV scheduling potential. https://doi.org/10.1371/journal.pone.0297855.g001 The scheduling potential of electric vehicles at time t in a large-scale electric vehicle system refers to the total potential of dispatchable charging and discharging load for participating electric vehicles in the region. Based on the current total power value at time t, the upper and lower limits of dispatchable charging and discharging load at time t+1 can be predicted. (8)(9)(10) In the expression, represents the upper limit of scheduling, represents the lower limit of scheduling, represents the charging or discharging power of electric vehicles, represents the additional EV load that must be added at time t+1 for scheduling, represents the EV load added only for charging scheduling, and represents the EV load reduced only for discharging scheduling. Multi-objective optimization scheduling model Parameters for optimization. For a given scheduling period, let xs denote the decision variables for all scheduling time slots θT in the model. For time slot t, the decision variables xs are shown in Table 1. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 1. Decision variables. https://doi.org/10.1371/journal.pone.0297855.t001 Objective function. Objective function for economic objective. The economic objective [19–21] function primarily considers the operational costs Cag of the regional power grid, including the cost of exchanging power with the external grid Cgrid, the cost of wind power generation Cw, the cost of photovoltaic power generation Cpv, the cost of thermal power generation Cf, and the cost of energy storage system charging and discharging Cb. (11) (12) (13) (14) (15) (16) In Eqs 12–16, represents the exchange power price between the regional power grid and external power grid, is the generation cost per unit of wind power generation, is the generation cost per unit of photovoltaic power generation, is the generation cost per unit of thermal power generation, is the cost per unit of energy storage device for charging or discharging,. Objective function for peak shaving and valley filling. (17) In the equation, Vload represents the variance of the load curve, represents the power of in-network electricity consumption, and represents the average load of electricity consumption. Objective function for new energy integration. Due to the limitations of transmission lines and substations, there are constraints on the integration of new energy. The expression for the objective function for new energy integration is as follows: (18) In the equation, ΔPx represents the curtailed wind and photovoltaic power generation at a certain time,represents the maximum power generation of the photovoltaic units at time t, represents the maximum power generation of the wind power units at time t. Constraints. For ∀t∈θT, the variables should satisfy the following constraints. 1) Constraints on power balance [21]. (19) In the equation, represents the conventional electricity load at time t. 2) Constraints on electric vehicle dispatch potential. (20) In the equation, represents the lower limit power of electric vehicle charging or discharging scheduling during time period t, and represents the upper limit power of electric vehicle charging or discharging scheduling during time period t. 3) Constraints on upper and lower power limit for each device [21]. (21) In the equation, Pgrid_min represents the lower limit of power exchanged with the external power grid, Pgrid_max represents the upper limit of power exchanged with the external power grid, Pev_min represents the lower limit of electric vehicle charging and discharging power, Pev_max represents the upper limit of electric vehicle charging and discharging power, Pb_min represents the lower limit of battery charging and discharging power, Pb_max represents the upper limit of battery charging and discharging power, Pf_min represents the lower limit of thermal power generation, Pf_max represents the upper limit of thermal power generation, Pw_max represents the upper limit of wind power generation, and Ppv_max represents the upper limit of photovoltaic power generation. 4) Ramp rate constraint [22] for thermal power generation units. (22) In the equation, Pdown_max represents the maximum ramp rate constraint for thermal power generation units in the downward direction, and Pup_max represents the maximum ramp rate constraint for thermal power generation units in the upward direction. 5) State of charge constraint [20, 21] for energy storage devices (23) In the equation, SOCT represents the final value of SOC for energy storage devices, SOCinit represents the initial value of SOC for energy storage devices, SOCb_min represents the lower limit of SOC for energy storage devices, and SOCb_max represents the upper limit of SOC for energy storage devices. Interval multi-objective optimization scheduling model Interval optimization theory. The multi-objective optimization [23, 24] problem with intervals can be explicitly expressed as follows. (24) In the equation, x is the decision variable, Ω is the decision space, c is the interval vector, gj(x,c)≥qj represents the interval inequality constraint, and hk(x,c) = rk represents the interval equality constraint, is one of the objective functions, as it contains an interval vector c, and its value is also an interval number. Description of uncertainty factors. The uncertainty and random fluctuation range of wind power generation [25, 26], photovoltaic power generation [26, 27], and electrical load power are described by the interval numbers [Pw], [Ppv], and [PImload], respectively. Through probabilistic analysis, the aforementioned uncertain factors randomly fluctuate around their predicted values, and these fluctuations are symmetrical. (25)(26)(27) Outer optimization model. The outer optimization model [28, 29] is usually used to determine the range of decision variables, with the goal of finding the optimal range of decision variables that enables the inner optimization model to optimize within this range and obtain the optimal solution. The objective function can be represented as follows. (28) The constraints of the outer optimization model are basically the same as those in the ’Constraints’ section, and the constraints involving the number of intervals are as follows, for ∀t ∈θT: 1) Power balance. (29) 2) Power constraints for each device. (30) Inner optimization model. The inner optimization model [28, 29] optimizes within the range of decision variables x* s determined by the given outer optimization model to achieve optimization objectives. (31)(32)(33) Model solving. This paper uses an improved NSGA-II algorithm to solve the problem of outer optimization model solving. By introducing interval credibility and interval overlap, it judges the situation where individuals meet the constraint conditions and compares individuals with the same rank to calculate individual crowding distance. For the interval numbers [30, 31] q1 = [a1,b1], q2 = [a2,b2] ∈ I(R), the expression of interval credibility is as follows, where w1 = w(q1) and w2 = w(q2): (34) Eq (34) is recorded as the interval confidence level of q1≥q2. For an individual x, if its interval confidence level satisfying the j-th constraint is ξj, then: (35) Correspondingly, the degree of violation Lj of the j-th constraint for individual x is: (36) By setting the confidence threshold ξ* j, an individual x is considered feasible if its interval confidence level for a certain constraint is greater than or equal to this threshold, otherwise, it is considered infeasible. The dominance relation is determined by comparing the dominance relationships between feasible solutions and infeasible solutions. For evolutionary individuals x1 and x2 with the same rank value, the intersection of their objective function values can be represented by fi(x1,c)∩fi (x2,c). The interval overlap degree between x1 and x2 is expressed as: (37) Model for temporal and spatial distribution characteristics of electric vehicle load The model has been modified based on the probability distribution fitting of electric vehicle load time characteristics and spatial characteristics [4] using electric vehicle daily travel data from a certain district in Chongqing. The resulting model is as follows. Parking time. Based on the statistical data, the approximation of parking duration is estimated to follow a log-normal distribution [4, 17], with a mean of μp = 2.8 and a variance of σp = 0.95. The probability density function is given by the following equation: (1) Arrival/departure time. The start time of charging or discharging for electric vehicles at a charging station can be approximated as the entry time, which follows a normal distribution with a mean of μs = 17.5 and a variance of σs = 4.1. The probability density function of electric vehicle arrival time is defined by the given equation. (2) Electric vehicle departure time from the station is given by the following equation: (3) Daily driving mileage. The initial State of Charge (SOC) of an electric vehicle is determined by its daily driving mileage. By fitting data, it has been found that the daily driving mileage of electric vehicles follows a log-normal distribution with a mean value of μr = 3.4 and a variance of σr = 0.98. The probability density function for this distribution is given as follows: (4) Planned charging power. In order to ensure that the EVs have the desired SOC by the expected departure time, it is necessary to calculate the planned charging power Pevplan at time t, which can only participate in scheduling if it falls within the upper and lower limits of the EV charging power. (5) In the equation, SOCq represents the expected SOC of the user when leaving the station, represents the current SOC, Ce represents the rated capacity, tq represents the user’s expected departure time, and td represents the current time. Initial SOC. (6) In the equation, SOCinit represents the initial SOC of the electric vehicle upon arrival at the charging station, dr represents the daily travel distance of the electric vehicle, and de represents the rated travel distance of the electric vehicle. Analysis of dispatch potential for electric vehicles Charging behavior of a single electric vehicle follows the parallelogram law [4, 17]. The charging or discharging power for each time period can be determined based on the actual scheduling needs between the upper and lower limits of scheduling, and the maximum waiting time is a measure of whether the electric vehicle can accept discharging scheduling at time t (i.e., the electric vehicle must start charging when this time is reached). In the Fig 1, Cev represents the scheduling capacity. (7) In the expression, Cevq represents the desired electric quantity, and Cevs represents the current electric quantity. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 1. Schematic diagram of EV scheduling potential. https://doi.org/10.1371/journal.pone.0297855.g001 The scheduling potential of electric vehicles at time t in a large-scale electric vehicle system refers to the total potential of dispatchable charging and discharging load for participating electric vehicles in the region. Based on the current total power value at time t, the upper and lower limits of dispatchable charging and discharging load at time t+1 can be predicted. (8)(9)(10) In the expression, represents the upper limit of scheduling, represents the lower limit of scheduling, represents the charging or discharging power of electric vehicles, represents the additional EV load that must be added at time t+1 for scheduling, represents the EV load added only for charging scheduling, and represents the EV load reduced only for discharging scheduling. Multi-objective optimization scheduling model Parameters for optimization. For a given scheduling period, let xs denote the decision variables for all scheduling time slots θT in the model. For time slot t, the decision variables xs are shown in Table 1. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 1. Decision variables. https://doi.org/10.1371/journal.pone.0297855.t001 Objective function. Objective function for economic objective. The economic objective [19–21] function primarily considers the operational costs Cag of the regional power grid, including the cost of exchanging power with the external grid Cgrid, the cost of wind power generation Cw, the cost of photovoltaic power generation Cpv, the cost of thermal power generation Cf, and the cost of energy storage system charging and discharging Cb. (11) (12) (13) (14) (15) (16) In Eqs 12–16, represents the exchange power price between the regional power grid and external power grid, is the generation cost per unit of wind power generation, is the generation cost per unit of photovoltaic power generation, is the generation cost per unit of thermal power generation, is the cost per unit of energy storage device for charging or discharging,. Objective function for peak shaving and valley filling. (17) In the equation, Vload represents the variance of the load curve, represents the power of in-network electricity consumption, and represents the average load of electricity consumption. Objective function for new energy integration. Due to the limitations of transmission lines and substations, there are constraints on the integration of new energy. The expression for the objective function for new energy integration is as follows: (18) In the equation, ΔPx represents the curtailed wind and photovoltaic power generation at a certain time,represents the maximum power generation of the photovoltaic units at time t, represents the maximum power generation of the wind power units at time t. Constraints. For ∀t∈θT, the variables should satisfy the following constraints. 1) Constraints on power balance [21]. (19) In the equation, represents the conventional electricity load at time t. 2) Constraints on electric vehicle dispatch potential. (20) In the equation, represents the lower limit power of electric vehicle charging or discharging scheduling during time period t, and represents the upper limit power of electric vehicle charging or discharging scheduling during time period t. 3) Constraints on upper and lower power limit for each device [21]. (21) In the equation, Pgrid_min represents the lower limit of power exchanged with the external power grid, Pgrid_max represents the upper limit of power exchanged with the external power grid, Pev_min represents the lower limit of electric vehicle charging and discharging power, Pev_max represents the upper limit of electric vehicle charging and discharging power, Pb_min represents the lower limit of battery charging and discharging power, Pb_max represents the upper limit of battery charging and discharging power, Pf_min represents the lower limit of thermal power generation, Pf_max represents the upper limit of thermal power generation, Pw_max represents the upper limit of wind power generation, and Ppv_max represents the upper limit of photovoltaic power generation. 4) Ramp rate constraint [22] for thermal power generation units. (22) In the equation, Pdown_max represents the maximum ramp rate constraint for thermal power generation units in the downward direction, and Pup_max represents the maximum ramp rate constraint for thermal power generation units in the upward direction. 5) State of charge constraint [20, 21] for energy storage devices (23) In the equation, SOCT represents the final value of SOC for energy storage devices, SOCinit represents the initial value of SOC for energy storage devices, SOCb_min represents the lower limit of SOC for energy storage devices, and SOCb_max represents the upper limit of SOC for energy storage devices. Parameters for optimization. For a given scheduling period, let xs denote the decision variables for all scheduling time slots θT in the model. For time slot t, the decision variables xs are shown in Table 1. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 1. Decision variables. https://doi.org/10.1371/journal.pone.0297855.t001 Objective function. Objective function for economic objective. The economic objective [19–21] function primarily considers the operational costs Cag of the regional power grid, including the cost of exchanging power with the external grid Cgrid, the cost of wind power generation Cw, the cost of photovoltaic power generation Cpv, the cost of thermal power generation Cf, and the cost of energy storage system charging and discharging Cb. (11) (12) (13) (14) (15) (16) In Eqs 12–16, represents the exchange power price between the regional power grid and external power grid, is the generation cost per unit of wind power generation, is the generation cost per unit of photovoltaic power generation, is the generation cost per unit of thermal power generation, is the cost per unit of energy storage device for charging or discharging,. Objective function for peak shaving and valley filling. (17) In the equation, Vload represents the variance of the load curve, represents the power of in-network electricity consumption, and represents the average load of electricity consumption. Objective function for new energy integration. Due to the limitations of transmission lines and substations, there are constraints on the integration of new energy. The expression for the objective function for new energy integration is as follows: (18) In the equation, ΔPx represents the curtailed wind and photovoltaic power generation at a certain time,represents the maximum power generation of the photovoltaic units at time t, represents the maximum power generation of the wind power units at time t. Constraints. For ∀t∈θT, the variables should satisfy the following constraints. 1) Constraints on power balance [21]. (19) In the equation, represents the conventional electricity load at time t. 2) Constraints on electric vehicle dispatch potential. (20) In the equation, represents the lower limit power of electric vehicle charging or discharging scheduling during time period t, and represents the upper limit power of electric vehicle charging or discharging scheduling during time period t. 3) Constraints on upper and lower power limit for each device [21]. (21) In the equation, Pgrid_min represents the lower limit of power exchanged with the external power grid, Pgrid_max represents the upper limit of power exchanged with the external power grid, Pev_min represents the lower limit of electric vehicle charging and discharging power, Pev_max represents the upper limit of electric vehicle charging and discharging power, Pb_min represents the lower limit of battery charging and discharging power, Pb_max represents the upper limit of battery charging and discharging power, Pf_min represents the lower limit of thermal power generation, Pf_max represents the upper limit of thermal power generation, Pw_max represents the upper limit of wind power generation, and Ppv_max represents the upper limit of photovoltaic power generation. 4) Ramp rate constraint [22] for thermal power generation units. (22) In the equation, Pdown_max represents the maximum ramp rate constraint for thermal power generation units in the downward direction, and Pup_max represents the maximum ramp rate constraint for thermal power generation units in the upward direction. 5) State of charge constraint [20, 21] for energy storage devices (23) In the equation, SOCT represents the final value of SOC for energy storage devices, SOCinit represents the initial value of SOC for energy storage devices, SOCb_min represents the lower limit of SOC for energy storage devices, and SOCb_max represents the upper limit of SOC for energy storage devices. Interval multi-objective optimization scheduling model Interval optimization theory. The multi-objective optimization [23, 24] problem with intervals can be explicitly expressed as follows. (24) In the equation, x is the decision variable, Ω is the decision space, c is the interval vector, gj(x,c)≥qj represents the interval inequality constraint, and hk(x,c) = rk represents the interval equality constraint, is one of the objective functions, as it contains an interval vector c, and its value is also an interval number. Description of uncertainty factors. The uncertainty and random fluctuation range of wind power generation [25, 26], photovoltaic power generation [26, 27], and electrical load power are described by the interval numbers [Pw], [Ppv], and [PImload], respectively. Through probabilistic analysis, the aforementioned uncertain factors randomly fluctuate around their predicted values, and these fluctuations are symmetrical. (25)(26)(27) Outer optimization model. The outer optimization model [28, 29] is usually used to determine the range of decision variables, with the goal of finding the optimal range of decision variables that enables the inner optimization model to optimize within this range and obtain the optimal solution. The objective function can be represented as follows. (28) The constraints of the outer optimization model are basically the same as those in the ’Constraints’ section, and the constraints involving the number of intervals are as follows, for ∀t ∈θT: 1) Power balance. (29) 2) Power constraints for each device. (30) Inner optimization model. The inner optimization model [28, 29] optimizes within the range of decision variables x* s determined by the given outer optimization model to achieve optimization objectives. (31)(32)(33) Model solving. This paper uses an improved NSGA-II algorithm to solve the problem of outer optimization model solving. By introducing interval credibility and interval overlap, it judges the situation where individuals meet the constraint conditions and compares individuals with the same rank to calculate individual crowding distance. For the interval numbers [30, 31] q1 = [a1,b1], q2 = [a2,b2] ∈ I(R), the expression of interval credibility is as follows, where w1 = w(q1) and w2 = w(q2): (34) Eq (34) is recorded as the interval confidence level of q1≥q2. For an individual x, if its interval confidence level satisfying the j-th constraint is ξj, then: (35) Correspondingly, the degree of violation Lj of the j-th constraint for individual x is: (36) By setting the confidence threshold ξ* j, an individual x is considered feasible if its interval confidence level for a certain constraint is greater than or equal to this threshold, otherwise, it is considered infeasible. The dominance relation is determined by comparing the dominance relationships between feasible solutions and infeasible solutions. For evolutionary individuals x1 and x2 with the same rank value, the intersection of their objective function values can be represented by fi(x1,c)∩fi (x2,c). The interval overlap degree between x1 and x2 is expressed as: (37) Interval optimization theory. The multi-objective optimization [23, 24] problem with intervals can be explicitly expressed as follows. (24) In the equation, x is the decision variable, Ω is the decision space, c is the interval vector, gj(x,c)≥qj represents the interval inequality constraint, and hk(x,c) = rk represents the interval equality constraint, is one of the objective functions, as it contains an interval vector c, and its value is also an interval number. Description of uncertainty factors. The uncertainty and random fluctuation range of wind power generation [25, 26], photovoltaic power generation [26, 27], and electrical load power are described by the interval numbers [Pw], [Ppv], and [PImload], respectively. Through probabilistic analysis, the aforementioned uncertain factors randomly fluctuate around their predicted values, and these fluctuations are symmetrical. (25)(26)(27) Outer optimization model. The outer optimization model [28, 29] is usually used to determine the range of decision variables, with the goal of finding the optimal range of decision variables that enables the inner optimization model to optimize within this range and obtain the optimal solution. The objective function can be represented as follows. (28) The constraints of the outer optimization model are basically the same as those in the ’Constraints’ section, and the constraints involving the number of intervals are as follows, for ∀t ∈θT: 1) Power balance. (29) 2) Power constraints for each device. (30) Inner optimization model. The inner optimization model [28, 29] optimizes within the range of decision variables x* s determined by the given outer optimization model to achieve optimization objectives. (31)(32)(33) Model solving. This paper uses an improved NSGA-II algorithm to solve the problem of outer optimization model solving. By introducing interval credibility and interval overlap, it judges the situation where individuals meet the constraint conditions and compares individuals with the same rank to calculate individual crowding distance. For the interval numbers [30, 31] q1 = [a1,b1], q2 = [a2,b2] ∈ I(R), the expression of interval credibility is as follows, where w1 = w(q1) and w2 = w(q2): (34) Eq (34) is recorded as the interval confidence level of q1≥q2. For an individual x, if its interval confidence level satisfying the j-th constraint is ξj, then: (35) Correspondingly, the degree of violation Lj of the j-th constraint for individual x is: (36) By setting the confidence threshold ξ* j, an individual x is considered feasible if its interval confidence level for a certain constraint is greater than or equal to this threshold, otherwise, it is considered infeasible. The dominance relation is determined by comparing the dominance relationships between feasible solutions and infeasible solutions. For evolutionary individuals x1 and x2 with the same rank value, the intersection of their objective function values can be represented by fi(x1,c)∩fi (x2,c). The interval overlap degree between x1 and x2 is expressed as: (37) Results and discussion The regional power grid in a certain urban area of a city was selected as the research object, and a case study was conducted on the energy optimization and scheduling of the regional power grid. The 224-hour day was divided into 24 scheduling periods, with each hour serving as a scheduling period. In the case study, the wind power unit in the regional power grid was configured as 30 MW, the photovoltaic unit was 15 MW, the thermal power unit was 65 MW, the energy storage equipment was 10 MWh, and approximately 10,000 electric vehicles participated in the scheduling. The exchange power limit of the regional power grid’s purchase and sale of electricity to the external power grid is ±60 MW. The NSGA-II algorithm is set with the population size of Npop = 200, the maximum number of iterations of gmax = 10000, the crossover probability of Pc = 0.9, and the mutation probability of Pm = 0.2. EVs participating in the scheduling In the Pareto frontier, the distribution of objective function values is shown in Fig 2, where each point represents a Pareto optimal solution. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 2. The pareto frontier with the participation of electric vehicles. https://doi.org/10.1371/journal.pone.0297855.g002 Three optimization and scheduling schemes were selected for analysis based on their economic viability, high renewable energy absorption capacity, and effective peak-shaving and valley-filling effects. The most user-demand-oriented optimization and scheduling scheme was determined using the AHP decision-making analysis method. The objective function values under different scheduling schemes are shown in Table 2. Clearly, due to the lower cost of renewable energy generation compared to thermal power generation, the economic index of the power grid and the peak-shaving and valley-filling indexes and renewable energy absorption capacity index exhibit a trend of trade-off. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 2. The objective function values of scheduling schemes considering extreme scenarios. https://doi.org/10.1371/journal.pone.0297855.t002 Plot the three optimization and scheduling schemes that are the most economically viable, have the best peak-shaving and valley-filling effects, and have the highest renewable energy absorption capacity. 1) The best economic dispatch results. The power curves of the most economical scheduling scheme are shown in Figs 3 and 4. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 3. New energy generation power and EV power. https://doi.org/10.1371/journal.pone.0297855.g003 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 4. Thermal power generation and load power curve. https://doi.org/10.1371/journal.pone.0297855.g004 In the most economically efficient scheduling plan, during periods of high electricity prices and peak usage, the regional power grid needs to sell electricity to the large grid to meet the demand, resulting in an increase in power generation from wind, solar, and thermal sources. However, due to the limitations of power transmission capacity, some of the power generated from wind and solar sources may be wasted, resulting in a waste of energy resources. Conversely, during periods of low electricity prices and low usage, the regional power grid needs to buy more electricity from the large grid, resulting in a decrease in power generation from wind, solar, and thermal sources. As energy storage devices have lower participation costs in scheduling, they are used more during this time period, while the participation power of EVs decreases. This may cause a decrease in the system’s ability to balance peak and off-peak demand, as the scheduling capacity of energy storage devices is limited and may not be able to meet the demands during peak periods. 2) The results of the optimal peak-shaving and valley-filling scheme. The power curves of the optimal scheduling scheme for peak shaving and valley filling are shown in Figs 5 and 6. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 5. New energy generation power and EV power. https://doi.org/10.1371/journal.pone.0297855.g005 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 6. Thermal power generation and load power curve. https://doi.org/10.1371/journal.pone.0297855.g006 For the optimal peak-shaving and valley-filling scheme, in order to make the load curve as smooth as possible and reduce the peak-to-valley difference, EVs and energy storage devices are charged at high power during periods of low electricity demand and discharged at high power during periods of high electricity demand. As both thermal power generation and energy storage systems can control their power output, they can adjust their generation capacity flexibly to provide additional electricity during peak load periods and reduce the load during off-peak periods. However, this will increase the involvement of thermal power generation in load regulation, which will lead to a decrease in economic goals and a lower capacity for integrating new energy sources. 3) The results of the scheme with the highest capacity for accommodating new energy consumption. The power curves of the scheduling scheme with the highest new energy consumption capacity are shown in Figs 7 and 8. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 7. New energy generation power and EV power. https://doi.org/10.1371/journal.pone.0297855.g007 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 8. Thermal power generation and load power curve. https://doi.org/10.1371/journal.pone.0297855.g008 For the optimal strategy aimed at the highest new energy consumption capacity, the best approach is to maintain relatively high power generation of wind and photovoltaic power units at each scheduling time. However, since the output of wind and photovoltaic power generation fluctuates due to environmental factors, the peak-shaving and valley-filling effect of this approach is not ideal. In contrast, thermal power generation with higher output during the low valley of wind and photovoltaic power generation can compensate for the insufficient output of wind and photovoltaic power generation, while reducing output during the peak of wind power generation to avoid wasting resources. In addition, the power output of EVs participating in scheduling should have higher discharge power and lower charging power during the low valley of wind and photovoltaic power generation, and lower discharge power and higher charging power during the peak of wind and photovoltaic power generation. Based on the importance of three decision factors, a pairwise comparison judgment matrix A is constructed. After consistency check and hierarchical ranking, the normalized weight values of each index are W = [0.1047 0.6370 0.2583]. The optimal scheduling scheme is selected. The variable values of the optimal scheduling scheme are shown in Figs 9 and 10: Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 9. Power output of thermal power units and load power. https://doi.org/10.1371/journal.pone.0297855.g009 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 10. Grid connected power of new energy generation. https://doi.org/10.1371/journal.pone.0297855.g010 EVs do not participate in scheduling When electric vehicles do not participate in scheduling, we believe that they start charging when they enter the station, until the charging is completed or the electric vehicle leaves the station. Establishing a regional power grid scheduling model based on electric vehicles as regular loads that do not participate in scheduling, with the goal of power grid economics and new energy consumption as the objective functions. The decision variables include thermal power generation, photovoltaic power generation, and wind power generation, with the optimal Pareto solution set and frontier obtained. The optimal dispatching scheme is obtained based on the AHP. The comparison of dispatching results between the cases where EVs participate and do not participate in scheduling is shown in Figs 11 and 12. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 11. Thermal power, load curve. https://doi.org/10.1371/journal.pone.0297855.g011 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 12. New energy generation power. https://doi.org/10.1371/journal.pone.0297855.g012 By comparing the power generation of thermal power units, it is found that the participation of electric vehicles in scheduling reduces the power generation during peak electricity consumption, while increasing the power during off peak electricity consumption, ensuring the long-term stable operation of thermal power plants. Comparing the load curve, it is found that when electric vehicles participate in scheduling, the fluctuation rate of power load decreases by 38.21% compared to when electric vehicles do not participate in scheduling. Comparing the output power of wind turbines and photovoltaic power generation units, it is found that when electric vehicles participate in scheduling, the abandonment rate of wind and photovoltaic power generation decreases by 31.24% compared to when electric vehicles do not participate in scheduling. The optimized scheduling results obtained from the participation of EVs in the scheduling can be compared with the simulation results without the participation of EVs. It is evident that EVs have significant advantages in balancing electricity loads and peak shaving. The fluctuation rate of electricity load is reduced by 38.21%. The participation of EVs in the scheduling reduces the output of thermal power generation during peak electricity consumption, while increasing its output during off-peak consumption, ensuring the long-term stable operation of the power plants. The abandonment rate of wind and PV power is reduced by 31.24%. Interval multi-objective optimization scheduling model solution Configure the NSGA-II algorithm with the population size of Npop = 100, the maximum number of iterations of gmax = 10000, the crossover probability of Pc = 0.9, and the mutation probability of Pm = 0.1. Based on historical data analysis, the fluctuation range of predicted values for wind power, photovoltaic power generation, and conventional power load obtained through MATLAB simulation are shown in Figs 13–15. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 13. Wind power prediction. https://doi.org/10.1371/journal.pone.0297855.g013 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 14. Photovoltaic power prediction. https://doi.org/10.1371/journal.pone.0297855.g014 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 15. Conventional electrical load power prediction. https://doi.org/10.1371/journal.pone.0297855.g015 Based on the above data and the improved NSGA-II algorithm, the interval multi-objective optimization model was solved to obtain the Pareto solution set in the objective functions, and the scatter plot of the median values in the objective function was drawn as shown in Fig 16. The distribution of each objective function value is shown in Figs 17–19. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 16. Scatter plot of median values. Incorporating uncertainty into interval multi-objective optimization models leads to a more reliable and robust set of solutions for the objective function values. https://doi.org/10.1371/journal.pone.0297855.g016 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 17. Economic objective. https://doi.org/10.1371/journal.pone.0297855.g017 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 18. Peak shaving and valley filling objective. https://doi.org/10.1371/journal.pone.0297855.g018 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 19. New energy consumption and integration objective. https://doi.org/10.1371/journal.pone.0297855.g019 Based on the AHP method for multi-objective decision-making, the optimal scheduling scheme was determined. The results are shown in Table 3. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 3. The objective function values of scheduling schemes considering extreme scenarios. https://doi.org/10.1371/journal.pone.0297855.t003 Compare the scheduling results of deterministic optimization model and interval optimization model for the most economical scheduling scheme. Assuming the most economical scheme is applied to the actual power grid, considering the impact of uncertainty factors, mainly due to the smoothing of load power fluctuations by thermal power units. When the output of thermal power units is subject to ramping constraints and power upper and lower limits, the load power fluctuations will be smoothed by the transmission power between the external power grid. The objective function values under different fluctuation coefficients are shown in Tables 4 and 5. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 4. The objective function value of the most economical scheme of the deterministic optimization model. https://doi.org/10.1371/journal.pone.0297855.t004 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 5. The objective function value of the most economical solution for interval optimization model. https://doi.org/10.1371/journal.pone.0297855.t005 Compared with the most economical solution obtained by the deterministic optimization model, the economic objective function interval of the most economical solution obtained by the interval optimization model has a larger median and a smaller interval width. The interval median and width of the peak shaving and valley filling objective function are both smaller. The interval median of the new energy consumption objective function is smaller, but the interval width is larger. This indicates that the introduction of uncertain factors into the interval optimization model enables the most economical solution obtained to pursue better economic performance while also considering the peak shaving and valley filling objectives and the new energy consumption objectives. This further proves that the solution obtained by the interval optimization model can better balance multiple objectives when extreme scenarios are considered, and has stronger robustness. EVs participating in the scheduling In the Pareto frontier, the distribution of objective function values is shown in Fig 2, where each point represents a Pareto optimal solution. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 2. The pareto frontier with the participation of electric vehicles. https://doi.org/10.1371/journal.pone.0297855.g002 Three optimization and scheduling schemes were selected for analysis based on their economic viability, high renewable energy absorption capacity, and effective peak-shaving and valley-filling effects. The most user-demand-oriented optimization and scheduling scheme was determined using the AHP decision-making analysis method. The objective function values under different scheduling schemes are shown in Table 2. Clearly, due to the lower cost of renewable energy generation compared to thermal power generation, the economic index of the power grid and the peak-shaving and valley-filling indexes and renewable energy absorption capacity index exhibit a trend of trade-off. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 2. The objective function values of scheduling schemes considering extreme scenarios. https://doi.org/10.1371/journal.pone.0297855.t002 Plot the three optimization and scheduling schemes that are the most economically viable, have the best peak-shaving and valley-filling effects, and have the highest renewable energy absorption capacity. 1) The best economic dispatch results. The power curves of the most economical scheduling scheme are shown in Figs 3 and 4. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 3. New energy generation power and EV power. https://doi.org/10.1371/journal.pone.0297855.g003 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 4. Thermal power generation and load power curve. https://doi.org/10.1371/journal.pone.0297855.g004 In the most economically efficient scheduling plan, during periods of high electricity prices and peak usage, the regional power grid needs to sell electricity to the large grid to meet the demand, resulting in an increase in power generation from wind, solar, and thermal sources. However, due to the limitations of power transmission capacity, some of the power generated from wind and solar sources may be wasted, resulting in a waste of energy resources. Conversely, during periods of low electricity prices and low usage, the regional power grid needs to buy more electricity from the large grid, resulting in a decrease in power generation from wind, solar, and thermal sources. As energy storage devices have lower participation costs in scheduling, they are used more during this time period, while the participation power of EVs decreases. This may cause a decrease in the system’s ability to balance peak and off-peak demand, as the scheduling capacity of energy storage devices is limited and may not be able to meet the demands during peak periods. 2) The results of the optimal peak-shaving and valley-filling scheme. The power curves of the optimal scheduling scheme for peak shaving and valley filling are shown in Figs 5 and 6. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 5. New energy generation power and EV power. https://doi.org/10.1371/journal.pone.0297855.g005 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 6. Thermal power generation and load power curve. https://doi.org/10.1371/journal.pone.0297855.g006 For the optimal peak-shaving and valley-filling scheme, in order to make the load curve as smooth as possible and reduce the peak-to-valley difference, EVs and energy storage devices are charged at high power during periods of low electricity demand and discharged at high power during periods of high electricity demand. As both thermal power generation and energy storage systems can control their power output, they can adjust their generation capacity flexibly to provide additional electricity during peak load periods and reduce the load during off-peak periods. However, this will increase the involvement of thermal power generation in load regulation, which will lead to a decrease in economic goals and a lower capacity for integrating new energy sources. 3) The results of the scheme with the highest capacity for accommodating new energy consumption. The power curves of the scheduling scheme with the highest new energy consumption capacity are shown in Figs 7 and 8. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 7. New energy generation power and EV power. https://doi.org/10.1371/journal.pone.0297855.g007 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 8. Thermal power generation and load power curve. https://doi.org/10.1371/journal.pone.0297855.g008 For the optimal strategy aimed at the highest new energy consumption capacity, the best approach is to maintain relatively high power generation of wind and photovoltaic power units at each scheduling time. However, since the output of wind and photovoltaic power generation fluctuates due to environmental factors, the peak-shaving and valley-filling effect of this approach is not ideal. In contrast, thermal power generation with higher output during the low valley of wind and photovoltaic power generation can compensate for the insufficient output of wind and photovoltaic power generation, while reducing output during the peak of wind power generation to avoid wasting resources. In addition, the power output of EVs participating in scheduling should have higher discharge power and lower charging power during the low valley of wind and photovoltaic power generation, and lower discharge power and higher charging power during the peak of wind and photovoltaic power generation. Based on the importance of three decision factors, a pairwise comparison judgment matrix A is constructed. After consistency check and hierarchical ranking, the normalized weight values of each index are W = [0.1047 0.6370 0.2583]. The optimal scheduling scheme is selected. The variable values of the optimal scheduling scheme are shown in Figs 9 and 10: Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 9. Power output of thermal power units and load power. https://doi.org/10.1371/journal.pone.0297855.g009 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 10. Grid connected power of new energy generation. https://doi.org/10.1371/journal.pone.0297855.g010 EVs do not participate in scheduling When electric vehicles do not participate in scheduling, we believe that they start charging when they enter the station, until the charging is completed or the electric vehicle leaves the station. Establishing a regional power grid scheduling model based on electric vehicles as regular loads that do not participate in scheduling, with the goal of power grid economics and new energy consumption as the objective functions. The decision variables include thermal power generation, photovoltaic power generation, and wind power generation, with the optimal Pareto solution set and frontier obtained. The optimal dispatching scheme is obtained based on the AHP. The comparison of dispatching results between the cases where EVs participate and do not participate in scheduling is shown in Figs 11 and 12. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 11. Thermal power, load curve. https://doi.org/10.1371/journal.pone.0297855.g011 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 12. New energy generation power. https://doi.org/10.1371/journal.pone.0297855.g012 By comparing the power generation of thermal power units, it is found that the participation of electric vehicles in scheduling reduces the power generation during peak electricity consumption, while increasing the power during off peak electricity consumption, ensuring the long-term stable operation of thermal power plants. Comparing the load curve, it is found that when electric vehicles participate in scheduling, the fluctuation rate of power load decreases by 38.21% compared to when electric vehicles do not participate in scheduling. Comparing the output power of wind turbines and photovoltaic power generation units, it is found that when electric vehicles participate in scheduling, the abandonment rate of wind and photovoltaic power generation decreases by 31.24% compared to when electric vehicles do not participate in scheduling. The optimized scheduling results obtained from the participation of EVs in the scheduling can be compared with the simulation results without the participation of EVs. It is evident that EVs have significant advantages in balancing electricity loads and peak shaving. The fluctuation rate of electricity load is reduced by 38.21%. The participation of EVs in the scheduling reduces the output of thermal power generation during peak electricity consumption, while increasing its output during off-peak consumption, ensuring the long-term stable operation of the power plants. The abandonment rate of wind and PV power is reduced by 31.24%. Interval multi-objective optimization scheduling model solution Configure the NSGA-II algorithm with the population size of Npop = 100, the maximum number of iterations of gmax = 10000, the crossover probability of Pc = 0.9, and the mutation probability of Pm = 0.1. Based on historical data analysis, the fluctuation range of predicted values for wind power, photovoltaic power generation, and conventional power load obtained through MATLAB simulation are shown in Figs 13–15. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 13. Wind power prediction. https://doi.org/10.1371/journal.pone.0297855.g013 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 14. Photovoltaic power prediction. https://doi.org/10.1371/journal.pone.0297855.g014 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 15. Conventional electrical load power prediction. https://doi.org/10.1371/journal.pone.0297855.g015 Based on the above data and the improved NSGA-II algorithm, the interval multi-objective optimization model was solved to obtain the Pareto solution set in the objective functions, and the scatter plot of the median values in the objective function was drawn as shown in Fig 16. The distribution of each objective function value is shown in Figs 17–19. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 16. Scatter plot of median values. Incorporating uncertainty into interval multi-objective optimization models leads to a more reliable and robust set of solutions for the objective function values. https://doi.org/10.1371/journal.pone.0297855.g016 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 17. Economic objective. https://doi.org/10.1371/journal.pone.0297855.g017 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 18. Peak shaving and valley filling objective. https://doi.org/10.1371/journal.pone.0297855.g018 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 19. New energy consumption and integration objective. https://doi.org/10.1371/journal.pone.0297855.g019 Based on the AHP method for multi-objective decision-making, the optimal scheduling scheme was determined. The results are shown in Table 3. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 3. The objective function values of scheduling schemes considering extreme scenarios. https://doi.org/10.1371/journal.pone.0297855.t003 Compare the scheduling results of deterministic optimization model and interval optimization model for the most economical scheduling scheme. Assuming the most economical scheme is applied to the actual power grid, considering the impact of uncertainty factors, mainly due to the smoothing of load power fluctuations by thermal power units. When the output of thermal power units is subject to ramping constraints and power upper and lower limits, the load power fluctuations will be smoothed by the transmission power between the external power grid. The objective function values under different fluctuation coefficients are shown in Tables 4 and 5. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 4. The objective function value of the most economical scheme of the deterministic optimization model. https://doi.org/10.1371/journal.pone.0297855.t004 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 5. The objective function value of the most economical solution for interval optimization model. https://doi.org/10.1371/journal.pone.0297855.t005 Compared with the most economical solution obtained by the deterministic optimization model, the economic objective function interval of the most economical solution obtained by the interval optimization model has a larger median and a smaller interval width. The interval median and width of the peak shaving and valley filling objective function are both smaller. The interval median of the new energy consumption objective function is smaller, but the interval width is larger. This indicates that the introduction of uncertain factors into the interval optimization model enables the most economical solution obtained to pursue better economic performance while also considering the peak shaving and valley filling objectives and the new energy consumption objectives. This further proves that the solution obtained by the interval optimization model can better balance multiple objectives when extreme scenarios are considered, and has stronger robustness. Conclusion The charging behavior of EVs is characterized by flexibility and uncertainty, and a large number of EVs with uncontrolled charging connected to the grid greatly increases the pressure of the grid. Electric vehicles have dual functions of energy storage and energy supply, and can be used as power grid energy storage devices when the vehicle is idle, which provides convenience for power grid peak regulation. In addition, electric vehicles can also be charged and discharged through the grid to achieve two-way flow of energy, thereby further improving the utilization efficiency of electric energy. As a dispatching resource, EVs can be scheduled for charging and discharging based on time-of-use pricing incentives, and the following conclusions can be drawn: Firstly, compared to simply exchanging electricity with the external power grid as a smoothing method, using electric vehicles as scheduling resources guides electric vehicle users to charge during low electricity consumption periods and discharge during peak electricity consumption periods, helping the regional power grid suppress load fluctuations and further reducing volatility. Secondly, EV aggregators can group multiple EVs into an aggregate according to the needs of the grid, so as to achieve the effect of load regulation and peak cutting and valley filling on the power grid. This approach can not only improve the dispatching capacity of the power grid, but also provide economic incentives for EV owners, promoting the popularization of EVs. Finally, the paper considered the impact of uncertainty factors on the dispatching process, and introduced the concept of interval number to make the model more realistic. The scheduling scheme obtained by interval multi-objective scheduling model is more robust and practical than that obtained by multi-objective scheduling scheme without uncertainty factors. Supporting information S1 Dataset. https://doi.org/10.1371/journal.pone.0297855.s001 (XLSX) S1 File. https://doi.org/10.1371/journal.pone.0297855.s002 (7Z) Acknowledgments The author would like to thank State Grid Henan Electric Power Company for providing the data, as well as the laboratory teachers for their assistance in algorithm and writing. TI - Optimal dispatching of regional power grid considering vehicle network interaction JF - PLoS ONE DO - 10.1371/journal.pone.0297855 DA - 2024-07-16 UR - https://www.deepdyve.com/lp/public-library-of-science-plos-journal/optimal-dispatching-of-regional-power-grid-considering-vehicle-network-oo2zNoh8MB SP - e0297855 VL - 19 IS - 7 DP - DeepDyve ER -