基于合作博弈的综合能源系统利益分配优化调度（Matlab代码实现）

 💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码、数据

---

💥1 概述

摘要:"碳达峰、碳中和"重大战略决策的提出，进一步的推动了新能源的发展，同时风能与光能具有不稳定性，因此构建综合能源系统对新能源的消纳具有重要意义。本文对综合能源系统的合作运行进行了研究，在考虑系统运营成本、碳排放、利益分配的基础上，首先构建综合能源系统框架模型，再对系统内三个主体构建合作联盟，采用改进的NSGA-Il算法（或者cplex求解器）求解．并提出一种改进的Shapley值法对综合能源系统联盟中各成员进行利益分配。最后通过算例进行模拟，结果证明所提模型能够有效降低系统总用能成本提高收益的同时减少碳排放量，改进的Shapley值法可对联盟收益分配起到很好的效果。
关键词:合作博弈;综合能源系统;碳排放;利益分配;
1引言
近年来全球能源危机正在逐渐加剧，二氧化碳的大量排放对环境造成了严重的破坏。十四五规划明确提出实现"碳达峰、碳中和"重大战略决策的首要途径便是构建绿色能源体系，进一步推动了新能源的发展。
风能和太阳能作为两大重要的新能源，具有绿色、可再生、储量丰富等天然优势。但受自然气象变化等因素影响，风能和太阳能存在不稳定、不连续的特点，随着发电侧集中式新能源的大规模接入，电力系统出力难以快速跟随用电负荷，其安全稳定运行面临巨大挑战。同时目前我国正逐步减少对风光发电的补贴，积极推进风光发电平价上网，这将对风光发电的效益造成一定的影响。所以构建综合能源系统协同优化、互济互补的运行模式为上述问题提供了一种解决思路,且已成为当下的研究热点。
综合能源系统(integrated energy system. IES)是一种将电、热、气、氢等多种异质能源相互耦合的供能系统."横向多能互补、纵向源-网―荷-储协调"的特性可显著提高风光消纳、降低系统运行成本.减少二氧化碳排放，增加联盟收益降低风险。

📚2 运行结果

部分代码：

P_wt=value(P_wt);     %风机消纳功率
P_pv=value(P_pv);      %光伏消纳功率
V_gas=value(V_gas);     %P2G产生天然气体积
Q_co2=value(Q_co2);     %制气厂消耗二氧化碳的量
P_H2=value(P_H2);      %电解池产生的氢气量
P_P2G=value(P_P2G);     %P2G消耗电功率量
P_h2cha=value(P_h2cha);   %储氢罐充气功率
P_h2dis=value(P_h2dis);   %储氢罐放气功率
H_h2cha=value(H_h2cha);   %储氢罐充气变量
H_h2dis=value(H_h2dis);   %储氢罐放气变量
P_H2total=value(P_H2total);   
E_h=value(E_h);       %储氢罐放气功率
P_CHPe=value(P_CHPe);     %CHP发电功率
P_CHPh=value(P_CHPh);     %CHP发热功率
P_CHPg=value(P_CHPg);     %CHP耗气功率
Q_co2out=value(Q_co2out);   %CHP燃烧产生的co2
P_be=value(P_be);       %购电功率
P_bg=value(P_bg);       %购气功率
V_CB=value(V_CB);       %燃气锅炉进气量
Q_CB=value(Q_CB);       %燃气锅炉发热量
P_MRh=value(P_MRh);       %甲烷反应器输入氢功率
P_MRg=value(P_MRg);       %甲烷反应器输出天然气功率
F=value(F);

%% 绘图
figure(1);
plot(load_e,'linewidth',1.5);
hold on
plot(load_h,'-*','linewidth',1.5);
hold on
plot(load_g,'--','linewidth',1.5);
hold on
plot(load_H,'-^','linewidth',1.5);
hold on
axis([1,24,0,1000]);
 title('典型日负荷曲线');
 legend('电负荷','热负荷','气负荷','氢负荷');
 xlabel('时刻');
 ylabel('功率/kW');
 grid on

figure(2);
plot(P_wt,'r--','linewidth',1.5);
hold on
plot(P_pv,'b-*','linewidth',1.5);
hold on
axis([1,24,0,1000]);
 legend('风机出力','光伏出力');
 xlabel('时刻');
 ylabel('功率/kW');
 grid on
 
figure(3);
stairs(C_e,'r--','linewidth',1.5);
hold on
stairs(C_g,'b-*','linewidth',1.5);
hold on
axis([1,24,0,4]);
 legend('购电电价','购气价格');
 xlabel('时刻');
 ylabel('功率/kW');
 grid on
 
 
figure(4)
F11=zeros(2,24);
F11(1,:)=load_e;
F11(2,:)=P_P2G;
AA=bar(F11(:,:)','stack');
set(AA(1),'FaceColor',[0.1 0.5 0.9]);
set(AA(2),'FaceColor',[0.9 0.1 0.5]);
hold on
FF=zeros(4,24);
FF(1,:)=P_wt;
FF(2,:)=P_pv;
FF(3,:)=P_CHPe;
FF(4,:)=P_be;
bar(-FF(:,:)','stack')
axis([0.5,24.5,-1500,1500]);
xlabel('时刻/h');
ylabel('电功率/kW'); 
legend('电负荷','P2G消耗电量','风机出力','光伏出力','CHP发电功率','购电功率');
title('电负荷平衡');
grid on

figure(5)
F22=zeros(3,24);
F22(1,:)=load_g;
F22(2,:)=P_CHPg;
F22(3,:)=V_CB;
BB=bar(F22(:,:)','stack');
set(BB(1),'FaceColor',[0.1 0.5 0.9]);
set(BB(2),'FaceColor',[0.9 0.1 0.5]);
set(BB(2),'FaceColor',[0.2 0.4 0.2]);
hold on
FFF=zeros(2,24);
FFF(1,:)=P_MRg;
FFF(2,:)=P_bg;
bar(-FFF(:,:)','stack')
axis([0.5,24.5,-600,600]);
xlabel('时刻/h');
ylabel('气功率/kW'); 
legend('气负荷','CHP耗气量','燃气锅炉进气量','甲烷反应器输出天然气','购气功率');
title('气负荷平衡');
grid on
%% 合作博弈处理

F1=sc1();
F2=sc2();
F3=sc3();
F12=sc12();
F13=sc13();
F23=sc23();

s=[1,0,0
    0,1,0
    0,0,1
    1,1,0
    1,0,1
    0,1,1
    1,1,1
];
v=[F1 F2 F3 F12 F13 F23 F];
abcd=shapley1(s,v);

clc
disp(['********************']);

disp(['shapley分配前的收益为：']);

disp(['   新能源电厂收益为 ',num2str(F1),'元']);

disp(['   燃气热电厂收益为 ',num2str(F2),'元']);

disp(['   电制气厂收益为 ',num2str(F3),'元']);

disp(['********************']);
disp(['********************']);
disp(['shapley分配后的收益为：']);

disp(['   新能源电厂收益为 ',num2str(abcd(1)),'元']);

disp(['   燃气热电厂收益为 ',num2str(abcd(2)),'元']);

disp(['   电制气厂收益为 ',num2str(abcd(3)),'元']);
disp(['********************']);

🎉3 参考文献

文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。

[1]王永利,刘振,薛露等.基于多主体博弈的综合能源系统运行优化方法[J].控制理论与应用,2022,39(03):499-508.

[2]祝荣,任永峰,孟庆天等.基于合作博弈的综合能源系统电-热-气协同优化运行策略[J].太阳能学报,2022,43(04):20-29.DOI:10.19912/j.0254-0096.tynxb.2022-0112

🌈4 Matlab代码、数据