計測震度の算出方法

　計測震度は、震度計内部で以下のようなディジタル処理によって計算されます。
2000年１０月６日に発生した鳥取県西部地震の米子市（計測震度＝５．１）を例に示します。

1. ディジタル加速度記録３成分（水平動２成分、上下動１成分）（図１）のそれぞれの フーリエ変換（図２）を求める。
2. 地震波の周期による影響を補正するフィルター（図３）を掛ける。
3. 逆フーリエ変換を行い、時刻歴の波形（図４）にもどす。
4. 得られたフィルター処理済みの３成分の波形をベクトル的に合成をする。
5. ベクトル波形（図５）の絶対値がある値 a 以上となる時間の合計を計算したとき、これがちょうど 0.3秒となるような a を求める。この例ではａ＝１２７．８５ｇａｌとなる。
6. 5.で求めたaを、　I = 2 log a + 0.94　により計測震度 Ｉ を計算する。計算されたIの小数第３位を四捨五入し、小数第２位を切り捨てたものを計測震度とする。

地震情報などにより発表される震度階級は、観測点における揺れの強さの程度を数値化した計測震度から換算されるものです。（表１）

表１　気象庁震度階級表

震度階級	計　測　震　度	震度階級	計　測　震　度
０	0.5未満	５弱	4.5以上5.0未満
１	0.5以上1.5未満	５強	5.0以上5.5未満
２	1.5以上2.5未満	６弱	5.5以上6.0未満
３	2.5以上3.5未満	６強	6.0以上6.5未満
４	3.5以上4.5未満	７	6.5以上

　図１　オリジナルの加速度波形



　図２　オリジナル加速度波形の周波数スペクトル



　図３　震度計算のためのフィルター特性



　図４　フィルター補正後の加速度波形



　図５　フィルター後の３成分合成加速度