发明名称 一种基于量化的广域电力系统输出反馈控制方法
摘要 一种基于量化的广域电力系统输出反馈控制方法,步骤如下:1)在广域测量系统当中,信息的采集是通过PMU来完成的,PMU以固定的采样周期对电力系统状态量进行测量,建立以T为采样周期的闭环电力系统离散化状态空间模型。2)在实际系统中,系统的状态难以全部测得。3)采用Lyapunov稳定性理论及线性矩阵不等式方法对得到的电力系统模型进行稳定性分析,使得满足不等式(15),然后设计控制器增益分析验证。本发明能够保证电力系统在带有量化器情况下的稳定性,有效缓解了通信网络拥塞问题,减少了在通信网络上的建设费用,同时节约数据存储中心的存储空间,有利于系统的扩展,增加经济效益。
申请公布号 CN103414245B 申请公布日期 2015.10.28
申请号 CN201310220056.3 申请日期 2013.06.04
申请人 浙江工业大学 发明人 俞立;白丽叶;刘安东;张文安;孙辉;凌荣耀
分类号 H02J13/00(2006.01)I;G06F17/50(2006.01)I 主分类号 H02J13/00(2006.01)I
代理机构 杭州天正专利事务所有限公司 33201 代理人 王兵;黄美娟
主权项 一种基于量化的广域电力系统输出反馈控制方法,其特征在于:所述控制方法包括以下步骤:1)在广域测量系统当中,信息的采集是通过PMU来完成的,PMU以固定的采样周期对电力系统状态量进行测量,以T为采样周期的闭环电力系统离散化状态空间模型为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open = '{' close = ''><mtable><mtr><mtd><mi>x</mi><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo><mo>=</mo><mi>A</mi><mi>x</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>+</mo><mi>B</mi><mi>u</mi><mo>(</mo><mi>k</mi><mo>)</mo></mtd></mtr><mtr><mtd><mi>y</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>=</mo><mi>C</mi><mi>x</mi><mo>(</mo><mi>k</mi><mo>)</mo></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000747948040000011.GIF" wi="529" he="157" /></maths>其中,k表示第k个传输周期kT;A,B表示电机的常数矩阵;C为系统的输出矩阵;x(k+1)表示广域电力系统在第k+1周期的状态变量;u(k)为系统第k周期的输入变量;y(k)表示系统在第k周期的输出变量;2)首先根据输出信号对系统的状态进行观测,然后设计基于观测器的输出反馈控制器对系统进行控制,并在电力系统闭环系统中加入量化器,对PMU的输出量进行量化,采用对数量化器:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open = '' close = ''><mtable><mtr><mtd><mrow><msub><mi>y</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mi>f</mi><mrow><mo>(</mo><mi>y</mi><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><msub><mi>u</mi><mi>i</mi></msub><mo>,</mo></mrow></mtd><mtd><mrow><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><mi>&delta;</mi></mrow></mfrac><msub><mi>u</mi><mi>i</mi></msub><mo>&lt;</mo><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&le;</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>-</mo><mi>&delta;</mi></mrow></mfrac><msub><mi>u</mi><mi>i</mi></msub><mo>;</mo><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&gt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mi>f</mi><mrow><mo>(</mo><mo>-</mo><mi>y</mi><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo><mo>,</mo></mrow></mtd><mtd><mrow><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000747948040000012.GIF" wi="1124" he="317" /></maths>相应的量化等级定义为:U={±u<sub>(i)</sub>,u<sub>(i)</sub>=ρ<sup>i</sup>u<sub>(0)</sub>,i=±1,±2,…}∪{±u<sub>(0)</sub>}∪{0}0<ρ<1,u<sub>(0)</sub>>0其中,y<sub>q</sub>(k)为量化后的输出变量,u<sub>(i)</sub>为第i级量化水平、u<sub>(0)</sub>为初始量化水平,ρ为量化密度,<img file="FDA0000747948040000021.GIF" wi="229" he="140" />带有量化器的观测器结构为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mi>A</mi><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>B</mi><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>L</mi><mo>(</mo><msub><mi>y</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mover><mi>y</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow>]]></math><img file="FDA0000747948040000022.GIF" wi="871" he="94" /></maths>其中,<img file="FDA0000747948040000023.GIF" wi="168" he="71" />为k+1时刻系统状态的估计值、L为观测器增益、<img file="FDA0000747948040000024.GIF" wi="111" he="71" />为观测器的输出变量;3)采用Lyapunov稳定性理论及线性矩阵不等式方法可得到保证整个广域电力系统在量化误差存在的情况下保持稳态运行的充分条件:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfenced open = '[' close = ']'><mtable><mtr><mtd><mrow><mo>-</mo><mi>P</mi><mo>+</mo><msup><mi>&epsiv;&delta;</mi><mn>2</mn></msup><msup><mi>C</mi><mi>T</mi></msup><mi>C</mi></mrow></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mi>Q</mi></mrow></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd></mtr><mtr><mtd><mrow><mi>P</mi><mi>A</mi><mo>+</mo><mi>Y</mi></mrow></mtd><mtd><mrow><mo>-</mo><mi>Y</mi></mrow></mtd><mtd><mrow><mo>-</mo><mi>P</mi></mrow></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mi>Q</mi><mi>A</mi><mo>+</mo><mi>Z</mi><mi>C</mi></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mi>Q</mi></mrow></mtd><mtd><mo>*</mo></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><msup><mi>Z</mi><mi>T</mi></msup></mtd><mtd><mrow><mo>-</mo><mi>&epsiv;</mi><mi>I</mi></mrow></mtd></mtr></mtable></mfenced><mo>&lt;</mo><mn>0</mn></mrow>]]></math><img file="FDA0000747948040000025.GIF" wi="1009" he="387" /></maths>其中,P为Lyapunov函数的状态权重矩阵,Q为Lyapunov函数的误差权重矩阵,ε为1维参数变量,矩阵Y为Y=PBK,矩阵Z为Z=QL,通过Matlab中的LMI工具箱可解出矩阵Y和Z,则所求的控制器增益为K=(PB)<sup>‑1</sup>Y,再将得到的控制器加入系统中,使系统保持稳定。
地址 310014 浙江省杭州市下城区潮王路18号