发明名称 一种电力系统相量频率测量方法
摘要 一种电力系统相量频率测量方法,对原始的定时间间隔离散化采样值序列先经过有限长冲击响应数字低通滤波器进行低通滤波,以消除原始输入信号中二次及二次以上谐波分量的干扰;对从第i个采样值开始的前推2T<sub>0</sub>时间内的采样值序列分别计算实部<img file="35714dest_path_image001.GIF" wi="76" he="26" />和虚部<img file="704592dest_path_image002.GIF" wi="79" he="31" />,由公式<img file="697956dest_path_image003.GIF" wi="125" he="44" />计算得出相位<img file="389969dest_path_image004.GIF" wi="63" he="26" />;根据公式<img file="392560dest_path_image005.GIF" wi="297" he="41" />计算得出修正后的当前的高精度频率值<img file="94107dest_path_image006.GIF" wi="72" he="28" />。本发明不受电源基准电平波动的干扰;当接入同一块A/D转换器的多路电力系统相量的频率不一致时也能适用,且不会持续受过去计算值的影响;同时精度也满足工程上的实际需求。
申请公布号 CN102236048A 申请公布日期 2011.11.09
申请号 CN201110162398.5 申请日期 2011.06.16
申请人 国网电力科学研究院;南京南瑞集团公司 发明人 任建锋;叶振风;宣筱青;李雪明;曹一中;李祝昆
分类号 G01R23/02(2006.01)I 主分类号 G01R23/02(2006.01)I
代理机构 南京知识律师事务所 32207 代理人 汪旭东
主权项 1.一种电力系统相量频率测量方法,其特征在于,包括以下步骤:1)对原始的定时间间隔离散化采样值序列先经过有限长冲击响应数字低通滤波器进行低通滤波,以消除原始输入信号中二次及二次以上谐波分量的干扰;2)对从第i个采样值开始的前推2T<sub>0</sub>时间内的采样值序列分别计算实部<img file="20770DEST_PATH_IMAGE001.GIF" wi="76" he="62" />和虚部<img file="31451DEST_PATH_IMAGE002.GIF" wi="79" he="62" />,由公式<img file="999407DEST_PATH_IMAGE003.GIF" wi="125" he="62" />计算得出相位<img file="273701DEST_PATH_IMAGE004.GIF" wi="63" he="62" />;3)重复步骤1)、2),可依次得到<img file="985305DEST_PATH_IMAGE005.GIF" wi="99" he="62" />,<img file="53755DEST_PATH_IMAGE006.GIF" wi="99" he="62" />,<img file="192613DEST_PATH_IMAGE007.GIF" wi="99" he="62" />,…;4)假设当前采样值是第p个采样值,则其对应的相位为<img file="697412DEST_PATH_IMAGE008.GIF" wi="69" he="62" />,由公式<img file="212707DEST_PATH_IMAGE009.GIF" wi="293" he="62" />计算得到初次频率值;5)对初次频率值f(2N,p)进行四舍五入处理得到round(f(2N,p));给定一个事先人为设定的频率测量上限值<img file="197981DEST_PATH_IMAGE010.GIF" wi="23" he="62" />和下限值<img file="445422DEST_PATH_IMAGE011.GIF" wi="21" he="62" />(<img file="250567DEST_PATH_IMAGE010.GIF" wi="23" he="62" />和<img file="490925DEST_PATH_IMAGE011.GIF" wi="21" he="62" />均为整数),若round(f(2N,p))<img file="330705DEST_PATH_IMAGE012.GIF" wi="13" he="62" />[<img file="811365DEST_PATH_IMAGE013.GIF" wi="49" he="62" />],则认为计算出来的初次频率值不可信,同时将计算出来的初次频率值清零;否则根据round(f(2N,p))值查找事先人为设定的修正表得到round(f(2N,p))值所对应的v值;修正表在[<img file="307068DEST_PATH_IMAGE013.GIF" wi="49" he="62" />]频率范围内和整数频率一一对应,修正表是这样得到的:给定原始信号序列,<img file="898586DEST_PATH_IMAGE014.GIF" wi="188" he="62" />,此处频率f给定为整数且<img file="858452DEST_PATH_IMAGE015.GIF" wi="89" he="62" />,j=0,1,2,...,10N-1;在原始信号序列中每次取2N个连续的值作为一组采样值序列,由上述步骤1~4计算得到其初次频率值序列f(2N,j);对<img file="696964DEST_PATH_IMAGE016.GIF" wi="88" he="62" />的初次频率值序列f(2N,j),找出使关系式<img file="476701DEST_PATH_IMAGE017.GIF" wi="209" he="62" />(<img file="871910DEST_PATH_IMAGE018.GIF" wi="11" he="62" />是一个人为设定的阀值)成立的所有j值,记为s(j),对所有相邻的两个s(j)值取其差值的绝对值,并对上述所有差值的绝对值取平均值后再四舍五入取整,即可得到该频率f所对应的v值;改变频率f的值,则可一一求出在<img file="623966DEST_PATH_IMAGE019.GIF" wi="60" he="62" />范围内,每个整数频率对应的v值,从而得到该修正表;6)根据公式<img file="446428DEST_PATH_IMAGE020.GIF" wi="297" he="62" />计算得出修正后的当前的高精度频率值<img file="713461DEST_PATH_IMAGE021.GIF" wi="72" he="62" />。
地址 210003 江苏省南京市鼓楼区南瑞路8号