发明名称 一种基于导纳矩阵快速修正的电网潮流计算方法
摘要 本发明涉及一种基于导纳矩阵快速修正的电网潮流计算方法,利用电网现有广域测量系统实现电网各变电站电压、电流相量的同步测量,获取电网同步数据;若发生电网单一线路跳闸,通过低维数矩阵求逆计算中间矩阵和变量矩阵,进而计算线路跳闸后电网各变电站的电压变化量,计算线路跳闸后电网各条线路的电流变化量,然后利用计算得到电网各条线路电流变化量与线路跳闸前电网各条线路的电流进行叠加得到线路跳闸后电网各条线路流过的电流。本发明方法计算速度快,整个计算过程可在200ms时间内完成,能为电力系统稳定分析控制提供实时电网基础数据,增强电网安全运行的稳定性。
申请公布号 CN104065066A 申请公布日期 2014.09.24
申请号 CN201410339280.9 申请日期 2014.07.16
申请人 国家电网公司;国网福建省电力有限公司;国网福建省电力有限公司检修分公司 发明人 曾惠敏
分类号 H02J3/00(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 厦门市首创君合专利事务所有限公司 35204 代理人 张松亭
主权项 一种基于导纳矩阵快速修正的电网潮流计算方法,其特征在于,包括如下依序步骤:1)利用电网现有广域测量系统实现电网各变电站电压、电流相量的同步测量,实时获取电网同步数据,若线路ij两端的断路器跳开后,进入步骤2);其中,线路ij连接i变电站和j变电站;i=1,…,n;j=1,…,n;n为电网变电站总数目;2)计算线路ij两端的断路器跳开后,电网各变电站的电压变化量矩阵<img file="FDA0000538642410000018.GIF" wi="101" he="66" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>&Delta;</mi><mover><mi>U</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mo>-</mo><msup><mi>Y</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>M</mi><msup><mrow><mo>(</mo><mi>E</mi><mo>+</mo><msup><mi>Y</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>M</mi><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>J</mi><mover><mi>U</mi><mo>&CenterDot;</mo></mover></mrow>]]></math><img file="FDA0000538642410000011.GIF" wi="592" he="100" /></maths>其中,Y为电网节点导纳矩阵;M=[0 … 1 … ‑1 … 0]<sup>T</sup>为n×1维矩阵,M矩阵中第i行第1列元素为1,M矩阵中第j行第1列元素为‑1,其他元素为0;i=1,…,n;j=1,…,n;n为电网变电站总数目;E为n×n阶单位矩阵;J=[0 … ‑y<sub>ij</sub> … y<sub>ij</sub> … 0]为1×n维矩阵,J矩阵中第1行第i列元素为‑y<sub>ij</sub>,J矩阵中第1行第j列元素为y<sub>ij</sub>,其他元素为0;<img file="FDA0000538642410000012.GIF" wi="54" he="67" />为线路ij两端的断路器跳开前,广域测量系统测量到的电网各变电站电压相量矩阵;y<sub>ij</sub>为线路ij的线路导纳;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>&Delta;</mi><mover><mi>U</mi><mo>&CenterDot;</mo></mover><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&Delta;</mi><msub><mover><mi>U</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub></mtd><mtd><mi>&Delta;</mi><msub><mover><mi>U</mi><mo>&CenterDot;</mo></mover><mn>2</mn></msub></mtd><mtd><mi>&Delta;</mi><msub><mover><mi>U</mi><mo>&CenterDot;</mo></mover><mn>3</mn></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mi>&Delta;</mi><msub><mover><mi>U</mi><mo>&CenterDot;</mo></mover><mrow><mi>n</mi><mo>-</mo><mn>2</mn></mrow></msub></mtd><mtd><mi>&Delta;</mi><msub><mover><mi>U</mi><mo>&CenterDot;</mo></mover><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><mi>&Delta;</mi><msub><mover><mi>U</mi><mo>&CenterDot;</mo></mover><mi>n</mi></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup></mrow>]]></math><img file="FDA0000538642410000013.GIF" wi="956" he="86" /></maths>为n×1维矩阵,<img file="FDA0000538642410000014.GIF" wi="540" he="71" /><img file="FDA0000538642410000015.GIF" wi="238" he="71" />分别为线路ij两端的断路器跳开后,编号为第1、2、3、…、n‑2、n‑1、n变电站的电压变化量;3)计算线路ij两端的断路器跳开后线路lk的电流<img file="FDA0000538642410000016.GIF" wi="424" he="134" />其中,n为电网变电站总数目;l=1,…,n;k=1,…,n;z<sub>lk</sub>为线路lk的线路阻抗;<img file="FDA0000538642410000017.GIF" wi="60" he="79" />为线路ij两端的断路器跳开前,广域测量系统测量到线路lk上流过的电流;H=[0 … 1 … ‑1 … 0]<sup>T</sup>为n×1维矩阵,H矩阵中第l行第1列元素为1,H矩阵中第k行第1列元素为‑1,其他元素都为0;线路lk连接l变电站和k变电站;线路lk与线路ij为不同线路。
地址 100031 北京市西城区西长安街86号