A fault location method and device based on berlincourt differential flow and a storage medium

Abstract

The application discloses a fault location method based on Bergeron differential flow, comprising the following steps: acquiring voltage sampling values and current sampling values of M ends and N ends of a line; performing phase-mode transformation on the voltage sampling values and current sampling values of the two ends to obtain mode voltage and current of the two ends; calculating the mode Bergeron differential flow of the two ends according to the mode voltage and current of the two ends, wherein the Bergeron differential flow is the difference between an actual value of the current end and a calculated value of the current end obtained by using a Bergeron model of the opposite end; calculating the time t m and t n when the mode Bergeron differential flow of the two ends is equal to a preset current threshold value respectively; calculating the difference between the time t m and t n , and combining the traveling wave transmission speed to calculate the distance from the fault point to the line end. The method has clear principles, simple algorithm, does not need to identify the traveling wave head, and is easy to implement in engineering.

Description

TECHNICAL FIELD

[0001] The application belongs to the field of power system relay protection, and particularly relates to a fault distance measurement method and device based on Bergeron differential current and a storage medium. BACKGROUND

[0002] High-voltage transmission lines are the lifelines of power systems, and bear the heavy responsibility of transmitting electric energy. At the same time, it is the place where faults occur most frequently in the system, and it is extremely difficult to find. Therefore, quickly and accurately finding the fault point after a fault not only is important for timely repairing the line and ensuring power supply reliability, but also is important for the safe and stable and economic operation of the power system.

[0003] Fault distance measurement methods are divided into fault analysis method and traveling wave method in principle, and can be divided into single-end method and double-end method according to the source of distance measurement information. The fault analysis method is a general method for calculating the distance between the fault point and the distance measurement point according to the system parameters and the voltage and current series of the distance measurement point when a fault occurs on the transmission line, and then analyzing and calculating the distance measurement equation to obtain the distance between the fault point and the distance measurement point. It is not only suitable for single-end distance measurement, but also suitable for double-end distance measurement, and the difference between the two is that the former uses the voltage and current information at one end of the transmission line, and the latter uses the voltage and current information at both ends of the line. The traveling wave method is a distance measurement method based on the traveling wave theory. The traveling wave method is also divided into single-end method and double-end method. The former is realized by using the transient traveling wave quantity detected at one end of the line, and the latter is a fault distance measurement method realized by using the quantities at both ends of the line. The fault analysis method and the traveling wave method based on single-end electrical quantities are limited in distance measurement accuracy due to the limited fault information and uncertain system structure. The distance measurement method based on double-end electrical quantities has more comprehensive fault information and can obtain higher distance measurement accuracy. The traveling wave distance measurement method needs to identify the traveling wave front, and it is difficult to accurately identify the traveling wave front after the traveling wave front decays through the line. The current identification method is relatively complex, and it is necessary to improve the traveling wave method to make the algorithm simpler and easier to implement in engineering. SUMMARY

[0004] The purpose of the present application is to provide a fault distance measurement method and device based on Bergeron differential current and a storage medium, which calculates the time difference of the Bergeron differential current at both ends of the line reaching the same threshold, determines the time difference of the fault traveling wave reaching both ends of the line, and determines the fault point position in combination with the wave speed and the line length. The distance measurement method has clear principles, simple algorithm, does not need to identify the traveling wave front, and is easy to implement in engineering.

[0005] In order to achieve the above purpose, the solution of the present application is:

[0006] The first aspect of the present application provides a fault distance measurement method based on Bergeron differential current, comprising the following steps:

[0007] Obtaining voltage sampling values and current sampling values at M ends and N ends at both ends of the line;

[0008] performing phase-to-magnitude transformation on the voltage sampling value and the current sampling value of the two ends to obtain a magnitude voltage and a magnitude current of the two ends respectively;

[0009] calculating a magnitude Bergeron difference flow of the two ends according to the magnitude voltage and the magnitude current of the two ends, the Bergeron difference flow being a difference between an actual value of the two ends and a calculated value of the two ends calculated by using a Bergeron model of a counter end;

[0010] calculating times t m and t n at which the magnitude Bergeron difference flow of the two ends is equal to a preset current threshold value respectively; m and t n ;

[0011] calculating a difference value of the times t m and t n , and calculating a distance from a fault point to a line end according to a traveling wave transmission speed.

[0012] Preferably, the phase-to-magnitude transformation uses a 3-order conversion matrix.

[0013] Preferably, the magnitude voltage and the magnitude current are any one of a 0-magnitude voltage and a 0-magnitude current, a 1-magnitude voltage and a 1-magnitude current, and a 2-magnitude voltage and a 2-magnitude current.

[0014] Preferably, a calculation formula of the magnitude voltage and the magnitude current is as follows:

[0015]

[0016] wherein S -1 is a phase-to-magnitude transformation matrix, u ma , u mb , u mc are voltage sampling values of an M end, i ma , i mb , i mc are current sampling values of the M end, u na , u nb , u nc are voltage sampling values of an N end, i na , i nb , i nc are current sampling values of the N end, u m0 , u m1 , u m2 are 0-magnitude voltage, 1-magnitude voltage and 2-magnitude voltage values of the M end respectively, i m0 , i m1 , i m2 are 0-magnitude current, 1-magnitude current and 2-magnitude current values of the M end respectively, u n0 , u n1 , u n2 are 0-magnitude voltage, 1-magnitude voltage and 2-magnitude voltage values of the N end respectively, i n0 , in1 , i n2 are the 0-mode, 1-mode and 2-mode current values at the N terminal, respectively, and T represents the transpose of the matrix.

[0017] Preferably, the modulus Baryon difference flow calculation formulae of the M terminal and the N terminal are respectively

[0018]

[0019]

[0020] Wherein: τ mφ is the propagation time of the φ (φ = 0, 1, 2) mode traveling wave between the M terminal protection installation and the fault point, τ nφ is the propagation time of the φ mode traveling wave between the N terminal protection installation and the fault point, i mφ (t) is the φ mode current value of the M terminal at time t, u mφ (t) is the φ mode voltage value of the M terminal at time t, u nφ (t-τ φ ) is the φ mode voltage value of the N terminal at time (t-τ φ ), i nφ (t-τ φ ) is the φ mode current value of the N terminal at time (t-τ φ ), i nφ (t) is the φ mode current value of the N terminal at time t, i nφ (t) is the φ mode voltage value of the N terminal at time t, u mφ (t-τ φ ) is the φ mode voltage value of the M terminal at time (t-τ φ ), i mφ (t-τ φ ) is the φ mode current value of the M terminal at time (t-τ φ );

[0021] Z Cφ is the φ mode wave impedance of the line, τ φ is the φ mode traveling wave propagation time of the whole line, and the calculation formula is:

[0022]

[0023] In the formula: D L is the whole length of the line, v φ is the φ mode wave speed, L φ and C φ are the φ mode inductance and capacitance per unit length of the line, respectively.

[0024] Preferably, the preset current threshold value is in the range of 0.1-0.3 times the rated secondary current.

[0025] Preferably, the distance D of the fault point to the line end M m The calculation formula is:

[0026]

[0027] In the formula, D L is the total length of the line, v φ is the wave velocity of φ (φ = 0, 1, 2).

[0028] Preferably, the distance D of the fault point to the line end N n The calculation formula is:

[0029]

[0030] In the formula, D L is the total length of the line, v φ is the wave velocity of φ (φ = 0, 1, 2).

[0031] In a second aspect, the application provides a fault distance measuring device based on Bergeron differential flow, comprising:

[0032] A collection unit is configured to acquire voltage sampling values and current sampling values at two ends M and N of a line.

[0033] A phase-mode conversion unit is configured to perform phase-mode conversion on the voltage sampling values and current sampling values at the two ends to obtain mode voltage and current at the two ends, respectively.

[0034] A Bergeron differential flow calculation unit is configured to calculate mode quantity Bergeron differential flow at the two ends according to the mode voltage and current at the two ends, wherein the Bergeron differential flow is the difference between an actual value at the current end and a calculated value at the current end calculated by using a Bergeron model at the opposite end.

[0035] A time difference calculation unit is configured to calculate times t m and t n at which the mode quantity Bergeron differential flow at the two ends is equal to a preset current threshold value, respectively.

[0036] A fault distance calculation unit is configured to calculate the difference between the times t m and t n , and calculate the distance of a fault point to a line end in combination with a wave transmission speed.

[0037] In a third aspect, the application provides a computer readable storage medium having a processor program stored thereon, wherein the processor program is configured to execute the method of any one of claims 1 to 8.

[0038] Advantages:

[0039] This application proposes a fault location method based on Berylone differential current. By calculating the time difference between the Berylone differential currents on both sides of the line reaching the same threshold, the time difference between the arrival of the fault traveling wave on both sides of the line is determined. This, combined with wave velocity and line length, allows for the determination of the fault location. This location method has a clear principle, a simple algorithm, does not require identification of the traveling wave front, and is easy to implement in engineering. Attached Figure Description

[0040] Figure 1 This is a flowchart of a fault location method based on Berylone differential current provided in an embodiment of this application.

[0041] Figure 2 This is a schematic diagram of a fault location device based on Berylone differential current provided in an embodiment of this application. Detailed Implementation

[0042] The technical solution and beneficial effects of the present invention will be described in detail below with reference to the accompanying drawings.

[0043] like Figure 1 The figure shown is an embodiment of a fault location method based on Berylone differential current provided in this application, which includes the following steps:

[0044] S110: Acquires voltage and current sampling values ​​at both ends of the line, M and N. Typically, a fiber optic channel is used to synchronously sample the voltage and current at both ends of the line.

[0045] S120: Perform phase-mode transformation on the voltage and current sample values ​​at both ends to obtain the mode voltage and current at both ends, respectively. The phase-mode transformation uses a 3rd order transformation matrix. There are three types of mode voltage and current: 0-mode voltage and current, 1-mode voltage and current, and 2-mode voltage and current; any one of these three can be selected here.

[0046] Two preferred 3rd-order transformation matrices are listed below.

[0047]

[0048] Preferably, the formulas for calculating the mode voltage and mode current are as follows:

[0049]

[0050] Wherein: S -1 Let u be the phase mode transformation matrix. ma u mb u mc The voltage sample value at terminal M, i ma i mb i mc The current sample value at terminal M, u na u nb unc is the voltage sampling value of the N terminal, i na is the current sampling value of the N terminal, u nb is the current sampling value of the N terminal, u nc is the current sampling value of the N terminal, u m0 is the 0-mode voltage value of the M terminal, i m1 is the 1-mode voltage value of the M terminal, i m2 is the 2-mode voltage value of the M terminal, i m0 is the 0-mode current value of the M terminal, i m1 is the 1-mode current value of the M terminal, i m2 is the 2-mode current value of the M terminal, i n0 is the 0-mode voltage value of the N terminal, i n1 is the 1-mode voltage value of the N terminal, i n2 is the 2-mode voltage value of the N terminal, i n0 is the 0-mode current value of the N terminal, i n1 is the 1-mode current value of the N terminal, i n2 is the 2-mode current value of the N terminal, i

[0051] S130: calculating the modulus Berering difference flow of the two terminals according to the modulus voltage and current of the two terminals, the Berering difference flow being the difference between the actual value of the local terminal and the calculated value of the local terminal calculated by the Berering model of the opposite terminal.

[0052] Preferably, the modulus Berering difference flow calculation formulas of the M terminal and the N terminal are respectively

[0053]

[0054]

[0055] wherein: τ mφ is the propagation time of the φ (φ = 0, 1, 2) mode traveling wave between the protection installation of the M terminal and the fault point, τ nφ is the propagation time of the φ mode traveling wave between the protection installation of the N terminal and the fault point, i mφ (t) is the φ mode current value of the M terminal at time t, u mφ (t) is the φ mode voltage value of the M terminal at time t, u nφ (t-τ φ ) is the φ mode voltage value of the N terminal at time (t-τ φ ), i nφ (t-τ φ ) is the φ mode current value of the N terminal at time (t-τ φ ), i nφ (t) is the φ mode current value of the N terminal at time t, i nφ (t) is the φ mode voltage value of the N terminal at time t, u mφ (t-τ φ ) is the φ mode voltage value of the M terminal at time (t-τ φ ), imφ (t-τ φ ) is the φ-mode current value at time t-τ φ ;

[0056] Z Cφ is the φ-mode wave impedance of the line, τ φ is the φ-mode wave propagation time of the full length of the line, and the calculation formula is:

[0057]

[0058] In the formula, D L is the full length of the line, v φ is the φ-mode wave speed, L φ and C φ are the φ-mode inductance and capacitance per unit length of the line, respectively.

[0059] S140: Calculate the time t m and t n at which the modal Bererion difference flows at both ends are equal to the preset current threshold value.

[0060] Preferably, the preset current threshold value is in the range of 0.1-0.3 times the rated secondary current.

[0061] S150: Calculate the difference between the times t m and t n , and calculate the distance from the fault point to the line end in combination with the wave transmission speed. Here, the distance from the fault point to either end M or N of the line can be calculated.

[0062] The distance D m from the fault point to the line end M is calculated by the formula:

[0063]

[0064] In the formula, D L is the full length of the line, and v φ is the φ(φ=0,1,2) mode wave speed.

[0065] The distance D n from the fault point to the line end N is calculated by the formula:

[0066]

[0067] In the formula, D L is the full length of the line, and v φ is the φ(φ=0,1,2) mode wave speed.

[0068] A specific embodiment of a fault location method based on the Bererion difference flow is described below using the 2-mode voltage and current. The method includes the following steps:

[0069] S210: Obtain the voltage and current samples at both ends by using the fiber channel, the sample voltage and sample current at M end are u ma , u mb , u mc and i ma , i mb , i mc , the sample voltage and sample current at N end are u na , u nb , u nc and i na , i nb , i nc .

[0070] S220: Perform phase-mode transformation on the sample voltage and sample current at both ends to obtain the 2-mode voltage and current. The phase-mode transformation matrix used in this embodiment is

[0071]

[0072] The calculation formula for calculating the mode voltage and mode current by using the sample voltage and sample current is

[0073]

[0074] Wherein: S -1 is the phase-mode transformation matrix, u ma , u mb , u mc are the voltage sample values at M end, i ma , i mb , i mc are the current sample values at M end, u na , u nb , u nc are the voltage sample values at N end, i na , i nb , i nc are the current sample values at N end, u m0 , u m1 , u m2 are respectively the 0-mode voltage, 1-mode voltage and 2-mode voltage values at M end, i m0 , i m1 , i m2 are respectively the 0-mode current, 1-mode current and 2-mode current values at M end, u n0 , u n1 , u n2 are respectively the 0-mode voltage, 1-mode voltage and 2-mode voltage values at N end, i n0 , i n1 , i n2 are respectively the 0-mode current, 1-mode current and 2-mode current values at N end, and T represents the transpose of the matrix.

[0075] S230: Calculate the 2-mode Bergeron differential current of the two ends according to the 2-mode voltage and current of the two ends respectively. The 2-mode Bergeron differential current calculation formula of the M end and the N end is respectively

[0076]

[0077]

[0078] Wherein: τ m2 is the propagation time of the 2-mode traveling wave between the M-end protection installation and the fault point, τ n2 is the propagation time of the 2-mode traveling wave between the N-end protection installation and the fault point, i m2 (t) is the 2-mode current value of the M end at time t, u m2 (t) is the 2-mode voltage value of the M end at time t, u n2 (t-τ2) is the 2-mode voltage value of the N end at time (t-τ2), i n2 (t-τ2) is the 2-mode current value of the N end at time (t-τ2), i n2 (t) is the 2-mode current value of the N end at time t, i n2 (t) is the 2-mode voltage value of the N end at time t, u m2 (t-τ2) is the 2-mode voltage value of the M end at time (t-τ2), i m2 (t-τ2) is the 2-mode current value of the M end at time (t-τ2);

[0079] Z c2 is the 2-mode wave impedance of the line, τ2 is the 2-mode traveling wave propagation time of the whole length of the line, and the calculation formula is

[0080]

[0081] In the formula, D L is the whole length of the line, v2 is the 2-mode wave speed, L2 and C2 are the 2-mode inductance and capacitance per unit length of the line respectively.

[0082] S240: Set a preset current threshold value k, calculate the time t m and t n when the Bergeron differential current of the two ends is equal to k. K is an adjustable threshold, which is set to 0.1 times the rated secondary current in this embodiment.

[0083] S250: Calculate the distance D m from the fault point to the M end. The calculation formula of D m is

[0084]

[0085] Wherein, D L is the whole length of the line, v2 is the 2-mode wave speed.

[0086] After the above scheme, the time difference of the traveling wave reaching the two ends of the line is determined by calculating the time when the two-sided Bergeron differential flow is equal to the fixed threshold, and then the distance from the fault point to the M end can be obtained by calculation. The distance measurement method has clear principle, simple algorithm, does not need to identify the traveling wave head, and is easy to implement in engineering.

[0087] As shown in Figure 2 The application provides a fault distance measurement device 1000 based on Bergeron differential flow, which comprises a collection unit 1001, a phase-mode conversion unit 1002, a Bergeron differential flow calculation unit 1003, a time difference calculation unit 1004 and a fault distance calculation unit 1005 connected in sequence.

[0088] The collection unit 1001 is used to obtain voltage sampling values and current sampling values at the two ends M and N of the line;

[0089] The phase-mode conversion unit 1002 is used to perform phase-mode conversion on the voltage sampling values and current sampling values at the two ends to obtain the mode voltage and current at the two ends, respectively.

[0090] The Bergeron differential flow calculation unit 1003 is used to calculate the modulus Bergeron differential flow at the two ends according to the mode voltage and current at the two ends, wherein the Bergeron differential flow is the difference between the actual value at the local end and the calculated value at the local end calculated by using the Bergeron model at the opposite end;

[0091] The time difference calculation unit 1004 is used to calculate the time t m and t n when the modulus Bergeron differential flow at the two ends is equal to a preset current threshold value, respectively.

[0092] The fault distance calculation unit 1005 is used to calculate the difference between the times t m and t n , and calculate the distance from the fault point to the end of the line in combination with the traveling wave transmission speed.

[0093] In the phase-mode conversion unit 1002, the phase-mode conversion adopts a 3-order conversion matrix.

[0094] Preferably, the mode voltage and current are any one of 0-mode voltage and current, 1-mode voltage and current, and 2-mode voltage and current.

[0095] Preferably, the calculation formula for calculating the mode voltage and current in the phase-mode conversion unit 1002 is as follows:

[0096]

[0097] Wherein: S -1 is a phase-mode conversion matrix, u ma , u mb , umc is the voltage sampling value of the M terminal, i ma is the voltage sampling value of the M terminal, i mb is the voltage sampling value of the M terminal, i mc is the current sampling value of the M terminal, u na is the current sampling value of the M terminal, u nb is the current sampling value of the M terminal, u nc is the voltage sampling value of the N terminal, i na is the voltage sampling value of the N terminal, i nb is the voltage sampling value of the N terminal, i nc is the current sampling value of the N terminal, u m0 is the current sampling value of the N terminal, u m1 is the current sampling value of the N terminal, u m2 is the 0-mode voltage value of the M terminal, i m0 is the 0-mode voltage value of the M terminal, i m1 is the 0-mode voltage value of the M terminal, i m2 is the 0-mode current value of the M terminal, u n0 is the 0-mode current value of the M terminal, u n1 is the 0-mode current value of the M terminal, u n2 is the 0-mode voltage value of the N terminal, i n0 is the 0-mode voltage value of the N terminal, i n1 is the 0-mode voltage value of the N terminal, i n2 is the 0-mode current value of the N terminal, T represents the transpose of a matrix.

[0098] Preferably, the modulus Berge differential current calculation formula of the M terminal and the N terminal in the Berge differential current calculation unit 1003 is respectively

[0099]

[0100]

[0101] wherein: τ mφ is the propagation time of the φ (φ = 0, 1, 2) mode traveling wave between the protection installation of the M terminal and the fault point, τ nφ is the propagation time of the φ mode traveling wave between the protection installation of the N terminal and the fault point, i mφ (t) is the φ mode current value of the M terminal at time t, u mφ (t) is the φ mode voltage value of the M terminal at time t, u nφ (t-τ φ ) is the φ mode voltage value of the N terminal at time (t-τ φ ), i nφ (t-τ φ ) is the φ mode current value of the N terminal at time (t-τ φ ), i nφ (t) is the φ mode current value of the N terminal at time t, i nφ (t) is the φ mode voltage value of the N terminal at time t, u mφ (t-τ φ ) is the φ mode voltage value of the M terminal at time (t-τφ ) the φ-mode voltage value at time t-τ mφ (t-τ φ ) is the φ-mode current value at time t-τ φ

[0102] Z Cφ is the φ-mode wave impedance of the line, τ φ is the φ-mode wave propagation time of the whole length of the line, and the calculation formula is:

[0103]

[0104] In the formula, D L is the whole length of the line, v φ is the φ-mode wave speed, L φ and C φ are the φ-mode inductance and capacitance per unit length of the line, respectively.

[0105] Preferably, the preset current threshold value in the time difference calculation unit 1004 is in the range of 0.1-0.3 times the rated secondary current.

[0106] Preferably, the fault distance calculation unit 1005 calculates the distance D m from the fault point to the line end M, and the calculation formula is:

[0107]

[0108] In the formula, D L is the whole length of the line, v φ is the φ-mode wave speed (φ=0, 1, 2).

[0109] Preferably, the fault distance calculation unit 1005 calculates the distance D n from the fault point to the line end N, and the calculation formula is:

[0110]

[0111] In the formula, D L is the whole length of the line, v φ is the φ-mode wave speed (φ=0, 1, 2).

[0112] The embodiment of the application further provides a computer readable storage medium, which stores a processor program, wherein the processor program is used for executing the fault distance measurement method based on the Bererion differential flow.

[0113] ​Those skilled in the art will appreciate that embodiments of the application can be provided as methods or as computer program products. The application can take the form of a computer program product embodied in one or more computer-usable storage media (including, but not limited to, disk memory and optical memory) having computer-usable program code embodied therein.

[0114] The present application is described in reference to the flowchart illustrations and / or block diagrams according to the embodiments of the application. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart illustrations and / or block diagrams block or blocks. Figure 1 one or more functions specified in the flowchart illustrations and / or block diagrams block or blocks. Figure 1 one or more functions specified in the flowchart illustrations and / or block diagrams block or blocks.

[0115] These computer program instructions can also be stored in a computer- readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instructions which implement the function specified in the flowchart illustrations and / or block diagrams block or blocks. Figure 1 one or more functions specified in the flowchart illustrations and / or block diagrams block or blocks. Figure 1 one or more functions specified in the flowchart illustrations and / or block diagrams block or blocks.

[0116] These computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the function specified in the flowchart illustrations and / or block diagrams block or blocks. Figure 1 one or more functions specified in the flowchart illustrations and / or block diagrams block or blocks. Figure 1 one or more functions specified in the flowchart illustrations and / or block diagrams block or blocks.

[0117] The above embodiments are only illustrative of the technical idea of the present application, and cannot limit the protection scope of the present application. Any modification made according to the technical idea of the present application on the basis of the technical solution falls within the protection scope of the present application.

Claims

1. A fault location method based on the differential flow of the Berrilon, characterized by, The method comprises the following steps: obtaining synchronous voltage sample values and current sample values of M end and N end at both ends of a line; performing phase-to-magnitude transformation on the voltage sample values and current sample values at both ends to obtain magnitude voltage and current at both ends respectively; the phase-to-magnitude transformation adopts a 3-order conversion matrix; calculating magnitude Bergeron differential flow at both ends according to the magnitude voltage and current at both ends; the Bergeron differential flow is the difference between an actual value at the end and a calculated value at the end calculated by adopting a Bergeron model at the opposite end; The time t when the modulus Barylon difference flow at both ends is equal to the preset current threshold value m and t n ; the preset current threshold value ranges from 0.1 to 0.3 times the rated secondary current. The time t is calculated m The difference between t n and t is calculated, in combination with the wave transmission speed, to calculate the distance from the fault point to the line end.

2. The method for fault location based on the Berlincourt differential current as claimed in claim 1, wherein: the magnitude voltage and current are any one of 0-magnitude voltage and current, 1-magnitude voltage and current, and 2-magnitude voltage and current.

3. The fault location method based on Berylone differential current as described in claim 1, characterized in that: The calculation formula of the magnitude voltage and current is: wherein: S -1 is a phase modulation transformation matrix, u ma , u mb , u mc is a voltage sampling value of the M terminal, i ma , i mb , i mc is a current sampling value of the M terminal, u na , u nb , u nc is a voltage sampling value of the N terminal, i na , i nb , i nc is a current sampling value of the N terminal, u m0 , u m1 , u m2 are respectively 0-mode voltage, 1-mode voltage and 2-mode voltage values of the M terminal, i m0 , i m1 , i m2 are respectively 0-mode current, 1-mode current and 2-mode current values of the M terminal, u n0 , u n1 , u n2 are respectively 0-mode voltage, 1-mode voltage and 2-mode voltage values of the N terminal, i n0 , i n1 , i n2 are respectively 0-mode current, 1-mode current and 2-mode current values of the N terminal, and T represents a transpose of a matrix.

4. The fault location method based on Bergeron differential flow according to claim 1, characterized in that, the calculation formula of the magnitude Bergeron differential flow at M end and N end is wherein: τ mφ is the propagation time of the φ (φ = 0, 1, 2) mode between the point of fault and the installation of the protection at the M end, τ nφ is the propagation time of the φ mode between the point of fault and the installation of the protection at the N end, i mφ (t) is the φ mode current value at the M end at time t, u mφ (t) is the φ mode voltage value at the M end at time t, u nφ (t - τ φ ) is the φ mode voltage value at the N end at time (t - τ φ ), i nφ (t - τ φ ) is the φ mode current value at the N end at time (t - τ φ ), i nφ (t) is the φ mode current value at the N end at time t, i nφ (t) is the φ mode voltage value at the N end at time t, u mφ (t - τ φ ) is the φ mode voltage value at the M end at time (t - τ φ ), i mφ (t - τ φ ) is the φ mode current value at the M end at time (t - τ φ ). Z Cφ φ mode wave impedance of the line, τ φ φ mode time of flight of the line, calculated as: where D L is the full length of the line, v φ is the modal wave speed, L φ and C φ are the per-unit length line inductance and capacitance, respectively.

5. The fault location method based on Berylone differential current as described in claim 1, characterized in that: Distance D of the fault point to the line end M m The calculation formula is: In the formula, D L is the full length of the line, v φ is the phase (φ = 0, 1, 2) wave velocity.

6. The method of claim 1, wherein: Distance D of the fault point to the line end N n The calculation formula is: where D L is the full length of the line, v φ is the phase (φ = 0, 1, 2) mode wave velocity.

7. A fault location device based on the Barelon differential flow, characterized by, comprise: a collection unit configured to obtain synchronous voltage sample values and current sample values of M end and N end at both ends of a line; a phase-to-magnitude transformation unit configured to perform phase-to-magnitude transformation on the voltage sample values and current sample values at both ends to obtain magnitude voltage and current at both ends respectively; the phase-to-magnitude transformation adopts a 3-order conversion matrix; a Bergeron differential flow calculation unit configured to calculate magnitude Bergeron differential flow at both ends according to the magnitude voltage and current at both ends; the Bergeron differential flow is the difference between an actual value at the end and a calculated value at the end calculated by adopting a Bergeron model at the opposite end; A time difference calculation unit is configured to calculate a time t at which the modulus Berreman difference flows at both ends are respectively equal to a preset current threshold value m and t n ; the preset current threshold value is in a range of 0.1-0.3 times of a rated secondary current. a fault distance calculation unit for calculating the difference between the times t m and t n in combination with the wave propagation transmission speed to calculate the distance of the fault point to the line end.

8. A computer readable storage medium having stored thereon a processor program, wherein, the processor program is configured to execute the method according to any one of claims 1 to 6.

Images