发明名称 基于佳点集量子粒子群算法的微网优化方法
摘要 一种基于佳点集量子粒子群算法的微网优化方法,包括:建立微网多目标优化模型,包括微网多目标优化的目标函数公式和约束条件公式;使用佳点集改进量子粒子群算法;使用佳点集量子粒子群算法对微网多目标优化模型求解。本发明根据最大负荷确定微型燃气轮机的容量,再优化间歇式分布式电源以及储能系统的容量。求解过程中采用佳点集量子粒子群算法,保证了优化结果同时具有良好的选择性和指导性。内置的基于滤波的调度策略能充分发挥微型燃气轮机和储能系统的特性,不但拥有频谱分析法的优势,还具有更好的实时性,能够计算得到微网中分布式电源、储能和微型燃气轮机的位置与容量,实现孤立微网中经济性和污染物排放的整体优化。
申请公布号 CN105958482A 申请公布日期 2016.09.21
申请号 CN201610382257.7 申请日期 2016.05.31
申请人 天津天大求实电力新技术股份有限公司;国网天津市电力公司;国家电网公司 发明人 申刚;张岩;尚德华;杨毅;张源超;庄剑;于建成;项添春;王旭东;丁一;戚艳
分类号 H02J3/00(2006.01)I;H02J3/38(2006.01)I;G06F17/50(2006.01)I;G06Q10/06(2012.01)I;G06Q50/06(2012.01)I;G06N3/00(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 天津市北洋有限责任专利代理事务所 12201 代理人 杜文茹
主权项 一种基于佳点集量子粒子群算法的微网优化方法,其特征在于,如下步骤:1)建立微网多目标优化模型,包括微网多目标优化的目标函数公式和约束条件公式,其中,所述的微网多目标优化的目标函数公式包括:(1)经济性目标公式:min C<sub>s</sub>=∑C<sub>i</sub>                           (1)<img file="FDA0001005613380000011.GIF" wi="1078" he="126" />C<sub>i,j</sub>=C<sub>i,j,Cap</sub>+C<sub>i,j,Rep</sub>+C<sub>i,j,OM</sub>+C<sub>i,j,Sal</sub>             (3)式中,C<sub>s</sub>为微网系统全寿命周期n年的净现值,C<sub>i</sub>为微网第i个元件的净现值,C<sub>i,j</sub>为第i个元件在第j年的现金流,所述的元件包括风力发电机、光伏、储能和微型燃气轮机;C<sub>i,j,Cap</sub>为初建成本,C<sub>i,j,Rep</sub>为更新成本,C<sub>i,j,OM</sub>为运行维护成本,C<sub>i,j,Sal</sub>为残值;(2)污染物排放目标公式:<img file="FDA0001005613380000012.GIF" wi="1078" he="79" />式中,V<sub>CO2</sub>为年CO<sub>2</sub>排放体积,V<sub>Gas</sub>为年天然气耗量,σ为天然气排放系数,即单位体积的天然气燃烧产生的二氧化碳量;所述的微网多目标优化的约束条件公式,包括:(3)孤立微网运行的等式约束有:电功率平衡约束公式,P<sub>Load</sub>=∑P<sub>PV,i</sub>+∑P<sub>WT,i</sub>+∑P<sub>MT,i</sub>+P<sub>ES</sub>                  (10)和储能充放电量平衡约束公式E<sub>charge</sub>=E<sub>discharge</sub>                          (11)储能充放电量平衡约束公式是指在一个评估时间段内,储能充放电量要保证相等以便储能的循环利用;式中P<sub>Load</sub>、∑P<sub>PV,i</sub>、∑P<sub>WT,i</sub>、∑P<sub>MT,i</sub>、P<sub>ES</sub>分别对应为微网中总负荷功率、总光伏功率、总风机功率、总微型燃气轮机功率以及储能功率;(4)不等式约束公式,包括有:分布式电源约束公式、储能设备的配置数目约束公式、储能的寿命约束公式和储能的荷电状态约束公式,其中,所述分布式电源约束公式:0<N<sub>DG</sub><N<sub>DG,max</sub>                            (12)所述储能设备的配置数目约束公式:0<N<sub>ES</sub><N<sub>ES,max</sub>                              (13)所述的储能寿命约束公式:L<sub>Bat</sub>>L<sub>Limit</sub>                                  (14)所述的储能的荷电状态约束公式:S<sub>min</sub><S<S<sub>max</sub>                              (15)式中,N<sub>DG,max</sub>和N<sub>ES,max</sub>为对给定容量的分布式电源和储能所预设的最大配置数目,N<sub>DG</sub>和N<sub>ES</sub>分别为对给定容量的分布式电源和储能所预设的实际数目,L<sub>Limit</sub>和L<sub>Bat</sub>分别为预期的储能运行寿命实际值和最小值,S为储能的荷电状态,S<sub>min</sub>和S<sub>max</sub>分别为储能的荷电状态的最小和最大值。2)使用佳点集改进量子粒子群算法,包括:(1)更新量子位幅角增量更新后的量子位幅角增量公式如下:θ<sub>ij</sub>(t+1)=θ<sub>ij</sub>(t)+△θ<sub>ij</sub>(t+1)                (18)△θ<sub>ij</sub>(t+1)=w△θ<sub>ij</sub>(t)+c<sub>1</sub>r<sub>1</sub>(△θ<sub>1</sub>)+c<sub>2</sub>r<sub>2</sub>(△θ<sub>g</sub>)        (19)其中:c<sub>1</sub>和c<sub>2</sub>分别为个体和全局的学习因子,r<sub>1</sub>和r<sub>2</sub>为区间[0,1]内的随机数,w为混沌时间序列数映射到[0.1,0.9]区间上的数值,θ<sub>ij</sub>为量子比特的相位,Δθ<sub>1</sub>为当前个体与个体之间的角度差,Δθ<sub>g</sub>为当前个体与全局最优之间的角度差,Δθ<sub>1</sub>与Δθ<sub>g</sub>的公式如下<img file="FDA0001005613380000021.GIF" wi="1182" he="246" /><img file="FDA0001005613380000022.GIF" wi="1182" he="254" />式中,θ<sub>ilj</sub>为当前代数最优量子比特的相位,θ<sub>gj</sub>为全局最优量子比特的相位。(2)佳点集交叉操作选取两个粒子θ<sub>i</sub>和θ<sub>j</sub>作为父代,设θ<sub>i</sub>=θ<sub>i1</sub>,θ<sub>i2</sub>,…,θ<sub>is</sub>,θ<sub>j</sub>=θ<sub>j1</sub>,θ<sub>j2</sub>,…,θ<sub>js</sub>,θ<sub>i</sub>和θ<sub>j</sub>共同确定了一个有界闭区间<img file="FDA0001005613380000023.GIF" wi="251" he="61" /><img file="FDA0001005613380000024.GIF" wi="918" he="151" />a<sub>k</sub>、b<sub>k</sub>分别为父代两个粒子最小和最大值组合,s为父代粒子个数D是R上的超长方体,即<img file="FDA0001005613380000025.GIF" wi="886" he="118" />式中,U([a,b])为子代个体的距离空间,父代个体通过佳点集交叉产生子代个体,子 代个体通过切割父代个体确定的超长方体获得新的基因片段,然后重组产生;(3)变异操作使用量子非门实现变异操作,过程如下:<img file="FDA0001005613380000031.GIF" wi="901" he="143" />令变异概率为p<sub>m</sub>,每个粒子在0~1之间设定一个随机数N<sub>mdi</sub>,若N<sub>mdi</sub>&lt;p<sub>m</sub>,则用量子非门兑换两个概率幅,该粒子的自身最优位置和转向角仍保持不变,N<sub>mdi</sub>≥p<sub>m</sub>,则不进行变异操作;3)使用佳点集量子粒子群算法对微网多目标优化模型求解,包括:(1)初始化佳点集量子粒子群算法的有关参数,包括种群规模、变量个数、迭代次数、解空间范围;(2)对量子位幅角进行初始化,生成光伏容量、风力发电机容量、储能容量和微型燃气轮机容量的值;(3)将每个粒子带入微网多目标优化模型的目标函数,计算得到微网全寿命周期成本和污染物排放量状态变量值;(4)对每个粒子的初始位置进行评价,计算出每个粒子位置的适应值,若粒子目前的位置优于自身记忆的最优位置,则用目前位置替换;若目前全局最优位置优于到目前为止所搜索到的最优位置,则用全局最优位置替换;(5)根据步骤2)第(2)步中更新后的量子位幅角增量公式更新粒子位置;(6)对粒子进行交叉变异操作;(7)判断是否达到最大代数,如果达到,则计算结束,否则将第(6)步进行交叉变异操作得到的粒子重新进行第(3)~第(7)步的操作。
地址 300384 天津市南开区华苑产业园区海泰发展六道6号海泰绿色产业基地J-408