具体问题是Applied Energy期刊文章《Combined heating and cooling networks with waste heat recovery based on energy hub concept 》中的能源枢纽问题(该问题是多目标问题,这里主要是复现这篇文章)
from gurobipy import *
import random
import numpy as np
try:
# 天然气锅炉的转换率(天然气->热负荷)
eta_fh_b = 0.85
# 热泵的转换率(电->热负荷;(eta_eh_hp-1)为电->冷负荷)
eta_eh_hp = 6
# 制冷机组的转换率(电->冷负荷)
eta_eh_ch = 4
# 电锅炉的转换率(电->热负荷)
eta_eh_eh = 0.95
# 热泵的价格(热泵的容量,每1kw需要的钱)
phi_hp = 230
# 电的价格(每一度(1kwh)需要的钱)
phi_e = 0.0327
# 天然气的价格(每m³的天然气价格,产生一度电的天然气的价格)
phi_f = 0.016123
# 电排放二氧化碳的比例(每度电排放的二氧化碳量)
lambda_e = 0.02072
# 天然气排放二氧化碳的比例(产生一度电的天然气排放的二氧化碳量)
lambda_f = 0.17644
# 热负荷函数的振幅
A_h = 50000
# 冷负荷函数的振幅
A_c = 50000