CN103761349B - 一种基于风电机组概率同调性的风电场等值建模方法 - Google Patents
一种基于风电机组概率同调性的风电场等值建模方法 Download PDFInfo
- Publication number
- CN103761349B CN103761349B CN201310323508.0A CN201310323508A CN103761349B CN 103761349 B CN103761349 B CN 103761349B CN 201310323508 A CN201310323508 A CN 201310323508A CN 103761349 B CN103761349 B CN 103761349B
- Authority
- CN
- China
- Prior art keywords
- wind
- turbine set
- wind energy
- energy turbine
- equivalent
- 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
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)
Abstract
本发明涉及一种基于风电机组概率同调性的风电场等值建模方法,本发明基于电力***实际运行中各类短路故障发生的概率和风电场中风资源状况的统计结果,通过时域仿真采集不同短路故障类型和风速组合条件下风力发电机转速在暂态过程中的动态变化情况,并依据此信息使用two‑step分类法对风电机组聚合分类。在汇总所有风电机组分群结果的基础上,根据机群划分结果出现的概率确定风电场概率等值模型,并对风力发电机和集电***等参数进行等值。本发明提出的风电场等值建模方法充分考虑了不同短路故障和机组运行工况对同调机群划分结果的影响,有效避免了工程应用中对同调机群变化的实时监测环节,提高了风电场等值模型的工程应用价值。
Description
技术领域
本发明涉及一种对风电场并网外特性的风电场等值建模方法。
背景技术
目前,适用于离线暂态稳定性分析的风电场等值模型大多基于电力***同调等值思想建立。传统电力***动态等值中,发电单元间的同调性可依据暂态过程中发电机转子摇摆曲线提供的信息进行判别,并且研究指出,机组的同调性主要与三个因素有关:1)故障扰动类型;2)发电机组运行工况;3)网络拓扑结构。文献[1]考虑了运行工况对同调性的影响,依据风电机组的输入风速进行机群划分。但是,风速具有随机性和波动性,基于风速差异的机组分群结果会随输入风速大小或风向的改变而变化。文献[2]虽然考虑了风速大小和风向随机波动的特性,但在实际应用中仍然需要动态识别同调机群来调用相应的风电场等值模型。文献[3]依据反映桨距角动作情况的故障前机组有功功率、机端电压和风速进行机群划分,文献[4]依据风电机组故障初始时刻的转速划分机群,但是都没有考虑故障扰动因素对机组同调性的影响。文献[5]在基于风电场内所有风电机组机端短路电压全都为零的假设条件下,仅通过近似计算得到的风力发电机故障切除时刻这一采样点下的转速对定速风电机组进行分群,难以准确、全面反映风电机组在暂态过程中的电压跌落差异和动态运行特性。在风电场等值建模研究过程中,考虑的基本上都是风电场公共节点处发生三相对称短路故障的情况,很少对单相短路或相间短路情况进行研究,这与实际电力***运行中故障类型发生概率的大小恰好相反。因此,需要更深入***地研究风电机组在不同故障类型下的暂态响应特性以及它们之间的差异。
参考文献:
(1)Fernandez L M,Saenz J R,Jurado F.Dynamic models of wind farms with fixed speed wind turbines[J].Renewable Energy,2006,31(8):1203-1230。
(2)曹娜,于群.一种考虑风电场输入风速和风向随机波动的风电场等值方法: 中国,CN102012956A[P].2011.04.13。
(3)苏勋文.基于桨距角动作情况的双馈机组风电场动态等值建模方法:中国,CN101937483A[P],2011.01.05。
(4)张保会,王小立,郝治国.基于转速分群的双馈机组风电场等值建模***及方法:中国,CN102760191A[P]。2012.10.31。
(5)苏勋文.风电场动态等值建模方法研究[D]。北京:华北电力大学,2010。
发明内容
本发明为了提高风电场等值模型在工程应用中的实用性和通用性,避免对风电机组同调性的动态识别环节。本发明提出了针对不同故障类型和风电场全工况条件下的风电场等值建模方法。
本发明是通过以下技术方案实现的:
一种基于风电机组概率同调性的风电场等值建模方法,包括以下步骤:
(1)建立了风电场风资源统计库和电力***短路故障库;
(2)基于风电场的网络拓扑结构和参数值,在仿真软件下搭建风电场详细模型;
(3)读取风资源统计库中的风速信息,利用尾流效应计算推导出各单台风电机组的输入风速,完成动态仿真前模型的潮流初始化;
(4)从电力***故障库中选取某一电网侧短路故障类型,采集时域仿真过程中单台风电机组在故障初始时刻、故障切除时刻、故障切除后0.2s和0.4s时刻上的发电机转子角速度,形成风力发电机的转子角速度向量,作为同调机组划分的标准;
(5)依据发电机转子角速度向量,使用two-step分类法先将风电机组压缩为子聚类集合,再利用层级聚类方法将子聚类进一步合并为更大的聚类;每个机群用一台等值机组表征,并且根据选取的风速和故障类型发生概率确定该组合条件下的风电场等值模型概率;
(6)重复上述第(3)步至第(5)步,完成风资源统计库和故障类型库中所有组合情况下的风电机组同调机群划分,进而根据不同组合条件下机组同调性的概率确定风电场等值概率模型;
(7)根据风电场等值概率模型计算风速的等值;风力发电机参数的等值;风电 场内部集电***线路的等值;风电机组机端无功补偿电容的等值。
一种基于风电机组概率同调性的风电场等值建模方法,所述的仿真软件为DIgSILENT PowerFactory仿真软件。
一种基于风电机组概率同调性的风电场等值建模方法,所述步骤(3)中利用尾流效应计算推导出各单台风电机组输入风速的计算方法利用以下公式计算得出,
式中:vj(t)为作用在任意台风力机上的输入风速;vj0为没有经过任何塔影效应作用在第j台风力机的输入风速;vw-k(xkj)为考虑机组间尾流效应时第k台风力机作用在第j台风力机上的尾流风速;βk=(Ashad-jk/Arot-j)表示在第j台风力机处第k台风力机的投影面积与第j台风力机面积的比值;n为风力机总台数。
一种基于风电机组概率同调性的风电场等值建模方法,所述步骤(5)中根据组合条件下输入风速的统计概率fi wind和故障类型发生概率fj fault确定此次机组分群结果的概率fij wtg,即:
fij wtg=fi wind×fj fault
其中,风速的统计概率fi wind为风速大小和风向概率的乘积。
一种基于风电机组概率同调性的风电场等值建模方法,步骤(7)中所述计算风速的等值采用以下方法计算得出:
式中:A为等效风力机的扫风面积;Ai、cpi、vi分别为机群中第i台机组的风力机扫风面积、风能利用系数和输入风速;m为机群中风电机组台数;下标eq表示等值模型的参数。
一种基于风电机组概率同调性的风电场等值建模方法,步骤(7)中所述计算风力发电机参数的等值采用以下方法计算得出:
利用参数变换方法计算等值风力发电机的参数,计算公式为:
其中虚拟阻抗Zfi可通过下式计算:
式中:n为机群中风电机组台数;Zs,Zm,Zr,Zf分别为发电机的定子阻抗、激磁阻抗、转子阻抗和虚拟阻抗;S、P分别为风电机组的额定功率和有功功率;ST为机端变压器的额定容量;H、K和D分别为惯性时间常数、轴系刚度系数和轴系阻尼系数;下标i表示第i台风电机组;下标eq表示等值模型的参数。
一种基于风电机组概率同调性的风电场等值建模方法,步骤(7)中所述计算风电场内部集电***线路的等值采用以下方法计算得出:
线路阻抗的并联化分解从线路末端向首段逐步进行,即首先从Zl2开始。Zl2,1、Zl2,2的计算公式为:
修正已经完成并联化变换的风电机组等值阻抗,即:
Zn′=Zn+Zli,n
依据修正后的风电机组支路等值阻抗进行下一个线路阻抗的分解,在纯并联结构 下按照以下公式计算线路的等值阻抗,
式中:Pi为第i台风电机组额定容量;Zl′i′为纯并联结构下第i台风电机组支路的线路阻抗;m为等值风电机组台数。
一种基于风电机组概率同调性的风电场等值建模方法,步骤(7)中所述风电机组机端无功补偿电容的等值采用以下方法计算得出:
通过引进变参数电容对风电机组无功补偿电容进行等值,由于风电机组吸收的无功功率与有功功率之间存在非线性函数关系,即
Qe=-Q0-Q1Pe-Q2Pe 2=f(Pe)
式中:Q0、Q1、Q2可由实验测得,在保证等值模型有功功率与等值前相等的前提下,等值前后风电场中机组吸收的无功功率不相等,通过引入变参数电容来等效由于等值前后风场内机组吸收不同的无功功率而导致的无功补偿差,变参数电容值为;
其中:
Pei为第i台风电机组的有功功率;m为等值风电机组台数,最终的补偿电容为:
式中:Ci为第i台风电机组的机端补偿电容值。
本发明提出的基于风电机组概率同调性的风电场等值建模方法根据故障初始时刻、故障切除时刻以及故障切除后0.2s和0.4s时风电机组的转速进行机群划分,完全计及了故障扰动类型、风电机组运行工况和网络拓扑结构这三个因素对机组同调性的影响;根据风电场实际风资源信息和不同故障类型统计概率确定的风电场概率模型符合工程应用要求;通过概率统计建立的风电场等值模型很好 的解决了由于风速随机性和波动性引起的机组同调性不确定的问题,增加了风电场等值模型的可信度。
本发明建立了风电场风资源统计库和电力***短路故障库。对风电场全年的风速状况进行统计,包括风速大小和风向信息,如图3所示,通过对风资源情况的统计可以直观的确定风电机组运行于各种工况下的概率。电力***短路故障库的建立基于电力***的运行经验。电力***运行过程中,各类短路发生的次数在短路总次数中所占的比例不同,通过合理推算确定各类短路故障发生的概率,如表2所示。记录时域仿真过程中风力发电机转子在多个采样点上的角速度大小,包括故障前、故障切除时刻和故障切除后各个时间段上的值,并且以向量的形式构建机群划分标准。与仅采用故障初始或故障切除时刻发电机转子角速度的标准相比,该机群划分标准可以更全面、准确的反映风力发电机转子角速度在暂态过程中的变化,提高了风电机组机群划分精度。聚类算法的选取。本发明选用的two-step分类法第一步完成简单数据处理,将原始输入数据压缩为可管理的子聚类集合,第二步使用层级聚类方法将子聚类进一步合并为更大的聚类。Two-step法能够为训练数据自动估计最佳聚类数,解决了传统k-means聚类算法中聚类数无法自动寻优以及分群结果对初始聚类中心敏感的问题,因此可进一步提高风电机组机群划分精度。风电场内部概率同调机组的确立。汇总各种短路故障和风电场全工况组合条件下的风电机组机群划分结果,通过对所有分群结果的概率统计确定概率同调机群,包括概率统计结果最大的机群数和每个机群中包含的风电机组。本发明中使用基于网络变换前后风电机组机端电压不变为原则的集电***参数等值方法,通过对集电***网络纯并联化处理,更准确的反映了线路阻抗对风电机组机端电压的影响。与基于等值前后线路功率损耗相等原则的等值方法相比,其更好的表征了由于电气距离对风电机组行为相关性的影响,适用于风电场机组聚合等值建模时不同馈线上风电机组聚合的情况。依据风电机组故障初始时刻、故障切除时刻以及故障切除后的转速向量作为同调机群的分群指标,一是转速比风速更能准确表征风电机组的暂态特性,二是转速向量也可以更完整的反映机组的暂态特性。所述方法综合了电力***中多种不同故障类型对机组同调性的影响,建立的风电场等值模型不再需要对风电机组间的同调性进行实时监测。
表1风电场电气参数表
表2短路故障概率统计表
表3风速7m/s、风向180。工况时不同故障类型条件下机组分群结果
表4风电场机组分群结果概率统计表
表5风电机组输入风速表
表6风电机组转速表
表7基于two-step分类法的分群结果
附图说明
图1风电场概率等值建模流程图。
图2风电场接入IEEE14节点***图。
图3风电场风资源统计图。
图4多台电机参数变换法等值示意图。
图5风电机组串并联混合结构下的集电***等效电路图
图6风电机组纯并联网络变换示意图
图7风电场概率等值模型图
图8风电场PCC点处有功功率变化曲线图
图9风电场PCC点处无功功率变化曲线图
图10风电场PCC点处有功功率变化曲线图
图11风电场PCC点处无功功率变化曲线图
具体实施方式
参见附图1-11所示。本发明的实施方式包括两部分内容,一是风电机组概率同调性的研究,二是风电场概率等值模型参数计算。
步骤a:风电机组概率同调性研究,具体包括以下内容:
第一步;根据气象数据惯例,将00-3600的风向均分为16个区域,风速的间隔步长取为1m/s,对风电场的风资源状况进行概率统计,统计结果如图3所示。提取单台风电机组和集电***线路参数。
第二步:基于风电场的网络拓扑结构和参数值,在DIgSILENT PowerFactory仿真软件下搭建风电场详细模型,风电场公共节点接入IEEE14节点测试***的14号母线,如图2所示。其中,风电机组的机械传动***采用计及轴系刚度和阻尼系数的二质量块模型,即考虑在故障过程中轴系能量的变化及其对风力发电机转速的影响。
第三步:读取风资源统计库中的风速信息,利用尾流效应计算推导出各单台风电机组的输入风速,完成动态仿真前模型的潮流初始化。考虑尾流效应条件下作用在任意台风电机组的输入风速为:
式中:vj(t)为作用在任意台风力机上的输入风速;vj0为没有经过任何塔影效应作用在第j台风力机的输入风速;vw-k(xxj)为考虑机组间尾流效应时第k台风力机作用在第j台风力机上的尾流风速;βk=(Ashad-jk/Arot-j)表示在第j台风力 机处第k台风力机的投影面积与第j台风力机面积的比值;n为风力机总台数。第四步:从电力***故障库中选取某一故障类型,进行暂态时域仿真,观察采集每台风电机组在故障初始时刻、故障切除时刻、故障切除后0.2s和0.4s时刻对应的发电机转子角速度ω0、ω1、ω3和ω4,形成风电机组转子角速度向量{ω1,ω2,ω3,ω4},作为风电机组机群划分的标准。
第五步:基于风电机组转子角速度向量,使用two-step分类法先将风电机组压缩为子聚类集合,再利用层级聚类方法将子聚类进一步合并为更大的聚类。根据该组合条件下输入风速的统计概率fi wind和故障类型发生概率fj fault确定此次机组分群结果的概率fij wtg,即:
fij wtg=fi wind×fj fault
其中,风速的统计概率fi wind为风速大小和风向概率的乘积。
第六步:重复上述第三步至第五步,完成所有输入风速和故障类型组合条件下的风电机组机群划分。
对经过two-step法聚合分类的风电机组分群结果进行统计分析,本发明中提出的依据故障初始时刻、故障切除时刻以及故障切除后0.2s和0.4s风电机组转速构成的机群分类指标可以计及故障扰动类型对风电机组同调性的影响,如表3所示。可以看出,虽然风电机组的初始运行工况相同,但不同类型故障的扰动导致了不同的风电机组机群划分结果。
按照以上六步对two-step法聚合分类后的风电机组分群结果进行统计分析,总共有1344组分群结果,其中有108个不同的组别,包括不同的机群数或相同机群数下包含有不同的风电机组。依据图3中风电场全年的风资源统计信息和表2中电力***运行中不同故障类型发生的概率统计,得到如表4所示的风电机组分群结果概率统计值,其中机群划分结果A出现的概率最大,为依据本发明中风电场等值建模方法确定的概率等值模型。
步骤b:风电场概率等值模型参数计算,其中包括以下内容:
第一步:根据公式二计算等效风速。
式中:A为等效风力机的扫风面积;Ai、cpi、vi分别为机群中第i台机组的风力机扫风面积、风能利用系数和输入风速;m为机群中风电机组台数;下标eq表示等值模型的参数。
第二步:利用参数变换方法计算等值风力发电机的参数,其变换过程见图4,计算公式为:
其中虚拟阻抗Zfi可通过下式计算:
式中:n为机群中风电机组台数;Zs,Zm,Zr,Zf分别为发电机的定子阻抗、激磁阻抗、转子阻抗和虚拟阻抗;S、P分别为风电机组的额定功率和有功功率;ST为机端变压器的额定容量;H、K和D分别为惯性时间常数、轴系刚度系数和轴系阻尼系数;下标i表示第i台风电机组;下标eq表示等值模型的参数。
第三步:基于网络变换前后风电机组机端电压不变的原则对线路参数进行等值,将风电机组串并联混合形式的集电***网络结构转换成纯并联结构。以风电场中第一条馈线为例,变换前的等效电路如图5所示,纯并联网络变换的方法即将线路的阻抗分解成多个阻抗串联到风电机组支路中,如图6所示。
线路阻抗的并联化分解从线路末端向首段逐步进行,即首先从Zl2开始。Zl2,1、 Zl2,2的计算公式为:
修正已经完成并联化变换的风电机组等值阻抗,即:
Zn′=Zn+Zli,n 公式六
依据修正后的风电机组支路等值阻抗进行下一个线路阻抗的分解。在纯并联结构下按照公式七计算线路的等值阻抗。
式中:Pi为第i台风电机组额定容量;Zli″为纯并联结构下第i台风电机组支路的线路阻抗;m为等值风电机组台数。
第四步;通过引进变参数电容对风电机组无功补偿电容进行等值。由于风电机组吸收的无功功率与有功功率之间存在非线性函数关系,即
Qe=-Q0-Q1Pe-Q2Pe 2=f(Pe) 公式八
式中:Q0、Q1、Q2可由实验测得。在保证等值模型有功功率与等值前相等的前提下,等值前后风电场中机组吸收的无功功率不相等。通过引入变参数电容来等效由于等值前后风场内机组吸收不同的无功功率而导致的无功补偿差,变参数电容值为;
其中:
Pei为第i台风电机组的有功功率;m为等值风电机组台数。因此,最终的补偿电容为:
式中:Ci为第i台风电机组的机端补偿电容值。
为验证本发明中依据风电机组转速向量作为分群指标的有效性,在第0.2s时风电场PCC点处一回输电线路发生三相短路,150ms后故障切除,风电机组的输入风速如表5所示,表6中为采集得到的机组转速。分别选取转速向量、故障初始时刻转速和故障切除时刻转速作为机组分群指标,利用two-step法对风电机组聚合分类,结果如表7所示。图8、9给出了风电场详细模型、基于转速向量建立的两机等值模型、基于故障初始时刻转速建立的两机等值模型以及基于故障切除时刻转速建立的三机等值模型在风电场PCC点处的有功功率、无功功率动态响应过程。从图8、9中可以看出,基于本发明中转速向量分群指标建立的风电场等值模型的有功功率、无功功率与详细模型几乎相同,误差小于基于故障初始时刻转速或故障切除时刻转速分群指标建立的风电场等值模型,说明本发明提出的依据转速向量进行机群划分的方法更加准确。
当风电场输入风速大小5m/s、风向315。时,风电场PCC点处0.2s发生单相接地短路,150ms后故障切除,基于转速向量的风电机组分群结果为表4中C所示。从图10、11可以看出,风电场概率等值模型的有功功率、无功功率与详细模型保持一致,误差大小符合工程应用要求,说明本发明建立的风电场概率等值模型即使在分群结果不同的运行工况和故障类型组合条件下,仍然具有较高的模型精度。
Claims (8)
1.一种基于风电机组概率同调性的风电场等值建模方法,其特征在于包括以下步骤:
(1)建立了风电场风资源统计库和电力***短路故障库;
(2)基于风电场的网络拓扑结构和参数值,在仿真软件下搭建风电场详细模型;
(3)读取风资源统计库中的风速信息,利用尾流效应计算推导出各单台风电机组的输入风速,完成动态仿真前模型的潮流初始化;
(4)从电力***故障库中选取某一电网侧短路故障类型,采集时域仿真过程中单台风电机组在故障初始时刻、故障切除时刻、故障切除后0.2s和0.4s时刻上的发电机转子角速度,形成风力发电机的转子角速度向量,作为同调机组划分的标准;
(5)依据发电机转子角速度向量,使用two-step分类法先将风电机组压缩为子聚类集合,再利用层级聚类方法将子聚类进一步合并为更大的聚类;每个机群用一台等值机组表征,并且根据选取的风速和故障类型发生概率的乘积确定该条件下的风电场等值模型概率;
(6)重复上述第(3)步至第(5)步,完成风资源统计库和故障类型库中所有组合情况下的风电机组同调机群划分,进而根据不同组合条件下机组同调性的概率确定风电场等值概率模型;
(7)根据风电场等值概率模型计算风速的等值;风力发电机参数的等值;风电场内部集电***线路的等值;风电机组机端无功补偿电容的等值。
2.根据权利要求1所述的一种基于风电机组概率同调性的风电场等值建模方法,其特征在于所述的仿真软件为DIgSILENT PowerFactory仿真软件。
3.根据权利要求1所述的一种基于风电机组概率同调性的风电场等值建模方法,其特征在于:所述步骤(3)中利用尾流效应计算推导出各单台风电机组输入风速的计算方法利用以下公式计算得出,
式中:vj(t)为作用在任意台风力机上的输入风速;vj0(t)为没有经过任何塔影效应作用在第j台风力机的输入风速;vw-k(xkj,t)为考虑机组间尾流效应时第k台风力机作用在第j台风力机上的尾流风速;βk=(Ashad-jk/Arot-j)表示在第j台风力机处第k台风力机的投影面积Ashad-jk与第j台风力机面积Arot-j的比值;n为风力机总台数。
4.根据权利要求1所述的一种基于风电机组概率同调性的风电场等值建模方法,其特征在于:所述步骤(5)中根据组合条件下输入风速的统计概率fi wind和故障类型发生概率fj fault确定此次机组分群结果的概率fij wtg,即:
fij wtg=fi wind×fj fault
其中,风速的统计概率fi wind为风速大小和风向概率的乘积。
5.根据权利要求1所述的一种基于风电机组概率同调性的风电场等值建模方法,其特征在于,步骤(7)中所述风速的等值采用以下方法计算得出:
式中:A为等效风力机的扫风面积;Ai、cpi、vi分别为机群中第i台机组的风力机扫风面积、风能利用系数和输入风速;m为机群中风电机组台数;下标eq表示等值模型的参数;veq为等值风力机的输入风速,cp_eq为等值风力机的风能利用系数。
6.根据权利要求1所述的一种基于风电机组概率同调性的风电场等值建模方法,其特征在于,步骤(7)中所述风力发电机参数的等值采用以下方法计算得出:利用参数变换方法计算等值风力发电机的参数,计算公式为:
其中虚拟阻抗Zfi可通过下式计算:
式中:n为机群中风电机组台数;Zsi,Zmi,Zri,Zfi分别为第i台发电机的定子阻抗、激磁阻抗、转子阻抗和虚拟阻抗;Si、Pi分别为第i台风电机组的额定功率和有功功率;STi为第i台风电机组端口变压器的额定容量;Hi、Ki和Di分别为第i台风电机组的惯性时间常数、轴系刚度系数和轴系阻尼系数;下标i表示第i台风电机组;下标eq表示等值模型的参数。
7.根据权利要求1所述的一种基于风电机组概率同调性的风电场等值建模方法,其特征在于,步骤(7)中所述风电场内部集电***线路的等值采用以下方法计算得出:
线路阻抗的并联化分解从线路末端向首段逐步进行,即首先从Zl2开始;Zl2,1、Zl2,2的计算公式为:
修正已经完成并联化变换的风电机组等值阻抗,即:
Z′n=Zn+Zli,n
依据修正后的风电机组支路等值阻抗进行下一个线路阻抗的分解,在纯并联结构下按照以下公式计算线路的等值阻抗,
式中:Pi为第i台风电机组额定容量;Z″li为纯并联结构下第i台风电机组支路的线路阻抗;m为等值风电机组台数;Z1为第1台风电机组的阻抗,Z2为第2台风电机组的阻抗;Zn为第n台风电机组的阻抗;Zl2为第2台风电机组支路的线路阻抗;Zli,n表示由第i台风电机组支路的线路阻抗分解成的多个阻抗串联到第n台风电机组支路上的阻抗。
8.根据权利要求1所述的一种基于风电机组概率同调性的风电场等值建模方法,其特征在于,步骤(7)中所述风电机组机端无功补偿电容的等值采用以下方法计算得出:
通过引进变参数电容对风电机组无功补偿电容进行等值,由于风电机组吸收的无功功率与有功功率之间存在非线性函数关系,即
式中:Q0、Q1、Q2可由实验测得,为非线性函数中无特殊含义的三个不同的常数;Pe为风电机组吸收的有功功率,Qe和f(Pe)为风电机组吸收的无功功率;在保证等值模型有功功率与等值前相等的前提下,等值前后风电场中机组吸收的无功功率不相等,通过引入变参数电容来等效由于等值前后风场内机组吸收不同的无功功率而导致的无功补偿差,变参数电容值为;
其中:
Pei为第i台风电机组的有功功率;Qei和f(Pei)为第i台风电机组吸收的无功功率;
为等值后风电机组吸收的无功功率;f0为频率,U为变参数电容的电压;
m为等值风电机组台数,最终的补偿电容为:
式中:Ci为第i台风电机组的机端补偿电容值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310323508.0A CN103761349B (zh) | 2013-07-29 | 2013-07-29 | 一种基于风电机组概率同调性的风电场等值建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310323508.0A CN103761349B (zh) | 2013-07-29 | 2013-07-29 | 一种基于风电机组概率同调性的风电场等值建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103761349A CN103761349A (zh) | 2014-04-30 |
CN103761349B true CN103761349B (zh) | 2016-08-24 |
Family
ID=50528585
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310323508.0A Expired - Fee Related CN103761349B (zh) | 2013-07-29 | 2013-07-29 | 一种基于风电机组概率同调性的风电场等值建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103761349B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104504285A (zh) * | 2015-01-06 | 2015-04-08 | 合肥工业大学 | 一种计及撬棒保护的双馈风电场等值建模方法 |
CN104795812B (zh) * | 2015-04-07 | 2017-04-26 | 中国南方电网有限责任公司电网技术研究中心 | 一种考虑定子阻尼效应并采用变量代换的慢同调分区法 |
CN105260553A (zh) * | 2015-10-26 | 2016-01-20 | 中国南方电网有限责任公司电网技术研究中心 | 考虑电压动态特性下的同调机组判别方法 |
CN106991229B (zh) * | 2017-03-30 | 2020-08-21 | 华北电力大学 | 一种针对复杂拓扑的风电场等值建模方法 |
CN108446515B (zh) * | 2018-03-30 | 2020-05-22 | 河海大学 | 一种基于双馈风力发电机组短路电流特征分析的风电场等值方法 |
CN108804742A (zh) * | 2018-04-08 | 2018-11-13 | 西南交通大学 | 一种大规模风电场等值模型及其构建方法和应用 |
CN109086527B (zh) * | 2018-08-03 | 2022-10-04 | 合肥工业大学 | 一种基于风电机组运行状态的实用化等值建模方法 |
CN109816261A (zh) * | 2019-01-31 | 2019-05-28 | 广东电网有限责任公司电网规划研究中心 | 海上半直驱式风电场等值方法、***和装置 |
CN110210170B (zh) * | 2019-06-19 | 2022-10-18 | 河北工业大学 | 一种大规模风电机群等值小信号模型建模方法 |
CN115062563B (zh) * | 2022-08-18 | 2022-11-18 | 南京理工大学 | 基于空间变化的风力机三维尾流风速计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101764561A (zh) * | 2008-12-25 | 2010-06-30 | 上海电气自动化设计研究所有限公司 | 一种风力发电电网频率的仿真电路及其方法 |
CN101882167A (zh) * | 2009-11-25 | 2010-11-10 | 中国电力科学研究院 | 一种大规模风电集中接入电网的风电场等值建模方法 |
CN102592026A (zh) * | 2012-01-12 | 2012-07-18 | 甘肃省电力公司风电技术中心 | 一种双馈风机等效模拟的仿真建模方法 |
CN103034763A (zh) * | 2012-12-12 | 2013-04-10 | 上海市电力公司 | 大容量风电场风机聚合模型建立和仿真方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008155830A1 (ja) * | 2007-06-20 | 2008-12-24 | Fujitsu Limited | 回路シミュレーションモデル生成装置、回路シミュレーションモデル生成プログラム、回路シミュレーションモデル生成方法及び回路シミュレーション装置 |
CN103282904B (zh) * | 2010-12-29 | 2016-01-20 | 株式会社村田制作所 | 等效电路模型,程序以及存储介质 |
-
2013
- 2013-07-29 CN CN201310323508.0A patent/CN103761349B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101764561A (zh) * | 2008-12-25 | 2010-06-30 | 上海电气自动化设计研究所有限公司 | 一种风力发电电网频率的仿真电路及其方法 |
CN101882167A (zh) * | 2009-11-25 | 2010-11-10 | 中国电力科学研究院 | 一种大规模风电集中接入电网的风电场等值建模方法 |
CN102592026A (zh) * | 2012-01-12 | 2012-07-18 | 甘肃省电力公司风电技术中心 | 一种双馈风机等效模拟的仿真建模方法 |
CN103034763A (zh) * | 2012-12-12 | 2013-04-10 | 上海市电力公司 | 大容量风电场风机聚合模型建立和仿真方法 |
Non-Patent Citations (2)
Title |
---|
"基于同调机群识别的动态等值方法的工程应用";杨靖萍等;《电网技术》;20050930;第29卷(第17期);第68-71页 * |
蔡国伟等."基于两步聚类法的受扰机群同调性分析".《中国电力》.2013,第46卷(第2期),第46-49页. * |
Also Published As
Publication number | Publication date |
---|---|
CN103761349A (zh) | 2014-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103761349B (zh) | 一种基于风电机组概率同调性的风电场等值建模方法 | |
Zou et al. | A fuzzy clustering algorithm-based dynamic equivalent modeling method for wind farm with DFIG | |
CN101937483B (zh) | 基于桨距角动作情况的双馈机组风电场动态等值建模方法 | |
CN102760191B (zh) | 基于转速分群的双馈机组风电场等值建模***及方法 | |
CN102136036B (zh) | 适于电力***小扰动稳定分析的双馈风电场等值建模方法 | |
CN103094920B (zh) | 一种直驱式风电机组风电场等值方法 | |
Li et al. | Transient stability analysis of wind turbines with induction generators considering blades and shaft flexibility | |
CN105678033B (zh) | 一种适用于电磁暂态仿真的风电场等值建模方法 | |
CN104036073B (zh) | 适于有功功率特性分析的双馈风电场动态等值建模方法 | |
CN108667005B (zh) | 一种计及新能源影响的电网静动态结合脆弱性评估方法 | |
CN103558768B (zh) | 一种基于风电场内风速分布特性的等值建模方法 | |
Xu et al. | Influence of different flexible drive train models on the transient responses of DFIG wind turbine | |
CN106383947B (zh) | 风电场集电网络动态等值参数的快速获取方法 | |
CN102013702B (zh) | 外部电网故障时并网风电场的动态等值方法 | |
Al-Bayati et al. | Aggregated models of wind farms: Current methods and future trends | |
CN102034008B (zh) | 风速波动下定速风电机组风电场的变参数电容补偿等值建模方法 | |
CN105243230B (zh) | 一种计及卸荷电路的直驱风电场的等值建模方法 | |
CN105576654B (zh) | 直驱风电场等值方法和*** | |
CN104820741A (zh) | 兼顾风场分散性与机组差异性的风电场动态等值方法 | |
CN104504285A (zh) | 一种计及撬棒保护的双馈风电场等值建模方法 | |
CN106410862B (zh) | 基于有功恢复斜率校正的风电场单机等值方法 | |
CN103986189A (zh) | 一种储能型双馈风电场简化模型建模方法 | |
CN104242336B (zh) | 基于发电单元相互作用系数的新能源电站等效建模方法 | |
CN111027179A (zh) | 一种计及辅助调频服务的双馈风电场的等值建模方法 | |
CN104617578A (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: 20160824 Termination date: 20190729 |
|
CF01 | Termination of patent right due to non-payment of annual fee |