廖圣瑄，馬曉明，韓中合，賈海坤

（1.華北電力大學(xué)能源動(dòng)力與機(jī)械工程學(xué)院，河北 保定 071003; 2.中國(guó)電力科學(xué)研究院有限公司，北京 100192）

0 引 言

近年來(lái)，隨著風(fēng)電行業(yè)的發(fā)展，風(fēng)電開(kāi)發(fā)的重心由陸地向海上轉(zhuǎn)移，風(fēng)電場(chǎng)建造規(guī)模不斷擴(kuò)大，風(fēng)電機(jī)組運(yùn)行過(guò)程中的安全性問(wèn)題也越發(fā)嚴(yán)重。風(fēng)電機(jī)組運(yùn)行過(guò)程中對(duì)載荷進(jìn)行評(píng)估可以保障機(jī)組運(yùn)行的安全性和經(jīng)濟(jì)性[1]。為此有必要對(duì)風(fēng)電機(jī)組主要部位的受載情況進(jìn)行分析，以保障機(jī)組正常運(yùn)行。風(fēng)電機(jī)組載荷主要來(lái)自于重力和空氣動(dòng)力等，載荷狀況會(huì)對(duì)風(fēng)電機(jī)組運(yùn)行造成一定影響。目前對(duì)風(fēng)電機(jī)組正常運(yùn)行的載荷情況進(jìn)行分析時(shí)，通常使用實(shí)驗(yàn)測(cè)試或者使用FAST和GH bladed等載荷計(jì)算軟件進(jìn)行仿真分析[2]。以GH bladed為例，該軟件計(jì)算載荷時(shí)主要通過(guò)實(shí)驗(yàn)日的環(huán)境參數(shù)和機(jī)組參數(shù)建立動(dòng)態(tài)數(shù)學(xué)模型，對(duì)風(fēng)機(jī)組正常運(yùn)行的載荷情況進(jìn)行仿真分析[3]，實(shí)際計(jì)算過(guò)程中存在計(jì)算量大、耗時(shí)長(zhǎng)的問(wèn)題，不適用于短期風(fēng)電機(jī)組載荷評(píng)估。周士棟等[4]通過(guò)建立BP神經(jīng)網(wǎng)絡(luò)對(duì)風(fēng)電機(jī)組部分載荷進(jìn)行預(yù)測(cè)；秦斌等[5]通過(guò)建立極限學(xué)習(xí)機(jī)對(duì)風(fēng)電機(jī)組葉根載荷進(jìn)行辨識(shí)。然而單一模型的神經(jīng)網(wǎng)絡(luò)在預(yù)測(cè)過(guò)程中往往存在陷入局部最優(yōu)解和誤差較大的問(wèn)題。鑒于實(shí)測(cè)和仿真過(guò)程中存在成本高、耗時(shí)長(zhǎng)的問(wèn)題，本文提出了基于組合算法的風(fēng)電機(jī)組輪轂載荷預(yù)測(cè)方法，在機(jī)組輪轂載荷預(yù)測(cè)時(shí)將單一神經(jīng)網(wǎng)絡(luò)進(jìn)行優(yōu)化，使用遺傳算法[6]優(yōu)化BP神經(jīng)網(wǎng)絡(luò)[7]，使用粒子群算法[8]優(yōu)化極限學(xué)習(xí)機(jī)[9]，再使用定權(quán)值輸入[10]和GRNN[11]融合的方式分別將上述兩個(gè)模型組合，建立新的模型進(jìn)行載荷預(yù)測(cè)。隨后以某3.4 MW海上風(fēng)電機(jī)組為例，按照IEC標(biāo)準(zhǔn)對(duì)機(jī)組載荷進(jìn)行實(shí)測(cè)，利用實(shí)測(cè)數(shù)據(jù)基于神經(jīng)網(wǎng)絡(luò)原理建立了GA-BP、PSO-ELM、BP-ELM-GRNN模型對(duì)輪轂載荷進(jìn)行預(yù)測(cè)。研究結(jié)果可為風(fēng)電機(jī)組輪轂載荷預(yù)測(cè)提供參考。

1 輪轂載荷預(yù)測(cè)組合模型

1.1 遺傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)建模

BP神經(jīng)網(wǎng)絡(luò)是一種多層前饋網(wǎng)絡(luò)，目前在數(shù)據(jù)預(yù)測(cè)領(lǐng)域運(yùn)用十分廣泛。訓(xùn)練BP神經(jīng)網(wǎng)絡(luò)時(shí)需不斷調(diào)整網(wǎng)絡(luò)的權(quán)值和閾值，其隱含層層數(shù)一般按經(jīng)驗(yàn)取值，隱含層神經(jīng)元個(gè)數(shù)計(jì)算公式如下：

式中：h——隱含層神經(jīng)元個(gè)數(shù)；

n——輸入層數(shù)；

m——輸出層數(shù)；

a——隨機(jī)取0～10內(nèi)的常數(shù)。

本文在訓(xùn)練網(wǎng)絡(luò)時(shí)選取最大訓(xùn)練次數(shù)為1 000次，目標(biāo)誤差設(shè)置為0.001，學(xué)習(xí)效率設(shè)置為0.1。在選用訓(xùn)練次數(shù)時(shí)，比較一般訓(xùn)練次數(shù)為500次和最大訓(xùn)練次數(shù)為1 000次時(shí)的相對(duì)誤差均值和決定系數(shù)如表1所示。由表可知，選用最大訓(xùn)練次數(shù)為1 000時(shí)精確度更高。

表1 不同訓(xùn)練次數(shù)下的誤差對(duì)比

遺傳算法（genetic algorithm，GA）按照自然選擇規(guī)律進(jìn)行，通過(guò)選擇、交叉、變異步驟逐次進(jìn)行篩選，以達(dá)到最優(yōu)適應(yīng)度。遺傳算法迭代計(jì)算過(guò)程中，具有更優(yōu)適應(yīng)度的種群取代原有種群，多次更新種群以達(dá)到目標(biāo)解。

遺傳算法對(duì)BP神經(jīng)網(wǎng)絡(luò)進(jìn)行優(yōu)化的過(guò)程中主要通過(guò)遺傳算法對(duì)BP神經(jīng)網(wǎng)絡(luò)的初始權(quán)值和閾值進(jìn)行修正。修正時(shí)將風(fēng)電機(jī)組的實(shí)測(cè)載荷與神經(jīng)網(wǎng)絡(luò)的預(yù)測(cè)載荷之間的誤差作為遺傳算法的適應(yīng)度，選取15次種群迭代后誤差最小值對(duì)應(yīng)的權(quán)值和閾值代入BP神經(jīng)網(wǎng)絡(luò)。X方向上彎矩在不同隱含層神經(jīng)元個(gè)數(shù)時(shí)對(duì)應(yīng)的適應(yīng)度曲線(xiàn)如圖1所示，圖1中各隱含層神經(jīng)元個(gè)數(shù)不同進(jìn)化代數(shù)對(duì)應(yīng)的適應(yīng)度為進(jìn)行5次試驗(yàn)計(jì)算后所取的平均值。由圖可以看出，在進(jìn)化代數(shù)為15時(shí)曲線(xiàn)逐漸平緩，每次迭代進(jìn)化需要6～10 min，且迭代時(shí)間隨著輸入數(shù)據(jù)增加而增加，選取過(guò)多的種群代數(shù)時(shí)增加迭代次數(shù)會(huì)導(dǎo)致計(jì)算冗長(zhǎng)。

圖1 不同隱含層節(jié)點(diǎn)適應(yīng)度曲線(xiàn)

相比于單一BP神經(jīng)網(wǎng)絡(luò)而言，利用GA優(yōu)化后的網(wǎng)絡(luò)在風(fēng)電機(jī)組輪轂載荷預(yù)測(cè)時(shí)精準(zhǔn)性更高。本文在參數(shù)設(shè)計(jì)時(shí)遺傳算法選擇進(jìn)化代數(shù)為15，種群規(guī)模為70，交叉概率為0.75，變異概率為0.05。GA-BP模型訓(xùn)練網(wǎng)絡(luò)時(shí)選取最大訓(xùn)練次數(shù)為1 000次，目標(biāo)誤差設(shè)置為0.001，學(xué)習(xí)效率設(shè)置為0.1。GA-BP預(yù)測(cè)模型計(jì)算流程見(jiàn)圖2。

圖2 GA-BP模型流程

1.2 粒子群算法優(yōu)化的極限學(xué)習(xí)機(jī)

極限學(xué)習(xí)機(jī)（extreme learning machines，ELM）相比于BP神經(jīng)網(wǎng)絡(luò)只有一層隱含層，在設(shè)定參數(shù)時(shí)隨機(jī)設(shè)定兩層之間的權(quán)值以及隱含層的閾值，與BP神經(jīng)網(wǎng)絡(luò)相比計(jì)算過(guò)程更快，只需要設(shè)置隱含層神經(jīng)元個(gè)數(shù)。極限學(xué)習(xí)機(jī)相比于支持向量機(jī)而言誤差相仿，但極限學(xué)習(xí)機(jī)建模更簡(jiǎn)單，計(jì)算速度更快。

粒子群算法（particle swarm optimization，PSO）最早起源于模擬自然界中鳥(niǎo)群捕食的過(guò)程。該算法通過(guò)迭代計(jì)算逐漸逼近最優(yōu)解，在迭代過(guò)程中，粒子不斷更新自身速度和位置，朝著自身最優(yōu)解和種群最優(yōu)解接近。粒子群算法進(jìn)行中粒子更新后的速度 Vid和位置 Xid如下：

其中c1、c2為學(xué)習(xí)因子；為慣性因子，通過(guò)調(diào)整的取值對(duì)全局網(wǎng)絡(luò)進(jìn)行修正；Pid為第i個(gè)變量的d維位置；Pgd為全局最優(yōu)解的d維位置；random（0,1）為0～1之間的任意數(shù)。

本文使用粒子群算法（PSO）優(yōu)化極限學(xué)習(xí)機(jī)（ELM）時(shí)，優(yōu)化其初始輸入的權(quán)值和閾值。PSO優(yōu)化ELM的過(guò)程中，先初始化粒子群算法的學(xué)習(xí)因子、慣性因子、進(jìn)化代數(shù)等參數(shù)，再將風(fēng)電機(jī)組輪轂的實(shí)際載荷與預(yù)測(cè)載荷之間的均方誤差作為粒子群算法的適應(yīng)度進(jìn)行迭代尋優(yōu)。X方向上彎矩不同學(xué)習(xí)因子對(duì)應(yīng)的適應(yīng)度曲線(xiàn)如圖3所示，其中各學(xué)習(xí)因子不同進(jìn)化代數(shù)所對(duì)應(yīng)的適應(yīng)度為進(jìn)行5次試驗(yàn)所取的平均值。

圖3 不同學(xué)習(xí)因子適應(yīng)度曲線(xiàn)

優(yōu)化過(guò)程中每個(gè)粒子需要將其當(dāng)前的適應(yīng)度與自身最優(yōu)適應(yīng)度、種群最優(yōu)適應(yīng)度進(jìn)行對(duì)比，不斷更新自身位置。本文計(jì)算時(shí)設(shè)置進(jìn)化代數(shù)為200，種群規(guī)模為80，初始隱層神經(jīng)元個(gè)數(shù)為40，取c1=c2=2.05，PSO迭代尋優(yōu)后將取得的各層之間的權(quán)值和初始隱層的閾值代入ELM中重塑網(wǎng)絡(luò)進(jìn)行建模預(yù)測(cè)。PSO-ELM預(yù)測(cè)模型計(jì)算流程見(jiàn)圖4。

圖4 PSO-ELM模型流程圖

1.3 定權(quán)值組合模型

在風(fēng)電機(jī)組載荷預(yù)測(cè)過(guò)程中，由于不同預(yù)測(cè)模型網(wǎng)絡(luò)結(jié)構(gòu)不同、預(yù)測(cè)方式有所差異，部分預(yù)測(cè)值與實(shí)測(cè)值誤差較大，對(duì)此提出建立定權(quán)值預(yù)測(cè)模型。定權(quán)值預(yù)測(cè)模型結(jié)合多個(gè)模型進(jìn)行組合預(yù)測(cè)，每個(gè)模型所占權(quán)重為固定值，模型之間相互組合以補(bǔ)償彼此誤差。組合模型公式如下：

式中：Td——定權(quán)組合模型預(yù)測(cè)結(jié)果，N；

TBPout——優(yōu)化后BP模型預(yù)測(cè)值，N；

TELMout——優(yōu)化后ELM模型預(yù)測(cè)值，N；

x——固定權(quán)值。

1.4 廣義回歸融合模型

在組合模型預(yù)測(cè)過(guò)程中還可以采用神經(jīng)網(wǎng)絡(luò)融合的方法，采用廣義回歸神經(jīng)網(wǎng)（general regression neural network，GRNN）對(duì)BP神經(jīng)網(wǎng)絡(luò)和極限學(xué)習(xí)機(jī)進(jìn)行融合。相比于傳統(tǒng)徑向基神經(jīng)網(wǎng)絡(luò)而言，GRNN同樣通過(guò)非線(xiàn)性回歸分析進(jìn)行建模，但其映射能力和學(xué)習(xí)效率都要優(yōu)于RBF。相比于傳統(tǒng)神經(jīng)網(wǎng)絡(luò)而言，GRNN沒(méi)有隱含層到輸出層這個(gè)環(huán)節(jié)，而是引入了模式層和求和層，簡(jiǎn)化結(jié)構(gòu)，加快了運(yùn)算速度。GRNN函數(shù)求和層與輸出值的維數(shù)有關(guān)，其模式層傳遞函數(shù)fi如下：

式中：σ——光滑因子；

Y——輸入值；

Yi——第i個(gè)神經(jīng)元的學(xué)習(xí)樣本。

本文利用廣義回歸神經(jīng)網(wǎng)絡(luò)對(duì)遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)和粒子群算法優(yōu)化后的極限學(xué)習(xí)機(jī)進(jìn)行融合時(shí)，首先將從風(fēng)電機(jī)組獲取的原始實(shí)驗(yàn)數(shù)據(jù)進(jìn)行處理，得到12 360組有效數(shù)據(jù)，將其分為三部分，分別有11 000組、1 260組、100組數(shù)據(jù)。將第一部分11 000組數(shù)據(jù)作為訓(xùn)練集1，其中的輪轂風(fēng)速等參數(shù)作為輸入集1，輪轂載荷作為輸出集1，分別代入GA-BP和PSO-ELM網(wǎng)絡(luò)中進(jìn)行訓(xùn)練。將第二部分1 260組數(shù)據(jù)和第三部分100組數(shù)據(jù)作為測(cè)試集1，代入訓(xùn)練完成后的BP網(wǎng)絡(luò)和ELM中進(jìn)行預(yù)測(cè)，分別保存其載荷預(yù)測(cè)值。將建立好的GA-BP和PSO-ELM模型中預(yù)測(cè)輸出的第二部分和第三部分載荷預(yù)測(cè)結(jié)果作為新的樣本，代入GRNN中作為輸入集2，仍然選擇其機(jī)組實(shí)測(cè)載荷作為輸出集2。對(duì)于新建立的GRNN模型，將第二部分各網(wǎng)絡(luò)預(yù)測(cè)結(jié)果和實(shí)測(cè)載荷作為訓(xùn)練集2，將第三部分各項(xiàng)數(shù)據(jù)作為測(cè)試集2進(jìn)行訓(xùn)練。訓(xùn)練完成后的BP-ELM-GRNN組合模型用于風(fēng)電機(jī)組輪轂載荷預(yù)測(cè)時(shí)誤差相對(duì)較低。

風(fēng)電機(jī)組通過(guò)實(shí)驗(yàn)獲得的輸入?yún)?shù)維數(shù)不定，將所得參數(shù)代入組合模型中，經(jīng)過(guò)GA-BP和PSOELM初步處理后的一次預(yù)測(cè)結(jié)果成為新的二維參數(shù)。將其代入GRNN模型中進(jìn)行預(yù)測(cè)可以降低輸入維數(shù)。GRNN用于融合神經(jīng)網(wǎng)絡(luò)過(guò)程見(jiàn)圖5。

圖5 GRNN融合神經(jīng)網(wǎng)絡(luò)過(guò)程

本文采用相對(duì)誤差均值Er1，絕對(duì)誤差百分比Er2和決定系數(shù)r2對(duì)模型的誤差進(jìn)行評(píng)測(cè)。其計(jì)算公式如下：

式中：n——預(yù)測(cè)輪轂載荷個(gè)數(shù)；

ei——各點(diǎn)預(yù)測(cè)誤差絕對(duì)值；

Ri——第i個(gè)點(diǎn)的實(shí)際測(cè)試值；

Ti——載荷預(yù)測(cè)模型的預(yù)測(cè)值；

Ti,f——模型擬合值；

2 工程概況

選取某3.4 MW海上風(fēng)電機(jī)組進(jìn)行載荷實(shí)測(cè)，其輪轂高度為95 m，額定風(fēng)速為10.2 m/s，葉片長(zhǎng)度為76.2 m，機(jī)組設(shè)計(jì)按照IEC61400-1-2005規(guī)范進(jìn)行。本文建立輪轂旋轉(zhuǎn)主軸坐標(biāo)系進(jìn)行預(yù)測(cè)，選取坐標(biāo)系原點(diǎn)位于輪轂中心，坐標(biāo)系隨著風(fēng)輪一起旋轉(zhuǎn)，坐標(biāo)系如圖6所示。圖6中XR，YR,ZR為坐標(biāo)軸，MiR為i方向上的彎矩，F(xiàn)iR為i方向上的受力。本文研究來(lái)自該風(fēng)電機(jī)組正常發(fā)電工況下連續(xù)運(yùn)行10 min所得到的風(fēng)機(jī)參數(shù)和載荷狀況，剔除異常值和缺失值后保留12 460組數(shù)據(jù)用于評(píng)估風(fēng)電機(jī)組輪轂短期受載情況。本次實(shí)驗(yàn)采用SCADA系統(tǒng)對(duì)風(fēng)電機(jī)組輪轂參數(shù)進(jìn)行收集，將輪轂轉(zhuǎn)速傳感器、偏航角傳感器等數(shù)據(jù)采集裝置布置在機(jī)組相應(yīng)部位。對(duì)于風(fēng)速等環(huán)境參數(shù)，采用氣象數(shù)據(jù)采集系統(tǒng)進(jìn)行收集。

圖6 風(fēng)電機(jī)組輪轂旋轉(zhuǎn)主軸坐標(biāo)系

為了得到較為精準(zhǔn)的載荷預(yù)估模型，本文將輪轂風(fēng)速、偏航偏差等基本參數(shù)現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)與輪轂載荷進(jìn)行關(guān)聯(lián)性分析后，選取合適的輸入集建立模型，將取自機(jī)組的參數(shù)設(shè)為變量X，輪轂載荷設(shè)為變量Y。目前常用的關(guān)聯(lián)性分析方法主要是皮爾遜相關(guān)系數(shù)法和斯皮爾曼相關(guān)系數(shù)法。本文計(jì)算了各項(xiàng)風(fēng)機(jī)參數(shù)與載荷間的皮爾遜相關(guān)系數(shù)，其計(jì)算公式如下：

其中E為數(shù)學(xué)期望。

在統(tǒng)計(jì)學(xué)中，一般通過(guò)相關(guān)系數(shù)來(lái)形容兩變量之間的相關(guān)程度，皮爾遜系數(shù)越高，相關(guān)性越強(qiáng)，當(dāng)皮爾遜系數(shù)在0.2以下時(shí)，一般認(rèn)為兩組數(shù)據(jù)間不相關(guān)。風(fēng)電機(jī)組主要參數(shù)相關(guān)于載荷的皮爾遜系數(shù)如表2所示。分析表中數(shù)據(jù)，選取關(guān)聯(lián)度較高的系數(shù)進(jìn)行建模，本文選取輪轂風(fēng)速、轉(zhuǎn)子轉(zhuǎn)速和平均俯仰角作為輸入集，MXR作為輸出集進(jìn)行建模。選取輪轂風(fēng)速和平均俯仰角作為輸入集，F(xiàn)XR作為輸出集建立模型。

表2 輪轂載荷的皮爾遜相關(guān)系數(shù)

3 結(jié)果分析

以輪轂載荷X方向上的彎矩MX建模預(yù)測(cè)時(shí)，其GA-BP、PSO-ELM、BP-ELM-GRNN預(yù)測(cè)模型的預(yù)測(cè)值與實(shí)際值對(duì)比如圖7～圖9所示。通過(guò)BPELM-GRNN建模預(yù)測(cè)輪轂載荷X方向上的受力FX與實(shí)際值對(duì)比如圖10所示。圖11～圖14為各類(lèi)模型用于載荷預(yù)測(cè)時(shí)預(yù)測(cè)結(jié)果與實(shí)測(cè)值之間的殘差圖。

圖7 GA-BP預(yù)測(cè)彎矩與實(shí)測(cè)值對(duì)比

圖8 PSO-ELM預(yù)測(cè)彎矩與實(shí)測(cè)值對(duì)比

圖9 GRNN預(yù)測(cè)彎矩與實(shí)測(cè)值對(duì)比

圖10 GRNN預(yù)測(cè)受力與實(shí)測(cè)值對(duì)比

圖11 風(fēng)電機(jī)組輪轂旋轉(zhuǎn)主軸坐標(biāo)系

圖12 PSO-ELM預(yù)測(cè)彎矩殘差圖

圖13 GRNN預(yù)測(cè)彎矩殘差圖

圖14 GRNN預(yù)測(cè)受力殘差圖

由圖7～圖9可以看出，風(fēng)電機(jī)組輪轂X方向上的彎矩在較大時(shí)，各學(xué)習(xí)器預(yù)測(cè)結(jié)果基本與真實(shí)值相同。對(duì)比各圖發(fā)現(xiàn)單一預(yù)測(cè)模型在預(yù)測(cè)過(guò)程中相比于優(yōu)化后的模型，由于其初始權(quán)值和閾值未達(dá)到最優(yōu)，預(yù)測(cè)誤差相對(duì)較大。同時(shí)單一預(yù)測(cè)模型在網(wǎng)絡(luò)建立過(guò)程中可能陷入局部最優(yōu)解，影響實(shí)驗(yàn)結(jié)果。由圖11～圖14的殘差圖可以看出，使用各類(lèi)模型進(jìn)行預(yù)測(cè)時(shí)的擬合效果較好，用于機(jī)組載荷預(yù)測(cè)的精確度較高。由圖7～圖14可以看出，采用GRNN對(duì)BP、ELM組合可解決單一預(yù)測(cè)機(jī)在預(yù)測(cè)過(guò)程中出現(xiàn)的問(wèn)題，達(dá)到很好的擬合效果。BP、ELM和GRNN都屬于強(qiáng)學(xué)習(xí)器，其預(yù)測(cè)性能整體而言準(zhǔn)確性較高。各類(lèi)學(xué)習(xí)器在風(fēng)電機(jī)組輪轂載荷預(yù)測(cè)過(guò)程中誤差都比較小，在風(fēng)電機(jī)組早期載荷評(píng)價(jià)過(guò)程中使用可以達(dá)到目標(biāo)精度。通過(guò)神經(jīng)網(wǎng)絡(luò)融合的方式在一定程度上改善了單個(gè)學(xué)習(xí)器的局限性，提高了預(yù)測(cè)的準(zhǔn)確性。各類(lèi)學(xué)習(xí)器預(yù)測(cè)誤差和真實(shí)度如表3所示。

表3 不同算法預(yù)測(cè)結(jié)果

由表可知，BP-ELM-GRNN組合模型訓(xùn)練風(fēng)電機(jī)組輪轂載荷效果相比于單個(gè)網(wǎng)絡(luò)而言，誤差有所減小。采用GRNN對(duì)遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)和粒子群算法優(yōu)化后的極限學(xué)習(xí)機(jī)進(jìn)行融合后預(yù)測(cè)精度得到提高，輪轂彎矩MX的相對(duì)誤差均值降低到0.044 4，決定系數(shù)r2提高到0.879 6，輪轂受力FX的相對(duì)誤差均值降低為0.048 3，決定系數(shù)r2達(dá)到0.801 7。定權(quán)值組合模型模擬效果與GRNN組合算法預(yù)測(cè)結(jié)果相差不大，同樣起到優(yōu)化單一預(yù)測(cè)模型的效果。由此可以看出，使用神經(jīng)網(wǎng)絡(luò)組合模型可在一定程度上提高風(fēng)電機(jī)組輪轂載荷狀況預(yù)測(cè)精度。

4 結(jié)束語(yǔ)

本文以某3.4 MW海上風(fēng)電機(jī)組為例，研究了基于GA-BP、PSO-ELM、BP-ELM-GRNN的風(fēng)電機(jī)組輪轂載荷預(yù)測(cè)方法。建立多個(gè)預(yù)測(cè)模型進(jìn)行對(duì)比，得到以下結(jié)論：

1）本文從實(shí)測(cè)所得的機(jī)組參數(shù)中利用相關(guān)性分析找到主要輸入?yún)?shù)，提出用組合算法對(duì)風(fēng)電機(jī)組輪轂載荷進(jìn)行建模預(yù)測(cè)，使用GRNN組合算法預(yù)測(cè)時(shí)輪轂X方向上彎矩MX的相對(duì)誤差均值為0.044 4，決定系數(shù)r2為0.879 6，其X方向上的受力FX的相對(duì)誤差均值為0.048 3，決定系數(shù)r2為0.801 7。

2）本文分別采用遺傳算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)的初始權(quán)值和閾值，采用粒子群算法優(yōu)化極限學(xué)習(xí)機(jī)，使得優(yōu)化后的模型預(yù)測(cè)誤差相比于原模型有所減小，提高了預(yù)測(cè)的準(zhǔn)確性。

3）與傳統(tǒng)的神經(jīng)網(wǎng)絡(luò)相比，本文提出的組合模型預(yù)測(cè)方法有效結(jié)合了各個(gè)網(wǎng)絡(luò)的優(yōu)勢(shì)，彌補(bǔ)了各單一模型在預(yù)測(cè)時(shí)可能存在的不足。采用組合算法預(yù)測(cè)的方式提高了預(yù)測(cè)時(shí)的精度，保障了機(jī)組運(yùn)行的安全性和經(jīng)濟(jì)性。

本文利用BP-ELM-GRNN組合算法進(jìn)行機(jī)組輪轂載荷預(yù)測(cè)還存在以下預(yù)期：1）對(duì)于取自NWP的數(shù)據(jù)可以增加輸入?yún)?shù)，如地形地貌參數(shù)、大氣參數(shù)等；2）考慮針對(duì)尾流進(jìn)行分析建模；3）研究從單個(gè)機(jī)組到整個(gè)風(fēng)電場(chǎng)同一時(shí)間段的載荷狀況。