Assessment of Voltage Imbalance Improvement and Power Loss Reduction in Residential Distribution Systems in Taiwan

In this study, the non-dominated sorting genetic algorithm II (NSGA-II) is used to optimize the annual phase arrangement of distribution transformers connected to primary feeders to improve three-phase imbalance and reduce power loss. Based on the data of advanced metering infrastructure (AMI), a quasi-real-time ZIP load model and typical sample distribution systems in Taiwan are constructed. The equivalent circuit models and solution algorithms for typical distribution systems in Taiwan are built using the commercial software package MATLAB. A series of simulations, analyses, comparisons, and explorations is executed. Finally, the quantitative evaluation results for improving the voltage imbalance and reducing the power loss are summarized. For the series of studies, the percentage reductions in (1) total power imbalance TSI, (2) total line loss TLL, (3) average voltage drop AVD, (4) total voltage imbalance factors for zero/negative sequences Td0/Td2, and (5) neutral current of the main transformer ILCO are up to 45.48%, 4.06%, 16.61%, 63.99%, 21.33%, and 88.01%, respectively. The results obtained in this study can be applied for energy saving and can aid the authorities to implement sustainable development policies in Taiwan.


Background
The worldwide energy demand is continuously increasing. Meeting this demand requires careful work from all parts of power generation, power distribution, and user integration. However, many problems still persist in distribution systems, such as energy loss and voltage imbalance [1], where the complexity of loading conditions is among the many causative factors. In Taiwan, the complexity of the load behavior is also a problem. Taiwan Power Company (Taipower), a state-owned electric power industry, has adopted three-phase three-wire (3Φ3w) and single-phase three-wire (1Φ3w) in its power distribution systems to service three-phase and single-phase loads, respectively. Since the loads are simultaneous, currently, three-phase four-wire (3Φ4w) is preferably deployed. However, it may pose a new problem of three-phase current and voltage imbalance in some feeders.
The system operating conditions are regarded as time-based load behaviors that vary. The optimal phase arrangement of distribution transformers can be used to solve the imbalance problem. Although it can improve voltage and current imbalances and reduce power loss in a practical setting, an extended evaluation is needed.

Literature Review
Energy efficiency and saving are the main concerns in distribution systems because power losses caused by voltage imbalance occur unpredictably. Hence, many researchers and engineers have focused on conservation voltage regulation (CVR). In accordance with studies were heavily dependent on the specific time. Therefore, the non-dominated sorting genetic algorithm (NSGA-II) has been widely adopted [19][20][21] due to its flexibility and efficiency. There have been several successful applications of NSGA-II in recent years. For example, in one study [22], NSGA-II was associated with a discrete event simulation for safety system design and multi-objective maintenance problems. The simulation was performed using different encoding methods, and the results showed various accuracies. In [23], NSGA-II was employed with the adaptive neuro-fuzzy interference system (AN-FIS) to solve the optimization problems of short-term electric power demand forecasting because of its superior forecasting accuracy. Moreover, NSGA-II was implemented in [24] to analyze the energy consumption of a deep-sea self-sustaining profile buoy (DSPB). The results show that NSGA-II can lead to good performance accuracy and timeliness. In addition, a comparative study of multi-objective optimization (MOO) was investigated by Li et al. [25]. The performances of MOOs, such as multi-objective particle swarm optimization (MOPSO), multi-objective genetic algorithm (MOGA), and multi-objective differential evolution (MODE), are compared with NSGA-II. NSGA-II is one of the most commonly used MOOs. The results show that NSGA-II has a better average performance than the others. Moreover, NSGA-II also shows good performance in terms of computational time. Gao et al. [26] used an improved NSGA-II algorithm to predict the load power of a hybrid ship under load uncertainty. The results indicate that the improved NSGA-II can solve the optimization problems with a high running speed and good convergence of the solution set. In another study [27], NSGA-II was selected to optimize the rotor tooth profile of a gerotor pump. The results demonstrate that NSGA-II provides a good solution convergence set at a high computational running speed. Furthermore, Xu et al. [28] used NSGA-II to calculate the performance and optimize the configuration of the annular radiator using heat transfer unit simulation. This NSGA-II application indicates that penalty functions and Lagrange multipliers are not necessary in this method, and they offer faster computational times than the traditional Pareto ranking method.

