一种基于星载三频毫米波雷达的云雨微物理参数反演方法技术

本发明专利技术公开了一种基于星载三频毫米波雷达的云雨微物理参数反演方法。方法为：首先统计分析云雨粒子时空分布特征并进行归纳，用广义Gamma分布函数表征云雨滴粒子分布，并进行参数化；然后基于贝叶斯估计理论的思想建立反演模型，预设参数的先验分布，将星载W、Ka和Ku三频毫米波雷达反射率因子输入物理模型，输入时空匹配的地基毫米波雷达数据作为调整因子，进行后验，经过迭代计算，不断修正先验概率分布，使代价函数最小，输出最优反演结果；最后对各个距离库内的反演结果进行线性化处理，得到整条廓线内反演结果。本发明专利技术时空分辨率高，成本低，反演结果精细，能够同时反演云雨粒子的实际分布。

An inversion method of cloud and rain microphysical parameters based on space borne three frequency millimeter wave radar

【技术实现步骤摘要】
一种基于星载三频毫米波雷达的云雨微物理参数反演方法
本专利技术属于遥感
，特别是一种基于星载三频毫米波雷达的云雨微物理参数反演方法。
技术介绍
云雨微物理参数的准确测量是气象准确预报的必要保证，是分析各种天气现象发生的关键决策辅助，对于人工影响天气作业也起到至关重要的作用。对云雨的探测主要分为直接测量和遥感探测两种。其中，直接测量方式包括升放无线电探空仪和机载实验测量云雨滴谱等；遥感探测方式是使用各个波段的雷达、激光雷达和辐射计等遥感仪器。其中，投放探空仪和机载实验测量实施物力和人力成本较大，所得结果较为精确，但是时空分辨率较低。用雷达进行探测能够连续不断地获得时空分辨率较高的回波数据，并可以基于相关特定算法来反演云雨物理参数。目前，地基探测主要利用天气雷达，以S波段和C波段雷达为主进行对云雨物理参数的探测。相比较厘米波而言，毫米波雷达波长更接近云等小粒子尺度，且其较窄的波束宽度也更利于探测云和弱降水的结构特征。目前，已经研发出W波段(94GHz)测云雷达和Ka(35GHz)/Ku(13.6GHz)波段测雨雷达，并已投入使用，如目前在轨的CloudSat星载云廓线雷达和GPM星载测雨雷达。无论是测云雷达还是测雨雷达，现有的反演算法的基础都只是针对单种类型云或雨粒子，用简易化的粒子谱函数来拟合分布，并不能精细化地同时反映云雨粒子的实际分布。由于星载W，Ka和Ku三频毫米波雷达同时包括了W波段和Ka/Ku波段，所以在同时获取云雨微物理参数方面具有较大的优势。目前尚无针对星载W，Ka和Ku三频毫米波雷达反演云雨微物理参数的方法，因此，亟需寻找一种星载W，Ka和Ku三频毫米波雷达反演云雨微物理参数的方法。
技术实现思路
本专利技术的目的在于提供一种时空分辨率高、反演结果精细、能够同时反演云雨粒子的实际分布的基于星载W、Ka和Ku三频毫米波雷达的云雨微物理参数反演方法。实现本专利技术目的的技术解决方案为：一种基于星载三频毫米波雷达的云雨微物理参数反演方法，包括以下步骤：步骤1、基于地基毫米波雷达和机载毫米波雷达的云雨探测数据，统计分析云雨滴粒子时空分布特征，拟合典型天气状况下的云雨滴谱函数并进行归纳，典型天气状况包括对流云降水、层状云降水；用广义Gamma分布函数表征云雨滴粒子分布，并进行参数化；步骤2、对星载三频毫米波雷达数据和地基毫米波雷达数据进行时空匹配；步骤3、建立贝叶斯估计反演模型，对广义Gamma分布谱函数中的任意一个谱参数预设先验概率分布，将星载W、Ka和Ku三频毫米波雷达反射率因子输入物理模型，得到其它三个谱参数；步骤4、将地基毫米波雷达的时空匹配数据作为调整因子，输入贝叶斯估计反演模型进行后验，然后进行迭代计算，以修正先验概率分布，使代价函数最小，最后输出最优反演结果；步骤5、进行云相态识别，计算在不同相态中每个距离库的云雨微物理参数，并进行线性化处理，得到整个廓线云雨微物理参数的反演结果；步骤6、分析比对典型天气过程中的机载毫米波雷达云雨探测数据，将步骤4中反演得到的云雨微物理参数与所述云雨探测数据进行比对，计算偏差以验证结果的有效性。进一步地，步骤1中所述的用广义Gamma分布函数表征云雨滴粒子分布，并进行参数化，公式如下：n(r)dr＝arμexp(-brν)dr式中，n(r)为数密度分布函数，r为云雨滴粒子半径；N为粒子总数浓度，因式子中含有Γ函数，故称为广义Γ分布函数；该函数a是截距参数、b是尺度参数、μ是形状参数、ν是功率参数，四个参数均为正数，调节四个参数的值，拟合各种特征谱型的云雨滴谱。进一步地，步骤3所述的建立贝叶斯估计反演模型，对广义Gamma分布谱函数中的任意一个谱参数预设先验概率分布，将星载W、Ka和Ku三频毫米波雷达反射率因子输入物理模型，得到其它三个谱参数，具体如下：步骤3.1、建立贝叶斯估计反演模型为：式中，X为a、b、μ和ν四个参数中的一个参数，设定为已知量，服从预设的先验概率分布；Y为测量向量，是在X服从先验概率分布条件下，输入星载W、Ka和Ku三频毫米波雷达的三个频段雷达反射率因子；Ppr表示预设的先验概率密度函数，Pf表示条件概率密度函数；步骤3.2、联立方程组求得其余三个谱参数：式中Ze1、Ze2和Ze3分别为星载W、Ka和Ku三频毫米波雷达的三个频段经过衰减订正的雷达反射率因子值，λ1、λ2和λ3分别为星载W、Ka和Ku三频毫米波雷达的三个频段工作波长，σb1、σb2和σb3分别为星载W、Ka和Ku三频毫米波雷达的三个频段测得的后向散射截面；K为与复折射指数m有关的常数，D为粒子直径。进一步地，步骤4中所述的进行迭代计算，以修正先验概率分布，使代价函数最小，具体如下：迭代公式为：式中，n为设定的廓线内的云距离库个数，x1和xn分别为在云低和云高距离库内的先验状态参数，a1和an、b1和bn、ν1和νn分别为在先验条件下云低和云高距离库内求得的其余三个谱参数；反演算法在代价函数Φ取得最小值的条件下求得最优解，Φ为与X差分以及与Y差分的权重之和，其中，为迭代计算后的状态参数，为迭代计算后的物理模型输出的谱参数，Sx为先验状态参数向量的协方差矩阵，为物理模型输出的谱参数误差协方差矩阵。进一步地，步骤5中所述的进行线性化处理，公式如下：式中，同名像点α为线性化因子，T为温度，m为常数，根据云层相态信息进行取值，N1,N2...Nn分别为各个有云距离库内经过反演模型得到的云雨微物理参数，包括粒子数浓度、水含量、水路径、有效粒子半径、消光系数和雨强。本专利技术与现有技术相比，其显著优点在于：(1)使用雷达不间断获得回波数据，时空分辨率高；(2)经过迭代计算不断修正先验概率分布，使代价函数最小，反演结果精细，能够同时反演云雨粒子的实际分布；(3)利用获得的W，Ka和Ku三频的雷达回波数据，基于相关算法来反演云雨微物理参数，成本低。附图说明图1是本专利技术基于星载三频毫米波雷达的云雨微物理参数反演方法的流程示意图。图2是本专利技术实施例中云雨微物理参数的反演结果曲线图。具体实施方式本专利技术基于星载三频毫米波雷达的云雨微物理参数反演方法，包括以下步骤：步骤1、基于地基毫米波雷达和机载毫米波雷达的云雨探测数据，统计分析云雨滴粒子时空分布特征，拟合典型天气状况下的云雨滴谱函数并进行归纳，典型天气状况包括对流云降水、层状云降水；用广义Gamma分布函数表征云雨滴粒子分布，并进行参数化；步骤2、对星载三频毫米波雷达数据和地基毫米波雷达数据进行时空匹配；步骤3、建立贝叶斯估计反演模型，对广义Gamma分布谱函数中的任意一个谱参数预设先验概率分布，将星载W、Ka和Ku三频毫米波雷达反射率因子输入物理模型，得到其它三个谱参数；步骤4、将地基毫米波雷达的时空匹配数据作为调整因子，输入贝叶斯估计反演模型进行后验，然后进行迭代计算，以修正先验概率分布，使代价函数最小，最后输出最优反演结果；步骤5、进行云相态识别，计算在不同相态中每个距离库的云雨微物理参数，并进行线性化处理，得到整个廓线云雨微物理参数的反演结果；步骤6、分析比对典型天气过程中的机载毫米波雷达云雨探测数据，将步骤4中反演得到的云雨微物理参数与所述云雨探测数据进行比对，计算偏差以验证结果的有效性。进一步地，步骤1中所述的用广义Gamma分本文档来自技高网...

【技术保护点】
1.一种基于星载三频毫米波雷达的云雨微物理参数反演方法，其特征在于，包括以下步骤：步骤1、基于地基毫米波雷达和机载毫米波雷达的云雨探测数据，统计分析云雨滴粒子时空分布特征，拟合典型天气状况下的云雨滴谱函数并进行归纳，典型天气状况包括对流云降水、层状云降水；用广义Gamma分布函数表征云雨滴粒子分布，并进行参数化；步骤2、对星载三频毫米波雷达数据和地基毫米波雷达数据进行时空匹配；步骤3、建立贝叶斯估计反演模型，对广义Gamma分布谱函数中的任意一个谱参数预设先验概率分布，将星载W、Ka和Ku三频毫米波雷达反射率因子输入物理模型，得到其它三个谱参数；步骤4、将地基毫米波雷达的时空匹配数据作为调整因子，输入贝叶斯估计反演模型进行后验，然后进行迭代计算，以修正先验概率分布，使代价函数最小，最后输出最优反演结果；步骤5、进行云相态识别，计算在不同相态中每个距离库的云雨微物理参数，并进行线性化处理，得到整个廓线云雨微物理参数的反演结果；步骤6、分析比对典型天气过程中的机载毫米波雷达云雨探测数据，将步骤4中反演得到的云雨微物理参数与所述云雨探测数据进行比对，计算偏差以验证结果的有效性。

