发明名称 考虑实时健康状态的电力变压器可靠性评估方法
摘要 本发明公开了一种考虑实时健康状态的电力变压器可靠性评估的方法。本发明包括如下步骤:步骤1、确定当前时刻、最后状态监测时刻、最后状态转移点、未来考察时间点、电力变压器状态逗留时间的初始分布等条件;步骤2、使用接受拒绝采样方法生成随机数,即每个状态的状态逗留时间;步骤3、计算本次采样样本下的可靠度与平均剩余寿命;步骤4、根据方差系数判断收敛条件是否满足;若不满足,返回第2步继续新的采样;若满足,计算可靠度函数与平均剩余寿命,评估完毕。与现有方法相比,本发明提出的模型中假设变压器处于周期监测中,但是健康状态的转移可以发生在任何时刻,这比现有的大多数模型更加符合实际。本发明可靠、易行,便于推广。
申请公布号 CN106569052A 申请公布日期 2017.04.19
申请号 CN201610890465.8 申请日期 2016.10.11
申请人 国网湖北省电力公司;浙江大学;国家电网公司 发明人 李大虎;孙建波;楼贤嗣;何宇斌;郭创新
分类号 G01R31/00(2006.01)I 主分类号 G01R31/00(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 叶志坚
主权项 考虑实时健康状态的电力变压器可靠性评估方法,其特征在于包括如下步骤:步骤1、确定条件,包括确定当前时刻、最后状态监测时刻、最后状态转移点、未来考察时间点、电力变压器状态逗留时间的初始分布;所述的当前时刻是进行电力变压器停运率计算的时刻t<sub>0</sub>;所述的最后状态监测时刻是最后一次状态监测时刻Δ<sub>m</sub>;所述的未来考察时间点是计算电力变压器停运率的目标时刻;所述的电力变压器状态逗留时间的初始分布分为两种情况;第一种情况为当前时刻t<sub>0</sub>是状态监测时刻,则当前时刻的健康状态已知,变压器健康状态通过在线监测或离线检测直接获取;在这种情况下,t<sub>0</sub>是状态监测时刻,假设t<sub>0</sub>时刻电力变压器的状态Z(t<sub>0</sub>)=j,则状态j的逗留时间X<sub>j</sub>需要满足X<sub>j</sub>>t<sub>0</sub>‑T<sub>j‑1</sub>;其中T<sub>j‑1</sub>为上一次变压器健康状态改变的时刻;状态j的逗留时间须按照条件分布进行采样,之后的逗留时间分布g(y)为原始分布;第二种情况为t<sub>0</sub>不是状态监测时刻,则当前时刻的健康状态未知,此时对应为:变压器健康状态在之前状态监测时刻已经得到,下一个在线监测或离线检测时刻尚未达到;在这种情况下,t<sub>0</sub>不是状态监测时刻,则距离t<sub>0</sub>最近的上一次状态监测时刻为Δ<sub>m</sub>,则状态j的逗留时间X<sub>j</sub>需满足X<sub>j</sub>>Δ<sub>m</sub>‑T<sub>j‑1</sub>;状态j的逗留时间分布为条件分布,之后的逗留时间分布g(y)为原始分布;此时,t<sub>0</sub>时刻不一定处于状态j,由实际采样样本确定;步骤2、通过接受拒绝采样方法生成每个状态的状态逗留时间;所述的接受拒绝采样方法步骤如下:假设G(y)是待采样的分布,F(y)是一个已知的能够直接进行采样的简单分布,则:①产生分布为F(y)的随机变量y;②生成一个(0,1)均匀分布的随机数U;③如果<img file="FDA0001128691940000021.GIF" wi="223" he="110" />则接收y作为分布G(y)的一个样本;否则,返回第①步;其中,系数c应确保<img file="FDA0001128691940000022.GIF" wi="346" he="111" />步骤3、计算本次采样样本下的可靠度与平均剩余寿命;所述的可靠度计算方法如下:若已知变压器在t<sub>0</sub>时刻正常运行,则根据可靠性数学理论可知,在未来时刻t的可靠度函数为<maths num="0001"><math><![CDATA[<mrow><mi>R</mi><mrow><mo>(</mo><mi>t</mi><mo>|</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><mi>P</mi><mrow><mo>(</mo><mi>T</mi><mo>&GreaterEqual;</mo><mi>t</mi><mo>|</mo><mi>T</mi><mo>&gt;</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><msubsup><mo>&Integral;</mo><msub><mi>t</mi><mn>0</mn></msub><mi>t</mi></msubsup><mi>&lambda;</mi><mo>(</mo><mrow><mi>s</mi><mo>,</mo><mi>Z</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow><mo>)</mo><mi>d</mi><mi>s</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001128691940000023.GIF" wi="845" he="103" /></maths>其中λ(s,Z(s))=λ<sub>0</sub>(s)·Ψ(Z(s))(2)式中,T是故障时刻,λ<sub>0</sub>(s)是基准函数,用来描述设备基本的老化过程,模型中采用了威布尔分布的基准函数<img file="FDA0001128691940000024.GIF" wi="144" he="116" />Ψ(Z(s))是连接函数,用来量化协变量Z(s)对强迫停运率的影响,模型中法采用了指数分布的连接函数e<sup>γZ(s)</sup>;Z(s)称为时刻s的协变量状态,即变压器实时健康状态,由变压器DGA的周期性监测数据得出;设备强迫停运率λ的函数表达式随Z(s)变化而变化,当变压器健康状态Z(s)发生状态转移时,强迫停运率λ(s,Z(s))发生跃变,因此可靠度函数(1)能够看作一个分段函数,分段的时间节点即为状态转移时刻;在之前采样得到的确定样本下,根据公式(1)能够计算出当前样本下未来任意时刻的可靠度函数;所述的平均剩余寿命计算方法如下:平均剩余寿命评估了设备从当前时刻起,在未来时间内的期望无故障工作时间,可通过下式计算<maths num="0002"><math><![CDATA[<mrow><mi>M</mi><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><mi>E</mi><mrow><mo>(</mo><mi>T</mi><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>|</mo><mi>T</mi><mo>&gt;</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><msubsup><mo>&Integral;</mo><msub><mi>t</mi><mn>0</mn></msub><mi>&infin;</mi></msubsup><mi>R</mi><mrow><mo>(</mo><mi>t</mi><mo>|</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mi>d</mi><mi>t</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001128691940000025.GIF" wi="823" he="108" /></maths>步骤4、根据方差系数判断收敛条件是否满足;若不满足,返回步骤2继续新的采样;若满足,计算可靠度函数与平均剩余寿命,评估完毕;所述的方差系数计算公式为:<maths num="0003"><math><![CDATA[<mrow><mi>c</mi><mi>v</mi><mo>=</mo><mfrac><msqrt><mrow><mfrac><mn>1</mn><mrow><mi>n</mi><mi>s</mi><mrow><mo>(</mo><mi>n</mi><mi>s</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><msubsup><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>n</mi><mi>s</mi></mrow></msubsup><mrow><mo>(</mo><msub><mi>C</mi><mi>k</mi></msub><mo>-</mo><mover><mi>C</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow></mrow></msqrt><mover><mi>C</mi><mo>&OverBar;</mo></mover></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001128691940000031.GIF" wi="677" he="151" /></maths>其中,ns指代采样的样本数量;C<sub>k</sub>指代第k个样本的平均剩余寿命;<img file="FDA0001128691940000032.GIF" wi="42" he="62" />指代平均剩余寿命的期望;若cv<0.05,则认为整体计算过程收敛;若不满足,则返回步骤2继续新的采样,进入下一次计算,直到最终结果满足收敛标准;变压器最终的可靠度函数与平均剩余寿命是所有样本下的均值。
地址 430077 湖北省武汉市洪山区徐东大街197号