Optimal Day-Ahead Scheduling of a Smart Distribution Grid Considering Reactive Power Capability of Distributed Generation

Abstract: In the traditional paradigm, large power plants provide active and reactive power required for the transmission system and the distribution network purchases grid power from it. However, with more and more distributed energy resources (DERs) connected at distribution levels, it is necessary to schedule DERs to meet their demand and participate in the electricity markets at the distribution level in the near future. This paper proposes a comprehensive operational scheduling model to be used in the distribution management system (DMS). The model aims to determine optimal decisions on active elements of the network, distributed generations (DGs), and responsive loads (RLs), seeking to minimize the day-ahead composite economic cost of the distribution network. For more detailed simulation, the composite cost includes the aspects of the operation cost, emission cost, and transmission loss cost of the network. Additionally, the DMS effectively utilizes the reactive power support capabilities of wind and solar power integrated in the distribution, which is usually neglected in previous works. The optimization procedure is formulated as a nonlinear combinatorial problem and solved with a modified differential evolution algorithm. A modified 33-bus distribution network is employed to validate the satisfactory performance of the proposed methodology.


Introduction
Global warming, energy crises, and ongoing serious environmental pollution have promoted the development of renewable energy greatly.With the benefit of low pollution and economic cost, they are encouraged to be used instead of fossil fuels.However, the high penetration of distributed energy resources (DERs) would pose a more complex operational situation for the distribution network [1].Under the circumstances, it is significant to assure the reliable and secure operation of the distribution; additionally, the reasonable utilization of DERs is crucial to the sustainable development of renewable energies.
Over these years, the smart grid has been presented as a good notion to deal with the new problems accompanied by the integration of DERs.As the rapid development of intelligent electronic devices and information communication technology, the concept of the smart grid has been extended to the distribution level.In 2008, at the Conference International des Grands Reseaux Electriques (CIGRE), the active distribution network (ADN) was first put forward [2].Different from the traditional network characterized with passivity and single-phase, the DGs and demand-side management programs, such as RLs connected in the ADN are online-controllable.The distribution management system (DMS), Energies 2016, 9, 311 2 of 17 as the core of secure and economical operation of the distribution network, is responsible for the optimal operation scheme in the ADN, considering the participation of DERs.
Relevant research has been reported in recent years.A day-ahead operational scheduling paradigm for a distribution network has been presented in [3] where DGs are judiciously taken into consideration for their goodness factors.Results showed significant benefits in terms of reduced losses and cost savings.The authors in [4] have proposed a DMS algorithm to manage an ADN.The algorithm targets the minimization of the cost of system operation through the coordination of flexible network topology with the management of energy resources.In [5], a short-term energy resource scheduling methodology for a smart grid has been developed considering the intensive penetration of DGs and demand-response programs.The methodology aims at minimizing the operation costs and involves different anticipation times: day-ahead, hour-ahead, and 5-min-ahead.The authors in [6] have presented a comprehensive day-ahead operation model that includes network technical constraints and can participate in both energy and reserve markets.Considering the minimization of cost of all integrated DERs and the minimization of the voltage magnitude difference in all buses, a multi-objective paradigm has been proposed to determine the optimal solution of resource scheduling for a smart grid in [7].The results have shown the importance of taking into account the reactive scheduling in the distribution network.
Although remarkable achievements have been obtained, the previous works generally ignore the pollutant emission and reactive power support capabilities of wind and solar power in the optimal resource scheduling of the smart grid.As atmospheric pollution and global warming have become one of the main critical environmental issues, the electrical industry, most of which use fossil energy as fuel, should make a contribution to environmental protection [8].Thus, the emission ought to be considered in the distribution scheduling.Meanwhile, the authors in [9] have demonstrated the steady-state reactive power loading capability of doubly-fed induction generators (DFIG)-based wind turbines.The authors in [10] have achieved a significant reduction in losses by reasonable utilization of the wind farm reactive power capability.In [11], the reactive capability limits of different renewable DGs have been discussed and incorporated in the distribution system expansion planning strategy.A scheme for a reactive power market for DGs in the distribution network is presented in [12], where DGs are not only used to deliver active power, but also encouraged to take part in the reactive power support.Therefore, in the optimal scheduling of smart distribution grid, it is sensible to take advantage of the DGs' reactive power capability, especially wind turbines (WT) and photovoltaics (PV), which are usually neglected.
Based on the preceding discussions, an intelligent scheduling program for DMS is presented in this paper.The proposed DMS targets to seek the optimal control of available active elements through the minimization of composite economic cost of distribution network.The composite cost consists of three components.The first is the operation cost, including the cost of purchasing power from the whole market, the cost of DG dispatches, and the cost of RL participations; the second is the emission cost, including pollutant cost of purchasing power and DG power; the third is the cost of network losses.Furthermore, the reactive power capabilities of WT and PV are analyzed; then they are wisely scheduled to settle the optimal reactive power dispatch of the ADN, as well as the conventional DG.The technical constraints of the whole system and the DERs integrated to the network are accommodated in the AC power flow fashion.A modified differential evolution algorithm (MDE) is adopted to resolve the proposed program.The obtained results are investigated to assess the effect of proposed DMS on the technical and economical improvements.
The rest of this paper is organized as follows.After the introductory section, Section 2 addresses the reactive power capabilities of different DG systems integrated in the distribution system.Section 3 presents the proposed scheduling framework.Section 4 focuses on the optimization technique.Section 5 carries out several test cases and discusses the numerical results.Finally, Section 6 concludes the paper.

Reactive Power Capabilities of Different DG Systems
In the paper, the wind systems based on doubly-fed induction generators (DFIGs), PV systems are based on voltage source inverters (VSI), and diesel generators are considered as the DERs that are integrated to the distribution.The static reactive power capabilities of the different DERs are discussed as follows.

Wind System
DFIG can adjust its reactive power injection or absorption and maintain the maximum active power output at the same time.The maximum reactive power injection of the DFIG-based wind system is mainly limited by the maximum rotor current I r n , active power output P WT n,st and slip s n [13], which is expressed as: where Q WTin n,st is the maximum reactive power injection limit of the WT at bus n over a system state st; V s n is the stator voltage of the WT at bus n; χ s n is the stator reactance; χ m n is the magnetizing reactance.The maximum reactive power absorption Q WTab n,st is determined by the steady-state stability limits [13], which is expressed as:

PV System
The active power output of the VSI-based PV system depends on the light intensity, and the reactive power output can be independently controlled by the power angle and modulation index [14].The maximum reactive injection Q PVin n,st and absorption Q PVab n,st is expressed as: where S ipv n,st is the rated apparent power of the inverter; P PV n,st is the active power injection of the PV system.

Diesel Generator
A diesel generator (DE) is based on the traditional synchronous machine.The reactive power capability is controlled by the field excitation while maintaining stable active power output.The maximum reactive power injection Q DEin1 n,st and absorption Q DEab n,st are limited by the maximum stator current I s n [15], which can be calculated by: where V n,st is the system voltage at bus n; P DE n,st is the active power injection of DE.The maximum reactive injection Q DEin2 n,st is also limited by the maximum rotor voltage E r n [15].
where X d n is the direct axis component of synchronous reactance.

The Proposed Scheduling Framework
In this section, a mathematical model for DMS is presented to seek the optimal scheduling of the purchasing power from the whole market, and make decisions on DERs options including DEs, WTs, PVs, and RLs.The proposed scheme schedules the available control variables of ADN in a 24 h day-ahead time horizon.

Objection Function
In addition to the economic benefit, the impacts of environment pollution and technology revenue are also considered in the scheduling framework.Thus, the composite economic cost of distribution network Fpx, uq including operation cost F1px, uq, emission cost F2px, uq, and network loss cost F3px, uq are proposed to be the objective function.That is: where x and u respectively denote the vectors of control and dependent variables.In Equation ( 7), V is the vector of bus voltages; Ñ S f is the vector of apparent power flow of feeder f.
x " r 3.1.1.Operation Cost F1(x, u) The distribution operation cost is composed of several terms.The first term is the cost of purchasing power from the day-ahead wholesale market; the second term is DG start-up and shut-down costs; the third term is DG fuel consumption cost and maintenance cost; the fourth term is the maintenance cost for WTs and PVs; and the fifth term is cost of RLs participations.WTs and PVs are powered by renewable energy, so its fuel costs can be ignored.Therefore, the operation cost is expressed as: The wholesale power is considered to be thermal power.Carbon dioxide (CO2), sulfur dioxide (SO2), and nitrogen oxide (NOx) are the main pollutants accompanied with the purchasing power and DGs power.WTs and PVs are deemed to be non-polluting.Thus, the emission cost is: where M is set of pollutant species; α j is emission coefficient of pollutant j, $/kg; E j is emission rate of pollutant j for thermal power, kg/MW; F i,j is emission rate of pollutant j for ith DG, kg/MW.In Equation (10), the first term is emission cost for the purchasing power from whole market; the second term is emission cost for DGs dispatches.
3.1.3.Network Loss Cost F3(x, u) The active power loss of distribution network can be calculated by power flow.The cost is: where ρ loss,t is unit cost of active power loss at time interval t, $/MW; P loss,t is active power loss of distribution at time interval t.

Steady-State Security Constraints
In these constraints, P D , Q D are, respectively, active and reactive power loads in each bus; V i , V j are voltage amplitudes of buses i and j; N node is set of system buses; N O is the set of system buses, excluding bus i; Y ij is the admittance between bus i and j; δ i , δ j are voltage phases of buses i and j; θ ij is phase of Y ij .S Grid is the substation apparent power capacity; the superscripts "max" and "min" respectively denote the corresponding upper and lower limits.RR is ramping capability of DG.PF RL is constant power factor for RLs.N f is set of distribution feeders.In Equations ( 16), (18), and (20), the reactive power limits of WTs, PVs, and DGs should be associated with their steady-state characteristics as discussed in Section 2.

Optimization Technique
The proposed day-ahead optimal scheduling framework is a nonlinear combinatorial optimization problem.In this paper, a modified differential evolution (MDE) algorithm is presented to solve it.The differential evolution algorithm (DE) was first put forward by Storn and Price [16].Its principle is simple and it is easy to be understood and implemented with fewer controlled parameters.DE has been used to solve different engineering optimization problems and achieved good effects [17][18][19][20].However, the control parameters in DE have a great effect on the algorithm performance and their appropriate values must be obtained through numerous tests.

Original Differential Evolution Algorithm
DE is a kind of population-based algorithm, which can remember the best solutions and share the information within the population.Supposing the population size is N, the ith individual at Gth iteration can be expressed as X i,G pi " 1, 2, ¨¨¨, Nq.Then a new population is generated through mutation, crossover, and selection operation [16].
For each individual X i,G , three different individual vectors are randomly selected from the current iteration G; then the mutation vector V i,G`1 can be produced by: where F is a scaling factor, which scales the difference vector pX r2,G ´Xr3,G q.The trial vector U i,G`1 is generated by crossover operator between the target vector X i,G and the mutation vector V i,G`1 .
where j P r1, 2, ¨¨¨, Ds, D is the dimension of the individual vector; rand(j) is a uniform random number within [0,1]; i rand is a random integer within r1, 2, ¨¨¨, Ds.CR is the crossover rate.After mutation, the selection operator is performed to select the next generation, which is expressed as: In original DE, all individuals in the population use the same F and CR, and their values remain unchanged in the evolution.Additionally, if the optimization problem is different, the values of F and CR should be reset.In order to make the parameters selection independent of optimization problems and improve the robustness of DE, a parameter self-learning strategy is proposed.In the strategy, the scaling factor F and crossover rate CR participates in individual coding and each individual X i,G corresponds to one F i,G and CR i,G .When initialized, F i,0 and CR i,0 are randomly generated within a given range.
In the evolution process, if an individual has not yielded one better child in successive m iterations, such as m = 5, it can be considered that the current control parameters of the individual is not appropriate, and should be reset; conversely, if an individual has yielded one or more better children in each successive m iterations, it can be considered that the current control parameters of the individual is appropriate, and should be retained.Thus, it is defined that in each successive n iterations, such as n = 10, the individual that has yielded the most, better, children is the "local optimal individual" and its corresponding control parameters are the "local optimal parameters".Furthermore, in the current Gth iteration, the individual that has yielded the most, better, child individuals is the "global optimal individual" and its corresponding control parameters are the "global optimal parameters" in the current population.The inappropriate individual control parameters should learn from the "local optimal parameters" and "global optimal parameters" in the current population when reset.Thus, a self-learning parameters modification strategy is proposed.
where k is the number of the individual whose control parameters need to be reset; C k,l and C l are, respectively, the cumulative improvement times of the individual X k and the "local optimal individual" in successive n iterations; C k,G and C g are, respectively, the cumulative improvement times of the individual X k and the "global optimal individual" in the current Gth iteration; F l and CR l are the "local optimal parameters"; F g and CR g are the "global optimal parameters"; w is the weight factor which is equal to 0.5 here; p1 ´Ck,l C l `1 q and p1 ´Ck,G C g `1 q are named as the learning step.It can be seen that, for an individual, the smaller the cumulative improvement times is, the larger the learning step is, namely learning more from the "local optimal parameters" and "global optimal parameters".Therefore, memory and information exchange functions are added to the proposed parameters modification strategy and the strategy has the ability of self-learning.

Boundary Handling
In the evolution process, if one dimension of an individual X i " rx 1 , x 2 ¨¨¨, x D s exceeds its corresponding feasible regionrx min , x max s, it is generally reset to the boundary, namely x i " x max or x i " x min .However, many individuals that gather around the boundary are harmful for the algorithm convergence; thus, a boundary handing method is presented.
If x i ą x max : where h is a user-specified constant.Through the method, the individuals which beyond the limits are distributed in a feasible region instead of on the boundary, so it increases the population diversity to a certain extent and accelerates the algorithm convergence.

Case Studies
This section is devoted to explore the efficiency of the proposed day-ahead optimal scheduling framework through the modified IEEE 33-bus distribution system [21,22].Figure 1 shows the single line diagram of the 33-bus radial distribution system.In the network, three diesel-based DGs are supplemented at buses 7, 13, and 25 whose characteristics are provided in Table 1.Additionally, the system includes two 2 MW DFIG-based WTs and two 2 MW VSI-based PVs at buses 14, 17, 2, 4 and 31, respectively.The loads at buses 24 and 30 are considered as RLs with constant PF which can be curtailed up to 10 percent of their actual demand.The contract price for 1 MW curtailed by RLs is determined as $115/MW.Bus 1 is the only substation connected to the wholesale market.
where h is a user-specified constant.Through the method, the individuals which beyond the limits are distributed in a feasible region instead of on the boundary, so it increases the population diversity to a certain extent and accelerates the algorithm convergence.

Case Studies
This section is devoted to explore the efficiency of the proposed day-ahead optimal scheduling framework through the modified IEEE 33-bus distribution system [21,22].Figure 1 shows the single line diagram of the 33-bus radial distribution system.In the network, three diesel-based DGs are supplemented at buses 7, 13, and 25 whose characteristics are provided in Table 1.Additionally, the system includes two 2 MW DFIG-based WTs and two 2 MW VSI-based PVs at buses 14, 17, 2, 4 and 31, respectively.The loads at buses 24 and 30 are considered as RLs with constant PF which can be curtailed up to 10 percent of their actual demand.The contract price for 1 MW curtailed by RLs is determined as $115/MW.Bus 1 is the only substation connected to the wholesale market.In the test cases, the main assumption is that the forecast of WT power, PV power, and load profiles are fairly accurate; thus, the prediction errors are considered negligible.Figure 2 shows the network day-ahead predicted total load and the scaled wholesale market prices are demonstrated in Figure 3. Supposing that two WTs and two PVs are, respectively, faced with similar geographical conditions, the power generations by WTs and PVs are shown in Figure 4. Table 2 lists the maintenance cost of DERs.Emission quantity and emission cost for thermal power and DGs power are provided in Table 3 [23,24].The cost coefficient of network loss is considered to be equal to the electricity price.In the test cases, the main assumption is that the forecast of WT power, PV power, and load profiles are fairly accurate; thus, the prediction errors are considered negligible.Figure 2 shows the network day-ahead predicted total load and the scaled wholesale market prices are demonstrated in Figure 3. Supposing that two WTs and two PVs are, respectively, faced with similar geographical conditions, the power generations by WTs and PVs are shown in Figure 4. Table 2 lists the maintenance cost of DERs.Emission quantity and emission cost for thermal power and DGs power are provided in Table 3 [23,24].The cost coefficient of network loss is considered to be equal to the electricity price.
Energies 2016, 9, 311 9 of 17 here.First, DMS determines the optimal solution of the available active elements including the reactive power support capabilities of WTs and PVs.Second, WTs and PVs are evoked only to output active power without the reactive power support capabilities.In the simulations, the proposed algorithm (MDE) is adopted to solve the problem.In the algorithm, the population size N and maximum number of iterations max G are, respectively, set to 100 and 200 for all cases; when initialized, the scaling factor ( ,0 i F ) and crossover rate ( ,0 CR i ) are randomly generated within the given range (0,1].here.First, DMS determines the optimal solution of the available active elements including the reactive power support capabilities of WTs and PVs.Second, WTs and PVs are evoked only to output active power without the reactive power support capabilities.In the simulations, the proposed algorithm (MDE) is adopted to solve the problem.In the algorithm, the population size N and maximum number of iterations max G are, respectively, set to 100 and 200 for all cases; when initialized, the scaling factor ( ,0 i F ) and crossover rate ( ,0 CR i ) are randomly generated within the given range (0,1].In order to examine the efficiency of the proposed scheduling framework, two scenarios are devised here.First, DMS determines the optimal solution of the available active elements including the reactive power support capabilities of WTs and PVs.Second, WTs and PVs are evoked only to output active power without the reactive power support capabilities.In the simulations, the proposed algorithm (MDE) is adopted to solve the problem.In the algorithm, the population size N and maximum number of iterations max G are, respectively, set to 100 and 200 for all cases; when initialized, the scaling factor ( ,0 i F ) and crossover rate ( ,0 CR i ) are randomly generated within the given range (0,1].In order to examine the efficiency of the proposed scheduling framework, two scenarios are devised here.First, DMS determines the optimal solution of the available active elements including the reactive power support capabilities of WTs and PVs.Second, WTs and PVs are evoked only to output active power without the reactive power support capabilities.In the simulations, the proposed algorithm (MDE) is adopted to solve the problem.In the algorithm, the population size N and maximum number of iterations Gmax are, respectively, set to 100 and 200 for all cases; when initialized, the scaling factor (F i,0 ) and crossover rate (CR i,0 ) are randomly generated within the given range (0,1].

Considering the Reactive Power Support Capabilities of WTs and PVs
The proposed scheduling model seek the optimal scheduling of active elements by minimizing the economic costs of ADN.In order to assess the effects of different cost functions on the scheduling results, four cases are tested.
Case 1: the operation cost is considered as the objective function; Case 2: the emission cost is considered as the objective function; Case 3: the network loss cost is considered as the objective function; and Case 4: the composite cost including operation cost, emission cost, and network loss cost is considered as the objective function.
The energy resource scheduling for Cases 1-4 are demonstrated in Figures 5-8.It can be seen from Figure 5 that, in Case 1, DGs and RLs are made on during peak hours 12-20 to meet the network load when their costs are lower than the market electricity price; during hours 1-11 and 21-24, the network load demand is less and the market electricity price is low, DGs and RLs are switched off.The distribution purchases all of its requirements from the wholesale market.As the emission cost of DGs' power are less than the purchased power (thermal power) and RLs' participation can decrease the load demand to reduce the emissions; meanwhile, the dispatches of DGs and RLs can also reduce the network loss with the alteration of the network power flow distribution; hence, Figures 6 and 7 shows that DGs and RLs are switched on during all the hours in Cases 2-3.Case 4 takes the composite cost as objective function and notable changes can be observed from Figure 8.Compared with Case 1, RLs and DGs are, respectively, scheduled from hours 10 and 11; their dispatch time continues to hour 21.In these hours, due to the participation of DERs, the operation cost may increase, but the emission and active power loss will drop largely.Thus, the composite cost would be a minimum.Furthermore, in comparison to Case 1, the RLs' participation in Case 4 are significantly more.
Table 4 depicts the optimal hourly operation costs, emission costs, and active power loss costs of the distribution network for Cases 1-4.As it is seen, the total operation cost in Case 1 ($73,229.59),emission cost in Case 2 ($12,879.09),and network loss cost in Case 3 ($983.94)are, respectively, minimum for the individual cost in Cases 1-4.The corresponding individual costs in Case 4 are all more than them; but the composite cost in Case 4 is the minimum in Cases 1-4, which is $89,415.89.It can be observed that the total operation cost in Case 4 is $214.64 higher than Case 1; but with respect to the emission cost and network loss cost, Case 4 is, respectively, $654.26 and $297.03lower than Case 1.Thus, the composite cost for Case 4 is $736.65 lower than Case 1.Compared with Cases 2 and 3, due to the dispatches of DERs in off-peak hours for Cases 2-3, namely hours 1-11 and 21-24, the operation costs are too high; hence, the composite cost for Case 4 is, respectively, $7,542.85 and $7,569.25 lower than Cases 2 and 3.
The distribution purchases all of its requirements from the wholesale market.As the emission cost of DGs' power are less than the purchased power (thermal power) and RLs' participation can decrease the load demand to reduce the emissions; meanwhile, the dispatches of DGs and RLs can also reduce the network loss with the alteration of the network power flow distribution; hence, Figures 6 and 7 shows that DGs and RLs are switched on during all the hours in Cases 2-3.Case 4 takes the composite cost as objective function and notable changes can be observed from Figure 8.Compared with Case 1, RLs and DGs are, respectively, scheduled from hours 10 and 11; their dispatch time continues to hour 21.In these hours, due to the participation of DERs, the operation cost may increase, but the emission and active power loss will drop largely.Thus, the composite cost would be a minimum.Furthermore, in comparison to Case 1, the RLs' participation in Case 4 are significantly more.Cases 1-4 demonstrate the superiority and effectiveness of the proposed scheduling framework.The presented algorithm MDE have solved the nonlinear optimization problem effectively.In order to further illustrate the advantage of MDE, Particle Swarm Optimization (PSO) [25,26], DE [16], and MDE are used, respectively, to solve Case 4 ten times.In PSO and DE, the population size and maximum number of iterations are all set to 100 and 200.The acceleration constants c1, c2 are set to 2; the maximum and minimum inertia weight factor ω are set, respectively, to 0.9 and 0.4 in PSO.In DE, the scaling factor (F) and crossover rate (CR) are set, respectively, to 0.2 and 0.6.The statistical result is shown in Table 5.It can be seen that MDE obtains the values of composite cost and its standard deviation.Thus, MDE has the better performance on optimizing capacity and algorithm robustness.

Non-Considering The Reactive Power Support Capabilities of WTs and PVs
To assess the effect of reactive power support capabilities of WTs and PVs on the scheduling result, Cases 5-8 are added.They are similar to Cases 1-4, but in the Cases 5-8, the WTs and PVs are evoked only to output active power without the reactive power support capabilities.6.As it is seen, the total composite cost in Case 8 ($90,400.01)is also the minimum compared with Cases 5-7.Thus, in Cases 5-8, which do not consider the reactive power support capabilities of WTs and PVs, the proposed composite cost objective function is also effective and achieves the minimum composite economic cost of the distribution network.Comparing Tables 4  and 6 four comparisons between Cases 1 and 5, Cases 2 and 6, Cases 3 and 7, and Cases 4 and 8 are developed, it can be found that considering the reactive power support capabilities of WTs and PVs, the total daily composite cost in Cases 1-4 are $90,152.54,$96,958.74,$96,985.14 and $89,415.89which are, respectively, $934.77,$833.78,$866.66 and $984.12lower than Cases 5-8.Therefore, the conclusion can be drawn that implementation of the reactive power support capabilities of WTs and PVs can effectively decrease the economic cost of the network.To further illustrate the effect of reactive power support capabilities of WTs and PVs, Case 4 and Case 8 are used for analysis.In Cases 4 and 8, the proposed composite cost is both considered as the objective function.Figure 13 explores the imported total reactive power from the main grid in the optimal day-ahead scheduling framework.It can be observed that because of the reactive power support capabilities of WTs and PVs, the total reactive power injected from the main grid in Case 4 is much less than Case 8.Meanwhile, as Figure 14 shows, the total apparent power injected from the main grid in Case 4 is also less than Case 8. Additionally, the network loss is usually concerned when the reactive power of distribution is optimally dispatched.Figure 15 depicts the hourly network losses for Cases 4 and 8.In Case 8, the main reactive power source is the wholesale market, but in Case 4, WTs and PVs integrated into ADN are also used as reactive power sources to accommodate the network demand.Thus, the reactive power regulation means Case 4 is more flexible.Accordingly, as it is seen, the total network losses for Case 4 are evidently less than Case 8.  To further illustrate the effect of reactive power support capabilities of WTs and PVs, Case 4 and Case 8 are used for analysis.In Cases 4 and 8, the proposed composite cost is both considered as the objective function.Figure 13 explores the imported total reactive power from the main grid in the optimal day-ahead scheduling framework.It can be observed that because of the reactive power support capabilities of WTs and PVs, the total reactive power injected from the main grid in Case 4 is much less than Case 8.Meanwhile, as Figure 14 shows, the total apparent power injected from the main grid in Case 4 is also less than Case 8. Additionally, the network loss is usually concerned when the reactive power of distribution is optimally dispatched.Figure 15 depicts the hourly network losses for Cases 4 and 8.In Case 8, the main reactive power source is the wholesale market, but in Case 4, WTs and PVs integrated into ADN are also used as reactive power sources to accommodate the network demand.Thus, the reactive power regulation means Case 4 is more flexible.Accordingly, as it is seen, the total network losses for Case 4 are evidently less than Case 8.   To further illustrate the effect of reactive power support capabilities of WTs and PVs, Case 4 and Case 8 are used for analysis.In Cases 4 and 8, the proposed composite cost is both considered as the objective function.Figure 13 explores the imported total reactive power from the main grid in the optimal day-ahead scheduling framework.It can be observed that because of the reactive power support capabilities of WTs and PVs, the total reactive power injected from the main grid in Case 4 is much less than Case 8.Meanwhile, as Figure 14 shows, the total apparent power injected from the main grid in Case 4 is also less than Case 8. Additionally, the network loss is usually concerned when the reactive power of distribution is optimally dispatched.Figure 15 depicts the hourly network losses for Cases 4 and 8.In Case 8, the main reactive power source is the wholesale market, but in Case 4, WTs and PVs integrated into ADN are also used as reactive power sources to accommodate the network demand.Thus, the reactive power regulation means Case 4 is more flexible.Accordingly, as it is seen, the total network losses for Case 4 are evidently less than Case 8.

Conclusions
This paper presented an intelligent day-ahead scheduling program for a smart distribution grid.In the program, a composite economic cost objective function is established which involves the operation cost, emission cost, and network loss cost; additionally, the reactive power support capabilities of WTs and PVs are considered in the optimal energy resources scheduling.Through simulation studies, it can be observed that considering the composite cost of ADN, the total daily operation cost increases, but the emission cost and network loss cost decline more; thus, the total daily composite cost would be the minimum.Meanwhile, implementation of the reactive power support capabilities of WTs and PVs has reduced the imported reactive power from external grid and improved the economic benefit of ADN by decreasing the total daily composite cost.The total network loss of the system is decreased as well.Results certify the validity of the proposed methodology in the optimal scheduling of a smart distribution grid.

Conclusions
This paper presented an intelligent day-ahead scheduling program for a smart distribution grid.In the program, a composite economic cost objective function is established which involves the operation cost, emission cost, and network loss cost; additionally, the reactive power support capabilities of WTs and PVs are considered in the optimal energy resources scheduling.Through simulation studies, it can be observed that considering the composite cost of ADN, the total daily operation cost increases, but the emission cost and network loss cost decline more; thus, the total daily composite cost would be the minimum.Meanwhile, implementation of the reactive power support capabilities of WTs and PVs has reduced the imported reactive power from external grid and improved the economic benefit of ADN by decreasing the total daily composite cost.The total network loss of the system is decreased as well.Results certify the validity of the proposed methodology in the optimal scheduling of a smart distribution grid.

Figure 1 .
Figure 1.Single line diagram of the IEEE-33 bus test system.

Table 1 .
DGs' characteristics and cost data.

Figure 1 .
Figure 1.Single line diagram of the IEEE-33 bus test system.

Figure 4 .
Figure 4. Forecasted power output of WTs and PVs.

Figure 5 .
Figure 5. Energy resources scheduling results for Case 1.

Figure 7 .
Figure 7. Energy resources scheduling result for Case 3.

Figure 8 .Figure 6 . 17 Figure 6 .
Figure 8. Energy resources scheduling result for Case 4. Table 4 depicts the optimal hourly operation costs, emission costs, and active power loss costs of the distribution network for Cases 1-4.As it is seen, the total operation cost in Case 1 ($73,229.59),emission cost in Case 2 ($12,879.09),and network loss cost in Case 3 ($983.94)are, respectively, minimum for the individual cost in Cases 1-4.The corresponding individual costs in Case 4 are all more than them; but the composite cost in Case 4 is the minimum in Cases 1-4, which is $89,415.89.It can be observed that the total operation cost in Case 4 is $214.64 higher than Case 1; but with respect to the emission cost and network loss cost, Case 4 is, respectively, $654.26 and $297.03lower than

Figure 7 .
Figure 7. Energy resources scheduling result for Case 3.

Figure 8 .Figure 7 .
Figure 8. Energy resources scheduling result for Case 4. Table 4 depicts the optimal hourly operation costs, emission costs, and active power loss costs of the distribution network for Cases 1-4.As it is seen, the total operation cost in Case 1 ($73,229.59),emission cost in Case 2 ($12,879.09),and network loss cost in Case 3 ($983.94)are, respectively, minimum for the individual cost in Cases 1-4.The corresponding individual costs in Case 4 are all more than them; but the composite cost in Case 4 is the minimum in Cases 1-4, which is $89,415.89.It can be observed that the total operation cost in Case 4 is $214.64 higher than Case 1; but with respect to the emission cost and network loss cost, Case 4 is, respectively, $654.26 and $297.03lower than

Figure 6 .
Figure 6.Energy resources scheduling results for Case 2.

Figure 7 .
Figure 7. Energy resources scheduling result for Case 3.

Figure 8 .Figure 8 .
Figure 8. Energy resources scheduling result for Case 4. Table 4 depicts the optimal hourly operation costs, emission costs, and active power loss costs of the distribution network for Cases 1-4.As it is seen, the total operation cost in Case 1 ($73,229.59),emission cost in Case 2 ($12,879.09),and network loss cost in Case 3 ($983.94)are, respectively, minimum for the individual cost in Cases 1-4.The corresponding individual costs in Case 4 are all more than them; but the composite cost in Case 4 is the minimum in Cases 1-4, which is $89,415.89.It can be observed that the total operation cost in Case 4 is $214.64 higher than Case 1; but with respect

Case 5 :
Figures 9-12 show the energy resource scheduling for Cases 5-8.It can been seen that in Case 5, DGs and RLs are, respectively, switched on during peak hours 12-20 and 11-20; in Cases 6 and 7, DGs and RLs are switched on during all hours, as well as in Cases 2-3; in Case 8, DGs and RLs are, respectively, scheduled during hours 11-21 and 9-22.The optimal economic costs associated with Cases 5-8 are depicted in Table6.As it is seen, the total composite cost in Case 8 ($90,400.01)is also the minimum compared with Cases 5-7.Thus, in Cases 5-8, which do not consider the reactive power support capabilities of WTs and PVs, the proposed composite cost objective function is also effective and achieves the minimum composite economic cost of the distribution network.Comparing Tables4 and 6four comparisons between Cases 1 and 5, Cases 2 and 6, Cases 3 and 7, and Cases 4 and 8 are developed, it can be found that considering the reactive power support capabilities of WTs and PVs, the total daily composite cost in Cases 1-4 are $90,152.54,$96,958.74,$96,985.14 and $89,415.89which are, respectively, $934.77,$833.78,$866.66 and $984.12lower than Cases 5-8.Therefore, the conclusion can be drawn that implementation of the reactive power support capabilities of WTs and PVs can effectively decrease the economic cost of the network.

Figure 9 .
Figure 9. Energy resources scheduling results for Case 5.

Figure 10 .
Figure 10.Energy resources scheduling results for Case 6.

Figure 11 .
Figure 11.Energy resources scheduling results for Case 7.

Figure 10 .
Figure 10.Energy resources scheduling results for Case 6.

Figure 11 .
Figure 11.Energy resources scheduling results for Case 7.

Figure 10 .
Figure 10.Energy resources scheduling results for Case 6.

Figure 11 .
Figure 11.Energy resources scheduling results for Case 7.

Figure 10 .
Figure 10.Energy resources scheduling results for Case 6.

Figure 11 .
Figure 11.Energy resources scheduling results for Case 7.

Figure 13 .
Figure 13.Imported total reactive power from the main grid in Cases 4 and 8.

Figure 14 .Figure 13 .
Figure 14.Imported total apparent power from main grid in Cases 4 and 8.

Figure 13 .
Figure 13.Imported total reactive power from the main grid in Cases 4 and 8.

Figure 14 .
Figure 14.Imported total apparent power from main grid in Cases 4 and 8.

Figure 14 .
Figure 14.Imported total apparent power from main grid in Cases 4 and 8.

Figure 15 .
Figure 15.Active power losses of the network in Cases 4 and 8.

Figure 15 .
Figure 15.Active power losses of the network in Cases 4 and 8.
vectors of the active and reactive power from the whole market; Ñ are index and set of time intervals, respectively; N S , N DG , N WT , N PV and N RL are respectively sets of substations, DGs, WTs, PVs, and RLs.ρ Gird is the day-ahead wholesale electricity price; SU DG , SD DG are DG start-up and shut-down costs; X, Z are binary variables for DG start-up and shut-down status; a, b, c are fuel cost function coefficients of DGs; ρ DG , ρ WT , ρ PV are maintenance cost coefficients of DGs, WTs, and PVs; P WT , P PV are active power output of WTs and PVs; ρ RL is the contract price of

Table 1 .
DGs' characteristics and cost data.

Table 3 .
Emission characteristics and cost data.

Table 5 .
Statistical result of different algorithms.