地震震度計算

以地震觀測資料計算出中央氣象局（CWB）、日本氣象廳（JMA）、福建省地震局（FJEA）所使用的震度。

中央氣象局地震震度計算與繪圖的整合版本請參考 https://github.com/chemars/Seismic-Intensity

準備步驟

1. 安裝 GNU Octave（計算程式適用於 Windows-64 版）。
2. 下載所有檔案，並另外下載觀測資料，觀測資料需包括三軸加速度訊號（單位為 gal）及取樣率（單位為 Hz），把所有檔案放在同一目錄。
3. 開啟 GNU Octave，切換到檔案所在的目錄。

計算

計算程式皆為函數（function），呼叫格式為函數名稱(觀測資料檔案名稱,分隔符號,要忽略的列數,南北向行的位置,東西向行的位置,上下向行的位置,取樣率)，cwb2000無取樣率參數。cwb2020輸出為字串，其餘輸出為數值。

• cwb2000對應的震度為中央氣象局於2000年8月1日公告修訂的震度，分為0級至7級。
• cwb2020對應的震度為中央氣象局於2020年1月1日起施行的震度，分為0級、1級、2級、3級、4級、5弱、5強、6弱、6強、7級。
• jma對應的震度為日本氣象廳的計測震度，帶有一位小數，程式碼修改自Rで計測震度の算出。
• fjea對應的震度為福建省地震局的地震儀器烈度，分為1度至12度，程式碼修改自Rで計測震度の算出。

另有震度轉換程式（亦為函數），呼叫格式為intensity_scale(震度種類,震度數值)，輸出為字串。

• jma是將日本氣象廳的計測震度轉換為震度階級，日本氣象廳震度階級分為0、1、2、3、4、5弱、5強、6弱、6強、7。

範例

地震觀測資料來源：中央氣象局 個案地震報導、日本氣象廳 強震観測データ

以中央氣象局的檔案為例，可使用下列方式計算出震度（jma 可更改為 fjea）

>> cwb2000("T073001.263.txt","",11,3,4,2)
ans =  6
>> jma("T073001.263.txt","",11,3,4,2,200)
ans =  5.2000
>> cwb2020("T073001.263.txt","",11,3,4,2,200)
ans = 6-

以日本氣象廳的檔案為例，可使用下列方式計算出震度（jma 可更改為 fjea）

>> cwb2000("AA06EA01.csv",",",7,1,2,3)
ans =  6
>> jma("AA06EA01.csv",",",7,1,2,3,100)
ans =  5.1000
>> cwb2020("AA06EA01.csv",",",7,1,2,3,100)
ans = 5-

以日本氣象廳的計測震度為例，可使用下列方式轉換為震度階級

>> intensity_scale("jma",5.1)
ans = 5+

聲明

程式計算、轉換的方式與官方實際做法有所差異，得出的結果僅供參考，實際震度應以官方發布的為準。jma計算出的數值與日本氣象廳公布的計測震度可能有0.1的誤差。fjea沒有福建省地震局的資料可供驗證，誤差大小未知。

授權條款

MIT 授權條款

Rで計測震度の算出內的程式碼的授權條款為 CC0（相關說明）

參考資料

• Rで計測震度の算出，GitHub repository
• 交通部中央氣象局有感地震報告發布作業要點（中華民國103年12月22日中象地字第1030015380號函修正）
• 交通部中央氣象局有感地震報告發布作業要點（中華民國108年12月2日中象地字第1080016036號函修正）
• 計測震度の算出方法
• 気象庁震度階級の解説
• 地震仪器烈度表