CN110018448B - 基于任意取向微波链路的双偏振雷达衰减订正方法 - Google Patents
基于任意取向微波链路的双偏振雷达衰减订正方法 Download PDFInfo
- Publication number
- CN110018448B CN110018448B CN201910253715.0A CN201910253715A CN110018448B CN 110018448 B CN110018448 B CN 110018448B CN 201910253715 A CN201910253715 A CN 201910253715A CN 110018448 B CN110018448 B CN 110018448B
- Authority
- CN
- China
- Prior art keywords
- radar
- attenuation
- attenuation coefficient
- link
- path
- 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.)
- Active
Links
- 238000012937 correction Methods 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000002310 reflectometry Methods 0.000 claims description 19
- 238000001556 precipitation Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 2
- 150000001875 compounds Chemical class 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000010287 polarization Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information 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
本发明公开了一种基于任意取向微波链路的双偏振天气雷达强度衰减订正方法,由微波链路实测的路径衰减计算平均衰减系数;设定衰减系数与比差分相位之比的初始值α,根据天气雷达观测量计算链路路径平均衰减系数;不断调整衰减系数与比差分相位之比α,以雷达计算的链路路径平均衰减系数与微波链路计算的路径平均衰减系数绝对差异为目标函数,当目标函数最小时,α值为最优;采用此最优α值对雷达强度进行衰减订正。本发明能够有效实现双偏振天气雷达强度的衰减订正,增强了丰富的已建微波链路资源的可用性,本发明的适用性和鲁棒性较强。
Description
技术领域
本发明属于大气探测技术,具体为基于任意取向微波链路的双偏振雷达衰减订正方法。
背景技术
天气雷达以其较高的时空分辨率和准确性而成为当前降水监测的常用手段。较短波长(如X波段)天气雷达因具备天线尺寸小、发射功率低、时空分辨率高、机动性好等优势,成为区域降水监测的重要途径。常规天气雷达主要依靠反射率因子与雨强的关系定量估测降水。然而,较短波长雷达受降水衰减影响较大,导致反射率因子降低,所以在使用较短波长雷达反射率因子进行定量降水估测时,必须首先对反射率因子实施衰减订正,以提高降水估测的准确性。
对常规天气雷达来说,常用的衰减订正方法有解析法、迭代法、逐库法及其近似方法等,但由于真实衰减未知以及衰减系数与反射率因子关系的不确定,常会导致订正不稳定问题。针对双偏振雷达的衰减订正,基于典型雷达频率范围内差分传播相位与衰减总量基本呈线性比例关系这一特点,研究者们提出了降水廓线(ZPHI)算法和自一致(SC)方法,此类方法对该线性比例和雷达自身观测误差非常敏感。
利用微波链路能够得到链路雨致衰减的特性,出现了使用微波链路的路径衰减作为参考,实现雷达基数据衰减订正的方法。这些方法主要存在以下问题:一是大多要求所使用的微波链路沿雷达径向,而由于已建设微波链路相对于雷达是任意取向的,故该要求导致大量已建实际微波链路无法应用。二是使用沿雷达径向链路辅助雷达衰减订正,仅能订正链路所处径向的雷达数据,方法的效率不高。三是仅利用了常规雷达参量,没有充分发挥双偏振雷达能够得到偏振参量的优势。
发明内容
本发明的目的在于提出了一种基于任意取向微波链路的双偏振雷达衰减订正方法。
实现本发明的技术解决方案为:一种基于任意取向微波链路的双偏振雷达衰减订正方法,具体步骤为:
步骤1、根据微波链路路径总衰减观测值,计算雷达路径平均衰减系数;
步骤2、设定衰减系数与比差分相位之比α的初始值,根据雷达观测量计算链路路径平均衰减系数;
步骤3、计算目标函数,所述目标函数为雷达路径平均衰减系数和链路路径平均衰减系数的绝对差异;
步骤4、调整衰减系数与比差分相位之比α的值,重复执行步骤2~3,将目标函数最小时的衰减系数与比差分相位之比α作为最优参数值;
步骤5、根据最优α值以及衰减系数计算公式确定链路横跨方位范围内的每个雷达距离库的最优衰减系数AH(r;αopt),实现雷达反射率因子的衰减订正。
本发明与现有技术相比,其显著优点为:本发明将任意取向链路测值引入雷达衰减订正中,能够有效减小衰减对雷达强度带来的误差;本发明订正区域明显增大,由单条径向拓展到链路跨越的方位范围;突破了现有方法中对链路沿雷达径向的限制,使得广泛建设的实际微波链路的可用性大为增加;有效发挥了双偏振雷达得到偏振参量的优势,且本发明适用性增强,便于推广应用。
下面结合附图对本发明做进一步详细的描述。
附图说明
图1是本发明的流程图。
图2是本发明中系数α的优化过程示意图。
图3为本发明效果示意图。
具体实施方式
如图1所示,一种基于任意取向微波链路的双偏振雷达衰减订正方法,结合雷达和微波链路观测值,采用一个迭代过程对天气雷达衰减系数与比差分相位之比α进行最优估计,并使用最优α值进行衰减订正。在实际工程中,如果链路信号频率与雷达电磁波频率不同,将链路频率下的衰减转换为雷达频率下的衰减后,再采用本发明进行雷达衰减订正,具体步骤如下:
步骤1、链路计算雷达路径平均衰减系数:通过比较接收电平和基准电平得到微波链路路径总衰减观测值AttenuML,进而由计算路径平均衰减系数AML(dB/km),如下式:
式中,LML(km)为链路长度,AML为雷达路径平均衰减系数。
在判断微波链路基准电平时,干(晴)湿(雨)时段可以通过分析链路接收信号的功率谱来判断,本方法中,天气雷达和处于微波链路路径上的雨量计也可以作为干/湿时刻判断的指示器。如果谱分析方法的结果、或沿链路路径至少一个雷达距离库、或雨量计显示某时刻有降水,即认定该时刻为降水时刻,则基准电平为最后“干”时刻的接收电平,链路路径衰减即可通过比较接收电平和基准电平得到。
其中:
上式中,r0为第一个有降水的距离库,r1为有降水的最远距离库,r0<r<r1,Z'H(r)为实测的雷达反射率因子,b为雷达距离r处的衰减系数AH(r)[单位:dB/km]和实测的雷达反射率因子Z'H(r)[单位:mm6/m3]之间的关系AH(r)=a[ZH(r)]b中的指数系数,α为衰减系数与比差分相位之比,为距离r1和r2处差分相位的差异。其中,比差分相位是水平和垂直偏振电磁波的传播常数之间的差值,在均匀介质中,可通过两个距离位置的差分相移计算得到。
步骤3、计算目标函数:目标函数δA定义为雷达路径平均衰减系数和链路路径平均衰减系数的绝对差异:
步骤4、迭代确定最优值:如图2所示,根据已有研究结论确定α的取值范围内,设定一个微小步长,不断递增调整α值,重复上述步骤2、3,使δA最小的α值即为最优α;
步骤5、雷达反射率因子衰减订正:根据最优α值和式(3)和(4),计算得到链路横跨方位范围内的每个雷达距离库的最优衰减系数AH(r;αopt),使用下式实现雷达反射率因子的衰减订正:
其中,ZH和Z'H分别为真实(未经衰减)和实测(经衰减)的雷达反射率因子[单位:mm6/m3]。
图3显示了衰减订正前、后X波段双偏振雷达的强度回波以及对应位置处S波段天气雷达的强度分布。其中虚线为衰减订正前的X波段双偏振雷达某径向的强度分布廓线,实线为衰减订正后的径向强度廓线,“+”号线为附近S波段天气雷达在对应位置处的强度廓线。可以看到,订正前,因降水衰减影响,X波段天气雷达反射率因子显著低于长波长S波段雷达反射率因子,远距离处更加明显。而订正后,X波段雷达反射率因子大幅增强,订正量随距离增加而增加。特别是远距离处的降水回波,在订正前强度很弱,而在订正后得以恢复。可见,该订正方法在较大程度上缓解了衰减效应导致的反射率因子误差。
本发明将任意取向微波链路的路径衰减值用于天气雷达被链路横跨区域数据的衰减订正,能够有效实现双偏振天气雷达强度的衰减订正,增强了丰富的已建微波链路资源的可用性,本发明的适用性和鲁棒性较强。
Claims (3)
1.一种基于任意取向微波链路的双偏振雷达衰减订正方法,其特征在于,具体步骤为:
步骤1、根据微波链路路径总衰减观测值,计算雷达路径平均衰减系数,公式为:
式中,AttenuML为微波链路路径总衰减观测值,LML(km)为链路长度;
步骤2、设定衰减系数与比差分相位之比α的初始值,根据雷达观测量计算链路路径平均衰减系数,公式为:
步骤3、计算目标函数,所述目标函数为雷达路径平均衰减系数和链路路径平均衰减系数的绝对差异,具体为:
步骤4、调整衰减系数与比差分相位之比α的值,重复执行步骤2~3,将目标函数最小时的衰减系数与比差分相位之比α作为最优参数值;
步骤5、根据最优α值以及衰减系数计算公式确定链路横跨方位范围内的每个雷达距离库的最优衰减系数AH(r;αopt),实现雷达反射率因子的衰减订正。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910253715.0A CN110018448B (zh) | 2019-03-30 | 2019-03-30 | 基于任意取向微波链路的双偏振雷达衰减订正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910253715.0A CN110018448B (zh) | 2019-03-30 | 2019-03-30 | 基于任意取向微波链路的双偏振雷达衰减订正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110018448A CN110018448A (zh) | 2019-07-16 |
CN110018448B true CN110018448B (zh) | 2023-03-28 |
Family
ID=67190251
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910253715.0A Active CN110018448B (zh) | 2019-03-30 | 2019-03-30 | 基于任意取向微波链路的双偏振雷达衰减订正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110018448B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112213727B (zh) * | 2020-10-15 | 2024-01-02 | 国家卫星气象中心(国家空间天气监测预警中心) | 一种基于主被动微波联合探测的星载雷达的降水订正方法 |
CN113740934B (zh) * | 2021-08-18 | 2022-12-09 | 浙江省大气探测技术保障中心 | 一种基于s波段双偏振天气雷达降水估测方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104656163B (zh) * | 2015-02-04 | 2017-04-05 | 中国人民解放军理工大学 | 一种基于大数据的移动通信网降水分布及动态测量方法 |
CN106547037B (zh) * | 2016-10-21 | 2019-01-01 | 中国人民解放军理工大学 | 双频双偏振微波链路联合的雨滴谱反演及降水类型识别方法 |
CN106324580B (zh) * | 2016-11-02 | 2018-11-06 | 中国人民解放军理工大学 | 一种基于微波链路网的雷达回波衰减订正方法 |
CN107843884B (zh) * | 2017-09-13 | 2021-09-14 | 成都信息工程大学 | 基于双偏振雷达观测提高雷暴天气预警预报准确度的方法 |
CN107942305B (zh) * | 2017-10-11 | 2021-04-09 | 南京大学 | 双偏振雷达***初始差分相位的在线标定方法 |
-
2019
- 2019-03-30 CN CN201910253715.0A patent/CN110018448B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110018448A (zh) | 2019-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106569208B (zh) | 基于杂波图的机场跑道异物检测方法 | |
CN110488393B (zh) | X波段双偏振天气雷达定量测量降水方法及*** | |
EP2278353B1 (en) | Weather radar apparatus and rainfall rate calculation method | |
CN105445701B (zh) | Ddma‑mimo雷达目标的单脉冲角度估计方法 | |
KR100442009B1 (ko) | 이중편광레이다를이용한강수율결정방법및그방법을이용하는기상레이다 | |
CN110018448B (zh) | 基于任意取向微波链路的双偏振雷达衰减订正方法 | |
CN107576963B (zh) | 基于粒子滤波的双偏振雷达差分传播相移的估计方法 | |
CN107942305B (zh) | 双偏振雷达***初始差分相位的在线标定方法 | |
CN104360329A (zh) | 一种全数字阵列相控阵天气雷达的强度定标方法 | |
CN109085555B (zh) | 一种基于多通道特征的假目标抑制方法 | |
CN110749871B (zh) | 双偏振天气雷达的参量估算方法 | |
CN108717181B (zh) | 一种新型室外场静态rcs前置联动标校测试方法 | |
CN110082729B (zh) | 天气雷达补盲平滑处理***的处理方法 | |
CN111751797B (zh) | 一种基于方位角的高频地波雷达一阶和二阶回波谱分界线确定方法 | |
CN110618410B (zh) | 一种用于毫米波雷达的角度校准与计算方法 | |
NL2026449B1 (en) | Abnormal grain condition detection method based on radio tomographic imaging | |
CN116719023A (zh) | Gpr能量聚焦和分辨率联合约束下的介质介电常数估计方法 | |
KR100942689B1 (ko) | 편파파라메터 차등위상을 이용한 편파레이더 강수감쇠 보정방법 | |
CN111065934B (zh) | 利用天线装置测量仰角和/或方位角的方法 | |
CN113835085B (zh) | 一种基于复杂地形补偿的雷达快速测高方法 | |
Bringi et al. | Measurements and inferences of raindrop canting angles | |
CN111541015B (zh) | 一种改善天线角度分辨率的方法及天线 | |
CN113655454A (zh) | 基于毫米波雷达的太赫兹测云雷达反射率因子标校方法 | |
CN118011358B (zh) | 一种双波段雷达雨滴谱反演方法及装置 | |
CN112630732B (zh) | 一种基于isl约束的抗射频干扰设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |