X波段固态双极化天气雷达的金属球定标方法

文档序号:10510792阅读:597来源:国知局
X波段固态双极化天气雷达的金属球定标方法
【专利摘要】本发明提供了一种X波段固态双极化天气雷达的金属球定标方法,包括:在雷达宽脉冲探测范围和窄脉冲探测范围内分别选择一测试点;在每个测试点处,放置一金属球;针对每个金属球,利用其GPS信息与雷达经纬度和高度信息计算对应雷达观测仰角、方位角和所在库数;设置定标观测模式,获取观测数据,定标观测模式包括雷达观测方位的水平变化范围、和雷达观测仰角的变化范围和步进角度,观测数据包括雷达反射率因子观测值和差分反射率因子观测值;以及将雷达反射率因子观测值和差分反射率因子观测值分别与反射率因子理论值和差分反射率因子理论值相比较,以二者之间的差值作为雷达的测量误差,使用该测量误差修正雷达在相应探测范围内的探测结果。
【专利说明】
X波段固态双极化天气雷达的金属球定标方法
技术领域
[0001] 本发明涉及天气雷达,更具体地,涉及一种X波段固态双极化天气雷达的金属球定 标方法。
【背景技术】
[0002] 天气雷达在现代气象综合观测系统中起着非常重要的作用,对中小尺度天气系统 有很好的监测能力,为短临天气预报提供了直观的资料。随着双极化技术的应用,固态发射 机技术的成熟,天气雷达的性能不断改进。双极化天气雷达不仅能够得到雷达反射率因子 (Z)、径向速度(V)、速度谱宽(W),还可以得到差分反射率(Z DR)、零滞后互相关系数(P0)、差 分传播相移(Φ?Ρ)、单位差分传播相移(KDP)和线性退极化比(L DR)五个探测变量。探测变量 的增加为定量估测降水、降水粒子识别等资料应用提供更多的信息和判断依据,但同时要 求这些探测变量有较高的准确性。
[0003] 随着雷达气象学研究的发展,出现了不同的雷达定标技术。对回波强度的定标技 术主要包括:比较利用频移反射器得到发生频率偏移后对应的雷达回波;利用球载或机载 的标准物体产生的雷达回波进行定标,由于标准物的不稳定性会引起观测误差;在垂直方 向上,利用雨滴谱仪得到的雨滴谱分布计算出回波强度值,与雷达观测值相比较进行标定, 观测数据的格点匹配可能会引起误差。差分反射率的定标方法主要包括:信号源法适合标 校接收有源通道;太阳法只能标定雷达接收通路产生的Z DR系统误差而并未包括雷达发射通 路引入的ZDR系统误差,而且该方法只能在每天有限的时段内实施;小雨滴法能够标校整部 雷达的系统偏差,局限性在于标定过程必须在雨强很小的层状云降水条件下完成,难以满 足业务运行的要求。

【发明内容】

[0004] 本发明提供了一种X波段固态双极化天气雷达的金属球定标方法,包括:在雷达宽 脉冲探测范围和窄脉冲探测范围内分别选择一测试点;在每个测试点处,放置一金属球;针 对每个金属球,利用其GPS信息与雷达经炜度和高度信息计算对应雷达观测仰角、方位角和 所在库数;设置定标观测模式,获取观测数据,定标观测模式包括雷达观测方位的水平变化 范围、和雷达观测仰角的变化范围和步进角度,观测数据包括雷达反射率因子观测值和差 分反射率因子观测值;以及将雷达反射率因子观测值和差分反射率因子观测值分别与反射 率因子理论值和差分反射率因子理论值相比较,以二者之间的差值作为雷达的测量误差, 使用该测量误差修正雷达在相应探测范围内的探测结果。
【附图说明】
[0005] 为了更完整地理解本发明及其优势,现在将参考结合附图的以下描述,其中:
[0006] 图1示出了使用该定标方法对X波段固态双极化天气雷达进行定标的示意图。
[0007] 图2示出了根据本发明实施例的定标观测模式示意图。
[0008] 图3示出了根据本发明的实施例获取观测数据的示意图。
【具体实施方式】
[0009] 根据结合附图对本发明示例性实施例的以下详细描述,本发明的其它方面、优势 和突出特征对于本领域技术人员将变得显而易见。
[0010] 在本发明中,术语"包括"和"含有"及其派生词意为包括而非限制;术语"或"是包 含性的,意为和/或。
[0011] 在本说明书中,下述用于描述本发明原理的各种实施例只是说明,不应该以任何 方式解释为限制发明的范围。参照附图的下述描述用于帮助全面理解由权利要求及其等同 物限定的本发明的示例性实施例。下述描述包括多种具体细节来帮助理解,但这些细节应 认为仅仅是示例性的。因此,本领域普通技术人员应认识到,在不背离本发明的范围和精神 的情况下,可以对本文中描述的实施例进行多种改变和修改。此外,为了清楚和简洁起见, 省略了公知功能和结构的描述。此外,贯穿附图,相同参考数字用于相似功能和操作。
[0012] 发明人注意到,采用金属球作为标准物体对雷达进行定标有如下优点:相对于雷 达波长,选择金属球半径使得后向散射位于几何光学散射区域,则雷达后向散射面积仅取 决于金属球半径;金属球的Zdr值接近于零,由金属球圆率引起的Zdr值小于0. ldB;利用金属 球可以同时完成回波强度和差分反射率的定标。
[0013] 根据本发明的实施例,提供了一种X波段固态双极化天气雷达的金属球定标方法。 图1示出了使用该定标方法对X波段固态双极化天气雷达1进行定标的示意图。
[0014] 如图1所示,在X波段固态双极化天气雷达1的宽脉冲探测范围和窄脉冲探测范围 内分别选择一个测试点。例如,根据本发明的实施例,如雷达1的脉冲宽度为40ys,则可以以 到雷达1的距离6km为宽脉冲探测范围和窄脉冲探测范围之间的边界,并选取3km和7km位置 作为两个测试点。在晴好、静风(风速最好小于lm/s)的天气条件下,通过固定在绞车2上的 汽艇3,将金属球4悬挂于空中测试点位置处。汽艇3上悬挂全球定位系统GPS装置,可以实时 获取金属球的GPS信息,包括金属球所在的经炜度信息和高度信息。根据金属球GPS信息与 雷达站点信息计算出金属球所在仰角。
[0015] 根据本发明的实施例,利用金属球的全球定位系统GPS信息与雷达的经炜度和高 度信息计算对应的雷达观测仰角、方位角和所在库数包括:根据金属球所在的经炜度 (Lonl,Latl)和高度(hi)信息与雷达所在的经炜度(Lon2,Lat2)和高度(h2)信息,计算出雷 达观测仰角A1、方位角Az和金属球所在距离库数Num,计算公式为:
[0016] L = sin(Latl )*sin(Lat2)*cos(Lon2_Lonl )+cos(Latl )*cos(Lat2)
[0019] Az = arctan2(cos(Lat2)*sin(Lon2-Lonl),cos(Latl)*sin(Lat2)-sin(Latl)* cos(Lat2)*cos(Lon2_Lonl))
[0021] 其中,Re为地球半径,取值为6371km,c为光速,τ为雷达脉冲宽度。
[0022] 然后,设置雷达的定标观测模式,取得观测数据。所述定标观测模式包括雷达观测 方位的水平变化范围和步进角度、以及雷达观测仰角的变化范围的步进角度,所述观测数 据包括雷达的反射率因子观测值和差分反射率观测值。
[0023] 图2示出了根据本发明实施例的定标观测模式示意图。所述定标观测模式是指方 位上以金属球所在方位为中心角的变化范围,该变化范围的角度大于等于30°且小于等于 60°,并且仰角上以金属球所在仰角为中心,11层仰角,步进为小于等于0.2°。根据本发明的 实施例,雷达每完成一次观测,雷达实时控制程序就生成一个以时间、站点、文件形式等命 名的二进制数据文件,文件格式采用中国气象局气象探测中心下发的格式,该数据文件可 以存储站址信息、雷达性能参数、本次探测的观测参数、本次探测的回波数据等数据。
[0024] 图3示出了根据本发明的实施例获取观测数据的示意图。如图3所示,获取观测数 据可以包括:选择与金属球方位角相差1.5°以内、仰角相差0.4°以内、距离库数相差5个以 内的雷达观测数据中反射率因子的最大值,即为金属球对应的观测值,根据该值所在对应 的距离库数和层数得到差分反射率因子观测值。
[0025] 然后,将雷达的反射率因子观测值和差分反射率观测值分别与反射率因子理论值 和差分反射率理论值相比较,以二者之间的差值作为雷达的测量误差,根据该测量误差值 对雷达进行定标。例如,如果以反射率因子理论值减去反射率因子观测值所得差值作为反 射率因子测量误差,则在实际测量中,将反射率因子观测值与反射率因子测量误差相加得 到反射率因子修正值。如果以差分反射率因子理论值减去差分反射率因子观测值所得差值 作为差分反射率因子测量误差,贝U在实际测量中,将差分反射率因子观测值与差分反射率 因子测量误差相加得到差分反射率因子修正值。
[0026]所述反射率因子理论值Ζ是:
[0028] 其中:
[0029] λ为雷达波长,θ、φ为雷达水平和垂直方向上的3dB波束宽度,h为雷达径向分辨率, 即h = cT,c为光速,τ为雷达脉冲宽度,k为金属球的介电常数,r为金属球半径,R为金属球与 雷达间的距离,差分反射率因子的理论值为零。
[0030] 根据本发明的实施例,在测试之前,可以利用信号源产生与雷达相同频率且固定 能量的信号Pr,从雷达接收机测试端口贯入,经过高频放大和混频器、中频放大器和信号处 理机,在显示器上显示对应的雷达反射率因子Z,信号强度每次递增5dB,记录相应的Z值,统 计并分析Pr与Z的对应关系,检查是否为线性关系,以判断雷达接收机是否正常工作。
[0031] 根据本发明的实施例,在测试之前,利用频谱仪测试雷达发射机的脉冲持续时间, 脉冲形状,发射机峰值功率随时间的变化,检验其工作的稳定性。
[0032] 本发明解决了 X波段固态双极化天气雷达由于自身因素引起探测变量误差的难 题,实现对反射率因子和差分反射率同时定标,可以提高观测数据的质量,提高定量估测降 水的精度和降水粒子类型识别的准确性。本发明可广泛地应用于双极化天气雷达的定标, 特别是固态双极化天气雷达的定标。
[0033] 尽管已经参照本发明的特定示例性实施例示出并描述了本发明,但是本领域技术 人员应该理解,在不背离所附权利要求及其等同物限定的本发明的精神和范围的情况下, 可以对本发明进行形式和细节上的多种改变。因此,本发明的范围不应该限于上述实施例, 而是应该不仅由所附权利要求来进行确定,还由所附权利要求的等同物来进行限定。
【主权项】
1. 一种X波段固态双极化天气雷达的金属球定标方法,包括: 在雷达宽脉冲探测范围和窄脉冲探测范围内分别选择一个测试点; 在每个测试点处,放置一个金属球; 针对每个金属球,利用金属球的全球定位系统GPS信息与雷达的经炜度和高度信息计 算对应的雷达观测仰角、方位角和所在库数; 设置定标观测模式,获取观测数据,所述定标观测模式包括雷达观测方位的水平变化 范围、以及雷达观测仰角的变化范围和步进角度,所述观测数据包括雷达的反射率因子观 测值和差分反射率因子观测值;以及 将雷达反射率因子观测值和差分反射率因子观测值分别与反射率因子理论值和差分 反射率因子理论值相比较,以二者之间的差值作为雷达的测量误差,使用该测量误差修正 雷达在相应探测范围内的探测结果。2. 根据权利要求1所述的方法,其特征在于,还包括:在测试之前,利用信号源产生与雷 达相同频率且固定能量的信号Pr,从雷达接收机测试端口贯入,经过高频放大和混频器、中 频放大器和信号处理机,在显示器上显示对应的雷达反射率因子Z,信号强度每次递增5dB, 记录相应的Z值,统计并分析Pr与Z的对应关系,检查是否为线性关系,以判断雷达接收机是 否正常工作。3. 根据权利要求1所述的方法,其特征在于,还包括:在测试之前,利用频谱仪测试雷达 发射机的脉冲持续时间,脉冲形状,发射机峰值功率随时间的变化,检验其工作的稳定性。4. 根据权利要求1所述的方法,其特征在于:对于采用脉冲压缩技术的固态双极化雷 达,分别在3km、7km位置选取两个测试点。5. 根据权利要求1所述的方法,其特征在于,利用金属球的全球定位系统GPS信息与雷 达的经炜度和高度信息计算对应的雷达观测仰角、方位角和所在库数包括: 根据金属球所在的经炜度和高度信息与雷达所在的经炜度和高度信息计算出雷达观 测仰角、方位角及金属球所在距离库数。6. 根据权利要求1所述的方法,其特征在于,所述定标观测模式是指: 方位上以金属球所在方位为中心角的变化范围,该变化范围的角度大于等于30°且小 于等于60°;以及 仰角上以金属球所在仰角为中心,11层仰角,步进为小于等于〇. 2°。7. 根据权利要求1所述的方法,其特征在于,获取观测数据包括: 提取与金属球方位角相差1.5°以内、仰角相差0.4°以内、距离库数相差5个以内的雷达 观测数据中反射率因子的最大值,即为金属球对应的观测值,根据该值所在对应的距离库 数和层数得到差分反射率因子观测值8. 根据权利要求1所述的方法,其特征在于: 所述反射率因子理论值Z是:其中: λ为雷达波长,θ、φ分别为雷达水平和垂直方向上的3dB波束宽度,h为雷达径向分辨率, 即h = cT,c为光速,τ为雷达脉冲宽度,k为金属球的介电常数,r为金属球半径,R为金属球与 雷达间的距离, 差分反射率因子的理论值为零。
【文档编号】G01S7/40GK105866751SQ201610165887
【公开日】2016年8月17日
【申请日】2016年3月22日
【发明人】陈洪滨, 李兆明, 段树, 毕永恒, 陆建兵, 李东洋, 吴海军, 赖悦
【申请人】中国科学院大气物理研究所, 南京恩瑞特实业有限公司
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1