"地震動の初動の到達時刻と該地震動の最大振幅となる時刻との時間差Tを求め、この時間差Tに基づき、あらかじめ統計的に求めておいたlogT=a×Mba,bは統計解析によって決定される定数、Mはマグニチュードからマグニチュードを推定することを特徴とする地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法。地震動の初動の到達時刻と該地震動の最大振幅となる時刻との時間差Tを求め、この時間差Tに基づき、あらかじめ統計的に求めておいたlogT=a×Mbc×Ra,b,cは統計解析によって決定される定数、Mはマグニチュード、Rは震央距離もしくは震源距離からマグニチュードを推定することを特徴とする地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法。請求項1又は2記載の地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がP波であり、前記時間差TとしてP波初動の到達時刻と該P波の最大振幅となる時刻との時間差TPを用いることを特徴とする地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法。請求項1又は2記載の地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がS波であり、前記時間差TとしてS波初動の到達時刻と該S波の最大振幅となる時刻との時間差TSを用いることを特徴とする地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法。地震動の初動の到達時刻と該地震動の最大振幅となる時刻との時間差Tを求め、この時間差Tに基づき、イベント毎のTの代表値平均値or中央値を求め、あらかじめ統計的に求めておいたlogT=a×Mba,bは統計解析によって決定される定数、Mはマグニチュードからマグニチュードを推定することを特徴とする地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法。地震動の初動の到達時刻と該地震動の最大振幅となる時刻との時間差Tを求め、この時間差Tに基づき、イベント毎のTの代表値平均値or中央値を求め、あらかじめ統計的に求めておいたlogT=a×Mbc×Ra,b,cは統計解析によって決定される定数、Mはマグニチュード、Rは震央距離もしくは震源距離からマグニチュードを推定することを特徴とする地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法。請求項5又は6記載の地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がP波であり、前記時間差TとしてP波初動の到達時刻と該P波の最大振幅となる時刻との時間差TPを用いることを特徴とする最大振幅値出現までの時間を利用した簡易マグニチュード推定方法。請求項5又は6記載の地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がS波であり、前記時間差TとしてS波初動の到達時刻と該S波の最大振幅となる時刻との時間差TSを用いることを特徴とする最大振幅値出現までの時間を利用した簡易マグニチュード推定方法。従来の地震のマグニチュードの推定方法では、それぞれ以下のような問題があった。振幅の大きさを用いる方法気象庁マグニチュードMJ、実体波マグニチュードmB、表面波マグニチュードMSなどでは、大地震の場合、マグニチュードの値が飽和してしまう。また、長周期波形のスペクトル解析から求める方法モーメントマグニチュードMWでは、時間がかかる。更に、地震動の継続時間を用いるものも、時間がかかる上に、地震計特性や観測点周辺の環境に大きく依存してしまう。本発明は、上記状況に鑑みて、大地震でもマグニチュードの値が飽和することなく、即時的に安定したマグニチュードを推定することができる、地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法を提供することを目的とする。本発明は、上記目的を達成するために、1地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、地震動の初動の到達時刻とこの地震動の最大振幅となる時刻との時間差Tを求め、この時間差Tに基づき、あらかじめ統計的に求めておいたlogT=a×Mba,bは統計解析によって決定される定数,Mはマグニチュードからマグニチュードを推定することを特徴とする。2地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、地震動の初動の到達時刻とこの地震動の最大振幅となる時刻との時間差Tを求め、この時間差Tに基づき、あらかじめ統計的に求めておいたlogT=a×Mbc×Ra,b,cは統計解析によって決定される定数、Mはマグニチュード、Rは震央距離もしくは震源距離からマグニチュードを推定することを特徴とする。3上記1又は2記載の地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がP波であり、前記時間差TとしてP波初動の到達時刻とこのP波の最大振幅となる時刻との時間差TPを用いることを特徴とする。4上記1又は2記載の地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がS波であり、前記時間差TとしてS波初動の到達時刻とこのS波の最大振幅となる時刻との時間差TSを用いることを特徴とする。5地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、地震動の初動の到達時刻とこの地震動の最大振幅となる時刻との時間差Tを求め、この時間差Tに基づき、イベント毎のTの代表値平均値or中央値を求め、あらかじめ統計的に求めておいたlogT=a×Mba,bは統計解析によって決定される定数、Mはマグニチュードからマグニチュードを推定することを特徴とする。6地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、地震動の初動の到達時刻とこの地震動の最大振幅となる時刻との時間差Tを求め、この時間差Tに基づき、イベント毎のTの代表値平均値or中央値を求め、あらかじめ統計的に求めておいたlogT=a×Mbc×Ra,b,cは統計解析によって決定される定数、Mはマグニチュード、Rは震央距離もしくは震源距離からマグニチュードを推定することを特徴とする。7上記5又は6記載の地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がP波であり、前記時間差TとしてP波初動の到達時刻とこのP波の最大振幅となる時刻との時間差TPを用いることを特徴とする。8上記5又は6記載の地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、前記地震動がS波であり、前記時間差TとしてS波初動の到達時刻とこのS波の最大振幅となる時刻との時間差TSを用いることを特徴とする。本発明によれば、大地震でもマグニチュードの値が飽和することなく、即時的に安定したマグニチュードを推定することができる。地震の最大振幅値出現までの時間を利用した簡易マグニチュード推定方法において、観測される地震動のうち最大値を示す成分がアスペリティ付近で生成されることに起因して、地震動の初動の到達時刻とこの地震動の最大振幅となる時刻との時間差Tを求め、この時間差Tに基づき、あらかじめ統計的に求めておいたlogT=a×Mba,bは統計解析によって決定される定数,Mはマグニチュードからマグニチュードを推定する。以下、本発明の実施の形態について詳細に説明する。まず、上記した非特許文献1によれば、破壊開始点震源がアスペリティasperity:強震動生成域の外側に存在することが有意に多いことが指摘されている。また、上記した非特許文献2によれば、より大きなアスペリティを破壊するためには、そこから離れた強度のより弱い場所から破壊が徐々に進行していく必要があることが指摘されている。したがって、観測される地震動のうち最大値を示す成分がアスペリティ付近で生成されたものであることを考慮すれば、上記の指摘より、地震波の初動の到達から地震波の最大振幅値が出現するまでの時間はマグニチュードにほぼ対応するといえる。本発明では、これに着目し、P波またはS波の最大振幅値となる時刻の差を利用してマグニチュードを推定する。図1は本発明のマグニチュード推定方法を説明するための概念図である。図1において、P波初動の到達時刻をTpr、P波初動の到達時刻TprからS波初動の到達までをP相、P相での最大振幅の時刻をTpmax、Tpmax-TprをTpとする。また、S波初動の到達時刻Tsr、S波初動の到達時刻Tsr以降をS相、S相での最大振幅の時刻をTsmax、Tsmax-TsrをTsとする。図1に示した、P波初動到達からP相での最大振幅を示す時刻の差TpとS波初動到達からS相での最大振幅を示す時刻の差Tsとが、それぞれマグニチュードMとの高い相関を持つので、これらをあらかじめ統計的に求めておいた関数でフィッティングしてマグニチュードを推定する。1第1のマグニチュード推定方法まず、P波初動の到達時刻Tprを決定し、次いで、P相で最大振幅となる時刻Tpmaxを決定し、P波初動の到達時刻TprとP相で最大振幅となる時刻Tpmaxの差をとり、求めた値をTpとする。そこで、あらかじめ統計的に求めておいたlogTp=a×Mba,bは統計解析によって決定される定数からMを推定する。2第2のマグニチュード推定方法まず、S波初動の到達時刻Tsrを決定し、次いで、S相で最大振幅となる時刻Tsmaxを決定し、S波初動の到達時刻TsrとS相で最大振幅となる時刻Tsmaxの差をとり、求めた値をTsとする。そこで、あらかじめ統計的に求めておいたlogTs=a×Mba,bは統計解析によって決定される定数からMを推定する。3第3のマグニチュード推定方法まず、P波初動の到達時刻Tprを決定し、次いで、P相で最大振幅となる時刻Tpmaxを決定し、P波初動の到達時刻TprとP相で最大振幅となる時刻Tpmaxの差をとり、求めた値をTpとする。そこで、あらかじめ統計的に求めておいたlogTp=a×Mbc×Ra,b,cは統計解析によって決定される定数、Rは震央距離もしくは震源距離からMを推定する。4第4のマグニチュード推定方法まず、S波初動の到達時刻Tsrを決定し、次いで、S相で最大振幅となる時刻Tsmaxを決定し、S波初動の到達時刻TsrとS相で最大振幅となる時刻Tsmaxの差をとり、求めた値をTsとする。そこで、あらかじめ統計的に求めておいたlogTs=a×Mbc×Ra,b,cは統計解析によって決定される定数、Rは震央距離もしくは震源距離からMを推定する。5第5のマグニチュード推定方法まず、P波初動の到達時刻Tprを決定し、次いで、P相で最大振幅となる時刻Tpmaxを決定し、P波初動の到達時刻TprとP相で最大振幅となる時刻Tpmaxの差をとり、求めた値をTpとする。そこで、イベント毎のTpの代表値平均値or中央値を求める。そして、あらかじめ統計的に求めておいたlogTp=a×Mba,bは統計解析によって決定される定数からMを推定する。6第6のマグニチュード推定方法まず、S波初動の到達時刻Tsrを決定し、次いで、S相で最大振幅となる時刻Tsmaxを決定し、S波初動の到達時刻TsrとS相で最大振幅となる時刻Tsmaxの差をとり、求めた値をTsとする。そこで、イベント毎のTsの代表値平均値or中央値を求める。そして、あらかじめ統計的に求めておいたlogTs=a×Mba,bは統計解析によって決定される定数からMを推定する図3参照。7第7のマグニチュード推定方法まず、P波初動の到達時刻Tprを決定し、次いで、P相で最大振幅となる時刻Tpmaxを決定し、P波初動の到達時刻TprとP相で最大振幅となる時刻Tpmaxの差をとり、求めた値をTpとする。そこで、イベント毎のTpの代表値平均値or中央値を求める。そして、あらかじめ統計的に求めておいたlogTp=a×Mbc×Ra,b,cは統計解析によって決定される定数、Rは震央距離もしくは震源距離からMを推定する。8第8のマグニチュード推定方法まず、S波初動の到達時刻Tsrを決定し、次いで、S相で最大振幅となる時刻Tsmaxを決定し、S波初動の到達時刻TsrとS相で最大振幅となる時刻Tsmaxの差をとり、求めた値をTsとする。そこで、イベント毎のTsの代表値平均値or中央値を求める。そして、あらかじめ統計的に求めておいたlogTs=a×Mbc×Ra,b,cは統計解析によって決定される定数、Rは震央距離もしくは震源距離からMを推定する。図2は本発明にかかるS波初動到達からS波最大振幅出現までの時間TsとマグニチュードMとの関係を示す図であり、横軸にマグニチュードM、縦軸にlogTsを示している。ここでは、K-NETのマグニチュード5以上のデータを用いて計算している。また、●は、各M毎のlogTsの平均値、エラーバーaは標準偏差を表す。回帰分布によって決定された直線logT=0.352×M-1.53は、上記のlogTsの平均値を用いて求めている。なお、実際の解析では、表面波の影響を少なくするため、ハイパスフィルター等の処理を行うこともある。図3は本発明にかかる真のMMW:モーメントマグニチュードと推定されたMの関係上記第6のマグニチュード推定方法に基づくを示す図であり、真のMと推定されたMの差のRMSRootMeanSquare,二乗平均平方根は0.52となり、十分な精度を持つことがわかった。なお、本発明は上記実施例に限定されるものではなく、本発明の趣旨に基づき種々の変形が可能であり、これらを本発明の範囲から排除するものではない。"