Multiple Chaos Synchronization System for Power Quality Classification in a Power System

This document proposes multiple chaos synchronization (CS) systems for power quality (PQ) disturbances classification in a power system. Chen-Lee based CS systems use multiple detectors to track the dynamic errors between the normal signal and the disturbance signal, including power harmonics, voltage fluctuation phenomena, and voltage interruptions. Multiple detectors are used to monitor the dynamic errors between the master system and the slave system and are used to construct the feature patterns from time-domain signals. The maximum likelihood method (MLM), as a classifier, performs a comparison of the patterns of the features in the database. The proposed method can adapt itself without the need for adjustment of parameters or iterative computation. For a sample power system, the test results showed accurate discrimination, good robustness, and faster processing time for the detection of PQ disturbances.


Introduction
The increasing use of electronic equipment can cause electromagnetic disturbances and degrade voltage and current quality. Signal processing, in the time-domain, may have timevarying characteristics [1,2]. For signal analysis applications, Fast Fourier transformation (FFT), short-time Fourier transformation (STFT), and wavelet transformation (WT) based methods [3][4][5][6] have provided promising results in PQ studies. However, the size chosen for the time-window affects both the frequency and time resolution when using a frequencydomain method. Because of the short length of the window at high frequencies and the long length of the window at low frequencies, the frequency-domain method is not suitable for nonstationary signals. The WT, time-frequency domain method, uses cascaded low-pass or high-pass filters and down/upsampling operations to generalize several frequency bands. Specific features render it suitable for the classification of different patterns at specific coefficients in a trial procedure with wavelet decomposition.
Frequency analysis can verify the spectrum characteristics for the classification of PQ disturbances, but the transformation process needs time and requires much memory; so for, online applications, the number of samples is limited.
To identify a broad range of features in the time-domain, a dynamic chaos system is a nonlinear deterministic system that constructs mathematical models for systems that are characterized by random signals. This system has some remarkable characteristics, such as heightened sensitivity to initial conditions, broad spectrums of Fourier transforms, and the fractal properties of motion in phase space [7,8]. It has been extensively studied in many engineering applications, such as adaptive control systems, biomedical signal processing, fluid mechanics, secure communication, and information processing [9][10][11][12][13][14]. The system consists of a master system (MS) and a slave system (SS). The trajectory of a SS can automatically track the trajectory of a MS. Let MS be a reference, so the dynamic errors are defined as the differences between the SS and the MS and are used to construct various patterns for the classification of 2 The Scientific World Journal PQ disturbances. Therefore, multiple Chen-Lee systems are composed of multiple CS detectors that extract the feature patterns from periodic and nonperiodic signals.
Pattern recognition is important in a detection system, so the artificial intelligence method is designed to classify different types of disturbances, such as multilayer neural networks, wavelet neural networks (WNNs), and support vector machines (SVMs) [4-6, 15, 16]. A classifier uses the maximum likelihood method (MLM) [17][18][19] to distinguish between different patterns found in the database, with reference to the threshold of rejection. It performs a comparison of the templates of stored patterns and selects the maximum matching likelihood of the pattern above the threshold. It is simpler and less sensitive to the distortion of patterns than traditional methods. It also has a flexible pattern mechanism, with changing training patterns in a dynamic modeling system, to further heighten the accuracy.
For a 14-bus system, the tests demonstrated computational efficiency and accurate recognition. The organization of this paper is as follows: Section 2 addresses the mathematical methods including multiple CS systems and maximum likelihood method and Section 3 describes the implementation of CS system. Sections 4, 5, and 6 present the construction of the multiple CS detectors, experimental results, and conclusions, which show the efficiency of the proposed model.

Multiple Chaos Synchronization (CS) System
In PQ studies, the disturbances encountered in the process may be time-varying characteristics and voltage fluctuation phenomena, occurring in transient, short-duration, or long-duration variations, involving sag, swell, and harmonic events. The distortion level is investigated and characterized using long-duration recording and statistical events. In many practical applications, FFT and STFT have been applied to analyze transient or short-time fluctuation phenomena. However, when the operating states are partially or totally unavailable, the length chosen for the time-window affects both the frequency and time resolution, so the application of state tracking is restrictive. Therefore, discrete-time, chaotic system is proposed for the tracking of dynamic behaviors, using a computer. With the definition of the master system as a reference, a slave input was designed to synchronously track the dynamic errors, with respect to the reference model. The Chen-Lee system, with known and unknown parameters, shows good symmetrical behavior and has been used in CS applications [20][21][22]. Generally, it consists of a master system (MS) and a slave system (SS) and can be described as Then, subtract (2) from (1), and the dynamic-error system can be expressed by According to previous studies [20,23], the system (4) acts as a chaotic attractor that must satisfy the system parameters , , and , for the following specific condition: > 0, < 0, < 0, 0 < < − ( + ) .
The scaling parameters , , and are nonzero constants that control the chaotic motion of the system (4), and phase trajectories indicate that chaotic motion can be controlled. The MS and SS are synchronization such that 1 = 0, 2 = 0, and 3 = 0. For time series analysis, this technique can be used to examine various dynamic errors from the bilateral signals or patterns, such as periodic signals or aperiodic signals in the time-domain.
For computational analysis, the discrete CS systems are used to describe the dynamic rules and time is represented by the integer . Let the error variables be 1  and ; let be a sequence of sampling data obtained from known signal or training patterns; let be a sequence of sampling data obtained from unknown signal or testing patterns; and is the number of sampling points in one periodic signal. So the general form of dynamic-error system is defined as ] ] ] .
The Scientific World Journal 3 A Chen-Lee based CS detector is used to show that the tracking dynamic errors are zero or have bounded ranges. The dynamic errors Φ 1 , Φ 2 , and Φ 3 can be calculated using system (6) with system parameters , , and , subject to the constricted condition of (5). System (6) can be used to track the dynamic errors and identify the relevant features of the unknown signals.
The multiple CS systems consist of two stages, as shown in Figure 1. In first stage, it specifies the fundamental signal and unknown signal , cycle by cycle. Then, the dynamicerror equation Φ 2 1 ∈ real number can be used to track the dynamic errors between the starting point, = 1, and the ending point, = − 2, and then testing pattern Φ 2 1 connects to each SS. In second stage, we can systematically create training patterns Φ 2 ( ) as where = 1, 2, 3, . . . , and is the number of PQ disturbances, including harmonics (har), voltage sags (sa), voltage swells (sw), sags or swells involving harmonics (sah and swh), or momentary interruptions (int). With reference to (7), the input of each MS is derived from the training patterns , as shown in Figure 1.
For each disturbance, we can store training patterns in the database. Input testing pattern Φ 2 1 , error patterns Φ 2 can be obtained using dynamic-error equation (6). In this study, we choose system parameters = 2, = −3, and = −2 to produce a stable CS attractor. For classes of classification, the likelihood function , = 1, 2, 3, . . . , , is defined as [18,19] The covariance 2 for each class is defined as where the number of classes is , the number of training patterns in class is , and is the weighted value for each class. More training patterns of a class can be continuously stored in the database. For classes classification, the dimension of the vector can be reduced from -dimension to -dimension ( > ). The likelihood function of the same class will be slightly added in a classifier. Index gives the likelihood that a testing pattern Φ 2 1 belongs to training pattern Φ 2 by its value and classification results. The final classification is obtained by choosing the class with the highest probability estimate as * = arg max For all of the total likelihoods relating to a hypothesis, the maximum likelihood * is identified as the correct likelihood.

Implementation of CS System
A dynamic-error system (4) can be expressed using the integrators and four-fundamental arithmetical operations, including adders, subtractors, and multipliers. We use the dynamic-error equations to implement the CS system, using nonlinear electronic circuits [21,24], as shown in Figure 2. Using the analog electronic analysis, it can simulate the chaotic Chen-Lee system as an electronic oscillator circuit. The nonlinear circuit uses operational amplifiers (OPAs), resistors, and capacitors. The system parameters (amplifier gains) can be adjusted by  Figure 1: Architecture of multiple chaos synchronization (CS) systems.

Multiplier
We can estimate the values of resistors 5 , 7 , 15 , 17 , 23 , and 24 by using (14) and satisfying specific condition (5). A nonlinear electronic circuit is a chaotic attractor, so the default initial conditions are zero. The phase diagrams can be plotted for the time interval 1/60 seconds. The time scale is where is the number of sampling points.  For the system parameters = 2, = −3, and = −2 and the initial conditions as 1 = 2 = 3 = 0, and 1 = 2 = 3 = 0, the phase diagrams of the simulation experimental results are shown in Figure 3. One of the two periodic signals is the ideal voltage signal with 60 Hz-fundamental frequency and the other is the voltage sag signal, as shown in Figure 3(a). Voltage sags are 70% to 90% reductions in the rated voltage with short-duration variations. The sampling rate is = 1.44 kHz, and the number of sample points is 24 ( = 24, 12 points in the positive half-cycle, and 12 points in the negative half-cycle). The Chen-Lee based CS detector acts to track the dynamic errors between the fundamental and sag signals. The dynamic errors Φ 1 , Φ 2 , and Φ 3 reflect the chaotic behaviors and control the motions in a bounded region with the specific parameters , , and . As can be seen, from the phase diagram in Figure 3(b), the dynamic errors increase as voltage sag increases, and the trajectories of the chaotic motions increases. Each trajectory is bounded within an elliptical region, Φ 1 : ±1.0, Φ 2 : ±2.0 and Φ 3 : ±1.0, as shown in Figure 3(c). The dynamic errors increase as the distortions gradually become pronounced, especially the dynamic error Φ 2 . The same phenomena occur in other PQ disturbances. Therefore, the dynamic error Φ 2 was used to construct various patterns for the classification of PQ disturbances. (i) Voltage sag (dip) is a sudden voltage drop of 10%∼90% in magnitude. It often lasts for one-half cycle to 1 minute. When the voltage drops by 30% or more, we consider the system status to be severe, such as a fault in a power system or switching of heavy loads.

Multiple CS Detectors Construction
(ii) Voltage swell is a voltage rise above 110% of the normal voltage for the duration of between one-half cycle and 1 minute. It has been known to cause a single lineto-ground fault on the system, a failure of marginal components in electronic equipment, and large loads or capacitors are removed.
(iii) Interruptions can be classified as momentary or sustained voltage interruptions. An interruption occurs when the supply voltage drops to less than 90% of the normal voltage for the duration of one-half cycle to 5 minutes. These events can be the result of faults, equipment failures, or control malfunctions.
(iv) Harmonics are fundamental voltages that have integer multiples of the frequency. Total harmonic distortion (THD) levels exceed the standard [1,2], due to the nonlinear characteristics of power electronic equipment and loads on the power system. A 14-bus system was used for the test example, as shown in Figure 4. The system has 5 generator buses, 15 lines, and 5 transformers. There are 8 harmonic sources in Zone 1 (Buses 6, 11, 12, and 13) and Zone 2 (Buses 4, 7, 9, and 10), including power rectifiers or converters. Table 1 Figure 4: One-line diagram of the 14-bus power system. A sampling rate, = 1.44kHz, was used and the number of sampling points was 24 ( = 24, 12 points in the positive half-cycle, and 12 points in the negative halfcycle). According to (6) with = 2, = −3, and = −2, the CS system is a chaotic attractor. Dynamic errors Φ 1 , Φ 2 , and Φ 3 can be estimated using normal voltages and distorted voltages. Their motion trajectories are in a bounded region within ±2, and show good symmetrical behavior, as shown in Figure 5. These elliptical trajectories follow the socalled butterfly patterns. The motion trajectories of normal voltages describe small ellipses and increase as the variation in voltage magnitudes gradually becomes more pronounced. In addition, the patterns reveal saw tooth and polygonal characteristics for voltage sag or voltage swell involving harmonics.  The scatter diagram of magnitude-duration shows that for Taiwan's science parks and high-tech customers, most magnitudes were distributed between 0.70∼0.90 per unit (p.u.) and 1.10∼1.30 p.u. within 10∼30 cycles [2]. Therefore, this study considers the short-and long-duration voltage fluctuation phenomena. The 61 sets of training patterns for multiple CS detectors at each observation bus, include the following events. (v) voltage interruption (int): voltage magnitude less than 10% of nominal with 3 sets of training patterns, = 3.
A database was constructed for the 61 training patterns ( = 61 in this study) to validate the proposed method. These patterns can be classified into 7 categories. We systematically created template patterns Φ( ), = 1, 2, 3, . . . , 61, such that we had 61 sets of comparative patterns for the MLM based classifier. The final probability (Φ/Φ( )) is between 0 and 1, found using (8). Equation (10) was used to find the highest probability. For 7 categories classification, a threshold of rejection was used to confirm the decision. This can be defined as [25] The index * = arg max { ℎ , , ℎ , , ℎ , , } is maximum one and is up to the threshold value judge . It is as high as possible and approaches to 1. The threshold judge as (17) provides a degree of confidence for the classification of PQ disturbances.

