# Energy Storage – PVS

**Energy Storage System**

For now, the scope of this feature is limited to AC-coupled storage systems. In practice, this means that the storage is connected to the PV system at the MV level, after the MV transformers and associated cabling. The storage system (shaded in blue below) consists of: batteries, inverters and MV transformers and cabling.

*Figure 1 – Single-line diagram of an AC-coupled storage system*

*Nodal positions 1-7 are defined for reference. Power is reported for certain nodes in the downloadable nodal data.*

**Inputs**

sdc_{np} |
Nameplate Energy Capacity |
[Wh] |

s_{ussoc} |
Energy Capacity Factor |
unitless |

sdc_{us} |
Initial Usable Energy Capacity |
[Wh] |

s_{caldeg} |
Energy Capacity Calendar Degradation |
[%/year] |

s_{cycdeg} |
Energy Capacity Cycling Degradation |
[%/cycle] |

dc_{rte} |
Initial Roundtrip DC Efficiency |
[%] |

rte_{caldeg} |
Roundtrip DC Efficiency Calendar Degradation |
[%/year] |

rte_{cycdeg} |
Roundtrip DC Efficiency Cycling Degradation |
[%/cycle] |

sac |
Inverter Real Power |
[W] |

eff_{sinv} |
Inverter Efficiency |
[%] |

P_{Loss,MV} |
Storage MV Transformer Losses |
[W] |

P_{Loss,HV} |
Plant HV Transformer Losses |
[W] |

P_{Loss,T} |
Transmission Line Losses |
[W] |

s_{auxconst} |
HVAC Constant Loss |
[W/MWh] |

s_{auxop} |
HVAC Operating Loss |
[W/MW] |

ic |
Power Output Limit |
[W] |

P_{avail} |
Plant Power with Availability Loss Applied |
[W] |

P_{AC,out,PV} |
PV MV Transformer Output |
[W] |

TP |
Target Period Hour |
{User True/False selection} |

DM |
Dispatch Model |
{User selection} |

Δt |
Time step |
[s] |

**Outputs**

P_{AC,out,S} |
Net Output (high-side Storage MV Transformers) AC power |
[W] |

**Algorithm**

__Built-In Dispatch Algorithms__

The dispatch logic considers the power produced by the PV plant, the interconnect capacity, and the State of Charge of the ESS, in addition to the value of energy output to the grid in each hour throughout the day. The logic takes into account energy losses and equipment power limitations.

- First, the ESS State of Charge at the beginning of the day is determined (initialized at full usable State of Charge). The PV power output profile for the day is collected. A target period of high-priority hours during the day is identified.

a. The charging algorithms offered are listed below:

i. **Interconnect Excess Algorithm**: The storage system will be charged if and only if the loss included PV output capability exceeds the Power Output Limit. It will discharge energy at its maximum capability during the target period window.

ii. **Energy Available Algorithm**: The storage system will be charged as and when PV energy becomes available with the sole objective of charging it to full capacity. Once the storage system is fully charged, the energy from the PV system will then be re-directed to the grid. The storage system will discharge energy at maximum capability during the target period window.

b. The maximum power by which the storage is charged is limited by the total power rating of the storage inverters, and the PV power available.

c. The maximum energy which can be stored is limited by the Usable Energy Capacity.

d. Both Usable Energy Capacity (energy the storage can hold) and Roundtrip Efficiency degrade as a function of time and # of cycles (where # of cycles is cumulative energy charged/usable energy capacity).

e. Roundtrip Efficiency is only applied when charging the storage (not when discharging), so the State of Charge shows the available dischargable DC energy.

f. Inverters associated with the storage system are modeled with flat efficiencies.

g. Transformers associated with the storage system are modeled as for the PV system, with an auxiliary consumption defined by a no-load loss and a full-load loss.

2. Next, the power that is required to meet the interconnect capacity not met by PV output during target hours is calculated, where the maximum power to be discharged is limited by the ESS power rating.

3. The ESS then discharges the amount of power required to meet the interconnect limit during each target hour, limited by the ESS power rating, until the ESS depletes the stored energy.

__Custom Dispatch__

The custom dispatch option allows the user more control of the storage dispatch.

The user defines a charge and discharge target profile by identifying a “Charge” or “Discharge” command and power targets as a fraction of the total inverter rated capacity, for an index corresponding to each time step in the prediction weather file. Indexes without a command will not result in any charge or discharge.

The target power is defined at the input to the storage MV transformers, and will be limited by the total power rating of the storage inverters, the State of Charge of the storage, the total PV generated, or the total power rating of the MV transformers, as relevant.

The portion of commanded power that exceeds the limitation is not diverted as commanded.

However, a command to discharge a level of power that results in the total discharged and PV generated power exceeding the interconnection capacity is followed but results in the excess power being clipped.

**Algorithm**

a) Determine the amount of power exceeding the interconnect from just the PV system. For the Interconnect Excess algorithm, this will determine how much power to divert from the PV system to charge the storage system.

PVExceedingInterconnect_{6 }= max(**P _{avail}** –

**ic**, 0)

b) Determine the excess power available at Node 4 for charging the storage system after losses.

AvailableExcesstoCharge_{4 }= min(**sac**,(PVExceedingInterconnect_{6 + }**P _{Loss,HV }+ P_{Loss,T} – P_{Loss,MV}**) *

**eff**)

_{sinv}c) Determine the power at Node 4 that the storage system has capacity to charge in addition to the excess power. Evaluate discharge limits. Calculate the maximum amount of power that can be discharged from the storage system after considering the PV output and Power Output Limit

AvailableCapacitytoChargefromNonExcess_{4} = min(**sac, P _{AC,out,PV }– P_{Loss,MV}) * eff_{sinv}**

d) Evaluate discharge limits. Calculate the maximum amount of power that can be discharged from the storage system after considering the PV output and interconnect limit.

AvailableCapacitytoDischarge_{6,LGIALimit} = max(**ic **–** P _{avail}**, 0)

e) Calculate the maximum amount of power that can be discharged from the storage system, translated to the interconnect (Node 6), after considering the storage inverter limits

AvailableCapacitytoDischarge_{6,StorageLimitt} = **sac – P _{Loss,MV}**(full load)

**– P**–

_{Loss,HV }**P**

_{Loss,T}The maximum power possible to discharge through the plant is the lesser of the two above constraints.

AvailableCapacitytoDischarge_{6} = min(AvailableCapacitytoDischarge_{6,LGIALimit}, AvailableCapacitytoDischarge_{6,StorageLimitt})

Calculated at Node 2:

AvailableCapacitytoDischarge_{2} = AvailableCapacitytoDischarge_{6} ** + P _{Loss,HV }** +

**P**

_{Loss,T}Calculated at Node 3, with MV transformer losses applicable to the power:

AvailableCapacitytoDischarge_{3} = AvailableCapacitytoDischarge_{2} ** + P _{Loss,MV}( **AvailableCapacitytoDischarge

_{2})

f) Determine the maximum power that the system has capacity to discharge.

AvailableCapacitytoDischarge_{4} = min(AvailableCapacitytoDischarge_{3}, **sac**)/** eff _{sinv}**

g) Determine the charge/discharge schedule.

Evaluate whether the time step is a Target Period for discharge.

If **TP**=True, Discharge=1

Identify hours in which power can charge based on the selected model type.

If **DM** = LGIAExcess and **P _{avail}**>

**ic**,

Charge=1 and AvailableCapacitytoCharge = AvailableExcesstoCharge_{4}

If **DM** = EnergyAvailable and AvailableCapacitytoChargefromNonExcess_{4} >0 and if **TP**=False,

Charge=1 and AvailableCapacitytoCharge = AvailableCapacitytoChargefromNonExcess_{4}

h) Calculate the maximum usable energy capacity after degradation. For each time step, reference from the previous time step: the previous System Age, assuming the it is 0 years in the first time step; the previous Battery AC Discharge at Node 3, assuming it is 0MW in the first time step; the previous Battery Maximum Capacity, assuming the Battery Maximum Capacity in the first time step is set at **sac**; and the previous Usable Energy Capacity Cycling-Induced Degradation, assuming it is 0% in the first time step.

UsableEnergyCap_CycleDeg = (prevBatteryACDischarge_{3}/** eff _{sinv}**)/prevBatteryMaxCapacity *

**s**+ prevUsableEnergyCap_CycleDeg

_{cycdeg}SystemAge = Δt+prevSystemAge

UsableEnergyCap_CalendarDeg = **s _{caldeg}** * SystemAge

BatteryMaxCapacity = **sdc _{us}*** (1- (UsableEnergyCap_CycleDeg + UsableEnergyCap_CalendarDeg))

i) Calculate the losses due to DC Roundtrip Efficiency. For each time step, reference: the previous Battery DC Roundtrip Efficiency, assuming in the first time step it is set at **dc _{rte}**; and the previous Battery Roundtrip Efficiency Cycling Degradation, assuming it is 0% in the first time step.

BatteryRTE_CycleDeg = (prevBatteryACDischarge_{3}/** eff _{sinv}**)/prevBatteryDCRTE *

**rte**+ prevBatteryRTE_CycleDeg

_{cycdeg}BatteryRTE_CalendarDeg = **rte _{caldeg}** * SystemAge

BatteryDCRTE = **dc _{rte}** * (1- (BatteryRTE_CycleDeg + BatteryRTE_CalendarDeg))

j) Find the State of Charge of the storage system at each time step. For each time step, reference the previous Battery SOC, assuming the Battery SOC in the first time step is set at **sdc _{us}**.

If Charge=1,

BatterySOC = min(BatteryMaxCapacity, prevBatterySOC+ AvailableCapacitytoCharge * BatteryDCRTE)

If Discharge=1,

BatterySOC =max(0, prevBatterySOC+ AvailableCapacitytoDischarge_{4})

Else,

BatterySOC = prevBatterySOC

** **

k) Calculate the power throughput on the high side of the storage system transformers.

If BatterySOC – prevBatterySOC > 0,

BatteryDCPower_{4} = (BatterySOC – prevBatterySOC)/ BatteryDCRTE

BatteryACCharge_{3} = BatteryDCPower_{4}/**eff _{sinv}**

**P _{AC,out,S} = **BatteryACCharge

_{3}+

**P**

_{Loss,MV}Else,

BatteryDCPower_{4} = (BatterySOC – prevBatterySOC)

BatteryACDischarge_{3} = BatteryDCPower_{4*}**eff _{sinv}**

**P _{AC,out,S} = **BatteryACDischarge

_{3}–

**P**

_{Loss,MV}