CN108595861A - 基于psasp直驱风电机组简化建模及参数辨识方法 - Google Patents
基于psasp直驱风电机组简化建模及参数辨识方法 Download PDFInfo
- Publication number
- CN108595861A CN108595861A CN201810405039.XA CN201810405039A CN108595861A CN 108595861 A CN108595861 A CN 108595861A CN 201810405039 A CN201810405039 A CN 201810405039A CN 108595861 A CN108595861 A CN 108595861A
- Authority
- CN
- China
- Prior art keywords
- psasp
- wind
- protection
- voltage
- current
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000001052 transient effect Effects 0.000 claims abstract description 35
- 238000006243 chemical reaction Methods 0.000 claims abstract description 28
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 9
- 230000002068 genetic effect Effects 0.000 claims abstract description 9
- 238000004088 simulation Methods 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000005611 electricity Effects 0.000 claims description 5
- 230000009897 systematic effect Effects 0.000 claims description 4
- 230000009471 action Effects 0.000 claims description 3
- 230000006978 adaptation Effects 0.000 claims description 3
- 238000009434 installation Methods 0.000 claims description 3
- 230000014759 maintenance of location Effects 0.000 claims description 3
- 230000001681 protective effect Effects 0.000 claims description 3
- 238000011835 investigation Methods 0.000 claims description 2
- 230000000630 rising effect Effects 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 7
- 238000011084 recovery Methods 0.000 description 7
- 239000013589 supplement Substances 0.000 description 6
- 238000011160 research Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 108700039691 Genetic Promoter Regions Proteins 0.000 description 1
- 241000233805 Phoenix Species 0.000 description 1
- BQCADISMDOOEFD-UHFFFAOYSA-N Silver Chemical compound [Ag] BQCADISMDOOEFD-UHFFFAOYSA-N 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000008929 regeneration Effects 0.000 description 1
- 238000011069 regeneration method Methods 0.000 description 1
- 229910052709 silver Inorganic materials 0.000 description 1
- 239000004332 silver Substances 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Biophysics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Economics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- General Health & Medical Sciences (AREA)
- Computational Linguistics (AREA)
- Computer Hardware Design (AREA)
- Biomedical Technology (AREA)
- Genetics & Genomics (AREA)
- Data Mining & Analysis (AREA)
- Physiology (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Geometry (AREA)
- Artificial Intelligence (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明公开了一种基于PSASP直驱风电机组简化建模及参数辨识方法。本发明使用PSASP/UD模型搭建D‑PMSG稳态潮流模型和D‑PMSG暂态简化模型;所述的稳态潮流模型忽略发电机及变流器的损耗,其为将风速转换成风机输出机械功率的风速‑功率转换模块;所述的暂态简化模型包括功率转换模块、输出电流计算模块、网侧逆变器运行控制模块、故障穿越及保护控制模块和参数辨识模块。本发明使用PSASP/UD搭建D‑PMSG稳态和暂态简化模型,克服了PSASP自带的直驱永磁模型不能表述实际直驱永磁风电机组运行特性的问题,仿真结果与实测数据更加接近,并且在低压故障穿越及保护部分运行效果比现有模型更加理想,并用遗传算法实现了PSASP中重要参数的合理辨识。
Description
技术领域
本发明属于PSASP风机建模技术领域,具体地说是一种基于PSASP直驱风电机组简化建模及参数辨识方法。
背景技术
直驱永磁风机电磁暂态模型虽然可以精确反映风电机组运行特性,但是结构复杂、计算量大,不适合大规模风电场暂态稳定性研究,因此需要根据研究目标对暂态模型进行简化。
文献1:CONROY J,WATSON R.Aggregate modelling of wind farms containingfull-converter wind turbine generators with permanent magnet synchronousmachines:transient stability studies[J].IET Renewable Power Generation,2009,3(1):39-52;文献2:ELLIS A,KAZACHKOV Y,MULJADI E,et al.Description and technicalspecifications for generic WTG models a status report[C]//2011IEEE/PES PowerSystems Conference and Exposition(PSCE).Phoenix,Arizona,US,2011:1-8。文献[1-2]提出在暂态稳定分析的时间尺度内,PMSG的机侧元件没有动态响应,PMSG的网侧响应主要由机组的网侧变流器及其控制决定,因此可以不对机侧元件及其控制器详细建模,而是仅考虑它们的运行状态和控制效果。文献[1-2]提出将D-PMSG等效成受控电压源或受控电流源,为D-PMSG的简化工作提出了有效的思路。
文献3:徐力、王刚、侍乔明等人.直驱永磁风电机组电压暂态简化建模及仿真[J].电机与控制应用,2015,42(9):47-51;文献4:雷虹云、郑超、汪宁渤等人.基于变频器受控源模拟的双馈风机等效仿真[J].中国电力,2012,45(6):87-91。文献[3-4]以文献[1-2]提出的简化模型基础上,根据风电机组低压穿越过程中的跌落及无功动态响应对机组进一步简化,通过模拟低压穿越时风电机组的动态特性验证了简化模型的有效性。文献5:刘忠义、刘崇茹、李庚银.适用于暂态稳定分析的直驱永磁风机建模研究[J].电网与清洁能源,2016,32(2):96-102,则是将D-PMSG两种简化模型对比,得出在大规模风电集群并网时,受控电流源型PMSG比受控电压源模型更能保证暂态稳定分析正确性的结论。然而上述文献主要只考虑了发电机和机侧变流器的简化等效,在以下三个方面略有不足:
1.在风机简化过程中,忽略了风机转速达到转速额定值,但风速未达到风速额定值的风机恒转速运行情况。在实际工程中,这种情况普遍出现,需要考虑。
2.在低压穿越过程中,未能在无功限流环节的基础上考虑短路电流的设置,无法保证系统在低压过程中的稳定时;同时没有考虑有功在恢复过程出现的过有功现象对系统稳定性的影响。
3.系统参数的获取不够准确,未能有效模拟简化模型与实际系统之间的关系,从而证明D-PMSG简化模型的可行性。
发明内容
本发明针对上述现有文献的不足,重点关注适用于实际工程的D-PMSG建模方法,同时考虑到大规模直驱永磁风电系统机电暂态仿真时采用PSASP自带的直驱永磁模型不能表述实际直流输电特性的问题,研究适用于不同时间尺度的直驱永磁风机模型,利用PSASP自定义建模环境搭建直驱永磁风机模型。
为此,本发明采用如下的技术方案:基于PSASP直驱风电机组简化建模及参数辨识方法,其使用PSASP/UD模型搭建稳态潮流模型和D-PMSG暂态简化模型;
所述的稳态潮流模型忽略发电机及变流器的损耗,其为将风速转换成风机输出机械功率的风速-功率转换模块;所述的暂态简化模型包括功率转换模块、输出电流计算模块、网侧逆变器运行控制模块、故障穿越及保护控制模块和参数辨识模块;
用PSASP/UD模型搭建D-PMSG暂态简化模型的步骤包括:
1)根据功率转换模块获取电网有功和无功,同时完成PSASP/UD初始化设置,获取初始电流;
2)根据网侧逆变器运行控制模块完成有功无功解耦控制,并设置有功无功限流模式,正常情况下采用有功限流模式,低压穿越时采用无功限流模式;
3)获取dq坐标系下电流分量,经过故障穿越及保护控制模块,完成低压穿越故障情况下短路电流的设置及故障恢复过程的有功坡道设置,再根据故障时间设置频率和电压二级保护控制;
4)根据输出电流计算模块完成PSASP/UD到电网系统的交互工作,将dq坐标系下的电流转换到xy坐标系下,同时计算电网经过D-PMSG暂态简化模型后的实际功率,并作为反馈量反馈到网侧逆变器控制模块,进一步完善D-PMSG暂态简化模型中的模型参数。
5)为了达到系统参数辨识的目的,本发明按照PSASP/UD上的逻辑控制模块在另一个电力系统仿真软件Matlab/Simulink中搭建D-PMSG暂态简化模型,通过遗传算法(GA)完成参数辨识工作。
作为上述技术方案的补充,所述的稳态潮流模型研究不同风速和转速对输出功率的影响时,考虑MPPT最优控制理论,风机的输出功率为:
式中,P是风机输出的有功功率;V是风速;R为风力机叶片半径;Vstart-up是风机运行的起始风速;Kmppt是风机最优功率输出的比例系数;λopt是最优叶尖比;Cpmax是最大风能转换系数;Vs是当风机转速达到最大转速所对应的风速;Vref是额定风速;Cp为风能转换系数;ω*为额定转速;λ是尖速比;ρ是空气密度;θ为桨距角;P*为额定功率;β是中间变量。
作为上述技术方案的补充,所述的功率转换模块:在D-PMSG稳态潮流模型的基础上,增加一阶轴系方程获得风机转速,然后根据MPPT控制得到最大功率输出,再除以转速得到电磁转矩并反馈到轴系方程中,由此跟踪风力机最优转速,一阶轴系方程如下:
Jeqdωw/dt=Tm-Te-Bmωw,
其中:Tm、Te分别是风力机机械转矩和用磁发电机电磁转矩;Bm是风电机组的阻尼系数;Jeq为PMSG等效转动惯量,ωw是风机转速。
作为上述技术方案的补充,所述的输出电流计算模块用于完成UD模型到PSASP系统的交互工作,
初始化过程功率到电流的转换公式:
电流从dq坐标到xy坐标转换计算公式:
式中,ITR、ITI、ITR0、ITI0分别是xy坐标系下电流的实轴和虚轴分量及它们的初值;VT0、PG0、QG0分别是母线电压、有功和无功初值;VT是母线电压值;Angel和θ都是电压相角。
作为上述技术方案的补充,所述的网侧逆变器运行控制模块包括有功无功解耦控制及限流环节,在网侧扰动引起的电网低电压期间,有功无功的值由网侧输出电流和电网电压决定,在dq解耦下点功率方程为:
Ps=UdsIds
Qs=-UdsIqs,
式中,Ps、Uds、Ids、Qs、Iqs分别是dq坐标系下有功、d轴电压、d轴电流、无功、q轴电流;
由此看出,通过调节d轴电流分量Ids能直接控制网侧输出的有功,通过调节q轴电流分量Iqs能直接调节网侧输出的无功,从而实现解耦控制;
限流环节主要是对电路进行保护,共有两种模式,当系统正常运行时,系统有功对系统稳定性影响更大,采用有功优先模式;当低电压穿越时,系统无功的波动对系统影响更大,采用无功优先模式。
作为上述技术方案的补充,所述的故障穿越及保护控制模块中,低压穿越的过/欠电压、过/欠频保护策略为:
保护信号flag初始值为0,当下列任一事件发生时,flag为1,保护装置动作,不再输出电流,系统退出运行:若ULVP2≤U≤ULVP1且持续事件超过tLVP1,一级欠压保护;若U<ULVP2且持续时间超过tLVP2,二级欠压保护;若UHVP2≤U≤UHVP1且持续时间超过tHVP1时,一级过压保护;若U<UHVP2且持续时间超过tHVP2,二级过压保护;若fLf2≤f≤fLf1且持续时间超过tLf1时,一级过频保护;若f<fLf2且持续时间超过tLf2,二级过频保护;若fHf2≤f≤fHf1且持续时间超过tHf1时,一级过频保护;若f<fHf2且持续时间超过tHf2,二级过频保护;
U、f分别为母线电压和系统频率,ULVP2、ULVP1、tLVP1、tLVP2分别是欠压保护的上、下电压阀值和一、二级欠压保护的时间;UHVP2、UHVP1、tHVP1、tHVP2分别是过压保护的上、下电压阀值和一、二级过压保护的时间;fLf2、fLf1、tLf1、tLf2分别是欠频保护的上、下电压阀值和一、二级欠压保护的时间;fHf2、fHf1、tHf1、tHf2分别是过频保护的上、下电压阀值和一、二级过频保护的时间。
作为上述技术方案的补充,具体的参数辨识步骤如下:
1)省去限流环节等对参数识别结果影响不大的部分,在matlab/simulink中搭建D-PMSG暂态简化模型;
2)设定遗传算法的运行参数,初始化种群;
3)读取真实系统数据,输入激励信号,进行仿真计算;设置目标函数,目标函数为低压穿越时D-PMSG实测有功无功与当前有功无功的差值;根据输出误差计算个体适应度;
4)对种群进行轮盘赌、交叉、变异、精英保留策略,生成下一代群落;
5)判断终止条件是否满足,满足结束,不满足的跳转到步骤3)。
本发明具有的有益效果如下:本发明使用PSASP/UD搭建D-PMSG稳态和暂态简化模型,克服了PSASP自带的直驱永磁模型不能表述实际直驱永磁风电机组运行特性的问题,仿真结果与实测数据更加接近,并且在低压故障穿越及保护部分运行效果比现有模型更加理想,并用遗传算法实现了PSASP中重要参数的合理辨识。
附图说明
图1为本发明实施例中稳态风机转速-风速-输出功率特性曲线图;
图2为本发明实施例中暂态简化模型整体结构图(V、P、Q、Uterm分别为风速、风电场向电网输送的有功、风电场向电网输送的无功和母线电压;Pord、Qord为功率转换模型得到的有功和无功;Ip_cmd、Iq_cmd为网侧运行控制模块输出的电流直轴和交轴分量,Ip、Iq分别是故障穿越及保护控制模块输出的电流分量直轴和交轴分量);
图3为本发明实施例中使用PSASP/UD搭建暂态简化模型整体流程图;
图4为本发明实施例中网侧逆变器运行控制模块示意图(P、Q为电网有功和无功,Pord、Qord为功率转换模型得到的有功和无功;IP_cmd、Iq_cmd为网侧运行控制模块输出的电流直轴和交轴分量,Ipmax、Ipmin、Iqmax、Iqmin分别为直轴电流最大值、直轴电流最小值、交轴电流最大值、交轴电流最小值);
图5为本发明实施例中低压穿越及保护控制模块示意图(IP_cmd、Iq_cmd为网侧运行控制模块输出的电流直轴和交轴分量;Ip、Iq分别是故障穿越及保护控制模块输出的电流分量直轴和交轴分量;U、f、flag、k分别为电网电压、系统频率、保护信号和短路电流比例系数;Ip_FRT、Iq0_LV、Iqmax_LV分别为电网电压低于二级低压值时的直轴低压穿越电流设置值、交轴低压时无功电流初值、交轴低压时无功电流最大值;UHV1、ULV1分别是二级、一级低压保护阀值);
图6为本发明实施例中基于遗传算法参数辨识流程图;
图7为本发明应用例中三级九节点系统图(其中,GEN1-230、GEN2-230、GEN3-230、STNA-230、STNB-230、STNC-230都是母线名称);
图8为本发明应用例中matlab/simulink详细模型低压穿越测试图(图8a为有功低压穿越图,图8b为无功低压穿越图);
图9为本发明应用例中实测数据与PSASP/UD模型测试数据对比图(图9a为有功对比图,图9b为无功对比图);
图10为本发明应用例中不考虑低压穿越故障恢复坡道设置示意图。
具体实施方式
下面结合说明书附图和具体实施方式对本发明作进一步说明。
实施例
(1)本发明主要分为暂态模型和稳态模型的PSASP/UD搭建。如图1所示的稳态风机转速-风速-输出功率特性曲线图,主要是阐述D-PMSG在风机启动区、Cp恒定区、转速恒定区和功率恒定区时风速和转速与风机输出功率之间的关系,便于理解D-PMSG功率转换时的特点。
(2)图2主要表示的是使用构建D-PMSG暂态简化模型整体结构图,主要由功率转换模块、输出电流计算模块、网侧逆变器运行控制模块、故障穿越及保护控制模块、参数辨识模块5个模块组成。图3为使用PSASP/UD模型搭建D-PMSG暂态简化模型的系统流程图,实施步骤为:
1)根据功率转换模块获取电网有功和无功,同时完成PSASP/UD初始化设置获取初始电流;
2)根据网侧逆变器运行控制模块完成有功无功解耦控制,并设置有功无功限流模式,正常情况下采用有功限流模式,低压穿越时采用无功限流模式。具体的解耦过程和限流方法如图4所示。
3)获取dq坐标系下电流分量,经过故障穿越及保护控制模块,完成低压穿越故障情况下短路电流的设置及故障恢复过程的有功坡道设置,再根据故障时间设置频率和电压二级保护控制。具体的短路电流设置方法和故障保护策略如图5所示。
4)根据输出电流计算模块完成PSASP/UD到电网系统的交互工作,将dq坐标系下的电流转换到xy坐标系下,同时计算电网经过D-PMSG暂态简化模型后的实际功率,并作为反馈量反馈到网侧逆变器控制模块,进一步完善D-PMSG暂态简化模型中的模型参数。
5)为了达到系统参数辨识的目的,本发明按照PSASP/UD上的逻辑控制模块在另一个电力系统仿真软件Matlab/Simulink中搭建D-PMSG暂态简化模型,通过遗传算法(GA)完成参数辨识工作。详细的参数辨识流程图如图6所示。
下面对上述的五个模块进行详细说明。
1)功率转换模块:主要在潮流稳态模型的基础上,增加一阶轴系方程获得风机转速,然后根据MPPT控制得到最大功率输出,再除以转速得到电磁转矩并反馈到轴系方程中去,由此跟踪风力机最优转速,一阶轴系方程如下:
Jeqdωw/dt=Tm-Te-Bmωw,
其中:Tm、Te分别是风力机机械转矩和用磁发电机电磁转矩;Bm是机组的阻尼系数;Jeq为PMSG等效转动惯量;ωw是风机转速。
2)输出电流计算模块:主要完成UD模型到PSASP系统的交互工作。由于PSASP/UD模型只能通过修改xy轴坐标系下电流的实部和虚部来完成对系统的作用,因此在建立UD模型和PSASP系统联系时需要将功率量转换成电流量,在计算控制保护部分输出电流时需要将dq坐标系下电流转换到xy坐标系下。
初始化过程功率到电流的转换:
电流从dq坐标到xy坐标转换计算公式:
式中,ITR、ITI、ITR0、ITI0分别是xy坐标系下电流的实轴和虚轴分量及它们的初值;VT0、PG0、QG0分别是母线电压、有功和无功初值;VT是母线电压值;Angel和θ都是电压相角。
3)网侧逆变器运行控制模块:主要采用的是用受控电流源来简化变流器。D-PMSG与传统同步机不同,在暂态稳定分析的时间尺度内,网侧的扰动对机侧元件没有影响,D-PMSG的网侧响应主要由机组的网侧变流器及其控制决定。因此研究以网侧低电压穿越特性作为直驱风机暂态特性的问题,可以忽略机侧变流器及发电机,将变流器简化成受控电流源。
网侧逆变器运行控制模块主要由有功无功解耦控制及限流环节组成。在网侧扰动引起的电网低电压期间,有功无功的值由网侧输出电流和电网电压决定,在dq解耦下点功率方程为:
Ps=UdsIds
Qs=-UdsIqs,
式中,Ps、Uds、Ids、Qs、Iqs分别是dq坐标系下有功、d轴电压、d轴电流、无功、q轴电流;
由此可以看出,通过调节d轴电流分量Ids可以直接控制网侧输出的有功,通过调节q轴电流分量Iqs可以直接调节网侧输出的无功,从而实现解耦控制。
限流环节主要是为了防止出现过电流的情况,对电路进行保护,共有2种模式。当系统正常运行时,系统有功对系统稳定性影响更大,因此采用有功优先模式。当低电压穿越时,系统无功的波动对系统影响更大,采用无功优先模式。
式中,Iq_cmd为网侧运行控制模块输出的电流直轴和交轴分量,Imax、Ipmax、Ipmin、Iqmax、Iqmin分别为限制的最大电流、直轴电流最大值、直轴电流最小值、交轴电流最大值、交轴电流最小值。
4)低压穿越及故障保护模块:处理PMSG低电压穿越的方法是当运行点发生故障引起电压过低时,需要给系统提供一定的有功和无功支撑。在PSASP/UD中的表示为给有功电流和无功电流一个定值,防止系统震荡失稳。同时为了防止有功变化过快,形成过有功的情况,应在故障恢复过程中设置有功恢复速率,使得有功功率随着坡道逐步上升。
低压穿越的过/欠电压、过/欠频保护策略为:
保护信号flag初始值为0,当下列任一事件发生时,flag为1,保护装置动作,不再输出电流,系统退出运行:若ULVP2≤U≤ULVP1且持续事件超过tLVP1,一级欠压保护;若U<ULVP2且持续时间超过tLVP2,二级欠压保护;若UHVP2≤U≤UHVP1且持续时间超过tHVP1时,一级过压保护;若U<UHVP2且持续时间超过tHVP2,二级过压保护;若fLf2≤f≤fLf1且持续时间超过tLf1时,一级过频保护;若f<fLf2且持续时间超过tLf2,二级过频保护;若fHf2≤f≤fHf1且持续时间超过tHf1时,一级过频保护;若f<fHf2且持续时间超过tHf2,二级过频保护。
U、f分别为母线电压和系统频率,ULVP2、ULVP1、tLVP1、tLVP2分别是欠压保护的上、下电压阀值和一、二级欠压保护的时间;UHVP2、UHVP1、tHVP1、tHVP2分别是过压保护的上、下电压阀值和一、二级过压保护的时间;fLf2、fLf1、tLf1、tLf2分别是欠频保护的上、下电压阀值和一、二级欠压保护的时间;fHf2、fHf1、tHf1、tHf2分别是过频保护的上、下电压阀值和一、二级过频保护的时间。
5)参数辨识模块。具体的参数辨识步骤是:
1.省去限流环节等对参数识别结果影响不大的部分,在matlab/simulink中搭建D-PMSG简化模型。
2.设定GA(遗传算法)的运行参数,初始化种群。
3.读取真实系统数据,输入激励信号,进行仿真计算。设置目标函数,目标函数为低压穿越时D-PMSG实测有功无功与当前有功无功的差值。根据输出误差计算个体适应度。
4.对种群进行轮盘赌、交叉、变异、精英保留策略,生成下一代群落。
5.判断终止条件是否满足,满足的话结束,不满足的跳转到步骤3。
应用例
应用例采用如图7所示的三机九节点系统,首先基于matlab/simulink自带的D-PMSG详细模型,设置仿真2秒时发生三相短路故障,详细模型仿真波形如图8所示。图8显示无功也在故障恢复时突变到近0.3左右,大体上与测试数据吻合,而有功在低压穿越过后有功一度达到了0.24左右的过功率状态,过了0.6s后基本恢复到原有的系统状态。然后将近海风电场实测数据与本发明基于PSASP/UD搭建的简化模型测试数据做对比,对比图为图9,对比图中实测数据与简化模型测试波形图基本接近,证明了参数辨识工作确实很好的完成了预期的目标。最后去除故障穿越模式下短路电流限制和有功故障恢复坡道设计,图10显示不考虑低压穿越故障恢复坡道的有功示意图,从图中可以看出在故障恢复时产生了明显了过有功情况,这在实际工程中对系统稳定性有着很大的影响并有可能因此产生系统震荡,更加证明了增加故障恢复坡道设计的必要性。所以得出结论:本发明设计的基于PSASP/UD的D-PMSG简化模型正确合理,适用于实际工程的运用。
以上所述实施例仅表达了本发明的实施方式,并不能因此理解为对本发明保护范围的限制,也并非对本发明的结构作任何形式上的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进,这些都属于本发明的保护范围。
Claims (7)
1.基于PSASP直驱风电机组简化建模及参数辨识方法,其特征在于,使用PSASP/UD模型搭建D-PMSG稳态潮流模型和D-PMSG暂态简化模型;
所述的D-PMSG稳态潮流模型忽略发电机及变流器的损耗,其为将风速转换成风机输出机械功率的风速-功率转换模块;所述的暂态简化模型包括功率转换模块、输出电流计算模块、网侧逆变器运行控制模块、故障穿越及保护控制模块和参数辨识模块;
用PSASP/UD模型搭建D-PMSG暂态简化模型的步骤包括:
1)根据功率转换模块获取电网有功和无功,同时完成PSASP/UD初始化设置,获取初始电流;
2)根据网侧逆变器运行控制模块完成有功无功解耦控制,并设置有功无功限流模式,正常情况下采用有功限流模式,低压穿越时采用无功限流模式;
3)获取dq坐标系下电流分量,经过故障穿越及保护控制模块,完成低压穿越故障情况下短路电流的设置及故障恢复过程的有功坡道设置,再根据故障时间设置频率和电压二级保护控制;
4)根据输出电流计算模块完成PSASP/UD到电网系统的交互工作,将dq坐标系下的电流转换到xy坐标系下,同时计算电网经过D-PMSG暂态简化模型后的实际功率,并作为反馈量反馈到网侧逆变器控制模块,进一步完善D-PMSG暂态简化模型中的模型参数。
5)为了达到系统参数辨识的目的,本发明按照PSASP/UD上的逻辑控制模块在另一个电力系统仿真软件Matlab/Simulink中搭建D-PMSG暂态简化模型,通过遗传算法(GA)完成参数辨识工作。
2.根据权利要求1所述的基于PSASP直驱风电机组简化建模及参数辨识方法,其特征在于,所述的稳态潮流模型研究不同风速和转速对输出功率的影响时,考虑MPPT最优控制理论,风机的输出功率为:
式中,P是风机输出的有功功率;V是风速;R为风力机叶片半径;Vstart-up是风机运行的起始风速;Kmppt是风机最优功率输出的比例系数;λopt是最优叶尖比;Cpmax是最大风能转换系数;Vs是当风机转速达到最大转速所对应的风速;Vref是额定风速;Cp为风能转换系数;ω*为额定转速;λ是尖速比;ρ是空气密度;θ为桨距角;P*为额定功率;β是中间变量。
3.根据权利要求1或2所述的基于PSASP直驱风电机组简化建模及参数辨识方法,其特征在于,所述的功率转换模块:在D-PMSG稳态潮流模型的基础上,增加一阶轴系方程获得风机转速,然后根据MPPT控制得到最大功率输出,再除以转速得到电磁转矩并反馈到轴系方程中,由此跟踪风力机最优转速,一阶轴系方程如下:
Jeqdωw/dt=Tm-Te-Bmωw,
其中:Tm、Te分别是风力机机械转矩和用磁发电机电磁转矩;Bm是风电机组的阻尼系数;Jeq为PMSG等效转动惯量,ωw是风机转速。
4.根据权利要求1或2所述的基于PSASP直驱风电机组简化建模及参数辨识方法,其特征在于,所述的输出电流计算模块用于完成UD模型到PSASP系统的交互工作,
初始化过程功率到电流的转换公式:
电流从dq坐标到xy坐标转换计算公式:
式中,ITR、ITI、ITR0、ITI0分别是xy坐标系下电流的实轴和虚轴分量及它们的初值;VT0、PG0、QG0分别是母线电压、有功和无功初值;VT是母线电压值;Angel和θ都是电压相角。
5.根据权利要求1或2所述的基于PSASP直驱风电机组简化建模及参数辨识方法,其特征在于,所述的网侧逆变器运行控制模块包括有功无功解耦控制及限流环节,在网侧扰动引起的电网低电压期间,有功无功的值由网侧输出电流和电网电压决定,在dq解耦下点功率方程为:
式中,Ps、Uds、Ids、Qs、Iqs分别是dq坐标系下有功、d轴电压、d轴电流、无功、q轴电流;
由此看出,通过调节d轴电流分量Ids能直接控制网侧输出的有功,通过调节q轴电流分量Iqs能直接调节网侧输出的无功,从而实现解耦控制;
限流环节主要是对电路进行保护,共有两种模式,当系统正常运行时,系统有功对系统稳定性影响更大,采用有功优先模式;当低电压穿越时,系统无功的波动对系统影响更大,采用无功优先模式。
6.根据权利要求1或2所述的基于PSASP直驱风电机组简化建模及参数辨识方法,其特征在于,所述的故障穿越及保护控制模块中,低压穿越的过/欠电压、过/欠频保护策略为:
保护信号flag初始值为0,当下列任一事件发生时,flag为1,保护装置动作,不再输出电流,系统退出运行:若ULVP2≤U≤ULVP1且持续事件超过tLVP1,一级欠压保护;若U<ULVP2且持续时间超过tLVP2,二级欠压保护;若UHVP2≤U≤UHVP1且持续时间超过tHVP1时,一级过压保护;若U<UHVP2且持续时间超过tHVP2,二级过压保护;若fLf2≤f≤fLf1且持续时间超过tLf1时,一级过频保护;若f<fLf2且持续时间超过tLf2,二级过频保护;若fHf2≤f≤fHf1且持续时间超过tHf1时,一级过频保护;若f<fHf2且持续时间超过tHf2,二级过频保护;
U、f分别为母线电压和系统频率,ULVP2、ULVP1、tLVP1、tLVP2分别是欠压保护的上、下电压阀值和一、二级欠压保护的时间;UHVP2、UHVP1、tHVP1、tHVP2分别是过压保护的上、下电压阀值和一、二级过压保护的时间;fLf2、fLf1、tLf1、tLf2分别是欠频保护的上、下电压阀值和一、二级欠压保护的时间;fHf2、fHf1、tHf1、tHf2分别是过频保护的上、下电压阀值和一、二级过频保护的时间。
7.根据权利要求1或2所述的基于PSASP直驱风电机组简化建模及参数辨识方法,其特征在于,具体的参数辨识步骤如下:
1)省去限流环节等对参数识别结果影响不大的部分,在matlab/simulink中搭建D-PMSG暂态简化模型;
2)设定遗传算法的运行参数,初始化种群;
3)读取真实系统数据,输入激励信号,进行仿真计算;设置目标函数,目标函数为低压穿越时D-PMSG实测有功无功与当前有功无功的差值;根据输出误差计算个体适应度;
4)对种群进行轮盘赌、交叉、变异、精英保留策略,生成下一代群落;
5)判断终止条件是否满足,满足结束,不满足的跳转到步骤3)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810405039.XA CN108595861B (zh) | 2018-04-28 | 2018-04-28 | 基于psasp直驱风电机组简化建模及参数辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810405039.XA CN108595861B (zh) | 2018-04-28 | 2018-04-28 | 基于psasp直驱风电机组简化建模及参数辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108595861A true CN108595861A (zh) | 2018-09-28 |
CN108595861B CN108595861B (zh) | 2019-08-20 |
Family
ID=63619406
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810405039.XA Active CN108595861B (zh) | 2018-04-28 | 2018-04-28 | 基于psasp直驱风电机组简化建模及参数辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108595861B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109241664A (zh) * | 2018-09-29 | 2019-01-18 | 贵州电网有限责任公司 | 一种关于风力发电机雷电电磁暂态特性的分析计算方法 |
CN110165703A (zh) * | 2019-05-29 | 2019-08-23 | 国网内蒙古东部电力有限公司 | 一种永磁直驱风机电磁暂态仿真的初始化方法及系统 |
CN110188466A (zh) * | 2019-05-30 | 2019-08-30 | 国网内蒙古东部电力有限公司 | 一种双馈风机电磁暂态仿真的初始化方法及系统 |
CN111725840A (zh) * | 2020-06-29 | 2020-09-29 | 浙江大学 | 一种直驱风电机组控制器参数辨识方法 |
CN111884267A (zh) * | 2020-08-12 | 2020-11-03 | 华北电力大学(保定) | 一种d-pmsg机组频率自适应附加惯性控制方法 |
CN112327217A (zh) * | 2020-10-13 | 2021-02-05 | 上海电力大学 | 基于变压器稳态和暂态运行的绕组漏电感在线辨识方法 |
CN112886634A (zh) * | 2021-01-14 | 2021-06-01 | 国网浙江省电力有限公司电力科学研究院 | 一种风电场中长期控制建模方法以及等效模型 |
CN112910006A (zh) * | 2021-03-19 | 2021-06-04 | 哈尔滨工业大学 | 一种直驱风电机组通用电磁暂态建模方法 |
CN112968464A (zh) * | 2019-12-13 | 2021-06-15 | 新疆金风科技股份有限公司 | 永磁直驱风力发电机组高低电压连续穿越控制方法及系统 |
CN114123336A (zh) * | 2021-11-30 | 2022-03-01 | 国网山东省电力公司德州供电公司 | 一种直驱风机变流器的控制参数辨识方法及系统 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102592026A (zh) * | 2012-01-12 | 2012-07-18 | 甘肃省电力公司风电技术中心 | 一种双馈风机等效模拟的仿真建模方法 |
CN102664427A (zh) * | 2012-05-25 | 2012-09-12 | 华北电力大学(保定) | 永磁直驱风电机组低电压穿越时有功和无功协调控制方法 |
CN102999675A (zh) * | 2012-12-12 | 2013-03-27 | 上海市电力公司 | 双馈变速恒频风电机组系统电磁暂态仿真方法 |
CN103034761A (zh) * | 2012-12-12 | 2013-04-10 | 上海市电力公司 | 双馈变速恒频风电机组系统机电暂态仿真方法 |
CN103034764A (zh) * | 2012-12-12 | 2013-04-10 | 上海市电力公司 | 双馈变速恒频风电机组系统建模与仿真方法 |
CN103094921A (zh) * | 2013-01-06 | 2013-05-08 | 宁夏电力公司电力科学研究院 | 直驱式风力发电机组机电暂态建模方法及机电暂态模型 |
CN103955572A (zh) * | 2014-04-23 | 2014-07-30 | 国家电网公司 | 一种双馈式风力发电机组机电暂态模型的建模方法 |
CN105591402A (zh) * | 2016-01-28 | 2016-05-18 | 云南电网有限责任公司电力科学研究院 | 一种直驱永磁风电机组的建模与仿真方法和装置 |
CN107453389A (zh) * | 2016-05-31 | 2017-12-08 | 中国电力科学研究院 | 一种组串式光伏逆变器机电暂态模型建模方法 |
-
2018
- 2018-04-28 CN CN201810405039.XA patent/CN108595861B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102592026A (zh) * | 2012-01-12 | 2012-07-18 | 甘肃省电力公司风电技术中心 | 一种双馈风机等效模拟的仿真建模方法 |
CN102664427A (zh) * | 2012-05-25 | 2012-09-12 | 华北电力大学(保定) | 永磁直驱风电机组低电压穿越时有功和无功协调控制方法 |
CN102999675A (zh) * | 2012-12-12 | 2013-03-27 | 上海市电力公司 | 双馈变速恒频风电机组系统电磁暂态仿真方法 |
CN103034761A (zh) * | 2012-12-12 | 2013-04-10 | 上海市电力公司 | 双馈变速恒频风电机组系统机电暂态仿真方法 |
CN103034764A (zh) * | 2012-12-12 | 2013-04-10 | 上海市电力公司 | 双馈变速恒频风电机组系统建模与仿真方法 |
CN103094921A (zh) * | 2013-01-06 | 2013-05-08 | 宁夏电力公司电力科学研究院 | 直驱式风力发电机组机电暂态建模方法及机电暂态模型 |
CN103955572A (zh) * | 2014-04-23 | 2014-07-30 | 国家电网公司 | 一种双馈式风力发电机组机电暂态模型的建模方法 |
CN105591402A (zh) * | 2016-01-28 | 2016-05-18 | 云南电网有限责任公司电力科学研究院 | 一种直驱永磁风电机组的建模与仿真方法和装置 |
CN107453389A (zh) * | 2016-05-31 | 2017-12-08 | 中国电力科学研究院 | 一种组串式光伏逆变器机电暂态模型建模方法 |
Non-Patent Citations (8)
Title |
---|
SHUHUI LI 等: "Optimal and Direct-Current Vector Control of Direct-Driven PMSG Wind Turbines", 《IEEE TRANSACTIONS ON POWER ELECTRONICS》 * |
夏玥 等: "基于直驱式永磁同步发电机组的风电场动态建模", 《电网技术》 * |
徐力 等: "直驱永磁风电机组电压暂态简化建模及仿真", 《电机与控制应用》 * |
朱自强 等: "基于PSASP保护自定义建模的电网连锁故障仿真", 《合肥工业大学学报自然科学版》 * |
李和明 等: "永磁直驱风电机组低电压穿越时的有功和无功协调控制", 《电工技术学报》 * |
杨越 等: "基于PSCAD的永磁直驱风电机组简化模型研究", 《水电能源科学》 * |
王慧 等: "基于PSO算法的定速风电机组三质块传动系统模型参数辨识", 《河海大学学报(自然科学版)》 * |
许德操 等: "基于PSASP的直流系统用户自定义建模", 《电力系统电动化》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109241664A (zh) * | 2018-09-29 | 2019-01-18 | 贵州电网有限责任公司 | 一种关于风力发电机雷电电磁暂态特性的分析计算方法 |
CN109241664B (zh) * | 2018-09-29 | 2021-05-25 | 贵州电网有限责任公司 | 一种关于风力发电机雷电电磁暂态特性的分析计算方法 |
CN110165703B (zh) * | 2019-05-29 | 2020-09-29 | 国网内蒙古东部电力有限公司 | 一种永磁直驱风机电磁暂态仿真的初始化方法及系统 |
CN110165703A (zh) * | 2019-05-29 | 2019-08-23 | 国网内蒙古东部电力有限公司 | 一种永磁直驱风机电磁暂态仿真的初始化方法及系统 |
CN110188466A (zh) * | 2019-05-30 | 2019-08-30 | 国网内蒙古东部电力有限公司 | 一种双馈风机电磁暂态仿真的初始化方法及系统 |
CN110188466B (zh) * | 2019-05-30 | 2020-09-29 | 国网内蒙古东部电力有限公司 | 一种双馈风机电磁暂态仿真的初始化方法及系统 |
CN112968464A (zh) * | 2019-12-13 | 2021-06-15 | 新疆金风科技股份有限公司 | 永磁直驱风力发电机组高低电压连续穿越控制方法及系统 |
US11677345B2 (en) | 2019-12-13 | 2023-06-13 | Xinjiang Goldwind Science & Technology Co., Ltd. | Control method and system for continuous high and low voltage ride through of permanent-magnet direct-drive wind-driven generator set |
CN112968464B (zh) * | 2019-12-13 | 2022-12-09 | 新疆金风科技股份有限公司 | 永磁直驱风力发电机组高低电压连续穿越控制方法及系统 |
CN111725840A (zh) * | 2020-06-29 | 2020-09-29 | 浙江大学 | 一种直驱风电机组控制器参数辨识方法 |
CN111725840B (zh) * | 2020-06-29 | 2021-09-21 | 浙江大学 | 一种直驱风电机组控制器参数辨识方法 |
CN111884267A (zh) * | 2020-08-12 | 2020-11-03 | 华北电力大学(保定) | 一种d-pmsg机组频率自适应附加惯性控制方法 |
CN112327217B (zh) * | 2020-10-13 | 2023-02-10 | 上海电力大学 | 基于变压器稳态和暂态运行的绕组漏电感在线辨识方法 |
CN112327217A (zh) * | 2020-10-13 | 2021-02-05 | 上海电力大学 | 基于变压器稳态和暂态运行的绕组漏电感在线辨识方法 |
CN112886634A (zh) * | 2021-01-14 | 2021-06-01 | 国网浙江省电力有限公司电力科学研究院 | 一种风电场中长期控制建模方法以及等效模型 |
CN112910006A (zh) * | 2021-03-19 | 2021-06-04 | 哈尔滨工业大学 | 一种直驱风电机组通用电磁暂态建模方法 |
CN114123336A (zh) * | 2021-11-30 | 2022-03-01 | 国网山东省电力公司德州供电公司 | 一种直驱风机变流器的控制参数辨识方法及系统 |
CN114123336B (zh) * | 2021-11-30 | 2023-06-16 | 国网山东省电力公司德州供电公司 | 一种直驱风机变流器的控制参数辨识方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN108595861B (zh) | 2019-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108595861B (zh) | 基于psasp直驱风电机组简化建模及参数辨识方法 | |
CN102799722B (zh) | 一种风电场低电压穿越能力仿真验证方法 | |
Conroy et al. | Aggregate modelling of wind farms containing full-converter wind turbine generators with permanent magnet synchronous machines: transient stability studies | |
Fadaeinedjad et al. | The impact of tower shadow, yaw error, and wind shears on power quality in a wind–diesel system | |
Wang et al. | Comprehensive modeling and parameter identification of wind farms based on wide-area measurement systems | |
CN103558768B (zh) | 一种基于风电场内风速分布特性的等值建模方法 | |
CN103955572A (zh) | 一种双馈式风力发电机组机电暂态模型的建模方法 | |
Mutlu et al. | Power quality analysis of wind farm connected to Alaçatı substation in Turkey | |
Li et al. | Dynamic aggregation of doubly-fed induction generators (DFIGs) for stability analysis of wind power systems | |
Li et al. | DFIG-based wind farm electromagnetic dynamic model and impact on protection relay of transmission network | |
CN116090175B (zh) | 一种用于新能源配电网系统的等值建模求解方法及装置 | |
Altimania | Modeling of doubly-fed induction generators connected to distribution system based on eMEGASim® real-time digital simulator | |
Seyedi | Evaluation of the DFIG wind turbine built-in model in PSS/E | |
Qi et al. | Dynamic equivalent modeling of direct-drive PMSG wind farms based on the transient active power response characteristics | |
Li et al. | Structure Preserving Aggregation Method for Doubly-Fed Induction Generators in Wind Power Conversion | |
Karthik et al. | Steady-state initialization of doubly fed induction generator based wind turbine considering grid side filter | |
Gevorgian et al. | Wgrid-49 GMLC project report: Understanding the role of short-term energy storage and large motor loads for active power controls by wind power | |
Keung et al. | Generic models of wind turbines for power system stability studies | |
Kanabar et al. | Rotor speed stability analysis of constant speed wind turbine generators | |
Li et al. | Wind farm electromagnetic dynamic model and outgoing line protection relay RTDS testing | |
Ye et al. | Simulation and Analysis of PMSG-based Wind Energy Conversion System using Different Coverter Models | |
Ju et al. | Simplified modeling of Directly Driven Wind Turbine with Permanent Magnet Synchronous Generator based on PSASP/UD | |
Mi et al. | Study on voltage stability of power grid with large scale wind farm interconnected | |
Kabashi et al. | Wind Grid Code requirements regarding connection and operation of Wind Turbine in Kosovo | |
Rosyadi et al. | Novel control design of variable speed PM wind generator considering grid code requirement |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |