风电

WRF 动力降尺度:复杂地形风电的气象增强

· 南京运梦科技算法团队 · 评审 算法负责人

WRF 动力降尺度:复杂地形风电的气象增强 封面

如果你在山区、丘陵或沿海台地上做风电选址,迟早会撞上同一堵墙:全球再分析和粗网格预报的水平分辨率往往在几十公里量级,而决定一台风机出力的山脊加速、背风涡流、河谷急流,尺度只有几公里甚至几百米。把粗网格的风"硬塞"到复杂地形里,结果常常是山脊低估、河谷高估,整场的年发电量预估随之失真。动力降尺度(dynamical downscaling)正是为弥合这道鸿沟而生——它用一个真正求解大气运动方程的中尺度模式,把粗网格的大尺度环流"翻译"成贴合地形的精细风场。在这条技术路线上,WRF(Weather Research and Forecasting Model)是被使用最广、文献最厚的工具。本文带工程师视角,把它讲清楚,也讲明白它的边界。

关键要点

  • 动力降尺度用 WRF 这类中尺度模式重新求解大气方程,把几十公里的粗网格风场细化到公里乃至百米级,重建插值得不到的山脊加速与谷地通道结构。
  • WRF 是全球使用最广的大气模式之一:按 2017 年 BAMS 综述统计,累计注册超过 36,000 人次、覆盖 162 个国家/地区。
  • 嵌套加密常用约 3:1 的相邻网格比例,可从几十公里逐级降到 1 公里以内;复杂地形配合大涡模拟可一路嵌套到 30 米量级。
  • 降尺度收益高度依赖地形位置:地形开阔的山脊点收益明显,背风坡与河谷站点的负偏差反而可能被放大;分辨率最高也不等于统计分数最优。
  • 水平分辨率落在约 1 公里到 100 米的"灰区"是物理盲区,PBL 参数化与 LES 假设都不成立,结果须谨慎解读。

背景与定位

WRF 不是一篇论文,而是一整套社区数值天气预报系统,权威出处主要有两份。其一是综述性文章 The Weather Research and Forecasting Model: Overview, System Efforts, and Future Directions,作者为 Powers、Klemp、Skamarock 等人,2017 年发表于 Bulletin of the American Meteorological Society(BAMS,美国气象学会会刊),卷 98、第 8 期、第 1717–1737 页。其二是技术细节文档 A Description of the Advanced Research WRF Model Version 4.3,作者为 Skamarock、Klemp、Dudhia 等人,2021 年由美国国家大气研究中心(NCAR)以技术报告 NCAR/TN-556+STR 发布,系统描述了 ARW(Advanced Research WRF)求解器的动力核、物理参数化选项、初始化、边界条件与网格嵌套机制。

为什么说它权威?一是历史与体量:WRF 自 2000 年首次公开发布,按 2017 年那篇 BAMS 综述的统计,全球累计注册已超过 36,000 人次、覆盖 162 个国家/地区,是被引用与使用最广泛的大气模式之一。二是它"既科研又业务"的双重定位——同一套代码既支撑学术研究,也支撑各国气象与能源机构的实际预报。对风电从业者而言,这意味着 WRF 的降尺度方法不是小众实验,而是有海量同行实践与可对标口径的成熟路线。

方法 / 它做了什么

动力降尺度的核心思想,可以用一句话概括:用粗网格数据提供"边界和初始条件",让一个高分辨率的物理模式在小区域内自己把细节算出来。

具体到 WRF,工作流大致是这样的。先取一份粗网格的大尺度场——通常是全球再分析或全球预报——作为模式最外层网格的初始场和侧边界条件(lateral boundary conditions)。然后在内部布置一层或多层逐级加密的嵌套网格(nesting),常见的相邻网格分辨率比例约为 3:1,例如 27 公里 → 9 公里 → 3 公里 → 1 公里逐层细化。每加密一层,模式都会用更精细的地形高程、地表类型与粗糙度数据重新刻画下垫面,并在每个网格点上数值积分大气的动量、热力、水汽等控制方程。

关键在于,降尺度产出的精细风场不是对粗网格的简单插值,而是模式根据真实地形重新求解出来的——山脊处的气流抬升与加速、背风侧的减速与分离、谷地里的通道效应,这些插值永远变不出来的结构,动力模式能够在物理上重建。这也是它相对"统计降尺度/插值"的根本优势:它解的是物理,而不是拟合的相关关系。

对风电最相关的是轮毂高度风场。WRF 内置多种行星边界层(PBL)参数化方案与地表层方案,用于刻画近地面湍流交换,从而给出 10 米、100 米乃至更高层次的风速风向廓线——而 100 米附近正是现代陆上风机轮毂的典型高度区间。

关键结论

  • WRF 是全球使用最广的大气模式之一。 按 2017 年 BAMS 综述统计,累计注册超过 36,000 人次、覆盖 162 个国家/地区,自 2000 年首发以来持续迭代;ARW 求解器的权威技术规格见 2021 年 NCAR/TN-556+STR(版本 4.3)技术报告。这是判断"该路线是否成熟可靠"的硬指标。
  • 嵌套加密能把粗网格细化到公里乃至百米级。 典型相邻网格比例约 3:1,可从几十公里逐级降到 1 公里以内;在复杂地形多尺度试验中(如 2017 年葡萄牙 Perdigão 双山脊野外观测),WRF 配合大涡模拟可一路嵌套到 30 米量级的微尺度网格。需要强调:能加密到的最细网格尺度,决定了所能"解析"的最小湍流涡尺度(即可解析湍流的下限)——网格越细,能直接解析的湍流涡越小。
  • 降尺度对地形暴露良好的山脊点收益明显,对背风/谷地点反而可能恶化。 多项复杂地形评估指出,地形开阔的山脊站点经降尺度后更贴合实测,而处于背风坡或河谷通道里的站点,负偏差可能被进一步放大。换言之,降尺度的增益高度依赖站点的地形位置,不能一概而论。
  • 最高分辨率不等于最优统计分数。 在复杂地形评估中常见一个反直觉现象:分辨率最高的模拟在 MAE、RMSE、偏差这类传统指标上未必最好,但它对局地地形效应的"形态"重建得更真实。对工程应用而言,这意味着选分辨率要看用途——刻画空间分布与回测口径不同。
  • 存在"灰区"(terra incognita)这一物理盲区。 当水平分辨率落在约 1 公里到 100 米之间时,传统一维 PBL 方案"次网格湍流全部参数化"的假设被破坏,而 LES 所需的"湍流主要被解析"也尚未满足——这一过渡带是当前中—微尺度耦合公认的难点,在该口径下风场对湍流方案的选择高度敏感,结果须谨慎解读。

对新能源 / 运梦平台的意义

把上述结论翻译成新能源工程语言,有几条很实在的用法。

风资源评估与选址。 平坦地形上,再分析往往已经够用;但在山区、丘陵、沿海台地这类复杂地形里,动力降尺度能还原山脊加速、谷地通道等关键结构,对机位排布(micro-siting)和单机年发电量预估的价值最大。务必结合上一条结论审视站点:山脊机位可信度高,背风/谷地机位要叠加实测校验,别被"高分辨率"三个字麻痹。

风电功率预测。 降尺度思路同样适用于短期功率预测——用一份业务化的全球预报驱动高分辨率模式,得到贴合场区地形的轮毂高度风,再经功率曲线映射为出力。这条链路对地形复杂、爬坡事件频繁的场站尤其关键。

对待不确定性的纪律。 灰区问题与地形依赖的偏差提醒我们:降尺度不是"分辨率越高越准"的银弹。稳妥做法是把降尺度风场当作物理一致的高质量背景场,再用实测和统计后处理做偏差订正,并以多源对比量化不确定性。运梦气象 API 把多套数据源放在同一接口下,正是为这种"历史回测—预报上线—多源交叉验证"的闭环服务。

在运梦气象 API 上手

要把降尺度的思路落到代码,第一步是拿到对的字段。运梦气象 API 通过 downloadSync 接口,按 dataSourceIdfields 拉取某点、某时间段的变量。

针对复杂地形风电,核心字段是 u100 / v100(100 米高度的纬向、经向风分量)以及派生的 ws / wd(风速、风向)——轮毂高度风速是出力建模的主力变量。若要补充背景气象,可加 tas / sp / hurs(气温、气压、相对湿度)辅助稳定度与空气密度修正。

下面是一个面向复杂地形风资源回测的最小请求示例,历史回测把 dataSourceId 设为 era5

{
  "dataSourceId": "era5",
  "lat": 26.85,
  "lon": 100.23,
  "stime": "2024-01-01 00:00",
  "etime": "2024-12-31 23:00",
  "fields": ["u100", "v100", "ws", "wd", "sp", "tas"],
  "timezone": "8"
}

拿到整年的小时级 100 米风序列后,就可以按本文思路推进:先评估再分析在该地形位置的可信度,对山脊与谷地机位分别核验偏差,再用功率曲线把风速映射为出力并与实测发电量对照。若要把历史回测延伸到上线预报,把 dataSourceId 换成 ger(德国气象局预报源)即可复用同一套字段与代码框架;如需对比 NASA 或 zg1 等其他源做交叉验证,只需替换 dataSourceId

字段口径与单位以官方文档为准:完整字段表见 /docs/weather/data-elements/,接口参数见 /docs/weather/api-reference/,ERA5 数据源说明见 /products/weather-api/era5/,风电功率预测能力见 /products/wind-forecast/

常见问题

WRF 动力降尺度和统计降尺度有什么区别? 动力降尺度用 WRF 这样真正求解大气运动方程的物理模式,根据真实地形重新算出精细风场;统计降尺度只是拟合粗网格与局地的相关关系。前者能重建山脊加速、背风分离、谷地通道这些插值永远变不出来的结构。

WRF 能把分辨率降到多细?复杂地形风电要用多少分辨率? 嵌套加密常用约 3:1 的相邻网格比例,可从几十公里逐级降到 1 公里以内,配合大涡模拟在复杂地形可嵌套到 30 米量级。但分辨率最高不等于统计分数最优,选分辨率要看用途——刻画空间分布与回测口径不同。

为什么降尺度后有些站点反而变差了? 降尺度的增益高度依赖站点的地形位置。地形开阔的山脊机位经降尺度后更贴合实测、可信度高;而处于背风坡或河谷通道里的站点,负偏差可能被进一步放大,这类机位要叠加实测校验。

什么是 WRF 降尺度的"灰区"问题? 当水平分辨率落在约 1 公里到 100 米之间时,传统一维 PBL 方案"次网格湍流全部参数化"的假设被破坏,而 LES 所需的"湍流主要被解析"也尚未满足。这一过渡带是中—微尺度耦合公认的难点,风场对湍流方案高度敏感,结果须谨慎解读。

在运梦气象 API 上做复杂地形风电要拉哪些字段? 核心字段是 100 米高度的 u100 / v100 风分量及派生的 ws / wd(风速、风向),轮毂高度风速是出力建模的主力变量;可加 tas / sp / hurs 辅助稳定度与空气密度修正。历史回测把 dataSourceId 设为 era5,上线预报换成 ger(德国气象局预报源)即可复用同一套字段。

引用与原文

Powers, J. G., Klemp, J. B., Skamarock, W. C., Davis, C. A., Dudhia, J., Gill, D. O., et al. (2017). The Weather Research and Forecasting Model: Overview, System Efforts, and Future Directions. Bulletin of the American Meteorological Society, 98(8), 1717–1737. DOI: 10.1175/BAMS-D-15-00308.1

Skamarock, W. C., Klemp, J. B., Dudhia, J., Gill, D. O., Liu, Z., Berner, J., et al. (2021). A Description of the Advanced Research WRF Model Version 4.3 (NCAR Technical Note NCAR/TN-556+STR). National Center for Atmospheric Research. DOI: 10.5065/1dfh-6p97