CN114152840A - LCC-MMC hybrid direct current transmission line fault distance measurement method and system - Google Patents

LCC-MMC hybrid direct current transmission line fault distance measurement method and system Download PDF

Info

Publication number
CN114152840A
CN114152840A CN202111437314.4A CN202111437314A CN114152840A CN 114152840 A CN114152840 A CN 114152840A CN 202111437314 A CN202111437314 A CN 202111437314A CN 114152840 A CN114152840 A CN 114152840A
Authority
CN
China
Prior art keywords
frequency
signals
calculating
lcc
direct current
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.)
Granted
Application number
CN202111437314.4A
Other languages
Chinese (zh)
Other versions
CN114152840B (en
Inventor
束洪春
田鑫萃
杨竞及
代月
唐玉涛
董俊
单节杉
韩一鸣
曹璞璘
张广斌
安娜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Kunming University of Science and Technology
Original Assignee
Kunming University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Kunming University of Science and Technology filed Critical Kunming University of Science and Technology
Priority to CN202111437314.4A priority Critical patent/CN114152840B/en
Publication of CN114152840A publication Critical patent/CN114152840A/en
Application granted granted Critical
Publication of CN114152840B publication Critical patent/CN114152840B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Locating Faults (AREA)

Abstract

The invention relates to a fault location method and a system for an LCC-MMC hybrid direct-current transmission line, and belongs to the technical field of relay protection of a power system. The method comprises the steps that traveling wave coupling box current signals are respectively collected from a rectification side (LCC) and an inversion side (MMC) of a hybrid direct current power transmission system (LCC-MMC), and ranging signals are respectively constructed; respectively calculating frequency spectrums of ranging signals at two ends, and screening out equal-interval frequency spectrums reflecting fault positions; respectively calculating frequency difference values of frequency spectrums at two ends, calculating a ratio of the frequency difference values at the two ends, and calculating a ranging estimation result k%, (1-k)%, by using the ratio; and respectively calculating double-end check signals, comparing whether the difference of the double-end check signals is smaller than a preset check threshold, if so, determining the fault distance measurement result to be max [ k%, (1-k)% ], and if not, determining the fault distance measurement result to be min [ k%, (1-k)% ]. The invention does not need to identify the properties of the traveling wave head, calibrate the arrival time of the wave, synchronize the clocks at the two ends and is not influenced by the steepness of the traveling wave head, the missing of the wave head, the data defect of the traveling wave and the wave speed.

Description

LCC-MMC hybrid direct current transmission line fault distance measurement method and system
Technical Field
The invention relates to a fault location method and a system for an LCC-MMC hybrid direct-current transmission line, and belongs to the technical field of relay protection of a power system.
Background
The traditional high-voltage direct-current transmission (LCC-HVDC) technology based on the phase control converter tends to be mature, has the advantages of large transmission capacity, low cost, mature operation technology and the like, is not influenced by power transmission distance, and solves the problem of stable power angle in an alternating-current system, so the phase control converter (LCC) is an important means for solving the interconnection between power long-distance transmission and a power grid. However, when the LCC is used as an inverter, there are the following problems: firstly, there is the risk of commutation failure, secondly the transmission voltage level is lower, at the bottom of the harmonic number of times that the interchange side formed, need install the wave filter, thirdly LCC transverter absorbs a large amount of reactive power, requires very high to reactive power compensation arrangement. Compared with the LCC converter, the flexible direct current based on the full-control device IGBT and the PWM modulation technology has no commutation failure and small harmonic wave, can be connected with a weak alternating current and passive system and can be used for black start after power grid faults, has the remarkable characteristics and advantages in many aspects, and becomes an important technical means for solving the problem of large-scale direct current power feed limitation outside a receiving end system receiving area. By integrating the advantages of LCC and MMC (modular multilevel converter), the LCC is used as a rectifying side, the MMC is used as an inverting side, and the formed hybrid high-voltage direct-current transmission system (LCC-MMC) gradually becomes the mainstream mode of direct-current remote transmission engineering.
The fault location technology of the LCC-MMC hybrid direct-current transmission line can be divided into a traveling wave method, a natural frequency method and a fault analysis method in principle. The traveling wave method is characterized in that a traveling wave method is used for calibrating a wave head of a fault traveling wave reaching a measuring end by using a signal processing method, and calculating a fault distance according to a constraint relation between traveling wave reaching time and the fault distance, and the key point is precise calibration of the traveling wave head; the fault analysis method realizes fault positioning by optimizing and solving according to the mathematical relation between the line parameters and the electrical quantities; the fault location based on the natural frequency method mainly utilizes the mathematical relationship between the main frequency component of fault traveling wave and the fault distance to realize fault location, and the location precision is influenced by the accurate extraction of main frequency, frequency spectrum aliasing, frequency spectrum leakage, system boundary conditions and the like, and has not been practically applied.
Disclosure of Invention
The invention aims to solve the technical problem of providing a fault location method and a fault location system for an LCC-MMC hybrid direct current transmission line, which are used for solving the problems in the prior art in the aspect of fault location.
The technical scheme of the invention is as follows: a fault location method and a system for an LCC-MMC hybrid direct current transmission line can construct a location equation by respectively utilizing data of a rectifying side and data of an inverting side, and compared with a traveling wave location method, the method does not need to accurately detect and calibrate a wave head, double-end communication does not need synchronous time synchronization, and the anti-noise interference capability and the anti-transition resistance capability are stronger; compared with the natural frequency method, the mixed direct current system has different system physical boundaries at two ends, inconsistent wave velocity, no error of system control response to the ranging result, no influence of the accuracy of the natural frequency main frequency value, and elimination of frequency spectrum aliasing caused by fault point refraction and reflection.
The method comprises the following specific steps:
step 1: and respectively acquiring current signals of the traveling wave coupling box from the rectifying side and the inverting side of the hybrid direct-current power transmission system, and respectively constructing ranging signals.
Step 2: and respectively calculating the frequency spectrums of the ranging signals at the two ends, and screening out the equal-interval frequency spectrums reflecting the fault positions.
Step 3: and respectively calculating frequency difference values of frequency spectrums at two ends, calculating the ratio of the frequency difference values at the two ends, and calculating the distance measurement estimation result k%, (1-k)%, by using the ratio.
Step 4: and respectively calculating double-end check signals, comparing whether the difference of the double-end check signals is smaller than a preset check threshold, if so, determining the fault distance measurement result to be max [ k%, (1-k)% ], and if not, determining the fault distance measurement result to be min [ k%, (1-k)% ].
The Step1 is specifically as follows:
step1.1: and respectively acquiring current signals of the traveling wave coupling box from the rectifying side and the inverting side of the hybrid direct-current power transmission system.
Step1.2: decoupling a direct current line through a Karenbauer transformation matrix, and solving polar space modulus current signals of an rectifying side and an inverting side as ranging signals.
The Step2 is specifically as follows:
step2.1: and respectively calculating frequency spectrums of ranging signals on the rectifying side and the inverting side.
Step2.2: and judging whether the frequency spectrums of the frequency domain signals at the two ends are at equal intervals, if not, executing Step2.3, and if so, respectively calculating the frequency difference of the signals at the two ends.
Step2.3: and respectively calculating power transformation signals of the ranging signals at the two ends of the direct current transmission line, and returning to Step2.1.
The power transformation refers to an odd power.
The Step3 is specifically as follows:
step3.1: frequency values of spectral peak positions of the double-end equal interval frequency spectrums are respectively extracted, any two adjacent frequency values are subjected to difference, and frequency differences corresponding to the frequency spectrums at the two ends are obtained.
Step3.2: the ratio omega of the double-end frequency difference is obtained according to the formula (1)MNM denotes a rectification side of the double-ended dc system, and N denotes an inversion side of the double-ended dc system.
Figure RE-RE-RE-GDA0003495101950000021
In the formula (1), Δ ωM、ΔωNRespectively representing the angular frequency difference, delta f calculated by the M terminal and the N terminalM、ΔfNRespectively representing the frequency difference obtained by the calculation of the M end and the N end, wherein the double-end frequency difference and the angular frequency difference respectively satisfy the relational expressions of the expressions (2) and (3).
Figure RE-RE-RE-GDA0003495101950000022
Figure RE-RE-RE-GDA0003495101950000023
Step3.3: using the ratio omega of the frequency differences of both endsMNAnd (4) obtaining a ranging estimation result k% by the constructed ranging formula (4).
Figure RE-RE-RE-GDA0003495101950000031
In the formula (4), xfThe distance between the fault point and the end M is shown, k% represents the percentage of the distance between the fault point and the end M to the total length of the line, and l represents the total length of the direct current line.
The Step3.2 can also be used for calculating the ratio omega of the double-end frequency difference according to the formula (5)NM。
Figure RE-RE-RE-GDA0003495101950000032
Figure RE-RE-RE-GDA0003495101950000033
The Step3.3 can also be realized by utilizing the ratio omega of double-end frequency differenceNMThe constructed distance measurement formula (6) obtains the distance measurement result (1-k)%.
The Step4 is specifically as follows:
step4.1: the absolute value of the difference between the ratio of the two-terminal frequency differences and 1 is calculated.
Step4.2: judging whether the absolute value of the difference between the double-end frequency difference ratio and 1 is larger than a preset threshold value, if so, determining that the fault distance measurement result is k%, otherwise, executing Step4.3
Step4.3: respectively calculating the change rate K between three adjacent sampling points of the wave head detected by the first two endsMAnd KN。
Step4.4: comparison KM-KNAnd whether the value is smaller than a preset check threshold value or not, if so, executing Step4.5, and otherwise, executing Step4.6.
Step4.5: the fault location result is max [ k%, (1-k)% ].
Ste4.6: the fault location result is min [ k%, (1-k)% ].
A LCC-MMC mixed direct current transmission line fault distance measuring system comprises:
and the signal acquisition module is used for respectively acquiring current signals of the traveling wave coupling box from the rectifying side and the inverting side of the hybrid direct-current power transmission system.
And the signal preprocessing module is used for constructing a ranging signal.
And the frequency domain signal analysis module is used for respectively calculating the frequency spectrums of the ranging signals at the two ends and screening out the equal interval frequency spectrums reflecting the fault positions.
And the frequency domain signal processing module is used for calculating the equal interval frequency difference value of the frequency domain signal and calculating the ratio of the frequency difference values at two ends.
And the fault distance measurement module is used for calculating and recording fault distance measurement results.
And the signal checking module is used for calculating the checking signal.
The signal acquisition module includes:
and the voltage and current transmitting unit is used for converting the voltage and current signals on the secondary side of the transformer into signals acquired by the traveling wave device A/D.
And the analog quantity acquisition unit is used for converting the voltage and current analog quantity signals into digital signals.
And the data storage unit is used for naming the wave recording data files according to the time stamps and storing the wave recording data files in the local memory.
And the abrupt change starting unit is used for judging whether the waveform abrupt change is larger than a set starting threshold value or not, and storing the voltage and current signals into a recording data file if the waveform abrupt change is larger than the set starting threshold value.
The signal preprocessing module comprises:
and the numerical value calculation unit 1 is used for decoupling the direct current line by using a Karenbauer transformation matrix to obtain a polar space current modulus traveling wave transient signal.
And the numerical value calculating unit 2 is used for calculating power transformation signals of the ranging signals at the two ends of the direct current transmission line respectively.
The frequency domain signal analysis module comprises:
and the time-frequency transformation unit is used for extracting frequency domain signals of the related signals at the two ends of the direct current transmission line.
And the frequency spectrum calculating unit is used for calculating the frequency spectrum of the frequency domain signals at the two ends.
And the frequency spectrum judging unit is used for judging whether the frequency spectrums of the frequency domain signals at the two ends are at equal intervals.
The frequency domain signal processing module includes:
and the frequency difference calculating unit is used for extracting the frequency values of the spectral peak positions of the double-end equally-spaced frequency spectrums, and performing difference on any two adjacent frequency values to obtain the frequency differences corresponding to the frequency spectrums at the two ends.
A ratio calculation unit for calculating a ratio omega of the two-terminal frequency differencesMN、ΩNMM denotes one end of the double-ended dc system, and N denotes the other end of the double-ended dc system.
The signal verification module includes:
and the checking signal calculating unit 1 is used for calculating the absolute value of the difference between the double-ended frequency difference ratio and 1.
A check signal calculating unit 2 for calculating the change rate K between three adjacent sampling points of the first detected wave head at both endsMAnd KN。
And the check signal comparison unit 1 is used for comparing whether the absolute value of the difference between the double-end frequency difference ratio and 1 is greater than a preset threshold value.
A check signal comparison unit 2 for comparing KM-KNWhether or not toIs less than the preset check threshold.
The invention has the beneficial effects that:
1. the method is not influenced by different system boundary conditions of the LCC side and the MMC side, and avoids errors caused by unreliable traveling wave head identification and inaccurate wave arrival time calibration. .
2. According to the invention, a ranging equation is constructed by using the frequency difference ratio obtained at the LCC side and the MMC side, the ranging precision is not influenced by wave head distortion, waveform defect and head wave head loss, the influence of wave velocity attenuation is mathematically reduced, and the method does not depend on double-end synchronous time setting and has better robustness on weak fault modes such as high-resistance fault and the like.
3. The invention can be realized by utilizing the traveling wave signals collected by the traveling wave distance measuring device and the signals collected by the fault recorder, is not influenced by interference signals, channel noise and channel synchronism, is used as supplement and assistance of the existing traveling wave distance measuring algorithm, reduces the distance measuring result to the number of the tower, does not participate in distance measuring calculation of the line length, and has the precision not influenced by the actual change of the line length.
Drawings
FIG. 1 is a schematic diagram of a system according to an embodiment;
FIG. 2 is a schematic diagram of a system boundary on the rectifying side according to the first embodiment;
FIG. 3 is a schematic diagram of a current ranging signal obtained by a corresponding M-port traveling-wave coupling box according to the first embodiment;
FIG. 4 is a schematic diagram of a current ranging signal obtained by a corresponding N-terminal traveling-wave coupling box according to the first embodiment;
FIG. 5 is a M-terminal spectrum according to an embodiment;
FIG. 6 shows an N-terminal spectrum according to an embodiment;
fig. 7 is a schematic flow chart of a fault location method for an LCC-MMC hybrid dc transmission line according to a second embodiment of the present invention;
FIG. 8 is a flowchart illustrating the detailed steps of Step1 according to the second embodiment;
FIG. 9 is a flowchart illustrating the detailed steps of Step2 according to the second embodiment;
FIG. 10 is a flowchart illustrating the detailed steps of Step3 according to the second embodiment;
FIG. 11 is a flowchart illustrating the detailed steps of Step4 according to the second embodiment;
fig. 12 is a functional block diagram of a system corresponding to the LCC-MMC hybrid dc transmission line fault ranging of the second embodiment.
Detailed Description
The invention is further described with reference to the following drawings and detailed description.
As shown in FIG. 1, the +/-800 kV true bipolar double-ended hybrid direct current transmission system (LCC-MMC) is built under a PSCAD/EMTDC environment. The rectification side is an M end, an LCC converter is adopted, the inversion side is an N end, and an MMC converter is adopted. Each pole of the rectifying side is formed by connecting double 12 pulse current converters in series. The direct current transmission line is a six-split conductor, the line length is 1500km, and 400mH smoothing reactors are arranged on two sides of the line.
Fig. 2 is a schematic diagram of the physical boundary formed by the smoothing reactor and the dc filter. As shown in fig. 2, the characteristic harmonic order on the dc side is Kp, where p is the pulse number of the inverter and K is any positive integer, so that 12/24/36 triple-tuned filters are selected as the dc filter. The parameters are as follows: l is 400mH, L1＝39.09mH、L2＝26.06mH、L3＝19.545mH、L4＝34.75mH、 C1＝0.9μF、C2＝0.9μF、C3＝1.8μF、C40.675 μ F. In the figure B1～B4、D1、D2The lightning arresters are respectively a direct current filter, a smoothing reactor and a direct current bus. The smoothing reactor can be regarded as an open circuit within the traveling wave frequency range. The boundary of the inversion side is the smoothing reactor.
Example 1: suppose that a unipolar metallic short fault occurs at 400km from the M end, with a sampling rate of 200 kHz. The specific implementation steps are as follows:
step 1: collecting current signals of a traveling wave coupling box from a rectifying side (LCC) and an inverting side (MMC) of a hybrid direct current power transmission system (LCC-MMC) respectively, and constructing distance measurement signals respectively.
Step1.1: and respectively acquiring current signals of the traveling wave coupling box from the rectifying side and the inverting side of the hybrid direct-current power transmission system.
Step1.2: decoupling a direct current line through a Karenbauer transformation matrix, and solving polar space modulus current signals on the rectifying side and the inverting side as ranging signals, as shown in fig. 3 and 4.
Step 2: and respectively carrying out Fourier transform on the ranging signals at the two ends to obtain corresponding frequency spectrums, and screening out the equidistant frequency spectrums reflecting the fault positions.
Step2.1: the frequency spectrums of the ranging signals on the rectifying side and the inverting side are calculated respectively as shown in fig. 5 and 6.
Step2.2: and judging whether the frequency spectrums of the frequency domain signals at the two ends are at equal intervals, if not, executing Step2.3, and if so, respectively calculating the frequency difference of the signals at the two ends.
Step2.3: and respectively calculating power transformation signals of the ranging signals at the two ends of the direct current transmission line, and returning to Step2.1.
The power transformation refers to an odd power.
Step 3: respectively calculating frequency difference values of frequency spectrums at two ends, calculating a ratio of the frequency difference values at the two ends, and calculating a ranging estimation result k%, (1-k)%, by using the ratio, specifically comprising:
step3.1: frequency values ω of spectral peak positions of the double-ended equally spaced spectrum are extracted respectively, as indicated by circles in fig. 4 and 5.
And (5) according to the formulas (7) and (8), carrying out difference on any two adjacent frequency values to obtain the frequency difference corresponding to the frequency spectrums at the two ends.
ΔωM＝ωM(k2+1)-ωM(k2) (7)
ΔωN＝ωN(k2+1)-ωN(k2) (8)
The calculation result is obtained by calculation according to the formulas (2), (3), (7) and (8): Δ fM＝372.5Hz，ΔfN＝135.45Hz。 k 21,2, … (length (ω) -1), length (ω) represents the length of the sequence of spectral peak frequencies ω.
Step3.2: the ratio omega of the double-end frequency difference is obtained according to the formula (1)MNM denotes a rectification side (LCC) of the double-ended DC system, and N denotes an inverter side of the double-ended DC system(MMC)。
Figure RE-RE-RE-GDA0003495101950000061
In the formula (1), Δ ωM、ΔωNRespectively representing the angular frequency difference, delta f calculated by the M terminal and the N terminalM、ΔfNRespectively representing the frequency difference obtained by the calculation of the M end and the N end, wherein the double-end frequency difference and the angular frequency difference respectively satisfy the relational expressions of the expressions (2) and (3).
Figure RE-RE-RE-GDA0003495101950000062
Figure RE-RE-RE-GDA0003495101950000063
The ratio omega of the double-end frequency difference is obtained according to the formula (1)MN＝2.75。
Step3.3 uses the ratio omega of the two-terminal frequency differencesMNAnd (4) obtaining a ranging estimation result k% by the constructed ranging formula (4).
Figure RE-RE-RE-GDA0003495101950000064
In the formula (4), xfThe distance between the fault point and the end M is shown, k% represents the percentage of the distance between the fault point and the end M to the total length of the line, and l represents the total length of the direct current line. k% ═ 26.47%.
It is also possible to pass the ratio Ω of the two-terminal frequency differencesNMTo obtain a fault location result, specifically:
omega calculated from the formula (5)NMUsing the ratio omega of the frequency differences of the two endsNMThe constructed distance measurement formula (6) gave a distance measurement result of (1-k)%, 73.53%.
Figure RE-RE-RE-GDA0003495101950000071
Figure RE-RE-RE-GDA0003495101950000072
Step 4: and respectively calculating double-end check signals, comparing whether the difference of the double-end check signals is smaller than a preset check threshold, if so, determining the fault distance measurement result to be max [ k%, (1-k)% ], and if not, determining the fault distance measurement result to be min [ k%, (1-k)% ].
Step4.1: calculating the absolute value of the difference between the ratio of the two end frequency differences and 1, abs (omega)MN-1)＝1.75。
Step4.2: and judging whether the absolute value of the difference between the double-end frequency difference ratio and 1 is larger than a preset threshold value, setting the absolute value to be 0.05, and if so, setting the fault distance measurement result to be that k% is 26.47%.
Example 2: as shown in fig. 7, a method for locating a fault of an LCC-MMC hybrid dc transmission line includes the following steps:
step 1: collecting current signals of a traveling wave coupling box from a rectifying side (LCC) and an inverting side (MMC) of a hybrid direct current power transmission system (LCC-MMC) respectively, and constructing distance measurement signals respectively.
Step 2: and respectively calculating the frequency spectrums of the ranging signals at the two ends, and screening out the equal-interval frequency spectrums reflecting the fault positions. .
Step 3: and respectively calculating frequency difference values of frequency spectrums at two ends, calculating the ratio of the frequency difference values at the two ends, and calculating the distance measurement estimation result k%, (1-k)%, by using the ratio.
Step 4: and respectively calculating double-end check signals, comparing whether the difference of the double-end check signals is smaller than a preset check threshold, if so, determining the fault distance measurement result to be max [ k%, (1-k)% ], and if not, determining the fault distance measurement result to be min [ k%, (1-k)% ].
As shown in fig. 8, Step1 is specifically:
step1.1: and respectively acquiring current signals of the traveling wave coupling box from the rectifying side and the inverting side of the hybrid direct-current power transmission system.
Step1.2: decoupling a direct current line through a Karenbauer transformation matrix, and solving polar space modulus current signals of an rectifying side and an inverting side as ranging signals.
As shown in fig. 9, Step2 is specifically:
step2.1: and respectively calculating frequency spectrums of ranging signals on the rectifying side and the inverting side.
Step2.2: and judging whether the frequency spectrums of the frequency domain signals at the two ends are at equal intervals, if not, executing Step2.3, and if so, respectively calculating the frequency difference of the signals at the two ends.
Step2.3: and respectively calculating power transformation signals of the ranging signals at the two ends of the direct current transmission line, and returning to Step2.1.
The power transformation refers to an odd power.
As shown in fig. 10, Step3 is specifically:
step3.1: frequency values of spectral peak positions of the double-end equal interval frequency spectrums are respectively extracted, any two adjacent frequency values are subjected to difference, and frequency differences corresponding to the frequency spectrums at the two ends are obtained.
Step3.2: the ratio omega of the double-end frequency difference is obtained according to the formula (1)MNM denotes a rectification side (LCC) of the double-ended dc system, and N denotes an inversion side (MMC) of the double-ended dc system.
Figure RE-RE-RE-GDA0003495101950000081
In the formula (1), Δ ωM、ΔωNRespectively representing the angular frequency difference, delta f calculated by the M terminal and the N terminalM、ΔfNRespectively representing the frequency difference obtained by the calculation of the M end and the N end, wherein the double-end frequency difference and the angular frequency difference respectively satisfy the relational expressions of the expressions (2) and (3).
Figure RE-RE-RE-GDA0003495101950000082
Figure RE-RE-RE-GDA0003495101950000083
Step3.3: using the ratio omega of the frequency differences of both endsMNAnd (4) obtaining a ranging estimation result k% by the constructed ranging formula (4).
Figure RE-RE-RE-GDA0003495101950000084
In the formula (4), xfThe distance between the fault point and the end M is shown, k% represents the percentage of the distance between the fault point and the end M to the total length of the line, and l represents the total length of the direct current line.
The ratio of the double-end frequency difference is obtained in Step3.2, and omega calculated by the formula (5) is also includedNMUsing the ratio omega of the frequency differences of the two endsNMThe constructed distance measurement formula (6) obtains the distance measurement result (1-k)%.
Figure RE-RE-RE-GDA0003495101950000085
Figure RE-RE-RE-GDA0003495101950000086
As shown in fig. 11, Step4 is specifically:
step4.1: the absolute value of the difference between the ratio of the two-terminal frequency differences and 1 is calculated.
Step4.2: judging whether the absolute value of the difference between the double-end frequency difference ratio and 1 is larger than a preset threshold value, if so, determining that the fault distance measurement result is k%, otherwise, executing Step4.3
Step4.3: respectively calculating the change rate K between three adjacent sampling points of the wave head detected by the first two endsMAnd KN。
Step4.4: comparison KM-KNAnd whether the value is smaller than a preset check threshold value or not, if so, executing Step4.5, and otherwise, executing Step4.6.
Step4.5, the failure ranging result is max [ k%, (1-k)% ].
Step4.6, the fault location result is min [ k%, (1-k)% ].
As shown in fig. 12, an LCC-MMC hybrid dc transmission line fault location system includes:
the signal acquisition module 201 is configured to acquire traveling wave coupling box current signals from a rectification side (LCC) and an inversion side (MMC) of a hybrid direct current power transmission system (LCC-MMC), respectively.
A signal preprocessing module 202 for constructing ranging signals.
And the frequency domain signal analysis module 203 is used for respectively calculating frequency spectrums of the ranging signals at the two ends and screening out the equal interval frequency spectrums reflecting the fault positions.
The frequency domain signal processing module 204 is configured to calculate an equally spaced frequency difference value of the frequency domain signal, and calculate a ratio of the two end frequency difference values.
And the fault distance measurement module 205 is used for calculating and recording a fault distance measurement result.
And a signal checking module 206 for calculating a checking signal.
The signal acquisition module 201 includes:
and the voltage current transmission unit 2011 is used for converting a voltage current signal on the secondary side of the transformer into a signal acquired by the traveling wave device A/D.
The analog acquisition unit 2012 is configured to convert the voltage-current analog signal into a digital signal.
And the data storage unit 2013 is used for naming the wave recording data file according to the time mark and storing the wave recording data file in a local memory.
And the abrupt change starting unit 2014 is used for judging whether the waveform abrupt change is larger than a set starting threshold value or not, and if so, storing the voltage and current signals into a wave recording data file.
The signal preprocessing module 202 includes:
the numerical calculation unit 2021 is configured to decouple the dc line by using a Karenbauer transformation matrix to obtain a polar space current modulus traveling wave transient signal.
A numerical calculation unit 2022 for calculating power transformation signals for respectively calculating distance measurement signals at two ends of the DC transmission line
The frequency domain signal analysis module 203 comprises:
the time-frequency transform unit 2031 is configured to extract frequency domain signals of the signals associated with the two ends of the dc transmission line.
A spectrum calculating unit 2032 for calculating the spectrum of the two-end frequency domain signal.
Spectrum determining section 2033 configured to determine whether or not the frequency spectrums of the two-end frequency domain signals are at equal intervals.
The frequency domain signal processing module 204 includes:
the frequency difference calculating unit 2041 is configured to extract frequency values at spectral peak positions of the double-end equally-spaced frequency spectrums, and perform a difference between any two adjacent frequency values to obtain equal frequency differences corresponding to the frequency spectrums at the two ends.
A ratio calculation unit 2042 for calculating a ratio Ω of the two-terminal frequency differenceMN、ΩNMM denotes one end of the double-ended dc system, and N denotes the other end of the double-ended dc system.
The signal verification module 206 includes:
a check signal calculation unit 2061 for calculating an absolute value of a difference between the double-ended frequency difference ratio and 1.
A check signal calculation unit 2062 for calculating the change rate K between three adjacent sampling points of the wave head detected by the first two endsMAnd KN。
A check signal comparing unit 2063 for comparing whether the absolute value of the difference between the double-ended frequency difference ratio and 1 is greater than a preset threshold value
A check signal comparison unit 2064 for comparing KM-KNWhether the value is less than a preset check threshold value.
While the present invention has been described in detail with reference to the embodiments shown in the drawings, the present invention is not limited to the embodiments, and various changes can be made without departing from the spirit and scope of the present invention.

Claims (12)

1. A fault location method for an LCC-MMC hybrid direct current transmission line is characterized by comprising the following steps:
step 1: respectively acquiring current signals of a traveling wave coupling box from a rectifying side and an inverting side of the hybrid direct-current power transmission system, and respectively constructing ranging signals;
step 2: respectively calculating frequency spectrums of ranging signals at two ends, and screening out equal-interval frequency spectrums reflecting fault positions;
step 3: respectively calculating frequency difference values of frequency spectrums at two ends, calculating a ratio of the frequency difference values at the two ends, and calculating a ranging estimation result k%, (1-k)%, by using the ratio;
step 4: and respectively calculating double-end check signals, comparing whether the difference of the double-end check signals is smaller than a preset check threshold, if so, determining the fault distance measurement result to be max [ k%, (1-k)% ], and if not, determining the fault distance measurement result to be min [ k%, (1-k)% ].
2. The LCC-MMC hybrid direct current transmission line fault location method of claim 1, wherein Step1 specifically is:
step1.1: respectively collecting current signals of a traveling wave coupling box from a rectifying side and an inverting side of the hybrid direct-current power transmission system;
step1.2: decoupling a direct current line through a Karenbauer transformation matrix, and solving polar space modulus current signals of an rectifying side and an inverting side as ranging signals.
3. The LCC-MMC hybrid direct current transmission line fault location method of claim 1, wherein Step2 specifically is:
step2.1: respectively calculating frequency spectrums of ranging signals on a rectifying side and a inverting side;
step2.2: judging whether the frequency spectrums of the frequency domain signals at the two ends are at equal intervals, if not, executing Step2.3, and if so, respectively calculating the frequency difference of the signals at the two ends;
step2.3: and respectively calculating power transformation signals of the ranging signals at the two ends of the direct current transmission line, and returning to Step2.1.
4. The LCC-MMC hybrid direct current transmission line fault location method of claim 1, wherein Step3 specifically is:
step3.1: respectively extracting frequency values of spectral peak positions of double-end equal-interval frequency spectrums, and carrying out difference on any two adjacent frequency values to obtain corresponding equal frequency differences of the frequency spectrums at the two ends;
step3.2: according to formula (1)Finding the ratio omega of the frequency differences of both endsMNM represents the rectification side of the double-end direct-current system, and N represents the inversion side of the double-end direct-current system;
Figure FDA0003381877400000011
in the formula (1), Δ ωM、ΔωNRespectively representing the angular frequency difference, delta f calculated by the M terminal and the N terminalM、ΔfNRespectively representing the frequency difference obtained by the calculation of the M end and the N end, wherein the double-end frequency difference and the angular frequency difference respectively satisfy the relational expressions of the expressions (2) and (3);
Figure FDA0003381877400000012
Figure FDA0003381877400000021
step3.3: using the ratio omega of the frequency differences of both endsMNThe constructed distance measurement formula (4) obtains a distance measurement estimation result k%;
Figure FDA0003381877400000022
in the formula (4), xfThe distance between the fault point and the end M is shown, k% represents the percentage of the distance between the fault point and the end M to the total length of the line, and l represents the total length of the direct current line.
5. The LCC-MMC hybrid direct current transmission line fault location method of claim 4, characterized in that: the Step3.2 can also be used for calculating the ratio omega of the double-end frequency difference according to the formula (5)NM；
Figure FDA0003381877400000023
Figure FDA0003381877400000024
The Step3.3 can also be realized by utilizing the ratio omega of double-end frequency differenceNMThe constructed distance measurement formula (6) obtains the distance measurement result (1-k)%.
6. The LCC-MMC hybrid direct current transmission line fault location method of claim 1, wherein Step4 specifically is:
step4.1: calculating the absolute value of the difference between the ratio of the two end frequency differences and 1;
step4.2: judging whether the absolute value of the difference between the double-end frequency difference ratio and 1 is larger than a preset threshold value, if so, determining that the fault distance measurement result is k%, otherwise, executing Step4.3
Step4.3: respectively calculating the change rate K between three adjacent sampling points of the wave head detected by the first two endsMAnd KN；
Step4.4: comparison KM-KNWhether the current value is smaller than a preset check threshold value or not is judged, if yes, Step4.5 is executed, and if not, Step4.6 is executed;
step4.5: the fault location result is max [ k%, (1-k)% ];
ste4.6: the fault location result is min [ k%, (1-k)% ].
7. The utility model provides a mixed direct current transmission line trouble ranging system of LCC-MMC which characterized in that includes:
the signal acquisition module is used for respectively acquiring current signals of the traveling wave coupling box from the rectifying side and the inverting side of the hybrid direct-current power transmission system;
the signal preprocessing module is used for constructing a ranging signal;
the frequency domain signal analysis module is used for respectively calculating frequency spectrums of the ranging signals at the two ends and screening out the equal interval frequency spectrums reflecting the fault positions;
the frequency domain signal processing module is used for calculating the equal interval frequency difference value of the frequency domain signal and calculating the ratio of the frequency difference values at the two ends;
the fault distance measurement module is used for calculating and recording a fault distance measurement result;
and the signal checking module is used for calculating the checking signal.
8. The LCC-MMC hybrid direct current transmission line fault ranging system of claim 7, wherein said signal acquisition module comprises:
the voltage and current transmission unit is used for converting a voltage and current signal on the secondary side of the transformer into a signal acquired by the traveling wave device A/D;
the analog quantity acquisition unit is used for converting the voltage and current analog quantity signals into digital signals;
the data storage unit is used for naming the wave recording data files according to the time stamps and storing the wave recording data files in a local memory;
and the abrupt change starting unit is used for judging whether the waveform abrupt change is larger than a set starting threshold value or not, and storing the voltage and current signals into a recording data file if the waveform abrupt change is larger than the set starting threshold value.
9. The LCC-MMC hybrid direct current transmission line fault ranging system of claim 7, wherein said signal pre-processing module comprises:
the numerical calculation unit 1 is used for decoupling the direct current line by using a Karenbauer transformation matrix and solving polar space current modulus traveling wave transient signals;
and the numerical value calculating unit 2 is used for calculating power transformation signals of the ranging signals at the two ends of the direct current transmission line respectively.
10. The LCC-MMC hybrid direct current transmission line fault ranging system of claim 7, wherein said frequency domain signal analysis module comprises:
the time-frequency transformation unit is used for extracting frequency domain signals of related signals at two ends of the direct current transmission line;
the frequency spectrum calculating unit is used for calculating the frequency spectrum of the frequency domain signals at the two ends;
and the frequency spectrum judging unit is used for judging whether the frequency spectrums of the frequency domain signals at the two ends are at equal intervals.
11. The LCC-MMC hybrid direct current transmission line fault ranging system of claim 7, wherein said frequency domain signal processing module comprises:
the frequency difference calculating unit is used for extracting frequency values of spectral peak positions of the double-end equal-interval frequency spectrums, and performing difference on any two adjacent frequency values to obtain the corresponding equal frequency differences of the frequency spectrums at the two ends;
a ratio calculation unit for calculating a ratio omega of the two-terminal frequency differencesMN、ΩNMM denotes one end of the double-ended dc system, and N denotes the other end of the double-ended dc system.
12. The LCC-MMC hybrid direct current transmission line fault ranging system of claim 7, wherein said signal verification module comprises:
the check signal calculating unit 1 is used for calculating the absolute value of the difference between the double-end frequency difference ratio and 1;
a check signal calculating unit 2 for calculating the change rate K between three adjacent sampling points of the first detected wave head at both endsMAnd KN；
The check signal comparison unit 1 is used for comparing whether the absolute value of the difference between the double-end frequency difference ratio and 1 is greater than a preset threshold value or not;
a check signal comparison unit 2 for comparing KM-KNWhether the value is less than a preset check threshold value.
CN202111437314.4A 2021-11-29 2021-11-29 LCC-MMC hybrid direct current transmission line fault distance measurement method and system Active CN114152840B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111437314.4A CN114152840B (en) 2021-11-29 2021-11-29 LCC-MMC hybrid direct current transmission line fault distance measurement method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111437314.4A CN114152840B (en) 2021-11-29 2021-11-29 LCC-MMC hybrid direct current transmission line fault distance measurement method and system

Publications (2)

Publication Number Publication Date
CN114152840A true CN114152840A (en) 2022-03-08
CN114152840B CN114152840B (en) 2022-10-25

Family

ID=80455261

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111437314.4A Active CN114152840B (en) 2021-11-29 2021-11-29 LCC-MMC hybrid direct current transmission line fault distance measurement method and system

Country Status (1)

Country Link
CN (1) CN114152840B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114518506A (en) * 2021-11-05 2022-05-20 昆明理工大学 Method and system for detecting fault of direct-current transmission line
CN115267428A (en) * 2022-07-20 2022-11-01 福州大学 LCC-MMC single-pole grounding fault positioning method based on VMD-ET feature selection
CN117388642A (en) * 2023-12-11 2024-01-12 昆明理工大学 Grounding electrode line fault location method and system based on single-ended frequency difference

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105738759A (en) * 2014-12-12 2016-07-06 国家电网公司 Transient recording data-based direct-current power transmission line fault locating method
CN106019076A (en) * 2016-05-16 2016-10-12 西安理工大学 Fault distance detecting method for high-voltage DC transmission line
CN110361627A (en) * 2019-06-24 2019-10-22 昆明理工大学 A kind of single-ended traveling wave fault location method based on MMC-HVDC
CN112083272A (en) * 2020-08-18 2020-12-15 昆明理工大学 10kV cable fault location method based on sheath current traveling wave natural frequency difference method
CN112083271A (en) * 2020-08-18 2020-12-15 昆明理工大学 10kV cable online distance measurement method based on sheath current traveling wave time-frequency composite analysis
CN113471940A (en) * 2021-07-02 2021-10-01 天津大学 MMC-HVDC protection method based on voltage high-low frequency ratio difference
CN113589101A (en) * 2021-07-30 2021-11-02 华北水利水电大学 Fault location method and system suitable for direct-current power transmission system

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105738759A (en) * 2014-12-12 2016-07-06 国家电网公司 Transient recording data-based direct-current power transmission line fault locating method
CN106019076A (en) * 2016-05-16 2016-10-12 西安理工大学 Fault distance detecting method for high-voltage DC transmission line
CN110361627A (en) * 2019-06-24 2019-10-22 昆明理工大学 A kind of single-ended traveling wave fault location method based on MMC-HVDC
CN112083272A (en) * 2020-08-18 2020-12-15 昆明理工大学 10kV cable fault location method based on sheath current traveling wave natural frequency difference method
CN112083271A (en) * 2020-08-18 2020-12-15 昆明理工大学 10kV cable online distance measurement method based on sheath current traveling wave time-frequency composite analysis
CN113471940A (en) * 2021-07-02 2021-10-01 天津大学 MMC-HVDC protection method based on voltage high-low frequency ratio difference
CN113589101A (en) * 2021-07-30 2021-11-02 华北水利水电大学 Fault location method and system suitable for direct-current power transmission system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
WEI CHEN ET AL.: "Travelling wave fault location approach for MMC-HVDC transmission line based on frequency modification algorithm", 《INTERNATIONAL JOURNAL OF ELECTRICAL POWER & ENERGY SYSTEMS》 *
陈仕龙等: "一种特高压直流输电线路神经网络双端故障测距新方法", 《电工技术学报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114518506A (en) * 2021-11-05 2022-05-20 昆明理工大学 Method and system for detecting fault of direct-current transmission line
CN114518506B (en) * 2021-11-05 2022-07-22 昆明理工大学 Method and system for detecting fault of direct-current transmission line
CN115267428A (en) * 2022-07-20 2022-11-01 福州大学 LCC-MMC single-pole grounding fault positioning method based on VMD-ET feature selection
CN115267428B (en) * 2022-07-20 2024-06-11 福州大学 LCC-MMC monopole ground fault positioning method based on VMD-ET feature selection
CN117388642A (en) * 2023-12-11 2024-01-12 昆明理工大学 Grounding electrode line fault location method and system based on single-ended frequency difference
CN117388642B (en) * 2023-12-11 2024-04-09 昆明理工大学 Grounding electrode line fault location method and system based on single-ended frequency difference

Also Published As

Publication number Publication date
CN114152840B (en) 2022-10-25

Similar Documents

Publication Publication Date Title
CN114152840B (en) LCC-MMC hybrid direct current transmission line fault distance measurement method and system
CN108107321B (en) Fault waveform comparison method for power system
CN104597376B (en) A kind of high voltage direct current transmission line fault location method considering actual measurement velocity of wave
CN108344923B (en) High-adaptability power transmission line fault location method and system
CN110687399B (en) Method for judging waveform fault starting time of power distribution network fault indicator
CN102967779B (en) Identifying method of distribution parameters of transmission line
CN115015686A (en) LCC-VSC hybrid high-voltage direct-current transmission line fault distance measurement method and system
CN116338525B (en) Wind power alternating current outgoing line fault location method and system
CN112986775B (en) Active calibration method for partial discharge pulse phase
CN113945806B (en) Single-ended fault location method and system for flexible direct current transmission line
CN105785124A (en) Method for measuring harmonics and interharmonics of electric power system through spectrum estimation and cross correlation
CN115097253A (en) MMC-HVDC direct current transmission line fault distance measurement method and system
CN108414838A (en) A kind of inverter parallel system line impedance measurement method
CN114578184B (en) Direct-current transmission line double-end traveling wave frequency difference ratio fault location method and system
CN115047284A (en) Fault distance measuring method and system for high-voltage direct-current transmission line
CN115015687A (en) Four-end annular flexible direct-current power grid fault location method and system
CN114200246B (en) Fault location method and system for direct-current transmission line
CN114966194A (en) ADC-based low-frequency metering three-phase electric energy meter
CN111999597B (en) Traveling wave fault positioning device of hybrid power transmission line
CN114460412B (en) Direct-current transmission line single-end fault location method and system
CN114002559B (en) Flexible direct current transmission line traveling wave double-end distance measurement method and system
CN108152583B (en) Method and device for separating fundamental wave and harmonic component of leakage current
CN114184896B (en) Single-ended fault location method and system for hybrid direct-current transmission line
CN110018401B (en) Distribution line single-phase earth fault positioning method
CN113702762A (en) Power distribution network single-phase earth fault distance measurement method using zero sequence information quantity

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