JP7072531B2 - Anomaly detection device and anomaly detection method - Google Patents

Anomaly detection device and anomaly detection method Download PDF

Info

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
Application number
JP2019044798A
Other languages
Japanese (ja)
Other versions
JP2020149208A (en
Inventor
修一 西納
徹 矢崎
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
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 Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP2019044798A priority Critical patent/JP7072531B2/en
Priority to CN201911211451.9A priority patent/CN111693794B/en
Publication of JP2020149208A publication Critical patent/JP2020149208A/en
Application granted granted Critical
Publication of JP7072531B2 publication Critical patent/JP7072531B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements 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.

特開2009-245228号公報Japanese Unexamined Patent Publication No. 2009-245228

上述の特許文献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.

第一の実施形態における所要時間異常検出部の例の概要を示す図である。It is a figure which shows the outline of the example of the required time abnormality detection part in 1st Embodiment. 標準時系列モデル数値パラメータ記憶部の初期確率表の例を示す図である。It is a figure which shows the example of the initial probability table of the standard time series model numerical parameter storage part. 標準時系列モデル数値パラメータ記憶部の出力確率表の例を示す図である。It is a figure which shows the example of the output probability table of the standard time series model numerical parameter storage part. 標準時系列モデル数値パラメータ記憶部の遷移確率表の例を示す図である。It is a figure which shows the example of the transition probability table of the standard time series model numerical parameter storage part. 標準遷移時間モデル数値パラメータ記憶部の遷移時間表の例を示す図である。It is a figure which shows the example of the transition time table of the standard transition time model numerical parameter storage part. 異常度計算部の異常度記憶表の例を示す図である。It is a figure which shows the example of the abnormality degree memory table of the abnormality degree calculation part. 平滑化処理の効果の例を示す図である。It is a figure which shows the example of the effect of a smoothing process. 異常度出力画面(ピークリスト)の例を示す図である。It is a figure which shows the example of the abnormality degree output screen (peak list). 異常度出力画面(遅延時点推定)の例を示す図である。It is a figure which shows the example of the abnormality degree output screen (delay time point estimation). 異常度出力画面(遅延度合)の例を示す図である。It is a figure which shows the example of the abnormality degree output screen (delay degree degree). 第二の実施形態における所要時間異常検出部の例の概要を示す図である。It is a figure which shows the outline of the example of the required time abnormality detection part in the 2nd Embodiment. 第三の実施形態における所要時間異常検出部の例の概要を示す図である。It is a figure which shows the outline of the example of the required time abnormality detection part in 3rd Embodiment. 整合演算の結果の例を示す図である。It is a figure which shows the example of the result of the matching operation. 第四の実施形態における所要時間異常検出システムの例を示す図である。It is a figure which shows the example of the required time abnormality detection system in 4th Embodiment. 第四の実施形態における作業異常検出システムの構成例を示す図である。It is a figure which shows the configuration example of the work abnormality detection system in 4th Embodiment. 第五の実施形態における所要時間異常検出システムの例を示す図である。It is a figure which shows the example of the required time abnormality detection system in the 5th Embodiment. 第五の実施形態における骨格認識処理の例を示す図である。It is a figure which shows the example of the skeleton recognition processing in the 5th Embodiment. 第六の実施形態における所要時間異常検出システムの例を示す図である。It is a figure which shows the example of the required time abnormality detection system in the 6th Embodiment. 第七の実施形態における異常度出力部の構成例を示す図である。It is a figure which shows the structural example of the abnormality degree output part in 7th Embodiment. 異常度出力画面(作業対照)の例を示す図である。It is a figure which shows the example of the abnormality degree output screen (work control).

以下の実施形態においては便宜上その必要があるときは、複数のセクションまたは実施の形態に分割して説明するが、特に明示した場合を除き、それらはお互いに無関係なものではなく、一方は他方の一部または全部の変形例、詳細、補足説明等の関係にある。 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 data input unit 103 are connected in order of elapsed time from the time when the data points are obtained. The variable representing the time when the data point is obtained is defined as the time stamp, and the variable representing the value measured at the data point is defined as the data value. Data values and time stamps are arranged in chronological order as a data value series and a time stamp series, and it is possible to retain and express time series data by these two. However, it is not essential to use a time stamp series in the representation of time series data, and if the data points are sampled according to a certain rule, it is possible to keep the sampling rule instead of the time stamp series. It can be calculated from. The data values include not only scalars but also vectors and images.

また、モデル記憶部102には事前に所定の単一のまたは複数の基準となる時系列データを状態モデル化した情報がプリセットされているものとする。プリセットの方法については、第二の実施形態に示すような方法により実現可能である。 Further, it is assumed that the model storage unit 102 is preset with information obtained by modeling a state of a predetermined single or a plurality of reference time-series data in advance. The preset method can be realized by the method shown in the second embodiment.

図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 state 1021A and the initial probability 1021B in association with each other. That is, it can be said that the probability that the state is the initial state is stored for each state.

図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 state 1022A, the mean value 1022B, and the variance 1022C in association with each other. That is, it can be said that the distribution of values that can be taken in each state is stored for each state.

図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 original state 1023A and the transition destination state 1023B in association with each other. That is, it can be said that the probability that the state transitions to another state is stored for each state.

図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 state 1024A, the mean value 1024B, and the variance 1024C in association with each other. That is, it can be said that the distribution of possible values for the duration of the state is stored for each state.

所要時間異常検出部108は、時系列データ入力部103から時系列データを受け付けて、異常度出力部104へ異常度の情報を受け渡す。時系列データ入力部103には、異常検出の対象となる時系列データが入力される。所要時間異常検出部108は、時系列データをモデル記憶部102にプリセットされたモデル情報との整合演算を行い、異常度を特定して異常度出力部104へ受け渡す。異常度出力部104は、異常度の情報を用いて表示、あるいは記憶領域にその異常度の情報を記録する等の出力処理を行う。異常度出力部104は、このとき、特定の基準よりも異常度が大きい場合は警告を出力するようにしても良い。 The required time abnormality detection unit 108 receives the time-series data from the time-series data input unit 103, and passes the abnormality degree information to the abnormality degree output unit 104. Time-series data to be detected as an abnormality is input to the time-series data input unit 103. The required time abnormality detection unit 108 performs a matching operation with the model information preset in the model storage unit 102 for the time series data, identifies the abnormality degree, and passes it to the abnormality degree output unit 104. The abnormality degree output unit 104 performs output processing such as displaying using the abnormality degree information or recording the abnormality degree information in a storage area. At this time, the abnormality degree output unit 104 may output a warning when the abnormality degree is larger than a specific reference.

所要時間異常検出部108には、異常度計算部100と、状態推定部101と、モデル記憶部102と、タイムスタンプ取得部105と、が含まれる。 The required time abnormality detection unit 108 includes an abnormality degree calculation unit 100, a state estimation unit 101, a model storage unit 102, and a time stamp acquisition unit 105.

タイムスタンプ取得部105は、時系列データ入力部103から受け付けた対象時系列データのタイムスタンプ系列を取得し、異常度計算部100へ送る。 The time stamp acquisition unit 105 acquires the time stamp series of the target time series data received from the time series data input unit 103 and sends it to the abnormality degree calculation unit 100.

状態推定部101は、モデル記憶部102より読み出した標準時系列モデル数値パラメータ(図2の初期確率表1021と、図3の出力確率表1022と、図4の遷移確率表1023と、を含む)を用いて、時系列データ入力部103が受け付けた時系列データに対して、標準時系列モデルによる状態推定を行い、各データ点の状態変数を推定し、推定した状態変数について、異常度計算部100へ送る。なお、状態変数の推定には、状態推定部101は標準時系列モデルと対応した状態推定手法を用いる。例えば、標準時期系列モデルとして隠れマルコフモデルを用いる場合、状態推定部101は、ビタビアルゴリズムを用いて状態推定を行う。 The state estimation unit 101 contains standard time-series model numerical parameters read from the model storage unit 102 (including the initial probability table 1021 in FIG. 2, the output probability table 1022 in FIG. 3, and the transition probability table 1023 in FIG. 4). For the time-series data received by the time-series data input unit 103, state estimation is performed by a standard time-series model, state variables of each data point are estimated, and the estimated state variables are sent to the abnormality degree calculation unit 100. send. The state estimation unit 101 uses a state estimation method corresponding to the standard time series model to estimate the state variables. For example, when a hidden Markov model is used as a standard time series model, the state estimation unit 101 performs state estimation using a Viterbi algorithm.

また、状態推定部101は、複数の標準時系列モデル数値パラメータによる状態推定を行った結果を統合して状態変数を推定する。例えば、状態推定部101は、複数の推定された状態変数の平均値や中央値を状態変数として推定してもよいし、平均値と分散を求めて確率として推定するようにしてもよい。 Further, the state estimation unit 101 integrates the results of state estimation using a plurality of standard time series model numerical parameters to estimate state variables. For example, the state estimation unit 101 may estimate the mean value or the median value of a plurality of estimated state variables as state variables, or may obtain the mean value and the variance and estimate them as probabilities.

異常度計算部100は、モデル記憶部102から標準遷移時間モデル数値パラメータを読み出し、状態推定部101から送られてきた各データ点の状態変数と、タイムスタンプ取得部105から受け取ったタイムスタンプ系列と、を用いて異常度を算出し、異常度記憶表1001に記憶させる。 The anomaly degree calculation unit 100 reads out the standard transition time model numerical parameters from the model storage unit 102, the state variables of each data point sent from the state estimation unit 101, and the time stamp series received from the time stamp acquisition unit 105. , Is calculated and stored in the abnormality degree storage table 1001.

異常度計算部100は、対象時系列データの異常度を算出するデータ点Aからデータ点Bまでの対象区間について、データ点Aからデータ点Bに遷移するのに所要する遷移時間T1をタイムスタンプ系列から求め、対象時系列データに対する状態推定の結果から対象区間に含まれる状態を取得する。さらに、異常度計算部100は、標準遷移時間モデルを参照して得た各状態の遷移時間分布もしくはその代表値(平均値や最頻値など)から、区間全体にわたる遷移時間の分布である区間遷移時間分布もしくはその代表値T2を導出または近似させる。異常度計算部100は、遷移時間T1と代表値T2、または遷移時間T1と区間遷移時間分布の差異を、異常度として算出する。 The abnormality degree calculation unit 100 time stamps the transition time T1 required for transitioning from the data point A to the data point B for the target section from the data point A to the data point B for calculating the abnormality degree of the target time series data. Obtain from the series and acquire the state included in the target section from the result of state estimation for the target time series data. Further, the anomaly degree calculation unit 100 is a section which is a transition time distribution over the entire section from the transition time distribution of each state obtained by referring to the standard transition time model or its representative value (mean value, mode value, etc.). Derivation or approximation of the transition time distribution or its representative value T2. The abnormality degree calculation unit 100 calculates the difference between the transition time T1 and the representative value T2, or the transition time T1 and the interval transition time distribution as the abnormality degree.

区間遷移時間分布の近似方法として、異常度計算部100は、サンプリングによる方法を採用する。代表値の近似方法としては、異常度計算部100は、代表値が平均値であるとき累積和をとる方法を採用する。 As a method for approximating the interval transition time distribution, the anomaly degree calculation unit 100 adopts a sampling method. As a method for approximating the representative value, the abnormality degree calculation unit 100 adopts a method of taking a cumulative sum when the representative value is an average value.

区間遷移時間分布と遷移時間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 degree calculation unit 100 may perform the abnormality degree calculation unit 100. It is also possible to calculate the degree of abnormality as a series by repeating the calculation of the degree of abnormality for two points at sufficiently short time intervals on the target time series data.

図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 time point 1001A, which is information for specifying a time point, and an abnormality degree 1001B associated with the time point 1001A.

また、異常度計算部100は、算出した異常度の系列に対して、平滑化などの事後処理を行うようにしてもよい。平滑化のアルゴリズムとしては、移動平均法等の公知の方法を採用する。 Further, the abnormality degree calculation unit 100 may perform post-processing such as smoothing on the calculated abnormality degree series. As the smoothing algorithm, a known method such as a moving average method is adopted.

図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 degree calculation unit 100 stores the calculated abnormality degree series or the smoothed information in the abnormality degree storage table 1001 and passes it to the abnormality degree output unit 104.

図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 degree output unit 104 by graphing a series of abnormality degrees or information obtained by smoothing the series. The abnormality degree output screen (peak list) 1041 has a delay degree display area 1041A for displaying a delay degree graph showing the magnitude of the delay degree (abnormality degree) generated for each product, and a delay degree display area 1041A for each product in descending order of the degree of delay. It includes a peak list 1041B that displays the time when the peak value of the degree is taken, and a confirmation button 1041C that accepts an instruction to end the screen display. This screen allows the user to compare the peak delays of each product. Not limited to this, the abnormality degree output unit 104 may create and output the following output screen.

図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 degree output unit 104 by graphing a series of abnormality degrees or information obtained by smoothing the series. On the abnormality degree output screen (delay time estimation) 1042, a delay degree display area 1042A for displaying a delay degree graph showing the magnitude of the delay degree (abnormality degree) generated in the product, and a time point when a particularly large delay degree is detected. (Peak) includes a comment 1042B that displays the occurrence of a delay, and a confirmation button 1042C that accepts an instruction to end the screen display. From this screen, the user can know the time when a large delay occurs in the target product and get a hint for improving the work content.

図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 degree output unit 104 by graphing a series of abnormality degrees or information obtained by smoothing the series. On the abnormality degree output screen (delay degree) 1043, a delay degree display area 1043A for displaying a delay degree graph showing the magnitude of the delay degree (abnormality degree) generated in the product, a time zone in which the delay degree is detected, and the like thereof are displayed. A comment 1043B for displaying the amount of delay in time and a confirmation button 1043C for receiving an instruction to end the screen display are included. From this screen, the user can know the time when the delay occurs continuously for the target product and can obtain a hint for improving the work content.

なお、異常度出力画面(遅延度合)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 degree calculation unit 100 specifies the delay amount by calculating the difference (T1-T2) between the transition time T1 and the representative value T2. By doing so, the delay amount can be specified, and the abnormality degree output screen (delay degree degree) 1043 can be created. In addition, the user can also grasp the amount of delay by creating a graph for the amount of abnormality and using geometric information (specifically, an area indicating the amount of delay) on the graph.

[第二の実施形態]図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 abnormality detection unit 108 ′ according to the second embodiment is basically the same as the required time abnormality detection unit 108 according to the first embodiment, but there are some differences. Hereinafter, the difference from the required time abnormality detection unit according to the first embodiment will be mainly described.

第二の実施形態においては、モデル記憶部102に記憶される標準モデルを、時系列データ入力部103から得た単一あるいは複数の標準時系列データから学習可能としている。 In the second embodiment, the standard model stored in the model storage unit 102 can be learned from a single or a plurality of standard time series data obtained from the time series data input unit 103.

第二の実施形態にかかる所要時間異常検出部108´には、さらに、標準時系列モデル推定部106と、標準遷移時間モデル推定部107と、が含まれる。 The required time abnormality detection unit 108 ′ according to the second embodiment further includes a standard time series model estimation unit 106 and a standard transition time model estimation unit 107.

時系列データ入力部103から入力された標準時系列データは、標準時系列モデル推定部106と、状態推定部101と、タイムスタンプ取得部105と、に受け渡される。タイムスタンプ取得部105は、タイムスタンプ系列を算出もしくは取得し、標準遷移時間モデル推定部107に受け渡す。 The standard time series data input from the time series data input unit 103 is passed to the standard time series model estimation unit 106, the state estimation unit 101, and the time stamp acquisition unit 105. The time stamp acquisition unit 105 calculates or acquires a time stamp series and passes it to the standard transition time model estimation unit 107.

標準時系列モデル推定部106では、受け渡された単一または複数の標準時系列データを用いて、標準時系列モデルの数値パラメータを単一または複数推定してモデル記憶部102へ記憶させる。標準時系列モデルの学習は、隠れマルコフモデルや確率的オートマトンといった機械学習の枠組み等、AI(人工知能)等を用いて数値パラメータを推定するなどする。例えば隠れマルコフモデルを用いた場合、初期確率、遷移確率、出力確率が数値パラメータとして推定され、モデル記憶部102は、推定された数値パラメータ(初期確率表1021、出力確率表1022、遷移確率表1023)を標準時系列データごとに格納する。 The standard time-series model estimation unit 106 estimates a single or a plurality of numerical parameters of the standard time-series model using the passed single or a plurality of standard time-series data and stores them in the model storage unit 102. The learning of the standard time series model involves estimating numerical parameters using AI (artificial intelligence), etc., such as a machine learning framework such as a hidden Markov model or a probabilistic automaton. For example, when a hidden Markov model is used, the initial probability, the transition probability, and the output probability are estimated as numerical parameters, and the model storage unit 102 uses the estimated numerical parameters (initial probability table 1021, output probability table 1022, transition probability table 1023). ) Is stored for each standard time series data.

状態推定部101では、モデル記憶部102より読み出した標準時系列データの数値パラメータを用いて、時系列データ入力部103から入力された標準時系列データのセットに対して、標準時系列モデルを用いた状態推定を行い、時系列データの各データ点の状態変数を推定する。推定した状態変数について、状態推定部101は、標準遷移時間モデル推定部107に受け渡す。状態変数の推定には、標準時系列モデルと対応した状態推定手法を用いる。例えば隠れマルコフモデルであれば、ビタビアルゴリズムを用いて状態推定が可能である。また、複数の標準時系列データの数値パラメータを用いる場合、複数の状態が推定されるため、状態推定部101はその分布を求めるようにしてもよい。 The state estimation unit 101 uses the numerical parameters of the standard time series data read from the model storage unit 102 to estimate the state using the standard time series model for the set of standard time series data input from the time series data input unit 103. To estimate the state variable of each data point of the time series data. The state estimation unit 101 passes the estimated state variables to the standard transition time model estimation unit 107. For the estimation of state variables, the state estimation method corresponding to the standard time series model is used. For example, in the case of a hidden Markov model, the state can be estimated using the Viterbi algorithm. Further, when a plurality of numerical parameters of standard time series data are used, since a plurality of states are estimated, the state estimation unit 101 may obtain the distribution thereof.

標準遷移時間モデル推定部107は、状態推定部101から各データ点の状態変数を受け取り、タイムスタンプ取得部105から各点のタイムスタンプを受け取る。状態とタイムスタンプとを対応づけることで、標準遷移時間モデル推定部107は、対象の時系列データにおいて特定の状態が始まってから次の状態へ遷移するまでの遷移時間を算出する。 The standard transition time model estimation unit 107 receives a state variable of each data point from the state estimation unit 101, and receives a time stamp of each point from the time stamp acquisition unit 105. By associating the state with the time stamp, the standard transition time model estimation unit 107 calculates the transition time from the start of a specific state to the transition to the next state in the target time series data.

そして、標準遷移時間モデル推定部107は、各データ点の状態変数とタイムスタンプとの複数の標準遷移時間モデルにおけるすべてのペアに対して遷移時間の算出を行うことで、各状態についての遷移時間の分布を求め、モデル記憶部102の標準遷移時間モデル(遷移時間表1024)として記憶させる。なお、遷移時間の分布については、標準遷移時間モデル推定部107は、ガウス分布等で近似させてそのパラメータ(平均と分散など)を代表値として記憶する。これに限られず、すべての算出した遷移時間をそのまま保存しても良い。 Then, the standard transition time model estimation unit 107 calculates the transition time for all the pairs in the plurality of standard transition time models of the state variable and the time stamp of each data point, so that the transition time for each state is calculated. Is obtained and stored as a standard transition time model (transition time table 1024) of the model storage unit 102. The transition time distribution is approximated by the standard transition time model estimation unit 107 with a Gaussian distribution or the like, and the parameters (mean, variance, etc.) are stored as representative values. Not limited to this, all calculated transition times may be saved as they are.

このようにすることで、モデル記憶部102に記憶される単一または複数の標準モデルを、時系列データ入力部103から得た標準時系列データを用いて学習することができる。 By doing so, it is possible to learn a single or a plurality of standard models stored in the model storage unit 102 by using the standard time series data obtained from the time series data input unit 103.

[第三の実施形態]図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 abnormality detection unit 108 ″ according to the third embodiment is basically the same as the required time abnormality detection unit 108 according to the second embodiment, but there are some differences. Hereinafter, the difference from the required time abnormality detection unit according to the second embodiment will be mainly described.

第三の実施形態においては、所要時間異常検出部108´´では、単一の(一連の)標準時系列データを受け取り標準モデルとすることで、標準時系列モデル推定部106および標準遷移時間モデル推定部107を不要とする例である。 In the third embodiment, the required time abnormality detection unit 108 ″ receives a single (series) standard time series data and uses it as a standard model, so that the standard time series model estimation unit 106 and the standard transition time model estimation unit 106 ′. This is an example in which 107 is unnecessary.

第三の実施形態にかかる所要時間異常検出部108´´には、異常度計算部100と、状態推定部101と、モデル記憶部102と、タイムスタンプ取得部105と、が含まれる。 The required time abnormality detection unit 108 ″ according to the third embodiment includes an abnormality degree calculation unit 100, a state estimation unit 101, a model storage unit 102, and a time stamp acquisition unit 105.

時系列データ入力部103は、単一の標準時系列データの入力を受け付ける。時系列データ入力部103は、データ値を時系列に並べたデータ値系列はモデル記憶部102に、タイムスタンプ系列はタイムスタンプ取得部105に、それぞれ送る。 The time series data input unit 103 accepts input of a single standard time series data. The time-series data input unit 103 sends the data value series in which the data values are arranged in time series to the model storage unit 102, and the time stamp series to the time stamp acquisition unit 105.

タイムスタンプ取得部105は、入力された標準時系列データのタイムスタンプ系列を取得し、モデル記憶部102に受け渡す。 The time stamp acquisition unit 105 acquires the time stamp series of the input standard time series data and passes it to the model storage unit 102.

そして、モデル記憶部102は、時系列データ入力部103から受け渡された標準時系列データのデータ値系列と、タイムスタンプ取得部105から受け渡されたタイムスタンプ系列と、を記憶する。これが、単一の標準モデルの学習の流れとなる。次に、異常度計算時の流れを説明する。 Then, the model storage unit 102 stores the data value series of the standard time series data passed from the time series data input unit 103 and the time stamp series passed from the time stamp acquisition unit 105. This is the learning flow for a single standard model. Next, the flow at the time of calculating the degree of abnormality will be described.

時系列データ入力部103において異常検出の対象となる対象時系列データを受け付けると、時系列データ入力部103は、データ値系列は状態推定部101に、タイムスタンプ系列もしくはその算出に必要な情報はタイムスタンプ取得部105にそれぞれ送る。タイムスタンプ取得部105では、対象時系列データのタイムスタンプ系列を取得もしくは算出し、異常度計算部100へ受け渡す。 When the time-series data input unit 103 receives the target time-series data to be detected as an abnormality, the time-series data input unit 103 sends the data value series to the state estimation unit 101, and the time stamp series or the information necessary for its calculation. It is sent to the time stamp acquisition unit 105 respectively. The time stamp acquisition unit 105 acquires or calculates the time stamp series of the target time series data and passes it to the abnormality degree calculation unit 100.

状態推定部101は、受け渡された対象時系列データのデータ値系列と、モデル記憶部102から読み出した標準時系列データのデータ値系列について整合演算を行い、整合演算の結果を異常度計算部100へ受け渡す。 The state estimation unit 101 performs a matching operation on the data value series of the passed target time series data and the data value series of the standard time series data read from the model storage unit 102, and the result of the matching calculation is obtained by the abnormality degree calculation unit 100. Hand over to.

図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 state estimation unit 101 adopts and specifies a known algorithm such as DTW.

なお、異常度計算部100による異常度の算出に関しては、基本的に第一の実施形態および第二の実施形態と同様である。第三の実施形態の構成においては、標準時系列モデルの数値パラメータとして単一の標準時系列データを用いるため、第二の実施形態に示した複数の標準時系列データを用いてモデリングを行う標準時系列モデル推定部106が不要となる。また、標準時系列データの各データ点を異なる状態として扱い、整合演算によって対象時系列の各データ点を対応付けることで状態推定を実現しているため、各状態の標準遷移時間は単一の値であり、タイムスタンプ系列の隣り合ったタイムスタンプの差分と同一となる。そのため、標準遷移時間モデルとしてはタイムスタンプ系列をそのまま記憶すればよく、標準遷移時間モデル推定部107についても不要である。したがって、第三の実施形態に係る所要時間異常検出部108´´は、その構成を簡便なものとすることができる。 The calculation of the degree of abnormality by the abnormality degree calculation unit 100 is basically the same as that of the first embodiment and the second embodiment. In the configuration of the third embodiment, since a single standard time series data is used as the numerical parameters of the standard time series model, the standard time series model estimation is performed by modeling using the plurality of standard time series data shown in the second embodiment. The unit 106 becomes unnecessary. In addition, since each data point of the standard time series data is treated as a different state and the state estimation is realized by associating each data point of the target time series by the matching operation, the standard transition time of each state is a single value. Yes, it is the same as the difference between adjacent time stamps in the time stamp series. Therefore, the time stamp series may be stored as it is as the standard transition time model, and the standard transition time model estimation unit 107 is also unnecessary. Therefore, the required time abnormality detection unit 108 ″ according to the third embodiment can be simplified in its configuration.

[第四の実施形態]図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 abnormality detection system 205 that detects an abnormality in manual work by a wearable sensor (a detachable sensor that acquires the amount of change in position). be.

作業者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 sensor receiving unit 204, the information processing terminal 203 inputs the sensor data to the work abnormality detection system 205.

図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 sensor receiving unit 204 to the feature quantity extraction unit 206, and the feature quantity extraction unit 206 converts the sensor data into time-series data composed of feature quantities suitable for abnormality detection. For example, if the sensor data is missing, the feature amount extraction unit 206 performs complementary processing and removes noise contained in the sensor data with a filter such as a low-pass filter, so that the accuracy of abnormality detection can be expected to be improved. .. Further, the feature amount extraction unit 206 obtains the feature amount by Fourier transforming the sensor data.

特徴量抽出部206により特徴量が抽出された時系列データは、系列区間切り出し部207へ受け渡され、切り出し制御部210からの制御信号を元に、繰り返し作業における一回の作業に対応する区間が切り出されて、所要時間異常検出部108へ受け渡される。 The time-series data from which the feature amount has been extracted by the feature amount extraction unit 206 is passed to the series section cutout unit 207, and based on the control signal from the cutout control unit 210, the section corresponding to one operation in the repetitive work. Is cut out and handed over to the required time abnormality detection unit 108.

所要時間異常検出部108は、上述の第一、第二、第三の実施形態に示したように、受け取った時系列データに対して異常度を算出し、異常度出力部104へ算出結果を送る。異常度出力部104は、異常度を記録し、ディスプレイ等に表示する。このようにすることで、例えば、画面に表示された異常度を確認することで、正常に作業が行われているかを判断することが可能となる。特定の基準よりも異常度が大きい場合は警告を出力し、作業者や監督者などが対処を容易となるように支援するようにしても良い。 The required time abnormality detection unit 108 calculates the degree of abnormality with respect to the received time-series data as shown in the first, second, and third embodiments described above, and outputs the calculation result to the abnormality degree output unit 104. send. The abnormality degree output unit 104 records the abnormality degree and displays it on a display or the like. By doing so, for example, by checking the degree of abnormality displayed on the screen, it becomes possible to determine whether the work is being performed normally. If the degree of abnormality is greater than a specific standard, a warning may be output to assist workers and supervisors in making it easier to deal with.

[第五の実施形態]図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 abnormality detection system 205 that detects abnormalities in manual work by a camera.

作業者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 abnormality detection system 205.

情報処理端末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 abnormality detection system 205 according to the present embodiment is basically the same as the work abnormality detection system 205 according to the fourth embodiment, but there are some differences. Hereinafter, the differences from the work abnormality detection system 205 according to the fourth embodiment will be mainly described.

特徴量抽出部206は、動画データ内の被写体の関節位置を特徴量として抽出する。具体的には、特徴量抽出部206は、動画データを解析して骨格認識を行う。例えば、骨格認識アルゴリズムを用いて、特徴量抽出部206は、骨格を認識する。 The feature amount extraction unit 206 extracts the joint position of the subject in the moving image data as the feature amount. Specifically, the feature amount extraction unit 206 analyzes the moving image data and performs skeleton recognition. For example, using a skeleton recognition algorithm, the feature amount extraction unit 206 recognizes a skeleton.

図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 section cutting unit 207 is the same as that of the fourth embodiment.

第五の実施形態によれば、例えば、画面に表示された異常度を確認することで、正常な速度で作業が行われているかを判断することが可能となる。特定の基準よりも異常度が大きい場合は警告を出力し、作業者や監督者などがカメラ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 abnormality detection system 205 that detects the abnormality of the machine tool work by the current sensor.

工作機械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 current sensor 502 sends the measured current value to the information processing terminal 503 via wired communication or wireless communication.

情報処理端末503は、電流センサ502から電流値の時系列データを受け取ると、作業異常検出システム205へセンサデータを入力する。 When the information processing terminal 503 receives the time-series data of the current value from the current sensor 502, the information processing terminal 503 inputs the sensor data to the work abnormality detection system 205.

本実施形態に係る作業異常検出システム205は、基本的に第四の実施形態に係る作業異常検出システム205と同様であるが、一部相違がある。以下、第四の実施形態にかかる作業異常検出システム205との相違を中心に説明する。 The work abnormality detection system 205 according to the present embodiment is basically the same as the work abnormality detection system 205 according to the fourth embodiment, but there are some differences. Hereinafter, the differences from the work abnormality detection system 205 according to the fourth embodiment will be mainly described.

特徴量抽出部206は、電流値の時系列データを対象として、センサデータを異常検出に適した特徴量からなる時系列データに変換する。例えば、センサデータに含まれるノイズをローパスフィルタ等のフィルタで除去することで、異常検出の精度向上が期待できる。例えば、特徴量抽出部206は、センサデータに欠落がある場合には、補完処理を行い、センサデータに含まれるノイズをローパスフィルタ等のフィルタで除去することで、異常検出の精度向上が期待できる。また、特徴量抽出部206は、センサデータをフーリエ変換することで特徴量を得る。 The feature amount extraction unit 206 converts the sensor data into time-series data composed of feature amounts suitable for abnormality detection for the time-series data of the current value. For example, by removing the noise contained in the sensor data with a filter such as a low-pass filter, the accuracy of abnormality detection can be expected to be improved. For example, if the sensor data is missing, the feature amount extraction unit 206 performs complementary processing and removes noise contained in the sensor data with a filter such as a low-pass filter, so that the accuracy of abnormality detection can be expected to be improved. .. Further, the feature amount extraction unit 206 obtains the feature amount by Fourier transforming the sensor data.

特徴量抽出された時系列データが系列区間切り出し部207へ受け渡された後の処理は、第四の実施形態と同様である。 The processing after the time-series data extracted by the feature amount is passed to the series section cutting unit 207 is the same as that of the fourth embodiment.

第六の実施形態によれば、例えば、画面に表示された異常度を確認することで、正常な速度で工作機械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 degree output unit 104 holds an identifier of an event that occurs on the standard time series data, and time information about the time and duration of the event. The standard identification information 402 is provided, and the standard identification information 402 is connected to the abnormality degree display unit 401. The standard identification information 402 includes, for example, in the example of work in a factory as shown in the fourth to sixth embodiments, the order of work in a series of work and the standard work time of each are predetermined. In many cases, information is stored in advance based on this.

所要時間異常検出部108から特定時間の異常度についての情報が渡されると、異常度表示部401は、標準識別情報402を参照し、標準時系列データ上において該当時間に発生する事象の識別子を検索し、異常度と識別子を対応付けて表示画面を作成する。 When the information about the abnormality degree of the specific time is passed from the required time abnormality detection unit 108, the abnormality degree display unit 401 refers to the standard identification information 402 and searches for the identifier of the event occurring at the corresponding time on the standard time series data. Then, create a display screen by associating the degree of abnormality with the identifier.

図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 degree output unit 104 by graphing a series of abnormality degrees or information obtained by smoothing the series. On the abnormality degree output screen (work control) 1044, a delay degree display area 1044A for displaying a delay degree graph showing the magnitude of the delay degree (abnormality degree) generated in the product and an identifier display for displaying a work identifier for each time point are displayed. An area 1044B and a confirmation button 1044C for receiving an instruction to end the screen display are included. From this screen, the user can know the work at the time when a large delay occurs for the target product and get a hint for improving the work content.

