JPH057674B2 - - Google Patents
Info
- Publication number
- JPH057674B2 JPH057674B2 JP58070599A JP7059983A JPH057674B2 JP H057674 B2 JPH057674 B2 JP H057674B2 JP 58070599 A JP58070599 A JP 58070599A JP 7059983 A JP7059983 A JP 7059983A JP H057674 B2 JPH057674 B2 JP H057674B2
- Authority
- JP
- Japan
- Prior art keywords
- earthquake
- time
- occurrence
- exceeded
- comparison circuit
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/01—Measuring or predicting earthquakes
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Description
(1) 発明の技術分野
本発明は地震波発生の時刻、即ち、地震波発生
の時点（time point）を精度高く判定する方式に
関する。 (2) 従来技術と問題点 従来から地震発生時刻を判定するため、複数の
地震波観測テレメータ装置における観測信号をデ
ータセンタに集め、移動平均法により地震波発生
の時刻を判定していた。すなわち第１図に示すよ
うに地震波を電圧レベルに変換して観測すると
き、１秒間程度の短時間a0〜a40平均の電圧レベ
ルaLと、40秒間程度の長時間b40〜b0平均値bLと
を比較し、 bL×α＜aL となるときの短時間a40〜a0のいずれかの開始時
刻を地震波発生の時刻と判定する。ただしαは１
より大きな任意の定数とし、起こり易い地震の大
きさに対応させて選定する。第１図ではa40、a39
……と推移し、a0において所定のスレシヨルドレ
ベルを超えたとき、長時間平均値bLは急激な変化
を起こさないので、短時間a0における平均値の方
が急に大きく変化する。この場合、長時間とする
長さ、αの大きさについて、理論的な妥当値を求
めることが困難であり、経験により適当な値を使
用していた。気象変化に対応してαを定めること
が良いという検討もされている。また、時刻ａに
おいて、地震波と同一周波数帯域でノイズが発生
したときデータセンタは地震発生と判断するの
で、実際の地震とノイズとを区別できずに、地震
発生時刻を誤るという欠点があつた。 (3) 発明の目的 本発明の目的は前述の欠点を改善し、複数のデ
ータ処理を併せて行うことにより、精度高く地震
波発生の時刻を判定する方式を提供することにあ
る。 (4) 発明の構成 前述の目的を達成するための本発明の構成は、
地震波観測テレメータ装置からの観測信号をデー
タセンタに集め、該観測信号の変化から地震発生
の判定を行う方式において、前記観測信号を常時
高速フーリエ変換し、地震波特有の周波数範囲で
のスペクトルの積算値を平常時の積算値と比較す
る第１比較回路と、前記観測信号を移動平均演算
し平常時の値と比較する第２比較回路と、を具備
し、前記第１比較回路で規定レベルを超えたこと
を検出し、更に後記第２比較回路においても規定
レベルを超えたことを検出したとき地震発生と判
定し、その何れの比較回路によるか予め定めてお
いた比較回路により、規定レベルを超えたと判定
した時刻を地震発生時刻とする。 (5) 発明の実施例 第２図は本発明の一実施例の構成を示す図で、
地震計EQとそのテレメータ装置TMからの観測
信号はデータセンタDTCに集められる。データ
センタDTCではテレメータ受信装置TRVにおい
てＡ／Ｄ変換され、得られたデイジタル信号につ
いて判定を行う。デイジタル信号を前述の移動平
均法による移動平均処理回路TDと、本発明によ
る高速フーリエ変換処理回路FFTに入力する。
高速フーリエ変換処理は或る時刻のデイジタル信
号について、第３図Ａ，Ｂに示すように、周波数
対振幅値を求めることである。即ち地震波の特有
の周波数範囲W1〜W2について、第３図Ａに示す
前述のようなスレシヨルドレベルを超える信号が
入力したときに、振幅積分値Ａ１と、第３図Ｂに
示す所謂平常時の雑音レベルにおける或る時刻t1
の振幅積分値Ａ２とを比較する。そして Ａ２×β＜Ａ１ （ただしβは１より大きな任意の定数とする） となつたときの時刻t0を地震波発生の時刻と判定
する。中央処理装置CPUは、高速フーリエ変換
処理を行うと共に、前記積算値Ａ１，Ａ２の比較
を行う第１比較回路として動作する。 更に中央処理装置CPUは従来の装置と同様に
移動平均法による判定を行うため、第２比較回路
としても動作する。 その結果中央処理装置CPUは前述と同様の移
動平均法により地震波発生を検出し、更にFFT
法により地震波発生を検出したとき、真の地震波
発生と判定し、この時地震発生時刻は移動平均法
により求められた時刻又はFFT法により求めた
時刻のうち、例えば後者のFFT法により求めた
時刻とするように予め定めておく。若し、移動平
均法によつて地震発生と判断したが、FFT法に
よつては発生としなければ、外来雑音などによる
こととして、中央処理装置はその後の処理を行わ
ない。更にFFTにより演算するとき、地震波特
有の帯域におけるノイズを常時調べているから、
地震波特有の帯域における観測信号について前記
ノイズ信号を差し引いて後、演算処理を行うこと
で、地震波に混入している雑音の影響を取り除く
ことが出来る。 (6) 発明の効果 このようにして本発明によると、高速フーリエ
変換による判定を従来の移動平均による判定と併
用することにより、異なる地点からの地震波を判
定するときに、途中の地質の状態や、気象条件に
基づく誤差を減少し、精度の高い地震波発生時刻
を求めることができる。また突発的な外来雑音が
到来し、バツクグランドノイズが変化していると
きでも、それらの影響が容易に取り除けるため、
時刻観測が正確となり、地震発生地点をより正確
に求めることができる。
の時点（time point）を精度高く判定する方式に
関する。 (2) 従来技術と問題点 従来から地震発生時刻を判定するため、複数の
地震波観測テレメータ装置における観測信号をデ
ータセンタに集め、移動平均法により地震波発生
の時刻を判定していた。すなわち第１図に示すよ
うに地震波を電圧レベルに変換して観測すると
き、１秒間程度の短時間a0〜a40平均の電圧レベ
ルaLと、40秒間程度の長時間b40〜b0平均値bLと
を比較し、 bL×α＜aL となるときの短時間a40〜a0のいずれかの開始時
刻を地震波発生の時刻と判定する。ただしαは１
より大きな任意の定数とし、起こり易い地震の大
きさに対応させて選定する。第１図ではa40、a39
……と推移し、a0において所定のスレシヨルドレ
ベルを超えたとき、長時間平均値bLは急激な変化
を起こさないので、短時間a0における平均値の方
が急に大きく変化する。この場合、長時間とする
長さ、αの大きさについて、理論的な妥当値を求
めることが困難であり、経験により適当な値を使
用していた。気象変化に対応してαを定めること
が良いという検討もされている。また、時刻ａに
おいて、地震波と同一周波数帯域でノイズが発生
したときデータセンタは地震発生と判断するの
で、実際の地震とノイズとを区別できずに、地震
発生時刻を誤るという欠点があつた。 (3) 発明の目的 本発明の目的は前述の欠点を改善し、複数のデ
ータ処理を併せて行うことにより、精度高く地震
波発生の時刻を判定する方式を提供することにあ
る。 (4) 発明の構成 前述の目的を達成するための本発明の構成は、
地震波観測テレメータ装置からの観測信号をデー
タセンタに集め、該観測信号の変化から地震発生
の判定を行う方式において、前記観測信号を常時
高速フーリエ変換し、地震波特有の周波数範囲で
のスペクトルの積算値を平常時の積算値と比較す
る第１比較回路と、前記観測信号を移動平均演算
し平常時の値と比較する第２比較回路と、を具備
し、前記第１比較回路で規定レベルを超えたこと
を検出し、更に後記第２比較回路においても規定
レベルを超えたことを検出したとき地震発生と判
定し、その何れの比較回路によるか予め定めてお
いた比較回路により、規定レベルを超えたと判定
した時刻を地震発生時刻とする。 (5) 発明の実施例 第２図は本発明の一実施例の構成を示す図で、
地震計EQとそのテレメータ装置TMからの観測
信号はデータセンタDTCに集められる。データ
センタDTCではテレメータ受信装置TRVにおい
てＡ／Ｄ変換され、得られたデイジタル信号につ
いて判定を行う。デイジタル信号を前述の移動平
均法による移動平均処理回路TDと、本発明によ
る高速フーリエ変換処理回路FFTに入力する。
高速フーリエ変換処理は或る時刻のデイジタル信
号について、第３図Ａ，Ｂに示すように、周波数
対振幅値を求めることである。即ち地震波の特有
の周波数範囲W1〜W2について、第３図Ａに示す
前述のようなスレシヨルドレベルを超える信号が
入力したときに、振幅積分値Ａ１と、第３図Ｂに
示す所謂平常時の雑音レベルにおける或る時刻t1
の振幅積分値Ａ２とを比較する。そして Ａ２×β＜Ａ１ （ただしβは１より大きな任意の定数とする） となつたときの時刻t0を地震波発生の時刻と判定
する。中央処理装置CPUは、高速フーリエ変換
処理を行うと共に、前記積算値Ａ１，Ａ２の比較
を行う第１比較回路として動作する。 更に中央処理装置CPUは従来の装置と同様に
移動平均法による判定を行うため、第２比較回路
としても動作する。 その結果中央処理装置CPUは前述と同様の移
動平均法により地震波発生を検出し、更にFFT
法により地震波発生を検出したとき、真の地震波
発生と判定し、この時地震発生時刻は移動平均法
により求められた時刻又はFFT法により求めた
時刻のうち、例えば後者のFFT法により求めた
時刻とするように予め定めておく。若し、移動平
均法によつて地震発生と判断したが、FFT法に
よつては発生としなければ、外来雑音などによる
こととして、中央処理装置はその後の処理を行わ
ない。更にFFTにより演算するとき、地震波特
有の帯域におけるノイズを常時調べているから、
地震波特有の帯域における観測信号について前記
ノイズ信号を差し引いて後、演算処理を行うこと
で、地震波に混入している雑音の影響を取り除く
ことが出来る。 (6) 発明の効果 このようにして本発明によると、高速フーリエ
変換による判定を従来の移動平均による判定と併
用することにより、異なる地点からの地震波を判
定するときに、途中の地質の状態や、気象条件に
基づく誤差を減少し、精度の高い地震波発生時刻
を求めることができる。また突発的な外来雑音が
到来し、バツクグランドノイズが変化していると
きでも、それらの影響が容易に取り除けるため、
時刻観測が正確となり、地震発生地点をより正確
に求めることができる。
第１図は地震波を移動平均法により判定する説
明図、第２図は本発明の一実施例の構成を示す
図、第３図は高速フーリエ変換処理の説明図であ
る。 EQ……地震計、TM……テレメータ装置、
DTC……データセンタ、TRV……移動平均法処
理回路、FFT……高速フーリエ変換処理回路、
CPU……中央処理装置。
明図、第２図は本発明の一実施例の構成を示す
図、第３図は高速フーリエ変換処理の説明図であ
る。 EQ……地震計、TM……テレメータ装置、
DTC……データセンタ、TRV……移動平均法処
理回路、FFT……高速フーリエ変換処理回路、
CPU……中央処理装置。
Claims (1)
- 【特許請求の範囲】 １ 地震波観測テレメータ装置からの観測信号を
データセンタに集め、該観測信号の変化から地震
発生の判定を行う方式において、 前記観測信号を常時高速フーリエ変換し、地震
波特有の周波数範囲でのスペクトルの積算値を平
常時の積算値と比較する第１比較回路と、 前記観測信号を移動平均演算し平常時の値と比
較する第２比較回路と、を具備し、 前記第１比較回路で規定レベルを超えたことを
検出し、更に後記第２比較回路においても規定レ
ベルを超えたことを検出したとき地震発生と判定
し、 その何れの比較回路によるか予め定めておいた
比較回路により、規定レベルを超えたと判定した
時刻を地震波発生時刻とすること を特徴とする地震発生判定方式。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP58070599A JPS59195178A (ja) | 1983-04-21 | 1983-04-21 | 地震判定方式 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP58070599A JPS59195178A (ja) | 1983-04-21 | 1983-04-21 | 地震判定方式 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS59195178A JPS59195178A (ja) | 1984-11-06 |
| JPH057674B2 true JPH057674B2 (ja) | 1993-01-29 |
Family
ID=13436186
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP58070599A Granted JPS59195178A (ja) | 1983-04-21 | 1983-04-21 | 地震判定方式 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS59195178A (ja) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62103596A (ja) * | 1985-10-31 | 1987-05-14 | Radio Res Lab | 電位差又は電波のダイナミツク・スペクトルによる地震予報法 |
| US5890094A (en) * | 1995-07-05 | 1999-03-30 | Geoforschungszentrum Potsdam | Earthquake monitoring |
| WO2004021298A2 (en) * | 2002-08-30 | 2004-03-11 | Seismic Warning Systems, Inc. | Sensor apparatus and method for detecting earthquake generated p- waves and generating a responsive control signal |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS53115281A (en) * | 1977-03-18 | 1978-10-07 | Oki Electric Ind Co Ltd | Earthquake discrimination system |
| JPS6026472B2 (ja) * | 1980-05-28 | 1985-06-24 | 富士通株式会社 | 地震判定方式 |
| JPS608468B2 (ja) * | 1980-09-25 | 1985-03-02 | 日本国有鉄道 | 地震早期検知警報システム |
-
1983
- 1983-04-21 JP JP58070599A patent/JPS59195178A/ja active Granted
Also Published As
| Publication number | Publication date |
|---|---|
| JPS59195178A (ja) | 1984-11-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10386510B2 (en) | Earthquake Detection System and Method | |
| Sereno Jr et al. | Seismic detection capability at NORESS and implications for the detection threshold of a hypothetical network in the Soviet Union | |
| Tappert et al. | New full‐wave approximation for ocean acoustic travel time predictions | |
| MaMullen et al. | A simple rising-edge detector for time-of-arrival estimation | |
| JPH057674B2 (ja) | ||
| JP4465509B2 (ja) | 震央距離推定装置、震央距離推定システム及び震央距離推定方法 | |
| Wenz | Some periodic variations in low‐frequency acoustic ambient noise levels in the ocean | |
| CN112834774B (zh) | 一种阈值自适应转速信号处理系统及其方法 | |
| JP2024057455A (ja) | 震源距離推定装置、震源距離推定プログラム及び震源距離推定方法 | |
| US3717812A (en) | Real time analysis of waves | |
| AU2024264565B2 (en) | Method and device for detecting earthquakes | |
| JPH0682067B2 (ja) | 地震検出装置 | |
| JPH07181249A (ja) | 超音波距離測定装置 | |
| RU2796528C1 (ru) | Способ отбраковки шумов при магнитотеллурическом зондировании | |
| CN115372920B (zh) | 一种雷达结露补偿方法、装置、设备及存储介质 | |
| JPH0263194B2 (ja) | ||
| JPS587538A (ja) | ノツキング検出装置 | |
| KR100236234B1 (ko) | 차량용 쇽 업소버의 레그 판별 방법 | |
| Gupta et al. | Micropulsations in the Pc 3 to Pc 4 period range at four Canadian observatories | |
| KR101034537B1 (ko) | 지진신호 감지 및 보정장치 및 그 구동 방법 | |
| JPS5621059A (en) | Method and device for determining position in ae measurement | |
| JPS6095376A (ja) | 津波予測装置 | |
| JPH0311751Y2 (ja) | ||
| JPH01133198A (ja) | 地震警報方法 | |
| CN121049669A (zh) | 一种局部放电信号源定位方法及相关装置 |