CN102012956B - 一种考虑风电场输入风速和风向随机波动的风电场等值方法 - Google Patents
一种考虑风电场输入风速和风向随机波动的风电场等值方法 Download PDFInfo
- Publication number
- CN102012956B CN102012956B CN201010576662.5A CN201010576662A CN102012956B CN 102012956 B CN102012956 B CN 102012956B CN 201010576662 A CN201010576662 A CN 201010576662A CN 102012956 B CN102012956 B CN 102012956B
- Authority
- CN
- China
- Prior art keywords
- wind
- equivalent
- energy turbine
- turbine set
- wind turbines
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 230000000694 effects Effects 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims description 28
- 238000006243 chemical reaction Methods 0.000 claims description 16
- 230000003068 static effect Effects 0.000 claims description 7
- 238000006116 polymerization reaction Methods 0.000 claims description 6
- 230000001052 transient effect Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000013016 damping Methods 0.000 claims description 4
- 230000000717 retained effect Effects 0.000 claims description 4
- 238000002347 injection Methods 0.000 claims description 2
- 239000007924 injection Substances 0.000 claims description 2
- 230000005611 electricity Effects 0.000 abstract description 6
- 238000004458 analytical method Methods 0.000 abstract description 5
- 238000004519 manufacturing process Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000005684 electric field Effects 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000005624 perturbation theories Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- 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
- Control Of Eletrric Generators (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了一种考虑风电场输入风速、风向随机波动、风电机组间尾流效应及风电场内电网络等因素建立风电场随机模型的方法。技术方案:考虑风电场风速和风向随机波动对风电机组分组;根据等值前后风电场与电网之间功率交换不变的原则对风电场功率等值;根据等值前后风电场动态特性不变的原则,对风电机组的发电机、轴系和控制***等模型参数等值;通过合理地选择等效风电机组基准值,使得等值风电机组模型的详细程度与单台风电机组的相同。利用此方法建立的模型进行风电场并网分析,克服了不考虑风电场风速、风向随机波动、风电机组间尾流影响以及电网络等因素带来的缺点;为电网调度合理安排生产计划提供了较为准确的风电数据。
Description
技术领域
本发明涉及一种并网风电场的等值方法,特别涉及一种采用恒速风电机组的风电场,考虑风电场输入风速和风向随机波动的风电场等值方法,用于大型风电场并网分析。
背景技术
随着我国近年来风力发电的快速发展,进行风电场接入电力***的研究和评估、开展***及关键节点接纳风电能力的研究是当前风电发展中迫切需要解决的重要问题。从电力***的角度来看,对风电场进行研究所关心的不是风电场内部每台风电机组的特性而是风电场作为一个整体的动态特性以及对电力***的影响,在风电场接入电力***的分析中不可能也没有必要把风电场内每台风电机组都作为一个单独元件列入仿真程序中进行分析,随着风电场规模的增大,这一特点越来越明显。但是,风电场不同于常规的发电厂,风电场是由大量分散布置的风电机组群组成的,风电场内每台风电机组的输入风速除了都随着风电场来风风速和来风风向的波动而波动外,风电机组间的尾流效应还可能使得不同安装地点的风电机组输入风速具有明显差异,造成同一时刻风电场内风电机组的运行状态不完全相同。因此,在风电场并网研究中需对风电场进行合理简化以建立满足分析要求的动态模型。
外部电网遭受故障和风速波动两种情况下风电场的动态特性及其对电力***的影响是风电场并网分析的两个重要方面。在这两种情况下对风电场输出特性研究的侧重点不同,因此风电场的等值方法也就不同。本发明主要研究用于风速波动风电场动态特性研究的风电场动态模型的等值方法。
建立风电场的等值模型就是考虑影响风电场特性的因素合理简化风电场,包括对风电机组合理分组、确定每组风电机组等值输入风速、等值风电机组静态和动态模型建立及参数确定以及风电场内部电网简化等几个方面。
目前,现有的风电机组分组方法有:
(1)把风电场内所有风电机组归为一组。
(2)连结于同一条集电线路上风电机组归为一组。
(3)对于风电机组排列布置规则的风电场,按照风电机组安装位置进行分组。如海上风电场,假设与风电场来风风向垂直的每排风电机组输入风速相同,将每排风电机组归为一组并等值成一台等值风电机组。
(4)当风电机组排列布置不规则时,根据发电机的机械暂态数学模型的特征根对双馈风电机组分组。
现有风电场等值方法有:
(1)按奇异摄动理论进行风电场模型降阶,将描述风电机组模型的状态变量分解为快变化和慢变化两部分,通过忽略快变化的状态变量,得到用慢变化状态变量表示的降阶风电场模型。
(2)对输入风速不同的风电机组等值时,风电机组的发电机和风力机等值方法不同。按照容量加权法将每台风电机组的发电机等值成一台等值发电机,等值发电机的容量为所有发电机容量之和,等值发电机模型的详细程度与单台发电机模型相同;但是在对风电机组的风力机部分进行等值时,不是把输入风速不同的风力机等值成一台等值风力机,而是采用输入风速不同的多台简化风力机表示。
(3)按照容量加权法将归为一组的风电机组等值成与单台风电机组模型详细程度相同的等值风电机组,其额定容量为每台风电机组额定容量之和。
等值风电机组参数的确定方法有:
(1)采用遗传算法优化的风电机组参数。利用遗传算法的全局寻优功能,优选风电场等值模型参数,最终可找到一组用优化参数表示的风电场等值模型。不过,它没有考虑风电场风速随时间的变化和风电场中作用到各个风电机组的风速分布,而是将风电场中的风速分布看作一系列稳态过程来得到整个风电场的风速分布。
(2)采用基于风电场输出特性的风电场动态等值模型参数优化解法。兼顾了风电场机组众多、风能波动较强的特点,将风电场输出特性引入风电等值参数求解过程,将风电场等值为一台风电机组。
(3)按容量加权方法对风电机组参数进行等值计算。利用容量加权法对恒速风电机组和双馈变速风电机组模型参数进行等值计算。
(4)把风电机组的阻抗和风电机组间连接线路的阻抗进行串并联计算,求得阻抗值作为风电场的等值阻抗。
风电场内部线路长度随着风电场装机容量的增加而增加,它将影响风电场的输出特性。现有处理风电场内部线路的方法有:
(1)在某些小型风电场建模中,忽略风电场内部线路,只是将风电机组进行等值。
(2)利用等值损耗模型对连结于同一集电线路上风电机组间的线路进行简化,得出等值风电机组与风电场升压变之间等值线路的阻抗。但是当输入风速相同的风电机组连接于不同集电线路时,利用此方法不能计算等值风电机组与风电场升压变之间等值线路的阻抗。
综上所述,在风电场接入电力***的分析研究中,国内外研究机构建立了用于风电场电能质量、风电场动态特性以及可靠性等研究的风电场动态模型。但是在风电场建模过程中仍存在一些问题,如:
(1)没有考虑风电场输入风速和风向随机变化对风电机组分组的影响。风电场输入风速、风向和风电机组间尾流影响直接影响每台风电机组的输入风速,并且随风速和风向的随机波动而变化,风电场内风电机组的分组也应是随机变化的,不是固定分组;
(2)在建立风电场模型时,没有考虑风电场内部电网络、风电机组间尾流效应等影响风电场的输出特性因素。风电机组类型、内部电网络、风电机组间尾流效应以及风电场的地势地貌等因素都同时影响风电场的输出特性,在对风电场建模时应同时考虑这些因素。
结合目前风电场简化模型研究现状以及我国风力发电的发展来看,进一步深入研究风速和风向的随机波动、风电机组间尾流效应和风电场内部线路对风电场输出特性的影响,并建立风电场等值模型,把风电场建模方法应用到实际风电场的并网研究中是非常必要的。
发明内容
本发明的目的是考虑风电场输入风速和风向的随机波动、风电机组间的尾流效应及风电场内部线路等因素对风电场输出特性的影响,提出建立风电场等值模型的方法。包括风电机组分组方法;每组风电机组母线的等值;风电场内部电网络的化简;风电场功率的等值;风电机组的发电机、轴系和控制***等模型及其参数等值。利用此方法可建立风电场的随机模型,研究风电场输入风速和风向随机波动时风电场输出特性,为电力***合理调度提供风电场运行数据。
本发明为实现上述目的,采用的技术方案是:
1.考虑风电场输入风速和风向随机波动对风电机组分组并计算每组风电机组的等值风速
对于实际运行的风电场来说,场内每台风电机组输入风速、风电机组分组都随风电场输入风速和风向的随机波动而变化,因此,要想对正常运行的风电场进行仿真就需要建立随风速和风向随机波动的风电场随机模型。针对这一问题,考虑风电场输入风速(变化范围为4m/s~25m/s)和风向(变化范围为0°~360°)变化范围比较宽,并且是随机波动的特点,提出了利用相关系数对风电机组分组的方法,并可以在已知风电场输入风速的情况下方便地求出每组风电机组的等值输入风速;利用此方法就可以建立风电场的随机模型。
风电机组相关系数的计算框图如图1所示。通过对风电场每个风速和风向上风电机组相关系数的计算可以得到如图2所示的3维相关系数矩阵。3维系数矩阵的行表示风向,其行数依赖于所选风向的间隔;矩阵的列对应于风电机组的台数;系数矩阵的第三个轴对应于风电场的输入风速,其列数依赖于所选风速的间隔;系数矩阵内每个元素CI是风电机组的相关系数。
对于给定的风电场输入风速和风向,查询图2得出对应的相关系数,然后把相关系数相等的风电机组归为一组,建立风电场模型。当研究风速和风向随机波动风电场输出特性时,可通过识别风电场输入风速和风向来调用风电场相应等值模型,然后对输入风速和风向随机波动的风电场进行仿真,程序框图如图3所示。
在已知风电场输入风速和风向的情况下,利用风电机组相关系数可以计算出每台风电机组的等值输入风速ve_i:
式中:ve_i为第i组风电机组等值风速;vinf为风电场输入风速;vstep为风速步长;CIi为第i组风电机组的相关系数;sgn()为符号函数;NOT()为非函数。
风电机组分组完成后,下面就需要对风电场进行等值。风电场等值包括静态等值和动态等值两部分。
2.风电场的静态等值
风电场的静态等值包括风电机组出口母线等值、网络化简和风电场功率等值三部分。
2.1风电机组母线等值
把风电场内相关系数相等的风电机组母线集合为{C},与{C}相联的母线的集合为{R},相关系数相等的风电机组母线聚合就是用一个等值母线e代替这一组母线{C},如图4所示,{R}表示与{C}相联的母线集合。
相关系数相等的风电机组母线的等值是将被合并母线{C}通过复变比的理想移相变压器连接到等值母线e,然后将移相变压器化为非移相变压器,再通过网络化简的方法消去母线{C},如图5所示。
等值母线的电压幅值和相角δe等于被合并母线电压幅值和相角的平均值:
其中:m为被等值的风电机组台数,为第k台风电机组出口电压,δk为第k台风电机组出口电压的相角。
2.2网络化简
网络化简主要是通过对***的节点导纳矩阵进行消去运算,消去不需要保留的节点。网络化简的关键是使等值后稳态潮流保持等值前的状态,而动态误差尽可能小。
按照网络化简要求,可将原风电场网络划分为三个部分:
(1)内部***(即风电场内每组风电机组的等值母线e1,e2…eN)
该部分是需要保留的部分,其节点集合用{I}表示。
(2)边界节点(即风电场的等值风电机组的连接母线P)
该部分也是需要保留的部分,其节点集合用{B}表示。
(3)外部***(即风电场内部除e1,e2…eN和P节点以外的所有节点)
该部分是需要消去的部分,其节点集合用{D}表示。
将原型***节点导纳方程按照I,B,D集合用分块矩阵形式表示如下:
消去外部***节点电压和注入电流,则化简后***方程为
其中:
简化后***节点导纳方程只含内部***和边界节点的电流和电压。
2.3风电场功率的等值
稳态时风电场视在功率Se等于所有风电机组视在功率Sij之和:
式中,j为风电机组的组数,i为风电机组的台数。
风电场有功功率Pe等于风电场内每台风电机组Pij向电网提供有功功率之和:
风电场内恒速风电机组的接线如图6所示,风电场等值风电机组从电网中吸收的无功功率Qe为:
其中,Qi-1,j为风电场内每台风电机组从电网中吸收的无功功率;ΔQleq为风电场内部网络所消耗的无功功率,i,j分别决定风电机组在风电场内位置,N为风电场内风电机组的组数,n为每组风电机组中风电机组的台数。
风电场等值风电机组在未加入无功补偿时吸收的无功功率Qg,e为:
其中,Xg,e、Rg,e为风电场等值风电机组的等值电抗和电阻,Ug,e为等值风电机组出口电压。
为保证等值前后风电场与电网之间的无功功率交换保持不变,则等值风电机组补偿的无功功率ΔQg为:
ΔQg=Qe-Qg,e(11)
等值电容器的阻抗XCg为:
3.风电场的动态等值
传统的发电机聚合方法是在频域中拟合发电机各个环节的传递函数,通过迭代优化的方法使等值机和同调机群具有相近的频域特性进行聚合。本文采用保留发电机方程系数矩阵结构不变的时域聚合法。该方法是在时域中进行,直接计算等值机参数,无需迭代,计算速度快。
对风电机组动态等值时假定如下:
(1)被等值的风电机组有相同的输入风速;
(2)被等值的风电机组通过理想变压器移到等值母线e上;
(3)等值风电机组机械功率和电磁功率与被等值风电机组总机械功率和总电磁功率相等。
恒速风电机组中感应发电机采用三阶实用模型,用发电机暂态电抗表示的第i台风电机组的电流方程为:
上式可简化为:
ii=Ai(ui-E′i)(14)
3.1恒速风电机组等值发电机暂态阻抗的计算
等值发电机也有与式(14)类似的电流方程式:
ie=Ae(ue-E′e)(15)
式中,ie为流入等值风电机组的电流;ue为等值风电机组的出口电压;E′e为等值风电机组的内电势;Ae为等值风电机组电流方程中的矩阵,
根据假设(3),等值发电机电流等于被等值的发电机电流之和;风电机组输入风速(相关系数)相等(或相似)的风电机组出口电压和内电势相同(或相似);等值母线e的电压等于被合并母线电压的平均值,即等值母线电压ue=ui,等值风电机组内电势E′e=E′i。由此可得:
经过求逆变换得出等值感应发电机的等值阻抗为:
3.2等值感应发电机暂态时间常数的计算
等值感应发电机的暂态时间常数是通过对感应发电机状态变量的等值来求取。由矩阵表示的第i台感应发电机转子微分方程为:
式中,上式也可写为:
相应地,等值发电机也有与式(19)类似的转子微分方程:
式中:ie为流入等值风电机组的电流矩阵;E′e为等值感应发电机内电势矩阵;为等值感应发电机内电势导数矩阵;Ce、De为等值感应发电机转子微分方程中的矩阵。由于且流入m台风电机组的等值机组电流等于流入每台风电机组电流之和,所以式(20)变为:
比较(20)和(21)两式可得:
Ci=Ce,Di=mDe(22)
则等值感应发电机的暂态时间常数T′d0e为:
T′d0e=T′d0(23)
3.3等值风电机组机械输入功率的等值
等值风电机组的等值机械功率PTeq为:
其中:PTi为第i台风力机从风中捕获的风能;ρ为空气密度;β为叶片桨距角;λ称为叶尖速比;Cpi为风力机的功率系数:vi为风电机组的等值输入风速。
3.4风电机组机械参数的等值
风电机组的机械运动方程为:
式中,Jt、Jg分别为风力机和发电机轴的转动惯量,ωt、ωg分别为风力机和发电机转速,Tt为作用在风力机上的转矩,Te为发电机输出的电磁转矩,Tm为风力机轴的输出转矩,K为轴的刚度系数,D为轴的阻尼系数,θtg为扭转角。
风电机组等值准则是等值机的机械功率和电磁功率分别等于被等值风电机组机械功率和电磁功率之和。由于风电机组的输入风速相等,则发电机具有相同的转速,将被等值风电机组的转子运动方程累加得等值机的转子运动方程为:
即可得出:
若在以等值风电机组的额定电压和等值风电机组的额定功率作为基准值的标幺值***中,等值风电机组转动惯量、刚度系数和阻尼系数的标幺值与单台风电机组的相等,即:
Jte=Jt;Jge=Jg;De=D;Ke=K(28)
3.5风电机组控制***的等值
在以等值风电机组额定容量和额定电压作为基准值的标么制***下,等值风电机组控制模型的参数也与单台风电机组的参数相同,只是把控制***中的相关功率变为等值风电机组的功率。同样保护***模型也与单台风电机组的相同。
本发明的有益效果体现在:利用本发明提出的风电场方法,克服了以往不考虑风电场来风风速、来风风向变化、风电机组间尾流影响和风电场内部线路等影响而建立风电场模型的缺点;利用此方法可方便建立风电场输入风速和风向变化时建立风电场的随机模型,进一步决定风电场的出力,研究并网风电场的输出特性及对电网运行特性的影响,同时为电网调度部门合理安排生产计划、减少***备用容量提供了较为准确的风电数据。
附图表说明
图1风电机组相关系数的计算框图;
图2每个风速和风向上风电机组3维相关系数矩阵;
图3风电场输入风速、风向随机波动时风电场仿真框图;
图4相关系数相等的风电机组母线聚合;
图5等值网络图;
图6风电场内恒速风电机组连接
图7风电场内风电机组的排列布置;
图8风向γ=45°时风电场等值模型;
图9风电场来风风速;
图10风向为γ=45°时风电场内每台风电机组的输入风速;
图11风电场的有功功率和出口电压的变化情况
具体实施方案
下面利用附图和实施例子对本发明进一步描述。本发明提出的考虑风电场来风风速、来风风向变化、风电机组间尾流影响和风电场内部线路等建立风电场随机模型的方法,解决传统的把风电场内所有风电机组归为一组或根据风电机组排列分组进行等值无法处理的问题,提高了风电场动态特性研究的准确性。具体实施方案如下:
1.风电场由16台容量为1.5MW的恒速风电机组组成如图7所示,叶轮直径为70m。风电机组分4排,每排风电机组中相邻两台风电机组之间的距离为叶轮直径的7倍,相邻两排风电机组间的距离也为叶轮直径的7倍。每排风电机组通过距离为1km、型号为LGJ-185的架空线路连接到风电场主变,每排中相邻两台风电机组间线路采用型号为LGJ-185架空线。
2.利用风电机组相关系数计算方法,可算出风电机组的相关系数,并把它存储为三维矩阵。通过计算比较可得,在一个风速间隔内可以用这个间隔的平均风速所计算出的相关系数作为3维系数矩阵中系数。
3.任取其中六个风向计算风电场内每台风电机组的输入风速,并把相关系数相等的风电机组归为一组,风电机组的分组情况如表1所示。
4.选取任一风向建立风电场模型,并验证建模方法的合理性。风向γ=45°时风电场等值电路如图8所示。每个等值风电机组的参数如表2所示。
5.风电场来风风速如图9所示,通过风电机组输入风速计算方法,得出风电场来风风向为γ=45°时风电场内每台风电机组的输入风速如图10所示。分别采用风电场详细模型和等值模型进行仿真,可得风电场的有功功率和出口电压的变化情况如图11所示。从图11中可以看出两种模型下风电场的动态特性非常相近,也就是说风电场内风电机组的参数等值和风电场内部电网的等值是较为可行的。
表1风电机组的分组情况
表2等值风电机组参数
Rse(p.u.) | X′e(p.u.) | Ke(p.u.) | Jte(p.u.) | Jge(p.u.) | De(pu) | |
WT1e | 0.001 | 0.156 | 2.1 | 56.4 | 10.5 | 0.105 |
WT2e | 0.0014 | 0.218 | 1.5 | 42 | 9 | 0.09 |
WT3e | 0.0023 | 0.364 | 0.9 | 24.6 | 4.5 | 0.045 |
WT4e | 0.007 | 1.092 | 0.3 | 8.2 | 1.5 | 0.015 |
Claims (6)
1.一种考虑风速和风向随机波动的风电场等值方法,其特征在于:通过风电场输入风速和风向的随机波动、风电机组间尾流效应、风电机组排列布置及风电场内部馈线,建立风电场等值模型;
通过对风电场每个风速和风向上风电机组相关系数的计算可以得到3维相关系数矩阵;3维系数矩阵的行表示风向,其行数依赖于所选风向的间隔;矩阵的列对应于风电机组的台数;系数矩阵的第三个轴对应于风电场的输入风速,其列数依赖于所选风速的间隔;系数矩阵内每个元素CI是风电机组的相关系数;
对于给定的风电场输入风速和风向,查询3维相关系数矩阵的相关系数,然后把相关系数相等的风电机组归为一组,建立风电场模型;当风速和风向随机波动风电场输出特性时,可通过识别风电场输入风速和风向来调用风电场相应等值模型,然后对输入风速和风向随机波动的风电场进行仿真;
在已知风电场输入风速和风向的情况下,利用风电机组相关系数可以计算出每台风电机组的等值输入风速ve_i:
其中,ve_i为第i组风电机组等值风速;vinf为风电场输入风速;vstep为风速步长;CIi为第i组风电机组的相关系数;sgn()为符号函数;NOT()为非函数;
风电机组分组完成后,需要对风电场进行等值;风电场等值包括静态等值和动态等值两部分;
风电场的静态等值步骤:
风电场的静态等值包括风电机组出口母线等值、网络化简和风电场功率等值三部分;
风电机组母线等值步骤:
把风电场内相关系数相等的风电机组母线集合为{C},与{C}相联的母线的集合为{R},相关系数相等的风电机组母线聚合就是用一个等值母线e代替这一组母线{C},;
相关系数相等的风电机组母线的等值是将被合并母线{C}通过复变比的理想移相变压器连接到等值母线e,然后将移相变压器化为非移相变压器,再通过网络化简的方法消去母线{C};
等值母线的电压幅值和相角δe等于被合并母线电压幅值和相角的平均值:
其中,m为被等值的风电机组台数,为第k台风电机组出口电压,δk为第k台风电机组出口电压的相角;
网络化简步骤:
网络化简主要是通过对***的节点导纳矩阵进行消去运算,消去不需要保留的节点;网络化简的关键是使等值后稳态潮流保持等值前的状态,使动态误差尽可能小;
可将原风电场网络划分为三个部分:
内部***,即风电场内每组风电机组的等值母线e1,e2…eN,该部分为保留的部分,其节点集合用{I}表示;
边界节点,即风电场的等值风电机组的连接母线P,该部分为保留的部分,其节点集合用{B}表示;
外部***,即风电场内部除e1,e2…eN和P节点以外的所有节点;该部分为消去的部分,其节点集合用{D}表示;
将原风电场网络节点导纳方程按照I,B,D集合用分块矩阵形式表示如下:
消去外部***节点电压和注入电流,则化简后***方程为:
其中:
简化后***节点导纳方程只含内部***和边界节点的电流和电压;
风电场功率的等值步骤:
稳态时风电场视在功率Se等于所有风电机组视在功率Sij之和:
其中,j为风电机组的组数,i为风电机组的台数;
风电场有功功率Pe等于风电场内每台风电机组Pij向电网提供有功功率之和:
其中,Qi-1,j为风电场内每台风电机组从电网中吸收的无功功率;ΔQleq为风电场内部网络所消耗的无功功率,i,j分别决定风电机组在风电场内位置,N为风电场内风电机组的组数,n为每组风电机组中风电机组的台数;
风电场等值风电机组在未加入无功补偿时吸收的无功功率Qg,e为:
其中,Xg,e、Rg,e为风电场等值风电机组的等值电抗和电阻,Ug,e为等值风电机组出口电压;
为保证等值前后风电场与电网之间的无功功率交换保持不变,则等值风电机组补偿的无功功率ΔQg为:
ΔQg=Qe-Qg,e(11)等值电容器的阻抗XCg为:
风电场的动态等值步骤:
对风电机组动态等值时设定如下:
被等值的风电机组有相同的输入风速;
被等值的风电机组通过理想变压器移到等值母线e上;
等值风电机组机械功率和电磁功率与被等值风电机组总机械功率和总电磁功率相等;
恒速风电机组中感应发电机采用三阶实用模型,用发电机暂态电抗表示的第i台风电机组的电流方程为:
上式可简化为:
ii=Ai(ui-E′i)(14)
恒速风电机组等值发电机暂态阻抗的计算
等值发电机与式(14)类似的电流方程式:
ie=Ae(ue-E′e)(15)
其中,ie为流入等值风电机组的电流;ue为等值风电机组的出口电压;E′e为等值风电机组的内电势;Ae为等值风电机组电流方程中的矩阵,
根据上述设定,等值发电机电流等于被等值的发电机电流之和;风电机组输入风速相等或类似的风电机组出口电压和内电势相同;等值母线e的电压等于被合并母线电压的平均值,即等值母线电压ue=ui,等值风电机组内电势E′e=E′i;由此可得:
等值感应发电机暂态时间常数的计算步骤:
等值感应发电机的暂态时间常数是通过对感应发电机状态变量的等值来求取;矩阵表示的第i台感应发电机转子微分方程为:
式中, 上式也可写为:
相应地,等值发电机有与式(19)类似的转子微分方程:
其中:ie为流入等值风电机组的电流矩阵;E′e为等值感应发电机内电势矩阵;为等值感应发电机内电势导数矩阵;Ce、De为等值感应发电机转子微分方程中的矩阵。由于且流入m台风电机组的等值机组电流等于流入每台风电机组电流之和,所以式(20)变为:
比较式(20)和式(21)两式可得:
Ci=Ce,Di=mDe(22)
则等值感应发电机的暂态时间常数T′d0e为:
T′d0e=T′d0(23)
等值风电机组机械输入功率的等值
等值风电机组的等值机械功率PTeq为:
其中:PTi为第i台风力机从风中捕获的风能;ρ为空气密度;β为叶片桨距角;λ称为叶尖速比;Cpi为风力机的功率系数:vi为风电机组的等值输入风速;
风电机组机械参数的等值步骤:
风电机组的机械运动方程为:
其中,Jt、Jg分别为风力机和发电机轴的转动惯量,ωt、ωg分别为风力机和发电机转速,Tt为作用在风力机上的转矩,Te为发电机输出的电磁转矩,Tm为风力机轴的输出转矩,K为轴的刚度系数,D为轴的阻尼系数,θtg为扭转角;
风电机组等值准则是等值机的机械功率和电磁功率分别等于被等值风电机组机械功率和电磁功率之和;由于风电机组的输入风速相等,则发电机具有相同的转速,将被等值风电机组的转子运动方程累加得等值机的转子运动方程为:
即可得出:
若在以等值风电机组的额定电压和等值风电机组的额定功率作为基准值的标幺值***中,等值风电机组转动惯量、刚度系数和阻尼系数的标幺值与单台风电机组的相等,即:
Jte=Jt;Jge=Jg;De=D;Ke=K(28)。
2.根据权利要求1所述的一种考虑风速和风向随机波动的风电场等值方法,其特征在于:根据风电场输入风速和风向随机波动,把相关系数相等的风电机组归为一组,得出了计算每组风电机组的等效风速的方法。
3.根据权利要求1所述的一种考虑风速和风向随机波动的风电场等值方法,其特征在于:得出了每组等值风电机组母线的聚合方法。
4.根据权利要求1所述的一种考虑风速和风向随机波动的风电场等值方法,其特征在于:根据恒功率等效定则,得出由恒速风电机组组成风电场的有功功率和无功功率的等值方法。
5.根据权利要求1所述的一种考虑风速和风向随机波动的风电场等值方法,其特征在于:得出了恒速风电机组风电场等值机组的暂态电抗和暂态时间常数等值方法。
6.根据权利要求1所述的一种考虑风速和风向随机波动的风电场等值方法,其特征在于:得出了恒速风电机组风电场等值机组输入机械功率的等值方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010576662.5A CN102012956B (zh) | 2010-11-30 | 2010-11-30 | 一种考虑风电场输入风速和风向随机波动的风电场等值方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010576662.5A CN102012956B (zh) | 2010-11-30 | 2010-11-30 | 一种考虑风电场输入风速和风向随机波动的风电场等值方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102012956A CN102012956A (zh) | 2011-04-13 |
CN102012956B true CN102012956B (zh) | 2016-01-20 |
Family
ID=43843128
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201010576662.5A Expired - Fee Related CN102012956B (zh) | 2010-11-30 | 2010-11-30 | 一种考虑风电场输入风速和风向随机波动的风电场等值方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102012956B (zh) |
Families Citing this family (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102254092B (zh) * | 2011-06-16 | 2013-02-27 | 国网电力科学研究院 | 一种含双馈风电机组的大型风电场动态等值方法 |
CN102270256B (zh) * | 2011-06-30 | 2013-04-03 | 内蒙古电力勘测设计院 | 平坦地区风场规划方法 |
CN102235314B (zh) * | 2011-06-30 | 2013-03-06 | 内蒙古电力勘测设计院 | 复杂地形下提高风能利用的方法 |
CN102609590B (zh) * | 2012-02-16 | 2013-08-07 | 中国科学院寒区旱区环境与工程研究所 | 风电场群布局方法 |
CN103093027B (zh) * | 2012-12-06 | 2015-09-16 | 广东电网公司电力科学研究院 | 基于双馈风电场等值模型的电力***分析方法 |
EP2757257B1 (en) * | 2013-01-22 | 2017-11-22 | GE Renewable Technologies | Methods and arrangements for controlling power generators |
CN103277251B (zh) * | 2013-05-24 | 2015-02-04 | 长沙理工大学 | 一种风电场的控制方法及*** |
CN103401236B (zh) * | 2013-07-30 | 2015-12-09 | 华北电力大学 | 一种基于风电场流动相关性的风电场机组分组方法 |
CN103617308B (zh) * | 2013-10-30 | 2016-06-08 | 河海大学 | 一种风电场频域等效模型的构建方法 |
CN103956767B (zh) * | 2014-02-21 | 2016-04-20 | 国家电网公司 | 一种考虑尾流效应的风电场并网稳定性分析方法 |
CN103886178A (zh) * | 2014-02-25 | 2014-06-25 | 国家电网公司 | 一种永磁直驱机组风电场聚合模型参数等值方法 |
CN103853888B (zh) * | 2014-03-10 | 2017-01-25 | 南方电网科学研究院有限责任公司 | 一种适用于无功补偿的海上风电***等值建模方法 |
US9551322B2 (en) | 2014-04-29 | 2017-01-24 | General Electric Company | Systems and methods for optimizing operation of a wind farm |
CN105320825A (zh) * | 2014-06-19 | 2016-02-10 | 国网山西省电力公司电力科学研究院 | 基于简化模型的双馈机组风电场动态等值建模***及方法 |
CN104166777A (zh) * | 2014-08-29 | 2014-11-26 | 重庆大学 | 计及多重相关性的风速矢量数据模拟生成方法 |
CN105631175B (zh) * | 2014-10-27 | 2018-10-23 | 国家电网公司 | 一种特高压线路操作电压等值风速的修正方法 |
CN106844795A (zh) * | 2015-12-03 | 2017-06-13 | 甘肃省电力公司风电技术中心 | 一种风电场内部风向场确认方法 |
CN105447255B (zh) * | 2015-12-03 | 2018-06-12 | 哈尔滨工业大学 | 双馈风电场动态等值方法、建模方法及*** |
US10385829B2 (en) | 2016-05-11 | 2019-08-20 | General Electric Company | System and method for validating optimization of a wind farm |
CN106951614A (zh) * | 2017-03-09 | 2017-07-14 | 南京理工大学 | 风力机模拟器的考虑阻尼补偿的兆瓦级风机模拟方法 |
CN107194097B (zh) * | 2017-05-27 | 2021-01-12 | 中国大唐集团科学技术研究院有限公司 | 基于风电场气动模拟和风速风向数据的分析方法 |
CN108879799B (zh) * | 2018-08-10 | 2021-10-29 | 三一重能有限公司 | 风电场有功功率控制方法、装置及电子设备 |
CN110492438A (zh) * | 2019-07-10 | 2019-11-22 | 山东科技大学 | 一种基于瞬时频率差的风电场主变压器保护方法 |
CN110968958B (zh) * | 2019-12-11 | 2020-12-01 | 贵州电网有限责任公司 | 一种基于单机等值与选择模态分析的风电场等值建模方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101414317A (zh) * | 2008-11-27 | 2009-04-22 | 东北电力大学 | 一种处理风电场静态功率等值分散性问题的等效风速法 |
CN101852172A (zh) * | 2010-03-09 | 2010-10-06 | 山东科技大学 | 一种考虑风电场内尾流影响的风电机组输入风速计算方法 |
CN101882167A (zh) * | 2009-11-25 | 2010-11-10 | 中国电力科学研究院 | 一种大规模风电集中接入电网的风电场等值建模方法 |
-
2010
- 2010-11-30 CN CN201010576662.5A patent/CN102012956B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101414317A (zh) * | 2008-11-27 | 2009-04-22 | 东北电力大学 | 一种处理风电场静态功率等值分散性问题的等效风速法 |
CN101882167A (zh) * | 2009-11-25 | 2010-11-10 | 中国电力科学研究院 | 一种大规模风电集中接入电网的风电场等值建模方法 |
CN101852172A (zh) * | 2010-03-09 | 2010-10-06 | 山东科技大学 | 一种考虑风电场内尾流影响的风电机组输入风速计算方法 |
Non-Patent Citations (3)
Title |
---|
潮流计算中风电场的等值;曹娜 等;《电网技术》;20060531;第30卷(第9期);全文 * |
风电场动态分析中风速模型的建立及应用;曹娜 等;《中国电机工程学报》;20071231;第27卷(第36期);全文 * |
风速波动时风电场动态特性分析;曹娜 等;《太阳能学报》;20090430;第30卷(第4期);引言第2-4段、第1节第1段、第2节、第3.1-3.3节、图2、4、6、7 * |
Also Published As
Publication number | Publication date |
---|---|
CN102012956A (zh) | 2011-04-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102012956B (zh) | 一种考虑风电场输入风速和风向随机波动的风电场等值方法 | |
CN103887815B (zh) | 基于运行数据的风电场参数辨识及动态等值方法 | |
Hansen et al. | Dynamic modelling of wind farm grid interaction | |
CN101882167B (zh) | 一种大规模风电集中接入电网的风电场等值建模方法 | |
CN102013702B (zh) | 外部电网故障时并网风电场的动态等值方法 | |
CN102760191B (zh) | 基于转速分群的双馈机组风电场等值建模***及方法 | |
CN103094920B (zh) | 一种直驱式风电机组风电场等值方法 | |
CN104269867B (zh) | 一种节点扰动功率转移分布均衡度分析方法 | |
CN103138256A (zh) | 一种新能源电力消纳全景分析***及方法 | |
CN105703396B (zh) | 一种双馈风力发电场多机表征等值建模方法 | |
CN106774276A (zh) | 风电场自动发电控制***测试平台 | |
CN105069236B (zh) | 考虑风电场节点空间相关性的广义负荷联合概率建模方法 | |
Wang et al. | Dynamic equivalent modeling for wind farms with DFIGs using the artificial bee colony with K-means algorithm | |
CN105762816B (zh) | 一种抑制次同步相互作用的方法 | |
CN102034008B (zh) | 风速波动下定速风电机组风电场的变参数电容补偿等值建模方法 | |
Wu et al. | Development of an equivalent wind farm model for frequency regulation | |
CN104820741A (zh) | 兼顾风场分散性与机组差异性的风电场动态等值方法 | |
CN106777499A (zh) | 一种双馈异步风力发电机组的整机动态建模方法 | |
CN104617578B (zh) | 一种含风电场电力***的可用输电能力的获取方法 | |
CN107769227B (zh) | 一种适用于次同步研究的风电场等值建模方法 | |
CN103853888A (zh) | 一种适用于无功补偿的海上风电***等值建模方法 | |
CN110571794A (zh) | 一种适用于双馈风机风电场的暂态模型等值计算方法 | |
CN104037806B (zh) | 一种基于风力发电机组基本模型的电力***潮流计算方法 | |
Ni et al. | Cross‐gramian‐based dynamic equivalence of wind farms | |
CN109657380A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160120 |