Method Article

Power Distribution Cable Defect Localization Technology Based on the Maximum Entropy Spectral Method

DOI:

10.3791/70701

June 5th, 2026

In This Article

Summary

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

This protocol presents a maximum-entropy spectral-estimation-based method for precise defect localization in long power distribution cables using frequency-domain reflectometry. It details high-frequency cable modeling, reflection coefficient spectrum analysis, spectral extrapolation, and validation via simulations and experiments, achieving superior accuracy compared to conventional FFT approaches.

Abstract

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

As cable laying distances in modern urban distribution networks continue to expand, the accuracy of fault location via frequency-domain reflectometry (FDR) diminishes markedly with increasing length. To improve precision for long cables, this study proposes a maximum-entropy spectral-estimation-based method. A high-frequency distributed parameter model for distribution cables is developed to quantitatively assess the impacts of insulation aging and mechanical damage on per-unit-length capacitance and inductance. Drawing on transmission line theory, the relationship between the reflection coefficient spectrum and the defect position is derived, along with a step-frequency optimization criterion to enhance spectral resolution. To mitigate spectral leakage and limited resolution in fast Fourier transform (FFT)-based frequency-domain analysis, the maximum entropy spectral method is employed for high-fidelity spectrum estimation, thereby elevating fault distance localization accuracy. Simulations reveal that, relative to conventional FFT, the proposed method reduces location errors by 2–4.5 times in typical long-cable cases. Experimental results confirm a relative error below 0.25% for the maximum entropy approach, surpassing the conventional method's error under 0.55%, thus validating its efficacy and superiority in engineering practice.

Introduction

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

With the increasing demand for power supply reliability and safety in modern power systems, cables have been widely adopted in urban distribution networks. Among them, cross-linked polyethylene (XLPE) cables are extensively utilized due to the author’s outstanding advantages1,2. Distribution cables are typically installed underground, where the operating environment is harsh. As service time increases, combined with load fluctuations and overloading, cable aging becomes more pronounced. Consequently, localized defects such as water trees, electrical trees, partial discharges, and insulation moisture may gradually occur. If not addressed in time, these defects may evolve into severe failures, eventually leading to power outages in the grid3,4. Therefore, defect identification and fault location of distribution cables have become essential components of power grid operation and maintenance.

The fault location methods for distribution cables can be broadly categorized into two types: online and offline approaches. Online methods are mainly represented by traveling-wave location, partial discharge (PD) monitoring, and the common-mode leakage current method. Traveling-wave location identifies fault positions by analyzing transient traveling waves generated at the fault. However, it faces several challenges, including difficulties in accurately extracting wavefronts, severe superposition of fault features with nonlinear distribution, and inconvenient installation. More importantly, this method is not applicable for locating aging-related defects5. PD-based online monitoring detects defect signals generated inside the cable and enables real-time assessment of operating conditions. Nevertheless, it still encounters technological bottlenecks in recognizing aging defects and performing long-term monitoring. Problems such as signal interference, noise suppression, and complex data analysis often lead to misjudgments in practical applications6,7,8. The common-mode leakage current method determines the position and severity of aging defects by measuring variations in common-mode leakage currents at the cable terminals. This technique offers non-intrusive detection with high sensitivity. However, its anti-interference capability is relatively weak, as high-frequency differential-mode currents, such as those generated during variable-frequency equipment operation, can significantly reduce measurement accuracy. In addition, online monitoring of energized cables must address signal isolation in high-voltage environments, which increases equipment costs and maintenance complexity9,10,11.

In terms of offline methods, several approaches have been proposed, including time domain reflectometry (TDR)12, frequency domain reflectometry (FDR)13, and time-frequency domain reflectometry (TFDR)14 Among them, TDR suffers from low energy in the high-frequency components of the injected pulse signal and is highly susceptible to cable attenuation effects and noise interference, resulting in limited location accuracy and poor capability in detecting incipient defects at early stages of cable faults15. TFDR, by injecting Gaussian chirp pulses, mitigates the deficiency of insufficient high-frequency components in TDR16. However, TFDR requires complex signal modulation circuits and extraction algorithms, and its high-frequency modulated signals attenuate rapidly, making defect detection in long cables challenging17. FDR, on the other hand, injects swept-frequency signals with equal power distribution. It contains the richest high-frequency components, achieves higher location accuracy, and demonstrates greater sensitivity compared with TDR, which has led to its rapid development. Depending on the type of acquired signal, FDR can be further categorized into broadband impedance spectrum (BIS) methods18,19,20 and reflection coefficient spectrum (RCS) methods21,22,23,24.

In 2003, one group first proposed a cable local defect localization method based on the input reflection coefficient spectrum of cables25. However, this approach did not account for the frequency-dependent effects of distributed parameters, which may lead to misjudgments. Later, an inverse fast Fourier transform (IFFT) was applied to extract information from the reflection coefficient spectrum at the cable input, enabling effective localization of local defects such as thermal aging, radiation aging, and extrusion deformation26. It was also observed that the higher the upper frequency of the sweep signal, the higher the spatial resolution of frequency-domain reflectometry (FDR). Nevertheless, when applied to long cables, the upper frequency of the sweep signal decreases significantly, resulting in reduced localization accuracy.

The core of the FDR method lies in performing the FFT on the acquired reflection signals. However, due to the inherent spectral leakage problem of the FFT, researchers have proposed several improvements. A study introduced the use of a distance window combined with a Kaiser window function in the discrete Fourier transform (DFT) to extract the defect components from the real part of the reflection coefficient spectrum, thereby enabling the localization of loosened copper shielding defects in cables27. A Kaiser window-based DFT28 was applied to the real part of the input reflection coefficient spectrum, achieving fault-type identification for open-circuit, short-circuit, high-resistance, and low-resistance conditions. A study proposed applying a Blackman window-based DFT to the imaginary part of the cable input impedance for the detection and localization of moisture in intermediate joints29. An interpolation algorithm using the Hanning self-convolution window in the DFT was introduced to estimate the frequency and phase of periodic components in the real part of the reflection coefficient spectrum30, which enabled polarity determination of cable impedance anomalies. Subsequently, improvements were achieved by incorporating Nuttall windows, fourth-order three-term Nuttall windows, second-order Nuttall self-convolution windows, and three-point interpolated DFT algorithms31. Overall, existing studies have mainly focused on window function design or interpolation algorithm refinement, which have enhanced defect localization accuracy to a certain extent.

From the perspective of the FFT principle, in frequency-domain reflectometry (FDR), the higher the frequency of the injected signal, the better the measurement accuracy. In theory, if the injection frequency approaches infinity, the cable defect can be located with perfect precision. However, the transmission of ultra-high-frequency signals along cables experiences severe attenuation, making them unsuitable for medium- and long-length cables32,33. Therefore, there exists an inherent trade-off between the maximum injection frequency in FDR and the achievable defect location accuracy. To address this issue, this paper proposes a cable defect localization technique based on the maximum entropy spectral method. The method predicts the high-frequency portion of the reflection coefficient spectrum from the measured low-frequency data according to the maximum entropy criterion, and then computes the maximum entropy spectrum(MES) of the reflection coefficient to determine the defect location. Compared with conventional FFT-based FDR methods, which are limited by spectral leakage and resolution constraints, especially in long-distance cable applications, the proposed maximum entropy spectral method provides higher spectral resolution by effectively extrapolating high-frequency information from limited measured data. Unlike window-function-based improvements that mainly suppress sidelobes without fundamentally enhancing resolution, the MES approach improves both peak sharpness and localization accuracy. In addition, compared with interpolation-based or parametric spectral estimation methods, the MES method offers a better balance between resolution and robustness, making it particularly suitable for detecting incipient defects in long distribution cables. This approach improves the localization accuracy, reducing the relative error from 0.55% in conventional methods to 0.25%, while producing smoother localization curves that facilitate the identification of defect peaks.

Access restricted. Please log in or start a trial to view this content.

Protocol

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

NOTE: Perform all measurements on the test cable in an offline state (de-energized). Use a vector network analyzer (VNA) or equivalent frequency-domain reflectometry instrument capable of swept-frequency signal injection and reflection coefficient acquisition. Ensure the cable is open-circuited at the far end unless otherwise specified. Minimize connection reflections by using short, direct coaxial connections and conductive tape for secure contact. Refer to the Table of Materials for recommended equipment and software.

1. Prepare the cable sample and measurement setup

  1. Select the power distribution cable for testing (e.g., YJLV22-8.7/10 kV-3×70 XLPE cable or equivalent).
  2. Measure and record the total cable length and known defect positions (if artificially introduced).
  3. Measure the electromagnetic wave propagation velocity in the cable.
    1. Connect the cable to the instrument at one end and leave the far end open.
    2. Acquire the reflection coefficient spectrum over a suitable frequency band.
    3. Identify the oscillation period corresponding to the full cable length (round-trip).
    4. Calculate velocity as Wave velocity equation \(v=2 \times L \times f\); physics formula; educational use.   (1)
  4. Prepare the cable head connection.
    1. Strip the cable end to expose the core and shield.
    2. Use double-sided conductive copper tape to securely attach the instrument leads to the core and shield.
    3. Minimize branch lead length to reduce head-end reflections.

2. Configure the frequency sweep parameters

  1. Determine the frequency step (Δf) using the selection criterion.
    1. Use the cable-end oscillation period Equation for wave travel time, Tend=2L/v, in kinematics study, symbol representation.    (2)
    2. Set   Frequency change equation Δf = Tend/200  to ensure at least 200 sampling points per cycle.
    3. The frequency step Δf determines the maximum observable distance range:Maximum distance formula \(d_{max}=\frac{v}{2\Delta f}\), key in signal processing analysis.   (3)
    4. To ensure that the full cable length L can be resolved without ambiguity:
      Frequency uncertainty equation Δf<v/2L describes wavelength limits in optical systems.    (4)
    5. In this study, a more stringent condition is adopted to improve resolution: Frequency stability equation, Δf<1/200T_end, critical for signal processing analysis.    (5)
  2. Set the frequency band.
    1. Set the lower limit to 100 kHz (or instrument minimum).
    2. Set the upper limit to the maximum attainable without excessive attenuation (typically 30–50 MHz for cables <200 m; reduce for longer cables).
    3. The upper frequency limit fmax directly affects spatial resolution:Distance resolution inequality, formula: Δd < v/2f_max, used in signal processing analysis.    (6)
    4. A higher fmax leads to finer resolution, but is constrained by cable attenuation.
  3. Calculate the required number of sampling pointsFrequency calculation; equation N=(f_upper-f_lower)/(Δf+1); mathematical formula.    (7)
  4. Configure the instrument with the calculated frequency band, step, and points (e.g., 100 kHz–50 MHz, ~5,000 points for high resolution).

3. Acquire the reflection coefficient spectrum

  1. Calibrate the instrument (open-short-load calibration recommended for VNA).
  2. Connect the prepared cable head to the instrument port.
  3. Inject swept-frequency sinusoidal signals and measure the input reflection coefficient Γ(ω) at each discrete frequency point.
    1. The measured reflection coefficient can be expressed as:Frequency-dependent reflection coefficient formula, Γ(ω)=Γ₀·e⁻ʲ²βd, used in wave analysis.    (8)
    2. where: Γd: reflection coefficient at the defect, Beta ratio formula β=ω/v, shown as an equation relevant in physics studies.: phase constant, d: defect distance.
    3. This relationship shows that the phase of the reflection coefficient varies linearly with frequency, forming the basis for distance estimation.
  4. Record the complex reflection coefficient spectrum (magnitude and phase, or real and imaginary parts).
  5. Save the data as frequency vs. real/imaginary reflection coefficient.

4. Process the data using the conventional FFT method (for comparison)

  1. Extract the imaginary part of the reflection coefficient spectrum.
  2. Apply a window function (e.g., Blackman window) to reduce sidelobe leakage.
  3. Perform inverse fast Fourier transform (IFFT) on the windowed imaginary part.
    1. The inverse Fourier transform is introduced as:Inverse Fourier transform equation, x(t)=F⁻¹{Γ(ω)}, mathematical concept, signal processing.    (9)
  4. The defect distance is calculated as:
    Kinematic equation, d=vt/2, solving uniform acceleration; physics formula illustration.    (10)
    1. where t is the time corresponding to the peak position after the inverse Fourier transform.
  5. Identify defect peaks (excluding head-end and terminal reflections) in the reflection coefficient spectrum.

5. Apply the maximum entropy spectral method for defect localization

  1. Load the acquired reflection coefficient spectrum data.
  2. Detrend and center the imaginary part sequence.
  3. Select the autoregressive model order
    M ≈ N/3 equation, illustrating static equilibrium principles in a mathematical context.    (11)
    1. where N is the number of data points.
  4. Implement the Burg recursion algorithm.
    1. Initialize forward and backward prediction errors.
    2. Recursively compute reflection coefficients and autoregressive parameters while minimizing forward and backward error power.
    3. Extrapolate the autocorrelation sequence based on the maximum entropy criterion.
  5. Compute the power spectrum from the extrapolated sequence.
  6. Identify peaks in the maximum entropy power spectrum.
  7. Convert peak positions to defect distances using Equation for distance in spectroscopy, d = (peak position × v) / 2, concept in wave analysis.  (12)
  8. Compare peak clarity and localization error with the conventional FFT results.

6. Interpret the results

  1. Distinguish positive peaks (successful defect detection) from spurious peaks or suboptimal traces.
  2. Calculate relative localization error as Equation for calculating percentage error in measurements.   (13)
  3. Confirm success by relative error <0.25% and a clear, smooth peak without significant sidelobes.

Access restricted. Please log in or start a trial to view this content.

Results

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

The reflection coefficient spectrum and defect localization results were obtained from both simulations and experiments using the maximum entropy spectral method and compared with conventional FFT-based approaches. A successful defect localization is characterized by a smooth power spectrum with a sharp, prominent peak at the exact defect position, minimal spurious peaks (pseudopeaks caused by sidelobes or leakage), and a low relative error (typically <0.3%). In contrast, suboptimal results from conventional FFT metho...

Access restricted. Please log in or start a trial to view this content.

Discussion

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

In this study, the authors present a detailed protocol for localizing defects in power distribution cables using frequency-domain reflectometry (FDR) enhanced by the maximum entropy spectral (MES) method. Previous studies have often relied on conventional fast Fourier transform (FFT)-based analysis, which struggles with spectral leakage, limited resolution in long cables, and spurious peaks that complicate defect identification34. Additionally, researchers without advanced signal processing expert...

Access restricted. Please log in or start a trial to view this content.

Disclosures

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

The authors have no conflicts of interest to declare that are relevant to the content of this article.

Acknowledgements

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

This work was supported by the Science and Technology Project of Guangxi Power Grid Co., Ltd. under Project No. 040100KC23110005.

Access restricted. Please log in or start a trial to view this content.

Materials

List of materials used in this article
NameCompanyCatalog NumberComments
Conductive copper tape, double-sidedNot specifiedN/AUsed for connecting test leads to the cable core and shield to minimize head-end reflection interference
https://www.shcenxin.com/metaltape/93.html
YJLV22-8.7/10 kV-3×70 medium-voltage cableNot specifiedYJLV22-8.7/10 kV-3×70Cross-linked polyethylene (XLPE) aluminum core steel tape armored PVC sheathed power cable; total length 171 m; used as test object with externally induced defect
http://www.xmdl518.com/jspn01-Products-6234877/
10 kV cable defect localization deviceNot specifiedNot specifiedFrequency-domain reflectometry instrument; frequency band 100 kHz–50 MHz; configured with 5,001 sampling points for reflection coefficient spectrum acquisition
https://www.163.com/dy/article/K6OHFJOO0552X57P.html

Reprints and Permissions

Request permission to reuse the text or figures of this JoVE article

Request Permission

Tags

Power Cable DefectsFault LocalizationMaximum Entropy MethodSpectral EstimationFrequency Domain ReflectometryTransmission Line TheoryReflection Coefficient SpectrumStep Frequency OptimizationInsulation AgingMechanical Damage

Related Articles