发明名称 一种基于智能优化技术的含新能源配电网状态估计方法
摘要 本发明公开了一种基于智能优化技术的含新能源配电网状态估计方法,对非线性、离散配电网的运行状态进行在线实时的估计,根据输入的REDG有功功率和负载功率的量测值实时的估计出REDGs的有功功率和负载功率的平均输出值。本发明方法将智能优化算法中量子进化算法QIEA全局寻优能力强、收敛速度快和贪婪随机自适应算法GRASP局部搜索能力强的优点相结合,给出一种混合量子进化算法(mQIEA),能够根据输入的REDG有功功率和负载功率的量测值实时的估计出REDGs的有功功率和负载功率的平均输出值,便于系统有效地监控系统运行状态。
申请公布号 CN103208797A 申请公布日期 2013.07.17
申请号 CN201310085442.6 申请日期 2013.03.18
申请人 西南交通大学 发明人 张葛祥;赵俊博
分类号 H02J3/00(2006.01)I;G06N3/12(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 成都信博专利代理有限责任公司 51200 代理人 张澎
主权项 1.一种基于智能优化技术的含新能源配电网状态估计方法,对非线性、离散配电网的运行状态进行在线实时的估计,根据输入的REDG有功功率和负载功率的量测值实时的估计出REDGs的有功功率和负载功率的平均输出值,具体作法包括:步骤1:读取系统的网络数据,包括网络结构和线路阻抗、量测值、REDGs有功功率和负载功率平均值以及它们的标准差;通过量测数据和负载率来初始化每个负载值;通过每个REDG的输出功率平均值来初始化每个REDG的有功功率;步骤2:通过前推回代法计算配电网的潮流,并获得每个变量的上下限约束;步骤3:产生初始种群Q(t),<img file="FDA00002929735200011.GIF" wi="390" he="79" />(t=1),n为种群规模,t为进化代数,<img file="FDA00002929735200012.GIF" wi="47" he="65" />(i=1,2,...,n)为种群中第t代的第i个个体,其形式如下所示<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>q</mi><mi>i</mi><mi>t</mi></msubsup><mo>=</mo><mfenced open='[' close='|'><mtable><mtr><mtd><msubsup><mi>&alpha;</mi><mrow><mi>i</mi><mn>1</mn></mrow><mi>t</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>&beta;</mi><mrow><mi>i</mi><mn>1</mn></mrow><mi>t</mi></msubsup></mtd></mtr></mtable></mfenced><mfenced open='' close='|'><mtable><mtr><mtd><msubsup><mi>&alpha;</mi><mrow><mi>i</mi><mn>2</mn></mrow><mi>t</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>&beta;</mi><mrow><mi>i</mi><mn>2</mn></mrow><mi>t</mi></msubsup></mtd></mtr></mtable></mfenced><mfenced open='' close=''><mtable><mtr><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd></mtr></mtable></mfenced><mfenced open='|' close=']'><mtable><mtr><mtd><msubsup><mi>&alpha;</mi><mi>im</mi><mi>t</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>&beta;</mi><mi>im</mi><mi>t</mi></msubsup></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>式中,<img file="FDA00002929735200014.GIF" wi="174" he="85" />(j=1,2,...,m)是一个量子比特,并且满足归一化条件<img file="FDA00002929735200015.GIF" wi="279" he="86" />m为量子比特的个数,即一个量子染色体的长度,在初始化时,将全部n条染色体的2mn个概率幅,即所有的表示状态|0&gt;概率幅的<img file="FDA00002929735200016.GIF" wi="53" he="77" />和表示状态|1&gt;概率幅的<img file="FDA00002929735200017.GIF" wi="55" he="74" />(i=1,2,...,n;j=1,2,...,m)均置为<img file="FDA00002929735200018.GIF" wi="128" he="98" />即<maths num="0002"><![CDATA[<math><mrow><mi>Q</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><mn>1</mn><mo>/</mo><msqrt><mn>2</mn></msqrt></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><mn>1</mn><mo>/</mo><msqrt><mn>2</mn></msqrt></mtd></mtr><mtr><mtd><mi>M</mi></mtd><mtd><mi>O</mi></mtd><mtd><mi>M</mi></mtd></mtr><mtr><mtd><mn>1</mn><mo>/</mo><msqrt><mn>2</mn></msqrt></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><mn>1</mn><mo>/</mo><msqrt><mn>2</mn></msqrt></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>它表示在t=1时,每条染色体以相同的概率处于所有可能状态的线性叠加m=hl,h是REDGs和负载的总个数,l是每个REDG的有功功率和负载功率量子比特的长度;步骤4:将有约束的配电网状态估计问题转化为无约束的配电网状态估计问题;在此步骤中,采用罚函数将配电网状态估计问题中任何违反约束条件的值纳入惩罚性因子,分别将外点法和内点法用于处理等式和不等式约束,从而配电网状态估计的目标函数转换成一个无约束增广目标函数:<maths num="0003"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mi>F</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mn>1</mn><msqrt><msup><mi>r</mi><mi>t</mi></msup></msqrt></mfrac><msup><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mn>1</mn></msub></munderover><msub><mi>J</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mi>r</mi><mi>t</mi></msup><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mn>2</mn></msub></munderover><mfrac><mn>1</mn><mrow><msub><mi>g</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></mfrac></mtd></mtr><mtr><mtd><msub><mi>J</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2,3</mn><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>N</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>g</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2,3</mn><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>N</mi><mn>2</mn></msub></mtd></mtr></mtable></mfenced></math>]]></maths>其中,f(x)是配电网状态估计的目标函数值,N<sub>1</sub>和N<sub>2</sub>分别为等式和不等式约束的数量,J<sub>j</sub>(x)和g<sub>j</sub>(x)分别为等式和不等式约束;通常惩罚函数的初始值都被设置为1,本发明中罚函数为<maths num="0004"><![CDATA[<math><mrow><msup><mi>r</mi><mi>t</mi></msup><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>)</mo></mrow><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msup><mo>;</mo></mrow></math>]]></maths>步骤5:对种群Q(t)进行观测,产生观测态值<img file="FDA00002929735200022.GIF" wi="410" he="82" />其中<img file="FDA00002929735200023.GIF" wi="47" he="61" />(i=1,2,...n)是由<img file="FDA00002929735200024.GIF" wi="43" he="63" />(i=1,2,...n)观测产生的一个长度为m的二进制串;观测过程如下:如果random[0,1)&lt;|α|<sup>2</sup>,观测值就为0,否则为1,其中,α为任一量子比特的概率幅;步骤6:通过步骤4中的目标函数来评价P(t);步骤7:将P(t)中的最优个体存储到最优个体的集合s(t);步骤8:采用贪婪随机自适应算法GRASP对s(t)进行局部搜索,产生新的种群R(t);步骤9:通过步骤4中的目标函数来评价局部搜索所产生的个体集合R(t);步骤10:将R(t)中的最优个体存储到s(t);步骤11:量子门更新;步骤12:终止条件判断,若满足终止条件,则向智能电网能量管理系统或配电网实时状态监控系统输出REDG有功功率和负载功率平均值的估计值,若不满足,则令t←t+1,<maths num="0005"><![CDATA[<math><mrow><msup><mi>r</mi><mi>t</mi></msup><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>)</mo></mrow><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msup><mo>,</mo></mrow></math>]]></maths>并转步骤5;
地址 610031 四川省成都市二环路北一段111号西南交通大学科技处