Abstract

Traditional distribution networks are transforming into active distribution networks (ADNs) with the advancement of distributed generators (DGs). Break faults that lead to voltage and current imbalances and fluctuations threat the safety of sensitive power electronic equipment in ADNs. However, locating break faults in ADNs remains a challenge under the influence of DG fault polymorphism. This paper proposes a new method to locate break faults by monitoring information of DG current based on the observability of ADNs. The DG equivalent model is established based on the DG output characteristics under the presence of a single-phase break fault. The characteristics of the fault current contributed by DG are analyzed. An identification method of break faults is proposed based on the variations in DG current. Thus, a fault-locating matrix algorithm is proposed combined with current information and location information of DGs. The simulation results demonstrate that the proposed method can locate break faults quickly and accurately. The method is not affected by topology and changes of DG output and load.

1. Introduction

The scale of distributed generator (DG) access to the distribution network is expanding along with the development and increasing application of renewable energies. Apart from solving energy shortage and environmental pollution problems, DGs have a significant role in protecting and controlling distribution networks [1, 2]. DGs have significantly changed the fault characteristics of distribution networks, thereby introducing challenges in fault identification and location. Traditional distribution networks have been transformed into active distribution networks (ADNs) with a high proportion of DGs. ADNs can actively manage DGs, controllable loads, and switches and flexibly adjust the network topology, thereby providing a new means for protection and control.

The existing fault location methods for distribution networks mainly include the impedance method [3], traveling wave method [4], active injection method [5], matrix algorithm [6], and artificial intelligence algorithm [7]. Many scholars have recently examined these methods while considering the influence of DGs. The study in [8] proposed a method that compensates for fault impedance and errors associated with DGs for distance relays yet requires a multipoint simultaneous measurement. The work in [9] proposed a fault location method based on voltage sag measurements that reduces voltage measurement requirements and is proven suitable for various types of DGs. The study in [10] proposed a heuristic algorithm based on bus currents and voltages that employs a pattern search method to estimate the optimal fault point. However, the reliability of such algorithm is affected by the inaccuracy of feeder parameters and the DG equivalent model. On the basis of the block-sparse Bayesian learning method, the study in [11] formulated the fault location problem as the estimation of a block-sparse fault injection current signal that does not require a DG equivalent yet relies on μPMUs. However, most of the existing fault location methods are aimed at short-circuit faults, and only few of these methods consider break faults. Given that the equivalent circuit and electrical characteristics of short-circuit and break faults completely differ, short-circuit fault location methods cannot be used to identify and isolate break faults.

Distribution lines are susceptible to lightning strikes, external forces, electrical effects, etc., resulting in break faults. With the increasingly complex topology of distribution networks, break faults have been frequently reported in recent years. Break faults lead to a phase-deficient operation and generate negative sequence components that affect power quality and threaten personal safety. ADNs contain many sensitive power electronic equipment, which can be easily damaged by severe voltage and current fluctuations caused by break faults [12]. At present, break faults are mainly identified by detecting abnormal voltages in substations and checked by manual line patrolling [13]. However, this method cannot be easily adapted to the modern distribution networks given its time consumption, poor economy, and low reliability.

Only few studies have examined break faults in distribution networks. For instance, the authors of [14, 15] analyzed the voltage and current characteristics of break faults in a small current grounding system. The study in [16] easily realized break fault line selection given that the negative sequence current of a fault line is much larger than that of a nonfault line. Meanwhile, the work in [17] used the decrease in fault phase current and the phase difference of nonfault phase currents to realize break fault line selection. The study in [18] proposed a break fault location method based on the correlation coefficient of the negative sequence voltage and current. The study in [19] used the zero sequence voltage-amplitude difference between the power and load sides to locate the break fault section. The study in [20] realized break fault detection and location by using intelligent algorithms. However, DGs are very sensitive to the terminal voltage. The voltage drops and imbalances caused by break faults can lead to variations in DG output, which inevitably change the fault characteristics of the distribution network. However, the influence of DGs is not considered in existing break fault location methods, and only few studies have examined the DG output characteristics and equivalent modeling methods.

Given their weak characteristics, accurately locating ADN break faults under the randomness and polymorphism of DGs presents a challenge. With the excellent controllability of ADNs [21], this paper proposes a break fault location method for ADNs based on DG monitoring. Unlike traditional methods that take DG fault output as disturbance in fault identification, this paper uses rich DG fault output information to quickly and accurately locate break faults. The rest of the paper is organized as follows: Section 2 analyzes the DG current characteristics of the break fault in a distribution network and establishes an equivalent model. On the basis of the ADN break fault compound sequence network, Section 3 deduces the analytical expressions of the DG current when a break fault occurs upstream and downstream of the DG. Section 4 constructs a fault identification criterion based on the characteristics of the variation in DG fault current. In Section 5, a new break fault location matrix algorithm for ADNs is proposed. Section 6 presents the simulation results of the proposed method. Conclusions are drawn in Section 7.

The main contributions of this paper are as follows: the equivalent model of DG under the break faults in the distribution network is established. The break fault characteristics of ADNs with consideration of the impact of DGs are obtained. A new method for break fault identification based on DG current is proposed. The method is not affected by the changes of ADN topology, load, and output of DG and has high reliability. A new matrix algorithm for break fault location according to DG current is proposed. The fault location can be realized by simple matrix calculation without mathematical complexity and does not rely on external devices and is highly feasible.

2. Equivalent DG Model under the Presence of Break Faults

DGs can be divided into the asynchronous motor interface and the inverter interface based on the ways of grid connection. The asynchronous motor interface DGs mainly include doubly fed wind turbines, whereas the inverter interface DGs mainly include photovoltaic power, direct-drive wind turbine, microturbine, and energy storage. The stator of the doubly fed wind turbine is directly connected to the power grid, and the rotor is connected by a back-to-back converter. Given the limited exchange capacity of the grid-side converter (GSC), the GSC current is much lower than the stator current. Therefore, the GSC current can be ignored [22]. In a normal operation, the doubly fed wind turbine output tracks the reference values and the d- and q-axis components of the stator current can be expressed aswhere and denote the d- and q-axis components of the stator current, respectively, is the terminal voltage amplitude, and and are the reference active and reactive power values, respectively.

The three-phase voltage becomes nonsymmetrical when a break fault occurs. The positive sequence voltage drops, and negative and zero sequence voltages are generated. The drop in the positive sequence voltage is usually small under the presence of a single-phase break fault. By relying on the millisecond-level response and regulation capability of converters, DGs quickly enter the fault steady state [23]. In case of a large drop in the positive sequence voltage, DG protection or control can quickly act to avoid fault impact and then exit within a few tens of milliseconds to restore DG control performance and enter the steady state. To avoid the negative sequence current from affecting the stability, dynamic performance, and safety of power electronic equipment, DGs are generally equipped with negative sequence suppression control [24, 25]. In case the break fault is in a steady state, the negative sequence current of DGs is 0.

The rotor-side converter of the doubly fed wind turbine generally adopts a stator-flux-oriented vector control, which means that the q-axis component of the stator flux is 0. According to the mathematical model of a doubly fed wind turbine [26], the d- and q-axis positive sequence components of the rotor current in the case of a steady-state break fault are formulated aswhere and are the stator and mutual inductances, respectively, is the positive sequence component of the stator flux under a break fault, and is the positive sequence component of the terminal voltage under a break fault. Subscripts + and − denote the positive and negative sequence components, respectively.

Equation (2) is then substituted into the stator voltage and flux linkage equations of the doubly fed wind turbine. The d- and q-axis positive sequence components of the stator current in the break fault steady state are

The inverter interface DG is connected to the power grid by the converter. The positive and negative sequence components of the inverter output are independently decoupled. Under a negative sequence suppression control, the power equation can be expressed as [27]where and are the d- and q-axis components of terminal voltage and and are the d- and q-axis components of the inverter output current, respectively.

The inverter interface DG adopts vector control based on the grid voltage orientation, which means that the q-axis component of the terminal voltage is 0. The current of the inverter interface DG in the break fault steady state can be expressed as

Equations (2)–(5) show that the DG output current is determined by the terminal voltage and the active and reactive power reference values in the break fault steady state. These active and reactive power reference values are usually adjusted according to the preset strategy. Given the fast response of converters, these reference values can be treated as constant. DG can be equivalent to a positive sequence current source controlled by the terminal voltage in the break fault.

3. Analysis of Active Distribution Network Break Faults

The fault phase current is 0, and the nonfault phase voltage is continuous at the fracture under the single-phase break fault. The boundary condition can be set as follows using the symmetrical component method:where , , and denote the positive, negative, and zero sequence currents at the fault point, respectively, and , , and are the positive, negative, and zero sequence voltage differences across the fracture.

When the break fault occurs upstream of DG, by ignoring the influence of line impedance, the composite sequence network can be established as shown in Figure 1(a). denotes the system voltage of ADN, denotes the bus current, and denotes the DG current in the case of upstream break fault, which is computed as the ratio of complex power to terminal voltage , where . According to Kirchhoff’s law, the relationship between voltage and current is expressed aswhere Z2 is the equivalent impedance of the negative sequence network and Z0 is the equivalent impedance of the zero sequence network that can be computed aswhere ZT1 is the DG step-up transformer impedance, ZS1 and ZS2 are the positive and negative sequence system equivalent impedances, ZiM1 and ZiM2 are the positive and negative sequence load impedances upstream of the fault point, ZiN1 and ZiN2 are the positive and negative sequence load impedances downstream of the fault point, is the neutral grounding impedance, CiM is the capacitance upstream of the fault point, and CiN is the capacitance downstream of the fault point.

From (7), can be expressed as

When a break fault occurs downstream of DG, according to (6), the composite sequence network can be established as shown in Figure 1(b). and denote the DG current and terminal voltages in case of a downstream break fault. The relationship between voltage and current is expressed aswhere is the parallel equivalent impedance and can be expressed as

From (11), can be expressed as

Relative to normal operations, the DG current changes in the presence of a break fault. Combining (10) and (13) reveals that the DG break fault current is related to DG power, system impedance, system voltage, and load impedance. When the DG power SDG, system impedance ZS, and voltage are determined, the DG break fault current is mainly affected by load impedance. The DG break fault current changes at different fault locations. According to (10) to (13), decreases as increases, whereas increases along with . Meanwhile, decreases as the distance between the fault location and bus increases. The characteristics of the variation in DG break fault current along with fault location are shown in Figure 2. When the break fault occurs upstream of DG, a closer distance between the fault location and the head of the feeder corresponds to a smaller and a greater DG break fault current. Meanwhile, a closer distance between the fault location and the DG grid-connected point corresponds to a larger and a smaller DG break fault current. When the break fault occurs downstream of DG, a closer distance between the fault location and the DG grid-connected point corresponds to a smaller and DG break fault current. Similarly, a closer distance between the fault location and the end of the feeder corresponds to a larger and DG break fault current.

4. Break Fault Identification Method

According to (1), DG can be equivalent to a voltage-controlled current source in normal operations. In this case, the relationship between voltage and current can be expressed aswhere and are the DG current and terminal voltages in the normal operation, respectively.

From (14), can be expressed aswhere is the total positive sequence load impedance.

Combining (10), (13), and (15) reveals that when a break fault occurs upstream of DG, the variation in DG current relative to that under normal operations is formulated as

When a break fault occurs downstream of DG, the variation in DG current relative to that under normal operations is formulated aswhere, , , and .

Combined with (17), when , the equation can be solved aswhere St is the zero-crossing power.

Combined with (20), when , the equation can be solved aswhere Sr is the power at the extreme point.

When a break fault occurs upstream of DG, the DG current may be higher or lower than that under normal operations. Equation (17) shows that the DG current variation is mainly determined by DG power and load impedance. The DG current increases when but decreases when . Given that and are much smaller than , always holds, that is, the power at the extreme point is less than that at the zero-crossing point. The curve of the DG current variation when a break fault occurs upstream of DG is shown in Figure 3. When , increases along with power. When , decreases along with increasing power and crosses zero at .

Combined with (18), when , this equation can be solved as

Let , which can be solved as

The load impedance of distribution networks is often significantly greater than the system impedance, and the capacitive reactance of feeders is significantly greater than the load impedance. Therefore, equations (12) and (16) can be expressed as

If , the ratio of to can be expressed as

When a break fault occurs downstream of DG, equation (24) shows that a larger k corresponds to a smaller deviation between and , whereas a smaller deviation between and corresponds to a smaller variation in the DG current. Combined with (18), given that , can be obtained. In other words, equations (21) and (22) are always true. Therefore, the DG current variation always increases along with SDG when a break fault occurs downstream of DG as shown in Figure 3. is always smaller than .

When the DG output current increases or decreases beyond the setting value, a break fault occurs upstream of DG. The setting value is determined by the maximum DG current variation in the case of a break fault downstream of DG. Therefore, the fault identification criterion can be expressed aswhere is the DG current variation and is the setting value of DG current variation that can be calculated aswhere is the DG output current at time t, is the DG output current at time t-1, is the reliability factor that ranges between 1.1 and 1.2, is the error margin, is the minimum DG current in the case of a break fault downstream of DG, is the DG output current under normal operations, is the maximum value of , is the minimum value of , is the maximum value of that is equal to the positive sequence load impedance at the end of the feeder, and k is the ratio of negative sequence impedance to positive sequence impedance of load.

When and , one may obtain . In other words, the DG output current remains unchanged under the presence of a break fault downstream of DG. Given the influence of actual measurement errors and to ensure the sensitivity of the break fault identification criterion, the error margin in (27) can be set between 2 A and 5 A.

5. Break Fault Location Algorithm for Active Distribution Networks

The presence of multiple DGs on the feeder is equivalent to having multiple controlled sources connected in parallel upstream or downstream of the fault point. Given that the structure of the composite sequence network is unchanged, the characteristics of the DG break fault current and the break fault identification criteria remain valid. As shown in Figure 4, each DG grid-connected point is defined as a node and is represented by ①, ②, ③, ④, and ⑤. Meanwhile, the section between the node and bus or between two nodes is defined as the area and is denoted by (1), (2), (3), (4), and (5). Therefore, the network topology matrix D can be constructed as

The elements of the network topology matrix D are defined as

The elements of the fault information matrix F are defined as

The fault location matrix P is defined as

Matrix D represents the positional relationship between areas and nodes, whereas matrix F represents the information of nodes downstream of the break fault. Therefore, matrix P can reflect the upstream, downstream, and fault point areas of a break fault. Each row element of matrix P presents the fault information of the corresponding area. If all elements in the row of matrix P are 0, then this row presents information on the upstream area of the break fault. If the row elements contain −1 and 1 and if the sum of elements is 0, then this row presents information on the downstream area of the break fault. If the row elements only contain 0 and 1 and if the sum of these elements is 1, then this row presents information on the break fault area.

As shown in Figure 4, if a break fault occurs in area (2), then the fault information matrix F can be built as

The fault location matrix P can be calculated as

Equation (36) shows that areas (1), (4), and (5) are the upstream areas of the break fault, area (3) is the downstream area of the break fault, and area (2) is the break fault area.

The flow of the ADN break fault location algorithm is described as follows: First, the DGs are numbered, and the network topology matrix D is formed. Second, the output current of each DG is monitored at each time. Third, the setting value and actual variation in the DG current are calculated to check whether the break fault identification criterion is satisfied. Fourth, the DG information of each node is integrated to form the fault information matrix F. Finally, the fault location matrix P is calculated, and the break fault area is determined according to the characteristics of each row of matrix P.

6. Cases

To verify the reliability of the proposed method, a 10 kV ADN model was built in Matlab/Simulink as shown in Topology A in Figure 5. The distribution network has 5 feeders, and the positive sequence parameters are , , and . Meanwhile, the zero sequence parameters are , , and . The neutral point is grounded with a small resistance Rd of 10 Ω. The system equivalent impedance is 0.5 Ω, and the DG step-up transformer impedance is 1 Ω. The loads are kept static, and the parameters of loads and DGs are shown in Table 1.

All DGs and areas in the system are numbered according to the sequence of feeders and the distance between the node and bus. The nodes on feeder 1 are labeled ①, ②, and ③, which constitute areas (1), (2), and (3), respectively. The nodes on feeder 2 are labeled ④ and ⑤, which constitute areas (4) and (5), respectively. The nodes on feeders 3 and 4 are labeled ⑥ and ⑦, respectively, which constitute areas (6) and (7). The nodes on feeder 5 are labeled ⑧ and ⑨, which constitute areas (8) and (9), respectively.

Following the definition of (30), the network topology matrix D is constructed as

By combining 2731 with the parameters shown in Table 1, the setting values of the DG current variation in each node are calculated as shown in Table 2, where Krel is 1.1 and IE is 2 A.

When , a single-phase break fault occurs at 3 km of feeder 1. The variation in the DG output current before and after the fault is presented in Figure 6, which shows that the amplitude of the DG output current at nodes ② and ③ significantly increases, whereas the other nodes remain unchanged after the occurrence of the break fault. The DG output current of each node is recorded in Table 3. The DG current variation before and after the occurrence of a break fault at nodes ① and ④ to ⑨ ranges from approximately 0.1 A to 0.2 A. The DG current variation at nodes ② and ③ is about tens of Amperes. A greater DG power corresponds to a greater variation. The DG current variation of each node is compared with the corresponding setting value according to Tables 2 and 3. Nodes ① and ④ to ⑨ satisfy , whereas nodes ② and ③ satisfy . By combining (25) and (33), the fault information matrix F can be constructed as

According to (34), the fault location matrix P can be calculated as

Given that the elements in rows 1 and 4 to 9 of matrix P are all 0, areas (1), (4), (5), (6), (7), (8), and (9) represent the upstream area of the break fault. The elements in row 3 contain −1 and 1, and the sum of these elements is 0. Therefore, area (3) is selected as the downstream area of the break fault. The elements in row 2 only contain 0 and 1, and the sum of these elements is 1. Therefore, area (2) is selected as the break fault area; that is, a break fault occurs in the section between nodes ② and ③.

The resulting DG current variation before and after the fault in different fault locations is shown in Table 4. As can be seen from the table, when the fault is located 1 or 2 km away from the bus, the DG current variation of nodes ④ to ⑨ is very small and negative, whereas that of nodes ①, ②, and ③ is positive and satisfies the break fault identification criterion. Area (1) can be selected as the break fault area based on the proposed matrix algorithm. Similarly, when the fault location is 3, 4, or 5 km away from the bus, area (2) can be selected as the break fault area. When the fault location is 6, 7, or 8 km away from the bus, area (3) can be selected as the break fault area. Therefore, the proposed method can reliably locate break faults at different locations and areas.

The St corresponding to each node is calculated using (19) and shown in Table 5. The SDG of each node in Table 1 is less than the calculated St. The analysis in Section 3 reveals that regardless of the break fault location, the output current of the DG downstream of the fault point always increases. The simulation results in Table 3 are consistent with the results of the theoretical analysis.

To verify the effectiveness of the proposed method under different DG power values, only the SDG of node ② is changed. The DG current variation is shown in Table 6. When the SDG of node ② is less than the corresponding St (5.6 MW), the DG output current increases after the occurrence of a break fault. When the SDG of node ② is greater than the corresponding St, the DG output current decreases. As the SDG increases, the DG current variation initially increases and then decreases. These simulation results are consistent with those of the theoretical analysis. Obviously, under various DG power values, node ② always satisfies the break fault identification criterion. Therefore, the fault information matrix F and fault location matrix P do not change, the fault area is still area (2), and the fault location is deemed accurate.

In order to verify the scalability of the proposed method, the structure of the ADN model is changed to Topology B as shown in Figure 5. The matrix of Topology B is constructed as

When a break fault occurs in area (3), the DG current variation of nodes ① to ⑨ are −0.1, −0.1, 7.7, 35.7, −0.1, −0.1, −0.2, and −0.1 A, respectively. Thus, the fault information matrix can be constructed as

Then, the fault location matrix can be calculated as

At present, phase voltage rise and phase current sag are usually used to identify break faults in substations. In order to compare with the proposed method, load parameters of feeder 1 in Table 1 are reset. Loads of 1, 1.5, 2, and 0.5 MW are connected to 1, 4, 6, and 10 km of feeder 1, respectively. When a single-phase break fault occurs at 7 km of feeder 1, the phase current variation of feeder 1 is −3.8 A. The phase current is only reduced by 2%. The phase current decreases slightly, which is difficult to distinguish from normal operating conditions such as load switching. However, in the proposed method, the output current variation of each DG is −0.1, −0.1, 2.5, −0.1, −0.1, −0.2, and −0.1 A, which can accurately reflect that the fault is located between node ② and node ③. Compared with the existing method, the advantage of the proposed method is that it is not affected by load size and distribution and has high reliability.

7. Conclusions

This paper proposes a new locating method of break faults in ADNs based on DG monitoring. The paper studies the DG output current variation characteristics when break faults occur upstream or downstream of DG. Based on the matrices constructed by the network topology and DG fault information, a break fault location algorithm is proposed. The simulation verifies that the method is not affected by the changes of DG power, network topology, load, and fault location. The method has a good application effect in the ADN with a high proportion of DG. In addition, the algorithm is simple to calculate, highly feasible, and does not rely on external monitoring devices.

Data Availability

Data sharing not applicable to this article as no datasets were generated or analyzed during the current study.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China under Grant U1866603.