Abstract

Network reconfiguration is an effective approach to reduce the power losses in distribution system. Recent studies have shown that the reconfiguration problem considering load profiles can give a significant improvement on the distribution network performance. This work proposes a novel method to determine the optimal daily configuration based on variable photovoltaic (PV) generation output and the load profile data. A good combination and coordination between these two varying data may give the lowest power loss in the system. Gravitational Search Algorithm (GSA) is applied to determine the optimum tie switches positions for 33-Bus distribution system. GSA based proposed method is also compared with Evolutionary Programming (EP) to examine the effectiveness of GSA algorithm. Obtained results show that the proposed optimal daily configuration method is able to improve the distribution network performance in term of its power loss reduction, number of switching minimization and voltage profile improvement.

1. Introduction

Over the last few decades, the power loss of distribution system has raised attention of researchers to increase the efficiency of the existing network. Different researchers have proposed different approaches for minimization of power system losses [1]. One of the common approaches is through network reconfiguration.

Network reconfiguration is defined as altering the topological structures of the distribution feeders, by changing the position of tie and sectionalizing switches; however, under normal operation, medium voltage distribution networks operate in radial manner [2, 3]. During normal condition, network reconfiguration is used to improve the power quality, maintain the voltage stability, and reduce power loss, while for fault condition it is used for service restoration. Once the source of fault is removed by opening all the switches surrounding the affected portions, the service can be restored within a very short time.

In the past, various algorithms have been proposed for network reconfiguration. Merlin and Back were the pioneers in proposing network reconfiguration as a tool for power loss reduction [4] in 1975. The network reconfiguration problem was solved by closing all tie switches in the system (resulting in mesh network); then the switches are successively opened, one at a time until a new radial network with minimum power losses is obtained. Later on, Shirmohammadi and Hong improved the Merlin method by avoiding the approximation made in [4]. Shirmohammadi theorem starts by closing all tie switches which are then opened by anther based on optimum power flow results [5]. This method was also found efficient in terms of computation time as compared to Merlin method. Reference [2] introduced the “feeder exchange approach” and calculated the loss reduction due to load transfer between two feeders. A year later, Baran and Wu improved the branch exchange method [2] by proposing two approximated power flow equations and solved the problem of loss minimization and loss reduction as an integer programming problem [6]. In the same year, Lin and Chin [7] have considered service restoration, in addition to loss minimization in solving network reconfiguration problem. Zhu et al. performed the network reconfiguration based on modified heuristic solution and experience system operation rules [8]. In a latest research, Mohd Zin et al. in [9] have used the heuristic search algorithm to find the minimum power losses as an objective function, introduced by the author in [6]. The update principle of heuristic search is made by finding the branch having minimal current, till the optimum solution is obtained.

Most of the methods discussed above require exhaustive search mechanism by considering all possible solutions from a predetermined set and subsequently pick the best one. However, such methods are considered as nonefficient in terms of computation time and space requirement. Intelligent, greedy, and nature observed heuristic techniques have also been proposed in the literature and they are commonly known as Artificial Intelligent (AI) methods [10]. Artificial Intelligent methods are a special class of heuristic search methods [10]. Intelligent based optimization methods have also been utilized in finding optimum tie switch combinations in network reconfiguration problem. Authors have used Genetic Algorithms (GA) [11], fuzzy [1214], neural network [15, 16], fuzzy-GA [17], Particle Swarm Optimization (PSO) [18], Matroid Theory [19], Hybrid Evolutionary algorithm [20], Ant Colony Search (ACS) algorithm [21, 22], Harmony Search Algorithm [23], and Bacterial Foraging Algorithm [24]. Das [12] and Savier and Das [25] have proposed fuzzy based multiobjective approach for loss reduction and considered voltage limit and line limits as well in fuzzy set. Penalty based multiobjective-single-fitness approaches have been utilized in [21], to combine power losses with system constraints including voltage limits and line limit. The multiobjective problem is solved using ACS algorithm.

Besides network reconfiguration, Distributed Generation (DG) such as minihydro, photovoltaic, and wind is installed in distribution system to solve the power loss problem due to voltage drop. However, the placement and size of DG are essential as they have a great impact on the power loss [26]. Many researchers have proposed different methods for optimum DG location and size separately [2628] or simultaneously [29, 30]. However, all the aforementioned works have considered a constant load and fixed DG output. In practical, load pattern and DG generation output dynamically change with time. Therefore, some studies have also considered a variable load pattern in their researches. For example, a gradual approaching method is proposed in [31] to solve dynamic reconfiguration. In [32], researchers solved the network reconfiguration problem with daily variable load implementation as well. These studies have proven that the network performance was improved when the time-varying load profiles were taken into consideration.

Apart from that, there are few works on network reconfiguration with the presence of DG considering variable load profile. For example, network reconfiguration with the presence of wind and photovoltaic generation for total daily power loss reduction using improved genetic algorithm is proposed in [33]. Network reconfiguration with photovoltaic generation for power loss reduction considering the DG penetration and three different load levels is presented in [34].

In this work, a novel method is proposed to solve the network reconfiguration problem based on daily PV generation output and variable load profiles. Gravitational Search Algorithm (GSA) proposed in [35] is used in this paper to determine the best daily system configuration. The proposed method can be viewed in two parts. In the first part, hourly optimal configuration considering the PV generation and load profile is carried out using GSA. In the second part, a selective approach is used to determine the best configuration that gives the lowest power loss for the whole day (24 hrs) or optimal daily configuration. The proposed method is implemented on 33-bus distribution system. The proposed method results are also compared with Evolutionary Programming (EP) to show the effectiveness of GSA.

This paper is organized as follows. In Section 2, mathematical model for network reconfiguration is described. The proposed method including application of GSA and selection approach is provided in Section 3. Selection of PV module is presented in Section 4 while load modelling and PV generation are presented in Section 5. In Section 6, analysis of the performance of the proposed method is discussed. Finally, the conclusions are presented in Section 7.

2. Mathematical Model for Distribution Network Reconfiguration Problem

Network reconfiguration is used to minimize the distribution system power loss. The new configuration selected must fulfil the constraints and maintain the radial structure of the distribution system. In this study, the main objective is to gain minimum power loss of the distribution system. Thus, the objective function can be formulated as follows:where is total line in the system. is real active current of line . is resistance of line . is line number. is the topology status of line (1 = close, 0 = open).

is the total hour considered in the time frame and is the current time considered. In this work, (1) is used to analyse the power loss for every hour () of a day with the load data of hour () and PV generation of hour ().

The power system constraints that are also being considered in the study are as follows.

(a) Voltage Bus Constraint. The value of the bus voltage at each node must be within the acceptable limits to maintain the power quality. For example, in this work, the minimum voltage () used is 0.95 p.u. and maximum voltage () is 1.05 p.u.:

(b) Radial Configuration Constraint. The radial configuration of the network must be maintained after the reconfiguration process. In order to ensure this requirement, graph theory is used to determine the radiality of the network. In this work, a MATLAB function is used to check the radiality of a network as follows:where is the distribution system. If the network is radial, equals 1 (true); else it is 0 (false). A radial distribution system must connect all the buses and it contains no loop. In order to maintain the radial structure of the distribution system, a network with loop form, one of the switches from that loop must be opened to retain the radial structure.

There are several methods for radial power flow analysis such as backward-forward method. However, in this work the proposed method uses Newton-Raphson load flow for power flow analysis to calculate the power loss. This load flow is selected to avoid fussy requirements on radial structure and prevent complicated computation [36]. It can handle the DG in the calculation more effectively as compared to other distribution load flow methods. Furthermore, this method provides high convergence with low storage required which can be easily applied with any metaheuristic optimization method.

The real power of bus , , reactive power of bus , , the differences in real power, , the differences in reactive power, , rectangular Newton-Raphson, and power loss can be computed aswhere and are voltage magnitude of bus and bus , respectively, and are magnitude and angel of element in the bus admittance matrix, respectively, and are voltage angel of bus and bus , respectively, and are specified real and reactive power at bus , respectively, and are real and reactive power of bus , respectively, and is line resistance between bus and bus .

3. Proposed Method

The proposed method for finding the best configuration based on daily PV generation and load profiles is presented in this section. The proposed algorithm consists of two main parts, as shown in Figure 1. In the first part, a day is divided into hours, where equals 24 hours in this work. Using GSA and graph theory, optimal configuration is obtained for each hour while the distribution network is maintained to be radial. Hence, a total of effective configurations are obtained for the network. In part two, the optimal daily configuration is obtained with a selective approach.

Part 1: Finding Effective Configuration Using Gravitational Search Algorithm (GSA). Gravitational Search Algorithm (GSA) is first introduced by Rashedi et al. in [35]. GSA is the newest stochastic search algorithm which is based on Newtonian laws and mass interactions. In this algorithm, the population individuals are referred to as masses and their performances are measured by their position masses. Each mass will have four particulars: its position, its inertial mass, its active gravitational mass, and passive gravitational mass. The position of the mass represented a solution while its gravitational and inertial masses are corresponding to the fitness function. According to Newtonian laws, all these objects will attract each other due to the gravity force. Due to this force, all these objects will move towards the object with heavier mass. In other words, heavy masses equaled good solutions and they move slower than the lighter masses that equaled bad solutions. In this way, the exploitation step of the algorithm is guaranteed.

The details of the proposed method to reduce power loss are based on the following steps.

Step 1 (input data and parameters). Input data are inserted in the program. This includes bus data, line data, PV output, base MVA, base voltage, predefined voltage bus range, radial configuration constraint, maximum iteration, and accuracy. The parameter needed to be set in GSA is the number of masses ().

Step 2 (generating initial masses). The initial population is determined by selecting random switches to be opened in the distribution network to form the masses. The number of switches to be opened is . The set of switches to be opened is written as follows:where represents the position of th mass in the th dimension or th component in th mass. , , and are the switches to be opened in th dimension. Only the set of switches to be opened that satisfy all the constraints will be generated as the mass. At the end of this step, all the elements in the masses for network reconfiguration are then converted to whole number by rounding the decimal number to the nearest positive number.

Step 3 (evaluate fitness of each mass). In the fitness function, power loss is calculated in the current mass population using Newton-Raphson load flow analysis.

Step 4 (update gravitational constant, best mass, worst mass, and inertia mass of each iteration). In order to control the searching accuracy, the gravitational constant, , is initialized at the beginning and is reduced with iteration. Hence, gravitational constant is a function of initial value of gravitational constant, , and iteration, iter, as follows:where is a user specified constant, iter is the current iteration, and max_iter is the total number of iterations.
The active gravitational mass, , passive gravitational mass, , and inertial mass of mass , , are computed using fitness evaluation. According to Newton’s law and law of motion, a heavier mass has higher attractions and move more slowly. Hence, in GSA, a heavier mass is represented as good solution and the pattern of movement is represented by the explorations. The inertia mass, , is updated as follows by assuming all masses are equal:where represents the power loss of the mass at iteration, iter, and and represent the strongest and weakest masses with respect to the lowest and highest power losses in current iteration. The mass of current iteration is represented as . For minimization problem, the and are defined as

Step 5 (calculation of the total force in different directions). The gravitational force, , of mass due to mass at current iteration, iter, can be computed as follows:where is the inertial mass of the mass , is the inertial mass of mass , is a small constant, and is the Euclidian distance between and masses specified as follows:

Step 6 (calculation of acceleration and velocity). The acceleration of the mass at current iteration, iter, in th dimension, and is defined as follows:where is the total force that acts on mass of th dimension and is calculated as follows:The random number between interval , is introduced in GSA. In order to compromise between exploration and exploitation in this algorithm, the exploration must fade out and exploitation must fade in with lapse of iterations. In other words, all masses apply force to each other in the beginning and only one mass applies force to others in the end of the algorithm. Based on the aforementioned concept, , a function of iteration which is the set of first masses with the lowest power loss and biggest mass, is introduced to this algorithm. , the initial value of , is set at the beginning and decreased with iterations. Thus, is decreased linearly with iterations as well. In this way, a smaller value of will result in less interaction between the masses by gravitational force. Hence, the movement and computational of power loss are reduced and lead to convergence [37]. The next velocity of a mass is given as follows:

Step 7 (update masses’ position). The next position of a mass can be expressed as follows:

Step 8 (convergence). Steps 3 to 8 are repeated until the stopping criterion is reached.

Step 9 (matrix of feasible configurations). A matrix consists of feasible configuration and its power loss is created. Each configuration consists of switch to be opened which depends on the chosen test system. Configuration is the best configuration with the minimum power loss at hour . For example, configuration 1 is the best configuration at 0100 and so on as illustrated in

Part 2: A Selective Approach for Optimal Daily Configuration. In the steps above, number of configurations will be obtained for total time, . Since most of the PV generations and the load of certain hour in total time, , are approximately the same, most of the obtained configurations in (22) are similar. The similar configurations, , are removed from the matrix. Hence, only the nonsimilar effective configurations are left in the matrix, . The first nonsimilar effective configuration in the matrix represents the set switch 1, and so forth, until the last nonsimilar effective configuration, and set switch . Hence, a matrix consists of all nonsimilar effective configurations and is obtained as illustrated inEach set switch obtained in Part 2 is used to calculate the total power loss for a day using Newton-Raphson load flow analysis. Then, the total daily power losses of the distribution system for all sets of the switch are obtained and the switch set with the minimum total daily power loss is selected as the optimum configuration of the day, as illustrated inwhere optimal daily configuration = min(TotDaily).

4. Selection of PV Module

The PV generation in this study is dependent on the solar irradiance and ambient temperature of Malaysian Meteorological Department for year 2008 at Kuantan site. The capacity factor (CF) is calculated using four PV modules from [38] based on the solar irradiance mentioned as above for a day. As shown in Figure 2, module D has the highest CF (0.2223); therefore, module D is selected for this study. The PV output for 10000 units and each time interval are calculated using (8) to (12). The PV generation is shown in Figure 3:where is cell temperature [°C] during current time, is ambient temperature [°C], is voltage temperature coefficient of the PV module [V/°C], is current temperature of the PV module [A/°C], is nominal operating temperature of cell of the PV module [°C], FF(module) is fill factor of the PV module, is short circuit current of the PV module [A], is open circuit voltage of the PV module [V], is current at maximum power point of the PV module [A], is voltage at maximum power point of the PV module [V], is output generation of PV module during current time, and is average solar irradiance during current time.

5. Load Modelling and PV Generation

The load profiles were received from Tenaga Nasional Berhad (TNB) [39] as shown in Figure 3. This study provides hourly peak load per unit of the daily peak load.

6. Analysis of Proposed Method

The test system for this study consists of 12.66 kV 33-bus radial distribution. The single line diagram of the system is shown in Figure 2. The system consists of three DGs, five labelled tie switches, 32 unlabelled sectionalizing switches, and 33 labelled buses. The line between bus 1 and bus 2 is referred to as switch 1, line between bus 2 and bus 3 is referred to as switch 2, and so forth until switch 32. The five tie switches are referred to as switches 33, 34, 35, 36, and 37 as shown in Figure 4. The line and load data are taken from [6]. For analysis purpose, the DGs are assumed to be installed at buses 6, 18, and 22. However, the locations can be changed based on the actual site data or using other placement methods such as in [29]. The parameters of GSA used in the simulation are set to , , and = sum of all active power in the system and . The simulations are carried out on a computer with Intel Duo Core, 3.0 GHz, and 3.0 GB RAM.

The equals 100 MVA and all the calculations are performed per unit. In the simulation of the network, four scenarios were considered to analyse the capability of the proposed method. They are as follows.

Scenario 1. The system is without reconfiguration and variation of the load and PV generation are considered.

Scenario 2. The system is with optimal hourly reconfiguration and variation of the load and PV generation are considered.

Scenario 3. The system is with optimal daily configuration obtained from the proposed method.

Scenario 4. The system is with optimal daily configuration obtained by considering 24 demand levels and PV generation as whole.

The results of a matrix consist of feasible configurations and their power loss obtained in Part 1 of the study is shown in Table 4. The highlighted configurations in Table 4 represent repeated configurations that will be removed in selection approach in Part 2 to obtain the nonsimilar effective configurations matrix in Part 2. The remaining are the non-imilar effective configurations. The results of nonsimilar effective configurations matrix in Part 2 of the study are shown in Table 5.

Then, for each set of the switch in Table 5, their total daily power loss is obtained by carrying out Newton-Raphson load flow analysis. The results of the total power loss evaluation for all sets switch obtained in Part 2 are shown in Table 1.

From Table 1, the total daily power loss of set switch 19 which is represented with is the lowest among the 20 sets switch. Even set switch 16 with the lowest power loss at 1300 does not yield the lowest total daily power losses. Hence, set switch 19 is chosen as the optimal daily configuration which is the best configuration for the day. The real results of the four scenarios are summarised in Table 2. The proposed method is carried out using EP for these four scenarios as well as for comparison as shown in Table 2. EP is selected for comparison since it is a common and well-established optimization method. The optimal daily configuration obtained in Scenario 4 is 37, 10, 7, 13, and 17 with GSA and 6, 8, 17, 10, and 5 with EP.

