在机载雷达系统设计中,雷达的空域覆盖是表征雷达性能的一个重要指标。它能清晰、形象地描述雷达在不同空域的最大作用距离,协助优化雷达的系统设计。然而,准确的仿真机载雷达的空域覆盖范围是非常困难的,机载雷达的探测距离相对较远,电磁波在远距离传输过程中受影响的因素很多,例如大气密度随高度引起的折射率变化、不同类型的天气状况、不同类型的地表对电磁波的反射、散射等[1-3]。这些均会影响电磁波的传输从而影响雷达的空域覆盖,特别是地表面反射的影响。直达波和反射波的叠加使雷达的空域覆盖范围在空域上产生起伏,对雷达的探测影响很大。因此,研究机载雷达的多径效应,对雷达系统的设计和优化具有指导意义。
电磁波在均匀介质中沿直线传播的,然而实际地球的大气层并非均匀介质,通常可以把地球的大气层分为3个部分:对流层、平流层和电离层[4]。对流层的高度约为0~30 km,空气的密度随高度的升高而降低,因此,电磁波在传输的过程中会产生折射(向下弯曲)现象。在本文的仿真中,采用等效地球半径(4/3倍的地球半径)的方法研究对流层大气折射对电波传播的影响,平流层和电离层的影响暂不考虑。
电磁波在大气中传播,由于空气的密度随高度的升高而降低,所以实际的电波并非沿直线传播。为了满足电波的直线传播条件,通常采用等效地球半径的模型[4]。在建立空域覆盖范围的研究模型时,首先确定直达波,然后根据直达波计算反射波的幅值和相位,最后将这两种电波进行矢量叠加,确定电磁波能量的空间分布,进而确定雷达的空域覆盖范围。
图1给出了雷达空域覆盖范围的示意图,图中Rd为直达波,R1和R2为反射波,直达波和反射波在空间进行矢量叠加,决定了雷达的空域覆盖范围[5]。实际的雷达探测过程中,反射波不仅仅只有一条,但是经过多次反射形成的反射波能量相对较小,在研究多径效应影响时,通常只考虑反射最强的一条反射波。
图1 机载预警雷达多径对空域覆盖的影响
研究多径效应对雷达空域覆盖范围的影响中,需要计算出反射波的幅值和相位。因此,仿真过程中两个关键参数有待确定:1)反射点的位置;2)地面的反射系数。由反射点的位置能计算出R1和R2的长度,进一步确定反射波的相位。由地面的反射系数可以计算出反射波的幅值。
如果雷达和目标的距离比较近,地面可以近似看成平面,利用镜像原理可以比较容易求出反射点的位置。但是,在机载条件下,雷达探测距离较远,此时地球表面曲率对反射的影响就凸显出来,如图1所示。镜像原理在此处依然适用,但是反射面为反射点的切平面。根据文献[2],散射点的位置由如下的三次方程求解得到:
(1)
式中,hr为雷达到地面的高度,ht为目标位置点到地面的高度,re为地球的等效半径,r为雷达和目标之间对应地球表面的弧长,r1为雷达和反射点之间对应地球表面的弧长。
散射点位置可以解出为
r1=r/2-psin(ξ/3)
(2)
式中,
根据弧长r1和r2确定对应地心角Φ1和Φ2。由图中的三角关系可以计算出反射波的路径长度R1,R2和直达波路径长度Rd:
(3)
(4)
(5)
因此,两条路径的路径差为
ΔR=R1+R2-Rd
(6)
此时,两条路径之间的相位差为
(7)
散射波幅值的大小由反射系数确定。机载条件下地球表面的散射不仅需要考虑镜向反射,同时还需要考虑粗糙表面的散射和地球曲率引起的扩散,因此总散射系数可以表示为
Γt=Γ(H,V)×Sr×D
(8)
式中:Γ(H,V)是平面的镜向反射系数,反射系数和电磁波的极化有关,H为水平极化,V为垂直极化;D是扩散因子,由于此时的反射面是凸起的球面,根据射线光学[6],电磁波将被发射,发散的程度跟地球的半径和入射的擦地角等因素相关;Sr是粗糙面对平面反射的影响系数,由地面的随机起伏导致的漫反射,跟电磁波的波长和地面的起伏强度等因素相关,为了计算出地面的散射系数,需要分别计算这三个部分。
平面的镜向反射系数如下所示,式(9)为水平极化的镜向反射系数,式(10)为垂直极化的镜向反射系数[6]。
(9)
(10)
式中,Ψg为擦地角,ε为地面的介电常数。
实际的地球表面是随机粗糙面,如图2所示。图中给出了平面和不同程度粗糙面反射、散射差异的示意图。反射面的粗糙程度越大,散射越强。粗糙面对反射系数的影响用参数Sr表示[7]:
(11)
式中:σ为粗糙面的均方根高度,表明地球表面的起伏程度;Ψg为擦地角;λ为入射电磁波的波长,波长越短,地面散射越大。
图2 镜面反射和粗糙面反射差异示意图
由于地球表面是球形的,电磁波在发生反射的同时也会由于曲率而出现发散的现象,如图3所示。发散因子D的表达式如下所示[7]:
(12)
图3 电波传播发射示意图
地球表面的总反射系数Γt由上述三个参数确定后,可以计算出反射波的幅值。由直达波和反射波的矢量叠加,可以计算出多径效应对雷达空域覆盖范围的影响。通常采用方向图传播因子F表示多径效应的影响。F的表达式为
(13)
式中,Ed为直达波的电场,Ei为反射波的电场。
雷达探测距离受反射波的影响可以修正为[7]
(14)
图4~图6仿真了发射频率为500 MHz时不同高度、不同极化的雷达空域覆盖。仿真时设置地面的海拔高度均为1 000 m,机载雷达的海拔高度分别为4 000 m、5 000 m和6 000 m。天线的俯仰向波束宽度为10°,地球表面的起伏均方根高度为1 m。考虑多径效应时,两种极化的空域覆盖结果如图所示。其中,HH极化为电磁波的电场方向平行于水平面,VV极化为电磁波的电场方向垂直于水平面。
(a) HH极化
(b) VV极化
图4 机载雷达海拔高度4 000 m
(a) HH极化
(b) VV极化
图5 机载雷达海拔高度5 000 m
(a) HH极化
(b) VV极化
图6 机载雷达海拔高度6 000 m
由仿真结果可以看出,受多径效应的影响,雷达的空域覆盖在不同高度范围呈现增强和减弱。这种效应与直达波和反射波的相位差有关,当相位差为半波长的偶数倍时,干涉相长,雷达的探测距离增大;当相位差为半波长的奇数倍时,干涉相消,雷达的探测距离减小。当相位差为其他值时,雷达的探测距离介于最大和最小值之间[8-9]。所以,雷达的空域覆盖范围在空间分布中呈现强、弱起伏分布。
比较不同极化雷达的空域覆盖时,发现HH极化雷达的最大探测距离大于VV极化雷达,主要原因是雷达的最大探测距离由直达波和反射波相干叠加所决定。根据式(9)和式(10),HH极化的镜向反射系数大于VV极化的镜像反射系数,所以,叠加后HH极化电波的总能量大于VV极化。另一方面,VV极化的多径效应相对于HH极化小,测高性能较好。
当机载雷达处于不同海拔高度时,雷达空域覆盖的高度随海拔高度的升高而减小。原因是随着雷达高度的升高,直达波和反射波叠加时,反射波的传输距离增大,传输路径损耗也随之增加,探测距离将有所减小。多径效应的影响另一个表现是空域覆盖的空间起伏上,雷达的海拔高度越高,空域覆盖的起伏特性越明显(多径效应的影响越大)。主要原因是当雷达的海拔高度增大,波束照射到地面的面积也增大,将会有更多的反射波与直达波矢量叠加影响雷达的空域覆盖。
(a) HH极化
(b) VV极化
图7 机载雷达发射频率500 MHz
(a) HH极化
(b) VV极化
图8 机载雷达发射频率1.2 GHz
(a) HH极化
(b) VV极化
图9 机载雷达发射频率3.0 GHz
图7~图9仿真了雷达发射频率分别为500 MHz(P波段)、1.2 GHz(L波段)和3 GHz(S波段)时雷达HH极化和VV极化的空域覆盖范围。仿真时,假定雷达天线阵面的物理口径一致,天线阵子按照半波长的间隔排布。P波段的波束宽度为10°,根据理论公式折算出L波段和S波段的天线阵面波束宽带分别为4.2°和1.7°。地面的海拔高度均为1 000 m,起伏的均方根高度为1 m。机载雷达的海拔高度均为6 000 m。
由上述仿真结果,可以发现随着发射频率的增大,雷达多径效应对空域覆盖范围的影响逐渐减小。当发射频率为3 GHz时,多径对雷达空域覆盖的影响已经很小,如图所示,主要是两方面的原因:
1) 雷达的发射频率增大后,首先波束宽度变窄,被反射的波束范围也相对较小,形成的反射波较少。所以,多径效应对空域覆盖范围的影响随之减少。
2) 当雷达的发射频率增大(波长变小)时,地球的粗糙面散射增强,镜面反射减小(公式(11))。多径效应中镜面反射是影响雷达空域覆盖的主要因素,所以,多径效应对空域覆盖范围的影响也随之减少。
同时,由仿真结果可以看出随着雷达发射频率的增大,雷达的最大探测距离也随之增加(对于某RCS确定的目标,如20 m2目标),主要是天线波束变窄,天线的增益变大。
进一步,对地基雷达的空域覆盖范围进行仿真并与机载雷达进行比较。仿真时地面的海拔高度设置为1 000 m,雷达一般会放置在地势较高的阵地,假设雷达与地面的高度为50 m,雷达仰角设置为2°,波束的3 dB宽带为1.7°(与上图机载雷达发射频率3.0 GHz波形参数相当),仿真结果如图10所示。
(a) HH极化
(b) VV极化
图10 地基雷达发射频率3.0 GHz
根据仿真结果可知,雷达的空域覆盖范围与机载雷达的变化规律基本一致,HH极化的多径效应要强于VV极化;与机载雷达有所不同的是,多径效应对空域覆盖范围的影响主要针对空域的下半部分,如图10所示。主要是因为地基雷达离地面的距离很近,反射波的擦地角将会很小,所以,多径反射波中的大部分只能影响空域的下半部分。
(a) HH极化
(b) VV极化
图11 发射频率500 MHz海面起伏均方根高度0 m
(a) HH极化
(b) VV极化
图12 发射频率500 MHz海面起伏均方根高度2 m
(a) HH极化
(b) VV极化
图13 发射频率500 MHz海面起伏均方根高度5 m
图11~图19仿真了机载雷达在海面时空域覆盖范围随雷达发射频率和海面高度起伏的变化规律。图11~图13仿真了雷达发射频率为500 MHz,海面起伏均方根高度分别为0 m、2 m和5 m时的空域覆盖范围。与地面相比,海面的空域覆盖范围受电磁波极化影响更大,HH极化时,上半空域雷达的最大作用距离明显大于VV极化,主要原因是海面的介电常数大于地面,同时存在较大的虚部。介电常数模值影响反射系数的模值,虚部影响反射相位。同时,VV极化的反射系数受介电常数的影响大于HH极化(式(9)和式(10)所示)。所以两种极化的空域覆盖范围具有差异并且海面的差异大于地面[10]。
海面起伏均方根高度对应不同的海浪起伏,均方根高度越大,表示海面的海浪也越大,此时海面对电磁波的散射增强,平面的镜向反射减弱[11]。如图13所示,在空域覆盖的上半部分,当海面的均方根高度为5 m时空域覆盖范围没有起伏,说明此时雷达反射波中大部分被散射,与直达波叠加的反射波很小。而在雷达空域覆盖的下半部分,由于反射波的擦地角减小,散射特性减弱(式(11)),空域覆盖范围又出现了起伏特性。
(a) HH极化
(b) VV极化
图14 发射频率1.2 GHz海面起伏均方根高度0 m
(a) HH极化
(b) VV极化
图15 发射频率1.2 GHz海面起伏均方根高度2 m
(a) HH极化
(b) VV极化
图16 发射频率1.2 GHz海面起伏均方根高度5 m
图14~图16仿真了L波段的空域覆盖结果,此时雷达的波束宽度为4.2°。由于波束宽带较小,雷达照射到海面的范围也比较小,被反射的电磁波也较少,因此不同极化条件下空域覆盖范围的差异性比P波段的小。海面起伏均方根高度对空域覆盖范围的影响与P波段的规律基本一致。
(a) HH极化
(b) VV极化
图17 发射频率3.0 GHz海面起伏均方根高度0 m
(a) HH极化
(b) VV极化
图18 发射频率3.0 GHz海面起伏均方根高度2 m
(a) HH极化
(b) VV极化
图19 发射频率3.0 GHz海面起伏均方根高度5 m
图17~图19仿真了S波段机载雷达在海面上的空域覆盖范围,雷达的波束宽度为1.7°。仿真结果与雷达在地面上空的一致,几乎没有多径效应的影响,并且海面起伏均方根高度越大,受多径效应的影响越小。
上文仿真、分析了不同发射频率、不同高度、不同极化和不同地面类型雷达空域范围受多径效应的影响。P波段空域覆盖的起伏特性最明显,L波段次之,S波段受影响很小。在设计雷达系统时,考虑到多径效应的影响可以采用两种或者多种波形相互配合,增强雷达探测性能。利用一个波形相长干涉的特性增大了雷达的探测距离,同时采用另一个波形对相消干涉的空域部分进行空域补偿。
针对雷达远距离探测的威力和测高精度的问题,由于VV极化多径效应较HH极化弱,可以根据极化间的差异特性,采用两种极化协同的方法,采用HH极化雷达进行远距离探测、采用VV极化雷达进行高精度测高,两部雷达相互协同,增强雷达的探测性能。
本文基于电波传播理论仿真了雷达空域覆盖范围受多径效应的影响,研究了雷达处于不同海拔高度、不同发射频率和不同地表类型雷达空域覆盖范围的变化规律。仿真结果和理论分析发现:HH的反射系数大于VV极化,所以HH极化时雷达最大探测距离大于VV极化。机载雷达海拔高度越大,路径损耗越大,雷达最大探测距离随海拔高度的增大而减小。机载雷达所处的海拔高度越高,雷达波束的地面投影越大,多径效应引起的空域覆盖范围的起伏特性越显著。发射频率越高,波束越窄,地面镜向反射系数越小(散射系数越大),地面反射波对直达波的影响越小。海面的介电常数远大于地面,所以不同极化之间的差异较大。海面的起伏均方根高度越大,镜向反射系数越小(散射系数越大),反射波对直达波的影响越小。最后,分析了不同发射频率受多径效应影响对空域覆盖范围的影响,采用低频和高频相互配合的方法,增强雷达的探测性能。同时,针对不同极化雷达的差异性较大的特点,分析了采用极化间协同,提高雷达探测威力的同时增强雷达的测高精度。
[1] 肖景明,王元坤.电波传播工程计算[M].西安:西安电子科技大学,1989.
[2] SKOLNIK M I.雷达手册 [M]. 2版. 北京:电子工业出版社,2003.
[3] 丁鹭飞,耿富录,陈建春.雷达原理[M].5版. 北京:电子工业出版社,2014.
[4] MAHAFZA B R.雷达系统分析与设计[M].3版. 北京:电子工业出版社,2016.
[5] 陆鹏程,江胜利,同伟.米波雷达应对多路径效应的系统设计[J].雷达科学与技术,2018,16(3):256-260.
LU Pengcheng, JIANG Shengli, TONG Wei. VHF Radar Designing Analysis Dealing with Multipath Effect[J]. Radar Science and Technology, 2018, 16(3):256-260. (in Chinese)
[6] 汪茂光.几何绕射理论[M].2版. 西安:西安电子科技大学,1994.
[7] MAHAFZA B R.Radar Systems Analysis and Design Using Matlab[M].Alabama:COLSA Corporation Huntsville, 2000.
[8] 汪安戈,胡国平,周豪. 雷达多径效应抑制技术分析及展望[J].火力与指挥控制,2019,44(5):12-16.
[9] 史鉴,钟文. 多路径效应对侦察测向的影响[J].舰船电子对抗,2019,42(3):26-29.
[10] 李海军,王天然. 海面环境对导弹雷达多径效应的影响及仿真分析[J].兵器装备工程学报,2021,42(1):85-89.
[11] 邓恩民,刘冬利,隋波,等. 复合多径效应下舰载雷达探测性能分析[J].指挥控制与仿真,2021,43(3):40-44.
刘 波 男,1974年5月出生,博士,高级工程师,中国电子学会高级会员,主要从事雷达系统论证和试验方面的研究工作。