# Power Load Event Detection and Classification Based on Edge Symbol Analysis and Support Vector Machine

Authors: Lei Jiang; Jiaming Li; Suhuai Luo; Sam West; Glenn Platt
Journal: Applied Computational Intelligence and Soft Computing
 (2012)

---

## Abstract

Energy signature analysis of power appliance is the core of nonintrusive load monitoring (NILM) where the detailed data of the appliances used in houses are obtained by analyzing changes in the voltage and current. This paper focuses on developing an automatic power load event detection and appliance classification based on machine learning. In power load event detection, the paper presents a new transient detection algorithm. By turn-on and turn-off transient waveforms analysis, it can accurately detect the edge point when a device is switched on or switched off. The proposed load classification technique can identify different power appliances with improved recognition accuracy and computational speed. The load classification method is composed of two processes including frequency feature analysis and support vector machine. The experimental results indicated that the incorporation of the new edge detection and turn-on and turn-off transient signature analysis into NILM revealed more information than traditional NILM methods. The load classification method has achieved more than ninety percent recognition rate.


---

## Body

## 1. Introduction

Nowadays, multifarious classification techniques are widely used in many signal-processing areas, such as face recognition, road traffic analyze, medical image, weather forecasting, and so forth [1, 2]. It has been considered whether this kind of techniques could serve the purposes of power load monitoring, which is a process for obtaining what appliances are used in the house as well as their individual energy consumption by analyzing changes in the voltage and current. In particular, information about the operating conditions of consumers (such as what appliances are operated and what state they are in) is useful for demand estimation and prediction. In addition, the operating conditions of electrical appliances in modern society clearly reflect consumers’ lifestyles and behavior patterns [3]. Research is being performed on the use of such information for circuit diagnoses, safety confirmation systems for elderly persons living alone [4], demand management, and optimization [5–10], and for other uses [11].There are mainly two classes of approaches in power monitoring algorithm, including intrusive appliance load monitoring (IALM) and nonintrusive appliance load monitoring (NIALM) [12]. In order to measure consumption, IALM distributes direct sensors at each device or appliance. Although conceptually straightforward and potentially highly accurate, direct sensing is often expensive due to time consuming installation and the requirement for one sensor for each device or appliance. In response to limitations with the direct sensing approach, researchers have explored methods to infer energy usage via a single sensor.In contrast to the direct sensing methods, the standard NIALM configuration includes a sole sensor set to measure current and voltage, as well as a processing algorithm for determining the status of various devices [13]. It was firstly studied by Hart of MIT in the early 1990s with funding from the Electric Power Research Institute [12]. Since then, different strategies for NIALM have been developed. One of the important approaches on classification is to apply support vector machine (SVM) [14–18] in the nonintrusive appliance load monitoring. It is imperative for power system research field to evaluate the SVM on this task from a practical point of view. As the preprocess for finding out the features for SVM, transient detection also plays an important role in NIALM.In order to decompose the total loads into their components, Bijker et al. employed different power level, or named step change for the purpose [19]. Chang et al. used DWT and a new method coreless HCT to detect the power events [20]. Nonetheless, they only supposed to detect the events of appliances with complicated features, such as motor, microwave oven, and thyristor rectifier. In 2002, Onoda et al. used the data proposed in [21] for the purpose of estimating the state of household electric appliances. They compared different types of SVMs obtained by choosing different kernels. They reported results of polynomial kernels, radial basis function kernels, and sigmoid kernels. All results for the three different kernels achieved almost same error rates. However, in the estimation of the state of household electric appliances, the results for the three different kernels achieved different error rates. They also compared different capacity of SVMs obtained by choosing different regularization constants and parameters of kernels experimentally. The results showed that the capacity control is as important as the choice of kernel functions. Kadouche and his colleagues later presented [22] their ongoing work on the house occupant prediction issue based on daily life habits in smart houses. Most of their works were based on supervised learning technical. They used SVM to build behavior classification model for learning the user’s habits, analyzed the publicly available dataset from the Washington State University Smart Apartment Test-bed. Particularly, they evaluated the grooming, having breakfast, and bed to toilet activities [23]. Their experimental results showed that the user can be recognized with a high precision which means that each user has his own way to perform activities. As future work, the users’ patterns which allow a person to be discriminated and recognized among a group, performing multiactivities in the same environment without using intrusive technologies were being studied.Grinblat and his fellows presented a new method for generating adaptive classifiers and capable of learning concepts that change with time, which is the time-adaptive support vector machine (TA-SVM) [24]. The basic idea of TA-SVM is to use a sequence of classifiers, each one is appropriate for a small time window but, in contrast to other proposals, learning all the hyperplanes in a global way. Starting from the solution of independent SVMs, they showed that the addition of a new term in the cost function (which penalizes the diversity between consecutive classifiers) produces in fact a coupling of the sequence. Once coupled, the set of SVMs acts as a single adaptive classifier. They evaluated different aspects of the TA-SVM using artificial drifting problems. In particular, they showed that changing the number of classifiers and the coupling constant can effectively regularize the sequence of classifiers. They compared TA-SVM with other state-of-the-art methods in three different settings: estimation, prediction, and extrapolation, including problems with small datasets, high-dimensional input spaces, and noise. TA-SVM showed in all cases to be equivalent to or better than the other methods. Even for the most unfavorable situation for TA-SVM, that is, the sudden changes of the dataset, their new method showed a very good performance. However, the limitation of this method may be reflected in the long drifting time and the requirement of additional hardware. A more recent research of Liang’s group [25], combined various features including current waveform (CW), active/reactive power (PQ), harmonics (HAR), instantaneous admittance waveform (IAW), instantaneous power waveform (IPW), eigenvalues (EIG), and switching transient waveform (STW). The results of their research provided a higher degree of recognition precision, but the algorithm requires mountains of work on collecting and processing appliance signatures.In this paper, the data collected in real world is used and clearly analyzed what are the important issues for applying SVM to power system research field. The remainder is organized as below. Section2 gives the details of the new transient event detection process named ESA and the reasons why we choose SVM as our method and the main process of load classification. The details of data collecting and preprocessing along with the experiment results are shown in Section 3. Finally, conclusions are made in Section 4.


## 2. Power Load Events Detection and Classification

It is known that although NIALM based on different techniques, it has several common principles [26].(i)
Load features classification: specific appliance features, or signatures, need to be selected and mathematically characterized.(ii)
Mechanism: a hardware installation (sensor and data acquisition system) that can detect the selected features is required.(iii)
Decomposition: a mathematical algorithm detects the features in the overall signal and output.Based on this theory, Figure1 describes the process of this research. The mechanism collects the whole operation of a circuit, then, the appropriate algorithms draw out all the electrical events in the circumstance. After successfully seeking the point “load x switch on,” it is feasible to classify load x, afterwards seek the point “load x+1 switch on” excluding faked edges generated by previous loads. The same procedure subsequently goes for searching switch off events.Figure 1
The load events detection procedure.


### 2.1. Power Load Events Detection Using Edge Symbol Analysis Method

#### 2.1.1. Background

Methods of edge detection have been widely studied in many areas, primarily developed in the discrete cases and have been confined to slightly noisy signals. The regularization or smoothing and optimal approaches of Canny [27] have led to several efficient continuous operators for noisy and blurred signals. Other advanced methods that consider the Canny criterion have been developed to deal with noise, uneven illumination, and image contrast. In a marginal way, discrete approaches for regularization have been developed and have improved results by considering the discrete nature of the signals [28, 29].Introducing nonlinearity into the global filtering process, as in noisy edge detection, is a marginal yet efficient method of obtaining good performance. Generally, nonlinear filtering is used in a preliminary regularization stage. Pitas and Venetsanopoulos proposed a class of nonlinear filters that reject additive and impulse noise, while preserving the edges [30]. More recently, Benazza-Benyahia et al. introduced a nonlinear filter bank leading to a multiresolution approximation of the signal in which discontinuities are well preserved [31]. A nonlinear filter for edge enhancement, using a morphological filter, has been proposed by Schulze [32]. The author showed that local variation analysis allows enhancing edges corrupted by multiplicative noise. Hwang and Haddad presented an integrated nonlinear edge-detection-based denoising scheme [33]. A thresholded derivative is computed from two half filters (median for impulse noise, mean for Gaussian noise, and min-max for uniform noise) and edge detection is used to select the second filtering stage, that is, mean for noise or median for edge points. In the scheme, edge detection could be considered as a by-product and the optimal performance is obtained only when the correct first filter is selected according to the noise statistic. Based on these theories, Laligant and his colleague propose to obtain both noise reduction and edge detection by a one-stage nonlinear derivative scheme [34]. The scheme, which consists of combining two polarized differences, yields significant improvements in signal-to-noise ratio without using regularization or increasing the computational requirements.


#### 2.1.2. The Approach

It is known that three kinds of slopes can be found in all power waveforms, including positive, negative, and constant. So the edge symbol analysis method (ESA) can find the position of the edge point according to the sign of the signal area. In this way, the edge point will be located after the transition if the slope is positive. Otherwise, when it is negative, the edge point will be located before the slope. Here, for the sake of detecting these two kinds of edges, two detector filtersL+ and L- have to be introduced. Since L is chosen as an antisymmetric linear filter, which gives the localization changes and this shifted pixel localization depending on edge’s orientation. These filters have the relationship described in (1).
(1)L-(S)=-L+(S-1).
With the simple threshold H, we use the combination of L+ and L- as below to detect edge points of the original input signal I. 										(2)O+(S)=H[L+(S)I(S)],O-(S)=-H[-L-(S)I(S)].


### 2.2. SVM-Based Power Load Classification

The proposed power classification method is composed of two major stages, including: first a frequency feature analysis is applied on current signal, and then a trained classifier based on SVM is applied to identify different appliances. The following sections give the details of the approach.


#### 2.2.1. Current Feature Analysis

To save computational resources and improve performance, only current signals were used for frequency analysis with fewer sample points [35]. For example, according to previous experiments, the accuracy rate of the classification is almost the same or even higher when employing 500 frequency signal points as features instead of 2000 unprocessed current signal points.After collecting the data of single device and extracting the current signals, a short-time fast Fourier transform (FFT) of the signal is performed as in (3). Where ωN=e(2πi)/N is an nth root of the unity. 										(3)X(k)=∑j=1Nx(j)ωN(j-1)(k-1),x(j)=(1N)∑k=1NX(k)ωN-(j-1)(k-1),
and an example is shown as Figure 2.Figure 2
Load current signal of LCD monitor in time and frequency domains.


#### 2.2.2. SVM-Based Classification

SVM, as one of the method of NIALM system, is relatively insensitive to the number of data points and the classification complexity does not depend on the dimensionality of the feature space [36]. Therefore, it can potentially learn a larger set of patterns and be able to scale better than some other methods. Once the data are classified into two classes, a suitable optimizing algorithm can be used if necessary for further feature identification, depending on the application.Support vector machine is a training algorithm for learning classification and regression rules from data SVM was first introduced by Cortes and Vapnik [37] in the 1990s for classification and have recently become an area of intense research owing to developments in the techniques and theory coupled with extensions to regression and density estimation. It is based on the structural risk minimization principle which incorporates capacity control to prevent overfitting. It is a partial solution to the bias-variance trade-off dilemma [38]. It has been widely used in different areas. As a classifier, it gives a set of training examples, each marked as belonging to one of two categories. The SVM training algorithm builds a model that predicts whether a new example falls into one category or the other. Intuitively, an SVM model is a representation of the examples as points in space, mapped so that the examples of the separate categories are divided by a clear gap that is as wide as possible [39]. New examples are then mapped into that same space and predicted to belong to a category based on which side of the gap they fall on.SVM delivers state-of-the-art performance in real-world applications such as text categorization, hand-written character recognition, biosequences analysis, image classification, and so forth. It is now established as one of the standard tools for machine learning and data mining. The SVM decision function is defined as (4)
(4)f(y)=∑i=1NαiK(xi,y)+b.
Here, y is the unclassified tested vector, xi is the support vectors and αi is their weights, and b is a constant bias. K(x,y) is the kernel function introduced into SVM to solve the nonlinear problems by performing implicit mapping into a high-dimensional feature space [40].Consider the problem of separating the set of training vectors belonging to two separate classes [41]. 										(5)D={(x1,y1),…,(xl,yl)},x∈Rn,y∈{-1,1}
with a hyperplane,
(6)〈ω,x〉+b=0.
The kernel function chosen in our algorithm is Gaussian
(7)k(xi,xj)=exp(-γ∥xi-xj∥2),forγ>0,
where scalar γ is identical for all coordinates.


## 2.1. Power Load Events Detection Using Edge Symbol Analysis Method

### 2.1.1. Background

Methods of edge detection have been widely studied in many areas, primarily developed in the discrete cases and have been confined to slightly noisy signals. The regularization or smoothing and optimal approaches of Canny [27] have led to several efficient continuous operators for noisy and blurred signals. Other advanced methods that consider the Canny criterion have been developed to deal with noise, uneven illumination, and image contrast. In a marginal way, discrete approaches for regularization have been developed and have improved results by considering the discrete nature of the signals [28, 29].Introducing nonlinearity into the global filtering process, as in noisy edge detection, is a marginal yet efficient method of obtaining good performance. Generally, nonlinear filtering is used in a preliminary regularization stage. Pitas and Venetsanopoulos proposed a class of nonlinear filters that reject additive and impulse noise, while preserving the edges [30]. More recently, Benazza-Benyahia et al. introduced a nonlinear filter bank leading to a multiresolution approximation of the signal in which discontinuities are well preserved [31]. A nonlinear filter for edge enhancement, using a morphological filter, has been proposed by Schulze [32]. The author showed that local variation analysis allows enhancing edges corrupted by multiplicative noise. Hwang and Haddad presented an integrated nonlinear edge-detection-based denoising scheme [33]. A thresholded derivative is computed from two half filters (median for impulse noise, mean for Gaussian noise, and min-max for uniform noise) and edge detection is used to select the second filtering stage, that is, mean for noise or median for edge points. In the scheme, edge detection could be considered as a by-product and the optimal performance is obtained only when the correct first filter is selected according to the noise statistic. Based on these theories, Laligant and his colleague propose to obtain both noise reduction and edge detection by a one-stage nonlinear derivative scheme [34]. The scheme, which consists of combining two polarized differences, yields significant improvements in signal-to-noise ratio without using regularization or increasing the computational requirements.


### 2.1.2. The Approach

It is known that three kinds of slopes can be found in all power waveforms, including positive, negative, and constant. So the edge symbol analysis method (ESA) can find the position of the edge point according to the sign of the signal area. In this way, the edge point will be located after the transition if the slope is positive. Otherwise, when it is negative, the edge point will be located before the slope. Here, for the sake of detecting these two kinds of edges, two detector filtersL+ and L- have to be introduced. Since L is chosen as an antisymmetric linear filter, which gives the localization changes and this shifted pixel localization depending on edge’s orientation. These filters have the relationship described in (1).
(1)L-(S)=-L+(S-1).
With the simple threshold H, we use the combination of L+ and L- as below to detect edge points of the original input signal I. 										(2)O+(S)=H[L+(S)I(S)],O-(S)=-H[-L-(S)I(S)].


## 2.1.1. Background

Methods of edge detection have been widely studied in many areas, primarily developed in the discrete cases and have been confined to slightly noisy signals. The regularization or smoothing and optimal approaches of Canny [27] have led to several efficient continuous operators for noisy and blurred signals. Other advanced methods that consider the Canny criterion have been developed to deal with noise, uneven illumination, and image contrast. In a marginal way, discrete approaches for regularization have been developed and have improved results by considering the discrete nature of the signals [28, 29].Introducing nonlinearity into the global filtering process, as in noisy edge detection, is a marginal yet efficient method of obtaining good performance. Generally, nonlinear filtering is used in a preliminary regularization stage. Pitas and Venetsanopoulos proposed a class of nonlinear filters that reject additive and impulse noise, while preserving the edges [30]. More recently, Benazza-Benyahia et al. introduced a nonlinear filter bank leading to a multiresolution approximation of the signal in which discontinuities are well preserved [31]. A nonlinear filter for edge enhancement, using a morphological filter, has been proposed by Schulze [32]. The author showed that local variation analysis allows enhancing edges corrupted by multiplicative noise. Hwang and Haddad presented an integrated nonlinear edge-detection-based denoising scheme [33]. A thresholded derivative is computed from two half filters (median for impulse noise, mean for Gaussian noise, and min-max for uniform noise) and edge detection is used to select the second filtering stage, that is, mean for noise or median for edge points. In the scheme, edge detection could be considered as a by-product and the optimal performance is obtained only when the correct first filter is selected according to the noise statistic. Based on these theories, Laligant and his colleague propose to obtain both noise reduction and edge detection by a one-stage nonlinear derivative scheme [34]. The scheme, which consists of combining two polarized differences, yields significant improvements in signal-to-noise ratio without using regularization or increasing the computational requirements.


## 2.1.2. The Approach

It is known that three kinds of slopes can be found in all power waveforms, including positive, negative, and constant. So the edge symbol analysis method (ESA) can find the position of the edge point according to the sign of the signal area. In this way, the edge point will be located after the transition if the slope is positive. Otherwise, when it is negative, the edge point will be located before the slope. Here, for the sake of detecting these two kinds of edges, two detector filtersL+ and L- have to be introduced. Since L is chosen as an antisymmetric linear filter, which gives the localization changes and this shifted pixel localization depending on edge’s orientation. These filters have the relationship described in (1).
(1)L-(S)=-L+(S-1).
With the simple threshold H, we use the combination of L+ and L- as below to detect edge points of the original input signal I. 										(2)O+(S)=H[L+(S)I(S)],O-(S)=-H[-L-(S)I(S)].


## 2.2. SVM-Based Power Load Classification

The proposed power classification method is composed of two major stages, including: first a frequency feature analysis is applied on current signal, and then a trained classifier based on SVM is applied to identify different appliances. The following sections give the details of the approach.


### 2.2.1. Current Feature Analysis

To save computational resources and improve performance, only current signals were used for frequency analysis with fewer sample points [35]. For example, according to previous experiments, the accuracy rate of the classification is almost the same or even higher when employing 500 frequency signal points as features instead of 2000 unprocessed current signal points.After collecting the data of single device and extracting the current signals, a short-time fast Fourier transform (FFT) of the signal is performed as in (3). Where ωN=e(2πi)/N is an nth root of the unity. 										(3)X(k)=∑j=1Nx(j)ωN(j-1)(k-1),x(j)=(1N)∑k=1NX(k)ωN-(j-1)(k-1),
and an example is shown as Figure 2.Figure 2
Load current signal of LCD monitor in time and frequency domains.


### 2.2.2. SVM-Based Classification

SVM, as one of the method of NIALM system, is relatively insensitive to the number of data points and the classification complexity does not depend on the dimensionality of the feature space [36]. Therefore, it can potentially learn a larger set of patterns and be able to scale better than some other methods. Once the data are classified into two classes, a suitable optimizing algorithm can be used if necessary for further feature identification, depending on the application.Support vector machine is a training algorithm for learning classification and regression rules from data SVM was first introduced by Cortes and Vapnik [37] in the 1990s for classification and have recently become an area of intense research owing to developments in the techniques and theory coupled with extensions to regression and density estimation. It is based on the structural risk minimization principle which incorporates capacity control to prevent overfitting. It is a partial solution to the bias-variance trade-off dilemma [38]. It has been widely used in different areas. As a classifier, it gives a set of training examples, each marked as belonging to one of two categories. The SVM training algorithm builds a model that predicts whether a new example falls into one category or the other. Intuitively, an SVM model is a representation of the examples as points in space, mapped so that the examples of the separate categories are divided by a clear gap that is as wide as possible [39]. New examples are then mapped into that same space and predicted to belong to a category based on which side of the gap they fall on.SVM delivers state-of-the-art performance in real-world applications such as text categorization, hand-written character recognition, biosequences analysis, image classification, and so forth. It is now established as one of the standard tools for machine learning and data mining. The SVM decision function is defined as (4)
(4)f(y)=∑i=1NαiK(xi,y)+b.
Here, y is the unclassified tested vector, xi is the support vectors and αi is their weights, and b is a constant bias. K(x,y) is the kernel function introduced into SVM to solve the nonlinear problems by performing implicit mapping into a high-dimensional feature space [40].Consider the problem of separating the set of training vectors belonging to two separate classes [41]. 										(5)D={(x1,y1),…,(xl,yl)},x∈Rn,y∈{-1,1}
with a hyperplane,
(6)〈ω,x〉+b=0.
The kernel function chosen in our algorithm is Gaussian
(7)k(xi,xj)=exp(-γ∥xi-xj∥2),forγ>0,
where scalar γ is identical for all coordinates.


## 2.2.1. Current Feature Analysis

To save computational resources and improve performance, only current signals were used for frequency analysis with fewer sample points [35]. For example, according to previous experiments, the accuracy rate of the classification is almost the same or even higher when employing 500 frequency signal points as features instead of 2000 unprocessed current signal points.After collecting the data of single device and extracting the current signals, a short-time fast Fourier transform (FFT) of the signal is performed as in (3). Where ωN=e(2πi)/N is an nth root of the unity. 										(3)X(k)=∑j=1Nx(j)ωN(j-1)(k-1),x(j)=(1N)∑k=1NX(k)ωN-(j-1)(k-1),
and an example is shown as Figure 2.Figure 2
Load current signal of LCD monitor in time and frequency domains.


## 2.2.2. SVM-Based Classification

SVM, as one of the method of NIALM system, is relatively insensitive to the number of data points and the classification complexity does not depend on the dimensionality of the feature space [36]. Therefore, it can potentially learn a larger set of patterns and be able to scale better than some other methods. Once the data are classified into two classes, a suitable optimizing algorithm can be used if necessary for further feature identification, depending on the application.Support vector machine is a training algorithm for learning classification and regression rules from data SVM was first introduced by Cortes and Vapnik [37] in the 1990s for classification and have recently become an area of intense research owing to developments in the techniques and theory coupled with extensions to regression and density estimation. It is based on the structural risk minimization principle which incorporates capacity control to prevent overfitting. It is a partial solution to the bias-variance trade-off dilemma [38]. It has been widely used in different areas. As a classifier, it gives a set of training examples, each marked as belonging to one of two categories. The SVM training algorithm builds a model that predicts whether a new example falls into one category or the other. Intuitively, an SVM model is a representation of the examples as points in space, mapped so that the examples of the separate categories are divided by a clear gap that is as wide as possible [39]. New examples are then mapped into that same space and predicted to belong to a category based on which side of the gap they fall on.SVM delivers state-of-the-art performance in real-world applications such as text categorization, hand-written character recognition, biosequences analysis, image classification, and so forth. It is now established as one of the standard tools for machine learning and data mining. The SVM decision function is defined as (4)
(4)f(y)=∑i=1NαiK(xi,y)+b.
Here, y is the unclassified tested vector, xi is the support vectors and αi is their weights, and b is a constant bias. K(x,y) is the kernel function introduced into SVM to solve the nonlinear problems by performing implicit mapping into a high-dimensional feature space [40].Consider the problem of separating the set of training vectors belonging to two separate classes [41]. 										(5)D={(x1,y1),…,(xl,yl)},x∈Rn,y∈{-1,1}
with a hyperplane,
(6)〈ω,x〉+b=0.
The kernel function chosen in our algorithm is Gaussian
(7)k(xi,xj)=exp(-γ∥xi-xj∥2),forγ>0,
where scalar γ is identical for all coordinates.


