发明名称 一种水轮机调节系统的参数辨识方法
摘要 本发明公开了一种水轮机调节系统的参数辨识方法,用于在水轮机调节系统中对系统参数进行有效辨识。根据水轮机调节系统建立仿真模型,然后依据该仿真模型建立采用实际系统输出与辨识系统输出的加权误差平方和作为参数辨识的目标函数,运用本发明设计的优选方法求解目标函数得到最优控制参数。本发明设计的水轮机调节系统参数辨识方法,采用一种新型启发式优化算法优化目标函数,可以搜索到更小的目标函数值,能得到更精确的辨识参数。更精确的辨识参数使得对应的辨识系统与实际系统输出吻合较好。
申请公布号 CN105298734A 申请公布日期 2016.02.03
申请号 CN201510759863.1 申请日期 2015.11.10
申请人 华中科技大学 发明人 李超顺;董伟;毛翼丰;张楠;罗萌;王文潇
分类号 F03B15/00(2006.01)I 主分类号 F03B15/00(2006.01)I
代理机构 华中科技大学专利中心 42201 代理人 曹葆青
主权项 一种水轮机调节系统参数辨识方法,其特征在于,所述方法包括如下步骤:步骤(1):建立水轮机调节系统仿真模型,确定待辨识参数:所述的水轮机调节系统包括PID控制器、电液随动系统、引水系统、水轮机、发电机及负载;在所述水轮机调节系统中,当频率给定发生改变时,PID控制器根据机组频率偏差产生调节控制信号驱动电液随动系统,改变水轮机导叶开度,水轮机进口流量随之发生改变,在水轮机导叶开度改变瞬间,引水系统中发生水锤现象,导致水轮机蜗壳压力发生变化;水轮机进口流量和蜗壳压力的变化,会使水轮机力矩发生改变,从而使水轮机力矩与发电机的负载阻力矩产生差值,发电机转速随之发生变化,转速改变同步地调节了频率大小,达到调节机组频率的目的;需要辨识的参数向量为<img file="FDA0000843227010000011.GIF" wi="684" he="102" />其中k<sub>p</sub>,k<sub>i</sub>,k<sub>d</sub>为PID控制器参数,T<sub>y1</sub>和T<sub>y</sub>为电液随动系统的两个时间常数,h<sub>w</sub>为引水管道特征常数,T<sub>r</sub>为水锤反射时间常数,T<sub>a</sub>为机组惯性时间常数,e<sub>n</sub>为发电机与水轮机自调节系数之差;步骤(2):采集水轮机调节系统动态过程数据:对水轮机调节系统进行频率扰动试验,采集系统动态过程数据,动态过程数据包括导叶开度y<sub>t</sub>,机组频率x;步骤(3):建立水轮机调节系统辨识目标函数:采用实际系统输出与辨识系统输出的加权误差平方和作为参数辨识的目标函数,目标函数定义如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>min</mi><mi> </mi><msub><mi>f</mi><mrow><mi>W</mi><mi>M</mi><mi>S</mi><mi>E</mi></mrow></msub><mrow><mo>(</mo><mover><mi>&theta;</mi><mo>^</mo></mover><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>s</mi></msub></munderover><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>s</mi></munderover><msub><mi>w</mi><mi>j</mi></msub><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>j</mi></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>-</mo><msub><mover><mi>y</mi><mo>^</mo></mover><mi>j</mi></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FDA0000843227010000012.GIF" wi="796" he="151" /></maths>其中N<sub>S</sub>为系统输出采样点数,s为系统输出个数,w=[w<sub>1</sub>,w<sub>2</sub>]为权重,<img file="FDA0000843227010000021.GIF" wi="658" he="110" />是待辨识参数,y<sub>j</sub>(k)∈{y<sub>t</sub>(k),x(k)}为实际系统输出;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>y</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&Element;</mo><mo>{</mo><msub><mover><mi>y</mi><mo>^</mo></mover><mi>t</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>}</mo></mrow>]]></math><img file="FDA0000843227010000022.GIF" wi="494" he="70" /></maths>为辨识系统输出;步骤(4):运用启发式优化算法求解水轮机调节系统参数辨识的目标函数,得到待辨识系统参数。
地址 430074 湖北省武汉市洪山区珞喻路1037号