的配电网重构》以考虑网损成本、弃风弃光成本和开关操作惩罚成本的综合成本最小为目标针对配电网重构模型的非凸性引入中间变量并对其进行二阶锥松弛构建混合整数凸规划模型采用改进的

IEEE33

节点配电网进行算例仿真分析了需求响应措施和清洁能源渗透率对配电网重构结果的影响。

该程序复现效果和出图较好详见程序结果部分注释清楚方便学习

注意该程序运行环境为matlabmosek需要各位同学下载并安装mosek求解器通过官网可以申请学术许可可免费使用365天。

目标函数

目标函数为配电网综合运行成本最小其中考虑了网损成本、弃风弃光成本以及分段开关操作惩罚成本。

重要约束条件

常规的功率平衡、节点电压电流等约束不再赘述重点分析一下网络结构约束和需求响应约束。

网络结构约束

该网络结构约束是采用虚拟潮流方式之前有几个重构代码也是采用虚拟潮流形式参考的是《A

New

仔细观察不难发现上面的模型是下面的简洁版在不考虑分布式电源节点对网络切割情况下两者是等价的。

经验证见结果图最后一张该种约束方式下能够保证网络的连通性和辐射性。

需求响应约束

在配电网中采用需求响应策略可以在降低负荷峰谷差的同时减少配电网运行的综合成本提高配电网运行的经济性和可靠性。

2部分代码

渗透率P_DGsum(mpc.bus(:,3))*F/mpc.baseMVA/5;

光伏最大有功P_loadmpc.bus(:,3)/mpc.baseMVA*P_L0;

有功负荷Q_loadmpc.bus(:,4)/mpc.baseMVA*P_L0;

支路功率最大值r_ijmpc.branch(:,3)*ones(1,T);

线路电阻x_ijmpc.branch(:,4)*ones(1,T);

线路电抗wind[9

储能接入位置Umax[1;1.06*1.06*ones(32,1)];

电压上限的平方Umin[1;0.94*0.94*ones(32,1)];

充电功率上限0.2MWP_dis_max0.2/mpc.baseMVA;

放电功率上限0.2MWE_min0.15/mpc.baseMVA;

储能容量下限0.15MWhE_max0.8/mpc.baseMVA;

储能容量上限0.8MWhn_ch0.9;

初始荷电状态为0.3MWhQ_CB_st0.15/mpc.baseMVA;

弹性系数c13;

分段开关操作惩罚成本系数15元/次rhozeros(1,24);

谷时电价rho00.35;

电源有功功率最大值Q_g_max10/mpc.baseMVA;

电源无功功率最大值branch_to_nodezeros(nb,nl);

流入节点的支路branch_from_nodezeros(nb,nl);

k1:nl

branch_to_node(mpc.branch(k,2),k)1;

%举例说明k1,流入节点2是支路1同时流出节点1的是支路1同理k2,流入节点3且流出节点2的是支路2这一步建立支路和节点的连接关系

branch_from_node(mpc.branch(k,1),k)1;end

储能荷电状态N_CBintvar(1);

1.潮流约束m_ij(1-alpha_ij)*M*ones(1,T);

Constraints

P_g-P_curbranch_to_node*P_ij-branch_to_node*(I_ij.*r_ij)-branch_from_node*P_ij

0];Constraints

Q_g-Q_loadbranch_to_node*Q_ij-branch_to_node*(I_ij.*x_ij)-branch_from_node*Q_ij

0];Constraints

[Constraints,U_i(mpc.branch(:,1),:)-U_i(mpc.branch(:,2),:)

m_ij

[Constraints,U_i(mpc.branch(:,1),:)-U_i(mpc.branch(:,2),:)

-m_ij

I_ij(k,t)-U_i(mpc.branch(k,1),t)],I_ij(k,t)U_i(mpc.branch(k,1),t))];

endendConstraints

P_g_dot-P_L_dotbranch_to_node*P_ij_dot-branch_from_node*P_ij_dot

0];

n_ch*P_ch(:,1)-1/n_dis*P_dis(:,1)E0];Constraints

[Constraints

n_ch*P_ch(:,t)-1/n_dis*P_dis(:,t)E_ESS(:,t-1)];

3程序结果