# Patent

## Title
HVAC system control integrated with demand response, on-site energy storage system and on-site energy generation system

## Abstract

A thermal behavior model of a building may be constructed based on time series data. Based on the constructed thermal behavior model, forecasted zone temperature and energy usage for a next control time period may be predicted. An objective function may be constructed based on at least a dynamically priced grid energy cost, occupant comfort matrix, and one or more of: energy storage system cost and associated operational cost, and energy generation system cost and associated green house emission cost and associated operational cost. Constraints may be constructed based on at least the forecasted zone temperature values and energy usage values for the next control time period. A control profile of the HVAC system and sourcing decision of energy load of the HVAC system may be determined simultaneously based on the objective function and the plurality of constraints.

## Background
FIELD
The present application relates generally to buildings and energy consumption, and more particularly to control systems and controlling Heating Ventilation and Air Conditioning (HVAC) systems.
BACKGROUND
Buildings such as commercial and residential buildings consume a large amount of energy, and approximately half of the energy consumption is taken up by Heating Ventilation and Air Conditioning (HVAC) systems in the buildings. A traditional approach to optimized HVAC control in a building is to compute an optimal control profile (e.g., a set point profile of each heating and/or cooling zone, supply flow rate or supply temperature of air handling unit (AHU)) that minimizes the total energy consumption while satisfying thermal comfort (e.g., zone temperature and humidity) and physical limitations of HVAC equipment (e.g., supply temperature, supply flow rate of AHU). Such traditional approaches typically assume that the electricity price is constant throughput the day.
Another known approach develops an HVAC control method that minimizes the total energy costs considering demand response signal (e.g., dynamic, time of day pricing of electricity) while satisfying thermal comfort, and an example of such method is a commercial product called BuildingIQ. However, that approach does not decide how the load of HVAC system resulting from the optimized control is optimally sourced through energy supplies, e.g., grid electricity with demand response, on-site stored electricity (e.g., lead acid battery) and on-site generated electricity (e.g., combined heat and power (CHP) generator).
BRIEF SUMMARY
A method of controlling an HVAC system, in one aspect, may comprise receiving a time series data of at least a zone temperature associated with a building zone, energy usage associated with the building zone, building operations associated with the building zone and ambient temperature. At least some of the time series data may be received from a building automation controller that is connected to one or more sensors and/or one or more meters in a building. The method may also comprise constructing a thermal behavior model of the building zone based on the time series data. The method may further comprise predicting based on the constructed thermal behavior model at least forecasted zone temperature values and energy usage values for a next control time period. The method may also comprise constructing an objective function based on at least a dynamically priced grid energy cost, occupant comfort matrix, and one or more of: energy storage system cost and associated operational cost, and energy generation system cost and associated green house emission cost and associated operational cost. The method may further comprise constructing a plurality of constraints based on at least the forecasted zone temperature values and energy usage values for the next control time period. The method may also comprise determining a control profile of the HVAC system and sourcing decision of energy load of the HVAC system simultaneously based on the objective function and the plurality of constraints. The method may further comprise transmitting the control profile to the building automation controller for controlling the HVAC system in accordance with the control profile. The method described herein may be used for controlling HVAC systems in commercial buildings, residential building and other types of buildings such as hospital, schools, warehouses, airports and hotels.
A system for controlling an HVAC system, in one aspect, may comprise at least one processor. A processor may be operable to receive a time series data of at least a zone temperature associated with a building zone, energy usage associated with the building zone, building operations associated with the building zone and ambient temperature. The processor may be operable to receive at least some of the time series data from a building automation controller that is connected to one or more sensors and/or one or more meters in a building. The processor may be further operable to construct a thermal behavior model of the building zone based on the time series data. The processor may be further operable to predict based on the constructed thermal behavior model at least forecasted zone temperature values and energy usage values for a next control time period. The processor may be further operable to construct an objective function based on at least a dynamically priced grid energy cost, occupant comfort matrix, and one or more of (i) energy storage system cost and associated operational cost, and (ii) energy generation system cost and associated green house emission cost and associated operational cost. The processor may be further operable to construct a plurality of constraints based on at least the forecasted zone temperature values and energy usage values for the next control time period. The processor may be further operable to determine a control profile of the HVAC system and sourcing decision of energy load of the HVAC system simultaneously based on the objective function and the plurality of constraints. The processor may be further operable to transmit the control profile to the building automation controller for controlling the HVAC system in accordance with the control profile. The processor may be used for controlling HVAC systems in commercial buildings, residential building and other types of buildings such as hospital, schools, warehouses, airports and hotels.
A computer readable storage medium storing a program of instructions executable by a machine to perform one or more methods described herein also may be provided.
Further features as well as the structure and operation of various embodiments are described in detail below with reference to the accompanying drawings. In the drawings, like reference numbers indicate identical or functionally similar elements.



BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS
FIG. 1 is a diagram illustrating an optimized demand response control of HVAC system with on-site energy storage system and on-site energy generation system in one embodiment of the present disclosure.
FIG. 2 is a diagram illustrating an optimized demand response control of HVAC system in another embodiment of the present disclosure.
FIG. 3 is a diagram illustrating an optimized demand response control of HVAC system with on-site energy storage system in one embodiment of the present disclosure.
FIG. 4 is a diagram illustrating an optimized demand response control of HVAC system with on-site energy generation system in one embodiment of the present disclosure.
FIG. 5 is a diagram illustrating an overview of an HVAC system control integrated with demand response, on-site energy storage system and/or on-site energy generation system in one embodiment of the present disclosure.
FIG. 6 is a diagram illustrating a control flow for controlling HVAC system by integrating it with a demand response, on-site energy storage system and/or on-site energy storage system in one embodiment of the present disclosure.
FIG. 7 is a diagram illustrating a structure of a non-linear autoregressive with external input (NARX) neural network used in one embodiment of the present disclosure to model thermal behavior.
FIG. 8 is a flow diagram illustrating a method of controlling an HVAC system in one embodiment of the present disclosure.
FIG. 9 illustrates a schematic of an example computer or processing system that may implement an HVAC control system in one embodiment of the present disclosure.



DETAILED DESCRIPTION
A method and computer system may be presented that in one embodiment compute control profiles of an HVAC system as well as how the power (load) needed by the HVAC system is optimally sourced through optimized combination of grid purchased energy with demand response, on-site stored energy and on-site generated energy. In one aspect, the control profiles are optimal. The control profile is then communicated to the controller (e.g., a Building Automation System) to control the HVAC system. In one embodiment, the method includes a Model Predictive Control (MPC) framework for a building HVAC control that reduces energy costs of HVAC operations proactively using an optimal control method and considering the dynamic price of grid purchased electricity, on-site stored electricity, and/or on-site generated electricity. The method and/or computer system can serve as an energy reduction and demand response tool that not only optimize the energy costs in buildings but also reduce energy production, stabilize energy supply (grid), and promote smart grid.
In one embodiment, the method and/or computer system of the present disclosure may compute optimal control profile of HVAC system that minimizes the total energy costs and green house gas (GHG) emission considering, e.g., (1) demand response signal, (2) on-site energy storage system (3) on-site energy generation system, while at the same time, satisfying thermal comfort (e.g., zone temperature and humidity), physical limitations of HVAC equipment and physical limitations of one or more energy storage systems (ESS) and one or more energy generators. For example, the method and/or computer system of the present disclosure may optimally compute how much electricity to generate and/or store and when or at what time (e.g., when cost or price is low) and compute how much of the stored or generated electricity or grid purchased electricity to use and when or at what time for operations of an HVAC system. For instance, in addition to using the grid purchased electricity which may be priced dynamically, e.g., based on demand and supply, the method and/or computer system of the present disclosure may also determine how much of the stored and/or generated electricity to use and when to use the stored and/or generated electricity. The method and/or computer system of the present disclosure may power the HVAC system from the optimal combination of grid electricity, on-site stored electricity and on-site generated electricity.
In one aspect, using a predictive model that predicts HVAC energy demand at various times, weather conditions and operational situations in a building (e.g., whether the building or a zone of the building is in operation or not), the method and/or computer system of the present disclosure in one embodiment compute how energy demand should be changed as well as how much energy from the storage should be used, how much energy to store simultaneously and how much energy to generate on-site and use it for HVAC operation.
In another aspect, the method and/or computer system of the present disclosure in one embodiment source the load of optimized HVAC control profile (control variable values over time) through an optimized combination of on-site energy storage (e.g., lead acid battery etc.) or on-site generation (e.g., CHP generator) along with demand response responded grid purchased energy.
The method and/or computer system of the present disclosure in one embodiment determine how the load of an HVAC system resulting from an optimized control is optimally sourced through energy supplies, e.g., grid electricity with demand response, on-site stored electricity (e.g., lead acid battery) and on-site generated electricity (e.g., CHP generator). In another aspect, the method and/or computer system of the present disclosure in one embodiment may allow for generating and storing energy at anytime (e.g., when price or cost is low) and using it at anytime for operations of HVAC system.
A Model Predictive Control (MPC) method in the present disclosure in one embodiment may model or describe the thermal behavior of zones (i.e., spaces) in a building by using two sets of equations. The one equation is the equation of state variables:




  
    
      
        
          
            d
            
              x
              ⁡
              
                (
                t
                )
              
            
          
          
            d
            t
          
        
        =
        
          f
          ⁡
          
            (
            
              
                x
                ⁡
                
                  (
                  t
                  )
                
              
              ,
              
                u
                ⁡
                
                  (
                  t
                  )
                
              
            
            )
          
        
      
    
    
      
        (
        1
        )
      
    
  




The thermal behavior can also be represented by algebraic expression as

xt=f(xt−1,xt−2, . . . ,ut,ut−1,ut−2, . . . )  (2)

Here, xt is the value of state variable x at time t, and ut is the value of control variables u at time t. The equations (1) and (2) describe that the state variable values for a time period of interest depend on the state variable values in the previous time periods and control variable values at the time of interest as well as the previous time periods. The state variable and control variable are vectors, and include multiple components. Examples of the state variable include zone temperature of zone z at time t, Tt,zzone, day of the week (DOW), time of the day (TOD) and ambient temperature, Ttamb. Examples of the control variables include Tt,zsp, set point of zone z and time t, Tt,zsys, supply temperature of air handling unit (AHU) at zone z, and {dot over (m)}t,zsys, supply flow rate of air handling unit (AHU) at zone z. For the example of state variable Tt,zzone and control variable Tt,zsp, the example state variable equation (2) becomes

Tt,zzone=f(Tt−1,zzone,Tt−2,zzone, . . . ,xt,xt−1,xt−2, . . . ,Tt,zsp,Tt−1,zsp,Tt−2,zsp, . . . )  (3)

where x are vectors of other state variables such as day of the week (DOW), time of the day (TOD) and ambient temperature, Ttamb. Equations (1)-(3) are referred to as equations of system state.

Another set of equations that describe the thermal phenomena of zones in a building is the equations of system output:

yt=h(yt−1,yt−2 . . . ,xt,xt−1,xt−2, . . . ,ut,ut−1,ut−2, . . . )  (4)

Here, yt is the value of system output variable y at time t. An example of system output variable is Pt,zHVAC, power consumption of HVAC system for zone z at time t. For the example of state variable Tt,zzone and control variable Tt,zsp, output variable can be Pt,zHVAC, and the example system output equation (3) becomes

Pt,zHVAC=h(Pt−1,zHVAC,Pt−2,zHVAC, . . . ,Tt,zzone,Tt−1,zzone,Tt−2,zzone, . . . ,xt,xt−1,xt−2, . . . ,Tt,zsp,Tt−1,zsp,Tt−2,zsp, . . . )  (5)

An example method for modeling the thermal behavior (equation 2 and 3) of the building zone is artificial neural network modeling. As HVAC system behavior is dynamic and non-linear, a non-linear autoregressive with external input (NARX) data driven model can capture its properties and states. NARX is a feed-forward time delay neural network, which maps input data to an output, with additional external input (see FIG. 7). The NARX network includes three layers: input layer, hidden layer and output layer. However, deep NARX network, with multiple hidden layers can also be considered. The choice of the activation function influences the optimization scheme. The methodology of the present disclosure may consider three types of activation functions: hyperbolic tangent, symmetric saturated liner transfer function and hard-limit transfer function. Choosing the first transfer function (hyperbolic tangent) results in Mixed-Integer Non-Linear Programming (MINLP) problem, while choosing the other two function (symmetric saturated liner transfer function and hard-limit transfer function) results in Mixed-Integer Linear Programming (MILP) problem. The latter is easier to solve. In one embodiment of the present disclosure, the neural network is trained on historic time-series data (e.g., few weeks' time series data). The entire dataset for neural network training may be randomly divided into three contiguous blocks: training, validation and testing.
This network results in the following algebraic equation:

y(t)=TF2[W3·TF1{W1·(x(t),x(t−1), . . . ,x(t−d1))+W2·(y(t−1), . . . ,y(t−d2))+b1}+b2]  (6)

where W1, W2, W3 are weight matrices, b1, b2 are biases, d1, d2 are network time delays and TF1, TF2 are transfer functions, of which TF2 may be chosen to be linear. x(t) is the network input (e.g., an array of input) at time t and y(t) is the network output (e.g., zone temperature or power).

In one embodiment of the present disclosure, an optimal control modeling may be described by the following equations:
Compute the Optimal Control Profile, u*(τ)




  
    
      
        
          
            min
            ⁢
            
              
            
            ⁢
            
              g
              ⁡
              
                (
                
                  x
                  ,
                  u
                
                )
              
            
          
          =
          
            
              ∫
              
                t
                k
              
              
                t
                
                  k
                  +
                  T
                
              
            
            ⁢
            
              
                φ
                ⁡
                
                  (
                  
                    
                      x
                      ⁡
                      
                        (
                        τ
                        )
                      
                    
                    ,
                    
                      u
                      ⁡
                      
                        (
                        τ
                        )
                      
                    
                  
                  )
                
              
              ⁢
              
                
              
              ⁢
              
                d
                τ
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            φ
            ⁡
            
              (
              
                
                  x
                  ⁡
                  
                    (
                    τ
                    )
                  
                
                ,
                
                  u
                  ⁡
                  
                    (
                    τ
                    )
                  
                
              
              )
            
          
          :=
          
            
              
                α
                1
              
              ⁡
              
                (
                EnergyCost
                )
              
            
            +
            
              
                α
                2
              
              ⁡
              
                (
                GHG_Cost
                )
              
            
            +
            
              
                α
                3
              
              ⁡
              
                (
                Comfort
                )
              
            
            +
            ...
          
        
        ⁢
        
          

        
        ⁢
        
          
            
              s
              .
              t
              .
              
                
                  d
                  
                    x
                    ⁡
                    
                      (
                      τ
                      )
                    
                  
                
                
                  d
                  t
                
              
            
            =
            
              f
              ⁡
              
                (
                
                  x
                  ,
                  u
                
                )
              
            
          
          ,
          
            
              or
              ⁢
              
                
              
              ⁢
              
                f
                ⁡
                
                  (
                  
                    x
                    ,
                    u
                  
                  )
                
              
            
            =
            0
          
        
        ⁢
        
          

        
        ⁢
        
          
            
              h
              ⁡
              
                (
                
                  x
                  ,
                  u
                
                )
              
            
            =
            0
          
          ,
          
            

          
          ⁢
          
            
              x
              L
            
            ≤
            x
            ≤
            
              x
              H
            
          
          ,
          
            

          
          ⁢
          
            
              u
              L
            
            ≤
            u
            ≤
            
              u
              H
            
          
          ,
        
      
    
    
      
        (
        7
        )
      
    
  




other constraints (physical limitations).
In the above equation (7), L represents lower bound, H represents upper bound, g represents an overall objective function, φ represents an objective function value at certain time, τ represents t or time, f represents an equation of system state, h represents an equation of system output (e.g., power), k represents an arbitrary index value to indicate time.
A method for modeling the optimal control problem in one embodiment may be formulated as an objective function and a set of constraints. The following describes the notations, the objective function and the set of constraints in one embodiment:
Indices and sets:
tεT: Time periods
zεZ: Zones
Decision Variables:
pt Power (electricity) (for HVAC) purchased from grid at time t [kW]
gt Power of CHP generated by generator at time t [kW]
stin Power charged by ESS at time t [kW]
stout Power discharged from ESS at time t [kW]
stac Total amount of energy stored in the ESS at time t [kWh]
Pt,zHVAC Total power used by HVAC system at time t for zone z [kW]
ut Binary variable, on/off status of generator at time t
vt Binary variable, on/off status of ESS at time t
Tt,zZone Zone temperature at time t [° C.]
Tt,zsp Zone set point at time t [° C.] (Control Variable)
Parameters:
Cte Cost of grid purchased electricity at time t [$/kW for one hour]:
Cg Cost of natural gas purchased [$/GJ or $/kWh]:
Ge GHG emission rate of grid purchased electricity (source) at time t [kg CO2e/kWh]
Gg GHG emission rate of natural gas purchased (source) [kg CO2e/kWh]
λs Efficiency of energy charged to ESS
λd Efficiency of energy discharged from ESS
λsd Rate of energy self-discharged from ESS [%]
λg Efficiency of CHP co-gen
α1 Weight of cost in the objective function: e.g., 1
α2 Weight of emission in the objective function: e.g., 5
α3 Weight of comfort in the objective function: e.g., 2
α4 Weight of EGS operational cost in the objective function: e.g., 1
α5 Weight of ESS operational cost in the objective function: e.g., 1
gmax, gmin Maximum and minimum generation amount [kW]
sin,max, sin,min Maximum and minimum storage (charge) rate in unit time [kW]
sout,max, sout,max Maximum and minimum discharge rate in unit time [kW]
sac,max, sac,min Maximum and minimum total storage amount [kWh]
δ Time unit (to convert power (kW) to energy (kWh), [hr]
rtup Max ramp-up rate of generator [kW/hr]
rtdown Max ramp-down rate of generator [kW/hr]
Tt,zzone* Target zone temperature of time t in zone z [° C.]
Stg Startup cost (capital) of generator
Sts Startup cost (capital) of storage system
Objective Function:




  
    
      
        
          min
          
            
              p
              t
            
            ,
            
              s
              t
              in
            
            ,
            
              s
              t
              out
            
            ,
            
              g
              t
            
            ,
            
              T
              
                t
                ,
                z
              
              sp
            
          
        
        ⁢
        
          
            ∑
            
              t
              ∈
              T
            
          
          ⁢
          
            
          
          ⁢
          
            [
            
               
              
                
                  
                    α
                    1
                  
                  ⁢
                  
                    {
                    
                      
                        
                          C
                          t
                          e
                        
                        ⁡
                        
                          (
                          
                            
                              p
                              t
                            
                            +
                            
                              
                                s
                                in
                              
                              
                                λ
                                s
                              
                            
                          
                          )
                        
                      
                      +
                      
                        
                          C
                          g
                        
                        ⁢
                        
                          
                            g
                            t
                          
                          
                            λ
                            g
                          
                        
                      
                    
                    }
                  
                
                +
                
                  
                    α
                    2
                  
                  ⁢
                  
                    {
                    
                      
                        
                          G
                          e
                        
                        ⁡
                        
                          (
                          
                            
                              p
                              t
                            
                            +
                            
                              
                                s
                                in
                              
                              
                                λ
                                s
                              
                            
                          
                          )
                        
                      
                      +
                      
                         
                        
                          
                             
                            
                              
                                G
                                g
                              
                              ⁢
                              
                                
                                  g
                                  t
                                
                                
                                  λ
                                  g
                                
                              
                            
                            }
                          
                          +
                          
                            
                              α
                              3
                            
                            ⁢
                            
                              
                                ∑
                                
                                  z
                                  ∈
                                  Z
                                
                              
                              ⁢
                              
                                
                                  
                                    
                                    
                                      
                                        T
                                        
                                          t
                                          ,
                                          z
                                        
                                        
                                          zone
                                          *
                                        
                                      
                                      -
                                      
                                        f
                                        ⁡
                                        
                                          (
                                          ·
                                          )
                                        
                                      
                                    
                                    
                                  
                                  2
                                
                                2
                              
                            
                          
                          +
                          
                            
                              α
                              4
                            
                            ⁢
                            
                              S
                              t
                              g
                            
                            ⁢
                            
                              μ
                              t
                            
                          
                          +
                          
                            
                              α
                              5
                            
                            ⁢
                            
                              S
                              t
                              s
                            
                            ⁢
                            
                              v
                              t
                            
                          
                        
                        ]
                      
                    
                  
                
              
            
          
        
      
    
    
      
        (
        8
        )
      
    
  




Constraints:




  
    
      
        
          
            
              
                
                  T
                  
                    t
                    ,
                    z
                  
                  zone
                
                =
                
                  f
                  (
                  
                    
                      T
                      
                        
                          t
                          -
                          1
                        
                        ,
                        z
                      
                      zone
                    
                    ,
                    
                      T
                      
                        
                          t
                          -
                          2
                        
                        ,
                        z
                      
                      zone
                    
                    ,
                    ...
                    ⁢
                    
                      
                    
                    ,
                    
                      
                        x
                        _
                      
                      
                        t
                        -
                        1
                      
                    
                    ,
                    
                      
                        x
                        _
                      
                      
                        t
                        -
                        2
                      
                    
                    ,
                    ...
                    ⁢
                    
                      
                    
                    ,
                    
                      T
                      
                        t
                        ,
                        z
                      
                      sp
                    
                    ,
                    
                      T
                      
                        
                          t
                          -
                          1
                        
                        ,
                        z
                      
                      sp
                    
                    ,
                    
                      T
                      
                        
                          t
                          -
                          2
                        
                        ,
                        z
                      
                      sp
                    
                    ,
                    ...
                  
                  )
                
              
              ,
            
          
          
            
              
                ∀
                
                  t
                  ∈
                  T
                
              
              ,
              
                z
                ∈
                Z
              
            
          
        
      
    
    
      
        (
        9
        )
      
    
  
  
    
      
        
          
            
              
                
                  P
                  
                    t
                    ,
                    z
                  
                  HVAC
                
                =
                
                  h
                  (
                  
                    
                      P
                      
                        
                          t
                          -
                          1
                        
                        ,
                        z
                      
                      HVAC
                    
                    ,
                    
                      P
                      
                        
                          t
                          -
                          2
                        
                        ,
                        z
                      
                      HVAC
                    
                    ,
                    ...
                    ⁢
                    
                      
                    
                    ,
                    
                      T
                      
                        t
                        ,
                        z
                      
                      zone
                    
                    ,
                    
                      T
                      
                        
                          t
                          -
                          1
                        
                        ,
                        z
                      
                      zone
                    
                    ,
                    
                      T
                      
                        
                          t
                          -
                          2
                        
                        ,
                        z
                      
                      zone
                    
                    ,
                    ...
                    ⁢
                    
                      
                    
                    ,
                    
                      
                        x
                        _
                      
                      
                        t
                        -
                        1
                      
                    
                    ,
                    
                      
                        x
                        _
                      
                      
                        t
                        -
                        2
                      
                    
                    ,
                    ...
                    ⁢
                    
                      
                    
                    ,
                    
                      T
                      
                        t
                        ,
                        z
                      
                      sp
                    
                    ,
                    
                      T
                      
                        
                          t
                          -
                          1
                        
                        ,
                        z
                      
                      sp
                    
                    ,
                    
                      T
                      
                        
                          t
                          -
                          2
                        
                        ,
                        z
                      
                      sp
                    
                    ,
                    ...
                  
                  )
                
              
              ,
            
          
          
            
              
                ∀
                
                  t
                  ∈
                  T
                
              
              ,
              
                z
                ∈
                Z
              
            
          
        
      
    
    
      
        (
        10
        )
      
    
  
  
    
      
        
          
            
              
                
                  ∑
                  
                    z
                    ∈
                    Z
                  
                
                ⁢
                
                  P
                  
                    t
                    ,
                    z
                  
                  HVAC
                
              
              =
              
                p
                t
              
            
          
          
            
              (
              
                example
                ⁢
                
                  
                
                ⁢
                shown
                ⁢
                
                  
                
                ⁢
                in
                ⁢
                
                  
                
                ⁢
                
                  FIG
                  .
                  
                    
                  
                  ⁢
                  2
                
              
              )
            
          
          
            
              ∀
              
                t
                ∈
                
                  {
                  
                    1
                    ,
                    ...
                    ⁢
                    
                      
                    
                    ,
                    T
                  
                  }
                
              
            
          
        
      
    
    
      
        (
        11
        )
      
    
  
  
    
      
        
          
            
              
                
                  ∑
                  
                    z
                    ∈
                    Z
                  
                
                ⁢
                
                  P
                  
                    t
                    ,
                    z
                  
                  HVAC
                
              
              =
              
                
                  p
                  t
                
                +
                
                  
                    s
                    t
                    out
                  
                  ⁢
                  
                    λ
                    d
                  
                
              
            
          
          
            
              (
              
                
              
              ⁢
              
                example
                ⁢
                
                  
                
                ⁢
                shown
                ⁢
                
                  
                
                ⁢
                in
                ⁢
                
                  
                
                ⁢
                
                  FIG
                  .
                  
                    
                  
                  ⁢
                  3
                
              
              )
            
          
          
            
              ∀
              
                t
                ∈
                
                  {
                  
                    1
                    ,
                    ...
                    ⁢
                    
                      
                    
                    ,
                    T
                  
                  }
                
              
            
          
        
      
    
    
      
        (
        12
        )
      
    
  
  
    
      
        
          
            
              
                
                  ∑
                  
                    z
                    ∈
                    Z
                  
                
                ⁢
                
                  P
                  
                    t
                    ,
                    z
                  
                  HVAC
                
              
              =
              
                
                  p
                  t
                
                +
                
                  g
                  t
                
              
            
          
          
            
              (
              
                example
                ⁢
                
                  
                
                ⁢
                shown
                ⁢
                
                  
                
                ⁢
                in
                ⁢
                
                  
                
                ⁢
                
                  FIG
                  .
                  
                    
                  
                  ⁢
                  4
                
              
              )
            
          
          
            
              ∀
              
                t
                ∈
                
                  {
                  
                    1
                    ,
                    ...
                    ⁢
                    
                      
                    
                    ,
                    T
                  
                  }
                
              
            
          
        
      
    
    
      
        (
        13
        )
      
    
  
  
    
      
        
          
            
              
                
                  ∑
                  
                    z
                    ∈
                    Z
                  
                
                ⁢
                
                  P
                  
                    t
                    ,
                    z
                  
                  HVAC
                
              
              =
              
                
                  p
                  t
                
                +
                
                  
                    s
                    t
                    out
                  
                  ⁢
                  
                    λ
                    d
                  
                
                +
                
                  g
                  t
                
              
            
          
          
            
              (
              
                example
                ⁢
                
                  
                
                ⁢
                shown
                ⁢
                
                  
                
                ⁢
                in
                ⁢
                
                  
                
                ⁢
                
                  FIG
                  .
                  
                    
                  
                  ⁢
                  1
                
              
              )
            
          
          
            
              ∀
              
                t
                ∈
                
                  {
                  
                    1
                    ,
                    ...
                    ⁢
                    
                      
                    
                    ,
                    T
                  
                  }
                
              
            
          
        
      
    
    
      
        (
        14
        )
      
    
  
  
    
      
        
          
            
              
                T
                
                  t
                  ,
                  z
                
                
                  zone
                  ,
                  L
                
              
              ≤
              
                T
                
                  t
                  ,
                  z
                
                zone
              
              ≤
              
                T
                
                  t
                  ,
                  z
                
                
                  zone
                  ,
                  H
                
              
            
          
          
            
              ∀
              
                t
                ∈
                
                  {
                  
                    
                      t
                      L
                    
                    ≤
                    t
                    ≤
                    
                      t
                      H
                    
                  
                  }
                
              
            
          
        
      
    
    
      
        (
        15
        )
      
    
  




where Tt,zzone,L is a lower bound for Tt,zzone and Tt,zzone,H is an upper bound for Tt,zzone, and tL is a lower boundary of building operating time (e.g., hour) and tH is an upper boundary of building operating time (e.g., hour)
Co-gen (for scenario 3 and 4)

gmin·ut≦gt≦gmax·ut  (16)

−δrtdown≦gt+1−gt≦δrtup  (17)

utε{0,1},gt≧0,  (18)

ESS (for Scenario 2 and 4)

stac=[st−1ac+δ(stin−stout)](1−λtsd)  (19)

δstout≦stac  (20)

sin,min·vt≦stin≦sin,max·vt  (21)

sout,min·vt≦stout≦sout,max·vt  (22)

sac,min≦stac≦sac,max  (23)

vtε{0,1},stin≧0,stout≧0  (24)

In this optimal control method, the method and/or computer system of the present disclosure in one embodiment determine a profile of a control variable, e.g., Tt,zsp, set point of zone z and time t for a future period (e.g., next 24 hours) that minimizes total energy costs of operating HVAC system subject to thermal behavior (Eq. 9 and 10) inside a zone in a building, energy balance (Eq. 11, 12, 13 and 14) in a zone, and physical limitations of equipments, e.g., on-site energy storage equipment (Eq. 19, 20, 21, 22, 23 and 24) and on-site energy generation equipment (Eq. 16, 17 and 18). The control profile can also be of Tt,zsys, supply temperature of air handling unit (AHU) at zone z, and {dot over (m)}t,zsys, supply flow rate of air handling unit (AHU) at zone z. The total cost (i.e., objective function) includes but is not limited to the costs related to energy usage, greenhouse gas emission deviation from comfort temperature range for building occupants and operational costs of on-site operating energy storage system and on-site energy generation system.
The method and/or computer system of the present disclosure in one embodiment may determine a profile of a variable (e.g., control variable) for one or more zones. A profile may be determined for each of multiple zones, if multiple zones are considered, simultaneously using an objective function and associated constraints. For instance, multiple profiles may be determined simultaneously.
Constraints of equations (9) and (10) describe thermal behavior inside a zone in a building. Constraints of equations (11)-(14) describe energy balance, i.e., the total used energy is equal to the sum of the sourced energy.
Constraints of equations (16)-(18) are associated with energy generation. The constraint of equation (16) describes that energy generation amount is bounded by a minimum and a maximum value. The minimum and maximum values may be given, e.g., by a manufacturer of an energy generator. The constraint of equation (17) describes that a change is the amount of generated energy from a time period to the next time period is bounded by ramp-up and ramp-down values. The ramp-up and ramp-down values may be given, e.g., by a manufacturer of an energy generator. The constraint of equation (18) describes that the status of a generator is on or off and that the quantity of energy generated is a positive number.
Constraints of equations (19)-(24) are associated with energy storage system. The constraint of equation (19) describes energy balance of the energy store: e.g., the total energy accumulated at time t is accounted by the energy accumulated at previous time, energy charged or stored at time t, and energy discharged at time t considering self-discharge efficiency factor. The constraint of equation (20) describes that the amount of energy discharged is not greater than the amount that is accumulated. The constraint of equation (21) describes that the amount of energy charged is bounded by minimum and maximum rates. The constraint of equation (22) describes that energy discharge capacity is bounded by minimum and maximum values. The constraint of equation (23) describes that charged or accumulated energy is bounded by minimum and maximum values. The constraint of equation (23) describes that the status of an energy storage system is on or off and that the quantity of energy stored and discharged is a positive number. The minimum and maximum values may be given by a manufacturer of an energy storage system.
In one embodiment, the optimal control problem above may be formulated as a MINLP problem, e.g., because of the non-linear expression of the thermal behavior modeled (Eq. 9 and 10), which includes sigmoid activation functions, in artificial neural network. Solving the MINLP problem is usually difficult, and one way to simplify the problem is to approximate the MINLP to MILP by linearizing the sigmoid function. In this method, Eq. 9 and 10 are reformulated to a MILP, which may be solved using a variety of methods, depending on its size. One way to solve the MILP problem is with the IBM-ILOG CPLEX solver.
FIG. 1 is a diagram illustrating an optimized demand response control of HVAC system with on-site energy storage system and on-site energy generation system in one embodiment of the present disclosure. The total energy used by an HVAC system at time t (102) may be expressed as a function of control variables (e.g., set point temperature) and energy supplied to the HVAC system at time t (108). The energy supplied to the HVAC system 108 is expressed in terms of decision variables. Examples of decision variables may include, power purchased from grid at time t, power charge rate, power discharge rate and power generated at time t. For example, power may include purchased grid electricity 104, which may be directly used for an HVAC system 106. Purchased grid electricity may be also stored or accumulated as shown at 112. The amount of stored energy may depend on charge efficiency 110 and self-discharge rate of an energy storage system. The stored electricity may be discharged subject to discharge efficiency 114 of the energy storage system, to be supplied as part of the energy supplied to HVAC 108. Energy supplied to HVAC 108 also may be sourced from generated energy. For example, another source of power, natural gas 116 may be converted to power 120 subject to the efficiency factor of a generator 118 and supplied to the HVAC system as part of the energy supplied to HVAC 108.
A methodology of the present disclosure in one embodiment shown in FIG. 1 may determine simultaneously both the amounts of energy sources to be supplied and the control variables (e.g., temperatures to set to) over a future period of time that optimize the cost of supplying the energy or power and the desired comfort level for that future period of time. For example, using the objective function of equation (8), the methodology of the present disclosure in one embodiment finds optimal values for control variable profile u(t) (e.g., set point, Tt,zsp) and decision variables p(t), stin (charge rate), stout (discharge rate) and g(t), (energy generation) that minimize the objective function of equation (8) subject to a set of constraints, e.g., described above in equations (9)-(10), (14), (16)-(24).
In the objective function of equation (8), the first term represents a cost of buying power (e.g., gas and electricity), the second term represents a gas emission cost, the third terms represents occupant comfort matrix, e.g., expressed as a deviation from the desired temperature level, and the fourth and fifth terms represents an operational cost for running energy storage system and energy generation system.
FIG. 2 is a diagram illustrating an optimized demand response control of HVAC system in another embodiment of the present disclosure. The total energy used by an HVAC system at time t (202) may be expressed as a function of control variables (e.g., set point temperature) and energy supplied to the HVAC system at time t (208). The energy supplied to the HVAC system 208 is expressed in terms of decision variables. In this embodiment, examples of decision variables may include, power purchased from grid at time t. For example, power may include purchased grid electricity 204, which may be used for an HVAC system 206. Cost of grid purchased electricity may depend on demand response signal at time t, for example, as shown at 210.
A methodology of the present disclosure in one embodiment shown in FIG. 2 may determine simultaneously both the amount of energy source to be supplied and the control variables (e.g., set point temperatures) over a future period of time that optimize the cost of supplying and the desired comfort level for that future period of time. For example, using the objective function:




  
    
      
        
          min
          
            
              p
              t
            
            ,
            
              T
              
                t
                ,
                z
              
              sp
            
          
        
        ⁢
        
          
            ∑
            
              t
              ∈
              T
            
          
          ⁢
          
             
            
              [
              
                
                  
                    α
                    1
                  
                  ⁢
                  
                    C
                    t
                    e
                  
                  ⁢
                  
                    p
                    t
                  
                
                +
                
                  
                    α
                    2
                  
                  ⁢
                  
                    
                      ∑
                      
                        z
                        ∈
                        Z
                      
                    
                    ⁢
                    
                      
                        
                          
                          
                            
                              T
                              
                                t
                                ,
                                z
                              
                              
                                zone
                                *
                              
                            
                            -
                            
                              f
                              ⁡
                              
                                (
                                ·
                                )
                              
                            
                          
                          
                        
                        2
                      
                      2
                    
                  
                
              
              ]
            
          
        
      
    
    
      
        (
        25
        )
      
    
  




the methodology of the present disclosure in one embodiment finds optimal values for control variable profile u(t) (e.g., set point, Tt,zsp) and decision variable p(t) (e.g., purchased grid electricity) that minimize the objective function subject to a set of constraints, e.g., described above in equations (9)-(10) and (11). In the object function of equation (25), the first term represents a cost of buying power (e.g., electricity) and the second term represents occupant comfort matrix, e.g., expressed as a deviation from the desired temperature level.

FIG. 3 is a diagram illustrating an optimized demand response control of HVAC system with on-site energy storage system in one embodiment of the present disclosure. The total energy used by an HVAC system at time t (302) may be expressed as a function of control variables (e.g., set point temperature) and energy supplied to the HVAC system at time t (308). The energy supplied to the HVAC system 108 is expressed in terms of decision variables. In this embodiment, examples of decision variables may include, power purchased from grid at time t, power storage system's charge rate and discharge rate. For example, power may include purchased grid electricity 304, which may be directly used for an HVAC system 306. Purchased grid electricity may be also stored as shown at 312. The amount of stored energy may depend on charge efficiency 110 of the power or energy storage system. The stored electricity may be discharged subject to discharge efficiency 114 of the power or energy storage system, to be supplied as part of the energy supplied to HVAC 308.
A methodology of the present disclosure in one embodiment shown in FIG. 3 may determine simultaneously both the amounts of energy sources to be supplied and the control variables (e.g., temperatures to set to) over a future period of time that optimize the cost of supplying and desired comfort level for that future period of time. For example, using the objective function:




  
    
      
        
          min
          
            
              p
              t
            
            ,
            
              s
              t
              in
            
            ,
            
              s
              t
              out
            
            ,
            
              T
              
                t
                ,
                z
              
              sp
            
          
        
        ⁢
        
          
            ∑
            
              t
              ∈
              T
            
          
          ⁢
          
            
              [
              
                
                  
                    α
                    1
                  
                  ⁢
                  
                    
                      C
                      t
                      e
                    
                    ⁡
                    
                      (
                      
                        
                          p
                          t
                        
                        +
                        
                          
                            s
                            in
                          
                          
                            λ
                            s
                          
                        
                      
                      )
                    
                  
                
                +
              
               
            
            ⁢
            
               
              
                 
                
                   
                  
                    
                      
                        α
                        3
                      
                      ⁢
                      
                        
                          ∑
                          
                            z
                            ∈
                            Z
                          
                        
                        ⁢
                        
                          
                            
                              
                              
                                
                                  T
                                  
                                    t
                                    ,
                                    z
                                  
                                  
                                    zone
                                    *
                                  
                                
                                -
                                
                                  f
                                  ⁡
                                  
                                    (
                                    ·
                                    )
                                  
                                
                              
                              
                            
                            2
                          
                          2
                        
                      
                    
                    +
                    
                      
                        α
                        5
                      
                      ⁢
                      
                        S
                        t
                        s
                      
                      ⁢
                      
                        v
                        t
                      
                    
                  
                  ]
                
              
            
          
        
      
    
    
      
        (
        26
        )
      
    
  




the methodology of the present disclosure in one embodiment finds optimal values for control variable profile u(t) (e.g., set point, Tt,zsp) and decision variables p(t), stin (charge rate) and stout (discharge rate) that minimize the objective function subject to a set of constraints, e.g., described above in equations (9)-(10), (12), (19)-(24). In the object function of equation (26), the first term represents a cost of buying power (e.g., electricity), the second term represents occupant comfort matrix, e.g., expressed as a deviation from the desired temperature level, and the third term represents an operational cost of running an energy storage system.

FIG. 4 is a diagram illustrating an optimized demand response control of HVAC system with on-site energy generation system in one embodiment of the present disclosure. In this embodiment, the total energy used by an HVAC system at time t (402) may be expressed as a function of control variables (e.g., set point temperature) and energy supplied to the HVAC system at time t (408). The energy supplied to the HVAC system 408 is expressed in terms of decision variables. Examples of decision variables in this embodiment may include, power purchased from grid at time t and power generated at time t. For example, power may include purchased grid electricity 404, which may be directly used for an HVAC system 406. Energy supplied to HVAC 408 also may be sourced from generated energy. For example, another source of power, natural gas 410 may be converted to power 414 subject to an efficiency factor of a generator 412 and supplied to the HVAC system as part of the energy supplied to HVAC 408.
A methodology of the present disclosure in one embodiment shown in FIG. 4 may determine simultaneously both the amounts of energy sources to be supplied and the control variables (e.g., set point temperatures) over a future period of time that optimize the cost of supplying energy and the desired comfort level for that future period of time. For example, using the objective function:




  
    
      
        
          min
          
            
              p
              t
            
            ,
            
              g
              t
            
            ,
            
              T
              
                t
                ,
                z
              
              sp
            
          
        
        ⁢
        
          
            ∑
            
              t
              ∈
              T
            
          
          ⁢
          
            
          
          ⁢
          
             
            
              
                
                  [
                   
                
                ⁢
                
                  
                    α
                    1
                  
                  (
                  
                    
                      
                        C
                        t
                        e
                      
                      ⁢
                      
                        p
                        t
                      
                    
                    +
                    
                      
                        C
                        g
                      
                      ⁢
                      
                        
                          g
                          t
                        
                        
                          λ
                          g
                        
                      
                    
                  
                  )
                
              
              +
              
                 
                
                   
                  
                    
                      
                        α
                        2
                      
                      (
                      
                        
                          
                            G
                            e
                          
                          ⁢
                          
                            p
                            t
                          
                        
                        +
                        
                          
                            G
                            g
                          
                          ⁢
                          
                            
                              g
                              t
                            
                            
                              λ
                              g
                            
                          
                        
                      
                      )
                    
                    +
                    
                      
                        α
                        3
                      
                      ⁢
                      
                        
                          ∑
                          
                            z
                            ∈
                            Z
                          
                        
                        ⁢
                        
                          
                            
                              
                              
                                
                                  T
                                  
                                    t
                                    ,
                                    z
                                  
                                  
                                    zone
                                    *
                                  
                                
                                -
                                
                                  f
                                  ⁡
                                  
                                    (
                                    ·
                                    )
                                  
                                
                              
                              
                            
                            2
                          
                          2
                        
                      
                    
                    +
                    
                      
                        α
                        4
                      
                      ⁢
                      
                        S
                        t
                        g
                      
                      ⁢
                      
                        μ
                        t
                      
                    
                  
                  ]
                
              
            
          
        
      
    
    
      
        (
        27
        )
      
    
  




the methodology of the present disclosure in one embodiment finds optimal values for control variable profile u(t) (e.g., set point, Tt,zsp) and decision variables p(t) and g(t) (energy generation amount), that minimize the objective function of equation (27) subject to a set of constraints, e.g., described above in equations (9)-(10), (13), (16)-(18).

In the object function of equation (27), the first term represents a cost of buying power (e.g., gas and electricity), the second term represents a gas emission cost, the third terms represents occupant comfort matrix, e.g., expressed as a deviation from the desired temperature level, and the fourth term represents an operational cost of running an energy generation system.
The objective functions described above in equations (8), (25)-(27) used set point temperature as an example of a control variable. Other examples of control variables that may be determined using the objective functions may include supply temperature and supply flow rate of AHU.
FIG. 5 is a diagram illustrating an overview of an HVAC system control integrated with demand response, on-site energy storage system and on-site energy generation system in one embodiment of the present disclosure. At 502, the thermal behavior models for zone temperature (e.g., described above with reference to equations (1)-(3) and formulated as a constraint in equation (9)) and power usage (e.g., described above with reference to equations (4)-(5) and formulated as a constraint in equation (10)) may be periodically calibrated in a continuous manner. Dynamic and historical data for system variables, control variables, power (HVAC) usage, weather data, and building operations data may be collected or received to calibrate the thermal model. Dynamic and historical data here refers to a time series data, e.g., data received over time, e.g., 15 minute interval data for past 4 weeks.
At 504, an objective function may be constructed, e.g., as a mixed-integer non-linear programming (MINLP) problem. The objective function in one embodiment of the present disclosure may take into consideration a grid electricity for HVAC cost, energy storage system (ESS) cost, energy generation system (EGS) cost, greenhouse gas (GHG) emission cost of ESS and ESG, and occupant comfort matrix (factors or variables associated with providing comfort to occupant, e.g., comfortable level of temperature). For example, the objective function may be formulated as equation (8), equation (25), equation (26) or equation (27).
At 506, the calibrated thermal models of 502 are utilized to generate predicted or forecasted values for zone temperature and power usage over time (e.g., over a control time period) given the values for zone temperature, weather data, and building operations data. The given values for zone temperature may include zone temperature data from previous time periods. The given values for weather data may include ambient temperature from previous time periods and/or the current time period. The given values for building operations data may include calendar data such a day of the week, time of the day, holiday that provide information as to whether the building is operating or not. As an example, at 506, for every interval of time (e.g., every 10 minutes) for the control time period (e.g., the next 24 hours), a zone temperature and power usage may be predicted.
At 508, MINLP model constraints are constructed based on the predicted values. In addition, the modeled constraints may include sourcing the predicted HVAC load (e.g., computed according to equation (10) at 506) as a sum of grid electricity, ESS discharging rate and EGS generation rate, e.g., as expressed in equations (11)-(14), for example, based on which objective function is being constructed. The modeled constraints may also include physical limitations of ESS, e.g., charging capacity, discharging capacity, ESS storage capacity, etc., e.g., as expressed in equations (19)-(24), for example, based on which objective function is being constructed. The modeled constraints may further include physical limitations of EGS, e.g., ramp-up, ramp-down, min/max generation power, etc., as expressed in equations (16)-(18), for example, based on which objective function is being constructed.
At 510, MINLP is solved using the objective function constructed at 504 and constraints constructed at 508 to determine a control variable profile (e.g., temperature set point at every interval of time over the control period for every zone being considered, e.g., Tt,zzone), energy to use (e.g., grid electricity to use, pt) at the corresponding every interval of time over the control period for every zone being considered. Based on the objective function, stored energy to use (e.g., ESS discharge rate, e.g., stin) at the corresponding every interval of time over the control period for every zone being considered, energy to store (e.g., ESS charge rate, e.g., stout), and/or energy to generate and use (e.g., generated electricity, e.g., gt) at the corresponding every interval of time over the control period for every zone being considered, may be also determined.
At 512, the control variable profile may be transmitted or fed into an HVAC controller, e.g., of a building automation system (BAS). Signals that include the data of the control variable profile, e.g., may be transmitted to the HVAC controller. The control variable profile, for example, may include a temperature value (temperature set point) for every interval of time that the HVAC controller may set in order to achieve the desired target zone temperature for that control period of time at the same time minimizing the various costs described in the objective function. For instance, responsive to receiving the control variable profile signals, the HVAC controller may set the temperatures at the corresponding time intervals.
In the scenario in which amount of energy to store (charge) and used (discharge) is computed, a signal to store and discharge energy may be transmitted to an energy storage system with the appropriate data for instructing or controlling the energy storage system accordingly. In the scenario in which energy is generated and used, a signal to generate energy and use the generated energy may be transmitted to an energy generation system with appropriate data for instructing or controlling the energy generation system accordingly.
In one aspect, the system shown in FIG. 5 may function as a feedback control system loop, in which the temperature set by the HVAC controller 512 affects the control and decision variables, which influence the thermal behavior of the zones in the building, and the results of the thermal behavior of the building are in turn used to recalibrate the thermal models at 502, which are used to predict another set of zone temperatures and power usage at 506, used as constraints in 508.
FIG. 6 is a diagram illustrating a control flow for controlling HVAC system by integrating it with a demand response, on-site energy storage system and/or on-site energy storage system in one embodiment of the present disclosure. BACnet is a communications protocol for building automation and control networks. It is an American Society of Heating, Refrigerating and Air-Conditioning Engineers (ASHRAE), American National Standards Institute (ANSI), and International Organization for Standardization (ISO) standard protocol. An MPC server 602 incorporates functionalities of an MPC model described above. For example, MPC server 602 may perform the methodologies described above with respect to reference items 502, 504, 506, 508 and 510 in FIG. 5, and may determine among others one or more control variables, e.g., temperature set point, supply temperature of AHU, supply flow rate of AHU, based on input data such as weather data and building operational schedule (e.g., described with respect to FIG. 5 at 506) and demand response signal (e.g., described with respect to FIG. 5 at 504). The MPC server 602 may be implemented on a processor including memory. The MPC server 602 sends the determined values, e.g., one or more of temperature set point, supply temperature of AHU, or supply flow rate of AHU to a controller (e.g., building automation system controller (BAS)) 604 for controlling a building's HVAC system (610, 612). The MPC server 602 may also receive data such as historic and/or current zone temperature data, past and/or current power usage data, past and/or current temperature set point data from the BAS controller 604. This data may be used by the MPC server 602 in recalibrating the thermal models, e.g., zone temperature and power usage models, e.g., shown at 502 in FIG. 5. The controller 604 may receive information such as the zone temperature data and power usage data from one or more sensors and/or one or more meters 606. For instance, the controller 604 may be connected to at least a sensor or a meter or both. The controller 604 may be connected to multiple of those sensors, meters or both.
Responsive to receiving one or more determined control variables from the MPC server 602, the BAS controller 604 may send signals to devices to control the devices according to the values of the control variables. For instance, the BAS controller 604 may send control signals indicating temperature set points (e.g., Tt,zsp) over time according to the determined temperature set point control profile, to a zone thermostat 608, to set the zone thermostat 608 at the specified temperature at specified time interval. If the control variable determined at the MPC server 602 includes supply temperature (e.g., Tt,zsys) of AHU or supply flow rate (e.g., {dot over (m)}t,zsys) of AHU, the BAS controller 604 may send control signals to an HVAC system, e.g., an AHU 610 or electric heat pump (EHP) with variable air volume (VAV) or a like device to blow air at the determined temperature or flow rate. An AHU 610 or a like device may be connected to a zone thermostat 608 for receiving temperature data. An AHU 610 or a like device may be also connected to a chiller or boiler 612 of the HVAC system.
FIG. 7 is a diagram illustrating a structure of a non-linear autoregressive with external input (NARX) neural network used in one embodiment of the present disclosure to model thermal behavior. An embodiment of the present disclosure models thermal behavior of a building zone by a nonlinear autoregressive neural network (NARNET). FIG. 7 shows a feed-forward time delay neural network, which maps input data to an output, with additional external input, and illustrates the algebraic equation (6). The methodology of the present disclosure in one embodiment may consider three types of activation functions (for example, shown at TF1 in the figure): hyperbolic tangent, symmetric saturated liner transfer function and hard-limit transfer function. TF2 is a linear function in one embodiment.
In one aspect, optimal control problem of the present disclosure may be formulated as a mixed-integer non-linear programming (MINLP) model. The optimal control problem may be solved to minimize the total energy costs of powering HVAC system and corresponding GHG emission and the cost of operating on-site energy storage system (ESS) and on-site energy generation (EGS) system, and satisfy thermal comfort of building occupants within the physical limitation of HVAC equipment and physical limitation of on-site energy storage system and on-site energy generators. The optimal control problem may consider dynamic demand response signal (dynamic time of the day pricing), on-site energy storage system and on-site energy generation system.
As described above, a methodology of the present disclosure in one embodiment may optimally decide control profiles of an HVAC system and how the energy or power (load) needed by the HVAC system is sourced. A methodology of the present disclosure in one embodiment may model an energy load (e.g., power usage) with a machine learning technique such as a neural network technique, and calibrate the model dynamically and continuously based on dynamic profiles of sensor data (and/or meter data). In one aspect, examples control variables in the present disclosure in an HVAC system may include thermostat set point, supply temperature of air handling unit (AHU), and supply flow of AHU. A methodology of the present disclosure in one embodiment may model a gas-fired co-generator in determining a control profile of an HVAC system. A methodology of the present disclosure in one embodiment may use an on-site generation of carbon-emission data (e.g., carbon generated by on-site gas-fired generators which produces electricity that is used to power the HVAC system) in determining a control profile of an HVAC system. In one aspect, a methodology of the present disclosure in one embodiment generates carbon-emission data associated with gas-fired on-site generated energy, and therefore need not depend on a manufacturer or utility company to provide such information.
A methodology of the present disclosure in one embodiment provides for a demand response based Model Predictive Control (MPC) framework that computes control profiles of an HVAC system that minimize costs, considering grid purchased energy at dynamic pricing, on-site stored energy and/or on-site generated energy, GHG emission and thermal comfort of occupants. FIG. 8 is a flow diagram illustrating a method of controlling an HVAC system in one embodiment of the present disclosure. At 802, a time series data is received. The time series data includes at least a zone temperature data associated with a building zone, energy usage data associated with the building zone, building operations data associated with the building zone and ambient temperature data. In one embodiment, at least some of the time series data are received from a building automation controller that is connected to one or more sensors and/or one or more meters in a building. For instance, the building automation controller may be connected to at least a sensor or a meter or both a sensor and a meter.
At 804, a thermal behavior model of the building zone is constructed based on the time series data. The thermal behavior model may include a zone temperature model and energy usage model. In one embodiment, the thermal behavior model is modeled by a nonlinear autoregressive neural network. At 806, forecasted zone temperature values and energy usage values for a next control time period is predicted based on the constructed thermal behavior model.
At 808, an objective function is constructed. In one embodiment, an objective function may be constructed based on a dynamically priced grid energy cost and occupant comfort matrix, and energy storage system cost and associated energy storage system operational cost. In another embodiment, an objective function may be constructed based on a dynamically priced grid energy cost and occupant comfort matrix, and energy generation system cost and associated green house emission cost and associated energy generation system operational cost. Yet in another embodiment, an objective function may be constructed based on a dynamically priced grid energy cost and occupant comfort matrix, and energy storage system cost and associated energy storage system operational cost, and energy generation system cost and associated green house emission cost and associated energy generation system operational cost.
At 810, a plurality of constraints is constructed based on at least the forecasted zone temperature values and energy usage values for the next control time period.
At 812, a control profile of HVAC system and sourcing decision of energy load of the HVAC system are determined simultaneously based on the objective function and the plurality of constraints. In another embodiment, a control profile of HVAC system, sourcing decision of energy load of the HVAC system and an energy storage decision may be determined simultaneously based on the objective function and the plurality of constraints.
At 814, the control profile may be transmitted to the building automation controller for controlling the HVAC system in accordance with the control profile data. In one embodiment, the control profile may include a series of temperature set points over the next control time period for the building automation controller to set one or more zone thermostats. In another embodiment, the control profile may include a series of supply temperature of AHU over the next control time period for the building automation controller to set an HVAC device. Yet in another embodiment, the control profile may include a series of supply flow rate of AHU over the next control time period for the building automation controller to set an HVAC device.
In one embodiment, the sourcing decision of energy load may include sourcing the energy load of the HVAC system from purchased energy, stored energy and generated energy. In another embodiment, the sourcing decision of energy load may include sourcing the energy load of the HVAC system from purchased energy and stored energy. Yet in another embodiment, the sourcing decision of energy load may include sourcing the energy load of the HVAC system from purchased energy and generated energy. The sourcing decision of energy load may include amount of energy to use from which source and at what time interval over the next control time period. The energy storage decision in one embodiment may include amounts of energy to store at which intervals of time over the next control time period. A sourcing decision associated with energy storage and discharge of stored energy may be transmitted to an energy storage system. A sourcing decision associated with energy generation and using generated energy may be transmitted to an energy generation system.
In one aspect, the control profile of HVAC system and the sourcing decision of energy load of the HVAC system may be determined for each of multiple building zones simultaneously based on the objective function and the plurality of constraints.
The method of the present disclosure in one embodiment may compute a cost effective control profile of an HVAC system and a sourcing decision of the energy load of the HVAC system that minimize the total energy costs considering demand response signals and one or more of on-site energy storage system and on-site energy generation system while satisfying thermal comfort of building occupants, a physical limitation of HVAC equipment, a physical limitation of energy storage system and a physical limitation energy generation system by modeling the dynamics of zone temperature and power consumption, adaptively calibrating the models using the dynamics of sensor data (and/or meter data), and solving the model-predictive optimal control problem using an objective function and a set of constraints. An objective function of the cost effective control model may include grid electricity (for operating HVAC) cost, one or more of on-site energy storage system cost and on-site energy generation system cost, greenhouse gas (GHG) emission cost for one or more of grid electricity and on-site energy generation fuel, and occupant comfort matrix. The cost effective control problem may be solved to compute the amount of electricity from grid to be used, the amount of energy discharged from energy storage system and/or the amount of energy charged to the on-site energy storage system and/or the amount of on-site generated electricity.
FIG. 9 illustrates a schematic of an example computer or processing system that may implement a HVAC control system in one embodiment of the present disclosure. The computer system is only one example of a suitable processing system and is not intended to suggest any limitation as to the scope of use or functionality of embodiments of the methodology described herein. The processing system shown may be operational with numerous other general purpose or special purpose computing system environments or configurations. Examples of well-known computing systems, environments, and/or configurations that may be suitable for use with the processing system shown in FIG. 9 may include, but are not limited to, personal computer systems, server computer systems, thin clients, thick clients, handheld or laptop devices, multiprocessor systems, microprocessor-based systems, set top boxes, programmable consumer electronics, network PCs, minicomputer systems, mainframe computer systems, and distributed cloud computing environments that include any of the above systems or devices, and the like.
The computer system may be described in the general context of computer system executable instructions, such as program modules, being executed by a computer system. Generally, program modules may include routines, programs, objects, components, logic, data structures, and so on that perform particular tasks or implement particular abstract data types. The computer system may be practiced in distributed cloud computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed cloud computing environment, program modules may be located in both local and remote computer system storage media including memory storage devices.
The components of computer system may include, but are not limited to, one or more processors or processing units 12, a system memory 16, and a bus 14 that couples various system components including system memory 16 to processor 12. The processor 12 may include a module 10 that performs the methods described herein. The module 10 may be programmed into the integrated circuits of the processor 12, or loaded from memory 16, storage device 18, or network 24 or combinations thereof.
Bus 14 may represent one or more of any of several types of bus structures, including a memory bus or memory controller, a peripheral bus, an accelerated graphics port, and a processor or local bus using any of a variety of bus architectures. By way of example, and not limitation, such architectures include Industry Standard Architecture (ISA) bus, Micro Channel Architecture (MCA) bus, Enhanced ISA (EISA) bus, Video Electronics Standards Association (VESA) local bus, and Peripheral Component Interconnects (PCI) bus.
Computer system may include a variety of computer system readable media. Such media may be any available media that is accessible by computer system, and it may include both volatile and non-volatile media, removable and non-removable media.
System memory 16 can include computer system readable media in the form of volatile memory, such as random access memory (RAM) and/or cache memory or others. Computer system may further include other removable/non-removable, volatile/non-volatile computer system storage media. By way of example only, storage system 18 can be provided for reading from and writing to a non-removable, non-volatile magnetic media (e.g., a “hard drive”). Although not shown, a magnetic disk drive for reading from and writing to a removable, non-volatile magnetic disk (e.g., a “floppy disk”), and an optical disk drive for reading from or writing to a removable, non-volatile optical disk such as a CD-ROM, DVD-ROM or other optical media can be provided. In such instances, each can be connected to bus 14 by one or more data media interfaces.
Computer system may also communicate with one or more external devices 26 such as a keyboard, a pointing device, a display 28, etc.; one or more devices that enable a user to interact with computer system; and/or any devices (e.g., network card, modem, etc.) that enable computer system to communicate with one or more other computing devices. Such communication can occur via Input/Output (I/O) interfaces 20.
Still yet, computer system can communicate with one or more networks 24 such as a local area network (LAN), a general wide area network (WAN), and/or a public network (e.g., the Internet) via network adapter 22. As depicted, network adapter 22 communicates with the other components of computer system via bus 14. It should be understood that although not shown, other hardware and/or software components could be used in conjunction with computer system. Examples include, but are not limited to: microcode, device drivers, redundant processing units, external disk drive arrays, RAID systems, tape drives, and data archival storage systems, etc.
The present invention may be a system, a method, and/or a computer program product. The computer program product may include a computer readable storage medium (or media) having computer readable program instructions thereon for causing a processor to carry out aspects of the present invention.
The computer readable storage medium can be a tangible device that can retain and store instructions for use by an instruction execution device. The computer readable storage medium may be, for example, but is not limited to, an electronic storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any suitable combination of the foregoing. A non-exhaustive list of more specific examples of the computer readable storage medium includes the following: a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), a static random access memory (SRAM), a portable compact disc read-only memory (CD-ROM), a digital versatile disk (DVD), a memory stick, a floppy disk, a mechanically encoded device such as punch-cards or raised structures in a groove having instructions recorded thereon, and any suitable combination of the foregoing. A computer readable storage medium, as used herein, is not to be construed as being transitory signals per se, such as radio waves or other freely propagating electromagnetic waves, electromagnetic waves propagating through a waveguide or other transmission media (e.g., light pulses passing through a fiber-optic cable), or electrical signals transmitted through a wire.
Computer readable program instructions described herein can be downloaded to respective computing/processing devices from a computer readable storage medium or to an external computer or external storage device via a network, for example, the Internet, a local area network, a wide area network and/or a wireless network. The network may comprise copper transmission cables, optical transmission fibers, wireless transmission, routers, firewalls, switches, gateway computers and/or edge servers. A network adapter card or network interface in each computing/processing device receives computer readable program instructions from the network and forwards the computer readable program instructions for storage in a computer readable storage medium within the respective computing/processing device.
Computer readable program instructions for carrying out operations of the present invention may be assembler instructions, instruction-set-architecture (ISA) instructions, machine instructions, machine dependent instructions, microcode, firmware instructions, state-setting data, or either source code or object code written in any combination of one or more programming languages, including an object oriented programming language such as Smalltalk, C++ or the like, and conventional procedural programming languages, such as the “C” programming language or similar programming languages. The computer readable program instructions may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the latter scenario, the remote computer may be connected to the user's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or the connection may be made to an external computer (for example, through the Internet using an Internet Service Provider). In some embodiments, electronic circuitry including, for example, programmable logic circuitry, field-programmable gate arrays (FPGA), or programmable logic arrays (PLA) may execute the computer readable program instructions by utilizing state information of the computer readable program instructions to personalize the electronic circuitry, in order to perform aspects of the present invention.
Aspects of the present invention are described herein with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and/or block diagrams, and combinations of blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer readable program instructions.
These computer readable program instructions may be provided to a processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks. These computer readable program instructions may also be stored in a computer readable storage medium that can direct a computer, a programmable data processing apparatus, and/or other devices to function in a particular manner, such that the computer readable storage medium having instructions stored therein comprises an article of manufacture including instructions which implement aspects of the function/act specified in the flowchart and/or block diagram block or blocks.
The computer readable program instructions may also be loaded onto a computer, other programmable data processing apparatus, or other device to cause a series of operational steps to be performed on the computer, other programmable apparatus or other device to produce a computer implemented process, such that the instructions which execute on the computer, other programmable apparatus, or other device implement the functions/acts specified in the flowchart and/or block diagram block or blocks.
The flowchart and block diagrams in the Figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer program products according to various embodiments of the present invention. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of instructions, which comprises one or more executable instructions for implementing the specified logical function(s). In some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems that perform the specified functions or acts or carry out combinations of special purpose hardware and computer instructions.
The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used herein, the singular forms “a”, “an” and “the” are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms “comprises” and/or “comprising,” when used in this specification, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof.
The corresponding structures, materials, acts, and equivalents of all means or step plus function elements, if any, in the claims below are intended to include any structure, material, or act for performing the function in combination with other claimed elements as specifically claimed. The description of the present invention has been presented for purposes of illustration and description, but is not intended to be exhaustive or limited to the invention in the form disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the invention. The embodiment was chosen and described in order to best explain the principles of the invention and the practical application, and to enable others of ordinary skill in the art to understand the invention for various embodiments with various modifications as are suited to the particular use contemplated.

## Claims

1. A method of controlling an HVAC system, comprising:
receiving a time series data of at least a zone temperature associated with a building zone, energy usage associated with the building zone, building operations associated with the building zone and ambient temperature, at least some of the time series data received from a building automation controller that is connected to at least a sensor or a meter or both in a building;
constructing a thermal behavior model of the building zone based on the time series data;
predicting based on the constructed thermal behavior model at least forecasted zone temperature values and energy usage values for a next control time period;
constructing an objective function based on at least a dynamically priced grid energy cost and associated greenhouse gas emission, occupant comfort matrix, and one or more of:
energy storage system cost and associated operational cost, and
energy generation system cost and associated green house emission cost and associated operational cost;

constructing a plurality of constraints based on at least the forecasted zone temperature values and energy usage values for the next control time period;
determining a control profile of the HVAC system and sourcing decision of energy load of the HVAC system simultaneously based on the objective function and the plurality of constraints; and
transmitting the control profile to the building automation controller for controlling the HVAC system in accordance with the control profile.

2. The method of claim 1, wherein the thermal behavior model comprises a zone temperature model and energy usage.

3. The method of claim 1, wherein the thermal behavior model is modeled by a nonlinear autoregressive neural network.

4. The method of claim 1, wherein the control profile comprises a series of temperature set points of a zone thermostat over the next control time period for the building automation controller to set one or more zone thermostats.

5. The method of claim 1, wherein the control profile comprises a series of supply temperature of an AHU over the next control time period for the building automation controller to set an HVAC device.

6. The method of claim 1, wherein the control profile comprises a series of supply flow rate of an AHU over the next control time period for the building automation controller to set an HVAC device.

7. The method of claim 1, wherein the sourcing decision of energy load comprises sourcing the energy load of the HVAC system from purchased energy, stored energy and generated energy, sourcing the energy load of the HVAC system from purchased energy and stored energy, or sourcing the energy load of the HVAC system from purchased energy and generated energy, wherein the sourcing decision of energy load comprises amount of energy to use from which source and at what time interval over the next control time period.

8. The method of claim 1, wherein the determining a control profile of HVAC system and sourcing decision of energy load of the HVAC system comprises determining the control profile of HVAC system, the sourcing decision of energy load of the HVAC system and an energy storage decision simultaneously based on the objective function and the plurality of constraints.

9. The method of claim 8, wherein the energy storage decision comprises amounts of energy to store at which intervals of time over the next control time period.

10. The method of claim 1, wherein the control profile of HVAC system and the sourcing decision of energy load of the HVAC system are determined for each of multiple building zones simultaneously based on the objective function and the plurality of constraints.

11. A system for controlling an HVAC system, comprising:
a processor;
the processor operable to receive a time series data of at least a zone temperature associated with a building zone, energy usage associated with the building zone, building operations associated with the building zone and ambient temperature,
the processor operable to receive at least some of the time series data from a building automation controller that is connected to at least a sensor or a meter or both in a building,
the processor further operable to construct a thermal behavior model of the building zone based on the time series data,
the processor further operable to predict based on the constructed thermal behavior model at least forecasted zone temperature values and energy usage values for a next control time period,
the processor further operable to construct an objective function based on at least a dynamically priced grid energy cost and associated greenhouse gas emission, occupant comfort matrix, and one or more of (i) energy storage system cost and associated operational cost, and (ii) energy generation system cost and associated green house emission cost and associated operational cost,
the processor further operable to construct a plurality of constraints based on at least the forecasted zone temperature values and energy usage values for the next control time period,
the processor further operable to determine a control profile of the HVAC system and sourcing decision of energy load of the HVAC system simultaneously based on the objective function and the plurality of constraints, and
the processor further operable to transmit the control profile to the building automation controller for controlling the HVAC system in accordance with the control profile.

12. The system of claim 11, wherein the thermal behavior model comprises a zone temperature model and energy usage.

13. The system of claim 11, wherein the thermal behavior model is modeled by a nonlinear autoregressive neural network.

14. The system of claim 11, wherein the control profile comprises:
a series of temperature set points over the next control time period for the building automation controller to set one or more zone thermostats;
a series of supply temperature of an air handling unit over the next control time period for the building automation controller to set an HVAC device; or
a series of supply flow rate of the air handling unit over the next control time period for the building automation controller to set an HVAC device; or
combinations thereof.

15. The system of claim 11, wherein the sourcing decision of energy load comprises sourcing the energy load of the HVAC system from purchased energy, stored energy and generated energy, sourcing the energy load of the HVAC system from purchased energy and stored energy, or sourcing the energy load of the HVAC system from purchased energy and generated energy, wherein the sourcing decision of energy load comprises amount of energy to use from which source and at what time interval over the next control time period.

16. The system of claim 11, wherein the processor further determines an energy storage decision simultaneously with the control profile of HVAC system and the sourcing decision of energy load of the HVAC system based on the objective function and the plurality of constraints.

17. The system of claim 16, wherein the energy storage decision comprises amounts of energy to store at which intervals of time over the next control time period.

18. The system of claim 11, wherein the control profile of HVAC system and the sourcing decision of energy load of the HVAC system are determined for each of multiple building zones simultaneously based on the objective function and the plurality of constraints.

19. A computer readable storage medium storing a program of instructions executable by a machine to perform a method of controlling an HVAC system, the method comprising:
receiving a time series data of at least a zone temperature associated with a building zone, energy usage associated with the building zone, building operations associated with the building zone and ambient temperature, at least some of the time series data received from a building automation controller that is connected to at least a sensor or a meter or both in a building;
constructing a thermal behavior model of the building zone based on the time series data;
predicting based on the constructed thermal behavior model at least forecasted zone temperature values and energy usage values for a next control time period;
constructing an objective function based on at least a dynamically priced grid energy cost and associated greenhouse gas emission, occupant comfort matrix, and one or more of:
energy storage system cost and associated operational cost, and
energy generation system cost and associated green house emission cost and associated operational cost;

constructing a plurality of constraints based on at least the forecasted zone temperature values and energy usage values for the next control time period;
determining a control profile of the HVAC system and sourcing decision of energy load of the HVAC system simultaneously based on the objective function and the plurality of constraints; and
transmitting the control profile to the building automation controller for controlling the HVAC system in accordance with the control profile.

20. The computer readable storage medium of claim 19, wherein the determining a control profile of HVAC system and sourcing decision of energy load of the HVAC system comprises determining the control profile of HVAC system, the sourcing decision of energy load of the HVAC system and an energy storage decision simultaneously based on the objective function and the plurality of constraints.