CN1114833C - 双极化雷达确定降水概率的方法及采用该方法的气象雷达 - Google Patents

双极化雷达确定降水概率的方法及采用该方法的气象雷达 Download PDF

Info

Publication number
CN1114833C
CN1114833C CN96199383A CN96199383A CN1114833C CN 1114833 C CN1114833 C CN 1114833C CN 96199383 A CN96199383 A CN 96199383A CN 96199383 A CN96199383 A CN 96199383A CN 1114833 C CN1114833 C CN 1114833C
Authority
CN
China
Prior art keywords
estimated value
radar
value
factor
decay
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.)
Expired - Fee Related
Application number
CN96199383A
Other languages
English (en)
Other versions
CN1206469A (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.)
Thales SA
Original Assignee
Thomson CSF SA
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 Thomson CSF SA filed Critical Thomson CSF SA
Publication of CN1206469A publication Critical patent/CN1206469A/zh
Application granted granted Critical
Publication of CN1114833C publication Critical patent/CN1114833C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

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
    • 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/024Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using polarisation effects
    • G01S7/025Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using polarisation effects involving the transmission of linearly polarised waves
    • 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)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及一种用双极化雷达确定降水概率的方法和一种使用该方法的气象雷达。在用雷达的两个极化通道H和V接收信号的基础上,本发明方法包括:首先确定单个极化通道的雷达射束方向上每个测距门i范围内的降水概率的估计值Ra(19,30到32),然后确定降水概率积分的第一估计值(I1),最后根据衰减反射因子(ZHa(r))确定未衰减差分反射因子(ZDRs(r))(35到38、47、48)。此外还要确定衰减差分反射因子的测量值(ZDRs(r))(30、40到42)并由此推导出降水概率的积分(I2)(43、44),然后将(I2)与第一估计值(I1)加以比较(45、46)以便当该方法只被应用于H极化通道时对该方法中的参数(α)进行调节。当α被调节以后就可得到未衰减反射因子的校准值和降水概率的校准值(38、47)。本发明适用于气象雷达。

Description

双极化雷达确定降水概率的方法及采用该方法的气象雷达
本发明涉及一种用脉冲雷达确定降水概率与水汽凝结体的反射因子的方法以及一种使用该方法的双极化气象雷达,其中,脉冲雷达的工作波长是被水汽凝结体衰减的波长。
在气象领域,特别是为了短期预报问题,需要观测降雨量和其它水汽凝结体。为此目的,许多年来人们一直在考虑用雷达来定量地确定降水概率,因为降水概率是这些现象的特征参量。通常采用的方法是基于用雷达射束来测量雨粒或水汽凝结体的反射率。然而,对大多数方法来说,水汽凝结体产生的散射会造成测量偏差,因此需要以衰减尽可能小的波长工作,即以相对长一些的波长工作。这增加了工作难度并增大了设备成本。
还提出另一个更具吸引力的方法:双波长衰减和单极化。这种方法的优点是不受雷达校准(calibration)的限制。不过,这种方法要获得需要的精度有一定的技术难度。这是因为在使两个波长的射束一致和使测量同步中会迂到一些问题。
现已进行的所有研究结果表明,用较短波长的雷达比较好,它使得可能选择体积更为紧凑、价格更加便宜的设备。不过,这会带来校准衰减问题,这些问题是非常棘手的,特别是当降雨量很大时,这些问题很难解决。
已经提出一种改进,包括利用双极化雷达来利用差分反射确定降水概率。这种类型的解决方法例如在E.GORGUCCI等人的文章“ARobust Estimator of Rainfall Rate Using Differential Reflectivity”(环境与海洋技术杂志,vol.11,1994年4月,第586-592页),或者在C.W.ULBRICH等人的“Assessment of the contribution of differentialpolarization to improved rainfall measurements”(无线电科学,vol.19,No.1,1984年1-2月,第49-57页)。这些解决方法改进了对于中等衰减得到的结果,但没有解决大降雨时的所有问题。
本发明涉及能够克服这些问题的方法,尽管采用会对穿过由水汽凝结体(雨、冰雹、云)构成的吸收介质的反射路径的电磁波产生相对较大衰减的短波长。
本发明基于利用双极化雷达测量两个极化分量 的反射率同时对衰减给采用的不同测量方法所测得的值造成的影响进行校准。
因此本发明提供一种用脉冲雷达来确定降水概率和水汽凝结体的反射因子的方法,其中,脉冲雷达的工作波长为被所述水汽凝结体衰减的波长,所述雷达是双极化雷达,它传送两串与各个极化通道相应的各种测距门的衰减反射因子[ZHa(r)、ZVa(r)]的测量值,该方法中:降水概率(R)是基于一个极化通道的反射因子(ZH)和未衰减的差分反射因子(ZDR)来估计的,所述方法的特征在于它包括下列步骤:
A)在雷达的每个测距门内,在测得的单极化反射因子的基础上用反射因子与降水概率的关系式和单极化衰减系数与降水概率的关系式来确定降水概率(Rs)的第一估计值;
B)在所述降水概率的第一估计值的基础上求得未衰减差分反射因子(ZDRs)的估计值;
C)从所述未衰减差分反射因子(ZDRs)的估计值和衰减后的差分反射因子的雷达测量值(ZDRa)推导出差分极化衰减估计值(ADP);和
D)通过反复调节反射因子与降水概率关系式中的一个参数(α)而同时使另一个参数(β)保持在一预定的常数值来设法使步骤A中估计出的单极化雷达的降水概率的积分[I1(r,R)]与由所述差分极化衰减的估计值(ADP)推导出的降水概率的积分[I2(r,R)]相等,以便能求出衰减反射因子的校准值[ZHs(r)]和降水概率[Rs(r)]。
本发明的显著优点在于它不受雷达校准的限制和不受由可能未探测到的云造成的衰减的限制,因为它是以差分测量为基础的。
本发明的另一个方面提供了一种以被水汽凝结体衰减的波长工作的双极化气象雷达,它用于确定降水概率和反射因子,所述雷达包括:
-双极化天线(1),它与提取与两个垂直的极化通道H和V相对应的模式的装置相关;
-用于向所述天线的两个通道发送脉冲的发射器;
-至少一个用于对两个垂直的极化通道的输出信号进行处理的接收通道;
所述雷达的特征在于它具有用来实现上述方法的数字信号处理设备。
借助于下列说明和附图可以更清楚地理解本发明的特征和优点,在附图中:
图1是按照本发明的气象雷达的示意图;和
图2是用于实现本发明方法的数字信号处理设备的方框图。
图1示出了双极化气象雷达的示意图。该雷达具有一对包括照明反射器(未示出)双极化初始源1的双极化天线。该双极化源与具有两个极化通道
Figure C9619938300122
的模式(正交模式)提取器2相关。脉冲发射器(b)向双极化源提供重复频率与所要求的范围相匹配的(例如,对150km范围来说,Fr=1KHz)脉冲。发射功率通过混合接头分布到天线的两个通道以便获得定向得与水平面成45°的极化,或者通过3dB耦合器5分布到天线的两个通道中以便获得圆极化。为了发射和接收两个相互垂直的极化的信号并根据它们的幅值对它们进行处理,可以事先将这两种发射模式设定为等效模式。圆极化波的第二个优点是天线反射波与发射器的发射波之间的比值是固定的。
提供了两个相同的接收通道,它们对由两个功率循环器3和4所抽取的信号进行处理。在传统的方法中,每个通道包括限制器7与8、放大器9与10、混合器11与12、中频放大器15与16、对数放大器17与18和模/数编码器19与20,其中混合器11、12用于使频率改变成中频,对数放大器17、18用于获得大的动态范围。经编码的信号被传送到根据本发明的数字信号处理设备21(该设备将在后面详细说明),然后数字信号处理设备输出校准后的反射因子Z和降水概率R信号作为其输出。混合器11和12还接收来自本地振荡器13的信号,本地振荡器13受自动频率控制电路14的控制。
在某些要求获得有关水汽凝结体性质的补充信息的气象应用场合,用另外的多普勒接收通道来测量被接收信号的相位更加有用。这些通道包括矢量解调电路22和25,一方面矢量解调电路的后面接有两个模/数编码器23、24,另一方面矢量解调电路25的后面接有两个模/数编码器26、27,用于输出正交编码的信号IH、QH(通道 )和Iv、Qv(通道 )。
本发明方法的基本原理如下:在单极化的基础上,例如
Figure C9619938300133
通道,距离为r处的雨粒(或其它水汽凝结体)的衰减反射因子ZHa(r)与衰减系数之间的关系式可用下式表示: Z Ha ( r ) = Z o ( r ) - 2 Δr Σ i = 1 n - 1 a i . . . . . . . . ( 1 )
此处Zo(r)为未衰减的反射因子,Δr为测距门的宽度,ai为第i个测距门的衰减系数,它以每单位长度dB为单位,Zha(r)和Zo(r)以dBZ为单位,也就是说等于10logZ(mm6m-3),其中r=nΔr。
衰减系数与降水概率R的关系用已知等式表示: a = k H R γ H . . . . . . . . . ( 2 )
此处KH和γH是极化通道 的参数,它们与温度和降水粒子分布(DSD)的关系不大,对应每个波长都可计算KH和γH
组合等式(1)和(2)可得出:
Zo(r)=ZHa(r)+2△rkHI1(r,R)  (3)及 I 1 ( r , R ) = Σ i = 1 n - 1 R i γ H . . . . . . . . ( 4 )
I1表示所选择的H极化通道的雷达射束方向上的降水概率的积分,也就是说,I1表示测距门n之前的n-1个测距门的累计衰减。
现在,己知雨粒反射因子Z与降水概率R之间的关系为:
              Z=αRβ……………(5)
此处α和β是参数。因此等式(4)可以用另一种等效形式表示: I 1 ( r , R ) = α - γ / β Σ i = 1 n - 1 Z i γ / β . . . . . . . . ( 4 ' )
参数α和β主要取决于DSD。在用迭带法测得的n个Za(r)值的基础上利用等式(3)和(5)中的四个参数K、γ、α和β的近似值可以计算出校准的估计值Z’(r),其中,迭带法是P.H.Hildebrand于1978年在“Iterative correction for attenuation of 5cm radar in rain”J.Appl.Meteor,17,第508-514页中提出的。正如该文章中所述的那样,DSD和温度的假定误差能大大减小这些衰减估计值,但这些减小量与由雷达校准误差所造成的减小量相比却很小。由于受雷达校准误差影响的参数是与Z和R有关的等式中的参数,因此I1(r,R)的主要不确定因素是由参数α和β的不确定因素所造成的。
为了判定这些不确定因素,利用使用该雷达的两个极化通道测得的差分测量结果。该差分反射因子ZDR由下式定义:
              ZDR=ZH-Zv……………(6)
此处Z以dBZ为单位,下标H和V对应于雷达的两个极化通道,如水平通道和垂直通道。
为了求得衰减波长,对应于H和V极化通道利用等式(1),然后从H项中减去V项得到:
              ZDRa(r)=ZDRs(r)-ZADP(r)……………(7)
此处ZDRa是衰减极化雷达测量值,ZDRs=ZOH-ZOV是仅由水滴形状造成的未衰减差分反射因子,ADP是电磁波路径上的差分极化衰减系数,即: A DP ( r ) = Δr Σ i = 1 n - 1 ( a iH - a iV ) . . . . . . . ( 8 )
对于假定的DSD,由水滴形状造成的差分反射因子ZDrs可以作为降水概率的函数计算:
              ZDrs=f(R)                  (9)
等式(8)中的衰减系数aiH和aiv可根据等式(2)对于H和V极化通道求出。假定DSD后,可以求出参数KH、KV、γH和γV。计算结果表明,差值aH-aV并不主要取决于DSD。对一给定的波长,KH与KV存在很大差异,但γH≈γV≈γ。
组合等式(2)和(8)可得到:
              ADP(r)=Δr(kH-kV)I2(r,R))  (10)
I 2 ( r , R ) = Σ i = 1 n - 1 R i r . . . . . . . ( 11 )
可注意到等式(11)可以写成类似于等式(4’)的形式,也就是说,作为Z的函数的形式。
I2(r,R)表示假定已正确地知道KH-KV时从差分衰减系数ADP推导出的降水概率的积分。由于ADP是从差分测量值推出的,因此ADP与雷达的校准无关。
由于存在雨滴DSD的自然变化和校准误差,因此在估计I1和I2时只有参数α和β是有效的。
按照本发明的方法中,建议将β保持为与气候平均值相等的恒定值并用差分衰减系数ADP来调整参数α。其中气候平均值与波长及降雨量的类型相对应。
有关文献中的许多测量结果证明对β的这种假定是有效的。
因此基本原理是:根据单个极化通道测得的一系列测量值确定I1,同时选择α的任一中间值;根据差分雷达测量值及利用I1估计未衰减的差分反射因子的计算值来确定I2;并将I1与I2加以比较以便近似时将调整α到使I1≈I2的值,这样就能求出未衰减反射因子的校准估计值Zi(r)和降水概率Ri(r)。
图2表示采用本方法的数字信号处理设备。图的左边部分涉及单个极化通道(H通道)的处理过程,而右边部分涉及差分处理过程。
来自H通道的编码器19的对数信号在步骤30中被重新转换成线性信号ZHa(r)。这些信号在雷达的N个脉冲上进行平均(步骤31),然后按照等式(5)基于参数α和β计算出降水概率估计值 Ra(步骤32),其中参数α和β是预定的,但α可任意选取。对每个测距门I来说都要进行上述计算。
步骤33包括根据等式(2)计算衰减系数 ai
然后根据等式(2)和(4)推导出值 I1(r,R)的估计值(步骤34)。
此外,在步骤36和37中根据单个极化通道H的一系列测量结果计算出衰减反射因子 ZHa(r),在步骤35中估计出等式(3)中的衰减系数 AH。在步骤38中应用等式(3)求和得到了H极化通道的衰减反射因子的估计值ZHs(r),然后在步骤47中,利用等式(5)和选定的参数α与β求出校准的降水概率 Rs(r)。
在这个估计值的基础上,可以在步骤48中用等式(9)计算出未衰减差分反射因子的估计值ZDRs(r)。
此外,根据步骤30和40传来的H通道和V通道的编码的信号可以得到比率Zha(r)/ZVa(r)(步骤41),给出将在步骤42中对N个脉冲作平均的衰减极化雷达测量值ZDRa(r)。该平均值与从步骤48传来的未衰减差分反射因子的估计值之间的差值构成了差分极化衰减系数的估计值 ADP(r)。在Δr和KH-KV值的基础上,步骤44输出由等式(10)得到的 I2(r,R)的估计值。然后在步骤45中对 I1和 I2加以比较。
如果 I1≠ I2,那么就在步骤46中使参数α在适当方向上改变Δα值并将新的α’值用到步骤32和47中重新进行计算。这些工作将重复进行直到 I1= I2为止,这时对应于给出α的正确值。然后就能求出未衰减反射因子的校准值Zi(r)= ZHs(r)和降水概率的校准值Ri(r)= Rs(r)(步骤49)。
在对本发明方法的基本原理加以说明的过程中,假定等式(9)中的ZDRs以及等式(10)中的KH、KV是为单个假定的DSD而计算的。这种假定是可以接受的,因为对给定的气候条件来说,平均的DSD是稳定的参数。因此,可以根据气候条件来选择特定的DSD。
此外,最好将测距门选择在反射因子较小的距离r处,因为这样会使ZDRs比较小,对由等式(9)得到的估计值造成的偏差也比较小。由于在不使I1、I2和α产生较大变化的情况下可以用另一个近似于r的且具有不同反射因子的测距门来代替在距离r处的测距门,因此可以更加方便、合理地选择特定的测距门。
可以在两个或三个相邻的测距门的基础上计算α,并且可以求出它们的平均值以便减小DSD分散的影响。
还要注意,所观测的辐射方向可以被分成几个部分,本发明方法可以同时应用于这几个部分。
上述方法对影响一个极化通道雷达所测得的反射因子或影响差分值的随机采样误差很不敏感。
现已说明了本发明方法的一个主要优点:这就是该方法与雷达的校准无关,与由未检测的云造成的衰减无关。
另外,另一种可能的用法是定性检测冰雹。现在对此加以说明:当雨粒中夹杂着冰雹时,由于冰雹大体上呈球形,它们对两个极化通道都会产生较高的反射因子并且ZH≈ZV,从而使差分反射因子基本上等于零。另一方面,会产生很大的与极化通道无关的衰减,因而使差分衰减系数ADP几乎为零。这样,冰雹的出现就改变了I1,但不会相应地改变I2,这与均匀雨粒的情况相反。因而能探测到α有相当大的变化,这可以用来作为冰雹的指示。
还应当注意,尽管上述说明是从利用差分极化衰减系数ADP的角度叙述的,但本发明方法还可以用来处理差分相位移KDP=KH-KV,KDP以传播角度Km为单位。
显然,所述例子不会对本发明产生限制。特别是,所述的气象雷达可以采用单个接收通道而不是所述的两个通道,两个极化通道的输入信号是采用切换方式输入的以便分时使用通道,输出信号则是用“不切换”方式输出的,只不过需要选择切换速率以便在切换期间极化通道H和V接收到的信号相互关联。
此外应注意可以用多普勒通道输出信号I和Q来计算每个通道H和V的信号幅值 I 2 + Q 2 ,这个幅值可以用来代替编码器19和20的输出信号以便不使用通道17、19与18、20。不过,这种做法存在这样的缺点:它不会象对数放大器那样提供动态范围。
最后,在图2所示的按照本发明的处理设备中,可以继续采用对数信号而不是将对数信号转换成线性信号。

Claims (15)

1、一种用脉冲雷达确定降水概率和水汽凝结体的反射因子的方法,其中脉冲雷达的工作波长为被所述水汽凝结体衰减的波长,所述雷达是双极化雷达,它传送两串与每个极化通道相应的各种测距门的衰减反射因子[ZHa(r)、ZVa(r)]的测量值,该方法中:降水概率(R)是基于一个极化通道的反射因子(ZH)和未衰减的差分反射因子(ZDR)来估计的,所述方法的特征在于它包括下列步骤:
A)在雷达的每个测距门内,在测得的单极化反射因子的基础上用反射因子与降水概率的关系式和单极化衰减系数与降水概率的关系式来确定降水概率(Rs)的第一估计值;
B)在所述降水概率的第一估计值的基础上求得未衰减差分反射因子(ZDRs)的估计值;
C)从所述未衰减差分反射因子(ZDRs)的估计值和衰减后的差分反射因子的雷达测量值(ZDRa)推导出差分极化衰减估计值(ADP);和
D)通过反复调节反射因子与降水概率关系式中的一个参数(α)而同时使另一个参数(β)保持在一预定的常数值来设法使步骤A中估计出的单极化雷达的降水概率的积分[I1(r,R)]与由所述差分极化衰减的估计值(ADP)推导出的降水概率的积分[I2(r,R)]相等,以便能求出衰减反射因子的校准值[ZHs(r)]和降水概率[Rs(r)]。
2、根据权利要求1所述的方法,其特征在于步骤A)包括下列步骤:
a)为了得到所述第一串测量值[ZHa(r)],根据任意设定参数(α、β)计算各个测距门(1到n-1)的降水概率的估计值(Ri);
b)从所述降水概率的估计值(Ri)推导出代表降雨所造成的衰减的第一估计值[I1(r,R)],所述第一估计值是为单个极化通道估算的降水概率的积分;
c)基于所述第一估计值[I1(r,R)]和第一串测量值[ZHa(r)]确定未衰减反射因子的估计值[ZHs(r)]和校准降水概率的估计值(Rs)。
3、根据权利要求1和2之一所述的方法,其特征在于步骤B)包括下列步骤:
d)基于所述校准降水概率的估计值(Rs)和给定降水粒子分布(DSD)计算未衰减差分反射因子的估计值[ZDRs(r)]。
4、根据权利要求1或2所述的方法,其特征在于所述步骤C)包括下列步骤:
e)基于所述两串衰减反射因子的测量值[ZHa(r)、ZVa(r)]计算衰减差分反射因子的测量值[ZDRa(r)];
f)根据衰减差分反射因子[的测量值ZDRa(r)]和在步骤d)中得到的未衰减差分反射因子的估计值[ZDRs(r)]推导出差分极化衰减的估计值(ADP)。
5、根据权利要求1或2所述的方法,其特征在于所述步骤D)包括下列步骤:
g)根据所述差分极化衰减(的估计值ADP)推导出代表降雨所造成的衰减的第二估计值[I2(r,R)],即降水概率的积分;
h)将所述第一[I1(r,R)]与表征降水衰减的所述第二[I2(r,R)]估计值比较;
i)如果步骤h)得出的比较结果表明所述第一估计值与所述第二估计值不一样,那么在步骤a)中改变至少一个所述参数(α、β)并重复步骤a)到h);
j)如果步骤h)得出的比较结果表明所述第一估计值与所述第二估计值相同,那么就求出所述未衰减反射因子的估计值[ZHs(r)]和所述校准降水概率的估计值(Rs)。
6、根据权利要求2所述的方法,其特征在于步骤a)包括:
a.1)在雷达发出N个脉冲的基础上求取所述第二串测量值[ZHa(r)]的平均值;
a.2)在任意设定参数(α、β)的基础上根据等式 R ‾ i = ( Z ‾ Ha ( r ) α ) 1 / β ,为每个测距门(i)计算降水概率的估计值[ Ri]。
7、根据权利要求2所述的方法,其特征在于步骤b)包括:
b.1)根据降水概率的估计值(Ri),用等式: ai/kH= Ri rH计算各个测距门的衰减系数(ai),此处,KH和γH是为某个温度和某个降水粒子分布预先设定的参数;和
b.2)根据这些衰减系数( ai),用等式: I ‾ 1 ( r , R ) = Σ i = 1 n - 1 a ‾ i / k H 推导出所述第一估计值( I1(r,R)),此处,n是从雷达到距离r处的测距门的数量。
8、根据权利要求2所述的方法,其特征在于步骤c)包括:
c.1)基于所述第一估计值(I1(r,R)),计算衰减值A,使得
A=2ΔrkH I1(r,R)
此处,Δr是测距门的宽度;
c.2)在雷达的n个脉冲上求取所述第一串测量值(ZHz(r))的平均值;
c.3)求所述衰减值 A的总和以及求衰减反射因子的平均值 ZHa(r)以便获得未衰减反射因子的估计值 ZHs(r);和
c.4)根据等式: R ‾ s = ( Z ‾ Hs ( r ) α ) 1 / β 计算校准降水概率的估计值 Rs
9、根据权利要求4所述的方法,其特征在于步骤e)包括:
e.1)求所述两串测量值(ZHa(r),ZVa(r))的比率以便获得雷达的衰减差分反射因子的值(ZDRa(r));
e.2)在雷达N个脉冲上求所述衰减差分反射因子的平均值。
10、根据权利要求5所述的方法,其特征在于步骤f)包括求取在步骤d)中得到的所述未衰减差分反射因子的估计值( ZDRs(r))与所述衰减差分反射因子的平均值( ZDRa(r))之间的差值以便获得雷达路径上差分极化衰减的估计值( ADP(r))。
11、根据权利要求5所述的方法,其特征在于步骤g)包括:根据所述差分极化衰减的估计值( ADP(r)),用等式 I ‾ 2 ( r , R ) = A ‾ DP ( r ) Δr ( k H - k V ) 计算所述第二估计值 I2(r,R),此处,KV是以与Kh同样的方式为第二串测量值(ZVa(r))预先设定的参数。
12、一种以被水汽凝结体衰减的波长工作的双极化气象雷达,它用于测定降水概率和反射因子,所述雷达包括:
-双极化天线(1),它与提取两个垂直的极化通道
Figure C9619938300062
Figure C9619938300063
的相应模式的装置(2,3,4)相关;
-用于向所述天线的两个通道发送脉冲的发射器(5,6);
-至少一个用于对两个垂直的极化通道的输出信号进行处理的接收通道(7,9,11,13,15,17,19;10,12,13,16,18,20);
所述雷达的特征在于它具有用来实现按照权利要求1或2所述的方法的数字信号处理设备(21),该数字信号处理设备(21)的输入端连接于该雷达中的模拟/数字编码器(19,21),该数字信号处理设备(21)的输出端连接该雷达的输出装置。
13、根据权利要求12所述的雷达,其特征在于,所述接收通道的输出部分包含用于增大动态范围的对数放大器(17、18),所述信号处理设备(21)的输入部分具有转换设备(30、36、40)以便恢复被接收的信号值。
14、根据权利要求12和13之一所述的雷达,其特征在于它还包括两个分别与两个极化接收通道的输出信号有关的多普勒接收通道(22到24;25到27);
15、根据权利要求12或13所述的雷达,其特征在于它具有单个的与受同步控制的输入和输出切换装置有关的接收通道,以便交替处理两个极化通道,其切换速率是选择得使在切换期间两个极化通道接收到的信号不是不相关联的。
CN96199383A 1995-12-26 1996-12-13 双极化雷达确定降水概率的方法及采用该方法的气象雷达 Expired - Fee Related CN1114833C (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR95/15483 1995-12-26
FR9515483A FR2742876B1 (fr) 1995-12-26 1995-12-26 Procede de determination du taux de precipitation par radar a double polarisation et radar meteorologique le mettant en oeuvre

Publications (2)

Publication Number Publication Date
CN1206469A CN1206469A (zh) 1999-01-27
CN1114833C true CN1114833C (zh) 2003-07-16

Family

ID=9485950

Family Applications (1)

Application Number Title Priority Date Filing Date
CN96199383A Expired - Fee Related CN1114833C (zh) 1995-12-26 1996-12-13 双极化雷达确定降水概率的方法及采用该方法的气象雷达

Country Status (9)

Country Link
US (1) US6061013A (zh)
EP (1) EP0870206B1 (zh)
JP (1) JP3720848B2 (zh)
KR (1) KR100442009B1 (zh)
CN (1) CN1114833C (zh)
CA (1) CA2242225C (zh)
DE (1) DE69604088T2 (zh)
FR (1) FR2742876B1 (zh)
WO (1) WO1997024629A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107643522A (zh) * 2017-09-19 2018-01-30 中国电子科技集团公司第三十八研究所 一种双频段双极化星载云雨测量雷达系统

Families Citing this family (76)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2742876B1 (fr) * 1995-12-26 1998-02-06 Thomson Csf Procede de determination du taux de precipitation par radar a double polarisation et radar meteorologique le mettant en oeuvre
US6164130A (en) * 1997-06-15 2000-12-26 Spherics Mess- Und Analysetechnik Gmbh Electromagnetic measuring system for meteorology
FR2774174B1 (fr) * 1998-01-23 2000-07-07 Centre Nat Rech Scient Technique pour l'estimation de pluie a partir d'un radar meteorologique a diversite de polarisation
US6563432B1 (en) * 2001-01-12 2003-05-13 Safegate International Ab Aircraft docking system and method with automatic checking of apron and detection of fog or snow
FR2827391B1 (fr) * 2001-07-11 2005-04-08 Centre Nat Rech Scient Technique pour l'estimation de pluie avec radar lorsque l'attenuation est negligeable
US6522972B2 (en) 2001-07-18 2003-02-18 Preston W Helms Method for determining an instantaneous unit hydrograph
US7049997B2 (en) * 2002-09-20 2006-05-23 Enterprise Electronics Corporation Simultaneous dual polarization radar system
US6803875B1 (en) * 2002-09-20 2004-10-12 Drs Weather Systems, Inc. Simulatneous dual polarization radar system
US6859163B2 (en) * 2002-09-20 2005-02-22 Drs Weather Systems, Inc. Simultaneous dual polarization radar system
JP3707469B2 (ja) * 2002-11-19 2005-10-19 三菱電機株式会社 二重偏波レーダ装置及びその信号処理方法
DE10305139B4 (de) * 2003-02-07 2006-02-16 Metek Metereologische Meßtechnik GmbH Verfahren und Vorrichtung zur Ermittlung von Niederschlagsarten in der Atmosphäre
OA13040A (en) * 2003-03-21 2006-11-10 Entpr Electronics Corp Simultaneous dual polarization radar system.
US7053813B1 (en) * 2004-04-22 2006-05-30 University Corporation For Atmospheric Research Radar system
US6879280B1 (en) * 2004-06-28 2005-04-12 Rockwell Collins, Inc. Vertical weather profile display system and method
JP2006184130A (ja) * 2004-12-27 2006-07-13 Tdk Corp レーダー装置
JP4762566B2 (ja) * 2005-02-15 2011-08-31 三菱電機株式会社 気象レーダ装置
FI117950B (fi) * 2005-05-23 2007-04-30 Vaisala Oyj Samanaikaisen kaksoispolarisoinnin tutkajärjestelmä, jossa on jalustaan asennettu vastaanotin
CN1967283B (zh) * 2005-11-18 2010-05-12 中国科学院物理研究所 一种用于气象雷达接收机前端的微波子系统
US7623064B2 (en) * 2005-12-06 2009-11-24 Arthur Robert Calderbank Instantaneous radar polarimetry
US7671785B1 (en) * 2005-12-15 2010-03-02 Baron Services, Inc. Dual mode weather and air surveillance radar system
US7551123B2 (en) * 2006-03-22 2009-06-23 Enterprise Electronics Corporation Phase shifted transmitted signals in a simultaneous dual polarization weather system
US7439899B2 (en) * 2006-03-22 2008-10-21 Enterprise Electronics Corporation Encoded transmitted signals in a simultaneous dual polarization weather system
US7369082B2 (en) * 2006-07-12 2008-05-06 Enterprise Electronics Corporation Method and apparatus implementing a scan strategy for automatic high power tube recovery
US7750573B2 (en) * 2006-07-12 2010-07-06 Enterprise Electronics Corporation Method and apparatus for automatic high power tube recovery
US7683828B2 (en) * 2006-07-12 2010-03-23 Enterprise Electronics Corporation System and method for measuring phase and power variance
DE102006054320A1 (de) * 2006-11-17 2008-05-21 Robert Bosch Gmbh Verfahren und Vorrichtung zum Erkennen von Niederschlag mittels Radar
JP4667426B2 (ja) * 2007-06-26 2011-04-13 三菱電機株式会社 気象レーダ装置
FR2925697B1 (fr) * 2007-12-21 2010-01-15 Thales Sa Procede de traitement d'une image radar.
DE112009000617B4 (de) * 2008-03-14 2021-03-11 Colorado State University Research Foundation Adaptive spezifische differenzielle Phase in Dual-Doppelpolarisationsradaren
EP2286511B1 (en) * 2008-05-07 2017-03-22 Colorado State University Research Foundation Dual-polarization radar processing system using time domain method
KR100921424B1 (ko) * 2008-11-13 2009-10-14 부경대학교 산학협력단 도플러 기상레이더의 관측 자료를 이용한 강수량 추정 방법
KR100942689B1 (ko) * 2008-11-13 2010-02-16 부경대학교 산학협력단 편파파라메터 차등위상을 이용한 편파레이더 강수감쇠 보정방법
JP5325361B2 (ja) * 2009-02-06 2013-10-23 古野電気株式会社 レーダ装置
US8264908B2 (en) * 2009-03-09 2012-09-11 Second Wind Systems, Inc. Method of detecting and compensating for precipitation in sodar systems
EP2278353B1 (en) * 2009-07-24 2018-02-14 Kabushiki Kaisha Toshiba Weather radar apparatus and rainfall rate calculation method
JP2011047742A (ja) * 2009-08-26 2011-03-10 Toshiba Corp 気象レーダシステムとその降水強度算出方法及びプログラム
JP2011027546A (ja) * 2009-07-24 2011-02-10 Toshiba Corp 気象レーダシステムとその降水量算出方法
JP5214562B2 (ja) * 2009-08-26 2013-06-19 株式会社東芝 気象レーダシステムとその降水強度算出方法及びプログラム
JP5491788B2 (ja) * 2009-07-24 2014-05-14 株式会社東芝 気象レーダシステムとその降水強度算出方法
JP2011047743A (ja) * 2009-08-26 2011-03-10 Toshiba Corp 気象レーダシステムとその降水強度算出方法及びプログラム
JP5135317B2 (ja) * 2009-11-04 2013-02-06 株式会社ホンダエレシス 車載レーダ装置、及びプログラム
US8988274B2 (en) * 2009-11-16 2015-03-24 The Board Of Regents Of The University Of Oklahoma Cylindrical polarimetric phased array radar
JP4832597B2 (ja) * 2011-01-11 2011-12-07 三菱電機株式会社 気象レーダ装置
US9261583B2 (en) * 2011-04-07 2016-02-16 Baron Services, Inc. Systems and methods for calibrating dual polarization radar systems
KR101221773B1 (ko) 2011-05-09 2013-01-11 경북대학교 산학협력단 이중편파 레이더를 이용한 기상 및 비기상 에코 분류 방법
US9182485B1 (en) * 2011-05-24 2015-11-10 Garmin International, Inc. Transmit/receive module for electronically steered weather radar
US8854250B2 (en) * 2011-06-28 2014-10-07 Vaisala Oyj Off-diagonal element echo power estimator for polarization weather radar
CN102662162B (zh) * 2012-02-16 2013-09-11 邓勇 X波段多普勒雷达降水估测方法
US9465106B2 (en) * 2012-03-27 2016-10-11 Colorado State University Research Foundation Robust attenuation correction system for radar reflectivity and differential reflectivity
KR101255736B1 (ko) 2012-09-21 2013-04-26 대한민국 단일편파 레이더를 이용한 기상 및 비기상 에코 분류 방법
CN102998670B (zh) * 2012-11-29 2014-03-26 北京无线电测量研究所 一种Ka频段固定指向双极化全固态毫米波云雷达
KR101294681B1 (ko) 2013-04-18 2013-08-08 (주) 웨더링크 기상 신호 처리장치 및 그 처리방법
JP6297037B2 (ja) * 2013-07-11 2018-03-20 古野電気株式会社 気象情報処理装置、気象レーダシステムおよび気象情報処理方法
KR101469860B1 (ko) * 2013-10-07 2014-12-08 부경대학교 산학협력단 기상 레이더 자료를 이용한 대류성 강우 특성 정보를 얻는 방법
KR101437587B1 (ko) * 2013-10-16 2014-11-03 부경대학교 산학협력단 정확한 강우량 추정을 위한 차등위상차 값의 보정방법
US10488556B2 (en) 2014-03-28 2019-11-26 Northeastern University System for multivariate climate change forecasting with uncertainty quantification
CN104330801B (zh) * 2014-11-15 2017-01-11 安徽四创电子股份有限公司 一种基于全数字阵列的有源相控阵天气雷达系统
US10026163B2 (en) 2015-02-25 2018-07-17 Cale Fallgatter Hydrometeor identification methods and systems
RU2648237C2 (ru) * 2015-04-27 2018-03-23 МИНИСТЕРСТВО ОБОРОНЫ РОССИЙСКОЙ ФЕДЕРАЦИИ ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ КАЗЁННОЕ ВОЕННОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ Военная академия ракетных войск стратегического назначения имени Петра Великого Способ измерения положения фронта волны
RU2611587C1 (ru) * 2015-12-23 2017-02-28 Игорь Владимирович Рябов Базовая станция дистанционного зондирования атмосферы
RU2655164C2 (ru) * 2016-07-04 2018-05-24 Федеральное государственное бюджетное военное образовательное учреждение высшего образования "Военно-космическая академия имени А.Ф.Можайского" Министерства обороны Российской Федерации Система для определения скорости распространения и направления прихода ионосферного возмущения
FR3061775B1 (fr) * 2017-01-09 2020-03-27 Novimet Procede d'estimation d'un taux de precipitations, notamment de pluie ou de neige
CN107024698A (zh) * 2017-03-16 2017-08-08 河海大学 基于降雨衰减模型的雨区合成孔径雷达成像仿真方法
KR101881418B1 (ko) * 2017-08-23 2018-07-24 영남대학교 산학협력단 기후환경을 고려한 라이다 장치 및 그의 출력조절 방법
CN107703487B (zh) * 2017-09-15 2019-10-08 北京无线电测量研究所 一种一体化的天气雷达双极化组件
RU2693842C1 (ru) * 2018-07-30 2019-07-05 Федеральное Государственное Бюджетное Учреждение "Институт Прикладной Геофизики Имени Академика Е.К. Федорова" (Фгбу "Ипг") Способ зондирования ионосферы и тропосферы
CN109884605B (zh) * 2019-02-02 2020-12-15 中国气象科学研究院 云雨对雷达信号的吸收衰减和米氏散射衰减的提取方法
CN111044984B (zh) * 2020-01-08 2023-07-18 福建省气象科学研究所 一种基于威力分布的雷达性能检测方法
CN111474606B (zh) * 2020-03-30 2022-04-05 中国人民解放军国防科技大学 基于双频双极化微波链路的降雨与水汽综合测量装置及方法
KR102221897B1 (ko) * 2020-08-31 2021-03-04 대한민국 실시간 레이더 반사도 보정오차 감시 방법, 이를 수행하기 위한 기록 매체 및 장치
RU2756977C1 (ru) * 2020-11-22 2021-10-07 Федеральное государственное бюджетное образовательное учреждение высшего образования "Поволжский государственный технологический университет" Радиокомплекс для метеорной и трансионосферной связи
CN113341420B (zh) * 2021-06-07 2023-07-18 中国气象局上海台风研究所(上海市气象科学研究所) 适用不同气候态降雨的双偏振气象雷达雨滴谱反演方法
CN113466868B (zh) * 2021-06-09 2022-03-11 中山大学 一种处理双极化雷达差分相位的方法及系统
KR102632005B1 (ko) * 2021-07-02 2024-02-01 대한민국 다중 지상 통신 기지국 마이크로파를 활용한 엑스밴드 단일 편파 레이더 반사도 보정 장치 및 그 방법
KR102632006B1 (ko) * 2021-07-02 2024-02-01 대한민국 다중 지상 통신 기지국 마이크로파를 활용한 엑스밴드 이중 편파 레이더 반사도 보정 장치 및 그 방법
CN115425389A (zh) * 2022-09-02 2022-12-02 大连海事大学 一种基于差分测距算法的自适应波束扫描天线系统

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS60100676U (ja) * 1984-10-11 1985-07-09 株式会社トキメック レ−ダ装置
JP2635679B2 (ja) * 1988-05-09 1997-07-30 株式会社東芝 気象レーダ装置
JP2660942B2 (ja) * 1991-02-25 1997-10-08 日本無線株式会社 気象観測用レーダにおける降水形態識別方法
FR2742876B1 (fr) * 1995-12-26 1998-02-06 Thomson Csf Procede de determination du taux de precipitation par radar a double polarisation et radar meteorologique le mettant en oeuvre

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107643522A (zh) * 2017-09-19 2018-01-30 中国电子科技集团公司第三十八研究所 一种双频段双极化星载云雨测量雷达系统
CN107643522B (zh) * 2017-09-19 2020-07-07 中国电子科技集团公司第三十八研究所 一种双频段双极化星载云雨测量雷达系统

Also Published As

Publication number Publication date
KR19990076841A (ko) 1999-10-25
DE69604088D1 (de) 1999-10-07
JP2000502794A (ja) 2000-03-07
JP3720848B2 (ja) 2005-11-30
CN1206469A (zh) 1999-01-27
EP0870206B1 (fr) 1999-09-01
EP0870206A1 (fr) 1998-10-14
FR2742876A1 (fr) 1997-06-27
WO1997024629A1 (fr) 1997-07-10
FR2742876B1 (fr) 1998-02-06
US6061013A (en) 2000-05-09
CA2242225A1 (fr) 1997-07-10
CA2242225C (fr) 2004-09-28
DE69604088T2 (de) 2000-01-20
KR100442009B1 (ko) 2004-10-02

Similar Documents

Publication Publication Date Title
CN1114833C (zh) 双极化雷达确定降水概率的方法及采用该方法的气象雷达
CN1078949C (zh) 减轻脉冲多普勒雷达中距离一多普勒模糊的设备和方法
US4851848A (en) Frequency agile synthetic aperture radar
JP3798247B2 (ja) 並列ドップラ処理を行うレーダ式の氷測深器
CN1227539C (zh) 用脉冲雷达测量雷达反射性和多谱勒偏移的方法和系统
Gray et al. Theory of incoherent-scatter measurements using compressed pulses
CN102798856B (zh) 一种基于小波包的24GHz LFMCW雷达测距方法
CN109061648B (zh) 基于频率分集的速度/距离解模糊雷达波形设计方法
CN103675788B (zh) 散射计回波信号中降水回波及海面后向散射的分离方法
CN108445477A (zh) 机场场面异物检测雷达的高精度测距方法
CN109358330B (zh) 一种双基前视sar线性时间同步误差估计方法
CN107942307B (zh) 调频连续波着陆雷达波束的地面入射角估计方法
CN110609263B (zh) 一种同时计算脉冲激光雷达目标回波时延和频偏的方法
DE19934212B4 (de) Verfahren und Vorrichtung zum Messen der Strömungsgeschwindigkeit eines Fluidstromes
CN116047447A (zh) 基于调频连续波雷达的目标距离、速度测算方法及装置
CN103995262B (zh) 用于起伏界面的mimo稀疏阵列的超声测量方法及系统
Zhu et al. Doppler effect and compensation in a Rotating Fanbeam Spaceborne Scatterometer
CN111948621A (zh) 一种极低信噪比下最佳稀疏域变换的lfm雷达信号处理方法
Develet Performance of a synthetic-aperture mapping radar system
Pant et al. Wind profiler and radio acoustic sounding system at IMD, Pune: Some preliminary results
Elhefhawy et al. Study of short-range synthetic aperture radar system
CN113608188B (zh) 一种基于天基tdm-mimo雷达空间碎片的测角方法及系统
Weiner et al. Dual-use waveforms for simultaneous automotive radar and vehicular communication
Yanyan et al. Technique of Doppler compensation for phase-coded signal pulse compression
Hagfors et al. Two-dimensional chirp mapping of delay and Doppler spread targets

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C19 Lapse of patent right due to non-payment of the annual fee
CF01 Termination of patent right due to non-payment of annual fee