Patent Description:
Nowadays, electrical fires account for a high proportion of fire accidents. An arc fault is one of the main causes of the electrical fires. The fault arc is usually caused by aging and breakage of insulation of lines and equipment, or poor electrical connections. When an arc fault occurs, the arc may release high heat which can easily cause fire. Based on the type of the arc fault, the arc fault includes a series arc fault or a parallel arc fault. In a case of the parallel arc fault, the current is large, and the conventional overcurrent protection device and short-circuit protection device may perform partially protection. In a case of the series arc fault, the current is abnormal and is lower than a protection threshold, and the conventional overcurrent protection device cannot effectively detect the fault arc and protect a circuit.

With the conventional arc detection method, a threshold is set for an extracted feature value to determine whether a half wave is an arc. Due to the diversity of loads in the actual electricity environment, different thresholds are set for different loads. Therefore, the conventional arc detection method cannot adapt to different load environments, and it is still required to improve the performance of the conventional arc detection method. The patent <CIT> discloses a low-voltage power distribution system series fault arc identification method based on all-phase deep learning, including: under a low-voltage alternating current system, carrying out current signal collection on different loads in a low-voltage loop; carrying out all-phase discrete Fourier transform on a collected current signal, carrying out full-phase spectrum characteristic quantity extraction of a load, and constructing an all-phase spectrum characteristic vector; constructing a deep learning neural network model based on Logistic regression, carrying out deep learning training on all-phase spectrum characteristic quantity. The patent <CIT> discloses an electrical load determination apparatus and method using a modified MFCC, in which an electrical load-specific feature vector showing a clear difference in an electrical load-specific current signal waveform can be efficiently extracted by using a predetermined number of filters of frequency bands that differ from the frequency bands used in a conventional MFCC algorithm, and processing time and size can be reduced by performing a modified MFCC algorithm. The patent <CIT> discloses an arc light fault identification device and method based on panoramic information. By analyzing the characteristic changes of a current in a time domain and a frequency domain when an arc fault occurs, appropriate characteristic values are extracted as input of a BP neural network, and the BP neural network is adopted to detect the arc fault. The patent <CIT> discloses a DC arc fault detection method based on the wavelet transform and the CNN (convolutional neural network), including: (<NUM>) collecting an arc current signal of a to-be-detected line, and performing high-pass filtering to filter the DC and the low-frequency noise signal in the original collected signal; (<NUM>) calculating a wavelet coefficient of the arc current signal of the to-be-detected line by using the wavelet transform; and (<NUM>) inputting the normalized wavelet coefficient of the arccurrent signal of the to-be-detected line into a trained CNN. The patent <CIT> discloses a DC arc detection method based on impedance characteristics and a DC arc detection system thereof. The method comprises the steps that step <NUM>, an arc resistance value is determined in advance according to a DC arc nonlinear resistance model; step <NUM>, the digital signals of the DC voltage and the DC current of a DC series power supply system are acquired according to the preset frequency; and step <NUM>, the mean resistance in the sampling period is calculated according to the digital signals. The patent <CIT> discloses a series fault arc detection method, including: classifying an arc generated in a low-voltage power supply and distribution system into a possible fault arc and a possible non-fault arc according to a current signal collected in the low-voltage power supply and distribution system; adopting a depth convolution neural network method to identify the possible faultarc and the possible non-fault arc and determining the non-fault arc or the fault arc.

In view of the defects of the conventional technology, a method for detecting a fault arc signal by using a convolutional neural network is provided according to the embodiments of the present disclosure. Compared with the conventional detection method in which a threshold is set, the method according to the present disclosure has higher accuracy and higher reliability.

In a first aspect, a method for detecting a fault arc signal by using a convolutional neural network is provided in an embodiment of the present disclosure. The method includes: S1, filtering, by using three band-pass filters with different pass-bands, a sampled current signal to extract time-frequency eigenvectors; S2, constructing a three-dimensional matrix based on the time-frequency eigenvectors; S3, constructing a two-dimensional convolutional neural network model, and training the two-dimensional convolutional neural network model to obtain a two-dimensional convolutional neural network model with best performance; S4, performing, by using the two-dimensional convolutional neural network model with the best performance, an online determination on the three-dimensional matrix to obtain an arc detection result, where an arc detection result of <NUM> indicates that no arcing occurs, and an arc detection result of <NUM> indicates that an arcing occurs; S5, counting the number of fault half-waves in an observation time period ΔT based on the arc detection result; and S6, comparing the number of the fault half-waves in the observation time period ΔT with a threshold, performing a tripping operation in a case that the number of the fault half-waves exceeds the threshold, and performing no operation in a case that the number of the fault half-waves does not exceed the threshold.

In step S1, an AD conversion is performed on the sampled current signal, and the sampled signal passes through the three band-pass filters with different pass-bands. The pass-bands respectively range from <NUM> to <NUM>, from <NUM> to <NUM>, and from <NUM> to <NUM>. <NUM> time-frequency eigenvectors are extracted from each of output signals of the filters, and each of the eigenvector is a <NUM>-dimensional eigenvector.

In the method, in performing time-frequency feature extraction on the filtered signal, segmentation processing is performed. A <NUM> half-wave is divided evenly into <NUM> segments, and a waveform of each of the segments is analyzed and processed. Extracted similar eigenvalues of each of the segments form a time series.

Specifically, each of the time-frequency eigenvectors includes a time-domain feature and a frequency-domain feature, and the time-domain feature includes a time dispersion, an amplitude dispersion, and the number of waveforms. In extracting the time-domain feature, for each of the waveforms, the waveform is preprocessed, a non-local extreme point of the waveform is eliminated, remaining sampling points are connected sequentially to obtain a new waveform, and the time-domain feature is extracted based on the new waveform.

The time dispersion is calculated by using the following equation: <MAT> where Ti represents a time interval between two adjacent minimum values.

The amplitude dispersion is calculated by using the following equation: <MAT> where each of a difference between VFH and VDF, a difference between VHJ and VFH and a difference between VHJ and VFH in the numerator represents an amplitude difference between two adjacent minimum points, and each of VDE, VFG and VMP in the denominator represents an amplitude difference between a minimum point and a maximum point adjacent to the minimum point.

The number of waveforms is calculated as follows. Assuming that y represents a sequence of the new waveform obtained by preprocessing the waveform, the number N of the waveforms is calculated by using the following equation: <MAT> where length(y) represents a sequence length of the preprocessed waveform, and <MAT> represents a rounding down operation.

The above calculation is performed on <NUM> segments of each of the half-waves, and extracted <NUM> eigenvalues chronologically form a <NUM>-dimensional vector.

In the method, a frequency-domain feature of the filtered signal is extracted by: step <NUM>, assuming that L represents a length of data in each of time period, and performing a <NUM>-point FFT operation on the data in each of the time period for <MAT> times; step <NUM>, selecting <NUM> feature frequency points from an FFT operation result, obtaining a vector based on M FFT operation results, and performing median filtering respectively on M-dimensional eigenvectors at the <NUM> frequency points; step <NUM>, summing the <NUM> filtered M-dimensional eigenvectors to obtain <NUM> eigenvalues; and step <NUM>, performing the above operations on <NUM> time periods in each of the half-waves, and obtaining a eigenvector based on of the eigenvalues in different time periods in a time sequence.

In the method, before performing processing by using the neural network, it is required to perform normalization on the extracted eigenvectors. The normalization is performed by using the following equation: <MAT> where x[n] represents an n-th element in an eigenvector, and x[n] represents an element after normalization; X represents an eigenvector before normalization; max(X) represents a maximum element in the eigenvector X; and min(X) represents a minimum element in the eigenvector X.

In the method, the eigenvectors are spliced. <NUM> eigenvectors are extracted from the output of each of the three filters, and the eigenvectors corresponding to each of the filters form a <NUM>*<NUM> matrix. Further, matrices corresponding to the three filters are stacked to form a <NUM>*<NUM>*<NUM> three-dimensional feature matrix, where the last dimension <NUM> represents the number of channels.

In an embodiment, a topology structure of the convolutional neural network used in the detection method mainly includes two convolution layers. A first convolution layer has three <NUM>*<NUM>*<NUM> convolution kernels, where the number <NUM> in <NUM>*<NUM>*<NUM> represents the number of channels. A second convolution layer has five <NUM>*<NUM>*<NUM> convolution kernels. The topology structure of the convolutional neural network further includes a MaxPooling2D pooling layer, a Dropout layer, and a Flatten layer that stretches a multi-dimensional output of a convolution layer to a one-dimensional vector. The one-dimensional vector is inputted to a fully connected layer. The fully connected layer includes two layers which respectively have <NUM> neurons and <NUM> neurons. An output layer has one neuron.

