JP3793423B2 - 評価対象の状態評価装置、方法、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体 - Google Patents
評価対象の状態評価装置、方法、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体 Download PDFInfo
- Publication number
- JP3793423B2 JP3793423B2 JP2001136221A JP2001136221A JP3793423B2 JP 3793423 B2 JP3793423 B2 JP 3793423B2 JP 2001136221 A JP2001136221 A JP 2001136221A JP 2001136221 A JP2001136221 A JP 2001136221A JP 3793423 B2 JP3793423 B2 JP 3793423B2
- Authority
- JP
- Japan
- Prior art keywords
- state
- evaluation
- attractor
- point
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Images
Landscapes
- Testing Or Calibration Of Command Recording Devices (AREA)
- Complex Calculations (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、評価対象の状態評価装置、方法、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体に関する。
【０００２】
【従来の技術】
従来は、ある評価対象の状態を評価する場合、評価対象の状態が反映された物理量（温度等）を計測し、その計測結果に基づいて評価対象の状態を推定していた。
【０００３】
【発明が解決しようとする課題】
しかしながら、計測結果に基づいて評価対象の状態を推定していくだけでは、評価対象の時間的変化を捉えにくいといった問題があった。
【０００４】
さらに、計測点が評価対象から離れているような場合には、評価対象物と計測点との間での減衰等により計測する物理量が変化してしまうため、評価対象の状態を正確に評価することができないといった問題があった。
【０００５】
本発明は上記のような点に鑑みてなされたものであり、評価対象の時間的変化を捉えることができ、更には、評価対象と計測点とが離れている場合でも、計測点において計測した信号から評価対象の位置における信号を求め、評価対象の状態を正確に評価できるようにすることを目的とする。
【０００６】
【課題を解決するための手段】
本発明の評価対象の状態評価装置は、評価対象から離れた計測点で計測した信号に基づいて、評価対象の状態を評価する評価対象の状態評価装置であって、上記計測点で計測した信号に基づいて、逆問題解析により上記評価対象の位置における信号を求める逆問題解析手段と、上記逆問題解析手段により求めた評価対象の位置における信号から所定の次元を有する再構成アトラクタを作成するアトラクタ作成手段と、上記アトラクタ作成手段により作成された再構成アトラクタ上の基準時刻での点を中心とする直径εの超球に存在する再構成アトラクタ上の近傍点が、上記基準時間から所定時間推移後での点を中心とする直径εの超球にいくつ存在するかの割合を求め、その割合に基づいて状態の評価を行う状態評価手段とを備えた点に特徴を有する。
【０００８】
また、本発明の評価対象の状態評価装置の他の特徴とするところは、上記状態評価手段は、上記割合に応じて評価対象の状態の安定性又は定常状態からの乖離度を評価する点にある。
【００１０】
本発明の評価対象の状態評価方法は、評価対象から離れた計測点で計測した信号に基づいて、評価対象の状態を評価する評価対象の状態評価方法であって、上記計測点で計測した信号に基づいて、逆問題解析により上記評価対象の位置における信号を求める逆問題解析処理を有し、上記逆問題解析処理により求めた評価対象の位置における信号から所定の次元を有する再構成アトラクタを作成するアトラクタ作成処理を有し、上記アトラクタ作成処理により作成された再構成アトラクタ上の基準時刻での点を中心とする直径εの超球に存在する再構成アトラクタ上の近傍点が、上記基準時間から所定時間推移後での点を中心とする直径εの超球にいくつ存在するかの割合を求め、その割合に基づいて状態の評価を行い、状態評価に応じて診断処理及び制御処理のうち少なくとも1つの処理を行う点に特徴を有する。
【００１２】
本発明のコンピュータプログラムは、評価対象から離れた計測点で計測した信号に基づいて、評価対象の状態を評価する処理をコンピュータに実行させるコンピュータプログラムであって、上記計測点で計測した信号に基づいて、逆問題解析により上記評価対象の位置における信号を求める逆問題解析処理をコンピュータに実行させ、上記逆問題解析処理により求めた評価対象の位置における信号から所定の次元を有する再構成アトラクタを作成するアトラクタ作成処理と、上記アトラクタ作成処理により作成された再構成アトラクタ上の基準時刻での点を中心とする直径εの超球に存在する再構成アトラクタ上の近傍点が、上記基準時間から所定時間推移後での点を中心とする直径εの超球にいくつ存在するかの割合を求め、その割合に基づいて状態の評価を行う状態評価処理とをコンピュータに実行させる点に特徴を有する。
【００１４】
本発明のコンピュータ読み取り可能な記憶媒体は、上記コンピュータプログラムを格納した点に特徴を有する。
【００１５】
上記のようにした本発明においては、所定の次元を有する再構成アトラクタを作成し、再構成アトラクタ上の基準時刻での点の周囲に存在する再構成アトラクタ上の近傍点が、上記基準時間から所定時間推移後での点の周囲にいくつ存在するかの割合を求めることによって、評価対象の状態を時間的な変化を捉えて評価することができる。
【００１６】
【発明の実施の形態】
以下、図面を参照して、本発明の評価対象の状態評価装置、方法、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体の実施の形態について説明する。
【００１７】
図１には、本実施の形態の状態評価装置の構成を示す。同図において、１３０は計測部であり、評価対象１１０から所定の距離を有する計測点において、評価対象１１０から出力され、抵抗１２０により減衰された状態信号を計測する。上記評価対象１１０から出力される状態信号は、評価対象１１０の状態が反映された信号であり、熱、信号、音等がある。すなわち、計測部１３０は、計測点における評価対象１１０の状態が反映された物理量等を計測する。抵抗１２０は、評価対象１１０から出力される状態信号を減衰させるものである。例えば、上記状態信号が熱に関する信号であれば、大きな熱容量を有するものや熱伝導度の低いもの等の伝熱抵抗となるものである。
【００１８】
１４０は状態評価装置であり、上記計測部１３０における状態信号の計測結果（以下、「計測点における計測値」と称する）に基づいて、評価対象１１０の外表面又は外表面上の点（以下、「評価面」と称する）における状態信号を求めて、評価対象の状態を評価する。状態評価装置は、逆問題解析部１４１、アトラクタ作成部１４２、記憶部１４３、状態評価部１４４を含んで構成される。
【００１９】
逆問題解析部１４１は、計測点における計測値及び抵抗１２０の物性値等を用いて逆問題解析を行い、評価対象１１０の評価面における状態信号を求める。すなわち、逆問題解析部１４１は、仮に評価対象１１０の位置で計測したと仮定した場合に計測される状態信号の算出を行う。
【００２０】
上記逆問題解析では、例えば、評価対象１１０、抵抗１２０、及び計測部１３０を含む系を対象にした所定の方程式（偏微分方程式等）と、評価対象１１０の評価面における状態信号の仮定値とを用いて、計測点における状態信号を算出する。その算出した状態信号と、計測部１３０により実際に計測された計測点における計測値との誤差が所定の値より小さくなるように、評価対象１１０の評価面における状態信号の仮定値を修正し、計測点における状態信号の算出を繰り返す。そして、算出した状態信号と、実際に計測された計測点における計測値との誤差が所定の値より小さくなったときの仮定値を、評価対象１１０の評価面における状態信号とする。
【００２１】
また、例えば、下記の数１に示す式（１）、（２）に基づいて、評価対象１１０の評価面における状態信号を算出する。
【００２２】
【数１】
【００２３】
上記式（１）は非定常方程式であり、例えば、熱伝導であれば熱伝導方程式である。式（１）に対して所定の演算等を施すと、式（２）に示すような積分境界方程式になる。式（２）において、Ｇは共役方程式の解、ｕはスカラー量（例えば、熱伝導であれば温度等）、∂ｕ／∂ｎはスカラー勾配（例えば、熱伝導であれば熱流束）である。
【００２４】
上記式（２）において、左辺は評価対象１１０の評価面に関する積分であり、右辺は所定の既知境界面、例えば計測点を含む面に関する積分である。したがって、計測点における計測値に基づいて、式（２）の右辺の値が求められ、その求めた値から式（２）の左辺の評価面におけるスカラー勾配∂ｕ／∂ｎが求められる。さらに、上述のようにして得られた評価面におけるスカラー勾配∂ｕ／∂ｎを、計測点における計測値を境界条件として解くことにより、評価面のスカラー値ｕを算出する。このスカラー値ｕが、評価面における状態信号である。
【００２５】
アトラクタ作成部１４２は、逆問題解析部１４１により算出された時系列の評価面における状態信号に基づいて、アトラクタと呼ばれる軌道を再構成する。まず、アトラクタ作成部１４２は、逆問題解析部１４１により算出された時系列の評価面における状態信号から、対象とする現象の２倍以上の次元ｍを持つ遅延ベクトルｖ（ｔ）＝（ｕ（ｔ），ｕ（ｔ＋τ），ｕ（ｔ＋２τ），…，ｕ（ｔ＋（ｍ−１）τ））を作成する。
【００２６】
ここで、上記遅延ベクトルｖ（ｔ）において、ｕ（Ｔ）は時刻Ｔにおける評価面における状態信号の値、τは時間遅れ間隔である。例えば、上記逆問題解析部１４１が、上述した式（１）、（２）を用いて、評価面における状態信号を算出した場合には、上記ｕ（Ｔ）は評価面のスカラー値ｕである。なお、上記ｕ（Ｔ）として、評価面におけるスカラー勾配∂ｕ／∂ｎを用いてもよい。
【００２７】
次に、アトラクタ作成部１４２は、上記作成した遅延ベクトルｖ（ｔ）を所定の次元を有する位相空間に写像する。この写像した遅延ベクトルｖ（ｔ）の時間推移による軌道を作成することによりアトラクタを再構成する。なお、以下の説明では、この再構成したアトラクタを「再構成アトラクタ」と称する。
【００２８】
１４３は記憶部であり、逆問題解析部１４１、アトラクタ作成部１４２、状態評価部１４４での処理に用いるデータを記憶するものである。記憶部１４３には、計測部１３０により計測された状態信号の計測値、逆問題解析部１４１により算出された評価面における状態信号、アトラクタ作成部１４２により作成された再構成アトラクタ等が記憶される。
【００２９】
１４４は状態評価部であり、上記アトラクタ作成部１４２により得られた結果に基づいて、評価対象の状態を評価する。評価対象の状態の挙動をΔｔの時間スケールで観測したときに、時間発展の様子が決定論的、すなわちある法則性に支配されて推移するようにみえるということは、図２に示すように、再構成された軌道群の近接した部分がΔｔ後に同じように近傍した部分に移されることを意味する。そこで、状態評価部１４４では、再構成アトラクタ上の現時刻点ｘ（ｔ）を中心とする直径εの超球を考え、そこに含まれる再構成アトラクタ上の近傍点が、Δｔ秒後の点ｘ（ｔ＋Δｔ）を中心とする直径εの超球にいくつ存在するかの割合を指標（以下、「定常性評価指標」と称する）として定義する。
【００３０】
具体的には、定常性評価指標は、定常性評価指標＝（Δｔ時刻後に生き残った近傍点数）／（時刻ｔにおける近傍点数）により表される。この定常性評価指標は、時系列変化の法則性依存度を表し、法則性依存度が大きいほど１に近づき、完全なランダム状態、例えば白色ノイズでは０となる。図２においては、定常性評価指標は３／５＝０．６となる。
【００３１】
以下、図３のフローチャートに基づいて、評価対象の状態評価のための処理動作について説明する。まず、逆問題解析部１４１において、計測部１３０から供給される計測値を取り込み（ステップＳ３０１）、その計測値を用いて上述したように逆問題解析処理を行い、評価対象１１０の評価面における状態信号を求める（ステップＳ３０２）。
【００３２】
図４には、参考として、ある事例において、２種類の計測点における計測値（同図（Ａ）に示す４１、４２）と、これら２種類の計測点における計測値を用いて逆問題解析により求めた評価面における状態信号（同図（Ｂ）に示す４４、４５）とを示す。同図（Ａ）に示すように、計測点における計測値は、抵抗２０により減衰されるために滑らかな状態で検出されるのに対して、同図（Ｂ）に示すように、評価面における状態信号は、鋭い変化を示す。特に、同図（Ａ）の領域４３では、計測点における計測値変化は小さいが、同図（Ｂ）で上記領域４３に対応する領域４６では、状態信号は大きな変化を示している。
【００３３】
次に、アトラクタ作成部１４２では、上記ステップＳ３０２で算出された評価面における状態信号に基づいて、上述したようにアトラクタを再構成する（ステップＳ３０３）。すなわち、逆問題解析部１４１により算出された時系列の評価面における状態信号から、対象とする現象の２倍以上の次元ｍを持つ遅延ベクトルｖ（ｔ）を作成する。そして、その遅延ベクトルｖ（ｔ）を所定の次元を有する位相空間に写像し、遅延ベクトルｖ（ｔ）の時間推移による軌道を作成することによりアトラクタを再構成する。
【００３４】
図５には、参考として、図４で述べた事例において、図４（Ｂ）に示す評価面における状態信号に基づいて遅延ベクトルｖ（ｔ）を作成し、その遅延ベクトルｖ（ｔ）を（ｘ（ｔ），ｘ（ｔ＋７τ），ｘ（ｔ＋（ｍ−１）τ））を座標としてプロットすることにより得られた再構成アトラクタを示す。同図において、軌道が密な領域６１は安定な状態を示す軌道領域であり、軌道が疎な領域６２、６３は不安定な状態を示す軌道領域である。
【００３５】
そして、状態評価部１４４では、再構成アトラクタ上の現時刻点ｘ（ｔ）を中心とする直径εの超球を考え、上述したように、直径εの超球内に存在する再構成アトラクタ上の近傍点（過去の軌道上の点）が、Δｔ秒後の点ｘ（ｔ＋Δｔ）を中心とする直径εの超球内にいくつ存在しているかという定常性評価指標を求める。そして、定常性評価指標が１に近ければ、評価対象は安定状態にあり、０に近ければ、評価対象は不安定な状態にあると評価する（ステップＳ３０４）。
【００３７】
以下、具体例を挙げて説明する。例えば、高炉操業において炉底温度状態（評価対象の状態）を評価するために、炉底温度状態を直接計測することはできないため、高炉炉底煉瓦に熱電対を埋め込んでおき、その熱電対で測定された温度から炉底温度状態を評価する場合を説明する。
【００３８】
図６には、経過時間と、高炉炉底煉瓦に埋め込まれた熱電対（「測定部３０」に相当）により計測された温度との関係を示す。同図の矢印Ａに示すように、何らかの原因により熱電対温度の急激な上昇が発生している。
【００３９】
図７には、上記逆問題解析を行わずに、熱電対温度に基づいてアトラクタを再構成し（図示せず）、定常性評価指標を求めた結果を示す。同図の矢印Ｂに示すように、図６の矢印Ａに示した変化の前に、定常性評価指標の低下、すなわち不安定な状態が観察される。これは、現在の状態が別の状態に移行しようとする予兆を表していると考えられる。
【００４０】
図８には、図６、７で説明したのと同じ高炉における、経過時間と、逆問題解析により求めた炉底熱流束との関係を示す図である。このように逆問題解析を行うことにより、炉壁に埋め込まれた熱電対（「測定部３０」に相当）により測定された温度（「計測点における計測値」に相当）に基づいて、炉底熱流束（「評価面における状態信号」に相当）を求めることが可能となる。この関係からも、同図の矢印Ａに示すように、何らかの原因により熱電対温度の急激な上昇が発生していることがわかる。
【００４１】
図９には、上記のようにして求めた炉底熱流束に基づいてアトラクタを再構成し（図示せず）、定常性評価指標を求めた結果を示す。同図の矢印Ｂに示すように、図８の矢印Ａに示した変化の前に、定常性評価指標に大きな乱れと値の低下、すなわち不安定な状態が観察される。これは、図７と同様、現在の状態が別の状態に移行しようとする予兆を表していると考えられる。また、図７と比較して、逆問題解析を介在させることにより、この予兆が鮮明に現れており、予兆の検出能が向上していることがわかる。
【００４２】
以上述べたように本実施の形態によれば、評価対象から離れた計測点における計測値に基づいて逆問題解析を行うことにより、抵抗１２０の影響を取り除いて、評価面における状態信号を求めることができる。そして、求められた評価面における状態信号の振る舞いを示すアトラクタを再構成し、ある法則性に支配されて推移するか否かを表す定常性評価指標を求めることにより、評価対象の状態の時間的変化を捉えて正確に評価することができる。
【００４３】
なお、定常性評価指標を求めるための超球の直径εは、評価対象に応じて定めればよく、例えば、評価を厳しく行う場合には、超球の直径εを小さくすればよい。
【００４４】
また、本実施の形態では、評価対象の状態を評価するものとしているが、さらに評価対象の将来の状態を予測したり、評価対象の状態を制御したりするのに用いてもよい。例えば、定常性評価指標が１に近い値となるように、評価対象に対して制御を行うようにすればよい。
【００４５】
（その他の実施の形態）
なお、上記実施の形態の状態評価装置は、コンピュータのＣＰＵ或いはＭＰＵ、ＲＡＭ、ＲＯＭ等により構成されるものであり、ＲＡＭやＲＯＭに格納されたプログラムが動作することによって各機能が実現される。この場合、上記機能を果たすようにコンピュータを動作させるプログラム自体、さらには当該プログラムを供給するための手段、例えばかかるプログラムを格納した記録媒体は本発明を構成する。記録媒体としては、例えばフレキシブルディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、磁気テープ、不揮発性のメモリカード、ＲＯＭ等を用いることができる。
【００４６】
また、コンピュータが供給されたプログラムコードを実行することにより、上述の実施の形態の機能が実現されるだけでなく、そのプログラムコードがコンピュータにおいて稼働しているＯＳ（オペレーティングシステム）或いは他のアプリケーションソフト等と共同して上述の実施の形態の機能が実現される場合にもかかるプログラムコードは本発明の実施の形態に含まれる。
【００４７】
また、本発明をネットワーク環境で利用すべく、全部或いは一部のプログラムが他のコンピュータで実行されるようになっていてもよい。例えば、画面入力処理は、遠隔端末コンピュータで行われ、各種判断、ログ記録等は他のセンターコンピュータ等で行われるようにしてもよい。上記実施の形態において示した各部の形状及び構造は、何れも本発明を実施するにあたっての具体化のほんの一例を示したものに過ぎず、これらによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその精神、又はその主要な特徴から逸脱することなく、様々な形で実施することができる。
【００４８】
【発明の効果】
以上述べたように本発明によれば、評価面における状態信号の振る舞いを示すアトラクタを再構成し、ある法則性に支配されて推移するか否かを表す割合を求めるようにしたので、評価対象の状態の時間的変化を捉えて評価することができる。更に、計測点が評価対象から離れているような場合に、計測点における計測値に基づいて逆問題解析を行うようにすれば、抵抗の影響を取り除いて評価面における状態信号を求めることができ、評価対象の状態を正確に評価することが可能となる。
【図面の簡単な説明】
【図１】本実施の形態の状態評価装置の構成を示すブロック図である。
【図２】アトラクタ近傍点の時間推移を説明するための図である。
【図３】状態評価のための処理動作を示すフローチャートである。
【図４】ある事例における２種類の計測点における計測値と、これら２種類の計測点における計測値を用いて逆問題解析により求めた評価面における状態信号とを示す図である。
【図５】ある事例における再構成アトラクタを示す図である。
【図６】ある高炉における経過時間と、高炉炉底煉瓦に埋め込まれた熱電対により計測された温度との関係を示す図である。
【図７】熱電対温度に基づいてアトラクタを再構成し、定常性評価指標を求めた結果を示す図である。
【図８】ある高炉における経過時間と、逆問題解析により求めた炉底熱流束との関係を示す図である。
【図９】逆問題解析により求めた炉底熱流束に基づいてアトラクタを再構成し、定常性評価標を求めた結果を示す図である。
【符号の説明】
１４０ 状態評価装置
１４１ 逆問題解析部
１４２ アトラクタ作成部
１４３ 記憶部
１４４ 状態評価部
【発明の属する技術分野】
本発明は、評価対象の状態評価装置、方法、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体に関する。
【０００２】
【従来の技術】
従来は、ある評価対象の状態を評価する場合、評価対象の状態が反映された物理量（温度等）を計測し、その計測結果に基づいて評価対象の状態を推定していた。
【０００３】
【発明が解決しようとする課題】
しかしながら、計測結果に基づいて評価対象の状態を推定していくだけでは、評価対象の時間的変化を捉えにくいといった問題があった。
【０００４】
さらに、計測点が評価対象から離れているような場合には、評価対象物と計測点との間での減衰等により計測する物理量が変化してしまうため、評価対象の状態を正確に評価することができないといった問題があった。
【０００５】
本発明は上記のような点に鑑みてなされたものであり、評価対象の時間的変化を捉えることができ、更には、評価対象と計測点とが離れている場合でも、計測点において計測した信号から評価対象の位置における信号を求め、評価対象の状態を正確に評価できるようにすることを目的とする。
【０００６】
【課題を解決するための手段】
本発明の評価対象の状態評価装置は、評価対象から離れた計測点で計測した信号に基づいて、評価対象の状態を評価する評価対象の状態評価装置であって、上記計測点で計測した信号に基づいて、逆問題解析により上記評価対象の位置における信号を求める逆問題解析手段と、上記逆問題解析手段により求めた評価対象の位置における信号から所定の次元を有する再構成アトラクタを作成するアトラクタ作成手段と、上記アトラクタ作成手段により作成された再構成アトラクタ上の基準時刻での点を中心とする直径εの超球に存在する再構成アトラクタ上の近傍点が、上記基準時間から所定時間推移後での点を中心とする直径εの超球にいくつ存在するかの割合を求め、その割合に基づいて状態の評価を行う状態評価手段とを備えた点に特徴を有する。
【０００８】
また、本発明の評価対象の状態評価装置の他の特徴とするところは、上記状態評価手段は、上記割合に応じて評価対象の状態の安定性又は定常状態からの乖離度を評価する点にある。
【００１０】
本発明の評価対象の状態評価方法は、評価対象から離れた計測点で計測した信号に基づいて、評価対象の状態を評価する評価対象の状態評価方法であって、上記計測点で計測した信号に基づいて、逆問題解析により上記評価対象の位置における信号を求める逆問題解析処理を有し、上記逆問題解析処理により求めた評価対象の位置における信号から所定の次元を有する再構成アトラクタを作成するアトラクタ作成処理を有し、上記アトラクタ作成処理により作成された再構成アトラクタ上の基準時刻での点を中心とする直径εの超球に存在する再構成アトラクタ上の近傍点が、上記基準時間から所定時間推移後での点を中心とする直径εの超球にいくつ存在するかの割合を求め、その割合に基づいて状態の評価を行い、状態評価に応じて診断処理及び制御処理のうち少なくとも1つの処理を行う点に特徴を有する。
【００１２】
本発明のコンピュータプログラムは、評価対象から離れた計測点で計測した信号に基づいて、評価対象の状態を評価する処理をコンピュータに実行させるコンピュータプログラムであって、上記計測点で計測した信号に基づいて、逆問題解析により上記評価対象の位置における信号を求める逆問題解析処理をコンピュータに実行させ、上記逆問題解析処理により求めた評価対象の位置における信号から所定の次元を有する再構成アトラクタを作成するアトラクタ作成処理と、上記アトラクタ作成処理により作成された再構成アトラクタ上の基準時刻での点を中心とする直径εの超球に存在する再構成アトラクタ上の近傍点が、上記基準時間から所定時間推移後での点を中心とする直径εの超球にいくつ存在するかの割合を求め、その割合に基づいて状態の評価を行う状態評価処理とをコンピュータに実行させる点に特徴を有する。
【００１４】
本発明のコンピュータ読み取り可能な記憶媒体は、上記コンピュータプログラムを格納した点に特徴を有する。
【００１５】
上記のようにした本発明においては、所定の次元を有する再構成アトラクタを作成し、再構成アトラクタ上の基準時刻での点の周囲に存在する再構成アトラクタ上の近傍点が、上記基準時間から所定時間推移後での点の周囲にいくつ存在するかの割合を求めることによって、評価対象の状態を時間的な変化を捉えて評価することができる。
【００１６】
【発明の実施の形態】
以下、図面を参照して、本発明の評価対象の状態評価装置、方法、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体の実施の形態について説明する。
【００１７】
図１には、本実施の形態の状態評価装置の構成を示す。同図において、１３０は計測部であり、評価対象１１０から所定の距離を有する計測点において、評価対象１１０から出力され、抵抗１２０により減衰された状態信号を計測する。上記評価対象１１０から出力される状態信号は、評価対象１１０の状態が反映された信号であり、熱、信号、音等がある。すなわち、計測部１３０は、計測点における評価対象１１０の状態が反映された物理量等を計測する。抵抗１２０は、評価対象１１０から出力される状態信号を減衰させるものである。例えば、上記状態信号が熱に関する信号であれば、大きな熱容量を有するものや熱伝導度の低いもの等の伝熱抵抗となるものである。
【００１８】
１４０は状態評価装置であり、上記計測部１３０における状態信号の計測結果（以下、「計測点における計測値」と称する）に基づいて、評価対象１１０の外表面又は外表面上の点（以下、「評価面」と称する）における状態信号を求めて、評価対象の状態を評価する。状態評価装置は、逆問題解析部１４１、アトラクタ作成部１４２、記憶部１４３、状態評価部１４４を含んで構成される。
【００１９】
逆問題解析部１４１は、計測点における計測値及び抵抗１２０の物性値等を用いて逆問題解析を行い、評価対象１１０の評価面における状態信号を求める。すなわち、逆問題解析部１４１は、仮に評価対象１１０の位置で計測したと仮定した場合に計測される状態信号の算出を行う。
【００２０】
上記逆問題解析では、例えば、評価対象１１０、抵抗１２０、及び計測部１３０を含む系を対象にした所定の方程式（偏微分方程式等）と、評価対象１１０の評価面における状態信号の仮定値とを用いて、計測点における状態信号を算出する。その算出した状態信号と、計測部１３０により実際に計測された計測点における計測値との誤差が所定の値より小さくなるように、評価対象１１０の評価面における状態信号の仮定値を修正し、計測点における状態信号の算出を繰り返す。そして、算出した状態信号と、実際に計測された計測点における計測値との誤差が所定の値より小さくなったときの仮定値を、評価対象１１０の評価面における状態信号とする。
【００２１】
また、例えば、下記の数１に示す式（１）、（２）に基づいて、評価対象１１０の評価面における状態信号を算出する。
【００２２】
【数１】
上記式（１）は非定常方程式であり、例えば、熱伝導であれば熱伝導方程式である。式（１）に対して所定の演算等を施すと、式（２）に示すような積分境界方程式になる。式（２）において、Ｇは共役方程式の解、ｕはスカラー量（例えば、熱伝導であれば温度等）、∂ｕ／∂ｎはスカラー勾配（例えば、熱伝導であれば熱流束）である。
【００２４】
上記式（２）において、左辺は評価対象１１０の評価面に関する積分であり、右辺は所定の既知境界面、例えば計測点を含む面に関する積分である。したがって、計測点における計測値に基づいて、式（２）の右辺の値が求められ、その求めた値から式（２）の左辺の評価面におけるスカラー勾配∂ｕ／∂ｎが求められる。さらに、上述のようにして得られた評価面におけるスカラー勾配∂ｕ／∂ｎを、計測点における計測値を境界条件として解くことにより、評価面のスカラー値ｕを算出する。このスカラー値ｕが、評価面における状態信号である。
【００２５】
アトラクタ作成部１４２は、逆問題解析部１４１により算出された時系列の評価面における状態信号に基づいて、アトラクタと呼ばれる軌道を再構成する。まず、アトラクタ作成部１４２は、逆問題解析部１４１により算出された時系列の評価面における状態信号から、対象とする現象の２倍以上の次元ｍを持つ遅延ベクトルｖ（ｔ）＝（ｕ（ｔ），ｕ（ｔ＋τ），ｕ（ｔ＋２τ），…，ｕ（ｔ＋（ｍ−１）τ））を作成する。
【００２６】
ここで、上記遅延ベクトルｖ（ｔ）において、ｕ（Ｔ）は時刻Ｔにおける評価面における状態信号の値、τは時間遅れ間隔である。例えば、上記逆問題解析部１４１が、上述した式（１）、（２）を用いて、評価面における状態信号を算出した場合には、上記ｕ（Ｔ）は評価面のスカラー値ｕである。なお、上記ｕ（Ｔ）として、評価面におけるスカラー勾配∂ｕ／∂ｎを用いてもよい。
【００２７】
次に、アトラクタ作成部１４２は、上記作成した遅延ベクトルｖ（ｔ）を所定の次元を有する位相空間に写像する。この写像した遅延ベクトルｖ（ｔ）の時間推移による軌道を作成することによりアトラクタを再構成する。なお、以下の説明では、この再構成したアトラクタを「再構成アトラクタ」と称する。
【００２８】
１４３は記憶部であり、逆問題解析部１４１、アトラクタ作成部１４２、状態評価部１４４での処理に用いるデータを記憶するものである。記憶部１４３には、計測部１３０により計測された状態信号の計測値、逆問題解析部１４１により算出された評価面における状態信号、アトラクタ作成部１４２により作成された再構成アトラクタ等が記憶される。
【００２９】
１４４は状態評価部であり、上記アトラクタ作成部１４２により得られた結果に基づいて、評価対象の状態を評価する。評価対象の状態の挙動をΔｔの時間スケールで観測したときに、時間発展の様子が決定論的、すなわちある法則性に支配されて推移するようにみえるということは、図２に示すように、再構成された軌道群の近接した部分がΔｔ後に同じように近傍した部分に移されることを意味する。そこで、状態評価部１４４では、再構成アトラクタ上の現時刻点ｘ（ｔ）を中心とする直径εの超球を考え、そこに含まれる再構成アトラクタ上の近傍点が、Δｔ秒後の点ｘ（ｔ＋Δｔ）を中心とする直径εの超球にいくつ存在するかの割合を指標（以下、「定常性評価指標」と称する）として定義する。
【００３０】
具体的には、定常性評価指標は、定常性評価指標＝（Δｔ時刻後に生き残った近傍点数）／（時刻ｔにおける近傍点数）により表される。この定常性評価指標は、時系列変化の法則性依存度を表し、法則性依存度が大きいほど１に近づき、完全なランダム状態、例えば白色ノイズでは０となる。図２においては、定常性評価指標は３／５＝０．６となる。
【００３１】
以下、図３のフローチャートに基づいて、評価対象の状態評価のための処理動作について説明する。まず、逆問題解析部１４１において、計測部１３０から供給される計測値を取り込み（ステップＳ３０１）、その計測値を用いて上述したように逆問題解析処理を行い、評価対象１１０の評価面における状態信号を求める（ステップＳ３０２）。
【００３２】
図４には、参考として、ある事例において、２種類の計測点における計測値（同図（Ａ）に示す４１、４２）と、これら２種類の計測点における計測値を用いて逆問題解析により求めた評価面における状態信号（同図（Ｂ）に示す４４、４５）とを示す。同図（Ａ）に示すように、計測点における計測値は、抵抗２０により減衰されるために滑らかな状態で検出されるのに対して、同図（Ｂ）に示すように、評価面における状態信号は、鋭い変化を示す。特に、同図（Ａ）の領域４３では、計測点における計測値変化は小さいが、同図（Ｂ）で上記領域４３に対応する領域４６では、状態信号は大きな変化を示している。
【００３３】
次に、アトラクタ作成部１４２では、上記ステップＳ３０２で算出された評価面における状態信号に基づいて、上述したようにアトラクタを再構成する（ステップＳ３０３）。すなわち、逆問題解析部１４１により算出された時系列の評価面における状態信号から、対象とする現象の２倍以上の次元ｍを持つ遅延ベクトルｖ（ｔ）を作成する。そして、その遅延ベクトルｖ（ｔ）を所定の次元を有する位相空間に写像し、遅延ベクトルｖ（ｔ）の時間推移による軌道を作成することによりアトラクタを再構成する。
【００３４】
図５には、参考として、図４で述べた事例において、図４（Ｂ）に示す評価面における状態信号に基づいて遅延ベクトルｖ（ｔ）を作成し、その遅延ベクトルｖ（ｔ）を（ｘ（ｔ），ｘ（ｔ＋７τ），ｘ（ｔ＋（ｍ−１）τ））を座標としてプロットすることにより得られた再構成アトラクタを示す。同図において、軌道が密な領域６１は安定な状態を示す軌道領域であり、軌道が疎な領域６２、６３は不安定な状態を示す軌道領域である。
【００３５】
そして、状態評価部１４４では、再構成アトラクタ上の現時刻点ｘ（ｔ）を中心とする直径εの超球を考え、上述したように、直径εの超球内に存在する再構成アトラクタ上の近傍点（過去の軌道上の点）が、Δｔ秒後の点ｘ（ｔ＋Δｔ）を中心とする直径εの超球内にいくつ存在しているかという定常性評価指標を求める。そして、定常性評価指標が１に近ければ、評価対象は安定状態にあり、０に近ければ、評価対象は不安定な状態にあると評価する（ステップＳ３０４）。
【００３７】
以下、具体例を挙げて説明する。例えば、高炉操業において炉底温度状態（評価対象の状態）を評価するために、炉底温度状態を直接計測することはできないため、高炉炉底煉瓦に熱電対を埋め込んでおき、その熱電対で測定された温度から炉底温度状態を評価する場合を説明する。
【００３８】
図６には、経過時間と、高炉炉底煉瓦に埋め込まれた熱電対（「測定部３０」に相当）により計測された温度との関係を示す。同図の矢印Ａに示すように、何らかの原因により熱電対温度の急激な上昇が発生している。
【００３９】
図７には、上記逆問題解析を行わずに、熱電対温度に基づいてアトラクタを再構成し（図示せず）、定常性評価指標を求めた結果を示す。同図の矢印Ｂに示すように、図６の矢印Ａに示した変化の前に、定常性評価指標の低下、すなわち不安定な状態が観察される。これは、現在の状態が別の状態に移行しようとする予兆を表していると考えられる。
【００４０】
図８には、図６、７で説明したのと同じ高炉における、経過時間と、逆問題解析により求めた炉底熱流束との関係を示す図である。このように逆問題解析を行うことにより、炉壁に埋め込まれた熱電対（「測定部３０」に相当）により測定された温度（「計測点における計測値」に相当）に基づいて、炉底熱流束（「評価面における状態信号」に相当）を求めることが可能となる。この関係からも、同図の矢印Ａに示すように、何らかの原因により熱電対温度の急激な上昇が発生していることがわかる。
【００４１】
図９には、上記のようにして求めた炉底熱流束に基づいてアトラクタを再構成し（図示せず）、定常性評価指標を求めた結果を示す。同図の矢印Ｂに示すように、図８の矢印Ａに示した変化の前に、定常性評価指標に大きな乱れと値の低下、すなわち不安定な状態が観察される。これは、図７と同様、現在の状態が別の状態に移行しようとする予兆を表していると考えられる。また、図７と比較して、逆問題解析を介在させることにより、この予兆が鮮明に現れており、予兆の検出能が向上していることがわかる。
【００４２】
以上述べたように本実施の形態によれば、評価対象から離れた計測点における計測値に基づいて逆問題解析を行うことにより、抵抗１２０の影響を取り除いて、評価面における状態信号を求めることができる。そして、求められた評価面における状態信号の振る舞いを示すアトラクタを再構成し、ある法則性に支配されて推移するか否かを表す定常性評価指標を求めることにより、評価対象の状態の時間的変化を捉えて正確に評価することができる。
【００４３】
なお、定常性評価指標を求めるための超球の直径εは、評価対象に応じて定めればよく、例えば、評価を厳しく行う場合には、超球の直径εを小さくすればよい。
【００４４】
また、本実施の形態では、評価対象の状態を評価するものとしているが、さらに評価対象の将来の状態を予測したり、評価対象の状態を制御したりするのに用いてもよい。例えば、定常性評価指標が１に近い値となるように、評価対象に対して制御を行うようにすればよい。
【００４５】
（その他の実施の形態）
なお、上記実施の形態の状態評価装置は、コンピュータのＣＰＵ或いはＭＰＵ、ＲＡＭ、ＲＯＭ等により構成されるものであり、ＲＡＭやＲＯＭに格納されたプログラムが動作することによって各機能が実現される。この場合、上記機能を果たすようにコンピュータを動作させるプログラム自体、さらには当該プログラムを供給するための手段、例えばかかるプログラムを格納した記録媒体は本発明を構成する。記録媒体としては、例えばフレキシブルディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、磁気テープ、不揮発性のメモリカード、ＲＯＭ等を用いることができる。
【００４６】
また、コンピュータが供給されたプログラムコードを実行することにより、上述の実施の形態の機能が実現されるだけでなく、そのプログラムコードがコンピュータにおいて稼働しているＯＳ（オペレーティングシステム）或いは他のアプリケーションソフト等と共同して上述の実施の形態の機能が実現される場合にもかかるプログラムコードは本発明の実施の形態に含まれる。
【００４７】
また、本発明をネットワーク環境で利用すべく、全部或いは一部のプログラムが他のコンピュータで実行されるようになっていてもよい。例えば、画面入力処理は、遠隔端末コンピュータで行われ、各種判断、ログ記録等は他のセンターコンピュータ等で行われるようにしてもよい。上記実施の形態において示した各部の形状及び構造は、何れも本発明を実施するにあたっての具体化のほんの一例を示したものに過ぎず、これらによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその精神、又はその主要な特徴から逸脱することなく、様々な形で実施することができる。
【００４８】
【発明の効果】
以上述べたように本発明によれば、評価面における状態信号の振る舞いを示すアトラクタを再構成し、ある法則性に支配されて推移するか否かを表す割合を求めるようにしたので、評価対象の状態の時間的変化を捉えて評価することができる。更に、計測点が評価対象から離れているような場合に、計測点における計測値に基づいて逆問題解析を行うようにすれば、抵抗の影響を取り除いて評価面における状態信号を求めることができ、評価対象の状態を正確に評価することが可能となる。
【図面の簡単な説明】
【図１】本実施の形態の状態評価装置の構成を示すブロック図である。
【図２】アトラクタ近傍点の時間推移を説明するための図である。
【図３】状態評価のための処理動作を示すフローチャートである。
【図４】ある事例における２種類の計測点における計測値と、これら２種類の計測点における計測値を用いて逆問題解析により求めた評価面における状態信号とを示す図である。
【図５】ある事例における再構成アトラクタを示す図である。
【図６】ある高炉における経過時間と、高炉炉底煉瓦に埋め込まれた熱電対により計測された温度との関係を示す図である。
【図７】熱電対温度に基づいてアトラクタを再構成し、定常性評価指標を求めた結果を示す図である。
【図８】ある高炉における経過時間と、逆問題解析により求めた炉底熱流束との関係を示す図である。
【図９】逆問題解析により求めた炉底熱流束に基づいてアトラクタを再構成し、定常性評価標を求めた結果を示す図である。
【符号の説明】
１４０ 状態評価装置
１４１ 逆問題解析部
１４２ アトラクタ作成部
１４３ 記憶部
１４４ 状態評価部
Claims (5)
- 評価対象から離れた計測点で計測した信号に基づいて、評価対象の状態を評価する評価対象の状態評価装置であって、
上記計測点で計測した信号に基づいて、逆問題解析により上記評価対象の位置における信号を求める逆問題解析手段と、
上記逆問題解析手段により求めた評価対象の位置における信号から所定の次元を有する再構成アトラクタを作成するアトラクタ作成手段と、
上記アトラクタ作成手段により作成された再構成アトラクタ上の基準時刻での点を中心とする直径εの超球に存在する再構成アトラクタ上の近傍点が、上記基準時間から所定時間推移後での点を中心とする直径εの超球にいくつ存在するかの割合を求め、その割合に基づいて状態の評価を行う状態評価手段とを備えたことを特徴とする評価対象の状態評価装置。 - 上記状態評価手段は、上記割合に応じて評価対象の状態の安定性又は定常状態からの乖離度を評価することを特徴とする請求項１に記載の評価対象の状態評価装置。
- 評価対象から離れた計測点で計測した信号に基づいて、評価対象の状態を評価する評価対象の状態評価方法であって、
上記計測点で計測した信号に基づいて、逆問題解析により上記評価対象の位置における信号を求める逆問題解析処理と、
上記逆問題解析処理により求めた評価対象の位置における信号から所定の次元を有する再構成アトラクタを作成するアトラクタ作成処理とを有し、
上記アトラクタ作成処理により作成された再構成アトラクタ上の基準時刻での点を中心とする直径εの超球に存在する再構成アトラクタ上の近傍点が、上記基準時間から所定時間推移後での点を中心とする直径εの超球にいくつ存在するかの割合を求め、その割合に基づいて状態の評価を行い、状態評価に応じて診断処理及び制御処理のうち少なくとも１つの処理を行うことを特徴とする評価対象の状態評価方法。 - 評価対象から離れた計測点で計測した信号に基づいて、評価対象の状態を評価する処理をコンピュータに実行させるコンピュータプログラムであって、
上記計測点で計測した信号に基づいて、逆問題解析により上記評価対象の位置における信号を求める逆問題解析処理と、
上記逆問題解析処理により求めた評価対象の位置における信号から所定の次元を有する再構成アトラクタを作成するアトラクタ作成処理と、
上記アトラクタ作成処理により作成された再構成アトラクタ上の基準時刻での点を中心とする直径εの超球に存在する再構成アトラクタ上の近傍点が、上記基準時間から所定時間推移後での点を中心とする直径εの超球にいくつ存在するかの割合を求め、その割合に基づいて状態の評価を行う状態評価処理とをコンピュータに実行させることを特徴とするコンピュータプログラム。 - 請求項４に記載のコンピュータプログラムを格納したことを特徴とするコンピュータ読み取り可能な記憶媒体。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001136221A JP3793423B2 (ja) | 2001-05-07 | 2001-05-07 | 評価対象の状態評価装置、方法、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001136221A JP3793423B2 (ja) | 2001-05-07 | 2001-05-07 | 評価対象の状態評価装置、方法、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2002333346A JP2002333346A (ja) | 2002-11-22 |
| JP3793423B2 true JP3793423B2 (ja) | 2006-07-05 |
Family
ID=18983544
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001136221A Expired - Lifetime JP3793423B2 (ja) | 2001-05-07 | 2001-05-07 | 評価対象の状態評価装置、方法、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3793423B2 (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002143997A (ja) * | 2000-11-10 | 2002-05-21 | Nippon Steel Corp | 鋳型内鋳片の状態検知装置、方法、及びコンピュータ読み取り可能な記憶媒体 |
| JP2002346715A (ja) * | 2001-05-18 | 2002-12-04 | Nippon Steel Corp | 鋳型内鋳片の状態評価装置、方法、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体 |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6616960B2 (ja) * | 2015-04-20 | 2019-12-04 | フクダ電子株式会社 | 生体信号処理装置およびその制御方法 |
-
2001
- 2001-05-07 JP JP2001136221A patent/JP3793423B2/ja not_active Expired - Lifetime
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002143997A (ja) * | 2000-11-10 | 2002-05-21 | Nippon Steel Corp | 鋳型内鋳片の状態検知装置、方法、及びコンピュータ読み取り可能な記憶媒体 |
| JP2002346715A (ja) * | 2001-05-18 | 2002-12-04 | Nippon Steel Corp | 鋳型内鋳片の状態評価装置、方法、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2002333346A (ja) | 2002-11-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107943874B (zh) | 知识图谱处理方法、装置、计算机设备及存储介质 | |
| JP3728198B2 (ja) | 評価対象の状態検出装置、状態検出方法およびコンピュータ読み取り可能な記録媒体 | |
| JP3793423B2 (ja) | 評価対象の状態評価装置、方法、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体 | |
| US20190178814A1 (en) | State assessing device, state assessing method, and storage medium for storing program | |
| WO2014033874A1 (ja) | 測定位置指示装置、測定位置指示方法 | |
| JP6981526B2 (ja) | システム同定装置、システム同定方法及びコンピュータプログラム | |
| JP7322332B2 (ja) | 測定機器 | |
| JP2021012078A (ja) | 温度分布推定装置及び温度分布推定方法 | |
| JP5845794B2 (ja) | 移動軌跡補間装置、移動軌跡補間方法、及びプログラム | |
| JP7036209B2 (ja) | 診断装置、診断方法、及びプログラム | |
| JP4074443B2 (ja) | 鋳型内鋳片の状態評価装置、方法、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体 | |
| JPWO2018198267A1 (ja) | 因果関係学習方法、プログラム、装置および異常分析システム | |
| CN109642827B (zh) | 温度推算系统、温度推算方法及记录有温度推算程序的记录媒体 | |
| JPWO2017077609A1 (ja) | 構造解析方法、及び構造解析プログラム | |
| JP7254262B2 (ja) | 作業推定装置、作業推定方法、及び、作業推定プログラム | |
| JP4041279B2 (ja) | 鋳型内鋳片の状態検知装置、方法、及びコンピュータ読み取り可能な記憶媒体 | |
| JP7297240B2 (ja) | ユーザ状態推定装置 | |
| JP7105447B2 (ja) | 形状推定装置、形状推定方法及びプログラム | |
| JP2002328914A (ja) | 状態類似度検出装置、状態類似度検出方法、記録媒体及びコンピュータプログラム | |
| JP4050871B2 (ja) | 高炉炉内状態の判断装置、判断方法および記録媒体 | |
| JPS60249069A (ja) | ステップ応答測定装置 | |
| JP6022623B2 (ja) | 熱貫流率推定システム、方法、およびプログラム、ならびに、熱貫流率試験装置 | |
| JPWO2018131311A1 (ja) | センシングシステム、センサノード装置、センサ測定値処理方法及びプログラム | |
| WO2022264342A1 (ja) | 深さ推定装置、深さ推定方法、及び深さ推定プログラム | |
| Adhikari et al. | SpiroMilli: Bringing Ad-hoc Spirometry to 5G Devices |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050214 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050222 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050425 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20060328 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20060407 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 3793423 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090414 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100414 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110414 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120414 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130414 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130414 Year of fee payment: 7 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130414 Year of fee payment: 7 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130414 Year of fee payment: 7 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130414 Year of fee payment: 7 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140414 Year of fee payment: 8 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |