Skip to content

Coop battery implementation #137

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
GillisHommen merged 6 commits into from
Mar 31, 2025
Merged

Coop battery implementation #137

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
6 commits
Select commit Hold shift + click to select a range
f95d9e6
Eerste setup coop batterij
AteZenmo Mar 21, 2025
4be5401
grid battery is put on false again if no coop
AteZenmo Mar 21, 2025
094fdb4
activation of coop battery should be coming from a trigger in the int…
AteZenmo Mar 26, 2025
0cbad9c
Merge branch 'main' into coopBattery
AteZenmo Mar 27, 2025
609dc89
Merge branch 'main' into coopBattery
AteZenmo Mar 27, 2025
323962b
Merge branch 'main' into coopBattery
AteZenmo Mar 28, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 4 additions & 0 deletions Zero_engine.alpx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -973,6 +973,10 @@
<Id>1725879557284</Id>
<Name><![CDATA[SELF_CONSUMPTION]]></Name>
</Option>
<Option>
<Id>1742557314771</Id>
<Name><![CDATA[BALANCE_COOP]]></Name>
</Option>
<Option>
<Id>1734001947459</Id>
<Name><![CDATA[OFF]]></Name>
Expand Down
16 changes: 16 additions & 0 deletions _alp/Agents/EnergyCoop/Code/Functions.java
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1186,6 +1186,7 @@
f_calculateKPIs();
}

f_connectCoopBattery();
/*ALCODEEND*/}

double f_getGroupContractDeliveryCapacity_kW()
Expand Down Expand Up @@ -1664,3 +1665,18 @@ HashSet<GridConnection> f_getAllChildCustomerGridConnections_recursion(HashSet<G

/*ALCODEEND*/}

double f_connectCoopBattery()
{/*ALCODESTART::1742569887460*/
GCGridBattery coopBattery = findFirst(energyModel.GridBatteries, bat -> bat.p_batteryOperationMode == OL_BatteryOperationMode.BALANCE_COOP);

if(coopBattery != null){
//Reset previous state
coopBattery.v_previousPowerElectricity_kW = 0;

//Connect to coop
coopBattery.c_parentCoops.add(this);
c_memberGridConnections.add(coopBattery);
v_liveAssetsMetaData.hasBattery = true;
}
/*ALCODEEND*/}

16 changes: 16 additions & 0 deletions _alp/Agents/EnergyCoop/Code/Functions.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -565,4 +565,20 @@
<ShowLabel>true</ShowLabel>
<Body xmlns:al="http://anylogic.com"/>
</Function>
<Function AccessType="default" StaticFunction="false">
<ReturnModificator>VOID</ReturnModificator>
<ReturnType>double</ReturnType>
<Id>1742569887460</Id>
<Name><![CDATA[f_connectCoopBattery]]></Name>
<X>890</X>
<Y>160</Y>
<Label>
<X>10</X>
<Y>0</Y>
</Label>
<PublicFlag>false</PublicFlag>
<PresentationFlag>true</PresentationFlag>
<ShowLabel>true</ShowLabel>
<Body xmlns:al="http://anylogic.com"/>
</Function>
</Functions>
3 changes: 3 additions & 0 deletions _alp/Agents/EnergyModel/Code/Functions.java
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1773,6 +1773,9 @@ EnergyCoop f_removeEnergyCoop(EnergyCoop energyCoop)
// Removing this coop from the list of coops in the GC
for (GridConnection GC : energyCoop.f_getAllChildMemberGridConnections()) {
GC.c_parentCoops.remove(energyCoop);
if(GC instanceof GCGridBattery && GC.p_batteryOperationMode == OL_BatteryOperationMode.BALANCE_COOP){
GC.f_setActive(false);
}
}

// Remove energyCoop from pop_energyCoops.
Expand Down
55 changes: 55 additions & 0 deletions _alp/Agents/GCGridBattery/Code/Functions.java
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -68,8 +68,11 @@ else if(EA.j_ea instanceof J_EAStorageElectric && ((J_EAStorageElectric)EA.j_ea)
f_batteryManagementPrice(v_batterySOC_fr);
} else if (p_batteryOperationMode == OL_BatteryOperationMode.NODAL_PRICING){
f_batteryManagementPriceGrid(v_batterySOC_fr);
} else if (p_batteryOperationMode == OL_BatteryOperationMode.BALANCE_COOP){
f_batteryManagementBalanceCOOP(v_batterySOC_fr);
}


p_batteryAsset.f_updateAllFlows(p_batteryAsset.v_powerFraction_fr);
//J_FlowsMap flowsMap = flowsPair.getFirst();
//v_batteryPowerElectric_kW = flowsMap.get(OL_EnergyCarriers.ELECTRICITY);
Expand Down Expand Up @@ -342,3 +345,55 @@ else if (WTP_discharge_eurpkWh < currentElectricityPriceCharge_eurpkWh) {

/*ALCODEEND*/}

double f_batteryManagementBalanceCOOP(double batterySOC)
{/*ALCODESTART::1742557292690*/

if(c_parentCoops.isEmpty()){
return;
}

EnergyCoop parentCoop = c_parentCoops.get(0);

if (p_batteryAsset.getStorageCapacity_kWh() != 0){
double CoopConnectionCapacityDelivery_kW = parentCoop.v_liveConnectionMetaData.contractedDeliveryCapacity_kW * 0.9; //10% reduction, to accord for time step delay
double CoopConnectionCapacityFeedin_kW = parentCoop.v_liveConnectionMetaData.contractedDeliveryCapacity_kW * 0.9; //10% reduction, to accord for time step delay
double v_previousPowerBattery_kW = v_previousPowerElectricity_kW;// Assumes battery is only asset on gridconnection!! p_batteryAsset.electricityConsumption_kW-p_batteryAsset.electricityProduction_kW;
double currentCoopElectricityBalance_kW = parentCoop.v_liveData.data_liveElectricityBalance_kW.getY(parentCoop.v_liveData.data_liveElectricityBalance_kW.size() - 1) - v_previousPowerBattery_kW;

double availableChargePower_kW = CoopConnectionCapacityDelivery_kW - currentCoopElectricityBalance_kW; // Max battery charging power within grid capacity
double availableDischargePower_kW = CoopConnectionCapacityFeedin_kW + currentCoopElectricityBalance_kW; // Max discharging power within grid capacity
//Tuning parameters
double FeedbackGain_kWpSOC = 3 * p_batteryAsset.getCapacityElectric_kW(); // How strongly to aim for SOC setpoint
double FeedforwardGain_kWpKw = 0.1; // Feedforward based on current surpluss in Coop
double chargeOffset_kW = 0; // Charging 'bias', basically increases SOC setpoint slightly during the whole day.
double chargeSetpoint_kW = 0;

//Congestion prevention
if (availableChargePower_kW < 0) { //Delivery side
p_batteryAsset.v_powerFraction_fr = max(-1, availableChargePower_kW / p_batteryAsset.getCapacityElectric_kW());
return;
}
if (availableDischargePower_kW < 0) { //Feedin side
p_batteryAsset.v_powerFraction_fr = min(1, -availableDischargePower_kW / p_batteryAsset.getCapacityElectric_kW());
return;
}

if (energyModel.v_currentSolarPowerNormalized_r > 0.1) {
if (currentCoopElectricityBalance_kW < 0) {
p_batteryAsset.v_powerFraction_fr = max(-1, min(1, -currentCoopElectricityBalance_kW / p_batteryAsset.getCapacityElectric_kW()));
}
}
else {
double expectedWind_kWh = parentCoop.v_liveAssetsMetaData.totalInstalledWindPower_kW * energyModel.v_WindYieldForecast_fr * energyModel.p_forecastTime_h;
double expectedSolar_kWh = parentCoop.v_liveAssetsMetaData.totalInstalledPVPower_kW * energyModel.v_SolarYieldForecast_fr * energyModel.p_forecastTime_h;
double incomingPower_fr = (expectedSolar_kWh + expectedWind_kWh) / p_batteryAsset.getStorageCapacity_kWh();
double SOC_setp_fr = 1 - incomingPower_fr;

chargeSetpoint_kW = FeedbackGain_kWpSOC*(SOC_setp_fr - p_batteryAsset.getCurrentStateOfCharge());
chargeSetpoint_kW = min(max(chargeSetpoint_kW, availableDischargePower_kW),availableChargePower_kW); // Don't allow too much (dis)charging!
p_batteryAsset.v_powerFraction_fr = max(-1,min(1, chargeSetpoint_kW / p_batteryAsset.getCapacityElectric_kW()));
}
}

/*ALCODEEND*/}

24 changes: 24 additions & 0 deletions _alp/Agents/GCGridBattery/Code/Functions.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -164,4 +164,28 @@
</Parameter>
<Body xmlns:al="http://anylogic.com"/>
</Function>
<Function AccessType="public" StaticFunction="false">
<ReturnModificator>VOID</ReturnModificator>
<ReturnType>double</ReturnType>
<Id>1742557292690</Id>
<Name><![CDATA[f_batteryManagementBalanceCOOP]]></Name>
<Description><![CDATA[Regeling voor balancering van Coop-load.

@author: Gillis Hommen
@since: 28-12-2022]]></Description>
<X>1160</X>
<Y>750</Y>
<Label>
<X>10</X>
<Y>0</Y>
</Label>
<PublicFlag>false</PublicFlag>
<PresentationFlag>true</PresentationFlag>
<ShowLabel>true</ShowLabel>
<Parameter>
<Name><![CDATA[batterySOC]]></Name>
<Type><![CDATA[double]]></Type>
</Parameter>
<Body xmlns:al="http://anylogic.com"/>
</Function>
</Functions>
12 changes: 12 additions & 0 deletions _alp/Agents/GridConnection/Code/Functions.java
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2471,6 +2471,12 @@ else if (flowsMap.get(EC) > 0){
energyModel.v_liveAssetsMetaData.totalInstalledWindPower_kW -= v_liveAssetsMetaData.totalInstalledWindPower_kW;
energyModel.v_liveAssetsMetaData.totalInstalledBatteryStorageCapacity_MWh -= v_liveAssetsMetaData.totalInstalledBatteryStorageCapacity_MWh;

for(EnergyCoop coop : c_parentCoops){
coop.v_liveAssetsMetaData.totalInstalledPVPower_kW -= v_liveAssetsMetaData.totalInstalledPVPower_kW;
coop.v_liveAssetsMetaData.totalInstalledWindPower_kW -= v_liveAssetsMetaData.totalInstalledWindPower_kW;
coop.v_liveAssetsMetaData.totalInstalledBatteryStorageCapacity_MWh -= v_liveAssetsMetaData.totalInstalledBatteryStorageCapacity_MWh;
}

// Reset Connection Capacity to default
f_nfatoSetConnectionCapacity(true);

Expand Down Expand Up @@ -2523,6 +2529,12 @@ else if (flowsMap.get(EC) > 0){
energyModel.v_liveAssetsMetaData.totalInstalledPVPower_kW += v_liveAssetsMetaData.totalInstalledPVPower_kW;
energyModel.v_liveAssetsMetaData.totalInstalledWindPower_kW += v_liveAssetsMetaData.totalInstalledWindPower_kW;
energyModel.v_liveAssetsMetaData.totalInstalledBatteryStorageCapacity_MWh += v_liveAssetsMetaData.totalInstalledBatteryStorageCapacity_MWh;

for(EnergyCoop coop : c_parentCoops){
coop.v_liveAssetsMetaData.totalInstalledPVPower_kW += v_liveAssetsMetaData.totalInstalledPVPower_kW;
coop.v_liveAssetsMetaData.totalInstalledWindPower_kW += v_liveAssetsMetaData.totalInstalledWindPower_kW;
coop.v_liveAssetsMetaData.totalInstalledBatteryStorageCapacity_MWh += v_liveAssetsMetaData.totalInstalledBatteryStorageCapacity_MWh;
}
}

//Update the 'isActive' variable
Expand Down