JP7072531B2 - Anomaly detection device and anomaly detection method - Google Patents
Anomaly detection device and anomaly detection method Download PDFInfo
- Publication number
- JP7072531B2 JP7072531B2 JP2019044798A JP2019044798A JP7072531B2 JP 7072531 B2 JP7072531 B2 JP 7072531B2 JP 2019044798 A JP2019044798 A JP 2019044798A JP 2019044798 A JP2019044798 A JP 2019044798A JP 7072531 B2 JP7072531 B2 JP 7072531B2
- Authority
- JP
- Japan
- Prior art keywords
- time
- abnormality
- degree
- standard
- state
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
Description
本発明は、異常検出装置および異常検出方法に関する。 The present invention relates to an abnormality detection device and an abnormality detection method.
特許文献1には、「前記標準時系列データ算出ステップで算出された標準時系列データを一軸に配置し前記異常検出対象データ取得ステップで取得された対象時系列データを他軸に配置したDTW表記において該標準時系列データと該対象時系列データとの間の距離を算出する距離算出ステップと、前記距離算出ステップで算出された距離が該時系列データ中で予め規定された基準よりも大きいものを異常として検出する異常データ点検出ステップとを有する異常検出方法」が開示されている。 In Patent Document 1, "in a DTW notation in which the standard time series data calculated in the standard time series data calculation step is arranged on one axis and the target time series data acquired in the abnormality detection target data acquisition step is arranged on the other axis. An abnormality is defined as a distance calculation step for calculating the distance between the standard time series data and the target time series data, and a distance calculated in the distance calculation step larger than a predetermined reference in the time series data. An abnormality detection method including an abnormality data point detection step to be detected ”is disclosed.
上述の特許文献1は、時系列データの異常を検出するために、システムがDTW(Dynamic Time Warping:動的時間伸縮法)によって標準時系列データと対象時系列データを時間伸縮することで双方のデータ間の整合点を求め、距離を算出して合計することで異常の度合いを測り、異常か否かの判定を行っている。これにより、時間的なずれを度外視して比較することが可能となっている。その反面、所要時間異常(ある動作の遅延等)を検出するのが困難となっている。 In the above-mentioned Patent Document 1, in order to detect an abnormality in time-series data, the system expands and contracts the standard time-series data and the target time-series data by DTW (Dynamic Time Warping). The degree of abnormality is measured by finding the matching points between the two, calculating the distance, and summing up, and determining whether or not it is abnormal. This makes it possible to make comparisons without regard to time lag. On the other hand, it is difficult to detect an abnormality in the required time (delay of a certain operation, etc.).
本発明の目的は、時系列データ上で、遅延や前倒し等の所要時間異常を考慮した異常検出を可能とする技術を提供することを目的とする。 An object of the present invention is to provide a technique capable of detecting anomalies on time-series data in consideration of required time anomalies such as delay and advancement.
本願は、上記課題の少なくとも一部を解決する手段を複数含んでいるが、その例を挙げるならば、以下のとおりである。 The present application includes a plurality of means for solving at least a part of the above problems, and examples thereof are as follows.
本発明の一態様は、異常検出装置であって、基準となる時系列データを状態モデル化した標準時系列モデル数値パラメータと、各状態の持続時間をモデル化した標準遷移時間モデル数値パラメータと、を記憶するモデル記憶部と、異常検出の対象となる時系列データの入力を受け付ける時系列データ入力部と、上記異常検出の対象となる時系列データに対して、上記標準時系列モデル数値パラメータによる状態推定を行い、各データ点の状態変数を推定する状態推定部と、上記標準遷移時間モデル数値パラメータと、上記状態推定部が推定した上記各データ点の状態変数と、を用いて上記各データ点のうち所定の状態から別の状態への遷移に要する所要時間の異常度を算出する異常度計算部と、上記異常度を経過時間ごとの遅延度としてグラフ化して出力する異常度出力部と、を備える。 One aspect of the present invention is an abnormality detection device, which comprises a standard time-series model numerical parameter in which a reference time-series data is modeled, and a standard transition time model numerical parameter in which the duration of each state is modeled. State estimation by the standard time series model numerical parameters for the model storage unit to be stored, the time series data input unit that accepts the input of the time series data to be the target of abnormality detection, and the time series data to be the target of the abnormality detection. The state estimation unit that estimates the state variable of each data point, the standard transition time model numerical parameter, and the state variable of each data point estimated by the state estimation unit are used to determine each data point. Of these, an abnormality degree calculation unit that calculates the degree of abnormality of the time required for transition from a predetermined state to another state, and an abnormality degree output unit that outputs the above abnormality degree as a graph as a delay degree for each elapsed time. Be prepared.
本発明によれば、時系列データ上で、所要時間異常を考慮した異常検出を可能とする技術を提供することができる。 INDUSTRIAL APPLICABILITY According to the present invention, it is possible to provide a technique capable of detecting an abnormality in consideration of a required time abnormality on time-series data.
上記した以外の課題、構成及び効果は、以下の実施形態の説明により明らかにされる。 Issues, configurations and effects other than those described above will be clarified by the following description of the embodiments.
以下の実施形態においては便宜上その必要があるときは、複数のセクションまたは実施の形態に分割して説明するが、特に明示した場合を除き、それらはお互いに無関係なものではなく、一方は他方の一部または全部の変形例、詳細、補足説明等の関係にある。 In the following embodiments, when necessary for convenience, the description will be divided into a plurality of sections or embodiments, but unless otherwise specified, they are not unrelated to each other, and one is the other. It is related to some or all modifications, details, supplementary explanations, etc.
また、以下の実施形態において、要素の数等(個数、数値、量、範囲等を含む)に言及する場合、特に明示した場合および原理的に明らかに特定の数に限定される場合等を除き、その特定の数に限定されるものではなく、特定の数以上でも以下でもよい。 Further, in the following embodiments, when the number of elements (including the number, numerical value, quantity, range, etc.) is referred to, except when explicitly stated or when the number is clearly limited to a specific number in principle. , The number is not limited to the specific number, and may be more than or less than the specific number.
さらに、以下の実施形態において、その構成要素(要素ステップ等も含む)は、特に明示した場合および原理的に明らかに必須であると考えられる場合等を除き、必ずしも必須のものではないことは言うまでもない。 Furthermore, it goes without saying that, in the following embodiments, the components (including element steps and the like) are not necessarily essential unless otherwise specified or clearly considered to be essential in principle. stomach.
同様に、以下の実施形態において、構成要素等の形状、位置関係等に言及するときは特に明示した場合および原理的に明らかにそうではないと考えられる場合等を除き、実質的にその形状等に近似または類似するもの等を含むものとする。このことは、上記数値および範囲についても同様である。 Similarly, in the following embodiments, when the shape, positional relationship, etc. of the constituent elements are referred to, the shapes, etc. It shall include those that are close to or similar to. This also applies to the above numerical values and ranges.
また、実施形態を説明するための全図において、同一の部材には原則として同一の符号を付し、その繰り返しの説明は省略する。以下、本発明の各実施形態について図面を用いて説明する。 Further, in all the drawings for explaining the embodiment, the same members are in principle the same reference numerals, and the repeated description thereof will be omitted. Hereinafter, each embodiment of the present invention will be described with reference to the drawings.
DTW(Dynamic Time Warping)は、整合演算という上位概念で表すことができる。本実施形態では、時系列データの時間伸縮を行い、整合点を算出する手法を整合演算としている。整合演算は時間伸縮を伴い、標準時系列データと対象時系列データにおける変化の速度が異なる場合でも、整合点を求めることができる。 DTW (Dynamic Time Warping) can be expressed by a superordinate concept of matching operation. In the present embodiment, the method of calculating the matching point by expanding and contracting the time series data with time is defined as the matching operation. The matching operation involves time expansion and contraction, and even if the rate of change in the standard time series data and the target time series data is different, the matching point can be obtained.
一般に、標準時系列データと整合演算はそれぞれ、標準時系列モデルと状態推定という上位概念で表すことができる。ここでは、状態推定は時系列データを入力として数値パラメータを用いて各データ点に対応する状態変数を推定する演算方法であり、用いられる数値パラメータと演算方法を標準時系列モデルとしている。整合演算は標準時系列データの各データ点を異なる状態として扱い、与えられた対象時系列データの各データ点をそのどれかに整合させるという意味で、状態推定の一種として考えることができる。しかし、遷移時間に遅延等の異常が見られる場合にこれを異常として検出することはできない。 In general, standard time series data and matching operations can be represented by the superordinate concepts of standard time series model and state estimation, respectively. Here, the state estimation is a calculation method in which the state variables corresponding to each data point are estimated by using the time series data as an input and the numerical parameters, and the numerical parameters and the calculation method used are used as a standard time series model. The matching operation can be considered as a kind of state estimation in the sense that each data point of the standard time series data is treated as a different state and each data point of the given target time series data is matched with one of them. However, when an abnormality such as a delay is observed in the transition time, this cannot be detected as an abnormality.
そこで、本発明に係る実施形態においては、新たに、各状態の標準的な遷移時間を表現した標準遷移時間モデルを導入する。対象時系列データの異常度を算出する任意のデータ点Aからデータ点Bまでの区間について、対象時系列データに対する状態推定の結果から取得した対象区間に含まれる状態と、標準遷移時間モデルを参照して得た各状態の遷移時間分布から、区間についての遷移時間分布もしくはその代表値を算出する。対象時系列データにおけるデータ点Aからデータ点Bまでの区間遷移時間の区間遷移時間分布への不適合の度合い、もしくは遷移時間と区間遷移時間分布の代表値と差異を算出して異常度とする。 Therefore, in the embodiment of the present invention, a standard transition time model expressing the standard transition time of each state is newly introduced. For the section from any data point A to data point B for calculating the degree of abnormality of the target time series data, refer to the state included in the target section obtained from the result of state estimation for the target time series data and the standard transition time model. From the transition time distribution of each state obtained in the above, the transition time distribution for the section or its representative value is calculated. The degree of nonconformity of the interval transition time from the data point A to the data point B in the target time series data to the interval transition time distribution, or the difference between the transition time and the representative value of the interval transition time distribution is calculated and used as the degree of abnormality.
[第一の実施形態]図1は、第一の実施形態における所要時間異常検出部の例の概要を示す図である。ここでは、時系列データ入力部103によって入力を受け付ける時系列データは、データ点が得られた時点から経過時間順に複数連なっているものとする。データ点が得られた時点を表す変数をタイムスタンプ、データ点で測定された値を表す変数をデータ値とする。データ値とタイムスタンプをそれぞれ時間順に並べたものをデータ値系列、タイムスタンプ系列とし、この二つによって時系列データを保持・表現することが可能である。ただし、時系列データの表現において、タイムスタンプ系列を用いることは必須ではなく、データ点が一定の規則にしたがってサンプリングされている場合は、タイムスタンプ系列の代わりにサンプリング規則を保持しておけば後から算出可能である。また、データ値についてはスカラーだけでなくベクトルや画像なども含まれる。
[First Embodiment] FIG. 1 is a diagram showing an outline of an example of a required time abnormality detection unit in the first embodiment. Here, it is assumed that a plurality of time-series data for which input is received by the time-series
また、モデル記憶部102には事前に所定の単一のまたは複数の基準となる時系列データを状態モデル化した情報がプリセットされているものとする。プリセットの方法については、第二の実施形態に示すような方法により実現可能である。
Further, it is assumed that the
図2は、標準時系列モデル数値パラメータ記憶部の初期確率表の例を示す図である。初期確率表1021は、単一の時系列データにおける標準時系列モデルの初期状態の確率を示す表である。初期確率表1021には、状態1021Aと、初期確率1021Bと、が対応付けて含まれている。すなわち、状態ごとにその状態が初期状態である確率が格納されているといえる。
FIG. 2 is a diagram showing an example of an initial probability table of the standard time series model numerical parameter storage unit. The initial probability table 1021 is a table showing the probabilities of the initial state of the standard time series model in a single time series data. The initial probability table 1021 includes the
図3は、標準時系列モデル数値パラメータ記憶部の出力確率表の例を示す図である。出力確率表1022は、単一の時系列データにおける標準時系列モデルの出力確率を示す表である。出力確率表1022には、状態1022Aと、平均値1022Bと、分散1022Cと、が対応付けて含まれている。すなわち、状態ごとにその状態で取りうる値の分布が格納されているといえる。
FIG. 3 is a diagram showing an example of an output probability table of the standard time series model numerical parameter storage unit. The output probability table 1022 is a table showing the output probabilities of the standard time series model in a single time series data. The output probability table 1022 includes the
図4は、標準時系列モデル数値パラメータ記憶部の遷移確率表の例を示す図である。遷移確率表1023は、単一の時系列データにおける標準時系列モデルの状態間の遷移の確率を示す表である。遷移確率表1023には、元の状態1023Aと、遷移先の状態1023Bと、の組み合わせについて遷移する確率が対応付けて含まれている。すなわち、状態ごとにその状態が別の状態に遷移する確率が格納されているといえる。
FIG. 4 is a diagram showing an example of a transition probability table of the standard time series model numerical parameter storage unit. The transition probability table 1023 is a table showing the probability of transition between the states of the standard time series model in a single time series data. The transition probability table 1023 includes the probability of transition for the combination of the
図5は、標準遷移時間モデル数値パラメータ記憶部の遷移時間表の例を示す図である。遷移時間表1024は、単一の時系列データにおける標準時系列モデルのある状態ごとにその状態が持続する持続時間の分布をモデル化して示す表である。遷移時間表1024には、状態1024Aと、平均値1024Bと、分散1024Cと、が対応付けて含まれている。すなわち、状態ごとにその状態が持続する時間の取りうる値の分布が格納されているといえる。
FIG. 5 is a diagram showing an example of a transition time table of the standard transition time model numerical parameter storage unit. The transition time table 1024 is a table that models and shows the distribution of the duration of a standard time series model for each state in a single time series data. The transition time table 1024 includes the
所要時間異常検出部108は、時系列データ入力部103から時系列データを受け付けて、異常度出力部104へ異常度の情報を受け渡す。時系列データ入力部103には、異常検出の対象となる時系列データが入力される。所要時間異常検出部108は、時系列データをモデル記憶部102にプリセットされたモデル情報との整合演算を行い、異常度を特定して異常度出力部104へ受け渡す。異常度出力部104は、異常度の情報を用いて表示、あるいは記憶領域にその異常度の情報を記録する等の出力処理を行う。異常度出力部104は、このとき、特定の基準よりも異常度が大きい場合は警告を出力するようにしても良い。
The required time
所要時間異常検出部108には、異常度計算部100と、状態推定部101と、モデル記憶部102と、タイムスタンプ取得部105と、が含まれる。
The required time
タイムスタンプ取得部105は、時系列データ入力部103から受け付けた対象時系列データのタイムスタンプ系列を取得し、異常度計算部100へ送る。
The time
状態推定部101は、モデル記憶部102より読み出した標準時系列モデル数値パラメータ(図2の初期確率表1021と、図3の出力確率表1022と、図4の遷移確率表1023と、を含む)を用いて、時系列データ入力部103が受け付けた時系列データに対して、標準時系列モデルによる状態推定を行い、各データ点の状態変数を推定し、推定した状態変数について、異常度計算部100へ送る。なお、状態変数の推定には、状態推定部101は標準時系列モデルと対応した状態推定手法を用いる。例えば、標準時期系列モデルとして隠れマルコフモデルを用いる場合、状態推定部101は、ビタビアルゴリズムを用いて状態推定を行う。
The
また、状態推定部101は、複数の標準時系列モデル数値パラメータによる状態推定を行った結果を統合して状態変数を推定する。例えば、状態推定部101は、複数の推定された状態変数の平均値や中央値を状態変数として推定してもよいし、平均値と分散を求めて確率として推定するようにしてもよい。
Further, the
異常度計算部100は、モデル記憶部102から標準遷移時間モデル数値パラメータを読み出し、状態推定部101から送られてきた各データ点の状態変数と、タイムスタンプ取得部105から受け取ったタイムスタンプ系列と、を用いて異常度を算出し、異常度記憶表1001に記憶させる。
The anomaly
異常度計算部100は、対象時系列データの異常度を算出するデータ点Aからデータ点Bまでの対象区間について、データ点Aからデータ点Bに遷移するのに所要する遷移時間T1をタイムスタンプ系列から求め、対象時系列データに対する状態推定の結果から対象区間に含まれる状態を取得する。さらに、異常度計算部100は、標準遷移時間モデルを参照して得た各状態の遷移時間分布もしくはその代表値(平均値や最頻値など)から、区間全体にわたる遷移時間の分布である区間遷移時間分布もしくはその代表値T2を導出または近似させる。異常度計算部100は、遷移時間T1と代表値T2、または遷移時間T1と区間遷移時間分布の差異を、異常度として算出する。
The abnormality
区間遷移時間分布の近似方法として、異常度計算部100は、サンプリングによる方法を採用する。代表値の近似方法としては、異常度計算部100は、代表値が平均値であるとき累積和をとる方法を採用する。
As a method for approximating the interval transition time distribution, the anomaly
区間遷移時間分布と遷移時間T1の差異を算出する方法としては、遷移時間T1が区間遷移時間分布にどのくらい適合しているかを尤度等によって評価し、逆数や正負反転等を行ったものを異常度として用いる方法を採用する。 As a method of calculating the difference between the interval transition time distribution and the transition time T1, the degree to which the transition time T1 fits the interval transition time distribution is evaluated by the likelihood, etc., and the reciprocal or positive / negative inversion is abnormal. The method used as a degree is adopted.
また、遷移時間T1と代表値T2の差異を算出する方法としては、例えば、遷移時間T1を代表値T2で除す方法を採用する。この方法では、除した結果の商である異常度が「1」より大きくなった場合には、対象時系列データにおいてデータ点Aからデータ点Bへ遷移する間に遅れが生じたか、もしくは標準時系列データとは異なる事象がデータ点Aからデータ点Bの間に挿入された可能性が高いことを示す。つまり、作業動作が遅延したか、作業動作の間に別の作業動作が行われた可能性が高いことを示す。逆に、異常度が「1」より小さくなった場合には、対象時系列データにおいてデータ点Aからデータ点Bへの変化が早まったか、もしくは標準時系列データでは存在した事象がデータ点Aからデータ点Bの間で欠落した可能性が高いことを示す。つまり、作業動作が通常よりも前倒して行われたか、作業動作が省略された可能性が高いことを示す。 Further, as a method of calculating the difference between the transition time T1 and the representative value T2, for example, a method of dividing the transition time T1 by the representative value T2 is adopted. In this method, when the degree of abnormality, which is the quotient of the result of division, becomes larger than "1", is there a delay in the transition from data point A to data point B in the target time series data, or is there a delay in the transition? It indicates that an event different from the data is likely to be inserted between the data point A and the data point B. In other words, it is highly possible that the work operation was delayed or another work operation was performed during the work operation. On the contrary, when the degree of abnormality becomes smaller than "1", the change from the data point A to the data point B is accelerated in the target time series data, or the event existing in the standard time series data is data from the data point A. It indicates that there is a high possibility that the data is missing between the points B. In other words, it is highly possible that the work operation was performed earlier than usual, or that the work operation was omitted.
また、上記では任意の2点間(データ点Aとデータ点B間)について異常度(遅延度あるいは前倒し度)を算出する例を示したが、これに限られず、異常度計算部100は、対象時系列データ上の十分に短い時間間隔にある2点についての異常度計算を繰り返すことで、異常度を系列として算出することも可能である。
Further, in the above, an example of calculating the degree of abnormality (delay degree or advancement degree) between arbitrary two points (between data point A and data point B) is shown, but the present invention is not limited to this, and the abnormality
図6は、異常度計算部の異常度記憶表の例を示す図である。異常度記憶表1001には、時点を特定する情報である時点1001Aと、これに対応付けられた異常度1001Bと、が含まれる。
FIG. 6 is a diagram showing an example of an abnormality degree storage table of the abnormality degree calculation unit. The anomaly degree storage table 1001 includes a
また、異常度計算部100は、算出した異常度の系列に対して、平滑化などの事後処理を行うようにしてもよい。平滑化のアルゴリズムとしては、移動平均法等の公知の方法を採用する。
Further, the abnormality
図7は、平滑化処理の効果の例を示す図である。図7では、平滑化前を点線で、平滑化後を実線で表している。例えば時系列データにノイズ等の影響があった場合、異常度のばらつきが大きくなる。このような場合であっても、平滑化を行うことで、図7に示すようにばらつきを抑えた異常度の情報を得ることが可能となる。 FIG. 7 is a diagram showing an example of the effect of the smoothing process. In FIG. 7, the line before smoothing is represented by a dotted line, and the line after smoothing is represented by a solid line. For example, when the time-series data is affected by noise or the like, the variation in the degree of abnormality becomes large. Even in such a case, by performing smoothing, it is possible to obtain information on the degree of abnormality with suppressed variation as shown in FIG. 7.
異常度計算部100は、算出した異常度の系列またはこれを平滑化した情報を異常度記憶表1001に記憶させるとともに、異常度出力部104へ受け渡す。
The abnormality
図8は、異常度出力画面(ピークリスト)の例を示す図である。異常度出力画面(ピークリスト)1041は、異常度出力部104が異常度の系列またはこれを平滑化した情報をグラフ化して作成する画面情報である。異常度出力画面(ピークリスト)1041には、製品ごとに発生した遅延度(異常度)の大きさを示す遅延度グラフを表示する遅延度表示領域1041Aと、製品ごとに遅延度の大きい順に遅延度のピーク値を取った時点を表示するピークリスト1041Bと、画面表示を終える指示を受け付ける確認ボタン1041Cと、が含まれる。この画面により、利用者は製品毎の遅延度のピークを比較することが可能となる。なお、これに限られず、異常度出力部104は、以下のような出力画面を作成して出力するものであってもよい。
FIG. 8 is a diagram showing an example of an abnormality degree output screen (peak list). The abnormality degree output screen (peak list) 1041 is screen information created by the abnormality
図9は、異常度出力画面(遅延時点推定)の例を示す図である。異常度出力画面(遅延時点推定)1042は、異常度出力部104が異常度の系列またはこれを平滑化した情報をグラフ化して作成する画面情報である。異常度出力画面(遅延時点推定)1042には、製品に発生した遅延度(異常度)の大きさを示す遅延度グラフを表示する遅延度表示領域1042Aと、特に遅延度が大きく検出された時点(ピーク)に遅延発生を表示するコメント1042Bと、画面表示を終える指示を受け付ける確認ボタン1042Cと、が含まれる。この画面により、利用者は対象の製品について、大きな遅延が発生した時点を知り、作業内容を改善するヒントを得ることができる。
FIG. 9 is a diagram showing an example of an abnormality degree output screen (delay time point estimation). The abnormality degree output screen (delay time point estimation) 1042 is screen information created by the abnormality
図10は、異常度出力画面(遅延度合)の例を示す図である。異常度出力画面(遅延度合)1043は、異常度出力部104が異常度の系列またはこれを平滑化した情報をグラフ化して作成する画面情報である。異常度出力画面(遅延度合)1043には、製品に発生した遅延度(異常度)の大きさを示す遅延度グラフを表示する遅延度表示領域1043Aと、遅延度が検出された時間帯とその時間における遅延量を表示するコメント1043Bと、画面表示を終える指示を受け付ける確認ボタン1043Cと、が含まれる。この画面により、利用者は対象の製品について、連続的に遅延が発生した時点を知り、作業内容を改善するヒントを得ることができる。
FIG. 10 is a diagram showing an example of an abnormality degree output screen (delay degree). The abnormality degree output screen (delay degree) 1043 is screen information created by the abnormality
なお、異常度出力画面(遅延度合)1043を作成するためには、遅延度が検出された時間帯における遅延量を特定する必要がある。この特定処理では、異常度計算部100が遷移時間T1と代表値T2との差(T1-T2)を算出することにより遅延量を特定する。このようにすることで、遅延量を特定可能となり異常度出力画面(遅延度合)1043を作成することが可能となる。また、利用者も、グラフが異常度の量について作成され、グラフ上の幾何的情報(具体的には、遅延度の量を示す面積)により遅延度の量を把握することができる。
In order to create the abnormality degree output screen (delay degree) 1043, it is necessary to specify the delay amount in the time zone in which the delay degree is detected. In this specifying process, the abnormality
[第二の実施形態]図11は、第二の実施形態における所要時間異常検出部の例の概要を示す図である。第二の実施形態に係る所要時間異常検出部108´は、基本的に第一の実施形態にかかる所要時間異常検出部108と同様であるが、一部相違がある。以下、第一の実施形態にかかる所要時間異常検出部との相違を中心に説明する。
[Second Embodiment] FIG. 11 is a diagram showing an outline of an example of a required time abnormality detection unit in the second embodiment. The required time
第二の実施形態においては、モデル記憶部102に記憶される標準モデルを、時系列データ入力部103から得た単一あるいは複数の標準時系列データから学習可能としている。
In the second embodiment, the standard model stored in the
第二の実施形態にかかる所要時間異常検出部108´には、さらに、標準時系列モデル推定部106と、標準遷移時間モデル推定部107と、が含まれる。
The required time
時系列データ入力部103から入力された標準時系列データは、標準時系列モデル推定部106と、状態推定部101と、タイムスタンプ取得部105と、に受け渡される。タイムスタンプ取得部105は、タイムスタンプ系列を算出もしくは取得し、標準遷移時間モデル推定部107に受け渡す。
The standard time series data input from the time series
標準時系列モデル推定部106では、受け渡された単一または複数の標準時系列データを用いて、標準時系列モデルの数値パラメータを単一または複数推定してモデル記憶部102へ記憶させる。標準時系列モデルの学習は、隠れマルコフモデルや確率的オートマトンといった機械学習の枠組み等、AI(人工知能)等を用いて数値パラメータを推定するなどする。例えば隠れマルコフモデルを用いた場合、初期確率、遷移確率、出力確率が数値パラメータとして推定され、モデル記憶部102は、推定された数値パラメータ(初期確率表1021、出力確率表1022、遷移確率表1023)を標準時系列データごとに格納する。
The standard time-series
状態推定部101では、モデル記憶部102より読み出した標準時系列データの数値パラメータを用いて、時系列データ入力部103から入力された標準時系列データのセットに対して、標準時系列モデルを用いた状態推定を行い、時系列データの各データ点の状態変数を推定する。推定した状態変数について、状態推定部101は、標準遷移時間モデル推定部107に受け渡す。状態変数の推定には、標準時系列モデルと対応した状態推定手法を用いる。例えば隠れマルコフモデルであれば、ビタビアルゴリズムを用いて状態推定が可能である。また、複数の標準時系列データの数値パラメータを用いる場合、複数の状態が推定されるため、状態推定部101はその分布を求めるようにしてもよい。
The
標準遷移時間モデル推定部107は、状態推定部101から各データ点の状態変数を受け取り、タイムスタンプ取得部105から各点のタイムスタンプを受け取る。状態とタイムスタンプとを対応づけることで、標準遷移時間モデル推定部107は、対象の時系列データにおいて特定の状態が始まってから次の状態へ遷移するまでの遷移時間を算出する。
The standard transition time
そして、標準遷移時間モデル推定部107は、各データ点の状態変数とタイムスタンプとの複数の標準遷移時間モデルにおけるすべてのペアに対して遷移時間の算出を行うことで、各状態についての遷移時間の分布を求め、モデル記憶部102の標準遷移時間モデル(遷移時間表1024)として記憶させる。なお、遷移時間の分布については、標準遷移時間モデル推定部107は、ガウス分布等で近似させてそのパラメータ(平均と分散など)を代表値として記憶する。これに限られず、すべての算出した遷移時間をそのまま保存しても良い。
Then, the standard transition time
このようにすることで、モデル記憶部102に記憶される単一または複数の標準モデルを、時系列データ入力部103から得た標準時系列データを用いて学習することができる。
By doing so, it is possible to learn a single or a plurality of standard models stored in the
[第三の実施形態]図12は、第三の実施形態における所要時間異常検出部の例の概要を示す図である。第三の実施形態に係る所要時間異常検出部108´´は、基本的に第二の実施形態にかかる所要時間異常検出部108と同様であるが、一部相違がある。以下、第二の実施形態にかかる所要時間異常検出部との相違を中心に説明する。
[Third Embodiment] FIG. 12 is a diagram showing an outline of an example of a required time abnormality detection unit in the third embodiment. The required time
第三の実施形態においては、所要時間異常検出部108´´では、単一の(一連の)標準時系列データを受け取り標準モデルとすることで、標準時系列モデル推定部106および標準遷移時間モデル推定部107を不要とする例である。
In the third embodiment, the required time
第三の実施形態にかかる所要時間異常検出部108´´には、異常度計算部100と、状態推定部101と、モデル記憶部102と、タイムスタンプ取得部105と、が含まれる。
The required time
時系列データ入力部103は、単一の標準時系列データの入力を受け付ける。時系列データ入力部103は、データ値を時系列に並べたデータ値系列はモデル記憶部102に、タイムスタンプ系列はタイムスタンプ取得部105に、それぞれ送る。
The time series
タイムスタンプ取得部105は、入力された標準時系列データのタイムスタンプ系列を取得し、モデル記憶部102に受け渡す。
The time
そして、モデル記憶部102は、時系列データ入力部103から受け渡された標準時系列データのデータ値系列と、タイムスタンプ取得部105から受け渡されたタイムスタンプ系列と、を記憶する。これが、単一の標準モデルの学習の流れとなる。次に、異常度計算時の流れを説明する。
Then, the
時系列データ入力部103において異常検出の対象となる対象時系列データを受け付けると、時系列データ入力部103は、データ値系列は状態推定部101に、タイムスタンプ系列もしくはその算出に必要な情報はタイムスタンプ取得部105にそれぞれ送る。タイムスタンプ取得部105では、対象時系列データのタイムスタンプ系列を取得もしくは算出し、異常度計算部100へ受け渡す。
When the time-series
状態推定部101は、受け渡された対象時系列データのデータ値系列と、モデル記憶部102から読み出した標準時系列データのデータ値系列について整合演算を行い、整合演算の結果を異常度計算部100へ受け渡す。
The
図13は、整合演算の結果の例を示す図である。この結果は、第三の実施形態に限定されるものではなく、第一の実施形態から後述する第七の実施形態まで全ての実施形態において同様である。実線で示した対象時系列データのデータ値系列と、標準時系列データのデータ値系列について、整合する点を斜線で結んでいる。この結びつく点を特定するのが、整合演算であり、状態推定部101は、例えばDTW等の既知のアルゴリズムを採用して特定する。
FIG. 13 is a diagram showing an example of the result of the matching operation. This result is not limited to the third embodiment, and is the same in all embodiments from the first embodiment to the seventh embodiment described later. The matching points of the data value series of the target time series data shown by the solid line and the data value series of the standard time series data are connected by diagonal lines. It is the matching operation that specifies the points to be connected, and the
なお、異常度計算部100による異常度の算出に関しては、基本的に第一の実施形態および第二の実施形態と同様である。第三の実施形態の構成においては、標準時系列モデルの数値パラメータとして単一の標準時系列データを用いるため、第二の実施形態に示した複数の標準時系列データを用いてモデリングを行う標準時系列モデル推定部106が不要となる。また、標準時系列データの各データ点を異なる状態として扱い、整合演算によって対象時系列の各データ点を対応付けることで状態推定を実現しているため、各状態の標準遷移時間は単一の値であり、タイムスタンプ系列の隣り合ったタイムスタンプの差分と同一となる。そのため、標準遷移時間モデルとしてはタイムスタンプ系列をそのまま記憶すればよく、標準遷移時間モデル推定部107についても不要である。したがって、第三の実施形態に係る所要時間異常検出部108´´は、その構成を簡便なものとすることができる。
The calculation of the degree of abnormality by the abnormality
[第四の実施形態]図14は、第四の実施形態における所要時間異常検出システムの例を示す図である。第四の実施形態においては、ウェアラブルセンサ(位置の変化量を取得する着脱可能なセンサ)による人手作業の異常検出を行う作業異常検出システム205を対象として、所要時間異常検出部を適用する例である。
[Fourth Embodiment] FIG. 14 is a diagram showing an example of a required time abnormality detection system in the fourth embodiment. In the fourth embodiment, an example in which a required time abnormality detection unit is applied to a work
作業者201は、ウェアラブルセンサ202を身体に取り付けて、繰り返し作業を行っているものとする。ウェアラブルセンサ202にて計測した作業者201の動きのセンサデータを情報処理端末203へ無線通信を介して送る。動きを計測するウェアラブルセンサ202としては、加速度センサやジャイロセンサなどの位置の変化量を取得するセンサを用いることができる。情報処理端末203は、センサ受信部204でセンサデータを受け取ると、作業異常検出システム205へセンサデータを入力する。
It is assumed that the worker 201 attaches the wearable sensor 202 to the body and repeatedly performs the work. The sensor data of the movement of the worker 201 measured by the wearable sensor 202 is sent to the information processing terminal 203 via wireless communication. As the wearable sensor 202 for measuring the movement, a sensor such as an acceleration sensor or a gyro sensor that acquires the amount of change in position can be used. When the information processing terminal 203 receives the sensor data at the
図15は、第四の実施形態における作業異常検出システムの構成例を示す図である。センサ受信部204からセンサデータ(時系列データ)が特徴量抽出部206に受け渡され、特徴量抽出部206は、センサデータを異常検出に適した特徴量からなる時系列データに変換する。例えば、特徴量抽出部206は、センサデータに欠落がある場合には、補完処理を行い、センサデータに含まれるノイズをローパスフィルタ等のフィルタで除去することで、異常検出の精度向上が期待できる。また、特徴量抽出部206は、センサデータをフーリエ変換することで特徴量を得る。
FIG. 15 is a diagram showing a configuration example of the work abnormality detection system according to the fourth embodiment. Sensor data (time-series data) is passed from the
特徴量抽出部206により特徴量が抽出された時系列データは、系列区間切り出し部207へ受け渡され、切り出し制御部210からの制御信号を元に、繰り返し作業における一回の作業に対応する区間が切り出されて、所要時間異常検出部108へ受け渡される。
The time-series data from which the feature amount has been extracted by the feature
所要時間異常検出部108は、上述の第一、第二、第三の実施形態に示したように、受け取った時系列データに対して異常度を算出し、異常度出力部104へ算出結果を送る。異常度出力部104は、異常度を記録し、ディスプレイ等に表示する。このようにすることで、例えば、画面に表示された異常度を確認することで、正常に作業が行われているかを判断することが可能となる。特定の基準よりも異常度が大きい場合は警告を出力し、作業者や監督者などが対処を容易となるように支援するようにしても良い。
The required time
[第五の実施形態]図16は、第五の実施形態における所要時間異常検出システムの例を示す図である。第五の実施形態においては、カメラによる人手作業の異常検出を行う作業異常検出システム205を対象として、所要時間異常検出部を適用する例である。
[Fifth Embodiment] FIG. 16 is a diagram showing an example of a required time abnormality detection system in the fifth embodiment. In the fifth embodiment, the required time abnormality detection unit is applied to the work
作業者301は、繰り返し作業を行っているものとする。カメラ302にて撮影した作業者301の動画データを情報処理端末303へ有線通信または無線通信を介して送る。カメラ302としては、可視光を撮影するカメラに限られず、赤外線等の不可視光を撮影するものであってもよく、時系列に静止画を並べて動画とするものであれば足りる。 It is assumed that the worker 301 is performing the work repeatedly. The video data of the worker 301 taken by the camera 302 is sent to the information processing terminal 303 via wired communication or wireless communication. The camera 302 is not limited to a camera that captures visible light, and may capture invisible light such as infrared rays, and any camera that arranges still images in chronological order to form a moving image is sufficient.
情報処理端末303は、カメラ302から動画データを受け取ると、作業異常検出システム205へセンサデータを入力する。
When the information processing terminal 303 receives the moving image data from the camera 302, the information processing terminal 303 inputs the sensor data to the work
情報処理端末303は、動画データを、フレームデータが時間順に並んだ時系列データとみなすことができる。本実施形態に係る作業異常検出システム205は、基本的に第四の実施形態に係る作業異常検出システム205と同様であるが、一部相違がある。以下、第四の実施形態にかかる作業異常検出システム205との相違を中心に説明する。
The information processing terminal 303 can regard the moving image data as time-series data in which frame data are arranged in chronological order. The work
特徴量抽出部206は、動画データ内の被写体の関節位置を特徴量として抽出する。具体的には、特徴量抽出部206は、動画データを解析して骨格認識を行う。例えば、骨格認識アルゴリズムを用いて、特徴量抽出部206は、骨格を認識する。
The feature
図17は、第五の実施形態における骨格認識処理の例を示す図である。骨格認識処理においては、作業者301の映った動画データの各フレーム310に対して骨格認識を実施し、作業者の関節位置を抽出する(フレーム311)。関節位置を特徴量とすることで、作業者301の動きを反映した時系列データへ変換することができ、異常検出が容易となる。 FIG. 17 is a diagram showing an example of the skeleton recognition process in the fifth embodiment. In the skeleton recognition process, skeleton recognition is performed for each frame 310 of the moving image data of the worker 301, and the joint position of the worker is extracted (frame 311). By using the joint position as a feature amount, it is possible to convert it into time-series data that reflects the movement of the worker 301, and it becomes easy to detect an abnormality.
特徴量抽出された時系列データが系列区間切り出し部207へ受け渡された後の処理は、第四の実施形態と同様である。
The processing after the time-series data extracted by the feature amount is passed to the series
第五の実施形態によれば、例えば、画面に表示された異常度を確認することで、正常な速度で作業が行われているかを判断することが可能となる。特定の基準よりも異常度が大きい場合は警告を出力し、作業者や監督者などがカメラ302の撮影した映像で動作を確認しながら対処を容易となるように支援するようにしても良い。 According to the fifth embodiment, for example, by checking the degree of abnormality displayed on the screen, it is possible to determine whether the work is being performed at a normal speed. If the degree of abnormality is larger than a specific standard, a warning may be output to assist the operator, the supervisor, or the like to check the operation with the image taken by the camera 302 and facilitate the countermeasure.
[第六の実施形態]図18は、第六の実施形態における所要時間異常検出システムの例を示す図である。第六の実施形態においては、電流センサによる工作機械作業の異常検出を行う作業異常検出システム205を対象として、所要時間異常検出部を適用する例である。
[Sixth Embodiment] FIG. 18 is a diagram showing an example of a required time abnormality detection system in the sixth embodiment. In the sixth embodiment, the required time abnormality detection unit is applied to the work
工作機械501は、切削加工を行う機械であり、主軸電力線の電流値の大きさは、切削・アイドリング・停止という工作機械501の動作状態に応じたものである。電流センサ502は、計測した電流値を、情報処理端末503へ有線通信または無線通信を介して送る。
The machine tool 501 is a machine that performs cutting, and the magnitude of the current value of the spindle power line corresponds to the operating state of the machine tool 501 such as cutting, idling, and stopping. The
情報処理端末503は、電流センサ502から電流値の時系列データを受け取ると、作業異常検出システム205へセンサデータを入力する。
When the information processing terminal 503 receives the time-series data of the current value from the
本実施形態に係る作業異常検出システム205は、基本的に第四の実施形態に係る作業異常検出システム205と同様であるが、一部相違がある。以下、第四の実施形態にかかる作業異常検出システム205との相違を中心に説明する。
The work
特徴量抽出部206は、電流値の時系列データを対象として、センサデータを異常検出に適した特徴量からなる時系列データに変換する。例えば、センサデータに含まれるノイズをローパスフィルタ等のフィルタで除去することで、異常検出の精度向上が期待できる。例えば、特徴量抽出部206は、センサデータに欠落がある場合には、補完処理を行い、センサデータに含まれるノイズをローパスフィルタ等のフィルタで除去することで、異常検出の精度向上が期待できる。また、特徴量抽出部206は、センサデータをフーリエ変換することで特徴量を得る。
The feature
特徴量抽出された時系列データが系列区間切り出し部207へ受け渡された後の処理は、第四の実施形態と同様である。
The processing after the time-series data extracted by the feature amount is passed to the series
第六の実施形態によれば、例えば、画面に表示された異常度を確認することで、正常な速度で工作機械501が作業動作を行っているか否かを判断することが可能となる。特定の基準よりも異常度が大きい場合は警告を出力し、工作機械の修理や状態維持作業の手配等の対処を容易となるように支援するようにしても良い。 According to the sixth embodiment, for example, by checking the degree of abnormality displayed on the screen, it is possible to determine whether or not the machine tool 501 is performing the work operation at a normal speed. If the degree of abnormality is greater than a specific standard, a warning may be output to assist in facilitating measures such as repairing machine tools and arranging state maintenance work.
[第七の実施形態]図19は、第七の実施形態における異常度出力部の構成例を示す図である。第七の実施形態においては、異常度の出力を拡張するものである。第七の実施形態は、第一から第六の実施形態に示したいずれかの実施形態に係る所要時間異常検出システムと基本的に同様の構成を備えるものであるが、一部相違がある。以下、その相違について説明する。 [Seventh Embodiment] FIG. 19 is a diagram showing a configuration example of an abnormality degree output unit in the seventh embodiment. In the seventh embodiment, the output of the degree of anomaly is expanded. The seventh embodiment has basically the same configuration as the required time abnormality detection system according to any one of the first to sixth embodiments, but there are some differences. The difference will be described below.
図19に示すように、第七の実施形態においては、異常度出力部104は、標準時系列データ上で発生する事象の識別子と、該事象の発生時点および継続期間についての時間情報と、を保持した標準識別情報402を備え、標準識別情報402は、異常度表示部401と連結している。標準識別情報402は、例えば、第四の実施形態~第六の実施形態に示したような工場における作業の例では、一連の作業における作業の順番とそれぞれの標準作業時間が予め定められていることが多く、これに基づき予め情報が格納されている。
As shown in FIG. 19, in the seventh embodiment, the anomaly
所要時間異常検出部108から特定時間の異常度についての情報が渡されると、異常度表示部401は、標準識別情報402を参照し、標準時系列データ上において該当時間に発生する事象の識別子を検索し、異常度と識別子を対応付けて表示画面を作成する。
When the information about the abnormality degree of the specific time is passed from the required time
図20は、異常度出力画面(作業対照)の例を示す図である。異常度出力画面(作業対照)1044は、異常度出力部104が異常度の系列またはこれを平滑化した情報をグラフ化して作成する画面情報である。異常度出力画面(作業対照)1044には、製品に発生した遅延度(異常度)の大きさを示す遅延度グラフを表示する遅延度表示領域1044Aと、時点ごとの作業識別子を表示する識別子表示領域1044Bと、画面表示を終える指示を受け付ける確認ボタン1044Cと、が含まれる。この画面により、利用者は対象の製品について、大きな遅延が発生した時点の作業を知り、作業内容を改善するヒントを得ることができる。
FIG. 20 is a diagram showing an example of an abnormality degree output screen (work control). The abnormality degree output screen (work control) 1044 is screen information created by the abnormality
ここで、標準識別情報402が保持する識別子に対する時間情報は、あくまで標準時のものであり、対象時系列データにおいて時間異常が頻発した場合は、その誤差が無視できないほど大きい可能性がある。そこで、異常度表示部401は、異常度を元に対象時系列データにおける標準の場合との時間のずれを算出し、時間情報を補正した上で識別子を用いるようにしても良い。具体的には、異常度計算部100において、開始時刻から特定経過時刻までの区間を局所区間に区切り、各局所区間ABの異常度をT1およびT2の差分で求めた場合(T1、T2の定義は第一の実施形態と同様)、この累積値を算出することで、特定時間までの標準時に対する時間のずれを算出することができ、これを用いて時間情報を補正した上で表示に用いることができる。以上が、第七の実施形態である。
Here, the time information for the identifier held by the
なお、本発明は上記した実施形態に限定されるものではなく、様々な変形例が含まれる。例えば、上記した実施例は本発明を分かりやすく説明するために詳細に説明したものであり、必ずしも説明した全ての構成を備えるものに限定されるものではない。また、ある実施例の構成の一部を他の実施例の構成に置き換えることが可能であり、また、ある実施例の構成に他の実施例の構成を加えることも可能である。 The present invention is not limited to the above-described embodiment, and includes various modifications. For example, the above-described embodiment has been described in detail in order to explain the present invention in an easy-to-understand manner, and is not necessarily limited to the one including all the described configurations. Further, it is possible to replace a part of the configuration of one embodiment with the configuration of another embodiment, and it is also possible to add the configuration of another embodiment to the configuration of one embodiment.
また、各実施例の構成の一部について、他の構成の追加・削除・置換をすることが可能である。また、上記の各構成、機能、処理部、処理手段等は、それらの一部又は全部を、例えば集積回路で設計する等によりハードウェアで実現してもよい。また、上記の各構成、機能等は、プロセッサがそれぞれの機能を実現するプログラムを解釈し、実行することによりソフトウェアで実現してもよい。各機能を実現するプログラム、テーブル、ファイル等の情報は、メモリや、ハードディスク、SSD(Solid State Drive)等の記録装置、または、ICカード、SDカード、DVD(Digital Versatile Disc)等の記録媒体に置くことができる。 Further, it is possible to add / delete / replace a part of the configuration of each embodiment with another configuration. Further, each of the above configurations, functions, processing units, processing means and the like may be realized by hardware by designing a part or all of them by, for example, an integrated circuit. Further, each of the above configurations, functions, and the like may be realized by software by the processor interpreting and executing a program that realizes each function. Information such as programs, tables, and files that realize each function is stored in a memory, a hard disk, a recording device such as an SSD (Solid State Drive), or a recording medium such as an IC card, an SD card, or a DVD (Digital Versaille Disc). Can be placed.
また、制御線や情報線は説明上必要と考えられるものを示しており、製品上必ずしも全ての制御線や情報線を示しているとは限らない。実際には殆ど全ての構成が相互に接続されていると考えてもよい。 In addition, the control lines and information lines indicate what is considered necessary for explanation, and do not necessarily indicate all the control lines and information lines in the product. In practice, it can be considered that almost all configurations are interconnected.
本発明は、異常検出部、異常検出装置及び異常検出方法に限られず、異常検出システム、コンピュータ読み取り可能なプログラム、画像処理回路などの様々な態様で提供できる。 The present invention is not limited to the abnormality detection unit, the abnormality detection device, and the abnormality detection method, and can be provided in various forms such as an abnormality detection system, a computer-readable program, and an image processing circuit.
100・・・異常度計算部、101・・・状態推定部、102・・・モデル記憶部、103・・・時系列データ入力部、104・・・異常度出力部、105…タイムスタンプ取得部、108…所要時間異常検出部。 100 ... Abnormality calculation unit, 101 ... State estimation unit, 102 ... Model storage unit, 103 ... Time series data input unit, 104 ... Abnormality output unit, 105 ... Time stamp acquisition unit , 108 ... Required time abnormality detection unit.
Claims (13)
異常検出の対象となる時系列データの入力を受け付ける時系列データ入力部と、
前記異常検出の対象となる時系列データに対して、前記標準時系列モデル数値パラメータによる状態推定を行い、各データ点の状態変数を推定する状態推定部と、
前記標準遷移時間モデル数値パラメータと、前記状態推定部が推定した前記各データ点の状態変数と、を用いて前記各データ点のうち所定の状態から別の状態への遷移に要する所要時間の異常度を算出する異常度計算部と、
前記異常度を経過時間ごとの遅延度としてグラフ化して出力する異常度出力部と、
を備えることを特徴とする異常検出装置。 A model storage unit that stores standard time-series model numerical parameters that model the reference time-series data and standard transition time model numerical parameters that model the duration of each state.
A time-series data input unit that accepts input of time-series data that is the target of abnormality detection,
A state estimation unit that estimates the state of the time-series data to be detected by the standard time-series model numerical parameters and estimates the state variables of each data point.
Abnormality of time required for transition from a predetermined state to another state of each data point using the standard transition time model numerical parameter and the state variable of each data point estimated by the state estimation unit. Anomalous degree calculation unit that calculates the degree, and
An abnormality degree output unit that graphs and outputs the abnormality degree as a delay degree for each elapsed time, and
Anomaly detection device characterized by being equipped with.
前記異常検出の対象となる時系列データからタイムスタンプを取得するタイムスタンプ取得部を備え、
前記異常度計算部は、前記標準遷移時間モデル数値パラメータと、前記状態推定部が推定した前記各データ点の状態変数と、前記タイムスタンプと、を用いて各状態から別の状態への遷移に要する所要時間の異常度を算出する、
ことを特徴とする異常検出装置。 The abnormality detection device according to claim 1.
It is equipped with a time stamp acquisition unit that acquires a time stamp from the time series data that is the target of the abnormality detection.
The anomaly degree calculation unit makes a transition from each state to another state by using the standard transition time model numerical parameter, the state variable of each data point estimated by the state estimation unit, and the time stamp. Calculate the degree of abnormality of the required time,
Anomaly detection device characterized by this.
前記標準時系列モデル数値パラメータと、前記標準遷移時間モデル数値パラメータとは、前記モデル記憶部にそれぞれ複数記憶されており、
前記状態推定部は、複数の前記標準時系列モデル数値パラメータによる状態推定を行った結果を統合して前記状態変数を推定する、
ことを特徴とする異常検出装置。 The abnormality detection device according to claim 1.
A plurality of the standard time series model numerical parameters and the standard transition time model numerical parameters are stored in the model storage unit.
The state estimation unit estimates the state variables by integrating the results of state estimation using a plurality of standard time series model numerical parameters.
Anomaly detection device characterized by this.
受け渡された単一または複数の標準時系列データを用いて、標準時系列モデルの数値パラメータを単一または複数推定して前記標準時系列モデル数値パラメータとして前記モデル記憶部へ記憶させる標準時系列モデル推定部、
を備えることを特徴とする異常検出装置。 The abnormality detection device according to claim 1.
A standard time series model estimation unit that estimates single or multiple numerical parameters of a standard time series model using the passed single or multiple standard time series data and stores them in the model storage unit as the standard time series model numerical parameters.
Anomaly detection device characterized by being equipped with.
受け渡された単一または複数の標準時系列データにおいて特定の状態が始まってから次の状態へ遷移するまでの遷移時間を算出して前記標準遷移時間モデル数値パラメータとして前記モデル記憶部へ記憶させる標準遷移時間モデル推定部、
を備えることを特徴とする異常検出装置。 The abnormality detection device according to claim 1.
A standard that calculates the transition time from the start of a specific state to the transition to the next state in the passed single or multiple standard time series data and stores it in the model storage unit as the standard transition time model numerical parameter. Transition time model estimation unit,
Anomaly detection device characterized by being equipped with.
前記時系列データ入力部は、位置の変化量を取得するセンサから出力される時系列データの入力を受け付ける、
ことを特徴とする異常検出装置。 The abnormality detection device according to claim 1.
The time-series data input unit receives the input of time-series data output from the sensor that acquires the amount of change in position.
Anomaly detection device characterized by this.
動画データを解析して骨格認識を行う特徴量抽出部を備え、
前記時系列データ入力部は、カメラから出力される動画データの入力を受け付け、
前記状態推定部は、前記特徴量抽出部により前記動画データを解析した結果得られる骨格の関節位置の変化を時系列データとして状態を推定する、
ことを特徴とする異常検出装置。 The abnormality detection device according to claim 1.
Equipped with a feature amount extraction unit that analyzes video data and recognizes the skeleton.
The time-series data input unit receives the input of video data output from the camera and receives the input.
The state estimation unit estimates the state using the change in the joint position of the skeleton obtained as a result of analyzing the moving image data by the feature amount extraction unit as time series data.
Anomaly detection device characterized by this.
前記時系列データ入力部は、電流値を取得するセンサから出力される時系列データの入力を受け付ける、
ことを特徴とする異常検出装置。 The abnormality detection device according to claim 1.
The time-series data input unit receives an input of time-series data output from a sensor that acquires a current value.
Anomaly detection device characterized by this.
前記異常度出力部は、前記基準となる時系列データ上で発生する事象の識別子と、該事象の発生時点および継続期間についての時間情報と、を含む標準識別情報を予め記憶しており、
前記グラフ化において前記所定の状態に対応した前記事象の識別子を対応付ける、
ことを特徴とする異常検出装置。 The abnormality detection device according to claim 1.
The abnormality degree output unit stores in advance standard identification information including an identifier of an event that occurs on the reference time-series data and time information about the time and duration of the event.
In the graphing, the identifier of the event corresponding to the predetermined state is associated.
Anomaly detection device characterized by this.
前記異常度出力部は、前記グラフに対応付けて、前記遅延度のピーク値が大きくなる経過時間を示す、
ことを特徴とする異常検出装置。 The abnormality detection device according to claim 1.
The abnormality degree output unit indicates an elapsed time in which the peak value of the delay degree becomes large in association with the graph.
Anomaly detection device characterized by this.
前記異常度出力部は、前記グラフの遅延度のピーク値が大きくなる経過時間の位置に、遅延発生のコメントを示す、
ことを特徴とする異常検出装置。 The abnormality detection device according to claim 1.
The abnormality degree output unit indicates a comment of the occurrence of the delay at the position of the elapsed time when the peak value of the delay degree of the graph becomes large.
Anomaly detection device characterized by this.
前記異常度出力部は、前記グラフを前記異常度の量について作成し、前記グラフ上の幾何的情報により遅延度の量を示す、
ことを特徴とする異常検出装置。 The abnormality detection device according to claim 1.
The anomaly degree output unit creates the graph for the amount of the anomaly degree, and indicates the amount of the delay degree by the geometric information on the graph.
Anomaly detection device characterized by this.
前記情報処理装置は、
基準となる時系列データを状態モデル化した標準時系列モデル数値パラメータと、各状態の持続時間をモデル化した標準遷移時間モデル数値パラメータと、を記憶するモデル記憶部と、制御部と、を備え、
前記制御部は、
異常検出の対象となる時系列データの入力を受け付ける時系列データ入力ステップと、
前記異常検出の対象となる時系列データに対して、前記標準時系列モデル数値パラメータによる状態推定を行い、各データ点の状態変数を推定する状態推定ステップと、
前記標準遷移時間モデル数値パラメータと、前記状態推定ステップで推定した前記各データ点の状態変数と、を用いて前記各データ点のうち所定の状態から別の状態への遷移に要する所要時間の異常度を算出する異常度計算ステップと、
前記異常度を経過時間ごとの遅延度としてグラフ化して出力する異常度出力ステップと、
を実施することを特徴とする異常検出方法。 It is an abnormality detection method using an information processing device.
The information processing device is
It is equipped with a model storage unit and a control unit that store standard time-series model numerical parameters that model the reference time-series data and standard transition time model numerical parameters that model the duration of each state.
The control unit
A time-series data input step that accepts input of time-series data that is the target of anomaly detection,
A state estimation step for estimating the state variable of each data point by performing state estimation with the standard time series model numerical parameters for the time series data to be the target of the abnormality detection, and
Abnormality of time required for transition from a predetermined state to another state of each data point using the standard transition time model numerical parameter and the state variable of each data point estimated in the state estimation step. Abnormality calculation step to calculate the degree and
An abnormality degree output step that graphs and outputs the abnormality degree as a delay degree for each elapsed time, and
Anomaly detection method characterized by carrying out.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019044798A JP7072531B2 (en) | 2019-03-12 | 2019-03-12 | Anomaly detection device and anomaly detection method |
CN201911211451.9A CN111693794B (en) | 2019-03-12 | 2019-11-29 | Abnormality detection device and abnormality detection method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019044798A JP7072531B2 (en) | 2019-03-12 | 2019-03-12 | Anomaly detection device and anomaly detection method |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2020149208A JP2020149208A (en) | 2020-09-17 |
JP7072531B2 true JP7072531B2 (en) | 2022-05-20 |
Family
ID=72429967
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2019044798A Active JP7072531B2 (en) | 2019-03-12 | 2019-03-12 | Anomaly detection device and anomaly detection method |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP7072531B2 (en) |
CN (1) | CN111693794B (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113296990B (en) * | 2020-09-30 | 2022-06-24 | 阿里云计算有限公司 | Method and device for recognizing abnormity of time sequence data |
JP7162312B2 (en) * | 2021-03-04 | 2022-10-28 | 株式会社ハイテックシステム | Sequence execution status notification device and program |
CN113645096A (en) * | 2021-08-11 | 2021-11-12 | 四川华腾国盛科技有限公司 | Building intelligent engineering detection system |
CN116610482B (en) * | 2023-07-18 | 2023-10-17 | 山东理工大学 | Intelligent monitoring method for operation state of electrical equipment |
CN116738151B (en) * | 2023-08-09 | 2023-11-17 | 广东电网有限责任公司广州供电局 | Enterprise energy consumption carbon emission monitoring accounting method and system based on big data |
CN116826977B (en) * | 2023-08-28 | 2023-11-21 | 青岛恒源高新电气有限公司 | Intelligent management system for optical storage direct-flexible micro-grid |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007257184A (en) | 2006-03-22 | 2007-10-04 | Nec Corp | Obstacle factor estimation system, method and program |
JP2018163574A (en) | 2017-03-27 | 2018-10-18 | サクサ株式会社 | Log management device and program for log management |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09198420A (en) * | 1996-01-22 | 1997-07-31 | Toshiba Corp | Method and device for predicting delay time |
JPH10103308A (en) * | 1996-09-24 | 1998-04-21 | Mitsubishi Heavy Ind Ltd | Opening controlling method for air actuating valve and system thereof |
JP3996428B2 (en) * | 2001-12-25 | 2007-10-24 | 松下電器産業株式会社 | Abnormality detection device and abnormality detection system |
EP1676992A4 (en) * | 2003-09-24 | 2014-12-10 | A & D Co Ltd | Multi-signal analysis device |
KR101137535B1 (en) * | 2006-05-26 | 2012-04-20 | 가부시키가이샤 어드밴티스트 | Testing apparatus and test module |
KR100925370B1 (en) * | 2007-12-21 | 2009-11-09 | 주식회사 하이닉스반도체 | Data input device |
JP5140476B2 (en) * | 2008-03-31 | 2013-02-06 | アズビル株式会社 | Abnormality detection method and abnormality detection device |
JP5301310B2 (en) * | 2009-02-17 | 2013-09-25 | 株式会社日立製作所 | Anomaly detection method and anomaly detection system |
CN101865974B (en) * | 2009-04-20 | 2013-06-12 | 普诚科技股份有限公司 | Logic tester and method for simultaneously measuring delay time of a plurality of devices to be measured |
CN105652120A (en) * | 2015-12-31 | 2016-06-08 | 国网重庆潼南区供电有限责任公司 | Power transformer fault detection method and detection system |
CN107133343B (en) * | 2017-05-19 | 2018-04-13 | 哈工大大数据产业有限公司 | Big data abnormal state detection method and device based on time series approximate match |
CN109032829B (en) * | 2018-07-23 | 2020-12-08 | 腾讯科技(深圳)有限公司 | Data anomaly detection method and device, computer equipment and storage medium |
-
2019
- 2019-03-12 JP JP2019044798A patent/JP7072531B2/en active Active
- 2019-11-29 CN CN201911211451.9A patent/CN111693794B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007257184A (en) | 2006-03-22 | 2007-10-04 | Nec Corp | Obstacle factor estimation system, method and program |
JP2018163574A (en) | 2017-03-27 | 2018-10-18 | サクサ株式会社 | Log management device and program for log management |
Also Published As
Publication number | Publication date |
---|---|
JP2020149208A (en) | 2020-09-17 |
CN111693794B (en) | 2022-08-30 |
CN111693794A (en) | 2020-09-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP7072531B2 (en) | Anomaly detection device and anomaly detection method | |
JP5045770B2 (en) | Process analysis system | |
US20160231738A1 (en) | Information processing apparatus and analysis method | |
JP6672712B2 (en) | Abnormal work detection system and abnormal work detection method | |
JP7254546B2 (en) | Information processing device, information processing method and program | |
KR20200110404A (en) | Work motion analysis system and work motion analysis method | |
CN110480678A (en) | A kind of industrial robot collision checking method | |
JP2009205282A (en) | Motion analysis method, motion analysis device, and motion evaluation device using the analysis device | |
CN114365161A (en) | Work element analysis device and work element analysis method | |
EP3336746A1 (en) | System and method of video content filtering | |
JP7139987B2 (en) | Process information acquisition system, process information acquisition method, and process information acquisition program | |
US10996235B2 (en) | System and method for cycle duration measurement in repeated activity sequences | |
US20220404792A1 (en) | Method, system and computer program product for monitoring a shaping process | |
JP2020140365A (en) | Product quality defect prediction system | |
JP6247777B2 (en) | Abnormality diagnosis apparatus and abnormality diagnosis method | |
CN114648809A (en) | Computer implemented process monitoring method, apparatus, system and recording medium | |
JP2004078338A (en) | Method and system for evaluating computer performance | |
CN117296067A (en) | Machine learning device, classifying device, and control device | |
JP7449248B2 (en) | Support devices, support methods and support programs | |
JP7368189B2 (en) | Analysis equipment | |
KR100324654B1 (en) | Method for Measuring of Operation Time | |
JP2020008932A (en) | Image processing device, image processing method, and program | |
CN109844739B (en) | Method for pattern recognition in multiple signals | |
JP7311820B1 (en) | Abnormality determination method, abnormality determination device, and program | |
JP7235094B2 (en) | Abnormal work detection system, abnormal work detection method, and abnormal work detection program |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20210419 |
|
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: 20220426 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20220510 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 7072531 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |