Multiobjective Synergistic Scheduling Optimization Model for Wind Power and Plug-In Hybrid Electric Vehicles under Different Grid-Connected Modes

In order to promote grid’s wind power absorptive capacity and to overcome the adverse impacts of wind power on the stable operation of power system, this paper establishes benefit contrastive analysis models of wind power and plug-in hybrid electric vehicles (PHEVs) under the optimization goal of minimum coal consumption and pollutant emission considering multigrid connected modes. Then, a two-step adaptive solving algorithm is put forward to get the optimal system operation scheme with the highest membership degree based on the improved ε constraints method and fuzzy decision theory. Thirdly, the IEEE36 nodes 10-unit system is used as the simulation system. Finally, the sensitive analysis for PHEV’s grid connected number ismade.The result shows the proposed algorithm is feasible and effective to solve the model. PHEV’s grid connection could achieve load shifting effect and promote wind power grid connection. Especially, the optimization goals reach the optimum in fully optimal charging mode. As PHEV’s number increases, both abandoned wind and thermal power generation cost would decrease and the peak and valley difference of load curve would gradually be reduced.


Introduction
With the promotion of the national energy conservation policy, implementing energy-saving generation scheduling in power system has been already put in the process.This means renewable energies are facing excellent development opportunities.With wind power's scale increasing, in 2013, wind power newly increased installed capacity was 16093.3MW and the gross installed capacity reached 92038.49MW.Both of the two indicator values rank first in the world.Currently, China is planning to build 8 million level wind power bases.The gross installed capacity of wind power will reach one hundred million in 2020 MW [1].
However, because of wind power output volatility, power sources away from the load centers, uncoordinated grid construction, and other factors, wind power absorptive problem has become the bottleneck of wind power in China.According to the statistical data released by The National Energy Board, in 2012, the total abandoned wind power reaches 20TW⋅h in China.The average equivalent utilization hours of wind power units are 1890 h and less than that in 2011.In Mongolia East, Jilin, Qinghai, and other wind resource-rich regions, the average equivalent utilization hours of wind power units are 1400 h and the abandoned wind power problem is very serious.To achieve the target, that renewable energy proportion in primary energy reaches 15% [2], abandoned wind problem should be solved mainly by improving the wind power absorptive capacity and efficiency of the grid.
There are already many researches focused on improving the wind power absorptive problem.Optimization methods are mainly aimed at optimizing wind power's backup service.They could be divided into four categories, namely, wind power-AGC collaborative grid connection [3][4][5][6], wind power and adjustable units collaborative grid connection [4][5][6][7], wind power and energy storage devices collaborative grid connection [8,9], and wind power and demand side response collaborative grid connection [10][11][12][13].Although all 2 Mathematical Problems in Engineering the optimization methods could bring some benefits, they also have certain deficiencies.For example, the adjustment capability of backup service is weak, and backup service cost is too high [14,15], which promote power system to search for new backup service for wind power.
Plug-in hybrid electric vehicles (PHEVs) have chargeable and dischargeable capacity, which makes PHEV have the potential to be the backup service for wind power grid connection.PHEV's orderly grid connection can bring multiple benefits such as achieving energy conservation by using electricity rather than oil, smoothing the load curve by controlling charging and discharging behavior, and promoting wind power grid connection.However, if PHEV's grid connection behavior is not under control, its arbitrary charging and discharging behavior would exacerbate the volatility of system load, which would decrease wind power consumption capacity.To solve these problems, wind power-PHEV synergistic scheduling optimization method needs to be studied.
Researches on wind power-PHEV synergistic scheduling optimization are mainly focused on two aspects, namely, model building and solution algorithm.As for model building, literature [16][17][18] established PEHV-thermal power synergistic scheduling optimization model and taxonomic studies on 4 kinds of charging modes' influence on the optimal combination of thermal power.In those literatures, there are only thermal units in the assumed system.Literature [19,20], respectively, established PHEV-wind power synergistic scheduling stochastic optimization model.Literature [21] put forward a PHEV-wind power synergistic scheduling optimization model in multiple time scales.And using the model, the literature studied wind power-PHEV synergistic scheduling's influence on the grid under freecharging mode.However, literature [19][20][21] assumes that PHEVs are in the free-charging mode, ignoring the elasticity relationships between user's electricity consumption behavior and electricity price.
In terms of solution algorithm, minimizing pollution emission becomes another optimization objective of wind power-PHEV synergistic scheduling under energy conservation requirements.This makes the optimization problems become multiobjective problems.Multiobjective problems need suitable solution algorithms to get the optimal solution sets.Literature [22,23] and literature [24], respectively, use POS, the simulated annealing, differential evolution, and other intelligent algorithms to solve unit commitment problems in wind power-PHEV synergistic scheduling.Conversional solution algorithms generally have difficulties in determining parameters and transforming constraints, which make their optimization levels not high [25,26].Intelligent algorithms can avoid these problems.But when individual extremes do not meet the multiobjective planning principles, intelligent algorithms are easy to fall into the local extreme point, limiting the search abilities [27][28][29].
Based on the analysis above, to overcome the deficiencies of the study on wind power-PHEV synergistic scheduling, this paper mainly contains 7 sections.Section 2 analyzed PHEV's charging and discharging characteristics and put forward the PHEV charging and discharging power output model.Section 3 simulates the uncertainty of wind power and constructs wind power scenarios simulation and scenario reduction strategy.Section 4 considered system demand and supply balance and units' output and backup service constraints, established wind power-PHEV synergistic scheduling models, respectively, and set the minimum coal cost and pollution emission cost as the optimization objectives during power generation.Section 5 puts forward a two-step solution algorithm based on the improved  constraints and fuzzy decision theory.Section 6 introduced the solving progress of the solution algorithm.Section 7 used the improved IEEE36 points system as the simulation system to analyze the optimization results in different charging and discharging modes and then made a sensitive analysis for the grid connection PHEV's number.

Charging Modes.
In terms of charging model, there are some research results.PHEV's charging modes can be divided from two aspects, namely, ordered or disordered charging, distinguish time periods, or continue charging.In this way, its charging modes can be classified into 3 kinds as follows: the no-control charging mode, the continuous charging mode, and the delayed charging mode [29].We can compare these three charging modes with the fully optimized charging mode.Their charging load curves are as shown in Figures 1,  2, and 3. Their characteristics can be described as follows.
(1) The no-control charging mode: in this mode users' charging behaviors are not influenced by anything.They just follow their own needs and habits.So PHEV's charging load would have strong volatility.Therefore PHEV's grid connection would bring some bad effects to electric system's short term operation stability.However, in the long term operation aspect, electric system could grasp users' charging habit and formulate corresponding measures to overcome the bad effects.
(2) The delayed charging mode: being different from the no-control charging mode, this charging mode applies time of use (TOU) and other control means to lead users' charging behaviors.Users are trend to adjust their charging behavior according to the load curve.This charging mode is usually considered to be a practical scenario, because in this mode PHEV mainly charge in the valley period and discharge in the peak period.So PHEV could play an important role in load shifting, smoothing the load curve.
(3) The continuous charging mode: in the no-control charging mode or delayed charging mode, PHEV is supposed to charge and discharge only once a day.That can be true for individual PHEV, but it is not so practical for public PHEV.In the continuous charging mode, as long as PHEV is not in operation, it can charge or discharge.So in this mode user's charging and discharging behavior would not be limited.PHEV's grid connection would increase, which means gasoline consumption and CO 2 emission would decrease.However, users' casual charging and discharging behavior would bring some uncontrolled effects to electric system's operation.For example, the volatility of load curve may increase.This is not conducive to wind power grid connection.
(4) The fully optimal charging mode: based on the delayed charging mode, the fully optimal charging mode takes measure to control users' charging behavior directly.In this mode, PHEV's charging or discharging would absolutely be in accordance with power system's supply and demand situation.If TOU is rational in the continuous charging mode, its optimization result would be almost the same with the fully optimal mode.So the load curve of the fully optimal charging mode is similar to the continuous mode.

Discharging
Modes.PHEV's discharging mode is influenced by battery type, capacity, related parameters, discharging cycle, discharging loop power, and other factors.Because it is hardly to decide discharging parameters, this paper uses existing research results [11] as data basis of discharging load distribution.
PHEV is used as a kind of backup service for wind power grid connection, in a certain extent; it is determined by controlling and estimating PHEV's discharging capacity.So the uncertainty of PHEV's discharging time and number brings considerable difficulties to accurately estimate the discharge capacity.Literature [31] has already studied PHEV's battery type, capacity, discharge cycle, discharge power, and other parameters' influence on PHEV's discharging behavior.This paper uses its result to make further study.Setting the maximum discharge capacity to be 9.6 kW, PHEV's discharging time and capacity 11 curve are shown in   Mathematical Problems in Engineering 2.3.PHEV Charging and Discharging Power Model.PHEVs need storage battery as their energy storage component to achieve charging and discharging behavior.Load change of storage battery can be tracked by energy storage controller [2].When it comes to the maximum capacity of the storage battery, the energy storage controller would stop storage battery from charging to ensure the battery life and operation safety.And when it comes to the minimum capacity, the energy storage controller would stop storage battery from discharging.Consider where  max and  min are the upper and lower limit of the storage battery capacity;   is the storage electricity of PHEV at time .
(1) When PHEV is in the charging status, (2) When PHEV is in the discharging status, where    and    are the storage electricity at time  of PHEV  and PHEV ;   is the self-consumed electricity during its discharging period of PHEV ;   is the self-consumed electricity during its charging period of PHEV ;    and    are discharging and charging power at time  of PHEV  and PHEV .
Generally, PHEV's charging or discharging power should not exceed 20% of the maximum capacity of the storage battery [25].The constraints can be described as

Simulate the Uncertainty of Wind Power
Wind power output is limited by the income-wind velocity.However, if the income-wind velocity is lower than the cutin wind velocity or is higher than the cut-out wind velocity, the wind farm will not generate power.The relationship between output power and the income-wind velocity could be expressed as where  * , is the actual available output of wind power unit  at time ;   is the rated output power of wind power units; V , is the cut-in speed; V , is the cut-out speed; V , is the rated wind speed; V  is the actual speed at time .

Wind Power Scenarios Simulation.
In the actual scheduling progress, wind power output can be determined through two ways, namely, forecast and simulation.In terms of wind power forecast method, there are classical forecast method [32,33], modern forecast method [34][35][36], and intelligent forecast method [37][38][39].The development of intelligent forecast method, especially, largely promoted wind power forecast technology.However, the practical applications show that the forecasting accuracy of these methods is still not able to reach the requirements of system scheduling.Therefore, this paper does not forecast wind farm's output power but uses scenario analysis method to simulate each power generation scenario.Assume the random variable of wind power output is  , ; its composition is detailed as follows: where   , is the forecast power of wind power;  , is the forecast error of wind power forecast.The forecast error is assumed to follow the normal distribution  , ∼ [0,  2 , ]; then  , could be though to follows the normal distribution . In order to gain the wind power scenarios, interval method is used to simulate wind power output.Divide wind power output into several intervals.Set the value of a point in the interval as wind power output expectation.When the number of intervals is sufficient, the forecast value could be regarded as the real output.The details are shown in Figure 1.
According to Figure 5, in each interval, there are three statues, namely, high, normal, and low.The expectation of wind power output in each statue is   , , where  = 1, 2, 3 are, respectively, corresponding with three statues.The probabilities of each output statue are    ; then the wind power output set in each scenario is  = {  , |  = 1, 2, . . ., } and the scenario probability of wind power output is  = ∏ ∈    .Figure 6 shows the wind power output forecast simulation scene scenario of three statues.introduces the Kantorovich distance [41], and set minimize the Kantorovich distance between the initial scenario and the reduced scenario as the optimize objective.

Wind Power Scenario Reduction
Assume  = [  ] × is the initial wind power scenario set and  = [ ξ ] × Ñ is the reduced scenario set.Then the Kantorovich distance between the initial scenario set  and the reduced scenario set  can be defined as Then the Kantorovich distance would be Define   to be equal to the summation of the scenario's occurrence probability in the initial scenario and the closest deleted scenario.This is detailed as () ∈ arg min (  , ξ ) , ∀ ∈ .
Then build wind farm scenario reduction optimization method based on (7)- (12).The optimized scenario reduction method can be expressed as min {∑ ∈   min ∉   (  , ξ ) :  ⊂ {1, 2, . . ., Ñ} ,  =  − Ñ} .(13) According to (13), the scenario reduction mechanism would influence the reduction result directly.To set an appropriate scenario deleted number, this paper puts forward the maximum reduction strategy: where ( 14) is to ensure the similar degree of the reduced scenario set and the initial scenario set in the required range.
The scenario reduction model consists of ( 7)- (14).To solve this model, this paper uses the multistage heuristic algorithm, referring to literature [42].

Wind Power-PHEV Synergistic Scheduling Optimization Model
where  is the quantity of power generation units;  is the total optimization period and, in this paper,  = 24;   is a 0-1 variable and when   = 1 it means that unit  is in operation at time  and   = 0 means unit  is shutdown at time ;   is the startup-shutdown cost of unit  at time ;   is the output of unit  at time ;   (  ) is the generation coal consumption cost function of unit  at time .Consider where   ,   ,   are coal cost coefficients of unit  determined by the regression of historical data of power generation.Consider where    is the cold-start cost of unit ;  ℎ  is the hot-start cost of unit ;  min , is the allowed minimum downtime;  off  is the continuous downtime of unit  at time ;   , is the cold-start time of unit ;  off  is the summation of minimum downtime and cold-start time of unit .

Minimize Pollutant Emission Cost Objective Function.
Pollutant emission cost of power generation can be calculated by the least squares method using historical pollutant emission data.The objective function is where  is the number of pollutant types; in this paper  = 3 and  = 1, 2, 3, respectively, means that the pollutant is CO 2 , SO where  is the number of all PHEVs;  is the number of PHEVs in charging; () is the load demand before PHEVs' grid connection at time ;  is the number of wind power units;   is the self-consumption rate of unit ;  , is the output of wind power at time ;   is the power consumption rate of the wind farm.

Thermal Units Output Constraints.
Thermal power output constraints include power generation upper and lower constraints, power climbing constraints, and the minimum startup and downtime constraints, as described in formulas ( 14)- (17).Consider where  max  and  min  are the upper and lower power generation limitation of thermal unit .Consider where Δ +  and Δ −  are the upper and lower power climbing limitation of thermal unit .Consider where  on  is the minimum startup time of thermal unit ;  on ,−1 is the continuous operation time of thermal unit  at time  − 1.Consider where  off  is the minimum downtime of thermal unit ;  off ,−1 is the continuous downtime of thermal unit  at time  − 1.

Wind Power Output Constraints.
Consider wherein  * , is the upper output limitation of wind power units at time .

PHEV Charging and Discharging Constraints. Consider
where   and   are, respectively, the number of charging and discharging PHEVs and are, respectively, the maximum number of chargeable and dischargeable PHEVs.

System Spinning Reserve Capacity. Consider the following:
(1) During the charging period, system's spinning reserve capacity is where   is the spinning reserve capacity before wind power's grid connection;  (1)   ( , ) and  (2)   ( , ) are, respectively, increased upper and lower limitation of spinning reserve capacity after wind power's grid connection.

Contracts Linearization.
To facilitate solving, we need to do linearization of formula (10) and formula (11).
(1) Initial state constraints are as follows: where  is the number of operating thermal units at initial state and if there is no initial state, then  = 0;  is the number of thermal units in shutdown statue at initial state.Set  0 as the period of operating thermal units at initial state; then (2) Startup and shutdown constraints are as follows:

Two-Step Adaptive Solving Algorithm
5.1.The Improved  Constraints.When solving multiobjective models, if a solution  can make all objective functions   () ( = 1, 2, . . ., ) achieve optimum, we could define  as an absolutely optimal solution.However, general contradictions exist between objectives and there is no absolutely optimal solution, but a set of optimal solutions, called Pareto optimal solution set [21].
This paper selects the improved constraints method to calculate multiobjective optimization problems.And then it uses the determination method in literature [43] to select feasible solutions.That means if   ∈  meets   () >   (  ),   () >   (  ) and  ̸ = , then   is a feasible solution.The constraints method has many advantages.But there are the following two big problems: (1) this method can only be used in the range of feasible solutions, which makes optimization results easy to fall into partial optimum; (2) results may not satisfy Pareto optimal solution's feasibility and are nondominant.According to literature [44,45], lexicographic optimization [33] and enhanced  constraints method can overcome these problems [29].
However, the enhanced  constraints method does not take different important degrees of each objective into consideration during its solving progress.This paper combines the enhanced  constraints method, lexicographic optimization method, and weight coefficient to get an improved  constraints method.And by defining an iteration parameter,    , we can get the feasible solution for the multiobjective problem, as shown in formula (31).Consider where   is the direction of objective ;   equals −1 means that objective  needs to be minimized, while   equals +1 means that objective  needs to be maximized.To avoid objectives' scale expand problems we pull     /  into objective functions.  is the weight of objective ;   is the residual variables of constraints;   is the range of objective function, which is determined by the decision attribute table and its determination progress is detailed in literature [46].Different from conventional weighting method, the improved  constraints method regards weights as the optimization variable.And the iteration parameter would update the weights during the solving progress.

Fuzzy Decision Theory.
To meet decision-makers' demand, this paper chooses fuzzy decision method to calculate membership degrees of the Pareto optimal solution set.Define membership function    , which stands for the optimal degree of objective  in the Pareto optimization scheme .
Then the membership degree function of Pareto optimization scheme  can be described as where   is the weight of objective  in (35);  is the number of elements in Pareto solution set;    and    , respectively, stand for the value and membership degrees of objective  and in Pareto optimization scheme .The membership functions of the optimization objective are shown in Figure 7.

Solving Progress
Based on the two-step adaptive solving algorithm, we can analyze different wind power-PHEV effectiveness in multigrid connected mode.Its solving steps are as follows.
(1) Input the original data of the model.
(3) Based on step (1) and step (2), the value range of objective function  2 () can be determined.And the improved  constraints method would help to solve the model and get the optimal solution sets.(4) Based on step (3), formula (36) could get the Pareto solutions and membership degrees of the objective functions.Choosing the solution set with the highest membership degree can get the satisfactory solution for the optimization model.Figure 8 shows the details.Table 1: The input-output tables of objective function.

Simulation
Objective function coal consumption and pollutant emission parameters are listed in Table 3.The emission parameters are referenced from paper [28].Based on wind power scenario simulation method, MATLAB is used to simulate 100 wind power output scenarios.Then the scenarios are reduced according to ( 7)-( 14).Finally we got 20 basic scenarios and wind power output, use the average output value of wind power in 20 basic scenarios as its available output power.The details are listed in Table 2. Load demand of a typical load day is also listed in Table 2.
Currently, there are three main types of PHEV, namely, BEVs, V2G, Triple-VG2 (equivalent to three V2G), and PCEV.Compared with other types, V2G has the advantage of being rechargeable and dischargeable, which makes it have better prospects.Therefore, this paper chooses 50000 V2G cars and studies V2G cars grid connection influence on wind power consumption.Assume that their average charging power is 1.8 kW and the maximum charging power is 2.4 kW.Charging period lasts for 6 hours and total charged electricity is 10.8 kWh [11].The number of PHEVs in charging period is equal to that in the discharging period.PHEV's power discharging consumption occupies 5.6% of its total electricity capacity, which means it can discharge 10.2 kWh.The discharging behavior is limited in 6 hours.

Simulation Results.
The simulation has been implemented in GAMS optimization software using CPLEX 11.0 linear solver from ILOG solver.The CPU time required for solving the problem for different case studies with a VAIO E series laptop computer powered by core i3 processor and 2 GB of RAM was less than 10 s.Firstly, we verified the validity and applicability of the solving algorithm in Section 7.2.1 to ensure the reliability of the result.And then the optimization schemes for different charging modes are compared.

Verification of Algorithm's Validity and Applicability.
This paper sets the no-control charging mode as the basic simulation scenario and chooses improved multiobjective nondominated sorting genetic algorithm-II (NSGA-II) [47] to solve the wind power-PHEV synergistic scheduling optimization model, where the maximum iterations time of NSGA-II algorithm is 10000, population size is 100, genetic operations crossover probability is 0.95, and mutation probability is 0.05.Results comparing the two algorithms are shown in Table 4.
According to Table 4, compared with the NSGA-II algorithm, the algorithm put forward by this paper could get a better optimal solution with lower iteration number and solving time.The grid-connected electricity of wind power is also more than that in the NSGA-II algorithm.
Analyzing the solution of the algorithm put forward by this paper can verify its applicability of solving wind power-PHEV synergistic scheduling optimization model under multiobjective functions.Figure 9 shows charging and discharging load distribution in the no-control charging mode.PHEV's charging time is concentrated in low load periods (1:00-5:00, 16:00-17:00, and 23:00-24:00).Its discharging time is concentrated in peak load periods (9:00-12:00 and 19:00-21:00).So the optimization result shows that PHEVs tend to  32) is applied to get the optimal result with the highest membership degree.
(1) Thermal Units' Output Structure.After PHEV's grid connection, thermal units' total output decreased.The output structure is significantly optimized.For example, 7# and 8# units no longer generate power and 6# unit only generates power in the unrestricted charging mode.And 3# and 4# units generate less while units 1# and 2# generate more.In terms of the output structure, consider the following.
(1) Compared with the no-control charging mode, in the other 3 charging modes, 1# unit achieved full capacity operation and the large capacity generators' output increased, namely, 2# and 3# units.
(2) In the continuous charging mode, charge and discharge times are not limited.Therefore, PHEV's grid connection is the maximum, which makes thermal units' output decrease.But for units with small installed capacity, like 4# and 5# units, their output would be more than that in the delayed charging mode or the fully optimal charging mode.
(3) In the fully optimal charging mode, thermal units' output allocation is determined by their installed capacity.That means units with big installed capacity would output much more than the small ones.In this way, the output structure comes to the optimal.Table 5 shows output allocation of thermal power units in four modes.
According to Table 6, thermal units' output before and after PHEV grid connection can be compared.The output allocation in the fully optimal charging mode is described as Figure 10 and that before PHEV grid connection is described as Figure 11.From these two figures, we can see the characteristics of the output allocation in the fully optimal charging mode as follows.
(1) Thermal units' output decreased obviously and the output curve is relatively smooth.
(2) In section of the base load, unit 1# and unit 2# supply the basis load demand.
(3) In section of the waist load, units 3# and 6# supply the waist load demand before PHEV grid connection and that changed to be unit 3# and unit 4# in the fully optimal charging mode.
(4) In section of the peak load, units 7# and 8# supply the peak load demand before PHEV grid connection.And that changed to be unit 5# in the fully optimal charging mode.
(2) Wind Power Grid Connection Situation.Before PHEV's grid connection, abandoned wind was 1511.39MW.But after PHEV's grid connection, abandoned wind decreased obviously.It is because in delayed charging mode or fully optimal charging mode, users' charging behaviors are influenced by the load curve.That makes load curve gentler and provide a stronger backup for wind power.In these two charging modes, the abandoned winds are, respectively, 233 MW⋅h and 26 MW⋅h.Wind power output in 4 charging modes is shown in Figure 12.
(3) Economic and Environmental Benefits.PHEV's economic and environmental benefits in 4 modes are shown in Table 7.After its grid connection, both economic and environmental benefits increased.In the fully optimal charging mode, thermal units' startup and shutdown costs decreased by 50%, coal consumption cost decreased by 14%, emission of CO 2 , SO 2 , and NO  , respectively, decreased by 7%, 6%, and 6%, and abandoned wind decreased to 25.64 MWh.PHEV's grid connection would bring, obviously, economic and environmental benefits.To achieve the maximum benefit, electric enterprise should take relative incentive mechanisms to guide PHEV users' rational charging and discharging.

Sensitive Analysis. As government's support to develop
PHEV is increasing, PHEV's number grows.This means the charging and discharging behaviors of PHEV would influence power system's load curve directly.To verify the influence of PHEV's grid connection on power system's load curve economic benefit and environmental benefit, this paper sets 25, 50, 75, and 1000 thousand PHEVs in succession to do system simulation.The results are shown in Figures 13 and  14 and peak and valley differences of different PHEV number are listed in Table 8.
As the number of PHEVs increases, system's load curve becomes gentler, which means more obvious load shifting effect.When PHEV number is 25 thousand, the peak load is 2700 MW and the valley load is 1700 MW.And when PHEV  number is 100 thousand, the peak and valley loads changed to 2545 MW and 1830 MW.Their difference decreased by 29%.
According to Figure 14, as the number of PHEVs increases, abandoned wind decreases.When PHEV number is 100 thousand, abandoned wind is 72.5 MW, 77% less than that of 315.4 MW when PHEV number is 25 thousand.Furthermore, part of the PHEV would use wind power to charge at the valley time and discharge at the peak time.Thermal power's generation cost decreased when the number of PHEVs increased.

Conclusions
To analyze PHEV's grid connection benefits, this paper established a wind power-PHEV synergistic scheduling optimization model and put forward a two-step adaptive solving algorithm based on improved  constraints theory and fuzzy decision theory.According to the simulation results, there are some conclusions as follows.
(1) This paper uses the improved  constraints method to get Pareto optimal solution set and fuzzy decision theory to select a satisfactory solution with the highest membership degree.This two-step adaptive algorithm can be more suitable to solve the multiobjective wind power-PHEV synergistic scheduling optimization model.
(2) PHEV's grid connection can help in load shifting, decrease abandoned wind and power generation cost, and obviously bring economic and environmental benefits.And in delayed charging pattern and fully optimal charging pattern, system would get the optimal benefits.
(3) Compared with the no-control charging mode, the delayed charging mode has the maximum grid connection electricity amount, the least thermal power generation, and little abandoned wind.However, since users' casually charging behavior, PHEV grid connection can bring big impact on system's safe and stable operation.So the abandoned wind in this mode is higher than that in the delayed charging mode or fully optimal charging mode.
(4) In the fully optimal charging mode, abandoned wind comes to be the least.That is because in this mode users would adjust their charging behavior according to the load curve.Then the load shifting effect would be increased and PHEV grid connection can help wind power grid connection, increasing economic and environment benefits.
(5) Based on the conclusions above, to make full use of PHEV grid connection, we need to commence from two aspects.On one hand, use incentive policies to promote the development of electric vehicles.And on the other hand, use reasonable demand side electricity to guide users' charging behavior, which would be an important research direction on wind power-PHEV synergistic scheduling.

Figure 4 .
In the figure, almost all PHEVs can finish their discharging work in 6 hours.
Strategy.The basic concept of scenario reduction is comparing a scenario with other scenarios and removing the closest one.And the bigger the scenario number, the bigger the workload of scenario reduction.To overcome this problem, this paper

Figure 7 :
Figure 7: The corresponding membership function of optimization objective.

2 𝐹 min 2
Remark: * indicates that the objective function is the target to solve.

Figure 9 :
Figure 9: Electric vehicle charging and discharging time distributions in the no-control charging mode.

Figure 10 :Figure 11 :Figure 12 :
Figure 10: Thermal units' output allocation in the fully optimal charging mode.

Figure 14 :
Figure 14: Economic and environmental benefits comparison.
2 , and NO  ;    ,    , and    are pollutant emission parameters of unit .
1 (  ) can be expressed as piecewise functions.When  , ∈ [ min ] into  segments, so that quadratic functions

Table 1 .
(2)Based on step (1), objective functions  1 () and  2 () The objectives  1 () and  2 () are, respectively, set as a singleobjective optimization model and it could get objective values under different optimization objectives.The results form a decision attribute table of the objectives, as shown in

Table 2 :
Wind power output and power load demand.

Table 3 :
Coefficients of coal-fired power units.

Table 4 :
Results comparing the two algorithms.

Table 5 :
Best Pareto-optimal solution in different modes.

Table 6 :
Output allocation of thermal power unit in four modes (MW).

Table 7 :
System optimization results contrast in different modes.

Table 8 :
Sensitive analysis result on PHEV number.