发明名称 基于STATCOM与SVC的电能质量调节器协调控制方法
摘要 本发明公开了一种基于STATCOM与SVC的电能质量调节器协调控制方法。在单机无穷大电力系统中,根据发电机、STATCOM与SVC的动态模型,推导出了装有STATCOM与SVC的单机无穷大电力系统线性化状态方程,并建立了基于STATCOM与SVC的单机无穷大电力系统的扩展Phillips-Heffron模型。在此基础上,建立了STATCOM与SVC的控制系统。由此推导出线性定常组合控制系统。本发明在能够有效减小两者控制器之间的交互影响,从而保证了两个控制装置能够安全稳定运行,使它们能够高效率地相互配合补偿电网无功功率并维持电压稳定,提高电网的稳定性和电能质量。
申请公布号 CN101183790A 申请公布日期 2008.05.21
申请号 CN200710192662.3 申请日期 2007.12.19
申请人 湖南大学 发明人 罗安;彭双剑;赵伟;帅智康;戴亚培;浣威
分类号 H02J3/00(2006.01);H02J3/16(2006.01);H02J3/18(2006.01) 主分类号 H02J3/00(2006.01)
代理机构 长沙正奇专利事务所有限责任公司 代理人 马强
主权项 1.一种基于STATCOM与SVC的电能质量调节器协调控制方法,其特征在于,依据装有STATCOM与SVC的单机无穷大电力系统实物图建立了基于STATCOM与SVC的无穷大电力系统模型,在单机无穷大电力系统中,根据发电机、STATCOM与SVC的动态模型:<maths num="0001"><![CDATA[<math><mfenced open='' close='}'><mtable><mtr><mtd><mover><mi>&delta;</mi><mo>.</mo></mover><mo>=</mo><msub><mi>&omega;</mi><mn>0</mn></msub><mi>&omega;</mi><mo>,</mo></mtd></mtr><mtr><mtd><mover><mi>&omega;</mi><mo>.</mo></mover><mo>=</mo><mrow><mo>(</mo><msub><mi>P</mi><mi>m</mi></msub><mo>-</mo><msub><mi>P</mi><mi>e</mi></msub><mo>-</mo><mi>D&omega;</mi><mo>)</mo></mrow><mo>/</mo><mi>M</mi><mo>,</mo></mtd></mtr><mtr><mtd><msubsup><mi>E</mi><mi>q</mi><mo>&prime;</mo></msubsup><mo>=</mo><mrow><mo>(</mo><mo>-</mo><msub><mi>E</mi><mi>q</mi></msub><mo>+</mo><msub><mi>E</mi><mi>fd</mi></msub><mo>)</mo></mrow><mo>/</mo><msubsup><mi>T</mi><mrow><mi>d</mi><mn>0</mn></mrow><mo>&prime;</mo></msubsup><mo>,</mo></mtd></mtr><mtr><mtd><mover><msub><mi>E</mi><mi>fd</mi></msub><mo>.</mo></mover><mo>=</mo><mo>[</mo><mo>-</mo><msub><mi>E</mi><mi>fd</mi></msub><mo>+</mo><msub><mi>K</mi><mi>A</mi></msub><mrow><mo>(</mo><msub><mi>V</mi><mrow><mi>t</mi><mn>0</mn></mrow></msub><mo>-</mo><msub><mi>V</mi><mi>t</mi></msub><mo>)</mo></mrow><mo>]</mo><mo>/</mo><msub><mi>T</mi><mi>A</mi></msub></mtd></mtr></mtable></mfenced></math>]]></maths>式中:δ是发电机转子摇摆角(rad);ω为发电机的转子角速度(rad/s);ω<sub>0</sub>是同步角速度,ω<sub>0</sub>=2πf<sub>0</sub>=314rad/s;P<sub>e</sub>是发电机的输出电磁功率;P<sub>m</sub>是发电机的输入机械功率;D是发电机阻尼参数;M是发电机惯性时间常数;E<sub>q</sub>是发电机q轴电势;E<sub>q</sub>′是发电机q轴暂态电势;T<sub>d0</sub>′是发电机d轴暂态短路时间常数;V<sub>b</sub>是无穷大母线电压;x<sub>d</sub>是发电机d轴电抗;x<sub>d</sub>′是发电机d轴暂态电抗;且,P<sub>e</sub>=E<sub>q</sub>′I<sub>tq</sub>+(X<sub>q</sub>-X<sub>d</sub>′)I<sub>td</sub>I<sub>tq</sub>,V<sub>tq</sub>=E<sub>q</sub>′-X<sub>d</sub>′I<sub>td</sub>,E<sub>q</sub>=E<sub>q</sub>′+(X<sub>d</sub>-X<sub>d</sub>′)I<sub>td</sub>,V<sub>t</sub><sup>2</sup>=V<sub>td</sub><sup>2</sup>+V<sub>tq</sub><sup>2</sup>,V<sub>tq</sub>=X<sub>q</sub>I<sub>tq</sub>.导出装有STATCOM与SVC的单机无穷大电力系统线性化状态方程:<img file="S2007101926623C00021.GIF" wi="1310" he="581" /><maths num="0002"><![CDATA[<math><mrow><mo>+</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>-</mo><mfrac><msub><mi>K</mi><mi>pa</mi></msub><mi>M</mi></mfrac></mtd><mtd><mo>-</mo><mfrac><msub><mi>K</mi><mi>pc</mi></msub><mi>M</mi></mfrac></mtd><mtd><mo>-</mo><mfrac><msub><mi>K</mi><mi>p&psi;</mi></msub><mi>M</mi></mfrac></mtd></mtr><mtr><mtd><mo>-</mo><mfrac><msub><mi>K</mi><mi>qa</mi></msub><msubsup><mi>T</mi><mrow><mi>d</mi><mn>0</mn></mrow><mo>&prime;</mo></msubsup></mfrac></mtd><mtd><mo>-</mo><mfrac><msub><mi>K</mi><mi>qc</mi></msub><msubsup><mi>T</mi><mrow><mi>d</mi><mn>0</mn></mrow><mo>&prime;</mo></msubsup></mfrac></mtd><mtd><mo>-</mo><mfrac><msub><mi>K</mi><mi>q&psi;</mi></msub><msubsup><mi>T</mi><mrow><mi>d</mi><mn>0</mn></mrow><mo>&prime;</mo></msubsup></mfrac></mtd></mtr><mtr><mtd><mo>-</mo><mfrac><mrow><msub><mi>K</mi><mi>A</mi></msub><msub><mi>K</mi><mi>va</mi></msub></mrow><msub><mi>T</mi><mi>A</mi></msub></mfrac></mtd><mtd><mo>-</mo><mfrac><mrow><msub><mi>K</mi><mi>A</mi></msub><msub><mi>K</mi><mi>vc</mi></msub></mrow><msub><mi>T</mi><mi>A</mi></msub></mfrac></mtd><mtd><mo>-</mo><mfrac><mrow><msub><mi>K</mi><mi>A</mi></msub><msub><mi>K</mi><mi>v&psi;</mi></msub></mrow><msub><mi>T</mi><mi>A</mi></msub></mfrac></mtd></mtr><mtr><mtd><msub><mi>K</mi><mi>da</mi></msub></mtd><mtd><msub><mi>K</mi><mi>dc</mi></msub></mtd><mtd><msub><mi>K</mi><mi>d&psi;</mi></msub></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&Delta;&alpha;</mi></mtd></mtr><mtr><mtd><mi>&Delta;c</mi></mtd></mtr><mtr><mtd><mi>&Delta;&psi;</mi></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>并建立了基于STATCOM与SVC的单机无穷大电力系统的扩展Phillips-Heffron模型,在扩展Phillips-Heffron的基础上,建立了STATCOM与SVC的控制系统,该控制系统是一个三输入、三输出系统,对每一对输入、输出设计一个单输入单输出的控制器,由此推导出线性定常组合控制系统:<img file="S2007101926623C00023.GIF" wi="400" he="181" />式中,<maths num="0003"><![CDATA[<math><mrow><mi>A</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>&omega;</mi><mn>0</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>-</mo><mfrac><msub><mi>K</mi><mn>1</mn></msub><mi>M</mi></mfrac></mtd><mtd><mo>-</mo><mfrac><mi>D</mi><mi>M</mi></mfrac></mtd><mtd><mo>-</mo><mfrac><msub><mi>K</mi><mn>2</mn></msub><mi>M</mi></mfrac></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mfrac><msub><mi>K</mi><mi>pDC</mi></msub><mi>M</mi></mfrac></mtd></mtr><mtr><mtd><mo>-</mo><mfrac><msub><mi>K</mi><mn>4</mn></msub><msubsup><mi>T</mi><mrow><mi>d</mi><mn>0</mn></mrow><mo>&prime;</mo></msubsup></mfrac></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mfrac><msub><mi>K</mi><mn>3</mn></msub><msubsup><mi>T</mi><mrow><mi>d</mi><mn>0</mn></mrow><mo>&prime;</mo></msubsup></mfrac></mtd><mtd><mfrac><mn>1</mn><msubsup><mi>T</mi><mrow><mi>d</mi><mn>0</mn></mrow><mo>&prime;</mo></msubsup></mfrac></mtd><mtd><mo>-</mo><mfrac><msub><mi>K</mi><mi>qDC</mi></msub><msubsup><mi>T</mi><mrow><mi>d</mi><mn>0</mn></mrow><mo>&prime;</mo></msubsup></mfrac></mtd></mtr><mtr><mtd><mo>-</mo><mfrac><mrow><msub><mi>K</mi><mi>A</mi></msub><msub><mi>K</mi><mn>5</mn></msub></mrow><msub><mi>T</mi><mi>A</mi></msub></mfrac></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mfrac><mrow><msub><mi>K</mi><mi>A</mi></msub><msub><mi>K</mi><mn>6</mn></msub></mrow><msub><mi>T</mi><mi>A</mi></msub></mfrac></mtd><mtd><mo>-</mo><mfrac><mn>1</mn><msub><mi>T</mi><mi>A</mi></msub></mfrac></mtd><mtd><mo>-</mo><mfrac><mrow><msub><mi>K</mi><mi>A</mi></msub><msub><mi>K</mi><mi>vDC</mi></msub></mrow><msub><mi>T</mi><mi>A</mi></msub></mfrac></mtd></mtr><mtr><mtd><msub><mi>K</mi><mn>7</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>K</mi><mn>8</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>K</mi><mi>dDC</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mi>B</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>-</mo><mfrac><msub><mi>K</mi><mi>p&alpha;</mi></msub><mi>M</mi></mfrac></mtd><mtd><mo>-</mo><mfrac><msub><mi>K</mi><mi>pc</mi></msub><mi>M</mi></mfrac></mtd><mtd><mo>-</mo><mfrac><msub><mi>K</mi><mi>p&psi;</mi></msub><mi>M</mi></mfrac></mtd></mtr><mtr><mtd><mo>-</mo><mfrac><msub><mi>K</mi><mi>q&alpha;</mi></msub><msubsup><mi>T</mi><mrow><mi>d</mi><mn>0</mn></mrow><mo>&prime;</mo></msubsup></mfrac></mtd><mtd><mo>-</mo><mfrac><msub><mi>K</mi><mi>qc</mi></msub><msubsup><mi>T</mi><mrow><mi>d</mi><mn>0</mn></mrow><mo>&prime;</mo></msubsup></mfrac></mtd><mtd><mo>-</mo><mfrac><msub><mi>K</mi><mi>q&psi;</mi></msub><msubsup><mi>T</mi><mrow><mi>d</mi><mn>0</mn></mrow><mo>&prime;</mo></msubsup></mfrac></mtd></mtr><mtr><mtd><mo>-</mo><mfrac><mrow><msub><mi>K</mi><mi>A</mi></msub><msub><mi>K</mi><mi>v&alpha;</mi></msub></mrow><msub><mi>T</mi><mi>A</mi></msub></mfrac></mtd><mtd><mo>-</mo><mfrac><mrow><msub><mi>K</mi><mi>A</mi></msub><msub><mi>K</mi><mi>vc</mi></msub></mrow><msub><mi>T</mi><mi>A</mi></msub></mfrac></mtd><mtd><mo>-</mo><mfrac><mrow><msub><mi>K</mi><mi>A</mi></msub><msub><mi>K</mi><mi>v&psi;</mi></msub></mrow><msub><mi>T</mi><mi>A</mi></msub></mfrac></mtd></mtr><mtr><mtd><msub><mi>K</mi><mi>d&alpha;</mi></msub></mtd><mtd><msub><mi>K</mi><mi>dc</mi></msub></mtd><mtd><msub><mi>K</mi><mi>d&psi;</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msup><mi>C</mi><mo>&prime;</mo></msup><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn><mo>-</mo><msub><mi>B</mi><mn>6</mn></msub><msub><mi>K</mi><mi>DCP</mi></msub></mtd><mtd><mn>1</mn><mo>-</mo><msub><mi>B</mi><mn>5</mn></msub><msub><mi>K</mi><mi>ACP</mi></msub></mtd><mtd><mn>1</mn><mo>-</mo><msub><mi>B</mi><mn>3</mn></msub><msub><mi>K</mi><mi>&alpha;P</mi></msub></mtd></mtr><mtr><mtd><mn>1</mn><mo>-</mo><msub><mi>D</mi><mn>6</mn></msub><msub><mi>K</mi><mi>DCP</mi></msub></mtd><mtd><mn>1</mn><mo>-</mo><msub><mi>D</mi><mn>5</mn></msub><msub><mi>K</mi><mi>ACP</mi></msub></mtd><mtd><mn>1</mn><mo>-</mo><msub><mi>D</mi><mn>3</mn></msub><msub><mi>K</mi><mi>&alpha;P</mi></msub></mtd></mtr><mtr><mtd><mn>1</mn><mo>-</mo><msub><mi>C</mi><mn>6</mn></msub><msub><mi>K</mi><mi>DCP</mi></msub></mtd><mtd><mn>1</mn><mo>-</mo><msub><mi>C</mi><mn>5</mn></msub><msub><mi>K</mi><mi>ACP</mi></msub></mtd><mtd><mn>1</mn><mo>-</mo><msub><mi>C</mi><mn>3</mn></msub><msub><mi>K</mi><mi>&alpha;P</mi></msub></mtd></mtr></mtable></mfenced><mrow><mo>-</mo><mn>1</mn></mrow></msup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>B</mi><mn>1</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>B</mi><mn>2</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>B</mi><mn>4</mn></msub></mtd></mtr><mtr><mtd><msub><mi>D</mi><mn>1</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>D</mi><mn>2</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>D</mi><mn>4</mn></msub></mtd></mtr><mtr><mtd><msub><mi>C</mi><mn>1</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>C</mi><mn>2</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>C</mi><mn>4</mn></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>X=[Δδ Δω ΔE<sub>q</sub>′ ΔE<sub>fd</sub> ΔV<sub>DC</sub>]<sup>T</sup>,U=[Δα Δc Δψ]<sup>T</sup>,Y=[ΔV<sub>DC</sub> ΔV<sub>l</sub> ΔV<sub>s</sub>]<sup>T</sup>;线性定常组合系统(X(0)=X<sub>0</sub>)在给定信息结构S的约束下,求静态输出反馈控制规律U(t)=KY(t),其中K为具有特定信息S的常数矩阵,使得线性二次型LQ性能指标<maths num="0006"><![CDATA[<math><mrow><mi>J</mi><mo>=</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mo>&infin;</mo></msubsup><mrow><mo>(</mo><msup><mi>X</mi><mi>T</mi></msup><mi>QX</mi><mo>+</mo><msup><mi>U</mi><mi>T</mi></msup><mi>RU</mi><mo>)</mo></mrow><mi>dt</mi></mrow></math>]]></maths>为极小,其中Q,R为相应维数的对称正半定矩阵,如果初始状态X<sub>0</sub>是随机变量,且E<sub>0</sub>=E(X<sub>0</sub>X<sub>0</sub><sup>T</sup>),其中E(·)是期望运算,考虑目标函数的期望值<maths num="0007"><![CDATA[<math><mrow><mover><mi>J</mi><mo>^</mo></mover><mo>=</mo><mi>E</mi><mrow><mo>(</mo><mi>J</mi><mo>)</mo></mrow><mo>=</mo><mi>tr</mi><mrow><mo>(</mo><mi>P</mi><msub><mi>E</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>为简化计算取E<sub>0</sub>=I(单位阵),则<img file="S2007101926623C00034.GIF" wi="44" he="73" />取得极小值的必要条件是(P,K,V)满足S信息结构下推广的Levine-Athans方程组:<maths num="0008"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><mrow><mo>&PartialD;</mo><mover><mi>J</mi><mo>^</mo></mover></mrow><mrow><mo>&PartialD;</mo><mi>K</mi></mrow></mfrac><mo>=</mo><mn>2</mn><mrow><mo>(</mo><mi>RKCV</mi><msup><mi>C</mi><mi>T</mi></msup><mo>+</mo><msup><mi>B</mi><mi>T</mi></msup><mi>PV</mi><msup><mi>C</mi><mi>T</mi></msup><mo>)</mo></mrow><mi>s</mi><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><mi>P</mi><mrow><mo>(</mo><mi>A</mi><mo>+</mo><mi>BKC</mi><mo>)</mo></mrow><mo>+</mo><msup><mrow><mo>(</mo><mi>A</mi><mo>+</mo><mi>BKC</mi><mo>)</mo></mrow><mi>T</mi></msup><mi>P</mi><mo>+</mo><mi>Q</mi><mo>+</mo><msup><mi>C</mi><mi>T</mi></msup><msup><mi>K</mi><mi>T</mi></msup><mi>RKC</mi><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><mrow><mo>(</mo><mi>A</mi><mo>+</mo><mi>BKC</mi><mo>)</mo></mrow><mi>V</mi><mo>+</mo><msup><mrow><mo>(</mo><mi>A</mi><mo>+</mo><mi>BKC</mi><mo>)</mo></mrow><mi>T</mi></msup><mo>+</mo><msub><mi>E</mi><mn>0</mn></msub><mo>=</mo><mn>0</mn></mtd></mtr></mtable></mfenced></math>]]></maths>上述Levine-Athans方程组采用梯度算法和直接迭代法求解,对于基于STATCOM与SVC的单机无穷大电力系统线性状态方程,由于各局部控制器均采用本地输出反馈来实现,相应的信息结构矩阵退化为分块对角阵;对单机无穷大电力系统在典型运行点获得相应的小范围线性化状态空间模型,选取LQ性能指标中的权系数矩阵Q,得到协调控制器的静态反馈增益K。
地址 410082湖南省长沙市岳麓山