发明名称 一种故障录波数据故障起点检测方法
摘要 本发明涉及一种故障录波数据故障起点检测方法,尤其是涉及一种利用新型突变量电流检测和最大似然估计检测相结合的故障起点双重检测方法。本发明根据CONTRADE标准对故障录波文件解码,得到电压、电流采样序列,利用新型突变量电流检测方法得到初始突变点,然后用最大似然估计检测对初始突变点进行二次检测,如果满足条件,则该突变点即为故障起点;如果不满足条件,则对初始突变点的后一个点进行最大似然估计检测,直到检测出满足条件的突变点为止。该方法能够准确检测出采样序列的故障起点,检测灵敏度高,几乎不受系统频率变化的影响,而且可以有效抵制噪声的影响,故障起点检测精度高,结果可靠。
申请公布号 CN105510745A 申请公布日期 2016.04.20
申请号 CN201510988411.0 申请日期 2015.12.24
申请人 武汉大学 发明人 龚庆武;魏东;占劲松;林燕贞;梁志远;杨安琪
分类号 G01R31/00(2006.01)I 主分类号 G01R31/00(2006.01)I
代理机构 武汉科皓知识产权代理事务所(特殊普通合伙) 42222 代理人 严彦
主权项 一种故障录波数据故障起点检测方法,其特征在于:包括以下步骤,步骤1、提取故障录波数据的采样序列;步骤2、检测得到初始突变点,记为K,实现方式如下,计算突变量电流,公式如下,<math><![CDATA[<mrow><mi>&Delta;</mi><mi>i</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>N</mi><mn>1</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>|</mo><mi>i</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mi>i</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><msub><mi>N</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>|</mo><mo>,</mo></mrow></mtd><mtd><mrow><mi>k</mi><mo>=</mo><msub><mi>N</mi><mn>1</mn></msub><mo>+</mo><mn>1</mn><mo>,</mo><msub><mi>N</mi><mn>1</mn></msub><mo>+</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mn>2</mn><msub><mi>N</mi><mn>1</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>|</mo><mo>|</mo><mi>i</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mi>i</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><msub><mi>N</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>|</mo><mo>-</mo><mo>|</mo><mi>i</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><msub><mi>N</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><mi>i</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>2</mn><msub><mi>N</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>|</mo><mo>|</mo><mo>,</mo></mrow></mtd><mtd><mrow><mi>k</mi><mo>=</mo><mn>2</mn><msub><mi>N</mi><mn>1</mn></msub><mo>+</mo><mn>1</mn><mo>,</mo><mn>2</mn><msub><mi>N</mi><mn>1</mn></msub><mo>+</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>N</mi></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000889293430000011.GIF" wi="1382" he="247" /></maths>其中,N为采样点总数,N<sub>1</sub>是一个工频周期的采样点数,Δi(k)表示在k点对应采样时刻的突变量电流;i(k)表示在k点对应采样时刻的测量电流采样值;i(k‑N<sub>1</sub>)表示在k点对应采样时刻前一周期的采样电流值,i(k‑2N<sub>1</sub>)表示在k点对应采样时刻前二周期的采样电流值;将满足以下两个判据的第一个点k作为被检测故障数据的初始突变起点,记为K,判决1,突变量电流占同相正常电流的百分比大于整定值ξ,记作|Ψ(k)|>ξ,其中突变量电流占同相正常电流的百分比Ψ(k)=Δi(k)/i(k‑N<sub>1</sub>);判决2,时间窗中满足判据(1)的元素个数大于等于整定值β,记作card(A<sub>k</sub>)≥β,其中,集合A<sub>k</sub>={Ψ(n)||Ψ(n)|>ξ,k≤n<k+α},card(A<sub>k</sub>)为集合A<sub>k</sub>的元素的个数,α为时间窗宽度;步骤3、根据最大似然估计检测方法确定故障起点,实现方式如下,设定用z(k)描述“k是故障起点”这一事件可能性的大小,z(k)的值介于0和1之间;最大似然检测判据设为,如果z(k)>0.5,则k点是故障起点;如果z(k)≤0.5,则k不是故障起点;执行以下过程,(1)令当前的点d=K,即先对初始突变点K进行验证;(2)对当前的点d进行验证,根据最大似然检测判据,针对当前的点d,若z(d)>0.5,参数Xn取值为1,则当前的点d是故障起点,进入步骤4;若z(d)≤0.5,参数Xn取值为0,d不是故障起点,继续执行(3);(3)令d=d+1,返回执行(2)对下一个采样点进行验证,直到检测出满足条件的点为止,进入步骤4;步骤4、根据步骤3所确定故障起点,取故障起点之前的采样序列为故障前的电气量,故障起点及以后的采样序列为故障后的电气量。
地址 430072 湖北省武汉市武昌区珞珈山武汉大学