双参云微物理方案伽马雨滴谱函数高精度快速求解方法技术

本发明专利技术涉及一种双参云微物理方案伽马雨滴谱函数高精度快速求解方法，包括：计算M0/M3的比值c；建立μ‑Λ线性相关关系函数；对M0和M3分贝值进行样本分类以构建查找表，通过查表法查出a、b值；通过μ‑Λ线性相关关系和M0/M3的比值构建求解谱形参数的一元三次方程；代入a、b和M0/M3值以求出该谱形参数的解析解；剔除不合理谱形参数解，得到唯一合理实数解；计算得到截距参数N0以构建适用于双参云微物理方案的Γ雨滴谱三参数。采用本发明专利技术所述方法可以快速求解得到高精度的伽马雨滴谱函数三参数，使得雨滴谱分布特征更接近于华南地区雨滴谱观测事实，相比于现有方案计算结果的误差大幅度降低，计算值接近于实际观测结果，实现简单，计算速度快，结果精度高。

A high precision and fast method for solving the gamma raindrop spectrum function of two parameter cloud microphysics scheme

【技术实现步骤摘要】
双参云微物理方案伽马雨滴谱函数高精度快速求解方法
本专利技术涉及数值天气预报
，尤其涉及一种双参云微物理方案伽马雨滴谱函数高精度快速求解方法。
技术介绍
数值天气预报是现代天气预报的主要技术手段和工具，因此发展精细数值天气预报技术是提高强降水预报能力的必经之路。随着计算能力持续快速地提升，数值天气预报模式空间分辨率不断提高，现已到达云尺度/对流尺度(～1km)。此时，传统的积云参数化方案已不再适用，显式处理云降水必须需要使用云微物理方案。数值天气预报模式中云降水微物理过程通过吸收/释放潜热，与模式动力及辐射特性之间存在复杂非线性相互作用，此时云降水微物理过程变得非常重要，并直接决定降水系统的发生发展和演变。因此，在进行高分辨率数值模拟时，云微物理方案模拟精度直接影响模式预报精度。当前全球和区域业务数值天气预报模式对华南强降水的预报能力明显低于平均水平，因此，精确描述云降水微物理过程是提升华南地区强降水预报水平的关键所在。云微物理参数化方案计算效率高，被广泛用于业务精细数值预报模式。数值模式中云微物理参数化方案选用谱分布函数N(D)来描述水凝物的谱分布特征。云微物理参数化方案的计算效率取决于与微物理过程方程相关的PSD阶矩，而与水凝物的物种类本身并没有太大关系。在考虑空气密度和水凝物密度后，水凝物数浓度(Nx)等价于0阶距(M0)，水凝物质量含量(Qx)正比于3阶距(M3)。随着观测技术不断进步，对水凝物谱型认识深入，越来越多研究表明，相比于指数分布，伽马(Gamma，Γ)分布函数中的三个变量：截距参数(N0)、谱形参数μ和斜率参数Λ能更好地描述各类水凝物谱变化情形，且上述变量的描述精度明显高于指数分布函数，适用性更广。当谱形参数μ为正值时，曲线向上弯曲；当谱形参数μ为负值时，曲线向下弯曲；当谱形参数μ为0时，曲线蜕化为指数函数分布。目前常见的云微物理参数化方案主要有单参方案(one-momentschemes)和双参方案(two-momentschemes)两种。单参方案只预报水凝物质量含量(等价于M3)，固定N0和μ的值，Λ的值由水凝物质量含量求解得到，水凝物谱型的演变仅由斜率参数Λ来确定。常见单参方案主要有：Lin(Linetal.,1983)、WSM6(HongandLim,2006)、Goddard(Taoetal.,1989)、SBU-YLin(LinandColle,2011)等。双参方案同时预报水凝物质量含量(等价于M3)和数浓度(等价于M0)，固定μ的值(一般设为0)，谱型参数N0和Λ可以方程联立求解得出，水凝物谱型由N0和Λ两个参数共同决定。常见双参方案主要有：Thompson(Thompsonetal.,2008)、Morrison(Morrisonetal.,2005)、WDM6(LimandHong,2010)、NSSL(Manselletal.,2010)、Milbrandt(MilbrandtandYau,2005)等。由于云降水微物理过程与水凝物谱分布特征有关，谱分布函数参数描述的精度直接影响云微物理方案的模拟精度，因此，水凝物谱分布也是决定数值模式双偏振雷达资料同化观测算子误差的重要因素，尤其是谱形参数，其能够直接影响雷达反射率因子和差分反射率模拟精度。现有单参云微物理方案只能预报截距参数N0，而双参云微物理方案也只能预报截距参数N0和斜率参数Λ，都不能预报谱形参数μ。目前尚无成熟有效预报或诊断谱形参数的云微物理方案。华南地区对流活动频繁，致灾性强降水频发，常导致城市内涝、山洪、山体滑坡等灾害。华南地区大量雨滴谱观测事实表明，谱形参数可正可负，变化范围在-5～20区间，极少为0的情形。而当前单参方案和双参方案通常假设伽马雨滴谱形状参数为常数(对于雨滴一般设为0)，使其模拟水凝物谱分布特征与观测事实严重不符，存在较大误差。这不仅导致现有云微物理方案在模拟华南强降水时有明显不适用性，也使得进行双偏振雷达资料同化时效果退化，是亟待解决的关键技术瓶颈。大量雨滴谱观测事实表明，伽马雨滴谱分布函数三个参数之间并非完全独立。谱形参数μ与斜率参数Λ之间并非独立变量，两者呈现很好的正相关关系。μ-Λ关系主要与云微物理过程有关，但其还可能受气候、地理位置、降水类型等因素影响，导致这种关系随地域和气候差异性的不同而存在明显的差异性。因此，需要基于本地区大样本高精度的雨滴谱观测资料，才能获取本地区具有代表性的μ-Λ关系。μ-Λ相关关系非常重要，它使Γ雨滴谱分布函数由三个独立参数减少为两个。理论上通过雨水的质量含量和数浓度信息，能得到Γ雨滴谱三参数(N0、μ、Λ)。然而现有参考文献中，都采用抛物线函数Λ＝aμ2+bμ+c拟合μ-Λ关系，这使得联立方程时f(μ)成为一个一元六次函数。一元六次方程难有解析解。而采用迭代法计算解时，经常出现来回震荡，难以收敛的情况，且最多存在6个解。基于双参方案预报的雨水质量含量和数浓度信息，采用现有抛物线函数μ-Λ关系，理论上是能求解Γ雨滴谱三参数。一元六次方程f(μ)本身难以有数学解析解。如果采用数值逼近方法，不仅计算量巨大，且容易出现震荡情况，导致不收敛情形。与此同时，当存在多个解时(最多6个)，无法高效剔除无效解和不合理解。因此采用抛物线函数拟合μ-Λ关系，难以构建适用于双参云微物理方案的快速得到高精度伽马雨滴谱分布参数求解方案。
技术实现思路
为此，本专利技术提供一种双参云微物理方案伽马雨滴谱函数高精度快速求解方法，用以克服现有技术中将谱形参数设为固定值，导致计算结果严重偏离观测事实，精度差的问题。为实现上述目的，本专利技术提供一种双参云微物理方案伽马雨滴谱函数高精度快速求解方法，包括：步骤1：进行双参云微物理方案预报，并将预报中的雨水数浓度Nr和物质含量Qr转化为谱分布0阶距参量M0和谱分布3阶距参量M3，分别将M0值和M3值转化为分贝值，并计算得到M0/M3的比值c；步骤2：基于华南地区大量雨滴谱观测得到的伽马雨滴谱分布函数特点，建立谱形参数μ与斜率参数Λ的线性关系函数μ＝aΛ+b；步骤3：对华南地区进行雨滴谱观测，并基于观测事实，采用M0和M3分贝值进行样本分类以构建所述μ-Λ线性相关关系的a、b值查找表，通过查表法将μ-Λ线性相关关系中的a、b值查出；步骤4：通过采用谱形参数μ与斜率参数Λ线性相关关系μ＝aΛ+b和M0/M3的比值c，构建求解谱形参数的一元三次方程f(μ)＝Aμ3+Bμ2+Cμ+D；步骤5：将a、b和M0/M3值代入已构建谱形参数的一元三次方程f(μ)，并求出该谱形参数的实数解析解；步骤6：通过谱形参数μ与斜率参数Λ线性相关关系，计算得到斜率参数值Λ，并剔除不合理解；步骤7：基于华南地区的雨滴谱观测事实，建立华南地区的雷达反射率因子Z与雨强R的相关关系，通过Z-R相关关系以得出唯一合理谱形参数解μ；步骤8：求解后，使用μ、Λ、M0和M3值计算得到截距参数N0以构建适用于双参云微物理方案的Γ雨滴谱三参数(N0、μ、Λ)以描述各类水凝物谱本文档来自技高网...

【技术保护点】
1.一种双参云微物理方案伽马雨滴谱函数高精度快速求解方法，其特征在于，包括：/n步骤1：进行双参云微物理方案数值天气预报，并将预报中的雨水数浓度Nr和物质含量Qr转化为谱分布0阶距参量M0和谱分布3阶距参量M3，分别将M0值和M3值转化为分贝值，并计算得到M0/M3的比值c；/n步骤2：基于华南地区大量雨滴谱观测资料获取的伽马雨滴谱分布函数特点，建立谱形参数μ与斜率参数Λ的线性关系函数μ＝aΛ+b；/n步骤3：对华南地区进行雨滴谱观测，并基于大量观测事实，采用M0和M3分贝值进行样本分类以构建所述μ-Λ线性相关关系的a、b值查找表，通过查表法将μ-Λ线性相关关系中的a、b值查出；/n步骤4：通过采用谱形参数μ与斜率参数Λ线性相关关系μ＝aΛ+b和M0/M3的比值c，构建求解谱形参数的一元三次方程f(μ)＝Aμ

【技术特征摘要】
1.一种双参云微物理方案伽马雨滴谱函数高精度快速求解方法，其特征在于，包括：
步骤1：进行双参云微物理方案数值天气预报，并将预报中的雨水数浓度Nr和物质含量Qr转化为谱分布0阶距参量M0和谱分布3阶距参量M3，分别将M0值和M3值转化为分贝值，并计算得到M0/M3的比值c；
步骤2：基于华南地区大量雨滴谱观测资料获取的伽马雨滴谱分布函数特点，建立谱形参数μ与斜率参数Λ的线性关系函数μ＝aΛ+b；
步骤3：对华南地区进行雨滴谱观测，并基于大量观测事实，采用M0和M3分贝值进行样本分类以构建所述μ-Λ线性相关关系的a、b值查找表，通过查表法将μ-Λ线性相关关系中的a、b值查出；
步骤4：通过采用谱形参数μ与斜率参数Λ线性相关关系μ＝aΛ+b和M0/M3的比值c，构建求解谱形参数的一元三次方程f(μ)＝Aμ3+Bμ2+Cμ+D；
步骤5：将a、b和M0/M3值代入已构建谱形参数的一元三次方程f(μ)，并求出该谱形参数的实数解；
步骤6：通过谱形参数μ与斜率参数Λ线性相关关系，计算得到斜率参数值Λ，并剔除不合理解；
步骤7：基于华南地区的雨滴谱观测事实，建立华南地区的雷达反射率因子Z与雨强R的相关关系，通过Z-R相关关系以得出唯一合理谱形参数解μ；
步骤8：求解后，使用μ、Λ、M0和M3值计算得到截距参数N0以构建适用于双参云微物理方案的Γ雨滴谱三参数(N0、μ、Λ)以描述各类水凝物谱变化情形。


2.根据权利要求1所述的双参云微物理方案伽...

【专利技术属性】
技术研发人员：刘显通，万齐林，胡胜，肖辉，黎慧琦，郑腾飞，夏丰，冯璐，
申请(专利权)人：中国气象局广州热带海洋气象研究所，
类型：发明
国别省市：广东;44