Providing a New Multiobjective Two-Layer Approach for Developing Service Restoration of a Smart Distribution System by Islanding of Faulty Area

One of the essential capabilities of a smart distribution network is to improve network restoration performance using the postfault islanding method. Islanding of the faulty area can be done ofine and online. Online islanding will decrease load shedding and operation cost. In this study, a novel two-step mathematical method for system restoration after the fault is presented. A new mathematical model for the optimal arrangement of the system for the faulty area in the frst layer is proposed. In this layer, the main objective is to decrease the distribution system’s load shedding and operational costs. In this regard, after the fault event, the boundary of the islanded MGs is determined. Ten, in the second layer, the problem of unit commitment in the smart distribution network is addressed. In addition to the load shedding, optimal planning of energy storage systems (ESSs) and nondispatchable distributed generation (DG) resource rescheduling are also determined in this layer. Te important advantages of the proposed approach are low execution time and operational costs. A demand response (DR) program has also been used for optimal system restoration. Solving the problem using the multiobjective method with the epsilon-constraint method is another goal of the paper, which simultaneously minimizes the cost and the emissions of the smart distribution network. Te proposed model has been tested on an IEEE 33-bus system. Better performance of the proposed model compared to the techniques in the literature has been proven.


Introduction
To develop renewable resources and reduce concerns about the problem of environmental pollution and the increasing use of the smart grid, government incentives aimed at developing distributed generations and implementing demand response programs have a signifcant efect.Some smart grid goals include the following: participation of distributed generation resources at the distribution network level, the use of electrical storage resources to participate in DR programs, responding to real-time loads and prices, selfhealing, and smart protection.
One of the important capabilities of smart grids is service restoration in the event of a fault in the power systems.Optimal restoration is carried out with the aim of minimizing load shedding.Also, all system constraints should be considered.Te system's ability to quickly detect faults, take necessary actions to reduce the adverse efects of faults, and quickly restore the system to a stable operating state is defned as self-healing [1].Te structure of the smart grid, while implementing the self-healing strategy, is categorized into smart transmission grid, smart distribution grid, and smart MGs [2].Te method of diferential evolution algorithm is chosen in [3] to optimize the performance of network MGs in self-healing mode.Te evolving search space method for new confgurations is performed by combining graph theory and a heuristic method [4].To restore service to the maximum customers afected due to the contingency, a fuzzy multiagent system was used in the distribution system [5].In [6], an adaptive restoration decision support system based on a two-step MILP problem has been used to address the challenging online restoration problem.In [7], a self-healing control strategy is proposed that includes fault detection, fault localization, faulted area isolation, and power restoration in the electrical distribution system.In [8], a voltage control method is presented based on a comprehensive architectural model consisting of several distributed photovoltaic resources and management systems that support network reconfguration after the selfhealing operation.A self-healing approach is presented in [9], which evaluates the robustness of networked MGs during the islanded mode.Network fexibility and resilience are essential in risky events such as bad weather conditions and cyberattacks.A mixed-integer linear optimization method is proposed to achieve self-healing operation to minimize load shedding [10].Concepts related to the selfhealing capability of smart grids include transmission, distribution, MGs, transient stability, and cyberattack [11].
In many studies, the process of energy distribution system management has been investigated.Te optimal model for self-healing management in active distribution systems, including electric vehicles (EVs), renewable energy sources, and DR programs, is presented in [10].In [12], with the aim of day-ahead planning of a smart distribution system, a comprehensive operating method has been made for normal and emergency conditions.In [13], the daily optimal scheduling problem of networked MGs using a metaheuristic algorithm under uncertainties of renewable energy systems and loads is investigated in a proposed energy management system.In [14], the authors have studied the energy consumption management of a rural MG in the optimization framework based on a new stochastic method to balance generation and demand.
Recently, DR programs have been converted into an essential interest for researchers.A variety of loads, such as industrial, commercial, and residential, in the DR programs, can help to introduce them as fexible and smart loads [15].Te authors in [16] have investigated the efect of DR programming simultaneously with the multi-MG-based operation of smart distribution systems.A step-wise DR program has been used in the energy management system for an isolated structure of networked MGs.Scenario-based analysis has dealt with the uncertainty of renewable energy sources and loads [17].In [18], an ESS, an MG consisting of a combined heating and cooling system, photovoltaic power generation, and the response load are implemented to check the optimal scheduling process of these units.In [19], the efectiveness of DR programs and local load energy supply by microgrid based on renewable resources has been studied.Te authors in [20] present a method to evaluate the way of real-time energy management in smart grids.A multiobjective approach is used to implement a DR program considering the power grid's load factor and residential users' energy costs [21].In [22], an energy management toolbox for the buildings connected to the electric grid has been modeled.
Various studies have dealt with faults in the self-healing mode and focused on island building.In [23], the efect of EVs on improving the self-healing properties in an islanding mode in the power grid is evaluated.A real-time self-healing scheme with measurement and model-based algorithms has been proposed [24] to deal with severe power system disturbances.
An MG is a structure that can be connected or disconnected from the main grid.MGs are usually switched to the islanded mode due to economic problems, maintenance purposes, and network faults.When an outage occurs, MGs can operate in islanded mode or grid-connected mode.MGs downstream of the grid can act as a standalone source for their customers by disconnecting from the network [25].In [26], an optimal self-healing method for MG islanding under diferent scenarios has been investigated, an islanded MG optimization scheduling focused on frequency adjustment of units, and DR is proposed.A novel concept of conventional droop control is developed as a fast droop controller in [27] cooperation with a modern frequency controller to ensure the reliability of the microgrid system.In [28], an optimal integral proportional control method is used in an island microgrid.Te diferent categories for intentional islanding of MG at the point of standard coupling have been analyzed in [29].In the self-healing scheme for the reconfguration of smart distribution networks with the presence of DGs, a heuristic algorithm with a tree structure is used [30].An energy management system with convex relaxation has been developed in [31] for an islanded MG that optimizes its operating cost.An optimal method for energy management is used in [32], which increases the reliability of an MG structure and considers the proft of small network owners.A comprehensive study of diferent MG control methods, such as hierarchical control of the islanded MGs, has been carried out in [33].
Tere are various methods to optimize the energy management of MGs by considering the technical constraints, which can be mentioned using mathematical approaches or heuristic optimization techniques.A mathematical approach for enhancing the resiliency of the MG energy management system, based on islanded mode, has been developed [34].In [35], a mathematical model for renewable energy penetration in an island MG, in the presence of DR and ES, is proposed to mitigate the imbalance of supply and demand due to the intermittent nature of renewable production.Te authors in [36] provide an overview of optimization and mathematical modeling of the DR algorithm and its implementation at various levels in the smart grid.Te study in [37] proposes a performance analysis of power system parameters, including the voltage, frequency, real and reactive power, and phase angle, to detect the island using mathematical morphology.In [38], a new strategy for MG protection based on mathematical morphology is introduced.Te proposed scheme was investigated for islanded systems and grid connected with loop and radial confgurations.In [39], a two-step method is used that, in the frst layer, solves a mixed-integer linear programming (MILP) problem and, in the second layer, solves a nonlinear programming (NLP) problem.In [40], a twostage stochastic program has been used, which optimizes the output power of WTs and PVs and load consumption in the frst stage and adjusts production in the second stage with 2 International Transactions on Electrical Energy Systems suitable scenarios based on the output power of renewable sources.In [41,42], a two-step method is used.In the frst step, graph theory integrated with the binary particle swarm optimization (PSO) algorithm is utilized to determine the optimal arrangement of the distribution system.In the second step of a nonlinear model, the unit commitment problem is solved mathematically [41].Te models presented in [39][40][41][42] have a very high execution time due to their nonlinearity, and it is practically impossible to use these models in large systems.On the other hand, due to intelligent algorithms, achieving the optimal solution will not be guaranteed.However, in the present study, a new linearized two-step method is used.Te proposed model is solved using mathematical solvers.Te proposed model also uses a load-shedding tool integrated with DR schemes, while, in [39][40][41][42], DR-based methods were not used.Te online island construction method retrieves the service in the fault zone in the proposed scheme.In the online islanding manner, the boundaries of island MGs and the number of MGs are determined optimally after the fault in a part of the distribution system.After the fault occurs, the fault line is isolated.Ten, the optimal island MG is determined by closing the tie switches and opening switches in the faulty area.In addition to creating an optimal island operation, the radial condition of the system is maintained.Terefore, several tie switches may be closed, and several switches may be opened.Of course, upstream network switches are not in planning optimal MG confguration in the faulty area.Because the faulty area will be connected to the upstream network, the system will no longer act as an island.While in the proposed method, the performance of the system in the faulty area is considered based on the formation of one/multiple islands.Terefore, in the frst layer, islanded MGs will be formed via the formation of optimal MGs using a linear mathematical model.In the second layer, the unit commitment problem of smart distribution systems is solved by the optimal arrangement formed in the frst step.At this stage, the unit commitment problem integrated with load-shedding tools and DR schemes is solved mathematically using a linear model.One of the advantages highlighted by the proposed method is using an optimal approach that can signifcantly decrease the execution time of the problem and achieve an optimal solution.Also, in the proposed model, all possible tools, such as smart load shedding, DR, various DG resources, and ESS, have been used for optimal restoration.
Another issue that we will discuss in this study is emissions.In multiobjective modeling, the increase in investors' profts, along with the reduction of environmental pollutants, has become more prominent.Terefore, a solution that only maximizes the total proft or minimizes the total scheduling costs may not be suitable for power grids alone.Tus, the multiobjective framework is proposed for solving the problem mentioned above, considering emissions.Te multiobjective method is used in articles, such as [43,44], to maximize profts and minimize emissions.Te authors in [43] investigate the role of the renewable-based VPP in maximizing proft and minimizing emissions in a two-objective manner.Te authors in [44] solve the bilevel problem with the augmented epsilon-constraint method in a biobjective way, which minimizes the emission of virtual power plant units and maximizes the proft.In the epsilonconstraint method, the primary objective function is designated the leading objective, while the second to n-th objectives are constrained to a specifc maximum value.Altering this value can result in multiple solutions, some of which may not exhibit efciency.It should be noted that none of the papers devoted to solving the problem of service restoration of a smart distribution network has incorporated multiobjective issues.
Te main components of the article are described as follows: (i) Providing a novel mathematical model for solving the reconfguration problem of a smart distribution system after a fault.
In literature, intelligent algorithms, graph theory, and nonlinear models have mainly been used in the frst layer.However, this paper uses a new linear model for reconfguration after the fault.Te proposed model leads to a decrease in the solution time and optimal solutions.(ii) Solving the unit commitment problem in a smart distribution system with smart load shedding and DR tools.
In the second layer, the unit commitment problem is solved by DG rescheduling and responsive loads.
Te proposed linear model leads to an optimal response at low solution times like the frst layer.Execution time is essential in the restoration problem.(iii) Creating optimal islanded MGs in the faulty area.
Another advantage of the proposed method is the creation of optimal MGs after the fault; while, in the methods with the islanding approach in the literature, MGs are formed ofine.Online formation of MGs is essential because, according to the available resources and the fault location, the boundaries of MGs are formed, which will lead to a reduction in load shedding of the distribution systems.(iv) Providing mathematical modeling and new objective functions for the restoration problem of the smart distribution system problem.(v) A multiobjective model is introduced to minimize the cost and emission in the smart distribution network based on the epsilon-constraint method.
In the following, two-layer linear mathematical modeling is used.Te formulation of the model is described in Section 3. Section 4 is related to case studies and results, and fnally, in the last section, the conclusion is presented.

