Source code for sdom.models.formulations_hydro

from pyexpat import model
from pyomo.environ import Set, Param, value, NonNegativeReals
from .models_utils import add_generation_variables, add_alpha_and_ts_parameters, add_budget_parameter, add_upper_bound_paramenters, add_lower_bound_paramenters, generic_budget_rule
from pyomo.core import Var, Constraint, Expression

from ..constants import VALID_HYDRO_FORMULATIONS_TO_BUDGET_MAP, MONTHLY_BUDGET_HOURS_AGGREGATION, DAILY_BUDGET_HOURS_AGGREGATION
from ..io_manager import get_formulation


####################################################################################|
# ----------------------------------- Parameters -----------------------------------|
####################################################################################|


[docs]
def add_large_hydro_parameters(model, data: dict):
    add_alpha_and_ts_parameters(model.hydro, model.h, data, "AlphaLargHy", "large_hydro_data", "LargeHydro")
    formulation = get_formulation(data, component='hydro')
    add_budget_parameter(model.hydro, formulation, VALID_HYDRO_FORMULATIONS_TO_BUDGET_MAP)




[docs]
def add_large_hydro_bound_parameters(model, data: dict):
    # Time-series parameter data initialization
    add_upper_bound_paramenters(model.hydro, model.h, data, "large_hydro_max", "LargeHydro")
    add_lower_bound_paramenters(model.hydro, model.h, data, "large_hydro_min", "LargeHydro")




####################################################################################|
# ------------------------------------ Variables -----------------------------------|
####################################################################################|


[docs]
def add_hydro_variables(model):
    add_generation_variables(model.hydro, model.h, domain=NonNegativeReals, initialize=0)


####################################################################################|
# ----------------------------------- Constraints ----------------------------------|
####################################################################################|


[docs]
def add_hydro_run_of_river_constraints(model, data: dict):
    model.hydro.run_of_river_constraint = Constraint(model.h, rule=lambda m,h: m.generation[h] == m.alpha * m.ts_parameter[h] )
    return




[docs]
def add_hydro_budget_constraints(model, data: dict):
    
    model.hydro.upper_bound_ts_constraint = Constraint(model.h, rule=lambda m,h: m.generation[h] <= m.alpha * m.ts_parameter_upper_bound[h] )
    model.hydro.lower_bound_ts_constraint = Constraint(model.h, rule=lambda m,h: m.generation[h] >= m.alpha * m.ts_parameter_lower_bound[h] )

    model.hydro.budget_constraint = Constraint(model.hydro.budget_set, rule = generic_budget_rule )
    
    return