Pumped_Storage.py 文件源码

def obj_func_cont(xx, e_g, e_p, g, rho, Q_g, Q_p, head_g, head_p, optimizing = True):
    H_T = int(price_duration.Frequency.max()) # total duration (100%)
    x1 = np.arange(0,xx)
    y1 = f(x1)
    x2 = np.arange(H_T-xx,H_T)
    y2 = f(x2)
    Power_Revenue = np.trapz(y1, x1, dx=0.1, axis = -1)*e_g*rho*g*Q_g*head_g/(10**6)
    Pumping_Cost = np.trapz(y2, x2, dx=0.1, axis = -1)/e_p*rho*g*Q_p*head_p/(10**6)   
    z = Power_Revenue - Pumping_Cost # profit
    return -z if optimizing else z

# objective function to maximize - discrete