发明名称 一种汽轮机调速系统控制参数整定方法
摘要 一种汽轮机调速系统控制参数整定方法,包括以下步骤:1)建立包括了电网模型时的单机无穷大系统的汽轮机调速系统的数学模型及传递函数;2)当控制参数K<sub>i</sub>和K<sub>p</sub>在0~10范围内发生变化时,不断求出给定参数下系统传递函数特征根的数值解;特征根表达式为:s=A+Bj;3)再求出每组特征根特征参数ξ的大小,并获得其中值最小的ξ;4)绘制每组特征参数ξ最小值随控制参数K<sub>i</sub>和K<sub>p</sub>变化趋势图;5)根据特征参数ξ的正负来判断系统的稳定与否,从图中获取使系统稳定的控制参数的取值范围,采用试凑法进行微调,整定出真正实用于机组的最优控制参数。本发明对于进一步规范汽轮机控制系统参数设置,减少系统低频振荡的发生具有十分重要的意义。
申请公布号 CN104063584B 申请公布日期 2017.02.08
申请号 CN201410254317.8 申请日期 2014.06.10
申请人 国家电网公司;湖北省电力公司电力科学研究院;华中科技大学 发明人 李阳海;潘剑;刘魏然;杨涛;高伟;黄树红
分类号 G06F19/00(2011.01)I;G05B13/04(2006.01)I 主分类号 G06F19/00(2011.01)I
代理机构 武汉楚天专利事务所 42113 代理人 雷速
主权项 一种汽轮机调速系统控制参数整定方法,其特征在于包括如下步骤:步骤1:建立包括了电网模型时的单机无穷大系统的汽轮机调速系统的数学模型,利用现有工况系统中相应的对象模型识别模块得到被控对象的数学模型以及各参数,各参数之间的传递函数亦即包括电网模型时的汽轮机调速系统的数学模型的传递函数,如下式所示:<maths num="0001"><math><![CDATA[<mrow><mo>-</mo><msub><mi>&Delta;P</mi><mi>m</mi></msub><mo>=</mo><mrow><mo>(</mo><mfrac><mn>1</mn><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>T</mi><mn>1</mn></msub><mi>s</mi><mo>)</mo><mi>&delta;</mi></mrow></mfrac><mo>(</mo><mrow><mn>1</mn><mo>+</mo><msub><mi>G</mi><mrow><mi>p</mi><mi>i</mi><mi>d</mi></mrow></msub></mrow><mo>)</mo><mo>+</mo><mfrac><mn>1</mn><mi>s</mi></mfrac><msub><mi>G</mi><mn>1</mn></msub><msub><mi>G</mi><mrow><mi>p</mi><mi>i</mi><mi>d</mi></mrow></msub><mo>)</mo></mrow><msub><mi>G</mi><mi>e</mi></msub><msub><mi>G</mi><mi>t</mi></msub><mi>&Delta;</mi><mi>&omega;</mi><mo>;</mo></mrow>]]></math><img file="FDA0001157247820000011.GIF" wi="846" he="119" /></maths><maths num="0002"><math><![CDATA[<mrow><mo>(</mo><msub><mi>&Delta;P</mi><mi>m</mi></msub><mo>-</mo><msub><mi>&Delta;P</mi><mi>e</mi></msub><mo>)</mo><mfrac><mn>1</mn><mrow><msub><mi>T</mi><mi>&sigma;</mi></msub><mi>s</mi></mrow></mfrac><mo>=</mo><mi>&Delta;</mi><mi>&omega;</mi><mo>;</mo></mrow>]]></math><img file="FDA0001157247820000012.GIF" wi="397" he="117" /></maths><maths num="0003"><math><![CDATA[<mrow><msub><mi>&Delta;P</mi><mi>e</mi></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>K</mi><mn>1</mn></msub><mo>-</mo><mfrac><mrow><msub><mi>K</mi><mn>2</mn></msub><msub><mi>K</mi><mn>3</mn></msub><mo>&lsqb;</mo><msub><mi>K</mi><mn>4</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>sT</mi><mi>R</mi></msub><mo>+</mo><msub><mi>K</mi><mn>5</mn></msub><msub><mi>K</mi><mi>A</mi></msub><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mrow><msup><mi>s</mi><mn>2</mn></msup><msub><mi>T</mi><mn>3</mn></msub><msub><mi>T</mi><mi>R</mi></msub><mo>+</mo><mi>s</mi><mrow><mo>(</mo><msub><mi>T</mi><mn>3</mn></msub><mo>+</mo><msub><mi>T</mi><mi>R</mi></msub><mo>)</mo></mrow><mo>+</mo><mn>1</mn><mo>+</mo><msub><mi>K</mi><mn>3</mn></msub><msub><mi>K</mi><mn>6</mn></msub><msub><mi>K</mi><mi>A</mi></msub></mrow></mfrac><mo>)</mo></mrow><mi>&Delta;</mi><mi>&delta;</mi><mo>;</mo></mrow>]]></math><img file="FDA0001157247820000013.GIF" wi="854" he="119" /></maths><maths num="0004"><math><![CDATA[<mrow><mo>&DoubleRightArrow;</mo><mfrac><mrow><msub><mi>&Delta;P</mi><mi>m</mi></msub></mrow><mrow><msub><mi>&Delta;&omega;</mi><mi>r</mi></msub></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mi>G</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>G</mi><mrow><mi>p</mi><mi>i</mi><mi>d</mi></mrow></msub><mo>)</mo></mrow><msub><mi>G</mi><mi>e</mi></msub><msub><mi>G</mi><mi>t</mi></msub></mrow><mrow><mn>1</mn><mo>+</mo><msub><mi>G</mi><mn>3</mn></msub><msub><mi>G</mi><mn>4</mn></msub></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0001157247820000014.GIF" wi="494" he="127" /></maths>上述公式中K<sub>1</sub>,K<sub>2</sub>,K<sub>3</sub>,K<sub>4</sub>,K<sub>5</sub>,K<sub>6</sub>分别为比例系数,K<sub>A</sub>为励磁机比例系数,s为经拉氏变换后的微分算子,δ为转速的调差系数,Δω为转速角速度偏量,Δω<sub>r</sub>为转速扰动值,Δδ为发电机功角偏差,ΔP<sub>m</sub>为机械功率增量,ΔP<sub>e</sub>为电磁功率增量,T<sub>1</sub>表示转速变送器时间常数,T<sub>3</sub>为励磁回路时间常数,T<sub>R</sub>为电压传感器时间常数,T<sub>σ</sub>为汽轮机转子的时间常数,G<sub>1</sub>为单机无穷大系统的传递函数,G<sub>2</sub>为控制系统中的转速变送器的一阶惯性传递函数,G<sub>e</sub>为电液伺服器的传递函数,G<sub>t</sub>为串联组合、单再热汽轮机的传递函数,G<sub>pid</sub>为控制器的传递函数,G<sub>3</sub>为汽轮机调速系统地传递函数,G<sub>4</sub>是为简化形式而建立的没有相应的物理含义;其中G<sub>1</sub>、G<sub>2</sub>、G<sub>3</sub>、G<sub>4</sub>定义为:<maths num="0005"><math><![CDATA[<mrow><msub><mi>G</mi><mn>1</mn></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>K</mi><mn>1</mn></msub><mo>-</mo><mfrac><mrow><msub><mi>K</mi><mn>2</mn></msub><msub><mi>K</mi><mn>3</mn></msub><mo>&lsqb;</mo><msub><mi>K</mi><mn>4</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>sT</mi><mi>R</mi></msub><mo>+</mo><msub><mi>K</mi><mn>5</mn></msub><msub><mi>K</mi><mi>A</mi></msub><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mrow><msup><mi>s</mi><mn>2</mn></msup><msub><mi>T</mi><mn>3</mn></msub><msub><mi>T</mi><mi>R</mi></msub><mo>+</mo><mi>s</mi><mrow><mo>(</mo><msub><mi>T</mi><mn>3</mn></msub><mo>+</mo><msub><mi>T</mi><mi>R</mi></msub><mo>)</mo></mrow><mo>+</mo><mn>1</mn><mo>+</mo><msub><mi>K</mi><mn>3</mn></msub><msub><mi>K</mi><mn>6</mn></msub><msub><mi>K</mi><mi>A</mi></msub></mrow></mfrac><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0001157247820000015.GIF" wi="772" he="111" /></maths><maths num="0006"><math><![CDATA[<mrow><msub><mi>G</mi><mn>2</mn></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>T</mi><mn>1</mn></msub><mi>s</mi><mo>)</mo></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0001157247820000016.GIF" wi="246" he="109" /></maths><maths num="0007"><math><![CDATA[<mrow><msub><mi>G</mi><mn>3</mn></msub><mo>=</mo><mrow><mo>(</mo><mfrac><mn>1</mn><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>T</mi><mn>1</mn></msub><mi>s</mi><mo>)</mo><mi>&delta;</mi></mrow></mfrac><mo>(</mo><mrow><mn>1</mn><mo>+</mo><msub><mi>G</mi><mrow><mi>p</mi><mi>i</mi><mi>d</mi></mrow></msub></mrow><mo>)</mo><mo>+</mo><mfrac><mn>1</mn><mi>s</mi></mfrac><msub><mi>G</mi><mn>1</mn></msub><msub><mi>G</mi><mrow><mi>p</mi><mi>i</mi><mi>d</mi></mrow></msub><mo>)</mo></mrow><msub><mi>G</mi><mi>e</mi></msub><msub><mi>G</mi><mi>t</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA0001157247820000017.GIF" wi="726" he="111" /></maths><maths num="0008"><math><![CDATA[<mrow><msub><mi>G</mi><mn>4</mn></msub><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>T</mi><mi>&sigma;</mi></msub><mi>s</mi><mo>+</mo><msub><mi>G</mi><mn>1</mn></msub><mfrac><mn>1</mn><mi>s</mi></mfrac></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0001157247820000018.GIF" wi="286" he="151" /></maths>步骤2:对控制器的传递函数<img file="FDA0001157247820000021.GIF" wi="322" he="115" />中的调速系统控制参数K<sub>p</sub>和K<sub>i</sub>在0~10范围内设置多组不同的数组,分别求出各组数组参数下被控对象的传递函数的特征根,通过求取数值方法求取特征根的数值解,进而求取调速系统的控制参数的取值范围,式中K<sub>i</sub>为调速系统积分控制参数、K<sub>p</sub>为比例控制参数、K<sub>d</sub>为微分控制参数,s为微分算子;分析被控对象的传递函数,得到被控对象传递函数的特征根可为实轴上的非零根以及共轭复数根,用下式表示:s=A+Bj,式中A和B分别为特征根的实部和虚部;对于二阶振荡环节,其特征根为一对共轭复数根,表达式如下所示:<maths num="0009"><math><![CDATA[<mrow><msub><mi>s</mi><mrow><mn>1</mn><mo>,</mo><mn>2</mn></mrow></msub><mo>=</mo><mo>-</mo><msub><mi>&xi;&omega;</mi><mi>n</mi></msub><mo>&PlusMinus;</mo><msub><mi>j&omega;</mi><mi>n</mi></msub><msqrt><mrow><mn>1</mn><mo>-</mo><msup><mi>&xi;</mi><mn>2</mn></msup></mrow></msqrt><mo>,</mo></mrow>]]></math><img file="FDA0001157247820000022.GIF" wi="595" he="86" /></maths>其中特征参数ξ为阻尼比,ω<sub>n</sub>为无阻尼振荡频率;步骤3:求出每组特征根的特征参数ξ的大小,并获得其中值最小的ξ,由上述公式可知,A=‑ξω<sub>n</sub>,<maths num="0010"><math><![CDATA[<mrow><mi>B</mi><mo>=</mo><mo>&PlusMinus;</mo><msub><mi>j&omega;</mi><mi>n</mi></msub><msqrt><mrow><mn>1</mn><mo>-</mo><msup><mi>&xi;</mi><mn>2</mn></msup></mrow></msqrt><mo>,</mo></mrow>]]></math><img file="FDA0001157247820000023.GIF" wi="446" he="79" /></maths>将上述方程中A、B当作已知,联立方程求得阻尼比<maths num="0011"><math><![CDATA[<mrow><mi>&xi;</mi><mo>=</mo><mfrac><mrow><mo>-</mo><mi>A</mi></mrow><msqrt><mrow><msup><mi>A</mi><mn>2</mn></msup><mo>+</mo><msup><mi>B</mi><mn>2</mn></msup></mrow></msqrt></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0001157247820000024.GIF" wi="350" he="150" /></maths>由上式知ξ小于零时,系统特征根实部不为负,二阶系统不稳定,二阶系统的阻尼比可以反映二阶系统的稳定与否,在分析考虑电网模型的汽轮机调速系统是否稳定时,对每个特征根求取其特征参数ξ;步骤4:根据步骤3中求得的调速系统中各组参数对应的特征根的最小特征参数ξ绘制特征参数ξ随调速器的参数K<sub>p</sub>、K<sub>i</sub>变化的曲线,当该特征根为共轭复根时,ξ为其二阶环节的阻尼比;当该特征根为非零实根时,ξ则取±1,若此实根小于零则ξ为1,若此实根大于零则ξ为‑1;当系统稳定时,则系统全部特征根均具有负实部,同时可以推出所有特征根的特征参数ξ均取正值,由此可根据被控对象所呈现的特征参数ξ随各组控制参数的变化图像确定使系统稳定的控制参数的取值范围;步骤5:根据前述步骤1至3所建立的汽轮机调速系统的数学模型以及步骤4绘制的变化趋势图,采用试凑法进行微调,整定出真正实用于机组的最优控制参数。
地址 100031 北京市西城区西长安街86号
您可能感兴趣的专利