基于Simulink的风电场全系统建模与并网仿真平台

资 源 简 介

本项目利用MATLAB及其电力系统仿真工具箱（SimPowerSystems）构建高保真的风电场仿真模型，重点针对目前主流的双馈感应风力发电机（DFIG）或永磁直驱风力发电机（PMSG）进行深度建模与分析。系统功能包括但不限于：构建模拟自然界真实风况的风速模型（包含基本风、阵风、渐变风及随机噪声风）；建立包含空气动力学特性、传动系统及桨距角控制的风力机模型，实现宽风速范围内的最大功率点跟踪（MPPT）；设计基于电压定向矢量控制（VOC）的转子侧与网侧变流器控制策略，实现有功功率与无功功率的解耦控制及直流母线电压的稳定调节；模拟多机组汇集并网过程，分析风电场对电网电能质量的影响；设置电网故障场景（如三相短路、电压跌落），引入Crowbar保护电路或斩波电阻，仿真验证风电机组的低电压穿越（LVRT）能力及故障恢复特性。该项目旨在为风电并网控制策略优化、电网稳定性评估及故障特性分析提供可靠的数字化仿真环境。

详 情 说 明

风电场全系统建模与并网仿真平台 (DFIG机组)

项目介绍

本项目是一个基于MATLAB纯代码实现的风电场全系统仿真平台。不同于传统的Simulink图形化建模，本项目主要通过编写MATLAB脚本（m文件），利用离散化Euler法手动求解微分方程，从而实现对双馈感应风力发电机（DFIG）机组的深度建模与仿真。

该平台不仅包含高保真的物理模型（空气动力学、传动系统、电机电气动态），还完整实现了矢量控制策略（RSC/GSC）以及电网故障下的低电压穿越（LVRT）逻辑。项目旨在为风电并网控制策略研究、MPPT算法验证及故障特性分析提供一个透明、可编程的数字化仿真环境，特别适用于需要深入理解底层算法逻辑的科研与教学场景。

功能特性

• 纯代码仿真引擎：摆脱Simulink模块库依赖，通过离散时间步长（100us）迭代求解全系统状态，运行效率高且算法逻辑完全透明。
• 多维度风速模拟：内置组合风速模型，涵盖基本恒定风、阵风（余弦函数拟合）及随机白噪声风，模拟真实自然环境。
• 高精度气动模型：实现基于叶尖速比（$lambda$）和桨距角（$beta$）的非线性$C_p$曲线计算，动态求解气动转矩。
• 双馈电机（DFIG）电气建模：基于dq坐标系下的定子磁链定向模型，包含定子/转子磁链、电流及电磁转矩的动态计算。
• 完整的矢量控制系统：

* 机侧变流器（RSC）：实现最大功率点跟踪（MPPT）及有功/无功解耦控制。 * 网侧变流器（GSC）：负责维持直流母线电压稳定。

• 电网故障与保护仿真：支持模拟电网电压跌落（三相短路场景），内嵌Crowbar（撬棒）保护逻辑及LVRT期间的无功电流支撑策略。
• 数据记录与分析：自动记录全过程的风速、功率、转速、电压电流及故障状态，便于后续THD分析及波形绘图。

系统要求

• MATLAB：推荐 R2016a 及以上版本。
• 工具箱：本项目核心逻辑基于基础MATLAB编写，无需额外安装SimPowerSystems或其他专用工具箱即可运行主程序逻辑。

使用方法

1. 直接运行主程序函数。
2. 程序将初始化全局参数（1.5MW DFIG机组数据）。
3. 执行时间步长为 $100mu s$、总时长为 3.0s 的时域仿真。
4. 仿真过程中会自动计算并更新各个子系统的状态变量。
5. 仿真结束后，所有关键数据将保存在 Data 结构体中，可直接用于绘图或分析。

核心算法与代码实现逻辑

主程序采用单一入口函数结构，并在内部通过大循环进行时步迭代。以下是各模块的实际实现细节：

1. 全局参数配置

代码定义了基于1.5MW DFIG机组的详细物理参数：

• 空气动力学：空气密度 $1.225 kg/m^3$，叶轮半径 38.5m。
• 电机参数：极对数 $p=2$，定转子电阻、漏感及互感参数均已折算并定义。
• 控制参数：预设了桨距角、RSC（电流环、转速环、功率环）及GSC（电压环、电流环）的PI控制增益。

2. 复杂风速模型构建

在仿真循环中，程序实时合成风速信号 $v_{wind}$：

• 基础风：设定为 $11.0 m/s$。
• 阵风模型：在仿真时间 $t in [1.0, 1.5]$ 秒期间，叠加余弦形式的阵风分量。
• 随机风：全时段叠加高斯白噪声，模拟风速的随机波动。

3. 空气动力学与传动系统

• 计算叶尖速比（Lambda）：根据实时机械转速和风速计算。
• $C_p$ 曲线拟合：采用经典经验公式，利用指数函数和多项式组合，根据 Lambda 和桨距角 Beta 实时计算风能利用系数 $C_p$。
• 转矩计算：基于 $C_p$ 和风速的三次方计算气动转矩 $T_{aero}$，作为系统输入的机械功率源。

4. 电网模型与故障逻辑

• 电网电压生成：采用Thevenin等效电路，基于理想PLL生成电网相位。
• 故障模拟（LVRT）：

* 在 $1.5s$ 至 $1.65s$ 期间，模拟电网电压跌落故障，跌落深度由 Fault.DipDepth (20%) 控制。 * Crowbar 逻辑：当检测到故障（电压跌落）时，自动置位 Fault.Crowbar 标志位，用于闭锁RSC部分控制功能或改变控制策略。

5. DFIG 电气动态模型 (dq坐标系)

• 磁链-电流解耦：利用定子自感、转子自感及互感构成的电感矩阵行列式，将磁链状态变量（$psi_{ds}, psi_{qs}, psi_{dr}, psi_{qr}$）实时转换为定转子电流（$i_{ds}, i_{qs}, i_{dr}, i_{qr}$）。
• 电磁转矩：基于 $i_{qs}i_{dr} - i_{ds}i_{qr}$ 的叉乘关系计算实时电磁转矩 $T_e$。

6. RSC 矢量控制策略 (定子磁链定向)

控制器采用双闭环结构：

• 外环（转速/功率）：

* 计算MPPT参考转速：$omega_{ref} = v_{wind} times 7.0 / R$。 * 正常模式：速度误差经过PI调节生成转矩电流指令 $i_{qr}^*$，控制有功功率。 * 故障模式（LVRT）：根据电压跌落深度，注入无功电流 $Q_{ref}$ 以支撑电网电压，优先满足无功需求。

• 内环（电流）：

* 包含前馈解耦项，补偿由转速差（$w_{slip}$）引起的耦合电压。 * PI调节器输出转子电压参考值 $v_{dr}, v_{qr}$。 * 实施变流器最大输出电压限幅。

7. GSC 控制策略 (电压定向)

• 直流母线电压控制：GSC的外环负责维持 $V_{dc}$ 稳定在 1150V。
• 通过PI控制器根据母线电压误差计算网侧d轴电流指令 $i_{dg}^*$，实现有功功率平衡。

8. 数据采集

程序在每一步仿真中不仅更新状态，还将以下关键数据存入数组以便分析：

• 风速与 $C_p$ 值
• 电磁转矩与机械转速
• 电网三相电压与电流
• 直流母线电压
• 故障状态标志位