A Three-Stage Optimal Approach for Power System Economic Dispatch Considering Microgrids

Wei-Tzer Huang 1,*, Kai-Chao Yao 1, Chun-Ching Wu 1, Yung-Ruei Chang 2, Yih-Der Lee 2 and Yuan-Hsiang Ho 2 1 Department of Industrial Education and Technology, National Changhua University of Education, Changhua 500, Taiwan; kcyao@cc.ncue.edu.tw (K.-C.Y.); wugch001@yahoo.com.tw (C.-C.W.) 2 The Institute of Nuclear Energy Research, Taoyuan City 325, Taiwan; raymond@iner.gov.tw (Y.-R.C.); ydlee@iner.gov.tw (Y.-D.L.); twingo_ho@iner.gov.tw (Y.-H.H.) * Correspondence: vichuangl@cc.ncue.edu.tw; Tel.: +886-939-828-628


Introduction
Energy efficiency, economic dispatch (ED), and greenhouse gas emissions are major issues associated with utilities.Although electrical energy has been used for over one hundred years since the first electric network was established in 1882 at the Pearl Street Station in New York by Thomson Edison [1], the development of power systems has evolved into generation, transmission, and distribution with a vertical system architecture.In a vertical system structure, a generation system is composed of major large centralized thermal power plants, minor large reservoir hydroelectric power plants, and a small portion of renewable energy generation.Carbon dioxide and other greenhouse gas emissions caused by thermal power plants exacerbate global warming because of the usage of carbon-based fuels in such plants.The electrical energy produced is delivered through transmission and distribution systems until it reaches customers.Electrical energy losses during transmission are considerable.Such losses not only affect the environment, but also reduce the overall system efficiency.Consequently, the traditional large centralized power system architecture is impractical in terms of high efficiency, energy savings, and environmental protection.Conversely, Distributed Energy Resources (DERs) composed of distributed generations (DGs) and energy storage devices [2] are incorporated Energies 2016, 9, 976 2 of 18 in local high-or low-voltage distribution networks to satisfy the load demand.As such, DERs can reduce the losses caused by power delivery, and renewable energy-based wind turbine, solar power, or relatively clean units, such as fuel cells and microturbines, can be utilized to reduce greenhouse gas emissions and environmental impact.
DERs interconnected into an existing power grid to form an area grid, which can be operated in grid-tied or autonomous modes, are called MGs [3].The power grid integrated with information and communication technologies is a smart grid.The development of low-carbon central generation and local energy networks has emerged as a new trend.Europe, the USA, Japan, and Canada have been actively involved in pilot system construction and related research.For instance, the European SmartGrids Technology Platform for future power grids began its work in 2005.This platform aims to formulate and promote a vision for the development of electricity networks in Europe in 2020 and beyond [4].Future power systems will be composed of low-emission power plants, solar energy, wind turbines, biomass, ocean wave energy, hydroelectric, microturbines, gas turbines, fuel cells, and energy storage devices.As a result, the operation of the future electricity network will be shared between central and distributed generators.Thus, the control of distributed generators can be aggregated to form MGs or virtual power plants (VPPs) to facilitate their integration in physical systems and in markets.Overall, the planning, design, control, and operation of MGs are essential, especially the on-line management or optimal dispatch of various types of DG in an MG under grid-tied and autonomous operation modes.The optimal dispatch of future power grids integrated with MGs as VPPs for safety, economic, and energy loss reduction planning and operation is also relevant in power system research.
ED problems in traditional power systems should be solved in terms of the total minimum operation costs by adjusting the output of each generator to satisfy the load demand and related operational constraints.Since 1920, several power engineering experts and scholars have conducted relevant research.ED problems are based on circuit and optimal theories.The primary task aims to obtain real-time system states, and the power flow algorithm is applied to determine the exact bus voltages and line flows of power systems via traditional AC power flow method [5][6][7][8].Although solutions generated by such methods are precise, they require time-consuming processes and iterations and exhibit a convergence problem.Newton-Raphson method and Gauss-Seidel method are not applicable for the real-time calculation of power flows.The methods of network sensitivity factors, such as generation shift distribution factor [9], generalized generation shift distribution factor [10], Z-bus distribution factor [11], and Jacobian-based distribution factor (JBDF) [12], have been proposed to improve the computation speed and to eliminate the need for iterations.The application of these approaches overcomes the weakness of Newton and Gauss-based power flow algorithms in ED.Common optimal algorithms in ED and unit commitment are based on Lagrange multiplier-based methods [13,14] and meta-heuristic approaches [15].Soroudi and Rabiee [16] proposed a new model to solve the multi-area dynamic ED problem by considering the uncertainties in wind power generation, energy prices, and system demands.Optimality condition decomposition is also utilized to verify the proposed approach for the real-time operation of practical power systems.Niu et al. [17] presented an improved quantum-behaved particle swarm optimization (QPSO) to solve the ED problems with valve-point effects and multiple fuel options.This method connects the basic QPSO [18] with a selective probability operator to preserve the fast convergence characteristics and to reduce suffering from premature convergence.Cheng and Zhang [19] combined a dynamic economic dispatch model and an improved particle swarm optimization (PSO) approach to solve the ED problem considering the uncertain nature of wind generation.Li et al. [20] proposed a Quasi-Newton method to update the multiplier, together with an initialization strategy for the multipliers and the approximation matrix.It can be used in dynamic economic dispatch to large-scale integration of renewable generation in power systems.Zaman et al. [21] presented a self-adaptive differential evolution and a real-coded genetic algorithm to solve the dynamic dispatch problem, simulation results demonstrated that the approach is with high solution quality and reliability.Liu et al. [22] introduced an inexact mix-integer two-stage programming model to address low-carbon energy systems in terms of costs, energy security, and carbon dioxide emission.In a previous study [23], fruit fly algorithm is applied to solve an optimal dispatch problem in power systems with a VPP.This algorithm is a nonlinear mixed-integer programming with inter-temporal constraints, and the proposed approach can manage the optimal dispatch strategy among DGs in a VPP and a trading plan with electricity market.Although, these algorithms provide solution quality, convergence speed, and robustness depending on their applications.However, all the related studies mentioned above seldom consider the distribution substation-level MGs in ED problems in modern power systems.This is vital because of the inclusion of MGs in power systems will significantly affect the power dispatch due to the large magnitude of the import and export power flows.Therefore, in this study, an innovative three-stage approach is proposed for the ED of modern power systems comprising MGs as VPPs.The ED model that is derived by the network sensitivity factors in main power grid is with high efficiency and without iteration problem after load changes; besides, the improved direct search method (IDSM), which is with easy and fast characteristics is proposed for optimal dispatch in MGs.Finally, the power exchanges between the main power grid and MGs are modeled by linear incremental model for solving the entire system ED problem by linear programming method.Consequently, this new approach can be utilized for smart grid applications.This paper is divided into five sections.Section 1 introduces the background and objectives of this study.Section 2 describes the traditional Lagrange multiplier-based ED and optimal dispatch in MGs and then derives the proposed three-stage ED model.Section 3 explains the solution procedure of the proposed model.Section 4 discusses the simulation results.Section 5 presents our conclusions.

Problem Description and Formula Derivation
MGs play important roles as VPPs in power system operations.Therefore, power systems will be composed of the main power grid and several MGs.The generation dispatch corresponding to load demand will become highly complicated in the near future.In this section, the ED problem is divided into four parts, namely, system structure, ED model of the main power grid, ED model of the MG, and integrated ED model of the power system.

System Structure
Figure 1 illustrates the system structure of a power system that consists of the main power grid and N distribution substation-level MGs.In terms of transmission systems, the distribution networks with DERs and loads can be regarded as MGs at the primary side of the main transformer as the point of common coupling.According to a multi-area ED model, power systems can be partitioned into the main power grid and a MG in each area.The MG in each area is considered a VPP, which is composed of diesel engine generators, fuel cells, microturbine generators, wind turbine generators, photovoltaic systems, and loads.

ED of the Main Power Grid
The ED in the main power grid can be described as follows: subject to: Energies 2016, 9, 976 4 of 18 where ƒ i and P Gi represent the cost function and power generation of the ith thermal generator unit, respectively; P D denotes the total power demand in the main power grid; P L is the transmission loss; and P m is the real power flow in the ith transmission line.Using Lagrange multiplier and adding equality constraints, we can rewrite the optimal problem of ED as follows [24]: On the basis of Equation ( 5), we let partial L to P Gi be 0, as shown in Equation ( 6): Using Equation ( 6), we derive Equation (7): we let: equality constraints, where the generation is equal to the total demands plus losses, can be obtained, as shown in Equation (2).In Equation (7), ∂P L /∂P Gi is the incremental transmission loss (ITL), and Equation ( 7) is rearranged as: we rewritten Equation (9) as: where Pf i is equal to 1/(1 − ITL i ) corresponding to the penalty factor.Although ED problems can be solved by artificial intelligence-based optimal algorithms, the Lagrange multiplier-based methods remain the most reliable, most accurate, and fastest in practical system applications.On the basis of Equations ( 9) and (10), we observe that the solution kernel of ED depends on the rapid calculation of ITL i , Pf i , and P L , and the first task is to derive the mathematical model of P L to rapidly compute ITL i and Pf i .The formula derivation is described in the following sections.

ED in a Microgrid
The DERs included in MGs are shown in Figure 1.In terms of the multi-area ED concept, the large thermal power plants, transmission systems, and distribution substations in case of the Taipower system regarded as equivalent loads can be formed as the main power grid, and partial distribution stations whose distribution networks with DERs and loads are developed as MGs.Consequently, the primary task of the integrated ED problem in modern power systems is the optimal dispatch in MGs.Thus, the first step is to establish the cost function in an MG.
Similar to large thermal units, the fuel cost of a diesel engine generator is the quadratic polynomial function of its actual power output.However, the fuel cost of a fuel cell generation system is the linear function of its actual power output and efficiency.Therefore, the aggregate cost function of the MG is combined with nonlinear and linear functions.This problem can be solved by direct or indirect search algorithms, where the direct search algorithm called nongradient or zeroth-order method only considers objective function and the indirect search algorithm should consider the partial derivatives of objective function.In this study, an IDSM is proposed to solve the optimal dispatch problem in MGs effectively.In general, a direct search method (DSM) [25] is unsuitable for large-scale systems.The first step of this search approach is to set the incremental cost (IC) value (∆P).The unit with the minimum IC value is then selected to increase the generation until the IC of each unit is equal and thus satisfies the power balance equation.The weakness of this search method is the start point, that is, it usually starts from 0. As such, the searching time becomes lengthy.The proposed IDSM can improve the limitation of the DSM.The first step is to obtain the ICs of all units, including fuel cell and microturbine with linear costs, as shown in Equations ( 11) and ( 12): where f FCi denotes the fuel cost of the ith fuel cell; β i and P FCi represent the cost coefficient and power generation of the ith fuel cell unit, respectively.f MTi , γ i , and P MTi are the fuel cost, cost coefficient, and power generation of the ith microturbine unit, respectively.NFC and NMT represent the unit number of fuel cell and microturbine in an MG, respectively.IC can be derived as follows: The ICs of fuel cell and microturbine are constant, as shown in Equations ( 13) and ( 14).Therefore, they can be sequentially arranged as follows: where IC 1 < IC 2 < ... < IC N , and N is the summation of N FC and N MT in MG-A.The fuel cost of diesel engine generator can be expressed as follows: where f DGi denotes the fuel cost of the ith diesel engine unit; a DGi , b DGi , and c DGi are the cost coefficients.NDG represents the unit number of the diesel engine unit in an MG.The IC can be derived as follows: substituting the IC 1 of Equation (15) in Equation ( 17) results in: the initial generation, P 0 DGi , is started, as shown in Equation ( 18).P DGi is gradually adjusted according to the proposed IDSM method until this procedure satisfies the power balance equation, as shown in Equation ( 19): where P DA is the total power demand in MG-A, P DA is the real power loss in MG-A, P PV is the power generation of photovoltaic in MG-A, and P WT is the power generation of wind turbine in MG-A.

ED Integrated with Main Power Grid and MGs
In this section, the integrated model of the main power grid and MGs is derived on the basis of Sections 2.2 and 2.3.According to the multi-area dispatch concept, the completed ED problem involves three stages: Stage III: integrated ED of main power grid and MGs The ED in stage III is based on those in stages I and II.The detailed derivation in each stage to establish the linear increment model is described as follows.

Thermal Unit in Main Power Grid
The linear cost model of the thermal power unit in the main power grid is expressed as Equation (20).The coefficient of each unit depends on the power plant type:

ED Integrated with Main Power Grid and MGs
In this section, the integrated model of the main power grid and MGs is derived on the basis of Sections 2.2 and 2.3.According to the multi-area dispatch concept, the completed ED problem involves three stages: Stage III: integrated ED of main power grid and MGs The ED in stage III is based on those in stages I and II.The detailed derivation in each stage to establish the linear increment model is described as follows.

Thermal Unit in Main Power Grid
The linear cost model of the thermal power unit in the main power grid is expressed as Equation (20).The coefficient of each unit depends on the power plant type: Energies 2016, 9, 976 7 of 18

NonRenewable Unit in MGs
The linear cost models of diesel engine, fuel cell, and microturbine units in an MG are expressed as Equations ( 21)- (23), respectively:

Linear Programing Model of ED
The linearized dispatch model is established in stage III on the basis of the linear cost model described in Sections 2.4.1 and 2.4.2 and is solved by the linear programming method.The objective function is expressed as follows: subject to: where P eA is the power interchange between the main power grid and MG-A, and positive P eA means that power is fed into MG-A.These equations can be expressed as a standard linear model as follows: subject to: each extreme point is gradually searched in accordance with the linear programming method.
The optimal power generation of each unit is derived, as shown in Equations ( 34)-(37).The optimal ED solution of power systems considering MGs as VPPs can be obtained: Energies 2016, 9, 976 8 of 18 P MTi = P 0 MTi + ∆P MTi . (37)

Solution Procedure
In this section, the detailed solution algorithm and model of the main power grid and MGs are explained.The proposed approach can be applied in multi-area complex power systems.The solution procedure is described as follows.

Stage I
The JBDF [12]-based Lagrange multiplier method is proposed to solve the ED in the main power grid.A power flow program should be executed to solve the ED problem in each iteration by the traditional Lagrange multiplier method.However, this procedure is time consuming and may result in divergence condition.The key point is to obtain the ITL and penalty factor (Pf ) efficiently after the load demand and generation are changed.The JBDF-based formula is derived to obtain ITL and Pf in the ED procedure instead of running a power flow program and to overcome the weakness of ITL and Pf derivation after load demand changes.

Models of Transmission Line Loss and Line Flow
P m and Q m represent the line m flows of real and reactive powers, and the transmission line loss P L is expressed as follows: where R m denotes the resistance of line m, and NL is the number of transmission lines.The real and reactive power line flows can be expressed as [12]: in the ED procedure, the real power is generally only adjustable using generators corresponding to load demand changes, and the reactive power is kept constant.Consequently, Equations (39) and (40) are simplified as Equations ( 41) and (42), respectively.∆P Gi is the change in the real power of the ith generator from the base case: where F p and F q are the distribution factors of JBDF, which are expressed as Equations ( 44) and (45), respectively: where ∂ V p /∂P i , ∂ V q /∂P i , ∂δ p /∂P i and ∂δ q /∂P i represent the magnitude and phase angle of the bus voltage sensitive to real power injections.These terms can be calculated in the inverse Jacobian matrix of the base case power flow solution.
3.1.2.Formula of ITL and Penalty Factor The partial derivatives of Equation ( 38) and the JBDF-based ITL can be derived as Equation ( 46).Pf can be calculated by 1/(1 − ITL i ):  2) and ( 7), the model of JBDF-based ED conditions is rewritten as follows: extending Equation ( 48) to all units results in (NG + 1) simultaneous nonlinear equations in a matrix form as follows: Equation ( 49) can be solved by the Newton-based algorithm, where J 1 , J 2 , J 3 , and J 4 are the sub-matrices of Jacobian matrix, as explained in the following section.
The diagonal and off-diagonal elements of J 1 are expressed as follows: the partial derivatives of ITL are shown in Equations ( 52) and (53): the elements of J 2 are: the elements of J 3 are: the elements of J 4 are: the model of the JBDF-based ED in the main power grid can be derived according to these equations.The solution flow chart is shown in Figure 2. The solution is obtained when the convergence conditions, ∆P and ∆λ, are less than a specified accuracy, i.e., ε.On the contrary, it is divergence, the iterations will reach to maximum iteration number and then stop the iteration procedure.It is worth noting that the Newton-based algorithm is sensitive to the initial guess value.In general, the physical meaning of λ is the system incremental cost of delivered power; therefore, the initial value could be guessed in a reasonable range according to the average cost per megawatt for avoiding the unreasonable solution or divergence.
the elements of J3 are: the elements of J4 are: the model of the JBDF-based ED in the main power grid can be derived according to these equations.The solution flow chart is shown in Figure 2. The solution is obtained when the convergence conditions, P and   , are less than a specified accuracy, i.e.,  .On the contrary, it is divergence, the iterations will reach to maximum iteration number and then stop the iteration procedure.It is worth noting that the Newton-based algorithm is sensitive to the initial guess value.
In general, the physical meaning of  is the system incremental cost of delivered power; therefore, the initial value could be guessed in a reasonable range according to the average cost per megawatt for avoiding the unreasonable solution or divergence.

Start
Input system parameters for ED, such as bus ,line data, fuel cost function, cost coefficient, etc.
And then Establish real and reactive power JBDF,i.e.Fp, Kp, Fq andKq.

Yes
Output optimal dispatch results in stage I.

Stage II
The DERs are divided into nonrenewable and renewable units.The cost functions of nonrenewable units and the prediction of the available power of renewable units are essential for optimal dispatch in an MG.The characteristics of the DERs are described as follows.
For the nonrenewable units, the fuel cost of a diesel engine generator is expressed as Equation ( 16), and the total fuel cost of diesel engine units in MG-A are shown in Equation (57).Molten carbonate fuel cells and solid oxide fuel cells are used in this study.The fuel cost of the fuel cell unit is shown in Equation (11), and the total fuel cost of fuel cell units in MG-A is shown in Equation (58).Similar to the cost function of the fuel cell unit, the fuel cost of the fuel cell unit is shown in Equation ( 12), and the total fuel cost of fuel cell units in MG-A is shown in Equation (59): The power output of the wind turbine unit is expressed as Equation ( 60), where C p is the power coefficient, ν is the air density, V w is the wind velocity, and A is the rotor swept area.The total power generation of wind power units in MG-A is shown in Equation (61).The power output of the photovoltaic module can be expressed as Equation ( 62), with the measured solar irradiation (W/m 2 ) and ambient temperature ( • C), where P stc is the maximum power at the standard test condition, G ing is the incident irradiance, G stc is the irradiance at the standard test condition (1000 W/m 2 ), k is the temperature coefficient of power, T c is the cell temperature, and T r is the reference temperature.The total power generation of photovoltaic units in MG-A is shown in Equation ( 63).Therefore, the renewable units are regarded as negative power loads in optimal dispatch in MGs: C pi νV 3 w A (61) On the basis of the cost function and power generation of each unit, we can express the objective function for the optimal dispatch problem in MG-A as follows: subject to: P min DGi ≤ P DGi ≤ P max DGi , (66) Finally, the IDSM is proposed to solve the optimal dispatch in an MG, and the solution flow chart is shown in Figure 3. (68) Finally, the IDSM is proposed to solve the optimal dispatch in an MG, and the solution flow chart is shown in Figure 3.

Start
Input system parameters for optimal dispatch in the MG, such as bus data ,line data, fuel cost function, cost coefficient, available power of renewable units, etc.. Select the minimum value, IC1, and substituting IC1 into (18) to obtain .
Increasing real power of the minimum value in (15), and adjust PDG of diesel engine unit as well as its IC.

Satisfy (65)?
Output optimal dispatch results in stage II.

Yes
No Using ( 13) and ( 14) to compute IC of the unit with linear cost function, and then establish IC matrix, (15).15) and ( 17) of all units and their limits of (66)~(68).

Stage III
In stage III, the linear incremental model of the entire power system is established on the basis of the optimal dispatch solutions of the two previous stages.Considering the power interchange of MG-A, PeA, into the main power grid, we can obtain the linearly optimized model of the main power grid and MG-A.The linearly optimized model in stage III can be derived by increasing N number of MGs in the power system, as shown in Equation (69): According to the linear model and satisfying the related constraints, the optimal dispatch generations are solved by linear programming.The optimal solution of the entire power system with MGs as VPPs can also be found.

Stage III
In stage III, the linear incremental model of the entire power system is established on the basis of the optimal dispatch solutions of the two previous stages.Considering the power interchange of MG-A, P eA , into the main power grid, we can obtain the linearly optimized model of the main power grid and MG-A.The linearly optimized model in stage III can be derived by increasing N number of MGs in the power system, as shown in Equation (69): According to the linear model and satisfying the related constraints, the optimal dispatch generations are solved by linear programming.The optimal solution of the entire power system with MGs as VPPs can also be found.

Numerical Results and Discussion
An effective ED program is coded in the Matlab environment according to the proposed three-stage model and solution algorithms.The proposed approach is first used to evaluate the ED problem of the original IEEE 14-bus test system and to verify the accuracy and effectiveness of the proposed approach, and the related parameters for power flow and ED are obtained from previous studies [25,26].In this scenario, the MG is assumed to be operated autonomously and therefore not considered in ED in the main power grid for the comparison of the traditional Lagrange multiplier method and the proposed approach based on the known results.The numerical results (Table 1) demonstrate that the developed ED program is valid.The ED problem of the main power grid considering MGs is then solved by using the three-stage ED program.The IEEE 14-bus test system and an MG (Figures 4 and 5) are modified as a sample system to verify the feasibility and accuracy of the proposed approach.The equivalent lumped daily load demand curves of the main power grid are illustrated in Figure 6.The equivalent daily load demand curves that correspond to different load types at each bus are shown in Figure 7.The purchased/sold electricity cost of MG-A are depicted in Figure 8.This cost comprises the peak-load (0.2392 $USD/kW•h), half-peak-load (0.0894 $USD/kW•h), and off-peak-load (0.0447 $USD/kW•h) pricing from the main power grid during the summer season.On the basis of the three-stage ED algorithm and constraints, we show the numerical results of the optimal dispatch of each generator in the main power grid in Figure 9, and the optimal power interchange between the main power grid and MG-A is shown in Figure 10.The numerical results of the generation cost per kWh and the total generation cost of the sample system are illustrated in Figure 11a,b, respectively.The numerical results also demonstrate that the bus voltage, current flow, and power generation limits of each unit are consistent with the constraints.The executing time of the proposed approach is around 25% of the traditional Lagrange multiplier method in the IEEE 14-bus test system with a MG; furthermore, there is no divergence problem for power flow solution of the proposed approach because of no iteration after load demand changes.Therefore, the performance of both computation time and reliability of the proposed approach is superior to those of the traditional Lagrange multiplier method.Hence, the proposed approach exhibits high efficiency and rapidly solves the ED of power systems considering MGs as VPPs.

Conclusions
In this study, effective and fast three-stage ED model and approach are proposed to solve the problem of modern power systems considering MGs as VPPs.Using the traditional Lagrange multiplier method, we derive the JBDF-based ED algorithm for stage I dispatch according to the network sensitivity factors in the main power grid, and the IDSM for stage II is developed for the optimal dispatch in an MG.The linear incremental model and solution algorithm for stage III are obtained and integrated with stages I and II.For the algorithms proposed in this paper, the power systems are divided into the main power grid and area MGs.The ED problem is solved on the basis of the concept of multi-area ED by using the proposed three-stage approach.Besides, the traditional methods of ED required iterations in each solution procedure; nevertheless the proposed approach only need to execute one iteration procedure by Newton algorithm as base case and then without any iteration after load demand changes.It is absolutely verified that not only the computing speed but also the convergence are superior to other methods that need iterations after load demand changes in real-time dispatch.The numerical results of the test system demonstrate that the proposed approach exhibit high accuracy and efficiency.Therefore, this approach is helpful for optimal dispatch because of the large amounts of DERs interconnected into distribution systems.The proposed approach can also be utilized to solve the area power dispatch problems because of the distribution systems formed as MGs, which are regarded as VPPs for future smart grid applications.

Figure 1 .
Figure 1.Structure of a power system considering N distribution substation-level MGs.

Figure 2 .
Figure 2. Solution flow chart of the JBDF-based ED in the main power grid.Figure 2. Solution flow chart of the JBDF-based ED in the main power grid.

Figure 2 .
Figure 2. Solution flow chart of the JBDF-based ED in the main power grid.Figure 2. Solution flow chart of the JBDF-based ED in the main power grid.

Figure 3 .
Figure 3. Solution flow chart of IDSM in an MG.

Figure 3 .
Figure 3. Solution flow chart of IDSM in an MG.

Figure 8 .
Figure 8. Electricity purchased/sold cost of MG-A.

Figure 9 .
Figure 9. Numerical results for the optimal generation of each unit in the main power grid.

Figure 9 .
Figure 9. Numerical results for the optimal generation of each unit in the main power grid.

Figure 9 .
Figure 9. Numerical results for the optimal generation of each unit in the main power grid.

Figure 9 .
Figure 9. Numerical results for the optimal generation of each unit in the main power grid.

Figure 10 .Figure 11 .Figure 11 .
Figure 10.Power interchange between the main power grid and MG-A.

Table 1 .
Comparison of the simulation results of the traditional Lagrange multiplier method and the proposed approach for the ED of the original IEEE 14-bus test system.