发明名称 一种基于贝叶斯网络的变压器故障诊断分析方法
摘要 本发明公开了一种基于贝叶斯网络的变压器故障诊断分析方法,具体步骤为:步骤一:确定变压器油色谱的属性变量Y={Y<sub>1</sub>,.Y<sub>2</sub>,Y<sub>3</sub>…Y<sub>n</sub>}y<sub>i</sub>是Y<sub>i</sub>的取值和故障类型变量D={D<sub>1</sub>,D<sub>2</sub>,D<sub>3</sub>…D<sub>m</sub>}d<sub>m</sub>为D<sub>m</sub>的取值:以变压器油色谱的特征气体作为变压器油色谱的属性变量,以变压器出现的故障类型作为故障类型变量;步骤二:根据变压器油色谱的属性变量和故障类型变量确立贝叶斯网络分类模型及网络的结构参数和概率参数;步骤三:利用连接树算法确定变压器的故障类型。
申请公布号 CN103197177B 申请公布日期 2015.09.23
申请号 CN201310090961.1 申请日期 2013.03.20
申请人 山东电力集团公司济宁供电公司;国家电网公司 发明人 王彦良;王继文;陈晓红;张凡华;赵洪振;郑超;冯维华;王森;王卓;邓凸;王宏
分类号 G01R31/00(2006.01)I 主分类号 G01R31/00(2006.01)I
代理机构 济南圣达知识产权代理有限公司 37221 代理人 张勇
主权项 一种基于贝叶斯网络的变压器故障诊断分析方法,其特征是,具体步骤为:步骤一:确定变压器油色谱的属性变量Y={Y<sub>1</sub>,Y<sub>2</sub>,Y<sub>3</sub>,…,Y<sub>n</sub>}y<sub>i</sub>是Y<sub>i</sub>的取值和故障类型变量D={D<sub>1</sub>,D<sub>2</sub>,D<sub>3</sub>,…,D<sub>m</sub>}d<sub>m</sub>为D<sub>m</sub>的取值:以变压器油色谱的特征气体作为变压器油色谱的属性变量,以变压器出现的故障类型作为故障类型变量;步骤二:根据变压器油色谱的属性变量和故障类型变量确立贝叶斯网络分类模型及网络的结构参数和概率参数;步骤三:利用连接树算法确定变压器的故障类型;所述步骤二的具体步骤为:(2‑1)确定为建立网络模型有关的变量及其解释;(2‑2)建立一个表示条件独立断言的有向无环图,根据概率乘法公式有<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><mi>Y</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Pi;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mi>P</mi><mrow><mo>(</mo><msub><mi>Y</mi><mi>i</mi></msub><mo>/</mo><msub><mi>y</mi><mn>1</mn></msub><mo>,</mo><msub><mi>y</mi><mn>2</mn></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000725439640000011.GIF" wi="1139" he="197" /></maths>用P<sub>bm</sub>表示Y<sub>i</sub>的父节点集,则<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><msub><mi>Y</mi><mn>1</mn></msub><mo>,</mo><msub><mi>Y</mi><mn>2</mn></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>Y</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>Y</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>=</mo><munder><mi>&Pi;</mi><mi>i</mi></munder><mi>P</mi><mrow><mo>(</mo><msub><mi>Y</mi><mi>i</mi></msub><mo>/</mo><msub><mi>P</mi><mi>bm</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000725439640000012.GIF" wi="1274" he="158" /></maths>其中,n表示变压器油色谱属性变量的个数,i=1,2,…,n;(2‑3)为每一个属性变量Y<sub>i</sub>的父节点集的各个状态指派局部概率分布P(Y<sub>i</sub>/P<sub>bm</sub>);所述步骤三的具体步骤为:(3‑1)收集变压器属性变量的发生情况,根据贝叶斯定理,I<sub>i</sub>=(y<sub>1</sub>,y<sub>2</sub>,…,y<sub>n</sub>)属于故障类别变量D<sub>m</sub>的概率为<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><msub><mi>D</mi><mi>j</mi></msub><mo>/</mo><msub><mi>y</mi><mn>1</mn></msub><mo>,</mo><msub><mi>y</mi><mn>2</mn></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>y</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>P</mi><mrow><mo>(</mo><msub><mi>y</mi><mn>1</mn></msub><mo>,</mo><msub><mi>y</mi><mn>2</mn></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>y</mi><mi>n</mi></msub><mo>/</mo><msub><mi>D</mi><mi>j</mi></msub><mo>)</mo></mrow><mi>P</mi><mrow><mo>(</mo><msub><mi>D</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow><mrow><mi>P</mi><mrow><mo>(</mo><msub><mi>y</mi><mn>1</mn></msub><mo>,</mo><msub><mi>y</mi><mn>2</mn></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>y</mi><mi>n</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><mi>&beta;P</mi><mrow><mo>(</mo><msub><mi>D</mi><mi>j</mi></msub><mo>)</mo></mrow><mi>P</mi><mrow><mo>(</mo><msub><mi>y</mi><mn>1</mn></msub><mo>,</mo><msub><mi>y</mi><mn>2</mn></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>y</mi><mi>n</mi></msub><mo>/</mo><msub><mi>D</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000725439640000013.GIF" wi="1641" he="163" /></maths>式中:β是正则化因子,一般取1;P(D<sub>j</sub>)是D<sub>j</sub>的先验概率;P(D<sub>j</sub>/y<sub>1</sub>,y<sub>2</sub>,…y<sub>n</sub>)是D<sub>j</sub>的后验概率;后验概率反映了样本数据对D<sub>j</sub>的影响,将公式(3)表示为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><msub><mi>D</mi><mi>j</mi></msub><mo>/</mo><msub><mi>y</mi><mn>1</mn></msub><mo>,</mo><msub><mi>y</mi><mn>2</mn></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>y</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>&beta;P</mi><mrow><mo>(</mo><msub><mi>D</mi><mi>j</mi></msub><mo>)</mo></mrow><munderover><mi>&Pi;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mi>P</mi><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>/</mo><msub><mi>y</mi><mn>1</mn></msub><mo>,</mo><msub><mi>y</mi><mn>2</mn></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>D</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000725439640000021.GIF" wi="1240" he="136" /></maths>(3‑2)选择没有进行的但已经发生的具有最高概率的属性变量,并分析其父节点的状态;(3‑3)选择没有进行的且具有最高概率的父节点,根据其概率信息计算其对于故障的绝对概率信息,若该节点还有父节点,则按照步骤(3‑2)的方法执行直到无父节点为止;(3‑4)根据贝叶斯最大后验准则,公式(4)表示为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><msub><mi>D</mi><mi>j</mi></msub><mo>/</mo><msub><mi>y</mi><mn>1</mn></msub><mo>,</mo><msub><mi>y</mi><mn>2</mn></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>y</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>&beta;P</mi><mrow><mo>(</mo><msub><mi>D</mi><mi>j</mi></msub><mo>)</mo></mrow><munderover><mi>&Pi;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mi>P</mi><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>/</mo><msub><mi>D</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000725439640000022.GIF" wi="944" he="165" /></maths>根据公式(5),确定I<sub>i</sub>=(y<sub>1</sub>,y<sub>2</sub>,…,y<sub>n</sub>)属于故障类别变量D<sub>m</sub>的概率;(3‑5)根据概率信息确定故障原因,按故障原因概率大小输出结果,由最大概率解对变压器进行检修,排除故障。
地址 272129 山东省济宁市市中区火炬28号