基于Simulink搭建的双馈风力发电(DFIG)模型

1、双馈风力发电机

  双馈异步风力发电机(DFIG)是一种绕线式感应发电机,是变速恒频风力发电机组的核心部分。DFIG主要由两大部分组成,分别是冷却系统和本体:本体的组成是定子、转子以及轴承,而冷却系统则分为空空冷、水冷和空水冷三类结构。
  DFIG为何称为“双馈”发电机:
  (1)发电机其中的定子和转子侧都参与到了电网输电的过程中;
  (2)实现发电机的有功功率以及无功功率进行单独的调节,对电网实施无功补偿;
  (3)DFIG同时囊括了异步和同步发电机的特点。

2、双馈异步风力发电机的数学模型

  DFIG风电系统由电气子系统以及机械子系统两大核心构成。机械子系统是由齿轮箱、风机以及部分组件组成的;而电气子系统由网侧和机侧这两部分组成。网侧包含网侧变换器、电网等,机侧包含DFIG发电机以及机侧变换器。DFIG系统结构如下图所示:

unity3d风力发电机 风力发电机模型怎么做_数学模型


  首先在对DFIG进行建模之前,我们需先完成DFIG正方向惯例的参数设定, 本文选取电压降落的方向作为电流的正方向。为了达到数学模型的简化,可以假定以下几点:

  (1)不考虑空间谐波的情况,设发电机定、转子三相绕组对称分布杂在空间上,它们之间相差120电角度,气隙磁场呈正弦规律分布。

  (2)不考虑铁耗以及磁路的非线性,假定绕组的自感和互感是与磁路工作点有关的恒定值,也不考虑频率和温度变化对DFIG风电系统参数的影响。

  (3)以定子侧为参考,将转子的参数均归算至定子侧,定、转子绕组归算后的有用匝数相同。

  DFIG仿真模型如下所示:

unity3d风力发电机 风力发电机模型怎么做_数学模型_02


  带转子电压补偿的双馈发电机控制器模型:

unity3d风力发电机 风力发电机模型怎么做_unity3d风力发电机_03


  风力发电模型:

unity3d风力发电机 风力发电机模型怎么做_Simulink_04


  DFIG+Wind组合模型如下图所示:

unity3d风力发电机 风力发电机模型怎么做_Simulink_05

3、仿真分析

3.1、清空变量

clear; clc;

3.2、DFIG参数设置

f_rated = 60;         % 额定频率 Hz
w_syn = 2*pi*f_rated; % 同步电转速 rad/s
V_ll_rated = 690;     % 额定线电压 V
p = 6;                % 极点数
s = 0.01;             % 滑移率
J = 70;               % 转动惯量 kg*m^2
R_s = 2e-3;           % 定子下标s Ohm
R_r = 1.5e-3;         % 转子下标r Ohm
X_ls = 50e-3;         
X_lr = 47e-3;
X_m = 860e-3;         
L_ls = X_ls/w_syn;    % H
L_lr = X_lr/w_syn;
L_m  = X_m /w_syn;
L_s  = L_ls+L_m;
L_r  = L_lr+L_m;
tau_r = L_r/R_r;      % 转子绕组时间常数
sigma = 1-L_m^2/(L_s*L_r);

3.3、DFIG初始条件设置

% 额定(满)负载下的转子转速
w_mech_rated = (1-s)*w_syn/(p/2); 
% 额定定子电流 A
I_s_rated = V_ll_rated/sqrt(3) / (R_s + j*X_ls + j*X_m*(R_r/s+j*X_lr)/(j*X_m+R_r/s+j*X_lr));
% 额定转子电流 A
I_r_rated = -I_s_rated*j*X_m/(j*X_m+R_r/s+j*X_lr);
% 额定扭矩
T_em_rated = (3*abs(I_r_rated)^2*R_r*(1-s)/s) / w_mech_rated;
T_load_rated = T_em_rated;
% dq 域中的电压
V_sd_rated = V_ll_rated;
V_rd_rated = 0;
V_rq_rated = 0;
% dq 域中的定子电流
I_sd_rated = sqrt(3)*real(I_s_rated);
I_sq_rated = sqrt(3)*imag(I_s_rated);
% dq 域中的转子电流
I_rd_rated = sqrt(3)*real(I_r_rated);
I_rq_rated = sqrt(3)*imag(I_r_rated);
% Rated Qs
Q_s_rated = (V_sd_rated)^2/(w_syn*L_s)+(L_m/L_s)*V_sd_rated*I_rq_rated;
% 额定流量
fl_sd_rated = L_s*I_sd_rated + L_m*I_rd_rated;
fl_sq_rated = L_s*I_sq_rated + L_m*I_rq_rated;
fl_rd_rated = L_m*I_sd_rated + L_r*I_rd_rated;
fl_rq_rated = L_m*I_sq_rated + L_r*I_rq_rated;

3.4、风机参数设置

A = 3904;            % 面积 m^2
R = 70.5/2;          % 转子半径 m
J_turb = 2.4*10^6;   % 转动惯量 kg * m^2
rho = 1.2;           % 空气密度
v_wind = [12; 9; 6]; % 风速 m/s

3.5、结果分析

unity3d风力发电机 风力发电机模型怎么做_matlab_06


unity3d风力发电机 风力发电机模型怎么做_Simulink_07


unity3d风力发电机 风力发电机模型怎么做_数学模型_08


  由上图可以看出,搭建的Simulink模型能很好的显示出DFIG风电机组工作在理想电网电压工况状态下。