Experimental Results
The proposed CS based detector was designed using a PC Pentium-IV (2.4 GHz with 480 MB RAM) and MATLAB software. Studies of harmonic and voltage fluctuation demonstrate the effectiveness of the proposed method.

Voltage Sag and Swell Involving Harmonics.
Suppose that there are multiple harmonic sources at Bus6, Bus9, Bus12, and Bus13, and the total harmonic voltage distortion of each bus is greater than 2.5%. If the voltage sags caused by heavy motor loads are at Bus12, then multiple harmonic sources are at the Bus6, Bus9, and Bus13. Figure 6(a) shows the harmonic voltages and voltage sag within about 30 cycles in the timedomain. When the voltage sag suddenly reaches 86∼88% of the rated voltage, the multiple CS detectors acts to track the dynamic errors between the known pattern Φ 2 ( ) and unknown pattern Φ 2 . As shown in Figure 6(b), the critical dynamic errors at the beginning and ending of events were noted to monitor the voltage sags involving harmonics. As can be seen in Figure 6(c), dynamic errors Φ 1 , Φ 2 , and Φ 3 are estimated using normal voltage and distorted voltages.
The various morphologies can also be used to observe PQ disturbances. For 50 detection cycles, the MLM based classifier performed a comparison of the patterns in the database. If pattern Φ 2 was similar to any patterns Φ 2 ( ), = 1, 2, 3, . . . , 61, the dynamic errors approached zero. Figure 7 shows the maximum matching likelihood between different patterns. The average output was computed using (8) , ℎ , , }, the final outputs indicate the "har" and "sah, " respectively. This confirms that the proposed method has a higher accuracy for PQ disturbances.
When large loads are suddenly removed at Bus 12, voltage swells rise to 116∼118% and 112∼114% of the rated voltage. The scenario shown in Figure 6

Performance Test and Discussion.
In this study, computer simulation was used to change traced disturbances. Voltage data were produced for varying voltage magnitudes and a variation in harmonic components varying from −50% to +50%, THD % ≥ 2.5%. Figure 8(a) shows the average probability versus voltage magnitude variance for 350 untrained data including nor, sa, sah, sw, swh, and int events. With training patterns for sa and sah events having a specific sag range between 0.7 and 0.9 per unit, voltage magnitudes between 0.50 and 0.90 per unit were identified as "sa" and "sah" events. The results can be also observed for "sw" and "swh" events between 1.10 and 1.50 per unit. Voltage interruptions were gradually identified for a voltage magnitude below 0.20 per unit. With the 61 sets of training patterns, the proposed method could work in a dynamic environment, with voltage fluctuation with harmonics. The experimental results show that the proposed method has high confidence in the classification of all events, using the rejection threshold judge .
To compare the proposed method with the traditional method, a multilayer neural network was also designed as a classifier, using a least mean square (LMS) algorithm. However, the choices of the network architectures, initial network parameters, learning rates, and convergent condition affect learning performances [6]. As the number of training patterns increases (high-dimensional pattern space), the training process and classification efficiency become a bottleneck using the gradient methods, steepest descent methods, or optimization methods. The proposed method was used to determine the maximum matching likelihood. Although the performance of the MLM based method is affected by covariance, values are automatically adjusted by the weighted value and training patterns for various voltage fluctuation phenomena, as shown in Figure 8(b). The proposed method can adapt itself in real-time system applications without the need for parameter assignment or iterative computation.

Conclusion
The multiple CS detectors for power quality disturbances classification have been examined in this paper. The Chen-Lee based chaos system was used to track the dynamic errors between the normal signal and the distorted signals in the time-domain. The butterfly patterns reveal various morphological motions for harmonics, voltage fluctuation phenomena, and multiple disturbances. These motions can also point out the beginning and ending of disturbance occurrences. Subsequent, multiple CS detectors act to track dynamic errors between feature pattern and training patterns. The maximum likelihood method (MLM) then used these feature patterns to classify PQ disturbances. It can automatically adjust the parameters and could be suitable for application in a dynamical environment. To ensure the continuity service and higher reliability, the proposed method also provides a promising result for fault detection in intelligent microgrid, including short-circuit and line-to-ground faults. Intelligent electronic devices, such as digital protective relays and digital fault recorders, can be implemented through analog   electronic circuits, microprocessor, and communication links for further use in protection system.