Reinforcement Learning Applied to an Electric Water Heater: From Theory to Practice

2015·Arxiv

Abstract

Abstract

Electric water heaters have the ability to store energy in their water buffer without impacting the comfort of the end user. This feature makes them a prime candidate for residential demand response. However, the stochastic and non-linear dynamics of electric water heaters, makes it challenging to harness their flexibility. Driven by this challenge, this paper formulates the underlying sequential decision-making problem as a Markov decision process and uses techniques from reinforcement learning. Specifically, we apply an auto-encoder network to find a compact feature representation of the sensor measurements, which helps to mitigate the curse of dimensionality. A wellknown batch reinforcement learning technique, fitted Q-iteration, is used to find a control policy, given this feature representation. In a simulation-based experiment using an electric water heater with 50 temperature sensors, the proposed method was able to achieve good policies much faster than when using the full state information. In a lab experiment, we apply fitted Q-iteration to an electric water heater with eight temperature sensors. Further reducing the state vector did not improve the results of fitted Q-iteration. The results of the lab experiment, spanning 40 days, indicate that compared to a thermostat controller, the presented approach was able to reduce the total cost of energy consumption of the electric water heater by 15%.

Index Terms—Auto-encoder network, demand response, electric water heater, fitted Q-iteration, reinforcement learning.

I. INTRODUCTION

THE share of renewable energy sources is expected toreach 25% of the global power generation portfolio by 2020 [1]. The intermittent and stochastic nature of most renewable energy sources, however, makes it challenging to integrate these sources into the power grid. Successful integration of these sources requires flexibility on the demand side through demand response programs. Demand response programs enable end users with flexible loads to adapt their consumption profile in response to an external signal. A prominent example of flexible loads are electric water heaters with a hot water storage tank [2], [3]. These loads have the ability to store energy in their water buffer without impacting the comfort level of the end user. In addition to having significant flexibility, electric water heaters can consume about 2 MWh per year for a household with a daily hot water demand of 100 liters [4]. As a result, electric water heaters are a prime candidate for residential demand response programs.

F. Ruelens and R. Belmans are with the Department of Electrical Engineering, KU Leuven/EnergyVille, Belgium (frederik.ruelens@esat.kuleuven.be).

B. J. Claessens is with the Energy Department of Vito/EnergyVille, Belgium (bert.claessens@vito.be).

S. Quaiyum is with the Department of Electrical Engineering, Uppsala University, Sweden.

B. De Schutter and R. Babuˇska are with the Delft Center for Systems and Control, Delft University of Technology, The Netherlands.

Previously, the flexibility offered by electric water heaters has been used for frequency control [5], local voltage control [6], and energy arbitrage [7], [8]. Amongst others, two prominent control paradigms in the demand response literature on electric water heaters are model-based approaches and reinforcement learning.

Perhaps the most researched control paradigm applied to demand response are model-based approaches, such as Model Predictive Control (MPC) [7], [9], [8]. Most MPC strategies use a gray-box model, based on general expert knowledge of the underlying system dynamics, requiring a system iden-tification step. Given this mathematical model, an optimal control action can be found by solving a receding horizon problem [10]. In general, the implementation of MPC consists of several critical steps, namely, selecting accurate models, estimating the model parameters, estimating the state of the system, and forecasting of the exogenous variables. All these steps make MPC an expensive technique, the cost of which needs to be balanced out by the possible financial gains [11]. Moreover, possible model errors resulting from an inaccurate model or forecast, can effect the stability of the MPC controller [12], [13].

In contrast to MPC, Reinforcement Learning (RL) techniques [14] do not require expert knowledge and consider their environment as a black-box. RL techniques enable an agent to learn a control policy by interacting with its environment, without the need to use modeling and system identification techniques. In [15], Ernst et al. state that the trade-off in applying MPC and RL mainly depends on the quality of the expert knowledge about the system dynamics that could be exploited in the context of system identification. In most residential demand response applications, however, expert knowledge about the system dynamics or future disturbances might be unavailable or might be too expensive to obtain relative to the expected financial gain. In this context, RL techniques are an excellent candidate to build a general purpose agent that can be applied to any demand response application.

This paper proposes a learning agent that minimizes the cost of energy consumption of an electric water heater. The agent measures the state of its environment through a set of sensors that are connected along the water buffer of the electric water heater. However, learning in a high-dimension state space can significantly impact the learning rate of the RL algorithm. This is known as the “curse of dimensionality”. A popular approach to mitigate its effects is to reduce the dimensionality of the state space during a pre-processing step. Inspired by the work of [16], this paper applies an auto-encoder network to reduce the dimensionality of the state vector. By so doing, this paper makes following contributions: (1) we demonstrate how a well-established RL technique, fitted Q-iteration, can be combined with an auto-encoder network to minimize the cost of energy consumption of an electric water heater; (2) in a simulation-based experiment, we assess the performance of different state representations and batch sizes; (3) we successfully apply an RL agent to an electric water heater in the lab (Fig. 1).

The remainder of this paper is organized as follows. Section II gives a non-exhaustive literature overview of RL related to demand response. Section III maps the considered demand response problem to a Markov decision process. Section IV describes how an auto-encoder network can be used to find a low-dimensional state representation, followed by a description of the fitted Q-iteration algorithm in Section V. Section VI presents the results of the simulation-based experiments and Section VII presents the results of the lab experiment. Finally, Section VIII draws conclusions and discusses future work.

II. REINFORCEMENT LEARNING

This section gives a non-exhaustive overview of recent developments related to Reinforcement Learning (RL) and demand response. Perhaps the most widely used model-free RL technique applied to a demand response setting is standard Q-learning [17], [18], [19], [20]. After each interaction with the environment, Q-learning uses temporal difference learning [14] to update its state-action value function or Qfunction. A major drawback of Q-learning is that the given observation is discarded after each update. As a result, more interactions are needed to spread already known information through the state space. This inefficient use of information limits the application of Q-learning to real-world applications.

In contrast to Q-learning, batch RL techniques [21], [22], [23] are more data efficient, since they store and reuse past interactions. As a result, batch RL techniques require less interactions with their environment, which makes them more practical for real-world applications, such as demand response. Perhaps the most popular batch RL technique which has been applied to a wide range of applications [16], [24], [25] is fitted Q-iteration developed by Ernst et al. [21]. Fitted Q-iteration iteratively estimates the Q-function given a fixed batch of past interactions. An online version that uses a neural network, neural fitted Q-iteration, has been proposed by Riedmiller et al. in [22]. Finally, an interesting alternative is to combine experience replay to an incremental RL technique such as Qlearning or SARSA [26]. In [27], the authors demonstrate how fitted Q-iteration can be used to control a cluster of electric water heaters. The results indicate that fitted Q-iteration was able to reduce the cost of energy consumption of a cluster of 100 electric water heaters after a learning period of 40 days. In addition, [28] shows how fitted Q-iteration can be extended to reduce the cost of energy consumption of a heat-pump thermostat given that a forecast of the outside temperature is provided.

A promising alternative to the previously mentioned model-free techniques are model-based or model-assisted RL techniques. For example, the authors of [29] present a model-based

Fig. 1. Setup of the electric water heater used during the lab experiment. Eight temperature sensors were placed under the insulation material of the buffer tank.

policy search method that learns a Gaussian process to model uncertainties. In addition, inspired by [30], the authors of [31] demonstrate how a model-assisted batch RL technique can be applied to control a building heating system.

III. PROBLEM FORMULATION

The aim of this work is to develop a controller or agent that minimizes the cost of energy consumption of an electric water heater, given an external price profile. This price profile is provided to the agent at the start of each day. The agent can measure the temperature of the water buffer through a set of temperature sensors that are connected along the hull of the buffer tank. Following the approach presented in [28], the electric water heater is equipped with a backup controller that overrules the control action from the agent when the safety or comfort constraints of the end user are violated. A challenge in developing such an agent is that the dynamics of the electric water heater, the future hot water demand and the settings of the backup controller are unknown to the agent. To overcome this challenge, this paper leverages on the previous work of [16], [21], [28] and applies techniques from RL.

A. Markov decision process framework

To apply RL, this paper formulates the underlying sequential decision-making problem of the learning agent as a Markov decision process formulation. The Markov decision process formulation is defined by its d-dimensional state space , its action space , its stochastic discrete-time transition function f, and its cost function . The optimization horizon is considered finite, comprising steps, where at each discrete time step k, the state evolves following:

where a random disturbance is drawn from a conditional distribution is the control action and X the state. Associated with each state transition, a cost is given by:

The goal of the learning agent is to find an optimal control policy that minimizes the expected T-stage return for any state in the state space. The expected T-stage return starting from and following a policy h is defined as follows:

where E denotes the expectation operator over all possible stochastic realizations.

