发明名称 一种基于自动微分的含统一潮流控制器的最优潮流算法
摘要 本发明公布了一种基于自动微分的含统一潮流控制器的最优潮流算法。随着统一潮流控制器(UPFC)的灵活控制功能受到越来越多的关注,其工程化应用方兴未艾,研究计及UPFC的电力系统最优潮流算法(OPF)模型成为一个热点。UPFC灵活多样的控制功能使算法编程实现的工作量变得很大,自动微分技术(AD)技术可以自动生成内点法计算过程中需要用到的雅可比矩阵和海森矩阵,从而提高开发效率,减少编程工作量。本发明首先建立了含UPFC的电力系统OPF模型,然后采用AD技术和内点法相结合的方法进行求解,算例结果表明,该方法的收敛性强,适应性好,实现简单。
申请公布号 CN106253288A 申请公布日期 2016.12.21
申请号 CN201610567278.6 申请日期 2016.07.19
申请人 河海大学 发明人 孙永辉;钱臻;何天雨;卫志农
分类号 H02J3/06(2006.01)I 主分类号 H02J3/06(2006.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 刘渊
主权项 一种基于自动微分的含统一潮流控制器的最优潮流算法,其特征在于,所述方法在计算机中依次按以下步骤实现:(1)获得电力系统的网络参数信息;(2)程序初始化,包括:对算法中的状态总变量x设置初值、对拉格朗日乘子y和松弛变量u和w设置初值、设置迭代计数器k=0,设置最大迭代次数K<sub>max</sub>、收敛精度ε、形成节点导纳矩阵Y;(3)计算互补间隙Gap,判断其是否满足精度要求,若满足,则输出最优解,结束循环,否则,继续;(4)利用自动微分技术计算目标函数f(x)、等式约束h(x)、不等式约束<img file="FDA0001053185230000011.GIF" wi="122" he="83" />的一阶导数雅可比矩阵<img file="FDA0001053185230000012.GIF" wi="443" he="79" />以及目标函数、等式约束、不等式约束的二阶导数海森矩阵<img file="FDA0001053185230000013.GIF" wi="699" he="79" />然后计算各常数项L′<sub>x</sub>、L<sub>y</sub>、L<sub>w</sub>、<img file="FDA0001053185230000014.GIF" wi="84" he="70" />并根据以下方程求解总变量x、拉格朗日乘子y和松弛变量u、w的增量Δx、Δy、Δu、Δw;(5)确定原始变量和对偶变量的迭代步长:<maths num="0001"><math><![CDATA[<mrow><msub><mi>&alpha;</mi><mi>p</mi></msub><mo>=</mo><mn>0.9995</mn><mi>m</mi><mi>i</mi><mi>n</mi><mo>{</mo><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mi>i</mi></munder><mrow><mo>(</mo><mfrac><mrow><mo>-</mo><msub><mi>u</mi><mi>i</mi></msub></mrow><mrow><msub><mi>&Delta;u</mi><mi>i</mi></msub></mrow></mfrac><msub><mo>|</mo><mrow><msub><mi>&Delta;u</mi><mi>i</mi></msub><mo>&lt;</mo><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>,</mo><mn>1</mn><mo>}</mo></mrow>]]></math><img file="FDA0001053185230000015.GIF" wi="806" he="206" /></maths><maths num="0002"><math><![CDATA[<mrow><msub><mi>&alpha;</mi><mi>d</mi></msub><mo>=</mo><mn>0.9995</mn><mi>m</mi><mi>i</mi><mi>n</mi><mo>{</mo><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mi>i</mi></munder><mrow><mo>(</mo><mfrac><mrow><mo>-</mo><msub><mi>w</mi><mi>i</mi></msub></mrow><mrow><msub><mi>&Delta;w</mi><mi>i</mi></msub></mrow></mfrac><msub><mo>|</mo><mrow><msub><mi>&Delta;w</mi><mi>i</mi></msub><mo>&gt;</mo><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>,</mo><mn>1</mn><mo>}</mo></mrow>]]></math><img file="FDA0001053185230000016.GIF" wi="821" he="207" /></maths>其中:u<sub>i</sub>、Δu<sub>i</sub>分别表示松弛变量向量u及其增量Δu的第i个元素,w<sub>i</sub>、Δw<sub>i</sub>分别表示松弛变量向量w及其增量Δw的第i个元素,i=1,2...,r;(6)按照下式更新所有变量和拉格朗日乘子:<maths num="0003"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><msup><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></msup></mtd></mtr><mtr><mtd><msup><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></msup></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msup><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup></mtd></mtr><mtr><mtd><msup><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup></mtd></mtr></mtable></mfenced><mo>+</mo><msub><mi>&alpha;</mi><mi>p</mi></msub><mfenced open = "[" close = "]"><mtable><mtr><mtd><mi>&Delta;</mi><mi>x</mi></mtd></mtr><mtr><mtd><mi>&Delta;</mi><mi>u</mi></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001053185230000017.GIF" wi="590" he="184" /></maths><maths num="0004"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><msup><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></msup></mtd></mtr><mtr><mtd><msup><mi>w</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></msup></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msup><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup></mtd></mtr><mtr><mtd><msup><mi>w</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup></mtd></mtr></mtable></mfenced><mo>+</mo><msub><mi>&alpha;</mi><mi>d</mi></msub><mfenced open = "[" close = "]"><mtable><mtr><mtd><mi>&Delta;</mi><mi>y</mi></mtd></mtr><mtr><mtd><mi>&Delta;</mi><mi>w</mi></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001053185230000018.GIF" wi="612" he="183" /></maths>(7)判断迭代次数是否小于最大迭代次数K<sub>max</sub>,若是,则令迭代次数加1,返回步骤(3),否则,输出“计算不收敛”,结束程序。
地址 210024 江苏省南京市江宁区佛城西路8号