TDLAS技術用于燃燒場氣體溫度和濃度重建研究

殷可為，胥 頔，張 龍，楊富榮

?

TDLAS技術用于燃燒場氣體溫度和濃度重建研究

殷可為，胥 頔，張 龍，楊富榮

( 中國空氣動力研究與發展中心，四川綿陽 621000 )

基于可調諧半導體激光吸收光譜(Tunable Diode Laser Absorption Spectroscopy, TDLAS)技術及修正型代數迭代算法(Algebraic Reconstruction Technique, ART)，選擇頻率為7 153.7 cm-1和7 154.3 cm-1的H2O吸收線作為測溫譜線，測量了CH4/Air預混平面火焰爐燃燒場溫度和H2O氣濃度分布。仿真采用1 000 pixels×1 000 pixels個像素點描述方形模型區域，假設40條光束從兩個方向正交通過重建區域，通過引入初始分布并經修正型ART算法計算和三次多項式插值處理，溫度場重建偏差在4.5%以內，H2O氣濃度場重建偏差在4%以內。實驗采用導軌和轉臺實現探測光對穩態燃燒場兩個方向上的平行掃描，共獲取24路光束吸收信號，經修正型ART算法重建和三次多項式插值處理得到整個火焰爐燃燒區域的溫度場和H2O氣濃度場分布。將溫度場重建結果與熱電偶的測量結果進行對比，表明該方法能夠有效實現燃燒場溫度二維重建。

可調諧半導體激光吸收光譜；代數迭代算法；溫度場；濃度場；二維重建

0 引 言

燃燒流場參數測量是許多特種風洞試驗以及高超聲速推進系統研究的重要組成部分，需要對流場溫度、物質濃度、流場速度等參數開展定量測量，以便深刻理解特種風洞的復雜流動機理，優化高超聲速推進系統設計。因此，高精度的燃燒流場測量數據對于了解和預測一些特種風洞及高超聲速推進系統的性能至關重要。

如熱電偶等一系列的傳統接觸式測量設備，不僅響應速度慢、靈敏度低而且會擾亂燃燒場。此外，很多傳感器也存在無法長時間在高溫、高壓惡劣環境下工作的問題。而TDLAS技術可以進行氣體溫度、組分濃度、氣流流速的非接觸實時測量，設備重量輕，占用空間小，環境適應能力強，而且抗干擾能力強，測量數據精度高，完全擁有在線實時測量能力。然而，TDLAS的不足之處在于它是一種視線測量技術[1-4]，即測得的光譜吸收信號是沿測量路徑的積分值，反演的溫度和氣體濃度值則是該光路上的平均值，不能反映待測區域的場分布情況。

近年來國內外正興起基于TDLAS技術和計算機斷層掃描技術(Computed Tomography, CT)的場分布二維重建研究。它的計算過程是先將待測區域切分成很多網格，然后采用多條光路從不同方位同時測量待測區域，或者采用單條光路從不同方位掃描待測區域，經過重構算法反演，計算出每個網格內的待測值，從而實現待測區域二維高空間分辨率、高精度實時測量[5-7]。如今TDLAS傳感器已廣泛應用于航空發動機進氣道[8-9]、超然沖壓發動機隔離段/燃燒室[10-11]等流場診斷領域。在科學研究中有助于分析反應機理，對理論模型和數值模型進行驗證。在工程應用中有助于提高燃燒系統安全性，提升發動機性能。

本文基于TDLAS技術，采用波長掃描直接吸收的方式，通過修正型ART算法對燃燒場溫度和H2O氣濃度二維分布開展了仿真與實驗研究。在文獻[7]的基礎上，通過引入初始分布并改進迭代算法，大大提高了重建的網格數量，取得了更高分辨率的重建結果，有助于在風洞試驗中的應用。

1 測量原理

1.1 可調諧激光吸收光譜技術溫度濃度測量原理

TDLAS技術的理論基礎是Beer-Lambert定律，假設一束頻率為的激光束穿過待測區域，被均勻介質氣體吸收，光路上的積分吸光度可表示：

例如H2O分子的配分函數多項式系數如表1所示。

表1 H2O配分函數多項式系數

TDLAS雙線測溫法需要挑選同種氣體分子的兩條吸收線，將兩個吸收線強值相除消掉配分函數，得到不同吸收譜線強度的比值為

則待測區域溫度可以通過式(2)和式(4)表示為

在得到溫度參數后，一般選擇吸收較大的譜線來計算介質氣體組分濃度：

1.2 修正型ART重建算法

進行重建之前，將探測區域劃分成=×個網格，共有=2條光線相互垂直交叉穿過待測區域，如圖1所示。

圖1 待測區域離散網格圖

當路徑上介質氣體分布不均勻時，則式(1)變為