A more convenient way to characterize a policy is by using a state-action value function or Q-function:

which indicates the cumulative return starting from state x and by taking action u and following h thereafter.

The optimal Q-function corresponds the best Q-function that can be obtained by any policy:

Starting from an optimal Q-function for every state-action pair, the optimal policy is calculated as follows:

where satisfies the Bellman optimality equation [32]:

Following the notation introduced in [28], the next three paragraphs give a description of the state, the action, and the cost function tailored to an electric water heater.

B. Observable state vector

The observable state vector of an electric water heater contains a time-related component and a controllable component. The time-related component describes the part of the state related to timing, which is relevant for the dynamics of the system. Specifically, the tap water demand of the end user is considered to have diurnal and weekly patterns. As such, the time-related component contains the day of the week and the quarter in the day. The controllable component represents physical state information that is measured locally and is influenced by the control action. The controllable component contains the temperature measurements of the sensors that are connected along the hull of the storage tank. The observable state vector is given by:

where is the current day of the week, is the quarter in the day and denotes the temperature measurements of sensor i at time step k.

C. Control action

The learning agent can control the heating element of the electric water heater with a binary control action , where 0 indicates off and 1 on. However, the backup mechanism, which enacts the comfort and safety constraints of the end user, can overrule this control action of the learning agent. The function maps the control action to a physical action according to:

where the vector defines the safety and user-defined comfort settings of the backup controller. In order to have a generic approach we assume that the logic of the backup controller is unknown to the learning agent. However, the learning agent can measure the physical action enforced by the backup controller (see Fig. 2), which is required to calculate the cost.

The logic of the backup controller of the electric water heater is defined as:

where is the electrical power rating of the heating element, is the current state of charge and and are the upper and lower bounds for the state of charge. A detailed description of how the state of charge is calculated can by found in [3].

D. Cost function

At the start of each optimization period , the learning agent receives a price vector for the next T time steps. At each time step, the agent receives a cost according to:

where is the electricity price during time step k, and the length of one control period.

IV. BATCH OF FOUR-TUPLES

Generally, batch RL techniques estimate the Q-function based on a batch of four-tuples . This paper, however, considers the following batch of four-tuples:

where for each l, the next state , and the physical action have been obtained as a result of taking control action in the state . Note that, F does not include the observed cost , since the cost depends on the price vector that is provided to the learning agent at the start of each day.

As defined by (8), contains all temperature measurements of the sensors connected to the hull of the water buffer. Learning in a high-dimensional state space requires more observations from the environment to estimate the Q-function, as more tuples are needed to cover the state-action space. This is known as the “curse of dimensionality”. This curse becomes even more pronounced in practical applications,

where each observation corresponds to a “real” interation with the environment.

A pre-processing step can be used to find a compact and more efficient representation of the state space and can help to converge to a good policy much faster [33]. A popular technique to find a compact representation is to use a handcrafted feature vector based on insights of the considered control problem [34]. Alternative approaches that do not require prior knowledge are unsupervised feature learning algorithms, such as auto-encoders [16] or a principal component analysis [33].

As illustrated in Fig. 2, this paper demonstrates how an auto-encoder can be used to find a compact representation of the sensory input data. An auto-encoder network is a neural network that maps its output back to its input. By selecting a lower number of neurons in the middle hidden layer than in the input layer p < d, the auto-encoder can be used to reduce the dimensionality of the input data.

The reduced feature vector can be calculated as follows:

where w and b denote the weights and the biases that connect the input layer with the middle hidden layer of the auto-encoder network. The function is an encoder function and maps the observed state vector to the feature vector . To train the weights of the auto-encoder, a conjugate gradient descent algorithm is used as presented in [35].

In the next section, fitted Q-iteration is used to find a policy that maps every feature vector to a control action using batch R:

which consists of feature vectors with a dimensionality p.

Since we apply the auto-encoder on the input data of the supervised learning algorithm, we assume that all input data is equally important. As such, it is possible that we ignore low-variance yet potentially useful components during the learning process. A possible route of future work would be to add a regularization term to the regression algorithm of the supervised learning algorithm to prevent the risk of overfitting without the risk of ignoring potentially important data.

Fig. 2. Setup of the simulation-based experiment. An auto-encoder network is used to find a compact representation of the temperature measurements.

V. FITTED Q-ITERATION

This section describes the learning algorithm and the ex- ploration strategy of the agent based on the batch of feature vectors R presented in the previous section.

