风电

欧洲风图集与 WAsP 方法 — 风资源评估的奠基之作

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

欧洲风图集与 WAsP 方法 — 风资源评估的奠基之作 封面

做过风电前期的人都遇到过同一个尴尬:测风塔立在 A 点,机位却要排到几公里外的 B 点;A 点的 50 米实测风速曲线漂亮,可 B 点的山脊、林带、村落和坡度完全不一样,直接把 A 点的风速搬过去显然不靠谱。更现实的问题是——一座测风塔只能覆盖一小片地形,而一个风电场动辄几十平方公里,靠"加塔"把每个机位都测一遍既不经济也不现实。

那么,凭一两座塔的实测数据,如何科学地外推到整个场区的每一个机位?这正是 1989 年《欧洲风图集》(European Wind Atlas)要回答的核心问题。这本由丹麦 Risø 国家实验室的 Ib Troen 与 Erik Lundtang Petersen 撰写的报告,提出了一套被称为"风图集方法"(wind atlas methodology)的思路,并把它固化进了 WAsP(Wind Atlas Analysis and Application Program)软件。三十多年过去,今天几乎所有主流风资源评估软件的产能计算内核,都能追溯到这套方法学。

本文梳理《欧洲风图集》的物理框架与工程逻辑,并说明在现代实践里,如何用运梦气象 API 的 ERA5 长序列(dataSourceId=era5)作为"长期参考风气候"的输入,让这套经典方法在数据可得性上重新焕发活力。

关键要点

  • 《欧洲风图集》由丹麦 Risø 国家实验室的 Troen 与 Petersen 于 1989 年完成,奠定了沿用至今的 WAsP 风资源评估方法学。
  • 核心是"风图集方法":把站点实测风做"去地形化"(剥离当地粗糙度、地形与障碍物影响),还原为"广义风气候"(generalized wind climate),再应用到目标点(重新引入目标点的地形效应),实现风资源的区域外推。
  • 关键物理链条包括地转拖曳定律(geostrophic drag law,连接高空地转风与近地面风)、粗糙度长度 z0 及粗糙度变化、地形(orography)抬升、障碍物尾流。
  • 输出是各风向扇区的威布尔分布参数与风功率密度,直接服务于选址与产能估算。
  • 局限在于:对复杂地形(山地)和非中性大气采用了线性化假设,存在误差,后续业界发展了 CFD 类方法做补强。

背景与定位

《欧洲风图集》并非一篇期刊论文,而是一份体量厚重的技术报告/专著(ISBN 87-550-1482-8),由欧洲共同体委托、Risø 国家实验室在丹麦罗斯基勒(Roskilde)牵头完成。它的两位作者 Ib Troen 与 Erik Lundtang Petersen,都是大气边界层与风能领域的资深研究者。这份工作把分散在欧洲各地气象站的长期观测,整理成可供工程使用的区域风资源图,本身就是一项浩大的数据与方法学工程。

它为什么重要?因为在此之前,风资源评估很大程度上依赖经验外推和局部测风,缺乏一套既有物理基础、又能在工程上落地的统一框架。《欧洲风图集》提供的不只是"地图",更是"如何从一个测点推算到另一个点"的可计算方法——这正是它真正的遗产。WAsP 软件把这套方法封装成了工程师手里的工具,使得"一塔评估一片场区"从经验判断变成了有据可依的标准化流程。

从定位上看,它处在"观测数据"与"工程决策"之间的桥梁位置。上游是气象站、测风塔的实测时间序列;下游是机位排布、年发电量(AEP)估算、投资决策。《欧洲风图集》的方法学正是把上游数据转译成下游可用结论的那套规则,因此被称为现代风电选址与产能评估软件的方法学源头并不为过。

方法 / 它做了什么

风图集方法的精髓,可以用"两步走"来概括:先"上行"(generalization,广义化),再"下行"(application,应用)。

第一步"上行"是去地形化。测风塔记录到的风,是当地粗糙度、地形起伏、周边障碍物共同作用后的"被污染"结果。方法把这些当地效应逐一剥离:用粗糙度长度 z0 及其在不同风向扇区上的变化来描述地表摩擦,用地形模型描述山丘与坡度带来的加速或减速,用障碍物模型描述建筑、林带在其尾流区造成的风速亏损。剥离之后,得到的就是一组"广义风气候"——可以理解为:如果把测站所在地变成一块理想的、平坦均匀粗糙度的开阔地,在若干标准高度和标准粗糙度下应有的风气候。这个广义风气候不再依赖某一个具体测点的地貌,因而具备了区域代表性。

支撑这一步的关键物理是地转拖曳定律(geostrophic drag law)。它把近地面摩擦层的风与高空自由大气中的地转风联系起来。其背后的核心假设是:在一片不太大的区域内,高空的地转风气候是基本一致的——也就是说,几公里范围内"天上吹的风"差不多,差异主要来自地面附近的摩擦、地形与障碍物。正因如此,只要把当地效应去掉、还原到高空,再换一个地点把当地效应加回去,就能完成跨点外推。

第二步"下行"是应用。拿到广义风气候后,针对目标机位重新引入它自己的粗糙度、地形与障碍物——相当于把"理想地块"重新"做旧"成目标点真实的地貌环境,从而预测出目标点轮毂高度处的风况。整个过程对每个风向扇区分别处理,因为不同来风方向看到的上风向地表和地形完全不同。

最终的输出形态,是分扇区的威布尔分布参数(尺度参数与形状参数)以及风功率密度。有了威布尔分布,结合机组功率曲线就能估算年发电量;分扇区信息则用于机位排布和尾流分析。这套"实测→广义风气候→目标点"的链路,正是后世风资源软件的标准工作流。

关键结论

  1. 风资源的跨点外推可以建立在清晰的物理框架上:通过去地形化得到广义风气候、再在目标点重新引入地形效应,而非简单的经验搬运。这是《欧洲风图集》方法学的核心贡献。
  2. 地转拖曳定律是连接高空地转风与近地面风的桥梁;"区域内高空地转风气候基本一致"是支撑跨点外推的关键假设。
  3. 当地效应需被分解为粗糙度(z0 及其扇区变化)、地形抬升、障碍物尾流三类分别建模,且需按风向扇区逐一处理。
  4. 方法的标准输出是分扇区威布尔分布参数与风功率密度,可直接对接产能评估与选址。
  5. 该方法对复杂地形(山地)和非中性大气采用了线性化假设,因而在陡峭、复杂地形条件下存在系统性误差;这一局限催生了后续以 CFD 为代表的非线性流场方法作为补强(CFD 补强属业界后续发展,非本报告内容)。

需要强调:上述第 1–5 点中所涉及的物理机制与方法论判断,均来自《欧洲风图集》本身;而具体到某地形下"误差有多大"的量化结论,原报告并未给出统一数字,工程实践中也高度依赖站点条件,不应臆造百分比。

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

风图集方法最大的现实约束,从来不是物理模型,而是"长期参考风气候"的可得性。方法要求有一份能代表长期气候、且时间序列足够长的参考风数据,作为广义化的输入与测风塔短期数据做长期订正(MCP)的锚点。传统做法是找附近的长期气象站,但气象站位置、迁站历史、观测高度往往不理想,离场区又远。

现代实践的转变在于:越来越多的项目改用再分析数据作为长期参考。运梦气象 API 的 ERA5(dataSourceId=era5)提供数十年的逐时长序列,覆盖全球、网格规整,正好补上"长期参考风气候"这块短板。取 100 米高度风分量 u100/v100,以及风速 ws、风向 wd,就能为某一片场区构造长期的、分扇区的风气候特征,供广义化与 MCP 长期订正使用。

具体到运梦平台的工作流,可以这样组织:用 ERA5 长序列刻画区域的长期风气候与年际波动,作为风图集方法里的"参考"输入;测风塔的短期实测则负责捕捉局地真实风况;两者通过长期订正(MCP)结合,再走"广义化—应用"链路外推到各机位。这样既保留了实测的局地准确性,又用 ERA5 弥补了实测序列偏短、难以代表长期气候的不足。

需要明确边界的是:再分析数据的网格分辨率有限,本身并不能替代复杂地形下的精细流场刻画。在山地等复杂场景,仍应叠加更精细的地形与流场方法(如 CFD)来订正 WAsP 线性模型的偏差。运梦 API 的角色是提供可靠、规整、长序列的长期参考输入,而非取代选址软件的流场计算。预报环节则使用德国气象局源(dataSourceId=ger),与用于资源评估的历史再分析在用途上严格区分开。

在运梦气象 API 上手

下面给出一个最小可用的取数示例:拉取目标场区附近一个网格点的 ERA5 长序列,字段取 100 米风分量 u100/v100 与风速风向 ws/wd,用于构造长期参考风气候、做广义化与 MCP 长期订正的锚点。把 stime/etime 拉到尽可能长的区间(示例仅取一段说明格式),即可得到逐时序列。

{
  "dataSourceId": "era5",
  "lat": 40.25,
  "lon": 116.18,
  "stime": "2014-01-01 00:00",
  "etime": "2023-12-31 23:00",
  "fields": ["u100", "v100", "ws", "wd"],
  "timezone": "8"
}

拿到序列后,先用 u100/v100 计算合成风速与风向(也可直接用返回的 ws/wd),按风向扇区统计威布尔分布参数,得到该网格点的长期风气候;再把它与测风塔短期实测做相关性建模(MCP),将测风塔序列订正到长期水平,最后送入风图集方法的"广义化—应用"链路完成机位外推。若需要短期运行预报,则另行调用德国气象局源(dataSourceId=ger)。

常见问题

问:ERA5 能直接替代测风塔吗? 答:不能。ERA5 的强项是"长期、规整、全球覆盖",适合做长期参考风气候和长期订正的锚点;但它的网格分辨率有限,无法刻画机位级的局地地形与障碍物效应。正确用法是 ERA5 提供长期背景、测风塔提供局地真实,两者通过 MCP 结合。

问:为什么要先"去地形化"再外推,而不直接搬运测风塔风速? 答:因为测风塔风速里混入了它所在地的粗糙度、地形和障碍物影响,这些影响在目标机位完全不同。去地形化先把这些当地效应剥离、还原到广义风气候,再在目标点把目标点自己的地形效应加回去,才能保证外推的物理一致性。

问:WAsP 方法在山地会失准吗? 答:会。该方法对复杂地形和非中性大气采用线性化假设,在陡峭、复杂地形下存在系统性偏差。这正是业界后续引入 CFD 类非线性流场方法做补强的原因。运梦 API 不替代流场计算,只负责提供可靠的长期参考输入。

问:应该取哪些字段? 答:风资源评估推荐 u100/v100(100 米风分量)配合 ws/wd(风速、风向)。100 米接近现代机组轮毂高度,u100/v100 便于按扇区合成与统计威布尔分布。如需 10 米风可另取 uas/vas。

问:资源评估和功率预报用同一个数据源吗? 答:不是。历史资源评估用 ERA5 再分析(dataSourceId=era5)长序列;短期运行功率预报用德国气象局源(dataSourceId=ger)。两者用途不同,不应混用。

引用与原文

Troen, I., & Petersen, E. L. (1989). European Wind Atlas. Risø National Laboratory, Roskilde, Denmark. ISBN 87-550-1482-8.

该报告无 DOI;引用时以书目信息(ISBN 87-550-1482-8)为准。