1. 項(xiàng)目概述同步電機(jī)與構(gòu)網(wǎng)型變流器的頻率穩(wěn)定性研究在電力系統(tǒng)向高比例可再生能源轉(zhuǎn)型的背景下構(gòu)網(wǎng)型變流器Grid-Forming Converter正逐步取代傳統(tǒng)同步電機(jī)成為電網(wǎng)頻率支撐的主力。這個(gè)項(xiàng)目通過Simulink搭建了同步電機(jī)與構(gòu)網(wǎng)型變流器的混合系統(tǒng)模型研究兩者在頻率動(dòng)態(tài)響應(yīng)中的交互機(jī)制。我將在本文分享完整的建模方法、關(guān)鍵參數(shù)設(shè)計(jì)原則以及通過Matlab代碼實(shí)現(xiàn)的自動(dòng)化仿真流程。提示本文所有模型文件與代碼已打包文末提供下載方式。建議使用Matlab R2021b及以上版本運(yùn)行。2. 核心模型構(gòu)建與參數(shù)設(shè)計(jì)2.1 同步電機(jī)建模要點(diǎn)同步電機(jī)的Simulink模型采用經(jīng)典六階模型關(guān)鍵參數(shù)包括慣性時(shí)間常數(shù)H4.5s典型火力發(fā)電機(jī)組取值阻尼系數(shù)D2.0 pud/q軸同步電抗Xd1.8 pu, Xq1.7 pu勵(lì)磁系統(tǒng)采用IEEE Type-AC1A標(biāo)準(zhǔn)模型% 同步電機(jī)參數(shù)初始化代碼示例 syn_machine sprintf([ H 4.5; D 2.0; Xd 1.8; Xq 1.7;\n... Tdo 8.0; Tqo 0.4; Xd 0.3; Xq 0.3;]);2.2 構(gòu)網(wǎng)型變流器控制策略構(gòu)網(wǎng)型變流器采用虛擬同步機(jī)(VSG)控制架構(gòu)核心控制環(huán)路由以下部分組成虛擬慣量環(huán)節(jié)模擬同步機(jī)轉(zhuǎn)子運(yùn)動(dòng)方程J_vsg 0.5 kg·m2等效慣性時(shí)間常數(shù)約2sD_vsg 15 N·m·s/rad電壓電流雙環(huán)控制電流環(huán)帶寬500Hz電壓環(huán)帶寬50Hz功率同步環(huán)下垂系數(shù)Rp0.05 pu濾波時(shí)間常數(shù)Tp0.1s% VSG控制參數(shù)設(shè)置 J_vsg 0.5; % 虛擬轉(zhuǎn)動(dòng)慣量 D_vsg 15; % 虛擬阻尼系數(shù) Rp 0.05; % 有功下垂系數(shù) Tp 0.1; % 功率測(cè)量濾波器時(shí)間常數(shù)3. 系統(tǒng)級(jí)仿真搭建3.1 主電路拓?fù)湓O(shè)計(jì)系統(tǒng)采用如圖1所示的測(cè)試架構(gòu)同步電機(jī)額定容量100MVA構(gòu)網(wǎng)型變流器容量50MVA公共連接點(diǎn)(PCC)電壓230kV負(fù)載擾動(dòng)設(shè)置階躍變化±10%負(fù)荷3.2 仿真步長與求解器選擇為保證數(shù)值穩(wěn)定性與仿真效率變流器開關(guān)頻率2kHz → 仿真步長取1e-6s機(jī)電暫態(tài)過程采用ode23tb變步長求解器總仿真時(shí)間10s含1s初始化注意在Model Configuration Parameters中需勾選Algebraic loop選項(xiàng)以處理控制環(huán)路的代數(shù)環(huán)問題。4. 頻率穩(wěn)定性分析案例4.1 負(fù)荷階躍擾動(dòng)測(cè)試設(shè)置0.5s時(shí)突增10%負(fù)荷觀測(cè)系統(tǒng)頻率響應(yīng)純同步機(jī)系統(tǒng)最大頻差-0.35Hz穩(wěn)定時(shí)間8.2s混合系統(tǒng)最大頻差-0.28Hz穩(wěn)定時(shí)間5.6s純VSG系統(tǒng)最大頻差-0.15Hz穩(wěn)定時(shí)間3.1s% 頻率偏差計(jì)算代碼 [peaks, locs] findpeaks(abs(freq_deviation)); max_deviation max(peaks); settling_time time(find(abs(freq_deviation)0.01,1)) - 0.5;4.2 慣量支撐效果對(duì)比通過改變VSG虛擬慣量參數(shù)J_vsg觀測(cè)系統(tǒng)頻率變化率(RoCoF)J_vsg [kg·m2]初始RoCoF [Hz/s]最大頻差 [Hz]0.1-1.25-0.420.5-0.83-0.281.0-0.61-0.215. 自動(dòng)化仿真流程實(shí)現(xiàn)5.1 批量仿真腳本設(shè)計(jì)通過Matlab腳本實(shí)現(xiàn)參數(shù)掃描與結(jié)果自動(dòng)導(dǎo)出% 批量仿真示例 param_range linspace(0.1, 1.0, 10); % J_vsg參數(shù)掃描范圍 results cell(length(param_range),3); for i 1:length(param_range) set_param(GFM_Model/VSG_Control/J_vsg, Value, num2str(param_range(i))); simout sim(GFM_Model); % 結(jié)果提取 results{i,1} param_range(i); results{i,2} max(abs(simout.freq_deviation.Data)); results{i,3} simout.settling_time.Data(end); end5.2 結(jié)果可視化技巧使用Matlab圖形對(duì)象句柄實(shí)現(xiàn)專業(yè)級(jí)繪圖figure(Position, [100 100 800 400]) subplot(1,2,1) plot(time, freq_deviation, LineWidth,1.5) grid on; xlabel(Time (s)); ylabel(Frequency Deviation (Hz)) subplot(1,2,2) bar(cell2mat(results(:,1)), cell2mat(results(:,2))) xlabel(Virtual Inertia J_{vsg}); ylabel(Max Frequency Deviation)6. 工程經(jīng)驗(yàn)與問題排查6.1 常見仿真報(bào)錯(cuò)處理代數(shù)環(huán)問題現(xiàn)象仿真報(bào)錯(cuò)Algebraic loop contains...解決方案在反饋環(huán)路中加入單位延遲(1/z)模塊數(shù)值振蕩現(xiàn)象高頻小幅振蕩調(diào)整方法減小仿真步長或增加濾波器時(shí)間常數(shù)初始化失敗檢查所有Powergui模塊配置確認(rèn)電機(jī)初始滑差設(shè)置為06.2 模型加速技巧將連續(xù)系統(tǒng)模塊替換為離散版本使用Accelerator模式運(yùn)行對(duì)非關(guān)鍵子系統(tǒng)啟用Atomic Subsystem7. 模型文件與代碼下載完整項(xiàng)目包含主仿真模型 GFM_Model.slx參數(shù)初始化腳本 init_parameters.m批量仿真程序 batch_simulation.m結(jié)果處理工具 result_analysis.m點(diǎn)擊下載項(xiàng)目文件包我在實(shí)際調(diào)試中發(fā)現(xiàn)當(dāng)VSG容量超過系統(tǒng)總?cè)萘?0%時(shí)需要特別注意電壓環(huán)參數(shù)調(diào)整否則會(huì)導(dǎo)致PCC點(diǎn)電壓波動(dòng)加劇。建議按照IEC 62909標(biāo)準(zhǔn)進(jìn)行阻抗穩(wěn)定性校驗(yàn)后再確定最終控制參數(shù)。