一种基于主被动微波联合探测的星载雷达的降水订正方法与流程

本发明涉及卫星雷达系统遥感数据处理领域，特别涉及一种基于主被动微波联合探测的星载雷达的降水订正方法。

背景技术：

海洋和气象卫星上载有的微波雷达系统比如星载微波散射计、微波高度计、降水测量雷达等，能够提供覆盖全球、高分辨率、高频次和长周期的中小尺度表面风场、海浪场、垂直降雨率廓线等观测资料，广泛应用于灾害性天气预报预警、海洋气象相互作用、气候和数值模式资料同化等气象海洋科学研究领域。在轨运行的我国海洋二号(hy-2)卫星上的散射计和高度计已在应用领域发挥了极大价值，我国风云三号(fy-3)卫星上也即将装载散射计和降水测量雷达系统，预期将提供全球气象服务。

微波雷达在进行地气系统探测时，雷达回波受到降水云的影响产生衰减和体散射，为了准确获取目标的散射回波信息，需要对原始回波进行准确的降水检测和降水订正。降水订正是星载雷达观测数据预处理的重要步骤，也是决定雷达定量反演产品精度的关键因素之一。

目前星载微波雷达降水订正方案主要包括基于雷达单独观测和主被动联合观测两种方法。基于雷达单独观测的回波衰减订正方法一般用于降水测量雷达系统，包括通用法和表面反射技术法，基于雷达反射率因子与衰减系数的关系(k-z关系)，计算每层的衰减量，根据全路径衰减等约束条件进行迭代，重构衰减系数，该方法降水衰减的订正精度受到k-z关系、约束条件和方法本身的影响产生较大的计算误差，而且只能用于有大气垂直分辨能力的雷达如降水测量雷达，难以用于散射计这类观测地表的雷达。

基于主被动联合观测的方法一般利用同步观测的微波辐射计，这种方法用在散射计和高度计领域较多，美国seawinds散射计和星载高度计topex/poseidon的业务大气校正算法均利用微波辐射计amsr和tmr的多通道亮温估算ku波段透过率，海洋卫星散射计和高度计的业务算法也是使用相似算法，该校正算法基于提前建立的降雨模型、通过米散射计算雷达频段的降水订正参数(路径总衰减和总体散射)与辐射计透过率、亮温等观测的相关关系，根据被动观测亮温直接对雷达回波进行降水订正，或者基于辐射计反演的降水廓线计算降水订正参数。该方法依赖于多种理想假设精度受到影响，另外最大的问题是只适用于具有相近频段的雷达和辐射计系统，应用于ku频段散射计(13.5ghz)，只能同步搭载相近频段(18.7ghz和37ghz)的微波成像类辐射计才能应用，受卫星载荷规划的约束较大，当由于卫星平台限制无法搭载同步观测的微波成像仪时(例如风云三号e星)，该方法失效。同时，该方法在降水测量雷达上难以应用，有些学者基于结合辐射传输的联合迭代求解方法尝试应用，但校正精度不高。

技术实现要素：

本发明的目的在于克服上述技术缺陷，基于主被动微波联合探测，提出了一种适用于气象卫星和海洋卫星的雷达观测数据的通用降水订正方法，该方法适用于星载散射计、高度计、降水测量雷达和云雷达等雷达系统，当平台上载有同步观测的被动微波辐射计包括微波成像类仪器或微波探测类仪器时均可适用，不再受到卫星平台上载荷规划的限制。

为实现上述目的，本发明提出了一种基于主被动微波联合探测的星载雷达的降水订正方法，所述方法包括：

利用数值预报背景场历史数据和观测历史数据，对降水存在条件下主动雷达观测数据和被动辐射计观测数据的偏差进行订正；

对订正后的主动雷达观测数据、订正后的被动辐射计观测数据和数值预报背景场大气三维空间分布信息进行时空匹配；

以数值预报背景场大气三维空间分布信息作为先验信息，利用订正后的主动雷达观测数据和订正后的被动辐射计观测数据计算得到最优分析大气云雨廓线；

根据最优分析大气云雨廓线计算主动雷达的相关降水订正系数。

作为上述方法的一种改进，所述方法还包括：通过被动辐射计观测数据进行降水检测的步骤，具体包括：

对于微波成像仪器，根据85ghz实际观测与估算的亮温差构建散射指数si85v：

si85v＝f-tb85v

其中，tb19v、tb22v、tb85v为19、22、85ghz通道的垂直极化观测亮温；f为中间参数；a、b、c和d均为系数；当si85v>10k，则标识为降水区；

对于微波探测类仪器，计算窗区通道亮温差siwin：

siwin＝tb89-tb150

其中，tb89为89ghz通道的观测亮温；tb150为150ghz通道的观测亮温；

如果3k<siwin<10k，则标识为层状降水区，如果siwin>10k则标识为对流降水区。

作为上述方法的一种改进，所述利用数值预报背景场历史数据和观测历史数据，对降水存在条件下主动雷达观测数据和被动辐射计观测数据的偏差进行订正；具体包括：

利用数值预报背景场历史数据、辐射传输模式和观测历史数据，选择南北纬60°以内的海洋晴空稳定目标区域，计算全球观测-模拟偏差，按照纬度带进行划分，进行扫描偏差统计，计算扫描偏差订正项；

选用1000～300hpa厚度，200～50hpa厚度，表面温度，水汽总量作为预报因子，根据气团偏差和预报因子间的线性关系，利用最小二乘法，计算气团偏差订正系数；最后根据下式计算偏差订正后的观测值：

其中，j为通道号，为订正后观测，为观测值，为扫描偏差订正项，为气团偏差订正项，xi为预报因子，aij和bj为气团订正系数，和θ分别为观测方位角和观测入射角。

作为上述方法的一种改进，所述数值预报背景场大气三维空间分布信息包括全球温度、湿度、压力垂直廓线、云和降水廓线。

作为上述方法的一种改进，所述对订正后的主动雷达观测数据、订正后的被动辐射计观测数据和数值预报背景场大气三维空间分布信息进行时空匹配；具体包括：

根据主动观测的像元中心位置和时间信息以及被动观测的像元中心位置和时间信息，基于时间和空间阈值将被动观测位置匹配到主动观测位置；

采用线性插值或二次项插值法，将数值预报背景场大气三维空间分布信息经时空插值到主动雷达观测时刻和位置点；

对主动雷达的测量数据，采用线性插值法进行空间垂直方向的插值，到雷达反射率因子垂直逐层位置。

作为上述方法的一种改进，所述以数值预报背景场大气三维空间分布信息作为先验信息，利用订正后的主动雷达观测数据和订正后的被动辐射计观测数据计算得到最优分析大气云雨廓线，具体包括：

对每个观测像元，选择观测位置附近的20×20背景场格点作为窗口，窗口内的数值模式大气廓线作为先验信息；

通过对数据库廓线的求和来计算反演廓线的后验概率分布，积分公式如下：

其中，x^a为最优分析大气云雨廓线，xj为窗口内的大气廓线，y0为多通道主被动观测量组成的矢量；y(xj)为主被动观测算子，即主被动的辐射传输模式；被动辐射传输模式选用逐线模式或快速模式，主动辐射传输模式根据雷达类型自主选取；r为观测误差协方差，通过统计偏差订正后观测-模拟偏差的均方根误差获得，a为归一化因子，n为窗口内的廓线数目。

作为上述方法的一种改进，所述根据最优分析大气云雨廓线计算主动雷达的相关降水订正系数，具体包括：

从最优分析大气云雨廓线中获取云水、云冰、降雪和降雨混合比，经过主动辐射传输过程计算雷达回波的降水订正系数；所述降水订正系数包括逐层衰减系数、路径总衰减和路径总体散射的一个或多个；

利用最优分析大气云雨廓线中的温度廓线和湿度廓线，通过主动辐射传输过程计算得到气体吸收订正系数，完成气体订正。

作为上述方法的一种改进，当星载雷达为散射计时，所述从最优分析大气云雨廓线中获取云水、云冰、降雪和降雨混合比，经过主动辐射传输过程计算雷达回波的降水订正系数；具体包括：

微波散射计ku波段信号受降雨的影响包括大气中的雨水粒子造成雷达信号衰减，以及雨水粒子的体散射使信号增强两种效应；根据辐射传输理论，受降雨影响的后向散射回波σ0′为：

σ0′＝σ0kpath+ηpath

其中，路径总衰减kpath为：

路径总体散射ηpath为：

其中，z为层高，h为总层高，θ为入射角；k(z)为衰减效率，η(z)为后向散射效率；后向散射效率和衰减效率通过滴谱、单个粒子的标准化衰减截面和散射截面计算得出；σ0为晴空海面的后向散射回波；

假设降水粒子滴谱满足m-p分布，根据米散射理论、滴谱假设和降雨混合比廓线，基于上述公式计算得到路径总衰减kpath和路径总体散射ηpath这两项降雨订正系数。

本发明的优势在于：

1、考虑到现有的雷达数据降水订正方法存在精度不高、受平台遥感仪器配置制约等缺点，本发明提出一种气象和海洋卫星雷达系统遥感数据降水订正的通用方法，适用于微波散射计、高度计、降水测量雷达和云雷达等雷达系统，且对平台上载有同步观测的被动微波成像类仪器或微波探测类仪器时均可适用，不再受到卫星平台上载荷规划的限制；

2、本发明的方法利用实时预报的降水背景信息和主被动仪器观测信息，通过贝叶斯后验概率方法精确构建降水客观廓线，通过辐射传输过程计算降水订正系数，能够保证降水订正的精度；

3、本发明的方法能够借助背景场中的物理约束同时对降水、云水和大气温湿度廓线进行反演，因此在降水订正的同时也可以同步进行云水、氧气和水汽的大气订正，这也是其它方法所不具备的优势。

附图说明

图1为本发明的基于主被动微波联合探测的星载雷达的降水订正方法的流程图。

具体实施方式

下面结合附图对本发明的技术方案进行详细说明。

本发明的方法通过引入降水垂直廓线的背景场作为先验信息，以主被动微波联合探测数据作为约束，利用主被动辐射传输理论和贝叶斯方法计算得到最优降水垂直廓线，最后通过主动辐射传输理论计算雷达回波的逐层衰减系数、路径总衰减和路径总体散射等相关订正系数。本发明的方法能够结合降水背景信息和观测信息，通过贝叶斯方法精确构建降水客观廓线，通过辐射传输计算降水订正系数能够保证降水订正精度，为雷达回波的后端反演和应用提供精度更高的输入数据，同时该方法能够借助背景场中的物理约束同时对降水、云水和大气温湿度廓线进行反演，因此在降水订正的同时也可以同步进行云水、氧气和水汽的大气订正，这也是现有技术的其它方法所不具备的。

本发明提出了一种气象和海洋卫星雷达系统遥感数据降水订正的通用方法，适用于多种雷达系统和被动微波仪器；本发明的方法考虑到背景场中云和降水场往往存在空间上的错位问题，在观测位置附近设置窗口，利用窗口内的廓线，计算贝叶斯后验概率，精确构建降水客观廓线，是一种新型的降水反演和订正的思路和方法；

如图1所示，本发明提出了一种基于主被动微波联合探测的星载雷达的降水订正的方法，该方法包括：

步骤1)获取数值预报背景场大气三维空间分布信息和主被动微波观测资料；

由于降水反演需要背景廓线中包含云和降水信息，对背景场数据的选择有特殊要求。如果对历史雷达数据进行重处理的降水订正，或对时效性没有很高的要求，可使用的开源背景场数据包括era-5再分析数据集等，从数据集中提取温压湿廓线、云水成物(云水、云冰、降雪和降雨)廓线、地表参数；如果处理实时数据，背景场数据选择业务运行的全球和区域模式，如中国气象局数值预报中心的grapes模式、wrf中尺度模式实时模拟的结果作为预报的背景场，从数据中提取的云水成物廓线根据模式的微物理方案略有差别。

步骤2)降水检测；

使用被动微波仪器进行降水检测，如果检测为降水，执行降水订正方案，检测为无降水，流程结束。

使用散射指数法识别降水区，对于微波成像仪器，根据85ghz实际观测与估算的亮温差构建散射指数，定义如下：

si85v＝f-tb85v

其中tb19v、tb22v、tb85v为19、22、85ghz通道的垂直极化观测亮温。如果si85v>10k，则标识为降水区。

对于微波探测类仪器，采用窗区si，为窗区通道亮温差(siwin＝tb89-tb150)，其中，tb89为89ghz通道的观测亮温；tb150为150ghz通道的观测亮温；

如果3k<siwin<10k则标识为层状降水区，如果siwin>10k则标识为对流降水区。

步骤3)主被动观测偏差订正；

分为两部分：与卫星扫描角有关的扫描偏差订正和与气团相关的气团偏差订正。

卫星扫描角有关的扫描偏差订正：

利用数值预报背景场历史数据、辐射传输模式和观测历史数据，选择南北纬60°以内的海洋晴空稳定目标区域，计算全球观测-模拟偏差，按照纬度带进行划分，进行扫描偏差统计，计算扫描偏差订正项。

与气团相关的气团偏差订正：

选用1000～300hpa厚度，200～50hpa厚度，表面温度，水汽总量作为预报因子，根据气团偏差和预报因子间的线性关系，利用最小二乘法，计算气团偏差订正系数。最后根据下式计算偏差订正后的观测值：

其中为订正后观测，为扫描偏差订正项，为气团偏差订正项，xi为预报因子，j为通道，aij和bj为气团订正系数，θ分别为观测方位角和入射角。

步骤4)主被动观测和背景场的时空匹配；

由于使用的主被动观测是共平台准同步观测，首先需要根据主被动观测的像元中心位置和时间信息，根据时间和空间阈值将被动观测匹配到主动观测位置；第二步将背景场三维空间的大气和地表信息经时间插值到雷达观测时刻和位置点，时间上可以选择线性插值或二次项插值方法，对降水测量雷达，还需要进行空间垂直方向的插值，到雷达反射率因子垂直逐层位置，方法采用线性插值法。

步骤5)贝叶斯反演；

贝叶斯反演的优势在于对先验信息的利用，利用先验概率密度分布和条件概率密度分布两类信息和贝叶斯公式得到后验概率密度分布，实际在应用中使用先验数据库来进行外部约束，对每个观测像元，选择观测位置附近的20×20背景场格点作为窗口，窗口内的数值模式大气廓线作为先验信息。

通过对数据库廓线的求和来计算反演廓线的后验概率分布，积分公式如下：

x^a为最优分析大气云雨廓线，xj为窗口内的大气廓线，y0为多通道主被动观测量组成的矢量，以微波成像仪和散射计为例，包含微波成像仪十通道观测亮温和散射计两通道雷达后向散射系数。y(xj)为主被动观测算子，即主被动的辐射传输模式，为节省数据处理时间，被动辐射传输模式选用逐线模式、快速模式等，主动辐射传输模式可根据雷达类型自主选取。r为观测误差协方差，可以通过统计偏差订正后观测-模拟偏差的均方根误差获得。a为归一化因子，n为窗口内的廓线数目。

步骤6)降水订正。

根据步骤5)得到的最优分析的大气廓线中的云水、云冰、降雪、降雨混合比廓线，经过主动辐射传输过程计算雷达回波的逐层衰减系数、路径总衰减和路径总体散射等相关订正参数。最后将降水订正参数订正到雷达回波上，完成降水订正。同时可以利用分析廓线中的温度、湿度廓线，通过主动辐射传输过程计算得到气体吸收订正系数，完成气体订正。

以星载微波散射计雷达系统为例，计算的具体过程如下：

微波散射计ku波段信号受降雨的影响包括大气中的雨水粒子造成雷达信号衰减，以及雨水粒子的体散射使信号增强两种效应。根据辐射传输理论，受降雨影响的雷达海面回波可以推导为：

σ0′＝σ0kpath+ηpath

其中衰减项定义为：

体散射项定义为：

σ0为晴空海面后向散射系数，后向散射效率η和衰减效率k可以通过滴谱、单个粒子的标准化衰减截面和散射截面计算得出。z为层高，h为总层高，θ为入射角。假设降水粒子滴谱满足m-p分布，根据米散射理论、滴谱假设、步骤5)得到的降雨混合比廓线，基于以上公式可以计算得到衰减项和体散射项这两项降雨订正系数；同理可计算降雪的订正系数。

最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。