In performing calculation in a convolution layer, an input matrix is multiplied by corresponding points in a convolution kernel of the convolution layer, then products are summed, and then a bias value is added. The calculation is performed by using the following equation: <MAT> where K represents the number of channels, M represents the number of rows of a convolution kernel in each of the channels, N represents the number of columns of the convolution kernel in each of the channels, yn represents a convolution output result, bn represents a direct-current bias in a linear operation, ai, j, k represents a weighting coefficient in the linear operation, and xi, j, k represents an originally inputted feature element or an output result of a previous convolution layer.

Assuming that an original matrix is an A*B*K matrix and the convolution operation is performed with a stride of <NUM>, then a (A-M+<NUM>)*(B-N+<NUM>) result is outputted after the above convolution operation. The number of the channels of the convolution operation result is determined by the number of the convolution kernels.

Dimensionality reduction is performed on the convolution result through the pooling layer. In the method, a MaxPooling2D pooling layer is used to extract a maximum value from data in an region in a channel. The pooling operation is performed by using the following equation: <MAT> where xi,j represents an elements in a region covered by a pooling window, and an represents a pooling output result representing a maximum xi,j in the channel.

In an embodiment, the convolution layers and the fully connected layers adopt a ReLu activation function, which is expressed as: <MAT> and the output layer adopts a sigmoid activation function, which is expressed as: <MAT> The output layer outputs a probability value ranging from <NUM> and <NUM> by using the sigmoid activation function. A half-wave is determined as a normal half-wave or as a fault arc based on the probability value. For example, if <NUM> is set as a threshold, the half-wave is determined as an arc in a case that the probability value is greater than <NUM>, and the half-wave is determined as a normal half-wave in a case that the probability value is less than <NUM>.

In an embodiment, in the detection method, before the neural network model is used for determination, it is required to train the neural network model offline based on training data to obtain and save a model with best performance. Then, online determination is performed on the obtained feature matrix by using the trained model.

In an embodiment, in the detection method, experimental data labeled as normal data and experimental data labeled as arcing data are required for training the neural network. In collecting arcing experimental data in the laboratory, there may be a case in which data labeled as arcing data is collected and no arcing occurs, thus it is required to eliminate the mislabeled data in this case. In a series arc experiment, a voltage across an arc generator or a voltage across a carbonized cable is measured, then the mislabeled data is eliminated based on the measured voltage, that is, the data labeled as arcing data while no arcing occurring is eliminated. In a parallel arc experiment, a current at a position at which an arc occurs in a cable is measured, and then the data labeled as arcing data while no arcing occurring is eliminated based on the measured current.

With the detection method, it is determined whether a <NUM> half-wave is an arc. In a practical arc detection, it is required to comprehensively process a half-wave in an observation time period ΔT to determine whether an arcing occurs in this time period. Counting is performed by using the following equation: <MAT> where yi represents a determination result of an i-th half-wave in the observation time period, a determination result equal to <NUM> indicates that the half-wave is a normal half-wave, and a determination result of equal to <NUM> indicates that an arcing occurs. The number of fault half-waves in the observation time period is obtained by summing vectors, and then is compared with a fault half-wave threshold corresponding to ΔT to determine whether to perform tripping operation.

Compared with the conventional method in which a single eigenvalue is obtained and then the single eigenvalue is compared with a threshold to determine a half-wave is a fault arc half-wave, with the method based on a convolutional neural network according to the present disclosure, a higher accuracy and higher reliability can be achieved in identifying a fault arc half-wave, and adaptability can be achieved in performing training for different load conditions.

In order for a clearer illustration of technical solutions in the embodiments of the present disclosure or the conventional technology, accompanying drawings used in the description of the embodiments or the conventional technology are briefly introduced below.

The technical solutions in the embodiments of the present disclosure are clearly and completely described below with reference to the accompanying drawings used in the embodiments of the present disclosure. Apparently, the described embodiments are a part, rather than all, of the embodiments of the present disclosure. Any other embodiments obtained by those of ordinary skill in the art from the embodiments of the present disclosure without any creative effort shall fall within the protection scope of the present disclosure.

