0   引　言

对流层散射通信是一种利用对流层不均匀体对电波的散射作用来实现超视距通信（over-the-horizon communication, OHC）的方式. 由于这种通信方式具有传输容量大、通信距离远、保密性好、抗干扰和抗截获能力强等优点，自投入使用以来在军事和民用通信领域中始终占有重要的地位 ^[1-4].

对流层中经常存在着湍流运动，参与湍流运动的大气以一种随机的速度做不规则运动，形成一个个涡旋. 涡旋在入射电波的照射下变成偶极子，将电磁能量再辐射出去，即散射^[4]8. 散射的电波具有明显的衰落特性，并且快衰落包络大致服从瑞利分布，因此在大多数文献中通常将对流层散射信道作为瑞利衰落信道进行研究. 其中，对于宽带的对流层散射通信，信道被看作频率选择性衰落的瑞利信道，常用Kailath提出的抽头延迟线^[5]（tapped delay line,TDL）模型和Sunde提出的Sunde模型^[6]进行建模；而对于窄带的对流层散射通信，信道则被看作平坦衰落的瑞利信道，采用Clark和Jakes模型进行建模^[7].

但上述模型均为统计性模型，没有考虑具体信道的物理性质，而是根据实测数据的统计规律，从数学上作抽象的研究，因此难以根据实际环境因素的改变获取相应的信道传递函数. 此外，在非视距传输系统中，如果源自某一个散射体路径的信号功率特别强，信号的衰落会服从莱斯分布^[8]，而上述的统计性模型难以同时考虑信号的多种分布.

相较于统计性模型，确定性模型可基于环境特征进行精确的电磁建模，可靠地实现特定场景中的电波传播预测^[9]. 对流层散射的确定性模型方法主要有射线追踪（ray tracing，RT）法和抛物方程法（parabolic equation method, PEM）. RT能够描述电磁波的传播路径轨迹，便于计算多径时延. Dinc等^[10]提出了基于RT的多输入多输出（multiple-input multiple-output，MIMO）对流层散射信道模型，并计算出信道的功率时延谱. S. Zhang等 ^[11]考虑了射线的弯曲、折射效应，提出用三维RT方法来估计散射信道的功率损耗和时延. 但这些方法均不能同时表征复杂的气象和地理环境，并且射线路径的考虑过于理想化，且考虑的传输路径有限.

PEM是适用于近轴传播的全波算法，能够对辐射方向图进行建模，并准确地计算复杂大气和地理环境对电波传播的影响，体现电波的折射、反射、绕射、散射和衰减效应^[12]. 结合分步傅里叶变换（split-step Fourier transform, SSFT），能快速求解大区域复杂环境中的电波传播问题，成为了目前对流层电波传播模拟的主流算法. 为了描述对流层不均匀体对电波的散射作用，Hitney ^[13]将折射率随机扰动因子引入抛物方程，提出了对流层散射PEM. 李磊^[14]提出了一种计算散射抛物方程最大传播角和最大高度的方法，对散射PEM进行了改进，使计算同时具备高效性和准确性. 上述有关散射PEM的研究大多停留在预测传输损耗方面，鲜有研究散射信号特性，并将其用于通信性能评估方面的报道.

为此，本文采用对流层散射PEM，构建了考虑环境因素的对流层散射通信信道模型. 并进一步推导了在QPSK调制下，散射信道的误码率公式. 为验证散射PEM描述对流层散射信道的正确性，利用蒙特卡洛模拟方法，统计计算出对流层散射信道传输函数的概率密度分布. 最后，分析了不同通信距离下，不同气象气候条件对对流层散射信道误码率的影响.

1   湍流散射对电波传播的影响

1.1   湍流大气结构模型

现有理论中，湍流非相干散射被认为是造成对流层散射传播的主要原因. 该理论认为，对流层散射传播源于对流层中的湍流运动. 由于对流层中大气温度、湿度和气压等环境因素的改变，导致大气流动速度随机变化，形成一个个涡旋. 每一个涡旋都是一个介电常数局部不均匀体，因此在湍流运动中，折射率是随机起伏的，通常用大气折射率结构常数$C_{{n}}^2$来描述湍流起伏的强弱.

根据局地均匀各向同性的平稳湍流理论，得到大气折射率结构常数的计算公式为

式中：$a^{2}$为常数，在实验室测得其值约为2.8^[15]；$\alpha^{'}$为湍流扩散系数的比率，该值变化不大，可近似取为1；${L_0}$为湍流外尺度，根据实验资料，其量级约为几十米^[4]，在本文仿真中取10 m^[16]；大气修正折射率

$M{'^{}}$为折射率垂直梯度，

式中：$P$为大气压强，单位hPa；T为绝对温度，单位K；$e$为水汽压；$\Theta$为位温，单位K；$h$为海拔高度，单位m；$q$为比湿. 各参数表达式如下：

式中，$U$为相对湿度. 在实际情况中，除了大气压强基本只随高度递减外，对流层大气的温度、湿度均随时间、地点和高度的变化而变化. 但假定大气是静止和干净的情况下，对流层气象参数随高度呈现一定的规律性，能粗略反映某一特定地区多年的大气年平均情况. 下面给出某一地区温度$T$和湿度$U$关于高度$h$的变化关系 ^[17-18]：

式中：${T_0}$，${U_0}$分别为近地面的绝对温度和相对湿度. 对流层大气压强随高度的变化关系可近似表示为^[14]22

式中，${P_0}$为近地面的压强.

为了便于后续的分析，本文根据给定的近地面温度、湿度、压强，采用式(7)~(9)这三个经验公式来推导其垂直剖面数据. 若是有真实气象垂直剖面数据，也同样可用本文模型进行分析. 将气象数据带入式(1)~(6)可得大气折射率结构常数$C_{{n}}^2$. 再根据式(10)得出由湍流运动造成的折射率的随机扰动值^[13]：

式中：$r$为$[0,1]$之间均匀分布的随机数；${S _{n}}(k)$为折射率功率谱；$k = 2{\text{π}}/{L_{\text{s}}}$，${L_{\text{s}}}$为对散射场有贡献的不均匀性尺度^[4]12.

1.2   对流层散射PEM

本文采用PEM对传播环境进行电磁建模，设时谐因子为${{\text{e}}^{ - {\text{i}}\omega t}}$，假设波函数$\psi$与y无关，$\psi (x,{\textit{z}})$可用二维Helmholtz方程描述：

对于沿$x$轴正向传播的电磁波，简化场$u$可以表示为

将式(13)带入二维Helmholtz方程可推导出Feit-Fleck型宽角抛物方程：

式中：$u\left( {x,{\textit{z}}} \right)$为$(x,{\textit{z}})$处的简化场；${k_0}$为自由空间波数；$n$为媒质的折射指数. 引入SSFT技术，可得相应的递推公式^[19]：

式中：$\mathfrak{J}$和${\mathfrak{J}^{ - 1}}$分别为傅里叶变换和逆变换；${{\text{e}}^{{\text{i}}{k_0}\Delta x(n - 2)}}$为折射指数项，反映了空间媒质对电磁波的影响；${{\text{e}}^{{\text{i}}\Delta x\sqrt {k_0^2 - {p^2}} }}$为绕射项，反映了传播路径上障碍物对电波的绕射效应，其中$p = {k_0}\sin\; \theta$为角谱域变量，$\theta$为电磁波与水平方向上的夹角. 当给定初始场分布$u\left( {{x_0},{\textit{z}}} \right)$后，即可运用式(15)得到下一步进处的场分布$u\left( {{x_0} + \Delta x,{\textit{z}}} \right)$，从而迭代求解出整个计算空间中的场.

在处理湍流运动对电波传播造成的影响时，根据实际测得的近地面温度${T_0}$、近地面的相对湿度值${U_0}$和近地面压强${P_0}$，从而得到实时的大气折射率结构常数$C_{{n}}^2$，再利用式(10)和式(11)求得湍流运动导致的随机折射率扰动${n_{\text{f}}}$. 将湍流运动产生的折射率随机扰动因子${n_{\text{f}}}$和大气修正折射率剖面${M}$相结合，即可形成反应湍流运动影响的折射指数剖面，即

将${n_{\text{c}}}$带入式(15)即可分析湍流散射对电波传播的影响.

2   对流层散射通信

2.1   对流层散射通信信道传输函数

PEM的计算过程相当于求解单频点信号经过线性时不变系统的传输函数^[20-21]. 因此，可以用对流层散射PEM直接将对流层散射信道的传输函数$H$计算出来.

在PEM中，收发天线功率之比可用§1.2中求得的简化场$u(x,{\textit{z}})$来表示^[19]，表达式为

式中：${P_{\text{r}}}$为接收功率；${P_{\text{t}}}$为视轴等效全向辐射功率；$\lambda$为波长. 式(17)并没有考虑到天线增益. 这里，我们将散射信道传输函数$H(\omega ,x,{\textit{z}})$模的平方定义为考虑天线增益后收发天线功率之比，即

式中：${G_{\text{t}}}$和${G_{\text{r}}}$分别为收发天线增益. 对于散射信道传输函数$H(\omega ,x,{\textit{z}})$，有

2.2   对流层散射通信平均误码率

为实现高速率、高可靠性的对流层散射通信，需要对系统的通信性能进行分析. 误码率是研究系统通信性能的一个重要指标，研究不同条件下系统误码率性能的变化，有助于实现可靠、高质量的散射通信. 目前以AN/TRC-170为代表的很多散射通信站均采用BPSK/QPSK调制，对于QPSK调制，在只考虑高斯白噪声影响的情况下，系统的误比特率可表示为^[1]

式中：${E_{\text{b}}}$为平均比特能量；$\dfrac{{{N_0}}}{2}$为噪声功率谱密度；$Q(x)$函数为标准正态分布的右尾函数，

经对流层散射信道$H$后，有效比特信噪比为$\dfrac{{{{\left| H \right|}^2}{E_{\text{b}}}}}{{{N_0}}}$，因此误码率可以表示为^[22]

令频谱效率$\dfrac{{{R_{\text{b}}}}}{W} = 2$，其中${R_{\text{b}}}$为比特率，$W$为带宽，可将式(22)改写为

式中，N为噪声功率. QPSK系统对应的在对流层散射信道下的平均误码率为

式中，${f_{{H}}}(H)$为信道传输函数的概率密度分布，可运用蒙特卡洛模拟方法，通过对流层散射PEM多次计算统计得到.

3   仿真与分析

3.1   模型验证

为验证模型描述对流层散射信道传输损耗$({P_{\text{r}}}/{P_{\text{t}}})$的结果，比较PEM与ITU-R P.617-5模型^[23]计算结果. ITU模型中采用的是平均折射率与折射率梯度，为便于比较，本次仿真中PEM与ITU-R P.617-5模型均采用标准大气条件，即${N_0} = 315$ N单位，${\text{d}}N/ {\text{d}}h = - 39$ N单位/km. 发射天线高度${h_{\text{t}}}$，接收天线高度${h_{\text{r}}}$，3 dB波瓣宽度${W_{\rm{b}}}$，天线仰角，频率$f$如表1所示. 图1给出了不同距离下PEM计算的中值结果和ITU-R P.617-5模型不同百分比时间内未超过的损耗结果.

表  1  天线参数设置

Tab.  1  Antenna parameter setting

ht/m	hr/m	${W_{\rm{b} } }$/(°)	仰角/(°)	f/GHz	极化方式
20	20	1	0	3	水平极化

下载: 导出CSV 

| 显示表格

图  1  PEM与ITU-R模型路径损耗对比

Fig.  1  Comparison of path loss between PEM and ITU-R model

下载: 全尺寸图片 幻灯片

由于PEM计算的是损耗的瞬时分布，因此利用PEM对相同传播条件下的传输链路进行1000次模拟取中值，得到的结果应接近ITU模型中50%时间内不超过的平均年中值传输损耗，图1中的结果也确实与50%曲线最为吻合. 另外，图1中传播距离较近时，PEM模型与ITU模型出现较大差异. 这是因为此时无线电波主要以视距(line-of-sight, LoS)、绕射等传播方式为主，而本文所用的仅是ITU-R P.617-5模型中与对流层散射相关的部分，因此在较近距离时，ITU模型并不适用.

为了验证模型描述散射信道衰落特性的结果，采用文中提到的PEM进行多次计算，得到信道$H$的分布，再将其拆分成慢衰落信道${H_{\text{L}}}$与快衰落信道${H_{\text{S}}}$的结合. 即将信号看成是先通过${H_{\text{L}}}$经历固定的衰减和相移，再通过${H_{\text{S}}}$经历随机的衰落和相移.

由于慢衰落描述信号电平中值的长期变化，将慢衰落传输函数${H_{\text{L}}}$模的平方定义为对流层散射PEM计算所得的接收功率中值${P_{\text{m}}}$与视轴等效全向辐射功率${P_{\text{t}}}$之比，即

慢衰落传输函数${H_{\text{L}}}$的相位与其对应的接收场相同，快衰落的信道传输函数${H_{\text{S}}}$可由下式求出：

运用PEM求得的简化场$u(x,{\textit{z}})$是个包含幅度信息和相位信息的复数. 在水平极化情况下，只有${E_{{y}}}$为非零的电场分量，波函数$\psi (x,{\textit{z}}) = {E_{{y}}}(x,{\textit{z}})$，因此取其相位为接收场的相位，$\psi (x,{\textit{z}})$可以表示为^[24]

设近地面温度${T_0} = 300\;{\text{K}}$，近地面压强${P_0} = 1\;013.25\;{\text{hPa}}$，近地面相对湿度${U_0} = 80\text{%}$，天线参数设置如表1所示. 最大LoS距离可由公式${d_{{\text{LoS}}}} < 4.12\sqrt {{h_{\text{t}}}} + 4.12\sqrt {{h_{\text{r}}}}$计算，根据表1算得最大LoS距离约为37 km. 分别计算通信距离远离LoS距离即$d = 300$ km和接近LoS距离即$d = 60$ km处的快衰落信道传输函数概率密度分布，计算高度在5000 m以内. 图2给出了上述参数下${H_{\rm{S}}}$幅值和相位的概率密度函数(probability density function, PDF).

从图2可以看出，当通信距离较远时，收发天线之间没有直射波路径，接收信号为各个幅度、相位彼此无关的路径分量的叠加. 快衰落信道包络满足瑞利分布，相位在$[ - {\text{π, π}} ]$内均匀分布.

图  2  快衰落信道分布特性

Fig.  2  Distribution characteristics of fast fading channels

下载: 全尺寸图片 幻灯片

当通信距离变短到接近LoS距离时，源自某一散射体路径的信号功率会增强，接收端的信号可以看作功率较强的主信号与服从瑞利分布的多径信号分量的和. 其幅度服从莱斯分布，相位则趋向于分布在0附近.

综上所述，在路径损耗方面PEM结果与ITU模型较为吻合；在衰落特性方面，得到的快衰落包络和相位分布与众多理论研究和链路试验^[1-4,25]中描述的一致.

3.2   对流层散射超视距通信仿真

利用§1.1中的湍流模型，可以给出在高度2 km处大气折射率结构常数$C_{{n}}^2$随绝对温度、相对湿度等气象气候条件的变化，结果如图3所示. 文献[26]中提到，气压对$C_{{n}}^2$的影响很小，可以忽略. 因此仿真没有考虑气压变化对$C_{{n}}^2$的影响，近地面压强取值与§3.1中一致.

图  3  温度、湿度对$C_{{n}}^2$的影响

Fig.  3  Influence of temperature and humidity on $C_{{n}}^2$

下载: 全尺寸图片 幻灯片

由图3可以看到，随着近地面温度和近地面相对湿度的增加，在2 km高空中的湍流活动也逐渐变强.

对流层散射信号的慢衰落产生于气象气候条件的变化^[4]56，气象气候条件决定大气折射率结构常数$C_{{n}}^2$，而文献[16]指出大气折射率结构常数与对流层传输损耗之间有着强的相关性，因此气象气候条件也影响着对流层传输损耗. 利用对流层散射PEM计算出成都7月平均气候条件和12月平均气候条件下的传输损耗，结果如图4所示，天线参数设置同表1.

图  4  成都冬季和夏季的路径损耗对比

Fig.  4  Comparison of path loss between winter and summer in Chengdu

下载: 全尺寸图片 幻灯片

从图4可以看出，7月份路径损耗比12月小，这与文献[1]中所描述的我国散射信号电平夏季最高、冬季最低是一致的. 文献[14]中提到，对流层散射传播慢衰落的季节变化可能是由于折射率结构常数变化引起的，而折射率结构常数因地域、季节而异. 在本文的模型中，其变化取决于温度和湿度，因此温度和湿度较高的7月份的路径损耗较小. 由此可见对流层散射PEM能考虑气象环境的变化对对流层散射传播的影响.

3.3   误码率计算

通过调节发射天线功率，并将接收机处的信噪比转化为比特信噪比，得出了不同比特信噪比下PEM仿真误码率与理论误码率的对比结果. 如图5所示，天线参数设置同表1.

图  5  误码率仿真结果与理论结果对比

Fig.  5  Comparison of simulation results and theoretical results

下载: 全尺寸图片 幻灯片

d=300 km时，信号包络服从瑞利分布，PEM仿真结果与理论瑞利误码率曲线吻合较好. 此时误码率对信噪比的变化不敏感，通过增加信噪比改善误码率的效果不显著，应该采用分集技术等抗衰落技术来改善误码率. d=60 km时，信号包络服从莱斯分布，PEM仿真结果与理论莱斯误码率曲线吻合较好. 此时误码率随着信噪比提高快速下降，通过增加发射天线功率，提高信噪比，能有效改善通信质量. 在莱斯分布中，莱斯因子可由PEM计算出的快衰落信道传输函数的概率密度分布拟合得到，环境因素和距离的变化会对应得到不同的莱斯因子. 这说明了PEM能根据环境因素和通信距离的变化综合考虑对流层散射信道中信号的多种分布，从而计算出相应的误码率.

按照前文所述方法，对对流层散射信道的误码率随气象气候条件的变化进行分析，相关系统参数设置参考了文献[2]中的工程示例，如表2所示，其余参数设置同表1. 图6给出了不同通信距离和不同湿度环境下误码率随近地面温度的变化.

表  2  系统参数设置

Tab.  2  System parameter setting

d/km	Pt/W	N0/( dBm/Hz)	W/MHz	Gt/dB	Gr/dB
300	1000	−174	1	25	25
60	100	−174	1	5	5

下载: 导出CSV 

| 显示表格

图  6  不同通信距离下温度、湿度对误码率的影响

Fig.  6  Relationship between BER and climatic conditions at different communication distances

下载: 全尺寸图片 幻灯片

d=300 km时，信道服从瑞利分布，误码率随温度和湿度的变化较为平缓. 此时温度和湿度的改变对散射通信性能的影响较小，这表明温度和湿度的升高没有影响信道状态的概率密度分布，仅仅只是减小了电波的传输损耗，提高了接收端的信噪比. d=60 km时，信道服从莱斯分布，误码率随温度和湿度的变化较为剧烈，且变化的速度在逐渐加快. 此时温度和湿度的改变对散射通信的影响较大，这表明温度和湿度的变化不仅会影响信噪比，还会影响折射率梯度和湍流强度，从而使莱斯因子发生改变.

4   结　论

目前对于散射信道的研究大多基于瑞利分布的假设构建信道模型，没有考虑莱斯分布的情况及环境因素对信道的影响，在建模完备性上有一定局限. 对此，本文构建了对流层散射信道模型，推导了QPSK调制下散射信道误码率表达式，并对实际传播场景下的散射通信链路进行电磁建模和仿真计算，分析了不同气象环境因素对误码率的影响. 结果表明，较近通信距离下误码率随温度、湿度变化较为剧烈；通信距离较远时，环境对误码率的影响较小；相同通信距离下，高温和潮湿气候有利于超视距通信.