发明名称 一种配电网无功优化方法
摘要 本发明公开一种配电网无功优化方法,包括以下几个步骤:(1)建立风机的数学模型;(2)初始化节点电压;(3)建立无功规划数学模型,转化成多目标优化问题;(4)输入原始数据;(5)形成初始粒子群,初始化粒子速度和位置;(6)潮流计算,得到粒子适应值和当前最优值;(7)用小生境权重飞行时间法修正粒子速度和位置;(8)计算整个种群适应度,得到适应度值,并更新当前最优解pBest<sub>id</sub>;(9)对个体粒子中越界点的控制变量进行调整,修正越界状态变量;;(10)判断是否满足终止条件,若满足则输出结果,否则返回(6)。
申请公布号 CN102856917B 申请公布日期 2014.11.19
申请号 CN201210269539.8 申请日期 2012.07.31
申请人 上海交通大学;吉林省电力有限公司延边供电公司;国家电网公司 发明人 王昕;李立学;郑益慧;王希;张忠保;邵凤鹏;柳杨
分类号 H02J3/18(2006.01)I 主分类号 H02J3/18(2006.01)I
代理机构 上海汉声知识产权代理有限公司 31236 代理人 胡晶
主权项 一种配电网无功优化方法,其基于小生境权重飞行时间粒子群优化算法,其特征在于,包括以下步骤:(1)建立风机的数学模型,给定风电场风速,并在此风速下确定风机输出的功率,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>P</mi><mi>wfarm</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mi>V</mi><mo>></mo><msub><mi>V</mi><mi>co</mi></msub><mi>orV</mi><mo>&lt;</mo><msub><mi>V</mi><mi>ci</mi></msub></mtd></mtr><mtr><mtd><mfrac><mrow><msub><mi>P</mi><mi>r</mi></msub><mrow><mo>(</mo><msubsup><mi>V</mi><mi>w</mi><mn>3</mn></msubsup><mo>-</mo><msubsup><mi>V</mi><mi>ci</mi><mn>3</mn></msubsup><mo>)</mo></mrow></mrow><mrow><mo>(</mo><msubsup><mi>V</mi><mi>r</mi><mn>3</mn></msubsup><mo>-</mo><msubsup><mi>V</mi><mi>ci</mi><mn>3</mn></msubsup><mo>)</mo></mrow></mfrac></mtd><mtd><msub><mi>V</mi><mi>ci</mi></msub><mo>&le;</mo><mi>V</mi><mo>&le;</mo><msub><mi>V</mi><mi>r</mi></msub></mtd></mtr><mtr><mtd><msub><mi>P</mi><mi>r</mi></msub></mtd><mtd><msub><mi>V</mi><mi>r</mi></msub><mo>&le;</mo><mi>V</mi><mo>&le;</mo><msub><mi>V</mi><mi>co</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000537337270000011.GIF" wi="1172" he="439" /></maths>其中,V为风速,V<sub>co</sub>为切除风速;V<sub>ci</sub>为切入风速;V<sub>r</sub>为风机组额定风速;P<sub>r</sub>为风机的额定功率;V<sub>w</sub>为风机轮毂高度的风速;(2)在含风电场的电力系统中形成节点导纳矩阵,设定各节点电压初值;(3)建立无功规划的数学模型,将此数学模型看作是一个多目标优化问题,其中包括目标函数、功率约束方程;(4)获取系统节点信息和支路信息,获取控制变量的个数及各自的取值范围;(5)将无功功率的数学模型的解看成一个粒子,形成初始种群,获取粒子群的群体规模等参数,设置最大迭代次数,随机产生全部粒子及其自身的初始位置和初始速度;(6)对粒子群中的每个粒子带入进行潮流计算,计算群中每个粒子的网损,得到各粒子的适应度值和当前最优解pBest<sub>id</sub>;(7)运用小生境权重飞行时间粒子群优化算法,调整粒子的速度<img file="FDA0000537337270000012.GIF" wi="93" he="83" />和位置<img file="FDA0000537337270000013.GIF" wi="124" he="83" />计算得到新的适应度值,并更新最优解pBest<sub>ij</sub>;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>V</mi><mi>ij</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mi>w</mi><msubsup><mi>V</mi><mi>ij</mi><mi>t</mi></msubsup><mo>+</mo><msub><mi>c</mi><mn>1</mn></msub><msub><mi>r</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>pBest</mi><mi>ij</mi></msub><mo>-</mo><msubsup><mi>X</mi><mi>IJ</mi><mi>T</mi></msubsup><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000537337270000014.GIF" wi="907" he="95" /></maths>w=w<sub>0</sub>‑Δw·t/t<sub>max</sub>;<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>T</mi><mi>t</mi></msub><mo>=</mo><msub><mi>T</mi><mn>0</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mfrac><mi>t</mi><msub><mi>t</mi><mi>max</mi></msub></mfrac><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000537337270000021.GIF" wi="495" he="169" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>X</mi><mi>ij</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mi>X</mi><mi>ij</mi><mi>t</mi></msubsup><mo>+</mo><msubsup><mi>V</mi><mi>ij</mi><mi>t</mi></msubsup><msub><mi>T</mi><mi>t</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA0000537337270000022.GIF" wi="550" he="101" /></maths>式中:t为当前迭代次数,c<sub>1</sub>为学习因子;r<sub>1</sub>是介于[0,1]的随机数;V<sub>ij</sub>∈[‑V<sub>max</sub>,V<sub>max</sub>],V<sub>max</sub>为粒子进化过程中得最大速度;w<sub>0</sub>为权重系数的初始值;Δw权重系数变化率;t<sub>max</sub>为最大迭代次数;T<sub>t</sub>表示第t代粒子的飞行时间;T<sub>0</sub>表示粒子最长飞行时间;t<sub>0</sub>为比例系数,起调节作用;(8)对个体粒子中越界点的控制变量进行调整,修正越界状态变量;(9)若粒子群迭代次数已经达到最大迭代次数或粒子群迄今为止搜索到的最优位置满足预定的最小适应阈值,则执行步骤(10);若不满足,则回到步骤(6);(10)停止迭代,输出最优解和无功配置。
地址 200240 上海市闵行区东川路800号