def minimize_cost(self, c_blocks=(), r_blocks=()):
    """ Builds objective function that minimises the total costs.

    Costs included are:
                        opex_var,
                        opex_fix,
                        curtailment_costs (dispatch sources),
                        annualised capex (investment components)

    Parameters:
    ------------
    self : pyomo model instance
    c_blocks : pyomo blocks containing components that are included in
               cost terms of objective function
    r_blocks : pyomo blocks containing components that are included in revenue
               terms of objective function
    """
    expr = 0
    c_blocks = (str(transformer.Simple), str(transformer.CHP), str(source.FixedSource))
    r_blocks = ()

    blocks = [
        block
        for block in self.block_data_objects(active=True)
        if not isinstance(block, solph.optimization_model.OptimizationModel)
    ]

    for block in blocks:
        if block.name in c_blocks:
            if block.name == str(transformer.Storage):
                ref = "capacity"
                expr += objexpr.add_opex_var(self, self.simple_storage, ref="input")
            else:
                ref = "output"
            # variable costs

            expr += objexpr.add_opex_var(self, block, ref="output")
            # fix costs
            expr += objexpr.add_opex_fix(self, block, ref=ref)
            # investment costs
            if block.optimization_options.get("investment", False):
                expr += objexpr.add_capex(self, block, ref=ref)
            # revenues
        if block.name in r_blocks:
            expr += objexpr.add_revenues(self, block, ref="output")

    # costs for dispatchable sources
    if hasattr(self, str(source.DispatchSource)):
        expr += objexpr.add_curtailment_costs(self, getattr(self, str(source.DispatchSource)))

    if getattr(self, str(Bus)).shortage_uids:
        expr += objexpr.add_shortage_slack_costs(self, block)
    # artificial costs for excess or shortage
    if getattr(self, str(Bus)).excess_uids:
        expr += objexpr.add_excess_slack_costs(self, block)

    self.objective = po.Objective(expr=expr)
Exemplo n.º 2
0
def minimize_cost(self, cost_objects=None, revenue_objects=None):
    """ Builds objective function that minimises the total costs.

    Costs included are:
                        opex_var,
                        opex_fix,
                        curtailment_costs (dispatch sources),
                        annualised capex (investment components)

    Parameters
    ----------
    self : pyomo model instance
    cost_blocks : array like
       list containing classes of objects that are included in
               cost terms of objective function
    revenue_blocks : array like
       list containing classes of objects that are included in revenue
               terms of objective function
    """
    expr = 0
    # TODO: Fix naming
    c_blocks = cost_objects
    r_blocks = revenue_objects

    if cost_objects is None:
        c_blocks = [str(transformer.Simple),
                    str(transformer.TwoInputsOneOutput),
                    str(transformer.VariableEfficiencyCHP),
                    str(transformer.SimpleExtractionCHP),
                    str(transformer.Storage),
                    str(transformer.CHP),
                    str(source.FixedSource),
                    str(source.Commodity)]
    if revenue_objects is None:
        r_blocks = []

    blocks = [block for block in self.block_data_objects(active=True)
              if not isinstance(block,
                                solph.optimization_model.OptimizationModel)]

    for block in blocks:
        if block.name in c_blocks:
            if block.name == str(transformer.Storage):
                ref = 'capacity'
                expr += objexpr.add_opex_var(self,
                                             getattr(self,
                                                     str(transformer.Storage)),
                                             ref='input')
            else:
                ref = 'output'
            # variable costs
            expr += objexpr.add_opex_var(self, block, ref='output')
            # fix costs
            if block != str(source.Commodity):
                expr += objexpr.add_opex_fix(self, block, ref=ref)
            # investment costs
            if block.optimization_options.get('investment', False):
                expr += objexpr.add_capex(self, block, ref=ref)
            if hasattr(block, 'z_start'):
                expr += objexpr.add_startup_costs(self, block)

            # revenues
        if block.name in r_blocks:
            expr += objexpr.add_revenues(self, block, ref='output')

    # costs for dispatchable sources
    if hasattr(self, str(source.DispatchSource)):
        expr += \
            objexpr.add_curtailment_costs(self,
                                          getattr(self,
                                                  str(source.DispatchSource)))

    # artificial costs for excess or shortage
    if hasattr(self, str(ExcessSlack)):
        expr += objexpr.add_excess_slack_costs(self,
                                               getattr(self, str(ExcessSlack)))
    if hasattr(self, str(ShortageSlack)):
        expr += objexpr.add_shortage_slack_costs(self,
                                                 getattr(self, str(ShortageSlack)))


    self.objective = po.Objective(expr=expr)
 def objective_function_expressions():
     objfuncexprs.add_opex_var(om, block, ref='output')
     objfuncexprs.add_opex_fix(om, block, ref='output')
     objfuncexprs.add_curtailment_costs(om, block)
 def objective_function_expressions():
     objfuncexprs.add_opex_var(om, block, ref='output')
     objfuncexprs.add_opex_fix(om, block, ref='capacity')
 def objective_function_expressions():
     objfuncexprs.add_opex_var(om, block, ref='output')
     objfuncexprs.add_opex_fix(om, block, ref='output')
     objfuncexprs.add_input_costs(om, block)
     objfuncexprs.add_revenues(om, block)