发明名称 一种用于机械储能系统的控制方法
摘要 本发明公开了一种用于机械储能系统的控制方法,步骤如下:1)将一个飞轮储能单元与一台风力发电机并联,或与其他多个飞轮储能单元一起组成飞轮储能矩阵;2)转子位置传感器采集飞轮电机的转速;电流传感器采集直流母线电流信号;3)计算转速误差的采样值;4)应用神经元自适应PID控制器作为转速控制器,计算其输出控制量;5)计算直流母线电流误差信号;6)将电流误差信号输入电流控制器,用于控制PWM波的输出,PWM驱动双向变流器控制飞轮电机加速或减速运行。该控制方法从根本上不同于其它的针对三相永磁同步式飞轮驱动电机设计的控制方法,无需进行矢量控制,大大提高了飞轮储能系统的转速控制精度和系统的鲁棒性。
申请公布号 CN103248292B 申请公布日期 2015.11.04
申请号 CN201310201718.2 申请日期 2013.05.27
申请人 重庆大学 发明人 王磊;宋永端;杜晓强;陈柳;王玉兴
分类号 H02P6/16(2006.01)I;H02P6/06(2006.01)I;H02J3/30(2006.01)I 主分类号 H02P6/16(2006.01)I
代理机构 重庆博凯知识产权代理有限公司 50212 代理人 穆祥维
主权项 一种用于机械储能系统的控制方法,其特征在于,该方法包括如下步骤:1)将一个飞轮储能单元与一台风力发电机并联在一条直流母线上,组成一个发电单元;或者将多个飞轮储能单元并联在一条直流母线上,组成飞轮储能矩阵系统,飞轮储能矩阵的直流母线通过双向变流器与风电场的交流母线连接;2)采用转子位置传感器采集飞轮电机的转速,电机实际输出转速为y(k);采用电流传感器采集飞轮储能单元输出的直流母线电流信号I<sub>dc</sub>;3)计算转速误差的采样值e(k):e(k)=y(k)‑r(k),其中,r(k)为飞轮参考转速;4)应用神经元自适应PID控制器作为转速控制器,计算其控制量u(k):先计算变化量△u(k):△u(k)=w<sub>1</sub>[e(k)‑e(k‑1)]+w<sub>2</sub>e(k)+w<sub>3</sub>[e(k)‑2e(k‑1)+e(k‑2)];e(k)、e(k‑1)和e(k‑2)为转速误差的采样值,k、k‑1和k‑2分别为采样时刻;w<sub>1</sub>为转速误差的比例权值,w<sub>2</sub>为转速误差的积分权值,w<sub>3</sub>为转速误差的微分权值;再由△u(k)=u(k)‑u(k‑1)计算出控制量u(k),并将其作为直流母线电流的参考输出信号<img file="FDA0000777221890000011.GIF" wi="74" he="82" />输入电流内环控制;其中,转速误差的比例权值w<sub>1</sub>、积分权值w<sub>2</sub>和微分权值w<sub>3</sub>能够通过自主学习过程自适应的调节,控制器的比例、积分、微分状态变量x<sub>i</sub>,i=1,2,3表示为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mo>{</mo><mrow><mtable><mtr><mtd><msub><mi>x</mi><mn>1</mn></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>=</mo><mi>&Delta;</mi><mi>e</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>=</mo><mi>e</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>-</mo><mi>e</mi><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>x</mi><mn>2</mn></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>=</mo><mi>e</mi><mo>(</mo><mi>k</mi><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>x</mi><mn>3</mn></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>=</mo><msup><mi>&Delta;</mi><mn>2</mn></msup><mi>e</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>=</mo><mi>e</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>-</mo><mn>2</mn><mi>e</mi><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo><mo>+</mo><mi>e</mi><mo>(</mo><mi>k</mi><mo>-</mo><mn>2</mn><mo>)</mo></mtd></mtr></mtable><mo>;</mo></mrow></mrow>]]></math><img file="FDA0000777221890000012.GIF" wi="1094" he="293" /></maths>通过调节权值w<sub>i</sub>,为了使系统误差e(k)趋于0,引入一个二次型性能指标函数J(k):<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>J</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mrow><mo>&lsqb;</mo><mi>r</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>e</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000777221890000021.GIF" wi="902" he="171" /></maths>权值的自适应调节目标是J(k)沿着负梯度方向变化,即使误差e(k)变小;J(k)与w<sub>i</sub>(k)的关系可表示为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>&Delta;w</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><msub><mi>&eta;</mi><mi>i</mi></msub><mo>&lsqb;</mo><mo>&part;</mo><mi>J/</mi><mo>&part;</mo><msub><mi>w</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>=</mo><mo>-</mo><msub><mi>&eta;</mi><mi>i</mi></msub><mo>&lsqb;</mo><mo>&part;</mo><mi>J/</mi><mo>&part;</mo><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>&lsqb;</mo><mo>&part;</mo><mi>y</mi><mi>/</mi><mo>&part;</mo><msub><mi>w</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><msub><mi>&eta;</mi><mi>i</mi></msub><mi>e</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&lsqb;</mo><mo>&part;</mo><mi>y</mi><mi>/</mi><mo>&part;</mo><mi>&Delta;</mi><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>&lsqb;</mo><mo>&part;</mo><mi>&Delta;</mi><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>/</mo><mo>&part;</mo><msub><mi>w</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow></mtd></mtr></mtable><mo>;</mo></mrow>]]></math><img file="FDA0000777221890000022.GIF" wi="1453" he="230" /></maths>其中,△w<sub>i</sub>(k)=w<sub>i</sub>(k+1)‑w<sub>i</sub>(k);η<sub>i</sub>,η<sub>i</sub>&gt;0,i=1,2,3分别表示比例、积分、微分环节的学习率;在神经网络有监督的Hebb学习算法下,J(k)与w<sub>i</sub>(k)的关系式重写为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>&Delta;w</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&eta;</mi><mi>i</mi></msub><mi>e</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>x</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&lsqb;</mo><mo>&part;</mo><mi>y</mi><mo>/</mo><mo>&part;</mo><mi>&Delta;</mi><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>3</mn><mo>;</mo></mrow>]]></math><img file="FDA0000777221890000023.GIF" wi="1217" he="118" /></maths>为了简化学习率η<sub>i</sub>初始值的选择,引入K作为学习率的系数,K&gt;0;因此变化量△u(k)的方程被重写为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>&Delta;</mi><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mi>K</mi><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>3</mn></munderover><msubsup><mi>w</mi><mi>i</mi><mn>0</mn></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>x</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000777221890000024.GIF" wi="832" he="229" /></maths>其中,<img file="FDA0000777221890000025.GIF" wi="201" he="122" />为归一化权值,<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mi>w</mi><mi>i</mi><mn>0</mn></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>w</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>/</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>3</mn></munderover><mo>|</mo><msub><mi>w</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo><mo>;</mo></mrow>]]></math><img file="FDA0000777221890000026.GIF" wi="786" he="220" /></maths>在这种情况下,三个变量<img file="FDA0000777221890000027.GIF" wi="243" he="107" />i=1,2,3随着<img file="FDA0000777221890000028.GIF" wi="207" he="110" />i=1,2,3的自适应调节而改变;5)计算电流误差信号e<sub>i</sub>:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>e</mi><mi>i</mi></msub><mo>=</mo><msubsup><mi>I</mi><mrow><mi>d</mi><mi>c</mi></mrow><mo>*</mo></msubsup><mo>-</mo><msub><mi>I</mi><mrow><mi>d</mi><mi>c</mi></mrow></msub><mo>,</mo></mrow>]]></math><img file="FDA0000777221890000029.GIF" wi="346" he="90" /></maths>式中:<img file="FDA00007772218900000210.GIF" wi="78" he="83" />为直流母线电流的参考输出信号,<img file="FDA00007772218900000211.GIF" wi="419" he="93" />I<sub>dc</sub>为飞轮储能单元输出的直流母线电流信号;6)将电流误差信号e<sub>i</sub>输入电流控制器,电流控制器的输出信号d<sup>*</sup>用于控制PWM波的输出,PWM驱动双向变流器控制飞轮电机加速或减速运行。
地址 400044 重庆市沙坪坝区沙正街174号