This invention relates to an apparatts for determination of a ground fault. The invention also relates to a power distribution network arid an associated method.
According to a first aspect of the invention an apparatus for determining a phase associated with a ground fault on a multi-phase power distribution network is provided, the apparatus configured to, based on determination of a spike fault, determine which of a plurality of RMS values of phase-to-phase fault currents is the minimum, said minimum RMS value indicative of a particular phase of the multi-phase power distribution network associated with the ground fault;
a. wherein the apparatus is configured to determine the occurrence of the spike fault based on a plurality of sample sets obtained at a predetermined sample rate and a residual current value, iR, the sample sets comprising, at a pafticular sample time, the phase current of each phase of the multi-phase power distribution network and the residual current value, iR, comprising a sum of said phase currents, said apparatus configured to determine an RMS value of the residual current, IR(n), as follows;
b. where k is an index value of the sum, n designates a recent sample of a sliding window, and N is the sliding window size in terms of the number of samples, wherein the occurrence of the spike fault is determined based on at least IR(n) being above a predetermined residual current threshold (IRset) for less than a predetermined spike time threshold (IRtime); and
c. wherein the apparatus is configured to determine which of the plurality of RMS values of phase-to-phase fault currents Δixy(n) is the minimum based on phase-to-phase currents ixy(n) comprising the difference between respective phase currents of the samples of each sample set, and the phase-to-phase fault currents Δixy(n) comprising the difference between the phase-to-phase currents for each phase at the beginning and end of the sliding window, the apparatus configured to determine the RMS values of the phase-to-phase fault current ΔIxy(n), as follows;
where x and y designate the two phases used to determine each phase-to-phase current.
In one or more examples, the apparatus is configured for use with a three phase power distribution network comprising phases a, b and c and wherein the RMS values of the phase to phase fault currents that is the minimum is indicative of a particular phase of the multi-phase power distribution network on which a ground fault is present in that if the minimum ΔImin comprises;
ΔIbc, then the ground fault is associated with phase a;
ΔIca, then the ground fault is associated with phase b;
ΔIab, then the ground fault is associated with phase c.
In one or more examples, the apparatus is configured to, in determination of the occurrence of the spike fault, determine a rate of change of the RMS value of the residual current and whether or not said rate of change is greater than a rate of change threshold;
a. the occurrence of the spike fault determined based on at least IR(n) being above the predetermined residual current threshold for less than the predetermined spike time threshold and the rate of change being greater than the rate of change threshold at least during part of the time IR(n) is greater than the predetermined residual current threshold.
In one or more examples, the apparatus is configured to determine whether the fault is of a forward fault direction;
a. based on the RMS values of the phase-to-phase fault current ΔIxy(n) for each combination of phases and on an RMS value of the residual fault current ΔIR(n), the residual fault current comprising the difference between the residual current derived from one of the sample sets at the beginning of the sliding window and derived from one of the sample sets at the end of the sliding window and said RMS value of the residual fault current ΔIR(n) comprising;
b. the apparatus configured to identify a minimum of the RMS values of the phase-to-phase fault current ΔImin and a maximum of the RMS values of the phase-to-phase fault current ΔImax, the apparatus configured to determine a forward fault if;
c. ΔImin is less than a first predetermined fraction of the RMS value of the residual fault current; and d. ΔImax is greater than second predetermined fraction of the RMS value of the residual fault current.
In one or more examples, the first predetermined fraction is equal to or less than 0.2 or 0.15 or 0.1. In one or more examples, the first predetermined fraction is between 0.1 and 0.2 or 0.05 and 0.25. In one or more examples, the second predetermined fraction is equal to or greater than 0.8 or 0.85 or 0.9. In one or more examples, the second predetermined fraction is between 0.8 and 0.9 or 0.7 and 0.95.
In one or more examples, the apparatus is configured to raise a fault alarm based on the determination of a phase associated with a ground fault, the raising of a fault alarm dependent on;
a. the occurrence of the spike fault and the identification of a particular phase of the multi-phase power distribution network associated with the ground fault during said spike fault, at least a predetermined, alarm, number of times in a predetermined alarm period.
In one or more examples, the raising of the fault alarm is further dependent on the determination of a fault direction according to claim 3 for each occurrence of the spike fault.
In one or more examples, the sample rate of the sample sets is in the range of 2-3 kHz, such as 2400 Hz.
In one or more examples, the sliding window size corresponds to a single period of the fundamental frequency of the multi-phase distribution network. Thus, if the sampling rate is 2400 Hz and the fundamental frequency is 50 Hz, then N=2400/50=48 samples. Accordingly, the sliding window would encompass N consecutive samples from the most recently acquired sample.
In one or more examples, the predetermined residual current threshold is less than 10% or less than 5% of the rated current. The rated current may be the current applied to the multi-phase power distribution network, such as from a secondary side of a current transformer connected to the power distribution network.
In one or more examples, the predetermined residual current threshold may be less than 100 mA, or less than 50 mA.
In one or more examples, the predetermined spike time threshold is based on one cycle of the fundamental frequency of the multi-phase distribution network. Thus, for example, for a 50 Hz network, the spike time threshold may be 1/50*1000=20 ms or for a 60 Hz network, the spike time threshold may be 1/60*1000=16.67 ms.
According to a second aspect of the invention a power distribution network including the apparatus for determining a phase associated with a ground fault according to the first aspect is provided.
According to a third aspect of the invention a method for determining a phase associated with a ground fault on a multi-phase power distribution network is provided, the method comprising, based on determination of a spike fault, determining which of a plurality of RMS values of phase-to-phase fault currents is the minimum, said minimum RMS value indicative of a particular phase of the multi-phase power distribution network associated with the ground fault;
a. determining the occurrence of the spike fault based on a plurality of sample sets obtained at a predetermined sample rate and a residual current value, iR, the sample sets comprising, at a particular sample time, the phase current of each phase of the multi-phase power distribution network and the residual current value, iR, comprising a sum of said phase currents, said apparatus configured to determine an RMS value of the residual current, IR(n), as follows;
where k is an index value of the sum, n designates a recent sample of a sliding window, and N is the sliding window size in terms of the number of samples, wherein the occurrence of the spike fault is determined based on at least IR(n) being above a predetermined residual current threshold (IRset) for less than a predetermined spike time threshold (IRtime); and
b. determining which of the plurality of RMS values of phase-to-phase fault currents Δixy(n) is the minimum based on phase-to-phase currents ixy(n) comprising the difference between respective phase currents of the samples of each sample set, and the phase-to-phase fault currents Δixy(n) comprising the difference between the phase-to-phase currents for each phase at the beginning and end of the sliding window, the apparatus configured to determine the RMS values of the phase-to-phase fault current ΔIxy(n), as follows;
where x and y designate the two phases used to determine each phase-to-phase current.
According to a fourth aspect a computer program or a computer program product comprising computer program code which, when executed on a processor having memory associated therewith is configured to perform the method of the third aspect is provided.
There now follows, by way of example only, a detailed description of embodiments of the invention with reference to the following figures, in which:
Ground faults may occur in a power distribution network when one or more of the phases inadvertently connect to ground, such as due to damage to cabling that forms the network. The detection of ground faults in power distribution networks is difficult, particularly when the ground fault is intermittent. When the ground fault on an electrical cable of the power distribution network is at an initial stage, the electrical insulation may not be fully broken. This may result in an intermittent ground fault with an extra high fault impedance. This kind of fault is very hard to detect and it is especially difficult to determine the fault direction.
The apparatus and methods of the invention may advantageously provide for detection of ultra-high impedance ground faults that are in the above-described initial stage. The apparatus of the invention may also be advantageous in that determination of such a fault is based only on measurements of currents. As such, voltages measurements are not required (although may be used in some implementations). The apparatus may provide for, without any additional special hardware, fault direction detection and/or intermittent high impedance fault detection, which may be used to predict the future occurrence of a complete ground fault with high fault current. Thus, the apparatus may advantageously provide for fault direction detection using only measured currents and may provide for intermittent ultra-high impedance ground fault detection, as will be described below.
The apparatus 100 includes a spike fault element 102 for determining the occurrence of a spike fault based on a plurality of the sample sets, as will be described in more detail below. A spike fault flag signal 103 (termed DDB_ST) is generated by the spike fault element 102 on detection of a spike fault. The apparatus further comprises a faulty phase detection element 104 configured to determine which of a plurality of RMS values of phase-to-phase fault currents is the minimum, said minimum RMS value indicative of a particular phase of the multi-phase power distribution network associated with the ground fault (the operation of which will be discussed in more detail below). The faulty phase detection element 104 provides for generation of a phase indicator signal comprising one of a phase A fault signal 105 (termed DDB_FTA), phase B fault signal 106 (termed DDB_FTB) or phase C fault signal 107 (termed DDB_FTC) depending on which phase is determined to have the ground fault. The signal 103 and one of the signals 105-107 may provide for the ultimate output of the apparatus 100.
The apparatus 100, in this examples includes several optional components, which may increase its accuracy and/or functionality. in particular, the apparatus 100 includes a direction element 108 configured to determine the fault direction. A fault direction signal 110 (termed DDB_DIR) is provided by the direction element 108 indicative of the determined direction of the fault, which may be indicative of a forward fault direction. The apparatus 100 may also include an alarm element 111 for providing a fault alarm indication signal (collectively 112, 113, 114) to a further apparatus, which may provide for shutting down of the affected part of the power distribution network for repairs/inspection. The alarm element 105 may provide an indication of the phase to which the alarm relates as its fault alarm indication signal and may thus comprise any one of a phase A alarm 112 (termed DDB_Alarm_A), a phase B alarm 113 (termed DDB_Alarm_B) and a phase C alarm 114 (termed DDCB_Alarm_C).
The apparatus 100 is further shown as receiving a plurality of predetermined thresholds or settings at 115. The settings include a predetermined residual current threshold 116, a predetermined alarm time, Tset, 117, a predetermined, alarm, number of times Nset 118, and CT ratio 119. The CT ratio may comprise the ratio of turns of a primary winding to a secondary winding of a current transformer connected to the power distribution network.
In general, the spike fault element 102 is configured to, based on current samples of each phase, discriminate between an intermittent/transient “spike” or a permanent ground fault. The spike fault element 102 uses a root mean squared (RMS) value of a plurality of residual current values covering a recent time window, each residual current value comprising the sum of the individual phase currents at a particular time, and comparison thereof with predetermined thresholds to determine the occurrence of a spike fault. In some examples a comparison of a derivative of the RMS value to a threshold is additionally used.
In general, the direction element 108 is configured to compare the magnitudes of RMS values derived from phase-to-phase current differences against thresholds to determine a fault direction.
In general, the faulty phase detection element 104 is configured to compare the magnitudes of RMS values derived from phase-to-phase fault currents to determine a minimum indicative of the faulty phase.
In general, the alarm element 111 is configured to count the number of the determined fault spikes, dependent on a faulty phase being identified at the time of the fault spike, over a predetermined sliding time window and raise an alarm if the count is greater than a predetermined threshold.
The spike fault element 102 is configured to perform the following steps. The spike fault element is configured to determine, at block 202, a residual current value for each set of phase current samples. The residual current value ik is determined by the sum of the individual phase currents iA, iB, iC at a particular sample time (equation 1);
i
R
=i
A
+i
B
+i
C (1)
where, iR is residual current sample; iA is current sample of phase A; iB is current sample of phase B; and iC is current sample of phase C.
In other examples, the calculation of the residual current value is performed for a subset of the received phase current samples, e.g. at a lower rate than the sample rate of the current samples.
An RMS value of the residual current is then calculated, at block 203, using a plurality of residual current values determined over a recent period defined by a window size N extending from a current time. Equation 2 shows the calculation of the RMS value;
where n is the most recent index and thus represents its associated value at a recent time, e.g. IR(n) may represent the residual current value derived from the most recent sample set at n; N is the sliding window size in terms of the number of samples used to form the RMS value i.e. corresponding to the number samples sets of phase currents, and k represents an index of the sum function.
The window size is, in this example, set to one cycle of fundamental frequency of the multi-phase power distribution network. For example, if the sampling rate of the phase currents is 2400 Hz and the fundamental frequency of the network is 50 Hz, then N=2400/50=48 samples.
Thus, each RMS value of equation 2 may have an overlap of N-1 samples with a subsequent RMS value in terms of the residual current values used. In other examples, the RMS values may be calculate with a different or no overlap.
The determination of the occurrence of a fault spike may be based on a number of parameters. In some examples, the magnitude of the RMS value may be compared to a threshold and the length of time the RMS value exceeds the threshold may be used to determine the occurrence of a fault spike. In this example a further optional parameter is used, namely a derivative of the RMS value, as will be described below.
Each RMS value IR(n) is compared to the predetermined residual current threshold (IRset) at block 204. If the value of IR(n) is greater than IRset then a flag, DDB_IR is raised at 205. IR(n) being greater than IRset may be indicative of an electrical disturbance and the raising of the flag provides an indication that this parameter is satisfied. The ultimate decision on whether a fault is present may require multiple parameters to be satisfied simultaneously.
The predetermined residual current threshold IRset may be determined based on a rated current in a secondary winding of a current transformer connected to the power transmission network. Typically, the rated secondary current of a current transformer is 1A. For example, the predetermined residual current threshold IRset may be less than 10% or less than 5% or about 0.05 of said rated “secondary” current. For example, if said rated current is 1A, then the threshold IRset may be 50 mA (i.e. 5%).
Block 206 provides for determining the length of time that the IR(n) value is greater than IRset. Block 206 thus provides for discrimination between a transient spike and a permanent ground fault. In this example, block 206 acts on flag DDB_IR, although in other examples it may use the output of block 203, the threshold 201 and a timer. Block 206 thus measures the temporal width of the flag signal DDB_IR at block 207 and compares it, at block 208 to a predetermined spike time threshold (IRtime) 209. The predetermined spike time threshold, in this example, is set to one cycle of the fundamental frequency. Thus, if the temporal width of the DDB_IR flag is longer than IRtime=20 ms (for 50 Hz system) or for a 60 Hz system IRtime=1/60*1000=16.67 ms, then it indicates that the fault is not a spike fault.
The block 208 generates a “true” condition when the temporal width exceeds the IRtime threshold and thus the block 206 includes a logic NOT block 210 such that a flag signal 211 is indicative of the DDB_IR flag persisting for less than the IRtime threshold and thus comprising a “spike”.
As mentioned above, a logic AND function based on the flag 205 and flag 211 may provide for generation of the spike fault flag Signal 103. However, this example shows the RMS values from block 203 provided to a derivative block 212, which may determine a derivative with respect to time of two or more RMS values, i.e. a rate of change of the RMS value of the residual current.
Block 213 compares the derivative of the RMS value to a rate of change threshold 214. If the derivative of the RMS value is greater than the rate of change threshold a flag DDB_deltaR 215 is raised. The flag 215 is provided to dwell for a dwell time, by block 216, on the condition of block 213 being satisfied. Thus, the flag 215 will be provided for the dwell time despite the condition at block 213 being satisfied for less than the dwell time. In this example, the dwell time is 10ms. The dwell time may be related to the fundamental frequency, such as less than the fundamental frequency or less than or about 0.5*fundamental frequency.
The rate of change threshold, in this example, is based on the predetermined residual current threshold IRset, although in other embodiments it may be independent thereof. In particular, the rate of change threshold, in this example, comprises 0.5*IRset. Block 220 provides for scaling adjustment of the IRset value to form the rate of change threshold 214 used by the block 213.
The DDB_IR flag 205 is combined with the DDB_deltaR flag 215 by a logic AND block 217 to generate a DDB_FT flag signal 218.
A further logic AND block 221 is provided to combine the DDB_FT flag signal 218 and the flag signal 211 from the transient/permanent ground fault discrimination block 206 to provide the ultimate output of the element 102, spike fault flag signal DDB_ST 103. A delay block 222 may be used to delay the presentation of the DDB_FT flag signal 218 to the block 221 by a delay time. The delay time, in this example, comprises 20 ms, or some other value related to the fundamental frequency or less than the fundaniental period.
It will be appreciated that other parameters may be used, which may generate respective flags, on which a presumption of the existence of a spike fault may be made by element 102.
The fault direction element 108 receives, at 300, the sample set of the current samples of each phase 101a, 101b, 101c. A block 301 is configured to determine phase-to-phase currents ixy(n) comprising the difference between respective phase currents of the samples of each sample set and, subsequently, phase-to-phase fault currents Δixy(n) comprising the difference between the phase-to-phase currents for each phase at the beginning and end of the sliding window.
Thus block 301 is configured to use equation 3 to determine the phase-to-phase currents ixy(n), as follow;
i
ab(n)=ia(n)−ib(n)
i
bc(n)=ib(n)−ic(n)
i
ca(n)=ic(n)−ia(n) (3)
Thus, the difference between the sampled phase currents in each sample set are determined by equation 3 which is then used to derive the phase-to-phase fault currents Δixy(n). The phase-to-phase fault currents comprise the difference between one of the phase-to-phase currents at a first time and the same phase-to-phase current at a second time. The first and second times, in this example, comprise the beginning and end of the sliding window described above wherein the sliding window covers one cycle of rated frequency (in this example). Thus, the phase-to-phase fault currents are determined by equation 4;
where n represents the associated term at the present time and N is the sliding window size. The sliding window size, as mentioned above may comprise one period of the fundamental frequency and thus N is number of samples per cycle of fundamental frequency. If the sampling rate of the phase current is 2400 Hz and the finidamental frequency is 50 Hz, N=2400/50=48.
Equation 4 additionally shows the determination of the residual fault current ΔiR(n) from the different between the residual current at time n and the residual current at time n-N.
The apparatus is configured to determine, at block 302, the RMS values of the phase-to-phase fault current ΔIxy(n), as follows by general equation 5;
where x and y designate the two phases used to determine each phase-to-phase current.
Thus, in this three phase example, the RMS value of the three phase-to-phase fault currents together with RMS value of the residual fault current are calculated by more specific equations 6-9:
Block 303 is configured to determine which of the RMS values of the three phase-to-phase fault currents determined by equations 6-9 is the minimum and which is the maximum. Accordingly, block 303 uses the equations 10 and 11 below;
ΔImin=min{ΔIab, ΔIbc, ΔIca} (10)
ΔImax=max{ΔIab, ΔIbc, ΔIca} (11)
The block 303 provides the values of the minimum and maximum for threshold comparison. The minimum RMS value ΔImin is provided to comparison block 304. The maximum RMS value ΔImax is provided to comparison block 305. The comparison blocks 304, 305 use thresholds based on the RMS value of the residual fault current. Block 306 determines a first predetermined fraction of the RMS value of the residual fault current and block 307 determines a second predetermined fraction of the RMS value of the residual fault current. In this example, the first predetermined fraction comprises 0.1 and the second predetermined fraction comprises 0.9. However it will be appreciated that other thresholds may be used to suit the characteristics of the network or the desired confidence of the determined direction.
A logic AND block 308 combines the output of the comparison blocks 304 and 305. Thus, the combination of the blocks 304-307 provide for determination of whether or not there is a forward fault by comparing the minimum RMS value of the RMS values of the three phase-to-phase fault component current with 0.1 times the RMS value of the residual fault current, and by comparing the maximum value with 0.9 times the RMS value of the residual fault current. That is, if
ΔImin<0.1ΔIR AND ΔImax>0.9ΔIR (12)
then a forward fault is presumed and the fault direction flag signal 110 is raised to show a forward fault has been identified.
Thus, as described above, the block 108 is configured to determine whether the fault is of a forward fault direction based on the RMS values (calculated at block 302) of the phase-to-phase fault current ΔIxy(n) (calculated at block 101) for each combination of phases and on an RMS value of the residual fault current ΔIxy(n) (calculated at block 302), the residual fault current comprising the difference between the residual current of one of the sample sets at the beginning of the sliding window and one of the sample sets at the end of the sliding window (calculated at block 301). Further, the block 108 is configured to identify a minimum of the RMS values of the phase-to-phase fault current ΔImin and a maximum of the RMS values of the phase-to-phase fault current ΔImax at block 303 and determine a forward fault at blocks 304, 305 if; ΔImin is less than a first predetermined fraction of the RMS value of the residual fault current (calculated at block 304); and ΔImax is greater than second predetermined fraction of the RMS value of the residual fault current (calculated at block 305).
Returning to
The faulty phase detection element 104 determines the phase associated with the ground fault from the minimum RMS value as follows;
If ΔImin=ΔIbc, then element 104 determines that the faulty phase is phase A and provides phase A fault signal 105.
If ΔImin=ΔIca, then element 104 determines that the faulty phase is phase B, and provides phase B fault signal 106.
If ΔImin=ΔIab, then element 104 determines that the faulty phase is phase C, and provides phase C fault signal 107.
In some examples, the output of the apparatus 100 may be a signal indicative of the phase that has ground fault based on the output 103 of the spike fault determination element 102 and the output 105, 106, 107 of the faulty phase detection element 104. Thus, output from the elements 102 and 104 at the same time may provide for output of the apparatus 100.
In this example, the apparatus includes the alarm element 111 which receives the output from the faulty phase detection element 104, spike fault element 102 and direction element 108 and provides for further processing thereof before raising an alarm. The alarm may comprise the output of the apparatus 100.
In particular, the apparatus is configured to determine a forward direction spike fault of phase A when flag signals DDB_ST, DDB_DIR and DDB_FTA are received by the block 600. The apparatus is configured to determine a forward direction spike fault of phase B when flag signals DDB_ST, DDB_DIR and DDB_FTB are received by the block 600. Finally, the apparatus is configured to determine a forward direction spike fault of phase C when flag signals DDB_ST, DDB_DIR and DDB_FTC are received by the block 600. These forward direction spike fault signals for each phase may form the output of the apparatus 100.
However, in the present example, an alarm signal is not generated until a predetermined number of fault flag signals are received within a predetermined recent amount of time. A counter block 602 is configured to count the number of fault flag signals 601 received within a recent time period defined by a predetermined alarm period, Tset, shown at 603. The output 604 of the counter thus comprises an integer of the number of tithes the fault flag signal is generated in a sliding time window. A comparison block 605 compares the output 604 with a predetermined “alarm” number of times, Nset, shown at 606. If the output 604 exceeds the Nset threshold then an alarm signal 607 is generated. Given that the arrangement of
The predetermined values Nset arid Tset, received by the apparatus 100 at 117 and 118, may be user defined.
The functionality of the components discussed above in relation to
Number | Date | Country | Kind |
---|---|---|---|
16275030.1 | Feb 2016 | EP | regional |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/EP2017/053322 | 2/14/2017 | WO | 00 |