|
|
Parsivel雨滴譜儀及其在氣象領域的應用 |
| 來源：科技論文 發布日期：2009-03-17 16:55:07 |
|
濮江平 趙國強 蔡定軍 姚展予 袁冬梅 呂 梅 |
|
(解放軍理工大學氣象學院,南京211101;河南省氣象局,鄭州450003;
江西省氣象局,南昌330046;中國氣象科學研究院,北京100081) |
摘 要
降水粒子特性是大氣運動和云內微物理過程的綜合結果,在云降水物理及人工影響天氣領域有著重要的意義。傳統的測量方法不適合對大量數據分析尋找規律,德國OTT公司的Parsivel“激光降水粒子譜測量系統可以較好解決自動測量難題。該儀器是以激光測量為基礎的粒子測量傳感器,采用平行激光束和光電管陣列結合,當有降水粒子穿越采樣空間時,自動記錄遮擋物的寬度,通過穿越時間計算降水粒子的尺度和速度,根據各種參數的綜合信息對降水粒子進行分類,并能夠以數字形式顯示瞬時降水強度、降水粒子總數、累積降水量、降水時的能見度和雷達反射因子,以圖形方式顯示降水粒子尺度譜、速度譜、降水粒子分類且自動生成天氣現象代碼,實現天氣現象的自動識別。
激光降水粒子譜儀主要用于氣象水文觀測。在雷達氣象學領域可用于Z /R關系的擬合修正,比傳統的用雨量筒觀測數據擬合效果好得多;由降水粒子譜儀測量雨滴的降落速度,可以對天氣雷達垂直向上測量的粒子徑向速度譜進行校正。
人工影響天氣的效果檢驗一直是一個難題,自然降水粒子譜分布形式與人工催化以后的降水粒子譜型理論上應當具有較大的差別,人工增雨作業降水滴譜變化物理響應和降水強度時間變化響應都有明顯的區別。如果能夠實時檢測到這些差別,就能夠充分說明人工催化的有效性。未來如果能夠進行聯網觀測記錄區域性降水、降雹,就有充分證據表明人影作業的有效性,在定量化作業效果評估以及災害損失評估等方面應用潛力巨大。
利用該儀器已經對一年的自然降水過程進行了連續觀測,并將所獲得的降水粒子譜、雨滴濃度值隨時間變化狀況與衛星反演的云頂有效粒子半徑時間變化趨勢進行了對比,發現有較好的一致性。
關鍵詞: 降水強度;降水量;雷達反射率;降水粒子測量
中圖分類號: TH765. 6 + 6
文獻標識碼: A
文章編號: 1673 - 7148 (2007) 02 - 0003 - 06
引 言
降水粒子特性是云內動力過程和微物理過程的綜合結果,在云降水物理研究,特別是人工影響天氣領域有著重要的意義。降水量測量歷史悠久,雨滴譜研究也由來已久。反映降水粒子譜的特性參數, 主要有雨滴的直徑大小、組成分布、降落速度和動能。Marshall and Palmer ( 1948)[ 1 ]提出了降水粒子譜分布, Gunn and Kinzer (1949)[ 2 ]得到了降水粒子速度譜分布。雨滴直徑的大小,決定了降落雨滴的質量和速度的大小,從而決定了雨滴所具有動能的大小,所以測量雨滴直徑是首要的。目前常見的測定方法有方格法、快速攝影法、浸入法、面粉球法、濾紙色斑法等。傳統的測量方法工作量大,勞動強度高,無法自動完成測量分類工作,不適合對大量數據分析尋找規律。Knollenberg(1970)[ 3 ]研制了光學陣列探頭(OAP) ,用于測量單個雨滴粒子尺度,但Ill2 ingsworth and Stevens (1987)[ 4 ]認為,用OAP測量速度不夠理想, Baumgardner ( 1993)[ 5 ]和Jameson and Kostinski的一系列文章研究了降水粒子的特性 ( 2000, 2001, 2002)[ 6-8 ]。本文介紹一個測量降水粒子尺度和粒子下落速度的觀測系統???德國OTT 公司的Parsivel激光降水粒子譜測量系統(Laser2 based op tical Disdrometer for simultaneous measure2 ment of Particle Size and Velocity of all liquid and sol2 id p recip itation) ,可以較好解決自動測量難題。
Parsivel降水粒子譜儀是以激光為基礎的新一代光學粒子測量器及氣象傳感器,可同時測量降水中所有液體和固體粒子的尺度和速度,采用的是直接測量體制,并對降水粒子進行分類。它可滿足氣象、水文的傳感探測器需求,并達到了世界氣象組織和美國國家氣象局相關技術規定的要求。在所有天氣狀況下的降水過程, Parsivel都提供了完整的圖像記錄,并有準確的降水類型和強度報告,而且在惡劣戶外天氣中仍有很好的表現。它的一體化加熱裝置可以把冰凍天氣對儀器表面的負面影響降到最低,并可在任何天氣狀況下進行作業。
1 Parsivel激光降水粒子譜儀介紹
1. 1 主要技術指標
Parsivel降水粒子譜儀采用平行激光束為采樣空間,光電管陣列為接受傳感器,當有降水粒子(無論固態還是液態甚至沙塵)穿越采樣空間時,自動記錄遮擋物的寬度和穿越時間,從而計算降水粒子的尺度和速度。儀器的主要技術指標見表1,室外傳感器部分見圖1。
|
表1 Parsivel主要技術指標 |
|
|
|
|
|
圖1 Parsivel室外傳感器外形 |
1. 2 儀器功能
Parsivel降水粒子譜儀是用以激光測量為基礎的光學傳感器,通過測量遮擋物的寬度和通過時間來計算降水粒子的尺度和速度。傳感器的發射部分產生一水平平直的激光束帶,接收傳感器把它轉化成電子信號。當降水粒子通過光束帶下落時信號將改變,儀器自動記錄被遮擋的觀點傳感器個數和時間。由Parsivel快速運算記錄每一個降水粒子的尺度和速度,并通過一臺裝有其配套軟件ASDO的電腦對降水粒子譜進行分析統計,將結果以圖表數據形式進行顯示,進而得出累積參數。同時ASDO還能夠對室外傳感器進行實時監控,包括激光光束信號強度、環境溫度、電源狀況、加熱電流甚至鏡頭污染需要維護等信息。ASDO有兩種顯示模式:在線模式用于采集室外傳感器數據,離線模式用于歷史數據回放。兩者都能夠用二維或三維圖形表示粒子尺度譜與速度譜分布。在在線模式下能夠將傳感器數據記錄并傳到電腦上,數據采集時間間隔可以人為設定,最短采樣周期為10 s,最長為2 h,然后保存在數據庫中。在離線模式下可以回放歷史數據。
ASDO計算統計顯示軟件的界面可同時顯示降水粒子尺度譜、速度譜;同時可顯示降水粒子分類, 可分為無降水、毛毛雨、雨(小雨、雨、大雨) 、冰粒子、冰雹(小冰雹、冰雹、大冰雹) 、雪(小雪、雪、大雪) 、混合型降水等類型,并且自動生成天氣現象代碼、氣候代碼和航空氣象代碼,實現自動識別天氣現象;另外還能夠以數字形式顯示瞬時降水強度、粒子總數以及累積降水量,可代替雨量筒和翻斗式雨量計;通過粒子譜可以計算降水時的能見度和雷達反射因子。
Parsivel降水粒子譜儀主要用于氣象水文觀測。其瞬時降雨強度測量能力有助于早期洪水預測預警;其粒子譜測量能力在人工影響天氣效果檢驗方面有著重要的應用潛力;其粒子濃度測量能力及計算雷達反射因子能力在氣象雷達定標和雷達Z /R 關系計算方面有著諸多應用。這些應用還有待于進一步開發。可以預計,降水粒子譜儀在常規氣象觀測、航空氣象保障、公路氣象保障、洪水預測預警、天氣雷達等領域將得到廣泛應用。
1. 3 工作原理
Parsivel降水粒子譜儀的傳感器由發射機、接收機和控制、運算、存儲等電路部分組成。所產生的高穩定度水平平行激光束帶和穩定的接受光電陣列是傳感器的核心部件。高速運算和數據存儲是測量的基礎。高速響應光學接收光電陣列傳感器把遮蔽物圖像轉化成電子信號,通過即時測量每一個遮擋物的寬度和通過時間來計算降水粒子的尺度和速度。當降水粒子通過光束帶下落時信號將改變,儀器自動記錄被遮擋的觀點傳感器個數和時間。由Parsivel快速運算記錄每一個降水粒子的尺度和速度 (見圖2) 。
粒子尺度的測量:降水粒子下降過程中經過探測儀器的采樣空間時,粒子的直徑a由被遮擋的接收二極管d個數確定。
粒子速度的測量:粒子的下降速度V 由接收二極管d被遮擋的時間t和粒子直徑a計算確定,即V = f ( d, t) 。假設粒子尺度均勻非變形, 則有Vt = f ( a) 。
由于計算速度足夠快,系統能夠響應所有經過的遮蔽物,從而可以統計所有降水粒子顆粒的尺度與速度。但是如果有兩個甚至多個降水粒子同時到達采樣面時,就會產生重疊誤差,無論是對尺度測量還是對速度的計算,都會產生很大的影響。
2 降水粒子譜分析
2. 1 譜的定義
在云降水物理學和大氣氣溶膠研究領域,研究對象具有明顯的多分散性。為了全面描述分散系的物理特性,應對質粒尺度的分布特征進行定量描述。通常采用質粒濃度隨尺度的分布來定量描述顆粒物的物理特性,這種尺度分布常稱為譜或譜分布。譜分布有兩種形式,即離散分布和連續分布。實際測量數據總是離散的,為了數學處理和表述的方便,常轉化成連續性函數形式。
2. 2 數據分析
Parsivel降水粒子譜儀測量的數據共有32個尺度測量通道和32個速度測量通道。其中粒子尺度測量32個通道對應的數據范圍為0. 2～25 mm,粒子速度測量通道數據范圍對應著0. 2～20 m / s (見表2) 。每一次采樣間隔內的粒子譜測量數據都有 32 ×32 = 1 024個。以2006年3月份以來的幾次降水資料為例進行觀測分析,發現實際的數據中并不能涵蓋所有的1024個通道,常見的只有0. 2～8 mm、 0. 2～14. 4 m / s。
2. 3 Parsivel的觀測個例分析
2006年3月11日降雨過程是一次典型的鋒面過境持續的降水過程,整個降水過程從3 月11 日 12: 30 一直持續到3 月12 日01: 30,總降水量為 7. 67 mm,在10 s間隔的采樣時間段內,最大降水強度為7. 498 mm /h,雨滴濃度最大值0. 1630個?cm- 2?s- 1 (見圖3) 。
從圖3a中分析可見,整個降水過程分為5個峰值過程。其中,中間過程降水強度和雨滴濃度都比較大,而起始和結束的降水都是強度小濃度較大,后面幾個過程也有類似的現象,尤其是降水結束前雨滴濃度都比較大,但降水量卻比較小。圖3a中的幾個瞬時降水強度高值區對應著圖3b中累積降水量躍升臺階。這種時間高分辨觀測數據用傳統的雨量筒是無法獲得的。圖中的彩色塊狀圖就是降水粒子譜的分布,其中的擬合彎曲線為Hobbs在美國Cascade計劃中從大量觀測數據中分析獲得的降水擬合曲線。圖3c是整個降水過程雨滴譜隨時間的變化狀況。
2. 4 Parsivel的觀測數據的變換
ASDO保存數據格式以32 ×32順序排列1024 個數據形式,如果需要按照常規的尺度與濃度柱狀圖關系排列,需要對原始數據進行格式變換。
先從Parsivel獲得的二維數據,每個時間對應的都有1024個數據,橫坐標是尺度(mm) ,縱坐標是速度(m / s) ,然后把各個尺度分檔對應的一列速度數據進行累加,這樣縱坐標就成了質粒數濃度,與橫坐標尺度配合就可以初步畫出對應的譜分布,利用Excel可做出柱狀譜分布圖(圖4) 。
3 Parsivel的應用及其優缺點
Parsivel降水粒子譜儀能夠實現無人值守,自動記錄,維護保養方便,并可在任何天氣條件下作業; 降水粒子的尺度和速度二維、三維分布的全面分析, 顯示界面直觀易懂;實現非接觸式測量,對原始場無干擾;工藝上實現防雨滴飛濺干擾測量精度;可靠地識別WMO天氣現象,自動生成電碼;自控電加熱系統能夠適應凍雨環境工作。因此,它在常規氣象觀測和水文觀測洪水預警方面被廣泛應用,在天氣雷達探測和人工影響天氣領域有著潛在的應用價值。
3. 1 在天氣雷達探測和雷達氣象學領域潛在的應用
目前氣象雷達學領域常用的Z /R關系都是用雨量筒觀測數據進行訂正(數據是每小時的降水量),雷達反射因子Z定義為單位體積內所有粒子直徑D的6次方之和,即
雨強R 定義為單位時間內落到單位面積上的降水質量,所以雨強與落到單位面積上總的雨滴質量 (與雨滴直徑D 的3次方成正比)和雨滴下落速度有關。可見,反射因子Z和雨強R 都和雨滴直徑有關,則Z /R 關系可用
表示。雷達探測時間上的高分辨率與雨量筒的低反應效率之間很難比較,而利用Parsivel測得的高時間分辨率降水粒子譜數據和瞬時降水強度數據可以擬合較為準確的雷達Z /R 關系。
Parsivel降水粒子譜儀能夠實際測量雨滴的降落速度,可以對天氣雷達垂直向上測量的徑向速度譜進行校正。雷達的后向散射功率與降水粒子譜直接相關,而與降水量的相關性不是唯一的,同樣波長的雷達探測小球形粒子時,粒子的半徑愈大,后向散射截面愈大。降水量不變而降水粒子譜變化時,雷達回波功率肯定發生變化。因此雷達探測降水量和降水強度急需降水譜資料,這在雷達校驗定標方面作用潛力巨大。
3. 2 人工影響天氣領域潛在的應用
人工影響天氣的效果檢驗一直是一個難題,目前主要是統計檢驗和物理檢驗。統計檢驗需要大量試驗樣本,而且隨機化試驗在我國難以實現;物理檢驗目前主要是雷達觀測演變狀況和飛機直接觀測云內微物理參數,這些方法不是直接測量定量化困難, 就是成本昂貴。自然降水粒子譜分布形式與人工催化以后的降水粒子譜形式理論上應當具有較大的差別,人工增雨作業降水滴譜變化物理響應和降水強度時間變化響應都有明顯的區別。如果能夠實時檢測到這些差別,就能夠充分說明人工催化的有效性。同時無論是自然降水還是人工催化降水的強度、時間變化能夠被精確記錄,也有益于分析作業效果。
在人工防雹作業效果檢驗方面, Parsivel降水粒子譜儀能夠記錄降雹的最大尺度、持續時間、雹譜形式的變化及降雹能量計算,在人工抑制冰雹作業中對降雹譜變化、降水與降雹的區分和時間分布及能量計算,無論是效果檢驗還是雹災損失估算都具有很高的利用價值。未來如果能夠進行聯網觀測記錄區域性降水、降雹,就有充分證據表明人影作業的有效性及其效果評估定量化,應用潛力巨大。
3. 3 Parsivel降水粒子譜儀的缺點
任何觀測儀器都有其弱點, Parsivel降水粒子譜儀也同樣有其弱點,主要是采樣空間的確定和采樣周期的選擇。采樣空間受到工藝技術的限制。采樣寬度受到接受光電陣列的限制(目前是30 mm寬) , 采樣長度受到激光強度(目前是180 mm長)和光電陣列的靈敏度限制。采樣面積太小時,由于小雨雨滴密度很低,采樣缺乏代表性;采樣面積大時,對于儀器的穩定性有影響。此外,在大雨時候,雨滴重疊時的測量難題無法解決,樣本數太大超過計數器上限將飽和溢出。
采樣周期的選擇應當與降水類型相匹配,小雨采樣周期應當長一些,以保證有足夠的樣本數進行統計;暴雨時應當選擇較短的采樣周期,以免計數器飽和產生溢出計算錯誤。根據已經探測到的數個降水過程分析,強暴雨時選擇10 ～15 s;毛毛雨時選 120～180 s;正常待機時綜合選擇60 s為宜。
參考文獻
[ 1 ]Marshall J S, PalmerW M. The distribution of raindrops with size [ J ]. Meteorol, 1948, 5: 165 - 166.
[ 2 ]Gunn R, Kinzer G D. The terminal velocity of fall forwater drop lets in stagnant air[ J ]. Meteorol, 1949, 6: 243 - 248.
[ 3 ]Knollenberg R G. The op tical array: An alternative to scattering and extinction for airborne particle size determination [ J ]. Appi Met, 1970, 9: 86 - 103.
[ 4 ] Illingsworth A J, Stevens C J. An optical disdrometer for measurement of raind rop-size spectra in windy conditions [ J ]. Atmos Oceanic Tech, 1987, 4: 411 - 421.
[ 5 ]BaumgardnerD B B, Weaver K. A technique for themeasurement of cloud structure on centimeter scales[ J ]. AtmosOceanic Tech, 1993, 10: 557 - 565.
[ 6 ] Jameson A R, Kostinski S B. Fluctuation properties of precipitation. Part V I: Observations of hyperfine clustering and drop size distribution structures in three-dimensional rain [ J ]. Atmos Sci, 2000, 57: 373 - 388.
[ 7 ] Jameson A R, Kostinski S B. What is a raind rop size distribution [ J ]. Bull AmerMet Soc, 2001, 82: 1169 - 1177.
[ 8 ] Jameson A R, Kostinski S B. When is rain steady [ J ]. Appl Met, 2002, 41: 83 - 90. | |
|
|
|
|
|
