[發明專利]一種電力系統廣域測量系統通信時延的實時仿真方法

有效

• 申請號: 201210434116.7
• 申請日: 2012-11-01
• 公開（公告）號: CN102970180A
• 公開（公告）日: 2013-03-13
• 發明（設計）人: 黎雄;孫元章;張放;程林;劉志雄;曾敏;王舉;徐箭;林成;趙維興;唐建興;王國松;馬覃峰
• 申請（專利權）人: 武漢大學;貴州電網公司電力調度控制中心
• 主分類號: H04L12/26
• 分類號: H04L12/26;H04L12/24;H02J13/00
• 代理公司: 武漢科皓知識產權代理事務所(特殊普通合伙) 42222
• 代理人: 嚴彥

1.一種電力系統廣域測量系統通信時延的實時仿真方法，其特征在于：定義變量j的初始值為0；設關注的WAMS量測有K個，對應有K組中間輔助變量，每組中間輔助變量包括相應的測量值緩沖數據序列、通信時延緩沖數據序列、時延數據序列和WAMS數據鏈；所述測量值緩沖數據序列是用于存放當前仿真時刻起至之前最大仿真時延dT_max段內的測量值的序列，所述通信時延緩沖數據序列是用于存放與測量值緩沖數據序列中量測值的對應時延數據的序列，所述時延數據序列是用于存放待仿真的時延數據的序列，所述WAMS數據鏈是根據時標信息進行同步的帶時標WAMS測量數據的序列；

每當完成一個仿真步長的仿真計算后，得到多個量測的瞬時值，開始當前仿真步長的通信時延仿真，通信時延仿真包括以下步驟，

步驟1，令j＝j+1；

步驟2，根據限定處理時延的最大循環次數N值對關注的K個WAMS量測分別執行小步長子流程，所述小步長子流程以實時仿真系統的仿真步長dt實時輸出量測受通信時延影響的值；

其中，dT為WAMS測量步長，仿真的最大時延為dT_max；

步驟3，輸出經過小步長子流程處理后的K個WAMS量測的受通信網絡時延影響的值；步驟4，判斷是否到達一次WAMS測量周期，判斷條件為j等于M，如果不相等，則沒有到達WAMS測量周期，進行下一次仿真，再次執行步驟1；如果相等，則到達一次WAMS測量周期，繼續執行步驟5；

其中，dt為實時仿真系統的仿真步長，dT為WAMS測量步長；

步驟5，根據限定處理時延的最大循環次數N值對關注的K個WAMS量測分別執行大步長子流程，所述大步長子流程用于以WAMS測量步長dT輸出量測受通信時延影響的值以及相應的時標；

步驟6，根據經過大步長子流程處理后的K個WAMS量測的對應的WAMS數據鏈，協調處理K個受時延影響的量測數據；

步驟7，結束當前仿真步長的通信時延仿真。

2.根據權利要求1所述電力系統廣域測量系統通信時延的實時仿真方法，其特征在于：步驟2中，將關注的WAMS量測依次標記為第k個量測，k=1，...,K；

對第k個量測執行的小步長子流程如下：

步驟a，初始化變量i=0；

步驟b，判斷i是否小于N，若i大于等于N，則結束小步長子流程，進入步驟3；若i小于N則繼續進行步驟c；

步驟c，逐一判斷通信時延緩沖數據序列中的第i個受時延影響的測量值是否能用；具體方法為，設通信時延緩沖數據序列為td_buf，從通信時延緩沖數據序列td_buf中讀取第k個量測的第i個通信時延緩沖數據td_buf(k)(i)，記為t，表達式為，

t=td_buf(k)(i)

然后，判斷t的值是否與預設的特殊網絡狀態對應的數值之一相等，若沒有相等的情況，則執行步驟d；若有相等的情況，則進行針對該特殊狀態的相應處理，i=i+1，并返回步驟a；

步驟d，判斷第i個受時延影響的測量值是否能用，包括先用t值減相應的時延，表達式為，

t=t-i×dT-j×dt

若t>0，則該時刻受時延影響的測量值不能用，i=i+1，并返回步驟b；若t<0，則執行步驟e；

步驟e，計算第k個量測受時延影響后的實時測量值yk_out為，

yk_out=y_buf(k)(i)

其中，y_buf為測量值緩沖數據序列，y_buf(k)(i)為第k個量測的第i個測量值緩沖數據；

步驟f，結束小步長子流程。

3.根據權利要求1或2所述電力系統廣域測量系統通信時延的實時仿真方法，其特征在于：步驟5中，將關注的WAMS量測依次標記為第k個量測，k=1，...,K；

對第k個量測執行的大步長子流程如下：

步驟a，初始化變量i=0；

步驟b，判斷i是否小于N，若i大于等于N，則結束大步長子流程，進入步驟6；若i小于N則繼續進行步驟c；

步驟c，逐一判斷通信時延緩沖數據序列中的第i個受時延影響的測量值是否能用；具體方法為，從通信時延緩沖數據序列td_buf中讀取第k個量測的第i個通信時延緩沖數據td_buf(k)(i)，記為t，表達式為，

t=td_buf(k)(i)

然后，判斷t的值是否與預設的特殊網絡狀態對應的數值之一相等，若沒有相等的情況，則執行步驟d；若有相等的情況，則進行針對該特殊狀態的相應處理，則i=i+1，并返回步驟a；

步驟d，判斷第i個受時延影響的測量值是否能用，包括先用t值減相應的時延，表達式為，

t=t-i×dT

若t>0，則該時刻受時延影響的測量值不能用，i=i+1，并返回步驟b；若t<0，則執行步驟e；

步驟e，計算第k個量測受時延影響后，用于保存的測量值yk_d為，

yk_d=y_buf(k)(i)

其中，y_buf為測量值緩沖數據序列，y_buf(k)(i)為第k個量測的第i個測量值緩沖數據；

步驟f，更新測量值緩沖數據序列、通信時延緩沖數據序列如下，

y_buf(k)(M-1)=y_buf(k)(M-2);

...

y_buf(k)(2)=y_buf(k)(1);

y_buf(k)(1)=y_buf(k)(0);

y_buf(k)(0)=yk;

及

td_buf(k)(M-1)=td_buf(k)(M-2);

...

td_buf(k)(2)=td_buf(k)(1);

td_buf(k)(1)=td_buf(k)(0);

td_buf(k)(0)=td_new(k);

其中，yk為實時仿真系統中第k個量測的實時量測值，td_new(k)為依次從自第k個量測的相應時延數據序列中取的一個數據；

步驟g，結束大步長子流程。