Skip to content

Districtheating_merit_order_and_losses #301

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

Open
AteZenmo wants to merge 2 commits into
base: main
Choose a base branch
from
Open

Districtheating_merit_order_and_losses #301

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
d7ae3fd
Districtheating_merit_order_and_losses
AteZenmo Jun 17, 2026
82f7804
bugfix: should ofcourse only call energyBalance for active GC distric…
AteZenmo Jun 23, 2026
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 @@ -2229,6 +2229,10 @@ LARGE_CONNECTION (grootverbruik, > 3x80)
<Name><![CDATA[J_BackupGeneratorManagementExternalSetpoint]]></Name>
<Folder>1776078970513</Folder>
</JavaClass>
<JavaClass>
<Id>1781710977201</Id>
<Name><![CDATA[GISUtil]]></Name>
</JavaClass>
</JavaClasses>
<RequiredLibraryReference>
<LibraryName>com.anylogic.libraries.modules.markup_descriptors</LibraryName>
Expand Down
14 changes: 11 additions & 3 deletions _alp/Agents/EnergyModel/Code/Functions.java
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -159,14 +159,22 @@
v_batteryStoredEnergy_kWh = 0;

if (b_parallelizeGridConnections) {
c_gridConnections.parallelStream().forEach(gc -> gc.f_calculateEnergyBalance(p_timeVariables, v_isRapidRun));
c_gridConnections.parallelStream().forEach(gc -> {if(!(gc instanceof GCDistrictHeating)){ gc.f_calculateEnergyBalance(p_timeVariables, v_isRapidRun);}});
}
else {
for(GridConnection gc : c_gridConnections) {
gc.f_calculateEnergyBalance(p_timeVariables, v_isRapidRun);
if(!(gc instanceof GCDistrictHeating)){
gc.f_calculateEnergyBalance(p_timeVariables, v_isRapidRun);
}
}
}

for(GCDistrictHeating GCDH : DistrictHeatingSystems){
if(GCDH.f_isActive()){
GCDH.f_calculateEnergyBalance(p_timeVariables, v_isRapidRun);
}
}

//OLD LOCATION OF AGGREGATOR CALL, NECESSARY IF MANAGEMENT IS WITHOUT TIME STEP DELAY

for(GridConnection gc : c_gridConnections) { // Can't do this in parallel due to different threads writing to the same values!
Expand Down Expand Up @@ -502,7 +510,7 @@
double f_gridNodeRecursiveAdd(GridNode GN)
{/*ALCODESTART::1716886716306*/
c_gridNodeExecutionList.add(GN);
for (GridNode GNchild : GN.c_connectedGridNodes) {
for (GridNode GNchild : GN.f_getConnectedGridNodes()) {
f_gridNodeRecursiveAdd(GNchild);
}

Expand Down
6 changes: 6 additions & 0 deletions _alp/Agents/EnergyModel/EmbeddedObjects.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -90,6 +90,12 @@
<Parameter>
<Name><![CDATA[p_profileType]]></Name>
</Parameter>
<Parameter>
<Name><![CDATA[p_heatLossFactorDistribution_kW]]></Name>
</Parameter>
<Parameter>
<Name><![CDATA[p_heatLossFactorTransport_kWpm]]></Name>
</Parameter>
</Parameters>
<ReplicationFlag>true</ReplicationFlag>
<Replication Class="CodeValue">
Expand Down
12 changes: 6 additions & 6 deletions _alp/Agents/GridNode/AOC.GridNode.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -115,7 +115,7 @@ if (p_energyType.equals(OL_EnergyCarriers.HEAT)) {
<Id>1705919330576</Id>
<Name><![CDATA[data_liveLoad_kW]]></Name>
<X>50</X>
<Y>850</Y>
<Y>880</Y>
<Label>
<X>15</X>
<Y>0</Y>
Expand Down Expand Up @@ -144,7 +144,7 @@ if (p_energyType.equals(OL_EnergyCarriers.HEAT)) {
<Id>1712668669449</Id>
<Name><![CDATA[data_liveCapacityDemand_kW]]></Name>
<X>50</X>
<Y>890</Y>
<Y>920</Y>
<Label>
<X>15</X>
<Y>0</Y>
Expand Down Expand Up @@ -173,7 +173,7 @@ if (p_energyType.equals(OL_EnergyCarriers.HEAT)) {
<Id>1712673776914</Id>
<Name><![CDATA[data_liveCapacitySupply_kW]]></Name>
<X>50</X>
<Y>870</Y>
<Y>900</Y>
<Label>
<X>15</X>
<Y>0</Y>
Expand Down Expand Up @@ -202,7 +202,7 @@ if (p_energyType.equals(OL_EnergyCarriers.HEAT)) {
<Id>1715167370828</Id>
<Name><![CDATA[data_totalLoad_kW]]></Name>
<X>50</X>
<Y>920</Y>
<Y>950</Y>
<Label>
<X>15</X>
<Y>0</Y>
Expand Down Expand Up @@ -231,7 +231,7 @@ if (p_energyType.equals(OL_EnergyCarriers.HEAT)) {
<Id>1715167392629</Id>
<Name><![CDATA[data_summerWeekLoad_kW]]></Name>
<X>50</X>
<Y>940</Y>
<Y>970</Y>
<Label>
<X>15</X>
<Y>0</Y>
Expand Down Expand Up @@ -260,7 +260,7 @@ if (p_energyType.equals(OL_EnergyCarriers.HEAT)) {
<Id>1715167408643</Id>
<Name><![CDATA[data_winterWeekLoad_kW]]></Name>
<X>50</X>
<Y>960</Y>
<Y>990</Y>
<Label>
<X>15</X>
<Y>0</Y>
Expand Down
159 changes: 61 additions & 98 deletions _alp/Agents/GridNode/Code/Functions.java
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -54,28 +54,6 @@
}
}*/

/*ALCODEEND*/}

double f_sumLoads()
{/*ALCODESTART::1660122738707*/
v_currentLoad_kW = 0;

// determine the net energy flows from all subconnections by nodetype

for( GridNode GN : c_connectedGridNodes ) {
v_currentLoad_kW += GN.v_currentLoad_kW;
}

for( GridConnection GC : c_connectedGridConnections) {
v_currentLoad_kW += GC.fm_currentBalanceFlows_kW.get(p_energyCarrier);
}


/*if( p_energyType == OL_EnergyCarriers.ELECTRICITY ){
v_electricLoadRatioExclBattery = v_currentLoadElectricity_kW / p_capacity_kW;
}*/


/*ALCODEEND*/}

double f_nodeMetering(J_TimeVariables timeVariables,J_TimeParameters timeParameters,boolean isRapidRun)
Expand Down Expand Up @@ -179,62 +157,6 @@ else if (p_energyCarrier == OL_EnergyCarriers.HEAT) {
}
/*ALCODEEND*/}

double f_instantiateEnergyAssets()
{/*ALCODESTART::1676387466304*/
//traceln("asset " + p_energyAssetList);

if( p_energyAssetList != null) {
for( JsonNode l : p_energyAssetList ) {
OL_EACategories assetCategory = OL_EACategories.valueOf(l.required( "category" ).textValue());
switch( assetCategory ) {
case STORAGE:
//EnergyAsset storageAsset = main.add_pop_energyAssets();
//storageAsset.set_p_parentAgentID( this.p_gridNodeID );
OL_EnergyAssetType assetType = ( OL_EnergyAssetType.valueOf(l.required( "type" ).textValue()));
//storageAsset.set_p_defaultEnergyAssetPresetName(l.required( "name" ).textValue());
if( assetType == OL_EnergyAssetType.STORAGE_HEAT){
double capacityHeat_kW = l.path( "capacityHeat_kW").doubleValue();
double lossFactor_WpK = l.path( "lossFactor_WpK" ).doubleValue();
double heatCapacity_JpK = l.path( "heatCapacity_JpK" ).doubleValue();
double minTemperature_degC = l.path( "minTemp_degC" ).asDouble(40.0); // provide default values
double maxTemperature_degC = l.path( "maxTemp_degC" ).asDouble(90.0);
double setTemperature_degC = l.path( "setTemp_degC" ).asDouble(60.0);
double initialTemperature_degC = l.path( "initialTemperature_degC" ).doubleValue();
String ambientTempType2 = l.path( "ambientTempType" ).textValue();

// minTemperature_degC = 35; // TO DELETE TEMP FIX
// setTemperature_degC = uniform_discr(38, 45); // TO DELETE TEMP FIX
// initialTemperature_degC = uniform_discr(70,80) * 1.0; // TO DELETE TEMP FIX
// heatCapacity_JpK = heatCapacity_JpK * uniform(0.7, 1.3); // TO DELETE TEMP FIX
traceln("heatstorage asset initialisation check! minTemp = "+minTemperature_degC+", maxTemperature_degC = "+maxTemperature_degC+", setTemp_degC = "+ setTemperature_degC+", initialTemperature_degC = "+initialTemperature_degC);

//traceln("Heat Storage init: minTemperature = "+minTemperature_degC+", maxTemperature_degC = "+maxTemperature_degC+", setTemperature_degC = "+setTemperature_degC+", initialTemperature_degC = "+initialTemperature_degC);


p_transportBuffer = new J_EAStorageHeat((Agent)this, OL_EAStorageTypes.HEATBUFFER, capacityHeat_kW, lossFactor_WpK, energyModel.p_timeStep_h, initialTemperature_degC, minTemperature_degC, maxTemperature_degC, setTemperature_degC, heatCapacity_JpK, ambientTempType2);
//J_EAStorageHeat(Agent parentAgent, OL_EAStorageTypes heatStorageType, double capacityHeat_kW, double lossFactor_WpK, double timestep_h, double initialTemperature_degC, double minTemperature_degC, double maxTemperature_degC, double setTemperature_degC, double heatCapacity_JpK, String ambientTempType ) {
//p_transportBuffer = storageAsset.j_ea;
p_transportBuffer.updateAmbientTemperature( energyModel.p_undergroundTemperature_degC );
if(heatCapacity_JpK > 0 & capacityHeat_kW > 0) {
b_transportBufferValid = true;
}
}
else{
traceln("F_instantiateEnergyAssets: ERROR, storage asset type not available");
}
//storageAsset.f_connectToParentNode( this );
//main.c_storageAssets.add(storageAsset);

break;
default:
traceln("not a valid energy asset category." + assetCategory);
break;
}
}
//traceln("GridConnection "+this.p_gridNodeID+" has finished initializing its energyAssets!");
}
/*ALCODEEND*/}

double f_calculateEnergyBalance(J_TimeVariables timeVariables,J_TimeParameters timeParameters,boolean isRapidRun)
{/*ALCODESTART::1688370981599*/
f_sumLoads();
Expand Down Expand Up @@ -302,16 +224,6 @@ else if (p_energyCarrier == OL_EnergyCarriers.HEAT) {
f_getCurrentChargingInformation(timeParameters);
/*ALCODEEND*/}

double f_addGridBatteryLoad()
{/*ALCODESTART::1688372319365*/
for( Agent a : subConnections.getConnections() ) {
if ( a instanceof GCGridBattery){
v_currentLoadElectricity_kW += ((GCGridBattery)a).v_currentPowerElectricity_kW;
}
}

/*ALCODEEND*/}

double f_resetStates(J_TimeParameters timeParameters)
{/*ALCODESTART::1698919552330*/
// Current status
Expand Down Expand Up @@ -530,38 +442,54 @@ J_LoadDurationCurves f_getDuurkrommes(J_TimeParameters timeParameters)

DataSet f_getPeakImportWeekDataSet(J_TimeParameters timeParameters)
{/*ALCODESTART::1754461654982*/
double[] elecBalance_kW = acc_annualElectricityBalance_kW.getTimeSeries_kW();
ZeroAccumulator acc_annualBalance_kW;
if(p_energyCarrier == OL_EnergyCarriers.ELECTRICITY){
acc_annualBalance_kW = acc_annualElectricityBalance_kW;
}
else{
acc_annualBalance_kW = acc_annualHeatBalance_kW;
}

double[] balance_kW = acc_annualBalance_kW.getTimeSeries_kW();

Integer maxIndex = 0; // index with peak import
for (int i = 1; i < elecBalance_kW.length; i++) {
if (elecBalance_kW[i] > elecBalance_kW[maxIndex]) {
for (int i = 1; i < balance_kW.length; i++) {
if (balance_kW[i] > balance_kW[maxIndex]) {
maxIndex = i;
}
}

double peakTime_h = maxIndex*timeParameters.getTimeStep_h();
double duration_h = acc_annualElectricityBalance_kW.getDuration();
double duration_h = acc_annualBalance_kW.getDuration();
double accStartTime_h = min(duration_h-7*24,max(0,peakTime_h - 3.5*24));
DataSet ds = acc_annualElectricityBalance_kW.getDataSet(timeParameters.getRunStartTime_h(), accStartTime_h, accStartTime_h+24*7);
DataSet ds = acc_annualBalance_kW.getDataSet(timeParameters.getRunStartTime_h(), accStartTime_h, accStartTime_h+24*7);

return ds;
/*ALCODEEND*/}

DataSet f_getPeakExportWeekDataSet(J_TimeParameters timeParameters)
{/*ALCODESTART::1754462048041*/
double[] elecBalance_kW = acc_annualElectricityBalance_kW.getTimeSeries_kW();
ZeroAccumulator acc_annualBalance_kW;
if(p_energyCarrier == OL_EnergyCarriers.ELECTRICITY){
acc_annualBalance_kW = acc_annualElectricityBalance_kW;
}
else{
acc_annualBalance_kW = acc_annualHeatBalance_kW;
}

double[] balance_kW = acc_annualBalance_kW.getTimeSeries_kW();

Integer minIndex = 0; // index with peak import
for (int i = 1; i < elecBalance_kW.length; i++) {
if (elecBalance_kW[i] < elecBalance_kW[minIndex]) {
for (int i = 1; i < balance_kW.length; i++) {
if (balance_kW[i] < balance_kW[minIndex]) {
minIndex = i;
}
}

double peakTime_h = minIndex*timeParameters.getTimeStep_h();
double duration_h = acc_annualElectricityBalance_kW.getDuration();
double duration_h = acc_annualBalance_kW.getDuration();
double accStartTime_h = min(duration_h-7*24,max(0,peakTime_h - 3.5*24));
DataSet ds = acc_annualElectricityBalance_kW.getDataSet(timeParameters.getRunStartTime_h(), accStartTime_h, accStartTime_h+24*7);
DataSet ds = acc_annualBalance_kW.getDataSet(timeParameters.getRunStartTime_h(), accStartTime_h, accStartTime_h+24*7);

return ds;
/*ALCODEEND*/}
Expand Down Expand Up @@ -601,3 +529,38 @@ int f_getCurrentNumberOfChargeRequestsBalancingThisGN()
return v_currentChargingPowerBalancingThisGN_kW;
/*ALCODEEND*/}

double f_calculateHeatNodeTransportLoss_kW(double latConnectedGN,double lonConnectedGN)
{/*ALCODESTART::1781699944841*/
double totalDist_m = GISUtil.calculateManhattanDistance_m(p_latitude, p_longitude, latConnectedGN, lonConnectedGN);
double heatLoss_kW = totalDist_m * p_heatLossFactorTransport_kWpm;

return heatLoss_kW;
/*ALCODEEND*/}

double f_sumLoads()
{/*ALCODESTART::1781700026205*/
v_currentLoad_kW = 0;
v_currentLoss_kW = 0;

// determine the net energy flows from all subconnections by nodetype

for( GridNode connectedGN : c_connectedGridNodes ) {
v_currentLoad_kW += connectedGN.v_currentLoad_kW;
if (p_energyCarrier == OL_EnergyCarriers.HEAT && p_parentNodeID != null) {
double heatLoss_kW = 0;
heatLoss_kW = f_calculateHeatNodeTransportLoss_kW(connectedGN.p_latitude, connectedGN.p_longitude);
v_currentLoad_kW += heatLoss_kW;
v_currentLoss_kW += heatLoss_kW;
}
}

for( GridConnection GC : c_connectedGridConnections) {
v_currentLoad_kW += GC.fm_currentBalanceFlows_kW.get(p_energyCarrier);
if (p_energyCarrier == OL_EnergyCarriers.HEAT && p_parentNodeID != null) {
double heatLoss_kW = p_heatLossFactorDistribution_kW;
v_currentLoad_kW += heatLoss_kW;
v_currentLoss_kW += heatLoss_kW;
}
}
/*ALCODEEND*/}

Loading