陳林菲,卞真穩(wěn)
(安徽四創(chuàng)電子股份有限公司,合肥 230088)
地基微波輻射計通過被動接收一定距離外大氣水汽和氧氣分子的輻射信息,應(yīng)用反演的方法獲取大氣參數(shù)的一種探測技術(shù)和方法。微波輻射計是被動無源探空設(shè)備,無源探測設(shè)備是大氣探測的重要組成,與常規(guī)的地基探空雷達相比,成本低,功耗低,可靠性高,且無電磁污染;與探空氣球相比,具有可以進行連續(xù)不間斷觀測的優(yōu)點,能夠?qū)ζ溆^測覆蓋區(qū)域內(nèi)的氣象要素進行常年不間斷監(jiān)測。可應(yīng)用于探測大氣云水含量、水汽、降水和大氣成分、大氣溫度廓線以及濕度廓線等重要大氣參數(shù)。同時在氣候預(yù)報模式、在數(shù)值天氣、人工影響天氣、航空飛行保障以及災(zāi)害性天氣監(jiān)測等許多領(lǐng)域都有著重要的應(yīng)用。
大氣微波遙感不但能夠彌補普通常規(guī)探測資料的不足和探測手段的應(yīng)用局限性,而且還能夠獲得很多有利于研究大氣物理結(jié)構(gòu)等方面新穎而有價值的氣象觀測資料。各種形式的星載輻射計在太陽觀測和地球大氣觀測領(lǐng)域都有著廣泛的應(yīng)用,并得到了廣泛的認可。但是星載輻射計對于對流層以下的近地面氣象目標的觀測能力受到地球大氣層的嚴重阻礙,同時氣象衛(wèi)星對于某地觀測的時空分辨率也非常低,因而對于短時天氣變化趨勢的預(yù)測力不從心。
地基微波輻射計具有對流層大氣廓線長期連續(xù)探測、功耗低、體積小、可靠性高以及可以實現(xiàn)常年無人值守等優(yōu)點,是無線探空氣球的有效補充。目前,美歐發(fā)達國家率先將微波輻射計原理用于地面氣象探測領(lǐng)域,成功開發(fā)出專用的地基微波輻射計產(chǎn)品,并成功轉(zhuǎn)入商業(yè)化。在整機系統(tǒng)上來說,我國早在上世紀七、八十年代就由北京大學(xué)完成輻射計樣機的研制,然而目前國內(nèi)在該領(lǐng)域的工程實踐與國外先進國家還是存在一定的差距。因此,本文通過繼承安徽四創(chuàng)電子股份有限公司早期開發(fā)的雙通道毫米波輻射計的工程設(shè)計經(jīng)驗和借鑒國外先進的設(shè)計理念,開展地基微波輻射計的國產(chǎn)化系統(tǒng)設(shè)計開發(fā)工作。
微波輻射計主要由天饋分系統(tǒng)、接收分系統(tǒng)、伺服分系統(tǒng)、數(shù)據(jù)采集與轉(zhuǎn)換、監(jiān)控與恒溫控制分系統(tǒng)、標定分系統(tǒng)、終端數(shù)據(jù)處理分系統(tǒng)和電源分系統(tǒng)。其系統(tǒng)組成框圖如圖1所示。
圖1 多通道毫米波輻射計系統(tǒng)組成框圖
氧分子在2.52mm及5mm(頻率60GHz屬V波段)波段有強吸收帶。如圖2所示,水汽分子在1.35cm及1.6mm波段(頻率22GHz-32GHz屬K波段),有強吸收帶。據(jù)基爾霍夫定律,這些波段吸收強,輻射也強。上層大氣發(fā)射該波段輻射向著地面的某一方向傳輸時,它一方面會受到所經(jīng)氣層的削弱,另一方面,所經(jīng)氣層也會按自身溫度放射該波段的熱輻射,僅是當各氣層溫度不同時,所放射的該波段的輻射強度也不同。地面接收到的該波段的輻射強度,與該氣層中的溫度、壓力(密度)分布有關(guān),我們就是根據(jù)地面上接收到的該波段輻射情況,來反演出大氣溫度、壓力層結(jié)情況的。式(1)是地基微波大氣遙感的基本方程,其右邊第一項是大氣層向下的輻射貢獻,第二項是宇宙背景輻射,一般Tcλ=2.7K,有時可以忽略。
對于非降水云的有云大氣,式中αλ為有云大氣的吸收系數(shù),它包括三項即 αλ=αH2O+αO2+αc,αH2O,αO2和 αc分別為水汽、氧氣和云的吸收系數(shù)。對降水云層αλ為雨云大氣的吸收系數(shù),它包括四項即
式中αr為雨的吸收系數(shù)。所以,地基遙測到輻射亮溫為頻率λ和天頂角θ,大氣的溫度T(z),壓強P(z)和水汽密度ρv(z)參數(shù)以及云雨等大氣信息的函數(shù),即知道以上參數(shù)就可以計算出亮溫的值,這個過程稱為正向或前向問題。與此相反的,用測量到的幾個頻率v或天頂角θ上的幾個值作為反演算法的輸入變量而得到感興趣的大氣參數(shù)變量的過程稱為反演過程。根據(jù)上述探測原理,微波輻射計應(yīng)能通過探測大氣亮溫,結(jié)合使用物理方法和統(tǒng)計方法建立反演方法,能長期、自動、連續(xù)的提供0km~10km的溫濕度廓線,并進而反演得出大氣柱積分水汽量、大氣柱積分云水含水量等多種應(yīng)用產(chǎn)品,滿足氣象業(yè)務(wù)和科研的實際使用需求。
圖2 典型中緯度大氣在兩個不同海拔高度下的吸收光譜
表1 系統(tǒng)主要設(shè)計指標
1.4.1 天饋系統(tǒng)設(shè)計要求
理想的天線是在天線觀測主波束內(nèi)有著固定的增益,而在主波束之外增益為零。實際的天線不僅在主波束內(nèi)有增益,在主波束外也存在增益。同時在主波束內(nèi)天線的波束形狀也不是理想的銳截止,往往存在一定的寬度。在工程上通常定義波束的-3dB點對應(yīng)的波束寬度為天線的波束寬度,即從增益最大點下降-3dB對應(yīng)的波束寬度,我們稱之為主波束。主波束之外的稱之為副瓣。天線的設(shè)計主要關(guān)注的就是波束寬度和效率。天線的主要形式有三種:喇叭天線,相控陣天線,反射面天線。其中,反射面天線是一種常用的天線類型。喇叭天線增益作為一種低增益設(shè)備一般不會應(yīng)用在較遠距離的遙感系統(tǒng)中;較為常見的是其應(yīng)用在實驗室標定、機載系統(tǒng)以及作為反射面天線的饋源。相控陣天線由于其造價高,相對損耗大,毫米波天線設(shè)計一般也不會選用。反射面天線的類型有前饋型、卡塞格倫天線(后饋式)和偏饋型。前饋型天線會導(dǎo)致系統(tǒng)有較長的波導(dǎo)走線,在毫米波波段損耗會比較大。卡塞格倫天線由于反射副面的遮擋,天線的效率會降低,副瓣抬升嚴重。偏饋型反射面天線沒有副面的遮擋,是微波輻射計的理想選擇,同時饋源直接和輻射計的接收機相連減少了系統(tǒng)的損耗;饋源采用波紋饋電喇叭,波紋饋電喇叭的反射損失小,結(jié)構(gòu)緊湊,可以提供一個較寬的帶寬、低的交叉極化電平和旋轉(zhuǎn)對稱的波束形式。
因此,微波輻射計的天饋系統(tǒng)采用偏饋型型反射面天線和波紋饋電喇叭的組合,實現(xiàn)高增益、低副瓣和旋轉(zhuǎn)對稱。具體設(shè)計上天線在輻射方向的投影直徑為250mm的圓,饋源偏饋角90°。為保證天線副瓣穿過太陽時亮溫誤差小于0.2°,拋物面反射面天線和波紋饋電喇叭的副瓣電平必須低于-30dB。同時考慮到設(shè)備的空間分辨率和結(jié)構(gòu)的緊湊性,系統(tǒng)將-3dB波束寬度定義在3°左右。天線的參數(shù)如表2所示:
表2 天線的參數(shù)
1.4.2 接收機系統(tǒng)靈敏度
天線輸出端口的負載是傳輸線和接收機,傳輸線使總能量損失,定義傳輸線輸入端功率與輸出端功率之比為損耗因子L。傳輸線自身熱輻射使功率增加(設(shè)傳輸線的溫度與天線的物理溫度相同)。接收機可等效為一個無噪聲接收機和一個位于接收機輸入端溫度為Tr的熱噪聲源。噪聲源輸出功率為kTrB。即接收機輸入端口總功率為
其中,
微波輻射計系統(tǒng)的靈敏度可以由式(4)確定。
式中,Tsys由系統(tǒng)損耗天線溫度和接收機噪聲溫度組成;ι為系統(tǒng)的積分時間。
設(shè)系統(tǒng)損耗天線溫度為230k,接收機噪聲溫度K通道為400K,V通道為700K。傳輸線損耗因子L為1.01,天線物理溫度為300K,圖3為K通道接收機靈敏度與積分時間仿真結(jié)果,圖4為V通道接收機靈敏度與積分時間仿真結(jié)果。1.4.3 反演算法選擇
圖3 K通道接收機靈敏度與積分時間的關(guān)系
圖4 V通道接收機靈敏度與積分時間的關(guān)系
溫度、濕度、云液態(tài)水廓線等的反演問題,即把地基微波輻射計所測亮溫值和外部傳感器所測參數(shù)代入網(wǎng)絡(luò)反推溫度、濕度、云液態(tài)水廓線。依據(jù)輻射計的系統(tǒng)設(shè)計,反演使用的主要方法來有:牛頓迭代統(tǒng)計回歸法;貝葉斯最大似然法;經(jīng)驗正交函數(shù)展開法;Smith迭代;線性迭代回歸法;Chedin等提出的改進的初始化反演法;神經(jīng)網(wǎng)絡(luò)反演算法。這些反演方法基本可以分為統(tǒng)計反演方法和非統(tǒng)計反演方法。統(tǒng)計反演方法是統(tǒng)計大量歷史探空資料,建立起亮溫值與反演大氣參數(shù)之間的統(tǒng)計關(guān)系。統(tǒng)計反演方法獲得的方程固定,計算速度快,但反演精度不夠;非統(tǒng)計反演是根據(jù)大氣輻射傳輸方程進行求解,解法復(fù)雜。綜合考慮我國微波輻射計溫濕廓線反演方法的實際情況和工程的可實現(xiàn)性。神經(jīng)網(wǎng)絡(luò)反演算法是一種常見且比較成熟的非線性統(tǒng)計反演方法,具有很多傳統(tǒng)方法不具備的優(yōu)點,其中最大的一個優(yōu)點是理論上可以逼近任意復(fù)雜的非線性關(guān)系,且不需要專門設(shè)計特別復(fù)雜反演算法。
四創(chuàng)公司地基多通道毫米波輻射計與2017年5月18日前,完成公司內(nèi)測試和調(diào)試工作;2017年5月18日到達安慶國家基本氣象站,截至2017年10月,共獲得168組探測數(shù)據(jù)。其中5月~6月50組數(shù)據(jù),7月依據(jù)前期數(shù)據(jù)對設(shè)備的算法進行改進,無有效數(shù)據(jù);8月份58組數(shù)據(jù),8月17、18日設(shè)備跳閘,未收集到數(shù)據(jù);9月份60組數(shù)據(jù)。初步完成8月和9月份118組輻射計反演的溫濕度廓線與安慶場站的GTS1型數(shù)字式探空儀獲取的探空數(shù)據(jù)比對分析。探空氣球受氣象條件的影響,數(shù)據(jù)在時間和空間上存在不規(guī)律的變化,而微波輻射計獲取的數(shù)據(jù)是在固定的觀測點上垂直上空56個高度層上的數(shù)據(jù)。為了使探空資料和和微波輻射計的資料可比,本文采用將探空資料的溫度和相對濕度數(shù)據(jù)進行線性插值,以得到與微波輻射計數(shù)據(jù)相同高度上的溫度和相對濕度值。分析中采用相關(guān)系數(shù)R,平均偏差Tbias和均方根誤差Tstd分別描述微波輻射計資料與探空資料的趨勢吻合度、系統(tǒng)偏差和觀測值精度,其計算公式為:
式中,TS為探空觀測值;TR為微波輻射計觀測值;N為觀測樣本數(shù)。
通過對8月份和9月份微波輻射計獲取的數(shù)據(jù)與探測氣球獲的數(shù)據(jù)分析我們得出:微波輻射計和探空氣球獲取的溫度數(shù)據(jù)相關(guān)性達到99%,平均氣溫垂直遞減率為0.57K/100m,溫度的平均誤差是-1.33°,溫度的標準差為1.87°;濕度數(shù)據(jù)的相關(guān)性為67%,相對濕度的平均誤差是-9.28%,標準差為21.14%。
針對國內(nèi)外目前地基微波輻射計的研究和開發(fā)現(xiàn)狀,本文初步提出了一種國產(chǎn)14通道地基微波輻射計的系統(tǒng)設(shè)計方法,通過對關(guān)鍵技術(shù)參數(shù)的計算和仿真分析,論證了系統(tǒng)工程化實現(xiàn)的可行性。系統(tǒng)采用多通道濾波器組并行接收技術(shù)和系統(tǒng)自動標校技術(shù),實現(xiàn)多通道毫米波輻射計系統(tǒng)對大氣溫濕廓線同時高精度探測;同時針對大氣溫濕廓線的非線性特點,終端信號處理采用BP神經(jīng)網(wǎng)絡(luò)反演算法實現(xiàn)對溫濕廓線的反演。經(jīng)過在安慶國家基本氣象站的近5個月的連續(xù)觀測試驗,驗證了設(shè)備的穩(wěn)定性和可靠性,通過與安慶場站的探空數(shù)據(jù)對比分析,初步驗證了設(shè)備的探測能力。但同時也可以看到,微波輻射計獲取的數(shù)據(jù)與探空獲取的數(shù)據(jù)還存在著一定的不一致性,后期仍需在應(yīng)用過程中通過不斷地提升設(shè)備性能、優(yōu)化反演算法來進一步完善設(shè)備。