結合MODIS云產品的黑龍江省2015年人工增雨潛力研究

摘要[目的]研究2015年黑龍江省人工增雨潛力的時空分布。[方法]利用再分析資料計算降水轉化率，并通過結合MODIS衛星云產品，從人工增雨的微物理條件出發進一步研究了黑龍江省2015年春季和秋季的增雨潛力。[結果]在春季的松嫩平原大部分地區以及秋季的松嫩平原與三江平原大部分地區，降水量和降水轉化率均較小，存在一定的增雨潛力；春季，松嫩平原北部和三江平原中部存在適宜進行暖云催化增雨的區域；秋季，從松嫩平原北部至三江平原西側存在許多適宜暖云催化增雨的區域，同時在松嫩平原北部與三江平原西側的局部地區存在適宜冷云催化增雨的區域。[結論]該研究可為科學合理開展人工增雨規劃及作業提供依據。

關鍵詞人工增雨潛力；降水轉化率；MODIS；云微物理；黑龍江省

中圖分類號S423+.9文獻標識碼A文章編號0517-6611（2017）11-0171-05

Abstract[Objective] The research aimed to study the temporal and spatial distribution of artificial precipitation potentiality in Heilongjiang Province in 2015.[Method] The precipitation conversion rate was calculated by using the reanalysis data. The potential of rainfall in spring and autumn of Heilongjiang Province in 2015 was further studied by combining the MODIS satellite cloud product and the microphysical conditions of artificial rainfall.[Result] In most parts of Songnen Plain in spring and in Songnen Plain and most of the Sanjiang Plain in autumn， the precipitation and conversion rate of precipitation were small， and there was some potential for increasing rainfall.In spring， there were areas suitable for warm cloudy catalytic enhancement precipitation in the northern part of Songnen Plain and the middle part of the Sanjiang Plain.In autumn， there were many areas suitable for warm cloudy catalytic enhancement precipitation from the northern part of Songnen Plain to the western part of Sanjiang Plain. At the same time， there were areas suitable for cold cloudy catalytic enhancement precipitation in the northern part of Songnen Plain and the western part of Sanjiang Plain.[Conclusion] The study provides a basis for scientific and reasonable planning and operation of artificial rainfall.

Key wordsArtificial precipitation potentiality；Precipitation conversion rate；MODIS；Cloud microphysics；Heilongjiang Province

人工增雨是開發空中水資源、解決水資源緊缺與保障生產的重要手段之一，對人工增雨潛力的充分研究是科學合理開展人工增雨規劃及作業的必要前提。黑龍江省作為我國重要的商品糧基地，其水資源的重要性不言而喻，因此開展針對黑龍江省的人工增雨潛力研究是非常必要的。

針對云水資源開發和增雨潛力的研究有多種途徑，如齊彥斌等[1]對東北冷渦云系開展的飛機穿云觀測試驗、洪延超等[2]利用數值模式對層狀云增雨潛力開展的理論研究。而在國內針對區域性增雨潛力的研究中，更多的是采用基于觀測數據的統計變量來表征增雨潛力的時空分布，如以CWR-PEP評估法統計量[3-4]或降水轉化率[5-7]的大小來表示增雨潛力，針對黑龍江省也曾有過這類研究[8]。但這2種統計變量均存在一定的局限性。其中CWR-PEP法統計量通常主要基于理論估計的10%～15%增雨效率和地面降水計算得到，降水轉化率則通常以整層大氣可降水量和地面降水的比值計算得到，然而就目前切實可行的業務化增雨手段而言，人工增雨作業仍然是在適宜的條件下以人為改變云水的微物理特性及其變化過程為主要增雨手段[9]，事實上，在不同的云微物理條件下，既無法保證一定有10%～15%的增雨效率，也無法確保整層大氣可降水量一定可以轉化為地面降水。因此，在研究區域性的增雨潛力時對該地區的云微物理特性進行討論是非常有必要的。

