JP5705692B2

JP5705692B2 - 地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法

Info

Publication number: JP5705692B2
Application number: JP2011210668A
Authority: JP
Inventors: 俊太野田; 俊六山本; 将宏是永; 直泰岩田; 賀章伊藤
Original assignee: Railway Technical Research Institute
Current assignee: Railway Technical Research Institute
Priority date: 2011-09-27
Filing date: 2011-09-27
Publication date: 2015-04-22
Anticipated expiration: 2031-09-27
Also published as: JP2013072685A

Description

本発明は、地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法に関するものである。

これまでに、地震のマグニチュードを推定するための方法として様々な方法が提案されている。

例えば、実体波（Ｐ波・Ｓ波）や表面波などの振幅の大きさを用いるもの（下記特許文献１参照）、長周期波形のスペクトル解析から求めるもの、地震動の継続時間を用いるものなどがある。

特開２００２−２７７５５７号公報

Ｍａｉ，Ｐ．Ｍ.,Ｐ．Ｓｐｕｄｉｃｈ，ａｎｄＪ．Ｂｏａｔｗｒｉｇｈｔ，ＨｙｐｏｃｅｎｔｅｒＬｏｃａｔｉｏｎｓｉｎＦｉｎｉｔｅ−ＳｏｕｒｃｅＲｕｐｔｕｒｅＭｏｄｅｌｓ，Ｂｕｌｌ．Ｓｅｉｓｍ．Ｓｏｃ．Ａｍ., Ｖｏｌ．９５，Ｎｏ．３，ｐｐ．９６５−９８０，２００５Ｉｉｏ，Ｙ., Ｅａｒｔｈｑｕａｋｅｎｕｃｌｅａｔｉｏｎｐｒｏｃｅｓｓ−Ｄｏｅｓｔｈｅｉｎｉｔｉａｔｉｏｎｏｆｅａｒｔｈｑｕａｋｅｒｕｐｔｕｒｅｋｎｏｗａｂｏｕｔｉｔｓｔｅｒｍｉｎａｔｉｏｎ？，Ｅａｒｔｈｑｕａｋｅｓ，Ｔｓｕｎａｍｉｓ，ａｎｄＶｏｌｃａｎｏｅｓｉｎＥｎｃｙｃｌｏｐｅｄｉａｏｆＣｏｍｐｌｅｘｉｔｙａｎｄＳｙｓｔｅｍｓＳｃｉｅｎｃｅ，ｐｐ．２５３８−２５５５，２００９

従来の地震のマグニチュードの推定方法では、それぞれ以下のような問題があった。

振幅の大きさを用いる方法（気象庁マグニチュードＭ_J、実体波マグニチュードｍ_B、表面波マグニチュードＭ_Sなど）では、大地震の場合、マグニチュードの値が飽和してしまう。

また、長周期波形のスペクトル解析から求める方法（モーメントマグニチュードＭ_W）では、時間がかかる。

更に、地震動の継続時間を用いるものも、時間がかかる上に、地震計特性や観測点周辺の環境に大きく依存してしまう。

本発明は、上記状況に鑑みて、大地震でもマグニチュードの値が飽和することなく、即時的に安定したマグニチュードを推定することができる、地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法を提供することを目的とする。

本発明は、上記目的を達成するために、
〔１〕地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法であって、地震動の初動の到達時刻とこの地震動の最大振幅となる時刻との時間差Ｔを求め、この時間差Ｔに基づき、あらかじめ統計的に求めておいたｌｏｇＴ＝ａ×Ｍ＋ｂ（ａ，ｂは統計解析によって決定される定数，Ｍはマグニチュード）からマグニチュードを推定する地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がＰ波であり、前記時間差ＴとしてＰ波初動の到達時刻と該Ｐ波の最大振幅となる時刻との時間差Ｔ _P を用いることを特徴とする。

〔２〕地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法であって、地震動の初動の到達時刻とこの地震動の最大振幅となる時刻との時間差Ｔを求め、この時間差Ｔに基づき、あらかじめ統計的に求めておいたｌｏｇＴ＝ａ×Ｍ＋ｂ＋ｃ×Ｒ（ａ，ｂ，ｃは統計解析によって決定される定数、Ｍはマグニチュード、Ｒは震央距離もしくは震源距離）からマグニチュードを推定する地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がＰ波であり、前記時間差ＴとしてＰ波初動の到達時刻と該Ｐ波の最大振幅となる時刻との時間差Ｔ _P を用いることを特徴とする。

〔３〕地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法であって、地震動の初動の到達時刻とこの地震動の最大振幅となる時刻との時間差Ｔを求め、この時間差Ｔに基づき、あらかじめ統計的に求めておいたｌｏｇＴ＝ａ×Ｍ＋ｂ（ａ，ｂは統計解析によって決定される定数，Ｍはマグニチュード）からマグニチュードを推定する地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がＳ波であり、前記時間差ＴとしてＳ波初動の到達時刻とこのＳ波の最大振幅となる時刻との時間差Ｔ _S を用いることを特徴とする。

〔４〕地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法であって、地震動の初動の到達時刻とこの地震動の最大振幅となる時刻との時間差Ｔを求め、この時間差Ｔに基づき、あらかじめ統計的に求めておいたｌｏｇＴ＝ａ×Ｍ＋ｂ＋ｃ×Ｒ（ａ，ｂ，ｃは統計解析によって決定される定数、Ｍはマグニチュード、Ｒは震央距離もしくは震源距離）からマグニチュードを推定する地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がＳ波であり、前記時間差ＴとしてＳ波初動の到達時刻とこのＳ波の最大振幅となる時刻との時間差Ｔ _S を用いることを特徴とする。

〔５〕地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法であって、地震動の初動の到達時刻とこの地震動の最大振幅となる時刻との時間差Ｔを求め、この時間差Ｔに基づき、イベント毎のＴの代表値（平均値ｏｒ中央値）を求め、あらかじめ統計的に求めておいたｌｏｇＴ＝ａ×Ｍ＋ｂ（ａ，ｂは統計解析によって決定される定数、Ｍはマグニチュード）からマグニチュードを推定する地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がＰ波であり、前記時間差ＴとしてＰ波初動の到達時刻とこのＰ波の最大振幅となる時刻との時間差Ｔ _P を用いることを特徴とする。

〔６〕地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法であって、地震動の初動の到達時刻とこの地震動の最大振幅となる時刻との時間差Ｔを求め、この時間差Ｔに基づき、イベント毎のＴの代表値（平均値ｏｒ中央値）を求め、あらかじめ統計的に求めておいたｌｏｇＴ＝ａ×Ｍ＋ｂ＋ｃ×Ｒ（ａ，ｂ，ｃは統計解析によって決定される定数、Ｍはマグニチュード、Ｒは震央距離もしくは震源距離）からマグニチュードを推定する地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がＰ波であり、前記時間差ＴとしてＰ波初動の到達時刻とこのＰ波の最大振幅となる時刻との時間差Ｔ _P を用いることを特徴とする。

〔７〕地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法であって、地震動の初動の到達時刻とこの地震動の最大振幅となる時刻との時間差Ｔを求め、この時間差Ｔに基づき、イベント毎のＴの代表値（平均値ｏｒ中央値）を求め、あらかじめ統計的に求めておいたｌｏｇＴ＝ａ×Ｍ＋ｂ（ａ，ｂは統計解析によって決定される定数、Ｍはマグニチュード）からマグニチュードを推定する地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がＳ波であり、前記時間差ＴとしてＳ波初動の到達時刻とこのＳ波の最大振幅となる時刻との時間差Ｔ _S を用いることを特徴とする。

〔８〕地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法であって、地震動の初動の到達時刻とこの地震動の最大振幅となる時刻との時間差Ｔを求め、この時間差Ｔに基づき、イベント毎のＴの代表値（平均値ｏｒ中央値）を求め、あらかじめ統計的に求めておいたｌｏｇＴ＝ａ×Ｍ＋ｂ＋ｃ×Ｒ（ａ，ｂ，ｃは統計解析によって決定される定数、Ｍはマグニチュード、Ｒは震央距離もしくは震源距離）からマグニチュードを推定する地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がＳ波であり、前記時間差ＴとしてＳ波初動の到達時刻とこのＳ波の最大振幅となる時刻との時間差Ｔ _S を用いることを特徴とする。

本発明によれば、大地震でもマグニチュードの値が飽和することなく、即時的に安定したマグニチュードを推定することができる。

本発明のマグニチュード推定方法を説明するための概念図である。本発明にかかるＳ波初動到達からＳ波最大振幅出現までの時間Ｔ_sとマグニチュードＭとの関係を示す図である。本発明にかかる真のＭと推定されたＭの関係（上記第６のマグニチュード推定方法に基づく）を示す図である。

本発明にかかる地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法は、観測される地震動のうち最大値を示す成分がアスペリティ付近で生成されることに起因して、地震動の初動の到達時刻とこの地震動の最大振幅となる時刻との時間差Ｔを求め、この時間差Ｔに基づき、あらかじめ統計的に求めておいたｌｏｇＴ＝ａ×Ｍ＋ｂ（ａ，ｂは統計解析によって決定される定数，Ｍはマグニチュード）からマグニチュードを推定する地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がＰ波であり、前記時間差ＴとしてＰ波初動の到達時刻と該Ｐ波の最大振幅となる時刻との時間差Ｔ _P を用いる。

