发明名称 一种燃气轮机燃烧室振源位置的确定方法
摘要 一种燃气轮机燃烧室振源位置的确定方法,该方法包括:在燃烧室多个位置安装多路动态压力传感器;多路传感器与采集系统的不同数据采集通道相连;各路传感器同时测量设定时长的动态压力信号并通过数据采集通道同步记录和存储信号;将各路信号进行傅里叶展开,分解为不同频率的谐波信号叠加;以某路传感器测量的信号中含有的某阶频率的谐波信号为基准,分别与其它几路传感器测量信号中含有的同阶频率的谐波信号做相关分析,然后记录各相关函数中第一次出现尖峰所对应的时刻,这些时刻即为某阶频率的谐波信号的振源点传到各路传感器的时差;最后根据波的传播原理确定出不同频率的谐波信号的振源点,这些振源点的集合即为振源区域在燃烧室内的位置。
申请公布号 CN103389213A 申请公布日期 2013.11.13
申请号 CN201310339976.7 申请日期 2013.08.06
申请人 北京华清燃气轮机与煤气化联合循环工程技术有限公司 发明人 谢法;井文明
分类号 G01M15/14(2006.01)I 主分类号 G01M15/14(2006.01)I
代理机构 北京鸿元知识产权代理有限公司 11327 代理人 邸更岩
主权项 1.一种燃气轮机燃烧室振源位置的确定方法,其特征是所述方法包括如下步骤:1)沿气流方向,在燃烧室下游不同横截面上安装N路动态压力传感器,其中N大于等于5;每路动态压力传感器的受感部位都深入到燃烧室内;N路动态压力传感器分别与采集系统的不同数据采集通道相连;2)N路动态压力传感器同时测量设定时长t的动态压力信号,并通过不同的数据采集通道同步记录和存储动态压力信号;3)将每路动态压力信号进行傅里叶展开,分解为不同频率的谐波信号叠加,即:a)将各路动态压力传感器测量的动态压力信号分别记为S<sub>1</sub>(t),S<sub>2</sub>(t)、……,S<sub>j</sub>(t),S<sub>j+1</sub>(t)……S<sub>N</sub>(t),t∈[τ<sub>1</sub>,τ<sub>2</sub>];S<sub>j</sub>(t)为第j路动态压力信号,N为传感器个数,τ<sub>1</sub>为开始记录数据的时刻,τ<sub>2</sub>为结束记录数据的时刻;b)将各路动态压力信号进行傅里叶展开,分解为不同频率的谐波信号叠加,其傅里叶展开的表达式为:<maths num="0001"><![CDATA[<math><mrow><msub><mi>S</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>i</mi></munderover><msub><mi>G</mi><mi>jn</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>n</mi></msub><mo>&CenterDot;</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>t∈[τ<sub>1</sub>,τ<sub>2</sub>],j∈[1,N],n∈[1,i],其中,n为频率阶数;f<sub>n</sub>为第n阶频率;G<sub>jn</sub>(f<sub>n</sub>·t)为第j路动态压力传感器测量的信号中含有的第n阶频率谐波信号;i为动态压力传感器测量的信号傅里叶展开所取的最高阶数,i取2~10中的任意一个整数值;4)以某路传感器测量的信号中含有的某阶频率的谐波信号为基准,分别与其它几路传感器测量的信号中含有的同阶频率的谐波信号做相关函数分析,然后记录各相关函数中第一次出现尖峰所对应的时刻,这些时刻即为某阶频率的谐波信号的振源点传到各路动态压力传感器的时差,即:a)以某路动态压力信号S<sub>j</sub>(t)中含有的第n阶频率的谐波信号G<sub>jn</sub>(f<sub>n</sub>·t)为基准,分别与其它几路动态压力信号中含有的第n阶频率的谐波信号做相关函数分析,从而获得各路动态压力传感器测量的信号中含有的第n阶频率的谐波信号之间的相关函数<img file="FDA00003628160900012.GIF" wi="164" he="84" /><img file="FDA00003628160900013.GIF" wi="158" he="81" />……,<img file="FDA00003628160900014.GIF" wi="166" he="80" />……,,<img file="FDA00003628160900015.GIF" wi="172" he="80" />相关函数的表达式为:<maths num="0002"><![CDATA[<math><mrow><msub><mi>R</mi><mrow><msub><mi>G</mi><mi>jn</mi></msub><msub><mi>G</mi><mi>kn</mi></msub></mrow></msub><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><msub><mi>&tau;</mi><mn>1</mn></msub></mrow><msub><mi>&tau;</mi><mn>2</mn></msub></munderover><msub><mi>G</mi><mi>jn</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>&CenterDot;</mo><mi>t</mi><mo>)</mo></mrow><msub><mi>G</mi><mi>kn</mi></msub><mo>[</mo><msub><mi>f</mi><mi>n</mi></msub><mo>&CenterDot;</mo><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mi>m</mi><mo>)</mo></mrow><mo>]</mo></mrow></math>]]></maths>k∈[1,N]且k≠j,m∈[0,τ<sub>2</sub>-τ<sub>1</sub>],n∈[1,i]b)分别记录各相关函数<img file="FDA00003628160900023.GIF" wi="148" he="74" /><img file="FDA00003628160900024.GIF" wi="144" he="75" />……,<img file="FDA00003628160900025.GIF" wi="154" he="73" />……,<img file="FDA00003628160900026.GIF" wi="132" he="80" />中首次出现尖峰所对应的时刻,该时刻即为对应的第n阶频率的谐波信号的振源点传到各个动态压力传感器的时差,分别记为Δt<sub>j1,n</sub>,Δt<sub>j2,n</sub>,……,Δt<sub>jk,n</sub>,……,Δt<sub>jN,n</sub>;Δt<sub>jk,n</sub>;5)确定振源区域所在的具体位置a)以第j路动态压力传感器受感部中心位置为坐标原点建立空间三维笛卡尔坐标系,设第k路动态压力传感器与第j路动态压力传感器在三维坐标轴方向的距离为a<sub>jk</sub>,b<sub>jk</sub>,c<sub>jk</sub>;第n阶频率的谐波信号的振源点的坐标为P<sub>n</sub>(x<sub>n</sub>,y<sub>n</sub>,z<sub>n</sub>),传播的平均传播速度为C<sub>n</sub>;b)根据波的传播原理,可得到如下方程组:<maths num="0003"><![CDATA[<math><mrow><mfenced open='' close='}'><mtable><mtr><mtd><msqrt><msubsup><mi>x</mi><mi>n</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>y</mi><mi>n</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>Z</mi><mi>n</mi><mn>2</mn></msubsup></msqrt><mo>-</mo><msqrt><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>n</mi></msub><mo>-</mo><msub><mi>a</mi><mrow><mi>j</mi><mn>1</mn></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>n</mi></msub><mo>-</mo><msub><mi>b</mi><mrow><mi>j</mi><mn>1</mn></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>z</mi><mi>n</mi></msub><mo>-</mo><msub><mi>c</mi><mrow><mi>j</mi><mn>1</mn></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>-</mo><msub><mi>C</mi><mi>n</mi></msub><msub><mi>&Delta;t</mi><mrow><mi>j</mi><mn>1</mn><mo>,</mo><mi>n</mi></mrow></msub><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><msqrt><msubsup><mi>x</mi><mi>n</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>y</mi><mi>n</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>z</mi><mi>n</mi><mn>2</mn></msubsup></msqrt><mo>-</mo><msqrt><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>n</mi></msub><mo>-</mo><msub><mi>a</mi><mrow><mi>j</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>n</mi></msub><mo>-</mo><msub><mi>b</mi><mrow><mi>j</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>z</mi><mi>n</mi></msub><mo>-</mo><msub><mi>c</mi><mrow><mi>j</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>-</mo><msub><mi>C</mi><mi>n</mi></msub><msub><mi>&Delta;t</mi><mrow><mi>j</mi><mn>2</mn><mo>,</mo><mi>n</mi></mrow></msub><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msqrt><msubsup><mi>x</mi><mi>n</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>y</mi><mi>n</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>z</mi><mi>n</mi><mn>2</mn></msubsup></msqrt><mo>-</mo><msqrt><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>n</mi></msub><mo>-</mo><msub><mi>a</mi><mi>jk</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>n</mi></msub><mo>-</mo><msub><mi>b</mi><mi>jk</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>z</mi><mi>n</mi></msub><mo>-</mo><msub><mi>c</mi><mi>jk</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>-</mo><msub><mi>C</mi><mi>n</mi></msub><msub><mi>&Delta;t</mi><mrow><mi>jk</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msqrt><msubsup><mi>x</mi><mi>n</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>y</mi><mi>n</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>z</mi><mi>n</mi><mn>2</mn></msubsup></msqrt><mo>-</mo><msqrt><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>n</mi></msub><mo>-</mo><msub><mi>a</mi><mi>jN</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>n</mi></msub><mo>-</mo><msub><mi>b</mi><mi>jN</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>z</mi><mi>n</mi></msub><mo>-</mo><msub><mi>c</mi><mi>jN</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>-</mo><msub><mi>C</mi><mi>n</mi></msub><msub><mi>&Delta;t</mi><mrow><mi>jN</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>=</mo><mn>0</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>方程中,k∈[1,N]且k≠j;a<sub>jk</sub>,b<sub>jk</sub>,c<sub>jk</sub>是已知的;c)求解方程组(1)所得的x<sub>n</sub>,y<sub>n</sub>,z<sub>n</sub>的值即为第n阶频率的谐波信号的振源点P<sub>n</sub>的坐标值,n∈[1,i]在空间三维笛卡尔坐标系的位置,这些由不同阶频率的谐波信号确定的振源点的集合即为振源区域在燃烧室内的位置。
地址 100084 北京市海淀区中关村东路1号院8号楼1001号