JPH0263194B2

JPH0263194B2 -

Info

Publication number: JPH0263194B2
Application number: JP60154131A
Authority: JP
Inventors: Yutaka Nakamura; Kenji Tomita; Tamio Hashimoto
Original assignee: Railway Technical Research Institute; Oki Electric Industry Co Ltd
Current assignee: Railway Technical Research Institute; Oki Electric Industry Co Ltd
Priority date: 1985-07-15
Filing date: 1985-07-15
Publication date: 1990-12-27
Also published as: JPS6215482A

Description

【発明の詳細な説明】（技術分野）本発明は、地震波の観測点から震源までの距離
（震源距離）をリアルタイムで時々刻々と計算す
る震源距離算出装置である。

（従来技術）従来、震源距離を算出する方法として、３点以
上の複数の観測点で地震波の到来時刻を検出し
て、震源を求めてから震源距離を算出する方法、
１観測点でＰ波とＳ波の到来時刻の差（Ｐ〜Ｓ時
間）を検出し、Ｐ〜Ｓ時間が震源距離に比例する
ことを利用して、震源距離を算出する方法があ
る。いずれも、地震波の到来から数10秒ないし数
分経過しないと、震源距離を知ることができな
い。

（発明の目的）本発明の目的は、固有周期が異なるようにして
得た観測出力から、地震の震源距離を算出する装
置を提供することである。

（発明の構成）本発明は、固有周期の異なる２種の地震計によ
る地震観測出力を１組とし、固有周期が最も長い
地震計が含まれる組と固有周期が最も短い地震計
が含まれる組の２組の観測出力を入力し、各組の
観測出力の指数平滑値を算出する手段と、各組の
地震計の観測波形の前記指数平滑値と各組の地震
計の固有周期とにより、高低２つの周波数帯域に
おける卓越周波数とその成分振幅をそれぞれ算出
する卓越周波数算出手段と、前記２組の卓越周波
数算出手段より得た２つの成分振幅の比の対数と
前記２組の卓越周波数算出手段より得た２つの卓
越周波数の差との比を算出する手段とを有して観
測地点の震源距離を得ることを特徴とする震源距
離算出装置である。

（第１の実施例）震源付近からはあらゆる周波数成分の波動が均
一に周辺に放射していると考えることができる。
しかし、地震波動が伝播するとき、一般に高い周
波数の振動ほど減衰しやすいため、地震波動中の
高い周波数成分の振幅と低い周波数成分の振幅と
の比は波動の伝播距離が大きくなるほど小さくな
る。したがつて、ある地点で観測される地震波動
中の二つの卓越振動数の振動の振動数と成分振幅
値の組合せにより、その地点が震源からどのくら
い離れているか推定することができる。また、ひ
とつの地震を固有周期の異なるいくつかの地震計
で同時に観測すると、それぞれの地震計はそれぞ
れの固有周期に近い振動に強く感応するため、そ
れらの出力波形は第１図のように異なつたものと
なるが、地震計の固有周期があらかじめわかつて
いれば、地震動の卓越周期や振幅をこれらの地震
計の出力振幅から推定することができる。本発明
は、固有周期の異なる２種の地震計による地震波
観測出力を１組として、２つの振動数帯域に属す
る２組の地震波観測出力に基づいて、それぞれの
帯域での卓越振動数と振幅を算出し、これらを用
いて震源距離を推定しようとするものである。

第２図は第１の実施例を示すブロツク図であ
る。１，２，３，４は同一地点の地動を検出する
センサでそれぞれの固有周波数は、たとえば第１
のセンサ１は0.1Hz、第２のセンサ２は0.5Hz、第
３のセンサ３は１Hz、第４のセンサ４は５Hzのよ
うに異なる。５１，５２，５３，５４は前記セン
サ１，２，３，４よりの検出信号を増幅するアン
プ、６１，６２，６３，６４はバツフアアンプ、
８はセンサ１，２，３，４からの検出信号より地
動のデータを処理して、震源距離を算出する処理
装置である。処理装置８は標本化手段１１１，１
１２，１１３，１１４、指数平滑化手段１２１，
１２２，１２３，１２４、卓越周期算出手段１３
１，１３２、震源距離算出手段１４、震源距離出
力手段１５により構成される。第３図は実施例の
フローチヤートである。

センサ１〜４は設置されている地点の地動を常
時検出し、これを電気的信号に変換して処理装置
８に送出している。処理装置８は、センサ１，
２，３，４で検出されてアンプ５１，５２，５
３，５４、バツフアアンプ６１，６２，６３，６
４を介して絶えず送られてくる地動の情報を取り
こんでいる。

標本化手段１１は、前記各センサ１，２，３，
４の情報を所定の時間間隔（たとえば1/50〜1/15
０）で標本化し、その過去何回分かのサンプリン
グ情報の平均値より入力データの直流分であるオ
フセツトレベルの算出を行う。このオフセツトレ
ベルは時々刻々得られるサンプリング情報によつ
て絶えず更新されている。また、センサ１，２，
３，４によつて検出されて時々刻々送られてくる
４つの地動情報のサンプリング値からオフセツト
レベルを除去した値x₁（ｔ）、x₂（ｔ）、x₃（ｔ）、x₄
（ｔ）を算出する。サンプリング値をx_si（ｔ）、オ
フセツトレベルをx_piとすると、x_i（ｔ）＝x_si（ｔ）
−x_pi（ｉ＝１、２、３、４）である。

指数平滑化手段１２１，１２２，１２３，１２
４は前記のx_i（ｔ）の平滑化を行う。すなわち、
指数平滑値x_ai（ｔ）は次式で算出する。

x_ai（ｔ）＝x_ai（ｔ−１）×α_i＋x_i ²（ｔ） ……(1) α_iは0.9程度の定数である。またx_ai（ｔ−１）は
時刻ｔに対し１サンプル前の指数平滑値である。

本平滑法により、信号に重畳したノイズ成分の
影響を低く抑えて、振動波形の包絡線の変動をリ
アルタイムに把握することができる。

因みにこの(1)式の意味は、以下のことを示すこ
とは自明である。

x_ai（ｔ）＝x_ai（ｔ−１）×α_i＋x_i ²（ｔ）＝（x_ai（
ｔ−２）×α_i＋x_i ²（ｔ−１））×α_i＋x₁ ²（ｔ）＝x_ai（ｔ−２）×α_i ²＋x_i ²（ｔ−１）×α_i＋x₁（
ｔ）＝（x_ai（ｔ−３）×α_i＋x_i ²（ｔ−２））×α_i ²
＋x_i ²（ｔ−１）×
α_i＋x_i ²（ｔ）＝x_ai（ｔ−３）×α_i ³＋x_i ²（ｔ−２）×α_i ²＋x_i ²
（ｔ−１）×α_i＋x_i ²（ｔ）＝…＋x_i ²（ｔ−ｎ）×α_i ⁿ＋…＋x_i ²（ｔ−３） ×α_i ³＋x_i ²（ｔ−２）×α_i ²＋x_i ²（ｔ−１）×α
_i＋x_i ²（ｔ）すなわち、(1)式は、現在の瞬時自乗振幅の影響
はそのままで、１ステツプ過去に遡る毎にその時
の瞬時自乗振幅の影響度がα_i倍になるような瞬時
自乗振幅の累積を表わしている。もし、x_i ²（ｔ）
が一定であれば、(1)式はこの一定値に初期値１公
比α_iの等比級の和を乗じたものになる。もし、α_i
が１未満であればこの等比級数の和は有限とな
り、１／（１−α_i）、となる。したがつて、(1)式
に（１−α_i）を乗じたものは自乗振幅値を平滑化
した包絡振幅をあらわす。このような平滑は過去
の瞬時振幅の影響が指数関数的に減少するため、
指数平滑と呼んでいる。このような平滑法は電気
の平滑回路をデイジタル的に模したものであり、
電気では平滑回路は一般的に用いられている。

本平滑法により、信号に重畳したノイズ成分の
影響を低く抑えて、振動波形の包絡線の変動をリ
アルタイムに把握することができる。なお、(2)式
などにみられるx_iなどの振幅値は、式の展開から
も明らかなように瞬時瞬間の振幅値ではなく振動
波形の包絡線の変動で表わされる振幅値である。
また、それらの比を算出するものであるから包絡
振幅値が定数倍されていても問題はない。そこ
で、計算を簡単にするためにも包絡振幅値を算出
する平滑法を定義するにあたつて、（１−α_i）を
乗じなかつた。

１台のセンサの固有周期₁Hz、出力信号をx₁
（ｔ）、他の１台の固有周波数₂Hz、出力信号をx₂
（ｔ）とすると、卓越周波数（ｔ）は、次式で計
算する。

