发明名称 一种基于故障暂态行波衰减分量的电网故障区段定位方法
摘要 本发明公开了一种基于故障暂态行波衰减分量的电网故障区段定位方法。该方法在电网已配置测量装置情况下,通过对故障后信号的小波变换获取小波系数模极大值作为表征初始行波幅值的特征量;特征量最大的测量点对应于故障线路的一端,再利用故障行波经最短路径传播时特征量的衰减特性确定故障线路的另一端,判断故障区段。该区段定位方法准确度高,对测量装置的同步性没有特殊要求,具有较好地经济性和较强的实用价值。
申请公布号 CN106443347A 申请公布日期 2017.02.22
申请号 CN201610881840.2 申请日期 2016.10.09
申请人 中国矿业大学 发明人 梁睿;薛雪;杨智;王飞;彭楠;刘成磊
分类号 G01R31/08(2006.01)I 主分类号 G01R31/08(2006.01)I
代理机构 南京瑞弘专利商标事务所(普通合伙) 32249 代理人 杨晓玲
主权项 一种基于故障暂态行波衰减分量的电网故障区段定位方法,其特征在于,包括如下步骤:1)对于拓扑结构确定的电网,假设共有N个测量点,分别标记为m<sub>1</sub>,m<sub>2</sub>,…,m<sub>N</sub>,并构成集合M={m<sub>1</sub>,m<sub>2</sub>,…,m<sub>N</sub>},根据拓扑结构得到每个测量点的相邻测量点数目分别为j<sub>1</sub>,j<sub>2</sub>,…,j<sub>N</sub>,并构成集合J={j<sub>1</sub>,j<sub>2</sub>,…,j<sub>N</sub>};2)故障发生后,对各测量点获取的线模电压分别进行小波变换,提取小波系数模极大值表征该测量点的行波幅值,各测量点的行波幅值分别为f<sub>1</sub>,f<sub>2</sub>,…f<sub>N</sub>,并构成集合F={f<sub>1</sub>,f<sub>2</sub>,…f<sub>N</sub>};测量点m<sub>a</sub>的相邻测量点行波幅值记为集合{f<sub>a1</sub>,f<sub>a2</sub>,...f<sub>aja</sub>},a∈{1,2,...N},得到各测量点的相邻测量点行波幅值矩阵L表示为:<img file="FDA0001127093990000011.GIF" wi="582" he="326" />其中,Z为集合J中的最大值,L中的空元素用数值零来替代;3)步骤2)所获得最大行波幅值对应的测量点确定为故障的一端,该测量点记为m<sub>i</sub>,并得到测量点m<sub>i</sub>的P个相邻测量点分别记为b<sub>1</sub>,b<sub>2</sub>,…b<sub>P</sub>,并构成集合B={b<sub>1</sub>,b<sub>2</sub>,…b<sub>P</sub>};4)选择集合B中一个相邻测量点b<sub>j</sub>,相邻测量点b<sub>j</sub>对应的行波幅值记为f<sub>bj</sub>,j∈{1,2,...P};获取b<sub>j</sub>到一个测量点m<sub>n</sub>的最短路径所经过的节点,所经过的节点根据b<sub>j</sub>到测量点m<sub>n</sub>依次经过节点的序号排列为集合D={m<sub>x</sub>,…m<sub>y</sub>},其中m<sub>x</sub>为所经过的第一个节点的序号,m<sub>y</sub>为所经过的最后一个节点的序号,n∈{1,2,...N};若所述集合D中没有测量点m<sub>i</sub>,集合D中各节点对应的相邻测量点数目集合为JD={j<sub>x</sub>,…j<sub>y</sub>},则集合D中每个节点对应的折射系数为d={d<sub>mx</sub>,…d<sub>my</sub>}={2/(j<sub>x</sub>+1),…2/(j<sub>y</sub>+1)},由b<sub>j</sub>折算到测量点m<sub>n</sub>的行波幅值为f<sub>nj</sub>’=f<sub>bj</sub>×d<sub>mx</sub>×…×d<sub>my</sub>;若所述集合D中包含有测量点m<sub>i</sub>,删除集合D中m<sub>i</sub>之前的节点,所述集合D更新为D={m<sub>i</sub>,…m<sub>y</sub>},则由b<sub>j</sub>折算到测量点m<sub>i</sub>的行波幅值为f<sub>nj</sub>’=f<sub>mi</sub>×…×d<sub>my</sub>,f<sub>mi</sub>为测量点m<sub>i</sub>对应的行波幅值;初始n=1,j=1,遍历集合B中所有相邻测量点,得到由P个相邻测量点折算到所有测量点的行波幅值矩阵为G为:<img file="FDA0001127093990000021.GIF" wi="613" he="303" />5)将测量m<sub>a</sub>的行波幅值分别与其每一个相邻测量点行波幅值作差,得到集合<img file="FDA0001127093990000025.GIF" wi="597" he="87" />a∈{1,2,...N},遍历所述集合F中每个测量点得到矩阵△F:<img file="FDA0001127093990000022.GIF" wi="917" he="319" />根据行波幅值在传播过程中的衰减特性,矩阵ΔF中各元素作差后数值的正负表示行波的传播方向;6)对折算出来的行波幅值矩阵G中的第一列做如下处理,得到折算出来的传播方向矩阵△G<sub>1</sub>为:<img file="FDA0001127093990000023.GIF" wi="989" he="327" />用符号函数sgn(x)对矩阵△G<sub>1</sub>和△F进行标准化处理:<maths num="0001"><math><![CDATA[<mrow><mi>sgn</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mrow><mi>x</mi><mo>&gt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>x</mi><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mn>1</mn></mrow></mtd><mtd><mrow><mi>x</mi><mo>&lt;</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001127093990000024.GIF" wi="407" he="222" /></maths>其相同元素个数表征两个矩阵的相似程度,并记相同元素个数为s<sub>1</sub>,依次对G中的各列做上述处理,得到相同元素个数集合为S={s<sub>1</sub>,s<sub>2</sub>,…s<sub>p</sub>},集合S中的最大值对应的测量点即为故障线路的另一端,最终确定故障区段。
地址 221116 江苏省徐州市铜山区大学路中国矿业大学科研院