ここで、標準識別情報402が保持する識別子に対する時間情報は、あくまで標準時のものであり、対象時系列データにおいて時間異常が頻発した場合は、その誤差が無視できないほど大きい可能性がある。そこで、異常度表示部401は、異常度を元に対象時系列データにおける標準の場合との時間のずれを算出し、時間情報を補正した上で識別子を用いるようにしても良い。具体的には、異常度計算部100において、開始時刻から特定経過時刻までの区間を局所区間に区切り、各局所区間ABの異常度をT1およびT2の差分で求めた場合(T1、T2の定義は第一の実施形態と同様)、この累積値を算出することで、特定時間までの標準時に対する時間のずれを算出することができ、これを用いて時間情報を補正した上で表示に用いることができる。以上が、第七の実施形態である。 Here, the time information for the identifier held by the standard identification information 402 is only for standard time, and if a time abnormality occurs frequently in the target time series data, the error may be so large that it cannot be ignored. Therefore, the abnormality degree display unit 401 may calculate the time lag from the standard case in the target time series data based on the abnormality degree, correct the time information, and then use the identifier. Specifically, in the abnormality degree calculation unit 100, the section from the start time to the specific elapsed time is divided into local sections, and the abnormality degree of each local section AB is obtained by the difference between T1 and T2 (definition of T1 and T2). Is the same as in the first embodiment), by calculating this cumulative value, it is possible to calculate the time lag with respect to the standard time up to a specific time, and use this to correct the time information before using it for display. Can be done. The above is the seventh embodiment.

なお、本発明は上記した実施形態に限定されるものではなく、様々な変形例が含まれる。例えば、上記した実施例は本発明を分かりやすく説明するために詳細に説明したものであり、必ずしも説明した全ての構成を備えるものに限定されるものではない。また、ある実施例の構成の一部を他の実施例の構成に置き換えることが可能であり、また、ある実施例の構成に他の実施例の構成を加えることも可能である。 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.
請求項1に記載の異常検出装置であって、
前記異常検出の対象となる時系列データからタイムスタンプを取得するタイムスタンプ取得部を備え、
前記異常度計算部は、前記標準遷移時間モデル数値パラメータと、前記状態推定部が推定した前記各データ点の状態変数と、前記タイムスタンプと、を用いて各状態から別の状態への遷移に要する所要時間の異常度を算出する、
ことを特徴とする異常検出装置。
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.
請求項1に記載の異常検出装置であって、
前記標準時系列モデル数値パラメータと、前記標準遷移時間モデル数値パラメータとは、前記モデル記憶部にそれぞれ複数記憶されており、
前記状態推定部は、複数の前記標準時系列モデル数値パラメータによる状態推定を行った結果を統合して前記状態変数を推定する、
ことを特徴とする異常検出装置。
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.
請求項1に記載の異常検出装置であって、
受け渡された単一または複数の標準時系列データを用いて、標準時系列モデルの数値パラメータを単一または複数推定して前記標準時系列モデル数値パラメータとして前記モデル記憶部へ記憶させる標準時系列モデル推定部、
を備えることを特徴とする異常検出装置。
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.
請求項1に記載の異常検出装置であって、
受け渡された単一または複数の標準時系列データにおいて特定の状態が始まってから次の状態へ遷移するまでの遷移時間を算出して前記標準遷移時間モデル数値パラメータとして前記モデル記憶部へ記憶させる標準遷移時間モデル推定部、
を備えることを特徴とする異常検出装置。
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.
請求項1に記載の異常検出装置であって、
前記時系列データ入力部は、位置の変化量を取得するセンサから出力される時系列データの入力を受け付ける、
ことを特徴とする異常検出装置。
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.
請求項1に記載の異常検出装置であって、
動画データを解析して骨格認識を行う特徴量抽出部を備え、
前記時系列データ入力部は、カメラから出力される動画データの入力を受け付け、
前記状態推定部は、前記特徴量抽出部により前記動画データを解析した結果得られる骨格の関節位置の変化を時系列データとして状態を推定する、
ことを特徴とする異常検出装置。
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.
請求項1に記載の異常検出装置であって、
前記時系列データ入力部は、電流値を取得するセンサから出力される時系列データの入力を受け付ける、
ことを特徴とする異常検出装置。
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.
請求項1に記載の異常検出装置であって、
前記異常度出力部は、前記基準となる時系列データ上で発生する事象の識別子と、該事象の発生時点および継続期間についての時間情報と、を含む標準識別情報を予め記憶しており、
前記グラフ化において前記所定の状態に対応した前記事象の識別子を対応付ける、
ことを特徴とする異常検出装置。
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.
請求項1に記載の異常検出装置であって、
前記異常度出力部は、前記グラフに対応付けて、前記遅延度のピーク値が大きくなる経過時間を示す、
ことを特徴とする異常検出装置。
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.
請求項1に記載の異常検出装置であって、
前記異常度出力部は、前記グラフの遅延度のピーク値が大きくなる経過時間の位置に、遅延発生のコメントを示す、
ことを特徴とする異常検出装置。
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.
請求項1に記載の異常検出装置であって、
前記異常度出力部は、前記グラフを前記異常度の量について作成し、前記グラフ上の幾何的情報により遅延度の量を示す、
ことを特徴とする異常検出装置。
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.
JP2019044798A 2019-03-12 2019-03-12 Anomaly detection device and anomaly detection method Active JP7072531B2 (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
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