发明名称 基于改进人工鱼群混合优化算法的大规模风电并网无功电压优化方法
摘要 本发明公开了基于改进人工鱼群混合优化算法的大规模风电并网无功电压优化方法,包括:调度选定大规模风电基地的自动化管控系统,基于自动化管控系统中的能量管理系统,读取电网的原始数据信息;基于所得原始数据信息,获取用于潮流计算的有效数据,形成全网的有效数据集;基于预先确定的优化方式,对所得有效数据集进行筛选,获取基于全网数据的计算数集;调用改进人工鱼群混合优化算法,基于筛选所得计算数集,对风电并网无功电压进行无功优化处理,获取无功优化处理结果。本发明所述方法,可以克服现有技术中数据处理量大、收敛性差与稳定性不好等缺陷,以实现数据处理量小、收敛性好与稳定性好的优点。
申请公布号 CN102684207B 申请公布日期 2014.07.09
申请号 CN201210162209.9 申请日期 2012.05.23
申请人 甘肃省电力公司电力科学研究院;华北电力大学;国家电网公司 发明人 郑伟;刘文颖;王维洲;谢昶;韩旭杉;徐鹏;马超;李波;杨勇;金娜;温志伟;杜珣;安亮亮;李亚龙;卢甜甜;梁才;周海洋;门德月;但扬清;王久成;李扬;刘景延;王佳明;杨斌;曹俊龙;梁纪峰;杜波;吴耀昊;钟佳辰;文晶;葛润东;吴晓丹;杨楠;刘茜;邢晶;赵子兰
分类号 H02J3/18(2006.01)I;H02J3/38(2006.01)I;H02J13/00(2006.01)I 主分类号 H02J3/18(2006.01)I
代理机构 北京中恒高博知识产权代理有限公司 11249 代理人 宋敏
主权项 基于改进人工鱼群混合优化算法的大规模风电并网无功电压优化方法,其特征在于,包括:a、调度运行人员通过调度主站的能量管理系统(EMS),读取大规模风电的原始数据信息;所述原始数据信息,至少包括用作电网静态分析数据的电网静态参数与约束条件、以及用作反映电网实时运行状态计算数据的电网实时数据;b、基于步骤a读取的原始数据信息,获取用于潮流计算的有效数据,并形成全网的有效数据集;c、基于预先确定的优化方式,对步骤b所得有效数据集进行筛选,获取基于全网数据的计算数集;d、调用改进人工鱼群混合优化算法,基于步骤c筛选所得计算数集,对风电并网无功电压进行无功优化处理,获取无功优化处理结果;在步骤c中,所述预先确定的优化方式,至少包括常规风电并网的无功优化方式与大规模风电并网的无功优化方式;在步骤d中,所述对风电并网无功电压进行无功优化处理的操作,具体包括:建立至少包含风电间歇性和随机性因素的无功电压优化控制模型;基于该无功电压优化控制模型,利用遗传算法在原始数据信息的初始数据中形成用于进行全网潮流计算的有效数据集后,在所述有效数据集的待优化区域中,设立用于记录最优人工鱼个体状态的公告板;每条人工鱼在行动一次后,将表示自身当前状态的函数值与公告板进行比较,当表示自身当前状态的函数值优于公告板时,则用自身当前状态取代公告板,并将所述待优化区域的公告板从整个电网中划分出来,形成对风电并网无功电压无功优化的最终计算数据;主要叙述建立包含风电间歇性和随机性因素的无功电压优化控制模型的操作,具体包括:建立改进人工鱼群混合优化算法的目标函数:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mi>min</mi><mrow><mo>(</mo><msub><mi>P</mi><mi>LOSS</mi></msub><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><mfrac><msub><mi>&Delta;V</mi><mi>i</mi></msub><msub><mi>V</mi><mi>M</mi></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><mfrac><msub><mi>&Delta;Q</mi><mi>i</mi></msub><msub><mi>Q</mi><mi>M</mi></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000454625280000011.GIF" wi="1469" he="246" /></maths>上述目标函数的等式约束条件为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>P</mi><mi>i</mi></msub><mo>=</mo><msub><mi>V</mi><mi>i</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>V</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>G</mi><mi>ij</mi></msub><mi>cos</mi><msub><mi>&delta;</mi><mi>ij</mi></msub><mo>+</mo><msub><mi>B</mi><mi>ij</mi></msub><mi>sin</mi><msub><mi>&delta;</mi><mi>ij</mi></msub><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000454625280000021.GIF" wi="1538" he="326" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>i</mi></msub><mo>=</mo><msub><mi>V</mi><mi>i</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>V</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>G</mi><mi>ij</mi></msub><mi>sin</mi><msub><mi>&delta;</mi><mi>ij</mi></msub><mo>-</mo><msub><mi>B</mi><mi>ij</mi></msub><mi>cos</mi><msub><mi>&delta;</mi><mi>ij</mi></msub><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000454625280000022.GIF" wi="1550" he="326" /></maths>上述目标函数的不等式约束中控制变量约束条件为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>T</mi><mrow><mi>i</mi><mi>min</mi></mrow></msub><mo>&le;</mo><msub><mi>T</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>T</mi><mrow><mi>i</mi><mi>max</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>Q</mi><mrow><mi>ci</mi><mi>min</mi></mrow></msub><mo>&le;</mo><msub><mi>Q</mi><mi>ci</mi></msub><mo>&le;</mo><msub><mi>Q</mi><mrow><mi>ci</mi><mi>max</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>V</mi><mi>Gimi</mi></msub><mo>&le;</mo><msub><mi>V</mi><mi>Gi</mi></msub><mo>&le;</mo><msub><mi>V</mi><mrow><mi>Gi</mi><mi>max</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>V</mi><mrow><mi>i</mi><mi>min</mi></mrow></msub><mo>&le;</mo><msub><mi>V</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>V</mi><mi>imaz</mi></msub></mtd></mtr><mtr><mtd><msub><mi>Q</mi><mrow><mi>Gi</mi><mi>min</mi></mrow></msub><mo>&le;</mo><msub><mi>Q</mi><mi>Gi</mi></msub><mo>&le;</mo><msub><mi>Q</mi><mrow><mi>Gi</mi><mi>max</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000454625280000023.GIF" wi="1041" he="829" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>P</mi><mi>LOSS</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>V</mi><mi>i</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>V</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>G</mi><mi>ij</mi></msub><mi>cos</mi><msub><mi>&delta;</mi><mi>ij</mi></msub><mo>+</mo><msub><mi>B</mi><mi>ij</mi></msub><mi>sin</mi><msub><mi>&delta;</mi><mi>ij</mi></msub><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000454625280000024.GIF" wi="1537" he="269" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>&Delta;V</mi><mi>i</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>V</mi><mi>i</mi></msub><mo>-</mo><msub><mi>V</mi><mrow><mi>i</mi><mi>max</mi></mrow></msub></mtd><mtd><mrow><mo>(</mo><msub><mi>V</mi><mi>i</mi></msub><mo>></mo><msub><mi>V</mi><mrow><mi>i</mi><mi>max</mi></mrow></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>(</mo><msub><mi>V</mi><mrow><mi>i</mi><mi>min</mi></mrow></msub><mo>&le;</mo><msub><mi>V</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>V</mi><mrow><mi>i</mi><mi>max</mi></mrow></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>V</mi><mrow><mi>i</mi><mi>min</mi></mrow></msub><mo>-</mo><msub><mi>V</mi><mi>i</mi></msub></mtd><mtd><mrow><mo>(</mo><msub><mi>V</mi><mi>i</mi></msub><mo>></mo><msub><mi>V</mi><mrow><mi>i</mi><mi>min</mi></mrow></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000454625280000025.GIF" wi="1421" he="393" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>&Delta;Q</mi><mi>i</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>Q</mi><mi>i</mi></msub><mo>-</mo><msub><mi>Q</mi><mrow><mi>i</mi><mi>max</mi></mrow></msub></mtd><mtd><mrow><mo>(</mo><msub><mi>Q</mi><mi>i</mi></msub><mo>></mo><msub><mi>Q</mi><mrow><mi>i</mi><mi>max</mi></mrow></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>(</mo><msub><mi>Q</mi><mrow><mi>i</mi><mi>min</mi></mrow></msub><mo>&le;</mo><msub><mi>Q</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>Q</mi><mrow><mi>i</mi><mi>max</mi></mrow></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>Q</mi><mrow><mi>i</mi><mi>min</mi></mrow></msub><mo>-</mo><msub><mi>Q</mi><mi>i</mi></msub></mtd><mtd><mrow><mo>(</mo><msub><mi>Q</mi><mi>i</mi></msub><mo>></mo><msub><mi>Q</mi><mrow><mi>i</mi><mi>min</mi></mrow></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000454625280000026.GIF" wi="1657" he="363" /></maths>V<sub>M</sub>=V<sub>imax</sub>‑V<sub>imin</sub>;Q<sub>M</sub>=Q<sub>imax</sub>‑Q<sub>imin</sub>;其中,V<sub>i</sub>为线路首端节点i的电压幅值,V<sub>j</sub>为线路末端节点j的电压幅值,δ<sub>ij</sub>为线路首末端节点i、j的电压相角差,T<sub>i</sub>为可调变压器分接头的位置,Q<sub>ci</sub>为容性无功补偿容量,V<sub>Gi</sub>为发电机机端电压,Q<sub>Gi</sub>为发电机的无功出力,P<sub>LOSS</sub>电网的有功损耗,G<sub>ij</sub>是线路(i,j)的电导,B<sub>ij</sub>是线路(i,j)的电纳,N是系统的节点数,P<sub>i</sub>是节点注入的有功功率,Q<sub>i</sub>是节点注入的无功功率,T<sub>imax</sub>是变压器分接头的最高档位,T<sub>imin</sub>是变压器分接头最低档位,Q<sub>cimax</sub>和Q<sub>cimin</sub>分别是节点无功补偿量的上限和下限,V<sub>Gimin</sub>和V<sub>Gimax</sub>分别是发电机机端电压的上和下限,V<sub>imin</sub>和V<sub>imax</sub>分别是节点电压的上限和下限,Q<sub>Gimin</sub>和Q<sub>Gimax</sub>分别是发电机无功出力的上限和下限。
地址 730050 甘肃省兰州市七里河区西津东路648号