发明名称 危险废弃物热解焚烧炉内温度场分布非接触测量方法
摘要 本发明公开了一种危险废弃物热解焚烧炉内温度场分布非接触测量方法。它是利用焚烧炉内高温烟气的红外热辐射特性,首先对焚烧系统进行三维网格划分建立辐射投影方程,然后通过在不同位置安装红外辐射能传感器得到红外辐射投影图像,并从红外辐射能图像中提取不同波长下的辐射投影值,最后结合焚烧系统内三维温度分布的重建算法,计算危废焚烧系统中回转窑和二燃室内温度的分布数据。本发明提出的基于红外辐射能图像的温度场分布重建方法,是一种非接触式在线温度测量方法。它以危废焚烧系统内高温烟气的红外辐射传递特性为基础,通过采集的红外辐射能图像,在线重建焚烧系统内温度的三维分布。
申请公布号 CN101413826A 申请公布日期 2009.04.22
申请号 CN200810162524.5 申请日期 2008.12.01
申请人 浙江大学 发明人 严建华;黄群星;王飞;池涌;岑可法
分类号 G01J5/00(2006.01)I 主分类号 G01J5/00(2006.01)I
代理机构 杭州求是专利事务所有限公司 代理人 张法高
主权项 1.一种危险废弃物热解焚烧炉内温度场分布的非接触测量方法,其特征在于包括如下步骤:1)在回转窑(1)的头部(3)、回转窑(1)的尾部(4)和二燃室(2)的顶部(5)安装红外辐射能传感器,传感器响应波长范围1~7μm,以非接触方式获取红外辐射能图像;2)根据回转窑和二燃室结构尺寸进行X,Y,Z三个方向的网格划分,对所有网格进行编号为1…N;3)对传感器视场角内每条投影射线进行编号为1…M,并采用光线跟踪算法计算每条投影射线在每个网格内的穿越距离,建立如下的系数矩阵:<img file="A200810162524C00021.GIF" wi="541" he="341" />其中Δs<sub>i,k</sub>为第i条投影射线穿越第k个网格的距离;4)从红外辐射能传感器输出的红外辐射能图像中按照编号提取每条射线在两个不同红外波长λ<sub>1</sub>,λ<sub>2</sub>下的辐射能投影值,并构造向量:<maths num="0001"><![CDATA[<math><mrow><msub><mi>G</mi><msub><mi>&lambda;</mi><mn>1</mn></msub></msub><mo>=</mo><msup><mrow><mfenced open='[' close=']' separators=','><mtable><mtr><mtd><msub><mi>g</mi><mrow><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>,</mo><mn>1</mn></mrow></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>g</mi><mrow><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>,</mo><mi>i</mi></mrow></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>g</mi><mrow><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>,</mo><mi>M</mi></mrow></msub></mtd></mtr></mtable></mfenced></mrow><mi>T</mi></msup></mrow></math>]]></maths>和<maths num="0002"><![CDATA[<math><mrow><msub><mi>G</mi><msub><mi>&lambda;</mi><mn>2</mn></msub></msub><mo>=</mo><msup><mrow><mfenced open='[' close=']' separators=','><mtable><mtr><mtd><msub><mi>g</mi><mrow><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>,</mo><mn>1</mn></mrow></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>g</mi><mrow><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>,</mo><mi>i</mi></mrow></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>g</mi><mrow><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>,</mo><mi>M</mi></mrow></msub></mtd></mtr></mtable></mfenced></mrow><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths><img file="A200810162524C00024.GIF" wi="69" he="42" />和<img file="A200810162524C00025.GIF" wi="72" he="43" />为第i个条射线在波长λ<sub>1</sub>和λ<sub>2</sub>处的辐射投影强度;5)采用滤波反投影或正则化重建算法,求解如下的辐射投影方程,得到每个网格内两个波长下的红外发射辐射强度:G=H·S式中:<maths num="0003"><![CDATA[<math><mrow><msub><mi>G</mi><msub><mi>&lambda;</mi><mn>1</mn></msub></msub><mo>=</mo><msup><mrow><mfenced open='[' close=']' separators=','><mtable><mtr><mtd><msub><mi>g</mi><mrow><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>,</mo><mn>1</mn></mrow></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>g</mi><mrow><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>,</mo><mi>i</mi></mrow></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>g</mi><mrow><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>,</mo><mi>M</mi></mrow></msub></mtd></mtr></mtable></mfenced></mrow><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths>或者<maths num="0004"><![CDATA[<math><mrow><msub><mi>G</mi><msub><mi>&lambda;</mi><mn>2</mn></msub></msub><mo>=</mo><msup><mrow><mfenced open='[' close=']' separators=','><mtable><mtr><mtd><msub><mi>g</mi><mrow><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>,</mo><mn>1</mn></mrow></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>g</mi><mrow><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>,</mo><mi>i</mi></mrow></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>g</mi><mrow><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>,</mo><mi>M</mi></mrow></msub></mtd></mtr></mtable></mfenced></mrow><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths>S为网格穿越距离系数矩阵,<maths num="0005"><![CDATA[<math><mrow><msub><mi>H</mi><msub><mi>&lambda;</mi><mn>1</mn></msub></msub><mo>=</mo><msup><mrow><mfenced open='[' close=']' separators=','><mtable><mtr><mtd><msub><mi>H</mi><mrow><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>,</mo><mn>1</mn></mrow></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>H</mi><mrow><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>,</mo><mi>k</mi></mrow></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>H</mi><mrow><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>,</mo><mi>N</mi></mrow></msub></mtd></mtr></mtable></mfenced></mrow><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths>或者<maths num="0006"><![CDATA[<math><mrow><msub><mi>H</mi><msub><mi>&lambda;</mi><mn>2</mn></msub></msub><mo>=</mo><msup><mrow><mfenced open='[' close=']' separators=','><mtable><mtr><mtd><msub><mi>H</mi><mrow><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>,</mo><mn>1</mn></mrow></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>H</mi><mrow><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>,</mo><mi>k</mi></mrow></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>H</mi><mrow><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>,</mo><mi>N</mi></mrow></msub></mtd></mtr></mtable></mfenced></mrow><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths><img file="A200810162524C000210.GIF" wi="88" he="54" />和<img file="A200810162524C000211.GIF" wi="91" he="55" />为第k个网格内烟气在波长λ<sub>1</sub>和λ<sub>2</sub>处的发射辐射强度;6)根据每个网格内两个不同波长下的烟气发射辐射强度,采用如下的双波长温度计算公式得到系统各网格内的温度值:<maths num="0007"><![CDATA[<math><mrow><msub><mi>T</mi><mi>k</mi></msub><mo>=</mo><mfrac><mrow><mn>0.014388</mn><mrow><mo>(</mo><mfrac><mn>1</mn><msub><mi>&lambda;</mi><mn>2</mn></msub></mfrac><mo>-</mo><mfrac><mn>1</mn><msub><mi>&lambda;</mi><mn>1</mn></msub></mfrac><mo>)</mo></mrow></mrow><mrow><mi>ln</mi><mfrac><msub><mi>H</mi><mrow><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>,</mo><mi>k</mi></mrow></msub><msub><mi>H</mi><mrow><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>,</mo><mi>k</mi></mrow></msub></mfrac><mo>-</mo><mn>5</mn><mi>ln</mi><mfrac><msub><mi>&lambda;</mi><mn>2</mn></msub><msub><mi>&lambda;</mi><mn>1</mn></msub></mfrac></mrow></mfrac></mrow></math>]]></maths>式中T<sub>k</sub>为第k个网格内的温度值;7)采用OpenGL三维图形可视化技术,在线显示回转窑(1)和二燃室(2)内的所有网格内的温度值。
地址 310027浙江省杭州市浙大路38号