(57)【要約】【課題】構造物の取付状態の劣化を、検査員よる定期検査を必要とせずに、定量的に診断する。【解決手段】センサヘッドから取得した加速度情報に基づく特徴量として、傾き情報および固有振動数情報の少なくとも一方の情報を抽出する特徴量抽出手段（Ｓ２０１）と、特徴量抽出手段により抽出された特徴量に関して、正常状態における基準用確率密度分布と、劣化診断時における診断用確率密度分布を算出する確率密度分布算出手段（Ｓ２０４）と、確率密度分布算出手段で算出された基準用確率密度分布と診断用確率密度分布との間に所定量以上の有意差を有する分布が存在する場合には、構造物の取付状態に異常が発生していると判定する判定手段（Ｓ２０５、Ｓ２０６）とを備える。【選択図】図２
2014228472.tif 000002
【特許請求の範囲】
【請求項１】
固定面に取り付けられ、劣化診断対象である構造物に設置されたセンサヘッドから出力される前記構造物の加速度情報に基づいて、前記構造物の取付状態が正常か否かを診断するセンサコントローラを備えた構造物劣化診断システムにおいて、前記センサコントローラにより実行される構造物劣化診断システムであって、
前記センサヘッドから取得した前記加速度情報に基づく特徴量として、前記構造物の傾き情報および固有振動数情報の少なくとも一方の情報を抽出する特徴量抽出手段と、
前記特徴量抽出手段により抽出された前記特徴量に関して、正常状態における基準用確率密度分布と、劣化診断時における診断用確率密度分布を算出する確率密度分布算出手段と、
前記確率密度分布算出手段で算出された前記基準用確率密度分布と前記診断用確率密度分布との間に所定量以上の有意差を有する分布が存在する場合には、前記構造物の取付状態に異常が発生していると判定する判定手段と
を備えた構造物劣化診断システム。
【請求項２】
請求項１に記載の構造物劣化診断システムにおいて、
前記センサヘッドからは、３軸の加速度情報が出力され、
前記特徴量抽出手段は、３軸分の傾き情報および３軸分の固有振動数情報の少なくとも一方の情報を抽出する
構造物劣化診断システム。
【請求項３】
請求項２に記載の構造物劣化診断システムにおいて、
前記センサコントローラにおいて、前記傾き情報を前記特徴量として抽出することで劣化診断を行う場合には、
前記特徴量抽出手段は、前記センサヘッドから取得した前記３軸の加速度情報のそれぞれについて低周波成分を抽出することで、３軸に対する前記傾き情報を取得し、
前記確率密度分布算出手段は、前記特徴量抽出手段で取得されたそれぞれの軸の前記傾き情報をデシベル化し、所定のデシベル単位に基づいて、学習時および劣化診断時の確率密度分布を算出し、
前記判定手段は、３軸のそれぞれについて、劣化診断時の確率密度分布を学習時の確率密度分布で割ることで確率密度比を算出し、算出した３軸の前記確率密度比のユークリッド距離としてベクトル長を求め、劣化診断のための評価値とし、前記評価値が許容傾き基準値を超える場合には、前記構造物の傾きが許容できないレベルに達していると判断する
構造物劣化診断システム。
【請求項４】
請求項２に記載の構造物劣化診断システムにおいて、
前記センサコントローラにおいて、前記固有振動数情報を前記特徴量として抽出することで劣化診断を行う場合には、
前記特徴量抽出手段は、前記センサヘッドから取得した前記３軸の加速度情報に関して、いずれかの軸で基準レベルを越えた加速度データが得られた時点を含む一定数の加速度データを抽出することで、固有振動数情報を取得し、
前記確率密度分布算出手段は、前記特徴量抽出手段で取得された前記固有振動数情報について周波数解析を行ってパワースペクトルを算出し、前記パワースペクトルが最も大きい周波数に相当する第１固有振動数と、前記パワースペクトルが２番目に大きい周波数に相当する第２固有振動数を算出し、前記第１固有振動数および前記第２固有振動数による２次元の直交座標として各軸の学習時および劣化診断時の２次元の確率密度分布を算出し、
前記判定手段は、３軸のそれぞれについて、劣化診断時の２次元の確率密度分布を学習時の２次元の確率密度分布で割ることで２次元の確率密度比を算出し、劣化診断のための評価値とし、３軸のうちの少なくともいずれか１つの軸において、前記評価値が許容振動基準値を超える場合には、前記構造物の振動が許容できないレベルに達していると判断する
構造物劣化診断システム。
【請求項５】
請求項１から４のいずれか１項に記載の構造物劣化診断システムにおいて、
前記特徴量抽出手段により抽出された前記特徴量に対して、前記センサヘッドの取り付け方向を修正するキャリブレーション処理を行う誤差修正手段
をさらに備え、
前記確率密度分布算出手段は、前記誤差修正手段により前記キャリブレーション処理が施された後の特徴量に関して、前記基準用確率密度分布および診断用確率密度分布を算出する
構造物劣化診断システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、例えば、トンネル内の吊り下げ構造物、あるいは張り出し構造物のような構造物の取付状態の劣化を、長期にわたってモニタすることのできる構造物劣化診断システムに関する。
【背景技術】
【０００２】
構造物の劣化状態を検出する方法としては、検査員による定期検査により、目視あるいは何らかの計器を用いて行われることが主流であった。また、劣化診断対象であるトンネルなどの構造物に経年的に発生する亀裂に関して、定量的な検査を、簡単かつ迅速に行う従来技術が開示されている（例えば、特許文献１参照）。
【０００３】
この特許文献１では、紫外線または青色系可視光などの励起光によって発光する蛍光色素を、劣化診断対象である構造物にあらかじめ混入させている。そして、この構造物に紫外線または青色系可視光などを発光する光源を照射し、目視あるいはＣＣＤカメラ等による撮像画像の解析処理により、亀裂の発生を定量的に判断している。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１３−８３４９３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、従来技術には、以下のような課題がある。
特許文献１では、定量的な劣化診断を可能にしてはいるものの、あくまでも、検査員による定期検査を基本としている。さらに、劣化診断対象の構造物に対して、蛍光色素をあらかじめ混入させておく必要があった。
【０００６】
一方、近年では、構造物の劣化診断を定期検査よりも短い周期で、検査員を介さずに無人で行うことのできる劣化診断システムが望まれている。また、トンネル等においては、トンネル自身の壁面や天井の劣化状態以外にも、トンネル内の吊り下げ構造物、あるいは張り出し構造物のような構造物の取付状態の劣化を、定量的に診断する必要性がある。さらに、新規の構造物だけでなく、既存の構造物に対しても、容易に対応できることが望まれる。
【０００７】
本発明は、前記のような課題を解決するためになされたものであり、構造物の取付状態の劣化を、検査員よる定期検査を必要とせずに、定量的に診断することのできる構造物劣化診断システムを得ることを目的とする。
【課題を解決するための手段】
【０００８】
本発明に係る構造物劣化診断システムは、固定面に取り付けられ、劣化診断対象である構造物に設置されたセンサヘッドから出力される構造物の加速度情報に基づいて、構造物の取付状態が正常か否かを診断するセンサコントローラを備えた構造物劣化診断システムにおいて、センサコントローラにより実行される構造物劣化診断システムであって、センサヘッドから取得した加速度情報に基づく特徴量として、構造物の傾き情報および固有振動数情報の少なくとも一方の情報を抽出する特徴量抽出手段と、特徴量抽出手段により抽出された特徴量に関して、正常状態における基準用確率密度分布と、劣化診断時における診断用確率密度分布を算出する確率密度分布算出手段と、確率密度分布算出手段で算出された基準用確率密度分布と診断用確率密度分布との間に所定量以上の有意差を有する分布が存在する場合には、構造物の取付状態に異常が発生していると判定する判定手段とを備えるものである。
【発明の効果】
【０００９】
本発明によれば、加速度情報に基づく特徴量から算出された確率密度分布の変化度合を検出することで、定量的な劣化診断を可能とすることにより、構造物の取付状態の劣化を、検査員よる定期検査を必要とせずに、定量的に診断することのできる構造物劣化診断システムを得ることができる。
【図面の簡単な説明】
【００１０】
【図１】本発明の実施の形態１における構造物劣化診断システムの構成図である。
【図２】本発明の実施の形態１における第１の劣化診断部による劣化診断処理の流れを示すフローチャートである。
【図３】本発明の実施の形態１における３軸分の傾き情報の収集に関する具体例を示した図である。
【図４】本発明の実施の形態１における３軸分の傾き情報に関する確率密度分布の具体例を示した図である。
【図５】本発明の実施の形態１における、ある１軸の傾き情報に関する確率密度分布、確率密度比、および評価値の具体例を示した図である。
【図６】本発明の実施の形態１における第２の劣化診断部による劣化診断処理の流れを示すフローチャートである。
【図７】本発明の実施の形態１における３軸分の振動数情報の収集に関する具体例を示した図である。
【図８】本発明の実施の形態１における３軸分のパワースペクトルの算出結果を示した図である。
【図９】本発明の実施の形態１における、ある１軸について算出された学習時の確率密度分布、劣化診断時の確率密度分布、および両者の比である確率密度比の具体例を示した図である。
【発明を実施するための形態】
【００１１】
以下、本発明の構造物劣化診断システムの好適な実施の形態につき、図面を用いて説明する。
本発明は、劣化診断対象である構造物に取り付けられたセンサから得られる加速度情報に基づいて、構造物の傾きあるいは固有振動数に関する特徴量の確率密度分布を求め、正常時における確率密度分布と劣化診断時における確率密度分布との間に有意差がある場合には、取付状態に何らかの劣化が発生していると判断することを技術的特徴としている。
【００１２】
実施の形態１．
図１は、本発明の実施の形態１における構造物劣化診断システムの構成図である。本実施の形態１における構造物劣化診断システムは、センサコントローラ１０、およびＮ個（Ｎは、２以上の整数）のセンサヘッド２０（１）〜２０（Ｎ）を備えて構成されている。
【００１３】
なお、本実施の形態１における構造物劣化診断システムでは、最低限、１個のセンサヘッド２０を設けておけば、劣化診断を実施することが可能である。また、複数用いる場合のＮ個のセンサヘッドのそれぞれの機能は、全て共通である。そこで、以下の説明では、それぞれのセンサヘッドを区別する必要がない場合には、（１）〜（Ｎ）の添字を用いずに、単にセンサヘッド２０と記載する。
【００１４】
Ｎ個のセンサヘッド２０のそれぞれは、センサ部２１と、加速度情報出力部２２を有しており、劣化診断対象である構造物の異なる位置に設置されている。ここで、構造物の具体例としては、トンネル内の吊り下げ構造物、あるいは張り出し構造物が挙げられ、本発明の構造物劣化診断システムによって、長期にわたって構造物の取付状態の劣化診断が行われることとなる。
【００１５】
このような構造物としては、トンネル内に設置されるジェットファン、情報板、標識等が対象とされ、トンネルを形成する躯体のコンクリートから吊り下げまたは張り出される重量物が挙げられる。そして、構造物に対してセンサヘッド２０が好ましくは２個設置され、各々を構造物に対して両端となる部分に設置することで、構造物の変位をいずれかが大きく捕らえることとなる。
【００１６】
そのため、Ｎ個のセンサヘッド２０は、構造物の変位を大きく捕らえるために、設置面の方向に対して、各々が離れる位置に設置される。さらに、センサヘッド２０は、交換を容易にするため、２個等の複数個を２組としてもよく、１組ごとで交換することにより、監視状態を継続することができる。
【００１７】
センサ部２１は、例えば、薄膜の水晶振動子を用い、応答性に優れ、測定範囲がＤＣ〜数十Ｈｚ程度の加速度を測定可能な３軸加速度センサである。このように、加速度センサに３軸のセンサを用いることにより、水平出しが不要となり、傾きや振動の方向に関わらず、センサ出力を行うことができる。したがって、水平出しを行うなど、傾きの方向などが特定できる場合には、２軸や１軸であってもよい。
【００１８】
また、加速度情報出力部２２は、センサヘッドの設置箇所における構造物の３軸の加速度に関するアナログ信号を、所定のサンプリングレート（例えば、５０Ｈｚのサンプリングレート）でデジタル信号に変換し、加速度情報としてセンサコントローラ１０へ送信する。なお、加速度情報出力部２２は、一例として、ＣＡＮやＲＳ２３２Ｃ等の通信Ｉ／Ｆを介して、デジタル信号である加速度情報をセンサコントローラ１０に送信することができる。
【００１９】
また、センサコントローラ１０は、構造物近傍のＮ個のセンサヘッド２０に配線が届く位置に設置され、Ｎ個のセンサヘッド２０のそれぞれから取得した加速度情報に基づいて、構造物の取付状態が正常か否かを判断する。具体的には、センサコントローラ１０は、加速度情報に基づいて、傾き情報を第１の特徴量として抽出し、構造物の取付状態の劣化診断を行う第１の劣化診断部１１と、加速度情報に基づいて、固有振動数情報を第２の特徴量として抽出し、構造物の取付状態の劣化診断を行う第２の劣化診断部１２を有している。
【００２０】
また、センサコントローラ１０は、診断結果、あるいは診断に用いたデータ等を、図示しない制御装置に送信できるようにするために、例えば、イーサネット（登録商標）等の大容量通信Ｉ／Ｆを備えることができる。
【００２１】
構造物としてジェットファン等が設けられるトンネル設備の場合、詳細に示さないが、トンネル内の情報は、トンネル近傍に設けられた、いわゆる電気室内に設置される制御装置に集約され、さらに、監視者の居る遠方監視制御装置に移報されて、双方で情報が共有される。後述する劣化診断の処理に従い、劣化が生じていると判定された結果が、上記各制御装置に表示や警報音等で報知される。このような判定結果は、別途劣化診断システムの設備業者や点検業者のセンタ装置に送られてもよく、異常発生時に迅速な対応を可能にできる。
【００２２】
また、劣化を判定するレベルを次のように細分し、そのレベルに応じて対応を変えてもよい。例えば、劣化の判定レベルは、
レベル１：「定期点検等で詳細に点検する必要有り（確認レベル）」
レベル２：「現場を確認し、今後の改良計画等を検討する必要あり（計画レベル）」
レベル３：「速やかに改善を必要とする事態で、至急現場へ急行する必要有り（改善レベル）」
レベル４：「崩落を含む緊急事態（緊急レベル）」
とレベル分けし、レベル１、２は、各制御装置の盤面で簡単に表示するだけとし、レベル３、４の段階になって、詳細な警報表示や警報音鳴動等を行うようにしてもよい。
【００２３】
さらに、業者のセンタ装置で常時判定に用いるデータを収集し、別途、多面的に分析を行い、例えば、レベル３に達するまでの時期予測などを行ってもよい。トンネルは、遠方監視制御装置によって、常時監視員に監視されていることが多いが、そうでない場合もあり、日常の管理設備が十分でない場合は、業者のセンタ装置を介して、緊急時の一次対応として現場を確認し、状況を報告する対応を、業者の係員に行わせることも可能である。
【００２４】
そして、第１の劣化診断部１１および第２の劣化診断部１２は、いずれも、監視対象である構造物の正常時における確率密度分布と、劣化診断時における確率密度分布とを比較し、有意差が生じたことで構造物の劣化を検出する点では共通している。ただし、傾き情報を特徴量として劣化診断を行う第１の劣化診断部１１と、固有振動数情報を特徴量として劣化診断を行う第２の劣化診断部１２とでは、具体的な信号処理方法が異なっており、以下において、詳細に説明する。
【００２５】
（１）第１の劣化診断部１１による傾き情報に基づく劣化診断の処理の流れについて
図２は、本発明の実施の形態１における第１の劣化診断部１１による劣化診断処理の流れを示すフローチャートである。この図２のフローチャートに基づいて、傾き情報に基づいた第１の劣化診断部１１による具体的な劣化診断処理について説明する。
【００２６】
まず始めに、ステップＳ２０１において、第１の劣化診断部１１は、センサヘッド２０内の加速度情報出力部２２から取得した加速度情報に基づいて、傾き情報の抽出を行う。具体的には、第１の劣化診断部１１は、例えば、５０ＨｚのサンプリングレートでＸ、Ｙ、Ｚ軸の加速度情報を取得する。そして、第１の劣化診断部１１は、この加速度情報を用いて、各軸毎にディジタルフィルタによるローパスフィルタ処理を行う。ローパスフィルタのカットオフ周波数は、例えば、０．１Ｈｚ程度とする。
【００２７】
ここで、ＤＣ成分を含む３軸加速度の低周波成分は、傾き情報と等価である。そして、ローパスフィルタ処理された傾き情報は、急速に変化する高周波成分を含まない。このため、劣化検出のための傾きデータのサンプリングレートを遅くしても、情報は欠損しない。具体的には、カットオフ周波数の２倍程度のサンプリングレートがあれば、標本化定理を満たすことができ、０．１Ｈｚのカットオフ周波数に対しては０．２Ｈｚのサンプリング周波数、すなわち、５秒間隔のサンプリングでよい。
【００２８】
図３は、本発明の実施の形態１における３軸分の傾き情報の収集に関する具体例を示した図である。図３に示すように、第１の劣化診断部１１は、サンプリング間隔ΔＴごとに、３軸分の傾き情報を収集することとなる。
【００２９】
次に、第１の劣化診断部１１は、ステップＳ２０２において、センサヘッド２０の取り付け方向の修正を行う。具体的には、第１の劣化診断部１１は、センサヘッド２０の取り付け方向のキャリブレーションデータ（アフィン変換の行列の要素）と、先のステップＳ２０１で得られたサンプリングデータを行列演算することで、それぞれのセンサヘッド２０の座標軸を一致させるように、取り付け方向の誤差修正を行う。
【００３０】
このステップＳ２０２は、センサヘッド２０の取付方向に依存せず座標軸を揃えることが目的であり、センサヘッド２０が３軸で方向に関わらず傾きや振動を検出できることから、センサ出力をそのまま用いる場合には、なくてもよい。
【００３１】
次に、第１の劣化診断部１１は、ステップＳ２０３において、センサヘッド２０の取り付け方向の誤差修正後の傾き情報に対して、対数化処理（デシベル化）を行う。微小な傾き情報の計測のために、例えば、１／１００万［Ｇ］単位の分解能を元にして、後述する確率密度分布を得ようとすれば、そのために必要となる記憶領域は、膨大となる。そこで、本実施の形態１の第１の劣化診断部１１は、加速度を対数化（デシベル化）し、デシベルデータに基づいた分類により確率密度分布を求めることで、効率のよい量子化を行っている。
【００３２】
次に、第１の劣化診断部１１は、ステップＳ２０４において、傾き情報に基づく劣化診断を行うために、確率密度分布の生成処理を行う。第１の劣化診断部１１は、例えば、先のステップＳ２０３でデシベル化された傾き情報を用いて、０．５デシベル単位で累積度数分布を作成し、これに局所平滑化処理を施すことで、確率密度分布を生成する。なお、平滑化処理としては、例えば、ガウス関数を適用することができる。
【００３３】
なお、設置時の初期学習データとしての確率密度分布は、例えば、１ケ月程度のデータを用いて生成される。一方、劣化診断時の現在のデータとして作成する確率密度分布のデータ数は、例えば、過去３時間〜２４時間程度のデータに基づいて生成される。
【００３４】
また、基準データとなる正常時における学習データとしては、必ずしも設置時に取得したデータを継続して使用する必要はない。長期にわたって劣化診断を行うため、劣化診断を行う過程で、継続的に学習データを更新することも可能である。そして、劣化監視時に継続学習データとして採用する確率密度分布は、例えば、過去数時間〜数週間程度のデータに基づいて生成される。
【００３５】
図４は、本発明の実施の形態１における３軸分の傾き情報に関する確率密度分布の具体例を示した図である。各軸の１次元の確率密度分布は、傾きなしの中点に対して、正方向の傾きおよび負方向の傾きが大きくなるに従って、中点から両側に広がるような分布として表される。
【００３６】
次に、第１の劣化診断部１１は、ステップＳ２０５において、劣化診断時に求めた各軸の確率密度分布（診断時データに相当）を、学習時に求めた各軸の確率密度分布（基準データに相当）で割ることで、両者の比率を計算し、それぞれの軸について確率密度比の分布データを生成する。なお、基準データとしては、初期設定時のデータ以外に、継続学習データにより更新された学習データを採用することもできる。
【００３７】
図５は、本発明の実施の形態１における、ある１軸の傾き情報に関する確率密度分布、確率密度比、および評価値の具体例を示した図である。図５（ａ）、（ｂ）に示すように、確率密度分布のピークがずれた際には、そのずれたピーク位置を中心に、１よりも大きな確率密度比の分布が発生することとなる。
【００３８】
さらに、第１の劣化診断部１１は、各軸について算出された傾きに関する確率密度比をベクトルに見立てて、ベクトル長を計算する。具体的には、第１の劣化診断部１１は、図５（ｃ）に示すように、３軸の確率密度比のユークリッド距離としてベクトル長を求め、劣化診断のための評価値とする。
【００３９】
そして、最後に、ステップＳ２０６において、第１の劣化診断部１１は、確率密度比のベクトル長が、劣化判定の基準値未満であるならば正常とし、基準値を超える分布データが存在する場合には、劣化が生じていると判定する。
【００４０】
以上の内容を整理すると、傾き情報に関する確率密度分布を利用して、構造物の取付状態の劣化診断を実施するに当たっては、以下の処理を行うことを特徴としている。
［特徴１］加速度情報の低周波成分を抽出することで、傾き情報を取得する。
［特徴２］傾き情報のサンプリングレートは、情報が欠損しない程度に遅くでき、例えば、５秒間隔でデータ収集を行うことができる。
［特徴３］取り付け方向の修正がなされた傾きデータをデシベル化することで、記憶容量を削減した上で、効率的な量子化を行っている。
［特徴４］所定のデシベル単位に基づいて、学習時および劣化診断時の確率密度分布を算出し、３軸の確率密度比のユークリッド距離としてベクトル長を求め、劣化診断のための評価値としている。
［特徴５］基準値を超える評価値がある場合には、センサヘッド２０が設置された部分で、構造物の傾きが許容できないレベルに達していると判断し、劣化状態を定量的に判断する。
【００４１】
（２）第２の劣化診断部１２による固有振動数情報に基づく劣化診断の処理の流れについて
次に、第２の劣化診断部１２による具体的な処理内容を説明する。
図６は、本発明の実施の形態１における第２の劣化診断部１２による劣化診断処理の流れを示すフローチャートである。この図６のフローチャートに基づいて、固有振動数情報に基づいた第２の劣化診断部１２による具体的な劣化診断処理について説明する。
【００４２】
まず始めに、ステップＳ６０１において、第２の劣化診断部１２は、センサヘッド２０内の加速度情報出力部２２から取得した加速度情報に基づいて、振動数情報の抽出を行う。センサヘッド２０からは、例えば、５０ＨｚのサンプリングレートでＸ、Ｙ、Ｚ軸の加速度情報が出力されている。しかしながら、常に大きな振動が発生している訳ではない。一方、周波数解析を行うに当たっては、高いサンプリングレートが必要であり、常時、データを取得するためには、膨大な記憶領域を要する。
【００４３】
そこで、本実施の形態１における第２の劣化診断部１２は、取得した３軸の加速度情報のうち、少なくとも１軸において、あらかじめ設定した基準レベルを超えた加速度が得られた場合には、過去から未来に亘るデータを２ｎの一定数（例えば５１２個）取得している。
【００４４】
図７は、本発明の実施の形態１における３軸分の振動数情報の収集に関する具体例を示した図である。図７においては、Ｘ軸において測定された加速度データが、基準レベルを超えた際に、その時点をイベントトリガとして、３軸分それぞれに過去データと未来データからなる一定数のデータを記憶する場合を例示している。
【００４５】
次に、第２の劣化診断部１２は、ステップＳ６０２において、センサヘッド２０の取り付け方向の修正を行う。具体的には、第２の劣化診断部１２は、センサヘッド２０の取り付け方向のキャリブレーションデータ（アフィン変換の行列の要素）と、先のステップＳ６０１で得られた一定数のデータのそれぞれを行列演算することで、それぞれのセンサヘッド２０の座標軸を一致させるように、取り付け方向の誤差修正を行う。
【００４６】
次に、第２の劣化診断部１２は、ステップＳ６０３において、誤差修正された３軸データに対し、ＦＦＴ処理によりパワースペクトルを算出し、固有振動数を求める。図８は、本発明の実施の形態１における３軸分のパワースペクトルの算出結果を示した図である。加速度データには、直流成分（傾き成分）が含まれる。そこで、第２の劣化診断部１２は、パワースペクトルから低周波成分を取り除き、残ったスペクトルの中から最大の極大値（図８中のｆ１ｘ、ｆ１ｙ、ｆ１ｚに相当）、および第２の極大値（図８中のｆ２ｘ、ｆ２ｙ、ｆ２ｚに相当）を求め、これをそれぞれ、各軸の第１固有振動数、第２固有振動数とする。
【００４７】
次に、第２の劣化診断部１２は、ステップＳ６０４において、第１固有振動数と第２固有振動数をそれぞれ直交軸として、２次元の累積度数分布を作成し、これに局所平滑化処理を施すことで、２次元の確率密度分布を生成する。具体的には、第２の劣化診断部１２は、例えば、先のステップＳ６０１でのイベントトリガごとに、先のステップＳ６０３で得られた第１固有振動数および第２固有振動数で規定される２次元の累積度数分布を作成し、この累積度数分布に局所平滑化処理を施すことで、確率密度分布を生成する。なお、平滑化処理としては、例えば、ガウス関数を適用することができる。
【００４８】
なお、設置時の初期学習データとしての確率密度分布は、例えば、１ケ月程度のデータを用いて生成される。一方、劣化診断時の現在のデータとして作成する確率密度分布のデータ数は、例えば、過去３時間〜２４時間程度のデータに基づいて生成される。
【００４９】
また、基準データとなる正常時における学習データとしては、必ずしも設置時に取得したデータを継続して使用する必要はない。長期にわたって劣化診断を行うため、劣化診断を行う過程で、継続的に学習データを更新することも可能である。劣化監視時に継続学習データとして採用する確率密度分布は、例えば、過去数時間〜数週間程度のデータに基づいて生成される。
【００５０】
次に、第２の劣化診断部１２は、ステップＳ６０５において、劣化診断時に求めた各軸の２次元確率密度分布（診断時データに相当）を、学習時に求めた各軸の２次元確率密度分布（基準データに相当）で割ることで、両者の比率を計算し、それぞれの軸について確率密度比の２次元分布データを生成する。なお、基準データとしては、初期設定時のデータ以外に、継続学習データにより更新された学習データを採用することもできる。
【００５１】
図９は、本発明の実施の形態１における、ある１軸としてＸ軸について算出された学習時の２次元確率密度分布、劣化診断時の２次元確率密度分布、および両者の比である２次元の確率密度比の具体例を示した図である。図９に示したように、固有振動数に関する確率密度分布は、第１固有振動数および第２固有振動数で規定される２次元のデータとして生成されるが、確率密度比の考え方は、第１の劣化診断部１１による傾きデータに関する１次元の確率密度比の考え方と同じである。第２の劣化診断部１２は、各軸について算出された２つの固有振動数に関する２次元の確率密度比を、劣化診断のための評価値とする。
【００５２】
そして、最後に、ステップＳ６０６において、第２の劣化診断部１２は、３軸の確率密度比が、ともに劣化判定の基準値未満であるならば正常とし、いずれか１つでも基準値を超える確率密度比の２次元分布が存在する場合には、劣化が生じていると判定する。
【００５３】
以上の内容を整理すると、固有振動数情報に関する確率密度分布を利用して、構造物の取付状態の劣化診断を実施するに当たっては、以下の処理を行うことを特徴としている。
［特徴１］加速度情報に関して、いずれかの軸で基準レベルを越えた加速度データが得られた時点を含む一定数の加速度データを抽出することで、固有振動数情報を取得する。
［特徴２］固有振動数情報は、その後の周波数解析に用いられるため、センサヘッド２０からの出力レートをダウンサンプリングすることは、適切でない。そこで、いずれかの軸で基準レベルを越えた加速度データが得られた時点をイベントトリガとして、過去データと未来データからなる一定数の加速度データを抽出して固有振動数情報を記憶させることで、記憶容量の低減を図っている。
［特徴３］取り付け方向の修正がなされた固有振動数情報について周波数解析を行ってパワースペクトルを算出し、パワースペクトルから低周波成分を取り除いた残りのスペクトルの中から、最大の極大値および第２の極大値を求めることで、各軸の第１固有振動数、第２固有振動数を算出する。
［特徴４］イベントトリガごとに算出された第１固有振動数、第２固有振動数による２次元の直交座標として作成した各軸の累積度数分布から、学習時および劣化診断時の２次元の確率密度分布を算出し、３軸の確率密度比を、劣化診断のための評価値とする。
［特徴５］基準値を超える評価値がある場合には、センサヘッド２０が設置された部分で、構造物の振動が許容できないレベルに達していると判断し、劣化状態を定量的に判断する。
【００５４】
以上のように、実施の形態１によれば、劣化診断対象の構造物から得られた加速度情報に基づいて、傾きに関する第１の特徴量および固有振動数に関する第２の特徴量を抽出している。そして、それぞれの特徴量に関して、正常時の基準データに相当する学習時の確率密度分布と、劣化診断時の測定結果に基づく確率密度分布との比較により、有意差が検出された場合には、劣化が発生していると判断している。この結果、構造物の取付状態の劣化診断を、長期にわたって定量的に実施することを可能としている。
【００５５】
なお、上述した実施の形態１では、傾きに関する第１の特徴量に基づく劣化診断と、固有振動数に関する第２の特徴量に基づく劣化診断について説明したが、これら２つの診断は、いずれか１つのみを行うことによっても、劣化診断を長期にわたって定量的に実施することが可能である。
【００５６】
また、上述した実施の形態１による劣化診断は、最小限の構成として、センサヘッドを１個用いた場合にも、劣化診断が可能である。ただし、センサヘッドを設置した箇所では劣化が発生していなくても、他の場所で劣化が発生しているおそれはある。そこで、センサヘッドを複数箇所に設置し、いずれかのセンサヘッドで劣化状態が検出されたときに、構造物の取付劣化が発生したと判断することで、検出精度の向上が期待できる。
【００５７】
さらに、２つのセンサヘッドのデータを活用できる場合には、個々の診断結果に加え、傾きに関しては、２つのセンサヘッドの傾きの差として、振動数に関しては、２つのセンサヘッドの振動の位相差として、確率密度分布に基づく劣化診断を行うことができ、さらなる検出精度の向上が期待できる。なお、振動数に関しては、位相差として行うだけでなく、振幅差として同様に劣化診断を行うこともできる。
【００５８】
また、上述した実施の形態１においては、トンネル内の吊り下げ構造物、あるいは張り出し構造物を劣化診断対象の一例として挙げたが、本発明は、これに限定されない。固定面に取り付けられ、経年的に取付状態が変化してしまうおそれのある構造物であれば、長期にわたって定量的に劣化診断を行うことができる。また、既存の構造物に対して、センサヘッドを後付けすることによっても、センサヘッドの設置以降において、構造物の取付状態の経年的変化を、定量的に診断することができる。
【００５９】
また、上述した実施の形態１においては、正常時の基準データとして、初期段階で学習する場合と、劣化診断時における学習により更新する場合について説明したが、本発明は、これに限定されない。正常であることを判断するための確率密度分布は、唯一である必要はなく、例えば、時間帯毎に個別の基準データを設ける、あるいは構造物に発生する事象毎に別個の基準データを設けることもできる。さらに、経年変化を考慮して、初期段階での確率密度分布と、劣化診断時の学習により得られた新たな確率密度分布を併用することも可能である。
【符号の説明】
【００６０】
１０  センサコントローラ、１１  第１の劣化診断部、１２  第２の劣化診断部、２０  センサヘッド、２１  センサ部、２２  加速度情報出力部。
【図１】
2014228472.tif 000003
【図２】
2014228472.tif 000004
【図３】
2014228472.tif 000005
【図４】
2014228472.tif 000006
【図５】
2014228472.tif 000007
【図６】
2014228472.tif 000008
【図７】
2014228472.tif 000009
【図８】
2014228472.tif 000010
【図９】
2014228472.tif 000011