Nothing Special   »   [go: up one dir, main page]

CN112345895B - Series arc fault detection method based on discriminant analysis strategy - Google Patents

Series arc fault detection method based on discriminant analysis strategy Download PDF

Info

Publication number
CN112345895B
CN112345895B CN202011168112.XA CN202011168112A CN112345895B CN 112345895 B CN112345895 B CN 112345895B CN 202011168112 A CN202011168112 A CN 202011168112A CN 112345895 B CN112345895 B CN 112345895B
Authority
CN
China
Prior art keywords
matrix
current data
vector
factor
arc fault
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
CN202011168112.XA
Other languages
Chinese (zh)
Other versions
CN112345895A (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.)
Ningbo Lxing Polytron Technologies Inc
Original Assignee
Ningbo Lxing Polytron Technologies Inc
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 Ningbo Lxing Polytron Technologies Inc filed Critical Ningbo Lxing Polytron Technologies Inc
Priority to CN202011168112.XA priority Critical patent/CN112345895B/en
Publication of CN112345895A publication Critical patent/CN112345895A/en
Application granted granted Critical
Publication of CN112345895B publication Critical patent/CN112345895B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/12Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing
    • G01R31/1227Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing of components, parts or materials
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R19/00Arrangements for measuring currents or voltages or for indicating presence or sign thereof
    • G01R19/0092Arrangements for measuring currents or voltages or for indicating presence or sign thereof measuring current only

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Relating To Insulation (AREA)

Abstract

The invention discloses a series arc fault detection method based on a discriminant analysis strategy, which is characterized in that an arc generating circuit is detected through a current sensor and a data acquisition card, current data in the arc generating circuit is acquired through a computer, characteristic transformation is implemented, a current characteristic matrix is constructed and standardized processing is implemented to obtain a new characteristic matrix, a Fisher discriminant analysis algorithm is utilized to obtain a classification transformation matrix, and the classification characteristic matrix and a residual error matrix are calculated. The collected current data are constructed into a current data vector, the mean value, the standard deviation kurtosis, the skewness, the root mean square, the peak factor, the shape factor, the pulse factor, the edge factor and the maximum logarithm are obtained through feature transformation, a feature vector is constructed, and the feature vector is subjected to standardization processing to obtain a detection index. The method has the advantages that only current data under different loads when the circuit normally operates are needed, current data under an arc fault state do not need to be simulated or actually sampled, and negative effects on the whole circuit are small.

Description

Series arc fault detection method based on discriminant analysis strategy
Technical Field
The invention relates to an arc fault detection method, in particular to a series arc fault detection method based on a discriminant analysis strategy.
Background
With the trend of increasing complexity and diversification of electric equipment, the problem of electric safety is a non-negligible problem. According to the statistics of the national fire department, the fire caused by the problems of short circuit, overload, electric leakage and the like of electrical equipment is the main cause of fire. The current electric equipment is generally provided with circuit units such as a circuit breaker, a fuse, electric leakage protection and the like, so that the problems of short circuit, overload and electric leakage are solved respectively. In addition to these three typical and common types of electrical equipment faults, an arc fault refers to an electrical fault in which an unintended arc occurs in a live line. According to the connection relationship between the electric Arc and the circuit when the electric Arc Fault occurs, the electric Arc Fault can be divided into Series Arc Fault (SAF), parallel Arc Fault (PAF), ground Arc Fault (GAF) and the like. The parallel arc fault and the grounding arc fault are equivalent to a short-circuit fault due to large current during fault, and a circuit breaker or a fuse can be directly operated to play a protection role. When a series arc fault occurs, the load current is smaller than the normal working current due to the self impedance, and the circuit breaker and the fuse cannot take corresponding actions. Therefore, series arc fault detection is a weak link in current power supply and distribution system fault detection.
Generally, the task of detecting the series arc fault mainly lies in successfully detecting whether the fault occurs or not, and the diagnosis of the type of the arc fault is performed after corresponding protective measures are taken. There are basically two main types of strategies for implementing series arc fault detection: firstly, the arc fault detection is implemented by installing a physical sensing device to acquire physical characteristic signals in real time; second, data-driven arc fault detection is implemented using current or voltage signals in the circuit. The former requires installation of hardware equipment, increases the investment cost of the equipment, and is limited by the installation location space. In contrast, the latter belongs to software detection, and online real-time monitoring of series arc faults can be realized only by using a proper data mining algorithm. Therefore, implementing data-driven series arc fault detection is currently the main direction of research.
It is noted that under normal operating conditions, series arcs can exhibit different characteristics of change in current in the circuit due to different loads and complex operating conditions. The traditional implementation of arc fault detection by using data is mainly based on the idea of multi-class classification, such as the use of decision trees, neural networks, support vector machines, bayesian classifiers, and the like. In addition, in consideration of the current change characteristics, the current signal is subjected to characteristic transformation such as time domain characteristic analysis and frequency domain characteristic analysis, and then classification is performed by using the transformed characteristics, thereby diagnosing the series arc fault. From this point of view, the feature selection plays a key role, that is, it needs to select a proper feature to improve the precision effect of classification, which derives a new problem. The data-driven serial arc fault detection is implemented, the characteristics are analyzed and converted through a time domain or a frequency domain, the multi-working-condition characteristics of a serial circuit need to be considered, the characteristics are extracted through the converted characteristics, the characteristic components are fully utilized, and the arc fault is timely and effectively detected.
Disclosure of Invention
The invention aims to solve the main technical problems that: the purpose of series arc fault detection can be achieved only according to current data under different load conditions. Specifically, the method firstly carries out time domain and frequency domain feature analysis on current data under different load conditions under normal working conditions, then converts the time-frequency domain features into feature signals and residual error signals for classification by utilizing discriminant analysis, and finally carries out fault detection on the classification feature signals and the residual error feature signals respectively.
The technical scheme adopted by the method for solving the problems is as follows: a series arc fault detection method based on a discriminant analysis strategy comprises the following steps:
step 1: the method comprises the following steps of detecting an arc generating circuit through a current sensor and a data acquisition card, acquiring current data in the arc generating circuit through a computer, and dividing the acquired current data into three types: resistive load current data
Figure GDA0003939444150000021
Inductive load current data
Figure GDA0003939444150000022
And non-linear load current data
Figure GDA0003939444150000023
Wherein,
Figure GDA0003939444150000024
represents 1 XN 1 A vector of real numbers of the dimensions,
Figure GDA0003939444150000025
represents 1 XN 2 The vector of real numbers of the dimension(s),
Figure GDA0003939444150000026
represents 1 XN 3 Real number vector of dimension, N 1 ,N 2 ,N 3 Respectively representing the number of resistive load current data, the number of inductive load current data and the number of nonlinear load current data;
step 2: performing characteristic transformation on the three types of current data obtained in the step 1 to obtain a resistive load current characteristic matrix
Figure GDA0003939444150000027
Inductive load current signature matrix
Figure GDA0003939444150000028
And a non-linear load current signature matrix
Figure GDA0003939444150000029
Wherein,
Figure GDA00039394441500000210
represents n 1 A real number matrix of x 10 dimensions,
Figure GDA00039394441500000211
represents n 3 X 10 dimensional real number matrix, n 1 ,n 2 ,n 3 Respectively represents Z 1 ,Z 2 ,Z 3 The number of the row vectors of (a),
and step 3: constructing a current feature matrix Z = [ Z ] 1 T ,Z 2 T ,Z 3 T ] T And standardizing each column vector in Z to obtain a new feature matrix
Figure GDA00039394441500000212
Wherein the upper symbol T represents the transpose of a matrix or vector, n = n 1 +n 2 +n 3 ,R n×10 A real number matrix representing n × 10 dimensions;
and 4, step 4: discrimination by FisherThe analysis algorithm obtains a classification transformation matrix W belonging to R 10×2 And calculating a classification feature matrix
Figure GDA00039394441500000213
Sum residual matrix
Figure GDA00039394441500000214
Wherein the load matrix
Figure GDA00039394441500000215
And 5: and (3) carrying out singular value decomposition on the residual matrix E: e = UAV T And calculating a transformation matrix Q = VA -1 (ii) a Wherein U and V represent unitary matrix of singular value decomposition, the diagonal elements in diagonal matrix A are composed of singular values;
step 6: according to the formula psi = diag { UU T Calculating a detection index vector psi and determining a maximum value of psi as psi 0
And 7: the current data collected at the sampling moment are formed into a current data vector x t Where t represents the latest sampling instant, x is obtained by feature transformation t Mean value m of t Standard deviation of delta t Degree of kurtosis alpha t Deviation of degree beta t Root mean square upsilon t Crest factor gamma t Shape factor η t Pulse factor e t Edge factor f t Maximum logarithm of g t (ii) a And constructing a feature vector z t =[m t ,δ t ,α t ,β t ,υ t ,γ t ,η t ,e t ,f t ,g t ]∈R 1×10 For feature vector z t Carrying out standardization processing to obtain new characteristic vector
Figure GDA00039394441500000216
And 8: using formulas
Figure GDA00039394441500000217
Computing residual vector E t And detecting the index psi t =E t E t T
And step 9: and judging whether the step 8 meets the condition: psi t ≤ψ 0 (ii) a If yes, the series arc fault does not occur, and the steps 7 to 9 are repeated to continue to detect the arc fault of the arc generating circuit; if not, detecting the series arc fault.
Preferably, the specific implementation method for performing feature transformation on the three types of current data obtained in step 1 in step 2 is as follows:
setting the window length as L and the moving step length as D,
step 2.1: initializing j =1;
step 2.2: mixing X 1 The current interval vector q is composed of the (j-1) × L +1 to the (j × L) th elements j ∈R 1×L (ii) a Wherein R is 1×L A real number vector representing 1 × L dimensions;
step 2.3: respectively calculating mean values m according to formulas (1) to (R) j Standard deviation of delta j Degree of kurtosis alpha j Deviation of degree beta j Root mean square upsilon j Crest factor gamma j Shape factor η j Pulse factor e j Edge factor f j Maximum logarithm of g j
Figure GDA0003939444150000031
Figure GDA0003939444150000032
Figure GDA0003939444150000033
Figure GDA0003939444150000034
Figure GDA0003939444150000035
Figure GDA0003939444150000036
η j =υ j /m j
f j =q j (max)/m j
Figure GDA0003939444150000037
g j =log 10 (q j (max)) ⑩
Step 2.4: constructing a feature vector z j =[m j ,δ j ,α j ,β j ,υ j ,γ j ,η j ,e j ,f j ,g j ]∈R 1×10 And then judging whether the conditions are met: (j + 1). Times.L > N 1 (ii) a If not, after j = j +1 is set, returning to the step (2.2); if yes, all the obtained eigenvectors form a resistive load current characteristic matrix
Figure GDA0003939444150000038
Preferably, X in step 2.2 is 1 Is replaced by X 2 、X 3 And respectively obtaining the inductive load current characteristic matrixes according to the implementation processes of the step 2.1 to the step 2.4
Figure GDA0003939444150000039
And a non-linear load current signature matrix
Figure GDA00039394441500000310
Preferably, the latest sampling time is set as t, and the current data I acquired at the latest sampling time in the step 7 is set as t And the first L-1 thereofCurrent data I at the sampling instants t-1 ,I t-2 ,…,I t-L+1 Built into a current data vector x t =[I t ,I t-1 ,…,I t-L+1 ]∈R 1×L (ii) a Coupling X in step 2.2 1 Is replaced by x t Obtaining x t Mean value m of t Standard deviation of delta t Degree of kurtosis alpha t Degree of deviation beta t Root mean square upsilon t Crest factor gamma t Shape factor η t Pulse factor e t Edge factor f t Maximum logarithm of g t
By carrying out the steps described above, the advantages of the method of the invention are introduced as follows:
firstly, the method processes the multi-load characteristic of current data by discriminant analysis, thereby detecting whether a fault occurs by using residual errors; secondly, the method can judge whether the circuit generates the electric arc only by utilizing the current data under different loads when the circuit normally operates without simulating or actually sampling the current data under the electric arc fault state, thereby trying to reduce the influence on the circuit caused by the traditional data acquisition through equipment and realizing the purpose of minimizing the negative influence on the whole circuit; finally, the method provided by the invention uses the time domain characteristics of the current signal to carry out data mining, and is simpler and more practical.
Drawings
FIG. 1 is a schematic flow chart of the method of the present invention.
Detailed Description
The present invention will be described in detail below with reference to the accompanying drawings and specific embodiments.
As shown in fig. 1, the invention discloses a series arc fault detection method based on a discriminant analysis strategy. A specific embodiment of the method of the present invention is further set forth below in conjunction with fig. 1.
And constructing an arc generating circuit, namely connecting resistive, inductive and nonlinear loads with an arc generator and a current sensor respectively to form the arc generating circuit.
Step 1: by current sensors and data acquisitionThe card detects the arc generating circuit, and acquires current data in the arc generating circuit through the computer, and divides the acquired current data into three types: resistive load current data
Figure GDA0003939444150000041
Inductive load current data
Figure GDA0003939444150000042
And non-linear load current data
Figure GDA0003939444150000043
Wherein,
Figure GDA0003939444150000044
represents 1 XN 1 The vector of real numbers of the dimension(s),
Figure GDA0003939444150000045
represents 1 XN 2 Real number vector of dimension, N 1 ,N 2 ,N 3 Respectively representing the data number of resistive load current, inductive load current and nonlinear load current;
and 2, step: setting the window length as L and the moving step length as D, and performing characteristic transformation on the three types of current data to obtain a resistive load current characteristic matrix
Figure GDA0003939444150000046
Inductive load current signature matrix
Figure GDA0003939444150000047
And a non-linear load current signature matrix
Figure GDA0003939444150000048
Wherein,
Figure GDA0003939444150000049
represents n 1 A real number matrix of x 10 dimensions,
Figure GDA00039394441500000410
represents n 3 X 10 dimensional real number matrix, n 1 ,n 2 ,n 3 Respectively represents Z 1 ,Z 2 ,Z 3 The specific implementation process of the number of row vectors of (2) is as follows,
step 2.1: initializing j =1;
step 2.2: mixing X 1 The current interval vector q is composed of the (j-1) × L +1 to the (j × L) th elements j ∈R 1×L (ii) a Wherein R is 1×L A real number vector representing 1 × L dimensions;
step 2.3: calculating mean value m according to above formula (1) to formula (R) j Standard deviation of delta j Degree of kurtosis alpha j Deviation of degree beta j Root mean square upsilon j Crest factor gamma j Shape factor η j Pulse factor e j Edge factor f j Maximum logarithm of g j
Step 2.4: constructing a feature vector z j =[m j ,δ j ,α j ,β j ,υ j ,γ j ,η j ,e j ,f j ,g j ]∈R 1×10 And then judging whether the conditions are met: (j + 1). Times.L > N 1 (ii) a If not, after j = j +1 is set, returning to the step (2.2); if yes, all the obtained eigenvectors form a resistive load current characteristic matrix
Figure GDA0003939444150000051
Mixing X in the step 2.2 1 Is replaced by X 2 、X 3 And according to the implementation process from step 2.1 to step 2.4, the inductive load current characteristic matrixes can be respectively obtained
Figure GDA0003939444150000052
And a non-linear load current signature matrix
Figure GDA0003939444150000053
By standardizationProcessing to obtain a new feature matrix
Figure GDA0003939444150000054
Where the upper index T represents the transpose of the matrix or vector, n = n 1 +n 2 +n 3 ,R n×10 Representing a matrix of real numbers of dimension n x 10.
Step 4, obtaining a classification transformation matrix W belonging to R by utilizing a Fisher discriminant analysis algorithm 10×2 And calculating a classification feature matrix
Figure GDA0003939444150000055
Sum residual matrix
Figure GDA0003939444150000056
Wherein the load matrix
Figure GDA0003939444150000057
The specific implementation process of the Fisher discriminant analysis algorithm is as follows:
step 4.1: new feature matrix
Figure GDA0003939444150000058
Front n in (1) 1 The row vectors of the rows form a matrix
Figure GDA0003939444150000059
N of (1) 1 +1 line to nth line 1 +n 2 The row vectors of the rows form a matrix
Figure GDA00039394441500000510
After n in (1) 3 The row vectors of the rows form a matrix
Figure GDA00039394441500000511
Step 4.2: separately computing matrices
Figure GDA00039394441500000512
And
Figure GDA00039394441500000513
line mean vector mu of 1 ,μ 2 ,μ 3 The calculation method of the line mean vector is as follows: adding the row vectors in the matrix, and dividing the sum by the number of the row vectors;
step 4.3: the covariance matrix C between classes is calculated according to the formula 0 And within-class covariance matrix C 1
Figure GDA00039394441500000514
Figure GDA00039394441500000515
In the above formula, the first and second carbon atoms are,
Figure GDA00039394441500000516
representation matrix
Figure GDA00039394441500000517
The row vector of the j-th row in (1), i ∈ {1,2,3}, j ∈ {1,2, …, n i };
Step 4.4: calculation matrix C = C 1 -1 C 0 Maximum two eigenvalues λ 1 And λ 2 And its corresponding feature vector w 1 And w 2 Reconstructing a classification transformation matrix W = [ W = [) 1 ,w 2 ];
And 5: and (3) performing singular value decomposition on the residual matrix E: e = UAV T And calculating a transformation matrix Q = VA -1 (ii) a Wherein U and V represent unitary matrix of singular value decomposition, the diagonal elements in diagonal matrix A are composed of singular values;
step 6: according to the formula psi = diag { UU T Calculating a detection index vector psi and determining a maximum value of psi as psi 0 (ii) a Wherein diag { } denotes an operation of converting an element on a diagonal of the matrix in braces into a column vector;
and 7: current data I acquired at the latest sampling moment t And first L-1 sampling timeCurrent data of the moment I t-1 ,I t-2, …,I t-L+1 Built into a current data vector x t =[I t ,I t-1 ,…,I t-L+1 ]∈R 1×L (ii) a Where t represents the latest sampling instant. Mixing X in the step 2.2 1 Is replaced by x t Respectively calculate to obtain x t Mean value m of t Standard deviation of delta t Degree of kurtosis alpha t Deviation of degree beta t Root mean square upsilon t Crest factor gamma t Shape factor eta t Pulse factor e t Edge factor f t Maximum logarithm of g t (ii) a And constructing a feature vector z t =[m t ,δ t ,α t ,β t ,υ t ,γ t ,η t ,e t ,f t ,g t ]∈R 1×10 And carrying out the same standardization processing as the step 3 to obtain a new feature vector
Figure GDA0003939444150000061
And 8: using formulas
Figure GDA0003939444150000062
Computing residual vector E t And detecting the index psi t =E t E t T
And step 9: psi is judged in step 8 t Whether the condition is satisfied: psi t ≤ψ 0 (ii) a If yes, the series arc fault does not occur, and the steps 7 to 9 are repeated to continue to detect the arc fault of the arc generating circuit; if not, the series arc fault is detected, so that the detection of the arc generating circuit is realized, measures are convenient to take to break the circuit, and the generation of serious accidents is avoided.
The above-mentioned embodiments are only preferred embodiments of the present invention, and the description thereof is more specific and detailed, but not to be construed as limiting the scope of the present invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the inventive concept, which falls within the scope of the present invention.

Claims (4)

1. A series arc fault detection method based on a discriminant analysis strategy is characterized by comprising the following steps:
step 1: the method comprises the following steps of detecting an arc generating circuit through a current sensor and a data acquisition card, acquiring current data in the arc generating circuit through a computer, and dividing the acquired current data into three types: resistive load current data
Figure FDA0003939444140000011
Inductive load current data
Figure FDA0003939444140000012
And non-linear load current data
Figure FDA0003939444140000013
Wherein,
Figure FDA0003939444140000014
represents 1 XN 1 A vector of real numbers of the dimensions,
Figure FDA0003939444140000015
represents 1 XN 2 The vector of real numbers of the dimension(s),
Figure FDA0003939444140000016
represents 1 XN 3 Real number vector of dimension, N 1 ,N 2 ,N 3 Respectively representing the number of resistive load current data, the number of inductive load current data and the number of nonlinear load current data;
step 2: performing characteristic transformation on the three types of current data obtained in the step 1 to obtain a resistive load current characteristic matrix
Figure FDA00039394441400000113
Inductive load current signature matrix
Figure FDA00039394441400000114
And a non-linear load current signature matrix
Figure FDA00039394441400000115
Wherein,
Figure FDA00039394441400000116
represents n 1 A real number matrix of x 10 dimensions,
Figure FDA00039394441400000117
represents n 3 X 10 dimensional real number matrix, n 1 ,n 2 ,n 3 Respectively represents Z 1 ,Z 2 ,Z 3 The number of the row vectors;
and step 3: constructing a current characteristic matrix Z = [ Z ] 1 T ,Z 2 T ,Z 3 T ] T And standardizing each column vector in Z to obtain a new feature matrix
Figure FDA0003939444140000017
Wherein the upper symbol T represents the transpose of a matrix or vector, n = n 1 +n 2 +n 3 ,R n×10 A real number matrix representing n × 10 dimensions;
and 4, step 4: obtaining a classification transformation matrix W epsilon R by utilizing a Fisher discriminant analysis algorithm 10×2 And calculating a classification feature matrix
Figure FDA0003939444140000018
Sum residual matrix
Figure FDA0003939444140000019
Wherein the load matrix
Figure FDA00039394441400000110
And 5: and (3) carrying out singular value decomposition on the residual matrix E: e = UAV T And calculating a transformation matrix Q = VA -1 (ii) a Wherein U and V represent unitary matrix of singular value decomposition, the diagonal elements in diagonal matrix A are composed of singular values;
step 6: according to the formula psi = diag { UU T Calculating a detection index vector ψ, and determining a maximum value of ψ as ψ 0
And 7: the current data collected at the sampling moment are formed into a current data vector x t Where t represents the latest sampling instant, x is obtained by feature transformation t Mean value m of t Standard deviation of delta t Degree of kurtosis alpha t Deviation of degree beta t Root mean square upsilon t Crest factor gamma t Shape factor eta t Pulse factor e t Edge factor f t Maximum logarithm of g t (ii) a And constructing a feature vector z t =[m t ,δ t ,α t ,β t ,υ t ,γ t ,η t ,e t ,f t ,g t ]∈R 1×10 For feature vector z t Carrying out standardization processing to obtain new characteristic vector
Figure FDA00039394441400000111
And 8: using formulas
Figure FDA00039394441400000112
Computing residual vector E t And detecting the index psi t =E t E t T
And step 9: and judging whether the step 8 meets the condition: psi t ≤ψ 0 (ii) a If yes, the series arc fault does not occur, and the steps 7 to 9 are repeated to continue to detect the arc fault of the arc generating circuit; if not, detecting the series arc fault.
2. The method for detecting the series arc fault based on the discriminant analysis strategy as claimed in claim 1, wherein the specific implementation method for performing the feature transformation on the three types of current data obtained in step 1 in step 2 is as follows:
setting the window length as L and the moving step length as D,
step 2.1: initializing j =1;
step 2.2: mixing X 1 The current interval vector q is composed of the (j-1) × L +1 to the (j × L) th elements j ∈R 1×L (ii) a Wherein R is 1 ×L A real number vector representing 1 × L dimensions;
step 2.3: respectively calculating mean values m according to formulas (1) to (R) j Standard deviation delta j Degree of kurtosis alpha j Deviation of degree beta j Root mean square upsilon j Crest factor gamma j Shape factor η j Pulse factor e j Edge factor f j Maximum logarithm of g j
Figure FDA0003939444140000021
Figure FDA0003939444140000022
Figure FDA0003939444140000023
Figure FDA0003939444140000024
Figure FDA0003939444140000025
Figure FDA0003939444140000026
η j =υ j /m j
f j =q j (max)/m j
Figure FDA0003939444140000027
g j =log 10 (q j (max)) ⑩
Step 2.4: constructing a feature vector z j =[m j ,δ j ,α j ,β j ,υ j ,γ j ,η j ,e j ,f j ,g j ]∈R 1×10 And then judging whether the conditions are met: (j + 1). Times.L > N 1 (ii) a If not, after j = j +1 is set, returning to the step (2.2); if yes, all the obtained eigenvectors form a resistive load current characteristic matrix
Figure FDA0003939444140000028
3. The tandem arc fault detection method based on discriminant analysis strategy as claimed in claim 2, wherein X in step 2.2 1 Is replaced by X 2 、X 3 And respectively obtaining the inductive load current characteristic matrixes according to the implementation processes of the step 2.1 to the step 2.4
Figure FDA0003939444140000029
And a non-linear load current signature matrix
Figure FDA00039394441400000210
4. The method of claim 2, wherein the setting is based on a discriminant analysis strategy for series arc fault detectionThe latest sampling time is t, and the current data I acquired at the latest sampling time in the step 7 t And current data I of the first L-1 sampling moments t-1 ,I t-2 ,…,I t-L+1 Built into a current data vector x t =[I t ,I t-1 ,…,I t-L+1 ]∈R 1×L (ii) a Mixing X in step 2.2 1 Is replaced by x t Obtaining x t Mean value m of t Standard deviation of delta t Degree of kurtosis alpha t Deviation of degree beta t Root mean square upsilon t Crest factor gamma t Shape factor η t Pulse factor e t Edge factor f t Maximum logarithm of g t
CN202011168112.XA 2020-10-28 2020-10-28 Series arc fault detection method based on discriminant analysis strategy Active CN112345895B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011168112.XA CN112345895B (en) 2020-10-28 2020-10-28 Series arc fault detection method based on discriminant analysis strategy

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011168112.XA CN112345895B (en) 2020-10-28 2020-10-28 Series arc fault detection method based on discriminant analysis strategy

Publications (2)

Publication Number Publication Date
CN112345895A CN112345895A (en) 2021-02-09
CN112345895B true CN112345895B (en) 2023-02-10

Family

ID=74359266

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011168112.XA Active CN112345895B (en) 2020-10-28 2020-10-28 Series arc fault detection method based on discriminant analysis strategy

Country Status (1)

Country Link
CN (1) CN112345895B (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5459675A (en) * 1992-01-29 1995-10-17 Arch Development Corporation System for monitoring an industrial process and determining sensor status
US7366622B1 (en) * 2005-10-17 2008-04-29 X-L Synergy Arc fault identification using model reference estimation
CN103384805A (en) * 2011-01-31 2013-11-06 东北大学 Operation fault detection device for electric arc furnace and method thereof
CN106529079A (en) * 2016-11-29 2017-03-22 上海电机学院 Chemical process failure detection method based on failure-dependent principal component space
CN108181893A (en) * 2017-12-15 2018-06-19 宁波大学 A kind of fault detection method based on PCA-KDR
CN108490923A (en) * 2018-04-28 2018-09-04 南京航空航天大学 The design method of small fault detection and positioning for electric traction system
CN110221590A (en) * 2019-05-17 2019-09-10 华中科技大学 A kind of industrial process Multiple faults diagnosis approach based on discriminant analysis
CN110222553A (en) * 2019-03-29 2019-09-10 宁波大学 A kind of recognition methods again of the Multi-shot pedestrian based on rarefaction representation
CN111458599A (en) * 2020-04-16 2020-07-28 福州大学 Series arc fault detection method based on one-dimensional convolutional neural network

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8817431B2 (en) * 2009-12-18 2014-08-26 True-Safe Technologies, Inc. System and integrated method for a parallel and series arc fault circuit interrupter
CN106597260B (en) * 2016-12-29 2020-04-03 合肥工业大学 Analog circuit fault diagnosis method based on continuous wavelet analysis and ELM network

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5459675A (en) * 1992-01-29 1995-10-17 Arch Development Corporation System for monitoring an industrial process and determining sensor status
US7366622B1 (en) * 2005-10-17 2008-04-29 X-L Synergy Arc fault identification using model reference estimation
CN103384805A (en) * 2011-01-31 2013-11-06 东北大学 Operation fault detection device for electric arc furnace and method thereof
CN106529079A (en) * 2016-11-29 2017-03-22 上海电机学院 Chemical process failure detection method based on failure-dependent principal component space
CN108181893A (en) * 2017-12-15 2018-06-19 宁波大学 A kind of fault detection method based on PCA-KDR
CN108490923A (en) * 2018-04-28 2018-09-04 南京航空航天大学 The design method of small fault detection and positioning for electric traction system
CN110222553A (en) * 2019-03-29 2019-09-10 宁波大学 A kind of recognition methods again of the Multi-shot pedestrian based on rarefaction representation
CN110221590A (en) * 2019-05-17 2019-09-10 华中科技大学 A kind of industrial process Multiple faults diagnosis approach based on discriminant analysis
CN111458599A (en) * 2020-04-16 2020-07-28 福州大学 Series arc fault detection method based on one-dimensional convolutional neural network

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
SVD在串联故障电弧检测及选相中的应用;王泽伦等;《中国安全生产科学技术》;20200930(第09期);162-167 *
基于峰值指标和脉冲指标的串联电弧故障检测研究;刘德军等;《电器与能效管理技术》;20170830(第16期);11-15+43 *

Also Published As

Publication number Publication date
CN112345895A (en) 2021-02-09

Similar Documents

Publication Publication Date Title
Sarwar et al. High impedance fault detection and isolation in power distribution networks using support vector machines
CN107340456B (en) Power distribution network operating condition intelligent identification Method based on multiple features analysis
CN109670553B (en) Photovoltaic array fault diagnosis method based on adaptive neural fuzzy inference system
Ahmadi et al. A new method for detecting series arc fault in photovoltaic systems based on the blind-source separation
CN104181460B (en) A kind of permanent magnetic actuator vacuum circuit breaker method for diagnosing faults based on Multi-information acquisition
CN104914847B (en) Industrial process method for diagnosing faults based on direction core offset minimum binary
JP5049675B2 (en) Distribution line accident cause estimation system, method, and program
CN111679158A (en) Power distribution network fault identification method based on synchronous measurement data similarity
CN109061414A (en) Photovoltaic system DC Line Fault arc method for measuring
Hashizume et al. Fault Detection of Combinational Circuits Based on Supply Current.
Xing et al. A deep learning approach for series DC arc fault diagnosing and real-time circuit behavior predicting
CN117272143A (en) Power transmission line fault identification method and device based on gram angle field and residual error network
Moloi et al. Power distribution fault diagnostic method based on machine learning technique
CN112345895B (en) Series arc fault detection method based on discriminant analysis strategy
CN111999591A (en) Method for identifying abnormal state of primary equipment of power distribution network
CN113376474A (en) Neural network fault arc identification system and method based on generalized S transformation
CN108919041B (en) Transformer winding state online monitoring method based on cluster analysis
Singh et al. A novel methodology for identifying cross-country faults in series-compensated double circuit transmission lines
CN112444758B (en) Intelligent power distribution network line fault diagnosis and classification evaluation method
CN108509732A (en) Appraisal procedure, terminal device and the storage medium of steam turbine fault severity level
CN112858813A (en) Assessment method for lightning arrester characteristic distortion caused by high and low temperature factors
CN103514458A (en) Sensor fault distinguishing method based on combination of error correction codes and support vector machine
WO2020087128A1 (en) Solar farm fault detection and diagnosis
CN114692690A (en) Fault arc identification method based on multi-feature fusion support vector machine
CN112906489A (en) Transformer machine winding looseness identification method based on LMD and permutation entropy

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