近年來，衛星遙感技術的發展為云水資源時空分布及其微物理特性的研究提供了極大的便利[10-12]，宋松濤等[13]利用國際衛星云計劃（ISCCP）中衛星反演的云水路徑數據（大氣柱中單位面積液態或固態水含量的總和）研究了西北地區的云水資源時空分布，林丹[14]利用MODIS衛星云產品分析了西南地區的云水特征，其中MODIS衛星3級云產品中提供并區分了液相和冰相2種相態的云水路徑、粒子平均尺度以及粒子尺度標準差等可以反映云微物理特性的反演量。筆者在利用降水轉化率對增雨潛力進行初步討論的基礎上，利用上述MODIS衛星云產品所表征的云微物理特性對增雨條件做出進一步分析。

1資料與方法

選取歐洲中心ERA-Interim再分析數據集中0.5°×05°分辨率的2015年逐月的月累計降水和月均氣柱內總含水量來表征黑龍江省2015年的降水及水汽條件情況，并選取MODIS衛星3級產品中的月平均云產品（1°×1° 格點產品），利用其中的液相和冰相云水的云水路徑（反映含水量）、粒子平均尺度、粒子尺度的標準差（反映粒子尺度均勻性）來表征黑龍江省2015年云水部分微物理特性的時空分布。

在計算降水轉化率時，由于氣柱內總含水量和整層大氣可降水量在定義上是相同的，因此可按公式（1）計算出降水轉化率。需要指出的是，降水轉化率本身并非是一個與實際物理過程相對應的計算量，首先其分子是瞬時狀態量而其分母是累計量，其次其不能體現水汽及降水系統的輸送影響，但如果對其進行較長時間的平均，則其至少可以在一定程度上反映水汽向降水轉化在空間分布上的差異。因此文中討論的對象除地面降水外均為各個變量的季節平均量。

降水轉化率=日均整層大氣可降水量1地面24 h累計降水（1）

在分析云微物理特性與增雨潛力的關系時，不同類型的人工增雨作業所需要的云微物理條件不同，因此在利用MODIS云產品分析增雨潛力時應更有針對性地對液相和冰相的含水量、粒子平均尺度等反演產品與不同類型增雨條件的關系進行討論。除催化釋放潛熱所帶來的動力效應外，目前按照增雨催化劑種類可大致將增雨類型分為暖云催化和冷云催化[9]。其中暖云催化主要針對有一定液態含水量且云水粒徑較為均勻的云，這類云中由于云滴的碰并過程很弱導致云滴增長很慢，這時通過人工播撒吸濕性鹽核，可使云中迅速形成大滴并加速碰并過程，進而促進降水的產生。而冷云催化主要是通過播撒人工冰核或制冷劑使過冷水層中迅速產生人工冰晶，根據冰晶效應（貝吉龍）的原理，過冷云中水面飽和水汽壓大于冰面飽和水汽壓，導致云中的過冷水不斷消耗，并使冰晶不斷增長，進而加速啟動冷云降水過程，其必要條件是有一定的過冷水和缺乏冰晶。但MODIS云產品并不能在垂直方向上分辨云水分布，即不能給出過冷水含量，因此這里只能利用已知試驗進行推論，根據封閉系統中的冰水轉化規律[9]，冰晶和水滴的數密度比越小即冰晶相對水滴數量越少時，冰晶最終的尺度及所需增長時間均越大或越長，由此推論到衛星觀測及反演的情況，當有一定過冷水含量但冰晶數量少時，應能觀測到冰相水含量較少但冰相粒子尺度不均勻的情況。需要指出的是，上述變量雖然暫時無法找到精確的閾值進行定量判別，但在進行季節平均后，可以根據各個變量在空間分布上的相對大小進行分析。綜上所述，將MODIS云產品能夠反映的適宜增雨的微物理條件歸納為表1。

22015年黑龍江省降水量及轉化率

2.1年降水量概況黑龍江省山地與平原交錯分布（圖1a），其中西南的松嫩平原和東部的三江平原為主要的農牧產區，從西北延伸至中部的小興安嶺以及南部的張廣才嶺和老爺嶺山區則為主要的林業產區。從全年累計降水（圖1b）來看，黑龍江省2015年降水量為400～850 mm，其中從西北延伸至中南部的山區降水相對充沛，局地有700 mm以上的大值中心，而平原地區降水相對較少，三江平原西北部分地區及松嫩平原的大部分地區當年降水均不足550 mm，松嫩平原南部部分地區當年降水不足450 mm。

（b）從黑龍江省各季節的月均降水量（圖2）來看，2015年黑龍江省降水的季節差異很明顯，全省境內夏季月均降水量均在90 mm以上，而冬季（1、2、12月）全省月均降水量不足20 mm，春秋兩季降水介于冬夏兩季之間且整體偏少。其中春季松嫩平原大部分地區和秋季自松嫩平原到三江平原的大部分地區降水量均相對較少，月均降水量均小于30 mm。因此針對黑龍江省西南和東部2個平原地區的增雨潛力應重點討論。

2.2各季節的降水轉化率從圖3可看出，2015年黑龍江省各個季節的降水轉化率為0.05～0.30。在春季與秋季降水相對較少的2個平原地帶，降水轉化率與全省相比也相對較小，轉化率低值中心不足0.10，相比于其他地區0.15以上的降水轉化率，上述平原地帶在春秋兩季仍然具有相當的增雨潛力。鑒于冬夏兩季降水轉化率空間分布相對均勻且當年夏季降水充沛，冬季降水極少，在下文中主要針對黑龍江省春季和秋季上述相應平原地區的云微物理特征及增雨條件展開討論。

45卷11期趙麗斌結合MODIS云產品的黑龍江省2015年人工增雨潛力研究32015年黑龍江省云水的微物理特征及增雨條件分析

3.1春季云微物理特征由當年春季平均的液態云水含量（圖4a）可見，松嫩平原自西南向東北有一由少到多的液態云水含量的變化梯度（80～140 g/m2），其與該區春季降水量及降水轉化率均存在一定程度上的空間趨勢一致性，其西南部的轉化率和液態云水含量同時較小，表明水汽并沒有大量在此凝結為云水，因而春季松嫩平原并非大部分地區都適宜

進行暖云催化。而在松嫩平原北部（47°～48° N、125°～126° E），存在具有一定液態云水含量（120 g/m2左右）、液相云粒子尺度較小（11.5 μm左右）、液相粒子尺度標準差較小（4.5 μm以下）的區域（圖4a、c、e），在一定程度上符合表1所示的增雨微物理條件，表明春季松嫩平原北部存在適宜進行暖云催化的區域。而在三江平原中部同樣也存在液相粒子尺度較小且標準差較小的區域，且該區域的液態云水含量較松嫩平原更大（140 g/m2以上），因此春季三江平原中部雖然降水轉化率已達到全省平均（0.15左右），但仍然有繼續通過暖云催化增雨的潛力。而從春季云水的冰相微物理特性（圖4b、d、f）來看，并沒有明顯符合表1中所述的具有一定液水含量、冰相水含量較少且冰相粒子尺度標準差大的區域，因而并未發現春季全省境內有適合冷云催化的區域。

3.2秋季云微物理特征從秋季平均的液態云水含量（圖5a）來看，黑龍江省南部從西向東同樣存在著與降水量及降水轉化率小值區相對應的液態云水含量較小的區域，即并非整個南部地區均存在適宜增雨的云水條件，但在松嫩平原北部、三江平原東西兩側以及黑龍江省中部到南部的山區仍然存在具有一定液態云水含量（120 g/m2及以上）的區域，且從松嫩平原北部至三江平原西側（47°～48° N、125°～130° E）以及三江平原東部偏南區域存在著許多同時符合液態云水粒子尺度較小且液態云水粒子尺度標準差較小的局地性區域（圖5c、e），表明秋季這些區域符合暖元催化增雨的微物理條件。而從秋季平均的冰相微物理特征（圖5b、d、f）來看，在上述具有一定液態云水含量的松嫩平原北部和三江平原西側冰相云水含量較少（小于170 g/m2）且冰相粒子尺度標準差較大（大于8.5 μm），推斷這些區域很可能是有一定過冷水而缺乏冰晶的區域，因此很可能符合冷云催化的微物理條件并具有冷云催化增雨潛力。

4小結

利用歐洲中心再分析數據并結合MODIS衛星云產品，先后從降水轉化率和云水的微物理特征入手，對2015年黑龍江省不同季節的增雨潛力進行了研究，主要結論如下：

（1）在春季的松嫩平原大部分地區以及秋季的松嫩平原與三江平原大部分地區，降水相對全省而言均偏少，同時這些地區的降水轉化率也較小，表明在這些相對干旱的平原地帶仍然存在一定的增雨潛力。

（2）根據MODIS云產品中與云微物理特性有關的參量分析，在春季，松嫩平原北部和三江平原中部存在適宜進行暖云催化增雨的區域，但并沒有適合冷云催化增雨的區域。在秋季，從松嫩平原北部至三江平原西側存在許多適宜暖云催化增雨的區域，同時在松嫩平原北部與三江平原西側的局部地區存在適宜冷云催化增雨的區域。

參考文獻

[1] 齊彥斌，郭學良，金德鎮.一次東北冷渦中對流云帶的宏微物理結構探測研究[J].大氣科學，2007，31（4）：621-634.

[2] 洪延超，周非非.層狀云系人工增雨潛力評估研究[J].大氣科學，2006，30（5）：913-926.

[3] 陳祺，尹憲志，張龍，等.甘肅省空中云水資源分布及人工增雨潛力評估[J].現代農業科技，2015（1）：243-244.

[4] 韓輝邦，馬玉巖，康曉燕，等.應用CWR-PEP方法評估2012年青海省夏秋季人工增雨潛力[J].青海氣象，2015（1）：34-37.

[5] 趙玲，安沙舟，楊蓮梅，等.1976-2007年烏魯木齊可降水量及其降水轉化率[J].干旱區研究，2010，27（3）：433-437.

[6] 王維佳，周嵬，何奇瑾，等.飛機作業區域降水轉化率特征[J].安徽農業科學，2010，38（2）：18372-18373，18468.

[7] 代娟，黃建華，王華榮，等.襄樊市空中云水資源分布及人工增雨潛力研究[J].暴雨災害，2009，28（1）：79-83.

[8] 牛忠清，張苗苗，安英玉，等.2013年黑龍江省云水資源人工增雨潛力評估[J].安徽農業科學，2015，43（3）：251-253.

[9] 盛裴軒，毛節泰，李建國，等.大氣物理學[M].北京：北京大學出版社，2003：330-331，372-376.

[10] 劉貴華，余興，師春香，等.FY-3A/VIRR反演云微物理特征及與TERRA/MODIS反演結果的比較[J].高原氣象，2011，30（2）：461-470.

[11] 趙艷風，王東海，尹金方.基于CloudSat資料的青藏高原地區云微物理特征分析[J].熱帶氣象學報，2014，30（2）：239-248.

[12] 王洪強，陳勇航，彭寬軍，等.基于Aqua衛星總云量資料分析山區云水資源[J].自然資源學報，2011，26（1）：89-96.

[13] 宋松濤，張武，陳艷，等.中國西北地區近20年云水路徑時空分布特征[J].蘭州大學學報（自然科學版），2013，49（6）：787-793，798.

[14] 林丹.利用MODIS衛星產品分析西南地區云水特征[J].氣象科技，2015，43（1）：138-144.