基于非線性最小二乘法的多光譜擬合程序： 應用于12CH4譜線參數分析

馬宏亮，鄭健捷，劉 強，錢仙妹，朱文越

1. 中國科學院合肥物質科學研究院安徽光學精密機械研究所中國科學院大氣光學重點實驗室，安徽 合肥 230031 2. 安慶師范大學電子工程與智能制造學院，安徽 安慶 246133 3. 先進激光技術安徽省重點實驗室，安徽 合肥 230037 4. 中國科學技術大學研究生院科學島分院，安徽 合肥 230026

引 言

準確測量大氣中甲烷氣體濃度的時空分布對于研究溫室效應以及氣候變化等具有十分重要的意義[1-2]。鑒于此，人們發(fā)展了一系列的技術手段來實現這一目的，如基于星載或者地基的甲烷氣體濃度廓線探測系統[3]。由這些探測系統可以直接觀測獲得多組大氣透過率光譜，然后，利用已知的甲烷氣體分子譜線參數(譜線的位置、強度、碰撞展寬以及它們與溫度的依賴關系等)可以獲得所需要的濃度廓線信息。為了準確獲得這些譜線參數，需要在已知實驗條件(濃度，溫度，壓強以及分子種類的混合比等)的情況下，多次掃描同一波段范圍得到多組實驗測量光譜，然后利用特定的程序擬合這些光譜反演獲得。對于一次只能處理一張光譜的擬合程序稱為單光譜擬合程序。采用單光譜擬合程序逐張擬合上述光譜，并把得到的數據進行平均或者線性擬合即可獲得相應的譜線參數。然而，單光譜擬合程序處理實驗光譜最大的缺點是非常耗時，同時，容易引起擬合過程中的誤差傳遞[4]。為克服以上問題，國際上多個研究小組發(fā)展出了多光譜擬合程序[4-6]。然而，上述多光譜擬合程序都是針對分析由傅里葉光譜儀得到的光譜開發(fā)的，目前很少有針對分析由TDLAS所獲得的光譜而開發(fā)的此類程序。

本文利用最小二乘法開發(fā)了一款適用于處理TDLAS所測光譜的多光譜擬合程序，文中對該程序的功能做了詳細的介紹； 然后基于此程序分析了12CH4分子在2 958～2 959 cm-1波段內多個不同溫度下的吸收光譜，獲得了6條譜線的空氣展寬系數，并與之前采用單光譜擬合程序獲得的相關數值作了比較，驗證了程序的穩(wěn)定性。

1 多光譜擬合程序

1.1 基于非線性最小二乘法的光譜擬合程序

非線性最小二乘擬合(NLSF)技術廣泛用于實驗室或實際大氣測量光譜的分析工作[5]。當處理給定的單張實驗光譜時，非線性最小二乘方法擬合函數形式為

(1)

如果同時處理多組給定的實驗光譜時，只需對式(1)稍作修改即可獲得多光譜非線性最小二乘擬合(NLSF)函數形式為

(2)

1.2 多光譜擬合程序的結構和功能

多光譜擬合程序的結構如圖1所示，分為輸入、主程序、輸出三個功能不同的模塊，下面分別就每個模塊的功能和用法作相關介紹。

1.2.1 輸入模塊

如圖1所示，輸入模塊一共包括3個部分。方框1包含需要擬合的實驗光譜數據，一般為txt格式文件，每個文件對應一張實驗光譜，里面包含兩列ASCII原始實驗數據，第一列為采樣點的頻率，第二列為采樣點的信號強度(即待擬合數據)，文件的第一行需要包括采樣點數。擬合前需要給每張光譜文件命名，把所有文件與執(zhí)行程序放在一個文件夾里。方框2為待擬合光譜文件名列表，里面包含方框1中需要擬合的文件名稱(即方框1中命名的文件)。方框3為擬合模型的原始輸入參數，各參數的名稱和意義如表1所示。

圖1 多光譜擬合程序流程圖Fig.1 Simplified flow diagram of the multispectrum fitting program

表1 模型需要輸入的原始參數及其意義Table 1 The initial parameter names and meanings of the model

1.2.2 主程序模塊

如圖1所示，左邊為主程序的數據處理過程，其詳細功能如表2所示。

表2 主程序數據處理過程及意義Table 2 The data processing and significance of the main program

如圖1所示，右邊為主程序的執(zhí)行過程，開始擬合時，程序會依次提示如表3所示的5步操作。表3給出了各操作步驟的名稱和意義。

表3 主程序執(zhí)行步驟及各步驟意義Table 3 The execution steps and significance of the main program

1.2.3 輸出模塊

如圖1所示，輸出模塊共有3部分。方框4為方框1中單張光譜經過擬合之后相對應的輸出數據。方框5為采用全局擬合并且歸一化后得到的所有光譜的透過率譜圖。方框6為對所有光譜擬合后得到的譜線參數，如基線、譜線位置、多普勒展寬、壓力展寬、吸收面積、擬合殘差。

2 結果與討論

本文利用多光譜擬合程序處理了2 958～2 959 cm-1波數內12CH4分子6條躍遷譜線5組不同溫度下的空氣(Air)展寬實驗光譜用以驗證其性能。測量這些光譜的有關實驗設備細節(jié)和詳細的壓強、溫度等實驗條件可以參考文獻[7]。

圖2(a)為采用多光譜擬合程序同時對4張溫度為173 K, 壓力為18.158 kPa, Air+12CH4的吸收光譜進行Voigt線型擬合所獲得的全局擬合輸出結果, 擬合過程中多普勒展寬固定在其理論值，圖2(b)為對應的擬合殘差。由圖2(b)中的擬合殘差可以看出多光譜擬合程序可以較好的完成多張光譜數據同時擬合的任務。

圖2 采用多光譜擬合程序同時擬合4張光譜(a)： 測量(散點線)和計算(實線)得到的光譜；(b)： 相應的擬合殘差Fig.2 Example of multispectrum fitting program operation of with four spectra(a)： The measured (in scatter line) and calculated (in solid line)line profiles with a multi-peak fitting program； (b)： The corresponding fitting residuals

表4 多光譜擬合程序獲得的12CH4空氣展寬系數及其與文獻[7]相應數值的對比結果Table 4 The air-broadening coefficients of 12CH4 derived from the multispectrum fitting program and the difference between them and those from Ref. [7]

3 結 論

基于最小二乘法和Levenberg-Marquardt迭代算法開發(fā)了一款可以同時擬合多張光譜并獲得一套光譜參數的多光譜擬合程序，該程序適用于處理由TDLAS所獲得的吸收光譜。利用該程序處理了2 958～2 959 cm-1波數內12CH4分子6條躍遷譜線的實驗光譜，獲得了這些譜線相應的空氣展寬系數。通過比較發(fā)現，利用多光譜擬合程序獲得的空氣展寬系數與之前文獻報道的結果整體符合較好，說明該程序具有較好的可靠性。與單光譜擬合程序相比，該程序可以極大的節(jié)省實驗數據處理時間，并且可以減小數據處理過程中的誤差傳遞。

致謝： 作者感謝中國科學技術大學劉安雯副教授在多光譜擬合程序編寫以及論文撰寫過程中給予的指導。