基于可調諧二極管激光技術利用小波去噪在2.008μm波段對δ13CO2的研究?

牛明生 王貴師

1)(曲阜師范大學激光研究所,山東省激光偏光與信息技術重點實驗室,曲阜 273165)

2)(中國科學院安徽光學精密機械研究所大氣物理化學研究二室,合肥 230031)

基于可調諧二極管激光技術利用小波去噪在2.008μm波段對δ13CO2的研究?

牛明生1)2)?王貴師2)?

1)(曲阜師范大學激光研究所,山東省激光偏光與信息技術重點實驗室,曲阜 273165)

2)(中國科學院安徽光學精密機械研究所大氣物理化學研究二室,合肥 230031)

(2016年7月29日收到;2016年10月27日收到修改稿)

利用分布反饋式激光器和小型多通池建立了性能穩定的δ13CO2測量系統.基于可調諧二極管激光技術在2.008μm波段研究了幾種常見的小波評價方法、評價能力與適用性,選出最佳小波函數Haar作為小波基進行分層.在最優層上采用VisuShrink閾值函數對δ13CO2測量中的去噪效果和測量精度進行了研究.在相同實驗條件下,對去噪前后δ13CO2的測量結果進行了比較,然后從理論上分析了去噪前后測量結果不一致的原因,確定利用小波去噪對測量結果的精確性.結果表明,利用小波去噪對δ13CO2的測量精度比不用小波去噪時提高了7.3倍.

同位素測量,新型多通池,小波去噪,測量精度

1 引 言

二氧化碳(CO2)作為重要的溫室氣體之一,在控制全球氣候變化、理解生物地球化學循環與大氣碳收支方面起著關鍵作用[1,2],在環境科學和生態系統科學領域,自然界中的能量轉換過程大多會伴隨產生CO2.近50多年來,大氣中的CO2濃度一直都在增長,到2015年,全球大氣中CO2濃度已經達到了400.83 ppm[3].隨著人口增長和經濟高速發展,目前中國排放的CO2世界第一,并且在近期還可能繼續增加[4,5].為有效掌握我國CO2排放源的分布情況,制定有效的減排政策,為我國政府開展氣候變化談判和環境外交提供重要的支撐,我國對CO2排放源的研究更加迫切.CO2穩定同位素是CO2源和匯的理想示蹤劑.大氣、陸地和海洋中的物理、化學和生物交換過程的不同,導致大氣、陸地和海洋中特殊的同位素特征,可以提供參與CO2碳循環的源和匯的獨特信息[6],確定各CO2源和匯對大氣CO2濃度的貢獻.因此,為了更好地研究CO2源匯和輸送規律,了解我國碳循環收支狀況,研究CO2碳同位素比值高精度監測方法,實現大氣中CO2碳同位素比值的高精度、連續測量便顯得尤為重要.

傳統的同位素測量技術是測量精度較高的同位素比率質譜技術(IRMS),但是基于該技術制造的質譜儀結構復雜、體積龐大、成本高、需要專門的人進行維護、不能區分同質分子,測量過程程序復雜,耗時量大,很難用于實時在線連續檢測,不利于發展成便攜式儀器[7,8].近年來發展的光譜技術克服了傳統的IRMS的限制.光譜技術利用了這樣的事實,即分子的不同同位素有其自己的轉動-振動紅外光譜,根據其獨特的紅外吸收特征,可以把它與其他分子區分開來.利用光譜技術測量δ13CO2,只要選擇合適的吸收線,可以完全避開同一元素的不同同位素譜線或其他氣體譜線的干擾,且樣品不需要復雜的前期處理,使得該技術很容易發展成便攜式、可實時在線探測的儀器.應用光譜技術對同位素進行測量時,測量精度的提高一直是同位素測量中的挑戰.目前對提高同位素豐度測量精度的研究技術主要是通過設計新的吸收池或進行新的調腔技術[9-14].在這種情況下,除了系統不穩定性、樣品處理和注入吸收池的過程引入誤差外,信噪比往往是限制測量精度的主要因素.因此,在激光吸收光譜技術中,除了減小系統的固有噪聲外,可以通過選擇較強的吸收線對或高精度的光學腔來提高吸收信號強度.另一種提高信噪比的方法是通過Allan方差來優化激光器的平均掃描次數以減小檢測帶寬,使系統達到最佳平均時間.盡管這些方法在一定程度上能夠獲得較高的測量精度,但是在同位素豐度變化較快,尤其在高靈敏度和高精度的實時在線檢測中,如對呼吸氣體同位素比的檢測,系統的響應速度跟不上[15].要提高δ13CO2測量精度,需要新的去除噪聲的方法.小波變換方法在痕量氣體吸收光譜測量信號處理方面逐漸得到了應用[16,17],吸收光譜中小波變換的應用都是基于經驗方面的某一小波用來去除信號中的噪聲,以提高測量過程的信噪比.在同位素測量中,還沒有就多個小波變換進行比較,擇優應用的研究.

