1、数值预报基础：大气运动方程、数值模式分类、全球与区域模式区别
各位同事，咱们今天聊聊数值预报的底子。说实话，我刚开始接触风电调度那会儿，觉得数值预报就是个黑盒子——数据进去，结果出来，准不准全看运气。后来干久了才发现，不懂底层原理，你连预报为啥不准都说不清楚。
这一节，咱们把地基打牢。我会尽量用大白话讲，但该有的专业深度一点不少。
1.1 大气运动方程——预报的物理根基
数值预报到底在算什么？说白了，就是解一组方程。这组方程描述的是大气怎么运动、怎么变化。
核心方程有七个，我列个表，大家有个印象就行：
| 方程名称 | 物理含义 | 在风电调度中的意义 |
|---|---|---|
| 运动方程（纳维-斯托克斯） | 描述风速、风向的变化 | 直接决定风机出力预测 |
| 连续方程 | 质量守恒，空气不凭空消失 | 影响风场辐合辐散判断 |
| 热力学方程 | 温度变化与能量收支 | 影响大气稳定度，进而影响风速 |
| 水汽方程 | 水汽的输送与相变 | 影响雾、降水，间接影响风功率 |
| 状态方程 | 气压、温度、密度的关系 | 基础参数转换 |
嗯，这里要注意——这些方程是非线性的，没有解析解。什么意思？就是你不能像解一元二次方程那样直接写出答案。只能靠计算机一点一点算，这就是数值预报名字的由来。
我个人习惯：每次拿到新的预报产品，先看它用的动力核心是什么。不同的动力框架，对复杂地形的处理能力差别很大。比如我在西北项目上遇到过，同一个初始场，用谱方法和用有限差分法算出来的风速，在山区能差2-3m/s。
1.2 数值模式分类——从空间维度看
数值模式怎么分类？我一般从两个维度看：空间维度和时间维度。
从空间维度，主流分三类：
- 静力平衡模式：假设垂直方向气压梯度力与重力平衡。适合大尺度，水平分辨率>10km。计算快，但处理对流天气不行。
- 非静力平衡模式：不假设垂直平衡，能模拟对流、湍流。现在的主流，比如WRF、ECMWF的IFS都是这类。
- 大涡模拟（LES）：直接解析大尺度湍流涡旋。分辨率极高（<100m），但计算量巨大，目前只用于科研或风电场微观选址。
从时间维度，分确定性预报和集合预报。这个后面章节会细讲，今天先记住：确定性预报给一个答案，集合预报给一堆答案，然后告诉你概率。
避坑指南：我曾经在海上风电项目上，直接用静力平衡模式做72小时预报，结果台风路径偏了200公里。后来换成非静力模式，误差缩小到50公里以内。所以，做风电调度，至少要用非静力平衡模式。
1.3 全球模式与区域模式的区别
这个问题，我每次培训都会被问到。咱们直接上对比表：
| 对比项 | 全球模式 | 区域模式 |
|---|---|---|
| 覆盖范围 | 整个地球 | 特定区域（如中国、华北） |
| 典型分辨率 | 9-25km（ECMWF、GFS） | 1-9km（WRF、COSMO） |
| 计算资源 | 超级计算机，全球共享 | 本地或云端，相对灵活 |
| 边界条件 | 自洽，不需要外部输入 | 需要全球模式提供边界 |
| 预报时效 | 通常7-16天 | 通常48-72小时 |
| 对风电调度的价值 | 中期趋势判断（3天以上） | 短期功率预测（0-72小时） |
你想想看，全球模式就像一张世界地图，能看到大尺度天气系统——比如西风带波动、副热带高压。但你想看某个风电场门口那个小山包对风的影响？全球模式根本看不见。
区域模式就不一样了。它把全球模式的粗结果拿过来，在局部区域用更细的网格重新算一遍。分辨率高了，地形、海岸线、城市热岛这些细节就能体现出来。
我的经验：在内蒙古做风电调度时，我习惯用ECMWF的全球模式看大趋势，然后用WRF区域模式做72小时内的精细化预报。两者结合，准确率能提升15-20%。但要注意——区域模式的边界每6小时更新一次，如果全球模式本身跑偏了，区域模式也会跟着偏。这叫"垃圾进，垃圾出"。
1.4 知识体系框架
下面这张图，是我自己总结的数值预报知识体系。每次带新人，我都先让他们看这张图：
这张图的核心逻辑是：物理方程是根，数值方法是干，模式分类是枝，风电应用是果。根不牢，什么都白搭。
重要提醒：不要以为有了高分辨率区域模式，全球模式就没用了。我见过不少调度员，只看区域模式的结果，结果遇到大尺度天气系统调整（比如副高突然北跳），区域模式完全没报出来，因为它的边界条件来自全球模式，而全球模式在边界更新之前就已经变了。
所以我的建议是：全球模式看趋势，区域模式看细节，两者缺一不可。
好了，这一节的内容就这些。记住三个关键词：方程、分类、区别。下一节咱们会深入讲数据同化——就是怎么把观测数据塞进模式里，让预报更准。这个环节，我在实际项目中踩过不少坑，到时候跟大家好好聊聊。