It is to be understood that, when used in this specification and the appended claims, terms "comprise" and "include" (and variants thereof) indicate existence of described features, entireties, steps, operations, elements and/or components, and do not exclude existence or addition of one or more of other features, entireties, steps, operations, elements, components, and/or combinations thereof.

According to the conventional technology, when a fault arc occurs, a current signal in the circuit may be significantly distorted, and a voltage signal is similar to a normal voltage signal. Therefore, in the present disclosure, a current signal is sampled, then AD conversion is performed on the current signal, and then the current signal is analyzed to obtain features of various arc signals.

In addition, for the problems of the conventional arc detection method described in the background technology, with the solutions according to the present disclosure, an analog-to-digital conversion is performed on a sampled current signal, and then filtering is performed by using three filters with different pass-bands. For each of outputted half-wave signals after filtering, time-domain eigenvectors and frequency-domain eigenvectors of the half-wave are extracted. Eigenvectors corresponding to an output of a same filter are spliced to obtain a two-dimensional matrix. Feature matrices corresponding to the three filters are stacked to obtain a three-dimensional feature matrix. A two-class processing is performed on the three-dimensional feature matrix by using a two-dimensional convolutional neural network, and it is determined whether an arc occurs in the half-wave based on an outputted probability value. The number of half-waves of fault arcs occurring in an observation time period ΔT is counted, and is compared with a preset threshold. A tripping operation is performed in a case that the number of the half-waves of the fault arcs occurring in the observation time period ΔT exceeds the preset threshold, and no operation is performed in a case that the number of the half-waves of the fault arcs occurring in the observation time period ΔT does not exceed the preset threshold.

Hereinafter, the method for detecting a fault arc according to the present disclosure is described with reference to <FIG>.

<FIG> shows a flowchart of a method for detecting an arc based on a convolutional neural network according to the present disclosure. First, an analog-to-digital conversion is performed on a sampled current signal, and then the signal after the analog-to-digital conversion is filtered by three band-pass filters with different pass-bands. The pass-bands respectively range from <NUM> to <NUM>, from <NUM> to <NUM>, and from <NUM> to <NUM>. Each of filtered half-waves with a time length of <NUM> is equally divided into <NUM> segments, and an arc eigenvalue of a high-frequency signal in each of the segments is extracted respectively in a time domain and in a frequency domain. Eigenvalues of a same type of the <NUM> segments are arranged chronologically to form a <NUM>-dimensional eigenvector. Then, eigenvectors of different types form a feature matrix. The feature matrix is processed by using a multi-channel two-dimensional convolutional neural network. Before using the neural network for online determination, it is required to train the neural network to obtain and save an optimal model. The number of fault half-waves in an observation time period ΔT is counted and then is compared with a preset threshold to determine whether to perform a tripping operation.

The above process is described in detail below.

In detecting an arc signal, since the arc signal is non-stationary, each of the half-wave signals is processed by time segments. Based on this idea, in a preferred embodiment, a <NUM> half-wave is divided into <NUM> segments, a time-domain feature and a frequency-domain feature of each of the segments are extracted, and eigenvalues extracted from the <NUM> segments are arranged chronologically to form a <NUM>-dimensional eigenvector.

It is required to perform waveform preprocessing based on the time dispersion, the amplitude dispersion, and the number of waveforms in the time-domain feature. The waveform preprocessing is performed by: for each of the original waveforms, eliminating a non-local extreme point of the waveform, and connecting remaining sampling points in sequence to obtain a new waveform. The new waveform includes only local maximum points and local minimum points in the original waveform. <FIG> shows waveforms after eliminating non-extreme points based on the above description.

In performing time-domain feature analysis, the calculation of the time dispersion is shown in <FIG>. Based on the preprocessed waveform, the calculation is performed by dividing a sum of absolute values of time differences between adjacent waveforms by a sum of the time periods. The calculation is performed by using the following equation: <MAT>.

As shown in <FIG>, in the above equation, Ti represents a time interval between two adjacent minimum values, and represents a time length of a waveform unit.

