JPS59195178A - Earthquake deciding method - Google Patents

Earthquake deciding method

Info

Publication number
JPS59195178A
JPS59195178A JP58070599A JP7059983A JPS59195178A JP S59195178 A JPS59195178 A JP S59195178A JP 58070599 A JP58070599 A JP 58070599A JP 7059983 A JP7059983 A JP 7059983A JP S59195178 A JPS59195178 A JP S59195178A
Authority
JP
Japan
Prior art keywords
time
earthquake
moving average
decision
processing 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.)
Granted
Application number
JP58070599A
Other languages
Japanese (ja)
Other versions
JPH057674B2 (en
Inventor
Kiyoko Miyashita
宮下 清子
Tadamitsu Ryu
忠光 龍
Hideaki Arisaka
有坂 英明
Nobuaki Hashioka
橋岡 伸明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujitsu Ltd
Original Assignee
Fujitsu Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Fujitsu Ltd filed Critical Fujitsu Ltd
Priority to JP58070599A priority Critical patent/JPS59195178A/en
Publication of JPS59195178A publication Critical patent/JPS59195178A/en
Publication of JPH057674B2 publication Critical patent/JPH057674B2/ja
Granted legal-status Critical Current

Links

Classifications

    • G01V1/01

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

PURPOSE:To enable decision of the time when an earthquake wave arises with high accuracy by making use of a decision by a high speed Fourier transformation in combination with a decision by a moving average. CONSTITUTION:The observation signals from a seimsograph EQ and a telemeter device TM are collected in a center DTC and is subjected to an A/D conversion by a receiver TRV, from which the signal is supplied to a moving average processing circuit TD and a high speed Fourier transformation processing circuit FFT. An integrated amplitude value A1 and an integrated amplitude value A2 at a noise level are compared when the signal exceeding a threshold level with respect to the frequency level W1-W2 specific to the earthquake wave is inputted are compared and the time t0 when (A2Xbeta)<A1 is attained is decided as the time when the earthquake arises.

Description

【発明の詳細な説明】 (1)発明の技術分野 本発明は地震波発生の時刻を精度高(判定する方式に関
する。
DETAILED DESCRIPTION OF THE INVENTION (1) Technical Field of the Invention The present invention relates to a method for determining the time of seismic wave generation with high precision.

(2)従来技術と問題点 従来から地震発止時刻を1′1j定するため、複数の地
震波観測テレメータ装置におりJる観測信号をデータセ
ンタに集め、移動平均法により地震波発生の時刻を判定
していた。すなわち第1図に示すように地震波を電圧レ
ベルに変換して観測するとき、1秒間程度の短時間ao
−ar46平均の電圧レベルa1と、40秒間程度の長
時間す、+。〜ba平均値尻とを比較し、 bt−xα<aL となるときの短時間aw ””’ aoのいずれかの開
始時刻を地震波発生の時刻と判定する。ただしαはlよ
り大きな任意の定数とし、起こり易い地震の大きさに対
応させて選定する。第1図でばaQ、o+  a、’)
・7−−−−−−−と推移し、ao において所定のス
レショルドレベルを超えたとき、長時間平均値bL は
急激な変化を起こさないので、短時間a。における平均
値の方が急に太き(変化する。この場合、長時間とする
長さ、αの大きさについて、理論的な妥当値を求めるこ
とが困難であり、経験により適当な値を使用していた。
(2) Conventional technology and problems Conventionally, in order to determine the time of earthquake onset, observation signals from multiple seismic wave observation telemeter devices are collected at a data center, and the time of earthquake wave occurrence is determined using the moving average method. Was. In other words, as shown in Figure 1, when converting seismic waves to voltage levels and observing them, ao
-Ar46 average voltage level a1 and a long time of about 40 seconds, +. ~ba is compared with the end of the average value, and the start time of any of the short periods aw ""' ao when bt-xα<aL is determined to be the time of seismic wave generation. However, α is an arbitrary constant larger than l, and is selected in accordance with the magnitude of earthquakes that are likely to occur. In Figure 1, aQ, o+ a,')
・When ao exceeds a predetermined threshold level, the long-term average value bL does not change rapidly, so the short-term a. The average value of is suddenly thicker (changes). In this case, it is difficult to obtain theoretically appropriate values for the length of time and the size of α, and appropriate values are used based on experience. Was.

気象変化に対応してαを定めることが良いという検討も
されている。したがって短時間ao の開始時刻と言っ
ても精度が高くないという欠点があった。
It is also being considered that it would be better to determine α in response to weather changes. Therefore, there is a drawback that the accuracy of the start time of the short time ao is not high.

(3)発明の目的 本発明の目的は前述の欠点を改善し、後攻のデータ処理
を併せて行うごとにより、精度高(地震波発生の時刻を
判定する方式を提供することにある。
(3) Purpose of the Invention The purpose of the present invention is to improve the above-mentioned drawbacks and to provide a method for determining the time of seismic wave generation with high accuracy by also performing subsequent data processing.

4)発明の構成 前述の目的を達成するための本発明の構成は、地震波観
測テレメータ装置からの観測信号をデータセンタに集め
、地震発生の判定を行・う地震判定方式において、前記
観測信号について、常時高速フーリエ変換を行いその判
定結果と1.移動平均による判定結果とを併用し、両者
で地震発生を検出した時、地震波発生の時刻を求めるこ
とである。
4) Structure of the Invention The structure of the present invention to achieve the above-mentioned object is to collect observation signals from a seismic wave observation telemeter device in a data center and determine whether or not an earthquake has occurred. , always perform fast Fourier transform and the judgment results and 1. The method is to use both the determination results based on the moving average and determine the time of seismic wave generation when both detect the occurrence of an earthquake.

(5ン発明の実施例 第2図は本発明の一実施例の構成を示す図で、地震計E
Qとそのテレメータ装置TMからの観/l!11信号は
データセンタDTCに集められる。データセンタDTC
ではテレメータ受信装置TRVにおいてA/D変換され
、得られたディジタル信号について判定を行う。ディジ
タル信号を前述の移動平均法による移動平均処理回路T
Dと、本発明による高速フーリエ変換処理回路FFTに
入力する。
(Embodiment 5 of the Invention Figure 2 is a diagram showing the configuration of an embodiment of the invention.
View from Q and its telemeter device TM/l! 11 signals are collected at the data center DTC. Data center DTC
Now, the telemeter receiving device TRV performs A/D conversion, and a determination is made on the obtained digital signal. Moving average processing circuit T using the above-mentioned moving average method
D and is input to the fast Fourier transform processing circuit FFT according to the present invention.

高速フーリエ変換処理は成る時刻のテイジクル信号につ
いて、第3図A、Bに示すように、周波数対振幅値を求
めることである。即ぢ地震波の特許の周波数範囲W1〜
W2について、第3図Aに示す前述のようなスレショル
ドレベルを超える信号が入力したときに、振幅積分値A
lと1、第3図Bに示す所謂平常時の雑音レベルにおけ
る成る時刻t1の振幅積分値A2とを比較する。そして
A2Xβ<Al (ただしβは1より大きな任意の定数とする)とな9た
ときの時刻10を地震波発生の時刻と判定する。中央処
理装置CPUは前述と同様の移動平均法により地震波発
生を検出し、更にFFT法により地震波発生を検出した
とき、真の地震波発生と判定し、この時地震発生時刻は
移動平均法により求めた時刻又はFFT法により求めた
時刻のいずれを用いても良い。
The fast Fourier transform process is to obtain the frequency versus amplitude value for the tagged signal at the given time, as shown in FIGS. 3A and 3B. Ikuji seismic wave patent frequency range W1 ~
Regarding W2, when a signal exceeding the threshold level as shown in FIG. 3A is input, the amplitude integral value A
1 and the amplitude integral value A2 at time t1 at the so-called normal noise level shown in FIG. 3B are compared. Then, the time 10 when A2Xβ<Al (where β is an arbitrary constant greater than 1) is determined to be the time of seismic wave generation. When the central processing unit CPU detects the occurrence of seismic waves using the same moving average method as described above and further detects the occurrence of seismic waves using the FFT method, it determines that a true seismic wave has occurred, and at this time, the time of earthquake occurrence is determined using the moving average method. Either the time or the time determined by the FFT method may be used.

(6)発明の効果 このようにして本発明によると、高速フーリエ変換によ
る判定を従来の移動平均による判定と併用することによ
り、異なる地点からの地震波を判定するときに、途中の
地質の状態や、気象条件に基づく誤差を減少し、精度の
高い地震波発生時刻を求めることができ、したがって地
震発生地点をより正確に求めることができる。
(6) Effects of the Invention As described above, according to the present invention, by using determination by fast Fourier transform together with determination by conventional moving average, when determining seismic waves from different points, geological conditions along the way can be , it is possible to reduce errors based on weather conditions, to obtain highly accurate seismic wave generation times, and therefore to more accurately determine earthquake occurrence points.

【図面の簡単な説明】[Brief explanation of drawings]

第1図は地震波を移動平均法により判定する説明図、 第2図は本発明の一実施例の構成を示す図第3図は高速
フーリエ変換処理の説明図である。 ’T”RV−一移動平均法処理回路 F F T−高速フーリエ変換処理回路CP U−中央
処理装置 特許出願人    冨士通株式会社 代理人     弁理士 鈴木栄祐 ^ A 第1図
FIG. 1 is an explanatory diagram for determining seismic waves using the moving average method. FIG. 2 is an illustration of the configuration of an embodiment of the present invention. FIG. 3 is an explanatory diagram of fast Fourier transform processing. 'T'RV - Moving average processing circuit F F T - Fast Fourier transform processing circuit CP U - Central processing unit Patent applicant Fujitsu Co., Ltd. Agent Patent attorney Eisuke Suzuki A Figure 1

Claims (1)

【特許請求の範囲】[Claims] 地震波観測テレメータ装置からの観測信号をデータセン
タに集め、地震発生の判定を行う地震判定方式において
、前記観測信号に・ついて、常時高速フーリエ変換を行
いその判定結果と、移動平均による判定結果とを併用し
、両者で地震検出を行った時地震波発生の時刻を求める
ことを特徴とする地震判定方式。
In an earthquake determination method that collects observation signals from seismic wave observation telemeter equipment in a data center and determines the occurrence of an earthquake, the observation signals are constantly subjected to fast Fourier transform and the determination results are combined with the determination results based on a moving average. An earthquake determination method characterized by being used together to determine the time of seismic wave generation when earthquakes are detected by both.
JP58070599A 1983-04-21 1983-04-21 Earthquake deciding method Granted JPS59195178A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP58070599A JPS59195178A (en) 1983-04-21 1983-04-21 Earthquake deciding method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP58070599A JPS59195178A (en) 1983-04-21 1983-04-21 Earthquake deciding method

Publications (2)

Publication Number Publication Date
JPS59195178A true JPS59195178A (en) 1984-11-06
JPH057674B2 JPH057674B2 (en) 1993-01-29

Family

ID=13436186

Family Applications (1)

Application Number Title Priority Date Filing Date
JP58070599A Granted JPS59195178A (en) 1983-04-21 1983-04-21 Earthquake deciding method

Country Status (1)

Country Link
JP (1) JPS59195178A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62103596A (en) * 1985-10-31 1987-05-14 Radio Res Lab Seismic forecast by potential difference or dynamic spectrum of radio wave
WO1997002499A1 (en) * 1995-07-05 1997-01-23 Geoforschungszentrum Potsdam Earthquake monitoring
EP1540616A4 (en) * 2002-08-30 2006-06-07 Seismic Warning Systems Inc Sensor apparatus and method for detecting earthquake generated p- waves and generating a responsive control signal

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS53115281A (en) * 1977-03-18 1978-10-07 Oki Electric Ind Co Ltd Earthquake discrimination system
JPS56166481A (en) * 1980-05-28 1981-12-21 Fujitsu Ltd Deciding method of earthquake
JPS5757273A (en) * 1980-09-25 1982-04-06 Japanese National Railways<Jnr> Early sensing and warning system for earthquake

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS53115281A (en) * 1977-03-18 1978-10-07 Oki Electric Ind Co Ltd Earthquake discrimination system
JPS56166481A (en) * 1980-05-28 1981-12-21 Fujitsu Ltd Deciding method of earthquake
JPS5757273A (en) * 1980-09-25 1982-04-06 Japanese National Railways<Jnr> Early sensing and warning system for earthquake

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62103596A (en) * 1985-10-31 1987-05-14 Radio Res Lab Seismic forecast by potential difference or dynamic spectrum of radio wave
WO1997002499A1 (en) * 1995-07-05 1997-01-23 Geoforschungszentrum Potsdam Earthquake monitoring
DE19680510B4 (en) * 1995-07-05 2004-07-22 Geoforschungszentrum Potsdam Seismic monitoring system for earthquake forecasting - evaluates detected seismic parameters to locate minor tremors with spatial and chronological quantification of latter for earthquake prediction
EP1540616A4 (en) * 2002-08-30 2006-06-07 Seismic Warning Systems Inc Sensor apparatus and method for detecting earthquake generated p- waves and generating a responsive control signal

Also Published As

Publication number Publication date
JPH057674B2 (en) 1993-01-29

Similar Documents

Publication Publication Date Title
EP0153787B1 (en) System of analyzing human speech
CN101881628A (en) Detecting method of weak periodic signal based on chaotic system and wavelet threshold denoising
US5257309A (en) Dual tone multifrequency signal detection and identification methods and apparatus
US20190376874A1 (en) A method of estimating the number of modes for the sparse component analysis based modal identification
CN105865611B (en) A kind of method and device adjusting fiber-optic vibration detection threshold value
CN113674763B (en) Method, system, device and storage medium for identifying whistle by utilizing line spectrum characteristics
CN109377982B (en) Effective voice obtaining method
JPS59195178A (en) Earthquake deciding method
GB2049935A (en) Apparatus for monitoring vibrations in rotating machinery
CN115857013A (en) Method for calculating self-noise of seismometer by using improved welch method
CN107548007B (en) Detection method and device of audio signal acquisition equipment
US9384641B2 (en) Glass breakage detection system and method
CN114201991A (en) Partial discharge signal detection method and system based on ultrasonic sensor array
Kanai et al. Detection and discrimination of flaws in ball bearings by vibration analysis
CN113723207A (en) Method for detecting mutation signal based on histogram distance
US20060100828A1 (en) Impulse event separating apparatus and method
Diamant Computationally efficient calculations of target performance of the normalized matched filter detector for hydrocoustic signals
Wang et al. Dual Gaussian attenuation model of ultrasonic echo and its parameter estimation
JP2603758B2 (en) Cavitation noise detection method and device
KR940017487A (en) Micro Frequency Shift Detection Method
CN113109705B (en) GIS mechanical resonance spectrum analysis method based on sensitivity analysis of different resonance points
CN114173273B (en) Microphone array detection method, related device and readable storage medium
JP3241486B2 (en) Measurement seismic intensity meter
Miyano Estimating noise level in dynamical behavior for detecting mechanical damage
Ma et al. A Fast SNR-based Vibration Events Detection Algorithm for AETA Geo-acoustic Data