Load Signature Formulation for Non-Intrusive Load Monitoring Based on Current Measurements
Abstract
:1. Introduction
2. Measurements and Load Signature Formulation
2.1. Measurement Setup
2.2. Measured Appliances
2.3. Load Signatures Formulation
- LSi is the load signature of appliance i;
- indices 50, 150 and 250 represent the fundamental nominal current, the 3rd and the 5th harmonic currents respectively;
- indices and with and denote the number of values utilized for the respective part of the LS.
- Which one of these values should be considered for the formulation of the first part of the LS in (1) since all of them describe the operation of this appliance under steady-state? A simple and quick approach could be the mean value. The problem here is that the higher the variation range, the less representative the mean value would be. This could greatly hinder the performance of a NILM algorithm and the identification accuracy of the appliance.
- If the consideration of a single value yields inefficient LS, then how many values should be utilized in order to ensure that the operation of the appliance is captured in most of the possible operational modes? Τhe answer in this question defines the number of a, b and c indices of the LS.
- Each appliance is measured for a time period of 5 min under a recording frequency of 1 sample/s as described in Section 2.1. Therefore, a data series with approximately 300 current values for each harmonic order (i.e., at 50 Hz, 150 Hz and 250 Hz) for each appliance are stored in the database. For most of the typical appliances in a residence the time period of 5 min can be considered adequate, since it captures the typical residential usage. For those with multiple operational modes, e.g., washing machine, all of these different operational modes should be measured, for the LS formulation.
- The standard deviation (SD) for each appliance is computed: SD50, SD150 and SD250 respectively for the three examined frequencies.
- A threshold (th) is defined for each SD in order to identify if one or more values should be utilized for the formulation of the respective part of the LS.
- This threshold is taken as follows: SD50(th) = SD150(th) = SD250(th) = 0.02. The value of the threshold (th) has been selected after several trials since this specific value has provided relatively short load signatures (i.e., with relatively few representative currents values) but efficient enough for load identification in the disaggregation mode of the proposed methodology. This threshold value is proposed as the upper limit regarding the data processing towards the LSs formulation.
- For each appliance i the following rules are applied:
- If SD50i ≤ SD50(th), then compute the mean value (μ50i) for the data in this data series and formulate the first part of the LSi as follows: I50i = μ50i. Obviously in this case the value of index a is equal to 1, a = 1.
- If SD50i > SD50(th), then reorder the data in the data series in descending order. Afterwards, divide top-down the data in a (a = 1, …, z) non-overlapping sequential data groups in order to ensure that for each one SD50ia ≤ SD50(th).
- For each of these a groups, compute the mean value (μ50ia). Formulate the first part of the LSi as follows: I50i = μ50ia, …, μ50iz.
- Apply steps 5a–5c to the data series for 150 Hz and 250 Hz respectively under the corresponding SD threshold that is defined in step 4. This obtains the values of indices b and c.
- Store the formulated LSs for the i appliances and form the LS database for this residence.
2.4. Load Signatures of Agreggated Measurements
- Given a measurement of the total instantaneous current at the main feeding panel of the residence at time t, i.e., It, how could we identify the operating appliances in time t? Thus, the challenge here is to efficiently disaggregate the measured value to its components parts of certain loads.
- Apply FFT to the measured instantaneous current values within time t
- Determine the current amplitudes for frequencies of 50 Hz, 150 Hz and 250 Hz
- The recorded values at time t is the mean of the current amplitudes
3. Proposed Methodology towards NILM Implementation
3.1. Disaggregation and NILM Algorithm
3.2. Disaggregation Results
- Each appliance displays a relatively constant operational mode with smooth variations regarding the recorded current values.
- Due to the high fundamental current of each appliance (i.e., I50), when they operate simultaneously the variations of the aggregated fundamental current values could be higher than the current of small rated appliances. In this case the algorithm may incorrectly identify not operating appliances. Usually, such small rated appliances (e.g., TV or laptop) are not linear loads with significant 3rd and 5th order harmonic currents. Hence, in this case the OOHCs are expected to contribute by excluding these appliances from the identification procedure.
3.3. Results Evaluation
- The methodology for the LSs formulation could provide few and still representative current values that adequately cover the steady state operation of the appliances. A sensitivity analysis about the predefined SD threshold for each frequency could provide the optimal number of utilized current values for the formulation of the three parts of each LS.
- The proposed NILM scheme could be considered suitable for near real-time load identification. A NILM scheme with such high successful identification resolution could yield a detailed disaggregation of the consumption behavior of a residence and is highly appreciated by the retail energy providers. For example, the more detailed the knowledge of the consumption behavior of the customers the more efficiently demand response schemes can be designed and implemented.
- The identification approach performs almost flawlessly for combinations that include high consuming appliances without significant harmonic content. The latter is quite evident in comb#1 and comb#2 since Equation (2) is valid for the short LSs.
- For combinations with many appliances that present significant harmonic currents, the efficiency of Equation (2) is limited when only the harmonic current amplitudes are considered. In this case, the phase angle of each harmonic current should be also considered (using the fundamental voltage phase angle as the angular reference) because the aggregation should refer to vectors and not just amplitudes. For example, comb#4 includes five appliances that all present harmonic behavior. The summation of the harmonic current amplitudes shows high deviations from the considered recorded aggregated value of the combination. The latter explains the poor identification rate of the Resistive-heater appliance, since the algorithm identifies the Coffee-maker and Hair-dryer appliances instead. The problem here is that the contribution of the 5th harmonic current of the Resistive-heater is not identified in the aggregated recorded value. The TV appliance is not identified due to the same reason as well. Based on measurements in [29] the phase angle between the 5th harmonic currents of an LCD TV and a desktop PC is approximately 330°, thus the amplitudes should be almost subtracted concerning the aggregated respective amplitude under simultaneous operation.
4. Conclusions
Acknowledgments
Author Contributions
Conflicts of Interest
References
- Paudyal, S.; Canizares, C.A.; Bhattacharya, K. Optimal Operation of Distribution Feeders in Smart Grids. IEEE Trans. Ind. Electron. 2011, 58, 4495–4503. [Google Scholar] [CrossRef]
- Bouhouras, A.S.; Andreou, G.T.; Labridis, D.P.; Bakirtzis, A.G. Selective automation upgrade in distribution networks towards a smarter grid. IEEE Trans. Smart Grid 2010, 1, 278–285. [Google Scholar] [CrossRef]
- Mehrshad, M.; Tafti, A.D.; Effatnejad, R. Demand side Management in the Smart Grid Based on Energy Consumption Scheduling by NSGA-II. Int. J. Eng. Pract. Res. 2013, 2, 197–200. [Google Scholar]
- Chang, H.H.; Lin, C.L.; Lee, J.K. Load identification in nonintrusive load monitoring using steady-state and turn-on transient energy algorithms. In Proceedings of the 2010 14th International Conference on Computer Supported Cooperative Work in Design (CSCWD), Shanghai, China, 14–16 April 2010; pp. 27–32. [Google Scholar]
- Liang, J.; Ng, S.K.K.; Kendall, G.; Cheng, J.W.M. Load signature studypart I: Basic concept, structure, and methodology. IEEE Trans. Power Deliv. 2010, 25, 551–560. [Google Scholar] [CrossRef]
- Hart, G.W. Nonintrusive Appliance Load Monitoring. Proc. IEEE 1992, 80, 1870–1891. [Google Scholar] [CrossRef]
- Laughman, C.; Lee, K.; Cox, R.; Shaw, S.; Leeb, S.; Norford, L.; Armstrong, P. Power signature analysis. IEEE Power Energy Mag. 2003, 1, 56–63. [Google Scholar] [CrossRef]
- Lin, Y.-H.; Tsai, M.-S. A novel feature extraction method for the development of nonintrusive load monitoring system based on BP-ANN. In Proceedings of the 2010 International Symposium on Computer, Communication, Control and Automation, Tainan, Taiwan, 5–7 May 2010; Volume 2, pp. 215–218. [Google Scholar]
- Yang, H.T.; Chang, H.H.; Lin, C.L. Design a neural network for features selection in non-intrusive monitoring of industrial electrical loads. In Proceedings of the 2007 11th International Conference on Computer Supported Cooperative Work in Design (CSCWD), Melbourne, Australia, 26–28 April 2007; pp. 1022–1027. [Google Scholar]
- Ng, S.K.K.; Liang, J.; Cheng, J.W.M. Automatic appliance load signature identification by statistical clustering. In Proceedings of the 8th International Conference on Advances in Power System Control, Operation and Management (APSCOM 2009), Hong Kong, China, 8–11 November 2009; pp. 1–6. [Google Scholar]
- Chang, H.H.; Chen, K.; Tsai, Y.P.; Lee, W.J. A New Measurement Method for Power Signatures of Nonintrusive Demand A New Measurement Method for Power Signatures of Nonintrusive Demand Monitoring and Load Identification. IEEE Trans. Ind. Appl. 2012, 48, 764–771. [Google Scholar] [CrossRef]
- Du, Y.; Du, L.; Lu, B.; Harley, R.; Habetler, T. A review of identification and monitoring methods for electric loads in commercial and residential buildings. In Proceedings of the 2010 IEEE Energy Conversion Congress and Exposition (ECCE), Atlanta, GA, USA, 12–16 September 2010; pp. 4527–4533. [Google Scholar]
- Shao, H.; Tech, V.; Marwah, M. A Temporal Motif Mining Approach to Unsupervised Energy Disaggregation. In Proceedings of the 1st International Workshop on Non-Intrusive Load Monitoring, Pittsburgh, PA, USA, 7 May 2012; pp. 1–2. [Google Scholar]
- Li, J.; Yang, H. The investigation of residential load identification technologies. In Proceedings of the 2012 Asia-Pacific Power and Energy Engineering Conference (APPEEC), Shanghai, China, 27–29 March 2012; pp. 2–5. [Google Scholar]
- Srinivasan, D.; Ng, W.S.; Liew, A.C. Neural-network-based signature recognition for harmonic source identification. IEEE Trans. Power Deliv. 2006, 21, 398–405. [Google Scholar] [CrossRef]
- Ehrhardt-martinez, A.K.; Donnelly, K.A. Advanced Metering Initiatives and Residential Feedback Programs: A Meta-Review for Household Electricity-Saving Opportunities. Energy 2010, 123, 128. [Google Scholar]
- Chang, H.H. Non-intrusive demand monitoring and load identification for energy management systems based on transient feature analyses. Energies 2012, 5, 4569–4589. [Google Scholar] [CrossRef]
- Bacurau, R. Experimental Investigation on the Load Signature Parameters for Non-Intrusive Load Monitoring. Prz. Elektrotech. 2015, 1, 88–92. [Google Scholar] [CrossRef]
- Hong, Y.Y.; Chou, J.H. Nonintrusive energy monitoring for microgrids using hybrid self-organizing feature-mapping networks. Energies 2012, 5, 2578–2593. [Google Scholar] [CrossRef]
- Zoha, A.; Gluhak, A.; Imran, M.A.; Rajasegarar, S. Non-intrusive Load Monitoring approaches for disaggregated energy sensing: A survey. Sensors 2012, 12, 16838–16866. [Google Scholar] [CrossRef] [PubMed]
- Huang, S.J.; Hsieh, C.T.; Kuo, L.C.; Lin, C.W.; Chang, C.W.; Fang, S.A. Classification of home appliance electricity consumption using power signature and harmonic features. Proceedings of 2011 IEEE Ninth International Conference on Power Electronics and Drive Systems, Singapore, 5–8 December 2011; pp. 596–599. [Google Scholar]
- Bouhouras, A.S.; Milioudis, A.N.; Andreou, G.T.; Labridis, D.P. Load signatures improvement through the determination of a spectral distribution coefficient for load identification. In Proceedings of the 2012 9th International Conference on the European Energy Market (EEM), Florence, Italy, 10–12 May 2012; pp. 1–6. [Google Scholar]
- Bouhouras, A.S.; Milioudis, A.N.; Labridis, D.P. Development of distinct load signatures for higher efficiency of NILM algorithms. Electr. Power Syst. Res. 2014, 117, 163–171. [Google Scholar] [CrossRef]
- Semwal, S.; Shah, G.; Prasad, R.S. Identification residential appliance using NIALM. In Proceedings of the 2014 IEEE International Conference on Power Electronics, Drives and Energy Systems, Mumbai, India, 16–19 December 2014. [Google Scholar]
- Srinivasan, G.; Anandan, C.; Jain, S.A.K.; Ahmed, S.S.; Vijayaraghavan, V. Low-cost Non-Intrusive Device Identification System. In Proceedings of the 2016 IEEE Circuits and Systems Conference (DCAS), Arlington, TX, USA, 10 October 2016; pp. 4–7. [Google Scholar]
- Feng, C.Y.; Hoe, H.M.; Abdullah, P. Tracing of Energy Consumption by Using Harmonic Current. In Proceedings of the 2013 IEEE Student Conference on Research and Development (SCOReD), Putrajaya, Malaysia, 16–17 December 2013; pp. 16–17. [Google Scholar]
- Tran, T.T.; Lee, G.-D.; Pham, T.X.; Kim, G.-J.; Van Dang, C.; Kim, J.-W.; Kang, B. Identification of In-Home Appliances Through Analysis of Current Consumption. In Proceedings of the 10th International Conference on Ubiquitous Information Management and Communication, Danang, Vietnam, 4–6 January 2016. [Google Scholar]
- Gil-de-castro, A.; Moreno-muñoz, A. Street Lamps Aggregation Analysis through IEC 61000-3-6 Approach. In Proceedings of the 22nd International Conference on Electricity Distribution, Stockholm, Sweden, 10–13 June 2013; pp. 10–13. [Google Scholar]
- Nassif, A.B.; Yong, J.; Xu, W.; Chung, C.Y. Indices for comparative assessment of the harmonic effect of different home appliances. Eur. Trans. Electr. Power 2014, 23, 638–654. [Google Scholar] [CrossRef]
Measured Appliances | ||||
PC | Hair-dryer (hot) | Coffee-maker | Vacuum | Electric iron |
TV | Resistive-heater | Toaster | Refrigerator | Blender |
Measured Appliance Combinations | ||||
1 | Hair-dryer (hot) | Vacuum | Resistive-heater | |
2 | Coffee-maker | Electric iron | Resistive-heater | Toaster |
3 | Coffee-maker | Toaster | Refrigerator | Blender |
4 | PC | Electric iron | TV | Resistive-heater |
Appliance | Recording Frequency 1 Value/S Recording Time Approximately 5 min | ||
---|---|---|---|
SD50 | SD150 | SD250 | |
PC-desktop | 0.210 | 0.173 | 0.110 |
Hair dryer (hot) | 0.031 | 0.005 | 0.011 |
Coffee maker 1 | 0.020 | * | * |
Vacuum | 0.086 | 0.046 | 0.013 |
Electric iron 1 | 0.227 | * | 0.007 |
TV | 0.010 | 0.006 | 0.007 |
Resistive Heater 1 | 0.039 | * | 0.020 |
Toaster 1 | 0.009 | * | * |
PC-laptop | 0.165 | 0.107 | 0.032 |
Refrigerator | 0.016 | * | * |
Blender | 0.011 | * | * |
LS Database for the Examined Residence | |||||||
---|---|---|---|---|---|---|---|
Appliance | 1st Part for 50 Hz Current Amplitudes [A] | 2nd Part for 150 Hz Current Amplitudes [A] | 3rd Part for 250 Hz Current Amplitudes [A] | Indices Values a-b-c | |||
PC (desktop) | 1.122 | 1.154 | 1.078 | 1.022 | 0.818 | 0.739 | 10-10-9 |
1.082 | 1.014 | 0.946 | 0.891 | 0.703 | 0.654 | ||
0.940 | 0.865 | 0.856 | 0.774 | 0.579 | 0.516 | ||
0.800 | 0.736 | 0.702 | 0.661 | 0.416 | 0.384 | ||
0.630 | 0.588 | 0.563 | 0.522 | 0.346 | |||
Hair dryer (hot) | 11.328 | 11.273 | 0.084 | 0.180 0.147 | 4-1-2 | ||
11.162 | 11.136 | ||||||
Coffee maker | 4.668 | 4.570 | negligible | negligible | 4-0-0 | ||
4.425 | 4.389 | ||||||
Vacuum | 6.588 | 6.514 | 1.281 | 1.228 | 0.123 0.100 | 8-5-2 | |
6.432 | 6.363 | 1.190 | 1.150 | ||||
6.300 | 6.225 | 1.134 | |||||
6.182 | 6.009 | ||||||
Electric iron | 12.841 | 12.754 | negligible | 0.186 | 5-0-1 | ||
12.687 | 12.598 | ||||||
12.280 | |||||||
TV | 0.100 | 0.175 | 0.142 | 0.167 | 0.150 | 1-2-2 | |
Resistive Heater | 14.940 | 14.805 | negligible | 0.234 0.202 | 4-0-2 | ||
14.693 | 14.607 | ||||||
Toaster | 2.721 | 2.670 | negligible | negligible | 2-0-0 | ||
Refrigerator | 0.940 | 0.860 | negligible | negligible | 2-0-0 | ||
Blender | 0.483 | 0.465 | negligible | negligible | 2-0-0 |
comb#1 | ||||
Appliance | Time Activated (s) | Correctly Identified (s) | NILM Performance | Min and Max Value of f |
Vacuum | 157 | 157 | 100% | 0.044 0.423 |
Hair dryer (hot) | 109 | 108 | 99% | |
Resistive heater | 52 | 49 | 94% | |
comb#2 | ||||
Resistive heater | 215 | 208 | 97% | 0.010 0.233 |
Coffee maker | 78 | 78 | 100% | |
Electric iron | 19 | 19 | 100% | |
Toaster | 56 | 56 | 100% | |
2 times (at seconds 105 and 174) the algorithm identifies hairdryer-discarded. | ||||
comb#3 | ||||
Refrigerator | 243 | 195 | 80% | 0.053 0.152 |
Coffee maker | 97 | 97 | 100% | |
Blender | 135 | 66 | 49% | |
Toaster | 83 | 83 | 100% | |
failure to identify mixer is several cases—no irrelevant appliances proposed by the algorithm. | ||||
comb#4 | ||||
PC (desktop) | 300 | 163 | 54% | 0.015 0.850 |
TV | 240 | 0 | 0% | |
Vacuum | 192 | 192 | 100% | |
Resistive heater | 132 | 0 | 0% | |
Electric iron | 24 | 24 | 100% | |
failure to identify TV (low current at 50 Hz falls within the variation ranges); failure to identify resistive heater due to the fact that Equation (2) was not valid for the aggregated current values of the 3rd and 5th harmonic currents. |
© 2017 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
Share and Cite
Bouhouras, A.S.; Gkaidatzis, P.A.; Chatzisavvas, K.C.; Panagiotou, E.; Poulakis, N.; Christoforidis, G.C. Load Signature Formulation for Non-Intrusive Load Monitoring Based on Current Measurements. Energies 2017, 10, 538. https://doi.org/10.3390/en10040538
Bouhouras AS, Gkaidatzis PA, Chatzisavvas KC, Panagiotou E, Poulakis N, Christoforidis GC. Load Signature Formulation for Non-Intrusive Load Monitoring Based on Current Measurements. Energies. 2017; 10(4):538. https://doi.org/10.3390/en10040538
Chicago/Turabian StyleBouhouras, Aggelos S., Paschalis A. Gkaidatzis, Konstantinos C. Chatzisavvas, Evangelos Panagiotou, Nikolaos Poulakis, and Georgios C. Christoforidis. 2017. "Load Signature Formulation for Non-Intrusive Load Monitoring Based on Current Measurements" Energies 10, no. 4: 538. https://doi.org/10.3390/en10040538
APA StyleBouhouras, A. S., Gkaidatzis, P. A., Chatzisavvas, K. C., Panagiotou, E., Poulakis, N., & Christoforidis, G. C. (2017). Load Signature Formulation for Non-Intrusive Load Monitoring Based on Current Measurements. Energies, 10(4), 538. https://doi.org/10.3390/en10040538