JP5313706B2 - 地震動波形推定方法 - Google Patents
地震動波形推定方法 Download PDFInfo
- Publication number
- JP5313706B2 JP5313706B2 JP2009015348A JP2009015348A JP5313706B2 JP 5313706 B2 JP5313706 B2 JP 5313706B2 JP 2009015348 A JP2009015348 A JP 2009015348A JP 2009015348 A JP2009015348 A JP 2009015348A JP 5313706 B2 JP5313706 B2 JP 5313706B2
- Authority
- JP
- Japan
- Prior art keywords
- earthquake
- waveform
- filter
- seismic
- motion
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 230000005540 biological transmission Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 5
- 230000035939 shock Effects 0.000 description 4
- 238000012795 verification Methods 0.000 description 4
- 230000002265 prevention Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Description
そして、震源から数Km以上離れた地点においては、S波に起因する地震動の主要動が到達するまでに数秒から数十秒の余裕があることから、上記リアルタイム地震情報を利用することにより、地震発生後数秒程度で、発生した地震のマグニチュードや震源位置等に関する情報を受け取って、地震被害の発生を防止することが可能になっている。
また、請求項2に記載の発明は、請求項1に記載の発明において、上記第1のフィルタが、観測された上記地震動波形からリアルタイムで周波数帯の地震動波形を抽出可能なデジタルフィルタまたはアナログフィルタであり、かつ上記第2のフィルタが、高速フーリエ変換機能(Fast Fourier Transform:以下FFT)を用いたバンドパスフィルタであることを特徴とするものである。
そこで、新たな地震が発生した際に、上記地震発生予測地域で観測された地震動波形をリアルタイムに上記地震動予測フィルタに入力することにより、想定点において対象となる周波数帯の地震動波形をリアルタイムで、かつ高い精度で推定することができる。
先ず、今後地震が起きると推定される震源地またはその近傍の地震発生予測地域と、当該地震発生予測地域から離れた特定の想定点を特定する。そして、上記震源地において過去に発生した地震時に観測された上記地震発生予測地域および想定点での地震動波形に基づいて、上記地震発生予測地域と想定点との間の伝達関数を同定する。
例えば、超高層ビルや免震構造物の制御に利用する場合は、これらの構造物の卓越周波数が存在する0.5Hz〜0.16Hz(2秒〜6秒)程度の範囲を含む周波数帯を対象とすればよい。また、1Hzの卓越周波数を持つ機器の場合は、1Hzを含む前後の周波数帯(例えば0.5HZ〜2hz)程度を対象とすればよい。
図1に示すように、地震動予測フィルタを作成するための過去の地震として、2004年紀伊半島南東沖地震を選択した。この地震では、前震、本震および余震が発生している。そこで、前震の観測記録から地震動予測フィルタを作成し、本震と余震の想定点における地震動を推定した。
図3は、式(3)によって示される上記4次のバンドパスフィルタの周波数特性を示すものであり、位相特性からも周波数成分に位相遅れが生じることが判る。
そして、このバンドパスフィルタに、地震発生予測地域A1の地震動波形を通して、さらに1回積分することにより、上記対象となる周波数帯の速度波形を抽出した。
図4は、上記フィルタ処理によって抽出された地震発生予測地域A1、想定点A3および想定点A4の上記対象となる周波数帯の地震動波形を示すものである。
ここで、地震発生予測地域A1と想定点A3と間の距離は、52kmであり、地震発生予測地域A1と想定点A4との間の距離は、78kmである。そこで、A1−A3間の伝達関数を同定する場合には、想定点A3で観測された地震動波形の開始時刻を10秒間進めた。また、A1−A4間の伝達関数を同定する場合には、想定点A4で観測された地震動波形の開始時間を15.4秒進めた。これを換言すれば、想定点A3に対して得られる地震動波形は、当該地震動が到達する10秒前に、想定点A4に対して得られる地震動波形は、当該地震動が到達する15.4秒前に、それぞれ推定可能である事を示している。
そして、得られた伝達関数に、地震発生予測地域A1の地震動波形のフィルタ処理に用いた4次のバンドパスフィルタの周波数特性を掛け合わせて、地震動予測フィルタを作成した。
また、式(1)、(2)において、A、B、C、Dの行列が同定された伝達関数の状態方程式は、式(6)、(7)に示すとおりである。
そして、地震発生時に、上記式(8)、(9)に、地震発生予測地域で観測された地震動波形χ(i)を直接入力して、想定点A3(A4)の対象となる周波数帯の推定地震動波形を直接計算することができる。この際に、地震動波形χ(i)から推定地震動波形y(i)は、地震動データのサンプリング周期(0.01秒程度)1ステップ分の遅れで計算することができる。
A3、A4 想定点
Claims (3)
- 過去の地震時に地震発生予測地域で観測された地震動波形と上記地震発生予測地域から離れた特定の想定点において観測された地震動波形とに基づいて、新たな地震発生時に上記想定点における地震動波形をリアルタイムに推定する地震動波形推定方法において、
上記過去の地震時における上記地震動波形を得るに際して、上記地震発生予測地域で観測された上記地震動波形から、リアルタイムで第1のフィルタを介して上記想定点において上記新たな地震発生時に到達する地震動の大きさを判断するための対象となる周波数帯の地震動波形を抽出し、上記想定点で観測された上記地震動波形から当該地震動波形の形状および位相にずれを生じない第2のフィルタを介して上記周波数帯の地震動波形を抽出し、これら地震発生予測地域および想定点で得られた上記周波数帯の地震動波形から上記地震発生予測地域と上記想定点との間の伝達関数を求めた後に、上記第1のフィルタの周波数特性に当該伝達関数を掛け合わせた周波数伝達特性を有する地震動予測フィルタを作成し、
上記新たな地震発生時に、上記地震発生予測地域で観測された上記地震動波形をリアルタイムに上記地震動予測フィルタに入力することにより、上記想定点における上記周波数帯の地震動波形をリアルタイムで推定することを特徴とする地震動波形推定方法。 - 上記第1のフィルタは、観測された上記地震動波形からリアルタイムで周波数帯の地震動波形を抽出可能なデジタルフィルタまたはアナログフィルタであり、かつ上記第2のフィルタは、高速フーリエ変換機能を用いたバンドパスフィルタであることを特徴とする請求項1に記載の地震動波形推定方法。
- 予め、一箇所の上記想定点と複数の上記地震発生予測地域との間に対して、各々対応する複数の上記地震動予測フィルタを作成し、上記新たな地震発生時に、最も近接した上記地震発生予測地域と上記想定点との間の上記地震動予測フィルタを用いことを特徴とする請求項1または2に記載の地震動波形推定方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009015348A JP5313706B2 (ja) | 2009-01-27 | 2009-01-27 | 地震動波形推定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009015348A JP5313706B2 (ja) | 2009-01-27 | 2009-01-27 | 地震動波形推定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2010175274A JP2010175274A (ja) | 2010-08-12 |
JP5313706B2 true JP5313706B2 (ja) | 2013-10-09 |
Family
ID=42706405
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2009015348A Active JP5313706B2 (ja) | 2009-01-27 | 2009-01-27 | 地震動波形推定方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5313706B2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105676286A (zh) * | 2016-01-25 | 2016-06-15 | 西南交通大学 | 一种用于地震预警系统的实时地震震级估计方法 |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5806031B2 (ja) * | 2011-07-29 | 2015-11-10 | 株式会社ミツトヨ | 長周期振動感知器、及び長周期振動感知器の出力値補正方法 |
JP6177628B2 (ja) * | 2013-08-29 | 2017-08-09 | 公益財団法人鉄道総合技術研究所 | 地中地震動を利用した即時地震動予測方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4385948B2 (ja) * | 2005-01-11 | 2009-12-16 | 大成建設株式会社 | リアルタイム地震情報を利用したリアルタイム地震応答波形推定方法 |
JP4758759B2 (ja) * | 2005-12-27 | 2011-08-31 | 株式会社ミツトヨ | 地震災害防止システム |
JP4229337B2 (ja) * | 2006-09-28 | 2009-02-25 | 独立行政法人防災科学技術研究所 | 計測震度概算装置、それを用いた計測震度概算システム及び計測震度概算方法 |
-
2009
- 2009-01-27 JP JP2009015348A patent/JP5313706B2/ja active Active
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105676286A (zh) * | 2016-01-25 | 2016-06-15 | 西南交通大学 | 一种用于地震预警系统的实时地震震级估计方法 |
CN105676286B (zh) * | 2016-01-25 | 2017-11-14 | 西南交通大学 | 一种用于地震预警系统的实时地震震级估计方法 |
Also Published As
Publication number | Publication date |
---|---|
JP2010175274A (ja) | 2010-08-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Rainieri et al. | Automated output-only dynamic identification of civil engineering structures | |
A Hasan et al. | Enhanced frequency domain decomposition algorithm: a review of a recent development for unbiased damping ratio estimates | |
Rainieri et al. | Robust output-only modal identification and monitoring of buildings in the presence of dynamic interactions for rapid post-earthquake emergency management | |
US20150330950A1 (en) | Structural fatigue crack monitoring system and method | |
JP4745170B2 (ja) | 漏水検出装置および漏水検出方法 | |
KR102209749B1 (ko) | 지진가속도 시간이력 생성 방법 및 장치 | |
JP6429042B2 (ja) | 地盤特性を自動校正する現地型地震早期警報システム及び関連方法 | |
KR101086448B1 (ko) | 지진파형의 시간영역 부지응답특성 변환을 통한 인접부지의 지진동세기추정시스템 및 그 추정방법 | |
JP5313706B2 (ja) | 地震動波形推定方法 | |
KR101724152B1 (ko) | 금속 이물질의 충격 위치 추정을 위한 분산 선도 추출 방법 | |
TWI809293B (zh) | 地震觀測裝置、地震觀測方法以及記錄地震觀測程式的記錄媒體 | |
JP4385948B2 (ja) | リアルタイム地震情報を利用したリアルタイム地震応答波形推定方法 | |
Chen et al. | An empirical evolutionary magnitude estimation for early warning of earthquakes | |
Paul et al. | Phase space interrogation of the empirical response modes for seismically excited structures | |
Liu et al. | Lower-order modal parameters identification for offshore jacket platform using reconstructed responses to a sea test | |
KR101923166B1 (ko) | 지진배경잡음을 이용하여 추출된 지진신호의 진폭 크기 보정방법 | |
JP2006275696A (ja) | 早期地震諸元推定方法及びそのシステム | |
JP5770051B2 (ja) | 大地震の破壊領域の即時的推定方法 | |
JP2018200313A (ja) | 建物への入力地震動の評価方法 | |
JP5042872B2 (ja) | 高架橋のリアルタイム地震被害推定方法及びその装置 | |
JP6401003B2 (ja) | 単独観測点の地震波を用いた近距離地震に対する早期警報方法 | |
JP6559579B2 (ja) | 観測地震の距離パラメータを推定する方法ならびに距離パラメータ推定プログラムおよび距離パラメータ推定プログラムを記録したコンピュータ読み取り可能な記録媒体 | |
JP5760268B2 (ja) | 地震動継続時間予測システム | |
JP2009186385A (ja) | 高架橋の揺れ方によるリアルタイム地震被害推定方法及びその装置 | |
JP6559578B2 (ja) | 地震の主要動の最大振幅を推定する方法ならびに最大振幅推定プログラムおよび最大振幅推定プログラムを記録したコンピュータ読み取り可能な記録媒体 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20111018 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20121120 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20130117 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20130219 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20130702 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20130704 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 Ref document number: 5313706 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |