孫晨 馮玉濤 傅頔 張亞飛 李娟 劉學(xué)斌

1) (中國科學(xué)院西安光學(xué)精密機(jī)械研究所, 光譜成像技術(shù)重點(diǎn)實(shí)驗(yàn)室, 西安 710119)

2) (中國科學(xué)院大學(xué), 北京 100049)

多普勒差分干涉儀通過測量干涉圖相位變化反演大氣風(fēng)速.干涉圖信噪比是工程應(yīng)用中影響測風(fēng)干涉儀相位反演精度的核心指標(biāo)之一, 基于觀測數(shù)據(jù)分析多普勒差分干涉儀中相位不確定度與干涉圖信噪比之間的定量關(guān)系建立解析模型, 對儀器設(shè)計、性能評價及風(fēng)場數(shù)據(jù)應(yīng)用具有重要意義.本文基于傅里葉變換光譜學(xué)噪聲傳遞理論和多普勒差分干涉儀相位反演模型, 建立干涉圖信噪比與相位不確定度傳遞的理論模型.利用多組不同信噪比仿真及實(shí)驗(yàn)實(shí)測數(shù)據(jù), 將相位反演模型計算的統(tǒng)計結(jié)果與本文提出的理論模型分析結(jié)果進(jìn)行比較驗(yàn)證, 結(jié)果表明, 多組實(shí)驗(yàn)數(shù)據(jù)平均殘差小于0.07 mrad, 兩者的一致性較好, 建立的模型為多普勒差分干涉儀效能評估、干涉儀設(shè)計優(yōu)化提供理論依據(jù).

1 引 言

大氣風(fēng)場是描述空間環(huán)境特征的重要物理參數(shù), 是空間天氣與空間環(huán)境學(xué)研究領(lǐng)域的熱點(diǎn)課題, 也是空間科學(xué)中最有應(yīng)用價值的研究內(nèi)容之一[1?3].星載測風(fēng)干涉儀獲取中高層大氣風(fēng)場廓線的長期觀測數(shù)據(jù), 對空間天氣預(yù)報、航天器發(fā)射、運(yùn)行和返回軌道參數(shù)計算、無線電通信、空間科學(xué)實(shí)驗(yàn)等都具有重要意義, 為研究上下層大氣之間的能量傳輸、大氣成分變化以及大氣波現(xiàn)象提供重要理論依據(jù)[4?6].

多普勒差分干涉光譜技術(shù)是由Englert團(tuán)隊(duì)在2006年提出的一種新型的中高層大氣風(fēng)場探測技術(shù)[6?9], 多普勒差分干涉儀是通過探測氣輝譜線相位變化量反演多普勒頻移, 進(jìn)而計算大氣風(fēng)速,相位反演精度決定風(fēng)速測量精度[10?15].在實(shí)際應(yīng)用中, 儀器采集到的原始干涉圖會存在光電轉(zhuǎn)換、探測器暗電流、讀出電路帶來的各種隨機(jī)噪聲, 這些噪聲會隨相位反演模型傳入相位, 導(dǎo)致反演相位隨機(jī)誤差.在儀器工程設(shè)計中通常采用原始干涉圖信噪比作為衡量儀器信號質(zhì)量的重要指標(biāo), 是系統(tǒng)基本光學(xué)、電子學(xué)參數(shù)確定的依據(jù)[16?20].因此, 研究多普勒差分干涉儀干涉圖信噪比與反演相位不確定度之間的傳遞關(guān)系, 建立定量數(shù)理模型對指導(dǎo)儀器工程設(shè)計是非常重要的.

傳統(tǒng)傅里葉變換光譜儀以反演光譜為目的, 已建立了從原始干涉圖到復(fù)原光譜噪聲傳遞關(guān)系, 但由于多普勒差分干涉儀目標(biāo)是測量多普勒頻移帶來的相位漂移, 故已有的傅里葉變換型光譜儀數(shù)據(jù)反演的噪聲傳遞理論模型不能直接用于衡量多普勒差分干涉儀相位不確定度與原始干涉圖信噪比的關(guān)系.Luo[17]分析了基于對稱式的空間外差干涉光譜儀中以光子噪聲為主要噪聲源時的干涉域信噪比表達(dá)式, 提出空域與頻域信噪比的數(shù)理模型,但無法說明原始干涉圖信噪比與相位不確定度之間的關(guān)系.Pritt[19]針對傳統(tǒng)的傅里葉變換光譜儀,建立原始干涉圖噪聲與反演光譜噪聲之間的傳遞模型, 但未進(jìn)一步證明其與復(fù)干涉圖噪聲和相位噪聲之間的傳遞關(guān)系.

本文基于傅里葉變換光譜學(xué)基礎(chǔ)理論, 根據(jù)多普勒差分干涉儀的相位反演方法, 結(jié)合原始干涉圖表達(dá)式、系統(tǒng)噪聲模型及相位反演模型, 建立原始干涉圖信噪比與相位不確定度傳遞的理論模型, 并利用實(shí)驗(yàn)室研制的多普勒差分干涉儀驗(yàn)證其正確性, 分析誤差來源, 為儀器的工程設(shè)計與性能評價提供重要理論基礎(chǔ).

2 理論推導(dǎo)分析

2.1 多普勒差分干涉儀原理

多普勒差分干涉儀是將邁克爾遜干涉儀中兩干涉臂末端的平面反射鏡替換成閃耀光柵, 在一個干涉臂中加入非對稱量 ?d, 增大光程差, 提高干涉儀風(fēng)速測量靈敏度.如圖1所示.當(dāng)任意單色光譜輸入B(σ) 進(jìn)入儀器入瞳, 經(jīng)過準(zhǔn)直鏡、濾光片和分束器后形成兩束光, 并分別在兩個干涉臂末端的閃耀光柵處發(fā)生衍射, 經(jīng)過閃耀光柵后返回, 兩束光的波陣面在分束器處相遇, 經(jīng)過條紋成像系統(tǒng)L2、L3, 在探測器上形成的干涉條紋I滿足[21]:

圖1 多普勒差分干涉儀原理圖Fig.1.Schematic of the Doppler asymmetric spatial heterodyne interferometer.

式中,σ為入射光波數(shù),σL為 Littrow 波數(shù),θL為Littrow角,x為探測器上位置坐標(biāo), 定義探測器中心x=0 , ?d為兩干涉臂非對稱量.

2.2 相位反演模型

多普勒差分干涉儀通過探測干涉圖的相位變化反演氣輝譜線的多普勒頻移, 相位反演模型如圖2所示, 按照參考文獻(xiàn)[6], 將式(1)簡化為

式中,j為通過濾光片進(jìn)入儀器的譜線,Sj為譜線強(qiáng)度系數(shù),Ej(x) 為第j根譜線的包絡(luò)函數(shù),kj≡4(σj?σL)tanθL為差分頻率,?j≡4π(σj?σL)?d為附加相位項(xiàng), δφj為多普勒頻移引起的相移項(xiàng).對(2)式傅里葉變換可得空間頻率在 +kj和 ?kj有特征峰的反演光譜S, 用窗函數(shù)分離S中的一個特征峰 ( +kj或者 ?kj), 其余的所有元素均賦 0, 傅里葉逆變換后可得:

此時可計算任一采樣點(diǎn)對應(yīng)的相位值:

圖2 相位反演模型流程圖Fig.2.Flow chart of phase retrieval method.

2.3 原始干涉圖信噪比與相位不確定度理論推導(dǎo)分析

在多普勒差分干涉儀成像過程中, 噪聲主要來源于電路噪聲(讀出噪聲、暗電流噪聲)和散粒噪聲中的光子噪聲(屬于高斯白噪聲).因此疊加噪聲的干涉圖I′可以表示為

式中εN表示電路讀出噪聲、暗電流噪聲和光子噪聲引起的原始干涉圖強(qiáng)度的隨機(jī)變化, 由隨機(jī)噪聲特性可知, 數(shù)學(xué)期望為0, 且相鄰空間位置不相關(guān).假設(shè)原始干涉圖以光子噪聲為主要噪聲源, 光子噪聲方差與探測器光電轉(zhuǎn)換后生成的信號對應(yīng)電子數(shù)相等[16].分別對(5)式中I′、I進(jìn)行傅里葉變換、窗函數(shù)單頻提取、傅里葉逆變換后得、Ic, 結(jié)合隨機(jī)信號分析方法, 用(4)式求得、Ic對應(yīng)的反演相位 P hase′、 P hase , 兩相位差 ? Phase 即為原始干涉圖中隨機(jī)噪聲引入的相位噪聲:

由于

因此, 用反正切函數(shù)的一階泰勒展開式求解(6)式,即:

求 ? Phase 標(biāo)準(zhǔn)差得到由原始干涉圖光子噪聲引入的相位不確定度:

式中,n為窗函數(shù)尺寸,N為干涉圖采樣數(shù), 將(10)式代入(9)式中可得:

原始干涉圖任一采樣點(diǎn)信噪比表達(dá)式為

式中εN和分別對應(yīng)任一采樣點(diǎn)的光子噪聲和原始干涉圖信號.將(12)式代入(11)式可得原始干涉圖信噪比與相位不確定度之間的解析表達(dá)式:

3 計算機(jī)仿真分析

圖3 信噪比與相位不確定度仿真結(jié)果及殘差 (a)為相位反演模型與理論模型計算的相位不確定度; (b)為兩種計算結(jié)果殘差Fig.3.Simulation results of retrieved phase uncertainty and residual with different SNR: (a) Phase uncertainty of two models; (b) the differences between two results.

為充分驗(yàn)證推導(dǎo)的理論模型的正確性, 本節(jié)保證仿真干涉圖與實(shí)驗(yàn)室實(shí)測干涉圖有完全相同的邊界條件, 使用計算機(jī)仿真生成信噪比與實(shí)測信噪比相同的干涉圖, 每組數(shù)據(jù) 1000 幀.首先, 采用第2節(jié)中的相位反演模型, 去除原始干涉圖的低頻基線, 傅里葉變換用長度為5的矩形窗函數(shù)對反演譜單頻提取, 經(jīng)傅里葉逆變換后, 反演每幀干涉圖相位, 取第512個采樣點(diǎn)相位統(tǒng)計標(biāo)準(zhǔn)差作為每組干涉圖的相位不確定度; 其次, 再用本文推導(dǎo)的(13)式, 根據(jù)每幀干涉圖信噪比計算相位不確定度, 統(tǒng)計平均后得到每組數(shù)據(jù)對應(yīng)的結(jié)果.兩種方法計算結(jié)果如圖3所示, 其中右圖為兩方法計算結(jié)果殘差, 平均殘差為 0.0239 mrad.造成殘差的原因主要是在理論推導(dǎo)中為了方便計算, 用反正切函數(shù)的泰勒一階展開式對相位方差近似計算.從圖中可看出, 隨著信噪比的提高, 相位不確定度與殘差均減小, 兩種方法的計算結(jié)果一致性較好.

4 實(shí)驗(yàn)室測量分析

為了進(jìn)一步驗(yàn)證本文推導(dǎo)的原始干涉圖信噪比與相位不確定度理論模型正確性, 結(jié)合實(shí)驗(yàn)室自研的多普勒差分干涉儀(主要參數(shù)參見表1)搭建實(shí)驗(yàn)平臺, 實(shí)驗(yàn)器材包括多普勒差分干涉儀、波長為632.8 nm的穩(wěn)頻激光器、分光鏡、衰減片、準(zhǔn)直鏡, 如圖4 所示.

圖4 實(shí)驗(yàn)平臺Fig.4.Photograph of laboratory setup.

為了提高系統(tǒng)的穩(wěn)定性, 降低相位熱漂移導(dǎo)致的系統(tǒng)性偏差, 對多普勒差分干涉儀進(jìn)行被動隔熱, 并對干涉儀組件主動溫控, 精度為 0.1 ℃, 探測器工作溫度為–75 ℃, 并且在采集記錄數(shù)據(jù)之前先使整個系統(tǒng)工作數(shù)小時到達(dá)穩(wěn)定狀態(tài), 圖5為積分時間為0.5 s時的原始干涉圖.

實(shí)驗(yàn)中用兩種方法分別采集多組不同信噪比干涉圖數(shù)據(jù), 第一種方法使用同一衰減片, 調(diào)整積分時間控制原始干涉圖強(qiáng)度, 獲得不同信噪比干涉圖; 第二種方法使用相同的積分時間, 調(diào)整衰減片獲得不同信噪比干涉圖.這樣設(shè)計主要有以下兩種原因: 首先, 驗(yàn)證理論模型對不同噪聲類型是否具有普適性; 其次, 為了使得實(shí)測干涉圖信噪比覆蓋更廣.由文獻(xiàn)[16,17]可知, 實(shí)測干涉圖噪聲主要包含光子噪聲、暗電流噪聲、讀出噪聲, 因此當(dāng)改變積分時間時, 三種噪聲均發(fā)生改變, 當(dāng)調(diào)整衰減片強(qiáng)度時, 僅光子噪聲發(fā)生改變.使用第一種方法選擇兩組不同透過率的衰減片, 其中第一組衰減片,積分時間起始為 0.2 s, 間隔為 0.2 s; 第二組衰減片, 積分時間起始為 0.5 s, 間隔為 0.5 s.每組衰減片實(shí)驗(yàn)中各采集9組數(shù)據(jù), 每組100幀.由于實(shí)驗(yàn)中沒有精確測量入射光強(qiáng), 因此選擇多幀平均方法計算實(shí)測干涉圖信噪比: 將每組100幀干涉圖平均得到平均干涉圖, 每幀干涉圖與其偏差作為每幀干涉圖噪聲, 進(jìn)而計算干涉圖信噪比.

圖5 原始干涉圖Fig.5.Interferogram.

表1 多普勒差分干涉儀主要參數(shù)Table 1.Principle components.

在用相位反演模型和本文推導(dǎo)的理論模型計算之前, 對原始干涉圖預(yù)處理(去除沖擊噪聲和暗電平).與仿真計算方法相同, 首先用第2節(jié)中的方法計算每幀干涉圖相位, 得出每組干涉圖的相位不確定度, 整個計算方法與前文計算仿真數(shù)據(jù)的數(shù)據(jù)處理方法相同.其次, 用本文推導(dǎo)的理論模型計算每幀干涉圖的理論相位不確定度, 再對每組結(jié)果統(tǒng)計平均得到每組干涉圖的相位不確定度.圖6分別為兩組衰減片對應(yīng)的不同信噪比干涉圖用相位反演模型與本文推導(dǎo)的理論模型計算結(jié)果及殘差.圖7為相同積分時間, 不同衰減片對應(yīng)的實(shí)測數(shù)據(jù)計算結(jié)果及殘差.

圖6中實(shí)線圈, 第3組數(shù)據(jù)為積分時間0.2 s的干涉圖, 平均信噪比為73.488倍, 反演模型計算的不確定度為1.2213 mrad, 理論模型計算結(jié)果為1.0552 mrad, 仿真生成干涉圖信噪比為72.8292倍, 相位反演模型計算結(jié)果為 1.0396 mrad, 理論模型計算結(jié)果為1.0075 mrad.殘差偏大的原因是由于實(shí)測干涉圖相位漂移較大, 圖8中黑色為第3組實(shí)測數(shù)據(jù)反演相位差, 紅色為仿真數(shù)據(jù)相同信噪比反演相位差, 從圖中可以明顯看出, 實(shí)測數(shù)據(jù)反演相位差的方差大于仿真數(shù)據(jù).

圖6中虛線圈, 第4組數(shù)據(jù)為積分時間1.5 s的干涉圖, 平均信噪比為74.286倍, 反演模型計算的不確定度為1.1006 mrad, 理論模型計算結(jié)果為0.8663 mrad, 仿真生成的干涉圖信噪比為72.8292倍時, 相位反演模型計算結(jié)果為 1.0396 mrad, 理論模型計算結(jié)果為 1.0075 mrad.分析發(fā)現(xiàn), 整組實(shí)測數(shù)據(jù)的信噪比波動較大, 如圖9所示, 由于波動較大, 因此在平均計算時出現(xiàn)偏差, 造成計算結(jié)果殘差較大.

圖6 變積分時間實(shí)測數(shù)據(jù)計算結(jié)果及殘差 (a)為相位反演模型與理論模型計算的相位不確定度; (b)為兩種計算結(jié)果殘差Fig.6.Experimental results of retrieved phase uncertainty and residual with same attenuation coefficient and different integration time: (a) Phase uncertainty of two models; (b) the differences between two results.

圖7 變衰減片實(shí)測數(shù)據(jù)計算結(jié)果及殘差 (a)為相位反演模型與理論模型計算的相位不確定度; (b)為兩種計算結(jié)果殘差Fig.7.Experimental results of retrieved phase uncertainty and residual with different attenuation coefficients: (a) Phase uncertainty of two models; (b) the differences between two results.

圖8 相同信噪比實(shí)測數(shù)據(jù)與仿真數(shù)據(jù)對比分析Fig.8.Comparison and analysis of measured and simulated data with the same SNR.

圖9 積分時間 1.5 s實(shí)測干涉圖信噪比Fig.9.SNR of 1.5 s integration time interferograms.

分析圖6結(jié)果可得, 由于積分時間相對較短,因此在使用不同方法采集干涉圖對最終的計算結(jié)果影響甚小.由圖6、圖7結(jié)果可得, 本文推導(dǎo)的理論模型與相位反演模型計算得到的相位不確定度一致性較好.圖6中兩組不同衰減片計算結(jié)果平均殘差分別為0.0634 mrad, 圖7計算結(jié)果平均殘差為 0.0573 mrad, 對應(yīng)風(fēng)速誤差分別為 0.0373 m/s,0.0337 m/s.(以波長為 632.8 nm的激光作為輸入光源, 實(shí)驗(yàn)室自研的多普勒差分干涉儀基礎(chǔ)光程差為 50 mm, 1 m/s的風(fēng)速變化量對應(yīng)的相位變化量為1.7 mrad).兩種方法的計算結(jié)果的一致性表明本文推導(dǎo)的理論模型正確.

對比圖6、圖7可發(fā)現(xiàn), 在相同條件下, 兩種計算方法的殘差主要有以下幾點(diǎn): 首先, 在實(shí)際測量中, 隨著積分時間的加長, 使得整組干涉圖的采集時間加長, 信噪比波動范圍變大, 如圖6中當(dāng)信噪比大于 200 倍時, 積分時間為 4 s和 4.5 s的兩組實(shí)測數(shù)據(jù), 兩種方法的計算結(jié)果殘差變大; 其次,實(shí)驗(yàn)中對多普勒差分干涉儀中的成像系統(tǒng)被動隔熱, 并未主動溫控, 可能會對實(shí)測干涉圖造成相位漂移.后續(xù)研究中可將成像系統(tǒng)放大率帶來的相位漂移修正后, 進(jìn)一步提高反演相位精度.

5 結(jié) 論

基于多普勒差分干涉方程, 結(jié)合相位反演模型, 對干涉圖進(jìn)行傅里葉變換、窗函數(shù)單頻提取、逆傅里葉變換、相位計算、反正切函數(shù)的一階泰勒展開, 分步解析多普勒差分干涉儀中相位不確定度與干涉圖信噪比之間的定量關(guān)系.分別用多組不同信噪比的仿真與實(shí)測數(shù)據(jù)驗(yàn)證本文推導(dǎo)的理論模型, 與相位反演模型計算結(jié)果的高度一致性驗(yàn)證了本文提出模型的正確性.本文成果在工程設(shè)計中結(jié)合目標(biāo)源輻射特性分析、光學(xué)設(shè)計參數(shù)及電子學(xué)參數(shù)能夠通過確定原始干涉圖信噪比, 進(jìn)而得出相位不確定度, 實(shí)現(xiàn)對多普勒差分干涉儀全鏈路觀測能力分析的目的.除此之外基于觀測數(shù)據(jù)的信噪比分析, 還能夠?qū)崿F(xiàn)對實(shí)測觀測數(shù)據(jù)的誤差限分析, 為數(shù)據(jù)后續(xù)的科學(xué)應(yīng)用提供不確定度輸入.

感謝中國科學(xué)院西安光學(xué)精密機(jī)械研究所馮玉濤研究員和傅頔博士的討論.