基于GPU加速的大電網(wǎng)N-1故障掃描批量計(jì)算方法
發(fā)布時(shí)間:2021-02-04 08:16
隨著電網(wǎng)規(guī)模的不斷擴(kuò)大,從各種可能的設(shè)備開(kāi)斷情況中篩選出嚴(yán)重故障集成為N-1安全校驗(yàn)的重要耗時(shí)部分。為了加速大電網(wǎng)N-1安全校驗(yàn)的故障篩選,提出了一種基于中央處理器-圖形處理器(CPU-GPU)異構(gòu)計(jì)算框架的實(shí)時(shí)N-1故障掃描批量計(jì)算方法?紤]到不同工況下的計(jì)算存在粗粒度并行性,進(jìn)一步挖掘計(jì)算中的細(xì)粒度并行性是提高計(jì)算效率的有效途徑。提出了同時(shí)考慮斷線(xiàn)故障和發(fā)電機(jī)開(kāi)斷故障這2種預(yù)想事故下的細(xì)粒度并行計(jì)算方法,并設(shè)計(jì)了關(guān)鍵計(jì)算步驟的核函數(shù)。增加考慮了網(wǎng)絡(luò)拓?fù)渲械囊葡嗥?使得計(jì)算精度更高,通過(guò)與IEEE標(biāo)準(zhǔn)算例和歐洲真實(shí)電網(wǎng)算例對(duì)比,驗(yàn)證了各工況下批量計(jì)算方法的正確性,并取得了顯著的加速效果。
【文章來(lái)源】:電力自動(dòng)化設(shè)備. 2020,40(08)北大核心
【文章頁(yè)數(shù)】:10 頁(yè)
【部分圖文】:
圖1 雙層并行算法架構(gòu)
以下2種情形將導(dǎo)致等效功率注入列向量Pdc發(fā)生改變:(1)斷線(xiàn)支路中含有移相器;(2)發(fā)電機(jī)開(kāi)斷故障。創(chuàng)建m+s個(gè)線(xiàn)程,前m個(gè)線(xiàn)程處理斷線(xiàn)故障,后s個(gè)線(xiàn)程處理發(fā)電機(jī)開(kāi)斷故障。針對(duì)情形(1),假設(shè)支路ij開(kāi)斷,根據(jù)上文分析Pdc中節(jié)點(diǎn)i對(duì)應(yīng)元素和節(jié)點(diǎn)j對(duì)應(yīng)元素將發(fā)生改變,節(jié)點(diǎn)i處修改變量表示為ΔPidc=-θshift(ij)/(kijxij),節(jié)點(diǎn)j處修改變量與節(jié)點(diǎn)i處修改變量互為相反數(shù),即ΔPjdc=-ΔPidc,若支路中不含移相器,則θshift(ij)=0;針對(duì)情形(2),假設(shè)節(jié)點(diǎn)a處的發(fā)電機(jī)開(kāi)斷,節(jié)點(diǎn)a處的注入功率將發(fā)生改變,進(jìn)而改變Padc。圖2為潮流方程右端項(xiàng)的并行計(jì)算過(guò)程。圖中,i(1)和j(1)分別為支路1對(duì)應(yīng)的流出和流入節(jié)點(diǎn);a(1)為發(fā)電機(jī)1對(duì)應(yīng)的節(jié)點(diǎn)編號(hào);(Pdci(1))1為場(chǎng)景1中的向量Pdc中節(jié)點(diǎn)i(1)對(duì)應(yīng)的元素,其余同理。當(dāng)線(xiàn)程編號(hào)索引小于m時(shí),由線(xiàn)程編號(hào)索引到斷線(xiàn)支路編號(hào),根據(jù)對(duì)應(yīng)的支路編號(hào)將等效功率注入列向量的第i個(gè)元素和第j個(gè)元素進(jìn)行修改;當(dāng)線(xiàn)程編號(hào)索引為m到m+s時(shí),由發(fā)電機(jī)編號(hào)索引對(duì)應(yīng)的節(jié)點(diǎn)編號(hào),對(duì)節(jié)點(diǎn)注入功率進(jìn)行修改。上述設(shè)計(jì)過(guò)程中不同的計(jì)算指令將分配到不同的線(xiàn)程束中,GPU并行處理此類(lèi)運(yùn)算,且并行度很高。
直流潮流模型平衡節(jié)點(diǎn)的相角為0,因此,對(duì)于2.5節(jié)中求得的電壓相角向量,需要在平衡節(jié)點(diǎn)的位置加零元素,其核函數(shù)設(shè)計(jì)思路如圖3所示。偽代碼見(jiàn)附錄A中表A2,令x為2.5節(jié)求得的電壓相角向量數(shù)組,新數(shù)組用x1表示。創(chuàng)建包含(m+s)×(n+1)個(gè)元素的一維數(shù)組,采用cudaMemset函數(shù)將其全部置0;然后創(chuàng)建(m+s)×n個(gè)線(xiàn)程,每個(gè)線(xiàn)程處理1個(gè)元素。CUDA異構(gòu)平臺(tái)采用單指令多線(xiàn)程架構(gòu)管理和執(zhí)行線(xiàn)程,每32個(gè)線(xiàn)程為1個(gè)線(xiàn)程束,線(xiàn)程束中所有線(xiàn)程同時(shí)執(zhí)行相同的指令,由于線(xiàn)程并不能根據(jù)批量處理規(guī)模完美分配,此處對(duì)線(xiàn)程編號(hào)設(shè)置判斷條件如下:編號(hào)小于(m+s)×n的線(xiàn)程執(zhí)行任務(wù),多余的線(xiàn)程被過(guò)濾。此時(shí),同1個(gè)線(xiàn)程束中依然執(zhí)行相同的指令,不會(huì)導(dǎo)致線(xiàn)程束的分支,保持很高的并行性。同時(shí)上述算法中內(nèi)存連續(xù)讀取,保證同一線(xiàn)程束下內(nèi)存的合并訪問(wèn),使算法具有高效性。
【參考文獻(xiàn)】:
期刊論文
[1]電力系統(tǒng)靜態(tài)安全最優(yōu)潮流并行計(jì)算方法[J]. 陽(yáng)育德,陶琢,劉輝,覃智君. 電力自動(dòng)化設(shè)備. 2019(01)
[2]一種GPU-CPU異構(gòu)運(yùn)算框架加速的實(shí)時(shí)N-1交流潮流計(jì)算方法[J]. 唐坤杰,董樹(shù)鋒,宋永華. 中國(guó)電機(jī)工程學(xué)報(bào). 2018(15)
[3]匹配GPU編程架構(gòu)的直流故障篩選算法設(shè)計(jì)[J]. 溫柏堅(jiān),謝恩彥,劉明波. 電力自動(dòng)化設(shè)備. 2017(05)
[4]基于網(wǎng)絡(luò)分塊算法的靜態(tài)安全快速計(jì)算方法[J]. 王方雨,劉文穎,田浩,梁琛,鄭晶晶. 電力自動(dòng)化設(shè)備. 2017(04)
[5]電力系統(tǒng)靜態(tài)安全混合控制方法[J]. 傅旭. 電力自動(dòng)化設(shè)備. 2017(01)
[6]基于補(bǔ)償法的預(yù)防性靜態(tài)安全控制[J]. 陽(yáng)育德,馮彥維,韋化. 電力自動(dòng)化設(shè)備. 2015(12)
[7]一種求解機(jī)組組合問(wèn)題的內(nèi)點(diǎn)半定規(guī)劃GPU并行算法[J]. 張寧宇,高山,趙欣. 電力自動(dòng)化設(shè)備. 2013(07)
[8]基于并行計(jì)算和數(shù)據(jù)復(fù)用的快速靜態(tài)安全校核技術(shù)[J]. 李峰,李虎成,於益軍,馮樹(shù)海. 電力系統(tǒng)自動(dòng)化. 2013(14)
本文編號(hào):3018019
【文章來(lái)源】:電力自動(dòng)化設(shè)備. 2020,40(08)北大核心
【文章頁(yè)數(shù)】:10 頁(yè)
【部分圖文】:
圖1 雙層并行算法架構(gòu)
以下2種情形將導(dǎo)致等效功率注入列向量Pdc發(fā)生改變:(1)斷線(xiàn)支路中含有移相器;(2)發(fā)電機(jī)開(kāi)斷故障。創(chuàng)建m+s個(gè)線(xiàn)程,前m個(gè)線(xiàn)程處理斷線(xiàn)故障,后s個(gè)線(xiàn)程處理發(fā)電機(jī)開(kāi)斷故障。針對(duì)情形(1),假設(shè)支路ij開(kāi)斷,根據(jù)上文分析Pdc中節(jié)點(diǎn)i對(duì)應(yīng)元素和節(jié)點(diǎn)j對(duì)應(yīng)元素將發(fā)生改變,節(jié)點(diǎn)i處修改變量表示為ΔPidc=-θshift(ij)/(kijxij),節(jié)點(diǎn)j處修改變量與節(jié)點(diǎn)i處修改變量互為相反數(shù),即ΔPjdc=-ΔPidc,若支路中不含移相器,則θshift(ij)=0;針對(duì)情形(2),假設(shè)節(jié)點(diǎn)a處的發(fā)電機(jī)開(kāi)斷,節(jié)點(diǎn)a處的注入功率將發(fā)生改變,進(jìn)而改變Padc。圖2為潮流方程右端項(xiàng)的并行計(jì)算過(guò)程。圖中,i(1)和j(1)分別為支路1對(duì)應(yīng)的流出和流入節(jié)點(diǎn);a(1)為發(fā)電機(jī)1對(duì)應(yīng)的節(jié)點(diǎn)編號(hào);(Pdci(1))1為場(chǎng)景1中的向量Pdc中節(jié)點(diǎn)i(1)對(duì)應(yīng)的元素,其余同理。當(dāng)線(xiàn)程編號(hào)索引小于m時(shí),由線(xiàn)程編號(hào)索引到斷線(xiàn)支路編號(hào),根據(jù)對(duì)應(yīng)的支路編號(hào)將等效功率注入列向量的第i個(gè)元素和第j個(gè)元素進(jìn)行修改;當(dāng)線(xiàn)程編號(hào)索引為m到m+s時(shí),由發(fā)電機(jī)編號(hào)索引對(duì)應(yīng)的節(jié)點(diǎn)編號(hào),對(duì)節(jié)點(diǎn)注入功率進(jìn)行修改。上述設(shè)計(jì)過(guò)程中不同的計(jì)算指令將分配到不同的線(xiàn)程束中,GPU并行處理此類(lèi)運(yùn)算,且并行度很高。
直流潮流模型平衡節(jié)點(diǎn)的相角為0,因此,對(duì)于2.5節(jié)中求得的電壓相角向量,需要在平衡節(jié)點(diǎn)的位置加零元素,其核函數(shù)設(shè)計(jì)思路如圖3所示。偽代碼見(jiàn)附錄A中表A2,令x為2.5節(jié)求得的電壓相角向量數(shù)組,新數(shù)組用x1表示。創(chuàng)建包含(m+s)×(n+1)個(gè)元素的一維數(shù)組,采用cudaMemset函數(shù)將其全部置0;然后創(chuàng)建(m+s)×n個(gè)線(xiàn)程,每個(gè)線(xiàn)程處理1個(gè)元素。CUDA異構(gòu)平臺(tái)采用單指令多線(xiàn)程架構(gòu)管理和執(zhí)行線(xiàn)程,每32個(gè)線(xiàn)程為1個(gè)線(xiàn)程束,線(xiàn)程束中所有線(xiàn)程同時(shí)執(zhí)行相同的指令,由于線(xiàn)程并不能根據(jù)批量處理規(guī)模完美分配,此處對(duì)線(xiàn)程編號(hào)設(shè)置判斷條件如下:編號(hào)小于(m+s)×n的線(xiàn)程執(zhí)行任務(wù),多余的線(xiàn)程被過(guò)濾。此時(shí),同1個(gè)線(xiàn)程束中依然執(zhí)行相同的指令,不會(huì)導(dǎo)致線(xiàn)程束的分支,保持很高的并行性。同時(shí)上述算法中內(nèi)存連續(xù)讀取,保證同一線(xiàn)程束下內(nèi)存的合并訪問(wèn),使算法具有高效性。
【參考文獻(xiàn)】:
期刊論文
[1]電力系統(tǒng)靜態(tài)安全最優(yōu)潮流并行計(jì)算方法[J]. 陽(yáng)育德,陶琢,劉輝,覃智君. 電力自動(dòng)化設(shè)備. 2019(01)
[2]一種GPU-CPU異構(gòu)運(yùn)算框架加速的實(shí)時(shí)N-1交流潮流計(jì)算方法[J]. 唐坤杰,董樹(shù)鋒,宋永華. 中國(guó)電機(jī)工程學(xué)報(bào). 2018(15)
[3]匹配GPU編程架構(gòu)的直流故障篩選算法設(shè)計(jì)[J]. 溫柏堅(jiān),謝恩彥,劉明波. 電力自動(dòng)化設(shè)備. 2017(05)
[4]基于網(wǎng)絡(luò)分塊算法的靜態(tài)安全快速計(jì)算方法[J]. 王方雨,劉文穎,田浩,梁琛,鄭晶晶. 電力自動(dòng)化設(shè)備. 2017(04)
[5]電力系統(tǒng)靜態(tài)安全混合控制方法[J]. 傅旭. 電力自動(dòng)化設(shè)備. 2017(01)
[6]基于補(bǔ)償法的預(yù)防性靜態(tài)安全控制[J]. 陽(yáng)育德,馮彥維,韋化. 電力自動(dòng)化設(shè)備. 2015(12)
[7]一種求解機(jī)組組合問(wèn)題的內(nèi)點(diǎn)半定規(guī)劃GPU并行算法[J]. 張寧宇,高山,趙欣. 電力自動(dòng)化設(shè)備. 2013(07)
[8]基于并行計(jì)算和數(shù)據(jù)復(fù)用的快速靜態(tài)安全校核技術(shù)[J]. 李峰,李虎成,於益軍,馮樹(shù)海. 電力系統(tǒng)自動(dòng)化. 2013(14)
本文編號(hào):3018019
本文鏈接:http://sikaile.net/kejilunwen/dianlilw/3018019.html
最近更新
教材專(zhuān)著