CN106894950A - 一种基于风电机组传动特性的功率特性验证方法 - Google Patents
一种基于风电机组传动特性的功率特性验证方法 Download PDFInfo
- Publication number
- CN106894950A CN106894950A CN201710138372.4A CN201710138372A CN106894950A CN 106894950 A CN106894950 A CN 106894950A CN 201710138372 A CN201710138372 A CN 201710138372A CN 106894950 A CN106894950 A CN 106894950A
- Authority
- CN
- China
- Prior art keywords
- power
- output
- wind
- curve
- current transformer
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 57
- 238000012795 verification Methods 0.000 title claims abstract description 17
- 238000010977 unit operation Methods 0.000 claims abstract description 6
- 238000013461 design Methods 0.000 claims abstract description 5
- 230000008878 coupling Effects 0.000 claims description 19
- 238000010168 coupling process Methods 0.000 claims description 19
- 238000005859 coupling reaction Methods 0.000 claims description 19
- 230000008569 process Effects 0.000 claims description 13
- 238000005259 measurement Methods 0.000 claims description 12
- 230000001360 synchronised effect Effects 0.000 claims description 9
- 239000013256 coordination polymer Substances 0.000 claims description 7
- 238000009795 derivation Methods 0.000 claims description 7
- 238000010408 sweeping Methods 0.000 claims description 3
- 230000002159 abnormal 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
- 235000008434 ginseng Nutrition 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 6
- 238000011156 evaluation Methods 0.000 abstract description 4
- 238000011160 research Methods 0.000 abstract description 4
- 238000012360 testing method Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000005211 surface analysis Methods 0.000 description 1
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D7/00—Controlling wind motors
- F03D7/02—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor
- F03D7/028—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor controlling wind motor output power
-
- 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/72—Wind turbines with rotation axis in wind direction
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Wind Motors (AREA)
Abstract
本发明涉及一种基于风电机组传动特性的功率特性验证方法,通过获取一定时间段的机组运行数据,构建能量传输模型进行功率特性验证,本发明为以机组能量传输过程为研究对象,基于各部件能量传输过程中的损耗,建立一个通用的能量传输模型,计算推导风速,拟合实际功率曲线,将实际曲线与机组设计时的理论曲线进行对比,通过一定的方法定量计算两条曲线的差异,形成一整套机组功率曲线验证方法和技术路线,为机组性能评价提供了一个评价准确、可行性高、实用性强的方法。
Description
技术领域
本发明涉及一种风电机组功率评价方法,特别是一种基于风电机组传动特性的功率特性验证方法。
背景技术
目前国内外有关研究机构已经开展了风电机组功率特性测试的相关研究,2005年IEC发布了IEC 61400-12-1标准,规定了风力发电机组功率特性测试的相关标准,包括测试准备、测试程序、后期数据处理等,2012年该标准被引入中国并作为国家标准施行,标准号为GB/T 18451.2-2012《风力发电机组功率特性测试》。由于标准的相关实施条件要求较高,比如要求在机组前某些特定扇区、位置树立测风塔连续测量一定时段的风资源,很难满足现场需求。
目前国内在进行功率曲线验证时,较通用的是采用一机一验证的办法,具体就是用超声波或激光雷达测风仪参照IEC和行业的有关标准测定轮毂前一定距离的来风风速,同时读取机组相同时段的功率,经过一定数据处理后进行曲线拟合,同时参考机组历史数据进行一定的修正。该方法具有简单灵活、实施条件简单的优点,但其不足也很明显,比如:
1)未考虑风电机组空气密度的影响,曲线不够准确;
2)用历史数据修正功率曲线,使曲线带有机组运行特性,降低可信度和推广度;
3)一机一测定增加工作量,加大了工作量,不利于大面积推广;
4)在具体测量手段、数据处理方面难以实现完全共识。
发明内容
本发明所要解决的技术问题是提供一种以机组能量传输过程为研究对象,基于各部件能量传输过程中的损耗,建立一个通用的能量传输模型,计算推导风速,拟合实际功率曲线,将实际曲线与机组设计时的理论曲线进行对比,通过一定的方法定量计算两条曲线的差异,形成一整套机组功率曲线验证方法和技术路线,为机组性能评价提供了一个评价准确、可行性高、实用性强的方法。
本发明解决上述技术问题所采用的技术方案是:
一种基于风电机组传动特性的功率特性验证方法,包括如下步骤:
S1)获取一定时间段的机组运行数据,包括机组出力、发电机转速、叶轮转速和桨距角;
S2)能量传输模型构建,所述能量传输模型构建依次包括风轮捕捉机械功率的计算、推导风速的计算和模型修正;
S21)风轮捕捉机械功率的计算:根据变流器、发电机、齿轮箱的效率曲线,寻取某一功率输出情况下,各传动部件的工作效率,得到风轮捕捉的机械功率;
S22)推导风速的计算:计算出风轮捕捉的机械功率后,利用风能利用系数和风速之间关系,推导实际空气密度下功率对应的风速,再将该风速换算到标准大气密度条件下;
S23)模型修正:架设测风仪器,同步获取测定的风速V2和机组出力P0,根据步骤S21和S22计算机组出力P0下的推导风速V1,比较V1和V2差异并调整传输模型,直至差异缩小到设定范围,完成能量传输模型构建;
S3)功率特性验证,所述功率特性验证包括实际功率曲线拟合、偏差计算和机组功率评价;
S31)实际功率曲线拟合:以传输模型为基础,输入P0得到V1,拟合V1-P0的对应实际功率曲线;
S32)偏差计算:将实际功率曲线与机组设计时的理论功率曲线放在一个坐标轴下,计算曲线的偏差e;
S33)根据偏差e对机组功率进行评价。
进一步地,所述步骤S1中还包括如下步骤:
S11)对采集的数据源进行甄别,利用网格法筛选机组正常运行状态时刻的有功功率。
作为优选,所述步骤S1中还包括如下步骤:
S12)所述机组运行数据还包括温度和压强,通过四分位法将温度和压强中的异常数据剔除,利用空气密度与温度和压强的函数关系求出空气密度。
进一步地,所述步骤S21中通过折半查找法确定变流器的效率和其输入功率,变流器的输入功率作为发电机的输出功率,通过折半查找法确定发电机的效率和其输入功率,发电机的输入功率作为主轴联轴器的输出功率,通过折半查找法确定主轴联轴器的效率和输入功率,主轴联轴器的输入功率作为齿轮箱的输出功率,通过折半查找法确定齿轮箱的效率和其输入功率,齿轮箱的输入功率即为风轮捕捉的机械功率。
作为优选,其特征在于:
所述步骤S21的具体计算过程如下:
网侧变流器输出功率为Pr wc、发电机定子输出功率为Ps、机侧变流器输出功率为Pr jc,Pr jc=sPs,其中s为转差率、机组出口功率为PCK,PCK=Pr+Ps、发电机输出功率为PG、联轴器输出功率为PLZQ、齿轮箱输出功率为PCLX、叶轮输出功率为P,齿轮箱、联轴器、发电机、变频器效率参数为η1、η2、η3、η4;
机组出口功率与发电机定子、网侧变流器输出功率之间的关系可由下面的公式表示:PCK=Pr wc+Ps;
发电机组有两种运行方式,当转子转速小于同步转速时s<0,转子侧馈入来自电网的功率,此时Pr jc<0;转子转速大于同步转速时s≥0,转子向电网馈出功率,此时Pr jc≥0;
若Pr jc≥0,根据发电机转速ω,求得转差率s,再假设Ps 0,机侧变流器输入功率为sPs,根据变流器运行效率曲线η4 0=η(sPs 0)计算出网侧变流器输出功率Pr wc=η(sPs 0)·sPs 0,则机组出力比较实测机组出力PCK与二者之间的差距,差距过大重新假设Ps 0,采用相同的步骤计算直至PCK二者接近到设定范围,求得Ps,发电机输出功率为PG=Ps+sPs;
若Pr jc<0,先假设Ps 0,机侧变流器输出至发电机转子的功率为sPs 0,假设网侧变流器输入功率Pr wc0,根据变流器运行效率曲线η4 0=η(Pr wc0),返回计算机侧变流器输入功率sPs 0′,比较sPs 0、sPs 0′差异,差异设定范围,重新假设Pr wc0,直至在设定范围内,迭代计算出η0以及网侧变流器输入功率Pr wc0后,计算机组出力比较实测PCK、差异,差异过大,重新假设Ps 0,过程与Pr jc≥0相同,求得Ps,发电机输出功率为PG=Ps+sPs;
再推算发电机输入机械功即联轴器输出功率PLZQ,假设联轴器输出功率初始值为PLZQ 0,求得发电机此出力下的效率则发电机输出功率比较迭代与已知的PG之间的差异,差异过大重新假设PLZQ 0,直至在设定范围内;根据同样的方法,推算出齿轮箱输出功率为PCLX、叶轮输出功率P。
作为优选,其特征在于:
所述步骤S22的具体计算过程如下:
机械功率P与风能捕捉系数CP和推导风速V1有如下关系:
机械功率P与风能捕捉系数CP和推导风速V1有如下关系:其中,P为风轮捕捉的机械功率、CP为风能捕捉系数、ρ为空气密度、A为扫风面积、V1为推导风速、θ为错风角;
CP与λ和β有关,λ与V1有关,将P转换为V1的唯一函数,计算出V1。
作为优选,其特征在于:
所述步骤S23的具体修正过程如下:
测风仪器数据覆盖从机组切入至满功率风速段;根据传输模型计算推导风速V1,比较V1、V2,计算相对误差为ε=(V1-V2)/V2,计算其均方误差为n为运行数据总组数,σ大于设定值时,在齿轮箱效率基础上加上一个调整系数α,即齿轮箱传动效率为α·η2,重新进行步骤S21至S23直至σ小于等于设定值。
作为优选,其特征在于:所述步骤S31中通过多项式拟合方法绘制出V1-P0的对应实际功率曲线。
作为优选,其特征在于:所述步骤S32中以切入风速为起始,切出风速为终止,1m/s的风速为步长,对比每点理论曲线和实际曲线的差值,计算曲线各点的平均差值e。
本发明同现有技术相比具有以下优点及效果:本发明以机组能量传输过程为研究对象,基于各部件能量传输过程中的损耗,建立一个通用的能量传输模型,计算推导风速,拟合实际功率曲线,将实际曲线与机组设计时的理论曲线进行对比,通过一定的方法定量计算两条曲线的差异,形成一整套机组功率曲线验证方法和技术路线,为机组性能评价提供了一个评价准确、可行性高、实用性强的方法。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为风电机组功率转化及流向示意图。
图2为本发明的流程图。
具体实施方式
下面结合实施例对本发明做进一步的详细说明,以下实施例是对本发明的解释而本发明并不局限于以下实施例。
实施例1:
本实施例包括如下步骤:
步骤101:按照一定比例均匀抽取对象机组,从SCADA系统中导出对象机组的风速、功率、压力、温度、发动机转速、错风角、变桨角度和机组运行状态以及测风塔导出的气压、温度、湿度等数据。根据需要确定抽取比例,一般的对于同型号、同配置、同期投役的风电机组,随机抽取一台机组作为研究机组;导出SCADA运行数据及测风数据,要求覆盖从切入风速至满功率风速所有运行段。
步骤102:通过四分位法将温度和压强中的异常数据剔除,利用空气密度与温度和压强的函数关系求出某时刻的空气密度;利用网格法将机组状态中的正常运行状态时刻的有功功率筛选出来。
所谓四分位法是一种数据分组方法,将全部数据按由小到大的顺序排列,分成四等份,处于前25%的数据称为第一四分位数,处于25%~50%的数据称为第二四分位数,处于50%~75%的数据称为第三四分位数,最后25%的数据称为第四分位数。
步骤201:风轮捕捉机械功率的计算:首先确定变流器的效率和其输入功率,再确定发电机的效率和其输入功率,以此类推,确定主轴联轴器、齿轮箱的效率和其输入功率,齿轮箱的输入功率即为风轮捕捉到的机械功率。计算过程的关键是通过部件输出功率迭代计算其效率和输入功率,最终得出齿轮箱的输入功率,具体过程如下:
设网侧变流器输出功率为Pr wc(定经网侧变流器流进电网方向为功率正向,反之为负向)、发电机定子输出功率为Ps、机侧变流器输出功率为Pr jc(Pr jc=sPs,其中s为转差率,方向定义与网侧变流器相同)、机组出口功率为PCK(PCK=Pr+Ps)、发电机输出功率为PG、联轴器输出功率为PLZQ、齿轮箱输出功率为PCLX、叶轮输出功率为P,设齿轮箱、联轴器、发电机、变频器效率参数为η1、η2、η3、η4。
机组出口功率与发电机定子、网侧变流器输出功率之间的关系可由下面的公式表示:
PCK=Pr wc+Ps (1)
双馈发电机组有两种运行方式,当转子转速小于同步转速时(s<0),转子侧馈入来自电网的功率,此时Pr jc<0;转子转速大于同步转速时(s≥0),转子向电网馈出功率,此时Pr jc≥0。以下阐述两种不同运行方式下发电机输出功率的计算方法。
若Pr jc≥0,通过SCADA系统,获取某一时刻一组运行数据,数据类型如步骤101所述。先根据发电机转速ω,求得转差率s,再假设Ps 0,则转子侧对应的输出功率(即机侧变流器输入功率)为sPs,根据变流器运行效率曲线η4 0=η(sPs 0)计算出网侧变流器输出功率Pr wc=η(sPs 0)·sPs 0,则机组出力比较实测机组出力PCK与二者之间的差距,差距过大重新假设Ps 0,采用相同的步骤计算直至PCK二者接近到接受范围内(边界根据具体情况确定,一般可定为0.1‰)。求得Ps,发电机转子输出功率也求得,则发电机输出功率为PG=Ps+sPs。
若Pr jc<0,同样的,先假设Ps 0,求得机侧变流器输出至发电机转子的功率为sPs 0,再假设网侧变流器输入功率Pr wc0,根据变流器运行效率曲线η4 0=η(Pr wc0),返回计算机侧变流器输入功率sPs 0′,比较sPs 0、sPs 0′差异,差异超过边界(一般可定为0.1‰),重新假设Pr wc0,直至在边界范围内,迭代计算出η0以及网侧变流器输入功率Pr wc0后,计算机组出力比较实测PCK、差异,差异过大,重新假设Ps 0,过程与Pr jc≥0相同。求得Ps,发电机转子输出功率也求得,则发电机输出功率为PG=Ps+sPs。
再推算发电机输入机械功即联轴器输出功率PLZQ,假设联轴器输出功率初始值为PLZQ 0,求得发电机此出力下的效率则发电机输出功率比较迭代与已知的PG之间的差异,差异过大重新假设PLZQ 0,直至在边界范围内(一般定0.1‰)。
根据同样的方法,推算出齿轮箱输出功率为PCLX、叶轮输出功率P。
步骤202:推导风速的计算:由于机舱对风与来风风向不是实时一致的,实际叶轮输出的机械功率P与风能捕捉系数CP和推导风速V1有如下关系:
其中,P为风轮捕捉的机械功率、CP为风能捕捉系数、ρ为空气密度、A为扫风面积、V1为推导风速、θ为错风角。
CP与λ和β有关,而λ又与V1有关,因此在做一定数学处理后将P转换为V1的唯一函数,从而计算出V1。计算过程如下:
1)先初始假设一个V1′,计算相应的λ′,再根据机组叶片的实时桨距角,按照公式(3)计算基于当前V1′下的CP′。
2)将CP′带入公式(2)中,计算V1″,比较初始假设值V1′和计算输出值V1″,差距大于设定值(一般定0.1‰),重新假设V1′,计算V1″,再计算二者差值,直至其小于设定值,得出满足条件的V1′,定其为对应于机械功率P的自由来流风速。
由于风速与密度有关,同样的机械功,密度越大,轮毂前风速越小,因此应将实际空气密度带入公式对风速进行修正。
步骤203:模型修正:架设激光雷达测风仪等测风仪器,同步获取雷达测风仪测定的风速V2和机组出力P0,要求数据覆盖从机组切入至满功率风速段。对获取的数据进行一定处理,根据传输模型计算V1,比较V1、V2,计算相对误差为ε=(V1-V2)/V2,计算其均方误差为(n为运行数据总组数),σ过大时分析传输模型的原因,调整传输模型,直至差异缩小到可接受范围,完成能量传输模型的构建。调整过程如下,在齿轮箱效率基础上加上一个调整系数α,即齿轮箱传动效率为α·η2,重新进行步骤201至203直至满足结束条件为止(一般的,结束条件定σ<0.1%)。
步骤301:功率曲线拟合:以传输模型为基础,输入P0得到V1,通过多项式拟合方法绘制出V1-P0拟合曲线。
步骤302:偏差计算:将实际曲线与理论功率曲线放在一个坐标轴下,以切入风速为起始,切出风速为终止,1m/s的风速为步长,对比每点理论曲线和实际曲线的差值,计算曲线各点的平均差值e。
步骤303:差异评价:将e按三个等级进行分级,分界点分别是ε1、ε2,ε1=0.02,ε2=0.05。定e≤ε1机组功率特性优秀、ε1<e≤ε2合格、e>ε2不合格。对于不合格机组,可以按照先软件后硬件的思路,从控制系统控制变量调节特性着手,再分析机组叶片、齿轮箱散热等方面分析曲线偏低的原因。
此外,需要说明的是,本说明书中所描述的具体实施例,其零、部件的形状、所取名称等可以不同。凡依本发明专利构思所述的构造、特征及原理所做的等效或简单变化,均包括于本发明专利的保护范围内。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,只要不偏离本发明的结构或者超越本权利要求书所定义的范围,均应属于本发明的保护范围。
Claims (9)
1.一种基于风电机组传动特性的功率特性验证方法,包括如下步骤:
S1)获取一定时间段的机组运行数据,包括机组出力、发电机转速、叶轮转速和桨距角;
S2)能量传输模型构建,所述能量传输模型构建依次包括风轮捕捉机械功率的计算、推导风速的计算和模型修正;
S21)风轮捕捉机械功率的计算:根据变流器、发电机、齿轮箱的效率曲线,寻取某一功率输出情况下,各传动部件的工作效率,得到风轮捕捉的机械功率;
S22)推导风速的计算:计算出风轮捕捉的机械功率后,利用风能利用系数和风速之间关系,推导实际空气密度下功率对应的风速,再将该风速换算到标准大气密度条件下;
S23)模型修正:架设测风仪器,同步获取测定的风速V2和机组出力P0,根据步骤S21和S22计算机组出力P0下的推导风速V1,比较V1和V2差异并调整传输模型,直至差异缩小到设定范围,完成能量传输模型构建;
S3)功率特性验证,所述功率特性验证包括实际功率曲线拟合、偏差计算和机组功率评价;
S31)实际功率曲线拟合:以传输模型为基础,输入P0得到V1,拟合V1-P0的对应实际功率曲线;
S32)偏差计算:将实际功率曲线与机组设计时的理论功率曲线放在一个坐标轴下,计算曲线的偏差e;
S33)根据偏差e对机组功率进行评价。
2.根据权利要求1所述的基于风电机组传动特性的功率特性验证方法,其特征在于:所述步骤S1中还包括如下步骤:
S11)对采集的数据源进行甄别,利用网格法筛选机组正常运行状态时刻的有功功率。
3.根据权利要求2所述的基于风电机组传动特性的功率特性验证方法,其特征在于,所述步骤S1中还包括如下步骤:
S12)所述机组运行数据还包括温度和压强,通过四分位法将温度和压强中的异常数据剔除,利用空气密度与温度和压强的函数关系求出空气密度。
4.根据权利要求1所述的基于风电机组传动特性的功率特性验证方法,其特征在于:所述步骤S21中通过折半查找法确定变流器的效率和其输入功率,变流器的输入功率作为发电机的输出功率,通过折半查找法确定发电机的效率和其输入功率,发电机的输入功率作为主轴联轴器的输出功率,通过折半查找法确定主轴联轴器的效率和输入功率,主轴联轴器的输入功率作为齿轮箱的输出功率,通过折半查找法确定齿轮箱的效率和其输入功率,齿轮箱的输入功率即为风轮捕捉的机械功率。
5.根据权利要求4所述的基于风电机组传动特性的功率特性验证方法,其特征在于:
所述步骤S21的具体计算过程如下:
网侧变流器输出功率为Pr wc、发电机定子输出功率为Ps、机侧变流器输出功率为Pr jc,Pr jc=sPs,其中s为转差率、机组出口功率为PCK,PCK=Pr+Ps、发电机输出功率为PG、联轴器输出功率为PLZQ、齿轮箱输出功率为PCLX、叶轮输出功率为P,齿轮箱、联轴器、发电机、变频器效率参数为η1、η2、η3、η4;
机组出口功率与发电机定子、网侧变流器输出功率之间的关系可由下面的公式表示:PCK=Pr wc+Ps;
发电机组有两种运行方式,当转子转速小于同步转速时s<0,转子侧馈入来自电网的功率,此时Pr jc<0;转子转速大于同步转速时s≥0,转子向电网馈出功率,此时Pr jc≥0;
若Pr jc≥0,根据发电机转速ω,求得转差率s,再假设Ps 0,机侧变流器输入功率为sPs,根据变流器运行效率曲线η4 0=η(sPs 0)计算出网侧变流器输出功率Pr wc=η(sPs 0)·sPs 0,则机组出力比较实测机组出力PCK与二者之间的差距,差距过大重新假设Ps 0,采用相同的步骤计算直至PCK二者接近到设定范围,求得Ps,发电机输出功率为PG=Ps+sPs;
若Pr jc<0,先假设Ps 0,机侧变流器输出至发电机转子的功率为sPs 0,假设网侧变流器输入功率Pr wc0,根据变流器运行效率曲线η4 0=η(Pr wc0),返回计算机侧变流器输入功率sPs 0′,比较sPs 0、sPs 0′差异,差异设定范围,重新假设Pr wc0,直至在设定范围内,迭代计算出η0以及网侧变流器输入功率Pr wc0后,计算机组出力比较实测PCK、差异,差异过大,重新假设Ps 0,过程与Pr jc≥0相同,求得Ps,发电机输出功率为PG=Ps+sPs;
再推算发电机输入机械功即联轴器输出功率PLZQ,假设联轴器输出功率初始值为PLZQ 0,求得发电机此出力下的效率则发电机输出功率比较迭代与已知的PG之间的差异,差异过大重新假设直至在设定范围内;根据同样的方法,推算出齿轮箱输出功率为PCLX、叶轮输出功率P。
6.根据权利要求4所述的基于风电机组传动特性的功率特性验证方法,其特征在于:
所述步骤S22的具体计算过程如下:
机械功率P与风能捕捉系数CP和推导风速V1有如下关系:
机械功率P与风能捕捉系数CP和推导风速V1有如下关系:其中,P为风轮捕捉的机械功率、CP为风能捕捉系数、ρ为空气密度、A为扫风面积、V1为推导风速、θ为错风角;
CP与λ和β有关,λ与V1有关,将P转换为V1的唯一函数,计算出V1。
7.根据权利要求4所述的基于风电机组传动特性的功率特性验证方法,其特征在于:
所述步骤S23的具体修正过程如下:
测风仪器数据覆盖从机组切入至满功率风速段;根据传输模型计算推导风速V1,比较V1、V2,计算相对误差为ε=(V1-V2)/V2,计算其均方误差为n为运行数据总组数,σ大于设定值时,在齿轮箱效率基础上加上一个调整系数α,即齿轮箱传动效率为α·η2,重新进行步骤S21至S23直至σ小于等于设定值。
8.根据权利要求1所述的基于风电机组传动特性的功率特性验证方法,其特征在于:所述步骤S31中通过多项式拟合方法绘制出V1-P0的对应实际功率曲线。
9.根据权利要求1所述的基于风电机组传动特性的功率特性验证方法,其特征在于:所述步骤S32中以切入风速为起始,切出风速为终止,1m/s的风速为步长,对比每点理论曲线和实际曲线的差值,计算曲线各点的平均差值e。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710138372.4A CN106894950B (zh) | 2017-03-09 | 2017-03-09 | 一种基于风电机组传动特性的功率特性验证方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710138372.4A CN106894950B (zh) | 2017-03-09 | 2017-03-09 | 一种基于风电机组传动特性的功率特性验证方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106894950A true CN106894950A (zh) | 2017-06-27 |
CN106894950B CN106894950B (zh) | 2019-02-12 |
Family
ID=59185914
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710138372.4A Active CN106894950B (zh) | 2017-03-09 | 2017-03-09 | 一种基于风电机组传动特性的功率特性验证方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106894950B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107330183A (zh) * | 2017-06-29 | 2017-11-07 | 华北电力大学 | 一种基于运行数据的风电利用率计算方法 |
CN108397358A (zh) * | 2018-01-10 | 2018-08-14 | 内蒙古久和能源装备有限公司 | 一种风力发电机组多点效率检测方法和系统 |
CN111159885A (zh) * | 2019-12-27 | 2020-05-15 | 大唐(赤峰)新能源有限公司 | 一种基于曲线拟合的风机参数离散率异常分析方法 |
CN112131703A (zh) * | 2020-08-13 | 2020-12-25 | 国网新源张家口风光储示范电站有限公司 | 一种基于matlab功能模块的风机运行特性分析系统 |
CN114219299A (zh) * | 2021-12-16 | 2022-03-22 | 重庆大学 | 考虑区间不确定性的多风轮风电机组传动效率评估方法 |
CN114992047A (zh) * | 2022-07-13 | 2022-09-02 | 华电电力科学研究院有限公司 | 一种风力发电机组控制方法及相关组件 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1744058A1 (en) * | 2004-05-07 | 2007-01-17 | Mitsubishi Denki Kabushiki Kaisha | Wind power generation evaluation system and prediction control service system for wind power generator |
CN101545953A (zh) * | 2009-04-24 | 2009-09-30 | 重庆大学 | 小型风力发电机性能测试方法及其专用测试平台 |
CN103699804A (zh) * | 2013-12-31 | 2014-04-02 | 西安热工研究院有限公司 | 一种实时监测风电机组运行时的功率出力偏差的方法 |
CN104156575A (zh) * | 2014-07-28 | 2014-11-19 | 国家电网公司 | 基于测风塔数据外推法的风电场理论功率计算方法 |
-
2017
- 2017-03-09 CN CN201710138372.4A patent/CN106894950B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1744058A1 (en) * | 2004-05-07 | 2007-01-17 | Mitsubishi Denki Kabushiki Kaisha | Wind power generation evaluation system and prediction control service system for wind power generator |
CN101545953A (zh) * | 2009-04-24 | 2009-09-30 | 重庆大学 | 小型风力发电机性能测试方法及其专用测试平台 |
CN103699804A (zh) * | 2013-12-31 | 2014-04-02 | 西安热工研究院有限公司 | 一种实时监测风电机组运行时的功率出力偏差的方法 |
CN104156575A (zh) * | 2014-07-28 | 2014-11-19 | 国家电网公司 | 基于测风塔数据外推法的风电场理论功率计算方法 |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107330183A (zh) * | 2017-06-29 | 2017-11-07 | 华北电力大学 | 一种基于运行数据的风电利用率计算方法 |
CN107330183B (zh) * | 2017-06-29 | 2020-04-17 | 华北电力大学 | 一种基于运行数据的风电利用率计算方法 |
CN108397358A (zh) * | 2018-01-10 | 2018-08-14 | 内蒙古久和能源装备有限公司 | 一种风力发电机组多点效率检测方法和系统 |
CN111159885A (zh) * | 2019-12-27 | 2020-05-15 | 大唐(赤峰)新能源有限公司 | 一种基于曲线拟合的风机参数离散率异常分析方法 |
CN111159885B (zh) * | 2019-12-27 | 2024-01-02 | 大唐(赤峰)新能源有限公司 | 一种基于曲线拟合的风机参数离散率异常分析方法 |
CN112131703A (zh) * | 2020-08-13 | 2020-12-25 | 国网新源张家口风光储示范电站有限公司 | 一种基于matlab功能模块的风机运行特性分析系统 |
CN114219299A (zh) * | 2021-12-16 | 2022-03-22 | 重庆大学 | 考虑区间不确定性的多风轮风电机组传动效率评估方法 |
CN114219299B (zh) * | 2021-12-16 | 2024-04-19 | 重庆大学 | 考虑区间不确定性的多风轮风电机组传动效率评估方法 |
CN114992047A (zh) * | 2022-07-13 | 2022-09-02 | 华电电力科学研究院有限公司 | 一种风力发电机组控制方法及相关组件 |
CN114992047B (zh) * | 2022-07-13 | 2024-05-28 | 华电电力科学研究院有限公司 | 一种风力发电机组控制方法及相关组件 |
Also Published As
Publication number | Publication date |
---|---|
CN106894950B (zh) | 2019-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106894950A (zh) | 一种基于风电机组传动特性的功率特性验证方法 | |
CN104794347B (zh) | 一种风电机组机舱传递函数分区拟合方法 | |
CN106150904A (zh) | 一种风力发电机组偏航系统控制性能优化方法及系统 | |
CN105649875B (zh) | 风力发电机组的变桨控制方法及装置 | |
CN110761958B (zh) | 风力发电机组的叶片失速诊断方法及装置 | |
Odabaee et al. | CFD simulation of a supercritical carbon dioxide radial-inflow turbine, comparing the results of using real gas equation of estate and real gas property file | |
CN107002636A (zh) | 用于估计风速,包括计算针对叶片扭转调节的桨距角的方法 | |
CN103825298B (zh) | 一种定速风电机组驱动系统参数辨识方法 | |
CN103219725A (zh) | 一种基于实时运行数据的风电场等值建模方法 | |
CN105298748A (zh) | 提高风电机组发电效率的方法 | |
CN110808608A (zh) | 规模化新能源参与受端电网调频调压能力评估方法及系统 | |
CN108869174A (zh) | 一种非线性建模的风力发电机叶片固有频率工况补偿方法 | |
CN104888396A (zh) | 一种检测消防水泵在快速启动过程中瞬态故障的方法 | |
CN111412114B (zh) | 一种基于定子电流包络谱的风电机组叶轮不平衡检测方法 | |
CN111311021A (zh) | 一种风电场理论功率预测方法、装置、设备和存储介质 | |
CN115549139A (zh) | 一种新能源发电和负荷混合模型辨识建模方法 | |
CN108278184A (zh) | 基于经验模态分解的风电机组叶轮不平衡监测方法 | |
CN106894949B (zh) | 基于环境因素的风机功率信号反馈方法 | |
CN107339194B (zh) | 最佳叶尖速比控制方法、数据处理方法及装置、系统 | |
CN104454381A (zh) | 一种优化风力机风能捕获过程的功率增发的方法 | |
CN111756039B (zh) | 基于概率统计的新能源电力系统惯量估计方法 | |
CN108691730A (zh) | 风力发电机系统及风能数据误差修正方法和偏航控制方法 | |
Hiron et al. | The Performance of a Three-blades Fish-ridge Turbine in an Oscillating Water Column System for Low Waves | |
CN112600199A (zh) | 一种基于转子动能的风电机组等效转动惯量在线评估方法 | |
Bao et al. | A data-driven approach for identification and compensation of wind turbine inherent yaw misalignment |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |