CN113762602A - 一种风电场短期风速预测方法 - Google Patents
一种风电场短期风速预测方法 Download PDFInfo
- Publication number
- CN113762602A CN113762602A CN202110929023.0A CN202110929023A CN113762602A CN 113762602 A CN113762602 A CN 113762602A CN 202110929023 A CN202110929023 A CN 202110929023A CN 113762602 A CN113762602 A CN 113762602A
- Authority
- CN
- China
- Prior art keywords
- individual
- wind speed
- sub
- individuals
- population
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 61
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 72
- 238000013528 artificial neural network Methods 0.000 claims abstract description 56
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 40
- 238000012360 testing method Methods 0.000 claims abstract description 26
- 238000005457 optimization Methods 0.000 claims abstract description 24
- 238000012549 training Methods 0.000 claims abstract description 24
- 238000005259 measurement Methods 0.000 claims abstract description 6
- 239000013598 vector Substances 0.000 claims description 27
- 230000008569 process Effects 0.000 claims description 14
- 230000007246 mechanism Effects 0.000 claims description 12
- 230000035772 mutation Effects 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 210000000349 chromosome Anatomy 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 6
- 108090000623 proteins and genes Proteins 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000006872 improvement Effects 0.000 claims description 5
- 210000002569 neuron Anatomy 0.000 claims description 5
- 238000013459 approach Methods 0.000 claims description 4
- 230000009286 beneficial effect Effects 0.000 claims description 4
- 230000001174 ascending effect Effects 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 125000004122 cyclic group Chemical group 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 claims description 3
- 230000004083 survival effect Effects 0.000 claims description 3
- 230000001360 synchronised effect Effects 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract 1
- 238000012163 sequencing technique Methods 0.000 description 3
- 238000010835 comparative analysis Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 101001095088 Homo sapiens Melanoma antigen preferentially expressed in tumors Proteins 0.000 description 1
- 102100037020 Melanoma antigen preferentially expressed in tumors Human genes 0.000 description 1
- 101710097421 WD repeat and HMG-box DNA-binding protein 1 Proteins 0.000 description 1
- 102100029469 WD repeat and HMG-box DNA-binding protein 1 Human genes 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009396 hybridization Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
Images
Classifications
-
- 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
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/086—Learning methods using evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- 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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Abstract
本发明公开了一种风电场短期风速预测方法,包括以下步骤:S1:获取风电场测风塔基本信息;S2:对于一定时段的风速数据,利用小波分解工具进行多层分解;S3:基于改进差分混合蛙跳‑神经网络算法建立测风数据预测模型,将经小波分解算法得到的各分量分为训练样本与测试样本,分别经所建模型进行模型训练与神经网络权值与阈值参数优化;S4:利用拟合样本数据,利用所建模型分层进行预测,后进行预测风速数据重构,与拟合数据进行对比,检验风速预测模型可靠性。本发明结构设计合理,采用人工神经网络算法进行风速数据的预测,更好的利用了人工神经网络算法鲁棒性、自学习、自适应特性,具有更好的风速数据预测精度。
Description
技术领域
本发明涉及风速预测技术领域,具体涉及一种风电场短期风速预测方法。
背景技术
随着近几年新能源行业的大力发展,风电装机容量占比不断提高,截止2019年,全国风力发电并网容量达到21005万千瓦,增长14%,占全部容量的10.45%。风力发电出力具有随机性强的特点,随着并网容量占比的不断提升,所需要的调峰调频容量也随之提高,近年来,电网实施的“两个细则”考核,提出了对风力发电的功率预测要求,实现对风电场出力的精准预测,风速的预测是关键。通过对短期风速的精准预测,进而得到准确的功率预测值,减少风电场应对“两个细则”要求的考核,提高电网整体运行水平和优化调度能力。
风速预测所用的算法主要由时间序列、支持向量机、小波神经网络、人工神经网络等,单一使用某种算法进行风速预测时,误差常达到25%~40%,预测精度低,因此通常使用组合算法进行风速预测,以提高预测精度,为此,我们提出一种风电场短期风速预测方法。
发明内容
本发明的目的在于克服现有技术中存在的上述问题,提供一种风电场短期风速预测方法,采用小波分解、混合蛙跳算法与人工神经网络进行组合,通过风速数据的小波多层分解,降低分辨率,改善风速量的波动性,利用混合蛙跳算法的遍历性强、不易陷入局部最优的特点,改进人工神经网络权值与阈值的寻优,通过对各层风速分解量进行神经网络权值与阈值优化,分别预测,以达到提高风速预测精度的目的。
为实现上述技术目的,达到上述技术效果,本发明是通过以下技术方案实现:
一种风电场短期风速预测方法,包括以下步骤:
S1:获取风电场测风塔基本信息;
S2:对于一定时段的风速数据,利用小波分解工具进行多层分解;
S3:基于改进差分混合蛙跳-神经网络算法建立测风数据预测模型,将经小波分解算法得到的各分量分为训练样本与测试样本,分别经所建模型进行模型训练与神经网络权值与阈值参数优化;
S4:利用拟合样本数据,利用所建模型分层进行预测,后进行预测风速数据重构,与拟合数据进行对比,检验风速预测模型可靠性。
优选地,本发明所述的一种风电场短期风速预测方法,所述步骤S1具体过程包括:
S101:收集风电场装机容量、场址、风机选型测风塔位置、实测各通道风速风向、场地气象环境数据。
优选地,本发明所述的一种风电场短期风速预测方法,所述步骤S2具体过程包括:
S201:利用小波分解算法进行风速数据的多层分解,得到高频分量与低频分量,选择合适的小波分解函数与分解层数P;
S202:基于选择的小波分解函数与确定的分解层数,计算得到低频分量Va,高频分量Vd1,Vd2...Vd(P-1)。
优选地,本发明所述的一种风电场短期风速预测方法,所述步骤S3具体过程包括:
S301:确定神经网络为三层拓扑结构,输入层、隐含层、输出层神经元数量分别为M、N、Q个,设置人工神经网络训练迭代次数T,学习率Lr,最小目标误差Er;
S302:确定染色体,包括输入隐含层权值w1、隐含阈值b1、隐含输出层权值w2、输出阈值b2,其基因数量分别为M*N、N、N*Q、Q,共同组成一条染色体,基因总数为D;
S303:确定混合蛙跳算法中蛙群个体数量P,子种群数量L,子种群迭代寻优次数S,全局混合迭代寻优次数G,每个子种群中蛙群个体数量为F=P/L,个体Xi=(x1,x2,...,xD),i∈P,种群X=(X1,X2,...,XP),个体Xi按照随机概率在[-R,R]间生成;
S304:计算每个个体的适应度值,按照升序或者降序统一排序,进行子种群个体排序,将第一只青蛙分配到第一个子种群中,第二个青蛙分配到第二个子种群中,以此类推,将所有个体划分为个子种群,形成子种群个体序列Xi=(xi1,xi2,...,xiF),i∈L,进行各子种群内部局部搜索更新,待满足子种群更新条件后,所有个体混合,按照适应度值进行排序分配,进行子种群更新寻优,周而复始,完成全局循环更新后输出最优个体;
S305:将最优个体中权值与阈值分别赋予人工神经网络结构,完成人工神经网络拓扑结构建模。
优选地,本发明所述的一种风电场短期风速预测方法,所述步骤S304具体过程包括:
S3041:将S2中经小波分解算法得到的高频与低频分量作为训练数据与测试数据,参与模型的训练进化;计算种群个体Xi的适应度值公式为:
式中:Vture为测风塔时段实测风速值,fnet(Xi)为人工神经网络同时段预测风速值,J∈R,为该时段内风速数据样本个数,f(Xi)越小表示个体Xi越优秀;
S3042:子种群内部搜索更新机制,采用差分进化算法(Difference Evolution,DE)进行改进,替换原有的子种群搜索更新机制,DE中特有的进化策略使得算法具有针对较优个体的记忆能力,促使较差个体迅速向最优个体靠拢,具有较强的全局收敛性和鲁棒性;
具体包括变异、交叉、选择三种操作。
优选地,本发明所述的一种风电场短期风速预测方法,所述步骤S3042具体包括以下内容:
S3042a:变异操作,通过把子种群中随机三个个体向量中两个向量的加权差向量加到第三个个体向量上得到变异个体,公式如下:
S3042b:由于差分算法在改进混合蛙跳算法中充当局部搜索的角色,因此,对上式提出改进措施,如下式所示:
式(1-3)中F为介于[0,2]的实数型常量因子,参数F的大小与差分矢量对变异个体的影响强弱有着紧密关系,当取值F较大时,生成的变异个体与当前个体差异较大,有助于保持种群的多样性和跳出局部最优,加快全局收敛;当F取值较小时,生成的变异个体与当前个体差异较小,使得DE围绕最优解精细化搜索,有利于局部搜索和加速收敛;基于F的这一特性,设计F动态变化策略,使得F随着混合蛙跳算法全局迭代次数的逐渐增加而逐渐减小,如下所示:
式中:Fmax为F最大取值,Fmin为F最小取值,GEN为混合蛙跳算法全局迭代次数,i=1,2,...,GEN;
S3042c:交叉操作,通过将目标个体与变异个体按照一定概率交换向量元素的方式产生试验个体,以此增加种群个体的多样性;
本文采用二项式交叉方法实现交叉操作,根据独立的Bernoulli实验确定变异个体Vi(t)和目标个体xi(t)中不同部位上元素进入到试验个体ui中,依据下面公式生成试验个体ui={ui1,ui2,...,uiD};
公式如下:
式中:CR∈[0,1],为杂交概率,控制变异个体和目标个体对试验个体的影响程度,rand1是[0,1]之间的随机数,rand2是[1,D]之间的随机整数,j=rand2保证uij中至少从变异个体中继承一个元素,确保试验个体与变异个体的差异性;
交叉概率CR的取值大小对DE的影响与参数F类似,当CR去取值较大时,交叉个体与变异个体差异较大,有助于算法未知空间的拓展;当CR取值较小时,交叉个体与变异个体差异较小,当算法迭代后期时,较小的CR利于保证具有良好向量因子的个体免遭过多的破坏,加快算法收敛;
采用如下方法使得CR随进化代数线性变化:
式中:CRmax为CR最大取值,CRmin为CR最小取值,其余同上;
S3042d:选择操作,基于适者生存的贪婪选择机制,将执行变异、交叉操作生产的个体,与子种群当前适应度值最差个体进行比较,若新个体适应度值优于当前最差个体,则进行替换,开始下一次子种群迭代更新;否则,重新进行变异与交叉操作,重复上述操作,进行适应度值比较,优选最佳个体进入下一次子种群迭代更新,若重复N次后仍未产生较佳个体,则随机选择变异、交叉个体进入下一次子种群迭代更新。
优选地,本发明所述的一种风电场短期风速预测方法,所述步骤S4具体内容包括:
S401:以低频分量Va为例,经模型进行神经网络初始权值与阈值寻优计算,完成模型构建,输入测试;
基于已完成权值、阈值参数优化,并完成建模的人工神经网络结构,利用拟合数据进行风速数据短期预测;
S402:将拟合样本数据输入风速预测模型,经过运算后得到输出样本数据;
S403:以此类推,分别完成剩余高频分量Vd1,Vd2...Vd(P-1)的输入与预测数据输出,将输出的各高频分量与低频分量重构成风速数据,将输出样本数据与同期实测风速数据进行对比分析,计算并评价拟合效果;采用平均绝对百分比误差、均方根误差指标进行计算,公式如下:
与现有技术相比,本发明的有益效果如下:
第一:采用人工神经网络算法进行风速数据的预测,更好的利用了人工神经网络算法鲁棒性、自学习、自适应特性,具有更好的风速数据预测精度。
第二:基于人工神经网络自学习易陷入局部最优的问题,引入混合蛙跳算法进行输入层、隐含层、输出层的权值与阈值寻优,较好的解决了人工神经网络易陷入局部最优的问题,提高了人工神经网络模型训练的效率与拟合精度。
第三:将差分算法引入混合蛙跳算法子种群迭代寻优过程,并对变异、交叉策略进行适当改进,特有的进化策略使得算法具有针对较优个体的记忆能力,促使较差个体迅速向最优个体靠拢,使其具有自适应、较强的全局收敛能力和鲁棒性的特点,使得新算法具有较优的收敛速度与精度。
当然,实施本发明的任一产品并不一定需要同时达到以上所述的所有优点。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明的基于混合蛙跳差分算法改进的人工神经网络算法流程图;
图2为本发明的DE算法选择操作流程图;
图3为本发明的原始风速数据;
图4为本发明的原始风速数据小波分解分量;
图5为本发明的预测数据与实测数据拟合与相对误差过程线。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
实施例1
请参阅图1-5所示,一种风电场短期风速预测方法,包括以下步骤:
S1:获取风电场测风塔基本信息
具体过程包括:S101:收集风电场装机容量、场址、风机选型测风塔位置、实测各通道风速风向、场地气象环境数据;
S2:对于一定时段的风速数据,利用小波分解工具进行多层分解;
具体过程包括:
S201:利用小波分解算法进行风速数据的多层分解,得到高频分量与低频分量,选择合适的小波分解函数与分解层数P;
S202:基于选择的小波分解函数与确定的分解层数,计算得到低频分量Va,高频分量Vd1,Vd2...Vd(P-1);
S3:基于改进差分混合蛙跳-神经网络算法建立测风数据预测模型,将经小波分解算法得到的各分量分为训练样本与测试样本,分别经所建模型进行模型训练与神经网络权值与阈值参数优化;
具体过程包括:
S301:确定神经网络为三层拓扑结构,输入层、隐含层、输出层神经元数量分别为M、N、Q个,设置人工神经网络训练迭代次数T,学习率Lr,最小目标误差Er;
S302:确定染色体,包括输入隐含层权值w1、隐含阈值b1、隐含输出层权值w2、输出阈值b2,其基因数量分别为M*N、N、N*Q、Q,共同组成一条染色体,基因总数为D;
S303:确定混合蛙跳算法中蛙群个体数量P,子种群数量L,子种群迭代寻优次数S,全局混合迭代寻优次数G,每个子种群中蛙群个体数量为F=P/L,个体Xi=(x1,x2,...,xD),i∈P,种群X=(X1,X2,...,XP),个体Xi按照随机概率在[-R,R]间生成;
S304:计算每个个体的适应度值,按照升序或者降序统一排序,进行子种群个体排序,将第一只青蛙分配到第一个子种群中,第二个青蛙分配到第二个子种群中,以此类推,将所有个体划分为个子种群,形成子种群个体序列Xi=(xi1,xi2,...,xiF),i∈L,进行各子种群内部局部搜索更新,待满足子种群更新条件后,所有个体混合,按照适应度值进行排序分配,进行子种群更新寻优,周而复始,完成全局循环更新后输出最优个体;
步骤S304具体过程包括:
S3041:将S2中经小波分解算法得到的高频与低频分量作为训练数据与测试数据,参与模型的训练进化;计算种群个体Xi的适应度值公式为:
式中:Vture为测风塔时段实测风速值,fnet(Xi)为人工神经网络同时段预测风速值,J∈R,为该时段内风速数据样本个数,f(Xi)越小表示个体Xi越优秀;
S3042:子种群内部搜索更新机制,采用差分进化算法(Difference Evolution,DE)进行改进,替换原有的子种群搜索更新机制,DE中特有的进化策略使得算法具有针对较优个体的记忆能力,促使较差个体迅速向最优个体靠拢,具有较强的全局收敛性和鲁棒性;
具体包括变异、交叉、选择三种操作:
具体包括:
S3042a:变异操作,通过把子种群中随机三个个体向量中两个向量的加权差向量加到第三个个体向量上得到变异个体,公式如下:
S3042b:由于差分算法在改进混合蛙跳算法中充当局部搜索的角色,因此,对上式提出改进措施,如下式所示:
式(1-3)中F为介于[0,2]的实数型常量因子,参数F的大小与差分矢量对变异个体的影响强弱有着紧密关系,当取值F较大时,生成的变异个体与当前个体差异较大,有助于保持种群的多样性和跳出局部最优,加快全局收敛;当F取值较小时,生成的变异个体与当前个体差异较小,使得DE围绕最优解精细化搜索,有利于局部搜索和加速收敛;基于F的这一特性,设计F动态变化策略,使得F随着混合蛙跳算法全局迭代次数的逐渐增加而逐渐减小,如下所示:
式中:Fmax为F最大取值,Fmin为F最小取值,GEN为混合蛙跳算法全局迭代次数,i=1,2,...,GEN;
S3042c:交叉操作,通过将目标个体与变异个体按照一定概率交换向量元素的方式产生试验个体,以此增加种群个体的多样性;
本文采用二项式交叉方法实现交叉操作,根据独立的Bernoulli实验确定变异个体Vi(t)和目标个体xi(t)中不同部位上元素进入到试验个体ui中,依据下面公式生成试验个体ui={ui1,ui2,...,uiD};
公式如下:
式中:CR∈[0,1],为杂交概率,控制变异个体和目标个体对试验个体的影响程度,rand1是[0,1]之间的随机数,rand2是[1,D]之间的随机整数,j=rand2保证uij中至少从变异个体中继承一个元素,确保试验个体与变异个体的差异性;
交叉概率CR的取值大小对DE的影响与参数F类似,当CR去取值较大时,交叉个体与变异个体差异较大,有助于算法未知空间的拓展;当CR取值较小时,交叉个体与变异个体差异较小,当算法迭代后期时,较小的CR利于保证具有良好向量因子的个体免遭过多的破坏,加快算法收敛;
采用如下方法使得CR随进化代数线性变化:
式中:CRmax为CR最大取值,CRmin为CR最小取值,其余同上;
S3042d:选择操作,基于适者生存的贪婪选择机制,将执行变异、交叉操作生产的个体,与子种群当前适应度值最差个体进行比较,若新个体适应度值优于当前最差个体,则进行替换,开始下一次子种群迭代更新;否则,重新进行变异与交叉操作,重复上述操作,进行适应度值比较,优选最佳个体进入下一次子种群迭代更新,若重复N次后仍未产生较佳个体,则随机选择变异、交叉个体进入下一次子种群迭代更新;
S305:将最优个体中权值与阈值分别赋予人工神经网络结构,完成人工神经网络拓扑结构建模;
S4:利用拟合样本数据,利用所建模型分层进行预测,后进行预测风速数据重构,与拟合数据进行对比,检验风速预测模型可靠性:
步骤S4具体过程包括:
S401:以低频分量Va为例,经模型进行神经网络初始权值与阈值寻优计算,完成模型构建,输入测试;
基于已完成权值、阈值参数优化,并完成建模的人工神经网络结构,利用拟合数据进行风速数据短期预测;
S402:将拟合样本数据输入风速预测模型,经过运算后得到输出样本数据;
S403:以此类推,分别完成剩余高频分量Vd1,Vd2...Vd(P-1)的输入与预测数据输出,将输出的各高频分量与低频分量重构成风速数据,将输出样本数据与同期实测风速数据进行对比分析,计算并评价拟合效果;采用平均绝对百分比误差、均方根误差指标进行计算,公式如下:
实施例2
请参阅图1-5所示,一种风电场短期风速预测方法,包括以下步骤:
S1:获取风电场测风塔基本信息,采用激光雷达进行测风,获取该测风塔2019-6-27,0:00至2019-6-27,23:50,90m高度处实测风速数据,采样频率为10min,数据点共380个,其中前350个数据点为训练样本,剩余30个数据点为拟合样本。
S2:选择小波分解算法db3,将样本数据分解为低频分量Va与高频分量Vd1,Vd2,Vd3。
S3:基于混合蛙跳算法,优化人工神经网络参数设置,引入差分算法改进混合蛙跳算法子种群寻优机制,建立测风数据插补模型;
具体的,步骤S3包括以下内容:
S301:确定人工神经网络隐含层神经元数量,经过反复试验,隐含层神经元数量为5,人工神经网络训练迭代次数T=100,学习率Lr=0.1,最小目标误差Er=0.0001;
S302:输入隐含层权值w1、隐含阈值b1、隐含输出层权值w2、输出阈值b2分别为1、5、5、1,蛙群个体向量包括16个元素,构成个体向量;
S303:设置混合蛙跳算法蛙群个体数量为P=100,子种群数量L=10,子种群迭代寻优次数S=20,全局混合迭代寻优次数G=100,每个子种群青蛙个体F=10,按照随机概率在(-1,1)之间生成;
S304:首先对全部种群个体进行初始适应度值计算按照适应度值由小到大进行排序,将第一只青蛙个体分配至第一个子种群,将第二只青蛙个体分配至第二个子种群,以此类推,共将100只青蛙个体分配至10个子种群,每个子种群中个体数量为10,进行子种群内部寻优迭代;
具体的,步骤S304包括以下内容:
a、采用差分进化算法代替原有的混合蛙跳算法子种群迭代寻优机制,包括变异、交叉和选择三种操作。设置差分缩放矢量因子F∈[0.6,0.8],杂交概率因子CR∈[0.3,0.45],F与CR随着全局迭代次数的增加逐渐减小;
b、在10个子种群内部分别按照变异、交叉、选择操作进行新个体产生与适应度值比较,优选优秀青蛙个体更新子种群;
c、当所有子种群内部寻优迭代全部完成后,混合所有青蛙个体,重新计算适应度值、排序与分配,开始第二轮全局迭代,当满足全局迭代次数时,结束程序;
S305:将最优青蛙个体中权值与阈值赋予人工神经网络结构,完成人工神经网络拓扑结构建模。
S4:完成人工神经网络结构参数设置,拟合样本数据进行短期风速数据预测,经模型运算后,将输出样本数据与同期测风塔实测数据进行对比分析;
具体的,S4步骤主要包括:
d、以低频分量Va为例,将输入数据进行归一化处理,基于步骤S3进行模型计算,生成初始蛙群,按照混合差分算法进行迭代寻优,完成子种群与全局寻优,输出最佳适应度值个体;
e、将最佳适应度值个体赋值给已构建好的三层神经网络,作为初始权值与阈值,完成神经网络训练;
f、基于已训练好的神经网络,将经归一化处理后的低频分量Va作为神经神经网络输入数据,得到预测数据Vd'1;
g、按照步骤S401~S403,完成剩余分解分量的神经网络预测,分别得到Va'1,Vd'2,Vd'3,将所有预测分量进行小波重构,得到预测风速数据;
h、利用式(1-7)、式(1-8)计算拟合精度指标εMAPE,εRMSE,评价算法模型短期风速预测效果。
实施例3
现已某风电场激光雷达测风数据为例,进行实施方式解析,按照上述具体实施方式中模型与参数进行设置,将381个数据样本分为训练样本(1~350)与拟合样本(351~381个),具体见图3,经小波分解算法db3分解为低频分量与3个高频分量,具体见图4,按照上述步骤分别输入模型进行神经网络初始权值与阈值寻优,待种群全局寻优结束,将最佳个体元素作为神经网络初始权值与阈值,进行神经网络训练,基于已训练好的神经网络完成分解分量的数据预测,将各分量预测数据进行小波重构,与实测拟合样本数据进行对比分析,具体见图5,计算绝对误差εMAPE为8.54%,均方差εRMSE为0.52,效果较好。
在本说明书的描述中,参考术语“一个实施例”、“示例”、“具体示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
以上公开的本发明优选实施例只是用于帮助阐述本发明。优选实施例并没有详尽叙述所有的细节,也不限制该发明仅为的具体实施方式。显然,根据本说明书的内容,可作很多的修改和变化。本说明书选取并具体描述这些实施例,是为了更好地解释本发明的原理和实际应用,从而使所属技术领域技术人员能很好地理解和利用本发明。本发明仅受权利要求书及其全部范围和等效物的限制。
Claims (7)
1.一种风电场短期风速预测方法,其特征在于:包括以下步骤:
S1:获取风电场测风塔基本信息;
S2:对于一定时段的风速数据,利用小波分解工具进行多层分解;
S3:基于改进差分混合蛙跳-神经网络算法建立测风数据预测模型,将经小波分解算法得到的各分量分为训练样本与测试样本,分别经所建模型进行模型训练与神经网络权值与阈值参数优化;
S4:利用拟合样本数据,利用所建模型分层进行预测,后进行预测风速数据重构,与拟合数据进行对比,检验风速预测模型可靠性。
2.根据权利要求1所述的一种风电场短期风速预测方法,其特征在于:所述步骤S1具体过程包括:
S101:收集风电场装机容量、场址、风机选型测风塔位置、实测各通道风速风向、场地气象环境数据。
3.根据权利要求1所述的一种风电场短期风速预测方法,其特征在于:所述步骤S2具体过程包括:
S201:利用小波分解算法进行风速数据的多层分解,得到高频分量与低频分量,选择合适的小波分解函数与分解层数P;
S202:基于选择的小波分解函数与确定的分解层数,计算得到低频分量Va,高频分量Vd1,Vd2...Vd(P-1)。
4.根据权利要求1所述的一种风电场短期风速预测方法,其特征在于:所述步骤S3具体过程包括:
S301:确定神经网络为三层拓扑结构,输入层、隐含层、输出层神经元数量分别为M、N、Q个,设置人工神经网络训练迭代次数T,学习率Lr,最小目标误差Er;
S302:确定染色体,包括输入隐含层权值w1、隐含阈值b1、隐含输出层权值w2、输出阈值b2,其基因数量分别为M*N、N、N*Q、Q,共同组成一条染色体,基因总数为D;
S303:确定混合蛙跳算法中蛙群个体数量P,子种群数量L,子种群迭代寻优次数S,全局混合迭代寻优次数G,每个子种群中蛙群个体数量为F=P/L,个体Xi=(x1,x2,...,xD),i∈P,种群X=(X1,X2,...,XP),个体Xi按照随机概率在[-R,R]间生成;
S304:计算每个个体的适应度值,按照升序或者降序统一排序,进行子种群个体排序,将第一只青蛙分配到第一个子种群中,第二个青蛙分配到第二个子种群中,以此类推,将所有个体划分为个子种群,形成子种群个体序列Xi=(xi1,xi2,...,xiF),i∈L,进行各子种群内部局部搜索更新,待满足子种群更新条件后,所有个体混合,按照适应度值进行排序分配,进行子种群更新寻优,周而复始,完成全局循环更新后输出最优个体;
S305:将最优个体中权值与阈值分别赋予人工神经网络结构,完成人工神经网络拓扑结构建模。
5.根据权利要求4所述的一种风电场短期风速预测方法,其特征在于:所述步骤S304具体过程包括:
S3041:将S2中经小波分解算法得到的高频与低频分量作为训练数据与测试数据,参与模型的训练进化;计算种群个体Xi的适应度值公式为:
式中:Vture为测风塔时段实测风速值,fnet(Xi)为人工神经网络同时段预测风速值,J∈R,为该时段内风速数据样本个数,f(Xi)越小表示个体Xi越优秀;
S3042:子种群内部搜索更新机制,采用差分进化算法(DifferenceEvolution,DE)进行改进,替换原有的子种群搜索更新机制,DE中特有的进化策略使得算法具有针对较优个体的记忆能力,促使较差个体迅速向最优个体靠拢,具有较强的全局收敛性和鲁棒性;
具体包括变异、交叉、选择三种操作。
6.根据权利要求5所述的一种风电场短期风速预测方法,其特征在于:所述步骤S3042具体包括:
S3042a:变异操作,通过把子种群中随机三个个体向量中两个向量的加权差向量加到第三个个体向量上得到变异个体,公式如下:
S3042b:由于差分算法在改进混合蛙跳算法中充当局部搜索的角色,因此,对上式提出改进措施,如下式所示:
式(1-3)中F为介于[0,2]的实数型常量因子,参数F的大小与差分矢量对变异个体的影响强弱有着紧密关系,当取值F较大时,生成的变异个体与当前个体差异较大,有助于保持种群的多样性和跳出局部最优,加快全局收敛;当F取值较小时,生成的变异个体与当前个体差异较小,使得DE围绕最优解精细化搜索,有利于局部搜索和加速收敛;基于F的这一特性,设计F动态变化策略,使得F随着混合蛙跳算法全局迭代次数的逐渐增加而逐渐减小,如下所示:
式中:Fmax为F最大取值,Fmin为F最小取值,GEN为混合蛙跳算法全局迭代次数,i=1,2,...,GEN;
S3042c:交叉操作,通过将目标个体与变异个体按照一定概率交换向量元素的方式产生试验个体,以此增加种群个体的多样性;
本文采用二项式交叉方法实现交叉操作,根据独立的Bernoulli实验确定变异个体Vi(t)和目标个体xi(t)中不同部位上元素进入到试验个体ui中,依据下面公式生成试验个体ui={ui1,ui2,...,uiD};
公式如下:
式中:CR∈[0,1],为杂交概率,控制变异个体和目标个体对试验个体的影响程度,rand1是[0,1]之间的随机数,rand2是[1,D]之间的随机整数,j=rand2保证uij中至少从变异个体中继承一个元素,确保试验个体与变异个体的差异性;
交叉概率CR的取值大小对DE的影响与参数F类似,当CR去取值较大时,交叉个体与变异个体差异较大,有助于算法未知空间的拓展;当CR取值较小时,交叉个体与变异个体差异较小,当算法迭代后期时,较小的CR利于保证具有良好向量因子的个体免遭过多的破坏,加快算法收敛;
采用如下方法使得CR随进化代数线性变化:
式中:CRmax为CR最大取值,CRmin为CR最小取值,其余同上;
S3042d:选择操作,基于适者生存的贪婪选择机制,将执行变异、交叉操作生产的个体,与子种群当前适应度值最差个体进行比较,若新个体适应度值优于当前最差个体,则进行替换,开始下一次子种群迭代更新;否则,重新进行变异与交叉操作,重复上述操作,进行适应度值比较,优选最佳个体进入下一次子种群迭代更新,若重复N次后仍未产生较佳个体,则随机选择变异、交叉个体进入下一次子种群迭代更新。
7.根据权利要求1所述的一种风电场短期风速预测方法,其特征在于:所述步骤S4具体过程包括:
S401:以低频分量Va为例,经模型进行神经网络初始权值与阈值寻优计算,完成模型构建,输入测试;
基于已完成权值、阈值参数优化,并完成建模的人工神经网络结构,利用拟合数据进行风速数据短期预测;
S402:将拟合样本数据输入风速预测模型,经过运算后得到输出样本数据;
S403:以此类推,分别完成剩余高频分量Vd1,Vd2...Vd(P-1)的输入与预测数据输出,将输出的各高频分量与低频分量重构成风速数据,将输出样本数据与同期实测风速数据进行对比分析,计算并评价拟合效果;采用平均绝对百分比误差、均方根误差指标进行计算,公式如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110929023.0A CN113762602A (zh) | 2021-08-13 | 2021-08-13 | 一种风电场短期风速预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110929023.0A CN113762602A (zh) | 2021-08-13 | 2021-08-13 | 一种风电场短期风速预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113762602A true CN113762602A (zh) | 2021-12-07 |
Family
ID=78789201
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110929023.0A Pending CN113762602A (zh) | 2021-08-13 | 2021-08-13 | 一种风电场短期风速预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113762602A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116405355A (zh) * | 2023-04-21 | 2023-07-07 | 重庆邮电大学 | 一种基于模式选择的信号解调方法 |
CN117277314A (zh) * | 2023-11-21 | 2023-12-22 | 深圳航天科创泛在电气有限公司 | 风电功率预测方法、装置、电子设备及可读存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108256683A (zh) * | 2018-01-16 | 2018-07-06 | 云南电网有限责任公司电力科学研究院 | 一种基于改进bp神经网络的风电场短期功率预测方法 |
CN109615142A (zh) * | 2018-12-18 | 2019-04-12 | 中国能源建设集团江苏省电力设计院有限公司 | 一种基于小波分析的风电场风速组合预测方法 |
-
2021
- 2021-08-13 CN CN202110929023.0A patent/CN113762602A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108256683A (zh) * | 2018-01-16 | 2018-07-06 | 云南电网有限责任公司电力科学研究院 | 一种基于改进bp神经网络的风电场短期功率预测方法 |
CN109615142A (zh) * | 2018-12-18 | 2019-04-12 | 中国能源建设集团江苏省电力设计院有限公司 | 一种基于小波分析的风电场风速组合预测方法 |
Non-Patent Citations (3)
Title |
---|
方国华等: "梯级水库生态调度多目标混合蛙跳差分算法研究", 水资源与水工程学报, vol. 28, no. 1, pages 69 - 73 * |
王红亚等: "混合蛙跳算法优化神经网络的同步电机转子故障检测", 微电机, vol. 50, no. 9, pages 22 - 26 * |
肖迁等: "基于改进的小波-BP神经网络的风速和风电功率预测", 电力系统保护与控制, vol. 42, no. 15, pages 80 - 86 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116405355A (zh) * | 2023-04-21 | 2023-07-07 | 重庆邮电大学 | 一种基于模式选择的信号解调方法 |
CN117277314A (zh) * | 2023-11-21 | 2023-12-22 | 深圳航天科创泛在电气有限公司 | 风电功率预测方法、装置、电子设备及可读存储介质 |
CN117277314B (zh) * | 2023-11-21 | 2024-03-08 | 深圳航天科创泛在电气有限公司 | 风电功率预测方法、装置、电子设备及可读存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110363344B (zh) | 基于miv-gp算法优化bp神经网络的概率积分参数预测方法 | |
CN112348271A (zh) | 基于vmd-ipso-gru的短期光伏功率预测方法 | |
CN111860982A (zh) | 一种基于vmd-fcm-gru的风电场短期风电功率预测方法 | |
CN104715096B (zh) | Bp神经网络预测二肽模型多极展开属性计算方法 | |
CN113762602A (zh) | 一种风电场短期风速预测方法 | |
CN107346459B (zh) | 一种基于遗传算法改进的多模式污染物集成预报方法 | |
CN107506861B (zh) | 一种短期风速预测方法、装置及系统 | |
CN102594909A (zh) | 基于共邻矩阵谱信息的多目标社区检测方法 | |
CN104408562A (zh) | 一种基于bp神经网络的光伏系统发电效率综合评估方法 | |
CN110807554A (zh) | 基于风电/光伏经典场景集的生成方法及系统 | |
CN108399470B (zh) | 一种基于多示例遗传神经网络的室内pm2.5预测方法 | |
CN112884236B (zh) | 一种基于vdm分解与lstm改进的短期负荷预测方法及系统 | |
CN113468817A (zh) | 一种基于igoa优化elm的超短期风电功率预测方法 | |
CN112069656A (zh) | 基于lssvm-nsgaii耐久性混凝土配合比多目标优化的方法 | |
Zhao et al. | Construction cost prediction based on genetic algorithm and BIM | |
CN104732067A (zh) | 一种面向流程对象的工业过程建模预测方法 | |
CN112765902A (zh) | 基于TentFWA-GD的RBF神经网络软测量建模方法及其应用 | |
CN112367130A (zh) | 一种无线网络信号传输强度计算方法及计算机存储介质 | |
Gui et al. | Wireless sensor network fault sensor recognition algorithm based on MM* diagnostic model | |
CN108717506A (zh) | 一种预测焦炭热态强度的方法 | |
CN109408896B (zh) | 一种污水厌氧处理产气量多元智能实时监控方法 | |
CN113553768B (zh) | 模型数据混合驱动的电网可靠性快速计算方法及装置 | |
Zhu | Research on adaptive combined wind speed prediction for each season based on improved gray relational analysis | |
CN116542382A (zh) | 基于混合优化算法的污水处理溶解氧浓度预测方法 | |
CN111310974A (zh) | 一种基于ga-elm的短期需水预测方法 |
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 |