In performing time-domain feature analysis as shown in <FIG>, the amplitude dispersion is calculated by: dividing a sum of absolute values of amplitude differences of adjacent waveforms by a sum of amplitudes in the segment. The calculation is performed by using the following equation: <MAT> where each of a difference between VFH and VDF, a difference between VHJ and VFH and a difference between VHJ and VFH in the numerator represents an amplitude difference between two adjacent minimum points, and each of VDE, VFG and VMP in the denominator represents an amplitude difference between a minimum point and a maximum point adjacent to the minimum point.

In an embodiment, it is assumed that y represents a sequence of the new waveform obtained by preprocessing the waveform, then the number N of the waveforms is calculated by using the following equation: <MAT> where length(y) represents a sequence length of the preprocessed waveform, and <MAT> represents a rounding down operation.

The frequency-domain feature is extracted from the filtered signal. Same frequency-domain processing is performed on each of signals outputted from the three filters. The FFT eigenvalues are extracted by performing the following steps <NUM> to <NUM>.

In step <NUM>, data of each of the <NUM> half-waves filtered by different sub-band filters is divided into <NUM> segments.

In step <NUM>, a <NUM>-point FFT transform is performed on data in each of the segments. Assuming that L represents a length of the data in each of the segments, the <NUM>-point FFT transform is performed on the data in each of the segments for <MAT> times.

In step <NUM>, <NUM> frequency channels are selected from FFT operation results corresponding to two pass-bands, and FFT transform values of M same frequency points in a segment form an M-dimensional eigenvector.

In step <NUM>, median filtering is performed on M-dimensional eigenvectors respectively corresponding to <NUM> frequency points to obtain <NUM> median filtered eigenvectors.

In step <NUM>, the median filtered eigenvectors corresponding to <NUM> frequency points are summed according to the frequency points to obtain eigenvalues corresponding to <NUM> frequency points in each of the segments.

In step <NUM>, the above operations are performed on each of the <NUM> segments of each of the half-waves, and eigenvalues of the <NUM> segments in a same FFT channel form an eigenvector. The eigenvectors corresponding to the <NUM> frequency points form a <NUM>*<NUM> feature matrix.

Described above are eigenvalues designed and adopted in a preferred embodiment of the present disclosure, and the eigenvalues to be processed by the method for detecting an arc based on the multi-channel two-dimensional convolutional neural network are not limited to the eigenvalues mentioned above.

With the time-domain feature analysis and the frequency-domain feature analysis, multiple eigenvectors may be obtained. Before processing the eigenvectors by using the neural network, it is required to perform normalization on the eigenvectors to eliminate the influence of the dimensions of different eigenvalues. Since each of the half-waves is divided into <NUM> segments and the obtained eigenvectors are <NUM>-dimensional, normalization is performed on each of the eigenvectors by using the following equation: <MAT> where x[n] represents an n-th element in the eigenvector, and x[n] represents an element after normalization; X represents the eigenvector; max(X) represents an element with a maximum value in the eigenvector X; and min(X) represents an element with a minimum value in the eigenvector X.

In an embodiment, based on the image processing method by using a convolutional neural network, the eigenvectors after normalization are spliced to obtain a feature matrix in the detection method by using the neural network according to the present disclosure. Two-dimensional feature matrices corresponding to different filters are similar to different channels in an image. For each of signals outputted from the filters, <NUM> time-domain eigenvectors and <NUM> eigenvectors are extracted, each of the half-waves is divided into <NUM> segments, and the eigenvectors are spliced to obtain a <NUM>*<NUM> feature matrix. The two-dimensional feature matrices corresponding to the three filters may be stacked. As shown in <FIG>, three <NUM>*<NUM> matrices are stacked to obtain a <NUM>*<NUM>*<NUM> three-dimensional matrix, where <NUM> indicates the number of channels of the feature matrix.

In the embodiment, a topology structure of the neural network is shown in <FIG>. With reference to the process of the neural network shown in <FIG>, the topology structure of the neural network is briefly described in the following steps <NUM> to <NUM>.

In step <NUM>, a <NUM>*<NUM>*<NUM> three-dimensional feature matrix corresponding to each of half-waves is inputted through an input layer, and then is processed by two convolution layers.

In step <NUM>, a first convolution layer has three <NUM>*<NUM>*<NUM> convolution kernels, where the number <NUM> in <NUM>*<NUM>*<NUM> indicates that the number of the convolution kernel is same as the number of channels of the inputted feature matrix. Each of the convolution kernels outputs a <NUM>*<NUM> result. The three convolution kernels of the first convolution layer correspond to three channels. Thus, the first convolution layer outputs a <NUM>*<NUM>*<NUM> result.

In step <NUM>, a first pooling layer having a <NUM>*<NUM> pooling window performs dimension reduction on the output of the first convolution layer to output a <NUM>*<NUM>*<NUM> result.

In step <NUM>, the output of the first pooling layer is inputted to a second convolution layer which has five <NUM>*<NUM>*<NUM> convolution kernels, and a <NUM>*<NUM>*<NUM> result is outputted.

In step <NUM>, a pooling layer having a <NUM>*<NUM> pooling window performs dimension reduction on the result outputted from the second convolution layer, and output a <NUM>*<NUM>*<NUM> result.

In step <NUM>, a Flatten layer stretches the three-dimensional matrix to obtain a one-dimensional vector including <NUM> elements.

In step <NUM>, the one-dimensional vector is inputted to a fully connected layer having <NUM> neurons, and then is inputted to a fully connected layer having <NUM> neurons, and then is inputted to an output layer having one neuron.

In step <NUM>, after the neuron in the output layer perform processing, the output layer outputs a probability value for performing two-class processing to determine whether an arc occurs or no arc occurs.

The Dropout layer in <FIG> is mainly used in training to reduce overfitting in training.

In the embodiment, a multi-channel two-dimensional convolution operation is performed. In performing the multi-channel two-dimensional convolution operation, a two-dimensional convolution operation is performed for each of channels, then a convolution result for each of the channels are summed, and then a bias value is added. The calculation is performed by using the following equation: <MAT> where K represents the number of channels, M represents the number of rows of a convolution kernel in each of the channels, N represents the number of columns of the convolution kernel in each of the channels, yn represents a convolution output result, bn represents a direct-current bias in a linear operation, ai, j, k represents a weighting coefficient in the linear operation, and xi, j, k represents an originally inputted feature element or an output result of a previous convolution layer.

In the embodiment, the convolution operation is performed with a stride of <NUM>. As shown in <FIG>, in performing a next convolution operation, a row sliding is performed on the matrix inputted to the convolution layer according to the stride. First, the row is fixed, and a column sliding is performed until sliding to the end of the column, and then a row sliding is performed along the direction of the row according to the stride. It is assumed that an original matrix is an A*B*K matrix, where K represents the number of channels in the data matrix. A convolution operation with an M*N*K convolution kernel is performed, then a (A-M+<NUM>)*(B-N+<NUM>) result is outputted. The number of the channels of the convolution operation result is determined by the number of the convolution kernels.

In the embodiment, dimensionality reduction is performed on the convolution result by using a pooling layer, a MaxPooling2D pooling layer is used. As shown in <FIG>, it is assumed that a convolution result has a <NUM>*<NUM> channel matrix and a <NUM>*<NUM> pooling window is used in the pooling process, the pooling process is performed by using the following equations: <MAT>.

Thus, the pooling process outputs a <NUM>*<NUM> result. During the pooling process, adjacent pooling operation windows do not overlap with each other.

In the embodiment, the outputs by the neurons in the fully connected layers and the output layer are obtained by using the following calculation: <MAT> where yn represents an output of the fully connected layer or the output layer after performing a linear operation, ai represents a weighting coefficient for operations in the fully connected layer or the output layer, xi represents an input to the fully connected layer or the output layer, and bn represents a direct current bias in the linear operation.

N is <NUM> for the calculation of the neuron in a first fully connected layer. N is <NUM> for the calculation of the neuron in the second fully connected layer. N is <NUM> for the calculation of the neuron in the output layer.

In the embodiment, the convolution layers and the fully connected layers adopt a ReLu activation function, which is expressed as: <MAT> where x represents a weighted sum result after convolution operations or a weighted sum result after processing by the fully connected layer.

In the embodiment, the output layer adopts a sigmoid function, which is expressed as: <MAT> where x represents a weighted sum result of a last fully connected layer. A result outputted by the output layer ranges from <NUM> and <NUM>, which represents a probability of a classification result being <NUM> or <NUM>.