（ｔ）² ＝｛（₂ ²×x₂（ｔ）−₁ ²×x₁（ｔ）／₂ ²×x₁ ²
（ｔ）−₁ ²×x₂ ²（ｔ）｝^1/2・₁ ₂ センサ１，２の固有周波数をω₁、ω₂、減衰定
数をh₁、h₂とし、地表振動の真の加速度波形をｙ
とすると(2)、(3)式で示す運動方式に従つて出力波
形が電気信号として出力される。

x¨₁＋2h₁ω₁x〓₁＋ω₁ ²x₁＝2h₁ω₁y¨ (2) x¨₂＋2h₂ω₂x〓₂＋ω₂ ²x₂＝2h₂ω₂y¨ (3) この場合、センサ１，２は振子の振動がそのま
ま電気信号に変えられて出力するものであつて、
x₁、x₂が相対変位に、x〓₁、x〓₂が相対速度に、x¨₁、
x¨₂が相対加速度に対応した出力波形である。

センサ１，２の出力値は種々のものをとること
が可能であるが、本実施例では変位対応のx₁、x₂
を出力するものを例示している。

ｙ＝ae^j〓^t（ａ＞０） (4) x₁＝x₁e^j(〓^t+〓¹⁾（x₁＞０） (5) x₂＝x₂e^j(〓^t+〓²⁾（x₂＞０） (6) とおくと、(2)式より（ω₁ ²−ω²＋j2h₁ω₁ω）x₁e^j(〓^t+〓¹⁾ ＝−2h₁ω₁aω²e^j〓^t ∴x₁＝−2h₁ω₁aω²／ω₁−ω²＋j2h₁ω₁ω・e^-j〓¹(7) x₁は速度に比例する出力であり、実数（正）で
あるから、同様にして、入力速度振幅はa〓、これを未知数Ｘとおく、
また入力振動数はω、これも未知数と考え、Ｘ、
ωをx₁、x₂から導く。

(8)式、(9)式より、 x₁ ²｛（ω₁ ²−ω²）²＋（2h₁ω₁ω）²｝＝4h₁ ²ω₁ ²ω²X²(10) x₂ ²｛（ω₂ ²−ω²）²＋（2h₁ω₁ω）²｝＝4h₂ ²ω₂ ²ω²X²(11) さらに(10)、〓式より 4h₂ ²ω₂ ²x₁ ²｛（ω₁ ²−ω²）²＋（2h₁ω₁ω）²｝−4h₁ ²ω₁ ²x₂ ²｛（ω₂ ²−ω²）²＋（2h₂ω₂ω）²｝＝０（h₂ ²ω₂ ²x₁ ²−h₁ ²ω₁ ²x₂ ²）ω⁴＋（−2h₂ ²x₁ ²＋4h₁ ²h₂ ²x₁ ² ＋2h₁ ²x₂ ²−4h₂ ²h₁ ²x₂ ²）ω₁ ²ω₂ ²ω²＋h₂ ²ω₂ ²x₁ ²ω
₁ ⁴−h₁ ²ω₁ ²x₂ ²ω₂ ⁴＝０（h₂ ²ω₂ ²x₁ ²−h₁ ²ω₁ ²x₂ ²）ω⁴＋｛−2h₂ ²x₁（１−2h₁ ²）＋2h₁ ²x₂ ²（１ −2h₂ ²）｝ω₁ ²ω₂ ²ω²＋h₂ ²ω₂ ²x₁ ²−h₁ ²ωx₂ ²ω₂ ⁴
＝０(12) ここで、Ａ＝h₂ ²ω₂ ²x₁ ²−h₁ ²ω₁ ²x₂ ² （13）Ｂ＝｛−h₂ ²x₁ ²（１−2h₁ ²）＋h₁ ²x₂ ²（１−2h₂ ²）
｝ω₁ ²ω₂ ²（14）Ｃ＝h₂ ²ω₂ ²x₁ ²ω₁ ⁴−h₁ ²ω₁ ²x₂ ²ω₂ ⁴ （15）とおけば、(12)式は Aω⁴＋2Bω²＋Ｃ＝０（16）となる。（16）式をω²についてω²＞０を考慮して
解けばとなる。ω²＞０であるためには B²−AC＞０かつAC＜０。

(10)、(11)式より x₁ ²＝（2h₁ω₁ωX）²／（ω₁ ²−ω²）²＋（2h₁ω₁ω
）²（18） x₂ ²＝（2h₂ω₂ωX）²／（ω₂ ²−ω²）²＋（2h₂ω₂ω
）²（19）（18）、（19）式を（13）式に代入してＡ＝4h₁ ²h₂ ²ω₁ ²ω₂ ²ω²X²／（ω₁ ²−ω²）²＋（2h₁ω₁ω）²−4h₁ ²h₂ ²ω₁ ²ω₂ ²ω²X²／（ω₂ ²−ω²）²＋（2h₂ω₂ω）² ＝4h₁ ²h₂ ²ω₁ ²ω₂ ²ωX²／｛（ω₁ ²−ω²）²＋（2h₁ω
₁ω）²｝｛（ω₂ ²−ω）²｝＋（2h₂ω₂ω）²｝・｛（ω₂ ²＋ω₁ ²）（ω₂ ²−ω₁ ²）−２（ω₂ ²−ω₁ ²）ω²＋４（h₂ ²ω₂ ²−h₁ ²ω₁ ²）ω²｝（20）（18）、（19）式を（15）式に代入してＣ＝h₂ ²ω₂ ²・（2h₁ω₁ωX）²／（ω₁ ²−ω²）²＋（2h₁ω₁ω）²・ω₁ ⁴−h₁ ²ω₁ ²・（2h₂ω₂ωX）²／（ω₂ ²−
ω²）²＋（2h₂ω₂ω）²・ω₂ ⁴ ＝4h₁h₂ ²ω₁ ²ω₂ ²ω²X²／｛（ω₁ ²−ω²）²＋（2h₁ω
₁ω）²｝｛（ω₂ ²−ω²）²＋（2h₂ω₂ω）²｝・〔ω₁ ⁴｛（ω₂ ²−ω²）²＋（2h₂ω₂ω）²−ω₂ ⁴｛
（ω₁ ²−ω²）²＋（2h₁ω₁ω）²｝〕（21）（ω₁ ²−ω²）²＋（2h₁ω₁ω）²＞０、（ω₂ ²−ω²）²＋２（h₂ω₂ω）²＞０ ω₂ ²＞ω₁ ²であるので、（20）式×（21）式＝AC
＜０であるためには、｛（ω₂ ²＋ω₁ ²）（ω₂ ²−ω₁ ²） −２（ω₂ ²−ω₁ ²）ω² ＋４（h₂ ²ω₂ ²−h₁ ²ω₁ ²）ω²｝・〔ω₁ ⁴｛（ω₂ ²−ω²）²＋（2h₂ω₂ω）²｝ −ω₂ ⁴｛（ω₁ ²−ω²）²＋（2h₁ω₁ω）²｝〕＜０であればよい。上式左辺を整理すると、 ω²〔（ω₂ ²＋ω₁ ²）（ω₂ ²−ω₁ ²）＋2ω²｛（１−2h₁ ²）ω₁ ²−（１−2h₂ ²）ω₂ ²）｝〕・〔ω²（ω₁ ² ＋ω₂ ²）（ω₁ ²−ω₂ ²）−2ω₁ ²ω₂ ²｛（１−2h₂ ²）
ω₁ ²−（１−2h₁ ²）ω₂ ²｝〕＜０（22）したがつて（１−2h₁ ²）ω₁ ²−（１−2h₂ ²）ω₂ ²０（１−2h₂ ²）ω₁ ²−（１−2₁ ²）ω₂ ²０であれば、（22）式を満足する。

h₁＝h₂＝１／√２に選べばＡ＝１／２（ω₂ ²x₁ ²−ω₁ ²x₂ ²）Ｂ＝０Ｃ＝１／２ω₁ ²ω₂ ²（x₁ ²ω₁ ²−x₂ ²ω₂ ²）よつて（14）式はとなり、この（23）式が卓越周期を与える。すな
わち ω＝｛（x₂ ²ω₂ ²−x₁ ²ω₁ ²）／（x₁ ²ω₂ ²−x₂ ²ω₁ ²）
｝^1/4（ω₁ω₂）^1/2（24）＝｛（x₂ ² ₂ ²−x₁ ² ₁ ²）／（x₁ ² ₂ ²−x₂ ² ₁ ²）
｝^1/4（₁ ₂）^1/2 一方、Ｘについては、(10)式より X²＝x₁ ²｛（ω₁ ²−ω²）²＋（2h₁ω₁ω）²｝／（4h₁ ²ω₁ ²ω²）（25）この式にh₁＝１／√２を代入して X²＝x₁（ω₁ ⁴＋ω⁴）／2ω₁ ²ω² （26）さらに、（17）式のω²を代入すると従つて第１の卓越周期算出手段１３１は第１、第２の
指数平滑手段１２１，１２２の出力信号x_a1（ｔ）
とx_a2（ｔ）より低周波の卓越周期_L（ｔ）と振幅
a_L（ｔ）を算出する。第２の卓越周期算出手段１
３２は第３、第４の指数平滑手段１２３，１２４
の出力信号x_a3（ｔ）とx_a4（ｔ）より高周波の卓越
周期_H（ｔ）と振幅a_H（ｔ）を算出する。ここで、
低周波側の卓越周期算出手段１３１には、第１の
センサ１からの信号を入力しており、高周波側の
卓越周期算出手段１３２には第４のセンサ４から
の信号に入力している。第１のセンサ１の固有周
期は0.1Hz、第４のセンサ４の固有周期は５Hzで
あり、第１から第４のセンサのうちれぞれ最小、
最大値を有しており、このため第１の卓越周期算
出手段１３１は低周波側となり、第２の卓越周期
算出手段１３２は高周波側となる。

電源距離算出手段１４は、第１及び第２の卓越
周期算出手段１３１，１３２からの卓越周波数_L
（ｔ）、_H（ｔ）と振幅a_L（ｔ）、a_H（ｔ）より、震
源
距離Ｒ（ｔ）を次式で算出する。

Ｒ（ｔ）＝Ｑ・V_p・log（a_L（ｔ）／a_H（ｔ））／π（l
oge）（_H−_L）（29）ここで、Ｑは波動の減衰を示す定数で、たとえ
ば150を採用する。V_pは伝播経路上のＰ波の平均
速度を示し、たとえば6.5Km／秒を採用する。π
は円周率である。

（29）式で示す震源距離の推定式は次のように
して求められる。

震源からＲKmのときの振幅をａ、R_OKmのとき
の振幅をa_Oとすると、ａは次の式で示される。

ａ＝（a_O／√_O）e^-〓^R （30） αは定数で、α＝2πh／λ＝ωh／V_p＝2πh／V_p λ＝V_p （30）式の対数をとると log ａ＝log a_O−αRlog ｅ−0.5logR_O ＝log a_O−2πhR／V_plog ｅ−１／２logR_O ここで、２台のセンサで観測した振幅をa_L、a_H
とすると log a_L＝log a_O −2πh（log ｅ）Ｒ／V_p・_L−１／２logR_O （31） log a_H＝log a_O −2πh（loge）Ｒ／V_p・_H−１／２logR_O （32）（31）、（32）式よりa_Oを消去すると、 log a_L−log a_H＝2πh（log ｅ）Ｒ／V_p（_H−_L）（33） ∴Ｒ＝V_plog（a_L／a_H）／2πh（log ｅ）・（_H−_L
）＝QV_plog（a_L／a_H）／π（log ｅ）（_H−_L）（34
）但し、_H＞_L、a_L＞a_Hである。

震源距離出力手段１５は、震源距離算出手段１
４で算出したその時刻の震源距離Ｒ（ｔ）を外部
に出力する。

地震は断層が破壊するときの振動であり、断層
が大きければ破壊に時間を要する。また、震源は
現在破壊が進行しているところであり、これは常
に移動している。したがつて、正確には震源距離
は断層運動が継続している間変動するものである
が、観測点と震央の離れが断層の大きさより十分
大きいと、推定される震源距離はあまり変動しな
いが、断層の大きさに比べて近いところに観測点
があると、震源距離は大きく変動することが考え
られる。したがつて、むしろ本発明は従来見逃し
ていたこのような変化をも促えることができ、本
発明の目的は十二分に達成される。

（第２の実施例）第２の実施例は、センサを３台用意しておき、
固有周期の低いセンサの出力を第１の標本化手段
１２１に、固有周期の高いセンサの出力を第４の
標本化手段１２４に、また前記２つの固有周期に
対し中間の固有周期を有するセンサの出力を第２
及び第３の標本化手段１２２及び１２３に、各々
入力する。このようにして、第１の実施例よりセ
ンサを１台少なくして構成することができる。

（第３の実施例）第４図に、第３の実施例のブロツク図を示す。
第４図において、５はアンプ、６はバツフアアン
プ、７１，７２，７３は特性変換装置であり、他
は第２図のものと同様のものである。

第３の実施例においては、センサを１台設置し
ており、当該センサ１の信号の特性を変換する特
性変換装置７１，７２，７３を、バツフアアンプ
６と第２、第３、第４の標本化手段１１２，１１
３，１１４との間に設置している。特性変換装置
７１，７２，７３はその出力を、センサ１の固有
周期とは異なつた固有周期を持つたセンサから出
力されたかのように、標本化手段１１２，１１
３，１１４に与えるものであり、例えば第１の実
施例と同様にセンサ１の固有周期を0.1Hzとし、
固有周期に関し第１の特性変換装置７１を0.5Hz、
第２の特性変換装装置７２を１Hz、第３の特性変
換装置７３を５Hzであるように標本化手段１１
１，１１２，１１３及び１１４の各々に入力す
る。この特性変換装置７１，７２，７３は、特開
昭56−46479（特願昭54−122559）として公開され
ているものである。すなわち、第５図に示すよう
に、振動波形が電気信号で入力され、任意に設定
可能な周波数特性と制動特性を有し、これらに基
づいてセンサをシユミレートする仮想地震計機能
と、センサと仮想地震計機能の形式、特性によつ
て決定される定数を前記センサの出力に掛けあわ
せて、これを仮想地震計機能へ入力する機能と、
仮想地震計機能のシユミレート結果によつて出力
される各種の出力信号にセンサの周波数特性と制
動特性から導かれる定数を掛け、これらを加え合
わせる機能とを備えた装置とによつて構成される
ものである。第５図において、５０１はセンサ、
５０２はｎ倍器、５０３は仮想地震計、５０４は
加算器である。このような構成によつて得られた
特性変換装装置７１，７２，７３の出力とバツフ
アアンプ６の出力とを各々標本化手段１１１，１
１２，１１３，１１４に送出することにより、以
下第１の実施例と同様にして、震源距離を得るこ
とができる。

また、第１又は第２の特性変換装置７１又は７
２のいずれか一方、例えば第２の特性変換装置７
２を省略して、第１の特性変換装置７１の出力を
第２及び第３の標本化手段１１２，１１３の共通
の入力とすることにより、第２の実施例と同様に
震源距離を得ることができる。

（発明の効果）本発明によれば、地震波の震源距離を、１地点
の観測によつて地震波の到来直後に算出すること
ができ、従来のように多数の地震計を広範囲に配
置することが不要になるという利点がある。ま
た、リアルタイムで震源距離を知ることができる
ので、地震の警報システムの判定材料として利用
することができる。

【図面の簡単な説明】

第１図は、一の地震を異なる固有周期を持つ地
震計で同時に観測したときの観測例を示す波形
図、第２図は第１の実施例を示すブロツク図、第
３図は実施例の処理の概略フローチヤート、第４
図は第３の実施例を示すブロツク図、第５図は特
性変換装置を説明するブロツク図である。１，２，３，４……センサ、５１，５２，５
３，５４……アンプ、６１，６２，６３，６４…
…バツフアアンプ、７１，７２，７３……特性変
換装置、８……処理装置、１１１，１１２，１１
３，１１４……標本化手段、１２１，１２２，１
２３，１２４……指数平滑化手段、１３１，１３
２……卓越周期算出手段、１４……震源距離算出
手段、１５……震源距離出力手段。

Claims

【特許請求の範囲】

１固有周期の異なる２種の地震計による地震観
測出力を１組とし、固有周期が最も長い地震計が
含まれる組と固有周期が最も短い地震計が含まれ
る組の２組の観測出力を入力し、各組の観測出力
の指数平滑値を算出する手段と、各組の地震計の
観測波形の前記指数平滑値と各組の地震計の固有
周期とにより、高低２つの周波数帯域における卓
越周波数とその成分振幅をそれぞれ算出する卓越
周波数算出手段と、前記２組の卓越周波数算出手
段より得た２つの成分振幅の比の対数と前記２組
の卓越周波数算出手段より得た２つの卓越周波数
の差との比を算出する手段とを有して観測地点の
震源距離を得ることを特徴とする震源距離算出装
置。