馬維光,周曉彬,曹振松,許非,田建飛,周月婷,劉建鑫,趙剛*

(1 山西大學激光光譜研究所量子光學與光量子器件國家重點實驗室,山西 太原 030006;2 中國科學院合肥物質(zhì)科學研究院安徽光學精密機械研究所,中國科學院大氣光學重點實驗室,安徽 合肥 230031;3 山西大學極端光學協(xié)同創(chuàng)新中心,山西 太原 030006)

0 引言

根據(jù)歐盟哥白尼氣候變化服務中心報告,2020 年成為近十年最高溫年份。相較于工業(yè)革命之前,全球年平均氣溫上升了1.2°C,并且上升率逐漸增加。二氧化碳(CO2)是最主要的溫室氣體,對全球變暖的貢獻達到了60%以上。化石燃料的燃燒造成全球CO2平均濃度相較于工業(yè)革命前增加了50%。因此,對大氣中CO2氣體濃度進行實時、高靈敏的監(jiān)測有助于了解CO2的排放源、濃度變化趨勢以及與惡劣氣候的關系等。

傳統(tǒng)的氣體檢測方式包括電化學式[1]、接觸燃燒[2]、熱導式[3]等,它們具有靈敏度低、需要采樣分析、響應速度慢等缺點。而激光吸收光譜技術(LAS)基于光與氣體分子相互作用,可以實時在線、高靈敏地響應氣體濃度的變化。直接吸收光譜技術是最基本的LAS,但是由于受到相對強度噪聲的影響,其最小可探測光強衰減率通常只能達到10-3。為了提升探測靈敏度,通常采用兩種方法:第一種是使用調(diào)制技術抑制1/f噪聲的影響,基于此發(fā)展了波長調(diào)制光譜技術[4]和頻率調(diào)制光譜技術[5];第二種是通過使用多通道池[6,7]或者光學腔[8-10]增長激光與氣體介質(zhì)作用長度,從而增強吸收信號。而腔衰蕩光譜(CRDS)技術是腔增強光譜技術(CEAS)的一種,它通過測量光子在腔內(nèi)的壽命獲取腔內(nèi)的吸收損耗[10,11]。因此,CRDS 不僅借助光學腔增強吸收信號,并且由于通過測量時間信號來反演吸收信息,避免了激光強度噪聲的影響,從而可以獲得比光頻鎖定CEAS 更高的探測靈敏度。

1988 年O’Keefe 和Deacon 首次提出CRDS[12],他們使用脈沖光作為光源,觀測到了腔衰蕩信號,對可見光區(qū)域的氧分子吸收光譜進行了測量。1997 年,Romanini 等實現(xiàn)了基于連續(xù)激光源的CRDS(CW-CRDS)[13],由于連續(xù)激光源具有更好的重復性和穩(wěn)定性,CW-CRDS 可以獲得更好的探測靈敏度,在此基礎上,他們探測了乙炔在570 nm 處的吸收線,噪聲等效吸收達到了10-9cm-1。

目前,世界上有多個公司和機構都進行了基于CRDS 的氣體分析儀的開發(fā),以美國的Picarro 公司和Tigeroptics 公司的產(chǎn)品為典型代表。Picarro 公司精準控制系統(tǒng)溫度(穩(wěn)定性優(yōu)于20 mK)和腔內(nèi)壓強(控制精度優(yōu)于±50.65 Pa)的同時,還基于標準具開發(fā)了一種高精度的波長監(jiān)視器,很好地解決了無法準確獲取激光二極管輸出波長的問題。得益于以上措施,他們的CRDS 氣體分析儀體積比檢測極限達到了10-12量級[14]。而Tigeropitcs 的儀器以AOM 為光開關,獲得了1.0 nmol/mol 的最低檢測極限[15]。這兩家公司的分析儀結(jié)構緊湊,檢測靈敏度和穩(wěn)定性也很高,國內(nèi)相關產(chǎn)品與之相比還存在一定差距。

在這樣的環(huán)境背景下,發(fā)展國產(chǎn)高性能CRDS 氣體分析儀具有重要的現(xiàn)實意義。本文詳細介紹了基于CW-CRDS 的氣體分析儀集成化研制過程,包括理論分析、光路和電路設計以及處理程序;然后基于該系統(tǒng),實現(xiàn)了對不同濃度CO2氣體的探測;最后通過評估空腔衰蕩時間,評估了系統(tǒng)的探測極限。

1 理 論

在CW-CRDS 技術中,當激光頻率與光腔共振時激光會高效耦合入光腔,如突然關斷激光,腔透射光強會隨時間衰減。當腔內(nèi)存在氣體吸收時,定義激光往返一周損耗因子為

式中:R為腔鏡光強反射率;L為腔長;α 為腔內(nèi)氣體單通吸收系數(shù);αt為光強單通衰減系數(shù),可表示為

由此可得,腔內(nèi)光強隨時間變化的表達式為

式中:tr=2L/c為光在腔內(nèi)往返一周的時間,I0為入射光強。對于腔內(nèi)無吸收的情況,此時α=0,腔內(nèi)光強為

對于高精細度的腔,腔鏡具有高的反射率,可做近似

因此腔透射光強可表示為

式中:τ0為空腔衰蕩時間,由腔長L、光速c、腔鏡反射率R決定,即

當有腔內(nèi)存在吸收時,吸收系數(shù)α 是光學頻率ν 的函數(shù),衰蕩時間τ(ν)滿足

由此可知,腔內(nèi)氣體樣品的吸收越強,測得的衰蕩時間就越短。綜合(7)、(8)式可得氣體吸收系數(shù)

通過測量不同激光波長時的吸收系數(shù)即可獲得完整的吸收線型。由Beer-Lambert 定律可知,α(ν)可以表示為

式中:S為吸收線強度,單位是cm-2·Pa-1;φ(ν)為面積歸一化的線型函數(shù);P為腔內(nèi)氣壓,單位為Pa;x為待測氣體的濃度。

2 研究方案與裝置

2.1 光路結(jié)構

圖1 為CW-CRDS 氣體分析裝置結(jié)構圖,系統(tǒng)采用DFB 激光器(NTT,NLK1L5EAAA) 作為激光源,中心波長為1578 nm,可以覆蓋CO2位于6359.97 cm-1處的吸收線,吸收線強度為1.78×10-23cm-1/(mol·cm-2)[16]。激光器的輸出頻率隨著溫度和電流的變化而改變,通過自制三角波掃描電路產(chǎn)生變化的電壓信號后接入激光器電流源掃描激光的輸入電流,可以實現(xiàn)對激光器輸出波長的連續(xù)調(diào)諧。輸出激光先通過一個光纖隔離器(Thorlabs,IO-G-1550-APC)以防止反射光對激光器產(chǎn)生光學反饋,后通過一個3 dB 的法蘭式固定增益衰減器將光功率衰減到5 mW 以下,滿足半導體光放大器(SOA)(北京科揚光電技術有限公司,KY-PLM-15-M-FA)的輸入要求。衰減后的光從SOA 的入射端輸入,放大后的激光從出射端輸出,通過一個光纖準直器(OZ Optics,HPUCO-T,3A-1530/1585-P-3.9AS)將光纖中的光場轉(zhuǎn)化為空間光場。在系統(tǒng)中,SOA 主要有兩個作用:一是作為光學開關,觸發(fā)衰蕩事件;二是對激光功率進行放大,放大系數(shù)為13.2 dB,用來提升信噪比。

圖1 CW-CRDS 氣體分析裝置結(jié)構圖Fig.1 Schematic diagram of CW-CRDS gas analyzer

實驗中使用的Fabry-Perot(F-P)光學腔由兩個高反鏡組成,反射率皆為99.95%,對應腔的精細度約為6500。兩個高反鏡(長春新產(chǎn)業(yè)光電技術有限公司)分別粘在低膨脹腔體兩端。腔長為40 cm,對應自由光譜區(qū)(FSR)約為380 MHz。使用一個透鏡組完成激光到高精細度F-P 腔的模式匹配。通過掃描激光器電流來調(diào)諧激光頻率,當激光頻率恰好滿足腔的駐波條件時,激光與腔達到共振,在腔內(nèi)建立起強共振光場。腔透射光射入一個雪崩探測器(Thorlabs,APD110C/M),用于測量腔衰蕩信號。

2.2 電路結(jié)構

裝置的電路部分主要由激光器驅(qū)動和信號處理兩部分組成。激光器驅(qū)動電路包括三角波掃描和激光器溫控電流源電路。激光器的溫控電流源基于Wavelength 公司的FL500 (電流工作范圍設置為0～250 mA,電流噪聲為3 μA)和WTC3243,電路采用了芯片的標準引用。如圖2(a)的三角波掃描電路輸出一個三角波模擬信號,其中紅色虛線內(nèi)為三角波發(fā)生模塊,用于產(chǎn)生三角波信號,FREQ 電位器用于調(diào)節(jié)三角波頻率,AMP 電位器用于調(diào)節(jié)三角波幅度；橙色虛線內(nèi)為觸發(fā)信號發(fā)生模塊,用于產(chǎn)生與三角波信號同步的方波信號;綠色虛線內(nèi)為偏置調(diào)節(jié)模塊,OFFSET 電位器用于調(diào)節(jié)三角波信號的偏置;藍色虛線內(nèi)是一個加法器,將偏置電壓加到三角波信號上。產(chǎn)生的三角波信號送入激光器電流驅(qū)動模塊,用于調(diào)諧激光器頻率,本系統(tǒng)激光頻率掃描范圍約為10 GHz,掃描頻率約為3 Hz。

信號處理部分主要包括腔模信號采集,觸發(fā)衰蕩事件以及擬合腔衰蕩信號。探測器采集到腔模信號后分為兩路輸出:一路輸入數(shù)據(jù)采集卡(National Instruments,PCI-6115)后利用計算機進行數(shù)據(jù)處理,采樣率設置為10 MHz;另一路送入自制的閾值電路,原理如圖2(b)所示,用于控制SOA 開斷來觸發(fā)衰蕩事件。為了防止兩路信號相互干擾,使用一個跟隨器電路將信號分離。閾值電路的核心芯片是lm555,當輸入信號沒有達到設定閾值時,閾值電路輸出5 V 高電平,SOA 將光功率放大后輸出;當輸入信號達到閾值后,閾值電路立刻輸出一個2 μs 的低脈沖信號,SOA 關斷激光,產(chǎn)生衰蕩信號,采集完衰蕩信號后,閾值電路輸出自動回復到高電平,從而開啟SOA。

圖2 (a)三角波掃描電路及(b)閾值電路原理圖Fig.2 Schematic diagram of(a)triangular wave scanning circuit and(b)threshold circuit

2.3 機械結(jié)構

如圖3(a)所示,整個系統(tǒng)被集成在一個700 mm×300 mm×185 mm 的定制機箱中,由一個線性電源模塊供電,可以提供一路±15 V、一路±10 V 和三路+5 V 低噪聲電壓,分別驅(qū)動探測器、三角波電路、激光電流源溫控電路[見圖3(c)]、閾值電路和SOA。腔體[見圖3(b)](大連齊維理化器械有限公司)材料為超低熱膨脹系數(shù)的殷鋼,側(cè)面設有出氣和進氣口,通過管路與裝置外殼連接用于對環(huán)境空氣進行實時探測。系統(tǒng)中的鏡桿、支架、電路支撐盤均由硬鋁加工而成并固定在底板上。系統(tǒng)的信號輸出分別是探測器信號輸出(來自探測器,接入數(shù)據(jù)采集卡進行處理)和觸發(fā)信號TTL 輸出(來自三角波掃描電路的方波輸出,用于觸發(fā)數(shù)據(jù)采集卡采集信號)。

圖3 (a)CW-CRDS 氣體分析裝置內(nèi)部結(jié)構圖;(b)F-P 腔實物圖;(c)激光驅(qū)動電路實物圖Fig.3 (a)Mechanical structure diagram of CW-CRDS gas analyzer;(b)Picture of F-P cavity;(c)Picture of laser driver circuit

3 測量結(jié)果

采集卡采集到的信號送入計算機由Labview 程序處理。為了提升對光譜的采樣速率,在完成10 GHz激光頻率掃描后對透射腔模進行逐個處理。通過腔模尋峰,程序可以獲取每一個腔衰蕩信號觸發(fā)的位置,而后使用e 指數(shù)函數(shù)對衰蕩信號進行擬合,結(jié)果如圖4 所示,圖中圓點為衰蕩信號,曲線為擬合曲線。為了減小噪聲影響,本系統(tǒng)選取了衰蕩時間8 倍時長作為采樣總長度,并且考慮到閾值電路和SOA 的響應時間,去除了衰蕩起始處的前8 個數(shù)據(jù)點。

圖4 (a)衰蕩信號及擬合曲線;(b)擬合殘差Fig.4 (a)Ring-down signal and fitting curve;(b)The fitting residual

將一個掃描范圍內(nèi)所有的衰蕩事件進行擬合后,確定出衰蕩時間τ 與激光頻率的關系,結(jié)果如圖5(a)所示,其中橫坐標為相對激光頻率,利用腔模間隔與自由光譜區(qū)的關系來定標,可以明顯觀察到由于氣體的吸收,衰蕩時間變小,并且越靠近吸收線中心,衰蕩時間越小。利用(4)式可以將衰蕩時間直接轉(zhuǎn)換為吸收系數(shù)α,因此該技術是一種免定標光譜測量技術,得到的結(jié)果如圖5(b)所示。由于腔內(nèi)氣壓約為9.3219×104Pa(太原本地氣壓),線型函數(shù)采用面積歸一化洛倫茲函數(shù)來近似表示[17],基于(10)式對圖5(b)中的吸收線型進行擬合,如圖中紅線所示,以此可確定出CO2氣體的濃度。圖5(c)顯示了兩者的殘差,可見峰值殘差小于0.3%。

圖5 (a)衰蕩時間隨相對頻率的變化曲線;(b)吸收系數(shù)變化曲線及擬合曲線;(c)擬合殘差Fig.5 (a)Curve of ring-down time versus relative frequency;(b)Absorption coefficient extracted from ringdown time and its fitting curve;(c)The fitting residual

為了驗證系統(tǒng)濃度測量的準確性,對體積占比分別為2×10-4、3×10-4、4×10-4、5×10-4、6×10-4的CO2標準氣進行了實時測量,濃度測量結(jié)果如圖6(a)所示,其中每個濃度測量500 s,基于該測量可以獲得每種濃度的測量平均值和標準偏差。圖6(b)為測量濃度的平均值與標氣濃度值的對應關系,其中黑色點為測量濃度均值,紅色線為線性擬合結(jié)果,誤差條為相應濃度的標準偏差。通過線性擬合,可見其斜率為0.992,R2為0.99963。

圖6 (a)不同濃度CO2 標準氣測量結(jié)果;(b)測量濃度的平均值與標氣濃度值的對應關系Fig.6 (a)Measurement results of different concentrations of CO2 standard gas;(b)Correspondence between the average value of the measured concentration and the calibration gas concentration value

為了確定本套系統(tǒng)的探測極限,對空腔衰蕩時間進行了10 min 的連續(xù)測量,結(jié)果如圖7 所示。由圖可見空腔衰蕩時間的平均值約為2.818 μs,標準偏差Δτ 為0.00917 μs,將該數(shù)據(jù)代入

圖7 空腔衰蕩時間Fig.7 Cavity ring-down time without absorption

后可得系統(tǒng)的探測極限為3.85×10-8cm-1[18],對應CO2氣體的最小可探測體積比為4×10-6。

4 結(jié)論

基于CRDS 技術研發(fā)了一套高集成化的CO2氣體分析裝置,實現(xiàn)了CO2濃度的長期、實時監(jiān)測。詳細介紹了所研發(fā)裝置的主要器件及信號處理過程。本裝置作為原型機很好地完成了室內(nèi)CO2氣體實時監(jiān)測的任務,獲得了4×10-8cm-1的探測極限,對應可探測CO2體積比為4×10-6。其探測極限主要受到電子噪聲的限制,有待進一步改進。