JP5832644B2 - 殊にガスタービンまたは風力タービンのような技術システムのデータドリブンモデルを計算機支援で形成する方法 - Google Patents

殊にガスタービンまたは風力タービンのような技術システムのデータドリブンモデルを計算機支援で形成する方法 Download PDF

Info

Publication number
JP5832644B2
JP5832644B2 JP2014513206A JP2014513206A JP5832644B2 JP 5832644 B2 JP5832644 B2 JP 5832644B2 JP 2014513206 A JP2014513206 A JP 2014513206A JP 2014513206 A JP2014513206 A JP 2014513206A JP 5832644 B2 JP5832644 B2 JP 5832644B2
Authority
JP
Japan
Prior art keywords
data
model
technical system
method characterized
loop control
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.)
Active
Application number
JP2014513206A
Other languages
English (en)
Other versions
JP2014525063A (ja
Inventor
デュル ズィークムント
デュル ズィークムント
ヘンチェル アレクサンダー
ヘンチェル アレクサンダー
シュテアツィング フォルクマー
シュテアツィング フォルクマー
ウートルフト シュテフェン
ウートルフト シュテフェン
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Siemens AG
Original Assignee
Siemens AG
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Siemens AG filed Critical Siemens AG
Publication of JP2014525063A publication Critical patent/JP2014525063A/ja
Application granted granted Critical
Publication of JP5832644B2 publication Critical patent/JP5832644B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00Machine learning
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03DWIND MOTORS
    • F03D17/00Monitoring or testing of wind motors, e.g. diagnostics
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03DWIND MOTORS
    • F03D7/00Controlling wind motors 
    • F03D7/02Controlling wind motors  the wind motors having rotation axis substantially parallel to the air flow entering the rotor
    • F03D7/04Automatic control; Regulation
    • F03D7/042Automatic control; Regulation by means of an electrical or electronic controller
    • F03D7/043Automatic control; Regulation by means of an electrical or electronic controller characterised by the type of control logic
    • F03D7/046Automatic control; Regulation by means of an electrical or electronic controller characterised by the type of control logic with learning or adaptive control, e.g. self-tuning, fuzzy logic or neural network
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B17/00Systems involving the use of models or simulators of said systems
    • G05B17/02Systems involving the use of models or simulators of said systems electric
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B23/00Testing or monitoring of control systems or parts thereof
    • G05B23/02Electric testing or monitoring
    • G05B23/0205Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
    • G05B23/0218Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
    • G05B23/0224Process history based detection method, e.g. whereby history implies the availability of large amounts of data
    • G05B23/024Quantitative history assessment, e.g. mathematical relationships between available data; Functions therefor; Principal component analysis [PCA]; Partial least square [PLS]; Statistical classifiers, e.g. Bayesian networks, linear regression or correlation analysis; Neural networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/048Activation functions
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/044Recurrent networks, e.g. Hopfield networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/084Backpropagation, e.g. using gradient descent
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N7/00Computing arrangements based on specific mathematical models
    • G06N7/01Probabilistic graphical models, e.g. probabilistic networks
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/72Wind turbines with rotation axis in wind direction

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Software Systems (AREA)
  • Mathematical Physics (AREA)
  • General Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Automation & Control Theory (AREA)
  • Chemical & Material Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Computing Systems (AREA)
  • Mechanical Engineering (AREA)
  • Combustion & Propulsion (AREA)
  • Sustainable Development (AREA)
  • Sustainable Energy (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Computational Linguistics (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Fuzzy Systems (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Feedback Control In General (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Testing And Monitoring For Control Systems (AREA)
  • Wind Motors (AREA)

Description

本発明は、殊にガスタービンまたは風力タービンのような技術システムのデータドリブンモデルを計算機支援で形成する方法ならびに相応するコンピュータプログラム製品に関する。
技術システムを自動で動作させるためには、この技術システムの動作パラメタを表す相応するトレーニングデータに基づいて学習することが可能なデータドリブンモデルを使用することができる。このように学習されたデータドリブンモデルは、上記の技術システムの挙動を反映しようするものであり、引き続き、相応する技術システムの監視モデルおよび/または制御モデルを形成するために使用することができる。つぎにこれらの監視モデルないしは制御モデルは、実際の動作パラメタに基づいてこのシステムの相応する操作量を調整するかまたはこの技術システムのエラー動作状態を自動的に確定するため、この技術システムの実際の動作に使用することができる。
一般的にはデータドリブンモデルを学習する上記のトレーニングデータは、技術システムの動作パラメタの状態空間内に均一に分散されてはいない。一般的には通常動作を表す技術システムの動作領域には、この技術システムの動作の境界領域におけるデータよりも、ないしはこの技術システムのエラー動作状態に対するデータよりも格段に多くのデータが存在する。この問題は殊に、損傷を回避するためにふつう境界領域では動作されないガスタービンおよび風力タービンを制御する際に発生する。したがってガスタービンないしは風力タービンの頻繁に利用される動作領域においては、極めて大量の有効なトレーニングデータセットが存在するが、その情報内容は大部分が冗長的である。これに対してまれに運転が行われる動作モードに対してはわずかなデータしか存在しないのである。
したがってデータドリブンモデルを学習するために使用されるトレーニングデータに対し、データ密度の相異なる複数の領域が存在するのである。一般的にはこのデータ密度は、相応するトレーニングデータを選択する際、およびこれに基づいて、相応する技術システムのモデルを形成する際には考慮されない。これにより、上記のデータドリブンモデルに基づいて計算される相応する制御モデルおよび監視モデルは、技術システムの動作時に誤ったないしは最適以下の制御動作または誤った警告を出力する結果にもなり得るのである。したがって例えば、ガスタービンないしは風力タービンの動作時に発生しかつ上記のデータドリブンモデルのトレーニングフェーズに含まれていなかった寒空イベントが、誤った警報をトリガすることにもなり得るのである。また典型的な動作状態に対するトレーニングデータだけで学習されるデータドリブンモデルは、多くのモデルエラーを生じさせ得るため、これに基づく閉ループ制御モデルないしは開ループ制御モデルは、技術システムの実際動作時に、例えばガスタービンおいて供給されるガスの量ないしは案内羽根の調整などの、相応する操作量を最適に調整しないのである。このことは技術システムの出力ないしは効率ないしは寿命の低下に結び付き得るのである。
従来、技術システムのデータドリブンモデルを形成する際には、異常値またはエラーを取り除いた後、トレーニングデータのすべてのデータセットを使用するか、またはこれらのデータセットをランダムに選択している。すべてのトレーニングデータを考慮することにより、極めて大量かつ常時増大し続けるトレーニングデータセットが形成され、このトレーニングデータセットは、相応するデータドリブンモデルの形成を遅くし、データドリブンモデルを形成する1つまたは複数の計算機の資源問題(例えば、メインメモリサイズ)に結び付き得るのである。トレーニングデータセットのランダムな選択に単純に制限することは、データドリブンモデルの計算を確かに高速化するが、技術システムについての関連する情報が失われるという危険性が高まる。この問題は、技術システムの動的特性を表すデータドリブンモデルにおいて、モデルを学習する際に、時間的に関連するデータセットからなる所定の1つの集合が、つねに共通に考慮されることによって増大するのである。
本発明の課題は、一層量少ない計算資源によって技術システムの良好なデータドリブンモデルを学習する、計算機支援によるデータドリブンモデルの形成する方法を提供することである。
この課題は、独立請求項によって解決される。本発明の発展形態は従属請求項に定められている。
本発明による方法は、技術システムのデータドリブンモデルおよび殊にガスタービンまたは風力タービンのデータドリブンモデルを計算機支援で形成するために使用される。このデータドリブンモデルは有利には、技術システムにおける自動監視作業および/または予測作業および/または制御作業(MPC)に使用され、複数のトレーニングデータに基づいて形成される。これらのトレーニングデータには多くのデータセットが含まれており、これらのデータセットはそれぞれ、技術システムの複数の動作パラメタを表す。複数の動作パラメタとは、技術システムの測定量のことであり、また場合によっては操作量のことでもある。ここでこれらの測定量は、技術システムにおいてまたはその周囲環境において適切なセンサによって求められたものかまたはセンサ値から導出された量のことである。操作量は、技術システムの動作に影響を及ぼすために調整可能ないしは変更可能である、技術システムの動作パラメタのことである。
本発明による方法は繰り返して実行され、ここでは所定の回数の繰り返しステップに対してそれぞれ、トレーニングデータからならデータセットによって上記のデータドリブンモデルおよび1つ密度推定器が学習される。この場合に上記の繰り返しの終了後に学習されたデータドリブンモデルが、本発明による方法によって形成されるデータドリブンモデルである。上記の(学習された)密度推定器は、トレーニングデータからなる複数のデータセットに対して1つずつの信頼度尺度を出力し、この信頼度尺度は、トレーニングデータからなる各データセットと、別の複数のデータセットとの類似度が大きければ大きいほど高い。上記のデータセットの類似度は、有利な実施形態において、これらのデータセットによって表されるデータ点の空間における相応の間隔尺度によって表される。上記の密度推定器は、その信頼度尺度を介して、相応するデータセットに対し、上記のデータ空間において密度が高いか否か、ひいては類似した上記のトレーニングデータによるデータセットが多くあるか否かを出力する。上で今まさに説明した密度推定器を実現可能な方法は、従来技術から十分に公知であり、したがって詳しく説明しない。有利な1つの実施形態では、密度推定器として、ニューラルクラウドに基づく方法を使用する。有利には刊行物[1]に記載された方法を使用する。この刊行物の開示内容部全体は、参照によって本発明の内容にするものとする。
本発明による方法において繰り返しステップによって学習したデータドリブンモデルは、上記のトレーニングデータのデータセットを再構成するために使用される。ここでは、この再構成されるデータセットと、上記のデータセットによる実際のデータセットとの間の違いを表すモデルエラーが発生する。質的に良好なデータドリブンモデルにおいては、このモデルエラーを可能な限りに小さくしたい。ここでこのデータドリブンモデルは、殊に、トレーニングデータの各データセットの動作パラメタの第1の部分が、このデータドリブンモデルの入力値を表すのに対して、このデータドリブンモデルの出力値が、各データセットに含まれている動作パラメタの第2の部分になるように構成される。この場合に上記のモデルエラーは、上記のデータドリブンモデルを介して出力される動作パラメタの第2の部分と、トレーニングデータの相応するデータセットの動作パラメタの相応する第2の部分との差分によって得られる。
本発明による方法の枠内では、各繰り返しステップにおいた学習した上記の密度推定器およびデータドリブンモデルを用いて、つぎの繰り返しステップにおける学習のために上記のトレーニングデータからなるデータセットを選択するかまたは重み付けを行う。ここでは、トレーニングデータからなる信頼度尺度が低いデータセットおよびモデルエラーの大きいデータセットは、より前に選択されるかまたは大きく重み付けされる。すなわち、モデルエラーが同じ場合、別のデータセットよりも低い信頼度尺度を有する、トレーニングデータからなるデータセットは、上記の別のデータセットよりも前に選択されるかないしは大きく重み付けされるのである。また信頼度尺度が同じ場合、別のデータセットよりも大きなモデルエラーを有するデータセットは、上記の別のデータセットよりも有利に選択されるかないしは大きく重み付けされるのである。
本発明は、高い度合いで「新規性」を有するトレーニングデータ、すなわちその信頼度尺度が小さくかつその結果として類似のデータセットが少ないトレーニングデータを使用する場合に技術システムのデータドリブンモデルが極めて効率的に学習されるという知識に基づいている。同様に学習法に対して、モデルがまだエラーを有する、すなわちまだ大きなモデルエラーが発生するトレーニングデータの領域が関心の対象になる。したがって学習の新たな繰り返しステップにおいて信頼度尺度の低いないしはモデルエラーの大きなデータセットがより大きく考慮される適切な関係を上記のようにモデル化することにより、技術システムの良好なデータドリブンモデルを効率的に計算して形成できる学習が達成されるのである。
本発明による方法によれば、大きなトレーニングデータ集合内に存在する関係を適切に抽出することができ、またこの関係を技術システムの、殊にガスタービンまたは風力タービンの特性をモデル化するために利用することができる。相応する技術システムの良好なモデル化は、トレーニングデータのデータセットがわずかな動作パラメタしか含まないか、ないしはトレーニングデータの集合が少ないデータセットにわたって分散される場合であっても可能になる。本発明による方法では、トレーニングデータの重要な情報に関連するデータセットないしは領域だけを利用するため、相応するこのモデルは、高速かつ少ない計算資源で形成することができる。さらに本発明による方法により、形成されるデータドリブンモデルのパーセンタイルパフォーマンスが最適化される。このことは、上記のデータドリブンモデルが技術システムの計算機支援による制御または監視の枠内で使用される場合に殊に有利である。
本発明による方法はさらに、引き続いて技術システムの開ループ制御モデルおよび/または閉ループ制御モデルないしは監視モデルを計算する際に組み込み可能な学習型の密度推定器が、本発明によって提供されるという利点を有する。相応する開ループ制御モデルないしは監視モデルに使用されるこの密度推定器は、技術システムの実際の動作パラメタの重要性についての情報を供給する。殊に操作量の変化に対し、相応する報告の必要な閾値または最大値を動的に適合させることである。これにより、技術システムおよび殊にガスタービンまたは風力タービンをリアルタイムで監視する際のエラーアラームを回避することができる。さらに、形成されるデータドリブンモデルのオーソライズを、上記の信頼度尺度に依存させ、ひいては技術システムの実際の動作パラメタの相応するデータセットのデータ密度に依存させることができる。すなわち、データ密度が十分である場合にのみ、上記のデータドリブンモデルないしは技術システムの動作を受け入れるないしは継続するのである。これにより、例えば、ガスタービンの燃焼動的特性が大きすぎることによって技術システムの損傷に結び付き得る、技術システムの操作量設定を回避することができる。
本発明による方法によって形成されるデータドリブンモデルを用いれば、殊に技術システムの動的(すなわち時間に依存する)特性もモデル化することができる。この場合、データセットには、同じ時点の動作パラメタだけでなく、データドリブンモデルのアーキテクチャに相応し、関連している時点の一連の動作パラメタが含まれる。
本発明の方法によれば、データドリブンモデルとして有利にはニューラルネットワークが、また殊にリカレントニューラルネットワークが学習される。このようなリカレントニューラルネットワークを用いれば、相応するトレーニングデータに基づいて技術システムの特性を極めて良好に再構成することができる。上のニューラルネットワークの学習は、有利にはそれ自体公知のエラー・バックプロパブゲーション法によって行われる。殊に有利な実施形態では、ニューラルネットワークとして、刊行物DE 10 2010 011 221 A1に記載された、マルコフ決定過程抽出ネットワークとも称されるネットワークを使用する。
本発明による方法は、ガスタービンのモデル化に殊に良好に適している。データセットにおいて技術システムの動作パラメタとして考慮される、ガスタービンの適切な測定値には、以下の1つまたは複数のパラメタが含まれる。すなわち、
− ガスタービンが動作する周囲環境における温度、
− ガスタービンが動作する周囲環境における大気湿度、
− ガスタービンが動作する周囲環境における気圧、
− 殊にパーセンタイルパフォーマンス値のような前記ガスタービンの性能、
− 前記ガスタービンの燃料良度、
− 例えば窒素酸化物および/または一酸化炭素の放出のようなガスタービンの有害物質放出、
− ガスタービンの1つまたは複数のタービン羽根の温度、
− ガスタービンの排ガス温度、
− ガスタービンの燃焼動的特性
が含まれるのである。
ここで上で挙げた燃焼動的特性という概念は、適切なセンサで検出した、ガスタービンにおける圧力変動の振幅値に関するものである。
これに対し、ガスタービンについては、上記のデータセットにおける動作パラメタとして以下の1つまたは複数の操作量、すなわち、
− ガスタービンにおける1つまたは複数の案内羽根の位置、
− ガスタービンに供給されるガスの量、
− ガスタービンの排ガス温度に対する目標値
が殊に有利であることが判明している。
ここで注意すべきであるのは、ガスタービンでは一般的に、適切な制御器によってタービンの排ガス温度に対する所望値を調整できることである。この値は、上述の目標値に相応する。
本発明の方法の別の有利な使用分野は、風力タービン用のデータドリブンモデルの形成である。ここでは、上記のデータセットにおける動作パラメタの測定量として、有利には以下の1つまたは複数の量、すなわち、
− 風力タービンによって形成される電力、
− 風力タービンが動作される周囲環境における風速、
− 風力タービンの回転数、
− 風力タービンの回転翼における曲げ力および/または撓みおよび/または応力変動、
− 風力タービンの塔の振動
が使用される。
風力タービンの上記の動作パラメタの操作量とし、上記のデータセットにおいて有利には、以下の1つまたは複数のパラメタ、すなわち、
− 風力タービンの前記回転翼の迎え角、
− 回転数と、風力タービンによって形成される電力との間の関係を表す、風力タービンの特性マップ
を使用する。
殊に有利な別の実施形態では、本発明による方法によって形成したデータドリブンモデルを用い、1つまたは複数の最適化判定基準に基づいて、技術システムの動作に対する開ループ制御モデルおよび/または閉ループ制御モデルを計算する。計算したこの開ループ制御モデルおよび/または閉ループ制御モデルにより、技術システムの動作パラメタの1つまたは複数の測定量に依存して、この技術システムの動作パラメタの1つまたは複数の操作量の変化が得られ、これによって技術システムの実際動作において、上記の操作量を適切に設定することができる。
ガスタービンの上記の開ループ制御モデルおよび/または閉ループ制御モデルが計算されると、有利にも上記の1つまたは複数の最適化判定基準に、可能な限りに少ない有害物質放出および/または可能な限りに少なくとも燃焼動的特性が含まれ、または可能な限りに高い効率(熱消費率)または可能な限りに高い出力も含まれる。これにより、ガスタービンの環境的に健全な動作ないしはタービンの長い寿命が保証されるのである。
風力タービンに対する上記の開ループ制御モデルおよび/または閉ループ制御モデルが計算されると、上記の1つまたは複数の最適化判定基準には殊に、風力タービンの回転翼における可能な限りに小さい曲げ力および/または可能な限りに高い風力タービンの効率が含まれる。これにより、風力タービンの最適な動作が保証され、ないしはここでも長い寿命が保証されるのである。
殊に有利な別の実施形態では、上記の開ループ制御モデルおよび/または閉ループ制御モデルを構成して、上記の学習した密度推定器により、技術システムの動作時に発生する動作パラメタからなる複数のデータセットの各信頼度尺度が求められているように上記の開ループ制御および/または閉ループ制御モデルを構成する。ここでは、相応する信頼度尺度が、あらかじめ設定した閾値を下回る場合、上記の動作パラメタの1つまたは複数の操作量の上記の変更を制限し、および/または、上記のデータドリブンモデルに関連する機能の動作および/または技術システムの動作を中断するかまたは最初からスタートさせない。これにより、上記の技術システムの動作パラメタが、十分に高いデータ密度を有する領域内にある場合にのみ、この技術システムの動作をオーソライズするために、本発明の方法の枠内で学習した上記の密度推定器を使用することも可能である。
本発明による方法の別の実施形態では、形成したデータドリブンモデルにより、監視モデルを計算する。ここでこの監視モデルは、技術システムの動作時に発生する1つまたは複数の動作パラメタが、形成したデータドリブンモデルを介して再構成される1つまたは複数の相応する動作パラメタから、あらかじめ設定した量よりも偏差する場合、警告通知および/またはアラーム通知を出力し、および/または、この技術システムの自動化されたセーフティ手段をトリガする。有利にはこのような監視モデルにおいて上記の学習した密度推定器を適当に考慮することもできる。この場合には、技術システムの動作時に発生する動作パラメタからなるデータセットに対する各信頼度尺度が、上記の学習された密度推定器によって求められるようにこの監視モデルを構成し、ここで上記のあらかじめ設定した量は、上記のモデルの入力量に対する信頼度尺度が小さければ小さいほど、大きくなるように設定され、このあらかじめ設定した量に基づいて警告通知が形成される。これにより、データ密度の低いデータセットに対する、報告の必要な閾値が上がる。場合によっては、各信頼度尺度が、あらかじめ設定した閾値よりも小さい場合、データ密度が小さすぎるために目下のところは警告通知を出力することができず、ひいては監視できないことの通知を出力することも可能である。
殊に有利な実施形態では上記の開ループ制御モデルおよび/または閉ループ制御モデルおよび/または監視モデルは、強化学習に基づき、および/または、モデル予測制御に基づき、および/または最適化方法に基づいて計算される。強化学習ならびにモデル予測制御の方法は、従来技術から十分に公知である。相応する最適化方法も公知である。例えば、この最適化方法は、上記の開ループ制御モデルないしは閉ループ制御モデルないしは監視モデルにしたがい、相応に設定される最適値を求めるグラジエントベースまたは確率的なサーチ方法に基づいて行うことができる。殊に有利な別の実施形態では、例えばリカレントニューラルネットワークのようなニューラルネットワークに基づいて上記の開ループ制御モデルおよび/または閉ループ制御モデルおよび/または監視モデルを学習する。
相応の開ループ制御モデルおよび/または閉ループ制御モデルないしは監視モデルを計算するための技術としては、例えば、刊行物DE 10 2007 001 025 A1(Recurrent-Control-Neural-Network リカレントコントロールニューラルネットワーク)ないしはDE 10 2007 017 259 A1(Neural-Rewoards-Regression ニューラル報酬回帰)ないしはDE 10 2007 042 440 B3(Policy-Gradient-Neural-Reward-Regression ポリシィ勾配ニューラル報酬回帰)に記載されている方法が適している。相応する開ループ制御モデルおよび/または閉ループ制御モデルないしは監視モデルを計算するため、従来技術から公知のNeural-Fitted-Q-Iteration(刊行物[2]を参照されたい)も使用可能である。
場合によっては開ループ制御モデルおよび/または閉ループ制御モデルないしは監視モデルをニューラルネットワークによって学習する際にはレギュラリゼーション法を使用することができる。例えば、各データセットに対するエラーに加えて、上記のニューラルネットワークにおけるいわゆる重み付けに対するペナルティ項(英語ではweight penalty term)を追加し、これによってデータ密度の低い領域において一層線形であり、ひいては一層良好な外挿性を有するモデルを形成することができる。またそれ自体公知のInverse Kurtosis評価関数によって選択された、上記のニューラルネットワークの複数の重み付けをカットオフし、これによってデータ密度の小さいデータ領域において上記の方法を安定化させる。このことは、殊にガスタービンまたは風力タービンに対して開ループ制御モデルおよび/または閉ループ制御モデルを計算する際に有利である。
有利な実施形態では、上記のデータドリブンモデルを形成する際、上記のデータセットの重み付けないしは選択と、信頼度尺度ないしはモデルエラーとの間の関係を関数によってモデル化する。ここでは、信頼度尺度およびモデルエラーに依存する関数を定める。ここでこの関数は、データセットの上記の信頼度尺度およびモデルエラーに基づき、1つのデータセットを選択するか否か、ないしはこのデータセットの上記の重み付けをどの程度にするかを確定する。上記の関数は、適用事例に応じてさまざまに構成することが可能である。このような関数の具体的な決定は、当業者の行い得る範囲内にある。例えば、上記の関数は、上記の相応する関係を、階段によって不連続にモデル化するかまたは適当に連続にモデル化することができる。
本発明による方法の別の有利な実施形態では、上記のデータドリブンモデルを形成する際、上記のデータセットを選択して、このデータセットが、つぎの繰り返しステップに対し、上記の複数のデータセットに対応する確率に基づき、標本抽出によって取り出されるようにする。すなわち、上記のデータセットには、各信頼度尺度ないしはモデルエラーに基づいて適当に確率が割り当てられ、つぎにこれらがサンプリングの枠内で考慮されるのである。上記の確率の割り当ては、例えば、上で説明した関数によって行うことができる。
本発明による方法の別の実施形態では、上記のトレーニングデータからなるデータセットに、上記の密度推定器およびデータドリブンモデルを学習する際に考慮される重み付けを対応付ける。ここで各繰り返しステップにおいて学習した密度推定器またはデータドリブンモデルにより、低い信頼度尺度および大きなモデルエラーを有する、上記のトレーニングデータからなるデータセットがより大きく重み付けされるように、上記の重み付けを更新する。これらの重み付けは、種々異なる仕方で上記の密度推定器ないしはデータドリブンモデルよって考慮することが可能である。本発明の1つの変化形態において、上記の重み付けは、上記の密度推定器およびデータドリブンモデルによって処理される確率を表す。例えば、上記の密度推定器および/またはデータドリブンモデルにより、それらの確率に基づき、学習のためにトレーニングデータからなるデータセットを選択することができる。これらの確率は、学習の枠内において、例えば、上記の密度推定器またはデータドリブンモデルによって使用される確率的勾配降下法において適切に考慮することも可能である。
本発明による方法は、種々異なる仕方で初期化することができる。1つの実施形態では、例えばランダムに決定される、上記のトレーニングデータからなるデータセットの部分集合によって上記の方法を初期化する。また適切に決定した初期の重み付けを有する、トレーニングデータからなるすべてのデータセットによって上記の方法を初期化することも可能である。
本発明による方法の枠内で求めた上記のモデルエラーは、例えば、上記のデータドリブンモデルによって再構成したデータセット、ないしはこのデータセットの再構成された動作パラメタと、上記のトレーニングデータの相応するデータセット、ないしはこのデータセットの相応する動作パラメタとの偏差の絶対値によって表すことができる。また上記のモデルエラーは、相対的な偏差を表すこともでき、この相対的な偏差では、上記のデータセットのすべてのモデルエラーの標準偏差を基準にして上記の偏差の絶対値が求められ、これが各繰り返しステップにおいて再構成されるのである。
本発明による方法において出力される信頼度尺度も、例えばゼロと1との間の区間に適当に正規化され、ここで値ゼロは信頼度がない(すなわちデータ密度が低い)ことおよび値1は信頼度が高い(すなわちデータ密度)ことを表す。
本発明による方法の枠内で実行される繰り返しは、相異なる判定基準に基づいて終了させることが可能である。有利な実施形態では、つぎの繰り返しステップに対して選択されるかまた重み付けされたデータセットのモデルエラーの少なくとも1つのあらかじめ設定した割合が、あらかじめ設定した閾値を下回るか、および/または、つぎの繰り返しステップに対して選択されるかまたは重み付けされたデータセットの信頼度尺度の少なくとも1つのあらかじめ設定した割合が、あらかじめ設定した閾値を上回る場合に、上記の繰り返しを終了する。これによって保証されるのは、上記の学習したデータドリブンモデルが所定の最小品質に達した場合に上記の方法を終了することである。ここで上記のあらかじめ設定した割合は、例えば、すべてのモデルエラーないしは信頼度尺度の95%または99%以上とすることが可能である。
本発明による方法の別の実施形態では、つぎの繰り返しステップに対して選択および重み付けされたデータセットの平均モデルエラーが、先行する繰り返しステップに対してこれ以上は小さくならない、および/またはつぎの繰り返しステップに対して選択および重み付けされたデータセットのすべての信頼度尺度が、あらかじめ設定した閾値を上回る場合に、上記の繰り返しステップを終了する。
データドリブンモデルを形成する上記の方法に加え、本発明はさらに技術システムを動作させる方法に関しており、ここでは上記の本発明による方法を用いて、ないしは、本発明のよる方法の1つまたは複数の有利な変化形態を用いてデータドリブンモデルを形成するか、ないしは形成して引き続き、この技術システムの動作を解析するため、および/またはこれに影響(殊に直接の影響)を及ぼすために使用する。
したがってこの方法は、本発明によるこの方法は、技術システムの実際動作の枠内で、本発明によって形成した上記のデータドリブンモデルを使用するのである。ここでは有利には、開ループ制御モデルおよび/または閉ループ制御モデルないしは監視モデルを計算する変形形態を使用して上記のデータドリブンモデルを形成する。この場合に技術システムは、この開ループ制御モデルおよび/または閉ループ制御モデルおよび/または監視モデルにより、開ループ制御および/または閉ループ制御および/または監視される。
上で説明した方法に加え、本発明はさらに機械読み出し可能媒体に記憶されたプログラムコードを有するコンピュータプログラム製品に関しており、このコンピュータプログラムコードは、これが計算機で実行される場合に本発明による上記の方法、ないしは本発明による上記の方法の1つまたは複数の有利な変形形態を実行するためのものである。
以下では添付の図面に基づき、本発明の実施例を詳しく説明する。
本発明による方法の第1実施例の概略図である。 本発明による方法の第2実施例の概略図である。
以下で説明する本発明による方法の実施形態は、複数のデータセットを含むトレーニングデータを使用して、技術システムのデータドリブンモデルを形成する方法に使用される。これらの実施形態によって形成されるデータドリブンモデルは有利には、上で説明したように、技術システムに対する相応の開ループ制御モデルおよび/または閉ループ制御モデルないしは監視モデルを計算するために使用可能である。続いてこれらの開ループ制御モデルおよび/または閉ループ制御モデルないしは監視モデルは、同様に上で説明したように、技術システムの実際動作において使用することができる。以下で説明する実施形態の殊に有利な適用事例は、ガスタービンないしは風力タービン用のデータドリブンモデルの形成である。
ここで説明する実施形態において個々のデータセットは、技術システムの複数の動作パラメタを含むベクトルである。これらのような動作パラメタの例は、ガスタービンおよび風力タービンに対して上で挙げた。1つの実施形態において上記のベクトルの動作パラメタは、技術システムの時間的に連続する測定量および操作量を表すため、データドリブンモデルにより、技術システムの動的特性が学習される。上記のトレーニングデータにおける動作パラメタは、技術システムの実際動作において検出され、データベースに記憶されている。
本発明による方法の目的は、関連する情報を上記のトレーニングデータから抽出し、また冗長な情報を取り除くことであるため、上記のトレーニングデータによって学習されるデータドリブンモデルの品質が改善される。ここでは殊に、学習されるモデルのいわゆるパーセンタイルパフォーマンスを改善する。すなわち、技術システムの、データセットの少ない動作領域においてデータドリブンモデルの品質を高めようとしているのである。
図1には、本発明による方法の1つの実施形態が示されており、ここでは、相応するデータベースに格納されているトレーニングデータTDが処理される。すでに述べたように、これらのトレーニングデータは、技術システムの動作パラメタの形態のデータセットに関連しており、場合によっては事前にフィルタリングが行われているため、データドリブンモデルに関連するデータセットの所定の動作パラメタだけが、ないしは上記の方法における所定のデータセットだけが引き続いて処理される。図1の実施形態は、データドリブンモデルを学習するために、トレーニングデータTDの部分集合をサンプリング(すなわち標本抽出)によって適当に選択するというアイデアに基づいている。この方法を初期化する際にはまず、ステップRS(RS=Random Sampling)においてトレーニングデータの部分集合をランダムに選択する。この部分集合は、図1においてRD(RD=Reduced Dataset)と記されている。この部分集合により、ニューラルネットワークないしはニューラルモデルNMの形態のデータドリブンモデルを学習する。ここではこのために、例えば強化学習のようなそれ自体公知の学習方法を使用する。
これに並行して上記の部分集合により、密度推定器ないしはデータエンカプスレータDEを学習する。ここで説明している実施形態では、データエンカプスレータとして、刊行物[1]から公知のニューラルクラウドに基づく密度推定器を使用する。この方式に基づく密度推定器は公知であり、したがってここでは詳しく説明しない。この密度推定器の目標は、相応するデータセットないしはデータ点に対し、上記のトレーニングデータの別のデータセットが、観察中のデータセットにどの程度類似しているかという点についてその「新規性」を求めることである。ここでこの類似度は、データセット相互の距離によって表される。この距離は、データ空間におけるデータ点のユークリッド距離によって与えられる。刊行物[1]の密度推定器によりまず、min-max正規化法に基づき、個々のデータセットないしはデータ点の正規化を行う。引き続き、「拡張K平均」アルゴリズムとも称される従来技術から公知のK平均化アルゴリズムの変化形態に基づいて上記のデータ点のクラスリングを行う。このアルゴリズムにより、データ点の空間内に複数のクラスタ中心が得られる。つぎのステップでは、これらの中心にわたってガウスベルを配置し、引き続いてここれらのガウスベルの正規化を行う。つぎのこのアルゴリズムの最終形成物として、学習されたデータエンカプスレータが得られる。このデータエンカプスレータにより、各データセットに対して信頼度が求められ、この信頼度は、このデータセットと、トレーニングデータからなるデータセットとの類似度を表し、また上記の刊行物[1]では観察対象の技術システムに対する故障確率として評価される。
上記のニューラルネットワークNMの学習が終了した後、モデルエラーを除いて上記のトレーニングデータからなる複数のデータセットを再構成することができる。ここでは、データセットの動作パラメタの所定の割合が、上記のニューラルネットワークの入力ベクトルとして処理される。このニューラルネットワークはつぎに、相応するデータセットの別の動作パラメタを出力ベクトルとして供給する。上記のトレーニングデータによる実際の出力ベクトルと、上記のニューラルネットワークを介して計算した出力ベクトルとの間の差分により、相応するモデルエラーが得られる。
図1の実施形態では、データエンカプスレータDEの相応する信頼度CMと、上記のニューラルネットワークNMの相応するモデルエラーMEとを使用して、重み付きサンプリングWSを実行する。この重み付きサンプリングを行うため、トレーニングデータTDからなるすべてのデータセットに対してそれぞれ、上記の学習したニューラルネットワークによって得られるモデルエラーと、上記のデータエンカプスレータによって得られる信頼度とを求める。これらの値の組み合わせ毎に、相応するデータセットに対して確率を求める。これらの確率は、すべてのデータセットに対して1に正規化される。これらの確率は有利には、モデルエラーが同じ場合には、信頼度の低いデータセットが、信頼度の高いデータセットよりも高い確率を有するように決定される。さらに上記の確率は、信頼度が同じ場合にモデルエラーの大きなデータセットが、モデルエラーの小さいデータセットよりも大きな確率を有するように決定される。引き続いてこれらの確率によって上記のサンプリングWSが行われる。ここでは、このサンプリングにより、信頼度が低いことによって新規性が高いかないしはこのデータドリブンモデルによっては良好には再構成できないデータセットが上記のトレーニングデータからあらかじめ取り出される。
標本抽出WSを実行した後、データセットの新たな部分集合RDが得られ、この部分集合は、つぎの繰り返しステップにおいてデータエンカプスレータDEおよびニューラルネットワークNMをさらに学習するために使用可能である。ここでも引き続き、新たに重み付けしたサンプリングWSが行われる。しかしながらこのために、このサンプリングに対する相応の確率は新たに、再度学習したデータエンカプスレータないしは再度学習したニューラルネットワークによって決定される。最終的に上記の方法は、新たに形成したデータセットの部分集合により、つぎの繰り返しステップにおいて継続されるのである。
データエンカプスレータないしはニューラルネットワークの上記の繰り返しの学習は、適当な停止判定基準に基づいて終了される。この停止断判定基準は、新たに形成した部分集合のすべてのデータセットに対して各モデルエラーがあらかじめ設定した閾値よりも小さく、かつ、各信頼度があらかじめ定めた引き位置よりも大きい場合に上記の方法が終了するように決定される。この方法の最終結果として、引き続いて技術システムの開ループ制御モデルおよび/または閉ループモデルないしは監視モデルを計算するのに使用可能な学習されたニューラルネットワークが得られるのである。
図2には、本発明による方法の第2の実施形態が示されている。図1とは異なり、上記のデータエンカプスレータDEおよびニューラルネットワークNMを学習する際には(場合によって事前にフィルタリングした)トレーニングデータTDのすべてのデータセットを使用する。しかしながらこれらのデータセットは、上記の学習したデータエンカプスレータないしは学習したニューラルネットワークを使用して適切に重み付けされる。図2に示した方法は、はじめに適当に初期化され、ここではこの初期化の枠内で、観察対象のN個のデータセットに対する相応する重み付けw1,w2…,wNが上記のトレーニングデータにしたがって求められる。これらの重み付けは、例えば、ランダムにまたは値1/Nで初期化することができる。ここで説明している実施形態においてこれらの重み付けは、データエンカプスレータDEまたはニューラルネットNMによって考慮される確率を表す。
図1の実施形態とは異なり、ここではすべてのトレーニングデータが事前の選択なしに考慮される。上記の重み付けによって表される確率により、刊行物[1]の密度推定器に相応するデータエンカプスレータは、相応するデータセットを選択する。上記の確率は、データエンカプスレータにおいて上記の正規化の際に組み込むことができ、ここではこれらの確率にしたがい、上記の正規化の際に考慮されるデータセットが選択される。ここではまた上で説明したクラスタリングの際に上記の確率を使用することも可能である。ここでは上記の確率により、クラスタ中心をサーチする際のどのデータセットを使用するのかを選択することができる。またここでは、上記のデータエンカプスレータにおいて使用されるクラスタ中心の距離計算を行う際に、これらの確率を使用することができる。ここでも、上記の距離を計算する際にどのデータセットを考慮するのかが、上記の確率によって決定される。上記の正規化、クラスタリングおよび距離の計算の枠内では、上記の相応に使用されるデータセットは、一度だけ求めるかないかはステップ毎に新たに別個に求めることができる。
さらに上記の重み付きにより、上記のニューラルネットワークの学習のために考慮されるかないしニューラルネットワークの学習の枠内で考慮される相応のデータが選択される。例えば上記の学習の際には上記のデータセットの選択に基づき、公知の確率的勾配降下法の枠内で上記のグラジエントを求めることができる。
上記のデータエンカプスレータないしはニューラルネットワークを学習した後、更新ステップUPが行われ、ここでは上記の学習したエンカプスレータおよび学習したニューラルネットワークに基づき、最初に初期化した重み付けが新たに計算される。これは、図1に類似してつぎのようにすることによって行われる。すなわち、上記のトレーニングデータからなるすべてのデータセットに対し、上記のデータエンカプスレータの相応する信頼度および上記のニューラルネットワークの相応するモデルエラーを求め、これにより、考察対象のデータセットの信頼度が小さければ小さいほど、ないしは考察対象のデータセットのモデルエラーが大きければ大きいほど確率が高くなるように、重み付けないしは確率を導出されることによって行われるのである。ここでは、例えば、1つの関数に基づく関係のような相応の関係を決定することは当業者がなし得ることである。例えば、相応する関数によって線形関係をモデル化することができる。
図2の方法は、ステップUPにおいて更新された重み付けにより、つぎの繰り返しステップにおいて継続される。この繰り返しステップでは、上記のデータエンカプスレータおよびニューラルネットワークはこれらの新たな重み付けを考慮して再度学習され、引き続き、学習したデータエンカプスレータないしは学習したニューラルネットワークにより、上記の重み付けが新たに更新される。この繰り返しの方法は、適切な中断判定基準によって終了し、この中断判定基準は、図1に類似して、各データセットに対して発生するデータエンカプスレータの信頼度が、あらかじめ設定した閾値を上回るか、ないしは各データセットに対して発生するニューラルネットワークのモデルエラーが、あらかじめ設定した閾値を下回る場合、上記の方法が終了するように決定される。このようにした学習したニューラルネットワークはここでも、引き続き、相応する技術システムを制御するため、ないしは監視するためのモデルの計算に使用することができる。
本発明の上で説明した実施形態は、一連の利点を有する。殊にトレーニングデータの大きな集合から、データセットのベースにある技術システムを良好に表すことができるデータセットを抽出することができる。ここでは上記のデータセットから部分集合を適当に選択することにより、ないしは上記のデータセットを重み付けることにより、少ないデータセットによって強力な領域が考慮されるデータドリブンモデルを学習することができる。したがって、通常動作からは逸脱する、技術システムの動作状態が一層強力に考慮されるのである。
本発明による方法によれば、関連するデータセットを上記のように選択することにより、トレーニングデータに基づき、相応するデータドリブンモデルを高速に学習することができる。さらに学習したこのデータドリブンモデルは、良好なパーセンタイルパフォーマンスを有する。なぜならば有利にも、データの少ない複数の領域において学習されるからである。良好なパーセンタイルパフォーマンスは、技術システムの動作において、エラーを正しく予測するために重要である。
参考文献
[1] B. Lang et al., Neural Clouds for Monitoring of Complex Systems, Optical Memory and Neural Networks (Information Optics), 2008, Vol. 17, No.3, p.183-192
[2] M. Riedmiller, Neural Fitted Q Iteration - First Experiences with a Data Efficient Neural Reinforcement Learning Method, Proc. of the European Conf. on Machine Learning, 2005, p.317 - 328

Claims (28)

  1. 自動監視作業および/または診断作業および/または制御作業に使用するために、ガスタービンまたは風力タービンの技術システムのデータドリブンモデル(NM)を計算機支援で形成する方法であって、
    前記技術システムの複数の動作パラメタをそれぞれ表す複数のデータセットを含むトレーニングデータに基づいて前記データドリブンモデルを形成する、方法において、
    − 所定数の繰り返しステップに対し、前記トレーニングデータ(TD)からなるデータセットを用いて、前記データドリブンモデル(NM)および密度推定器(DE)をそれぞれ学習させ、ただし、
    前記密度推定器(DE)は、前記トレーニングデータからなる複数のデータセットに対して1つずつの信頼度尺度(CM)を出力し、ここで当該信頼度尺度(CM)は、各データセットと、前記トレーニングデータ(TD)からなる別のデータセットとの類似度が高ければ高いほど大きく、
    前記データドリブンモデル(NM)は、それぞれ1つのモデルエラー(ME)を有する前記各トレーニングデータのデータセットを再構成し、
    − 各繰り返しステップにおいて学習した前記データドリブンモデル(NM)および前記密度推定器(DE)を用いて、つぎの繰り返しステップにおける学習のため、前記トレーニングデータ(TD)からなるデータセットを選択するかまたは重み付けし、その際に低い信頼度尺度(CM)、および、大きなモデルエラー(ME)を有する、前記トレーニングデータ(TD)からなるデータセットがより前に選択されるかまたは大きく重み付けされるようにする、
    ことを特徴とする方法。
  2. 請求項1に記載の方法において、
    データドリブンモデル(NM)として、ニューラルネットワークを学習する、
    ことを特徴とする方法。
  3. 請求項1または2に記載の方法において、
    前記技術システムの前記動作パラメタには、ガスタービンの以下の1つまたは複数の測定量、すなわち、
    − 前記ガスタービンが動作する周囲環境における温度、
    − 前記ガスタービンが動作する周囲環境における大気湿度、
    − 前記ガスタービンが動作する周囲環境における気圧、
    前記ガスタービンの性能、
    − 前記ガスタービンの燃料良度、
    前記ガスタービンの有害物質放出、
    − 前記ガスタービンの1つまたは複数のタービン羽根の温度、
    − 前記ガスタービンの排ガス温度、
    − 前記ガスタービンの燃焼動的特性
    が含まれる、
    ことを特徴とする方法。
  4. 請求項1から3までのいずれか1項に記載の方法において、
    前記技術システムの前記動作パラメタには、ガスタービンの以下の1つまたは複数の操作量、すなわち、
    − 前記ガスタービンにおける1つまたは複数の案内羽根の位置、
    − 前記ガスタービンに供給されるガスの量、
    − 前記ガスタービンの排ガス温度に対する目標値
    が含まれる、
    ことを特徴とする方法。
  5. 請求項1から4までのいずれか1項に記載の方法において、
    前記技術システムの前記動作パラメタには、風力タービンの以下の1つまたは複数の測定量、すなわち、
    − 前記風力タービンによって形成される電力、
    − 前記風力タービンが動作する周囲環境における風速、
    − 前記風力タービンの回転数、
    − 前記風力タービンの回転翼における曲げ力および/または撓みおよび/または応力変動、
    − 前記風力タービンの塔の振動
    が含まれる、
    ことを特徴とする方法。
  6. 請求項1から5までのいずれか1項に記載の方法において、
    前記技術システムの前記動作パラメタには、風力タービンの以下の1つまたは複数の操作量、すなわち、
    − 前記風力タービンの回転翼の迎え角、
    − 回転数と、前記風力タービンによって形成される電力との間の関係を表す、前記風力タービンの特性マップ
    が含まれる、
    ことを特徴とする方法。
  7. 請求項1から6までのいずれか1項に記載の方法において、
    形成された前記データドリブンモデル(NM)を用い、1つまたは複数の最適化判定基準に基づいて、前記技術システムの動作に対する開ループ制御モデルおよび/または閉ループ制御モデルを計算し、
    前記開ループ制御モデルおよび/または閉ループ制御モデルにより、前記技術システムの前記動作パラメタの1つまたは複数の測定量に依存して、前記技術システムの前記動作パラメタの1つまたは複数の操作量の変化を得る、
    ことを特徴とする方法。
  8. 請求項7に記載の方法において、
    ガスタービンの前記開ループ制御モデルおよび/または閉ループ制御モデルを計算し、
    前記1つまたは複数の最適化判定基準には、有害物質放出および/または燃焼動的特性および/または効率および/または出力が含まれる、
    ことを特徴とする方法。
  9. 請求項7または8に記載の方法において、
    風力タービンの前記開ループ制御モデルおよび/または閉ループ制御モデルを形成し、
    前記1つまたは複数の最適化判定基準には、前記風力タービンの効率および/または前記風力タービンの回転翼における曲げ力および/または撓みおよび/または応力変動が含まれる、
    ことを特徴とする方法。
  10. 請求項7から9までのいずれか1項に記載の方法において、
    前記学習した密度推定器(DE)により、前記技術システムの前記動作時に発生する動作パラメタからなる複数のデータセットの前記各信頼度尺度(CM)が求められるように前記開ループ制御モデルおよび/または閉ループ制御モデルを構成し、
    相応する前記信頼度尺度(CM)が、あらかじめ設定した閾値を下回る場合、前記動作パラメタの1つまたは複数の操作量の前記変化を制限し、および/または、前記データドリブンモデルに関連する機能の動作および/または前記技術システムの動作を中断するかまたはスタートさせない、
    ことを特徴とする方法。
  11. 請求項1から10までのいずれか1項に記載の方法において、
    形成した前記データドリブンモデルによって監視モデルを計算し、
    該監視モデルは、前記技術システムの動作時に発生する1つまたは複数の動作パラメタが、形成した前記データドリブンモデル(NM)を介して再構成される1つまたは複数の相応する動作パラメタから、あらかじめ設定した量よりも偏差する場合、警告通知および/またはアラーム通知を出力し、および/または、前記技術システムの自動化されたセーフティ手段をトリガする、
    ことを特徴とする方法。
  12. 請求項11に記載の方法において、
    前記技術システムの動作時に発生する動作パラメタからなるデータセットに対する各信頼度尺度(CM)が、前記学習された密度推定器によって求められるように前記監視モデルを構成し、
    前記あらかじめ設定した量は、前記各信頼度尺度(CM)が小さければ小さいほど、大きくなるように設定され、前記あらかじめ設定した量に基づいて警告通知が形成される、
    ことを特徴とする方法。
  13. 請求項7から10までのいずれか1項に記載の方法において、
    強化学習に基づき、および/または、モデル予想制御に基づき、および/または最適化方法に基づき、前記開ループ制御モデルおよび/または閉ループ制御モデルを計算する、
    ことを特徴とする方法。
  14. 請求項7から10までのいずれか1項に記載の方法において、
    ニューラルネットワークに基づいて前記開ループ制御モデルおよび/または閉ループ制御モデルを学習する、
    ことを特徴とする方法。
  15. 請求項1から14までのいずれか1項に記載の方法において、
    前記信頼度尺度(CM)および前記モデルエラー(ME)に依存する関数を介して、前記トレーニングデータ(TD)からなる前記データセットの前記選択または重み付けを行う、
    ことを特徴とする方法。
  16. 請求項1から15までのいずれか1項に記載の方法において、
    つぎの繰り返しステップに対し、前記トレーニングデータからなる前記データセットが、前記データセットに対応づけられる確率に基づき、標本抽出によって取り出されるように前記トレーニングデータからなるデータセットを選択する、
    ことを特徴とする方法。
  17. 請求項1から16までのいずれか1項に記載の方法において、
    前記トレーニングデータ(TD)からなる前記データセットに、前記密度推定器(DE)および前記データドリブンモデル(NM)を学習する際に考慮される重み付け(w1,w2,…,wN)を対応付け、
    各繰り返しステップにおいて学習した前記密度推定器(DE)および前記データドリブンモデル(NM)により、低い信頼度尺度(CM)および高いモデルエラー(ME)を有する、前記トレーニングデータ(TD)からなるデータセットが、より大きく重み付けされるように、前記重み付け(w1,w2,…,wN)を更新する、
    ことを特徴とする方法。
  18. 請求項17に記載の方法において、
    前記重み付け(w)は、前記密度推定器(DE)および前記データドリブンモデル(NM)によって処理される確率を表す、
    ことを特徴とする方法。
  19. 請求項18に記載の方法において、
    前記密度推定器(DE)および/または前記データドリブンモデル(NM)により、前記確率に基づき、前記学習のためまたは前記学習の枠内で、前記トレーニングデータ(TD)からなるデータセットを選択する、
    ことを特徴とする方法。
  20. 請求項1から19までのいずれか1項に記載の方法において、
    前記トレーニングデータ(TD)からなるデータセットの部分集合(RD)により、または初期の重み付け(w)を有する、前記トレーニングデータ(TD)からなるすべてのデータセットによって前記方法を初期化する、
    ことを特徴とする方法。
  21. 請求項1から20までのいずれか1項に記載の方法において、
    前記密度推定器(DE)は、ニューラルクラウドアルゴリズムに基づく、
    ことを特徴とする方法。
  22. 請求項1から21までのいずれか1項に記載の方法において、
    つぎの繰り返しステップに対して選択されるかまたは重み付けされたデータセットの前記モデルエラー(ME)の少なくとも1つのあらかじめ設定した割合が、あらかじめ設定した閾値を下回る場合、および/またはつぎの繰り返しステップに対して選択されるかまたは重み付けされたデータセットの前記信頼度尺度(CM)の少なくとも1つのあらかじめ設定した割合が、あらかじめ設定した閾値を上回る場合に、前記繰り返しを終了する、
    ことを特徴とする方法。
  23. 技術システムを動作させる方法において、
    請求項1から22のいずれか1項に記載の方法によってデータドリブンモデルを形成し、引き続いて前記データドリブンモデルを使用し、前記技術システムの動作を解析す、および/または、当該技術システムの動作に影響を及ぼす、
    ことを特徴とする方法。
  24. 請求項7から10までのいずれか1項と組み合わせた請求項23に記載の方法において、
    前記開ループ制御モデルおよび/または閉ループ制御モデルによって前記技術システムを閉ループ制御および/または開ループ制御する
    ことを特徴とする方法。
  25. 請求項11または12に記載の方法において、
    強化学習に基づき、および/または、モデル予想制御に基づき、および/または最適化方法に基づき、前記監視モデルを計算する、
    ことを特徴とする方法
  26. 請求項11または12または25に記載の方法において、
    ニューラルネットワークに基づいて前記監視モデルを学習する、
    ことを特徴とする方法
  27. 請求項11または12または25または26と組み合わせた請求項23に記載の方法において、
    前記技術システムを前記監視モデルによって監視する、
    ことを特徴とする方法
  28. コンピュータプログラムであって、
    請求項1から27までのいずれか1項に記載の方法を
    計算機に実行させる
    ことを特徴とするコンピュータプログラム
JP2014513206A 2011-06-03 2012-06-01 殊にガスタービンまたは風力タービンのような技術システムのデータドリブンモデルを計算機支援で形成する方法 Active JP5832644B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
DE102011076936 2011-06-03
DE102011076936.6 2011-06-03
PCT/EP2012/060400 WO2012164075A2 (de) 2011-06-03 2012-06-01 Verfahren zur rechnergestützten generierung eines datengetriebenen modells eines technischen systems, insbesondere einer gasturbine oder windturbine

Publications (2)

Publication Number Publication Date
JP2014525063A JP2014525063A (ja) 2014-09-25
JP5832644B2 true JP5832644B2 (ja) 2015-12-16

Family

ID=46275807

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014513206A Active JP5832644B2 (ja) 2011-06-03 2012-06-01 殊にガスタービンまたは風力タービンのような技術システムのデータドリブンモデルを計算機支援で形成する方法

Country Status (5)

Country Link
US (1) US9466032B2 (ja)
EP (1) EP2697695B1 (ja)
JP (1) JP5832644B2 (ja)
CN (1) CN103597413B (ja)
WO (1) WO2012164075A2 (ja)

Families Citing this family (54)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9015093B1 (en) 2010-10-26 2015-04-21 Michael Lamport Commons Intelligent control with hierarchical stacked neural networks
US8775341B1 (en) 2010-10-26 2014-07-08 Michael Lamport Commons Intelligent control with hierarchical stacked neural networks
EP2706422B1 (de) * 2012-09-11 2016-07-27 Siemens Aktiengesellschaft Verfahren zur rechnergestützten Überwachung des Betriebs eines technischen Systems, insbesondere einer elektrischen Energieerzeugungsanlage
EP2971699B8 (en) 2013-03-15 2020-01-15 Rolls-Royce Corporation Lifing and performance optimization limit management for turbine engine
DE102013205356B4 (de) * 2013-03-26 2016-07-07 Siemens Aktiengesellschaft Verfahren zur rechnergestützten Steuerung und/oder Regelung eines technischen Systems
ES2568052T3 (es) * 2013-09-27 2016-04-27 Deutsche Telekom Ag Procedimiento y sistema para la evaluación de valores medidos registrados de un sistema
WO2015100632A1 (zh) 2013-12-31 2015-07-09 西门子公司 燃气轮机运行参数检测故障的诊断方法
JP6710039B2 (ja) * 2015-10-28 2020-06-17 三菱日立パワーシステムズ株式会社 計画装置、計画方法およびプログラム
US10551818B2 (en) * 2015-11-25 2020-02-04 United Technologies Corporation Fault detection methods and systems
FR3046265B1 (fr) * 2015-12-29 2018-08-10 Thales Systeme de surveillance d'une installation industrielle ; procedes de configuration et de surveillance associes
US10430716B2 (en) * 2016-02-10 2019-10-01 Ground Rounds, Inc. Data driven featurization and modeling
EP3260700A1 (en) * 2016-06-22 2017-12-27 Siemens Aktiengesellschaft A computer-based method and a prediction device for predicting operating data of a wind turbine
DE102017122928A1 (de) * 2016-10-11 2018-01-18 FEV Europe GmbH Verfahren zum Bestimmen eines Effizienzwertes einer Turbine mit einer variablen Geometrie
JP6815184B2 (ja) * 2016-12-13 2021-01-20 株式会社東芝 情報処理装置、情報処理方法、および情報処理プログラム
CN106704103B (zh) * 2017-01-05 2019-12-17 华北电力大学 一种基于叶片参数自学习的风电机组功率曲线获取方法
US9881430B1 (en) 2017-02-22 2018-01-30 General Electric Company Digital twin system for a cooling system
JP2018178900A (ja) * 2017-04-18 2018-11-15 株式会社日立製作所 風力発電システム
US10922627B2 (en) 2017-06-15 2021-02-16 Microsoft Technology Licensing, Llc Determining a course of action based on aggregated data
US10503580B2 (en) 2017-06-15 2019-12-10 Microsoft Technology Licensing, Llc Determining a likelihood of a resource experiencing a problem based on telemetry data
US10805317B2 (en) 2017-06-15 2020-10-13 Microsoft Technology Licensing, Llc Implementing network security measures in response to a detected cyber attack
US11062226B2 (en) 2017-06-15 2021-07-13 Microsoft Technology Licensing, Llc Determining a likelihood of a user interaction with a content element
KR102025145B1 (ko) * 2017-09-01 2019-09-25 두산중공업 주식회사 플랜트 데이터 예측 장치 및 방법
EP3489773A1 (de) * 2017-11-24 2019-05-29 Siemens Aktiengesellschaft Verfahren zum rechnergestützten steuern eines technischen systems, insbesondere einer energieerzeugungsanlage
JP7006364B2 (ja) * 2018-02-22 2022-01-24 富士通株式会社 強化学習プログラム、強化学習方法、および強化学習装置
JP7090440B2 (ja) * 2018-03-16 2022-06-24 株式会社Lixil 制御システム、学習装置、制御装置、及び制御方法
JP2019157841A (ja) * 2018-03-16 2019-09-19 株式会社Lixil 制御装置、及び制御方法
JP7090439B2 (ja) * 2018-03-16 2022-06-24 株式会社Lixil 制御システム、学習装置、制御装置、及び制御方法
JP7090445B2 (ja) * 2018-03-20 2022-06-24 株式会社Lixil 制御システム、学習装置、制御装置、及び制御方法
JP7090443B2 (ja) * 2018-03-20 2022-06-24 株式会社Lixil 制御システム、学習装置、及び制御方法
JP7090442B2 (ja) * 2018-03-20 2022-06-24 株式会社Lixil 制御システム、学習装置、及び制御方法
JP7090444B2 (ja) * 2018-03-20 2022-06-24 株式会社Lixil 制御システム、学習装置、及び制御方法
US20200219027A1 (en) * 2018-06-15 2020-07-09 Dalian University Of Technology Knowledge transfer-based modeling method for blast furnace gas scheduling systems
EP3588211A1 (de) * 2018-06-27 2020-01-01 Siemens Aktiengesellschaft Steuereinrichtung zum steuern eines technischen systems und verfahren zum konfigurieren der steuereinrichtung
WO2020142542A1 (en) * 2018-12-31 2020-07-09 Sentient Science Corporation Methods and systems for predicting risk of observable damage in wind turbine gearbox components
US11892003B2 (en) * 2019-03-06 2024-02-06 General Electric Company Application of machine learning to process high-frequency sensor signals of a turbine engine
JP7113958B2 (ja) * 2019-03-11 2022-08-05 三菱電機株式会社 運転支援装置及び運転支援方法
JP7345744B2 (ja) * 2019-07-08 2023-09-19 株式会社野村総合研究所 データ処理装置
EP4024295A4 (en) * 2019-09-30 2023-05-10 Siemens Aktiengesellschaft OPTIMIZATION METHOD AND DEVICE FOR INTEGRATED POWER SYSTEM AND COMPUTER READABLE STORAGE MEDIUM
EP4024145A4 (en) * 2019-09-30 2023-05-31 Siemens Aktiengesellschaft NONLINEAR MODEL MODELING METHOD, DEVICE, AND STORAGE MEDIA
EP3816885A1 (en) * 2019-10-28 2021-05-05 Siemens Gamesa Renewable Energy Innovation & Technology, S.L. Method for computer-implemented forecasting of wind phenomena with impact on a wind turbine
JP7381861B2 (ja) 2019-11-29 2023-11-16 富士通株式会社 複合モデル生成プログラム及び方法、並びに情報処理システム
EP3828408A1 (en) * 2019-11-29 2021-06-02 Siemens Gamesa Renewable Energy A/S A method and an apparatus for computer-implemented monitoring of a wind turbine
EP3839666B1 (de) 2019-12-18 2023-10-11 Siemens Aktiengesellschaft Steuereinrichtung zum steuern eines technischen systems und verfahren zum konfigurieren der steuereinrichtung
US11060504B1 (en) * 2020-02-07 2021-07-13 General Electric Company Systems and methods for continuous machine learning based control of wind turbines
CN113807000A (zh) * 2020-07-23 2021-12-17 广东毓秀科技有限公司 基于数据驱动的ups电池剩余寿命预测方法
CN111739134B (zh) * 2020-07-30 2023-05-26 腾讯科技(深圳)有限公司 虚拟角色的模型处理方法、装置及可读存储介质
CN111965981B (zh) * 2020-09-07 2022-02-22 厦门大学 一种航空发动机强化学习控制方法及系统
DE102021202287B4 (de) 2021-03-09 2023-03-23 Robert Bosch Gesellschaft mit beschränkter Haftung Verfahren und Vorrichtung zum Betreiben eines Kraftstoffeinspritzventils mithilfe maschineller Lernverfahren
DE102021202289A1 (de) 2021-03-09 2022-09-15 Robert Bosch Gesellschaft mit beschränkter Haftung Computerimplementiertes Verfahren und Vorrichtung zum Trainieren eines datenbasierten Zeitpunktbestimmungsmodells für die Bestimmung eines Öffnungs- oder Schließzeitpunkts eines Einspritzventils mithilfe maschineller Lernverfahren
DE102021113547A1 (de) * 2021-05-26 2022-12-01 Rwe Renewables Gmbh Verfahren zum Trainieren eines Machine-Learning-Modells verwendbar zum Bestimmen einer Restnutzungsdauer einer Windkraftanlage
US11649804B2 (en) 2021-06-07 2023-05-16 General Electric Renovables Espana, S.L. Systems and methods for controlling a wind turbine
US11860060B2 (en) 2022-04-05 2024-01-02 Rtx Corporation Integrally bladed rotor analysis and repair systems and methods
US20230314341A1 (en) * 2022-04-05 2023-10-05 Raytheon Technologies Corporation Data transfer systems and methods for analysis inputs
CN115493740B (zh) * 2022-11-14 2023-02-28 长江勘测规划设计研究有限责任公司 利用机器视觉测量水轮机内部流道压力脉动的方法及系统

Family Cites Families (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5485545A (en) * 1991-06-20 1996-01-16 Mitsubishi Denki Kabushiki Kaisha Control method using neural networks and a voltage/reactive-power controller for a power system using the control method
US5335291A (en) 1991-09-20 1994-08-02 Massachusetts Institute Of Technology Method and apparatus for pattern mapping system with self-reliability check
US5598075A (en) * 1995-09-13 1997-01-28 Industrial Technology Research Institute Servo control method and apparatus for discharging machine
GB0016561D0 (en) * 2000-07-05 2000-08-23 Rolls Royce Plc Health monitoring
US7321809B2 (en) * 2003-12-30 2008-01-22 The Boeing Company Methods and systems for analyzing engine unbalance conditions
US7216061B2 (en) * 2004-08-25 2007-05-08 Siemens Corporate Research, Inc. Apparatus and methods for detecting system faults using hidden process drivers
US8768664B2 (en) * 2005-03-18 2014-07-01 CMC Solutions, LLC. Predictive emissions monitoring using a statistical hybrid model
US7421348B2 (en) 2005-03-18 2008-09-02 Swanson Brian G Predictive emissions monitoring method
US8065022B2 (en) 2005-09-06 2011-11-22 General Electric Company Methods and systems for neural network modeling of turbine components
US7660636B2 (en) * 2006-01-04 2010-02-09 Accelerated Care Plus Corp. Electrical stimulation device and method for the treatment of dysphagia
DE102007001025B4 (de) * 2007-01-02 2008-11-20 Siemens Ag Verfahren zur rechnergestützten Steuerung und/oder Regelung eines technischen Systems
DE102007017259B4 (de) 2007-04-12 2009-04-09 Siemens Ag Verfahren zur rechnergestützten Steuerung und/oder Regelung eines technischen Systems
DE102007042440B3 (de) 2007-09-06 2009-01-29 Siemens Ag Verfahren zur rechnergestützten Steuerung und/oder Regelung eines technischen Systems
US8194938B2 (en) * 2009-06-02 2012-06-05 George Mason Intellectual Properties, Inc. Face authentication using recognition-by-parts, boosting, and transduction
CN101695190B (zh) 2009-10-20 2011-07-27 北京航空航天大学 一种基于神经网络的三维无线传感器网络节点自定位方法
JP5027859B2 (ja) 2009-10-26 2012-09-19 パナソニック デバイスSunx株式会社 信号識別方法および信号識別装置
CN101701845B (zh) 2009-11-04 2011-06-01 西安理工大学 一种机车车轮运行状态的识别方法
AU2011205557A1 (en) * 2010-01-14 2012-09-06 Venture Gain LLC Multivariate residual-based health index for human health monitoring
DE102010011221B4 (de) 2010-03-12 2013-11-14 Siemens Aktiengesellschaft Verfahren zur rechnergestützten Steuerung und/oder Regelung eines technischen Systems
CN101819411B (zh) 2010-03-17 2011-06-15 燕山大学 一种基于gpu的改进加权关联规则的设备故障预警与诊断方法
US20120072029A1 (en) * 2010-09-20 2012-03-22 Heatvu Inc. Intelligent system and method for detecting and diagnosing faults in heating, ventilating and air conditioning (hvac) equipment

Also Published As

Publication number Publication date
CN103597413A (zh) 2014-02-19
EP2697695B1 (de) 2019-05-08
RU2013157156A (ru) 2015-07-20
US20140100703A1 (en) 2014-04-10
WO2012164075A3 (de) 2013-02-21
JP2014525063A (ja) 2014-09-25
CN103597413B (zh) 2017-01-18
EP2697695A2 (de) 2014-02-19
US9466032B2 (en) 2016-10-11
WO2012164075A2 (de) 2012-12-06

Similar Documents

Publication Publication Date Title
JP5832644B2 (ja) 殊にガスタービンまたは風力タービンのような技術システムのデータドリブンモデルを計算機支援で形成する方法
US20210383041A1 (en) In-situ thermodynamic model training
Nian et al. A review on reinforcement learning: Introduction and applications in industrial process control
US20220146993A1 (en) Machine learning method and machine learning device for learning fault conditions, and fault prediction device and fault prediction system including the machine learning device
US11568310B2 (en) Apparatus for generating temperature prediction model and method for providing simulation environment
JP6513015B2 (ja) 機械の動作を制御する方法、および機械の動作を反復的に制御する制御システム
JP6896432B2 (ja) 故障予知方法、故障予知装置および故障予知プログラム
US11035346B2 (en) Method and device for predicting wind turbine fault
CN101923662B (zh) 信息处理设备、信息处理方法以及程序
US11604934B2 (en) Failure prediction using gradient-based sensor identification
JP7379833B2 (ja) 強化学習方法、強化学習プログラム、および強化学習システム
US8290885B2 (en) Information processing apparatus, information processing method, and computer program
CN110168220A (zh) 一种评估风力涡轮机发电机性能的方法和系统
EP3842635A1 (en) Operating a wind turbine with sensors implemented by a trained machine learning model
JP4028384B2 (ja) エージェント学習装置、方法、プログラム
Toubakh et al. Hybrid dynamic data-driven approach for drift-like fault detection in wind turbines
EP3985460B1 (en) Controller, control system and method to control a facility and to generate a control model by learning
US20230266721A1 (en) Method for configuring a control agent for a technical system, and control device
US20220108132A1 (en) Method and system for an adversarial training using meta-learned initialization
Ding et al. Diffusion world model
CN113906353A (zh) 用于对风力涡轮机的阻力系数进行计算机实现确定的方法
Toubakh et al. Advanced data mining approach for wind turbines fault prediction
CN111356959B (zh) 用于计算机辅助地控制技术系统的方法
Manrique Piramanrique et al. Data driven fault detection and isolation: a wind turbine scenario
US20240176311A1 (en) Method and apparatus for performing optimal control based on dynamic model

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150105

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150401

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: 20150928

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20151027

R150 Certificate of patent or registration of utility model

Ref document number: 5832644

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250