太赫茲雷達散射截面的仿真與時域光譜測量

張旭濤 闕肖峰 蔡禾 孫金海 張景 李糧生 劉永強

1)(電磁散射重點實驗室,北京 1000854)

2)(電子科技大學電子工程學院,成都 610054)

1 引 言

太赫茲波(頻率0.3—10 THz,波長1—30 μm)是位于微波與紅外之間的一段特殊電磁波段,長期以來,由于缺乏有效的太赫茲波發射和探測裝置而一直處于未開發狀態,從而形成了電磁波研究中的“空隙”,近年來,隨著技術的發展,這一波段已經在醫療衛生、航天航空、基礎科研等領域有了廣泛的應用.太赫茲在目標特性領域的主要優勢有: 相對于激光,太赫茲波的穿透力更強,能夠穿透煙霧、沙塵等,可用于復雜環境; 相對于微波,帶寬更大,帶來更高的目標分辨率,高分辨圖像能夠獲得目標更豐富的特征信息,更好地獲取目標的精細結構和微動特征; 在反隱身方面有極大的應用潛力[1,2].

目前,美、歐等國已建立了多個微波倍頻太赫茲波目標特性實驗室.其中比較典型的是美國馬薩諸塞大學羅威爾分校(University of Massachusetts Lowell)所屬的亞毫米波技術實驗室(Submillimeter-Wave Technology Laboratory,STL),研究人員已經建立了多套連續太赫茲波實驗系統(例如: 0.24,0.35,0.52,0.585,1.56 THz等頻段)[3,4],并對這些頻率的雷達散射截面(radar cross section,RCS)等目標特性進行了深入研究.國外基于時域光譜技術的太赫茲頻段目標特性研究同樣起步較早,2000年,美國俄克拉荷馬州立大學的研究人員利用其搭建的太赫茲時域光譜雷達系統對3 mm直徑氧化鋁圓柱體的散射進行了研究,在0.2—1.4 THz范圍內,RCS的測量結果與PO理論解具有非常好的一致性[5].2010年,丹麥科技大學基于時域光譜測量系統獲得了1.1 THz等頻點處的1︰150縮比F-16飛機模型的360°全方位向RCS結果,并對縮比模型進行了成像[6].2012年,德國太赫茲中心的研究人員通過光纖耦合太赫茲時域光譜系統測量了200 GHz等頻點處金屬球、圓柱等定標體RCS,與理論解比較測量值達到了較小的誤差,并測量了1︰250縮比旋風飛機模型的180°方位向RCS[7].

近來,國內國防科技大學對太赫茲目標特性進行了相關的研究[8]; 中國工程物理研究院于2011年采用自主研制的太赫茲器件構建了國內首個0.14 THz高分辨率逆合成孔徑雷達成像實驗系統[9],2013年又完成了670 GHz全固態成像雷達[10].

當前國內的研究成果多基于太赫茲低頻段連續波系統,而對基于太赫茲時域光譜技術的RCS測量一直缺少準確的定量研究,時域光譜技術目前能夠覆蓋非常大的頻率范圍(2 THz以上),在如此寬的頻域范圍內要對目標的RCS進行準確的測量具有非常大的難度.本文基于太赫茲時域光譜技術,搭建了時域光譜散射測量系統,具備0.1—1.6 THz目標RCS測量及成像能力.在系統性能測試的基礎上[11],測量了幾種典型定標體的RCS,與理論結果誤差在3 dB以內,驗證了系統的測量能力.然后測量了幾種復雜目標的RCS,并將表面積分方程法與多層快速多極子算法相結合對目標的雷達散射截面進行了數值計算,經測量與仿真結果對比,0.2—1.6 THz范圍內誤差均在3 dB以內.本文的工作解決了時域光譜寬頻段(0.2—1.6 THz)RCS定量測量的難題,對后續目標RCS測量等夯實了基礎; 測量結果與理論仿真結果具有非常好的一致性,驗證了測量與數值計算方法的準確性.

2 實驗系統

基于太赫茲時域光譜技術建設了一套太赫茲時域光譜目標散射測量系統,該緊縮場系統可以實現遠場測試環境,系統后向散射角為6.5°,具有目標準后向散射信號的測量能力,測量系統如圖1所示.太赫茲波由光電導天線產生后經緊縮場系統擴束準直后入射到目標處,目標散射的太赫茲波經收集緊縮場光路進入接收天線,接收天線采集到的數據由鎖相放大器處理最終得到原始時域回波數據.有機玻璃密閉罩將整個太赫茲波傳輸光路密封在濕度可控的環境中.測量時,系統濕度控制在3%以下,以最大限度地減少太赫茲波傳輸中的空氣損耗.系統搭建完成后,我們對系統的各項性能進行了測試,其中對目標區太赫茲波強度的測量結果如圖2所示,圖中顯示了在垂直于入射波方向的二維平面上的太赫茲波強度分布,呈明顯的高斯分布,由此可確定靜區的范圍.

圖1 太赫茲時域光譜散射測量系統光路圖Fig.1.Scattering measurement system of terahertz time domain.

上述區域場強峰值變化較大,場強均勻區域較小,對于直徑超出該區域的的金屬球仍然可測得其準確的RCS.原因如下: 測試所用的金屬球直徑在厘米量級,電尺寸超過20,處于光學區,此時金屬球的散射中爬行波可忽略,主要為鏡面反射,而金屬球鏡面反射區域小于場強均勻區域,保證了測量結果的準確.

3 目標RCS仿真與測量

3.1 定標體

采用相對定標法測量目標的RCS,公式如下:

式中,σcal(f)為定標體在頻率f處的理論RCS,Pobj(f)和Pcal(f)分別為目標和定標體時域回波測量值經傅里葉變換后在頻點f處的回波功率值.實際計算時,由于時域光譜系統的頻譜分辨率較低,本文測試時將信號頻譜分辨率統一處理為10 GHz,為消除個別頻點數據的偶然性,(1)式中頻點f處的功率由以該頻點為中心頻率、帶寬40 GHz范圍內的功率的平均值得到.

圖2 目標區太赫茲波場強分布圖(單位: μV)Fig.2.Distribution of the field intensity(unit: μV).

金屬球、圓盤、圓柱的RCS理論精確解分別由(2),(5)及(6)式得出[12].

=jn(x)+iyn(x)為第一類球漢克爾函數;jn(x)為 第一類球貝塞爾函數;yn(x)為第二類球貝塞爾函數.

式中r是圓盤半徑,φ是入射角.

式中a,l分別是圓柱地面半徑及柱高,φ是入射角.

水資源管理水平得到提高。國家最嚴格水資源管理制度試點建設通過水利部中期評估。水資源論證納入政府投資項目聯合審批，規模以上地下水用戶全部實現水量在線監測。推行再生水、淡化海水、生態水和地下水“三增一減”管理，年利用再生水1.4億t、淡化海水3 000萬t，中心城區年生態補水1.73億m3，壓采深層地下水1 200萬m3。順利通過國家節水型城市復查，完成水功能區跨界斷面核定，全市所有水功能區納入水質監測。全年用水總量控制在24億m3以內，萬元GDP用水量降至18 m3，萬元工業增加值用水量降至8.3 m3，重要江河湖泊水功能區達標率提高到21%。

此外,為實現導體目標太赫茲頻段散射特性的數值仿真,本文采用表面積分方程法(surface integral equation method,SIE)結合阻抗邊界條件(impedance boundary condition,IBC)建模技術進行分析[13].

在太赫茲頻段導體目標的電導率不再趨于無限大,而是隨頻率變化.導體目標可以認為是有耗介質,表面波阻抗增加.此時導體目標邊界上滿足阻抗邊界條件,其表面電場與磁場的關系為:

其中η0表示自由空間波阻抗,η0為相對表面阻抗,為導體表面單位外法向矢量.

根據等效原理J=,M=,推導出表面電流和磁流的關系:

從而建立表面積分方程進行求解.

在太赫茲頻段采用數值法進行電磁建模時往往面臨電大尺寸問題,本文采用多層快速多極子算法(multilevel fast multipole algorithm,MLFMA)對矩陣方程加速計算[14,15].MLFMA利用幾何分組,將子散射體(基函數)之間的耦合關系分為近區和遠區關系.近區組內子散射體的耦合通過嚴格的矩量法(moment of method,MoM)直接計算;遠區組內的子散射體耦合關系則利用加法定理將格林函數在譜域內進行展開:

其中i和j為場源基函數編號,m和m'為場源分組編號,k為自由空間波數,為轉移因子.最后通過聚合-轉移-配置過程來實現.將傳統MoM的計算量和存儲量從量級減少為O(NlogN)量級,其中N為數值計算未知量數.MLFMA采用迭代算法求解矩陣方程,應用了廣義最小殘差方法(generalized minimal residual,GMRES)并結合預條件技術加速收斂.GMRES每次迭代只需進行一次矩陣矢量相乘運算,具有較穩健的收斂性能; 預條件則利用MLFMA的附近組阻抗矩陣元素構造預條件矩陣,能夠顯著提高迭代求解的計算效率.

圖3 直徑50 mm金屬球RCSFig.3.RCS of the sphere with the diameter of 50 mm.

圖4 直徑20.9 mm金屬圓盤RCSFig.4.RCS of the disk with thediameterof 20.9 mm.

