留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

基于扰动后暂态量的戴维南等值参数辨识方法

上一篇

下一篇

肖俊, 司徒友, 袁炜灯, 等. 基于扰动后暂态量的戴维南等值参数辨识方法[J]. 西南师范大学学报(自然科学版), 2020, 45(6): 113-119. doi: 10.13718/j.cnki.xsxb.2020.06.017
引用本文: 肖俊, 司徒友, 袁炜灯, 等. 基于扰动后暂态量的戴维南等值参数辨识方法[J]. 西南师范大学学报(自然科学版), 2020, 45(6): 113-119. doi: 10.13718/j.cnki.xsxb.2020.06.017
Jun XIAO, You SITU, Wei-deng YUAN, et al. Identification Method for Thevenin Equivalent Parameters under Transient Disturbance Data[J]. Journal of Southwest China Normal University(Natural Science Edition), 2020, 45(6): 113-119. doi: 10.13718/j.cnki.xsxb.2020.06.017
Citation: Jun XIAO, You SITU, Wei-deng YUAN, et al. Identification Method for Thevenin Equivalent Parameters under Transient Disturbance Data[J]. Journal of Southwest China Normal University(Natural Science Edition), 2020, 45(6): 113-119. doi: 10.13718/j.cnki.xsxb.2020.06.017

基于扰动后暂态量的戴维南等值参数辨识方法

  • 基金项目: 广东电网有限责任公司科技项目GDKJXM20162533(031900KK52160081)
详细信息
    作者简介:

    肖俊(1983-), 男, 工程师, 博士, 主要从事电力系统计算分析、稳定与控制、交直流系统、调度运行、检修策划等方面的研究 .

  • 中图分类号: TM74

Identification Method for Thevenin Equivalent Parameters under Transient Disturbance Data

  • 摘要: 针对现有戴维南等值参数辨识算法对等值系统发生扰动时适应能力不强的情况,提出了基于系统扰动后暂态量的戴维南等值参数辨识方法.首先对等值系统扰动下传统等值方法无法辨识戴维南等值参数进行了分析,然后阐明了基于扰动后暂态电压电流信号辨识等值参数的实现原理,利用最小二乘法求出端口电流电压约束关系下微分方程的系数,对比等值网络稳态时的等值参数,从而获得该等值网络暂态状态的戴维南等值参数.通过对简单电路网络和IEEE9节点系统进行仿真分析,结果验证了此方法在等值系统扰动后对戴维南等值参数进行跟踪辨识的可行性和有效性.
  • 加载中
  • 图 1  戴维南等值示意图

    图 2  电力系统等值

    图 3  简单电路网络模型

    图 4  3机9节点系统

    图 5  稳态实测曲线与等值系统仿真曲线

    表 1  简单电路网络模型参数

    参数 R1 R2 R3 R4 R5
    参数值 40 Ω 40 Ω 10 Ω 20 Ω 20 Ω
    Rn1 Rn2 L1 L2 L3 C
    5 Ω 10 Ω 0.5 H 0.5 H 0.2 H 1e-5F
    下载: 导出CSV

    表 2  网络N2扰动下辨识结果(S1开, S2由开到闭)

    对应的网络侧 R L/H Eeq/V
    网络N1 24.998 0.500 220.728
    网络N2 -24.998 -0.500
    下载: 导出CSV

    表 3  网络N1扰动下辨识结果(S2开, S1由开到闭)

    对应的网络侧 R L/H Eeq/V
    网络N1 -41.730 -0.238
    网络N2 41.730 0.238 35.074
    下载: 导出CSV

    表 4  负荷侧扰动下的系统等值参数辨识结果

    等值系统 R X/H Eeq/105V
    系统侧 A相 12.026 1 0.167 5 1.423 3
    B相 11.112 8 0.167 3 1.424 5
    C相 11.526 9 0.165 4 1.421 3
    负荷侧 A相 -12.026 1 -0.167 5
    B相 -11.112 8 -0.167 3
    C相 -11.526 9 -0.165 4
    下载: 导出CSV

    表 5  系统侧扰动下的负荷等值参数辨识结果

    等值系统 R X/H Eeq/105V
    系统侧 A相 -469.017 4 -0.512 4
    B相 -468.413 0 -0.516 2
    C相 -469.215 2 -0.517 1
    负荷侧 A相 469.017 4 0.512 4 0.003 9
    B相 468.413 0 0.516 2 0.004 0
    C相 469.215 2 0.517 1 0.003 9
    下载: 导出CSV
  • [1] 汤涌, 林伟芳, 孙华东, 等.基于戴维南等值跟踪的电压失稳和功角失稳的判别方法[J].中国电机工程学报, 2009, 29(25):1-6. doi: 10.3321/j.issn:0258-8013.2009.25.001
    [2] 李来福, 柳焯.基于戴维南等值参数的紧急态势分析[J].电网技术, 2008, 32(21):64-67.
    [3] 刘道伟, 宋墩文, 王红印, 等.基于WAMS和EMS的电压稳定在线评估系统[J].电网技术, 2014, 38(7):1934-1938. doi: http://d.old.wanfangdata.com.cn/Periodical/dwjs201407033
    [4] 汤涌, 顾卓远.区域电网远端故障引发互联系统稳定破坏现象探究[J].电网技术, 2014, 38(1):154-159. doi: http://d.old.wanfangdata.com.cn/Periodical/dwjs201401024
    [5] 李连伟, 吴政球, 钟浩, 等.基于节点戴维南等值的静态电压稳定裕度快速求解[J].中国电机工程学报, 2010, 30(4):79-83. doi: http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgdjgcxb201004012
    [6] 林越, 孟祥玮, 刘廷章, 等.计入交易偏好的多微网电力系统纳什均衡分析[J].西南师范大学学报(自然科学版), 2017, 42(10):100-105. doi: http://xbgjxt.swu.edu.cn/article/doi/10.13718/j.cnki.xsxb.2017.10.018
    [7] 颜伟, 文一宇, 余娟, 等.基于戴维南等值的静态电压稳定广域切负荷控制策略[J].电网技术, 2011, 35(8):88-92. doi: http://d.old.wanfangdata.com.cn/Periodical/dwjs201108016
    [8] 马世英, 仲悟之, 汤涌, 等.基于多信息源的电压安全全过程预警及防控系统[J].电网技术, 2013, 37(9):2459-2465. doi: http://d.old.wanfangdata.com.cn/Periodical/dwjs201309013
    [9] 李来福, 柳进, 于继来, 等.节点戴维南等值参数在线跟踪简捷算法[J].中国电机工程学报, 2006, 26(10):40-44. doi: 10.3321/j.issn:0258-8013.2006.10.008
    [10] VU K, BEGOVIC M, NOVOSEL D, et al. Use of Local Measurements to Estimate Voltage-Stability Margin[J]. IEEE Transaction on Power Systems, 1999, 14(3):1029-1035. doi: 10.1109/59.780916
    [11] 王漪, 柳焯.基于戴维南等值的系统参数跟踪估计[J].电网技术, 2000, 24(11):28-30. doi: 10.3321/j.issn:1000-3673.2000.11.006
    [12] 廖国栋, 王晓茹.电力系统戴维南等值参数辨识的不确定性模型[J].中国电机工程学报, 2008, 28(28):74-79. doi: 10.3321/j.issn:0258-8013.2008.28.013
    [13] 李来福, 于继来, 柳焯.戴维南等值参数的漂移问题研究[J].中国电机工程学报, 2005, 25(20):1-5. doi: 10.3321/j.issn:0258-8013.2005.20.001
    [14] 汤涌, 孙华东, 易俊, 等.基于全微分的戴维南等值参数跟踪算法[J].中国电机工程学报, 2009, 29(13):48-53. doi: 10.3321/j.issn:0258-8013.2009.13.008
    [15] 孙华东, 陈树勇, 常晓鹏, 等.适应大扰动的电力系统戴维南等值参数跟踪算法[J].中国电机工程学报, 2012, 32(32):126-132. doi: http://d.old.wanfangdata.com.cn/Periodical/zgdjgcxb201222019
    [16] 李卫星, 牟晓明, 李志民.电力系统戴维南等值参数的解析与思考[J].中国电机工程学报, 2012, 32(增刊):28-34. doi: http://d.old.wanfangdata.com.cn/Periodical/zgdjgcxb2012z1005
    [17] 罗华伟, 吴政球, 戴庆华, 等.电网戴维南等值参数的快速计算[J].中国电机工程学报, 2009, 29(1):35-39. doi: 10.3321/j.issn:0258-8013.2009.01.006
    [18] 李卓艺, 刘俊勇, 刘友波, 等.系统持续扰动下的电力系统戴维南等值参数在线辨识[J].中国电机工程学报, 2015, 35(8):1900-1908.
    [19] 罗建, 徐馨, 杨桦, 等.微小扰动状况下的戴维南等值参数辨识方法[J].中国电机工程学报, 2014, 34(S):61-66. doi: http://d.old.wanfangdata.com.cn/Periodical/zgdjgcxb2014z1009
  • 加载中
图( 5) 表( 5)
计量
  • 文章访问数:  1218
  • HTML全文浏览数:  1218
  • PDF下载数:  125
  • 施引文献:  0
出版历程
  • 收稿日期:  2019-09-22
  • 刊出日期:  2020-06-20

基于扰动后暂态量的戴维南等值参数辨识方法

    作者简介: 肖俊(1983-), 男, 工程师, 博士, 主要从事电力系统计算分析、稳定与控制、交直流系统、调度运行、检修策划等方面的研究
  • 广东电网有限责任公司东莞供电局, 广东 东莞 523008
基金项目:  广东电网有限责任公司科技项目GDKJXM20162533(031900KK52160081)

摘要: 针对现有戴维南等值参数辨识算法对等值系统发生扰动时适应能力不强的情况,提出了基于系统扰动后暂态量的戴维南等值参数辨识方法.首先对等值系统扰动下传统等值方法无法辨识戴维南等值参数进行了分析,然后阐明了基于扰动后暂态电压电流信号辨识等值参数的实现原理,利用最小二乘法求出端口电流电压约束关系下微分方程的系数,对比等值网络稳态时的等值参数,从而获得该等值网络暂态状态的戴维南等值参数.通过对简单电路网络和IEEE9节点系统进行仿真分析,结果验证了此方法在等值系统扰动后对戴维南等值参数进行跟踪辨识的可行性和有效性.

English Abstract

  • 戴维南等值方法通过简化电力系统, 提高了安全稳定评估的效率[1-8], 而戴维南等值参数的实时辨识是基于戴维南等值研究及应用的关键.现有的戴维南等值参数实时辨识方法通常就等值网络的负荷侧扰动状况进行研究[9-13], 而对于等值网络的系统侧扰动状况研究不多.

    对于等值网络的扰动状况, 由于前后两个运行状况下的等值网络参数的变化, 那些基于两个运行状况或多个运行状况下的量测量进行戴维南等值参数实时辨识的方法存在辨识困难, 显然利用一个运行状况下的量测量进行戴维南等值参数的实时辨识可以满足等值网络系统扰动状况下的参数辨识要求.

    现有的利用一个运行状况下的量测量进行戴维南等值参数实时辨识的方法主要是基于全微分的戴维南等值参数跟踪算法[14-18], 在等值网络等值参数初值确定的条件下, 该算法在原理上可以跟踪等值网络参数的变化, 但该算法对系统小扰动的适应能力不强, 等值网络参数的渐进变化将影响算法的辨识效果.

    为改进对系统小扰动的适应能力, 文献[19]提出了一种小扰动状况下的戴维南等值参数辨识的方法, 该方法在系统和负荷扰动状况下都能跟踪辨识等值网络的戴维南等值参数.虽然文献[19]给出了戴维南等值参数辨识方法的基本原理, 但未就等值网络系统扰动状况下的戴维南等值参数辨识实现原理加以阐明, 为此本文就等值网络系统扰动状况下的戴维南等值参数辨识实现原理进行阐明, 最后通过仿真验证了辨识方法的准确性和有效性.

  • 若从某一负荷向系统看进去, 任意k时刻的实际电力系统可以等值为一个含有电压源和阻抗部分, 与外部负荷构成图 1所示的简单两节点系统, 其数学模型可以描述为

    其中, $\dot E$SZS分别为某等值系统在k时的戴维南等值电势和等值阻抗, $\dot U$k时刻的负荷母线电压, $\dot I$为负荷电流, PQ分别为负荷侧的有功功率、无功功率.等值母线处电压$\dot U$和电流$\dot I$与负荷ZL的关系如下

  • 在等值母线处对左侧(系统侧)进行等值, 假定相邻时刻戴维南等值参数未发生太大变化, 则系统戴维南等值阻抗可由下式计算:

    式中, 上标表示相邻时刻等值母线处电压电流的测量信息.

    若扰动发生在系统侧, 负荷侧参数保持不变, 假定系统侧等值阻抗ZS变为ZS+ΔZS, 等值电势$\dot E$S保持不变, 将式(2)带入式(3)则有

    由上式可知, 若扰动仅存在于系统侧, 则辨识得到的等值参数为负荷侧的等值阻抗.因此, 等值系统的系统侧发生扰动, 传统戴维南等值方法无法准确跟踪等值参数的变化.

    文献[19]的小扰动状况下的戴维南等值参数辨识方法利用暂态电量信号(当前运行状况下的量测量与前一运行状况下的量测量的差值)对等值网络参数进行辨识, 对于等值网络系统侧扰动而言, 由于扰动发生前后网络的参数不一致, 导致未知参数的个数增加, 无法对等值参数进行有效的辨识.等值网络系统侧扰动后的参数不变, 利用扰动后暂态过程中的量测量可有效辨识当前运行状况下的等值网络参数.

  • 利用扰动后的量测信息辨识网络参数, 从而得到当前运行状况下结构和参数唯一确定的等值网络, 电力系统等值过程如图 2所示.

    根据电路理论, 等值网络的电压电流约束关系可用K阶微分方程表示:

    式中, v(t)和i(t)分别表示等值网络端口的电压和电流信号, v(t)(k)i(t)(k)分别表示电压和电流信号的k阶微分求导, vs(t)表示等值网络的等效电源, vs(t)(k)表示等效电源的k阶微分求导, 微分方程的阶次K和系数akbkck(k=0, …, K)决定了等值网络的特性.

    在系统实际运行中, 等值网络系统侧发生扰动的情况较为普遍, 等值网络系统侧扰动变化, 其参数会随之改变, 导致高阶微分方程中各系数改变, 其表征的网络特性也随之发生改变.等值网络系统侧发生扰动后的结构和参数保持不变, 对应当前运行状况下的等值参数, 高阶微分方程表征的等值网络的特性不变, 因而可利用扰动后的暂态过程中的电量信息辨识等值网络参数.

    扰动后网络端口的电压电流信号为周期分量和暂态分量的叠加, 即

    式中, iT(t)、vT(t)分别表示端口暂态电压、电流信号的周期分量, iNT(t)、vNT(t)表示相应暂态信号的暂态分量, vsT(t)表示暂态过程中的等效电源.

    通过提取网络端口电压电流的暂态分量信号, 利用电压电流的暂态分量信号辨识等值网络当前运行状况下的等值参数.将式(6)带入式(5)中得到

    取暂态过程中电压电流的相邻半个周期的信号, 可得:

    由式(7)和式(8)可得:

    式中, $\mathit{\Delta }{v_{NT}}(t) = {v_{NT}}\left( {t + \frac{T}{2}} \right) + {v_{NT}}(t), \mathit{\Delta }{i_{NT}}(t) = {i_{NT}}\left( {t + \frac{T}{2}} \right) + {i_{NT}}(t)$分别表示电压和电流暂态分量信号.

    式(9)表示的高阶微分方程能够有效表征等值网络的特性.根据网络的实际特性建立相应的等值网络模型, 利用端口提取的暂态过程中的电压电流以及扰动后的稳态信号求解高阶微分方程系数akbk, 通过建立高阶微分方程系数与等值模型参数之间的关系, 即可求出等值模型的参数.

    利用暂态分量信号辨识等值参数, 具体方法如下:令ΔvNT=Vmsin(ωt+φv), ΔiNT=Imsin(ωt+φi), 带入微分方程(9)中.

    K为偶数时, 令N=K/2, 并计算

    K为奇数时, 令N=(K-1)/2, 并计算

    由此可得:

    上式中左边部分提取$\sqrt {A_0^2 + A_1^2} $, 并令$\sin \theta = \frac{{{A_1}}}{{\sqrt {A_0^2 + A_1^2} }}, \cos \theta = \frac{{{A_0}}}{{\sqrt {A_0^2 + A_1^2} }}$, 推导可得:

    对上式坐标进行平移变换, 可得

    将上式进一步整理可得:

    将上式改写为$\mathit{\Delta }v = {R_{eq}}\mathit{\Delta }i + {L_{eq}}\frac{{{\rm{d}}\mathit{\Delta }i}}{{{\rm{d}}t}}$的形式, 得

    即可得到等值网络的等值电阻和等值电抗为

    利用端口提取的暂态分量信号进行戴维南参数辨识, 能实时跟踪等值网络的参数变化, 得到的结果能进一步用于电力系统计算和分析, 使该等值系统更贴近于并反映系统的实际运行状况.

  • 为了说明本文方法正确可行, 通过建立简单电路网络模型进行验证, 电路网络模型如图 3所示.

    图 3所示的简单电路网络中, Us1=220 V, Us2=35 V, 各支路阻抗参数如表 1所示, 将电路网络分为等值网络N1和外部网络N2.利用理想断路器S1和S2动作分别模拟等值网络N1的扰动和外部网络N2的扰动, 由此产生小扰动信号.仿真时间为0.2 s, 在t=0.1 s时理想断路器断开线路, 记录检测点电压电流信号, 采样频率设置为10 kHz. S1断开、闭合状态下网络N1的等值电阻电感真实值分别为:25 Ω, 0.5 H及20 Ω, 0.5 H. S2断开、闭合状态下网络N2的等值电阻电感真实值分别为:41.724 4 Ω、0.238 4 H及36.724 4 Ω, 0.238 4 H.

    保持等值网络N1的断路器S1处于断开状态, 断开或闭合外部网络N2的S2, 模拟等值系统外部网络的扰动, 利用扰动状况下的暂态信号辨识戴维南等值参数, 结果如表 2所示.类似地, 保持负荷侧网络的断路器S2处于断开状态, 断开或闭合等值网络系统侧的S1, 模拟等值系统系统侧网络的扰动, 辨识的戴维南等值参数如表 3所示.

    表 2表 3可知, 基于扰动后暂态分量能准确辨识戴维南等值参数.

  • 图 4所示, 在IEEE9节点系统进行戴维南等值, 负荷设置如下:Bus5所接负荷为恒阻抗负荷和感应电动机的组合, Bus7和Bus8处的负荷大小不变、功率因数恒定, 通过调整Bus5负荷变化、系统结构变化进行辨识来验证本方法的准确性.

    针对图 4所示的系统, 分别施以如下的扰动进行戴维南等值参数辨识:

    1) 外部负荷变化:对Bus5处负荷减小1%作为外部扰动;

    2) 等值系统网络结构变化:改变Bus4与Bus5间线路L45的长度作为等值系统扰动.

    以上情况的仿真时间均为1.8 s, 扰动开始时刻均为1.5 s, 在Bus5处测量电压电流信号, 采样频率为10 kHz, 利用扰动下的三相电压电流暂态信号分别辨识戴维南等值参数, 辨识结果如表 4表 5所示.

    表 4表 5的戴维南等值参数辨识结果可知, 扰动后的A、B、C三相暂态电流电压信号辨识的等值参数非常接近, 说明利用扰动后暂态分量能够有效辨识戴维南等值参数.针对等值系统网络结构变化引起的等值参数变化, 该方法能够利用一个运行状态下的量测数据实时辨识等值系统的戴维南参数.

    为了验证戴维南等值参数辨识结果的准确性, 对比稳态下Bus5处原系统实测曲线与等值系统仿真曲线.限于篇幅, 仅给出A相的电压电流, 如图 5所示, 其中, (a)和(b)为等值系统负荷变化时的电压电流稳态响应曲线, (c)和(d)为等值网络系统侧网络结构变化时的电压电流稳态响应曲线.

    从以上仿真结果可以看出, 在等值系统负荷或网络结构变化时, 稳态响应曲线与实测曲线差别很小, 参数辨识的结果接近真值, 说明本文采用的戴维南等值参数辨识方法在等值系统内外扰动下均能准确辨识系统参数, 有效地解决了等值系统结构和参数变化的辨识问题.

  • 本文提出了一种基于扰动后暂态分量信号的戴维南等值参数辨识方法, 通过建立能够表征等值网络实际特性的电路网络模型, 利用扰动后暂态电压电流信号求解出与之匹配的高阶微分方程系数, 然后等值网络稳态时参数与微分方程系数相关联, 从而求出等值系统的实时参数.简单电路网络和IEEE9节点系统的仿真结果表明, 利用端口提取的暂态分量信号能够跟踪等值网络的实时变化, 有效解决了等值网络扰动下等值参数的准确跟踪问题.

参考文献 (19)

目录

/

返回文章
返回