主权项 |
一种基于故障电阻非线性识别的高阻接地故障检测方法,其特征在于,包括以下步骤:步骤一、对被监测馈线的相电压和零序电流进行采样,获得一个工频周波的采样值序列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>>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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>R<sub>set</sub>,则判断为发生了高阻接地故障,其中R<sub>set</sub>为一整定常数,取值范围为1.2‑1.5;步骤六、每0.02秒重复一次步骤一至步骤五,得到每个工频周波是否发生高阻接地故障的结果;如果高阻接地故障持续时间超过阈值T<sub>set</sub>,则确定为发生了稳态高阻接地故障;其中T<sub>set</sub>取值为1秒。 |