发明名称 一种基于暂态电压安全裕度的超实时重合时序整定方法
摘要 本发明涉及一种基于暂态电压安全裕度的超实时重合时序整定方法,属电力系统稳定与控制技术领域。当交流输电线路故障时,根据将广域信息采集系统提供的故障引起的断路器跳闸后的节点电压和电流等实测数据,采用最小二乘法在线辨识导纳阵,修改辨识出的导纳阵分别模拟线路首端重合和末端重合,采用改进欧拉法预测两种重合时序下的节点电压变化轨迹,并计算暂态电压安全裕度,实现对电压变化轨迹的量化评估,选取暂态电压安全裕度较大者对应重合时序为优化重合方案。本方法利用最小二乘法和改进欧拉法预测节点电压变化轨迹,优化重合方案,在提高预测精度的同时,解决了传统“离线计算、在线匹配”方法得到的重合时序方案运行方式和故障位置失配问题。
申请公布号 CN102403720A 申请公布日期 2012.04.04
申请号 CN201110375495.2 申请日期 2011.11.23
申请人 昆明理工大学 发明人 孙士云;束洪春;孙向飞;董俊
分类号 H02J3/00(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 代理人
主权项 1.一种基于暂态电压安全裕度的超实时重合时序整定方法,其特征在于具体步骤如下:1.1通过WAMS(wide area measurement system)广域信息采集系统,采集扰动后节点电压和电流实测数据,然后根据式(1)网络节点方程和式(2)导纳阵最小二乘辨识值方程,用最小二乘法辨识法求取导纳阵<img file="173760DEST_PATH_IMAGE001.GIF" wi="17" he="18" />;<img file="887638DEST_PATH_IMAGE002.GIF" wi="49" he="18" />(1)<img file="667375DEST_PATH_IMAGE003.GIF" wi="105" he="33" />(2)式中,<img file="875633DEST_PATH_IMAGE004.GIF" wi="17" he="18" />、<img file="690006DEST_PATH_IMAGE005.GIF" wi="13" he="17" />为节点电压、电流向量,<img file="574785DEST_PATH_IMAGE006.GIF" wi="16" he="29" />为导纳阵<img file="841818DEST_PATH_IMAGE001.GIF" wi="17" he="18" />的最小二乘辨识值;1.2采用下述方法,模拟线路首端投入重合闸,修改导纳阵<b><i>Y</i></b>中故障线路对应元素,使故障线路首端断线状况消失,得到导纳阵<b><i>Y</i></b><sub>1</sub>:若故障线路首端断路器对应母线侧和线路侧的节点编号分别为<i>p、q,</i>按式(3)修改导纳阵中<i>p、q</i>点对应自导纳和互导纳:<img file="40719DEST_PATH_IMAGE007.GIF" wi="109" he="84" />(3)式中,<img file="519717DEST_PATH_IMAGE008.GIF" wi="25" he="26" />、<img file="513081DEST_PATH_IMAGE009.GIF" wi="25" he="26" />为节点<i>p、q</i>的自导纳,<img file="329727DEST_PATH_IMAGE010.GIF" wi="25" he="26" />为节点<i>p、q</i>间的互导纳,<img file="332318DEST_PATH_IMAGE011.GIF" wi="20" he="25" />为模拟<i>p、q</i>点间的故障消失而采用的导纳,其标幺值为9999 pu;若线路发生单相故障并采用单相重合闸时,按上述方法修改正、负、零三序导纳阵中对应元素;1.3将通过WAMS广域信息采集系统采集的第<img file="855703DEST_PATH_IMAGE012.GIF" wi="14" he="16" />个时步的实测发电机功角<img file="833018DEST_PATH_IMAGE013.GIF" wi="20" he="25" />、角速度<img file="136960DEST_PATH_IMAGE014.GIF" wi="21" he="25" />和电磁功率<img file="677663DEST_PATH_IMAGE015.GIF" wi="24" he="25" />代入式(4)转子运动方程,利用式(5)改进欧拉法求得功角和角速度第<img file="55555DEST_PATH_IMAGE016.GIF" wi="33" he="20" />个时步的预测值<img file="469350DEST_PATH_IMAGE017.GIF" wi="28" he="26" />和<img file="932692DEST_PATH_IMAGE018.GIF" wi="30" he="26" />;<img file="339402DEST_PATH_IMAGE019.GIF" wi="168" he="86" />(4)式中,<img file="571801DEST_PATH_IMAGE020.GIF" wi="16" he="20" />、<img file="343448DEST_PATH_IMAGE021.GIF" wi="17" he="16" />为发电机功角、角速度,<img file="107135DEST_PATH_IMAGE022.GIF" wi="21" he="25" />、<img file="989641DEST_PATH_IMAGE023.GIF" wi="18" he="25" />为发电机机械、电磁功率,<img file="404442DEST_PATH_IMAGE024.GIF" wi="9" he="17" />为时间,<img file="346990DEST_PATH_IMAGE025.GIF" wi="22" he="18" />为惯性时间常数;<img file="657361DEST_PATH_IMAGE026.GIF" wi="202" he="140" />(5)式中,<img file="343557DEST_PATH_IMAGE013.GIF" wi="20" he="25" />、<img file="550548DEST_PATH_IMAGE014.GIF" wi="21" he="25" />为第<i>n</i>个时步的实测功角、角速度,<img file="477046DEST_PATH_IMAGE017.GIF" wi="28" he="26" />、<img file="402277DEST_PATH_IMAGE018.GIF" wi="30" he="26" />为第<img file="688902DEST_PATH_IMAGE016.GIF" wi="33" he="20" />个时步的功角、角速度预测值,<i>h</i>为采样时间间隔;1.4按式(6)计算功角和角速度增量:<img file="15978DEST_PATH_IMAGE027.GIF" wi="109" he="26" />,<img file="847799DEST_PATH_IMAGE028.GIF" wi="114" he="26" />(6)1.5基于导纳阵<b><i>Y</i></b><sub>1</sub>得到雅可比矩阵元素,将功角和角速度增量代入式(7-1)和(7-2),分别求得负荷节点的电压幅值、相角和电磁功率增量<img file="260326DEST_PATH_IMAGE029.GIF" wi="49" he="26" />、<img file="350641DEST_PATH_IMAGE030.GIF" wi="45" he="26" />、<img file="532224DEST_PATH_IMAGE031.GIF" wi="42" he="26" />;<img file="987476DEST_PATH_IMAGE032.GIF" wi="330" he="52" />(7-1)<img file="434769DEST_PATH_IMAGE033.GIF" wi="390" he="52" />(7-2)式中,<img file="880DEST_PATH_IMAGE034.GIF" wi="23" he="19" />、<img file="364865DEST_PATH_IMAGE035.GIF" wi="27" he="19" />为发电机和负荷节点的有功、无功增量,△<i>δ</i>为发电机的功角增量,<img file="991018DEST_PATH_IMAGE036.GIF" wi="26" he="19" />、<img file="922678DEST_PATH_IMAGE037.GIF" wi="29" he="18" />为负荷节点相角、电压幅值增量,<i>H</i><sub><i>11</i></sub><i>、N</i><sub><i>11</i></sub><i>、H</i><sub><i>12</i></sub><i>、N</i><sub><i>12</i></sub><i>、J</i><sub><i>11</i></sub><i>、L</i><sub><i>11</i></sub><i>、J</i><sub><i>12</i></sub><i>、L</i><sub><i>12</i></sub><i>、H</i><sub><i>21</i></sub><i>、N</i><sub><i>21</i></sub><i>、H</i><sub><i>22</i></sub><i>、N</i><sub><i>22</i></sub><i>、J</i><sub><i>21</i></sub><i>、L</i><sub><i>21</i></sub><i>、J</i><sub><i>22</i></sub><i>、L</i><sub><i>22</i></sub>为雅可比矩阵元素,<img file="292479DEST_PATH_IMAGE038.GIF" wi="98" he="25" /><img file="448654DEST_PATH_IMAGE039.GIF" wi="45" he="25" />为静态负荷模型下电压特性系数,<img file="308026DEST_PATH_IMAGE040.GIF" wi="22" he="25" />、<img file="916862DEST_PATH_IMAGE041.GIF" wi="18" he="25" />、<img file="637824DEST_PATH_IMAGE042.GIF" wi="21" he="25" />分别为负荷节点电压、有功功率和无功功率;<img file="648505DEST_PATH_IMAGE040.GIF" wi="22" he="25" />、<img file="678778DEST_PATH_IMAGE041.GIF" wi="18" he="25" />、<img file="774910DEST_PATH_IMAGE042.GIF" wi="21" he="25" />与电压特性系数之间的关系式为:<img file="299564DEST_PATH_IMAGE043.GIF" wi="148" he="54" />;1.6按式(8)计算第<img file="164751DEST_PATH_IMAGE016.GIF" wi="33" he="20" />个时步负荷节点的电压幅值<img file="303609DEST_PATH_IMAGE044.GIF" wi="38" he="26" />、相角<img file="949354DEST_PATH_IMAGE045.GIF" wi="34" he="26" />和电磁功率<img file="464649DEST_PATH_IMAGE046.GIF" wi="33" he="26" />的预测值:<img file="997392DEST_PATH_IMAGE047.GIF" wi="137" he="26" />,<img file="307151DEST_PATH_IMAGE048.GIF" wi="126" he="26" />,<img file="440192DEST_PATH_IMAGE049.GIF" wi="120" he="26" />(8)1.7将第<img file="493599DEST_PATH_IMAGE016.GIF" wi="33" he="20" />个时步电磁功率的预测值<img file="143498DEST_PATH_IMAGE046.GIF" wi="33" he="26" />代入式(4),求得<img file="624158DEST_PATH_IMAGE050.GIF" wi="52" he="26" />;用第<img file="182179DEST_PATH_IMAGE016.GIF" wi="33" he="20" />个时步负荷节点的电压幅值和相角的预测值<img file="836014DEST_PATH_IMAGE044.GIF" wi="38" he="26" />、<img file="530300DEST_PATH_IMAGE045.GIF" wi="34" he="26" />替换第<img file="260490DEST_PATH_IMAGE012.GIF" wi="14" he="16" />个时步电压幅值和相角值<img file="40227DEST_PATH_IMAGE051.GIF" wi="29" he="26" />、<img file="435436DEST_PATH_IMAGE052.GIF" wi="25" he="26" />,修正雅可比矩阵矩阵中的元素<i>H</i><sub><i>11</i></sub><i>、N</i><sub><i>11</i></sub><i>、H</i><sub><i>12</i></sub><i>、N</i><sub><i>12</i></sub><i>、J</i><sub><i>11</i></sub><i>、L</i><sub><i>11</i></sub><i>、J</i><sub><i>12</i></sub><i>、L</i><sub><i>12</i></sub><i>、H</i><sub><i>21</i></sub><i>、N</i><sub><i>21</i></sub><i>、H</i><sub><i>22</i></sub><i>、N</i><sub><i>22</i></sub><i>、J</i><sub><i>21</i></sub><i>、L</i><sub><i>21</i></sub><i>、J</i><sub><i>22</i></sub><i>、L</i><sub><i>22</i></sub>;1.8将<img file="312126DEST_PATH_IMAGE050.GIF" wi="52" he="26" />代入式(5),求得第<img file="134588DEST_PATH_IMAGE016.GIF" wi="33" he="20" />个时步的功角和角速度的校正值<img file="214671DEST_PATH_IMAGE053.GIF" wi="28" he="26" />和<img file="413571DEST_PATH_IMAGE054.GIF" wi="30" he="26" />;1.9重复步骤1.3-1.8,直至预测和校正过程结束,得到负荷节点的电压预测轨迹和发电机节点功角预测轨迹;1.10利用式(9)和(10),对负荷节点电压变化轨迹进行量化评估,得到线路首端重合时的暂态电压安全裕度<img file="144767DEST_PATH_IMAGE055.GIF" wi="24" he="25" />:<img file="138130DEST_PATH_IMAGE056.GIF" wi="18" he="25" />=min(<img file="892460DEST_PATH_IMAGE057.GIF" wi="21" he="25" />,<img file="442521DEST_PATH_IMAGE058.GIF" wi="26" he="26" />)                            (9)<img file="965906DEST_PATH_IMAGE059.GIF" wi="185" he="73" />(10)式中,<img file="458067DEST_PATH_IMAGE060.GIF" wi="20" he="25" />为暂态电压安全裕度,<img file="699693DEST_PATH_IMAGE061.GIF" wi="22" he="25" />为暂态电压稳定裕度,<img file="784936DEST_PATH_IMAGE062.GIF" wi="28" he="26" />为暂态电压偏移可接受性裕度,<i>H</i>为感应电动机的惯性时间常数(秒),<i>s</i>为感应电动机的转差(p.u.),<img file="428407DEST_PATH_IMAGE063.GIF" wi="26" he="25" />为感应电动机的机械功率(W或p.u.),<img file="763573DEST_PATH_IMAGE064.GIF" wi="20" he="25" />为母线的电压偏移门槛值(V或p.u.),<img file="554812DEST_PATH_IMAGE065.GIF" wi="20" he="25" />为超越偏移门槛值的允许持续时间(秒),<i>k</i>为把临界电压偏移时间换算成电压的折算因子;1.11采用下述方法,模拟线路末端投入重合闸,修改导纳阵<img file="633626DEST_PATH_IMAGE001.GIF" wi="17" he="18" />,得到导纳阵<b><i>Y</i></b><sub>2</sub>:若故障线路末端断路器对应母线侧和线路侧的节点编号分别为<i>p’、q’,</i>则导纳阵中<i>p’、q’</i>点对应自导纳和互导纳按式(11)修改:<img file="944653DEST_PATH_IMAGE066.GIF" wi="118" he="92" />(11)式中,<img file="716300DEST_PATH_IMAGE067.GIF" wi="30" he="29" />、<img file="729255DEST_PATH_IMAGE068.GIF" wi="29" he="29" />为节点<i>p’、q’</i>的自导纳,<img file="611761DEST_PATH_IMAGE069.GIF" wi="30" he="29" />为节点<i>p’、q’</i>间的互导纳,<img file="964245DEST_PATH_IMAGE011.GIF" wi="20" he="25" />为模拟<i>p’、q’</i>点间的故障消失而采用的导纳,其标幺值为9999 pu;若线路发生单相故障并采用单相重合闸时,则按上述方法同时修改正、负、零三序导纳阵中对应元素;重复步骤1.3-1.10,步骤1.5中改为基于导纳阵<b><i>Y</i></b><sub>2</sub>计算得到雅可比矩阵元素,求得线路末端重合时的暂态电压安全裕度<img file="719842DEST_PATH_IMAGE070.GIF" wi="21" he="25" />;1.12比较线路首、末端重合时的暂态电压安全裕度<img file="157777DEST_PATH_IMAGE055.GIF" wi="24" he="25" />,<img file="640711DEST_PATH_IMAGE070.GIF" wi="21" he="25" />,选取较大者对应的重合时序作为输出的优化重合时序方案;若<img file="113280DEST_PATH_IMAGE055.GIF" wi="24" he="25" />&gt;<img file="226730DEST_PATH_IMAGE070.GIF" wi="21" he="25" />,由线路首端首先重合;若<img file="965010DEST_PATH_IMAGE070.GIF" wi="21" he="25" />&gt;<img file="189318DEST_PATH_IMAGE055.GIF" wi="24" he="25" />,则由线路末端首先重合。
地址 650093 云南省昆明市五华区学府路253号