Next Article in Journal
Assessment of the Internal Catalyst Efficiency in a Diesel Engine of a Vehicle under the Conditions Simulating Real Driving
Previous Article in Journal
Densified Laser-Induced Graphene for Flexible Microsupercapacitors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploration of DG Placement Strategy of Microgrids via FMFO Algorithm: Considering Increasing Power Demand and Diverse DG Combinations

1
Academic Affairs Office, Jiujiang University, Jiujiang 332005, China
2
School of Computer, Information and Communication Engineering, Kunsan National University, Gunsan 54150, Korea
3
Department of Electrical and Electronics Engineering, P.A. College of Engineering and Technology, Pollachi 642003, Tamilnadu, India
*
Author to whom correspondence should be addressed.
Energies 2020, 13(24), 6568; https://doi.org/10.3390/en13246568
Submission received: 25 September 2020 / Revised: 24 November 2020 / Accepted: 10 December 2020 / Published: 13 December 2020
(This article belongs to the Section A1: Smart Grids and Microgrids)

Abstract

:
Distributed energy resource (DER) has been widely deployed, and distributed generation (DG) can complement the distribution system. Favorable DG deployment provides the grid-connected microgrid (MG) with stable voltage and reduces emission and power generation costs. DGs are considered distributed feeders, and MG is required to be operated under the optimal state. Reconfiguration is a practical approach to optimizing resource allocation. The optimal global solution is obtained via optimization algorithms. In this paper, three objectives are defined, namely, minimization of economic cost (ECC), emission cost (EMC), and voltage deviation (VD). Consequently, a fuzzy moth-flame optimization (FMFO) algorithm is proposed to coordinate the interests of multiple objectives. Moreover, the simulation is conducted based on the standard IEEE 33-bus radial distribution system (RDS), under which the impact of deployment of various DG type and quantity on the MG is explored. In particular, diverse DG combinations are tried under the increasing power demand, and a high-stable voltage strategy is proposed to meet the specific demands. The simulation results reveal that: (1) the DG type has a significant impact on ECC and EMC; (2) penetration level of DG shows a positive-like relationship with the MG stability; and (3) the proposed FMFO algorithm exhibits an efficient performance in convergence.

1. Introduction

Energy plays a vital role in social development. The massive consumption of traditional energy resources, represented by fossil fuels, has given rise to multiple problems, such as climate change [1] and environmental pollution [2]. With the rapid development of distributed energy resource (DER) technology, the demand for DER is also increasing dynamically, promoting profound changes in planning, configuration, operation, and control of the power system. Currently, DERs are accepted and deployed near houses as a part of the local power supply. A DER is a small-scale energy resource, making the traditional power generation mode undergo a significant transformation. In the traditional operation model, power is supplied by one or more centralized large-scale power plants, which provide power for long-distance loads through high-voltage transmission lines. The DER generation units, also called distributed generations (DGs), are integrated into the power network to provide power compensation for the power system [3]. The DER consists of controllable units (e.g., diesel engine (DE), microturbine (MT), and fuel cell (FC)) and uncontrollable units (e.g., wind turbine (WT) and photovoltaic (PV) cell) [4]. It is more flexible and faster to deploy the DGs compare with traditional power generators. Moreover, DGs’ deployment can bring higher reliability to the power system, provide higher power quality to the loads, reduce transmission loss in the distribution system, decrease power generation cost, and limit environmental pollution.
As the power system is deepening reform, traditional power systems are developing into smart grids (SGs). Microgrids (MGs) are considered cost-effective, environment-friendly, and user-satisfied solutions to the SG. The MG is a small-scale power generation and distribution system, which generally includes multiple components such as power converters (PCs), energy storage systems (ESSs), transformers, switchgears, transmission lines, and power consumers (loads) [5]. There are two operation modes of MG: grid-connected mode and islanded mode [6], respectively. However, considering unpredictable changes in the MG (such as transmission line failures, considerable fluctuations in loads, or DG faults) and the power distribution system’s economic and environmental factors, a practical solution is conducting MG reconfiguration [7]. The MG reconfiguration improves the system reliability, reduces the power generation cost, limits the emission pollution, reduces the power transmission loss (PTL), and stabilizes the voltage deviation (VD) in MGs, keeping the MG operated under the optimal state.
Various management strategies, control technologies, and optimization methods related to MG have emerged in recent years to provide reliable technical support for MGs. In [8], the progress of intelligent distribution technology is reviewed concisely, and the emerging technology and future development trend are determined to support the active management of the distribution network. Different MG frameworks are proposed in [9,10,11,12] to support energy management, power distribution, DG placement, and renewable energy sources. Hierarchical structures of MG control methods are introduced in [13,14,15,16,17] to handle the system control, energy storage, power distribution, and market behaviors.
Investigation and study on MG reconfiguration are not sufficient [18] and still facing numerous critical challenges. The MG reconfiguration is regarded as a multiobjective optimization process that requires the high performance of optimization algorithms. Considering DG deployment, [19,20,21] conducted the reconfiguration of the power systems to reduce PTL by applying novel optimization methods. Novel intelligent algorithms [22,23,24,25] are springing out to provide options to the optimization problems. The intelligent algorithms are high-efficient in solving single-objective optimization problems but mostly cannot independently deal with multiobjective optimization problems. [26,27,28,29,30] proposed hybrid methods, which combined conventional optimization algorithms with multiobjective approaches, such as stochastic method, Pareto dominance method, and fuzzy decision-making theory, to provide opportunities for multiobjective optimization solutions.
The construction of the power system infrastructure is a significant issue related to the livelihood of people. As the future of the power system, the MG presents an evident value in research and application. In the present studies, the research on multiobjective optimization algorithms is not sufficient; and few works have explored the impact of DG deployment on optimizing the MG system, especially on DG deployment strategies under increasing power demand.
This paper focuses on obtaining and analyzing the optimal solution for reconfigurable MGs. Three objectives: minimization of economic cost (ECC), minimization of emission cost (EMC), and minimization of voltage deviation (VD) are defined to describe multiple requirements of the MG, and a multiobjective optimization algorithm is in need for coordinating the global interests of multiple objectives. To supplement the shortcomings of the existing multiobjective optimization solutions, a novel FMFO algorithm, hybridized the nature-inspired Moth-flame optimizer (MFO) and heuristic fuzzy method, is proposed to solve the multiobjective problem defined in this paper. While considering the integration of MG and DG, the deployment strategy is crucial to achieving the MG voltage stability, especially in the high-stable VD required scenarios. Thus, the analysis and comparison of diverse DG combinations are necessary for developing the DG deployment strategy. Furthermore, the employment of three types of controllable DER generation units, DE, MT, and FC, are considered to explore the impact of diverse DG combinations on optimal reconfiguration results, such as the effect on load demand, bus voltage, and branch line current.
The main contributions of this paper are:
  • Uncontrollable DER generation units, DE, MT, and FC, are studied, and the corresponding models of ECC, EMC, and VD are established in mathematical expression patterns, respectively;
  • An FMFO algorithm is proposed to solve the multiobjective optimization problem and acquire the optimal global solution in MG reconfiguration;
  • An in-depth exploration of the potential impact of DG type and quantity on the power distribution system’s performance is conducted by employing of penetration level of DG;
  • The step-by-step changing power demands of prosumers are considered to study the optimization strategy under critical situations;
  • A method to realize the high-stable voltage was studied that can provide a DG selection strategy for those who require high-quality power.
The remaining parts of this paper are organized as follows: Section 2 introduces the formation of mathematical models, including all relevant objective functions and the corresponding constraints. Section 3 presents the primary research process and the proposed hybrid optimization algorithm. Section 4 describes several related cases and scenarios, and the results obtained in the cases are analyzed. Finally, Section 5 gives a conclusion to this paper.

2. Problem Formulation

2.1. Objective Functions

Three objectives, minimization of ECC, minimization of EMC, and minimization of VD, are defined based on the desired changes in the MG. The relevant objective functions can be optimized through the novel intelligent algorithms while conducting system reconfiguration, and the mathematical models for the three objectives are established.

2.1.1. Minimization of ECC

In the test system of this paper, there are two primary sources of the power supply: one is from the utility grid (substation), and the other is from DER generation units (DGs). In this case, the total power generation cost can only be produced by these two types of energy supply. It should be equal to the sum of DGs’ generation cost and the utility grid’s generation cost. Due to the commercial nature of power generation, the power generation cost is considered the ECC. Therefore, the total ECC can be defined as:
T ECC = G ECC + DE ECC + MT ECC + FC ECC
where T ECC is the total ECC, G ECC , DE ECC , MT ECC , and FC ECC represent the ECC caused by utility grid, DE, MT, and FC, respectively. Furthermore, the G ECC can be calculated as:
G ECC = P g r i d × K g r i d
where P g r i d is the active power supplied from the utility grid, and K g r i d is the coefficient for calculating the ECC of the utility grid. For the DGs, the calculation of the ECC depends on the type and quantity of DGs. This paper defines that the three types of DGs are DE, MT, and FC, respectively. The ECC of DE is expressed as:
DE ECC = i = 1 n DE a i + b i P DE i + c i P DE i 2
where P DE i is the active power generated by the i-th DE, n DE is the DE quantity, and a i , b i , c i are the coefficients of the ECC of DE. Similarly, the ECC of MT is expressed as:
MT ECC = j = 1 n MT c ng P MT j η MT j Q LHV
where P MT j is the active power generated by the j-th MT, n MT is MT quantity, c ng is the price of the natural gas, η MT j is the operation efficiency of the j-th MT, and Q LHV is the low-hot value of natural gas. The subsequent ECC of FC is expressed as:
FC ECC = k = 1 n FC c ng P FC k η FC k Q LHV
where P FC k is the power generated by the k-th FC, n FC is the FC quantity, and η FC k is the operation efficiency of the k-th FC.

2.1.2. Minimization of EMC

The power generation may lead to varying degrees of environmental pollution according to the variety of power resources. Among the substances emitted, CO 2 , SO 2 and NO x are considered to be the chief components that are the culprit for causing damage to the environment. For quantifying the degree of emissions, the EMC function is established as follows:
T EMC = G EMC + D E EMC + M T EMC + F C EMC
G EMC = ( CO 2 g r i d + NO x g r i d + SO 2 g r i d ) × P g r i d
DE EMC = i = 1 n DE ( CO 2 DE i + NO x DE i + SO 2 DE i ) × P DE i
MT EMC = j = 1 n MT ( CO 2 MT j + NO x MT j + SO 2 MT j ) × P MT j
FC EMC = k = 1 n FC ( CO 2 FC k + NO x FC k + SO 2 FC k ) × P FC k
where T EMC is the total EMC, G EMC , DE EMC , MT EMC , and FC EMC are the EMC of the utility grid, DE, MT, and FC, respectively. P DE i , P MT j , and P FC k represent the active power generated by the i-th DE, the j-th MT, and the k-th FC, respectively.

2.1.3. Minimization of VD

Voltage stability is a crucial index to reflect power quality. Controlling the bus VD within the tolerable range plays an essential role in improving the reliability of the distribution system [31]. VD refers to the difference between the real-time voltage and the reference voltage of all load buses, which can be expressed as:
V D = i = 1 N B |   V i V ref   |
where V i is the voltage of the i-th bus, V ref is the reference voltage, and N B is the number of load buses.

2.2. Constraints

When reconfiguring the MG, all possibilities should be considered to set the necessary constraints to meet the actual situation. According to the characteristics, this work divides the constraints into equality constraints and inequality constraints.

2.2.1. Equality Constraint

In this work, equation constraints are mainly discussing the power balance in the distribution system, and the total input power should be equal to the sum of all loads and transmission loss. The power balance formula is as follows:
{   P g r i d + P DE + P MT + P FC = P l o s s + P l o a d   Q g r i d + Q DE + Q MT + Q FC = Q l o s s + Q l o a d
{   P DE = i = 1 n DE P DE i ,   P MT = j = 1 n MT P M T j ,   P FC = k = 1 n FC P FC k     Q DE = i = 1 n DE Q DE i ,   Q MT = j = 1 n MT Q MT j ,   Q FC = k = 1 n F C Q FC k
where P g r i d , P DE , P MT , and P FC are the active power generated by the grid, DE, MT, and FC, respectively. Q g r i d , Q DE , Q MT , and Q FC stand for the reactive power generated by the grid, DE, MT, and FC, respectively. P l o s s and Q l o s s are the active and reactive power losses, P l o a d and Q l o a d are active and reactive load demands.

2.2.2. Inequality Constraint

Inequality constraints are a set of conditional constraints that consider all possibilities (such as the boundaries of variables) and can only be expressed with inequalities. Constraints on power generation operation:
max ( P i min , P i 0 D R i ) P i min ( P i max , P i 0 + U R i )
P i { P i m i n P i P i , 1 l P i , j 1 u P i P i , j l P i , n i u P i P i m a x
where P i is the current output power, and P i 0 is the previous output power, U R i is the up-ramp limit of the i-th generator (MW/time-period), and D R i is the down-ramp limit of the i-th generator (MW/time-period). Due to the physical characteristics of power equipment, the constraints of line flows are:
P i m i n P i P i m a x 1 i N G
Q i m i n Q i Q i m a x 1 i N G
V i m i n V i V i m a x 1 i N B
| S k | S k m a x 1 k N E
where N G is the number of generator buses, and N E is the number of line branches.
The constraints of DG Power factor:
0.8 P F 1
Each kind of dispatchable DG’s power handling capacity is limited because of its current and voltage control loops. Therefore, in any case, the power P gen , i , s and Q gen , i , s generated by DG follow the droop relationship:
P gen , i , s P gen ,   i , m a x , i N droop
Q gen , i , s Q gen ,   i , m a x , i N droop
The constraints of minimum in-service and out-service time:
{   T in , m t , s T in , m m i n       T out , m t , s T out , m m i n
The constraints of storage battery:
E SB t , s = {   E SB t 1 , s P SB t , s Δ t η ch P SB t , s < 0 E SB t 1 , s P SB t , s Δ t / η dis P SB t , s 0
P SB , ch m a x P SB t , s P SB , dis m a x
E SB m i n E SB t , s E SB m a x

3. Hybrid FMFO Algorithm

Intelligent optimization algorithms with reliable performance and fast convergence came into being and have impressed many researchers and are successfully applied in many fields. Compared with traditional algorithms, intelligent optimization algorithms have many advantages and higher efficiency in optimizing nonlinear or multi-extreme complex functions. For acquiring the optimal DG placement in the reconfigurable MGs, this paper proposes a hybrid multiobjective optimization algorithm which applied the fuzzy method to the MFO algorithm.

3.1. MFO Algorithm

The MFO algorithm [25] is a novel nature-inspired algorithm first proposed by Seyedali Mirjalili in 2015. The algorithm sets a fixed number of moths to wander in the variable space to search for the flame, the optimal solution of the objective function. The execution process of the MFO algorithm is summarized and described in Algorithm 1.
Algorithm 1 MFO Algorithm
Energies 13 06568 i001

3.2. Heuristic Fuzzy

The MG optimization involves three objectives: minimization of ECC, EMC, and VD. The scheme of heuristic fuzzy combined with the MFO algorithm is proposed to balance the benefit of multiple objectives. The fuzzy membership functions are employed for each objective function with upper and lower limits to enjoy high flexibility for making decisions. The developed membership functions can be modified to change the contribution of each objective function of the problem based on their requirements. The objective functions involved are described in detail below.

3.2.1. Membership Function for ECC

The total power generation cost of the solution should be minimized. Therefore, for a given load level E C C , the ratio of the total power generation cost before and after optimization of the distribution system is defined as follows:
E C C = E C C a f t e r E C C b e f o r e
In the equation, for lower E C C values, the solution can get more benefits. The membership function corresponding to E C C minimization is expressed as:
μ E C C = {   1 E C C E C C m i n E C C m a x E C C   E C C m a x E C C m i n E C C m i n < E C C < E C C m a x 0 E C C E C C m a x
In this work, E C C min and E C C max equal to 0.5 and 1, respectively. The value of E C C min refers to the minimum acceptable amount of saving to 50%.

3.2.2. Membership Function for EMC

The same as E C C , the EMC should be minimized, and the E M C is defined as follows:
E M C = E M C a f t e r E M C b e f o r e
The membership function corresponding to the minimization of E M C is expressed as,
μ E M C = {   1 E M C E M C m i n E M C m a x E M C E M C m a x E M C m i n E M C m i n < E M C < E M C m a x 0 E M C E M C m a x
where E M C min and E M C max are equal to 0.5 and 1, respectively. The value E M C min refers to the minimum acceptable amount of saving to 50%.

3.2.3. Membership Function for VD

Among the power flow constraints, bus voltage constraints are more sensitive when external equipment is connected to the system. Therefore, minimization of VD was considered as one of the objectives, along with the reduction of total PTL. The minimum bus voltage (MBV) of the system was obtained from the load flow equation for each new population of MFO. This membership function determines the deviation of the V D of the new population ( V D a f t e r ) to the initial V D ( V D b e f o r e ). The membership function for the population can be expressed as follows:
μ V D = {   1.0 f o r   V D x V D m a x V D m a x V D x V D m a x V D m i n f o r   V D m i n < V D x < V D m a x 0.0 f o r   V D x V D m a x
where V D m i n and V D m a x are considered as 0.05 and 0.1, respectively. The V D m i n reflects the maximum voltage of network buses is <1.05 (or the minimum voltage of network buses is >0.95), then μ V D = 1 . The value of V D x is calculated and expressed below:
V D x = V D a f t e r V D b e f o r e
From the Equations (28), (30), and (31), the fuzzified values for the E C C , E M C , and V D can be obtained as μ E C C , μ E M C , and μ V D , respectively. In order to obtain a compromise solution among these three fuzzified values, the fuzzy decision set is introduced, which is expressed as
μ D = m i n ( μ E C C , μ E M C , μ V D )
The equation is considered as the fitness function for the MFO algorithm. As a result, each population of the MFO produces a compromised solution among the three objective functions.

4. Simulation and Discussion

The simulations are prepared in a desktop computer installed with Microsoft Windows 10 Enterprise Edition as the operating system for coding and testing. The CPU of the computer is Intel Core (TM) i5-4660, equipped with an 8G Samsung DDR3-1600 as the memory. In terms of software, the Java NetBeans IDE 7.2 is adopted as a development tool to simulate the standard IEEE 33-bus RDS, which integrated with MGs deployed according to the loops. Figure 1 shows the schematic diagram of the test system. Since there are five loops in the standard IEEE 33-bus RDS, the test system maximizes the deployment of the hypothetical MGs in these loops to ensure that the MGs can generally work in the grid-connected mode and are also able to operate in the islanded mode under the specific conditions [32].
The primary purpose of the simulation is to obtain the optimal global solution for optimization and explore the optimal DG placement strategy. For evaluating the performance of the FMFO algorithm, the simulations are conducted for the case study, and numerous results are gathered for analysis.

4.1. Initialization and Settings

In the standard IEEE 33-bus RDS, the initial basic kV and MVA of the test system are set to 12.66 kV and 100 MVA, and the basic active power and reactive power demands are 3.715 MW and 2.3 MVAR, respectively. For facilitating the research, the simulations take the basic load as a reference. It assumes that the MBV limit is 0.9 per unit (pu) and that each line has a maximum branch current capacity of 400 A (excluding lines connected to the feeder). There are 33 buses in the test system (from 0 to 32). In the initial state, the normally open (NO) switches are S33, S34, S35, S36, and S37, and the normally closed (NC) switches are from S1 to S32. Table 1 shows the initial state of the lines and buses in the test system.
This investigation divides the primary parameters or variables involved in the optimization process into two categories: one is called generation-parameter, which plays a crucial role in energy generation, and the values can be collected from the corresponding objective models. Another category is named optimization-variable, which plays a decisive role in the simulation results, and the values are the optimal locations of DGs and the open switches. Table 2 and Table 3 list the primary generation-parameters and corresponding values in our simulation.
According to the characteristics, this study further divides the optimization-variable into the location-optimization variable, type-optimization variable, and quantity-optimization variable. The location-optimization variables refer to the locations of the switches and DGs, the type-optimization variables describe the types of DGs, and the quantity-optimization variables define the number of active DGs. The simulations also restricted each MG to one open switch to keep the radial structure. By changing the location of the open switches, the topology of the power distribution system is changed. Table 4 lists the boundaries of the location-optimization variable values.

4.2. Simulation and Case Study

Simulations are conducted in the test system, which has integrated the standard IEEE 33-bus RDS and five MGs. For exploring the optimal DG placement strategy, DE, MT, and FC are deployed with holistic consideration of the ECC, EMC, MBV, and PTL. It is assumed that all DGs operate at rated power, where the rated power of DE, MT, and FC are set to 0.4 MW/h, 0.6 MW/h, and 0.5 MW/h, respectively. Then, the changes of crucial parameters in reconfiguration were analyzed, and the obtained performance feedback is discussed from the perspective of multiple DG combinations. Through analyzing the results, the changes and their trends can be found to assist in formulating the DG placement strategy in different cases.

4.2.1. Case 1: Simple DG Deployment

In this phase, the number of DE, MT, and FC deployed in the test system is restricted to 1, but the types of DGs are not limited in simulation. The results of system reconfiguration shown in Table 5 are effectively obtained due to the low computational complexity. For better comparing the changes brought by reconfiguration after deploying DG(s), a simulation without DG deployment is conducted as a reference. The open switches were changed to the optimal locations after reconfiguration.
Since a DE deployed in the MGs, the EMC, MBV, and PTL are improved in evidence, while ECC is increased on the contrary due to the higher generation cost of DE. However, Table 5 discloses the apparent improvements in all aspects if any MG is deployed with an MT or FC. The ECC, EMC, and PTL were efficiently reduced, the MBV is significantly optimized, and the MBV trended to be more stable with the increasing usage of DGs. Different types of DG bring different changes to MG, which is closely related to their physical and electrical characteristics. Due to the ECC, DE has a negative contribution to reducing ECC. Therefore, the ECC of DG should be fully considered in the selection of MG before the deployment of DG. Compared with DE, FC and MT play a positive role in ECC and EMC, which can be considered while selecting DGs.

4.2.2. Case 2: Effect of DG Quantity

Most of the changes caused by DG deployment are positive, but the qualitative or quantitative measurement of the impact of DG type and quantity on the test system is still an issue. To make up for the deficiencies, Case 2 gradually increased the number of different types of DG from 1 to maximum and monitored the changes in system parameters. It is assumed that the DGs deployed in the test system are active and running at the rated power. Owing to the limitation of the total load demand, DGs cannot be deployed endlessly. According to the generation capacity of different DG types, the total DG generation should be less than or equal to the sum of load demands and line transmission losses. The maximum number of DE, MT, and FC are restricted to 9, 6, and 7, respectively. Moreover, the penetration level (PL) [36] of DG is employed to measure the impact of DG deployment on MG optimization. The PL of DGs P L D G is calculated as:
P L G D = P DE + P MT + P FC P l o a d + P l o s s × 100 %
In this phase, the simulation is carried out for probing the impact of increasing DG number on total ECC. Figure 2 shows the impact on ECC, EMC, MBV, and PTL in reconfiguration by increasing DE, MT, and FC quantity, limiting the maximum number of DG to 6 for comparison. Figure 2a indicates that the increase of DE has the opposite effect on reducing ECC, while with the increase in the number of MT and FC, the total ECC shows a linear-like slow downward trend. In the aspect of EMC shown in Figure 2b, with the increase of three types of DG, the EMC reduced in varying degrees, and the performance of FC is outstanding, followed by MT, and the improvement of DE is the lowest. In terms of voltage stability, it can be found in Figure 2c that there is a significant correlation between EMC and DG rated power. Finally, in terms of line transmission loss, pointed in Figure 2d, various types of DGs brought varying degrees of improvement, with FC still performs best among the three types.
The utilize of PL quantified the impact of DG with different generation capacities on the MG system targets and lessened the difference caused by distinct DG rated power. From the perspective of DG PL, as shown in Table 6, Table 7 and Table 8, the increasing number of DE, MT, and FC leads to an increase in DG penetration levels. In contrast, ECC, EMC, PTL, and MBV all show divergent trends and degrees on the change. Moreover, the increase of PL makes DG play a more significant role in MGs and profoundly impacts the optimal fitness values of essential objectives and the MG configurations.
The highest PL of DE, MT, and FC in tables are 94.64%, 94.56%, and 92.01%. Comparing the corresponding results of the maximum PL, the FC manifests the overall advantages over the other two types of DGs: the lower ECC, fewer EMC, less PTL, and more stable bus voltage.

4.2.3. Case 3: Diverse DG Combinations

This phase fixed the DG quantity to 5 in the simulation, based on which diverse combinations of DE, MT, and FC are considered to study the potential impact in system optimization on the essential parameters. Moreover, both manual mode and automatic mode are implemented while coding, where the manual mode is to simulate system reconfiguration in the case of specified DG type and quantity. In contrast, the automatic mode can automatically choose an optimal combination for system reconfiguration according to the given DG quantity. As listed in Table 9, the reconfiguration results of six given combinations and an auto-select combination show some different performance between two modes and all combinations. On account of the significant advantages of FC presented in Section 4.2.2, the automatic mode tends to consider the high performance of FC and the generation capacity of MT in choosing DGs, thus obtaining better configuration than the other six fixed combinations in terms of synthesis.

4.2.4. Case 4: Under Changing Power Demand

Power demand is usually not fixed but continually changing. MGs need to deal with sudden high-power demand, so it is necessary to consider the DG placement and system reconfiguration under various complicated conditions. Based on the forecasts of power demand and consumption in [37], this phase adjusted each bus’s energy demand within an acceptable range to imitate the changes in real power demands. For investigating the significance of system reconfiguration and DG placement under the changing power demand, the simulation increases the total power demand from 100% of the reference power demand to 200%, by 10% for each step.
Moreover, the simulation performed three scenarios based on the changing power demands (shown in Table 10). Scenario 1 simulates with a fixed DG combination, in which the number of DE, MT, and FC are set to 2 separately. Scenario 2 conducted an auto-select mode and set the DG quantity to 6, which decided the six FCs as the optimal combination. For obtaining more results, this simulation maximizes the deployment of DGs in Scenario 3 and lets the proposed strategy decide the DG quantity and the optimal DG combination. The comparison of the scenarios is shown in Figure 3.
The results obtained in the three scenarios are shown in Table 11, Table 12 and Table 13. It is significant that with the increase of power demand from reference quantity to 200%, the total ECC, total EMC, and PTL in the three scenarios increase to varying degrees; only the voltage stability shows a fluctuating decrease. Scenario 2 shows more significant advantages from a scenario perspective than Scenario 1 in terms of ECC and EMC. In terms of PTL and voltage stability, Scenario 2 performs slightly better than Scenario 1 in the reference power demand’s vicinity. As the power demand increases to 200%, the PTL and MBV of Scenario 2 are getting higher than those in Scenario 1. Scenario 3 performs significantly better than Scenario 1 and Scenario 2 in EMC, ECC, PTL, and voltage stability. Moreover, as power demand increases, this advantage becomes more apparent.

4.3. Analysis and Discussion

The significance of MG optimization is to achieve the optimal allocation of resources. DER plays a prominent role in regulating the peak energy demands, considering both economic and environmental protection. Deploying DGs in the reconfiguration of MG can solve many problems encountered in the development of SG.

4.3.1. Evaluation of Algorithm

This paper proposed an FMFO algorithm to solve the multiobjective issue. As shown in Figure 4, the improvement of system optimization by employing the FMFO algorithm is evident in the bus voltage and line current in MGs. The MBV was improved to 0.9378 (pu) after system reconfiguration without deploying DG. While deploying a DE, the MBV has reached 0.9503 (pu), even reached 0.9661 (pu) while deploying a simple DG combination. Simultaneously, the line current is optimized with the appropriate utilization of DG(s) or DG combination; each line’s current is significantly decreased.
In order to reflect the convergence of the proposed algorithm, our simulation randomly selected ten results in the case of without DG deployment, with DE deployment, with MT deployment, and with FC deployment, respectively. Figure 5 shows the convergence of the proposed algorithm on ECC that the optimal solutions are obtained within a maximum of 20 iterations without DG deployment and within 44 iterations among the 30 runs while DG is deployed. Moreover, a fuzzy genetic algorithm (FGA) and fuzzy particle swarm optimization (FPSO) are implemented in the simulation to compare the convergence with the proposed FMFO on the PTL. As shown in Figure 6, 20 runs of FGA, FPSO, and FMFO are conducted separately to calculate the average PTL. The proposed FMFO algorithm shows advantages in both simulations deployed without DG and deployed with an MT. In the simulations without DG deployed, the FMFO, FPSO, and FGA obtained the optimal solutions at the average iteration numbers of 9.7, 14.9, and 34.2, respectively; and in the simulations with an MT deployed, the FMFO, FPSO, and FGA obtained the optimal solutions at the average iteration numbers of 13.6, 18.5, and 46.3, respectively.
The effect of the DG penetration level on improving PTL is usually positive [36]. As shown in Table 14, this simulation obtained better results in which the PTL is 62.18 kW under the auto-select mode, in the case where the DG penetration level was 47.65%. Compared with the optimal results in [38], the PTL obtained in this paper was reduced by 7.35%; meanwhile, the MBV was still kept at a considerable level.

4.3.2. Data Analysis and DG Placement Strategy

The deployment of DG(s) leads to a significant impact on MG optimization, as shown in Table 5. Compared with no DG deployed, employing a DE has increased the ECC from $169.6002/h to $297.4400/h by 75.38%. In contrast, the deployment of an MT and FC have decreased the corresponding ECC by 2.98% and 3.80%. In terms of EMC, the applying of a DE, MT, and FC has reduced the total EMC from 7982.8072 lb/h to 7608.5086 lb/h, 7536.1043 lb/h, and 7332.4768 lb/h by 3.60%, 4.52%, and 7.10%, respectively. Similarly, by applying a DE, MT, and FC, the PTL decreased by 22.15%, 30.09%, and 26.43%, and the MBV increased from 0.9378 pu to 0.9503 pu, 0.9522 pu, and 0.9450 pu, respectively. When DE, MT, and FC are combined, the effects of different DGs on system optimization are superimposed.
Through analyzing the impact of DG type on system optimization and the superposition effect of DG quantity, the simulation results where DG quantity increases from 1 to the maximum are shown in Table 6, Table 7 and Table 8. In corresponding Figure 2, Figure 3 and Figure 4, the ECC curves and EMC curves present a quasi-linear decrease with the increase of DG quantity except for DE; meanwhile, the PTL shows an exponent-like reduction, and the MBV shows an irregular improvement. The changing trends of ECC and EMC rely mainly on the characteristics of DGs.
As further exploring the reconfiguration under diverse DG combinations, Table 9 shows that the EMC, ECC, PTL, and the MBV improved to varying degrees under various DG combinations. Compared with Table 5, the MBV is more stable, illustrating the importance of DG quantity increment to voltage stability. However, it is worth noticing that the MBV value of auto-select mode is high, reaching 0.9739 pu, due to its higher DG penetration level among all DG combinations. From the perspective of MBV stability, the performance of an appropriate multitype DG combination is slightly better than that of the optimal single-type DG combination in case that the DG quantity is fixed.
Thus, if the high voltage stability is required, it is needed to maximize the utilization of DGs, and then choose the appropriate multitype DG combination to meet the different power replenishment needs of buses.
Compare Table 11 with Table 12, when the power demand increases from 100% to 200%, the PTL and MBV of the combination with two DEs, two MTs, and two FCs are 252.1833 kW and 0.9325 pu; the combination of auto-select mode are 243.3098 kW and 0.9396 pu. The optimal DG combination showed more considerable advantages in PTL and MBV with increased power demand. Nevertheless, the ECC and EMC of auto-select mode are $447.5775/h and 13,709.5232 lb/h that is much better than that of the fixed 6 DGs $580.3721/h and 13,777.4195 lb/h. With the increase of power demand in the distribution system, the fixed DG quantity makes the penetration level decreased. Table 13 shows relatively stable DG penetration levels by maximizing DGs’ deployment under the auto-select mode. Compared with Table 11 and Table 12, the ECC, EMC, PTL, and MBV in Table 13 are comprehensively improved due to the higher DG penetration level. The ECC, EMC, PTL, and MBV have also been improved to $261.6486/h, 8267.5435 lb/h, 159.0155 kW, and 0.9607 pu under the 200% power demand.
As a result of the bus voltage’s reliable stability in the simulation, the FMFO algorithm paid more attention to energy saving and emission reduction while ensuring high-quality power supplying. Maximizing the DG penetration level can effectively optimize the power distribution system’s output, but the appropriate DG combination is necessary to improve system parameters further. The fundamental cause is that making the current average and at a lower level will stable the bus voltage and minimize the energy transmission loss.

5. Conclusions

This paper conducted the system optimization of MGs based on the standard IEEE 33-bus RDS and proposed an FMFO algorithm to solve the multiobjective optimization problem. The proposed FMFO algorithm showed an efficient performance on convergence comparing with classical FGA and FPSO. The obtained optimal solution via MG reconfiguration has considered integrating with DGs and achieved minimum PTL, ECC, EMC, and VD comparing with existing research. Moreover, our simulations have considered various DG combinations to explore the effect of DG type and quantity on MG reconfiguration results. Furthermore, we have investigated the results in four step-by-step cases and gained the optimal solution for each. The optimal DG placement and combination were explored under the increasing power demand, which formed a DG deployment strategy for decision-making. The simulation results reveal that: the proposed FMFO algorithm is a reliable supplement for multiobjective optimization approaches in MG reconfiguration; it requires fewer iterations to obtain the optimal solution than other optimizers; appropriate DG selection and deployment strategy is conducive to obtaining the optimal global solution, especially under the critical power demand conditions.

Author Contributions

Conceptualization, X.W. and I.-h.R.; methodology, T.S. and H.K.; investigation, X.W. and T.S.; writing—original draft preparation, X.W.; writing—review and editing, T.S. and H.K.; funding acquisition: I.-h.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This work was supported by the Institute for Information and Communications Technology Promotion (IITP) grant funded by the Korea government (MSIT) (No. 2018-0-00508), Development of blockchain-based embedded devices and platform for MG security and operational efficiency. This work was supported in part by the Human Resources Development Program (Grant No. 20194010201800) of the Korea Institute of Energy Technology Evaluation and Planning (KETEP) grants funded by the Korean government (Ministry of Trade, Industry, and Energy).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
DEDiesel engine
DERDistributed energy resource
DGDistributed generation
ECCEconomic cost
EMCEmission cost
ESSEnergy storage systems
FCFuel cell
FMFOFuzzy moth-flame optimization
MGmicrogrid
MTmicroturbine
MBVMinimum bus voltage
MFOMoth-flame optimizer
NCnormally closed
NOnormally open
PLPenetration level
PVphotovoltaic
PCPower converter
PTLPower transmission loss
RDSRadial distribution system
SGSmart grid
VDVoltage deviation
WTWind turbine

References

  1. Höök, M.; Tang, X. Depletion of fossil fuels and anthropogenic climate change—A Review. Energy Policy 2013, 52, 797–809. [Google Scholar] [CrossRef] [Green Version]
  2. Withagen, C. Pollution and exhaustibility of fossil fuels. Resour. Energy Econ. 1994, 16, 235–242. [Google Scholar] [CrossRef]
  3. Singh, M.; Khadkikar, V.; Chandra, A.; Varma, R.K. Grid interconnection of renewable energy sources at the distribution level with power-quality improvement features. IEEE Trans. Power Deliv. 2011, 26, 307–315. [Google Scholar] [CrossRef]
  4. Li, Y.; Wang, P.; Gooi, H.B.; Ye, J.; Wu, L. Multi-Objective optimal dispatch of microgrid under uncertainties via interval optimization. IEEE Trans. Smart Grid 2019, 10, 2046–2058. [Google Scholar] [CrossRef]
  5. Meng, L.; Sanseverino, E.R.; Luna, A.; Dragicevic, T.; Vasquez, J.C.; Guerrero, J.M. Microgrid supervisory controllers and energy management systems: A Literature Review. Renew. Sustain. Energy Rev. 2016, 60, 1263–1273. [Google Scholar] [CrossRef]
  6. Parhizi, S.; Lotfi, H.; Khodaei, A.; Bahramirad, S. State of the art in research on microgrids: A Review. IEEE Access 2015, 3, 890–925. [Google Scholar] [CrossRef]
  7. Kavousi-Fard, A.; Khodaei, A.; Bahramirad, S. Improved efficiency, enhanced reliability and reduced cost: The transition from static microgrids to reconfigurable microgrids. Electr. J. 2016, 29, 22–27. [Google Scholar] [CrossRef]
  8. Zhao, J.; Wang, C.; Zhao, B.; Lin, F.; Zhou, Q.; Wang, Y. A review of active management for distribution networks: Current status and future development trends. Electr. Power Compon. Syst. 2014, 42, 280–293. [Google Scholar] [CrossRef]
  9. Niknam, T.; Kavousi-Fard, A. Optimal energy management of smart renewable micro-grids in the reconfigurable systems using adaptive harmony search algorithm. Int. J. Bio-Inspired Comput. 2016, 6, 184–194. [Google Scholar] [CrossRef]
  10. Gazijahani, F.S.; Salehi, J. Stochastic multi-objective framework for optimal dynamic planning of interconnected microgrids. IET Renew. Power Gener. 2017, 11, 1749–1759. [Google Scholar] [CrossRef]
  11. Wang, Z.; Chen, B.; Wang, J.; Kim, J.; Begovic, M.M. Robust optimization based optimal dg placement in microgrids. IEEE Trans. Smart Grid 2014, 5, 2173–2182. [Google Scholar] [CrossRef]
  12. Khorram-Nia, R.; Bahmani-Firouzi, B.; Simab, M. Optimal switching in reconfigurable microgrids considering electric vehicles and renewable energy sources. J. Renew. Sustain. Energy 2018, 10, 1–18. [Google Scholar] [CrossRef]
  13. Olivares, D.E.; Mehrizi-Sani, A.; Etemadi, A.H.; Cañizares, C.A.; Iravani, R.; Kazerani, M.; Hajimiragha, A.H.; Gomis-Bellmunt, O.; Saeedifard, M.; Palma-Behnke, R.; et al. Trends in microgrid control. IEEE Trans. Smart Grid 2014, 5, 1905–1919. [Google Scholar] [CrossRef]
  14. Wu, X.; Wang, X.; Qu, C. A Hierarchical framework for generation scheduling of microgrids. IEEE Trans. Power Deliv. 2014, 29, 2448–2457. [Google Scholar] [CrossRef]
  15. Palizban, O.; Kauhaniemi, K.; Guerrero, J.M. Microgrids in active network management—Part I: Hierarchical control, energy storage, virtual power plants, and market participation. Renew. Sustain. Energy Rev. 2014, 36, 428–439. [Google Scholar] [CrossRef] [Green Version]
  16. Bidram, A.; Davoudi, A. Hierarchical structure of microgrids control system. IEEE Trans. Smart Grid 2012, 3, 1963–1976. [Google Scholar] [CrossRef]
  17. Guerrero, J.M.; Vasquez, J.C.; Matas, J.; de Vicuna, L.G.; Castilla, M. Hierarchical control of droop-controlled AC and DC microgrids—A general approach toward standardization. IEEE Trans. Ind. Electron. 2011, 58, 158–172. [Google Scholar] [CrossRef]
  18. Shariatzadeh, F.; Vellaithurai, C.B.; Biswas, S.S.; Zamora, R.; Srivastava, A.K. Real-Time implementation of intelligent reconfiguration algorithm for microgrid. IEEE Trans. Sustain. Energy 2014, 5, 598–607. [Google Scholar] [CrossRef]
  19. Abd El-salam, M.F.; Beshr, E.; Eteiba, M.B. A new hybrid technique for minimizing power losses in a distribution system by optimal sizing and siting of distributed generators with network reconfiguration. Energies 2018, 11, 3351. [Google Scholar] [CrossRef] [Green Version]
  20. Badran, O.; Mekhilef, S.; Mokhlis, H.; Dahalan, W. Optimal reconfiguration of distribution system connected with distributed generations: A Review of different methodologies. Renew. Sustain. Energy Rev. 2017, 73, 854–867. [Google Scholar] [CrossRef]
  21. Jabbari-Sabet, R.; Moghaddas-Tafreshi, S.-M.; Mirhoseini, S.-S. Microgrid operation and management using probabilistic reconfiguration and unit commitment. Int. J. Electr. Power Energy Syst. 2016, 75, 328–336. [Google Scholar] [CrossRef]
  22. Eskandar, H.; Sadollah, A.; Bahreininejad, A.; Hamdi, M. Water cycle algorithm–A novel metaheuristic optimization method for solving constrained engineering optimization problems. Comput. Struct. 2012, 110–111, 151–166. [Google Scholar] [CrossRef]
  23. Mirjalili, S.; Lewis, A. The whale optimization algorithm. Adv. Eng. Softw. 2016, 95, 51–67. [Google Scholar] [CrossRef]
  24. Mirjalili, S. The ant lion optimizer. Adv. Eng. Softw. 2015, 83, 80–98. [Google Scholar] [CrossRef]
  25. Mirjalili, S. Moth-flame optimization algorithm: A novel nature-inspired heuristic paradigm. Knowl. Based Syst. 2015, 89, 228–249. [Google Scholar] [CrossRef]
  26. Sedighizadeh, M.; Bakhtiary, R. Optimal multi-objective reconfiguration and capacitor placement of distribution systems with the Hybrid Big Bang–Big Crunch algorithm in the fuzzy framework. Ain Shams Eng. J. 2016, 7, 113–129. [Google Scholar] [CrossRef] [Green Version]
  27. Jafari, A.; Ganjeh Ganjehlou, H.; Baghal Darbandi, F.; Mohammadi-Ivatloo, B.; Abapour, M. Dynamic and multi-objective reconfiguration of distribution network using a novel hybrid algorithm with parallel processing capability. Appl. Soft Comput. 2020, 90, 106146. [Google Scholar] [CrossRef]
  28. Chen, Q.; Wang, W.; Wang, H.; Wu, J.; Li, X.; Lan, J. A social beetle swarm algorithm based on grey target decision-making for a multi-objective distribution network reconfiguration considering partition of time intervals. IEEE Access 2020, 8, 204987–205013. [Google Scholar] [CrossRef]
  29. Souifi, H.; Kahouli, O.; Hadj Abdallah, H. Multi-objective distribution network reconfiguration optimization problem. Electr. Eng. 2019, 101, 45–55. [Google Scholar] [CrossRef]
  30. Kianmehr, E.; Nikkhah, S.; Rabiee, A. Multi-objective stochastic model for joint optimal allocation of DG units and network reconfiguration from DG owner’s and DisCo’s perspectives. Renew. Energy 2019, 132, 471–485. [Google Scholar] [CrossRef]
  31. Li, H.; Li, F.; Xu, Y.; Rizy, D.T.; Kueck, J.D. Adaptive voltage control with distributed energy resources: Algorithm, theoretical analysis, simulation, and field test verification. IEEE Trans. Power Syst. 2010, 25, 1638–1647. [Google Scholar] [CrossRef]
  32. Ganjian-Aboukheili, M.; Shahabi, M.; Shafiee, Q.; Guerrero, J.M. Seamless transition of microgrids operation from grid-connected to islanded mode. IEEE Trans. Smart Grid 2020, 11, 2106–2114. [Google Scholar] [CrossRef]
  33. Li, P.; Xu, D.; Zhou, Z.; Lee, W.-J.; Zhao, B. Stochastic optimal operation of microgrid based on chaotic binary particle swarm optimization. IEEE Trans. Smart Grid 2016, 7, 66–73. [Google Scholar] [CrossRef]
  34. El-Ela, A.A.A.; El-Sehiemy, R.A.; Abbas, A.S. Optimal placement and sizing of distributed generation and capacitor banks in distribution systems using water cycle algorithm. IEEE Syst. J. 2018, 12, 3629–3636. [Google Scholar] [CrossRef]
  35. Kefayat, M.; Lashkar Ara, A.; Nabavi Niaki, S.A. A hybrid of ant colony optimization and artificial bee colony algorithm for probabilistic optimal placement and sizing of distributed energy resources. Energy Convers. Manag. 2015, 92, 149–161. [Google Scholar] [CrossRef]
  36. Srivastava, A.K.; Kumar, A.A.; Schulz, N.N. Impact of distributed generations with energy storage devices on the electric grid. IEEE Syst. J. 2012, 6, 110–117. [Google Scholar] [CrossRef]
  37. Eseye, A.T.; Lehtonen, M.; Tukia, T.; Uimonen, S.; John Millar, R. Machine learning based integrated feature selection approach for improved electricity demand forecasting in decentralized energy systems. IEEE Access 2019, 7, 91463–91475. [Google Scholar] [CrossRef]
  38. Mohamed Imran, A.; Kowsalya, M.; Kothari, D.P. A novel integration technique for optimal network reconfiguration and distributed generation placement in power distribution networks. Int. J. Electr. Power Energy Syst. 2014, 63, 461–472. [Google Scholar] [CrossRef]
  39. Saonerkar, A.K.; Bagde, B.Y. Optimized DG placement in radial distribution system with reconfiguration and capacitor placement using genetic algorithm. In Proceedings of the 2014 IEEE International Conference on Advanced Communications, Control and Computing Technologies, Ramanathapuram, India, 8–10 May 2014; pp. 1077–1083. [Google Scholar]
  40. Rao, R.S.; Ravindra, K.; Satish, K.; Narasimham, S.V.L. Power loss minimization in distribution system using network reconfiguration in the presence of distributed generation. IEEE Trans. Power Syst. 2013, 28, 317–325. [Google Scholar] [CrossRef]
  41. Rajaram, R.; Sathish Kumar, K.; Rajasekar, N. Power system reconfiguration in a radial distribution network for reducing losses and to improve voltage profile using modified plant growth simulation algorithm with Distributed Generation (DG). Energy Rep. 2015, 1, 116–122. [Google Scholar] [CrossRef] [Green Version]
Figure 1. A single-line diagram of standard IEEE 33-bus radial distribution system (RDS).
Figure 1. A single-line diagram of standard IEEE 33-bus radial distribution system (RDS).
Energies 13 06568 g001
Figure 2. ECC, emission cost (EMC), minimum bus voltage (MBV), and power transmission loss (PTL) changing curve with the increasing number of DGs.
Figure 2. ECC, emission cost (EMC), minimum bus voltage (MBV), and power transmission loss (PTL) changing curve with the increasing number of DGs.
Energies 13 06568 g002
Figure 3. Comparison of ECC, EMC, PTL, and MBV under changing power demand in three scenarios.
Figure 3. Comparison of ECC, EMC, PTL, and MBV under changing power demand in three scenarios.
Energies 13 06568 g003
Figure 4. Effect of DG deployment on bus voltage and line current of system reconfiguration.
Figure 4. Effect of DG deployment on bus voltage and line current of system reconfiguration.
Energies 13 06568 g004
Figure 5. The convergence of the proposed algorithm on ECC.
Figure 5. The convergence of the proposed algorithm on ECC.
Energies 13 06568 g005
Figure 6. Comparison of convergence of fuzzy genetic algorithm (FGA), fuzzy particle swarm optimization (FPSO), and fuzzy moth-flame optimization (FMFO) on power transmission loss (PTL).
Figure 6. Comparison of convergence of fuzzy genetic algorithm (FGA), fuzzy particle swarm optimization (FPSO), and fuzzy moth-flame optimization (FMFO) on power transmission loss (PTL).
Energies 13 06568 g006
Table 1. The initial state of the distribution test system.
Table 1. The initial state of the distribution test system.
MGNC SwitchesNO SwitchBuses
1S4, S5, S6, S7, S19, S20S33{3,4,5,6,7,18,19,20}
2S8, S9, S10, S11, S21S35{8,9,10,11,21}
3S12, S13, S14S34{12,13,14}
4S23, S24, S25, S26, S27, S28S37{22,23,24,25,26,27,28}
5S15, S16, S17, S29, S30, S31, S32S36{15,16,17,29,30,31,32}
Table 2. Primary parameters of economic cost (ECC) models [33].
Table 2. Primary parameters of economic cost (ECC) models [33].
ParameterValueParameterValue
K g r i d 0.044 $/kWh c ng 0.3
a i 6.0 η MT j 0.8
b i 0.012 Q LHV 9.7 kWh/m3
c i 0.00085 η FC k 0.9
Table 3. Emission coefficients related to grid and distributed energy resources (DERs) [34,35].
Table 3. Emission coefficients related to grid and distributed energy resources (DERs) [34,35].
Emission TypeEmission Factors (lb/MWh)
GridDEMTFC
CO22031149415961078
SO211.60.0080.0080.006
NOx5.061.150.440.03
Table 4. Boundaries of location-optimization variable values.
Table 4. Boundaries of location-optimization variable values.
VariableLower BoundaryUpper BoundaryVariableLower BoundaryUpper Boundary
LOS1 117LOS4 117
LOS2 116LOS5 118
LOS3 114LDGi 2132
1 LOS: location of the open switch; 2 LDGi: location of the i-th distributed generation (DG).
Table 5. Reconfiguration results for single DG and corresponding combinations.
Table 5. Reconfiguration results for single DG and corresponding combinations.
DG CombinationOpen SwitchesDG LocationECC ($/h)EMC (lb/h)Ploss(kW)MBV (pu)
--7,9,14,32,37--169.60027892.8072139.54970.9378 (31)
DE7,9,14,31,3731297.44007608.5086108.63740.9503 (30)
MT7,9,14,30,3731164.54867536.104397.56240.9522 (30)
FC7,9,14,30,3731163.15967332.4768102.66930.9450 (30)
DE/MT7,10,14,28,3129/17293.08487284.208782.47460.9619 (30)
DE/FC7,9,14,28,3129/32291.61117076.645685.65950.9603 (31)
MT/FC7,10,14,28,3117/29158.90707012.958078.84140.9637 (31)
DE/MT/FC7,9,14,28,3011/31/28287.59226767.996767.13990.9584 (30)
Auto-select (3 FCs)7,11,14,28,307,28,31156.58366922.201462.18130.9661 (30)
Table 6. Reconfiguration results of the increasing diesel engine (DE) deployment.
Table 6. Reconfiguration results of the increasing diesel engine (DE) deployment.
Quantity of DEPL (%)ECC ($/h)EMC (lb/h)Ploss (kW)MBV (pu)
--0.0000%169.60027892.8072139.54970.9378 (31)
110.1722%297.44007608.5086108.63740.9503 (30)
220.5395%425.81857349.272789.96490.9542 (32)
331.0433%554.37267098.216875.28710.9577 (30)
441.6333%683.07776854.181764.03800.9653 (32)
552.2608%811.92326616.683355.98120.9633 (30)
662.8975%940.87616384.183850.36560.9640 (31)
773.5551%1069.87676153.905545.83480.9709 (31)
884.1063%1199.03365930.896944.85430.9699 (31)
994.6364%1328.21875709.203544.51610.9677 (30)
Table 7. Reconfiguration results of the increasing microturbine (MT) deployment.
Table 7. Reconfiguration results of the increasing microturbine (MT) deployment.
Quantity of MTPL (%)ECC ($/h)EMC (lb/h)Ploss (kW)MBV (pu)
--0.0000%169.60027892.8072139.54970.9378 (31)
115.3448%164.54867536.104397.56240.9522 (30)
231.0358%160.38487220.715075.75110.9637 (31)
346.8828%156.58366922.201462.18130.9661 (30)
462.8166%152.96786632.314052.82420.9667 (31)
578.7368%149.53306350.853447.58240.9745 (31)
694.5581%146.26336077.074146.09180.9781 (30)
Table 8. Reconfiguration results of the increasing MT deployment.
Table 8. Reconfiguration results of the increasing MT deployment.
Quantity
of FC
PL (%)ECC ($/h)EMC (lb/h)Ploss (kW)MBV (pu)
--0.0000%169.60027892.8072139.54970.9378 (31)
112.7540%163.15967332.4768102.66930.9450 (30)
225.7797%157.43286805.367982.01310.9604 (31)
338.9600%151.97866290.943067.55130.9596 (30)
452.2406%146.68425783.950956.71940.9679 (31)
565.5179%141.58725286.151350.37670.9687 (30)
678.7867%136.59334793.145246.37500.9670 (30)
792.0122%131.68944304.331844.42090.9819 (31)
Table 9. Reconfiguration results under various combinations of 5 DGs.
Table 9. Reconfiguration results under various combinations of 5 DGs.
DG Combination (DE/MT/FC)Open SwitchesDG Location (DE/MT/FC)PG 1 (MW)ECC ($/h)EMC (lb/h)Ploss (kW)MBV 2 (pu)BN 3PL (%)
DGGridDGGridTotalDGGridTotal
1/1/36,9,14,31,3710/24/7,29,322.51.265221.54255.661277.2033172.9862590.3255763.31150.0170.97223165.5302%
1/2/211,28,30,33,3413/7,24/23,312.61.164227.55651.217278.7733591.8372383.5185975.35549.0210.96433068.1870%
1/3/16,9,14,30,3711/7,24,32/282.71.064233.57046.796280.3654010.6882177.7576188.44448.5340.96583070.8277%
2/1/210,28,31,33,3413,32/7/24,292.41.366351.16060.093411.2543232.0312796.6106028.64150.7590.97023162.8845%
2/2/111,28,30,33,3414,28/7,24/312.51.265357.17455.653412.8263650.8822589.9406240.82249.8290.96583065.5367%
3/1/110,28,31,33,3414,29,32/24/72.31.467480.77864.547545.3253291.0763003.8806294.95751.9820.97173160.2258%
auto-select9,14,27,30,33--/7,23,28,31/112.90.862109.96637.947147.9124370.4931765.9506136.44447.4240.97393076.1185%
1 PG: power generation; 2 MBV: minimal bus voltage; 3 BN: bus number.
Table 10. The changing power demands of all buses from 100% to 200%.
Table 10. The changing power demands of all buses from 100% to 200%.
Bus Number100% (MW)110% (MW)120% (MW)130% (MW)140% (MW)150% (MW)160% (MW)170% (MW)180% (MW)190% (MW)200% (MW)
PQPQPQPQPQPQPQPQPQPQPQ
10.10000.06000.14080.09290.06400.04220.17090.10250.16100.10140.25330.15200.18240.11160.18340.11000.15000.09450.17440.11100.30860.1944
20.09000.04000.10560.05160.15200.07430.13370.06540.13300.06210.12000.05600.15680.07110.16690.07420.20400.09520.20750.09780.17140.0800
30.12000.08000.09680.07100.12800.09390.11890.08720.14700.10290.14000.09800.17600.11970.20190.13460.19800.13860.18610.13150.18290.1280
40.06000.03000.11440.06290.16800.09240.23030.11510.14470.07600.20440.10730.21970.11210.24520.12260.27000.14180.21240.11260.24950.1310
50.06000.02000.13200.04840.08800.03230.11140.04090.11200.03920.13330.04670.08960.03050.05560.01850.10200.03570.18710.06610.20000.0700
60.20000.10000.15840.08710.18400.10120.10400.05720.21470.11280.14440.07580.20050.10230.21020.10500.23400.12290.29820.15800.22670.1190
70.20000.10000.18480.09240.13600.07480.12630.06950.22870.12010.21780.11430.20690.10550.20400.10200.19800.10400.28650.15180.30670.1610
80.06000.02000.09680.03550.20000.07330.08910.03270.09570.03350.00440.00160.08850.03010.07420.02480.19200.06720.23680.08370.06670.0233
90.06000.02000.07920.02900.11400.04180.23030.08440.13530.04740.24060.08420.24110.08190.29470.09820.25350.08870.13790.04870.25760.0902
100.04500.03000.10560.07740.15200.11150.18570.12380.11780.08250.15220.10660.16590.11280.17410.11610.21900.15330.19630.13870.19900.1393
110.06000.03500.07920.05080.09600.06160.15600.10010.11200.06860.16670.10210.16640.09900.19160.11180.18000.11030.13450.08310.19430.1190
120.06000.03500.09680.06210.10400.06670.14110.09060.11200.06860.14670.08980.14080.08380.14630.08530.16200.09920.15980.09880.18860.1155
130.12000.08000.14080.10330.16800.11720.20060.13370.18430.12900.21560.15090.23150.15740.25350.16900.26400.18480.24850.17560.27620.1933
140.06000.01000.10560.01940.23200.04250.16340.03000.12130.02120.06890.01210.15790.02680.16480.02750.27000.04730.25920.04580.12760.0223
150.06000.02000.11440.04190.17600.06450.15690.05750.12160.04460.11000.03850.14600.04960.14350.04780.21680.07590.22790.08050.16860.0590
160.06000.02000.07920.02900.10400.03810.16340.05990.11430.04190.16890.05910.17390.05920.20190.06730.19200.06720.13930.04920.19620.0687
170.09000.04000.10560.05160.15200.07440.10400.05080.12370.06050.08440.03940.12690.05750.12570.05590.18000.08400.21140.09960.14100.0658
180.09000.04000.18480.09030.12000.05870.15600.07630.16100.07510.19330.09020.13120.05950.08650.03850.14400.06720.26020.12260.28570.1333
190.09000.04000.14960.07310.07200.03520.16340.07990.15400.07180.23330.10890.16000.07250.14630.06500.13800.06440.18900.08910.29710.1387
200.09000.04000.13200.05870.08800.04300.09660.04720.12830.06090.13560.06330.10350.04690.07830.03480.10800.05040.19190.09040.20190.0942
210.09000.04000.11440.05590.16800.08210.15600.07630.14230.06640.13560.06330.17390.07880.18550.08240.22800.10640.22510.10600.19050.0889
220.09000.05000.12320.07530.10400.06360.13370.08170.13770.08030.16440.09590.14610.08280.14220.07900.15600.09100.19000.11190.22100.1289
230.42000.20000.14960.07840.15200.07960.15600.08170.38270.18220.37780.17990.46930.22800.61200.29140.39000.19500.28060.14160.42100.2105
240.42000.20000.18480.09680.11200.05870.14110.07390.38730.18440.40220.20110.43310.21030.54190.25810.32400.16200.28840.14560.46480.2324
250.06000.02500.07920.03630.20000.09160.18570.08510.12130.05310.11560.05060.19630.08340.23280.09700.28200.12340.20660.09120.15050.0658
260.06000.02500.11440.05240.08800.04030.12630.05790.11200.04900.14670.06420.11520.04900.10100.04210.12600.05510.16760.07400.20000.0875
270.06000.02000.13200.04840.12000.04400.11890.04360.11430.04000.11560.04040.09710.03300.06590.02200.13200.04620.20950.07400.18480.0647
280.12000.07000.14960.09600.19200.12320.16340.10240.17500.10210.15330.09230.18880.11230.18960.11060.24600.15070.27960.17290.22860.1400
290.20000.60000.17050.51150.21600.61560.21540.64630.25290.75860.25420.76250.30520.91990.34771.03570.33980.99890.31901.01040.32860.9973
300.15000.07000.16720.08580.18400.08590.11890.06100.18670.08710.13110.06120.16210.07720.14630.06830.21000.10290.30010.14850.22100.1083
310.21000.10000.18480.08800.13600.07120.17090.08140.24970.11890.27780.13230.26130.12690.28020.13340.24000.11950.28160.14210.35810.1790
320.06000.04000.11440.07650.08800.06450.14110.09410.11670.07780.16440.10960.13010.08850.12160.08110.13800.09660.16560.11710.21520.1507
Total3.71502.30004.08652.53004.45802.76004.82952.99005.20103.22005.57253.45005.94403.68006.31553.91006.68704.14007.05854.37007.43004.6000
Table 11. Scenario 1: system reconfiguration results of DG combinations under 6 DGs and increasing power demand.
Table 11. Scenario 1: system reconfiguration results of DG combinations under 6 DGs and increasing power demand.
PD (%)Open SwitchesDG Locations (DE/MT/FC)PG 1 (MW)ECC ($/h)EMC (lb/h)Ploss (kW)MBV 2 (pu)BN 3PL (%)
DGGridDGGridTotalDGGridTotal
1007,9,14,28,307,11/23,32/6,283.00.7602374.356033.4484407.80444189.91556.61095746.510945.19010.96863079.7832%
1107,9,14,28,327,31/11,26/17,243.01.1360374.356049.9828424.33884189.92326.08596515.985949.47270.97561372.5343%
1207,10,14,28,365,32/8,29/11,153.01.5267374.356067.1750441.53104189.93126.17227316.072268.70470.96571766.2734%
1309,14,28,31,339,13/24,26/3,143.01.9066374.356083.8918458.24784189.93904.13218094.032177.13100.96643161.1417%
1407,9,14,28,3214,32/9,28/26,303.02.3010374.3560101.2425475.59854189.94711.59588901.495899.96590.96191356.5935%
1509,14,28,31,3313,29/15,24/11,323.02.6842374.3560118.1055492.46154189.95496.36119686.2611111.71570.96023152.7777%
1609,14,28,31,3326,31/12,15/24,293.03.0905374.3560135.9813510.33734189.96328.259610,518.1596146.48360.95003149.2572%
1709,14,28,30,334,31/8,12/17,283.03.4959374.3560153.8207528.17674189.97158.464611,348.3646180.42440.94553046.1828%
1807,10,14,27,329,30/12,29/14,173.03.8918374.3560171.2403545.59634189.97969.135112,159.0351204.82540.94053243.5298%
1907,9,14,28,3614,30/12,17/24,323.04.2900374.3560188.7605563.11654189.98784.486212,974.3862231.51210.94271341.1522%
2007,9,14,27,3215,17/10,24/13,33.04.6822374.3560206.0161580.37214189.99587.519513,777.4195252.18330.93253139.0514%
1 PG: power generation; 2 MBV: minimal bus voltage; 3 BN: bus number.
Table 12. Scenario 2: system reconfiguration results of auto-select DG combination under 6 DGs and increasing power demand.
Table 12. Scenario 2: system reconfiguration results of auto-select DG combination under 6 DGs and increasing power demand.
PD 1 (%)Open
Switches
DG Locations (DE/MT/FC)PG 2 (MW)ECC ($/h)EMC (lb/h)Ploss (kW)MBV 3
(pu)
BN 4PL (%)
DGGridDGGridTotalDGGridTotal
1007,10,14,28,30--/6/8,21,24,28,323.10.6620109.106529.1268138.23333652.95881355.49395008.452746.97220.96773082.4036%
1109,14,28,32,33--/11,26/5,24,30,323.20.9362115.120341.1917156.31204071.80961916.96905988.778649.67540.97901377.3662%
1207,9,14,28,36--/5,10,16/7,24,323.31.2242121.134053.8630174.99704490.66042506.66156997.321966.15910.96871372.9417%
1307,10,14,28,36--/8,12,25,28/15,323.41.5036127.147866.1597193.30744909.51123078.92027988.431474.12860.96701769.3364%
1407,9,14,28,32--/8,12,25,31/16,283.41.8955127.147883.3999210.54774909.51123881.24278790.753994.45270.96463064.2060%
1509,13,17,28,33--/6,9,15,24,29/173.52.1795133.161595.8995229.06105328.3624462.94569791.3076107.03450.96363161.6248%
1607,10,13,17,28--/8,15,24,31/11,293.42.6859127.1478118.1781245.32594909.51125499.739710,409.2509141.86570.95301755.8672%
1707,11,28,31,3415/24,32/9,13,293.13.3938244.7381149.3270394.06514130.85486949.338311,080.1931178.29500.93443147.7379%
1807,10,14,28,3614/12,16,29/25,313.23.6872250.7519162.2360412.98794549.70567550.094912,099.8005200.18190.94191746.4631%
1907,11,14,28,3613/9,16,29/25,323.24.0825250.7519179.6307430.38264549.70568359.603812,909.3094224.01560.93341743.9409%
2007,11,14,28,3616/12,14,28/31,323.24.4733250.7519196.8256447.57754549.70569159.817613,709.5232243.30980.93961041.7030%
1 PD: power demand; 2 PG: power generation; 3 MBV: minimal bus voltage; 4 BN: bus number.
Table 13. Scenario 3: system reconfiguration results of auto-select DG combination under maximum DG quantity and increasing power demand.
Table 13. Scenario 3: system reconfiguration results of auto-select DG combination under maximum DG quantity and increasing power demand.
PD 1 (%)Open
Switches
DG Locations
(DE/MT/FC)
PG 2 (MW)ECC ($/h)EMC (lb/h)Ploss (kW)MBV 3 (pu)BN 4PL (%)
DGGridDGGridTotalDGGridTotal
1008,14,28,31,3323,30/1,5,7,11,243.70.0603132.30242.6546134.95704610.8276123.53894734.366545.33170.97473198.3956%
11010,14,28,32,3325/2,7,12,15,20,24,304.10.0341143.47081.5014144.97224730.994869.87084800.865647.62230.97803299.1746%
1209,14,28,31,33--/1,5,10,14,17,21,24,26,304.50.0195154.63920.8582155.49744851.16239.94094891.102961.50560.97173199.5684%
1309,14,28,31,336,10,23/13,15,19,25,29,324.80.0961172.68044.2303176.91076107.7144196.86806304.582466.64290.97883198.0364%
1409,14,28,31,3320,24/3,4,13,15,23,26,28,325.20.0830183.84883.6520187.50086227.8816169.95396397.835581.99910.97393198.4289%
1509,14,27,30,337/1,2,6,10,12,17,20,24,28,315.60.0601195.01722.6438197.66096348.0488123.03496471.083787.58560.97393098.9384%
1609,14,28,31,3315,23,13,8,19,3,29,24,2,9,31,256.00.0500206.18562.2006208.38626468.216102.41046570.6264106.01340.97003099.1733%
1707,9,14,27,3112,23,28,29/4,7,10,16,17,19,23,256.40.0429230.24051.8885232.12918143.619287.88738231.5065127.42090.96493199.3338%
1807,10,13,28,363,11,29/1,6,9,14,16,20,23,24,25,306.80.0279241.40891.2267242.63568263.786457.08768320.8740140.87940.95221799.5917%
1907,9,14,28,3211,31/3,4,6,7,13,14,18,19,23,26,28,327.20.0076252.57730.3324252.90988383.953615.47078399.4243149.05530.95681799.8952%
20011,28,32,33,341,2,5,7,9,12,15,17,19,21,23,24,27,29,317.50.0890257.73203.9167261.64868085.27182.27358267.5435159.01550.96073298.8270%
1 PD: power demand; 2 PG: power generation; 3 MBV: minimal bus voltage; 4 BN: bus number.
Table 14. Comparison of reconfiguration results under the deployment of 3 DGs.
Table 14. Comparison of reconfiguration results under the deployment of 3 DGs.
Reconfiguration without DGLiterature ResultsProposed
Ref. [39]Ref. [40]Ref. [41]Ref. [38]
Open switches7,9,14,32,377,10,14,28,327,10,14,28,327,10,14,28,317,10,14,28,327,11,14,28,30
Size of DG (MW)--0.5258 (32)0.5258 (32)0.6311 (18)0.5315 (18)0.6 (7)
(Bus number)--0.5586 (31)0.5586 (31)0.5568 (32)0.6158 (29)0.6 (28)
--0.5840 (33)0.5840 (33)0.5986 (22)0.5367 (32)0.6 (31)
PTL (kW)139.5573.7773.0572.2367.1162.18
MBV (pu)0.93780.97000.97000.97240.9713 (14)0.9661 (30)
Penetration level of DG--44.02%44.15%47.17%45.33%47.65%
Considered objectivesPTLPTLPTLPTLPTL and voltage stabilityPTL, ECC, EMC, and VD
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, X.; Srinivasan, T.; Kim, H.; Ra, I.-h. Exploration of DG Placement Strategy of Microgrids via FMFO Algorithm: Considering Increasing Power Demand and Diverse DG Combinations. Energies 2020, 13, 6568. https://doi.org/10.3390/en13246568

AMA Style

Wang X, Srinivasan T, Kim H, Ra I-h. Exploration of DG Placement Strategy of Microgrids via FMFO Algorithm: Considering Increasing Power Demand and Diverse DG Combinations. Energies. 2020; 13(24):6568. https://doi.org/10.3390/en13246568

Chicago/Turabian Style

Wang, Xiancheng, Thiruvenkadam Srinivasan, Hyuntae Kim, and In-ho Ra. 2020. "Exploration of DG Placement Strategy of Microgrids via FMFO Algorithm: Considering Increasing Power Demand and Diverse DG Combinations" Energies 13, no. 24: 6568. https://doi.org/10.3390/en13246568

APA Style

Wang, X., Srinivasan, T., Kim, H., & Ra, I. -h. (2020). Exploration of DG Placement Strategy of Microgrids via FMFO Algorithm: Considering Increasing Power Demand and Diverse DG Combinations. Energies, 13(24), 6568. https://doi.org/10.3390/en13246568

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop