JP2003043154A

JP2003043154A - 振動波形弁別方法および加速度データ出力装置並びに液状化検知装置

Info

Publication number: JP2003043154A
Application number: JP2001226067A
Authority: JP
Inventors: Hiroyuki Furukawa; 洋之古川; Hikari Takubo; 光田久保; Takashi Yanada; 貴簗田; Yoshihisa Shimizu; 善久清水; Wataru Nakayama; 渉中山; Kenichi Koganemaru; 健一小金丸
Original assignee: Tokyo Gas Co Ltd; Azbil Corp
Current assignee: Tokyo Gas Co Ltd; Azbil Corp
Priority date: 2001-07-26
Filing date: 2001-07-26
Publication date: 2003-02-13
Anticipated expiration: 2021-07-26
Also published as: JP3652629B2

Abstract

(57)【要約】【課題】本来の入力信号（加速度信号）にノイズ成分
が重畳した場合に、ノイズ成分をリアルタイムに弁別す
ることができる振動波形弁別方法、およびこの方法を用
いて、ノイズ成分の影響のない加速度データを出力でき
る加速度データ出力装置を提供する。【解決手段】加速度素子１０１に入力された加速度信
号に対し、ＣＰＵ１０２の振動波形弁別手段１０２Ａに
より正負の最大値を求め、最大値の絶対値の比が所定の
範囲にあるか否かを判断する。この比が所定の範囲にな
く、この状態が少なくとも所定の時間継続した場合に
は、入力信号はノイズであると判断する。この比が所定
の範囲内である場合または比が所定の範囲にない状態が
所定の時間継続しなかった場合には、入力信号が地震波
形であると判断し、出力手段１０２Ｂにおいて最大加速
度Ａ_maX、ＳＩ値を算出し、出力する。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は、振動波形とノイズ
成分とを弁別するための振動波形弁別方法、および加速
度素子から出力される加速度信号を基にして地震動の被
害推定等のための加速度データを出力する加速度データ
出力装置、並びにこの加速度データ出力装置から出力さ
れる加速度データを基に、地震動の大きさの指標である
ＳＩ（Spectrum Intensity）値を演算出力し、地震動に
より液状化現象が発生したか否かを検知するための液状
化検知装置に関する。

【０００２】

【従来の技術】従来、この種のＳＩ値演算出力装置や液
状化検知装置など、加速度センサからの出力を基に液状
化判定等の演算判定処理を行い処理結果を出力する機器
においては、加速度素子が出力する信号に加速度信号以
外のノイズ信号が重畳していることがある。この種の装
置では、ノイズ信号を含んだ信号を用いて演算処理を行
うと、適正な演算出力や液状化判定が得られない場合が
ある。

【０００３】例えば、加速度センサのアナログ信号処理
過程において整流増幅回路等が存在する場合には、外部
からの電磁波ノイズが侵入すると、ノイズ自身は交流成
分であるにもかかわらず、整流増幅回路によって片側
（正負の一方）のみの成分を有する信号として出力され
る。

【０００４】図１０は、実際の静電容量式の加速度セン
サが数十ＭＨｚの電磁波ノイズに曝されたときに得られ
た出力波形である。出力波形は、ノイズが注入されてい
る間、ノイズ未注入の出力レベルに対して片側（単一）
方向へドリフトしている。これは、内部処理回路中のあ
る部分で注入ノイズが積分されているためと考えられ
る。図１１は実際に屋外で観測された電磁波ノイズと思
われる出力波形であり、電磁波を発生しながら移動体
（例えばトラックなどの車両）が機器周辺を通過する際
に、移動体が電波を放射しながら機器に近づき、遠ざか
ることによってこのような山形（三角波に近い）波形が
観測されると考えられる。

【０００５】また、別の例を挙げるならば、地震動を計
測すべく設計された低周波応答用の加速度センサに対し
て衝撃が加えられたときにも同様の現象が発生する。図
１２（Ａ），（Ｂ）は低周波加速度検出用に設計された
圧電式加速度センサに衝撃を加えたときの出力波形であ
り、衝撃により大きな電荷が一瞬にしてデバイス両端に
発生し処理回路側の入力インピーダンスが周波検出用に
非常に高く設計されているために出力が片側に振り切
れ、もとの出力レベルに戻るまでに非常に長い時間を必
要としていることがわかる。なお、同図（Ｂ）は同図
（Ａ）の部分拡大図である。また、図１３は容量式加速
度センサに衝撃を加えたときの出力波形であるが、この
ようなインパルス状の波形が観測される場合がある。

【０００６】

【発明が解決しようとする課題】従来、このようなノイ
ズの影響を排除するための技術としては、信号波形がゼ
ロ（基軸）を通過する時間間隔を計測し、これが期待さ
れる周波数帯域の振動周期を越える場合に、この信号は
ノイズであると判定する方法が採用されている。しかし
ながら、この方法では、図１１に示したノイズ信号や図
１３に示したノイズ信号に対しての効果は期待できない
という問題があった。

【０００７】本発明はかかる問題点に鑑みてなされたも
ので、その第１の目的は、本来の入力信号（加速度信
号）にノイズ成分が重畳した場合に、ノイズ成分を弁別
することができる振動波形弁別方法、およびこの方法を
用いて、ノイズ成分の影響のない加速度データを出力で
きる加速度データ出力装置を提供することにある。

【０００８】本発明の第２の目的は、上記加速度データ
出力装置を用いてＳＩ値演算や液状化現象を精度良く検
知することができる液状化検知装置を提供することにあ
る。

【０００９】

【発明を解決するための手段】請求項１記載の振動波形
弁別方法は、振動波形として加速度信号を出力する加速
度素子と、加速度信号を基にした加速度データ、若しく
は前記加速度データに演算処理を行なった結果を出力す
る出力手段とを含む加速度データ出力装置において、加
速度素子から出力される加速度信号に対して、その正負
の最大値の絶対値の比が所定の範囲内にあるか否かを検
出し、加速度信号の正負の最大値の絶対値の比が所定の
範囲内にない場合には、加速度信号が入力値として期待
されている信号ではないと判断するものである。

【００１０】請求項２記載の振動波形弁別方法は、請求
項１に記載の方法において、加速度信号の正負の最大値
の絶対値の比が所定の範囲内にない状態が所定の時間幅
以上にわたって継続した場合に、加速度信号が入力値と
して期待されている信号ではないと判断するものであ
る。

【００１１】請求項３記載の振動波形弁別方法は、請求
項１または２に記載の方法において、加速度信号を時系
列的な連続波形として処理し、所定の時間幅毎に前記加
速度信号の正負の最大値の絶対値の比を算出するもので
ある。

【００１２】請求項４記載の振動波形弁別方法は、請求
項１乃至４のいずれか１項に記載の方法において、加速
度信号に基線補正を行った後、その正負の最大値の絶対
値の比を算出するものである。

【００１３】また、請求項５記載の加速度データ出力装
置は、加速度信号を振動波形として出力する加速度素子
と、この加速度素子から出力される加速度信号を基にし
て基にした加速度データ、若しくは加速度データに演算
処理を行なった結果を出力する出力手段と、請求項１乃
至４のいずれか１項に記載の振動波形弁別方法を用いて
判断を行う振動波形弁別手段とを備えている。

【００１４】請求項６記載の加速度データ出力装置は、
請求項５記載の出力装置において、出力手段は、振動波
形弁別手段が、加速度信号は入力値として期待されてい
る信号ではないと判断した場合に、加速度データの演算
処理または出力を停止するように構成したものである。

【００１５】請求項７記載の加速度データ出力装置は、
請求項５または６記載の出力装置において、出力手段
は、地震波による振動波形を加速度信号としてその最大
加速度およびＳＩ値を出力するように構成したものであ
る。

【００１６】更に、請求項８記載の液状化検知装置は、
請求項５乃至７のいずれか１項に記載の加速度データ出
力装置と、この加速度データ出力装置から出力された加
速度データを基に、地震動により液状化現象が発生した
か否かを判定する液状化判定装置とを備えるように構成
したものである。

【００１７】

【発明の実施の形態】以下、本発明の実施の形態につい
て図面を参照して詳細に説明する。

【００１８】具体的な実施の形態の説明の前に、本発明
に係る振動波形弁別方法の基本的な原理について説明す
る。

【００１９】一般に、前述の説明（図１０乃至図１３）
からも明らかなように、加速度センサの電磁ノイズに対
する挙動は、加速度信号において正（＋）方向へのＤＣ
（直流）バイアス的なドリフトとなって観測されるが、
これに対し、実際の地震波による加速度波形は、必ず正
（＋）側と負（−）側に成分を有するＡＣ（交流）的な
挙動を示す。このようなことから、本発明においては、
加速度素子の出力が、ＤＣ的な挙動を示しているか、Ａ
Ｃ的な挙動を示しているかによって、地震の振動波とノ
イズとを弁別するものである。

【００２０】これを図５（Ａ），（Ｂ）の波形図を参照
して説明する。図５（Ａ）は走行車両により電磁波ノイ
ズが注入された場合の典型的な加速度素子の出力を、同
図（Ｂ）は実際の地震波が入力された場合の加速度素子
の出力波形をそれぞれ示している。図からも明確なよう
に、ノイズ波の場合と実際の地震波の場合では基線に対
する出力波形の正側方向の最大値α+ と負側方向の最大
値α- の絶対値の比（｜α+ ｜／｜α- ｜）が明らかに
異なっている。本発明の振動波形弁別方法は、これに着
眼したものであり、次式（１）に示したように、この絶
対値の比が、所定の範囲内にあるか否かによってノイズ
と実際の地震波とを識別する方法である。すなわち、式
（１）が成立する（真）ならば地震による加速度信号、
成立しない（否）場合にはノイズと判定するものであ
る。なお、ｒ（以下、ノイズレシオという）は、ｒ＞１
であり、例えばｒ＝５とする。

【００２１】１／ｒ＜｜α+ ｜／｜α- ｜＜ｒ …（１）

【００２２】本発明の振動波形弁別方法の基本的な考え
方は上述の通りであるが、加速度センサの実機での動作
を考慮した場合、次のような問題が生じる。すなわち、
上記方法は、図５（Ａ），（Ｂ）に示した波形からも明
らかなように、加速度素子の出力パターンを、１ノイズ
波全体あるいは１地震波全体を判定対象にした場合の考
え方である。しかし、実機においては、加速度素子の出
力は刻々と変化する量として連続的に得られるものであ
り、この出力を基にリアルタイムに判定を行い、かつ誤
判定動作を生じないような応用が必要となる。

【００２３】従って、本発明の振動波形弁別方法を実機
に適用する場合には、連続する時系列信号として得られ
る加速度素子の出力値を用いてノイズの有無の判定を２
段階に行う必要があり、これによりノイズ判定をリアル
タイムに行うことができると共に、誤判定の防止を図る
ことができる。

【００２４】次に、この方法によるリアルタイムな判定
動作について説明する。図６はノイズ信号、図７は地震
波による加速度信号に対する一例である。このうち、図
６（Ａ）はノイズ波形、同図（Ｂ）はノイズ波形の絶対
値の比｜α+ ｜／｜α- ｜についてそれぞれ経時変化を
示したものである。このようなノイズ波が演算処理部
（ＣＰＵ）へ加速度信号として時系列に沿って入力され
るとき、（Ａ）に示したように、時刻ｔa-x において正
側の最大出力量α+ はα（ax+)であり、負側の最大出力
量α- はα（ax-)である。この時点では、｜α（ax+)｜
と｜α（ax-)｜の比（｜α（ax+)｜／｜α（ax-)｜）
は、（Ｂ）に示したように式（１）の範囲内であり、ノ
イズとは判定されない。更に時刻が経過し、時刻ta-yに
なると、α+の値はα（ay+)に更新され、｜α（ay+)｜
／｜α（ax-)｜がノイズレシオ“ｒ”を超えるためにノ
イズと判定される。なお、以後のα+ の値は、ピーク値
が得られるまで増大し、より大きなピークが入力されて
こない限り、この値が最大値として固定される。その
間、α- は、α（ax-)のまま、あるいは同程度の値とな
っている。よって、｜α+ ｜／｜α- ｜は、（Ｂ）のよ
うに“ｒ”を超えて増加した後に信号波形がピークに達
した時点でほぼ一定に固定される。従って、このような
波形では、一旦ノイズと判定されれば｜α+ ｜／｜α-
｜の値が保存されている限り、それ以降の波形もノイズ
であると判断されることになる。

【００２５】図７（Ａ）は実際の地震による加速度信
号、同図（Ｂ）はその信号波形の絶対値の比｜α+ ｜／
｜α- ｜の経時変化をそれぞれ示している。このような
加速度信号がＣＰＵに時系列として入力されるとき、
（Ａ）に示したように時刻ｔb-xにおいてα+ はα（bx
+)であり、α- はα（bx-)である。この時点では、
（Ｂ）に示したように｜α（bx+)｜／｜α（bx-)｜は式
（１）の範囲内であり、ノイズと判定されない。ところ
が、時刻tb-yとなると、α+ はα（by+)に更新される。
このとき、｜α（by+)｜／｜α（bx-)｜はノイズレシオ
“ｒ”を超えてしまう。更に、時刻tb-zでは、α- はα
（bz-)に更新され、｜α（by+)｜／｜α（bz-)｜はノイ
ズレシオ“ｒ”以下となる。その後、時刻tb-y' では、
α+ はα（by'+)に更新され、｜α（by'+) ｜／｜α（b
z-)｜はノイズレシオ“ｒ”を再度超えることになる。
時刻tb-z' では、α- はα（bz'-) に更新され、この時
点で｜α（by+)｜／｜α（bz'-) ｜はノイズレシオ
“ｒ”以下となる。

【００２６】このように、先の基本的な考え方をそのま
まリアルタイムにノイズ判定を行う連続した時系列に適
用すると、地震波の振動状態によってはノイズと誤認さ
れる場合がある。それは出力波形の包絡線が基線に対し
大きく変動する場合であり、特に、無振動状態から振動
が始まる地震の初期振動において頻繁に現れることがみ
てとれる。それでも、時間が経過して時刻tb-z' となる
とα- はα（bz'-) に更新され、この時刻で｜α（by+)
｜／｜α（bz-)｜は式（１）の範囲内となる。これ以
後、｜α+ ｜／｜α- ｜は式（１）の範囲内にあり、正
常に地震波と判定される。従って、本発明の振動波弁別
方法を実機に適用する場合には、地震動の初期振動に生
じ易いこの種の誤認を回避するため、先の基本的な方法
に対し下記の論理動作を付加する必要がある。

【００２７】判定時間幅“Ｔｒ”を設定し、加速度出力
の正側方向の最大値α+ と負側方向の最大値α- の絶対
値の比がノイズレシオｒを超えた状態が“Ｔｒ”時間以
上継続した場合に、入力信号はノイズであると判定す
る。ここで、判定時間幅“Ｔｒ”は、地震動による加速
度波形の基本周期変動を包含する大きさ（例えば２秒）
に設定する。

【００２８】本発明の振動波弁別方法では、この論理動
作を付加することにより、図６に示したノイズ波形では
ノイズレシオ“ｒ”を超えた時刻ta-yから“Ｔｒ”経っ
た時刻ta-jになったときにノイズ判定が行われる。ま
た、図７に示した加速度波形では、時刻tb-yから“Ｔ
ｒ”後の時刻tb-j、および時刻tb-y' から“Ｔｒ”後の
時刻tb-j' のどちらの時点においても、｜α+ ｜／｜α
- ｜は式（１）の範囲内に収まり、誤ってノイズと判定
されることはない。

【００２９】次に、本発明の振動波弁別方法において、
加速度波形の正側の最大値α+ および負側の最大値α-
の絶対値を具体的に求める方法について説明する。その
方法としては、図８に示した固定幅時間窓法および図９
に示したフレックス時間窓法がある。なお、ここでは、
加速度の最大値の絶対値｜α+ ｜，｜α- ｜は、速度応
答の最大値Ｓｖに対応している。

【００３０】まず、固定幅時間窓法について説明する。
この固定幅時間窓法は、図８（Ａ）〜（Ｃ）に示したよ
うに、単純に固定幅の時間窓を演算間隔Δｔずつ移動さ
せる方法である。ここに、図８（Ａ）は、地震波のある
軸方向の加速度に対する速度応答ｖ（ｔ）を表してい
る。Ｓｖは速度応答の最大値を意味するが、固定幅の時
間窓を現時刻からｎ秒前までと設定し、時間窓内のｖ
（ｔ）を使用してＳｖを求める。演算間隔はΔｔで規定
する。図８（Ｂ）に示したように、ｎ＝４〔ｓ〕，Δｔ
＝１〔ｓ〕の条件下で、現時刻ｔ５としたときのＳｖ
（ｔ５）を求める場合には、ｔ２〜ｔ５の時間窓にな
り、ｖ（ｔ２），ｖ（ｔ３），ｖ（ｔ４），ｖ（ｔ５）
の最大値がＳｖ（ｔ５）になる。次の時刻ｔ６のＳｖ
（ｔ６）を求める場合には、ｔ３〜ｔ６の時間窓となる
ので、ｖ（ｔ２）は忘却され、代わりにｖ（ｔ６）が記
憶されてＳｖ（ｔ６）が求まる。図８（Ｃ）は、求めら
れたＳｖの時刻変化を表している。

【００３１】この固定幅時間窓法では、時間窓内の各時
刻でのｖ（ｔ）を全て記憶しておく必要がある。実際に
は、地震波の加速度信号からＳＩ値を求めるため、Δｔ
＝０．０１〔ｓ〕とする必要があり、上記例に対して信
号データを格納するメモリは１００倍となる。また、Ｓ
ｖを求める度に時間窓内の全てのｖ（ｔ）を比較する必
要があり、演算処理部での演算負荷が重くなるという問
題がある。

【００３２】フレックス時間窓法は、上述の固定幅時間
窓法の問題を解決するための方法である。以下、図９
（Ａ）〜（Ｄ）を参照してこの方法について説明する。

【００３３】このフレックス時間窓法では、使用する時
間窓を、１つの固定幅の記憶時間窓Ｗ1 と隣接して伸び
縮みする演算時間窓Ｗ0 との合計２つの時間窓で構成す
る。演算時間窓Ｗ0 はサンプリング時間とともに伸びて
ゆき、その間、サンプリングの度にメモリＭ0 に格納さ
れた時間窓Ｗ0 におけるｖ（ｔ）の最大値が更新され
る。なお、演算時間窓Ｗ0 が規定の幅となると、メモリ
Ｍ0 の値はメモリＭ1 にコピーされてメモリＭ0 はゼロ
クリアされる。

【００３４】図９（Ａ）は、加速度素子の出力におい
て、ある固有周期の速度応答ｖ（ｔ）の特性を表してい
る。図９（Ｂ）は、２つの時間窓Ｗ1 ，Ｗ0 で構成され
る時間窓の幅を各サンプリング時間毎（ｔ＝ｔ０，ｔ
１，ｔ２，…）に示している。図９（Ｃ）は、記憶時間
窓Ｗ1 ，演算時間窓Ｗ0 のそれぞれに対応してｖ（ｔ）
の最大値を格納するメモリＭ1 ，Ｍ0 を表している。ま
た、図９（Ｄ）は、Ｓｖの経時変化を表している。

【００３５】演算は時刻ｔ０から開始し、このとき全変
数は０とする。各時刻毎にｖ（ｔ）をメモリＭ0 の値と
比較し、ｖ（ｔ）の方が大きいならばＭ0 の値をｖ
（ｔ）で更新する。時刻ｔ４でＷ0 の値はＷ１としてＭ
１にコピーされ、Ｗ0 の値はゼロクリアされる。このと
き、Ｓｖ（ｔ４）はｔ＝ｔ０〜ｔ４までのｖ（ｔ）の最
大値となる。時刻の経過と共に時間窓の幅は大きくな
り、時刻ｔ７でＷ1 とＷ0 を合わせた時間窓の幅は最大
（ｔ０〜ｔ７）になる。

【００３６】時刻ｔ８でＷ0 の値をＷ1 の値へのコピー
した後、Ｗ0 の値のゼロクリア処理が行われ、以後同様
の処理が行われる。以上の演算手法を用いることによ
り、メモリには時間窓Ｗ0 ，Ｗ1 に対応したｖ（ｔ）の
最大値を２つ格納するだけで済み、これらを比較するこ
とでその時点でのＳｖが求められる。従って、固定時間
窓法のように時間窓内全てのｖ（ｔ−ｎ）〜ｖ（ｔ）を
記憶する必要がなくなりメモリ消費量を大幅に節約する
ことができる。また、Ｓｖを求めるに当たり、フレック
ス時間窓法ではＭ0 の値と時刻ｔ毎のｖ（ｔ）、および
Ｍ1 の値の比較を行なうだけでよいので、固定時間窓法
のように、時間窓内全てのｖ（ｔ−ｎ）〜ｖ（ｔ）同士
を比較して最大値を求める必要がなくなり、そのため演
算処理部の処理負荷が大幅に軽減される。これにより、
フレックス時間窓法を用いて加速度の最大値の絶対値｜
α+ ｜，｜α- ｜の演算をリアルタイムに行うことがで
きる。

【００３７】なお、これらの方法では、時間窓あたりに
格納可能なデータ数と窓数によって決定される所定時間
後には、保持されていた最大加速度α+ ，α- の値は、
自動的に消去される。また、α+ ，α- として一旦保持
された値が以降のノイズ検出に影響することを防ぐため
に、｜α+ ｜／｜α- ｜が所定期間以上にわたって所定
の閾値以下である場合には現時点でノイズは入力されて
いないものと判断し、過去のものとなったα+ ，α- の
値をクリアするようにしてもよい。

【００３８】次に、以上のような基本原理に基づいた振
動波形弁別方法を適用した具体的な実施の形態について
説明する。

【００３９】図１は本発明の一実施の形態に係る液状化
検知装置１０のブロック構成を表すものである。この液
状化検知装置１０は、加速度データ出力装置１００およ
び液状化判定装置２００を備えている。

【００４０】加速度データ出力装置１００は、例えばＳ
Ｉセンサと称されるものである。ＳＩセンサは、地震に
よる振動波形（加速度波形）を基に最大加速度Ａ_maxを
検出すると共に、加速度波形の速度応答を基にしてＳＩ
値を得るものである。最大加速度Ａ_maxは、例えば、地
震波形が基線を横切る時間（ゼロクロス時間）間隔を基
にして推定されるゼロクロス周期などの固有周期Ｔ内に
おける加速度の絶対値の最大値である。ＳＩ値は、地震
による一般的な建物の揺れの程度を数値化したものであ
り、地震の加速度を入力として応答解析を行うことによ
り求められた特定の固有周波数範囲における平均値であ
る。このＳＩ値は、具体的には、以下の式（２）により
求めることができ、過去の観測結果から最大加速度より
も、震度や地震被害の有無との相関関係が高いことが知
られている。

【００４１】

【数１】

【００４２】この加速度データ出力装置１００は、地震
等によって加速度信号を出力する加速度素子１０１と、
この加速度素子１０１からの加速度信号を受けて、演算
処理を行うＣＰＵ（Central Processor Unit；中央演算
装置) １０２とを有している。

【００４３】加速度素子１０１としては例えば静電容量
型のものがある。この静電容量型の加速度素子は、固定
の上部電極と、重りのついた可動の下部電極とを対向配
置した構造を有し、加速度に応じて重りと共に下部電極
が動くようになっている。このときの上部電極と下部電
極との間の静電容量の変化を読み取ることにより、加速
度の方向および大きさを検出することができる。なお、
加速度素子１０１としては、その他、ピエゾ抵抗型ある
いは圧電型などの半導体加速度素子であってもよい。

【００４４】ＣＰＵ１０２は、図示しないＲＯＭ（Read
Only Memory）等に格納されたプログラムに従って演算
処理を実行するものであり、やはり図示はしないが、演
算のためのデータを格納するメモリ（Ｍ1 ，Ｍ0 ）を有
している。このＣＰＵ１０２は、入力値として期待され
ている地震による加速度信号とノイズ信号とを弁別する
ための振動波形弁別手段１０２Ａ、および加速度素子１
０１から出力される加速度信号を基にして加速度データ
を出力するための出力手段１０２Ｂとしての機能を備え
ている。

【００４５】振動波形弁別手段１０２Ａは、加速度素子
１０１から出力される加速度信号の正負の最大値の絶対
値の比が所定の範囲内にあるか否かを検出し、その絶対
値の比が、前述の式（１）に示した所定の範囲になく、
かつその状態が所定の時間（判定時間幅Ｔｒ）以上継続
した場合には、計測されている信号は、入力値として期
待されている加速度信号ではない、すなわち、前述のよ
うなノイズが発生したと判断するものである。

【００４６】出力手段１０２Ｂは、例えば、加速度素子
１０１から出力される加速度信号を基に最大加速度Ａ
_maXを検出すると共にＳＩ値を求め、これら最大加速度
Ａ_maXおよびＳＩ値を加速度データとして出力するもの
である。

【００４７】また、本実施の形態では、加速度素子１０
１とＣＰＵ１０２との間には基線補正用のフィルタ１０
３、増幅器１０３およびＡ／Ｄ（アナログ／ディジタ
ル）変換器１０５が介挿されている。基線補正用のフィ
ルタ１０３は、加速度素子１０１から出力された加速度
信号のＤＣバイアス成分を零に補正するもので、これに
より振動波形弁別手段１０２Ａにおいて加速度信号の正
負それぞれの最大値が求められる。増幅器１０４は、フ
ィルタ１０３から出力された加速度信号を所定のレベル
に増幅するものであり、Ａ／Ｄ変換器１０５は増幅され
た加速度信号（アナログ信号）をディジタル信号に変換
してＣＰＵ１０２へ出力するものである。

【００４８】液状化判定装置２００は、加速度データ出
力装置１００から出力された最大加速度Ａ_maXおよびＳ
Ｉ値を基にして液状化現象が発生したか否かの判断を行
うものである。具体的には、例えば、これら最大加速度
Ａ_maXおよびＳＩ値に加え、地震波形の固有周期Ｔおよ
び地盤変位Ｄ_maXを求め、これら４つのパラメータが予
め定めた閾値以上である場合に、液状化現象が発生した
と判断する。なお、地盤変位Ｄ_maXは例えば、最大加速
度Ａ_maXおよびＳＩ値を用いて次式（３）によって近似
的に算出することができる。

【００４９】Ｄ_maX＝２×（ＳＩ）²／Ａ_maX…（３）

【００５０】なお、このような構成の液状化検知装置１
では、加速度データ出力装置１００において、衝撃や電
磁波の影響によって、フィルタ１０３、増幅器１０４、
Ａ／Ｄ変換器１０５等に対してノイズ信号が入力され、
このノイズ信号がＣＰＵ１０２に入力してしまうことが
ある。そのため、従来では、入力として期待されている
実際の地震波による加速度信号に基づく正確なＳＩ値を
出力することができず、液状化判定装置において、液状
化現象の判定を精度良く行うことができなかったことは
前述の通りである。これに対して、本実施の形態では、
以下に説明するように、加速度データ出力装置１００の
振動波形弁別手段１０２Ａによってノイズの影響をなく
すことが可能になる。

【００５１】次に、図２〜図４の流れ図を参照して、本
実施の形態の液状化検知装置１０の作用について説明す
る。

【００５２】〔正常判定処理〕まず、地震が発生する
と、加速度データ出力装置１００においては、加速度素
子１０１から加速度信号が出力され、この加速度信号は
フィルタ１０３により基線補正がなされ、増幅器１０４
により所定のレベルに増幅された後、Ａ／Ｄ変換器１０
５においてディジタル信号に変換される。このディジタ
ル信号に変換された加速度信号がＣＰＵ１０２へ入力さ
れる（ステップＳ１００）。加速度信号が入力される
と、ＣＰＵ１０２では、振動波形弁別手段１０２Ａによ
り、前述のフレックス時間窓法によって、その加速度信
号の正負の最大値を求める（ステップＳ１０１）。この
時点で“ノイズ確定状態”にない場合（ステップＳ１０
２；Ｎ）には、これら最大値の絶対値の比が（１）式で
示した所定の範囲にあるか否かを判断する（ステップＳ
１０３）。“ノイズ確定状態”とは、後述するように、
先に入力された一連の加速度信号を基にして、既にノイ
ズと判定されている状態をいう。

【００５３】最大値の絶対値の比が（１）式で示した範
囲である場合（ステップＳ１０３；Ｙ）には、入力信号
が入力値として期待されている加速度信号であると判断
し、“出力一時停止状態”をクリアし（ステップＳ１０
４）、出力手段１０２Ｂにおいて最大加速度Ａ_maXを算
出し（ステップＳ１０５）、ＳＩ値を算出する（ステッ
プＳ１０６）。続いて、現在“出力一時停止状態”であ
るか否かを判断するが（ステップ１０７）、いまは“出
力一時停止状態”にない（Ｎ）ため、出力手段１０２Ｂ
は、算出した最大加速度Ａ_maXおよびＳＩ値を出力する
（ステップＳ１０８）。更に、液状化判定装置２００で
は、これら最大加速度Ａ_maXおよびＳＩ値に加え、地震
波形の固有周期Ｔおよび地盤変位Ｄ_maXを算出し、これ
らの値に基づいて液状化現象が発生したか否かの判定を
行う（ステップＳ１０９）。

【００５４】〔ノイズ検出処理〕一方、最大値の絶対値
の比が（１）式で示した範囲にない場合（ステップＳ１
０３；Ｎ）に、更にこの状態が判定時間Ｔｒ継続してい
ない場合（ステップＳ１１０；Ｎ）は、入力信号が入力
値として期待されている加速度信号であるか否かの判定
ができない。この場合には、“出力一時停止状態”とし
たうえで（ステップＳ１１１）、最大加速度Ａ_maX，Ｓ
Ｉ値を算出する（ステップＳ１０５，Ｓ１０６）。続い
て“出力一時停止状態”であるか否かを判断するが（ス
テップ１０７）、いまは“出力一時停止状態”である
（Ｙ）から、最大加速度Ａ_maX，ＳＩ値の演算結果を出
力することなく入力信号待ち状態に戻る。なお、このの
ちに、信号がノイズによるものではないと判断された場
合には、“出力一時停止状態”はクリアされ（ステップ
Ｓ１０４）、最大加速度Ａ_maX，ＳＩ値を算出して（ス
テップＳ１０５，Ｓ１０６）、これらの値に加え、クリ
アされる以前に求められていた最大加速度Ａ_maXおよび
ＳＩ値を出力する（ステップＳ１０８）。

【００５５】入力された信号がノイズ信号である場合
は、上記ステップＳ１０３からステップＳ１０７を繰り
返した後に、最大値の絶対値の比が（１）式で示した範
囲にない状態（ステップＳ１０３；Ｎ）が判定時間Ｔｒ
継続する（ステップＳ１１０；Ｙ）ことになる。このと
き、ＣＰＵ１０２は、入力信号が入力値として期待され
ている加速度信号ではなくノイズによるものと確定し
て、現在の状態を“ノイズ確定状態”とし（ステップＳ
１１２）、Ａ_maX，ＳＩ値の演算および出力を停止し
（ステップＳ１１３）、入力信号待ち状態に戻る。

【００５６】〔ノイズクリア処理〕こうして、一旦“ノ
イズ確定状態”となった後は、次の入力信号についてフ
レックス時間窓法により、その加速度信号の正負の最大
値を求め（ステップＳ１０１）、“ノイズ確定状態”と
判断され（ステップＳ１０２；Ｙ）、一連のノイズクリ
ア処理を行う。すなわち、加速度信号の正負の最大値の
絶対値の比が（１）式で示した所定の範囲にあるか否か
を判断し（ステップＳ１１４）、絶対値の比が所定の範
囲にない場合（Ｎ）は、再び入力信号待ち状態に戻る。
絶対値の比が所定の範囲にある場合（Ｙ）には、ＳＩ値
の演算などに用いられる変数としてＣＰＵ１０２に格納
されているデータを初期化し（ステップＳ１１５）、
“ノイズ確定状態”をクリアして（ステップＳ１１
６）、再び入力信号待ち状態に戻る。なお、前述したよ
うに、フレックス時間窓法を用いて最大加速度α+ ，α
- を求める場合には、後にノイズと判定された信号の最
大加速度α+ ，α- は、フレックス時間窓法の時間窓あ
たりに格納可能なデータ数と窓数によって決定される所
定時間後には、保持されていた最大加速度α+ ，α- の
値は、自動的に消去されるため、ノイズ信号が入力され
なくなって所定時間が経過すると必ず“ノイズ確定状
態”はクリアされることになる。以上の判定動作はリア
ルタイムで連続的に行われる。

【００５７】以上のように本実施の形態の液状化検知装
置１０では、加速度データ出力装置１００において、加
速度信号の正負の最大値の絶対値の比が所定の範囲内に
あるか否かを検出し、加速度信号の正負の最大値の絶対
値の比が所定の範囲内にない状態が所定の時間継続した
場合に、計測されている入力信号はノイズであると判断
し、このときにはＳＩ値の演算または出力を停止するよ
うにしたので、液状化判定装置２００に対してノイズの
影響のないＳＩ値を出力することができる。従って、液
状化判定装置２００での判定の信頼性が向上し、地震の
被害の推定が容易になる。

【００５８】以上実施の形態を挙げて本発明を説明した
が、本発明は上記実施の形態に限定されるものではな
く、種々変形が可能である。例えば、上記実施の形態で
は、加速度データ出力装置１００から最大加速度Ａ_maX
およびＳＩ値を出力するようにしたが、いずれか一方を
出力するような構成のものとしてもよい。

【００５９】また、液状化判定装置２００の判定方法は
上記実施の形態のものに限らず、その他の判定方法を採
用してもよい。

【００６０】

【発明の効果】以上説明したように請求項１乃至５のい
ずれか１に記載の振動波形弁別方法によれば、加速度信
号の正負の最大値の絶対値の比が所定の範囲内にあるか
否かを検出し、加速度信号の正負の最大値の絶対値の比
が所定の範囲内にない場合には、計測されている入力信
号をノイズと判断するようにしたので、容易にノイズ発
生の有無を検出することができる。

【００６１】特に、請求項２記載の振動波形弁別方法に
よれば、請求項１記載の方法において、加速度信号の正
負の最大値の絶対値の比が所定の範囲内にない状態が所
定の時間継続した場合に、この加速度信号をノイズと判
断するようにしたので、より精度良くノイズの発生の有
無を検出することができる。

【００６２】また、請求項６乃至８のいずれか１に記載
の加速度データ出力装置によれば、上記振動波形弁別方
法を適用し、ノイズ成分がある場合には加速度データの
演算または出力を停止するようにしたので、ノイズの影
響のない加速度データを出力することが可能となる。

【００６３】更に、請求項９に記載の液状化検知装置に
よれば、上記加速度データ出力装置を用いるようにした
ので、ノイズの影響のない加速度データを基に、液状化
現象発生の有無の判定の信頼性が向上し、地震の被害の
推定が容易になるという効果を奏する。

【図面の簡単な説明】

【図１】本発明の一実施の形態に係る液状化検知装置の
構成を表すブロック図である。

【図２】図１に示した液状化検知装置の動作を説明する
ための流れ図である。

【図３】図２に続く流れ図である。

【図４】図２に続く流れ図である。

【図５】ノイズ波形の特性（Ａ）および地震による加速
度波形の特性（Ｂ）の違いを説明するための波形図であ
る。

【図６】本発明に係る振動波弁別方法を説明するための
波形図であり、（Ａ）はノイズ波形、（Ｂ）はこのノイ
ズ波形における正負の最大値の絶対値の比の変化状態を
説明するための図である。

【図７】本発明に係る振動波弁別方法を説明するための
波形図であり、（Ａ）は地震波形、（Ｂ）はこの地震波
形における正負の最大値の絶対値の比の変化状態を説明
するための図である。

【図８】本発明に係る振動波弁別方法において、加速度
の最大値を検出するための方法（固定幅時間窓法）を説
明するための図である。

【図９】同じく、加速度の最大値を検出するための他の
方法（フレックス時間窓法）を説明するための図であ
る。

【図１０】ノイズの一例を説明するための波形図であ
る。

【図１１】ノイズの他の例を説明するための波形図であ
る。

【図１２】ノイズの更に他の例を説明するための波形図
である。

【図１３】ノイズの更に他の例を説明するための波形図
である。

【符号の説明】

１０…液状化検知装置、１００…加速度データ出力装
置、１０１…加速度素子、１０２…ＣＰＵ、１０２Ａ…
振動波弁別手段、１０２Ｂ…出力手段、１０３…フィル
タ、１０４…増幅器、１０５…Ａ／Ｄ変換器、２００…
液状化判定装置

───────────────────────────────────────────────────── フロントページの続き (72)発明者田久保光東京都渋谷区渋谷２丁目12番19号株式会社山武内 (72)発明者簗田貴東京都渋谷区渋谷２丁目12番19号株式会社山武内 (72)発明者清水善久東京都港区海岸一丁目５番20号東京瓦斯株式会社内 (72)発明者中山渉東京都港区海岸一丁目５番20号東京瓦斯株式会社内 (72)発明者小金丸健一東京都港区海岸一丁目５番20号東京瓦斯株式会社内

Claims

【特許請求の範囲】

【請求項１】振動波形として加速度信号を出力する加
速度素子と、前記加速度信号を基にした加速度データ、
若しくは前記加速度データに演算処理を行なった結果を
出力する出力手段とを含む加速度データ出力装置におい
て、前記加速度素子から出力される加速度信号に対して、そ
の正負の最大値の絶対値の比が所定の範囲内にあるか否
かを検出し、前記加速度信号の正負の最大値の絶対値の
比が所定の範囲内にない場合には、前記加速度信号が入
力値として期待されている信号ではないと判断すること
を特徴とする振動波形弁別方法。
【請求項２】前記加速度信号の正負の最大値の絶対値
の比が所定の範囲内にない状態が所定の時間幅以上にわ
たって継続した場合に、前記加速度信号が入力値として
期待されている信号ではないと判断することを特徴とす
る請求項１記載の振動波形弁別方法。
【請求項３】前記加速度信号を時系列的な連続波形と
して処理し、所定の時間幅毎に前記加速度信号の正負の
最大値の絶対値の比を算出することを特徴とする請求項
１または２に記載の振動波形弁別方法。
【請求項４】前記加速度信号に基線補正を行った後、
その正負の最大値の絶対値の比を算出することを特徴と
する請求項１乃至３のいずれか１項に記載の振動波形弁
別方法。
【請求項５】加速度信号を振動波形として出力する加
速度素子と、この加速度素子から出力される加速度信号を基にして基
にした加速度データ、若しくは前記加速度データに演算
処理を行なった結果を出力する出力手段と、請求項１乃至４のいずれか１項に記載の振動波形弁別方
法を用いて判断を行う振動波形弁別手段とを備えたこと
を特徴とする加速度データ出力装置。
【請求項６】前記出力手段は、前記振動波形弁別手段
が、前記加速度信号は入力値として期待されている信号
ではないと判断した場合に、加速度データの演算処理ま
たは出力を停止することを特徴とする請求項５に記載の
加速度データ出力装置。
【請求項７】前記出力手段は、地震波による振動波形
を加速度信号としてその最大加速度およびＳＩ値を出力
することを特徴とする請求項５または６に記載の加速度
データ出力装置。
【請求項８】請求項５乃至７のいずれか１項に記載の
加速度データ出力装置と、この加速度データ出力装置から出力された加速度データ
を基に、地震動により液状化現象が発生したか否かを判
定する液状化判定装置とを備えたことを特徴とする液状
化検知装置。