CN115905201A - 一种风机限电损失计算方法 - Google Patents
一种风机限电损失计算方法 Download PDFInfo
- Publication number
- CN115905201A CN115905201A CN202211522353.9A CN202211522353A CN115905201A CN 115905201 A CN115905201 A CN 115905201A CN 202211522353 A CN202211522353 A CN 202211522353A CN 115905201 A CN115905201 A CN 115905201A
- Authority
- CN
- China
- Prior art keywords
- fan
- power
- agc
- time
- calculating
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 claims abstract description 26
- 238000007781 pre-processing Methods 0.000 claims abstract description 7
- 238000010248 power generation Methods 0.000 claims description 11
- 210000001061 forehead Anatomy 0.000 claims description 5
- 230000005611 electricity Effects 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 239000000843 powder Substances 0.000 claims description 2
- 238000012423 maintenance Methods 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 230000005856 abnormality Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002864 sequence alignment Methods 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
Images
Landscapes
- Wind Motors (AREA)
Abstract
本发明公开了一种风机限电损失计算方法,风电场SCADA系统实时数据和历史采集和存储、数据预处理、基于时间窗的滑动平均风速计算、风机限功率实时状态判定、风机理论发电量计算、风机及风场限功率损失电量计算。该方法利用SCADA采集的模拟量帮助更加准确性的计算和分析风机限功率状态,可更加广泛的应用于大量存量老旧风机无准确限功率状态标志的风电企业的限功率损失电量计算,并提升常规风电企业的限功率损失计算准确性;该方法还可以结合风电场升压站AGC信号准确计算风电场的AGC限功率损失电量,为风电场精细化限电损失管理提供依据。
Description
技术领域
本发明属于电力控制领域,具体涉及一种风机限电损失计算方法。
背景技术
风电能源企业在日常运维和管理过程中,需要实时的统计风电场整体及各台风机的运行损失电量,以针对性的对问题风机进行检修和维护,最大程度的提升风电场发电效能。其中,准确计算风机运行过程中的限功率状态的电量损失是个难点,也非常重要,其可以为风电场精确运维管理提供巨大帮助,可使得风电企业更加准确的开展人员绩效考核,精确的计算非自身因素导致的电量损失以免除自身责任等。
目前风电企业的限电损失计算大多依赖风机的状态码进行限功率状态判定,其计算的准确性极度依赖风机状态码的准确性。实际运行过程中,很多风机缺乏限功率状态码或者虽然能够提供限功率状态码,但是其计算准确性不高,计算策略偏向保守导致状态误判。此外,目前大多风机都无法提供更加精细化的AGC限功率状态码,导致无法计算得到风场的AGC限功率损失电量,从而无法帮助风电企业实现更加精细化的风场运行管理和运维。
因此,需要设计开发一种新型风电企业风机限电损失计算方法,提升常规风电企业的限功率损失计算准确性,并且能够准确计算风电场的AGC限功率损失电量,为风电场精细化运维提供依据。
公开号为CN114297251A的专利发明了一种风电机组损失电量计算方法,其中对限电损失的计算依赖风机的AGC信号和自限电信号来进行状态判定。但是该方法极度依赖风机PLC给出的限功率状态码的准确性,实际应用中,几乎很少有风机可以给出AGC限功率状态码,很多老旧的风机甚至连基本的限功率状态码都没有,因此无法基于状态码去判定风机的限功率状态乃至进一步判定AGC限功率状态。此外,该方案中,风机的理论功率根据风机实际功率曲线计算得到,而准确的应发功率的获取还依赖准确的风速数据。在大多时候处于的非稳定工况下,传统的10min平均风速的准确性较低,且数据获取周期10min较长,导致时效性也较差,因此基于此平均风速得到的理论功率的准确性也较差,从而导致损失电量计算误差较大。
发明内容
为解决现有技术的问题,本发明提供一种风机限电损失计算方法,首先分别针对风电场内各风机,执行如下步骤S1至步骤S6,获得风机的有功功率差值、以及限功率运行状态,进而获得各风机分别所对应的有功功率差值、以及限功率运行状态;然后执行步骤S7至S10,完成所有风机总的AGC自限电损失电量统计;
步骤S1、采集风机状态码、瞬时风速、有功功率、有功发电量、风机桨叶角对应历史时长的数据量,获得风机各指定类型数据对应历史时长的连续数据量;
步骤S2、对风机各指定类型数据对应历史时长的连续数据量进行预处理并更新,获得预设固定时间间隔的离散数据,进一步获得风机各指定类型数据对应历史时长的离散数据量;
步骤S3、根据离散数据的瞬时风速序列、预设时间窗,结合预设时间窗滑动时长,对该瞬时风速序列进行划分,获得各个时间窗对应的瞬时风速序列,并计算各个时间窗对应的瞬时风速序列的平均速度;进一步获得滑动时间窗平均风速序列;
步骤S4、基于风机各指定类型数据对应历史时长的连续数据量包括风机桨叶角、风机有功功率的历史数据,按预设条件进行预处理,并绘制风机实际功率曲线,再根据风机对应的滑动时间窗平均风速序列,获得风机的理论有功功率值;根据风机实时数据、风机理论有功功率值,计算其有功功率差值;
步骤S5、根据风机状态码,将风机分类为:第一类风机、第二类风机,之后执行步骤S6;
步骤S6、根据风机状态码、实时有功功率、风机桨叶角,以及有功功率差值,分别按照预设条件判定第一类风机、以及第二类风机的限功率运行状态;
步骤S7、根据步骤S6获得的第一类风机、第二类风机中各个风机的限功率状态、以及各个风机的理论有功功率值,计算各个风机的理论发电量,进一步获得所有风机总限功率损失电量;
步骤S8、采集风场升压站AGC设定值信号和各个风机有功功率数据,计算获得风场升压站AGC设定值与风机总有功功率的差值e,并按预设条件判断风场AGC限功率状态;
步骤S9、在风场AGC限功率状态下,并且在目标时间段内,针对步骤S6获得的限功率状态风机,获得AGC限功率状态下风机时间段集合,然后对该集合内的风机计算其AGC限功率损失电量,并求和获得总AGC限功率损失电量。
进一步地,前述的S1中,利用风电场SCADA系统采集风机的数据。
进一步地,前述的步骤S2中,包括如下子步骤:
S201、对风机各指定类型数据对应历史时长的连续数据量存储与数据库,并对超过预设范围的数据进行剔除,并且进行数据更新;
S202、采用临近估计法将更新的数据进行字段统一处理,包括分别将时间戳、测点名称、测点值、测点所属风机名称、测点所属风场名称统一。
进一步地,前述的步骤S3具体为:
根据预设时间窗ΔT,对瞬时风速序列{Vi,i=1,2...a…m....b...n},按如下公式计算时间窗瞬时风速序列的滑动平均风速:
V平均(T-1)=∑(V1,...,Vm)/m;
其中,Vm为T1+ΔT时刻的风速,a,b,m为时间值,n为总时长。
进一步地,前述的步骤S4包括如下子步骤:
S401、按如下公式计算各个时间窗对应的瞬时风速序列的平均速度:
V平均(T-i)=∑(Va,...,Vb)/(b-a+1),
S402、按如下公式计算风机滑动平均风速序列:
{V平均(T-X),X=1,2..i...N},
其中,N为时间窗序号。
进一步地,前述的步骤S5包括如下子步骤:
S501、将风机滑动平均风速序列{V平均(T-X),X=1,2..i...N}代入风机实际功率曲线P=f(V),计算风机的理论有功功率值P理论;
S502、根据风机实时有功功率值P实际,以及P理论,按如下公式计算风机有功功率差值:
ΔP=P理论-P实际。
进一步地,前述的步骤S7包括如下子步骤:
S701、根据风机理论有功功率P理论,按如下公式计算各个风机的理论发电量:
Q理论=∫P理论dt;
S702、根据风机的限功率状态时间段集合{TXi,i=1,2...n},得到每个限功率时间段Txi的开始时间和结束时间,计算时间段TXi的风机实际发电量Q实际,如下式:
Q实际=μ*(Q末-Q初),
其中,Q初为时间段Txi开始时间的有功电量,Q末为时间段Txi结束时间的有功电量,μ为风机预设电量斜率参数;
S703、按如下公式计算n个限功率时间段Txi的风机总限功率损失电量Q总损:
Q总损=∑Q总损i=∑(Q理论i-Q实际i)=∑(∑Q理论i5/10min-∑Q实际i5/10min),
其中,i=1,2...n。
进一步地,前述的步骤S8包括如下子步骤:
S801、计算风场所有风机总有功功率如下式:
P总=∑Pi,
其中,Pi为第i个风机的有功功率,i=1,2...N风机;风机数为N风机;
S802、按如下公式计算风场所有风机总装机容量:
P额总=∑Pi额定,其中,Pi额定为第i个风机的额定有功功率,i=1,2...N风机;
S803、如下公式计算风场总有功功率的差值:
e=PAGC-P总。
进一步地,前述的步骤S9包括如下子步骤:
S901、根据AGC限功率损失电量QAGC损i,按如下公式计算风场的总的AGC限功率损失电量:
QAGC总损=∑QAGC损i,i=1,2...,Nx;
S902、按如下公式计算所有风机总的自限电损失电量:
Q自限电总损=Q总损-QAGC总损。
附图说明
图1是风机和风场限功率损失电量计算方法流程图。
图2是滑动时间窗示意图。
图3是风场AGC限功率损失电量计算流程图。
具体实施方式
为了更了解本发明的技术内容,特举具体实施例并配合所附图式说明如下。
在本发明中参照附图来描述本发明的各方面,附图中示出了许多说明性实施例。本发明的实施例不局限于附图所述。应当理解,本发明通过上面介绍的多种构思和实施例,以及下面详细描述的构思和实施方式中的任意一种来实现,这是因为本发明所公开的构思和实施例并不限于任何实施方式。另外,本发明公开的一些方面可以单独使用,或者与本发明公开的其他方面的任何适当组合来使用。
如图1所示,一种风机限电损失计算方法包括如下步骤:
步骤S1、数据采集与存储。风电场SCADA采集风机实时数据包括:风机状态码、瞬时风速、有功功率、有功发电量、风机桨叶角等,采集过去一年10min平均风速、10min平均有功功率、桨叶角、风机状态码的历史数据,获得风机各指定类型数据对应历史时长的连续数据量;可以基于标准TCP协议开发历史数据接口将上述采集到的历史数据存入历史数据库,也可以开发实时数据接口实时采集上述实时数据,并存入实时数据库。
步骤S2、数据的更新与预处理。实时数据库保持数据为过去30天,超过时效的可进行删除。通过数据转存模块,定期的(如每隔24h)将实时数据库中的数据复制保存到历史数据库中。历史数据库保持数据为过去一年,历史数据库对超过一年的历史数据进行删除。
(1)异常值剔除:将采集到的数据进行合理性判定,对于显著超过合理范围的数据进行丢弃,异常判定规则如下:
风机状态码:为空值、乱码、或超出合理范围的异常标识符;
风速(m/s):小于0或大于50;
桨叶角(°):小于0或大于95;
有功功率(kW):大于150%P额定或小于-30%P额定;其中P额定为风机的额定有功功率;
有功发电量(kWh):小于0或大于1.2*P额定*Tq;其中Tq为电量计算时长,单位h;
(2)时间序列对齐:采集到的各个测点的时间戳是各不相同的,需要按照整数分钟时刻进行调整对齐,得到固定间隔时间(例如1min)的整齐的所有数据值。为了简化计算过程,模拟量和开关量的时间对齐计算方法均采用临近估计法。即每个标准时刻的值,采用离其时间最近的历史值进行等效替代。
本发明采用了基于时间窗滑动平均法求取平均风速,如图2所示:
步骤S3、根据离散数据的瞬时风速序列、预设时间窗,结合预设时间窗滑动时长,对该瞬时风速序列进行划分,获得各个时间窗对应的瞬时风速序列,并计算各个时间窗对应的瞬时风速序列的平均速度;进一步获得滑动时间窗平均风速序列;时间窗宽度和时间窗间隔可灵活配置,设时间窗宽度ΔT配置为5min,时间窗间隔L配置为1min。具体算法如下:
采集的瞬时风速值按照时间序列分布为{Vi,i=1,2...m...n},其中T1时刻为时间窗ΔT的起始时间,其对应的风速为V1,T1+ΔT时刻的风速为Vm。则第一个时间窗(T-1)的滑动平均风速为:
V平均(T-1)=∑(V1,...,Vm)/m,
将时间窗往后平移时间间隔L,得到第二时间窗(T-2),其起始时间为T1+L,对应的风速为Va(1<a<m),T1+L+ΔT时刻的风速为Vb(m<b)。则第二时间窗(T-2)的滑动平均风速为:
V平均(T-2)=∑(Va,...,Vb)/(b-a+1)
则各个时间窗对应的瞬时风速序列的平均速度:
V平均(T-i)=∑(Va,...,Vb)/(b-a+1),
以此不断重复,继而得到了基于时间窗的滑动平均风速序列{V平均(T-X),X=1,2...N},N为时间窗序号。
与原始瞬时风速比,该基于时间窗的滑动平均风速时间隔为L(如1min),相比传统的10min平均风速计算法,其时效性提高了10倍。此外,该方法通过时间窗内的多瞬时风速点进行平均计算,也可实现常规10min平均风速计算对风速随机误差的降低效果。
需要注意的是,对于部分极端工况,风速的突变和跳跃可能会非常大,在剔除异常值之后,对有效数据进行上述基于时间窗的滑动平均风速计算,仍然可能得到随机误差较大的风速数据。此时可以通过略微调整时间窗的宽度和时间间隔进行计算,但需要注意,时间窗的宽度不能太大,否则会增大时域范围的累积误差。此时,可采用二阶时间窗滑动平均计算法,即基于第一次滑动平均计算后的风速数据基础上进行二次滑动平均计算,得到最终的结果可有效消除较大的随机误差。
将基于时间窗滑动平均法计算得到的平均风速结合功率曲线得到的理论功率与常规10min平均风速结合功率曲线得到的理论功率对比,比较基准为风机状态为正常运行时的实际功率趋势,发现基于时间窗滑动平均法计算的功率差值为基于常规10min平均法计算的功率差值的1/2~2/3,因此可确认能够显著提升理论功率的计算准确性。
步骤S4、基于风机各指定类型数据对应历史时长的连续数据量包括风机桨叶角、风机有功功率、风速的历史数据,按预设条件进行预处理,并绘制风机实际功率曲线,再根据风机对应的滑动时间窗平均风速序列,获得风机的理论有功功率值;根据风机实时数据、风机理论有功功率值,计算其有功功率差值;具体为:
首先获取历史数据库中的10min平均风速、桨叶角和10min平均有功功率数据,将10min平均有功功率<0.5kW、10min平均风速小于3m/s以及桨叶角大于3°的历史时间段所有数据剔除,将剩下的10min平均风速和10min平均有功功率数据拟合绘制成风机实际功率曲线P=f(V);
然后将风机滑动平均风速序列{V平均(T-X),X=1,2..i...N}代入风机实际功率曲线P=f(V),计算风机的理论有功功率值P理论;
最后根据风机实时有功功率值P实际,以及P理论,按如下公式计算风机有功功率差值:
ΔP=P理论-P实际。
步骤S5、根据风机状态码是否存在限功率运行状态,将风机分类为:第一类风机、第二类风机,随后执行步骤S6;
步骤S6、根据风机状态码、实时有功功率、风机桨叶角,以及有功功率差值,分别按照预设条件判定第一类风机、第二类风机的限功率运行状态;
第一类风机限功率运行状态判定条件为:
(1)风机状态码解析为正常运行;
(2)实际有功功率>0.5kW;
(3)风机桨叶角>3°;(定桨风机无该条件)
(4)理论和实际有功功率差值ΔP>k*P额定,其中,k为可调系数,一般可设置为0.05~0.15之间,可根据实际情况取最佳值;
(5)持续时长>3min;(可灵活配置,一般不超过10min)
以上5个条件须全部同时满足(定桨风机无第3个满足条件),可判定第一类风机为限功率运行状态。
上述计算频率为1min。
第二类风机限功率运行状态判定条件为:
条件1:风机状态码解析为限功率运行;
条件2:
(1)风机状态码解析为正常运行;
(2)实际有功功率>0.5kW;
(3)风机桨叶角>3°;(定桨风机无该条件)
(4)理论和实际有功功率差值ΔP>k*P额定,其中,P额定为风机的额定有功功率;k为可调系数,一般可设置为0.05~0.15之间,可根据实际情况取最佳值;
(5)持续时长>3min;(可灵活配置,一般不超过10min)
若条件1和条件2满足其中一个,可判定第二类风机为限功率运行状态。其中,条件2满足指的是条件2中所有条件同时满足。
上述计算频率为1min。
步骤S7、根据步骤S6获得的第一类风机、第二类风机中各个风机的限功率状态、以及各个风机的理论有功功率值,计算各个风机的理论发电量,进一步获得所有风机总限功率损失电量;
S701、根据风机理论有功功率P理论,按如下公式计算各个风机的理论发电量:
Q理论=∫P理论dt;
S702、根据风机实际发电量和理论发电量以及限功率运行状态时间段,计算限功率损失电量,步骤如下:
1)确定目标时间段内的所有风机为限功率状态的时间段集合{TXi,i=1,2...n},并得到每个限功率时间段TXi的开始时间和结束时间。
2)计算目标时间段内所有的限功率时间段Txi的风机实际发电量Q实际,可根据风机的有功电量测点终值与初始值的差值计算得到。计算公式如下:Q实际=μ*(Q末-Q初),其中,Q初为时间段TXi开始时间的有功电量,Q末为时间段Txi结束时间的有功电量,μ为风机预设电量斜率参数,一般为预设固定值。
实际计算过程中,可每间隔5min或10min计算一次该时段内的实际发电量,各限功率时间段Txi内的总的风机实际发电量为每5min或每10min的实际发电量之和。(非整数时段的实际发电量可根据每个限功率时间段Txi的开始时间和结束时间进行相应的折算得到)
3)计算目标时间段内所有的n个限功率时间段Txi的风机总限功率损失电量Q总损,计算公式如下:
Q总损=∑Q总损i=∑(Q理论i-Q实际i)=∑(∑Q理论i5/10min-∑Q实际i5/10min)
其中,i=1,2...n。
如图3所示,风场AGC限功率损失电量计算流程图:
步骤S8、采集风场升压站AGC设定值信号和各个风机有功功率数据,计算获得风场升压站AGC设定值与风机总有功功率的差值e,并按预设条件判断风场AGC限功率状态;具体为:根据采集的风场升压站AGC设定值信号和风机有功功率数据,首先需要计算风场升压站AGC设定值和风场总有功功率的差值e,计算频率为1min。
风场风机数为N风机,风场升压站AGC设定值为PAGC,风场总有功功率为P总,风场所有风机总装机容量为P额总,则:
P总=∑Pi,其中,Pi为第i个风机的有功功率,i=1,2...N风机;
P额总=∑Pi额定,其中,Pi额定为第i个风机的额定有功功率,i=1,2...N风 机;
e=PAGC-P总。
风场AGC限功率状态判定逻辑为以下两个条件同时成立:
e<h*P总;其中:h为系数,一般可取值0.02~0.05;
PAGC<P额总;
实际运行时,风场AGC限功率状态判定模块会根据上述状态满足/不满足的持续时间进行判定,假设t0时刻上述条件成立,分别延迟Δt和2Δt后进行第二次和第三次判定(Δt可配置为1min,一般不大于5min),若连续三次判定条件成立,则认定此次风场AGC限功率状态判定条件有效,记录开始时间t0。
此后每隔Δt进行一次判定,若在t1时刻判定上述条件不成立,则需要间隔Δt后进行二次判定,若连续两次判定条件均不成立,则判定风场AGC限功率状态结束,结束时间为t1。否则,风场AGC限功率状态持续不变。遍历风场所有风机,寻找风场AGC限功率时间段[t0,t1]内的处于限功率状态的风机,如存在风机限功率状态时间段包含在[t0,t1]中,则认定该时间段内的风机限功率状态为AGC限功率状态,不在[t0,t1]时间段内的风机限功率状态为风机自限电状态。统计目标时间段范围(如过去一天)的所有上述风场AGC限功率时间段,可以记录所有处于AGC限功率状态的风机状态时间段集合。
步骤S9、在风场AGC限功率状态下,并且在目标时间段内,针对步骤S6获得的限功率状态风机,获得AGC限功率状态下风机时间段集合,然后对该集合内的风机计算其AGC限功率损失电量,并求和获得总AGC限功率损失电量。具体为:
设统计目标时间段内所有风场AGC限功率时间段内出现AGC限功率状态的风机数为Nx,首先计算Nx个风机中的每个风机处于AGC限功率状态时间段内的损失电量QAGC损i,然后对所有Nx个风机的AGC限功率损失电量进行汇总求和,根据AGC限功率损失电量QAGC损i,在此基础上,统计目标时间段内所有风机总的限功率损失电量之和减去风场的总的AGC限功率损失电量即可得到风场所有风机总的自限电损失电量Q自限电总损,即:
Q自限电总损=Q总损-QAGC总损。
虽然本发明已以较佳实施例阐述如上,然其并非用以限定本发明。本发明所属技术领域中具有通常知识者,在不脱离本发明的精神和范围内,当可作各种的更动与润饰。因此,本发明的保护范围当视权利要求书所界定者为准。
Claims (9)
1.一种风机限电损失计算方法,其特征在于,首先分别针对风电场内各风机,执行如下步骤S1至步骤S6,获得风机的有功功率差值、以及限功率运行状态,然后执行步骤S7至S9,完成所有风机总的AGC自限电损失电量统计;
步骤S1、采集风机状态码、瞬时风速、有功功率、有功发电量、风机桨叶角对应历史时长的数据量,获得风机各指定类型数据对应历史时长的连续数据量;
步骤S2、对风机各指定类型数据对应历史时长的连续数据量进行预处理并更新,获得预设固定时间间隔的离散数据,进一步获得风机各指定类型数据对应历史时长的离散数据量;
步骤S3、根据离散数据的瞬时风速序列、预设时间窗,结合预设时间窗滑动时长,对该瞬时风速序列进行划分,获得各个时间窗对应的瞬时风速序列,并计算各个时间窗对应的瞬时风速序列的平均速度;进一步获得滑动时间窗平均风速序列;
步骤S4、基于风机各指定类型数据对应历史时长的连续数据量包括风机桨叶角、风机有功功率、风速的历史数据,按预设条件进行预处理,并绘制风机实际功率曲线,再根据风机对应的滑动时间窗平均风速序列,获得风机的理论有功功率值;根据风机实时数据、风机理论有功功率值,计算其有功功率差值;
步骤S5、根据风机状态码,将风机分类为:第一类风机、第二类风机,随后执行步骤S6;
步骤S6、根据风机状态码、实时有功功率、风机桨叶角,以及有功功率差值,分别按照预设条件判定第一类风机、以及第二类风机的限功率运行状态;
步骤S7、根据步骤S6获得的第一类风机、第二类风机的限功率状态、以及各个风机的理论有功功率值,计算各个风机的理论发电量,进一步获得所有风机总限功率损失电量;
步骤S8、采集风场升压站AGC设定值信号和各个风机有功功率数据,计算获得风场升压站AGC设定值与风机总有功功率的差值e,并按预设条件判断风场AGC限功率状态;
步骤S9、在风场AGC限功率状态下,并且在目标时间段内,针对步骤S6获得的限功率状态风机,获得AGC限功率状态下风机时间段集合,然后对该集合内的风机计算其AGC限功率损失电量,并求和获得总AGC限功率损失电量。
2.根据权利要求1所述的一种风机限电损失计算方法,其特征在于,步骤S1中,利用风电场SCADA系统采集风机的数据。
3.根据权利要求1所述的一种风机限电损失计算方法,其特征在于,步骤S2中,包括如下子步骤:
S201、对风机各指定类型数据对应历史时长的连续数据量存储与数据库,并对超过预设范围的数据进行剔除,并且进行数据更新;
S202、采用临近估计法将更新的数据进行字段统一处理,包括分别将时间戳、测点名称、测点值、测点所属风机名称、测点所属风场名称统一。
4.根据权利要求1所述的一种风机限电损失计算方法,其特征在于,所述步骤S3具体为:
根据预设时间窗ΔT,对瞬时风速序列{Vi,i=1,2...a…m....b...n},按如下公式计算时间窗瞬时风速序列的滑动平均风速:
V平均(T-1)=∑(V1,...,Vm)/m;
其中,Vm为T1+ΔT时刻的风速,a,b,m为时间值,n为总时长。
5.根据权利要求4所述的一种风机限电损失计算方法,其特征在于,步骤S4包括如下子步骤:
S401、按如下公式计算各个时间窗对应的瞬时风速序列的平均速度:
V平均(T-i)=∑(Va,...,Vb)/(b-a+1),
S402、按如下公式计算风机滑动平均风速序列:
{V平均(T-X),X=1,2..i...N},
其中,N为时间窗序号。
6.根据权利要求5所述的一种风机限电损失计算方法,其特征在于,步骤S5包括如下子步骤:
S501、将风机滑动平均风速序列{V平均(T-X),X=1,2..i...N}代入风机实际功率曲线P=f(V),计算风机的理论有功功率值P理论;
S502、根据风机实时有功功率值P实际,以及P理论,按如下公式计算风机有功功率差值:
ΔP=P理论-P实际。
7.根据权利要求6所述的一种风机限电损失计算方法,其特征在于,步骤S7包括如下子步骤:
S701、根据风机理论有功功率P理论,按如下公式计算各个风机的理论发电量:
Q理论=∫P理论dt;
S702、根据风机的限功率状态时间段集合{TXi,i=1,2...n},得到每个限功率时间段Txi的开始时间和结束时间,计算时间段TXi的风机实际发电量Q实际,如下式:
Q实际=μ*(Q末-Q初),
其中,Q初为时间段Txi开始时间的有功电量,Q末为时间段Txi结束时间的有功电量,μ为风机预设电量斜率参数;
S703、按如下公式计算n个限功率时间段Txi的风机总限功率损失电量Q总损:
Q总损=∑Q总损i=∑(Q理论i-Q实际i)=∑(∑Q理论i5/10min-∑Q实际i5/10min),
其中,i=1,2...n。
8.根据权利要求7所述的一种风机限电损失计算方法,其特征在于,步骤S8包括如下子步骤:
S801、计算风场所有风机总有功功率如下式:
P总=∑Pi,
其中,Pi为第i个风机的有功功率,i=1,2...N风机;风机数为N风机;
S802、按如下公式计算风场所有风机总装机容量:
P额总=∑Pi额定,其中,Pi额定为第i个风机的额定有功功率,i=1,2...N风机;
S803、如下公式计算风场总有功功率的差值:
e=PAGC-P总,其中,PAGC为风场升压站AGC设定值。
9.根据权利要求8所述的一种风机限电损失计算方法,其特征在于,步骤S9包括如下子步骤:
S901、根据AGC限功率损失电量QAGC损i,按如下公式计算风场的总的AGC限功率损失电量:
QAGC总损=∑QAGC损i,i=1,2...,Nx;
S902、按如下公式计算所有风机总的自限电损失电量:
Q自限电总损=Q总损-QAGC总损。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211522353.9A CN115905201B (zh) | 2022-11-30 | 2022-11-30 | 一种风机限电损失计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211522353.9A CN115905201B (zh) | 2022-11-30 | 2022-11-30 | 一种风机限电损失计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115905201A true CN115905201A (zh) | 2023-04-04 |
CN115905201B CN115905201B (zh) | 2023-11-17 |
Family
ID=86475924
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211522353.9A Active CN115905201B (zh) | 2022-11-30 | 2022-11-30 | 一种风机限电损失计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115905201B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108336761A (zh) * | 2018-04-03 | 2018-07-27 | 北京金风科创风电设备有限公司 | 风电场的功率控制方法、装置、系统和计算机设备 |
CN108631362A (zh) * | 2018-05-11 | 2018-10-09 | 北京中宜天成科技有限公司 | 一种风电场功率控制方法 |
CN112231361A (zh) * | 2020-10-28 | 2021-01-15 | 西安热工研究院有限公司 | 一种基于风机运行数据的风电项目发电量评估方法 |
KR20220028756A (ko) * | 2020-08-31 | 2022-03-08 | 인천대학교 산학협력단 | 풍력 터빈 발전기의 전력 편차를 최소화하는 풍력 발전기 시스템 |
CN114297251A (zh) * | 2021-11-08 | 2022-04-08 | 湖南大唐先一科技有限公司 | 一种风电机组损失电量计算方法 |
-
2022
- 2022-11-30 CN CN202211522353.9A patent/CN115905201B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108336761A (zh) * | 2018-04-03 | 2018-07-27 | 北京金风科创风电设备有限公司 | 风电场的功率控制方法、装置、系统和计算机设备 |
CN108631362A (zh) * | 2018-05-11 | 2018-10-09 | 北京中宜天成科技有限公司 | 一种风电场功率控制方法 |
KR20220028756A (ko) * | 2020-08-31 | 2022-03-08 | 인천대학교 산학협력단 | 풍력 터빈 발전기의 전력 편차를 최소화하는 풍력 발전기 시스템 |
CN112231361A (zh) * | 2020-10-28 | 2021-01-15 | 西安热工研究院有限公司 | 一种基于风机运行数据的风电项目发电量评估方法 |
CN114297251A (zh) * | 2021-11-08 | 2022-04-08 | 湖南大唐先一科技有限公司 | 一种风电机组损失电量计算方法 |
Non-Patent Citations (2)
Title |
---|
GULSHAN SHARMA等: "Automatic Generation Control (AGC) of Wind Power System: An Least Squares-Support Vector Machine (LS-SVM) Radial Basis Function(RBF) Kernel Approach", 《ELECTRIC POWER COMPONENTS AND SYSTEMS》, vol. 46, no. 14 * |
辛燕;李景才;任旦元;周民强;: "基于IEC信息分类的风电机组可利用率分析与提升方法", 山东大学学报(工学版), no. 06 * |
Also Published As
Publication number | Publication date |
---|---|
CN115905201B (zh) | 2023-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107885959B (zh) | 一种基于置信等效功率曲线带的风电建模及性能评估方法 | |
CN106570790B (zh) | 一种计及风速数据分段特性的风电场出力数据修复方法 | |
CN111522808B (zh) | 一种风电机组异常运行数据处理方法 | |
CN112966395B (zh) | 一种风电机组偏航系统静态误差的计算方法 | |
CN117150419B (zh) | 基于数据分析的储能电源运行环境检测系统 | |
CN115293372A (zh) | 基于多维度、多参量数值分析的光伏组串故障诊断方法 | |
CN113626767A (zh) | 一种风机限电识别方法、系统、设备和存储介质 | |
CN115456304A (zh) | 考虑台风影响的海上风电场可靠性指标计算方法及装置 | |
CN105591387A (zh) | 计及风电不确定性的基于安全域的电力系统安全监测方法 | |
CN110555784B (zh) | 基于门控循环神经网络风速修正的风电功率预测方法 | |
CN116089777A (zh) | 一种基于信息智能匹配的新能源智能结算方法及系统 | |
CN108471137B (zh) | 一种风功率预测中风速功率概率性映射方法 | |
CN115905201A (zh) | 一种风机限电损失计算方法 | |
CN111984906B (zh) | 考虑光伏与负荷时序相关性的概率潮流快速计算方法 | |
CN110363338B (zh) | 基于cos衰减性原子分解和混沌理论的风电功率超短期预测方法 | |
CN116706940A (zh) | 一种计及不同频率响应特性数据范围的惯量评估优化方法 | |
CN113991700B (zh) | 基于历史数据的一次调频优化决策方法及装置 | |
CN115577807A (zh) | 基于pmu数据驱动的轨道交通电力系统故障事件检测方法 | |
CN115170347A (zh) | 一种基于运行风功率曲线的阵列式风电场低效风电机组挖掘方法 | |
CN113554203B (zh) | 基于高维网格化和LightGBM的风电预测方法及装置 | |
CN111794921B (zh) | 一种基于迁移成分分析的陆上风电机组叶片结冰诊断方法 | |
CN113378102A (zh) | 一种用于短期负荷预测数据缺失预处理方法、介质及应用 | |
CN106251238B (zh) | 风电场建模序列离散化步长的选取及模型误差分析方法 | |
CN112561734A (zh) | 一种台区线损分析系统及其分析方法 | |
CN105242143B (zh) | 基于多时段精密测量单元数据状态估计不良数据修正方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |