JPS59195178A - Earthquake deciding method - Google Patents
Earthquake deciding methodInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims description 14
- 238000006243 chemical reaction Methods 0.000 abstract description 2
- 230000009466 transformation Effects 0.000 abstract 2
- 238000010586 diagram Methods 0.000 description 3
- 239000013256 coordination polymer Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
Classifications
-
- G01V1/01—
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
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.
第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)
タに集め、地震発生の判定を行う地震判定方式において
、前記観測信号に・ついて、常時高速フーリエ変換を行
いその判定結果と、移動平均による判定結果とを併用し
、両者で地震検出を行った時地震波発生の時刻を求める
ことを特徴とする地震判定方式。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.
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)
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)
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 |
-
1983
- 1983-04-21 JP JP58070599A patent/JPS59195178A/en active Granted
Patent Citations (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 |
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)
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 |