Skip to content
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

Make the time step (k) dependent on the representative period (rp) #199

Merged
merged 3 commits into from
Oct 25, 2023
Merged

Make the time step (k) dependent on the representative period (rp) #199

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
494e5cd
Add new columns to the rp file
gnawin Oct 25, 2023
0f78626
Make k dependent of rp
gnawin Oct 25, 2023
f78e4d7
Fix error by explicit define rp for a variable
gnawin Oct 25, 2023
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
2 changes: 1 addition & 1 deletion src/io.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -24,7 +24,7 @@ function create_parameters_and_sets_from_file(input_folder::AbstractString)
assets_conversion = assets_data_df[assets_data_df.type.=="conversion", :].name #conversion assets in the energy system
assets_investment = assets_data_df[assets_data_df.investable.==true, :].name #assets with investment method in the energy system
rep_periods = unique(assets_profiles_df.rep_period_id) #representative periods
time_steps = unique(assets_profiles_df.time_step) #time steps in the RP (e.g., hours)
time_steps = Dict(row.id => 1:row.num_time_steps for row in eachrow(rep_period_df)) #time steps in the RP (e.g., hours), that are dependent on RP

# Parameters for system
rep_weight = Dict((row.id) => row.weight for row in eachrow(rep_period_df)) #representative period weight [h]
Expand Down
32 changes: 16 additions & 16 deletions src/model.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -26,10 +26,10 @@ function create_model(graph, params, sets; verbose = false, write_lp_file = fals
set_attribute(model, "output_flag", verbose)

# Variables
@variable(model, flow[F, RP, K]) #flow from asset a to asset aa [MW]
@variable(model, flow[F, rp ∈ RP, K[rp]]) #flow from asset a to asset aa [MW]
@variable(model, 0 ≤ assets_investment[Ai], Int) #number of installed asset units [N]
@variable(model, 0 ≤ flows_investment[Fi], Int)
@variable(model, 0 ≤ storage_level[As, RP, K])
@variable(model, 0 ≤ storage_level[As, rp ∈ RP, K[rp]])

# TODO: Fix storage_level[As, RP, 0] = 0

Expand All @@ -56,7 +56,7 @@ function create_model(graph, params, sets; verbose = false, write_lp_file = fals
model,
sum(
params.rep_weight[rp] * params.flows_variable_cost[f] * flow[f, rp, k] for
f ∈ F, rp ∈ RP, k ∈ K
f ∈ F, rp ∈ RP, k ∈ K[rp]
)
)

Expand All @@ -73,7 +73,7 @@ function create_model(graph, params, sets; verbose = false, write_lp_file = fals
# - consumer balance equation
@constraint(
model,
c_consumer_balance[a ∈ Ac, rp ∈ RP, k ∈ K],
c_consumer_balance[a ∈ Ac, rp ∈ RP, k ∈ K[rp]],
sum(flow[f, rp, k] for f ∈ F if f[2] == a) -
sum(flow[f, rp, k] for f ∈ F if f[1] == a) ==
get(params.assets_profile, (a, rp, k), 1.0) * params.peak_demand[a]
Expand All @@ -84,7 +84,7 @@ function create_model(graph, params, sets; verbose = false, write_lp_file = fals
# TODO: Fix the initial storage_level
@constraint(
model,
c_storage_balance[a ∈ As, rp ∈ RP, k ∈ K],
c_storage_balance[a ∈ As, rp ∈ RP, k ∈ K[rp]],
storage_level[a, rp, k] ==
(k ≥ 2 ? storage_level[a, rp, k-1] : 0.0) +
sum(flow[f, rp, k] * params.flows_efficiency[f] for f ∈ F if f[2] == a) -
Expand All @@ -94,15 +94,15 @@ function create_model(graph, params, sets; verbose = false, write_lp_file = fals
# - hub balance equation
@constraint(
model,
c_hub_balance[a ∈ Ah, rp ∈ RP, k ∈ K],
c_hub_balance[a ∈ Ah, rp ∈ RP, k ∈ K[rp]],
sum(flow[f, rp, k] for f ∈ F if f[2] == a) ==
sum(flow[f, rp, k] for f ∈ F if f[1] == a)
)

# - conversion balance equation
@constraint(
model,
c_conversion_balance[a ∈ Acv, rp ∈ RP, k ∈ K],
c_conversion_balance[a ∈ Acv, rp ∈ RP, k ∈ K[rp]],
sum(flow[f, rp, k] * params.flows_efficiency[f] for f ∈ F if f[2] == a) ==
sum(flow[f, rp, k] / params.flows_efficiency[f] for f ∈ F if f[1] == a)
)
Expand All @@ -111,7 +111,7 @@ function create_model(graph, params, sets; verbose = false, write_lp_file = fals
# - overall output flows
@constraint(
model,
c_overall_output_flows[a ∈ Acv∪As∪Ap, rp ∈ RP, k ∈ K],
c_overall_output_flows[a ∈ Acv∪As∪Ap, rp ∈ RP, k ∈ K[rp]],
sum(flow[f, rp, k] for f ∈ F if f[1] == a) ≤
get(params.assets_profile, (a, rp, k), 1.0) * (
params.assets_init_capacity[a] +
Expand All @@ -122,7 +122,7 @@ function create_model(graph, params, sets; verbose = false, write_lp_file = fals
# - overall input flows
@constraint(
model,
c_overall_input_flows[a ∈ As, rp ∈ RP, k ∈ K],
c_overall_input_flows[a ∈ As, rp ∈ RP, k ∈ K[rp]],
sum(flow[f, rp, k] for f ∈ F if f[2] == a) ≤
get(params.assets_profile, (a, rp, k), 1.0) * (
params.assets_init_capacity[a] +
Expand All @@ -137,7 +137,7 @@ function create_model(graph, params, sets; verbose = false, write_lp_file = fals
a ∈ A,
f ∈ F,
rp ∈ RP,
k ∈ K;
k ∈ K[rp];
!(a ∈ Ah ∪ Ac) && f[1] == a && f ∉ Ft,
],
flow[f, rp, k] ≤
Expand All @@ -150,35 +150,35 @@ function create_model(graph, params, sets; verbose = false, write_lp_file = fals
# Constraints that define a lower bound for flows that are not transport assets
@constraint(
model,
c_lower_bound_asset_flow[f ∈ F, rp ∈ RP, k ∈ K; f ∉ Ft],
c_lower_bound_asset_flow[f ∈ F, rp ∈ RP, k ∈ K[rp]; f ∉ Ft],
flow[f, rp, k] ≥ 0
)

# Constraints that define bounds for a transport flow Ft
@expression(
model,
e_upper_bound_transport_flow[f ∈ F, rp ∈ RP, k ∈ K],
e_upper_bound_transport_flow[f ∈ F, rp ∈ RP, k ∈ K[rp]],
get(params.flows_profile, (f, rp, k), 1.0) * (
params.flows_init_capacity[f] +
(f ∈ Fi ? (params.flows_export_capacity[f] * flows_investment[f]) : 0.0)
)
)
@constraint(
model,
c_transport_flow_upper_bound[f ∈ Ft, rp ∈ RP, k ∈ K],
c_transport_flow_upper_bound[f ∈ Ft, rp ∈ RP, k ∈ K[rp]],
flow[f, rp, k] ≤ e_upper_bound_transport_flow[f, rp, k]
)
@expression(
model,
e_lower_bound_transport_flow[f ∈ F, rp ∈ RP, k ∈ K],
e_lower_bound_transport_flow[f ∈ F, rp ∈ RP, k ∈ K[rp]],
get(params.flows_profile, (f, rp, k), 1.0) * (
params.flows_init_capacity[f] +
(f ∈ Fi ? (params.flows_import_capacity[f] * flows_investment[f]) : 0.0)
)
)
@constraint(
model,
c_transport_flow_lower_bound[f ∈ Ft, rp ∈ RP, k ∈ K],
c_transport_flow_lower_bound[f ∈ Ft, rp ∈ RP, k ∈ K[rp]],
flow[f, rp, k] ≥ -e_lower_bound_transport_flow[f, rp, k]
)

Expand All @@ -193,7 +193,7 @@ function create_model(graph, params, sets; verbose = false, write_lp_file = fals
)
@constraint(
model,
upper_bound_for_storage_level[a ∈ As, rp ∈ RP, k ∈ K],
upper_bound_for_storage_level[a ∈ As, rp ∈ RP, k ∈ K[rp]],
storage_level[a, rp, k] ≤
params.initial_storage_capacity[a] + (a ∈ Ai ? energy_limit[a] : 0.0)
)
Expand Down