Automatic fault inception and clearance detection using superimposed quantities and RMS change analysis.
In Detego
Accurate fault timing is the foundation of every COMTRADE recording analysis. The fault inception time determines when to measure fault quantities (phasors, impedance, sequence components), while the clearance time determines the total fault duration and relay operating time. Detego automatically detects both events using a two-method approach: superimposed quantities as the primary detector, with RMS change as a fallback.
The detection algorithm processes all current and voltage channels independently for inception detection, taking the earliest inception across all channels. For clearance detection, the algorithm first restricts to current channels only (voltage recovery dynamics produce artifacts), then classifies the fault type at inception + 1 cycle to identify the faulted phases, and uses only faulted-phase current channels for clearance timing. All clearance results are refined using a Stationary Wavelet Transform (SWT) to pinpoint the exact breaker-operation transient.
Fault inception is the instant when the power system transitions from its pre-fault steady state to the faulted condition. Two detection methods are used in priority order.
The superimposed delta quantity isolates the change component by subtracting the signal value from one cycle earlier. Under steady state, this difference is zero. At fault inception, the sudden change in current produces a spike in the delta signal. The reference value is obtained via timestamp-based linear interpolation, ensuring correct operation with multi-rate COMTRADE segments where the sample rate changes mid-recording.
Interpolated Superimposed Delta
Where
The detection algorithm computes the RMS of the delta signal over a half-cycle sliding window, then applies a guard period, noise-adaptive threshold, and persistence check:
Inception Threshold
Where
The dual-component threshold ensures sensitivity on both clean and noisy recordings.
Multi-Rate Support
For de-energized channels (pre-fault RMS below 0.01) or recordings shorter than the guard period, the algorithm falls back to detecting step changes in the RMS magnitude:
The RMS fallback is inherently slower than the superimposed method because the sliding RMS window requires approximately one full cycle to transition from the pre-fault to fault-level value.
Channels with negligible pre-fault signal (below baseline threshold) are automatically routed to the RMS fallback method. This prevents the superimposed delta threshold from collapsing to the noise floor and producing false triggers on inactive channels such as neutral CTs with no ground current or spare VT channels that carry only noise.
The diagram below illustrates both detection methods working on the same fault event. The superimposed delta (red) produces an immediate spike at inception, while the RMS magnitude (blue) takes approximately one cycle to fully transition to the fault-level value.
Fault inception detection using two methods. The superimposed delta (red) spikes within a single sample of the fault, while the RMS magnitude (blue) takes approximately one cycle to fully transition. Superimposed detection is the primary method; RMS is the fallback.
The complete inception detection pipeline processes the raw signal through three stages, each refining the signal to isolate the fault transition from noise and steady-state oscillation:
Fault clearance detection identifies when the circuit breaker has interrupted the fault current and the system has returned to a stable post-fault condition. Two methods are used: superimposed delta as the primary, with RMS transition as a fallback. The result is then refined using wavelet analysis to pinpoint the exact breaker-operation transient.
Clearance detection preferentially uses current channels only. Voltage channels are excluded from clearance aggregation because voltage recovery dynamics (RLC oscillations, tap-changer adjustments, recovery inrush) produce spurious late transient events that corrupt the fault window. If the recording contains no current channels (voltage-only recording), all channels are used as a fallback.
The clearance algorithm uses the same interpolated superimposed delta signal as inception detection. The key insight: the clearance transient produces the last significant delta spike in the recording.
When the superimposed delta approach does not find a clearance (e.g., recordings that end during the fault), the algorithm falls back to RMS-based detection:
After the primary or fallback method produces an approximate clearance time, a Stationary Wavelet Transform (SWT) refinement step pinpoints the exact breaker-operation transient. This refinement is applied to every detected clearance event, not only fallback cases.
The SWT computes level-1 detail coefficients using a Daubechies db4 high-pass filter. Unlike the standard DWT, the SWT does not downsample — each input sample maps 1:1 to an output coefficient. This makes the transform shift-invariant: a one-sample shift in the input produces a one-sample shift in the output. Detail coefficients spike sharply at discontinuities (fault inception, breaker operation) and are near-zero during smooth sinusoidal regions.
The refinement algorithm:
Safety Contract
A minimum fault duration of 10 ms (half a cycle at 50 Hz) is enforced. Events shorter than this threshold are rejected as noise or switching transients rather than true faults. This prevents the algorithm from detecting brief power system events (capacitor switching, load tap changes) as faults.
When per-channel clearance detection fails on all channels (e.g., on recordings where load current dominates fault current), the algorithm re-runs clearance detection anchored at the globally-detected inception time rather than each channel's individual inception.
This handles cases where the fault signature is difficult to see on current channels but voltage channels clearly show the inception. On load-dominated faults, per-channel inception detection on current channels may latch onto the breaker opening instead of the true fault start, leaving no room to find clearance afterward. By anchoring the search at the global inception (typically derived from voltage channels), the algorithm can detect the breaker opening as a large current drop from the correct starting point.
The fallback only fires when the primary path returned no valid clearance on any channel, so it cannot cause regressions on recordings where per-channel detection already works.
Power systems often experience auto-reclose sequences: the breaker trips to clear a fault, waits a dead time, then recloses. If the fault was transient (e.g., a lightning flashover), the reclose succeeds. If the fault is permanent, a second fault occurs and the breaker trips again. This produces multiple fault events in a single COMTRADE recording.
The multi-event detector uses a two-phase architecture: the first event is detected using the proven single-event inception primitive combined with a forward-scanning clearance detector, and subsequent events use delta scanning for inception with the same forward clearance approach.
The forward clearance detector scans half-cycle RMS values forward from the inception point, using a running-peak tracker and multi-phase validation to distinguish real clearances from mid-fault current dips:
After the first event clears, the algorithm scans the superimposed delta for additional inception spikes separated by quiet gaps, then applies the same forward clearance detection for each subsequent event.
Auto-reclose scenario with two fault events. Event 1 inception uses the proven superimposed delta detector; clearance uses forward RMS with running-peak tracking (orange/green dashes). Event 2 inception is found via delta scan after the dead-time quiet gap. Each event is numbered independently.
Three guards prevent false events from noise, ring-down transients, and mid-fault current variations:
When a single fault is detected, annotations are labelled as before: "Inception (IA)" and "Clearance (IA)". When multiple events are detected on the same channel, they are numbered: "Inception #1 (IA)", "Clearance #1 (IA)", "Inception #2 (IA)", "Clearance #2 (IA)", and so on. An event with no clearance (fault persisting to the end of the recording) will only have an inception marker.
When fault classification scans across the fault window, the scan region is tightened to exclude the first cycle after inception (DC offset settling) and the last cycle before clearance (breaker transient noise). Phasor filters produce spurious negative-sequence current during breaker operation (Kasztenny, SEL, 2019), which would cause false classification changes at the edges of the fault window.
A guard ensures the end of the scan window never precedes the start, so the scan degrades gracefully on very short faults.
Short recordings and de-energized channels fall back to single-event detection, returning at most one inception/clearance pair.
AI Analysis Integration
The X/R Ratio computed signal derives the system X/R from the exponential decay of the DC offset in fault current: , where is the DC decay time constant. The DC offset is measured per half-cycle using , and an exponential fit extracts .
Without bounding, the exponential fit would run from the DC peak all the way to the end of the recording. After the fault clears, the current drops to near-zero noise, which inflates the fitted time constant and causes the X/R ratio to rise artificially. This produces misleading values in the post-fault region.
X/R ratio with and without fault window bounding. The unbounded fit (red dashed) includes post-clearance noise, inflating the time constant and causing X/R to rise artificially. The bounded computation (green solid) drops to zero at clearance, preventing misleading values.
The X/R computation now uses multi-event fault detection internally to identify fault windows. For each detected fault event:
All samples outside any fault window output 0, clearly indicating that no meaningful X/R measurement is available in those regions. When multiple fault events exist (auto-reclose), each event gets its own independent X/R calculation.
No Fault Detected
The fault detection algorithm exposes several configurable parameters via the global settings system.
| Parameter | Default | Description |
|---|---|---|
| superimposedThreshold | Configurable | Signal-proportional threshold factor (fraction of pre-fault RMS) |
| guardCycles | Configurable | Number of cycles to skip at recording start (noise floor measurement) |
| persistenceCount | Configurable | Consecutive half-cycle windows required to confirm trigger |
| noiseMultiplier | Configurable | Noise floor scaling factor (~30 dB margin above quantization noise) |
| inceptionThreshold | Configurable | RMS fallback: relative change threshold for inception detection |
| clearanceThreshold | Configurable | Stability threshold for clearance confirmation |
| transitionDropThreshold | Configurable | RMS fallback: fraction of peak below which clearance transition is detected |
| minFaultDurationMs | Configurable | Minimum fault duration in ms (events shorter than this are rejected) |
Accurate fault inception and clearance timing enables several critical protection engineering assessments:
Timing Accuracy