参照元
ホーム>各種データ・資料>地震の活動状況>強震観測>用語の解説>計測震度の算出方法
　計測震度は、震度計内部で以下のようなディジタル処理によって計算されます。
2000年１０月６日に発生した鳥取県西部地震の米子市（計測震度＝５．１）を例に示します。
　1.ディジタル加速度記録３成分（水平動２成分、上下動１成分）（図１）のそれぞれのフーリエ変換（図２）を求める。
　2.地震波の周期による影響を補正するフィルター（図３）を掛ける。
　3.逆フーリエ変換を行い、時刻歴の波形（図４）にもどす。
　4.得られたフィルター処理済みの３成分の波形をベクトル的に合成をする。
　5.ベクトル波形（図５）の絶対値がある値a以上となる時間の合計を計算したとき、これがちょうど0.3秒となるようなaを求める。この例ではａ＝１２７．８５ｇａｌとなる。
　6. 5.で求めたaを、　I = 2 log a + 0.94　により計測震度 Ｉ を計算する。計算されたIの小数第３位を四捨五入し、小数第２位を切り捨てたものを計測震度とする。

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


元の加速度波形の図
　図１　オリジナルの加速度波形
http://www.data.jma.go.jp/svd/eqev/data/kyoshin/kaisetsu/image/orig_wave.png

オリジナル加速度波形の周波数スペクトル
　図２　オリジナル加速度波形の周波数スペクトル
http://www.data.jma.go.jp/svd/eqev/data/kyoshin/kaisetsu/image/orig_spec.png

震度計算のためのフィルター特性
　図３　震度計算のためのフィルター特性
http://www.data.jma.go.jp/svd/eqev/data/kyoshin/kaisetsu/image/filter_s.png

フィルター補正後の加速度波形
　図４　フィルター補正後の加速度波形
http://www.data.jma.go.jp/svd/eqev/data/kyoshin/kaisetsu/image/filter_wave.png

フィルター後の３成分合成加速度
　図５　フィルター後の３成分合成加速度
http://www.data.jma.go.jp/svd/eqev/data/kyoshin/kaisetsu/image/gosei.png