CN105629168B - 一种基于功率区间的风电机组闪变确定方法 - Google Patents
一种基于功率区间的风电机组闪变确定方法 Download PDFInfo
- Publication number
- CN105629168B CN105629168B CN201410637797.6A CN201410637797A CN105629168B CN 105629168 B CN105629168 B CN 105629168B CN 201410637797 A CN201410637797 A CN 201410637797A CN 105629168 B CN105629168 B CN 105629168B
- Authority
- CN
- China
- Prior art keywords
- power
- flickering
- coefficient
- wind turbines
- power interval
- 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 31
- 238000012360 testing method Methods 0.000 claims abstract description 11
- 230000001186 cumulative effect Effects 0.000 claims description 12
- 230000005611 electricity Effects 0.000 claims description 7
- 230000003247 decreasing effect Effects 0.000 claims description 6
- 230000001052 transient effect Effects 0.000 claims description 5
- 230000007423 decrease Effects 0.000 claims description 3
- 210000003205 muscle Anatomy 0.000 claims 2
- 238000013459 approach Methods 0.000 abstract description 2
- 238000002474 experimental method Methods 0.000 abstract 1
- 238000005259 measurement Methods 0.000 description 11
- 238000004364 calculation method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000010977 unit operation Methods 0.000 description 1
Landscapes
- Wind Motors (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提供了一种基于功率区间的风电机组闪变确定方法，该方法包括以下步骤：I、确定测试点；II、根据风电机组的运行状态在不同的功率区间内进行数据采集；III、根据风电机组年平均输出功率和电网阻抗角计算闪变值。该方法采用功率分区间的方法进行风电机组的闪变计算，不需单独采集风速，有效的降低了试验成本，提高了数据的利用率，减少了试验过程中的不确定度来源。
Description
技术领域
本发明涉及一种新能源接入与控制领域的方法，具体讲涉及一种基于功率区间的风电机组闪变确定方法。
背景技术
风力发电是风能利用的重要形式，是可再生能源开发中技术最成熟的、最具有规模开发前景的发电方式之一。风能是可再生、无污染、能量大、前景广的能源，大力发展清洁能源是世界各国的战略选择。近年来我国风力发电发展势头迅猛，风电装机每年翻番，大规模并网风电机组的运行对电网及用户的影响也日趋重要。
理想的电力系统应以恒定的频率(50Hz)和正弦的波形，按照规定的电压水平对用户供电。在三相交流电力系统中，各相的电压和电流应处于幅值大小相等、相位互差120°的对称状态。但由于电力系统中的发电机、变压器、线路等系统元件的参数并不是理想线性或对称的，再加上目前调控手段的不完善、负荷性质各异、负荷变化的随机性、运行操作以及各种故障等原因，这种理想的对称状态实际上并不存在，因此就产生了电能质量的概念。电能质量的定义为：导致用电设备故障或不能正常工作的电压、电流或频率偏差，闪变就是电能质量的重要内容之一。
对于风电机组而言，由于风具有随机波动的特性，而风电机组的输出功率与风速的三次方成正比且与风向有关，因此其输出功率也是随机波动的。风电机组自身的固有特性，如塔影效应、偏航误差、风剪切作用和叶片重力不平衡等，也会对风电机组的输出功率造成影响，引起机组输出功率的周期性波动。这种不稳定的功率注入电网，尤其是在大规模风电并网的情况下，会对电网的稳定性和电能质量造成影响，引起电压的波动和闪变问题。因此，电压波动和闪变是风电机组引起的主要的电能质量问题。
但是，现有技术中风电机组的闪变计算方法为根据风速进行区间统计后计算得出相应的闪变系数，但是此方法在测试期间需要额外增加测风的设备；计算方面，由于风速只作为风电机组的参考量，并不参与控制，按照风速区间统计方法进行统计计算，不能真实的反映出风电机组在不同运行状态及运行环境下的闪变值。
因此，需要提供一种改进的风电机组闪变确定方法。
发明内容
为克服上述现有技术的不足，本发明提供一种基于功率区间的风电机组闪变确定方法。
实现上述目的所采用的解决方案为：
一种基于功率区间的风电机组闪变确定方法，其改进之处在于：所述方法包括以下步骤：
I、确定测试点；
II、根据风电机组的运行状态在不同的功率区间内进行数据采集；
III、根据风电机组年平均输出功率和电网阻抗角计算闪变值。
进一步的，所述步骤I中，选取风电机组出口变压器低电压侧为所述测量点。
进一步的，所述步骤II中，所述功率区间根据输出的有功功率或额定功率的值进行分区；
所述功率区间包括(-0.05～0.1]、(0.1～0.2]、(0.2～0.3]、(0.3～0.4]、(0.4～0.5]、(0.5～0.6]、(0.6～0.7]、(0.7～0.8]、(0.8～0.9]、(0.9～1.05]。
进一步的，所述步骤II包括以下步骤：
S201、从三相瞬时电压和三相瞬时电流中确定各所述功率区间内相电压和相电流，确定有功功率和无功功率；
S202、根据有功功率确定功率区间；
S203、获取不同所述功率区间10min时间序列的所述测量点的三相瞬时电压和三相瞬时电流。
进一步的，所述步骤III中，计算闪变包括以下步骤：
S301、构建虚拟电网，根据所述虚拟电网确定模拟电压的瞬时值ufic(t)；
S302、获取各功率区间内10min时间序列的有功功率，确定各有功功率对应的闪变扰动系数Pst；
S303、按下式确定各功率区间的闪变系数：
其中，Pst为闪变扰动系数，Sk,fic/Sn为短路比；
S304、按下式确定功率区间的加权系数：
Pua为年平均输出功率，Pi为第i个功率区间内功率，Ni为第i个功率区间内测量闪变系数的实际出现数目，N为闪变系数在各功率区间的总个数；
S305、闪变系数加权累计分布概率统计。
进一步的，所述步骤S306中，闪变系数加权累计分布概率统计：
S3051、将所述闪变值降序排序，确定最大值对应的分布概率为100％，之后的闪变系数对应的加权累计分布顺序递减，递减值根据各功率区间对应的加权系数确定，确定各功率区间加权系数递减确定值为：
其中，wi为功率区间的加权系数，Ni为第i个功率区间内测量闪变系数的实际出现数目，Nb为第i个功率区间的闪变系数的递减个数；
S3052、根据各闪变系数对应的测量值的加权累计分布，确定最接近95％分布概率的闪变系数值为不同年平均输出功率及电网阻抗角下的闪变值c(ψk,Pua)。
进一步的，所述步骤III中，通过获取不同阻抗角下的闪变扰动系数，确定不同的闪变系数；
所述阻抗角分别取30°、50°、70°和85°。
进一步的，所述步骤III中，在不同的所述年平均输出功率下，确定不同的闪变系数；
所述年平均输出功率分别取0.2pu、0.5pu、0.7pu和1.0pu。
与现有技术相比，本发明具有以下优异效果：
1、本发明提供的方法采用功率分区间的方法进行风电机组的闪变计算，不需单独采集风速，有效的降低了试验成本。
2、本发明提供的方法根据实时获取的数据进行计算，减少部分计算过程，提高数据的精度，从而减少了不确定度的来源，提高计算的准确度；优化了现有技术的方法不能真实的反映出风电机组在不同运行状态及运行环境下的闪变值的缺陷。
附图说明
图1为风电机组闪变测试采集点示意图；
图2为本实施例中风电机组闪变计算流程图；
图3为本实施例中闪变仪计算流程图；
图4为本实施例中电网阻抗角k分别为30°，50°，70°和85°时Pst值随功率的变化图。
具体实施方式
下面结合附图对本发明的具体实施方式做进一步的详细说明。
本发明提供一种基于功率区间的风电机组闪变确定方法，该方法包括以下步骤：
步骤一、确定测试点；
步骤二、根据风电机组的运行状态在不同的功率区间内进行数据采集；
步骤三、根据风电机组年平均输出功率和电网阻抗角计算闪变值。
步骤一、确定测量点。
为准确获得风电机组在不同运行工况下的闪变值，选用风电机组出口变压器低压侧为所述测量点。
如图1所示，图1为风电机组闪变测试采集点示意图；获得风电机组出口变压器低压侧三相电流ILV和风电机组出口变压器低压侧三相电压ULV。
步骤二、根据风电机组的运行状态在不同的功率区间内进行数据采集。
数据采集需要确定风电机组的工作状态。
测试条件为：被测风电机组应已并网运行一段时间；测试期间，风电机组应可在不同工况下连续运行；为更真实的反映风电机组的闪变，测量点应选择在风电机组出口变压器低压侧。
在风电机组连续运行状态下测量闪变系数时，风电机组的无功功率输出应尽可能为零。
数据采集包括以下步骤：
获取测量点的三相瞬时电压和三相瞬时电流，确定基波正序分量的有功功率和无功功率。根据无功功率判断数据是否有效，并根据有功功率判断采集数据所处区间。具体方法如下：
①、测量获得相电压和相电流，确定单基波周期T内基波分量的傅里叶系数。A相相电压Ua的计算等式如下式：
其中，f1为基波频率。
②、按下式确定所述基波正序分量的电压矢量分量及电流矢量分量：
③、进一步的，按下式确定基波正序分量的有功功率P1+和无功功率Q1+：
以上以A相为例，B、C相相似。
④、在不同的功率区间内进行数据采集。
首先，确定功率区间。本发明中，对有功功率进行分区。分区如下表1所示，表1为有功功率区间表。
表1有功功率区间
然后，风电机组连续运行时，在风电机组出口变压器低压侧测量三相瞬时电流和三相瞬时电压；测量数据时应满足以下要求：
(1)、在每个功率区间内，至少采集5个10min时间序列的三相瞬时电压和电流测量数据。此处功率为10min平均值。
(2)、处理测量数据时应剔除风电机组连续运行状态以外的采集数据。
三、计算风电机组闪变。如图2所示，图2为本实施例中风电机组闪变计算流程图；具体包括以下步骤：
1、构建虚拟电网，根据所述虚拟电网确定模拟电压的瞬时值ufic(t)。方法如下：
虚拟电网由一个瞬时值为u0(t)的理想相电压源和电阻Rfic及电感Lfic串联组成的电网阻抗表示。风电机组用电流发生器im(t)表示，im(t)为机组线电流的测量瞬时值。该模型中模拟电压的瞬时值ufic(t)为：
其中理想电压源u0(t)可以通过各种不同的方式产生。但必须满足理想电压的两个特性：
①、理想电压不应有任何波动，即电压闪变为零；
②、u0(t)应与测量电压的基波具有相同的电气角αm(t)。
电阻Rfic和电感Lfic根据不同的电网阻抗角ψk计算得出，计算不同阻抗角下的瞬时值，从而确定不同电网阻抗角下不同的闪变扰动系数。
上述电网阻抗角ψk分别取30°、50°、70°和85°。
2、确定闪变扰动系数。
运用闪变仪(该闪变仪符合IEC 61000-4-15标准要求)计算每个功率区间段内的各10min数据的闪变扰动系数Pst。如图3所示，图3为本实施例中闪变仪计算流程图。
3、确定闪变系数。包括以下步骤
①、根据步骤2中，获得虚拟电网中每10min时间序列的闪变扰动系数Pst。
②、确定短路比。本实施例中，采样IEC 61400-21:2008标准推荐采用20～50之间的短路比Sk,fic/Sn。
③、按下式确定各闪变扰动系数值对应的闪变系数。
其中，Pst为闪变扰动系数，Sk,fic/Sn为短路比。
4、确定功率区间的加权系数。
采用功率分区间方法进行闪变计算，在不考虑机组运行故障或限功率运行的情况下，风电机组输出功率与风电机组平均风速相关，即风速越大，输出功率越大。
按下式确定风电机组输出功率由<0.1pu至>0.9pu之间每个0.1pu功率区间的加权系数wi：
参考风电机组风速年平均分布，风电机组输出功率假定也满足瑞利分布，根据假定年平均输出功率为Pua(所述年平均输出功率取0.2pu、0.5pu、0.7pu和1.0pu)，Pi为第i个功率区间内功率，所述功率取10min时间序列的三相瞬时电压和电流的计算出各有功功率的平均值；Ni为第i个功率区间内测量闪变系数的实际出现数目，N为闪变系数在各功率区间的总个数。
5、闪变系数加权累计分布概率统计
所有每10min采集数据计算得到的闪变值，按照由大到小的顺序依次排列，最大值对应的分布概率为100％，之后的闪变系数对应的加权累计分布顺序递减，递减值根据各功率区间对应的加权系数确定，各功率区间加权系数递减确定值为：
wi为第i个功率区间的加权系数，Ni为第i个功率区间内测量闪变系数的实际出现数目，Nb表示第i个功率区间内闪变系数的递减个数。
根据各闪变系数对应的测量值的加权累计分布，确定最接近95％分布概率的闪变系数值，认为该值为不同年平均输出功率及电网阻抗角下的闪变系数c(ψk,Pua)。
确定的多个闪变系数，分别是不同年平均输出功率及电网阻抗角下的闪变系数。
按照风电机组的年平均输出功率Pua分别为0.2pu、0.5pu、0.7pu和1.0pu时，电网阻抗角ψk分别为30°、50°、70°和85°的条件下，测试风电机组的闪变系数c(ψk,Pua)，取结果取累积分布概率为95％时各条件下对应的闪变系数值。如下表2所示，表2为一具体实施例用功率分区间方法进行闪变计算出的结果。
表2基于功率区间方法计算的闪变系数
本实施例中，运用上述方法测试闪变，采集风电机组出口变压器低压侧三相电压、三相电流，通过采集数据得到在不同电网阻抗角下Pst值随功率的变化图，如附图4所示。
最后应当说明的是:以上实施例仅用于说明本申请的技术方案而非对其保护范围的限制,尽管参照上述实施例对本申请进行了详细的说明,所属领域的普通技术人员应当理解:本领域技术人员阅读本申请后依然可对申请的具体实施方式进行种种变更、修改或者等同替换，但这些变更、修改或者等同替换，均在申请待批的权利要求保护范围之内。
Claims (6)
1.一种基于功率区间的风电机组闪变确定方法，其特征在于：所述方法包括以下步骤：
I、确定测试点；
II、根据风电机组的运行状态在不同的功率区间内进行数据采集；
III、根据风电机组年平均输出功率和电网阻抗角计算闪变值；
所述步骤III中，计算闪变包括以下步骤：
S301、构建虚拟电网，根据所述虚拟电网确定模拟电压的瞬时值ufic(t)；
S302、获取各功率区间内10min时间序列的有功功率，确定各有功功率对应的闪变扰动系数Pst；
S303、按下式确定各功率区间的闪变系数：
其中，Pst为闪变扰动系数，Sk,fic/Sn为短路比；
S304、按下式确定功率区间的加权系数：
Pua为年平均输出功率，Pi为第i个功率区间内功率，Ni为第i个功率区间内测量闪变系数的实际出现数目，N为闪变系数在各功率区间的总个数；
S305、闪变系数加权累计分布概率统计；
所述步骤S305中，闪变系数加权累计分布概率统计包括以下步骤：
S3051、将所述闪变值降序排序，确定最大值对应的分布概率为100％，之后的闪变系数对应的加权累计分布顺序递减，递减值根据各功率区间对应的加权系数确定，确定各功率区间加权系数递减确定值为：
其中，wi为功率区间的加权系数，Ni为第i个功率区间内测量闪变系数的实际出现数目，Nb为第i个功率区间的闪变系数的递减个数；
S3052、根据各闪变系数对应的测量值的加权累计分布，确定最接近95％分布概率的闪变系数值为不同年平均输出功率及电网阻抗角下的闪变值c(ψk,Pua)；
其中，ψk为电网阻抗角。
2.如权利要求1所述的方法，其特征在于：所述步骤I中，选取风电机组出口变压器低电压侧为所述测试点。
3.如权利要求1所述的方法，其特征在于：所述步骤II中，所述功率区间根据输出的有功功率或额定功率的值进行分区；
所述功率区间包括(-0.05～0.1]、(0.1～0.2]、(0.2～0.3]、(0.3～0.4]、(0.4～0.5]、(0.5～0.6]、(0.6～0.7]、(0.7～0.8]、(0.8～0.9]、(0.9～1.05]。
4.如权利要求1所述的方法，其特征在于：所述步骤II包括以下步骤：
S201、从三相瞬时电压和三相瞬时电流中确定各所述功率区间内相电压和相电流，确定有功功率和无功功率；
S202、根据有功功率确定采集数据所在的功率区间；
S203、获取不同所述功率区间10min时间序列的所述测试点的三相瞬时电压和三相瞬时电流。
5.如权利要求1所述的方法，其特征在于：所述步骤III中，通过获取不同阻抗角下的闪变扰动系数，确定不同的闪变系数；
所述阻抗角分别取30°、50°、70°和85°。
6.如权利要求1所述的方法，其特征在于：所述步骤III中，在不同的所述年平均输出功率下，确定不同的闪变系数；
所述年平均输出功率分别取0.2pu、0.5pu、0.7pu和1.0pu。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410637797.6A CN105629168B (zh) | 2014-11-06 | 2014-11-06 | 一种基于功率区间的风电机组闪变确定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410637797.6A CN105629168B (zh) | 2014-11-06 | 2014-11-06 | 一种基于功率区间的风电机组闪变确定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105629168A CN105629168A (zh) | 2016-06-01 |
| CN105629168B true CN105629168B (zh) | 2018-11-13 |
Family
ID=56044286
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410637797.6A Active CN105629168B (zh) | 2014-11-06 | 2014-11-06 | 一种基于功率区间的风电机组闪变确定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105629168B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106339567B (zh) * | 2016-11-01 | 2023-05-26 | 中国电力科学研究院 | 一种基于系统阻抗的风电场闪变传递计算方法及其装置 |
| CN110161299A (zh) * | 2019-04-23 | 2019-08-23 | 中国电力科学研究院有限公司 | 一种风电场短时闪变值确定方法及装置 |
| CN110346656B (zh) * | 2019-05-21 | 2023-08-11 | 中国电力科学研究院有限公司 | 一种风电机组三相电压不平衡度确定方法及装置 |
| CN116545023B (zh) * | 2023-07-06 | 2023-11-14 | 中国电力科学研究院有限公司 | 一种风电机组并网点闪变特性的仿真验证方法及装置 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102749555A (zh) * | 2012-06-20 | 2012-10-24 | 广东电网公司电力科学研究院 | 一种供电系统中闪变源的定位方法 |
| CN103323738A (zh) * | 2012-03-20 | 2013-09-25 | 华北电力科学研究院有限责任公司 | 一种多馈入系统闪变源定位方法 |
| CN103926456A (zh) * | 2014-04-10 | 2014-07-16 | 宁波恒力达科技有限公司 | 基于改进fft的闪变值计算方法及逆变器 |
-
2014
- 2014-11-06 CN CN201410637797.6A patent/CN105629168B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103323738A (zh) * | 2012-03-20 | 2013-09-25 | 华北电力科学研究院有限责任公司 | 一种多馈入系统闪变源定位方法 |
| CN102749555A (zh) * | 2012-06-20 | 2012-10-24 | 广东电网公司电力科学研究院 | 一种供电系统中闪变源的定位方法 |
| CN103926456A (zh) * | 2014-04-10 | 2014-07-16 | 宁波恒力达科技有限公司 | 基于改进fft的闪变值计算方法及逆变器 |
Non-Patent Citations (1)
| Title |
|---|
| 风电场电能质量测试方法与实践;陈兵;《百度文库》;20140514;第17页闪变测试部分 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105629168A (zh) | 2016-06-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104268332B (zh) | 一种基于区域互联电网频率特性的高频切机仿真配置方法 | |
| Liu et al. | A dynamic equivalent model for DFIG-based wind farms | |
| CN104362678B (zh) | 一种风电场接入对局部电网影响应对措施的分析方法 | |
| CN105629168B (zh) | 一种基于功率区间的风电机组闪变确定方法 | |
| CN110346656A (zh) | 一种风电机组三相电压不平衡度确定方法及装置 | |
| CN104617578B (zh) | 一种含风电场电力系统的可用输电能力的获取方法 | |
| CN104502775A (zh) | 一种电能质量对变压器综合能耗影响的定量分析方法 | |
| Elizondo et al. | Model reduction, validation, and calibration of wind power plants for dynamic studies | |
| CN106406272A (zh) | 一种风电场中静止无功发生器的控制器性能测试方法 | |
| CN106610446A (zh) | 一种基于功率区间自适应划分的风电机组谐波电流确定方法 | |
| CN109412163A (zh) | 一种分散式风电接入配电网的精确潮流计算方法 | |
| Xie et al. | Online parameter determination based adaptive single‐phase reclosing scheme for wind‐powered outgoing lines with shunt reactors | |
| Karegar et al. | Islanding detection of wind farms by THD | |
| CN104850751B (zh) | 一种电流质量评估方法 | |
| CN104577983B (zh) | 风电场主变压器纵联差动保护适应性分析方法 | |
| CN102928644B (zh) | 风力发电机组并网点电压闪变检测方法及装置 | |
| Wu et al. | Research on improvement of line loss algorithm based on three-phase unbalance degree | |
| CN102957166B (zh) | 一种基于轨迹灵敏度的风电配置比例快速计算方法 | |
| Zhan et al. | Oscillatory stability analysis for wind power systems based on operating-condition dependent impedance model | |
| CN112217220B (zh) | 新能源送出多端柔性直流电网安全稳定防御方法及装置 | |
| Li et al. | The optimization of undervoltage load shedding strategy | |
| Khoa et al. | Voltage Estimation of Distribution Network Integrated with Distributed Generation Using Remote Measurement Data | |
| CN105355122A (zh) | 一种虚负荷模拟组网发生装置及发生方法 | |
| CN105207249A (zh) | 一种风电场恒定功率因数控制系统及其控制方法 | |
| CN106610471B (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 |