发明名称 变压器差动保护的励磁涌流模糊识别方法
摘要 本发明公开了一种变压器差动保护的励磁涌流模糊识别方法,基于模糊理论的2可加模糊测度法,能够做到各种涌流判据的优势互补,综合波形间断角判据、波形对称判据、低电压判据,构成判据集,通过加权系数来调整各个判据对励磁涌流判据的影响。即使某一个判据不完善而出错,只会导致模糊量隶属度函数偏移,而不会导致最终结果的错误,解决了判据单一而使得错误动作相对高频率发生的问题,提高了变压器保护的性能。
申请公布号 CN103545789B 申请公布日期 2016.12.28
申请号 CN201310377571.2 申请日期 2013.08.26
申请人 江苏科技大学 发明人 曾庆军;卢梦宁;王彪
分类号 H02H7/045(2006.01)I;G01R19/00(2006.01)I 主分类号 H02H7/045(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 楼高潮
主权项 一种变压器差动保护的励磁涌流模糊识别方法,其特征在于,该方法包含下列步骤:1)选择波形间断角判据、波形对称判据、低电压判据构成判据集;①波形间断角判据隶属函数的选取令:<maths num="0001"><math><![CDATA[<mrow><msub><mi>K</mi><mi>i</mi></msub><mo>=</mo><mfrac><msub><mi>S</mi><mi>i</mi></msub><mrow><mi>S</mi><mo>/</mo><mn>2</mn></mrow></mfrac></mrow>]]></math><img file="FDA0001004714920000011.GIF" wi="190" he="135" /></maths>其中,S为电流波形一个周期的面积,S<sub>i</sub>为半个周期面积;通过移动半周期数据窗得n/2个面积值,再根据模糊贴进度原理,运用海明贴进度,建立判据:<maths num="0002"><math><![CDATA[<mrow><mi>N</mi><mo>=</mo><mn>1</mn><mo>-</mo><mfrac><mn>1</mn><mrow><mi>n</mi><mo>/</mo><mn>2</mn></mrow></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>n</mi><mo>/</mo><mn>2</mn></mrow></munderover><mrow><mo>|</mo><mrow><mn>1</mn><mo>-</mo><msub><mi>K</mi><mi>i</mi></msub></mrow><mo>|</mo></mrow></mrow>]]></math><img file="FDA0001004714920000012.GIF" wi="428" he="135" /></maths>其中n为电流波形一周期采样点数,故障时N值与1接近,涌流时远离1;选取梯形隶属函数,表达式为:<maths num="0003"><math><![CDATA[<mrow><msub><mi>f</mi><mn>1</mn></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mn>0.........................</mn><mi>N</mi><mo>&lt;</mo><msub><mi>c</mi><mn>1</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mn>1</mn><mrow><msub><mi>c</mi><mn>2</mn></msub><mo>-</mo><msub><mi>c</mi><mn>1</mn></msub></mrow></mfrac><mrow><mo>(</mo><mi>N</mi><mo>-</mo><msub><mi>c</mi><mn>1</mn></msub><mo>)</mo></mrow><mn>......</mn><msub><mi>c</mi><mn>1</mn></msub><mo>&lt;</mo><mi>N</mi><mo>&lt;</mo><msub><mi>c</mi><mn>2</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><mn>1.........................</mn><mi>N</mi><mo>&gt;</mo><msub><mi>c</mi><mn>2</mn></msub></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001004714920000013.GIF" wi="702" he="302" /></maths>其中c<sub>1</sub>=0.15,c<sub>2</sub>=0.24;②波形对称判据隶属函数的选取波形对称判据是将流入差动继电器的差流进行微分,将微分后的差流的前半波与后半波做对称比较来区分故障电流和励磁涌流,该对称度定义式如下:<maths num="0004"><math><![CDATA[<mrow><msub><mi>K</mi><mrow><mi>s</mi><mi>y</mi><mi>m</mi></mrow></msub><mo>=</mo><mrow><mo>|</mo><mfrac><mrow><msubsup><mi>I</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>+</mo><msup><msub><mi>I</mi><mrow><mi>i</mi><mo>+</mo><mi>M</mi><mo>/</mo><mn>2</mn></mrow></msub><mo>&prime;</mo></msup></mrow><mrow><msubsup><mi>I</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>-</mo><msup><msub><mi>I</mi><mrow><mi>i</mi><mo>+</mo><mi>M</mi><mo>/</mo><mn>2</mn></mrow></msub><mo>&prime;</mo></msup></mrow></mfrac><mo>|</mo></mrow><mo>&le;</mo><mi>K</mi></mrow>]]></math><img file="FDA0001004714920000014.GIF" wi="444" he="164" /></maths>其中,I<sub>i</sub>'为差电流导数前半波第i点的采样值,i为前半波任一采样点,I<sub>i+M/2</sub>'为差电流导数后半波对应第i点的采样值,M为一个周期的采样点数,k为比较阈值取0.85;当i点满足上式时,就称其为对称,否则为不对称;现取一周期采样点数M为20,选择正态形隶属函数,则具体表达式为:<maths num="0005"><math><![CDATA[<mrow><msub><mi>f</mi><mn>2</mn></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msup><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>a</mi><mo>)</mo></mrow><mn>2</mn></msup><mi>b</mi></mfrac></mrow></msup><mn>..............4</mn><mo>&le;</mo><mi>n</mi><mo>&le;</mo><mn>10</mn></mrow></mtd></mtr><mtr><mtd><mrow><mn>0...........................</mn><mi>n</mi><mo>&lt;</mo><mn>4</mn></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001004714920000015.GIF" wi="629" he="190" /></maths>其中,n为不满足上述对称度定义公式的采样点数,a=10,b=51.9;③低电压判据隶属函数的选取隶属函数选取梯形隶属函数为:<maths num="0006"><math><![CDATA[<mrow><msub><mi>f</mi><mn>3</mn></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mn>1...........................</mn><mfrac><mi>U</mi><msub><mi>U</mi><mi>N</mi></msub></mfrac><mo>&lt;</mo><msub><mi>c</mi><mn>1</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mn>1</mn><mrow><msub><mi>c</mi><mn>1</mn></msub><mo>-</mo><msub><mi>c</mi><mn>2</mn></msub></mrow></mfrac><mrow><mo>(</mo><mfrac><mi>U</mi><msub><mi>U</mi><mi>N</mi></msub></mfrac><mo>-</mo><msub><mi>c</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>......</mn><msub><mi>L</mi><mn>1</mn></msub><mo>&lt;</mo><mfrac><mi>U</mi><msub><mi>U</mi><mi>N</mi></msub></mfrac><mo>&lt;</mo><msub><mi>c</mi><mn>2</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><mn>0...........................</mn><mfrac><mi>U</mi><msub><mi>U</mi><mi>N</mi></msub></mfrac><mo>&gt;</mo><msub><mi>c</mi><mn>2</mn></msub></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001004714920000021.GIF" wi="782" he="455" /></maths>其中,U为变压器端电压,U<sub>N</sub>为额定电压,取c<sub>1</sub>=0.2,c<sub>2</sub>=0.7;2)选取N组差电流数据,平均分成两组数值:N<sub>1</sub>组为励磁涌流电流,N<sub>2</sub>组为故障电流,N组差电流代入各判据的隶属函数表达式得函数值f<sub>i</sub>,并代入公式:<maths num="0007"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mo>(</mo><mi>c</mi><mo>)</mo></mrow></mtd><mtd><mrow><mo>&Integral;</mo><mi>f</mi><mi>d</mi><mi>&mu;</mi><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>-</mo><msub><mi>f</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mi>&mu;</mi><mrow><mo>(</mo><mo>{</mo><msub><mi>a</mi><mi>j</mi></msub><mo>/</mo><msub><mi>f</mi><mi>j</mi></msub><mo>&GreaterEqual;</mo><msub><mi>f</mi><mi>i</mi></msub><mo>}</mo><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001004714920000022.GIF" wi="854" he="127" /></maths>其中,设判据集A=(a<sub>1</sub>,a<sub>2</sub>,……,a<sub>m</sub>),f是定义在有限集合A上的实值非负函数,令f<sub>i</sub>表示f(a<sub>i</sub>),a<sub>i</sub>是判据集A=(a<sub>1</sub>,a<sub>2</sub>,……,a<sub>m</sub>)中的子集或元素;f<sub>j</sub>表示f(a<sub>j</sub>),a<sub>j</sub>是判据集A=(a<sub>1</sub>,a<sub>2</sub>,……,a<sub>m</sub>)中的子集或元素;μ为定义在P(A)上的模糊测度,其中f<sub>i</sub>≥f<sub>i‑1</sub>,f<sub>0</sub>=0;μ({a<sub>j</sub>/f<sub>j</sub>≥f<sub>i</sub>})含义是μ(a<sub>j</sub>),表示判据子集a<sub>j</sub>的重要性或者权重;3)求解下述非线性优化模型从而得单点集和两点集上的2可加模糊测度,非线性优化模型:<maths num="0008"><math><![CDATA[<mrow><mi>m</mi><mi>a</mi><mi>x</mi><mfrac><msup><mrow><mo>(</mo><msub><mover><mi>y</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mo>-</mo><msub><mover><mi>y</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mn>1</mn></msub></munderover><msup><mrow><mo>(</mo><msub><mi>y</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>-</mo><msub><mover><mi>y</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mn>2</mn></msub></munderover><msup><mrow><mo>(</mo><msub><mi>y</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mo>-</mo><msub><mover><mi>y</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac></mrow>]]></math><img file="FDA0001004714920000023.GIF" wi="677" he="215" /></maths><maths num="0009"><math><![CDATA[<mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>y</mi><mrow><mn>2</mn><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>y</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>+</mo><mi>&Delta;</mi><mi>y</mi><mo>&le;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><munder><mo>&Sigma;</mo><mrow><msub><mi>a</mi><mi>j</mi></msub><mo>&Element;</mo><mi>k</mi></mrow></munder><msub><mi>&mu;</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>-</mo><munder><mo>&Sigma;</mo><mrow><msub><mi>a</mi><mi>j</mi></msub><mo>&Element;</mo><mi>k</mi></mrow></munder><msub><mi>&mu;</mi><mi>j</mi></msub><mo>-</mo><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>2</mn><mo>)</mo></mrow><msub><mi>&mu;</mi><mi>i</mi></msub><mo>&GreaterEqual;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mo>&ForAll;</mo><msub><mi>a</mi><mi>i</mi></msub><mo>&Element;</mo><mi>A</mi><mo>,</mo><mi>K</mi><mo>&Subset;</mo><mi>A</mi><mo>/</mo><msub><mi>a</mi><mi>i</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mi>&Delta;</mi><mi>y</mi><mo>&gt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><munder><mo>&Sigma;</mo><mrow><mo>{</mo><msub><mi>a</mi><mi>i</mi></msub><mo>,</mo><msub><mi>a</mi><mi>j</mi></msub><mo>}</mo><mo>&Subset;</mo><mi>A</mi></mrow></munder><msub><mi>&mu;</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>-</mo><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>2</mn><mo>)</mo></mrow><munder><mo>&Sigma;</mo><mrow><msub><mi>a</mi><mi>i</mi></msub><mo>&Element;</mo><mi>A</mi></mrow></munder><msub><mi>&mu;</mi><mi>i</mi></msub><mo>=</mo><mn>1</mn></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001004714920000024.GIF" wi="670" he="494" /></maths>其中,Δy表示变化量;μ<sub>i</sub>表示μ(a<sub>i</sub>);μ<sub>j</sub>表示μ(a<sub>j</sub>);μ<sub>ij</sub>表示μ(a<sub>i</sub>,a<sub>j</sub>),y<sub>1i</sub>、y<sub>2i</sub>为第i组励磁涌流和故障电流的模糊积分值,全部励磁涌流和故障电流的模糊积分结果的平均值:<img file="FDA0001004714920000031.GIF" wi="350" he="143" />j=1,2通过公式<img file="FDA0001004714920000032.GIF" wi="638" he="118" />计算全部子集上的测度值;4)将变压器差电流对各判据的隶属函数值及各判据的重要度代入公式:<img file="FDA0001004714920000033.GIF" wi="861" he="123" />得模糊积分,比较该值与定值0.7的大小;5)前者大,则为励磁涌流;反之,为故障电流。
地址 212003 江苏省镇江市梦溪路2号