## 3. Experiments

### 3.1. Data Collecting and Preprocessing

All the data are collected in Commonwealth Scientific and Industrial Research Organization (CSIRO) and a residential house in Newcastle, NSW, Australia. Collecting details are shown in Table1.Table 1
Data collecting details.
Load
Description
Dremel
High speed electric motor (a Dremel is a handheld tool for grinding, engraving, etc.). Has built-in speed control which modifies the current sine wave in all except full speed settings.
Fan heater
Small fan heater with fan, medium heat, and high heat modes selected using a slide switch
Heat gun
2-fanspeed digital temperature controlled heat gun
Kettle
A normal kettle
Laptop charger
Dell laptop charger for a D630 laptop
LCD monitor
A Dell LCD monitor
Incandescent light
Incandescent light bulb in small desk lamp
Fluorescent light
Circular fluorescent tube light from soldering station magnifying light stand
Mini PC
Dream PC (labeled PGWY1)
PC
A normal personal computer
Pedestal fan
A cheap $15 3-speed plastic pedestal fanFigure3 shows the segmentation process in this experiment. Each signal of load has three phrases: on phrase, steady state, and off phrase.Figure 3
Process of signal segmentation.


### 3.2. Experiment Result and Discussion

#### 3.2.1. Implementation of ESA

In Figure4, the overlapped data of Kettle and Toaster are shown as an example. Figure 4(a) is the original overlapped data of kettle and oven. The red figures in Figure 4(b) are the output of edge detection and the original data. The on and off edges are represented in red color. To clearly observe the resultant edge points, the edge points in Figure 4(b) (i.e., the red curve) are redrawn in Figure 4(c).The original data and edge analysis of Kettle + Oven.
(a)
(b)
(c)The example of overlapping kettle with microwave oven (Figure5(a)) and PC (Figure 5(b)) is shown in Figure 5. It can easily find the edges (as the yellow arrows shown in Figure 5) by repeating the process that mentioned above, and then relocate the transient points in the original data. When the transient features are recognized as “microwave oven on” or “PC on” by SVM, the algorithm ignores the next “Off” and “On” events because they are fake symbols. It is thus able to obtain the full process of the whole event.Examples of power signal edge analysis.
(a)
Kettle + Microwave oven
(b)
Kettle + PCData of 2 loads combinations from four typical appliances are employed for testing in this experiment. The algorithm runs on 10 different data for each group. As for a single case, the performance time consumption is lower than 4 seconds. The result of this process is shown in Table2. Practically, all events are accurately detected by ESA except the Load2-On, whereas the accuracy is still higher than 80%.Table 2
Result of ESA on four typical loads.
Overlapped load appliances groups
Accuracy of recognizing transient events (switch on and switch off) %
Load 1 On
Load 2 On
Load 2 Off
Load 1 Off
Kettle + Oven
100
100
100
100
Oven + Kettle
100
100
100
100
Kettle + Microwave oven
100
80
100
100
Microwave oven + Kettle
100
90
100
100
Kettle + PC
100
100
100
100
PC + Kettle
100
90
100
100
Oven + Microwave oven
100
100
100
100
Microwave oven + Oven
100
80
100
100
Oven + PC
100
100
100
100
PC + Oven
100
80
100
100


#### 3.2.2. Implementation of Load Classification

The recognition process on each phrase is running with Gaussian as the function kernel. Using the segmented steady-state power data for the test, some selected features are shown in Figures6 and 7, where Figure 6 shows the training unprocessed current data in time domain, and Figure 7 shows the training current data in frequency domain. The algorithm runs 70 times on the data collected from the 11 different loads. As for a single case, the performance time consumption is lower than 2 seconds. The recognition results are shown in Table 3.Table 3
Results of 11 loads classification.
Load appliance
Accuracy of features (500-currence)
Unprocessed current data
Current harmonic data
Dremel
100%
100%
Fan heater
97%
100%
Heater gun
97%
100%
Kettle
90%
97%
Laptop
100%
100%
LCD
100%
100%
Incandescent light
97%
97%
Fluorescent light
100%
100%
Mini PC
100%
100%
PC
100%
100%
Pedestal fan
100%
100%Figure 6
Examples of training set of unprocessed current data in time domain (steady state).Figure 7
Examples of training set of current data in frequency domain (steady state).It can be seen that all the recognition accuracy rates are higher than 90%. It shows the ideal recognition performance of our algorithm. Moreover, despite the training processes always take several minutes, it avoid a mess of pre-work on processing signatures. The testing processes are much quicker than previous methods because of employing comparatively a small quantity of points from load frequency features. However, we would like to indicate that the number of appliance is limited and no transient event is involved. The performance may come down when more devices and events are used in experiments.


## 3.1. Data Collecting and Preprocessing

All the data are collected in Commonwealth Scientific and Industrial Research Organization (CSIRO) and a residential house in Newcastle, NSW, Australia. Collecting details are shown in Table1.Table 1
Data collecting details.
Load
Description
Dremel
High speed electric motor (a Dremel is a handheld tool for grinding, engraving, etc.). Has built-in speed control which modifies the current sine wave in all except full speed settings.
Fan heater
Small fan heater with fan, medium heat, and high heat modes selected using a slide switch
Heat gun
2-fanspeed digital temperature controlled heat gun
Kettle
A normal kettle
Laptop charger
Dell laptop charger for a D630 laptop
LCD monitor
A Dell LCD monitor
Incandescent light
Incandescent light bulb in small desk lamp
Fluorescent light
Circular fluorescent tube light from soldering station magnifying light stand
Mini PC
Dream PC (labeled PGWY1)
PC
A normal personal computer
Pedestal fan
A cheap $15 3-speed plastic pedestal fanFigure3 shows the segmentation process in this experiment. Each signal of load has three phrases: on phrase, steady state, and off phrase.Figure 3
Process of signal segmentation.


## 3.2. Experiment Result and Discussion

### 3.2.1. Implementation of ESA

In Figure4, the overlapped data of Kettle and Toaster are shown as an example. Figure 4(a) is the original overlapped data of kettle and oven. The red figures in Figure 4(b) are the output of edge detection and the original data. The on and off edges are represented in red color. To clearly observe the resultant edge points, the edge points in Figure 4(b) (i.e., the red curve) are redrawn in Figure 4(c).The original data and edge analysis of Kettle + Oven.
(a)
(b)
(c)The example of overlapping kettle with microwave oven (Figure5(a)) and PC (Figure 5(b)) is shown in Figure 5. It can easily find the edges (as the yellow arrows shown in Figure 5) by repeating the process that mentioned above, and then relocate the transient points in the original data. When the transient features are recognized as “microwave oven on” or “PC on” by SVM, the algorithm ignores the next “Off” and “On” events because they are fake symbols. It is thus able to obtain the full process of the whole event.Examples of power signal edge analysis.
(a)
Kettle + Microwave oven
(b)
Kettle + PCData of 2 loads combinations from four typical appliances are employed for testing in this experiment. The algorithm runs on 10 different data for each group. As for a single case, the performance time consumption is lower than 4 seconds. The result of this process is shown in Table2. Practically, all events are accurately detected by ESA except the Load2-On, whereas the accuracy is still higher than 80%.Table 2
Result of ESA on four typical loads.
Overlapped load appliances groups
Accuracy of recognizing transient events (switch on and switch off) %
Load 1 On
Load 2 On
Load 2 Off
Load 1 Off
Kettle + Oven
100
100
100
100
Oven + Kettle
100
100
100
100
Kettle + Microwave oven
100
80
100
100
Microwave oven + Kettle
100
90
100
100
Kettle + PC
100
100
100
100
PC + Kettle
100
90
100
100
Oven + Microwave oven
100
100
100
100
Microwave oven + Oven
100
80
100
100
Oven + PC
100
100
100
100
PC + Oven
100
80
100
100


### 3.2.2. Implementation of Load Classification

The recognition process on each phrase is running with Gaussian as the function kernel. Using the segmented steady-state power data for the test, some selected features are shown in Figures6 and 7, where Figure 6 shows the training unprocessed current data in time domain, and Figure 7 shows the training current data in frequency domain. The algorithm runs 70 times on the data collected from the 11 different loads. As for a single case, the performance time consumption is lower than 2 seconds. The recognition results are shown in Table 3.Table 3
Results of 11 loads classification.
Load appliance
Accuracy of features (500-currence)
Unprocessed current data
Current harmonic data
Dremel
100%
100%
Fan heater
97%
100%
Heater gun
97%
100%
Kettle
90%
97%
Laptop
100%
100%
LCD
100%
100%
Incandescent light
97%
97%
Fluorescent light
100%
100%
Mini PC
100%
100%
PC
100%
100%
Pedestal fan
100%
100%Figure 6
Examples of training set of unprocessed current data in time domain (steady state).Figure 7
Examples of training set of current data in frequency domain (steady state).It can be seen that all the recognition accuracy rates are higher than 90%. It shows the ideal recognition performance of our algorithm. Moreover, despite the training processes always take several minutes, it avoid a mess of pre-work on processing signatures. The testing processes are much quicker than previous methods because of employing comparatively a small quantity of points from load frequency features. However, we would like to indicate that the number of appliance is limited and no transient event is involved. The performance may come down when more devices and events are used in experiments.


## 3.2.1. Implementation of ESA

In Figure4, the overlapped data of Kettle and Toaster are shown as an example. Figure 4(a) is the original overlapped data of kettle and oven. The red figures in Figure 4(b) are the output of edge detection and the original data. The on and off edges are represented in red color. To clearly observe the resultant edge points, the edge points in Figure 4(b) (i.e., the red curve) are redrawn in Figure 4(c).The original data and edge analysis of Kettle + Oven.
(a)
(b)
(c)The example of overlapping kettle with microwave oven (Figure5(a)) and PC (Figure 5(b)) is shown in Figure 5. It can easily find the edges (as the yellow arrows shown in Figure 5) by repeating the process that mentioned above, and then relocate the transient points in the original data. When the transient features are recognized as “microwave oven on” or “PC on” by SVM, the algorithm ignores the next “Off” and “On” events because they are fake symbols. It is thus able to obtain the full process of the whole event.Examples of power signal edge analysis.
(a)
Kettle + Microwave oven
(b)
Kettle + PCData of 2 loads combinations from four typical appliances are employed for testing in this experiment. The algorithm runs on 10 different data for each group. As for a single case, the performance time consumption is lower than 4 seconds. The result of this process is shown in Table2. Practically, all events are accurately detected by ESA except the Load2-On, whereas the accuracy is still higher than 80%.Table 2
Result of ESA on four typical loads.
Overlapped load appliances groups
Accuracy of recognizing transient events (switch on and switch off) %
Load 1 On
Load 2 On
Load 2 Off
Load 1 Off
Kettle + Oven
100
100
100
100
Oven + Kettle
100
100
100
100
Kettle + Microwave oven
100
80
100
100
Microwave oven + Kettle
100
90
100
100
Kettle + PC
100
100
100
100
PC + Kettle
100
90
100
100
Oven + Microwave oven
100
100
100
100
Microwave oven + Oven
100
80
100
100
Oven + PC
100
100
100
100
PC + Oven
100
80
100
100


## 3.2.2. Implementation of Load Classification

The recognition process on each phrase is running with Gaussian as the function kernel. Using the segmented steady-state power data for the test, some selected features are shown in Figures6 and 7, where Figure 6 shows the training unprocessed current data in time domain, and Figure 7 shows the training current data in frequency domain. The algorithm runs 70 times on the data collected from the 11 different loads. As for a single case, the performance time consumption is lower than 2 seconds. The recognition results are shown in Table 3.Table 3
Results of 11 loads classification.
Load appliance
Accuracy of features (500-currence)
Unprocessed current data
Current harmonic data
Dremel
100%
100%
Fan heater
97%
100%
Heater gun
97%
100%
Kettle
90%
97%
Laptop
100%
100%
LCD
100%
100%
Incandescent light
97%
97%
Fluorescent light
100%
100%
Mini PC
100%
100%
PC
100%
100%
Pedestal fan
100%
100%Figure 6
Examples of training set of unprocessed current data in time domain (steady state).Figure 7
Examples of training set of current data in frequency domain (steady state).It can be seen that all the recognition accuracy rates are higher than 90%. It shows the ideal recognition performance of our algorithm. Moreover, despite the training processes always take several minutes, it avoid a mess of pre-work on processing signatures. The testing processes are much quicker than previous methods because of employing comparatively a small quantity of points from load frequency features. However, we would like to indicate that the number of appliance is limited and no transient event is involved. The performance may come down when more devices and events are used in experiments.


## 4. Conclusion

This paper has proposed one approach of nonintrusive appliance load monitoring for electrical consumption managing. This approach can automatically detect the switch-on and switch-off events of domestic appliances and classify different appliances using load features and advanced algorithms, therefore, monitoring the house power consumption of individual devices. The new transient detection algorithm, in combination with turn-on and turn-off transient waveforms analysis, is developed to detect the mutative power events. The load classification technique, which employs support vector machine to recognize different appliances, is capable of identifying kinds of power loads with improved recognition accuracy and computational speed of NILM results. A recognition rate of higher than 90% has been achieved in recognizing various electrical devices. Our experiments have shown that recognition approaches based on support vector machines will be a trend for the process of residential load monitoring.



---

*Source: 742461-2012-11-19.xml*