发明名称 一种基于故障电阻非线性识别的高阻接地故障检测方法
摘要 一种基于故障电阻非线性识别的高阻接地故障检测方法,采集被监测馈线的相电压和零序电流瞬时值,分别计算各相电压与零序电流的相关系数,如果某相电压与零序电流的相关系数大于阈值R,则用最小二乘法分段线性拟合该相电压与零序电流构成的伏安特性曲线,进而根据各段拟合直线的斜率计算得到有效表征高阻接地故障特征的系数——故障电阻非线性系数,通过比较故障电阻非线性系数与阈值即可判断是否发生疑似高阻接地故障;连续检测疑似高阻接地故障的持续时间,如果超过延时,则最终判断为发生了高阻接地故障,本发明比基于谐波的检测方法灵敏度更高,比单纯基于时域量的检测方法具有更好的抗噪声能力。
申请公布号 CN103257302B 申请公布日期 2015.04.15
申请号 CN201310175572.9 申请日期 2013.05.13
申请人 清华大学 发明人 王宾;耿建昭;董新洲
分类号 G01R31/02(2006.01)I 主分类号 G01R31/02(2006.01)I
代理机构 西安智大知识产权代理事务所 61215 代理人 贾玉健
主权项 一种基于故障电阻非线性识别的高阻接地故障检测方法,其特征在于,包括以下步骤:步骤一、对被监测馈线的相电压和零序电流进行采样,获得一个工频周波的采样值序列u<sub>a</sub>(n)、u<sub>b</sub>(n)、u<sub>c</sub>(n)、i<sub>0</sub>(n),n=1,2,……60;步骤二、分别计算各相电压与零序电流的相关系数r<sub>x0</sub>,x=a,b,c,表示三相,如果存在r<sub>x0</sub>&gt;r<sub>set</sub>,r<sub>set</sub>为一整定常数,取值为0.966;则判断x相为疑似故障相,且令u<sub>f</sub>(n)=u<sub>x</sub>(n),n=1,2,……60;否则判断没有发生疑似高阻接地故障且重复步骤一,r<sub>x0</sub>的计算方法为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>r</mi><mrow><mi>x</mi><mn>0</mn></mrow></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mn>60</mn></munderover><msub><mi>i</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>u</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mn>60</mn></munderover><msubsup><mi>i</mi><mn>0</mn><mn>2</mn></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>u</mi><mi>x</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msqrt></mfrac></mrow>]]></math><img file="FDA0000654172640000011.GIF" wi="561" he="308" /></maths>步骤三、找到序列u<sub>f</sub>(n)的最大值点,记作u<sub>f</sub>(max),分别构建两组数据序列:x<sub>0</sub>(i)=i<sub>0</sub>(mod(max+42+i,60)),y<sub>0</sub>(i)=u<sub>f</sub>(mod(max+42+i,60)),i=1,2,……6x<sub>1</sub>(j)=i<sub>0</sub>(mod(max+55+j,60)),y<sub>1</sub>(j)=u<sub>f</sub>(mod(max+55+j,60)),j=1,2,……10其中mod(x,y)表示求x除以y的余数;步骤四、用最小二乘法分别线性拟合序列x<sub>0</sub>(i),y<sub>0</sub>(i)以及序列x<sub>1</sub>(j),y<sub>1</sub>(j),计算拟合直线的斜率k<sub>1</sub>、k<sub>2</sub>,并计算故障电阻非线性系数R,具体计算方法为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>k</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><mn>6</mn><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>6</mn></munderover><msub><mi>x</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><msub><mi>y</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>6</mn></munderover><msub><mi>x</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>6</mn></munderover><msub><mi>y</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow><mrow><mn>6</mn><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>6</mn></munderover><msubsup><mi>x</mi><mn>0</mn><mn>2</mn></msubsup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>6</mn></munderover><msub><mi>x</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>6</mn></munderover><msub><mi>x</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000654172640000012.GIF" wi="818" he="299" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>k</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><mn>10</mn><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>10</mn></munderover><msub><mi>x</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><msub><mi>y</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>10</mn></munderover><msub><mi>x</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>10</mn></munderover><msub><mi>y</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></mrow><mrow><mn>10</mn><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>10</mn></munderover><msubsup><mi>x</mi><mn>1</mn><mn>2</mn></msubsup><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>10</mn></munderover><msub><mi>x</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>10</mn></munderover><msub><mi>x</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000654172640000013.GIF" wi="874" he="307" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>R</mi><mo>=</mo><mfrac><msub><mi>k</mi><mn>1</mn></msub><msub><mi>k</mi><mn>2</mn></msub></mfrac></mrow>]]></math><img file="FDA0000654172640000021.GIF" wi="151" he="138" /></maths>步骤五、如果R&gt;R<sub>set</sub>,则判断为发生了高阻接地故障,其中R<sub>set</sub>为一整定常数,取值范围为1.2‑1.5;步骤六、每0.02秒重复一次步骤一至步骤五,得到每个工频周波是否发生高阻接地故障的结果;如果高阻接地故障持续时间超过阈值T<sub>set</sub>,则确定为发生了稳态高阻接地故障;其中T<sub>set</sub>取值为1秒。
地址 100084 北京市海淀区100084信箱82分箱清华大学专利办公室