基于激光吸收光譜的多用途氣體檢測系統設計

常岐海

（西藏民族學院，陜西 威陽 712082）

0 引 言

基于調制半導體激光吸收光譜 （tunable diode laser absorption spectrum，TDLAS）方法的氣體檢測技術，相對于傳統半導體、化學等傳感器，具有選擇性好、響應迅速和靈敏度高等特點，成為氣體檢測技術的主要研究發展方向之一。TDLAS技術不同于傳統傳感器之處在于可利用激光遙感形成非置入現場的自然目標短程遙感傳感器，避免高危氣體現場檢測人員的人身傷害[1]。

TDLAS方法利用發射波長調制激光穿越氣體，經光電探測器接收，分析接收光電信號變化測算激光照射路徑內的氣體含量。其中波長調制是實現高靈敏度檢測的重要環節。TDLAS波長調制檢測方法通過對激光進行波長掃描和光強調制，對接收光信號進行相敏檢測。通過提取光吸收信號的高次諧波幅度或線形作為反映氣體含量的特征。TDLAS技術利用氣體指紋吸收譜線進行氣體存在判定及測量，可使復雜背景氣體組份的干擾大大降低。該技術在使用前需要對檢測系統定標，即形成檢測譜線及相應的氣體含量數據，以便系統定量檢測[2-3]。

本文詳細介紹了TDALS方法的各種處理算法，并利用LabWindows/CVI開發了一套基于TDLAS方法的氣體檢測定標和測量的測試實驗系統。

1 TDLAS處理算法

圖 1（a）、（b）、（c）分別給出了 TDLAS 方法的氣體吸收光譜原始、一次諧波和二次諧波信號。由于TDLAS方法中對光譜信號的高次諧波信號的提取近似于對原始信號的求導。從原始調制吸收信號提取一次諧波信號就去除了線性本底，二次諧波信號的提取則不僅去除了線性本底噪聲，并去除了二次非線性本底噪聲。高次諧波信號隨階數的增加，信號幅度顯著下降[4]。根據探測靈敏度和信噪比，實際多選用二次諧波信號作為判定和測量氣體含量的特征信號。TDLAS處理算法可分為諧波信號最小二乘法和諧波比值法。

1.1 諧波信號最小二乘法

TDLAS方法中，可以將光電探測器接收的信號分為定標信號（Cal）、本底背景信號（Bgr）和現場檢測信號（Amb）。根據最小二乘法[5]，可以定義為

其中 Sn=f （Bgrn，Caln），N 為數據點個數。 函數 f一般認為是定標信號和本底背景信號的線性組合，如Sn=Bgr+c·Cal，其中c為待定系數。

現場檢測信號減去本底背景信號，然后參考定標信號，并使ξ最小化，求取的待定系數c即為現場檢測氣體的含量。根據基于最小二乘法的線性回歸使 ξ最小化，由式（1）可得

為了確定c，對式（2）求取ξ對于c的微分

1.2 諧波比值法

根據Beer-Lambert定律，調制激光輸出光強Ⅰ0與輸入光強Ⅰ和氣體濃度之間的關系[6]可表示為

式中：η——強度調制因子；

α（f）——氣體吸收系數；

f0——中心激光頻率；

fm——掃描頻率的幅度；

ω——電流調制的角頻率；

C——氣體濃度；

L——吸收路徑的長度。

在 α（f）CL＜＜1 和 η<<1 條件下，α（f）使用洛倫茨譜線線形，式（6）可簡化為

其中 x=fm/δf， δf為吸收譜線的半高半寬。 式（7）做傅里葉級數展開，可得一次諧波和二次諧波的表達式為 Ⅰ1=Ⅰ0η 和 Ⅰ2=-2kα0CLⅠ0，其中 k 為 x 的函數。一般選擇x=2.2，這樣可以使二次諧波的中心幅值最大，提高信噪比。二次諧波與一次諧波的比值為

圖1 TDLAS方法的氣體吸收光譜原始、一次諧波和二次諧波信號

該比值與吸收氣體含量成正比。式（8）比值消除了Ⅰ0，即消除本底噪聲的影響，提高了定量精度。此外，還有一種比值算法[7]，直接對氣體吸收信號做快速傅里葉變換，取二次諧波頻率對應的譜強度與一次諧波頻率對應的譜強度做比值，該比值在弱吸收條件下也與吸收氣體含量成正比。該算法是建立在吸收信號的諧波線形在任意頻率處的幅值都與氣體含量成正比基礎上。實際這也說明不需要一定選擇偶次諧波的中心頻率處幅值作為定標，選擇高階奇或偶次諧波非譜線頻率中心的其他給定頻率幅值一樣可以進行氣體定標和檢測。

上述兩種算法有其適用的范圍。對于最小二乘法，一般用于氣體吸收池或封閉空間檢測環境下的高精度痕量檢測。對于開放環境條件下，由于采用的線性回歸方法，其對背景噪聲的影響極其敏感，特別是在背景光多變的環境下，檢測量值會隨背景光的起伏而變化。為了在開放環境下使用最小二乘法，需要使用平衡探測接收消除背景光的起伏對該算法的影響。比值法相對于最小二乘法檢測精度低，其處理速度快，適用于快速現場檢測或高速移動檢測。比值法需對信號進行相關檢測，以判別是否是氣體的吸收信號。特別是利用比值法進行微弱氣體含量檢測，必須進行信號的相關性判別[8]。

