梁 浩，李國剛，徐艷國

(南京電子技術(shù)研究所,南京 210039)

?

飛行高度對(duì)雷達(dá)對(duì)海性能影響性分析評(píng)估

梁 浩，李國剛，徐艷國

(南京電子技術(shù)研究所,南京 210039)

回顧了點(diǎn)目標(biāo)和面雜波的經(jīng)典雷達(dá)方程和接收機(jī)噪聲基本原理，介紹了三種經(jīng)典海雜波單位面積反向散射率經(jīng)驗(yàn)?zāi)Ｐ?。使用?jīng)驗(yàn)?zāi)Ｐ头治隽瞬煌w行高度、特定海情條件海雜波分布特性。使用實(shí)際試飛數(shù)據(jù)印證了實(shí)錄機(jī)載海雜波與經(jīng)驗(yàn)?zāi)Ｐ偷内厔?shì)一致性。結(jié)合海雜波背景中目標(biāo)檢測(cè)的特點(diǎn)，提出了一種快速評(píng)估飛行高度對(duì)海雜波背景下目標(biāo)檢測(cè)的影響性分析評(píng)估方法。最終給出了機(jī)載海面監(jiān)視雷達(dá)使用高度的選擇建議。

海面監(jiān)視雷達(dá)；海雜波；飛行高度；海情海況

0 引 言

機(jī)載海面監(jiān)視雷達(dá)目標(biāo)檢測(cè)面臨的主要背景是海雜波，海雜波特性與諸多因素相關(guān)，包括：海情海況、風(fēng)速、風(fēng)向、雷達(dá)照射方向、浪涌傳播方向、載機(jī)飛行高度等[1-2]。上述因素中，對(duì)載機(jī)飛行高度影響海雜波的關(guān)注較少，然而，其影響又客觀存在，不容忽視。如何在相應(yīng)海情海況和風(fēng)向條件下，應(yīng)用本文提出的快速分析評(píng)估方法，在特定海情條件下，選擇合適的飛行高度，獲得最佳目標(biāo)檢測(cè)效果，是本文研究和分析的重點(diǎn)內(nèi)容。

1 海面目標(biāo)檢測(cè)原理

由于實(shí)際工作環(huán)境中，海情海況變化莫測(cè)，遇到完全風(fēng)平浪靜的理想海面概率較低，海面回波(海雜波)廣泛存在。海面監(jiān)視雷達(dá)在海面目標(biāo)檢測(cè)過程中，無法回避海雜波存在性問題。

隨著探測(cè)距離的延伸，海雜波將逐漸衰減，并最終淹沒在接收機(jī)噪聲背景中。由于檢測(cè)背景不同，因此，必須將目標(biāo)檢測(cè)的信號(hào)處理分為雜波區(qū)和噪聲區(qū)進(jìn)行考慮。

在海雜波背景中，信號(hào)處理從回波中檢測(cè)目標(biāo)需要對(duì)抗海雜波回波，關(guān)注指標(biāo)應(yīng)為信雜比(SCR)，即目標(biāo)回波功率與雜波回波功率的差異程度。

在純?cè)肼暠尘爸?，信?hào)處理從回波中檢測(cè)目標(biāo)需要對(duì)抗接收機(jī)噪聲，關(guān)注指標(biāo)為信噪比(SNR)，即目標(biāo)回波功率與接收機(jī)背景噪聲功率的差異程度。

1.1 點(diǎn)目標(biāo)雷達(dá)方程

點(diǎn)目標(biāo)雷達(dá)回波方程，給出了特定距離確定雷達(dá)反向散射截面積(RCS)的點(diǎn)目標(biāo)回波功率理論值的計(jì)算公式[3]

(1)

式中：Pr為雷達(dá)接收目標(biāo)回波功率；Pt為雷達(dá)發(fā)射峰值功率；G為雷達(dá)天線增益；σT為目標(biāo)平均RCS值；D為脈沖壓縮比；λ為雷達(dá)工作波長(zhǎng)；R為目標(biāo)存在距離；Ls為各類系統(tǒng)損耗和傳輸損耗總和。

1.2 面雜波雷達(dá)方程

面雜波雷達(dá)回波方程給出了特定距離面雜波(地物雜波或海面雜波)的回波功率理論值計(jì)算公式[3]

(2)

式中：Pr為雷達(dá)接收面雜波回波功率；Pt為雷達(dá)發(fā)射峰值功率；G為雷達(dá)天線增益；c為光速；τ為脈沖壓縮后脈沖時(shí)寬；θaz為方位波束寬度；σ0為單位面積反向散射因子；D為脈沖壓縮比；λ為雷達(dá)工作波長(zhǎng)；θel為雷達(dá)波束擦地角；R為雜波單元距離；Ls為各類系統(tǒng)損耗和傳輸損耗總和。

1.3 接收前端噪聲功率

雷達(dá)接收機(jī)的射頻前端的熱噪聲是系統(tǒng)噪聲的重要組成部分，其功率理論計(jì)算公式[3]為

PN=k×T0×B×Fn

(3)

式中：k為玻耳茲曼常數(shù)；T0為基準(zhǔn)溫度290K；kT0=4×10-21W/Hz；B為接收機(jī)工作帶寬；Fn為考慮非理想接收機(jī)(實(shí)際工程化接收機(jī))所引入的額外噪聲影響，引入的噪聲系數(shù)。

可使用的簡(jiǎn)便計(jì)算公式[4]為

PN=-114+Fn+10lgB

(4)

式中：PN單位為dBm；B單位為MHz。

1.4 目標(biāo)信雜比和信噪比

在系統(tǒng)處理能力可承載目標(biāo)和虛警數(shù)量條件下，通過設(shè)定不同的檢測(cè)門限(差異程度比較)，完成目標(biāo)CFAR檢測(cè)處理。

在雜波背景中，目標(biāo)信雜比可表述為

SCR=PTr-PCr

(5)

式中：PTr為點(diǎn)目標(biāo)回波功率，單位為dBm，可通過式(1)進(jìn)行計(jì)算；PCr為海雜波回波功率，單位為dBm，可通過式(2)進(jìn)行計(jì)算。

在噪聲背景中，目標(biāo)信噪比可表述為

SNR=PTr-PN

(6)

式中：PTr為點(diǎn)目標(biāo)回波功率，單位為dBm，可通過式(2)進(jìn)行計(jì)算；PN為接收機(jī)噪聲背景功率，單位為dBm，可通過式(3)進(jìn)行計(jì)算。

鑒于海雜波分布的拖尾特性，在信號(hào)檢測(cè)過程中，通過設(shè)定適當(dāng)?shù)腃FAR門限，抑制短時(shí)相關(guān)(2 ms左右)雜波虛警，容忍“海尖峰”長(zhǎng)時(shí)相關(guān)(2 s左右)雜波虛警，再通過掃描間相關(guān)的處理方法對(duì)“海尖峰”類虛警進(jìn)行抑制，具體實(shí)施方法可參考相關(guān)文獻(xiàn)[5]。

2 海雜波單位面積反向散射率σ0建模

海雜波建模模型復(fù)雜，目前尚沒有完整的海洋電磁學(xué)理論進(jìn)行可靠支撐，已經(jīng)獲得較廣泛認(rèn)可的模型，都是基于大量實(shí)際回波數(shù)據(jù)的規(guī)律總結(jié)和參數(shù)提煉。本文計(jì)算采用了被廣泛采納認(rèn)可的三種海雜波單位面積反射率σ0模型：GIT模型、TSC模型和HYB模型。這里沒有選用SIT模型，因?yàn)樵撃Ｐ椭唤o出了迎風(fēng)和側(cè)風(fēng)向的海雜波σ0建模方法。

通過比較三者計(jì)算結(jié)果進(jìn)行相互驗(yàn)證，最終以三者的均值作為參與系統(tǒng)理論計(jì)算的海雜波單位面積反射率σ0參考值，以期獲得最貼近實(shí)際環(huán)境的結(jié)果，理論計(jì)算所使用的海情海況描述相關(guān)參數(shù)見表1。

