[發(fā)明專利]一種電力系統(tǒng)諧波檢測器部署方法

有效

• 申請?zhí)枺?/td>: 201710316267.5
• 申請日: 2017-05-08
• 公開（公告）號: CN106990287B
• 公開（公告）日: 2019-07-19
• 發(fā)明（設計）人: 王榮杰;周海峰;蔣德松
• 申請（專利權）人: 集美大學;中國水產科學研究院漁業(yè)機械儀器研究所
• 主分類號: G01R23/16
• 分類號: G01R23/16
• 代理公司: 廣州恒華智信知識產權代理事務所(普通合伙) 44299
• 代理人: 周長智

1.一種電力系統(tǒng)諧波檢測器部署方法，其特征在于，包括如下三大步驟：

步驟1：初始化；

步驟2：優(yōu)化尋優(yōu)；

步驟3：選出F(l)最小的θ(l)，如果θ(l)中的元素為1，則相應的母線部署諧波檢測器；如果θ(l)中的元素為0，則相應的母線不部署；其中：

步驟1分為如下步驟完成：

步驟1.1：給定部署諧波檢測器的期望冗余值R；

步驟1.2：初始化迭代次數k＝1，最大迭代次數k_max，k_limit，沒有得到更優(yōu)的解的記錄k_count(l)，蜂群的規(guī)模N_S和待優(yōu)化解的維數N；隨機產生二元優(yōu)化解θ(l,d)的初始值，l＝1,2,…,N_S；d＝1,2,…,N；θ(l)＝[θ(l,1),θ(l,2),…,θ(l,N)]；

步驟1.3：判斷θ(l)元素中1的個數是否超過n個，如果超過，則重新隨機產生θ(l)；

步驟1.4：根據θ(l)的元素判斷各條母線是否有部署諧波檢測器，在此基礎上再根據母線的諧波檢測器能否檢測到其它母線的諧波信息，構造矩陣X；然后由式(1)計算每個θ(l)的代價函數F(l)；

F(l)＝ω₁J₁+ω₂J₂ (1)

式(1)中，ω₁和ω₂為權值系數；J₁＝A^TA，J₂＝(R-XA)^T(R-XA)，J₁越小表示部署的諧波檢測器數量越少，J₂越小表示部署的諧波檢測器所能獲得信息冗余度越接近期望冗余值R；

其中，步驟2分為如下步驟完成：

步驟2.1：工蜂尋優(yōu)階段，分為以下具體步驟：

步驟2.1.1：將l設置為1；

步驟2.1.2：如果l≤N_S,將執(zhí)行步驟2.1.3，否則跳至步驟2.2；

步驟2.1.3：

隨機產生l相鄰的序號i和j，i≠j≠l；

選出F(l)最小的θ(l)，記為θ_best；

產生一個N維的零列向量，記為zero；

對θ(l)、θ(i)、θ(j)、θ_best和zero兩兩進行隨機配對；

對配對后的θ(l)、θ(i)、θ(j)、θ_best和zero兩兩互換它們的部分元素；

由式(1)計算θ(l)、θ(i)、θ(j)、θ_best和zero的代價函數，并從中選出代價函數最小的，作為θ_EB(l)及相應地賦值予F_EB(l)；

如果F_EB(l)<F(l)，則將F_EB(l)賦值予F(l)，θ_EB(l)賦值予θ(l)；否則，k_count(l)＝k_count(l)+1；跳至步驟2.1.2；

步驟2.2：根據式(2)計算每個θ(l)相應代價函數的概率；

步驟2.3：工蜂尋優(yōu)階段，分為以下具體步驟：

步驟2.3.1：將l設置為1；

步驟2.3.2：如果l≤N_S,將執(zhí)行步驟2.3.3，否則跳至步驟2.4；

步驟2.3.3：如果p(l)<rand[0 1],則執(zhí)行步驟2.3.4，否則跳至步驟2.4；

步驟2.3.4：

隨機產生l相鄰的序號i和j，i≠j≠l；

由式(3)計算θ_OB(l)；

θ_OB(l)＝θ(l)+rand[0 1][θ(i)-θ(j)] (3)

如果θ_OB(l)中的元素等于或大于0.5，則把它置為1；否則，置為0；

對θ(i)和θ(j)隨機地互換它們的部分元素；

由式(1)計算θ_OB(l)、θ(i)和θ(j)的代價函數，并從中選出代價函數最小的，更新θ_OB(l)及其F_OB(l)；

如果F_OB(l)<F(l)，則將F_OB(l)賦值予F(l)，θ_OB(l)賦值予θ(l)；否則，k_count(l)＝k_count(l)+1；跳至步驟2.3.2；

步驟2.4：偵察蜂尋優(yōu)階段，分為以下具體步驟：

步驟2.4.1：將l設置為1；

步驟2.4.2：如果l≤N_S,將執(zhí)行步驟2.4.3，否則跳至步驟2.5；

步驟2.4.3：如果k_count(l)≥k_limit,則執(zhí)行步驟2.4.4，否則跳至步驟2.5；

步驟2.4.4：

隨機產生一個二元的N維列向量，記為θ_SB；如果θ_SB元素中1的個數超過n個，則重新隨機產生；

由式(1)計算θ_SB的代價函數，記為F_SB；

如果F_SB<F(l)，則將F_SB賦值予F(l)，θ_SB賦值予θ(l)，k_count(l)＝0；跳至步驟2.5；

步驟2.5：k＝k+1；

步驟2.6：如果k大于k_max，則跳至步驟3；否則，跳至步驟2.1.3。