JP2849297B2 - 震度予測システム - Google Patents
震度予測システムInfo
- Publication number
- JP2849297B2 JP2849297B2 JP358393A JP358393A JP2849297B2 JP 2849297 B2 JP2849297 B2 JP 2849297B2 JP 358393 A JP358393 A JP 358393A JP 358393 A JP358393 A JP 358393A JP 2849297 B2 JP2849297 B2 JP 2849297B2
- Authority
- JP
- Japan
- Prior art keywords
- point
- seismic intensity
- seismic
- observation
- prediction
- 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.)
- Expired - Lifetime
Links
- 230000033001 locomotion Effects 0.000 claims description 39
- 230000001133 acceleration Effects 0.000 claims description 15
- 230000003321 amplification Effects 0.000 claims description 9
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 9
- 238000012546 transfer Methods 0.000 claims description 9
- 238000004891 communication Methods 0.000 claims description 7
- 230000002238 attenuated effect Effects 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 9
- 238000011156 evaluation Methods 0.000 description 3
- 238000000034 method Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Landscapes
- Buildings Adapted To Withstand Abnormal External Influences (AREA)
- Geophysics And Detection Of Objects (AREA)
- Emergency Alarm Devices (AREA)
Description
【０００１】
【産業上の利用分野】この発明は震源から距離を隔てた
地点での震度を地震波が到達する以前に予測する、震度
予測システムに関するものである。
地点での震度を地震波が到達する以前に予測する、震度
予測システムに関するものである。
【０００２】
【従来技術及び発明が解決しようとする課題】地震動の
大きさは気象台や測候所等の特定の場所で体感され、発
生後に気象庁から発表される、揺れの程度を表す震度と
して評価されているが、実際には震源までの距離や地質
等によって揺れの程度が異なることから、発表される震
度の値には震源から距離を隔てた任意の地点での揺れの
大きさは反映されず、震度はあくまでも特定の地点で発
生した地震の結果を伝達する意味しか持たない。
大きさは気象台や測候所等の特定の場所で体感され、発
生後に気象庁から発表される、揺れの程度を表す震度と
して評価されているが、実際には震源までの距離や地質
等によって揺れの程度が異なることから、発表される震
度の値には震源から距離を隔てた任意の地点での揺れの
大きさは反映されず、震度はあくまでも特定の地点で発
生した地震の結果を伝達する意味しか持たない。
【０００３】地震の発生前にその時期や地点，規模を予
知することはあらゆる面での被害を最小にくい止めるた
めに重要な課題であるが、現状では各種の前兆現象を総
合しても地震学的な見地に基づく予知は難しく、その内
容は不確かである。
知することはあらゆる面での被害を最小にくい止めるた
めに重要な課題であるが、現状では各種の前兆現象を総
合しても地震学的な見地に基づく予知は難しく、その内
容は不確かである。
【０００４】多数の地域に分散して設置された地震計
（強震計）からは地震動の加速度，速度，変位の他、ス
ペクトル強度等の特性が直接得られ、震源までに距離が
ある場合にはこの特性値から地震波の到達以前に任意の
地点の震度を予測することが可能であるが、これらの物
理量は専門的であるため一般的な感覚から揺れの程度を
予想する材料にはなり難く、現在のところ、任意の地点
での震度を予測することは行われていない。
（強震計）からは地震動の加速度，速度，変位の他、ス
ペクトル強度等の特性が直接得られ、震源までに距離が
ある場合にはこの特性値から地震波の到達以前に任意の
地点の震度を予測することが可能であるが、これらの物
理量は専門的であるため一般的な感覚から揺れの程度を
予想する材料にはなり難く、現在のところ、任意の地点
での震度を予測することは行われていない。
【０００５】また地盤の揺れは地質や地形によって異な
り、建物の種別や階数等によっても揺れの大きさと振動
数が異なるため、震度予測の信頼度を高めるには認識が
容易で、異なる条件毎に個別に対応することが望まれ
る。
り、建物の種別や階数等によっても揺れの大きさと振動
数が異なるため、震度予測の信頼度を高めるには認識が
容易で、異なる条件毎に個別に対応することが望まれ
る。
【０００６】この発明はこうした背景を踏まえてなされ
たもので、任意の地点での震度を震源距離等に対応して
予測するシステムを新たに提案しようとするものであ
る。
たもので、任意の地点での震度を震源距離等に対応して
予測するシステムを新たに提案しようとするものであ
る。
【０００７】
【課題を解決するための手段】震度は加速度，または速
度と強い相関があることから、本発明では震源付近の観
測地点における地震波のＰ波到着からＳ波到着までの初
動部分の加速度波形，または速度波形から必要な情報を
算出し、この算出値から、震源からの距離に応じて揺れ
が減衰する式を用いて観測地点から隔てた、目的とする
予測地点での震度を評価することにより震源からの距離
に対応した適切な震度を地震波の到達前に予測し、その
予測地点での地震動の規模の認識を容易にする。
度と強い相関があることから、本発明では震源付近の観
測地点における地震波のＰ波到着からＳ波到着までの初
動部分の加速度波形，または速度波形から必要な情報を
算出し、この算出値から、震源からの距離に応じて揺れ
が減衰する式を用いて観測地点から隔てた、目的とする
予測地点での震度を評価することにより震源からの距離
に対応した適切な震度を地震波の到達前に予測し、その
予測地点での地震動の規模の認識を容易にする。
【０００８】初動部分の加速度波形，または速度波形か
らは観測地点の初動の震度と観測地点の震源距離，及び
目的とする予測地点の震源距離が算出され、これらの算
出値から距離減衰の回帰式を用いて予測地点での震度が
評価され、予測された震度は予測地点へ地震波が到達す
る前に伝達される。震度は人間が体感する揺れや震害の
程度を表す尺度として一般的に知られていることから、
容易に地震に関する情報として認識される。
らは観測地点の初動の震度と観測地点の震源距離，及び
目的とする予測地点の震源距離が算出され、これらの算
出値から距離減衰の回帰式を用いて予測地点での震度が
評価され、予測された震度は予測地点へ地震波が到達す
る前に伝達される。震度は人間が体感する揺れや震害の
程度を表す尺度として一般的に知られていることから、
容易に地震に関する情報として認識される。
【０００９】震源付近の観測地点で地震の初動を処理
し、必要な情報のみを算出して予測地点の震度を評価
し、これを予測地点へ伝達することにより地震波が予測
地点へ到達するまでの間の、短時間での予測が可能にな
る。地震波の到達以前に、10秒程度後に起こる地震の発
生が予告されることにより迅速に地震到来に備えた行動
に移ることが可能となり、また初動部分波形から得られ
た算出値から震源距離に応じて減衰する式を用いて予測
震度を評価することにより実測値に近い予測値が得ら
れ、実際の揺れに応じた事前の対応が可能になる。
し、必要な情報のみを算出して予測地点の震度を評価
し、これを予測地点へ伝達することにより地震波が予測
地点へ到達するまでの間の、短時間での予測が可能にな
る。地震波の到達以前に、10秒程度後に起こる地震の発
生が予告されることにより迅速に地震到来に備えた行動
に移ることが可能となり、また初動部分波形から得られ
た算出値から震源距離に応じて減衰する式を用いて予測
震度を評価することにより実測値に近い予測値が得ら
れ、実際の揺れに応じた事前の対応が可能になる。
【００１０】更に初動部分波形からの予測地点の震度の
予測に引き続き、観測地点へのＳ波到着後の主要動部分
波形から観測地点の震度を算出し、この値と既に求まっ
ている観測地点の震源距離，及び目的とする予測地点の
震源距離から距離減衰の回帰式を用いて地震波到達後の
予測地点の震度を予測することにより時々刻々と変化す
る震度をリアルタイムに更新し、より正確な予測を行
う。
予測に引き続き、観測地点へのＳ波到着後の主要動部分
波形から観測地点の震度を算出し、この値と既に求まっ
ている観測地点の震源距離，及び目的とする予測地点の
震源距離から距離減衰の回帰式を用いて地震波到達後の
予測地点の震度を予測することにより時々刻々と変化す
る震度をリアルタイムに更新し、より正確な予測を行
う。
【００１１】更にまた、震源付近の観測地点における地
震波のＳ波到着後の主要動部分の加速度波形及び速度波
形から地震動の周波数特性を算出する一方、目的とする
予測地点の、建物が建つ局所的な地盤地形の増幅特性を
作成し、両特性から建物に入力する地震動の周波数特性
を算出し、この周波数特性と建物の伝達関数を用いるこ
とにより個々の建物の、階数毎の震度を算出する。
震波のＳ波到着後の主要動部分の加速度波形及び速度波
形から地震動の周波数特性を算出する一方、目的とする
予測地点の、建物が建つ局所的な地盤地形の増幅特性を
作成し、両特性から建物に入力する地震動の周波数特性
を算出し、この周波数特性と建物の伝達関数を用いるこ
とにより個々の建物の、階数毎の震度を算出する。
【００１２】
【実施例】以下本発明を一実施例を示す図面に基づいて
説明する。
説明する。
【００１３】この発明は図１に示すように震源付近の観
測地点で地震波の初動部分波形を検出，処理し、この初
動部分波形から観測地点の初動の震度と観測地点の震源
距離，及び目的とする予測地点の震源距離を算出し、地
震波の予測地点への到達前にこれらの算出値から震源距
離に応じて揺れが減衰する式を用いて予測地点の震度を
予測するシステムである。
測地点で地震波の初動部分波形を検出，処理し、この初
動部分波形から観測地点の初動の震度と観測地点の震源
距離，及び目的とする予測地点の震源距離を算出し、地
震波の予測地点への到達前にこれらの算出値から震源距
離に応じて揺れが減衰する式を用いて予測地点の震度を
予測するシステムである。
【００１４】このシステムは図２に示すように観測地点
で地震波を処理する観測地点装置１と、観測地点装置１
で得られたデータから予測地点の震度を予測する中央観
測センター装置２と、中央観測センター装置２からの予
測震度を受信し、地震波の到来に備えた態勢を取る予測
地点装置３とから構成される。観測地点装置１は基本的
に地震計と地震波を処理するマイクロコンピュータ及び
中央観測センター装置２にデータを送信する送信装置と
から、中央観測センター装置２は受信装置とデータを処
理するホストコンピュータ及び送信装置とから、予測地
点装置３は受信装置からそれぞれ成り立っている。
で地震波を処理する観測地点装置１と、観測地点装置１
で得られたデータから予測地点の震度を予測する中央観
測センター装置２と、中央観測センター装置２からの予
測震度を受信し、地震波の到来に備えた態勢を取る予測
地点装置３とから構成される。観測地点装置１は基本的
に地震計と地震波を処理するマイクロコンピュータ及び
中央観測センター装置２にデータを送信する送信装置と
から、中央観測センター装置２は受信装置とデータを処
理するホストコンピュータ及び送信装置とから、予測地
点装置３は受信装置からそれぞれ成り立っている。
【００１５】観測地点装置１では地震計によって地震波
の水平２成分と上下成分の３成分の観測が行われ、観測
された波形は信号増幅器によって増幅された後、デジタ
ル・アナログ変換器によってデジタル信号に変換され
る。ここでマイクロコンピュータによって観測地点に最
初に到達する地震波のＰ波（初動部分）及び後続のＳ波
（主要動部分）が検知され、Ｐ波の加速度波形，または
速度波形から初動部分の振幅，位相，継続時間，水平動
と上下動の比等の特性が算出され、これらの特性値を用
いて震源の方向，震源深さ，観測地点での震度と、観測
地点の震源距離及び予測地点の震源距離が算定される。
の水平２成分と上下成分の３成分の観測が行われ、観測
された波形は信号増幅器によって増幅された後、デジタ
ル・アナログ変換器によってデジタル信号に変換され
る。ここでマイクロコンピュータによって観測地点に最
初に到達する地震波のＰ波（初動部分）及び後続のＳ波
（主要動部分）が検知され、Ｐ波の加速度波形，または
速度波形から初動部分の振幅，位相，継続時間，水平動
と上下動の比等の特性が算出され、これらの特性値を用
いて震源の方向，震源深さ，観測地点での震度と、観測
地点の震源距離及び予測地点の震源距離が算定される。
【００１６】マイクロコンピュータによって算出された
データは送信装置から通信回線を利用して中央観測セン
ター装置２へ送られる。
データは送信装置から通信回線を利用して中央観測セン
ター装置２へ送られる。
【００１７】中央観測センター装置２の受信装置で受信
されたデータはホストコンピュータに入力され、ホスト
コンピュータによって距離減衰の回帰式を用いて目的と
する予測地点の震度、更には建物の増幅特性を用いて建
物のある階の震度が評価される。
されたデータはホストコンピュータに入力され、ホスト
コンピュータによって距離減衰の回帰式を用いて目的と
する予測地点の震度、更には建物の増幅特性を用いて建
物のある階の震度が評価される。
【００１８】地震波の予測地点への到達前の震度の予測
には観測地点へのＰ波到着からＳ波到着までの初動部分
の加速度波形，または速度波形が用いられ、これから予
測地点の震度が算出される。Ｓ波の加速度波形，または
速度波形は予測地点への地震波到達後の、より正確な震
度の予測に利用される。
には観測地点へのＰ波到着からＳ波到着までの初動部分
の加速度波形，または速度波形が用いられ、これから予
測地点の震度が算出される。Ｓ波の加速度波形，または
速度波形は予測地点への地震波到達後の、より正確な震
度の予測に利用される。
【００１９】ここでは予測地点の震度が震源距離の対数
値に従って線形に減衰すると仮定して評価し、予測の内
容に応じて２段階に評価する方法を示す。
値に従って線形に減衰すると仮定して評価し、予測の内
容に応じて２段階に評価する方法を示す。
【００２０】予測地点の震度は以下の式から算出され
る。この式は予測地点で揺れが発生する10秒程度前に必
要な、緊急の予測に使用される。ここで用いられる観測
地点の震源距離と予測地点の震源距離は上記の通り、Ｐ
波の波形から観測地点装置１において既に求まってい
る。
る。この式は予測地点で揺れが発生する10秒程度前に必
要な、緊急の予測に使用される。ここで用いられる観測
地点の震源距離と予測地点の震源距離は上記の通り、Ｐ
波の波形から観測地点装置１において既に求まってい
る。
【００２１】予測震度＝c0＋（観測地点の初動の震度）
＋c1×log(観測地点の震源距離) −c2×log(予測地点の
震源距離) （c0，c1，c2は場所に応じた定数）上式の観
測地点の初動の震度は初動の波形を用い、２×log(初動
の最大振幅の絶対値) ＋c3より得られる。（c3は加速度
波形の場合は0.7 , 速度波形の場合は2.5 ）ここでは最
初に地震動を検知した地点、すなわち震源に最も近い観
測地点で観測された初動の波形が用いられる。
＋c1×log(観測地点の震源距離) −c2×log(予測地点の
震源距離) （c0，c1，c2は場所に応じた定数）上式の観
測地点の初動の震度は初動の波形を用い、２×log(初動
の最大振幅の絶対値) ＋c3より得られる。（c3は加速度
波形の場合は0.7 , 速度波形の場合は2.5 ）ここでは最
初に地震動を検知した地点、すなわち震源に最も近い観
測地点で観測された初動の波形が用いられる。
【００２２】より正確な予測を行う場合は上式による予
測地点における初動の震度の算出後、引き続いて以下の
式によりＳ波到着後の主要動波形を用い、時々刻々と変
化する震度を更新して評価する。
測地点における初動の震度の算出後、引き続いて以下の
式によりＳ波到着後の主要動波形を用い、時々刻々と変
化する震度を更新して評価する。
【００２３】予測震度＝d0＋（観測地点の震度）＋d1×
log(観測地点の震源距離) −d2×log(予測地点の震源距
離) （d0，d1，d2は場所に応じた定数）上式の観測地点
の震度は主要動波形を用い、２×log(初動の最大振幅の
絶対値) ＋c3より得られる。上記の各定数には観測記録
を分析した結果から得られた数値が用いられる。ここで
も最初に地震動を検知した観測地点で観測された主要動
の波形が用いられる。
log(観測地点の震源距離) −d2×log(予測地点の震源距
離) （d0，d1，d2は場所に応じた定数）上式の観測地点
の震度は主要動波形を用い、２×log(初動の最大振幅の
絶対値) ＋c3より得られる。上記の各定数には観測記録
を分析した結果から得られた数値が用いられる。ここで
も最初に地震動を検知した観測地点で観測された主要動
の波形が用いられる。
【００２４】建物のある階の震度は観測地点における地
震波のＳ波到着後の主要動部分の加速度波形及び速度波
形から得られる、地震波の周波数特性を表す卓越周波数
等のパラメータと、建物の増幅特性との関係から建物に
入力する地震動の周波数特性を評価した後、建物の伝達
関数のデータを加えて算出される。
震波のＳ波到着後の主要動部分の加速度波形及び速度波
形から得られる、地震波の周波数特性を表す卓越周波数
等のパラメータと、建物の増幅特性との関係から建物に
入力する地震動の周波数特性を評価した後、建物の伝達
関数のデータを加えて算出される。
【００２５】震源付近の観測地点で観測された、図３に
示すＳ波の主要動波形からはまず観測地点装置１のマイ
クロコンピュータによって地震動の周波数特性Ｓ0 が図
４に示すように分析，抽出され、この周波数特性Ｓ0 か
ら短時間で効率のよい評価とホストコンピュータへの伝
達を行うために、マイクロコンピュータでは周波数特性
Ｓ0 として卓越周波数Ｆd と周波数帯域幅Ｂの２つのパ
ラメータが算定され、共に通信回線を通じて中央観測セ
ンター装置２へ送信される。
示すＳ波の主要動波形からはまず観測地点装置１のマイ
クロコンピュータによって地震動の周波数特性Ｓ0 が図
４に示すように分析，抽出され、この周波数特性Ｓ0 か
ら短時間で効率のよい評価とホストコンピュータへの伝
達を行うために、マイクロコンピュータでは周波数特性
Ｓ0 として卓越周波数Ｆd と周波数帯域幅Ｂの２つのパ
ラメータが算定され、共に通信回線を通じて中央観測セ
ンター装置２へ送信される。
【００２６】観測地点の地震計では地震動の加速度と速
度の双方が計測され、図３に示すようにＳ波到着後の、
時間的余裕から判断される何秒間かの加速度波形と速度
波形を用い、卓越周波数Ｆd は最大加速度Ａmax と最大
速度Ｖmax の比からＦd ＝Ａmax ／Ｖmax ／（２π）と
して、周波数帯域幅Ｂは卓越周波数Ｆd と零交差率νd
（１秒当たりの零線を負から正方向に交差する回数）か
らＢ＝νd −Ｆd としてそれぞれ求められ、ホストコン
ピュータへ伝達される。
度の双方が計測され、図３に示すようにＳ波到着後の、
時間的余裕から判断される何秒間かの加速度波形と速度
波形を用い、卓越周波数Ｆd は最大加速度Ａmax と最大
速度Ｖmax の比からＦd ＝Ａmax ／Ｖmax ／（２π）と
して、周波数帯域幅Ｂは卓越周波数Ｆd と零交差率νd
（１秒当たりの零線を負から正方向に交差する回数）か
らＢ＝νd −Ｆd としてそれぞれ求められ、ホストコン
ピュータへ伝達される。
【００２７】ホストコンピュータでは図５に示すように
卓越周波数Ｆd と周波数帯域幅Ｂを用いて周波数特性Ｓ
0 が周波数特性Ｓ1 としてモデル化される一方、図６に
示すように予測地点の地盤や地形のデータに基き、建物
が建つ局所的な地盤地形の増幅特性Ｓ2 が作成され、図
７に示す建物の伝達関数Ｓ4 と共に前もってホストコン
ピュータに格納されている。
卓越周波数Ｆd と周波数帯域幅Ｂを用いて周波数特性Ｓ
0 が周波数特性Ｓ1 としてモデル化される一方、図６に
示すように予測地点の地盤や地形のデータに基き、建物
が建つ局所的な地盤地形の増幅特性Ｓ2 が作成され、図
７に示す建物の伝達関数Ｓ4 と共に前もってホストコン
ピュータに格納されている。
【００２８】建物に入力する地震動の周波数特性Ｓ3 は
モデル化された周波数特性Ｓ1 と増幅特性Ｓ2 を用い、
以下の式により評価される。
モデル化された周波数特性Ｓ1 と増幅特性Ｓ2 を用い、
以下の式により評価される。
【００２９】 Ｓi (f) ＝Ｓ1 (f) ×Ｓ2 (f) σ2 ＝２π∫Ｓi (f) df Ｓ3 (f) ＝Ｓi (f) ×10I-0.7 ／（γ2 σ2 ） ここでＩは前記の式で求められた建物地点での予測震
度、γはピークファクター、Ｓは全てパワースペクトル
で、振動数ｆの関数である。
度、γはピークファクター、Ｓは全てパワースペクトル
で、振動数ｆの関数である。
【００３０】この周波数特性Ｓ3 と前もって格納してお
いた建物の伝達関数Ｓ4 (f) のデータを用いて必要な階
の予測震度ＩS が以下の式により算定される。
いた建物の伝達関数Ｓ4 (f) のデータを用いて必要な階
の予測震度ＩS が以下の式により算定される。
【００３１】 ＳS (f) ＝Ｓ3 (f) ×Ｓ4 (f) ＩS ＝２log 〔γ×｛２π∫ＳS (f) df｝1/2 〕＋0.7 ∫はいずれも０から∞までの定積分であり、図８，図９
はそれぞれ建物への入力波の周波数特性Ｓ3 と建物の応
答周波数特性ＳS を示している。
はそれぞれ建物への入力波の周波数特性Ｓ3 と建物の応
答周波数特性ＳS を示している。
【００３２】上記の通り、中央観測センター装置２で算
出された予測震度は通信回線により予測地点装置３へ送
られ、予測地点装置３では予測地点における予測震度や
建物の階毎の予測震度に応じて警報を発する他、建物内
に設置された制震装置の抵抗力等を調整する等、数秒後
に起こる地震に備えた対応が取られる。地震の発生前に
警報を受けることにより建物内では火気類を消す等、危
険回避の行動に余裕を持って移行することが可能とな
る。
出された予測震度は通信回線により予測地点装置３へ送
られ、予測地点装置３では予測地点における予測震度や
建物の階毎の予測震度に応じて警報を発する他、建物内
に設置された制震装置の抵抗力等を調整する等、数秒後
に起こる地震に備えた対応が取られる。地震の発生前に
警報を受けることにより建物内では火気類を消す等、危
険回避の行動に余裕を持って移行することが可能とな
る。
【００３３】
【発明の効果】この発明は以上の通りであり、震源付近
の観測地点における地震波のＰ波到着からＳ波到着まで
の初動部分波形から必要な情報を算出し、この算出値か
ら、震源からの距離に応じて減衰する式を用いて観測地
点から隔てた、目的とする予測地点での震度を評価する
ものであるため震源からの距離に対応した適切な震度を
短時間で、地震波の到達前に予測することができ、また
揺れの程度を震度として表現するため容易に予測地点で
の地震動の規模の認識することが可能であり、事前に地
震到来に備えた態勢を取ることが可能である。特に初動
部分波形から得られた算出値から震源距離に応じて減衰
する式を用いて予測震度を評価するため実測値に近い予
測値が得られ、実際の揺れに応じた対応が可能になる。
の観測地点における地震波のＰ波到着からＳ波到着まで
の初動部分波形から必要な情報を算出し、この算出値か
ら、震源からの距離に応じて減衰する式を用いて観測地
点から隔てた、目的とする予測地点での震度を評価する
ものであるため震源からの距離に対応した適切な震度を
短時間で、地震波の到達前に予測することができ、また
揺れの程度を震度として表現するため容易に予測地点で
の地震動の規模の認識することが可能であり、事前に地
震到来に備えた態勢を取ることが可能である。特に初動
部分波形から得られた算出値から震源距離に応じて減衰
する式を用いて予測震度を評価するため実測値に近い予
測値が得られ、実際の揺れに応じた対応が可能になる。
【００３４】更に初動部分波形からの予測地点の震度の
予測に引き続き、観測地点へのＳ波到着後の主要動部分
波形から観測地点の初動の震度と観測地点の震源距離，
及び目的とする予測地点の震源距離を算出し、これらの
算出値から距離減衰の回帰式を用いて地震波到達後の予
測地点の震度を予測することにより時々刻々と変化する
震度を更新する。
予測に引き続き、観測地点へのＳ波到着後の主要動部分
波形から観測地点の初動の震度と観測地点の震源距離，
及び目的とする予測地点の震源距離を算出し、これらの
算出値から距離減衰の回帰式を用いて地震波到達後の予
測地点の震度を予測することにより時々刻々と変化する
震度を更新する。
【００３５】そして震源付近の観測地点における地震波
のＳ波到着後の主要動部分の加速度波形及び速度波形か
ら地震動の周波数特性を算出する一方、目的とする予測
地点の地盤や地形に基づく建物の増幅特性を作成してお
き、両特性から建物に入力する地震動の周波数特性を算
出し、この周波数特性と建物の伝達関数を用いることに
よりより正確な個々の建物の、階数毎の震度を予測する
ことができる。
のＳ波到着後の主要動部分の加速度波形及び速度波形か
ら地震動の周波数特性を算出する一方、目的とする予測
地点の地盤や地形に基づく建物の増幅特性を作成してお
き、両特性から建物に入力する地震動の周波数特性を算
出し、この周波数特性と建物の伝達関数を用いることに
よりより正確な個々の建物の、階数毎の震度を予測する
ことができる。
【図１】震度予測システムを示した概念図である。
【図２】震度予測システムの流れを示したブロック図で
ある。
ある。
【図３】観測地点での地震波を示した波形図である。
【図４】図３のＳ波の初動部分を抽出した周波数特性図
である。
である。
【図５】図４の波形をモデル化した周波数特性図であ
る。
る。
【図６】建物の増幅特性を示した周波数特性図である。
【図７】建物の伝達関数を示した周波数特性図である。
【図８】建物に入力する地震波を示した周波数特性図で
ある。
ある。
【図９】入力する地震波による建物の応答を示した周波
数特性図である。
数特性図である。
１……観測地点装置、２……中央観測センター装置、３
……予測地点装置。
……予測地点装置。
フロントページの続き (72)発明者 宮村 正光 東京都港区元赤坂１丁目２番７号 鹿島 建設株式会社内 (72)発明者 諸井 孝文 東京都港区元赤坂１丁目２番７号 鹿島 建設株式会社内 (72)発明者 山中 浩明 東京都港区元赤坂１丁目２番７号 鹿島 建設株式会社内 (72)発明者 大保 直人 東京都調布市飛田給２丁目19番１号 鹿 島建設株式会社技術研究所内 (72)発明者 高橋 克也 東京都調布市飛田給２丁目19番１号 鹿 島建設株式会社技術研究所内 (56)参考文献 特開 昭59−99277（ＪＰ，Ａ) 特開 昭62−268479（ＪＰ，Ａ) 特開 平３−115830（ＪＰ，Ａ) 実開 昭64−24268（ＪＰ，Ｕ) 日本建築設計管理協会連合会編著「建 築知識別冊・ＱＡシリーズ１新耐震設計 法」（株）建知出版発行，1982年10月10 日，21頁乃至31頁
Claims (1)
- 【請求項１】 震源付近の複数の観測地点に、地震計と
マイクロコンピュータと送信装置を備えた観測地点装置
を配設し、 中央観測地点に前記観測地点装置と通信回線で結ぶ、受
信装置とホストコンピュータと送信装置を備えた中央観
測センター装置を配設し、 震度予測地点に前記中央観測センター装置と通信回線で
結ぶ、受信装置を備えた予測地点装置よりなる震度予測
システムにおいて、 観測地点装置の地震計で地震を検知し、 その情報をマイクロコンピュータに入力し、地震波の初
動部分波形を検出し、 この初動部分波形から観測地点の初動の震度と、観測地
点の震源距離と、予測地点の震源距離を算出し、 この情報を送信装置によって、通信回線を介して中央観
測センター装置の受信装置が受信し、 中央観測センター装置に配設したホストコンピュータに
入力し、 地震波が予測地点へ到達する前に、震源距離に応じて揺
れが減衰する式を用いて、予測地点及び建物の周波数特
性を用いて建物の予測階の震度を予測するシステムであ
り、 観測地点装置のマイクロコンピュータで、観測地点へＳ
波到着後の主要動部分の加速度波形，または速度波形か
ら観測地点の震度と地震動の周波数特性を算出し、 前記観測地点へＳ波到着後の震度と、 前記観測地点の震源距離，及び目的とする予測地点の震
源距離から震源距離に応じて揺れが減衰する式を用い
て、中央観測センターのホストコンピュータで地震波到
達後の予測地点の震度を算出し、 予め目的とする予測地点の、建物が建つ局所的な地盤地
形の増幅特性を作成して、これをホストコンピュータに
格納しておき、前記マイクロコンピュータで得た周波数
特性とこの増幅特性から、建物に入力する地震動の周波
数特性を算出し、この周波数特性と建物の伝達関数と前
記ホストコンピュータで得た震度から建物の任意の階の
震度を算出し、 中央観測センター装置に配設した送信装置によって、通
信回線を介して予測した震度に関する情報を、予測地点
装置に配設した受信装置に送信するように構成したこと
を特徴とする震度予測システム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP358393A JP2849297B2 (ja) | 1993-01-12 | 1993-01-12 | 震度予測システム |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP358393A JP2849297B2 (ja) | 1993-01-12 | 1993-01-12 | 震度予測システム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH06324160A JPH06324160A (ja) | 1994-11-25 |
| JP2849297B2 true JP2849297B2 (ja) | 1999-01-20 |
Family
ID=11561479
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP358393A Expired - Lifetime JP2849297B2 (ja) | 1993-01-12 | 1993-01-12 | 震度予測システム |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2849297B2 (ja) |
Families Citing this family (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006045885A (ja) * | 2004-08-04 | 2006-02-16 | Taisei Corp | リアルタイム地震情報を利用した制振制御方法および制振制御システム |
| JP4506625B2 (ja) * | 2005-07-07 | 2010-07-21 | 大成建設株式会社 | リアルタイム地震情報を利用した地震動の予測システム |
| JP4491399B2 (ja) * | 2005-10-13 | 2010-06-30 | Ｏｋｉセミコンダクタ株式会社 | 地震防災システム |
| JP4771129B2 (ja) * | 2005-11-24 | 2011-09-14 | 清水建設株式会社 | 生産施設の地震リスク評価システム |
| JP5041727B2 (ja) * | 2006-04-28 | 2012-10-03 | 東芝エレベータ株式会社 | エレベータの地震監視制御装置及びその地震監視制御システム |
| JP2008249486A (ja) * | 2007-03-30 | 2008-10-16 | Railway Technical Res Inst | Ｓ波を用いたリアルタイム地震動分布推定方法 |
| JP5128973B2 (ja) * | 2007-03-30 | 2013-01-23 | 公益財団法人鉄道総合技術研究所 | 早期地震諸元推定情報の検証による警報の解除方法 |
| JP5064946B2 (ja) * | 2007-09-11 | 2012-10-31 | 東海旅客鉄道株式会社 | 警報用予測震度算出装置、地震警報報知システム |
| JP2013200284A (ja) * | 2012-03-26 | 2013-10-03 | Shimizu Corp | 地震検知装置 |
| JP5798545B2 (ja) * | 2012-11-02 | 2015-10-21 | 株式会社武田エンジニアリング・コンサルタント | 大地震と巨大地震の予知方法、予知装置、予知プログラム及び記録媒体 |
| TWI580991B (zh) * | 2013-02-25 | 2017-05-01 | Tokai Ryokaku Tetsudo Kk | Earthquake prediction device |
| JP6384093B2 (ja) * | 2014-04-01 | 2018-09-05 | 富士電機株式会社 | 地震計測装置 |
| JP6286264B2 (ja) * | 2014-04-08 | 2018-02-28 | 株式会社Ｎｔｔファシリティーズ | 構造物検証システム、構造物検証装置、構造物検証プログラム |
| CN109521221B (zh) * | 2018-12-10 | 2020-09-29 | 东北大学 | 一种钻爆法施工硬岩隧道微震波波速实时获取方法 |
| WO2022018789A1 (ja) * | 2020-07-20 | 2022-01-27 | 日本電信電話株式会社 | 予測方法、学習方法、予測装置、学習装置、予測プログラム及び学習プログラム |
| CN114200512A (zh) * | 2021-12-14 | 2022-03-18 | 四川旷谷信息工程有限公司 | 一种用于铁路施工期关键工点的地震烈度预警方法及系统 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0447498A (ja) * | 1990-06-13 | 1992-02-17 | Fujita Corp | 地震発生の早期警報装置 |
-
1993
- 1993-01-12 JP JP358393A patent/JP2849297B2/ja not_active Expired - Lifetime
Non-Patent Citations (1)
| Title |
|---|
| 日本建築設計管理協会連合会編著「建築知識別冊・ＱＡシリーズ１新耐震設計法」（株）建知出版発行，１９８２年１０月１０日，２１頁乃至３１頁 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH06324160A (ja) | 1994-11-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2849297B2 (ja) | 震度予測システム | |
| CN102262246B (zh) | 海上地震数据采集方法和装置 | |
| Li et al. | Seismic guided waves trapped in the fault zone of the Landers, California, earthquake of 1992 | |
| US5490062A (en) | Real-time neural network earthquake profile predictor | |
| Lam et al. | Generation of synthetic earthquake accelerograms using seismological modelling: a review | |
| US4562557A (en) | Method and apparatus for determining acoustic wave parameters from acoustic well logging waveforms | |
| KR102044041B1 (ko) | 지진측정장치 및 지진측정방법 | |
| KR101281803B1 (ko) | 지하 구조 영상화를 위해 라플라스 도메인에서의 파형 역산을 이용한 속도 분석 방법 | |
| KR102027252B1 (ko) | 임의 지진구역 대상 무작위 생성 지진학습데이터를 이용한 지진여부판별 방법 및 이에 의한 진도예측 시스템 | |
| KR102209749B1 (ko) | 지진가속도 시간이력 생성 방법 및 장치 | |
| US11789173B1 (en) | Real-time microseismic magnitude calculation method and device based on deep learning | |
| JP2007071707A (ja) | リアルタイム地震情報を利用した地震動強さの予測方法および防災システム | |
| US20190227184A1 (en) | Gauge length optimization for signal preservation and gauge length processing for distributed vibration sensing | |
| CN113189644B (zh) | 一种微震震源定位方法及系统 | |
| CN110164071B (zh) | 一种安防系统 | |
| KR101923166B1 (ko) | 지진배경잡음을 이용하여 추출된 지진신호의 진폭 크기 보정방법 | |
| Iezzi et al. | Seismo-acoustic characterization of Mount Cleveland volcano explosions | |
| Chandler et al. | An attenuation model for distant earthquakes | |
| RU2300122C1 (ru) | Способ дистанционного определения параметров инфразвукового сигнала вблизи неопознанного источника сигнала | |
| KR102196534B1 (ko) | Ｍｅｍｓ 네트워크를 이용한 실시간 지진동 감시 및 분포 지도 작성을 위한 장치 및 방법 | |
| Herrmann et al. | Surface wave inversion for shear wave velocity | |
| US7440854B2 (en) | Density and velocity based assessment method and apparatus | |
| CN107167230B (zh) | 利用建筑结构振动衰减率评估建筑物二次辐射噪声的方法 | |
| CN104570145B (zh) | 一种光纤传感生命探测设备 | |
| CN114200512A (zh) | 一种用于铁路施工期关键工点的地震烈度预警方法及系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 19970114 |