主权项 |
一种基于测点加权值的燃气轮机燃烧系统在线监测方法,其特征在于:所述基于测点加权值的燃气轮机燃烧系统在线监测方法通过以下步骤实现:步骤一、在燃气轮机的透平出口周向均匀地布置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>[</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>]</mo><mo>,</mo><msub><mi>T</mi><mn>2</mn></msub><mo>=</mo><mo>[</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>]</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>[</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>]</mo><mo>,</mo><msub><mi>T</mi><mi>n</mi></msub><mo>=</mo><mo>[</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>]</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>Σ</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>‾</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>‾</mo></mover><mo>)</mo></mrow></mrow><mrow><munderover><mo>Σ</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>‾</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>‾</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>‾</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>]的范围,则说明机组发生故障。 |