吳思利,王 輝,鄭世超,曾朝陽,陳 翔,孫家正,馮利鵬
(上海衛(wèi)星工程研究所,上海 201109)
應(yīng)用多普勒信息的合成孔徑雷達(dá)(Synthetic Aperture Radar,SAR)不能直接提供目標(biāo)速度信息。經(jīng)過合成孔徑成像過程,不具有徑向速度的靜態(tài)目標(biāo)會在方位向上聚焦到零多普勒位置[1],如果直接利用多普勒分析方法得到速度信息則會極大損失頻譜進(jìn)而降低場景分辨率。作為合成孔徑原理的重要應(yīng)用,順軌干涉SAR 技術(shù)能夠利用具有時間差的完整分辨率二維相干數(shù)據(jù)得到場景圖像的多普勒偏移,并進(jìn)一步提取對應(yīng)的徑向速度,已廣泛應(yīng)用到包括洋流測速、地面動目標(biāo)檢測在內(nèi)的豐富情景中。順軌干涉(Along-Track Interferometry,ATI)技術(shù)由GOLDSTEIN 和ZEBKER 于1987 年提出,并首次在1989 年完成機載實驗[2],此后眾多機載系統(tǒng)[3]、航天飛機[4]和衛(wèi)星進(jìn)行了ATI 實驗或業(yè)務(wù)運行,獲得了優(yōu)于0.2 m/s 的測速效果。對ATI 的理論分析, ROMEISER 分析并建立了ATI-SAR 結(jié)合海洋表面流場的成像機理[5];有 學(xué)者[6-9]提出了針對調(diào)頻連續(xù)波ATI 系統(tǒng)的殘留相位差估計方法以及斜視情況下的精度分析方法等。
根據(jù)波束寬度與天線尺寸的對應(yīng)關(guān)系,單通道SAR 的順軌干涉無法滿足高增益和寬幅同時實現(xiàn),這一點對于弱散射特性目標(biāo)和大工作視角格外重要。SUESS 等[10]提出面向高分寬幅(High Resolu?tion Wide Swath,HRWS)的距離向多通道數(shù)字波束合成?掃描接收技術(shù)(Digital Beam Forming Scan-on-Receive,DBF-SCORE)可以有效解決這一問題。發(fā)射天線的距離向孔徑?jīng)Q定了單脈沖下的場景大小,減小天線高度獲得大的照射面積,再利用多個接收孔徑信號波束合成,獲得理想的觀測幅寬與信噪比。本文分析了DBF-SAR 在順軌干涉模式下的工作性能與誤差影響,通過指標(biāo)仿真實驗驗證了系統(tǒng)對寬幅場景的測速效果,彌補了單通道系統(tǒng)對海洋、河流等弱后向散射場景測速精度不高的缺點。
Ka 波段SAR 系統(tǒng),其應(yīng)用波長較短,與傳統(tǒng)SAR 系統(tǒng)相比實現(xiàn)相同的幅寬時只需要較小的天線尺寸,具有高分辨率、小型化、輕型化等特點[11-13]。但是由于Ka 波段系統(tǒng)損耗較高,如果采用傳統(tǒng)的單發(fā)單收工作方式,為滿足系統(tǒng)增益要求就需要非常大的峰值發(fā)射功率,這在工程實現(xiàn)上非常困難。要想在Ka 波段獲取高性能的實時處理,需要從多通道技術(shù)上尋找新的突破。
數(shù)字波束形成-掃描接收技術(shù)在發(fā)射端用一個小尺寸天線照射整個幅寬,如圖1 所示,接收端采用多通道接收天線等效獲得一個高增益窄波束依快時間沿距離向掃描,可以有效提高系統(tǒng)增益,獲得更低的距離模糊、更高的信噪比和更好的輻射分辨率,同時大幅度減小了數(shù)據(jù)下傳量。相比于單發(fā)單收系統(tǒng),基于數(shù)字波束合成(Digital Beam Form?ing,DBF)的Ka 波段SAR 系統(tǒng)對天線功率密度和硬件的要求降低,有利于降低成本,縮短研發(fā)周期,實現(xiàn)小型化。
圖1 距離向多通道合成結(jié)構(gòu)Fig.1 Structure of multi-channel DBF in range
忽略慢時間,只考慮快時間t、孔徑n接收信號sn(t)、加權(quán)系數(shù)ωn(t)、接收波束中心方向θ(t)、天線法線方向β、單通道距離向孔徑長度d,總的信號為各通道信號之和:
在t時刻(如圖2 所示),波束指向α(t)=θ(t)-β??鞎r間對應(yīng)的接收波束方向為理論波達(dá)角方向:
式中:Re為地球半徑;c為光速。
圖2 時變波達(dá)角Fig.2 Time-variant direction of arrival
彌補不同通道傳播歷程相位差的加權(quán)系數(shù)ωn(t)為
式中:λ為波長天線;dn為第n個接收通道相對天線中心的距離,
加權(quán)求和就相當(dāng)于形成一個時變的高增益窄波束在一個脈沖重復(fù)周期內(nèi)掃描整個測繪帶,跟蹤接收地面回波。通過DBF 過程,理想均衡通道回波信號經(jīng)過相干累加,而隨機噪聲期望不變,理論信噪比改善可以達(dá)到10 log10N。
在單通道原始回波信噪比為20 dB 下,多通道合成前后接收天線方向圖等效示意圖如圖3 所示,多通道合成與否脈沖壓縮信號信噪比改善情況示意圖如圖4 所示,信噪比改善值接近6 dB,符合理論計算。
圖3 DBF 接收前后增益Fig.3 Receive-gain before and after DBF
如圖5 所示,假設(shè)雷達(dá)平臺具有飛行高度H和速度Vp,配備順軌分布接收天線兩組并且相位中心A1和A2間距沿航跡長B,分別在距離向上具有多個接收通道。設(shè)定坐標(biāo)系xyz分別對應(yīng)順軌、交軌和高度向。
目標(biāo)相對雷達(dá)位置位于入射角θ、斜距r0處,并具有三向速度(vx,vy,vz),可以得到目標(biāo)徑向速度vr為
在理想情況下,目標(biāo)由速度帶來的干涉相位φv為
ATI 模糊速度是帶來2π rad 干涉相位的目標(biāo)速度:
式中:B為基線長度;Vp為雷達(dá)平臺運動速度。
系統(tǒng)理論上能不模糊觀測的速度范圍為-vamb/2 <vr<vamb/2。
對沿航跡干涉性能的主要限制因素進(jìn)行表述,并量化分析對于DBF-SAR 的影響。干涉流程之前需要分別成像,得到兩幅聚焦圖像。
1)沿航跡速度。如果速度分量vx大到不匹配方位向參考方程φa,目標(biāo)就會在方位向上散焦。使用二階近似,φa為
式中:ta為方位時間。
忽略徑向速度分量,目標(biāo)的相位歷程為
則在合成孔徑末端的殘余相位為
式中:Ta為合成孔徑時間,等于。因為飛行載具速度遠(yuǎn)大于地面目標(biāo)速度,如在800 km 軌道高度衛(wèi)星平臺運動速度可超過7 000 m/s,存在Vp?vx。要求φr小于π/2,得到對沿航跡速度約束|vx| ≤。
2)垂直航跡速度。目標(biāo)的徑向速度vr也會引起圖像散焦,并且在相位歷程上會附加額外的線性項,進(jìn)而改變多普勒中心頻率,引起動目標(biāo)方位向頻譜與雜波的分離。用QDop表示目標(biāo)方位向頻譜的偏移量:
目標(biāo)帶寬只有1-QDop的部分會被處理,所以方位向分辨率會以系數(shù)1/(1-QDop)降低,目標(biāo)單像素的信噪比會以系數(shù)(1-QDop)2降低。
3)距離徙動。由于徑向速度帶來的額外距離徙動也會造成動目標(biāo)在距離向上散焦,當(dāng)目標(biāo)的方位頻譜相對于PRF 混疊,相鄰脈沖間未補償?shù)木嚯x徙動量Δρn為
式中:nDop為目標(biāo)多普勒,相對于處理器多普勒偏移量PRF的整數(shù)倍。在合成孔徑時間內(nèi)的總額外距離徙動為
因此,快速目標(biāo)會難以識別。
徑向測速精度可以通過微分得到誤差子項貢獻(xiàn)量,其平方和即可表征整體干涉測速能力[14]。定義xyz坐標(biāo)系下x軸航向的一發(fā)多收DBF-ATI 系統(tǒng),A1和A2表述順軌兩組多通道接收天線且具有沿航向基線B,a為軌道高度H與地球半徑之和,θ為對目標(biāo)點觀測下視角。測速精度為
式中:誤差源σVP、σB、σΦ12、σBy、σBz、σh、σa、σr0分別為飛行器速度、沿航跡基線長度、干涉相位、基線y軸、基線z軸、觀測場景高度、軌道高度以及觀測斜距的測量不確定性,各誤差貢獻(xiàn)因子分別為
在多視數(shù)為NL時,干涉相位[15]可以表示為
式 中:s1、s2分別為A1和A2多通道 合成后 單視復(fù)信號。
干涉相位標(biāo)準(zhǔn)差是相干系數(shù)γ的函數(shù),為
總相干系數(shù)可以表述為子相干系數(shù)乘積:
式中:γtgt、γdop、γpro、γasr分別為目標(biāo)信噪比相干系數(shù)、空間相干系數(shù)、多普勒相干系數(shù)、處理相干系數(shù)、模糊相干系數(shù)。其中,目標(biāo)信噪受到距離向多通道影響較大,主要有通道固定相位誤差、空間去相干、波束指向誤差和波達(dá)角誤差4 種因素降低信噪比,影響合成信號質(zhì)量,進(jìn)而降低順軌干涉效果。
考慮同場景SAR 順軌圖像組,對應(yīng)分辨單元圖像值a1和a2分別為
式中:s1和s2為目標(biāo)信號;n1和n2為應(yīng)噪聲分量。并且,具有相差φv,
在理想條件下對兩組順軌接收通道完全一致,n1和n2為獨立噪聲,則對于該分辨單元的相干系數(shù)為
去除非相關(guān)交叉項后得到
化簡為
熱噪是非相干噪聲的主要元素。在實際應(yīng)用時,多通道合成存在多種誤差形式,以及固有的去相干效應(yīng)主要存在固定相位誤差、波達(dá)角估計誤差以及通道固有空間分布帶來的去相干。雷達(dá)截面積能夠表征觀測目標(biāo)的可探測性,在SAR 常用的觀測入射角30°下,水面一類場景的歸一化雷達(dá)截面積會比陸地場景弱至少30 dB。定義雷達(dá)回波中觀測目標(biāo)的雷達(dá)截面積為σtgt,雷達(dá)系統(tǒng)的等效后向散射系數(shù)為σ0,則得到取對數(shù)的沿視角目標(biāo)信噪比為
1)通道固定相位誤差。設(shè)第n個通道信號回波信號,經(jīng)過快時間補償后為不再受回波視角調(diào)制:
式中:sref(t)為參考通道信號;φn為各通道相對于參考通道的殘余固定相位誤差;Kr、t0為斜距r0雙程時延,則求和就能得到DBF 合成信號為
可以發(fā)現(xiàn),DBF 接收系統(tǒng)的各通道信號存在固定相位誤差時,不會直接影響單通道脈沖壓縮帶來的能量累積,但多通道合成后會在每個對應(yīng)距離單元形成對消。定義固定相位增益損耗為
2)波達(dá)角估計誤差。實際在地面具有高度起伏時,通過理想平地地球幾何獲取波達(dá)角,估計具有不可忽略的誤差,定義SCORE 法指向的方向為θs(t),真實具有高程h的目標(biāo)回波方向為θ0,則角估計誤差為
式中:t0為回波真實雙程時延,造成回波并沒有按照理想狀態(tài)被多通道合成最大增益接收,而是因為地形帶來了增益損失;G0為接收天線等效方向圖,定義損失量為
如圖6 所示,3 000 m 的地形高度誤差就能引起近端0.6°以上的波達(dá)角估計誤差,并相應(yīng)地引起5 dB 的增益損失。
3)空間去相干。在理想航跡下,通??梢院雎匝睾桔E干涉系統(tǒng)下空間去相干的影響,但對于多通道系統(tǒng),由于距離向多通道帶來的入射角差異,會造成雷達(dá)回波固有的信號差異,降低在空間尺度上的相干性。參考交軌干涉空間去相干系數(shù)[16],定義多通道下空間基線增益損失:
式中:ρy為距離分辨率,則通道間最大視角差異δθ約為。
則DBF-SAR 圖像信噪比為
系統(tǒng)參數(shù)見表1,分析多通道合成對ATI 測速性能影響。
圖6 地形對DBF 性能的影響Fig.6 The influence of topography on DBF performance
表1 仿真輸入Tab.1 Simulation inputs
應(yīng)用多通道合成后,干涉圖像對相干系數(shù)能夠獲得有效的提升。通過圖7 的單視仿真結(jié)果,可以得到雷達(dá)系統(tǒng)單通道數(shù)據(jù)在目標(biāo)歸一化散射截面積小于-20 dB 時相干系數(shù)會低于0.5,而此時四通道系統(tǒng)對應(yīng)值可以達(dá)到0.67。目標(biāo)歸一化散射截面積大于-10 dB 時,多通道相干系數(shù)提升極限小于0.1。
圖7 通道理想合成后相干系數(shù)Fig.7 Correlation coefficient by ideal multi-channels DBF
誤差源配置見表2,仿真得到多通道理想情況下的測速性能如圖8 所示。在當(dāng)前誤差分配下,可以發(fā)現(xiàn)干涉相位誤差為主要誤差量,單通道地距700 km 處誤差源貢獻(xiàn)量主要有相位誤差貢獻(xiàn)0.1 m/s、z軸基線誤差貢獻(xiàn) 量0.05 m/s、y軸基線誤差貢獻(xiàn)量0.05 m/s。
理想測速性能如圖9 所示。由圖9 可見:在地面弱散射場景下,幅寬內(nèi)單通道近端精度大于0.4 m/s;16 通道合成后,近端精度0.16 m/s 有效改善了分米級精度。強散射場景下,多通道測速精度集中在0.15 m/s 左右。
表2 誤差源輸入Tab.2 Error source inputs
圖8 測速誤差貢獻(xiàn)量Fig.8 Error contributions to velocity measurement
以不同的地面高程以及通道間最大相位誤差15°下進(jìn)行試驗,并考慮空間去相干可以得到如圖10所示的多通道合成后ATI 測速精度??梢园l(fā)現(xiàn)在地距較小時,高程誤差會嚴(yán)重影響測速精度,并且可能會出現(xiàn)實際系統(tǒng)性能不如單通道的情況。同理,當(dāng)通道間相位差過大,如150°時會發(fā)生性能不如單通道的情況。
圖9 理想測速性能Fig.9 Ideal performance of velocity measurement
圖10 非理想情景測速性能Fig.10 Non-ideal performance of velocity measurement
容易得出結(jié)論,DBF-SAR 面對弱后向散射場景可以提升系統(tǒng)的干涉測量性能。在具體實現(xiàn)方面,對多通道雷達(dá)系統(tǒng)相位誤差進(jìn)行實時或者后處理可以有效提高大多數(shù)情況下的系統(tǒng)順軌干涉性能。同時對于復(fù)雜地形環(huán)境,有必要進(jìn)行基于自適應(yīng)或者基于外部數(shù)字高程模型(Digital Elevation Model,DEM)數(shù)據(jù)的波達(dá)角估計,完成高精確度的快時間權(quán)校正。
作為順軌干涉應(yīng)用的有效提升手段,通過距離向多通道SAR 系統(tǒng)可以改善回波增益,進(jìn)而改善干涉相位質(zhì)量,提升測速能力。本文考慮了包括DBF-SCORE 原理、動目標(biāo)特性以及順軌干涉精度分析在內(nèi)的定量分析,并對多通道系統(tǒng)相位一致性、波達(dá)角估計以及空間去相干對合成后干涉效能的影響進(jìn)行了討論及仿真,為未來高分寬幅ATI 系統(tǒng)提供參考意義。