发明名称 基于多工况时域响应的PID水轮机调速器参数优化方法
摘要 本发明公开了一种基于多工况时域响应的PID水轮机调速器参数优化方法,包括以下步骤:(1)建立水轮机调节系统模型,具体包括:PID水轮机调速器模型,水轮机-引水管道模型和发电机模型;(2)通过现场试验和参数辨识的方式求得水轮机调节系统模型参数;(3)设计基于多工况时域响应的综合适应度函数,确定优化目标为:找到最优PID控制参数,使得综合适应度函数值达到最小;(4)利用智能优化算法求取最优控制参数。该方法求取的控制参数使得系统能够在不同的工况下都保持满意的动态特性,增强了系统鲁棒性。
申请公布号 CN103590969B 申请公布日期 2016.04.20
申请号 CN201310588803.9 申请日期 2013.11.20
申请人 华中科技大学 发明人 刘昌玉;何雪松;颜秋容;王湛;张恩博
分类号 F03B15/00(2006.01)I 主分类号 F03B15/00(2006.01)I
代理机构 华中科技大学专利中心 42201 代理人 朱仁玲
主权项 一种基于多工况时域响应的PID水轮机调速器参数优化方法,具体包括以下步骤:(1)在某一具体工况C<sub>m</sub>下,根据最大超调量M、稳定时间T以及振荡次数N这三个时域响应指标,对一组给定的PID控制参数K,定义单一工况适应度函数为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>f</mi><msub><mi>C</mi><mi>m</mi></msub></msub><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>a</mi><mn>1</mn></msub><mi>M</mi><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>a</mi><mn>2</mn></msub><mi>T</mi><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>a</mi><mn>3</mn></msub><mi>N</mi><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>a</mi><mn>4</mn></msub><msubsup><mo>&Integral;</mo><mn>0</mn><mrow><mo>+</mo><mi>&infin;</mi></mrow></msubsup><mo>|</mo><mi>e</mi><mrow><mo>(</mo><mi>K</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>|</mo><mi>d</mi><mi>t</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000843950700000011.GIF" wi="1518" he="102" /></maths>其中,系数a<sub>i</sub>用于各个指标之间的加权,i=1,2,3,4;函数e(·)表示频率偏差;K为PID控制器比例系数K<sub>p</sub>、积分系数K<sub>i</sub>和微分系数K<sub>d</sub>构成的向量,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>K</mi><mover><mo>=</mo><mi>&Delta;</mi></mover><mrow><mo>(</mo><msub><mi>K</mi><mi>p</mi></msub><mo>,</mo><msub><mi>K</mi><mi>i</mi></msub><mo>,</mo><msub><mi>K</mi><mi>d</mi></msub><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000843950700000012.GIF" wi="381" he="84" /></maths>t表示时间;(2)在多个工况C<sub>m</sub>(m=1,2,…,M)的情况下,定义PID控制参数K的综合适应度函数为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>F</mi><mi>i</mi><mi>t</mi><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow><mo>=</mo><munder><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mrow><mi>m</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>M</mi></mrow></munder><mo>{</mo><msub><mi>f</mi><msub><mi>C</mi><mi>m</mi></msub></msub><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000843950700000013.GIF" wi="1540" he="102" /></maths>(3)确定优化目标为:找到最优PID控制参数K<sup>*</sup>,使得综合适应度函数值达到最小,即:<img file="FDA0000843950700000014.GIF" wi="486" he="94" />
地址 430074 湖北省武汉市洪山区珞喻路1037号