基于SOE的随机配电网重构研究（Matlab代码实现）

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

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

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

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

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

目录

💥1 概述

📚2 运行结果

🎉3 文献来源

🌈4 Matlab代码实现

---

💥1 概述

配电网重新配置（DNR）的目的是确定配电网的最佳拓扑结构，这是减少电网功率损耗的有效措施。电力负荷需求和光伏（PV）输出是不确定的，并且随一天中的时间而变化，并且会影响最佳网络拓扑。单小时确定性 DNR 无法处理这种不确定性和可变性。因此，本文提出求解多小时随机DNR（SDNR）。现有的DNR求解方法要么不准确，要么过于耗时，因此无法求解大型分配网络的数小时SDNR。对此，提出了一种开关开换（SOE）方法。SOE 从关闭所有交换机的环形网络开始，由三个步骤组成。第一步是按顺序打开开关，直到所有环路都打开。第二步和第三步修改第一步中获得的分支的状态，以获得更好的径向拓扑。使用五个测试系统来验证SOE的准确性和快速求解速度，以及多小时SDNR相对于单小时确定性DNR的优越性。

📚2 运行结果

 

  

  

 部分代码：

%% ------------------------ Tabu algorithm ------------------------%%
% run the core program for each upstream branch connected to the idx_force_open
% idx_considered = [35 69]
% for iter = idx_considered
for iter = 1:length(branch_idx_focused)
    fprintf('iter=%d/%d\n', iter, length(branch_idx_focused));
    Branch = Branch0;
    Branch(branch_idx_focused(iter), :) = [];
    
    ff0 = Branch(:, 1);   ff = ff0;
    tt0 = Branch(:, 2);   tt = tt0;
    
    brch_idx_in_loop = brch_idx_in_loop0;
    idx_tmp = find(brch_idx_in_loop == branch_idx_focused(iter));
    if isempty(idx_tmp)
    else
        brch_idx_in_loop(idx_tmp) = [];
        brch_idx_in_loop(idx_tmp:end) = brch_idx_in_loop(idx_tmp:end)-1;
    end

    t1 = toc;
    %%------------------- core algorithm in Tabu loop--------------------%%
    if version_LODF
        [Branch] = decrease_reconfig_algo_LODF(Bus, Branch, brch_idx_in_loop, ...
            ff0, tt0, substation_node, n_bus, loss0, distancePara); %%%  core algorithm
    else
        [Branch] = decrease_reconfig_algo(Bus, Branch, brch_idx_in_loop, ff0, tt0, ...
            substation_node, n_bus, loss0); %%%  core algorithm
    end
    t2 = toc;    
    time_consumption.tabu(iter) = t2-t1;

    from_to = show_biograph_not_sorted(Branch(:, [1 2]), substation_node, ...
        show_biograph); %%% show figure, take time
    mpc = generate_mpc(Bus, Branch, n_bus);
    t1 = toc;
    res_pf = runpf(mpc, mpopt);
    t2 = toc;    
    losses = get_losses(res_pf.baseMVA, res_pf.bus, res_pf.branch);
    lossi = sum(real(losses)) % loss = 0.5364    
%     if res_pf.success==0
%         [Vbus, IB, loss_sweep_MW] = SweepPowerFlow_v2(mpc, 1000);
%         lossi = loss_sweep_MW*1000;
%     end
    loss_tabu(iter,1) = lossi;
    yij_dec = generate_yij_from_Branch(Branch, Branch0);

    % record Branch and loss
    Branch_loss_record.tabu(iter,1).Branch = Branch; 
    Branch_loss_record.tabu(iter,1).loss = lossi;
    
    [PQ, PV, REF, NONE, BUS_I, BUS_TYPE, PD, QD, GS, BS, BUS_AREA, VM, ...
      VA, BASE_KV, ZONE, VMAX, VMIN, LAM_P, LAM_Q, MU_VMAX, MU_VMIN] = idx_bus;
%     Vm = res_pf.bus(:, VM)';
%     Va = res_pf.bus(:, VA)';
%     ending_bus = find_ending_node(Branch, substation_node);
%     [ending_bus'; Vm(ending_bus)]; 
    
    %% ---------------------one open and one close---------------------%%
    if lossi < 1.06*loss0_dec % the result of tabu is not too bad        
        % prepare nodes_focused for one_open_one_close
        t1 = toc;
        [nodes_focused] = get_nodes_focused_o1c1( ...
            from_to, Branch, Branch0, substation_node, brch_idx_in_loop, ...
        n1_down_substation, n2_up_ending);

        loss_before_switch0 = lossi;
        if combine3 
            [record_o1c1_loss_dec, loss_after_switch_combine_two_o1c1, Branch_loss] = ...
                one_open_one_close_combine3(nodes_focused, Bus, Branch0, Branch, from_to, ...
                substation_node, n_bus, loss_before_switch0);
        else

🎉3 文献来源

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

[1]J. Zhan, W. Liu, C. Y. Chung and J. Yang, "Switch Opening and Exchange Method for Stochastic Distribution Network Reconfiguration," in IEEE Transactions on Smart Grid, vol. 11, no. 4, pp. 2995-3007, July 2020, doi: 10.1109/TSG.2020.2974922.

🌈4 Matlab代码实现