Modeling a Hybrid Microgrid Using Probabilistic Reconfiguration under System Uncertainties

<jats:p>A novel method for a day-ahead optimal operation of a hybrid microgrid system including fuel cells, photovoltaic arrays, a microturbine, and battery energy storage in order to fulfill the required load demand is presented in this paper. In the proposed system, the microgrid has access to the main utility grid in order to exchange power when required. Available municipal waste is utilized to produce the hydrogen required for running the fuel cells, and natural gas will be used as the backup source. In the proposed method, an energy scheduling is introduced to optimize the generating unit power outputs for the next day, as well as the power flow with the main grid, in order to minimize the operational costs and produced greenhouse gases emissions. The nature of renewable energies and electric power consumption is both intermittent and unpredictable, and the uncertainty related to the PV array power generation and power consumption has been considered in the next-day energy scheduling. In order to model uncertainties, some scenarios are produced according to Monte Carlo (MC) simulations, and microgrid optimal energy scheduling is analyzed under the generated scenarios. In addition, various scenarios created by MC simulations are applied in order to solve unit commitment (UC) problems. The microgrid’s day-ahead operation and emission costs are considered as the objective functions, and the particle swarm optimization algorithm is employed to solve the optimization problem. Overall, the proposed model is capable of minimizing the system costs, as well as the unfavorable influence of uncertainties on the microgrid’s profit, by generating different scenarios.</jats:p>

section_content: Energy generation development, increasing interests toward environmental aspects, and improving the power system reliability provide enough motivation to change the distribution systems' mode from inactive to active, and increase distributed energy resources' (DERs) penetration in the microgrid (MG) environment [1, 2]. A MG is described as a small-scale distribution grid that consists of diversified distributed generation units, energy storage systems (ESSs), and local controllable loads that typically can either be operated in islanded or grid-tied modes. A grid-connected MG is tied to the upstream utility grid through a point of common coupling (PCC), in order to exchange power. More than that, providing voltage support as an ancillary service to the main distribution network is another advantage of a MG system compared with conventional end-user units [3, 4]. Since a MG can operate autonomously, in the case of events such as brownout or blackout, it is capable of managing the operation of a localized power system [5]. 

The smart grid is a novel concept that has become apparent in recent years. It is an intelligent energy system that consists of a variety of advanced energy resources/loads such as smart appliances, energy efficient and renewable energy resources, smart meters, and electric cars [6] [7] [8]. The originality of a smart grid is that it can request to improve operations, maintenance and scheduling by assuring that each component of the ingenious electric grid have the ability of both 'talk' and 'listen'. The infrastructure of energy grids are growing to be more advanced in terms of power production, data sharing and management, automation, and communication with intricate and fully united network. Therefore, intelligent management applications and services should follow in a way that attains the system objectives corresponding to supply/demand balance, utility maximization, energy saving, and operational cost reduction [9] [10] [11]. 

Also, the concept of a demand response (DR) program has been implemented to reduce or shift the energy usage at peak times, which can reduce utility electric bills, and also provide a higher capacity factor and distribution grids security at the same time [12, 13]. Basically, DR is defined as behavioral changes in electric consumption implemented directly or indirectly by customers or end users from their regular usage/injection templates in response to determined signals. The establishment of time-differentiated pricing, such as time-of-use (TOU) and real-time pricing (RTP) and the incentive-based DR contracts that can be offered to the customers drives the paradigm shifting, in which the wholesale price fluctuations will be reflected to the end users and will motivate curtailment and responsive load shifting [14]. In emergency cases, a MG is often needed to reconfigure itself as soon as possible in order to support steady electricity to critical loads. Different types of operations, including topology switching, load shedding, and generation regulation, are discussed in the reconfiguration concept. Therefore, the problem of reconfiguration becomes a generic constrained nonlinear optimal problem [15]. 

This paper addresses optimal energy scheduling/management matters in a grid-tied MG reconfiguration as a fundamental problem in interactive smart grids operation, with consideration of load and renewable resources uncertainties. Previous research studies have discussed the scheduling of distributed generation (DG) operations and loads in MGs [2, 16]. Also, a three-stage optimal method for economic dispatch of the power system using sensitive factors, improved direct search method, and the incremental linear models is presented in [17]. However, the current study focuses on the central problem of uncertainty in the predicted load demand and renewable energy units in the presence of a battery storage system in a MG environment. MG operation and management approaches are discussed and analyzed that consider the intermittent nature of renewable energy resources by utilizing stochastic methods in order to solve the problem. A stochastic energy-planning modeling for a MG, including plug-in electric vehicles and intermittent renewable energy units, has been proposed in [18, 19] in order to minimize the power losses and operational cost. Participation of a typical MG in the electricity market under certain policies of power market using hierarchical scheduling and genetic algorithm-based methods in order to obtain the economic evaluation and reduce uncertainties of the proposed system are discussed in [16, 20]. Also, different multi-time scale control strategies that deploy electric vehicle demand flexibility in order to solve grid unbalancing and congestions [21], reduce demand charges, and effectively address renewable energies production uncertainties, are discussed [22]. 

In [23], a scenario-robust mixed integer linear programing was modeled to employ ensemble weather forecasts for the performance improvement of a hybrid MG comprising both renewable and conventional power sources. The model was exercised to quantify the privilege of utilizing ensemble weather forecasts, and the optimal performance of a hypothetical power grid, including wind turbines, is predicted by employing simulated scenarios of realistic weather forecasts based on data. The comprehensive financial dispatch problem of a MG employing four various optimization techniques comprising lambda logic, lambda iteration, direct search method (DSM), and particle swarm optimization (PSO) under different test systems and control areas is discussed in [24]. Also, [25] proposes and designs a centralized and decentralized demand-response-based multi-agent control Energies 2017, 10, 1430 3 of 17 and scheduling model, including load curtailment and a backup diesel power generation, in order to measure end-user satisfaction. 

A distributed algorithm using the classical symmetrical assignment problem is resolved in [26, 27], where energy planning is modeled as a resource allocation problem, and distribution techniques are also discussed for distributed allocation. Furthermore, in [28], a convex problem formulation is described where a dual decomposition is utilized to develop a distributed energy management strategy (EMS) in MG systems in order to fulfill the energy balance of supply/demand. 

Where the privacy constraints are associated with the linear programming approach and distributed algorithms are developed, a privacy-preserving energy planning technique in MGs can be found in [29]. Also, [30] focuses on the additive-increase and multiplicative-decrease algorithms are developed to find the optimal operation of DERs in a distributed procedure. 

In this paper, an economic dispatch (ED) problem (fuel cost and generated emission minimization) of a MG containing both dispatchable and non-dispatchable DERs is formulated and modeled with diverse operational constraints. An energy planning in MGs is formulated specifically as an optimal power flow (OPF) problem. Solving the OPF problem is a complex task due to the non-convex constraints of power flow. 

In summary, the main contributions of the paper can be expressed as follows:

• Among the formerly stated optimization techniques, PSO has been considerably applied in optimal operation planning problems, mostly because of its capability for population-based searches as well as its robustness, convergence speed, and simplicity. In this paper, the Monte Carlo (MC) simulation outcomes, along with the MG essential data, are considered as the input of a PSO algorithm, in order to find the day ahead optimal operation management of a MG, which has not been properly taken into account in the aforementioned research.
section_name: Introduction

section_content: The uncertainty of renewable energy resource and local load are considered, and the solution guarantees strength against all possible realizations of the modeled operational uncertainties.
section_name: •

section_content: The existing model applies a probabilistic reconfiguration and unit commitment (UC) concurrently to attain the optimal set points of the generating units of the MG, in addition to the MG optimal topology for the next-day electricity market. The probabilistic reconfiguration and UC for a MG optimal energy planning, with the integration of renewable energy and load demand, is a novel operation to estimate MG's profit in the uncertain situation. The proposed approach has the capability to be combined with uncertainty simulation models in order to obtain part of the required inputs. 

This paper is organized as follows: The intermittent behavior of a solar PV system and electric load consumption are considered as uncertainties of the system, which are modeled through probability density functions. In Section 2, the MG architecture and associated generating units are modeled. Proposed objective functions and related uncertainties modeling are also presented in that section. The system's technical and operational constraints are formulated is Section 3. Additionally, proposed optimal operation planning is carried out in Section 4. Simulation results and discussions on system behavior under different test conditions are presented in Section 4. Consequently, conclusions are given in the last section.
section_name: •

section_content: The application of each particular DER can create as many challenges as it may resolve; thus, a proper procedure to realize the emerging potential of DERs would address a system approach that considers generation and associated consumptions as a subsystem or a small-scale grid. Furthermore, the integration of DERs in a MG and the utilization of renewable energies in bulk amounts can enable related problems on economy, technology and environment to be attentively studied within the target system, and proper decisions to be made that improve the management of its exploitation. 

The organic municipal waste, such as food waste, agricultural waste, municipal solid waste (MSW), and municipal sewage sludge (bio-solids), have been studied in many research studies that show their high potential for bioenergy extraction [31]. In recent years, microbial fuel cells (MFCs) have been perused for the generation of electricity [32] through utilizing reactors and anaerobic reformers that use bio-waste and natural gas as the input fuel, as shown in Figure 1. The proposed MG in this research is composed of a microturbine (MT), a fuel cells (FC) unit, a solar photovoltaic (PV) system, and an battery storage system (BSS). Also, a communication network is designed among energy resources, local demand, and the MG central controller (MGCC) in order to send power references and exchange data. The control structure relies on a MGCC to deliver the energy planning, and manages the OPF in the system to fulfill the local needs. The distribution system operator (DSO) takes the MGCC as one subsection, which is qualified to control a cluster of suppliers and flexible loads locally and allow renewable energy-based resources (REBG) to supply their full advantageousness, while minimizing the emission of energy components' pollutants. The entire MG is connected to an upstream network through a PCC. Also, a smart meter device is located between the input energy and local load to record the consumption of electric energy in hourly intervals or less, and also to communicate that information back to the grid and MGCC for monitoring and billing. Smart meters are capable of bidirectional communication between the meter and the central system. The proposed MG configuration studied in this research is presented in Figure 1. 

Energies 2017, 10, 1430 4 of 17

MG in this research is composed of a microturbine (MT), a fuel cells (FC) unit, a solar photovoltaic (PV) system, and an battery storage system (BSS). Also, a communication network is designed among energy resources, local demand, and the MG central controller (MGCC) in order to send power references and exchange data. The control structure relies on a MGCC to deliver the energy planning, and manages the OPF in the system to fulfill the local needs. The distribution system operator (DSO) takes the MGCC as one subsection, which is qualified to control a cluster of suppliers and flexible loads locally and allow renewable energy-based resources (REBG) to supply their full advantageousness, while minimizing the emission of energy components' pollutants. The entire MG is connected to an upstream network through a PCC. Also, a smart meter device is located between the input energy and local load to record the consumption of electric energy in hourly intervals or less, and also to communicate that information back to the grid and MGCC for monitoring and billing. Smart meters are capable of bidirectional communication between the meter and the central system. The proposed MG configuration studied in this research is presented in Figure 1. In the following sections, the proposed model objective functions consisting of operational cost functions and emission reduction functions are presented, which will be solved and minimized through employing the PSO algorithm.
section_name: Microgrid Configuration

section_content: One of the optimization goals of the MG utilization is to minimize operational costs. The total cost of the system is calculated by adding up the costs of each generating unit as the following equation:

where cost ( ) FC t and cost ( ) MT t are the operational costs of the FC and MT at time t, respectively. Additionally, cost ( ) PV t, cost ( ) BSS t, and cost ( ) grid t stand for the cost of solar PV, battery storage, and the cost of purchasing electricity from the main grid at time t, respectively. The cost of each generating resource is calculated based on its corresponding output power per hour. The considered operation cost includes the cost of generating power, fuel costs, and the operation and maintenance (O&M) costs of energy units. 

The cost of generating energy by a MT can be calculated using the following equation: In the following sections, the proposed model objective functions consisting of operational cost functions and emission reduction functions are presented, which will be solved and minimized through employing the PSO algorithm.
section_name: First Objective Function: Operational Costs

section_content: One of the optimization goals of the MG utilization is to minimize operational costs. The total cost of the system is calculated by adding up the costs of each generating unit as the following equation:

where cost FC (t) and cost MT (t) are the operational costs of the FC and MT at time t, respectively. Additionally, cost PV (t), cost BSS (t), and cost grid (t) stand for the cost of solar PV, battery storage, and the cost of purchasing electricity from the main grid at time t, respectively. The cost of each generating resource is calculated based on its corresponding output power per hour. The considered operation cost includes the cost of generating power, fuel costs, and the operation and maintenance (O&M) costs of energy units. 

The cost of generating energy by a MT can be calculated using the following equation:

where B G,MT, P MT, and u MT are the MT generation costs in $/kW, the produced power in kW, and the on and off status of the MT, respectively. It is assumed that u MT is one when the turbine is on, and is zero when the turbine is off. Also, the microturbine start-up cost can be calculated by the below equation:

where χ MT is the cost of turning the MT on and off in $. In addition, the cost of O&M by utilizing a micro-gas turbine can be presented as follows:

In the above equation, γ MT stands for the O&M cost of a MT unit in $/kW. In general, the total operation cost of a MT system can be calculated as in Equation ( 5 ). 

Additionally, the operational cost of a FC unit is presented:

where B G,FC, P FC, and u FC are the FC generation cost in $/kW, produced power in kW, and the on and off status of the FC, respectively. The u FC is one when the fuel cell is on, and it is considered zero when the unit is off. 

The FC start-up cost is also similar to the MT start-up cost equation, and it can be shown in Equation (7). 

cost

where χ FC is the cost of turning the FC on and off in $. Similarly, the O&M cost of a FC unit can be obtained using Equation (8). 

where γ FC stands for the O&M cost of a FC unit in $/kW. Overall, the total operation cost of a FC system can be obtained as shown in Equation ( 9 ). 

cost

The cost of power production by a solar PV system is presented as following. 

where B G,PV, P PV and u PV stand for the PV generation cost in $/kW, produced power in kW, and the on and off status of the PV, respectively. Also, the O&M cost of a PV system can be shown as in Equation (11). cost O&M,PV (t) = γ PV • P PV (t) (11) where γ PV represents the O&M cost of a PV unit in $/kW. In general, the operation cost of a solar PV system can be calculated as follows:

The operation cost of a BSS system includes the cost of output power delivery by the unit and the O&M cost of the battery energy storage system. The delivery power cost at the BSS output can be calculated as follows:

cost

where P BSS, B BSS, and u BSS are the battery storage produced power in kW, the cost of power delivery at the BSS output in $/kW and the on and off status of the battery sets (which can be either one or zero), respectively. Moreover, the O&M cost of a BSS can be obtained using the below equation:

where γ BSS stands for the O&M cost of a battery storage unit in $/kW. Overall, the total operation cost of a BSS is presented as follows:

The cost of purchasing electrical power from the upstream external grid is shown in Equation ( 16 ). 

cost grid (t) = P buy-grid (t)

where P buy-grid and B buy-grid represent the amount of purchased power from the main grid in kW, and the power unit price of electricity in $/kW, respectively.
section_name: First Objective Function: Operational Costs

section_content: The amount of greenhouse gases (GHG) generated by energy resources is another factor that can be minimized in order to optimize the operation of a MG. Thus, the second objective function proposed in this paper is expressed as follows:

where emission FC (t), emission MT (t), and emission BSS (t) are the amount of emission produced by FC, MT, and BSS during the operation time, respectively. The pollution generated in the environment by the non-renewable resources are mainly three GHGs including carbon dioxide (CO 2 ), nitrogen dioxide (NO x ) and sulfur dioxide SO 2. The quantity of emissions produced by operating a FC unit, a MT, and BSS, respectively, can be calculated using the below equations [33] :

emission

emission

emission

where ξ stands for the amount of emissions generated by each conventional DER in kg/kWh, and can be expressed as follows:

In Equation ( 21 ), CO 2 (t), SO 2 (t), and NO x (t) are the amounts of carbon dioxide, nitrogen oxides, and sulfur dioxide emissions produced by units at the time t, respectively. 

Energies 2017, 10, 1430 7 of 17
section_name: Second Objective Function: Generated Emission

section_content: Renewable DG resources i. e., PV or wind turbine are non-dispatchable, and their output power is a function of the availability of the primitive sources (such as solar irradiance or wind speed etc. ). Therefore, the output of such units is assumed to be random variables because of the stochastic nature of renewable DERs, and no generation cost would be considered for their production during the operation period. In this paper, the uncertainty is modeled utilizing MC simulation, where for uncertain inputs, i. e. solar radiation and load demand, some scenarios are created. Subsequently, the system is studied under obtained stochastic scenarios as deterministic inputs. Thus, there are different statuses, which are analyzed by using different scenarios. In order to calculate scenario aggregation, the expected value (h) is employed. 

where ρ s is the probability of each scenario, h s is the value of variables in each scenario, and N s is number of scenarios. Also, the coefficient of variance (Cv x ) is defined as follows. 

where σ x and µ x are the standard deviation and mean value of parameter x, respectively. If Cv x is less than a specific tolerance, the result will be relatively acceptable. Since solar power output has an intermittent nature and the solar radiation has a high degree of uncertainty, the sun radiation needs to be presented as a probability function. Variation in solar radiation is dependent on diverse factors such as weather conditions, the time of the day, month and season; the panels' tilt angle, or the orientation of the PV solar cells to the sun radiation. The probability density function (PDF) of solar radiation is modeled using the beta density as the following [34] :

