(57)【要約】【課題】地震発生時に、想定点において必要とされる周波数帯の地震動波形を高い精度でリアルタイムに推定することができる地震動波形推定方法を提供する。【解決手段】地震発生予測地域Ａ１で観測された過去の地震動波形からリアルタイムで第１のフィルタを介して想定点において対象となる周波数帯の地震動波形を抽出し、想定点Ａ３（Ａ４）で観測された上記過去の地震動波形から第２のフィルタを介して周波数帯の地震動波形を抽出し、得られた周波数帯の地震動波形から地震発生予測地域と想定点との間の伝達関数を求め、第１のフィルタの周波数特性に当該伝達関数を掛け合わせた周波数伝達特性を有する地震動予測フィルタを作成して、地震発生時に地震発生予測地域で観測された地震動波形を地震動予測フィルタに入力することにより、想定点における上記周波数帯の地震動波形をリアルタイムで推定する。【選択図】図１
2010175274.tif 000002
【特許請求の範囲】
【請求項１】
過去の地震時に地震発生予測地域で観測された地震動波形と上記地震発生予測地域から離れた特定の想定点において観測された地震動波形とに基づいて、新たな地震発生時に上記想定点における地震動波形をリアルタイムに推定する地震動波形推定方法において、
上記過去の地震時における上記地震動波形を得るに際して、上記地震発生予測地域で観測された上記地震動波形から、リアルタイムで第１のフィルタを介して上記想定点において上記新たな地震発生時に到達する地震動の大きさを判断するための対象となる周波数帯の地震動波形を抽出し、上記想定点で観測された上記地震動波形から当該地震動波形の形状および位相にずれを生じない第２のフィルタを介して上記周波数帯の地震動波形を抽出し、これら地震発生予測地域および想定点で得られた上記周波数帯の地震動波形から上記地震発生予測地域と上記想定点との間の伝達関数を求めた後に、上記第１のフィルタの周波数特性に当該伝達関数を掛け合わせた周波数伝達特性を有する地震動予測フィルタを作成し、
上記新たな地震発生時に、上記地震発生予測地域で観測された上記地震動波形をリアルタイムに上記地震動予測フィルタに入力することにより、上記想定点における上記周波数帯の地震動波形をリアルタイムで推定することを特徴とする地震動波形推定方法。
【請求項２】
上記第１のフィルタは、観測された上記地震動波形からリアルタイムで周波数帯の地震動波形を抽出可能なデジタルフィルタまたはアナログフィルタであり、かつ上記第２のフィルタは、高速フーリエ変換機能を用いたバンドパスフィルタであることを特徴とする請求項１に記載の地震動波形推定方法。
【請求項３】
予め、一箇所の上記想定点と複数の上記地震発生予測地域との間に対して、各々対応する複数の上記地震動予測フィルタを作成し、上記新たな地震発生時に、最も近接した上記地震発生予測地域と上記想定点との間の上記地震動予測フィルタを用いことを特徴とする請求項１または２に記載の地震動波形推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、地震が発生した際に、その地震動が特定の想定点に到達する前に、当該想定点における地震動波形をリアルタイムに、かつ高い精度で推定することができる地震動波形推定方法に関するものである。
【背景技術】
【０００２】
現在、地震発生時に、震源の近くで検知されたＰ波によって当該地震に関するリアルタイム地震情報を即時的に得る各種のシステムが開発されている。
そして、震源から数Ｋｍ以上離れた地点においては、Ｓ波に起因する地震動の主要動が到達するまでに数秒から数十秒の余裕があることから、上記リアルタイム地震情報を利用することにより、地震発生後数秒程度で、発生した地震のマグニチュードや震源位置等に関する情報を受け取って、地震被害の発生を防止することが可能になっている。
【０００３】
しかしながら、現在運用されているリアルタイム地震情報は、発生した地震のマグニチュード、震源位置および震源深さ等の情報に限られているため、地震規模のリアルタイム推定は可能になったものの、発生した地震により、特定の想定地点において、地震動の変化、ピークの到来時期、あるいは継続時間といった地震動波形の詳細を地震到達前に知ることができないという問題点があった。
【０００４】
そこで、本発明者等は、先に下記特許文献１において、地震発生時に、その地震動が特定の想定点に到達する前に、当該想定点における地震動波形をリアルタイムに推定することが可能となる地震動波形推定方法を提案した。
【０００５】
この地震動波形推定方法は、予め複数の地震発生予測地域に対する当該予測地域と上記想定点との間の地震観測点と、上記想定点との間の地震動の伝達関数を求めておき、地震発生時に検知された震源位置に係るリアルタイム地震情報から上記当該震源位置に最も近い上記震動発生予測地域に対する上記伝達関数を抽出し、次いで上記観測点で得られた地震動波形の情報をリアルタイムに取得して、上記伝達関数に基づいて上記想定点における地震動波形をリアルタイムに推定することを特徴とするものである。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００６−１９４６１９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
ところで、一般に、このような伝達関数を用いて、地震発生時に地震発生予測地域で観測された地震動波形から想定点における地震動波形を推定する方法においては、上記地震発生予測地域において観測された地震動波形に、地震計における感度等の特性に依存した周波数帯の地震動成分が含まれている。
【０００８】
このため、上記伝達関数を用いて得られた上記想定点における推定の地震動波形には、地震発生時に、当該想定点において防災のための構造制御や機器制御を行う上で、判断の対象となる周波数帯に加えて、それ以外の必ずしも上記判断に必要ではない短周期や長周期の地震動成分やノイズが含まれてしまうという問題点があった。
【０００９】
本発明は、かかる事情に鑑みてなされたものであり、地震発生時に、その地震動が特定の想定点に到達する前に、当該想定点において必要とされる周波数帯の地震動波形を高い精度でリアルタイムに推定することができる地震動波形推定方法を提供することを課題とするものである。
【課題を解決するための手段】
【００１０】
上記課題を解決するため、請求項１に記載の発明は、過去の地震時に地震発生予測地域で観測された地震動波形と上記地震発生予測地域から離れた特定の想定点において観測された地震動波形とに基づいて、新たな地震発生時に上記想定点における地震動波形をリアルタイムに推定する地震動波形推定方法において、上記過去の地震時における上記地震動波形を得るに際して、上記地震発生予測地域で観測された上記地震動波形から、リアルタイムで第１のフィルタを介して上記想定点において上記新たな地震発生時に到達する地震動の大きさを判断するための対象となる周波数帯の地震動波形を抽出し、上記想定点で観測された上記地震動波形から当該地震動波形の形状および位相にずれを生じない第２のフィルタを介して上記周波数帯の地震動波形を抽出し、これら地震発生予測地域および想定点で得られた上記周波数帯の地震動波形から上記地震発生予測地域と上記想定点との間の伝達関数を求めた後に、上記第１のフィルタの周波数特性に当該伝達関数を掛け合わせた周波数伝達特性を有する地震動予測フィルタを作成し、上記新たな地震発生時に、上記地震発生予測地域で観測された上記地震動波形をリアルタイムに上記地震動予測フィルタに入力することにより、上記想定点における上記周波数帯の地震動波形をリアルタイムで推定することを特徴とするものである。
【００１１】
ここで、地震発生予測地域とは、各種の地震調査によって、今後地震が起きると推定される予想震源地およびその近傍の地域をいう。
また、請求項２に記載の発明は、請求項１に記載の発明において、上記第１のフィルタが、観測された上記地震動波形からリアルタイムで周波数帯の地震動波形を抽出可能なデジタルフィルタまたはアナログフィルタであり、かつ上記第２のフィルタが、高速フーリエ変換機能（Fast Fourier Transform：以下FFT）を用いたバンドパスフィルタであることを特徴とするものである。
【００１２】
さらに、請求項３に記載の発明は、請求項１または２に記載の発明において、予め、一箇所の上記想定点と複数の上記地震発生予測地域との間に対して、各々対応する複数の上記地震動予測フィルタを作成し、上記新たな地震発生時に、最も近接した上記地震発生予測地域と上記想定点との間の上記地震動予測フィルタを用いることを特徴とするものである。
【発明の効果】
【００１３】
請求項１〜３のいずれかに記載の発明においては、地震発生予測地域で観測された地震動波形から、リアルタイムで第１のフィルタを介して対象となる周波数帯の地震動波形を抽出する。すると、時々刻々と上記第１のフィルタによって抽出された地震動波形には、位相特性の変化が生じ、波形に時間遅れが生じたり、波形形状の山谷の継時変化特性が変化したりして、実際に観測された地震動波形とは異なった歪んだ特性が上記対象となる周波数帯成分に与えられてしまう。
【００１４】
他方、伝達関数を同定するための想定点で観測された地震動波形に対しては、リアルタイムでフィルタ処理を行う必要がないために、予め地震の全継続時間の時刻歴データを用いて、当該地震動波形の形状および位相にずれを生じない第２のフィルタを介して上記周波数帯の地震動波形を抽出する。
【００１５】
次いで、地震発生予測地域における地震動波形から第１のフィルタで抽出された対象とする周波数帯の地震動波形と、想定点における地震動波形から第２のフィルタで抽出された対象とする周波数帯の地震動波形を用いて、伝達関数を作成する。この伝達関数は、例えば上記特許文献１に開示された解析方法や、部分空間法などを用いて適当な次数の伝達関数モデルを同定することにより作成することができる。そして、第１のフィルタの周波数特性に、作成された上記伝達関数を掛け合わせた周波数伝達特性を有する地震動予測フィルタを作成する。
【００１６】
このようにして得られた地震動予測フィルタを用いることにより、第１のフィルタで抽出する際に生じた移動特性の変化を補正することができる。
そこで、新たな地震が発生した際に、上記地震発生予測地域で観測された地震動波形をリアルタイムに上記地震動予測フィルタに入力することにより、想定点において対象となる周波数帯の地震動波形をリアルタイムで、かつ高い精度で推定することができる。
【００１７】
ここで、請求項２に記載の発明のように、第１のフィルタとしては、観測された地震動波形からリアルタイムで周波数帯の地震動波形を抽出可能なデジタルフィルタまたはアナログフィルタが好適であり、第２のフィルタとしては、FFTを用いたバンドパスフィルタが好適である。
【００１８】
さらに、請求項３に記載の発明によれば、一箇所の想定点に対して、複数の地震発生予測地域を想定して、各々の地震発生予測地域との間の地震動予測フィルタを作成しておくことにより、新たな地震がいずれかの地震発生予測地域で発生した場合にも、確実に上記想定点において対象となる周波数帯の地震動波形をリアルタイムで推定することが可能になる。
【図面の簡単な説明】
【００１９】
【図１】本発明の一実施例において用いた過去の地震と地震発生予測地域および想定点の位置を示す地図である。
【図２】本実施例において観測された地震動波形を示すものであり、（ａ）は地震発生予測地域Ａ１における波形図、（ｂ）は想定点Ａ３おける波形図、（ｃ）は想定点Ａ４おける波形図である。
【図３】図１の地震発生予測地域で用いたフィルタの周波数特性を示す図である。
【図４】図２の地震動波形からフィルタによって抽出された地震動波形を示すものであり、（ａ）は地震発生予測地域Ａ１における波形図、（ｂ）は想定点Ａ３おける波形図、（ｃ）は想定点Ａ４おける波形図である。
【図５】本実施例における地震動予測フィルタの周波数特性を示す図である。
【図６】上記地震動予測フィルタを用いて想定点Ａ４における地震動を推定した結果を示すものであり、（ａ）は前震での検証結果、（ｂ）は本震での検証結果、（ｃ）は余震での検証結果である。
【発明を実施するための形態】
【００２０】
以下、本発明に係る地震動波形推定方法の一実施形態について説明する。
先ず、今後地震が起きると推定される震源地またはその近傍の地震発生予測地域と、当該地震発生予測地域から離れた特定の想定点を特定する。そして、上記震源地において過去に発生した地震時に観測された上記地震発生予測地域および想定点での地震動波形に基づいて、上記地震発生予測地域と想定点との間の伝達関数を同定する。
【００２１】
この際に、先ず上記地震発生予測地域で観測された地震動波形から、リアルタイムでデジタルフィルタまたはアナログフィルタ（第１のフィルタ）によって、上記想定点において新たな地震発生時に到達する地震動の大きさを判断するための対象となる周波数帯の地震動波形を抽出する。これにより、上記地震発生予測地域で観測された地震動波形に含まれている地震計の感度等に依存した周波数帯の地震動成分やノイズが除去される。
【００２２】
ちなみに、このようにして得られた上記周波数帯の地震動波形は、時々刻々上記フィルタから抽出されるために、波形に時間遅れ（位相のずれ）が生じたり、あるいは波形形状の継時変化特性が変化したりして、歪んだ周波数特性が周波数成分に与えられてしまっている。
【００２３】
他方、上記想定点で観測された地震動波形からも、フィルタを介して同様の周波数帯の地震動波形を抽出する。この際には、上記地震動波形について全継続時間の時刻歴データ（例えば、地震動が６０秒間続いた場合には、当該６０秒間の地震動波形）を用いることができるために、FFTを用いたバンドパスフィルタ（第２のフィルタ）を用いる。これにより、上記想定点で観測された地震動波形から、波形形状および位相にずれを生じない上記周波数帯の地震動波形を抽出することができる。
【００２４】
次いで、これらの地震発生予測地域および想定点で観測された地震動波形から、フィルタを介して抽出された対象とする周波数帯の地震動波形を用いて、上記地震発生予測地域と想定点との間の伝達関数を求める。この伝達関数は、例えば部分空間法を用いて同定することができる。この部分空間法は、伝達特性を下式（１）、（２）の状態方程式で表現し、Ａ、Ｂ、Ｃ、Ｄの行列を同定する方法である。
【００２５】
【数１】
2010175274.tif 000003
【００２６】
【数２】
2010175274.tif 000004
【００２７】
ここで、u(i)は、地震発生予測地域における時刻ステップiで観測された地震動波形値であり、y(i)は、想定点における時刻ステップiで観測された地震動波形値である。また、X(i)は、状態ベクトルであり、伝達特性のモデル次数と等しい次元を持つ。
【００２８】
次に、地震発生予測地域の地震動波形から対象となる周波数帯の地震動波形を抽出する際に用いたリアルタイムでデジタルフィルタまたはアナログフィルタ（第１のフィルタ）の周波数特性に得られた伝達関数を掛け合わせて地震動予測フィルタを作成する。
【００２９】
そして、地震が発生した時に、上記地震発生予測地域で観測された地震動波形を、上記地震動予測フィルタに入力することにより、想定点における上記周波数帯の地震動波形をリアルタイムで推定する。
【００３０】
以上のように、上記構成からなる地震動波形推定方法によれば、地震が発生した際に、上記地震発生予測地域で観測された地震動波形を上記地震動予測フィルタに入力することにより、想定点において対象となる周波数帯の地震動波形をリアルタイムで、かつ高い精度で推定することができる。
【００３１】
この結果、第１および第２のフィルタによって抽出する「想定点において対象となる周波数帯」を、利用する制御対象の卓越周波数によって適宜設定することにより、上記地震動波形推定方法を用いて推定した地震動波形に基づいて、より高度な機器制御や建物の制御を行うことが可能になる。
例えば、超高層ビルや免震構造物の制御に利用する場合は、これらの構造物の卓越周波数が存在する０．５Ｈｚ〜０．１６Ｈｚ（２秒〜６秒）程度の範囲を含む周波数帯を対象とすればよい。また、１Ｈｚの卓越周波数を持つ機器の場合は、１Ｈｚを含む前後の周波数帯（例えば０．５ＨＺ〜２ｈｚ）程度を対象とすればよい。
【実施例】
【００３２】
次に、図１〜図６に基づいて、本発明に係る地震動波形推定方法の実施例について説明する。
図１に示すように、地震動予測フィルタを作成するための過去の地震として、２００４年紀伊半島南東沖地震を選択した。この地震では、前震、本震および余震が発生している。そこで、前震の観測記録から地震動予測フィルタを作成し、本震と余震の想定点における地震動を推定した。
【００３３】
先ず、防災科学技術研究所の基盤強震観測網（KiK-net）により収集された観測データの観測点から、Ａ、Ｂ、Ｃと直線状の測線を選定し、震源に最も近い観測点A1(SZOH36)、B1(AICH06)、C1(MIEH05)を上記地震発生予測地域とした。そして、本実施例においては、Ａ測線に沿って地震発生予測地域Ａ１で観測された地震動から、想定点Ａ３あるいは想定点Ａ４に伝搬する地震動について検証した。
【００３４】
図２は、上記地震発生予測地域Ａ１、想定点Ａ３および想定点Ａ４において観測された加速度の地震動波形を示すものであり、各々の地震動波形には、短周期成分が含まれている。
【００３５】
そこで次に、これらの観測された地震動波形に、フィルタをかけた。上述したように、当該フィルタは、地震動を推定する際にはオンラインで作用させる必要があるため、地震発生予測地域Ａ１の地震動波形に対しては、アナログフィルタやデジタルフィルタを用いる必要がある。ちなみに、本実施例においては、その１例として、４次のバンドパスフィルタを使用した。
【００３６】
このバンドパスフィルタの周波数特性は、下式（３）によって示される。
【００３７】
【数３】
2010175274.tif 000005
【００３８】
ここで、ζd＝0.707、ωd＝2.0、ζh＝0.6、ωh＝0.1とした。
図３は、式（３）によって示される上記４次のバンドパスフィルタの周波数特性を示すものであり、位相特性からも周波数成分に位相遅れが生じることが判る。
そして、このバンドパスフィルタに、地震発生予測地域Ａ１の地震動波形を通して、さらに１回積分することにより、上記対象となる周波数帯の速度波形を抽出した。
【００３９】
これに対して、想定点Ａ３、Ａ４で観測された地震動波形については、FFTを用いたバンドパスフィルタによって、対象となる周波数帯の位相遅れのない地震動波形を抽出し、さらに１回積分することにより速度波形を作成した。ただし、各々の地震動波形はＰ波到達時刻を０として、本実施例では２００秒間の波形を抽出して、それぞれに対して上記フィルタ処理を行った。
図４は、上記フィルタ処理によって抽出された地震発生予測地域Ａ１、想定点Ａ３および想定点Ａ４の上記対象となる周波数帯の地震動波形を示すものである。
【００４０】
次いで、上記フィルタ処理によって抽出された上記周波数帯の地震動波形を用いて、Ａ１−Ａ３間、あるいはＡ１−Ａ４間の伝達関数を同定した。
ここで、地震発生予測地域Ａ１と想定点Ａ３と間の距離は、５２ｋｍであり、地震発生予測地域Ａ１と想定点Ａ４との間の距離は、７８ｋｍである。そこで、Ａ１−Ａ３間の伝達関数を同定する場合には、想定点Ａ３で観測された地震動波形の開始時刻を１０秒間進めた。また、Ａ１−Ａ４間の伝達関数を同定する場合には、想定点Ａ４で観測された地震動波形の開始時間を１５．４秒進めた。これを換言すれば、想定点Ａ３に対して得られる地震動波形は、当該地震動が到達する１０秒前に、想定点Ａ４に対して得られる地震動波形は、当該地震動が到達する１５．４秒前に、それぞれ推定可能である事を示している。
【００４１】
また、伝達関数の同定には、式（１）、（２）で示した部分空間法を用いた。なお、伝達関数のモデル次数は６とした。
そして、得られた伝達関数に、地震発生予測地域Ａ１の地震動波形のフィルタ処理に用いた４次のバンドパスフィルタの周波数特性を掛け合わせて、地震動予測フィルタを作成した。
【００４２】
すなわち、式（３）で示した上記バンドパスフィルタの周波数特性は、状態方程式として、下記式（４）、（５）で表すことができる。
【００４３】
【数４】
2010175274.tif 000006
【００４４】
【数５】
2010175274.tif 000007
【００４５】
ここで、χ(i)は、地震動予測地域Ａ１で観測された（または観測される）地震動波形、u(i)は、上記バンドパスフィルタの出力、Z(i)は、上記バンドパスフィルタの状態量であり、ＡFは４×４、ＢFは４×１、ＣFは１×４、ＤFは４×１の行列である。
また、式（１）、（２）において、Ａ、Ｂ、Ｃ、Ｄの行列が同定された伝達関数の状態方程式は、式（６）、（７）に示すとおりである。
【００４６】
【数６】
2010175274.tif 000008
【００４７】
【数７】
2010175274.tif 000009
【００４８】
ここで、u(i)は、地震発生予測地域Ａ１で観測された地震動波形に対する上記バンドパスフィルタからの出力であって、上記伝達関数への入力である。また、y(i)は、上記伝達関数からの出力すなわち、想定点３（４）における地震動の推定波形である。
【００４９】
そして、上記バンドパスフィルタの周波数特性と上記伝達関数を掛け合わせて地震動予測フィルタを作成するとは、上記式（４）〜式（７）を併合して、下式（８）、（９）を作成することに相当する。
【００５０】
【数８】
2010175274.tif 000010
【００５１】
【数９】
2010175274.tif 000011
【００５２】
図５は、このようにして得られた地震動予測フィルタの周波数特性を示すものである。
そして、地震発生時に、上記式（８）、（９）に、地震発生予測地域で観測された地震動波形χ(i)を直接入力して、想定点Ａ３（Ａ４）の対象となる周波数帯の推定地震動波形を直接計算することができる。この際に、地震動波形χ(i)から推定地震動波形y(i)は、地震動データのサンプリング周期（0.01秒程度）１ステップ分の遅れで計算することができる。
【００５３】
図６は、上記地震動予測フィルタを用いて、想定点Ａ４における地震動を推定した結果を示すものである。なお、本実施例においては、前震の観測波形を用いて地震動予測フィルタを作成したため、前震での検証結果と、同じフィルタを用いて本震と余震の波形推定を行った結果を示している。また、本震と余震についても想定点Ａ４の観測波形は，前震と同じ１５．４秒進めている。
【００５４】
図６から、本発明に係る地震動波形推定方法によれば、地震発生時に、上記地震発生予測地域Ａ１で観測された地震動波形をリアルタイムで上記地震動予測フィルタに入力することにより、想定点Ａ４において対象となる周波数帯の地震動波形をリアルタイムで、かつ高い精度で推定することができることが検証された。
【００５５】
なお、上記実施形態および実施例においては、過去に観測された地震動波形を用いて地震動予測フィルタを作成する場合についてのみ説明したが、本発明はこれに限定されるものではなく、地震動シミュレーションによって得られた地震動波形を用いて上記地震動予測フィルタを作成することも可能である。
【産業上の利用可能性】
【００５６】
地震発生時に、その地震動が特定の想定点に到達する前に、当該想定点における地震動波形をリアルタイムに、かつ高い精度で推定するために利用される。
【符号の説明】
【００５７】
Ａ１  地震発生予測地域
Ａ３、Ａ４  想定点
【図１】
2010175274.tif 000012
【図２】
2010175274.tif 000013
【図３】
2010175274.tif 000014
【図４】
2010175274.tif 000015
【図５】
2010175274.tif 000016
【図６】
2010175274.tif 000017