由激光準直器射出的第條光線在第個網格內的吸收長度L可由網格與待測區域的幾何關系得到。當網格被劃分得足夠小時，每一個網格內可視為具有相同的溫度和氣體濃度，光線經過待測區域內介質氣體的吸收，在接收端光強減弱，將式(7)離散化表示，則第條光線的吸收方程為

重建算法就是對式(8)進行代數迭代求解，其中A為第條光線穿過測量區域后的積分吸光度；為激光束編號；為網格編號；L代表第條光束穿過編號為網格的路徑長度；代表第網格內氣體吸收系數。式(8)可以進一步表示為

式中：為迭代次數；取為0.2。設第條光線計算出的積分吸光度與實際測量值的誤差e為

圖2 修正型ART算法流程圖

2 數值仿真

譜線選擇的基本原則可以概括為：譜線為1.3 μm至1.5 μm的近紅外光譜區；根據被測對象，譜線對有較強的吸收從而保證有較高的信噪比；盡量減少臨近譜線的干擾。本文選擇的兩條H2O吸收線頻率分別為7 153.7 cm-1和7 154.3 cm-1，計算數據來自于Hitran數據庫，兩條譜線線強比及測溫靈敏度隨溫度的變化關系如圖3所示。

圖3 譜線線強比及測溫靈敏度隨溫度的變化關系

由圖3可見，線強比和測溫靈敏度都隨著溫度的升高而降低，但溫度在1 000 K時，兩個參數的值均大于1，足以在平面火焰爐上用于研究燃燒場的溫度及H2O氣濃度分布。

假設待測區域為10 cm×10 cm的正方形區域，用1 000 pixels×1 000 pixels個像素點描述，在此區域內H2O氣組分濃度均勻分布(=0.05)，燃燒流場中溫度場分布如圖4所示。另外，在溫度場和濃度場重建前，將待測區域劃分為20×20個網格，并按照圖1所示的方式布置光路，總共有40條光束相互垂直交叉穿過待測區域。

根據修正型ART算法的流程圖，預先給定混合氣體總壓=1 atm，迭代步長=0.2。對于圖4所示的待測溫度場模型，如果將初始溫度場分布設為一常數(1 000 K)，而不對初始溫度場進行預估，經修正型ART算法計算及三次多項式插值處理后，重建的溫度場分布如圖5所示，可見重建結果已經嚴重失真。

圖4 待重建的溫度場模型

圖5 重建后的溫度場分布

出現這種情況的主要原因是投影光線數目遠遠小于未知數的數目，非齊次線性方程組式(9)屬于有無窮多解的情況，因此計算機會根據不同的初始分布計算出不同的結果。有兩種方法可以解決這個問題：一是增加光線數目和不同的光線投影角度，文獻[5-6]都對這個問題進行了討論；二是引入預估的初始分布進行干預。該方法并不要求初始分布能達到多高的精度，但要求初始分布具有與實際分布比較接近的外形。在燃燒流場測量實驗中，往往空間狹小且設備振動強烈，布置的光線數量極為有限，便可采用第二種方案，使用計算流體力學軟件獲得待測參數的初始分布情況[12-13]。

對于圖4所示的待測溫度場，設置了圖6所示的初始溫度場分布(預計峰值溫度是900 K，剩余區域溫度是600 K)和濃度場分布(=0.1)，從而可計算得到各個網格內的初始吸收系數的值。

圖6 預估的初始溫度場分布

圖7 重建結果及溫度絕對誤差分布

由圖7可見，這種方法對雙峰溫度場的重建效果較好，溫度的最大絕對誤差為40 K，最大相對誤差小于4.5%。H2O氣濃度的最大絕對誤差小于0.002，最大相對誤差小于4%。

3 測量實驗

實驗裝置如圖8所示，先將平面火焰爐放置在固定的導軌與轉臺之上，平面火焰爐兩側分別固定激光準直器和紅外探測器，對爐面不同位置吸收光譜的測量是通過精確移動和旋轉火焰爐來實現。光線布置方式仍然如圖1所示。首先點燃火焰爐，通過約20 min的燃燒，待火焰穩定后，在方向上測量了12路信號，每路信號之間移動了0.5 cm。信號發生器的掃描頻率設為400 Hz，每掃描200次輸出一次信號的平均值，每路信號的測量時間約為30 s，每條光線的積分吸光度取時間上的平均值。然后將平面火焰爐旋轉了90°，以同樣的方式在方向上也測量了12路信號，共計有24路信號相互垂直交叉穿過待測區域。

圖8 實驗裝置

實驗采用NEL公司NLK1E5EAAA可調諧激光二極管，中心波長在1 397.8 nm附近，一次掃描可獲得H2O分子的兩條吸收譜線(頻率為7 153.7 cm-1和7 154.3 cm-1)的信號。使用的平面火焰爐直徑為6 cm，燃燒氣體是甲烷和空氣的混合氣體，通過精確控制氣體流量，甲烷和空氣的配氣維持在1.5:15。光線離爐面高度為0.5 cm。

實驗得到吸光度信號如圖9所示，為了精確計算兩條譜線積分吸光度的值，需要進一步對吸光度信號進行voigt線型擬合[14]，并給出了擬合之后的殘差。

圖9 吸光度信號和voigt線型擬合結果

24路信號都會獲得上述吸光度曲線和voigt線型擬合曲線，根據擬合的結果可分別算得兩條譜線的積分吸光度。將48個積分吸光度的值導入修正型ART算法程序，初始溫度分布設為溫度1 000 K，組分濃度5%，半徑3 cm的圓。通過計算得到各個網格內的溫度值和氣體組分濃度值，如圖10所示。

再利用三次多項式插值法對各網格之間的溫度和氣體組分濃度進行插值計算，便可得到高分辨率的溫度場分布和濃度場分布，如圖11所示。

從圖11的結果可以看出，平面爐燃燒區域的最高溫度為1 204 K，且位于爐面的左下方，與理論情況不符(理論上最高溫度應位于爐面正中心)，其原因是平面爐的排氣孔排氣不均勻；另外，平面火焰爐在移動時也會產生誤差，對重建結果造成一定的影響。燃燒區域周圍的溫度在700 K左右，低于燃燒中心區域的溫度。H2O氣濃度受外界空氣濕度的影響(實驗當天的空氣濕度在30%～40%)，導致爐面一些邊緣區域的濃度達到8%~9%左右，高于中心區域5%的H2O氣濃度。

圖10 溫度場和濃度場重建結果

此外，采用熱電偶，沿著軸測量了溫度的分布情況。熱電偶的測溫誤差為5%，重建的結果和熱電偶測量的結果如圖12所示。

圖12 重建結果與熱電偶測量結果比較

由圖12可見，重建結果和熱電偶的測量結果存在一定的偏差，其中在火焰爐邊緣偏差較大，而火焰爐中心區域偏差較小。原因是在實驗進行的過程中，采用單光路掃描測量方式，測量過程并不是同時進行，因此測量結果會受到火焰不穩定性的影響，而爐面邊緣的火焰穩定性較差；另外，重建結果還會受到平面爐位移誤差的影響。

4 結 論

本文從理論上研究了將可調諧半導體激光吸收光譜技術應用于燃燒場溫度和氣體濃度二維分布重建，并在此基礎上編寫了修正型ART算法程序，通過引入初始分布，大大提高了重建質量。通過仿真表明該程序具有較高的重建精度和分辨率。另外，進行了實驗研究。通過火焰爐在導軌上的平移和在轉臺上的旋轉，獲得了24路光束的H2O氣吸收光譜信號，再使用該算法實現了對待測區域溫度場和H2O氣濃度場的重建。通過與熱電偶的測量結果的比較，發現它們趨勢基本一致。證明該重建方法能夠實現燃燒場溫度二維分布測量。下一步工作將搭建多光路測量系統并將該技術用于風洞試驗。

[1] ZHOU Xin，LIU Xiang，Jay B Jeffries，. Development of a sensor for temperature and water concentration in combustion gases using a single tunable diode laser [J].Measurement Science &Technology(S0957-0233)，2003，14：1459-1468.

[2] 李飛，余西龍，陳立紅，等. TDLAS測量甲烷/空氣預混平面火焰溫度和H2O濃度[J]. 實驗流體力學，2009，23(2)：40-44.

LI Fei，YU Xilong，CHEN Lihong，. Temperature and water vapour concentration measurements of CH4/air premixed flat flame based on TDLAS [J]. Journal of Experiments in Fluid Mechanics，2009，23(2)：40-44.

[3] 許振宇，劉文清，劉建國，等. 基于可調諧半導體激光器吸收光譜的溫度測量方法研究[J]. 物理學報，2012，61(23)：234204.

XU Zhenyu，LIU Wenqing，LIU Jiangou，. Temperature measurements based on tunable diode laser absorption spectroscopy [J]. Acta Physics Sinica，2012，61(23)：234204.

[4] 孫鵬帥，張志榮，崔小娟，等. 燃燒場內多路徑溫度與H2O濃度的在線檢測[J]. 中國激光，2015，42(9)：0915002.

SUN Pengshuai，ZHANG Zhirong，CUI Xiaojuan，. Multipath Real-Time measurement of Temperature and H2O Concentration for Combustion Diagnosis [J]. Chinese Journal of Lasers，2015，42(9)：0915002.

[5] 宋俊玲，洪延姬，王廣宇，等. 基于激光吸收光譜技術的燃燒場氣體溫度和濃度二維分布重建研究[J]. 物理學報，2012，61(24)：240702.

SONG Junling，HONG Yanji，WANG Guangyu，. Two-dimensional reconstructions of gas temperature and concentration in combustion based on tunable diode laser absorption spectroscopy [J]. Acta Physics Sinica，2012，61(24)：240702.

[6] 李金義，朱可柯，杜振輝，等. 一種基于TDLAS的高分辨率二維溫度場重建算法及數值仿真[J]. 實驗流體力學，2014，28(3)：63-71.

LI Jinyi，ZHU Keke，DU Zhenhui，. A high-resolution algorithm for 2D temperature reconstruction using TDLAS and numerical simulation [J]. Journal of Experiments in Fluid Mechanics，2014，28(3)：63-71.

[7] 夏暉暉，劉建國，許振宇，等. 基于代數迭代算法的燃燒火焰溫度場和氣體濃度場重建研究[J]. 光譜學與光譜分析，2015，35(10)：2697-2702.

XIA Huihui，LIU Jianguo，XU Zhenyu，. Reconstruction Research for Gas Concentration and Temperature of Flame Based on Algebraic Reconstruction Technique [J]. Spectroscopy and Spectral Analysis，2015，35(10)：2697-2702.

[8] Miller M F，Kessler W J，Allen M G. Diode laser-based air mass flux sensor for subsonic aero-propulsion inlets [J]. Applied Optics(S0003-6935)，1996，35(24)：4905-4912.

[9] Lyle K H，Jeffries J B，Hanson R K. Diode laser sensor for air mass flux based on oxygen absorption 2：non-uniform flow modeling and aero-engine tests [J]. AIAA Journal(S0001-1452)，2007，45(9)：2213-2223.

[10] Chang L S，Strand C L，Jeffries J B，. Supersonic massflux measurements via tunable diode laser absorption and nonuniform flow modeling [J]. AIAA Journal(S0001-1452)，2011，49(12)：2783-2791.

[11] Rieker G B，Jeffries J B，Hanson R K. Calibration-free wavelength modulation spectroscopy for measurements of gas temperature and concentration in harsh environments [J]. Applied Optics(S0003-6935)，2009，48(29)：5546-5560.

[12] Lindstrom C D，Jackson K R，Williams S，. Shock-Train Structure Resolved with Absorption Spectroscopy Part I System Design and Validation [J]. AIAA Journal(S0001-1452)，2009，47(10)：2368-2378.

[13] Lindstrom C D，Davis D，Williams S，. Shock-Train Structure Resolved with Absorption Spectroscopy Part II Analysis and CFD Comparison [J]. AIAA Journal(S0001-1452)，2009，47(10)：2379-2390.

[14] 尹增謙，武臣，宮琬鈺，等. Voigt線型函數及其最大值的研究[J]. 物理學報，2013，62(12)：123301.

YIN Zengqian，WU Chen，GONG Wanyu，. Voigt profile function and its maximum [J]. Acta Physics Sinica，2013，62(12)：123301.

2D Reconstruction for Gas Temperature and Concentration Based on TDLAS

YIN Kewei，XU Di，ZHANG Long，YANG Furong

( China Aerodynamics Research & Development Center, Mianyang 621000, Sichuan, China)

7 153.7cm-1and 7 154.3cm-1absorption light spectrums of H2O vapor were chosen to measure temperature distributions and H2O vapor concentration distributions of the combustion region of CH4/Air flat flame furnace based on Tunable Diode Absorption Spectroscopy (TDLAS) and modified Algebraic Reconstruction Technique (ART). In the numerical simulation, combustion region was described by 1 000 pixels×1 000 pixels. There were 40 laser rays cross the combustion region in orthogonal two directions. By introducing the initial distribution, we found the temperature field reconstruction relative error was less than 4.5% and H2O vapor concentration field reconstruction relative error was less than 4% based on modified ART algorithm and cubic polynomial interpolation processing. In the experiment, 24 laser rays were used to scan the target region of CH4/Air flat flame furnace in orthogonal two directions with rotary table and slide way. Temperature field and H2O vapor concentration field were obtained based on modified ART algorithm and cubic polynomial interpolation. A comparison of temperature reconstruction result and thermo couple result indicates the method is effective to achieve temperature field reconstruction of the combustion region.

tunable diode laser absorption spectroscopy; algebraic reconstruction technique; temperature field; concentration field; two dimensional reconstruction

1003-501X(2016)12-0020-08

O433.1

A

10.3969/j.issn.1003-501X.2016.12.004

2016-09-07；

2016-11-05

殷可為(1986-)，男(漢族)，重慶人。工程師，博士，主要研究工作是風洞非接觸測量技術。E-mail：ykw.258@163.com。