[發(fā)明專利]基于受擾軌跡動(dòng)態(tài)特性的數(shù)值積分提前終止方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201811277682.5
• 申請(qǐng)日: 2018-10-30
• 公開（公告）號(hào): CN109492286B
• 公開（公告）日: 2022-07-01
• 發(fā)明（設(shè)計(jì)）人: 薛禹勝;黃天罡;薛峰
• 申請(qǐng)（專利權(quán)）人: 南瑞集團(tuán)有限公司;國(guó)電南瑞科技股份有限公司;國(guó)網(wǎng)青海省電力公司;國(guó)家電網(wǎng)有限公司
• 主分類號(hào): G06F30/20
• 分類號(hào): G06F30/20
• 代理公司: 南京縱橫知識(shí)產(chǎn)權(quán)代理有限公司 32224
• 代理人: 董建林;范青青

1.基于受擾軌跡動(dòng)態(tài)特性的數(shù)值積分提前終止方法，其特征在于，所述方法包括：根據(jù)受擾軌跡信息，結(jié)合包括：失穩(wěn)模式易變度、受擾軌跡非同調(diào)度、最大轉(zhuǎn)子角間隙、同向相鄰擺次轉(zhuǎn)子角間隙信息相似度及振蕩衰減度在內(nèi)的動(dòng)態(tài)特性，判斷是否提前終止數(shù)值積分；

所述方法具體包括如下步驟：

對(duì)測(cè)試算例全集中的每個(gè)算例分別執(zhí)行下述步驟：

步驟A：求取算例的失穩(wěn)模式易變度K；

步驟B：采用小步長(zhǎng)逐步積分方法向前積分一步，直至到達(dá)主導(dǎo)映象首擺最遠(yuǎn)點(diǎn)或動(dòng)態(tài)鞍點(diǎn)；

步驟C：若采用小步長(zhǎng)逐步積分已達(dá)主導(dǎo)映象首擺最遠(yuǎn)點(diǎn)，則根據(jù)下述判定條件判定是否提前終止積分：

1)失穩(wěn)模式易變度K值小于第一閾值ε₁；

2)主導(dǎo)映象機(jī)械功率P_m大于電磁功率分量P_c；

若同時(shí)滿足上述兩個(gè)判定條件則提前終止積分；否則，執(zhí)行下一步；

步驟D：若采用小步長(zhǎng)逐步積分已達(dá)主導(dǎo)映象首擺動(dòng)態(tài)鞍點(diǎn)，則根據(jù)下述判定條件判定是否提前終止積分：

1)失穩(wěn)模式易變度K值小于第一閾值ε₁；

若滿足上述判定條件則提前終止積分；否則，執(zhí)行下一步；

步驟E：基于首擺暫態(tài)過程受擾軌跡，求取首擺非同調(diào)度μ₁；

步驟F：若采用小步長(zhǎng)逐步積分已達(dá)主導(dǎo)映象首擺最遠(yuǎn)點(diǎn)，則根據(jù)下述判定條件判定是否提前終止積分：

1)失穩(wěn)模式易變度K值小于第二閾值ε₂；

2)主導(dǎo)映象機(jī)械功率P_m大于電磁功率分量P_c；

3)首擺非同調(diào)度μ₁值小于第三閾值ε₃；

若同時(shí)滿足上述三個(gè)判定條件則提前終止積分；否則，執(zhí)行下一步；

步驟G：若采用小步長(zhǎng)逐步積分已達(dá)主導(dǎo)映象首擺動(dòng)態(tài)鞍點(diǎn)，則根據(jù)下述判定條件判定是否提前終止積分：

1)失穩(wěn)模式易變度K值小于第二閾值ε₂；

2)首擺非同調(diào)度μ₁值小于第三閾值ε₃；

若同時(shí)滿足上述兩個(gè)判定條件則提前終止積分；否則，執(zhí)行下一步；

步驟H：采用小步長(zhǎng)逐步積分方法向前積分一步，直至最大轉(zhuǎn)子角間隙大于預(yù)設(shè)的第四閾值ε₄，或到達(dá)主導(dǎo)映象第三擺最遠(yuǎn)點(diǎn)或動(dòng)態(tài)鞍點(diǎn)；

步驟I：若最大轉(zhuǎn)子角間隙已大于第四閾值ε₄，則提前終止積分；否則，執(zhí)行下一步；其中：δ_ij表示機(jī)組i轉(zhuǎn)子角同機(jī)組j轉(zhuǎn)子角的差值，即：δ_ij＝δ_i-δ_j；i表示領(lǐng)前群中第i個(gè)機(jī)組；j表示余下群中第j個(gè)機(jī)組；S表示領(lǐng)前群；A表示余下群；

步驟J：若采用小步長(zhǎng)逐步積分已達(dá)主導(dǎo)映象第三擺最遠(yuǎn)點(diǎn)或動(dòng)態(tài)鞍點(diǎn)，求取前三擺非同調(diào)度μ_1-3；若前三擺非同調(diào)度μ_1-3值小于第五閾值ε₅，則提前終止積分；否則，執(zhí)行下一步；

步驟K：采用小步長(zhǎng)逐步積分方法向前積分一步，直至最大轉(zhuǎn)子角間隙大于預(yù)設(shè)的第四閾值ε₄，或到達(dá)主導(dǎo)映象下一擺最遠(yuǎn)點(diǎn)；

步驟L：若最大轉(zhuǎn)子角間隙已大于第四閾值ε₄，則提前終止積分；否則，執(zhí)行下一步；

步驟M：若采用小步長(zhǎng)逐步積分已達(dá)主導(dǎo)映象下一擺最遠(yuǎn)點(diǎn)，從該擺次始，向后溯洄四個(gè)擺次，若已小于第四擺，執(zhí)行步驟K，否則執(zhí)行下一步；

步驟N：根據(jù)步驟M中向后溯洄四個(gè)擺次的暫態(tài)過程，求取正向及反向相鄰擺次轉(zhuǎn)子角間隙信息相似度，若相似度小于第六閾值ε₆，則提前終止積分；否則，執(zhí)行下一步；

步驟O：若從該擺次始，向后溯洄六個(gè)擺次，若已小于第四擺，則執(zhí)行步驟K，否則執(zhí)行下一步；

步驟P：根據(jù)步驟O中向后溯洄的六個(gè)擺次的暫態(tài)過程，求取連續(xù)兩個(gè)正向及反向相鄰擺次振蕩衰減度，若振蕩衰減度小于第七閾值ε₇，則提前終止積分；否則，執(zhí)行步驟K。

2.根據(jù)權(quán)利要求1所述的基于受擾軌跡動(dòng)態(tài)特性的數(shù)值積分提前終止方法，其特征在于，所述失穩(wěn)模式易變度K采用大步長(zhǎng)泰勒級(jí)數(shù)展開快速求取。