李國(guó)林，焦 月，馬 坤，鄢志丹，張雪娜，吳赟輝，張澤成

（中國(guó)石油大學(xué)（華東）控制科學(xué)與工程學(xué)院，山東青島266580）

1 引 言

煤、石油、天然氣等化石燃料的燃燒、動(dòng)植物的呼吸作用都會(huì)造成大量CO2的排放，而CO2是溫室效應(yīng)的“罪魁禍?zhǔn)住敝唬?］，因此，檢測(cè)空氣中CO2的含量對(duì)于監(jiān)測(cè)世界氣候的變化具有重要意義。檢測(cè)呼吸氣體中CO2的含量是識(shí)別肺動(dòng)脈栓塞的有效手段之一，中高檔監(jiān)護(hù)儀已將呼吸末CO2列為必備檢測(cè)功能［2］。此外，烯烴裂解爐會(huì)因?yàn)榭s合反應(yīng)產(chǎn)生沉積的焦炭顆粒，影響裂解效率，需要定期停爐清焦，通過(guò)分析燒焦過(guò)程的CO2含量監(jiān)測(cè)清焦進(jìn)展，可以快速反應(yīng)燒焦進(jìn)度，從而精準(zhǔn)地控制燒焦時(shí)間，提升工藝生產(chǎn)的經(jīng)濟(jì)效益［3］。

目前，檢測(cè)CO2的常用方法有電化學(xué)法［4-5，10］、非分光紅外光譜技術(shù)［5-6，10］以及可調(diào)諧半導(dǎo)體激光吸收光譜技術(shù)（Tunable Diode Laser Absorption Spectroscopy， TDLAS）［7-15］等 。Schwandt 等人用固態(tài)電化學(xué)傳感器檢測(cè)CO2，這種傳感器壽命短且響應(yīng)時(shí)間長(zhǎng)［4-5，10］。Vincent 等用非分光紅外光譜技術(shù)檢測(cè)呼吸氣體中的CO2，系統(tǒng)分辨率約為10×10-6，該方法容易受到光源強(qiáng)度變化和背景漂移帶來(lái)的影響，需要時(shí)常校準(zhǔn)頻率以及更換不同的濾光片［5-6，10］。TDLAS 方法常用空芯光纖和多反氣室，Xiong 等人用基于空芯光纖的CO2傳感器監(jiān)測(cè)呼吸過(guò)程中CO2濃度，相比于傳統(tǒng)的多反氣室，空芯光纖更容易產(chǎn)生干涉條紋及受到外界環(huán)境的干擾，Li 等用多反氣室和基于級(jí)聯(lián)積分器梳狀濾波器的檢測(cè)系統(tǒng)檢測(cè)CO2在1 580 nm 附近的濃度，靈敏度高、檢出限低、響應(yīng)速度快［7-15］。綜上可知，電化學(xué)傳感器通常含有有毒有害材料，壽命短，重復(fù)性差且響應(yīng)時(shí)間長(zhǎng)；基于濾光片的非分光紅外光譜技術(shù)熱光源的光譜寬，容易涵蓋多種氣體，干擾多，容易受到光源強(qiáng)度變化和背景漂移帶來(lái)的影響，需要時(shí)常對(duì)頻率進(jìn)行校準(zhǔn)以及使用不同的濾光片，氣體吸收光程不夠長(zhǎng)，檢測(cè)精度不如TDLAS 高，響應(yīng)速度慢，不適用于工業(yè)過(guò)程中的實(shí)時(shí)在線分析。TDLAS 不易受到外界環(huán)境干擾，不需要參比氣室，結(jié)構(gòu)簡(jiǎn)單穩(wěn)定；響應(yīng)速度快，可實(shí)時(shí)在線分析；基于多反的長(zhǎng)光程氣室光譜分辨率高，靈敏度高［4-15］。因此，TDLAS 在工業(yè)應(yīng)用中更具潛力。

機(jī)器學(xué)習(xí)算法被廣泛應(yīng)用于氣體濃度預(yù)測(cè)中。Chen 等用BP 神經(jīng)網(wǎng)絡(luò)建立了基于衛(wèi)星遙感數(shù)據(jù)的區(qū)域尺度污染物濃度預(yù)測(cè)模型，BP 神經(jīng)網(wǎng)絡(luò)具有良好的非線性擬合能力，但訓(xùn)練速度慢，局部最優(yōu)下降，不能滿足實(shí)時(shí)性和高精度的要求［16］。遺傳算法（Genetic Algorithm，GA）是一種解決最優(yōu)化的搜索啟發(fā)式算法，Liu 等利用遺傳算法優(yōu)化后的BP 神經(jīng)網(wǎng)絡(luò)檢測(cè)瓦斯的主要成分甲烷，設(shè)計(jì)了燃?xì)庑孤╊A(yù)警系統(tǒng)［17］，但是其相關(guān)系數(shù)最高僅為0.92。Sun 等將粒子群尋優(yōu)（Particle Swarm Optimization，PSO）算法與BP結(jié)合，對(duì)河北省二氧化碳排放的影響因素進(jìn)行了研究，結(jié)構(gòu)簡(jiǎn)單，但是計(jì)算速度慢［18］。Yang 等提出了一種基于LSSVM 的氮氧化物濃度實(shí)時(shí)動(dòng)態(tài)預(yù)測(cè)模型，但是動(dòng)態(tài)時(shí)間長(zhǎng)［19］。Li 等用PSO 與支持向量回歸（Support Vector Regression，SVR）結(jié)合進(jìn)行大氣污染物濃度的預(yù)測(cè)，其預(yù)測(cè)精度高于GA-SVR［20］。黃廣斌在極限學(xué)習(xí)機(jī)（Extreme Learning Machine，ELM）的基礎(chǔ)上引入核函數(shù)，提出核極限學(xué)習(xí)機(jī)（Kernel Extreme Learning Machine，KELM），簡(jiǎn)化了ELM 的參數(shù)。Su 等用KELM 預(yù)測(cè)臭氧濃度，KELM 學(xué)習(xí)速度快，方法優(yōu)于逐步回歸方法和反向傳播神經(jīng)網(wǎng)絡(luò)（BPNN）［21］。Ho 等提出一種基于PSO 算法的直覺(jué)模糊季節(jié)性回歸（IFSR）方法，提高了一氧化碳日濃度的預(yù)測(cè)性能［22］。KELM 具有結(jié)構(gòu)簡(jiǎn)單、計(jì)算速度快、精度高及實(shí)時(shí)性好等優(yōu)勢(shì)，PSO 算法能快速找到最優(yōu)的模型參數(shù)。本文將PSO 算法與KELM 相結(jié)合，克服了BP 神經(jīng)網(wǎng)絡(luò)、ELM 等算法響應(yīng)速度慢、結(jié)構(gòu)復(fù)雜、模型相關(guān)系數(shù)低等缺陷。

本文采用基于比爾朗伯定律的TDLAS 技術(shù)，設(shè)計(jì)了基于LabVIEW 的近紅外CO2濃度反演系統(tǒng)，將EMD 作為預(yù)處理算法與PSO-KELM濃度反演算法嵌入LabVIEW 中，在此基礎(chǔ)上開(kāi)展2%～20%的CO2濃度檢測(cè)實(shí)驗(yàn)。數(shù)據(jù)采集卡采集到的二次諧波信號(hào)（2f信號(hào)）經(jīng)EMD 去噪后信噪比明顯提高，最后用PSO-KELM 算法進(jìn)行濃度反演。該模型預(yù)測(cè)準(zhǔn)確度高、響應(yīng)速度快，可精準(zhǔn)測(cè)量清焦過(guò)程中的CO2含量。

2 算法原理及應(yīng)用

2.1 核極限學(xué)習(xí)機(jī)［21，23-24］

設(shè)樣本(Xi，ti)的個(gè)數(shù)為N，Xi為采集的第i條CO2的二次諧波信號(hào)，ti為對(duì)應(yīng)的CO2濃度，則包含L個(gè)隱含層神經(jīng)元（本次實(shí)驗(yàn)L=500）、激勵(lì)函數(shù)為g(x)的ELM 模型可以表示為：

其中：Wi為輸入權(quán)值，bi是隱含層神經(jīng)元的閾值，Wi和bi隨機(jī)產(chǎn)生，βi為輸出權(quán)重。ELM 網(wǎng)絡(luò)的目標(biāo)是使得模型CO2的預(yù)測(cè)濃度oj盡可能地接近CO2的設(shè)定濃度tj，即存在βi，Wi和bi使得：

用矩陣表示為：

其中：H是隱含層神經(jīng)元的輸出，β為輸出權(quán)重，T為目標(biāo)輸出濃度。