2 多用途測試實驗系統

圖2給出了測試系統的原理框圖。該測試系統封裝在一個工控電腦機箱中，對外預留激光輸出接口和光電信號輸入接口。外接的氣體吸收池或短程遙測氣體檢測單元僅需通過這兩個接口即可互聯，方便地實現氣體檢測實驗系統。

圖2 測試系統原理框圖

圖3給出了測試系統的軟件原理框圖。系統定標流程與該原理框圖一致。系統根據光學接口參數的設置進行不同光譜信號處理算法的選擇。

軟件測試系統采用LabWindows開發，采用多線程實現。多線程分為采集線程、處理線程、顯示線程和數據存儲線程。線程之間使用線程安全隊列TSQ實現數據傳輸。LabWindows提供了兩種在次線程中運行代碼的高級機制，分別是線程池和異步定時器。本測試系統采用線程池方式運行次線程，使用指針的指針進行TSQ的數據讀寫。

圖3 測試系統軟件流程圖

2.1 多線程的實現

實驗系統的多線程實現通過調用Utility Lib庫函數CmtScheduleThreadPoolFunctionAdv，把想要運行的次線程的函數名傳遞給上述函數，線程池調用這個函數在它的其中一個線程中運行。

Visual C的多線程設計線程只是在線程內運行，而LabWindows的多線程設計線程分為兩部分，一部分維護線程運行，另一部分使用回調函數執行線程處理[9]。具體在LabWindows中的實現，通過處理線程實例說明如下：

上述處理線程函數中，首先放置一個線程間傳遞數據的隊列，CmtInstallTSQCallback函數用于在線程中安裝線程安全隊列（TSQ）。While語句判定線程是否運行，ProcessSystemEvents函數進行運行線程消息抓取。線程結束后順序使用CmtUninstallTSQCallback和CmtFlushTSQ函數完成TSQ的釋放。ProcessSystemEvents函數實現事件的響應，防止線程陷入線程死循環和造成用戶界面無法響應操作。否則程序盡管編譯通過，但運行會出錯。

2.2 TSQ指針的指針數據傳遞

LabWindows平臺給出的TSQ多設計成FIFO，每次采集的數據都放置在TSQ中[10]。對于大數據量的采集工作，這種數據結構安排是不符合結構化的數據管理。本實驗測試系統利用指針的指針，每次采集的數據通過在堆上分配空間，然后將指針傳遞給TSQ實現大量分批次的數據傳遞。這樣不僅可以保證程序設計的明晰，而且保證數據處理的連續性。以處理線程回調函數為例：

要實現指針的指針數據傳遞，TSQ函數指針必須設置成指針的指針。采集的數據必須掛在預先定義的指針數組中，才可有效實現。ProcessArrayPtr為指針數組，readPtr和writePtr分別為讀寫不同的TSQ指針。調試結果顯示，指針數組的存儲地址一致，有效地實現了線程間指針數據的傳遞。

3 結束語

本文詳細地論述了基于調制半導體激光吸收光譜方法的氣體檢測原理和不同應用現場環境下處理算法。對不同算法進行詳細比較和應用過程分析，給出了不同算法處理過程的適用性結論。面對多種檢測形式，設計了一種接口單一的調制半導體激光吸收光譜多用途氣體實驗測試系統。在該測試系統軟件中，基于LabWindows開發環境，提出使用指針的指針在線程安全隊列中進行優化的數據結構傳遞方法，并詳細地給出了一種多線程采集處理的線程執行結構。該調制半導體激光吸收光譜多用途氣體實驗測試系統完成后，已用于科研項目中的氣體標定和測試。

[1]Uehara H T K，Kimura K.Real-time monitoring of environmental methane and other gases with semiconductor lasers： a review[J].Sens.Actuators B，1997（38）：136-140.

[2]Yamanaka T.Analysis of gas sensor transient response by visualizing instantaneous gas concentration using smoke[J].Sens.Actuators A，1998（69）：77-81.

[3]Werle P.High sensitivity gas analysis by mid and near infrared diode laser[C]∥Proc.5th Int.Symposium on gas analysis by tunable diode laser，1998：1-16.

[4]Bomse D S， Stanton A C， Silver J A.Frequency modulation and wavelength modulation spectroscopies:comparison of experimental methods using a lead-salt diode laser[J].Appl.Opt.，1999（31）：718-731.

[5]Silver J A.Frequency-modulation spectroscopy for trace species detection:theory and comparison among experimental methods[J].Appl Opt，1992（31）：707-717.

[6]Werle P W，Mazzinghi P，Amatoa F D，et al.Signal processing and calibration procedures for in situ diodelaser absorption spectroscopy[J].Spectrochimica Acta Part A，2004（60）：1685-1705.

[7]Hodgkinson J，Well B，Padgett M，et al.Modelling and interpretation of gas detection using remote laser pointers[J].Spectrochimica Acta Part A，2006（63）：929-939.

[8]Jin W，Xu Y Z，Demokan M S，et al.Investigation of interferometric noise in fiber-optic gas sensors with use of wavelength modulation spectroscopy[J].Appl Opt，1997，36（28）：7239-7246.

[9]National Instruments corporation.LabWindows/CVI Basic course manual[Z].2006：202-209.

[10]National Instruments corporation.LabWindows/CVI user manual[Z].2006：153-160.