🌬️ 风能风资源计算与评价

本笔记以风电场工程的实际推进流程为主线，逐步梳理并记录各阶段的计算逻辑、操作要点与规范要求。

---

一、第一步：导入风机、测风塔及地形图数据（空间校验）

在开展任何风速和能量计算之前，首要任务是在软件（如 WindPRO、WAsP、WindSim 等）或地理信息系统（如 Global Mapper、ArcGIS）中建立正确的三维空间场景，校验输入数据的坐标及空间代表性，避免“差之毫厘，谬以千里”。

1.1 输入数据清单与准备

在导入阶段，需要准备好以下空间地理数据：

1. 风电机组坐标表 ：包含初拟风电机组（WTG）的编号、X 坐标、Y 坐标、轮毂高度以及初拟机型，格式通常为 Excel、CSV 或 KML/KMZ 格式。
2. 测风塔基本信息表 ：记录测风塔的编号、经度、纬度、海拔高度、传感器安装层高以及地磁偏角，可参考规范 [NB/T 31147-2018 附录 A.2.1](../../references/技术规范/风电/风电场工程风能资源测量与评估技术规范（NB-T 31147-2018）.md#a2-基本信息)。
3. 地形高程数据（DEM） ：描述场区及周边地形起伏的数据。常用格式为 GeoTIFF（如 SRTM 分辨率或 ASTER GDEM 数据）或者由矢量地形图转出的格网/等高线数据（如 DXF、MAP 格式）。
4. 地表粗糙度/下垫面文件 ：反映森林、农田、草地等不同地表覆盖类型对风速的阻碍作用。

二、第二步：3D 模式可视化校验（地形匹配与遮挡核查）

在完成基础空间数据的二维导入后，第二步必须使用 3D 模式 （如在 WindPRO 中使用 3D-Hills、Google Earth 3D 视图，或者在 WindSim/WAsP 中的三维地形视图）进行全方位的立体视觉校验。这一步骤主要是检验测风塔位置的三维合理性，排除二维平面分析中无法识别的地形及遮挡陷阱。

2.1 地形贴合与海拔核准（排除高程漂移）

• 悬空与下埋检验 ：在 3D 场景中旋转视角，重点检查测风塔模型基座是否与 DEM 地面精准重合。若测风塔在 3D 场景中表现为“悬挂于半空”或“埋入山体内部”，通常说明导入数据时的高程基准（如不同投影高程系统、DEM 分辨率偏差）不匹配，或者坐标点定位有误。
• 高程一致性对齐 ：提取 3D 模型在地形格网上的基础海拔高度，与测风塔基本信息表中的 GPS 实测海拔进行对齐。对于复杂地形，若海拔高度差值较大，应以实测高程或大比例尺高精地形图为准，重新对齐 3D 场景中的基础海拔。

2.2 微地形代表性与风障分析

• 局部微地形判定 ：测风塔应当布设在能代表场区宏观风力特征的位置（如平缓的脊顶）。在 3D 模式下，可以直观判断测风塔是否陷于极其局部的“微地形挡风面”（如陡峭的断崖正下方、孤立山峰背后的涡流区）。
• 山体坡度检验 ：利用 3D 分析检查测风塔前侧的迎风坡度。若前侧迎风坡太陡（如坡度超过 ），会导致风发生急剧的气流分离，实测得出的风速和湍流会产生局部变异（畸变），无法代表该区域风电机组的正常风况。

2.3 三维障碍物及邻近遮挡核查

• 高大林木及构筑物遮挡 ：在 3D 场景中加入三维障碍物（如林带高度、附近房屋高度等）。
  • 结合规范，在 3D 模式下直观校验测风塔各高度的传感器（尤其是最高层及初拟轮毂高度处）是否被周围树冠或建筑物遮挡。
  • 若测风塔位于高大密集的树林内部，应结合 3D 空间结构检查顶部传感器是否高出森林冠层至少 10 m ，以获取不受树木湍流干扰的自由流风速。
• 与邻近风机的 wake 遮挡校验 ：如果在已建或拟建风机附近测风，可在 3D 模式下模拟各盛行风向的流线走向。核实测风塔是否直接处于邻近风电机组的主风轮扫掠 wake 区域，避免实测风速偏低和湍流偏高。

2.4 测风支臂（Boom）朝向与塔影效应规避核对

• 支臂朝向 3D 立体校对 ：根据场址的盛行风向，在 3D 模式下检查各个测风通道的水平均风支臂（Boom）的朝向设计。
  • 横梁与主风向宜成 （即支臂朝向主风向）。
  • 检查传感器支臂的夹角（如双风速通道设计），从三维视角确保最高层测风设备不受塔架自身的下风向湍流区（即避开 塔影效应 区域）干扰。

---

三、第三步：敏感性因素排查（三区三线与约束红线核查）

在完成空间布置和 3D 可视化地形校验后，第三步是进行敏感性因素排查。这主要是将拟布置的风机、测风塔、进场道路、集电线路及升压站等坐标数据，与国家及地方提供的 “三区三线” （及其他各类红线）空间地理矢量数据库进行叠加分析，排查是否存在法规或政策性的一票否决（颠覆性）因素。

3.1 “三区三线”空间红线套核

• 生态保护红线 ：
  • 严禁任何风电机组、升压站、道路及电缆桥架侵入 生态保护红线 边界。
  • 在 GIS 中将边界矢量数据与项目布局进行高精度重叠（Overlap）分析，对距离红线边界较近（如 以内）的风机，建议进行退让以防施工过程中意外越界。
• 永久基本农田 ：
  • 风电项目的永久征地（如风机塔筒基础、升压站、永久道路）原则上严禁占用 永久基本农田 。
  • 临时占地（如吊装广场、集电线路基坑）若必须占用，需在后期执行极严格的复垦及恢复手续，并需经过专题报批。排查阶段应尽可能避开，以规避合规风险。
• 城镇开发边界 ：
  • 核算风机与 城镇开发边界 及居民规划区的距离。若风机距离开发边界太近，不仅会面临噪音、光影闪烁（Shadow Flicker）等环保纠纷，也极易因后续城镇规划调整而面临被动拆除。

3.2 其他敏感限制性因素排查

除了“三区三线”，还需同步核查以下敏感性边界：

• 林地属性排查 ：区分防护林、国家级公益林、地方林地等。国家级公益林及防护林严禁建设风电。
• 文物与遗迹保护区 ：各级文物保护单位、风景名胜区、自然遗产地以及饮用水水源保护区均有相应的法律限制，严禁在核心区及缓冲区布置设备。
• 矿产压覆与地质灾害 ：核实风机基础是否压覆重要矿产资源（需征得矿权人同意或避开）；核实机位是否处于山体滑坡、泥石流等 地质灾害隐患点 。
• 军事与雷达影响 ：核实风电场是否处于军事演习区、雷达航路或微波通道的遮挡与干扰范围内，必要时需取得军委和气象部门的正式复函。

3.3 安全退让距离校核

将风机点与地图上的既有民用/公共基础设施进行距离量算，确保满足以下安全退让基准：

• 居民区/村庄 ：一般建议退让 1.5 倍至 2 倍 叶轮直径（或不小于 ），防止噪音和冰投射（Ice Throw）安全风险。
• 高压输电线路/铁路/高速公路 ：风机机位至高压线、铁路及高速公路边缘的距离，应不小于 1.5 倍 倒塔安全距离（即 ）。

---

四、第四步：测风原始数据读取与软件导入（NRG SymphoniePRO / Windographer）

在完成前期所有的空间位置及敏感性校验后，项目将正式进入测风数据分析阶段。针对测风塔下载的原始测风数据，首要工作是使用专用的测风管理与分析软件进行数据读取、格式解码与初步的时间序列展示。

4.1 NRG SymphoniePRO Desktop 原始数据处理

如果测风塔使用的是 NRG 记录仪（如 SymphoniePRO 系列），下载的原始文件通常为加密的二进制格式（如 .rld 文件）。

• 数据导入与解码 ：打开 NRG SymphoniePRO Desktop 软件，将原始的 .rld 二进制文件包导入数据库。软件会自动将其解码还原为明文的原始时间序列数据。
• 通道与标定核查（Channel Configuration） ：在软件中核对各通道（Channel）的传感器配置和标定参数（坡度 Slope 与截距 Offset）。确保风速仪和风向标的标定参数与仪器随附的出厂标定证书完全一致，避免由于公式参数配置错误导致读入的风速偏大或偏小。
• 文本数据导出 ：校验通道无误后，通过软件批量导出为标准的 10 分钟平均、标准差、极大极小风速文本文件（如 .csv 、 .txt 格式），供通用评估软件读取。

4.2 Windographer 数据集成与时间序列初判

对于通用的测风数据分析、清洗与定级，业内标准软件是 Windographer 。

• 多格式直接支持 ： Windographer 能够直接读取 NRG 的二进制文件（如 .rld ）或者解码后的文本 CSV 文件。导入时，软件会引导用户定义各个通道的高度、类型（风速、风向、温度、气压等）。
• 时间序列图可视化初判 ：读入数据后，第一步应当生成 时间序列图（Time Series Plot） 。通过横向 and 纵向滚动，初步审视数据流的整体连续性：
  • 检查是否存在大面积的“数据中断/数据缺测”段。
  • 检查传感器是否发生“结冰死线”（风速长时间为零或极小，风向标准差为零，且气温在 以下）。
  • 检查是否存在“超出物理范围”的野点数据（如突发的异常大风速，或气温突变）。

4.3 每日数据的批量导入与合并（解决单日分散文件问题）

由于测风数据下载时通常以天为单位生成独立文件，我们需要将其一次性批量合并打开：

• 在 Windographer 中批量拼接 ：
  1. 点击菜单 File -> Open ，或者将所有单日文件所在文件夹拖入软件。
  2. 在文件选择框中，使用键盘快捷键 Ctrl+A 全选所有天数的原始数据文件，或者直接在弹出的向导中选择整个数据文件夹。
  3. 点击打开，软件会自动识别各文件的起止时间并将其按时间线前后拼接，形成一条完整且连续的时间序列数据集。
  4. 如果后续有新增的一天或多天数据，可以在已打开的数据集上点击 File -> Append Data... （追加数据）进行增量合并。
• 在 NRG SymphoniePRO Desktop 中批量处理与导出 ：
  1. SymphoniePRO 软件基于数据库管理模式。切换到 Import （导入）界面。
  2. 点击 Add Folder （添加文件夹）按钮，直接选择存放所有单日 .rld 文件包的本地文件夹。
  3. 软件会批量自动解密解码并导入对应的测风塔数据库。
  4. 导入完成后，在主界面中选择对应的测风站点，进入 Export （导出）界面，设置您要导出的整个时段（例如整年或从第一天到最后一天），软件即可一键导出包含全部天数数据的完整文本文件。

4.4 避坑：Windographer 导入导出的精度损失问题（默认保留 3 位小数）

在实务操作中需特别注意： Windographer 在默认配置下，无论是软件内的数据显示，还是执行数据导出时，可能会将数据格式化并截断为 小数点后 3 位 （例如风速标准差、平均值截断为 ）。这对于一般的风速统计可能足够，但对于以下高精度计算会带来累积误差（精度损失）：

• 风切变指数拟合 ：在利用最小二乘法进行对数线性回归拟合时，小数位被截断会微幅改变回归直线的斜率（即风切变指数 ）。
• 威布尔分布拟合（Weibull parameters） ：拟合形状参数 与尺度参数 时，输入数据的截断误差会导致拟合曲线发生微调，进而影响估算发电量。
• 湍流强度（TI）计算 ：标准差与平均风速的轻微舍入误差，在除法计算中可能导致湍流强度数据产生不自然的阶梯状跳变。

🛠️ 解决方法 ：

1. 全局显示与存储精度调整 ：进入软件的 Tools -> Options （或 Preferences ），找到 Units & Decimals （单位与保留位数）选项，手动将风速、标准差、气压等核心通道的十进制保留位数（Decimal Places）调整为 4 位或 5 位 以上。
2. 数据导出精度控制 ：当在 Windographer 中清洗完数据并向外导出（如导出为文本或准备用于 WindPRO 的 .txt 文件）时，在导出设置弹窗中，务必将保留小数位数（Decimal Places）设置为“保留全部有效位数（All significant digits）”或手动指定一个较高的保留位数（如 5 位或 6 位 ），防止软件默认截断导致数据失真。

---

五、第四步（进阶方案）：基于 Python 自主编程处理测风数据

由于商业软件存在精度截断、闭源黑盒、难以进行大规模定制化批量处理等局限，采用 Python 自主编写代码来进行测风数据的读取、清洗、插补与参数计算，是极具实用价值的进阶选择。

5.1 自主编程处理的优势

• 无精度损失 ：使用 Python (numpy / pandas) 处理时，数据默认以双精度浮点数（ float64 ）存储和计算，规避了商业软件小数点后 3 位的强制截断。
• 流程高度自动化 ：对于每日生成的测风文件，用几行代码即可完成全自动的批量读取、时间排序、缺失值对齐与数据拼接。
• 算法完全透明 ：无论是风切变的对数线性拟合，还是威布尔参数的参数估计，公式和权重逻辑完全由代码掌控，可审计、可复现。
• 数据清洗高度定制 ：可以根据测风塔传感器实际的表现，自定义极复杂的关联条件逻辑（例如：结合气温与风速波动特征，开发自动判定凝冻结冰的 QC 算法）。

5.2 核心 Python 数据处理架构设计

自主编写脚本时，推荐的数据处理管道（Pipeline）架构设计如下：

import os
import glob
import pandas as pd
import numpy as np
from scipy.optimize import curve_fit

def batch_load_and_combine(folder_path):
    """
    1. 批量读取每日测风文件并按时间线拼接
    """
    # 匹配文件夹下所有的原始 CSV / TXT 文件
    file_list = sorted(glob.glob(os.path.join(folder_path, "*.csv")))
  
    df_list = []
    for file in file_list:
        # 根据 NRG 导出的格式设定 header 和 index_col
        temp_df = pd.read_csv(file, parse_dates=['Timestamp'], index_col='Timestamp')
        df_list.append(temp_df)
  
    # 沿时间轴纵向拼接数据，并根据时间排序去重
    combined_df = pd.concat(df_list).sort_index()
    return combined_df

def quality_control(df):
    """
    2. 数据合理性初判与数据清洗 (QC)
    按照规范标准，对异常点赋 NaN，不破坏时间索引的完整性
    """
    cleaned_df = df.copy()
  
    # 范围检验示例 (例如 10m 风速范围 0 ~ 40 m/s)
    cleaned_df.loc[(cleaned_df['WS_10m'] < 0) | (cleaned_df['WS_10m'] > 40), 'WS_10m'] = np.nan
  
    # 趋势检验示例 (例如风速连续 6 小时无任何变化判定为传感器故障)
    # 10 min 记录下 6 小时共 36 个 point
    is_stuck = cleaned_df['WS_10m'].rolling(window=36).std() == 0
    cleaned_df.loc[is_stuck, 'WS_10m'] = np.nan
  
    # 相关性检验示例 (例如高层风速小于低层风速超过合理阈值判定为反常)
    anomalous_shear = (cleaned_df['WS_hub'] - cleaned_df['WS_10m']) < -2.0
    cleaned_df.loc[anomalous_shear, ['WS_hub', 'WS_10m']] = np.nan
  
    return cleaned_df

5.3 实用工具：时间序列空缺行自动补全（fill_timeseries_gaps.py）

在实际测风过程中，由于记录仪断电、存储卡故障或传输中断，原始数据导出的文本中往往会产生“时间戳行缺失”的断档。这会导致时间序列不连续，无法直接进行步长滑动计算或插补。

为此，开发了专用的数据补全工具 fill_timeseries_gaps.py 。

1. 工作原理与核心逻辑

该脚本专门针对 Tab 键分隔 （Tab-delimited）的测风文本导出文件（如 SymphoniePRO 的导出文本）进行流式处理：

• 保持头部元数据不变 ：自动检索定位包含 Timestamp （可通过参数修改）的列标题行，并在其上方的所有非结构化元数据（Metadata/Header）原样保留，不破坏文件属性。
• 自动推断测风步长 ：若未指定频率，脚本会统计数据中出现频次最高的相邻时间戳差值，作为期望的步长（如自动识别为 10 min 或 1 min ）。
• 空位自动插行占位 ：按步长递增时间戳。当发现时间戳断档时，自动在缺损处插入对应的时间戳行，并为该行的所有数据列填充空值（保留 Tab 分隔符），确保输出序列的绝对连续。
• 兼容编码与换行符 ：自动识别文件的 BOM 头 （Byte Order Mark）与系统换行符（CRLF/LF），保证读写前后编码与排版无损。

2. 使用方法与 CLI 命令行参数

可以在命令行中通过以下参数调用该脚本：

python private/个人技术/fill_timeseries_gaps.py <输入测风文件路径> [选项]

常用选项说明 ：

• -o, --output ：指定输出文件路径。默认自动在输入文件名后追加 _continuous 后缀（如 site01_continuous.txt ）。
• --timestamp-column ：指定时间戳列的标题名称（默认值为 Timestamp ），用于定位数据起始行。
• --freq-minutes ：强制指定测风频率分钟数（如 10 或 1 ），不指定时则由系统自动推断。
• --report-only ：只打印出缺失时间戳的明细报告，不实际写入输出文件（可用于首期缺测诊断）。
• --force ：如果输出文件已存在，强制覆盖写入。

3. 补全后效果

处理后的数据集将形成完美连续的网格。缺失的时段在文本中表现为：

2026-06-18 00:00:00\t12.4\t312\t...\n
2026-06-18 00:10:00\t\t\t...\n  <-- 自动补全的空数据行
2026-06-18 00:20:00\t12.7\t315\t...\n

这使得后续在进行 Python 中 df.interpolate() 线性插补、 df.rolling() 滑动窗口计算或者 WAsP 转换时，拥有了齐整的时间网格基准。

---

六、第五步：测风数据标准化与字段清洗（过滤 Min、Max、Gust 等）

当多天原始测风文本成功合并并对齐时间线后，第五步是进行数据标准化处理。原始记录仪文件（如从 SymphoniePRO 导出的数据列）中通常包含大量冗余的统计列，例如各通道的 Min （最小风速）、 Max （最大风速）以及 Gust （3 秒极大阵风/极大风速）等。

在进行风特征参数统计、风切变拟合及 WAsP 计算等主流风资源评估流程时，我们通常只需要 平均值（Avg） 和 标准差（SD） 。过滤冗余列不仅可以使数据文件体积缩小、提高计算效率，还能统一数据接口，便于后续批量自动化脚本的开发。

6.1 过滤与标准化流程

1. 确定保留核心列 ：
   • 风速通道 ：保留平均风速（如 WS_100m） and 标准差（SD，如 WS_SD_100m，用于计算湍流强度）。
   • 风向通道 ：保留平均风向（如 WD_100m）。
   • 气象通道 ：保留平均气温（Temp_Avg）、平均大气压（Pres_Avg） and 平均相对湿度（RH_Avg）。
2. 剔除冗余字段 ：
   • 全局扫描列名，将所有带有 Min、Max、Gust 等辅助统计与诊断字段的列直接丢弃。

6.2 实务操作方法

根据项目条件与批量性，可以选择手动 Excel 过滤或脚本自动化处理：

1. 方法一：Excel 手动处理（当前采用）

对于单个或少量的合并文本，直接在 Excel 中操作非常直观：

1. 数据导入 ：以 Tab 键分隔 格式将合并后的测风文本（或 .csv）导入 Excel。
2. 列筛选与删除 ：整列选中包含 Min （最小值）、 Max （最大值）、 Gust （极大值）以及传感器工作电压等非统计核心列，右键选择“删除”。
3. 字段批量替换 ：使用快捷键 Ctrl + H 打开“查找和替换”对话框：
   • 查找内容输入： _Avg
   • 替换为：留空
   • 点击“全部替换”，一键将所有 WS_Avg_xxx 、 WD_Avg_xxx 重命名为 WS_xxx 和 WD_xxx 。
4. 保存导出 ：将处理后的表格另存为“文本文件（制表符分隔） *.txt ”或标准 *.csv 格式，完成标准化。

2. 方法二：Python 代码自动化（备选方案）

当测风塔数量较多或需要高频次迭代计算时，可以在您的 Python 处理流（Pipeline）中集成以下处理逻辑：

def standardize_measurement_data(df):
    """
    数据列标准化清洗：
    1. 过滤掉 Min, Max, Gust, Volt 等不需要用于评估计算的辅助统计列
    2. 规范化剩余的核心列名 (去掉 _Avg 标签，保留 _SD 标签)
    """
    # 1. 识别并过滤冗余列
    redundant_keywords = ['Min', 'Max', 'Gust', 'Volt', 'Batt']
    columns_to_drop = [
        col for col in df.columns 
        if any(keyword in col for keyword in redundant_keywords)
    ]
    df_filtered = df.drop(columns=columns_to_drop)
  
    # 2. 批量重命名核心列
    rename_rules = {}
    for col in df_filtered.columns:
        if '_Avg' in col:
            # 例如: WS_Avg_100m -> WS_100m, WD_Avg_100m -> WD_100m
            rename_rules[col] = col.replace('_Avg', '')
        elif '_SD_' in col and '_Avg' in col:
            # 兼容有些记录仪带有 WS_SD_Avg 的复杂命名
            rename_rules[col] = col.replace('_Avg', '')
          
    df_standardized = df_filtered.rename(columns=rename_rules)
    return df_standardized

---

七、第六步至第十一步：基于公司自主研发的在线电算平台进行云端计算处理

在完成前五步的测风数据补全与标准化清洗后，后续的复杂电算化工作（如合理性校验、插补修复、MCP 长系列订正、风切变及威布尔参数拟合、50 年一遇极端风速计算等）无需自行编写本地脚本，而是直接通过公司自主研发的 风资源在线电算化平台 进行直接电算处理。

7.1 平台接入与项目准备

• 电算平台网址 ： 风资源评估系统 （由公司开发，可直接在线进行电算处理）
• 登录账号 ： 18423509125 （潘刚的账户，密码需自行填写并保管）
• 操作流程 ：

1. 访问上述网址登录平台。
2. 登录后进入系统配置（System Configuration）或新建/导入对应的项目工作区。
3. 上传前述在 Excel 或 Python 中整理完成的、去除了 Min/Max/Gust 冗余列的标准化连续测风文本文件（ .txt 或 .csv 格式）。
4. 平台在云端会自动读取、建档并一键启动以下核心计算流程。

7.2 平台云端电算化计算逻辑与数学算法（技术底层）

虽然计算在云端平台自动运行，但我们需要掌握其底层执行的规范逻辑与数学公式：

1. 数据合理性检验（QC）

云端系统会对时间序列执行三项核心 QC：

• 范围检验 ：校验数据是否在物理区间内（如风速 ，温度 ）。
• 趋势检验 ：校验数据变幅。如风速或风向连续 6 小时 无任何变化（标准差为 0），判定为传感器死线，该段数据赋 NaN 丢弃。
• 相关性检验 ：核对相邻层高的风速差值是否合理，检查并剔除高低层异常倒挂的记录。

2. 数据修正与相关性插补

对于检验出的数据空洞（ NaN ），平台会自动采用线性回归算法修复：

• 优先选择同一测风塔其他有效高度的通道作为自变量 ，待插补层作为因变量 ，利用 最小二乘法 进行一元回归：

  （要求两个通道同期数据相关系数 ）。

• 本塔通道不可用时，系统会自动关联空间邻近测风塔的相同高度数据进行水平插补。

3. 长系列代表性订正（MCP）

系统将对现场实测完整年数据进行历史订正（MCP）：

• 关联参证气象站不低于 30 年 的历史长序列，利用 方差比法（Variance Ratio） 或 线性回归 反推代表年的风况。方差比法公式如下：

  （要求实测同期扇区相关系数 ）。

4. 轮毂高度风特征参数计算

系统基于代表年长系列计算以下微观特征：

• 风切变指数（）拟合 ：利用幂律公式对不同高度的实测大风段风速进行拟合求得 指数。
• 威布尔参数（Weibull parameters）拟合 ：采用极大似然估计（MLE）拟合形状参数 和尺度参数 ，确定风速概率密度分布。
• 特征湍流强度计算 ：计算每个风速区间内（以 步长）的湍流均值与标准差，通过下式求得特征值：

5. 50 年一遇极端风速计算

平台调取气象站历年最大风速记录，拟合 耿贝尔分布（Gumbel 分布） ：

设定重现期 （对应累积概率 ）反解得出 年一遇极大/最大风速，从而评估风电机组选型的安全等级。

---

八、第十二步：风资源分布图绘制及评估报告编制（流程终点）

在云端平台电算化结束并导出所有特征参数（各高度代表年平均风速、风功率密度、Weibull 参数、风向频率及风能方向分布、风切变、湍流和极端大风等）后，最终阶段的工作包含：

1. 资源微观分布图绘制 ：结合风资源计算结果（如利用 WAsP/WindSim 在流场中外推）与地形图，绘制轮毂高度风速及风功率密度的平面等值线分布图（微观选址分布图）。
2. 定级核对 ：根据规范对照表（如表 5.3.6-2）确定风资源等级（1 级至 7 级，或 D-1 至 D-3 级），评估风能可利用价值。
3. 机型选择与可研报告 ：结合生产商提供的机组性能曲线，在确定微观选型后，开展发电量折减及安全性核算。

---

九、风能利用理论与风能利用系数（Cp）计算基础

在进行风机微观选址、功率曲线拟合以及年发电量（AEP）敏感性分析时，风能利用系数 Cp 与风力发电机的空气动力学特性是评估机组发电效率的核心基础。

9.1 风能利用系数（Cp）的物理定义与公式

风能利用系数 Cp （Power Coefficient）表示风力发电机旋转的风轮从流过的风能中吸收并转化为机械轴功率的效率：

其中：

• —— 风轮实际获得的轴功率，单位为瓦特（W）；
• —— 空气密度，单位为千克/立方米（kg/m³）；
• —— 风轮扫风面积，即 （其中 为叶轮直径，单位为 m）；
• —— 上游来流的实际风速，单位为米/秒（m/s）。

9.2 贝兹极限（Betz Limit）与实际转换效率

• 贝兹极限 ：根据 1919 年德国物理学家贝兹（Albert Betz）基于一维动量理论的严格证明，任何风力发电机在气流不发生旋转的理想假设下，能够从风能中提取的最大理论极限效率为：

即任何风机最多只能转化大约 59.3% 的流经风能，这是风能利用的理论物理天花板。

• 实际机组效率 ：在实际工程中，受叶片三维流动损失、叶尖涡流损失、桨叶表面摩擦阻力、尾流旋转以及传动机构效率等因素的影响，实际机组的最大 Cp 值通常在 0.40 至 0.50 之间（目前主流先进机组的最大设计值约为 0.48 至 0.49 ）。

9.3 Cp 与叶尖速比（）及桨叶节距角（）的非线性关系

在实际运行中， Cp 的大小并不是一个常数，而是随着机组的运行状态（叶片转速、变桨角度和风速）动态变化的。它是 叶尖速比 （）和 桨叶节距角 （）的多元高度非线性函数。

1. 叶尖速比（Tip Speed Ratio, ）

叶尖速比表示风轮叶片尖端的线速度与上游未受干扰的来流风速之比：

其中：

• —— 风轮转动角速度，单位为弧度/秒（rad/s）；
• —— 风轮半径，单位为米（m）；
• —— 风轮发电机转速，单位为转/分钟（r/min）；
• —— 实际来流风速，单位为米/秒（m/s）。

2. 桨叶节距角（Pitch Angle, ）

桨叶节距角是叶片弦线与风轮旋转平面的夹角（单位为度， ）。

• 在额定风速以下，为了最大化捕获风能，通常将节距角保持在最佳设计值（一般为 0° 或微小的负角），此时系统主要通过调节风轮转速以保持最佳叶尖速比 ，使 Cp 处于最高点。
• 在额定风速以上，为了防止风轮和发电机过载，机组控制系统启动 变桨控制 ，主动增大节距角 （使叶片“顺风”），故意降低 Cp 效率，使发电机输出功率稳定在额定功率。

3. 典型数学拟合公式（仿真与建模常用）

在风电机组的仿真与控制系统设计中，常用以下半经验指数公式来模拟 Cp 的动态特性：

其中，内部变量 的表达式为：

系数常数通常取值为：

通过该公式，可以绘制出一条以 为横坐标、以 Cp 为纵坐标、以 为参变量的 Cp-λ-β 效率特性曲面 。微观选址软件（如 WT、WindPRO）也是基于这一特性来解算机组在过渡风速段的精细出力的。

---

← 返回目录 · ← 返回深度研究