地空路徑統(tǒng)計雨衰減建模的幾個問題

林樂科 趙振維 盧昌勝 張鑫 李娜

(中國電波傳播研究所 電波環(huán)境特性及?；夹g(shù)重點實驗室,青島 266107)

?

地空路徑統(tǒng)計雨衰減建模的幾個問題

林樂科 趙振維 盧昌勝 張鑫 李娜

(中國電波傳播研究所 電波環(huán)境特性及?；夹g(shù)重點實驗室,青島 266107)

地空路徑雨衰減統(tǒng)計特性是衛(wèi)星通信等地空無線電系統(tǒng)設(shè)計和運行評估的關(guān)鍵參數(shù). 針對地空路徑統(tǒng)計雨衰減的建模研究,分析了影響模型準(zhǔn)確度的幾個關(guān)鍵因素,它們包括雨滴譜分布、雨頂高度、降雨率在水平和垂直路徑的不均勻性等. 同時比較分析了現(xiàn)有主要模型計算的路徑調(diào)整因子隨仰角的變化、與頻率的關(guān)系等,討論了現(xiàn)有模型的局限性和可能的改進途徑. 通過上述對比研究,提出了地空路徑雨衰減統(tǒng)計建模向確定性模型方向發(fā)展的觀點.

雨衰減;衛(wèi)星通信;地空路徑;降雨率

DOI 10.13443/j.cjors.2015081001

引 言

互聯(lián)網(wǎng)和信息時代的發(fā)展使得地空無線電系統(tǒng)應(yīng)用的頻段向微波的高端和毫米波段擴展,以衛(wèi)星通信系統(tǒng)為例,從Ku、Ka到極高頻(Extremely High Frequency,EHF)頻段都開始逐步地應(yīng)用. 降雨衰減是導(dǎo)致信號中斷的首要傳播因素,也是10 GHz以上頻段地空傳播首要的考慮因素,它是由于雨滴對電磁波的吸收與散射作用產(chǎn)生的,通常雨衰減隨頻率和降雨率的增加而增大. 雨衰減統(tǒng)計特性是衛(wèi)星通信等地空無線電系統(tǒng)可靠性設(shè)計和運行評估的關(guān)鍵參數(shù). 長期以來,隨著雨衰減數(shù)據(jù)積累,雨衰減預(yù)報模型得到持續(xù)改進. 1997年Dissanayake等[1]基于對數(shù)正態(tài)分布由相似原理推導(dǎo)出DAH(Dissanayake-Allnutt-Haidara)模型,經(jīng)由國際電信聯(lián)盟無線電通信部(International Telecommunication Union-Radiocommunication sector,ITU-R)測試驗證其準(zhǔn)確度好于當(dāng)時其他模型,并被ITU-R P.618建議[2]采納,作為國際通用標(biāo)準(zhǔn)沿用至今. 但隨著國際電聯(lián)雨衰減數(shù)據(jù)庫積累數(shù)據(jù)的增多以及國際上進行的熱帶地區(qū)雨衰減測試等,逐步發(fā)現(xiàn)了該模型一些不合理的地方以及在熱帶地區(qū)預(yù)報偏差較大等問題,有待繼續(xù)改進. 2003年英國提出了全概率的地空路徑雨衰減模型[3],其同時提出的雨衰減率的計算公式修訂了ITU-R P.839建議書. 最近十年來,巴西、韓國等多個國家也向國際電聯(lián)提交了新的地空路徑雨衰減預(yù)測模型[4-5],國際電聯(lián)正在比較與測試中. 其中我們提出的基于降雨率調(diào)整因子新概念的全概率雨衰減模型[6-7](以下簡稱推薦模型)代表中國向國際電聯(lián)提交,其綜合性能在目前的比較中為最優(yōu). 由于降雨時空分布的復(fù)雜性,地空路徑雨衰減統(tǒng)計建模方法以半經(jīng)驗公式居多,除各國提交ITU的方法外,還有其他學(xué)者提出的方法,如Stutzman 等提出的simple model等[8];也有采用純經(jīng)驗方法的,例如,楊紅衛(wèi)等[9]也提出了基于神經(jīng)網(wǎng)絡(luò)計算雨衰減的方法,由于其沒有解析公式不方便進行比較,基本屬于純經(jīng)驗?zāi)Ｐ?其主要缺點是沒有顯性公式不方便進行比較,同時其模型參數(shù)完全依據(jù)實驗數(shù)據(jù),適用范圍有限,缺乏對不同頻率、不同降雨氣候區(qū)拓展的物理基礎(chǔ). ITU-R P.311建議書[10]對于雨衰減模型除了要求精度高以外,還要求預(yù)報模型有較好的物理基礎(chǔ)和簡單的形式. 此外本文考慮的是通用的雨衰減模型,對于只適用于局部氣候區(qū)域的模型不進行分析.

雨衰減建模涉及的主要因素包括雨衰減率的計算(其中最重要的因素為雨滴尺寸分布,即雨滴譜分布)、降雨率在水平和垂直路徑不均勻分布的等效處理、雨頂高度等,本文對于影響雨衰減建模精度的主要因素進行了討論,同時分析了幾個主要模型計算的路徑調(diào)整因子隨仰角的變化、與頻率的關(guān)系并討論了其合理性. 最后討論了地空路徑雨衰減建模研究的發(fā)展方向.

1 地空雨衰減建模的一般過程

圖1為地空路徑雨衰減示意圖.圖中, A為一般不產(chǎn)生衰減的冰晶層,B為雨頂,C為雨區(qū),D為地空路徑,其它主要參數(shù)如下:

hs: 地面站的平均海拔高度

θ: 路徑仰角

hR: 雨頂海拔高度

Ls: 電波信號穿越雨區(qū)的斜路徑長度

LG:Ls對應(yīng)的水平路徑

圖1 地空路徑雨衰減示意圖

地空路徑統(tǒng)計雨衰減一般可以由下式計算:

Ap=γp·LE.

(1)

式中:

Ap為給定時間概率p%的統(tǒng)計雨衰減;

LE為等效路徑長度(降雨率調(diào)整因子也可轉(zhuǎn)換為等效路徑長度),用以調(diào)整降雨率分布在水平路徑和垂直路徑的不均勻性;

γp為給定時間概率p%的雨衰減率,即假定電波穿過降雨率均勻分布的單位距離(通常為1 km)產(chǎn)生的衰減,影響雨衰減率最重要的因素是雨滴尺寸分布,即雨滴譜分布.

根據(jù)上面的分析可知,地空雨衰減建模的一般過程如下:首先采用一定的雨滴譜分布依據(jù)降雨率的統(tǒng)計分布計算雨衰減率,根據(jù)雨頂高度和路徑仰角確定穿越雨區(qū)的電波路徑長度,然后建立路徑調(diào)整因子(降雨率調(diào)整因子)模型,去調(diào)整降雨率分布在水平路徑和垂直路徑的不均勻性. 對于全概率模型,通過上述過程建立的式(1)即可獲得對應(yīng)時間概率的雨衰減;對于概率轉(zhuǎn)換模型,一般先采用式(1)計算給定時間概率(通常為0.01%)的雨衰減,然后建立概率轉(zhuǎn)換公式得到其他時間概率的雨衰減. 下面對地空雨衰減建模的主要因素進行分析,由于全概率模型是雨衰減統(tǒng)計模型的發(fā)展趨勢,對于概率轉(zhuǎn)換公式不展開討論.

2 地空雨衰減建模的主要因素

2.1 雨滴譜與雨衰減率

雨衰減率即電磁波通過1 km均勻降雨路徑后的衰減,它是雨衰減預(yù)報的基礎(chǔ). 雨衰減率與雨滴散射特性及雨滴譜分布有關(guān),由于降雨為稀疏介質(zhì),采用單次散射近似,則雨衰減率可由下式計算:

A=4.343×103∫Qt(D)N(D)dD dB/km.

(2)

式中:N(D)dD為雨滴譜,即dD直徑間隔的雨滴數(shù);Qt(D)為單個雨滴的消光截面. 雨滴消光截面可由雨滴散射特性獲得,雨滴散射特性隨電磁波頻率的不同和雨滴形狀的變化而采用不同的計算方法. 當(dāng)雨滴尺寸參數(shù)x=ka?1(k為傳播常數(shù),a為雨滴等效半徑)時,可以采用瑞利散射計算雨滴的散射特性,例如C波段降雨散射特性即可由瑞利散射計算. 當(dāng)x=ka?1時,雨滴散射特性可以采用幾何光學(xué)的方法計算. 當(dāng)x=ka接近1時,若雨滴為球形,可以由Mie理論獲得其散射特性的嚴(yán)格解;實際雨滴一般為非球形,通?？梢越茷楸鈾E球形,在這種情況下,需采用數(shù)值方法求得近似解,其中點匹配法和T矩陣法應(yīng)用最為廣泛.

雨滴譜分布是指對應(yīng)不同降雨率,不同尺寸的雨滴在空間的分布狀況. 雨滴譜分布是研究降雨特性、雷達氣象和無線電波傳播的重要參數(shù),得到了廣泛的測量和模式化研究. 雨滴譜在全世界范圍內(nèi)得到廣泛的測量,在不同地區(qū)測量數(shù)據(jù)的基礎(chǔ)上得到了大量雨滴尺寸分布模型,如:Laws-Parsons(L-P)分布、Marshall-Palmer(M-P)分布、Joss et al分布、對數(shù)正態(tài)分布、gamma分布等,其中L-P分布和負指數(shù)分布的M-P分布及Joss et al分布被廣泛使用,ITU-R的雨衰減率模式即是使用擬合的L-P分布得到的.

采用經(jīng)驗?zāi)Ｊ接嬎阌晁p率可以避免理論計算的繁瑣,同時應(yīng)用較為簡便. 雨衰減率計算中最為常用的經(jīng)驗?zāi)Ｊ綖?

γR=kRαdB/km.

(3)

式中:R為降雨率(mm/h);γR為對應(yīng)的雨衰減率;k與α為與頻率、仰角、極化角有關(guān)的參數(shù),其是通過選定一定的雨滴譜計算其對應(yīng)的雨衰減率然后擬合得到的.

從上面的分析可以看到,雨衰減率計算的關(guān)鍵因素是雨滴譜分布.目前的ITU-R P.838-3建議書中的雨衰減率模型是由英國專家提出的利用改進的L-P分布計算的雨衰減率擬合得到的. 其在擬合過程中利用了P.618建議書的統(tǒng)計雨衰減模式(以下簡稱ITU雨衰模型)和國際電聯(lián)地空雨衰減數(shù)據(jù)庫的數(shù)據(jù). 由于ITU雨衰模型本身也存在一定的誤差,必然影響雨衰減率模型的建立. 趙振維[11]通過對相近降雨率的雨滴尺寸分布求和平均及大雨滴截斷和對大雨滴的平滑處理,得到青島、廣州和新鄉(xiāng)三地區(qū)的雨滴譜分布. 利用負指數(shù)模型對測量數(shù)據(jù)的擬合結(jié)果表明,擬合值與M-P分布有較好的一致性,其計算的雨衰減率在毫米波段與838-3模式有較大誤差,與M-P分布結(jié)果比較接近. 圖2為不同雨滴譜雨衰減率的比較,可以看到在40 GHz以下頻段,除gamma分布以外,其他雨滴譜分布計算結(jié)果與實測雨滴譜計算結(jié)果接近;當(dāng)頻率大于40 GHz以后,gamma分布的結(jié)果逐漸與實測雨滴譜分布的結(jié)果接近,L-P分布、Weibull分布、對數(shù)正態(tài)分布的結(jié)果偏差逐漸增大,M-P分布雖在頻率的高端誤差增大,但總體上與實測雨滴譜的結(jié)果最吻合.

最近在熱帶和亞熱帶地區(qū)的測試[12-14]表明,熱帶和亞熱帶地區(qū)的雨滴譜分布計算得到的雨衰減率與838-3模式有較大差距. 在2015年ITU-R第三研究組的會議上,奧地利與歐洲太空局(European Space Agency,ESA)(簡稱歐空局)的專家聯(lián)合提出了利用本地雨滴譜分布計算雨衰減率的重要性[15],并推薦建立各地區(qū)雨滴譜分布的數(shù)據(jù)庫,其比較的數(shù)據(jù)包含溫帶地區(qū),說明即使在溫帶地區(qū)雨衰減率特征也不能由838-3建議書統(tǒng)一表征. 圖3給出了利用本地雨滴譜計算雨衰減率與838-3模式的誤差,可以看出來有較大的差異,其中Gan站的差異一般在2～3個dB以上.

圖2 不同雨滴譜雨衰減率的比較(水平極化)

圖3 本地雨滴譜雨衰減率與838-3模式的比較[15]

2.2 雨頂高度

雨頂高度是計算穿越雨區(qū)斜路徑長度的關(guān)鍵參數(shù). ITU-R最新的雨頂高度模式[17]中將-2 ℃等溫層高度作為雨頂高度,這實際上是考慮了融化層的效應(yīng). 其計算公式為:

hR=h0+0.36 km.

(4)

式中h0為0 ℃等溫層高度， ITU-R給出精度為1.5°×1.5°的0 ℃層高度柵格數(shù)據(jù)，其它任意地點0 ℃層高度數(shù)據(jù)由雙線性插值得到， 這里雨頂高度僅與地理位置有關(guān).

無線電氣象學(xué)上通常將降雨分為層狀云降雨(stratiform precipitation)和對流云降雨(convective precipitation)兩類. 層狀云降雨一般范圍較廣,雨強小,持續(xù)時間長,其雨頂高度可采用平均雨頂高度;對流云降雨范圍小,雨強大,持續(xù)時間短,由于強烈的垂直對流作用,雨頂高度往往伸展得很高,在-2 ℃層以上存在過冷雨滴,并且其雨頂高度與降雨率有關(guān). Stutzman[8]推薦了一種雨頂高度與降雨率的關(guān)系式:

(5)

該模式有一定的合理性,但由于缺乏實測數(shù)據(jù)與其對照,在地空傳播預(yù)測中尚未得到廣泛應(yīng)用.

2015年,法國、意大利和歐空局的專家向國際電聯(lián)提交了關(guān)于新的全球統(tǒng)計雨頂高度的提案[18],其中包含了逐月的雨頂高度統(tǒng)計.

從上面的分析可以看出,雨頂高度與降雨類型、降雨強度、季節(jié)、溫度等多種因素有關(guān),需要在雨衰減統(tǒng)計建模中認真考慮.

2.3 降雨在路徑上不均勻分布的處理

由于實際降雨是不均勻的,不僅空間上不均勻,時間上也不均勻,因此斜路徑上的統(tǒng)計降雨率與點降雨率統(tǒng)計是不相同的,所以用點降雨率統(tǒng)計預(yù)測斜路徑雨衰減時不能直接利用降雨路徑長度,這時需要一個等效路徑長度反映其實際關(guān)系,即:

(6)

式中: Rp為給定時間概率p%的統(tǒng)計不超過的降雨率; Le為等效路徑長度; Ls為經(jīng)過雨區(qū)的實際路徑長度； r為路徑調(diào)整因子. 等效路徑長度是指電磁波通過一段降雨率為R的均勻雨區(qū),其產(chǎn)生的衰減等于實際路徑衰減,則稱波經(jīng)過均勻雨區(qū)的斜路徑長度為等效路徑長度. 目前主要的雨衰減模式都采用這種概念. 我們[19]提出了降雨率調(diào)整因子的概念:

Ap=k(rpRp)αLs.

(7)

式中rp為給定時間概率p%的降雨率調(diào)整因子,可以看到其可以轉(zhuǎn)換為路徑調(diào)整因子:

(8)

降雨率空間分布上的不均勻不僅表現(xiàn)在水平方向,也表現(xiàn)在垂直方向,韓國的測試發(fā)現(xiàn)降雨率在垂直方向有較大的不均勻分布[20]. ITU雨衰模型的路徑調(diào)整因子由水平路徑調(diào)整因子和垂直路徑調(diào)整因子組成,也有的模型是將其綜合考慮成一個路徑調(diào)整因子來處理. 通常情況下,路徑調(diào)整因子(降雨率調(diào)整因子)的模型都是在假定路徑存在一定數(shù)量的雨胞且雨胞內(nèi)降雨率服從某一特定分布的前提下推導(dǎo)出來的,然后利用雨衰數(shù)據(jù)庫的數(shù)據(jù)擬合其參數(shù)得到的.

3 雨衰減模型的合理性分析

對于已建立的雨衰減模型,除比較其精度之外,還需要驗證其在物理意義上的合理性. 例如模型預(yù)報結(jié)果是否存在奇異值或不連續(xù)性,預(yù)報結(jié)果隨頻率或仰角等的變化是否與物理規(guī)律以及測試數(shù)據(jù)相符等. 下面我們以雨衰減隨仰角的變化為例進行分析.

通常人們認為,對于同一個站點仰角越低的傳播鏈路,由于電波穿越雨區(qū)的路徑較長,產(chǎn)生的雨衰減越大,隨著傳播仰角的增大,電波穿過雨區(qū)的距離減小,此時降雨造成的衰減也應(yīng)隨之降低. 也就是說同一站點隨仰角的增加降雨衰減應(yīng)單調(diào)遞減. 但是,2003年在巴西召開的ITU-R會議上,有專家指出ITU雨衰模型在仰角逐漸增大時雨衰減預(yù)報值先減小后增大[21]的問題,我們稱其為ITU雨衰模型隨仰角變化的奇異性現(xiàn)象. 國際上一直比較重視熱帶和亞熱帶地區(qū)的降雨衰減預(yù)報. 這些地區(qū)的地球站與同步衛(wèi)星的傳播仰角一般較大. 此外,對于利用非同步軌道衛(wèi)星進行通信、探測(如合成孔徑雷達)等功能的地空無線電系統(tǒng),由于衛(wèi)星位置的變動,經(jīng)常會發(fā)生高仰角的情況. 因此雨衰減在高仰角下的準(zhǔn)確預(yù)報對地空無線電系統(tǒng)的高性能工作非常有意義.

根據(jù)公式(6),等效路徑長度可以表示為

(9)

由前面的分析可知,降雨率調(diào)整因子可以轉(zhuǎn)換為路徑調(diào)整因子. 而由式(9)可知,對于給定降雨率和頻率的傳播鏈路來說,等效路徑長度決定著路徑降雨總衰減. 因此我們只須分析等效路徑長度隨仰角的變化即可. 利用式(9)和國際電聯(lián)地空雨衰數(shù)據(jù)庫的數(shù)據(jù)得到如圖4所示的等效路徑長度與傳播仰角的分布圖. 可以看到,等效路徑長度隨傳播仰角的增大而單調(diào)減小. 因此試驗數(shù)據(jù)并不具有隨傳播仰角的增大,降雨衰減先減小后增加的奇異性變化.

圖4 等效路徑長度與傳播仰角的變化關(guān)系

我們設(shè)定以下鏈路參數(shù)比較不同模型計算的雨衰隨仰角的變化:

雨衰年平均時間概率: 0.01%

0.01%時間不超過的降雨率: 88 mm/h

地球站的緯度: 22°N

地球站的經(jīng)度: 114°W

雨頂高度: 5.21 km

頻率: 30 GHz

極化: 水平

參與比較的模型為ITU雨衰模型和推薦模型. 比較結(jié)果如圖5所示. 由圖5和圖4對比可以看出,推薦模型對于給定的觀測點降雨率,等效路徑長度隨傳播仰角的增大而單調(diào)減小. 因此其不存在隨仰角變化的奇異性現(xiàn)象,而ITU雨衰模型當(dāng)仰角在60°左右時,路徑調(diào)整因子出現(xiàn)隨仰角增大而增大的奇異現(xiàn)象. 此外目前向國際電聯(lián)提交的其他國家的地空雨衰模型也都存在不同程度的這種奇異現(xiàn)象,限于篇幅,不再列舉.

(a) ITU雨衰模型(ITU model)

(b) 推薦模型(proposed model)圖5 不同模型雨衰減隨傳播仰角的變化

4 討論與結(jié)論

根據(jù)上面的分析可以看到,制約地空路徑雨衰減統(tǒng)計模型精度的主要因素包括:雨衰減率的計算其中關(guān)鍵是雨滴譜分布的準(zhǔn)確度;雨頂高度;準(zhǔn)確并且符合物理規(guī)律的路徑調(diào)整因子(降雨率調(diào)整因子的構(gòu)造). 就實驗觀測而言,目前的觀測手段要較以前更豐富,由于衛(wèi)星信標(biāo)頻段的擴展,地空路徑雨衰減的測試頻段也由Ku頻段逐步擴展到Ka頻段和EHF等毫米波頻段,對降雨的觀測由傳統(tǒng)的翻斗式雨量計改進為激光雨滴譜儀、微雨雷達、多普勒天氣雷達等,可以提供實際的雨胞分布、降雨在水平和垂直路徑的不均勻分布等多種信息,用于雨衰減率和路徑調(diào)整因子的建模. 由于降雨和雨衰監(jiān)測手段的日益齊全和豐富,將有助于推動地空路徑雨衰減統(tǒng)計模型向確定性模型方向發(fā)展(類似點對面?zhèn)鞑ヮA(yù)測的P.1812建議書[22]),例如,建立能反映各地不同雨滴譜分布的雨衰減率計算模式;發(fā)展能反映降雨氣候和隨降雨率大小變化的雨頂高度模型;建立能準(zhǔn)確反映降雨在水平和垂直路徑不均勻分布的路徑調(diào)整因子(降雨率調(diào)整因子);等等. 為便于工程應(yīng)用、模型改進和誤差分析,模型通常推薦采用顯式公式建立.

另外,ITU雨衰模型采用的是利用0.01%時間雨衰減換算其他時間概率雨衰減的方法,這種方法會帶來幾個問題:一是對于不同的降雨氣候區(qū)其0.01%時間雨衰減可能相近,但在其他不同時間概率可能有較大差別;二是其沒有考慮有雨和無雨概率分點的問題. 隨著降雨率統(tǒng)計的精細化和準(zhǔn)確化,全概率的雨衰減預(yù)測模型成為新的趨勢.

[1] DISSANAYAKE A, ALLNUTT J, HAIDARA F.A prediction model that combines rain attenuation and other propagation impairments along earth-satellite paths[J]. IEEE transaction on antennas propagation, 1997, 45(10): 1546-1558.

[2] Rec. ITU-R P.618-11. Propagation data and prediction methods required for the design of earth-space telecommunication systems[S]. Geneva: ITU, 2013.

[3] ITU-R Doc. 3M/29. Proposed modification to recommendation ITU-R P.618-8: an improved model of rain attenuation prediction along earth-space path[R]. Geneva: ITU, 2003.

[4] ITU-R Doc. 3M/27. Prediction of rain attenuation in slant paths and terrestrial links using the full rainfall rate distribution[R]. Geneva: ITU, 2008.

[5] ITU-R Doc. 3J/89. Proposed revision of the recommandation ITU-R P.618-11: a new approach for the effective path-length model for rain attenuation based on rain-cell characteristics: global model[R]. Geneva: ITU, 2014.

[6] ZHAO Z W, LU C S, LIN L K. ITU-R doc.3J/16 and 3M/31: modification to rainfall rate adjustment factor for modelling and prediction methods of rain attenuation statistics[R]. Geneva: ITU, 2012.

[7] 趙振維, 盧昌勝, 林樂科, 等. 地空鏈路雨衰減預(yù)測模式研究[J]. 電波科學(xué)學(xué)報, 2013, 28(3): 412-419.

ZHAO Z W, LU C S, LIN L K, et al. Prediction model of rain attenuation for the earth-space links[J].

Chinese journal of radio science, 2013, 28(3): 412-419. (in Chinese)

[8] STUTZMAN W L, DISHMAN W K. A simple model for the estimation of rain-induced attenuation along earth-space paths at millimeter wavelengths[J]. Radio science, 1982, 17(6): 1465-1476.

[9] 楊紅衛(wèi), 何晨, 諸鴻文, 等. 基于進化神經(jīng)網(wǎng)絡(luò)的地空路徑雨衰減模型[J]. 通信學(xué)報, 2001, 22(9): 20-27.

YANG H W, HE C, ZHU H W, et al. An earth-space rain attenuation model based on evolutionary neural network[J]. Journal of china institute of communications, 2001, 22(9): 20-27.(in Chinese)

[10]Rec. ITU-R P.311-14. acquisition, presentation and analysis of data in studies of tropospheric propagation[S]. Geneva: ITU, 2013.

[11]趙振維. 水凝物的電波傳播特性與遙感研究[D]. 西安:西安電子科技大學(xué), 2001.

ZHAO Z W. Study on radiowave propagation characteristics and remote sensing of hydrometeors[D]. Xi’an: Xidian University, 2001.(in Chinese)

[12]DAS S, MAITRA A, SHUKLA A K. Rain attenuation modeling in the 10-100 GHz frequency using drop size distributions for different climatic zones in tropical India[J]. Progress in electromagnetics research B, 2010, 25: 211-224.

[13]MANDEEP J S, NG Y Y, ABDULLAH H, et al. The study of rain specific attenuation for the prediction of satellite propagation in Malaysia[J]. Journal of infrared millimeter & terahertz waves, 2010, 31(6): 681-689.

[14]CHEN K S, CHU C Y, TZENG Y C. A semi-empirical model of rain attenuation at ka-band in northern Taiwan[J]. Progress in electromagnetics research M, 2011, 16: 213-223.

[15]ITU-R Doc.3J/131& 3M/225. Document on location dependent variations of drop size distributions[R]. Geneva: ITU, 2015.

[16]ITU-R Doc.3J/138& 3M/237. Improvement of the model of rain height in rec. ITU-R P.839-4 from ERA interim profile and surface data[R], Geneva: ITU, 2015.

[17]Rec. ITU-R P.839-4. rain height model for prediction methods[S]. Geneva: ITU, 2013.

[18]ITU-R Doc.3J/138&3M/237. Improvement of the model of rain height in Rec. ITU-R P.839-4 from ERA interim profile and surface data[R]. Geneva: ITU, 2015.

[19]趙振維, 盧昌勝, 林樂科, 等. 雨衰減預(yù)測中的降雨率調(diào)整因子概念[J]. 電波科學(xué)學(xué)報, 2013, 28(1): 33-38.

ZHAO Z W, LU C S, LIN L K, et al. The conception of rainfall rate adjustment factor for the rain attenuation prediction[J]. Chinese journal of radio science, 2013, 28(1): 33-38.(in Chinese)

[20]ITU-R Doc.3J/122&3M/215. Contribution to ITU-R study group 3 databanks: rain cell statistics from meteorological radar measurement data in Korea[R]. Geneva: ITU, 2015.

[21]ITU-R Doc. 3M/17. Comments on recommendation ITU-R P.618-7[R]. Geneva: ITU, 2003.

[22]Rec ITU-R P.1812-3. A path-specific propagation prediction method for point-to-area terrestrial services in the VHF and UHF bands[R]. Geneva: ITU, 2013.

林樂科 (1972-),男,山東人,研究員,博士,長期從事對流層傳播與環(huán)境特性研究.

趙振維 (1965-),男,河北人,研究員,博士,中國電子科技集團公司第二十二研究所總工程師,ITU-R SG3研究組副主席,長期從事電波環(huán)境及其傳播特性的研究和國際標(biāo)準(zhǔn)的制定.

盧昌勝 (1983-),男,甘肅人,高級工程師,西安電子科技大學(xué)博士研究生,主要從事毫米波、太赫茲波大氣傳播特性方面的研究.

Issues on modelling rain attenuation statistics on Earth-Space links

LIN Leke ZHAO Zhenwei LU Changsheng ZHANG Xin LI Na

(China Research Institute of Radiowave Propagation, National Key Laboratory ofElectromagneticEnvironment,Qingdao266107,China)

The rain attenuation statistics are key parameters for the proper planning and evaluating of earth-space radio systems such as satellite communication systems. The key factors affecting the accuracy of the rain attenuation prediction model on slant path are analyzed herein, which include rain drop size distribution(DSD), rain height, the non-uniformity of rain rate in horizontal and vertical direction, etc. Meanwhile, the dependencies of main rain attenuation models on elevation and frequency are compared and analyzed, and the limitations of the existing models and the approaches of improving the models potentially are discussed. Through the investigation of the above comparison, an opinion is put forward that the tendency of rain attenuation statistics modeling on earth-space path is to the path-specific model.

rain attenuation; satellite communication; earth-space path; rain rate

10.13443/j.cjors.2015081001

2015-08-10

TN011

A

1005-0388(2016)03-0616-07

林樂科, 趙振維, 盧昌勝, 等. 地空路徑統(tǒng)計雨衰減建模的幾個問題[J]. 電波科學(xué)學(xué)報,2016,31(3):616-622.

LIN L K, ZHAO Z W, LU C S, et al. Issues on modelling rain attenuation statistics on earth-space links[J]. Chinese journal of radio science,2016,31(3):616-622. (in Chinese). DOI: 10.13443/j.cjors.2015081001

聯(lián)系人： 林樂科 E-mail: llk22s@163.com