Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t)...

35
Smart refrigerators: a distribution-referred approach to decentralised control Simon Tindemans, Vincenzo Trovato, Michael Evans, Goran Strbac [email protected] Probabilistic methods for energy networks King’s College London, 14 March 2017

Transcript of Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t)...

Page 1: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Smartrefrigerators:adistribution-referredapproachtodecentralisedcontrol

SimonTindemans,VincenzoTrovato,MichaelEvans,GoranStrbac

[email protected]

ProbabilisticmethodsforenergynetworksKing’sCollegeLondon,14March2017

Page 2: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Outline

1. Opportunitiesforflexibleloadcontrol2. Distribution-referredcontrol3. Applicationtorefrigerators4. Systemusesofaggregatecontrol

Page 3: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Flexibilityinthefuturepowersystem

electricitygridin=out

CO2CO2

CO2

+-flexibility

+flexibility+flexibility

Page 4: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Theflexibilityspectrum

Flexiblegenerators

Gridscalestorage

IndustrialandcommercialDR

+-

residentialDR

10-100MW 1-1000MW 100kW- 1MW 10W– 1kW

source

individualmagnitude

number 100 10- 1000 1000- 10000 10sofmillions

‘Longtail’ofdemandresponse

Page 5: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Unlocking the longtail

Constraints/challenges• Flexibilityisaby-productofotheractivities,whichwill

affectserviceavailability.• Smallper-deviceflexibilitycontribution,socostand

attentionbudgetissmall• Privacyconcerns• Significantheterogeneity

Opportunity• Verylargenumberofdevices,sostatisticaldiversity

canbeanasset.

Page 6: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Flexiblerefrigeration

1000 2000 3000 4000 5000 6000 7000t [s]

20

40

60

80P [W]

0

Pon

P

0 1000 2000 3000 4000 5000 6000 7000t [s]

2

4

6

8

T [°C]

Tmin

Tmax

T

sourceofflexibility

Powerconsumption

Temperature

≅( )

Page 7: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Flexiblerefrigeration:from‘what’to‘how’

Theopportunity• Refrigeratorsrepresent5-15%of

systemload(est.2-3GWinGB)*• Loadshiftingfor~30minutesis

free*secondaryuse

0

20000

40000

60000

potentialzerocostflexibility!

Thechallenge• Maintaincoolingperformance: Secondaryuse(flexibility)shouldnot

compromisetheprimaryuse(cooling)ofdevices.• Robustnessandscalability:Relianceonreal-timecommunication

mayresultinbottlenecksandsinglepointsoffailure• Controllability:Ensuresufficientcontroloverpowerconsumption,

andavoidunwantedinteractions.

Page 8: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

DISTRIBUTION-REFERREDCONTROL

Page 9: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Controlrequirements

Aim:controltheaggregatepowerconsumptionofacollectionofappliances𝒜:𝑃$%$&' 𝑡 = ∑ 𝑃&(𝑡)�

&∈𝒜

Constraints:• Useonly aglobalcontrolsignalΠ(𝑡),suchthat𝑃$%$&' 𝑡 ≅ Π 𝑡 𝑃0$%$&'

• Applianceshavelimitedcontrollability• Discretepowerlevels• Localqualityofservice

𝑃$%$&' 𝑡 ≅1𝑃&(𝑡)�

&

= 𝑃0$%$&'×Π 𝑡Π(𝑡)

Page 10: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Complexpopulationconstraints

thermalmodels

states

(temperature,on/off*)

shouldnotswitchoff

shouldnotswitchon

x

x

x

Page 11: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Distribution-referredcontrol

𝐸 𝑃$%$&' 𝑡 = 𝑁𝐸 𝑃 𝑡 |Π7$

= 𝑁𝐸8 𝐸9(8) 𝑃(𝑆) 𝑀, Π7$ lawoftotalexpectation

= 𝑁 11𝑁𝐸9(>

?) 𝑃(𝑆) 𝑚&, Π7$�

&∈𝒜

= 1 𝐸9(>?) 𝑃(𝑆) 𝑚&, Π7$�

&∈𝒜

Assumptions• State𝑠 ofafridgewithmodel𝑀 isarandomrealisationof𝑆(𝑀)• ThedistributionofSissufficientlydiversetopermitvariationin𝑃(𝑆)• Independencebetweenappliances(O 𝑁CDE convergence)

𝐸9(>?) 𝑃(𝑆) 𝑚&, Π7$ = 𝑃0&Π 𝑡

= 1 𝑃0&Π 𝑡�

&∈𝒜

Controlconstraint

= 𝑃0$%$&'Π 𝑡

Page 12: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Distribution-referredcontrol

thermalmodels

states

(temperature,on/off*)

x

Howcanweoperatetheaggregatewithoutreal-timecommunication?

designcontrolstrategybasedonall statesforaspecificmodel

x

x

Page 13: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

DISTRIBUTION-REFERREDCONTROLFORREFRIGERATORS

• Decentralisedcontrolofthermostaticloadsforflexibledemandresponse.SimonTindemans,VincenzoTrovato,GoranStrbacIEEETransactionsonControlSystemsTechnology(2015)

• Nondisruptive decentralizedcontrolofthermalloadswithsecondorderthermalmodelsSimonTindemans,GoranStrbac2016IEEEPESGeneralMeeting,Boston.

Page 14: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Constructingapopulation

𝑑𝑇(𝑡)𝑑𝑡 = H −𝛼(𝑇(𝑡) − 𝑇%K)

−𝛼(𝑇(𝑡) − 𝑇&>LMNK$)(on)(off)

0 1000 2000 3000 4000 5000 6000 7000t [s]

2

4

6

8

T [°C]

Tmin

Tmax

T

cabinet𝑇R(𝑡)

[on/off]

𝑇S>LMNK$environment

Refrigeratordynamics

0 2 4 6 8T [°C]

0.05

0.10

0.15

0.20

0.25

0.30probability

Tmin TmaxT

Steadystatetemperaturedistribution

Page 15: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Distribution-referredcontrol

Chooseone-parameterfamilyoftransformations:• Compatiblewithtemperaturelimits• Sufficientlyslowtoberealisable• Distributionknownatalltimes

𝑣 𝑇, 𝑡 = 𝛼𝛽(𝑡)(𝑇 − 𝑇>&V) 0 2 4 6 8T [°C]

0.05

0.10

0.15

0.20

0.25

0.30probability

Tmin TmaxT

𝑑𝑇W(𝑡)𝑑𝑡 = −𝛼 𝑇W 𝑡 − 𝑇%XX + 𝑇%XX − 𝑇W0 Π(𝑡)

Relateaveragetemperatureandpower

DesiredpowerconsumptionΠ 𝑡

Populationaveragetemperature𝑇W 𝑡

one-to-oneTransformationofpopulationdensity(parametrisedby

𝑣 𝑇, 𝑡 )

one-to-many

continuityequation𝜕𝜕𝑡 𝑓 𝑇, 𝑡 = −

𝜕𝜕𝑇 𝑣 𝑇, 𝑡 𝑓 𝑇, 𝑡

Page 16: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Computeswitchingactions

2 3 4 5 6 7 8T [°C]

0.05

0.10

0.15

0.20

2 3 4 5 6 7 8T [°C]

0.05

2 3 4 5 6 7 8T [°C]

0.05

0.10

0.15

on

off

-4000 -2000 0 2000t [s]

2

4

6

8

T [°C]

Tlow

Thigh

θ

stochastic switching event

desiredcollectiveresponse

device-levelactions

2 3 4 5 6 7 8T [°C]

0.05

0.10

0.15

0.20

requiredflux

computeswitchingactions1. Deterministicswitching2. Stochasticswitching

Π(𝑡)

Page 17: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Controllerimplementation(summary)

2 3 4 5 6 7 8T [°C]

0.05

0.10

0.15

0.20

-4000 -2000 0 2000t [s]

2

4

6

8

T [°C]

Tlow

Thigh

θ

stochastic switching event

2. Constructahomogeneous ’virtualpopulation’withrandomtemperatures.

4. Determinedevice-specificactions,basedontheactualdevicetemperature

2 3 4 5 6 7 8T [°C]

0.05

0.10

0.15

0.20

3. Manipulatethe‘virtualpopulation’tocontrolits(virtual)powerconsumptioninlinewithΠ 𝑡 .

𝑑𝑇(𝑡)𝑑𝑡 = H −𝛼(𝑇(𝑡) − 𝑇%K)

−𝛼(𝑇(𝑡) − 𝑇&>LMNK$)(on)(off)

1. Eachapplianceknowsitsstate andmodel

Π(𝑡)

Page 18: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Controllerimplementation(summary)

2 3 4 5 6 7 8T [°C]

0.05

0.10

0.15

0.20

-4000 -2000 0 2000t [s]

2

4

6

8

T [°C]

Tlow

Thigh

θ

stochastic switching event

2. Constructahomogeneous ’virtualpopulation’withrandomtemperatures.

4. Determinedevice-specificactions,basedontheactualdevicetemperature

2 3 4 5 6 7 8T [°C]

0.05

0.10

0.15

0.20

3. Manipulatethe‘virtualpopulation’tocontrolits(virtual)powerconsumptioninlinewithΠ 𝑡 .

𝑑𝑇(𝑡)𝑑𝑡 = H −𝛼(𝑇(𝑡) − 𝑇%K)

−𝛼(𝑇(𝑡) − 𝑇&>LMNK$)(on)(off)

1. Eachapplianceknowsitsstate andmodel

Π(𝑡)Eachapplianceconsidersitselfasarandom

representativeofapopulation...

...andtakesactionsinlinewithpopulationobjectives

Page 19: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Aggregateresponse

0 1000 2000 3000t

Π(t)

Collectively,heterogeneousfridgestrackthereferencesignal𝛱 𝑡

N=1000domesticrefrigerators

Page 20: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Nextsteps:labtesting

• ModifyfridgewithtemperaturesensorsandRaspberryPitocontrolthecompressor

• ...• Profit!

Challenges• Howmanytemperaturemeasurementstoinclude?• Sensitivitytomodelmisspecification• Howtotestastochasticcontrollerinasmall-scalelabtrial?

Page 21: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

AGGREGATESERVICEPROVISION

• TheLeakyStorageModelforoptimalmulti-serviceallocationofthermostaticloads.VincenzoTrovato,SimonTindemans,GoranStrbacIETGeneration,Transmission&Distribution(2016)

• UnderstandingtheAggregateFlexibilityofThermostaticallyControlledLoadsVincenzoTrovato,SimonTindemans,GoranStrbacPowerTech 2017(submitted)

Page 22: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

High-levelapproach

Aggregator

AggregateloaddispatchModelofcollectiveflexibility

controlsigna

l

resourceavailability

flexibilityproducts

Markets

Page 23: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Theleakystorageunit

Six-parametermodeltodescribetheflexibilityofahomogeneouspopulation

𝑑𝑆(𝑡)𝑑𝑡 = −

𝑆 𝑡𝜏 + 𝑃(𝑡)

withconstraints:

𝑃>MK ≤ 𝑃 𝑡 ≤ 𝑃>&V𝑆>MK ≤ 𝑆(𝑡) ≤ 𝑆>&V

_ 𝑆 𝑡 𝑑𝑡 = 𝑆0`

0

P(t)

αS(t)

Smax

S0Smin

S(t)

preservethefood!

𝑑𝑇W(𝑡)𝑑𝑡 = −𝛼 𝑇W 𝑡 − 𝑇%XX + 𝑇%XX − 𝑇W0 Π(𝑡)

𝜏

Page 24: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Aggregationofleakystorageunits

Heterogeneousmodelsaremergedintoaconservativeenvelopeflexibilitymodel.

Themodelissufficientandlinear,foreasyembeddingindispatchmodels.

P(t)

αS(t)

Smax

S0Smin

S(t)

P(t)

αS(t)

Smax

S0Smin

S(t)

P(t)

αS(t)

Smax

S0Smin

S(t)

P(t)

αS(t)

Smax

S0Smin

S(t)

P(t)

αS(t)

Smax

S0Smin

S(t)

Page 25: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Propertiesofapplianceclasses

𝑑𝜎(𝑡)𝑑𝑡 = −

𝜎 𝑡�̂� + Π(𝑡)

Πc>MK ≤ Π 𝑡 ≤ Πc>&V𝜎d>MK ≤ 𝜎(𝑡) ≤ 𝜎d>&V

_ 𝜎 𝑡 𝑑𝑡 = 1`

0

(heterogeneous)

Page 26: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Suitabilityofapplianceforservicedelivery

a b c

a b c

Powermodulation≅ Πc>&V- Πc>MK

30-minserviceprovision Energyarbitrage≅ 𝜏 𝜎d>&V − 𝜎d>MK

Firstorderanalysisofsuitabilityfordifferentservices

Page 27: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

AggregatorAggregateloaddispatchModelofcollective

flexibility

controlsigna

l

resourceavailability

flexibilityproducts

Markets

P(t)

αS(t)

Smax

S0Smin

S(t)

Page 28: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Casestudy:Optimaluseofdifferentdeviceclasses

domestic commercial

OptimalserviceallocationinapproximateGBsystem• Energyarbitrage• Shortterm(30sec)frequency

services• Mediumterm(30min)

frequencyservices

Serviceallocationsreflectphysicalcharacteristics:• Slowthermaltimeconstants

aregoodforenergyarbitrage• Lowdutycyclesindomestic

appliancesleaveheadroomforhighfrequencyresponse.

Page 29: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Simulation:optimalallocationofflexibility59,524refrigerators(1MW);24-hourallocation

simulationofheterogeneous

refrigerators

Usingrefrigeratorstoprovideenergyarbitrageandfrequencyservices,makingoptimaluseofdeviceflexibility

deliveryoffrequencyresponse

Page 30: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Simulation:Clusteringofsimilarappliancesforbetterserviceallocation

Leakystoragemodelforflexibilityislimitedbylowestcommondenominator.

Flexibilitycanbeimprovedbycreatingclustersofappliancesaccordingtotheircapabilities.

Page 31: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

WRAPPINGUP

Page 32: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Thedemandresponsecontrolspectrum

Ourapproach:distribution-referredcontrol• Devicesaresemi-autonomous• Collectivegoals aresetcentrally• Actions aredecidedlocally,withreference

toexpectedgroupbehaviour

Indirectcontrolusingincentives

Decentralised actionsonthebasisofanon-localcontrolsignals.

UsefultaxonomyofindirectcontrolinHeussen etal.,IEEEPESISGTEurope2012

Goalsandactionsaredecidedcentrally,orinadistributedfashion

Directdispatchofflexibleresources

• Requiresreal-timecommunication• Limitedautonomy• Privacyconcerns

V • Controllability

Page 33: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

CommunicationrequirementsRobust‘semi-autonomous’operation

aggregator

responsedesigncycle

operationalcycle

1. Measuretemperature2. Updatemodel3. Switchon/off

realtime

aheadoftime

Significantchangesin:• thermalmodel• constraints

real-timecontrolsignalΠ(𝑡)

or

powerresponsemodel

Page 34: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Summaryandoutlook

Wehavedevelopedanend-to-endcontrolschemeforTCLsthatis• Nondisruptive: fridgesrespectlocalconstraintsatalltimesandarefree

torespondtoindividualcoolingrequirements• Decentralised:(semi-)autonomousoperationdoesnotrequirereal-time

commandandcontrolinfrastructure• Accurate:accuratecontroloveraggregatepowerconsumption,despite

on/offcharacteroffridgesandpopulationheterogeneity

Openquestionsandfurtherdevelopment• Robustness: Howsensitiveisthisschemetomodelmisspecification?• Thelimitoflargebutnotinfinitenumbers:meanfieldfeedbackeffects• Complexconstraints:lockout,temperaturevariations• Labtesting:firstsmallscale,thenlarge.• Optimaldispatchofleakystorageunits.

Page 35: Smart refrigerators: a distribution-referred … refrigerators: a distribution-referred ... (t) P(t) αS(t) Smax S0 Smin ... distributed fashion Direct dispatch of flexible

Wanttoknowmore?

[email protected]

• Decentralisedcontrolofthermostaticloadsforflexibledemandresponse.SimonTindemans,VincenzoTrovato,GoranStrbacIEEETransactionsonControlSystemsTechnology (2015)

• TheLeakyStorageModelforoptimalmulti-serviceallocationofthermostaticloads.VincenzoTrovato,SimonTindemans,GoranStrbacIETGeneration,Transmission&Distribution(2016)

• Nondisruptive decentralizedcontrolofthermalloadswithsecondorderthermalmodelsSimonTindemans,GoranStrbac2016IEEEPESGeneralMeeting,Boston.

• UnderstandingtheAggregateFlexibilityofThermostaticallyControlledLoadsVincenzoTrovato,SimonTindemans,GoranStrbacPowerTech 2017(submitted)