实时电价机制下交直流混合微网优化运行方法（Matlab代码实现）

👨‍🎓个人主页：研学社的博客 

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

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

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

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

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

---

💥1 概述

文献来源：

摘要：为了提高交直流混合微网经济效益,提出实时电价机制下的交直流混合微网优化运行方法。首先针对交直流混合微网的结构和电气特性建立其优化运行模型;其次针对微网中负荷构成建立了基于负荷分类的负荷水平对实时电价的响应模型。所提出的日前调度模型将交直流混合微网制定的网内实时电价作为控制变量,给出实时电价制定策略和交直流混合微网日前调度方案。对所建立的模型采用混沌粒子群算法求解。最后通过一个算例验证了所提出方法可以提高风光消纳率,增加交直流混合微网收益并且减小用户平均电费支出。

关键词：

交直流混合微网;优化运行;实时电价;需求侧响应;混沌粒子群算法;

📚2 运行结果

 部分代码“

function load=DR3(pload,detapr)
%pload=ones(1,24);
price=[0.4042,0.4042,0.4042,0.4042,0.4042,0.4042,0.4042,0.8084,0.8084,1.3339,1.3339,1.3339,1.3339,1.3339,1.3339,0.8084,0.8084,1.3339,1.3339,1.3339,1.3339,0.8084,0.4042,0.4042];
%detapr=[-0.1.*ones(1,12) 0.1.*ones(1,12)];
pl1=0.6.*pload;%刚性负荷
pl2=0.2.*pload;%Ⅰ类
pl3=0.2.*pload;%Ⅱ类
[row,colz]=find(detapr>0);%电价上升时段
[row,colf]=find(detapr<0);%电价下降时段
num_colz=size(colz,2);%电价上升时段数量
num_colf=size(colf,2);%%电价下降时段数量
T=24;
a=0.08;%死区阈值
%Ⅰ类负荷
linj=0;lini=0;
for i=1:num_colf
    linj=linj+abs(detapr(colf(i)));%公式（2）
end
for i=1:num_colz
    lini=lini+abs(detapr(colz(i)));
end
for t=1:T
    if ismember(t,colz)==1
        fdetap=fhzy(detapr(t),a);
        linp=0;
        for kj=1:num_colf
            linp=linp+fdetap*abs(detapr(colf(kj)))*pl2(t)/linj;
        end
         lam_pl2(t)=pl2(t)-linp;%公式（5）
       elseif ismember(t,colf)==1
        fdetap=fhzy(detapr(t),a);
        linp=0;
        for ki=1:num_colz
            linp=linp+fdetap*abs(detapr(colz(ki)))*pl2(t)/lini;
        end
         lam_pl2(t)=pl2(t)+linp;%公式（6）
    end
end
%Ⅱ类负荷
e0=-1;e1=0.8;e2=0.6;
Est=zeros(T,T);
for i=1:T
    for j=1:T
        if  i==j
            Est(i,j)=-1;
        elseif i==j+1
            Est(i,j)=0.8;
        elseif i+1==j
            Est(i,j)=0.6;
        else
            Est(i,j)=0;
        end
    end
end
for t=1:T
    lam_pl3(t)=pl3(t)+pl3(t)*(Est(t,:)*(detapr./price)');%公式（10）
end
load=pl1+lam_pl2+lam_pl3;

🎉3 文献来源

部分理论来源于网络，如有侵权请联系删除。

[1]陈安伟,华浩瑞,李鹏等.实时电价机制下交直流混合微网优化运行方法[J].电力系统保护与控制,2017,45(07):13-20.

🌈4 Matlab代码实现