1. Introduction
Residential customers are faced with diverse pricing schemes implemented by utilities to stimulate demand response [
1,
2,
3], and are hence faced with the subsequent problem on how to operate the household appliances to minimize the electricity payment under the premise of human comfort. The growing integration of household generation, mainly intermittent renewable generation, makes this problem even more complex. By integrating advanced automatic control and information and communication technologies, home energy management systems (HEMS) make it possible for residential customers to optimally manage their household appliances, i.e., to shift and curtail household loads to properly respond to the various pricing schemes [
4]. Under this background, HEMS have been attracting increasing attention in recent years.
There are many studies that model the household load scheduling problem as a deterministic optimization problem [
5,
6]. In [
7], a detailed home energy management system structure is developed to determine the optimal day-ahead appliances scheduling under hourly pricing and peak power-limiting. In [
8], in order to distinguish the energy consumption modes and corresponding status of different appliances, home appliances are assigned dynamic priorities, and a real-time household load priority scheduling algorithm, based on the allocated priority and renewable source availability prediction, is proposed to minimize the cost of energy consumption with customers’ comfort constraints. Aiming at minimizing the electricity bill and customers’ dissatisfaction at the same time, Soares et al. [
9] propose a multi-objective genetic algorithm to optimize the time allocation of domestic loads, and similar to the Non-dominated Sorting Genetic Algorithm II, some changes have been introduced in the proposed algorithm to adapt the physical characteristics of the load scheduling problem and improve the usability of results.
However, in these studies, many parameters, such as the time-varying prices, ambient temperature, customer behaviors, etc., are assumed to be obtained through day/hour/minute-ahead forecasting, which contains inevitable uncertainties that will undermine the feasibility and optimality of the load schedules obtained. In order to tackle this issue, many researchers have conducted a large number of works and made very large achievements. In general, these works are divided into the following two methods. One method is to enhance the accuracy of forecast parameters by improving the forecasting methods. In [
10], a new hybrid model that combines two well-known methods for short-term power forecasting of a grid-connected photovoltaic plant is introduced. And hourly forecasts of photovoltaic power output show a quite good accuracy and efficiency of the developed hybrid model. In [
11], Weron collects a variety of methods and ideas about electricity price forecasting. Although these methods have different strengths and weaknesses, it is proved that these methods have obtained varying degrees of success and forecasting precision. Another method is to make ultimate schedules to withstand the uncertainty through various kinds of uncertain optimizations. In [
12], a stochastic scheduling technique is used to handle the uncertainties contained in the energy consumption and runtime of household appliances. In [
13], typical uncertain parameters in the day-ahead temperature scheduling for air-conditioning are modeled by membership functions with fuzzy set theory. In [
14], normally distributed random variables are used to describe the uncertain parameters, and a chance constrained optimization model is then formulated to accommodate the uncertainties.
Although the aforementioned papers make great success to deal with these uncertain parameters in HEMS, there exists some limitations and shortcomings. Due to the complexity of impacting sources in the forecast model and randomness of human behaviors, some forecast parameters, especially like the ambient parameters and customer behavioral parameters, are still not accurate enough to represent the real environment. Meanwhile, most forecast methods need a large amount of historical data which is difficult or expensive for household customers to obtain. Especially, the stochastic programming accounts for a large part of the uncertain optimization methods. Since the characteristics of heuristic algorithms are widely used in the stochastic programming, it is inevitable that the results will be trapped in local optima and make the iteration rather time-consuming.
Differing form the above studies, this paper applies a robust optimization method to tackle the uncertainties in household load scheduling. In 1973, Soyster firstly gained the solutions under the worst situation of uncertain parameters through linear robust optimization [
15]. To decrease the conservatism of linear robust optimization, Ben-Tal et al. [
16] proposed the adjustable robust optimization counterpart. Bertsimas and Sim [
17,
18] extended the robust optimization by introducing an adjustable factor to reflect the different choices from the uncertain parameters set. The robust optimization proposed by Bertsimas has some advantages in handling the uncertainties in parameters: through modeling the uncertain parameters varying in a given uncertainty set, it avoids complex random variables which are subjected to the probability distribution or fuzzy membership function, and it maintains the linearity of the problem so that the global optimal solution can be obtained by a linear programming solver. Additionally, based on setting different values of the factor, it can flexibly control the conservative levels of the results. Based on these, quite a few works [
19,
20,
21] pay significant attention to robust optimization in HEMS. Thus, this robust method is applied in HEMS to solve the scheduling problem with uncertainties.
In this paper, a home energy management system equipped with photovoltaic generation and an energy storage device is established. A household load scheduling with uncertain parameters is formulated to minimize the electricity bill of customers and operates under the comfort constraints. Specifically, the main contributions of this paper can be summarized as follows:
Two typical examples of uncertain parameters, outdoor temperature and hot water demand, are modelled as uncertainty sets, based on which household load scheduling problem with uncertainties is formulated. For researching the uncertainties, the uncertain parameters are presented in the form of interval numbers.
A robust optimization method is applied to deal with the uncertainties in the comfort constraints. The robust counterpart transformation is the key component of robust optimization. Via deducing the robust counterpart, the original scheduling problem is transformed into a mixed integer linear programming problem, of which the global optimum can be found by mature tools, such as CPLEX. The proposed method avoids the time-consuming iterations in many other uncertain optimization methods.
An uncertainty analysis that quantifies the violation degree of the comfort constraints is designed and conducted to test the proposed method. The results show that the complete robust schedules is able to guarantee that the comfort constraints will not be violated. Moreover, schedules with different robust levels can be obtained to make a trade-off between the comfort violation and electricity payment.
This paper is organized as follows:
Section 2 introduces the mathematical model of the home energy management system, including the uncertain parameters, objective function, and constraints;
Section 3 introduces the robust optimization method and completes the derivation of the robust counterpart; and, finally, the simulation results are presented and analyzed in
Section 4 and
Section 5.
2. Mathematical Model
2.1. Uncertain Parameters
In household load scheduling, the outdoor temperature and hot water demand for the next day are representative uncertain parameters that are closely affected by the randomness in the environment and human behaviors and cannot be predicted precisely in practice. Thus, the deviation from the forecast value will lead to the ineffectiveness of the original optimal schedules especially for the loads associated with the uncertain parameters.
To formulate the uncertainties, it is assumed that the uncertain parameters vary in given intervals (i.e., the uncertainty sets) of which the boundaries are set around the predicted values. In practice, the uncertain parameters may take any value in the uncertainty set. Note that no complex probability distribution is assumed within the uncertain parameters.
According to the analysis in [
13] and [
22], the actual outdoor temperature will fluctuate around the forecast values. Thus, the uncertainty set for the outdoor temperature can be described as follows:
where
represents the outdoor temperature (which is uncertain during the scheduling stage) at the
ith time step of the next day.
and
are the forecast value and the maximum deviation value of the outdoor temperature considered at the
ith time step of the next day. Generally, the forecast values are obtained through the weather forecast issued by the local meteorological department. The maximum deviation considered can be confidently summarized by using sufficient historical data. If no sufficient data is available, it can be set as a certain percentage of the forecast value as an approximate estimation.
According to [
23], the uncertainty of hot water demand comes from random human behavior, such as an extra shower that will lead to extra usage of hot water. Therefore, the uncertainty set for the hot water demand can be described as follows:
where
represents the uncertain hot water demand at the
ith time step of the next day;
and
are the forecast value and maximum extra usage of the hot water demand at the
ith time step of the next day.
To respond to flexible electricity prices, the household load scheduling problem is formulated into a mathematical programming problem, of which the objective is to minimize the next-day electricity bill for a customer and the constraints are the operational and comfort limits of the appliances. The solution of the programming problem are the operational schedules for every appliances. The formulation of the problem is detailed below.
2.2. Objective Function
Considering the fact that customers with local generation may feed electricity back to the bulk grid, the objective function is formulated to minimize the electricity bill of the next day, which consists of the payment from buying electricity and the revenue from selling electricity. The objective function is expressed as:
In Equations (3) and (4), N is the total number of the time steps for the next day; is the time length of a single step; A is the set of all household loads; is the power consumption of the load at the ith time step; and represent the active power output of the energy storage device and photovoltaic generation respectively; and and denote the price of buying and selling electricity at the ith time step.
2.3. Constraints
The constraints include the operational limits of household loads and the comfort requirement of customers. For a residential home, the number of appliances can reach about 15 or more. Several appliances have similar operational characteristics so they are classified into the same category. In this paper, the household loads are classified into four main categories: the loads with uncertain parameters, uncontrollable loads, uninterruptible loads, and interruptible loads. Additionally, the energy storage device can be regarded as one special load in which the active power is negative while discharging. Furthermore, because the uncontrollable loads cannot be scheduled, it is directly modeled as a fixed load curve, which will be used in the objective function to calculate the electricity bill.
2.3.1. Loads with Uncertain Parameters
The loads with uncertain parameters widely exist in houses, such as the air conditioner (
AC), water heater (WH), refrigerator, etc. As stated in
Section 2.1, this paper takes the
AC and WH as the typical examples to be studied, which are associated the uncertainties in the outdoor temperature and hot water use.
As for an
AC, due to the heat exchange between the indoor and outdoor environment, the indoor temperature is closely related to the uncertain outdoor temperature, thus being an uncertain variable as well. This relationship can be described as the following Equation (5), which is derived based on the recursive formula presented in [
24]:
where
and
are the uncertain indoor and outdoor temperatures at the
ith time step, respectively;
represents the working status of the
AC (0—off status, 1—heating status and −1—cooling status); the constants
R and
C are the equivalent heat resistance and capacity;
represents the rated power of the
AC. It should be noted that
is the initial indoor temperature, which is given as an initial condition before scheduling.
The main constraint associated with the
AC is the thermal comfort constraint. That is, the uncertain indoor temperature has to be maintained in a preset range at each time step:
where
and
are the lower and upper boundaries of the preset comfort range. Equation (7) shows the relationship among the power consumption, the working status and the rated power of the
AC.
As for a WH, the decrease of the water temperature in tank results from the heat exchange with the environment and the hot water usage. Due to the excellent insulation, the heat loss in the process of the heat exchange can be ignored [
25]. Then the main cause of the temperature dropping is the hot water usage and the followed supplementary cold water. Based on the energy conservation law, the water temperature in tank can be calculated as follows:
where
is the uncertain water temperature in tank at the
ith time step;
is uncertain hot water usage;
is the heating power; the constant
is the transfer coefficient between J and kWh, taking the value of 3.6 × 10
6; the constant C is the specific heat capacity of water, being 4.2 × 10
3;
M is the mass of the water in full tank; and
is the initial water temperature in tank, which is also given before scheduling.
Similar to
AC, the constraint of the WH is also about thermal comfort:
where
and
are the lower and upper thresholds which are preset by customers. Equation (10) shows that the heating power of the WH is the decision variable to be optimized.
2.3.2. Interruptible Loads
Interruptible loads (
IL) are allowed to be stopped during their working time so that load shifting can be achieved. Electric vehicles (EVs) and clothes washers (CWs) are two examples of this category of load. The constraints for interruptible loads can be expressed as follows:
where
is the working status of interruptible loads (0—off, 1—on);
represents the allowable working period;
is the required working time to; and
is the rated power.
2.3.3. Uninterruptible Loads
Uninterruptible loads (
UIL) have very similar operational characteristic to interruptible loads and are also allowed to be shifted. However, uninterruptible loads, like a clothes dryer (CD) or a dishwasher (DW), cannot be interrupted once it starts working. Therefore, additional constraints should be proposed [
20].
Similar to interruptible loads, is the working status of uninterruptible loads (0—off, 1—on); represents the allowable working period; is the required working time; represents the rated power. Among the constraints, Equation (16) is the additional one for describing the uninterruptible feature.
2.3.4. Energy Storage Device
Household energy storage devices are usually batteries (e.g., lead-acid, Li-ion). Batteries charge and discharge to match the output of the local photovoltaic generation. The operational limits during the charging and discharging process are as follows [
26]:
where
and
represent the charging and discharging power;
and
represent the charging and discharging status (
= 1 and
= 0 indicate that the energy storage device is charging;
= 0 and
= 1 indicate that the energy storage device is discharging);
and
are the efficiency of charging and discharging;
and
are maximum power of charging and discharging;
represents self-discharge rate; Equations (18)–(22) constrain that the battery should not overcharge and overdischarge; and Equation (23) ensures that the battery capacity at the final time step of the day is no less than that of the initial time step at the beginning of the day.
3. Robust Optimization for Household Load Scheduling
3.1. Robust Optimization
In the view of mathematics, the above scheduling problem is a mixed integer programming (MIP) problem with uncertainties. A nominal mixed integer programming is given as Equation (24).
However, the existence of the uncertain parameters (some elements in of the Equation (24)) hinder the available methods like branch and bound method from solving the problem. Therefore, the robust optimization technique proposed by Bertsimas et al. is applied to make the uncertain programming problem solvable and to make the solutions robust.
In the robust optimization, if the element
in the matrix
A is uncertain, it should be replaced by
that takes value in
of which
and
are the mean value and range of the coefficient (the elements in
can be seen as
through some transformations [
12,
13]). A new parameter is then defined as
, so
is a bounded random variable (but with unknown distribution) which is within [−1,1].
In order to deal with the uncertain but bounded parameters, Bertsimas et al. also introduced a number that takes value in the interval , where for each constraint. The purpose of intruding this parameter is to adjust the robustness against the level of conservatism of the solution. Specifically, the number of the coefficients that are considered as uncertain equals to . For the former uncertain coefficients, the coefficients are considered to vary in the interval ; for the final uncertain coefficient, it is considered to vary in . The larger is, the more uncertainties will be considered during the optimization, and the more robust of the solution will be. When equals to , it means that the solution can satisfy all the situations regardless of what value the uncertain parameter takes.
3.2. Robust Counterpart Transformation
Based on the above definition and method, the original programming problem Equation (24) can be transformed into the following robust counterpart [
18]:
3.2.1. Robust Counterpart Transformation of the Uncertain Outdoor Temperature
When it comes to household load scheduling, the constraints with uncertain parameters need to be transformed into the robust counterpart. Firstly, we should transform the constraints of
AC into the matrix form. From the above paper, only one parameter is uncertain in each constraint and it does not appear in matrix
A. Thus, introducing a new variable
makes this uncertain parameter classified into matrix
A. To simplify the proof, set
. Then the matrix form for constraints is as follows:
Since
, then there exists
. The details are shown in Equation (28):
Based on the robust optimization method proposed by Bertsimas,
[0,1] and
. Thus, the original robust counterparts from Equations (25) and (26) can be transformed into the following formulas:
Then, going on the simply conversion of the max and min functions, the transformed formulas are:
Bringing the detailed elements in Equation (28) back into Equations (32) and (33), the final linear robust counterparts can be obtained:
where
is the parameter that control the solution robustness for the air conditioner, which takes a value in the range [0,1].
3.2.2. Robust Counterpart Transformation of Uncertain Water Demand
Similarly for the water heater, the linear robust counterparts for the constraints in Equations (8) and (9) with the uncertain water demand should be given. However, unlike the air conditioner, strong duality should be used to deduce the robust counterpart considering the multiple uncertainties (i.e., outdoor temperature and hot water use) involved. Except the uncertain parameters in matrix
A, there is another one uncertain parameter which do not appear in matrix
A in each constraint. Thus, introducing the new variable
makes this uncertain parameter classified into matrix
A. To simplify the proof, set
. Then the matrix form for constraints is as follows:
The functions with uncertain variables in the coefficient matrix are all monotonically decreasing function, and
. Thus, the elements in the coefficient matrix are decreasing, along with the increase of
. Thus,
and the details are shown in Equations (37) and (38):
Based on the robust optimization method proposed by Bertsimas, we can obtain
,
and
[0,
i]. Thus, the original robust counterparts (25–26) can be transformed into the following formulas:
Obviously, there is
in Equation (39). While Equation (40) cannot be directly solved which should apply the strong duality to transform. The final form of the linear robust counterparts is directly given as follows:
where
is the parameter that controls the solution robustness for the water heater, which takes a value in the range [0,
i].
3.3. Solution to the Robust Counterpart
The household load scheduling with uncertain parameters is finally transformed into a mixed integer linear programming (MILP) problem. In this MILP problem, the objective function is still the one before transformation, but the constraints consist of the new constraints as presented as Equations (34)–(35), (42)–(47), and the original constraints as presented by Equations (11)–(23). Many existing methods, tools, and commercial software are available to obtain the ultimate solution. Among the abundant methods, the heuristic-based evolutionary algorithms, like particle swarm optimization (PSO), the model predictive method, and the commercial software CPLEX are some representative ones.
However, because of numerous appliances, the MILP proposed in this paper contains too many decision variables and constraints which may cause the curse of dimensionality. In this situation, the heuristic algorithms are not suitable for solving. Given the powerful calculating capacity and various solvers, this paper adopts the commercial software CPLEX (Version 12.6.3.0, IBM, New York, NY, USA, 2016) to gain the optimal schedules. CPLEX is a commercial mathematical programming solver produced by IBM Corporation in the United States. It can solve various mathematical programming problems, including linear programming (LP), quadratic programming (QP), and mixed integer linear programming (MILP) stably and efficiently.
5. Conclusions
In this paper, the household load scheduling problem with uncertain parameters is studied. First of all, the scheduling problem is formulated as a mathematical programming problem which aims at minimizing the electricity bill under various constraints. Unlike the deterministic household load scheduling, the uncertain parameters, such as the outdoor temperature and extra hot water usage, are focused and modeled into uncertainty sets. After that, robust optimization is applied in the scheduling to deal with the uncertainties.
In the simulation results, the impact of uncertainties on the feasibility of the schedules is evaluated and analyzed first. The result indicates that the uncertain parameters may cause the infeasibility of the schedules derived from deterministic optimization. After that, the complete robust schedules are proposed and verified, which are capable of withstanding all the uncertainties, leading to no comfort violation. Finally, the economy and comfort of schedules with different robust levels are compared quantitatively. The proposed robust optimization method allows customers to make a trade-off between the economy and comfort, by choosing the schedules with different robust levels.