本文基于可調諧二極管激光技術(TDLAS),利用性能穩定的新型多通池,把小波去噪的方法引入到CO2同位素探測中.通過對小波去噪理論和實驗的研究對比,選出最佳小波函數并對其分層.在最優層數的基礎上采用適合的閾值對δ13CO2測量過程中的噪聲進行去除,以達到提高δ13CO2測量精度的目的,為今后TDLAS技術對大氣中δ13CO2的源匯進行高精度的實時在線檢測奠定了基礎.

2 實驗裝置與譜線選擇

2.1 實驗裝置

圖1為同位素測量的實驗儀器示意圖,圖1中內嵌圖是多通池通光后的實物照片,吸收池的兩反射鏡是曲率半徑為100 mm的鍍銀球面鏡,其通光孔徑為50 mm,兩鏡間的距離為12 cm,吸收池的容積為235.5 cm3.根據光斑多而不重合的原則,采取光線追跡的方法,通過理論模擬可知光在兩反射鏡間來回反射243次,其有效光程達到了29 m,經過精心的光束處理與入射光調節可以使多通池的干涉效應降至最低.Mass Controller用于控制池內氣體的流速,確保池內壓強分布均勻、穩定.Mass Controller前端三通分別連接了兩個閥門,與裝有水的瓶子相連的閥門用于緩解因氣流過大導致壓強超過設定值時的壓強.為防止水汽影響,在進氣口端連接了一個干燥器.實驗采用后端壓強控制,在氣流壓強超過設定值時,反饋給緩壓閥門,多余氣體從瓶中排出,使池內保持在壓強設定值.

圖1 同位素測量實驗儀器示意圖Fig.1.Schematic of the isotopic measurement instrument.

利用光譜對物質的相關特性進行測量,線強是溫度的函數,溫度對同位素豐度值的影響可以描述為[18,19]

式中Δδ是同位素比率變化,k為玻爾茲曼常數,T為測量樣品的絕對溫度,ΔT為測量樣品與參考氣體間的溫差,ΔE為同位素分子間吸收譜線的兩能級之間的能量差.由(1)式可知,要降低測量的δ值對溫度的依賴性,最好選取低態且能級盡量相近的吸收線對作為測量對象.

壓強是影響同位素測量精度的另一因素.池內的壓強不均勻或者有溫度梯度會引起同一譜線的線寬差異.理論上,如果不出現同位素譜線重疊的情況,且在相同濃度下,這種壓力差不會影響譜線的積分吸收值,也就不會影響測量的δ值.但實際測量時,同位素譜線可能在遠翼端重疊,這增加了對是否重疊的判斷難度.因而,為保障探測結果的精確性與準確性,實驗中要盡量消除吸收池的壓強差或壓強梯度.若吸收池內存在壓強差,需要通過下式對δ值進行修正[20]：

式中δ為測量后的修正值,δ?為實際測量值,Γsample為樣品氣體的線寬,Γref為參考氣體的線寬,γ是δ值的矯正系數,γ值可由實驗測量得到.

為保持池內壓強和溫度長期的均勻性和穩定性,吸收池放在溫控箱中.對池內壓強和溫度的控制分別采用美國生產的型號為640 B的MKS高精度壓力控制器(PC：Pressure Controller)和OMRON E5CN溫控儀,實驗參數通過計算機進行設定和連續檢測.為防止因氣體流速過大可能產生的湍流現象對測量結果的影響,采用北京七星華創公司生產的高精度的質量流量控制器對池內樣品流速進行控制.在400 sccm(1 sccm=1 mL/min)流量下,任意設定壓強和溫度分別為6.67 kPa和302.142 K,對吸收池進行控制測量,經過連續4 h的測量可知,流量變化范圍為±0.0235 sccm,標準偏差為0.0072 sccm,溫度和壓強的變化范圍分別為±11 mK和±3.2 Pa,標準偏差分別為4.33 mK和1.11 Pa,說明該同位素測量系統具有較好的溫度和壓強穩定性,完全能夠滿足同位素豐度測量的需要.

2.2 譜線選擇

同位素吸收線對的選取對于能否獲得較好的探測靈敏度與探測精度尤其重要.在利用單池進行同位素豐度測量時,一般要求同一元素的不同同位素吸收線對具有適當的間隔,與其他分子譜線不重疊,并且具有相近低態能級的兩條譜線要處于激光器一次調諧掃描范圍內,以便把溫度對測量結果的影響降至最低.因此,要想獲得較高的靈敏度和測量精度,進行譜線對選擇時通常要考慮以下因素：

1)選擇的同位素譜線對的線強應該足夠強而且應該相似或相差不大,較強的吸收線可以減小吸收池的體積,并能提高系統的信噪比,在實際測量時還能限制比爾-朗伯定律的內在非線性效應與探測器的潛在非線性效應;

2)所選同位素譜線對具有相近的低態能級,使同位素線豐度的溫度依賴性較小或溫度依賴性相同,這能夠確保吸收分子的線強隨溫度變化對不同同位素吸收線一致;

3)所選同位素吸收譜線對要處于激光器調諧掃描范圍內,小于兩個波數最佳,但是譜線對不能太近,以防掃描時不能分開,導致譜線擬合不準確;

4)在所選的同位素譜線對間隔內,不應有其他氣體分子和同位素的干擾,特別是水汽的干擾.

基于上述條件,根據HITRAN2012數據庫,我們選擇波數為4978.20475 cm-1(12CO2)和4978.02204 cm-1(13CO2)兩條吸收線.該吸收線對的頻率間隔為0.183 cm-1,吸收線強比為6.27,溫度依賴系數(Δδ/ΔT)為10.8‰ K-1,296 K時的空氣加寬系數分別為0.074和0.076 cm-1/atm.通過數據庫數據可知該吸收線對對測量條件的依賴性相對較弱,適宜于CO2穩定同位素比率的高精度測量.采用中心波長為2.008μm蝶形封裝的可調諧DFB激光器作為光源,該類激光器具有輸出功率穩定、可調諧范圍可達十幾個波數、激光控制電流與激光發射功率以及發射波長有很好的線性關系、激光發射線型的邊模抑制比高等優點.激光器的溫度和電流由ILX Lightwave生產的LDX-3724B激光控制器控制.激光器注入電流變化可以通過外部的信號發生器控制,激光器的溫度控制在33.0°C,中心電流在65 mA.實驗中選用鋸齒波信號,掃描的頻率為9 Hz,掃描電壓為1.4 V,掃描帶寬為0.47 cm-1.掃描區間涵蓋了12CO24978.20475 cm-1和13CO24978.02204 cm-1線對,而沒有其他干擾線存在.電流控制精度為0.01 mA,激光器的電流與波長的變化關系為0.017 cm-1/mA.對CO2吸收光譜進行三次多項式基線擬合,利用擬合所得的背景對光譜進行歸一化,通過計算可以得到δ13CO2.

3 結果與討論

3.1 小波去噪

同位素分析中,噪聲的存在常影響同位素豐度測量的準確度和精度,為降低噪聲提高系統的信噪比(SNR),實驗中利用Labview程序平臺編制了小波去噪程序,對一個信號進行小波包分析,采用多種小波包基,對信號進行分析對比,從中選擇最優的一種小波包基.小波分析是近十幾年來發展起來的一種數學理論和方法,在數學、物理、通訊、醫學、地質等領域獲得了極其廣泛的應用[21].目前已開發出了大量濾波技術,但只有少數,如Savitzky-Golay、傅里葉和卡爾曼濾波器被化學家廣泛使用[22].小波去噪在同位素探測中還沒有應用報道.由于小波變換具有多分辨率分析的特性,可將信號分解成一系列不同頻率的子帶,而將高頻信號與低頻信號分離,因而可實現信號的平滑和去噪.在氣體同位素豐度測量過程中,主要利用小波變換對信號進行濾波、去噪,這也是光譜分析數據處理中最基本的數據預處理環節之一.小波變換用到的函數(即小波函數)具有不惟一性,對同一個吸收信號用不同的小波函數進行分析有時結果相差很遠.正因為如此,小波函數的選取是小波分析在實際應用中的一個難點問題,也是小波分析研究的一個熱點問題.在實際應用中,我們往往通過經驗或者不斷的試驗來選擇最適合某一實際問題的小波基函數.比較常用的小波去噪方法有：模極大值去噪算法、相關性去噪算法和小波閾值去噪算法等.模極大值去噪算法采用的交替投影法,算法復雜,容易造成投影信號的偏差,難以在實際應用中對信號進行實時處理.相關性去噪算法利用小波系數在各尺度上的相關性,尤其是在信號邊緣附近,其相關性更加明顯,這種方法對信號幅度的大小不好把握;小波閾值算法比較簡單,運算量小,因此我們采用小波閾值去噪的方法對同位素測量過程中的信號進行去除.小波閾值去噪的基本思想是：染噪信號經小波變換后,信號產生的小波系數含有信號的重要信息,信號經小波分解后的小波系數較大,噪聲的小波系數較小,并且噪聲的小波系數要小于信號的小波系數,選取合適的閾值,大于閾值的小波系數被認為是由信號產生的,應予以保留,小于閾值的則認為是噪聲產生的,將其置為零,從而達到去噪的目的.基于小波閾值分析去噪的方法主要有三個步驟：1)計算含噪聲的正交小波變換,選擇合適的小波和小波分層數,將含噪信號進行小波分解,得到相應的小波分解系數,常用于去噪的小波函數有Harr小波函數[23,24]、dbn小波系[25]、symn小波系[26]、coif小波系[27]和bior小波系等[28],層數一般為3-5層;2)對分解得到的小波系數進行閾值處理,選擇適當的閾值對每一層小波系數進行量處理;3)進行小波逆變換.將經閾值處理過的小波系數重構,得到恢復的原始信號.

實驗中所采用基于Labview的小波去噪程序界面由參數選擇和光譜顯示兩部分構成.參數選擇面板中,每一個參數的選擇均采用下拉框的形式,使用方便快捷.光譜顯示面板中,不僅將不同小波去噪后的光譜顯示出來,還將去噪前后的光譜進行同框顯示,有利于對去噪效果進行實時監測和評估.同位素豐度值直接就能在界面上顯示并被采集.

利用小波進行去噪的過程中要涉及三種信號：真實信號si、原始染噪信號fi、去噪后信號gi.通常用信號的SNR與去噪后信號和真實信號的方均根誤差(RMSE)來衡量小波去噪的效果：

從(3)式可以看出,SNR并不是光譜學中能夠給出靈敏度的信噪比,我們命名為去噪信噪比,SNR越大RMSE越小去噪效果越好.以上評價標準必須在已知真實信號的前提下才能給出.然而,實驗采集的光譜,其對應的真實譜線很難獲得.因此,采用原始染噪信號相對于去噪后信號的余差的標準偏差來評價去噪效果：

由(5)式可知,在確保信號不失真的前提下,Residual-SD越大去噪效果越好,反之則差;若小波去噪沒產生任何作用,去噪后的信號與原始染噪信號會完全相同,則Residual-SD的值為0.Residual-SD對去噪結果的評價是不完備的,這是由于部分去噪后失真的信號也可能會具有較大的Residual-SD值.因而,完備的評價要同時考慮信號是否失真和Residual-SD值的大小這兩個方面的因素.在以下討論中,利用12CO2吸收峰的幅度值來判斷信號是否失真.

我們利用1000 ppm(1 ppm=1 mg/L)的CO2在室溫和40 Torr(1 Torr=133.322 Pa)的壓強下,對其吸收信號進行了測量.在測量過程中利用不同小波函數對其噪聲進行去除,如圖2所示,12CO2吸收峰的幅度值保持在0.423,根據(5)式計算得Haar小波去噪后的Residual-SD值最大為3.61985×10-5,bior2-2,coif,dbo2和sym2的Residual-SD分別為3.35733×10-5,3.46294×10-5,3.60727×10-5,3.60728×10-5. 計算結果表明,Haar小波去噪后的Residual-SD值最大,而且,12CO2吸收峰的幅度值也沒有失真.對圖2中的曲線比較也可以看出,各種小波都不同程度地濾除了噪聲,圖中矩形框內為不同小波去噪后光譜的明顯差異部分.在相同的實驗條件下,bior2-2小波函數去噪結果抑制奇異點效果較差,db02小波、coi fl小波和sym2小波去噪后曲線的平滑效果不好,相對于其他小波函數去噪效果,Haar小波去噪后曲線平滑效果要優于其他四種小波函數.因而,在這組小波函數中,Haar小波是對實驗采集光譜去噪的理想小波.

小波分解層數是小波去噪中的一個關鍵問題.理論上講,可以選取的最大分解尺度為J=|log2N|,||表示向下取整.但在實際中,一般取J為3-5.事實上,J越大,噪聲和信號表現的不同特征越明顯,越有利于信噪分離.然而,對重構來講,J越大則失真越明顯,即重構誤差越大,二者是相互矛盾的,必須選擇適當的J值,既能使信噪分離可辨又可避免重構失真.信號的最大分解尺度J應與其信噪比有關.如信噪比較大,即信號占有主要成分,則J取得稍小點就可以把噪聲分離;若信噪比較小,即噪聲占有可觀的比例,則要有效地抑制噪聲必須取較大J值.實驗表明,對于一般光譜信號,若信噪比大于20,則J取3,否則J取4或5.圖3是對選取的Haar基小波進行不同分解層數后去噪結果的對比.從圖3(a)中可以看出,隨小波分解層次增加,光譜的高頻噪聲部分被逐級去除,光譜的平滑度越來越好.但是從圖3(b)中可以看出,大于5層時,去噪后的吸收譜產生了明顯的失真,因此在進行分層時,要在最優的分解層數上進行噪聲去除.

圖2 不同小波函數的去噪結果對比Fig.2.The comparision of the denoising results with different wavelet function.

在同位素測量去噪時,確定了小波函數和最優的分解層數后,就應該選擇合適的小波閾值.不同的閾值選取將有不同的去噪效果.采用同一種小波基函數對信號進行去噪處理,那么閾值的選取將直接關系到去噪效果的優劣.如果閾值選取過小,則有一部分噪聲的小波系數將不能被置為零,這樣在去噪后的信號中保留了部分噪聲,此時信號去噪效果較差.相反地,如果閾值選取過大,則有一部分有用信號的小波系數被置零而隨噪聲一起被濾掉,使得去噪后的信號中部分有用信息丟失.目前使用的閾值分為全局閾值和局部閾值.全局閾值對各層所有小波系數或同一層內的小波系數都是統一的;局部閾值則是根據當前系數周圍的局部情況來確定適合的閾值,更具靈活性.光譜去噪時,采用局部閾值,典型的閾值估計方法有：VisuShrink(VISU)閾值、Sureshrink(SURE)閾值、Hybrid閾值和Minimax閾值.

圖3 選取Haar小波基對不同分解層數去噪結果的對比 (a)去噪后光譜平滑程度;(b)去噪前后光譜形狀的變化Fig.3.Comparison of different decomposition layers denoising results based on Haar wavelet basis：(a)The spectrum smoothness after denoising;(b)the Spectrum change before and after denoising.

VISU閾值是基于最小最大估計得出的最優閾值.閾值的選擇滿足式中σn為噪聲的均方差,N為信號的長度尺寸.這個閾值與信號的長度尺寸對數的平方根成正比.當N較大時,閾值趨向于所有小波系數置零,小波濾波器退化為低通濾波器,這在圖像去噪時造成圖像模糊,使重建圖像誤差增大.SURE閾值是一種軟件估計器,是專門針對軟閾值函數得出的結論,SURE閾值趨近于理想閾值.如果信噪比很低,SURE估計就有很大的噪聲,在這種情況下,就需要采用固定閾值形式而在高信噪比的情況下,基于SURE產生的閾值抑制噪聲的效果不明顯.此時,利用啟發函數自動在前兩種閾值選擇中選取一個較小者作為閾值.Minimax閾值是按照極大極小準則選取閾值,采用的是一種固定的閾值,它產生一個最小均方差的極值,在統計學上這種極值原理用于設計估計器.因為被去噪的信號可以看作與未知回歸函數的估計式相似,這種極值估計器可以在一個給定的函數中實現最大均方誤差最小化.Hybrid閾值也稱啟發式閾值是VISU閾值與SURE閾值的綜合.

圖4 采用不同閾值進行小波去噪結果的對比Fig.4.The comparision of the denoising results with different thresholds.

采用相同去噪條件,用以上討論的閾值法進行小波去噪.去噪后的效果如圖4所示.通過分析圖4可以明顯的看出,SURE閾值與Minimax閾值對光譜奇異點抑制效果較差,VISU閾值去噪后光譜的平滑效果優于Hybrid閾值的結果.因而,VISU閾值去噪結果最好.實驗中將采用VISU閾值對同位素測量時進行去噪.

利用上面討論的最佳去噪方法,對1000 ppm的CO2進行去噪前后的吸收信號進行對比,結果如圖5所示.通過計算對比發現,利用小波去噪后,系統信噪比是原始信號信噪比的3.85倍.因此在對CO2同位素豐度進行測量時,采用小波去噪方法對CO2同位素進行測量可以提高系統的信噪比.

圖5 吸收信號去噪前后對比Fig.5.Comparision of the absorption signals befor and after denosing.

3.2 結果與討論

利用小波去噪的方法,對濃度為400 ppm標準CO2氣體利用直接吸收法進行測量(這一濃度與實際大氣中CO2濃度基本一致).為了能更好地控制池內溫度,防止由于壓強過高使同位素吸收線展寬與其他吸收線重合,實驗在溫度為28°C和壓強為10 Torr的條件下進行,每隔2 min采集一次數據,共采集16次數據,結果如圖6所示.CO2穩定同位素的豐度值在-12.5‰左右,在沒有去噪時等效噪聲吸收系數為9.27×10-4,去噪后的等效噪聲吸收系數為1.12×10-4,同位素的測量精度前后相差7.3倍.

從圖6中可以得出,未去噪時的δ13CO2平均值為-12.50252‰,去噪后的δ13CO2平均值為-12.50032‰,其差值是未去噪時豐度的0.176‰.這是因為：在沒有進行小波濾噪時,由于噪聲對探測信號的影響,進而造成檢測靈敏度、檢測下限等性能的不準確,同時由于吸收信號的波動,其測量幅值發生變化,導致檢測穩定性下降.這就是為什么去噪前同位素豐度低的原因.而且,在測量過程中,尤其是氣體濃度很小時(這里的13CO2線強較小),各種干擾因素包括光源發光強度、環境參數和光路傳輸特性的變化以及檢測電路電子元器件的溫度漂移噪聲等的影響,對提取信號及其測量幅度帶來困難,進而導致測量精度不十分準確、穩定性差等問題.在處理具有噪聲的信號時,會把部分信號當成噪聲去掉,因此會丟失細節信號.小波閾值濾噪方法在最小均方誤差意義上有效并可達到很好的效果,它的主要理論依據是在Besov空間[29]的信號能量主要集中在幾個有限的系數中,而噪聲的能量卻分布于整個小波域中,經小波分解后信號的系數要大于噪聲的系數,因此采用閾值的辦法可以把信號的系數保留,而使大部分噪聲系數減小到零,使去噪后的測量結果更接近于真實值.

圖6 400 ppmδ13CO2去噪前后同位素豐度測量結果Fig.6.The measurement results of 400 ppmδ13CO2before and after denosing.

4 結 論

本文利用壓強和溫度都比較穩定、結構簡單、體積小的多通池TDLAS同位素測量系統,對小波去噪在同位素測量中的應用進行了研究,探索了小波去噪檢測方法及其對δ13CO2測量結果噪聲和精度的影響.在最優小波函數、小波分層和閾值函數下對400 ppm的δ13CO2進行了研究,并對小波去噪前后同位素測量精度進行了理論分析,相對于未進行去噪情況下的測量結果,利用測量的標準偏差計算[30,31],小波去噪的應用使δ13CO2的測量準確度更高,使測量精度從0.927‰提高到0.112‰,比不用小波去噪時提高了7.3倍.該實驗系統具有長時間的壓強和溫度穩定性,測得的δ13CO2精度高、速度快.小波去噪的應用對同位素豐度變化的情況,如呼吸氣體中δ13CO2的測量,能夠更加準確、快速地實時測量,具有廣闊的使用前景,為TDLAS在大氣δ13CO2監測等高精度需求領域中的應用奠定了基礎.

[1]Liu L X,Zhou L X,Vaughn B,Miller J B,Brand W A,Rothe M,Xia L J 2014J.Geophys.Res.Atmos.119 5602

[2]Mai B R,Deng X J,An X Q,Liu X T,Li F,Liu X 2014China Environmental Science34 1098(in Chinese)[麥博儒,鄧雪嬌,安興琴,劉顯通,李菲,劉霞 2014中國環境科學34 1098]

[3]Trend W 2016Nature531 281

[4]Quéré C L,Andres R J,Boden T,et al.2013Earth Syst.Sci.Data5 165

[5]Marland G 2012Nat.Clim.Change2 645

[6]Sturm P,Tuzson B,Henne S,Emmenegger L 2013Atmos.Meas.Tech.6 1659

[7]Brass M,R?ckmann T 2010Atmos.Meas.Tech.3 1707

[8]Zare R N,Kuramoto D S,Haase C,Tan S M,Crosson E R,Saad Nabil M R 2009PNAS106 10928

[9]McManus J B,Nelson D D,Zahniser M S 2015Opt.Express23 6569

[10]Sayres D S,Moyer E J,Hanisco T F,et al.2009Rev.Sci.Instrum.80 044102

[11]Wang C,Srivastava N,Jones B A,Rreese R B 2008Appl.Phys.B92 259

[12]Joseph F B,Todd B S,Max L 1992Appl.Opt.31 1921

[13]Li X X,Xu L,Gao M G,Tong J J,Jin L,Li S,Wei X L,Feng M C 2013Acta Phys.Sin.62 180203(in Chinese)[李相賢,徐亮,高閩光,童晶晶,金嶺,李勝,魏秀麗,馮明春2013物理學報62 180203]

[14]Pang J P,Wen X F,Sun X M 2016Sci.Total.Environ.539 322

[15]Bartlome R,Sigrist M W 2009Opt.Lett.34 866

[16]Li J S,Yu B L,Fischer H 2015Appl.Spectrosc.69 496

[17]Zheng C T,Ye W L,Huang J Q,Cao T S,Lv M,Dang J M,Wang Y D 2014Sens.Actuators B190 249

[18]Bergamaschi P,Schupp M,Harris G W 1994Appl.Opt.33 7704

[19]Zhang T W,Krooss B M 2001Geochim.Cosmochim.Acta65 2723

[20]Kerstel E R,Trigt R V,Dam N,Reuss J,Meijer H A J 1999Anal.Chem.71 5297

[21]Wang Y,Mo J Y 2005Spectrosc.Spect.Anal.25 124(in Chinese)[王瑛,莫金垣2005光譜學與光譜分析25 124]

[22]Wu T,Chen W D,Kerstel E,Fertein E,Gao X M,Koeth J,R??ner K,Brückner D 2010Opt.Lett.35 0146

[23]Yang Q 2011International Conference on Electronics&Optoelectronics3 129

[24]Luisier F,Vonesch C,Blu T,Unser M 2009IEEE International Symposium on Biomedical Imaging29 310

[25]Ma Y,Wang X Y,Yong H 2011CAC28 303(in Chinese)[馬毅,汪西原,雍慧 2011計算機與應用化學 28 303]

[26]Kumar H S,Pai P S,Sriram N S,Vijay G S 2013Procedia Engineering64 805

[27]Gradolewski D,Redlarski G 2014Comput.Biol.Med.52 119

[28]Naga R A,Chandralingam S,Anjaneyulu T,Satyanarayana K 2012Meas.Sci.Rev.12 46

[29]Xu J,Kawashima S 2015Arch.Ration.Mech.Anal.28 1

[30]Joseph F B,Todd B S,Max L 1992Appl.Opt.3 1921

[31]Werle P 2011Appl.Phys.B102 313

PACS:42.15.Eq,42.40.My,42.62.Fi,92.20.Xy DOI:10.7498/aps.66.024202

The research ofδ13CO2by use of wavelet de-noising at 2.008μm based on tunable diode laser absorption spectroscopy?

Niu Ming-Sheng1)2)?Wang Gui-Shi2)?

1)(Shandong Provincial Key Laboratory of Laser Polarization and Information Technology,Laser Institute,Qufu Normal University,Qufu 273165,China)

2)(Laboratory of Atmospheric Physico-Chemistry,Anhui Institute of Optics&Fine Mechanics,Chinese Academy of Sciences,Hefei 230031,China)

29 July 2016;revised manuscript

27 October 2016)

Development of optical isotope techniques has provided scientists with a set of powerful tools for investigating the sources and sink of atmospheric CO2.Here we describe a continuous,high precision,compact and portable carbon dioxide isotope ratio laser multi-pass cell spectrometer with a tunable distribute feedback laser at 2.008μm based on tunable diode laser absorption spectroscopy and,the spectrometer has good temperature and pressure stability.In order to deduce the noise,drift effect and background changes associated with low level signals,a superior signal processing technique of wavelet denoising,which possesses multi-level analytical resolutions both in time and frequency-domains,is introduced.After evaluating the method,evaluation ability and applicabilities of several common wavelet functions are analyzed and tested,the wavelet function of Haar is selected as an optimal wavelet basis function.Based on the analysis of the optimal decomposition level of Haa wavelet function,the VISU function is selected as an optimal wavelet threshold function.The denoising effect and measurement precision are evaluated by use of the VISU threshold function in the measurement process of carbon dioxide stable isotope ratio.The measurement results of carbon dioxide stable isotope ratio before and after suppressing the noises are compared in the same experiment conditions and,the inconsistent reasons of the measured results are theoretically analyzed.This technique allows the measurement of theδ-value for carbon dioxide isotopic ratios with a precision of-12.5‰and the measuremnt results show that the wavelet denoising measuring results have higher measurement accuracy,and the measurement precise of carbon dioxide isotope ratio is 7.3 times the original measurement results.The application of the wavelet denoising to the carbon dioxide isotope ratio measurement for the first time proves that the capability of the new near-infrared direct absorption technique to measure isotope ratio can permit high-frequency,near-continuous isotope measurement and obtain the high precision and accurate real-time stable isotope data directly in the field.This technique provides an important tool for studying the resource and sink of green house gases in the future.

isotope measurement,novel multi-pass cell,wavelet denoising,measurement precise

:42.15.Eq,42.40.My,42.62.Fi,92.20.Xy

10.7498/aps.66.024202

?國家自然科學基金(批準號：41405022)和曲阜師范大學博士啟動基金(批準號：20130760)資助的課題.

?通信作者.E-mail：nmsheng@163.com

?通信作者.E-mail：wulixi2004@126.com

*Project supported by the National Natural Science Foundation of China(Grant No.41405022)and the Qufu Normal University Fundation,China(Grant No.20130760).

?Corresponding author.E-mail：nmsheng@163.com

? Corresponding author.E-mail：wulixi2004@126.com