(57)【要約】【課題】従来、精度上の問題や矛盾があった回転機固定子コイルの余寿命評価を、Weibull確率分布から直接絶縁破壊確率を求めることにより、高精度かつ合理的に行える装置を提供する。【解決手段】確率分布関数は、その履歴における絶縁耐力の低下因子を全て含んでいるとともに、当該履歴におけるBath Tub曲線上の関数になっているため、これから直接絶縁破壊確率を求めることによって精度上の問題と矛盾を解決する。解析装置１０は管理値回転機と任意の回転機固定子コイルの分布関数, 確率密度関数および絶縁耐力を算出し、これらの関数から絶縁破壊確率を求め、この算出値と管理値の裕度から余寿命を評価する。なお、前者の関数と絶縁耐力は設定手段１２で予め設定しておく。また、補助的に絶縁耐力の裕度を併用することによって評価をより明確にする。これらを作図装置１４で表示する。【選択図】図３
2008089550.tif 000002
【特許請求の範囲】
【請求項１】
回転機固定子コイルの余寿命を評価する装置であって、
管理値に相当する絶縁耐力に至った回転機（以下、管理値回転機という）の固定子コイルの絶縁破壊電圧のデータ（ｖ１ｍ，ｖ２ｍ，…，ｖｎｍ）からＭ．Ｌ．Ｅ．（Ｍａｘｉｍｕｍ  Ｌｉｋｅｌｉｈｏｏｄ  Ｅｓｔｉｍａｔｉｏｎ；最尤法）によってワイブル（Ｗｅｉｂｕｌｌ）分布のｍｍ（形状母数）とηｍ（尺度母数）を求め、ｍｍ，ηｍの値からＧ（ｖ）（分布関数），ｇ（ｖ）（確率密度関数）を決定して、設定手段により予め設定しておき、
次に評価対象となる任意の回転機の絶縁破壊電圧のデータ（ｖ１ｊ，ｖ２ｊ，…，ｖｎｊ）とｎｊ（コイル総本数）を与え、前記と同様にｖ１ｊ，ｖ２ｊ，…，ｖｎｊから、Ｍ．Ｌ．Ｅ．によってｍｊ（形状母数）とηｊ（尺度母数）を求め、ｍｊ，ηｊの値からＦｊ（ｖ）（分布関数），ｆｊ（ｖ）（確率密度関数）を決定し、
以上によって得られたｆｊ（ｖ），Ｇ（ｖ）（またはＦｊ（ｖ），ｇ（ｖ））からＰｒ．（Ｖｆｊ＜Ｖｍ）（絶縁破壊確率，Ｖｆｊは任意の回転機固定子コイルの絶縁破壊電圧のＷｅｉｂｕｌｌ分布であり、Ｖｍは管理値回転機固定子コイルの絶縁破壊電圧のＷｅｉｂｕｌｌ分布である）をモンテカルロシミュレーションによって近似解を求め、この値と絶縁破壊確率の管理値Ｐｒ．（Ｖｆｊ＜Ｖｍ）＝０．５との確率値の差によって余寿命評価をおこなうアルゴリズムをもつ回転機固定子コイルの余寿命評価装置。
【請求項２】
請求項１に記載の回転機固定子コイルの余寿命評価装置であって、
管理値回転機の、前記ｍｍ，ηｍとｎｍ（管理値回転機のコイル総本数）から求めたθｍ（管理値回転機の絶縁耐力，ここでいう絶縁耐力とは絶縁破壊電圧の最小値で、確率的にはコイル総本数分の１の値である）を設定手段により予め設定しておき、
次に、評価対象となる任意の回転機の、前記ｍｊ，ηｊとｎｊ（任意の回転機のコイル総本数）からθｊ（任意の回転機の絶縁耐力）を求め、
前記の絶縁破壊確率の近似解と絶縁破壊確率の管理値Ｐｒ．（Ｖｆｊ＜Ｖｍ）＝０．５との確率値の差と、θｊと予め設定してあるθｍとの差を併用して、余寿命評価をおこなうアルゴリズムをもつ回転機固定子コイルの余寿命評価装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、回転機固定子コイルの余寿命評価装置に関するものである。
【背景技術】
【０００２】
従来の回転機固定子コイル，例えば発電機固定子コイルの余寿命評価は、発電機の運転履歴（起動回数，運転時間，負荷変動など）に対する絶縁耐力の関係を直線回帰モデルで表し、このモデル上でおこなっていた（例えば特許文献１を参照）。
【特許文献１】特開２０００−２１４２１２号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
しかし、運転履歴の正確な把握が困難である他、絶縁耐力に影響を及ぼす因子には、運転履歴の他に絶縁材の材質の不均一，製作時の工作のバラツキなどが考えられ、これらを含めた絶縁耐力の低下現象を正確に表現することは困難であった（以下では運転履歴の他に絶縁耐力の低下に寄与する因子の全てを単に履歴と表現する）。また、絶縁耐力は初期，一様，加速劣化の各領域を経て低下する（つまり、Bath Tub曲線状に低下する）と考えられているが、回帰モデルは直線モデルで、運転履歴の経過如何に関わらず絶縁耐力が一定値で低下するという矛盾点があった。
【０００４】
本発明は、このような点に鑑みてなしたものであり、ある履歴をもち、管理値に相当する絶縁耐力に至った回転機，つまり管理値回転機固定子コイルの絶縁破壊電圧（ｖ）のWeibull分布と、任意の回転機固定子コイルのｖのWeibull分布から管理値に到達する絶縁破壊確率（Ｐｒ．（Ｖｆｊ＜Ｖｍ））を直接求める手法を開発し、そのアルゴリズムに基づいた余寿命評価装置を提供することを目的とする。
【課題を解決するための手段】
【０００５】
上記課題を解決するために、請求項１に記載したとおりの回転機固定子コイルの余寿命評価装置を構成した。ここで、請求項１でいう絶縁破壊確率とは、固定子コイルの絶縁耐力が管理値に到達する確率値で、管理値は、現在、日本国内で用いられている２Ｅ＋１（ＫＶ）（Ｅは定格電圧）とした。なお、以下では絶縁破壊電圧（ｖ）および絶縁耐力（θ）の単位は、絶縁耐力の初期値に対する割合（P.U.；Per Unit）で表す。
【０００６】
管理値回転機固定子コイルのWeibull分布関数は、その履歴における分布関数となっているため、当該履歴におけるBath Tub曲線上における分布関数を表していることになる。一方、任意の回転機固定子コイルのWeibull分布関数も、その履歴における分布関数となっていることから同様のことがいえる。管理値は、任意の回転機の確率密度関数（ｆｊ（ｖ））と管理値回転機の確率密度関数（ｇ（ｖ））が等しいときにＰｒ．（Ｖｆｊ＜Ｖｍ）＝０．５となるから、この値に対する確率値の差，つまり裕度（ｌａｄ）が確率値からみた余寿命となる（図４参照）。
【０００７】
また、請求項２に記載したとおりの回転機固定子コイルの余寿命評価装置を構成した。余寿命評価はｌａｄによっておこなうが、確率値であるため、やや判りにくい。このため、θｊと予め設定してあるθｍとの値の差，つまり、この裕度（ｌｂｄ）を併用することによって、これを補うことができる。これによって、任意の回転機の管理値に対する裕度が、より明確になる。管理値を超過している場合の、その度合いは同様に、それぞれｌａｕ，ｌｂｕで表すことができる。
【発明の効果】
【０００８】
本発明によれば、分布関数には、そのBath Tub曲線上の履歴における絶縁耐力の低下因子を全て含んでいるため、精度上の問題点と矛盾点を解消できることになる。このため、高精度かつ合理的な余寿命評価をおこなうことができる。
【発明を実施するための最良の形態】
【０００９】
本発明を実施するための最良の形態について、図１〜図５を参照しながら説明する。なお、本形態では回転機としてタービン発電機を適用した例を説明する。
【００１０】
まず、過去に固定子コイルの巻替やサンプリングにともなって絶縁破壊試験を実施して取得した１２台の絶縁破壊電圧のデータを用いた。これらの発電機の定格容量は１７０〜４４２（ＭＶＡ），定格電圧は１２〜２２（ＫＶ），コイル総本数は９６〜１４４（本）の間に分布している。解析対象の発電機固定子コイルは、絶縁材がエポキシレジン，コイルエンド部の支持方式がリジット型を対象にした。管理値（θｍ）は、ここで解析の対象にした発電機では、０．２９４〜０．３７５（P.U.）となる。
【００１１】
ここで、本発明を完成するにあたり、最適な確率分布を相関係数（Ｒ）によって決定した。絶縁破壊現象を表す分布として、極値III型分布であるWeibull分布 , 極値Ｉ型分布であるグンベル分布（Gumbel分布），正規分布（Normal分布），および対数正規分布（Log-Normal分布）を選定し、各分布について、ｖ１ｊ，ｖ２ｊ，…，ｖｎｊとＦｊ（ｖ）との間の相関係数（Ｒ）を求めた結果、図１のようになった。
【００１２】
図１において、横軸は１２台のタービン発電機毎の符号（Ａ，Ｂ，…，Ｌ）を、縦軸は相関係数（Ｒ）を表し、各発電機ごとにWeibull，Gumbel，Normal，Log-Normalの各分布を一組とした棒グラフを描いている。この結果から、Weibull分布が他の分布に比べて相関度合いが高いといえる。一般的に破壊現象はWeibull分布に従う場合が多いといわれているが、絶縁破壊現象も、またWeibull分布に従うということができる。従って、本発明では、タービン発電機の余寿命評価にあたってWeibull分布を用いることとした。
【００１３】
図３には、タービン発電機固定子コイルの余寿命評価装置に相当する解析装置の構成例をブロック図で模式的に表す。図３に表す解析装置１０は、設定手段１２や作図装置１４などを有する。解析装置１０は、ＣＰＵがプログラムを実行することでデータの入出力や計算式の演算等を実現する。設定手段１２はWeibull分布関数（Ｇ（ｖ）），確率密度関数（ｇ（ｖ））および絶縁耐力（θｍ）を設定可能な手段、例えばＲＯＭやＲＡＭ等の記憶媒体が該当する。作図装置１４は表示器やプリンタ等の出力機器が該当する。
【００１４】
解析装置１０では、以下のような処理をおこなう。本装置内でのWeibull分布は２個のパラメーターをもつ２母数（ｍ；形状母数，η；尺度母数）Weibull分布として取り扱った。まず、管理値に相当する絶縁耐力に至ったタービン発電機（以下では「管理値発電機」という）の固定子コイルの絶縁破壊電圧のデータ（ｖ１ｍ，ｖ２ｍ，…，ｖｎｍ；ｎｍはコイル総本数）からｍｍ，ηｍをM.L.E.（Maximum Likelihood Estimation；最尤法）で求めた。管理値発電機の絶縁耐力（θｍ）はｍｍ，ηｍを用いて絶縁破壊電圧の最低値、つまり、コイル総本数（ｎｍ）分の１の確率値として定義し、次の式（１）で求める。
【数１】
2008089550.tif 000003
【００１５】
次に、評価対象となる任意のタービン発電機（以下では「任意の発電機」という）の各発電機固定子コイル毎に、絶縁破壊電圧のデータ（ｖ１ｊ，ｖ２ｊ，…，ｖｎｊ；ｎｊはコイル総本数）から、各発電機毎にｍｊ，ηｊをM.L.E.で求めた。任意の発電機の絶縁耐力（θｊ）は、ｍｊ，ηｊを用いて絶縁破壊電圧の最低値、つまり、コイル総本数（ｎｊ）分の１の確率値と定義し、次の式（２）で求める。
【数２】
2008089550.tif 000004
【００１６】
任意の発電機固定子コイルの絶縁破壊電圧のWeibull分布をＶｆｊ，管理値発電機固定子コイルの絶縁破壊電圧のWeibull分布をＶｍとすると、これらは互いに独立であるから、管理値に到達する絶縁破壊確率Ｐｒ．（Ｖｆｊ＜Ｖｍ）を次の式（３）で求める。なお、絶縁破壊確率の管理値はｆｊ（ｖ）＝ｇ（ｖ）のときにＰｒ．（Ｖｆｊ＜Ｖｍ）＝０．５となる。
【００１７】
【数３】
2008089550.tif 000005
【００１８】
ここで、ｆｊ（ｖ），ｇ（ｖ）は、それぞれ任意の発電機，管理値発電機固定子コイルの絶縁破壊電圧のWeibull確率密度関数であり、また、Ｆｊ（ｖ），Ｇ（ｖ）は、それぞれ任意の発電機，管理値発電機固定子コイルの絶縁破壊電圧のWeibull分布関数である。なお、ｖは絶縁破壊電圧である。これらはｍ，ηによって次の式（４）のように表される。
【数４】
2008089550.tif 000006
【００１９】
ここで、サヒックスのｊは、任意の発電機のｊ番目のタービン発電機を，サヒックスのｍは管理値発電機を表している。
【００２０】
Ｐｒ．（Ｖｆｊ＜Ｖｍ）の概念図を図２に示す。ここで、Ｖｆｊは任意の発電機（以下、ｊ番目の発電機という）固定子コイルの絶縁破壊電圧のWeibull分布であり、Ｖｍは管理値発電機固定子コイルの絶縁破壊電圧のWeibull分布である。Ｖｆｊは、ｍｊ（形状母数）とηｊ（尺度母数）で表される２母数Weibull分布であり、また、Ｖｍは、ｍｍ（形状母数）とηｍ（尺度母数）で表される２母数Weibull分布である。Ｖｍを表す２母数（ｍｍ，ηｍ）は、管理値である２Ｅ＋１（ＫＶ）に、ほぼ合致する値を示した２台のタービン発電機の絶縁破壊電圧のデータからM.L.E.で求めた平均値を採用し、ｍｍ＝７．２３，ηｍ＝０．６５８とした。なお、上記の式（３）は、ｆｊ（ｖ）＝ｇ（ｖ）のときに理論解が得られ、絶縁破壊確率の管理値Ｐｒ．（Ｖｆｊ＜Ｖｍ）＝０．５となる。しかし、ｆｊ（ｖ）≠ｇ（ｖ）のときは理論解が得られないため、モンテカルロシミュレーションによって近似解を求める。
【００２１】
上述したアルゴリズムによって開発した余寿命評価装置を図３に示す。本装置はｊ番目の発電機固定子コイルの絶縁破壊電圧のデータ（ｖ１ｊ，ｖ２ｊ，…，ｖｎｊ）とコイル総本数（ｎｊ）を与えることによって、絶縁破壊確率Ｐｒ．（Ｖｆｊ＜Ｖｍ）とθｊを得ることができ、これによって、ｊ番目の発電機のθｊに対するＰｒ．（Ｖｆｊ＜Ｖｍ）の値を求めることができる。ここで、Ｖｆｊはｊ番目の発電機固定子コイルの絶縁破壊電圧のWeibull分布を，また、Ｖｍは管理値発電機固定子コイルの絶縁破壊電圧のWeibull分布である。余寿命評価の例を図４に示す。絶縁破壊確率の管理値はＰｒ．（Ｖｆｊ＜Ｖｍ）＝０．５であるから、この値に対する確率値の差，つまり裕度（ｌａｄ）が確率値からみた余寿命となる。余寿命評価はｌａｄによっておこなうが、確率値であるため、やや判りにくい。このため、予め設定してあるθｍと算出値θｊとの値の差，つまり、この裕度（ｌｂｄ）を併用することによって、これを補うことができる。以上のことから、ｊ番目の発電機の管理値に対する裕度が明確になる。管理値を超過している場合の、その度合いは同様に、それぞれｌａｕ，ｌｂｕで表すことができる。
【００２２】
図３において、管理値発電機固定子コイルの絶縁破壊電圧のデータ（ｖ１ｍ，ｖ２ｍ，…，ｖｎｍ）は既知であるから、そのデータからWeibull分布の形状母数（ｍｍ），尺度母数（ηｍ）をM.L.E.によって算出し、それらの値（ｍｍ，ηｍ）から、Weibull分布関数（Ｇ（ｖ）），確率密度関数（ｇ（ｖ））を決定するとともに、コイル総本数（ｎｍ）から、コイル総本数（ｎｍ）分の１の確率値として絶縁耐力（θｍ）を求める。これらの値（ｍｍ，ηｍ，Ｇ（ｖ），ｇ（ｖ），θｍ）は、設定手段１２によって予め装置内に設定しておく。
【００２３】
一方、任意の発電機，つまりｊ番目の発電機固定子コイルの絶縁破壊電圧のデータ（ｖ１ｊ，ｖ２ｊ，…，ｖｎｊ）が得られたとき、装置にそのデータを与えることによって、Weibull分布の形状母数（ｍｊ）および尺度母数（ηｊ）をM.L.E.で求める。それらの値（ｍｊ，ηｊ）からWeibull分布関数（Ｆｊ（ｖ））および確率密度関数（ｆｊ（ｖ））を決定するとともに、コイル総本数（ｎｊ）を装置に与えることによって、絶縁耐力（θｊ）をコイル総本数（ｎｊ）分の１の確率値として求める。
【００２４】
以上の分布関数と確率密度関数から、絶縁破壊確率を前述した式（３）で求める。
【００２５】
ここで、Ｖｆｊはｊ番目の発電機固定子コイルの絶縁破壊電圧のWeibull分布であり、また、Ｖｍは管理値発電機固定子コイルの絶縁破壊電圧のWeibull分布である。Ｐｒ．（Ｖｆｊ＜Ｖｍ）を模式的に表すと図２のようになる。上記の式（５）はｆｊ（ｖ）＝ｇ（ｖ）のときに限って理論解が得られ、その値は０．５となる。しかし、それ以外の時には理論解が求まらないため、モンテカルロシミュレーションによって近似解を得る。
【００２６】
以上の結果から、θｊに対するＰｒ．（Ｖｆｊ＜Ｖｍ）の関係を作図装置１４で図示すれば、ｊ番目の発電機の絶縁破壊確率が求まることになる。これから求められる値と絶縁破壊確率の管理値Ｐｒ.（Ｖｆｊ＜Ｖｍ）＝０．５との確率値の差が余寿命である。また、同図中にθｍを併記すれば、これとθｊとの値の差も余寿命とみることができるため、これを併用することによって管理値に対する裕度が、より明確になる。前述の１２台の発電機の絶縁耐力（θｊ）に対する絶縁破壊確率値Ｐｒ．（Ｖｆｊ＜Ｖｍ）を「○」印でプロットすると図５のようになった。同図では、１２台分のタービン発電機の絶縁耐力に対する絶縁破壊確率値を、まとめて表示しているため、θｍが巾（ｄ）をもっている。この巾ｄは０．２９４〜０．３７５（P.U.）の範囲にある。
【００２７】
同図からθｊが０．６（P.U.）前後の値になると管理値に到達する確率値が増え始め、管理値に到達する絶縁耐力の値は、ほぼ０．３５（P.U.）になるという結果が得られた。なお、図４，５の縦軸のＰｒ．（Ｖｆｊ＜Ｖｍ）の値は０〜１となっているが、確率値であるため、０に限りなく近い０より大きい値〜１に限りなく近い１より小さい値を示している。
【００２８】
本装置は、あるタービン発電機固定子コイルの余寿命評価を、絶縁破壊電圧のデータとコイル総本数を与えることによって、管理値に達する度合い，あるいは管理値を超過する度合いが簡単に，また高精度かつ合理的に求められる。この度合いが裕度であり、この値によって任意の発電機固定子コイルの余寿命評価をおこなうことができる。
【００２９】
上述した実施の形態では、回転機としてタービン発電機を適用したが、他の発電機や高圧電動機の固定子コイルの余寿命評価装置に適用することが可能である。他の発電機としては、タービン発電機以外の同期発電機や誘導発電機などが該当する。これらの回転機固定子コイルについても、精度上の問題点と矛盾点を解消できることになり、高精度かつ合理的な余寿命評価をおこなうことができる。
【００３０】
本余寿命評価装置は、高精度かつ合理的な余寿命評価がおこなえるため、固定子コイルの巻替を最適な時期におこなうことができる。その結果、絶縁破壊の防止と巻替に伴う不要な支出を防ぎ、経済性の確保に寄与できる。特に絶縁破壊が地絡から短絡に移行した場合は、発電機やタービン等の機器が損壊する可能性もあるが、巻替のタイミングを的確に把握できるため、機器の損壊も防止できる。
【図面の簡単な説明】
【００３１】
【図１】確率分布の相関度合いを比較した図である。
【図２】絶縁破壊確率の概念を説明する図である。
【図３】余寿命評価装置の構成例を模式的に表すブロック図である。
【図４】余寿命評価を説明する図である。
【図５】絶縁耐力と絶縁破壊確率値との関係を説明する図である。
【符号の説明】
【００３２】
１０  解析装置（回転機固定子コイルの余寿命評価装置）
１２  設定手段
１４  作図装置
【図１】
2008089550.tif 000007
【図４】
2008089550.tif 000008
【図５】
2008089550.tif 000009
【図２】
2008089550.tif 000010
【図３】
2008089550.tif 000011