发明名称 智能电网广域负荷建模方法
摘要 本发明公开了一种智能电网广域负荷建模方法,将安装到户的智能电表按电压等级建立分层、分级的广域负荷建模体系构架,通过自动抄表系统将智能电表测量的功率信息和辨识出的负荷相关信息分层收集并集中起来,利用负荷建模算法对收集得到的负荷进行处理,建立起准确的负荷模型:利用各电压等级以及安装到户的智能电表建立广域负荷建模体系,通过自动抄表系统将电表的功率或负荷辨识结果等相关信息收集起来,综合利用已有的负荷建模算法对这些信息进行处理,实现电力系统负荷的广域建模,基于智能电表可实现负荷建模的长效机制以及覆盖最广大的范围。
申请公布号 CN102868160B 申请公布日期 2015.03.11
申请号 CN201210357915.9 申请日期 2012.09.24
申请人 南京南瑞继保电气有限公司;南京南瑞继保工程技术有限公司 发明人 王亮;李力;侯学勇;吕航;文继锋;陈松林
分类号 H02J3/00(2006.01)I;G06F17/50(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 代理人
主权项 智能电网广域负荷建模方法,其特征是步骤如下:将安装到户的智能电表按电压等级建立分层、分级的广域负荷建模体系构架,通过自动抄表系统将智能电表测量的功率信息和辨识出的负荷相关信息分层收集并集中起来,利用负荷建模算法对收集得到的负荷进行处理,建立起准确的负荷模型:1)广域负荷建模的系统构架:电力系统分为工业负荷、城市民用负荷、商业负荷、农业负荷及其它负荷,利用智能电表进行建模的工作需要按照负荷分区、分类的原则,有针对性的按电压等级逐步开展;2)广域负荷建模的方法:将统计综合法和总体测辨法二者相结合,形成基于智能电表输出数据的广域建模方法:首先,采用在线统计综合法,充分利用智能电表提供的功率数据,按电压等级以及电表的归属关系逐级向上统计,获得负荷模型和参数;同时,利用智能电表的电压、电流采样数据采用在线测辨法获得各电压等级负荷的模型和参数,最后,二种方法获得的辨识结果相互验证;具体建模策略如下:2‑1)在线测辨,智能电表作为在线测辨装置,对输入的电压、电流进行数据采集,然后确定负荷模型结构,最后根据实时采集的数据辨识出模型参数;对于负荷成分比较复杂的情况,采用总体测辨方法,无需知道用户的负荷组成及参数,获得较好的参数估计值,由于长期在线测量,根据各个时刻测量数据得到相应的负荷特性参数,从而较好的解决了负荷特性的时变特性,无需专门安装大量的在线测辨装置;2‑2)横向对比,对于同一电压等级的同类负荷得到的辨识结果要进行横向对比,对于明显超出误差范围的坏数据辨识应剔除,验证本层级模型结构、参数的合理性,若不合理,还将修正模型结构,再进行辨识;处于末端的智能电表要重点关注负荷的成分、比重以及时变特性而输出数据;2‑3)统计综合,充分利用智能电表进行功率统计,逐级向上实现负荷建模的在线统计综合;由于智能电表定时刷新的用电情况更准确、更广泛,因此能更好地实现负荷模型参数的动态、长效管理:负荷辨识算法的移植:采用“先动后静”的辨识方法,同时获得动态和静态参数;首先,采用动态模型去描述负荷;其次在电压及频率与初值偏离很小时,将静态模型线性化,动态模型的公式如下:<img file="FDA0000567241140000021.GIF" wi="928" he="170" /><img file="FDA0000567241140000022.GIF" wi="1572" he="193" />即:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&Delta;P</mi></mtd></mtr><mtr><mtd><mi>&Delta;Q</mi></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>dP</mi><mo>/</mo><mi>dU</mi></mtd><mtd><mi>dP</mi><mo>/</mo><mi>df</mi></mtd></mtr><mtr><mtd><mi>dQ</mi><mo>/</mo><mi>dU</mi></mtd><mtd><mi>dQ</mi><mo>/</mo><mi>df</mi></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&Delta;U</mi></mtd></mtr><mtr><mtd><mi>&Delta;f</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000567241140000023.GIF" wi="1563" he="158" /></maths>其中ΔU,Δf为电压和频率是输入变量,ΔP,ΔQ为有功和无功是输出变量;其次,在稳态条件下即电压与频率在小范围变化时动态模型退化为静态模型,其状态方程为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfrac><mi>dX</mi><mi>dt</mi></mfrac><mo>=</mo><mi>AX</mi><mo>+</mo><mi>BV</mi></mrow>]]></math><img file="FDA0000567241140000024.GIF" wi="354" he="126" /></maths>Y=CX+DV                                (3)其中:X:状态矢量Y:输出状态矢量V:控制矢量A:状态系数矩阵B:控制系数矩阵C:输出状态系数矩阵D:输出控制系数矩阵根据负荷的物理本质可知,动态模型中各变量导数为零即得静态模型,令:<img file="FDA0000567241140000025.GIF" wi="215" he="141" />可得:Y=[‑CA<sup>‑1</sup>B+D]V                             (4)其中:Y=[ΔP,ΔQ]<sup>T</sup>,V=[ΔU,Δf]<sup>T</sup>比较式(2)与(4)可得:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><mi>dP</mi><mo>/</mo><mi>dU</mi></mtd><mtd><mi>dP</mi><mo>/</mo><mi>df</mi></mtd></mtr><mtr><mtd><mi>dQ</mi><mo>/</mo><mi>dU</mi></mtd><mtd><mi>dQ</mi><mo>/</mo><mi>df</mi></mtd></mtr></mtable></mfenced><mo>=</mo><msup><mrow><mo>-</mo><mi>CA</mi></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>B</mi><mo>+</mo><mi>D</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000567241140000031.GIF" wi="1285" he="168" /></maths>而静态特征系数为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>p</mi><mi>u</mi></msub><mo>=</mo><mrow><mo>(</mo><mi>dP</mi><mo>/</mo><mi>dU</mi><mo>)</mo></mrow><mfrac><msub><mi>U</mi><mn>0</mn></msub><msub><mi>P</mi><mn>0</mn></msub></mfrac><mo>,</mo><msub><mi>p</mi><mi>f</mi></msub><mo>=</mo><mrow><mo>(</mo><mi>dP</mi><mo>/</mo><mi>df</mi><mo>)</mo></mrow><mfrac><msub><mi>f</mi><mn>0</mn></msub><msub><mi>P</mi><mn>0</mn></msub></mfrac></mrow>]]></math><img file="FDA0000567241140000032.GIF" wi="783" he="147" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>q</mi><mi>u</mi></msub><mo>=</mo><mrow><mo>(</mo><mi>dQ</mi><mo>/</mo><mi>dU</mi><mo>)</mo></mrow><mfrac><msub><mi>U</mi><mn>0</mn></msub><msub><mi>Q</mi><mn>0</mn></msub></mfrac><mo>,</mo><msub><mi>q</mi><mi>f</mi></msub><mo>=</mo><mrow><mo>(</mo><mi>dQ</mi><mo>/</mo><mi>df</mi><mo>)</mo></mrow><mfrac><msub><mi>f</mi><mn>0</mn></msub><msub><mi>Q</mi><mn>0</mn></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000567241140000033.GIF" wi="1389" he="159" /></maths>其中p<sub>u</sub>,p<sub>f</sub>,q<sub>u</sub>,q<sub>f</sub>为负荷静态特征系数,而dP/dU,dP/df,dQ/dU,dQ/df由动态模型求得,这样就从动态模型中得到了负荷的静态特征系数。
地址 211102 江苏省南京市江宁区苏源大道69号