[發(fā)明專利]一種包含隨機功率單元相關(guān)性問題的隨機潮流計算方法

有效

• 申請?zhí)枺?/td>: 201510040789.8
• 申請日: 2015-01-27
• 公開（公告）號: CN104699950B
• 公開（公告）日: 2018-02-02
• 發(fā)明（設(shè)計）人: 楊曉東;唐田;林章歲;李喜蘭;蔡霽霖;徐青山
• 申請（專利權(quán)）人: 國家電網(wǎng)公司;國網(wǎng)福建省電力有限公司;國網(wǎng)福建省電力有限公司經(jīng)濟技術(shù)研究院;東南大學(xué)
• 主分類號: G06F19/00
• 分類號: G06F19/00
• 代理公司: 福州元創(chuàng)專利商標(biāo)代理有限公司35100
• 代理人: 蔡學(xué)俊

1.一種包含隨機功率單元相關(guān)性問題的隨機潮流計算方法，其特征在于包括如下步驟：

步驟S1；對系統(tǒng)包含的所有隨機功率單元按照相關(guān)性和地理位置進行分組，每個隨機功率單元組視作一個隨機變量，記錄所述隨機變量的個數(shù)為n；

步驟S2；根據(jù)隨機功率單元的分組情況，確定簡化后的系統(tǒng)網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)和蒙特卡洛法的仿真次數(shù)m；

步驟S3；獲取步驟S1得到的各隨機變量在仿真周期內(nèi)的歷史時序數(shù)據(jù)，統(tǒng)計得到其累積分布函數(shù)；

步驟S4；利用改進聯(lián)合正態(tài)變換法，即JNT法，得到各隨機變量的采樣值，并形成一個n維隨機向量，重復(fù)m次；

步驟S5；利用蒙特卡洛法，利用m個n維隨機向量作為系統(tǒng)輸入，分別求解電力系統(tǒng)潮流方程，獲得經(jīng)典潮流問題的解的集合；

所述步驟S4具體包括以下步驟：

步驟S41：計算n個隨機變量x₁至x_n之間的等級相關(guān)系數(shù)矩陣R_r，所述的R_r為n維矩陣，并根據(jù)公式將R_r轉(zhuǎn)化為積差相關(guān)系數(shù)矩陣R，

其中r_r為矩陣R_r中的任一元素，r是矩陣R中與r_r處于對應(yīng)位置的元素；

步驟S42：由于所述步驟S41得到的R是正定矩陣，將R矩陣分解為R＝AA^T的形式；

步驟S43：對服從標(biāo)準(zhǔn)正態(tài)分布的任意一維隨機變量進行n次獨立的采樣，構(gòu)成各分量相互獨立的總體采樣向量并將所述通過進行正交變換得到所需的正態(tài)邊際分布域中的n維隨機向量

步驟S44：利用正態(tài)分布的分布函數(shù)Φ將轉(zhuǎn)換到均勻邊際分布域：

U→=Φ(N→)=(u1,u2,...,un)]]>

Φ函數(shù)采用離散數(shù)據(jù)對或離散數(shù)據(jù)點來保存，兩個保存的相鄰離散點之間的數(shù)據(jù)點采用線性插值的方法來近似獲取；

步驟S45：利用x₁至x_n各自的累積分布函數(shù)的反函數(shù)對所述中的各個對應(yīng)分量進行變換，得到實際邊際分布域中的隨機變量x₁至x_n的采樣值。

2.根據(jù)權(quán)利要求1所述的一種包含隨機功率單元相關(guān)性問題的隨機潮流計算方法，其特征在于：所述步驟S1具體為：所述的隨機功率單元分為隨機處理單元和負(fù)荷單元，將地理位置接近并且具有強相關(guān)性的隨機出力單元和負(fù)荷單元分別分組，將分好的每一個隨機功率單元組視作一個隨機變量，并假設(shè)各組內(nèi)的功率單元均為完全正相關(guān)，根據(jù)分組情況確定對應(yīng)的隨機出力單元變量和負(fù)荷變量。

3.根據(jù)權(quán)利要求1所述的一種包含隨機功率單元相關(guān)性問題的隨機潮流計算方法，其特征在于：所述步驟S3包括如下具體步驟：

步驟S31：將各隨機變量記為x₁至x_n，其歷史時序數(shù)據(jù)按時間順序編號，形成形式為(編號,功率數(shù)值)的數(shù)據(jù)對；

步驟S32：通過對步驟S31中所述數(shù)據(jù)對中的“編號”進行均勻分布采樣，獲取數(shù)據(jù)對的采樣點；

步驟S33：將步驟S32中采樣獲得的數(shù)據(jù)對按照“功率數(shù)值”從小到大排列，統(tǒng)計不同的“功率數(shù)值”的采樣值個數(shù)，并進行累加和歸一化，用以形成各隨機變量x₁至x_n的累積分布函數(shù)。