发明名称 一种基于直流潮流模型的支路开断灵敏度计算方法
摘要 本发明提供了一种源自直流法但保留了高阶项的支路开断灵敏度计算方法。该方法基于直流潮流模型,在分析传统的直流法求灵敏度的基础上提出了一种简便且计算结果精度更高的计算方法。该方法在求解灵敏度时,充分利用开断前后的状态信息,对高阶项进行了保留,在计算过程上并没有因此而增加计算量。而且,由于没有略去高阶项,误差来源只来源于所基于的模型上,计算精度得到了提高。由于计算简便且精度高,可以广泛应用于电力系统网络分析的各个应用领域,如EMS、DTS以及电力系统仿真和分析计算等领域,具有很好的推广应用价值。
申请公布号 CN102142682A 申请公布日期 2011.08.03
申请号 CN201010517291.3 申请日期 2010.10.18
申请人 中国电力科学研究院;东北电网有限公司 发明人 韩巍;李强;贾伟;孟令愚;易善军;杨宁;邵广惠;高德宾;刘家庆;李群;李忠满
分类号 H02J3/00(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 北京安博达知识产权代理有限公司 11271 代理人 徐国文
主权项 1.一种基于直流潮流模型的支路开断灵敏度计算方法,其特征在于包括以下步骤:(1)将电力系统物理节点及其连接关系通过对节点沿连接边搜索并归并,形成计算节点及其拓扑关系;电力系统中网络元件可以分为单端元件和双端元件,单端元件包括机组、负荷、电容器或电抗器,双端元件包括开关、线路和变压器,单端元件表征为一个物理节点,双端元件表征为两个物理节点及其连接它们的一条边,其关系可以表达为,图G=&lt;V,E&gt;,其中V={u<sub>1</sub>,……,u<sub>n</sub>},E={(u<sub>i</sub>,u<sub>j</sub>)|u<sub>i</sub>∈V,u<sub>j</sub>∈V,u<sub>i</sub>≠u<sub>j</sub>},通过依次对图G中节点V沿着边E搜索,将V中所有通过闭合开关相连的节点合并,形成新的图G<sub>1</sub>,其中G<sub>1</sub>=&lt;V<sub>1</sub>,E<sub>1</sub>&gt;,其中,V<sub>1</sub>表达计算点,E<sub>1</sub>表达计算点的连接关系;(2)依据计算节点及其拓扑连接关系生成电纳矩阵B<sub>0</sub>;其中n为计算节点数,其中无平衡节点对应行和对应列:<img file="FSA00000315560800011.GIF" wi="736" he="296" />B<sub>ij</sub>由下式定义,式中b<sub>ij</sub>为支路ij的电纳:<maths num="0001"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>B</mi><mi>ii</mi></msub><mo>=</mo><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow><mi>n</mi></munderover><msub><mi>b</mi><mi>ij</mi></msub></mtd></mtr><mtr><mtd><msub><mi>B</mi><mi>ij</mi></msub><mo>=</mo><msub><mi>b</mi><mi>ij</mi></msub></mtd></mtr></mtable></mfenced></math>]]></maths>(3)根据开断支路km及其支路的电纳值形成n-1维向量Y,其中n为计算节点数,Y<sub>i</sub>为向量第i个元素,b<sub>km</sub>为支路km的电纳值,如下所示:<maths num="0002"><![CDATA[<math><mrow><mi>Y</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>b</mi><mi>km</mi></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>-</mo><msub><mi>b</mi><mi>km</mi></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mfenced open='' close=''><mtable><mtr><mtd><mo>&LeftArrow;</mo><mi>k</mi></mtd></mtr><mtr><mtd></mtd></mtr><mtr><mtd></mtd></mtr><mtr><mtd></mtd></mtr><mtr><mtd><mo>&LeftArrow;</mo><mi>m</mi></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>(4)通过解方程X=B<sub>0</sub>Y,得到n-1维解向量X;(5)当计算支路ij对开断支路km的灵敏度时,将解向量X的第i、j、k、m处的向量元素代入如下公式,其中,b<sub>km</sub>为支路km的电纳值,b<sub>ij</sub>为支路ij的电纳值:<maths num="0003"><![CDATA[<math><mrow><msub><mi>SENS</mi><mrow><mi>ij</mi><mo>&RightArrow;</mo><mi>km</mi></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>b</mi><mi>ij</mi></msub><mrow><mo>(</mo><msub><mi>X</mi><mi>i</mi></msub><mo>-</mo><msub><mi>X</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>b</mi><mi>km</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>X</mi><mi>k</mi></msub><mo>+</mo><msub><mi>X</mi><mi>m</mi></msub><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths>
地址 100192 北京市海淀区清河小营东路15号