"実際の変圧器の構成に沿う回路モデルにおける共振周波数に対する前記回路モデルのパラメータを一つずつ個別に変化させたときの共振周波数の変化量を計算し、共振周波数の変化の程度が他と比べて大きいパラメータに関係する素子若しくは素子の特定部位に異常が発生していると判断することを特徴とする変圧器の異常診断方法。前記実際の変圧器の構成に沿う回路において電源と回路のインピーダンスとの間に抵抗を直列に追加挿入し、前記抵抗が十分に小さいときの状態方程式と前記抵抗が十分に大きいときの状態方程式とを導出し、これら状態方程式の係数行列の固有値を計算することによって前記実際の変圧器の構成に沿う回路モデルにおける共振周波数を計算することを特徴とする請求項1記載の変圧器の異常診断方法。実際の変圧器の構成に沿う回路モデルにおける共振周波数に対する前記回路モデルのパラメータを一つずつ個別に変化させたときの共振周波数の変化量を計算する手段と、共振周波数の変化の程度が他と比べて大きいパラメータに関係する素子若しくは素子の特定部位に異常が発生していると判断する手段とを有することを特徴とする変圧器の異常診断装置。前記実際の変圧器の構成に沿う回路において電源と回路のインピーダンスとの間に抵抗を直列に追加挿入し、前記抵抗が十分に小さいときの状態方程式と前記抵抗が十分に大きいときの状態方程式とを導出し、これら状態方程式の係数行列の固有値を計算することによって前記実際の変圧器の構成に沿う回路モデルにおける共振周波数を計算することを特徴とする請求項3記載の変圧器の異常診断装置。実際の変圧器の構成に沿う回路モデルにおける共振周波数に対する前記回路モデルのパラメータを一つずつ個別に変化させたときの共振周波数の変化量を計算する手段、共振周波数の変化の程度が他と比べて大きいパラメータに関係する素子若しくは素子の特定部位に異常が発生していると判断する手段としてコンピュータを機能させるための変圧器の異常診断プログラム。前記実際の変圧器の構成に沿う回路において電源と回路のインピーダンスとの間に抵抗を直列に追加挿入し、前記抵抗が十分に小さいときの状態方程式と前記抵抗が十分に大きいときの状態方程式とを導出し、これら状態方程式の係数行列の固有値を計算することによって前記実際の変圧器の構成に沿う回路モデルにおける共振周波数を計算することを特徴とする請求項5記載の変圧器の異常診断プログラム。しかしながら、周波数応答解析を用いた従来の診断では、変圧器における巻線の異常の発生を検出することはできても、検出した異常即ち共振周波数の変化がどのような巻線異常或いは鉄心異常を意味しているのか診断する言い換えると、特定することはできない。このため、周波数応答解析を用いての変圧器の異常の具体的な原因言い換えると、変圧器の異常の様相を診断する指針の検討が十分とは言えず、周波数応答解析を用いた従来の診断では有用性が十分であるとは言い難い。そこで、本発明は、変圧器における異常の具体的な原因異常の様相を特定することができる変圧器の異常診断方法、異常診断装置及び異常診断プログラムを提供することを目的とする。本発明者は、変圧器における異常の診断手法の検討を行う中で、計算に用いる回路モデルのパラメータを一つずつ個別に変化させて得られる共振周波数の変化の程度変化量変化率に基づいて変圧器における異常の具体的な原因異常の様相を特定できることを突き止めた。請求項1記載の変圧器の異常診断方法は、前記知見に基づくものであり、実際の変圧器の構成に沿う回路モデルにおける共振周波数に対する回路モデルのパラメータを一つずつ個別に変化させたときの共振周波数の変化量を計算し、共振周波数の変化の程度が他と比べて大きいパラメータに関係する素子若しくは素子の特定部位に異常が発生していると判断するようにしている。また、請求項3記載の変圧器の異常診断装置は、実際の変圧器の構成に沿う回路モデルにおける共振周波数に対する回路モデルのパラメータを一つずつ個別に変化させたときの共振周波数の変化量を計算する手段と、共振周波数の変化の程度が他と比べて大きいパラメータに関係する素子若しくは素子の特定部位に異常が発生していると判断する手段とを有するようにしている。また、請求項5記載の変圧器の異常診断プログラムは、実際の変圧器の構成に沿う回路モデルにおける共振周波数に対する回路モデルのパラメータを一つずつ個別に変化させたときの共振周波数の変化量を計算する手段、共振周波数の変化の程度が他と比べて大きいパラメータに関係する素子若しくは素子の特定部位に異常が発生していると判断する手段としてコンピュータを機能させるようにしている。これら変圧器の異常診断方法、異常診断装置及び異常診断プログラムによると、診断対象である実際の変圧器の構成に沿う回路であって回路モデルのパラメータを変化させていない回路の共振周波数に対する回路モデルのパラメータを一つずつ個別に変化させたときの共振周波数の変化量を計算し、共振周波数の変化の程度が他と比べて大きいパラメータを選定するようにしているので、変圧器において異常が発生している素子や素子の部位が特定される。なお、本発明は、周波数応答解析を実施可能な変圧器であれば、例えば電源変圧器や配電変圧器などの比較的大型の変圧器並びに柱上変圧器や更に小型の変圧器などあらゆるサイズ用途の変圧器を対象として適用可能である。また、請求項2,請求項4,請求項6記載の発明は、実際の変圧器の構成に沿う回路において電源と回路のインピーダンスとの間に抵抗を直列に追加挿入し、抵抗が十分に小さいときの状態方程式と抵抗が十分に大きいときの状態方程式とを導出し、これら状態方程式の係数行列の固有値を計算することによって実際の変圧器の構成に沿う回路モデルにおける共振周波数を計算するようにしている。この場合には、実際の変圧器の構成に沿う回路モデルに仮想の抵抗を挿入すると共に当該仮想の抵抗が十分に大きいときと十分に小さいときとの二種類の回路モデル及び状態方程式を用いることにより、伝達関数を計算することなく変圧器回路の共振周波数のみが計算されるので、実際の変圧器の構成に沿う回路モデルにおける共振周波数が非常に簡便に計算される。本発明の変圧器の異常診断方法、異常診断装置及び異常診断プログラムによれば、変圧器において異常が発生している素子や素子の部位を特定することができるので、変圧器の異常診断としての有用性の向上を図ることが可能になる。さらに、請求項2,請求項4,請求項6に記載の変圧器の異常診断方法、異常診断装置及び異常診断プログラムによれば、伝達関数を計算することなく変圧器回路の共振周波数のみを計算することができ、共振周波数を非常に簡便に計算することができるので、計算処理にかかる時間を大幅に短縮して変圧器の異常診断としての適用容易性と有用性との向上を図ることが可能になる。以下、本発明の構成を図面に示す実施の形態の一例に基づいて詳細に説明する。図1及び図2に、本発明の変圧器の異常診断方法、異常診断装置及び異常診断プログラムの実施形態の一例を示す。本実施形態の変圧器の異常診断方法は、変圧器の健全時と診断時との共振周波数を測定するS1と共にこれら共振周波数を比較してS2変化がある変圧器に対しS3:Yes、実際の変圧器の構成に沿う回路モデルにおける共振周波数を計算しS4,S5,S6、当該共振周波数に対する回路モデルのパラメータを一つずつ個別に変化させたときの共振周波数の変化量を計算しS7、共振周波数の変化の程度が他と比べて大きいパラメータを選定しS8、当該パラメータに関係する素子若しくは素子の特定部位に異常が発生していると判断するS9ようにしている。また、本実施形態の変圧器の異常診断方法は、実際の変圧器の構成に沿う回路において電源と回路のインピーダンスとの間に抵抗を直列に追加挿入しS4、抵抗が十分に小さいときの状態方程式と抵抗が十分に大きいときの状態方程式とを導出しS5、これら状態方程式の係数行列の固有値を計算することによって実際の変圧器の構成に沿う回路モデルにおける共振周波数を計算するS6ようにしている。そして、上述の変圧器の異常診断方法は、本発明の変圧器の異常診断装置として実現される。本実施形態の変圧器の異常診断装置は、実際の変圧器の構成に沿う回路モデルにおける共振周波数に対する回路モデルのパラメータを一つずつ個別に変化させたときの共振周波数の変化量を計算する手段と、共振周波数の変化の程度が他と比べて大きいパラメータに関係する素子若しくは素子の特定部位に異常が発生していると判断する手段とを有するようにしている。また、上述の変圧器の異常診断方法及び変圧器の異常診断装置は、本発明の変圧器の異常診断プログラムをコンピュータ上で実行することによっても実現される。本実施形態では、変圧器の異常診断プログラムをコンピュータ上で実行することによって変圧器の異常診断装置が実現されると共に変圧器の異常診断方法が実施される場合を例に挙げて説明する。変圧器の異常診断プログラム17を実行するためのコンピュータ10本実施形態では、変圧器の異常診断装置10でもあるの全体構成を図2に示す。この変圧器の異常診断装置10は、制御部11、記憶部12、入力部13、表示部14及びメモリ15を備え相互にバス等の信号回線によって接続されている。また、変圧器の異常診断装置10にはデータサーバ16がバス等の信号回線によって接続されており、その信号回線を介してデータや制御指令等の信号の送受信即ち出入力が相互に行われる。制御部11は記憶部12に記憶されている変圧器の異常診断プログラム17によって変圧器の異常診断装置10全体の制御並びに変圧器の異常診断処理に係る演算を行うものであり、例えばCPUである。記憶部12は少なくともデータやプログラムを記憶可能な装置であり、例えばハードディスクである。メモリ15は制御部11が種々の制御や演算を実行する際の作業領域であるメモリ空間となるものであり、例えばRAMである。入力部13は少なくとも作業者の命令を制御部11に与えるためのインターフェイスであり、例えばキーボードである。表示部14は制御部11の制御によって文字や図形等の描画表示を行うものであり、例えばディスプレイである。そして、記憶装置としてのデータサーバ16にアクセス可能なコンピュータである変圧器の異常診断装置10の制御部11には、変圧器の異常診断プログラム17を実行することにより、実際の変圧器の構成に沿う回路モデルにおける共振周波数を計算する手段回路モデル作成部11a,状態方程式導出部11b,共振周波数計算部11c、前記共振周波数に対する回路モデルのパラメータを一つずつ個別に変化させたときの共振周波数の変化量を計算する手段感度解析部11d、共振周波数の変化の程度が他と比べて大きいパラメータを選定する手段パラメータ選定部11e、前記パラメータに関係する素子若しくは素子の特定部位に異常が発生していると判断する手段異常様相診断部11fが構成される。また、本実施形態の変圧器の異常診断プログラム17を実行することにより、制御部11には、実際の変圧器の構成に沿う回路において電源と回路のインピーダンスとの間に抵抗を直列に追加挿入する手段回路モデル作成部11a、抵抗が十分に小さいときの状態方程式と抵抗が十分に大きいときの状態方程式とを導出する手段状態方程式導出部11b、これら状態方程式の係数行列の固有値を計算することによって実際の変圧器の構成に沿う回路モデルにおける共振周波数を計算する手段共振周波数計算部11cが構成される。本発明の変圧器の異常診断方法の実施に際しては、まず、異常診断対象の変圧器言い換えると、実機の共振周波数の測定が行われるS1。S1の処理における共振周波数の測定は、変圧器が健全である時における測定を予め行っておき、また、変圧器の診断時において行う。S1の処理における、異常診断対象の変圧器実機の共振周波数を測定する方法は、特定の方法に限られるものではなく、変圧器の共振周波数が測定解析可能であればいずれの方法でも良い。具体的には例えば周波数応答解析が用いられ得る。本実施形態では、例えばとして、健全時において第一共振周波数XkHz,第二共振周波数YkHz,第三共振周波数ZkHzの三つの共振周波数が測定され、診断時において第一共振周波数XkHz,第二共振周波数YkHz,第三共振周波数ZkHzの三つの共振周波数が測定されたとする。次に、健全時における共振周波数と診断時における共振周波数との比較が行われるS2。本実施形態では、第一共振周波数XとXとが比較され,第二共振周波数YとYとが比較され,第三共振周波数ZとZとが比較される。そして、健全時における共振周波数と診断時における共振周波数との間に変化がない場合S3:Noには、変圧器に異常は発生していないと判断されて異常の診断処理は終了する言い換えると、本発明による異常様相の診断は行われない;END。一方、健全時における共振周波数と診断時における共振周波数との間に変化がある場合S3:Yesには、変圧器に異常が発生していると判断されて本発明による異常様相の診断としてS4以降の処理が行われる。なお、変圧器に異常が発生していると判断するための健全時と診断時との間の共振周波数の変化の程度は、特定の値に限定されるものではなく、従来から用いられている判断基準を参考にして適宜設定される。ここで、健全時における共振周波数と診断時における共振周波数との間に変化があってS3:Yes本発明による異常様相の診断の対象になった変圧器については、変化が見られた共振周波数に関する健全時における共振周波数の値が共振周波数データベース18としてデータサーバ16に保存される。なお、本発明において共振周波数値データが蓄積されるデータベースやデータファイルが保存されるのは、データサーバ16に限られるものではなく、記憶部12に保存されるようにしても良いし、光記憶媒体等の各種記憶媒体や外部記憶装置に保存されるようにしても良い。本実施形態では、健全時における第一共振周波数XkHzと診断時における第一共振周波数XkHzとの間に変化が見られたとする。そして、本実施形態では、健全時における第一共振周波数XkHzが共振周波数データベース18としてデータサーバ16に保存される。なお、以下においては、健全時と診断時との間で変化が見られた共振周波数のことを着目共振周波数と呼ぶ。以上のS1からS3までの処理により、異常が発生している変圧器が抽出されると共に当該変圧器についての着目共振周波数の値データが整備される。これらS1からS3までの処理は、本発明による異常の様相の特定にとって前段階の処理であり、本発明の要部を構成するものではなく、従来若しくは新規のどのような方法によって行われても構わない。一方、S4以降は本発明の変圧器の異常診断の要部を構成する処理である。また、本実施形態の変圧器の異常診断プログラム17をコンピュータ上で実行することによって実現される変圧器の異常診断装置10はS4以降に対応する処理を実行する。異常の様相の特定に係る具体的な処理として、S3の処理において異常が発生していると判断されて本発明による異常様相の診断の対象になった変圧器についての回路モデルの作成が行われるS4。回路モデル等価回路とも呼ばれるの作成自体は従来の方式に従って作成される。一方で、本発明においては、後記のS5及びS6の処理の説明において詳述するように、解析対象である実際の変圧器の回路構成に抵抗を直列に追加挿入した回路モデルを用いる。本実施形態では、制御部11の回路モデル作成部11aが、異常様相の診断の対象になった変圧器についての回路モデルのデータを読み込み、若しくはデータの入力を受け、この実際の変圧器の構成に沿う回路モデルに対して電源1と回路のインピーダンスZTrとの間に抵抗Rを直列に追加挿入した回路モデル図3参照のデータをメモリ15に記憶させる。なお、異常様相の診断の対象になった変圧器についての回路モデルのデータの読み込みや入力は、特定の方式に限定されるものではなく、例えば、変圧器の異常診断プログラム17の処理形式に合わせたデータファイルとして予め整備されて記憶部12やデータサーバ16等の記憶装置に保存されたものを読み込むようにしても良いし、汎用モジュールやコンポーネントなどとして整備されているGUIを利用して入力として与えられるものを受けるようにしても良い。ここで、回路モデルが与えられれば、汎用されている回路解析ソフトウェアを用いて伝達関数を計算して共振周波数を抽出することが可能である。しかしながら、本発明においては共振周波数の計算を非常に多数回繰り返す必要があるので、伝達関数を計算して共振周波数を抽出する方法では計算処理に長時間を要する虞があって現実的ではない。そこで、本実施形態では、伝達関数を計算することなく、共振周波数のみを計算する。具体的には、本実施形態では、伝達関数を計算することなく変圧器回路の共振周波数のみを計算するため、状態方程式の導出を行いS5、共振周波数の計算を行うS6ようにする。回路の方程式を導出する定式化手法として状態変数法がある。状態変数法は、数式1-1で表される状態微分方程式と数式1-2で表される入出力状態方程式とで回路を表現する手法である。なお、状態変数法自体は周知の技術であるのでここでは詳細については省略する例えば、RARohrer著CircuitTheoryAnIntroductiontotheStateVariableApproach,MacGrawHillBookCompany,1971年、また、この文献の和訳として、斉藤正男,篠崎寿夫共訳回路理論状態変数解析入門,学献社,1973年を参照。そして、数式1-1における係数行列Aの固有値の虚部は回路の直列共振周波数になることが知られている例えば、前掲CircuitTheoryAnIntroductiontotheStateVariableApproach,回路理論状態変数解析入門を参照。次に、回路の並列共振周波数の計算について検討する。解析対象である実際の変圧器の回路構成を図3に示す回路モデルで表して状態微分方程式の係数行列Aから直列共振周波数が計算されるところ、回路のインピーダンスZTrを数式2で表すとこれは零点を計算していることと等価である。ここに、ZTr:回路のインピーダンス,s:ラプラス演算子,z:零点,p:極,添字m:零点の番号,添字n:極の番号をそれぞれ表す。ここで、図3に示すように電源1と回路のインピーダンスZTrとの間に抵抗Rを直列に追加挿入した回路モデルを考えると、このときの回路全体のインピーダンス[ZTrR]は数式3で表される。そして、数式3の零点に着目するために分子のみを考える。抵抗Rが十分に大きいとき、数式3の分子は数式4で表される。数式4から、回路に高抵抗を直列に挿入した場合の零点はもとの回路の極と一致することが分かる。すなわち、解析対象である実際の変圧器の構成に沿う回路において電源と回路のインピーダンスとの間に高抵抗を直列に追加挿入して導出した状態方程式から、もとの回路の並列共振周波数を計算できることが考察される。そして、本実施形態では、本発明者による上述の考察と変圧器異常診断手法の検討において得られた知見とに基づいて、まず、異常様相診断対象の実際の変圧器の回路モデルに仮想の抵抗Rを直列に追加挿入した回路モデルを用い、仮想抵抗Rが十分に小さいときの状態方程式と十分に大きいときの状態方程式とを導出してそれぞれの係数行列Aの固有値を計算する。ここで、係数行列Aの固有値は、回路モデルのパラメータの個数に対応して複数個計算される。そして、本発明者による変圧器異常診断手法の検討において得られた知見より、仮想抵抗Rが十分に小さいときと仮想抵抗Rが十分に大きいときとの間で虚部に変化がない固有値虚部は共振周波数は実際の回路言い換えると、実際の変圧器に対する測定においては共振周波数としては出現しないものであって変圧器の異常様相の診断には関係がないので、仮想抵抗Rが十分に小さいときと十分に大きいときとの間で虚部に変化があった固有値のみを用いて共振周波数を算出する。なお、仮想抵抗Rについて十分に小さい値と十分に大きい値とは、特定の値に限定されるものではないが、具体的には例えば十分に小さい値としてR=1mmΩを用いると共に十分に大きい値としてR=1010Ωを用いることが考えられる。本実施形態では、状態方程式の導出S5として、制御部11の状態方程式導出部11bが、S4の処理においてメモリ15に記憶された仮想抵抗Rを含む回路モデルのデータをメモリ15から読み込み、当該データを用いて仮想抵抗Rが十分に小さいときの状態方程式及び仮想抵抗Rが十分に大きいときの状態方程式を導出する。状態方程式の導出は、グラフ理論を用いてキルヒホッフの定理を適用することで行う。このグラフ理論を用いてのキルヒホッフの定理の適用による状態方程式の導出自体は周知の技術であるので例えば、RARohrer著CircuitTheoryAnIntroductiontotheStateVariableApproach,MacGrawHillBookCompany,1971年、また、この文献の和訳として、斉藤正男,篠崎寿夫共訳回路理論状態変数解析入門,学献社,1973年を参照、ここでは詳細については省略するが、概要は以下の通りである。状態方程式導出部11bが、上述の通り回路モデルのデータをメモリ15から読み込み、各素子両端に接点番号を与える。このとき、接続されている素子には同じ接点番号を与えることで回路の接続/非接続を表現する。なお、グラフ理論では、素子を枝という。回路網にグラフ理論を適用して木と補木とに分割する。なお、木はグラフの全ての接点を接続して閉路を含まない枝の集合であり、補木は全体から木を取り除いた残りである。そして、木に対して基本カットセットと呼ばれる枝の部分集合を抽出し、カットセットに含まれる枝の情報を電流が流れる方向も考慮して1,0,-1を要素としたカットセット-枝行列で表現し、このカットセット-枝行列から状態方程式を導出する。そして、状態方程式導出部11bは、仮想抵抗Rが十分に小さいときと十分に大きいときとの二種類の状態方程式をメモリ15に記憶させる。次に、共振周波数の計算S6として、制御部11の共振周波数計算部11cが、S5の処理においてメモリ15に記憶された仮想抵抗Rが十分に小さいときの状態方程式をメモリ15から読み込んで係数行列Aの固有値を計算すると共に、S5の処理においてメモリ15に記憶された仮想抵抗Rが十分に大きいときの状態方程式をメモリ15から読み込んで係数行列Aの固有値を計算する。そして、仮想抵抗Rが十分に小さいときと十分に大きいときとの間で虚部が同一でない固有値のみについて固有値の虚部から回路の共振周波数を算出する。そして、共振周波数計算部11cは、算出した共振周波数の値をメモリ15に記憶させる。さらに、共振周波数計算部11cは、共振周波数の計算過程において算定される固有値及び固有ベクトルをメモリ15に記憶させる。次に、回路モデルの各パラメータを個別に微小変化させたときそれぞれの共振周波数の変化の程度に基づいて共振周波数の感度解析が行われるS7。S7の処理は、着目共振周波数について行う。すなわち、回路モデルの各パラメータを個別に微小変化させたときそれぞれにおける着目共振周波数の変化の程度に基づいて感度解析を行う。本発明における回路モデルのパラメータとは、回路モデルを構成する種々の素子の性能や仕様などに関する値である。具体的には例えば、変圧器についての回路モデルの説明の例としての図4に示す回路モデルにおける、1次巻線,2次巻線を表すインダクタンスL1,L2や、その上下端に接続されている巻線間静電容量を二分割した静電容量C1,C4や、巻線の上下端に二分割して接続されている1次巻線の対地静電容量C3,C6上端の対地静電容量C3にはブッシングの静電容量を加えるようにしても良いや、巻線の上下端に二分割して接続されている2次巻線の対地静電容量C2,C5上端の対地静電容量C2にはブッシングの静電容量を加えるようにしても良いなどが挙げられる。なお、図4に示すものはあくまでも説明のための非常に簡単な例であり、変圧器についての回路モデルのパラメータは実際には数百程度になることが考えられる。そして、本発明における回路モデルのパラメータは、上記の例に限られるものではなく、異常様相診断対象の実際の変圧器についての回路モデルを構成する種々の素子の性能や仕様などに関する値が適宜抽出される。S7の処理においては、回路モデルのパラメータを一つずつ個別に微小に変化させ、パラメータを一つ微小に変化させた状態毎に共振周波数を計算する。なお、パラメータを変化させる程度は、特定の値に限定されるものではないが、具体的には例えば1~3程度にすることが考えられる。S7の処理における感度解析は、上述のように回路モデルのパラメータを一つずつ個別に変化させて計算される共振周波数と、もとの回路即ち、パラメータを変化させていない回路について計算される共振周波数本実施形態では、S6の処理において計算されるとを対比することによって行う。したがって、回路モデルのパラメータが例えば三百個である場合には共振周波数の計算を三百回行うことによってパラメータ毎の変化後の共振周波数を計算することができるものの、このように共振周波数を個別に繰り返し計算する方法では計算処理に長時間を要する虞がある。そこで、本実施形態では、S6の処理における共振周波数の計算過程で算定される固有値と固有ベクトルとを用いて回路モデルの各パラメータが微小変化したときの共振周波数の変化量を計算することによって感度解析を行う。具体的には、パラメータαが微小変化したときの固有値の変化量∂λ/∂αを表す数式5を用いて固有値感度解析を行う。ここに、λ:固有値虚部は共振周波数である,α:回路モデル中の着目するパラメータ,p:固有ベクトル,q:固有ベクトル,A:状態行列をそれぞれ表す。本実施形態では、制御部11の感度解析部11dが、S6の処理においてメモリ15に記憶された固有値及び固有ベクトルをメモリ15から読み込むと共に、回路モデルの各パラメータが微小変化したときの固有値虚部が共振周波数であるの変化量∂λ/∂αを数式5によって計算する。感度解析部11dは、また、データサーバ16内の共振周波数データベース18から着目共振周波数の値を読み込むと共に、S6の処理において計算されメモリ15に記憶された共振周波数であって着目共振周波数に対応する共振周波数の値をメモリ15から読み込む。そして、感度解析部11dは、着目共振周波数について、回路モデルのパラメータ毎に、もとの回路即ち、実際の変圧器の構成に沿う回路であって回路モデルのパラメータを変化させていない回路の共振周波数に対する回路モデルのパラメータを一つずつ個別に変化させたときの共振周波数の変化量を整理する。なお、回路モデルのパラメータの種類言い換えると、S7の処理において個別に微小変化させるパラメータの指定及びパラメータを変化させる程度変化量変化率は、変圧器の異常診断プログラム17内に予め規定しておくようにしても良いし、微小変化させるパラメータの種類やパラメータの変化量変化率の指定を要求するメッセージをS7の処理の段階で表示部14に表示させると共に入力部13を介して入力されたものを感度解析部11dに与えるようにしても良いし、微小変化させるパラメータの種類やパラメータの変化量変化率を記録したデータファイルを記憶部12やデータサーバ16等の記憶装置に予め保存しておくと共に当該データファイルから感度解析部11dが読み込むようにしても良い。そして、感度解析部11dは、着目共振周波数についての、回路モデルのパラメータ毎の周波数変化量をメモリ15に記憶させる。次に、着目共振周波数における周波数変化に対して微小変化の影響が大きいパラメータの選定が行われるS8。具体的には、着目共振周波数についての、S7の処理において得られた回路モデルのパラメータが一つずつ個別に微小変化したときの周波数の変化量を参照し、変化の程度が他と比べて大きいパラメータを選定する。変化の程度が他と比べて大きいと判断するための基準は、特定の基準に限定されるものではなく、全てのパラメータ毎の周波数の変化の程度の中において相対的に大きいか否かが判断される基準が適宜設定される。具体的には例えば、パラメータ毎の周波数の変化の程度の多くが01未満である場合において02~03を超える変化をしているパラメータについては、変化の程度が大きく、周波数変化に対して影響が大きいパラメータであると選定することなどが考えられる。本実施形態では、制御部11のパラメータ選定部11eが、S7の処理においてメモリ15に記憶された着目共振周波数についての回路モデルのパラメータ毎の周波数変化量をメモリ15から読み込む。そして、パラメータ選定部11eは、読み込んだ周波数変化量の中から、他と比べて変化の程度が大きいパラメータを選定する。なお、変化の程度が他と比べて大きいと判断するための基準は、変圧器の異常診断プログラム17内に予め規定しておくようにしても良いし、判断基準の指定を要求するメッセージをS8の処理の段階で表示部14に表示させると共に入力部13を介して入力されたものをパラメータ選定部11eに与えるようにしても良いし、全てのパラメータ毎の周波数の変化の程度に基づいて判断基準を設定する仕組み例えばサブルーチンを変圧器の異常診断プログラム17内に規定するようにしても良い。そして、パラメータ選定部11eは、選定したパラメータの種類をメモリ15に記憶させる。最後に、S8の処理の結果に基づいて異常診断対象の変圧器における異常様相の診断が行われるS9。本発明では、S8の処理において選定された着目共振周波数における周波数変化に対して微小変化の影響が大きいパラメータに関係する素子や素子の特定部位に異常が発生していると判断することにより、変圧器における異常様相異常発生箇所の診断を行う。具体的には例えば、S8の処理において1次巻線のインダクタンスL1が選定された場合には1次巻線に異常が発生していると判断する。なお、本発明においては、回路モデルのパラメータの設定の細かさに対応して異常様相異常発生箇所の特定も精緻になるので、回路モデルのパラメータの設定を細かくすることによって異常様相の診断の精度を向上させることが可能になる。例えば、上述の例では1次巻線のインダクタンスL1と1次巻線とを対応させているだけであるが、1次巻線に纏わる回路モデルパラメータの設定を細かくすることにより、1次巻線の上端部,中間部,下端部などのように異常が発生している部位の特定も可能になり得る。本発明による診断結果の具体的な内容としては、例えば、巻線の変形が生じている具体的な部位の特定や、巻線の位置ずれの様相左右のずれ,上下のずれ,スペーサが劣化によって収縮し全体が縮むなどして各ターンの隙間が小さくなるの特定や、鉄心の変形などが考えられる。ただし、本発明による診断結果の具体的な内容はこれらに限られるものではなく、S8の処理において選定されたパラメータの種類場合によっては複数のパラメータの組み合わせと実際に変圧器を確認することによって判明した異常様相との組み合わせの実績を蓄積することによって異常様相の内容が充実したり異常様相の特定の精度が向上したりし得る。S8の処理において選定されたパラメータの種類場合によっては複数のパラメータの組み合わせと異常様相異常発生箇所との対応は、変圧器の異常診断プログラム17内に予め規定しておくようにしても良いし、パラメータの種類組み合わせと異常様相異常発生箇所との組み合わせを定義したデータファイルを記憶部12やデータサーバ16等の記憶装置に予め保存しておくと共に当該データファイルから読み込むようにしても良い。本実施形態では、異常様相診断部11fがS8の処理において選定されてメモリ15に記憶されたパラメータの種類を読み込むと共に異常様相異常発生箇所を特定する。そして、異常様相診断部11fは、特定した異常様相異常発生箇所を、表示部14に表示させたり、診断結果データファイルとして記憶部12やデータサーバ16等の記憶装置に保存させたりする。そして、制御部11は異常診断対象の変圧器に対する変圧器の異常診断の処理を終了するEND。以上のように構成された本発明の変圧器の異常診断方法、異常診断装置及び異常診断プログラムによれば、診断対象である実際の変圧器の構成に沿う回路であって回路モデルのパラメータを変化させていない回路の共振周波数に対する回路モデルのパラメータを一つずつ個別に変化させたときの共振周波数の変化量を計算し、共振周波数の変化の程度が他と比べて大きいパラメータを選定することにより、異常が発生している素子や素子の部位を特定することができる。なお、上述の形態は本発明の好適な形態の一例ではあるがこれに限定されるものではなく、本発明の要旨を逸脱しない範囲において種々変形実施可能である。例えば、上述の実施形態では実際の変圧器の構成に沿う回路であって回路モデルのパラメータを変化させていない回路の共振周波数を、伝達関数を計算することなく変圧器回路の共振周波数のみを計算するために状態方程式の導出を行ってS5共振周波数の計算を行うS6ようにしているが、これに限られず、例えば汎用されている回路解析ソフトウェアを用いて伝達関数を計算して共振周波数を抽出するようにしても良い。また、上述の実施形態ではS7の処理において回路モデルの各パラメータを個別に微小変化させたときそれぞれの共振周波数の変化の程度を計算するために固有値の変化量∂λ/∂αを表す数式5を用いて固有値感度解析を行うようにしているが、これに限られず、例えば汎用されている回路解析ソフトウェアを用いて回路モデルのパラメータを一つ変化させた状態毎に共振周波数を計算するようにしても良いし、或いは、上述の実施形態のS5及びS6の処理と同様に状態方程式を導出して共振周波数を計算するようにしても良い。"