Pith. sign in

REVIEW 3 major objections 4 minor 21 references

Optimization-based Settingless Algorithm Combining Protection and Fault Identification

T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash

Pith's one-line read A settingless protection algorithm decides whether a power line is healthy or faulted by solving convex optimization problems that fit measured waveforms to both line models, and it locates and classifies the fault in the same step.

desk verdict A genuinely new convex model-selection approach to settingless line protection with strong simulation work, but the fault-identifiability claim is undercut by a real linear-algebra gap and the model-selection rule lacks a complexity penalty. read the letter →

arxiv 1908.03753 v1 pith:A3MSYNUD submitted 2019-08-10 eess.SY cs.SY

classification eess.SYcs.SY
keywords time-domainprotectionsettinglesslinefaultlocationtypeidentificationconvexoptimizationmedium-voltagelinesparameters
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper tries to establish that a medium-voltage line can be protected without any empirically tuned set-points, by treating protection as a model-selection problem: given a short window of synchronized voltage and current measurements from both ends, the algorithm asks which of two line models, healthy or internally faulted, fits the data better. Each fit is obtained by solving a small convex optimization problem, and the faulted model's optimal parameters give fault location, fault type, resistances, and inception time in the same calculation. This would matter because threshold-based relays need re-tuning as grids change and can misoperate with converter-based generation, while a settingless scheme unifies protection with fault analysis. The reported simulations back the claim: correct decisions in all tested normal, external-fault, and internal-fault scenarios, including noise, line-parameter uncertainty, and converter-based sources.

What carries the argument

The load-bearing object is the pair of line models in equations (1)-(4): a healthy lumped-parameter R-L line and the same line with a shunt fault represented by a $3\times 3$ fault-resistance matrix $Z_F$ parameterized by $R_a,R_b,R_c,R_g$ and location $\alpha$. Around these models the paper builds a model-selection procedure: for each hypothesis about the observation window it minimizes the mean squared mismatch of the model equations, which for fixed measurement matrices becomes a convex quadratic program $x^T H x + F^T x + d$ with bound constraints that can be solved to global optimality. This convexity is what lets one calculation both make the trip decision and estimate fault characteristics quickly enough for protection.

What would settle it

Run the algorithm on a large set of fault-free windows generated by the same grid model with Gaussian measurement noise near a 60 dB signal-to-noise ratio and count how often it declares an internal fault; because the faulted model has five free parameters, a nonzero false-trip rate would show that the raw error comparison needs a model-complexity penalty. A sharper test is to take a healthy-line window, deliberately mis-synchronize a few samples at one end, and check whether the faulted model absorbs the glitch and trips.

Watch

Extended reading notes

Core claim

The central claim is that a settingless time-domain unit protection algorithm can decide whether a medium-voltage line is healthy or internally faulted purely by comparing how well two physics-based models explain one window of synchronized measurements. The healthy model is the differential equation of a lumped R-L line; the faulted model is the same line with a fault at relative distance $\alpha$ and resistance parameters $R_a$, $R_b$, $R_c$, $R_g$. The algorithm enumerates $M+2$ hypotheses: normal operation over the whole window, faulted operation over the whole window, and $M$ mixtures in which a suspected fault-inception interval slides across the window. For each hypothesis it computes the smallest mean squared residual of the model equations by solving a convex quadratic program in the fault parameters, then selects the hypothesis with the smallest error as the true state, using a ratio-based tie-break when the healthy case and the first mixture case are numerically close. The paper reports 100% detection of all simulated internal faults, no false trips for normal operation or external faults, accurate fault-location and resistance estimates, correct inception-interval identification, and stable performance with converter-based generation for signal-to-noise ratios down to roughly 60 dB and line-parameter errors up to about 10%.

Load-bearing premise

The algorithm assumes that the model with the smaller average squared fitting error is the true state of the line, even though the faulted model has five extra adjustable parameters and can therefore fit noise and measurement glitches more easily.

Editorial extensions

If this is right

  • Because the same optimization problems yield both the trip decision and the fault parameters, a relay using the algorithm can send a trip signal and a fault report (location, type, resistances, inception interval) to operators at the same time.
  • Removing threshold set-points means the algorithm does not need to be re-tuned when grid conditions or fault characteristics change, which the paper shows by testing across different source types, noise levels, and line-parameter errors.
  • The algorithm's security against external faults and normal operation was maintained in all reported scenarios, so it can replace phasor-based relays that may be unreliable with converter-based distributed generation.
  • Within a 2 ms observation window and with communication-plus-propagation delays below 1 ms, the remaining computation budget of under 1 ms is achievable with dedicated hardware, so the approach is compatible with sub-cycle protection speeds.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • An immediate testable extension would be to add a model-complexity penalty (e.g., an information criterion) to the mean squared error comparison; the paper's Section III-D heuristic $a_1$ vs $a_2$ suggests the authors already saw the need to break ties at low error levels, and a principled penalty could remove the heuristic altogether.
  • The same model-selection machinery could be applied to transmission lines, but only if the lumped R-L model is replaced by a distributed-parameter model, which would make the optimization non-convex; the paper marks this as future work.
  • The sliding inception-interval design implies a natural refinement: once a faulted window is found, run the mixture cases again with a finer grid over the selected interval only, which should improve inception-time estimates without increasing the original computational budget.
  • A practical deployment constraint not explored in the paper is sensitivity to synchronization errors between the two line ends; because equations (1)-(4) assume time-aligned samples, the 100 kHz sampling and optical-transformer requirements could be paired with a synchronization-error test to see how much misalignment the convex fit tolerates.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The manuscript proposes a settingless time-domain unit protection algorithm for medium-voltage lines. The algorithm evaluates M+2 hypotheses for each observation window: a healthy line over the full window, a faulty line over the full window, and M mixed pre-/post-fault windows. Each faulty hypothesis is formulated as a small convex quadratic program in the fault resistances Ra, Rb, Rc, Rg and the fault location α, and the hypothesis with the smallest mean-squared error is selected. The optimized resistances and α are then used for fault type identification and fault location. The authors evaluate the algorithm on a Simulink implementation of a lumped R-L line model under normal operation, external faults, internal faults, converter-based generation, measurement noise, and line-parameter uncertainty, reporting 100% dependability and security over the tested scenarios and accurate fault-parameter estimates for SNR ≥ 60 dB and ±10% line-parameter errors.

Significance. The central idea—recasting protection as a model-selection problem solved by small convex programs with no empirical thresholds—is novel and potentially useful for MV time-domain protection. The authors also provide extensive simulation coverage of generation types and operating conditions, and the computational footprint appears compatible with relay time budgets. However, the paper's as-stated claim that the algorithm merges protection, fault location, and fault type identification is stronger than what the presented evidence supports: the fault-resistance identification problem is rank-deficient for common fault types, the model-selection rule has no complexity penalty, and the validation uses the same lumped R-L model that the algorithm assumes. These issues are load-bearing for the merged-functionality claim and need to be addressed in revision.

major comments (3)
  1. [Section III-C, Eq. (10)] The claim that H is positive definite is not correct for typical fault types. For a single-line-to-ground fault, KCL at the fault node gives I1p+I2p ≈ 0 for the two unfaulted phases, so in the regression matrix implied by W(RF, α) in Eq. (7) the columns corresponding to Rb and Rc are near zero; for phase-to-phase faults the column for the unfaulted phase vanishes. Consequently H is only positive semidefinite, the optimal Rb/Rc are non-unique, and a QP solver will return arbitrary values for those components. This undermines the claimed fault-type identification from RF and makes the resistance-error averages in Tables III and V uninterpretable for components whose true value is infinite or unidentifiable. The protection decision may still be reliable, but the fault-identification function needs to be reformulated, for example by solving separate structured models for each fault type or by reporting only identifiable parameters.
  2. [Section III-D] The case-selection rule compares raw mean-squared errors without penalizing the larger number of adjustable parameters of the faulty model. The healthy model has no fault parameters, whereas each faulty/mixed case has five free parameters (Ra, Rb, Rc, Rg, α); under noise the more flexible faulty model can fit the data spuriously. The heuristic a1 := ∆1/∆3 and a2 := ∆4/∆1 introduced in Section III-D is not justified by any statistical model-complexity criterion. This is the most likely reason that the security results in Section IV-F, Figure 6, degrade sharply below SNR = 60 dB. The authors should introduce a principled penalty (e.g., AIC/BIC or a chi-square test with proper degrees of freedom) or otherwise justify the model comparison; otherwise the 'settingless' claim that no thresholds are needed is not supported.
  3. [Section IV, simulation setup] The validation is circular in an important respect: the protected line in the Simulink model is represented with the same lumped series R-L equations, Eqs. (1)–(4), that the algorithm itself uses, with no shunt capacitance or distributed-parameter effects. The tests therefore demonstrate internal consistency between the algorithm and its own assumed model, not robustness to the dominant modeling errors of real MV lines. Since the final paragraph of the paper already notes that HV lines require more detailed modeling, the claims should be correspondingly limited for MV applications, or the algorithm should be tested on a higher-fidelity EMTP-type line model (e.g., frequency-dependent distributed-parameter) to establish that the model mismatch does not cause misoperation.
minor comments (4)
  1. [Section IV-A] The text states that xmax provides an upper bound on RF, but the simulation setup sets xmax := [∞, ∞, ∞, ∞, 1]^T. Please clarify whether a finite bound is actually used in the reported tests; with truly infinite bounds and rank-deficient H, the QP can have unbounded or arbitrary solutions.
  2. [Section IV-D, Table III] For K1 faults, Table II sets Rb = Rc = ∞, so the reported small mean resistance errors cannot be reproduced without stating how errors are computed for components whose true value is infinite; please specify whether the average is taken only over finite fault resistances.
  3. [Equations (3) and (4)] The current reference directions are not fully specified; in particular, the appearance of (I1+I2) in Eq. (7) presumes a particular sign convention at the fault node. Please state the directions explicitly so that the identifiability discussion is unambiguous.
  4. [References] Reference [1] contains garbled author names ('M. hrstrm', 'L. Sder', 'G. Andersson') and [15] says '72rd' instead of '72nd'; a careful copyedit of the reference list is needed.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the fault-characteristic outputs are the optimized fit variables of an inverse problem, and no load-bearing self-citation was found.

full rationale

The paper's derivation chain is an inverse-problem formulation rather than a prediction derived from assumptions that already contain the result. The healthy and faulted model mismatch matrices S and W are constructed directly from the line differential equations (1)-(4), and the optimization problems (9) and (14) minimize the mismatch between these model equations and the measured voltage and current samples. The outputs (R_F, alpha) are by definition the best-fit fault parameters, so reporting them as the estimated fault characteristics is the intended estimation task, not a prediction that was presupposed. The comparison of the mean squared errors Delta_m across the M+2 cases is a model-selection heuristic; it may be statistically under-penalized because the faulted model has extra adjustable parameters, and the assertion that H in (10) is positive definite is questionable for fault types where some unfaulted-phase columns vanish. However, these are correctness and robustness concerns, not circularity. The simulation validation generates data with the same lumped R-L model used by the algorithm, so it does not independently validate the model form, but this is a standard in-model-class test rather than a circular derivation. There are no load-bearing self-citations: references [11] and [12] are prior external work on setting-less protection, not the authors' own, and no uniqueness theorem or ansatz is imported from the authors' previous papers. Accordingly, no circular step can be exhibited, and the appropriate circularity score is 0.

Assumptions & free parameters 5 free parameters · 4 assumptions · 0 invented entities

The algorithm depends on several user-defined parameters (N, M, l, sampling rate, Rmax_F) and on modeling assumptions about the line, measurements, and model comparison. No new physical entities are introduced. The most fragile assumption is the unpenalized MSE comparison, which can bias the selection toward the faulted model in noisy conditions.

free parameters (5)
  • N (observation window length) = 200 samples at 100 kHz (2 ms)
    User-defined; controls trade-off between speed and accuracy. Selection of optimal N is stated as outside the paper's scope.
  • M (number of mixture cases) = 10
    User-defined; controls granularity of fault inception timing. Larger M increases accuracy but computational burden.
  • l (samples for derivative estimation) = 5
    User-defined; polynomial fit of second degree over l adjacent samples. Affects derivative quality and removal of samples near inception.
  • Sampling rate = 100 kHz
    Chosen to enable accurate derivative estimation within a few millisecond window; a hardware requirement, not an algorithm input.
  • Rmax_F (upper bound on fault resistances) = infinity (unbounded in simulation)
    Introduced in optimization constraints (9b) but set to infinity in the simulation, so the bound is not active.
assumptions (4)
  • domain assumption The protected line is accurately represented by a lumped R-L circuit without shunt capacitance or frequency-dependent parameters (equations 1-4).
    Invoked in Section II. The conclusion concedes that high-voltage transmission lines require more detailed modeling, but medium-voltage lines are claimed to fit this model without validation against more detailed line models.
  • domain assumption Measurements from both line ends are synchronized, unbiased, and have low noise, and current derivatives can be reliably estimated from sampled data.
    Section III-E lists hardware requirements: optical transformers, high sampling frequency, and low-latency communication. The algorithm's performance under realistic noise levels (SNR >= 60 dB) depends on this assumption.
  • ad hoc to paper The residual comparison across models with different numbers of parameters is a valid model-selection criterion without a complexity penalty.
    Section III-D chooses the model with the smallest mean squared error. The faulted model has five extra free parameters and can therefore fit noise more easily, yet no statistical penalty or significance test is applied.
  • domain assumption The fault model of equations (3)-(4) with resistance matrix (5) covers all internal faults of interest, including high-impedance faults.
    Section II defines the fault model. The simulation tests only the four fault types listed in Table II, all of which are consistent with this model structure.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Optimization-based Settingless Algorithm Combining Protection and Fault Identification." pith.science (2026). https://pith.science/paper/A3MSYNUD

@misc{pith2026190803753,
  author       = {Pith},
  title        = {Pith review of: Optimization-based Settingless Algorithm Combining Protection and Fault Identification},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/A3MSYNUD}},
  note         = {Machine review of arXiv:1908.03753}
}
read the original abstract

The demand for faster protection algorithms is growing due to the increasingly faster dynamics in the system. The majority of existing algorithms require empirically selected set-points, which may reduce sensitivity to internal faults and cause security problems. This paper addresses these challenges by proposing a settingless time-domain unit protection algorithm for medium-voltage lines. The main idea of the algorithm is to identify which model of a protected line, i.e. healthy or with an internal fault, is more consistent with the input measurements. This is done by solving a number of small-scale convex optimization problems, which at the same time determine the characteristics of an internal fault that best fit the measurements. Thus, the proposed algorithm merges protection, fault location and fault type identification functionalities. The algorithm's performance is extensively tested on a grid model in MATLAB Simulink for different types of generation and grid operating conditions. The results demonstrate that the algorithm can operate quickly and reliably, and accurately estimate fault characteristics even in the presence of noisy measurements and uncertain line parameters.

Figures

Figures reproduced from arXiv: 1908.03753 by the authors.

Figure 1
Figure 1. Model of a power line with/without an internal fault [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 3
Figure 3. Measured current and its numerically estimated derivatives for [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 4
Figure 4. Outline of second step of proposed algorithm [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figures from the paper (3 more)
Figure 5
Figure 5. Figure 5: Grid model for evaluation of proposed algorithm [PITH_FULL_IMAGE:figures/full_fig_p005_5.png]
Figure 6
Figure 6. Figure 6: Impact of measurement noise on algorithm performance [PITH_FULL_IMAGE:figures/full_fig_p007_6.png]
Figure 7
Figure 7. Figure 7: Impact of uncertainty in line parameters on algorithm performance [PITH_FULL_IMAGE:figures/full_fig_p008_7.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

21 extracted references · 21 canonical work pages

  1. [1]

    Evaluation of travelling wave based protection schemes for implementation in medium voltage distribution systems

    M. hrstrm, M. Geidl, L. Sder, and G. Andersson, “Evaluation of travelling wave based protection schemes for implementation in medium voltage distribution systems”, in 18th International Conference and Exhibition on Electricity Distribution , 2005, pp. 1–5

  2. [2]

    Speed of line protection - can we break free of phasor limitations?

    E. O. Schweitzer, B. Kasztenny, A. Guzmn, V . Skendzic, and M. V . Mynam, “Speed of line protection - can we break free of phasor limitations?”, in 68th Annual Conference for Protective Relay Engineers , 2015, pp. 448–461

  3. [3]

    Two-terminal traveling-wave-based transmission-line protection

    F. B. Costa, A. Monti, F. V . Lopes, K. M. Silva, P. Jamborsalamati, and A. Sadu, “Two-terminal traveling-wave-based transmission-line protection”, IEEE Transactions on Power Delivery , vol. 32, no. 3, pp. 1382–1393, 2017

  4. [4]

    A new differential protection of transmission line based on equivalent travelling wave

    L. Tang, X. Dong, S. Luo, S. Shi, and B. Wang, “A new differential protection of transmission line based on equivalent travelling wave”, IEEE Transactions on Power Delivery , vol. 32, no. 3, pp. 1359–1369, 2017

  5. [5]

    New time-domain line protection principles and implementation

    E. O. Schweitzer, B. Kasztenny, M. Mynam, A. Guzmn, and V . Skendzic, “New time-domain line protection principles and implementation”, 8th Annual Protection, Automation and Control World Conference , 2017

  6. [6]

    An ultra-high-speed directional relay based on correlation of incremental quantities

    A. Lei, X. Dong, and V . Terzija, “An ultra-high-speed directional relay based on correlation of incremental quantities”, IEEE Transactions on Power Delivery, vol. 33, no. 6, pp. 2726–2735, 2018

  7. [7]

    Transmission-line protection: A directional comparison scheme using the average of superimposed components

    S. M. Hashemi, M. T. Hagh, and H. Seyedi, “Transmission-line protection: A directional comparison scheme using the average of superimposed components”, IEEE Transactions on Power Delivery , vol. 28, no. 2, pp. 955–964, 2013

  8. [8]

    Study of ultra-high-speed protection of transmission lines using a directional comparison scheme of transient energy

    Y . Kong, B. Zhang, and Z. Hao, “Study of ultra-high-speed protection of transmission lines using a directional comparison scheme of transient energy”, IEEE Transactions on Power Delivery, vol. 30, no. 3, pp. 1317– 1322, 2015

Show all 21 references
  1. [9]

    Traveling- wave and incremental quantity directional elements speed up directional comparison protection schemes

    A. Guzmn, M. V . Mynam, V . Skendzic, and J. L. Eternod, “Traveling- wave and incremental quantity directional elements speed up directional comparison protection schemes”, 14th International Conference on Developments in Power System Protection , 2018

  2. [10]

    Practical setting considerations for protective relays that use incremen- tal quantities and traveling waves

    B. Kasztenny, A. Guzmn, N. Fischer, M. V . Mynam, and D. Taylor, “Practical setting considerations for protective relays that use incremen- tal quantities and traveling waves”, 43rd Annual Western Protective Relay Conference, 2016

  3. [11]

    Setting-less protection: Feasibility study

    A. P. S. Meliopoulos, G. J. Cokkinides, Z. Tan, S. Choi, Y . Lee, and P. Myrda, “Setting-less protection: Feasibility study”, in 46th Hawaii International Conference on System Sciences , 2013, pp. 2345–2353

  4. [12]

    Dynamic state estimation based protection on series compensated transmission lines

    Y . Liu, A. P. S. Meliopoulos, R. Fan, L. Sun, and Z. Tan, “Dynamic state estimation based protection on series compensated transmission lines”, IEEE Trans. on Power Delivery , vol. 32, no. 5, pp. 2199–2209, 2017

  5. [13]

    Estimation of transmission line parameters using multiple methods

    E. C. M. Costa and S. Kurokawa, “Estimation of transmission line parameters using multiple methods”, Transmission Distribution IET Generation, vol. 9, no. 16, pp. 2617–2624, 2015

  6. [14]

    Online optimal transmission line parameter estimation for relaying applications

    Y . Liao and M. Kezunovic, “Online optimal transmission line parameter estimation for relaying applications”, IEEE Transactions on Power Delivery, vol. 24, no. 1, pp. 96–102, 2009

  7. [15]

    Protection and testing considerations for iec 61850 sampled values-based distance and line current differential schemes

    S. Chase, E. Jessup, M. Silveira, J. Dong, and Q. Yang, “Protection and testing considerations for iec 61850 sampled values-based distance and line current differential schemes”, 72rd Annual Conference for Protective Relay Engineers , 2019

  8. [16]

    Wavelet transforms in power systems. I. General introduction to the wavelet transforms

    Chul Hwan Kim and Raj Aggarwal, “Wavelet transforms in power systems. I. General introduction to the wavelet transforms”, Power Engineering Journal, vol. 14, no. 2, pp. 81–87, 2000

  9. [17]

    Tu- torial on single-pole tripping and reclosing

    E. Godoy, A. Celaya, H. J. Altuve, N. Fischer, and A. Guzmn, “Tu- torial on single-pole tripping and reclosing”, in 39th Annual Western Protective Relay Conference , 2012, pp. 1–22

  10. [18]

    Optical current sensors eliminate CT saturation

    J. D. P. Hrabluik, “Optical current sensors eliminate CT saturation”, 2002, vol. 2, pp. 1478–1481 vol.2

  11. [19]

    Gurobi optimizer reference manual

    LLC Gurobi Optimization, “Gurobi optimizer reference manual”, 2019

  12. [20]

    Progress on protection strategies to mitigate the impact of renewable distributed generation on distribution systems

    M. Norshahrani, H. Mokhlis, A. H. Abu Bakar, J. J. Jamian, and S. Sukumar, “Progress on protection strategies to mitigate the impact of renewable distributed generation on distribution systems”, Energies, vol. 10, no. 11, pp. 1864, 2017

  13. [21]

    Optical current sensors for high power systems: A review

    Ricardo M. Silva, Hugo Martins, Ivo Nascimento, Jos M. Baptista, Antnio Lobo Ribeiro, Jos L. Santos, Pedro Jorge, and Orlando Frazo, “Optical current sensors for high power systems: A review”, Applied Sciences, vol. 2, no. 3, pp. 602–628, 2012

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.