一种LES超高分辨率TC大气流场数值模拟方法及系统

1.一种基于飞艇雷达与下投式探空数据同化驱动的LES超高分辨率TC大气流场数值模拟方法，其特征在于，包括以下步骤：S100. 在目标TC上空预设高度布设飞艇平台，利用搭载的多普勒雷达获取覆盖TC核心区的径向风观测数据，同时释放下投式探空仪，获取压力、温度、湿度及三维风场的垂直剖面观测数据；S200. 对飞艇雷达与下投探空观测数据进行时间同步、质量控制和空间配准，基于飞艇平台位置及雷达波束几何关系构建雷达径向风观测算子，基于探空轨迹构建剖面观测算子，建立观测量与模式变量之间的映射关系；S300. 在具有多重嵌套网格的TC数值模式框架中，采用数据同化方法将经预处理的飞艇雷达与下投探空观测数据与模式背景场融合，得到与嵌套网格分辨率匹配的三维分析场，并据此生成外层与内层网格的边界条件；S400. 在覆盖TC边界层及内核区域的最内层网格上配置LES求解器，将水平网格间距设定为百米量级以下并在垂直方向对边界层进行加密，以构建的三维分析场及其边界条件对LES嵌套域进行初始化；S500. 在LES嵌套域内积分非静力控制方程，采用亚网格尺度湍流参数化方案显式解析主导能量的大涡运动和部分小尺度湍流过程，获得包含边界层涡旋、小尺度对流生成、湍流结构、大涡演变在内的TC三维大气流场；S600. 对LES模拟结果进行后处理，输出TC三维流场、湍流动能、涡旋结构识别、边界层高度变化、对流生成机制分析结果。 2.根据权利要求1所述的方法，其特征在于，步骤S100中，飞艇平台布设于对流层顶附近的平流层下部，并在TC上空处于近似悬停或缓速巡航状态，其高度范围控制在16~22km；根据TC实时位置与最大风半径RMW确定飞艇与TC中心的水平距离DIS，使DIS位于0.8~1.5RMW区间内，确保雷达波束在扫描过程中穿越或覆盖TC最大风速带；同时按照预设扫描策略控制雷达在10°~60°的俯仰角范围内进行体积扫描，兼顾低层边界层结构与中高层内核风场信息的获取。 3.根据权利要求1或2所述的方法，其特征在于，步骤S100中，多普勒雷达在同一扫描周期内采用多仰角、多方位重复体积扫描模式，在雷达扫描期间，根据TC移动速度与飞艇平台相对位置动态调整扫描起始方位与扫描扇区宽度，使雷达观测扇区始终覆盖TC眼墙及近眼区；在释放下投式探空仪时，按照预设的径向和方位布点策略在TC内核及外围雨带不同方位半径上布设多条剖面线路，确保探空剖面在径向和方位上的采样密度。 4.根据权利要求1所述的方法，其特征在于，步骤S200中，在构建雷达径向风观测算子时，显式考虑雷达波束方位角、俯仰角以及探测点地理高度，根据三维风场分量与雷达视线方向之间的几何投影关系，将模式中的切向风、径向风与垂直速度分解到径向观测空间，同时引入降水粒子终端下落速度参数化，并将终端下落速度分量从径向速度中扣除以提高径向风与真实风场的一致性；在构建剖面观测算子时，基于下投式探空仪下落轨迹和时间戳，将模式格点场通过时间插值与空间插值映射到探空测点位置。 5.根据权利要求1所述的方法，其特征在于，步骤S200中，对原始雷达体积扫描数据进行回波品质控制和地物杂波剔除，利用反射率阈值、径向速度退混叠判别及邻域一致性检验去除异常观测；在径向风观测数据同化前，按照预设径向距离间隔和方位角间隔对高密度雷达观测进行空间聚合，形成径向和方位上具有代表性的超级观测，并为每个超级观测估计等效观测误差方差，以降低观测误差空间相关性和数据体量；对下投式探空观测数据则采用垂直平滑与异常值剔除相结合的方式进行预处理。 6.根据权利要求1所述的方法，其特征在于，步骤S300中，数据同化采用集合卡尔曼滤波方法，在外层与内层嵌套网格上分别维持不少于30个集合成员，通过对模式初始场和边界条件施加扰动构建流依赖背景误差协方差；在同化时段内，以雷达径向风观测与下投式探空剖面观测构成联合观测向量，分别利用已构建的雷达径向风观测算子和剖面观测算子对各集合成员进行观测空间投影，获得与各观测对应的模拟观测值，并基于实际观测值与模拟观测值之间的差异对各集合成员的模式状态进行加权修正，得到与嵌套网格分辨率相匹配的三维分析场。 7.根据权利要求1所述的方法，其特征在于，步骤S300中还包括多尺度局地化与权重分配过程，针对不同嵌套网格分别设置不同的协方差局地化半径，使外层网格在100～500km尺度上体现环境场平滑调整，内层网格在不大于100km尺度上突出TC内核结构的精细修正；同时根据雷达径向风与探空剖面在空间代表性与误差特性上的差异分别设置观测误差协方差矩阵，对TC内核环流敏感的雷达径向风赋予高权重，对TC热力层结特征敏感的温度、湿度剖面赋予高权重，在三维分析场中实现动力与热力要素的协同优化。 8.根据权利要求1所述的方法，其特征在于，步骤S400中，在最内层网格中将水平网格间距设定在50~100 m范围内，并在垂直方向上对0~2km的TC边界层区域进行加密布层，使近地面垂直分辨率达到10~40m量级；且在LES嵌套域内关闭行星边界层参数化方案，启用LES求解模式，并选用具有各向异性和非线性回散特性的亚网格尺度湍流参数化方案对未解析小尺度湍流通量进行闭合；同时通过单向或双向嵌套边界条件，将外层中尺度模式的风场、温度场和湿度场以时间插值形式传递给LES嵌套域。 9.根据权利要求1或8所述的方法，其特征在于，步骤S400中，利用同化生成的三维分析场通过三次样条插值方法向LES网格投影；对插值后的初始场进行动力平衡调整，通过数值滤波技术抑制高频重力波噪声，迭代求解非静力平衡方程组至残差收敛；针对边界层内的湍流初始场，基于同化生成的三维分析场的风切变和热力层结特征，叠加符合Kolmogorov谱特性的随机扰动，扰动幅度控制在平均流场的5~10 %范围内。 10.根据权利要求1所述的方法，其特征在于，步骤S500中，在LES嵌套域内对非静力控制方程进行显式时间积分，积分时长至少覆盖TC内核结构稳定发展的若干小时；积分过程中，分别计算解析大涡分量与亚网格尺度通量分量，得到垂直动量通量、湍流动能及其随高度的分布；同时在预定高度上提取垂直速度和涡度的瞬时场及时间平均场，识别并统计边界层条带状涡结构和类龙卷尺度涡旋的空间尺度、强度和出现频次。 11.根据权利要求1或10所述的方法，其特征在于，步骤S500中针对TC边界层物理过程的处理包括：海表通量参数化方案采用考虑海浪影响的COARE算法，根据10米高度风速和海气温差计算感热通量和潜热通量，海表粗糙度长度根据风速和海浪状态动态调整；边界层顶夹卷过程通过显式解析的湍流混合刻画，捕捉边界层顶逆温层的侵蚀和自由大气干空气的下传；辐射传输过程采用RRTMG方案，时间步长设定为5~10min。 12.根据权利要求1所述的方法，其特征在于，步骤S600中，对LES模拟的长时间序列数据进行统计分析，利用时间平均与扰动分解方法获得TC边界层内的平均环流场与湍流扰动场，将解析得到的边界层结构与观测统计结果进行对比，验证LES配置和同化方案的合理性；构建包括最大风速空间分布、极端阵风概率分布、涡旋结构高频出现区域在内的风险指示量场，并以网格化或矢量化形式输出。 13.一种基于飞艇雷达与下投式探空数据同化驱动的LES超高分辨率TC大气流场数值模拟系统，其特征在于，包括：飞艇平台协同观测模块，包括在目标TC上空布设的飞艇平台，控制搭载的多普勒雷达获取覆盖TC核心区的径向风观测数据，并控制下投式探空释放与接收装置获取压力、温度、湿度及三维风场的垂直剖面观测数据；数据预处理与观测算子构建模块，用于对飞艇雷达与下投探空观测数据进行时间同步、质量控制和空间配准，并基于飞艇平台位置、雷达波束几何关系及探空轨迹构建雷达径向风观测算子和剖面观测算子；数据同化与三维分析场构建模块，配置在TC数值模式框架中，用于调用各观测算子，将预处理后的飞艇雷达与下投探空观测数据与模式背景场融合，生成与模式各嵌套网格分辨率相匹配的三维分析场及相应边界条件；LES嵌套域配置模块，用于在覆盖TC边界层及内核区域的最内层网格配置LES求解器，将水平网格间距设定为百米量级以下并在垂直方向对边界层进行加密，并利用三维分析场及边界条件对LES嵌套域进行初始化；LES数值模拟模块，用于在LES嵌套域内积分非静力控制方程，采用亚网格尺度湍流参数化方案显式解析主导能量的大涡运动和部分小尺度湍流过程，获得TC超高分辨率三维大气流场；模拟结果输出模块，用于对LES数值模拟结果进行后处理，输出TC三维流场及湍流统计量、涡旋结构和极端风场诊断结果。