From Table 2, the total daily power loss is reduced by 50.36% in Scenario 2, 48.40% in Scenario 3, and 48.02% in Scenario 4 using GSA as optimizing tool and 44.77% in Scenario 2, 41.97% in Scenario 3, and 41.00% in Scenario 4 using EP as optimizing tool compared with Scenario 1 with same load profile and PV generation. The total daily power loss reduction using Scenario 2 is 1.96% (GSA) and 2.80% (EP) more compared with Scenario 3 and 2.34% (GSA) and 3.77% (EP) more compared with Scenario 4. However, no switching is required for Scenarios 3 and 4. For example, 24 times switching are requested in Scenario 2 while no switching is requested in Scenarios 3 and 4. Scenarios 3 and 4 far outweigh Scenario 2 with an acceptable total power loss reduction less than Scenario 2. Apart from that, Scenario 3 gives slightly higher total daily power loss reduction compared with Scenario 4. This is because more constraints need to be considered when 24 demand levels and PV generation are considered as whole. Hence, this limited the solution options.

The DG installation location and number are assumed in this work. Altering DG installation location and number will automatically effect the power loss. The performance of the proposed method is not affected since the algorithm still the same just altering the load data. In other words, more power loss can be reduced with optimal DG installation location and number for a distribution network. An optimal daily configuration still can be obtained with the proposed method. Similarly for different type of DG with their typical output generation.

The performance of Scenario 2 is reduced when changing the hour slots from 1 hour to 2 hours or 4 hours and so forth. In this case, lower total daily power loss is obtained in Scenario 3 compared to Scenario 2 if the hour slot is increased. The total daily power loss for Scenarios 2 and 3 and optimal daily configuration obtained from Scenario 3 for hour slots of one, two, and four are presented in Table 3. It can be noticed that the total daily power loss in Scenario 3 is less than Scenario 2 for two- and four-hour slots. However, these results are the optimal results since the load demand and PV generation are changing with time. Therefore, the analysis in this work is carried out based on hourly basis.

From this analysis, it clearly shows that the proposed method is effective to determine optimal daily configuration based on PV generation and load profiles for a day. Apart from that, GSA with selection approach gives the best solution compared with EP with selection approach. Figure 5 illustrates the optimal configuration obtained from the proposed method.

The voltage profiles for all scenarios are shown in Figure 6. All the scenarios show the similar shape of the voltage profiles with slight change in the magnitude. The proposed method also improves the overall voltage profiles of the distribution system besides generating lower power loss.

From Figure 6, all bus voltages in Scenario 1 are below the 0.95 p.u. voltage constraint in a day. Although all the bus voltages satisfy the minimum voltage constraint only from time 0200 to 1800, the overall voltage profiles are improved clearly for all the buses in Scenarios 2, 3, and 4. The voltage profiles at time 1800 for all the buses for Scenarios 1, 2, 3, and 4 are shown in Figure 7.

From Figure 7, it can be seen that, at time 1800, the voltage profiles for Scenarios 1, 2, 3, and 4 are almost the same for all buses. Scenarios 2, 3, and 4 improve the voltage profiles of all the buses and the voltages are within acceptable range especially for buses 28 to 33. This implies that reconfiguration considering PV generation and load profiles yields desired results of minimizing the power loss, reducing the number of switching and improving the voltage profiles. Meanwhile, Figure 8 shows the ability of the GSA with proposed method in finding the optimum configuration at times 0100, 0800, and 1800 which required less than 10 numbers of iterations to reach the optimal point. It is to note that the large the number of iterations, the longer computing time required as the iteration is counted based on the fixed looping in programming. Thus, in this work, the number of iterations is set to 100.

7. Conclusion

In this work, a novel method to determine the best configuration based on the daily PV generation and load profiles in the distribution network has been successfully proposed. Experimental results indicate that the proposed method is simple and effective to minimize total daily power loss with no switching required with the optimal daily configuration obtained. Hence, same method can be applied for a year. It also improves the voltage profiles. A 33-bus distribution system with three photovoltaic DGs, different sets of PV generation, and load profiles are successfully used to demonstrate the effectiveness of the proposed method. A comparison between GSA and EP with selection approach is carried out and the result proven GSA outperforms EP in terms of producing lower total daily power loss.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors thank the University of Malaya and Malaysian Ministry of Education (MOE) for supporting this work through High Impact Research Grant (HIR-MOHE D000004-16001) and Universiti Teknologi Malaysia (Q.J130000.2723.00K88).