JP3469509B2 - 振動強度の測定方法 - Google Patents

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】この発明は制御用地震計等に
利用する振動強度の測定方法に関するものである。

【０００２】

【従来の技術】従来、地震の強度に応じて各種のシステ
ムを制御し、被害の拡大や二次災害の発生を防止するた
めの装置として、制御用地震計が利用されている。この
制御用地震計としては、構造物の被害（損傷）との相関
が高いスペクトラム・インテンシティ−（以下、ＳＩ
値）という地震動の強度の尺度を用いて振動の強度を測
定し、これに基づいて判別信号などの制御信号を出力す
るものが提案されている（特開平１０−１２３２５８号
公報など）。

【０００３】地震動が発生すると、その加速度ｙ（ｔ）
の時間的な変化は例えば図９（ａ）のように観測され
る。このときの、ある構造物の速度応答ｖ（ｔ）および
最大速度応答Ｓｖ（ｔ）の時間的変化は図９（ｂ）のよ
うになる。なお、最大速度応答Ｓｖ（ｔ）とは、その時
刻までに発生した速度応答ｖ（ｔ）の最大値のことであ
る。また、このような振動のある時刻において固有周期
Ｔに対する最大速度応答Ｓｖ（ｔ）のスペクトル分布は
図１０のようになっており、この分布の平均値を求める
（すなわち図中の斜線部面積を積分区間で割る）ことに
より、その時刻のＳＩ値を得ることができる。なお、こ
の固有周期Ｔは構造物の大きさなどよって決まるもので
あり、主要な構造物は上記積分区間（０．１から２．
５）の範囲内の固有周期Ｔとなることが判っている。さ
らに詳細な点については、特開昭６２−１２８８４号公
報、特開昭６２−１２８８５号公報、特開昭６２−１２
８８６号公報、特開平６−２１４０４０号公報等の公知
文献を参考にされたい。

【０００４】さて、特開平１０−１２３２５８号公報に
開示されている地震動の測定方法を従来技術として説明
する。図１１に示す波形は、加速度検出部によって検出
された振動の加速度ｙ（ｔ）の時系列デ−タをプロット
して得られる波形である。図中Ｒは測定開始時刻を、Ｓ
は地震発生時刻を、Ｅは地震収束時刻を表す。この時系
列デ−タは極めて短い周期（例えば１０ｍ秒）でサンプ
リングされ、メモリ（ＲＡＭ）に記憶される。この時系
列デ−タからほぼリアルタイムでＳＩ値を得るために、
時刻Ｒから所定時間Ｔｐが経過した時点で、それまでに
得られた加速度ｙ（ｔ）の時系列デ−タをメモリから読
み出し、まず一旦各固有周期Ｔ毎の速度応答値ｖ（ｔ）
（下記式１）に変換し、次いで各固有周期Ｔ毎の最大速
度応答値Ｓｖ（ｔ）（下記式２）を求め、これらからＳ
Ｉ値（下記式３）を算出する。その後、所定の時間Δｔ
ｐが経過した時点で、時刻（Ｒ＋Δｔｐ）から時刻（Ｒ
＋Ｔｐ＋Δｔｐ）までの時系列デ−タをメモリから読み
出して同様の演算を行う。以下、時間Δｔｐが経過する
毎にＳＩ値を算出してゆく。

【０００５】なお、これらの式において、ｔ［ｓ］は現
在時刻、τは積分期間（時間）、Ｔは構造物の固有周
期、ｈは構造物の減衰定数、ωは固有角振動数（ω＝２
π／Ｔ［ｒａｄ／ｓ］）、ωｄは減衰固有角振動数（ω
ｄ＝ω×（１−ｈ² ）^1/2 ［ｒａｄ／ｓ］）である。

【０００６】

【数１】

【数２】

【数３】

【０００７】

【発明が解決しようとする課題】従来の振動強度の測定
方法は以上のような処理を実施するものなので、少なく
とも所定時間Ｔｐにおける加速度ｙ（ｔ）の時系列デ−
タを全てメモリに記憶させておかなければならず、精度
の高い演算結果を得ようとすると膨大なメモリ容量が必
要になってしまうという問題があった。例えば所定時間
Ｔｐを１０秒、サンプリング周期を１０μ秒とすると、
各デ−タ毎に１０００個もの値を記憶しておく必要があ
る。なお、大容量のメモリ装置も市販されてはいるが、
これを使用しても処理速度や価格の点で満足のゆくもの
とはならない。

【０００８】この発明は、演算速度や演算精度を損なう
ことなく、極めて少ないメモリ容量しか必要としない振
動強度の測定方法を得ることを目的とする。

【０００９】

【課題を解決するための手段】この発明に係る振動強度
の測定方法は、加速度検出素子からの加速度検出信号に
基づき速度応答値（ｖ（ｔ））を求め、この速度応答値
から最大速度応答値（Ｓｖ（ｔ））を求め、この最大速
度応答値から地震強度（ＳＩ）を求める振動強度の測定
方法において、所定の時間間隔毎に、速度応答値（ｖ
（ｔ））を記憶し、この記憶した値の所定個数（時間
窓）の中から最大値を求めてそれを最大速度応答値（Ｓ
ｖ（ｔ））とし、この最大速度応答値から地震強度値
（ＳＩ）を求めると共に、これらの演算を上記時間間隔
の整数倍ずつずらして反復的に行うものである。

【００１０】この発明に係る振動強度の測定方法は、加
速度検出素子からの加速度検出信号に基づき速度応答値
（ｖ（ｔ））を求め、この速度応答値から最大速度応答
値（Ｓｖ（ｔ））を求め、この最大速度応答値から地震
強度（ＳＩ）を求める振動強度の測定方法において、所
定の時間間隔毎に、速度応答値（ｖ（ｔ））を採取し、
第一記憶手段（Ｗ０）の記憶値と比較して採取値の方が
大きい場合には第一記憶手段の記憶を採取値に変更し、
第一記憶手段および第二記憶手段の記憶値のうち大きい
方を最大速度応答値（Ｓｖ（ｔ））として、この最大速
度応答値から地震強度値（ＳＩ）を求めるとともに、上
記時間間隔が所定個数に達する毎に、第二記憶手段の記
憶値を第一記憶手段の記憶値に変更した後第一記憶手段
の記憶値を「０」に変更するものである。

【００１１】この発明に係る振動強度の測定方法は、加
速度検出素子からの加速度検出信号に基づき速度応答値
（ｖ（ｔ））を求め、この速度応答値から最大速度応答
値（Ｓｖ（ｔ））を求め、この最大速度応答値から地震
強度（ＳＩ）を求める振動強度の測定方法において、所
定の時間間隔毎に、速度応答値（ｖ（ｔ））を採取し、
第一記憶手段（Ｗ０）の記憶値と比較して採取値の方が
大きい場合には第一記憶手段の記憶を採取値に変更し、
第一の記憶手段および他の複数の記憶手段（Ｗ１，Ｗ
２，Ｗ３）の記憶値のうち最も大きいものを最大速度応
答値（Ｓｖ（ｔ））として、この最大速度応答値から地
震強度値（ＳＩ）を求めるとともに、上記時間間隔が所
定個数に達する毎に、古い記憶値を持つ記憶手段の記憶
値をこれに次いで古い記憶値を持つ記憶手段の記憶値に
順次変更し、第二の記憶手段の記憶値を第一の記憶手段
の記憶値に変更した後当該第一の記憶手段の記憶値を
「０」に変更するものである。

【００１２】

【発明の実施の形態】以下、この発明の実施の一形態を
説明する。 実施の形態１．図１はこの発明のスペクトラム強度値演
算方法に基づく波形記録装置を適用したガス供給システ
ムを示すブロック図である。図において、１は防爆域を
備えたガス製造工場、２はガスタンク、３は高圧導管、
４は防爆域を備えたガバナステ−ション、５は緊急遮断
弁、６は中圧導管、７は防爆域に設置されたガバナ室、
８は低圧導管、９はガスを供給される家庭、１０はガス
を供給される工場である。また、１１は監視室、１２お
よび１３は緊急遮断指示、１４はガバナ、１５は緊急遮
断機能を有したガバナである。

【００１３】１６は３軸（Ｘ，Ｙ，Ｚ軸）計測が可能で
振動に応じて変化する検出信号を出力する半導体加速度
センサ、１７はこの検出信号に基づいて振動判別を行い
判別信号を出力する演算処理回路、１８は演算処理回路
１７から出力されるデジタル（２値）形式の判別信号、
１９は振り子などを利用して振動を検出する機械式地震
センサ、２０は判別信号１８と機械式地震センサ１９の
出力信号とが入力され、この両方共が振動検出である場
合に遮断信号を発生させる判定回路、２１はこの遮断信
号、２２はこの遮断信号２１により防爆域においてガス
流路を遮断するガバナである。なお、半導体加速度セン
サ１６および演算処理回路１７は防爆域内に設置される
ため防爆ケ−スに収容され、機械式地震センサ１９およ
び判定回路２０は非防爆域であるガバナ室外壁面に設置
された制御盤内に配設されている。

【００１４】なお、上記半導体加速度センサ１６として
は、例えば、特開平９−４３０６８号公報に記載された
静電容量の変化を利用したセンサおよびその信号処理回
路を用いることができる。同公報の従来の技術の欄にあ
る通り、固定基板と可撓基板との各対向面に電極を着設
して対向配置される静電容量素子を複数対設け、当該基
板面に平行なＸＹ平面を設定し、これと直交するＺ軸の
Ｘ，Ｙ，Ｚ軸３次元方向の位置の変化を、複数対の静電
容量素子間の静電容量変化に基づき各Ｘ，Ｙ，Ｚ軸方向
成分の検出を行うものである。そして、Ｘ軸方向の加速
度出力として、静電容量素子Ｃ２１とＣ２３との静電容
量差（Ｃ２１−Ｃ２３）、Ｙ軸方向の加速度出力とし
て、静電容量素子Ｃ２２とＣ２４との静電容量差（Ｃ２
２−Ｃ２４）、Ｚ軸方向の加速度出力として、静電容量
素子Ｃ２５の静電容量Ｃ２５あるいはＣ２１＋Ｃ２２＋
Ｃ２３＋Ｃ２４として検出することができる。他にも、
半導体加速度センサ１６としては、静電容量型加速度セ
ンサに限定されず、例えば、ピエゾ抵抗型あるいは圧電
型の加速度センサを使用することもできる。また、半導
体加速度センサ１６は３軸（Ｘ，Ｙ，Ｚ軸）計測ができ
るようにしたほうがよい。これにより、正確な地震動の
振動計測ができる。特に直下型地震では、Ｚ軸（鉛直）
方向が最初に振動することから、より重要な計測情報を
得ることができる。

【００１５】図２はこの発明の実施の形態１による演算
処理回路１７およびその周辺部材の詳細な構成を示すブ
ロック図である。図において、２３は半導体加速度セン
サ１６の検出信号とともに基準電圧回路２４および補正
用温度センサ２５の出力が入力され、これらの出力であ
る検出信号をΔｔ毎にサンプリングして加速度、温度、
基準電圧のデ−タを出力する（Ａ／Ｄ変換する）センサ
インタフェ−ス、２６はメモリ２７、リアルタイムクロ
ック２８、バックアップ回路２９、電源回路３０などが
接続されるとともに上記加速度、温度、基準電圧デ−タ
が入力され、各種振動の特性値を演算したり、デジタル
（２値）形式の判別信号１８を出力する中央処理装置
（ＣＰＵ）などの制御部であり、３１はこの判別信号１
８を外部に出力するデジタル出力インタフェ−スであ
る。

【００１６】制御部２６において、３２はＳＩ値演算部
３２ａによって算出された地震動の強度の一つの尺度で
あるＳＩ値と、震度演算部３２ｂによって算出された震
度と、加速度演算部３２ｃによって算出された加速度
と、変位量演算部３２ｄによって算出された変位量と、
速度演算部３２ｅによって算出された速度の各算出値の
中から任意の複数の算出値を選択する算出値選択手段で
ある。すなわち、この算出値選択手段３２は、各種演算
機能を必要に応じて適宜組み合わせることにより、設置
現場において最適な設定を提供することができる。ま
た、３３はこれらの算出値に基づいて振動判別を行い、
判別信号１８を出力するトリガ信号出力手段、３４は上
記リアルタイムクロック２８のクロック信号に基づいて
判別時間をカウントする第一タイマ、３５は上記リアル
タイムクロック２８のクロック信号に基づいて保持時間
をカウントする第二タイマである。

【００１７】また、３６はデジタル入力インタフェ−
ス、３７は演算処理結果をリレ−接点出力信号として出
力するリレ−インタフェ−ス、３８はアナログ出力イン
タフェ−ス、３９はロ−ダインタフェ−スである。この
リレ−インタフェ−ス３７は、例えば、直流２４ボルト
や交流１００ボルトなどの回路を直接、接続できるよう
に構成されており、したがって、アクチュエ−タやバル
ブ、大型の継電器などの各種の負荷を直接動作させるこ
とができるようになっているものである。また、このリ
レ−インタフェ−ス３７と図示しない接点側回路とは電
気的に分離して設けられている。なお、リレ−インタフ
ェ−ス３７は、リレ−だけでなく、フォトカプラやフォ
トモスリレ−などの分離絶縁型電気部品にも接続するこ
とができる。ロ−ダインタフェ−ス３９は、携帯型設定
器を接続するための図示しないロ−ダ接続端子を備え、
例えば特許２５２３０５３号公報に開示された手段によ
り算出値選択手段３２に算出値を選択させるための設定
信号などを与えるものである。

【００１８】なお、電源回路３０への外部電源からの配
線や各インタフェ−ス３１，３６，３７，３８，３９の
配線は、防爆ケ−スの外部に引き出されている。また、
判定信号１８はデジタル信号となっているのでノイズに
強くなっており、この配線をガバナ室７から制御盤まで
引き回したとしても誤動作し難くなっている。

【００１９】次に動作について説明する。センサインタ
フェ−スが例えば１０ｍｓ毎などの所定のサンプリング
周期毎に半導体加速度センサ１６の検出信号出力などを
サンプリングして加速度、温度、基準電圧のデ−タを出
力すると、制御部２６はこれらのデ−タに基づいて各種
振動の特性値を演算したり、デジタル（２値）形式の判
別信号１８を出力したりする。そして、例えばこの判別
信号１８が出力されると、判定回路２０は機械式地震セ
ンサ１９の出力も振動検出である場合には遮断信号２１
を発生させ、ガバナ２２はこの遮断信号２１に基づいて
防爆域においてガス流路を自動的に遮断する。これによ
り、地震で低圧導管８などに亀裂などが発生したとして
も、不要なガス漏れを防止することができる。

【００２０】なお、これとともに監視室１１から緊急遮
断指示１２，１３が出力されれば、ガバナ１４や緊急遮
断機能を有したガバナ１５においてもガス流路を遮断す
ることとなり、より確実に不要なガス漏れを防止するこ
とができる。

【００２１】次に、地震強度としてＳＩ値を測定する例
を、図３及び図４を参照しつつ以下に説明する。図３は
発明の実施の形態１のアルゴリズムを説明する模式図で
あり、図４は発明の実施の形態１の手順を示すフロ−チ
ャ−トである。なお、この図４のフロ−チャ−トは上記
半導体加速度センサ１６からの加速度が得られるサンプ
リングレ−トと同一周期毎に実施されるものである。但
し、サンプリングレ−トとこのフロ−チャ−トの演算周
期Δｔとは異なるものとすることも可能である。

【００２２】上記半導体加速度センサ１６からの出力は
時間的に変化する加速度信号であり、センサインタフェ
−ス２３によりデジタルデ−タに変換されて制御部２６
に取り込まれる（図４のＳＴ１１）。この加速度デ−タ
ｙ（ｔ）を受けたＳＩ値演算部３２ａは、極めて短い演
算間隔Δｔ（例えば１０ｍ秒）毎に速度応答値ｖ（ｔ）
に変換演算を行う（図４のＳＴ１２）。この変換演算は
上記式１に示すものであるが、実際には下記式群４に示
す直接積分法などを用いてデジタル的に処理する。これ
により毎回１ｍ秒以下で加速度値ｙ（ｔ）を速度応答値
ｖ（ｔ）に変換でき、図９（ｂ）に示すような速度応答
値ｖ（ｔ）の時系列デ−タを、実質的にリアルタイムに
得ることができる。

【００２３】

【数４】

【００２４】なお、固有周期Ｔを与える定数（数式１の
定数ω）はあらかじめ複数（例えば７個）定められてお
り、各定数毎に（すなわち各固有周期Ｔ毎に）速度応答
値ｖ（ｔ）の時系列デ−タが算出される。また、上記加
速度センサ１６は直交２軸方向（例えば東西方向及び南
北方向）の加速度信号を出力しており、演算により３６
０度任意の方向の加速度信号に変換することができるの
で、必要に応じて任意方向に関する速度応答値ｖ（ｔ）
の時系列デ−タを上記同様リアルタイムに演算すること
もできる。図３は、言うまでもないが、ある一つの固有
周期及び一つの方向に関して表現したものである。

【００２５】加速度応答ｖ（ｔ）から最大速度応答Ｓｖ
（ｔ）及びＳＩ値を求めるアルゴリズムを説明する。図
３（ａ）は速度応答ｖ（ｔ）の時間的変化を、図３
（ｂ）は次に述べる時間窓の時間的変化を、図３（ｃ）
は最大速度応答Ｓｖ（ｔ）の時間的変化を、それぞれ表
している。演算間隔Δｔのｎ倍（ｎは整数）だけ過去か
ら現時刻までの所定時間を時間窓と呼ぶことにする（図
３では簡単のためにｎ＝３の場合を例に図示している
が、実際にはｎ＝１０００程度が適用される）。速度応
答値ｖ（ｔ）用のメモリ領域には、この時間窓内での値
のみが記憶されるようになっている。

【００２６】ここで現時刻を仮にｔ５とすると、時間窓
はｔ２〜ｔ５となり、メモリ領域には４個の速度応答値
ｖ（ｔ２），ｖ（ｔ３），ｖ（ｔ４），ｖ（ｔ５）が記
憶されている（図４のＳＴ１３）。そして、この時刻が
ｎ（ここでは１）の倍数である場合には以下の処理を実
施する（図４のＳＴ１４）。これらの中で絶対値が最も
大きいものを、時刻ｔ５における最大速度応答値Ｓｖ
（ｔ５）として別のメモリ領域に格納する（図４のＳＴ
１５）。同時に、他の固有周期Ｔに関しても同様に最大
速度応答値Ｓｖ（ｔ）を選出する。そして全部（例では
７個）の最大速度応答値Ｓｖ（ｔ）から時刻ｔ５におけ
るＳＩ値（例では7個の平均値）を得る（図４のＳＴ１
６）。ちなみに時刻ｔ５からＳＩ値を得るまでの演算時
間は３ｍ秒程度である。

【００２７】次に、時刻ｔ５から演算間隔Δｔが経過し
て時刻ｔ６になると、時間窓はｔ３〜ｔ６となり、メモ
リ領域から速度応答値ｖ（ｔ２）が消去され、その代わ
りに新たな速度応答値ｖ（ｔ６）が格納される（図４の
ＳＴ１１〜ＳＴ１３）。すなわち時刻ｔ６においてメモ
リ領域に記憶されている速度応答値はｖ（ｔ３）、ｖ
（ｔ４）、ｖ（ｔ５）、ｖ（ｔ６）の４個となる。以
下、上記同様に時刻ｔ６における最大速度応答値Ｓｖ
（ｔ）とＳＩ値が算出される（図４のＳＴ１５、ＳＴ１
６）。

【００２８】前述の従来技術においては時間窓の範囲内
全ての加速度ｙ（ｔ）をメモリ２７に記憶しておいて後
で全ての演算を行うのに対し、この発明では加速度ｙ
（ｔ）から速度応答ｖ（ｔ）にリアルタイムで変換しそ
れを記憶するようにしたので、ＳＩ値の演算をするごと
に加速度ｙ（ｔ）から速度応答ｖ（ｔ）を算出する必要
がなく、速度応答ｖ（ｔ）の最大値Ｓｖ（ｔ）を求める
比較手順だけでよくなるため演算量が飛躍的に減少し、
制御部２６に用いるＣＰＵの演算負荷を大幅に減少させ
ることができた。

【００２９】以上のように、この実施の形態１によれ
ば、半導体加速度センサ１６からの加速度検出信号ｙ
（ｔ）に基づき速度応答値（ｖ（ｔ））を求め、この速
度応答値から最大速度応答値（Ｓｖ（ｔ））を求め、こ
の最大速度応答値から地震強度（ＳＩ）を求める振動強
度の測定方法において、半導体加速度センサ１６のサン
プリング間隔毎に、加速度検出信号ｙ（ｔ）に基づいて
演算された速度応答値（ｖ（ｔ））を記憶し、この記憶
した値の所定個数（時間窓）の中から最大値を求めてそ
れを最大速度応答値（Ｓｖ（ｔ））とし、この最大速度
応答値から地震強度値（ＳＩ）を求めると共に、これら
の演算を上記サンプリング間隔の整数倍（ここでは１
倍）ずつずらして反復的に行うので、従来のように地震
強度値（ＳＩ）を求める際に加速度検出信号ｙ（ｔ）に
基づき速度応答値（ｖ（ｔ））を求める必要がなくな
る。

【００３０】従って、地震強度値（ＳＩ）を求める度に
同じ加速度検出信号ｙ（ｔ）に基づいて速度応答値（ｖ
（ｔ））を求めることがなくなり（つまり同一の演算を
重複して実施することがなくなり）、その分ト−タルの
演算処理負荷を削減することができる。また、従来のも
のと比べた場合、地震強度値（ＳＩ）を求めるタイミン
グにおいてはそのタイミングにおける1つの速度応答値
演算を実施すればいいので、当該タイミングにおける瞬
時的な演算処理負荷を格段に削減することもできる。そ
れ故、演算速度や演算精度を損なうことなく演算負荷を
大幅に減少させることができる効果がある。

【００３１】実施の形態２．この実施の形態２は、上述
の実施の形態１に対して、さらにメモリ２７の使用量を
減少させたものであり、図５及び図６を用いて説明す
る。図５はこの発明の実施の形態２のアルゴリズムを説
明する模式図であり、図６はこの発明の実施の形態２の
手順を示すフロ−チャ−トである。なお、上記図６のフ
ロ−チャ−トは上記半導体加速度センサ１６からの加速
度が得られるサンプリングレ−トと同一周期毎に実施さ
れるものである。但し、サンプリングレ−トとこのフロ
−チャ−トの演算周期Δｔとは異なるものとすることも
可能である。また、速度応答ｖ（ｔ）をリアルタイムで
算出するところまでは実施の形態１と同様なので、ここ
では詳細説明を省略する（図６のＳＴ２４、ＳＴ２
５）。

【００３２】加速度応答ｖ（ｔ）から最大速度応答Ｓｖ
（ｔ）及びＳＩ値を求めるアルゴリズムを説明する。図
５（ａ）は速度応答ｖ（ｔ）の時間的変化を、図５
（ｂ）は時間窓の時間的変化を、図５（ｄ）は最大速度
応答Ｓｖ（ｔ）の時間的変化をそれぞれ表している。ま
た、図５（ｃ）は第一メモリ領域Ｗ０と第二メモリ領域
Ｗ２の使用状況を表している。

【００３３】この実施の形態２において特徴的な点は、
時間窓の長さが伸縮する点である。時刻ｔ１における時
間窓の長さはｔ０〜ｔ１、時刻ｔ２における時間窓はｔ
０〜ｔ２、というように演算間隔Δｔの経過と共に長さ
がΔｔ×ｍ（ｍ＝２ｎ−１、ｎは定数）に達するまで時
間窓が伸びてゆく（図５の例ではｎ＝４）。さらに時間
間隔Δｔが経過したとき（すなわち、時刻がｎの整数倍
に達したとき）、時間窓の長さは一旦Δｔ×ｎにまで縮
小される。そして再び演算間隔Δｔの経過と共に、時間
窓の長さはΔｔ×ｍの長さに達するまで伸びる。以後、
時間窓の長さはΔｔ×ｎとΔｔ×ｍとの間で伸縮を繰り
返す。

【００３４】さて、時刻ｔ０においては第一メモリ領域
Ｗ０及び第二メモリ領域Ｗ１には初期値０が格納されて
いる。演算時間間隔Δｔが経過して時刻ｔ１になると時
間窓はｔ０〜ｔ１となり（図６のＳＴ２１、ＳＴ２４、
ＳＴ２５）、第一メモリ領域Ｗ０の記憶値（両方とも
０）と速度応答値ｖ（ｔ１）の絶対値とが比較され、速
度応答値ｖ（ｔ１）の絶対値がＷ０よりも大きい場合に
は、第一のメモリ領域Ｗ０に速度応答値ｖ（ｔ１）の絶
対値がＷ０の最大速度応答値Ｓｖ（Ｗ０）として０に代
わって格納される（図６のＳＴ２６）。そしてＷ０とＷ
１との大きい方の値をその時刻における最大速度応答Ｓ
ｖ（ｔ１）とする（図６のＳＴ２７）。そして全部（例
では７個）の最大速度応答値Ｓｖ（ｔ）から時刻ｔ５に
おけるＳＩ値を得る（図６のＳＴ２８）。

【００３５】さらに演算時間間隔Δｔが経過して時刻ｔ
２になると時間窓はｔ０〜ｔ２となり（図６のＳＴ２
１、ＳＴ２４、ＳＴ２５）、第一メモリ領域Ｗ０の記憶
値（Ｓｖ（Ｗ０））および第二メモリ領域Ｗ１の記憶値
（０）と速度応答値ｖ（ｔ２）の絶対値とが比較され、
速度応答値ｖ（ｔ２）の絶対値が両者のいずれよりも大
きい場合には、第一のメモリ領域Ｗ０を速度応答値ｖ
（ｔ２）の絶対値で書き換え（図６のＳＴ２６）、最大
速度応答及びＳＩ値を求める（図６のＳＴ２７、ＳＴ２
８）。

【００３６】以下、上記同様に時刻ｔ３において手順を
繰り返して、時刻ｔ４に達すると（図６のＳＴ２１）、
第二メモリ領域Ｗ０の値を第一メモリ領域Ｗ１の最大速
度応答値Ｓｖ（ｔ４）に書き換える（図６のＳＴ２２）
と共に第一メモリ領域Ｗ０の値を０に書き換える（図６
のＳＴ２３）。また、この後上記同様に図６のＳＴ２４
からＳＴ２８の処理を実行してＳＩ値を算出する。

【００３７】さらに演算時間間隔Δｔが経過して時刻ｔ
５になると、時間窓はｔ０〜ｔ５となり（図６のＳＴ２
１、ＳＴ２４、ＳＴ２５）、第一メモリ領域Ｗ０の記憶
値（０）および第二メモリ領域Ｗ１の記憶値（Ｓｖ（ｔ
４））と速度応答値ｖ（ｔ５）の絶対値とが比較され、
速度応答値ｖ（ｔ５）の絶対値が両者のいずれよりも大
きい場合には、第一のメモリ領域Ｗ０に速度応答値ｖ
（ｔ５）の絶対値が０に代わって格納される（図６のＳ
Ｔ２６）。そして再び最大速度応答及びＳＩ値を求める
（図６のＳＴ２７、ＳＴ２８）。

【００３８】これらの手順が繰り返され、やがて時刻ｔ
８に達すると、第二メモリ領域Ｗ１の値を第一メモリ領
域Ｗ０の最大速度応答値Ｓｖ（ｔ８）に書き換えると共
に第一メモリ領域Ｗ０の値を０に書き換える（図６のＳ
Ｔ２２、ＳＴ２３）。このとき、今まで第二メモリ領域
Ｗ１に格納されていたＳｖ（ｔ４）が消失する（すなわ
ち時刻ｔ０〜ｔ４に至るまでの時間窓が消失する）。

【００３９】前述の発明の実施の形態１においては時間
窓分のデ−タを記憶するためにｎ個（例では５個）のメ
モリ領域が必要とされるが、この実施の形態２において
は二個のメモリ領域でこと足りるので、さらに少ないメ
モリ２７で演算処理が可能となる。

【００４０】以上のように、この実施の形態２によれ
ば、半導体加速度センサ１６からの加速度検出信号ｙ
（ｔ）に基づき速度応答値（ｖ（ｔ））を求め、この速
度応答値から最大速度応答値（Ｓｖ（ｔ））を求め、こ
の最大速度応答値から地震強度（ＳＩ）を求める振動強
度の測定方法において、半導体加速度センサ１６のサン
プリング間隔毎に、加速度検出信号ｙ（ｔ）に基づく速
度応答値（ｖ（ｔ））を採取し、第一記憶手段（Ｗ０）
の記憶値と比較して採取値の方が大きい場合には第一記
憶手段の記憶を採取値に変更し、第一の記憶手段および
他の複数の記憶手段（Ｗ１，Ｗ２，Ｗ３）の記憶値のう
ち最も大きいものを最大速度応答値（Ｓｖ（ｔ））とし
て、この最大速度応答値から地震強度値（ＳＩ）を求め
るとともに、上記サンプリング間隔が所定個数（ここで
は４）に達する毎に、古い記憶値を持つ記憶手段Ｗ１の
記憶値をこれに次いで古い記憶値を持つ記憶手段Ｗ０の
記憶値に順次変更し、第二の記憶手段Ｗ１の記憶値を第
一の記憶手段Ｗ０の記憶値に変更した後当該第一の記憶
手段Ｗ０の記憶値を「０」に変更するので、従来のよう
に地震強度値（ＳＩ）を求める際に加速度検出信号ｙ
（ｔ）に基づき全ての速度応答値（ｖ（ｔ））を求める
必要がなく、しかも、全ての速度応答値（ｖ（ｔ））に
ついて速度応答ｖ（ｔ）の最大値を求めるための比較手
順を実施する必要が無くなる。

【００４１】従って、地震強度値（ＳＩ）を求める度に
同じ加速度検出信号ｙ（ｔ）に基づいて速度応答値（ｖ
（ｔ））を求めることがなくなり（つまり同一の演算を
重複して実施することがなくなり）、且つ、全ての速度
応答値（ｖ（ｔ））について速度応答ｖ（ｔ）の最大値
を求める比較手順を実施する必要が無くなるので、その
分ト−タルの演算処理負荷を削減することができる。ま
た、実施の形態１と比べても、地震強度値（ＳＩ）を求
めるタイミングにおいては速度応答値同士の比較判断の
回数を半分にすることができるので、当該タイミングに
おける瞬時的な演算処理負荷を格段に削減することもで
きる。それ故、演算速度や演算精度を損なうことなく演
算負荷を大幅に減少させることができる効果がある。

【００４２】また、リアルタイムで比較判断を実施する
ことで記憶させる速度応答値ｖ（ｔ）を半分に削減して
いるので、演算速度や演算精度を損なうことなく極めて
少ないメモリ容量で振動強度を測定することができる効
果がある。

【００４３】この実施の形態２によれば、半導体加速度
センサ１６からの加速度検出信号ｙ（ｔ）に基づき速度
応答値（ｖ（ｔ））を求め、この速度応答値から最大速
度応答値（Ｓｖ（ｔ））を求め、この最大速度応答値か
ら地震強度（ＳＩ）を求める振動強度の測定方法におい
て、半導体加速度センサ１６のサンプリング間隔毎に、
加速度検出信号ｙ（ｔ）に基づく速度応答値（ｖ
（ｔ））を採取し、第一記憶手段（Ｗ０）の記憶値と比
較して採取値の方が大きい場合には第一記憶手段Ｗ０の
記憶を採取値に変更し、第一記憶手段Ｗ０および第二記
憶手段Ｗ１の記憶値のうち大きい方を最大速度応答値
（Ｓｖ（ｔ））として、この最大速度応答値から地震強
度値（ＳＩ）を求めるとともに、上記サンプリング間隔
が所定個数（ここでは４）に達する毎に、第二記憶手段
Ｗ１の記憶値を第一記憶手段Ｗ０の記憶値に変更した後
第一記憶手段Ｗ０の記憶値を「０」に変更するので、直
前の時間窓（所定数の時間間隔を合わせてなる時間帯の
こと）のデ−タから得た最大速度応答値と、現在の時間
窓の最大速度応答値との二個の値を記憶すればこと足り
るので、上述した場合に比べても極めて少ない最小のメ
モリ容量で、且つ、演算速度や演算精度を損なうことな
く振動強度を測定することができる効果がある。

【００４４】実施の形態３．発明の実施の形態３は、前
記発明の実施の形態１と実施の形態２との中間的な形態
の発明であり、図７および図８を参照して説明する。図
７は発明の実施の形態３のアルゴリズムを説明する模式
図であり、図８は発明の実施の形態３の手順を示すフロ
−チャ−トである。図７（ａ）は速度応答ｖ（ｔ）の時
間的変化を、図７（ｂ）は次に述べる時間窓の時間的変
化を、図７（ｃ）は時間窓の長さとそれに対応するメモ
リ２７の時間的変化を、図７（ｄ）は最大速度応答Ｓｖ
（ｔ）の時間的変化を、それぞれ表わしている。なお、
上記図８のフロ−チャ−トは上記半導体加速度センサ１
６からの加速度が得られるサンプリングレ−トと同一周
期毎に実施されるものである。但し、サンプリングレ−
トとこのフロ−チャ−トの演算周期Δｔとは異なるもの
とすることも可能である。この実施の形態３は、上記実
施の形態２においてメモリ領域が４つに増えたものに相
当し、その動作も実施の形態２に準ずるものなので、ｔ
０からｔ８まての動作説明を省略して以下簡略に説明す
る。

【００４５】現在時刻がｔ８が既に経過して時刻ｔ９に
至っていない時点であると仮定し、時間窓の長さはｔ２
〜ｔ８であり、第一のメモリ領域Ｗ０には速度応答ｖ
（ｔ８）が、第二のメモリ領域Ｗ１には速度応答ｖ（ｔ
６）とｖ（ｔ７）とのうちの大きいほうが、第三のメモ
リ領域Ｗ２には速度応答ｖ（ｔ４）とｖ（ｔ５）とのう
ちの大きいほうが、第四のメモリ領域Ｗ３には速度応答
ｖ（ｔ２）とｖ（ｔ３）とのうちの大きいほうが、それ
ぞれ記憶されているものとする。

【００４６】時刻ｔ９になると（図８のＳＴ３１）まず
加速度を読み込み（図８のＳＴ３６）、速度応答ｖ（ｔ
９）を算出し（図８のＳＴ３７）、この値と第一のメモ
リ領域Ｗ０の記憶値ｖ（ｔ８）とを比較してどちらか大
きいほうを第一のメモリ領域Ｗ０に上書きし（図８のＳ
Ｔ３８）、第一から第四のメモリ領域Ｗ０〜Ｗ３のそれ
ぞれの記憶値の中から最も大きなものを選択してこれを
最大速度応答Ｓｖ（ｔ９）として別の時系列メモリ領域
に記憶し（図８のＳＴ３９）、最大速度応答Ｓｖ（ｔ
９）値からＳＩ値を算出する（図８のＳＴ４０）。

【００４７】時刻がｔ１０になると（図８のＳＴ３
１）、第四のメモリ領域Ｗ３に第三のメモリ領域Ｗ２の
記憶値を上書きし（図８のＳＴ３２）、第三のメモリ領
域Ｗ２に第二のメモリ領域Ｗ１の記憶値を上書きし（図
８のＳＴ３３）、第二のメモリ領域Ｗ１に第一のメモリ
領域Ｗ０の記憶値を上書きし（図８のＳＴ３４）、第一
のメモリ領域Ｗ０に「０」を上書きする（図８のＳＴ３
５）ので、時刻ｔ２およびｔ３において得た情報は消失
する。次いで上記同様に図８のＳＴ３６からＳＴ４０の
処理を実行してＳＩ値を算出する。

【００４８】この実施の形態では、メモリ領域の数を４
つとしたが、もちろんこれに限るものではなく必要に応
じて増減することができる。また、一つのメモリ領域に
対応する演算間隔Δｔの数を２（すなわちｎ＝２）とし
て説明したが、実際にはｎ＝１０００程度にとることも
可能である。

【００４９】前述の発明の実施の形態１においては時間
窓分のデ−タを記憶するためにｎ個（例えば５個）のメ
モリ領域が必要とされるが、この実施の形態３において
は任意の複数個（例えば４個）のメモリ領域でこと足り
るので、さらに少ないメモリ２７で演算処理が可能とな
る。

【００５０】以上のように、この実施の形態３によれ
ば、半導体加速度センサ１６からの加速度検出信号ｙ
（ｔ）に基づき速度応答値（ｖ（ｔ））を求め、この速
度応答値から最大速度応答値（Ｓｖ（ｔ））を求め、こ
の最大速度応答値から地震強度（ＳＩ）を求める振動強
度の測定方法において、半導体加速度センサ１６のサン
プリング間隔毎に、加速度検出信号ｙ（ｔ）に基づく速
度応答値（ｖ（ｔ））を採取し、第一記憶手段Ｗ０の記
憶値と比較して採取値の方が大きい場合には第一記憶手
段の記憶を採取値に変更し、第一の記憶手段および他の
複数の記憶手段（Ｗ１，Ｗ２，Ｗ３）の記憶値のうち最
も大きいものを最大速度応答値（Ｓｖ（ｔ））として、
この最大速度応答値から地震強度値（ＳＩ）を求めると
ともに、上記サンプリング間隔が所定個数（ここでは
２）に達する毎に、古い記憶値を持つ記憶手段の記憶値
をこれに次いで古い記憶値を持つ記憶手段の記憶値に順
次変更し、第二の記憶手段Ｗ１の記憶値を第一の記憶手
段Ｗ０の記憶値に変更した後当該第一の記憶手段Ｗ０の
記憶値を「０」に変更するので、従来のように地震強度
値（ＳＩ）を求める際に加速度検出信号ｙ（ｔ）に基づ
き全ての速度応答値（ｖ（ｔ））を求める必要がなく、
しかも、全ての速度応答値（ｖ（ｔ））について速度応
答ｖ（ｔ）の最大値を求めるための比較手順を実施する
必要が無くなる。

【００５１】従って、地震強度値（ＳＩ）を求める度に
同じ加速度検出信号ｙ（ｔ）に基づいて速度応答値（ｖ
（ｔ））を求めることがなくなり（つまり同一の演算を
重複して実施することがなくなり）、且つ、全ての速度
応答値（ｖ（ｔ））について速度応答ｖ（ｔ）の最大値
を求める比較手順を実施する必要が無くなるので、その
分ト−タルの演算処理負荷を削減することができる。ま
た、実施の形態１と比べても、地震強度値（ＳＩ）を求
めるタイミングにおいては速度応答値同士の比較判断の
回数を少なくとも半分以下にすることができるので、当
該タイミングにおける瞬時的な演算処理負荷を格段に削
減することもできる。それ故、演算速度や演算精度を損
なうことなく演算負荷を大幅に減少させることができる
効果がある。

【００５２】また、リアルタイムで比較判断を実施する
ことで記憶させる速度応答値ｖ（ｔ）を削減しているの
で、演算速度や演算精度を損なうことなく極めて少ない
メモリ容量で振動強度を測定することができる効果があ
る。

【００５３】

【発明の効果】以上のように、請求項１に係る発明によ
れば、加速度検出素子からの加速度検出信号に基づき速
度応答値（ｖ（ｔ））を求め、この速度応答値から最大
速度応答値（Ｓｖ（ｔ））を求め、この最大速度応答値
から地震強度（ＳＩ）を求める振動強度の測定方法にお
いて、所定の時間間隔毎に、速度応答値（ｖ（ｔ））を
記憶し、この記憶した値の所定個数（時間窓）の中から
最大値を求めてそれを最大速度応答値（Ｓｖ（ｔ））と
し、この最大速度応答値から地震強度値（ＳＩ）を求め
ると共に、これらの演算を上記時間間隔の整数倍ずつず
らして反復的に行うので、従来のように地震強度値（Ｓ
Ｉ）を求める際に加速度検出信号に基づき速度応答値
（ｖ（ｔ））を求める必要がなくなる。

【００５４】従って、地震強度値（ＳＩ）を求める度に
同じ加速度検出信号に基づいて速度応答値（ｖ（ｔ））
を求めることがなくなり（つまり同一の演算を重複して
実施することがなくなり）、その分ト−タルの演算処理
負荷を削減することができる。また、従来のものと比べ
た場合、地震強度値（ＳＩ）を求めるタイミングにおい
てはそのタイミングにおける1つの速度応答値演算を実
施すればいいので、当該タイミングにおける瞬時的な演
算処理負荷を格段に削減することもできる。それ故、演
算速度や演算精度を損なうことなく演算負荷を大幅に減
少させることができる効果がある。

【００５５】請求項２および請求項３に係る発明によれ
ば、加速度検出素子からの加速度検出信号に基づき速度
応答値（ｖ（ｔ））を求め、この速度応答値から最大速
度応答値（Ｓｖ（ｔ））を求め、この最大速度応答値か
ら地震強度（ＳＩ）を求める振動強度の測定方法におい
て、所定の時間間隔毎に、速度応答値（ｖ（ｔ））を採
取し、第一記憶手段（Ｗ０）の記憶値と比較して採取値
の方が大きい場合には第一記憶手段の記憶を採取値に変
更し、第一の記憶手段および他の複数の記憶手段（Ｗ
１，Ｗ２，Ｗ３）の記憶値のうち最も大きいものを最大
速度応答値（Ｓｖ（ｔ））として、この最大速度応答値
から地震強度値（ＳＩ）を求めるとともに、上記時間間
隔が所定個数に達する毎に、古い記憶値を持つ記憶手段
の記憶値をこれに次いで古い記憶値を持つ記憶手段の記
憶値に順次変更し、第二の記憶手段の記憶値を第一の記
憶手段の記憶値に変更した後当該第一の記憶手段の記憶
値を「０」に変更するので、従来のように地震強度値
（ＳＩ）を求める際に加速度検出信号に基づき速度応答
値（ｖ（ｔ））を求める必要がなく、しかも、全ての速
度応答値（ｖ（ｔ））について速度応答ｖ（ｔ）の最大
値を求める比較手順を実施する必要が無くなる。

【００５６】従って、地震強度値（ＳＩ）を求める度に
同じ加速度検出信号に基づいて速度応答値（ｖ（ｔ））
を求めることがなくなり（つまり同一の演算を重複して
実施することがなくなり）、且つ、全ての速度応答値
（ｖ（ｔ））について速度応答ｖ（ｔ）の最大値を求め
る比較手順を実施する必要が無くなるので、その分ト−
タルの演算処理負荷を削減することができる。また、上
述したものと比べても、地震強度値（ＳＩ）を求めるタ
イミングにおいては速度応答値同士の比較判断の回数を
少なくとも半分以下にすることができるので、当該タイ
ミングにおける瞬時的な演算処理負荷を格段に削減する
こともできる。それ故、演算速度や演算精度を損なうこ
となく演算負荷を大幅に減少させることができる効果が
ある。

【００５７】また、リアルタイムで比較判断をして記憶
させる速度応答値を削減しているので、演算速度や演算
精度を損なうことなく極めて少ないメモリ容量で振動強
度を測定することができる効果がある。

【図面の簡単な説明】

【図１】この発明を適用したガス供給システムを示すブ
ロック図である。

【図２】この発明の実施の形態１による演算処理回路お
よびその周辺の詳細な構成を示すブロック図である。

【図３】この発明の実施の形態１のアルゴリズムを説明
する模式図である。

【図４】この発明の実施の形態１の手順を示すフロ−チ
ャ−トである。

【図５】この発明の実施の形態２のアルゴリズムを説明
する模式図である。

【図６】この発明の実施の形態２の手順を示すフロ−チ
ャ−トである。

【図７】この発明の実施の形態３のアルゴリズムを説明
する模式図である。

【図８】この発明の実施の形態３の手順を示すフロ−チ
ャ−トである。

【図９】ある固有周期における加速度、速度応答及び最
大速度応答の関係を示す波形図である。

【図１０】ある時点における最大速度応答のスペクトル
を示す波形図である。

【図１１】従来技術を示す波形図である。

【符号の説明】

１ ガス製造工場 ２ ガスタンク ３ 高圧導管 ４ ガバナステ−ション ５ 緊急遮断弁 ６ 中圧導管 ７ ガバナ室 ８ 低圧導管 ９ 家庭 １０ 工場 １１ 監視室 １２，１３ 緊急遮断指示 １４ ガバナ １５ 緊急遮断機能を有したガバナ １６ 半導体加速度センサ １７ 演算処理回路 １８ 判別信号 １９ 機械式地震センサ ２０ 判定回路 ２１ 遮断信号 ２２ ガバナ ２３ センサインタフェ−ス ２４ 基準電圧回路 ２５ 補正用温度センサ ２６ 制御部 ２７ メモリ ２８ リアルタイムクロック ２９ バックアップ回路 ３０ 電源回路 ３１ デジタル出力インタフェ−ス ３２ 算出値選択手段 ３２ａ ＳＩ値演算部 ３２ｂ 震度演算部 ３２ｃ 加速度演算部 ３２ｄ 変位量演算部 ３２ｅ 速度演算部 ３３ トリガ信号出力手段 ３４ 第一タイマ ３５ 第二タイマ ３６ デジタル入力インタフェ−ス ３７ リレ−インタフェ−ス ３８ アナログ出力インタフェ−ス ３９ ロ−ダインタフェ−ス

───────────────────────────────────────────────────── フロントページの続き (72)発明者 田久保 光 東京都渋谷区渋谷２丁目12番19号 株式 会社山武内 (72)発明者 清水 善久 埼玉県幸手市香日向２−26−３ (72)発明者 小金丸 健一 東京都豊島区東池袋１−48−６−908 (56)参考文献 特開 平10−123258（ＪＰ，Ａ) 特開 昭62−12884（ＪＰ，Ａ) 特開 平８−36062（ＪＰ，Ａ) 特開 平３−262930（ＪＰ，Ａ) 特開 平10−123257（ＪＰ，Ａ) 特開 昭61−234382（ＪＰ，Ａ) 特開 平３−262929（ＪＰ，Ａ) 特開2000−249771（ＪＰ，Ａ) 特開2000−249772（ＪＰ，Ａ) 特公 平４−35035（ＪＰ，Ｂ２) 特公 平６−35935（ＪＰ，Ｂ２) 特許2545504（ＪＰ，Ｂ２) 特許2800911（ＪＰ，Ｂ２) 特許2957572（ＪＰ，Ｂ１) 小金丸健一、清水善久、簗田貴、古川 洋之、田久保光，“新ＳＩセンサーのＳ Ｉ値計算方法”，第33回地盤工学研究発 表会−平成10年度発表講演集（２分冊の ２）−，日本，社団法人地盤工学会, 1998年 ５月25日，1233，ｐ．2471− 2472 Ｉｋｕｏ Ｔｏｗｈａｔａ， Ｊ． Ｋ． Ｐａｒｋ， Ｒ． Ｐ． Ｏｒｅ ｎｓｅ ａｎｄ Ｈｉｒｏｙｕｋｉ Ｋ ａｎｏ， ”Ｕｓｅ ｏｆ ｓｐｅｃｔ ｒｕｍ ｉｎｔｅｎｓｉｔｙ ｆｏｒ ｉｍｍｅｄｉａｔｅ ｄｅｔｅｃｔｉｏ ｎ ｏｆ ｓｕｂｓｏｉｌ ｌｉｑｕｅ ｆａｃｔｉｏｎ”，Ｓｏｉｌｓ ａｎｄ Ｆｏｕｎｄａｔｉｏｎｓ，日本，Ｊａ ｐａｎｅｓｅ Ｇｅｏｔｅｃｈｎｉｃａ ｌ Ｓｏｃｉｅｔｙ，1996年 ６月，第 36巻、第２号，ｐ．29−44 (58)調査した分野(Int.Cl.⁷，ＤＢ名) G01V 1/00 - 1/52 G01H 1/00 - 1/16 G01H 11/00 - 11/08 ＪＩＣＳＴファイル（ＪＯＩＳ)

Claims (3)

(57)【特許請求の範囲】
1. 【請求項１】 加速度検出素子からの加速度検出信号に
   基づき速度応答値を求め、この速度応答値から最大速度
   応答値を求め、この最大速度応答値から地震強度を求め
   る振動強度の測定方法において、 所定の時間間隔毎に、速度応答値を記憶し、この記憶し
   た値の所定個数の中から最大値を求めてそれを最大速度
   応答値とし、この最大速度応答値から地震強度値を求め
   ると共に、これらの演算を上記時間間隔の整数倍ずつず
   らして反復的に行うことを特徴とする振動強度の測定方
   法。
2. 【請求項２】 加速度検出素子からの加速度検出信号に
   基づき速度応答値を求め、この速度応答値から最大速度
   応答値を求め、この最大速度応答値から地震強度を求め
   る振動強度の測定方法において、 所定の時間間隔毎に、速度応答値を採取し、第一記憶手
   段の記憶値と比較して採取値の方が大きい場合には第一
   記憶手段の記憶を採取値に変更し、第一記憶手段および
   第二記憶手段の記憶値のうち大きい方を最大速度応答値
   として、この最大速度応答値から地震強度値を求めると
   ともに、 上記時間間隔が所定個数に達する毎に、第二記憶手段の
   記憶値を第一記憶手段の記憶値に変更した後第一記憶手
   段の記憶値を「０」に変更することを特徴とする振動強
   度の測定方法。
3. 【請求項３】 加速度検出素子からの加速度検出信号に
   基づき速度応答値を求め、この速度応答値から最大速度
   応答値を求め、この最大速度応答値から地震強度を求め
   る振動強度の測定方法において、 所定の時間間隔毎に、速度応答値を採取し、第一記憶手
   段の記憶値と比較して採取値の方が大きい場合には第一
   記憶手段の記憶を採取値に変更し、第一の記憶手段およ
   び他の複数の記憶手段の記憶値のうち最も大きいものを
   最大速度応答値として、この最大速度応答値から地震強
   度値を求めるとともに、 上記時間間隔が所定個数に達する毎に、古い記憶値を持
   つ記憶手段の記憶値をこれに次いで古い記憶値を持つ記
   憶手段の記憶値に順次変更し、第二の記憶手段の記憶値
   を第一の記憶手段の記憶値に変更した後当該第一の記憶
   手段の記憶値を「０」に変更することを特徴とする振動
   強度の測定方法。