发明名称 一种发电机励磁系统参数可辨识性判定方法
摘要 本发明公开了一种发电机励磁系统参数可辨识性判定方法,涉及发电机励磁系统参数可辨识性的判定。该方法通过构建参数时域灵敏度矩阵实现关联性参数的判断,将待辨识参数划分为良态参数集和病态参数集,取病态参数集中拟频域灵敏度之和最小的作为赋值代表,对赋值代表赋以经验值,并对其所对应的良态参数集进行参数辨识。优点是:有效地提高了发电机励磁系统参数辨识结果的准确性和稳定性,有助于提高电力系统安全稳定的运行水平,具有较高的工程实用价值。
申请公布号 CN104156504A 申请公布日期 2014.11.19
申请号 CN201410347835.4 申请日期 2014.07.21
申请人 国家电网公司;国网湖北省电力公司经济技术研究院;华中科技大学 发明人 彭昌勇;杜治;马蕊;徐敬友;文劲宇;杨东俊
分类号 G06F17/50(2006.01)I;G06F19/00(2011.01)I 主分类号 G06F17/50(2006.01)I
代理机构 武汉荆楚联合知识产权代理有限公司 42215 代理人 王健
主权项 一种发电机励磁系统参数可辨识性判定方法,其特征在于,该方法包含下列步骤:(1)确定励磁系统模型类型;(2)依据励磁系统模型的类型,分别计算待辨识参数的时域灵敏度矩阵和拟频域灵敏度的值;(3)求出时域灵敏度矩阵的秩,该秩的值R即为该励磁系统模型最多可辨识的参数的个数R;(4)在时域灵敏度矩阵中选取R个线性无关的列向量所对应的参数作为良态参数集,其余参数作为该良态参数集对应的病态参数集;列出时域灵敏度矩阵所有的良态参数集及其所对应的病态参数集;(5)计算每个病态参数集中参数的拟频域灵敏度之和,取和最小的作为赋值代表;(6)对赋值代表赋以经验值,并对其所对应的良态参数集进行参数辨识;其中,上述步骤(2)待辨识参数的拟频域灵敏度的值按下列方法计算:(2.1)对步骤(1)所确定励磁系统模型,给定扰动及参数参考值θ<sub>i</sub>,其中i=1,2,…,n,n为励磁系统模型参数的总个数,能仿真得到模型输出变量y的时域动态曲线;(2.2)对该时域输出变量y的自相关函数进行傅氏变换,计算功率谱密度,取其幅值自起始点至衰减到<img file="FDA0000541055100000015.GIF" wi="78" he="60" />倍时的频段区间f<sub>ef</sub>作为有效频率范围;(2.3)取Δθ<sub>i</sub>=0.01×θ<sub>i</sub>,给定参数θ<sub>i</sub>+Δθ<sub>i</sub>,维持其它参数和输入扰动不变,能仿真得到输出变量y的时域动态曲线,对其进行傅氏变换,得到时域曲线y在不同频率情况下的振荡幅值Y<sub>j</sub>((θ<sub>1</sub>,θ<sub>2</sub>,…,θ<sub>i</sub>+Δθ<sub>i</sub>,…,θ<sub>n</sub>),f<sub>j</sub>),其中f<sub>j</sub>代表时域曲线y所包含的第j个频率;(2.4)再给定参数θ<sub>i</sub>‑Δθ<sub>i</sub>,维持其他参数和输入扰动不变,重复步骤(2.3),获得该种情况下的Y<sub>j</sub>((θ<sub>1</sub>,θ<sub>2</sub>,…,θ<sub>i</sub>‑Δθ<sub>i</sub>,…,θ<sub>n</sub>),f<sub>j</sub>);(2.5)按照下式计算待辨识参数的拟频域灵敏度曲线<img file="FDA0000541055100000011.GIF" wi="80" he="70" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><mi>h</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>,</mo><msub><mi>f</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo></mtd></mtr><mtr><mtd><mo>[</mo><mo>|</mo><msub><mi>Y</mi><mi>j</mi></msub><mrow><mo>(</mo><mrow><mo>(</mo><msub><mi>&theta;</mi><mn>1</mn></msub><mo>,</mo><msub><mi>&theta;</mi><mn>2</mn></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>+</mo><mi>&Delta;</mi><msub><mi>&theta;</mi><mi>i</mi></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>&theta;</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>f</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>|</mo><mo>-</mo></mtd></mtr><mtr><mtd><mo>|</mo><msub><mi>Y</mi><mi>j</mi></msub><mrow><mo>(</mo><mrow><mo>(</mo><msub><mi>&theta;</mi><mn>1</mn></msub><mo>,</mo><msub><mi>&theta;</mi><mn>2</mn></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>-</mo><mi>&Delta;</mi><msub><mi>&theta;</mi><mi>i</mi></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>&theta;</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>f</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>|</mo><mo>]</mo><mo>/</mo><mn>2</mn><mi>&Delta;</mi><msub><mi>&theta;</mi><mi>i</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000541055100000012.GIF" wi="734" he="254" /></maths>(2.6)按照下式计算参数拟频域灵敏度曲线在有效频率范围f<sub>ef</sub>内的绝对值的平均值<img file="FDA0000541055100000013.GIF" wi="101" he="68" /><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>H</mi><msub><mi>&theta;</mi><mi>i</mi></msub></msub><mo>=</mo><mfrac><mn>1</mn><msub><mi>K</mi><mi>ef</mi></msub></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>K</mi><mi>ef</mi></msub></munderover><mo>|</mo><mi>h</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>,</mo><msub><mi>f</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>|</mo></mrow>]]></math><img file="FDA0000541055100000014.GIF" wi="417" he="143" /></maths>其中,f<sub>k</sub>为f<sub>j</sub>中第k个小于f<sub>ef</sub>的频率采样点,K<sub>ef</sub>为f<sub>k</sub>的总采样点数。
地址 100031 北京市西城区西长安街86号