LINE DOUBLE-END STEADY-STATE QUANTITY DISTANCE MEASURING METHOD AND SYSTEM BASED ON AMPLITUDE-COMPARISON PRINCIPLE

Abstract

A line double-end steady-state quantity distance measuring method and system based on an amplitude-comparison principle. According to the method and system, voltage values and current values of both sides of a line before and after a fault are collected (102), a voltage variable quantity and a current variable quantity of both sides of the line are calculated (103), and after a voltage phasor value and a current phasor value are determined according to the voltage variable quantity and the current variable quantity (104), the position of a short-circuit point is determined by performing iterative calculation on the voltage of the short-circuit point. The method is simple in principle, and can accurately recognize a fault point, achieving precise distance measurement of lines.

Description

CROSS-REFERENCE TO RELATED APPLICATION

This application claims priority to China patent application No. 202010053618.X, filed with China national intellectual property administration on Jan. 17, 2020, the disclosure of which is hereby incorporated by reference in its entirety.

TECHNICAL FIELD

The present disclosure relates to the field of relay protection, and relates to, for example, a line double-end steady-state quantity distance measuring method and system based on an amplitude-comparison principle.

BACKGROUND

In areas with diverse and complex topography, transmission lines are usually spread between mountains and rivers, and are often set in a manner of overhead lines. In addition, in response to a special situation of transmission lines crossing over ultra-wide waterways and straits, there has also been an ultra-high-voltage overhead line-cable hybrid line. With the complexity of a power transmission network, quickly and accurately determining a position of a fault point of a high-voltage line and timely solving the problems and eliminating potential safety hazards are of great significance to ensure the safety of a power system.

For detection of high-voltage transmission line faults, methods in related technologies have been proposed at home and abroad, such as an impedance method based on power frequency electrical quantities, a fault analysis method, and a traveling wave method based on transient-state components. However, since electrical parameters of cables are different from those of overhead transmission lines, the fault characteristics of the hybrid line are different from those of a conventional line to a certain extent. Therefore, the accuracy of a fault distance measuring method for an overhead transmission line in the related technology will be correspondingly affected, such that a new distance measuring method needs to be provided.

SUMMARY

The present disclosure provides a line double-end steady-state quantity distance measuring method and system based on an amplitude-comparison principle, which can solve the technical problem that there is a relatively large error in fault positioning due to the fact that the fault positioning carried out on a high-voltage overhead line-cable hybrid line based on power frequency electrical quantities is easily affected by external unstable factors.

The present disclosure provides a line double-end steady-state quantity distance measuring method based on an amplitude-comparison principle. The method includes:

at step 1, collecting voltage values and current values, after an overhead line-cable hybrid transmission line is faulted, of both sides of the line, and voltage values and current values of both sides of the line of one cycle before the line is faulted; here both sides of the line are respectively a side M and a side N;

at step 2, determining voltage value variations according to the collected voltage values of both sides of the line before and after the line is faulted, and determining current value variations according to the collected current values of both sides of the line before and after the line is faulted;

at step 3, calculating voltage phasor values of both sides of the line by performing Fourier transformation on the voltage value variations of both sides of the line, and calculating current phasor values of both sides of the line by performing the Fourier transformation on the current value variations of both sides of the line;

at step 4, calculating voltages Δ{dot over (U)}φ_Mxi and Δ{dot over (U)}_φNxi of a compensation point according to a distance x_i from the compensation point to the side M of the line after the line is faulted, the voltage phasor values and the current phasor values of both sides of the line, a length L₁ of an overhead line on the side M, a length L₃ of an overhead line on the side N, a length L₂ of a cable, a wave impedance Z_cT and a propagation coefficient γ_T of the overhead line, and a wave impedance Z_cC and a current propagation coefficient γ_C of the cable; here an initial value of i is 1, and φ is any phase in a three-phase circuit, φ=A,B,C;

at step 5, determining a distance x_i+1 from the compensation point to the side M of the line based on a set distance measuring model according to the voltages Δ{dot over (U)}_φMxi and Δ{dot over (U)}_φNxi of the compensation point; and

at step 6, setting i=i+1; determining that a distance measuring result is a distance x_N+1 from the compensation point to the side M of the line in case of i>R; and turning back to the step 4 in case of i≤R; here R is a number of iterations.

The present disclosure further provides a line double-end steady-state quantity distance measuring system based on an amplitude-comparison principle. The system includes:

a data acquisition unit, configured to collect voltage values and current values, after an overhead line-cable hybrid transmission line is faulted, of both sides of the line, and voltage values and current values of both sides of the line of one cycle before the line is faulted; here both sides of the line are respectively a side M and a side N;

a first calculation unit, configured to determine voltage value variations according to the collected voltage values of both sides of the line before and after the line is faulted, and determine current value variations according to the collected current values of both sides of the line before and after the line is faulted;

a second calculation unit, configured to calculate voltage phasor values of both sides of the line by performing Fourier transformation on the voltage value variations of both sides of the line, and calculate current phasor values of both sides of the line by performing the Fourier transformation on the current value variations of both sides of the line;

a third calculation unit, configured to calculate voltages Δ{dot over (U)}_φMxi and Δ{dot over (U)}_φNxi of a compensation point according to a distance x_i from the compensation point to the side M of the line after the line is faulted, the voltage phasor values and the current phasor values of both sides of the line, a length L₁ of an overhead line on the side M, a length L₃ of an overhead line on the side N, a length L₂ of a cable, a wave impedance Z_cT and a propagation coefficient γ_T of the overhead line, and a wave impedance Z_cC and a current propagation coefficient γ_C of the cable; here an initial value of i is 1, and φ is any phase in a three-phase circuit, φ=A,B,C; and

a result determination unit, configured to determine a distance x_i+1 from the compensation point to the side M of the line based on a set distance measuring model according to the voltages Δ{dot over (U)}_φMxi and Δ{dot over (U)}_φNxi of the compensation point, set i=i+1, determine that a distance measuring result is a distance x_N+1 from the compensation point to the side M of the line in case of i>R, and turn back to the third calculation unit in case of i≤R; here R is a number of iterations.

BRIEF DESCRIPTION OF THE DRAWINGS

The exemplary implementation modes of the present disclosure can be understood more completely by referring to the following accompanying drawings:

FIG. 1 is a flow chart of a line double-end steady-state quantity distance measuring method based on an amplitude-comparison principle according to an optional implementation mode of the present disclosure.

FIG. 2 is a schematic diagram of an overhead line-cable hybrid line according to an optional implementation mode of the present disclosure.

FIG. 3 is a schematic structural diagram of a line double-end steady-state quantity distance measuring system based on an amplitude-comparison principle according to an optional implementation mode of the present disclosure.

DETAILED DESCRIPTION

The exemplary implementation modes of the present disclosure will now be described with reference to the accompanying drawings. However, the present disclosure can be implemented in many different forms and is not limited to the embodiments described here. These embodiments are provided to disclose the present disclosure in detail and completely, and fully convey the scope of the present disclosure to those skilled in the art. The terms in the exemplary implementation modes illustrated in the accompanying drawings do not limit the present disclosure. In the drawings, the same units/elements use the same reference signs.

Unless otherwise specified, the terms (including scientific and technological terms) used herein have the usual meanings to those skilled in the art. In addition, it can be understood that the terms defined in commonly used dictionaries should be understood as having consistent meanings in the context of their related art, and should not be understood as idealized or overly formal meanings.

FIG. 1 is a flow chart of a line double-end steady-state quantity distance measuring method based on an amplitude-comparison principle according to an optional implementation mode of the present disclosure. As illustrated in FIG. 1, the line double-end steady-state quantity distance measuring method 100 based on the amplitude-comparison principle of the optional implementation mode starts from step 101.

At step 101, distance measuring parameters are set; a wave impedance Z_cT and a propagation coefficient γ_T of an overhead line are determined; and a wave impedance Z_cC and a current propagation coefficient γ_C of a cable are determined.

At step 102, voltage values and current values, after an overhead line-cable hybrid transmission line is faulted, of both sides of the line, and voltage values and current values of both sides of the line of one cycle before the line is faulted are collected. Both sides of the line are respectively a side M and a side N.

FIG. 2 is a schematic diagram of an overhead line-cable hybrid line according to an optional implementation mode of the present disclosure. As illustrated in FIG. 2, the overhead line-cable hybrid line is divided into a total of three sections, an overhead line section 1 close to the side M of the line, a cable section, and an overhead line section 2 close to the side N of the line. A connection point of the overhead line section 1 and the cable is M₁, and a connection point of the overhead line section 2 and the cable is N₁.

At step 103, voltage value variations are determined according to the collected voltage values of both sides of the line before and after the line is faulted, and current value variations are determined according to the collected current values of both sides of the line before and after the line is faulted.

At step 104, voltage phasor values of both sides of the line are calculated by performing Fourier transformation on the voltage value variations of both sides of the line, and current phasor values of both sides of the line are calculated by performing the Fourier transformation on the current value variations of both sides of the line.

At step 105, voltages Δ{dot over (U)}_φMxi and Δ{dot over (U)}_φNxi of a compensation point are calculated according to a distance x_i from the compensation point to the side M of the line after the line is faulted, the voltage phasor values and the current phasor values of both sides of the line, a length L₁ of the overhead line on the side M, a length L₃ of the overhead line on the side N, a length L₂ of the cable, the wave impedance Z_cT and the propagation coefficient γ_T of the overhead line, and the wave impedance Z_cC and the current propagation coefficient γ_C of the cable. An initial value of i is 1, and φ is any phase in a three-phase circuit, φ=A,B,C.

At step 106, a distance x_i+1 from the compensation point to the side M of the line is determined based on a set distance measuring model according to the voltages Δ{dot over (U)}_φMxi and Δ{dot over (U)}_φNxi of the compensation point.

At step 107, i=i+1 is set; it is determined that a distance measuring result is a distance x_N+1 from the compensation point to the side M of the line in case of i>R; and the step 105 is executed in case of i≤R.

Optionally, before collecting the voltage values and the current values of both sides of the line after the overhead line-cable hybrid transmission line is faulted, the method further includes the following steps: the distance measuring parameters are set, the wave impedance Z_cT and the propagation coefficient γ_T of the overhead line are determined, and the wave impedance Z_cC and the current propagation coefficient γ_C of the cable are determined. The distance measuring parameters include a length L of the transmission line, the length L₁ of the overhead line on the side M, the length L₃ of the overhead line on the side N, the length L₂ of the cable, the number R of iterations, and an initial distance x₁ from the compensation point to the side M.

Calculation formulas for the wave impedance Z_cT and the propagation coefficient γ_T of the overhead line, and the wave impedance Z_cC and the current propagation coefficient γ_C of the cable are respectively:

$\begin{array}{c}{Z}_{\mathrm{cT}}=\sqrt{\frac{{𝓏}_T}{y_T}},\\ {\gamma }_T=\sqrt{{𝓏}_T{y}_T},\\ Z_{\mathrm{cC}}=\sqrt{\frac{{𝓏}_C}{y_C}},\end{array}$ γ_C=√{square root over (z_Cy_C)};

where z_T is a unit impedance of the overhead line; y_T is a unit admittance of the overhead line; z_C is a unit impedance of the cable; and y_C is a unit admittance of the cable.

Optionally, calculation formulas for calculating the voltage value variations according to the collected voltage values of both sides of the line before and after the line is faulted, and calculating the current value variations according to the collected current values of both sides of the line before and after the line is faulted are:

Δu_φM(t_qd+kT+t)=u_φM(t_qd+kT+t)−u_φM(t_qd−T+t),

Δu_φN(t_qd+kT+t)=u_φN(t_qd+kT+t)−u_φN(t_qd−T+t),

Δi_φM(t_qd+kT+t)=i_φM(t_qd+kT+t)−i_φM(t_qd−T+t),

Δi_φN(t_qd+kT+t)=i_φN(t_qd+kT+t)−i_φN(t_qd−T+t).

where both sides of the line are respectively the side M and the side N; t_qd is a distance measuring start moment; T is a power frequency cycle, 0<t<T, k≥1, φ=A,B,C; u_φM(t_qd+kT+t) and u_φN(t_qd+kT+t)are respectively the collected voltage values of the side M and the side N of the line of phase φ after the fault occurs; u_φM(t_qd−T+t) and u_φN(t_qd−T+t) are respectively the collected voltage values of the side M and the side N of the line of phase φ of one cycle before the fault occurs; u_φM(t_qd+kT+t) and u_φN(t_qd+kT+t) are respectively the voltage variations of the side M and the side N of the line of phase φ after the phase occurs; i_φM(t_qd+kT+t) and i_φN(t_qd+kT+t) are respectively the collected current values of the side M and the side N of the line of phase φ after the fault occurs; i_φM(t_qd−T+t) and i_φN(t_qd−kT+t) are respectively the collected current values of the side M and the side N of the line of phase φ of one cycle before the fault occurs; and Δi_φM(t_qd+kT+t) and Δi_φN(t_qd+kT+t) are respectively the current variations of the side M and the side N of the line of phase φ after the fault occurs.

Optionally, calculation formulas for calculating the voltages Δ{dot over (U)}_φMxi and ΔU_φNxi of the compensation point according to the distance x_i from the compensation point to the side M of the line after the line is faulted, the voltage phasor values and the current phasor values of both sides of the line, the length L₁ of the overhead line on the side M, the length L₃ of the overhead line on the side N, the length L₂ of the cable, the wave impedance Z_cT and the propagation coefficient γ_T of the overhead line, and the wave impedance Z_cC and the current propagation coefficient γ_C of the cable include as follows.

When the compensation point is at the overhead line section on the side M, the distance from the compensation point to the side M of the line is x_i, and the voltages Δ{dot over (U)}_φMxi and Δ{dot over (U)}_φNxi of the compensation point are calculated by using following calculation formulas:

$\begin{array}{c}\{\begin{array}{c}\Delta {\stackrel{.}{I}}_{\phi N1}=\Delta {\stackrel{.}{I}}_{\phi N}\cosh \left({\gamma }_T{L}_3\right)-\frac{\Delta {\stackrel{.}{U}}_{\phi N}}{Z_{\mathrm{cT}}}\sinh \left({\gamma }_T{L}_3\right)\\ \Delta {\stackrel{.}{U}}_{\phi N1}=\Delta {\stackrel{.}{U}}_{\phi N}\cosh \left({\gamma }_T{L}_3\right)-\Delta {\stackrel{.}{I}}_{\phi N}{Z}_{\mathrm{cT}}\sinh \left({\gamma }_T{L}_3\right)\end{array},\\ \{\begin{array}{c}\Delta {\stackrel{.}{I}}_{\phi M1}=\Delta {\stackrel{.}{I}}_{\phi N1}\cosh \left({\gamma }_C{L}_2\right)-\frac{\Delta {\stackrel{.}{U}}_{\phi N1}}{Z_{\mathrm{cC}}}\sinh \left({\gamma }_C{L}_2\right)\\ \Delta {\stackrel{.}{U}}_{\phi M1}=\Delta {\stackrel{.}{U}}_{\phi N1}\cosh \left({\gamma }_C{L}_2\right)-\Delta {\stackrel{.}{I}}_{\phi N1}{Z}_{\mathrm{cC}}\sinh \left({\gamma }_C{L}_2\right)\end{array},\\ \{\begin{array}{c}\Delta {\stackrel{.}{U}}_{\phi \mathrm{Mxi}}=\Delta {\stackrel{.}{U}}_{\phi M}\cosh \left({\gamma }_{T}x\right)-\Delta {\stackrel{.}{I}}_{\phi M}{Z}_{\mathrm{cT}}\sinh \left({\gamma }_{T}x\right)\\ \Delta {\stackrel{.}{U}}_{\phi \mathrm{Nxi}}=\Delta {\stackrel{.}{U}}_{\phi M1}\cosh \left[{\gamma }_T\left(L_1-x\right)\right]-\Delta {\stackrel{.}{I}}_{\phi M1}{Z}_{\mathrm{cT}}\sinh \left[{\gamma }_T\left(L_1-x\right)\right]\end{array}.\end{array}$

When the compensation point is at the cable section, the distance from the compensation point to the side M of the line is x_i, and the voltages Δ{dot over (U)}_φMxi and Δ{dot over (U)}_φNxi of the compensation point are calculated by using following calculation formulas:

$\begin{array}{c}\{\begin{array}{c}\Delta {\stackrel{.}{I}}_{\phi N1}=\Delta {\stackrel{.}{I}}_{\phi N}\cosh \left({\gamma }_T{L}_3\right)-\frac{\Delta {\stackrel{.}{U}}_{\phi N}}{Z_{\mathrm{cT}}}\sinh \left({\gamma }_T{L}_3\right)\\ \Delta {\stackrel{.}{U}}_{\phi N1}=\Delta {\stackrel{.}{U}}_{\phi N}\cosh \left({\gamma }_T{L}_3\right)-\Delta {\stackrel{.}{I}}_{\phi N}{Z}_{\mathrm{cT}}\sinh \left({\gamma }_T{L}_3\right)\end{array},\\ \{\begin{array}{c}\Delta {\stackrel{.}{I}}_{\phi M1}=\Delta {\stackrel{.}{I}}_{\phi M}\cosh \left({\gamma }_T{L}_1\right)-\frac{\Delta {\stackrel{.}{U}}_{\phi M}}{Z_{\mathrm{cT}}}\sinh \left({\gamma }_T{L}_1\right)\\ \Delta {\stackrel{.}{U}}_{\phi M1}=\Delta {\stackrel{.}{U}}_{\phi M}\cosh \left({\gamma }_T{L}_1\right)-\Delta {\stackrel{.}{I}}_{\phi M}{Z}_{\mathrm{cT}}\sinh \left({\gamma }_T{L}_1\right)\end{array},\end{array}$ $\{\begin{array}{c}\Delta {\stackrel{.}{U}}_{\phi \mathrm{Mxi}}=\Delta {\stackrel{.}{U}}_{\phi M1}\cosh \left[{\gamma }_C\left(x-L_1\right)\right]-\Delta {\stackrel{.}{I}}_{\phi M1}{Z}_{\mathrm{cC}}\sinh \left[{\gamma }_C\left(x-L_1\right)\right]\\ \Delta {\stackrel{.}{U}}_{\phi \mathrm{Nxi}}=\Delta {\stackrel{.}{U}}_{\phi N1}\cosh \left[{\gamma }_C\left(L_1+L_2-x\right)\right]-\Delta {\stackrel{.}{I}}_{\phi N1}{Z}_{\mathrm{cC}}\sinh \left[{\gamma }_C\left(L_1+L_2-x\right)\right]\end{array}.$

When the compensation point is at the overhead line section on the side N, the distance from the compensation point to the side M of the line is x_i, and the voltages Δ{dot over (U)}_φMxi and Δ{dot over (U)}_φNxi of the compensation point are calculated by using following calculation formulas:

$\begin{array}{c}\{\begin{array}{c}\Delta {\stackrel{.}{I}}_{\phi M1}=\Delta {\stackrel{.}{I}}_{\phi M}\cosh \left({\gamma }_T{L}_1\right)-\frac{\Delta {\stackrel{.}{U}}_{\phi M}}{Z_{\mathrm{cT}}}\sinh \left({\gamma }_T{L}_1\right)\\ \Delta {\stackrel{.}{U}}_{\phi M1}=\Delta {\stackrel{.}{U}}_{\phi M}\cosh \left({\gamma }_T{L}_1\right)-\Delta {\stackrel{.}{I}}_{\phi M}{Z}_{\mathrm{cT}}\sinh \left({\gamma }_T{L}_1\right)\end{array},\\ \{\begin{array}{c}\Delta {\stackrel{.}{I}}_{\phi N1}=\Delta {\stackrel{.}{I}}_{\phi M1}\cosh \left({\gamma }_C{L}_2\right)-\frac{\Delta {\stackrel{.}{U}}_{\phi M1}}{Z_{\mathrm{cC}}}\sinh \left({\gamma }_C{L}_2\right)\\ \Delta {\stackrel{.}{U}}_{\phi N1}=\Delta {\stackrel{.}{U}}_{\phi M1}\cosh \left({\gamma }_C{L}_2\right)-\Delta {\stackrel{.}{I}}_{\phi M1}{Z}_{\mathrm{cC}}\sinh \left({\gamma }_C{L}_2\right)\end{array},\end{array}$ $\{\begin{array}{c}\Delta {\stackrel{.}{U}}_{\phi \mathrm{Mxi}}=\Delta {\stackrel{.}{U}}_{\phi N1}\cosh \left[{\gamma }_T\left(x-L_1-L_2\right)\right]-\Delta {\stackrel{.}{I}}_{\phi N1}{Z}_{\mathrm{cT}}\sinh \left[{\gamma }_T\left(x-L_1-L_2\right)\right]\\ \Delta {\stackrel{.}{U}}_{\phi \mathrm{Nxi}}=\Delta {\stackrel{.}{U}}_{\phi N}\cosh \left[{\gamma }_T\left(L_1+L_2+L_3-x\right)\right]-\Delta {\stackrel{.}{I}}_{\phi N}{Z}_{\mathrm{cT}}\sinh \left[{\gamma }_T\left(L_1+L_2+L_3-x\right)\right]\end{array};$

where M₁ and N₁ are respectively connections between the overhead line of the side M of the line and the cable as well as between the side N of the line and the cable; Δ{dot over (U)}_φM1 and Δ{dot over (U)}_φN1 are respectively the voltage phasor values at the position M₁ and the position N₁ of the phase φ of the transmission line; Δ{dot over (U)}_φM and Δ{dot over (U)}_φN are respectively the voltage phasor values of the side M and the side N of the phase φ of the transmission line; Δİ_φM1 and Δİ_φN1 are respectively the current phasor values at the position M₁ and the position N₁ of the phase φ of the transmission line; Δİ_φM and Δİ_φN are respectively the current phasor values of the side M and the side N of the phase φ of the transmission line; Δ{dot over (U)}_φMxi is the voltage of the compensation point calculated and determined according to the voltage phasor value and the current phasor value close to the side M of the compensation point; and Δ{dot over (U)}_φNxi is the voltage of the compensation point calculated and determined according to the voltage phasor value and the current phasor value close to the side N of the compensation point.

Optionally, in order to determine the distance x_i+1 from the compensation point to the side M of the line based on the set distance measuring model according to the voltages Δ{dot over (U)}_φMxi and Δ{dot over (U)}_φNxi of the compensation point, a calculation formula of the distance measuring model is:

$x_{i+1}=\{\begin{array}{cc}{x}_i-\frac{L}{2^{i+1}}& ❘\Delta {\stackrel{.}{U}}_{\phi \mathrm{Mxi}}❘\ge ❘\Delta {\stackrel{.}{U}}_{\phi \mathrm{Nxi}}❘\\ x_i+\frac{L}{2^{i+1}}& ❘\Delta {\stackrel{.}{U}}_{\phi \mathrm{Mxi}}❘<❘\Delta {\stackrel{.}{U}}_{\phi \mathrm{Nxi}}❘\end{array};$

where L is a length of the transmission line.

A simulation system is constructed by taking the line of FIG. 2 as an example Line fault points F1, F2, F3, F4, and F5 are away from the side M by 0 km, 2.3 km, 10.8 km, 19.3 km and 51.35 km, respectively. Phase-A metallic fault simulation and phase-AB metallic grounding fault simulation are carried out on different positions of the line. The phase A is simulated via a 100Ω transition resistor and distance measuring results and actual fault positions are as illustrated in Table 1.

TABLE 1
Comparison of distance measuring results and actual fault positions of different fault types
Fault point	F1	F2	F3	F4	F5
Actual position (km)	0	2.3	10.8	19.3	51.35
Phase-A metallic grounding distance	0.050	2.257	10.581	19.106	51.300
measuring result (km)
Phase-AB metallic grounding distance	0.050	2.257	10.681	19.306	51.200
measuring result (km)
Grounding distance measuring result of	0.050	2.257	10.781	19.306	51.300
phase A via the 100Ω transition resistor
(km)

It can be known from Table 1 that the distance measuring results of the fault points calculated according to the method of the present disclosure are close to actual results.

FIG. 3 is a schematic structural diagram of a line double-end steady-state quantity distance measuring system based on an amplitude-comparison principle according to an optional implementation mode of the present disclosure. As illustrated in FIG. 3, a line double-end steady-state quantity distance measuring system 300 based on an amplitude-comparison principle of the optional implementation mode may include: an initialization unit 301, a data acquisition unit 302, a first calculation unit 303, a second calculation unit 304, a third calculation unit 305 and a result determination unit 306.

The initialization unit 301 is configured to set distance measuring parameters, determine a wave impedance Z_cT and a propagation coefficient γ_T of an overhead line, and determine a wave impedance Z_cC and a current propagation coefficient γ_C of a cable.

The data acquisition unit 302 is configured to collect voltage values and current values, after an overhead line-cable hybrid transmission line is faulted, of both sides of the line, and voltage values and current values of both sides of the line of one cycle before the line is faulted, where both sides of the line are respectively a side M and a side N.

The first calculation unit 303 is configured to determine voltage value variations according to the collected voltage values of both sides of the line before and after the line is faulted, and determine current value variations according to the collected current values of both sides of the line before and after the line is faulted.

The second calculation unit 304 is configured to calculate voltage phasor values of both sides of the line by performing Fourier transformation on the voltage value variations of both sides of the line, and calculate current phasor values of both sides of the line by performing the Fourier transformation on the current value variations of both sides of the line.

The third calculation unit 305 is configured to calculate voltages Δ{dot over (U)}_φMxi and Δ{dot over (U)}_φNxi of a compensation point according to a distance x_i from the compensation point to the side M of the line after the line is faulted, the voltage phasor values and the current phasor values of both sides of the line, a length L₁ of the overhead line on the side M, a length L₃ of the overhead line on the side N, a length L₂ of the cable, the wave impedance Z_cT and the propagation coefficient γ_T of the overhead line, and the wave impedance Z_cC and the current propagation coefficient γ_C of the cable, where an initial value of i is 1, and φ is any phase in a three-phase circuit, φ=A,B,C.

The result determination unit 306 is configured to determine a distance x_i+1 from the compensation point to the side M of the line based on a set distance measuring model according to the voltages Δ{dot over (U)}_φMxi and Δ{dot over (U)}_φNxi of the compensation point, set i=i+1, determine that a distance measuring result is a distance x_N+1 from the compensation point to the side M of the line in case of i>R, and turn back to the third calculation unit 305 in case of i≤R.

Optionally, the initialization unit 301 is configured to set the distance measuring parameters, determine the wave impedance Z_cT and the propagation coefficient γ_T of the overhead line, and determine the wave impedance Z_cC and the current propagation coefficient γ_C of the cable. The distance measuring parameters include a length L of the transmission line, the length L₁ of the overhead line on the side M, the length L₃ of the overhead line on the side N, the length L₂ of the cable, the number R of iterations, and an initial distance x₁ from the compensation point to the side M. Calculation formulas for the wave impedance Z_cT and the propagation coefficient γ_T of the overhead line, and the wave impedance Z_cC and the current propagation coefficient γ_C of the cable are respectively:

$\begin{array}{c}{Z}_{\mathrm{cT}}=\sqrt{\frac{{𝓏}_T}{y_T}},\\ {\gamma }_T=\sqrt{{𝓏}_T{y}_T},\\ Z_{\mathrm{cC}}=\sqrt{\frac{{𝓏}_C}{y_C}},\end{array}$ γ_C=√{square root over (z_Cy_C)};

where z_T is a unit impedance of the overhead line; y_T is a unit admittance of the overhead line; z_C is a unit impedance of the cable; and y_C is a unit admittance of the cable.

Optionally, the first calculation unit 303 is configured to determine the voltage value variations according to the collected voltage values of both sides of the line before and after the line is faulted, and determine the current value variations according to the collected current values of both sides of the line before and after the line is faulted by using following calculation formulas:

Δu_φM(t_qd+kT+t)=u_φM(t_qd+kT+t)−u_φM(t_qd−T+t),

Δu_φN(t_qd+kT+t)=u_φN(t_qd+kT+t)−u_φN(t_qd−T+t),

Δi_φM(t_qd+kT+t)=i_φM(t_qd+kT+t)−i_φM(t_qd−T+t),

Δi_φN(t_qd+kT+t)=i_φN(t_qd+kT+t)−i_φN(t_qd−T+t).

where both sides of the line are respectively the side M and the side N; t_qd is a distance measuring start moment; T is a power frequency cycle, 0<t<T, k≥1, φ=A,B,C; u_φM(t_qd+kT+t) and u_φN(t_qd+kT+t) are respectively the collected voltage values of the side M and the side N of the line of phase φ after the fault occurs; u_φM(t_qd−T+t) and u_φN(t_qd−T+t) are respectively the collected voltage values of the side M and the side N of the line of phase φ of one cycle before the fault occurs; Δu_φM(t_qd+kT+t) and Δu_φN(t_qd+kT+t) are respectively the voltage variations of the side M and the side N of the line of phase φ after the phase occurs; i_φM(t_qd+kT+t) and i_φN(t_qd+kT+t) are respectively the collected current values of the side M and the side N of the line of phase φ after the fault occurs; i_φM(t_qd−T+t) and i_φN(t_qd−T+t) are respectively the collected current values of the side M and the side N of the line of phase φ of one cycle before the fault occurs; and Δi_φM(t_qd+kT+t) and Δi_φN(t_qd+kT+t) are respectively the current variations of the side M and the side N of the line of phase φ after the fault occurs.

Optionally, the third calculation unit 305 is configured to calculate the voltages Δ{dot over (U)}_φMxi and Δ{dot over (U)}_φNxi of the compensation point according to the distance x_i from the compensation point to the side M of the line after the line is faulted, the voltage phasor values and the current phasor values of both sides of the line, the length L₁ of the overhead line on the side M, the length L₃ of the overhead line on the side N, the length L₂ of the cable, the wave impedance Z_cT and the propagation coefficient γ_T of the overhead line, and the wave impedance Z_cC and the current propagation coefficient γ_C of the cable by using following calculation formulas.

When the compensation point is at the overhead line section on the side M, the distance from the compensation point to the side M of the line is x_i, and the voltages Δ{dot over (U)}_φMxi and Δ{dot over (U)}_φNxi of the compensation point are calculated by using following calculation formulas:

$\begin{array}{c}\{\begin{array}{c}\Delta {\stackrel{.}{I}}_{\phi N1}=\Delta {\stackrel{.}{I}}_{\phi N}\cosh \left({\gamma }_T{L}_3\right)-\frac{\Delta {\stackrel{.}{U}}_{\phi N}}{Z_{\mathrm{cT}}}\sinh \left({\gamma }_T{L}_3\right)\\ \Delta {\stackrel{.}{U}}_{\phi N1}=\Delta {\stackrel{.}{U}}_{\phi N}\cosh \left({\gamma }_T{L}_3\right)-\Delta {\stackrel{.}{I}}_{\phi N}{Z}_{\mathrm{cT}}\sinh \left({\gamma }_T{L}_3\right)\end{array},\\ \{\begin{array}{c}\Delta {\stackrel{.}{I}}_{\phi M1}=\Delta {\stackrel{.}{I}}_{\phi N1}\cosh \left({\gamma }_C{L}_2\right)-\frac{\Delta {\stackrel{.}{U}}_{\phi N1}}{Z_{\mathrm{cC}}}\sinh \left({\gamma }_C{L}_2\right)\\ \Delta {\stackrel{.}{U}}_{\phi M1}=\Delta {\stackrel{.}{U}}_{\phi N1}\cosh \left({\gamma }_C{L}_2\right)-\Delta {\stackrel{.}{I}}_{\phi N1}{Z}_{\mathrm{cC}}\sinh \left({\gamma }_C{L}_2\right)\end{array},\end{array}$ $\{\begin{array}{c}\Delta {\stackrel{.}{U}}_{\phi \mathrm{Mxi}}=\Delta {\stackrel{.}{U}}_{\phi M}\cosh \left({\gamma }_{T}x\right)-\Delta {\stackrel{.}{I}}_{\phi M}{Z}_{\mathrm{cT}}\sinh \left({\gamma }_{T}x\right)\\ \Delta {\stackrel{.}{U}}_{\phi \mathrm{Nxi}}=\Delta {\stackrel{.}{U}}_{\phi M1}\cosh \left[{\gamma }_T\left(L_1-x\right)\right]-\Delta {\stackrel{.}{I}}_{\phi M1}{Z}_{\mathrm{cT}}\sinh \left[{\gamma }_T\left(L_1-x\right)\right]\end{array}.$

When the compensation point is at the cable section, the distance from the compensation point to the side M of the line is x_i, and the voltages Δ{dot over (U)}_φMxi and Δ{dot over (U)}_φNxi of the compensation point are calculated by using following calculation formulas:

$\begin{array}{c}\{\begin{array}{c}\Delta {\stackrel{.}{I}}_{\phi N1}=\Delta {\stackrel{.}{I}}_{\phi N}\cosh \left({\gamma }_T{L}_3\right)-\frac{\Delta {\stackrel{.}{U}}_{\phi N}}{Z_{\mathrm{cT}}}\sinh \left({\gamma }_T{L}_3\right)\\ \Delta {\stackrel{.}{U}}_{\phi N1}=\Delta {\stackrel{.}{U}}_{\phi N}\cosh \left({\gamma }_T{L}_3\right)-\Delta {\stackrel{.}{I}}_{\phi N}{Z}_{\mathrm{cT}}\sinh \left({\gamma }_T{L}_3\right)\end{array},\\ \{\begin{array}{c}\Delta {\stackrel{.}{I}}_{\phi M1}=\Delta {\stackrel{.}{I}}_{\phi M}\cosh \left({\gamma }_T{L}_1\right)-\frac{\Delta {\stackrel{.}{U}}_{\phi M}}{Z_{\mathrm{cT}}}\sinh \left({\gamma }_T{L}_1\right)\\ \Delta {\stackrel{.}{U}}_{\phi M1}=\Delta {\stackrel{.}{U}}_{\phi M}\cosh \left({\gamma }_T{L}_1\right)-\Delta {\stackrel{.}{I}}_{\phi M}{Z}_{\mathrm{cT}}\sinh \left({\gamma }_T{L}_1\right)\end{array},\end{array}$

$\{\begin{array}{c}\Delta {\stackrel{.}{U}}_{\phi \mathrm{Mxi}}=\Delta {\stackrel{.}{U}}_{\phi M1}\cos h\left[{\gamma }_C\left(x-L_1\right)\right]-\Delta {\stackrel{.}{I}}_{\phi M1}{Z}_{\mathrm{cC}}\sin h\left[{\gamma }_C\left(x-L_1\right)\right]\\ \Delta {\stackrel{.}{U}}_{\phi \mathrm{Nxi}}=\Delta {\stackrel{.}{U}}_{\phi N1}\cos h\left[{\gamma }_C\left(L_1+L_2-x\right)\right]-\Delta {\stackrel{.}{I}}_{\phi N1}{Z}_{\mathrm{cC}}\sin h\left[{\gamma }_C\left(L_1+L_2-x\right)\right]\end{array}.$

When the compensation point is at the overhead line section on the side N, the distance from the compensation point to the side M of the line is x_i, and the voltages Δ{dot over (U)}_φMxi and Δ{dot over (U)}_φNxi of the compensation point are calculated:

$\begin{array}{c}\{\begin{array}{c}\Delta {\stackrel{.}{I}}_{\phi M1}=\Delta {\stackrel{.}{I}}_{\phi M}\cosh \left({\gamma }_T{L}_1\right)-\frac{\Delta {\stackrel{.}{U}}_{\phi M}}{Z_{\mathrm{cT}}}\sinh \left({\gamma }_T{L}_1\right)\\ \Delta {\stackrel{.}{U}}_{\phi M1}=\Delta {\stackrel{.}{U}}_{\phi M}\cosh \left({\gamma }_T{L}_1\right)-\Delta {\stackrel{.}{I}}_{\phi M}{Z}_{\mathrm{cT}}\sinh \left({\gamma }_T{L}_1\right)\end{array},\\ \{\begin{array}{c}\Delta {\stackrel{.}{I}}_{\phi N1}=\Delta {\stackrel{.}{I}}_{\phi M1}\cosh \left({\gamma }_C{L}_2\right)-\frac{\Delta {\stackrel{.}{U}}_{\phi M1}}{Z_{\mathrm{cC}}}\sinh \left({\gamma }_C{L}_2\right)\\ \Delta {\stackrel{.}{U}}_{\phi N1}=\Delta {\stackrel{.}{U}}_{\phi M1}\cosh \left({\gamma }_C{L}_2\right)-\Delta {\stackrel{.}{I}}_{\phi M1}{Z}_{\mathrm{cC}}\sinh \left({\gamma }_C{L}_2\right)\end{array},\end{array}$ $\{\begin{array}{c}\Delta {\stackrel{.}{U}}_{\phi \mathrm{Mxi}}=\Delta {\stackrel{.}{U}}_{\phi N1}\cosh \left[{\gamma }_T\left(x-L_1-L_2\right)\right]-\Delta {\stackrel{.}{I}}_{\phi N1}{Z}_{\mathrm{cT}}\sinh \left[{\gamma }_T\left(x-L_1-L_2\right)\right]\\ \Delta {\stackrel{.}{U}}_{\phi \mathrm{Nxi}}=\Delta {\stackrel{.}{U}}_{\phi N}\cosh \left[{\gamma }_T\left(L_1+L_2+L_3-x\right)\right]-\Delta {\stackrel{.}{I}}_{\phi N}{Z}_{\mathrm{cT}}\sinh \left[{\gamma }_T\left(L_1+L_2+L_3-x\right)\right]\end{array};$

where M₁ and N₁ are respectively connections between the overhead line of the side M of the line and the cable as well as between the side N of the line and the cable; Δ{dot over (U)}_φM1 and Δ{dot over (U)}_φN1 are respectively the voltage phasor values at the position M₁ and the position N₁ of the phase φ of the transmission line; ΔU_φM and Δ{dot over (U)}_φN are respectively the voltage phasor values of the side M and the side N of the phase φ of the transmission line; Δİ_φM1 and Δİ_φN1 are respectively the current phasor values at the position M₁ and the position N₁ of the phase φ of the transmission line; Δİ_φM and Δİ_φN are respectively the current phasor values of the side M and the side N of the phase φ of the transmission line; Δ{dot over (U)}_φMxi is the voltage of the compensation point calculated and determined according to the voltage phasor value and the current phasor value close to the side M of the compensation point; and Δ{dot over (U)}_φNxi is the voltage of the compensation point calculated and determined according to the voltage phasor value and the current phasor value close to the side N of the compensation point.

Optionally, the result determination unit 306 is configured to determine the distance x_i+1 from the compensation point to the side M of the line based on the set distance measuring model according to the voltages Δ{dot over (U)}_φMxi and Δ{dot over (U)}_φNxi of the compensation point by using a following calculation formula:

$x_{i+1}\{\begin{array}{cc}{x}_i-\frac{L}{2^{i+1}}& ❘\Delta {\stackrel{.}{U}}_{\phi \mathrm{Mxi}}❘\ge ❘\Delta {\stackrel{.}{U}}_{\phi \mathrm{Nxi}}❘\\ x_i+\frac{L}{2^{i+1}}& ❘\Delta {\stackrel{.}{U}}_{\phi \mathrm{Mxi}}❘<❘\Delta {\stackrel{.}{U}}_{\phi \mathrm{Nxi}}❘\end{array};$

where L is a length of the transmission line.

The distance measuring steps of the line double-end steady-state quantity distance measuring system based on the amplitude-comparison principle for the overhead line-cable hybrid line are the same as the steps used in the line double-end steady-state quantity distance measuring method based on the amplitude-comparison principle, and technical effects achieved are also the same, so descriptions thereof are omitted.

The present disclosure has been described by the above implementation modes. However, as is well known to those skilled in the art, as defined by the appended patent claims, other embodiments than those disclosed above in the present disclosure equally fall within the scope of the present disclosure.

Generally, all terms used in the claims are interpreted according to their ordinary meanings in the technical art, unless explicitly defined otherwise therein. All references to “a/the/the [apparatus, assembly, etc.]” are openly interpreted as at least one example of the apparatus, component, etc., unless explicitly stated otherwise. The steps of any method disclosed herein need not be run in the exact order disclosed, unless explicitly stated.

Those skilled in the art should understand that the embodiments of the present disclosure may be methods, systems or computer program products. Therefore, the present disclosure may adopt the form of a complete hardware embodiment, a complete software embodiment, or a software and hardware combination embodiment. In addition, the present disclosure may use the form of a computer program product implemented on one or multiple computer-sensitive storage media (including a magnetic disk memory, a compact disc read-only memory (CD-ROM), an optical memory and the like) including computer-sensitive program codes.

The present disclosure is described by referring to flowcharts and/or block diagrams of methods, devices (systems) and computer program products according to the embodiments of the present disclosure. It is worthwhile to note that computer program instructions can be used to implement each process and/or each block in the flowcharts and/or the block diagrams and a combination of a process and/or a block in the flowcharts and/or the block diagrams. These computer program instructions can be provided for a general-purpose computer, a dedicated computer, an embedded processor, or a processor of another programmable data processing device to generate a machine, so that the instructions executed by the computer or the processor of the another programmable data processing device generate an apparatus for implementing a specified function in one or more processes in the flowcharts and/or in one or more blocks in the block diagrams.

Alternatively, these computer program instructions can be stored in a computer-readable memory that can instruct the computer or the another programmable data processing device to work in a specific way, so that the instructions stored in the computer-readable memory generate an artifact that includes an instruction apparatus. The instruction apparatus implements a specified function in one or more processes in the flowcharts and/or in one or more blocks in the block diagrams.

Alternatively, these computer program instructions can be loaded onto the computer or the another programmable data processing device, so that a series of operations and steps are performed on the computer or the another programmable device, thereby generating computer-implemented processing. Therefore, the instructions executed on the computer or the another programmable device provide steps for implementing a specified function in one or more processes in the flowcharts and/or in one or more blocks in the block diagrams.

The line double-end steady-state quantity distance measuring method and system based on the amplitude-comparison principle provided by the technical solutions of the present disclosure determine the position of a short-circuit point by means of iteratively calculating voltages of the short-circuit point after collecting the voltage values and the current values of both sides of the line before and after the fault occurs, calculating the voltage variations and the current variations of both sides of the line, and determining the voltage phasor values and the current phasor values according to the voltage variations and the current variations. The method has a simple principle and can accurately identify a fault point and realize accurate distance measuring of the line.

Claims

1. A line double-end steady-state quantity distance measuring method based on an amplitude-comparison principle, comprising: at step 1, collecting voltage values and current values, after an overhead line-cable hybrid transmission line is faulted, of both sides of the line, and voltage values and current values of both sides of the line of one cycle before the line is faulted; wherein both sides of the line are respectively a side M and a side N;at step 2, determining voltage value variations according to the collected voltage values of both sides of the line before and after the line is faulted, and determining current value variations according to the collected current values of both sides of the line before and after the line is faulted;at step 3, calculating voltage phasor values of both sides of the line by performing Fourier transformation on the voltage value variations of both sides of the line, and calculating current phasor values of both sides of the line by performing the Fourier transformation on the current value variations of both sides of the line;at step 4, calculating voltages Δ{dot over (U)}φMxi and Δ{dot over (U)}φNxi of a compensation point according to a distance xi from the compensation point to the side M of the line after the line is faulted, the voltage phasor values and the current phasor values of both sides of the line, a length L1 of an overhead line on the side M, a length L3 of an overhead line on the side N, a length L2 of a cable, a wave impedance ZcT and a propagation coefficient γT of the overhead line, and a wave impedance ZcC and a current propagation coefficient γC of the cable; wherein an initial value of i is 1, and φ is any phase in a three-phase circuit, φ=A,B,C;at step 5, determining a distance xi+1 from the compensation point to the side M of the line based on a set distance measuring model according to the voltages Δ{dot over (U)}φMxi and Δ{dot over (U)}φNxi of the compensation point; andat step 6, setting i=i+1; determining that a distance measuring result is a distance xN+1 from the compensation point to the side M of the line in case of i>R; and turning back to the step 4 in case of i≤R; wherein R is a number of iterations.

2. The method of claim 1, wherein before collecting the voltage values and the current values of both sides of the line after the overhead line-cable hybrid transmission line is faulted, further comprising: setting distance measuring parameters, determining the wave impedance ZcT and the propagation coefficient γT of the overhead line, and determining the wave impedance ZcC and the current propagation coefficient γC of the cable; wherein the distance measuring parameters comprise a length L of the transmission line, the length L1 of the overhead line on the side M, the length L3 of the overhead line on the side N, the length L2 of the cable, the number R of iterations, and an initial distance x1 from the compensation point to the side M;wherein calculation formulas for the wave impedance ZcT and the propagation coefficient γT of the overhead line, and the wave impedance ZcC and the current propagation coefficient γC of the cable are respectively: ZcT=𝓏TyT,γT=𝓏T⁢yT,ZcC=𝓏CyC,γC=√{square root over (zCyC)};wherein zT is a unit impedance of the overhead line; yT is a unit admittance of the overhead line; zC is a unit impedance of the cable; and yC is a unit admittance of the cable.

3. The method of claim 1, wherein calculation formulas for calculating the voltage value variations according to the collected voltage values of both sides of the line before and after the line is faulted, and calculating the current value variations according to the collected current values of both sides of the line before and after the line is faulted are: ΔuφM(tqd+kT+t)=uφM(tqd+kT+t)−uφM(tqd−T+t),ΔuφN(tqd+kT+t)=uφN(tqd+kT+t)−uφN(tqd−T+t),ΔiφM(tqd+kT+t)=iφM(tqd+kT+t)−iφM(tqd−T+t),ΔiφN(tqd+kT+t)=iφN(tqd+kT+t)−iφN(tqd−T+t);wherein both sides of the line are respectively the side M and the side N; tqd is a distance measuring start moment; T is a power frequency cycle, 0<t<T , k≥1, φ=A,B,C; uφM(tqd+kT+t) and uφN(tqd+kT+t) are respectively the collected voltage values of the side M and the side N of the line of phase φ after the fault occurs; uφM(tqd−T+t) and uφN(tqd−T+t) are respectively the collected voltage values of the side M and the side N of the line of phase φ of one cycle before the fault occurs; ΔuφM(tqd+kT+t) and uΔφN(tqd+kT+t) are respectively the voltage variations of the side M and the side N of the line of phase φ after the phase occurs; iφM(tqd+kT+t) and iφN(tqd+kT+t) are respectively the collected current values of the side M and the side N of the line of phase φ after the fault occurs; iφM(tqd−T+t) and iφN(tqd−T+t) are respectively the collected current values of the side M and the side N of the line of phase φ of one cycle before the fault occurs; and ΔiφM(tqd+kT+t) and ΔiφN(tqd+kT+t) are respectively the current variations of the side M and the side N of the line of phase φ after the fault occurs.

4. The method of claim 2, wherein calculation formulas for calculating the voltages Δ{dot over (U)}φMxi and Δ{dot over (U)}φNxi of the compensation point according to the distance xi from the compensation point to the side M of the line after the line is faulted, the voltage phasor values and the current phasor values of both sides of the line, the length L1 of the overhead line on the side M, the length L3 of the overhead line on the side N, the length L2 of the cable, the wave impedance ZcT and the propagation coefficient γT of the overhead line, and the wave impedance ZcC and the current propagation coefficient γC of the cable comprises: when the compensation point is at an overhead line section on the side M, the distance from the compensation point to the side M of the line is xi, and the voltages Δ{dot over (U)}φMxi and Δ{dot over (U)}φNxi of the compensation point are calculated by using following calculation formulas:

5. The method of claim 2, wherein during determining the distance xi+1 from the compensation point to the side M of the line based on the set distance measuring model according to the voltages Δ{dot over (U)}φMxi and Δ{dot over (U)}φNxi of the compensation point, a calculation formula of the distance measuring model is:

6. A line double-end steady-state quantity distance measuring system based on an amplitude-comparison principle, comprising: a data acquisition component, configured to collect voltage values and current values, after an overhead line-cable hybrid transmission line is faulted, of both sides of the line, and voltage values and current values of both sides of the line of one cycle before the line is faulted; wherein both sides of the line are respectively a side M and a side N;a first calculator, configured to determine voltage value variations according to the collected voltage values of both sides of the line before and after the line is faulted, and determine current value variations according to the collected current values of both sides of the line before and after the line is faulted;a second calculator, configured to calculate voltage phasor values of both sides of the line by performing Fourier transformation on the voltage value variations of both sides of the line, and calculate current phasor values of both sides of the line by performing the Fourier transformation on the current value variations of both sides of the line;a third calculator, configured to calculate voltages Δ{dot over (U)}φMxi and Δ{dot over (U)}φNxi of a compensation point according to a distance xi from the compensation point to the side M of the line after the line is faulted, the voltage phasor values and the current phasor values of both sides of the line, a length L1 of an overhead line on the side M, a length L3 of an overhead line on the side N, a length L2 of a cable, a wave impedance ZcT and a propagation coefficient γT of the overhead line, and a wave impedance ZcC and a current propagation coefficient γC of the cable; wherein an initial value of i is 1, and φ is any phase in a three-phase circuit, φ=A,B,C; anda result determination component, configured to determine a distance xi+1 from the compensation point to the side M of the line based on a set distance measuring model according to the voltages Δ{dot over (U)}φMxi and Δ{dot over (U)}φNxi of the compensation point, set i=i+1, determine that a distance measuring result is a distance xN+1 from the compensation point to the side M of the line in case of i>R, and turn back to the third calculator in case of i≤R; wherein R is a number of iterations.

7. The system of claim 6, further comprising an initialization component; wherein the initialization component is configured to set distance measuring parameters, determine the wave impedance ZcT and the propagation coefficient γT of the overhead line, and determine the wave impedance ZcC an the current propagation coefficient γC of the cable; wherein the distance measuring parameters comprises a length L of the transmission line, the length L1 of the overhead line on the side M, the length L3 of the overhead line on the side N, the length L2 of the cable, the number R of iterations, and an initial distance x1 from the compensation point to the side M; wherein calculation formulas for the wave impedance ZcT and the propagation coefficient γT of the overhead line, and the wave impedance ZcC and the current propagation coefficient γC of the cable are respectively: ZcT=𝓏TyT,γT=𝓏T⁢yT,ZcC=𝓏CyC,γC=√{square root over (zCyC)};wherein zT is a unit impedance of the overhead line; yT is a unit admittance of the overhead line; zC is a unit impedance of the cable; and yC is a unit admittance of the cable.

8. The system of claim 6, wherein the first calculator is configured to use calculation formulas: ΔuφM(tqd+kT+t)=uφM(tqd+kT+t)−uφM(tqd−T+t),ΔuφN(tqd+kT+t)=uφN(tqd+kT+t)−uφN(tqd−T+t),ΔiφM(tqd+kT+t)=iφM(tqd+kT+t)−iφM(tqd−T+t),ΔiφN(tqd+kT+t)=iφN(tqd+kT+t)−iφN(tqd−T+t);wherein both sides of the line are respectively the side M and the side N; tqd is a distance measuring start moment; T is a power frequency cycle, 0<t<T, k≥1, φ=A,B,C; uφM(tqd+kT+t) and uφN(tqd+kT+t) are respectively the collected voltage values of the side M and the side N of the line of phase φ after the fault occurs; uφM(tqd−T+t) and uφN(tqd−T+t) are respectively the collected voltage values of the side M and the side N of the line of phase φ of one cycle before the fault occurs; ΔuφM(tqd+kT+t) and ΔuφN(tqd+kT+t) are respectively the voltage variations of the side M and the side N of the line of phase φ after the phase occurs; iφM(tqd+kT+t) and iφN(tqd+kT+t) are respectively the collected current values of the side M and the side N of the line of phase φ after the fault occurs; iφM(tqd−T+t) and iφN(tqd−T+t) are respectively the collected current values of the side M and the side N of the line of phase φ of one cycle before the fault occurs; and ΔiφM(tqd+kT+t) and ΔiφN(tqd+kT+t) are respectively the current variations of the side M and the side N of the line of phase φ after the fault occurs.

9. The system of claim 7, wherein the third calculator is configured to use calculation formulas: when the compensation point is at the overhead line section on the side M, the distance from the compensation point to the side M of the line is xi, and the voltages Δ{dot over (U)}φMxi and Δ{dot over (U)}φNxi of the compensation point are calculated by using following calculation formulas:

10. The system of claim 7, wherein a calculation formula, adopted by the result determination component, for the distance measuring model is: