技术特征：
1.一种顾及三维风场影响的雷达降雨估算误差校正方法，其特征在于，所述方法包括：获取雷达监测的降雨强度数据以及wrf模型的三维风场数据；对所述降雨强度数据进行分析，确定各雨滴的直径和初始位置；根据各雨滴的直径、初始位置以及三维风场数据，基于地面与雷达波之间的空间划分的多个三维子空间，模拟地面与雷达波之间的各三维子空间内的质量加权雨滴运动轨迹，直至模拟到各所述雨滴到达地面，获得在风漂移作用下各所述雨滴到达地面的最终实际位置；根据各所述雨滴到达地面的最终实际位置，获得雷达与雨量计之间的空间关系；根据所述空间关系下雷达和雨量计对的降雨值，计算雷达偏差的校正系数，对雷达降雨值进行校正。2.根据权利要求1所述的顾及三维风场影响的雷达降雨估算误差校正方法，其特征在于，所述对所述降雨强度数据进行分析，确定各雨滴的直径和初始位置的步骤，包括：按照wrf模型的水平分辨率将降雨强度数据划分成水平方网格，即若干个雷达像元；再根据wrf模型的垂直分辨率划分地面与雷达波之间的空间为若干个垂直层，最终地面与雷达波之间的空间体积被划分为若干三维子空间；根据降雨强度和雨滴直径的关系，得到各雷达像元相应雨滴的直径；以雷达像元的中心点为初始水平坐标（x，y），以该中心点对应的雷达波束高度为初始高度（），设置雨滴的初始位置（ x，y，）；式中：h
c
为该中心点对应的雷达波束高度；r
d
为该中心点到雷达测站的距离，r为地球半径；rea指雷达仰角；ir为折射率。3.根据权利要求1所述的顾及三维风场影响的雷达降雨估算误差校正方法，其特征在于，所述根据各雨滴的直径、初始位置以及三维风场数据，基于地面与雷达波之间的空间划分的多个三维子空间，模拟地面与雷达波之间的各三维子空间内的质量加权雨滴运动轨迹，直至模拟到各所述雨滴到达地面，获得在风漂移作用下各所述雨滴到达地面的最终实际位置的步骤，包括：根据雨滴的直径、初始位置以及三维风场数据，确定雨滴的位移关系，模拟地面与雷达波空间的每层雨滴的运动轨迹，所述运动轨迹通过雨滴的位移关系确定，雨滴的位移关系表达式为；
式中，x、y和是三维风向的三个相互正交的方向；u、v和w分别是风场的x、y和分量；m为雨滴质量，g为重力加速度；re为雷诺数；为空气密度；为水密度；为空气粘度；d为雨滴直径；c
d
为雨滴的阻力系数；在每个时间步长内进行轨迹模拟，在该时间步长的最后，得到雨滴的新位置和速度将用于下一个时间步长的轨迹模拟；将所有时间步长的x、y和方向的位移分别累加，反复迭代进行雨滴运动轨迹模拟，直到雨滴移出当前子空间，所述所有时间步长的x、y和方向的位移分别累加表示为：式中：dx
i,k
、dy
i,k
和是第i网格的第k层所有时间步长在x、y和方向上的位移累加量；dx
i,k,t
、dy
i,k,t
以及分别为第i网格的第k层以及时间步长t在x、y、方向上的位移；tn
k
指第k层的时间步长；将边界处的雨滴出口位置和速度作为下一子空间的输入；进入下一个子空间继续反复迭代进行雨滴运动轨迹模拟，直至模拟到所述雨滴到达地面，获得在风漂移作用下各所述雨滴到达地面的最终实际位置的坐标（x
i,r
，y
i,r
），该最终实际位置定义为雷达像元的地面映射点，其中，最终实际位置的坐标为：式中，（x
i,r
，y
i,r
）为网格i的坐标；（x
i,o
，y
i,o
）为雨滴的原始水平坐标；kn指层数。4.根据权利要求3所述的顾及三维风场影响的雷达降雨估算误差校正方法，其特征在于，根据各所述雨滴到达地面的最终实际位置，获得雷达与雨量计之间的空间关系的步骤，包括：针对每个雨量计，搜索每个雨量计周围的雷达像元的地面映射点，并从中选择空间最
邻近该雨量计的一个映射点，该映射点对应的雷达像元与该雨量计被视为一对经过校正的雷达-雨量计对，实现雷达-雨量计对的空间匹配。5.根据权利要求1所述的顾及三维风场影响的雷达降雨估算误差校正方法，其特征在于，根据所述空间关系下雷达和雨量计对的降雨值，计算雷达偏差的校正系数，对雷达降雨值进行校正的步骤，包括：根据所述空间关系下的雨量计的降雨值和雷达像元的降雨值计算雷达偏差的校正系数；式中，f是雷达偏差校正系数；是获得正确空间匹配的雷达-雨量计对的函数；g是雨量计的降雨值；r是雷达像元的降雨值；根据所述校正系数，校正雷达降雨值为：式中，是经过校正后的雷达降雨值；r是雷达像元的降雨值；f是雷达偏差校正系数。6.根据权利要求5所述的顾及三维风场影响的雷达降雨估算误差校正方法，其特征在于，所述雨量计的降雨值为校正后的雨量计雨量；获得所述校正后的雨量计雨量的方式为：根据雨滴运动轨迹的模拟和三维风场数据，计算雨滴和风速直径相应的风湍流误差；式中，w为风速，d为雨滴直径；e
num
为风湍流误差；、和是仅和风速相关的系数；通过所述风湍流误差在雨滴谱上积分，计算由于风效应而漏失的降雨量；式中：表示由于风效应而漏失的总降雨量；d
m
为雨滴质量加权直径；v为雨滴的最终下落速度；n(d，d
m
)为雨滴谱；校正后的雨量计雨量为雨量计初始测量值与漏失的降雨量之和。
技术总结
本申请涉及一种顾及三维风场影响的雷达降雨估算误差校正方法。该方法包括：获取雷达监测的降雨强度数据以及WRF模型的三维风场数据；确定各雨滴的直径和初始位置；根据各雨滴的直径、初始位置和三维风场数据，基于地面与雷达波之间的空间划分的多个三维子空间，模拟地面与雷达波之间的雨滴运动轨迹，获得在风漂移作用下各雨滴到达地面的最终实际位置；根据各雨滴到达地面的最终实际位置，获得雷达与雨量计之间的空间关系；根据该空间关系下雷达和雨量计对的降雨值，计算雷达偏差的校正系数对雷达降雨值进行校正，提高了雷达降雨估算的准确性。确性。确性。
技术研发人员：李文慧 戴强 杨琪祺 吴薇
受保护的技术使用者：南京师范大学
技术研发日：2022.04.02
技术公布日：2022/5/16