南大仙林X波段雷达_强回波区_相关系数变低的问题_2013年8月21日.docVIP

南大仙林X波段雷达在雷暴强回波区相关系数变低的问题 2013年8月21日 现象 8月9日，南京大学仙林校区的雷达观测发现四周有雷暴，但发现雷暴中强回波区域的相关系数不正常，只有0.7。 这部雷达的工作模式：X波段，重复频率2000Hz（即重复周期为500μs），不参差，天线转速为3转/分，相关脉冲数为128。 天线仰角为4.3度时的回波图（正北为雷暴）如下： 注意：上图的速度范围为-16~+16m/s，而不是-8~+8m/s! 图表 1 8月9日，雷达观测到的雷暴回波图（PPI） 然后，我们立即对这个雷暴进行了RHI扫描，如下： 注意：上图的速度范围为-16~+16m/s，而不是-8~+8m/s! 图表 2 8月9日，雷达观测到的雷暴回波图（RHI） 从上面的回波图中，我们发现的奇怪现象有： 雷暴中心区的相关系数特别低，只有0.7 雷暴中心区强度的起伏也比较大，仅仅相邻几个像素，回波强度起伏就有20dB（从绿色到红色） 差分反射率的起伏也比较大，而且夹杂着负的值 速度图和KDP的图基本正常 发现这种奇怪现象之后，我们在8月10日、8月19日、8月20日都用这部雷达进行了雷暴的观测，都有上述现象。 和别的雷达的对比 8月19日，用在溧水的S波段双偏振雷达（重复频率为1000Hz）对南京地区的雷暴进行了观测，发现雷暴的相关系数的值都是红色（在0.95以上），是正常的。 （能不能弄到当时溧水雷达的回波图？） 在此之前，溧水的雷达有问题，没有进行观测。 国外X波段雷达对雷暴的探测结果 根据目前我们看到的文献，国外别的波段的雷达，对雷暴的探测，没有发现过相关系数下降的现象。 找到了一张，反映的是雷暴中相关系数下降的现象，但作者认为是闪电造成的： （赵坤电脑上的PPT中有一张图） （在赵坤给的2012年的讲义中有一张闪电引起相关系数下降的图） 另外，国外有专家分析过，不同波长的雷达，对不同直径的降水粒子观测时，对dBZ、ZDR需进行修正的曲线，如下： （赵坤电脑上的PPT中有几张图） 但我们没有发现，对相关系数要不要进行修正的文献。 对IQ数据的初步分析结果 在进行探测的时候，我们同时将IQ数据也记录下来了。下面分析8月10日的IQ数据的结果。 当时，在终端软件上的显示画面如下： 图表 3 8月10日，雷达终端软件上显示的雷暴回波图 从上图可以看出，在方位=60度，距离=50km处，有雷暴。但在雷暴区，强度起伏大、相关系数低。 首先，对采集记录的IQ数据（文件为“500us_EL=4.3deg_161000.h5”，采用HDF5格式），用Matlab重新进行了信号处理（Matlab编的信号处理算法和在雷达上实时运行的信号处理软件的算法是一模一样的），并绘制出图像如下（仰角为4.3度）： 图表 4 Matlab绘制的根据IQ数据计算出的强度图 图表 5 Matlab绘制的根据IQ数据计算出的速度图 图表 6 Matlab绘制的根据IQ数据计算出的ZDR 图表 7 Matlab绘制的根据IQ数据计算出的相关系数 对比Matlab计算出的结果图和雷达实时终端软件上显示的结果，可以发现两者是一致的，在雷暴的强回波区，都有同样的奇怪现象。 然后，我们对回波进行FFT分析。首先，先看相关系数正确的回波区域的频谱分析结果，结果如下：（每个距离点表示96m） 图表 8 某个相关脉冲组128个触发脉冲的幅度图和局部放大的频谱图 图表 9相关系数随距离的变化曲线 图表 10 SNR随距离的变化曲线 从上图可以看出，距离单元从600~625的相关系数是正确的（距离点604~610由于进行了SNR的订正，造成相关系数略大于1），频谱也是合理的。 下面就开始分析相关系数比较低的地方的IQ结果，结果如下：（每个距离点表示96m） 图表 11 某个相关脉冲组128个触发脉冲的幅度图和局部放大的频谱图 图表 12相关系数随距离的变化曲线 图表 13 SNR随距离的变化曲线 从上图可以看出，重点看540~560的频谱，发现频谱既宽，而且有很强的点。 画出550距离点（该距离处的相关系数=0.779）的频谱的一维曲线，如下： 图表 14 相关系数=0.779的距离点的频谱曲线 画出580距离点（该距离处的相关系数=0.956）的频谱的一维曲线，如下： 图表 15 相关系数=0.956的距离点的频谱曲线 对比这两张频谱曲线，我们发现在强回波，相关系数降低的地方，其频谱既宽，而且完全不符合高斯型，有特别强的点。 我们怀疑的几种可能原因 对于X波段雷达，大雨滴处于瑞利散射和米散射的交界区 当电磁波传播遇到空气分子和云、雨质点时，入射的电磁波会从这些介质或粒子上向四面八方传播相同频率电磁波，称散射。 假定大雨滴的半径为3mm，X波段的波长为3cm，则x=2*pi*r/波