观潮 · 模式预报与气象 AI 观察
飓风 Laura 过境时的风机载荷：55.55 米网格算出的塔基弯矩系统性偏小，10 米网格才够
> 这篇是「Part I」。作者用 WRF 做中尺度、再用三档网格的大涡模拟（粗 55.55 米、中 25 米、细 10 米）把 2020 年飓风 Laura 过境时的边界层算出来，直接当成气动弹性模型的入流，驱动 IEA Wind 15 MW 参考风机在空转状态下的响应。结论对做台风区风资源评估的人很直接：粗网格把塔基弯矩系统性算小，细网格与大涡模拟能分辨的小尺度湍流才是关键；中间那档 25 米网格能把塔筒载荷估准，却估不准叶片。作者自己划的边界也很明确——只算了眼墙正面穿过风机位置这一个情形，只算了一个台风，且 Part II 才会做不同位置与设计标准的对比。
缺的是哪一块
风机已经装进了越来越多受热带气旋影响的地区，而有些地方热带气旋的重现期比风电场 25 年的寿命还短，灾害风险需要重新审视。现实中的破坏并不罕见：2003 年台风 Maemi 过境日本宫古岛，6 台风机里 3 台倒塌、另外 3 台叶片严重受损；2017 年飓风 Maria 东北象限扫过波多黎各 Punta Lima 风电场，13 台风机塔筒和叶片损毁，而离路径更远的另一个风电场毫发无损；中国的 Ragasa 和 Yagi 登陆时也造成过风机破坏。
问题出在设计规范用的风模型上。ASCE、API、IEC 这些标准给极端风廓线的方式不同（API 用对数律，ASCE 和 IEC 用幂律），但假设相同：近中性层结、只由地面切应力决定。它们用的湍流模型也不代表台风条件——Kaimal 谱基于平坦地形观测，假定单向、只随高度变化、平稳且水平均匀；台风边界层低层却是水平和垂直切变都很大、湍涡高度各向异性。Mann 均匀切变模型引入了均匀垂直切变，但仍处理不了水平切变和非零的平均垂直、横流分量，也可能抓不到强风暴里的低频能量。
缺的是台风条件下、风机轮毂高度上可信的平均风与湍流观测。实测在近地面和眼墙附近极其有限，所以只能靠高分辨率数值模拟补。
三档网格怎么搭起来的
中尺度部分用 WRF 两重嵌套：外层做谱逼近，内层随涡移动以提高风暴移动过程中的刻画；初始和边界强迫用 NCEP GFS 分析和业务海表温度海冰分析。垂直网格拉伸。
微尺度用离线耦合框架 MOSAIC：把随涡移动域每 30 秒一次的三维输出（2020 年 8 月 26 日 22:00 UTC 到 27 日 04:00 UTC）作为初边条件。粗档大涡模拟覆盖整个眼墙和远离中心的流场，用 WRF 跑，水平网格 55.55 米；中档和细档用 GPU 常驻的 FastEddy 跑，分别是 25 米和 10 米，时段为 26 日 22:30 UTC 到 27 日 04:00 UTC。物理过程上，粗档与中尺度一致（WSM6 微物理、RRTM 长波、Dudhia 短波），FastEddy 只做水汽与云水的转换、不含辐射；两者都用 Monin–Obukhov 相似理论算地表通量，动量交换系数按 Donelan 封顶，热和水汽通量按 Garratt 参数化。
档位 · 水平网格 · 模式 · 作用
粗 · 55.55 m · WRF · 覆盖整个眼墙与远场，提供初边条件
中 · 25 m · FastEddy · 分辨更小尺度湍流
细 · 10 m · FastEddy · 参考标准
湍流起转也被单独处理：粗档的一维能谱在初始化后 25 分钟收敛，因此前 30 分钟作为起转丢弃；细档里小尺度湍流在进入域内若干公里后才发展起来，而入流边界随风暴北移不断变化，所以作者把水平域取得足够大，保证从任一侧边界向内都有充分的发展距离。
风机侧用 OpenFAST 做 IEA Wind 15 MW 参考风机在单桩上的空转响应，启用 ElastoDyn、BeamDyn、AeroDyn、ServoDyn、InflowWind；叶片顺桨、变桨控制关闭、发电机在模拟时长内不启动，并假定有备用电源使偏航仍可工作。
飓风过境时风场怎么变
平均结构随时间变化很大：轮毂高度风速在风暴靠近时上升，进入风眼时下降，眼墙后缘过境时再上升，之后缓慢减小。由于大尺度结构不对称，眼墙前缘和后缘的平均风并不相同——最快的风出现在眼墙前缘，其次是后缘。
但垂直结构的变化要小得多。眼墙前缘靠近时，风速切变基本不变，风向切变（北半球常见的风随高度顺转）倾向于增大；进入风眼后风很静，幂律廓线不再能刻画。眼墙后缘过境后，风速切变和风向切变同时增大。
网格分辨率带来的差别主要在近地面：粗档的幂律指数平均是细档的两倍（在轮毂高度风速超过某阈值时），中档只略大一点；粗档和中档的风向顺转也平均更大。有意思的是，三档算出来的轮毂高度平均风速几乎一样——差别全在廓线形状和湍流上。
载荷结果：粗网格系统性偏小
塔基和叶根弯矩在眼墙附近最大、在风眼内最小，最大平均与瞬时载荷都出现在眼墙前缘过境时。一个值得注意的细节是：结构上的最大弯矩出现在平均风速峰值之前 2 分钟；眼墙后缘过境时的最大载荷则出现在平均轮毂高度风速峰值之后 1 分钟。
用粗档入流时，塔基弯矩在多数时刻被低估；靠近眼墙处，粗档给出的塔基和叶根载荷分布的分位数比细档小 20% 到 30%，且所有情形下这些分位数差异都统计显著（用非参数 bootstrap 判断置信区间是否重叠）。载荷分布两者都是右偏的，但细档入流的右尾更长，极端载荷更多。
中档的表现是分裂的：塔基载荷在眼墙过境时能估准，与细档的分位数彼此落在很小范围内；叶根载荷则在整个过境过程中都没估好，与细档的中位数和极端值差距明显。
为什么：激不激起结构自振频率
原因在频率域。载荷的功率谱显示，用更高分辨率入流时，高频部分的方差更大。粗档入流分辨不了小尺度湍流（对应更高频率），因此在这些频率上塔基和叶根弯矩的方差被严重低估。
更关键的是结构自振频率附近的窄带方差。粗档入流在塔筒和叶片自振频率附近几乎没有能量，激不起这些模态，于是整个过境期间的载荷都被低估；中档入流有足够的能量激起塔筒的自振频率（与细档的差异平均很小），但在叶片自振频率附近能量仍然不足，所以叶片载荷偏小。作者由此给出一个设计规则：网格要细到能分辨激起结构自振频率的那些小尺度湍流结构——给定结构自振频率、轮毂高度特征风速和模式的等效分辨率，可以反推出允许的最大网格间距。
这里还有个反直觉的结果：粗档入流在轮毂层上的风速切变和风向切变更大，而前人工作显示风向顺转会放大空转风机在飓风条件下的载荷，但粗档并没有因此给出更大的平均和最大载荷。作者的解读是，小尺度湍流的能量、它们在时空上的相关性，以及它们产生的不同平均廓线，才是空转风机载荷估计的关键。
原文图
图：中尺度到微尺度的嵌套域配置与 Laura 的近地面风速场。(a) 中尺度两重嵌套域与粗档大涡模拟域，(b) 粗档域内放大后的眼墙结构，观测路径与模拟路径逐点对照，(c) 细档大涡模拟域及位于域中心、与模拟中点风眼重合的虚拟风机位置。加密区就设在风暴中途经过的位置。
图：虚拟风机位置上平均风廓线与轮毂高度风速随风暴过境的变化。(a) 顺流向平均风速廓线，(b) 横流向分量廓线，两者均以 150 米轮毂高度为参考画出四个代表时刻，线色与 (c) 中的竖线对应；(c) 是三档网格算出的轮毂高度平均风速时间序列，眼墙前缘超过 70 米每秒的峰值、进入风眼时的骤降、以及眼墙后缘的次峰清晰可辨，三档曲线几乎重合——差别不在平均风速，而在廓线形状和湍流。
图：轮毂层上风速切变指数 α（上）与风向切变 Δφ（下）随时间的变化，四条竖线对应四个代表时刻。眼墙过境时段（约 2.5 到 3 小时）切变剧烈波动甚至变号，其余时段粗档曲线整体高于细档——幂律指数和风向顺转都偏大，差别主要来自轮毂高度以下的近地面风。
作者承认的限制
- 只算了眼墙正面穿过风机位置这一个情形。 作者明说热带气旋恰好从一台风机上方经过的可能性很低，不同距离和相对路径的响应放到 Part II。
- 只算了一个热带气旋。 相同强度的不同风暴结构可以不同，强度本身也会改变载荷，作者建议覆盖更广的风暴样本，并且随着分析更多风暴，网格敏感性也要一并重新检验。
- 粗网格未必永远不够。 对高强度热带气旋的眼墙，较粗的网格可能已经够用；但更弱的风暴或离眼墙更远的位置，特征风速更慢，可能需要更细的网格——作者明确说 10 米这一档也未必定够。
- 没算水动力载荷。 模拟聚焦对风条件的敏感性，单桩上的水动力载荷未纳入，也不使用 SubDyn，而是通过修改 ElastoDyn 输入让固支塔的自振频率和模态形状与单桩支撑的塔匹配。
- 两档微尺度的物理参数化并不完全一致。 WRF 与 FastEddy 在微物理和辐射上处理不同；作者的说法是小域尺度加上高频边界强迫保证了风暴结构在各模拟间一致。
- 合成湍流模型这条路线也有局限。 作者引用前人指出，合成湍流虽然常用，但可能代表不了实际能观测到的谱的多样性，需要仔细调参。
一句话判断
对台风区风机载荷分析，入流网格能细到什么程度不是精度偏好问题，而是能不能激起结构自振频率的物理门槛：25 米能估准塔筒、估不准叶片，要算叶片得往 10 米走。