以下、本発明の実施の形態について詳細に説明する。

まず、上記した非特許文献１によれば、破壊開始点（震源）がアスペリティ（ａｓｐｅｒｉｔｙ：強震動生成域）の外側に存在することが有意に多いことが指摘されている。

また、上記した非特許文献２によれば、より大きなアスペリティを破壊するためには、そこから離れた強度のより弱い場所から破壊が徐々に進行していく必要があることが指摘されている。したがって、観測される地震動のうち最大値を示す成分がアスペリティ付近で生成されたものであることを考慮すれば、上記の指摘より、地震波の初動の到達から地震波の最大振幅値が出現するまでの時間はマグニチュードにほぼ対応するといえる。

本発明では、これに着目し、Ｐ波またはＳ波の最大振幅値となる時刻の差を利用してマグニチュードを推定する。

図１は本発明のマグニチュード推定方法を説明するための概念図である。

図１において、Ｐ波初動の到達時刻をＴ_pr、Ｐ波初動の到達時刻Ｔ_prからＳ波初動の到達までをＰ相、Ｐ相での最大振幅の時刻をＴ_pmax、Ｔ_pmax−Ｔ_prをＴ_pとする。また、Ｓ波初動の到達時刻Ｔ_sr、Ｓ波初動の到達時刻Ｔ_sr以降をＳ相、Ｓ相での最大振幅の時刻をＴ_smax、Ｔ_smax−Ｔ_srをＴ_sとする。

図１に示した、Ｐ波初動到達からＰ相での最大振幅を示す時刻の差Ｔ_pとＳ波初動到達からＳ相での最大振幅を示す時刻の差Ｔ_sとが、それぞれマグニチュードＭとの高い相関を持つので、これらをあらかじめ統計的に求めておいた関数でフィッティングしてマグニチュードを推定する。
（１）第１のマグニチュード推定方法
まず、Ｐ波初動の到達時刻Ｔ_prを決定し、次いで、Ｐ相で最大振幅となる時刻Ｔ_pmaxを決定し、Ｐ波初動の到達時刻Ｔ_prとＰ相で最大振幅となる時刻Ｔ_pmaxの差をとり、求めた値をＴ_pとする。

そこで、あらかじめ統計的に求めておいたｌｏｇＴ_p＝ａ×Ｍ＋ｂ（ａ，ｂは統計解析によって決定される定数）からＭを推定する。

（２）第２のマグニチュード推定方法
まず、Ｓ波初動の到達時刻Ｔ_srを決定し、次いで、Ｓ相で最大振幅となる時刻Ｔ_smaxを決定し、Ｓ波初動の到達時刻Ｔ_srとＳ相で最大振幅となる時刻Ｔ_smaxの差をとり、求めた値をＴ_sとする。

そこで、あらかじめ統計的に求めておいたｌｏｇＴ_s＝ａ×Ｍ＋ｂ（ａ，ｂは統計解析によって決定される定数）からＭを推定する。

（３）第３のマグニチュード推定方法
まず、Ｐ波初動の到達時刻Ｔ_prを決定し、次いで、Ｐ相で最大振幅となる時刻Ｔ_pmaxを決定し、Ｐ波初動の到達時刻Ｔ_prとＰ相で最大振幅となる時刻Ｔ_pmaxの差をとり、求めた値をＴ_pとする。

そこで、あらかじめ統計的に求めておいたｌｏｇＴ_p＝ａ×Ｍ＋ｂ＋ｃ×Ｒ（ａ，ｂ，ｃは統計解析によって決定される定数、Ｒは震央距離もしくは震源距離）からＭを推定する。

（４）第４のマグニチュード推定方法
まず、Ｓ波初動の到達時刻Ｔ_srを決定し、次いで、Ｓ相で最大振幅となる時刻Ｔ_smaxを決定し、Ｓ波初動の到達時刻Ｔ_srとＳ相で最大振幅となる時刻Ｔ_smaxの差をとり、求めた値をＴ_sとする。

そこで、あらかじめ統計的に求めておいたｌｏｇＴ_s＝ａ×Ｍ＋ｂ＋ｃ×Ｒ（ａ，ｂ，ｃは統計解析によって決定される定数、Ｒは震央距離もしくは震源距離）からＭを推定する。

（５）第５のマグニチュード推定方法
まず、Ｐ波初動の到達時刻Ｔ_prを決定し、次いで、Ｐ相で最大振幅となる時刻Ｔ_pmaxを決定し、Ｐ波初動の到達時刻Ｔ_prとＰ相で最大振幅となる時刻Ｔ_pmaxの差をとり、求めた値をＴ_pとする。

そこで、イベント毎のＴ_pの代表値（平均値ｏｒ中央値）を求める。そして、あらかじめ統計的に求めておいたｌｏｇＴ_p＝ａ×Ｍ＋ｂ（ａ，ｂは統計解析によって決定される定数）からＭを推定する。

（６）第６のマグニチュード推定方法
まず、Ｓ波初動の到達時刻Ｔ_srを決定し、次いで、Ｓ相で最大振幅となる時刻Ｔ_smaxを決定し、Ｓ波初動の到達時刻Ｔ_srとＳ相で最大振幅となる時刻Ｔ_smaxの差をとり、求めた値をＴ_sとする。

そこで、イベント毎のＴ_sの代表値（平均値ｏｒ中央値）を求める。そして、あらかじめ統計的に求めておいたｌｏｇＴ_s＝ａ×Ｍ＋ｂ（ａ，ｂは統計解析によって決定される定数）からＭを推定する（図３参照）。

（７）第７のマグニチュード推定方法
まず、Ｐ波初動の到達時刻Ｔ_prを決定し、次いで、Ｐ相で最大振幅となる時刻Ｔ_pmaxを決定し、Ｐ波初動の到達時刻Ｔ_prとＰ相で最大振幅となる時刻Ｔ_pmaxの差をとり、求めた値をＴ_pとする。

そこで、イベント毎のＴ_pの代表値（平均値ｏｒ中央値）を求める。そして、あらかじめ統計的に求めておいたｌｏｇＴ_p＝ａ×Ｍ＋ｂ＋ｃ×Ｒ（ａ，ｂ，ｃは統計解析によって決定される定数、Ｒは震央距離もしくは震源距離）からＭを推定する。

（８）第８のマグニチュード推定方法
まず、Ｓ波初動の到達時刻Ｔ_srを決定し、次いで、Ｓ相で最大振幅となる時刻Ｔ_smaxを決定し、Ｓ波初動の到達時刻Ｔ_srとＳ相で最大振幅となる時刻Ｔ_smaxの差をとり、求めた値をＴ_sとする。

そこで、イベント毎のＴ_sの代表値（平均値ｏｒ中央値）を求める。そして、あらかじめ統計的に求めておいたｌｏｇＴ_s＝ａ×Ｍ＋ｂ＋ｃ×Ｒ（ａ，ｂ，ｃは統計解析によって決定される定数、Ｒは震央距離もしくは震源距離）からＭを推定する。

図２は本発明にかかるＳ波初動到達からＳ波最大振幅出現までの時間Ｔ_sとマグニチュードＭとの関係を示す図であり、横軸にマグニチュードＭ、縦軸にｌｏｇＴ_sを示している。ここでは、Ｋ−ＮＥＴのマグニチュード５以上のデータを用いて計算している。また、●は、各Ｍ毎のｌｏｇＴ_sの平均値、エラーバーａは標準偏差を表す。回帰分布によって決定された直線ｌｏｇＴ＝０．３５２×Ｍ−１．５３は、上記のｌｏｇＴ_sの平均値を用いて求めている。

なお、実際の解析では、表面波の影響を少なくするため、ハイパスフィルター等の処理を行うこともある。

図３は本発明にかかる真のＭ（Ｍ_W：モーメントマグニチュード）と推定されたＭの関係（上記第６のマグニチュード推定方法に基づく）を示す図であり、真のＭと推定されたＭの差のＲＭＳ（ＲｏｏｔＭｅａｎＳｑｕａｒｅ，二乗平均平方根）は０．５２となり、十分な精度を持つことがわかった。

なお、本発明は上記実施例に限定されるものではなく、本発明の趣旨に基づき種々の変形が可能であり、これらを本発明の範囲から排除するものではない。

本発明の地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法は、特に、大地震のマグニチュードの推定を即時的に安定して行うことができる簡易マグニチュード推定方法として利用することができる。

Claims

地震動の初動の到達時刻と該地震動の最大振幅となる時刻との時間差Ｔを求め、この時間差Ｔに基づき、あらかじめ統計的に求めておいたｌｏｇＴ＝ａ×Ｍ＋ｂ（ａ，ｂは統計解析によって決定される定数、Ｍはマグニチュード）からマグニチュードを推定する地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がＰ波であり、前記時間差ＴとしてＰ波初動の到達時刻と該Ｐ波の最大振幅となる時刻との時間差Ｔ _P を用いることを特徴とする地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法。
地震動の初動の到達時刻と該地震動の最大振幅となる時刻との時間差Ｔを求め、この時間差Ｔに基づき、あらかじめ統計的に求めておいたｌｏｇＴ＝ａ×Ｍ＋ｂ＋ｃ×Ｒ（ａ，ｂ，ｃは統計解析によって決定される定数、Ｍはマグニチュード、Ｒは震央距離もしくは震源距離）からマグニチュードを推定する地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がＰ波であり、前記時間差ＴとしてＰ波初動の到達時刻と該Ｐ波の最大振幅となる時刻との時間差Ｔ _P を用いることを特徴とする地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法。
地震動の初動の到達時刻と該地震動の最大振幅となる時刻との時間差Ｔを求め、この時間差Ｔに基づき、あらかじめ統計的に求めておいたｌｏｇＴ＝ａ×Ｍ＋ｂ（ａ，ｂは統計解析によって決定される定数、Ｍはマグニチュード）からマグニチュードを推定する地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がＳ波であり、前記時間差ＴとしてＳ波初動の到達時刻と該Ｓ波の最大振幅となる時刻との時間差Ｔ _S を用いることを特徴とする地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法。
地震動の初動の到達時刻と該地震動の最大振幅となる時刻との時間差Ｔを求め、この時間差Ｔに基づき、あらかじめ統計的に求めておいたｌｏｇＴ＝ａ×Ｍ＋ｂ＋ｃ×Ｒ（ａ，ｂ，ｃは統計解析によって決定される定数、Ｍはマグニチュード、Ｒは震央距離もしくは震源距離）からマグニチュードを推定する地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がＳ波であり、前記時間差ＴとしてＳ波初動の到達時刻と該Ｓ波の最大振幅となる時刻との時間差Ｔ _S を用いることを特徴とする地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法。
地震動の初動の到達時刻と該地震動の最大振幅となる時刻との時間差Ｔを求め、この時間差Ｔに基づき、イベント毎のＴの代表値（平均値ｏｒ中央値）を求め、あらかじめ統計的に求めておいたｌｏｇＴ＝ａ×Ｍ＋ｂ（ａ，ｂは統計解析によって決定される定数、Ｍはマグニチュード）からマグニチュードを推定する地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がＰ波であり、前記時間差ＴとしてＰ波初動の到達時刻と該Ｐ波の最大振幅となる時刻との時間差Ｔ _P を用いることを特徴とする地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法。
地震動の初動の到達時刻と該地震動の最大振幅となる時刻との時間差Ｔを求め、この時間差Ｔに基づき、イベント毎のＴの代表値（平均値ｏｒ中央値）を求め、あらかじめ統計的に求めておいたｌｏｇＴ＝ａ×Ｍ＋ｂ＋ｃ×Ｒ（ａ，ｂ，ｃは統計解析によって決定される定数、Ｍはマグニチュード、Ｒは震央距離もしくは震源距離）からマグニチュードを推定する地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がＰ波であり、前記時間差ＴとしてＰ波初動の到達時刻と該Ｐ波の最大振幅となる時刻との時間差Ｔ _P を用いることを特徴とする地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法。
地震動の初動の到達時刻と該地震動の最大振幅となる時刻との時間差Ｔを求め、この時間差Ｔに基づき、イベント毎のＴの代表値（平均値ｏｒ中央値）を求め、あらかじめ統計的に求めておいたｌｏｇＴ＝ａ×Ｍ＋ｂ（ａ，ｂは統計解析によって決定される定数、Ｍはマグニチュード）からマグニチュードを推定する地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がＳ波であり、前記時間差ＴとしてＳ波初動の到達時刻と該Ｓ波の最大振幅となる時刻との時間差Ｔ _S を用いることを特徴とする地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法。
地震動の初動の到達時刻と該地震動の最大振幅となる時刻との時間差Ｔを求め、この時間差Ｔに基づき、イベント毎のＴの代表値（平均値ｏｒ中央値）を求め、あらかじめ統計的に求めておいたｌｏｇＴ＝ａ×Ｍ＋ｂ＋ｃ×Ｒ（ａ，ｂ，ｃは統計解析によって決定される定数、Ｍはマグニチュード、Ｒは震央距離もしくは震源距離）からマグニチュードを推定する地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がＳ波であり、前記時間差ＴとしてＳ波初動の到達時刻と該Ｓ波の最大振幅となる時刻との時間差Ｔ _S を用いることを特徴とする地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法。