where S and Γ represent solar radiation in W/m 2 and the gamma distribution function, respectively. Also, α β, and β β are parameters of the beta probability distribution function in kW/m 2, and calculated utilizing the mean (µ s ) and standard deviation (σ s ) of solar irradiance s as following:

In addition, the behavioral patterns of various electric consumers lead to an alteration in the load demand profile in MG systems. By employing statistical analysis, the created alterations can be attained. Since electric consumption changes continually with a high level of uncertainty, in this paper, load demand is considered as a normal distribution function with a mean value and a standard deviation as follows [35]. 

where σ d and µ d are the mean and standard deviation of the load, respectively. Also, P d represents the power demand. In this paper, the uncertainty is modeled utilizing a MC simulation so that for uncertain inputs, i. e. solar radiation, temperature, and load demand, numerous scenarios are created. Therefore, the system is studied under achieved stochastic scenarios as deterministic inputs. Therefore, there are various statuses, which are analyzed by using distinctive scenarios. In fact, the MC technique is an approach that utilizes stochastic variables and their related PDFs to solve the uncertainty problems.
section_name: Uncertainties in Solar Radiation and Load

section_name: Modeling Constraints

section_content: The first and main constraint that is a required consideration in OPF is the power balance between supply and demand as following: P L (t) = P PV (t) + P MT (t) + P FC (t) + P ESS (t) + P grid (t) (28)
section_name: Power Balancing Constraint

section_content: In order to stabilize system operation, the power output of each dispatchable unit at time t should be between two max and min boundaries, as presented in Equation (29). Also, each conventional generating unit such as a MT has a minimum up/down time limit. Once the generating unit is switched on, it has to operate continually for a certain minimum time before switching it off again, as is expressed in Equations ( 30 ) and (31). Moreover, the number of stops and starts of each non-renewable DER should be less than a specific number, as indicated in Equation ( 32 ) [36]. 

where P min MT and P max MT are the upper and lower limits of MT output power in kW. Furthermore, T o f f t-1,j, T on t-1,j are the MT on and off time; additionally, u MT and u MT-1 denotes for the unit on/off status, which could be 0 or 1. Also, MUT MT /MDT MT are the minimum up/down time limits of a MG, ϑ MT,start/stop represent number of starts and stops of a MT, and N max is considered to be a certain number provided by manufacturers.
section_name: Micro-Turbine Constraints

section_content: Similarity, the constraints applied for a FC unit are presented in Equations ( 33 )- (36). 

where the parameters and variables definitions are similar to the previously mentioned MT unit.
section_name: Fuel Cell Constraints

section_content: The output power of the installed PV arrays in the studied MG during operational mode is restricted between upper and lower abounds as follows:

where P min PV and P max PV denote PV output power minimum and maximum limits, respectively. 

Energies 2017, 10, 1430 9 of 17
section_name: Photovoltaic Contraint

section_content: Applied constraints on BSS are formulated in this section. Technically, the power output of a battery set at time t is limited by lower and upper limits, as in Equation (38). Further, the battery state-of-charge (SOC) barricades the unit from overcharging or undercharging. Thus, placing constraints on the battery SOC can improve the ESS performance and lifetime. The corresponding constraint is formulated in Equation (39) so that the battery SOC should be within a given range at any hour. Furthermore, battery-ageing constraints are represented by a battery state-of-health (SOH) indicator, which is restricted by Equation (40). 

In Equation (38), ESS is in discharging mode when P ESS is positive, is in charging mode when P ESS is negative, and when P ESS has zero value, the ESS has no power generation. Likewise, SOC max, SOC min, and SOH min stand for a battery's maximum and minimum state of charge, and the unit's minimum state of health, respectively.
section_name: Battery Storage System Constraint

section_content: In the proposed model, the MG operates in on-grid mode; hence the MG is capable of injecting power into the main grid. Due to this power injection, the MG owner will take profit from the external gird. Power exchange between MG and the main utility grid also can be constrained. When the electivity that is purchased by the MG from the main grid P grid is positive, it has to be less than the grid maximum power. On the contrary, P grid is negative when the MG sells electricity to the utility, and it should be less than the grid minimum power, as follows:

P grid-sold (t) ≤ min(P grid ) (42)
section_name: Upstream Utility Grid Constraint

section_content: One of the most challenging issues that a MG operator faces is how to raise the profit. Two methods are available to achieve this purpose: unit UC and reconfiguration. UC technique has been studied in many research studies [37] ; however, the reconfiguration of this target has not been considered that much. Reconfiguration brings more benefits for the MG by changing the topology. Topology changes have advantages such as line allocation with high power transfer capacities and loss reduction. PV power and load consumption uncertainties are taken into account here. In the MG optimal operation, the decision vectors are considered to be the output power of MT (P MT ), the output power of FC (P FC ), the charge and discharge power of battery set (P BSS ), and the power exchange with the MG and external grid (P grid ) at each hour in a day. Thus, the system has a decision vector of four decision variables, and 96 decision variables for the following day that need to be determined. The decision matrix utilized in the modeling is presented in Equation (43).
section_name: Optimization Problem Formulation

section_content: PSO is introduced as a population-based search technique where each particle proceeds by suitably updating its corresponding velocity based on its own intelligence and the experience of neighboring particles. The search technique of PSO uses just a single velocity updating procedure, which includes both intensification and diversification strategies. The stepwise process of PSO is shown as follows.
section_name: System Methodology and PSO Algorithm

section_content: Allocate the total number of particles in the computational process. Set the initial velocity of the particle equal to zero. For the initial iteration t, the local memory scheme stores all of the values of initial particles called P best. The best value with minimum fitness among the set of initial particles is stored in the global memory arrangement called G best [38].
section_name: Particles Initialization

section_content: Each particle in a set is introduced by two vectors: velocity and position. For each movement at the next iteration t + 1, the particles move randomly within the search range by adding a new velocity to each of its present positions. The new velocity relating to the jth position of the ith particle and is calculated as follows:

where V i and x i are the velocity and position factor of particle i, respectively. Also, c 1 and c 2 are the acceleration constants in the range of [0, 2], and r 1 and r 2 are uniform random values in the range of [0, 1]. Additionally, k and t stand for the number of iterations and present iteration number, respectively. Also, ω represents the inertia weight factor.
section_name: Position Modification

section_content: The steps of the PSO algorithm can be described as following:

Step 1 Read system data and parameters of PSO: These include plant characteristics, iteration number, particle number, population size, cognitive and social parameters, and the constriction factor. 

Step 2 Initialization: The position and velocity of each particle is initialized randomly. 

Step 3 Fitness function evaluation: Using the objective function, the position of each particle is evaluated. The fitness function value (ft t+1 ) of each modified particle is calculated. 

Step 4 Memory update by determining P best and G best: If ft t+1 (i) < ft t (i), the ith particle related to the iteration t is replaced by the current ith particle in iteration t+1 in local memory P best (i). 

Afterwards, the best value with minimum fitness among the revised set of P best particles is stored in the global memory scheme as G best. 

Step 5 Update position and velocity: Using Equations ( 44 ) and (45) updates the velocity and position of all particles. 

Step 6 Termination point: When the maximum regeneration is reached, the process is terminated. 

If the termination process is not satisfied, it goes back to step 3. This process is repeated till the stopping criteria are reached. 

By implementing the proposed non-linear problem in the general algebraic modeling system (GAMS), the mixed integer linear program (MILP) problem is solved using Cplex 12 on a personal computer (2. 6 GHz with 8 GB of RAM). GAMS/Cplex is a GAMS solver that allows users to combine the high-level modeling capabilities of GAMS with the power of Cplex optimizers, which is a strong MILP solver and Cplex Callable Library based on the IBM ILOG Cplex.
section_name: Algorithm Procedure

section_content: In this section, the numerical and graphical results of proposed models implemented on a MG, including PV, battery, MT, and FC optimal output powers, are presented to demonstrate the feasibility of the proposed algorithm. In order to model the uncertainties in the next-day energy scheduling, first, MC simulation and hourly mean and standard deviation of solar radiation and load demand is used, and the stochastic scenarios are developed. Consequently, the expected values of the scenarios will be calculated. After meeting the convergence condition of obtained expected values, certain values of parameters related to solar radiation and MG load demand uncertainties are determined. 

The percentage values of rated output power of solar PV arrays that have been used in the system model to utilizing MC simulation in the presence of solar radiation at each hour are shown in Table 1. The MG average load demand results extracted from MC simulation over 24 h are presented in Table 2. This table has two peaks; the smaller peak occurs in the morning, and the larger peak arises in the afternoon around 6 pm. The maximum and minimum power production limits of the DERs, as well as the valid range of power exchange with upstream utility, are given in Table 3. The battery efficiency of charge is assumed to be 85%, and the efficiency of discharge is considered 90%. The system cost of start-up/shut-down in $, power generation cost in $/kWh, and the operation and maintenance cost in $/kWh of each generating unit are presented in Table 4. The real-time electricity market prices per each hour in $/kWh are mentioned in Table 5. The variation between the minimum and maximum price rates is around $1. 9, where the peak time occurs around noon, and the off-peak time happens during the nighttime. In addition, the amount of GHGs produced by distributed generation units in kg/MWh during the operation period of the MG are observed in Table 6. The MC simulation results, along with the MG required data, are considered as the input of PSO algorithm in order to find the next-day optimal operation scheduling of the MG. The output results of Energies 2017, 10, 1430 13 of 17 the PSO algorithm are the values of decision variables, which are the hourly power outputs of the MT and FC, battery charge and discharge power, and the amount of power exchange with the external grid. Also, the algorithm calculates the system total operation cost and the amount of generated emission by DERs in each hour and each day. Tables 7 and 8 show the MG hourly optimal scheduling results, system operation cost, and produced emissions over the next 24 h, respectively. Figure 2 demonstrates the daily energy consumption of the system vs. the load demand after PV generated power subtraction. It is clear from the results that the PV energy production during the daytime can smoothen the load profile, which is beneficial in terms of system cost and emission. The contribution of the remaining DERs will fulfill the MG base load. 24 6. 0000 25. 8402 0 -9. 5402 30. 0000

Figure 2 demonstrates the daily energy consumption of the system vs. the load demand after PV generated power subtraction. It is clear from the results that the PV energy production during the daytime can smoothen the load profile, which is beneficial in terms of system cost and emission. The contribution of the remaining DERs will fulfill the MG base load. The optimal power outputs of the MG's DGs, and also the power purchased/sold to the external grid in a typical test day, are shown in Figure 3. As included in the results, the MT and FC are on during the test day to meet the base load. The output of the FC and MT decreases during noontime, and around solar noon, when the PV generation is high. Additionally, there is no need to purchase power from the grid at that time, and there is even some excess generated energy that can be sold to the utility. Since electricity prices are low during the morning and very late at night, it is beneficial for the MG to purchase power from the external grid on its maximum boundary. BSS is in discharging mode when P BSS > 0 and charging mode when P BSS < 0; also, BSS has zero production when P BSS = 0. Thus, the BSS is in charging cycle when the electricity price is low at the beginning of the day, and it is in discharging cycle in the middle of the day when the electricity market price is high, and BSS is required to fulfill some part of the load. The BSS will go back to the initial condition at the end of the day. The optimal operation results in Figure 3 clearly demonstrate that the optimal system behavior and reaction to the various changes such as local power consumption, real-time electricity price, and weather condition. The optimal power outputs of the MG's DGs, and also the power purchased/sold to the external grid in a typical test day, are shown in Figure 3. As included in the results, the MT and FC are on during the test day to meet the base load. The output of the FC and MT decreases during noontime, and around solar noon, when the PV generation is high. Additionally, there is no need to purchase power from the grid at that time, and there is even some excess generated energy that can be sold to the utility. Since electricity prices are low during the morning and very late at night, it is beneficial for the MG to purchase power from the external grid on its maximum boundary. BSS is in discharging mode when PBSS > 0 and charging mode when PBSS < 0; also, BSS has zero production when PBSS = 0. Thus, the BSS is in charging cycle when the electricity price is low at the beginning of the day, and it is in discharging cycle in the middle of the day when the electricity market price is high, and BSS is required to fulfill some part of the load. The BSS will go back to the initial condition at the end of the day. The optimal operation results in Figure 3 clearly demonstrate that the optimal system behavior and reaction to the various changes such as local power consumption, real-time electricity price, and weather condition. Moreover, net system operation costs and produced emissions are presented in Figure 4. It is concluded from the figure that the system cost has its peak around solar noon, when all of the MG's DERs are on, and no power is injected from the grid. Although, the amount of generated emission Moreover, net system operation costs and produced emissions are presented in Figure 4. It is concluded from the figure that the system cost has its peak around solar noon, when all of the MG's DERs are on, and no power is injected from the grid. Although, the amount of generated emission has its off-peak at the same time, because a major part of the load is covered by PV at the mentioned period. In contrast, the MG's total cost is low at the beginning and end of the day, when the upstream grid meets the majority of the load, and the amount of generated GHGs is on its peak at the same time. 

Energies 2017, 10, 1430 15 of 17 has its off-peak at the same time, because a major part of the load is covered by PV at the mentioned period. In contrast, the MG's total cost is low at the beginning and end of the day, when the upstream grid meets the majority of the load, and the amount of generated GHGs is on its peak at the same time. Based on Table 8, a daily system total cost would be $1548. 8, and the total generated emission during that day would be 575. 6 kg. The optimal results can be compared with the case when the total power demand is purchased from the upstream utility with the electricity price rate, which is presented in Table 5. Also, the generated emission rate in kg per MWh by the utility grid can be extracted from Table 9. Utilizing all required data, the total operation cost when the whole load demand of the MG is fulfilled by the main grid is calculated and is equal to $2240. 6; also, the corresponding generated GHGs emission would be 1584 kg. As it is observed from the results, providing electric energy through the MG utilizing renewable and non-renewable DERs would be a cost-saving option. In addition to the economic benefits, operating a MG that includes renewable resources is an environmentally-friendly choice that would decrease the amount of GHGs in the environment. In this research, the generated emission operating a MG reduced the system emission by one third compared with the conventional power supply.
section_name: Simulation Results

section_content: In order to find the solutions of optimally DER operations in a MG, a novel optimization technique has been proposed in this paper, where the modeling goals were minimizing MG operation costs as well as maintenance and emission costs. The mathematical modeling and corresponding constraints of each MG's components were introduced in order to solve the proposed objective function. Utilizing the presented algorithm along with considering solar PV power generation and electrical load uncertainties, a day-ahead power generation scheduling of DERs was calculated. The contribution of the proposed model focuses on the adjustment of the uncertain parameters inside the approach, which are addressed by utilizing MC simulation in contrast to conventional techniques. 

The simulation results have shown that high renewable energy sources penetration will have a significant total effect on the grid operation, especially when emission reduction is one of the system objectives. However, the MG total operation cost rises correspondingly in the investigated time Based on Table 8, a daily system total cost would be $1548. 8, and the total generated emission during that day would be 575. 6 kg. The optimal results can be compared with the case when the total power demand is purchased from the upstream utility with the electricity price rate, which is presented in Table 5. Also, the generated emission rate in kg per MWh by the utility grid can be extracted from Table 9. Utilizing all required data, the total operation cost when the whole load demand of the MG is fulfilled by the main grid is calculated and is equal to $2240. 6; also, the corresponding generated GHGs emission would be 1584 kg. As it is observed from the results, providing electric energy through the MG utilizing renewable and non-renewable DERs would be a cost-saving option. In addition to the economic benefits, operating a MG that includes renewable resources is an environmentally-friendly choice that would decrease the amount of GHGs in the environment. In this research, the generated emission operating a MG reduced the system emission by one third compared with the conventional power supply.
section_name: Conclusions

section_content: In order to find the solutions of optimally DER operations in a MG, a novel optimization technique has been proposed in this paper, where the modeling goals were minimizing MG operation costs as well as maintenance and emission costs. The mathematical modeling and corresponding constraints of each MG's components were introduced in order to solve the proposed objective function. Utilizing the presented algorithm along with considering solar PV power generation and electrical load uncertainties, a day-ahead power generation scheduling of DERs was calculated. The contribution of the proposed model focuses on the adjustment of the uncertain parameters inside the approach, which are addressed by utilizing MC simulation in contrast to conventional techniques. 

The simulation results have shown that high renewable energy sources penetration will have a significant total effect on the grid operation, especially when emission reduction is one of the system objectives. However, the MG total operation cost rises correspondingly in the investigated time duration. Additionally, the proposed approach can facilitate the choice for the grid-connected MG operator to carry out a stable and economical operation at the planning decision stage for the next day. In consequence, it would be favorable for both system goals when a convenient means of energy exchange was provided between the MG and the upstream utility in the grid-tied mode.
section_name: Conclusions

section_name: Conflicts of Interest:

section_content: Author Contributions: Hadis Moradi, Amir Abtahi, and Ali Zilouchian conceived and developed the research problem; Hadis Moradi wrote the entire manuscript; Mahdi Esfahanian supervised the paper writing; Amir Abtahi and Ali Zilouchian participated in paper revisions.

section_content: The authors declare no conflict of interest.
section_name: Conflicts of Interest: