CN106451543B - 一种风电机组低电压穿越输出特性灵活仿真方法 - Google Patents
一种风电机组低电压穿越输出特性灵活仿真方法 Download PDFInfo
- Publication number
- CN106451543B CN106451543B CN201610936477.XA CN201610936477A CN106451543B CN 106451543 B CN106451543 B CN 106451543B CN 201610936477 A CN201610936477 A CN 201610936477A CN 106451543 B CN106451543 B CN 106451543B
- Authority
- CN
- China
- Prior art keywords
- active power
- active
- climb
- parameter
- value
- 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
- 230000009194 climbing Effects 0.000 claims description 21
- 238000013480 data collection Methods 0.000 claims description 4
- 238000012512 characterization method Methods 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims description 3
- 238000011084 recovery Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 2
- 230000000630 rising effect Effects 0.000 claims description 2
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 238000011217 control strategy Methods 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 238000000034 method Methods 0.000 abstract description 7
- 230000007547 defect Effects 0.000 abstract description 3
- 238000001514 detection method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 206010008190 Cerebrovascular accident Diseases 0.000 description 1
- 208000006011 Stroke Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000010354 integration 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
一种风电机组低电压穿越输出特性灵活仿真方法，属于电力***仿真技术领域，该方法基于双馈风机有功、无功独立控制，并且在有功控制模块中部分控制参数相对独立的特点，拟合出关键控制参数与低穿时期有功功率输出特性之间的线性对应关系，进一步利用任意机组在任意工况下的低电压穿越有功输出实测数据辨识出该控制参数，基于此可以完成对于其他任意工况下的低电压穿越有功功率输出特性的灵活仿真，本发明方法克服了传统的全工况低电压穿越测试耗时耗力的缺陷，是一种灵活性高、准确性好的风电机组低电压穿越输出特性仿真方法。
Description
技术领域
本发明属于电力***仿真技术领域，特别涉及一种风电机组低电压穿越输出特性灵活仿真方法。
背景技术
风电机组低电压穿越特性对并网电力***安全稳定运行有重要影响，因此低电压穿越特性测试是风机并网的必备检测项目。风电机组的全工况低电压穿越测试是对待检测风电机组最全面的测试评判，其中全工况测试指的是不同负荷状态、不同电压跌落程度下的低电压穿越能力检测。
对风电机组全面开展全工况低电压穿越测试存在以下三个缺点：1、单一风机需要进行多次测试才能确保其在各种工况下都具有低电压穿越能力，测试过程复杂、耗时长；2、每一次测试的负荷状态和电压跌落值都是离散的，只能得到测试点的低电压穿越情况；3、实际风电场中风机种类、型号众多，鉴于测试过程的复杂性，一般只能进行抽样测试，很难做到对全部机组的全覆盖测试。
但是，为了对有大型风电场并网的电力***进行安全性分析，必须对不同类型风机、不同运行工况下的低电压穿越输出特性进行模拟。
发明内容
为了克服传统全工况低电压穿越测试存在的缺陷，本发明的目的是提供一种风电机组低电压穿越输出特性灵活仿真方法，能够实现对不同工况下风电机组低电压穿越有功功率输出的灵活仿真。
为了达到上述目的，本发明的技术方案为：
一种风电机组低电压穿越输出特性灵活仿真方法，包括以下步骤：
步骤1：确定表征有功功率变化特性的特征量，包括：有功功率跌落深度P，有功功率快速爬升高度b和有功功率缓慢爬升速度a，
所述的有功功率跌落深度P，即有功功率从完全跌落到开始爬升的平均值P2与跌落前有功功率P1的比值，其中，完全跌落时刻指从该时刻起连续20个有功功率采样值与该时刻采样值之间的差值小于该采样值的10％；开始爬升时刻指从该时刻起连续20个有功功率采样值持续增加；
所述的有功功率快速爬升高度b，即有功功率在快速爬升期内的爬升量，其中，快速爬升期的时间起点指开始爬升的时刻，时间终点指从该时刻起连续20个有功功率采样值较其前一采样值的增量小于该时刻较前一采样值的增量；
所述的有功功率缓慢爬升速度a，即有功功率在缓慢爬升期内的爬升量Δa与爬升时间Δt的比值，其中，缓慢爬升期指的是快速爬升期结束至有功功率上升到95％P1的这一时期；
步骤2：对低电压穿越期间双馈风机有功功率控制第一阶段的策略进行修改，分别加入参数Km和Tramp将故障清除后阶段分为有功快速恢复和有功爬升两个阶段进行控制，并生成有功功率参考值以便后续生成相应的转子电流有功分量来实现对输出有功功率的调整；
低电压穿越期间的故障清除后时刻，首先在有功快速恢复阶段通过参数Km控制故障切除后有功快速上升的幅值：
Pref1＝K1(|u|-K2)+Km
然后在有功功率爬升阶段通过参数Tramp控制有功上升的速率：
式中：Pref1为有功功率参考值；u为故障期间发电机出口电压有效值；K1、K2为计算有功参考值时的调整参数，默认值分别为K1＝1.46，K2＝0.15；Pref1为故障穿越期间有功控制策略第一阶段的有功功率参考值；t为有功功率快速上升阶段结束时刻起到恢复到故障前数值所用时间；
步骤3：使用任意电力仿真软件，按照控制变量原则，改变控制参数Km和Tramp进行低电压穿越仿真，获得有功功率波形，提取有功功率波形中的三个特征量：有功功率跌落深度P、有功功率快速爬升高度b和有功功率缓慢爬升速度a，建立拟合数据集；
使用一阶多项式函数对控制参数Km与有功功率跌落深度P和有功功率快速爬升高度b之间的函数关系进行拟合；再对控制参数Tramp与有功功率跌落深度P、控制参数Km和有功功率缓慢上升速度a之间的函数关系进行拟合，即有：
式中：x1、x2、x3、y1、y2、y3分别为两个待拟合函数中的待定参数；
步骤4：提取实测有功功率波形的特征量，使用第三步中的公式进行辨识，得到控制参数Km和Tramp，利用Km和Tramp即可完成对于其他任意工况下的低电压穿越有功功率输出特性的灵活仿真。
本发明的特点及效果：
本发明方法快速准确地辨识出了风机低电压穿越有功功率控制部分的关键控制参数，实现了仿真模型在低电压穿越时期有功功率输出与实际机组较高程度的吻合，进而实现了对任意工况下风机低电压穿越有功输出特性的灵活仿真，大大提高了仿真的灵活性、准确度，弥补了现场低电压穿越测试的缺陷。
附图说明
图1为有功功率输出特征量提取示意图。
图2为实测数据与仿真数据对比图。
具体实施方式
本发明提出的风电机组低电压穿越输出特性灵活仿真方法，结合附图及实施例详细说明如下。
步骤1：参照图1，确定表征有功功率变化特性的特征量，包括：有功功率跌落深度P、有功功率快速爬升高度b、有功功率缓慢爬升速度a。
步骤2：对低电压穿越期间双馈风机有功功率控制第一阶段的策略进行修改，分别加入参数Km和Tramp将故障清除后阶段分为有功快速恢复和有功爬升两个阶段进行控制。
步骤3：按照控制变量原则，改变控制参数Km和Tramp进行低电压穿越仿真，获得有功功率波形。提取有功功率波形中的三个特征量：有功功率跌落深度P、有功功率快速爬升高度b和有功功率缓慢爬升速度a，建立拟合数据集如下：
表1拟合数据集
使用一阶多项式函数对控制参数Km与有功功率跌落深度P和有功功率快速爬升高度b之间的函数关系进行拟合；再对控制参数Tramp与有功功率跌落深度P、控制参数Km和有功功率缓慢上升速度a之间的函数关系进行拟合。即有：
Km＝f1(P,b)＝0.076+0.142×P+1.199×b
f1拟合结果的误差平方和小于2％，Km取不同值时f2拟合结果的误差平方和均小于0.1％，表明一阶多项式模型具有较高的拟合精度，可以满足参数辨识的需求。
步骤4：提取电压跌落20％情况下的有功实测波形的特征量，得到有功功率跌落深度P为0.089、有功功率快速爬升高度b为0.21、有功功率缓慢爬升速度a为0.6762。使用第三步中的公式进行辨识，将P和b带入到f1中可得控制参数Km为0.3268，选择f2中Km＝0.3时的对应的函数带入P和a得到控制参数Tramp为0.6603。利用得到的控制参数在仿真软件中进行低电压穿越的测试，将测试得到的有功曲线与实际的有功曲线进行对比，发现在误差允许的范围内两者匹配良好，达到了预期的参数辨识效果，具体请参见图2。因此，利用Km和Tramp即可完成对于其他任意工况下的低电压穿越有功功率输出特性的灵活仿真。
Claims (1)
1.一种风电机组低电压穿越输出特性灵活仿真方法，其特征在于，包括以下步骤：
步骤1：确定表征有功功率变化特性的特征量，包括：有功功率跌落深度P，有功功率快速爬升高度b和有功功率缓慢爬升速度a，
所述的有功功率跌落深度P，即有功功率从完全跌落到开始爬升的平均值P2与跌落前有功功率P1的比值，其中，完全跌落时刻指从该时刻起连续20个有功功率采样值与该时刻采样值之间的差值小于该时刻采样值的10％；开始爬升时刻指从该时刻起连续20个有功功率采样值持续增加；
所述的有功功率快速爬升高度b，即有功功率在快速爬升期内的爬升量，其中，快速爬升期的时间起点指开始爬升的时刻，时间终点指从该时刻起连续20个有功功率采样值较其前一采样值的增量小于该时刻较前一采样值的增量；
所述的有功功率缓慢爬升速度a，即有功功率在缓慢爬升期内的爬升量Δa与爬升时间Δt的比值，其中，缓慢爬升期指的是快速爬升期结束至有功功率上升到95％P1的这一时期；
步骤2：对低电压穿越期间双馈风机有功功率控制第一阶段的策略进行修改，分别加入参数Km和Tramp将故障清除后阶段分为有功快速恢复和有功爬升两个阶段进行控制，并生成有功功率参考值以便后续生成相应的转子电流有功分量来实现对输出有功功率的调整；
低电压穿越期间的故障清除后时刻，首先在有功快速恢复阶段通过参数Km控制故障切除后有功快速上升的幅值：
Pref1＝K1(|u|-K2)+Km
然后在有功功率爬升阶段通过参数Tramp控制有功上升的速率：
式中：u为故障期间发电机出口电压有效值；K1、K2为计算有功参考值时的调整参数，默认值分别为K1＝1.46，K2＝0.15；Pref1为故障穿越期间有功控制策略第一阶段的有功功率参考值；t为有功功率快速上升阶段结束时刻起到恢复到故障前数值所用时间；
步骤3：使用任意电力仿真软件，按照控制变量原则，改变控制参数Km和Tramp进行低电压穿越仿真，获得有功功率波形，提取有功功率波形中的三个特征量：有功功率跌落深度P、有功功率快速爬升高度b和有功功率缓慢爬升速度a，建立拟合数据集；
使用一阶多项式函数对控制参数Km与有功功率跌落深度P和有功功率快速爬升高度b之间的函数关系进行拟合；再对控制参数Tramp与有功功率跌落深度P、控制参数Km和有功功率缓慢上升速度a之间的函数关系进行拟合，即有：
式中：x1、x2、x3、y1、y2、y3分别为两个待拟合函数中的待定参数；
步骤4：提取实测有功功率波形的特征量，使用步骤3中的公式进行辨识，得到控制参数Km和Tramp，利用Km和Tramp即可完成对于其他任意工况下的低电压穿越有功功率输出特性的灵活仿真。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610936477.XA CN106451543B (zh) | 2016-11-01 | 2016-11-01 | 一种风电机组低电压穿越输出特性灵活仿真方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610936477.XA CN106451543B (zh) | 2016-11-01 | 2016-11-01 | 一种风电机组低电压穿越输出特性灵活仿真方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106451543A CN106451543A (zh) | 2017-02-22 |
| CN106451543B true CN106451543B (zh) | 2018-11-23 |
Family
ID=58178735
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610936477.XA Active CN106451543B (zh) | 2016-11-01 | 2016-11-01 | 一种风电机组低电压穿越输出特性灵活仿真方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106451543B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107302228B (zh) * | 2017-07-07 | 2019-11-08 | 中国石油大学(华东) | 一种基于dsp的低电压穿越区域通用判断方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101132107B1 (ko) * | 2010-09-29 | 2012-04-05 | 한국전력공사 | 분산전원이 연계된 전력계통의 전압/무효전력 제어 시스템 및 이를 위한 방법 |
| CN102510083A (zh) * | 2011-11-07 | 2012-06-20 | 清华大学 | 一种降低风电场低电压穿越容量需求的集成保护方法 |
| CN103606959A (zh) * | 2013-12-09 | 2014-02-26 | 华北电力大学 | 双馈型风力发电***低/高电压穿越的优化控制方法 |
| CN104242346A (zh) * | 2014-09-19 | 2014-12-24 | 许继集团有限公司 | 一种风电机组低电压穿越控制方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103141004B (zh) * | 2010-09-22 | 2016-12-07 | 东芝三菱电机产业***株式会社 | 电力转换装置 |
-
2016
- 2016-11-01 CN CN201610936477.XA patent/CN106451543B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101132107B1 (ko) * | 2010-09-29 | 2012-04-05 | 한국전력공사 | 분산전원이 연계된 전력계통의 전압/무효전력 제어 시스템 및 이를 위한 방법 |
| CN102510083A (zh) * | 2011-11-07 | 2012-06-20 | 清华大学 | 一种降低风电场低电压穿越容量需求的集成保护方法 |
| CN103606959A (zh) * | 2013-12-09 | 2014-02-26 | 华北电力大学 | 双馈型风力发电***低/高电压穿越的优化控制方法 |
| CN104242346A (zh) * | 2014-09-19 | 2014-12-24 | 许继集团有限公司 | 一种风电机组低电压穿越控制方法 |
Non-Patent Citations (2)
| Title |
|---|
| 基于实验数据辨识的双馈风力发电机组建模研究;王耀函;《中国优秀硕士学位论文全文数据库（电子期刊）工程科技Ⅱ辑》;20150215(第（2015年）02期);全文 * |
| 风电场电流保护动作特性等价低电压穿越容量性能分析;王宾等;《电力***自动化》;20111225;第35卷(第24期);第23-27页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106451543A (zh) | 2017-02-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107453401B (zh) | 一种双馈风力发电机参数辨识方法 | |
| CN102146812B (zh) | 电力***原动机及其调速器实测建模方法 | |
| CN106169747B (zh) | 一种双馈感应发电机参数辨识方法 | |
| CN106202815B (zh) | 基于有功响应的双馈型风电场动态等值方法 | |
| CN103093027B (zh) | 基于双馈风电场等值模型的电力***分析方法 | |
| CN103825521B (zh) | 一种双馈风电机组中驱动***及发电机参数的辨识方法 | |
| CN106058937B (zh) | 一种含直驱永磁风电场的配电网广义负荷建模方法 | |
| CN104392056B (zh) | 一种直驱风电机组建模方法及装置 | |
| CN103701389B (zh) | 一种双馈感应发电机控制器参数辨识方法 | |
| CN111293686A (zh) | 基于armax***辨识的电力***惯量实时评估方法 | |
| CN106383947B (zh) | 风电场集电网络动态等值参数的快速获取方法 | |
| CN102767472A (zh) | 一种风电机组功率曲线的提取方法 | |
| CN102830692B (zh) | 风电机组主控制***低电压穿越一致性测试方法 | |
| CN107947228B (zh) | 基于Markov理论的含风电电力***随机稳定性分析方法 | |
| CN105353304A (zh) | 一种风电机组电气模型低电压穿越特性验证方法 | |
| CN105119543A (zh) | 一种发电机组调速***远端线路甩负荷的控制方法及*** | |
| CN105134487A (zh) | 一种考虑湍流频率因素的风力机最大功率点跟踪控制方法 | |
| CN105576654B (zh) | 直驱风电场等值方法和*** | |
| CN104504285A (zh) | 一种计及撬棒保护的双馈风电场等值建模方法 | |
| CN107100795A (zh) | 一种低风速下风力发电机组mppt自适应控制方法 | |
| CN106410862B (zh) | 基于有功恢复斜率校正的风电场单机等值方法 | |
| CN106451543B (zh) | 一种风电机组低电压穿越输出特性灵活仿真方法 | |
| CN115392133A (zh) | 采用高斯混合模型的风电场最优聚类等值方法 | |
| CN107221933B (zh) | 一种概率潮流计算方法 | |
| CN101529347A (zh) | 识别有延时的调节对象的方法和装置、调节装置和计算机程序产品 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |