基于Simulink搭建的雙饋風力發電（DFIG）模型
1、雙饋風力發電機
雙饋異步風力發電機（DFIG）是一種繞線式感應發電機，是變速恒頻風力發電機組的核心部分。DFIG主要由兩大部分組成，分別是冷卻系統和本體：本體的組成是定子、轉子以及軸承，而冷卻系統則分為空空冷、水冷和空水冷三類結構。
DFIG為何稱為“雙饋”發電機：
（1）發電機其中的定子和轉子側都參與到了電網輸電的過程中；
（2）實現發電機的有功功率以及無功功率進行單獨的調節，對電網實施無功補償；
（3）DFIG同時囊括了異步和同步發電機的特點。
2、雙饋異步風力發電機的數學模型
DFIG風電系統由電氣子系統以及機械子系統兩大核心構成。機械子系統是由齒輪箱、風機以及部分組件組成的；而電氣子系統由網側和機側這兩部分組成。網側包含網側變換器、電網等，機側包含DFIG發電機以及機側變換器。DFIG系統結構如下圖所示：
首先在對DFIG進行建模之前，我們需先完成DFIG正方向慣例的參數設定，本文選取電壓降落的方向作為電流的正方向。為了達到數學模型的簡化，可以假定以下幾點：
（1）不考慮空間諧波的情況，設發電機定、轉子三相繞組對稱分布雜在空間上，它們之間相差120角度，氣隙磁場呈正弦規律分布。
（2）不考慮鐵耗以及磁路的非線性，假定繞組的自感和互感是與磁路工作點有關的恒定值，也不考慮頻率和溫度變化對DFIG風電系統參數的影響。
（3）以定子側為參考，將轉子的參數均歸算至定子側，定、轉子繞組歸算后的有用匝數相同。
DFIG仿真模型如下所示：
帶轉子電壓補償的雙饋發電機控制器模型：
風力發電模型：
DFIG+Wind組合模型如下圖所示：
3、仿真分析
%% 初始化
clear; clc;
%% 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（stator） Ohm
R_r = 1.5e-3; % 轉子下標r（rotor） 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);
%% 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;
%% 風機參數
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