《基于抛物方程的大气波导环境下电波传播的研究》.doc

基于抛物方程的大气波导环境下电波传播的研究 1 概述 第二次世界大战以来，雷达已在许多领域得到广泛的应用。人们在利用雷达进行探测时，经常会发现一些电磁波的异常传播现象。其中一种显著的现象是在一定的气象条件下，在大气边界层尤其是在近地层中传播的电磁波，受大气折射的影响，其传播轨迹弯向地面，当曲率超过地球表面曲率时，电磁波会部分地被陷获在一定厚度的大气薄层内，就像电磁波在金属波导管中传播一样，这种现象称为电磁波的大气波导传播，形成波导传播的大气薄层称为大气波导层。大气波导现象使得雷达有可能观测到数倍于雷达正常探测距离处的目标，实现所谓的超视距探测。 目前，国外在电波传播的抛物方程方法上做了许多研究。抛物方程（PE）方法最早是由Fock在1946年提出，但是直到1973年，随着计算机技术运用，Hardin和Tappert进一步提出了分步傅立叶方法,抛物方程方法才得以被运用到实际的工程计算中。与其它研究电波传播的方法相比，如几何光学、物理光学和简正模分析法等，抛物方程方法的限制是最小的。这种方法忽略了后向散射场，而且精确计算集中在近轴方向上。 本文介绍了大气折射的基本类型及其存在条件，给出了电波传播的抛物方程模型，并将其与双射线模型进行对比验证其可行性，然后利用该模型计算了标准大气、100m波导及100-150m波导环境下的传播损耗，最后讨论了大气波导对电波传播及雷达探测的影响。 2 原理 大气波导现象经常发生在海洋大气环境中。针对电磁波在大气波导中的异常传播现象，可用抛物方程方法进行分析计算其传播损耗。 2.1 大气折射与大气波导 1、大气折射的基本类型及其存在条件 影响大气环境中的电磁波传播特性的主要大气因子是大气折射率。对频率在1100GHz范围内的电磁波，大气折射率n或大气折射指数N(N单位)可表示为大气温度T单位：K、大气压力P单位：hPa和水汽压e单位：hPa的函数f(P,T,e)其关系由下式给出： （2.1） （2.2） 当电磁波传播距离很短时可近似认为地球表面为平面但若电磁波传播距离较长时就必须考虑地球曲率的影响此时为了将地球表面处理成平面通常使用进行了地球曲率订正的大气修正折射率m和大气修正折射指数又称大气折射指数模数M（M单位更为方便其表达式如下 式中R0=6.371×106m为平均地球半径，Z单位m）为地表以上的高度则简化为： （2.5） 将式分别对高度Z求导可得 （2.6） （2.7） 当大气折射指数垂直梯度单位m-1）dN/dZ 0时电磁波的传播轨迹将背着地球而凸起弯曲此时的大气为负折射当dNdZ= 0时电磁波的传播轨迹不发生弯曲沿直线传播此时的大气为零折射也称无折射当dNdZ0时电磁波的传播轨迹将凹着弯向地球此时的大气为正折射正折射包括正常折射也称标准折射、超折射、临界折射、陷获折射等表给出了各种折射类型的存在条件 dN/dZ（单位m-1） dM/dZ（单位m-1） 负折射 0 0.157 零折射（无折射） 0 0.157 正折射 0 0.157 正常折射（标准折射） -0.077~0 0.080~0.157 超折射 -0.157~-0.077 0~0.080 临界折射 -0.157 0 陷获折射 -0.157 0 2、陷获折射与大气波导 当dMdZ0(即dNdZ-0.157 m-1)时大气呈现陷获折射条件此时在大气中传播的在一定频率范围内的电磁波将部分地被陷获在大气波导层内传播 由式可知大气折射指数垂直梯度dNdZ或大气修正折射指数垂直梯度dMdZ与大气温、压、湿垂直梯度相关由于大气波导的存在条件是dMdZ0，而式?P/?Z0，即第二?T/?Z0且?e/?Z0，或者当?T/?Z和?e/?Z两项的综合贡献为负值且数值较大时，才有可能满足条件dM/dZ0，产生大气波导。 2.2 抛物方程（PE）模型 1、抛物方程方法的推导 PE方法是一种忽略电波后向传播，从而将Helmholtz方程近似为抛物线方程的计算方法。假设在推导过程中时谐因子取为e-jwt，在直角坐标系中，二维空间波ψ( x , z)满足Helmholtz方程[6]，即： （2.8） 式中，k=2π/λ为真空中的波数，λ为波长，n为介质的折射率，x为沿表面的水平距离，z为距离表面的高度。沿x轴正向传播的波函数为，代入式（2.8），根据近轴条件，得到标准抛物方程： （2.9） 2、分步傅里叶变换（SSFT） 目前PE的数值解法主要有分步傅立叶变换方法SSFT）、隐式有限差分法IFD）和有限元方法FEM）。相对于后两种方法SSFT在计算大场景、远距离、低仰角的对流层电波传播问题上具有计算方便特别是当大气折射率随传播距离和高度变化时、误差小和计算效率高的优点因此成为计算对流层电