Intelligent decision support algorithm for distribution system restoration

Distribution system is the means of revenue for electric utility. It needs to be restored at the earliest if any feeder or complete system is tripped out due to fault or any other cause. Further, uncertainty of the loads, result in variations in the distribution network’s parameters. Thus, an intelligent algorithm incorporating hybrid fuzzy-grey relation, which can take into account the uncertainties and compare the sequences is discussed to analyse and restore the distribution system. The simulation studies are carried out to show the utility of the method by ranking the restoration plans for a typical distribution system. This algorithm also meets the smart grid requirements in terms of an automated restoration plan for the partial/full blackout of network.

The outage of feeder needs to be isolated to restore the maximum possible consumer's load by restructuring the distribution network meeting the operational constraints and limits of feeder and equipments. In the present paper the criteria for selecting the restoration strategic plans are: 1. Minimizing the operation of switches. 2. Availability of supply to maximum consumer's load, and 3. Meeting the constraints and unconstrained criteria. Tsai et al. (1993), Berdandon et al. (2011), Feltes et al. (2006), Liu and Lin (2012), Chen et al. (2001); Lotfifard et al. (2011), Lim et al. (2006) and Nagata and Sasaki (2002) present the various methods for reconfiguration of networks after blackout or feeder outage. However, a limited work is available based on fuzzy-grey relation ranking method. Further, in most of these research works, constraint related to distribution transformer has not been included. Transformer is an important source of supply and hence we have incorporated the constraint of transformer capacity limit besides the constraints related to distribution network, like feeders and laterals capacity limitations. In order to compare and verify the accuracy of the study, a prototype distribution network of same structure (Chen et al. 2005) has been considered.

Problem definition
The objective of present study is to minimize the number of switching operations for restoration of power supply to all consumer's load and meeting the distribution system network constraints: 1. Let λ 1 (s) be the minimum number of switching operations. Let N sw be the maximum number of switches in the network. X is the position of switch during operartion. If a switch is opened from closed position or vice versa, the value of switch vector will be 1 and if there is no change in the status of switch, its value will be 0. 2. The maximum feeder's loading λ 2 (s) shall be within the rating of current capacity limitation, if λ 2 (s) ≤ maximum (feeder lines loading). Similarly the maximum lateral loading, λ 3 (s) is given by λ 3 (s) ≤ maximum (lateral lines loading). 3. Restoration of maximum loads at faulted area. 4. Minimum switching operations in the restoration plan. 5. The load in feeder, laterals and transformers should be balanced as far as possible, and the overloading of electrical equipment should be avoided. 6. The reconfigured distribution system should be nearly close to the original system, and the radial structure of distribution system should be maintained.

Mathematical model for reconfiguration
In order to design and develop the strategic plans for reconfiguration and implement them effectively, the following model of constraints and equalities are presented:

Number of switching operations
If λ 1 (S) defines the operations of the number of switches, then, here, X i is switch state vector given by [S 1 , S 2 , S 3 …, S NSW ], N sw = The total switches that can be operated in the network under consideration, X i = status of the switch. The conditions for the switch status are: X i = 1, if switch is opened from closed position or vice versa, X i = 0, if status of switch is not changed. Minimum number of switching operations indicates that the system will be more stability.

Maximum loading among backup feeders
The maximum loading, λ 2 (S) among supported feeder is given by Eq. (2): I FD i represents the current over the supported feeder FD i after switching operations. N FD defines the number of supported feeders. To meet the constraints criteria, λ 2 (s) shall be minimised. This objective function will give the most loaded backup feeder and by this we can have the remaining marginal load.

Maximum loading among backup laterals
Like loading criteria for feeders the supported laterals shall also meet the load criteria. This objective function will give the most loaded backup laterals. A lesser value of λ 3 (s) is preferred. λ 3 (s) is the capacity of supported laterals and LAT i is the load current over the laterals after switching operation and N LAT is the number of lateral branches. For techno-economic operation the λ 3 (S), Eq. (3) is desired to be minimized: where λ 3 (S) defines the supported laterals for maximum loading and I LAT i defines current over of the supported lateral LAT i after switching operation. N LAT defines the number of laterals in the distribution network. The load on the laterals should be minimum for the best operating conditions during restoration.

Unbalanced loading of feeders
The feeders as well as laterals shall have the balanced loading of feeders and laterals. It is an important feature for line loss reduction and voltage stability criteria. Thus, the load unbalancing index of feeders and laterals can be computed using Eqs. (4) and (5) respectively. (1) where, LV FD i is percentage load level of feeder FD i and LV FD is percentage refrence load level which is given by Eq. (5) In the above equation I FD i and IR FD i represents the load current and rated load current of feeder. In order to improve the performance of the system the unbalancing loading index shall be as minimised.

Unbalanced loading of laterals
Similarly, the lateral branches unbalance load index λ 5 (s) can be computed using equation: where, LV LAT i is percentage load level of lateral LAT i and LV LAT is percentage reference load level which is given by Eq. (7) as: In the above equation, I LAT i and IR LAT i represents the load current and rated load current of lateral respectively. This objective function is used to determine the degree of unbalanced loading of laterals, therefore, less value of λ 5 (s) is preferred.

Maximum loading among backup transformer
Transformer is the main source of power supply to feeders and laterals. Its maximum loading capacity and unbalanced loading index after the isolation of the fault needs to be computed and checked. These shall be as minimum as possible. The minimization of maximum loading of transformer due to supported feeders and laterals is desirable. Maximum loading of transformer, λ 6 (s) is computed by Eq. (8) as:

Unbalanced loading of transformer
The unbalanced loading index of transformer, λ 7 (s) is given by Eq. (9), where, where, LV TRSi is percentage load level of transformer TRS i and LV TRS is percentage reference load level which is given by Eq. (10) In the above equation I TRS i and IR TRS i represents the load current and rated load current of transformer respectively. This gives the degree of unbalance loading of transformer for the backup and the value of this function should be minimum.

Constraints
To further optimize the switching operation for the reconfiguration of distribution system, the following constraints shall have to be met: 1. Open switch operation have been complemented by closed switch operation. 2. I LATmin < I j < I LATmax 3. I FEEDERmin < I j < I FEEDERmax 4. I TRSmin < I j < I TRSmax

Fuzzy multi criteria evolution
The fuzzy logic method is a mathematical tool to make decision for vague and imprecise information in power system restoration problems (Chang 2008;Farahani et al. 2007;Gomes Flavio 2006;Gonzalez et al. 2012;Savier and Das 2007;Nagata et al. 1995;Pham et al. 2009;Wong and Lai 2000). Using fuzzy logic data base rules, a strategy with lesser switching operation and better load balance is achievable. Here, the linguistic terms like lesser, better etc. convey the vague nature of information. The restoration plan is considered more preferable if it involves fewer switching operations and better load balance. In restoration process, uncertainties arises when the feeder, lateral or transformer current is changing during the restoration process. These uncertainties can be taken into account using the fuzzy logic tool. It is based on rule-base (system operator experience), membership function of variable, and inference decision engine (IF THEN statement). One can consider the membership function of any type (triangular, sigmoid etc) but generally the fuzzy function is selected based on the nature of the problem. In the present study, three level (Low, Moderate and High) triangular fuzzy functions are considered to simplify the calculations during the restoration process, as shown in Fig. 1. "Appendix 2" presents the values of fuzzy membership function for different objective function. The fuzzy function transforms crisp value to fuzzy value which lie in the range {0, 1}. Then using rule-base and fuzzy inference decision procedure, the fuzzy value related to each defined variable/objective function is computed. These output fuzzy values, after inference, are transferred back to crisp values using de-fuzzification methods (either centre of gravity method or centre of area method). In the consequents the fuzzy sets are low, moderate and high, and can be crisply defined as 1, 0.5 and 0 respectively. Further inference is drawn by calculating the real value of objective function. By real value we can get the firing strength in IF-THEN rule and the weighted average. By using Eq. (11), the crisp de-fuzzification value is derived.
where, µ j and y j are firing strengths of anticedents and consequences; N R represents the number of fuzzy rules. The value of f i * represent the fitness degree of objective function 'f i ' for each restoration plan. For example, the number of switching operations performed for the restoration plan 5 is 5. The rules are as follows: The rule strength of R1, R2 and R3 will be 0.5, 0.5 and 0 respectively for the plan 5. After computing the rules, we use Eq. (11) to translate the rule results into real value by weighted average method. The singleton value for good, moderate and bad are 1, 0.5 and 0 respectively. The corresponding de-fuzzification value will be (0.5*1+0.5*0.5+0*0)/ (0.5+0.5) = 0.750.

Grey regression method
Grey relation theory is based upon the concept that available information is incomplete and/or unknown. It is data analysis technique to solve the multicriteria decision making (MCDM). Such problems (MCDM) are difficult to solve using fuzzy logic tools (Wong and Lai 2000;Zhang and Zhengeai 2008;Cheng et al. 1998;Deng 1982;Dong et al. 2003;Huang and Huang 1996;Liu and Forrest 2007;Deng 1989;Song et al. 2002;Tsai et al. 2003;Wong and Lai 1999;Chang and Yeh 2005;Huang et al. 2008;Lin et al. 2008). The coefficient of grey relation (GRC) of x i with respect to x 0 for kth term is given in Chen et al. (2005).
where, x 0 = (x 0 (1), x 0 (2), x 0 (3), …, x 0 (n)), o = 1, 2, 3, …, n, and The GRG between each comparative sequence x i and refrence sequence xo is derived from average value of GRC. The order of relation between comparative and reference sequences is given by Γ 0i. Higher value of Γ 0i means that the comparative sequence is more close to reference sequence than comparative sequence.
In the next stage of the grey analysis, the GRA is used to measure the preference degree for all feasible restoration plans. The various steps for the fuzzy grey approach for ranking the restoration plans and selecting the satisfactory plan is presented in the fuzzy grey relation.

Fuzzy-grey relation
In order to overcome the limitations of decision conditions related to grey relational method and fuzzy logic tool limitations, the two decision making tools with incomplete and vague information are fused together to form a hybrid fuzzy-grey relational tool (Lin et al. 2008;Basu and George 2014;Pereira Junior et al. 2014;Shahsavari et al. 2014;Liu et al. 2015). This tool overcomes the limitations in the two methods, and make the decision making more relevant and effective. Based on the minimum value of λ, the optimized objective function is decided among various alternatives. Optimized λ value is the best alternative switching operation sequence, loading on feeders, laterals and transformer. It provides the optimized solutions for decision making considering the various constraints and equalities. Choice of restoration plan is a type of multi criteria decision making problem which depends on all the objective functions and constraints considered. In this work we have tried to construct a measurement model via grey relational analysis to provide useful information and help system operator to make a right decision on the problem of service restoration. Figure 2 represents the flow chart for the entire restoration algorithm. Starting with the on-off status of the switches, the feasible restoration plans are generated. Objective function values are computed using Eqs. (1)-(10) for the feasible plans. Further fitness degree of the objective function is evaluated using the fuzzy multicriteria evaluation method. The grey regression analysis (GRA) method is used to calculate the preference index of the restoration plans. Based on the grey regression grades (GRG) the plans are ranked according to their preference order. The addition of objective function's minimization of unbalanced loading of transformer gives the stable restoration plans for the considered network. This improves the system reliability and stability, leading to the improved performance of the system.

Numerical application
In order to show the utility of fuzzy-grey relation method for reconfiguration of distribution network, a distributed transformer as a source of supply is added to the network configuration (Chen et al. 2005).
The distribution system of Taiwan Power Company is considered and presented in Fig. 3, which has main feeders YD 28 supplying power to LAT1, LAT2, LAT3, LAT4, LAT5, LAT6, LAT7, LAT8 and LAT9. Each lateral has its supporting lateral (LAT10, LAT11, LAT12, LAT13, LAT14, LAT15, LAT16 and LAT17). These alternating laterals are connected to the main lateral with the help of switches. When a fault occurs in the system the switches operate to restore the out of service area. "Appendix 1" presents the pre-fault load current of feeders. "Appendix 2" presents the maximum capacity of each feeder, laterals, and the loading of transformer. The switch state vector, X comprises of the main switches S W i (i = 1-9) and alternative switches S W j (j = 1-8). Switch Sw9 always remain closed since lateral 9 is not connected with any other supporting lateral.
To maintain the radial structure of the network, the switch open operation should be followed by switch close operation or vice versa. The switch pair for each restoration plan is given in Table 1 against each restoration plan based on the switching operation performed. The lateral, feeder and transformer loading is presented against the switching vector. All load values are in ampere. The maximum number of switching operations can be 8, thus, the total number of possible restoration plans are 2 8 = 256. The rated capacity of feeder, lateral, and transformer are assumed as 450, 100, and 800 A respectively before the fault condition. Table 1 lists all the feasible restoration plans and the load currents on supporting lines after the restoration of supply. The 22 feasible restoration plans are selected from 256 possible restoration plans. The maximum switching  Table 2 gives the values of all the objective functions. The corresponding data shown in Table 3 are obtained by fuzzy multi-criteria evaluation discussed in "Mathematical model for reconfiguration" section. This table gives us the values of reference sequence which is used further to calculate GRGs. the reference sequence selected for analysis is maximum value of fuzzy evaluation of objective functions: X 0 = (0.750, 0.2522, 0.3469,     Table 4. The grey relation grades computed by Chen et al. (2005) are reproduced also in this table. After the GRA steps are completed, the preference ranking of feasible plans are derived with the related analysis and the ranking is presented in Table 5. From the table we can see that the plan 5 is having the highest rank. The plans get their ranking modified because of the addition of transformer objective functions. The loading of all the transformer present in the distribution system is graphically presented in Fig. 4. Additional objective function increases the stability of the system as the unbalancing of the transformer load can also be handled. The new plans are compared with the earlier ones in Fig. 5. Higher the number of objective functions more is the stability of the system as more parameters are considered. The difference of old and new GRG is calculated in Table 4. Further, mean of difference of new and old GRG is calculated, which comes out to be 0.1092. This means that the system performance has increased by 1.09 % with new GRG. The restoration plans are ranked according to the new GRG values. If the current exceeds in a particular element, it could fail and the current would be shunted to other network element which eventually may fail also. Here all equipments are considered so if feeder gets overloaded, it can transfer load to lateral and it can shift the load to the transformer and system gets more stable. Entire restoration algorithm consists of the fuzzy evaluation and fuzzy grey multicriteria given by flow chart. The inclusion of transformer loading parameters in fuzzy analysis makes the system more stable as more parameters are restored. In previous work, feeder and laterals were restored but in this work, transformer is included and restored after the fault which increases the stability and reliability of the system. In this paper, the objective functions are considered so as to optimize the operation of switches and loading of feeders, laterals as well as on transformers and minimization of unbalanced loading index of feeders, laterals and transformer after switching operation. The constraints in the restoration process to be considered are: (a) maintain the radial structure of the network, (b) no overloaded equipment and (c) higher priority customers should always be supplied first. The result shows that inclusion of the transformer current limits has changed the ranking of plans. The system becomes more reliable with  the minimization of unbalanced loading of transformer current as there is no scope for exceeding the limits during restoration plans. Feasible plans which consider more objective functions make the system more reliable. The reliability and stability of the restored network has increased by 1.09 %. The priority customers can be supplied first on the basis of the preference index plans during partial blackout or full blackout. The consumer loads which are not energized may be fed by the supporting feeders in the neighbourhood via on-off switches. The result shows that restoration process is done using minimum number of switching operations. Safety and operability of transformer, laterals and feeders is taken up by maintaining the line currents within the operational limits of power system components. The topology of distribution system is maintained radial before and after implementing the restoration plans. A switch-opened operation is always followed by a switchclose operation, after every switching operation.

Conclusion
Simulation studies for restoration of distribution system are carried out considering multi-objective problem and fuzzy-grey algorithm. Transformer loading has been considered as an additional objective function in the optimization problem. The result shown is more stable and reliable because more are the objective functions, greater is the stability. Various strategies are derived based on rank. The best strategy is the one with the highest rank. The studies are useful for system operator in taking right decision during the restoration process. The rank of restoration plans is given which makes this method effective and very promising. The studies are useful for electric utilities/power distribution company to improve the customer services and revenue returns.

Discussion and future work
This research work has proposed an intelligent restoration algorithm for the distribution system using fuzzy grey combination. The computation implementation is done by adding the unbalance loading of transformer to the Taiwan Power Distribution network. The result shows that the stability of system has improved with the new additional objective functions. The fuzzy multicriteria evaluation gives the optimization values for various objective functions and using grey regression analysis the ranking of plans has been done. This can help in service restoration of priority based customers. The intelligent algorithm is capable of fulfilling the requirements of smart grid such as stability, automation and reliability. The power system automation enables rapid diagnosis and precise solutions to the particular network outages. The proposed algorithm can be useful in making the restoration process automatic with the predecided rank of plans.
Author's contribution All the research work and experimental analysis has been done by first author and layout and formatting of the manuscript has been done by co-authors. All authors read and approved the final manuscript.