CN109033027B - 一种高速阵风导致风电功率下爬坡事件的预测方法 - Google Patents

一种高速阵风导致风电功率下爬坡事件的预测方法 Download PDF

Info

Publication number
CN109033027B
CN109033027B CN201810894972.8A CN201810894972A CN109033027B CN 109033027 B CN109033027 B CN 109033027B CN 201810894972 A CN201810894972 A CN 201810894972A CN 109033027 B CN109033027 B CN 109033027B
Authority
CN
China
Prior art keywords
wind
speed
gust
wind power
operator
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
Application number
CN201810894972.8A
Other languages
English (en)
Other versions
CN109033027A (zh
Inventor
孙宏彬
潘欣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Baicheng Power Supply Co Of State Grid Jilin Electric Power Co ltd
Changchun Institute of Applied Chemistry of CAS
Original Assignee
Baicheng Power Supply Co Of State Grid Jilin Electric Power Co ltd
Changchun Institute of Applied Chemistry of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Baicheng Power Supply Co Of State Grid Jilin Electric Power Co ltd, Changchun Institute of Applied Chemistry of CAS filed Critical Baicheng Power Supply Co Of State Grid Jilin Electric Power Co ltd
Priority to CN201810894972.8A priority Critical patent/CN109033027B/zh
Publication of CN109033027A publication Critical patent/CN109033027A/zh
Application granted granted Critical
Publication of CN109033027B publication Critical patent/CN109033027B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations
    • G06F17/156—Correlation function computation including computation of convolution operations using a domain transform, e.g. Fourier transform, polynomial transform, number theoretic transform
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • Theoretical Computer Science (AREA)
  • Strategic Management (AREA)
  • Mathematical Physics (AREA)
  • Human Resources & Organizations (AREA)
  • Mathematical Analysis (AREA)
  • Data Mining & Analysis (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Health & Medical Sciences (AREA)
  • Computational Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Development Economics (AREA)
  • Quality & Reliability (AREA)
  • Public Health (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Computing Systems (AREA)
  • Operations Research (AREA)
  • Water Supply & Treatment (AREA)
  • Game Theory and Decision Science (AREA)
  • Algebra (AREA)
  • Primary Health Care (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Wind Motors (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供一种高速阵风导致风电功率下爬坡事件的预测方法，通过该方法在输入一个风电场的风速历史数据的同时，输入距离风电场最近的2个风力数据采集点的数据，通过这些数据构建高速阵风属性列表，进而进行一个风电场是否存在高速阵风导致风电功率下爬坡事件的情况进行预测。本发明不单纯依赖于一个风电场的数据，而且将距离风电场最近的2个风力数据采集点的数据也作为判断的辅助数据，解决了数据量高速阵风导致风电功率下爬坡事件量较少的问题；在风电功率下爬坡事件的预测方面具有较高的实用意义。

Description

一种高速阵风导致风电功率下爬坡事件的预测方法
技术领域：
本发明属于风电功率预测技术领域，尤其涉及一种高速阵风导致风电功率下爬坡事件的预测方法。
背景技术：
风电爬坡事件是指风电场或风电基地的风电功率输出在一个较短时间尺度内发生大幅度的变化。风电爬坡事件具有强波动性和强不确定性的特点，可造成系统中发/用电的严重不平衡，进而导致负荷损失。当风电场出现较高风速的阵风，使得风力涡轮机达到风速切出极限时，会导致典型的风电功率下爬坡事件，导致风电场有风但无输出功率的问题，对电网危害较大；准确的预测这种风电功率下爬坡事件，对于电网管理单位提前做好管理和应对准备，提高电网安全稳定运行、合理利用风电资源十分重要。
当前对于预测风电爬坡采用的手段主要是：收集风电爬坡事件，针对爬坡过程中的大气诸多因素构建样本,然后引入人工智能或统计学方法建立预测模型，随后利用预测模型对风电场的是否出现风功率爬坡可能进行预测。这些方法会遇到两个问题：（一）风电爬坡分为上爬坡和下爬坡事件类型，导致每种事件类型的引发原因又多种多样，当前方法试图通过一个模型或者算法将不同时间尺度、空间尺度和大气环境诸多因素整合在一起进行统一的预测，这种方式易于导致多种事件混淆，进而导致预测精度较低；（二）多数预测模型需要较多的训练爬坡事件来构建训练样本，风电场不可能选址在灾害频发、风力十分不稳定的地区，所以对于任意一个风电场，都难于保证获得足够多的样本来训练预测模型，可能导致产生的模型与特定的某几个爬坡事件数据耦合，预测能力较低。对于高速阵风导致风电功率下爬坡事件，目前较为实用的方法的较少。
因此需要提出一种方法，能够对一个风电场高速阵风导致风电功率下爬坡事件进行分析，并针对爬坡事件本身较少这一问题进行处理，使得风电场可以有效预测高速阵风导致风电功率下爬坡事件。
发明内容：
本发明提供一种高速阵风导致风电功率下爬坡事件的预测方法，解决了数据量高速阵风导致风电功率下爬坡事件量较少的问题；通过该方法在输入一个风电场的风速历史数据的同时，输入距离风电场最近的2个风力数据采集点的数据，通过这些数据构建高速阵风属性列表，进而进行一个风电场是否存在高速阵风导致风电功率下爬坡事件的情况进行预测。
本发明所述的一种高速阵风导致风电功率下爬坡事件的预测方法，包括以下步骤：
S1, 对于待预测风电场，输入其风电场风速历史数据H1，输入距离最近的2个风力数据采集点的风速历史数据D1,D2；分析起始时间点tstart，分析结束时间点tend；高速阵风导致风电功率下爬坡事件列表tlist;风电特征区间长度fdist;阵风预判区间长度zdist；计算采样间隔时间fgap；统计H1，获得H1的最大值H1Max，H1的最小值H1Min， H1的均值H1Mean，H1的标准差H1Stdv；
H1为待预测风电场的风速历史数据，为一个以时间为索引的数组，如H1[t]表示一个风电场在时间点t的风速值；
D1为风力数据采集点1的风速历史数据，风力数据采集点1为一个能够进行风力数据采集的地点，它距离待预测风电场直线距离最近，D1为一个以时间为索引的数组，如D1[t]表示一个在时间点t的风速值；
D2为风力数据采集点2的风速历史数据，风力数据采集点2为一个能够进行风力数据采集的地点，它距离待预测风电场直线距离比风力数据采集点1远，但是比其它风力数据采集点近，D2为一个以时间为索引的数组，如D2[t]表示一个在时间点t的风速值；
tlist为高速阵风导致风电功率下爬坡事件列表，tlist为一个列tlist={t1,t2,…,tn}，列表中的每一个元素对应一个时间点；
fdist为风电波动区间长度，该长度为描述和预测风电特征的时间长度值，其默认值为60，单位为分钟；
zdist为阵风预判区间长度，该长度为算法预判可能导致风电功率下爬坡事件的高速阵风出现的时间长度值，默认值为30，单位为分钟；
fgap为计算采样间隔时间，其默认值为10，单位为分钟；
S2，构建高速阵风特征算子WindOperator，算子可以抽取一个风速历史数据在一个时间点上的高速阵风特征；
S201，构建高速阵风特征算子WindOperator，该算子包含两个输入参数，第一个为一个以时间为索引的数组X1，第二个为特征抽取时间点tc；
S202，计算起始时间OStart=tc-fdist/2；计算结束时间OEnd=tc+fdist/2；
S203，风速数组WW=（从X1中以时间点OStart为起始，以OEnd为结束，以fgap为间隔取出的数据构成一个数组）；
S204，参考数量值NWW=(OEnd-OStart)/fgap+1；
S205，统计WW，获得WW的最大值WWMax,最小值WWMin,均值WWMean，标准差WWStdv；
S206，计算a1=(WWMean/H1Mean)，a2= (WW中高于WWMean的元素个数/NWW), a3=(WW的WWMax/H1Max)，a4=(WWStdv/H1Stdv)，a5=(WW中高于(WWMax-WWStdv/2)的元素个数/NWW),a6=(WW中高于(WWMean+WWStdv)的元素个数/NWW)；
S207, 构建高速阵风特征矢量v=(a1,a2,a3,a4,a5,a6)；
S208, 将高速阵风特征矢量v作为WindOperator的结果进行返回；
S3，基于高速阵风特征算子WindOperator，构建阵风判别算子DecisonOperator，算子可以对一个风速历史数据的一个时间点是否属于高速阵风进行判断；
S301, 构建阵风判别算子DecisonOperator,该算子包含两个输入参数，第一个为一个以时间为索引的数组DTX，第二个为特征抽取时间点dttc；
S302，高速阵风列表dtlist=空，高速阵风计数器dtcounter=1；
S303, dtnum=tlist的元素个数；
S304，dtnode=取出tlist的第dtcounter个元素,该元素为一个时间点数据；
S305，dtv=WindOperator(H1,dtnode),将H1和dtnode作为高速阵风特征算子WindOperator的输入，获得对应的高速阵风特征矢量；
S306，将dtv加入到高速阵风列表dtlist之中；
S307，dtcounter=dtcounter+1；
S308,如果dtcounter>dtnum则转到S307，否则转到S304；
S309，建立变量数组dty,dty包含dtnum个元素，dty每一个元素的初始值=1/dtnum；
S310, currentv=WindOperator(DTX,dttc),将DTX和dttc作为高速阵风特征算子的输入，获得对应的高速阵风特征矢量；
S310，将dty作为优化目标，优化求解以下公式的最小值，获得的结果存储到变量minvalue中；
S311，如果minvalue<(H1Stdv/2)则转到S312，否则转到S313；
S312，决策结果decision=TRUE，转到S314；
S313，决策结果decision=FALSE；
S314，将decision作DecisonOperator的结果返回；
S4，构建风速计算属性算子FeatureOperator，算子可以获取一个风速历史数据在一个时间点上的风速计算属性；
S401，构建风速计算属性算子FeatureOperator，该算子包含两个输入参数，第一个为一个以时间为索引的数组FX，第二个为特征抽取时间点ftc；
S402，计算属性起始时间FStart=ftc-fdist/2；计算属性结束时间FEnd=ftc+fdist/2；
S403，属性风速数组FWW=（从FX中以时间点OStart为起始，以OEnd为结束，以fgap为间隔取出的数据构成一个数组）；
S404，[ta1,ta2,ta3,ta4,ta5]=对FWW进行点数为5的傅里叶变换,并存储到5个变量之中；计算ta1=abs(ta1)/fdist,ta2=abs(ta2)/fdist,ta3=abs(ta3)/fdist，
ta4=abs(ta4)/fdist,ta5=abs(ta5)/fdist；
其中abs表示计算对应数的绝对值；
S405，对FWW进行统计获得FWW的均值为FWWMean，获得FWW最大值为FWWMax，获得FWW最小值FWMin；
S406，计算ta6=tanh（FWWMean/H1Max）, ta7=tanh(FWWMax/H1Max), ta8=tanh(FWMin/H1Max)；
S407，计算ta1=ta1,ta2=ta2-ta1,ta3=ta3-ta2,ta4=ta4-ta3,ta5=ta5-ta4；
S408，构建风速计算属性矢量Fwv=[ ta1,ta2,ta3,ta4,ta5,ta6,ta7,ta8]；
S409，将Fwv作为算子的结果返回；
S5,利用阵风判别算子DecisonOperator和风速计算属性算子FeatureOperator，从H1，D1，D2中抽取可用阵风特性，获得高速阵风属性列表highspeedlist；
S501，抽取时间currenttime=tstart；
S502，高速阵风属性列表highspeedlist=0；
S503，如果DecisonOperator(H1,currenttime)的结果为FALSE那么转到S505，,否则转到S504；
S504，计算FeatureOperator(H1,currenttime-zdist)获得的结果加入到highspeedlist之中；
S505，如果DecisonOperator(D1,currenttime)的结果为FALSE那么转到S507，否则转到S506；
S506，计算FeatureOperator(D1,currenttime-zdist)获得的结果加入到highspeedlist之中；
S507，如果DecisonOperator(D2,currenttime)的结果为FALSE那么转到S509,否则转到S508；
S508，计算FeatureOperator(D3,currenttime-zdist)获得的结果加入到highspeedlist之中；
S509，currenttime=currenttime+fgap；
S510，如果currenttime>tend则转到S511，否则转到S503；
S6，输入风电场新输入的风速数据数组FDCurrent，当前的时间点FDTime，利用风速计算属性算子FeatureOperator和高速阵风属性列表highspeedlist对当前时间点是否存在高速阵风导致风电功率下爬坡事件进行预测；
FDCurrent为待预测风电场的当前的风速数据数组，为一个以时间为索引的数组，如FDCurrent[t]表示一个风电场在时间点t的风速值；
S601，决策时间点DecisionPoint=FDTime-fdist/2；
S602, FDCV=计算FeatureOperator(FDCurrent,DecisionPoint)；
S603，FDCNum=highspeedlist中元素的个数；
S604，建立变量数组fdcy, fdcy包含FDCNum个元素，fdcy每一个元素的初始值=1/FDCNum；
S605，将fdcy作为优化目标，优化求解以下公式的最小值，获得的结果存储到变量fdcvalue中；
S606，如果fdcvalue<(H1Stdv/2)则转到S607，否则转到S608；
S607，预测当前风电场存在高速阵风导致风电功率下爬坡事件的风险，转到S609；
S608，预测当前风电场不存在高速阵风导致风电功率下爬坡事件的风险；
S609，预测过程结束。
本发明的积极效果在于：
提供一种高速阵风导致风电功率下爬坡事件的预测方法，通过该方法在输入一个风电场的风速历史数据的同时，输入距离风电场最近的2个风力数据采集点的数据，通过这些数据构建高速阵风属性列表，进而进行一个风电场是否存在高速阵风导致风电功率下爬坡事件的情况进行预测。本发明不单纯依赖于一个风电场的数据，而且将距离风电场最近的2个风力数据采集点的数据也作为判断的辅助数据，解决了数据量高速阵风导致风电功率下爬坡事件量较少的问题；
本发明针对于高速阵风导致风电功率下爬坡事件的方法，解决了多种事件混淆精度降低的问题；在风电功率下爬坡事件的预测方面具有较高的实用意义。
附图说明
图1为发明实施例2风电场风速历史数据H1图；
图2为发明实施例2风力数据采集点1的风速历史数据D1；
图3为发明实施例2风力数据采集点2的风速历史数据D2；
图4为发明实施例2当FDCurrent为该风电场在2016年4月1日至4月10日连续10天每天采样1440次的风速历史数据。
具体实施方式
通过以下实施例进一步举例描述本发明，并不以任何方式限制本发明，在不背离本发明的技术解决方案的前提下，对本发明所作的本领域普通技术人员容易实现的任何改动或改变都将落入本发明的权利要求范围之内。
实施例1
本发明一种高速阵风导致风电功率下爬坡事件的预测方法，包括以下步骤：
S1, 对于待预测风电场，输入其风电场风速历史数据H1，输入距离最近的2个风力数据采集点的风速历史数据D1,D2；分析起始时间点tstart，分析结束时间点tend；高速阵风导致风电功率下爬坡事件列表tlist;风电特征区间长度fdist;阵风预判区间长度zdist；计算采样间隔时间fgap；统计H1，获得H1的最大值H1Max，H1的最小值H1Min， H1的均值H1Mean，H1的标准差H1Stdv；
H1为待预测风电场的风速历史数据，为一个以时间为索引的数组，如H1[t]表示一个风电场在时间点t的风速值；
D1为风力数据采集点1的风速历史数据，风力数据采集点1为一个能够进行风力数据采集的地点，它距离待预测风电场直线距离最近，D1为一个以时间为索引的数组，如D1[t]表示一个在时间点t的风速值；
D2为风力数据采集点2的风速历史数据，风力数据采集点2为一个能够进行风力数据采集的地点，它距离待预测风电场直线距离比风力数据采集点1远，但是比其它风力数据采集点近，D2为一个以时间为索引的数组，如D2[t]表示一个在时间点t的风速值；
tlist为高速阵风导致风电功率下爬坡事件列表，tlist为一个列tlist={t1,t2,…,tn}，列表中的每一个元素对应一个时间点；
fdist为风电波动区间长度，该长度为描述和预测风电特征的时间长度值，其默认值为60，单位为分钟；
zdist为阵风预判区间长度，该长度为算法预判可能导致风电功率下爬坡事件的高速阵风出现的时间长度值，默认值为30，单位为分钟；
fgap为计算采样间隔时间，其默认值为10，单位为分钟；
S2，构建高速阵风特征算子WindOperator，算子可以抽取一个风速历史数据在一个时间点上的高速阵风特征；
S201，构建高速阵风特征算子WindOperator，该算子包含两个输入参数，第一个为一个以时间为索引的数组X1，第二个为特征抽取时间点tc；
S202，计算起始时间OStart=tc-fdist/2；计算结束时间OEnd=tc+fdist/2；
S203，风速数组WW=（从X1中以时间点OStart为起始，以OEnd为结束，以fgap为间隔取出的数据构成一个数组）；
S204，参考数量值NWW=(OEnd-OStart)/fgap+1；
S205，统计WW，获得WW的最大值WWMax,最小值WWMin,均值WWMean，标准差WWStdv；
S206，计算a1=(WWMean/H1Mean)，a2= (WW中高于WWMean的元素个数/NWW), a3=(WW的WWMax/H1Max)，a4=(WWStdv/H1Stdv)，a5=(WW中高于(WWMax-WWStdv/2)的元素个数/NWW),a6=(WW中高于(WWMean+WWStdv)的元素个数/NWW)；
S207, 构建高速阵风特征矢量v=(a1,a2,a3,a4,a5,a6)；
S208, 将高速阵风特征矢量v作为WindOperator的结果进行返回；
S3，基于高速阵风特征算子WindOperator，构建阵风判别算子DecisonOperator，算子可以对一个风速历史数据的一个时间点是否属于高速阵风进行判断；
S301, 构建阵风判别算子DecisonOperator,该算子包含两个输入参数，第一个为一个以时间为索引的数组DTX，第二个为特征抽取时间点dttc；
S302，高速阵风列表dtlist=空，高速阵风计数器dtcounter=1；
S303, dtnum=tlist的元素个数；
S304，dtnode=取出tlist的第dtcounter个元素,该元素为一个时间点数据；
S305，dtv=WindOperator(H1,dtnode),将H1和dtnode作为高速阵风特征算子WindOperator的输入，获得对应的高速阵风特征矢量；
S306，将dtv加入到高速阵风列表dtlist之中；
S307，dtcounter=dtcounter+1；
S308,如果dtcounter>dtnum则转到S307，否则转到S304；
S309，建立变量数组dty,dty包含dtnum个元素，dty每一个元素的初始值=1/dtnum；
S310, currentv=WindOperator(DTX,dttc),将DTX和dttc作为高速阵风特征算子的输入，获得对应的高速阵风特征矢量；
S310，将dty作为优化目标，优化求解以下公式的最小值，获得的结果存储到变量minvalue中；
S311，如果minvalue<(H1Stdv/2)则转到S312，否则转到S313；
S312，决策结果decision=TRUE，转到S314；
S313，决策结果decision=FALSE；
S314，将decision作DecisonOperator的结果返回；
S4，构建风速计算属性算子FeatureOperator，算子可以获取一个风速历史数据在一个时间点上的风速计算属性；
S401，构建风速计算属性算子FeatureOperator，该算子包含两个输入参数，第一个为一个以时间为索引的数组FX，第二个为特征抽取时间点ftc；
S402，计算属性起始时间FStart=ftc-fdist/2；计算属性结束时间FEnd=ftc+fdist/2；
S403，属性风速数组FWW=（从FX中以时间点OStart为起始，以OEnd为结束，以fgap为间隔取出的数据构成一个数组）；
S404，[ta1,ta2,ta3,ta4,ta5]=对FWW进行点数为5的傅里叶变换,并存储到5个变量之中；计算ta1=abs(ta1)/fdist,ta2=abs(ta2)/fdist,ta3=abs(ta3)/fdist，
ta4=abs(ta4)/fdist,ta5=abs(ta5)/fdist；
其中abs表示计算对应数的绝对值；
S405，对FWW进行统计获得FWW的均值为FWWMean，获得FWW最大值为FWWMax，获得FWW最小值FWMin；
S406，计算ta6=tanh（FWWMean/H1Max）, ta7=tanh(FWWMax/H1Max), ta8=tanh(FWMin/H1Max)；
S407，计算ta1=ta1,ta2=ta2-ta1,ta3=ta3-ta2,ta4=ta4-ta3,ta5=ta5-ta4；
S408，构建风速计算属性矢量Fwv=[ ta1,ta2,ta3,ta4,ta5,ta6,ta7,ta8]；
S409，将Fwv作为算子的结果返回；
S5,利用阵风判别算子DecisonOperator和风速计算属性算子FeatureOperator，从H1，D1，D2中抽取可用阵风特性，获得高速阵风属性列表highspeedlist；
S501，抽取时间currenttime=tstart；
S502，高速阵风属性列表highspeedlist=0；
S503，如果DecisonOperator(H1,currenttime)的结果为FALSE那么转到S505，,否则转到S504；
S504，计算FeatureOperator(H1,currenttime-zdist)获得的结果加入到highspeedlist之中；
S505，如果DecisonOperator(D1,currenttime)的结果为FALSE那么转到S507，否则转到S506；
S506，计算FeatureOperator(D1,currenttime-zdist)获得的结果加入到highspeedlist之中；
S507，如果DecisonOperator(D2,currenttime)的结果为FALSE那么转到S509,否则转到S508；
S508，计算FeatureOperator(D3,currenttime-zdist)获得的结果加入到highspeedlist之中；
S509，currenttime=currenttime+fgap；
S510，如果currenttime>tend则转到S511，否则转到S503；
S6，输入风电场新输入的风速数据数组FDCurrent，当前的时间点FDTime，利用风速计算属性算子FeatureOperator和高速阵风属性列表highspeedlist对当前时间点是否存在高速阵风导致风电功率下爬坡事件进行预测；
FDCurrent为待预测风电场的当前的风速数据数组，为一个以时间为索引的数组，如FDCurrent[t]表示一个风电场在时间点t的风速值；
S601，决策时间点DecisionPoint=FDTime-fdist/2；
S602，FDCV=计算FeatureOperator(FDCurrent,DecisionPoint)；
S603，FDCNum=highspeedlist中元素的个数；
S604，建立变量数组fdcy, fdcy包含FDCNum个元素，fdcy每一个元素的初始值=1/FDCNum；
S605，将fdcy作为优化目标，优化求解以下公式的最小值，获得的结果存储到变量fdcvalue中；
S606，如果fdcvalue<(H1Stdv/2)则转到S607，否则转到S608；
S607，预测当前风电场存在高速阵风导致风电功率下爬坡事件的风险，转到S609；
S608，预测当前风电场不存在高速阵风导致风电功率下爬坡事件的风险；
S609，预测过程结束。
实施例2
步骤S1，
输入风电场风速历史数据H1，该数据为某风电场2015年4月1日至4月10日连续10天每天采样1440次的风速历史数据（参见图1）；
输入风力数据采集点1的风速历史数据D1，该数据风力数据采集点1在2015年4月1日至4月10日连续10天每天采样1440次的风速历史数据（参见图2）；
风力数据采集点2的风速历史数据D2，风力数据采集点2的风速历史数据D1，该数据风力数据采集点1在2015年4月1日至4月10日连续10天每天采样1440次的风速历史数据（参见图3）；
分析起始时间点tstart=2015-4-1 0:0:0；
分析结束时间点tend=2015-4-1 23:59:59；
高速阵风导致风电功率下爬坡事件列表tlist={（2015-4-5 11:30:0）} 仅仅包含一次高速阵风导致的下爬坡事件；
fdist=60分钟；
zdist=30分钟；
fgap=10分钟；
在S1步骤对H1进行统计获得：
H1Max=30.9；
H1Min=4.3；
H1Mean=19.6；
H1Stdv=4.5；
S2，构建高速阵风特征算子WindOperator，算子可以抽取一个风速历史数据在一个时间点上的高速阵风特征；
S3，基于高速阵风特征算子WindOperator，构建阵风判别算子DecisonOperator，算子可以对一个风速历史数据的一个时间点是否属于高速阵风进行判断；
S4，构建风速计算属性算子FeatureOperator，算子可以获取一个风速历史数据在一个时间点上的风速计算属性；
S5，利用阵风判别算子DecisonOperator和风速计算属性算子FeatureOperator，从H1，D1，D2中抽取可用阵风特性，获得高速阵风属性列表highspeedlist；
该步骤计算结果highspeedlist包括以下内容：
{[0.31, 0.2, 0.11, 0.023, 0.17, 0.70, 0.83, 0.34]，
[0.27, 0.22, 0.17, 0.11, 0.19, 0.67, 0.87, 0.44]，
[0.21, 0.19, 0.09, 0.13, 0.17, 0.71, 0.81, 0.37]}；
S6，输入风电场新输入的风速数据数组FDCurrent，当前的时间点FDTime，利用风速计算属性算子FeatureOperator和高速阵风属性列表highspeedlist对当前时间点是否存在高速阵风导致风电功率下爬坡事件进行预测；
当FDCurrent为该风电场在2016年4月1日至4月10日连续10天每天采样1440次的风速历史数据（参见图4）；
测试的时间FDTime从2016-4-1 0:30:0为起始，以2016-4-1 23:30:0为终止，每次间隔10分钟，进行测试；测试的结果，在FDTime=2016-4-6 4:25:00时，FeatureOperator的输出结果为： [0.32, 0.19, 0.14, 0.14, 0.16, 0.75, 0.80, 0.40]；
步骤S6对应的输出结果为存在高速阵风导致风电功率下爬坡事件的风险，在FDTime为其他时间点时，步骤S6的输出结果为不存在高速阵风导致风电功率下爬坡事件的风险。
本发明的预测结果与其他方法的对比如下：
预测成功 漏掉的预测数 预测出现风险但未出现相关事件
本发明 19 0 2
支持向量机方法 8 11 34
决策树方法 10 9 23
结论：从表中可以看出，本发明可以成功预测所有的高速阵风导致风电功率下爬坡事件，相对支持向量机方法、决策树方法成功率较高，漏掉的预测数和预测出现风险但未出现相关事件的次数也明显优于传统方法。

Claims (1)

1.一种高速阵风导致风电功率下爬坡事件的预测方法，包括以下步骤：
S1，对于待预测风电场，风电场风速历史数据为H1，距离最近的2个风力数据采集点的风速历史数据为D1、D2；分析起始时间点tstart，分析结束时间点tend；高速阵风导致风电功率下爬坡事件列表tlist;风电特征区间长度fdist;阵风预判区间长度zdist；计算采样间隔时间fgap；统计H1，获得H1的最大值H1Max，H1的最小值H1Min， H1的均值H1Mean，H1的标准差H1Stdv；
H1为待预测风电场的风速历史数据，为一个以时间为索引的数组， H1[t]表示一个风电场在时间点t的风速值；
D1为风力数据采集点1的风速历史数据，风力数据采集点1为一个能够进行风力数据采集的地点，它距离待预测风电场直线距离最近，D1为风力数据采集点1的一个以时间为索引的数组，D1[t]表示一个在时间点t的风速值；
D2为风力数据采集点2的风速历史数据，风力数据采集点2为一个能够进行风力数据采集的地点，它距离待预测风电场直线距离比风力数据采集点1远，但是比其它风力数据采集点近，D2为风力数据采集点2的一个以时间为索引的数组，D2[t]表示一个在时间点t的风速值；
tlist为高速阵风导致风电功率下爬坡事件列表，tlist为一个列表，列表中每一个元素为一个时间点的数据；
fdist为风电波动区间长度，该长度为描述和预测风电特征的时间长度值，其默认值为60，单位为分钟；
zdist为阵风预判区间长度，该长度为算法预判可能导致风电功率下爬坡事件的高速阵风出现的时间长度值，默认值为30，单位为分钟；
fgap为计算采样间隔时间，其默认值为10，单位为分钟；
S2，构建高速阵风特征算子WindOperator，算子可以抽取一个风速历史数据在一个时间点上的高速阵风特征；
S201，构建高速阵风特征算子WindOperator，该算子包含两个输入参数，第一个为一个以时间为索引的数组X1，第二个为特征抽取时间点tc；
S202，计算起始时间OStart=tc-fdist/2；计算结束时间OEnd=tc+fdist/2；
S203，风速数组WW是从X1中以时间点OStart为起始，以OEnd为结束，以fgap为间隔取出的数据构成一个数组；
S204，参考数量值NWW=(OEnd-OStart)/fgap+1；
S205，统计WW，获得WW的最大值WWMax，最小值WWMin，均值WWMean，标准差WWStdv；
S206，计算a1=WWMean/H1Mean，a2=G1/NWW，a3=WWMax/H1Max，
a4=WWStdv/H1Stdv，a5=G2/NWW，a6=G3/NWW；
G1为WW中高于WWMean的元素个数；G2 为WW中高于(WWMax-WWStdv/2)的元素个数；G3为WW中高于(WWMean+WWStdv)的元素个数
S207，构建高速阵风特征矢量v=(a1，a2，a3，a4，a5，a6)；
S208，将高速阵风特征矢量v作为WindOperator的结果进行返回；
S3，基于高速阵风特征算子WindOperator，构建阵风判别算子DecisonOperator，算子可以对一个风速历史数据的一个时间点是否属于高速阵风进行判断；构建阵风判别算子DecisonOperator,该算子包含两个输入参数，第一个为一个以时间为索引的数组DTX，第二个为特征抽取时间点dttc；
S4，构建风速计算属性算子FeatureOperator，算子可以获取一个风速历史数据在一个时间点上的风速计算属性；
S401，构建风速计算属性算子FeatureOperator，该算子包含两个输入参数，第一个为一个以时间为索引的数组FX，第二个为特征抽取时间点ftc；
S402，计算属性起始时间FStart=ftc-fdist/2；计算属性结束时间FEnd=ftc+fdist/2；
S403，属性风速数组FWW为从FX中以时间点OStart为起始，以OEnd为结束，以fgap为间隔取出的数据构成一个数组；
S404，对FWW进行点数为5的傅里叶变换，并存储到5个变量ta1,ta2,ta3,ta4,ta5之中；计算ta1=abs(ta1)/fdist,ta2=abs(ta2)/fdist,ta3=abs(ta3)/fdist，ta4=abs(ta4)/fdist,ta5=abs(ta5)/fdist；
其中，abs表示计算对应数的绝对值；
S405，对FWW进行统计获得FWW的均值为FWWMean，获得FWW最大值为FWWMax，获得FWW最小值FWMin；
S406，计算ta6=tanh(FWWMean/H1Max)，ta7=tanh(FWWMax/H1Max)，
ta8=tanh(FWMin/H1Max)；
S407，计算ta1=ta1，ta2=ta2-ta1，ta3=ta3-ta2，ta4=ta4-ta3，ta5=ta5-ta4；
S408，构建风速计算属性矢量Fwv=[ ta1，ta2，ta3，ta4，ta5，ta6，ta7，ta8]；
S409，将Fwv作为算子的结果返回；
S5，利用阵风判别算子DecisonOperator和风速计算属性算子FeatureOperator，从H1，D1，D2中抽取可用阵风特性，获得高速阵风属性列表highspeedlist；
S501，抽取时间currenttime=tstart；
S502，高速阵风属性列表highspeedlist为空；
S503，如果DecisonOperator(H1，currenttime)的判断结果不属于高速阵风，转到S505，否则，转到S504；
S504，计算FeatureOperator(H1，C)获得的结果加入到highspeedlist之中；其中，C为currenttime-zdist
S505，如果DecisonOperator(D1，currenttime)的判断结果不属于高速阵风，转到S507，否则，转到S506；
S506，计算FeatureOperator(D1，C)获得的结果加入到highspeedlist之中；
S507，如果DecisonOperator(D2，currenttime)的判断结果不属于高速阵风，转到S509，否则，转到S508；
S508，计算FeatureOperator(D2，C)获得的结果加入到highspeedlist之中；
S509，currenttime=currenttime+fgap；
S510，如果currenttime>tend则结束，否则，转到S503；
S6，输入风电场新输入的风速数据数组FDCurrent，当前的时间点FDTime，利用风速计算属性算子FeatureOperator和高速阵风属性列表highspeedlist对当前时间点是否存在高速阵风导致风电功率下爬坡事件进行预测；
FDCurrent为待预测风电场的当前的风速数据数组，为一个以时间为索引的数组，FDCurrent[t]表示一个风电场在时间点t的风速值。
CN201810894972.8A 2018-08-08 2018-08-08 一种高速阵风导致风电功率下爬坡事件的预测方法 Active CN109033027B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810894972.8A CN109033027B (zh) 2018-08-08 2018-08-08 一种高速阵风导致风电功率下爬坡事件的预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810894972.8A CN109033027B (zh) 2018-08-08 2018-08-08 一种高速阵风导致风电功率下爬坡事件的预测方法

Publications (2)

Publication Number Publication Date
CN109033027A CN109033027A (zh) 2018-12-18
CN109033027B true CN109033027B (zh) 2022-07-29

Family

ID=64649975

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810894972.8A Active CN109033027B (zh) 2018-08-08 2018-08-08 一种高速阵风导致风电功率下爬坡事件的预测方法

Country Status (1)

Country Link
CN (1) CN109033027B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111582557A (zh) * 2020-04-20 2020-08-25 哈尔滨工业大学 基于变差函数的风电爬坡事件多级预警方法
CN111582751B (zh) * 2020-05-19 2022-06-14 国网吉林省电力有限公司 一种时间加权的购电风险预警方法
CN112270439B (zh) * 2020-10-28 2024-03-08 国能日新科技股份有限公司 超短期风电功率预测方法、装置、电子设备及存储介质
CN117196122B (zh) * 2023-11-02 2024-02-06 湖南赛能环测科技有限公司 一种基于风电功率爬坡时长的风电场调整方法及装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105068148A (zh) * 2015-07-14 2015-11-18 北京金风科创风电设备有限公司 风电场阵风预测方法和系统
CN106203693A (zh) * 2016-07-05 2016-12-07 华北电力大学 一种风电场输出功率爬坡事件预测的系统及方法
CN106933778A (zh) * 2017-01-22 2017-07-07 中国农业大学 一种基于爬坡事件特征识别的风电功率组合预测方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI476430B (zh) * 2010-08-23 2015-03-11 Inst Nuclear Energy Res Atomic Energy Council 具極端風速預測功能之風能預報方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105068148A (zh) * 2015-07-14 2015-11-18 北京金风科创风电设备有限公司 风电场阵风预测方法和系统
CN106203693A (zh) * 2016-07-05 2016-12-07 华北电力大学 一种风电场输出功率爬坡事件预测的系统及方法
CN106933778A (zh) * 2017-01-22 2017-07-07 中国农业大学 一种基于爬坡事件特征识别的风电功率组合预测方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
《Ramp Event Forecast Based Wind Power Ramp Control With Energy Storage System》;Yuzhong Gong 等;《IEEE Transactions on Power Systems》;20150630;全文 *
基于ASTFA-ELM的风电功率爬坡预测;黄麒元等;《电器与能效管理技术》;20170215(第03期);全文 *
基于小波变换和支持向量机的风电功率爬坡事件识别与预测;刘红柳等;《东北电力大学学报》;20161215(第06期);全文 *
基于气象背景选取邻近点的风电功率爬坡事件预测方法;欧阳庭辉等;《电网技术》;20151105(第11期);全文 *

Also Published As

Publication number Publication date
CN109033027A (zh) 2018-12-18

Similar Documents

Publication Publication Date Title
CN109033027B (zh) 一种高速阵风导致风电功率下爬坡事件的预测方法
US10312681B2 (en) Automatic network device electrical phase identification
CN102231216B (zh) 一种用于自助设备的现金管理方法和系统
CN108846555B (zh) 一种电力负荷大数据缺失值的高效精确填补法
US10598702B2 (en) Similarity detection of abnormal waveforms using PMU measurement
CN112712209B (zh) 水库入库流量预测方法、装置、计算机设备及存储介质
CN107437135B (zh) 一种新型储能选型方法
CN110580549A (zh) 一种考虑天气的区域短期能源功率预测方法及系统
EP3770423A1 (en) Wind turbine control method and device, controller, and control system
CN105846448B (zh) 一种基于随机矩阵理论的配电网无功补偿容量确定方法
CN109308589B (zh) 电网自动化数据质量监测方法、存储介质、终端设备和系统
CN113642827A (zh) 大坝监测数据的分析方法、装置、设备和存储介质
CN113988723A (zh) 一种基于用电数据异常分析的用户行为锁定方法和系统
CN117572837B (zh) 一种智慧电厂ai主动运维方法及系统
CN105469325B (zh) 一种确定火电机组负荷稳定状态的方法及系统
CN105743721A (zh) 数据上报方法、对上报数据进行处理的方法和装置
CN117277566B (zh) 基于大数据的电网数据分析电力调度系统及方法
CN106683001B (zh) 基于历史运行数据的火电机组辨识数据的选取方法
Zhang et al. A short-term user load forecasting with missing data
CN110020000B (zh) 判定异常风速数据的方法和装置
CN116226263A (zh) 一种可再生能源可视化智能管控制方法及系统
CN112596081B (zh) 一种基于关联性健康基线的健康状态监测的方法
CN110334308B (zh) 一种优化电源扫描参数测试中偏置条件的计算方法
CN113489005A (zh) 配电网潮流计算的配变负荷估算方法及系统
CN113554203A (zh) 基于高维网格化和LightGBM的风电预测方法及装置

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
TA01 Transfer of patent application right

Effective date of registration: 20220715

Address after: No. 395, KUANPING Road, Chaoyang District, Changchun City, Jilin Province

Applicant after: CHANGCHUN INSTITUTE OF TECHNOLOGY

Applicant after: BAICHENG POWER SUPPLY COMPANY OF STATE GRID JILIN ELECTRIC POWER Co.,Ltd.

Address before: No. 395, KUANPING Road, Chaoyang District, Changchun City, Jilin Province

Applicant before: CHANGCHUN INSTITUTE OF TECHNOLOGY

TA01 Transfer of patent application right