CN115508671B - Fault positioning method and system based on line voltage lowest point calculation - Google Patents
Fault positioning method and system based on line voltage lowest point calculation Download PDFInfo
- Publication number
- CN115508671B CN115508671B CN202211401819.XA CN202211401819A CN115508671B CN 115508671 B CN115508671 B CN 115508671B CN 202211401819 A CN202211401819 A CN 202211401819A CN 115508671 B CN115508671 B CN 115508671B
- Authority
- CN
- China
- Prior art keywords
- voltage
- calculation
- curve
- line
- point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
- G01R31/085—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution lines, e.g. overhead
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/088—Aspects of digital computing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
- Y04S10/52—Outage or fault management, e.g. fault detection or location
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Locating Faults (AREA)
Abstract
The invention discloses a fault positioning method and a fault positioning system based on voltage lowest point calculation along a line, wherein the method comprises the steps of performing voltage calculation along the line based on a Bergeron model, and performing complex integral calculation in a unit data window to obtain a voltage distribution curve; and calculating a voltage distribution curve, and performing fault location judgment on the line according to the calculation result of the voltage distribution curve. The invention is less influenced by fault transient process and line distributed capacitance, has higher fault positioning speed and precision, is hardly influenced by different fault moments, and greatly improves the fault location accuracy of the long-distance transmission line.
Description
Technical Field
The invention relates to the technical field of power system relay protection, in particular to a fault positioning method and system based on line voltage lowest point calculation.
Background
With the development of economy, the demand for long-distance and large-range power transmission is increased, and meanwhile, high-voltage power transmission lines are widely distributed, and the terrain of a crossing area is complex and changeable. The fault location is rapidly and accurately carried out after the transmission line has a fault, so that the fault is rapidly repaired, the stability of a power system is improved, and the fault location method is also very important for normal operation of national economy.
The single-ended impedance method for fault location only needs to measure information of one end of the power transmission line, and has the advantages of low requirement on hardware, easiness in implementation, stable algorithm and the like, so that the method is widely applied. The traditional fault location uses a line model based on centralized parameters, and when the fault location is applied to a longer power transmission line, the influence of line distributed capacitance becomes non-negligible, and a larger error is caused to the remote fault location precision. Therefore, it is very important to research a fault location model and algorithm suitable for long lines and improve the precision of fault location.
Disclosure of Invention
The invention aims to provide a fault location method and a fault location system based on the calculation of the lowest point of the voltage along the line, which are less influenced by the fault transient process and the distributed capacitance of the line, have higher fault location speed and precision, are hardly influenced by different fault moments, and greatly improve the fault location accuracy of a long-distance power transmission line
In order to achieve the above object, the present invention provides a fault location method based on the calculation of the lowest point of the line voltage, which comprises:
performing voltage calculation along the line based on the Bergeron model, and performing complex integral calculation in a unit data window to obtain a voltage distribution curve;
and calculating a voltage distribution curve, and performing fault location judgment on the line according to the calculation result of the voltage distribution curve.
Further, the obtaining of the voltage distribution curve includes,
determining the starting time of the complex integral calculation according to the full length of the line;
and determining a line fault phase through a phase selection element based on the initial moment, and calculating a voltage distribution curve of the corresponding fault phase through a Bergeron model according to the voltage and the current obtained by sampling.
Further, the starting time should satisfy the following condition:
wherein the content of the first and second substances,t 1 which is indicative of the time of the start,lin order to obtain the full length of the line,representing the mode wave velocity.
Further, the voltage distribution curve is:
wherein the content of the first and second substances,U inte (x) Which represents the curve of the voltage distribution,xindicating the distance, Δ, of any point on the line from the protective mountingTWhich represents the time between the sampling of the samples,urepresenting the instantaneous voltage at a point on the line,t 1 indicating a starting time within a unit data window,t 2 indicating the stop time within a unit data window,Nis shown in (t 2 -t 1 ) The number of voltage samples in the period;iindicating the number of calculations of the voltage profile.
Further, calculating a voltage distribution curve, and performing fault location judgment on the line according to the calculation result of the voltage distribution curve, including,
calculating the required line step length according to the Bergeron model, and determining the sampling rate and the interpolation mode of the line protection device;
calculating instantaneous values of step voltages of all points of a line from the initial time to the stopping time in a fault phase unit data window by using a Bergeron model;
and (4) calculating the integral value of the voltage of each point through the voltage distribution curve, comparing and judging to obtain the minimum point of the voltage integral value, and confirming the fault point.
Further, the calculating the required line step according to the berelon model includes,
the line step is obtained by:
wherein, the first and the second end of the pipe are connected with each other,represents a line step, <' > based on>Representing the velocity of the mode wave, ΔTRepresenting the sampling interval time.
Furthermore, the integral value of the voltage of each point is obtained through a voltage distribution curve for comparison and judgment, the minimum point of the voltage integral value is obtained, the fault point is confirmed, and the method comprises the following steps of,
selecting calculation interval time, and calculating a voltage distribution curve every other calculation interval time to obtain a complete complex voltage integral curve;
defining a comparison windowI.e. each calculation will have the latest successive->Carrying out cyclic comparison judgment on the integral value of the secondary complex voltage integral curve;
after calculating the initial time, continuously carrying out the operation of the complex voltage integral curveStarting the first judgment after the secondary calculation;
defining a time window for overall computationt set When the fault positioning judgment condition is not met and the current time is less than the time windowt set And then, along with the increase of the current moment, the complete cyclic judgment calculation of the complex voltage integral curve is continuously carried out.
Further, selecting calculation interval time, calculating a voltage distribution curve every other calculation interval time to obtain a complete complex voltage integral curve, including,
before each time of calculating the voltage distribution curve, obtaining the current time according to the calculation interval time and the initial time;
the current time is obtained by:
wherein the content of the first and second substances,tas the current time of day, the time of day,t 1 which indicates the starting moment in time of the day,nis shown asnThe sub-cycle calculates a complex voltage integral curve,Δ tthe table calculates the interval time.
Furthermore, after the initial time is calculated, the voltage integral curve is continuously repeatedStarting the first judgment after the second calculation, and judging the firstnUpon sub-calculation of a complex voltage integration curve>Whether or not it is greater than or equal ton:
If it isThen continue to return to executenCalculating a voltage integral curve repeatedly for +1 times and entering judgment cycle calculation again;
if it isIf so, judge thatnAnd whether a plurality of minimum value points exist in the complex voltage integral curve obtained by secondary calculation or not.
Further, ifIf so, it is judged asnWhether a plurality of minimum value points exist in the complex voltage integral curve obtained by the secondary calculation or not comprises,
if there are multiple minimum values in the said integrated curve, define the firstnThe minimum value point with the minimum integral value in the complex voltage integral curve obtained by the secondary calculation isThe next smallest extreme point of the integration value is->Judgment is made>And/or>Whether or not the integrated value of (b) satisfies a determination condition; wherein the content of the first and second substances,
if the current time of the complex voltage integral curve is still in the total calculation time windowt set If so, after a calculation interval time, continuing to execute the step onenCalculating a voltage integration curve for +1 times, and entering a judgment cycle again;
if the current time for calculating the complex voltage integral curve exceeds the total calculation time windowt set If yes, the loop is exited and the fault location is not possible;
when in useAnd &>If the integrated value of (a) meets the decision condition, then the comparison window is compared for consecutive @>And judging the minimum value point of the secondary composite voltage integral curve.
Further, the judgment is made by the following formulaAnd/or>Whether or not the integrated value of (b) satisfies the determination condition:
wherein, the first and the second end of the pipe are connected with each other,represents the minimum point in the duplicated voltage integration curve>Based on the value integrated in>Represents the next smallest extreme point in the duplicated voltage integration curve>Is greater than or equal to>And the setting value which excludes other minimum value point interference is judged.
Further, if<nIf so, it is judged asnWhether a plurality of minimum value points exist in the complex voltage integral curve obtained by the secondary calculation or not comprises,
if the plurality of minimum value points do not exist in the complex voltage integral curve, judging whether the complex voltage integral curve has one minimum value point or not; wherein the content of the first and second substances,
when the complex voltage integral curve has one and only one minimum value point, the continuous comparison window is comparedJudging minimum value points of the secondary repetition voltage integral curve;
when there is no minimum point on the complex voltage integral curve:
if the current time of the complex voltage integral curve is still in the total calculation time windowt set If so, after a calculation interval time, continuing to execute the step onenCalculating a voltage integration curve for +1 times, and entering a judgment cycle again;
if the current time for calculating the complex voltage integral curve exceeds the total calculation time windowt set The loop exits and the fault location is deemed to be impossible.
Further, the continuity in the comparison window is compared by the following formulaAnd (3) judging minimum value points of the secondary complex voltage integral curve:
wherein the content of the first and second substances,x n representing continuities within a comparison windowThe minimum point of the second-time-multiplexed voltage integral curve,x i is shown asiAs a result of the secondary fault location,ithe number of calculations for fault location.
Further, when the comparison window is continuousWhen the minimum value point of the secondary multiplexing voltage integral curve meets the condition of a judgment formula, determining that the calculated value of the fault distance is greater or less than the preset value>;
When continuous in the comparison windowWhen the minimum value point of the secondary complex voltage integral curve does not meet the condition of the judgment formula; wherein, the first and the second end of the pipe are connected with each other,
if the current time of the complex voltage integral curve is still in the total calculation time windowt set If so, after a calculation interval time, continuing to execute the step onenCalculating a voltage integration curve for +1 times, and entering a judgment cycle again;
if the current time for calculating the complex voltage integral curve exceeds the total calculation time windowt set And the calculation range of the output fault distance is the minimum value point of the complex voltage integral curve in the comparison time windowTo a maximum minimum point->。
The invention also provides a fault positioning system based on the calculation of the lowest point of the line voltage, which comprises a calculating unit and a judging unit,
the calculating unit is used for performing voltage calculation along the line based on the Bergeron model and performing complex integration in a unit data window to obtain a voltage distribution curve;
and the judging unit is used for calculating a voltage distribution curve and carrying out fault positioning judgment on the line according to the result of the voltage distribution curve.
Further, the calculating unit is configured to determine a starting time of the complex integral calculation according to the full length of the line;
the calculating unit is used for determining a fault phase through the phase selection element based on the initial moment, and calculating a voltage distribution curve of the corresponding fault phase through the Bergeron model according to the sampled voltage and current.
Further, the determination unit is configured to calculate a required line step length according to the berelon model, and determine a sampling rate and an interpolation mode of the line protection device;
calculating instantaneous values of step voltages of all points of a line from the initial time to the stopping time in a fault phase unit data window by using a Bergeron model;
and (4) calculating the integral value of the voltage of each point through the voltage distribution curve, comparing and judging to obtain the minimum point of the voltage integral value, and confirming the fault point.
The invention has the technical effects and advantages that: according to the invention, the voltage along the line is calculated based on the Bergeron model, and the complex integration is carried out in a unit data window, so that a complex voltage integration curve is obtained. Considering the different situations that the integral curve has only one minimum value point, a plurality of minimum value points and no minimum value point, the comprehensive analysis provides a fault positioning scheme based on the calculation of the lowest point of the voltage along the line. The invention is less influenced by fault transient process and line distributed capacitance, has higher fault positioning speed and precision, is hardly influenced by different fault moments, and greatly improves the fault location accuracy of the long-distance transmission line.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings required to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the description below are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a voltage distribution diagram illustrating an exemplary metallic single-phase ground fault in an embodiment of the present invention;
FIG. 2 is a flow chart of a method for fault location based on along-line voltage nadir calculation in an embodiment of the present invention;
FIG. 3 is a flowchart of a method for computing a fault in a fault location based on along-line voltage nadir computation in accordance with an embodiment of the present invention;
FIG. 4 is a partial flow chart of a calculation determination fault in a fault location method based on along-line voltage nadir calculation in accordance with an embodiment of the present invention;
FIG. 5 is a partial flowchart of a method for computing a fault in a fault location based on along-line voltage nadir computation in accordance with an embodiment of the present invention;
fig. 6 is a schematic structural diagram of a fault location system based on line voltage minimum point calculation according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
When breaking down on the transmission line, the bersulfuron model of whole transmission line is destroyed, can't calculate the whole line voltage distribution condition through single-ended bersulfuron model, nevertheless still satisfies the bersulfuron model between protection installation department and the fault point this moment, and accessible bersulfuron model calculates and seeks the fault point: for three-phase symmetrical short circuit and single-phase earth fault, the voltage value at the fault point of the short circuit earth phase is minimum; for a two-phase short circuit grounding fault, the phase voltage or the line voltage at the fault point of the short circuit grounding phase is the minimum value; for a two-phase short-circuit fault, the line voltage at the short-circuit phase fault point has a minimum value. Taking a metallic single-phase earth fault as an example, the voltage distribution is shown in fig. 1, and the equivalent power voltages of the power grids at two ends of the line are respectively shown asE m AndE n equivalent supply voltage acrossE m AndE n are respectively connected with system impedanceZ m AndZ n and the voltages at the two ends of the line are respectivelyU m AndU n if a fault at a certain point of the line is a fault pointFThen a point of failure in the lineFThe voltage of (2) is 0.
Therefore, after the protection is started after the line fault, after the fault phase is confirmed through the phase selection element, the corresponding voltage distribution is calculated through the Bergeron model, and then the fault point can be confirmed. However, due to various factors such as sampling noise, there may be a large error in calculating the line voltage distribution once.
In order to solve the defects of the prior art, the invention selects the voltage data of a certain data window to carry out the complex trapezoidal integration, and then the complex integral value of the line voltage is obtainedUinte(x)And comparing, searching for a fault point and enhancing the anti-interference capability.
The invention discloses a fault positioning method based on line voltage lowest point calculation, which comprises the following steps of:
calculating voltage along the line based on a Bergeron model, and performing complex integration in a unit data window to obtain a voltage distribution curve;
and calculating a voltage distribution curve, and performing fault location judgment on the line according to the result of the voltage distribution curve.
In particular to a method for preparing a high-performance nano-silver alloy,
and determining the step length required by calculation of the Bergeron equation according to the requirement of protection precision, thereby determining the sampling rate and the interpolation mode of the protection device. The line step is obtained by:
wherein the content of the first and second substances,X step which represents the step size of the line,representing the velocity of the mode wave, ΔTRepresenting the sampling interval time.
After the fault occurs, calculating the starting time to the stopping time in the fault phase unit data window by using a Bergeron model (t 2 -t 1 ) Step length of each point of the internal line (x step ,2x step ,3x step …) voltage.
It should be noted that, according to the characteristics of the beraprone model, the on-line routeStarting time evaluation line>At a voltage of (a) is actually requiredt 1 -x/v m ) And (t) 1 +x/v m ) Sampled data of time instants. For the sake of explanation, assume the wave velocity v m 3 x 105km/s, if the voltage at the position of 300km of the line 1ms after the fault is calculated, the voltage isx/v m = 1ms, the calculation also requires voltage and current sampling data at 0 time after the failure and 2ms after the failure. />
And after reading the voltage and current information at the protection installation position, carrying out complex voltage integral calculation positioning.
Firstly, determining the starting time of the complex integral calculation according to the full length of the line, wherein the starting time meets the following condition:
wherein, the first and the second end of the pipe are connected with each other,represents the starting time, l is the full length of the line,/>Representing the mode wave velocity.
The fault phase of the line is determined through the phase selection element, and after protection is started, the voltage distribution curve of the corresponding fault phase is calculated based on the voltage and the current obtained through samplingU inte (x) The calculation formula is as follows:
wherein the content of the first and second substances,U inte (x) A voltage profile is shown in the form of a voltage profile,xindicating the distance, Δ, of any point on the line from the protective mountingTWhich represents the time between the sampling of the samples,urepresenting the instantaneous voltage at a point on the line,t 1 indicating a starting time within a unit data window,t 2 indicating the stop time within a unit data window,Nis shown in (A)t 2 -t 1 ) The number of voltage samples in the period,iexpressed as the number of calculations of the voltage profile and finallyi=N-2。
Voltage distribution curve by equation (2)U inte (x) And (4) calculating and comparing the integral value of the voltage of each point, wherein the minimum point of the integral value of the voltage is the fault point. In particular to a method for preparing a high-performance nano-silver alloy,
next, n =0,n was set as the number of times the voltage distribution curve was calculated before the calculation. Selecting proper calculation interval timeAnd every other->Calculating a primary voltage distribution curve, and obtaining a complete complex voltage integral curve after n times of calculation;
before each voltage distribution curve is calculated, the current time is obtained according to the starting time and the calculation interval time, so that the complex voltage integral curve is judged in the subsequent cycle judgment, and the current time is obtained according to the following formula:
wherein the content of the first and second substances,tas the current time of day, the time of day,t 1 which indicates the starting moment in time of the day,nis shown asnThe sub-cycle calculates a complex voltage integral curve,Δ tthe table calculates the interval time.
Defining a comparison windowI.e. each calculation will have the latest successive->And carrying out cyclic comparison judgment on the integral value of the calculation result of the secondary complex voltage integral curve so as to enhance the reliability of fault positioning.
Therefore, at the calculation start timet 1 Then, the complex voltage integral curve is required to be continuously carried outThe first decision is started after the second calculation.
Defining a time window for overall computationt set When the fault positioning judgment condition is not met and the current momenttLess than the time windowt set At the current momenttThe calculation of the complete complex voltage integral curve continues.
Finally, the determination is performed according to the calculation result of the complete complex voltage integral curve, and the calculation and determination process of the complex voltage integral curve is shown in fig. 3, 4 and 5:
after calculating the starting time, obtaining the current time and the integrated curve of the complex voltageIs taken continuously>Starting the first judgment after the second calculation, judging the firstnUpon sub-calculation of a complex voltage integration curve>Whether or not it is greater than or equal ton(ii) a Wherein the content of the first and second substances,
if it isIf n is greater than or equal to n, the execution is continued to return ton+1 time multiplexing voltage integration curve->And calculating and entering the judgment loop again for calculation.
If it is<nIf so, it is judged asnThe recalculated voltage integration curve->Whether a plurality of minimum value points exist; the method comprises the following steps: />
If it is firstnThe integrated curve of the complex voltage obtained by the secondary calculationIf there are multiple minimum value points, defining the minimum value point of the integral value in the complex voltage integral curve as ^ greater than or equal to>The extreme point with the next smaller integrated value is->And the integral value corresponding to the two isAnd &>(ii) a Judgment->And/or>Whether or not the integrated value of (2) satisfies the determination condition.
Is judged by the following formulaBreak-offAnd/or>Whether or not the integrated value of (b) satisfies a determination condition:
wherein the content of the first and second substances,represents the minimum point in the duplicated voltage integration curve>Is greater than or equal to>Represents the next smallest extreme point in a complex voltage integration curve>Is greater than or equal to>Setting values for eliminating interference of other minimum value points are judged;
the smaller the value, i.e.. Sup>And/or>The greater the phase difference, at which a decision can be made>The greater the confidence probability for the fault distance. />The greater the value, i.e.>And/or>The smaller the phase difference is, the moment->The smaller the confidence probability for the distance to failure.
if the current time of the complex voltage integral curve is still in the total calculation time windowt set In, i.e.t<t set Then, after a calculation interval time, the execution is continuednCalculating a composite voltage integral curve +1 times, and entering a judgment cycle again;
if the current time for calculating the complex voltage integral curve exceeds the total calculation time windowt set I.e. byt≥t set The loop is exited and it is believed that fault location by the present method is not possible at this time due to transition resistance, etc.
When the temperature is higher than the set temperatureAnd &>When the integrated value of (4) meets a decision condition, then successive in the comparison window are compared>Judging minimum value points of the secondary complex voltage integral curve; wherein the content of the first and second substances,
comparing the window internal continuity byAnd (3) judging minimum value points of the secondary composite voltage integral curve:
wherein the content of the first and second substances,indicates a continuation in the comparison window>The minimum point of the second-order voltage integral curve,x i denotes the firstiAs a result of the secondary fault location, the fault location,ithe number of calculations for fault location.
When continuous in the comparison windowWhen the minimum value point of the secondary multiplexing voltage integral curve meets the condition of the judgment formula (5), determining that the fault distance calculation value is greater than or equal to>。
When continuous in the comparison windowWhen the minimum value point of the second-order composite voltage integral curve does not meet the condition of the judgment formula (5); wherein the content of the first and second substances,
if the current time of the complex voltage integral curve is still in the total calculation time windowt set In, i.e.t<t set If the calculation interval is one, the method continues to executenCalculating a voltage integration curve for +1 times, and entering a judgment cycle again;
if the current time for calculating the complex voltage integral curve exceeds the total calculation time windowt set I.e. byt≥t set If the calculated range of the output fault distance is the minimum value point of the complex voltage integral curve in the comparison time windowTo a maximum minimum point->。
If the complex voltage integral curve obtained by the nth calculation does not have a plurality of minimum value points, judging whether the complex voltage integral curve has one minimum value point or not; wherein, the first and the second end of the pipe are connected with each other,
when the complex voltage integral curve has one and only one minimum value point, the continuous comparison window is comparedJudging minimum value points of the secondary multiplexing voltage integral curve; wherein the content of the first and second substances,
the decision is made by equation (5):
when continuous in the comparison windowWhen the minimum value point of the secondary multiplexing voltage integral curve meets the condition of the judgment formula (5), determining that the fault distance calculation value is greater than or equal to>。
When continuous in the comparison windowMinimum point of secondary voltage integral curveWhen the condition of the decision formula (5) is not satisfied; wherein the content of the first and second substances,
if the current time of the complex voltage integral curve is still in the total calculation time windowt set In, i.e.t<t set Then, after a calculation interval time, the execution is continuednAnd calculating a voltage integral curve for +1 times, and entering a judgment loop again.
If the current time for calculating the complex voltage integral curve exceeds the total calculation time windowt set I.e. byt≥t set If the calculated range of the output fault distance is the minimum value point of the complex voltage integral curve in the comparison time windowTo a maximum minimum point->。
When there is no minimum point on the complex voltage integral curve:
if the current time of the complex voltage integral curve is still in the total calculation time windowt set In, i.e.t<t set Then, after a calculation interval time, the execution is continuednAnd calculating a voltage integral curve for +1 times, and entering a judgment loop again.
If the current time for calculating the complex voltage integral curve exceeds the total calculation time windowt set I.e. byt≥t set The loop is exited and it is believed that the fault location cannot be performed by the present method due to, for example, transition resistance.
The invention also discloses a fault location system based on the calculation of the lowest point of the voltage along the line, as shown in fig. 4, the system comprises a calculation unit and a determination unit,
and the calculating unit is used for calculating the voltage along the line based on the Bergeron model and performing complex integration in a unit data window to obtain a voltage distribution curve.
Specifically, the calculating unit is used for determining the starting moment of the complex integral calculation according to the full length of the line; the calculating unit is used for determining a fault phase through the phase selection element based on the initial moment, and calculating a voltage distribution curve of the corresponding fault phase through the Bergeron model according to the voltage and the current obtained through sampling.
And the judging unit is used for calculating a voltage distribution curve and carrying out fault positioning judgment on the line according to the result of the voltage distribution curve.
Specifically, the judging unit is used for calculating the required line step length according to the Bergeron model and determining the sampling rate and the interpolation mode of the line protection device;
the judging unit is used for calculating the instantaneous value of the step voltage of each point of the line from the starting time to the stopping time in the fault phase unit data window by using the Bergeron model;
and the judging unit is used for obtaining the integral value of the voltage at each point through the voltage distribution curve, comparing and judging to obtain the minimum point of the voltage integral value and confirm the fault point.
With regard to the system in the above embodiment, the specific manner in which each unit module performs the operation has been described in detail in the embodiment related to the method, and will not be elaborated here.
Finally, it should be noted that: although the present invention has been described in detail with reference to the foregoing embodiments, it will be apparent to those skilled in the art that modifications may be made to the embodiments or portions thereof without departing from the spirit and scope of the invention.
Claims (15)
1. A fault location method based on the calculation of the lowest point of voltage along the line is characterized by comprising the following steps,
performing voltage calculation along the line based on the Bergeron model, and performing complex integral calculation in a unit data window to obtain a voltage distribution curve;
calculating a voltage distribution curve, and performing fault location judgment on the line according to a calculation result of the voltage distribution curve;
wherein the obtaining of the voltage profile comprises the following operations:
determining the starting time of the complex integral calculation according to the full length of the line;
determining a line fault phase through a phase selection element based on the initial moment, and calculating a voltage distribution curve of the corresponding fault phase through a Bergeron model according to the voltage and the current obtained by sampling;
the voltage distribution curve is as follows:
wherein, the first and the second end of the pipe are connected with each other,Uinte(x)which represents the curve of the voltage distribution,xindicating the distance of any point on the line from the protective installation,ΔTwhich represents the time between the sampling of the samples,urepresenting the instantaneous voltage at a point on the line,t 1 indicating a starting time within a unit data window,t 2 indicating the stop time within a unit data window,Nis shown in(t 2 -t 1 )The number of voltage samples in the period;iindicating the number of calculations of the voltage profile.
2. The method of claim 1 wherein the fault location is based on line voltage nadir calculation,
the starting time should satisfy the following condition:
3. The fault location method based on the line voltage minimum point calculation according to claim 1 or 2, wherein a voltage distribution curve is calculated, and a fault location determination is made on a line according to the calculation result of the voltage distribution curve, comprising,
calculating the required line step length according to the Bergeron model, and determining the sampling rate and the interpolation mode of the line protection device;
calculating the instantaneous value of the step voltage of each point of the circuit from the initial moment to the stop moment in a unit data window of the fault phase by using a Bergeron model;
and (4) calculating the integral value of the voltage of each point through the voltage distribution curve, comparing and judging to obtain the minimum point of the voltage integral value, and confirming the fault point.
4. The fault location method based on the line voltage nadir calculation according to claim 3, wherein the calculating the required line step size according to the Bergeron model comprises,
the line step is obtained by:
5. The method for locating a fault based on the calculation of the lowest point of the line voltage according to claim 3, wherein the integral value of the voltages at each point is obtained through the voltage distribution curve for comparison and determination, the minimum point of the integral value of the voltages is obtained, and the fault point is identified, comprising,
selecting calculation interval time, and calculating a voltage distribution curve every other calculation interval time to obtain a complete complex voltage integral curve;
defining a comparison windowI.e. each calculation will be the latest successiveCarrying out cyclic comparison judgment on the integral value of the secondary composite voltage integral curve;
after calculating the initial time, continuously carrying out the operation of the complex voltage integral curveStarting the first judgment after the second calculation;
defining a time window for overall computationt set When the fault positioning judgment condition is not met and the current time is less than the time windowt set And then, along with the increase of the current moment, the complete cyclic judgment calculation of the complex voltage integral curve is continuously carried out.
6. The method of claim 5 wherein the calculation intervals are selected and the voltage profile is calculated every other calculation interval to obtain a complete and reproducible voltage integral curve, comprising,
before each time of calculating the voltage distribution curve, obtaining the current moment according to the calculation interval time and the initial moment;
the current time is obtained by:
wherein the content of the first and second substances,tis the current time of day and is,t 1 which indicates the starting moment in time of the day,nis shown asnThe sub-cycle calculates a complex voltage integral curve,Δtthe table calculates the interval time.
7. A method for fault location based on along-line voltage nadir calculation according to claim 5 or 6,
after calculating the initial time, the method continuously carries out the operation of the composite voltage integral curveStarting the first judgment after the second calculation, judging the firstnWhen the complex voltage integral curve is calculatedWhether or not it is greater than or equal ton:
If it isThen continue to return to executen+1Calculating a secondary complex voltage integral curve and entering judgment cycle calculation again;
8. The method as claimed in claim 7, wherein the fault location is based on the calculation of the lowest point of the line voltageIf so, it is judged asnWhether a plurality of minimum value points exist in the complex voltage integral curve obtained by the secondary calculation or not comprises,
if there are multiple minimum values in the said integrated curve, define the firstnComplex obtained by sub-calculationThe minimum point of the integration value in the integration curve of the voltage isThe second smallest extreme point of the integral value isJudgment ofAndwhether or not the integrated value of (b) satisfies a determination condition; wherein the content of the first and second substances,
if the current time of calculating the complex voltage integral curve is still in the total calculation time windowt set If so, after a calculation interval time, continuing to execute the step onenCalculating a voltage integration curve for +1 times, and entering a judgment cycle again;
if the current time for calculating the complex voltage integral curve exceeds the total calculation time windowt set If yes, the loop is exited and the fault location is not possible;
9. The method of claim 8 wherein the fault location is based on line voltage nadir calculation,
is judged by the following formulaAndwhether or not the integrated value of (b) satisfies a determination condition:
wherein the content of the first and second substances,representing the minimum point in the complex voltage integral curveThe value of the integral of (a) is,representing the next smallest extreme point in the complex voltage integral curveThe value of the integral of (a) is,and the setting value excluding other minimum value point interference is judged.
10. The method of claim 7 wherein the fault location is based on along-line voltage nadir calculationA process characterized in that saidIf so, judge thatnWhether a plurality of minimum value points exist in the complex voltage integral curve obtained by the secondary calculation or not comprises,
if the plurality of minimum value points do not exist in the complex voltage integral curve, judging whether the complex voltage integral curve has one minimum value point or not; wherein the content of the first and second substances,
when the said complex voltage integral curve has only one minimum value point, the comparison window is compared continuouslyJudging minimum value points of the secondary repetition voltage integral curve;
when there is no minimum point on the complex voltage integral curve:
if the current time of the complex voltage integral curve is still in the total calculation time windowt set If so, after a calculation interval time, continuing to execute the step onenCalculating a voltage integration curve for +1 times, and entering a judgment cycle again;
if the current time for calculating the complex voltage integral curve exceeds the total calculation time windowt set Then the loop is exited and the fault location is deemed to be impossible.
11. The fault location method based on the along-line voltage nadir calculation according to claim 8 or 10,
comparing the window internal continuity byAnd (3) judging minimum value points of the secondary composite voltage integral curve:
12. The method of claim 11 wherein the fault location is based on line voltage nadir calculation,
when continuous in the comparison windowWhen the minimum value point of the secondary complex voltage integral curve meets the condition of a judgment formula, determining the calculated value of the fault distance as;
When continuous in the comparison windowWhen the minimum value point of the secondary complex voltage integral curve does not meet the condition of the judgment formula; wherein the content of the first and second substances,
if the current time of the complex voltage integral curve is still in the total calculation time windowt set If so, after a calculation interval time, continuing to execute the step onenCalculating a voltage integration curve for +1 times, and entering a judgment cycle again;
13. A fault location system based on the calculation of the lowest point of the voltage along the line is characterized by comprising a calculation unit and a determination unit,
the calculating unit is used for performing voltage calculation along the line based on the Bergeron model and performing complex integration in a unit data window to obtain a voltage distribution curve;
and the judging unit is used for calculating a voltage distribution curve and carrying out fault positioning judgment on the line according to the result of the voltage distribution curve.
14. The fault location system based on along-line voltage nadir calculation of claim 13,
the calculating unit is used for determining the starting time of the complex integral calculation according to the full length of the line;
the calculating unit is used for determining a fault phase through the phase selection element based on the initial moment, and calculating a voltage distribution curve of the corresponding fault phase through the Bergeron model according to the voltage and the current obtained through sampling.
15. The fault location system based on along-line voltage nadir calculation of claim 13,
the judging unit is used for calculating the required line step length according to the Bergeron model and determining the sampling rate and the interpolation mode of the line protection device;
calculating the instantaneous value of the step voltage of each point of the circuit from the initial moment to the stop moment in a unit data window of the fault phase by using a Bergeron model;
and (4) calculating the integral value of the voltage of each point through the voltage distribution curve, comparing and judging to obtain the minimum point of the voltage integral value, and confirming the fault point.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211401819.XA CN115508671B (en) | 2022-11-10 | 2022-11-10 | Fault positioning method and system based on line voltage lowest point calculation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211401819.XA CN115508671B (en) | 2022-11-10 | 2022-11-10 | Fault positioning method and system based on line voltage lowest point calculation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115508671A CN115508671A (en) | 2022-12-23 |
CN115508671B true CN115508671B (en) | 2023-03-28 |
Family
ID=84513916
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211401819.XA Active CN115508671B (en) | 2022-11-10 | 2022-11-10 | Fault positioning method and system based on line voltage lowest point calculation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115508671B (en) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114371361A (en) * | 2021-11-01 | 2022-04-19 | 天津大学 | Power transmission line single-end fault location method based on multivariate model fusion analysis |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1804650A (en) * | 2006-01-24 | 2006-07-19 | 天津大学 | Distance measuring method based on Bergeron model for power transmission line |
CN101762774B (en) * | 2009-05-20 | 2011-10-19 | 中国南方电网有限责任公司超高压输电公司 | Method for identifying high voltage direct current transmission line fault location based on genetic algorithm parameter |
CN104092199A (en) * | 2014-07-25 | 2014-10-08 | 国家电网公司 | Line single-phase grounding voltage protection method based on distribution characteristics of voltage amplitudes along line |
CN109030957B (en) * | 2015-05-19 | 2020-12-25 | 江苏理工学院 | Dielectric loss measuring method |
CN114966301A (en) * | 2022-01-26 | 2022-08-30 | 昆明理工大学 | Fault distance measurement method and system for three-terminal hybrid high-voltage direct-current transmission line |
-
2022
- 2022-11-10 CN CN202211401819.XA patent/CN115508671B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114371361A (en) * | 2021-11-01 | 2022-04-19 | 天津大学 | Power transmission line single-end fault location method based on multivariate model fusion analysis |
Also Published As
Publication number | Publication date |
---|---|
CN115508671A (en) | 2022-12-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Korkali et al. | Traveling-wave-based fault-location technique for transmission grids via wide-area synchronized voltage measurements | |
CN108344923B (en) | High-adaptability power transmission line fault location method and system | |
CN109283430B (en) | Power distribution network fault location method based on voltage distribution principle | |
CN112946424B (en) | Method and device for accurately positioning fault | |
Silveira et al. | Transmission line fault location using two-terminal data without time synchronization | |
CN110927510A (en) | Frequency domain method for power transmission line double-end traveling wave fault location | |
CN108535597B (en) | Line model-based single-phase earth fault section positioning method | |
Taheri et al. | Single-end current-based algorithm for fault location in series capacitor compensated transmission lines | |
CN107247215B (en) | Distribution Network Failure population location algorithm based on Multipoint synchronous measurement data | |
CN115128395A (en) | Voltage sag source positioning method and device | |
CN111123027A (en) | Multi-end transmission line fault location method based on WAMS system | |
CN115508671B (en) | Fault positioning method and system based on line voltage lowest point calculation | |
CN107179476B (en) | Distribution network fault distance measurement method | |
CN111929611B (en) | Distribution line ground fault positioning method and system based on zero-mode current difference | |
CN113447758A (en) | Single-phase earth fault distance measuring method for multi-branch current collection circuit of wind power plant | |
CN115792507B (en) | Multiple lightning stroke distinguishing method and system based on monotonicity of short-time window slope | |
CN114878971B (en) | Power distribution network fault point positioning method, device, equipment and medium | |
CN107478941B (en) | Distribution network failure simulated annealing localization method based on Multipoint synchronous measurement data | |
CN113702762B (en) | Distribution network single-phase earth fault distance measurement method utilizing zero sequence information quantity | |
CN114371361A (en) | Power transmission line single-end fault location method based on multivariate model fusion analysis | |
CN114325226A (en) | High-frequency fault positioning method and system for single-end adaptive correction of power transmission line | |
CN114325240A (en) | Fault line identification method based on high-frequency fault information energy evaluation | |
CN114002542A (en) | Power frequency wide area information-based power distribution network fault positioning method and device | |
CN112731053A (en) | High-resistance grounding fault section positioning method for resonance grounding power distribution network | |
CN112505490A (en) | Power distribution network single-phase earth fault line selection method based on mutual difference absolute value sum |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |