第2章大气边界层基础

大气边界层基础：大气稳定度、风速廓线、湍流强度、地表粗糙度

各位同学好，我是老张。在风电场干了十几年气象模拟，今天咱们聊聊大气边界层。说白了，风电机组就活在这个层里——从地面往上大概几百米到一公里的范围。你想想看，叶片转不转、能发多少电，全看这层大气的"脸色"。

我个人习惯，讲边界层之前先问一个问题：为什么同一台风机，上午和下午的发电量能差30%？答案就藏在大气稳定度和风速廓线里。咱们一个一个来拆。

大气稳定度，说白了就是大气"想不想"上下翻腾。我遇到过不少刚入行的同事，觉得这概念太虚。其实它直接影响尾流恢复的速度——稳定层结下，尾流能飘出去十几公里都不散。

核心分类：

判断稳定度，我常用一个简单指标——理查森数 Ri。公式不复杂：

Ri = (g/θ) * (dθ/dz) / (du/dz)²

其中 g 是重力加速度，θ 是位温，u 是风速。Ri < 0 是不稳定，Ri > 0.25 是稳定，中间算中性。嗯，这里要注意：现场实测的 dθ/dz 经常跳变，我建议取10分钟滑动平均。

我的经验：在内蒙古某风场，白天不稳定条件下，尾流在5倍叶轮直径后就恢复了80%。但到了夜间稳定层结，同样距离只恢复了40%。这就是为什么夜间发电量"虚高"——前排吃风，后排喝汤。

风速不是均匀的。离地面越近，摩擦越大，风速越小。这个变化规律就是风速廓线。搞风资源评估，这是基本功。

最经典的模型是对数律：

u(z) = (u*/κ) * ln(z/z0)

u* 是摩擦速度，κ 是卡门常数（约0.4），z0 是地表粗糙度。说白了，风速随高度按对数增长。但我在海上风场发现，实际数据经常偏离对数律——这时候要用幂律：

u(z) = u_ref * (z/z_ref)^α

α 是风切变指数，一般在0.1到0.4之间。我建议你手头至少准备三组α值：白天不稳定用0.1，中性用0.15，夜间稳定用0.3。这样估算八九不离十。

避坑指南：我曾经在云南山地风场吃过亏——直接用平坦地形的α=0.14算轮毂高度风速，结果偏差了12%。山地、森林、城市，α值能差一倍。老老实实架测风塔，别偷懒。

湍流强度 TI，就是风速的"抖动"程度。公式很简单：

TI = σ_u / U

σ_u 是风速标准差，U 是平均风速。TI 越大，风机载荷越大，尾流混合越快。

我见过最极端的案例：某沿海风场，台风过境时 TI 飙到0.4，风机直接切出停机。正常情况 TI 在0.1到0.2之间。地表粗糙度越大，TI 越高——这是物理直觉，对吧？

实用技巧：做尾流模拟时，我习惯把 TI 分成两部分——背景湍流（由地表粗糙度决定）和风机 induced 湍流（由推力系数决定）。两者叠加才是真实的尾流湍流场。别直接用测风塔的 TI 去算，会低估。

地表粗糙度 z0，是个等效参数。它不直接等于草的高度或房子的高度，而是反映地面对气流的拖曳作用。我刚开始做这行时，总觉得 z0 是个"玄学"参数——其实不然。

确定 z0 的方法有三种：

我个人最推荐第二种。在项目现场，至少架3层测风仪（10m、50m、80m），用最小二乘法拟合对数律，z0 和 u* 一起解出来。我曾经在甘肃戈壁滩这么干过，算出来的 z0=0.005m，跟查表值完全吻合——那种成就感，嗯，你做了就知道。

关键影响：z0 每增加一个数量级，轮毂高度风速大约降低5-8%。对于100MW的风场，这意味着一两千万的年发电量差异。所以别小看这个"小"参数。

下面这张图，是我自己梳理的大气边界层核心逻辑。四个参数不是孤立的——稳定度影响廓线形状，粗糙度决定湍流强度，湍流强度又反过来改变稳定度。搞明白这个循环，尾流模拟就入门了。

好了，这一章的内容就到这儿。大气边界层是风电场气象模拟的根基——稳定度、廓线、湍流、粗糙度，这四个参数你吃透了，后面讲尾流模型、机舱传递函数、发电量折减，都会顺很多。下次咱们聊点更具体的：怎么用实测数据反推这些参数。