1. Introduction
Modern societies are highly dependent on the supply of electrical energy for their daily functioning. The per capita energy consumption has steadily grown, due to the increasing of industrialization and world population; so, fossil fuel and conventional energy resources might not be sufficient for meeting humanity’s energy needs in the medium term [
1]; in other words, the traditional electrical network will not be enough to meet the growing energy demand. For overcoming this issue, the massive implementation of renewable energy resources and efficient energy storage systems have been thought as the future of power generation [
2,
3]. Microgrids play a key role in this solution, since they facilitate the integration of renewable energy resources in distribution systems through distributed generation (DG) [
4]; nevertheless, such integration brings new issues regarding their control and operation [
5].
The fact that microgrids may exhibit different operational modes that are associated with the presence of DG units has become one of the main challenges in their protection coordination [
6]. DG units generate bi-directional power flows and variable short-circuit levels that modify the traditional behavior of distribution networks. According to [
7,
8], protection systems in microgrids must face two remarkable challenges: the first one corresponds to their inherent dynamic behavior, due to intermittent loads and generators, and the second one is related to their operating characteristics (grid connected or islanded). These two aspects tangle the problem of microgrid protection coordination; therefore, finding a protection scheme that guarantees speed, selectivity, and reliability has become one of the most challenging tasks in planning microgrid operation [
9].
Directional over-current relays (DOCRs) are widely used in distribution networks (DNs) and they are recognized as the main protection device due to their simplicity and low implementation cost [
10]. Nevertheless, non-radial flexible systems such as microgrids complicate the setting of DOCRs [
11,
12,
13]. Recently, researchers have focused their efforts on the implementation of optimization techniques for reducing the protection times of DOCRs, guaranteeing selectivity and reliability resorting to non-standard characteristics [
14,
15] and machine learning techniques [
16,
17,
18].
In [
14], an alternative approach to conventional protection coordination modeling was proposed. This is achieved by adding a new constraint to the coordination problem that takes into account the Plug Setting Multiplier (PSM), which is the ratio of the fault current that is seen by the relay to the pickup current. The tripping characteristic of over-current relays is limited by the PSM value. The approach of [
14] takes the limitation of conventional IEC tripping characteristics used nowadays in numerical relays into consideration.
The authors in [
10] proposed a setting group-based adaptive protection scheme to take different configurations of a microgrid that result from the connection and disconnection of DG into account. The work proposes different setting groups, limiting its number according to the capabilities of commercial relays; additionally, a new method based on integer linear programming and particle swarm optimization (PSO) is used to specify the adjustments of DOCRs. A similar approach to the one that was proposed in [
10] is implemented in [
16] and [
17] while using unsupervised machine learning techniques. In [
16], the authors exploit the facility to save multiple setting groups within digital DOCRs to adapt the setting of the relays to the current topology of the network. Because the number of setting groups is much lower than that of possible network topologies, a K-means algorithm is used to classify topologies into a number of clusters equal to the number of setting groups. Based on this work, the authors in [
17] proposed a technique for DOCRs coordination while using the Self-Organizing Map (SOM) clustering algorithm.
In [
18], the protection settings and network topology of a microgrid are modified through a rule-based adaptive protection scheme to enhance network reliability. The protection coordination is achieved by means of a machine-learning methodology that is based on a hybrid artificial neural network and a support vector machine model. In [
19], the authors proposed a protection coordination index to measure the impact of DG on DOCRs coordination. The index is determined by a two-phase non-linear programming optimization problem that considers variations of the maximum DG penetration level with changes in the protection coordination time interval. Authors in [
20] proposed a methodology for coordinating DOCRs that considers topology changes due to single outages of lines or generators. The optimal relay settings when considering the N-1 security criterion are achieved through the solution of a mixed integer nonlinear programming problem. In [
21], a multi-objective swarm optimization algorithm combined with a fuzzy decision-making tool is developed for the coordination of DOCRs in microgrids. The proposed approach aims to overcome the limitations of single-objective optimization algorithms that are applied to minimize the operation time of DOCRs coordination. In [
22], the authors proposed an online adaptive coordination protection scheme for DOCRs using intelligent electronic devices and a communication channel to obtain real-time information to update the configuration of the relays.
This paper proposes a new protection coordination approach of DOCRs in microgrids that integrates and complements the characteristics of the models that are presented in [
10,
14,
16,
17]. The proposed coordination model is applicable to adaptive protection schemes that use different settings groups and considers the current capabilities and limitations of DOCRs. The coordination problem is solved by a genetic algorithm, while several machine learning techniques carry out the clustering of operational scenarios. The contributions of this work are twofold: (1) it adopts the constraint that was proposed in [
14] regarding the PSM while taking into account the limitation of conventional IEC tripping characteristics currently used in numerical relays, and (2) it limits the number of clusters to that of the setting groups according to the capabilities of commercially available relays, as proposed in [
10,
16,
17].
These two features are available in normal operation of commercial relays; nonetheless, they are traditionally not considered in the optimization models of DOCRs. It is a fact that the number of relay setting groups is far less than the possible operating scenarios of a microgrid. Therefore, the proposed approach adopts machine learning techniques to intelligently group the operation scenarios according to the number of setting groups available in the relay. Additionally, a comparison of the performance of different machine learning techniques, which includes K-means and SOM techniques that were adopted in [
16,
17] is presented.
Table 1 shows the features of related research works in the field, evidencing the knowlege dgap in the existing literature that is filled by the proposed approach.
It is worth mentioning that these features are available in commercial relays, but are not considered in traditional protection coordination schemes of DOCRs. It is also a fact that the number of relay setting groups is far less than the possible operating scenarios of a microgrid. Therefore, the proposed approach adopts machine learning techniques to intelligently group operation scenarios, according to the number of setting groups available in the relay. Apart from the aforementioned contributions, this paper compares the performance of several Machine Learning techniques, as detailed in the results section. To summarize, a methodology is proposed for the coordination of DOCRs that jointly uses the K-means, self-Organizing Map, and hierarchical techniques to improve the operating times of the coordination scheme while guaranteeing network security.
4. Methodology
Initially, the microgrid under study is modeled and several simulations are carried out in order to characterize its behavior. In this step, information regarding short-circuit currents that are seen by the relays for a set of operative scenarios and faults is obtained. In this case, three-phase faults were considered in the middle of the lines. However, other types of faults and locations can be considered, which may result in different settings of the DOCRs. It is worth mentioning that every system operator determines the type of faults and locations to be evaluated beforehand in order to establish the protection coordination scheme. From the results of the aforementioned simulations, a cluster analysis is performed when considering several operative scenarios using the machine learning algorithms described in
Section 3. Subsequently, the optimal coordination of DOCRs is performed for each cluster that was obtained by the machine learning algorithms. A genetic algorithm (GA) is used in this step. Nonetheless, is worth mentioning that any other metaheuristic technique can be implemented in this step. Finally, while taking the results of the DOCRs coordination and the defined clusters into account, a heuristic is applied for improving the coordination.
Figure 4 depicts the flowchart of the proposed methodology. Each of the steps is described in detail below.
A short-circuit analysis is first performed in all operative scenarios to perform the clustering. In each operative scenario, each fault current that is seen by each relay of the system is taken, according to the established faults. Once it is established, the database with all of the fault currents that are seen by the relays, the clustering is carried out with each of the automatic learning techniques described in the document. In all of the automatic learning methods implemented, clustering is carried out from the fault currents that are seen by the relays in all operating scenarios. In the particular case of the K-means technique, the clustering is carried out, taking, as a starting point, a centroid that essentially depends on the fault currents seen by the relays.
Once the clusters are defined, in order to identify the coordination of each one, the previously established fault currents are considered and the coordination is evaluated while taking into account the fault current seen by each relay in all the scenarios found in the cluster, while using the GA. The above process is repeated for each cluster.
As already mentioned, the GA is in charge of coordinating each cluster based on the fault currents that are seen by each relay in the operating scenarios present in the cluster. In order to carry out this coordination, the initial individuals of the population are defined randomly, while considering the constraints of the decision variables. In this problem, the individuals are the variables that are to be adjusted in each relay of the system. The GA considers the operating times of all relays as the objective function. What is sought in this case is to reduce operating times, ensuring that the constraints that are presented in the formulation are met.
There are different stopping criteria in the steps of the methodology. In the general methodology, the first one is that there are no coordination violations between relays. The second one is that the maximum number of clusters allowed is reached (this is defined according to the characteristics of the commercial relay that is being used, depending on the maximum number of groups of settings that it has enabled). In the case of the GA, a maximum number of iterations was defined as the stopping criterion (2000 iterations for the test system used).
4.1. Simulation Process
The test microgrid is modeled using the Digsilent Power Factory software [
29], where a set of possible operating scenarios are configured. The configuration of scenarios is carried out taking into account combinations of different generation outputs and topological alternatives. Once each operation scenario is configured, short-circuit current calculations are performed for different fault locations and for each operating scenario. The output data of the simulation process are the short-circuit current values that are seen by each relay for each fault.
Figure 5 presents the flowchart of this stage of the methodology.
4.2. Clustering of Operating Scenarios
The proposed methodology seeks to identify different setting groups, limiting its number to the capabilities of commercial relays. The number of relay setting groups is far less than the possible operating scenarios of a microgrid; therefore, an intelligent grouping of scenarios that meets the limitations of commercial relays is necessary. The data regarding short-circuit currents of the set of operating scenarios is the characteristic evaluated to perform the clustering. The number of sets in which the operation scenarios are grouped is defined. Subsequently, the short-circuit currents obtained in Step 2 are used as input data for the machine learning algorithms. The clustering is performed according to the defined number of groups and short circuit results, while using the unsupervised learning technique that is described in
Section 3, as shown in
Figure 6. Each of these methods provides a possible distribution of the groups of operating scenarios.
4.3. Solution of the Protection Coordination Problem
The protection coordination problem is solved by means of a GA, as previously mentioned. It is worth to mention that in this step any other metaheuristic can be implemented; nonetheless, GAs have proven to be effective in addressing the coordination problem of DOCRs as outlined in [
24,
30,
31]. The short-circuit currents obtained in Step 1 and the groups obtained by the machine learning algorithms in Step 2 are used as input data of the GA. The DOCR coordination model presented in Equations (
1) to (
9) is solved for each group of settings provided by the machine learning algorithms (see
Figure 6).
The clustering carried out using the machine learning algorithms gives possible distributions of operative scenarios in the groups. The characteristics of each distribution must be evaluated in order to define which one guarantees the safest operating condition. The GA solves the DOCR coordination problem for each distribution corresponding to the operating scenarios. Once this process is carried out, the performance of each distribution is evaluated, verifying the coordination times between main and back up relays, as well as the operating times for each relay. Finally, if the performance of any of the distributions of the operating scenarios guarantees the enforcement of the coordination constraints, the process is finished; otherwise, it it proceed to step 4.
4.4. Heuristic Adjustment
An adjustment of the distributions is made using one of the techniques summarized in Section if the performance of the distributions of the operating scenarios does not comply with the constraints of the protection coordination problem
3. Each method provides a possible distribution of the operating scenarios. For performing such adjustment, the performance of each of the methods is analyzed.
The heuristic process is as follows: (1) the method with the best performance is selected, (2) the groups that present coordination problems are identified, (3) the group distributions of other methods with great similarity to groups of step 1 are detected; the coordination of such methods must be adequate, (4) the differences between the members of the selected groups in steps 2 and 3 are analyzed, (5) small modifications are made in the group that is selected in step 2, taking some elements form another group, (6) the DOCR coordination problem is solved for the groups of step 5, and (7) it is verified if the coordination performance of the new distributions complies with the selectivity constraint; if it is so, the process is terminated; otherwise, the second best performing method is selected and the procedure is repeated repeated through step 2.
Figure 7 shows an example of one of the proposed heuristic.
5. Tests and Results
The proposed approach was tested with a benchmark IEC microgrid that integrates DG.
Figure 8 depicts the test microgrid and its parameters can be consulted in [
32]. Also, the characteristics of the protection scheme can be found in [
14].
Table 2 presents the 16 operative scenarios (OSs) taken into account for this microgrid. Such OSs are generated by considering different generation and topological conditions. Note that the topology of the microgrid is modified while considering the operational states of switches CB-1 and CB-2; apart from this, the microgrid is able to operate connected or disconnected from the main network (grid-connected and islanded modes).
For the sake of clarity, all of the relays were labeled with numbers ranging from one to 15, preceded by the letter R located, as illustrated in
Figure 8. Standard IEEE 242 recommends that the coordination time between main relay and backup relays
should be equal to or greater than 0.2 s [
33]. In this case, a
of 0.3 s is considered for ensuring a greater margin of safety. A standard inverse IEC characteristic curve was adopted for the operating characteristic of the DOCRs, with parameters A and B of 0.14 and 0.02, respectively. This type of curve was adopted, since it is the most widely used in protection coordination schemes; nonetheless, any other characteristic curve may be used. It was also assumed that the relays of the microgrid are able to configure a maximum of four setting groups. The test system was implemented in DIgSILENT Power Factory and several faults were considered in the lines of the test network.
5.1. Results with No Clustering
The coordination problem was initially solved when considering a single set of parameters for the 16 operative scenarios described in
Table 2. The problem was solved using a GA. Results are summarized in
Table 3 that shows the total operating times and the number of cases where coordination is not guaranteed; that is, where the constraint given by Equation (
2) is not met. It is evident that coordination cannot be guaranteed in all cases since constraint (
2) is not fulfilled 13 times. From this test, it is concluded that it is extremely difficult (if not impossible) to guarantee security when considering a single set of parameters due to the high number of operating alternatives.
5.2. Results with Two Clusters
In this section, the DOCR coordination problem is solved by using the proposed methodology. The initial number of relay setting groups is defined. In this case, it starts with two clusters and the grouping of the OS that is presented in
Table 2 is carried out.
Table 4 details the OS belonging to each cluster, the operation time and number of violations of the coordination problem that was obtained with different methodologies.
In this case, all of the clustering methods provided the same distribution of OS. In cluster 1, there are no coordination problems; that is to say that selectivity is guaranteed; nonetheless, in cluster 2, there are four cases where coordination is not guaranteed. When comparing with the results that are presented in
Table 3, a reduction of the operating time and number of violations in the selectivity constraint is evident. Nonetheless, the adjustment with the heuristic technique does not present any additional improvement in this case. Therefore, it is proceeded to increase the number of clusters.
5.3. Results Considering Three Clusters
In this case, three clusters are considered and the OS presented in
Table 2 are grouped using the described methods. The distribution of the clusters and the obtained results are presented in
Table 5. Note that the Hierarchical single and average methods provide the same clustering, while the other methods result in different distributions of OS. When comparing with
Table 4, a reduction of the violations is evident; however, some cases still do not achieve coordination, since up to four violations are presented in some clusters. In this case, the method that obtained the best performance was K-means, while the ones with the worst performance were the hierarchical single and average methods. The adjustment using the heuristic technique did not present any additional improvement, so it is expected to increase the number of clusters.
5.4. Results Considering Four Clusters
Table 6 presents the results of the methodology considering four clusters. In this case, the hierarchical single and average methods provided the same distribution of OS, while the other methods found different clusters. When comparing with
Table 5, there is an improvement in all of the methods, except for K-means, which presents the worst performance, while the hierarchical Ward presents the best performance.
An adjustment of the clustering was made while using the heuristic that was proposed in
Section 4.4. The main idea with this heuristic is to improve the results of the clustering methods.
Table 7 and
Figure 9 present the adjustment results. It can be seen that the coordination is guaranteed, since constraint (
2) does not present any violation in any of the clusters. In this case, the optimal coordination is finally achieved and, therefore, no further clustering is required.
Table 8 presents the protection coordination settings for the test system. For each relay,
and
in each group are presented; furthermore, backup relay information is also presented.
For illustrative purposes,
Figure 10 depicts the coordination between a main and a backup relay in OS12 for the IEC microgrid when fault in line DL-5 takes place. In this case, the main and back up relays are RP2 and RB4, respectively. Note that the coordination between these two relays is evident.