Resilience-Oriented Framework for Microgrid Planning in Distribution Systems

: Recently, it has been suggested that microgrids (MGs) can improve the resilience of distribution systems. However, predictions about future faults are uncertain. This makes calculating the exact value of the beneﬁts of system resilience enhancement close to impossible at the time of MG planning. Therefore, this paper proposes a framework for MG planning, which focuses on resilience estimation. To consider the uncertainties of future failure events, the proposed method for estimating the resilience utilized the Monte Carlo simulation. In addition, an optimal scenario was estimated using a cost–beneﬁt analysis and constraints on the expected value of resilience enhancement. In the case study, an actual MG installation at D-university was evaluated to obtain the optimal MG planning scenario. The results show that the capacity and installation locations of the distributed generators (DGs) impact the resilience enhancement. The proposed method can effectively derive the optimal MG planning scenario by evaluating the possibility of future operations based on the segmentation of both the system conﬁguration and type of DG to improve the resilience of distribution systems. an optimal scenario. It shows that the NDG and ESS can contribute to improving the resilience in a distribution system when the ESS capacity is more than twice the PV capacity.


Introduction
Power system resiliency is defined as the ability to prepare for, adapt to, and recover rapidly from disruptive events [1]; it considers high-impact low-frequency events such as natural disasters and deliberate attacks [2]. Even a single occurrence of an event such as a natural disaster can cause critical problems in the power system. However, it is impossible to reduce the intensity of a natural disaster [3]. For this reason, the resilience enhancement of power systems is being increasingly researched.
In previous studies on resilience assessment, frameworks for resilience assessment with various goals were proposed to respond to natural disasters [4], including a specific disaster (typhoon) [5]. Moreover, the importance of a resilience metric for power systems has been highlighted [6]. Especially, rapid recovery is considered the most important feature of a resilience metric [7]. In this regard, microgrids (MGs), which operate the distributed generators (DGs), have been projected as a potential solution for supporting the rapid recovery of the resilience of a distribution system [8].
Existing studies on MG planning based on economic feasibility have considered the problem of DG sizing [9,10]. However, MG planning considering the resilience must involve the assessment process of the resilience of a distribution system to disruptive events. Hence, not only long-term planning but also short-term planning and operations should be reflected in the framework for MG planning [11].
In previous studies on resilience-oriented operation with MG, the categorization of DGs and sectionalization of systems were established as methods for rapidly recovering 2 of 17 from failures [12]. An optimization method based on the maximum allowable load loss was proposed for minimizing the investment costs [13].
In addition, a three-level model based on the multi-objective optimization problem was proposed in a previous study on resilience-oriented long-term planning with MG [14]. In the model, load profile uncertainty and some contingencies were considered, and the results proved a trade-off relationship between cost and resilience. In another study, multiple objective functions that maximize the critical load and minimize the capacity of dispatchable DGs were proposed for cost-effective MG planning [15]. A framework that incorporated resilience was developed for the economic assessment of MGs [16]. The cost-benefit analysis (CBA) in this framework comprised three steps and provided benefits of enhancement in the resilience and reliability through MG installation.
As mentioned in previous studies, to realize the success of long-term planning efforts, including MG installation, it is crucial to establish a framework that integrates various methods, such as that for system configuration, with the measurement criteria including that for the critical load. Nevertheless, owing to the uncertainty of events (type and intensity of disasters, number and location of failures, and time of failure occurrence and repair for each system), the uncertainty of various contingencies rather than the load profile uncertainty should be preferentially reviewed in the framework.
To address the problem of power distribution system reinforcement while avoiding the above issues, the resilience considering uncertainties was probabilistically evaluated using a Monte Carlo simulation (MCS) and the results were used as constraints for feasibility analysis [17]. To apply these prior studies to a real-world system, a more intuitive presentation of the effect of resilience enhancement is required for the system operator, and various types of components, such as DGs, should be applied in the system. Thus, this study proposes a resilience-oriented framework for MG planning. The framework consists of a probabilistic resilience estimation method, expected value of resilience enhancement (EoRE) estimation method, and CBA method. In the process, the post-contingency uncertainties are considered, and the various output characteristics of the DGs are reflected. As a case study, the MG planning of D-university was conducted using the proposed framework. The proposed method provides both the EoRE and CBA results, enabling the system operators in charge of MG planning or the MG operators to make an informed and reasonable judgment of the amount of resilience enhancement that should be set as a priority.

Framework of Microgrid Planning
MG planning is generally conducted based on the economic evaluation using load analysis, power generation forecasting, and system configuration [9,10]. In the case of MG planning considering the reliability, formulaic indices such as system average interruption duration index (SAIDI) and loss of load expectation (LOLE) are used [18]. In addition, the economic value of the reliability is estimated in terms of the outage costs per year according to the type of customers using the customer damage function [19]. However, owing to the uncertainty of contingencies in the future, determining a single value that reflects the benefit of resilience enhancement in MG planning is difficult. Furthermore, frameworks for MG planning including the generalized resilience indices are under development in various viewpoints.
With the objective of reducing the investment required for resilience enhancement in a contingency, this paper proposes a framework for MG planning, which comprises three steps, as shown in Figure 1. In step I, the various scenarios are established according to the purpose of the MG planner. In this step, the capacity as well as the type and location of the DG should be considered for each scenario, not only to maximize the economic benefits but also to improve the resilience of the system. The specific information of the established scenario is used to estimate the resilience enhancement in step II and cost-effectiveness in step III. In step II, the probabilistic evaluation of the resilience is conducted. In step III, the results of the probabilistic resilience are utilized to select the scenarios that satisfy the resilience constraints, and the optimal scenario is provided based on the commonly used CBA.

Estimation Method for Probabilistic Resilience
The probabilistic estimation of resilience, mentioned in step II of Figure 1, is illustrated in Figure 2 [17]. This method utilized MCS to derive the results of the probabilistic resilience, considering the post-contingency uncertainty in the development phase of MG planning. MCS is an inferential statistical technique that estimates an unknown value based on several simulations. It is used to derive a specific numerical value or probability distribution through the repeated random sampling of random variables. This study sets several random variables (e.g., location and occurrence time of failure and recovery time) for the event uncertainties to derive a probability distribution of the system's resilience performance.  In step I, the various scenarios are established according to the purpose of the MG planner. In this step, the capacity as well as the type and location of the DG should be considered for each scenario, not only to maximize the economic benefits but also to improve the resilience of the system. The specific information of the established scenario is used to estimate the resilience enhancement in step II and cost-effectiveness in step III. In step II, the probabilistic evaluation of the resilience is conducted. In step III, the results of the probabilistic resilience are utilized to select the scenarios that satisfy the resilience constraints, and the optimal scenario is provided based on the commonly used CBA.

Estimation Method for Probabilistic Resilience
The probabilistic estimation of resilience, mentioned in step II of Figure 1, is illustrated in Figure 2 [17]. This method utilized MCS to derive the results of the probabilistic resilience, considering the post-contingency uncertainty in the development phase of MG planning. MCS is an inferential statistical technique that estimates an unknown value based on several simulations. It is used to derive a specific numerical value or probability distribution through the repeated random sampling of random variables. This study sets several random variables (e.g., location and occurrence time of failure and recovery time) for the event uncertainties to derive a probability distribution of the system's resilience performance. In step I, the various scenarios are established according to the purpose of the MG planner. In this step, the capacity as well as the type and location of the DG should be considered for each scenario, not only to maximize the economic benefits but also to improve the resilience of the system. The specific information of the established scenario is used to estimate the resilience enhancement in step II and cost-effectiveness in step III. In step II, the probabilistic evaluation of the resilience is conducted. In step III, the results of the probabilistic resilience are utilized to select the scenarios that satisfy the resilience constraints, and the optimal scenario is provided based on the commonly used CBA.

