发明名称 计及分布式发电和储能单元的智能配网重构方法
摘要 本发明公开了一种计及分布式发电和储能单元的智能配网重构方法,首先建立配网重构的多目标多约束数学模型,对模型中的控制变量进行编码,然后采用改进的NSGA‑Ⅱ优化算法进行求解,解算过程中反复调用三相潮流计算模块,将NSGA‑Ⅱ中各代种群的个体进行解码,对解码后的网络重构方案进行三相潮流计算,进而得到各网络重构方案的适应度值,再反馈到NSGA‑Ⅱ进行优化处理。本发明所构建的多目标多约束重构模型能够同时计及分布式电源和储能单元的影响;所采用的多电源三相前推回推算法能够很好地考虑配网的三相不平衡特征,且计算速度快;所采用的带越限剔除策略的快速非支配排序遗传算法能够避免权重设置的盲目性,并能及时排除不可行解,计算效率高。
申请公布号 CN103310065B 申请公布日期 2016.08.10
申请号 CN201310256198.5 申请日期 2013.06.25
申请人 国家电网公司;山东电力集团公司电力科学研究院 发明人 邵志敏;李建修;李立生;张世栋;孙勇
分类号 G06F17/50(2006.01)I;G06Q50/06(2012.01)I 主分类号 G06F17/50(2006.01)I
代理机构 济南圣达知识产权代理有限公司 37221 代理人 张勇
主权项 一种计及分布式发电和储能单元的智能配网重构方法,其特征是,首先建立配网重构的多目标多约束数学模型,对模型中的控制变量进行编码,然后采用改进的NSGA‑Ⅱ优化算法进行求解,解算过程中反复调用三相潮流计算模块,将NSGA‑Ⅱ中各代种群的个体进行解码,对解码后的网络重构方案进行三相潮流计算,进而得到各网络重构方案的适应度值,再反馈到NSGA‑Ⅱ进行优化处理;所述一种分布式发电和储能单元的智能配网重构方法,具体步骤为:步骤一:建立智能配网重构的数学模型,模型中的目标函数有两个,分别为网络损耗最小数学模型和负荷均衡度最高数学模型,表达式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>minf</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><munder><mo>&Sigma;</mo><mrow><mi>i</mi><mo>&Element;</mo><mi>B</mi></mrow></munder><munder><mo>&Sigma;</mo><mrow><mi>m</mi><mo>&Element;</mo><mi>M</mi></mrow></munder><msub><mi>R</mi><mrow><mi>i</mi><mo>,</mo><mi>m</mi></mrow></msub><mfrac><mrow><mo>(</mo><msubsup><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>m</mi></mrow><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>Q</mi><mrow><mi>i</mi><mo>,</mo><mi>m</mi></mrow><mn>2</mn></msubsup><mo>)</mo></mrow><msubsup><mi>V</mi><mrow><mi>i</mi><mo>,</mo><mi>m</mi></mrow><mn>2</mn></msubsup></mfrac><mo>,</mo><mi>x</mi><mo>&Element;</mo><mi>X</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000926511370000011.GIF" wi="1515" he="170" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>minf</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><munder><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mrow><mi>i</mi><mo>&Element;</mo><mi>B</mi><mo>,</mo><mi>m</mi><mo>&Element;</mo><mi>M</mi></mrow></munder><mfrac><mrow><mo>|</mo><msub><mi>I</mi><mrow><mi>i</mi><mo>,</mo><mi>m</mi></mrow></msub><mo>|</mo></mrow><msub><mi>I</mi><mrow><mi>N</mi><mo>,</mo><mi>i</mi><mo>,</mo><mi>m</mi></mrow></msub></mfrac><mo>,</mo><mi>x</mi><mo>&Element;</mo><mi>X</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000926511370000012.GIF" wi="1339" he="167" /></maths>其中,f<sub>1</sub>为网络有功损耗,f<sub>2</sub>为支路最大负荷率;x为控制变量,即配网开关的状态组合,X为所有状态组合所构成的集合;M={A,B,C}为相序集合,B为配网支路集合,P<sub>i,m</sub>、Q<sub>i,m</sub>和I<sub>i,m</sub>分别为流过支路i的有功功率、无功功率和电流;R<sub>i,m</sub>为支路i的支路电阻;V<sub>i,m</sub>为支路i的末端电压,I<sub>N,i,m</sub>为支路i的额定电流;步骤二:确定智能配网重构的数学模型的约束条件;步骤三:对含分布式电源和储能单元的多电源配电网络,利用三相前推回推潮流算法将弱环网转化为辐射型网络,通过在断点两侧注入补偿量的方法来等效智能配网所含的环路的影响;同时根据各类分布式电源电网互联的接口方式以及它们的运行和控制方式,分别建立其在潮流计算中的等值模型,在各迭代步将分布式电源所在节点转换为PQ节点、PI节点或PV节点;步骤四:利用基于Pareto最优解理论的加入越限剔除策略的NSGA‑Ⅱ算法优化智能配网的重构;所述步骤三的具体步骤为:(3‑1)对智能配电网络的原始数据进行初始化,包括配网拓扑信息、元件参数、负荷数据;(3‑2)根据读入的数据计算各负荷节点的注入电流;(3‑3)确定分布式电源和储能单元的等值注入电流;(3‑4)判断网络中是否存在环网回路,若存在环网回路则向解环节点叠加注入电流;(3‑5)叠加储能单元及分布式电源的三相注入电流,进行三相前推回推潮流计算;(3‑6)从馈线末端节点开始进行前推计算,对支路电流进行求和,从而得到各条支路始端节点的三相电流;(3‑7)从馈线首端节点开始回推计算各支路末端节点的电压,同时对各节点的三相电压进行修正;(3‑8)修正PV节点注入的无功功率,然后判断PV节点的无功功率是否发生越界,若发生越界则将PV节点转化成PQ节点转入步骤(3‑2)重新进行计算,否则转入(3‑9);(3‑9)以节点相邻两次迭代电压幅值差是否满足预设精度为收敛条件,若满足条件则迭代收敛,转入(3‑10),否则重新修正所在节点的注入电流,然后转入步骤(3‑2)重新计算;(3‑10)计算结束,输出潮流计算结果;所述步骤二中的约束条件为:(2‑1)功率平衡约束,对智能配网任一支路i,需满足如下条件:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>S</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>=</mo><msub><mi>S</mi><mrow><mi>L</mi><mi>i</mi></mrow></msub><mo>+</mo><msub><mi>S</mi><mrow><mi>S</mi><mi>i</mi></mrow></msub><mo>-</mo><msub><mi>S</mi><mrow><mi>G</mi><mi>i</mi></mrow></msub><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msub><mi>S</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>,</mo><mi>i</mi><mo>&Element;</mo><mi>B</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000926511370000021.GIF" wi="1266" he="154" /></maths>S<sub>i</sub>=S′<sub>i</sub>+ΔS<sub>i</sub>,i∈B             (4)其中,S<sub>i</sub>和S′<sub>i</sub>分别为支路i首段和末端的复功率,S<sub>Li</sub>为节点i的负荷功率,S<sub>Gi</sub>和S<sub>Si</sub>分别为节点i所接DG的注入功率和SU的充放电功率,S<sub>i,j</sub>分别为节点i的第j条出支的首端复功率,ΔS<sub>i</sub>为支路i消耗的复功率,B为节点集合;(2‑2)支路电流约束,即各支路电流幅值应不大于最大允许载流量:|I<sub>i,m</sub>|≤I<sub>i,m,max</sub>,i∈B           (5)其中,I<sub>i,m,max</sub>为支路i第m相的最大允许电流;(2‑3)节点电压约束,即各节点的电压V<sub>i</sub>应满足:V<sub>i,min</sub>≤V<sub>i</sub>≤V<sub>i,max</sub>,i∈N      (6)其中,V<sub>i,min</sub>和V<sub>i,max</sub>分别为节点i的电压下限和上限,N为所有节点的集合;(2‑4)变压器容量约束:|S<sub>t</sub>|≤|S<sub>t,max</sub>|,t∈T     (7)其中,S<sub>t</sub>和S<sub>t,max</sub>分别为变压器t的实际功率和额定容量,T为所有变压器的集合;(2‑5)储能单元状态约束,包括功率约束、电压约束、注入电流约束;分别为:P<sub>gi,min</sub>≤P<sub>g,i</sub>≤P<sub>gi,max</sub>,i∈G<sub>P</sub>                   (8)Q<sub>gi,min</sub>≤Q<sub>g,i</sub>≤Q<sub>gi,max</sub>,i∈G<sub>Q</sub>          (9)V<sub>gi,min</sub>≤V<sub>g,i</sub>≤V<sub>gi,max</sub>,i∈G<sub>V</sub>         (10)I<sub>gi,min</sub>≤I<sub>g,i</sub>≤I<sub>gi,max</sub>,i∈G<sub>I</sub>         (11)其中,P<sub>g,i</sub>和Q<sub>g,i</sub>分别为第i个DG的有功和无功输出,P<sub>gi,max</sub>和P<sub>gi,min</sub>分别为其对应的有功功率上、下限,Q<sub>gi,max</sub>和Q<sub>gi,min</sub>分别为其对应的无功功率上、下限,V<sub>g,i</sub>为第i个DG的端电压,V<sub>gi,max</sub>和V<sub>gi,min</sub>分别为其对应的电压上、下限,I<sub>g,i</sub>为第i个DG的注入电流,I<sub>gi,max</sub>和I<sub>gi,min</sub>分别为其对应的输出电流上、下限,G<sub>P</sub>、G<sub>Q</sub>、G<sub>V</sub>和G<sub>I</sub>分别为有功可调、无功可调、电压可调和电流可调的可控DG集合;分布式电源充放电状态约束,储能单元i的充放电功率S<sub>Si</sub>应满足如下约束:S<sub>DN,i</sub>≤S<sub>Si</sub>≤S<sub>CN,i</sub>               (12)其中,S<sub>CN,i</sub>为储能单元i的额定充电功率,S<sub>DN,i</sub>为储能单元i的额定放电功率,其中,额定放电功率为负值;网络拓扑约束,重构后的网络g<sub>k</sub>应满足:g<sub>k</sub>∈G<sub>R</sub>         (13)其中,G<sub>R</sub>为所有可能的辐射状网络构成的集合;所述步骤四的具体步骤为:(4‑1)数据初始化,读入智能配电网络的基础参数和运行数据,同时设置遗传算法参数;设种群规模为W,最大进化代数为r<sub>max</sub>,收敛标准差σ<sub>a</sub>;(4‑2)置进化代数r为0,然后随机生成初始种群,采用二进制方式对配网重构模型进行编码,进化代数指种群整体进化的次数,进化前为第0代群体,每进化一次r增加1;(4‑3)对种群进行越限检查,将种群中的个体解码后,通过三相前推回推潮流计算,得到系统潮流分布状态,按智能配电网络重构数学模型的约束条件检查有无个体越限,若有则剔除越限个体,同时随机生成新个体进行补齐,直到不再有越限个体为止;(4‑4)按照智能配网重构的数学模型的两个目标函数确定函数值f<sub>1</sub>和f<sub>2</sub>,并进行快速Pareto非支配排序,确定个体所属的非支配层级,对处于同一层的个体,则通过拥挤距离的计算区分优劣;(4‑5)令r=r+1,通过选择算子、交叉算子、和变异算子的处理,得到下一代个体,重新进行越限剔除和个体评价;(4‑6)从第二代个体开始,经过越限剔除操作后,采用精英保留策略,将子代种群Cr和父代种群Pr进行合并,形成过渡种群,两代个体共同竞争,以个体所处的非支配层级和拥挤度进行排序,取前W个非支配个体作为下一代的父种群;(4‑7)满足如下收敛条件时,算法停止并输出最优解集;收敛条件:r≥r<sub>max</sub>或σ<sub>r</sub>≤σ<sub>a</sub>其中,σ<sub>r</sub>为第r代种群的适应度标准差。
地址 100031 北京市西城区西长安街86号