【技术特征摘要】
1.一种基于星载三频毫米波雷达的云雨微物理参数反演方法，其特征在于，包括以下步骤：步骤1、基于地基毫米波雷达和机载毫米波雷达的云雨探测数据，统计分析云雨滴粒子时空分布特征，拟合典型天气状况下的云雨滴谱函数并进行归纳，典型天气状况包括对流云降水、层状云降水；用广义Gamma分布函数表征云雨滴粒子分布，并进行参数化；步骤2、对星载三频毫米波雷达数据和地基毫米波雷达数据进行时空匹配；步骤3、建立贝叶斯估计反演模型，对广义Gamma分布谱函数中的任意一个谱参数预设先验概率分布，将星载W、Ka和Ku三频毫米波雷达反射率因子输入物理模型，得到其它三个谱参数；步骤4、将地基毫米波雷达的时空匹配数据作为调整因子，输入贝叶斯估计反演模型进行后验，然后进行迭代计算，以修正先验概率分布，使代价函数最小，最后输出最优反演结果；步骤5、进行云相态识别，计算在不同相态中每个距离库的云雨微物理参数，并进行线性化处理，得到整个廓线云雨微物理参数的反演结果；步骤6、分析比对典型天气过程中的机载毫米波雷达云雨探测数据，将步骤4中反演得到的云雨微物理参数与所述云雨探测数据进行比对，计算偏差以验证结果的有效性。2.根据权利要求1所述的基于星载三频毫米波雷达的云雨微物理参数反演方法，其特征在于，步骤1中所述的用广义Gamma分布函数表征云雨滴粒子分布，并进行参数化，公式如下：n(r)dr＝arμexp(-brν)dr式中，n(r)为数密度分布函数，r为云雨滴粒子半径；N为粒子总数浓度，因式子中含有Γ函数，故称为广义Γ分布函数；该函数a是截距参数、b是尺度参数、μ是形状参数、ν是功率参数，四个参数均为正数，调节四个参数的值，拟合各种特征谱型的云雨滴谱。3.根据权利要求1所述的基于星载三频毫米波雷达的云雨微物理参数反演方法，其特征在于，步骤3所述的建立贝叶斯估计反演模型，对广义Gamma分布谱函数中的任意一个谱参数预设先验概率分布，将星载W、Ka和...

【专利技术属性】
技术研发人员：胡雄，余茁夫，谭仲辉，胡申森，严卫，
申请(专利权)人：中国人民解放军国防科技大学，
类型：发明
国别省市：江苏,32