Estimation Method for Probabilistic Resilience
The probabilistic estimation of resilience, mentioned in step II of Figure 1, is illustrated in Figure 2 [17]. This method utilized MCS to derive the results of the probabilistic resilience, considering the post-contingency uncertainty in the development phase of MG planning. MCS is an inferential statistical technique that estimates an unknown value based on several simulations. It is used to derive a specific numerical value or probability distribution through the repeated random sampling of random variables. This study sets several random variables (e.g., location and occurrence time of failure and recovery time) for the event uncertainties to derive a probability distribution of the system's resilience performance.   As depicted in Figure 2, the probabilistic estimation process of resilience comprises three sub-steps of step II. In particular, the resilience curve is estimated in step II-1, considering the fault and recovery information for post-contingency. Thereafter, the resilience index (RI) is calculated in step II-2 for resiliency evaluation, and the RI results are fitted to the empirical distribution function (EDF) in step II-3. The MCS is completed when the EDF(RI) satisfies the convergence range ( ). The MCS was applied to all the MG planning scenarios (j T ).

Estimating the Resilience Curve (Step II-1)
In this study, the resilience curve was defined as the "time-varying performance of the system." An example of the estimation method used to establish the resilience curve is shown in Figure 3. As depicted in Figure 2, the probabilistic estimation process of resilience comprises three sub-steps of step II. In particular, the resilience curve is estimated in step II-1, considering the fault and recovery information for post-contingency. Thereafter, the resilience index (RI) is calculated in step II-2 for resiliency evaluation, and the RI results are fitted to the empirical distribution function (EDF) in step II-3. The MCS is completed when the EDF(RI) satisfies the convergence range ( ). The MCS was applied to all the MG planning scenarios ( ).

Estimating the Resilience Curve (Step II-1)
In this study, the resilience curve was defined as the "time-varying performance of the system." An example of the estimation method used to establish the resilience curve is shown in Figure 3.   Figure 3b illustrates an example to relate the fault situation in a distribution system to the resilience concept. In Figure 3b, the locations of the DGs were set to nodes 3 and 8 in the preparation process of long-term planning, and the failure locations were assumed to be lines 4-9 and   Figure 3b illustrates an example to relate the fault situation in a distribution system to the resilience concept. In Figure 3b, the locations of the DGs were set to nodes 3 and 8 in the preparation process of long-term planning, and the failure locations were assumed to be lines 4-9 and 3-4 in the withstand process. In the adapt process, the system was sectionalized via circuit breakers (CBs). DGs were contributed to supply power to each section. Figure 3c describes the resilience curve for Figure 3b. The diagram shows the system's resilience performance at each of the following time points: the moment of equipment failure at the time of the event (t 1 , t 3 ), the moment when the system is sectionalized via a CB (t 2 , t 4 ), the moment when power is supplied using the DG in each section (t 5 ), and the moment when the equipment is restored through repair or replacement (t 6 , t 7 ).
To estimate the resilience curve in this manner, the sectionalization of the system in event situations and a method to derive the nodes where power can be supplied using a DG are needed. This study proposes the resilience curve estimation process shown in Figure 4. 3-4 in the withstand process. In the adapt process, the system was sectionalized via circuit breakers (CBs). DGs were contributed to supply power to each section. Figure 3c describes the resilience curve for Figure 3b. The diagram shows the system's resilience performance at each of the following time points: the moment of equipment failure at the time of the event ( 1 , 3 ), the moment when the system is sectionalized via a CB ( 2 , 4 ), the moment when power is supplied using the DG in each section ( 5 ), and the moment when the equipment is restored through repair or replacement ( 6 , 7 ).
To estimate the resilience curve in this manner, the sectionalization of the system in event situations and a method to derive the nodes where power can be supplied using a DG are needed. This study proposes the resilience curve estimation process shown in Fig  In this study, the topology of a particular distribution system is defined as a set of nodes expressed in Equation (1): where represents the nth node, and the pre-failure topology of a single power system is a set that comprises all the nodes as elements.
When the events occur in a particular distribution system, the sections are separated via CBs and can be expressed as a subset of the topology as follows: where represents the sth section of the topology separated by an event.
In Figure 4, the fault condition consists of the number, type, and location of the failed equipment and the failure and recovery time for each piece of equipment. Based on this condition, the nodes ( ) included in each section are organized according to time (t), as depicted in Figure 5a. In this study, the topology χ of a particular distribution system is defined as a set of nodes expressed in Equation (1): where N n represents the nth node, and the pre-failure topology χ of a single power system is a set that comprises all the nodes as elements.
When the events occur in a particular distribution system, the sections are separated via CBs and can be expressed as a subset of the topology χ as follows: where χ S represents the sth section of the topology separated by an event.
In Figure 4, the fault condition consists of the number, type, and location of the failed equipment and the failure and recovery time for each piece of equipment. Based on this condition, the nodes (N n ) included in each section are organized according to time (t), as depicted in Figure 5a.
As depicted in Figure 5, the nodes (N n ) in each section (χ S (t)) are physically connected but power outages may occur owing to the loss of power supply. This means that the available nodes (ANs) in each section χ S (t) along with the time-variant section χ AN S (t) comprising the ANs, shown in Figure 5b, must be determined. As depicted in Figure 5, the nodes ( ) in each section ( ( )) are physically con nected but power outages may occur owing to the loss of power supply. This means tha the available nodes (ANs) in each section ( ) along with the time-variant section ( ) comprising the ANs, shown in Figure 5b, must be determined. Thus, the objective function was applied as expressed in Equation (3) with the con straints expressed in Equations (4)- (9). The proposed method evaluates the adequacy o the capability (available capacity) by using the peak load for each node and the capacit of each DG [20]. These equations are equally applied to each ( ) to select the optima ( ) containing the ANs.
Equation (3) represents the objective function used to determine the AN value tha maximizes the sum of the resilience performance (∑ ∈ ) for ( ). In this case, represents the performance metrics, which can be applied specifically according to th measurement criteria that the system operator intends to evaluate, such as the peak load of the AN, number of customers, and number of critical customers of AN.
In Equation (4), ∑ ∈ is the sum of the peak loads of the ANs and must be les than or equal to the sum of the capabilities ( + ∑ ∀ ) of the feeder and DG. Thus, the objective function was applied as expressed in Equation (3) with the constraints expressed in Equations (4)- (9). The proposed method evaluates the adequacy of the capability (available capacity) by using the peak load for each node and the capacity of each DG [20]. These equations are equally applied to each χ S (t) to select the optimal χ AN S (t) containing the ANs.
ω NDG = at NDG 24 (8) Equation (3) represents the objective function used to determine the AN value that maximizes the sum of the resilience performance ( ∑ N∈AN R ξ N ) for χ AN S (t). In this case, ξ represents the performance metrics, which can be applied specifically according to the measurement criteria that the system operator intends to evaluate, such as the peak load of the AN, number of customers, and number of critical customers of AN.
In Equation (4), ∑ N∈AN L N is the sum of the peak loads of the ANs and must be less than or equal to the sum of the capabilities (P f eeder + ∑ ∀DG P DG ) of the feeder and DG.
In Equation (4), if only the non-controllable DGs (NDGs) dependent on the natural environment are present in a specific section, e.g., photovoltaic (PV) and wind turbines (WT), they can cause an interruption in future operations. Equation (5)  characteristics of a DG divided into NDG and controllable DG (CDG), and the energy storage system (ESS) is included separately as it can charge as well as discharge power.
For Equation (6), the capability of each NDG (P NDG ) is calculated using ω NDG ·P NDG when the energy produced from NDG (E NDG ) [kWh] can be controlled using E ESS [kWh] during the time between the occurrence of disaster and recovery (T). However, if this is not feasible, P NDG is calculated as "0," assuming that curtailment occurs by manipulating the inverter for each NDG. Equation (7) depicts the evaluation of E NDG used in Equation (6). As expressed, E NDG can be calculated by using the weighting factor (ω NDG ) and time of occurrence of failure (T), where ω NDG is calculated using the average daily production time for NDG (at NDG ), as shown in Equation (8). In Equation (9), P ESS is used when the dischargeable capacity over time for ESS (E ESS /T) is less than the power conversion system capacity of ESS (P PCS ). Otherwise, P PCS is used.
R ξ T (t) from Equation (10) represents the time-varying performance of the system, whereas R ξ N (t) represents the time-varying performance of each node.

Calculating the Resilience Index (Step II-2)
The resilience index (RI ξ ) used in this study is expressed in Equation (11) [21].
In Equation (11), R ξ T (t) and R ξ T0 (t) represent the resilience curve assuming that a failure occurred and did not occur, respectively.

Fitting EDF for RI (Step II-3)
As mentioned above, RI ξ comprises the various sets (RI ξ,j ) expressed in Equation (12) according to the criteria for performance measurement and the scenarios of MG planning.
RI ξ,j = RI 1,ξ,j , . . . , RI i,ξ,j The set of RI is fitted as an EDF. This study expresses the probability density function (PDF) as PDF RI ξ,j . The cumulative distribution function (CDF) is calculated as shown in Equation (13) [17,22].

Optimal MG Planning Considering Resilience Based on Cost-Benefit Analysis
To determine the optimal MG planning scenario, this section presents the CBA and resilience constraints. As this study considers the resilience of the distribution system during the MG planning process, the costs and benefits in the CBA are applied from the perspectives of the utility and MG operators.

Cost-Benefit Analysis
In this study, the method for determining the optimal MG planning scenario is described based on the benefit-cost (B-C) ratio among the various methodologies for CBA. Generally, there is an investment value of B-C ratios greater than or equal to 1, and there is a decreased investment value when the B-C ratio is less than 1.
In this study, the scenario achieves the maximum B-C ratio, which is expressed as Equation (14).
The elements of B and C in Equation (14) are expressed in Equations (15)-(16); these are estimated from parameters such as life cycle and power output forecasting for each DG [9,10].
In Equation (15), the benefit (B) is categorized into the benefit from the perspective of the environment and those from the perspectives of the MG operators and utilities, and the calculations for each benefit category are expressed in Equations (17)- (19). On the contrary, the cost (C) presented in Equation (16) (20) and (21).

Resilience Constraint
As the constraints of power balancing and DG capacity were considered during resilience calculation, this section focuses on the resilience constraints.
In this study, the effect of resilience enhancement through MG installation is calculated as in Equation (22) below: where EoRE denotes the expected value of the resilience enhancement, and E RI ξ,j indicates the expected value of the RI for the respective MG installation scenarios. E RI ξ,Re f . denotes the comparison criteria for E RI ξ,j , which can be set considering the objectives, e.g., prior to MG implementation or increased DG capacity for MG. Using the CDF(RI) that was derived above, EoRE can be calculated as shown in Equation (23) below [23,24].
In Equation (23), CDF 1 (RI) denotes (CDF RI ξ,Re f . ), a function serving as the comparison criteria, and CDF 2 (RI) indicates the function of each scenario (CDF RI ξ,j ). An example of EoRE based on Equation (23) is presented in Figure 6 below. , ≥ In Equation (19), CR represents the constraint of resilience, which is the minimum allowable value of EoRE obtained through discussions between the utility and MG operators. Since the scenarios of the MG planning that satisfy the CR are provided, the proposed framework can avoid the calculation process of the benefit (B) for resilience enhancement in CBA.    Figure 6 shows that the proposed method allows the estimation of the difference from the expected value of the RI, thus providing intuitive results for the degree of resilience enhancement through MG installation.

System configuration
The resilience constraint is set as in Equation (24).
In Equation (24), CR represents the constraint of resilience, which is the minimum allowable value of EoRE obtained through discussions between the utility and MG operators. Since the scenarios of the MG planning that satisfy the CR are provided, the proposed framework can avoid the calculation process of the benefit (B) for resilience enhancement in CBA. Figure 7 is the system configuration at D-university in Korea with the MGs in operation. The information regarding the peak load for each node in the system (Figure 7) is organized in Table A1 (Appendix A). The resilience constraint is set as in Equation (19).

System Configuration
In Equation (19), CR represents the constraint of resilience, which is the minimum allowable value of EoRE obtained through discussions between the utility and MG operators. Since the scenarios of the MG planning that satisfy the CR are provided, the proposed framework can avoid the calculation process of the benefit (B) for resilience enhancement in CBA.

Scenarios of MG Planning (Step I)
The scenarios for MG planning, which consider the capacities of the DGs as well as their installation locations, are listed in Table 1. Table A1 (Appendix A) shows the DG capacity for each node in scenarios 10, 11, 15, and 16, assuming that the DGs in Table 1 are installed at different locations.  Table 2 shows the ratios calculated to minutely review the distribution of the DGs. Each ratio was calculated based on the following criteria: a : Ratio of DG capacity to peak load b : Ratio of controllable DGs (CHP) to peak load c : Ratio of controllable DGs (CHP) to total DGs In Tables 1 and 2, scenarios 1-3 assume the same DG capacity but different DG installation locations (nodes 2, 10, or 18). Scenarios 4-6 have the same total DG capacity as that in scenarios 1-3 but with different ratios of the DG composition. Scenarios 7-11 have the same total DG capacity as that installed in D-university. Scenario 10 represents the actual MG in D-university, and it is the calculation result of the DG capacity to maximize the economic feasibility without consideration of resilience. To ensure resilience enhancement, the total DG capacity in scenarios 12-16 is larger than that in scenarios 7-11, and a large ratio of the CDG to peak load is assumed.

Simulation Condition for Estimating the Probabilistic Resilience (Step II)
To consider the various post-contingency uncertainties, the simulation conditions were examined with the following limitations [25]: The unit price for each cost and benefit element is presented in Table A2 (Appendix A). Figure 8 shows the E-CDF result for each scenario of the three failed facilities. Table 3 shows the EoRE results for one to four failed facilities. The intermediate results (resilience curve) obtained for calculating the EoRE are presented in Figure A1 (Appendix A). The unit price for each cost and benefit element is presented in Table A2 (Appendix) Figure 8 shows the E-CDF result for each scenario of the three failed facilities. Table  3 shows the EoRE results for one to four failed facilities. The intermediate results (resilience curve) obtained for calculating the EoRE are presented in Figure A1 (Appendix).   In Figure 8, the E-CDF results for the scenarios with the same total capacity are similar, except for scenarios 12-16. This means that when a (the ratio of the DG capacity to peak load) in scenarios 1-11, as listed in Table 2, is small, the variation in the enhancement of resilience according to the installation locations of the DGs is negligible.

Resilience Enhancement
Meanwhile, although a of scenarios 7-11 is larger than a of scenarios 1-6 in Table 2, the enhancements in the resilience in scenarios 7-11 and 1-6 are similar. This is because the ratio of CDGs ( b , c ) in scenarios 7-11 in Table 2 is low, indicating that controllable DGs such as CHP contribute to enhancing the resilience more than the other DG types.
In Table 3, the variation in the EoRE according to changes in the number of failed facilities was significant for scenarios 12-16, which had relatively large DG capacities. Nevertheless, Scenario 12 showed a lower EoRE than scenarios 13-16 because the DGs in scenario 12 were located next to the feeder at node 2, as shown in Figure 7. These results indicate that when a failure occurs in the line connecting nodes 2 and 10 in the distribution system of Figure 7, it is more advantageous for the DGs to be installed at nodes 10 and 18 (scenarios 13 and 14) than at node 2 (scenario 12). The following resilience constraint was set to review the optimal scenario. The condition for the number of failed facilities was limited to three.

•
Resilience constraint (CR): 0.05 Figure 9 shows the EoRE and B-C ratio for all scenarios, as well as the baseline CR. Table 4 presents the B-C ratios for scenarios that satisfy the resilience constraint. Table 4 shows that scenario 8, where the capacity of the NDG was relatively large, was an optimal scenario. It shows that the NDG and ESS can contribute to improving the resilience in a distribution system when the ESS capacity is more than twice the PV capacity.
Energies 2022, 15, x FOR PEER REVIEW In Figure 8, the E-CDF results for the scenarios with the same total capacity ilar, except for scenarios 12-16. This means that when ⓐ (the ratio of the DG cap peak load) in scenarios 1-11, as listed in Table 2, is small, the variation in the enhan of resilience according to the installation locations of the DGs is negligible.
Meanwhile, although ⓐ of scenarios 7-11 is larger than ⓐ of scenarios 1-6 2, the enhancements in the resilience in scenarios 7-11 and 1-6 are similar. This is the ratio of CDGs (ⓑ, ⓒ) in scenarios 7-11 in Table 2 is low, indicating that con DGs such as CHP contribute to enhancing the resilience more than the other DG In Table 3, the variation in the EoRE according to changes in the number o facilities was significant for scenarios 12-16, which had relatively large DG ca Nevertheless, Scenario 12 showed a lower EoRE than scenarios 13-16 because the scenario 12 were located next to the feeder at node 2, as shown in Figure 7. Thes indicate that when a failure occurs in the line connecting nodes 2 and 10 in the dist system of Figure 7, it is more advantageous for the DGs to be installed at nodes 1 (scenarios 13 and 14) than at node 2 (scenario 12).

Optimal scenario based on the CBA and resilience constraint
The following resilience constraint was set to review the optimal scenario. T dition for the number of failed facilities was limited to three.

•
Resilience constraint ( ): 0.05 Figure 9 shows the EoRE and B-C ratio for all scenarios, as well as the base Table 4 presents the B-C ratios for scenarios that satisfy the resilience constraint.    However, when the reference value of the resilience constraint is set to be greater than 0.1 in Figure 9, scenario 15 (B-C ratio: 0.86) can be derived as an optimal result. This shows that the resilience (EoRE) enhancement and cost (B-C ratio) are in a trade-off relationship.
The proposed framework for MG planning can provide a cost-effective scenario by selecting scenarios that satisfy the conditions for enhancing resilience. The estimation method of probabilistic resilience included in this process provides the EoRE that considers not only the DG capacity and type but also the installation location.

Conclusions
This study developed a framework for MG planning considering the resilience of a distribution system as described in Section 2. In Section 3, the estimation method of the resilience curve considering the system configuration and output characteristics of the DG was presented, and the method for calculating the RI and EDF was described. In this process, the post-contingency uncertainties were considered based on the MCS. In Section 4, the CBA model was presented for various cost and benefit elements, and the EoRE was proposed to intuitively express the degree of resilience enhancement.
To consider the realistic situation of the distribution system, a case study was conducted at the D-university with an installed MG. The results in the case study showed that the proposed method for estimating the resilience could consider not only the capacity of the DGs but also their type and installation location. In addition, the proposed framework provided an optimal scenario that satisfied the resilience constraints (CR) and maximized the B-C ratio.
When the proposed framework is verified through various real systems, it could contribute to establishing a more resilient system using MG in long-term planning. Especially, the framework could be usefully utilized to establish MG of the distribution system that is oriented toward resilience enhancement, such as customers of hospitals and military and industrial facilities, because the scenarios of the MG planning that satisfy the resilience constraint are provided.