第3高 原 气 象0卷 第4期Vol.30 No.4
,2011年8月011Auust2gPLATEAU METEOROLOGY]():张志强,刘黎平.相控阵技术在天气雷达中的初步应用[J.高原气象,2011,3041102-1107
相控阵技术在天气雷达中的初步应用
1
张志强1-3, 刘黎平
(1.中国气象科学研究院国家灾害天气重点实验室,北京 100081;2中国科学院研究生院,北京 100049;
)3.国家气象信息中心,北京 100081
摘 要:由中国气象科学研究院国家灾害天气重点实验室与中国电子科技集团南京第14研究所联合研制的S波段相控阵天气雷达采用宽波束发射多波束接收,从而很大程度上缩短了雷达扫描周期,但是由于相控阵天气雷达其波束宽度与增益不是定值,而是随着仰角发生变化,同时相控阵天气雷达采用了脉冲压缩技术,因此需要对Probert-Jones提出的经典雷达气象方程中的波束宽度和增益参量进行订正,才能运用到对相控阵天气雷达反射率因子的估算中。本文通过对相控阵天气雷达波束宽度与增益特性的分析,在经典雷达气象方程的基础上给出了有关相控阵天气雷达的雷达反射率计算方法,分区域给出了详细的反射率因子估算参数,并在此基础上进行反射率的估算,为相控阵天气雷达的定标与定量测量提供理论基础。
关键词:相控阵气象雷达;雷达增益;波束宽度;雷达气象方程()文章编号:10000534201104110206P406--- 中图分类号:
文献标识码:A
1 引言
我国目前已经布网的新一代天气雷达大大提高]1-7
,但是由于了对灾害性天气的监测和预警能力[
新一代天气雷达受机械扫描体制的限制,造成扫描时间过长、采样层数过少等时空分辨率较低的缺陷。采用相控阵技术的天气雷达可以提高探测数据
8]
的时间分辨率,利用相控阵技术的天气雷达[可以
方程已经被广泛地运用在反射率算法中,这些算法用来从气象雷达采样信号中计算回波反射率。新一代天气雷达是采用机械扫描方式,而相控阵天气雷达采用的是相位扫描方式。对于相位扫描雷达来说其波束宽度和增益是随着仰角变化的,同时为了提高性能,相控阵天气雷达采用了脉冲压缩技术。这些内在的联系都影响着雷达的有效辐射功率、雷达可分辨单元的大小、雷达的平均返回功率以及反射率的估算。所以针对相控阵天气雷达对经典雷达气象方程的某些参数进行订正是很有必要的,由此可得到针对相控阵天气雷达的雷达气象方程参数。本文针对S波段相控阵天气雷达参数和扫描特性,将经典雷达气象方程应用到相控阵天气雷达中,得到相控阵雷达对云目标探测的雷达气象方程参数,为今后对该雷达的定标和定量测量提供必要的准备和理论依据。
使天气雷达的扫描时间从6min减少至几十秒钟,也可以对多个区域同时观测,从而大大地提高了对龙卷风、暴雨等灾害性天气的预警能力,并且缩短了预报时间。2006年,美国开始计划研制具有多功能的相控阵天气雷达(MutilefunctionPhased p
[]1
,MP,其主要任务是进行强RadarAR)9-1Array
天气的监测、空中交通管制和服务于国家安全等。
我国目前有关相控阵天气雷达的研究尚处于起步阶段,由中国气象科学研究院与中国电子科技集团南京第14研究所联合研制的S波段相控阵天气雷达正在进行出厂前的测试工作,如何实现相控阵雷达的定量测量,已成为当前首要解决的问题。由
[2]
在1Probert-Jones1962年提出的经典雷达气象
2 S波段相控阵天气雷达简介
由中国气象科学研究院研制的S波段相控阵天
气雷达天线采用波导窄边裂缝天线形式,以实现水平方向超低副瓣及在垂直方向形成相控阵多波束。
;改回日期:2010041920100601 收稿日期:----
()气象)项目“相控阵天气雷达系统关键技术研究”资助GYHY(QX)200762 基金项目:公益性行业(--
,男,北京人,博士生,主要从事天气雷达相关算法研究以及相关气象数据处理.1982-) 作者简介:张志强(
:ailzhiian82@cams.cma.ov.cn E-mqgg
4期张志强等:相控阵技术在天气雷达中的初步应用 1103
由32根窄边开缝的波导行馈组成,每根波导上开82个窄边缝,全阵共有32×82=2624个裂缝单元。
/馈线系统由移相器、TR组件和功分器等组成,收发共用32个组件。天线具体参数见表1。
表1 S波段相控阵参数
Table1 Parametersofhasedarraweatherradar py
2.7~3.0GHz 工作频率:
极化形式:水平极化
由32根窄边开缝的波导行馈组成,每根波导上开82个窄边缝,全阵共有32×82=2624个裂缝单元
1.57°±0.15° 波束宽度:水平:
垂直:接收3.12°±0.3°×4发射:8°±0.4°
,仰角方向同时4波束接收,覆盖空仰角扫描范围:0°5°~3
域为8°
)垂直波束宽度8G≥34.4dB(° 天线增益:宽波束发射:
)接收:垂直波束宽度3.G≥38.3dB(12°
最大副瓣电平:水平收/发SLL≤-35dB
垂直接收 SLL≤-25dB
为脉冲宽度。
相对于常规天气雷达,相控阵天气雷达的波束
形状是随仰角而发生变化的,但其方向图仍近似满足高斯分布;又由于相控阵雷达天线采用幅度加权的方法,将副瓣电平控制在-20dB以下,可以不用考虑副瓣的影响。因此Probert-Jones推导经典雷达气象方程时的假设对相控阵雷达仍然适用。经典雷达气象方程对相控阵天气雷达来说仍然适用。但是由于相控阵天气雷达的波束宽度和增益随着仰角发生变化,同时又由于脉压技术的运用,其脉压得益也是需要考虑的,因此需要对雷达气象方程中的波束宽度和增益量进行修正,才能运用到相
11]
。控阵天气雷达反射率的估算中[
3.1 相控阵雷达增益
对于相控阵天气雷达天线波束,由于采用相位
控制扫描,天线增益不再是固定不变,而是随着扫描角发生变化。取天线法向上的增益为G0,则任
13]
:意方向上的天线增益G[
夹角,0° 工作状态时,阵面竖起与水平方向成8
为了提高扫描速度,发射时采用1个8宽波束,接°收时采用4个3.窄波束接收。在水平方向采用12°机械扫描,这与S波段多普勒天气雷达一样;在垂直方向,采用相位控制扫描。相控阵天气雷达的这种扫描方式大大缩短了扫描时间,但是因为其扫描方式的不同,故在垂直方向上波束宽度不再是固定不变,而是随着仰角发生变化,增益也随着仰角而发生变化。同时,其收发波束的不同也造成收发增益的不同。因此需要在经典雷达气象方程的基础上,根据相控阵天气雷达的特征参数,得到适合相控阵天气雷达的反射率估算表达式。
0
G=G0cosθ,
λ
天线的发射增益为Gt,接收增益为Gr,则:
)
0
2
()2
)
G=Gcosθ,
)λ
2
0
Gcosθt=Gt0t,
λ
r
r0
r
2
()3()4
式中:GGt0为天线法向方向发射增益;r0为天线法向方向接收增益;θt为发射波束与天线法向方向夹角;θr为接收波束与天线法向方向夹角。
当波束指向垂直于天线阵面时,天线增益最
大,随着波束指向角偏离阵面法线的角度越大,天线增益越小,到了±4时,增益减少了1.5°51dB;到了±6时,增益减小了近3d0°B。天线增益的降低会使雷达天线在整个观察范围内的发射信号功率分配极不均匀,因此需要对天线增益进行修正。3.2 天线波束宽度
包括相控阵雷达在内的任何雷达增益有一个基本的表达式,即:
D=,
ΩB
其中:ΩB为波束的立体角。 由于ΩBD=4π为一常数,则:
3 相控阵天气雷达回波强度测量
robert-Jones 对于离散体目标P
典的雷达气象方程:
[12]
提出的经
4
1024ln22 r
,Z=232 2
)()PtcK×||λG(θτπββφ
()1()5
其中:PPG为天线发r为接收功率;t为发射功率;射增益;R为目标距离雷达的距离;λ为天线波长;为垂LLθa为雷达系统损耗;a为大气传播损耗;β
2
直波束宽度;φ为水平波束宽度;|K|=β
22,8-1
;m为复折射指数;c=3×10m·sτ2
m+2
θu0v0φ0
θuφv=cosθ
,
2
()6
式中:θθu、v分别为水平和垂直的波瓣宽度;u0、φ
1104
高 原 气 象 30卷
波束垂直方向与天线法线方向的夹角。3.5 相控阵天气雷达探测强度补偿
宽波束,° 由于相控阵天气雷达在发射时采用8
其同一个宽波束内能流密度相差是很大的,如果波束宽度内增益都按同一个值来估算反射率,必然造成反射率估算的误差。而且其接收时是采用4个窄波束进行接收,3.12°4个接收波束的特性也是不同的。上面推导的相控阵雷达气象方程只考虑了不同宽波束间的增益与波束宽度的变化情况,并没有考虑同一个宽波束中收发增益的变化情况,下面分情况讨论。把宽波束分成4个区域,如图1a所示,接收情况如图1b所示
。
θv0分别为天线法线方向水平和垂直的波瓣宽度;φ
为发射波束与天线法向的夹角。3.3 脉冲压缩
对于相控阵天气雷达,增加脉冲宽度可以增大作用距离,但是脉冲宽度的增加会使距离分辨率变差,影响测距精度。一个有效的方法就是采用脉冲
14]
。脉冲压缩的工作原理是:采用宽脉压缩技术[
冲发射,以提高发射机的平均功率,保证雷达的最大作用距离。接收时利用脉冲压缩技术,获得窄脉冲,从而提高距离分辨率。因而脉冲压缩技术是解决作用距离与距离分辨率之间矛盾的一种有效办法。
脉压技术采用相干累计技术,当采用宽脉冲发
射,脉冲宽度增加了N倍,返回信号的功率增加了
N2倍,等效于发射功率增加了N2倍,它相当于把一个扩大了N倍的宽脉冲分成N个小块连续发射,其取样体并没有随着脉冲宽度的增大而增大,仍和原窄脉冲的取样体积相同,因此其分辨率也没有发生减小。在此引入脉压得益GP来衡量由于采用脉冲压缩技术所带来的雷达性能的改善。脉冲宽度ω取脉冲压缩后的脉冲宽度时,脉压得益GP=
图1 天线波束示意图
()发射波束,()接收波束ab
Fi.1 Schematicdiaramofbeamofantenna. gg
())atransmitbeam,(breceivinbeam g
N2,相当于脉冲宽度减小了N倍,发射功率增加
了N2倍,则雷达性能提高了N倍。3.4 相控阵雷达的雷达气象方程
相控阵天气雷达的雷达气象方程是在上述讨论的基础上得到的,考虑到相控阵天气雷达在收发增益上的不同以及水平上采用机械扫描,垂直方向上采用相位控制扫描,分别考虑脉压得益、波束宽度与收发增益的变化以及收发端的系统损耗,将式()、()、()),整理得:带入式(3461
4
0
Z=2()32
λτπ|K|0c
下面先考虑天线的波束情况,阵面竖起与水平,对于发射夹角,其扫描范围为0方向成80°°5°~3宽波束中心有5个发射角:2,10,18,26,34,4个,,,发射角与天线法向方向所成夹角为-8°0°8°,,每个发射角又分为4个区域,因此把相控16°24°以2为区间分为1阵天气雷达扫描范围0°5°°8个~3,分别对其进行发射增益G区域(图2)t、接收增益补偿,具体见式()()
G8~10。r和波束宽度
θuvφ
]cosr
θu0v0φ)
LPra2
,()74Rff0
·
LLsstr
GpGtcosGrcosθθpt0t0r
]]
8-1
;式中:c=3×10m·sτ为脉压后的脉冲宽度;
B波束宽θu0、v0分别为接收波束水平和垂直的3dφ
度;PGP为脉压得益;Pt为发射功率;r为接收机
输出功率;θr为接收波束方向与天线法线的夹角;
Gθt为发射波束方向与天线法线的夹角;t0为天线法向方向发射增益;GRr0为天线法向方向接收增益;为目标距离雷达的距离;LLst、sr分别为雷达发射 和接收的系统损耗;La为大气传播损耗。
由于相控阵天气雷达在水平上是采用机械扫描,垂直上采用相位控制扫描,因此θθr、t分别为
图2 扫描区间示意图
Fi.2 Schematicdiaramofscaninterval gg
Gcosθt=Gt0t,
Gcosθr=Gr0r,()8
()9
4期张志强等:相控阵技术在天气雷达中的初步应用 1105
r
,=
θθuvu0v0φφ
()10
的宽波束,而同一个波° 由于发射波束是一个8
束内部增益的变化是很大的,假设在波束内部增益是按高斯递减的。如图1a所示,假设波束法向方向增益为G0,取同每个区域中线方向的增益为该区域的增益,则4个区域的增益分别为0.68G0,0.88G0,0.88G0,0.68G0.
,最终11) 将相控阵天气雷达各参数带入方程(得到相控阵雷达气象方程:
式中:分别为天线发射增益补偿角、天线θθθt、r、φ接收增益补偿角、天线波束宽度补偿角,3个角度都为与天线法线方向之间的夹角,由于天线只在垂直方向进行相位控制扫描,因此只考虑波束垂直方向与天线法向的夹角。
考虑到发射、接收以及波束宽度的波束补偿
)角,将式(化简得:7
4
1024ln20
Z=R2()32
λτπ|K|0c
2
Zd37+20loR+10loP0loLLsLagggBz=1r+1str
(/(0lo0loGp)-4-1ggff0)om+Cp .
()12
式中:R为距离;P10loLLsLa为gr为接收功率;str
)
LL)·11 .(/(cosPGGG]θ)ff)]
stsr p
a
r
t
r0t0
0
4
]θu0v0φ
t
表2 增益修正项随仰角变化
ainTable2 Chaneofcorrectionwithelevation gg
扫描区间/°补偿/dB 扫描区间/°补偿/dB
0~2 0.605 18~
20 0.605
2~4 0.585 20~22 1.755
4~6 1.691 22~24 1.788
6~8 1.681 24~26 0.705
8~10 0.556 26~28 0.745
10~12 0.556 28~30 1.915
12~4 1.681 30~32 1.705
14~16 1.691 32~34 0.915
16~180.58534~360.985
图3 2009年7月7日12:46相控阵天气雷达观测
())))未进行波束补偿,(进行波束补偿,(未进行波束补偿,(进行波束补偿,a1°b1°c2.5°d2.5°
())未进行波束补偿,(进行波束补偿e4°f4°
),Fi.3 Realdetectionofhasedarraweatherradarat12:46on7Jul2009.(a1°PPIbeforecomensation gpyyp
(),(),(),b1°PPIaftercomensationc2.5°PPIbeforecomensationd2.5°PPIaftercomensation ppp
(),()e4°PPIbeforecomensationf4°PPIaftercomensation pp
1106
高 原 气 象 30卷
加雷达的作用距离和保持高分辨率,采用了脉冲压缩技术。因此需要根据相控阵天气雷达的参数特性,以Probert-Jones经典雷达气象方程为基础,得到适合相控阵天气雷达的雷达估算参数,并以此为依据进行相控阵天气雷达的反射率估算。)24所联合研制的 ( 中国气象科学研究院与1
相控阵天气雷达,发射时采用1个8宽波束发射,°接收时采用4个3.的窄波束进行接收。在发射12°时宽波束内的增益情况是不同的,接收时4个窄波束的情况也是不同的。
)3 ( 根据相控阵雷达的扫描特性,在之前得出的经典雷达气象方程的基础上,给出了相控阵天气雷达的雷达气象方程简化式,并简单分析了不同区域相控阵天气雷达的探测差异,对相控阵天气雷达一次测试数据进行了补偿。
由于相控阵天气雷达还处于测试阶段,本文只是针对该雷达测量的初步应用,究竟如何衡量相控阵天气雷达探测数据的准确性,以及确定一种衡量相控阵天气雷达性能的方法是今后要解决的问题。参考文献
[]1 俞小鼎,王迎春,陈明轩,等.新一代天气雷达与强对流天气
]():预警[J.高原气象,2005,243753-761.
[]2OIQC技术在雷达反演VAD廓线资料退模 朱立娟,龚建东.
]():糊中的应用研究[J.高原气象,2006,255862-869.[]3X波段双偏振多普勒天气雷达 马学谦,董万胜,楚荣忠,等.
]):降雨估算试验[J.高原气象,2008,27(2382-391.[]4 杨毅,邱崇践,龚建东,等.同化多普勒雷达风资料的两种方
]():法比较[J.高原气象,2007,263247-256.
[]5 杨毅,邱崇践.利用多普勒雷达资料分析一次强降水的中尺
]():度流场[J.高原气象,2006,255925-931.
[]62008年8月1~2日滁州大暴雨双多普勒雷达三维风 周海光.
]():场反演的初步结果[J.高原气象,2009,2861422-1433.[]7 段云霞,邵爱梅,杨毅.利用双多普勒雷达资料对一次台风流
]():场结构分析[J.高原气象,2010,291187-196.[]8eslieRL.TheU.S.nationalweatherradartestbed(hased L p
):,Potentialimactonaviation[C].Prerints31stInarra -ppyternationalConferenceonRadarMeteorolo.AmericMeteo -gy,roloSoc2003:25-26.
[]9oulasEF,JamesFK,DusanSZ,etal.Proressreorton D ggp
)[thenationalweatherradartested(hasedarraC].21stIn - -py
(损耗补偿项;为脉压得益项,它是用来衡10loGp)g量脉冲压缩带来的雷达性能的改善。假设雷达的脉,此时的τ压比为128,则脉压得益项的值为42dBz/为脉压后的脉冲宽度,其值为宽脉冲宽度的1128,。因此实际的脉压得益值应为21dBzComp为增益修正项,是由波束偏离天线法向造成的增益改变以及同一个宽波束内部增益的改变组成。
8个区域,其增益修正项的值 把俯仰上分为1
如表2所示。
从得到的相控阵雷达气象方程知,在整个扫描区间的边缘,随着波束偏离天线法向,造成波束宽度的展宽及增益的减少,扫描强度相对于中心波束是偏弱的。同时,目标物距离雷达越远,造成取样体积的加大,导致扫描值明显偏弱,必然会造成一些小尺度的强对流回波探测的缺失。
3.6 相控阵天气雷达初步探测结果
针对上述分析,得到了一次相控阵天气雷达实际观测结果,图3为2009年7月7日12:46相控阵天气雷达在南京试验场的初步观测结果。由于雷达还在测试阶段,得到的只是一个测试数据。此次观测中,采用脉冲压缩,脉压比为128,雷达水平转,共取得了3个仰角的有效数动范围为135°15°~3
,,。由于相控阵天气雷达相据,分别为1°2.5°4°位控制扫描的特点,完成3层扫描只用了几十秒,大大缩短了扫描时间。根据波束偏离天线法向方向,,的不同对得到的数据进行补偿,其中图3ace分别为采用宽脉冲宽度,未考虑脉压得益以及波束补,,,偿计算1对应的反射率结果;图3°2.5°4°b,df
是采用脉压后的窄脉宽,考虑脉压得益、波束宽度,,和增益的变化计算1的反射率结果。可°2.5°4°
以看出,在补偿前,回波值的强度过强,考虑到脉冲压缩对雷达性能的影响,计算其脉压得益项为21。同时由于此次观测时天线水平放置,根据天dBz
线偏离法向的不同,其不同仰角的补偿情况是不同的。经过补偿后,回波值域在一个正常的范围,与当时观测的降水情况一致。此次观测仅仅是一次测试,获得的数据有限,该相控阵天气雷达探测的准确性还需要进一步研究。
4 小结
)1 ( 对于相控阵天气雷达,由于其采用相位
控制雷达扫描波束,其波束宽度不再像常规雷达一样保持不变,而是随着仰角发生变化的。由于其波束宽度的变化,雷达增益也发生变化。同时为了增
ternationalConferenceonIinteractiveInformationandProcess -):Am,inSstem(IIPS.SanDieoericMeteoroloSoc2005. gyg []D10oulasEF,JamesFK,DusanSZ,etal.Udateonthena -gp
()[tionalweatherradartestbedhasedarraC].23rdInterna - -pytionalConferenceonInteractiveInformationandProcessin g):,Sstem(IIPS.SanAntonioAmericMeteoroloSoc2007:15 y
4期
-20.
张志强等:相控阵技术在天气雷达中的初步应用 1107
,():JRoMeteorSoc1962,883485-495.Quart y
[]K,13norrJB.Weatherradareuationcorrectionforfreuenc qqy
aileandhasedarraradars.AerosaceandElectronicSs -gpypy [],():temsJ.IEEETrans2007,4331220-1227.
[]M14errillISkolnik.雷达手册[Z].北京:电子工业出版社,
2007:384-388.
[]11LeKD,PalmerRD,YuT Y,etal.Adativearrarocess -pyp
[informultiissionhasedarraradarC].33rdConference -m gpy :Am,RadarMeteorolo.AustraliaericMeteoroloSocon gy2007:3-12.
[]]12Probert-JonesJR.TheradareuationinmeteoroloJ. qgy[
PreliminarAlicationofPhasedArraTechnoloinWeatherRadar yppygy
1-31
ZHANGZhiIULiianin -q L -pg,g
(1.State KeLaboratoroSevere Weather,ChineseAcademoMeteoroloicalSciences,Beiin00081,China; y yf yf gjg1
2.GraduateUniversitoChineseAcademoSciences,Beiin00081,China; yf yf jg1
3.National MeteoroloicalInormationCentre,Beiin00081,China) gfjg1
:TbstractheStateKeLaboratorofSevereWeatherofCAMSandNaninResearchInstituteof A yyjg hasedElectronicsTechnolohasdeveloedthefirstSbandarraweatherofChina.Inordertoreduce - pgypy
,radarscanccletheSbandhasedarraweatherusesatransmittedwiderbeaminitchandnarrowerthe - ypyp
,inazimuth,andmultireceivedbeamstooverlatransmittedbeam.Howeverforusinthisscannin pgg ,b,smethodandulsecomressioneam widthandainischanedwithelevationotheProbertJones -ppggweatherradareuationshouldbecorrected.Thecharacteristicsofbeam widthandainareanalzed. qgy
,BasedontheProbertJonesweatherradareuationittheweatherradareuationofthisariveshased - -qqgpraweatherradarandtheestimationofreflectivitindifferentreions.Itatheoreticalarametersrovides yygpp foundationforthecalibrationandmeasurementofthehasedarraweatherradar. py
:;;ainKewordsPhasedarraweatherradarRadarBeam width;Weatherradareuation gyqy
第3高 原 气 象0卷 第4期Vol.30 No.4
,2011年8月011Auust2gPLATEAU METEOROLOGY]():张志强,刘黎平.相控阵技术在天气雷达中的初步应用[J.高原气象,2011,3041102-1107
相控阵技术在天气雷达中的初步应用
1
张志强1-3, 刘黎平
(1.中国气象科学研究院国家灾害天气重点实验室,北京 100081;2中国科学院研究生院,北京 100049;
)3.国家气象信息中心,北京 100081
摘 要:由中国气象科学研究院国家灾害天气重点实验室与中国电子科技集团南京第14研究所联合研制的S波段相控阵天气雷达采用宽波束发射多波束接收,从而很大程度上缩短了雷达扫描周期,但是由于相控阵天气雷达其波束宽度与增益不是定值,而是随着仰角发生变化,同时相控阵天气雷达采用了脉冲压缩技术,因此需要对Probert-Jones提出的经典雷达气象方程中的波束宽度和增益参量进行订正,才能运用到对相控阵天气雷达反射率因子的估算中。本文通过对相控阵天气雷达波束宽度与增益特性的分析,在经典雷达气象方程的基础上给出了有关相控阵天气雷达的雷达反射率计算方法,分区域给出了详细的反射率因子估算参数,并在此基础上进行反射率的估算,为相控阵天气雷达的定标与定量测量提供理论基础。
关键词:相控阵气象雷达;雷达增益;波束宽度;雷达气象方程()文章编号:10000534201104110206P406--- 中图分类号:
文献标识码:A
1 引言
我国目前已经布网的新一代天气雷达大大提高]1-7
,但是由于了对灾害性天气的监测和预警能力[
新一代天气雷达受机械扫描体制的限制,造成扫描时间过长、采样层数过少等时空分辨率较低的缺陷。采用相控阵技术的天气雷达可以提高探测数据
8]
的时间分辨率,利用相控阵技术的天气雷达[可以
方程已经被广泛地运用在反射率算法中,这些算法用来从气象雷达采样信号中计算回波反射率。新一代天气雷达是采用机械扫描方式,而相控阵天气雷达采用的是相位扫描方式。对于相位扫描雷达来说其波束宽度和增益是随着仰角变化的,同时为了提高性能,相控阵天气雷达采用了脉冲压缩技术。这些内在的联系都影响着雷达的有效辐射功率、雷达可分辨单元的大小、雷达的平均返回功率以及反射率的估算。所以针对相控阵天气雷达对经典雷达气象方程的某些参数进行订正是很有必要的,由此可得到针对相控阵天气雷达的雷达气象方程参数。本文针对S波段相控阵天气雷达参数和扫描特性,将经典雷达气象方程应用到相控阵天气雷达中,得到相控阵雷达对云目标探测的雷达气象方程参数,为今后对该雷达的定标和定量测量提供必要的准备和理论依据。
使天气雷达的扫描时间从6min减少至几十秒钟,也可以对多个区域同时观测,从而大大地提高了对龙卷风、暴雨等灾害性天气的预警能力,并且缩短了预报时间。2006年,美国开始计划研制具有多功能的相控阵天气雷达(MutilefunctionPhased p
[]1
,MP,其主要任务是进行强RadarAR)9-1Array
天气的监测、空中交通管制和服务于国家安全等。
我国目前有关相控阵天气雷达的研究尚处于起步阶段,由中国气象科学研究院与中国电子科技集团南京第14研究所联合研制的S波段相控阵天气雷达正在进行出厂前的测试工作,如何实现相控阵雷达的定量测量,已成为当前首要解决的问题。由
[2]
在1Probert-Jones1962年提出的经典雷达气象
2 S波段相控阵天气雷达简介
由中国气象科学研究院研制的S波段相控阵天
气雷达天线采用波导窄边裂缝天线形式,以实现水平方向超低副瓣及在垂直方向形成相控阵多波束。
;改回日期:2010041920100601 收稿日期:----
()气象)项目“相控阵天气雷达系统关键技术研究”资助GYHY(QX)200762 基金项目:公益性行业(--
,男,北京人,博士生,主要从事天气雷达相关算法研究以及相关气象数据处理.1982-) 作者简介:张志强(
:ailzhiian82@cams.cma.ov.cn E-mqgg
4期张志强等:相控阵技术在天气雷达中的初步应用 1103
由32根窄边开缝的波导行馈组成,每根波导上开82个窄边缝,全阵共有32×82=2624个裂缝单元。
/馈线系统由移相器、TR组件和功分器等组成,收发共用32个组件。天线具体参数见表1。
表1 S波段相控阵参数
Table1 Parametersofhasedarraweatherradar py
2.7~3.0GHz 工作频率:
极化形式:水平极化
由32根窄边开缝的波导行馈组成,每根波导上开82个窄边缝,全阵共有32×82=2624个裂缝单元
1.57°±0.15° 波束宽度:水平:
垂直:接收3.12°±0.3°×4发射:8°±0.4°
,仰角方向同时4波束接收,覆盖空仰角扫描范围:0°5°~3
域为8°
)垂直波束宽度8G≥34.4dB(° 天线增益:宽波束发射:
)接收:垂直波束宽度3.G≥38.3dB(12°
最大副瓣电平:水平收/发SLL≤-35dB
垂直接收 SLL≤-25dB
为脉冲宽度。
相对于常规天气雷达,相控阵天气雷达的波束
形状是随仰角而发生变化的,但其方向图仍近似满足高斯分布;又由于相控阵雷达天线采用幅度加权的方法,将副瓣电平控制在-20dB以下,可以不用考虑副瓣的影响。因此Probert-Jones推导经典雷达气象方程时的假设对相控阵雷达仍然适用。经典雷达气象方程对相控阵天气雷达来说仍然适用。但是由于相控阵天气雷达的波束宽度和增益随着仰角发生变化,同时又由于脉压技术的运用,其脉压得益也是需要考虑的,因此需要对雷达气象方程中的波束宽度和增益量进行修正,才能运用到相
11]
。控阵天气雷达反射率的估算中[
3.1 相控阵雷达增益
对于相控阵天气雷达天线波束,由于采用相位
控制扫描,天线增益不再是固定不变,而是随着扫描角发生变化。取天线法向上的增益为G0,则任
13]
:意方向上的天线增益G[
夹角,0° 工作状态时,阵面竖起与水平方向成8
为了提高扫描速度,发射时采用1个8宽波束,接°收时采用4个3.窄波束接收。在水平方向采用12°机械扫描,这与S波段多普勒天气雷达一样;在垂直方向,采用相位控制扫描。相控阵天气雷达的这种扫描方式大大缩短了扫描时间,但是因为其扫描方式的不同,故在垂直方向上波束宽度不再是固定不变,而是随着仰角发生变化,增益也随着仰角而发生变化。同时,其收发波束的不同也造成收发增益的不同。因此需要在经典雷达气象方程的基础上,根据相控阵天气雷达的特征参数,得到适合相控阵天气雷达的反射率估算表达式。
0
G=G0cosθ,
λ
天线的发射增益为Gt,接收增益为Gr,则:
)
0
2
()2
)
G=Gcosθ,
)λ
2
0
Gcosθt=Gt0t,
λ
r
r0
r
2
()3()4
式中:GGt0为天线法向方向发射增益;r0为天线法向方向接收增益;θt为发射波束与天线法向方向夹角;θr为接收波束与天线法向方向夹角。
当波束指向垂直于天线阵面时,天线增益最
大,随着波束指向角偏离阵面法线的角度越大,天线增益越小,到了±4时,增益减少了1.5°51dB;到了±6时,增益减小了近3d0°B。天线增益的降低会使雷达天线在整个观察范围内的发射信号功率分配极不均匀,因此需要对天线增益进行修正。3.2 天线波束宽度
包括相控阵雷达在内的任何雷达增益有一个基本的表达式,即:
D=,
ΩB
其中:ΩB为波束的立体角。 由于ΩBD=4π为一常数,则:
3 相控阵天气雷达回波强度测量
robert-Jones 对于离散体目标P
典的雷达气象方程:
[12]
提出的经
4
1024ln22 r
,Z=232 2
)()PtcK×||λG(θτπββφ
()1()5
其中:PPG为天线发r为接收功率;t为发射功率;射增益;R为目标距离雷达的距离;λ为天线波长;为垂LLθa为雷达系统损耗;a为大气传播损耗;β
2
直波束宽度;φ为水平波束宽度;|K|=β
22,8-1
;m为复折射指数;c=3×10m·sτ2
m+2
θu0v0φ0
θuφv=cosθ
,
2
()6
式中:θθu、v分别为水平和垂直的波瓣宽度;u0、φ
1104
高 原 气 象 30卷
波束垂直方向与天线法线方向的夹角。3.5 相控阵天气雷达探测强度补偿
宽波束,° 由于相控阵天气雷达在发射时采用8
其同一个宽波束内能流密度相差是很大的,如果波束宽度内增益都按同一个值来估算反射率,必然造成反射率估算的误差。而且其接收时是采用4个窄波束进行接收,3.12°4个接收波束的特性也是不同的。上面推导的相控阵雷达气象方程只考虑了不同宽波束间的增益与波束宽度的变化情况,并没有考虑同一个宽波束中收发增益的变化情况,下面分情况讨论。把宽波束分成4个区域,如图1a所示,接收情况如图1b所示
。
θv0分别为天线法线方向水平和垂直的波瓣宽度;φ
为发射波束与天线法向的夹角。3.3 脉冲压缩
对于相控阵天气雷达,增加脉冲宽度可以增大作用距离,但是脉冲宽度的增加会使距离分辨率变差,影响测距精度。一个有效的方法就是采用脉冲
14]
。脉冲压缩的工作原理是:采用宽脉压缩技术[
冲发射,以提高发射机的平均功率,保证雷达的最大作用距离。接收时利用脉冲压缩技术,获得窄脉冲,从而提高距离分辨率。因而脉冲压缩技术是解决作用距离与距离分辨率之间矛盾的一种有效办法。
脉压技术采用相干累计技术,当采用宽脉冲发
射,脉冲宽度增加了N倍,返回信号的功率增加了
N2倍,等效于发射功率增加了N2倍,它相当于把一个扩大了N倍的宽脉冲分成N个小块连续发射,其取样体并没有随着脉冲宽度的增大而增大,仍和原窄脉冲的取样体积相同,因此其分辨率也没有发生减小。在此引入脉压得益GP来衡量由于采用脉冲压缩技术所带来的雷达性能的改善。脉冲宽度ω取脉冲压缩后的脉冲宽度时,脉压得益GP=
图1 天线波束示意图
()发射波束,()接收波束ab
Fi.1 Schematicdiaramofbeamofantenna. gg
())atransmitbeam,(breceivinbeam g
N2,相当于脉冲宽度减小了N倍,发射功率增加
了N2倍,则雷达性能提高了N倍。3.4 相控阵雷达的雷达气象方程
相控阵天气雷达的雷达气象方程是在上述讨论的基础上得到的,考虑到相控阵天气雷达在收发增益上的不同以及水平上采用机械扫描,垂直方向上采用相位控制扫描,分别考虑脉压得益、波束宽度与收发增益的变化以及收发端的系统损耗,将式()、()、()),整理得:带入式(3461
4
0
Z=2()32
λτπ|K|0c
下面先考虑天线的波束情况,阵面竖起与水平,对于发射夹角,其扫描范围为0方向成80°°5°~3宽波束中心有5个发射角:2,10,18,26,34,4个,,,发射角与天线法向方向所成夹角为-8°0°8°,,每个发射角又分为4个区域,因此把相控16°24°以2为区间分为1阵天气雷达扫描范围0°5°°8个~3,分别对其进行发射增益G区域(图2)t、接收增益补偿,具体见式()()
G8~10。r和波束宽度
θuvφ
]cosr
θu0v0φ)
LPra2
,()74Rff0
·
LLsstr
GpGtcosGrcosθθpt0t0r
]]
8-1
;式中:c=3×10m·sτ为脉压后的脉冲宽度;
B波束宽θu0、v0分别为接收波束水平和垂直的3dφ
度;PGP为脉压得益;Pt为发射功率;r为接收机
输出功率;θr为接收波束方向与天线法线的夹角;
Gθt为发射波束方向与天线法线的夹角;t0为天线法向方向发射增益;GRr0为天线法向方向接收增益;为目标距离雷达的距离;LLst、sr分别为雷达发射 和接收的系统损耗;La为大气传播损耗。
由于相控阵天气雷达在水平上是采用机械扫描,垂直上采用相位控制扫描,因此θθr、t分别为
图2 扫描区间示意图
Fi.2 Schematicdiaramofscaninterval gg
Gcosθt=Gt0t,
Gcosθr=Gr0r,()8
()9
4期张志强等:相控阵技术在天气雷达中的初步应用 1105
r
,=
θθuvu0v0φφ
()10
的宽波束,而同一个波° 由于发射波束是一个8
束内部增益的变化是很大的,假设在波束内部增益是按高斯递减的。如图1a所示,假设波束法向方向增益为G0,取同每个区域中线方向的增益为该区域的增益,则4个区域的增益分别为0.68G0,0.88G0,0.88G0,0.68G0.
,最终11) 将相控阵天气雷达各参数带入方程(得到相控阵雷达气象方程:
式中:分别为天线发射增益补偿角、天线θθθt、r、φ接收增益补偿角、天线波束宽度补偿角,3个角度都为与天线法线方向之间的夹角,由于天线只在垂直方向进行相位控制扫描,因此只考虑波束垂直方向与天线法向的夹角。
考虑到发射、接收以及波束宽度的波束补偿
)角,将式(化简得:7
4
1024ln20
Z=R2()32
λτπ|K|0c
2
Zd37+20loR+10loP0loLLsLagggBz=1r+1str
(/(0lo0loGp)-4-1ggff0)om+Cp .
()12
式中:R为距离;P10loLLsLa为gr为接收功率;str
)
LL)·11 .(/(cosPGGG]θ)ff)]
stsr p
a
r
t
r0t0
0
4
]θu0v0φ
t
表2 增益修正项随仰角变化
ainTable2 Chaneofcorrectionwithelevation gg
扫描区间/°补偿/dB 扫描区间/°补偿/dB
0~2 0.605 18~
20 0.605
2~4 0.585 20~22 1.755
4~6 1.691 22~24 1.788
6~8 1.681 24~26 0.705
8~10 0.556 26~28 0.745
10~12 0.556 28~30 1.915
12~4 1.681 30~32 1.705
14~16 1.691 32~34 0.915
16~180.58534~360.985
图3 2009年7月7日12:46相控阵天气雷达观测
())))未进行波束补偿,(进行波束补偿,(未进行波束补偿,(进行波束补偿,a1°b1°c2.5°d2.5°
())未进行波束补偿,(进行波束补偿e4°f4°
),Fi.3 Realdetectionofhasedarraweatherradarat12:46on7Jul2009.(a1°PPIbeforecomensation gpyyp
(),(),(),b1°PPIaftercomensationc2.5°PPIbeforecomensationd2.5°PPIaftercomensation ppp
(),()e4°PPIbeforecomensationf4°PPIaftercomensation pp
1106
高 原 气 象 30卷
加雷达的作用距离和保持高分辨率,采用了脉冲压缩技术。因此需要根据相控阵天气雷达的参数特性,以Probert-Jones经典雷达气象方程为基础,得到适合相控阵天气雷达的雷达估算参数,并以此为依据进行相控阵天气雷达的反射率估算。)24所联合研制的 ( 中国气象科学研究院与1
相控阵天气雷达,发射时采用1个8宽波束发射,°接收时采用4个3.的窄波束进行接收。在发射12°时宽波束内的增益情况是不同的,接收时4个窄波束的情况也是不同的。
)3 ( 根据相控阵雷达的扫描特性,在之前得出的经典雷达气象方程的基础上,给出了相控阵天气雷达的雷达气象方程简化式,并简单分析了不同区域相控阵天气雷达的探测差异,对相控阵天气雷达一次测试数据进行了补偿。
由于相控阵天气雷达还处于测试阶段,本文只是针对该雷达测量的初步应用,究竟如何衡量相控阵天气雷达探测数据的准确性,以及确定一种衡量相控阵天气雷达性能的方法是今后要解决的问题。参考文献
[]1 俞小鼎,王迎春,陈明轩,等.新一代天气雷达与强对流天气
]():预警[J.高原气象,2005,243753-761.
[]2OIQC技术在雷达反演VAD廓线资料退模 朱立娟,龚建东.
]():糊中的应用研究[J.高原气象,2006,255862-869.[]3X波段双偏振多普勒天气雷达 马学谦,董万胜,楚荣忠,等.
]):降雨估算试验[J.高原气象,2008,27(2382-391.[]4 杨毅,邱崇践,龚建东,等.同化多普勒雷达风资料的两种方
]():法比较[J.高原气象,2007,263247-256.
[]5 杨毅,邱崇践.利用多普勒雷达资料分析一次强降水的中尺
]():度流场[J.高原气象,2006,255925-931.
[]62008年8月1~2日滁州大暴雨双多普勒雷达三维风 周海光.
]():场反演的初步结果[J.高原气象,2009,2861422-1433.[]7 段云霞,邵爱梅,杨毅.利用双多普勒雷达资料对一次台风流
]():场结构分析[J.高原气象,2010,291187-196.[]8eslieRL.TheU.S.nationalweatherradartestbed(hased L p
):,Potentialimactonaviation[C].Prerints31stInarra -ppyternationalConferenceonRadarMeteorolo.AmericMeteo -gy,roloSoc2003:25-26.
[]9oulasEF,JamesFK,DusanSZ,etal.Proressreorton D ggp
)[thenationalweatherradartested(hasedarraC].21stIn - -py
(损耗补偿项;为脉压得益项,它是用来衡10loGp)g量脉冲压缩带来的雷达性能的改善。假设雷达的脉,此时的τ压比为128,则脉压得益项的值为42dBz/为脉压后的脉冲宽度,其值为宽脉冲宽度的1128,。因此实际的脉压得益值应为21dBzComp为增益修正项,是由波束偏离天线法向造成的增益改变以及同一个宽波束内部增益的改变组成。
8个区域,其增益修正项的值 把俯仰上分为1
如表2所示。
从得到的相控阵雷达气象方程知,在整个扫描区间的边缘,随着波束偏离天线法向,造成波束宽度的展宽及增益的减少,扫描强度相对于中心波束是偏弱的。同时,目标物距离雷达越远,造成取样体积的加大,导致扫描值明显偏弱,必然会造成一些小尺度的强对流回波探测的缺失。
3.6 相控阵天气雷达初步探测结果
针对上述分析,得到了一次相控阵天气雷达实际观测结果,图3为2009年7月7日12:46相控阵天气雷达在南京试验场的初步观测结果。由于雷达还在测试阶段,得到的只是一个测试数据。此次观测中,采用脉冲压缩,脉压比为128,雷达水平转,共取得了3个仰角的有效数动范围为135°15°~3
,,。由于相控阵天气雷达相据,分别为1°2.5°4°位控制扫描的特点,完成3层扫描只用了几十秒,大大缩短了扫描时间。根据波束偏离天线法向方向,,的不同对得到的数据进行补偿,其中图3ace分别为采用宽脉冲宽度,未考虑脉压得益以及波束补,,,偿计算1对应的反射率结果;图3°2.5°4°b,df
是采用脉压后的窄脉宽,考虑脉压得益、波束宽度,,和增益的变化计算1的反射率结果。可°2.5°4°
以看出,在补偿前,回波值的强度过强,考虑到脉冲压缩对雷达性能的影响,计算其脉压得益项为21。同时由于此次观测时天线水平放置,根据天dBz
线偏离法向的不同,其不同仰角的补偿情况是不同的。经过补偿后,回波值域在一个正常的范围,与当时观测的降水情况一致。此次观测仅仅是一次测试,获得的数据有限,该相控阵天气雷达探测的准确性还需要进一步研究。
4 小结
)1 ( 对于相控阵天气雷达,由于其采用相位
控制雷达扫描波束,其波束宽度不再像常规雷达一样保持不变,而是随着仰角发生变化的。由于其波束宽度的变化,雷达增益也发生变化。同时为了增
ternationalConferenceonIinteractiveInformationandProcess -):Am,inSstem(IIPS.SanDieoericMeteoroloSoc2005. gyg []D10oulasEF,JamesFK,DusanSZ,etal.Udateonthena -gp
()[tionalweatherradartestbedhasedarraC].23rdInterna - -pytionalConferenceonInteractiveInformationandProcessin g):,Sstem(IIPS.SanAntonioAmericMeteoroloSoc2007:15 y
4期
-20.
张志强等:相控阵技术在天气雷达中的初步应用 1107
,():JRoMeteorSoc1962,883485-495.Quart y
[]K,13norrJB.Weatherradareuationcorrectionforfreuenc qqy
aileandhasedarraradars.AerosaceandElectronicSs -gpypy [],():temsJ.IEEETrans2007,4331220-1227.
[]M14errillISkolnik.雷达手册[Z].北京:电子工业出版社,
2007:384-388.
[]11LeKD,PalmerRD,YuT Y,etal.Adativearrarocess -pyp
[informultiissionhasedarraradarC].33rdConference -m gpy :Am,RadarMeteorolo.AustraliaericMeteoroloSocon gy2007:3-12.
[]]12Probert-JonesJR.TheradareuationinmeteoroloJ. qgy[
PreliminarAlicationofPhasedArraTechnoloinWeatherRadar yppygy
1-31
ZHANGZhiIULiianin -q L -pg,g
(1.State KeLaboratoroSevere Weather,ChineseAcademoMeteoroloicalSciences,Beiin00081,China; y yf yf gjg1
2.GraduateUniversitoChineseAcademoSciences,Beiin00081,China; yf yf jg1
3.National MeteoroloicalInormationCentre,Beiin00081,China) gfjg1
:TbstractheStateKeLaboratorofSevereWeatherofCAMSandNaninResearchInstituteof A yyjg hasedElectronicsTechnolohasdeveloedthefirstSbandarraweatherofChina.Inordertoreduce - pgypy
,radarscanccletheSbandhasedarraweatherusesatransmittedwiderbeaminitchandnarrowerthe - ypyp
,inazimuth,andmultireceivedbeamstooverlatransmittedbeam.Howeverforusinthisscannin pgg ,b,smethodandulsecomressioneam widthandainischanedwithelevationotheProbertJones -ppggweatherradareuationshouldbecorrected.Thecharacteristicsofbeam widthandainareanalzed. qgy
,BasedontheProbertJonesweatherradareuationittheweatherradareuationofthisariveshased - -qqgpraweatherradarandtheestimationofreflectivitindifferentreions.Itatheoreticalarametersrovides yygpp foundationforthecalibrationandmeasurementofthehasedarraweatherradar. py
:;;ainKewordsPhasedarraweatherradarRadarBeam width;Weatherradareuation gyqy