WRF模式中风力发电机拖曳参数化方案详解

WRF模式中风力发电机拖曳参数化方案详解

方案概述

WRF（Weather Research and Forecasting）模式中的风力发电机拖曳参数化方案是一种用于模拟风力发电场对大气边界层影响的先进方法。该方案由Fitch等人开发，基于Blahak等人的工作改进而来，能够更准确地表示风力涡轮机从大气中提取动能的过程。

核心原理

该参数化方案的关键创新点在于：

1. 使用涡轮机的推力系数（而非功率系数加估计损失因子）来计算从大气中提取的动能总量
2. 生成的湍流动能（TKE）是风速的函数
3. 通过表格形式引入功率曲线和推力系数，提高了方案的灵活性

使用方法

1. 激活参数化

在WRF的物理参数化namelist中设置：

windfarm_opt = 1

2. 准备风力涡轮机数据文件

需要创建windturbines.txt文件，每行定义一个涡轮机，包含三个参数：

1. 涡轮机纬度（度）
2. 涡轮机经度（度）
3. 涡轮机类型（整数）

示例：

55.574051 6.883480 1
55.569066 6.884697 1
30.000000 -77.000000 1

3. 涡轮机特性表

对于每种涡轮机类型（如类型1），需要创建对应的特性表文件wind-turbine-1.tbl，格式如下：

• 第一行：数据对数量N（整数）
• 第二行：4个涡轮机特性参数
  • 轮毂高度（米）
  • 转子直径（米）
  • 静态推力系数
  • 涡轮机额定功率（MW）
• 后续N行：风速、推力系数和功率输出（kW）数据对

示例：

23
80. 85. 0.1 2.0
3.    0.7     80.
4.    0.75   150.
...
25.   0.05   2000.0

技术要点

1. 分辨率建议：不建议在高于5倍转子直径的分辨率下运行此参数化方案，因为在这种高分辨率下需要考虑转子叶片的旋转效应（当前模型未实现）。

2. 湍流处理：该参数化与MYNN PBL方案配合使用，通过修改QKE场（代表2×TKE）来包含风力涡轮机产生的湍流动能。可以通过bl_mynn_tkeadvect开关控制QKE的平流。

3. 伪真实配置：测试模式下可以使用网格坐标（i,j）而非经纬度来定义涡轮机位置，需设置：

windfarm_ij = 1

并使用windturbines-ij.txt文件。

最佳实践建议

1. 数据来源：强烈建议使用涡轮机制造商提供的推力系数和功率曲线数据，而非使用模型自带的理想化数据。

2. 引用要求：使用该方案发表研究成果时，必须引用Fitch等人的原始论文（见下文参考文献）。

3. 应用范围：该方案特别适合研究中尺度风电场对局地和区域气候的影响。

参考文献

1. Fitch, A. C. et al. (2012): 关于风电场局地和中尺度影响的原始论文
2. Fitch, A. C. et al. (2013): 风电场在日循环周期中的中尺度影响研究
3. Fitch, A. C. et al. (2013): 气候模型中风电场参数化方法
4. Jimenez, P.A. et al.: 近海风电场尾流的中尺度建模研究

总结

WRF中的风力发电机拖曳参数化方案为研究风电场与大气相互作用提供了强大的工具。通过合理配置涡轮机参数和位置，研究人员可以准确模拟风电场对边界层动力学和能量平衡的影响，为风能资源评估和风电场规划提供科学依据。