第一章：风资源评估基础
各位好，我是老张。干风资源评估这行十几年了，今天咱们聊聊最基础、也最容易被忽视的东西——测风塔数据怎么处理、风切变和湍流强度怎么算、威布尔分布怎么拟合，还有年发电量怎么初步估算。
说实话，我见过太多人一上来就搞复杂的CFD模拟，结果基础数据都没处理好。嗯，这就像盖房子不打地基，早晚要出问题。
1.1 测风塔数据采集与处理
测风塔，说白了就是我们在风电场里竖的一根杆子，上面挂满风速仪、风向标。我习惯在10m、50m、70m、90m、100m这几个高度都装传感器。为什么？因为风切变需要多个高度的数据才能算准。
核心要点：测风塔数据采集的三大原则——完整性、代表性、准确性。
数据采集频率一般是10分钟一个点。一天144个点，一年52560个点。数据量不小，但真正干净的数据往往不到80%。
我在内蒙古一个项目上遇到过，测风塔数据突然少了两个月。后来一查，是传感器被冻住了。所以啊，数据预处理这步绝对不能省。
数据预处理步骤
- 完整性检查——看看有没有缺失时间段。缺失超过10%的月份，我建议直接剔除。
- 合理性检查——风速不能为负，不能超过60m/s。风向在0-360度之间。温度不能低于-50℃。
- 一致性检查——同一高度两个风速仪的数据差异不能超过5%。我遇到过差异20%的情况，结果是其中一个传感器坏了。
- 时间序列修正——把异常值标记出来，用插值法补全。
我的小技巧：处理数据时，我习惯先画个时间序列图。一眼就能看出哪些数据有问题。比看表格快多了。
1.2 风切变与湍流强度计算
风切变，说白了就是风速随高度变化的规律。你想想看，地面有摩擦，树啊房子啊都会挡风，所以越靠近地面风速越小。到了高空，风就顺畅了。
计算风切变，我常用两种方法：
方法一：幂律公式
V(z) = V(z_ref) × (z / z_ref)^α其中α就是风切变指数。我一般用最小二乘法拟合出α值。
方法二：对数公式
V(z) = (u*/k) × ln(z/z0)u*是摩擦速度，k是卡门常数（约0.4），z0是地表粗糙度。
| 地表类型 | 粗糙度z0 (m) | 典型α值 |
|---|---|---|
| 开阔水面 | 0.0002 | 0.10-0.12 |
| 草地 | 0.03 | 0.14-0.16 |
| 农田 | 0.1 | 0.18-0.22 |
| 森林 | 0.5-1.0 | 0.25-0.35 |
湍流强度呢？它反映的是风速的波动程度。计算公式很简单：
TI = σ / V_meanσ是风速标准差，V_mean是平均风速。TI值小于0.1算低湍流，0.1-0.2算中等，大于0.2算高湍流。
注意：我曾经在山区项目上吃过亏。地形复杂的地方，湍流强度能到0.3以上。这种地方装风机，叶片疲劳载荷会很大。一定要慎重。
1.3 威布尔分布拟合
为什么要拟合威布尔分布？因为我们需要知道一年里各个风速段出现的概率。说白了，就是给风速做个统计画像。
威布尔分布有两个参数：形状参数k和尺度参数c。k值一般在1.5-3.5之间。k越大，风速分布越集中。c值接近年平均风速。
拟合方法我常用两种：
方法一：最小二乘法
ln(-ln(1-F(v))) = k × ln(v) - k × ln(c)把数据点画在图上，拟合直线。斜率就是k，截距可以算出c。
方法二：矩估计法
V_mean = c × Γ(1 + 1/k)
σ² = c² × [Γ(1 + 2/k) - Γ²(1 + 1/k)]用迭代法解出k和c。
我的习惯：两种方法都算一遍，对比结果。如果差异超过5%，说明数据可能有问题。我会回头检查原始数据。
1.4 年发电量初步估算
年发电量估算，说白了就是算风机一年能发多少电。我习惯用三步法：
- 算风速频率分布——用威布尔分布算每个风速段出现的小时数。
- 查功率曲线——每个风速对应多少功率，风机厂家会给。
- 加权求和——风速频率 × 对应功率，累加起来就是年发电量。
公式长这样：
AEP = 8760 × ∫ f(v) × P(v) dvf(v)是风速概率密度函数，P(v)是功率曲线，8760是一年的小时数。
重要提醒：初步估算时，别忘了考虑折减系数。空气密度折减、尾流折减、湍流折减、停机维护折减……林林总总加起来，实际发电量大概只有理论值的75%-85%。
我在新疆一个项目上，初步估算年发电量是2.5亿度。结果实际运行下来只有2.1亿度。后来一查，是忽略了冬季低温导致的空气密度下降。嗯，这个教训我一直记着。
知识体系总览
下面这张图，是我自己总结的风资源评估基础框架。你看一眼，心里就有数了。
好了，这一章的内容就这些。数据采集处理是基础，风切变和湍流是核心参数，威布尔分布是统计工具，年发电量估算是最终目标。每一步都环环相扣，哪一步出了问题，后面的结果都不靠谱。