根據(jù)廣義逆的理論，其解為：

其中：H+是矩陣H的Moore-Penrose 廣義逆，且求得的的范數(shù)是最小并且唯一的。

ELM 的優(yōu)化目標(biāo)是最小化訓(xùn)練誤差和輸出權(quán)重的范數(shù)，其優(yōu)化函數(shù)為：

可持續(xù)發(fā)展指標(biāo)體系是一個(gè)由多方面指標(biāo)組成的復(fù)雜體系，指標(biāo)的選取應(yīng)該能夠適當(dāng)反映各地區(qū)在經(jīng)濟(jì)、社會(huì)、人口、資源、環(huán)境生態(tài)和科技教育等方面的基本狀況，并能滿足不同時(shí)期和不同地區(qū)對(duì)比分析的需要。為了保證指標(biāo)體系的合理性，在制定過(guò)程中遵循科學(xué)性、綜合性、代表性、可操作性、時(shí)間序列數(shù)據(jù)可比性等原則。

其中：ξi為訓(xùn)練誤差，C為誤差懲罰因子，用來(lái)權(quán)衡模型的訓(xùn)練誤差和輸出權(quán)重。

由KKT 條件，最終得到ELM 模型輸出權(quán)值：

在KELM 模型中，核矩陣Ω代替了ELM 中的隨機(jī)矩陣HHT，核矩陣可以表示為：

其中K(xi，xj)為Ωi，j的第i行、第j列光譜點(diǎn)。本次選取徑向基（Radial Basic Function，RBF）核函數(shù)，RBF 核函數(shù)為：

其中σ是核函數(shù)寬度。

則KELM 預(yù)測(cè)的輸出濃度表達(dá)式為：

2.2 基于粒子群尋優(yōu)的核極限學(xué)習(xí)機(jī)［25-26］

誤差懲罰因子C和核參數(shù)σ2是KELM 優(yōu)化的兩個(gè)關(guān)鍵參數(shù)，本文采用PSO 算法對(duì)C和σ2進(jìn)行優(yōu)化。PSO 算法是一種屬于群體智能的優(yōu)化算法，本次實(shí)驗(yàn)粒子維數(shù)D為2，粒子數(shù)為20，進(jìn)化次數(shù)設(shè)為30，xi為第i條光譜，根據(jù)目標(biāo)濃度計(jì)算xi的最小適應(yīng)度值得到：

其中：f(xi)為預(yù)測(cè)濃度，yi為目標(biāo)濃度。

樣本i的速度（即尋找最優(yōu)C和σ2的速度）為：

樣本第i個(gè)體經(jīng)歷過(guò)的最好的位置（最優(yōu)C和σ2的解）為：

群體經(jīng)歷過(guò)的最好位置（最優(yōu)C和σ2的解）為：

通常，在第d(1≤d≤D)維的位置變化限定在[xmin，d，xmax，d]內(nèi)，本次實(shí)驗(yàn)限定在［10，40 000］；速度變化范圍限定在[vmin，d，vmax，d]，本次實(shí)驗(yàn)限定在［10，1013］，即在迭代中若vid，xid超出了邊界，則該維的速度或位置被限制為該維的最大速度或邊界位置。第k次迭代后，第i個(gè)樣本的第d維帶有慣性權(quán)重因子和約束因子的速度和位置調(diào)整為：

其中：c1，c2為加速度常數(shù)，調(diào)節(jié)學(xué)習(xí)最大步長(zhǎng)，本次實(shí)驗(yàn)取默認(rèn)值2；r1，r2是兩個(gè)隨機(jī)函數(shù)，取值為［0，1］，以增加搜索隨機(jī)性；xk+1id，vk+1id為D維上第i個(gè)樣本進(jìn)行k+1 次迭代的速度和位置；ω為慣性權(quán)重因子，本次實(shí)驗(yàn)取默認(rèn)值0.9，調(diào)節(jié)對(duì)解空間的搜索范圍。

3 實(shí) 驗(yàn)

3.1 譜線選擇

燒焦過(guò)程中除了有CO2之外，還伴隨著CH4，H2O，CO 的干擾。據(jù)HIRTAN 數(shù)據(jù)庫(kù)，在濃度為0.000 1、溫度為300 K、壓力為1.01×105Pa、光程為20 m 時(shí)，CO2和CH4，H2O，CO 在1 579.4～1 580.4 nm 內(nèi)的吸收譜線如圖1 所示（彩圖見(jiàn)期刊電子版）。與CO2，CO 相比，CH4和H2O 在此波段的吸收很小，可以忽略不計(jì)，CO2在此波段內(nèi)有兩個(gè)很強(qiáng)的吸收峰，為了避免CO 氣體對(duì)CO2在1 579.75 nm 附近的強(qiáng)的吸收峰干擾，本次實(shí)驗(yàn)選取CO2在1 579.9～1 580.3 nm內(nèi)的吸收譜線。使用武漢六九傳感器公司提供的中心波長(zhǎng)為1 580 nm 的DFB 激光器（SN：19062063），用商用激光器驅(qū)動(dòng)模塊（上海昕虹光電，PN：M-DFB-GP03）和波長(zhǎng)計(jì)（HP 86120B）來(lái)檢測(cè)激光器的溫度調(diào)諧系數(shù)和電流調(diào)諧系數(shù)。首先，設(shè)定激光器的驅(qū)動(dòng)電流為50 mA，將溫度從19.5 ℃調(diào)節(jié)到34 ℃，分別記錄激光器的發(fā)光波長(zhǎng)，計(jì)算出溫度調(diào)諧系數(shù)為0.098 5 nm/℃；設(shè)定溫度為25 ℃，將驅(qū)動(dòng)電流從10 mA 調(diào)到60 mA，分別記錄激光器的發(fā)光波長(zhǎng)，計(jì)算出電流調(diào)諧系數(shù)為0.019 3 nm/mA。通過(guò)調(diào)節(jié)溫度和電流，激光器發(fā)光波長(zhǎng)的范圍可以覆蓋CO2吸收峰，并保持穩(wěn)定。

圖1 HITRAN 數(shù)據(jù)庫(kù)中CO2，H2O，CH4，CO 在1 580 nm 附近濃度為0.000 1、溫度為300 K、壓力為1.01×105、光程為20 m 時(shí)的吸收譜Fig.1 Absorption spetra of CO2，H2O，CH4 and CO in HITRAN database in range of 1 580 nm with concentration of 0.000 1，temperature of 300 K，pressure of 1.01×105 and optical path of 20 m

3.2 實(shí)驗(yàn)裝置

如圖2 所示，實(shí)驗(yàn)裝置主要包括自動(dòng)化配氣站、TDLAS 二氧化碳?xì)怏w測(cè)試系統(tǒng)、LabVIEW數(shù)據(jù)采集與濃度反演部分。自動(dòng)化配氣站對(duì)100%的純CO2和N2進(jìn)行配比，并采用松下公司生產(chǎn)的DP-101 系列氣體壓力傳感器和七星華創(chuàng)公司生產(chǎn)的精度優(yōu)于1%的數(shù)字式質(zhì)量流量計(jì)對(duì)配比后的氣體進(jìn)行控制，氣體送入長(zhǎng)光程氣室（Herriott 氣室，由武漢旭海光電公司提供，光程20 m）后經(jīng)洗氣瓶排出。 通過(guò)溫控模塊WTC3293 控制激光器溫度，將激光器的中心波長(zhǎng)調(diào)諧到1 580 nm 的波長(zhǎng)掃描范圍附近，溫度調(diào)諧是粗調(diào)方法，一旦確定后無(wú)需更改。電流控制模塊產(chǎn)生驅(qū)動(dòng)電流（5 Hz 的鋸齒信號(hào)疊加7.8 kHz 正弦信號(hào)）驅(qū)動(dòng)激光器，調(diào)整起始掃描電流和掃描電流范圍可以精細(xì)調(diào)整激光器的掃描波長(zhǎng)，使其發(fā)光波長(zhǎng)的范圍覆蓋目標(biāo)氣體的吸收峰。激光器輸出的近紅外光經(jīng)過(guò)聚焦準(zhǔn)直進(jìn)入氣室，在氣室中多次反射并被目標(biāo)氣體吸收，最終到達(dá)探測(cè)器。探測(cè)器選用InGaAs 探測(cè)器，響應(yīng)波長(zhǎng)為800～1 700 nm，將透射光信號(hào)轉(zhuǎn)換為電信號(hào)，探測(cè)器信號(hào)經(jīng)過(guò)鎖相放大器后被數(shù)據(jù)采集卡（美國(guó)國(guó)家儀器公司，SCB-68A）采集到自主開(kāi)發(fā)的LabVIEW 上位機(jī)中，通過(guò)數(shù)字信號(hào)預(yù)處理算法EMD 和濃度反演算法PSO-KELM 對(duì)采集到的數(shù)據(jù)進(jìn)行處理。

圖2 近紅外CO2濃度反演系統(tǒng)實(shí)驗(yàn)裝置Fig.2 Experimental setup of near infrared CO2 concentration inversion system

3.3 基于LabVIEW 的EMD 數(shù)據(jù)預(yù)處理算法設(shè)計(jì)

利用LabVIEW 進(jìn)行數(shù)據(jù)采集時(shí)，首先在NIDAQmx 生成數(shù)據(jù)采集任務(wù)，選擇采集通道及觸發(fā)方式。本次實(shí)驗(yàn)需采集TDLAS 測(cè)試系統(tǒng)傳出的光功率（DC）信號(hào)、二次諧波信號(hào)（2f信號(hào)）以及氣體壓力（P）信號(hào)（均為電壓信號(hào)），采樣點(diǎn)為1 k，采樣頻率為5 k，即光譜信號(hào)每個(gè)周期為200 ms；在STM32 微處理芯片的程序中，通過(guò)外部中斷處理函數(shù)，每隔200 μs 觸發(fā)一次定時(shí)器中斷，在定時(shí)器中斷函數(shù)中調(diào)用數(shù)模轉(zhuǎn)換器（Digital to Analog Convertor，DAC）模擬信號(hào)輸出函數(shù)，每200 μs 輸出一個(gè)模擬量，每1 000 個(gè)點(diǎn)為一個(gè)周期即完成了一個(gè)周期慢掃信號(hào)的輸出。當(dāng)計(jì)數(shù)輸出了1 000 個(gè)模擬量時(shí)，向指定管腳發(fā)出一個(gè)高電平信號(hào)，該信號(hào)即LabVIEW 程序采集的外部觸發(fā)信號(hào)，所以本實(shí)驗(yàn)選擇上升沿觸發(fā)方式。

由于經(jīng)數(shù)據(jù)線采集到的數(shù)據(jù)傳到LabVIEW后有很大的噪聲干擾，所以要對(duì)數(shù)據(jù)進(jìn)行預(yù)處理。本文程序中利用MATLAB 程序生成的組件對(duì)象模型（Component Object Model，COM）以Win32 動(dòng)態(tài)鏈接庫(kù)（DLL）的形式被LabVIEW 調(diào)用，實(shí)現(xiàn)對(duì)算法的跨平臺(tái)調(diào)用和對(duì)實(shí)時(shí)數(shù)據(jù)的分析。在軟件中嵌入了經(jīng)驗(yàn)?zāi)B(tài)分解（Empirical Mode Decomposition，EMD）算法，EMD 算法是一種自適應(yīng)時(shí)頻分析方法，它將原始信號(hào)分解為有限個(gè)固有模態(tài)函數(shù)（Intrinsic Mode Function，IMF）和一個(gè)余項(xiàng)，所分解出來(lái)的各IMF 分量包含了原始信號(hào)的不同時(shí)間尺度的局部特征信號(hào)，每個(gè)IMF 必須要滿足在整個(gè)信號(hào)上極值點(diǎn)的個(gè)數(shù)和過(guò)零點(diǎn)的個(gè)數(shù)相差不大于1 且在任意點(diǎn)處上下包絡(luò)的均值為0 的條件。首先，找到原始2f信號(hào)x(t)的所有極大值點(diǎn)與所有極小值點(diǎn)，通過(guò)三次樣條函數(shù)擬合出極大值包絡(luò)線與極小值包絡(luò)線，上下包絡(luò)線的均值作為原始2f信號(hào)的均值包絡(luò)m(t)，用原始2f信號(hào)x(t)減去均值包絡(luò)m(t)就得到一個(gè)去掉低頻的新信號(hào)xnew(t)。若新信號(hào)不滿足S≤a，則重復(fù)上述過(guò)程：

其中a為篩分門(mén)限，本次實(shí)驗(yàn)取0.3；若新信號(hào)滿足S＞a的條件，則用原始信號(hào)減去第一個(gè)IMF分量后再重復(fù)上述過(guò)程，直到余項(xiàng)r(t)是單調(diào)函數(shù)或常量時(shí)EMD 分解過(guò)程停止。通過(guò)EMD 算法分解系統(tǒng)噪聲和光學(xué)基線，將這幾項(xiàng)噪聲去除后重組有用的光譜信號(hào)，可提高光譜信噪比［27-28］。

LabVIEW 程序采用多次平均處理可事先去除一部分噪聲，利用while 循環(huán)框上的移位寄存器實(shí)現(xiàn)光功率信號(hào)（DC 信號(hào)）和二次諧波信號(hào)（2f信號(hào)）的累加，累加次數(shù)即為平均次數(shù)，本次實(shí)驗(yàn)的平均次數(shù)為10。因多次平均去噪效果有限，所以在2f信號(hào)后添加EMD 處理，并分別將原始信號(hào)raw2f與經(jīng)EMD 處理后的信號(hào)存儲(chǔ)到文件中。圖3（a）為單條raw2f信號(hào)，噪聲很強(qiáng)，信噪比為6.75 dB；圖3（b）為raw2f信號(hào)經(jīng)過(guò)EMD 分解后得到3 個(gè)噪聲分量和1 個(gè)res 分量，需要剔除掉；圖3（c）為raw2f信號(hào)去掉圖3（b）中分解出的信號(hào)后得到的重組信號(hào)，即2f信號(hào)，可以看出噪聲被有效去除，信號(hào)平滑了很多，信噪比為12.59 dB，比原始信號(hào)的信噪比提高了1.87 倍。

圖3 MATLAB 仿真結(jié)果Fig.3 MATLAB simulation result

4 實(shí)驗(yàn)結(jié)果與分析

4.1 步進(jìn)實(shí)驗(yàn)

為了進(jìn)一步驗(yàn)證EMD 算法的處理效果，用純氮?dú)馀c100%濃度的CO2進(jìn)行混合配比，分別開(kāi)展了0.03%～0.08%（步進(jìn)為0.01%）的6 組低濃度實(shí)驗(yàn)與2%～15%（步進(jìn)為1%）、18%，20% 的16 組高濃度實(shí)驗(yàn)，每條光譜是1 000 個(gè)點(diǎn)，去掉前后各250 個(gè)不穩(wěn)定的點(diǎn)，取中間500 個(gè)數(shù)據(jù)點(diǎn)進(jìn)行分析。圖4（a）為經(jīng)EMD 處理后的16組高濃度CO2光譜數(shù)據(jù)。每個(gè)濃度采30 多條光譜，取每個(gè)濃度譜峰的峰峰值進(jìn)行平均，并與CO2濃度進(jìn)行線性擬合，得到的結(jié)果如圖4（b）所示，線性相關(guān)系數(shù)R2=0.997 6。圖4（c）為經(jīng)EMD 處理后的6 組低濃度CO2光譜數(shù)據(jù)及其峰峰值與濃度的線性相關(guān)性曲線，相關(guān)系數(shù)R2=0.998 2。

圖4 經(jīng)EMD 處理后的CO2光譜數(shù)據(jù)及其峰峰值與濃度的線性關(guān)系Fig.4 CO2 spectral data after EMD processing and linear relationship between peak-to-peak value and concentration

4.2 PSO-KELM 與BP 神經(jīng)網(wǎng)絡(luò)和ELM 的算法對(duì)比

因本次主要是應(yīng)用于燒焦過(guò)程的CO2氣體濃度實(shí)時(shí)監(jiān)測(cè)，所以僅對(duì)高濃度（2%～20%）的CO2進(jìn)行算法分析。16 組經(jīng)過(guò)EMD 處理的高濃度CO2光譜共527 條，采用5 折校驗(yàn)法，422 條光譜作為訓(xùn)練集，105 條光譜作為測(cè)試集，通過(guò)PSO算法對(duì)RBF 核函數(shù)的核參數(shù)即核半徑σ2以及懲罰系數(shù)C進(jìn)行尋優(yōu)，發(fā)現(xiàn)當(dāng)C=4.397×103，σ2=4.947×1010時(shí)，對(duì)測(cè)試集105 個(gè)樣本的預(yù)測(cè)精度達(dá)到最高，此時(shí)精度為962×10-6。

將PSO-KELM（采用RBF 核函數(shù)）算法的濃度反演結(jié)果與PLS、CLS、BP 神經(jīng)網(wǎng)絡(luò)（激活函數(shù)為ReLU）及ELM（激活函數(shù)為ReLU）進(jìn)行對(duì)比，結(jié)果如表1 所示。PSO-KELM 反演結(jié)果的均方根誤差最小，線性相關(guān)系數(shù)最接近1。由此可知，PSO-KELM 的預(yù)測(cè)準(zhǔn)確度最高，均方根誤差比PLS 小了1.22 倍，比CLS 小了15.10 倍，比ELM 小了4.08 倍，比BP 神經(jīng)網(wǎng)絡(luò)小了6.53 倍。

表1 步進(jìn)測(cè)試CO2 氣體中BP 神經(jīng)網(wǎng)絡(luò)、ELM 和PSOKELM 算法的擬合精度對(duì)比Tab.1 Comparison of fitting accuracy of BP neural network，ELM and PSO-KELM algorithms for step test of CO2 gas

圖5 為BP 神經(jīng)網(wǎng)絡(luò)、ELM、PSO-KELM 3 種算法模型的擬合結(jié)果。

圖5 BP 神經(jīng)網(wǎng)絡(luò)、ELM 和PSO-KELM 算法的擬合結(jié)果對(duì)比Fig.5 Comparison of fitting results of BP neural network，ELM and PSO-KELM algorithms

4.3 穩(wěn)定性測(cè)試

為了測(cè)試檢測(cè)系統(tǒng)的穩(wěn)定性，在70 kPa 的壓力下通入10%的CO2，持續(xù)1.5 h，測(cè)得623 條光譜。如圖6（a）所示，測(cè)量濃度在9.940%～10.056%內(nèi)波動(dòng)，標(biāo)準(zhǔn)差為196.6×10-6。如圖6（b）所示，對(duì)測(cè)量結(jié)果進(jìn)行Allan 方差分析，積分時(shí)間為2 s 時(shí)，Allan 方差為90.9×10-6；隨著積分時(shí)間的增加，當(dāng)積分時(shí)間為16 s 時(shí)，Allan 方差為44.8×10-6，此時(shí)TDLAS 測(cè)試系統(tǒng)的穩(wěn)定性最好，達(dá)到了極限標(biāo)準(zhǔn)差。由此可見(jiàn)，TDLAS 測(cè)試系統(tǒng)具有高度穩(wěn)定性。

圖6 穩(wěn)定性測(cè)試結(jié)果與分析Fig.6 Stability test results and analysis

5 結(jié) 論

本文基于TDLAS 技術(shù)，以中心波長(zhǎng)為1 580 nm 的DFB 激光器（SN：19062063）作為光源，光程為20 m 的Herriott 氣室作為光與CO2作用的場(chǎng)所，開(kāi)展燒焦過(guò)程的CO2氣體檢測(cè)實(shí)驗(yàn)。用數(shù)據(jù)采集卡采集原始二次諧波信號(hào)，提出了將EMD 作為預(yù)處理算法嵌入到LabVIEW 上位機(jī)中，經(jīng)EMD 處理后的2f信號(hào)的信噪比提高了近2倍。將PSO-KELM 算法模型應(yīng)用于CO2氣體濃度反演，與CLS，PLS，BP 神經(jīng)網(wǎng)絡(luò)和ELM 算法做比較，其均方根誤差最小，線性相關(guān)系數(shù)最接近1。最后，在氣室壓力在70 kPa，CO2濃度為10%的條件下進(jìn)行1.5 h 的長(zhǎng)期穩(wěn)定性測(cè)試，測(cè)量濃度在9.940%～10.056%內(nèi)波動(dòng)，標(biāo)準(zhǔn)差為196.6×10-6。對(duì)測(cè)量結(jié)果進(jìn)行Allan 方差分析，積分時(shí)間為2 s 時(shí)，Allan 方差為90.9×10-6；隨著積分時(shí)間的增加，當(dāng)積分時(shí)間為16 s 時(shí)，Allan 方差為44.8×10-6，此時(shí)TDLAS 測(cè)試系統(tǒng)的穩(wěn)定性最好，達(dá)到了極限標(biāo)準(zhǔn)差。因此，該系統(tǒng)可以實(shí)現(xiàn)對(duì)燒焦過(guò)程中CO2氣體濃度的高靈敏度實(shí)時(shí)準(zhǔn)確監(jiān)測(cè)，提升工藝生產(chǎn)的經(jīng)濟(jì)效益。