基于空間監(jiān)視雷達散射回波的電離層電子密度探測方法

丁宗華，代連東，許正文，吳健，楊嵩

（1.中國電波傳播研究所，山東青島266107；2.電波環(huán)境特性及?；夹g重點實驗室，山東青島266107）

基于空間監(jiān)視雷達散射回波的電離層電子密度探測方法

丁宗華1,2，代連東1,2，許正文1,2，吳健1,2，楊嵩1,2

（1.中國電波傳播研究所，山東青島266107；2.電波環(huán)境特性及?；夹g重點實驗室，山東青島266107）

目的提高空間監(jiān)視雷達的電離層電波環(huán)境適應性。方法VHF-L波段大功率雷達信號經(jīng)過電離層傳播時會發(fā)生非相干散射，散射回波中包含電離層信息。結(jié)合雷達方程和電離層目標特征，基于電離層非相干散射原理，理論分析空間監(jiān)視雷達用于電離層電子密度非相干散射探測的條件（包括雷達參數(shù)設置和波束掃描方式），給出電離層回波功率、自相關函數(shù)和電子密度表達式，利用Matlab編程對某大功率空間監(jiān)視雷達原始數(shù)據(jù)處理得到電離層散射回波及電子密度。結(jié)果給出了空間監(jiān)視雷達用于電離層電子密度非相干散射探測的基本條件，實測數(shù)據(jù)發(fā)現(xiàn)電離層回波和電子密度符合電離層變化特征。結(jié)論利用VHF-L波段大功率雷達空間散射回波探測電離層電子密度是可行的，為空間監(jiān)視雷達電離層環(huán)境感知與傳播自適應修正提供了一種可能的新途徑。

雷達；散射回波；電離層；電子密度；非相干散射

KEYWORDS:radar;scatter echo;ionosphere;electron density;incoherent scatter

利用大功率雷達探測空間目標時，若雷達信號穿過電離層（地面以上90～1000 km的區(qū)域）傳播，電離層中的電子、離子等將對入射電磁波產(chǎn)生非相干散射。通過對這種空間非相干散射回波的提取和分析，可以反演雷達信號傳播路徑上的電離層電子密度等參數(shù)，從而可能實現(xiàn)一種電離層探測新手段，并且用于大功率空間監(jiān)視雷達的電離層電波傳播自適應修正。Gordon等[1]從理論上提出了利用大功率雷達探測電離層非相干散射散射回波和物理參數(shù)的設想，隨后Bowles等[2]利用美國Long Branch雷達（工作頻率為41 MHz，峰值功率為6 MW）進行了首次電離層非相干散射探測試驗，驗證了大功率雷達探測電離層的原理可行性，隨后電離層非相干散射技術得到迅速發(fā)展[3—5]，目前國際上已建設了10多部非相干散射雷達[6]。

作為目前最強大的電離層地基探測手段，大功率非相干散射雷達在空間目標監(jiān)測方面也具有重要應用價值[7—8]，同樣，反過來其他大功率空間監(jiān)視雷達也可能用于電離層非干散射探測。20世紀80年代初中國電波傳播研究所[5]利用原7010相控陣空間目標監(jiān)視雷達，通過研制專門的信號采集與處理終端，成功檢測出電離層非相干散射回波[9—10]，并計算了電離層電子密度、電子溫度、離子溫度等參數(shù)。由于受當時雷達技術、信號處理和數(shù)據(jù)反演方法的限制，在信號處理方法實時性、探測分辨率、反演參數(shù)精度等方面存在許多不足，也未說明利用大功率雷達開展電離層非相干散射探測的條件與要求，無法滿足大功率空間監(jiān)視雷達電離層探測與自適應修正的需求。文中重點分析了利用大功率雷達開展電離層非相干散射探測的條件并介紹了有關的電離層參數(shù)反演方法，最后給出了某大功率雷達實測數(shù)據(jù)初步分析結(jié)果。

1 大功率雷達的電離層非相干散射探測條件

1.1 電離層非相干散射方程

電離層非相干散射是指電離層等離子體隨機熱起伏對入射電磁波產(chǎn)生的散射，由于散射截面與等離子體電荷質(zhì)量比緊密相關，因此實際上主要是電離層中自由電子產(chǎn)生的非相干散射。由于電離層電子散射截面非常小，散射信號很弱，因此一般需要大功率雷達才能夠探測到。電離層非相干散射信號主要來自電離層自由電子散射，它彌散分布在雷達波束內(nèi)，屬于分布式軟目標。設電離層探測目標的散射體積為：

式中：Pt為發(fā)射功率；A為天線有效面積；L為系統(tǒng)損耗；τ為脈沖寬度。與常見的硬目標雷達方程不同，電離層非相干散射回波信號功率與距離的平方成反比。

1.2 非相干散射信號的可檢測性和測量誤差

假設在300 km高度處散射，電離層電子密度為1012/m3（對應最強的散射回波），脈沖寬度為32 μs，散射點處雷達波束橫向面積為1 km2，則單個距離門內(nèi)總電子數(shù)目為5×1021，總散射截面為5×10－7，雷達天線有效口徑面積為500 m2，雷達發(fā)射功率為1 MW。代入式（2）可得，達到天線處的空間散射回波功率為2.5×10－16W。

另一方面，設接收機帶寬為500 kHz，雷達系統(tǒng)噪聲溫度為200 K，則天線接收的噪聲功率為1.38×10－15W。由此可見，在不采取相干積累條件下，最強的電離層散射回波還遠低于噪聲。因此，非相干散射探測一般應采用大功率發(fā)射、大口徑天線接收、低系統(tǒng)噪聲溫度、非相干積累、復雜編碼和信號處理等技術才能檢測到電離層散射回波。

電離層非相干散射測量誤差主要與回波信噪比、脈沖積累系數(shù)等有關，表達式為[3]：

式中：n為積累次數(shù)；Pn為背景電磁噪聲功率；m為噪聲測量點數(shù)。

一般來說，電離層F2層電子密度大，散射信號強，信噪比較高。圖1給出了非相干散射功率測量誤差隨積累次數(shù)和信噪比的變化，橫坐標表示脈沖積累次數(shù)，縱坐標表示相對測量誤差，不同曲線對應不同信噪比。從圖1可見，非相干散射回波測量誤差對信噪比的依賴性很強。當信噪比低于0.1時，功率測量誤差很大。當信噪比為0.01時，即使積累1萬次（探測時間約100 s），功率測量誤差可達100%。另外當探測積累次數(shù)在2000次以上時，非相干散射功率測量誤差基本不隨積累次數(shù)增加而發(fā)生改善。因此當電離層電子密度很弱，散射信號信噪比很低時，測量誤差會很大，單純增加積累次數(shù)不一定能改善測量誤差，還需結(jié)合其他手段，如采用靈活的波形設計等。

圖1 非相干散射測量誤差與信噪比、脈沖積累數(shù)的關系

1.3 對雷達參數(shù)和波束掃描方式的要求

1.3.1 工作頻率

一般要求雷達波長遠大于等離子體德拜長度，而電離層等離子體德拜長度為厘米量級甚至更小，所以頻率通常應小于S波段。但頻率不能太低，頻率越低，則背景電磁噪聲偏大，同時天線和雷達發(fā)射機等硬件設備更復雜。低于VHF波段時背景噪聲過大，而且電離層散射譜譜寬很窄，導致散射信號信噪比很低，因此非相干散射探測的頻率一般為VHF～L波段。

1.3.2 發(fā)射功率

由于電離層散射信號非常微弱，為了得到滿意的信噪比，通常要求發(fā)射峰值功率在兆瓦量級。

1.3.3 天線及掃描方式

非相干散射信號非常微弱，為了提高對微弱信號的檢測能力，一般采用大口徑拋物面天線和天線陣，此時天線有效面積至少數(shù)百平米，天線增益約40 dB或更大。由于電離層不是空間均勻分布的，經(jīng)常會出現(xiàn)各種尺度的不規(guī)則結(jié)構(gòu)。為了提高對小尺度結(jié)構(gòu)探測的橫向分辨率，要求具有較窄的波束寬度。為了減小旁瓣信號的干擾，要求天線的旁瓣電平比主瓣低得多。此外為了快速調(diào)整波束指向，以實現(xiàn)電離層三維空間掃描探測，要求天線具有一定的波束指向控制能力。

對于大功率空間監(jiān)視雷達，當采用波束凝視探測模式時，由于波束指向固定的電離層空間區(qū)域，此時與電離層非相干散射探測相似，只需對固定指向的電離層區(qū)域進行連續(xù)數(shù)千次以上（10 s以上）積累，即可得到電離層非相干散射回波。

當采用跟蹤掃描模式時，由于雷達波束指向快速變化，不同指向的電離層回波不一定相關，因此對雷達波束掃描速度和脈沖積累次數(shù)具有一定限制。如果雷達波束掃描速度為1(°)/s，則不同掃描時間對應的不同高度處的掃描弧長見表1。由于電離層具有一定的水平空間均勻性（一般認為水平方向50 km范圍內(nèi)電離層狀態(tài)均勻，散射回波具有相關性），則掃描10 s時300 km處的弧長已超過50 km，無法滿足水平均勻性條件。另一方面由于電離層至少在1 min之內(nèi)滿足時間均勻性(在1 min之內(nèi)，電離層狀態(tài)基本不變)。因此若在1 min內(nèi)，雷達波束始位于同一相關區(qū)域，則可將此時間和相關區(qū)域內(nèi)的散射回波進行相干積累，以便獲取電離層非相干散射回波。

表1 不同掃描時間、不同高度處的掃描弧長

1.3.4 系統(tǒng)噪聲溫度

非相干散射探測除了采用高發(fā)射功率和高增益天線外，還應最大限度地降低系統(tǒng)噪聲溫度，通常要求小于300 K。

1.3.5 脈沖寬度與波形

為了得到較強的散射回波，一般要求脈沖寬度盡可能寬，至少為數(shù)百微秒。此時自相關函數(shù)時延范圍可達數(shù)百微秒，對應的非相干散射功率譜分辨率為千赫茲量級，滿足電離層非相干散射要求。脈沖寬度不能過寬，否則在同樣占空比條件下，脈沖重復周期太長，脈沖重復頻率太低，速度模糊嚴重。為了獲得高時寬帶寬積，一般采用相位編碼脈沖，可同時獲得較寬的脈寬寬度和較高的距離分辨率。

常見的硬目標不同，電離層為分布式軟目標。為了適應不同探測高度上對距離分辨率和功率譜分辨率的要求，同時為了盡可能壓縮旁瓣，提高信噪比，需采用復雜的編碼。常用的脈沖波形為交替碼、巴克碼、互補碼等。其中交替碼是非相干散射雷達專用相位調(diào)制碼。

2 信號處理與電離層反演方法

設雷達空間散射回波信號采樣數(shù)據(jù)為z，發(fā)射脈沖為m，則積累后的電離層回波功率為：

式中：?表示卷積運算；Pn為噪聲功率；i表示距離門（一般從100～700 km）的序號；N為相干積累次數(shù)；j表示積累次數(shù)的序號；cj為某次積累時的系數(shù)。

雷達采樣數(shù)據(jù)的自相關函數(shù)表達式[3]：

由非相干散射雷達方程得到電離層電子密度表達式為[3,9,11]：

式中：k與雷達系統(tǒng)參數(shù)，如發(fā)射峰值功率、脈沖寬度、天線增益和系統(tǒng)損耗等有關，一般需利用其他手段（如電離層垂直探測儀探測的最大電子密度）來標校；α為電離層等離子體德拜長度與雷達波長之比；Tr為電離層電子溫度與離子溫度之比。由于德拜長度為厘米量級甚至更小，而對VHF-L波段頻率來說，雷達波長大于15 cm，因此2α一般很小（約0.001），可以忽略。電離層的電子與離子溫度比隨著高度變化，一般約1～3，需利用散射功率譜反演的電子溫度和離子溫度計算獲得。若假定電子離子溫度比為常數(shù)（通常在夜間滿足）則式（4）可進一步簡化。

3 實測數(shù)據(jù)分析

2015年1月使用某P波段大功率空間環(huán)境監(jiān)測雷達[5]進行電離層非相干散射探測實驗。雷達參數(shù)：拋物面天線口徑為29 m、峰值功率為2 MW、平均功率為100 kW、工作頻率為500 MHz、天線增益為41 dB、系統(tǒng)損耗為3.5 dB，采用凝視探測模式，發(fā)射脈沖為13位巴克碼，脈沖重復周期為12 ms、脈寬為390 μs、碼元寬度為30 μs。圖2給出了該雷達空間散射回波的原始采樣，可見空間散射回波電平隨機變化，無顯著的電離層高度分層變化特征。利用1024次（約12.3 s）和10240次（約123 s）脈沖重復周期的數(shù)據(jù)分別進行相干積累，在普通臺式計算機上利用Matlab語言編程處理得到了電離層散射回波功率剖面和電子密度剖面，耗時為數(shù)十分鐘，如圖3和圖4所示?？梢钥闯?，電離層電子密度剖面具有先增加后降低的趨勢，在約250 km具有最大電子密度，此外在150～200 km具有疑似分層結(jié)構(gòu)，這些都符合電離層基本變化規(guī)律。相對于圖3，圖4的回波信噪比更高，電子密度變化曲線更光滑，精度相對更高。這里利用Matlab編程處理的時效性很低，為了滿足實時的電離層探測與修正的需求，可從以下幾個方面提高數(shù)據(jù)處理速度和實時反演：采用C語言等底層語言編程，因為Matlab編程執(zhí)行效率低，占用內(nèi)存高，耗時長；采用高性能的計算設備；進一步優(yōu)化數(shù)據(jù)處理算法和程序代碼。

圖2 大功率雷達空間散射回波的原始采樣

圖3 電離層散射回波功率和電子密度剖面（積累1024次）

圖4 電離層散射回波功率和電子密度剖面（積累10240次）

4 結(jié)語

詳細分析了利用大功率雷達開展電離層電子密度非相干散射探測的基本方法，并就雷達參數(shù)和天線掃描方式設置進行了討論，最后給出了某大功率雷達實際采集數(shù)據(jù)及其電離層回波功率剖面和電子密度剖面，初步驗證了以上方法的原理可行性。這對利用大功率空間監(jiān)視雷達開展電離層探測進而用于雷達目標定位的電波折射修正，從而提高目標定位精度提供了一種新思路，同時可能為電離層空間環(huán)境探測提供一種新手段，具有重要的科學意義和應用價值。

[1]GORDON W E.Incoherent Scatter of Radio Waves by Free Electrons with Applications to Space Exploration by Radar[J].Proc IRE,1958,46(11):1824—1829.

[2]BOWLES K L.Observation of Vertical Incidence Scatter from the Ionosphere at 41Mc/sec[J].Phys Rev Lett, 1958(1):454—455.

[3]EVANS J V.Theory and Practice of Ionosphere study by Thomson Scatter Radar[J].Proc IEEE,1969,57(4): 496—530.

[4]EVANS J V.Millstone Hill Thomson Scatter Results for 1965[J].Planetary Space Sci,1970,18(8):1225—1253.

[5]丁宗華,吳健,許正文,等.電離層非相干散射雷達探測技術應用展望[J].電波科學學報,2016,31(1):193—198.

[6]丁宗華,代連東,董明玉,等.非相干散射雷達進展:從傳統(tǒng)體制到EISCAT_3D[J].地球物理學進展,2014, 29(5):2376—2381.

[7]MANDEVILLE J C,RIBONI F,BLELLY P L. Interpretation of EISCAT RADAR data for Orbital Debris Study[J].Advances in Space Research,1995,16(11):29—33.

[8]丁宗華,許正文,趙振維.非相干散射雷達的空間碎片參數(shù)提取[J].現(xiàn)代雷達,2014,36(2):1—5.

[9]鄭傳青,景玉山,倪濟安,等.利用非相干散射雷達探測電離層[J].電子學報,1984,12(6):11—15.

[10]鄭傳青,朱太平,景玉山.利用R10雷達探測電離層電子濃度?電子和離子溫度分布[J].電波科學學報,1986, 1(1):11—19.

[11]丁宗華,魚浪,代連東,等.曲靖非相干散射雷達功率剖面的初步觀測與分析[J].地球物理學報,2014, 57(11):3564—3569.

Ionospheric Electron Density Measurement Based on Scattered Echo from the Space Surveillance Radar

DING Zong-hua1,2,DAI Lian-dong1,2,XU Zheng-wen1,2,WU Jian1,2,YANG S ong1,2
(1.China Research Institute of Radio-wave Propagation,Qingdao 266107,China 2.National Key Laboratory of Electro-Magnetic Environment,Qingdao 266107,China)

ObjectiveTo improve the environmental adaptation of ionospheric radio from space surveillance radar.MethodsWhen the VHF-L radio propagated through the ionosphere,the radar signal was incoherently scattered by the ionospheric free electron and the ionosphere information was included in the scattered echo.Combined with the radar equation and ionospheric target property,the radar parameter setting and beam scanning mode were investigated by the incoherent scatter mechanism of ionosphere.The ionospheric echo power,auto-correlation function and electron expression were introduced.And the raw sampling data from one space radar with high power was processed by using the Matlab routine and the ionospheric echo and electron density were obtained.ResultsConditions for the ionospheric incoherent scatter from the space surveillance radar was described.ionospheric echo and The electron density varied in accordance with the ionospheric characteristic.ConclusionIt is feasible to derive the ionospheric electron density from the scattered echo of space surveillance radar with high power in VHF-L band,which might be a new method for ionospheric environment measurement and propagation correction.

10.7643/issn.1672-9242.2017.07.005

TJ06；TN011

A

1672-9242(2017)07-0024-05

2017-03-30；

2017-04-22

裝備預研基金項目

丁宗華（1978—），男，湖北人，博士，主要研究方向為電離層空間環(huán)境探測。