👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆下载资源链接👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆
《《《《《《《《更多资源还请持续关注本专栏》》》》》》》
论文与完整源程序_电网论文源程序的博客-CSDN博客https://blog.csdn.net/liang674027206/category_12531414.html
以可再生能源制氢为特征的能源系统将是今后能源互联网建设的重要方向之一。该文针对风–光–氢多主体能源系统的合作运行展开研究。首先,考虑主体间的电能交易建立各主体的优化运行模型,然后基于纳什谈判理论建立风–光氢多主体合作运行模型,接着将其等效为联盟效益最大化和电能交易支付谈判两个子问题。为保护各主体隐私,运用交替方向乘子法提出上述两个子问题的分布式求解方法。最后通过算例验证所提合作运行模型以及分布式算法的有效性。仿真结果表明通过风–光–氢多主体的合作运行,可以较大幅度提高各主体的运行效益以及合作联盟的整体效益。此外,风光发电上网电价的降低,将促进风–光–氢多主体展开合作,以提升各自运行效益。
部分代码展示:
function [P_pv2h_2,P_wt2h_2,Obj_dzq]=Fun_dzq(value_P_wt2h_2,value_P_pv2h_2,rho_wt,rho_pv,lambda_wt,lambda_pv)
%电制氢主体dzq的分布式优化迭代模型
%% 决策变量初始化
P_pv2h_2=sdpvar(1,24); %光伏主体向氢气主体售电量(电制氢主体所期望的)
P_wt2h_2=sdpvar(1,24); %风电主体向氢气主体售电量(电制氢主体所期望的)
P_el=sdpvar(1,24); %产氢对应的耗电量
P_H2=sdpvar(1,24); %电制氢主体的产氢量
P_com=sdpvar(1,24); %氢气压缩机消耗电功率
m_com=sdpvar(1,24); %压缩机压缩氢气流量
P_hg=sdpvar(1,24); %电制氢主体从大电网购电量
E_bat=sdpvar(1,24); %储电容量
P_batc=sdpvar(1,24); %充电功率
P_batd=sdpvar(1,24); %放电功率
U_abs=binvar(1,24); %储电设备的放电状态位,取1时为放电,0为未放电
U_relea=binvar(1,24); %储电设备的充电状态位,取1时为充电,0为未充电
%% 导入氢能负荷和电网购电电价
L_H2=[13.1632653100000,12.6530612200000,13.6734693900000,10.6122449000000,16.1224489800000,27.9591836700000,49.7959183700000,51.0204081600000,47.5510204100000,46.5306122400000,47.1428571400000,48.3673469400000,49.1836734700000,51.2244898000000,62.4489795900000,69.5918367300000,56.3265306100000,42.6530612200000,33.0612244900000,26.1224489800000,20.6122449000000,17.9591836700000,16.3265306100000,14.4897959200000]; %单位为kg
p_hg=[0.3376*ones(1,7),0.5980*ones(1,4),0.8654*ones(1,3),0.5980*ones(1,4),0.8654*ones(1,4),0.3376*ones(1,2)];
%% 约束条件
C=[];
C=[C,
P_H2==0.019224*P_el, %式(15):产氢约束
0<=P_el<=5000, %式(16)中的最大功率约束
P_com==14.304*m_com*40*1.4*[(3)^(1-1/1.4)-1]/[3600*0.7*(1.4-1)], %式(17):压缩机电功率
0<=P_com<=3000, %式(18):压缩机耗电功率不等式约束 %式(19):压缩机内部气压连续约束(后面补充)
20<=P_H2<=35, %式(20):压缩机内部气压不等式约束
P_H2(1)==25, %式(20):初始气压 式(21)电储能模型(后面补充)
P_hg+P_wt2h_2+P_pv2h_2+P_batd==P_el+P_batc+P_com, %式(25):氢气主体功率平衡
0<=P_hg, %补充定义:氢气主体不能向电网售电,避免光电主体的售电收益被转移到氢气主体侧
];
for t=2:24
C=[C,
-1000<= P_el(t)-P_el(t-1)<=1000, %式(16):爬坡约束
P_H2(t)==P_H2(t-1)+(m_com-L_H2)*858, %式(19):压缩机内部气压连续约束
];
end
%式(21):储电约束
%储能电站荷电状态连续性约束
C=[C,E_bat(1)==500+0.95*P_batc(1)-P_batd(1)/0.96,]; %1时段约束
for t=2:24
C=[C,E_bat(t)==E_bat(t-1)+0.95*P_batc(t)-P_batd(t)/0.96,]; %储电设备容量变化约束
end
效果展示:
👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆下载资源链接👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆
《《《《《《《《更多资源还请持续关注本专栏》》》》》》》
论文与完整源程序_电网论文源程序的博客-CSDN博客https://blog.csdn.net/liang674027206/category_12531414.html