【需求响应】一种新的需求响应机制DR-VCG研究

简介: 【需求响应】一种新的需求响应机制DR-VCG研究

1 主要内容

该程序对应文章《Contract Design for Energy Demand Response》,电力系统需求响应(DR)用来调节用户对电能的需求,即在预测的需求高于电能供应时,希望通过需求响应减少用户用电,从而满足系统平衡。程序实现新的需求响应模型DR-VCG,该模型提供了灵活的用户参与DR过程合同,并且通过投标活动保证收益分配和价格计算的合理性。通过实例证实该方法的有效性,可靠性显著提升,总费用明显降低。该程序采用python编写。

2 部分代码

import grid
import agent
import contract
import matplotlib.pyplot as plt
import time
import statistics
import pandas as pd
from datetime import datetime
def main():
    for i in range(1):
        M = 1000
        number_of_agents = 200
        number_of_simulation_per_lanbda = 100
        generator_price_multiply = 1
        gamma = [1.0,1.166,1.333,1.5,1.666,1.833,2]
        df_columns = ['actuel expanse','gamma','M','actuel kWh reduced','Met the demand']
        row_data_df = pd.DataFrame(columns=df_columns)
        Fixed_cont_avg_cost = []
        Fixed_cont_avg_reliability = []
        Fixed_single_cont_avg_cost = []
        Fixed_single_cont_avg_reliability = []
        T_F_List_Fixed_cont_Met_the_demand = []
        Fixed_cont_Total_expense = []
        gamma_used = []
        for lb in gamma:
            Fixed_cont_reduce_list = []
            for i in range(number_of_simulation_per_lanbda):
                start = time.time()
                print('iteration:',i)
                Grid = grid.grid(M,lb)
                Grid.introduce_self()
                Agents = []
                for num in range(number_of_agents):
                    ag = agent.agent(num)
                    Agents.append(ag)
                Contracts = []
                for i in range(10, M+1, 10):
                    Contracts.append(contract.contract(i,0.3,0.5))
                single_contract = []
                single_contract.append(contract.contract(50,0.3,0.5))
                Grid.set_single_contract(single_contract)
                Grid.set_contract(Contracts)
                Grid.set_agents(Agents)
                Grid.send_contrects_to_agents()
                Grid.send_single_contrects_to_agents()
                for ag in Agents:
                    ag.Fixed_cont_bid_on_contract()
                for ag in Agents:
                    ag.Fixed_single_cont_bid_on_contract()
                Grid.Fixed_cont_get_bids_from_agent()
                Grid.Fixed_cont_generator_bids(price_multiply=generator_price_multiply)
                Grid.Fixed_cont_get_q_from_agent()
                Fixed_cont_sum_of_bids = Grid.knapsack(bids_type='Fixed_cont')
                Grid.Fixed_cont_pay_to_agents(Fixed_cont_sum_of_bids)
                Grid.Fixed_cont_reliability()
                Grid.Fixed_single_cont_get_bids_from_agent()
                Grid.Fixed_single_cont_get_q_from_agent()
                Fixed_cont_Total_expense.append(Grid.Fixed_cont_Total_expense_sum)
                Fixed_cont_reduce_list.append(Grid.Fixed_cont_reliability_sum_q)
                if Grid.Fixed_cont_reliability_sum_q >= Grid.M:
                    met_the_demand = 1
                else:
                    met_the_demand = 0
                T_F_List_Fixed_cont_Met_the_demand.append(met_the_demand)
                print('Fixed_cont- Met_the_demand: ', T_F_List_Fixed_cont_Met_the_demand)
                print('Fixed_cont- Total_expense: ', Fixed_cont_Total_expense)
                gamma_used.append(lb)
                row_data_df = row_data_df.append(pd.DataFrame(
                                                 {'actuel expanse':[Grid.Fixed_cont_Total_expense_sum],
                                                  'gamma':[lb],
                                                  'M': [M],
                                                  'actuel kWh reduced': [Grid.Fixed_cont_reliability_sum_q],
                                                  'Met the demand': [met_the_demand]}))
                end = time.time()
                print('iteration took:', (end - start), 'sec')
                print('-'*200)
            Fixed_cont_avg_cost.append(statistics.mean(Fixed_cont_Total_expense))
            if len(T_F_List_Fixed_cont_Met_the_demand) > 0:
                Fixed_cont_avg_reliability.append(T_F_List_Fixed_cont_Met_the_demand.count(True) / len(T_F_List_Fixed_cont_Met_the_demand))
            else:
                Fixed_cont_avg_reliability.append(0.0)
        filename = datetime.now().strftime('data/energy_demamd_row_data-%Y-%m-%d-%H-%M-%S.csv')
        row_data_df.to_csv(filename,index=False)
        graph_it(Fixed_cont_avg_reliability,
                 Fixed_single_cont_avg_reliability, Fixed_cont_avg_cost,
                 Fixed_single_cont_avg_cost)
def graph_it(Fixed_cont_avg_reliability =[],
             Fixed_single_cont_avg_reliability=[],
             Fixed_cont_avg_cost=[],
             Fixed_single_cont_avg_cost=[]):
    plt.rcParams["figure.figsize"] = (8, 8)
    fig, ax = plt.subplots()
    ax.plot(Fixed_cont_avg_reliability, Fixed_cont_avg_cost, color='blue',marker='o',label="fixed multiple cont")
    ax.plot(Fixed_single_cont_avg_reliability, Fixed_single_cont_avg_cost, color='black',marker='o', label="fixed single cont")
    ax.set(xlabel="Total_Reliability", ylabel="expenses ($)", title="(a)n= 400")
    fig.savefig("test.png")
if __name__ == "__main__":
    main()
    plt.show()


3 程序结果

原文结果图:

4 程序链接

关注+私信免费获取程序源码和资料!

相关文章
|
Java 数据库 Maven
根据请求中接收到的主动协商头字段,目标资源没有用户代理可以接受的当前表示
今天使用ssm框架搭建了一个项目,测试一个插入的接口结果返回406,但是数据库插入成功了,报错406是因为返回的结果接口无法解析。
855 0
|
1月前
|
存储 缓存 监控
|
2月前
|
Python
8. 如何解决 Tornado 检测到了有事件(events)被发送到一个已经关闭的流(stream)。在 Tornado 中,一个流代表一个请求或响应的数据流。这个警告可能意味着在请求处理的过程中,
8. 如何解决 Tornado 检测到了有事件(events)被发送到一个已经关闭的流(stream)。在 Tornado 中,一个流代表一个请求或响应的数据流。这个警告可能意味着在请求处理的过程中,
|
3月前
|
Web App开发 监控 安全
[译] 用 sendBeacon 发送分析信息的优点
[译] 用 sendBeacon 发送分析信息的优点
|
6月前
|
监控 安全 持续交付
【专栏】Webhook是服务器主动发送事件通知的机制,打破传统客户端轮询模式,实现数据实时高效传递。
【4月更文挑战第29天】Webhook是服务器主动发送事件通知的机制,打破传统客户端轮询模式,实现数据实时高效传递。常用于持续集成部署、第三方服务集成、实时数据同步和监控告警。具有实时性、高效性和灵活性优势,但也面临安全风险和调试挑战。理解并善用Webhook能提升系统性能,广泛应用于现代软件开发和集成。
409 0
|
6月前
|
存储 缓存 前端开发
如何实现设备组缓存的正确清除?——基于心跳请求和心跳响应的解决方案
如何实现设备组缓存的正确清除?——基于心跳请求和心跳响应的解决方案
73 0
|
6月前
|
缓存 调度
响应系统的作用与实现 (下)
响应系统的作用与实现
65 0
|
6月前
|
存储 JavaScript Java
响应系统的作用与实现(上)
响应系统的作用与实现
83 0
|
Python
【需求响应DR】一种新的需求响应机制DR-VCG研究(Python代码实现)
【需求响应DR】一种新的需求响应机制DR-VCG研究(Python代码实现)
|
JSON PHP 数据格式
响应 方式
响应 方式