GNSS-R海面風速反演技術研究

2020-05-23 06:37:24劉原華何孟然牛新亮

全球定位系統 2020年2期

劉原華,何孟然, 牛新亮

(1. 西安郵電大學通信與信息工程學院,陜西西安 710121;2. 中國空間技術研究院西安分院,陜西西安 710100)

0 引言

全球衛星導航系統反射信號技術(GNSS-R)是一種利用接收機接收物體表面的反射信號,并從反射信號提取所需特征的新興的地球遙感技術. Auber等[1]在1994年首次提出了在機載實驗中探測到了海面散射全球導航系統信號;Zavorotny等[2]在雙基雷達方程的基礎上利用Kirchhooff近似的幾何光學方法,建立了海面散射信號的時延-多普勒二維相關功率模型,即Z-V模型;之后美國NASA與科羅拉多大學聯合開展了一系列實驗,為海面風場的反演奠定了理論基礎[3-6]. 對于星載GNSS-R技術,國內在實測數據處理及地表參數反演的相關研究較少.

利用星載GNSS-R反演海面風速時,由于幾何因素的影響,功率校準后的雙基地雷達橫截面積不能直接用來反演風速. 本文以Z-V模型為理論基礎,對功率校準后的數據進行幾何校正,以獲得高精度的數據進行風速反演.

1 反演原理

GNSS-R散射信號相關功率是時延-多普勒的二維函數

(1)

2 幾何校正

星載接收機接收到信號的功率受溫度的影響,因此要將接收到的信號進行功率校準[7],得到雙基雷達橫截面積(BRCS)[8]

(2)

式中:Pg為接收到的散射信號的功率;C為接收到的信號總功率;CN為校準開關指向天線的噪聲功率;PB為儀器接收到的散射信號功率;Pr為儀器產生的熱噪聲功率;CB為校準開關指向黑體負載的噪聲功率.

圖1 反射信號場景圖

利用功率校準后的歸一化散射信號相關功率進行風速反演[9], 如圖1所示,與真實風速的光滑曲線相比,歸一化相關功率的曲線較為粗糙, 反演風速和實際風速相差較大,這是由于發射機發射功率、天線方向圖和幾何關系等的影響,為了得到高精度的數據,在估算地球物理參數之前,必須對非地表相關項進行校正.

GNSS-R產生的數據通常采用時延-多普勒相關功率圖(DDM)的形式,DDM圖像代表了在特定的時延和多普勒頻移下,海洋表面散射信號相關功率的分布情況. 由于DDM峰值功率受到鏡面反射點以外多種因素的影響,例如熱噪聲、斑點噪聲、幾何形狀和天線方向圖的不均性等,使用DDM峰值功率無法準確地預測鏡面反射點的位置. 反射區域的反射信號中路徑延遲最短的點定義為鏡面反射點[10],在鏡面反射點周圍構建一個大型的網格,使用DTU10 平均海面模型[11]對鏡面反射點進行校正. 為了提供合理的海面空間分辨率,反演海面風速時,通常僅使用DDM的近鏡面部分[12]. 圖2為功率校準后的DDM圖.

圖2功率校準后的DDM圖像

2.1 雙基雷達橫截面

使用時延-多普勒頻點上幾個變量的有效值來簡化等式(1),有效值包括時延和多普勒擴展函數Λ和S

(3)

由于DDM的散射截面積σ不是歸一化的,需要變換公式(3)將歸一化散射雷達橫截面積σ0替換為散射橫截面積σ,得到散射橫截面積σ表達式為

(4)

圖3 雙基雷達橫截面積

2.2 有效散射區和物理散射區

單個時延-多普勒區包含捕獲到的海洋表面一個或多個物理區域上的散射功率,每個DDM的物理面積可以表示為

(5)

GPS的模糊函數[8]增加了每個時延-多普勒頻段的有效面積,導致散射功率從幾何確定的物理散射區域散布到相鄰的時延-多普勒頻段中,改變了總的散射功率,DDM中的每個時延-多普勒頻段的物理面積和有效區域都會變化. 因此將每個時延-多普勒頻段的有效表面散射面積表示為模糊函數加權表面積分,即:

(6)

圖4 有效面積

2.3 歸一化雙基地雷達橫截面

將雙基雷達橫截面積和有效散射區的DDM相結合,計算出歸一化雙基雷達橫截面積σ0. 近鏡面部分的DDM功率的平均值稱為DDM平均值(DDMA),典型的DDMA是DDM圖像中以鏡面反射點為中心的5(多普勒)×3(時延)的區域. 由于衛星攜帶儀器中開環信號跟蹤器存在的誤差,這些實際的DDMA區域通常不會與儀器生成的DDM區域完全對齊,采用鏡面反射點周圍的DDM區域與DDMA邊緣區域加權組合的方式計算實際的DDMA總雷達散射截面積σ0. 計算DDMA雙基雷達的總橫截面積為