Aim and Contribution
In our previous work [16], a GA-based approach was proposed to optimize the phase arrangement of distribution transformers connected to a primary feeder to improve the three-phase imbalance and reduce power loss. However, the optimal phase arrangement between the distribution transformers and feeder is available for a specific time. Based on the above explanation, the contributions of this study are as follows: (1) An annual phase arrangement of distribution transformers connected to primary feeders for improving the three-phase imbalance and reducing power loss is proposed using NSGA-II. (2) The five objective functions used in this study are (1) total power imbalance, (2) total line loss, (3) average voltage drop, (4) total voltage imbalance factor, and (5) neutral current of the main transformer. (3) Seven typical distribution systems in Taipower were examined for potential CVR solutions on residential dwellings in Taiwan.

Paper Organization
This study is divided into five sections. Section 1 is the introduction, including the background, aim, and contributions. Section 2 explains the multi-objective optimization problems. In Section 3, NSGA-II is presented to improve the voltage imbalance and reduce power loss in residential distribution systems in Taiwan. In Section 4, a problem formulation with objectives and constraints is proposed. The results and discussion are presented in Section 5. Finally, Section 6 concludes the paper.

Multi-Objective Optimization Problems
A multi-objective optimization problem can be defined as, subject to, In practical problems, a clash between the objectives always occurs. Usually, it is difficult to optimize the parameters for all objectives simultaneously. Hence, Pareto optimality was used to address these problems. Pareto-optimal solutions are non-dominated solutions. All non-dominated solutions are merged to form the Pareto-optimal set. The Pareto front is defined as the set of objective vectors of the Pareto-optimal set, as shown in Figure 1. If the following conditions are satisfied, a solution x dominates y (x ≺ y), as explained below: A vector x * is a Pareto-optimal solution if f (x * ) is non-dominated by any f (x) for every x. subject to, = ( 1 , 2 ,⋯ , ) ∈ , = ( 1 , 2 ,⋯ , ) ∈ , In practical problems, a clash between the objectives always occurs ficult to optimize the parameters for all objectives simultaneously. Henc ity was used to address these problems. Pareto-optimal solutions are n lutions. All non-dominated solutions are merged to form the Pareto-op reto front is defined as the set of objective vectors of the Pareto-optima Figure 1. If the following conditions are satisfied, a solution x dominat plained below:

Non-Dominated Sorting Genetic Algorithm II
A single-objective GA-based method is modified to non-dominate each iteration as a population-based method. The GA has the ability to s solutions in the solution space with multiple boundaries for solving no tinuous, and multimodal problems. The mating procedure in the GA is for unexplored non-dominated solutions in the objective space. In gen optimization methods are used to find a single solution for a multi-obj the first type of method, a priori preference information and single-obje techniques are used. A multi-objective problem is a combination of th classical multi-objective optimization methods, one combines multi-o

Non-Dominated Sorting Genetic Algorithm II
A single-objective GA-based method is modified to non-dominated solutions during each iteration as a population-based method. The GA has the ability to search for different solutions in the solution space with multiple boundaries for solving non-convex, discontinuous, and multimodal problems. The mating procedure in the GA is expected to search for unexplored non-dominated solutions in the objective space. In general, two types of optimization methods are used to find a single solution for a multi-objective problem. In the first type of method, a priori preference information and single-objective optimization techniques are used. A multi-objective problem is a combination of these problems. For classical multi-objective optimization methods, one combines multi-objective functions that are transformed into a single objective function. The other method is to solve the Pareto optimal solution set. The Pareto optimal solution is a set of nondominated solutions. For real-world problems, the Pareto optimal solution is better than a weight-sum single-function solution. Decision makers can make choices based on the current situation.
Most multi-objective optimization algorithms focus on searching the Pareto front approaching the true Pareto front, as shown in Figure 2. In previous studies [19,20], NSGA-II required only a single parameter with easy modification and less computational complexity. The solutions can effectively approximate the true Pareto front.

Fast Non-Dominated Sorting Approach
The fast non-dominated sorting approach is a type of Pareto sorting method. In the Pareto sorting method, Pareto dominance is used to evaluate the fitness of each solution. The fast non-dominated sorting approach can be divided into two steps.

1.
Determination of dominance: Using (5), the dominance of all solutions can be determined. n p is the number of solutions that dominate solution p. S p represents a set of solutions dominated by solution p. If solution p dominates solution q, then solution q is added to the S p . Otherwise, n p = n p + 1. After the sorting procedure, the solutions with n p = 0 belong to the first non-dominated front (F 1 ), as shown in Figure 2.

2.
Sort of the non-dominated fronts: For each solution p with n p = 0, each member q of S p is searched, and n q = n q − 1. If any member q with n q = 0, these members belong to the second non-dominated front (F 2 ). The process continues until all the fronts are identified.

Fast Non-Dominated Sorting Approach
The fast non-dominated sorting approach is a type of Par Pareto sorting method, Pareto dominance is used to evaluate t The fast non-dominated sorting approach can be divided into

Diversity-Preserving Approach
Maintaining the diversity of solutions is an important co tive GA. In NSGA-II, crowding distance is used to maintain di nondominated solutions can be spread uniformly in the object In the crowding distance assignment, the population is so jective function value in ascending order. In (6) and (7), the bo maximum and minimum objective function values are regard The absolute normalized difference in the objective function v tions is regarded as the distance (see Figure 3). The crowding individual distances, as shown in (7)

Diversity-Preserving Approach
Maintaining the diversity of solutions is an important consideration in multi-objective GA. In NSGA-II, crowding distance is used to maintain diversity in a population. All nondominated solutions can be spread uniformly in the objective space.
In the crowding distance assignment, the population is sorted according to each objective function value in ascending order. In (6) and (7), the boundary solutions with the maximum and minimum objective function values are regarded as an infinite distance. The absolute normalized difference in the objective function values of two adjacent solutions is regarded as the distance (see Figure 3). The crowding distance is the sum of the individual distances, as shown in (7) below: cd(x) = ∑ k cd k (x).

Elitist-Preserving Approach
In single-objective optimization, the elitism strategy means that the best be able to survive in the next generation. However, non-dominated solution the population. In NSGA-II, the combined population is formed by parent population at the t-th generation is used to create a new populatio lection, crossover, and mutation. After non-dominated sorting, the non-dom ( 1 , 2 ,⋯ , ) is determined. In the elitist-preserving approach, all members selected for the new population +1 . The remaining members of the popula selected from the remaining non-dominated fronts according to their | 1 ∪ 2 ⋯ ∪ +1 | > , all members of 1 , 2 ,⋯ , are chosen for the pop Then, the first ( − | +1 |) members of the set +1 are chosen for the popu

Constraint Handling Approach
In practice, multi-objective problems are accompanied by multiple con constraint-handling approach is used for tournament selection. If solution x dominates solution y, if any of the following conditions are satisfied: Condition 1: Solution x is feasible and solution y is not. Any feasible solution has a better non-dominated rank than any infeas Regardless of the objective function values, the infeasible solutions can be el

Minimum Manhattan Distance
The Manhattan distance (also known as city block distance) between t an n-dimensional space is the sum of the absolute distances in each dimensi tion that minimizes the distance from the normalized ideal vector is the min hattan distance. After normalization, the ideal vector is denoted as,

Elitist-Preserving Approach
In single-objective optimization, the elitism strategy means that the best solution will be able to survive in the next generation. However, non-dominated solutions are elite in the population. In NSGA-II, the combined population R t is formed by P t and Q t . The parent population P t at the t-th generation is used to create a new population Q t via selection, crossover, and mutation. After non-dominated sorting, the non-dominated front (F 1 , F 2 , · · · , F r ) is determined. In the elitist-preserving approach, all members of set F 1 are selected for the new population P t+1 . The remaining members of the population P t+1 are selected from the remaining non-dominated fronts according to their ranking.
F r are chosen for the population P t+1 . Then, the first (N − |P t+1 |) members of the set F r+1 are chosen for the population P t+1 .

Constraint Handling Approach
In practice, multi-objective problems are accompanied by multiple constraints. The constrainthandling approach is used for tournament selection. If solution x constrained dominates solution y, if any of the following conditions are satisfied: Condition 1: Solution x is feasible and solution y is not. Condition 2: Solutions x and y are both infeasible, but solution x has a smaller constraint violation. Condition 3: Solution x and y are feasible and solution x dominates solution y.
Any feasible solution has a better non-dominated rank than any infeasible solution. Regardless of the objective function values, the infeasible solutions can be eliminated.

Minimum Manhattan Distance
The Manhattan distance (also known as city block distance) between two points in an n-dimensional space is the sum of the absolute distances in each dimension. The solution that minimizes the distance from the normalized ideal vector is the minimum Manhattan distance. After normalization, the ideal vector y opt is denoted as, In the non-dominated solution set, the minimum distance sum between the ideal vector y opt and the selected solution is defined as the minimum Manhattan distance, which can be calculated using, Mathematics 2021, 9, 3254 7 of 17 A flowchart of the proposed NSGA-II-based approach is presented in Figure 4.

Genes and Chromosomes
In this study, the control variables (discrete variables) in the NSGA II-based approach are the connection scheme of a distribution transformer tapped off a feeder. A binarycoded chromosome structure was formed, as shown in Figure 5. Six possible connection schemes were mapped to the six genes. Therefore, the length of the chromosome is the total number of distribution transformers in the target distribution system.

Genes and Chromosomes
In this study, the control variables (discrete variables) in the NSGA II-based approach are the connection scheme of a distribution transformer tapped off a feeder. A binary-coded chromosome structure was formed, as shown in Figure 5. Six possible connection schemes were mapped to the six genes. Therefore, the length of the chromosome is the total number of distribution transformers in the target distribution system.
For three-phase loads, six possible connection schemes were adopted, as shown in Figure 6. However, the individual loads in three-phase loads may be different because the loads at the secondary side of the distribution transformer may consist of three-phase loads, double-phase loads (phase-phase loads), and single-phase loads (phase-neutral loads), as shown in Figure 7. The possible connection schemes for various types of load points are listed in Table 1. For three-phase loads, six possible connection schemes were adopted, as shown in Figure 6. However, the individual loads in three-phase loads may be different because the loads at the secondary side of the distribution transformer may consist of three-phase loads, double-phase loads (phase-phase loads), and single-phase loads (phase-neutral loads), as shown in Figure 7. The possible connection schemes for various types of load points are listed in Table 1.

Objective Function
The five objective functions used in this study are (1) total power imbalance, (2) total line loss, (3) average voltage drop, (4) total voltage imbalance factor, and (5) neutral current of the main transformer.

Total Complex Power Imbalance
The imbalance of the three powers is an important indicator for evaluating the balance of power systems. The imbalance between the three powers can be evaluated as follows: For three-phase loads, six possible connection schemes were adopted, as shown in Figure 6. However, the individual loads in three-phase loads may be different because the loads at the secondary side of the distribution transformer may consist of three-phase loads, double-phase loads (phase-phase loads), and single-phase loads (phase-neutral loads), as shown in Figure 7. The possible connection schemes for various types of load points are listed in Table 1.

Objective Function
The five objective functions used in this study are (1) total power imbalance, (2) total line loss, (3) average voltage drop, (4) total voltage imbalance factor, and (5) neutral current of the main transformer.

Total Complex Power Imbalance
The imbalance of the three powers is an important indicator for evaluating the balance of power systems. The imbalance between the three powers can be evaluated as follows: For three-phase loads, six possible connection schemes were adopted, as shown in Figure 6. However, the individual loads in three-phase loads may be different because the loads at the secondary side of the distribution transformer may consist of three-phase loads, double-phase loads (phase-phase loads), and single-phase loads (phase-neutral loads), as shown in Figure 7. The possible connection schemes for various types of load points are listed in Table 1.

Objective Function
The five objective functions used in this study are (1) total power imbalance, (2) total line loss, (3) average voltage drop, (4) total voltage imbalance factor, and (5) neutral current of the main transformer.

Total Complex Power Imbalance
The imbalance of the three powers is an important indicator for evaluating the balance of power systems. The imbalance between the three powers can be evaluated as follows:

Objective Function
The five objective functions used in this study are (1) total power imbalance, (2) total line loss, (3) average voltage drop, (4) total voltage imbalance factor, and (5) neutral current of the main transformer.

Total Complex Power Imbalance
The imbalance of the three powers is an important indicator for evaluating the balance of power systems. The imbalance between the three powers can be evaluated as follows: where j is the j-th feeder segment, t is the t-th time interval, S 3∅ t,j is the average of the three-phase power, S p t,j is the power in the j-th feeder segment, p ∈ α p , and α p = {a, b, c}. If S t,j = 0, the powers in the jth feeder segment are balanced. The total power imbalance in a feeder at time t is evaluated by: where m is the total number of feeder segments in a feeder. The total power imbalance in a feeder during a designated period can be evaluated by, when time T is selected to be in days, weeks, months, or years, it should be used in 96, 672, 2880, or 35,040 time intervals, respectively.

Total Line Loss
From the perspective of power utility, decreasing the power loss of power systems is a vital issue. Therefore, minimizing the total line losses is regarded as an objective function, as follows: where r p j is the resistance of the j-th line segment at the p phase, and I p j is the current in the jth line segment at the p phase.

Average Voltage Drop
Generally, the greater the load balance, the lower the voltage drops in the feeders. To improve the efficiency of CVR, reducing voltage drops and improving voltage differences among the three phases are important objectives. The voltage drop at the t time interval is evaluated by, The average voltage drop is evaluated by, where n is the total number of load points on the feeder of interest, k is the load point k, t is the t-th time interval, V rated is the rated phase voltage, and V p t,k is the magnitude of the p-phase voltage.

Total Voltage Imbalance Factor
The total voltage imbalance factors for zero and negative sequences, labelled Td 0 and Td 2 are defined as follows: where k is the load point k, t is the t-th time interval, V t,k is the zero-sequence voltage, V (1) t,k is the positive-sequence voltage, and V (2) t,k is the negative-sequence voltage.

Neutral Current of Main Transformer
To avoid malfunction of the zero-sequence relay, the neutral current in a wye-connection transformer is restricted to a specified threshold. The neutral current is summarized by the three-phase currents of the transformer, as follows: where I p t is the p phase current at the t-th time interval.

Constraints
The ranges of operating voltages for various loads in the Chinese National Standard (CNS) are shown in Figure 8. In this study, the operating voltages are restricted to within 95-105% of the nominal voltage for hybrid loads. The operating voltages for all the buses are restricted by, , where k is the load point k, t is the t-th time interval, , (0) is the zero , (1) is the positive-sequence voltage, and , (2) is the negative-sequen

Neutral Current of Main Transformer
To avoid malfunction of the zero-sequence relay, the neutral cur nection transformer is restricted to a specified threshold. The neutral rized by the three-phase currents of the transformer, as follows: is the p phase current at the t-th time interval.

Constraints
The ranges of operating voltages for various loads in the Chinese (CNS) are shown in Figure 8

Sample System
In this study, according to the type of load and conductor, the distribution systems in Taiwan are shown in Figure 9. The sample distr be classified into seven cases: (1) rural-underground cable (Feeder A); line (Feeder B); (3) rural-mixed distribution of overhead primaries and erals (Feeder C); (4) urban-underground cable (Feeder D); (5) urban-ov E); (6) urban-mixed distribution of overhead primaries and undergrou F); and (7) total loads of feeding substation (Substation G). One-line di ple distribution systems are shown in Figure 10. Tables 2 and 3 show th transformer data for Feeder A, respectively. For weekdays and wee daily load curves obtained by load surveys for residential, commercial tion demands in the four seasons were used as load patterns. For the o

Sample System
In this study, according to the type of load and conductor, the sample residential distribution systems in Taiwan are shown in Figure 9. The sample distribution system can be classified into seven cases: (1) rural-underground cable (Feeder A); (2) rural-overhead line (Feeder B); (3) rural-mixed distribution of overhead primaries and underground laterals (Feeder C); (4) urban-underground cable (Feeder D); (5) urban-overhead line (Feeder E); (6) urban-mixed distribution of overhead primaries and underground laterals (Feeder F); and (7) total loads of feeding substation (Substation G). One-line diagrams of the sample distribution systems are shown in Figure 10. Tables 2 and 3 show the line segment and transformer data for Feeder A, respectively. For weekdays and weekends, the typical daily load curves obtained by load surveys for residential, commercial, and industrial section demands in the four seasons were used as load patterns. For the other sample cases, the network parameters proposed in [29] were used. The rearrangement of connection phases for a distribution transformer may require several hundred-dollar level fees. Compared with annual operating and maintenance costs, the rearrangement cost is relatively small. Furthermore, the optimal phase arrangement of distribution transformers connected to a primary feeder is an annual routine task for power utility. Not only power losses can be reduced, but also three-phase imbalance can be improved. Table 4 presents the system performance indices before rearranged structures for the typical residential distribution systems in Taiwan. atics 2021, 9, x FOR PEER REVIEW losses can be reduced, but also three-phase imbalance can be imp the system performance indices before rearranged structures for distribution systems in Taiwan.        The above-mentioned typical sample systems in Taipower were used to evaluate the performance of the proposed approach. In addition, the possible improvement benefits of voltage quality and power loss for typical residential distribution systems in Taiwan can be determined. In this section, the parameters for the proposed NSGA-II-based approach are as follows:  The above-mentioned typical sample systems in Taipower were used to evaluate the performance of the proposed approach. In addition, the possible improvement benefits of voltage quality and power loss for typical residential distribution systems in Taiwan can be determined. In this section, the parameters for the proposed NSGA-II-based approach are as follows: 1.
The service voltages for each typical sample distribution system are listed in Table 5.

Performance Study
In our previous work, a weight sum GA-based approach was used to optimize the phase arrangement of distribution transformers connected to primary feeders to improve three-phase imbalance and reduce power loss [16]. The previous work was extended to address the annual phase arrangement of distribution transformers for comparison purposes. The following are the parameters for the weighting-sum GA-based approach: 1.
Furthermore, multi-objective functions are transformed into a single objective function with weighting factors as follows: Subject to System performance indices after carrying out the GA-and NSGA-II-based optimization processes for typical residential distribution systems in Taiwan are presented in  Tables 6 and 7, respectively. Further comparison was performed between the system performance indices after carrying out the GA-and NSGA-II-based optimization processes. As presented in Table 8, the NSGA-II-based optimization process can provide a balance solution in most cases, especially for large-scale distribution networks with an amount of distribution transformers.

Comparison Study
For comparison, the initial chromosome and its corresponding genes were generated randomly. The chromosome consists of the genes of the connection schemes for loadtapped-off points. To evaluate the performance of the proposed NSGA-II-based approach for voltage imbalance improvement, the improvement rate (IR) was calculated using the following formula: where x be f ore and x a f ter are the electrical quantities before and after performing the NSGA-II-based optimization process, respectively. The improvement benefits of voltage quality and power loss for typical residential distribution systems in Taiwan are tabulated in Table 9. The influence and possibility of major factors influencing the improvement benefits for each case are as follows: 1.
Feeder A: In this case, the typical distribution system is operated with low-load density and short-length feeders. This residential distribution system is the radial feeding of a collective housing area. The sample system has only eight distribution transformers, and the connection schemes of the three distribution transformers cannot be adjusted. Therefore, the improvements in voltage quality and power loss are limited by up to 5%.

