发明名称 一种基于记忆的分散式直驱风电机组控制方法
摘要 本发明公开了一种基于记忆的分散式直驱风电机组控制方法,将记忆控制与分散式永磁直驱风电机组相结合,分散式的机组保证了风力机最佳布局,记忆控制为机组提供了高效的风能利用和可靠的变浆控制;从仿生学角度设计机组控制器,利用当前响应,联系过去响应,产生新的控制命令,并让控制器在运行中不断优化自身,提高控制精度和准确性;控制策略主要包括电磁转矩控制和变桨执行电流控制,即在低风速时保证最大功率跟踪,以便充分利用风能资源高效发电;在高风速时保证可靠变浆控制,限制超额功率吸收,稳定输出功率;整个系统高效、智能、可靠,具有良好的应用前景。
申请公布号 CN105626377A 申请公布日期 2016.06.01
申请号 CN201510988465.7 申请日期 2015.12.23
申请人 安徽理工大学 发明人 周文庆;周莉
分类号 F03D7/00(2006.01)I 主分类号 F03D7/00(2006.01)I
代理机构 代理人
主权项 一种基于记忆的分散式直驱风电机组控制方法,其特征在于将记忆控制与分散式永磁直驱风电机组相结合;从仿生学角度设计机组控制器算法,利用当前响应,联系过去响应,产生新的控制命令,并让控制器在运行中不断优化自身,提高控制精度和准确性;控制策略主要包括电磁转矩控制和变桨执行电流控制,即在低风速时保证最大功率跟踪,在高风速时保证可靠变浆控制。具体步骤如下:(1)根据变浆动力学特性,变浆过程方程可描述为:<maths num="0001"><math><![CDATA[<mrow><mover><mi>&beta;</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mi>&phi;</mi><mrow><mo>(</mo><mi>&beta;</mi><mo>)</mo></mrow><mo>+</mo><mi>N</mi><mi>i</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000888143560000011.GIF" wi="581" he="110" /></maths>其中<img file="FDA0000888143560000012.GIF" wi="48" he="87" />为变桨距角快慢程度的计量,φ为非线性函数;i为变浆执行电流;N为一常数;最大功率跟踪误差定义为:e=C‑C<sup>*</sup>     (2)其中C为实际风能利用系数,为桨距角β和叶尖速比λ的函数;C<sup>*</sup>为风能利用系数期望值;对式(2)进行微分如下:<maths num="0002"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mover><mi>e</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mover><mi>C</mi><mo>&CenterDot;</mo></mover><mo>-</mo><msup><mover><mi>C</mi><mo>&CenterDot;</mo></mover><mo>*</mo></msup><mo>=</mo><mfrac><mrow><mo>&part;</mo><mi>C</mi></mrow><mrow><mo>&part;</mo><mi>&beta;</mi></mrow></mfrac><mover><mi>&beta;</mi><mo>&CenterDot;</mo></mover><mo>+</mo><mfrac><mrow><mo>&part;</mo><mi>C</mi></mrow><mrow><mo>&part;</mo><mi>&lambda;</mi></mrow></mfrac><mover><mi>&lambda;</mi><mo>&CenterDot;</mo></mover><mo>-</mo><msup><mover><mi>C</mi><mo>&CenterDot;</mo></mover><mo>*</mo></msup><mo>=</mo><mfrac><mrow><mo>&part;</mo><mi>C</mi></mrow><mrow><mo>&part;</mo><mi>&beta;</mi></mrow></mfrac><mover><mi>&beta;</mi><mo>&CenterDot;</mo></mover><mo>+</mo><mfrac><mrow><mo>&part;</mo><mi>C</mi></mrow><mrow><mo>&part;</mo><mi>&lambda;</mi></mrow></mfrac><mrow><mo>(</mo><mfrac><mrow><msub><mover><mi>&omega;</mi><mo>&CenterDot;</mo></mover><mi>t</mi></msub><mi>R</mi><mi>V</mi><mo>-</mo><msub><mi>&omega;</mi><mi>t</mi></msub><mi>R</mi><mover><mi>V</mi><mo>&CenterDot;</mo></mover></mrow><msup><mi>V</mi><mn>2</mn></msup></mfrac><mo>)</mo></mrow><mo>-</mo><msup><mover><mi>C</mi><mo>&CenterDot;</mo></mover><mo>*</mo></msup></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mfrac><mrow><mo>&part;</mo><mi>C</mi></mrow><mrow><mo>&part;</mo><mi>&beta;</mi></mrow></mfrac><mover><mi>&beta;</mi><mo>&CenterDot;</mo></mover><mo>+</mo><mfrac><mrow><mo>&part;</mo><mi>C</mi></mrow><mrow><mo>&part;</mo><mi>&lambda;</mi></mrow></mfrac><mfrac><mrow><msub><mover><mi>&omega;</mi><mo>&CenterDot;</mo></mover><mi>t</mi></msub><mi>R</mi></mrow><mi>V</mi></mfrac><mo>-</mo><mfrac><mrow><mo>&part;</mo><mi>C</mi></mrow><mrow><mo>&part;</mo><mi>&lambda;</mi></mrow></mfrac><mfrac><mrow><msub><mi>&omega;</mi><mi>t</mi></msub><mi>R</mi><mover><mi>V</mi><mo>&CenterDot;</mo></mover></mrow><msup><mi>V</mi><mn>2</mn></msup></mfrac><mo>-</mo><msup><mover><mi>C</mi><mo>&CenterDot;</mo></mover><mo>*</mo></msup><mo>=</mo><mfrac><mrow><mo>&part;</mo><mi>C</mi></mrow><mrow><mo>&part;</mo><mi>&beta;</mi></mrow></mfrac><mover><mi>&beta;</mi><mo>&CenterDot;</mo></mover><mo>+</mo><mfrac><mrow><mo>&part;</mo><mi>C</mi></mrow><mrow><mo>&part;</mo><mi>&lambda;</mi></mrow></mfrac><mfrac><mrow><mi>R</mi><mrow><mo>(</mo><mrow><mfrac><mn>1</mn><mrow><mn>2</mn><mi>&lambda;</mi></mrow></mfrac><msub><mi>C&rho;A</mi><mi>w</mi></msub><msup><mi>RV</mi><mn>2</mn></msup><mo>-</mo><msub><mi>T</mi><mi>d</mi></msub><mo>-</mo><msub><mi>T</mi><mi>e</mi></msub></mrow><mo>)</mo></mrow></mrow><mrow><mi>J</mi><mi>V</mi></mrow></mfrac><mo>-</mo><mfrac><mrow><mo>&part;</mo><mi>C</mi></mrow><mrow><mo>&part;</mo><mi>&lambda;</mi></mrow></mfrac><mfrac><mrow><msub><mi>&omega;</mi><mi>t</mi></msub><mi>R</mi><mover><mi>V</mi><mo>&CenterDot;</mo></mover></mrow><msup><mi>V</mi><mn>2</mn></msup></mfrac><mo>-</mo><msup><mover><mi>C</mi><mo>&CenterDot;</mo></mover><mo>*</mo></msup></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mfrac><mrow><mo>&part;</mo><mi>C</mi></mrow><mrow><mo>&part;</mo><mi>&beta;</mi></mrow></mfrac><mover><mi>&beta;</mi><mo>&CenterDot;</mo></mover><mo>+</mo><mfrac><mrow><mo>&part;</mo><mi>C</mi></mrow><mrow><mo>&part;</mo><mi>&lambda;</mi></mrow></mfrac><mfrac><mrow><msub><mi>RC&rho;A</mi><mi>w</mi></msub><mi>R</mi><mi>V</mi></mrow><mrow><mn>2</mn><mi>&lambda;</mi><mi>J</mi></mrow></mfrac><mo>-</mo><mfrac><mrow><mo>&part;</mo><mi>C</mi></mrow><mrow><mo>&part;</mo><mi>&lambda;</mi></mrow></mfrac><mfrac><mrow><msub><mi>&omega;</mi><mi>t</mi></msub><mi>R</mi><mover><mi>V</mi><mo>&CenterDot;</mo></mover></mrow><msup><mi>V</mi><mn>2</mn></msup></mfrac><mo>-</mo><msup><mover><mi>C</mi><mo>&CenterDot;</mo></mover><mo>*</mo></msup><mo>-</mo><mfrac><mrow><mo>&part;</mo><mi>C</mi></mrow><mrow><mo>&part;</mo><mi>&lambda;</mi></mrow></mfrac><mfrac><mrow><mi>R</mi><mrow><mo>(</mo><mrow><msub><mi>T</mi><mi>d</mi></msub><mo>+</mo><msub><mi>T</mi><mi>e</mi></msub></mrow><mo>)</mo></mrow></mrow><mrow><mi>J</mi><mi>V</mi></mrow></mfrac></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000888143560000013.GIF" wi="1901" he="486" /></maths>把式(1)代入式(3),并拆分再配对得:<maths num="0003"><math><![CDATA[<mrow><mover><mi>e</mi><mo>&CenterDot;</mo></mover><mo>=</mo><msub><mi>f</mi><mn>1</mn></msub><msub><mi>T</mi><mi>e</mi></msub><mo>+</mo><msub><mi>f</mi><mn>2</mn></msub><mi>i</mi><mo>+</mo><msub><mi>f</mi><mn>3</mn></msub><mo>+</mo><msub><mi>f</mi><mn>4</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000888143560000014.GIF" wi="724" he="103" /></maths>其中:<maths num="0004"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>f</mi><mn>1</mn></msub><mo>=</mo><mo>-</mo><mfrac><mrow><mo>&part;</mo><mi>C</mi></mrow><mrow><mo>&part;</mo><mi>&lambda;</mi></mrow></mfrac><mfrac><mi>R</mi><mrow><mi>J</mi><mi>V</mi></mrow></mfrac><mo>,</mo><msub><mi>f</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><mo>&part;</mo><mi>C</mi></mrow><mrow><mo>&part;</mo><mi>&beta;</mi></mrow></mfrac><mi>N</mi><mo>,</mo><msub><mi>f</mi><mn>3</mn></msub><mo>=</mo><mo>-</mo><mfrac><mrow><mo>&part;</mo><mi>C</mi></mrow><mrow><mo>&part;</mo><mi>&lambda;</mi></mrow></mfrac><mfrac><mrow><msub><mi>RT</mi><mi>d</mi></msub></mrow><mrow><mi>J</mi><mi>V</mi></mrow></mfrac></mrow></mtd><mtd><mrow><msub><mi>f</mi><mn>4</mn></msub><mo>=</mo><mfrac><mrow><mo>&part;</mo><mi>C</mi></mrow><mrow><mo>&part;</mo><mi>&lambda;</mi></mrow></mfrac><mfrac><mrow><msub><mi>RC&rho;A</mi><mi>w</mi></msub><mi>R</mi><mi>V</mi></mrow><mrow><mn>2</mn><mi>&lambda;</mi><mi>J</mi></mrow></mfrac><mo>-</mo><mfrac><mrow><mo>&part;</mo><mi>C</mi></mrow><mrow><mo>&part;</mo><mi>&lambda;</mi></mrow></mfrac><mfrac><mrow><msub><mi>&omega;</mi><mi>t</mi></msub><mi>R</mi><mover><mi>V</mi><mo>&CenterDot;</mo></mover></mrow><msup><mi>V</mi><mn>2</mn></msup></mfrac><mo>-</mo><msup><mover><mi>C</mi><mo>&CenterDot;</mo></mover><mo>*</mo></msup><mo>+</mo><mfrac><mrow><mo>&part;</mo><mi>C</mi></mrow><mrow><mo>&part;</mo><mi>&beta;</mi></mrow></mfrac><mi>&phi;</mi></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000888143560000015.GIF" wi="1918" he="166" /></maths>式(4)表明,<img file="FDA0000888143560000016.GIF" wi="40" he="79" />是由变桨距执行电流i和电磁转矩T<sub>e</sub>共同决定的,为减小机组空气动力载荷,增强风电转换效率,采取控制策略:V<sub>in</sub>≤V<V<sub>R</sub>时,运行在低风速区,桨距角不变,进行电磁转矩控制;V<sub>R</sub>≤V≤V<sub>out</sub>时,运行在超额风速区,进行变桨距控制,以限值超额功率吸收;其中V<sub>in</sub>、V<sub>R</sub>、V<sub>out</sub>分别为切入风速、额定风速、切出风速。(2)仿照人类基于记忆的学习过程这种控制方法,利用过去的响应和反馈信息,结合当前不成熟的响应,不断优化控制,提高精度的策略,设计基于记忆的机组控制方法。(21)基于记忆原理构造数学模型如下:C(t)=(1‑σ(t))C<sub>NR</sub>+σ(t)C<sub>ER</sub>  (6)上式中,C(t)为当前响应行为,σ(t)为影响因子,取值为σ(t)∈[0,1],C<sub>NR</sub>和C<sub>ER</sub>分别为系统的固有响应和经验响应;固有响应对应的影响因子为1‑σ(t),经验响应对应为σ(t),它们共同影响系统的当前响应行为,影响比重由影响因子σ(t)所决定;其中,影响因子σ(t)设为如下形式函数:<maths num="0005"><math><![CDATA[<mrow><mi>&sigma;</mi><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>e</mi><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>-</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mi>c</mi><mrow><mo>|</mo><mi>e</mi><mo>|</mo><mo>+</mo><mi>&epsiv;</mi></mrow></mfrac><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000888143560000021.GIF" wi="795" he="175" /></maths>上式中c和ε均为大于零的常数,保证函数的收敛及分母不为零。(22)基于记忆的一阶控制器算法:式(4)改写如下:<maths num="0006"><math><![CDATA[<mrow><mover><mi>e</mi><mo>&CenterDot;</mo></mover><mo>=</mo><msub><mi>f</mi><mn>1</mn></msub><msub><mi>T</mi><mi>e</mi></msub><mo>+</mo><msub><mi>f</mi><mn>2</mn></msub><mi>i</mi><mo>+</mo><msub><mi>f</mi><mn>3</mn></msub><mo>+</mo><msub><mi>f</mi><mn>4</mn></msub><mo>=</mo><mi>L</mi><mo>+</mo><msub><mi>f</mi><mn>1</mn></msub><msub><mi>T</mi><mi>e</mi></msub><mo>+</mo><msub><mi>f</mi><mn>2</mn></msub><mi>i</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000888143560000022.GIF" wi="860" he="93" /></maths>其中L=f<sub>3</sub>+f<sub>4</sub>,为非线性项和扰动项;基于记忆的电磁转矩控制和变浆控制方法数学模型如下:<maths num="0007"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>T</mi><mrow><mi>e</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>=</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><mi>&sigma;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mo>)</mo></mrow><msub><mi>T</mi><mrow><mi>e</mi><mo>,</mo><mi>N</mi><mi>R</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>+</mo><mi>&sigma;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>T</mi><mrow><mi>e</mi><mo>,</mo><mi>E</mi><mi>R</mi><mo>,</mo><mi>k</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>T</mi><mrow><mi>e</mi><mo>,</mo><mi>N</mi><mi>R</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>=</mo><mo>-</mo><msub><mi>k</mi><mn>0</mn></msub><msub><mi>e</mi><mi>k</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>T</mi><mrow><mi>e</mi><mo>,</mo><mi>E</mi><mi>R</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>=</mo><msub><mi>w</mi><mn>0</mn></msub><msub><mi>i</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msubsup><mi>g</mi><mn>1</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msup><mi>T</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mrow><msub><mi>w</mi><mn>1</mn></msub><msub><mi>e</mi><mi>k</mi></msub><mo>+</mo><msub><mi>w</mi><mn>2</mn></msub><msub><mi>e</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000888143560000023.GIF" wi="1037" he="303" /></maths>式(9)为V<sub>in</sub>≤V<V<sub>R</sub>风速区间时的电磁转矩控制算法;<maths num="0008"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>i</mi><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&sigma;</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><msub><mi>i</mi><mrow><mi>N</mi><mi>R</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>+</mo><mi>&sigma;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>i</mi><mrow><mi>E</mi><mi>R</mi><mo>,</mo><mi>k</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>i</mi><mrow><mi>N</mi><mi>R</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>=</mo><mo>-</mo><msub><mi>k</mi><mn>0</mn></msub><msub><mi>e</mi><mi>k</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>i</mi><mrow><mi>E</mi><mi>R</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>=</mo><msub><mi>w</mi><mn>0</mn></msub><msub><mi>i</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msubsup><mi>g</mi><mn>2</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msup><mi>T</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>V</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>V</mi><mo>=</mo><msub><mi>w</mi><mn>1</mn></msub><msup><mi>T</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>e</mi><mi>k</mi></msub><mo>+</mo><msub><mi>w</mi><mn>2</mn></msub><msup><mi>T</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>e</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>g</mi><mrow><mn>1</mn><mi>k</mi></mrow></msub><msub><mi>T</mi><mrow><mi>e</mi><mi>k</mi></mrow></msub><mo>-</mo><msub><mi>g</mi><mrow><mn>1</mn><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msub><mi>T</mi><mrow><mi>e</mi><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000888143560000024.GIF" wi="1166" he="375" /></maths>上述式中:w<sub>0</sub>=1,w<sub>2</sub>=‑2,w<sub>3</sub>=1,w<sub>4</sub>=‑1,w<sub>5</sub>=1;式(10)为V<sub>R</sub>≤V≤V<sub>out</sub>风速区间时的变浆电流控制算法。
地址 232001 安徽省合肥市舜耕中路168号