劉清惓，施泓伊，楊 杰

（1.南京信息工程大學(xué) 電子與信息工程學(xué)院，江蘇 南京 210044；2.江蘇省大氣環(huán)境與裝備技術(shù)協(xié)同創(chuàng)新中心，江蘇 南京 210044；3.江蘇省氣象探測與信息處理重點實驗室，江蘇 南京 210044）

高空氣象探測通常利用氣象飛機、無線電探空儀、氣象雷達(dá)等儀器設(shè)備對各個高度上的氣溫、濕度進(jìn)行檢測，大氣溫度觀測數(shù)據(jù)的精確度對于天氣預(yù)報、大氣實時監(jiān)測、氣象災(zāi)害預(yù)測和預(yù)警至關(guān)重要。在高空氣溫測量中，探空溫度傳感器的測量精度會受到不同環(huán)境因素的影響，其中引起測量誤差的主要因素是太陽輻射。為了讓溫度傳感器測量元件表面有足夠的通風(fēng)量，設(shè)計者直接將溫度傳感器放置在空氣中，導(dǎo)致溫度傳感器吸收到的太陽輻射疊加在測溫元件上，從而造成太陽輻射誤差。探空氣溫資料顯示，1958—2017年氣溫在對流層（850～150 hPa）總體上升趨勢為0.03～0.15 ℃/（10 a），平流層下層下降趨勢為0.18 ℃/（10 a）；1979—2017年對流層升溫趨勢為0.1～0.25 ℃/（10 a），平流層下層下降趨勢為0.1 ℃/（10 a）。然而由于輻射誤差的存在，溫度傳感器在高空的輻射誤差最高可達(dá)3.3 ℃量級，難以滿足高空溫度觀測精度的需求。因此，探空溫度傳感器的測量精度需要進(jìn)一步提高。

根據(jù)國際氣象組織WMO 提供的資料顯示，各類電子溫度傳感器在太陽高度角為20°時，輻射誤差最高為1.8～5.1 ℃。目前國內(nèi)使用的GTS1 探空儀搭載的溫度傳感器選擇棒狀熱敏電阻，表面涂層對短波有較高的反射率，但長波吸收率達(dá)到90%，造成的輻射誤差約為3 K。我國的長峰探空儀搭載的珠狀溫度傳感器整體性能比較穩(wěn)定，其誤差在海拔高度15 km 以上，最高可達(dá)-1 ℃。綜上所述，傳統(tǒng)的溫度傳感器未能達(dá)到高空溫度觀測的高精度要求。

為了降低高空太陽輻射誤差，本文利用壓電陶瓷彎曲振動的特性，設(shè)計一種基于熱電偶的探空溫度傳感器。該傳感器能夠在氣流速度較低時增加溫度傳感器表面的對流。通過計算流體動力學(xué)（CFD）對其進(jìn)行仿真分析，使用支持向量機（Support Vector Machine，SVM）對仿真結(jié)果進(jìn)行誤差修正，并搭建模擬高空環(huán)境的實驗平臺驗證傳感器的測量精度。

1 流體動力學(xué)仿真分析

1.1 探空溫度傳感器探頭設(shè)計

微型熱電偶的尺寸小，雷諾系數(shù)較低，可以準(zhǔn)確地測量被測物體的溫度，相對于珠狀熱敏電阻，能有效降低太陽輻射測量誤差。本文使用的微型熱電偶直徑為0.13 mm，選擇穩(wěn)定性高、精度高的鉑電阻對熱電偶的冷端進(jìn)行補償。銅球表面反射率高，受太陽輻射影響低，因此把鉑電阻和熱電偶的冷端使用銅球裹住，可以有效減小鉑電阻的太陽輻射誤差。在高空中，當(dāng)傳感器升空速度較慢時，探空溫度傳感器進(jìn)行溫度檢測時受到的太陽輻射影響較大，因此本文利用壓電陶瓷的逆壓電效應(yīng)來避免這種情況。壓電陶瓷是鋯鈦酸鉛（PZT）化合物，將壓電陶瓷晶體接入一個與極化方向相同（或相反）的電場中，在電場的激勵下，晶體會沿極化方向（或相反方向）發(fā)生不同方向的形變，這種把電能轉(zhuǎn)變?yōu)闄C械能的現(xiàn)象稱作逆壓電效應(yīng)。將兩片長30 mm、寬10 mm、厚0.3 mm 的PZT 壓電陶瓷片分別黏合在一片長50 mm、寬10 mm、厚0.3 mm 銅片的正反面，以串聯(lián)的方式接入電源，制成一個壓電振子。在電場的激勵下，壓電陶瓷片會產(chǎn)生彎曲振動。將壓電陶瓷片連接在熱電偶的熱端，在高空中氣流速度較低的情況下，陶瓷片通電產(chǎn)生彎曲振動，增加熱電偶表面通風(fēng)量，提高散熱性，從而提高測量精度。熱電偶冷端結(jié)構(gòu)如圖1 所示。

圖1 熱電偶冷端結(jié)構(gòu)示意圖

1.2 建立物理模型及網(wǎng)格劃分

使用CFD 方法對探空溫度傳感器探頭部分建立模型，并在外部建立一個長方體充當(dāng)空氣域，探頭部分包括壓電陶瓷振子和熱電偶，其物理模型如圖2 所示。

圖2 探頭和外圍空氣域的物理模型

ICEM 擁有相當(dāng)靈活的網(wǎng)格劃分功能，文中使用ICEM 對壓電陶瓷和熱電偶以及長方體空氣域進(jìn)行非結(jié)構(gòu)化網(wǎng)格劃分。為減小誤差，對網(wǎng)格數(shù)100 萬～150 萬的模型進(jìn)行測驗。當(dāng)網(wǎng)格數(shù)為130 萬時，質(zhì)量最好，因此選取130 萬網(wǎng)格對傳感器進(jìn)行迭代計算。

1.3 CFD 仿真結(jié)果分析

使用CFD 仿真軟件Fluent 對探空溫度傳感器探頭進(jìn)行仿真分析及迭代計算。因模型中涉及到熱量傳遞，所以需打開能量方程；湍流選擇最常使用的k?epsilon模型；耦合速度與壓力的算法選用SIMPLE 算法；將空氣域氣流入口面和出口面分別設(shè)置為速度入口和壓力出口，墻面為流動墻面。對設(shè)定的各項條件進(jìn)行初始化及迭代計算。當(dāng)太陽輻射強度為1 000 W/m，海拔高度為16 km，氣流速度為3 m/s 時，溫度場分布如圖3 所示。

圖3 傳感器溫度場分布圖

由圖3 可知，左側(cè)入口面到右側(cè)出口面的溫度逐漸降低，計算得到熱電偶的升溫為0.124 K。

為了得到不同太陽輻射強度、海拔高度、氣流速度下的輻射誤差，利用CFD 方法對上述三個因素進(jìn)行仿真計算。海拔高度、太陽輻射強度、氣流速度變化范圍分別為10～32 km，200～1 000 W/m，1～6 m/s。仿真結(jié)果如圖4 所示。

圖4 氣流速度、海拔高度、太陽輻射強度與輻射誤差的關(guān)系

由圖4 可知：在輻射強度為200～1 000 W/m和氣流速度為1～6 m/s 的條件下，海拔高度和太陽輻射強度呈單調(diào)遞增關(guān)系，且輻射強度越大、海拔越高、氣流速度越小時，太陽輻射誤差越大；當(dāng)輻射強度為1 000 W/m、海拔高度為32 km、氣流速度為1 m/s 時，太陽輻射誤差最大為0.36 K；當(dāng)輻射強度為200 W/m、海拔高度為10 km、氣流速度為6 m/s 時，太陽輻射誤差最小為0.024 K。

2 基于SVM 誤差修正算法

本文采用SVM 對太陽輻射誤差數(shù)據(jù)進(jìn)行擬合，得到對應(yīng)的輻射誤差修正曲線。SVM 不過分依賴樣本集的數(shù)量，擁有較好的泛化性能。支持向量回歸SVR（Support Vactor Regression）是SVM 擬合曲線的重要方法，其基本思想是找到一個核心平面，讓所有的樣本點到該平面的距離最近，從而利用該平面進(jìn)行預(yù)測。

2.1 支持向量機回歸

在SVR 中，假設(shè)給定一個訓(xùn)練集合{(x,y) ,=1,2, }，其中x是訓(xùn)練集的輸入，y是訓(xùn)練集的輸出。首先將使用非線性映射映射()到一個更高維的空間，然后在此空間中構(gòu)建一個線性回歸模型，公式為：

式中：(,)是估計函數(shù)；和是模型的權(quán)向量和閾值；()是輸入向量。SVR 定義了一種不敏感損失函數(shù)，用來衡量回歸估計函數(shù)的正確性，其公式為：

SVM 回歸的目標(biāo)就是使用讓和(,)之間的差值越小，即如果和(,)之間的差值小于，那么損失值等于0。此時可將問題轉(zhuǎn)化為求解優(yōu)化目標(biāo)函數(shù)最小化范圍，其優(yōu)化公式為：

約束條件為：

式中：是常數(shù)；ξ,為松弛變量，表示訓(xùn)練集的誤差。利用拉格朗日函數(shù)解決約束最優(yōu)化問題，將其轉(zhuǎn)變?yōu)閷ε紗栴}，通過求解對偶最優(yōu)化問題得到線性回歸擬合函數(shù)為：

支持向量機回歸預(yù)測的步驟如下：

1）選取訓(xùn)練樣本的輸入、輸出。將太陽輻射強度、海拔高度、氣流速度作為SVM 的樣本輸入，太陽輻射誤差作為樣本輸出。

2）數(shù)據(jù)歸一化。將樣本數(shù)據(jù)映射到[-1，1]或[0，1]之間，消除誤差較大的數(shù)據(jù)，加快數(shù)據(jù)訓(xùn)練速度。

3）選用核函數(shù)類型。通常選用徑向基核函數(shù)，公式為：

4）選擇最佳參數(shù)和。用得到的值對訓(xùn)練集進(jìn)行訓(xùn)練獲得SVM 模型。

5）數(shù)據(jù)測試。用訓(xùn)練得到的模型進(jìn)行數(shù)據(jù)測試，并與實驗值進(jìn)行比較。

2.2 結(jié)果與分析

本文選用的程序庫是LibSVM。LibSVM 是2001年開發(fā)出來的一套SVM 的庫，運行速度快、運用靈活，便于對數(shù)據(jù)做分類或回歸。通過訓(xùn)練得到：=5.657；=0.088 4；=0.547 8；=-0.151；標(biāo)準(zhǔn)支持向量個數(shù)為163。將SVM 訓(xùn)練的結(jié)果與仿真結(jié)果進(jìn)行比較，如圖5所示。

圖5 仿真值與訓(xùn)練值對比

由圖5 可以看出仿真值與訓(xùn)練值基本接近。為了便于結(jié)果分析，將仿真值與訓(xùn)練值相減得出絕對誤差，如圖6 所示，誤差值介于[-0.003，0.003]。

圖6 仿真值與訓(xùn)練值的絕對誤差

3 實驗測試與分析

3.1 溫度傳感器電路設(shè)計

本文的探空溫度傳感器硬件系統(tǒng)主要包括高精度溫度采集模塊、主控制器、無線通信模塊、電源模塊四個部分。其中主控制器選用STM32F103C8T6 芯片；測溫電 路 選 用24 位∑Δ 模 擬 數(shù) 字 轉(zhuǎn) 換 器（Analog to Digital Converter）AD7794 芯片；無線通信采用ESP8266 WiFi模塊。溫度測量電路框圖如圖7 所示。

圖7 溫度測量電路框圖

3.2 實驗平臺搭建

為了驗證傳感器的測量精度，本文使用太陽模擬器和低氣壓風(fēng)洞搭建一個模擬高空氣壓和太陽光照狀態(tài)的實驗平臺，實驗裝置如圖8 所示。將傳感器放置在低氣壓風(fēng)洞石英玻璃腔內(nèi)，通過調(diào)節(jié)空氣進(jìn)氣口的氣流量來改變空氣流速，模擬探空溫度傳感器勻速升空時的氣流速度；真空腔體內(nèi)部氣壓可以提供不同海拔高度的大氣密度；使用CMP?10 型總輻射傳感器測量輻射強度，可以提供1 000 W/m輻射強度。

圖8 實驗裝置

3.3 實驗結(jié)果分析

將傳感器放置于實驗裝置中，對海拔高度、氣流速度、太陽輻射強度與太陽輻射誤差的關(guān)系進(jìn)行研究。選擇太陽輻射強度為600 W/m時的實驗數(shù)據(jù)進(jìn)行詳細(xì)說明，實驗數(shù)據(jù)如表1 所示。

表1 不同海拔高度、氣流速度下的輻射誤差 K

由表1 可知，氣流速度越小、海拔高度越高時，太陽輻射誤差越大；當(dāng)海拔高度為32 km，氣流速度為1 m/s時，輻射誤差最大為0.353 K。選取氣流速度為2 m/s時，將帶有壓電振子的傳感器與傳統(tǒng)熱電偶傳感器實驗測量得到的誤差數(shù)據(jù)進(jìn)行對比，如表2 所示。將實驗數(shù)據(jù)和算法預(yù)測數(shù)據(jù)進(jìn)行對比，如圖9 所示。

表2 使用壓電振子傳感器與傳統(tǒng)傳感器輻射誤差對比 K

圖9 實驗數(shù)據(jù)和預(yù)測數(shù)據(jù)對比

由表2 可知，當(dāng)太陽輻射強度為600 W/m、氣流速度為2 m/s 時，不使用壓電振子的傳感器輻射誤差最小為0.315 K，最大為0.673 K；使用壓電振子的傳感器輻射誤差最小僅0.120 K，最大為0.331 K。

由圖9 可知，實驗值與預(yù)測值的誤差較小。經(jīng)計算得到該傳感器平均測量誤差值為0.050 K，均方根誤差為0.051 K。

4 結(jié)論

為了提高氣流速度、減少太陽輻射誤差對熱電偶的影響，本文利用壓電陶瓷的逆壓電效應(yīng)設(shè)計一種基于熱電偶的探空溫度傳感器，通過CFD 仿真對不同海拔高度、氣流速度、太陽輻射強度進(jìn)行太陽輻射誤差計算，使用SVM 算法對誤差值進(jìn)行修正，并搭建一個模擬高空環(huán)境的實驗平臺對傳感器的測量精度進(jìn)行驗證，得出以下結(jié)論：

1）傳感器的太陽輻射誤差與海拔高度呈單調(diào)遞增關(guān)系，與氣流速度呈遞減關(guān)系，且輻射強度越大、海拔越高、氣流速度越小時，太陽輻射誤差越大。

2）通過SVM 算法對仿真值進(jìn)行訓(xùn)練得到的訓(xùn)練值與仿真值之間的誤差介于[-0.003，0.003]，擬合精度較高。

3）熱電偶傳感器的測量誤差平均值為0.050 K，均方根誤差為0.051 K，表明傳感器的測量精度較高，壓電陶瓷振子能有效降低太陽輻射誤差。

本文設(shè)計的溫度傳感器對比傳統(tǒng)溫度傳感器能有效降低太陽輻射誤差，但仍未考慮太陽高度角等因素的影響，后續(xù)會繼續(xù)研究并完善。