小波与曲波变换探地雷达数据去噪对比分析

探地雷达勘探工程目标及观测工程环境越加复杂给其精确的数据处理带来极大挑战,高效的探地雷达数据去噪算法是当前关注的重要研究领域.基于阈值去噪思想,小波变换和曲波变换去噪算法在探地雷达数据去噪应用中受到限制,有必要开展上述两种去噪算法适用性和实用性系统评价及改进.基于传统高阶相关统计阈值和块状复数域阈值函数,本文开展了小波变换及曲波变换去噪算法在合成含噪数据去噪效果对比分析;提出窗口高阶相关统计阈值小波变换去噪算法、探讨了块状复数域阈值函数取值变化对曲波变换去噪效果影响规律.通过对实测数据去噪分析,验证了窗口高阶相关统计阈值小波变换和估计块状复数域阈值函数曲波变换去噪算法的可行性及有效性.     

中国雷达强天气预警信息质量控制与集成应用

冰雹和龙卷等强对流天气容易产生严重气象灾害,而雷达可以为强对流天气提供高时效的监测信息。针对当前雷达业务系统仅提供单站预警,且强天气识别算法缺乏本地优化等问题,研究了全国雷达预警信息组网集成技术,为大范围监测中小尺度强对流天气提供及时高效的参考信息;通过对全国高空站3年0℃层高度资料的统计分析表明,雷达监测冰雹等算法的0℃层高度应在3~6km范围动态调整,提出了雷达强天气预警动态本地化等综合质量控制方法。通过质控前后统计分析表明,质量控制可以有效降低雷达预警的虚警率,虚警剔除率可达到95%以上;与雷电观测资料的检验结果显示,质控后的雷达预警信息与雷电分布具有较好的契合度。将基于上述组网集成与质量控制的综合技术应用于2016年6月23日盐城龙卷强天气监测,取得了较好的应用效果。     

利用AMDAR数据分析厦门、泉州双雷达三维风场反演误差

利用2016年全年进出厦门机场和晋江机场航班的AMDAR(Aircraft Meteorological DAta Relay)数据对双雷达反演风场进行检验,分析了厦门、泉州双雷达风场反演的总体误差,研究结果表明:1)验证了双雷达连线附近区域反演风场的平均绝对误差较大,发现了与两部雷达连线的夹角小于15°和大于165°的区域的反演结果的平均绝对误差明显大于误差的年平均值。2)对于反演风向而言,误差随着高度增加而减小。对于反演风速而言,高度在9 km以下的反演误差在5 m/s左右,而9 km以上的反演误差较小。3)剔除了双雷达连线附近区域(与两部雷达连线的夹角小于15°和大于165°的区域)和反射率因子小于5 dBZ以及等于100 dBZ(剔除非气象回波)的反演结果后,双雷达反演风场误差较小,相对于AMDAR数据的风向年平均绝对误差为29.4°,风速年平均绝对误差为3.28 m/s,总体误差在可接受范围之内,反演结果接近"真值",该反演方法稳定可靠。     

GPM卫星和地面雷达对江苏盐城龙卷风强降水估测的对比

为综合评估卫星和天气雷达在2016年6月23日盐城龙卷风期间的强降水过程的降水估测精度,以国家级雨量站观测数据为基准,结合相关系数(CC)、相对误差(RB)、均方根误差(RMSE)以及分级评分指标,利用S波段的天气雷达定量降雨估测产品(RQPE)和全球降水观测计划多卫星融合产品(IMERG_FRCal,IMERG_FRUncal,IMERG_ERCal)进行比较。结果表明,雷达和卫星的累积降水量与雨量站的空间相关性很强(相关系数大于0.9),基本上能捕捉到整个降水过程的空间分布。降水主要分布在江苏省北部,但卫星高估了江苏省东北部强降水中心的降水量;对于小时时序区域平均降水,卫星高估了降水,而雷达低估了累积降水量。综合降水中心区域分析,IMERG的强降水区域降水量与雨量站的时间序列的偏差显著;RQPE在降水峰值达到之前及峰值之后与地面雨量站的变化趋势基本一致,但对降雨量峰值有明显的偏低。RQPE能较为准确地在时间上捕捉到降雨强度的变化趋势,但对于大雨及暴雨的估测能力不佳;RQPE的POD、SCI值都远远高于IMERG, FAR也较小。IMERG几乎未能监测到强降水的发生。总体上,RQPE对此次龙卷风强降水量的估测表现优于3种IMERG产品,特别是在捕捉强降水区域的空间分布方面,但对于强降水的估测能力仍需进一步改善。     

Ka波段云雷达晴空回波垂直结构及变化特征

利用2017—2019年中国气象局大气探测试验基地Ka波段云雷达资料,结合地面自动气象站、激光云高仪资料,从强度、速度、线性退极化比以及晴空回波高度等方面,分析晴空回波垂直结构和变化特征。基于激光和微波对粒子半径和数密度散射的差异,区分云和晴空回波。结果表明:Ka波段云雷达探测到的晴空回波在边界层主要包含层状湍流回波和点状昆虫回波,且回波顶高在3000 m以内。晴空回波强度和高度具有明显的季节和日变化特征,冬季回波顶高较低,夏季回波顶高较高,与地面气温具有很好的相关性,每年的1,2,11,12月几乎没有晴空回波,而7月和8月回波顶平均高度最高。晴空回波反射率因子为-40~-15 dBZ,其中层状湍流回波反射率因子概率密度峰值处反射率因子为-35 dBZ,点状昆虫回波反射率因子概率密度峰值处反射率因子为-30 dBZ。晴空回波垂直移动速度为-1.5~+0.5 m·s-1,整体呈下沉运动。层状湍流回波线性退极化比较点状昆虫回波稍大,一般为-10~-5 dB,点状昆虫回波线性退极化比一般为-15~-8 dB。     

衰减夹层GPR模拟的时频域全局反射误差

全局反射误差分析是深入研究探地雷达（GPR）吸收边界条件吸收效率的有力工具.基于常规完全匹配层（PML）的标准交错网格有限差分算法必须满足严格的CFL条件限制,即在单位时间步长内,不容许电磁波传播的距离超过单元网格尺寸.为了提高主区域所有网格节点的计算效率,并有效地吸收传播后期出现的低频隐失波,提出基于非分裂递归卷积完全匹配层（UCPML）的旋转交错网格（RSG）GPR正演算法.该算法采用不同的网格交错策略,并在边界条件中引入了吸收低频隐失波的自由可变因子,使得该算法允许选取较大的时间步长,提高了计算效率,并且实现了对低频隐失波的高效吸收.本文首先给出了RSG差分格式,推导了基于UCPML的RSG差分更新方程,探讨了数值色散的稳定性条件,然后以绕射现象严重的衰减夹层数值模拟为例,分别从波场快照、单道波记录、时间域反射误差（TDRE）、频率域反射误差（FDRE）四个方面分析了UCPML与常规PML的全局反射误差,说明了基于UCPML和RSG的GPR正演算法能更有效地吸收低频隐失波.

     

汶川地震同震形变场的GPS和InSAR邻轨平滑校正与断层滑移精化反演

为克服InSAR观测汶川地震同震形变场的邻轨不连续问题,提出联合GPS观测值与邻轨平滑约束的同震位移校正方法,采用GPS观测形变去除PALSAR轨道误差引入的残留平地相位,基于形变平滑条件校正邻轨干涉相位的不连续性.ALOS/PALSAR干涉处理结果表明,校正后同震形变场的准确度与平滑性得以显著提高,InSAR高相干点残差达3.6 cm,校正后精度提高约60%,低相干点精度提高约40%,校正后形变场的邻接平滑因子标准差减小约33%,验证了轨道误差校正与邻轨平滑约束方法的准确性与可靠性.进一步基于弹性半空间位错模型的断层滑移反演结果表明,断层滑移主要分布于映秀、北川和青川地区,集中于地壳深度0~16 km范围,最大滑动量（位于北川县城）约为9.0 m,GPS反演模型残差为5.5 cm,InSAR反演模型残差达9.2 cm,InSAR反演精度约有30%的显著提高,由模型反演计算得到的地震矩为8.0469×1020 N·m.     

应用雷达回波检验人工增雨效果的个例分析

利用 713-C型雷达对 1次飞机增雨过程进行了跟踪探测 ,发现有利于人工催化的云区 ( 20dBz) ,经过人工催化作业后回波增强到 30dBz并产生增雨。也有一部分云区 (小于 20dBz)虽然也进行了人工催化 ,但无明显增雨效果