圖5 金屬圓柱RCS定標結果Fig.5.RCS of the cylinder with the diameter of 20 mm and height of 20 mm.

測量過程中,以直徑35 mm的金屬球作為定標體,測量直徑50 mm的金屬球、直徑20.9 mm的金屬圓盤以及底面直徑20 mm柱高20 mm的金屬圓柱,直徑50 mm的金屬球RCS測量結果及理論值如圖3所示,直徑20.9 mm的金屬圓盤測量結果如圖4所示,底面直徑20 mm柱高20 mm的金屬圓柱測量結果如圖5所示.目標由金屬鋁制作,表面粗糙度均小于0.3 μm.

3.2 復雜目標

在0.2—1.6 THz范圍內定標體的RCS測量誤差小于3 dB,由此驗證了系統及測量技術的有效性.為進一步測量復雜目標的RCS,我們制作了帶有細微結構的邊長30 mm的正方形復雜樣件,目標材質為金屬鋁,表面粗糙度均小于0.3 μm,如圖6和圖7所示.測量時,由于系統的雙站角為6.5°,因此目標姿態為垂直放置,入射角 3.25 °.對目標建模計算時,由于目標加工有一定隨機性,會造成模型與實物有一定的誤差.在建模時也對實物幾何尺寸進行了多次測量,盡量減小目標的加工帶來的誤差.計算時部分頻點處材料的復介電常數如表1所列,其余頻點處的數值由擬合得到.

普通定標體在太赫茲頻段的散射中心相對單一,主要為鏡面散射中心.而凹腔等結構的多次復雜散射是實際情況中普遍存在的,因此對此類目標的測量就更具現實意義.

圖6 貫穿孔結構平板實物Fig.6.Photos of the plate with through-holes.

圖7 螺釘結構平板實物Fig.7.RCS of the plate with screws.

表1 不同頻點處復介電常數Table 1. Complex permittivity in different frequency.

由于時域光譜系統時域(空域)的超高分辨力,我們能夠得到目標細微結構的回波,如圖8所示.與一維距離像相對應,目標的各處散射中心對應的回波在時域信號中清晰可辨,并可根據回波峰值的位置精確計算出目標不同散射面的幾何結構數據.兩類復雜目標的RCS理論計算及測量結果如圖9及圖10所示.復雜目標的RCS測量結果顯示0.2—1.6 THz范圍內測量與理論值誤差在3 dB以內,表明系統對復雜目標的RCS測試精度同樣在3 dB誤差內.

圖8 螺釘結構平板時域回波Fig.8.Time domain echo of the plate with screws.

圖9 貫穿孔結構平板RCS測量值與理論值Fig.9.RCS of the plate with through-holes.

圖10 螺釘結構平板RCS測量值與理論值Fig.10.RCS of the plate with screws.

4 結 論

太赫茲波長在亞毫米量級,而太赫茲時域光譜系統可以通過控制延遲線的掃描精度,從而在時域回波上獲得非常高的分辨率,達到微米量級,因此能夠識別目標的極細微結構.本文測試的帶有細微結構的復雜目標代表了軍事目標表面的螺釘(鉚釘)、小孔等結構,這正是太赫茲時域光譜技術在目標特性領域的重要應用.在對系統進行了充分的優化調試后,系統的動態范圍、線性度、底噪、頻率穩定性等均達到測試要求,降低了系統本身帶來的測試誤差.在放置目標時采用激光定位從而保證了定位誤差最小化.充分降低這些誤差才能保證最終的測量誤差達到最小.與微波、毫米波導體目標建模不同的是太赫茲頻段導體目標的電導率不再趨于無限大,而是隨頻率變化.導體目標可認為是有耗介質.本文基于時域光譜系統對幾種典型定標體及復雜目標在0.2—1.6 THz的RCS進行了測量,并采用表面積分方程法進行建模仿真,多層快速多極子算法加速計算得到了目標太赫茲頻段RCS的理論數值結果,測量結果與理論結果誤差在3 dB以內.3 dB的誤差來源主要有: 一是目標的定位誤差,由于太赫茲波的指向性非常好,導致目標回波對目標姿態極為敏感; 二是靜區的波束均勻性,當前對靜區波束的精確整形存在較大困難,導致靜區波束幅度起伏較大,增加了測量誤差.

當前時域光譜系統存在的難題在于,對1.6 THz以上的頻段存在測量不準的問題,主要原因是受限于當前太赫茲收發天線的功率,在高頻區動態范圍不足.另外,對于時域光譜系統,其不同頻率的太赫茲波波束寬度不同,頻率越高,波束寬度越小,且時域光譜系統的太赫茲波束為高斯線型,RCS測量可用波束范圍很小,如何將高斯波束整型成為更大的測量可用波束是今后的研究重點之一.