A. Fitted Q-iteration

Fitted Q-iteration iteratively builds a training set with all state-action pairs (z, u) in R as the input. The target values consist of the corresponding Q-values, based on the approximation of the Q-function of the previous iteration. In the first iteration (N = 1), the Q-values approximate the expected cost (line 5 in Algorithm 1). In the subsequent iterations, Q-values are updated using the Q-function of the previous iteration. As a result, Algorithm 1 needs T iterations until the Q-function contains all information about the future costs. Note that, the cost corresponding to each tuple is recalculated using price vector that is provided at the start of the day (line 4 in Algorithm 1). As a result, the algorithm can reuse past experiences to find a control policy for the next day. Following [21], Algorithm 1 applies an ensemble of extremely randomized trees as a regression algorithm to estimate the Qfunction. An empirical study of the accuracy and convergence properties of extremely randomized trees can be found in [21]. However, in principle, any regression algorithm, such as neural networks [23], [24], can be used to estimate the Q-function.

B. Boltzmann exploration

During the day, the learning agent uses a Boltzmann exploration strategy [36] and selects an action with the following probability:

where is the Boltzmann temperature at day d, is the Qfunction from Algorithm 1 and z is the current feature vector measured by the learning agent. If , the exploration will decrease and the policy will become greedier. Thus by starting with a high the exploration starts completely random,

3 5 10 25 50 75 100 150 200 2.6

Fig. 3. Simulation-based results of fitted Q-iteration using five state represen- tations and different batch sizes. The full state contains 50 temperature measurements. A non-linear dimensionality reduction with Auto-Encoder (AE) is used to find a compact representation of the temperature measurements. Each marker point represents the average result of 100 simulation runs.

however as decreases the policy directs itself to the most interesting state-action pairs. In the evaluation experiments, in (15) is linearly scaled between [0, 100] and the is set to 100 at the start of the experiment, which will result in an equal probability for all actions. The Boltzmann temperature is updated as follows , which increases the probability of selecting higher valued actions.

VI. SIMULATION-BASED RESULTS

This section describes the results of the simulation-based experiments, which use a non-linear stratified tank model with 50 temperature layers. A detailed description of the stratified tank model can be found in [3]. The specifications of the electric water heater are chosen in correspondence with the electric water heater used during the lab experiment (see Section VII). The simulated electric water heater has a power rating of 2.36kW and has a water buffer of 200 liter. The experiments use realistic hot water profiles with a mean daily consumption of 120 liter [37] and use price information from the Belgian day-ahead [38] and balancing market [39]. The learning agent can measure the temperature of the 50 temperature layers obtained with the simulation model.

The aim of the first simulation-based experiment is to find a compact state representation using an auto-encoder network and to assess the impact of the state representation on the performance of fitted Q-iteration. The second simulation-based experiment compares the result of fitted Q-iteration with the default thermostat controller.

A. Step 1: feature selection

This experiment compares the performance of fitted Qiteration combined with different feature representations for different fixed batch sizes. An auto-encoder (AE) network

Fig. 4. Cumulative energy cost of fitted Q-iteration with a non-linear dimensionality reduction and of the thermostat controller. Results for one year using day-ahead prices (a) and imbalance prices (b).

that reduces the original sensory input vector (50 temperature sensors) to 5 dimensions is denoted by AE 5. The simulations are repeated for 100 simulation days. The average cost of energy consumption of these 100 simulations is depicted in Fig. 3. As can be seen in Fig. 3, the performance of fitted Qiteration combined with a specific state representation depends on the number of tuples in the batch. For example for a batch size of 10 days, AE 3 results in a lower cost than AE 15, while after 75 days, AE 15 will result in a lower cost than AE 3. In addition, as can be seen from Fig. 3, AE 1 resulted in a relatively bad policy, independent of the batch size.

In general, it can be concluded that for a batch of limited size, fitted Q-iteration with a low-dimensional feature vector will outperform fitted Q-iteration using the full state information, i.e. 50 temperature measurements. By learning in a low-dimensional state space, it is possible to learn with a smaller and more efficient representation. As a result, the agent requires less observations to converge to a better control policy than when the full state information is used. In addition, as more observations will result in a more efficient coverage of the state-action space, it can be seen from Fig. 3 that the result of fitted Q-iteration with the full state improves significantly as the batch size increases. In the following subsection, we present the results of AE 5 in more detail. A method for selecting an appropriate feature representation during the learning process will be part of our future work (Section VIII).

B. Step 2: evaluation

Fig. 4 compares the total cost of energy consumption using fitted Q-iteration combined with AE 5 against the default thermostat controller for two relevant price profiles, i.e. day-ahead prices (top plot) and imbalance prices (bottom plot). The default thermostat controller enables the heating element

0 48 96 144 192 240 288 10

0 48 96 144 192 240 288 0

Fig. 5. Simulation-based results of a mature agent using fitted Q-iteration with a non-linear dimensionality reduction. a, Temperature profiles of the 50 simulation layers. b, Power consumption (black) and imbalance prices (gray).

when the state-of-charge drops below its minimum threshold and stays enabled until the state-of-charge reaches 100%. Note that in contrast to the learning agent, the default controller is agnostic about the price profile.

The experiment starts with an empty batch and the tuples of the current day are added to the given batch at the end of each day. At the start of each day, the auto-encoder is trained to find a batch of compact feature vectors, which are then used by fitted Q-iteration to estimate the Q-function for the next day. Online, the learning agent uses a Boltzmann exploration strategy with set to 10, which results in 10 days of exploration.

The results of the experiment indicate that fitted Q-iteration was able to reduce the total cost of energy consumption by 24% for the day-head prices and by 34% for the imbalance prices compared to the default strategy. Note that the imbalance prices are generally more volatile than the day-ahead prices, as they reflect real-time imbalances due to forecasting errors of renewable energy sources, such as wind and solar, which were not foreseen in the day-ahead market.

The temperature profiles of the simulation layers and power profiles of a “mature” learning agent (batch size of 100 days) for the day-ahead and imbalance prices are depicted in Fig. 5 and Fig. 6. It can be seen in the bottom plot of both figures that the mature learning agent reduces the cost of energy consumption by consuming energy during low price moments.

VII. LAB RESULTS

The aim of our lab experiment is to demonstrate that fitted Q-iteration can be successfully applied to minimize the cost of energy consumption of a real-world electric water heater.

A. Lab setup

The setup used in the lab experiment was part of a pilot project on residential demand response in Belgium [40], where

0 48 96 144 192 240 288 10

0 48 96 144 192 240 288 0

Fig. 6. Simulation-based results of a mature agent using fitted Q-iteration with a non-linear dimensionality reduction. a, Temperature profiles of the 50 simulation layers. b, Power consumption (black) and day-ahead prices (gray).

a cluster of 10 electric water heaters was used for direct load control. Fig. 1 shows the electric water heater used during the lab experiment. The electric water heater is a standard unit that was equipped with eight temperature sensors and a controllable power relay. A controllable valve connected to the outlet of the buffer tank is used to simulate the hot water demand of a household with a mean daily flow volume of 120 liter [37]. An Arduino prototyping platform with a JSON/RPC 2.0 interface is used to communicate with a computer in the lab1, which runs the learning agent that uses fitted Qiteration. Fitted Q-iteration is implemented in Python and Scikit-learn [41] is used to estimate the Q-function, using an ensemble of extremely randomized trees [21].

Similar as in the previous simulation-based experiments, it is assumed that the learning agent is provided with a deterministic external price profile for the following day. The learning agent uses a Boltzmann exploration strategy with set to 10, which results in 10 days of exploration. In order to compare the performance of the lab experiment with the performance of the simulation-based experiments, we used both day-ahead prices and imbalance prices.

B. Evaluation

The performance of the learning agent was evaluated using different feature vectors as presented in Section VI. The best performance, however, was obtained by including the eight temperature measurements in the observable state vector.

Using identical price profiles and tap water demands, Fig. 7 and Fig. 8 show the temperature measurements and the power profiles of the mature learning agent using imbalance prices and day-ahead prices. As can be seen, the learning agent was able to successfully minimize the cost of energy consumption by consuming during low price moments.

0 48 96 144 192 240 288 20

0 48 96 144 192 240 288 0

Time [15 min] 0 48 96 144 192 240 2880

Fig. 7. Lab results of a mature agent using fitted Q-iteration with the full state (eight temperature measurements). a, Measurements of the temperature sensors. b, Power consumption (black) and imbalance prices (gray).

Fig. 9 depicts the experimental results, spanning 40 days, of fitted Q-iteration and the default thermostat controller. The top plot of this figure indicates the cumulative costs of energy consumption and the bottom plot indicates the daily costs of energy consumption. After 40 days, fitted Q-iteration was able to reduce the cost of energy consumption by 15% compared to the default thermostat controller. Furthermore, by excluding the first ten exploration days, this reduction increases to 28%.

VIII. CONCLUSIONS AND FUTURE WORK

This paper has demonstrated how an auto-encoder network can be used in combination with a well-established batch reinforcement learning algorithm, called fitted Q-iteration, to reduce the cost of energy consumption of an electric water heater. The auto-encoder network was used to find a compact representation of the state vector. In a series of simulation-based experiments using an electric water heater with 50 temperature sensors, the proposed method was able to converge to good policies much faster than when using the full state information. Compared to a default thermostat controller, the presented approach has reduced the cost of energy consumption by 24% using day-ahead prices and by 34% using imbalance prices.

In a lab experiment, fitted Q-iteration has been successfully applied to an electric water heater with eight temperature sensors. A reduction of the state vector did not improve the performance of fitted Q-iteration. Compared to the thermostat controller, fitted Q-iteration was able to reduce the total cost of energy consumption by 15% within 40 days of operation.

Based on the results of both experiments the following four conclusions can be drawn: (1) learning in a compact feature space can improve the quality of the control policy when the number of observations is relatively small (25 days); (2) when the number of observations increases it is advisable to switch to higher state-space representation; (3) when only a limited number of temperature sensors is available, i.e. 1-10 sensors, it

0 48 96 144 192 240 288 20

0 48 96 144 192 240 288 0

Time [15 min] 0 48 96 144 192 240 28820

Fig. 8. Lab results of a mature agent using fitted Q-iteration with the full state (eight temperature measurements). a, Measurements of the temperature sensors. b, Power consumption (black) and day-ahead prices (gray).

0 5 10 15 20 25 30 35 40 0

Fig. 9. Lab results of the learning agent using fitted Q-iteration (dashed line) and of the default thermostat (solid line) during 40 days. a, Cumulative energy cost. b, Daily energy cost.

is recommended to use the full state vector; (4) when applied to a real-world stetting, fitted Q-iteration was able to obtain good control policies within 20 days of operation.

In our future research, we aim at developing a method for selecting an appropriate state representation during the learning process. A promising route is to construct experts, where each expert combines a learning algorithm with a different feature representation. A metric based on the performance of each expert, as presented in [42], could then be used to select the expert with the highest metric as described in [43].

ACKNOWLEDGMENT

The authors would like to thank Davy Geysen, Geert Jacobs, Koen Vanthournout, and Jef Verbeeck from Vito for providing us with the lab setup. This work was supported by a Ph.D. grant of the Institute for the Promotion of Innovation through Science and Technology in Flanders (IWT-Vlaanderen) and by Stable MultI-agent LEarnIng for neTworks (SMILE-IT).

REFERENCES

[1] F. Birol et al., “World Energy Outlook 2013: Renewable Energy Outlook, An annual report released by the International Energy Agency,” http://www.worldenergyoutlook.org/media/weowebsite/ 2013, Paris, France, [Online: accessed July 21, 2015].

[2] B. Hastings, “Ten years of operating experience with a remote controlled water heater load management system at detroit edison,” IEEE Trans. on Power Apparatus and Syst., no. 4, pp. 1437–1441, 1980.

[3] K. Vanthournout, R. D’hulst, D. Geysen, and G. Jacobs, “A smart domestic hot water buffer,” IEEE Trans. on Smart Grid, vol. 3, no. 4, pp. 2121–2127, Dec. 2012.

[4] U.S. Department of Energy, “Energy cost calculator for electric and gas water heaters,” http://energy.gov/eere/femp/ energy-cost-calculator-electric-and-gas-water-heaters-0{#}output, [Online: accessed November 10, 2015].

[5] R. Diao, S. Lu, M. Elizondo, E. Mayhorn, Y. Zhang, and N. Samaan, “Electric water heater modeling and control strategies for demand response,” in Proc. 2012 IEEE Power and Energy Society General Meeting,, pp. 1–8.

[6] S. Iacovella, K. Lemkens, F. Geth, P. Vingerhoets, G. Deconinck, R. D’Hulst, and K. Vanthournout, “Distributed voltage control mechanism in low-voltage distribution grid field test,” in Proc. 4th IEEE PES Innov. Smart Grid Technol. Conf. (ISGT Europe), Oct 2013, pp. 1–5.

[7] S. Koch, J. L. Mathieu, and D. S. Callaway, “Modeling and control of aggregated heterogeneous thermostatically controlled loads for ancillary services,” in Proc. 17th IEEE Power Sys. Comput. Conf. (PSCC), Stockholm, Sweden, Aug. 2011, pp. 1–7.

[8] J. Mathieu and D. Callaway, “State estimation and control of hetero- geneous thermostatically controlled loads for load following,” in Proc. 45th Hawaii Int. Conf. on System Science (HICSS), Maui, HI, Jan. 2012, pp. 2002–2011.

[9] F. Sossan, A. M. Kosek, S. Martinenas, M. Marinelli, and H. Bindner, “Scheduling of domestic water heater power demand for maximizing PV self-consumption using model predictive control,” in Proc. 4th IEEE PES Innov. Smart Grid Technol. Conf. (ISGT Europe), Oct 2013, pp. 1–5.

[10] E. F. Camacho and C. Bordons, Model Predictive Control, 2nd ed. London, UK: Springer London, 2004.

[11] J. Cigler, D. Gyalistras, J. ˇSirok`y, V. Tiet, and L. Ferkl, “Beyond theory: the challenge of implementing model predictive control in buildings,” in Proc. 11th REHVA World Congress (CLIMA), Czech Republic, Prague, 2013.

[12] Y. Ma, “Model predictive control for energy efficient buildings,” Ph.D. dissertation, University of California Berkeley, Mechanical Engineering, Berkeley, CA, 2012.

[13] M. Maasoumy, M. Razmara, M. Shahbakhti, and A. Sangiovanni Vincen- telli, “Selecting building predictive control based on model uncertainty,” in Proc. American Control Conference (ACC), Portland, OR, June 2014, pp. 404–411.

[14] R. S. Sutton and A. G. Barto, Reinforcement Learning: An Introduction. Cambridge, MA: MIT Press, 1998.

[15] D. Ernst, M. Glavic, F. Capitanescu, and L. Wehenkel, “Reinforcement learning versus model predictive control: a comparison on a power system problem,” IEEE Trans. Syst., Man, Cybern., Syst., vol. 39, no. 2, pp. 517–529, 2009.

[16] S. Lange and M. Riedmiller, “Deep auto-encoder neural networks in reinforcement learning,” in Proc. IEEE 2010 Int. Joint Conf. on Neural Networks (IJCNN), Barcelona, Spain, July 2010, pp. 1–8.

[17] E. C. Kara, M. Berges, B. Krogh, and S. Kar, “Using smart devices for system-level management and control in the smart grid: A reinforcement learning framework,” in Proc. 3rd IEEE Int. Conf. on Smart Grid Commun. (SmartGridComm), Tainan, Taiwan, Nov. 2012, pp. 85–90.

[18] G. P. Henze and J. Schoenmann, “Evaluation of reinforcement learning control for thermal energy storage systems,” HVAC&R Research, vol. 9, no. 3, pp. 259–275, 2003.

[19] Z. Wen, D. O’Neill, and H. Maei, “Optimal demand response using device-based reinforcement learning,” IEEE Trans. on Smart Grid, vol. 6, no. 5, pp. 2312–2324, Sept 2015.

[20] M. Gonz´alez, R. Luis Briones, and G. Andersson, “Optimal bidding of plug-in electric vehicles in a market-based control setup,” in Proc. 18th IEEE Power Sys. Comput. Conf. (PSCC), Wroclaw, Poland, 2014, pp. 1–7.

[21] D. Ernst, P. Geurts, and L. Wehenkel, “Tree-based batch mode reinforce- ment learning,” Journal of Machine Learning Research, pp. 503–556, 2005.

[22] M. Riedmiller, “Neural fitted Q-iteration–first experiences with a data efficient neural reinforcement learning method,” in Proc. 16th European Conference on Machine Learning (ECML), vol. 3720. Porto, Portugal: Springer, Oct. 2005, p. 317.

[23] V. Mnih, K. Kavukcuoglu, D. Silver, A. A. Rusu, J. Veness, M. G. Bellemare, A. Graves, M. Riedmiller, A. K. Fidjeland, G. Ostrovski et al., “Human-level control through deep reinforcement learning,” Nature, vol. 518, no. 7540, pp. 529–533, 2015.

[24] M. Riedmiller, T. Gabel, R. Hafner, and S. Lange, “Reinforcement learning for robot soccer,” Autonomous Robots, vol. 27, no. 1, pp. 55–73, 2009.

[25] R. Fonteneau, L. Wehenkel, and D. Ernst, “Variable selection for dynamic treatment regimes: a reinforcement learning approach,” in Proc. European Workshop on Reinforcement Learning (EWRL), Villeneuve d’Ascq, France, 2008.

[26] S. Adam, L. Busoniu, and R. Babuˇska, “Experience replay for real- time reinforcement learning control,” IEEE Trans. on Syst., Man, and Cybern., Part C: Applications and Reviews, vol. 42, no. 2, pp. 201–212, 2012.

[27] F. Ruelens, B. Claessens, S. Vandael, S. Iacovella, P. Vingerhoets, and R. Belmans, “Demand response of a heterogeneous cluster of electric water heaters using batch reinforcement learning,” in Proc. 18th IEEE Power Sys. Comput. Conf. (PSCC), Wrocław, Poland, Aug. 2014, pp. 1–8.

[28] F. Ruelens, B. Claessens, S. Vandael, B. De Schutter, R. Babuˇska, and R. Belmans, “Residential Demand Response Applications Using Batch Reinforcement Learning,” Submitted to IEEE Trans. on Smart Grid (http://arxiv.org/pdf/1504.02125.pdf), Apr. 2015.

[29] M. Deisenroth and C. E. Rasmussen, “PILCO: A model-based and data-efficient approach to policy search,” in Proceedings of the 28th International Conference on machine learning (ICML-11), 2011, pp. 465–472.

[30] T. Lampe and M. Riedmiller, “Approximate model-assisted neural fitted q-iteration,” in 2014 International Joint Conference on Neural Networks (IJCNN), July 2014, pp. 2698–2704.

[31] G. T. Costanzo, S. Iacovella, F. Ruelens, T. Leurs, and B. Claessens, “Experimental analysis of data-driven control for a building heating system,” CoRR, vol. abs/1507.03638, 2015. [Online]. Available: http://arxiv.org/abs/1507.03638

[32] R. Bellman, Dynamic Programming. New York, NY: Dover Publications, 1957.

[33] W. Curran, T. Brys, M. Taylor, and W. Smart, “Using PCA to efficiently represent state spaces,” in The 12th European Workshop on Reinforcement Learning (EWRL 2015), Lille, France, 2015.

[34] D. Bertsekas and J. Tsitsiklis, Neuro-Dynamic Programming. Nashua, NH: Athena Scientific, 1996.

[35] M. Scholz and R. Vig´ario, “Nonlinear PCA: a new hierarchical ap- proach.” in ESANN, 2002, pp. 439–444.

[36] L. P. Kaelbling, M. L. Littman, and A. W. Moore, “Reinforcement learning: A survey,” Journal of Artificial Intelligence Research, pp. 237– 285, 1996.

[37] U. Jordan and K. Vajen, “Realistic domestic hot-water profiles in different time scales: Report for the international energy agency, solar heating and cooling task (IEA-SHC),” Universit¨at Marburg, Marburg, Germany, Tech. Rep., 2001.

[38] “Belpex - Belgian power exchange,” http://www.belpex.be/, [Online: accessed March 21, 2015].

[39] “Elia - Belgian power exchange,” http://www.belpex.be/, [Online: accessed March 21, 2015].

[40] B. Dupont, P. Vingerhoets, P. Tant, K. Vanthournout, W. Cardinaels, T. De Rybel, E. Peeters, and R. Belmans, “LINEAR breakthrough project: Large-scale implementation of smart grid technologies in distribution grids,” in Proc. 3rd IEEE PES Innov. Smart Grid Technol. Conf. (ISGT Europe), Berlin, Germany, Oct. 2012, pp. 1–8.

[41] F. Pedregosa, G. Varoquaux, A. Gramfort, V. Michel, B. Thirion, O. Grisel, M. Blondel, P. Prettenhofer, R. Weiss, V. Dubourg et al., “Scikit-learn: Machine learning in Python,” The Journal of Machine Learning Research, vol. 12, pp. 2825–2830, 2011.

[42] R. Fonteneau, S. A. Murphy, L. Wehenkel, and D. Ernst, “Batch mode reinforcement learning based on the synthesis of artificial trajectories,” Annals of Operations Research, vol. 208, no. 1, pp. 383–416, 2013.

[43] M. Devaine, P. Gaillard, Y. Goude, and G. Stoltz, “Forecasting electricity consumption by aggregating specialized experts,” Machine Learning, vol. 90, no. 2, pp. 231–260, 2013.

designed for accessibility and to further open science