2.
Feeder B: In this case, the distribution system services low-load density areas with long-length feeders. This typical system feeds electricity demands for agriculture, animal husbandry, and aquaculture. Obviously, TS I , Td 0 and I LCO can be improved considerably in this case. The reduction in I LCO was up to 80%. However, AV D increased slightly because of the load balance effect. In other words, the phase current with a low load increases. Although the load demands of Feeder B are similar to those of Feeder A, the number of distribution transformers of feeder B is 7.5 times that of feeder A. Therefore, the improvement benefits in Feeder B are better than those in Feeder A.

3.
Feeder C: In this case, the typical distribution system feeds the demands for residential, commercial, agricultural, and industrial electric energy. Compared with other benefit indices, the improvement benefits of TL L and AV D are not significant. However, the improvement benefits were the best among all the cases. The larger the scale of the distribution systems, the greater the potential for improvement benefits in voltage quality and power loss. 4.
Feeder D: A high-reliability power service is required in a science and technology park. Usually, the behaviors of demand loads in a science and technology park are high-load density without the on-peak and off-peak periods. In this case, the typical distribution system is operated with a heavy load density and short-length feeders. The sample system has only eight distribution transformers, and the connection schemes of the distribution transformers cannot be adjusted significantly. The improvement benefits for Td 0 and Td 2 are in conflict. By suppressing the zero-sequence component of the distribution systems, Td 0 can be improved and I LCO can be reduced. Therefore, in the proposed multiple-criteria decision making, two objective functions can be improved.
That is, the zero-sequence component of the distribution systems has a high priority for improvement.

5.
Feeder E: In this case, a typical underground primary distribution system in an urban grid was investigated. Similar to Feeder D, the improvement benefits in Td 0 and Td 2 are conflicting in Feeder E. However, the number of distribution transformers of feeder E is greater than that of feeder D. More improvement benefits in voltage quality and power loss can be obtained in feeder E. 6.
Feeder F: In this case, a typical distribution system feeds the demands for residential and commercial electric energy. After evaluating a series of case studies, qualitative results were obtained. Mixed distribution feeders with a large number of transformers and long-distance transmission have great potential for improving TL L and AV D . 7.
Substation G: Substation G can be regarded as the summation of the above-mentioned six feeders. Therefore, as the individual feeders can be improved to balance as much as possible, substation G will achieve the best benefit, and its improvement benefit is less affected by the characteristics of the specific feeders. Table 9. Improvement benefits of system performance indices after carrying out the NSGA-II-based optimization process.

Conclusions
In this study, the problems caused by a three-phase voltage imbalance were first introduced. The principle of the NSGA-II-based approach was then briefly described. With the characteristics of a fast non-dominated sorting approach, a parameterless diversitypreserving approach, an elitist strategy, and a simple and effective constraint-handling approach, the proposed NSGA-II-based approach can be effectively used to evaluate and analyze the improvement benefits of typical residential distribution systems.
The system operation conditions may have a significant effect on the improvement benefit of the load balance strategy. In this study, the system service voltages for all load points satisfy the requirements of the rules in a year. The system service voltages in typical residential systems are set between 1.0 and 1.03 p.u. of nominal voltage. The larger the scale of the distribution systems, the more potential improvement benefits in voltage quality and power loss. Using the proposed method, the total power imbalance, total line loss, average voltage drop, total voltage imbalance factor, and neutral current of the main transformer can be improved considerably.
In this study, a set of typical daily load curves for residential, commercial, and industrial section demands in four seasons was used as load patterns to evaluate the annual phase arrangement of distribution transformers connected to primary feeders. Therefore, typical improvement benefits for typical residential distribution systems in Taiwan can be obtained. The results of this study can be applied for energy saving, and can aid the authorities in implementing sustainable development policies in Taiwan. In future work, the load demands recorded by supervisory control and data acquisition (SCADA) can be further used to evaluate the annual phase arrangement of distribution transformers for a specified residential distribution system.