1. Introduction
The application of photovoltaic technology is becoming more and more widespread in daily life, and the generation industry of photovoltaic power is at a stage of rapid development. In photovoltaic systems, including DC and alternating current (AC) systems, there have been considerable economic benefits from the wide use of various related equipment such as frequency converters, inverters, rectifiers, and charging piles. However, in the power grid, many power quality problems such as harmonic, noise, and distortion issues will be produced due to the imbalance, impact, and nonlinearity of the load in the circuit. The use of electrical appliances and the stability of the power supply network will be seriously affected by these power quality problems [
1,
2,
3]. Reference [
4] proposes an evaluation scheme for renewable solar energy, which is aimed at the risk–benefit assessment of historical and traditional buildings. In terms of the photovoltaic direct current, the identification and classification of disturbance in photovoltaic DC signals is the key and basis for power quality analysis, evaluation, and management [
5,
6], therefore the classification and identification of power quality disturbance has high research value. Reference [
5] proposed a comprehensive solution, which is a real-time analysis of power quality disturbances, and reference [
6] proposed a new type of framework for the recognition of complex power quality disturbances; this framework is based on the multiple fusion convolutional neural network (MFCNN), which could automatically extract and fuse signal features from multiple sources.
In the engineering field, in order to extract features from signals that contain noise, many time–frequency analysis methods [
7,
8,
9,
10] are used. Reference [
7] proposed a new method of emitter identification, which is based on spectral features and variational mode decomposition. In order to monitor family activities, reference [
8] proposed a feature extraction method of acoustic signals, which is based on the non-negative matrix factorization (NMF) classifier. Reference [
9] is a very classical algorithm, namely the empirical mode decomposition method (EMD), which was proposed by N. E. Huang et al. This method and the subsequent improved algorithm played a key part in signal feature extraction and signal decomposition. In order to extract time–frequency information features from vibration signals, wavelet packet transform was applied in reference [
10]. These adaptive time–frequency analysis methods for signal analysis can be roughly divided into two categories according to whether they are based on the Fourier transform. For example, the variational mode decomposition (VMD) algorithm is based on the Fourier transform. All of these methods realize signal decomposition by computing in the frequency domain. The essence of the VMD algorithm is an adaptive filter [
11], and the method was proposed by D. Zosso and K. Dragomiretskiy. VMD shows good performance [
12,
13,
14] compared with previous methods, especially in the analysis of complex, non-stationary signals, which has been proved in many references. Reference [
12] proposed a novel method of online chatter identification for milling processes. In order to make a prediction of short-term wind speed, reference [
13] combined linear and nonlinear prediction models and proposed a hybrid method based on variational mode decomposition. To detect the blade imbalance fault in the marine current turbine (MCT), reference [
14] proposed a novel VMD denoising method. The original signal was decomposed into several intrinsic mode functions, which was obtained by calculating the envelope of the extreme point; this is the essence of a method not based on the Fourier transform [
15], such as the ensemble empirical mode decomposition (EEMD), the local mode decomposition (LMD), and the empirical mode decomposition (EMD) methods. The EEMD algorithm was proposed by Z. Wu et al., and compared with EMD, EEMD improves the anti-noise ability of the EMD method by averaging the decomposition results and mixing different levels of white noise into the original signal [
16]. The VMD and EEMD methods are superior and effective [
17,
18], which can be shown by recent studies. Hence, the EEMD method was selected as a representative of the methods that are based on the non-Fourier transform, while the VMD method was selected as a representative of the methods that are based on the Fourier transform; these two algorithms were compared with the algorithm proposed in this paper.
Due to the incorporation and disconnection of linear loads in the circuit, an interference signal similar to a square wave signal is generated in the original photovoltaic DC signal, and this type of signal is generally sharp. Since signals with sharp corners, such as sawtooth signals and square wave signals, have frequency bands that are infinite and can be considered as wideband signals, the existing time–frequency analysis methods are expected show the Gibbs phenomenon when processing this kind of signal.
First, the essence of the algorithm based on the Fourier transform is multiscale adaptive filtering. However, after filtering, a series of interference may be caused at the breakpoint of the decomposition result; the characteristics of the sharp corner of the broadband signal disappearing or attenuating is known as the Gibbs phenomenon [
19]. For broadband signals, the use of an improved Hilbert–Huang transform (HHT) and the energy entropy algorithm [
20] to extract the features of alternating current (AC) signal square waves was proposed by K. He et al. Second, the essence of the methods that are based on the non-Fourier transform is the use of the interpolation function to calculate the envelope of the extreme point, and the division of the original signal into a number of intrinsic mode functions (IMFs) with a “smooth” narrow band. Therefore, when dealing with broadband components, errors will inevitably occur. Consequently, Y. Peng et al. proposed the Broadband Mode Decomposition (BMD) algorithm [
21,
22]. The innovative point of the BMD algorithm is to construct an association dictionary, which contains common broadband data such as narrowband data, square wave data, and sawtooth data. The sparse solution is then obtained by using the optimization method to search in the association dictionary. Compared with the previous signal decomposition algorithms, the BMD algorithm takes the adjustment differential operator as the optimal object; it can also extract signal features from complex signals, including narrowband signal features and broadband signal features. However, when applied to a wideband signal, the BMD algorithm may treat it as several narrowband components because its relative bandwidth of is not sufficiently small. Therefore, to denoise photovoltaic DC signals, a broadband mode decomposition (MBMD) method is proposed in this paper, which is based on the modulation differential operator. In order to make the decomposition result more accurate on the basis of the BMD algorithm, a high-frequency, single-frequency signal is multiplied by the MBMD so that the relative bandwidth of the effective wideband signal is converted into one that is sufficiently small and far lower than 1, and the wideband signal is then regarded as an approximate wideband signal.
Recently, J. Zheng et al. proposed the composite multiscale fuzzy entropy (CMFE) algorithm, which is based on the multiscale fuzzy entropy (MFE) algorithm. It has been proved to be suitable for processing non-stationary signals in engineering [
23,
24]. Therefore, it is applied to construct the feature vectors of the MBMD decomposition results. Then, the BP neural network algorithm is used to identify and classify the photovoltaic DC current disturbance.
The remaining sections of the paper are as follows.
Section 2 introduces the disturbance state in the photovoltaic DC current and establishes a photovoltaic DC mathematical model.
Section 3 introduces the problems in the identification of photovoltaic DC disturbances.
Section 4 describes the details of the MBMD algorithm.
Section 5 explains the other algorithms used.
Section 6 describes the simulation analysis of the feature extraction of the photovoltaic signal disturbance model.
Section 7 analyzes the collected experimental data sets of photovoltaic direct current signals.
Section 8 provides the conclusion. All the data in this paper come from the data collected by the photovoltaic experimental platform in the literature [
25].
2. Photovoltaic Electrical Signal Model
According to the analysis of the measured signals, the mathematical model of the DC signal is initially established, which consists of the DC signal, a harmonic signal superimposed on the DC, a distortion signal caused by the impact load, a noise signal caused by external interference, electromagnetic interference, etc.
2.1. Harmonic Signal
In normal operation, the current of the DC system contains abundant harmonic signals. The harmonic signal refers to the AC component superimposed on the DC. Therefore, the established harmonic signal model is as follows:
where,
is the amplitude of the DC current,
k is the amplitude coefficient of harmonics, and
n is the maximum number of harmonics.
2.2. Distortion Signal
Due to the impact, nonlinearity, and imbalance of the load, various disturbance components in the photovoltaic DC system increase. These disturbance components are distortion signals. For the photovoltaic DC system, there are two kinds of disturbances that have the greatest impact on photovoltaic power quality: current mutation and current pulse. The sudden change of current refers to the steep rise and fall of the current caused when the load switch in the system is opened and closed. For linear loads such as light bulbs, when the load switch is opened or closed, the current change is mostly straight up and down; for nonlinear loads such as motors, when the load switch is closed or opened, the current changes are characterized by slopes that rise and fall. This is determined by the characteristics of the nonlinear load.
Therefore, the linear load current mutation signal model is established as follows:
where
A is the amplitude coefficient of disturbance,
is the initial time of disturbance, and
is the end time of disturbance.
The nonlinear load current mutation signal model is the following:
where
and
are the slope of the signal,
and
are the amplitudes of signal change,
b is the duration of the signal rise or fall,
is the initial time of disturbance, and
is the end time of disturbance.
A current pulse is produced when faults such as a grounding fault or an automatic restart occur in the operation of the photovoltaic DC system, and refers to the current generally rising quickly and then falling quickly with a large amplitude change, but with a relatively short duration, generally 10–20 ms.
Therefore, the current pulse signal model is established as:
where
B is the amplitude coefficient of disturbance,
is the initial time of disturbance,
is the end time of disturbance, and
is the frequency of the pulse signal.
6. Simulation Analysis
The simulation signal of a photovoltaic DC current disturbance is established, as shown in
Figure 6. The second harmonic with an amplitude of 0.02 A and a frequency of 100 Hz is set in the signal, where disturbance
A is the disturbance generated by linear load access or connection; disturbance
B is the disturbance caused by nonlinear load access or connection; disturbance
C is the current pulse disturbance, and
is white noise with an SNR of 25. The signal model formula is as follows:
When analyzing broadband signals, the VMD algorithm is still affected by the Gibbs phenomenon because its essence is the multi-scale adaptive Wiener filtering.
Figure 7 shows the decomposition results of the photovoltaic DC simulation electrical signals by VMD.
It can be seen in
Figure 7a that the mutation information of the square wave signal in the
IMF1 component decomposed by VMD is lost, the Gibbs phenomenon appears, and the missing characteristic information appears separately in the remaining IMF components.
Figure 7d is the residual component.
To obtain the envelope of the extreme points, the EMD method and the EEMD method typically use interpolation functions, which are non-Fourier transform-based methods, thereby generating a smooth envelope. When dealing with the broadband component, since the IMF is obtained by calculating the envelope, the final decomposition results are also smoothed.
Figure 8 shows the decomposition result of the photovoltaic current simulation signal using EEMD. It can be seen that the obtained
IMF1 and
IMF2 are components with partial signal characteristics, the
IMF3 component is close to a sinusoidal signal, and
IMF4 is close to the DC component. Evidently, the extracted signal characteristics have changed.
It can be seen in
Figure 9 that the MBMD method is used to decompose the DC current signal, which can separate and extract the DC signal, pulse signal, distortion and harmonic signals, and the noise signal.
Table 1 lists the accuracy parameters of the separated IMF obtained by the three algorithms, including the accuracy parameters of the correlation coefficient
, the energy error
, and the time
T to make a more accurate comparison.
In order to compare the effect of signal decomposition under different signal-to-noise ratios, two additional simulation models are constructed here. Except for noise, the other components remain unchanged. White noise
with an SNR of 10 and white noise
with an SNR of 5 were added, respectively. The two signal models constructed are shown in
Figure 10 and
Figure 11.
Figure 12 and
Figure 13 show the VMD decomposition results of these two signal models. It can be seen that with the decrease of SNR, the components of the VMD containing noise increase, especially
IMF2 and
IMF3. Due to the increase of noise, the decomposed components also contain part of the noise signals, so that the characteristics of the square wave signal are not obvious.
Figure 14 and
Figure 15 are the decomposition results of EEMD. Some components obtained by EEMD are relatively smooth, and the increased noise components are reflected in
IMF1.
Figure 16 and
Figure 17 shows the decomposition results of MBMD.
Table 2 and
Table 3 show the relevant parameters (energy error
, correlation coefficient
, and time T), respectively. It can be seen that the decomposition effect of MBMD is the best under different SNRs.