JP7223839B2 - 異常検出および/または予知保全のためのコンピュータ実装方法、コンピュータプログラム製品およびシステム - Google Patents

異常検出および/または予知保全のためのコンピュータ実装方法、コンピュータプログラム製品およびシステム Download PDF

Info

Publication number
JP7223839B2
JP7223839B2 JP2021512922A JP2021512922A JP7223839B2 JP 7223839 B2 JP7223839 B2 JP 7223839B2 JP 2021512922 A JP2021512922 A JP 2021512922A JP 2021512922 A JP2021512922 A JP 2021512922A JP 7223839 B2 JP7223839 B2 JP 7223839B2
Authority
JP
Japan
Prior art keywords
neural network
deep neural
intermediate output
data
observations
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
JP2021512922A
Other languages
English (en)
Other versions
JPWO2020049087A5 (ja
JP2022500745A (ja
Inventor
リカルド シェーグレン
ヨハン トリグ
Original Assignee
ザルトリウス ステディム データ アナリティクス アーベー
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 ザルトリウス ステディム データ アナリティクス アーベー filed Critical ザルトリウス ステディム データ アナリティクス アーベー
Publication of JP2022500745A publication Critical patent/JP2022500745A/ja
Publication of JPWO2020049087A5 publication Critical patent/JPWO2020049087A5/ja
Application granted granted Critical
Publication of JP7223839B2 publication Critical patent/JP7223839B2/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
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • 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
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/213Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/22Matching criteria, e.g. proximity measures
    • 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
    • G06N3/0442Recurrent networks, e.g. Hopfield networks characterised by memory or gating, e.g. long short-term memory [LSTM] or gated recurrent units [GRU]
    • 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/045Combinations of networks
    • G06N3/0455Auto-encoder networks; Encoder-decoder 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/0464Convolutional networks [CNN, ConvNet]
    • 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/08Learning methods
    • 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/092Reinforcement learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/60Type of objects
    • G06V20/69Microscopic objects, e.g. biological cells or cellular parts
    • G06V20/698Matching; Classification
    • 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/088Non-supervised learning, e.g. competitive learning

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Mathematical Physics (AREA)
  • Computational Linguistics (AREA)
  • Software Systems (AREA)
  • Biophysics (AREA)
  • Computing Systems (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Multimedia (AREA)
  • Automation & Control Theory (AREA)
  • Image Analysis (AREA)
  • Testing And Monitoring For Control Systems (AREA)
  • Investigating Or Analysing Biological Materials (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Description

本出願は、異常検出および/または予知保全のためのコンピュータ実装方法、コンピュータプログラム製品およびシステムに関し、特に、構造化データまたは非構造化データにおいて外れ値の検出を使用することに関する。
予知保全、サイバーセキュリティ、不正防止など、様々なシステムにおける異常検出がますます普及している。このようなシステムでは、システムの適切な機能を損なう可能性のある異常を、確実かつ迅速に検出することが重要である。
例えば、予知保全技術は、機器(例えば、機械や機械部品)の破損を防ぐために、いつ保守が必要かを自動的に判断することを目的としている。機器の劣化を自動的に見積もることにより、定期的または時間に基づく保守と比較して、大幅なコスト削減を実現できる。対照的に、固定された日程で計画すると、保守が必要以上に早く実行され、過剰な費用が発生するリスクがある。また、機器が予想よりも早く故障し、致命的な障害やプロセスの停止につながる可能性があるというリスクがある。
予測分析の方法論は、大きく2つの異なる方法論に分けることができる。まず、いわゆる教師なしのモデリング/機械学習を使用して、外部ラベルを必要とせずに、調査対象のシステムの通常の状態をモデル化できる。教師なしモデリングは、「機器は正常に動作しているかどうか」という質問に答えることができるが、システムを調査するためのさらなる手がかりは提供しない。より広く適用されている別のアプローチは、教師ありのモデリング/機械学習を使用することである。通常、履歴データが収集され、過去に発生した障害がドメインエキスパートによって調査され、障害の種類に分類され、場合によっては障害を修正するためにアクションが実行される。次に、予測モデルをトレーニングして、発生する障害の時期と種類を予測し、システムを修正する方法に関するガイダンスを提供する。後者の場合、人工ニューラルネットワークに基づくモデルを使用して、保守が必要な時期を予測できる。
人工ニューラルネットワークモデルは、予知保全の多くの用途で使用されている。完全には網羅されていないが、リストとしては、次のものが含まれる:
-ニューラルネットワークモデルを使用し、測定された物理パラメータに基づいて、物理システムの状態を分類する。つまり、事前に決定された動作条件(良好、不良、未決定など)に従って、システムを分類する(例:特許文献1~4を参照)。具体的なシステム例は、次のものが含まれる:
o石炭燃焼プラント(例:特許文献5を参照);
o金型鋳造(例:特許文献6を参照);
o電気自動車用のバッテリ製造(例:特許文献7を参照);
o石油およびガスの地表機器(例:特許文献8を参照);
o電気モータ(例:特許文献9を参照);
oリカレントニューラルネットワークを使用した鉄道車両の保守(例:特許文献10を参照);
o産業プラントの制御弁(例:特許文献11を参照);
-ニューラルネットワークを使用して、機器の残存寿命を予測する(例:特許文献12、13を参照)。具体的なシステム例は、次のものが含まれる:
o振動に基づくフライス盤のカッターヘッド(例:特許文献14を参照);
・汚染物質が決定された制限内にあるかどうかを判断するために、製造工場の汚染物質のソフトセンサとして、ニューラルネットワークモデルを使用する(例:特許文献15を参照);
システム状態を分類するための方法例として、ニューラルネットワークモデルを使用する(例:特許文献16~20を参照);
その他の例には、以下が含まれるが、これらに限定されない。
・振動監視に基づく軸受の残存寿命の予測(例:非特許文献1を参照);
ニューラルネットワークを含む3つの統計的回帰法を使用した水力発電ユニットの故障予測(例:非特許文献2を参照);
・軸受の残存寿命の予測(例:非特許文献3~5を参照);
・運動電流の徴候に基づく、モータ負荷の分類(例:非特許文献6を参照);
・風力タービンのギアボックス軸受の早期故障検出(例:非特許文献7、8を参照);
・赤外線サーモグラフィの使用に基づく、電気機器の熱欠陥の検出(例:非特許文献9を参照);
・振動信号と音信号に基づく、電気モータ駆動のウォームギアの故障分類(例:非特許文献10を参照);
・荷重、速度、ヨー角の特性に基づく、車輪とレールの摩耗予測(例:非特許文献11を参照);
・集光型太陽光発電所の管欠陥検出のための吸収体パイプの温度予測(例:非特許文献12を参照);
異常検出のためのニューラルネットワークの他の用途には、例えば、ネットワーク侵入や他の種類の悪意のある攻撃などのサイバーセキュリティ脅威を検出することが含まれる。例えば、ネットワークトラフィックは、機械学習技術を使用してモデル化し得る。例として、悪意のあるアクションおよび/またはアクタを示すトラフィックパターンの微妙な変化を検出し、侵入として分類し得る。もう1つの例は、金融取引における不正、疑わしいパターン、または特定傾向を検出することである。例えば、履歴データを使用して、機械学習に基づく分類器をトレーニングして、トランザクションまたはアクタ(例えば、不正または信頼できる)を分類したり、傾向を識別したりできる。
上記用途に共通するのは、システムを特徴付けるために蓄積されたデータ(例えば、履歴データなど)を使用することと、正常動作と既知の種類の異常とを区別するために機械学習モデルを使用することである。機械学習モデルは、深層学習モデルであり得る。深層ニューラルネットワークを使用した深層学習は、深層ニューラルネットワークによって学習された強力な変換のおかげで、多くの用途で非常に普及している。深層ニューラルネットワークは、入力層、出力層、および入力層と出力層の間に複数の隠れ層を有する人工ニューラルネットワークであり得る。深層ニューラルネットワークを使用した深層学習モデルの用途の長いリストが、そのようなモデルが幅広い用途において教師あり予知保全と異常検出にうまく使用できることを示している。しかし、深層学習やニューラルネットワークを含むすべてのデータ駆動型モデリングの問題の1つは、モデルのトレーニング中に遭遇する観測値とは異なる観測値を処理できないことである。つまり、予測中に外れ値を処理することはできない。外れ値(外れ値の観測値)は、モデルのトレーニングに使用されるデータとは、大幅に異なる観測値であり得る。外れ値は、モデルが以前に遭遇したものとは大幅に異なるため、このような状況ではモデル動作が定義されていない、すなわち、予測を信頼することができない。外れ値が検出されたときのモデル動作が定義されていないので、以前に不明であった種類のエラーは、検出されないことがある。
ただし、自動運転車、診断アプリケーション、発電所制御システムなどの重要なシステムで深層学習を採用する場合は、未知の種類の異常にもフラグを立てることが極めて重要になり得る。さらに、深層学習システムによって提供される予測を信頼するべきでないことを、理解することが重要になり得る。言い換えれば、深層学習システムによって提供される予測が正確であるかどうかだけでなく、深層学習システムが、そもそも予測する必要があるかどうかを知る必要ことが重要になり得る。学習された表現の限界を理解し、新しい種類のデータまたは深層学習システムによって認識されないデータに遭遇したときを認識することによって、自律的な決定を改善することができる。
深層ニューラルネットワークが予測の不確実性を記述できるようにするために、多くの異なる方法が提案されている。これらの方法の多くは、予測中にポイント推論ではなくベイズ推論を実行するように深層ニューラルネットワークをトレーニングすることに基づいている。これは、単一ポイントを予測するのではなく、ネットワークが可能なポイント分布を予測するようにトレーニングされていることを意味する。直感的には、モデルで説明できない観測値は、分布が広いか、不確実性が大きいことになる。ベイズ予測を実行するための一般的な解決策の1つは、いわゆるモンテカルロドロップアウト(MCドロップアウト)であり、Y.GalおよびZ.Gharamaniによる「Dropout as a Bayesian approximation:Representing model uncertainty in deep learning」(機械学習の国際会議、2016年、第1050~1059頁)で、開示されている。MCドロップアウトは、「ドロップアウト」を使用してトレーニングされたニューラルネットワークに適用できる。これは、トレーニングフェーズ中のニューラルネットワーク内の一部のユニット(例、ニューロン、ノード)をドロップアウトする(または、言い換えると、無視する)ことで、ニューラルネットワークの過剰適合を減らすための正則化手法である。MCドロップアウトは、予測のためにいわゆるモンテカルロサンプリングを使用して、予測ごとに複数の推論を行い、予測時間中にネットワークニューロンをランダムに排除またはドロップアウトする。様々な予測は、ドロップアウトによって異なり、予測分布を表す。次に、結果の分布のエントロピーまたは分散を計算することにより、予測の不確実性を定量化できる。
ドロップアウトを使用してモンテカルロサンプルを生成する他の方法には、バッチ正規化パラメータに基づくサンプリング(M.Teye、H.Azizpour、およびK.Smithによる「Bayesian Uncertainty Estimation for Batch Normalized Deep Networks」、アーカイヴ180206455Stat、2018年2月)、アンサンブル内の様々なモデルからの予測(B.Lakshminarayanan、A.Pritzel、およびC.Blundellによる「Simple and scalable predictive uncertainty estimation using deep ensembles」、ニューラル情報処理システムの進歩、2017年、第6405~6416頁)、共有に基づくネットワーク内の複数の「予測ヘッド」(I.Osband、C.Blundell、A.Pritzel、およびB.Van Royによる「Deep exploration via bootstrapped DQN」、ニューラル情報処理システムの進歩、2016年、第4026~4034頁、E.Ilgらによる「Uncertainty Estimates for Optical Flow with Multi-Hypotheses Networks」、アーカイヴ180207095Cs、2018年2月)、通常点の重みではなく重み分布の変分推論(C.Blundell、J.Cornebise、K.Kavukcuoglu、およびD.Wierstraによる「Weight uncertainty in neural networks」、アーカイヴ・プレプリント、アーカイヴ150505424、2015年)、および、既存の重みからの分布のラプラス近似(H.Ritter、A.Botev、およびD.Barberによる「A Scalable Laplace Approximation for Neural Networks」、2018年2月)、を含む。
ただし、これらの方法は、使用されるモデルの種類を前提とし、ネットワークの構築またはトレーニングの方法を制約したり、および/または、予測ごとに複数の推論に依存したりする可能性がある。これにより、複数の推論を行うことが実用的でないリアルタイムシステム、および/または、制約に適合しない既存システムへの一般的な適用性が制限される可能性がある。
US2019/0137985A1 EP0633536A1 US5566092A JP5643387B2 CN109255440A KR101677358B1 CN106168799B US8676721B2 CN108921303A EP3179326A1 US2016/0071004A1 DE10161633A1 US7457785B1 US8781982B1 EP0712509B1 EP2871803B1 US8275642B2 FR2698704A1 US6298308B1 JP5043839B2
Wu S、Gebraeel N、Lawley MA、YihY.Aによる「Neural Network Integrated Decision Support System for Condition-Based Optimal Predictive Maintenance Policy」、IEEEシステム、人間、サイバネティックス上のトランザクション、Part Syst Hum、2007年3月、第37号(2)、第226~36頁 Lucifredi A、Mazzieri C、Rossi Mによる「APPLICATION OF MULTIREGRESSIVE LINEAR MODELS,DYNAMIC KRIGING MODELS AND NEURAL NETWORK MODELS TO PREDICTIVE MAINTENANCE OF HYDROELECTRIC POWER SYSTEMS」、機械システムと信号処理、2000年5月1日、第14号(3)、第471~94頁 Huang R、Xi L、Li X、Richard Liu C、Qiu H、Lee Jによる「Residual life predictions for ball bearings based on self-organizing map and back propagation neural network methods」、機械システムと信号処理、2007年1月1日、第21号(1)、第193~207頁 Shao Y、Nezu Kによる「Prognosis of remaining bearing life using neural networks」、英国機械学会の論文講演集、Part J Sys Control Eng、2000年5月1日、第214号(3)、第217~30頁 Vachtsevanos G、Wang Pによる「Fault prognosis using dynamic wavelet neural networks」、2001IEEEオートテストコン論文講演集IEEEシステム即応性技術会議(Cat No01CH37237)2001年、第857~70頁 Bansal D、Evans DJ、Jones Bによる「A real-time predictive maintenance system for machine systems」、機械ツールと製造の国際学術誌、2004年6月1日、第44号(7)、第759~66頁 Bangalore P,Letzgus S,Karlsson D,Patriksson Mによる「An artificial neural network-based condition monitoring method for wind turbines,with application to the monitoring of the gearbox」、風力エネルギー、2017年、第20号(8)、第1421~38頁 Bangalore P,Tjernberg LBによる「An Artificial Neural Network Approach for Early Fault Detection of Gearbox Bearings」、IEEE Transactions Smart Grid、2015年3月、第6号(2)、第980~7頁 Huda ASN,Taib Sによる「Application of infrared thermography for predictive/preventive maintenance of thermal defect in electrical equipment」、応用熱技術、2013年11月3日、第61号(2)、第220~7頁 Waqar T、Demetgul Mによる「Thermal analysis MLP neural network based fault diagnosis on worm gears」、メジャメント、2016年5月1日、86号、第56~66頁 Shebani A、Iwnicki Sによる「Prediction of wheel and rail wear under different contact conditions using artificial neural networks」、Wear、2018年7月15日、406~407号、第173~84頁 Jimenez AA、Munoz CQG、Marquez FPG、Zhang Lによる「Artificial Intelligence for Concentrated Solar Plant Maintenance Management」、Xu J、Hajiyev A、Nickel S、Gen M、editors、経営科学と工学管理に関する第10回国際会議の講演論文集、シュプリンガー、シンガポール、2017年、第125~34頁(インテリジェントシステムとコンピューティングの進歩) Saxena A、Goebel K、Simon D、Eklund Nによる「Damage propagation modeling for aircraft engine run-to-failure simulation」、2008年の予後と健全性管理に関する国際会議、2008年、第1~9頁 Gers FA、Schmidhuber J、Cummins Fによる「Learning to Forget:Continual Prediction with LSTM」、ニューロコンピューティング、1999年、第12号、第2451~2471頁 Hochreiter S,Schmidhuber Jによる「Long short-term memory」、ニューロコンピューティング、1997年、第9号(8)、第1735~1780頁 Kingma DP、Ba J.Adamによる「A Method for Stochastic Optimization」、アーカイヴ14126980Cs[Internet]、2014年12月22日[2018年7月3日引用]http://arxiv.org/abs/1412.6980から入手可能
本出願は、上記の課題の1つまたは複数に対処する。したがって、一態様によれば、課題は、既知の種類のイベントだけでなく、これまで知られていなかった未知の種類のイベントも処理できる、異常検出および/または予知保全のための改善された方法およびシステムに関する。さらなる課題は、深層ニューラルネットワークへの観測値入力に対して深層ニューラルネットワークによって与えられる予測の信頼性を効率的に評価することである。
本課題は、独立した請求項によって開示された特徴によって解決される。さらなる例示的な実施形態は、従属した請求項によって定義される。
一態様によれば、対象のエンティティ(例えば、監視されたエンティティ)における異常検出のためのコンピュータ実装方法が提供される。
この方法は、
新しい観測値を受信することであって、前記新しい観測値は、エンティティの少なくとも1つのパラメータを特徴付けることと、
複数の隠れ層を有する深層ニューラルネットワークに新しい観測値を入力し、深層ニューラルネットワークに入力できる可能性のある観測値を含むトレーニングデータセットを使用してトレーニングされることと、
受信した新しい観測値を深層ニューラルネットワークに入力することであって、深層ニューラルネットワークの複数の隠れ層の少なくとも1つから出力される中間出力値の第2のセットを取得することと、
記憶媒体に格納された潜在変数モデルを使用して、中間出力値の第2のセットを射影値の第2のセットにマッピングすることと、
潜在変数モデルおよび射影値の第二のセットに基づいて、受信した新しい観測値がトレーニングデータセットに関して、外れ値であるかどうかを判断することとを含み、
記憶媒体に格納された潜在変数モデルは、
深層ニューラルネットワークの複数の隠れ層のうちの前記1つから出力される中間出力値の第1のセットを取得することであって、中間出力値の第1のセットのそれぞれは、トレーニングデータセットの少なくとも一部に含まれる可能性のある観測値の異なる1つを入力することによって取得されることと、
中間出力値の第1のセットを使用して潜在変数モデルを構築/適合することであって、潜在変数モデルは、中間出力のセットの次元よりも低い次元を有する潜在変数モデルの部分空間内の射影値の第1のセットへ中間出力値の第1のセットのマッピングすることと、
によって構築される。
別の態様によれば、対象のエンティティにおける異常検出のためのコンピュータ実装方法が提供される。この方法は、
深層ニューラルネットワークおよび深層ニューラルネットワークをトレーニングにするために使用されるトレーニングデータセットの少なくとも一部を取得することであって、深層ニューラルネットワークは、複数の隠れ層を含み、トレーニングデータセットを使用してトレーニングされ、トレーニングデータセットは、深層ニューラルネットワークに入力できる複数の可能性のある観測値が含むことと、
複数の隠れ層のうちの少なくとも1つから出力される中間出力値の第1のセットを取得することであって、中間出力値の第1のセットのそれぞれは、トレーニングデータセットの少なくとも一部に含まれる可能性のある観測値の異なる1つを入力することによって取得することと、
中間出力値の第1のセットを使用して潜在変数モデルを構築(つまり適合)することであって、潜在変数モデルは、中間出力のセットの次元よりも低い次元を有する潜在変数モデルの部分空間内の射影値の第1のセットへの中間出力値の第1のセットのマッピングを提供することと、
深層ニューラルネットワークに入力される新しい観測値を受信することと、
受信した新しい観測値を深層ニューラルネットワークに入力することにより、深層ニューラルネットワークの複数の隠れ層のうちの前記少なくとも1つから出力される中間出力値の第2のセットを取得することと、
潜在変数モデルを使用して、中間出力値の第2のセットを射影値の第2のセットにマッピングすることと、
潜在変数モデルおよび射影値の第2のセットに基づいて、受信した観測値がトレーニングデータセットに関して外れ値であるかどうかを判断することと、
を含む。
受信した観測値が外れ値であると判断された場合、受信した観測値に対して深層ニューラルネットワークによって与えられた予測(例えば、出力)は、新しい種類の異常であると仮定され得、および/または、受信した観測値が外れ値ではないと判断された場合に、予測は信頼性が低いと、仮定され得る。したがって、上記の態様のいずれかによる方法は、
深層ニューラルネットワークによって、新しい観測値の予測を計算することと、
予測および新しい観測値が外れ値であるかどうかを判断に基づいて、エンティティに少なくとも1つの異常が発生したことを示す結果を判断することと、
をさらに含み得る。
さらに別の態様によれば、対象のエンティティの予知保全のためのコンピュータ実装方法が提供される。この方法は、
観測されるエンティティに対して実行された少なくとも1つのセンサ測定からのセンサデータに基づいて新しい観測値を取得することと、
本明細書に記載の様々な態様、例、および実施形態の方法に従って、エンティティの少なくとも1つの異常の発生を示す結果を判断することと、
少なくとも1つの異常の発生を示す結果に基づいて、修正処置を講じる必要があるかどうか、および、任意選択で講じるべき修正処置の種類を、判断することと、
を含む。
対象のエンティティは、任意の種類のエンティティ、例えば、機器(機械または機械部品など)、または、複数の機器を含むより複雑な技術システムであり得る。観測されるエンティティの非限定的な例は、バイオリアクタ、発酵槽、発電所、産業プラント、車両、石油、ガス、太陽エネルギーまたは風力エネルギー機器、モータ、バッテリ、制御弁、ギアなどである。観測されるエンティティはまた、コンピュータまたはコンピュータネットワークなどのコンピュータシステム、例えば、金融取引を可能にするコンピュータシステム/ネットワークであり得る。
異常は、対象のエンティティの通常操作および/または通常状態からの逸脱であり得る。例えば、異常は、センサデータ、ネットワークトラフィック量および/またはパターン、ユーザの行動および/または生体認証特性などの逸脱に関連し得る。
観測値は、観測されるエンティティの少なくとも1つのパラメータを特徴付けるまたは示す、任意のデータであり得る。例えば、データは、1つまたは複数の機器の特性の状態を測定する1つまたは複数のセンサによって取得されたデータ、ネットワークトラフィックを特徴付けるデータ(例えば、ネットワーク上の異なるコンピュータ間で転送されたデータパッケージのログ)、金融取引に関連するデータなどの他のトランザクションデータであり得る。測定は、定期的または実質的に連続的で、例えば、リアルタイム測定であり得る。特定のアクション(例えば、ユーザとコンピュータの相互作用)の結果として測定が開始される可能性もある。データはまた、異なる入力情報、例えば、異なるセンサからの測定の組合せまたは融合に基づくことができる。
深層ニューラルネットワークによって計算される予測には、異常が検出されたかどうか、または発生する可能性が高いかどうかの予測を含み得る。任意選択で、異常が検出された、または発生する可能性が高いと判断された場合、予測には、検出された、または発生する可能性のある異常の種類を含み得る。
予測は、対象のエンティティの状態または状況を示す1つまたは複数のパラメータに関連するか、またはそれを含み得る。1つまたは複数のパラメータは、物理的パラメータ(温度、速度など)、残存寿命、摩耗、少なくとも1つの機器の故障または損傷の確率、汚染または他の望ましくない状況のレベル、不正または攻撃などの確率を含み得る。
予測はまた、良好/不良/未決定、安全/安全でない/脆弱な、信頼できる/不正であるなど、所定の動作条件または基準に従って、対象のエンティティの状態または状況の分類に関連するか、またはそれを含み得る。
新しい観測値の予測と、新しい観測値が外れ値であるかどうかの判断の結果に基づいて、エンティティに少なくとも1つの異常が発生したことを示す結果が、判断または生成される。例えば、新しい観測値が外れ値ではないと判断された場合、予測は、少なくとも1つの異常の発生を示す結果(簡単には、異常検出結果)として、受け入れられ出力され得る。さらに、異常検出結果を判断するために、受け入れられた予測をさらに処理することが可能である。異常検出の結果に基づいて、さらなるアクション(修正処置など)が開始され得る。さらなるアクションは、特定のアプリケーションに依存し得る。例えば、予測管理の場合、予測/異常検出結果の種類に応じて、必要かどうか、および、任意選択で、どの修正保守アクションが必要かを判断し得る。サイバーセキュリティまたはトランザクション・セキュリティアプリケーションの場合、アプリケーションおよび予測/異常検出結果の種類に応じて、コンピュータ、ユーザ、アカウントなどがブロックされるか、または、その他の予防アクションを講じ得る。
例えば、予測/異常検出の結果が異常はないことを示している場合、または特定の種類の異常が検出された場合、観測されるデバイスまたはシステムの通常の動作、整合性、またはセキュリティに対する脅威ではないと見なされるので、さらなるアクションは、システムの操作および/または監視を継続し、ユーザ、コンピュータ、サーバなどへのアクセスを許可し、ユーザ、コンピュータ、サーバなどが、トランザクションを開始および/または完了することになり得る。予測/異常検出の結果に異常があることを示している場合、修正保守アクション、エンティティへのアクセスの拒否または制限などの修正アクションを講じ得る。修正処置は、予測よる指定通りに検出された異常の種類に依存する可能性があり、修正アクションの性質は、特定のアプリケーションに依存し得る。修正処置は、例えば、修正保守アクションのスケジューリング、対象のエンティティに関連する少なくとも1つの制御パラメータの変更、対象のエンティティまたはそれに含まれる少なくとも1つの機器の停止、または、対象のエンティティ(例えば、特定のコンピュータ、ユーザ、サーバなど)へのアクセスの拒否を含み得る。また、少なくとも1つのさらなるアクションは、観測されるシステムの状態または状況に関連する情報を表示すること、警告を発することなどを含み得る。少なくとも1つのさらなるアクションは、自動的に判断および/または開始され得るか、または、手動介入またはユーザからのさらなる入力を必要とし得る。
新しい観測値が外れ値であると判断された場合、異常検出結果は、新しい観測値が新しい種類または未知の種類の異常に関連していることを示し得る。異常検出結果には、外れ値を続行または破棄する前に、外れ値を手動で分析するようにユーザに警告することも含まれ得る。例えば、外れ値の観測値が監視システムの障害、例えばセンサ監視装置の故障によって引き起こされたと結論付けられた場合、システムオペレータは、例えば、故障したセンサを交換することによって、監視システム自体が修正されるように予定を立てることができる。外れ値の観測値は、履歴データでは観測されない珍しいまたは新しい異常、例えば、新しいまたは珍しい方法での機器の破損、コンピュータネットワークへの侵入を実行する新しい方法、または、不正なトランザクションを実行する新しい種類が、原因であり得る。このような場合、外れ値の観測値には、さらに調査した後、例えば手動調査などの新しい種類の異常クラスが割り当てられ、それを使用して、将来の監視で新しい異常種類を含む最新の異常検出モデルをトレーニングし得る。
そこに記載されている様々な実施形態および例において、監視あり予知保全の能力は、モデルトレーニング後に外れ値を検出する能力を予測深層ニューラルネットワークに与えることによって、教師なしモデリングと組み合わせることができる。このアプローチを使用すると、深層ニューラルネットワークを使用して異常と障害種類を分類すると同時に、これまで未知の異常や障害種類が発生した場合でも、定義された動作を行うことができる。
本明細書に記載の様々な実施形態および例では、深層ニューラルネットワークによって処理されるデータは、構造化データまたは非構造化データであり得る。非構造化データは、事前定義されたデータモデルがないか、または、事前定義された方法で編成されていないデータとして理解され得る。非構造化データは、内部構造を有する場合があるが、事前定義されたデータモデルまたはスキーマを介して構造化されていない。非構造化データの例には、画像、音声、ビデオ、センサデータ、およびテキストが含まれるが、これらに限定されない。データは、生データ(未処理データ)である場合もあれば、正規化、ノイズ除去、コーディング、コントラスト強調、エッジ強調、色または周波数変換、サイズ変更、トリミング、および/またはその他の既知のデータ処理操作などの、なんらかの前処理が行われたデータであり得る。
いくつかの特定の例では、深層ニューラルネットワークによって処理されるデータは、温度センサ、加速度センサ、圧力センサ、ひずみゲージ、力センサ、マイク、画像センサなどの1つまたは複数のセンサから取得されたセンサデータ(例えば、温度、圧力、電流、磁束など)であり得る。テストデータセットに含まれる可能性のある観測値および受信した観測値のそれぞれは、ある時点または一定期間中に1つまたは複数のセンサによって取得されたセンサデータを含み得る。
センサは、機器の状態が物理的測定に基づいて間接的に計算される、いわゆるソフトセンサとも呼ばれる。例示的なソフトセンサには、カルマンフィルタ、速度推定器、製品品質(例えば、食品品質)の推定器などが含まれる。センサは、音声または画像の形式で、非構造化データを収集する場合もある。データはまた、ネットワークトラフィックを特徴付けるデータ(例えば、ネットワーク上の異なるコンピュータ間で転送されたデータパッケージのログ)、または、金融取引に関連するデータなどの他の取引データであり得る。
本開示では、画像は、ピクセルの2次元配置を含み得る。各ピクセルは、少なくとも1つの値を含み得る。例えば、グレースケール画像のピクセルは、ピクセルの強度を示す1つの値を含み得る。カラー画像のピクセルは、RGB色空間などの色空間の座標を示す複数の値、例えば3つの値を含み得る。ピクセルはまた、異なる画像化モダリティの組合せ、例えば、位相差画像化および蛍光画像化などの2つの異なる画像化装置からの1つまたは複数の強度値の組合せを含み得る。画像が深層ニューラルネットワークによって処理される場合、テストデータセットに含まれる各観測値と受信された観測値は、画像を含み得る。
本開示において、「深層ニューラルネットワーク」という用語は、入力層、出力層、および入力層と出力層との間に提供される複数の隠れ層を有する人工ニューラルネットワークとして理解し得る。ニューラルネットワークの隠れ層は、「中間層」と呼ばれ得る。したがって、隠れ層からの出力は「中間出力」と呼ばれ得る。さらに、隠れ層からの出力は、複数の値を含み得、それらのそれぞれは、隠れ層に含まれるノードに対応し得る。したがって、本明細書で使用される「中間出力値のセット」という用語は、隠れ層のそれぞれのノードから出力される複数の値を含む、隠れ層の出力を示し得る。
使用される深層ニューラルネットワークの種類は特に限定されない。いくつかの例では、本開示における「深層ニューラルネットワーク」は、複数の隠れ層を有するフィードフォワードニューラルネットワークであり得る。フィードフォワードニューラルネットワークでは、ノード間の接続は、サイクルを形成しない。フィードフォワード深層ニューラルネットワークの特定の例は、視覚画像の分析に一般的に適用される畳み込みニューラルネットワーク(CNN)であり得る。
他のいくつかの例では、本開示における「深層ニューラルネットワーク」は、ノード間の接続がシーケンスに沿って有向グラフを形成するリカレントニューラルネットワークであり得る。リカレント深層ニューラルネットワークの特定の例は、シーケンシャルデータを処理できる長・短期記憶(LSTM)であり得る。
他の例示的な深層ニューラルネットワークは、リカーシブニューラルネットワークおよびトランスフォーマニューラルネットワークである。
本開示では、「潜在変数モデル」という用語は、観測可能な変数のセットを潜在変数のセットに関連付けるか、またはマッピングする統計モデルであり得る。本明細書に記載の様々な実施形態および例では、(第1または第2の)「中間出力値のセット」は、潜在変数モデルの観測可能な変数のセットと見なすことができる。さらに、本明細書に記載の様々な実施形態および例では、「射影値のセット」は、潜在変数モデルの潜在変数のセットと見なすことができる。
いくつかの例では、上記の態様による方法において、潜在変数モデルは、主成分分析に従って構築または適合され得る。
本開示では、PCAとも呼ばれる「主成分分析」という用語は、直交変換を使用して、相関する可能性のある変数の観測値のセットを、主成分と呼ばれる、線形無相関の変数の値のセットに変換する統計的手順を示すと理解することができる。
PCAの改良型には、スパースPCA(Jenatton,R.、Obozinski,G.およびBach,F.(2010年3月)、「Structured sparse principal component analysis」、人工知能と統計に関する第13回国際会議の講演論文集、(第366~373頁))が含まれるが、これらに限定されない。PCAの代替には、ランダムな射影(Dasgupta,S.(2000年6月)、「Experiments with random projection」、人工知能の不確実性に関する第16回会議の講演論文集(第143~151頁)、Morgan Kaufmann Publishers Inc.)、スパース・ランダム射影(Achlioptas,D.(2003年)、「Database-friendly random projections:Johnson-Lindenstrauss with binary coins」、コンピュータとシステムサイエンス学会誌、第66号(4)、第671~687頁)、非常にスパースなランダム射影(Li,P.、Hastie,T.J.およびChurch,K.W.(2006年8月)、「Very sparse random projections」、知識発見とデータマイニングに関する第12回ACM SIGKDD国際会議の講演論文集(第287~296頁)、ACM)、および、自己組織化マップ(Kohonen,T.(1998)、「The self-organizing map」、ニューロコンピューティング、第21号(1~3)、第1~6頁)、が含まれるが、これらに限定されない。
他のいくつかの例では、上記の態様による方法では、潜在変数モデルは、オートエンコーダを使用して構築または適合され得る。「オートエンコーダ」は、教師なしの方法で効率的なデータ・コーディングを学習するために使用される一種の人工ニューラルネットワークであり得る。オートエンコーダの1つの形態は、入力層、出力層、およびそれらを接続する1つまたは複数の隠れ層を有し、出力層が入力層と同じ数のノードを有し、独自の入力を再構築する目的を持つ、フィードフォワード、非リカレントニューラルネットワークであり得る。様々な種類のオートエンコーダには、変化型オートエンコーダ(Kingma,D.P.およびWelling,M.(2013年)、「Auto-encoding variational bayes」、アーカイヴプレプリント、アーカイヴ1312.6114)、および、スパースオートエンコーダ(Makhzani,A.およびFrey,B.(2013年)、「K-sparse autoencoders」、アーカイヴプレプリント、アーカイヴ1312.5663)が含まれるが、これらに限定されない。
潜在変数モデルを使用して、潜在変数モデルまでの新しい観測値の距離、言い換えれば潜在変数近似からの距離を決定し得る。決定された距離に基づいて、新しい観測値が外れ値であるかどうかに関する判断を行うことができる。受信した新しい観測値が外れ値であるかどうかに関する判断は、例えば、潜在変数モデルまでの新しい観測値の距離が閾値距離よりも大きいかどうかを判断することを含み得る。新しい観測値の計算された距離が閾値距離よりも大きい場合、新しい観測値が外れ値であると判断することができる。
潜在変数モデルを使用して、トレーニングデータセットの少なくとも一部の各観測値から潜在変数モデルまでの距離を計算し、決定された距離に基づいて閾値距離を決定することで、閾値距離を決定し得る。したがって、上記の態様のいずれかによる方法は、複数の距離に基づいて閾値距離を決定することをさらに含み得、そのそれぞれは、射影値の第1のセットの分布に関して、射影値の第1のセットの異なる1つについて計算される。
さらに、上記の態様による方法において、受信された観測値が外れ値であるかどうかを判断する前記ステップは、
射影値の第1のセットの分布に関して、射影値の第2のセットの距離を計算することと、
計算された距離が距離の閾値よりも大きい場合、受信した観測値がトレーニングデータセットに関して外れ値であると判断することと、を含み得る。
距離の閾値は、距離に基づいて決定することができ、距離のそれぞれは、射影値の第1のセットの分布に関して、射影値の第1のセットの異なる1つについて計算し得る。
距離指標は、潜在変数近似からの距離を定量化するのに適した任意の距離指標(すなわち、射影値の第1のセット)であり得る。例えば、距離は、残差平方和(RSS)、マハラノビス距離、局所外れ値因子、またはLOFであり得る(例えば、M.M.Breunig、H.-P.Kriegel、R.T.NgおよびJ.Sanderの「LOF:Identifying Density-based Local Outliers」、2000ACM SIGMODデータ管理の国際会議の講演論文集、ニューヨーク市、ニューヨーク州、USA、2000年、第93~104頁)。距離はまた、記述された距離の2つ以上を組み合わせることによって形成される複合距離指標に基づく複合距離であり得る。
決定された距離は、特にマハラノビス距離であり得る。本開示において、「マハラノビス距離」という用語は、点P(例えば、観測値に対応する)と分布Dとの間の距離の既知の尺度を示すと理解され得る。より具体的には、「マハラノビス距離」は、多変量確率分布の原点から、観測値がどの程度の標準偏差で離れているかを測定する。
決定された距離は、残差平方和(RSS)であり得る。例えば、上記の態様による方法では、受信した観測値が外れ値であるかどうかを判断する前記ステップは、
潜在変数モデルと射影値の第2のセットを使用して、中間出力値の第2のセットに対応する中間出力値の近似セットを決定することと、
中間出力値の第2のセットと中間出力値の近似セットの二乗近似残差を計算することと、
計算された二乗近似残差が二乗近似残差の閾値よりも大きい場合、受信した新しい観測値がトレーニングデータセットに関して外れ値であると判断することと、
を含み得る。
二乗近似残差の閾値は、二乗近似残差に基づいて決定し得、そのそれぞれは、中間出力値の第1のセットの異なる1つ、および中間出力値の第1のセットの前記1つに対応する中間出力値の近似セットについて計算し得る。
本開示において、「二乗近似残差」という用語は、残差平方の合計を示すものとして理解され得、ここで、残差は、観測値と潜在変数モデルによって提供される近似値との間の差である。
さらに、上記の態様による方法では、中間出力値の第1のセットを取得し、潜在変数モデルを構築するステップは、複数の隠れ層のうちの2つ以上に対して実行され得、
中間出力値の第2のセットを取得し、中間出力値の第2のセットを射影値の第2のセットにマッピングするステップは、複数の隠れ層のうちの前記2つ以上に関して実行され得、
受信した観測値が外れ値であるかどうかを判断するステップは、潜在変数モデルおよび、複数の隠れ層のうちの前記2つ以上に関して取得された射影値の第2のセットに基づいて実行され得る。
さらに、上記の態様による方法では、深層ニューラルネットワークの複数の隠れ層のうちの少なくとも1つから出力される中間出力値を取得するステップは、複数の隠れ層の少なくとも1つからの活性化を判断すること、および任意選択で、判断された活性化にグローバルプーリングを適用することを含み得る。活性化は、例えば、活性化ベクタの形態であり得る。
与えられた深層ニューラルネットワークのレイヤi(i=1、2、3、...)からの特定の観測値xの活性化ベクトルa(または活性化)は、レイヤiからの中間出力に対応し、観測値が所与のネットワークに入力されたときに、レイヤiのそれぞれのノードからの出力に対応する値を持つ成分を含み得る。深層ニューラルネットワークの様々なレイヤ上の活性化または活性化ベクトルのそれぞれは、入力データ(つまり、入力観測値)の特徴表現を提供し得る。言い換えると、活性化は、深層ニューラルネットワークの様々なレイヤで、入力データの変換された、または前処理された表現を提供し得る。
特に、トレーニングデータセットの観測値に基づいて、深層ニューラルネットワークの複数の隠れ層の少なくとも1つから出力される第1の中間出力値を取得するステップは、複数の隠れ層の少なくとも1つからの前記観測値の活性化を判断することと、および任意選択で、判断された活性化にグローバルプーリングを適用することと、を含み得る。同様に、新しい観測値に基づいて深層ニューラルネットワークの複数の隠れ層の少なくとも1つから出力される第2の中間出力値を取得するステップは、複数の隠れ層の少なくとも1つからの新しい観測値の活性化を判断すること、および任意選択で、判断された活性化にグローバルプーリングを適用することを含み得る。グローバルプーリングは、グローバル平均プーリングであり得る。
例えば、活性化が特徴マップからのものである場合、グローバルプーリング(例えば、グローバル平均プーリング)が各特徴マップ全体に適用され得る。潜在変数モデルは、その後、(任意選択でプールされた)活性化、つまり(任意選択でプールされた)トレーニングセット活性化、および深層ニューラルネットワークへの新しい観測値入力に対して計算された(任意選択でプールされた)活性化に、構築または適合され得る。
グローバル平均プーリング(GAP)は、ネットワーク活性化の特徴マップ表現をベクトル表現に変換する手順を指す(例えば、Lin、Min、Qiang ChenおよびShuicheng Yanによる「Network in network」、アーカイヴ・プレプリント、アーカイヴ1312.4400(2013年)を参照)。このステップでは、空間情報が破棄され、各要素が特定の特徴マップの平均的な活性化に対応するベクトルをもたらす。例えば、畳み込みニューラルネットワークの特定の隠れ層の場合、画像xの活性化は、次元wxhxcを有する可能性のある特徴マップAによって記述され、ここで、hは特徴マップのピクセル高さ、wは特徴マップのピクセル幅、および、cは現在の隠れ層内の特徴マップの数である。次に、Aの特徴マップのグローバル平均プーリング
Figure 0007223839000001
は次の式で与えられる。
Figure 0007223839000002
または
Figure 0007223839000003
例えば:
Figure 0007223839000004
次に、画像xのグローバル平均プール特徴マップAのベクトルは、次の式で与えられる。
Figure 0007223839000005
上記の式は、他の種類の観測値にも適用し得る。
さらに、上記の態様のいずれか1つによる方法は、深層ニューラルネットワークのトレーニング、またはトレーニングデータセットを使用してトレーニングされた深層ニューラルネットワークを微調整することを含み得る。
上記の態様による方法が、外れ値と非外れ値をどれだけうまく区別できるかは、例えば、受信者操作特性・曲線下面積(ROC-AUC)指標を使用して評価することができる。例えば、受信した観測値が外れ値であるかどうかを判断するために、マハラノビス距離および/または二乗近似残差および/または他の距離指標が計算される場合、マハラノビス距離および/または二乗近似残差および/または他の距離指標が、外れ値を非外れ値からどれだけうまく分離するかを比較することにより、ROC曲線を計算し得る。本方法の評価に使用できるその他の指標には、(1)適合率、これは、実際には外れ値であるシステムにより外れ値として識別された観測値の一部、(2)再現率、これは、システムによって外れ値として正常に識別されたすべての外れ値の一部、(3)F1値、これは、適合率と再現率の調和平均、(4)正解率、これは、すべての識別のうち、システムによる正しい識別の一部、が含まれるが、これらに限定されない。
さらに別の態様によれば、コンピュータプログラム製品が提供される。コンピュータプログラム製品は、コンピュータにロードされて実行されると、コンピュータに上記の態様のいずれかに従って方法を実行させるコンピュータ可読命令を含む。
さらに別の態様によれば、異常検出および/または予知保全のためのデータ分析のためのシステムが提供される。システムは、
対象のエンティティ(観測されるエンティティ)のパラメータを特徴付ける少なくとも1つの観測値を取得するためのデータ収集ユニットと、
深層ニューラルネットワークをトレーニングするために使用されるトレーニングデータセットを格納する記憶媒体であって、深層ニューラルネットワークは、複数の隠れ層を含み、トレーニングデータセットは、深層ニューラルネットワークに入力する可能性のある観測値を含む、記憶媒体と、
上記の態様および例のいずれか1つに従って方法を実行するように構成された少なくとも1つのプロセッサを含む予測分析ユニット(予測分析コンポーネント)と、
を含む。
例えば、プロセッサは、
トレーニングデータセットを使用してトレーニングされた深層ニューラルネットワークを取得し、
記憶媒体に格納されているトレーニングデータセットの少なくとも一部を取得し、
複数の隠れ層の少なくとも1つから出力される中間出力値の第1のセットを取得し、
ここで、中間出力値の第1のセットのそれぞれは、トレーニングデータセットの少なくとも一部に含まれる可能性のある観測値の異なる1つを入力することによって取得され、
中間出力値の第1のセットを使用して潜在変数モデルを構築/適合し、
ここで、潜在変数モデルは、中間出力のセットの次元よりも低い次元を有する部分空間内の射影値の第1のセットへの中間出力値の第1のセットのマッピングを提供し、
深層ニューラルネットワークに入力される新しい観測値を受診し、
受信した新しい観測値を深層ニューラルネットワークに入力することにより、複数の隠れ層のうちの前記少なくとも1つから出力される中間出力値の第2のセットを取得し、
潜在変数モデルを使用して、中間出力値の第2のセットを射影値の第2のセットにマップし、
潜在変数モデルおよび射影値の第2のセットに基づいて、受信した新しい観測値がトレーニングデータセットに関して外れ値であるかどうかを判断するように、
構成され得る。
少なくとも1つのプロセッサは、
深層ニューラルネットワークによって、新しい観測値の予測を計算し、
予測、および新しい観測値が外れ値であるかどうかの判断に基づいて、エンティティに少なくとも1つの異常が発生したことを示す結果を判断するように、
さらに構成され得る。
上記の態様によるシステムでは、潜在変数モデルは、主成分分析に従って、またはオートエンコーダを使用して、構築され得る。
上記の態様によるシステムにおいて、受信された観測値が外れ値であるかどうかを判断するために、プロセッサは、
射影値の第1のセットの分布に関して、射影値の第2のセットの距離を計算し、
計算された距離が距離の閾値よりも大きい場合、受信した新しい観測値がトレーニングデータセットに関して外れ値であると判断するように、
さらに構成され得る。
上記のように、距離は、マハラノビス距離、残差平方和、局所外れ値因子、複合距離、またはその他の適切な距離測定であり得る。
例えば、上記の態様によるシステムにおいて、受信された新しい観測値が外れ値であるかどうかを判断するために、少なくとも1つのプロセッサは、
潜在変数モデルと射影値の第2のセットを使用して、中間出力値の第2のセットに対応する中間出力値の近似セットを判断し、
中間出力値の第2のセットと中間出力値の近似セットの二乗近似残差を計算し、
計算された二乗近似残差が二乗近似残差の閾値よりも大きい場合、受信した新しい観測値がトレーニングデータセットに関して外れ値であると判断するように、
さらに構成され得る。
上記の態様によるシステムでは、少なくとも1つのプロセッサは、
中間出力値の第1のセットを取得し、複数の隠れ層のうちの2つ以上の潜在変数モデルを構築するステップを実行し、
中間出力値の第2のセットを取得し、中間出力値の第2のセットを、複数の隠れ層のうちの前記2つ以上に関する射影値の第2のセットにマッピングするステップを実行し、
潜在変数モデルと、複数の隠れ層のうちの前記2つ以上に関して取得された射影値の第2のセットに基づいて、受信した新しい観測値が外れ値であるかどうかを判断するステップを実行するように、
さらに構成され得る。
さらに、上記の態様によるシステムでは、少なくとも1つのプロセッサは、複数の隠れ層のうちの少なくとも1つからの活性化を判断し、任意選択で、判断された活性化にグローバルプーリングを適用するようにさらに構成され得る。したがって、(任意選択でプールされた)活性化は、第1/第2の中間出力値を構成し得る。
少なくとも1つのプロセッサは、汎用コンピュータ、専用コンピューティングユニット、サーバ、またはコンピュータネットワークの一部であり得る。
上記の態様によるシステムは、深層ニューラルネットワークのパラメータ(例えば、重み、ノード接続、フィルタなど)を格納するための記憶媒体をさらに含み得る。格納される特定のパラメータは、通常、使用する深層ニューラルネットワークの種類によって異なる。さらに、システムは、潜在変数モデル、距離、閾値距離などの潜在変数モデルを使用して得られたデータ、および/または他のデータを格納するための記憶媒体を備え得る。
上記の記憶媒体のうちの1つまたは複数は、収集されたデータが蓄積、処理、および/または分析されるデータストアの一部であり得る。データストアは、サーバ(例えば、オンプレミス・サーバやローカル・サーバなど)に実装することも、クラウドサービスなどの分散ネットワークサービスとして実装し得る。
さらに、システムは、状態を測定するための1つまたは複数のセンサ、または対象のエンティティまたは、対象のエンティティのコンポーネントの状態を特徴付ける少なくとも1つのパラメータを備え得る。測定されたパラメータは、温度、圧力、電流、磁束、ネットワークトラフィック、データログなどの物理的、電子的、またはその他のパラメータであり得る。少なくとも1つのセンサは、ソフトセンサであり得る。少なくとも1つのセンサは、定期的または実質的に連続的な測定を実行するように構成され得る。一例では、少なくとも1つのセンサは、リアルタイム測定を実行するように構成され得る。トリガ時に測定を開始するようにセンサを構成することもできる。例えば、測定は、ユーザとコンピュータとの相互作用などの特定のアクションの結果として開始され得る。一例では、複数のセンサが提供され得、深層ニューラルネットワークに入力された観測値は、異なる入力情報の組合せまたは融合、例えば、異なるセンサからの測定値に基づき得る。
さらに、システムは、1つまたは複数のセンサデータ収集ユニット、転送されたデータを格納するための記憶媒体、および/または予測分析ユニットの間のデータ転送を可能にするデータ通信ユニットまたは複数のラインを備え得る。一例では、1つまたは複数のセンサとデータ収集ユニットとの間のデータ転送は、無線であり得る。
システムは、予測、および/または少なくとも1つの機器および/または予測に基づくシステム全体の状態または状況に関する情報を伝達するように構成されたユーザインターフェース(例えば、グラフィカルおよび/またはオーディオおよび/または触覚ユーザインターフェース)をさらに含み得る。ユーザインターフェースは、少なくとも1つの機器の測定されたパラメータに関する情報を伝達する(例えば、表示する)ようにさらに構成され得る。ユーザインターフェースは、観測されるシステムの少なくとも1つの機器の少なくとも1つのパラメータを、ユーザが変更することができるように、さらに構成され得る。
予測分析ユニットは、例えば、いつ保守が必要とされるかを予測するために、収集されたデータを集約およびモデル化するようにさらに構成され得る。分析コンポーネントは、例えば、対象のエンティティ(例えば、機器、機械、コンピュータなど)で直接実行され、いわゆるエッジ分析を構成し得る。
上記の様々な様態および例によれば、深層ニューラルネットワークへのモデルトレーニング後に外れ値を検出することが可能である。深層ニューラルネットワークがどのように構築されているかについての仮定が少なくて済み、サンプリングが不要な場合があるため、トレーニングデータセットに関する入力データの外れ値を効率的に検出できる。さらに、上記の様々な態様および例によれば、すでに学習されたデータ表現を外れ値検出に使用することができる。個別の外れ値検出モデルとは対照的に、予測に使用されるのと同じ表現を外れ値の検出に使用できる。
上記の様々な態様および例による方法およびシステムのさらなる利点には、以下が含まれ得る。
-アーキテクチャに依存しない。つまり、あらゆる種類のニューラルネットワークモデルで使用できる。
-トレーニングに依存しない。計算には、ニューラルネットワークがトレーニングされた後のトレーニングデータのシングル・パスが含まれる。つまり、通常のトレーニング手順への干渉は最小限に抑えられる。
-タスクに依存しない。つまり、目前のタスクについて想定する必要がない。提案されたアプローチは、分類、セグメンテーション、回帰、強化学習などの様々なタスクに適用できる。これらは計算の点で効率的である。つまり、深層ニューラルネットワークモデルの推論と比較して、追加の計算はほとんど必要ない。
-外れ値の検出において、最先端の方法よりも優れている。
提案された技術は、あらゆる種類の入力、例えば、センサ入力、ネットワークトラフィック、金融取引トラフィックなどに適用可能であるという意味で、「汎用」技術である。上記の様々な態様および例による方法およびシステムの適用分野には、機械部品、機械、または複数の機械部品および/または機械を含むシステムの監視および/または予知保全が含まれるが、これらに限定されない。
特定の非限定的な例は、ベアリング、ギア、ホイール、制御バルブ、バッテリ、モータ(例えば、電気モータ)、タービン、車両(例えば、鉄道車両、電気自動車)、バイオリアクタ、発酵槽、コンピュータ、コンピュータネットワーク、レール、石油およびガスの地表機器、産業または発電所(太陽光発電所、石炭燃焼プラントなど)の制御および/または監視、産業プロセスの制御および/または監視(例えば、ダイキャスティング、バッテリ製造、バイオテクノロジープロセス、医薬品製造など)の、異常(例えば、故障)検出および/または予知保全に関する。
別の可能な用途は、サイバーセキュリティの脅威、例えば、ネットワーク侵入やその他の種類の悪意のある攻撃、を検出するためのものである。例えば、ネットワークトラフィックは、機械学習技術を使用してモデル化し得る。例として、悪意のあるアクションおよび/またはアクタを示すトラフィックパターンの微妙な変化を検出し、侵入として分類し得る。さらに別の可能な用途は、金融取引における不正、疑わしいパターン、または特定の傾向を検出するためのものである。例えば、履歴データを使用して、機械学習に基づく分類器をトレーニングして、トランザクションまたはアクタ(例えば、不正または信頼できる)を分類したり、傾向を識別したりできる。
本出願に記載されている主題は、方法として、またはシステムとして、おそらく1つまたは複数のコンピュータプログラム製品の形で実装することができる。本出願に記載されている主題は、データ信号または機械可読媒体に実装することができ、媒体は、CD-ROM、DVD-ROM、半導体メモリ、ハードディスクなどの、1つまたは複数の情報キャリアに組み込まれている。そのようなコンピュータプログラム製品は、データ処理装置に、本出願に記載されている1つまたは複数の操作を実行させることができる。
さらに、本出願に記載されている主題は、プロセッサ、およびプロセッサに結合されたメモリを含むシステムとして実装することもできる。メモリは、1つまたは複数のプログラムをエンコードして、プロセッサに、本出願で説明されている1つまたは複数の方法を実行させることができる。いくつかの例では、システムは、汎用コンピュータシステムであり得る。他の例では、システムは、組み込みシステムを含む専用コンピュータシステムであり得る。
1つまたは複数の実装の詳細は、以下の例示的な図面および説明に示されている。その他の特徴は、説明、図面、および特許請求の範囲から明らかになるであろう。しかしながら、実施形態が別々に説明されているとしても、異なる実施形態の単一の特徴を組み合わせてさらなる実施形態にすることができることを理解されたい。
本開示によるシステムを含むシステムの一例を示す。 画像を処理するための畳み込みニューラルネットワーク(CNN)の構成例を示す。 本開示によるシステムによって実行される例示的なプロセスのフローチャートを示す。 予知保全モデルに使用されるネットワークアーキテクチャの例、および、外れ値の検出のために抽出された特徴の例を示す。 予知保全のためのテストセットシーケンスでの外れ値の検出を示す。 PCA成分の数、およびカットオフパーセンタイルに応じた予知保全のための外れ値検出性能の評価を示す。 難読化されたパワーシェルスクリプトを検出するための予知保全モデルに使用されるネットワークアーキテクチャの例を示す。 本発明による外れ値距離指標のバイオリン図(図8A)、および、MCドロップアウトを使用した比較例によるバイオリン図(図8B)を示す。 本明細書で説明されるシステムの少なくとも一部を実装するために使用され得るコンピュータの例示的なハードウェア構成を示す。
以下の文章では、図面を参照して例の詳細な説明を行う。例には、様々な変更を加えることができることを理解されたい。特に、1つの例の1つまたは複数の要素を、組み合わせ、他の例で使用し、新しい例を形成することができる。
システム構成
図1は、本開示によるシステムを含むシステムの一例を示す。
図1に示されるシステムは、コンピューティングシステム1、ユーザクライアント20、制御システム30、およびサードパーティシステム40を含み得る。
コンピューティングシステム1は、本開示によると、データ分析のためのシステムを提供し得る。コンピューティングシステム1は、例えば、1つまたは複数の汎用コンピュータを使用して実装し得る。図1に示されるように、コンピューティングシステム1は、アプリケーション10およびデータ記憶装置12を備え得る。
アプリケーション10は、後でより詳細に説明するように、コンピューティングシステムの例示的なプロセスをコンピュータに実行させる命令を含む、ソフトウェアアプリケーションによって実装し得る。図1に示されるように、アプリケーション10は、深層ニューラルネットワーク100、予測110、外れ値指標112、およびインターフェース114を含み得る。アプリケーションは、予測分析コンポーネントまたはユニットによって提供または実行され得る。
深層ニューラルネットワーク(以下、「DNN」とも呼ばれる)100は、入力層、出力層、および入力層と出力層との間に複数の隠れ層を有する人工ニューラルネットワークであり得る。深層ニューラルネットワーク100は、例えば、センサデータ、コンピュータネットワークトラフィックデータ、トランザクションデータなどのような任意の種類のデータを処理するように、トレーニングされ得る。一例では、深層ニューラルネットワーク100は、深層ニューラルネットワーク100への可能な入力画像を含むトレーニングデータセットを使用して、それぞれのセンサによって取得された画像を処理するようにトレーニングされ得る。トレーニングデータセットは、アプリケーション10によってアクセス可能なデータ記憶装置12に格納され得る。図1の例では、深層ニューラルネットワーク100は、レイヤ102および外れ値検出モジュール104を含み得る。
レイヤ102は、入力層(図示せず)、出力層(図示せず)、および入力層と出力層との間に提供される複数の隠れ層(図示せず)を含み得る。
外れ値検出モジュール104は、深層ニューラルネットワーク100の複数の隠れ層のうちの少なくとも1つに接続され得、深層ニューラルネットワーク100をトレーニングするために使用されるトレーニングデータセットに関して観測値が外れ値であるかどうかを判断するように構成され得る。外れ値検出モジュール104によって実行されるプロセスの詳細は、後で説明される。
図2は、深層ニューラルネットワーク100の特定の例を示す。図2に示される例示的な深層ニューラルネットワークは、データ、例えば画像を分類するための畳み込みニューラルネットワーク(CNN)である。例えば、図2に示される例示的なCNNは、入力画像などの入力観測値をいくつかの異なるカテゴリに分類し得る。別の例では、図2に示される例示的なCNNは、入力観測値(例えば、画像)が特定のカテゴリに属するかどうかを判断し得る。特定の例では、図2に示される例示的なCNNは、入力画像(すなわち、入力観測値)として28×28ピクセルのグレースケール画像を受信するように構成され得る。入力画像の各ピクセルは、ピクセルの強度を示す値を含み得る。入力画像の各ピクセルの強度値は、例示的なCNNの入力層の入力ノードへの入力値と見なすことができる。図2に示される例示的なCNNは、4つの畳み込み層C1、C2、C3、C4と、2つの最大プール層MP1、MP2と、出力層に含まれるノードの活性化関数としてソフトマックス関数を使用する出力層とを、含む。
上記の例では28x28ピクセルの画像を使用しているが、理論的な観点からは、画像サイズに上限はない。ただし、画像サイズを大きくすると、計算量も増え得る。一般に、対象のアプリケーション用に選択された画像サイズは、必要な詳細レベルと計算要求との間のトレードオフを反映し得る。一般的な分類アプリケーションでは、画像の範囲は200x200ピクセルから300x300ピクセルであり得る。ただし、アプリケーションはこの範囲に限定されず、200x200ピクセルより小さい、または300x300ピクセルより大きい、場合がある。
例示的なCNNの畳み込み層C1に含まれる各ノードは、入力画像の特定の部分に適用されるサイズ3×3(ピクセル)のフィルタに対応し得、フィルタを入力画像の特定の部分に適用することから生じる値を出力し得る。図2の例では、畳み込み層C1の入力画像に32個のフィルタが適用されている。32個のフィルタのそれぞれは、入力画像の幅方向と高さ方向の両方に、Sピクセルのストライドでフィルタをスライドさせることにより、入力画像の全領域に適用し得る。入力画像上の32個のフィルタの位置ごとに、32個のフィルタに対応する32個のノードが畳み込み層C1に存在し得る。図2の例では、ストライドSは1ピクセルに設定し得る。したがって、畳み込み層C1の出力は、28×28×32の値を含むことができ、これは、28×28ピクセルの32個の画像と見なし得る。畳み込み層C1の出力は、畳み込み層C2に入力され得る。畳み込み層C2は、サイズ3×3の32個のフィルタを備えた畳み込み層C1に類似した構成を持ち得る。したがって、畳み込み層C2の出力はまた、28×28×32の値を含むことができ、これは、28×28ピクセルの32個の画像と見なし得る。
畳み込み層C2の出力は、20%のドロップアウトの影響を受け得る(図2のD1を参照)。言い換えると、畳み込み層C2の出力の値(ノードに対応)の20%が、ランダムに選択され無視され得る。
ドロップアウト操作D1の後、畳み込み層C2の出力は、第1の最大プール層MP1での最大プーリング操作によるダウン・サンプリングの対象となり得る。最大プーリング操作は、複数の入力値の中から最大値を選択し得る。図2の例では、最大プール層MP1は、畳み込み層C2から出力された28×28ピクセルの32個の画像のそれぞれに、2×2のサイズのフィルタを適用し、2ピクセルのストライドで適用され得る。これにより、対応するフィルタ内の対応する画像(畳み込み層C2からの出力)のピクセルの中で最大強度値をそれぞれが有する、14×14ピクセルを含む32個の出力画像を生成し得る。最大プーリング操作で使用される各フィルタは、最大プール層MP1のノードと見なし得る。
最大プール層MP1からの出力は、畳み込み層C1、C2の構成に類似した構成を有するが64個のフィルタを適用する、畳み込み層C3に提供され得る。畳み込み層C3の出力は、サイズ3×3の64個のフィルタを備えた畳み込み層C3と同様の構成を有する畳み込み層C4に提供され得る。さらに、畳み込み層C4の出力は、20%のドロップアウトの影響を受け(図2のD2を参照)、上述のように、最大プール層MP1での最大プール操作に類似した方法で、最大プール層MP2で最大プーリング操作の対象になり得る。最大プール層MP2の出力には、7×7×64の値が含まれ得る。
最後に、最大プール層MP2の出力は、softmax関数を使用して出力層に提供され得る。出力層は、入力画像が分類され得る1つまたは複数のグループ(またはカテゴリ)に対応する1つまたは複数の出力ノードを含み得る。この例は特定のパラメータ(例えば、フィルタの数、ドロップアウト率、畳み込み層の数、ストライドなど)に言及しているが、各パラメータの値の範囲は、本明細書で意図され、方法およびシステムはこれらの実施形態に限定されない。
いくつかの例において、図2に示される例示的なCNNに関して、最大プール層MP1およびMP2からの出力は、図2に示す例示的なCNNのトレーニングに使用されるトレーニングデータセットに関して、入力画像が外れ値であるかどうかを検出するために、外れ値検出モジュール104(図1)に提供され得る。例えば、外れ値検出モジュール104は、最大プール層MP1およびMP2のそれぞれについて、トレーニングデータセット内の可能な入力画像のために、それぞれの最大プール層MP1およびMP2からの出力を使用する潜在変数モデルを、構築または適合し得る。潜在変数モデルは、最大プール層MP1またはMP2の出力を、最大プール層MP1またはMP2の出力の次元よりも低い次元を持つ部分空間(潜在変数モデルの部分空間)の射影値のセットにマッピングし得る。構築または適合された潜在変数モデルは、入力画像(例えば、トレーニングデータセットに含まれていない新しく取得された入力画像)が、トレーニングデータセットに関して外れ値であるかどうかを判断するために使用され得る。より一般的には、外れ値検出モジュール104は、深層ニューラルネットワーク100の1つ(または複数)の隠れ層から出力を取得し、取得された出力を使用して潜在変数モデルを構築または適合させるように構成され得る。潜在変数モデルの構築または適合の詳細については、以下で説明する。
再び図1を参照すると、予測110は、深層ニューラルネットワーク100からの出力であり得る。いくつかの例では、予測110は、入力画像が深層ニューラルネットワーク100によってどのカテゴリに分類されるかを示し得る。
外れ値指標112は、入力画像が深層ニューラルネットワーク100をトレーニングするために使用されるトレーニングデータセットに関して、外れ値であるかどうかを決定するための指標を含み得る。外れ値指標112の詳細は、以下で説明する。
インターフェース114は、アプリケーション10がコンピューティングシステム1の外部に提供され得る様々なデバイスと通信するためのインターフェースであり得る。例えば、インターフェース114は、アプリケーション10によって生成された情報をそれらのデバイスに通信するように構成され得る。さらに、例えば、インターフェース114は、それらのデバイスからアプリケーション10に向けられた情報を受信するように構成され得る。
データ記憶装置12は、アプリケーション100によって使用されるデータを格納するように構成され得る。図1は、データ記憶装置12がコンピューティングシステム1の一部であることを示しているが、いくつかの例では、データ記憶装置12に記憶されたデータがアプリケーション10によりアクセス可能である限り、データ記憶装置12は、コンピューティングシステムの外部で提供されてもよい。
ユーザクライアント20は、コンピューティングシステム1に接続されたクライアントデバイスであり得る。ユーザクライアント20は、コンピューティングシステム1で実行される外れ値検出の予測および結果を使用することができるユーザアプリケーション22を含み得る。ユーザクライアント20の特定の例は、例えば、SSH(セキュアシェル)またはHTTP(ハイパーテキスト転送プロトコル)要求を使用して、計算サーバにリモート接続されたワークステーションであり得る。次に、DNN(深層ニューラルネットワーク)を、計算サーバ上のユーザが提供する入力に適用することができ、結果として得られる予測および外れ値指標を、ユーザクライアント20に返すことができる。ユーザクライアント20は、例えば、DNN予測を実行するように構成されたワークステーション上で、アプリケーション10を実行するコンピューティングシステム1と同じ物理デバイスの一部であり得る。
制御および/または分析システム30は、デバイスを制御し、および/またはコンピューティングシステム1で実行される外れ値検出の予測および結果を使用して、さらなるデータ分析を実行し得る。制御および/または分析システム30は、異常検出および/または予知保全のためのシステムを構成するか、またはその一部であり得る。制御および/または分析システム30の例は、機械部品または機械、産業プロセスまたはプラント、車両(自律走行車など)、コンピュータネットワーク、金融取引ユニットなどの制御および/または分析システム(異常検出または予知保全など)であり得る。
例えば、車両において、制御および/または分析システムは、車両操縦を制御するため、および/または、車両に含まれる1つまたは複数の機器(モータ、ギアなど)の状態、摩耗、および/または寿命に関連する表示を提供するために、外れ値検出の予測および結果を受け取り、および/または必要な修正保守アクションを受け取り得る。
制御および/または分析システム30は、入力を受信し、データ分析を実行し、アプリケーション固有のアクションに関する決定を行う、制御および/または分析システム30で実行されるプログラムクライアント32を含み得る。例えば、それは、観測されるエンティティ(例えば、機器またはシステムの一部)の保守に関して、および/または、観測されるエンティティ(例えば、機器またはシステムの一部)の制御に関して行われる。
本技術の特定の例は、バイオ医薬品製造プロセスの分析および/または制御システムへの入力材料の特徴付けに関する。そのような製造プロセスの所望の出力品質を確保するために、プロセスの制御パラメータは、入力材料の特性に従って調整され得る。さらに、例えば、センサ・温度・湿度など、調整器、栄養溶液、細胞培養物などの洗浄または交換で、修正保守操作が、決定および/または開始され得る。
データソースの異種の組合せを使用することで、構造化データと非構造化データの両方を含む入力材料を特徴付けることができる。データソースの例は、イベントレポート、画像、クロマトグラム、分光データ、化学分析などのテキストを含み得る。深層学習は、構造化データと非構造化データの両方から特定の出力への有用なマッピングを見つけることができるため、このような異種のデータの組合せを分析するのに適している。入力データが深層学習モデルのトレーニングに使用されたデータと一致しない場合の出力品質の低下や致命的な障害を回避するには、モデルのトレーニング後に外れ値を検出することが重要になり得る。入力データが外れ値として識別された場合、出力品質が保証されないため、入力データを製造プロセスへの入力として使用できない場合がある。
本技術の別の特定の例は、人間のコントローラによってプログラムされるのではなく、深層強化学習を使用してデータから学習された制御および監視システム(予知保全システムなど)にある。バイオ医薬品の製造プロセスでは、可能な制御パラメータが多数存在する可能性があり、高度に複雑なため、システムに存在するすべての相互作用とフィードバック・ループを見つけるのが難しい場合がある。代わりに、深層強化学習を使用して、データとシミュレーションから制御ポリシを学習できる。プロセスがモデルトレーニング中に遭遇した状態とは異なる状態に入るタイミングを検出するには、外れ値の検出が重要になり得る。そのような状態が外れ値として識別された場合、深層学習に基づく制御システムは、ユーザに警告し(例えば、ユーザに通知を送信することによって)、および/または出力品質の低下や致命的な障害を回避するために安全ルーチンを実行し得る。深層強化学習の使用例は、当技術分野に見出すことができる(例えば、Li,Y.、Wen,Y.、Guan,K.およびTao,D.による(2017年)「Transforming Cooling Optimization for Green Data Center via Deep Reinforcement Learning」、アーカイヴ・プレプリント、アーカイヴ1709.05077、Nguyen,P.およびTakashi,E.による(2018年)「Automating Water Purification Plant Operations Using Deep Deterministic Policy Gradient」、ワークショップの講演論文集(エンジニアリングシステムにおける安全上重要な深層学習のためのICMLワークショップで公開))。
深層学習モデルにおける外れ値検出のための従来技術の方法とは異なり、本方法は、モデルがどのように設計および/またはトレーニングされるかを制約せず、これらの方法を稼働中のモデルに追加することを可能にする。さらに、本技術はモデルトレーニングの後に使用されるので、本技術は、ライブアプリケーションに実装された、リアルタイムアプリケーションおよび/またはモデルで、使用できる。異常検出のアプリケーションは非常に高速であり得る。つまり、予測モデルは、高周波数で観測値を受信する。例えば、回転タービンの監視は10kHzで測定し得る。同様に、コンピュータネットワークまたは銀行システムを監視するシステムは、通常、そのようなシステムには非常に多くの同時アクタが存在するため、高頻度で観測値を受信し得る。対照的に、予測ごとに複数のフォワードパス/推論を使用する従来方法(例えば、予測が平均として使用され、標準偏差またはエントロピーが不確実性の尺度として使用される場合)は、待ち時間の増加につながる可能性があり、リアルタイムアプリケーションにとっては、問題となる。さらに、提案された技術は、アーキテクチャに依存せず、トレーニングに依存せず、タスクにも依存しない。
サードパーティシステム40は、コンピュータによって実装され得、サードパーティアプリケーション42を含み得る。サードパーティシステム40は、コンピューティングシステム1が属するエンティティとは異なるエンティティに属し得る。サードパーティアプリケーション42は、コンピューティングシステム1で実行された外れ値検出の予測および結果を使用し得る。サードパーティアプリケーション42の特定の例は、DNN予測およびDNN予測プロバイダから購入した外れ値検出の結果に依存するソフトウェアであり得る。
コンピューティングシステム1は、インターネットなどのネットワークを介して、ユーザクライアント20、制御システム30、およびサードパーティシステム40に接続され得る。いくつかの例では、ユーザクライアント20、制御システム30、およびサードパーティシステム40のうちの1つまたは複数は、イントラネットを介してコンピューティングシステム1に接続され得る。
部分最小二乗回帰を使用した予測-時間外れ値の検出
上記のように、図1に示されるコンピューティングシステム1の外れ値検出モジュール104は、深層ニューラルネットワーク100のトレーニングに使用されるトレーニングデータセットに関して、入力観測値(例えば、入力画像)が外れ値であるかどうかを判断するために、深層ニューラルネットワーク100(例えば、図2の最大プール層MP1またはMP2)の複数の隠れ層の少なくとも1つから出力される中間出力値を使用して、潜在変数モデルを構築または適合させるように構成され得る。中間出力値は、例えば、深層ニューラルネットワーク100の複数の隠れ層のうちの少なくとも1つからのそれぞれの入力画像の活性化であり得る。任意選択で、活性化が特徴マップからのものである場合、グローバルプーリング(例えば、グローバル平均プーリング)を各特徴マップ全体に適用し得る。
以下は、部分最小二乗回帰を使用する外れ値検出のための例示的な技術を提供し、これは、外れ値検出モジュール104がどのように外れ値を検出できるかについての理解を容易にし得る。
モデルの不確実性の概念は、予測時間の外れ値の検出の概念に関する。どちらの場合も、主な目的は、モデルがトレーニングされた観測値に対応しない観測値を見つけることであり得る。予測モデルは、製造プロセス監視における予測中に、分布外の観測値を検出するために長い間使用されてきた。この分野では、部分最小二乗(PLS)回帰などの線形潜在変数モデルが広く使用されている(例えば、P.GeladiおよびB.R.Kowalskiによる「Partial least-squares regression:a tutorial」、アナリティカ・ケミカ・アクタ、第185号、補足C、第1~17頁、1986年1月)。PLSモデルは、主成分分析(PCA)に類似するトレーニングデータの部分空間近似で見つけることができ(S.Wold、K.Esbensen、およびP.Geladiによる、「Principal component analysis」、ケモメトリクスと知的実験システム、第2巻、第1号、第37~52頁、1987年8月)、これは予測に使用される。次に、PLSモデルにより見つけられた部分空間に、新しい観測値を射影できる。したがって、部分空間内の距離と部分空間からの距離の両方を使用して、予測時間中に外れ値を見つけることができる。
予測時間とは、モデルが適合され、モデルがアプリケーションでの予測に使用された後の時間を指す(モデルが生成されるときのモデルの適合/トレーニングの時間は含まれない)。
単純で効果的ではあるが、上記のPLSモデルの原則を、深層学習システムに直接適用できない場合があることに注意すべきである。具体的には、PLSの予測の基礎は入力データ空間の部分空間であるため、新しい観測値の部分空間近似をトレーニングに使用されたものと簡単に比較できる。一方、深層学習モデルは通常、非線形変換の複数のレイヤに依存し、各変換が最終出力に寄与する。これは、表面上は異なる観測値が、空間のシーケンスを介して変換され、出力にマップされることを意味することがある。入力から出力への有用な変換を見つけるこの機能は、深層学習モデルが非構造化データに対して、うまく機能する理由であり得る。ただし、これらの変換の結果、この変換が、新しい観測値のモデルがトレーニングされたのと同じデータ分布に属しているかどうかを判断することを困難にする場合がある。深層学習モデルの知識の限界を判断するための明確な比較ポイントがない場合や、新しい観測値がそれらの限界内にあるかどうかを判断するための明確な方法がない場合がある。
PLS回帰は、製造プロセスの監視などで広く使用されている回帰モデルである。観測値のn行と特徴のp列を含むトレーニングデータ行列(X=[x...xで表記される)、応答のq列を含む対応するターゲット行列(Yで表記される)が与えられると、PLSモデルは両方のデータ行列を要約し得る。PCAと同様に、PLSは、XおよびYを近似する「潜在変数空間」と呼ばれる、等しい次元kのフルランクの部分空間を見つけ得る。PCAとは対照的に、PLSは、単一の部分空間の分散を最大化する代わりに、検出された部分空間間の共分散を最大化し得る。つまり、PLSは入力データを次のように近似し得る。
X=TP+E (1)
ここで、T=[t,... ,t]は、Xの列・空間のフルランク部分空間にまたがりYとの共分散を最大化するn×kの潜在変数行列であり、Pは、PLS負荷のp×k行列であり、Eは、近似残差のn×p行列であり得る。負荷には、TからXを近似するために使用される線形結合の重みで含み得る。Yの応答行列は、同様の方法で近似し得る。PLSモデルを計算して回帰に使用する方法の詳細については、以下で見つけられる。S.Wold、M.SjostromおよびL.Erikssonらによる「PLS-regression:a basic tool of chemometrics」、Chemom.Intell.Lab.Syst.、第58巻、第2号、第109~130頁、2001年10月。
PLSによって提供される入力データの近似値を使用して、予測中に遭遇した外れ値を検出できる。モデルの適合後に、次の方法でPLSを使用して検出された部分空間に新しい観測値を射影し得る。
new=XnewP (2)
これにより、新しい観測値を次のように近似できる場合がある。
Figure 0007223839000006
部分空間射影Tnewおよび近似の両方を使用して、適切な距離指標に基づいて、外れ値を検出し得る。予測時間中に外れ値を検出するために使用される2つの例示的な尺度は、潜在変数空間のマハラノビス距離と入力列空間の近似残差である。他の尺度は、複合距離、局所アウトライン因子(LOF)などである。
PLSを使用する場合、予測時間中に外れ値を検出する1つの方法は、マハラノビス距離を使用することである。マハラノビス距離は、多変量確率分布の原点から観測値までの標準偏差が、どの程度離れているかを測定する既知の統計距離である。直感的には、マハラノビス距離が大きい観測値は、与えられた確率分布の下では、可能性が低いということであり得る。潜在変数の確率分布は、潜在変数の平均がゼロであると仮定して、Tにおける潜在変数の共分散行列Cを使用して、パラメータ化し得る。次に、潜在変数空間の射影値tを備える観測値xのマハラノビス距離dは、次のように計算し得る。
Figure 0007223839000007
マハラノビス距離は、観測値が特定の分布の下にある可能性がどれほど低いかを測定できるが、マハラノビス距離は、観測値がその分布に属しているかどうかについての情報を提供しない場合がある。対照的に、近似残差は、新しい観測値が分布からどれだけ離れているかを簡単に測定できる。適合されたPLSモデルが与えられると、観測値xは、
Figure 0007223839000008
として近似し得る。二乗近似残差は、次の式で簡単に求め得る。
Figure 0007223839000009
ここで、xi,jは、観測値ベクトルのj番目の要素であり得る。直感的には、トレーニング分布以外の別の分布からの観測値の近似が失敗し、近似残差が大きくなるということであり得る。
深層ニューラルネットワークの予測時間外れ値検出
外れ値検出モジュール104は、深層ニューラルネットワーク100が説明できない可能性がある観測値を見つけるために、例えば、上記のように製造プロセス監視における予測時間外れ値の原理に基づいて、アドオン方法を実行し得る。深層ニューラルネットワーク100は、任意の種類のニューラルネットワーク、例えば、畳み込みニューラルネットワーク(CNN)、リカレントニューラルネットワーク、リカーシブニューラルネットワーク、トランスフォーマニューラルネットワークなどであり得る。
深層学習システムにおいて予測時間の外れ値を検出するために外れ値検出モジュール104によって実行される方法は、ニューラルネットワークが入力データを変換することによって機能し得るという事実に基づき得る。入力データが深層ニューラルネットワークを介して供給される場合、データの複数の中間表現が存在する可能性があり、(例えば、深層ニューラルネットワークが分類問題を解決するように構成されている場合に、入力データが分類されるグループの)予測に、中間表現が使用され得る。これらの中間表現の1つまたは複数は、予測中に外れ値を検出するため、および予測を実行するため、に使用され得る。
PLSと同様の方法で、予測中に外れ値を検出する可能性を提供するために、深層ニューラルネットワーク(例えば、フィードフォワードニューラルネットワーク)を、一連の非線形変換と見なすことができる。言い換えると、特定のネットワークのレイヤi(=1、2、3、4、....)からの観測値x(入力データ)の活性化ベクトルaは、ネストされた一連の変換によって次のように与えられる。
Figure 0007223839000010
ここで、f(k=1、2、...、i)は活性化関数であり得、W(k=1、2、...、i)は重み行列であり得る。観測値xの活性化ベクトルaは、所与のネットワークのレイヤiからの中間出力と見なすことができ、観測値xが所与のネットワークに入力されるときのレイヤiのそれぞれのノードからの出力に対応する要素値を含み得る。これらの活性化aのそれぞれは、入力データの特徴表現を提供し得る。重み行列は、通常、バックプロパゲーションによる教師ありトレーニングによって取得し得るが、活性化は、入力データの変換された、または前処理された表現を提供するだけであり得る。
変換をより詳細に説明するために、観測値xは、n次元の行ベクトルx=[x...x]であり得、ここで、nは、アプリケーションに依存する。xが単一チャネル画像の場合、nは、ピクセル単位の画像の長さにピクセル単位の画像の幅を掛けたものであり得、値x,x,...,xは画像のピクセル値である。第1のレイヤでは、xは、行列に第1の重み行列Wを掛けたものであり、n次元の線形射影
Figure 0007223839000011
を形成する。Wの次元は、n×nである。線形射影の後、活性化関数fが、
Figure 0007223839000012
に適用されて、活性化aが形成される。活性化関数fは非線形関数であり得る。活性化関数の一般的な選択肢には、とりわけ、正規化線形関数f(x)=max(0,x)、シグモイド関数f(x)=(1+e-x-1、ソフトマックス関数
Figure 0007223839000013
が含まれる。次に、活性化aは、行列に重み行列Wを掛けたものであり、結果の線形射影
Figure 0007223839000014
は、活性化関数を使用して変換される。重み行列を用いた行列乗算と活性化関数を使用した変換手順は、レイヤiの活性化aが得られるまで、i回繰り返される。
これらの中間表現を使用して、任意の深層ニューラルネットワーク(例えば、任意のフィードフォワードニューラルネットワーク)の外れ値を検出するためのプラグ・アンド・プレイ方式を提供し得る。トレーニングされたネットワークが与えられた場合、トレーニングデータは、レイヤiからの活性化行列A=[a1,i... an,iを使用して表すことができ、ここで、nは観測値の数である。Aは、フルランクではない場合があるため、近似Aとして使用できるフルランクの部分空間を見つけることが望ましい場合がある。いくつかの例では、PCAを使用して、データを低次元の部分空間にマッピングするために、線形射影を使用して部分空間を取得し得る。他のいくつかの例では、オートエンコーダ、ランダム射影、スパース主成分分析、および/または自己組織化マップなどの他の方法を使用して、データを低次元の部分空間にマッピングし得る。
PCAの使用例では、以下を取得するために、トレーニングデータの活性化の次元をm次の次元に減らすことができる。
,P=PCA(A) (7)
PLSと同様に、Tは、トレーニングセットの活性化の部分空間にまたがる潜在変数行列を示し得、Pは、PCA負荷を示し得る。例えば、行列形式のPCAは、次の最小二乗モデルを提供し得る。
Figure 0007223839000015
これは、上記の式(1)に類似していると見なし得る。
PLSと同じように、潜在変数の共分散行列と負荷行列を使用して、マハラノビス距離、近似残差、複合計量距離、局所アウトライン因子(LOF)などの適切な距離指標に基づいて、予測時間の外れ値を検出し得る。
トレーニングされたニューラルネットワークを使用して新しい観測値が予測されると、活性化Ai,newも抽出し得る。新しい活性化は、PLSを使用する場合と同じ方法で、トレーニングデータの活性化から見つかった部分空間に射影し得る。
A,new=Ai,new (9)
そして、新しい観測値の距離、例えばマハラノビス距離は、式(4)に従って計算し得る。
新しい活性化は、PCA負荷を使用して、次のように近似し得る。
Figure 0007223839000016
次の式(11)(式(5)に類似)に従って計算された近似残差を使用しても、PLSを使用する場合と同じ方法で、外れ値を検出し得る。
Figure 0007223839000017
ここで、pは、レイヤi内のノードの数を示し得る。
決定された距離はまた、記述された距離の2つ以上を組み合わせることによって形成される複合距離指標であり得る。1つの例は、ユークリッドノルムを次のように使用して、残差平方和の平方根によって与えられるマハラノビス距離および残差距離を組み合わせる方法である。
Figure 0007223839000018
上記は複合距離指標の一例にすぎず、他の多くの組合せも可能である。
上記の式、すなわち式(7)から(12)は、深層ニューラルネットワークからの活性化を適用するために、本明細書で提供されるような新しい方法で適合されている。
システム設定のプロセス
図3の左側の部分は、コンピューティングシステム1を設定するための例示的なプロセスのフローチャートを示す。図3に示される例示的なプロセスは、コンピューティングシステム1のアプリケーション10によって実行され得る。
ステップS10で、アプリケーション10は、データ記憶装置12からデータを入手し得る。データは、深層ニューラルネットワーク100をトレーニングするためのトレーニングデータセットの少なくとも一部であり得る。トレーニングデータセットは、深層ニューラルネットワーク100への可能な観測値入力を含み得る。トレーニングセットを構成する観測値は、センサデータ、コンピュータネットワークトラフィックデータ、金融取引データ、または観測された機器やシステムを特徴付けるその他の履歴データであり得る。例えば、深層ニューラルネットワーク100が図2に示されるようなCNNである場合、トレーニングデータセットは、可能性のある観測値として、CNNへの可能な入力画像を含み得る。
ステップS12で、アプリケーション10は、深層ニューラルネットワーク100がすでにトレーニングされているかどうかを判断し得る。ステップS12は、ステップS10の前、後、または並行して実行し得る。
深層ニューラルネットワーク100がまだトレーニングされていないと判断された場合(ステップS12で、いいえ)、プロセスはステップS14に進み得る。深層ニューラルネットワーク100がすでにトレーニングされていると決定された場合(ステップS12ではい)、プロセスはステップS16に進み得る。あるいは、ステップS16に進む前に、深層ニューラルネットワーク100を微調整し得る。
ステップS14で、アプリケーション10は、ステップS10で入手されたデータを使用して、深層ニューラルネットワーク100をトレーニングし得る。
ステップS16で、アプリケーション10は、深層ニューラルネットワーク100を使用してデータを変換し得る。ステップS17で、アプリケーション10は、深層ニューラルネットワーク100の複数の隠れ層のうちの少なくとも1つから出力される中間出力値(中間出力値の第1のセット)を取得し得る。具体的には、例えば、アプリケーション10は、トレーニングデータセット内の可能な入力画像の、上記の式(6)に従って、深層ニューラルネットワーク100の少なくとも1つのレイヤから活性化ベクトルを取得し得る。活性化とは、式(6)で説明されているように、観測値による活性化関数を使用した行列乗算と変換のシーケンスの結果を指す。任意選択で、活性化が特徴マップからのものである場合、各特徴マップ全体のグローバルプーリング、例えば平均プーリングを適用し得る。
ステップS18で、アプリケーション10は、例えば(任意選択でプールされた)トレーニングセットの活性化において、中間出力値に潜在変数モデルを適合させ得る。言い換えれば、アプリケーション10は、潜在変数モデルを構築し得る。例えば、アプリケーション10は、上記のようにPCAを使用して、潜在変数行列TおよびPCA負荷Pを取得し得る(例えば、式(7)および(8)を参照)。あるいは、アプリケーションは、(任意選択でプールされた)活性化ベクトルを入力として使用してオートエンコーダをトレーニングし得る。
ステップS20で、アプリケーション10は、潜在変数モデルを使用して活性化を射影し得る。例えば、アプリケーション10は、ステップS17で得られた活性化ベクトルから、潜在変数モデルを構築/適合させることによって見出された部分空間における射影値の対応するセット(例えば、T)を取得し得る。
ステップS22で、アプリケーション10は、潜在変数モデルを使用して、トレーニングセット画像の1つまたは複数の距離のセットを計算し得る。例えば、アプリケーションは、データの活性化の距離を計算し得る。任意の適切な距離指標を使用し得る。例えば、アプリケーション10は、上記の式(4)に従って、ステップS18で構築された潜在変数モデルに関して、ステップS17で取得された活性化ベクトルのそれぞれについてマハラノビス距離を計算し得る。追加的または代替的に、例えば、アプリケーション10は、上記の式(11)に従って、ステップS17で取得された各活性化ベクトルの二乗近似残差を計算し得る。
ステップS24で、アプリケーション10は、距離の閾値を決定し得る。追加的または代替的に、二乗近似残差の閾値を決定し得る。閾値は、後で、新しい観測値(例えば、入力画像)がトレーニングデータセットに関して外れ値であるかどうかを判断するために使用し得る。閾値を取得するために、ステップS22で計算された距離および/または二乗近似残差を使用し得る。例えば、閾値は、ステップS22で計算された距離(または二乗近似残差)のパーセンタイル、例えば95パーセンタイルで、あり得る。「95パーセンタイル」は単なる例であり、95より大きいまたは小さい値を、閾値を決定するためのパーセンタイルとして使用し得ることに注意すべきである。
システム設定のプロセスは、ステップS24の後に終了し得る。
外れ値検出のプロセス
図3の右側は、外れ値を検出するためにコンピューティングシステム1によって実行される例示的なプロセスを示している。
このプロセスは、深層ニューラルネットワークを導入することから始まり得る。例えば、深層ニューラルネットワークは、観測値を処理するために使用されるソフトウェアに組み込み得る。これは、監視対象の機器またはシステム自体の内部にある場合もあれば、例えば、サーバまたはクラウドサービスとして実現される分析用の独立型ソフトウェアモジュールとしての場合もある。
ステップS30で、アプリケーション10は、新しい観測値を受信し得る。例えば、深層ニューラルネットワーク100が図2に示すようにCNNである場合、CNNに入力される画像を、新たな観測値として受信し得る。ただし、観測値は画像に限定されず、観測される機器またはシステムを特徴付ける任意のデータであり得る。例えば、観測値は、1つまたは複数のセンサによって取得された特定の測定データであり得る。
ステップS32で、アプリケーション10は、深層ニューラルネットワークを使用して新しい観測値を変換し得る。変換には、例えば、新しい観測値の深層ニューラルネットワーク予測を計算することを含み得る。
その後、アプリケーションは、新しい観測値が外れ値であるかどうかを判断し得る。これは、次の方法で実行し得る。
ステップS33で、アプリケーション10は、深層ニューラルネットワーク100の複数の隠れ層のうちの少なくとも1つから出力される中間出力値(中間出力値の第2のセット)を取得し得る。例えば、アプリケーション10は、新しい観測値の、上記の式(6)に従って、深層ニューラルネットワーク100の少なくとも1つのレイヤから活性化ベクトルを取得し得る。
ステップS34で、アプリケーション10は、上記のようにステップS18で構築/適合された潜在変数モデルを使用して、決定された中間値(中間値の第2のセット)を射影し得る。例えば、アプリケーション10は、ステップS32で得られた活性化ベクトルから、潜在変数モデルを構築/適合させることによって見出された部分空間における対応する射影値のセット(例えば、TA,new)を取得し得る(式(9)を参照)。
ステップS36で、アプリケーション10は、上記のようにステップS18の潜在変数モデルを使用して、新しい観測値の距離を計算し得る。例えば、アプリケーション10は、上記の式(4)に従って、ステップS18で構築された潜在変数モデルに関して、ステップS32で取得された新しい観測値の活性化ベクトルのマハラノビス距離を計算し得る。追加的または代替的に、例えば、アプリケーション10は、上記の式(11)に従って、ステップS32で取得された新しい観測値の活性化ベクトルの二乗近似残差を計算し得る。上記のように、マハラノビス距離と二乗近似残差の代わりに、他の適切な距離指標を使用し得る。
ステップS38において、アプリケーション10は、ステップS36で計算された距離が、ステップS24で決定された閾値よりも大きいかどうかを判断し得る。
ステップS38で「はい」の場合、プロセスはステップS40に進むことができ、アプリケーション10は、新しい観測値が外れ値であると判断し得る。この場合、新しい観測値が外れ値であると判断されるため、システムはモデル予測を信頼できない予測として報告し得る。したがって、異常検出結果には、予測が新しい種類の異常に関連している、および/または信頼性の低い予測である、というデータを含み得る。プロセスは、ステップS40の後に終了し得る。あるいは、他のアプリケーション固有のアクションを実行し得る。例えば、深層ニューラルネットワークの予測は無視され、外れ値の観測値は破棄され得る。新しい種類の異常が検出されたという警告を生成し、任意選択で、外れ値の観測値および/または予測を保存して、さらに分析することもできる。
ステップS38でいいえの場合、プロセスはステップS42に進むことができ、アプリケーション10は、新しい観測値のために深層ニューラルネットワーク100によって行われた予測が信頼できると判断し得る。この場合、新しい観測値は外れ値ではないと判断されるため、システムはモデル予測を信頼できる予測として報告し得る。次に、予測は異常検出結果として出力し得る。プロセスは、ステップS42の後に終了し得る。あるいは、さらにアプリケーション固有のアクションを実行し得る。
予知保全の場合、アプリケーション固有のアクションは、有利には、実質的なエラーが実際に発生する前に、実行する必要のある修正保守アクションの種類を判断することに関し得る。判断された修正保守アクションに関する情報は、適切なユーザインターフェース(例えば、表示装置)を使用して、ユーザに伝達され得る。判断された修正保守アクションを、自動的にスケジュールおよび/または開始することも可能である。
サイバーセキュリティアプリケーションの場合、さらなるアプリケーション固有のアクションは、セキュリティ違反の種類を判断すること、または、特定のユーザ、コンピュータ、サーバ、ネットワークなどが信頼できないと見なされることなどであり得る。それぞれの情報は、適切なユーザインターフェース(例えば、表示装置)を使用して、ユーザに伝達され得る。さらに、例えば、特定のユーザ、コンピュータ、サーバなどへのアクセスを拒否するなどの修正処置は、自動的にスケジュールおよび/または開始し得る。
銀行業務において、さらなるアプリケーション固有のアクションは、電子取引、または、特定の参加者との電子取引が不正であると判断することに関連し得、さらに、適切なユーザインターフェース(例えば、表示装置)を使用して判断された情報をユーザに伝達することを含み得る。さらに、例えば、特定のユーザ、コンピュータ、サーバなどへのアクセスを拒否するなどの修正処置は、自動的にスケジュールおよび/または開始し得る。
本出願の実施形態および例では、深層ニューラルネットワークを使用する深層学習モデルが、それが何であれ、観測されるシステムを特徴付ける履歴データを使用してトレーニングされる。その後、深層学習モデルを使用して入力データ(新しい観測値)を処理し、正常な動作と既知の種類の異常を区別する。既知のアプローチとは対照的に、深層ニューラルネットワークを採用した深層学習モデルは、未知の種類の異常を識別することもできる。
実験
このセクションでは、例示的なコンピューティングシステム1によって実行された外れ値検出のための実験の結果が説明される。以下の実験では、深層ニューラルネットワーク100は、例えば、異常検出、特に予知保全のために、画像分類タスクを解決するように指示されている。
実施例1:ターボファンエンジンの予知保全
いくつかの実装形態では、CNNは、画像分析のための深層ニューラルネットワーク100として使用することができる。しかしながら、例示的なコンピューティングシステム1は、画像以外の種類のデータを分析する際に、外れ値検出を実行することもできる。例えば、例示的なコンピューティングシステム1による外れ値検出は、保守がいつ実行されるべきかを予測するために、稼働中の機器の状況を判断する予知保全タスクに対して実行され得る。このセクションでは、ターボファンエンジンの予知保全のための外れ値検出の実験結果を提供する。
この実験では、NASAエイムズのPrognostics Center of Excellenceによって提供されたターボファンエンジン劣化データセット(Saxena,A.、Goebel,K.による「Turbofan Engine Degradation Simulation Data Set」、NASAエイムズ予後診断データ・リポジトリ(http://ti.arc.nasa.gov/project/prognostic-data-repository)、2008年、NASAエイムズ研究センタ、モフェット・フィールド、カリフォルニア州)が使用された。データセットは、様々な動作条件と障害モードで、障害が発生するまで実行される模擬的なターボファンエンジンで構成されている。エンジンは、21個のセンサを使用して、経時的に監視された。使用されたセンサには、温度センサ、圧力センサ、ファンおよびコア速度センサ、冷却液漏れの測定値、および燃焼器の燃料・空気の比率の測定値が含まれていた。データセットの詳細は、例えば、Saxena,A.、Goebel,K.、Simon,D.およびEklund,N.による、「Damage propagation modeling for aircraft engine run-to-failure simulation」、予後および健全性管理、2008年、PHM2008、IEEE国際会議、第1~9頁、2008年10月で、見つけることができる。さらに、トレーニングデータセットの100個のエンジン、およびテストデータセットの100個のエンジンに対しての、3つの制御設定が記録されていた。課題は、致命的な障害を回避するために、障害が発生する15センサ・サイクル前に障害を予測することである。
障害を予測するために、長・短期記憶(LSTM)ブロック(Hochreiter、SeppおよびJurgen Schmidhuberによる「Long Short-Term Memory」、ニューロコンピューティング9、第8号、1997年、第1735~1780頁、Gers、Felix A.、Jurgen SchmidhuberおよびFred Cumminsによる「Learning to Forget:Continual Prediction with LSTM」、ニューロコンピューティング12、1999年、第2451~2471頁)に基づいて、ニューラルネットワークモデルがトレーニングされた。LSTMブロックは、リカレントニューラルネットワークのユニットとして理解でき、セル、入力ゲート、出力ゲート、および忘却ゲートで構成され得る。セルは、任意の時間間隔で値を「記憶する」ことができ、例えば、内部の「メモリ」を実装することができる。入力ゲート、出力ゲート、および忘却ゲートのそれぞれは、活性化関数を使用して加重和の活性化を計算する、ニューラルネットワークのノードと見なし得る。入力、出力、および忘却ゲートはセルに接続でき、LSTMブロックの接続を通過する値の流れの調整器と見なし得る。この実験では、LSTMモデルをコンピューティングシステム1の深層ニューラルネットワーク100として使用した。
図4は、この実験でターボファンエンジン劣化データセットを使用して障害を予測するためにトレーニングされたLSTMモデル(深層ニューラルネットワーク100の例)の構成例を示す。図4に示すLSTMモデルは、3つのLSTM層L1、L2、L3と、活性化関数としてシグモイド関数を持つ出力層で構成されている。LSTM層L1とL2はそれぞれ64個のLSTMブロックで構成され、LSTM層L3は50個のLSTMブロックで構成される。LSTM層L1、L2、L3は、深層ニューラルネットワーク100の隠れ層と見なし得る。図4に示すLSTMモデルをトレーニングするために、50サイクルの長さのスライディングウィンドウが入力として使用され、バイナリ応答がスライディングウィンドウの最後のサイクルが障害から15サイクル以内にあるかどうかを示した。21個のセンサ出力、3個の制御設定、および最初からの現在のサイクルが、故障予測の変数として使用された。すべての変数は0から1の範囲にスケーリングされ、トレーニングデータのスケーリングパラメータは、テストデータのスケーリングに使用された。
トレーニングの進行状況を監視するために、トレーニングデータセットの10%が検証に使用された。LSTMモデルは、Adamオプティマイザ(Kingma、Diederik P.およびJimmy Ba.による、「Adam:A Method for Stochastic Optimization」、アーカイヴ1412.6980[Cs]、2014年12月22日、http://arxiv.org/abs/1412.6980で入手可能[アクセス日:2018年7月25日])を使用して、検証セットのバイナリ・クロス・エントロピーの増加が停止するまで(8エポック)、トレーニングされた。結果として得られたモデルは、スライディングウィンドウ分類で、テストセットF1値が94.1%を達成した。予知保全モデルのテストセットの性能を次の表1に要約する。
Figure 0007223839000019
センサが機能しなくなる予測時間の外れ値を模擬実験するために、テストセット・エンジンの半分が外れ値として機能するようにランダムに選択された。これらのエンジンのそれぞれについて、センサがランダムに選択され、その出力はフルタイムシーケンスの途中からゼロに設定された。これらの外れ値を検出するための潜在変数モデルを提供するために、予知保全モデルのLSTM層L1、L2、L3から、特徴(例えば、中間出力)が抽出された(図4を参照)。最初の2つのLSTM層L1、L2では、各LSTMノードの時間全体の平均を使用して、各スライディングウィンドウを要約した。最後のLSTM層L3では、タイムシーケンスの最後のサイクルが使用された。次に、3つのLSTM層L1、L2、L3からの出力が、連結されて単一行列となった。
外れ値検出モジュールは、PCAに基づいて設定され、100成分のPCAモデルがすべてのトレーニングデータのLSTMに適合した。近似残差平方和は、すべてのトレーニングセットのスライディングウィンドウについて計算され、トレーニング残差平方和の99.9パーセンタイルが外れ値のカットオフとして設定された。次に、外れ値と残りのテストセットLSTM特徴がPCAモデルに射影され、図5からわかるように、すべてのスライディングウィンドウの近似残差平方和が計算された。したがって、外れ値の検出および管理を処理するために、予知保全モデルに加えて、本明細書に記載の方法およびシステムを使用し得る。予知保全モデルは、トレーニングされたタスクでは正確であるが、このモデルは、センサの破損(エンジンの破損ではない)などの予期しないイベントを処理できない。
図5は、予知保全のためのテストセットシーケンスでの外れ値検出の結果を示す。図5の左側は、センサ破損時と整合している1つのランダムなセンサ破損のシーケンスを示し(センサが破損したときを示す垂直の破線)、図5の右側は、中間時点で整合しているセンサ破損がないテストセットを示す。図5の上段の図は、外れ値検出の残差平方和を示す。水平の破線は、トレーニングセットの距離に基づく99.9パーセンタイルを示す。図5の中段の図は、タービンが15サイクル以内に破損する場合のモデル予測を示す。ここで、1は、エンジンが故障すると予測されることを意味する。図5の下段の図は、タービンが15サイクル以内に破損するかどうかを示すラベルを示す。ここで、1は、エンジンが故障することを意味する。外れ値は通常の変動とうまく区別された(この実験では適合率=再現率=1)。センサ破損に続いて、LSTMモデルは、エンジンの多くが破損することを予測するが、残差平方和は同時に増加する。テストセットの場合、モデルはエンジンの故障を正しく予測する一方で、残差平方和は増加しない。これは、現在の方法を使用することにより、信頼すべきでない予測と信頼すべき予測との区別することが可能であることを示している。また、本手法がフィードフォワードネットワークだけでなくリカレントネットワークにも適用可能であることも示している。
さらに、成分の数と残差平方和が、結果にどのように影響するかを調査した。図6からわかるように、成分の数は1から150まで変化し、外れ値検出の適合率および再現率は、トレーニングセットの残差平方和の95、99、および99.9パーセンタイルを、カットオフとして使用して報告された。図6は、PCA成分の数とカットオフパーセンタイルに応じた、予知保全のための外れ値検出性能の評価を示す。興味深いことに、すべてのカットオフで適合率および再現率は1が達成されたが、カットオフがより厳密になると、そのために必要な成分の数は増加する。すべてのカットオフについて、高い適合率と再現率を達成するには、説明された変動が非常に小さい成分が必要である(36成分のみで99.99%の累積R2に到達)。これは、外れ値と通常のテストセットシーケンスとの間で異なる変動が、小さく微妙な差異であり、非常に小さなPCA成分を補足する必要があることを示している。すべてのカットオフで、成分の数が多くなりすぎると、適合率が低下して過剰適合が観測された。ただし、より厳密なカットオフでは、過剰適合が発生する前に、より多くの成分が許容される。この実験では、トレーニングLSTM特徴の変動の99.999%を説明するPCAモデルが、99パーセンタイルカットオフおよび99.9パーセンタイルカットオフの両方で、適合率および再現率が1であり外れ値を検出するための優れた基盤を提供する。
実施例2:サイバーセキュリティの分野での異常検出
予知保全の分野において、本発明の他の特に効果的な使用例は、サイバーセキュリティの分野における使用である。サイバーセキュリティで深層学習を使用する1つの方法は、難読化されたパワーシェルスクリプト(Hendler D、Kels S、Rubin Aによる、「Detecting Malicious PowerShell Commands Using Deep Neural Networks」、コンピュータと通信セキュリティのアジア学会2018の講演論文集[インターネット]、ニューヨーク市、ニューヨーク州、USA:ACM、2018年[2019年8月29日引用]、第187~197頁、(ASIACCS’18);http://doi.acm.org/10.1145/3196494.3196511から入手可能)を、検出することある。
マイクロソフトのパワーシェルは、デフォルトでウインドウズコンピュータにインストールされるコマンドラインツールおよびプログラミング言語であり、システム管理者が、様々な操作タスクのために一般的に使用される。シマンテック社は、最近、サイバ犯罪者が、パワーシェルを攻撃の標的として使用することが増加していると報告した(「PowerShell Threats Grow Further and Operate in Plain Sight」[インターネット][2019年8月29日引用]https://www.symantec.com/blogs/threat-intelligence/powershell-threats-grow-further-and-operate-plain-sightから入手可能)。サイバ犯罪者が検出を回避するために使用できる方法の1つは、いわゆる難読化である。難読化とは、攻撃者のコードは、意図を自動的に検出するのが困難になるように書き直されるが、意図したとおりに実行されることを、意味する。つまり、サイバ攻撃を防ぐための多くの方法の内の1つは、難読化されたパワーシェルスクリプトを検出して削除することである。これは、難読化されたスクリプトの意図は検出できず、信頼できないからである。
本発明による外れ値検出が、既存の種類でのみトレーニングされた場合に、以前には見られなかった種類の難読化をどのように検出できるかを、小規模な実験におけるこの状況のシミュレーションが示している。
この目的のために、双方向LSTMに基づくバイナリ分類子を単一文字のシーケンスでトレーニングして、パワーシェルコーパス(Daniel Bohannonの難読化無効の調査用に集められた、https://github.com/danielbohannon/Revoke-Obfuscation)からの約40万のパワーシェルスクリプトを、難読化されているかどうかを分類できる。例示的なネットワークアーキテクチャを図7に示す。難読化されたスクリプトを構成するには、Daniel BohannonのInvoke-Obfuscation v1.8(https://github.com/danielbohannon/Invoke-Obfuscation)プログラム・ストリング/1、または、その他の(手動または自動、ランダムまたは疑似ランダムまたは決定論的)難読化技術が、通常のスクリプトで使用され、難読化されたスクリプトの新しいセットを作成することができる(以下、難読化1と呼ぶ)。デモンストレーションの目的で、モデルは、Adamオプティマイザ(Kingma DP、Ba J.Adamによる「A Method for Stochastic Optimization」、アーカイヴ14126980Cs[インターネット]、2014年12月22日[2018年7月3日引用]、http://arxiv.org/abs/1412.6980から入手可能)、を使用して、5エポックでトレーニングされ、学習率は、最初の3つのエポックで10-2、第4番目で10-3、第5番目で10-4であり、バイナリ・クロス・エントロピー損失を最小化して、90.4%のテストセット正解率であった。
検査すると、難読化1は非常に穏やかな形式の難読化であることがわかる。つまり、結果として得られる難読化されたスクリプトは、通常のスクリプトと区別するのが難しい場合が多く、モデルが外れ値の検出に使用する有益な特徴を学習したことを示すのに、達成されたテストセット正解率が合理的であり充分であるとの結果になる。
深層学習に基づく外れ値検出モデルを使用して新しい種類の難読化を検出する方法を説明するために、本発明による方法によって達成された結果を、MCドロップアウトによって達成された結果と比較する。外れ値を構成するために、難読化1スクリプトは、難読化実施圧縮プログラムを実行され、難読化されたスクリプトの別のセット(以下、難読化2と呼ぶ)を作成した。外れ値を検出するために、最初の双方向LSTM層から特徴が収集され、シーケンス全体で平均化され(図7参照)、元のパワーシェルスクリプトの変動の99%を補足するPCAモデルが適合された。通常のパワーシェルスクリプトからの任意の逸脱が検索されることがわかっているため、モデルトレーニング中に使用されるすべてのデータではなく、通常のスクリプトのみが使用される。次に、通常のスクリプト、難読化1、および難読化2、の残差平方和(図8Aに示す、Y軸の対数スケールに注意)が、それぞれ計算される。RSSの分布は、通常のスクリプトと難読化1スクリプトとの間で、明確に異なることに注意すべきである。つまり、教師ありの方法でトレーニングされたモデルに対して本発明の方法を使用すると、既知と未知の両方の種類の難読化されたスクリプトを確実に検出できる。
比較のために、10個のモンテカルロサンプルを使用したMCドロップアウトを実行し、予測の標準偏差を外れ値距離として報告する。これにより、通常のスクリプトと難読化されたスクリプトとを強力に分離する(図8Bを参照)。本発明による方法論とMCドロップアウトとの間の大きな違いは、推論時間である。この例では、MCドロップアウトは、未知の外れ値のより強力な識別力を示しているが、本発明による方法は、予測時間に最小限のオーバー・ヘッドを追加する。この例では、本発明による方法は、予測がGPUで加速されている場合でも、予測に必要な実測時間の約1/30を追加するが、外れ値の検出ではそうではない。対照的に、MCドロップアウトは、モンテカルロサンプルごとに1つの完全な予測を必要とする。つまり、この例では、MCドロップアウトは、本発明の方法と比較して300倍長くかかる。バッチ処理を使用して複数のモンテカルロサンプルを並行実行できるが、これにより、推論の実時間と予測に必要なハードウェアメモリとの間に、トレードオフが生じる。
さらなる例
異常検出および/または予知保全のためのデータ分析のための別の例示的な方法は、以下のステップを含む。
S101:データセット(観測値)(例えば、既知の障害または残存寿命を持つ機器のセンサデータの履歴)が提供される。提供された観測値は、トレーニングセットを構成する。
S102:深層ニューラルネットワークをトレーニングするか、事前にトレーニングされたニューラルネットワークを微調整する。深層ニューラルネットワークは、畳み込みニューラルネットワーク、リカレントニューラルネットワーク、リカーシブニューラルネットワーク、トランスフォーマニューラルネットワーク、またはその他の適切な種類のネットワークであり得る。
S103:深層ニューラルネットワークの1つまたは複数のレイヤを選択し、そのレイヤからトレーニングセットの観測値の活性化を計算する。
S103a:任意選択で、活性化が特徴マップからのものである場合(例えば、ニューラルネットワークが畳み込みニューラルネットワークの場合)、各特徴マップ全体にグローバルプーリング(例えば、グローバル平均プーリング)を適用する。
S104:(任意選択でプールされた)トレーニングセットの活性化に潜在変数モデルを適合する。
S105:潜在変数モデルを使用して、トレーニングセットの観測値の1つ以上の距離のセットを計算する。
S106:トレーニングセットの観測値の距離、例えば、トレーニングセットの距離の95パーセンタイル(つまり、トレーニングセットの観測値に対して取得された距離)に基づいて、閾値を選択する。
S107:深層ニューラルネットワークを導入する。言い換えると、異常を検出したり、障害や残存寿命を予測したりするために使用されるソフトウェアに深層ニューラルネットワークを組み込む。ソフトウェアは、監視対象の機器自体の内部、オンプレミスまたはクラウド環境に導入できる中央サーバ(例えば、中央予測サーバまたは異常検出サーバ)、および/または、コンピュータネットワーク上に、独立型ソフトウェアとして実装し得る。
S108:新しい観測値に対する深層ニューラルネットワーク予測を計算する。
S109:次の手順を実行して、特定の新しい観測値が外れ値の観測値であるかどうかを判断する。
S109a:S105で潜在変数モデルを適合するために使用されるレイヤの活性化を計算する。任意選択で、グローバルプーリングを適用し得る。
S109b:S109aからの活性化およびS105の潜在変数モデルを使用して、モデルまでの距離を計算する。
S109c:S109bからの距離を、S106からの閾値距離と比較する。
(i)S109cの距離がS106の閾値よりも大きい場合、新しい観測値は外れ値と見なされる。その後、アプリケーション固有のアクションが実行される。例えば、エンドユーザに、続行する前に外れ値を手動で分析するように警告し得る。S108からの深層ニューラルネットワーク予測は無視され、外れ値の観測値は破棄し得る。
(ii)S109cの距離がS106の閾値より大きくない場合、新しい観測値は外れ値とは見なされず、S108からの予測が受け入れられる。
例示的な予知保全システムは、以下のコンポーネントを含み得る。
予知保全システムは、少なくとも1つの機器などの対象のエンティティの状態を測定するためのデータ収集センサを含み得る。センサは、例えば、温度、圧力、電流、磁束などを測定する物理的センサを含み得る。センサは、物理的測定に基づいて機器の状態が間接的に計算される、いわゆるソフトセンサとも呼ばれ得る。例えば、ソフトセンサには、カルマンフィルタや速度推定器が含まれる。センサは、音声または画像の形式で構造化データまたは非構造化データを収集する場合もあり得る。
予知保全システムは、データ通信ネットワーク、および、例えばデータストアおよび/または他のユニットへの、センサからのデータ転送を可能にするそれぞれのデータ通信ユニットを、さらに備え得る。いくつかの例では、データ転送はリアルタイムで行われ得る。いくつかの例では、データ転送は無線転送を含み得る。
予知保全システムは、センサから収集されたデータ、および任意選択で他のデータが蓄積、処理、および/または分析される、データストア(例えば、データベース)をさらに含み得る。データストアは、オンプレミスサーバ、遠隔サーバ、または分散型クラウドサービスとして実装し得る。
予知保全システムは、本明細書に記載の方法を実行するように構成された予測分析コンポーネント(予測分析ユニット)をさらに含み得る。予測分析コンポーネントは、図1に示される制御および/または分析システム30を構成するか、またはその一部であり得る。特に、予測分析コンポーネントは、収集されたデータを集約およびモデル化して、保守が必要な時期、差し迫ったセキュリティ違反、不正なトランザクションまたはユーザなど、を予測するように構成し得る。予測分析コンポーネントは、観測される機器またはシステムで直接実行でき、いわゆるエッジ分析を構成し得る。予測分析コンポーネントは、サーバ上で実行することも、クラウドサービスとして実行し得る。
予知保全システムは、例えば、オペレータが機器の状態および/または起動に関する見識を得るために使用することができる、ユーザインターフェースをさらに任意選択で備え得る。ユーザインターフェースは、例えば、グラフィカルユーザインターフェースであり得る。
バリエーション
上記の例のいくつかでは、畳み込みニューラルネットワーク(CNN)が使用された。しかしながら、深層ニューラルネットワーク100は、畳み込みニューラルネットワークに限定されず、任意の種類のニューラルネットワーク、例えば、リカレントニューラルネットワーク、リカーシブニューラルネットワーク、トランスフォーマニューラルネットワークなどであり得る。
上記の例のいくつかでは、単一の隠れ層からの活性化が外れ値の検出を実行するために使用される。
他のいくつかの例では、2つ以上の隠れ層からの活性化を使用して外れ値の検出を実行し得る。例えば、外れ値の検出に使用するレイヤを選択する際の困難を回避するために、すべての(非表示の)レイヤからの測定値を、外れ値の検出に使用されているガウス・カーネル密度推定を使用して組み合わせることができる(例えば、L.J.Latecki、A.LazarevicおよびD.Pokrajacによる、「Outlier Detection with Kernel Density Functions」、パターン認識における機械学習とデータマイニング、2007年、第61~75頁、E.Schubert、A.ZimekおよびH.Kriegelによる、「Generalized Outlier Detection with Flexible Kernel Density Estimates」、データマイニングの2014SIAM国際会議の講演論文集、第0巻、産業および応用数学のための学会、2014年、第542~550頁)。カーネル密度推定は、トレーニングセットのマハラノビス距離と残差平方和に対して個別に計算できるが、すべてのレイヤを組み合わせる。各観測値の確率は、モンテカルロ積分を使用して、結果のカーネル密度関数で近似できる。外れ値の検出性能は、上記の実験と同じ方法で評価できる。
さらなる例では、複数の埋め込み(例えば、複数のレイヤからの活性化)を統合する場合、バギング・アプローチ(例えば、A.LazarevicおよびV.Kumarによる「Feature Bagging for Outlier Detection」、第11回ACM SIGKDDデータマイニングにおける知識発見国際会議の講演論文集、ニューヨーク市、ニューヨーク州、USA、2005年、第157~166頁を参照)を、直接カーネル密度の代わりに、使用し得る。
さらに、学習された表現、GMMワンクラスSMV、クラスタリングなどで他の方法論を使用できる場合がある。
さらに、本明細書に記載の様々な実施形態および例では、マハラノビス距離などの単純な測定値を使用する代わりに、局所外れ値因子またはLOFなどのより局所的な測定値(例えば、M.M.Breunig、H.-P.Kriegel、R.T.NgおよびJ.Sanderによる「LOF:Identifying Density-based Local Outliers」、2000ACM SIGMODデータ会議の国際会議の講演論文集、ニューヨーク市、ニューヨーク州、USA、2000年、第93~104頁)は、PCA空間に直接適用され得る。原則として、複合指標のモデル化された埋め込みまでの距離、平方和、局所外れ値因子など、任意の適切な距離指標、または距離指標の組合せ、が使用され得る。
上記の例のいくつかでは、特定の画像が観測値として(トレーニングセットの一部として、および予測が取得される新しい観測値として)、使用される。上記のように、故障、劣化、完全性またはセキュリティの違反、または観測されるエンティティの正常な状態または操作からの任意の他の逸脱(すなわち、異常)が判断さ得ることに基づいて、これらの画像の代わりに、例えば機器、機械部品の画像など任意の他の画像が使用され得る。
例えば、特定の特徴の存在、画像の破損のレベル、取得された画像の色の変化などの取得された画像の特性に基づいて、画像化されたオブジェクトおよび/または画像システムの異常を、判断し得る。画像化されたオブジェクトの異常には、例えば、障害(例えば、亀裂)、画像化されたオブジェクトの摩耗または劣化が含まれ得る。さらに、画像の欠陥または画像の破損の存在は、画像化システムの欠陥、不整合などを示し得る。例えば、画像は、画像書き込み中の書き込みプロセスの中断、画像キャプチャデバイスの集束システムの障害のために、少なくとも部分的に破損する可能性があり、そのような異常を検出すると、ユーザは個々に通知されるか、および/または、修正保守アクションが開始され得る。
さらに、画像の代わりに、他の種類の入力データ、例えば、1つまたは複数のセンサの測定に基づいて得られたデータ、コンピュータネットワークトラフィックデータ、トランザクションデータ(例えば、金融取引に関連する)などが、使用され得る。原則として、観測される機器またはシステムを特徴付ける任意の種類のデータが、処理され得る。次に、本明細書に記載の様々な実施形態および例による方法を使用して、正常な動作と既知の種類の異常とを区別するだけでなく、未知の種類の異常にもフラグを立てることができる。
さらに、本明細書に記載されるような方法およびシステムの様々な用途があり得る。可能な用途の1つは、位相差顕微鏡で得られた画像の異常を検出することである。位相コントラスト画像から核を識別する作業は困難であり、通常、核の位置を強調する追加の画像チャネルを提供するために、核に結合する蛍光マーカに依存する。深層CNNを適用することにより、細胞からバックグラウンドをセグメント化し、蛍光ラベルを使用せずに位相差画像から細胞の核を識別するモデルが、生成およびトレーニングされ得る。CNNが複数の細胞種類からの位相差画像でトレーニングされている場合、CNNは、異常と見なされる可能性のある、これまでに見たことがない新しい細胞種類も識別し得る。これまで見たことがない細胞種類が、CNNトレーニングに使用されるセル種類と大きく異なる場合を判断するために、予測時間外れ値検出が、適用され得る。新しい細胞種類からの位相差顕微鏡画像が外れ値として検出された場合、異常が検出された、および/またはCNNの細胞とバックグラウンドの分離、および核の検出が信頼できない、という警告を、ユーザが受信し得る。この場合、提供されるCNNシステムは、細胞をバックグラウンドから確実に分離し、新しい細胞種類の位相差画像から核を識別するために、CNNモデルを再トレーニングする必要があるという信号を送信し得る。
さらなる用途は、機械部品、機械、産業プラント、発電所などの様々な種類の技術システムでの異常検出である。
例えば、予知保全機器では、通常、物理センサを使用してリアルタイムで測定できる。次に、予測深層学習モデルを使用して、エラーが発生する前に、エラーを検出および分類できる。本明細書に記載の方法およびシステムは、既知のエラーに加えて、未知の種類のエラーを検出する可能性を提供し、未知の種類のエラーは、外れ値として確実に特定される。その結果に基づいて、機器の破損を防ぐために保守が必要な時期が、自動的に判断され得る。したがって、予知保全システムは、これまで未知であった障害種類が発生した場合でも、定義された動作を行い得る。したがって、大幅なコスト削減を達成でき、機器が予想よりも早く故障して、致命的な障害やプロセスの停止につながる可能性があるリスクを減らし得る。
もう1つの例は、サイバーセキュリティの技術分野であり、例えば、ネットワークへの侵入を検出する。例えば、深層学習を使用してネットワークトラフィック、つまりネットワーク上の異なるコンピュータ間で転送されるデータパッケージのログをモデル化することにより、ネットワーク内の悪意のある攻撃者を示すトラフィックパターンの微妙な変化が検出され、侵入として分類し得る。提案された方法とシステムを使用することにより、以前に見られなかったパターンがある場合にも、使用された深層学習モードは警告することができる。そのようなパターンは外れ値として確実に識別されるからである。
別の例は、金融取引(例えば、銀行)の分野である。金融コンピュータネットワークでは、経済取引は、例えば、いつ/どこで/どのくらいの金額がアクタ間で転送されるかによって、特徴付け得る。履歴データを使用して、深層学習に基づく分類器をトレーニングして、トランザクションまたはアクタを不正として分類し得る。本発明を使用して、新しい種類の不正なトランザクションが発生したとき、またはアクタが以前には見られなかった方法で行動したときに、深層学習モデルは、警告することができる。したがって、金融コンピュータネットワークのセキュリティを改善し得る。
上記の例では、履歴データを使用して、調査対象のシステムの特性を明らかにし、深層学習ニューラルネットワーク(深層学習モデル)をトレーニングし得る。次に、深層学習モデルを使用して、正常な動作と既知の種類の異常を区別し得る。さらに、提案された方法とシステムは、未知の種類の異常を識別する可能性が増える。
ハードウェア構成
図9は、コンピューティングシステム1の少なくとも一部を実装するために、および/または、上記の方法を実行するために使用され得るコンピュータの例示的なハードウェア構成を示す。例えば、コンピュータは、上記の予測分析コンポーネントを実装し得る。図9に示されるコンピュータ7は、中央処理装置(CPU)70、グラフィックス処理装置(GPU)88、システムメモリ72、ネットワークインターフェース74、ハードディスクドライブ(HDD)インターフェース76、外部ディスクドライブインターフェース78および入力/出力(I/O)インターフェース80を含む。コンピュータのこれらのコンポーネントは、システムバス82を介して互いに結合されている。CPU70は、システムメモリ72にアクセスすることによって、算術、論理、および/または制御操作を実行し得る。GPU88は、深層学習に必要なグラフィック画像および/または行列演算に関する処理を実行し得る。GPU88は、並列処理を実行し得る複数のプロセッサ(例えば、コア)を備えることができ、これは、コンピューティングシステム1のより高い性能になり得る。CPU70および/またはGPU88は、上記の例示的なデバイスおよび/またはシステムのプロセッサを実装し得る。しかしながら、いくつかの例では、コンピュータ7は必ずしもGPU88を備えているとは限らない。システムメモリ72は、CPU70と組み合わせて使用するための情報および/または命令を格納し得る。システムメモリ72は、ランダムアクセスメモリ(RAM)720および読み取り専用メモリ(ROM)722などの揮発性および不揮発性メモリを含み得る。起動時など、コンピュータ7内の要素間で情報を転送するのに役立つ基本ルーチンを含む基本入出力システム(BIOS)は、ROM722に格納され得る。システムバス82は、メモリバスまたはメモリコントローラ、周辺バス、および様々なバスアーキテクチャのいずれかを使用するローカルバスを含む、いくつかの種類のバス構造のいずれかであり得る。CPU70は、1つまたは複数の対応するインターフェース(図示せず)およびバス82を介して、1つまたは複数のセンサ(図示せず)にさらに接続し得る。センサは、温度、pH、圧力などを含むが、これらに限定されない、物理的状況または状態を測定し得る。さらに、センサは、画像化デバイス、マイクロフォン、スペクトル・センサなどを含むが、これらに限定されない、他の種類の測定または検出デバイスを含み得る。コントローラは、温度、磁束、攪拌などを含むが、これらに限定されない、物理的状況または状態を制御し得る。
コンピュータは、ネットワークを介して他のコンピュータおよび/またはデバイスと通信するためのネットワークインターフェース74を含み得る。
さらに、コンピュータは、ハードディスク(図示せず)からの読み取りおよび書き込みのためのハードディスクドライブ(HDD)84(図示せず)、および、リムーバブル・ディスク(図示せず)への読み取りまたは書き込みのための外部ディスクドライブ86を含み得る。リムーバブルディスクは、磁気ディスクドライブ用の磁気ディスク、または光ディスクドライブ用のCD-ROMなどの光ディスクであり得る。HDD84および外部ディスクドライブ86は、それぞれ、HDDインターフェース76および外部ディスクドライブインターフェース78によってシステムバス82に接続されている。ドライブとそれに関連するコンピュータ可読媒体は、汎用コンピュータ用のコンピュータ可読命令、データ構造、プログラムモジュール、およびその他のデータの不揮発性記憶を提供する。データ構造は、本明細書に記載されるように、例示的な方法およびその変形の実施に関連するデータを含み得る。関連するデータは、データベース、例えばリレーショナルデータベースまたはオブジェクトデータベースに編成し得る。
本明細書に記載の例示的な環境は、ハードディスク(図示せず)および外部ディスク(図示せず)を使用するが、コンピュータによってアクセス可能なデータを格納できる他の種類のコンピュータ可読媒体、磁気カセット、フラッシュ・メモリカード、デジタル・ビデオ・ディスク、ランダムアクセスメモリ、読み取り専用メモリなどのようなものも、例示的な動作環境で使用することができることは、当業者によって理解されるべきである。
いくつかのプログラムモジュールが、ハードディスク、外部ディスク、ROM722またはRAM720に格納され得る。そこには、オペレーティングシステム(図示せず)、1つまたは複数のアプリケーションプログラム7202、他のプログラムモジュール(図示せず)、およびプログラムデータ7204を含む。アプリケーションプログラムは、上記記載のような機能の少なくとも一部を含み得る。
コンピュータ7は、対応するI/Oインターフェース80aおよび80bならびにシステムバス82を介して、マウスおよび/またはキーボードなどの入力装置92および液晶ディスプレイなどの表示装置94が、接続され得る。コンピュータ7がタブレットコンピュータとして実装される場合、例えば、情報を表示し、入力を受信するタッチパネルは、対応するI/Oインターフェースおよびシステムバス82を介して、コンピュータ7に接続され得る。さらに、いくつかの例では、図9には示されていないが、コンピュータ7は、対応するI/Oインターフェースおよびシステムバス82を介して、プリンタおよび/またはカメラなどの画像化デバイスにさらに接続され得る。
図9に示されるようなコンピュータ7を使用する実装に加えて、またはその代替として、本明細書に記載の例示的な実施形態の機能の一部またはすべては、1つまたは複数のハードウェア回路として実装され得る。このようなハードウェア回路の例には、大規模集積回路(LSI)、縮小命令セット回路(RISC)、特定用途向け集積回路(ASIC)、およびフィールドプログラマブルゲートアレイ(FPGA)が含まれるが、これらに限定されない。
ASIC 特定用途向け集積回路
BIOS 基本入出力システム
CD ROM コンパクトディスク読み取り専用メモリ
CNN 畳み込みニューラルネットワーク
CPU 中央処理装置
DNN 深層ニューラルネットワーク
GPU グラフィックス処理装置
FPGA フィールドプログラマブルゲートアレイ
HDD ハードディスクドライブ
HTTP ハイパーテキスト転送プロトコル
I/O入力/出力
LSI 大規模集積回路
LOF 局所外れ値因子
MCドロップアウト モンテカルロドロップアウト
MNIST 修正されたアメリカ国立標準技術研究所
PCA 主成分分析
RAM ランダムアクセスメモリ
ReLu/Relu 正規化線形ユニット
RISC 縮小命令セット回路
ROC-AUC 受信者操作特性、曲線下面積
ROM 読み取り専用メモリ
RSS 残差平方和
SSH セキュアシェル

Claims (15)

  1. 対象のエンティティにおける異常検出のためのコンピュータ実装方法であって、
    前記コンピュータ実装方法は、
    新しい観測値を受信することであって、前記新しい観測値は、前記エンティティの少なくとも1つのパラメータを特徴付けることと、
    前記新しい観測値を深層ニューラルネットワーク(100)に入力することであって、前記深層ニューラルネットワーク(100)は、複数の隠れ層を有し、前記深層ニューラルネットワーク(100)に入力できる可能性のある観測値を含むトレーニングデータセットを使用してトレーニングされることと、
    前記受信した新しい観測値を前記深層ニューラルネットワーク(100)に入力することにより、前記深層ニューラルネットワーク(100)の前記複数の隠れ層のうちの少なくとも1つから出力される中間出力値の第2のセットを取得することと、
    記憶媒体に格納された潜在変数モデルを使用して、前記中間出力値の第2のセットを射影値の第2のセットにマッピングすることと、
    前記潜在変数モデルおよび前記射影値の第2のセットに基づいて、前記受信した新しい観測値が前記トレーニングデータセットに関して外れ値であるかどうかを判断することと、
    前記深層ニューラルネットワーク(100)によって、前記新しい観測値の予測を計算することと、
    前記予測および前記新しい観測値が外れ値であるかどうかの前記判断に基づいて、前記エンティティに少なくとも1つの前記異常が発生したことを示す結果を判断することとを含み、
    前記記憶媒体に格納された前記潜在変数モデルは、
    前記深層ニューラルネットワーク(100)の前記複数の隠れ層のうちの前記1つから出力される中間出力値の第1のセットを取得することであって、前記中間出力値の第1のセットのそれぞれは、前記トレーニングデータセットの少なくとも一部に含まれる前記可能性のある観測値の異なる1つを入力することによって取得されることと、
    前記中間出力値の第1のセットを使用して前記潜在変数モデルを構築することであって、前記潜在変数モデルは、前記中間出力の前記セットの次元よりも低い次元を有する前記潜在変数モデルの部分空間内の射影値の第1のセットへ前記中間出力値の第1のセットのマッピングを提供することと、によって構築される潜在変数モデルである、
    コンピュータ実装方法。
  2. データを処理するための前記深層ニューラルネットワーク(100)および前記深層ニューラルネットワーク(100)をトレーニングするために使用されるトレーニングデータセットの少なくとも一部を取得することであって、前記深層ニューラルネットワークは、複数の隠れ層を含み、前記トレーニングデータセットを使用してトレーニングされ、前記トレーニングデータセットは、前記深層ニューラルネットワークに入力できる可能性のある複数の観測値を含むことと、
    前記複数の隠れ層のうちの少なくとも1つから出力される中間出力値の第1のセットを取得することであって、前記中間出力値の第1のセットのそれぞれは、前記トレーニングデータセットの前記少なくとも一部に含まれる前記可能性のある観測値の異なる1つを入力することによって取得できることと、
    前記中間出力値の第1のセットを使用して潜在変数モデルを構築または適合することであって、前記潜在変数モデルは、前記中間出力の前記セットの次元よりも低い次元を有する部分空間内の射影値の第1のセットへの前記中間出力値の第1のセットのマッピングを提供することと、
    前記潜在変数モデルおよび前記射影値の第1のセットを、記憶媒体に格納することと、
    をさらに含む、請求項1に記載のコンピュータ実装方法。
  3. 前記受信した新しい観測値が外れ値であるかどうかを判断するステップは、
    前記射影値の第1のセットの分布に関して、前記射影値の第2のセットの距離を判断することと、
    前記計算された距離が前記判断された距離の閾値よりも大きい場合、前記受信した新しい観測値が前記トレーニングデータセットに関して外れ値であると判断することとを、
    含む、請求項1または2に記載の方法。
  4. 前記距離の前記閾値は、距離に基づいて決定され、前記距離のそれぞれは、前記射影値の第1のセットの前記分布に関して前記射影値の第1のセットの異なる1つに対して計算される、請求項3に記載の方法。
  5. 前記距離が、複合距離、残差平方和、マハラノビス距離、または局所外れ値因子のうちの1つである、請求項3または4に記載の方法。
  6. 前記受信した新しい観測値が外れ値であるかどうかを判断するステップは、
    前記潜在変数モデルおよび前記射影値の第2のセットを使用して、前記中間出力値の第2のセットに対応する中間出力値の近似セットを判断することと、
    前記中間出力値の第2のセットおよび前記中間出力値の近似セットに対して、二乗近似残差を計算することと、
    前記計算された二乗近似残差が前記二乗近似残差の閾値よりも大きい場合、前記受信した観測値が前記トレーニングデータセットに関して外れ値であると判断することと、
    を含む、請求項1から5のいずれか一項に記載の方法。
  7. 前記二乗近似残差の前記閾値は、二乗近似残差に基づいて判断され、前記二乗近似残差のそれぞれは、前記中間出力値の第1のセットの異なる1つ、および前記中間出力値の第1のセットの一つに対応する中間出力値の近似セットについて計算される、請求項6に記載の方法。
  8. 前記中間出力値の第1のセットを取得し、前記潜在変数モデルを構築するステップは、前記複数の隠れ層のうちの2つ以上に対して実行され、
    前記中間出力値の第2のセットを取得し、前記中間出力値の前記第2のセットを前記射影値の第2のセットにマッピングするステップは、前記複数の隠れ層のうちの前記2つ以上に関して実行され、
    前記受信した新しい観測値が外れ値であるかどうかを判断するステップは、前記潜在変数モデル、および前記複数の隠れ層の前記2つ以上に関して取得された前記射影値の第2のセットに基づいて実行される、
    請求項1から7のいずれか一項に記載の方法。
  9. 前記深層ニューラルネットワーク(100)の前記複数の隠れ層のうちの少なくとも1つから出力される前記中間出力値を取得するステップは、前記複数の隠れ層の前記の少なくとも1つからの活性化を判断することと、および、任意選択で、前記判断された活性化にグローバルプーリングを適用することと、を含む、請求項1から8のいずれか一項に記載の方法。
  10. 前記潜在変数モデルが、主成分分析に従って、またはオートエンコーダを使用して構築される、請求項1から9のいずれか一項に記載の方法。
  11. 前記観測値は、1つまたは複数のセンサによる物理的測定に基づいて取得される、請求項1から10のいずれか一項に記載の方法。
  12. エンティティの予知保全のためのコンピュータ実装方法であって、
    前記エンティティに対して実行された少なくとも1つのセンサ測定からのセンサデータに基づいて新しい観測値を取得し、
    請求項1から11のいずれか一項の方法に従って、前記エンティティの少なくとも1つの前記異常の発生を示す結果を判断し、
    少なくとも1つの前記異常の発生を示す前記結果に基づいて、修正処置を講じる必要があるかどうか、および、任意選択で講じるべき前記修正処置の種類を判断する、
    コンピュータ実装方法。
  13. コンピュータにロードされて実行されると、請求項1から12のいずれか一項に記載の方法を前記コンピュータに実行させる、コンピュータ可読命令を含むコンピュータプログラム製品。
  14. エンティティの異常検出または予知保全のためのシステムであって、
    前記エンティティの少なくとも1つのパラメータを特徴付ける1つまたは複数の観測値を取得するためのデータ収集ユニットと、
    深層ニューラルネットワーク(100)のトレーニングに使用されるトレーニングデータセットを格納する記憶媒体(12)であって、前記深層ニューラルネットワーク(100)は、複数の隠れ層を含み、前記トレーニングデータセットを使用してトレーニングされ、前記トレーニングデータセットは、前記深層ニューラルネットワーク(100)に入力できる可能性のある観測値を含む、記憶媒体(12)と、
    請求項1から12のいずれか一項に記載の方法を実行するように構成された少なくとも1つのプロセッサを含む予測分析ユニットと、
    を含む、システム。
  15. 前記エンティティの少なくとも1つのパラメータを測定するための少なくとも1つのセンサと、
    前記少なくとも1つのセンサと、前記データ収集ユニット、および/または前記少なくとも1つのプロセッサとの間のデータ転送のためのデータ通信ユニットと、
    をさらに含む、請求項14に記載のシステム。
JP2021512922A 2018-09-05 2019-09-05 異常検出および/または予知保全のためのコンピュータ実装方法、コンピュータプログラム製品およびシステム Active JP7223839B2 (ja)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
EP18192649.4A EP3620983B1 (en) 2018-09-05 2018-09-05 Computer-implemented method, computer program product and system for data analysis
EP18192649.4 2018-09-05
EP19180972.2A EP3620986A1 (en) 2018-09-05 2019-06-18 Computer-implemented method, computer program product and system for analysis of cell images
EP19180972.2 2019-06-18
PCT/EP2019/073670 WO2020049087A1 (en) 2018-09-05 2019-09-05 Computer-implemented method, computer program product and system for anomaly detection and/or predictive maintenance

Publications (3)

Publication Number Publication Date
JP2022500745A JP2022500745A (ja) 2022-01-04
JPWO2020049087A5 JPWO2020049087A5 (ja) 2022-10-07
JP7223839B2 true JP7223839B2 (ja) 2023-02-16

Family

ID=63517723

Family Applications (2)

Application Number Title Priority Date Filing Date
JP2021512922A Active JP7223839B2 (ja) 2018-09-05 2019-09-05 異常検出および/または予知保全のためのコンピュータ実装方法、コンピュータプログラム製品およびシステム
JP2021512919A Active JP7187681B2 (ja) 2018-09-05 2019-09-05 細胞画像の分析のためのコンピュータ実装方法、コンピュータプログラム製品およびシステム

Family Applications After (1)

Application Number Title Priority Date Filing Date
JP2021512919A Active JP7187681B2 (ja) 2018-09-05 2019-09-05 細胞画像の分析のためのコンピュータ実装方法、コンピュータプログラム製品およびシステム

Country Status (5)

Country Link
US (2) US12001935B2 (ja)
EP (4) EP3620983B1 (ja)
JP (2) JP7223839B2 (ja)
CN (3) CN112639834A (ja)
WO (3) WO2020049098A1 (ja)

Families Citing this family (102)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11132529B2 (en) * 2016-11-16 2021-09-28 Ventana Medical Systems, Inc. Convolutional neural networks for locating objects of interest in images of biological samples
US10923141B2 (en) 2018-08-06 2021-02-16 Spotify Ab Singing voice separation with deep u-net convolutional networks
US10977555B2 (en) * 2018-08-06 2021-04-13 Spotify Ab Automatic isolation of multiple instruments from musical mixtures
EP3620983B1 (en) * 2018-09-05 2023-10-25 Sartorius Stedim Data Analytics AB Computer-implemented method, computer program product and system for data analysis
US12001949B2 (en) 2018-09-05 2024-06-04 Sartorius Stedim Data Analytics Ab Computer-implemented method, computer program product and system for data analysis
KR102608470B1 (ko) * 2018-10-23 2023-12-01 삼성전자주식회사 데이터 인식 장치 및 방법과 트레이닝 장치 및 방법
US11520390B2 (en) * 2018-11-07 2022-12-06 Hewlett-Packard Development Company, L.P. Receiving thermal data and producing system thermal grades
US10885631B2 (en) * 2019-02-01 2021-01-05 Essen Instruments, Inc. Label-free cell segmentation using phase contrast and brightfield imaging
US11682069B2 (en) * 2019-08-05 2023-06-20 Intuit, Inc. Extending finite rank deep kernel learning to forecasting over long time horizons
WO2021040779A1 (en) * 2019-08-23 2021-03-04 Landmark Graphics Corporation Anomaly detection using hybrid autoencoder and gaussian process regression
US11409995B2 (en) * 2019-08-28 2022-08-09 Nxp B.V. Anomaly detection by classifying past behavior
US20210141355A1 (en) * 2019-11-07 2021-05-13 Global Energy Interconnection Research Institute Co. Ltd Systems and methods of autonomous line flow control in electric power systems
US20210182386A1 (en) * 2019-12-11 2021-06-17 Samsung Electronics Co., Ltd. Electronic apparatus that monitors a safety function and a controlling method thereof
JP7433876B2 (ja) * 2019-12-12 2024-02-20 株式会社メガチップス 自己組織化マップの学習装置、方法およびプログラム並びに状態判定装置
US11423365B2 (en) * 2020-02-17 2022-08-23 Mo Tecnologias, Llc Transaction card system having overdraft capability
US11973774B2 (en) * 2020-02-28 2024-04-30 Darktrace Holdings Limited Multi-stage anomaly detection for process chains in multi-host environments
US11494496B2 (en) * 2020-03-30 2022-11-08 International Business Machines Corporation Measuring overfitting of machine learning computer model and susceptibility to security threats
US20210312406A1 (en) * 2020-04-07 2021-10-07 Dgnss Solutions, Llc Artificial intelligence monitoring, negotiating, and trading agents for autonomous vehicles
US20230300354A1 (en) * 2020-06-25 2023-09-21 Telefonaktiebolaget Lm Ericsson (Publ) Method and System for Image Compressing and Coding with Deep Learning
CN111814653B (zh) * 2020-07-02 2024-04-05 苏州交驰人工智能研究院有限公司 一种视频中异常行为的检测方法、装置、设备及存储介质
US11475332B2 (en) * 2020-07-12 2022-10-18 International Business Machines Corporation Selecting forecasting models by machine learning based on analysis of model robustness
CN111970489B (zh) * 2020-08-05 2022-05-24 北京必可测科技股份有限公司 一种基于人机双向的智能监盘管理方法及***
CN112001482B (zh) * 2020-08-14 2024-05-24 佳都科技集团股份有限公司 振动预测及模型训练方法、装置、计算机设备和存储介质
EP3957712A1 (en) 2020-08-21 2022-02-23 Sartorius Stedim Data Analytics AB Automated foam detection
JP7318612B2 (ja) * 2020-08-27 2023-08-01 横河電機株式会社 監視装置、監視方法、および監視プログラム
SG10202008469RA (en) * 2020-09-01 2020-10-29 Ensign Infosecurity Pte Ltd A deep embedded self-taught learning system and method for detecting suspicious network behaviours
US11360155B2 (en) * 2020-09-03 2022-06-14 Volkswagen Aktiengesellschaft Battery state of health estimator
IL277424B1 (en) * 2020-09-16 2024-03-01 Israel Aerospace Ind Ltd Predicting events based on machine learning and engineering analysis tools
CN112541558A (zh) * 2020-09-18 2021-03-23 江苏师范大学 一种基于不完整数据的贝叶斯半监督鲁棒ppls软测量方法
CN112149912A (zh) * 2020-09-30 2020-12-29 中煤能源研究院有限责任公司 一种基于逆向传播神经网络的煤矿工作面矿压数据预测方法
US20220101068A1 (en) * 2020-09-30 2022-03-31 International Business Machines Corporation Outlier detection in a deep neural network using t-way feature combinations
US20220108097A1 (en) * 2020-10-05 2022-04-07 Rakuten, Inc. Dual encoder attention u-net
US20220121691A1 (en) * 2020-10-20 2022-04-21 Here Global B.V. Method, apparatus, and system for machine learning-based persistence filtering
CN112348839B (zh) * 2020-10-27 2024-03-15 重庆大学 一种基于深度学习的图像分割方法及***
CN112015749B (zh) * 2020-10-27 2021-02-19 支付宝(杭州)信息技术有限公司 基于隐私保护更新业务模型的方法、装置及***
CN112953924B (zh) * 2021-02-04 2022-10-21 西安电子科技大学 网络异常流量检测方法、***、存储介质、终端及应用
CN112861443B (zh) * 2021-03-11 2022-08-30 合肥工业大学 一种融入先验知识的深度学习故障诊断方法
CN113256507B (zh) * 2021-04-01 2023-11-21 南京信息工程大学 一种针对二进制流量数据生成图像的注意力增强方法
JP2022162903A (ja) * 2021-04-13 2022-10-25 パナソニックIpマネジメント株式会社 生産設備監視システムおよび生産設備監視方法
CN112967816B (zh) * 2021-04-26 2023-08-15 四川大学华西医院 一种急性胰腺炎器官衰竭预测方法、计算机设备和***
CN113255209B (zh) * 2021-04-27 2022-06-07 燕山大学 一种齿轮箱轴承剩余寿命的预测方法
US11893085B2 (en) 2021-05-04 2024-02-06 Ford Global Technologies, Llc Anomaly detection for deep neural networks
CN112987675B (zh) * 2021-05-06 2021-07-23 北京瑞莱智慧科技有限公司 一种异常检测的方法、装置、计算机设备和介质
CN113159225B (zh) * 2021-05-17 2022-04-19 中国石油大学(华东) 多变量工业过程故障分类方法
EP4347784A1 (en) * 2021-05-27 2024-04-10 Lynceus SAS Machine learning-based quality control of a culture for bioproduction
WO2022248935A1 (en) * 2021-05-27 2022-12-01 Lynceus Sas Machine learning-based quality control of a culture for bioproduction
CN113297677B (zh) * 2021-06-15 2023-09-22 中国航发沈阳发动机研究所 一种基于Kriging模型的飞机发动机空气***概率分析方法
CN113762333B (zh) * 2021-07-20 2023-02-28 广东省科学院智能制造研究所 一种基于双流联合密度估计的无监督异常检测方法和***
CN113596012A (zh) * 2021-07-26 2021-11-02 工银科技有限公司 攻击行为的识别方法、装置、设备、介质及程序产品
CN114740159B (zh) * 2022-04-14 2023-09-19 成都秦川物联网科技股份有限公司 一种天然气能量计量组分获取方法和物联网***
US11979697B2 (en) 2021-07-26 2024-05-07 Chengdu Qinchuan Iot Technology Co., Ltd. Methods and internet of things systems for obtaining natural gas energy metering component
US11630428B2 (en) * 2021-08-06 2023-04-18 Trimble Inc. Real-time analysis of vibration samples for operating environment classification and anomaly detection
CN113689040B (zh) * 2021-08-25 2024-03-15 长江存储科技有限责任公司 测量模型的控制方法和控制***及计算机可读介质
CN113836673A (zh) * 2021-09-06 2021-12-24 西安理工大学 一种基于信息熵的排水管网监测点的布置方法
CN113510234B (zh) * 2021-09-14 2022-01-07 深圳市信润富联数字科技有限公司 轮毂低压铸造的质量监控方法、装置及电子设备
US20230086049A1 (en) * 2021-09-21 2023-03-23 Abb Schweiz Ag Systems and methods of predicting the remaining useful life of industrial mechanical power transmission equipment using a machine learning model
US20230133652A1 (en) * 2021-10-29 2023-05-04 GE Grid GmbH Systems and methods for uncertainty prediction using machine learning
US20230136726A1 (en) * 2021-10-29 2023-05-04 Peter A. Chew Identifying Fringe Beliefs from Text
TWI809592B (zh) * 2021-12-08 2023-07-21 財團法人工業技術研究院 模型預測控制系統及其方法
JP7188821B1 (ja) 2021-12-17 2022-12-13 バイスリープロジェクツ株式会社 検査システム
CN114237128A (zh) * 2021-12-21 2022-03-25 华能澜沧江水电股份有限公司 基于趋势报警的水电站设备实时监控数据监控***及监控方法
CN114564877B (zh) * 2021-12-29 2024-05-10 西安交通大学 一种滚动轴承寿命预测方法、***、设备及可读存储介质
EP4206838A1 (en) * 2021-12-29 2023-07-05 Petkim Petrokimya Holding A.S. Forecasting and anomaly detection method for low density polyethylene autoclave reactor
CN114354534A (zh) * 2021-12-30 2022-04-15 中国航空油料有限责任公司 利用二元线性分类器建立航空煤油性质预测模型的方法
TWI829076B (zh) * 2022-01-22 2024-01-11 國立高雄科技大學 藉由機器學習輔助模擬壓力數值進行虛實校正的學習方法
CN114528874A (zh) * 2022-01-24 2022-05-24 五凌电力有限公司 水电机组的故障诊断方法、装置及电子设备
DE102022103844B3 (de) 2022-02-17 2023-06-22 Synsor.ai GmbH Verfahren zur Optimierung eines Produktionsprozesses auf Basis visueller Informationen und Vorrichtung zur Durchführung des Verfahrens
US12007832B2 (en) 2022-02-25 2024-06-11 Bank Of America Corporation Restoring a system by load switching to an alternative cloud instance and self healing
CN114925196B (zh) * 2022-03-01 2024-05-21 健康云(上海)数字科技有限公司 多层感知网络下糖尿病血检异常值辅助剔除方法
JP7444186B2 (ja) * 2022-03-22 2024-03-06 横河電機株式会社 モデル検証装置、モデル検証方法、および、モデル検証プログラム
CN114726614B (zh) * 2022-04-06 2023-05-23 广东工业大学 一种基于条件变分自编码器和SENet的XSS攻击检测方法
CN114692694B (zh) * 2022-04-11 2024-02-13 合肥工业大学 一种基于特征融合和集成聚类的设备故障诊断方法
CN114978586B (zh) * 2022-04-12 2023-07-04 东北电力大学 一种基于攻击基因的电网攻击检测方法、***和电子设备
CN114818904B (zh) * 2022-04-21 2024-03-15 桂林电子科技大学 基于Stack-GANs模型的风机故障检测方法及存储介质
DE102022205534A1 (de) * 2022-05-31 2023-11-30 Kuka Deutschland Gmbh Überwachen einer mehrachsigen Maschine mittels interpretierbarer Zeitreihenklassifikation
CN115147203B (zh) * 2022-06-08 2024-03-15 阿尔法时刻科技(深圳)有限公司 基于大数据的金融风险分析方法
CN115051864B (zh) * 2022-06-21 2024-02-27 郑州轻工业大学 基于pca-mf-wnn的网络安全态势要素提取方法及***
US20240005132A1 (en) * 2022-06-30 2024-01-04 Maplebear Inc. (Dba Instacart) Machine-learned neural network architectures for incremental lift predictions
WO2024077271A2 (en) * 2022-10-07 2024-04-11 National Resilience, Inc. Biosecure digital twin for cyber-physical anomaly detection and biological process modeling
CN115454781B (zh) * 2022-10-08 2023-05-16 杭银消费金融股份有限公司 基于企业架构***的数据可视化展现方法及***
WO2024086143A1 (en) * 2022-10-17 2024-04-25 Proteinea Inc. Native expansion of a sparse training dataset into a dense training dataset for supervised training of a synonymous variant sequence generator
CN116112193B (zh) * 2022-10-18 2023-07-28 贵州师范大学 一种基于深度学习的轻量级车载网络入侵检测方法
EP4362029A1 (en) 2022-10-24 2024-05-01 Sartorius Stedim Data Analytics AB Modelling of biological and biochemical assays
DE102022212232A1 (de) 2022-11-17 2024-05-23 Robert Bosch Gesellschaft mit beschränkter Haftung Verfahren und Vorrichtung zum Erstellen eines robusten Anomalieerkennungsmodells zum Erkennen einer aktuellen oder anstehenden Anomalie einer Gerätebatterie mithilfe erlernter Latent Space Embeddings
CN115859202B (zh) * 2022-11-24 2023-10-10 浙江邦盛科技股份有限公司 一种非平稳时序数据流场景下的异常检测方法及装置
CN116187540B (zh) * 2023-01-03 2023-09-12 国网江西省电力有限公司电力科学研究院 一种基于时空偏差修正的风电场站超短期功率预测方法
CN116089611B (zh) * 2023-01-13 2023-07-18 北京控制工程研究所 基于性能-故障关系图谱的航天器故障诊断方法及装置
CN116048031B (zh) * 2023-03-31 2023-08-04 克拉玛依市蓝润环保科技有限责任公司 用于石油助剂生产的控制***及其方法
CN116258579B (zh) * 2023-04-28 2023-08-04 成都新希望金融信息有限公司 用户信用评分模型的训练方法及用户信用评分方法
CN116382223B (zh) * 2023-06-02 2023-08-01 山东鲁能控制工程有限公司 一种基于dcs的火电机组监测***
CN117077062A (zh) * 2023-08-31 2023-11-17 木卫四(北京)科技有限公司 汽车指令异常检测方法及装置
CN116862908B (zh) * 2023-09-01 2023-11-24 菏泽学院 一种基于图像处理的金属齿轮啮合缺陷检测方法
CN116908808B (zh) * 2023-09-13 2023-12-01 南京国睿防务***有限公司 一种基于rtn的高分辨一维像目标识别方法
CN117151196B (zh) * 2023-10-26 2024-01-30 湘江实验室 一种基于逐层增量期望传播的可解释性故障诊断归因方法
CN117192063B (zh) * 2023-11-06 2024-03-15 山东大学 基于耦合卡尔曼滤波数据同化的水质预测方法及***
CN117709536B (zh) * 2023-12-18 2024-05-14 东北大学 一种深度递归随机配置网络工业过程精准预测方法和***
CN117436024B (zh) * 2023-12-19 2024-03-08 湖南翰文云机电设备有限公司 一种基于钻机运行数据分析的故障诊断方法及***
CN117435511B (zh) * 2023-12-20 2024-03-19 厦门市唯云网络科技有限公司 一种基于图神经网络的流量监测软件的测试方法、装置及存储介质
CN117710118A (zh) * 2023-12-28 2024-03-15 中保金服(深圳)科技有限公司 一种智能理赔分析方法及***
CN117557256B (zh) * 2024-01-12 2024-03-22 西南交通大学 一种高速铁路轨道平顺性多任务协调维护方法
CN117579400B (zh) * 2024-01-17 2024-03-29 国网四川省电力公司电力科学研究院 一种基于神经网络的工控***网络安全监测方法及***
CN117786467B (zh) * 2024-02-28 2024-04-30 上海交通大学四川研究院 基于自适应加点的飞机着陆风险预测的分类模型构建方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016132468A1 (ja) 2015-02-18 2016-08-25 株式会社日立製作所 データ評価方法および装置、故障診断方法および装置
CN107655850A (zh) 2016-07-25 2018-02-02 上海创和亿电子科技发展有限公司 基于近红外光谱的非线性建模方法及***

Family Cites Families (50)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS55128591A (en) 1979-03-29 1980-10-04 Oouchi Shinko Kagaku Kogyo Kk Manufacture of phosphorothiolate
FR2698704B1 (fr) 1992-11-30 1995-03-10 Didier Heckmann Procédé et dispositif de maintenance prédictive.
EP0633536A1 (en) 1993-07-09 1995-01-11 International Business Machines Corporation Diagnosis method and system for predictive machine maintenance
US5539638A (en) 1993-08-05 1996-07-23 Pavilion Technologies, Inc. Virtual emissions monitor for automobile
US5566092A (en) 1993-12-30 1996-10-15 Caterpillar Inc. Machine fault diagnostics system and method
WO1997006418A1 (en) * 1995-08-07 1997-02-20 Boehringer Mannheim Corporation Biological fluid analysis using distance outlier detection
US6298308B1 (en) 1999-05-20 2001-10-02 Reid Asset Management Company Diagnostic network with automated proactive local experts
US7457785B1 (en) 2000-08-25 2008-11-25 Battelle Memorial Institute Method and apparatus to predict the remaining service life of an operating system
DE10161633A1 (de) 2000-12-14 2002-08-08 Siemens Corp Res Inc Verfahren und Vorrichtung zum Bereitstellen einer virtuellen Altersschätzung für die Vorhersage der verbleibenden Lebensdauer eines Systems unter Verwendung eines neuronalen Netzes
US7958063B2 (en) * 2004-11-11 2011-06-07 Trustees Of Columbia University In The City Of New York Methods and systems for identifying and localizing objects based on features of the objects that are mapped to a vector
KR101322434B1 (ko) 2005-07-11 2013-10-28 브룩스 오토메이션 인코퍼레이티드 지능형 상태 감시 및 결함 진단 시스템
US8275642B2 (en) 2007-01-19 2012-09-25 International Business Machines Corporation System to improve predictive maintenance and warranty cost/price estimation
US8676721B2 (en) 2009-09-18 2014-03-18 Apo Offshore, Inc. Method, system and apparatus for intelligent management of oil and gas platform surface equipment
US8781982B1 (en) 2011-09-23 2014-07-15 Lockheed Martin Corporation System and method for estimating remaining useful life
US11074495B2 (en) * 2013-02-28 2021-07-27 Z Advanced Computing, Inc. (Zac) System and method for extremely efficient image and pattern recognition and artificial intelligence platform
US9785983B2 (en) * 2012-06-13 2017-10-10 Opera Solutions U.S.A., Llc System and method for detecting billing errors using predictive modeling
JP2013257251A (ja) 2012-06-14 2013-12-26 Internatl Business Mach Corp <Ibm> 異常検知方法、プログラム及びシステム
US9263036B1 (en) * 2012-11-29 2016-02-16 Google Inc. System and method for speech recognition using deep recurrent neural networks
US9892238B2 (en) 2013-06-07 2018-02-13 Scientific Design Company, Inc. System and method for monitoring a process
CA2870080C (en) 2013-11-08 2017-12-19 Accenture Global Services Limited Network node failure predictive system
US10824958B2 (en) 2014-08-26 2020-11-03 Google Llc Localized learning from a global model
GB2532075A (en) * 2014-11-10 2016-05-11 Lego As System and method for toy recognition and detection based on convolutional neural networks
US10229356B1 (en) * 2014-12-23 2019-03-12 Amazon Technologies, Inc. Error tolerant neural network model compression
KR101677358B1 (ko) 2015-01-28 2016-11-18 주식회사 엑센솔루션 인공 신경망 알고리즘을 이용한 주조 설비에 대한 고장 예지 시스템 및 방법
CN104881684A (zh) 2015-05-27 2015-09-02 天津大学 一种立体图像质量客观评价方法
EP3314541A1 (en) * 2015-06-26 2018-05-02 Sentiance NV Deriving movement behaviour from sensor data
US20170091615A1 (en) * 2015-09-28 2017-03-30 Siemens Aktiengesellschaft System and method for predicting power plant operational parameters utilizing artificial neural network deep learning methodologies
US20160071004A1 (en) 2015-10-23 2016-03-10 Sarkhoon and Qeshm LLC Method and system for predictive maintenance of control valves
US10831577B2 (en) 2015-12-01 2020-11-10 Preferred Networks, Inc. Abnormality detection system, abnormality detection method, abnormality detection program, and method for generating learned model
EP3179326B1 (en) 2015-12-10 2024-04-10 KNORR-BREMSE Systeme für Schienenfahrzeuge GmbH Railway vehicle maintenance system with modular recurrent neural networks for performing time series prediction
US11062229B1 (en) 2016-02-18 2021-07-13 Deepmind Technologies Limited Training latent variable machine learning models using multi-sample objectives
US20180284735A1 (en) 2016-05-09 2018-10-04 StrongForce IoT Portfolio 2016, LLC Methods and systems for industrial internet of things data collection in a network sensitive upstream oil and gas environment
CN106168799B (zh) 2016-06-30 2019-05-03 常伟 一种基于大数据机器学习进行电动汽车电池预测性维护的方法
US11042811B2 (en) * 2016-10-05 2021-06-22 D-Wave Systems Inc. Discrete variational auto-encoder systems and methods for machine learning using adiabatic quantum computers
US11321609B2 (en) 2016-10-19 2022-05-03 Samsung Electronics Co., Ltd Method and apparatus for neural network quantization
US10303979B2 (en) 2016-11-16 2019-05-28 Phenomic Ai Inc. System and method for classifying and segmenting microscopy images with deep multiple instance learning
KR101988956B1 (ko) 2016-12-06 2019-06-12 한국과학기술원 외부환경 변동을 고려한 적응적 다변량 통계적 공정 모니터링 기반 구조물 건전성 평가 방법
US10140574B2 (en) * 2016-12-31 2018-11-27 Via Alliance Semiconductor Co., Ltd Neural network unit with segmentable array width rotator and re-shapeable weight memory to match segment width to provide common weights to multiple rotator segments
WO2018192672A1 (en) 2017-04-19 2018-10-25 Siemens Healthcare Gmbh Target detection in latent space
US10997499B1 (en) * 2017-05-12 2021-05-04 Veritas Technologies Llc Systems and methods for file system metadata analytics
CN109255440B (zh) 2017-07-11 2020-08-04 上海有孚网络股份有限公司 基于递归神经网络(rnn)的对电力生产设备进行预测性维护的方法
CN107844795B (zh) * 2017-11-18 2018-09-04 中国人民解放军陆军工程大学 基于主成分分析的卷积神经网络特征提取方法
US11069062B2 (en) * 2017-11-27 2021-07-20 Deciphex Automated screening of histopathology tissue samples via analysis of a normal model
CN108445752B (zh) 2018-03-02 2021-08-17 北京工业大学 一种自适应选择深度特征的随机权神经网络集成建模方法
US10789703B2 (en) * 2018-03-19 2020-09-29 Kla-Tencor Corporation Semi-supervised anomaly detection in scanning electron microscope images
CN108921303A (zh) 2018-05-29 2018-11-30 青岛鹏海软件有限公司 工业电动机的故障诊断及预测性维护方法
EP3620983B1 (en) * 2018-09-05 2023-10-25 Sartorius Stedim Data Analytics AB Computer-implemented method, computer program product and system for data analysis
US12001949B2 (en) * 2018-09-05 2024-06-04 Sartorius Stedim Data Analytics Ab Computer-implemented method, computer program product and system for data analysis
WO2021026400A1 (en) * 2019-08-06 2021-02-11 Neuroenhancement Lab, LLC System and method for communicating brain activity to an imaging device
EP3926420B1 (en) 2020-06-16 2023-11-01 Robert Bosch GmbH Making time-series predictions of a computer-controlled system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016132468A1 (ja) 2015-02-18 2016-08-25 株式会社日立製作所 データ評価方法および装置、故障診断方法および装置
CN107655850A (zh) 2016-07-25 2018-02-02 上海创和亿电子科技发展有限公司 基于近红外光谱的非线性建模方法及***

Also Published As

Publication number Publication date
CN112655004B (zh) 2024-03-26
JP7187681B2 (ja) 2022-12-12
EP3620986A1 (en) 2020-03-11
EP4290412A2 (en) 2023-12-13
CN112639834A (zh) 2021-04-09
WO2020049087A1 (en) 2020-03-12
US12001935B2 (en) 2024-06-04
EP3847586A1 (en) 2021-07-14
WO2020049094A1 (en) 2020-03-12
CN112655003A (zh) 2021-04-13
US20210350113A1 (en) 2021-11-11
CN112655004A (zh) 2021-04-13
EP3620983A1 (en) 2020-03-11
JP2022500745A (ja) 2022-01-04
US20210334656A1 (en) 2021-10-28
EP3620983B1 (en) 2023-10-25
JP2022500744A (ja) 2022-01-04
WO2020049098A1 (en) 2020-03-12
EP4290412A3 (en) 2024-01-03

Similar Documents

Publication Publication Date Title
JP7223839B2 (ja) 異常検出および/または予知保全のためのコンピュータ実装方法、コンピュータプログラム製品およびシステム
US20200285997A1 (en) Near real-time detection and classification of machine anomalies using machine learning and artificial intelligence
US11921566B2 (en) Abnormality detection system, abnormality detection method, abnormality detection program, and method for generating learned model
US12001949B2 (en) Computer-implemented method, computer program product and system for data analysis
JP6811276B2 (ja) 多次元時系列におけるスパース・ニューラル・ネットワーク・ベース異常検出
Maya et al. dLSTM: a new approach for anomaly detection using deep learning with delayed prediction
US11494690B2 (en) Apparatus and method of high dimensional data analysis in real-time
US11381583B1 (en) Systems and methods for detecting anomalous system or network behavior
Lin et al. A data-driven fault diagnosis method using modified health index and deep neural networks of a rolling bearing
WO2022115419A1 (en) Method of detecting an anomaly in a system
Mustafa et al. An adaptive metaheuristic optimization approach for Tennessee Eastman process for an industrial fault tolerant control system
CN115427986A (zh) 用于从大容量、高速流式数据动态生成预测分析的算法学习引擎
CN110874601B (zh) 识别设备运行状态的方法、状态识别模型训练方法及装置
Yang et al. Reducing response delay in multivariate process monitoring by a stacked long-short term memory network and real-time contrasts
Briceno-Mena et al. Data mining and knowledge discovery in chemical processes: Effect of alternative processing techniques
Colace et al. Unsupervised Learning Techniques for Vibration-Based Structural Health Monitoring Systems Driven by Data: A General Overview
Finotti et al. Novelty detection on a laboratory benchmark slender structure using an unsupervised deep learning algorithm
Le et al. Unsupervised detecting anomalies in multivariate time series by Robust Convolutional LSTM Encoder–Decoder (RCLED)
Kun et al. Deep Multiscale Soft‐Threshold Support Vector Data Description for Enhanced Heavy‐Duty Gas Turbine Generator Sets’ Anomaly Detection
Song et al. Advancements in Bearing Remaining Useful Life Prediction Methods: A Comprehensive Review
Lin et al. Industrial data-driven modeling for imbalanced fault diagnosis
Ripon et al. Design & Implementation of Automatic Machine Condition Monitoring and Maintenance System in Limited Resource Situations
Shirshahi et al. Identification of propagation path and root cause of faults based on generative adversarial networks in industrial systems
Weng et al. Dirichlet probability navigated fault detection via key-group memory auto-encoder under non-stationary working conditions
Wei et al. Anomaly Detection Based on Regularized Vector Auto Regression in Thermal Power Plant

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210712

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20220628

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220705

A524 Written submission of copy of amendment under article 19 pct

Free format text: JAPANESE INTERMEDIATE CODE: A524

Effective date: 20220929

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230206

R150 Certificate of patent or registration of utility model

Ref document number: 7223839

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150