收稿日期：2022-07-25

基金項(xiàng)目：國(guó)家自然科學(xué)基金（62366012；61563012）；廣西自然科學(xué)基金（2021GXNSFAA220074）；廣西嵌入式技術(shù)與智能系統(tǒng)重點(diǎn)實(shí)驗(yàn)室基金

（2020-1-3）

通信作者：董明剛（1977—），男，博士、教授，主要從事智能計(jì)算、人工智能方面的研究。d2015mg@qq.com

DOI：10.19912/j.0254-0096.tynxb.2022-1106 文章編號(hào)：0254-0096（2023）11-0279-08

摘 要：針對(duì)現(xiàn)有的非均勻風(fēng)電場(chǎng)發(fā)電功率預(yù)測(cè)方法準(zhǔn)確度不高、可解釋性差等問(wèn)題，提出一種改進(jìn)的物理引導(dǎo)的圖神經(jīng)網(wǎng)絡(luò)（IPGNN）模型來(lái)預(yù)測(cè)非均勻風(fēng)電場(chǎng)的發(fā)電功率輸出。首先，構(gòu)建一個(gè)基于三維高斯尾流模型的物理信息基函數(shù)，可更準(zhǔn)確地反映非均勻風(fēng)電場(chǎng)中風(fēng)力機(jī)之間尾流相互作用的關(guān)系；其次，設(shè)計(jì)一組基于消息傳遞框架的圖神經(jīng)網(wǎng)絡(luò)更新策略，該策略在圖神經(jīng)網(wǎng)絡(luò)邊緣更新中結(jié)合注意力機(jī)制將物理信息基函數(shù)作為權(quán)重更新函數(shù)，可增強(qiáng)模型的可解釋性。在不同數(shù)量風(fēng)力機(jī)下的非均勻風(fēng)電場(chǎng)仿真實(shí)驗(yàn)表明，相比于基于一維尾流模型的PGNN模型和典型的數(shù)據(jù)驅(qū)動(dòng)模型，IPGNN模型均可獲得較好的預(yù)測(cè)效果，其中對(duì)具有20臺(tái)風(fēng)力機(jī)的非均勻風(fēng)電場(chǎng)，發(fā)電功率預(yù)測(cè)平均絕對(duì)誤差為3.92%，可認(rèn)為是一種有效的非均勻風(fēng)電場(chǎng)發(fā)電功率預(yù)測(cè)方法。

關(guān)鍵詞：深度學(xué)習(xí)；圖神經(jīng)網(wǎng)絡(luò)；風(fēng)電功率；風(fēng)電場(chǎng)；數(shù)值模擬

中圖分類號(hào)：TM614""""""""""""""""""""""""""""" 文獻(xiàn)標(biāo)志碼：A

0 引 言

在所有可持續(xù)能源中，風(fēng)能被認(rèn)為是最有前途的一種［1］。據(jù)全球風(fēng)能理事會(huì)（GWEC）2022年全球風(fēng)能報(bào)告指出，截至2021年底，全球風(fēng)電累計(jì)裝機(jī)已達(dá)837 GW，面對(duì)新冠肺炎疫情，風(fēng)電行業(yè)顯示出強(qiáng)勁韌性［2］。在風(fēng)電場(chǎng)設(shè)計(jì)中最需要關(guān)注的是風(fēng)力機(jī)之間產(chǎn)生的尾流效應(yīng)，上游風(fēng)力機(jī)從風(fēng)中獲取能量的同時(shí)在其下游形成尾流區(qū)，導(dǎo)致位于尾流區(qū)的下游風(fēng)力機(jī)的風(fēng)能減少。與無(wú)尾流效應(yīng)相比，尾流造成的損失通常能使發(fā)電量減少10%～20%［3］。因此，對(duì)風(fēng)電場(chǎng)預(yù)測(cè)其有效發(fā)電量并根據(jù)計(jì)算結(jié)果進(jìn)行合理優(yōu)化布局是風(fēng)電場(chǎng)規(guī)劃設(shè)計(jì)中的重要一環(huán)。

如今現(xiàn)實(shí)中所構(gòu)建的風(fēng)電場(chǎng)大多只由一種類型風(fēng)力機(jī)組成，只有小部分研究考慮了多種風(fēng)力機(jī)類型的選擇，這種風(fēng)電場(chǎng)被稱為非均勻風(fēng)電場(chǎng)。文獻(xiàn)［4］將不同輪轂高度風(fēng)力機(jī)組合起來(lái)設(shè)計(jì)了一種風(fēng)電場(chǎng)布局，發(fā)現(xiàn)與統(tǒng)一輪轂高度的風(fēng)電場(chǎng)相比，其發(fā)電功率輸出增加了13.5%，總成本卻降低了0.4%。非均勻風(fēng)電場(chǎng)也已被證實(shí)在地形上可受到更小限制，能充分利用不同高度的風(fēng)力資源，這使得非均勻風(fēng)電場(chǎng)將成為未來(lái)風(fēng)電場(chǎng)設(shè)計(jì)的首選。目前，為了模擬風(fēng)電場(chǎng)中風(fēng)力機(jī)之間的尾流效應(yīng)并研究其對(duì)發(fā)電功率的影響，研究方法大致分為兩大類：1）采用基于網(wǎng)格的模型，并使用基于搜索的優(yōu)化算法（遺傳算法、粒子群算法等）進(jìn)行求解；2）應(yīng)用數(shù)據(jù)驅(qū)動(dòng)模型（統(tǒng)計(jì)機(jī)器學(xué)習(xí)模型）來(lái)擬合風(fēng)電場(chǎng)風(fēng)流條件與發(fā)電功率之間的關(guān)系，建立相應(yīng)輸入與輸出模型。基于網(wǎng)格的模型是最早用于風(fēng)電場(chǎng)優(yōu)化的模型，由Mosetti等［5］在1994年提出并使用遺傳算法求解，隨后眾多研究學(xué)者致力于該領(lǐng)域的研究并提出粒子群算法、差分進(jìn)化算法等進(jìn)行求解［6-7］。使用基于搜索的求解方法有效地將其與工程尾流模型結(jié)合起來(lái)，但算法本身特性難以保證其結(jié)果的精確性，往往又需要大量數(shù)值校準(zhǔn)，通常十分耗時(shí)。數(shù)據(jù)驅(qū)動(dòng)的方法用于預(yù)測(cè)目標(biāo)風(fēng)電場(chǎng)在特定風(fēng)況與位置下風(fēng)力機(jī)發(fā)電量［8-9］，旨在通過(guò)基于各種統(tǒng)計(jì)模型構(gòu)建輸入和輸出之間的關(guān)系來(lái)描述觀察到的響應(yīng)，即機(jī)組布局和發(fā)電量之間的關(guān)系。這種方法可充分利用高性能圖形處理器（graphics processing unit，GPU）并行計(jì)算，一般可在較短時(shí)間內(nèi)完成模型的訓(xùn)練，但對(duì)數(shù)據(jù)的依賴較大，只適用于實(shí)際能獲取大量數(shù)據(jù)的風(fēng)電場(chǎng)，且模型的泛化能力和可解釋性并不理想。

為了解決上述方法所存在的計(jì)算成本昂貴、泛化能力差、可解釋性不足等問(wèn)題，Park等［10］提出一種物理引導(dǎo)的圖神經(jīng)網(wǎng)絡(luò)（physics-induced graph neural network，PGNN）解決均勻風(fēng)電場(chǎng)發(fā)電功率預(yù)測(cè)問(wèn)題，相比于傳統(tǒng)數(shù)據(jù)驅(qū)動(dòng)方法取得了更好的效果，但該方法僅適用于均勻風(fēng)電場(chǎng)，不能用于非均勻風(fēng)電場(chǎng)。本文對(duì)PGNN模型進(jìn)行改進(jìn)，利用三維高斯尾流模型作為物理信息基函數(shù)，并采用適用于非均勻風(fēng)電場(chǎng)的圖神經(jīng)網(wǎng)絡(luò)更新策略，提出一種改進(jìn)的物理引導(dǎo)的圖神經(jīng)網(wǎng)絡(luò)（improved physics-induced graph neural network，IPGNN）模型對(duì)多種風(fēng)況下的非均勻風(fēng)電場(chǎng)進(jìn)行功率預(yù)測(cè)，并通過(guò)與PGNN模型和數(shù)據(jù)驅(qū)動(dòng)方法對(duì)比來(lái)證實(shí)其性能優(yōu)勢(shì)。構(gòu)建一個(gè)基于三維高斯尾流模型的物理信息基函數(shù)，更準(zhǔn)確地反映非均勻風(fēng)電場(chǎng)中風(fēng)力機(jī)之間尾流相互作用關(guān)系，可有效提高預(yù)測(cè)精度；設(shè)計(jì)一組適用于非均勻風(fēng)電場(chǎng)的圖神經(jīng)網(wǎng)絡(luò)更新策略，將風(fēng)力機(jī)之間的相互尾流作用和注意力機(jī)制嵌入到圖神經(jīng)網(wǎng)絡(luò)更新過(guò)程中，使模型具有較強(qiáng)的可解釋性。實(shí)驗(yàn)結(jié)果表明，與PGNN模型和數(shù)據(jù)驅(qū)動(dòng)方法相比，該方法在多種風(fēng)況下的預(yù)測(cè)誤差均更小。

1 相關(guān)工作

1.1 預(yù)測(cè)發(fā)電量的工程模型

在風(fēng)電場(chǎng)中，風(fēng)力機(jī)發(fā)電量的主要影響是上游風(fēng)力機(jī)葉片旋轉(zhuǎn)產(chǎn)生的尾流和湍流對(duì)下游風(fēng)力機(jī)產(chǎn)生的不利影響，大多數(shù)方法采用工程尾流模型量化上游風(fēng)力機(jī)影響下的穩(wěn)態(tài)風(fēng)速來(lái)計(jì)算目標(biāo)風(fēng)力機(jī)的發(fā)電量，因此尾流模型的設(shè)計(jì)與構(gòu)建一直都是研究重點(diǎn)。目前風(fēng)電場(chǎng)布局中應(yīng)用最廣泛的是一維Jensen模型［11］，該模型在計(jì)算上較為簡(jiǎn)單但精度不高。后來(lái)陸續(xù)有二維尾流模型，不僅通過(guò)實(shí)驗(yàn)證明優(yōu)于一維模型，且成功地應(yīng)用于風(fēng)電場(chǎng)優(yōu)化問(wèn)題中［12］。但如果選擇安裝不同尺寸的風(fēng)力機(jī)，由于需將輪轂高度考慮進(jìn)模型中，因此一維與二維模型都不再適用。文獻(xiàn)［13］在二維Jensen模型基礎(chǔ)上提出一種三維尾流模型，增加了垂直風(fēng)向與風(fēng)切參考量，并通過(guò)多個(gè)實(shí)驗(yàn)證實(shí)其精確性。文獻(xiàn)［14］提出三維高斯尾流模型，證實(shí)其可用于非均勻風(fēng)電場(chǎng)布局優(yōu)化問(wèn)題且其預(yù)測(cè)與真實(shí)情況接近。圖1為三維高斯尾流模型示意圖，展示了尾流效應(yīng)對(duì)風(fēng)速的影響。

由美國(guó)國(guó)家可再生能源實(shí)驗(yàn)室（NREL）開(kāi)發(fā)并維護(hù)的用于風(fēng)電場(chǎng)控制建模、設(shè)計(jì)和分析的開(kāi)源框架FLORIS［15］在最新版本中改進(jìn)了針對(duì)大型風(fēng)力機(jī)陣列以及非均勻流場(chǎng)和二次尾流轉(zhuǎn)向效應(yīng)的研究方法，添加了用于非均勻風(fēng)電場(chǎng)研究的GCH（Gauss curl hybrid）［16］三維尾流模型，其考慮了更多的影響因素。這是FLORIS首次加入針對(duì)非均勻風(fēng)電場(chǎng)的模型，計(jì)算精度有很大提升，預(yù)測(cè)的平均絕對(duì)誤差最高可降低18.9%［17］。本文研究數(shù)據(jù)生成和網(wǎng)絡(luò)模型搭建中都使用該尾流模型，三維高斯尾流模型采用三維坐標(biāo)表示風(fēng)力機(jī)位置，可計(jì)算風(fēng)力機(jī)轉(zhuǎn)子中尾流任意點(diǎn)的流向速度虧損。風(fēng)力機(jī)轉(zhuǎn)子后面的流向速度[u（x，y，z）]的計(jì)算公式為：

[u（x，y，z）=U∞1-Cexp-（y-δ）22σ2y?exp-（z-zh）22σ2y]"" （1）

[C=1-1-（σy0σz0）C0（2-C0）σyσz]"" （2）

[C0=1-1-CT]" （3）

[σzD=kz（x-x0）D+σz0D，"" σz0D=12uRU∞+u0]"" （4）

[σyD=kz（x-x0）D+σy0D，"" σy0D=σz0Dcosγ]""" （5）

式中：[x、][y、][z]——在流向、橫向和垂直方向上的空間坐標(biāo)；[U∞]——自由氣流下風(fēng)速；[C]——尾流中心速度損失；[δ]——尾流偏轉(zhuǎn)；[σy]、[σz]——[y]、[z]方向尾流寬度，它們?nèi)Q于尾流膨脹率和推力系數(shù)[CT]；[zh]——輪轂高度；下標(biāo)0——該項(xiàng)在遠(yuǎn)尾流開(kāi)始時(shí)的初值；[D]——轉(zhuǎn)子直徑；[ky]、[kz]——尾流在橫向和縱向上的膨脹率，此處考慮由風(fēng)力機(jī)運(yùn)轉(zhuǎn)引起的額外湍流，以更準(zhǔn)確地計(jì)算尾流膨脹率；[uR]——轉(zhuǎn)子的速度；[u0]——尾流中最大虧損速度；[γ]——風(fēng)力機(jī)的偏航偏移角度；其他參數(shù)具體含義在文獻(xiàn)［17］中有詳細(xì)解釋。

1.2 PGNN模型

文獻(xiàn)［10］提出的PGNN模型結(jié)合了物理信息網(wǎng)絡(luò)和圖神經(jīng)網(wǎng)的優(yōu)勢(shì)，在均勻風(fēng)電場(chǎng)發(fā)電量上取得了不錯(cuò)的效果。PGNN將當(dāng)前風(fēng)電場(chǎng)配置用圖結(jié)構(gòu)表示，把風(fēng)力機(jī)看作圖的節(jié)點(diǎn)，把風(fēng)力機(jī)相互尾流影響關(guān)系看作圖的邊。為處理圖結(jié)構(gòu)表示的風(fēng)電場(chǎng)，設(shè)計(jì)一組更新策略，通過(guò)在風(fēng)電場(chǎng)圖的節(jié)點(diǎn)之間發(fā)送和更新消息來(lái)更新風(fēng)電場(chǎng)圖的邊、節(jié)點(diǎn)和全局特征，將輸入圖轉(zhuǎn)化為節(jié)點(diǎn)嵌入向量集；又通過(guò)全連接層將節(jié)點(diǎn)嵌入向量作為輸入來(lái)預(yù)測(cè)各風(fēng)力機(jī)的輸出功率，用于計(jì)算整個(gè)風(fēng)電場(chǎng)的發(fā)電量。其中，在消息傳遞、更新過(guò)程中，更新函數(shù)和過(guò)程受到尾流模型的函數(shù)形式和計(jì)算過(guò)程的啟發(fā)，把用于計(jì)算風(fēng)力機(jī)發(fā)電量的一維Park尾流模型做適當(dāng)變換，構(gòu)造一個(gè)帶有一組可訓(xùn)練參數(shù)的物理信息基函數(shù)，將其作為消息傳遞中的權(quán)重更新函數(shù)用于在更新過(guò)程中為各節(jié)點(diǎn)分配相應(yīng)權(quán)重，即調(diào)節(jié)各風(fēng)力機(jī)受到的尾流影響相互作用的“強(qiáng)度”，類似在物理上上游風(fēng)力機(jī)在工作時(shí)對(duì)下游風(fēng)力機(jī)所產(chǎn)生的尾流效應(yīng)。同時(shí)，為增強(qiáng)模型的泛化能力，PGNN采用隨機(jī)生成的風(fēng)電場(chǎng)及風(fēng)況作為訓(xùn)練數(shù)據(jù)，使得模型在多種風(fēng)況下均適用。PGNN量化了風(fēng)力機(jī)之間的尾流相互作用，使其在物理上具有較強(qiáng)的可解釋性。模擬實(shí)驗(yàn)表明，PGNN在多種風(fēng)電場(chǎng)條件下的效果均優(yōu)于數(shù)據(jù)驅(qū)動(dòng)模型，證實(shí)能夠用于實(shí)際風(fēng)電場(chǎng)發(fā)電量預(yù)測(cè)工作。

2 模型構(gòu)建

2.1 構(gòu)建物理信息基函數(shù)

PGNN模型未將輪轂高度考慮進(jìn)尾流影響中，對(duì)具有不同尺寸風(fēng)力機(jī)的風(fēng)電場(chǎng)并不適用。因此，本節(jié)采用1.1節(jié)提到的適用于非均勻風(fēng)電場(chǎng)的三維高斯尾流模型，做出相應(yīng)變換，構(gòu)造本文模型所需物理信息基函數(shù)。根據(jù)式（1），經(jīng)過(guò)尾流效應(yīng)后的風(fēng)速是由初始風(fēng)速減去經(jīng)過(guò)尾流效應(yīng)的損失得到的，因此在構(gòu)建物理信息基函數(shù)時(shí)只需考慮尾流效應(yīng)損失部分，即：

[L（x，y，z）=Cexp-（y-δ）22σy2?exp-（z-zh）22σz2]"""" （6）

本文中，各風(fēng)力機(jī)的轉(zhuǎn)子中心位置用三維坐標(biāo)表示，設(shè)定上游風(fēng)力機(jī)為[j]，下游風(fēng)力機(jī)為[i]，式（4）、式（6）中[（x-x0）]、[（y-δ）]、[（z-zh）]的具體形式則為[（xi-xj）]、[（yi-yj）]、[（zi-zj）]，為方便記錄和表示，分別記作[▽xij、][▽yij]和[▽zij]。此時(shí)[L（x，y，z）]的表達(dá)式變?yōu)椋?/p>

[L▽xij，▽yij，▽zij=Cexp-▽2yij2σy2?exp-▽2zij2σz2]""""" （7）

式（7）是一個(gè)帶有一組可訓(xùn)練參數(shù)的指數(shù)函數(shù)，然而指數(shù)函數(shù)在計(jì)算函數(shù)梯度時(shí)可能會(huì)發(fā)生梯度消失或者梯度爆炸，所以為了降低計(jì)算梯度時(shí)的不穩(wěn)定性，這里將其近似變?yōu)閮缂?jí)數(shù)形式：

[L▽xij，▽yij，▽zij≈Cn=05-▽2yij2σy2n/n！?n=05-▽2zij2σz2n/n！]"" （8）

此外，需增加一個(gè)約束來(lái)判斷風(fēng)力機(jī)是否受到尾流影響，以此決定尾流模型公式的使用。受到尾流影響必須滿足3個(gè)條件：1）該風(fēng)力機(jī)位于對(duì)其有尾流作用風(fēng)力機(jī)的下游；2）兩臺(tái)風(fēng)力機(jī)的偏移歐氏距離處在尾流模型的偏移范圍內(nèi)，即[δijlt;δmax]；3）風(fēng)力機(jī)偏航角不超過(guò)尾流影響的最大角度，即[γijlt;γmax]。其中偏移歐氏距離選用較大風(fēng)力機(jī)轉(zhuǎn)子的直徑作為基數(shù)。

2.2 GN層更新策略

GN層利用邊緣更新函數(shù)[Φe（?";θ1）]、節(jié)點(diǎn)更新函數(shù)[Φv（?";θ2）]、全局更新函數(shù)[Φu（?";θ3）]和權(quán)重更新函數(shù)[Φw（?";θ4）]其4個(gè)函數(shù)在消息傳遞過(guò)程中完成信息的更新。其中邊緣、節(jié)點(diǎn)和全局更新函數(shù)又分別包含各自的聚合函數(shù)[ρe（?）]、[ρv（?）]、[ρu（?）]，用于聚合所收到的信息（[e、][v、][u]分別代表邊緣特征、節(jié)點(diǎn)特征和全局特征）。權(quán)重更新函數(shù)采用2.1節(jié)中構(gòu)建的物理信息基函數(shù)。GN層更新過(guò)程如圖2所示。

2.2.1 風(fēng)電場(chǎng)的圖表示

有向圖[G=（V，E，g）]是包含節(jié)點(diǎn)特征、邊緣特征和全局特征的一個(gè)元組。設(shè)[V、][E]是節(jié)點(diǎn)和邊緣特征的任意無(wú)序子集，節(jié)點(diǎn)特征[vi∈V]是包含特定節(jié)點(diǎn)屬性的任意向量，邊緣特征[eij∈E]表示節(jié)點(diǎn)[j]到[i]的邊屬性。這種邊之間的關(guān)系稱為發(fā)送方[j]和接收方[i]。[I]是包含所有節(jié)點(diǎn)索引的索引集，[U]是包含所有邊索引[（i，j）]的索引集，[g]是描述圖[G]全局屬性的向量。

2.2.2 邊緣更新

邊緣更新函數(shù)[Φe（?";θ1）]根據(jù)式（9）更新圖的所有邊。它利用發(fā)送方節(jié)點(diǎn)特征[vj、]接收方節(jié)點(diǎn)特征[vi]以及兩個(gè)節(jié)點(diǎn)之間邊的特征[eij]生成更新后的邊的特征[eij′]。GN層利用權(quán)重計(jì)算函數(shù)[Φw（?";θ4）]根據(jù)風(fēng)力機(jī)的相對(duì)位置和風(fēng)力條件動(dòng)態(tài)調(diào)整[eij′]的權(quán)重。這里選擇風(fēng)力機(jī)三維坐標(biāo)[（x，y，z）]作為[Φw（?";θ4）]的輸入，[Φw（?";θ4）]計(jì)算權(quán)重[Wij]，用于得到更新后的邊[eij′]。利用尾流模型中的影響關(guān)系得到權(quán)重來(lái)調(diào)節(jié)風(fēng)力機(jī)之間影響的強(qiáng)度，這種計(jì)算在物理上有意義并具備較強(qiáng)的可解釋性。邊緣更新過(guò)程為：

[eij′←Φe（eij，vj，vi;θ1）Wij←Φw▽xij，▽yij，▽zij;θ4eij′←Wij×eij′]""nbsp; （9）

2.2.3 節(jié)點(diǎn)更新

節(jié)點(diǎn)更新函數(shù)[Φv（?";θ2）]使用當(dāng)前輸入節(jié)點(diǎn)特征[vi、]更新后的邊特征[eij′]和全局特征[g]更新節(jié)點(diǎn)特征。GN層使用邊聚合函數(shù)[ρe（?）]將邊緣特征集合減少為單個(gè)聚合的邊緣特征[eij′]。節(jié)點(diǎn)更新函數(shù)[Φv（?;θ2）]處理聚合后的邊緣特征[eij′]和當(dāng)前節(jié)點(diǎn)特征[vi]以產(chǎn)生更新后的節(jié)點(diǎn)特征[vi′]。邊緣聚合和節(jié)點(diǎn)更新過(guò)程為：

[ei′←ρe（Ei′），" Ei′=eij′?j∈Rivi′←Φv（ei′，vi，g;θ2）]""""" （10）

2.2.4 全局更新

全局更新函數(shù)[Φu（?";θ3）]用來(lái)更新全局特征。GN使用更新后的邊和節(jié)點(diǎn)特征來(lái)更新全局特征。為此，GN利用聚合函數(shù)[ρu（?）]和[ρv（?）]將所有邊緣特征[E′]和所有節(jié)點(diǎn)特征[V′]聚合成兩個(gè)向量[e′]和[v′]。全局特征更新過(guò)程為：

[E′={eij′|?（i，j）∈U}，" e′←ρu（E′）V′={ni′|?i∈I}，" e′←ρv（V′）g′←Φu（e′，v′，g;θ3）]" （11）

2.2.5 權(quán)重更新

模型中，邊緣更新函數(shù)[Φe（?";θ1）]、節(jié)點(diǎn)更新函數(shù)[Φv（?";θ2）]、全局更新函數(shù)[Φu（?";θ3）]都是多層感知機(jī)（MLP）。權(quán)重更新函數(shù)[Φw（?";θ4）]用于圖神經(jīng)網(wǎng)絡(luò)中節(jié)點(diǎn)之間邊緣信息的更新，表達(dá)式即為2.1節(jié)所構(gòu)建的物理信息基函數(shù)：

[Φw▽xij，▽yij，▽zij≈Cn=05-▽2yij2σy2nn！?n=05-▽2zij2σz2nn！]"""""" （12）

通過(guò)將這些更新函數(shù)依次應(yīng)用于圖的輸入，GN生成更新后的圖[G′=（V′，E′，g′）]。該方法通過(guò)優(yōu)化參數(shù)學(xué)習(xí)如何表示風(fēng)力機(jī)間的相互作用。使所有的更新函數(shù)共享參數(shù)，這種方法在不同風(fēng)電場(chǎng)布局和風(fēng)況下的輸出功率預(yù)測(cè)中具有更好的通用性。GN層完整計(jì)算過(guò)程為：

[algorithm Graph Network[（E，V，u）]

for[（i，j）]in[U]do

[eij′←Φe（eij，vj，vi;θ1）]"""" #更新邊緣特征

end for

for[i] in[I]do

for [j]in[Ri] do

[Wij←Φw（▽xij，▽yij，▽zij;θ4）] #計(jì)算權(quán)重

[eij′←Wij×eij′]"""" #分配權(quán)重

end for

let[Ei′=eji′?j∈Ri]

[ei′←ρe（Ei′）]""""""""" #聚合邊緣特征

[vi′←Φv（ei′，vi，g;θ2）]#更新節(jié)點(diǎn)特征

end for

let[E′={eij′?（i，j）∈U}]

[e′←ρu（E′）]"""""" #聚合、更新邊緣特征

let [V′={ni′?i∈I}]

[e′←ρv（V′）]"""""" #聚合、更新節(jié)點(diǎn)特征

[g′←Φu（e′，v′，g;θ3）]""" #更新全局特征

[G′←（E′，V′，g′）]"""""" #生成更新后的圖

end algorithm ]

2.3 IPGNN實(shí)現(xiàn)過(guò)程

IPGNN由圖網(wǎng)絡(luò)層和全連接層組成，其完整實(shí)現(xiàn)過(guò)程如圖3所示。每個(gè)GN層通過(guò)在風(fēng)電場(chǎng)圖的節(jié)點(diǎn)之間發(fā)送和更新消息來(lái)更新風(fēng)電場(chǎng)圖的邊、節(jié)點(diǎn)和全局特征以預(yù)測(cè)各風(fēng)力機(jī)的發(fā)電量。模型中使用3個(gè)GN層嵌入風(fēng)電場(chǎng)圖，將每個(gè)內(nèi)部節(jié)點(diǎn)、邊緣和全局特征設(shè)置為100維的向量。GN層的3個(gè)更新函數(shù)[Φe（?";θ1）]、[Φv（?";θ2）]、[Φu（?";θ3）]設(shè)計(jì)由528、528、256維隱藏層的MLP以及Leaky ReLU激活單元組成。使用一個(gè)全連接層將處理后的嵌入向量映射到所有風(fēng)力機(jī)的功率輸出（即功率輸出圖）。全連接層與GN層有相同的MLP架構(gòu)，輸出激活單元也為ReLU。經(jīng)過(guò)GN層計(jì)算，生成更新后的風(fēng)電場(chǎng)圖[G′]。全連接層將[G′]作為輸入，功率[P]作為輸出，其節(jié)點(diǎn)值為風(fēng)電場(chǎng)中單臺(tái)風(fēng)力機(jī)的輸出功率。輸出功率圖[Pi]采取歸一化策略，計(jì)算出的功率除以未經(jīng)過(guò)尾流影響所能產(chǎn)生的最大功率。所有風(fēng)力機(jī)的歸一化功率構(gòu)成輸出圖[Pi]的節(jié)點(diǎn)特征。優(yōu)化器選用Adam優(yōu)化器；學(xué)習(xí)率與Park模型相同，設(shè)置為0.001。

3 實(shí)驗(yàn)方案

3.1 風(fēng)力機(jī)類型選擇

在實(shí)驗(yàn)中，模擬生成的非均勻風(fēng)電場(chǎng)由兩種不同型號(hào)（分別為NREL 5 MW和IEA 10 MW）的風(fēng)力機(jī)組成，其的部分重要參數(shù)如表2所示。NREL 5 MW風(fēng)力機(jī)是目前應(yīng)用最為廣泛、行業(yè)標(biāo)準(zhǔn)制定最為完善的中型風(fēng)力機(jī)之一。IEA 10 MW風(fēng)力機(jī)是國(guó)際能源署（IEA）為滿足當(dāng)前對(duì)風(fēng)電機(jī)組發(fā)電功率不斷上升的需求而推出的一種大型風(fēng)力機(jī)，研究表明該風(fēng)力機(jī)是當(dāng)下及未來(lái)部署大型風(fēng)電場(chǎng)的較好選擇［18］。因此，本文選用上述兩種風(fēng)力機(jī)模擬非均勻風(fēng)電場(chǎng)進(jìn)行實(shí)驗(yàn)。

3.2 數(shù)據(jù)準(zhǔn)備

本文的數(shù)據(jù)生成需用到GraphFLORIS開(kāi)源框架。GraphFLORIS是繼承于FLORIS的可用于生成風(fēng)電場(chǎng)圖形數(shù)據(jù)的工具包，已被成功應(yīng)用于解決均勻風(fēng)電場(chǎng)發(fā)電功率預(yù)測(cè)問(wèn)題。而對(duì)于非均勻風(fēng)電場(chǎng)，這里需將其尾流模型變換為三維高斯尾流模型以生成IPGNN訓(xùn)練所需數(shù)據(jù)。訓(xùn)練數(shù)據(jù)是由成對(duì)的風(fēng)電場(chǎng)圖[Gi]和輸出功率圖[Pi]表示。輸入圖[gi]作為輸入變量是由隨機(jī)生成的風(fēng)電場(chǎng)和隨機(jī)風(fēng)況（風(fēng)速、風(fēng)向）組成的。初始風(fēng)速根據(jù)風(fēng)力機(jī)的特性使用均勻分布，設(shè)置風(fēng)速為[s～U（6.0" m/s，15.0" m/s）]，風(fēng)向?yàn)閇θ～U（0°，"360°）]。

3.3 訓(xùn)練目標(biāo)

使用目標(biāo)風(fēng)電場(chǎng)中所有風(fēng)力渦輪機(jī)的估計(jì)功率和目標(biāo)功率之間的均方誤差（mean-square error，MSE）作為IPGNN的損失函數(shù)，即為模型的訓(xùn)練目標(biāo)。

4 預(yù)測(cè)結(jié)果及分析

訓(xùn)練完成的風(fēng)電場(chǎng)功率預(yù)測(cè)模型用于以圖的形式在任意風(fēng)向[θ]和任意風(fēng)速[s]下對(duì)包含任意數(shù)量[n]的風(fēng)力機(jī)、任意布局的風(fēng)電場(chǎng)[l]估計(jì)所有風(fēng)力機(jī)的發(fā)電功率。為了驗(yàn)證模型性能，進(jìn)行隨機(jī)布局條件下風(fēng)速和風(fēng)向泛化實(shí)驗(yàn)。在評(píng)估時(shí)，使用MSE來(lái)評(píng)估模型性能。MSE主要用于評(píng)價(jià)預(yù)測(cè)精度，平均絕對(duì)百分比誤差（mean absolute percentage error，MAPE）主要用于評(píng)價(jià)結(jié)果好壞。包含預(yù)測(cè)功率[P′]和實(shí)際功率[P]的[MMSE（P，P′）]計(jì)算為：

[MMSE（P，P′）=1ni=1nP（i）-P（i）′2]"" （13）

式中：[P（i）]——第[i]個(gè)圖節(jié)點(diǎn)；[n]——節(jié)點(diǎn)數(shù)量，即風(fēng)力機(jī)數(shù)量。

為了減小一個(gè)風(fēng)電場(chǎng)估計(jì)精度造成的誤差影響，將采用10個(gè)風(fēng)電場(chǎng)布局的平均誤差作為估計(jì)結(jié)果誤差，其計(jì)算公式為：

[MMSE=j=1ljMMSE（P，P′）lj]"""" （14）

式中：[lj]——第[j]個(gè)風(fēng)電場(chǎng)，本文實(shí)驗(yàn)中設(shè)定[j=10]。

訓(xùn)練好的IPGNN在不同風(fēng)電場(chǎng)布局和不同風(fēng)速和風(fēng)向的情況下可很好地估計(jì)功率。測(cè)試中，把風(fēng)向每[10°]設(shè)置為一個(gè)測(cè)試單元，共設(shè)置36個(gè)風(fēng)向，風(fēng)速為[6.0～15.0" m/s]，之所以將風(fēng)速設(shè)置在這個(gè)區(qū)間內(nèi)，是因?yàn)轱L(fēng)力機(jī)在這個(gè)風(fēng)速區(qū)間內(nèi)可正常運(yùn)行，風(fēng)速過(guò)低與過(guò)高都可能導(dǎo)致風(fēng)力機(jī)停止工作。圖4用熱力圖的形式展示了估計(jì)結(jié)果誤差，設(shè)置風(fēng)力機(jī)數(shù)量為20（NREL 5 MW風(fēng)力機(jī)和IEA 10 MW風(fēng)力機(jī)各10臺(tái)）?？煽闯?，在風(fēng)速為[6.0～15.0 "m/s]的情況下，風(fēng)速越快經(jīng)過(guò)訓(xùn)練的網(wǎng)絡(luò)預(yù)測(cè)誤差越小。這種結(jié)果符合一般規(guī)律，因?yàn)楸M管風(fēng)速大的情況下湍流強(qiáng)度高，但氣流中已包含足夠的風(fēng)能來(lái)產(chǎn)生渦輪機(jī)的最大功率，所以對(duì)于高風(fēng)速情況估計(jì)誤差極小。對(duì)于低風(fēng)速的情況，部分風(fēng)向情況的誤差可能會(huì)略大，后續(xù)可通過(guò)調(diào)整模型參數(shù)等方式進(jìn)行控制，且這種誤差在可接受范圍內(nèi)。

5 對(duì)比實(shí)驗(yàn)

權(quán)重更新函數(shù)為每條邊分配連續(xù)的權(quán)重值，且傾向于將較大的權(quán)重放在處于強(qiáng)尾流相互作用下的兩臺(tái)風(fēng)力機(jī)之間構(gòu)成的邊。為了更加突出顯示其優(yōu)點(diǎn)，這里將其性能與另一個(gè)數(shù)據(jù)驅(qū)動(dòng)的網(wǎng)絡(luò)進(jìn)行對(duì)比，該網(wǎng)絡(luò)可稱為數(shù)據(jù)驅(qū)動(dòng)圖神經(jīng)網(wǎng)絡(luò)（DGNN）。DGNN利用512、512和256維的MLP以及Leaky ReLU激活函數(shù)來(lái)構(gòu)造[Φw（?";θ5）]，DGNN的[Φw（?";θ5）]輸出維度為1，與IPGNN類似，使用ReLU作為輸出激活函數(shù)，以產(chǎn)生與IPGNN的[Φw（?";θ4）]相同的輸出值范圍。為了更加顯著地對(duì)兩種網(wǎng)絡(luò)的預(yù)測(cè)結(jié)果量化比較，采用平均絕對(duì)百分比誤差進(jìn)行評(píng)估，計(jì)算公式為：

[MMAPE（P，P′）=100%ni=1nP（i）-P（i）′P（i）]"""" （15）

與式（15）類似，這里也采用10個(gè)風(fēng)電場(chǎng)布局的平均誤差作為估計(jì)結(jié)果誤差，計(jì)算公式為：

[MMAPE=j=1ljMMAPE（P，P′）lj]"""""" （16）

對(duì)比實(shí)驗(yàn)在多種數(shù)量風(fēng)力機(jī)場(chǎng)景中測(cè)試IPGNN和DGNN性能。同時(shí)為證實(shí)三維尾流模型在非均勻風(fēng)電場(chǎng)的作用，把采用一維尾流模型的Park模型用于預(yù)測(cè)非均勻風(fēng)電場(chǎng)發(fā)電功率并將結(jié)果與IPGNN對(duì)比。對(duì)比實(shí)驗(yàn)中的風(fēng)力機(jī)數(shù)量分別為20、40、60、80、100臺(tái)，這些數(shù)量基本可涵蓋小型到大中型風(fēng)電場(chǎng)規(guī)模，其中兩種類型風(fēng)力機(jī)的數(shù)量各為總數(shù)一半。風(fēng)向設(shè)置為90°，風(fēng)速設(shè)置為11.0 m/s，該風(fēng)速接近于兩種風(fēng)力機(jī)的額定風(fēng)速，在不經(jīng)過(guò)尾流作用下可實(shí)現(xiàn)最大化發(fā)電量。3種模型誤差結(jié)果如圖5所示。

從圖5可看到，對(duì)于各種數(shù)量風(fēng)力機(jī)的情況，IPGNN均比DGNN的[MMAPE]值小，其中在風(fēng)力機(jī)數(shù)量較?。ㄈ?0、20臺(tái)）的情況下，兩種模型的[MMAPE]值并不高且差別不大，分別為5.18%和3.92%，最小相差僅為1.26%。但隨著風(fēng)力機(jī)數(shù)量的增加，兩種模型的[MMAPE]值均有較大上升，但DGNN明顯比IPGNN高，其中80臺(tái)風(fēng)力機(jī)時(shí)最大相差為7.20%。而Park模型在各種風(fēng)力機(jī)數(shù)量下其[MMAPE]值均遠(yuǎn)大于另外兩種模型，這是因?yàn)橐痪S模型本身對(duì)非均勻風(fēng)電場(chǎng)并不適用，所以預(yù)測(cè)精度和可靠性不高，也證實(shí)三維模型的使用對(duì)非均勻風(fēng)電場(chǎng)功率預(yù)測(cè)有較好的適用性。

IPGNN雖然在各種數(shù)量下其[MMAPE]值都是3種模型中最小的，但在80臺(tái)和100臺(tái)風(fēng)力機(jī)下其值分別達(dá)到19.61%和30.32%，這表明模型對(duì)大型風(fēng)電場(chǎng)的預(yù)測(cè)效果并不十分理想。隨著風(fēng)力機(jī)數(shù)量的增加，各風(fēng)力機(jī)在網(wǎng)絡(luò)中傳遞更新次數(shù)成倍增加，這也導(dǎo)致原始誤差被放大多倍，造成最終預(yù)測(cè)誤差隨之增大。該實(shí)驗(yàn)結(jié)果也符合一般規(guī)律，即在風(fēng)力機(jī)數(shù)量少的情況下平均預(yù)測(cè)誤差較小，隨著風(fēng)力機(jī)數(shù)量的增大，平均預(yù)測(cè)誤差也增大。表3列出了在20臺(tái)風(fēng)力機(jī)情況下歸一化的3種模型的預(yù)測(cè)值與真實(shí)值。實(shí)驗(yàn)表明，IPGNN在估計(jì)風(fēng)電場(chǎng)功率性能方面具有較好的性能。

對(duì)于一個(gè)擁有50臺(tái)風(fēng)力機(jī)規(guī)模的風(fēng)電場(chǎng)，IPGNN需要約0.12 s來(lái)估計(jì)50臺(tái)風(fēng)力機(jī)的功率。雖然在評(píng)估時(shí)使用單線程運(yùn)行程序，但網(wǎng)絡(luò)模型在訓(xùn)練階段采用的分布式計(jì)算框架ray可使模型并行運(yùn)算以大量減少運(yùn)行時(shí)間。訓(xùn)練好的模型在功率估計(jì)階段的執(zhí)行時(shí)間大致隨風(fēng)力機(jī)數(shù)量呈二次方增加，時(shí)間復(fù)雜度約為[O（n2）]。該復(fù)雜度可根據(jù)圖的本身特點(diǎn)進(jìn)行定性分析：圖網(wǎng)絡(luò)中有[n]臺(tái)風(fēng)力機(jī)，即[n]個(gè)頂點(diǎn)，風(fēng)電場(chǎng)圖最多有[n2]條邊。而且編碼器中的前向傳播時(shí)間幾乎是恒定的，這里可忽略節(jié)點(diǎn)更新的計(jì)算復(fù)雜性，因?yàn)樗请S風(fēng)力機(jī)數(shù)量線性增加的，所以IPGNN時(shí)間復(fù)雜度約為[O（n2）]。

6 應(yīng) 用

從多數(shù)現(xiàn)有風(fēng)電場(chǎng)布局來(lái)看，風(fēng)電場(chǎng)選址在相對(duì)平坦地面時(shí)多數(shù)基于類似網(wǎng)格的規(guī)整布局。為了研究IPGNN在規(guī)整布局中的表現(xiàn)，模擬一個(gè)具有33臺(tái)風(fēng)力機(jī)的風(fēng)電場(chǎng)，其中NREL 5 MW風(fēng)力機(jī)24臺(tái)，IEA 10 MW風(fēng)力機(jī)9臺(tái)，如圖6所示。這種布局思路來(lái)自文獻(xiàn)［19］，在風(fēng)電場(chǎng)外圍布局一種小型風(fēng)力機(jī)，在內(nèi)側(cè)布局另一種相對(duì)大型風(fēng)力機(jī)。本文選擇NREL 5 MW風(fēng)力機(jī)在外側(cè)，因?yàn)槠漭嗇灨叨葹?0 m，而IEA 10 MW風(fēng)力機(jī)的輪轂高度為119 m，這種高度差使得部分垂直方向的風(fēng)不會(huì)受到尾流影響，從而使內(nèi)部風(fēng)力機(jī)所受尾流損失更低。在此風(fēng)電場(chǎng)布局中，風(fēng)力機(jī)間距為[5D]（[D]為渦輪轉(zhuǎn)子直徑）。風(fēng)速設(shè)置為11.0 m/s，風(fēng)向?yàn)?°～360°，增量為10°，然后將IPGNN和DGNN兩種模型的估計(jì)值與相應(yīng)目標(biāo)值進(jìn)行比較。

圖7為21號(hào)風(fēng)力機(jī)在風(fēng)速為10 m/s的條件下，各風(fēng)向下IPGNN和DGNN經(jīng)歷尾流效應(yīng)時(shí)IPGNN和DGNN的發(fā)電功率預(yù)測(cè)結(jié)果。從圖7可看到，在經(jīng)歷較強(qiáng)尾流效應(yīng)0°、90°、180°、270°時(shí)，IPGNN和DGNN都能較好地預(yù)測(cè)發(fā)電功率。但當(dāng)風(fēng)向偏移，風(fēng)力機(jī)所經(jīng)歷復(fù)雜尾流效應(yīng)時(shí)，DGNN往往對(duì)發(fā)電量估計(jì)過(guò)高，而IPGNN與目標(biāo)值接近，仍可較好地預(yù)測(cè)發(fā)電功率。因此，IPGNN相比DGNN不僅預(yù)測(cè)結(jié)果更準(zhǔn)確且其泛化能力也較強(qiáng)。

7 結(jié) 論

通過(guò)訓(xùn)練模型和進(jìn)行對(duì)比實(shí)驗(yàn)，可得三維高斯尾流模型作為物理信息基函數(shù)被用于圖神經(jīng)網(wǎng)絡(luò)模型，且其在非均勻風(fēng)電場(chǎng)中的效果優(yōu)于一維模型；所提出的更新策略能準(zhǔn)確地傳遞、更新非均勻風(fēng)電場(chǎng)中各風(fēng)力機(jī)受到的尾流影響程度，同時(shí)使得模型可解釋性較好。與數(shù)據(jù)驅(qū)動(dòng)方法相比，物理引導(dǎo)的圖神經(jīng)網(wǎng)絡(luò)模型可利用較少的數(shù)據(jù)對(duì)非均勻風(fēng)電場(chǎng)發(fā)電功率有更小的預(yù)測(cè)誤差，從而利于風(fēng)電場(chǎng)發(fā)電功率預(yù)測(cè)與優(yōu)化工作。

［參考文獻(xiàn)］

［1］"""" WISER R， LANTZ E， MAI T， et al. Wind vision： a new era for wind power in the United States［J］. The electricity journal， 2015， 28（9）： 120-132.

［2］"""" Global Wind Energy Council. GWEC global wind report［R］. Brussels， Belgium： Global Wind Energy Council， 2022. 1-148.

［3］"""" BARTHELMIE R J， HANSEN K， FRANDSEN S T， et al. Modelling and measuring flow and wind turbine wakes in large wind farms offshore［J］. Wind energy， 2009， 12（5）： 431-444.

［4］"""" CHEN Y， LI H， JIN K， et al. Wind farm layout optimization using genetic algorithm with different hub height""" wind""" turbines［J］.""" Energy""" conversion""" and management， 2013， 70： 56-65.

［5］"""" MOSETTI G， POLONI C， DIVIACCO B. Optimization of wind turbine positioning in large windfarms by means of a genetic"" algorithm［J］." Journal" of"" wind" engineering" and industrial aerodynamics， 1994， 51（1）： 105-116.

［6］"""" HOU P， HU W H， SOLTANI M， et al. Optimized placement of wind turbines in large-scale offshore wind farm using particle swarm optimization algorithm［J］. IEEE transactions on sustainable energy， 2015， 6（4）： 1272-1282.

［7］"""" WANG Y， LIU H， LONG H， et al. Differential evolution with a new encoding mechanism for optimizing wind farm layout［J］. IEEE transactions on industrial informatics， 2017， 14（3）： 1040-1054.

［8］"""" 袁桂麗， 吳振民， 劉驊騏， 等. 基于深度置信網(wǎng)絡(luò)的短期風(fēng)電功率預(yù)測(cè)［J］. 太陽(yáng)能學(xué)報(bào)， 2022， 43（2）： 451-457.

YUAN G L， WU Z M， LIU H Q， et al. Short-term wind power prediction based on deep belief network［J］. Acta energiae solaris sinica， 2022， 43（2）： 451-457.

［9］"""" 蘇向敬， 周汶鑫， 李超杰， 等. 基于雙重注意力LSTM神經(jīng)網(wǎng)絡(luò)的可解釋海上風(fēng)電出力預(yù)測(cè)［J］. 電力系統(tǒng)自動(dòng)化， 2022， 46（7）： 141-151.

SU X J， ZHOU W X， LI C J， et al. Interpretable offshore wind power output forecasting based on long short-term memory" neural" network" with" dual-stage" attention［J］. Automation of electric power systems， 2022， 46（7）： 141-151.

［10］""" PARK J Y， PARK J K. Physics-induced graph neural network： an application to wind-farm power estimation［J］. Energy， 2019， 187（1）： 115883-115897.

［11］""" JENSEN N O. A note on wind generator interaction［M］. Roskilde， Denmark： Riso National Laboratory， 1983.

［12］""" TAO S Y， XU Q S， FEIJóO A， et al. Bi-hierarchy optimization of a wind farm considering environmental impact［J］." IEEE" transactions" on" sustainable" energy， 2020， 11（4）： 2515-2524.

［13］""" 宋翌蕾， 田琳琳， 趙寧. 風(fēng)力機(jī)三維尾流模型的提出與校核［J］. 太陽(yáng)能學(xué)報(bào)， 2021， 42（2）： 129-135.

SONG Y L， TIAN L L， ZHAO N. Proposal and validation of a new 3D wake model for wind turbine［J］. Acta energiae solaris sinica， 2021， 42（2）： 129-135.

［14］""" SUN H Y， YANG H X. Numerical investigation of the average wind speed of a single wind turbine and development of a novel three-dimensional multiple wind turbine wake model［J］. Renewable energy， 2020， 147： 192-203.

［15］""" FLEMING P， KING J， BAY C J， et al. Overview of FLORIS" updates［J］." Journal" of" physics：" conference series， 2020， 1618（2）： 022028.

［16］""" KING J， FLEMING P， KING R， et al. Control-oriented model for secondary effects of wake steering［J］. Wind energy science， 2021， 6（3）： 701-714.

［17］""" FARRELL A， KING J， DRAXL C， et al. Design and analysis of a wake model for spatially heterogeneous flow［J］. Wind energy science， 2021， 6（3）： 737-758.

［18］""" GUO X J， ZHANG Y， YAN J T， et al. Integrated dynamics response analysis for IEA 10-MW spar floating offshore wind turbine［J］. Journal of marine science and engineering， 2022， 10（4）： 542-559.

［19］""" STANLEY A P J， NING A. Coupled wind turbine design and layout optimization with nonhomogeneous wind turbines［J］. Wind energy science， 2019， 4（1）： 99-114.

POWER PREDICTION METHOD OF NON-UNIFORM WIND FARM

LAYOUT BASED ON GRAPH NEURAL NETWORK

Liu Deao1，Dong Minggang1，2，Ye Wei1，Wang Yan1，Gan Guojun1

（1. School"of"Information"Science"and"Engineering，"Guilin"University"of"Technology，"Guilin"541004， China；

2. Guangxi"Key"Laboratory"of"Embedded"Technology"and"Intelligent"System，"Guilin"541004，"China）

Abstract：Aiming at the problems of low accuracy and poor interpretability of the existing power prediction methods of non-uniform wind farm layout， we propose an improved physics-induced graph neural network （IPGNN） model to predict the power output of non-uniform wind farms. Firstly， IPGNN constructs a physical information basis function based on three-dimensional Gaussian wake model， which can more accurately reflect the wake interaction between wind turbines in non-uniform wind farms； Secondly， IPGNN designs a set of graph neural network update strategies based on the message passing framework. In the graph neural network edge update strategy， combined with the attention mechanism， the physical information basis function is used as the weight update function， which enhances the interpretability of the model. The simulation experiments of non-uniform wind farms with different numbers of wind turbines show that IPGNN model has a good prediction effect compared with PGNN model based on one-dimensional wake model and typical data-driven model. For non-uniform wind farms with 20 wind turbines， the MAPE of power generation prediction is 3.92%. It is an effective method for power generation prediction of non-uniform wind farm layout.

Keywords：deep learning； graph neural network； wind power； wind farm layout； numerical simulation