表1 海情海況標(biāo)準(zhǔn)描述 m

2.1 GIT模型[6]

Georgia技術(shù)研究所建立了一個(gè)確定性、參數(shù)化的海雜波單位面積平均反向散射率模型，該模型包含擦地角、風(fēng)速、與風(fēng)速相關(guān)的平均浪高、雷達(dá)波束和風(fēng)向夾角關(guān)系、雷達(dá)工作頻率和收發(fā)極化形式在內(nèi)的一系列相關(guān)參數(shù)。

GIT模型包括海面后向散射機(jī)制在內(nèi)的多項(xiàng)理論模型為基礎(chǔ)，結(jié)合引入經(jīng)驗(yàn)因子而建立，是目前被廣泛所接受描述海雜波單位面積平均反向散射率最好的模型之一。GIT模型由以下幾個(gè)因素構(gòu)成：

(1) 干涉因子AI

σφ=(14.4λ+5.5)φhav/λ

(7)

(8)

式中：λ為雷達(dá)波長(zhǎng)；φ為波束擦地角；hav為平均浪高。

(2) 風(fēng)向因子Au

Au=exp[0.2cosΨ(1-2.8φ)(λ+0.02)-0.4]

(9)

式中：Ψ為波束指向與風(fēng)向間夾角。

(3) 風(fēng)速因子Aw

qw=1.1/(λ+0.02)0.4

(10)

Aw=[1.94Vw/(1+Vw/15.4)]qw

(11)

式中：Vw為海面風(fēng)速，若對(duì)于“充分發(fā)展”的海面，Vw和平均浪高h(yuǎn)av間存在下述關(guān)系

(12)

對(duì)于工作頻率在10 GHz以下的水平極化體制雷達(dá)，海雜波單位面積反射率σ0(單位：dBm2/m2)的GIT模型可描述為

σ0=10lg(3.9×10-6λφ0.4AIAuAw)

(13)

2.2 TSC模型[6]

TSC公司提供了一個(gè)可良好擬合Nathanson數(shù)據(jù)的模型。該模型與波束擦地角、Douglas海情海況等級(jí)、風(fēng)向角度、雷達(dá)工作頻率和收發(fā)極化形式相關(guān)。在TSC模型中，有效的風(fēng)速和浪高直接使用Douglas海情海況等級(jí)計(jì)算推導(dǎo)獲得。

在相對(duì)謹(jǐn)慎和保守的性能預(yù)測(cè)時(shí)，推薦使用TSC模型，因?yàn)橄啾菺IT模型而言，TSC模型更加接近平均條件(注：TSC模型針對(duì)小擦地角評(píng)估準(zhǔn)確性優(yōu)于中、大擦地角條件)。

TSC模型由以下幾個(gè)因素構(gòu)成：

(1) 小擦地角因子GA

σz=0.115S1.95

(14)

σα=14.9φ(σz+0.25)/λ

(15)

(16)

式中：λ為雷達(dá)波長(zhǎng)(單位：ft，1ft=0.304 8 m)；φ為擦地角(單位：rad)；S為海情海況等級(jí)。

(2) 風(fēng)速因子Gw

Vw=6.2S0.8

(17)

Q=φ0.6

(18)

A1=(1+(λ/0.03)3)0.1

(19)

A2=(1+(λ/0.1)3)0.1

(20)

A3=(1+(λ/0.3)3)Q/3

(21)

A4=1+0.35Q

(22)

A=2.63A1/(A2A3A4)

(23)

Cw=[(Vw+4.0)/15]A

(24)

式中：Vw為海面風(fēng)速(單位：節(jié))。

(3) 方向因子Gu

B1=exp(-φ/0.17)

(25)

B2=(λ2+0.005)0.2

(26)

Gu=exp(0.3cosΨB1/B2)

(27)

式中：Ψ為波束指向與風(fēng)向間夾角(單位：rad)。

對(duì)于水平極化體制雷達(dá)，海雜波單位面積反射率σ0(單位：m2/m2)的TSC模型可描述為

σ0=1.7×10-5φ0.5GuGwGA/(λ+0.05)1.8

(28)

2.3 HYB模型[6]

使用HYB模型估算海雜波平均單位面積方向散射率時(shí)，以5級(jí)海情、0.1°擦地角、垂直極化和迎風(fēng)觀測(cè)為參考基準(zhǔn)值，對(duì)于海況等級(jí)、波束擦地角、極化方式和波束與風(fēng)向夾角等變量，提供了約4 dB的動(dòng)態(tài)調(diào)節(jié)范圍。過渡擦地角取值源于GIT模型和Barton的分析結(jié)果。

HYB模型由以下幾個(gè)因素構(gòu)成：

(1) 參考反向散射率σ0(ref)

對(duì)于X波段的雷達(dá)，參考反向散射率σ0(ref)可描述為

σ0(ref)=24.4lg(f)-65.2

(29)

式中：f為雷達(dá)工作頻率。

(2) 擦地角調(diào)節(jié)因子Kg

對(duì)于擦地角調(diào)節(jié)因子，定義參考擦地角φr(單位：°)和過渡擦地角φt(單位：°)

φr=0.1°

(30)

φt=arcsin(0.063 2λ/σh)

(31)

式中：σh為浪高均方根值，其與海情海況等級(jí)S間存在關(guān)系如下

σh=0.031S2

(32)

對(duì)于φr≤φt≤ 30°的情況，擦地角調(diào)節(jié)因子Kg可定義為

(33)

對(duì)于φr>φt的情況，擦地角調(diào)節(jié)因子Kg可定義為

(34)

(3) 海情海況調(diào)節(jié)因子Ks

Ks=5(S-5)

(35)

(4) 極化形式調(diào)節(jié)因子Kp

Kp=1.1ln(hav+0.015)-1.1lnλ-

1.3ln(φ/57.3+0.000 1)-9.7

(36)

式中：hav為平均浪高；λ為雷達(dá)工作波長(zhǎng)。

海情海況等級(jí)S與平均浪高h(yuǎn)av間存在如下關(guān)系

hav=0.08S2

(37)

(5) 風(fēng)向調(diào)節(jié)因子Kd

Kd=[2+1.7lg(0.1/λ)](cosΨ-1)

(38)

式中：Ψ為相對(duì)于雷達(dá)波束照射方向的風(fēng)向，將迎風(fēng)照射時(shí)的風(fēng)向角度定義為0°。

對(duì)于水平極化體制X波段雷達(dá)，海雜波單位面積反射率σ0的HYB模型可描述為

σ0=σ0(ref)+Kg+Ks+Kp+Kd

(39)

2.4 海雜波單位面積反射率σ0模型共性分析

通過對(duì)GIT模型、TSC模型和HYB模型進(jìn)行仿真計(jì)算可以看出，本類經(jīng)驗(yàn)性模型體現(xiàn)的共同特征包括：

(1) 海雜波平均單位面積反向散射率隨海情海況等級(jí)、工作頻率和擦地角的增加而增加；

(2) 在小擦地角條件下，海雜波平均單位面積反向散射率與擦地角的依賴關(guān)系更加明顯；

(3) 海雜波平均單位面積反射率在迎風(fēng)向最大、順風(fēng)向最小，中間角度基本平滑過渡。

3 載機(jī)飛行高度對(duì)海面目標(biāo)探測(cè)影響性分析

通過上述分析可知，海雜波單位面積反向散射截面積與擦地角間存在明顯的強(qiáng)相關(guān)關(guān)系。因此，載機(jī)飛行高度將直接影響海雜波存在的延伸距離，也同步影響著雷達(dá)的海面目標(biāo)檢測(cè)性能。

3.1 載機(jī)飛行高度與雜波間相關(guān)趨勢(shì)分析

按表2所示參數(shù)進(jìn)行仿真，其仿真結(jié)果如圖1和圖2所示。仿真結(jié)果使用極坐標(biāo)顯示，其中，半徑表示相對(duì)強(qiáng)度；方位表示風(fēng)向，以真北順時(shí)針為參考，圖中所繪曲線均經(jīng)過海情250 m條件下的海雜波最小單位面積反射率σ0值歸一化處理。

表2 仿真參數(shù)設(shè)置

圖1 東南風(fēng)1級(jí)海情不同高度σ0仿真結(jié)果

圖2 東南風(fēng)3級(jí)海情不同高度的仿真結(jié)果

通過觀察可以看出，在相同海情條件下，載機(jī)飛行高度越高，海雜波單位面積反射率越強(qiáng)；在相同高度條件下，海情等級(jí)約高，海雜波單位面積反射率越強(qiáng)；在水平極化天線設(shè)計(jì)中，迎風(fēng)面海雜波單位面積反射率強(qiáng)，背風(fēng)面海雜波單位面積反射率弱，其間角度均勻過渡。

圖3和圖4為試驗(yàn)雷達(dá)同一架次試飛的實(shí)際回波視頻。從截圖可以看出，在相同風(fēng)向風(fēng)速和海情海況條件下，隨載機(jī)高度提升，雜波在迎風(fēng)向上的延伸現(xiàn)象非常明顯：在迎風(fēng)向上，當(dāng)載機(jī)高度為250 m左右時(shí)，海雜波延伸距離約25 km；當(dāng)載機(jī)高度為1 000 m左右時(shí)，海雜波延伸距離擴(kuò)展至50 km以上。

圖3 東南風(fēng)3級(jí)海情，250 m高度實(shí)際回波

圖4 東南風(fēng)3級(jí)海情，1 000 m高度實(shí)際回波

3.2 載機(jī)飛行高度對(duì)探測(cè)影響性分析方法

載機(jī)飛行高度對(duì)海面監(jiān)視雷達(dá)海面目標(biāo)探測(cè)性能影響性分析的步驟為：

(1) 設(shè)定仿真場(chǎng)景，包括：海洋環(huán)境參數(shù)、雷達(dá)設(shè)計(jì)參數(shù)、海面目標(biāo)參數(shù)、載機(jī)飛行高度參數(shù)、目標(biāo)檢測(cè)參數(shù)；

(2) 海雜波單位面積反向散射率計(jì)算；

(3) 海面目標(biāo)信雜比(SNR)計(jì)算；

(4) 海面目標(biāo)信噪比(CNR)計(jì)算；

(5) 海雜波雜噪比(CSR)計(jì)算；

(6) 區(qū)分雜噪分界距離，當(dāng)CSR> 1.0時(shí)，設(shè)定為雜波區(qū)，使用高門限檢測(cè)策略，抑制虛警概率；當(dāng)CSR≤1.0時(shí)，設(shè)定為噪聲區(qū)，使用低門限檢測(cè)策略，提高檢測(cè)概率；

(7) 對(duì)不同海情海況、不同飛行高度、不同目標(biāo)距離的條件，進(jìn)行獨(dú)立仿真，并通過繪圖形式進(jìn)行對(duì)比，判定目標(biāo)檢測(cè)性能。

3.3 載機(jī)飛行高度對(duì)探測(cè)影響性仿真結(jié)果

按表3所示參數(shù)進(jìn)行設(shè)置，通過仿真可以看出，3級(jí)海情，載機(jī)飛行高度200 m條件下，保守估計(jì)的雜噪分界距離約35 km，在雜波區(qū)內(nèi)，目標(biāo)SCR> 7 dB，在噪聲區(qū)內(nèi)，目標(biāo)SNR> 4 dB，因此，可獲得不小于50 km的探測(cè)能力，如圖5所示；載機(jī)飛行高度400 m條件下，保守估計(jì)的雜噪分界距離約44 km，在雜波區(qū)內(nèi)，部分距離段目標(biāo)SCR< 7 dB，在噪聲區(qū)內(nèi)，目標(biāo)SCR> 4 dB，因此，雷達(dá)在雜波區(qū)15 km～44 km段會(huì)失去目標(biāo)的穩(wěn)定檢測(cè)能力，如圖6所示。

表3 仿真參數(shù)設(shè)置

圖5 3級(jí)海情200 m高迎風(fēng)探測(cè)10 m2目標(biāo)能力

圖6 3級(jí)海情400 m高迎風(fēng)探測(cè)10 m2目標(biāo)能力

4 結(jié)束語

本文利用三種典型海雜波單位面積反向散射率經(jīng)驗(yàn)?zāi)Ｐ?，完成海雜波簡(jiǎn)單建模。在此基礎(chǔ)上，設(shè)定特定應(yīng)用場(chǎng)景，提出一種評(píng)估機(jī)載雷達(dá)工作高度與海面目標(biāo)檢測(cè)能力影響性的定量分析評(píng)估方法。

通過使用該方法，可在不同海情海況條件下，結(jié)合作戰(zhàn)任務(wù)，提前靈活規(guī)劃航線高度，實(shí)現(xiàn)雷達(dá)對(duì)海探測(cè)性能最優(yōu)化，對(duì)增進(jìn)任務(wù)完成效果具有明確指導(dǎo)意義。

[1] WARD K D,BARKER C J,WATTS.Maritime surveillance radar.Part 1:Radar scattering from the ocean surface[J].IEEE Proceedings of Radar & Signal Processing,1990,137(2):51-62.

[2] WARD K D,BAKER C J,WATTS.Maritime surveillance radar.Part 2:Detection performance prediction in sea clutter[J].IEEE Proceedings of Radar & Signal Processing,1990,137(2):63-72.

[3] SKOLNIK M I.雷達(dá)手冊(cè)[M].3版.北京:電子出版社,2010.

SKOLNIK M I.Radar handbook[M].3rd ed.Beijing:Publishing House of Electronics Industry,2010.

[4] 弋 穩(wěn).雷達(dá)接收機(jī)技術(shù)[M].北京:電子工業(yè)出版社,2005.

YI Wen.Radar receiver technology[M].Beijing:Publishing House of Electronics Industry,2005.

[5] SMITH J M,LOGAN R H.AN/APS-116 periscope detection radar[J].IEEE Transactions on Aerospace and Electronic Systems,1980,16(1):66-73.

[6] ANTIPOV I.Simulation of sea clutter returns[R].DSTO-TR-0679,AR-010-560.Salisbury South Australia:Electronic and Surveillance Research Laboratory,1998.

梁 浩 男，1983年生，碩士研究生。研究方向?yàn)橹鄙龣C(jī)機(jī)載雷達(dá)電訊總體技術(shù)。

李國剛 男，1983年生，碩士研究生。研究方向?yàn)橹鄙龣C(jī)機(jī)載雷達(dá)電訊總體技術(shù)。

徐艷國 男，1980年生，碩士研究生。研究方向?yàn)闄C(jī)載預(yù)警雷達(dá)電訊總體技術(shù)。

Analysis of the Influence of Flight Height on the Sea Surface Surveillance Radar Performance

LIANG Hao，LI Guogang，XU Yanguo

(Nanjing Research Institute of Electronics Technology,Nanjing 210039,China)

The radar equation of dot type target and surface clutter ,and the calculation of receiver front terminal noise power is reviewed，and then three normally used empirical models of mean sea clutter reflectivity are introduced in detail.The characteristic of sea clutter reflectivity under different flight height and specifically sea status is analyzed by means of introduced empirical models.The trend consistency between the practical radar sea clutter echo and empirical models is compared and verified.The characteristic of the target detection in sea clutter is taken into account,then a fast method to analysis and estimation the flight height influence on the radar performance is given.The suggestion of flight height selection to maximize the air-born surveillance radar performance is given in the end.

sea surface surveillance radar; sea clutter; flight height; sea status

??工程·

10.16592/j.cnki.1004-7859.2016.10.003

梁浩 Email：108754226@qq.com

2016-07-18

2016-09-20

TP971.1

A

1004-7859(2016)10-0008-05