薛正躍 李竣 劉笑海 王晶晶 高曉明 談圖?

1) (中國科學(xué)技術(shù)大學(xué)環(huán)境科學(xué)與光電技術(shù)學(xué)院,合肥 230031)

2) (中國科學(xué)院合肥物質(zhì)科學(xué)研究院安徽光學(xué)精密機(jī)械研究所,合肥 230031)

激光外差光譜探測由于其光譜分辨率高、體積小、重量輕等優(yōu)點(diǎn)近年來得到了快速的發(fā)展,可用于大氣溫室氣體垂直廓線測量和碳衛(wèi)星地面定標(biāo)等.本文報(bào)道了利用3.939 μm 帶間級聯(lián)激光器作為本振光源的測量大氣N2O 的激光外差系統(tǒng),自制高精度太陽跟蹤儀收集太陽光作為激光外差的信號光源,其跟蹤精度達(dá)到7 arcsec,激光外差系統(tǒng)的光譜分辨率達(dá)到0.004 cm—1,測量了合肥地區(qū)(31.902°N,117.167°E)大氣N2O 吸收光譜,得到2838.336 和2539.344 cm—1 兩個(gè)強(qiáng)吸收峰,并對吸收信號進(jìn)行波長標(biāo)定,得到了N2O 分子的整層大氣透過率譜,信噪比為93.將高分辨率光譜數(shù)據(jù)進(jìn)行歸一化處理和頻率校正,利用參考正向模型和最優(yōu)估計(jì)算法得到N2O 大氣整層濃度廓線,標(biāo)準(zhǔn)偏差體積分?jǐn)?shù)為0.000031 × 10—6—0.0026 × 10—6,對應(yīng)相對誤差范圍為0.009%—0.83%.研究結(jié)果表明,所搭建的激光外差系統(tǒng)能夠?qū)崿F(xiàn)對大氣中N2O 的吸收光譜測量以及對N2O 的廓線反演,為長期觀測大氣N2O 濃度提供保證.

1 引言

氧化亞氮(N2O)是大氣中僅次于CO2和CH4的一種重要的溫室氣體.全球大氣中氧化亞氮的體積分?jǐn)?shù)已從工業(yè)化前約270 × 10—9,增加到現(xiàn)在的340 × 10—9[1].氧化亞氮總排放量中超過三分之一是來源于人為排放[2],主要來自于農(nóng)業(yè),但是目前對N2O 的時(shí)空分布尚不明確.因此,持續(xù)準(zhǔn)確地檢測大氣中N2O 的時(shí)空分布對氣候變化的影響具有重要意義.

基于光譜技術(shù)的大氣探測方法主要有傅里葉變換紅外光譜技術(shù)[3,4]、激光外差光譜技術(shù)[5].激光外差光譜技術(shù)因其光譜分辨率高、空間分辨率高、光譜采樣時(shí)間短、體積小、靈敏度高等優(yōu)點(diǎn)被廣泛應(yīng)用,其利用窄線寬的本振激光與太陽光或其他光源得到攜帶整層大氣分子吸收信息的吸收光譜.Palmer 等[6]、Weidmann 等[7]和Tsai 等[8]利用外腔量子級聯(lián)激光器建立的中紅外激光外差測量系統(tǒng),實(shí)現(xiàn)了對大氣中CH4,CO2,H2O,O3和N2O 等多種溫室氣體的測量.美國NASA 的Wilson 等[9,10]利用小型的分布反饋式激光器建立了全光纖近紅外激光外差測量系統(tǒng),獲取了大氣中CH4與CO2的吸收光譜,同時(shí)對裝置進(jìn)行了集成化設(shè)計(jì).Rodin等[11]利用平衡探測器建立激光外差測量系統(tǒng)獲取了大氣中CH4與CO2的吸收光譜.國內(nèi),中國科學(xué)院安徽光學(xué)精密機(jī)械研究所也開展了以外腔量子級聯(lián)激光器與帶間級聯(lián)激光器為本振光源的自由空間結(jié)構(gòu)激光外差光譜技術(shù)的研究[12,13].

本文基于激光外差探測技術(shù)結(jié)合反演算法,設(shè)計(jì)了一套用于測量大氣N2O 濃度的激光外差系統(tǒng),其中使用3.939 μm 的帶間級聯(lián)激光器(interband cascade laser,ICL)作為光源,采用N2O 位于2538.336 cm—1的吸收譜線對大氣中的N2O 濃度進(jìn)行了探測.實(shí)驗(yàn)中采用自制高精度太陽跟蹤儀采集太陽光作為信號光源,利用ICL 作為本振光源,兩束空間光進(jìn)行光學(xué)混頻后,差頻信號經(jīng)過光電探測器、射頻電路及鎖相放大器,并最終由美國國家儀器(National Instruments,NI)采集卡采集得到高分辨率光譜數(shù)據(jù).將所采集的光譜數(shù)據(jù)進(jìn)行歸一化處理、頻率校正之后,結(jié)合參考正向模型和最優(yōu)估計(jì)算法反演得到整層大氣N2O 的濃度廓線.

2 實(shí)驗(yàn)原理

激光外差原理在文獻(xiàn)[14]中已有詳細(xì)描述,在此只做簡單描述.微弱的輸入信號與較強(qiáng)的本振(LO)信號混合在非線性電路元件中,差頻或中頻信號(IF)由濾波器選擇后,在窄帶固定頻率放大器中放大并檢測,中頻信號可以選擇落在可用放大裝置的頻譜范圍內(nèi).圖1 給出了顯示激光外差的原理示意圖.兩束光束疊加在一個(gè)混頻器上,然后經(jīng)過放大后檢測中頻信號.

圖1 激光外差原理圖Fig.1.Schematic diagram of laser heterodyne.

大多數(shù)光學(xué)探測器是電場平方律器件,響應(yīng)與光照強(qiáng)度成正比,假設(shè)兩個(gè)振幅為As和ALO、角頻率為ws和wLO且線性極化的完全對準(zhǔn)平面波入射到探測器上.探測器所產(chǎn)生的總光電流理論值為

經(jīng)過三角函數(shù)和差化積公式簡化為

式中,cos2wst,cos2wLOt,cos(wst+wLOt) 頻率遠(yuǎn)遠(yuǎn)大于光電探測器電路帶寬,所以積分輸出為0.如果ws和wLO足夠接近,ws-wLO在光電探測器的低通截止頻率之下.那么(2)式可以進(jìn)一步簡化為

其中,光電探測器輸出信號包含兩個(gè)部分:直流信號和交流信號.

直流信號表達(dá)式:

交流信號表達(dá)式:

直流信號反映了光電探測器對本振光和信號的直流響應(yīng),交流信號即是激光外差信號,如果本振光和信號光都是理想的單色光,那么交流部分的信號是一個(gè)頻率為兩束光光頻之差ws-wLO的正弦波.實(shí)際上,外差的功率表達(dá)式為

其中RL為光電探測器的負(fù)載電阻.

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

3.1 太陽跟蹤儀

激光外差技術(shù)可以在很小的儀器體積下獲得很高的光譜分辨率,光譜分辨率約為幾十兆赫茲.為了提高外差探測精度,過去主要是在光譜儀本身進(jìn)行改進(jìn),例如儀器線性函數(shù)(ILS)表征[15],在干涉圖上應(yīng)用直流信號(DC)校正、改進(jìn)了采樣精度和探測器非線性的表征.為了達(dá)到足夠的準(zhǔn)確度,必須考慮激光外差探測技術(shù)對信號光和本振光的空間模式有嚴(yán)格的要求,而太陽跟蹤儀為激光外差探測提供信號光,必須保證其提供的信號光具有足夠高的精度及穩(wěn)定度,才能滿足測量的要求.

太陽跟蹤儀跟蹤方式有基于太陽運(yùn)行軌跡方式和光電跟蹤兩種方式,基于太陽運(yùn)行軌跡方式在任何氣象條件下都能正常跟蹤,但是沒有形成環(huán)路反饋,精度有所欠缺,光電跟蹤方式利用太陽光斑進(jìn)行實(shí)時(shí)跟蹤,缺點(diǎn)是在陰天氣候無法工作.因此我們將兩種跟蹤方式結(jié)合起來提高跟蹤精度,首先利用已知天文數(shù)據(jù)算出太陽高度角h0和方位角A.

太陽高度角計(jì)算公式為

太陽方位角計(jì)算公式為

其中θ為赤緯角,β為緯度,γ為時(shí)角.上位機(jī)Labview 算法程序在計(jì)算獲得太陽高度角和方位角之后,得到與當(dāng)前兩角度差值,轉(zhuǎn)換成步進(jìn)電機(jī)所需要的轉(zhuǎn)動步數(shù),再由Labview 上位機(jī)程序控制電機(jī)轉(zhuǎn)動,實(shí)現(xiàn)太陽運(yùn)動軌跡跟蹤.上述跟蹤完成后,在圖像傳感器判斷是否接收到太陽光斑信息,如果沒有接收到太陽光斑信息,系統(tǒng)繼續(xù)運(yùn)行太陽運(yùn)動軌跡跟蹤程序,一直到圖像傳感器接收到太陽光斑信息,Labview 程序通過對太陽光斑進(jìn)行二值化和中值濾波處理,獲得太陽光斑質(zhì)心坐標(biāo),計(jì)算太陽光斑質(zhì)心坐標(biāo)與視場中心坐標(biāo)偏差,并轉(zhuǎn)換成步進(jìn)電機(jī)轉(zhuǎn)動步數(shù),用PID (proportionintegral-differential)控制算法實(shí)時(shí)不斷的調(diào)整電機(jī)轉(zhuǎn)動步數(shù).

太陽跟蹤儀部分的結(jié)構(gòu)裝置圖如圖2 所示,太陽光經(jīng)過兩面安裝在電機(jī)上的高反鏡反射,白片反射的太陽光穿過巴德膜之后打到CCD 光電探測器,穿過氟化鈣白片的光與本振激光拍頻產(chǎn)生外差信號.

如圖3 所示,太陽跟蹤儀在精跟蹤模式下,X軸的偏移量在[—0.3,0.3]像素之間,Y軸偏移量在[—0.5,0.5]像素之間,1 Pixel 約等于14 arcsec,所以整體跟蹤精度約為7 arcsec.

圖2(a) 太陽跟蹤儀內(nèi)部裝置;(b)精跟蹤模式下太陽光斑質(zhì)心偏移量Fig.2.(a) Internal device of suntracker;(b) sunlight spot centroid offset in precision tracking mode.

圖3 系統(tǒng)結(jié)構(gòu)原理圖(LIA,鎖相放大器;PC,電腦;DAQ,數(shù)據(jù)采集卡;FL1,2,聚焦透鏡1,2;BS1,2,分束器1,2)Fig.3.System structure schematic diagram.LIA,lock-in amplifier;PC,personal computer;DAQ,data acquisition card;FL1,2,focus lens 1,2;BS1,2,beam splitter 1,2.

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

實(shí)驗(yàn)中基于大氣中N2O 的吸收譜線,選取了Nanoplus 公司生產(chǎn)的3.9 μm ICL 作為本振光源,波數(shù)范圍為2537.5—2540.5 cm—1,激光器輸出約為1.8 mW,采用LDC501 激光控制器控制.在該波段N2O 大氣吸收占前七種大氣吸收分子的95%以上,其他氣體吸收的干擾較小.

如圖3 所示,使用上位機(jī)Labview 程序控制LDC501 進(jìn)行逐點(diǎn)步進(jìn)電流掃描,同時(shí)上位機(jī)控制程序控制太陽跟蹤儀實(shí)時(shí)采集太陽光,經(jīng)過斬波之后與激光在合束器上進(jìn)行拍頻,通過濾光片和聚焦透鏡后打到探測器上,探測器輸出的射頻信號經(jīng)過帶通濾波器、檢波器輸入到鎖相放大器,鎖相解調(diào)后,通過采集卡采集光譜數(shù)據(jù).其中,探測器的DC輸出用來監(jiān)測激光功率的穩(wěn)定性,透過半反鏡的部分本振激光用波長計(jì)進(jìn)行實(shí)時(shí)標(biāo)定.圖4 為激光外差實(shí)驗(yàn)整體結(jié)構(gòu)的實(shí)物圖.

4 實(shí)驗(yàn)結(jié)果與討論

4.1 數(shù)據(jù)分析

實(shí)驗(yàn)過程中通過頻譜密度圖(圖5 所示)選用了10—70 MHz 的帶通濾波器,避免了頻譜分布中噪聲較大的部分,在保證光譜分辨率的同時(shí)為了提高系統(tǒng)的信噪比,適當(dāng)增大了積分時(shí)間和濾波帶寬,系統(tǒng)參數(shù)設(shè)置如表1 所列.

實(shí)驗(yàn)中,本振激光使用逐點(diǎn)掃描方式控制激光器電流,系統(tǒng)的儀器函數(shù)主要由電信號響應(yīng)帶寬決定.太陽光作為一個(gè)寬帶光源用來進(jìn)行外差拍頻時(shí),在探測器的整個(gè)響應(yīng)帶寬(0—200 MHz)都有均勻的信號分布.實(shí)際測量時(shí),考慮電子噪聲和光譜分辨率的要求,需要選擇添加特定波段的帶通濾波器.使用頻譜分析儀(Agilent Technologies CXA Signal Analyzer N900A)采集并分析了系統(tǒng)的電子特性.測量結(jié)果如圖5(a)所示,其中采集了本振激光3.939 μm 的功率譜,并分別選擇加入45—98 MHz 以及10—70 MHz 帶通濾波器進(jìn)行采集,在100 MHz 左右有明顯的系統(tǒng)電子噪聲,最后選擇了噪聲較低的區(qū)域(10—70 MHz).

圖5(a) 射頻信號功率譜;(b) 儀器函數(shù)Fig.5.(a) RF signal power spectrum;(b) instrument functions.

根據(jù)激光外差原理以及頻譜儀采集的功率譜歸一化得到系統(tǒng)儀器函數(shù),如圖5(b)所示,系統(tǒng)實(shí)際雙邊帶寬為120 MHz,光譜分辨率為0.004 cm—1.

信噪比是激光外差探測系統(tǒng)的1 個(gè)重要參數(shù),可以反映系統(tǒng)的最小探測能力.儀器函數(shù)則可體現(xiàn)外差系統(tǒng)可分辨的最窄譜線寬度,所以信噪比和儀器函數(shù)是激光外差系統(tǒng)的兩個(gè)重要參數(shù),提高信噪比和獲得準(zhǔn)確的儀器函數(shù)也是反演氣體濃度廓線的必要工作.

對外差系統(tǒng)的信噪比(signal-to-noise ratio,SNR)S/N測量,1966 年Siegman[16]提出以黑體輻射光源為目標(biāo)光源的的外差系統(tǒng)的理論信噪比:

外差探測的大部分目標(biāo)光源都可以等效為一個(gè)特定溫度的黑體,所以這個(gè)信噪比公式可以簡單的估計(jì)外差裝置的信噪比等參數(shù),為激光外差系統(tǒng)的建立提供依據(jù).之后通過相關(guān)作者的研究發(fā)現(xiàn)通過鎖相放大等探測器技術(shù)可大大提高系統(tǒng)的信噪比,外差裝置的信噪比可以寫為

其中B是探測帶寬,τ是積分時(shí)間.

實(shí)驗(yàn)中首先將黑體溫度設(shè)置為1273 K,不經(jīng)過吸收直接與本振激光(3.939 μm)拍頻,記錄鎖相輸出的外差信號.實(shí)驗(yàn)其他參數(shù)保持不變,將黑體光完全遮擋,再次記錄輸出信號.兩次實(shí)驗(yàn)中,斬波頻率為75 Hz,濾波帶寬為10—70 MHz,鎖相積分時(shí)間為1 s.實(shí)驗(yàn)結(jié)果如圖6 所示,此狀態(tài)下的信噪比約為93.

圖6 外差裝置的信噪比測量Fig.6.SNR measurement of heterodyne device.

4.2 大氣透過率譜

大氣中氧化亞氮(N2O)的吸收光譜由激光外差系統(tǒng)進(jìn)行電流逐點(diǎn)步進(jìn)掃描獲得,如圖7 所示,激光逐點(diǎn)掃描范圍內(nèi)存在兩處N2O 強(qiáng)吸收,分別位于2838.336 和2539.344 cm—1處.實(shí)驗(yàn)中電流掃描間隔為0.1 mA,對應(yīng)波數(shù)變化小于0.004 cm—1,圖7(a)為激光外差系統(tǒng)測量的N2O 的外差信號,圖7(b)為ICL 波數(shù)與注入電流的關(guān)系曲線.實(shí)驗(yàn)地點(diǎn)位于31.902°N,117.167°E 的安徽省合肥市科學(xué)島.

圖7(a) N2O 吸收的外差信號;(b) 激光器波數(shù)與注入電流的關(guān)系Fig.7.(a) Heterodyne signal absorbed by N2O;(b) the relationship between laser wavenumber and injection current.

實(shí)驗(yàn)過程中,太陽光波動會對外差信號產(chǎn)生極大影響,因此實(shí)驗(yàn)過程中需要保持太陽光功率相對穩(wěn)定.對圖7 直接采集的外差信號進(jìn)行一次多項(xiàng)式擬合譜線,再做歸一化處理后得到氧化亞氮(N2O)的大氣透過率譜.由圖8 所示,實(shí)驗(yàn)測得的N2O 吸收光譜與Hirtran 數(shù)據(jù)庫的吸收光譜吻合.由測量結(jié)果(圖7(a))得到此系統(tǒng)信噪比(SNR)可以達(dá)到97,高于預(yù)期效果.實(shí)驗(yàn)所得兩個(gè)強(qiáng)吸收峰(2838.336 和2539.344 cm—1)可用于對本振激光進(jìn)行波長校正,為反演大氣氧化亞氮廓線和柱總濃度提供準(zhǔn)確的光譜數(shù)據(jù).

圖8 實(shí)驗(yàn)所得N2O 吸收光譜與模擬吸收光譜的對比Fig.8.Comparison of the N2O absorption spectrum obtained in the experiment and the simulation absorption spectrum.

4.3 反演計(jì)算

圖9 為激光外差光譜(laser heterodyne spectroscopy,LHR)數(shù)據(jù)反演的流程,由參考正向模型、實(shí)驗(yàn)數(shù)據(jù)處理和迭代反演等組成.首先讀取溫度、壓力、先驗(yàn)廓線及實(shí)驗(yàn)所得的光譜數(shù)據(jù),利用數(shù)據(jù)計(jì)算出初始透過率,然后對實(shí)驗(yàn)數(shù)據(jù)進(jìn)行預(yù)處理,包括去除太陽光波動干擾和背景信號、波長校準(zhǔn)[17]等,計(jì)算出初始Jacobian 矩陣和增益系數(shù),用于迭代計(jì)算.最后進(jìn)行循環(huán)迭代,這一步為反演算法的核心.迭代反演采用Rodger[18,19]提出的最佳估計(jì)方法(optimal estimation method,OEM).對于OEM 在基于LHR 大氣探測中的應(yīng)用,Weidmann 等[20]已經(jīng)做了詳細(xì)的介紹.在此僅簡要說明,測量數(shù)據(jù)與大氣狀態(tài)向量之間的關(guān)系被描述為

圖9 數(shù)據(jù)反演流程圖Fig.9.Data inversion flow chart.

其中ym是m維的測量向量;x是n維參考正向模型中狀態(tài)向量;ε是誤差向量,二階多項(xiàng)式基線(系數(shù)為a,b和c)表示目標(biāo)光譜區(qū)域中的非選擇性吸收和未校正的功率.

OEM 為最小化損耗函數(shù)χ2的迭代運(yùn)算過程,其中損耗函數(shù)表示為

其中,Sε是測量協(xié)方差矩陣,Sa是先驗(yàn)協(xié)方差矩陣,xa是先驗(yàn)N2O 廓線.迭代過程采用Levenberg-Marquardt (LM)方法:

其中K為雅可比矩陣 (或加權(quán)函數(shù)),γ為Levenberg-Marquardt (LM) 參數(shù).

在反演過程中,通過LM 算法迭代計(jì)算調(diào)用參考正向模型以最小化損耗函數(shù),最終反演出N2O的廓線.圖10(a)顯示了經(jīng)過預(yù)處理的激光外差實(shí)驗(yàn)數(shù)據(jù)和在2537.5—2539.5 cm—1的光譜區(qū)域內(nèi)的擬合光譜,擬合結(jié)果之間的殘差數(shù)據(jù)如圖10(b)所示,其隨機(jī)分布在±0.04 左右表明擬合效果良好.

圖10 LHR 數(shù)據(jù)反演結(jié)果 (a) 實(shí)驗(yàn)和擬合LHR 譜圖;(b) 殘差Fig.10.Inversion results of LHR data:(a) Experimental and fitted LHR spectrogram;(b) residual.

如圖11 所示,數(shù)據(jù)反演所采用的地表溫度、壓力數(shù)據(jù)來自于歐洲中期天氣預(yù)報(bào)中心(ECMWF)對于N2O 的先驗(yàn)廓線,使用了典型的中緯度冬天廓線數(shù)據(jù).大氣分層采用40 層,從地表到70 km高空,高度網(wǎng)格依次為1,2.5 和5 km.

圖11 溫度和壓力廓線Fig.11.Temperature and pressure profiles.

圖12 顯示了2019 年11 月20 日的連續(xù)測量反演結(jié)果,實(shí)驗(yàn)對中午12:00 左右的測量結(jié)果進(jìn)行了分析,這是因?yàn)闅W洲中期天氣預(yù)報(bào)中心(ECMWF)和中國氣象數(shù)據(jù)服務(wù)中心CMCC 提供的時(shí)間點(diǎn)為12:00,該段時(shí)間內(nèi)大氣參數(shù)是穩(wěn)定的.圖12 顯示了N2O 先驗(yàn)濃度廓線和反演結(jié)果所得濃度廓線,其中先驗(yàn)濃度廓線結(jié)果選取在以觀測點(diǎn) (合肥,31.9°N,117.166°E) 為中心,6°半徑范圍內(nèi)的數(shù)據(jù).在測量過程中,體積分?jǐn)?shù)標(biāo)準(zhǔn)偏差為0.000031 ×10—6—0.0026 × 10—6,對應(yīng)相對誤差范圍為0.009%—0.83%,與ECWMF 結(jié)果吻合較好,可以對大氣中的N2O 進(jìn)行長期觀測,進(jìn)而展開季節(jié)和年度的變化研究.

圖12 N2O 的先驗(yàn)和反演的垂直濃度分布圖(VMR 表示體積混合比)Fig.12.Prior and inversion vertical concentration profiles of N2O.VMR,volume mixing ratio.

5 結(jié)論

以3.939 μm 激光器作為本振光源,采用自制高精度太陽跟蹤儀采集太陽光信號,建立了一套測量整層大氣氧化亞氮(N2O)濃度廓線的中紅外激光外差探測系統(tǒng).利用光譜分辨率為0.004 cm—1的激光外差系統(tǒng)對實(shí)際大氣中氧化亞氮的透過率譜測量,得到信噪比為97 的吸收光譜.高分辨率大氣吸收光譜經(jīng)過信號預(yù)處理,利用參考正向模型以及最優(yōu)估計(jì)算法實(shí)現(xiàn)對整層N2O 濃度廓線的反演,體積分?jǐn)?shù)標(biāo)準(zhǔn)偏差為0.000031 × 10—6—0.0026 ×10—6,對應(yīng)相對誤差范圍為0.009%—0.83%.后續(xù)將對測量系統(tǒng)的噪聲、儀器函數(shù)進(jìn)行優(yōu)化處理,從而提高系統(tǒng)的信噪比,同時(shí)達(dá)到更高的分辨率,從而進(jìn)行長期觀測.