[go: up one dir, main page]

CN112564057B - Adaptive Reclosing Method of DC System Based on Main Frequency of Traveling Wave - Google Patents

Adaptive Reclosing Method of DC System Based on Main Frequency of Traveling Wave Download PDF

Info

Publication number
CN112564057B
CN112564057B CN202010597438.8A CN202010597438A CN112564057B CN 112564057 B CN112564057 B CN 112564057B CN 202010597438 A CN202010597438 A CN 202010597438A CN 112564057 B CN112564057 B CN 112564057B
Authority
CN
China
Prior art keywords
line
fault
traveling wave
frequency
time
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
Application number
CN202010597438.8A
Other languages
Chinese (zh)
Other versions
CN112564057A (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.)
Beijing Jiaotong University
State Grid Economic and Technological Research Institute
Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
Beijing Jiaotong University
State Grid Economic and Technological Research Institute
Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd
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 Beijing Jiaotong University, State Grid Economic and Technological Research Institute, Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd filed Critical Beijing Jiaotong University
Priority to CN202010597438.8A priority Critical patent/CN112564057B/en
Publication of CN112564057A publication Critical patent/CN112564057A/en
Application granted granted Critical
Publication of CN112564057B publication Critical patent/CN112564057B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02HEMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H7/00Emergency protective circuit arrangements specially adapted for specific types of electric machines or apparatus or for sectionalised protection of cable or line systems, and effecting automatic switching in the event of an undesired change from normal working conditions
    • H02H7/26Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured
    • H02H7/265Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured making use of travelling wave theory
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/081Locating faults in cables, transmission lines, or networks according to type of conductors
    • G01R31/085Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution lines, e.g. overhead
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/081Locating faults in cables, transmission lines, or networks according to type of conductors
    • G01R31/086Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/088Aspects of digital computing
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02HEMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H1/00Details of emergency protective circuit arrangements
    • H02H1/0007Details of emergency protective circuit arrangements concerning the detecting means
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02HEMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H3/00Emergency protective circuit arrangements for automatic disconnection directly responsive to an undesired change from normal electric working condition with or without subsequent reconnection ; integrated protection
    • H02H3/02Details
    • H02H3/06Details with automatic reconnection
    • H02H3/066Reconnection being a consequence of eliminating the fault which caused disconnection
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02HEMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H7/00Emergency protective circuit arrangements specially adapted for specific types of electric machines or apparatus or for sectionalised protection of cable or line systems, and effecting automatic switching in the event of an undesired change from normal working conditions
    • H02H7/26Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured
    • H02H7/268Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured for DC systems

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Emergency Protection Circuit Devices (AREA)

Abstract

本发明提供了一种基于行波主频率的直流系统自适应重合闸方法。该方法包括:在直流系统输电线路发生故障后,对线路的正、负极电压及正极电流进行采样,利用MUSIC算法计算出线路的线模电压行波的主频率,判断电压行波的主频率是否等于线路的固有特征频率的主频率,如果是,确定线路瞬时性故障,混合断路器重合闸;否则,利用STFT算法分别提取电压行波主频率和固有特征主频率对应的电压幅值,计算出特征信号幅值比值。根据特征信号幅值比值变化标定故障消失时刻,基于故障消失时刻确定混合断路器重合闸的具体时刻。本发明通过主动判断故障性质,避免了直流断路器盲目重合于故障,实现了直流断路器自适应重合闸。

Figure 202010597438

The invention provides an adaptive reclosing method for a DC system based on the dominant frequency of traveling waves. The method includes: sampling the positive and negative voltages and positive currents of the line after a fault occurs in the transmission line of the DC system, calculating the main frequency of the line-mode voltage traveling wave of the line by using the MUSIC algorithm, and judging whether the main frequency of the voltage traveling wave is not The main frequency equal to the natural characteristic frequency of the line, if so, determine the instantaneous fault of the line, and the hybrid circuit breaker recloses; otherwise, use the STFT algorithm to extract the voltage amplitude corresponding to the main frequency of the voltage traveling wave and the natural characteristic main frequency, and calculate Eigensignal amplitude ratio. The fault disappearance time is calibrated according to the change of the characteristic signal amplitude ratio, and the specific time of the hybrid circuit breaker reclosing is determined based on the fault disappearance time. By actively judging the nature of the fault, the invention avoids the blind reclosing of the DC circuit breaker on the fault, and realizes the self-adaptive reclosing of the DC circuit breaker.

Figure 202010597438

Description

Traveling wave dominant frequency-based direct current system self-adaptive reclosing method
Technical Field
The invention relates to the technical field of power system protection and control, in particular to a direct current system self-adaptive reclosing method based on traveling wave main frequency.
Background
The high-voltage direct-current transmission is used as a main line for trans-regional energy transmission and has the characteristics of long power transmission corridor, wide coverage area, large transmission capacity and the like. Considering economic factors, the direct current transmission mostly takes an overhead line as a channel for electric energy transmission. The overhead transmission line is exposed to the air, has severe working conditions, is an element with the highest fault probability in a direct current system, and mainly has transient faults. The research on the recovery scheme after the direct-current line fault is particularly important for guaranteeing the safety of equipment, improving the power supply reliability and maintaining the stability of the whole system.
When a line has a fault, a VSC-HVDC (voltage source converted based high voltage direct current) transmission system based on a voltage source converter has a converter outlet connected with a large capacitor in parallel and discharges to a fault point; even if the converter is locked, an alternating current system and a fault point can still form a loop through a freewheeling diode, and the fault cannot be cleared through a control system like a line communated converter based high voltage direct current (LCC-HVDC) system based on a grid commutation converter. Therefore, the flexible direct current system is additionally provided with the direct current circuit breakers at two ends of the line, the isolation of a fault line is realized, and the reliability of the operation of the system is improved.
The fault restarting strategy of the flexible direct current power transmission system based on the overhead line is as follows: converter lockout-circuit breaker trip- (fault nature identification I) -deionization- (fault nature identification II) -restart. The key of success or failure of system recovery is the identification of fault properties, and if the fault is an instantaneous fault, the circuit breaker is controlled to be switched on to establish system voltage; if the fault is a permanent fault, the converter is locked, and line maintenance is carried out. The concept of "adaptive reclosing" was proposed in the 80's of the 20 th century, with the main objective of avoiding reclosing under permanent faults. The relevant scholars further study deeply and clearly understand that the adaptive reclosing should also be able to capture the arc quenching time under transient faults to ensure the reclosing success in the voltage recovery phase. Currently, adaptive reclosure is generally applied to an alternating current power grid, but research on the technology in the field of direct current systems is just started.
At present, the fault recovery research methods of the flexible dc system in the prior art can be summarized into four categories: fault property identification schemes based on alternating current circuit breakers, based on converters, based on direct current circuit breakers, and based on characteristic electrical quantities.
The fault recovery scheme based on the alternating current circuit breaker mainly utilizes the matching of the alternating current circuit breaker and the direct current switch to realize fault isolation and recovery, the strategy can finish the isolation of a fault area without depending on communication, but the whole direct current system can be quitted from operation in a fault clearing stage, and the reclosing speed is slow.
The converter-based reclosing strategy may also be referred to as a restart strategy. The strategy mainly utilizes a converter topology with fault clearing capability, after the fault is isolated, the converter is operated in a specific mode, and the fault property is judged according to related electric quantities such as direct current voltage, direct current, bridge arm current and the like. However, the fault recovery scheme is limited by the structure of the converter, and when the fault recovery scheme is used for a multi-terminal direct-current power grid, the converter valve is locked to cause the whole direct-current system to exit, so that the power failure range is expanded.
Compared with the first two methods, the method has the advantages that fault lines can be selectively removed based on the fault clearing and system recovery strategies of the direct current circuit breakers, the power failure range after the fault is reduced to the maximum extent, the integral shutdown of the direct current system is avoided, single direct current line coincidence is achieved, however, the scheme needs to detect whether voltage is established to judge whether the fault disappears or not through multiple coincidence of the circuit breakers after the dissociation removing time is over, the pre-recognition of the fault property is lacked, and the coincidence strategy has certain blindness.
Disclosure of Invention
The embodiment of the invention provides a direct current system self-adaptive reclosing method based on traveling wave main frequency, which aims to overcome the problems in the prior art.
In order to achieve the purpose, the invention adopts the following technical scheme.
A direct current system self-adaptive reclosing method based on traveling wave main frequency comprises the following steps:
step 1, after a direct current system transmission line breaks down, continuously sampling positive and negative voltages and positive currents of the line, and calculating line mode voltage traveling waves of the line;
step 2, calculating the main frequency of the line mode voltage traveling wave by using the MUSIC algorithm, judging whether the main frequency of the line mode voltage traveling wave is equal to the main frequency of the inherent characteristic frequency of the line, if so, determining the instantaneous fault of the line, wherein the fault is at the midpoint of the line, the hybrid breaker of the line is reclosed, and the process is ended; otherwise, executing step 3;
and step 3: using short-time Fourier transform (STFT) algorithm to carry out linear mode voltage traveling wavePerforming time-frequency transformation, and extracting voltage amplitude U (f) corresponding to main frequency of the traveling wave of the line mode voltage and voltage amplitude U (f) corresponding to main frequency of the inherent characteristic frequencyL) According to said U (f) and U (f)L) Calculating the amplitude ratio of the characteristic signals:
and 4, calibrating the fault disappearance moment according to the change of the characteristic signal amplitude ratio, and determining the specific reclosing moment of the hybrid circuit breaker based on the fault disappearance moment.
Preferably, after the power transmission line of the dc system fails in step 1, continuously sampling the positive and negative voltages and the positive current of the line, and calculating the line mode voltage traveling wave of the line, includes:
at length L of DC system transmission line t1After the time is failed, the positive and negative voltages U of the line are applied according to a set time interval1p、U1nAnd positive electrode current I1pContinuously sampling, and calculating the line mode voltage traveling wave U of the line by using a phase-mode conversion formula1
Figure GDA0002903546060000031
Preferably, the line mode voltage traveling wave U is calculated by utilizing the MUSIC algorithm1Comprises:
after the mixed circuit breakers at two sides of the line are tripped, the positive current I is continuously detected according to a set time interval1pApplying a positive electrode current I1pThe time at which the current decays to 0.01 times the rated current is denoted as t3Intercepting line mode voltage traveling wave U in data window length1Data, the length of the data window is t-4L/v, v is the traveling wave propagation speed, and U in the length of the data window is calculated by utilizing the MUSIC algorithm1Power spectrum of data, maximum amplitude U in said data window1The minimum frequency corresponding to the wave crest is called as the line mode voltage traveling wave U1The main frequency f of (1).
Preferably, the step 3 of performing time-frequency transformation on the line-mode voltage traveling wave by using a short-time fourier transform (STFT) algorithm to extract the line-mode voltageVoltage amplitude U (f) corresponding to the main frequency of the traveling wave and voltage amplitude U (f) corresponding to the main frequency of the natural characteristic frequencyL) According to said U (f) and U (f)L) And calculating a characteristic signal amplitude ratio, including:
using STFT algorithm to measure U in data window length1Performing time-frequency transformation on the data to calculate U1Extracting the U of the line mode voltage traveling wave1The main frequency f of the natural characteristic frequency is extracted from the voltage amplitude U (f) corresponding to the main frequency fLCorresponding voltage amplitude U (f)L);
Defining the amplitude ratio of the characteristic signals as K:
Figure GDA0002903546060000041
preferably, the calibrating the fault disappearance time according to the change of the characteristic signal amplitude ratio in step 4, and determining the specific reclosing time of the hybrid circuit breaker based on the fault disappearance time includes:
at time t3Then, the anode current U is continuously monitored according to a set time interval1pAnd calculating the line mode voltage traveling wave U1The main frequency f and the characteristic signal amplitude ratio K;
judging the condition anode current U1pLess than 0.01 times of rated voltage, K<1, if yes, determining a permanent fault of the line, locking the hybrid circuit breaker, cutting off the fault line, and ending the process; otherwise, carrying out subsequent processing;
the judgment condition is t4Time of day t4>t3,K>1, and if so, determining that the line has a transient fault, the fault being at a point in the line, t4The moment is the disappearance moment of the fault, at t4And the hybrid breaker is reclosed at any moment.
According to the technical scheme provided by the embodiment of the invention, the direct-current circuit breaker is prevented from being blindly superposed on the fault by actively judging the fault property, the self-adaptive reclosing is realized, the reclosing success rate of the direct-current circuit breaker is improved, the closing reliability of the direct-current circuit breaker is improved, and the technical defect of blindly superposing of the direct-current circuit breaker in the prior art is overcome.
Additional aspects and advantages of the invention will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
Fig. 1 is a processing flow chart of a direct current system adaptive reclosing method based on traveling wave dominant frequency according to an embodiment of the present invention;
fig. 2 is a schematic diagram of a ± 200kV double-ended pseudo bipolar direct-current transmission system based on an MMC according to an embodiment of the present invention;
FIG. 3 is a schematic diagram showing the calculation results of the present invention applied to different faults, wherein FIG. 3(a) corresponds to the calculation result of the permanent fault of example 1, and FIG. 3(b) corresponds to the calculation result of the transient fault of example 2;
FIG. 4 is a diagram illustrating a calculation result of a dead zone detected at a midpoint of a line according to the present invention.
Detailed Description
Reference will now be made in detail to embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like reference numerals refer to the same or similar elements or elements having the same or similar function throughout. The embodiments described below with reference to the accompanying drawings are illustrative only for the purpose of explaining the present invention, and are not to be construed as limiting the present invention.
As used herein, the singular forms "a", "an", "the" and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms "comprises" and/or "comprising," when used in this specification, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof. It will be understood that when an element is referred to as being "connected" or "coupled" to another element, it can be directly connected or coupled to the other element or intervening elements may also be present. Further, "connected" or "coupled" as used herein may include wirelessly connected or coupled. As used herein, the term "and/or" includes any and all combinations of one or more of the associated listed items.
It will be understood by those skilled in the art that, unless otherwise defined, all terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. It will be further understood that terms, such as those defined in commonly used dictionaries, should be interpreted as having a meaning that is consistent with their meaning in the context of the prior art and will not be interpreted in an idealized or overly formal sense unless expressly so defined herein.
For the convenience of understanding the embodiments of the present invention, the following description will be further explained by taking several specific embodiments as examples in conjunction with the drawings, and the embodiments are not to be construed as limiting the embodiments of the present invention.
The invention provides a method for judging the bipolar short-circuit fault property of the flexible direct-current transmission line actively based on the discharge process of the line to the ground unit capacitor after the fault, aiming at the research requirement of the fault recovery of the flexible direct-current transmission system, and the self-adaptive reclosing is realized by calculating the amplitude ratio of the voltage traveling wave characteristic frequency signal and calibrating the fault disappearance moment.
The invention provides a direct current system self-adaptive reclosing method based on traveling wave main frequency. The main frequency change characteristic of the residual voltage traveling wave of the line after the insulation of the arrester is recovered is calculated, the characteristic frequency voltage signal amplitude ratio K is calculated, the fault disappearance moment is calibrated based on the change of the K value, and the circuit breaker is controlled to be reclosed.
For a DC line of length LWhen the circuit breaker breaks down, the circuit breakers on the two sides of the circuit trip, and the residual voltage traveling wave along the line is reflected between the circuit breakers on the two sides of the circuit and a fault point. Wherein, the breaker is equivalent to an open circuit, and the reflection coefficient of the traveling wave is + 1; and the fault point is short-circuited, and the traveling wave reflection coefficient is approximate to-1. Meter TL2L/v is the time length of the traveling wave going back and forth on the line; f. ofLv/2L is the main frequency of the inherent characteristic frequency corresponding to the line, v is the traveling wave propagation speed, and v is 3 × 108m/s。
The processing flow of the direct current system adaptive reclosing method based on the traveling wave main frequency provided by the embodiment of the invention is shown in fig. 1, and the method specifically comprises the following steps:
step 1: line t of length L1After the fault happens at any moment, the positive and negative voltages U of the line are applied1p、U1nAnd positive electrode current I1pContinuously sampling, and calculating line mode voltage traveling wave U by using a phase mode conversion formula shown in the specification1
Figure GDA0002903546060000071
Step 2: after the mixed circuit breakers positioned at two sides of the line are tripped, the tripping time of the mixed circuit breaker is recorded as t2. Continuously detecting the anode current I1pApplying a positive electrode current I1pThe time at which the current decays to 0.01 times the rated current is denoted as t3. Thereafter, the U in the time window is intercepted1And the length of the time window is t-4L/v.
Calculating U in time window by MUSIC (Multiple Signal Classification) algorithm1Power spectrum of data, U within data window1The frequency corresponding to the peak with the maximum amplitude is the minimum, and the minimum frequency is called as the line mode voltage traveling wave U1Is denoted as f, f is also called as fault feature main frequency f. After the line is in fault, the maximum value of the fault voltage traveling wave propagation period is 4L/v, and the maximum value can be used as a time window for FFT calculation.
And step 3: judging line mode voltage traveling wave U1Whether the main frequency of (A) is equal to the inherent characteristics of the lineDominant frequency f of the eigenfrequencyLIf so, determining the instantaneous fault of the line, enabling the fault to be at the middle point of the line, reclosing the hybrid circuit breaker, and ending the process; otherwise, step 4 is executed.
And 4, step 4: applying STFT (Short Time Fourier Transform) algorithm to U in Time window1Carrying out time-frequency transformation on the data, and extracting the traveling wave U of the outgoing line mode voltage1The voltage amplitude corresponding to the main frequency f is marked as U (f), and f is also extractedLCorresponding voltage amplitude, denoted as U (f)L) Defining a characteristic signal amplitude ratio K:
Figure GDA0002903546060000072
the STFT is one of common tools for analyzing signal time-frequency characteristics, has the characteristics of simple calculation and frequency scale linearity, and is an effective algorithm for extracting a specific frequency spectrum. The basic principle is that a section of long signal is divided into a limited number of short signal sections according to the length of a window, and fft (fast Fourier transform) is performed on each divided short signal section, so that the frequency components of each short signal section at different moments can be obtained. Thereafter. And arranging the obtained spectrogram of each short signal segment along a time axis, so as to obtain the time-varying relation of different frequency components of the original signal. The STFT can conveniently construct the change relation of each frequency signal contained in the signal along with time, and the fault characteristic frequency f and the line characteristic frequency f are already obtainedLUnder the condition (2), only the corresponding frequencies f and f are required to be marked in the time-frequency diagram obtained by STFTLThe signal amplitudes of (d) are U (f) and U (f)L)。
When the window of the STFT takes a larger value, a very high frequency resolution and a relatively low time resolution can be obtained, but considering that the self-adaptive reclosing strategy has no high requirement on time precision, the problem can be solved through proper time delay.
Taking into account the line mode voltage travelling wave U1Contains a plurality of frequency components, the frequency spectrum is very complex, and the algorithm is more pursuing the frequency resolution capability, so the Hanning window (Han) is selectedning) as a moving window at the time of STFT calculation, the window length was set to 2000.
And 5: time t3Then, the anode current U is continuously monitored according to a set time interval1pAnd calculating the line mode voltage traveling wave U1And judging the fault property according to the change rule of K by the main frequency f and the characteristic signal amplitude ratio K.
Judging the condition anode current U1pLess than 0.01 times of rated voltage, K<1, if yes, determining a permanent fault of the line, locking the hybrid circuit breaker, cutting off the fault line, and ending the process; otherwise, step 6 is executed.
Step 6: at t4Time (t)4>t3),K>1: then a line transient fault is determined, the fault being at a line midpoint, t4The moment is the disappearance moment of the fault, at t4And the hybrid breaker is reclosed at any moment.
Fig. 2 is a schematic diagram of a ± 200kV double-ended pseudo bipolar direct-current transmission system based on an MMC according to an embodiment of the present invention, where various parameters of the system are shown in table 1. Suppose the system is at t1When the time is 0s, a bipolar short-circuit fault occurs, and the direct current breaker is at t2Put into the arrester 3ms later. Wherein, the total length L of the line is 300kM, and the wave speed v of the travelling wave is 3 multiplied by 108m/s, main frequency f of line natural characteristic frequencyLv/2L 500 Hz; in addition, the arc resistance R of the fault branch is simulated by utilizing the Thomson arc model shown in the formula (1)arcWherein u isaIs the arc voltage; i.e. ifIs the arc current; i iss=10-18A is ionization current induced by external ionization (light, heat); c. C1,c2Is constant at 0.1MPa, c1≈2.86×106K·cm,c2≈7.47×104K·kV·cm。Rf20 Ω is the transition resistance of the fault branch; sampling frequency f 050 kHZ; STFT time frequency analysis selects Hanning window with window length of 2000.
Figure GDA0002903546060000091
TABLE 1 System simulation parameters
Figure GDA0002903546060000092
Example 1: distance to failure Lf60kM, permanent fault (no arc)
Step 1) failure at t1After 0ms, the positive and negative voltages U1p、U1nAnd current I1pContinuously sampling, and calculating line-mode voltage traveling wave U by using phase-mode conversion1
Step 2) lightning arrester at t23ms, positive electrode current I1PAt t3When the time is 6.2ms, the current is reduced to 0.01 times of rated current; with t3Intercepting data with the duration of t being 4L/v being 4ms for starting time, and calculating to obtain U by using MUSIC1The main frequency (frequency minimum) f is 1233 Hz.
Step 3) utilizing STFT to pair U1Performing time-frequency analysis, and extracting voltage amplitude signals U (f) corresponding to the fault characteristic main frequency f and the line characteristic main frequency fLCorresponding voltage amplitude signal U (f)L) And calculating the voltage signal amplitude ratio K.
Example 2: distance to failure Lf60kM, permanent fault (no arc)
Example 2: distance to failure Lf60kM, the nature of the fault is a transient fault (with arc)
Step 1) failure at t1After 0ms, the positive and negative voltages U1p、U1nAnd current I1pContinuously sampling, and calculating line-mode voltage traveling wave U by using phase-mode conversion1
Step 2) lightning arrester at t23ms, positive electrode current I1PAt t3When the time is 6.2ms, the current is reduced to 0.01 times of rated current; with t3Intercepting data with the duration of t being 4L/v being 4ms for starting time, and calculating to obtain U by using MUSIC1The main frequency (frequency minimum) f of (c) 1074 Hz.
Step 3) utilizing STFT to pair U1Performing time-frequencyAnalyzing and extracting voltage amplitude signals U (f) corresponding to the fault characteristic main frequency f and the line characteristic main frequency fLCorresponding voltage amplitude signal U (f)L) And calculating the voltage signal amplitude ratio K.
Fig. 3 shows the calculation results of the above-mentioned example, fig. 3(a) corresponds to the calculation result of the permanent fault of example 1, and fig. 3(b) corresponds to the calculation result of the transient fault of example 2.
As shown in FIG. 3(a), the upper portion of FIG. 3(a) is the DC line current I1pAnd fault branch current IfComparing the waveforms; the middle part of fig. 3(a) is a voltage U (f) corresponding to a fault main frequency of 1233Hz and a voltage U (f) of a line characteristic main frequency of 500HzL) Comparing the waveforms; the lower part of FIG. 3(a) is U (f) and U (f)L) Is measured.
The calculation result shows that K is at t3And then, the voltage is always less than 1, so that the detection criterion of permanent faults is met, the circuit breaker is locked, and a fault line is cut off.
As shown in FIG. 3(b), the upper portion of FIG. 3(b) is the DC line current I1pAnd fault branch current IfComparing the waveforms; the middle part of fig. 3(b) is the voltage U (f) corresponding to the fault main frequency 1074Hz and the voltage U (f) of the line characteristic main frequency 500HzL) Comparing the waveforms; the lower part of FIG. 3(b) is U (f) and U (f)L) Is measured. Wherein, t4Representing the actual arc current I resulting from the simulationfDisappearance time, t4' denotes the calculated arc current extinction time. Theoretically, t4=t4', but considering that STFT itself has a problem of insufficient time resolution, t4And t4There is some error between.
The calculation result shows that K is at t3Less than 1 after 6.2ms, but at t4' thereafter, K>1, the detection criterion of the transient fault is met, and the circuit breaker is reclosed.
Further analysis, when the fault occurs at the line midpoint (L)f150 kM/2), the characteristic frequency f is v/4L due to a faultfv/2L equals the natural characteristic frequency f of the lineLAt this time, U (f) is U (f)L) CalculatingThe resulting voltage amplitude ratio K is 1. Therefore, the midpoint of the line is the detection blind zone of the invention.
Example 3: distance to failure LfCompare permanent faults with transient faults at 150kM (line midpoint).
Step 1) failure at t1After 0ms, the positive and negative voltages U1p、U1nAnd current I1pContinuously sampling, and calculating line-mode voltage traveling wave U by using phase-mode conversion1
Step 2) lightning arrester at t23ms, positive electrode current I1PAt t3When the time is 6.7ms, the current is reduced to 0.01 times of rated current; with t3Intercepting data with the duration of t being 4L/v being 4ms for starting time, and calculating to obtain U by using MUSIC1The main frequencies (frequency minimum) of (f) 508Hz (permanent fault) and (f) 521Hz (transient fault).
Step 3) utilizing STFT to pair U1Performing time-frequency analysis, and extracting voltage amplitude signals U (f) corresponding to the fault characteristic main frequency f and the line characteristic main frequency fLCorresponding voltage amplitude signal U (f)L) And calculating the voltage signal amplitude ratio K.
FIG. 4 is a diagram illustrating a calculation result of a dead zone detected at a midpoint of a line according to the present invention. The calculation result shows that K is at t3Constantly equal to 1 after 6.7ms, it is difficult to distinguish the nature of the fault, but considering that the probability of the fault occurring exactly at the point in the line is low, it is treated as a transient fault.
In conclusion, the invention avoids the direct current breaker from being blindly superposed on the fault by actively judging the fault property, realizes self-adaptive reclosing, improves the success rate of the reclosing of the direct current breaker, improves the reliability of the closing of the direct current breaker and overcomes the technical defect of the blindly superposing of the direct current breaker in the prior art.
The method of the embodiment of the invention has simple and reliable principle, can accurately identify the bipolar short circuit fault property only by collecting the line voltage and the current magnitude after the fault, and has high identification degree and small detection blind area.
Those of ordinary skill in the art will understand that: the figures are merely schematic representations of one embodiment, and the blocks or flow diagrams in the figures are not necessarily required to practice the present invention.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for apparatus or system embodiments, since they are substantially similar to method embodiments, they are described in relative terms, as long as they are described in partial descriptions of method embodiments. The above-described embodiments of the apparatus and system are merely illustrative, and the units described as separate parts may or may not be physically separate, and the parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
The above description is only for the preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims (4)

1.一种基于行波主频率的直流系统自适应重合闸方法,其特征在于,包括:1. a direct current system adaptive reclosing method based on traveling wave dominant frequency, is characterized in that, comprises: 步骤1、在直流系统输电线路发生故障后,对线路的正、负极电压及正极电流进行持续采样,计算出线路的线模电压行波;Step 1. After the DC system transmission line fails, continuously sample the positive and negative voltages and positive current of the line, and calculate the line-mode voltage traveling wave of the line; 步骤2、利用MUSIC算法计算出线模电压行波的主频率,判断线模电压行波的主频率是否等于线路的固有特征频率的主频率,如果是,则确定线路瞬时性故障,故障在线路中点,线路的混合断路器重合闸,流程结束;否则,执行步骤3;Step 2. Use the MUSIC algorithm to calculate the main frequency of the line-mode voltage traveling wave, and judge whether the main frequency of the line-mode voltage traveling wave is equal to the main frequency of the natural characteristic frequency of the line. point, the hybrid circuit breaker of the line is reclosed, and the process ends; otherwise, go to step 3; 步骤3:利用短时傅里叶变换STFT算法对所述线模电压行波进行时频变换,提取出线模电压行波的主频率对应的电压幅值U(f)和固有特征频率的主频率对应的电压幅值U(fL),根据所述U(f)和U(fL)计算出特征信号幅值比值:Step 3: Use the short-time Fourier transform (STFT) algorithm to perform time-frequency transformation on the line-mode voltage traveling wave, and extract the voltage amplitude U(f) corresponding to the main frequency of the line-mode voltage traveling wave and the main frequency of the natural characteristic frequency. For the corresponding voltage amplitude U(f L ), the characteristic signal amplitude ratio is calculated according to the U(f) and U(f L ): 步骤4、根据所述特征信号幅值比值变化标定故障消失时刻,基于故障消失时刻确定混合断路器重合闸的具体时刻,包括:Step 4, calibrating the time of fault disappearance according to the change of the characteristic signal amplitude ratio, and determining the specific time of the hybrid circuit breaker reclosing based on the time of fault disappearance, including: 在时刻t3之后,按照设定的时间间隔持续监测正极电压U1p,并计算线模电压行波U1的主频率f,特征信号幅值比值K;其中t3为正极电流I1p衰减至0.01倍额定电流的时刻;After time t3 , the positive voltage U 1p is continuously monitored according to the set time interval, and the main frequency f of the line-mode voltage traveling wave U 1 and the characteristic signal amplitude ratio K are calculated; where t 3 is when the positive current I 1p decays to The moment of 0.01 times the rated current; 判断条件正极电压U1p小于0.01倍额定电压时,K<1是否成立,如果是,则确定线路永久性故障,所述混合断路器闭锁,切除故障线路,流程结束;否则,进行后续处理;When the positive voltage U 1p of the judgment condition is less than 0.01 times the rated voltage, whether K<1 is established, if so, it is determined that the line is permanently faulty, the hybrid circuit breaker is blocked, the faulty line is cut off, and the process ends; otherwise, follow-up processing is performed; 判断条件在t4时刻,t4>t3,K>1是否成立,如果是,则确定线路瞬时性故障,故障在线路中点,t4时刻为故障消失时刻,在t4时刻混合断路器重合闸。Judgment condition is at time t 4 , t 4 >t 3 , K>1 is established, if so, determine the line instantaneous fault, the fault is at the midpoint of the line, time t 4 is the time when the fault disappears, and the hybrid circuit breaker at time t 4 Reclosing. 2.根据权利要求1所述的方法,其特征在于,所述的步骤1的在直流系统输电线路发生故障后,对线路的正、负极电压及正极电流进行持续采样,计算出线路的线模电压行波,包括:2 . The method according to claim 1 , wherein in the step 1, after the DC system transmission line fails, the positive and negative voltages and the positive current of the line are continuously sampled, and the line mode of the line is calculated. 3 . Voltage traveling waves, including: 在长度为L的直流系统输电线路t1时刻发生故障后,按照设定的时间间隔对线路的正、负极电压U1p、U1n及正极电流I1p持续采样,利用相模变换公式计算出线路的线模电压行波U1After the fault occurs at the time t 1 of the DC transmission line of length L, the positive and negative voltages U 1p , U 1n and the positive current I 1p of the line are continuously sampled according to the set time interval, and the phase-mode transformation formula is used to calculate the line's line mode voltage traveling wave U 1 ;
Figure FDA0003340031430000021
Figure FDA0003340031430000021
3.根据权利要求2所述的方法,其特征在于,所述的利用MUSIC算法计算出线模电压行波U1的主频率,包括:3. method according to claim 2, is characterized in that, described utilizing MUSIC algorithm to calculate the main frequency of line-mode voltage traveling wave U 1 , comprising: 位于线路两侧的混合断路器跳闸后,按照设定的时间间隔持续检测正极电流I1p,将正极电流I1p衰减至0.01倍额定电流的时刻记为t3,截取数据窗长内的线模电压行波U1数据,所述数据窗的长度为t=4L/v,v为行波传播速度,利用MUSIC算法计算数据窗长内的U1数据的功率谱,所述数据窗内幅值最大的U1波峰对应的频率最小,称该最小频率为线模电压行波U1的主频率f。After the hybrid circuit breaker on both sides of the line trips, the positive current I 1p is continuously detected according to the set time interval, and the time when the positive current I 1p is attenuated to 0.01 times the rated current is recorded as t 3 , and the line mode within the data window length is intercepted. Voltage traveling wave U1 data, the length of the data window is t=4L/v, v is the traveling wave propagation speed, use the MUSIC algorithm to calculate the power spectrum of the U1 data within the data window length, and the amplitude value in the data window The frequency corresponding to the largest U 1 peak is the smallest, which is called the main frequency f of the line-mode voltage traveling wave U 1 . 4.根据权利要求3所述的方法,其特征在于,所述的步骤3的利用短时傅里叶变换STFT算法对所述线模电压行波进行时频变换,提取出线模电压行波的主频率对应的电压幅值U(f)和固有特征频率的主频率对应的电压幅值U(fL),根据所述U(f)和U(fL)计算出特征信号幅值比值,包括:4 . The method according to claim 3 , wherein in the step 3, the short-time Fourier transform (STFT) algorithm is used to perform time-frequency transformation on the line-mode voltage traveling wave, and the line-mode voltage traveling wave is extracted. 5 . The voltage amplitude U(f) corresponding to the main frequency and the voltage amplitude U(f L ) corresponding to the main frequency of the natural characteristic frequency, the characteristic signal amplitude ratio is calculated according to the U(f) and U(f L ), include: 利用STFT算法对数据窗长内的U1数据进行时频变换,计算U1的功率谱,提取出线模电压行波的U1的主频率f对应的电压幅值U(f),提取出固有特征频率的主频率fL对应的电压幅值U(fL);Use the STFT algorithm to perform time - frequency transformation on the U1 data within the data window length, calculate the power spectrum of U1, extract the voltage amplitude U( f ) corresponding to the main frequency f of U1 of the line-mode voltage traveling wave, and extract the inherent The voltage amplitude U(f L ) corresponding to the main frequency f L of the characteristic frequency; 定义特征信号幅值比值为K:Define the characteristic signal amplitude ratio as K:
Figure FDA0003340031430000022
Figure FDA0003340031430000022
CN202010597438.8A 2020-06-28 2020-06-28 Adaptive Reclosing Method of DC System Based on Main Frequency of Traveling Wave Active CN112564057B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010597438.8A CN112564057B (en) 2020-06-28 2020-06-28 Adaptive Reclosing Method of DC System Based on Main Frequency of Traveling Wave

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010597438.8A CN112564057B (en) 2020-06-28 2020-06-28 Adaptive Reclosing Method of DC System Based on Main Frequency of Traveling Wave

Publications (2)

Publication Number Publication Date
CN112564057A CN112564057A (en) 2021-03-26
CN112564057B true CN112564057B (en) 2022-02-15

Family

ID=75040869

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010597438.8A Active CN112564057B (en) 2020-06-28 2020-06-28 Adaptive Reclosing Method of DC System Based on Main Frequency of Traveling Wave

Country Status (1)

Country Link
CN (1) CN112564057B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113659541B (en) * 2021-07-23 2022-09-09 华中科技大学 Reclosing method and system for multi-terminal DC power grid based on waveform similarity matching
CN113625122A (en) * 2021-08-23 2021-11-09 全球能源互联网研究院有限公司 A method for identifying fault properties of pseudo-bipolar flexible DC systems

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS58127516A (en) * 1982-01-20 1983-07-29 三菱電機株式会社 Automatic reclosing circuit device
CN108551161A (en) * 2018-03-20 2018-09-18 西安交通大学 Utilize the transverter DC line fault guard method that actively injection and traveling wave couple
CN111313379A (en) * 2019-12-09 2020-06-19 国网江苏省电力有限公司镇江供电分公司 Disconnection protection method for comparing line voltage on two sides of line and spare power automatic switching

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2139087B1 (en) * 2008-06-25 2016-10-19 ABB Research LTD Circuit breaker with improved re-closing functionality
US8744637B1 (en) * 2013-09-12 2014-06-03 Deepak Maragal Methods and apparatus for detecting the clearance of fault in shunt reactor compensated transmission lines

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS58127516A (en) * 1982-01-20 1983-07-29 三菱電機株式会社 Automatic reclosing circuit device
CN108551161A (en) * 2018-03-20 2018-09-18 西安交通大学 Utilize the transverter DC line fault guard method that actively injection and traveling wave couple
CN111313379A (en) * 2019-12-09 2020-06-19 国网江苏省电力有限公司镇江供电分公司 Disconnection protection method for comparing line voltage on two sides of line and spare power automatic switching

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
HHT based adaptive re-closure method;Jiping Wang;《IEEE PES Innovative Smart Grid Technologies》;20120524;全文 *

Also Published As

Publication number Publication date
CN112564057A (en) 2021-03-26

Similar Documents

Publication Publication Date Title
CN107966633B (en) Method and system for rapidly judging single-phase earth fault line of power distribution network of power supply system
CN103344875B (en) Classification line selection method for single-phase earth fault of resonance earthing system
CN105119255B (en) Photovoltaic microgrid fault isolation method based on fault state
Zhang et al. A non-unit boundary protection of DC line for MMC-MTDC grids
CN106950459B (en) Distributed single-phase earth fault judgment method
CN110110856B (en) Power grid fault diagnosis method and system based on multi-source information deep reasoning
CN105548806A (en) Small current grounding fault line selection method based on transient traveling wave
Bo et al. Positional protection of transmission line using fault generated high frequency transient signals
CN103018627A (en) Adaptive fault type fault line detection method for non-effectively earthed system
CN111740391B (en) Transient high frequency power protection method and device for ring DC ship distribution network
CN112564057B (en) Adaptive Reclosing Method of DC System Based on Main Frequency of Traveling Wave
CN107543998B (en) DC side fault location system and method for multi-terminal flexible DC transmission system
CN109617026B (en) An ultra-high-speed protection method for ground fault of DC grid converter station
CN104237703B (en) Band shunt reactor circuit three-phase self adaptation reclosing nature of trouble recognition methods
CN108683160A (en) A kind of method for protecting transient state energy for direct current transmission line fault electric current
Dubey et al. Wavelet based energy function for symmetrical fault detection during power swing
CN111614065B (en) Line disconnection protection method based on line parameter identification
CN105425093B (en) A kind of short trouble early stage discrimination method based on voltage detecting
CN109061397A (en) Line fault area identification method
CN110568300B (en) A single-phase-to-ground fault identification method for distribution network based on multi-source information
CN102435896B (en) A fast identification method for intermittent ground faults in medium voltage power systems of ships
CN110045232A (en) A kind of system with non effectively earth ed neutral ground fault phase discrimination method
CN112436481B (en) Self-adaptive protection method for distribution line of neutral point non-effective grounding system
CN108493909A (en) The detection method of Distribution Network Failure based on Voltage Drop
CN108152680B (en) Method for detecting commutation failure of direct-current transmission

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