In the embodiment, a result outputted by an activation function of the output layer is classified based on a threshold of <NUM>, which is expressed as: <MAT> where x represents an output of a neuron in the output layer; and y represents a determination result of a half-wave, where that the half-wave is a normal half-wave in a case that y=<NUM>, and the half-wave is determined as a fault arc half-wave in a case that y=<NUM>.

In the detection method, before the neural network model is used for determination, it is required to train the neural network model offline based on training data to obtain and save a model with best performance. Then, online determination is performed on the obtained feature matrix by using the trained model. In collecting data in a laboratory, there may be a case in which data is labeled as arcing data while no arc occurs on an arc generator or a carbonized cable, thus it is required to clean and eliminate data before the data is provided to the neural network for training. According to the present disclosure, a voltage at a position at which a series fault arc occurs and a current at a position at which a parallel fault arc occurs are measured to determine whether the collected experimental data indicates an arcing occurs. <FIG> and <FIG> show the circuits.

In the series arc experiment, there are two cases in which there is no arcing. In one case, an iron rod in the arc generator is completely separated from a carbon rod in the arc generator, or two wires in the carbonized cable are separated from each other. In this case, the voltage across the arc generator or the carbonized cable is a standard line voltage, as shown by line c in <FIG>. In the other case, an iron rod in the arc generator is in complete contact with a carbon rod in the arc generator, or two wires in the carbonized cable are connected to each other. In this case, the voltage across the arc generator or the carbonized cable fluctuates within a small range around zero, as shown by line b in <FIG>. When an arc occurs, the voltage is lower than the standard line voltage and is seriously distorted, as shown by line a in <FIG>. Therefore, the experimental data mislabeled as arcing data may be eliminated based on a waveform of the measured voltage.

In the parallel arc experiment, there are two cases in which there is no arcing. In one case, two wires in the cable are separated from each other, and a current at a position at which the arc occurs is close to zero, as shown by line c in <FIG>. In the other case, two wires in the cable are short-circuited, and a current at a position at which the arc occurs is very large, as shown by line b in <FIG>. When an arc occurs, the current is less than a line conduction current and there is a flat shoulder feature which is an arc symbolic feature, as shown in line a in <FIG>. Due to the short acquisition time period for the parallel arc, an arc may be determined by manually checking whether the current has a flat shoulder feature.

In an embodiment, a breaking time of a circuit breaker varies with a current. Therefore, in addition to perform determination on the half-wave, it is required to perform determination on all half-waves in the observation time period ΔT by using the neural network, and it is determined whether to perform a tripping operation based on a determination result of the half-waves in the observation time period.

In an embodiment, the determination is performed by performing the following steps <NUM> to <NUM>.

In step <NUM>, an observation time period ΔT and a fault half-wave number threshold in the observation time period are determined by querying a table based on a calculated measurement current.

In step <NUM>, half-waves in the ΔT are detected and determined by using detection method based on the neural network, determination results are outputted, and a determination result vector is obtained based on the determination results.

In step <NUM>, elements in the determination result vector are summed up to obtain the number of fault half-waves in the observation time period. The calculation is performed by using the following equation: <MAT> where yi represents a determination result of an i-th half-wave in the observation time period, a determination result equal to <NUM> indicates that the half-wave is a normal half-wave, and a determination result of equal to <NUM> indicates that an arcing occurs; <MAT> represents the number of half-waves in the observation time period ΔT, and <MAT> presents a rounding up operation.

In step <NUM>, the number of the half-waves in the observation time period ΔT obtained above is compared with a threshold to determine whether to perform a tripping operation.

Claim 1:
A method for detecting a fault arc signal by using a convolutional neural network, comprising:
S1, filtering, by using three band-pass filters with different pass-bands, a sampled current signal to extract time-frequency eigenvectors;
S2, constructing a three-dimensional matrix based on the time-frequency eigenvectors;
S3, constructing a two-dimensional convolutional neural network model, and training the two-dimensional convolutional neural network model to obtain a two-dimensional convolutional neural network model with best performance; and
S4, performing, by using the two-dimensional convolutional neural network model with the best performance, an online determination on the three-dimensional matrix to obtain an arc detection result, wherein an arc detection result of <NUM> indicates that no arcing occurs, and an arc detection result of <NUM> indicates that an arcing occurs.