1. Introduction
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 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.
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.
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.
2. Microgrid Configuration
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.
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.
2.1. First Objective Function: Operational Costs
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
and
are the operational costs of the FC and MT at time
t, respectively. Additionally,
,
, and
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
BG,MT,
PMT, and
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
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
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,
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
BG,FC,
PFC, and
are the FC generation cost in
$/kW, produced power in kW, and the on and off status of the FC, respectively. The
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).
where
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
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).
The cost of power production by a solar PV system is presented as following.
where
BG,PV,
PPV and
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).
where
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:
where
PBSS,
BBSS, and
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
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).
where
Pbuy-grid and
Bbuy-grid represent the amount of purchased power from the main grid in kW, and the power unit price of electricity in
$/kW, respectively.
2.2. Second Objective Function: Generated Emission
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
,
, and
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]:
where
stands for the amount of emissions generated by each conventional DER in kg/kWh, and can be expressed as follows:
In Equation (21), CO2(t), SO2(t), and NOx(t) are the amounts of carbon dioxide, nitrogen oxides, and sulfur dioxide emissions produced by units at the time t, respectively.
2.3. Uncertainties in Solar Radiation and Load
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
is the probability of each scenario,
is the value of variables in each scenario, and
is number of scenarios. Also, the coefficient of variance (
Cvx) is defined as follows.
where σ
x and
μx are the standard deviation and mean value of parameter x, respectively. If
Cvx 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 (
) and standard deviation (
) 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
and
are the mean and standard deviation of the load, respectively. Also,
Pd 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.
6. Simulation Results
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 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.
Table 7 and
Table 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.
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 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.