发明名称 基于测点加权值的燃气轮机燃烧系统在线监测方法
摘要 基于测点加权值的燃气轮机燃烧系统在线监测方法,属于燃气轮机燃烧系统监测领域。现有的燃烧监测系统难以对燃烧状态变化趋势做出判断的问题。一种基于测点加权值的燃气轮机燃烧系统在线监测方法,在燃气轮机的透平出口周向均匀地布置n个温度测点,得到t<sub>m</sub>时段内正常运行的排温数据T<sub>i</sub>;增加T<sub>i</sub>与T<sub>1</sub>的相关因子α<sub>i,1</sub>,根据T<sub>i</sub>与T<sub>1</sub>的关系函数,分别得到T<sub>i</sub>的预测值;计算无故障温度测点1排温理论值T<sub>1</sub>';令温度测点1的理论值T<sub>1</sub>'与温度测点1的实测值T<sub>1</sub>之差为△T<sub>1</sub>满足均值为0、标准差为σ<sub>1</sub>的正态分布;通过△T<sub>1</sub>与范围[-3σ<sub>1</sub>,3σ<sub>1</sub>]的关系进行温度测点的监测。本发明实现燃气轮机排温的在线监测,充分利用排温各个测点之间的相关性,准确检测出异常演变过程。
申请公布号 CN105134386A 申请公布日期 2015.12.09
申请号 CN201510556992.0 申请日期 2015.09.02
申请人 哈尔滨工业大学 发明人 刘金福;刘娇;万杰;刘晟;李飞;于达仁
分类号 F02C9/00(2006.01)I 主分类号 F02C9/00(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 杨立超
主权项 一种基于测点加权值的燃气轮机燃烧系统在线监测方法,其特征在于:所述基于测点加权值的燃气轮机燃烧系统在线监测方法通过以下步骤实现:步骤一、在燃气轮机的透平出口周向均匀地布置n个温度测点,燃气轮机在t<sub>m</sub>时段内正常运行,得到t<sub>m</sub>时段内的排温数据T=[T<sub>1</sub>,T<sub>2</sub>,…,T<sub>i</sub>,…T<sub>n</sub>];其中,T<sub>i</sub>表示第i个温度测点在t<sub>m</sub>时段内各时间点测得的排温数据,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>T</mi><mn>1</mn></msub><mo>=</mo><mo>&lsqb;</mo><msub><mi>T</mi><mrow><mn>1</mn><msub><mi>t</mi><mn>1</mn></msub></mrow></msub><mo>,</mo><msub><mi>T</mi><mrow><mn>1</mn><msub><mi>t</mi><mn>2</mn></msub></mrow></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>T</mi><mrow><mn>1</mn><msub><mi>t</mi><mi>m</mi></msub></mrow></msub><mo>&rsqb;</mo><mo>,</mo><msub><mi>T</mi><mn>2</mn></msub><mo>=</mo><mo>&lsqb;</mo><msub><mi>T</mi><mrow><mn>2</mn><msub><mi>t</mi><mn>1</mn></msub></mrow></msub><mo>,</mo><msub><mi>T</mi><mrow><mn>2</mn><msub><mi>t</mi><mn>2</mn></msub></mrow></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>T</mi><mrow><mn>2</mn><msub><mi>t</mi><mi>m</mi></msub></mrow></msub><mo>&rsqb;</mo><mo>,</mo></mrow>]]></math><img file="FDA0000795424410000013.GIF" wi="1010" he="108" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>T</mi><mi>i</mi></msub><mo>=</mo><mo>&lsqb;</mo><msub><mi>T</mi><mrow><msub><mi>it</mi><mn>1</mn></msub></mrow></msub><mo>,</mo><msub><mi>T</mi><mrow><msub><mi>it</mi><mn>2</mn></msub></mrow></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>T</mi><mrow><msub><mi>it</mi><mi>m</mi></msub></mrow></msub><mo>&rsqb;</mo><mo>,</mo><msub><mi>T</mi><mi>n</mi></msub><mo>=</mo><mo>&lsqb;</mo><msub><mi>T</mi><mrow><msub><mi>nt</mi><mn>1</mn></msub></mrow></msub><mo>,</mo><msub><mi>T</mi><mrow><msub><mi>nt</mi><mn>2</mn></msub></mrow></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>T</mi><mrow><msub><mi>nt</mi><mi>m</mi></msub></mrow></msub><mo>&rsqb;</mo><mo>;</mo></mrow>]]></math><img file="FDA0000795424410000014.GIF" wi="969" he="103" /></maths>步骤二、利用最小二乘法分别得到排温数据T<sub>i</sub>与排温数据T<sub>1</sub>的关系函数,即:T<sub>1</sub>=f<sub>i,1</sub>(T<sub>i</sub>),i=2,3,4,…,n其中,f<sub>i,1</sub>表示排温数据T<sub>i</sub>与排温数据T<sub>1</sub>的关系函数;满足:T<sub>1</sub>=k<sub>i,1</sub>T<sub>i</sub>+b<sub>i,1</sub>,i=2,3,4,...,n,其中k<sub>i,1</sub>和b<sub>i,1</sub>可由最小二乘法得出,即:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>k</mi><mrow><mi>i</mi><mo>,</mo><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>a</mi><mo>=</mo><msub><mi>t</mi><mn>1</mn></msub></mrow><msub><mi>t</mi><mi>m</mi></msub></munderover><mrow><mo>(</mo><msub><mi>T</mi><mrow><mi>i</mi><mi>a</mi></mrow></msub><mo>-</mo><mover><msub><mi>T</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>T</mi><mrow><mn>1</mn><mi>a</mi></mrow></msub><mo>-</mo><mover><msub><mi>T</mi><mn>1</mn></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow></mrow><mrow><munderover><mo>&Sigma;</mo><mrow><mi>a</mi><mo>=</mo><msub><mi>t</mi><mn>1</mn></msub></mrow><msub><mi>t</mi><mi>m</mi></msub></munderover><msup><mrow><mo>(</mo><msub><mi>T</mi><mrow><mi>i</mi><mi>a</mi></mrow></msub><mo>-</mo><mover><msub><mi>T</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac></mrow>]]></math><img file="FDA0000795424410000011.GIF" wi="541" he="292" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>b</mi><mrow><mi>i</mi><mo>,</mo><mn>1</mn></mrow></msub><mo>=</mo><mover><msub><mi>T</mi><mn>1</mn></msub><mo>&OverBar;</mo></mover><mo>-</mo><msub><mi>k</mi><mrow><mi>i</mi><mo>,</mo><mn>1</mn></mrow></msub><mover><msub><mi>T</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>;</mo></mrow>]]></math><img file="FDA0000795424410000012.GIF" wi="317" he="90" /></maths>步骤三、利用皮尔逊积矩系数表示排温数据T<sub>i</sub>与排温数据T<sub>1</sub>之间的相关性,得到排温数据T<sub>i</sub>与排温数据T<sub>1</sub>的相关因子α<sub>i,1</sub>,并得出:T<sub>1</sub>与T<sub>2</sub>间的相关因子为α<sub>2,1</sub>,T<sub>1</sub>与T<sub>3</sub>间的相关因子为α<sub>3,1</sub>,T<sub>1</sub>与T<sub>n</sub>间的相关因子为α<sub>n,1</sub>;步骤四:根据步骤二得到的得到排温数据T<sub>i</sub>与排温数据T<sub>1</sub>的关系函数,分别将t时刻测得的排温数据T<sub>2</sub>至T<sub>n</sub>带入到关系函数:T<sub>1</sub>=f<sub>i,1</sub>(T<sub>i</sub>)中,分别得到排温数据T<sub>1</sub>的预测值T<sub>2,1</sub>,T<sub>3,1</sub>,…T<sub>i,1</sub>…,T<sub>n,1</sub>:T<sub>i,1</sub>=f<sub>i,1</sub>(T<sub>i</sub>),i=2,3,4,…,n其中,f<sub>i,1</sub>表示步骤二得出的排温数据T<sub>i</sub>与排温数据T<sub>1</sub>的关系函数;t时刻是指t<sub>1</sub>到t<sub>m</sub>中的一个时刻;步骤五:计算t时刻燃气轮机无故障时温度测点1的排温理论值T<sub>1</sub>';步骤六:定义各时刻机组无故障时温度测点1的理论值T<sub>1</sub>'与温度测点1的实测值T<sub>1</sub>之差为ΔT<sub>1</sub>:ΔT<sub>1</sub>=T<sub>1</sub>'‑T<sub>1</sub>;将燃气轮机t<sub>m</sub>时段内的排温数据T=[T<sub>1</sub>,T<sub>2</sub>,…,T<sub>i</sub>,…T<sub>n</sub>]代入步骤四和步骤五,得到ΔT<sub>1</sub>满足均值为0、标准差为σ<sub>1</sub>的正态分布,即ΔT<sub>1</sub>~N(0,σ<sub>1</sub>);步骤七:对燃气轮机进入运行阶段待测温度测点1的监测:将每个时刻测得的排温数据代入步骤四和步骤五,得到若机组无故障运行时温度测点1的理论值T<sub>1</sub>';监测每个时刻机组无故障时温度测点1的理论值T<sub>1</sub>'与温度测点1的实测值T<sub>1</sub>之差为ΔT<sub>1</sub>,若ΔT<sub>1</sub>在[‑3σ<sub>1</sub>,3σ<sub>1</sub>]的范围内,则说明机组无故障,若超出[‑3σ<sub>1</sub>,3σ<sub>1</sub>]的范围,则说明机组发生故障;步骤八:对燃气轮机进入运行阶段待测温度测点2至n的监测:重复步骤二至步骤六,分别得到若机组无故障时温度测点2至n的理论值T<sub>2</sub>',T<sub>3</sub>',…,T<sub>n</sub>',以及理论值与实测值之差ΔT<sub>2</sub>,ΔT<sub>3</sub>,…,ΔT<sub>n</sub>的标准差Δσ<sub>2</sub>,Δσ<sub>3</sub>,…,Δσ<sub>n</sub>;相应地,若ΔT<sub>2</sub>,ΔT<sub>3</sub>,…,ΔT<sub>n</sub>都分别在[‑3σ<sub>n</sub>,3σ<sub>n</sub>]的范围,则说明机组无故障,若超出[‑3σ<sub>1</sub>,3σ<sub>1</sub>]的范围,则说明机组发生故障。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号