Cable Fault Location in VSC-HVDC System Based on Improved Local Mean Decomposition

: Aiming at the problem of low positioning accuracy caused by modal aliasing and noise interference in DC cable fault location analysis of a VSC-HVDC system, a double-ended fault location method for ﬂexible DC cables based on improved local mean decomposition is proposed. Firstly, the local mean decomposition (LMD) is used to decompose the six-mode voltage signal to obtain the product function (PF) component; then, to overcome the problem that the instantaneous frequency function of the LMD is limited by the extreme value, the Hilbert transform is performed on the PF1 to obtain the instantaneous frequency curve, and the arrival time of the voltage traveling wave head is determined from the mutation information. Finally, the fault distance is obtained by using the principle of double-ended traveling wave fault location. Different fault conditions are simulated, analyzed, and compared with wavelet transform and Hilbert–Huang transform. The results show that the proposed method has a positioning error within 1%, and it is less affected by interference noise and transition resistance. The three-phase fully controlled bridge circuit the SPWM modulation control circuit are used. There is constant DC voltage and constant AC voltage control; the receiving end adopts constant active power and constant reactive power control; the rated DC voltage is set to 640 kV; the rated transmission active power is 1000 MW; and the reactive power is 0 Mvar. The single-core XLPE DC cable model adopts the armored frequency domain model, and the total cable length is 200 km.


Introduction
In the context of the development of the dual-carbon policy, the large-scale and effective use of renewable energy and distributed energy has promoted the development of flexible DC transmission projects. China has also successfully built flexible DC transmission projects such as Shanghai Nanhui, Zhejiang Zhoushan, and Fujian Xiamen [1]. It has broad application scenarios in terms of island power supply and asynchronous interconnection [2]. As an important link for reliable and continuous output of electric energy, transmission lines are mostly laid in cable trenches or directly buried underground, which is more reliable than overhead lines. The line cannot be directly followed, and the repair speed and excavation range put forward higher requirements for the measurement of fault distance.
At present, the available methods for distance measurement of DC transmission lines include the traveling wave method, fault analysis method, and natural frequency method. From the three dimensions of time, space, and frequency, the distance equation between fault distance and electrical quantity is constructed. The fault analysis method establishes the ranging equation system according to the system parameters and the fault recording data and solves the fault distance [3]. Reference [4] calculates the resistance value of the transition resistance according to the voltage distribution along the line and the current distribution along the line and uses the minimum variance of the transition resistance value at the fault point to locate the fault. Reference [5] is based on the timedomain model, considering the distribution parameters of HVDC transmission lines with infinitely small multi-order distances, and uses the particle swarm algorithm to search for the corresponding distance with the smallest voltage distribution difference along the entire line, which is the actual fault distance. However, the line model is required to be process of multiple iterations, which will affect the detection accuracy of the disturbance amplitude and frequency. The local mean decomposition method [25] obtains the product function by iterative division, and the process of obtaining the product function component has fewer iterations and less end effect, and so the degree of inward contamination of the instantaneous amplitude and frequency data is lesser. However, for the positioning of the arrival time of the traveling wave head, it is often necessary to observe the amplitude mutation point of the instantaneous frequency curve. The value of the frequency-modulated (FM) signal must be kept between −1 and 1 to be effective, but the corresponding value of the wave head arrival time is often far greater. If the valid range is required, it should be set to −1 or 1. The instantaneous frequency curve obtained after processing has no obvious change in the sudden change time and amplitude from its adjacent time value and amplitude, and so accurate positioning cannot be carried out.
Based on the above analysis, an improved local mean decomposition method (ILMD) is proposed. Aiming at the problem that the instantaneous frequency function in the general local mean decomposition method is limited by the extreme value and cannot be accurately located, the method of extracting the instantaneous frequency function after the Hilbert transform of the product function component is proposed to solve the problem. The first frequency mutation point of the frequency curve is used to determine the moment when the first wave head of the traveling wave arrives at the measurement point. Simulation verification proves that the proposed ranging method has high accuracy.

DC Cable Phase-Mode Conversion
This article takes a single-core XLPE cable with an armored layer as an example. The single-core XLPE cable is composed of a cable core, an insulating layer, and an armoring layer. Regardless of whether they are in normal operation or the fault stage, the positive and negative poles have a strong coupling degree, and the six groups of voltage and current interfere with each other. To eliminate the difference between these parameters, they must be decoupled first. The voltages of the core layer, insulation layer, and armor layer of the DC cable are, respectively, recorded as U MPc , U Mps , U Mpa and U Nnc , U Nns , U Nna , and the DC currents of them are, respectively, recorded as I MPc , I Mps , I Mpa and I Nnc , I Nns , I Nna . The specific reference direction is shown in Figure 1.   The wave process equation of the cable is: Converting R, L, G, and C into a diagonal matrix gives the phase-mode transfor- The wave process equation of the cable is: Converting R, L, G, and C into a diagonal matrix gives the phase-mode transformation matrix of the DC voltage [26].
where 0 and I are the second-order zero matrix and identity matrix, respectively, and K is the Karenbauer transformation matrix. Let Q = P −T . Q is the phase-mode transformation matrix of the cable current; then the wave process equation of the cable is transformed into Among them, the modulo voltage and current vectors after phase-mode transformation are shown in Equations (5) and (6). The six sets of voltage and current parameters in the bipolar cable are converted into six independent moduli.

Traveling Wave Transmission Characteristics
When a short-circuit fault occurs somewhere in the cable, the magnetic field energy and electric field energy around the fault point are transformed into each other, and then voltage and current traveling waves are generated to the lines at both ends. According to the superposition principle, it can be equivalent to the normal component and the fault component voltage at the fault point. The additional voltage source generates fault voltage traveling waves with rich frequencies to both ends of the cable. When the fault voltage traveling wave propagates in the DC cable, its own refraction occurs at the fault point, and the polarity of the traveling wave will be changed by reflection instead of refraction. The grid is shown in Figure 2. voltage and current traveling waves are generated to the lines at both ends. According to the superposition principle, it can be equivalent to the normal component and the fault component voltage at the fault point. The additional voltage source generates fault voltage traveling waves with rich frequencies to both ends of the cable. When the fault voltage traveling wave propagates in the DC cable, its own refraction occurs at the fault point, and the polarity of the traveling wave will be changed by reflection instead of refraction. The grid is shown in Figure 2. In the DC cable, the characteristic impedance (wave impedance) Z cm of each modulus and the modulus propagation coefficient γ m are calculated as follows [27]: R m , G m , L m , C m are the modulus resistance, modulus conductance, modulus inductance, and modulus capacitance per unit length of the line, respectively.
Combining the above two formulas, we obtain: In the formula, α m is the attenuation constant of each modulus, characterizing the attenuation characteristics of the traveling wave of each modulus; β m is the phase constant of each modulus, characterizing the phase change characteristics of the traveling wave of each modulus; w is the signal angular frequency; w= 2fπ.
Then, the formula for calculating the frequency-dependent wave velocity (propagation along the line) of each modulus of the DC cable is: The fault traveling wave signal can be decomposed into traveling wave components in different frequency segments. From Equations (9) and (10), it can be known that when different frequency components of traveling waves with different moduli propagate on the transmission line, the attenuation constant and propagation speed are different. This will cause traveling wave attenuation and dispersion, which will distort the traveling wave waveform. The higher the frequency of the traveling wave component, the more serious is the attenuation; that is, during the propagation process, the highest frequency component of the initial traveling wave head has the most serious attenuation. Under the attenuation of the initial wave head, the frequency band of the second traveling wave head will be reduced (the farther the distance is, the more the reduction will be); the speed of the traveling wave component and the singularity of the wave head will be weakened. The travel time of the component's traveling wave to the measurement point increases. The detection of the second traveling wave head will directly cause errors in the ranging accuracy, and so using the double-ended traveling wave method to identify the first traveling wave head at both ends of the line will have a higher ranging accuracy.
If the total length of the line is L, it is assumed that, at the time t f , a short-circuit fault occurs at a distance d from the cable to the rectifier side, and the time when the initial traveling wave of the fault reaches the rectifier side and the inverter side is t M and t N , respectively. The propagation velocity of the voltage traveling wave in the cable is v; then the distance from the fault point to the rectifier terminal is shown in Equation (11).
The relative positioning error is defined as Equation (12). Among them, L m represents the measured fault distance, L a represents the actual fault distance, and L 0 represents the full length of the line.
In this paper, the six-mode voltage is used as the ranging parameter; the propagation speed and attenuation constant of the six-mode voltage in the cable are comprehensively considered, and the fixed reference wave speed is selected as v = 1.978 × 10 8 m/s.

The Principle of LMD
The local mean decomposition signal processing process [28] is as follows.
(1) For an original signal, determine all local extreme points in the signal, and calculate the average value of any two adjacent local extreme points and the estimated value of the amplitude envelope.
(2) Connect the adjacent local mean points m i and m i+1 and the amplitude envelope estimates a i and a i+1 , and then smooth them with the moving average method to obtain the local mean function m 11 (t) and the amplitude envelope estimation function a 11 (t), respectively. The calculation formula of moving average processing is as follows.
Among them, y(i) is the original sequence of data, and 2n + 1 is the sliding width. When i < n, the span should be reduced accordingly, and the limit should not exceed the endpoint of the sequence.
(3) Separate the local mean function m 11 (t) from the original signal x(t) to obtain the signal h 11 (t).
(4) Divide h 11 (t) by the magnitude envelope estimation function a 11 (t) to obtain the FM signal s 11 (t).
Repeat the above steps s 11 (t). If the amplitude envelope estimation function a 12 (t) = 1 is obtained, this means that s 11 (t) is a pure FM signal. If a 12 (t) = 1, the above iterative Processes 2022, 10, 1655 7 of 17 process needs to be repeated until s 1n (t) is a pure FM signal, as shown in the following formula.
The conditions for iterative termination are as shown in Equation (20). In practical applications, in order to avoid excessive decomposition times, an iteration termination threshold ∆ is set.
(5) Multiply all the amplitude envelope estimation functions generated in the iterative process to obtain the instantaneous amplitude function a 1 (t).
(6) Multiply a 1 (t) and s 1n (t) to obtain the first product function component of the original signal x(t).
(7) After separating the first product function component from the original signal x(t), the remaining signal u 1 (t) is obtained, and the above steps are repeated with u 1 (t) as the new initial signal until u k (t) is a monotonic function. The calculation process is as follows.
The final original signal can be reconstructed from u k (t) and all product function components.

Hilbert Transform
Perform Hilbert transform on the product function components of the original signal x(t).
The corresponding analytical signal can be expressed by Equations (26) and (27).
Processes 2022, 10, 1655 Derivation of the phase function gives the instantaneous frequency function. Figure 3 is the flowchart of the fault location of DC cables based on the improved LMD algorithm. The input fault signal is first decomposed by local mean value; the obtained amplitude envelope signal A is used as the final instantaneous amplitude function, and the PF1 component is used for Hilbert transform to obtain the instantaneous frequency function. The sudden change point of the instantaneous frequency curve is the moment when the wave head reaches the measurement end, and then the fault distance is calculated.
Derivation of the phase function gives the instantaneous frequency function. Figure 3 is the flowchart of the fault location of DC cables based on the improved LMD algorithm. The input fault signal is first decomposed by local mean value; the obtained amplitude envelope signal A is used as the final instantaneous amplitude function, and the PF1 component is used for Hilbert transform to obtain the instantaneous frequency function. The sudden change point of the instantaneous frequency curve is the moment when the wave head reaches the measurement end, and then the fault distance is calculated.

Comparison of Mathematical Methods
There is a significant difference in the solution methods for the local mean functions EMD and LMD. EMD uses cubic spline interpolation for all extreme points to obtain the upper envelope and lower envelope of the original signal and then uses the average method to obtain the local mean function. For LMD to solve the local mean function, the

Comparison of Mathematical Methods
There is a significant difference in the solution methods for the local mean functions EMD and LMD. EMD uses cubic spline interpolation for all extreme points to obtain the upper envelope and lower envelope of the original signal and then uses the average method to obtain the local mean function. For LMD to solve the local mean function, the average of two adjacent extreme values is obtained, and the moving average algorithm is used to smooth it; we take a mathematical signal as an example and use LMD and EMD to decompose and then reconstruct. Observe Figures 4c and 5b; compared with the original signal, the frequency and amplitude after EMD resynthesis are quite different, and the decomposition result of LMD is more accurate.
The LMD decomposition results are as follows: average of two adjacent extreme values is obtained, and the moving average algorithm is used to smooth it; we take a mathematical signal as an example and use LMD and EMD to decompose and then reconstruct. Observe Figure 4c and Figure 5b; compared with the original signal, the frequency and amplitude after EMD resynthesis are quite different, and the decomposition result of LMD is more accurate.
The LMD decomposition results are as follows:  The EMD decomposition results are as follows: Processes 2022, 10, 1655 11 The EMD decomposition results are as follows:

Power Transmission Model Construction
At present, there are three types of HVDC transmission model: distributed param model, equivalent thermal circuit model, and finite element model. When the wavele of the electromagnetic wave is much smaller than the physical size of the correspon network, it is necessary to use the distributed parameter model to characterize the po network, and so the distributed parameter model is more used in practical applicat The VSC-HVDC simulation system model established by PSCAD is shown in Figu The AC system voltage is 50 Hz; the rated voltage is 220 kV; the rated DC voltage is kV; and the rated transmission power is 1000 MW. The three-phase fully controlled br circuit and the SPWM modulation control circuit are used. There is constant DC vo and constant AC voltage control; the receiving end adopts constant active power and stant reactive power control; the rated DC voltage is set to 640 kV; the rated transmis active power is 1000 MW; and the reactive power is 0 Mvar. The single-core XLPE cable model adopts the armored frequency domain model, and the total cable leng 200 km.

Power Transmission Model Construction
At present, there are three types of HVDC transmission model: distributed parameter model, equivalent thermal circuit model, and finite element model. When the wavelength of the electromagnetic wave is much smaller than the physical size of the corresponding network, it is necessary to use the distributed parameter model to characterize the power network, and so the distributed parameter model is more used in practical applications. The VSC-HVDC simulation system model established by PSCAD is shown in Figure 6. The AC system voltage is 50 Hz; the rated voltage is 220 kV; the rated DC voltage is ±320 kV; and the rated transmission power is 1000 MW. The three-phase fully controlled bridge circuit and the SPWM modulation control circuit are used. There is constant DC voltage and constant AC voltage control; the receiving end adopts constant active power and constant reactive power control; the rated DC voltage is set to 640 kV; the rated transmission active power is 1000 MW; and the reactive power is 0 Mvar. The single-core XLPE DC cable model adopts the armored frequency domain model, and the total cable length is 200 km. Among them, ① represents the AC power supply; ② represents the converter transformer; ③ represents the converter reactor; ④ represents the AC filter; ⑤ represents the DC side capacitor; and ⑥ represents the DC cable.

ILMD Method Validation
Using the above-mentioned building model, faults are set at different positions, and the improved LMD method is used for testing. Set the fault start time to 1.5 s, the fault duration to 0.02 s, and the sampling frequency to 0.1 MHz. Taking the metal ground fault 130 km away from the rectifier side as an example, take 1.45-1.55 s for data analysis. The corresponding fault current and fault voltage waveforms are shown in Figure 7, and the six-mode voltage waveform obtained by decoupling the rectifier side voltage is shown in Figure 8.   Among them, 1 represents the AC power supply; 2 represents the converter transformer; 3 represents the converter reactor; 4 represents the AC filter; 5 represents the DC side capacitor; and 6 represents the DC cable.

ILMD Method Validation
Using the above-mentioned building model, faults are set at different positions, and the improved LMD method is used for testing. Set the fault start time to 1.5 s, the fault duration to 0.02 s, and the sampling frequency to 0.1 MHz. Taking the metal ground fault 130 km away from the rectifier side as an example, take 1.45-1.55 s for data analysis. The corresponding fault current and fault voltage waveforms are shown in Figure 7, and the six-mode voltage waveform obtained by decoupling the rectifier side voltage is shown in Figure 8. Among them, ① represents the AC power supply; ② represents the converter transformer; ③ represents the converter reactor; ④ represents the AC filter; ⑤ represents the DC side capacitor; and ⑥ represents the DC cable.

ILMD Method Validation
Using the above-mentioned building model, faults are set at different positions, and the improved LMD method is used for testing. Set the fault start time to 1.5 s, the fault duration to 0.02 s, and the sampling frequency to 0.1 MHz. Taking the metal ground fault 130 km away from the rectifier side as an example, take 1.45-1.55 s for data analysis. The corresponding fault current and fault voltage waveforms are shown in Figure 7, and the six-mode voltage waveform obtained by decoupling the rectifier side voltage is shown in Figure 8.   Among them, ① represents the AC power supply; ② represents the converter transformer; ③ represents the converter reactor; ④ represents the AC filter; ⑤ represents the DC side capacitor; and ⑥ represents the DC cable.

ILMD Method Validation
Using the above-mentioned building model, faults are set at different positions, and the improved LMD method is used for testing. Set the fault start time to 1.5 s, the fault duration to 0.02 s, and the sampling frequency to 0.1 MHz. Taking the metal ground fault 130 km away from the rectifier side as an example, take 1.45-1.55 s for data analysis. The corresponding fault current and fault voltage waveforms are shown in Figure 7, and the six-mode voltage waveform obtained by decoupling the rectifier side voltage is shown in Figure 8.   The six-mode voltage is decomposed by LMD to obtain the high-frequency component PF1 as shown in Figure 9a, and the Hilbert analysis is performed to obtain its instantaneous frequency function curve, as shown in Figure 9b. The frequency mutation points of the Mterminal and N-terminal correspond to the moments of 1.50065 s and 1.50035 s, respectively. Putting them into Equation (11), the theoretical value of the fault distance is 129.6699 km; the error distance is 0.33 km, and the error percentage is 0.165%. The six-mode voltage is decomposed by LMD to obtain the high-frequency component PF1 as shown in Figure 9a, and the Hilbert analysis is performed to obtain its instantaneous frequency function curve, as shown in Figure 9b. The frequency mutation points of the M-terminal and N-terminal correspond to the moments of 1.50065 s and 1.50035 s, respectively. Putting them into Equation (11), the theoretical value of the fault distance is 129.6699 km; the error distance is 0.33 km, and the error percentage is 0.165%. The IMF component obtained by using HHT for the six-mode voltage component is shown in Figure 10a, and the instantaneous frequency of Hilbert transform of IMF1 is shown in Figure 10b. The frequency mutation points of the M-terminal and N-terminal appear at 1.50066 s and 1.50035 s, respectively. Putting them into Equation (11), the theoretical value of fault distance is 130.6590 km; the error distance is 0.659 km, and the error percentage is 0.3295%. Compared with HHT, ILMD has a high precision of 329 m, and, in the case of partial fault distance, the amplitude-frequency information, including the highfrequency feature quantity, is seriously polluted in the process of EMD decomposition, which makes it impossible to accurately locate. The IMF component obtained by using HHT for the six-mode voltage component is shown in Figure 10a, and the instantaneous frequency of Hilbert transform of IMF1 is shown in Figure 10b. The frequency mutation points of the M-terminal and N-terminal appear at 1.50066 s and 1.50035 s, respectively. Putting them into Equation (11), the theoretical value of fault distance is 130.6590 km; the error distance is 0.659 km, and the error percentage is 0.3295%. Compared with HHT, ILMD has a high precision of 329 m, and, in the case of partial fault distance, the amplitude-frequency information, including the high-frequency feature quantity, is seriously polluted in the process of EMD decomposition, which makes it impossible to accurately locate.
In the fault location problem, the initial traveling wave front is the singular point, and the wavelet mode maxima can be used to identify the singular point. In the fault data in Figure 6, the wavelet modulus maxima are obtained by discrete stationary wavelet transform. The wavelet basis function and the number of decomposition layers are Daubechies (db4) and three layers, respectively, and the modulus maxima of the voltage at both ends are obtained as shown in Figure 11. In the fault location problem, the initial traveling wave front is the singular point, and the wavelet mode maxima can be used to identify the singular point. In the fault data in Figure 6, the wavelet modulus maxima are obtained by discrete stationary wavelet transform. The wavelet basis function and the number of decomposition layers are Daubechies (db4) and three layers, respectively, and the modulus maxima of the voltage at both ends are obtained as shown in Figure 11.  In the fault location problem, the initial traveling wave front is the singular point, and the wavelet mode maxima can be used to identify the singular point. In the fault data in Figure 6, the wavelet modulus maxima are obtained by discrete stationary wavelet transform. The wavelet basis function and the number of decomposition layers are Daubechies (db4) and three layers, respectively, and the modulus maxima of the voltage at both ends are obtained as shown in Figure 11.  It can be seen from Figure 11 that the time when the initial traveling wave head reaches the M end and the N end corresponds to the 5069th and 5038th sampling points, respectively, and the corresponding time difference between the two sampling points is 0.00031 s. It can be calculated that the theoretical value of the fault distance is 130.6605 km; the error distance is 0.6605 km, and the relative error percentage is 0.3302%.
For metallic ground faults, the fault time is set to 1.5 s~1.52 s; the fault location is performed for different fault distance points; and the accuracy is compared with the methods using wavelet analysis and HHT. It can be obtained that the error of fault location using ILMD is significantly reduced compared with that of other methods. The corresponding experimental simulation data are shown in Table 1. Taking China as an example, according to the 'National Electric Power Dispatching System "Tenth Five-Year" Science and Technology Development Plan Outline' issued by State Grid, the comprehensive error should not exceed 1%. Due to technical and economic constraints, the ranging error cannot be too small. From a practical point of view, an absolute ranging error of less than 1 km can better meet the field requirements. The error in the ranging method proposed in this paper is up to more than 900 m under high impedance, and the error is smaller in most cases, which can meet the actual requirements of engineering. The fault time is set as 1.5 s~1.52 s, and it is grounded through different transition resistances under different fault distances. The ILMD distance measurement results are shown in Table 2. It can be seen from Table 2 that, at the same fault point, the error only increases by a few meters in the case of high-impedance grounding; that is, the result of fault location is hardly affected by the transition resistance.

Noise Immunity Test
Cables generally use transformers to collect fault traveling wave signals, but the detection of transformers is susceptible to noise interference. The influencing factors and sources of noise can be divided into external factors and internal factors. External factors include line construction, electrical equipment interference, etc., while internal factors generally originate from random fluctuations within the transformer itself.
Take the metallic ground fault 130 km away from the rectifier side as an example, and add noise interference with a signal-to-noise ratio (SNR) of 20 dB to the voltage traveling wave signal, as shown in Figure 12a. Gaussian white noise with signal-to-noise ratios of 10 dB, 15 dB, 25 dB, and 40 dB is added to the fault signal, and the instantaneous frequency is obtained by ILMD analysis, as shown in Figure 12b. The 2037th and 2067th sampling points appear at the M-terminal and N-terminal frequency mutation points, respectively, corresponding to 1.50065 s and 1.50034 s. Putting them into formula (11), the theoretical value of the fault distance can be obtained as 129.6699 km, and the error is 0.3301 km. A large number of fault points are added with different proportions of noise for testing, and the results are shown in Table 3. When the signal-to-noise ratio (SNR) is less than 10 dB, the instantaneous frequency mutation point caused by the arrival of the wave head will be submerged, and accurate detection cannot be performed. When the signal-to-noise ratio SNR ≥ 15 dB, the ranging accuracy is high.

Noise Immunity Test
Cables generally use transformers to collect fault traveling wave signals, but the detection of transformers is susceptible to noise interference. The influencing factors and sources of noise can be divided into external factors and internal factors. External factors include line construction, electrical equipment interference, etc., while internal factors generally originate from random fluctuations within the transformer itself.
Take the metallic ground fault 130 km away from the rectifier side as an example, and add noise interference with a signal-to-noise ratio (SNR) of 20 db to the voltage traveling wave signal, as shown in Figure 12a. Gaussian white noise with signal-to-noise ratios of 10 dB, 15 dB, 25 dB, and 40 dB is added to the fault signal, and the instantaneous frequency is obtained by ILMD analysis, as shown in Figure 12b. The 2037th and 2067th sampling points appear at the M-terminal and N-terminal frequency mutation points, respectively, corresponding to 1.50065 s and 1.50034 s. Putting them into formula (11), the theoretical value of the fault distance can be obtained as 129.6699 km, and the error is 0.3301 km. A large number of fault points are added with different proportions of noise for testing, and the results are shown in Table 3. When the signal-to-noise ratio (SNR) is less than 10 dB, the instantaneous frequency mutation point caused by the arrival of the wave head will be submerged, and accurate detection cannot be performed. When the signal-to-noise ratio SNR ≥ 15 dB, the ranging accuracy is high.

Conclusions
(1) This paper proposes a new fault location method for HVDC cable lines which combines LMD and Hilbert transform. It not only inherits the advantages of LMD in the decomposition process of the weak endpoint effect and of less pollution of amplitudefrequency information but also takes into account the ability of the Hilbert transform to completely retain the mutation information of the instantaneous frequency function. Through a large number of experimental verifications, the error of the proposed positioning method is within 1%; the positioning accuracy is high, and it is hardly affected by transition resistance and interference noise.
(2) Compared with HHT, the method in this paper has a small endpoint effect and less iteration processing, which avoids the large difference in the instantaneous frequency function extracted by HHT in some cases. Compared with wavelet transform, it overcomes the large positioning error in the case of high-resistance near-end faults, and the ranging accuracy meets the requirements of on-site operation and maintenance accuracy.