垂直大气湍流影响下的天气雷达回波信号优化检测方法

郝孝智，付亚平*，李 芳，李延辉

（1.山西省大气探测技术保障中心，山西 太原 030002；2.山西省气象灾害防御中心，山西 太原 030002）

0 引言

天气雷达是指一类专门用来监测、预警强对流天气的气象雷达，通常是脉冲雷达，主要工作机理是先在固定频段上发出一些保留时间较短的脉冲波，然后使用雨水、雪花等降水粒子散射脉冲波，最后利用接收散射来的回波脉冲确定气象目标的空间分布等［1］。气象雷达的波段范围一般为［1，10］cm，各波段衰减很小，能更好地进行暴雨等气象的监测。

在天气雷达的工作过程中，必须充分考虑到大气湍流的干扰效果［2］。大气湍流是一种大气移动方式，它可以明显增强大气中水汽、能量和物质之间的垂直、水平流动能力，而大气湍流的产生也可导致脉冲波在大气中的扩散。

关于天气雷达回波检测，沈妍琰等［3］通过贝叶斯分类器完成天气雷达回波信号检测，首先，通过统计得出先验概率；然后，分析回波特征获取似然函数；最后，利用贝叶斯分类器完成回波信号检测，该方法有效增强了天气雷达的数据质量，降低了因数据同化对回波造成的不利影响。文浩等［4］采用模糊逻辑对天气雷达回波进行检测，通过提取回波特征构建隶属函数，然后与设定阈值对比完成回波检测，该算法显著提升了回波检测的准确性。但上述方法存在当检测环境较为恶劣或者检测范围较大的情况下，天气雷达所接收的回波信号质量较差，极易被噪声所淹没，影响检测结果。基于此，本文提出垂直大气湍流影响下的天气雷达回波信号优化检测方法，提升天气雷达回波信号的质量。

1 天气雷达回波信号优化检测

1.1 天气雷达回波信号检测原理

天气雷达回波信号检测原理：通过发射机产生高频脉冲电磁波，经由定向天线发射检测脉冲、接收回波脉冲，脉冲传播过程中存在部分能量的损耗，通过降水粒子散射脉冲电磁波至按收机内，经按收机放大回波脉冲信号，最后在屏幕上显示气象目标物体的回波频率、强度等信息［5］。气象目标物体与天气雷达间的距离采用脉冲发射、接收的渡越时间测量，脉冲的传播速度为常数。

设置天气雷达测量气象目标物与雷达之间的距离时，其发射与完成回波接收的渡越时间为t，得到气象目标物与雷达之间的距离

式中：c为脉冲在大气中的运动速度。

在脉冲运动速度确定的情况下，所得测距的精度

式中：Δη，Δt分别为雷达距离信号、时域信号的分辨率。

1.2 垂直大气湍流对天气雷达回波信号传输特性的影响

垂直大气湍流对天气雷达回波信号的影响表现为天气雷达回波信号的大气折射率跟随时空的变化而任意变化［6］。而大气折射率的变化判断因素为温度。地表与大气中的热量均源自阳光照射，地表空气受热后开始上升并与高处的冷空气接触出现温度变化，形成垂直大气湍流。

为表现出天气雷达回波信号在大气内的传输特性，采用大气湍流折射率结构函数En表示垂直大气湍流特性，En定义为

式中：n(r1)，n(r2)分别为空间内r1，r2两点的折射率；r为二者的距离，即r=|r1-r2|；Li、Lo分别为垂直大气湍流的内、外尺度；为系综平均。

以全面描述垂直大气湍流性质和变化规律为目的，研究折射率结构函数En(r)的傅里叶变换结果——大气折射率功率频谱φn(e)（e为空间频率），常用的大气湍流谱模型Kolmogorov 谱、Von Karman谱分别描述为

式中：垂直大气湍流外尺度相应的空间频率通过eo描述，eo=2π/Lo；为描述垂直大气湍流起伏强度的大气折射率结构常数。

描述垂直大气湍流强度的另一指标为大气相干长度R，其是回波信号在垂直大气湍流内传播至距离发射机z位置时，脉冲横截面的外差距离。R与垂直大气湍流强度呈负相关关系，R越大则湍流强度越弱，回波信号畸变越少［7］。基于Kolmogorov 湍流原理，可得出R与的空间分布关系为

