CN105978023B - 基于波形松弛法含dfig电力***暂态稳定域边界求取方法 - Google Patents
基于波形松弛法含dfig电力***暂态稳定域边界求取方法 Download PDFInfo
- Publication number
- CN105978023B CN105978023B CN201610293078.6A CN201610293078A CN105978023B CN 105978023 B CN105978023 B CN 105978023B CN 201610293078 A CN201610293078 A CN 201610293078A CN 105978023 B CN105978023 B CN 105978023B
- Authority
- CN
- China
- Prior art keywords
- dfig
- driven generator
- wind
- power
- waveform relaxation
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 230000001052 transient effect Effects 0.000 title claims abstract description 15
- 238000011217 control strategy Methods 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 6
- 239000013598 vector Substances 0.000 claims description 5
- 238000005312 nonlinear dynamic Methods 0.000 claims description 2
- 238000000465 moulding Methods 0.000 claims 1
- 238000004422 calculation algorithm Methods 0.000 abstract description 8
- 238000004364 calculation method Methods 0.000 abstract description 5
- 230000005611 electricity Effects 0.000 abstract description 2
- 238000005183 dynamical system Methods 0.000 description 8
- 230000000750 progressive effect Effects 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 2
- 230000006641 stabilisation Effects 0.000 description 2
- 238000011105 stabilization Methods 0.000 description 2
- 230000000739 chaotic effect Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- H02J3/386—
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- 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
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Control Of Eletrric Generators (AREA)
- Wind Motors (AREA)
Abstract
本发明属于电力***暂态稳定域领域,考虑含双馈式异步风力发电机(DFIG)电力***在定功率控制策略运行时利用波形松弛算法刻画其稳定域边界。该方法针对风电场接入电网以后其稳定域计算的问题,依据非线性微分动力***稳定域计算的波形松弛算法,提出了用波形松弛算法计算双馈式风力发电机接入电网以后其稳定域计算的新算法,以双馈式风力发电机的单质量块轴系模型建立了目标函数,以这个目标函数作为非线性微分动力***,利用波形松弛算法的雅可比迭代格式来计算稳定域。本发明提供了一种计算DFIG接入电力***其在运行过程中的暂态稳定域计算的方法,为未来分布式风电大量接入电网其稳定域计算提供了相应的指导依据。
Description
技术领域
本发明属于电力***暂态稳定域领域,考虑含双馈式异步风力发电机(DFIG)电力***在定功率控制策略运行时利用波形松弛算法刻画其稳定域边界。
背景技术
风能是一种可再生能源,在能源日益枯竭环境日益污染严重的今天,特别是火力发电量占了我国总发电量的很大一部分比例,煤炭属于一次性能源,终有一天会面临着枯竭的局面,所以发展风能这样的清洁能源显得十分的有必要。风力发电在世界范围内取得了迅速发展。变速风力机***是一个非线性自治动力***,变速机的轴系特性、桨叶特性、桨距角变化、风速波动以及输出电磁功率的变化均会变速风力机的运行状态改变,这些方面是变速风力机的自身因素影响其稳定性。与变速风力机的自身的稳定性不同,变速风力机能否稳定还是受到变速风力机的控制策略、发电机输出功率的变化以及风速的突变都是引起变速风力机失稳的原因。目前风电机组的有功功率控制策略主要有两种:1、最大风功率捕获控制(maximum power point tracking),2、定功率控制。本专利选取的是变速风力机的转速和桨距角为状态变量,分析不同控制策略下变速风力机转速和桨距角的轨迹。MPPT控制的目标是调节变速风力机转速到最优转速,使风力机捕获风能的效率最高。变速风力机输出电磁功率大小的设定是为了配合调节变速风力机的转速,因此在MPPT运行过程中,变速风力机不会出现失去稳定的现象。定功率控制的目标是利用变速风力机的转动惯量,使风电机输出的电磁功率保持恒定,风电机组在运行过程中按照设定的目标功率输出电磁功率,变速风力机的转子将释放动能已维持输出电磁功率恒定;若变速风力机在运行过程中持续出现转子释放动能的状态,变速风力机将会出现失去稳定,因此在变速风力机状态空间中存在稳定的边界,若变速风力机的状态在边界以内,变速风力机将趋于稳定。
依据非线性动力***理论,稳定域边界应由其所有不稳定平衡点的稳定流形的并集构成,因此电力***的稳定域边界是由其稳定边界上所有一型不稳定平衡点处稳定流形的闭包之并集构成。波形松弛方法是微分方程数值解法,它将整个微分***分成若干个子***。各自独立进行仿真计算,子***之间的耦合关系通过子***中的一些波形信号的相互传递来实现,在***间联系较弱的情况下,比较容易得到收敛结果。通过对微分动力***的稳定域进行综合分析,得到微分动力***的一个渐进稳定平衡点稳定域的边界是由闭轨线或由不稳定平衡点以及该不稳定的稳定域流形构成。
发明内容
针对现有技术中计算含DFIG电力***暂态稳定域的方法中存在的问题,有必要提出一种新的计算含DFIG电力***暂态稳定域的新的方法。本发明专利构建DFIG的二阶模型,提出了基于波形松弛法含DFIG电力***暂态稳定域的计算。首先,依据变速风力发电机单质量块轴系模型建立微分动力***,写出微分动力***的雅可比迭代形式,以桨距角和转速为变量,找出微分动力***的所有平衡点,确定渐进稳定平衡点稳定域边界上的平衡点;其次,用微分动力***的反方向***确定原***渐进平衡点稳定域边界上平衡点的稳定流形;最后,渐进稳定平衡点稳定域的边界是由边界上的平衡点和该平衡点稳定流形的并集构成,从而计算变速风力机的稳定运行域。
为实现上述目的,本发明采取以下技术方案:
本发明解决上述问题采取的技术方案:
步骤1、变速风力发电机建模,模型的变量是DFIG的桨距角和转速,本发明主要是对风力发电机采用单质量块轴系模型,建立二阶轴系模型,得到***的渐进稳定平衡点xs,否者记为并求出在***在平衡点处的雅可比矩阵,求出雅可比矩阵具有正实部特征值且具有单位长度的不稳定特征向量,记为yi,求出这些具有单位长度的不稳定特征向量yi和平衡点的ε球边界的每一个交截(交截点是)。
步骤2、利用波形松弛法计算暂态稳定域边界。
从步骤(1)的得到的每个交截点开始,在DFIG模型反方向***用波形松弛方法计算积分轨线确定,求出原***渐进稳定平衡点稳定域边界上平衡点的稳定流形
步骤3、最后变速风力发电机的暂态稳定运行域是由步骤2确定的边界平衡点以及该平衡点点稳定流形构成。
附图说明
图1是本发明的算法总体流程图;
图2变速风力机的稳定域
具体实施方式
本发明包括以下步骤:
1、变速风力发电机建模以及确定波形松弛法的模型。
考虑到DFIG风力机单质量块模型中,其转子运动方程为:
变速风力机桨距角的运动方程为:
变速风力机捕获风能产生的机械转矩如下:
在以上三个式子中wr为变速风力机转速,J为变速风力机的转动惯量,Tm为变速风力机的机械转矩,Te为变速风力机输出的电磁功率,β桨距角,β0为变速风力机桨距角参考值,τβ为桨距角调节时间常数,Cp为风能捕获系数,ρ为空气密度,R为叶轮半径,Vwind为风速。变速风力机的风能捕获系数是风速、变速风力机转速和桨距角的函数,即:其中λ为叶尖速比,
DFIG的二阶轴系模型可表示为:
这个模型以变速风力机的桨距角β和转速wr为变量。这个二阶模型的雅可比迭代模型如下:
反方向***的雅可比迭代模型如下:
2、利用波形松弛算法计算暂态稳定域
2.1、找出微分***的所有平衡点
由单质量块轴系模型构成了一个所需的微分动力***,找出这个微分动力***的所有的平衡点,即Tm-Te=0与βo-β=0的所有解,并确定渐进平衡点稳定域边界上的平衡点。
2.2计算稳定流形
(1)利用DFIG单质量块模型,确定***的所有平衡点,零型平衡点记为xs,否则记为求出***(1-5)在处的雅可比矩阵。
(2)找出雅可比矩阵具有正实部特征值且具有单位长度的不稳定特征向量yi。
(3)找出这些具有单位长度的不稳定特征向量yi和平衡点的每一个交截和
(4)从这些交截点开始对***(1-6)利用波形松弛法计算积分轨线到某一步,如果曲线保留在这个ε球内,则进行下一步。否则用αε代替ε,相应的重复步骤(3)
(5)从这些交截点开始利用波形松弛方法计算积分轨线。
2.3渐进平衡稳定点xs稳定域的边界是由第二步确定的平衡点和该平衡点稳定流形的并集构成。渐进平衡稳定点的稳定域边界即时本发明所要求的含DFIG电力***暂态稳定域的边界。如图二所示,虚线左方为DFIG运行不稳定的区域,右方为DFIG运行的稳定区域。
Claims (2)
1.基于波形松弛法含DFIG电力***暂态稳定域边界求取方法,包括如下步骤:
1)、DFIG风力发电机建模;
2)、利用波形松弛法计算暂态稳定域边界;
在考虑DFIG风力发电机的单质量块模型中,选取DFIG风力发电机的转速和浆距角为状态变量,考虑定功率控制策略下,利用DFIG风力发电机的转动惯量,使DFIG风力发电机输出的电磁功率保持恒定,DFIG风力发电机在运行过程中按照设定的目标功率输出电磁功率,DFIG风力发电机的转子将释放动能已维持输出电磁功率恒定;
单质量块模型中,DFIG风力发电机的转子运动方程为:
DFIG风力发电机桨距角的运动方程为:
建立的DFIG风力发电机的模型如下:
建立模型的雅可比迭代造型:
建立模型的反方向***的波形松弛迭代模型:
2.根据权利要求l所述的基于波形松弛法含DFIG电力***暂态稳定域边界求取方法,其特征在于:所述步骤2)中依据非线性动力***理论,稳定域边界应由其所有不平衡点的稳定流形的并集构成,根据单质量块模型,得到***的渐进稳定平衡点,并求出***在平衡点处的雅可比矩阵,求出雅可比矩阵具有正实部特征值且具有单位长度的不稳定特征向量,求出这些具有单位长度的不稳定特征向量和平衡点的每一个交截点,在得到的每个交截点处开始利用单质量块模型的反方向***用波形松弛法计算积分轨线;最后DFIG风力发电机的稳定运行域是由以上的边界平衡点以及该平衡点的稳定流形构成。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610293078.6A CN105978023B (zh) | 2016-05-05 | 2016-05-05 | 基于波形松弛法含dfig电力***暂态稳定域边界求取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610293078.6A CN105978023B (zh) | 2016-05-05 | 2016-05-05 | 基于波形松弛法含dfig电力***暂态稳定域边界求取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105978023A CN105978023A (zh) | 2016-09-28 |
CN105978023B true CN105978023B (zh) | 2018-12-18 |
Family
ID=56992459
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610293078.6A Active CN105978023B (zh) | 2016-05-05 | 2016-05-05 | 基于波形松弛法含dfig电力***暂态稳定域边界求取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105978023B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104917201A (zh) * | 2015-06-16 | 2015-09-16 | 山东大学 | 模拟惯性与超速相结合的双馈风机有功频率控制器及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5952208B2 (ja) * | 2013-02-28 | 2016-07-13 | 東京瓦斯株式会社 | 給電システム、給電制御プログラムおよび給電制御方法 |
-
2016
- 2016-05-05 CN CN201610293078.6A patent/CN105978023B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104917201A (zh) * | 2015-06-16 | 2015-09-16 | 山东大学 | 模拟惯性与超速相结合的双馈风机有功频率控制器及方法 |
Non-Patent Citations (1)
Title |
---|
波形松弛法的电力***暂态稳定性并行仿真计算;林济铿等;《电工技术学报》;20061226;第21卷(第12期);第47-53和65页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105978023A (zh) | 2016-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kumar et al. | A review of conventional and advanced MPPT algorithms for wind energy systems | |
Shariatpanah et al. | A new model for PMSG-based wind turbine with yaw control | |
Musunuri et al. | Comprehensive review of wind energy maximum power extraction algorithms | |
Mansour et al. | A grid-connected variable-speed wind generator driving a fuzzy-controlled PMSG and associated to a flywheel energy storage system | |
Sabzevari et al. | MPPT control of wind turbines by direct adaptive fuzzy-PI controller and using ANN-PSO wind speed estimator | |
CN108488035B (zh) | 永磁直驱风力发电机组失速和变桨混合控制方法 | |
CN110206686A (zh) | 一种用于风力发电机组的自适应最大功率跟踪控制方法 | |
Mohammadi et al. | Flicker emission, voltage fluctuations, and mechanical loads for small-scale stall-and yaw-controlled wind turbines | |
CN104675629A (zh) | 一种变速风力发电机组的最大风能捕获方法 | |
CN102254092A (zh) | 一种含双馈风电机组的大型风电场动态等值方法 | |
Sarkar et al. | A study of MPPT schemes in PMSG based wind turbine system | |
Ciccarelli et al. | Sliding mode control based robust observer of aerodynamic torque for variable-speed wind turbines | |
Hosseini et al. | Control of pitch angle in wind turbine based on doubly fed induction generator using fuzzy logic method | |
Fadaeinedjad et al. | Using STATCOM to mitigate voltage fluctuations due to aerodynamic aspects of wind turbines | |
CN101854064B (zh) | 双馈感应风力发电***基于能量的机网侧联合控制算法 | |
Zhang et al. | Intelligent control for large-scale variable speed variable pitch wind turbines | |
Xia et al. | A new pitch control method for large scale wind turbine based on ADRC | |
CN105978023B (zh) | 基于波形松弛法含dfig电力***暂态稳定域边界求取方法 | |
Barambones | A robust wind turbine control using a Neural Network based wind speed estimator | |
CN107255062B (zh) | 一种自适应密度的风电机组转矩-转速控制方法 | |
Peng | Nonlinear feed forward pitch controller for wind turbine based on rotor's aerodynamic characteristic | |
Mohapatra et al. | Fractional-order sliding mode control based MPPT of a DFIG connected to a Variable Speed Wind Turbine | |
CN104863793A (zh) | 一种根据平均值触发风力发电机变桨动作指令的控制方法 | |
Fadaeinedjad et al. | Simulation of a wind turbine with doubly-fed induction machine using FAST and Simulink | |
Pan et al. | Multi-mode simulation modeling of full-size variable speed pumped storage system including the complete unit links |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |