CN114707550A - Method for identifying direct current series arc fault at inverter end of photovoltaic system - Google Patents

Method for identifying direct current series arc fault at inverter end of photovoltaic system Download PDF

Info

Publication number
CN114707550A
CN114707550A CN202210352668.7A CN202210352668A CN114707550A CN 114707550 A CN114707550 A CN 114707550A CN 202210352668 A CN202210352668 A CN 202210352668A CN 114707550 A CN114707550 A CN 114707550A
Authority
CN
China
Prior art keywords
arc fault
signal
formula
calculating
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.)
Pending
Application number
CN202210352668.7A
Other languages
Chinese (zh)
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.)
Tianjin University of Science and Technology
Original Assignee
Tianjin 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 Tianjin University of Science and Technology filed Critical Tianjin University of Science and Technology
Priority to CN202210352668.7A priority Critical patent/CN114707550A/en
Publication of CN114707550A publication Critical patent/CN114707550A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/08Feature extraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/12Classification; Matching

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Short-Circuits, Discontinuities, Leakage, Or Incorrect Line Connections (AREA)

Abstract

The invention relates to a method for identifying a direct current series arc fault at an inverter end of a photovoltaic system, and belongs to the technical field of direct current series arc fault detection of the photovoltaic system. The technical scheme of the invention extracts the arc fault current spectrogram characteristics by using the dimensionless data distribution statistics, and the dimensionless data distribution statistics can better reflect the spectrogram characteristics, unify the numerical range and better adapt to the starting working condition of the inverter. The grey correlation analysis has a strong ability to process random information. The method finds out the relevance of various characteristics and indexes to be analyzed and researched in different types of information through certain data processing, and can realize accurate and effective identification of the direct-current series arc fault under the starting working condition of the photovoltaic system inverter.

Description

Method for identifying direct current series arc fault at inverter end of photovoltaic system
Technical Field
The invention relates to a method for identifying a direct current series arc fault at an inverter end of a photovoltaic system, and belongs to the technical field of direct current series arc fault detection of the photovoltaic system.
Background
The traditional mode of power generation, of the extensive type, has posed a great number of environmental problems, with an ever-increasing population and an ever-expanding energy demand, going into the 21 st century. While the photovoltaic power generation technology benefits mankind, many faults, especially arc faults, seriously affect the safe operation of a photovoltaic power station. The existing protection methods mainly include a time domain protection method, a frequency domain protection method, a time-frequency domain protection method, and the like. The time domain protection method detects faults through signal waveform characteristics, and the technical implementation method is simple but has poor accuracy. And further excavating more arc fault characteristics for fault detection by the time domain and time-frequency domain protection method. However, in the existing arc fault research, attention is often paid to the comparison of fault characteristics and steady-state characteristics, and the influence of characteristics generated by current and voltage changes on the detection and protection of the direct-current arc fault under the normal working condition of the photovoltaic system is ignored. In actual environment, because the inverter internally comprises a nonlinear power electronic device, when the inverter works, particularly in the starting process, the sudden change characteristics generated by nonlinear fluctuation in the current rising process are easily confused with the current sudden change when an electric arc occurs, the existing time domain, frequency domain and time-frequency domain protection method is easily interfered, and the problem of misoperation of a direct current series arc fault detection device under the starting working condition of the inverter is not completely solved in the prior art. Therefore, the direct-current series arc fault is correctly identified under the starting working condition of the photovoltaic system inverter, the misoperation and the operation rejection probability are reduced, and the method has important significance for safe operation of the photovoltaic system.
Disclosure of Invention
The invention aims to solve the technical problem of misoperation of an arc fault detection method caused by current disturbance during the starting of an inverter, and provides a direct current series arc fault identification method which can be used for the inverter end close to a photovoltaic system, so that accurate and effective direct current series arc fault identification can be carried out under the starting working condition of the inverter.
In order to solve the problems, the technical scheme of the invention is as follows: provides a method for identifying a direct current series arc fault at an inverter end of a photovoltaic system, which sequentially comprises the following steps,
step 1: extracting arc fault current spectrogram characteristics by using dimensionless data distribution statistics mainly comprising kurtosis, skewness, peak factors, impact factors, margin factors and wave factors;
step 2: the arc fault data is classified using a grey correlation.
Further, step 1 comprises the following substeps:
step 1.1: calculating kurtosis ku of a signal spectrum1(x) The calculation formula is
Figure BSA0000270227870000011
In the formula, E { x4And E { x }2Are the fourth and second moments of the signal x;
step 1.2: calculating skewness skew (x) of signal spectrum, wherein the calculation formula is
Figure BSA0000270227870000012
Where μ is the expectation of signal x and σ is the standard deviation;
step 1.3: calculating the peak value factor C (x) of the signal spectrum by the formula
Figure BSA0000270227870000013
Wherein RMS (x) is the signal effective value;
step 1.4: calculating the impact factor I (x) of the signal frequency spectrum by the formula
Figure BSA0000270227870000014
Wherein the rectified mean value ARV (x) is the mean value of the absolute values of the signals and is calculated by
Figure BSA0000270227870000015
Step 1.5: calculating a margin factor M (x) of a signal spectrum by the formula
Figure BSA0000270227870000021
Where R (x) is the square root amplitude of the signal, is the square of the mean of the arithmetic square root, and is calculated as
Figure BSA0000270227870000022
Step 1.6: calculating the form factor F (x) of the signal spectrum by the formula
Figure BSA0000270227870000023
Further, step 2 comprises the following substeps:
step 2.1: carrying out dimensionless processing on the signal X, and recording the dimensionless quantity as XiFor xiAny one element x in (1)i(k) The computational expression of the dimensionless processing is
Figure BSA0000270227870000024
In the formula xi(1) Is the first data of the signal;
step 2.2: calculating the grey correlation coefficient by the formula
Figure BSA0000270227870000025
Wherein rho is a resolution coefficient and takes a value between 0 and 1,
Figure BSA0000270227870000026
referred to as the two-stage minimum difference,
Figure BSA0000270227870000027
referred to as the two-level maximum difference;
step 2.3: calculating the grey correlation degree according to the grey correlation coefficient, wherein the calculation formula is
Figure BSA0000270227870000028
In the formula of omegakIs the weight coefficient of the k index;
step 2.4: respectively calculating gray correlation degree p of the low frequency band and the high frequency band of the arc fault current data according to dimensionless data distribution statistics in the low frequency band of 1-20kHz and the high frequency band of 40-60kHz1
Step 2.5: respectively calculating gray correlation degrees p of the low frequency band and the high frequency band of current data when the inverter is started according to dimensionless data distribution statistics in the low frequency band of 1-20kHz and the high frequency band of 40-60kHz2
Step 2.6: when p is1≥p2Time of day is identified as an arc fault, p1<p2An arc-free fault is identified.
The technical scheme of the invention utilizes the dimensionless data distribution statistics to extract the arc fault current spectrogram characteristics, the dimensionless data distribution statistics can better reflect certain characteristics of the spectrogram, and simultaneously, the dimensionless data distribution statistics unifies the numerical range and can better adapt to the starting working condition of the inverter, so the method is used for extracting the arc fault spectrogram characteristics. The grey correlation analysis has a strong ability to process random information. The method finds out the relevance of various characteristics and indexes to be analyzed and researched in different types of information through certain data processing, and determines the type of data to be tested according to main contradiction, main characteristics and main influence factors. Since the gray correlation analysis can reflect the overall trend of the characteristics without knowing typical rules, the requirement on the size of the sample volume is not high, and the final analysis result generally coincides with the qualitative analysis, so that the method has wide practicability.
Drawings
FIG. 1 is an RLC equivalent oscillation model of DC series arc fault
FIG. 2 is a schematic diagram of over-damped, critical damped, and under-damped oscillation waveforms
FIG. 3 is a schematic diagram of frequency spectra of over-damped, critically damped, and under-damped oscillation waveforms
FIG. 4 is a schematic circuit diagram of a DC series arc fault detection device
FIG. 5 shows a fault current and AD converted voltage signal according to a first embodiment of the present invention
FIG. 6 is a frequency spectrum of a fault current data window according to a first embodiment of the present invention
FIG. 7 is a diagram of a fault current and AD converted voltage signal according to a second embodiment of the present invention
FIG. 8 is a frequency spectrum of a fault current data window according to a second embodiment of the present invention
FIG. 9 shows a fault current and an AD converted voltage signal according to a third embodiment of the present invention
FIG. 10 is a frequency spectrum of a fault current data window in a third embodiment of the present invention
FIG. 11 is a diagram of a fault current and an AD converted voltage signal according to a fourth embodiment of the present invention
FIG. 12 is a frequency spectrum of a fault current data window in a fourth embodiment of the present invention
Detailed Description
The present invention will now be described in further detail with reference to the accompanying drawings, which are simplified schematic drawings that illustrate, by way of illustration only, the basic structure of the invention and, therefore, show only the components that are relevant to the invention.
The arc may be equivalent to an RLC oscillator circuit as shown in fig. 1. When arcing begins, a localized, isolated capacitively charged spot begins to appear between the two contact points. Under the condition of sufficiently high charge density, a strong electric field close to the edge of a light spot can cause electric breakdown, arc discharge with the length equal to that of current is generated, and capacitance C appears at two ends of a fault breakpointsLocal isolation point of (1), arc self-inductance LarcAnd arc resistance RarcAre connected in series.
Suppose self-inductance LarcAnd a resistance RarcApproximately time independent over the life of the arc, the arc can be equated with Larc、RarcAnd a capacitor CsThe second-order oscillation equivalent model is formed by the calculation formula of a differential equation
Figure BSA0000270227870000031
Wherein α ═ Rarc/2LarcIs the rate of the inductive damping,
Figure BSA0000270227870000032
is the natural frequency of the circuit, and alpha and omega are all non-negative numbers;
the differential equation corresponds to a characteristic equation of p2+2αp+ω2=0;
For an over-damped transition state, the circuit i (t) general solution in the circuit can be expressed as
Figure BSA0000270227870000033
In the formula
Figure BSA0000270227870000035
Is the initial ramp rate of the arc current;
from the relationship between α and ω, it can be seen that when α > ω, i.e., α > ω
Figure BSA0000270227870000034
p has two different real numbers, and the current is an over-damping signal; when α ═ ω, p has two real roots which are equal, in this caseThe current is a critical damping signal; when alpha is less than omega, p has no real number, and the current is an underdamped signal;
FIG. 2 shows that the frequency spectrum voltage is converted into dimensionless normalized amplitude value through normalization processing, the bandwidths of over-damping, critical damping and under-damping are sequentially reduced, the arc current belongs to an over-damping signal, and three damped oscillation waveforms of over-damping, critical damping and under-damping and corresponding frequency spectrums are schematically shown in FIG. 3, wherein fsRepresenting the sampling frequency, where the over-damped resistance is 4 times the critical damping and the under-damping is 1/4 times the critical damping, so the arc current signal has a wider frequency band in the frequency domain.
The principle of the dc series arc fault detection device is shown in fig. 4, and mainly includes a current transformer, an impedance matching circuit, a signal amplifying and filtering circuit, a digital-to-analog conversion AD unit, and a controller. The AD uses SM73201 of TI company, and has double-end differential mode input, the double-end differential mode input has stronger noise suppression capability compared with single-end input, the maximum sampling frequency is 250kHz, and the sampling frequency of the oscilloscope is 250kHz as the characteristic frequency of the arc fault current can reach more than 100 kHz.
To facilitate understanding of the public, the technical solution of the present invention is further described in detail by the following specific examples.
The first embodiment is as follows:
when an arc fault occurs, a fault current IarcWhen transient change occurs, as shown by the dotted line in fig. 5, the operating voltage of the photovoltaic system is 600V, the arcing current is 14A, and data used by the arc fault detection device needs to be converted by the current transformer and the AD and then output as a voltage signal UarcThe converted voltage signal is shown in fig. 5(b), wherein a region I represents a signal waveform when the photovoltaic system normally operates, a dashed region II represents a signal waveform when an arc fault occurs, and a region III represents a signal waveform after the system is completely powered off. It can be seen from the figure that the arc fault signal represented within the red dashed line of region II has a transient abrupt change characteristic which may be used to identify an arc fault.
In order to facilitate analysis of the change condition of the fault current, 6 data windows in the region I and the region II before and after the fault are selected for Fourier analysis, 3 data windows selected in the region I are sequentially named as A1, A2 and A3 according to the time sequence, and 3 data windows selected in the region II are respectively named as A4, A5 and A6. Comparison of spectral amplitudes of a1-a6 is shown in fig. 6, the spectrograms of three groups of steady operating currents a1, a2 and A3 are (a), (b) and (c) in the graph, respectively, and the spectrums of arc fault data windows a4, a5 and a6 are shown as (d), (e) and (f) in the graph. By comparing the stable operation and the fault voltage signal frequency spectrogram, the frequency spectrum shapes and amplitudes of the A4, the A5 and the A6 are different from those of the A1, the A2 and the A3, and the amplitudes of the frequency spectrum shapes and amplitudes are different from those of the low frequency band of 1-20kHz and the high frequency band of 40-60 kHz.
The second embodiment:
the arc fault test working condition of the electric arc fault test of the graph 7 is different from that of the former group, the working voltage of a photovoltaic system is 600V, the arc discharge current is 10A, and the more important difference is that the frequency spectrum of the electric arc fault test working condition only contains fault characteristics between 1 kHz and 20 kHz. As can be seen from fig. 7(a), regions I and II also represent the photovoltaic system normal state and the arc fault state, respectively, and region III represents the system power-off state. It can be seen that the arc fault characteristics are substantially the same under different operating conditions, but the smaller the arc fault current, the larger the current transient amplitude, as shown in fig. 7(b), since the arc energy is small and the arc burning becomes unstable.
Three sets of data windows before and after the fault are selected from the regions I and II for analysis, wherein stable working data windows before the fault are named A7, A8 and A9 respectively, arc fault data windows are named A10, A11 and A12 respectively, the transient change amplitude of the A10 waveform is maximum, and the arc burning is seen to be gradually stabilized due to the transient change from A10 to A12. The frequency spectrum corresponding to the data window is shown in fig. 8, wherein the frequency spectrum change amplitude of a7-a9 representing the steady-state operating current is small, the characteristic frequency band of a10-a12 representing the arc fault is concentrated in the low frequency band of 1-20kHz, and particularly, the frequency spectrum of a10 is relatively remarkably improved and presents a broadband characteristic.
Example three:
fig. 9 is a waveform of current change at the time of starting operation of the inverter, and as can be seen from fig. 9(a), as time advances, the current first rises rapidly and then tends to be smooth. The working condition can be divided into 9 stages, the current has 3 jump processes in the areas II, IV and VI, the corresponding voltage signal in the graph 9(b) also has obvious transient shock sudden change, the jump can be mistaken for arc fault characteristics, the correct identification of the arc detection device to the fault is influenced, and false operation is caused.
From each of the four II, IV, VI and VIII regions 1 data window was selected for analysis, named M1, M2, M3 and M4, respectively, where the data window M4 was selected from the stationary current region VIII for comparison analysis with the transient current, and the total window length was 1024. FIG. 10 is a graph of spectra corresponding to four portions of the data window M1-M4, where the M1 spectral amplitude has a boost in the 0-60kHz band and the spectral shape has narrow band characteristics. The spectrum amplitude of M2 and M3 is increased in the range of 0-20kHz, and the shape of the part is close to that of an arc fault, so that the part is easy to interfere with the identification of the arc fault and also causes misoperation. The M4 part is in a stable phase with small current changes and therefore low spectral amplitude.
Example four:
fig. 11 is a current waveform under the same inverter starting condition, and is different from fig. 8 in that the current is gradually increased in the region II after the current starts to be gradually increased in the region I, and there is no multi-jump process until the current in the region III is stabilized after 10 seconds, so that the analysis can be divided into three typical regions.
One data window, named M5, M6, and M7, each 1024 in length, is selected from regions I, II and III, respectively. FIG. 12 is the frequency spectrum corresponding to M5-M7, and it can be seen from the figure that the frequency spectrum of M5 has a larger amplitude in the 20kHz band, while the frequency spectrum of M6 and M7 has a slightly increased amplitude in the 20-40kHz band. From spectrogram shape analysis, M5-M7 did not exhibit significant broadband characteristics.

Claims (9)

1. A method for identifying a direct current series arc fault at an inverter end of a photovoltaic system is characterized by extracting arc fault current spectrogram characteristics by using dimensionless data distribution statistics mainly comprising kurtosis, skewness, peak factors, impact factors, margin factors and wave form factors.
2. Calculating kurtosis ku of a signal spectrum1(x) The calculation formula is
Figure FSA0000270227860000011
In the formula, E { x4And E { x }2The fourth moment and the second moment of the signal x are calculated, and the skewness skew (x) of the signal frequency spectrum is calculated according to the formula
Figure FSA0000270227860000012
Where μ is the expectation of the signal x and σ is the standard deviation.
3. Calculating the peak value factor C (x) of the signal spectrum by the formula
Figure FSA0000270227860000013
Wherein RMS (x) is the effective value of the signal, calculating the impact factor I (x) of the signal spectrum
Figure FSA0000270227860000014
Wherein the rectified mean value ARV (x) is the mean value of the absolute values of the signals and is calculated by
Figure FSA0000270227860000015
4. Calculating a margin factor M (x) of a signal spectrum by the formula
Figure FSA0000270227860000016
Where R (x) is the square root amplitude of the signal, is the square of the mean of the arithmetic square root, and is calculated as
Figure FSA0000270227860000017
Calculating the form factor F (x) of the signal spectrum by the formula
Figure FSA0000270227860000018
5. The photovoltaic system inverter-side direct current series arc fault identification method of claim 1, wherein arc fault data is classified using a grey correlation.
6. Carrying out dimensionless processing on the signal X, and recording the dimensionless quantity as XiFor xiAny one element x in (1)i(k) The computational expression of the dimensionless processing is
Figure FSA0000270227860000019
In the formula xi(1) Is the first data of the signal.
7. Calculating the grey correlation coefficient by the formula
Figure FSA00002702278600000110
Wherein rho is a resolution coefficient and takes a value between 0 and 1,
Figure FSA00002702278600000111
referred to as the two-stage minimum difference,
Figure FSA00002702278600000112
called two-stage maximum difference, the grey correlation degree is calculated according to the grey correlation coefficient, and the calculation formula is
Figure FSA00002702278600000113
In the formula of omegakIs the weighting factor of the k index.
8. Respectively calculating gray correlation degree p of the low frequency band and the high frequency band of the arc fault current data according to dimensionless data distribution statistics in the low frequency band of 1-20kHz and the high frequency band of 40-60kHz1According to the invaluity in the low frequency band of 1-20kHz and the high frequency band of 40-60kHzCalculating gray correlation degree p between the current data and the current data at the start of the inverter according to the dimension data distribution statistics2
9. When p is1≥p2Time of day is identified as an arc fault, p1<p2An arc-free fault is identified.
CN202210352668.7A 2022-04-06 2022-04-06 Method for identifying direct current series arc fault at inverter end of photovoltaic system Pending CN114707550A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210352668.7A CN114707550A (en) 2022-04-06 2022-04-06 Method for identifying direct current series arc fault at inverter end of photovoltaic system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210352668.7A CN114707550A (en) 2022-04-06 2022-04-06 Method for identifying direct current series arc fault at inverter end of photovoltaic system

Publications (1)

Publication Number Publication Date
CN114707550A true CN114707550A (en) 2022-07-05

Family

ID=82172449

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210352668.7A Pending CN114707550A (en) 2022-04-06 2022-04-06 Method for identifying direct current series arc fault at inverter end of photovoltaic system

Country Status (1)

Country Link
CN (1) CN114707550A (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110501631A (en) * 2019-08-19 2019-11-26 重庆大学 A kind of online intermittent fault detection and diagnostic method
CN113552447A (en) * 2021-07-27 2021-10-26 上海电机学院 Series arc fault detection method based on random forest

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110501631A (en) * 2019-08-19 2019-11-26 重庆大学 A kind of online intermittent fault detection and diagnostic method
CN113552447A (en) * 2021-07-27 2021-10-26 上海电机学院 Series arc fault detection method based on random forest

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈进: "低压交流故障电弧的特征提取与识别", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》, 15 March 2022 (2022-03-15), pages 3 - 4 *

Similar Documents

Publication Publication Date Title
US8680873B2 (en) Method and apparatus for fault identification in a power tansmission line
Wang et al. Classification of power quality events using optimal time-frequency representations-part 1: Theory
Song et al. A negative selection algorithm-based identification framework for distribution network faults with high resistance
CN112730964A (en) Lightning overvoltage identification method based on overvoltage waveform characteristics
CN113191192A (en) Breaker fault detection method based on wavelet analysis and fuzzy neural network algorithm
CN109001602A (en) Shelf depreciation Severity method based on extreme learning machine algorithm
Wang et al. Stator single-line-to-ground fault protection for bus-connected powerformers based on S-transform and bagging ensemble learning
CN114707550A (en) Method for identifying direct current series arc fault at inverter end of photovoltaic system
Bonde et al. Review on detection and classification of underlying causes of power quality disturbances using signal processing and soft computing technique
CN112255495B (en) Micro-grid high-resistance fault detection method
CN111896890B (en) Micro-grid line fault diagnosis method and system based on Hilbert-Huang transform
Yao et al. Pattern recognition for partial discharge using multi-feature combination adaptive boost classification model
CN114113894A (en) Repetitive fault identification method based on natural frequency characteristics
Gao et al. Quantitative research on accumulative effect of transformer winding deformation and its influence degree based on time-frequency analysis of vibration signal
Gao et al. Internal overvoltage identification of distribution network via time-frequency atomic decomposition
CN106771850A (en) Distribution network fault line selection method based on cascade bistable system
Do et al. Wavelet packet-based passive islanding detection method for grid connected photovoltaic inverters
Wang et al. Faulty line detection for cross‐line same‐phase successive ground faults in distribution network based on transient characteristics
Wu et al. A Novel Method for Extracting Features of Partial Discharge Pulse Current Signal
CN115047240B (en) Transformer excitation surge current discriminating method by utilizing wavelet detail component change characteristics
CN113189446B (en) Low-current ground fault line selection method
CN113311219B (en) Power distribution network temporary overvoltage identification method
Shuai et al. Pattern recognition of partial discharge ultrahigh frequency signal based on similar matrix BSS and deep learning CNN
CN116500383B (en) Method and system for identifying high-resistance faults and switching disturbance based on Lissajous curves
CN114062833B (en) Single-phase earth fault insulation degradation identification method based on fuzzy inference system

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