式中：N为回波信号数；L为回波信号传播距离。

综上所述，当回波信号在垂直大气湍流内传播时，垂直大气湍流内折射率的起伏导致回波信号产生随机波动，并对应引发相关垂直大气湍流效应。

1.3 基于匹配滤波的天气雷达回波信号优化

受垂直大气湍流影响，回波信号质量下降，直接导致天气雷达回波信号检测精度的降低。为此提出通过分析滤波优化回波信号，匹配滤波器（WF）为接收待检测的回波信号时输出幅值最大、未接收时输出幅值最小的滤波器［8］，其输出为天气雷达接收的回波信号所具有的信噪比最大值。

将回波信号f(N)输入WF，与f(N)相随的噪声是高斯白噪声。此时按照最优接收原则可以得到f(N)波形相应的WF为

式中：υ为常数；τ为延迟时间；f*为伴有噪声的回波信号。

其离散函数表示为

式中：τ隶属于f(N)的自相关函数ρ(k-K)；k为WF长度。若k=K，则此时(k-K)存在最大值。当f(N)为实序列时，此时WF的冲激响应为

由式（9）可知，WF 的冲激响应为f(N)的镜像函数，得出WF转移函数εjwT的频率特性为

式中：DFT表示镜像函数；T为周期。

已知WF 转移函数取决于频谱函数ε-jkwT的复共轭S*，说明WF 与回波信号间具有最高功率转化的最优搭配，可由此求取输出信噪比最大值，通过WF可准确找出有效的回波信号。

1.4 基于模糊逻辑的天气雷达回波信号优化检测

以功率法为基础，通过模糊逻辑方法对利用匹配滤波器完成优化的天气雷达回波信号进行检测。由于特征分布在许多区域会发生重叠，单一特征参数的划分效度不够，遂选取多个可以有效描述回波特征的物理参数，依据参数的概率分布构建对应隶属函数，并设定回波信号的取值范围。通过联合多个参数可有效区分不同类型的回波信号。

首先在折射率因子内选取可以描述回波差异的特征参数，折射率因子的延展性OREF为

式中：划分效度函数

式中：Zj为方位角上某距离库的折射率因子；Val为有效检测值；AR为折射率因子的距离库数。

当前距离库回波功率的一致性DB为

式中：时间、空间一致性匹配函数为

式中：Oj为方位角上某距离库和天气雷达的距离。

局部折射率因子回波纹理的一致性FDBZ为

局部区间内邻近折射率因子的符号变化SPIN为

式中：折射率范围确认函数为

式中：Zthresh为两个距离库间折射率因子的变化临界值。

依据上述参数的概率分布构建对应隶属函数，研究不同特征参数差异。已知概率分布函数为量化某个距离库反射率因子值，是某种回波信号的置信度，所以，隶属函数形状大致类似于概率分布。按照特征参数OREF，DR，FDBZ及SPIN的概率分布，确定各参数的隶属函数。

针对OREF大于80%的回波，求取方位角上某距离库没有接受距离修改的Bj，然后运算该距离库所在回波最后10%的有效距离库的B的平均值，最后运算DB，将DR与所设阈值进行对比，低于阈值则判定此距离库为待检测的目标回波信号［9］。

构建包含不同类型天气雷达回波信号的数据库，通过统计数据库检测结果的临界成功指数（CSI）确定权重选择，具体表示为

选择权重时，需要应用先验约束条件缩小研究组数，先验约束条件包括：1）OREF，DB权重占比大于0；2）FDBZ，SPIN的权重占比按照回波信号在FDBZ与SPIN的概率分布重叠区域多少判断；3）单一决策的权重占比不应超过40%，各特征参数权重需控制在总权重的40%以内。

2 实验分析

以某台雷达站天气雷达为实验对象，参数详情见表1。发射机以同倍率发射天线的方式进行高频脉冲电磁波发射，同口径接收天线的方式进行回波信号检测，接收机的增益为40 dB，带宽为80 MHz。天气雷达所检测回波信号的数据集内共有8 000 条回波信号，信噪比区间为0～18 dB，共包括4 种调制类型（2PSK，4PSK，8PSK，16QAM），载波频率的分布集中在800～850 MHz，1 400～1 450 MHz 与 2 400～2 450 MHz范围内，回波信号长度为800个点。

表1 天气雷达参数详情Tab.1 Details of weather radar parameters

利用本文方法检测0.5°和1.5°两个仰角下的天气雷达回波信号，检测结果如图1 所示。

图1 不同仰角下的天气雷达回波信号检测结果Fig.1 Detection results of weather radar echo signals at different elevations

通过图1 可以看出，本文方法能有效检测0.5°和 1.5°两个仰角的天气雷达回波信号，所检测的范围相近，通过本文方法检测的不同仰角天气雷达回波信号没有存在孤岛，也没有在方位角周围留存任意天气雷达回波信号，原因是本文方法在检测天气雷达回波信号时，集合了多个特征参量，在某一特征参量无法满足阈值的情况下，其它特征向量可继续进行检测，不易遗落孤立天气雷达回波信号。这说明本文方法可有效检测不同仰角的天气雷达回波信号，检测效果较好。

分析本文方法检测天气雷达回波信号过程中，差异湍流外尺度跟随高度的变化结果如图2 所示。可以发现，随着高度的增长，应用本文方法检测的回波信号的垂直湍流外尺度逐渐接近平稳，在10～105之间波动，说明其惯性区的最大尺度较为稳定，可完整获取湍流强度和空间特征尺度的相关参数。

图2 垂直大气湍流外尺度随高度的变化Fig.2 Changes in the outer scale of vertical atmospheric turbulence with height

通过本文方法检测天气雷达回波信号过程中，研究温度变化条件下垂直大气湍流的折射率结构参数随高度的变化情况，如图3 所示。

图3 大气折射率结构参数随高度的变化Fig.3 Changes in atmospheric refractive index structural parameters with altitude

由于跟随高度提升变化最显著的指标是温度，其在垂直大气湍流环境下与高度呈反比例关系，温度的变化公式表示为

式中：I0为地表温度，℃；z为距离地面的高度，m。

根据图3 可以看出，与稳定大气环境下的折射率结构参数随高度的变化情况相比，垂直大气湍流折射率结构参数随高度的不断提升，以高度为200 m 为节点呈先下降后上升态势，在该节点垂直大气湍流折射率结构参数出现波谷，主要是因为温度受高度增长发生变化时，节点位置处的温度出现断点，该情况显著影响了垂直大气湍流折射率结构参数，但在高度为400 m 后呈微小程度增长且近似平稳。

结合图2、图3 可以看出，应用本文方法可有效改善垂直大气湍流对天气雷达回波信号的影响，从根本上优化垂直大气湍流外尺度与折射率结构参数受高度影响发生的变化，令其状态趋于平稳。应用本文方法对所采集天气雷达回波信号进行匹配滤波完成优化，结果如图4 所示。

图4 天气雷达回波信号优化结果Fig.4 Optimization results of weather radar echo signal

由图4 可以看出，初始采集的天气雷达回波信号波形存在抖动且抖动幅度较大，通过本文方法的匹配滤波进行滤波处理后，天气雷达回波信号波形均匀，且波动幅度小，优化效果较好。

对比原始天气雷达回波信号与本文方法检测的天气雷达回波信号，在差异天气雷达功率情况下的信噪比提升效果［10］如表2 所示。

表2 回波信号信噪比提升能力分析Tab.2 Analysis of signal to noise ratio improvement ability of echo signal

从表2 可以看出，本文方法检测的天气雷达回波信号较原始回波信号的信噪比存在较大程度的提升，两种天气雷达回波信号的信噪比均跟随天气雷达功率的减少不断降低，但利用本文方法检测的天气雷达回波信号的信噪比始终高于原始回波信号，且信噪比提升比例显著增加，说明本文方法的天气雷达回波信号检测结果较好，受垂直大气湍流影响程度较小，可有效扩充天气雷达在实际中的应用。

3 结论

本文提出垂直大气湍流影响下的天气雷达回波信号优化检测方法，通过匹配滤波优化垂直大气湍流影响下的天气雷达回波信号，可有效保留回波信号特征，增大输出信噪比，在检测天气雷达回波信号时帮助去除假目标，较好地提高后续检测精度，有效扩大了天气雷达回波信号检测的范围。