Proposed Approach
Te assumptions for optimizing the problem according to the study of references [19,22,23,45]  In this paper, a two-step algorithm is presented.In the frst step, the restoration of the distribution system is achieved by a new linear model.In the second step, the unit commitment problem of distribution systems is solved by the optimal arrangement created in the frst step.Ten, the load-shedding tool is used on responsive loads.Rescheduling of dispatchable DG resources is also carried out.If the condition of balance is not met through the reduction of responsive loads and rescheduling of DG resources, the load will be curtailed.Of course, load shedding with high priority will impose more costs on the operation of distribution systems.Te proposed solution algorithm is shown in Figure 1.

Formulation
Tis step presents a novel objective function to solve the restoration problem.In objective function (1), term 1 is the penalty for voltage deviation from the allowable amount.In this term, the values are in the form of per unit.Te c delv is the voltage deviation coefcient penalty in dollars.Term 2 is the proft from the sale of electricity to customers.Terms 3 to 6 show the cost of load shedding of the frst, second, and third priorities and load outage, respectively.Power fow in the lines of each MG is described in term 7. x(k.j) � 0 does not mean the line will be cut, but rather that the current on the line can be zero.In other words, x(k.j) � 0 means the line between the buses and is a candidate of boundary lines between two self-adequate MGs with minimized generationload imbalance [40].
Tis item minimizes the total number of connected lines that lead to island MGs.
In the frst term, the voltage may deviate higher or lower than the nominal voltage.Terefore, the absolute value must be used.If the absolute value of the objective function is not used, the minimum allowable voltage should be used.Since the absolute value causes the model to be nonlinear, the absolute value must be converted to a linear relationship.In the frst term of equation ( 2), the value xp(k, t) + xn(k, t) is substituted.Tere are two positive variables.Tese two variables can be converted from a nonlinear model to a linear model by modeling the value of voltage deviation with the help of equation (3).
Te following equation shows the allowable voltage value on each bus: 4 International Transactions on Electrical Energy Systems Te active power fow equation is expressed in equation ( 5).Tis relationship is nonlinear and makes the solver not reach the desired answer.For this reason, the continuation of the formulation of this relation will be linear.
Due to the nonlinearity of the power fow equations in ( 5) and ( 6), these two relations are linearized as follows.To linearize, the following assumptions must be met: Terefore, the voltage of a bus can be expressed as follows: Voltage changes must be within the allowable range (∆V min ≤ ∆V k ≤ ∆V max ).Terefore, according to equations ( 5)- (7), they change as follows: Due to the multiplication of the two variables, equations ( 8) and ( 9) are still nonlinear.Since numerical results ∆V k,t .θkj,t and ∆V j,t .θkj,t are expected to be very insignifcant.Terefore, these nonlinear terms are removed from equation (8), which are fnally equation ( 10) linear relation of the power fow equation.
Te transmission power of each line exists if the switch is closed (X k,j � 1) and should be zero if the switch is open (X k,j � 0).For this reason, by using the large M method and multiplying this binary variable by the injected active power in relation (9), this condition will be guaranteed.
Te active power balance equation can be seen in the following equation: Te battery storage at low-load times helps to reduce the load shedding on the distribution system by storing energy and discharging during peak times.In the below, the limit of discharge and charge is represented.
Te battery storage can only be discharged or charged by the following equation: Te battery charging mode is mentioned in the following equation: Te minimum and maximum battery charge statuses are expressed as follows: Te operating range related to the output power of MTi at time t is given in (15).pg (k,t) and pg (k,t) express the minimum and maximum time-dependent operating range, respectively.It is not necessary that two variables are equal to pg min k,t and pg max k,t .

International Transactions on Electrical Energy Systems
Te lower and upper limits of operation pg (k,t) and pg (k,t) are described in ( 16) and (18).Te up/down level constraints of MT ramps are shown in relations ( 17) and (19).
In ( 20) and ( 21), the costs related to the shutdown and start-up of MTs are stated.Tese costs are used in the planning of MTs and the objective function in (1) and (2).
Te on-of status of unit i at time t is determined by u i,t .Te shutdown and start-up of MTs are also represented by z i,t and y i,t .On and of constraints of MTs are described in (22) and (23), respectively.However, these two relations make the model to a nonlinear model.To linearize the minimum shutdown limit, linearized relations ( 24)-(26) instead of (22) are used.In (24), the shutdown and start-up of the MT are limited to the MTs on-of state in the last time.Equation (25) shows that the MT turns of or on at time t, and both states cannot occur simultaneously.Equation (26) also indicates that u i,t , y i,t , and z i,t are binary variables.

􏽘
Equations ( 31)-( 34) also indicate the minimum shutdown time of each MT.
Te load power in each bus (p D i,t ) consists of the total power of home appliances based on diferent priorities as follows (35).
Te task of the load collector is to combine the power of the devices of diferent priorities on diferent buses and, at any moment, give the total energy of the distribution network in each priority to the DSM controller.Equation (36) shows that the whole load shedding by the priority on each bus is less than the total power of the priority device 1 of all the houses on each bus.Equations ( 37) and ( 38) also show the load-shedding constraints for priority loads 2 and 3, respectively.Equation (39) shows the load outage of uncontrolled loads.
However, the unit commitment problem in the new system must be resolved after the restoration.At this stage, planning is carried out for dispatchable DG and ESS resources.Also, the amount of load shedding is determined at this step.
In the objective function of the second layer, all the practical objectives in modeling the problem of unit commitment problem are considered.In the objective function (40), the frst to third terms show the cost of production and shutdown and start-up of MTs.Term 4 shows the proft of 6 International Transactions on Electrical Energy Systems customers' sales.Terms 5 to 8 show the cost of load shedding based on priorities 1 to 3 and load outage.Te cost of buying electricity from the upstream network is stated in term 9, and the proft from the sale of electricity to the upstream network is stated in term 10.Te generation cost of MTs is considered in some studies as a quadratic function.However, in [42], a linear model has been used for the production cost function of MTs due to computational load reduction.Terefore, in the objective function, instead of this nonlinear term, a linear term is used in relation.
In equation ( 41), the active power fow is represented.However, here, x(k, j) is constant and indicates the status of the system lines resulting from the frst step.Terefore, this relation no longer needs to be linearized by equation (9).
Te second step equations are equations ( 8) and ( 10)-( 40) that index i should be used instead of index k.

Multiobjective Modeling.
Te emission of the smart distribution systems is formulated based on a linear equation (42).Te total production of each unit is multiplied by a fxed factor to obtain the emission value of the unit, which is shown in Table 1 [43].
In the problem of multiobjective mathematical programming, in most cases, there is a confict between the objectives, that is, the improvement of the result of one of the objectives leads to the deterioration of the results of the other objectives.Terefore, the main diference between multiobjective and single-objective problems is that there is no single optimal solution in multiobjective problems that optimize all objectives simultaneously [44].

Epsilon-Constraint Method.
Te epsilon-constraint method is described below to solve a multiobjective optimization problem.
In the two-objective problem of emission and cost of the units, acts as the main objective function, which is the system cost, and emission (f 2 (x)) act as a constraint.In addition, x is an array of decision variables, which in this problem are the generation power of units.By changing the e 2 , Pareto solutions to the problem are generated (in this problem, there are 10 Pareto solutions) [45].To solve the epsilonconstraint method, the problem is solved as a single objective by considering the main objective function (cost).In this case, cost is at its lowest amount, and the emission is at its highest amount.Ten, emission is solved as a single objective.As a result of this, the emission is at its lowest value, and the cost is at its highest value.Ten, to generate Pareto solutions, the model is solved as a single-objective problem by considering the cost as a main objective function and emission as a subobjective function (constraint) as shown in equation (43).By changing the amount of e 2 , Pareto solutions are generated.

Fuzzy Decision-Making Method. By accessing all the
Pareto solutions after solving the problem, the decisionmaker should choose one of the Pareto solutions according to the prioritization and as the fnal solution.To choose the best answer for the decision-maker, it is suggested to use the fuzzy method with a linear membership function [43].Te proposed fuzzy decision-making method and its membership functions are defned as (44) and (45), used for maximization and minimization, respectively.Te output of objective functions in case of worst and best answers is arranged in the order of the nadir point f SN n and ideal point f u n .Te optimality of the objective function is indicated by μ r n in the r-th Pareto solution number.μ r n is the membership function f n in the r-th Pareto solution number, and f r n is the expression of the value of the objective function f n in the rth Pareto solution number.Te general membership function of the r-th Pareto solution is called μ r , which is calculated according to (46), where ω n is the importance factor of the n-th objective function.Decision-maker is the responsible for the values of the importance coefcients.For example, if environmental issues are the decision-maker's top priority, f 1 will be given a lower value, and if fnancial issues are more important, f 1 will be given a higher value [43].Te best solution is the one with the highest amount of μ r .

International Transactions on Electrical Energy Systems
(1) Maximization: (2) Minimization: Both models are mixed-integer linear programming problems and are solved in GAMS software with the powerful CPLEX solver.

Case Studies
Te simulation is performed on an IEEE 33-bus system to validate the proposed scheme, as shown in Figure 2 [41].Te information about the ESS and DG resources used is shown in Table 2 [41].In the studied system, the power of each residential house is assumed to be 10 kW.Te smart devices inside each house are divided into adjustable, interruptible, and deferrable (shiftable) devices.For home appliances, uncontrollable appliances are considered 0.8 kW.Also, 3.5 kW of the total power of equipment with frst priority, 3.1 kW of equipment with second priority, and 2.6 kW of equipment with third priority are assumed for each house [45].Table 1 shows the MT information.Information on the wind, solar, and load coefcients and energy prices are shown in Figures 3 and 4, respectively.Table 3 shows the costs associated with the model parameters.Problems are solved in two case studies, single objective and multiobjective.
4.1.Te Single-Objective Solution to the Problem.In this case, the emission is not considered.A fault between buses 4 and 5 is considered for testing the proposed model in the frst layer.Figure 5 shows the optimal system arrangement after fault between buses 4 and 5 (frst layer output).Since the aim is to island the faulty area, connecting tie switches 33, 35, and 37 is impossible.Because by connecting these switches, the MG no longer acts as an island.Four island MGs are formed in the faulty area after isolating the line between buses 4 and 5. Te frst MG is shown in pink, the second MG in red, the third MG in green, and the fourth MG in blue.Tie switches 34 and 36 are closed.Switches 7, 10, 15, 16, and 31 are also open.As seen in Figure 6, the radiality condition of the faulty area is maintained.According to the optimal arrangement in layer 1, the unit commitment problem is implemented in the second layer.Te output results for the second layer are shown in Table 4. Te frst point in Table 4 is that the MTs connected to the upstream network (19,25) have no output.According to the energy price shown in Figure 7 at 19:00 and 20:00, it may be seen that the energy price in these two hours is less than 0.1 $/kWh.Since the generation cost of MTs is estimated at 0.1 $/kWh, the program prefers to purchase electricity from the upstream network instead of utilizing MTs to decrease operation costs.However, other MTs have outputs with the capacity specifed in Table 4. Te second point is that the MT 5 in MG 1 is charged (its battery) in this state, which is not logical (due to the emergency state).Te reason for charging battery 5 is that the minimum amount of battery 5 is 30 kW, while the initial charge status of battery 5 is 15 kW.Due to the minimum amount of the battery, in this situation, battery 5 reaches its minimum charge (30 kW) by charging 15.789 kW at 19:00.Te rest of the batteries are discharged until the minimum battery charge is reached.No interruptions occurred at 19:00 and 20:00.To reduce operating costs, the load shedding is 913,513 kW at 19:00 and 1012.724kW at 20:00.
Figures 5 and 7 show priorities 2 and 3 of load shedding per bus, respectively.Since priority 1 of load shedding is of the highest importance, and the highest cost of load shedding in the DR program, no priority 1 of load-shedding has been shed at 19:00 and 20:00.Priority 3 of load shedding has the lowest value due to its lower importance.Terefore,  7, the load shedding in these loads has occurred in the highest amount.According to Figure 5, priority 2 of load shedding is higher in diferent buses.In Figures 5 and 7, buses 24 and 25 have the highest loadshedding rates (because these two buses have a large load).Te proposed model is compared with the results of [41] to validate the results.Te execution time of the proposed method is lower than [41] (and all studies that have used the nonlinear model in this layer).Te technique used by [41] and most of the methods in the literature are not practically applicable to large systems.In contrast, the proposed scheme performs any system in the shortest possible time.Te second advantage of the proposed scheme is to reduce the   Table 3: Parameters for calculating the relevant costs [41,45].

Parameters
Value ($/MWh)   In contrast, in the proposed approach, the number of switches is 8. Te number of MGs formed is 4 and 2 for the proposed approach and the study in [41], respectively.Also, the total load shedding in the proposed approach in the selfhealing period is 104.726 less than the study in [41].Terefore, the proposed approach in this study has a considerable advantage over the existing methods in the literature.In order to validate the results of the proposed method, the comparison of the present study with the study in [41] is mentioned in Table 5.

Solving the Two-Objective Problem of Cost and Emission in the Smart Distribution System Using the Epsilon Constraint
Method.In order to solve the problem, the two-objective method has been used.Te frst objective function is considered the main objective function (f 1 ), and the second objective function, i.e. emission, is considered the problem constraint (f 2 ).Te problem is solved based on the equations mentioned earlier.As shown in Figure 8, the amount of emission and the cost are inversely related.As one increases, the other decreases.Tis is because in order to have less emission, it is necessary to use more units that are more expensive and produce less emission.
First, to solve the epsilon-constraint method, the problem is solved as a single objective by considering the cost.In this case, the cost is at its lowest (21,811.972$), and the emission is at its highest (1,790.768kg).Ten, the problem is solved as a single objective by only considering the emission.In this case, the emission is at its lowest value (186.472kg), and the cost is at its highest (189,082.811$).Ten, to generate Pareto solutions, the model is solved as a single-objective problem by considering the cost as a main objective function and emission as a subobjective function   (constraint).Te problem is solved this way, and 10 Pareto solutions are generated, as shown in Figure 9.
Te fuzzy method is used to select the best solution.On the other hand, this problem is solved in three diferent categories.
Category I: the importance factor of the cost objective function is considered to be double the emission one Category II: the importance factor of the emission objective function is considered to be double the cost Category III: the importance factor of the cost objective function is considered equal to the emission one According to the expressed categories and as shown in Figure 10, the best Pareto solution is obtained in each category.Te second Pareto solution is selected as the best one in category I. Te ninth and eighth Pareto solutions are chosen as the best solutions in category II and III, respectively.Figure 8 shows the cost and emission of these Pareto solutions, and Figure 10 shows their membership values.
Figure 9 shows the cost and emission cumulatively for 24 hours.

. Conclusion
Distribution system restoration is one of the steps of selfhealing smart distribution systems.After isolating the fault, restoration has two stages, including restoration and unit commitment.One of the challenges for researchers to solve the restoration problem in large-scale distribution systems is the considerable time to solve the problem.Tis study presents a new linear model for solving the postfault restoration problem, which signifcantly reduces the time and ensures the optimal solution due to mathematical solvers.In the frst layer, the faulty area is transformed into one or multi-islandoperating MGs.Optimal island-operating MGs after the fault is one of the essential advantages of the proposed scheme.
Furthermore, in the second layer, using the new arrangement in the frst step, the problem of unit commitment in a smart distribution system is solved by a linear model.Tis study uses ES and various dispatchable and nondispatchable generations, along with DR tools.Controllable loads are divided into adjustable loads and interruptible and shiftable loads to exploit the DR tool.Using DR schemes in self-healing mode reduces system operating costs.Comparing the presented numerical results with existing methods ensures the excellent performance of the proposed method.Another goal of this paper is to use multiobjective modeling and the epsilon constraint method to minimize the cost and emission of the system.Te results indicated that, in the multiobjective model, we should have a trade-of between cost and emission to have an appropriate solution among Pareto solutions.

Nomenclature
(i) Te voltage is always close to the nominal value.(ii) Te voltage angle diference is negligible.Terefore, sin θ k � θ k and cos θ k � 1 can be considered.
are mentioned as follows: International Transactions on Electrical Energy Systems (i) Te system under consideration is balanced.(ii) Reactive load distribution is omitted.(iii) Te penetration factor of microturbines (MTs) is considered to be one.Tis means that they cannot work in overload conditions.(iv) Tere are energy storage and dispatchable and nondispatchable DG in a smart distribution network.Storage is essential in smart distribution systems because these resources, due to their very high response speed, can overcome the uncertainty of nondispatchable DG resources.Another factor in a smart distribution network is the ability to run DR programs.In this paper, DR is used for the optimal restoration of distribution systems.In this paper, home devices are divided into two categories: uncontrollable devices and controllable devices.Controllable devices are divided into shiftable loads as the frst priority, interruptible loads as the second priority, and adjustable loads as the third priority.
Shiftable loads (frst priority), such as a microwave, can be transfered for use another time.When these devices start working, they should continue to work without stopping until their work cycle is completed.Interruptible loads (second priority), such as a vacuum cleaner, can also be interrupted at any time.Adjustable loads (third priority) such as a refrigerator can be reduced in emergency situations.

Table 4 :
Second layer results from fault between buses 4 and 5.
Figure 10: Variation of the total membership and cost and emission functions versus Pareto optimal solutions.12International Transactions on Electrical Energy Systems pg (i,t) /pg (i,t) : /∆V j,t : Te amount of voltage changes S i : Apparent power V i,t : Bus voltage ϕ i,t /λ i,t : Binary variable for discharge/charge ESS z i,t /y i,t : Binary variable for shutdown/start-up of MTs u i,t : Binary variable to indicate the up or down status of MTs Cs i,t : Binary variable to determine the buses that are omitted PLS4: Load outage λ t : Selling power to the upstream network μ t : Buying power from the upstream network SOC i,t : State of charge of ESS.
i /DT i : Minimum up-time/down-time of MTs