CN106019286A - 一种机载气象雷达多扫描多重频全空域气象目标探测方法 - Google Patents

一种机载气象雷达多扫描多重频全空域气象目标探测方法 Download PDF

Info

Publication number
CN106019286A
CN106019286A CN201610308458.2A CN201610308458A CN106019286A CN 106019286 A CN106019286 A CN 106019286A CN 201610308458 A CN201610308458 A CN 201610308458A CN 106019286 A CN106019286 A CN 106019286A
Authority
CN
China
Prior art keywords
spatial domain
radar
data
target
scanning
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201610308458.2A
Other languages
English (en)
Other versions
CN106019286B (zh
Inventor
石海杰
李京华
岳露
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201610308458.2A priority Critical patent/CN106019286B/zh
Publication of CN106019286A publication Critical patent/CN106019286A/zh
Application granted granted Critical
Publication of CN106019286B publication Critical patent/CN106019286B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/95Radar or analogous systems specially adapted for specific applications for meteorological use
    • G01S13/953Radar or analogous systems specially adapted for specific applications for meteorological use mounted on aircraft
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提供了一种机载气象雷达多扫描多重频全空域气象目标探测方法,涉及气象雷达信号处理技术领域,利用多扫描机载气象雷达,在垂直方向上采取低、中、高三波束的扫描方法,水平方向采取变密度扫描,航向及两侧区域采用较密集扫描,提高探测精度,远离航向区域采用适当稀疏扫描,提高扫描和后续数据处理的速度,对获取的三波束扫描回波信号进行合成,获取载机前方全空域气象目标分布情况,为飞行安全提供有力保障,本发明可提高机载气象雷达探测精度,增大探测距离,抑制地杂波干扰,获取全空域气象目标分布情况。

Description

一种机载气象雷达多扫描多重频全空域气象目标探测方法
技术领域
本发明涉及气象雷达信号处理技术领域,尤其是一种气象目标的探测方法。
背景技术
机载气象雷达是商业飞机强制安装的机载电子设备,能够实时探测飞机航路上的湍流、雷暴、冰雹等灾害性天气,也能探测起降阶段的风切变,可以称得上是飞机的“双眼”,对于保障飞行安全具有举足轻重的作用。随着我国越来越多航线的开通,商业航班更加繁忙,对于机载气象雷达的要求也越来越高。传统的单扫描机载气象雷达,由于受到雷达地平线的限制,只能探测到雷达波束与地球表面相切范围内的气象目标。同时受制于雷达波束宽度的限制,无法同时显示载机前方全空域气象目标分布情况。
为了能够探测全空域气象目标,传统方法是采用飞行员手动调节雷达波束仰角的方法。调低仰角,探测载机前下方近距离空域气象目标,调高仰角,探测远距离空域气象目标。该方法虽然可以获取到全空域的气象目标分布情况,但获取到的分布情况并不直观,需要飞行员根据经验在大脑中合成。这种方法过于依赖飞行员经验,增加了飞行员的工作复杂度,容易出错,甚至可能会对飞行员其它更重要的操作产生影响。由于无法抑制雷达地平线处的地杂波,该方法也不能探测更远距离气象目标分布情况。
多扫描技术是新型机载气象雷达的关键技术之一,其相关的数字信号处理技术是研制新型雷达的关键所在,现有的多扫描技术未能与多重频技术有效结合,虽然能够探测远距离空域气象目标,却没能实现近距离空域的高精度探测。与此同时,由于采用多扫描技术,数据量必然增加,将会影响处理速度。因此提出一种机载气象雷达多扫描多重频全空域气象目标探测方法,利用多扫描机载气象雷达,在垂直方向上采用低、中、高三波束进行扫描,探测载机前方近、远、中距离空域气象目标。三波束采用不同的脉冲重复频率,兼顾探测距离与探测精度。水平方向上采用变精度的扫描方法,载机航向及两侧区域上采用密集扫描,远离航向区域采用适当稀疏扫描,保证扫描精度同时尽量较少扫描和后续数据处理时间。对获取的三波束扫描回波信号进行合成,得到载机前方全空域气象目标分布情况,为飞行安全提供有力保障。
发明内容
为了克服现有技术的不足,本发明提供一种适用于机载气象雷达的多扫描多重频全空域气象目标探测方法,利用多扫描机载气象雷达,在垂直方向上采取低、中、高三波束的扫描方法,最低仰角波束采用高重复频率脉冲,提高探测精度,探测近距离空域范围内气象目标1;中仰角波束采用低重复频率脉冲,提高雷达无模糊作用距离,探测远距离空域范围内气象目标3;最高仰角波束采用中重复频率脉冲,探测中距离空域易受地物杂波干扰的气象目标2。水平方向采取变密度扫描,航向及两侧区域采用较密集扫描,提高探测精度,远离航向区域采用适当稀疏扫描,提高扫描和后续数据处理的速度。对获取的三波束扫描回波信号进行合成,获取载机前方全空域气象目标分布情况,为飞行安全提供有力保障。
本发明解决其技术问题所采用的技术方案包括下述步骤:
步骤一、根据载机高度和雷达波束垂直宽度确定低、中、高三波束的仰角:设θ为与地球表面相切波束的仰角,则θ=arccos[r/(r+h)],其中r为地球半径,h为载机高度,中仰角波束是与地球表面相切的波束,因此其仰角θ2=θ=arccos[r/(r+h)],最低仰角波束的仰角为最高仰角波束的仰角为其中为雷达波束垂直宽度,由雷达性能参数决定的;
步骤二、根据机载气象雷达性能参数确定近距离空域、中距离空域、远距离空域的范围,确定低、中、高三波束的脉冲重复频率,设Re为雷达探测距离,则近距离空域是指:中距离空域是指:远距离空域是指设低、中、高三波束对应的脉冲重复频率为f1、f2、f3,设临时频率变量则f2在PRF1与ftemp之间任取一固定值,f3=2f2,f1=2f3,其中c=3.0X108m/s,为大气中光速,PRF1为雷达性能参数中的脉冲重复频率最低值;
步骤三、调节天线仰角到θ1位置,即最低仰角波束的位置,设置脉冲重复频率为f1,完成一个周期即水平扫描角从0°到180°的水平扫描,存储扫描数据,存储的扫描数据为Data1=(R11,γ,Pr1),其中R1为雷达到目标点的斜距,R1取值范围在近距离空域范围之内,即θ1为最低仰角波束的仰角,γ为水平扫描角,取值范围为0到180°,Pr1为最低仰角波束回波功率;
步骤四、调节天线仰角到θ2位置,即中仰角波束的位置,设置脉冲重复频率为f2,完成一个周期即水平扫描角从0°到180°的水平扫描,存储扫描数据,存储的扫描数据为Data2=(R22,γ,Pr2),其中R2为雷达到目标点的斜距,R2的有效取值范围在近距离空域与远距离空域范围之内,即内,θ2为中仰角波束的仰角,Pr2为中仰角波束回波功率;
步骤五、调节天线仰角到θ3位置,即最高仰角波束的位置,设置脉冲重复频率为f3,完成一个周期即水平扫描角从0°到180°的水平扫描,存储扫描数据,存储的扫描数据为Data3=(R33,γ,Pr3),R3为雷达到目标点的斜距,R3的有效取值范围在近距离空域与中距离空域范围之内,即内,θ3为最高仰角波束的仰角,Pr3为最高仰角波束的回波功率;
步骤六、根据标准气象雷达方程和步骤三、四、五所存储的扫描数据Data1、Data2、Data3计算近、中、远距离空域的气象目标回波反射率因子:
标准气象雷达方程为其中,Z为气象雷达回波反射率因子,λ为雷达波长,R是雷达到目标的斜距,Pr为雷达回波功率,|K|2是由散射粒子介电性质确定的常数,G0为天线波束中心方向上的增益,Pt为雷达发射功率,经过计算后所得到的数据为,Data1_1=(R11,γ,Z1),Data2_1=(R22,γ,Z2),Data3_1=(R33,γ,Z3),Z1、Z2、Z3为回波的反射率因子;
步骤七、对步骤六计算得到的数据Data1_1、Data2_1、Data3_1做投影变换:
将所探测到的气象目标斜距投影到以地心为圆心、地心与载机连线为半径的圆弧上;
雷达波束与地球表面相切,设地球圆心为点O,气象目标位于点B,点A是载机所在位置,点A'是气象目标在载机航迹上的投影,则OA、OA'与地球表面和载机航迹分别垂直,则该气象目标点的斜距为r2=AB,对应的投影弧长为在ΔAOB中,设边长
d = O B = r 2 2 + ( r + h ) 2 - 2 r 2 ( r + h ) c o s ( 90 - θ )
则内角α=∠AOA'=arccos[((r+h)2+d2-r2 2)/2d(r+h)]
α同时也是l2对应的圆心角,可得
l2=2πα(r+h)/360
计算后得到的数据为Data1_2=(L1,γ,Z1),Data2_2=(L2,γ,Z2),Data3_2=(L3,γ,Z3),L1、L2、L3分别为R1、R2、R3在载机航迹上的投影;
步骤八、根据步骤六计算得到的回波反射率因子Z1、Z2、Z3对近距离空域气象目标强度进行合成:
近距离空域气象目标强度由低、中、高三波束回波反射率因子进行合成,其合成采用加权平均的方法,权值可由各波束自身强度确定,近距离空域气象目标强度得到数据Data_near=(Ln,γ,Z),其中Ln取值范围为近距离空域,即
步骤九、根据步骤六计算得到的目标回波反射率因子Z1、Z2、Z3对中距离空域气象目标回波强度进行合成:
中距离空域气象目标强度Z=Z3,得到数据Data_middle=(Lm,γ,Z),其中Lm取值范围为中距离空域,即
步骤十、根据步骤六计算得到的目标回波反射率因子Z1、Z2、Z3对远距离空域气象目标回波强度进行合成:
远距离空域气象目标强度Z=Z2,得到数据Data_far=(Lf,γ,Z),其中Lf取值范围为远距离空域,即
步骤十一、根据步骤八、步骤九和步骤十中得到的Data_near、Data_middle、Data_far获取载机前方全空域气象目标分布情况:
D a t a = ( L , γ , Z ) = D a t a _ n e a r L ∈ 0 ~ 1 4 R e D a t a _ m i d d l e L ∈ 1 4 R e ~ 1 2 R e D a t a _ f a r L ∈ 1 2 R e ~ R e
其中,γ为水平扫描角,取值范围为0到180°,Z为气象目标反射率因子,L是气象目标到载机的距离,Data就是最终获得的机载气象雷达探测数据。
本发明的有益效果是利用多扫描机载气象雷达,在垂直方向上采用低、中、高三波束进行扫描,探测载机前方近、远、中距离空域气象目标,其中低、中、高三波束采用不同的脉冲重复频率,兼顾探测距离与探测精度;在水平方向上采用变密度的扫描方法、载机航向及两侧区域采用密集扫描、远离航向区域采用适当稀疏扫描,保证扫描精度同时尽量减少扫描和后续数据处理时间;对获取的三波束扫描回波信号进行合成,得到载机前方全空域气象目标分布情况,为飞行安全提供有力保障。本方法可提高机载气象雷达探测精度,增大探测距离,抑制地杂波干扰,获取全空域气象目标分布情况。
附图说明
图1是本发明机载气象雷达多扫描多重频全空域气象目标探测方法示意图。
图2是本发明水平变密度扫描示意图。
图3是本发明波束仰角计算原理示意图。
图4是本发明斜距投影计算方法示意图。
图5是一种机载气象雷达多扫描多重频全空域气象目标探测方法扫描结果示意图。
其中,101-载机;102-地球表面;103-近距离空域;104-气象目标1;105-低仰角波束;106-中仰角波束;107-高仰角波束;108-中距离空域;109-气象目标2;110-地物;111-远距离空域;112-气象目标3,201-航向及两侧区域;202-远离航向区域;203-水平扫描波束,301-载机高度;302-与地表相切波束仰角;303-地球半径;304-与地表相切波束;401-载机所在位置点;402-气象目标在载机航迹上的投影点;403-气象目标位置点;404-载机航迹;405-波束;501-区域D;502-区域E;503-区域F;504-区域A;505-区域B;506-区域C。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
本发明为一种机载气象雷达多扫描多重频全空域气象目标探测方法,在垂直方向上采取多波束的扫描方法,提高探测精度,增大探测距离,有效抑制地杂波干扰;水平方向采取变密度扫描,提高扫描和后续数据处理的速度。对获取的三波束扫描回波信号进行合成,获取载机前方全空域气象目标分布情况,为飞行安全提供有力保障。
附图1为机载气象雷达多扫描多重频全空域气象目标探测方法示意图。利用多扫描机载气象雷达,在垂直方向上采取低、中、高三波束的扫描方法,最低仰角波束(105)采用高重复频率脉冲,提高探测精度,探测近距离空域(103)范围内气象目标1(104);中仰角波束(106)采用低重复频率脉冲,提高雷达无模糊作用距离,探测远距离空域(111)范围内气象目标3(112);最高仰角波束(107)采用中重复频率脉冲,探测中距离空域(108)易受地物(110)杂波干扰的气象目标2(109)。
附图2为水平变密度扫描示意图。水平方向采取变密度扫描,航向及两侧区域(201)采用较密集扫描,提高探测精度,远离航向区域(202)采用适当稀疏扫描,提高扫描和后续数据处理的速度。对获取的三波束扫描回波信号进行合成,获取载机前方全空域气象目标分布情况,为飞行安全提供有力保障。
本发明为一种机载气象雷达多扫描多重频全空域气象目标探测方法,包括下述步骤:
步骤一、根据载机高度和雷达波束垂直宽度确定低、中、高三波束的仰角:附图3为波束仰角计算原理示意图,设θ(302)为与地球表面相切波束(304)的仰角,则θ=arccos[r/(r+h)],其中r(303)为地球半径,h(301)为载机高度,中仰角波束(106)是与地球表面(102)相切的波束,因此其仰角θ2=θ=arccos[r/(r+h)],最低仰角波束(105)的仰角为最高仰角波束(107)的仰角为其中为雷达波束垂直宽度,由雷达性能参数决定的;
步骤二、根据机载气象雷达性能参数确定近距离空域(103)、中距离空域(108)、远距离空域(111)的范围,确定低、中、高三波束的脉冲重复频率。设Re为雷达探测距离,则近距离空域(103)是指:中距离空域(108)是指:远距离空域(111)是指:设低、中、高三波束对应的脉冲重复频率为:f1、f2、f3,设临时频率变量则f2可以在PRF1与ftemp之间任取一固定值,为了使探测距离有一定余量,应尽量将f2的取值靠近PRF1,f3=2f2,f1=2f3,其中c=3.0X108m/s,为大气中光速,PRF1为雷达性能参数中的脉冲重复频率最低值;
步骤三、调节天线仰角到θ1位置,即最低仰角波束(105)的位置,设置脉冲重复频率为f1,完成一个周期即水平扫描角从0°到180°的水平扫描,存储扫描数据,存储的扫描数据为Data1=(R11,γ,Pr1),其中R1为雷达到目标点的斜距,R1取值范围在近距离空域范围之内,即θ1为最低仰角波束(105)的仰角,γ为水平扫描角,取值范围为0到180°,Pr1为最低仰角波束(105)回波功率;
步骤四、调节天线仰角到θ2位置,即中仰角波束(106)的位置,设置脉冲重复频率为f2,完成一个周期即水平扫描角从0°到180°的水平扫描,存储扫描数据,存储的扫描数据为Data2=(R22,γ,Pr2),其中R2为雷达到目标点的斜距,R2的有效取值范围在近距离空域与远距离空域范围之内,即内,θ2为中仰角波束(106)的仰角,γ为水平扫描角,取值范围为0到180°,Pr2为中仰角波束(106)回波功率;
步骤五、调节天线仰角到θ3位置,即最高仰角波束(107)的位置,设置脉冲重复频率为f3,完成一个周期即水平扫描角从0°到180°的水平扫描,存储扫描数据,存储的扫描数据为Data3=(R33,γ,Pr3),R3为雷达到目标点的斜距,R3的有效取值范围在近距离空域与中距离空域范围之内,即内,θ3为最高仰角波束(107)的仰角,γ为水平扫描角,取值范围为0到180°,Pr3为最高仰角波束(107)的回波功率;
步骤六、根据标准气象雷达方程和步骤三、四、五所存储的扫描数据Data1、Data2、Data3计算近、中、远距离空域的气象目标回波反射率因子:
标准气象雷达方程为其中,Z为气象雷达回波反射率因子,λ为雷达波长,R是雷达到目标的斜距,Pr为雷达回波功率,|K|2是由散射粒子介电性质确定的常数,G0波束中心方向上的增益,Pt为雷达发射功率,经过计算后所得到的数据为,Data1_1=(R11,γ,Z1),Data2_1=(R22,γ,Z2),Data3_1=(R33,γ,Z3),Z1、Z2、Z3为回波的反射率因子;
步骤七、对经过步骤六计算后得到的数据Data1_1、Data2_1、Data3_1做投影变换:
将所探测到的气象目标斜距投影到以地心为圆心、地心与载机连线为半径的圆弧上,附图4为斜距投影计算方法示意图,雷达波束(405)与地球表面(102)相切,设地球圆心为点O,气象目标位于点B(403),点A(401)是载机(101)所在位置,点A'(402)是气象目标在载机航迹(404)上的投影,则OA、OA'与地球表面(102)和载机航迹(404)分别垂直,则该气象目标点的斜距为r2=AB,对应的投影弧长为在ΔAOB中,设边长
d = O B = r 2 2 + ( r + h ) 2 - 2 r 2 ( r + h ) c o s ( 90 - θ )
则内角α=∠AOA'=arccos[((r+h)2+d2-r2 2)/2d(r+h)]
α同时也是l2对应的圆心角,可得
l2=2πα(r+h)/360
计算后得到的数据为Data1_2=(L1,γ,Z1),Data2_2=(L2,γ,Z2),Data3_2=(L3,γ,Z3),L1、L2、L3分别为R1、R2、R3在载机航迹(404)上的投影;
步骤八、根据步骤六计算得到的回波反射率因子Z1、Z2、Z3对近距离空域气象目标强度进行合成:
近距离空域气象目标强度由低、中、高三波束回波反射率因子进行合成,其合成采用加权平均的方法,权值可由各波束自身强度确定,近距离空域气象目标强度得到数据Data_near=(Ln,γ,Z),其中Ln取值范围为近距离空域(103),即
步骤九、根据步骤六计算得到的目标回波反射率因子Z1、Z2、Z3对中距离空域气象目标回波强度进行合成:
中距离空域气象目标强度Z=Z3,得到数据Data_middle=(Lm,γ,Z),其中Lm取值范围为中距离空域(108),即
步骤十、根据步骤六计算得到的目标回波反射率因子Z1、Z2、Z3对远距离空域气象目标回波强度进行合成:
远距离空域气象目标强度Z=Z2,得到数据Data_far=(Lf,γ,Z),其中Lf取值范围为远距离空域(111),即
步骤十一、根据步骤八、步骤九和步骤十中得到的Data_near、Data_middle、Data_far获取载机前方全空域气象目标分布情况:
D a t a = ( L , γ , Z ) = D a t a _ n e a r L ∈ 0 ~ 1 4 R e D a t a _ m i d d l e L ∈ 1 4 R e ~ 1 2 R e D a t a _ f a r L ∈ 1 2 R e ~ R e
其中,γ为水平扫描角,取值范围为0到180°,Z为气象目标反射率因子,L是气象目标到载机的距离,Data就是最终获得的机载气象雷达探测数据。
一种机载气象雷达多扫描多重频全空域气象目标探测方法,其一种具体实施例为:某一机载气象雷达发射机功率150W,天线增益35dB,波束宽度3.5°,三波束脉冲重复频率分别为200Hz、400Hz、800Hz,脉冲宽度依次为5μs、10μs、20μs,距离分辨能力为750m、1500m、3000m。近距离空域范围取为0—175Km,中距离空域范围取为175Km—350Km,远距离空域范围为取350Km—600Km。水平扫描时,载机航向方向120°范围内采用1.5°间隔进行扫描,其余范围采用3.0°间隔扫描。
附图5为扫描结果示意图,区域A(504)内气象目标的距离分辨率为750m,角度分辨率为1.5°,该区域是对载机(101)威胁最大的区域,其回波数据是由三波束合成得到。区域B(505)内气象目标的距离分辨率为1500m,角度分辨率为1.5°,其回波数据主要来自最高仰角波束。区域C(506)内气象目标的距离分辨率为3000m,角度分辨率为1.5°,其回波数据主要来自中仰角波束。区域D(501)内气象目标的距离分辨率为750m,角度分辨率为3°,其回波数据是由三波束合成得到。区域E(502)内气象目标的距离分辨率为1500m,角度分辨率为3°,其回波数据主要来自最高仰角波束。区域F(503)内气象目标的距离分辨率为3000m,角度分辨率为3°,其回波数据主要来自中仰角波束。
本发明主要解决了机载气象雷达对航路上气象目标探测的问题,可以提高机载气象雷达探测精度,增大探测距离,提高处理速度,抑制地杂波干扰,获取全空域气象目标分布,为飞行安全提供有力保障。

Claims (1)

1.一种机载气象雷达多扫描多重频全空域气象目标探测的方法,其特征在于包括下述步骤:
步骤一、根据载机高度和雷达波束垂直宽度确定低、中、高三波束的仰角:设θ为与地球表面相切波束的仰角,则θ=arccos[r/(r+h)],其中r为地球半径,h为载机高度,中仰角波束是与地球表面相切的波束,因此其仰角θ2=θ=arccos[r/(r+h)],最低仰角波束的仰角为最高仰角波束的仰角为其中为雷达波束垂直宽度,由雷达性能参数决定的;
步骤二、根据机载气象雷达性能参数确定近距离空域、中距离空域、远距离空域的范围,确定低、中、高三波束的脉冲重复频率,设Re为雷达探测距离,则近距离空域是指:中距离空域是指:远距离空域是指设低、中、高三波束对应的脉冲重复频率为f1、f2、f3,设临时频率变量则f2在PRF1与ftemp之间任取一固定值,f3=2f2,f1=2f3,其中c=3.0X108m/s,为大气中光速,PRF1为雷达性能参数中的脉冲重复频率最低值;
步骤三、调节天线仰角到θ1位置,即最低仰角波束的位置,设置脉冲重复频率为f1,完成一个周期即水平扫描角从0°到180°的水平扫描,存储扫描数据,存储的扫描数据为Data1=(R11,γ,Pr1),其中R1为雷达到目标点的斜距,R1取值范围在近距离空域范围之内,即θ1为最低仰角波束的仰角,γ为水平扫描角,取值范围为0到180°,Pr1为最低仰角波束回波功率;
步骤四、调节天线仰角到θ2位置,即中仰角波束的位置,设置脉冲重复频率为f2,完成一个周期即水平扫描角从0°到180°的水平扫描,存储扫描数据,存储的扫描数据为Data2=(R22,γ,Pr2),其中R2为雷达到目标点的斜距,R2的有效取值范围在近距离空域与远距离空域范围之内,即内,θ2为中仰角波束的仰角,Pr2为中仰角波束回波功率;
步骤五、调节天线仰角到θ3位置,即最高仰角波束的位置,设置脉冲重复频率为f3,完成一个周期即水平扫描角从0°到180°的水平扫描,存储扫描数据,存储的扫描数据为Data3=(R33,γ,Pr3),R3为雷达到目标点的斜距,R3的有效取值范围在近距离空域与中距离空域范围之内,即内,θ3为最高仰角波束的仰角,Pr3为最高仰角波束的回波功率;
步骤六、根据标准气象雷达方程和步骤三、四、五所存储的扫描数据Data1、Data2、Data3计算近、中、远距离空域的气象目标回波反射率因子:
标准气象雷达方程为其中,Z为气象雷达回波反射率因子,λ为雷达波长,R是雷达到目标的斜距,Pr为雷达回波功率,|K|2是由散射粒子介电性质确定的常数,G0为天线波束中心方向上的增益,Pt为雷达发射功率,经过计算后所得到的数据为,Data1_1=(R11,γ,Z1),Data2_1=(R22,γ,Z2),Data3_1=(R33,γ,Z3),Z1、Z2、Z3为回波的反射率因子;
步骤七、对步骤六计算得到的数据Data1_1、Data2_1、Data3_1做投影变换:
将所探测到的气象目标斜距投影到以地心为圆心、地心与载机连线为半径的圆弧上;
雷达波束与地球表面相切,设地球圆心为点O,气象目标位于点B,点A是载机所在位置,点A'是气象目标在载机航迹上的投影,则OA、OA'与地球表面和载机航迹分别垂直,则该气象目标点的斜距为r2=AB,对应的投影弧长为在ΔAOB中,设边长
d = O B = r 2 2 + ( r + h ) 2 - 2 r 2 ( r + h ) c o s ( 90 - θ )
则内角α=∠AOA'=arccos[((r+h)2+d2-r2 2)/2d(r+h)]
α同时也是l2对应的圆心角,可得
l2=2πα(r+h)/360
计算后得到的数据为Data1_2=(L1,γ,Z1),Data2_2=(L2,γ,Z2),Data3_2=(L3,γ,Z3),L1、L2、L3分别为R1、R2、R3在载机航迹上的投影;
步骤八、根据步骤六计算得到的回波反射率因子Z1、Z2、Z3对近距离空域气象目标强度进行合成:
近距离空域气象目标强度由低、中、高三波束回波反射率因子进行合成,其合成采用加权平均的方法,权值可由各波束自身强度确定,近距离空域气象目标强度得到数据Data_near=(Ln,γ,Z),其中Ln取值范围为近距离空域,即
步骤九、根据步骤六计算得到的目标回波反射率因子Z1、Z2、Z3对中距离空域气象目标回波强度进行合成:
中距离空域气象目标强度Z=Z3,得到数据Data_middle=(Lm,γ,Z),其中Lm取值范围为中距离空域,即
步骤十、根据步骤六计算得到的目标回波反射率因子Z1、Z2、Z3对远距离空域气象目标回波强度进行合成:
远距离空域气象目标强度Z=Z2,得到数据Data_far=(Lf,γ,Z),其中Lf取值范围为远距离空域,即
步骤十一、根据步骤八、步骤九和步骤十中得到的Data_near、Data_middle、Data_far获取载机前方全空域气象目标分布情况:
D a t a = ( L , γ , Z ) = D a t a _ n e a r L ∈ 0 ~ 1 4 R e D a t a _ m i d d l e L ∈ 1 4 R e ~ 1 2 R e D a t a _ f a r L ∈ 1 2 R e ~ R e
其中,γ为水平扫描角,取值范围为0到180°,Z为气象目标反射率因子,L是气象目标到载机的距离,Data就是最终获得的机载气象雷达探测数据。
CN201610308458.2A 2016-05-11 2016-05-11 一种机载气象雷达多扫描多重频全空域气象目标探测方法 Expired - Fee Related CN106019286B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610308458.2A CN106019286B (zh) 2016-05-11 2016-05-11 一种机载气象雷达多扫描多重频全空域气象目标探测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610308458.2A CN106019286B (zh) 2016-05-11 2016-05-11 一种机载气象雷达多扫描多重频全空域气象目标探测方法

Publications (2)

Publication Number Publication Date
CN106019286A true CN106019286A (zh) 2016-10-12
CN106019286B CN106019286B (zh) 2018-06-08

Family

ID=57099362

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610308458.2A Expired - Fee Related CN106019286B (zh) 2016-05-11 2016-05-11 一种机载气象雷达多扫描多重频全空域气象目标探测方法

Country Status (1)

Country Link
CN (1) CN106019286B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106405513A (zh) * 2016-11-19 2017-02-15 陕西长岭电子科技有限责任公司 基于多层扫描对比的地杂波抑制方法
CN110426709A (zh) * 2019-07-02 2019-11-08 中国航空工业集团公司雷华电子技术研究所 一种机载气象雷达空域自适应扫描控制方法
CN110456315A (zh) * 2019-08-29 2019-11-15 西安电子工程研究所 一种基于位置预测的参差重频目标探测方法
CN110609286A (zh) * 2018-06-14 2019-12-24 邓勇 一种智能消除雷达地物杂波的处理方法
CN111264037A (zh) * 2017-10-23 2020-06-09 IPCom两合公司 减少飞行器引起的干扰
CN113687360A (zh) * 2021-10-26 2021-11-23 南京恩瑞特实业有限公司 一种适合民航空管气象保障的相控阵天气雷达扫描策略

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5103233A (en) * 1991-04-16 1992-04-07 General Electric Co. Radar system with elevation-responsive PRF control, beam multiplex control, and pulse integration control responsive to azimuth angle
CN1230692A (zh) * 1998-03-30 1999-10-06 三菱电机株式会社 气象观测雷达系统
US6590520B1 (en) * 2002-02-04 2003-07-08 Lockheed Martin Corporation Method and system for determining air turbulence using bi-static measurements
US20060121893A1 (en) * 2004-12-03 2006-06-08 Tillotson Brian J System for measuring turbulence remotely
US20140253370A1 (en) * 2013-03-06 2014-09-11 Kabushiki Kaisha Toshiba Weather radar apparatus, observation sequence generation method, and observation sequence generation program

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5103233A (en) * 1991-04-16 1992-04-07 General Electric Co. Radar system with elevation-responsive PRF control, beam multiplex control, and pulse integration control responsive to azimuth angle
CN1230692A (zh) * 1998-03-30 1999-10-06 三菱电机株式会社 气象观测雷达系统
US6590520B1 (en) * 2002-02-04 2003-07-08 Lockheed Martin Corporation Method and system for determining air turbulence using bi-static measurements
US20060121893A1 (en) * 2004-12-03 2006-06-08 Tillotson Brian J System for measuring turbulence remotely
US20140253370A1 (en) * 2013-03-06 2014-09-11 Kabushiki Kaisha Toshiba Weather radar apparatus, observation sequence generation method, and observation sequence generation program

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106405513A (zh) * 2016-11-19 2017-02-15 陕西长岭电子科技有限责任公司 基于多层扫描对比的地杂波抑制方法
CN106405513B (zh) * 2016-11-19 2018-12-21 陕西长岭电子科技有限责任公司 基于多层扫描对比的地杂波抑制方法
CN111264037A (zh) * 2017-10-23 2020-06-09 IPCom两合公司 减少飞行器引起的干扰
CN111264037B (zh) * 2017-10-23 2022-06-14 IPCom两合公司 减少飞行器引起的干扰
CN110609286A (zh) * 2018-06-14 2019-12-24 邓勇 一种智能消除雷达地物杂波的处理方法
CN110426709A (zh) * 2019-07-02 2019-11-08 中国航空工业集团公司雷华电子技术研究所 一种机载气象雷达空域自适应扫描控制方法
CN110456315A (zh) * 2019-08-29 2019-11-15 西安电子工程研究所 一种基于位置预测的参差重频目标探测方法
CN110456315B (zh) * 2019-08-29 2023-03-07 西安电子工程研究所 一种基于位置预测的参差重频目标探测方法
CN113687360A (zh) * 2021-10-26 2021-11-23 南京恩瑞特实业有限公司 一种适合民航空管气象保障的相控阵天气雷达扫描策略

Also Published As

Publication number Publication date
CN106019286B (zh) 2018-06-08

Similar Documents

Publication Publication Date Title
CN106019286A (zh) 一种机载气象雷达多扫描多重频全空域气象目标探测方法
Heinselman et al. High-temporal-resolution capabilities of the national weather radar testbed phased-array radar
Pichugina et al. Doppler lidar–based wind-profile measurement system for offshore wind-energy and other marine boundary layer applications
CN110095784B (zh) 一种复杂环境影响下的海洋-低层大气激光传输建模方法
CN108387885A (zh) 基于激光雷达探测的晴空条件下飞机尾流特征参数反演方法
EP2369361A1 (en) Radar system and method
CN104597430B (zh) 机载气象雷达低空风切变三维回波仿真方法
CN104215946B (zh) 一种天地波混合雷达回波谱仿真方法
May The organization of convection in the rainbands of Tropical Cyclone Laurence
CN104076338B (zh) 基于数字高程和数字地表覆盖的机载雷达杂波仿真方法
CN106501802A (zh) 高分辨率多维协同昆虫迁飞雷达测量仪
CN104035095A (zh) 基于空时最优处理器的低空风切变风速估计方法
CN102721967A (zh) 基于风场扰动类型的空中目标发现方法
CN103487798A (zh) 一种相控阵雷达测高方法
CN107783118A (zh) 基于毫米波雷达的固定翼无人机多目标防撞系统的防撞方法
CN110726980A (zh) 一种机载相控阵天气雷达对地表杂波的分析方法
CN107783114A (zh) 旋翼无人机远距离复杂环境防撞毫米波雷达信号处理系统及方法
CN107783128A (zh) 基于毫米波雷达的固定翼无人机多目标防撞系统
CN104199002B (zh) 抑制风力发电旋叶虚警目标的雷达信息处理方法
CN111781603A (zh) 一种机载气象雷达地杂波抑制方法
Karunarathna et al. Narrow bipolar pulse locations compared to thunderstorm radar echo structure
CN106019242A (zh) 一种天基双基地雷达飞行状态的配置方法
Norton et al. Boundary layer structure and decoupling from synoptic scale flow during NAMBLEX
CN103729485A (zh) 一种基于dem数据的宽带雷达相干杂波仿真方法
CN104914420A (zh) 基于多通道联合自适应处理的低空风切变风速估计方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180608

Termination date: 20210511

CF01 Termination of patent right due to non-payment of annual fee