丁 楓，劉清惓，楊 杰 ，陳高穎， 袁 宇，付川琪

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

0 引言

當太陽光線照射到溫度傳感器表面時，通過輻射吸收及熱傳導的方式造成的溫升將疊加在傳感器所測的環(huán)境溫度上，由此造成的測量誤差被稱為太陽輻射誤差[1]。目前常用的傳感器存在輻射誤差，E. Erell等[2]通過理論分析和實驗證明了太陽輻射誤差與太陽輻射強度成正相關，而與氣流速度負相關。因此，要利用防輻射溫度傳感器來減弱太陽輻射帶來的影響。根據(jù)不同的熱交換方式，防輻射溫度傳感器分為自然通風與強制通風。R. Nakamura等[3]通過實驗研究發(fā)現(xiàn)，不同型號的防輻射溫度傳感器在不同環(huán)境參數(shù)下的輻射誤差呈現(xiàn)差異，例如，USCRN(U. S. climate reference network)與MMTS(maximum and minimum temperature system)相比，夜間差異在0.2～0.4 K，白天的差異則多達1.0 K[4]。日本學者研究發(fā)現(xiàn)百葉箱的誤差和風速、太陽方向角等環(huán)境因素相關，引起的誤差約為1 K，還有學者發(fā)現(xiàn)誤差更大的實驗結(jié)果。因此，目前常用溫度傳感器難以滿足氣象觀測日益提高的精度要求。

為提高地表氣溫觀測精度，本文提出一種低輻射誤差的陣列式球形溫度傳感器，且無需安裝于百葉箱或防輻射罩內(nèi)，利用計算空氣動力學(CFD)方法計算出傳感器相互之間的輻射誤差比值，并采用極限學習機(ELM)算法獲得該比值與各環(huán)境因素的關系，為了驗證陣列式球形溫度傳感器的觀測精度，開展基于標準氣溫測量實驗的太陽模擬平臺數(shù)據(jù)測量實驗。

1 陣列式球形溫度傳感器設計

1.1 傳感器工作原理

陣列式球型傳感器包含3個具有不同反射率涂層的溫度傳感器探頭。在同一環(huán)境下，3個探頭的輻射誤差不同，均與太陽輻射強度、環(huán)境風速、下墊面反射率等環(huán)境因素有關。因此只需量化不同探頭之間的輻射誤差比值與各環(huán)境因素的關系，即可推導出真實環(huán)境溫度，傳感器結(jié)構(gòu)如圖1所示。

圖1 原理示意圖

實驗中，傳感器探頭A表面涂覆一層反射率約0.2的黑色涂層。探頭B表面無涂層，銅球反射率約0.5。探頭C表面涂覆一層反射率約0.2的白色涂層。無輻射環(huán)境下，3個探頭的溫度測量值均為Tx。當施加一定強度太陽輻射后，溫度傳感器A的溫度為Ta，輻射誤差為ΔTa；溫度傳感器B的溫度為Tb，輻射誤差為ΔTb；溫度傳感器C的溫度為Tc，輻射誤差為ΔTc。此時傳感器探頭之間的輻射誤差比值：

(1)

(2)

η(P,ν)，λ(P,ν)為關于P，ν的特定函數(shù)，P為太陽輻射強度，ν為風速。由此可推算出等式：

(3)

(4)

由上述計算等式變化可推算出大氣溫度真實值Tx：

(5)

(6)

(7)

將求得Tx1與Tx2的平均值作為真實環(huán)境溫度最終值Tx，其中式(5)與式(6)為輻射誤差比值換算公式。

1.2 傳感器及測量電路設計

本文設計了一種通過輻射誤差比值反推環(huán)境溫度真實值的陣列式球形溫度傳感器。利用空心銅球作為外殼，采用Pt100鉑電阻作為測溫探頭，并通過硅酮導熱膠使銅球與鉑電阻貼合。根據(jù)黑色、白色與原色3種涂層使球形溫度傳感器具有不同的太陽輻射吸收率，并通過支撐片每2.5 cm固定一種球型溫度傳感器探頭。陣列式球形溫度傳感器結(jié)構(gòu)圖如圖2所示。

圖2 陣列式球形溫度傳感器結(jié)構(gòu)圖

高精度測溫電路由STM32F103C8T6處理器、陣列式傳感器數(shù)據(jù)采集系統(tǒng)和串口傳輸系統(tǒng)組成，將陣列式球形溫度傳感器采集的數(shù)據(jù)經(jīng)過串口發(fā)送到上位機顯示與分析，具體設計框圖如圖3所示。

圖3 溫度測量電路系統(tǒng)設計框圖

本設計溫度采集系統(tǒng)中的傳感器探頭采用4線制接法[5]，運用比率法消除引線誤差，提高測量精度[6]。利用Isotech的鎵熔點室和水三相點室進行鉑電阻標定實驗，并利用Fluke1595A測量電橋?qū)φ鎸崪囟冗M行記錄，該測溫電橋儀器具有較好的穩(wěn)定性和準確性，電壓噪聲達到3 nV級別。

2 計算流體動力學仿真分析

2.1 傳感器物理模型建立

本文使用三維實體建模軟件Pro/E建立陣列式球型溫度傳感器的物理模型?？招你~球的直徑、內(nèi)壁直徑、開孔直徑分別為6 mm、5.5 mm、3 mm，將3個不同吸收率的傳感器小球，間隔6 mm放置，確保區(qū)域環(huán)境的一致性，傳感器基本材料參數(shù)見表1。

園林花卉養(yǎng)護與管理，直接影響園林工程整體質(zhì)量與綜合效益，故此管理人員一定要給予其一定重視，并編制相關的質(zhì)量管理標準。結(jié)合園林花卉工程具體養(yǎng)護管理項目，依照相關標準進行質(zhì)量檢驗，將基礎性養(yǎng)護設為始發(fā)點，實現(xiàn)對園林花卉養(yǎng)護管理工作質(zhì)量的有效管控。

表1 傳感器基本材料參數(shù)

然后，使用軟件ICEM對傳感器幾何模型進行網(wǎng)格劃分，利用長方體空氣域包裹3個溫度傳感器的方式，建立整個仿真模型。為滿足仿真的收斂性和準確性的要求，將長方形空氣域設置為120 mm×120 mm×360 mm，并分別將正方形面設置為速度入口與壓力出口，傳感器及空氣域網(wǎng)格。經(jīng)系統(tǒng)分析，該模型的網(wǎng)格質(zhì)量大于0.35，網(wǎng)格質(zhì)量基本滿足FLUENT仿真要求[7]。

2.2 流固耦合數(shù)值計算

本文采用CFD方法對傳感器進行仿真分析，首先利用熱流密度表示太陽輻射強度，其次采用能量方程計算輻射傳熱、對流換熱，并選用k-e模型和壓力與速度進行SIMPLE算法。針對動量、能量的求解，可以先選定一階迎風格式，再采用二階迎風方式。對于設置模型邊界條件，選擇空氣域的進風口為速度入口，出口為壓力出口，空氣域四周設置為滑移壁面，并將3個球型傳感器設置不同的表面吸收率，白色涂層設置為20%，原色涂層為50%，黑色涂層為80%。為模擬真實實驗環(huán)境的需求，本文設置不同的太陽輻射強度和風速進行仿真研究。

當太陽輻射強度為1 000 W/m2，風速為3 m/s時，陣列式球形傳感器的溫度場分布如圖4所示。上、中、下層分別為黑色涂層、原色涂層和白色涂層，在輻射與熱對流的影響下，誤差升溫分別為：1.67 K、1.37 K、1.04 K。

圖4 傳感器溫度場分布圖

2.3 仿真結(jié)果分析與研究

本文通過計算流體動力學仿真的方法，對不同風速和太陽輻射強度下的陣列式球形溫度傳感器進行輻射誤差比值分析。太陽輻射強度為100～1 050 W/m2，步長為50 W/m2；風速變化范圍為1～5 m/s，步長為1 m/s，黑色涂層傳感器探頭與原色涂層傳感器探頭溫差比值記為ΔTa/ΔTb；原色涂層傳感器探頭與白色涂層傳感器探頭溫差比值記為ΔTb/ΔTc，仿真數(shù)據(jù)結(jié)果如圖5所示。

(a)黑色與原色涂層仿真數(shù)據(jù)結(jié)果

(b)原色與白色涂層仿真數(shù)據(jù)結(jié)果圖5 仿真數(shù)據(jù)結(jié)果表

由上述可知，黑色涂層傳感器探頭與原色涂層傳感器探頭溫差比值ΔTa/ΔTb在2.8上下區(qū)間浮動時，數(shù)值較為穩(wěn)定，但部分區(qū)間有較為明顯的波動，在風速為2 m/s，太陽輻射強度為500 W/m2時，數(shù)值波動達到最大，波動值為0.059 5；原色涂層傳感器探頭與白色涂層傳感器探頭溫差比值ΔTb/ΔTc在1.6上下區(qū)間浮動，整體趨勢與ΔTa/ΔTb相同。在風速為5 m/s，太陽輻射強度為100 W/m2時，數(shù)值波動較大，波動值達到0.023 5。因此對于陣列式球形溫度傳感器中的輻射溫差比值，并不能單一的認為是一個常數(shù)，而是與環(huán)境因素相關的多元非線性函數(shù)。

2.4 輻射誤差比值計算

為解決溫度傳感器溫差比值波動問題以及CFD仿真僅能計算分析離散模型的問題，本文采用極限學習機(ELM)方法[8]，使用前饋型神經(jīng)網(wǎng)絡，獲得不同太陽輻射強度與風速下的溫差比值誤差訂正方程。利用CFD仿真得到的數(shù)據(jù)作為極限學習機神經(jīng)網(wǎng)絡的訓練組，將激活函數(shù)設定為sigmoid函數(shù)，通過計算輸出矩陣的廣義逆矩陣，通過最小二乘法取得連接權(quán)值，最終得出數(shù)學模型。本文將風速與太陽輻射強度作為輸入變量，陣列式球形溫度傳感器各測溫探頭之間的溫差比ΔTa/ΔTb、ΔTb/ΔTc作為2個輸出變量，分析得出輻射誤差修正方程的數(shù)學模型。采用ELM算法擬合生成的CFD仿真值與修正模型計算值的對比圖和誤差圖如圖6和圖7所示。

(a)ΔTa/ΔTb仿真值與計算值對比圖

(b)ΔTa/ΔTb仿真值與計算值誤差圖圖6 黑色與原色涂層修正模型結(jié)果

根據(jù)誤差對比圖可知，修正比值與仿真比值之間的均方根誤差為0.024 5 K，故該算法的回歸擬合精度良好。輸入環(huán)境因素中的太陽輻射強度數(shù)據(jù)與風速數(shù)值，可計算出球形溫度傳感器的溫差比值，并通過換算公式得到環(huán)境真實溫度值，進一步提高了陣列式球形溫度傳感器測量的準確性。

(a)ΔTb/ΔTc仿真值與計算值對比圖

(b)ΔTb/ΔTc仿真值與計算值誤差圖圖7 原色與白色涂層修正模型結(jié)果

3 實驗測試及分析

3.1 實驗平臺搭建

實驗過程中，利用低氣壓風洞與太陽模擬器搭建本實驗觀測平臺，將陣列式球形傳感器放置在低氣壓風洞中心處，保證區(qū)域環(huán)境的一致性與穩(wěn)定性。該低氣壓風洞包括2臺雙極旋片式氣泵、1根石英玻璃管和2個氣壓檢測儀，利用Fluke1595A測溫電橋測得的標準數(shù)據(jù)作為對照組，該測溫電橋的測量準確度在典型的測溫學比例(0.25～4.0)范圍內(nèi)，太陽模擬器采用全光譜射線模擬系統(tǒng)，系統(tǒng)內(nèi)部的球形短弧氙燈具有匹配度極佳的光譜，且輸出功率大、亮度高[9]，最大可模擬1 000 W/m2的輻射強度。使用CMP-10的測量結(jié)果作為輻射強度標準，確保實驗的真實性與可靠性，實驗平臺現(xiàn)場圖片如圖8所示。

圖8 太陽模擬實驗平臺

3.2 實驗驗證及數(shù)據(jù)分析

由表2可知，陣列式球型溫度傳感器中黑色涂層與原色涂層的輻射誤差比值ΔTa/ΔTb修正值與實驗差值的平均值為0.042；原涂層與白色涂層的輻射誤差比值ΔTb/ΔTc修正值與實驗差值的平均值為0.053。通過輻射誤差比值換算公式計算出當前環(huán)境真實溫度，與測溫電橋測得標準環(huán)境溫度值進行對比，對照表如表3所示，該傳感器溫度誤差的平均值為0.145 ℃，均方根誤差為0.16 ℃。

表2 輻射誤差比值修正值與實驗值對照表

表3 計算溫度與測量溫度對照表

4 結(jié)論

為降低太陽輻射對溫度傳感器的輻射升溫影響，本文設計了一種陣列式球形溫度傳感器，通過輻射誤差比值，計算得出環(huán)境真實溫度。采用計算流體動力學方法對該傳感器在不同太陽輻射強度與風速條件下進行仿真計算，得出輻射誤差比值，然后采用極限學習機方法對數(shù)據(jù)進行擬合處理，最后搭建太陽模擬實驗平臺對仿真比值進行實驗驗證分析，得出以下結(jié)論：

(1)陣列式球形溫度傳感器的輻射誤差比值是一條波動幅度較小的曲線，風速與太陽輻射強度對輻射誤差比值影響較小，但無法忽視環(huán)境因素對輻射誤差比值的影響。

(2)陣列式球形溫度傳感器的輻射誤差修正比值與仿真值比值之間的均方根誤差為0.024 5 K，說明該算法擬合效果良好，提高了陣列式球型溫度傳感器的測量精度。

(3)陣列式球形溫度傳感器的溫度測量誤差均值為0.145 K，均方根誤差為0.16 K，具有良好的測量精度，并符合大氣科學研究對地表低溫觀測的精度要求。

本文設計的陣列式球形溫度傳感器具有較好的測量精度，但未考慮海拔高度和散射輻射等因素的影響，在后續(xù)研究中會進一步完善傳感器設計和CFD模型，使其精度量級達到更高。