Patent ID: 12255456

DETAILED DESCRIPTION

Aspects of the present disclosure include demand response methodologies for PFR, including FFR, in response to under or over frequency events. Aspects of the present disclosure also include methods for controlling a fleet of DERs equipped for FFR and quantifying in real time an amount of primary frequency control capacity available in the fleet. Such methods can be used to provide new sources of PFR control and enable participation in ancillary markets, including FFR markets.

In some examples, a fleet of DERs are designed and configured to participate in primary frequency control by configuring the DERs to switch from a conventional packetized energy management (PEM) operating mode to a frequency-responsive operating mode. The methodologies disclosed herein are capable of supporting operation on large-scale transmission systems as well as smaller microgrid systems and other types of low-inertia grids. In some examples, a methodology based on a local control policy, for example, at an individual device level, that is responsive to measured grid frequency and the device's local states, such as the device's need for energy (such as temperature or state of charge (SOC) and/or a packet timer state that indicates an extent of completion of an energy packet. A device that is drawing power may be configured to immediately stop drawing power and in some examples, begin discharging power back to the grid in response to a deviation in the measured grid frequency from a setpoint value and/or a ROCOF exceeding threshold value and according to the local control policy. In some examples, a local control policy interrupts a charging or discharging packet based on a combination of a current energy state of the device (e.g., water temperature, state of charge (SOC), etc.) and the time since it last transitioned into an energy consuming or generating state and/or the time remaining in a current active energy packet. When the local control policy is aggregated across a fleet of DERs, it dynamically prioritizes which devices will have their packets interrupted first during any frequency event on the grid. A DER can be any electrical source or load connected to the grid. Non-limiting examples of DERs include electric water heaters, smart thermostats, air conditioners, pool pumps, batteries, electric vehicle chargers, and refrigerators.

An aggregator in communication with the fleet of DERs may be configured to provide a real time estimate of an amount of primary frequency control capacity available in the fleet. For example, an aggregator may be configured to provide real-time power-frequency response capability curves for a fleet of participating DERs, including participating DERs with known transition time (e.g., those operating under packetized energy management) and/or DERs that have communicated start times when they have transitioned to an ON state and/or time since OFF. In some examples, aggregators are configured to create a power deviation vs. frequency response map that can be used by grid operators to understand how much frequency response reserves are available from the PEM fleet in real-time. As will be appreciated by persons having ordinary skill in the art, such capability can be used to improve the reliability and resilience of a grid and be used to respond to a frequency event more quickly than conventional generators. Such an aggregated DER fleet with FFR capability can also participate in ancillary services markets, including existing, upcoming, and future FFR markets.

In some examples, a packet interruption threshold mechanism is utilized that turns PEM DERs such as batteries and thermostatically controlled loads (TCLs) configured for PEM ON (charge or discharge) or OFF (standby) as a function of a measured frequency and local device information, thereby resulting in a frequency-responsive device. In some examples, devices are equipped with a packet timer that is initiated when the device is turned ON and stopped when a packet is completed. The DER devices may be configured to interrupt a packet before the packet has been completed when a threshold is exceeded, where the threshold depends on locally measured frequency, timer status, and in some examples, the device's energy state (e.g., temperature, or state of charge). In some examples, the local control policy designed for all DERs can be broadcasted and updated as needed, for example, every day/week/month/year depending on, for example, grid planning, scenarios, and predicted grid, weather, or market conditions. The control policy and the resulting threshold determine how devices are automatically prioritized during grid frequency events such as when a large generator is suddenly taken offline due to an unexpected fault.

The methods and systems disclosed herein provide the ability to deliver and accurately estimate, in real time, the aggregate power deviation of an entire fleet of DERs even when the DERs are interacting in a bottom-up decentralized fashion, meaning the individual loads react individually and independently to their local measurements as opposed to a top-down central management control scheme. As opposed to top-down broadcast-based methods of load control, the estimate can be determined online and updated in real time through an aggregator that is aware of the total number of accepted packets and/or participating devices that are ON and have communicated the time when they transitioned to ON at any time, which can be used to accurately estimate the distribution of timers and, additionally, the energy state of the fleet. Armed with this information, at each instant of time, and before any frequency disturbance, a system operator can estimate the available frequency response reserves for primary frequency control that the registered controllable loads can deliver nearly instantly, e.g., in one to two seconds, to improve grid reliability and resilience.

Aspects of the present disclosure are also disclosed in the following publication, which is incorporated by reference herein in its entirety: H. Mavalizadeh, L. A. Duffaut Espinosa and M. R. Almassalkhi, “Decentralized Frequency Control using Packet-based Energy Coordination*,” 2020 IEEE International Conference on Communications, Control, and Computing Technologies for Smart Grids (SmartGridComm), Tempe, AZ, USA, 2020, pp. 1-7, doi: 10.1109/SmartGridComm47815.2020.9302972.

In a packet-based system, loads repeatedly consume energy packets of a predetermined duration and/or amount and cease consuming energy when the packet is complete even if the device still has a need for energy. For example, a water heater will temporarily stop drawing energy when a packet is complete even if the water temperature is below an upper temperature setpoint and then independently determine, according to a local control methodology, whether to make a request to an aggregator for another packet during the next communication interval. As another example, a battery may temporarily stop charging or discharging when a packet is complete even if the battery is not fully charged or discharged and then independently determine whether to make a request to an aggregator to charge or discharge another packet during the next communication epoch. In some examples, each demand-side energy consuming device is equipped with a packet timer that turns the device OFF when the timer reaches a predetermined energy packet time duration indicating the energy packet is complete. Packet-based systems can alternatively have a top-down centralized control architecture or a bottom-up decentralized control architecture. Examples of decentralized bottom-up PEM control methodologies and systems are disclosed in U.S. Pat. No. 10,256,631, titled Systems And Methods For Random-Access Power Management Using Packetization, US Patent Publication No. 2018/0082385, titled Methods And Systems For Secure Scheduling And Dispatching Synthetic Regulation Reserve From Distributed Energy Resources, and US Patent Publication No. 2019/0324415, titled Packetized Energy Management Control Systems and Methods of Using the Same, each of which are incorporated by reference herein in its entirety.

FIG.1illustrates an example of a system100for the distribution and consumption of electricity that is configured with demand side packet-based DERs102that are configured for fast frequency response. System100includes a plurality of DERs102that are physically connected to an electrical distribution network104of a utility106. Each DER102is communicatively coupled to a DER controller108that is configured to provide PEM control functionality to the DER. The DER controllers108are communicatively coupled with a cloud-based DER management system (DERMS)110that includes one or more aggregators112. DERMS110is configured to receive balancing reference signals SREF from, for example a grid operator114of utility106and coordinate the DERs102to track the balancing reference signal SREF. DERMS110may be configured to communicate with utility106to provide real time and forecasted availability of the DERMS and DERs102, for example with information on the availability of DERs102, and information relating to the provision of ancillary services to the utility including primary frequency response services. Utility106also receives information from distribution network104on the state of the distribution network such as such as voltage, frequency, and power flows.

In the illustrated example, for ease of illustration, two aggregators112are shown, with each aggregator in communication with a corresponding group116of two DERs102. In practice, system100can include any number of aggregators112, any number of groups116and any number of DERs102in each group. This is particularly true for bottom-up local PEM control methodologies, which are particularly well suited for scalability to millions of DERs without a corresponding increase in communication and computation overhead cost as compared to top-down centralized control schemes. However, the FFR methodologies disclosed herein may also be applied to systems designed for centralized top-down PEM control methodologies where there is a more centralized control over the energy consumption of the DERs102. DERs102may be grouped in any way, such as grouped by geographic region. For example, distribution network104may have any topology known in the art, such as transmission, distribution, feeder, and neighborhood levels. DERMS110may include a plurality of aggregators112for controlling the DERs connected to specific portions of the grid. Similarly, reference signals SREF may include reference signals for specific sections of the grid so that the aggregators112can separately control DERs102in different geographic regions so that the net energy load of the DERs located or otherwise electrically proximate each region tracks the corresponding grid reference signals for each region. In some examples, an aggregator112may be configured to control multiple regions of DERs102by comparing the requests for a specified region of DERs to a corresponding region-specific reference signal.

FIG.2is a functional block diagram of an example instantiation of one of aggregators112and one of DERs102. DER102can be virtually any type or electrical device including any of the types of DERs described herein. DER controller108may be communicatively coupled to DER102in any way, such as physically coupled to a corresponding DER, physically near or adjacent a DER and in wired or wireless communication, for example, located in the same building and connected on a local network, or implemented remotely as a cloud-based software as a service for remotely managing one or more DERs102.

DER102is configured to repeatedly consume and/or discharge energy in discrete quantities referred to herein as energy packets by drawing or discharging energy for a fixed period of time referred to herein as a control epoch202, also referred to herein as the energy packet time duration, which is stored in memory204. In one example, when a DER102begins to draw power, a local packet timer206is initiated and when the timer reaches the control epoch202value, DER102ceases drawing or discharging power. DER102is configured to repeatedly determine whether to draw or discharge power at regular intervals referred to herein as a communication epoch208, also referred to herein as the communication interval. The control epoch202and communication epoch208may be the same, in which case DER102may immediately determine whether to request another packet after a preceding packet ends. Or the communication epoch208may be greater or less than the control epoch202, and the communication epoch may vary according to whether the DER102previously determined to not request an energy packet or if aggregator112granted or denied the DER's prior request. In other examples, rather than requesting energy packets from an aggregator, DER102may be configured to repeatedly draw or discharge energy packets according to a duty cycle210that includes periods where the device is on followed by periods of time the device is off. More generally, aspects of the present disclosure may be applied to any DER in communication with an aggregator that communicates at least a start time to the aggregator when it transitions from an OFF or STANDBY state to an ACTIVE or ON state and that is configured to interrupt its active state and transition back to OFF or STANDBY in response to a timer value exceeding a calculated timer threshold due to a measured local frequency falling below a threshold value.

In some examples, one or both of the communication epoch208and control epoch202are configurable by a user of a particular DER, by aggregator112and/or by grid operator114. In some examples, memory204includes a plurality of values for communication epoch208and/or control epoch202and a processor212executing instructions according to a PEM application214stored in memory204may be configured to select a particular one of the plurality of values for the communication epoch and control epoch. Aggregator112may similarly include a plurality of values for communication epoch208and control epoch202stored in memory216and be configured to send instructions to one or more DERs102to change one or both of the communication epoch and control epoch.

Memory204of DER controller108may also store one or more probability profiles218that define a probability a DER102will make a request to consume or discharge power during a communication epoch. DER controller108and/or DER102may also include or be communicatively coupled with one or more sensors220for measuring one or more locally-sensed conditions, for example, sensors220may include one or more of a temperature sensor, a pressure sensor, a current sensor, a frequency sensor, among others and be configured to generate a signal indicating a locally measured real time temperature, pressure, revolution rate, state of charge, power line frequency or any other condition. Any type of sensor known in the art may be used. For example, an energy metering microchip may be used, such as Microchip Technology Inc.'s energy metering microchip model number ATM90E26, or other microchips or other products with frequency measurement capability.

PEM application214may include instructions for selecting a probability profile218and determining a probability value for a given communication epoch from, for example, the selected probability profile226and a locally-sensed condition from one or more sensors220. The probability value may be, for example, a probability that DER102will make a request to aggregator112to consume or discharge power during a given communication epoch. In other examples, DER102does not request permission from aggregator112and the probability value may be a probability that DER102will decide to consume or discharge power during a given communication epoch.

By way of non-limiting example of an implementation of such a probability value, the determined probability value may be a given value over a predetermined range, such as between zero and 1 and processor212may be configured to execute a random number generator that randomly generates a number between zero and 1. PEM application214may include instructions to request to consume or discharge energy when the random number generator value is less than the probability value and not make a request when the random number generator value is greater than the probability value. As will be appreciated by persons of ordinary skill in the art, other implementations may also be used.

DER controller108may also store a FFR application240in memory204that includes instructions for interrupting the DER's consumption or discharge of an energy packet before the packet timer206has expired when a deviation in the grid frequency measured by sensor220from a setpoint value and/or a ROCOF exceeds a threshold value. In one example, FFR application240may include instructions for interrupting an energy packet according to the following:

Interrupt⁢my⁢packet?={Yesif⁢ηn,[k]≥η⁡(Δ⁢f[k])NoelseEqn.(1)Wherein:ηn[k]=tn[k]δ∈[0,1]Eqn.(2)Eqn.(3)η⁡(Δ⁢f)={1,if⁢Δ⁢fdb<Δ⁢f[k]1-(Δ⁢fdb-Δ⁢f[k]Δ⁢fdb-Δ⁢fmax)⁢(1-ηmin),if⁢Δ⁢fmax≤Δ⁢f[k]≤Δ⁢fdbηmin,if⁢Δ⁢f[k]≤Δ⁢fmax.ηn[k] is a normalized timer state;ηn(Δf[k]) is a packet interruption threshold that is dependent on the locally measured AC powerline frequency measured locally frequency at time step k;tn[k] is the value of packet timer206at time step k;δ is the control epoch length or other packet length, such as the length of time a DER is drawing or discharging power according to a duty cycle210;Δf[k] is the absolute value of the deviation of the real time AC powerline frequency measured locally, such as at or proximate a DER from the frequency setpoint, which is typically 60 Hz in the United States;Δfabis the minimum deviation of the measured frequency from the target frequency (e.g., the upper part of the deadband for a frequency drop);Δfmaxis a maximum deviation of the measured frequency from the target frequency (e.g., the lower part of the deadband for a frequency drop); andηminis a minimum packet timer threshold value and is less than one (1) and greater than or equal to zero (0).

The above method in Eqn. (3) can be generalized to the following:

η⁡(Δ⁢f)={1,if⁢Δ⁢fdb<Δ⁢f[k]1-g⁢(Δ⁢f[k])⁢(1-ηmin),if⁢Δ⁢fmax≤Δ⁢f[k]≤Δ⁢fdbηmin,if⁢Δ⁢f[k]≤Δ⁢fmax.Eqn.(4)

Where g(Δf [k]) is a non-increasing function of frequency deviation (e.g., drop) that maps frequency deviation (from minimum deviation Δfabto maximum deviation Δfmax) to timer threshold (from 1.0 to ηmin) and is 0.0 when deviation is less than Δfaband 1.0 when deviation is more than Δfmax.

Clearly, the case of frequency rise (or over-frequency event) can be treated similarly by interrupting discharging (or generating) packets for DERs, which follows as above for positive (increasing) frequency deviation.

The foregoing equations 1-4 provide a dynamic packet interruption threshold, ηn(Δf [k]), that is a function of the locally measured frequency and is designed to be compared to the packet timer state of a single device. The result is that DERs102that have drawn or discharged more power because they are closer to the end of an energy packet or have otherwise been ON for a longer period of time, are interrupted first. This may be preferable to simply interrupting all DERs that are drawing or discharging power during a frequency event because, depending on the extent of participating DERs in the system, a sudden interruption of all of the DERs rather than just the ones with higher timer values could be too abrupt and cause a system instability. In addition, the controlled interruption of DERs provides exactly the mechanism to design, optimize, and quantity the MW/Hz sensitivities that define (synthetic) damping useful for FFR markets and grid operators. Additionally, prioritizing the interruption of DERs with greater packet timer values also supports quality of service because those DERs are more likely to have a lower immediate need for energy due to being at the later stage of an energy packet and/or having been on for a relatively longer period of time.

FIG.3graphically illustrates Equation 3 with a deviation of a locally measured frequency from a setpoint frequency, Δf(Hz), along the x axis and the packet timer interruption threshold, ηn(Δf[k]), along the y axis.FIG.3illustrates how the packet timer interruption threshold varies with measured frequency for both over frequency (Δf>0) and under frequency (Δf<0) events. In this example, the dead band frequency, Δfab, is defined. When the measured deviation of the locally measured frequency is within the dead band frequency, no packets are interrupted. When the deviation of the measured frequency exceeds the dead band frequency, the normalized packet timer interruption threshold, ηn(Δf[k]), falls below 1. In the illustrated example, the packet interruption threshold is defined as a linear equation with the packet timer interruption threshold decreasing linearly with respect to increasing frequency deviation, with the packet timer interruption threshold finally reaching a minimum value when the deviation in the measured frequency from the frequency setpoint reaching the maximum frequency threshold value, Δfmax. In other implementations, it can be a linear or nonlinear function of either frequency deviation or the rate of change in frequency deviation or any combination of the two.

In the illustrated example, the packet timer interruption threshold, ηn(Δf[k]) is symmetric about the y axis, with the same frequency dead band Δfab, maximum frequency deviation Δfmax, and slope of the threshold curve between the dead band and maximum values being the same for both under frequency and over frequency events. In other examples, each of the foregoing values may be different for under and over frequency.

In the illustrated example, the packet timer interruption threshold is a function of the deviation of the locally sensed frequency from the target frequency. In other implementations, the packet timer interruption threshold may be a function of other parameters as well, such as the DER's need for energy, and/or one or more other locally-sensed values for a particular DER that may indicate the DER's need for energy. For example, in the case of thermostatically controlled loads (TCL), the packet timer interruption threshold may be determined as a function of a deviation of a controlled temperature from a target temperature, such as a difference between a current temperature of a water heater, refrigerator, or room being heated or cooled by an HVAC system from the setpoint temperature. A battery's state of charge or other measure of the charge level of the battery may be used to determine the packet timer interruption threshold for a battery. In one example, the packets of DERs with a higher packet timer state and lower need for energy according to a locally sensed condition such as temperature or SOC may be prioritized for interruption prior to DERs with a higher need for energy that have the same packet timer state.

During an under-frequency event, DERs102may provide synthetic inertia from distribution networks or microgrids104by dynamically interrupting the energy packets of select proportions of the DERs that are drawing power from the grid in an active energy packet when the under-frequency event occurs. The DERs can be interrupted according to a local control law, such as a local control laws described herein, including interruption of packets according to timer state. For bidirectional DERs that are capable of both consuming power from the grid as well as discharging power to the grid, such as batteries, during an under-frequency event, any bidirectional DER that is drawing power and is interrupted may be controlled to immediately begin discharging power to the grid to further increase the synthetic damping provided by the DERs. Other bidirectional DERs that were not drawing power can also be controlled to immediately discharge power to the grid.

During an over-frequency event, DERs102may provide synthetic inertia to distribution network104by dynamically interrupting the energy packets of select proportions of the DERs that are discharging power to the grid in an active energy packet when the over-frequency event occurs. The DERs can be interrupted according to a local control law, such as a local control laws described herein, including interruption of packets according to timer state. For bidirectional DERs that are capable of both drawing power from the grid as well as discharging power to the grid, such as batteries, during an over frequency event, any bidirectional DER that is discharging power and is interrupted may be controlled to immediately begin drawing power to the grid to further increase the synthetic damping provided by the DERs. Other DERs, including bidirectional DERs that were not discharging power, can also be controlled to immediately begin drawing power from the grid.

For bi-directional DERs, by prioritizing for interruption the DERs in an active energy packet that have a higher timer state and are, therefore, nearing the end of the energy packet, and flipping those bi-directional DERs from drawing power to discharging power or vice versa, the DERs that will most likely provide the most damping will be interrupted and flipped first. This is because a discharging DER nearing the end of a energy packet will have recently discharged energy and therefore be more likely to be capable of drawing energy, and a charging DER nearing the end of a energy packet will have recently drawn and stored energy and therefore be more likely to be capable of discharging that recently stored energy.

FIG.4illustrates an example implementation of a FFR local control law as applied to 400,000 DERs102in the form of water heaters during an under-frequency event. Packet timer values are shown along the x axis and a locally sensed condition, here the temperature of the water in the water heaters, along the y axis, with each circle on the plot representing one of the 400,000 DERs. In this example, the control epoch (time duration of an energy packet) is three minutes. As would be expected, there is a general upward trend in water temperature as the value of the packet timer increases because the water heaters with a higher packet timer have been drawing energy to heat water for a longer period of time. At the point in time shown inFIG.4, the frequency has dropped below the dead band frequency such that the timer interruption threshold408is less than the energy packet length, resulting in all DERs with a packet timer value greater than timer threshold408being interrupted (interrupted packets402). If the measured frequency continues to fall, the timer interruption threshold408will move to the left and decrease resulting in a greater proportion of water heaters being interrupted.

Referring again toFIG.2, memory204of DER controller108may be configured to also store FFR parameters242in memory, which may include any of the FFR parameters disclosed herein including ΔfabΔfmaxand ηmin, among others. The foregoing parameters may be the same for both under-frequency and over-frequency or the memory may include one or more separate FFR parameters for over and under frequency events.

Memory216of aggregator112may include one or more software applications250including a PEM application252that includes instructions for controlling DERs102according to a PEM control scheme including any of the PEM control schemes disclosed herein. In embodiments where DERs102make a request to aggregator112to consume or discharge a energy packet, PEM application252may include instructions for approving or denying the request and also include instructions for tracking the accepted requests in memory216.

With the limited information of (1) the FFR local control law the DERs102are configured to execute, (2) the number of DER energy packets that begin at each time step, and (3) the time duration of the packets (also referred to herein in some examples as the control epoch), aggregator112can utilize damping and inertia calculation application256to calculate a real time estimate of the available synthetic damping (MW/s) and/or inertia (MW/s{circumflex over ( )}2) that is available from the fleet of DERs. Aggregator112may also be configured to provide real-time power-Frequency frequency-response capability curves for the fleet of PEM-enabled DERs and create a power deviation vs. frequency response map that can be used by grid operator114to understand how much frequency response reserves are available from the DER fleet. Aggregator112may also include a grid services bid application258which can include instructions for using the calculation of damping and inertia of the fleet in an optimization calculation to determine how much to bid day-ahead for primary frequency response services, which could also include consideration of expected participation in other grid services.

In one example, aggregator112is configured to maintain a timer histogram254in memory216that maintains a record of a number of DERs that are consuming or discharging a energy packet at any given time and the normalized packet timer state of each active DER. In one example, the normalized packet timer state is the packet timer value divided by the time duration of the energy packet. By tracking a normalized value, DERs102configured to operate with different-length energy packets can be considered together as a group. Referring again toFIG.4,FIG.4includes a conceptual graphical representation404of packet timer histogram254. In the illustrated example, aggregator112may be configured to break up a normalized control epoch into a plurality of bins406(only one labeled) and track the number of DERs102in each bin. With knowledge of the FFR local control law the DERs are configured to execute, including the packet timer interruption threshold value calculated by each DER, aggregator112can predict the number of active DERs that would automatically, independently, and instantaneously interrupt an active energy packet at any given point in time in response to a local power line frequency measurement.FIG.4also illustrates a temperature histogram410that maintains a status of temperatures of the fleet of water heaters. In some examples, the water heaters may be configured to communicate a water temperature to an aggregator while in others the aggregator may include instructions for calculating an estimate of the temperature histogram according to knowledge of the fleet of DERs and history of accepted or denied energy packets. Aggregator112may similarly be configured to maintain a SOC histogram as well as or alternatively a real time energy capacity distribution that quantifies a level of energy stored in by the fleet of DERs at any given point in time.

Aggregator112may be configured to track a packet timer status of DERs operating under a variety of different PEM control schemes. In some examples, one or more DERs102may be configured to independently determine when to consume power without requesting permission from aggregator112. Such DERs can send a notification to aggregator112indicating when it began a new energy packet and the time duration of the energy packet, which the aggregator can use to maintain the timer histogram254. The foregoing may include DERs configured to draw power according to a duty cycle, where the period of time in the duty cycle the DER is ON is considered a energy packet. In yet other examples, DERs that are not operating according to a PEM control scheme or a regular duty cycle may have a minimum ON time duration that may be used in place of a control epoch. In such cases, an aggregator may be configured with a default minimum ON time, which may be defined for a particular category of DER (such as air conditioner, water heater, battery charger, etc.) or the minimum ON time may be communicated from the DER to the aggregator. The DER may use the minimum on time to calculate a timer interruption threshold value. The DER102may be equipped with communication capability for transmitting a status signal to the aggregator for tracking a damping and inertia of the fleet.

Referring again toFIG.2, damping and inertia calculation application256may include instructions for calculating a real time estimate of an equivalent damping and inertia available from the fleet of DERs102according to the timer histogram254. In one example, the damping and inertia calculation application256includes instructions for performing a calculation according to the following:

⇒DPEM=Prate⁢ηmax⁢np⁢x_nomΔ⁢fdb-Δ⁢fmaxEqn.(5)

Wherein:DPEMis the damping coefficient (e.g., MW/mHz) for the fleet and describes the sensitivity of the fleet's packet interruption rate (change in power) as a function of change in frequency,Prateis power rating (kW) of each device and given by Aggregator,npis total number of timer bins based on resolution of sensors and desired response characteristics,Xnomis the average number of DERs that are ON in each bin.

Eqn. (5) provides a highly accurate estimate of the fleet's ability to respond to changes in frequency (MW/mHz) when the distribution of timers across the fleet is approximately uniform. If the fleet is participating in other grid services or the Aggregator(s) have otherwise been varying the control of the fleet an alternate formulation described in equations 6-8 described below may be used. When the fleet is in a stationary state where the distribution of timers is approximately uniform, the rate of response (change in MW vs. change in frequency) is constant, which explains why DPEMin Equation 5 is constant.

In the case where DERs are actively participating in other grid services (e.g., ancillary services or load shaping) and not in a stationary state, Eqn. (5) can be generalized to provide grid operators or a DER aggregator a real-time, power-frequency response map:

PPEM(Δ⁢f[k])=∑b=0β⁡(Δ⁢f[k])-1tB-b[k]⁢PB-brate[k]Eqn.(6)

Wherein:PPEMis the cumulative power reduction at time k achieved by the fleet for a given change in frequency Δf[k] between Δfmaxand Δfaband is 0 above Δfdb, wheretB-b[k] is the number of packets in coordinator's timer distribution bin B-b at time k, where B=δ/Δt is the total number of bins (and largest bin index corresponding with normalized timer states with value 1.0) in the coordinator's timer distribution and δ is the packet duration (seconds) and Δt is the timer distribution's bin width (seconds),PrateB-b[k] is the average height (power rating in kW) of charging packets in the coordinator's timer distribution bin B-b at time k, andB−(β(Δf[k])−1) is the smallest (last) bin in the timer interrupted by frequency drop (at nadir) and β(Δf[k]) is defined by the timer threshold and designates the bin to interrupt and is computed by the coordinator with the following:

β⁡(Δ⁢f[k])=δΔ⁢t⁢η⁡(Δ⁢f[k])Eqn.(7)

Wherein:β(Δf[k]) in Eqn (6) represents the bin number corresponding to the computed timer threshold. Rounding techniques may be used, if necessary, to match β with specific bin.

If the frequency exceeds Δfmaxthe power drop of the fleet saturates with the timer threshold at

PPEM(Δ⁢f[k])=∑b=0β⁡(Δ⁢fmax)-1tB-b[k]⁢PB-brate[k].Eqn.(8)

FIGS.5A-6Bconceptually illustrate a timer distribution and corresponding power deviation vs. frequency response map and associated damping coefficient for a fleet with a uniform timer distribution (FIGS.5A,5B) and for a fleet with a non-uniform timer distribution (FIGS.6A,6B).FIG.5Ashows a uniform timer distribution500which can be expected when a fleet of DERs is operating in steady state and following a steady control signal from an aggregator. The X axis of timer distribution500is a normalized control epoch divided into a plurality of bins502(only one labeled) and the y axis represents a height of each bin which represents a number of participating DERs having a timer value that falls within a corresponding bin. The choice of bin width is a design variable chosen by the aggregator and affects quantization of the estimate. It should not be chosen too large as that causes the Power-frequency curve504to become jagged. In some example applications, a bin width is about 5-10 second may be appropriate, with the optimum bin width depending on the size of fleet. For example, large fleets can employ smaller bin widths as more devices are participating and thus populating the timer bins502.

FIG.5Bis a power deviation vs. frequency response map504that illustrates a change in power (y axis) than can be expected from interrupting participating DERs as a function of a potential or expected deviation of frequency (x axis), with the slope of the curve equal to the damping. An aggregator, such as one of aggregators112or another computing device configured to perform a calculation with equation 5, may be used to generate a real time power deviation vs. frequency response map and synthetic damping illustrated inFIG.5B.

FIG.6Aconceptually illustrates an example of a non-uniform timer distribution600which can be expected when a fleet of DERs has not been operating in steady state and has been following a varying control signal from an aggregator, which may occur when the fleet of DERs are being used to provide additional ancillary services and/or after a grid disturbance event. The X axis of timer distribution600is a normalized control epoch divided into a plurality of bins602(only one labeled) and the y axis represents a number of participating DERs having a timer value that falls within a corresponding bin.FIG.6Bis a power deviation vs. frequency response map604that illustrates a change in power (y axis) than can be expected from interrupting participating DERs as a function of a potential or expected deviation of frequency (x axis), with the slope of the curve equal to the damping. An aggregator, such as one of aggregators112or another computing device configured to perform a calculation with equations 6-8 may be used to generate a real time power deviation vs. frequency response map and synthetic damping illustrated inFIG.6B. As conceptually illustrated inFIGS.6A and6B, the power deviation vs. frequency response map604can be expected to be non-linear and the resulting damping can be variable as a function of frequency deviation when the distribution of the fleet is non-uniform.

In some examples, aggregator112may be configured to broadcast updated values for the FFR parameters242to the fleet of DERs102at any time and in response to any of a variety of conditions. In one example, aggregator112may be configured to update one or more FFR parameter242according to a number of DERs that are active and the proportion of the total load of the DERs on the distribution network in order to optimize grid stability and availability of damping and inertia. For example, the minimum packet timer value, ηmin, the frequency dead band, Δfdb, and/or the slope of the curve of the packet interruption threshold, ηn(Δf [k]), may be adjusted.

Example System Architectures

System100may have any of a variety of architectures including various combinations of hardware and software to provide the various functionalities described herein. For example, processor234of each aggregator112, and processor212of each DER controller108may include any type of one or more processing units known in the art, including one or more data processors, image processors and/or central processing units, and can be separate components, or can be integrated in one or more integrated circuits, and can be coupled by one or more communication buses or signal lines. Memory216of aggregators112and memory204of DER controller108may similarly be of any of a variety of types and configuration and may be a non-transitory machine-readable storage media. Memory of various components in system100, such as memory216or204, may include high-speed random access memory and/or non-volatile memory, such as DRAM, SRAM, EDO RAM, one or more magnetic disk storage devices, one or more optical storage devices, and/or flash memory. Memory216and memory204can store a number of other applications and instructions in addition to those illustrated and described, such as operating system instructions and instructions for operating any hardware and peripherals operably coupled to the corresponding memory and computing device.

Various components of system100may be configured to exchange data and information via any means known in the art, including communication between aggregator112and other components of the system via communication interface236and communication between each DER102and other components of the system via communication interface260. Communication interfaces236and260may include wireless communication components, which can include antennas, radio frequency receivers and transmitters and/or optical (e.g., infrared) receivers and transmitters. The specific design and implementation can depend on the communication network(s) over which the particular computing device, e.g. particular aggregator112or DER102is intended to operate. By way of example, communication interfaces can include communication components designed to operate over a GSM network, a GPRS network, an EDGE network, a Wi-Fi or WiMax network, and a Bluetooth™ network.

Communication interfaces of system100, such as communication interfaces236and260, may include one or more network devices configured to communicate with one or more networks, for example, Local Area Network (LAN), Wide Area Network (WAN) or the Internet through a variety of connections including, but not limited to, standard telephone lines, LAN or WAN links (for example, 802.11, T1, T3, 56 kb, X.25), broadband connections (for example, ISDN, Frame Relay, ATM), wireless connections, or controller area network (CAN). Communication interfaces may also include a built-in network adapter, network interface card, PCMCIA network card, card bus network adapter, wireless network adapter, USB network adapter, modem or any other device suitable for interfacing a computing device of system100to any type of network capable of communication and performing the operations described herein. And the components of system100may be communicatively coupled by any means known in the art, e.g., local and/or remote communication utilizing any of a wide range of wired and/or wireless communications protocols, including, for example, a digital multiplexer (DMX) interface protocol; a Wi-Fi protocol; a Bluetooth protocol; a digital addressable lighting interface (DALI) protocol; a ZigBee protocol; a near field communication (NFC) protocol; a local area network (LAN)-based communication protocol; a cellular-based communication protocol; an Internet-based communication protocol; a satellite-based communication protocol; and/or a combination of any one or more thereof.

Aggregator112and DER controller108may be implemented in any type of computer system, such as a workstation, desktop computer, server, laptop, handheld computer, tablet computer, mobile computing device, server, one or more virtual machines, or other form of computing or telecommunications device that is capable of communication, has sufficient processor power and memory capacity to perform the operations described herein.

Example Methods of Operation

FIG.7illustrates an example method700of providing a frequency response to an electrical grid with a DER, such as one of DERs102, and may include, at block701, drawing or discharging energy packets, for example, according to any of the packetized energy control schemes described herein, and at block703, maintaining a timer that indicates a duration of time since the commencement of drawing or discharging energy or amount of time until a planned termination of the drawing or discharging of energy. In some examples, block703includes maintaining a packet timer that indicates a remaining amount of time before the completion of an energy packet. At block705, the method may include calculating a timer interruption threshold according to a measured frequency of the electrical grid. In one example, block705may include execution of FFR application240including performing calculations according to Equations 1-3 above. At block707, the method may include comparing the timer value to the calculated timer interruption threshold, and at block709, interrupting the drawing or discharging of energy in response to the timer exceeding the calculated timer interruption threshold. For example, a locally measured AC frequency of the electrical supply adjacent a DER may deviate from a frequency set point, e.g., 60 Hz, by an amount greater than a frequency dead band value, resulting in a decrease in the calculated timer interruption threshold such that the timer value suddenly exceeds the calculated timer interruption threshold. DER102may then independently and immediately interrupt the drawing or discharging of energy as part of a FFR local control policy.

FIG.8illustrates an example method800of providing frequency response services to an electrical grid with an aggregator, such as aggregator112in communication with a plurality of DERs, such as a plurality of DERs102. At block801, the method may include recording a start time when each DER begins to draw or discharge electrical energy from or to the electrical grid and at block803, maintaining a timer histogram, such as energy packet timer histogram254, that provides a status of timer states for the fleet of DERs. At block805, the method may include estimating an available synthetic damping (MW/s) and/or inertia (MW/s{circumflex over ( )}2) according to the timer histogram. In one example, block805may include execution of damping and inertia calculation application256, including performing calculations according to Equations 5-8 above. And at block807, providing frequency response services forecasting and bids to a utility according to the estimating.

The foregoing has been a detailed description of illustrative embodiments of the disclosure. It is noted that in the present specification and claims appended hereto, conjunctive language such as is used in the phrases “at least one of X, Y and Z” and “one or more of X, Y, and Z,” unless specifically stated or indicated otherwise, shall be taken to mean that each item in the conjunctive list can be present in any number exclusive of every other item in the list or in any number in combination with any or all other item(s) in the conjunctive list, each of which may also be present in any number. Applying this general rule, the conjunctive phrases in the foregoing examples in which the conjunctive list consists of X, Y, and Z shall each encompass: one or more of X; one or more of Y; one or more of Z; one or more of X and one or more of Y; one or more of Y and one or more of Z; one or more of X and one or more of Z; and one or more of X, one or more of Y and one or more of Z.

Various modifications and additions can be made without departing from the spirit and scope of this disclosure. Features of each of the various embodiments described above may be combined with features of other described embodiments as appropriate in order to provide a multiplicity of feature combinations in associated new embodiments. Furthermore, while the foregoing describes a number of separate embodiments, what has been described herein is merely illustrative of the application of the principles of the present disclosure. Additionally, although particular methods herein may be illustrated and/or described as being performed in a specific order, the ordering is highly variable within ordinary skill to achieve aspects of the present disclosure. Accordingly, this description is meant to be taken only by way of example, and not to otherwise limit the scope of this disclosure.