An Inventive Method for Eco-Efficient Operation of Home Energy Management Systems

A demand response (DR) based home energy management systems (HEMS) synergies with renewable energy sources (RESs) and energy storage systems (ESSs). In this work, a three-step simulation based posteriori method is proposed to develop a scheme for eco-efficient operation of HEMS. The proposed method provides the trade-off between the net cost of energy (CEnet) and the time-based discomfort (TBD) due to shifting of home appliances (HAs). At step-1, primary trade-offs for CEnet, TBD and minimal emissions TEMiss are generated through a heuristic method. This method takes into account photovoltaic availability, the state of charge, the related rates for the storage system, mixed shifting of HAs, inclining block rates, the sharing-based parallel operation of power sources, and selling of the renewable energy to the utility. The search has been driven through multi-objective genetic algorithm and Pareto based optimization. A filtration mechanism (based on the trends exhibited by TEMiss in consideration of CEnet and TBD) is devised to harness the trade-offs with minimal emissions. At step-2, a constraint filter based on the average value of TEMiss is used to filter out the trade-offs with extremely high values of TEMiss. At step-3, another constraint filter (made up of an average surface fit for TEMiss) is applied to screen out the trade-offs with marginally high values of TEMiss. The surface fit is developed using polynomial models for regression based on the least sum of squared errors. The selected solutions are classified for critical trade-off analysis to enable the consumer choice for the best options. Furthermore, simulations validate our proposed method in terms of aforementioned objectives.


Introduction
From the previous decades, the energy requirement has grown to a critical level; however, the generation units have not been maintained at a sufficient rate to manage this increasing demand.The balance between demand and generation is a vital requirement for stable power system operation.The problem to maintain this balance has conventionally been addressed in the past; utilities have upgraded their centralized generation units and transmission capabilities through some supply side management methodologies.However, during the previous decade, demand-side management (DSM) has become a substituent scheme to manage the increasing requirement of energy which focuses on the consumer side.The home energy management system (HEMS) is used to implement DSM in a home.Major approaches for HEMS operation include price-based demand response (DR), and DR synergized with renewable energy sources (RESs) and energy storage systems (ESSs) optimal dispatch (DRSREOD) [1].The DR-based HEMS operation schedules the consumer's loads by shifting them towards the off-peak periods.Such scheduling benefits the consumer with a minimized cost of energy (CE) based on the acceptable value of time-based discomfort (TBD) [2,3].The utility, on the other side, is benefited with a reduced cost of generation through a smoother demand profile.The DRSREOD-based HEMS operation schedules the load in coordination with the optimal dispatch of the power grid, renewable energy sources (RESs) and energy storage systems (ESSs).The operation of such HEMS introduces additional benefits by minimizing the electricity cost, minimizing high demands and permanent demands, increasing total cost minimization and empowering the selling of the extra power to the utility [4][5][6][7][8][9][10][11].The aforementioned HEMSs are modeled to optimize the objectives comprising the net CE (CEnet), consumer discomfort/inconvenience, and peak and permanent demands.The abbreviations and nomenclature are given in Tables 1 and 2, respectively.Furthermore, the rampant rise in green house gas (GHG) emissions, the consequent climate changes and the related environmental issues have raised serious concerns over the quality of the life on the earth.In order to mitigate the serious environmental issues, various proposals have been discussed for GHG emissions at the highest international forums to confine them.The Kyoto protocol of United Nations Framework Convention on climate change has been mutually signed by 192 countries all over the world which proposes a reduction in GHG emissions through selling of emission commodities [12].Such a trading sets penalties and quantitative limitations on emissions by polluters that may include utilities, independent microgrid (MGD) operators, and the prosumers having fossil fuel based generation deployed with DRSREODLDG-based HEMSs.
The aforementioned scenario has incentivized utilities to reduce not only the cost of generation of energy; however, also the supply-side emissions making use of RESs installed for DRSREOD-based HEMS.The research on HEMS now seems to focus on reducing the GHG emissions along with the other well-known objectives for CE, TBD, etc.In [13], a scheme for DR-based HEMS is presented.Non-critical house loads are shifted towards low demand periods for minimizing the daily bill of the generation-side and the supply-side emissions.It is validated that implementation of the DR program effectively reduces the cost of generation on the supply-side; however, the emission on this side is reduced only when peak demand is met by high emission fuels based peaking plants.The DRSREOD-based HEMS, on the other hand, through an optimal operation of HEMS devices, can easily be used for reducing the supply-side emissions along with the reductions in the CEnet for the consumer and cost of generation for the utility.In [14], the authors present a scheme for optimal scheduling of shiftable home appliances (SHAs) integrated with the optimal dispatch of an RES and an storage battery (SB).The objectives include reductions in CEnet, temperature based discomfort, peak load, and the GHG emissions.The supply-side emissions are computed using GHG emission factors (EFTs) for the energy mix adopted at different times of the day.The supply-side emissions are reduced through an optimal operation of local RESs and SBs during high emission times.
Furthermore, fossil fuel based DGs are integrated into MGDs to improve the self-healing structure and the flexibilty of the power supply.In [15], an operational scheme is developed for a stand-alone HEMS operation using PSO.The scheme is based on load shifting of SHAs, an optimal dispatch of a wind turbine (WTB), a DG, and an SB.The DG is operated at the rated power in order to improve its efficiency and to reduce emissions.The power from the grid, however, is not included in the modeling.
In [16], an optimal dispatch scheme for a PV unit, a WTB, an ESS, a DG, and the power grid to supply a fixed load profile in an MGD is computed using GA.Constraints for ESS charge/discharge rates, generator start/stop and supply capacity are taken into account.Net emission for the power supplied by the grid and local DGs are computed.be analyzed for the emissions released during the LDG operations.To implement an eco-efficient operation of DRSREODLDG-based HEMS, optimal trade-offs between CEnet, TBD and minimal GHG emissions (TEMiss) need to be computed.This research introduces a method to harness a diversified set of solutions to decision vector Tst and the related trade-offs for CEnet, TBD and minimal TEMiss for an eco-efficient HEMS operation.
The proposed method for an eco-efficient operation of DRSREODLDG-based HEMS is based on a three-step approach.In step-1, a set of primary solutions in terms of Tst and the related trade-offs for CEnet, TBD and minimal TEMiss are generated using Algorithm 1.The algorithm is based on a heuristic derived from our previous studies on HEMS in [11].The proposed heuristic takes into account PV availability, the state of charge, the related rates for the storage system and the similar functionality of the sources.To achieve maximum reduction in CEnet, SHAs are modeled for mixed scheduled (MS) as already validated in [11].This research formulates the trade-off parameters for: CEnet to include the cost of energy purchased from the grid, cost of energy sold to the grid and the cost of energy supplied by the LDG; TEMiss to include the energy supplied by the LDG during LSD hours, EFT based on the calorific value of the fuel, the consumption efficiency of the LDG and the related emission factors for GHGs; and TBD to include the delay in the starting times of delay scheduling (DS) type and advanced completion of the job of advanced scheduling (AS) type for HAs.The trade-off solutions obtained in step-1 are analyzed for TEMiss as related to the trade-offs between CEnet and TBD as shown in Figure 2. The plot in Figure 2 reveals a highly un-even relation between TEMiss and the related parameters for CEnet and TBD.This un-even trend for TEMiss is exploited to screen out/exclude a set of trade-offs with larger values of TEMiss using a constraint filtration mechanism as presented in Algorithm 2 (step-2 and step-3).In step-2, an average value based constraint filter (AVCF) for TEMiss is developed and applied to filter out the trade-offs with extremely high values of TEMiss.In step-3, average surface fits for TEMiss are developed in terms of CEnet and TBD using polynomial based regression.The most suitable polynomial is selected after cross-validation of 25 number of polynomial model that fits on their capabilities in order to reduce TEMiss and TBD, and to maximize the number of diverse trade-offs for CEnet and TBD.The average surface fit based constraint filter (ASCF) with the selected polynomial formulation is applied to screen out the trade-offs with even marginally higher values of TEMiss.
The solutions for an eco-efficient HEMS operation are thus achieved including diversified trade-offs for CEnet, TBD and minimal TEMiss.Followings are the novelties of this research:

•
An innovative method is proposed to harness diversified trade-offs between CEnet, TBD and minimal TEMiss for an eco-efficient operation of DRSREODLDG-based HEMS.

•
Trade-offs for such HEMS have rarely been computed by combining a multi-objective genetic algorithm or Pareto optimization (MOGA/PO) based heuristic and regression based constraint filtration.The polynomial model fit for regression is based on its capabilities to reduce the trade-off parameters for eco-efficient HEMS operation.

•
Most of the authors use the weighted sum method (WSMD) to handle multi-objectivity for similar problems.This research presents a diverse set of trade-offs that are critically analyzed to enable the consumer choosing the best option.

•
Trends exhibited by the trade-off parameters are analyzed based on vital factors affecting these parameters, e.g., loss of unused energy from the PV unit.

•
The proposed method validated to minimize the emissions from a local LDG for a DRSREODLDG-based HEMS; however, it is easily extendable to reduce the supply-side emissions as well.
The organization of the paper is as follows: Section 2 describes the related work whereas the system model is elaborated in Section 3. The problem formulation for eco-efficient operation of DRSREODLDG-based HEMS and the techniques used to solve this problem are presented in Section 4. The proposed Algorithm 1 to generate primary trade-offs for optimal HEMS operations, and Algorithm 2 to harness eco-efficient trade-offs through a constraint filtration mechanism are presented in Section 5.In Section 6, simulations are presented to demonstrate the validity of Algorithm 1 to generate schemes for DRSREODLDG-based HEMS operation in terms of Tst and the primary trade-offs between CEnet, TBD and TEMiss.The trends exhibited by the primary trade-offs are analyzed in detail and the bases for the selection of constraint filters including AVCF and ASCF are validated.Further simulations are presented to demonstrate the validity of Algorithm 2 to harness eco-efficient trade-off solutions using the optimal constraint filters.Conclusions and future work are discussed in Section 7.

Related Work
With the installation of smart grid technologies enabling DSM, a widespread deployment of DRand DRSREOD-based HEMSs has been carried out throughout the world in the past few years [17,18].The report in [17] has given an overview and boost to the RESs by forming policies among various contries.In [18], the current phase of Paris agreement has focused on developing a global approach, which limits the GHG emissions of all countries.In recent years, authors have presented various models and methods for the optimal operation of such systems [2][3][4][5][6][7][8].The objectives for optimal HEMS operation include minimizing CE, TBD, peak-to-average ratio (PAR) and peak/permanent demands [19][20][21].In [19], the authors have used different priorities to derive user comfort.Khan et al. [20] have used three different appliances to minimize CE, TBD, and PAR.Meta-heuristics approaches including optimal stopping rule are used as optimization schemes.Similarly, the authors in [21] have applied meta-heuristics approaches along with a hybrid approach to minimize CE.Furthermore, utilities owning energy deficient power networks in developing countries are subjecting their users for LSD to maintain energy demand and supply.In such power networks, consumers deploy a LSD-compensating DG in DRSREOD-based HEMS for ensuring the reliable distribution of the energy [11].The aforementioned objectives for optimal HEMS operation have been achieved using optimization techniques like linear programming (LP), mixed integer LP (MILP), advanced heuristics, etc.
Additionally, the issue regarding serious environmental concerns over the use of fossil fuels has been raised at international forums consistently in the past few decades.Recently, worldwide consensus has been reached to reduce the GHG emissions by selling them as commodities [12].Such trading sets quantitative limitations on the emissions made by polluters that may include utilities, independent MGD operators and the prosumers having local fossil fuel based generations.The present scenario based on the polluter pays principle has incentivized utilities to reduce not only the generation cost; however, the supply-side emissions as well while making use of the RESs installed for DRSREOD-based HEMSs [14,22,23].Furthermore, MGD operators having RESs, ESSs and DGs also include TEMiss as an objective in the optimal dispatch scheme for their systems [15,16,24].Furthermore, in energy-deficient power networks, DRSREODLDG-based HEMSs having LSD-compensating DGs are used to ensure an uninterrupted supply of power during LSD hours [11].The operation of LDG in such HEMSs, however, does accompany the release of emissions, which needs to be minimized.
The related work includes the recent research on models and methods to achieve important objectives for DR and DRSREOD-based HEMSs including reductions in TEMiss (supply-side), CEnet, and TBD; for MGDs including reductions in TEMiss and CEnet ; and for DRSREODLDG-based HEMS including reductions in TEMiss (local), CEnet and TBD.The recent research related to the proposed method for an eco-efficient operation of DRSREODLDG-based HEMS is summarized in Tables 3 and 4.

Emissions Reduction Using DR-Based HEMSs
Most of the research on DR-based HEMS has focused on objectives like CE, PAR, peak load and discomfort [2,3].Such systems have limited capabilities to play a role in the reduction of GHG emissions.In [13], a scheme for DR based HEMS is presented.Non-critical house loads are shifted towards off-peak hours to minimize the daily cost of generation and emissions for the supply-side.It is validated that implementation of DR program effectively reduces the cost of generation on the supply-side; however, the emissions on this side are reduced only when peak demand is met by peaking plants based on high emission fuels like coal, diesel, etc.

Emissions Reduction Using DRSREOD-Based HEMSs
Most of the models for DRSREOD-based HEMS presented in the recent past are based on optimal scheduling of SHAs integrated with the optimal dispatch of RESs and ESSs.HEMS problems for these models have been solved to reduce CE and discomfort for the consumer, and to minimize peak load/PAR and cost of generation for the utility [4][5][6][7][8].Recently, in the context of worldwide concerns over GHG emissions, authors have focused on the reduction in emission as an objective for DRSREOD-based HEMS.In [14], authors proposed a scheme for optimal scheduling of SHAs integrated with the optimal dispatch of RES, SB, and the utility.Major goals include reductions in CEnet, temperature based discomfort, peak load, and the supply-side emissions.Such emissions are computed using emission coefficients for the energy mix adopted by the utility during various times of the day.An optimal dispatch of local RESs and SBs results in the reduction of net supply-side emissions by supplying the load during high emission hours.MILP has been used to solve the model.In [22], an operating mechanism of major HAs including heating and cooling appliances integrated with the optimal dispatch of PV and SB is presented.The algorithm for real-time HEMS operation is based on user preferences, home occupancy, day ahead emissions and climate forecasts.The objectives for reduction in CE, electric consumption, TEMiss, and the peak demand are formulated.The net cost of emission includes carbon footprint of the customer from the grid electricity usage minus carbon reduction from injecting emission-free electricity from RES.In [23], a prosumer based algorithm is presented to maximize the sum of benefits to the consumer and the utility.The emission trading has been considered as a mean of mitigating this commodity.The utility is profited by reducing his carbon footprints while purchasing energy from locally installed RESs and ESSs during high emission times.The fitness function maximizes the welfare including consumption-based satisfaction and monetary benefits from RESs and ESSs to consumers and benefits of the reduced peak load, generating cost and emissions to the utility.A dynamic selling with dynamic buy-back pricing scheme is also proposed to implement the model.For scalability and user privacy, the problem is solved using Lagrange multipliers.The objectives in all of the above research are combined using the WSMD.

Emissions Reduction in MGDs
In MGDs, RESs and ESSs are integrated with DGs to enhance the quality and the reliability of the power supply.In [15], a solution for DRSREOD-based HEMS operations for a stand-alone home including WTB, DG, and SB is computed using PSO.The local fossil fueled DG is operated at rated power for an improved efficiency and reduced emissions.A separate objective function for emissions; however, is not included.An optimal dispatch for an MGD is computed in [16] using GA.Additional constraints for ESS charge/discharge rates, DG start/stop and supply capacity are considered.Total emission is computed using emission factors for the grid, power supplied from the local DG and the ESS.The model does not include load shifting while computing the dispatch for power sources.A method to compute an optimal dispatch of RESs and DGs for a MGD is presented in [24].The dispatch is based on costs of energy from WTB, PV and DG, EMiss and CE from/to main grid for a fixed load profile.The WTB and the PV are the preferred sources.The SB is discharged based on its SoC if local RESs are not able to meet the demand; else, the load is supplied through the economic dispatch of the DG, fuel cell (FC), SB and the grid.Non-critical loads are disconnected when local sources are insufficient.The DG is operated at rated power to minimize EMiss.DR based load shifting is not included.

Emissions Reduction in DRSREODLDG-Based HEMS
Energy-deficient power supply networks in developing countries are based on the compromises for the consumers to LSD in order to maintain the equilibrium between demand and generation of energy [10,11].While a number of consumers in developing countries are participating in DSM making use of DRSREOD-based HEMSs, LSD-compensating DGs are deployed in such HEMSs for ensuring the uninterrupted supply of electricity.An algorithm for optimum sizing of an LDG for DRSREOD-based HEMS was presented in our recent research [11]; however, such a DG does introduce emissions when operated during LSD hours.Based on the recent scenario for quantitative restrictions on carbon emissions, research on the optimized operation of DRSREODLDG-based HEMS focusing reduction in TEMiss looks pertinent.A simulation-based posteriori method for an eco-efficient operation of DRSREODLDG-based HEMS takes into account the trade-offs between CEnet, TBD, and minimal TEMmiss is proposed.A three-step approach is followed.At step-1, primary trade-off solutions for CEnet, TBD, and TEMmiss are generated using a heuristic proposed for an optimal operation of DRSREODLDG-based HEMS.The heuristic, which uses MOGA/PO to search optimal trade-offs, is detailed in Algorithm 1.At step-2, an AVCF is used to filter out the trade-offs with extremely high values of TEMiss, whereas an ASCF is used to screen out the trade-offs with marginally high values of TEMiss at step-3.The ASCF was developed using advanced regression techniques.The filtration mechanism including AVCF and ASCF used to harness eco-efficient trade-off solutions for CEnet, TBD, and minimal TEMmiss is detailed in Algorithm 2.

System Model
The architecture for DRSREODLDG-based HEMS is shown in Figure 1.The major components of such HEMSs include home appliances, renewable energy sources, an energy storage system, an LSD-compensating DG, a HEMS controller, a local communication network, and a smart meter for bidirectional interation between users and the power grid.The proposed optimal operation for such HEMS are based on DR synergized with the optimal dispatch scheme for RESs, ESSs and an LDG.The operating scheme takes into account the MS of SHAs, their combined corresponding functions of the PV unit, the SB and the utility, and the energy sold to the grid based on the parametric values of power vector from PV(P pv ), vectors of the state of charge (SoC), the maximum charge/discharge rates, and the tariff scheme.PV unit is the preferred source that is responsible for suppling the power to the scheduled appliances.The surplus PV energy is saved into the SB for utilizing the power during peak hours and it is sold to the utility for a monetary benefit.During the LSD hours, if SB is full and there is no energy demand than the excess energy from the PV unit is dissipated in a dummy load [9].The mentioned energy (shown by Pdl) represents a loss of the PV energy that could not be sold due to the unavailability of the main grid.The LDG is used for supplying the load in high demand periods which is contributing similar to PV unit and the SB to prevent the electricity blackouts.The operation of the LDG in such systems ensures an uninterrupted supply of power; however, such operation of the LDG accompanies the release of GHGs emissions as well.The problem for DRSREODLDG-based HEMS operation has been formulated as multi-objective-optimization (MOO) to minimize CEnet, TBD, and TEMiss.
Furthermore, according to [10], 31% and 21% energy is consumed in industrial and residential sectors, respectively.However, in this paper, we consider only residential area for implementation of our proposed scheme.Because our proposed schemes are based on load scheduling from ON-peak to OFF-peak hours, it is not possible in industrial or agriculture sectors to reduce electricity cost via load shifting due to production problems.Therefore, we consider only residential area for implementations.Moreover, there are a lack of sources in developing countries, so the energy management system is a huge opportunity for these countries.
A three-step simulation based posteriori method is proposed to provide trade-off solutions for an eco-efficient operation of DRSREODLDG-based HEMS.The method makes use of Algorithm 1 and Algorithm 2 to harness eco-efficient schemes for HEMS operation in terms of Tst and the related trade-offs for CEnet, TBD, and minimal TEMiss.At step-1, primary trade-off solutions for CEnet, TBD, and TEMiss are generated making use of Algorithm 1. Algorithm 1 is based on a MOGA/PO based heuristic proposed in this work.At step-2, the primary trade-off solutions are passed through an AVCF to filter out the trade-offs with extremely high and above average values of TEMiss.The filtrate is then passed through an ASCF to screen out the trade-offs with even the marginally higher values of TEMiss at step-3.The proposed filtration mechanism comprising AVCF and ASCF is detailed in Algorithm 2. The simulations to validate the method for harnessing the desired trade-offs for eco-efficient operation of DRSREODLDG-based HEMS are presented in Section 6.
Major components of the proposed model for DRSREODLDG-based HEMS are presented below.

Parameters for Scheduling
A scheduling resolution of 10 min.slot has been adopted.To formulate the HEMS operations, total time is sub-divided into 144 slots.While scheduling, each SHA is executed in the specific horizon for a specified number of slots.The proposed model for HEMS operation is dependant on a dynamic electric pricing signal, i.e., an IBR pricing signal, the PV panel, the SB, and the LDG.The control parameters for HEMS components are described in Section 4 on problem formulation.

HAs
Motivated from the literature [21,25], the HAs are classified into non-shiftable home appliances (NSHAs) and SHAs.NSHA, e.g., electric lamps, and fans, are working on required time slots and can not be opted for scheduling.SHAs are assumed to be scheduled towards the low demand hours and the PV harnessing hours for optimized HEMS operation.To achieve a maximized reduction in the cost of energy, shiftable appliances are modeled as AS and DS.Such classification enables more reduction in the cost of energy making use of enhanced flexibility in the appliances shifting and an increased direct usage of the PV energy from the PV unit [11].AS and DS type SHAs with the user prioritize settings and the NSHAs along with their forecasted load, used in the simulation section are described in Tables 5 and 6.The proposed method is generic in nature and is equally applicable to DRSREODLDG-based HEMS based on time-of-use pricing and Day-ahead pricing schemes.A 2-step ToU pricing tariff with respective IBR values and threshold power demand are used in simulations in Section 6.

RESs
Solar irradiation data as measured by the Pakistan Engineering Council in Islamabad has been applied for the simulations to validate the proposed system [11].It is possible to sell the surplus energy produced from the RESs [26].The PV system is used in the simulations and its parameters' configurations are given in Table 7.The energy generation profile of PV unit is displayed in Figure 3.The electricity bill of generations from local RESs has not been included in the model and such installations have been considered as a module of the current system [11].Efficiency of PV panels 15%

ESS
The SB and the inverter are important parts of the proposed DRSREODLDG-based HEMS.These components along with their specifications are given in Table 8.Net loss for the SB is initially supposed up to 20%; otherwise, it is considered for charging.

LDG
The consumers perform optimally sizing the LDGs according to their deficient load as [27].The specifications of the LDG to cope with the LSD are used in the simulation of proposed model as given in Table 9.The emission factor is computed as per Equation ( 8) using the pertinent data given in [28,29].The cost of energy for the LDG is according to the levelized cost of energy for such units given in [30].

Formulating DRSREODLDG-Based HEMS Problem and the Related Optimization Techniques
The contents of this section are inherently divided into two parts: (1) the problem formulation for HEMS optimization to generate optimal schedules of SHAs in terms of Tst and the primary trade-offs for CEnet, TBD and TEMiss along with the proposed techniques.
The problem to generate the primary trade-offs for CEnet, TBD and TEMiss is computed using the following input values: The problems for HEMS are formulated to compute the power requirement for whole of the scheduling horizon.With Tst as the decision vector, HEMS problem for the scheduling vector Pschd_sh is treated as MILP and computed using the following Equations.: where X(b, n) for the b th SHA is computed based on the following terms, The load vector Pload_nsh for NSHAs is added to Pschd_sh to calculate the final scheduled load vector (Pschd) using Equation (2): The problem is solved using a MOGA/PO based heuristic proposed in Algorithm 1 to obtain a decision vector, Tst, which optimizes the outcomes for trade-offs parameters by fulfilling the required constraints.

Objectives for DRSREODLDG-Based HEMS Problems
The main objective for DRSREODLDG-based HEMS optimization is to achieve the optimal trade-offs between the CEnet, TBD and TEMiss.To obtain the above-mentioned trade-offs, the problem for HEMS is formulated using SHA schedules by calculating Pschd paralelly, synergizing the scheduling with RES, ESS, LDG and power grid dispatch for N time intervals over a specified scheduling horizon.The LDG is integrated in the dispatch only for the LSD intervals to supply the load in collaboration with the PV unit and SB.The heuristic presented in Algorithm 1 computes the corresponding vectors for Pgd, Pgn, and Psold that are used to compute the objective functions/trade-off parameters.

Reduction of CEnet
The CEnet is computed using the Equation ( 5): Here, Pgd and PE are the power purchased from the main grid and its price in Cents/kWh.The terms Psold and PE f are the energy sold to the utility by the consumer and its feed-in price in Cents/kWh, respectively.Pgn and PEg are the energy supplied from the LDG and its levelized price in Cents/kWh, respectively.A factor Psold × CEMiss can be excluded from CEnet as a reward for reducing the supply-side emissions through the PV energy sold to the utility.

Minimization of TBD for the Consumer
We consider user discomfort in terms of average waiting time (TBD) of appliances.It simply means how much time a user will wait to switch ON any appliance.Moreover, the maximum average waiting time means maximum user discomfort and vice versa.For instance, if the average waiting time of all appliances is four hours, then user discomfort will increase by 0.4 because the user feels discomfort to wait.Moreover, in the case of unscheduled electricity consumption, users do not wait to turn ON their appliances.Electricity user can use any appliance at any time, so their waiting time is 0 and the user's comfort is maximum (no discomfort).The TBD formulation is given below: where STslot and ENslot indicate the users' time bounds flexibility for representing the starting and ending boundaries for SHAs working.LoT is a vector which is comprised of total length of operation time for eacg SHA for completing its execution.Tst is also a decision vector which consists of starting intervals of all SHAs, whereas k 1 is the number of SHAs designated as DS.
When the corresponding SHA starts its execution at STslot, the TBD(D) uses its initial value as 0, i.e., the start of the execution time is assigned from the users.When Tst(b) is equal to ENslot(b) − LoT(b) + 1, it obtains its maximum value at 1, i.e., late starting time for the SHA results in finishing of the task at the late alloted time ENslot(b).The boundaries for the selection of Tst should be considered feasibly which are calculated with the help of the next Equation: Due to the advanced completion of the jobs of AS type SHAs denoted by TBD(A), the average TBD is computed.It is calculated by taking the average of the normalized advance-completion times of all gadgets.This value is computed using the next equation: whereas, k 2 denotes the number of SHAs for AS type.
When the corresponding SHA completes its execution at ENslot(b), TBD(A) takes its initial value as 0, for example, when Tst(b) + LoT(b)-1 is equal to ENslot(b).When Tst(b) is equal to STslot(b), it gets its maximum value as 1, for example, using the finishing time of the appliance which is calculated with There are muitple appliances and some of them are categorized as AS and others are categorized as AS in MS.For minimizing the average TBD for total k appliances in MS mode, the objective function is defined as below: Minimize(TBD(D) + TBD(A)).
For a scheduling flexibility and better reduction in CEnet, TBD based on MS of SHAs as per Equation (7) has been opted for this model.

Reduction of TEMiss
Equation ( 8) computes the emissions' minimization ratio, i.e., TEMiss, from which the LDG is formulated: Here, EFT is the carbon emission factor (kg/kWh) and Pgn is the vector for the energy supplied by the LDG in kWh during the LSD hours [13,16].A value of 1.6 Lbs./kWh has been used for EFT for the LDG as per the data available in [28,29].

Techniques to Solve DRSREODLDG-Based HEMS Optimization Problems
The HEMS optimization is considered as a combinatorial optimization problem.Multiple HEMS problems are nonlinear, non-convex constrained, multi-dimensional in nature and they have a variety of the solutions available in literature.To solve these problems, both conventional and advanced heuristic optimization techniques are used.

Techniques to Handle Multi-Objectivity in HEMS Optimization
Multiple HEMS optimization problems faced in existing scenarios are multi-objective optimization (MOO); however, these problems have conflicting objectives.Minimization of CEnet is considered as a major aim in most existing studies on HEMS, while minimizing TBD is another significant aim in the consumer perspective.After the serious concerns over the environmental issues, the role of DRSREOD-based HEMS to reduce the emissions due to the local as well as due to the centralized generation by the utility has recently started to be investigated.Various methods have been used in the recent research to take into account important trade-offs between CEnet, TBD and TEMiss.The most widely used approach is the WSMD [12,13,15,[18][19][20].This is an a priori technique that converts MOO into single-objective optimization (SOO) in order to achieve one solution.These techniques do not give the feasible relations among the objectives to allow the consumer to choose among specific preferences.The methods may miss a number of good solutions for a specific user regardless of their preference standards.Analysis of the trade-offs among the above-mentioned aims is considered very pertinent, which enables consumers in decison-making from the set of the diverse set of optimized outcomes.The posteriori method known as Pareto-based multi-objective optimization is dependent on the Pareto dominance idea, which provides the diverse set of possible solutions for multiple objectives.
The concept of MOO problems for a HEMS using decision vector Tst and m objectives for Pareto-based optimization is computed as: minimizing the objective vector F(Tst) = [F 1 (Tst), F 2 (Tst), .., F m (Tst)], following the mentioned constraints.When Tst 1 is better than Tst 2 in any one of the given objectives and is not considered worse in any other than the solution Tst 1 , which is said to dominate another Tst 2 .The Pareto-optimal set is composed of the set of non-dominated solutions.The recently introduced MOGA includes features to implement Pareto optimization.Pareto-optimal sets providing optimal trade-offs between CEnet, TBD and TEMiss for a DRSREODLDG-based HEMS have been calculated in this work by using MOGA and the Pareto optimization features as described in Algorithm 1.One hundred of these aforementioned trade-offs computed through MOGA/PO based heuristic are processed to enhance eco-efficiency making use of a constrained filtration mechanism as discussed in Sections 4.4 and 4.5.

Constrained Filtration of Trade-Offs to HEMS Optimization
The trade-off solutions achieved for a multi-objective optimization problem can be passed through an adequately designed filtration mechanism in order to apply a constraint on any one or more of the specified trade-offs.Such filtration mechanism enables harnessing the trade-offs with enhanced efficiency.In this research, a filtration mechanism has been proposed to screen out the trade-offs with larger values of TEMiss as related to the trade-offs for CEnet and TBD.This mechanism comprises of an AVCF and an ASCF.AVCF makes use of an average value of the trade-off parameter TEMiss to filter out the trade-offs with above average and extremely high values of TEMiss.While ASCF takes into account an average surface fit for TEMiss in terms of CEnet and TBD to screen out the trade-offs with higher values of TEMiss.ASCF has been developed using polynomial based regression technique elaborated in Section 4.5.The formulation for AVCF and ASCF and their application to harness eco-efficient trade-off solutions are presented in Section 4.

Regression Based Surface Fitting Techniques to Develop ASCF
Regression models are used to establish a relation between the dependent and the independent variables in a set of data (z i , x i , y i ).In order to develop a surface for the variable z in terms of variables x and y, a regression model is fit to the set of input data.Major models for surface fitting include interpolant, lowess, and polynomial.A polynomial surface fit model has been used in this research for its flexibility in application to the input data.For polynomial surfaces, a general model is designated as Poly(kl), where k is the degree in x and l is the degree in y.The degree of x in each term will be less than or equal to k, and the degree of y in each term will be less than or equal to l.The maximum for the sum of k and l is m.The degree of the polynomial is the maximum of the values of k and l.A linear regression model (LRM) for i number of observations for the independent variables (x i , y i ) defines a curved surface for the dependent variable z i in a 3D-space.The LRM for the surface in terms of an order-m polynomial may be represented as: where 0 ≤ k + 1 ≤ m. z i is called a dependent variable or regressand, and x i and y i are called independent variable or regressors.The first term in Equation ( 9) is deterministic and represents the conditional mean of z i based on the given values of x i and y i .The second term u i , called the error term, is random in nature.The term is added or subtracted from the first term to realize the actual data.A(k, l) are called regression coefficients (RCs).In LRMs, they are assumed to be fixed numbers.The term linear in LRM refers to the linearity of the RCs.The fitting of a model is based on the estimates for the RCs.The estimation is carried out based on the minimization of the least squares of the error term called a least square method.Based on Equation ( 9), u i is the difference between the actual value of z i and the one obtained through the regression.For an optimal linear coefficient (LC) for surface fitting, the sum of the squared error term (SSE) to be minimized is given as follows: As u i is a function of RCs, the minimum value of the SSE may be computed by taking partial derivatives of the same with respect to each of the RCs and equating the expression to zero.Based on the estimated a(k, l), a sample model z is is formulated and the error term is also known as residuals, which is computed as e i = z i − z is .The regression coefficients a(k, l) are the estimators of A(k, l) and ei is the estimator of the error term u i .The numerical values taken by an estimator are called estimates.The least-squares solution to the problem is a vector a(k, l), which estimates the unknown vector of coefficients A(k, l).In present research, SSE has been used to estimate the model fit for an average surface for TEMiss in terms of CEnet and TBD.It is assumed that the observed data is of equal quality and thus has constant variance; however, the fit might be unduly influenced by the data of poor quality.Methods like weighted-least-squares regression are applied to reduce the influence of the low quality data on estimating the model fits [31,32].In present research, the use of AVCF to screen out the trade-offs with extremely high values of TEMiss inherently improves the data quality for the model fit for ASCF.
The model fit in this research is based on minimization of SSE that may be improved using methods like minimization of root mean square error and root mean square error of approximation.However, the improved well-fitting is of minimal value, if it is not based on the ideas from a theory validation point of view and in such cases an extensive cross-validation is required [33].Accordingly, various polynomial models fit, 25 in number, have been examined for their capabilities to reduce the average value of trade-off parameters for TEMiss and the number of diverse trade-offs available for CEnet and TBD after the application of filtration mechanism.

Algorithms for Eco-Efficient Trade-Offs for DRSREODLDG-Based HEMS
A three-step approach has been used to achieve eco-efficient trade-off solutions for DRSREODLDG-based HEMS.In step-1, schemes for optimal HEMS operation and the related trade-offs for TEMiss, CEnet, and TBD are computed using Algorithm 1.The trade-off solutions thus obtained are passed through a filtration mechanism to harness the ones with bare minimum TEMiss in terms of CEnet and TBD using Algorithm 2. The filtration mechanism is completed in two stages designated as step-2 and step-3.In step-2, an AVCF for TEMiss is developed and applied to the primary trade-offs to filter out the ones with extremely high and above average values of TEMiss.The remaining trade-offs are then passed to step-3 to screen out the trade-offs with even the marginally higher values of TEMiss.In step-3, an ASCF is used to filter out the trade-offs with TEMiss parameters residing above the average surface fit for TEMiss.Eco-efficient solutions including bare minimum TEMiss and diversified trade-offs for CEnet and TBD are thus achieved for DRSREODLDG-based HEMS operation.The followings algorithms have been proposed to harness the eco-efficient tradeoff solutions for DRSREODLDG-based HEMS: • Algorithm 1 to generate primary tradeoffs for DRSREODLDG-based HEMS (Step-1).• Algorithm 2 for filtration mechanism to harness eco-efficient tradeoffs for DRSREODLDG-based HEMS (Step-2 and step-3).

Algorithm 1 to Generate Operating Schemes and the Primary Tradeoffs for DRSREODLDG-Based HEMS (Step-1)
This algorithm computes a set of primary tradeoff solutions for optimized HEMS operation based on MS of SHAs synergized with the optimal dispatch of the PV system, the SB, the grid, and an LDG.The LDG supplies the load only during LSD hours in coordination with the PV unit and the SB.Tradeoffs for CEnet, TBD, and TEMiss are based on the underlying scheme for HEMS operation.At the start, vector Tst for SHAs is generated that is followed by the production of Pschd vector.The PV system is regarded as the preferred source to directly supply Pschd.The dispatch planning is mainly based on the excess PV energy in each slot denoted by Pres which is the arithmetic difference between Ppv and Pschd.Two main cases arise with regard to the relative values of these two quantities and in each case, SoC, the maximum charge/discharge rates, the grid status and the power from the LDG play major roles in the dispatch.In the first case, where excess PV energy is available, as shown on line 21, the energy is stored in the SB if SoC is less than its maximum value; otherwise, it is sold to the grid.However, during LSD hours, the excess energy that would be sold to the grid is instead supplied to a dummy load as shown on line 24.The SB charging state depends on the condition given on line 30.If a value other than Pres is computed, it indicates that either the maximum charge rate or the limiting value of SoC is restricting the complete storage of the excess PV energy in the SB.Hence, any excess energy left after charging the SB is sold to the grid, as shown on line 35.However, during LSD hours, the excess energy that would have been sold to the grid is instead supplied to a dummy load, as shown on line-33.In the second case, in which Ppv is less than or equal to Pschd, as shown on line 41, the PV energy is insufficient to completely supply the load.The residual energy, in this case, will be supplied from the grid if SoC is less than or equal to its minimum limit or from the SB otherwise.Moreover, the SB will still also not be discharged if cheap energy is available from the grid as given on line 42.However, during LSD hours, the LDG will supply the load in place of the grid, as shown on line 46.SB shall supply the load only during peak hours when the cost of energy is greater than a maximum price limit.The discharging state of the SB depends on the condition given on line 50.If the minimum computed value is equal to the maximum discharge rate or to the residual capacity of the SB before discharging to the minimum SoC, then one of these constraints is restricting the ability to supply the full load from the SB, and the remaining load must be supplied from the grid, as shown on lines 52 and 54.However, during LSD hour, the LDG will supply the remaining load in place of the grid as shown on line 57.For each slot in the scheduling horizon, one of the above two cases will hold, the vectors Ppv, Pgd, Pds, and Pgn will be computed for dispatch accordingly.Similarly, the loads for each slot is computed for Pschd, Pch, Pdl, and Psold.TEMiss is computed (applying EFT) for the net generation from LDG as shown on line 68.CEnet is computed by arithmetically adding CE (applying PE/IBR), cost of generation from LDG (applying PEg) and cost of energy sold to the grid (applying PE f ) as shown on line 69.The TBD values are computed on line 77 after adding TBD(A) and TBD(D) on lines 72 and 74.The values for the mentioned objective functions are computed for each MOGA iteration.The algorithm provides Pareto optimal sets comprising one hundred operating schemes for HEMS in terms of Tst and the related tradeoffs for CEnet, TBD and TEMiss.

Algorithm 2 for Filtration Mechanism to Harness Eco-Efficient Trade-Offs for DRSREODLDG-Based HEMS (Step-2 and Step-3)
The algorithm completes the filtration process in two steps as stated below: Step-2: An AVCF based on the average value of TEMiss is developed taking into account all of the primary trade-offs generated through Algorithm 1 as shown on line-2.The residuals for TEMiss (TEMiss_Resid_avg) for each solution are then computed as given on line-3.A trade-off solution with the value of TEMiss_Resid_avg less than 0 indicates an above average value for TEMiss, and all such trade-offs are filtered out at the step shown on line-4.The trade-off solutions with average (or less than average) TEMiss values are collected and forwarded to step-3 for further processing as shown on line-5.
Step-3: An ASCF based on the average surface fit (using polynomial-based regression) is developed making use of the trade-off solutions forwarded from step-2 as shown on line-9.The residuals for TEMiss (TEMiss_Resid_avgs) for each solution are then computed by taking the difference between the TEMiss and the average surface fit of TEMiss computed in terms of CEnet and TBD as shown on line-10.A trade-off solution with the value of TEMiss_Resid_avgs less than 0 indicates the TEMiss value greater than the respective value on the average surface fit, and all such trade-offs are filtered out at a step shown on line-11.The remaining trade-off solutions with the TEMiss values equal to (or less than) the respective values on the average surface fit are selected and declared final eco-efficient trade-offs for DRSREODLDG-based HEMS operation as shown on line-12.

Simulations for DRSREODLDG-Based HEMS Operation and the Filtration Mechanism to Harness Eco-Efficient Trade-Off Solutions
The simulations are conducted using MATLAB 2015 and are reported in Section 6.1 based on Algorithm 1.These results show the validity of MOGA/PO based heuristic for DRSREODLDG-based HEMS to compute operational schemes for SHAs in terms of vector Tst and the primary trade-offs for CEnet, TBD and TEMiss.The results of simulations enable analyzing the trends exhibited by the trade-off parameters taking into consideration vital factors affecting these parameters.The critical analysis of the primary trade-offs enabled designing a filtration mechanism to extract desired set of eco-efficient trade-off solutions with minimal TEMiss.The simulations reported in Section 6.2 are based on Algorithm 2. They demonstrated the validity of the filtration mechanism to harness eco-efficient trade-offs.Regression based polynomial formulations and the procedure to finalize the model fits for the proposed mechanism are also elaborated in Section 6.2.Simulations have been conducted for the following: -DRSREODLDG-based HEMS operation to compute primary trade-offs for HEMS ( based on Algorithm 1/step-1).
-Application of filtration mechanism to harness eco-efficient trade-offs for HEMS (based on Algorithm 2/step-2 and step-3).

Simulations for DRSREODLDG-Based HEMS Operation to Compute Primary Trade-Offs Using Algorithm 1
Simulations were performed to validate DRSREODLDG-based HEMS operation using Algorithm 1. Operating schemes for SHAs in terms of Tst and the primary trade-offs were computed.The trends exhibited by the trade-off parameters were analyzed.Critical analysis for validating the relation between the trade-off parameter: TEMiss and the trade-offs for CEnet, TBD, enabled designing a filtration mechanism required to harness the desired eco-efficient trade-off solutions with minimal TEMiss from a large set of primary trade-offs.
For the simulations, the 2-stage ToU pricing tariff with an IBR value of 1.4 are considered.This comprises of hourly price of 15 Cents/kWh during the peak hours from 7:00 p.m. to 11:00 p.m. (slot numbers 115-138) hours and a hourly price of 9 Cents/kWh for the whole day, as displayed in Figure 4.The IBR factor threshold is considered as the power demand of 2.4 kW.A feed-in tariff, PE f , valued at 0.7 times of PE was considered for the PV energy sold to the grid.The primary trade-offs for CEnet, TBD and TEMiss, generated through simulation for an optimal DRSREODLDG-based HEMS operation, are presented in Table 10.Due to space limitation, the related Tst vector is not shown in this table (however, it is presented with the final eco-efficient trade-offs in Table 13).The primary trade-offs are graphically shown in Figure 5.The trends exhibited by the trade-off parameters and the relationship between them has been analyzed to approach a filtration mechanism that enables harnessing trade-offs with diversified options for CEnet, TBD and minimal value of TEMiss.Refer to Table 10, each trade-off solution is related to a unique Tst.The decision vector Tst is generated through the MOGA based on the vectors STslot and ENslot.The vector Tst for each of the solution corresponds to a unique demand profile, Pschd.To supply this demand, a dispatch scheme for energy sources and ESS based on the parameters Ppv, Pgd, Pgn, Pch, and Pds is computed through the heuristic proposed in Algorithm 1. Preferably, the load is supplied from the PV unit.The extra energy from the PV unit is stored in the SB after supplying the load.The SB is discharged to supply the load during the peak hours for making use of the stored energy.The LDG supplies the load in coordination with the SB during LSD hours only.The excess PV energy is sold to the grid.This is designated as Psold after supplying the load and charging the SB.However, during the LSD hours, the excess energy from the PV ought to be dissipated into the dummy load viz.designated as Pdl.The PV, SB, and the charger system are considered a part of the existing infrastructure and their cost is not included in computation.

Time Slots
The trade-off parameter CEnet is based on the dispatch from various sources to supply the scheduled load and the energy sold to the utility according to Equation (3).The rates for energies including PE, PE f and PEg in different slots play vital role in the computation of CEnet.The loss of the harnessed PV energy due to the unavailability of the grid, given by Pdl, is another important factor affecting the value of CEnet.The parameter TEMiss primarily depends on the energy supplied by the LDG, Pgn, during LSD hours.The EFT for the LDG is also important while evaluating TEMiss.The TBD is based on the time shift of SHAs from their preferred times of operation and is computed using Equation (7).The relationships between the trade-off parameters for the primary trade-off solutions are graphically presented in Figures 2 and 6  The trends exhibited by the trade-off parameters comprising CEnet, TBD and TEMiss are analyzed in subsections below.The primary trade-offs with extreme values of the parameters have especially been investigated.

Trends for CEnet
The objective to minimize the CEnet is mainly based on the following factors: 1.
Maximized usage of the PV energy to supply the load directly: This avoids the loss of energy in the SB due to storage/re-use of the PV energy while supplying the load (a net loss of 20% has been assumed for the SB).The energy thus saved enables to reduce the demand from the grid and the LDG which ultimately results in a reduced value of CEnet.

2.
Maximized usage of the stored PV energy to supply the load during the peak hours: This reduces the energy to be supplied from the grid during the peak hours as well as from the LDG during the peak LSD hours, which results in a reduced value of CEnet.

3.
Selling of the extra PV energy to the utility: A direct usage of the energy from the PV unit is better than selling it to the utility as PE f is generally lesser than the PE (PE f is assumed as 70% of the PE).However, it is beneficial to sell the PV energy to the utility, if a surplus of it is available after supplying Pschd and the charging load.The above-mentioned factors enable reducing the CEnet parameter through an optimal use of the PV energy based on the PE, PE f , PEg and the SB efficiency.Other factors to reduce CEnet parameter include the followings: 4.
Load shifting towards the off-peak hours: The load left after being supplied from the PV and the SB unit should have been shifted towards off-peak hours.This shifting minimizes the CEnet based on the tariff PE.

5.
Load to be supplied by the LDG during LSD hours: The algorithm enables supply of the energy from the LDG during LSD hours.If more load is shifted towards LSD hours, LDG is required to supply that load in coordination with the PV/SB at a higher cost of energy (PEg), which results in an increased value of CEnet.

6.
Loss of the harnessed PV energy: The dummy load Pdl has been identified as a factor of vital importance for reducing CEnet.To investigate the variations in CEnet parameter based on the above-mentioned six factors, solution-1 and solution-100 with the maximum and the minimum values of CEnet are analyzed as case studies.The analysis is based on the related HEMS operation including the power profiles for the loads and the dispatch scheme for the power sources and the SB.Solution-1 shows a CEnet value of 52.87 Cents, the largest of all solutions.This largest value of CEnet may be analyzed based on the above-mentioned factors by focusing on the power profiles for this solution shown in Figure 8.  First, a very small portion of the load (Pschd) has been supplied directly from the PV energy that is available from time slot no.37. Some of the available PV energy has been used to charge the SB while most of the PV energy is sold to the utility at cheap rates (PE f equals 70% of PE).A part of the load, instead of being supplied directly from the PV unit, is shifted towards the off-peak slots and supplied from the grid at the off-peak time rate.This load thus has been supplied at a net 30% increased cost of the energy as compared to the cost of energy sold to the grid.Second, a load larger than the capacity of the SB is shifted towards the peak-time slots.An average load of 0.21 kWh is thus supplied from the grid during peak time slot nos.132-134.The CEnet could be reduced if the load exceeding the capacity of the SB was shifted towards off-peak time.Third, a net load of 0.348 kWh has been supplied from the LDG during LSD based slot nos.139-144 at a rate of PEg (viz.higher than PE).This load is based on NSHAs only and it can not be shifted.However, the LDG also supplies a load of 0.068 kWh during slot no.102 that may be shifted towards the grid/PV to reduce the CEnet.Fourth, the least of part of the the load has been shifted within the PV harnessed LSD hours starting from slot nos.61 and 97.Under this scenario, 1.87 kWh of the PV energy has been lost/dumped during slot nos.63-66 and slot nos.97-101.More load could be shifted towards the mentioned slots to minimize the loss of the harnessed energy from the PV and thus to reduce the CEnet.In brief, a load shifting resulted in a non-optimal use of the PV energy, a very large value of the Pdl and other aforementioned factors resulted in the largest value of CEnet for this solution.Solution-100, on the other hand, exhibits the lowest CEnet value of 26.22 Cents that is again based on the aforementioned factors.The lowest value of CEnet may again be analyzed by focusing on the corresponding power profiles for the solution as shown in Figure 9.  First, a larger portion of the load (Pschd), as compared to solution-1, has been supplied directly from the PV that is available from time slot no.37.The harnessed PV energy has been used to charge the SB as well as to supply the maximum of the load, while a smaller value of the PV energy is sold to the utility at cheap rates.Second, the remaining load viz.smaller as compared to solution-1 has been shifted towards the peak time slots so that the SB is able to supply most of the said load.Accordingly, an average load of 0.189 kWh is left to be supplied by the grid during the peak time slot nos.135-137, which is smaller as compared to the same load in solution-1.Third, the LDG supplies a total energy of 0.054 kWh during slot nos.100-101, which is smaller as compared to the same parameter in solution-1.Fourth, most of the load has been shifted towards the PV harnessed LSD hours and hence Pdl exhibits a minimal value 0.11 kWh.In brief, a load shifting enabling an optimal use of the PV energy minimized the value of Pdl and other aforementioned factors resulted in the lowest CEnet for this solution.Similarly, the solutions with an intermediate value of CEnet may also be validated by focusing the same above-mentioned factors affecting CEnet.

Trends for TBD
The value of TBD is based on the total time shifts of the SHAs from the preferred times (STslot or ENslot based on type of scheduling) provided by the consumers.It depends on the decision vector Tst and computed using Equation ( 7) through Algorithm 1.The simulations reveal an exponential relation between the CEnet and TBD as shown in Figure 10.The TBD increases exponentially while reducing the CEnet.The relationship between the CEnet and TBD is very important in the context of the consumer's welfare.The optimal solutions provide diverse choices to the consumer for trade-offs between CEnet and TBD.However, it has been observed that CEnet cannot be reduced beyond a specific value after the TBD reaches a knee-point value.A knee-point value of 0.48 for TBD may be realized from Figure 10.On the other hand, the relation between the TBD and TEMiss for DRSREODLDG-based HEMS is highly un-even as shown in Figure 11.Such relations are not possible to be defined using standard techniques.

Trends for TEMiss
The variation in TEMiss is analyzed based on the primary trade-offs presented in the Figure 6/Table 10. Figure 6 exhibits an extremely uneven variations in TEMiss as related to CEnet (and TBD), especially around the center of the data.The solution-23 with the largest, solution-27 with the smallest and solution-73 with moderate values of TEMiss are analyzed as case studies.
Solution-23 exhibits a TEMiss value of 1.55 Lbs., the largest of all solutions.The value of TEMiss parameter depends on the profile for Pgn parameter.The profile for this solution is analyzed by focusing on the power/emission profiles shown in Figure 12.The value of TEMiss mainly depends on the operation of the LDG during four number of LSD hours discussed as follows.The loads shifted in the first LSD hour (starting at slot no.61) and in the third LSD hours (the peak time hour starting at slot no.121) are completely supplied by the PV and the SB, respectively.Thus, in actuality, the LDG has to operate only during the second LSD hour (starting at slot no.97) and during the fourth LSD hour (starting at slot no.139) to supply the shifted load as neither the grid nor the SB is available to supply within these hours.During the fourth LSD hour, a fixed load made up of NSHAs is supplied by the LDG completely.As no other source is available to supply during this hour, the fixed load has been supplied by the LDG in all scenarios.Focusing the second LSD hour, PV is available to supply the shifted load; however, the demand exceeding the energy harnessed from the PV (named excess demand) is only supplied through the LDG.This excess demand to be supplied by the LDG during the second LSD hour combined with the fixed demand in the fourth LSD hour, in fact, determines the net value of TEMiss.A maximum shifting of the excess demand out of the second LSD hour results in the minimization of the TEMiss.For solution-23, a maximum excess demand supplied through the LDG during the second LSD hour resulted in a maximum TEMiss value of 1.55 Lb. for this solution.The CEnet parameter in this scenario assumes a near average value of 43.96 Cents that is based on the combined effect of the related parameters' values including: a PV energy loss of 1.09 kW; a supply of an average load of 0.2 kWh through the grid during peak time slot nos.132-134; and a maximum supply of 0.98 kWh of energy from the LDG at a higher cost of value (PE g ).Solution-27 exhibits a TEMiss value of 0.56 Lbs, the lowest in all solutions and the power profiles shown in Figure 13.The minimum value of TEMiss in this scenario is because of zero loading of LDG during the second LSD hour.On the other hand, the CEnet parameter shows a near average value of 43.57Cents that is very similar to the CEnet value in solution-23.The value is again based on the combined effect of the related parameters' values including: a PV energy loss of 1.75 kW; supply of an average load of 0.23 kW by the grid during the peak time slot nos.132-134; and a minimum supply of 0.35 kWh of energy from the LDG at a higher cost, PEg.Solution-73 shows a moderate TEMiss value of 1.20 Lbs.corresponding to the power profiles shown in Figure 14.The excess load during the second LSD hour has not been completely shifted out of this hour and so the same has been supplied through the LDG.The TEMiss for this solution, therefore, is higher as compared to its value for solution-27.A much lower CEnet of value 31.99Cents as compared to the value of CEnet in solution-27 is based on a more efficient shifting of the load and a smaller value of Pdl in solution-73.The relation between TEMiss parameter and the trade-offs for CEnet and TBD is analyzed based on the primary trade-offs (sorted on CEnet), presented in Table 10.The trade-offs are graphically shown in Figure 15.Based on the variations in TEMiss, the data may be divided into three classes.
Class-1, including solution nos.01-20 at the beginning of the data, class-2, including solution nos.21-73 around the center of the data, and class-3, including solution nos.74-100 at the end of the data.
Class-1 is characterized by the trade-offs with minimal values of TEMiss combined with maximal values of CEnet; and class-3 by the trade-offs with minimal values of both of the TEMiss and CEnet parameters.Whereas class-2 around the middle section of the data, including more than 50% of the trade-offs, exhibits a highly un-even/irregular trend for TEMiss as related to the trade-offs for CEnet and TBD, it includes an un-even distribution of the data with the minimal, average as well as extremely high values of the TEMiss.Such trends indicate the presence of numerous solutions with comparable values of the trade-offs for CEnet and TBD, however with large variations in the related values for TEMiss.Solutions-23 and 27, graphically shown as points A and B respectively in Figure 15, are an example of such large variation in the TEMiss parameter.For comparable values of (43.96, 0.22) and (43.57, 0.2) for CEnet and TDB, the solutions exhibit extremely varied values of 1.55 Lbs.(maximum of all solutions) and 0.56 Lbs.(minimum of all solutions) for TEMiss.Solution-69 and solution-72 shown as points C and D are another example of similar large variations in TEMiss.For comparable values of (32.37, 0.33) and (32.01, 0.35) for CEnet and TDB, the solutions show largely varied respective values of 0.65 Lbs. and 1.23 Lbs. for TEMiss.
Figure 15 reveals a large number of data points especially in class-2 exhibiting large variations in TEMiss with very small corresponding variation in the respective trade-off values of CEnet and TBD.The finding regarding the existence of a large number of multiple comparable trade-offs for CEnet and TBD with extremely varied values of TEMiss in the primary trade-offs was exploited to design a mechanism to harness eco-efficient trade-offs for DRSREODLDG-based HEMS operation.A filtration mechanism was proposed to screen out the trade-offs with larger values of TEMiss in order to harness eco-efficient trade-offs with minimal TEMiss and a set of diverse trade-offs for CEnet and TBD.The proposed mechanism, based on an average value constraint filter and an average surface based constraint filter, is presented in Algorithm 2.

Simulation for Filtration Using AVCF (Step-2)
This step includes the formulation and application of a constraint filter based on the average value of TEMiss for the primary trade-off solutions presented in Table 10.In the following are the software and hardware tools used to demonstrate the solution space, and to formulate and apply the filter to validate the AVCF based filtration:   Where average_sur f ace_ f it is the value of emission obtained through the average surface fit based polynomial for the respective CEnet and TBD trade-off.In addition, TEMiss_Resid_avgs is the decision element for the filter.The exclude option provided with the surface fit function has been used to screen out the trade-offs based on the formulation of the decision element.As per the formulation for TEMiss_Resid_avgs in this research, a trade-off solution with a negative value of the decision element TEMiss_Resid_avg indicates the average surface fit for TEMiss.The application of ASCF thus screened out the trade-offs with higher values of TEMiss lying above the average surface fit for TEMiss.(11) The proposed polynomial model is based on the coefficients (with 95% confidence bounds) as follows: The eco-efficient solutions harnessed after the application of Poly41 surface filter are graphically shown in Figure 17.The final set of trade-off solutions for eco-efficient operation of DRSREODLDG-based HEMS are tabulated as Table 13.

Pdl (kWh)
Ts1 Ts2 Ts3 Ts4 Ts5 Ts6 Ts7 Ts8 Ts9 Ts10 Ts11 Ts12 Ts13 Ts14  CTA given in Table 14 along with the respective scheduled times Tst in Table 13 enables consumer to select the best eco-efficient options according to his needs after consulting a diverse set of current optimized choices for CEnet, TBD and minimal TEMiss.

Figure 2 .
Figure 2. Un-even trends for TEMiss as related to CEnet and TBD.

Figure 5 .
Figure 5.Primary trade-off solutions with un-even surface for TEMiss generated through Algorithm 1. .

Figure 8 .
Figure 8. Power and emission profiles for DRSREODLDG-based HEMS operation for solution-1.

Figure 9 .
Figure 9. Power and emission profiles for DRSREODLDG-based HEMS operation for solution-100.

Figure 10 .
Figure 10.Relation between CEnet and TBD for a DRSREODLDG-based HEMS.

Figure 12 .
Figure 12.Power and emission profiles for DRSREODLDG-based HEMS operation for solution-23.

Figure 13 .
Figure 13.Power and emission profiles for DRSREODLDG-based HEMS operation for solution-27.

Figure 15 .
Figure 15.Variations in total GHG emissions during the scheduling horizon TEMiss with the related trade-offs for net cost of energy CEnet and time-based discomfort TBD.
application: • TEMiss_Resid_avg = average(TEMiss) − TEMiss,• Exclude = TEMiss_Resid_avg < 0,Where TEMiss_Resid_avg is the decision element for the filter.The exclude option provided with the surface fitting function can be used to screen out the trade-offs based on the formulation of the decision element.As per the formulation for TEMiss_Resid_avg, a trade-off solution with a negative value of the decision element TEMiss_Resid_avg indicates the above average value for TEMiss.The application of AVCF thus screens out the trade-offs with extremely high as well as above the average values of TEMiss.The function of the AVCF to screen out the un-desired trade-offs with larger values of TEMiss are graphically shown in Figure16.The selected solutions after the application of the AVCF are presented in Table11 .

Figure 16 .
Figure 16.Application of AVCF to screen out the trade-offs with larger TEMiss values.6.1.6.Simulation for Filtration Using ASCF (Step-3) This step includes the formulation and application of a constraint filter based on the average surface fit for TEMiss.The average surface fit for TEMiss in terms of CEnet and TBD is generated using polynomial based regression for the trade-offs achieved after the application of AVCF presented
LDG Figure 1.Architecture for demand response synergized with renewable energy sources and energy storage systems optimal dispatch integrated with load shedding-compensating dispatchable generator based home energy management system for a smart home [10].

Table 3 .
Related work of proposed method for eco-efficient DRSREODLDG-based HEMS operation.

Table 5 .
Shiftable home appliances and scheduling specifications.

Table 6 .
Non-shiftable home appliances considered for scheduling.

Table 8 .
SB and inverter specifications.

Table 14 .
Critical trade-off analysis for eco-efficient DRSREODLDG-BASED HEMS operation.