发明名称 一种基于变电站量测的负荷模型辨识方法
摘要 本发明公开了一种电力系统的负荷模型建模方法,首先对电网中的变电站进行同调机群划分,构建等值感应电动机的负荷模型,建立辨识准则函数,对负荷模型中待辨识的负荷参数通过输入电压和实际系统的输出响应采用粒子群算法进行寻优,统计不同站点的负荷模型,得到地区电网负荷模型参数库。本发明在辨识参数中引入初始有功功率比例系数K<sub>pm</sub>和额定初始负荷率系数M<sub>lf</sub>,消除负荷幅值的时变性带来的影响,使负荷模型更加准确;并设置了粒子群算法的速度调整因子和惯性因子,提高算法收敛精度。
申请公布号 CN103872678A 申请公布日期 2014.06.18
申请号 CN201410079486.2 申请日期 2014.03.06
申请人 国家电网公司;四川电力科学研究院;电子科技大学 发明人 丁理杰;王均;黄琦;魏巍;刘影;张周晶;郭巍;朱超
分类号 H02J3/00(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 成都行之专利代理事务所(普通合伙) 51220 代理人 温利平
主权项 1.一种基于变电站量测的负荷模型辨识方法,其特征在于,包括以下步骤:S1:对电网中的变电站根据电压等级进行划分,将每个变电站的等值负荷模型采用感应电动机并联静态负荷来表示;S2:等值感应电动机采用三阶感应电动机感应电动机负荷模型,α=[R<sub>s</sub>,X<sub>S</sub>,X<sub>M</sub>,R<sub>R</sub>,X<sub>R</sub>,H,A,B,k<sub>Pz</sub>,k<sub>Pi</sub>,k<sub>Qz</sub>,k<sub>Qi</sub>,K<sub>pm</sub>,M<sub>lf</sub>]<sup>T</sup>为负荷模型的独立待辨识参数向量,β=[e<sub>x0</sub>,e<sub>y0</sub>,s<sub>0</sub>,x<sub>s</sub>,K<sub>l</sub>]<sup>T</sup>为可以通过电机稳态条件和α求得的辨识参数向量,其中R<sub>S</sub>为定子电阻、X<sub>S</sub>为定子绕阻漏抗、X<sub>M</sub>为励磁电抗、R<sub>R</sub>为转子电阻、X<sub>R</sub>为转子漏抗、H为转子惯性时间常数、A、B为电机机械转矩特性参数,k<sub>Pz</sub>、k<sub>Pi</sub>、k<sub>Qz</sub>、k<sub>Qi</sub>为静态特性参数,e<sub>x0</sub>,e<sub>y0</sub>为电动机初始的暂态电压,s<sub>0</sub>为电动机的转差率,x<sub>s</sub>为定子和转子之间的同步电抗,K<sub>l</sub>为负载率;K<sub>pm</sub>表示初始有功功率比例系数,K<sub>pm</sub>=P<sub>0</sub>′/P,P<sub>0</sub>′为电动机初始有功功率,P为所测负荷点在暂态过程中所消耗的有功功率;M<sub>lf</sub>为额定初始负荷率系数,<img file="FDA0000473353530000011.GIF" wi="428" he="151" />S<sub>MB</sub>是感应电动机的额定容量,U<sub>B</sub>是负荷基准电压,U<sub>0</sub>是暂态过程中负荷母线电压初始值;S3:进行负荷模型参数辨识,具体步骤包括:S3.1:在实际变电站系统中输入电压激励U(k)并采样测量得到N+1个实际输出响应y(t)=[P(t),Q(t)]<sup>T</sup>,其中t表示测量样本序号,取值范围为0≤t≤N,P(t)为时刻t的负荷点在暂态过程中所消耗的有功功率值,Q(t)为时刻t的负荷点在暂态过程中所消耗的无功功率值;S3.2:将实际输出响应对应时刻的离散时间序列{u(0),u(1),…,u(N)}输入负荷模型,负荷模型的初始稳态条件为:初始电压U<sub>0</sub>=u(0),初始有功功率P<sub>0</sub>=P(0),初始无功功率Q<sub>0</sub>=Q(0);设定独立待辨识参数向量α的n个初始值,采用粒子群算法对以下目标函数寻优:<maths num="0001"><![CDATA[<math><mrow><mi>min</mi><mi>J</mi><mrow><mo>(</mo><mi>&alpha;</mi><mo>,</mo><mi>&beta;</mi><mo>)</mo></mrow><mo>=</mo><mi>min</mi><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>0</mn></mrow><mi>N</mi></munderover><msup><mrow><mo>[</mo><mi>y</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>y</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>]</mo></mrow><mi>T</mi></msup><mo>[</mo><mi>y</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>y</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>]</mo></mrow></math>]]></maths>其中,y<sub>m</sub>(t)=[P<sub>m</sub>(t),Q<sub>m</sub>(t)]<sup>T</sup>为输入u(t)时负荷模型得到的输出响应;粒子群算法中的速度调整因子a(m)按照以下公式确定:<maths num="0002"><![CDATA[<math><mrow><mi>a</mi><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>F</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>D</mi></munderover><msqrt><msub><mi>a</mi><mi>max</mi></msub><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>a</mi><mi>min</mi></msub><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow></msqrt></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msub><mi>a</mi><mi>max</mi></msub><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>D</mi></munderover><msup><mrow><mo>(</mo><msubsup><mi>n</mi><mi>ij</mi><mn>1</mn></msubsup><mo>-</mo><msubsup><mi>c</mi><mi>ij</mi><mn>1</mn></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>D</mi></munderover><msup><mrow><mo>(</mo><msubsup><mi>c</mi><mi>ij</mi><mn>1</mn></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msub><mi>a</mi><mi>min</mi></msub><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>D</mi></munderover><msup><mrow><mo>(</mo><msubsup><mi>n</mi><mi>ij</mi><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>c</mi><mi>ij</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>D</mi></munderover><msup><mrow><mo>(</mo><msubsup><mi>c</mi><mi>ij</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac></mrow></math>]]></maths>其中,F为粒子群中粒子个数,D为独立待辨识参数向量α的维度,<img file="FDA0000473353530000024.GIF" wi="170" he="82" />分别表示设置的各个参数的最大速度及最小速度对应的拟合数据,<img file="FDA0000473353530000025.GIF" wi="157" he="82" />分别表示设置的各个参数的最大速度及最小速度对应的测量数据计算值。S4:各等值感应电动机的负荷模型及其负荷参数辨识结果作为变电站挂接母线下等值负荷模型及负荷参数,将不同站点的负荷模型进行统计,得到地区电网负荷模型参数库。
地址 100031 北京市西城区西长安街86号