JPH0236196B2 - - Google Patents
Info
- Publication number
- JPH0236196B2 JPH0236196B2 JP58064560A JP6456083A JPH0236196B2 JP H0236196 B2 JPH0236196 B2 JP H0236196B2 JP 58064560 A JP58064560 A JP 58064560A JP 6456083 A JP6456083 A JP 6456083A JP H0236196 B2 JPH0236196 B2 JP H0236196B2
- Authority
- JP
- Japan
- Prior art keywords
- amplitude
- exponential smoothing
- vertical
- value
- sensor
- 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
- 238000009499 grossing Methods 0.000 claims description 22
- 238000001514 detection method Methods 0.000 description 10
- 238000005070 sampling Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 5
- 238000004891 communication Methods 0.000 description 3
- 239000002131 composite material Substances 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000000034 method Methods 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
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
（技術分野）
本発明は、地震を検出しその震央方位を即時に
計算する震央方位推定装置である。
（従来技術）
１地点で震央方位を推定するためには、地震波
が到来したときの３成分（上下方向、南北方向、
東西方向）の初期の振幅を利用することが一般的
である。第１図（宇津（1976）、地震学、共立出
版より）のように南北方向成分の初期振幅ANと
東西方向成分の初期振幅AWを合成し、上下動の
初動が下向きのときはANとAWの合成方向が震央
方位となり、上下動の初動が下向きのときは前記
合成方向とは反対の方向が震央方位となる。３成
分とも初動が明瞭なときには、この方法も有効で
あるが、不明瞭なときには、90゜又は180゜の誤り
が出る可能性がある。
（発明の目的）
本発明の目的は、地動を検出するセンサからの
信号に基づいて、初動振幅が小さくても安定して
震央方位を連続的に算出することができるように
した震央方位推定装置を提供することである。
（発明の構成）
一般に地表面近くの地震波の速度が遅くなるた
め、地震波の入射角が垂直に近くなる。このた
め、地震波のＰ波部分では上下動成分の振幅が大
きく、Ｓ波部分では水平動成分の振幅が大きいこ
とが多い。上下動成分と水平動成分の比をとれば
Ｐ波とＳ波の区別に利用できることが考えられ
る。
単純に上下動成分と水平動成分の瞬時値、移動
平均値を使つては、Ｐ波とＳ波の区別を容易に行
うことはできない。これはデータの揺らぎやばら
つきなどノイズ成分の影響によるものである。本
発明では、揺らぎやばらつきを抑えるとともに振
動波形の包絡線形状（いわゆる振幅の変動）を連
続的に把握するため、指数平滑値を使つて、上下
動振幅と水平動振幅の比を算出すること、さらに
は、前記比の上昇によりＰ波の到来であると区別
する。
さらに詳しく説明すると、本発明の第１の構成
は、上下方向、東西方向および南北方向の地動を
検出するセンサと、前記センサから検出された上
下方向成分と南北方向成分との振幅値の積の指数
平滑値を算出する第１の指数平滑値算出手段と、
前記センサから検出された上下方向成分と東西方
向成分との振幅値の積の指数平滑値を算出する第
２の指数平滑値算出手段と、前記第１の指数平滑
値算出手段の出力と第２の指数平滑値算出手段の
出力とを比較する手段とからなることを特徴とす
る震央方位推定装置である。前記第１および第２
の指数平滑値算出手段は上下方向と同一位相で振
動する水平方向成分を強調し、不規則なノイズ成
分を減少させるためのものである。前記第１およ
び第２の指数平滑値算出手段の出力を比較する手
段は南北方向と東西方向との振幅の比率をみるた
めのものであり、これにより振動の方向、すなわ
ち震央方位を推定することができる。また、本発
明の第２の構成は、前記第１の構成とともに、上
下方向の振幅と水平方向の振幅との振幅比を算出
する上下水平振幅比算出手段とを有することを特
徴とする震央方位推定装置である。これは、前記
のごとく前記第１の構成の前提であるＰ波を検出
するためのものである。
（実施例）
第２図は本発明の一実施例を示すブロツク図で
ある。１は地動を検出するセンサであり、上下方
向、東西方向及び南北方向の３成分の地動を検出
する形式のものである。２は前記センサ１よりの
検出信号を増幅するアンプ、３はバツフアアン
プ、４は前記センサ３からの検出信号より前記３
成分の地動を処理して震央を算出する制御処理装
置、５は時刻情報を発生する時計装置、６は前記
制御処理装置４が算出した震央方位を表示する表
示装置、７は前記制御処理装置４が算出した震央
方位を通信回線８を介して伝送する通信装置であ
る。制御処理装置４は指数平滑値算出手段１１
２，１１４、震央方位算出手段１１６、上下水平
振幅比算出手段１２、地震検出手段１３と標本化
手段１４により構成される。第４図は本発明によ
り行なわれる手順を示すフローチヤートである。
センサ１は常時設置されている地点の地動を検
出し、これを電気信号に変換して制御処理装置４
に送出している。
制御処理装置４は、センサ１で検出されてアン
プ２、バツフアアンプ３を介して絶えず送られて
くる地動の各成分の情報を取り込んでいる。標本
化手段１４は前記各成分の情報を所定の時間間隔
（例えば1/50〜1/500秒）で標本化し、その過去何
回分かのサンプリング情報の平均値より入力デー
タの直流分であるオフセツトレベルの算出を行
う。このオフセツトレベルは時々刻々得られるサ
ンプリング情報によつて絶えず更新されている。
また標本化手段１４は、センサ１によつて検出さ
れて時々刻々送られてくる上下方向、東西方向、
南北方向の３つの成分の地動情報のサンプリング
値からオフセツトレベルを除去した値x1（ｔ）、x2
（ｔ）、x3（ｔ）に基づいて、ノイズレベルNX1、
NX2、NX3と指数平滑値ax1（ｔ）、ax2（ｔ）、ax3
（ｔ）を算出する。ここで、サンプリング値をxsi
（ｔ）、オフセツトレベルをxpiとすると、ti（ｔ）＝
xsi（ｔ）−xpi（ｉ＝１、２、３）である。ノイズレ
ベルNXiは｜xi（ｔ）｜の平均値で、後述するよう
に地震波を検出するレベルＡの決定に利用する。
指数平滑値axi（ｔ）は次式で算出する。
axi（ｔ）＝axi（ｔ−１）・αi＋xi 2（ｔ）………(1
)
αiは0.9程度の定数である。またaxi（ｔ−１）は
時刻ｔに対し１サンプル前のaxiである。標本化
器１４はxi（ｔ）、axi（ｔ）及びNXiを指数平滑値
算出手段１１２，１１４と上下水平振幅比算出手
段１２に送出する。震央方位θ（ｔ）は次のよう
に計算され、過去３秒程度の計算結果を記憶す
る。まず、第１の指数平滑値算出手段１１２で
は、上下動の振幅値x1（ｔ）と南北動の振幅値x3
（ｔ）の積の指数平均値xUDNS（ｔ）を次の第(2)式
により計算する。
xUDNS（ｔ）＝xUDNS（ｔ−１）・α＋
x1（ｔ）・x3（ｔ）………(2)
また、第２の指数平滑算出手段１１４では、上
下動の振幅値x1（ｔ）と東西動の振幅値x2（ｔ）の
積の指数平均値xUDEW（ｔ）を次の第(3)式により計
算する。
xUDEW（ｔ）＝xUDEW（ｔ−１）・α＋
x1（ｔ）・x2（ｔ）………(3)
次いで、第(2)式と第(3)式の結果を震央方位算出
手段１１６により東から北回りの角度θ0（ｔ）を
θ0（ｔ）＝tan-1（xUDNS（ｔ）／xUD
EW（ｔ））………(4)
なる式で計算する。次にθ0（ｔ）は−90゜から90゜の
範囲にあるので、xUDNS（ｔ）とxUDEW（ｔ）の符号
より、次の第１表のように４象限に分ける（第３
図参照）。
(Technical Field) The present invention is an epicenter azimuth estimation device that detects an earthquake and instantly calculates its epicenter azimuth. (Prior art) In order to estimate the epicenter direction at one point, three components (vertical direction, north-south direction,
It is common to use the initial amplitude in the east-west direction). As shown in Figure 1 (from Utsu (1976), Seismology, Kyoritsu Shuppan), the initial amplitude A N of the north-south component and the initial amplitude A W of the east-west component are combined, and when the initial vertical motion is downward, A The composite direction of N and A W becomes the epicenter azimuth, and when the initial vertical movement is downward, the direction opposite to the composite direction becomes the epicenter azimuth. This method is effective when the initial movements of all three components are clear, but when they are unclear, there is a possibility of a 90° or 180° error. (Object of the Invention) The object of the present invention is to provide an epicenter azimuth estimation device that can stably and continuously calculate the epicenter azimuth even if the initial motion amplitude is small, based on signals from a sensor that detects ground motion. The goal is to provide the following. (Structure of the Invention) Generally, the velocity of seismic waves near the ground surface becomes slow, so the angle of incidence of the seismic waves becomes close to vertical. Therefore, the amplitude of the vertical motion component is often large in the P wave portion of the seismic wave, and the amplitude of the horizontal motion component is often large in the S wave portion. It is conceivable that the ratio of the vertical motion component to the horizontal motion component can be used to distinguish between P waves and S waves. P waves and S waves cannot be easily distinguished by simply using the instantaneous values and moving average values of the vertical motion component and horizontal motion component. This is due to the influence of noise components such as fluctuations and variations in data. In the present invention, in order to suppress fluctuations and variations and to continuously grasp the envelope shape of the vibration waveform (so-called amplitude fluctuation), the ratio of the vertical motion amplitude to the horizontal motion amplitude is calculated using an exponential smoothing value. , Furthermore, the arrival of a P wave is distinguished from the increase in the ratio. To explain in more detail, the first configuration of the present invention includes a sensor that detects ground motion in the vertical direction, east-west direction, and north-south direction, and a product of the amplitude values of the vertical component and the north-south component detected by the sensor. a first exponential smoothing value calculation means for calculating an exponential smoothing value;
a second exponentially smoothed value calculation means for calculating an exponentially smoothed value of the product of the amplitude values of the vertical component and the east-west direction component detected by the sensor; and an output of the first exponentially smoothed value calculating means and a second and means for comparing the output of the exponential smoothing value calculation means. Said first and second
The exponential smoothing value calculation means is for emphasizing the horizontal component vibrating in the same phase as the vertical direction and reducing irregular noise components. The means for comparing the outputs of the first and second exponential smoothing value calculating means is for checking the ratio of amplitudes in the north-south direction and the east-west direction, and from this, the direction of vibration, that is, the epicenter direction can be estimated. Can be done. Further, a second configuration of the present invention is characterized in that it has, in addition to the first configuration, vertical horizontal amplitude ratio calculation means for calculating an amplitude ratio between an amplitude in the vertical direction and an amplitude in the horizontal direction. It is an estimation device. This is for detecting P waves, which is the premise of the first configuration as described above. (Embodiment) FIG. 2 is a block diagram showing an embodiment of the present invention. A sensor 1 detects ground motion, and is of a type that detects three components of ground motion: vertical, east-west, and north-south directions. 2 is an amplifier for amplifying the detection signal from the sensor 1, 3 is a buffer amplifier, and 4 is the amplification signal for the detection signal from the sensor 3.
5 is a clock device that generates time information; 6 is a display device that displays the epicenter direction calculated by the control processing device 4; 7 is the control processing device 4; This is a communication device that transmits the calculated epicenter direction via a communication line 8. The control processing device 4 includes an exponential smoothing value calculation means 11
2, 114, an epicenter azimuth calculation means 116, an upper and lower horizontal amplitude ratio calculation means 12, an earthquake detection means 13, and a sampling means 14. FIG. 4 is a flowchart illustrating the steps performed in accordance with the present invention. The sensor 1 detects ground motion at the point where it is constantly installed, converts it into an electrical signal, and sends it to the control processing device 4.
It is being sent to The control processing device 4 takes in information on each component of ground motion detected by the sensor 1 and constantly sent via the amplifier 2 and buffer amplifier 3. The sampling means 14 samples the information of each component at a predetermined time interval (for example, 1/50 to 1/500 seconds), and calculates the DC component of the input data from the average value of the past sampling information. Calculate the set level. This offset level is constantly updated based on sampling information obtained from time to time.
In addition, the sampling means 14 collects data detected by the sensor 1 and sent from time to time in the vertical direction, east-west direction,
Values obtained by removing the offset level from the sampled values of ground motion information for the three components in the north-south direction x 1 (t), x 2
(t), x 3 (t), the noise level NX 1 ,
NX 2 , NX 3 and exponential smoothing values ax 1 (t), ax 2 (t), ax 3
(t) is calculated. Here, the sampling value is x si
(t), and the offset level is x pi , then t i (t)=
x si (t)−x pi (i=1, 2, 3). The noise level NX i is the average value of |x i (t)|, and is used to determine the level A at which seismic waves are detected, as will be described later.
The exponential smoothing value ax i (t) is calculated using the following formula. ax i (t)=ax i (t-1)・α i +x i 2 (t)……(1
) α i is a constant of about 0.9. Further, ax i (t-1) is ax i one sample before time t. The sampler 14 sends x i (t), ax i (t) and NX i to the exponential smoothing value calculation means 112, 114 and the upper and lower horizontal amplitude ratio calculation means 12. The epicenter azimuth θ(t) is calculated as follows, and the calculation results for about the past 3 seconds are stored. First, the first exponential smoothing value calculation means 112 calculates the vertical motion amplitude value x 1 (t) and the north-south motion amplitude value x 3
The exponential average value x UDNS (t) of the product of (t) is calculated using the following equation (2). x UDNS (t)=x UDNS (t-1)・α+
x 1 (t)・x 3 (t) (2) In addition, the second exponential smoothing calculation means 114 calculates the amplitude value x 1 (t) of vertical motion and the amplitude value x 2 (t) of east-west motion. The exponential average value of the product x UDEW (t) is calculated using the following equation (3). x UDEW (t)=x UDEW (t-1)・α+
x 1 (t)・x 2 (t)……(3) Next, the results of equations (2) and (3) are used to calculate the northward angle θ 0 (t) from the east using the epicenter direction calculation means 116. θ 0 (t)=tan -1 (x UDNS (t)/x UD
EW (t))......(4) Calculate using the formula. Next, since θ 0 (t) is in the range of -90° to 90 °, it is divided into four quadrants (the third
(see figure).
【表】
xUDNS（ｔ）が−、xUDEW（ｔ）が−のときについ
て（従来技術）より第１表を説明する。xUD（ｔ）
が−のときxNS（ｔ）とxEW（ｔ）はともに＋とな
る。従つて、震央は第象限となる。xUD（ｔ）が
＋のとき、xNS（ｔ）とxEW（ｔ）はともに−とな
る。xUD（ｔ）が＋のときには、xNS（ｔ）とxEW
（ｔ）示す方位の逆となるため、同様に第象限
となる。同様に、xUDNS（ｔ）とxUDEW（ｔ）の符号
について全ての組合せからその震央方位の象限を
説明できる。
さらにθ0（ｔ）を北から東回りの角度（０〜
360゜）に座標変換し、震央方位θ（ｔ）を決定す
る。
逆正接（tan-1）の計算では、その都度近似式
を解くよりはメモリに０〜90゜に相当する数値を
記憶しておき、それと比較することにより度単位
で求めることができ、時間の節約となる。
震央方位算出手段１１６はさらに、震央方位θ
（ｔ）について過去３秒程度の平均値（ｔ）を
計して、震央方位θ（ｔ）と平均値（ｔ）を地
震波検出手段１３に送出する。
震央方位θ（ｔ）の計算結果例として、第５図、
第６図、第７図を示す。地震波が到来するとθ
（ｔ）が急激に変動し、震央方位を示す。
一方、前述の上下水平振幅比算出手段１２は、
次式で上下動振幅と水平動振幅の振幅比（VH）
を計算し、過去３秒程度のVH（ｔ）を記憶する。
上下水平振幅比算出手段１２は、さらに過去３
秒程度の平均値（ｔ）を計算して、振幅比
VH（ｔ）と平均値（ｔ）を地震波検出手段１
３に送出する。
第５図の場合、地震波が到来するとVH（ｔ）
は上昇し、次いでＳ波が到来するとVH（ｔ）は
下降する。第７図の場合、Ｐ波の振幅が小さくて
検出できず、Ｓ波が到来したとき初めて地震波と
検出した場合でもVH（ｔ）が下降するのでそれ
がＳ波であることを識別することができる。
地震波検出手段１３は前述の上下方向のサンプ
リング値からオフセツトレベルを除去した値の｜
x1（ｔ）｜が連続してｎサンプル以上だけ地震波検
出レベルＡを超えたことにより地震波の到来を検
出する。地震波検出レベルＡはx1（ｔ）のノイズ
レベルNX1より決定される。地震波を検知する
と、当該検知のｍサンプル前に計算し、記憶して
おいた前記上下水平振幅比出手段１２からの振幅
比の平均値（tl）と地震波の検知時のｍサン
プル後に得た振幅比の平均値（tl+2n）とを比
較する。比較した結果（tl）が（tl+2n）よ
り大きいと検知し地震波はＰ波であると判定す
る。Ｐ波の検出後、θ（ｔ）の平均値を震央方位
として、表示装置６に表示するとともに通信制御
装置７に送出する。
（発明の効果）
本発明によれば、振幅の積の指数平滑値を算出
しているので正確な震央方位を推定でき、また正
確にＰ波であるかＳ波であるか判別することがで
きる。[Table] Table 1 will be explained from (prior art) when x UDNS (t) is - and x UDEW (t) is -. x UD (t)
When is -, x NS (t) and x EW (t) both become +. Therefore, the epicenter is in the fourth quadrant. When x UD (t) is +, x NS (t) and x EW (t) are both -. When x UD (t) is +, x NS (t) and x EW
(t) Since the orientation is opposite to that shown, it is also in the fourth quadrant. Similarly, the quadrant of the epicenter direction can be explained from all combinations of the signs of x UDNS (t) and x UDEW (t). Furthermore, θ 0 (t) is the angle from north to east (0~
360°) and determine the epicenter azimuth θ(t). When calculating the arctangent (tan -1 ), rather than solving an approximate equation each time, you can store numerical values corresponding to 0 to 90° in memory and compare them with that value to calculate the value in degrees. It saves money. The epicenter azimuth calculation means 116 further calculates the epicenter azimuth θ
(t), the average value (t) over the past 3 seconds or so is measured, and the epicenter azimuth θ(t) and the average value (t) are sent to the seismic wave detection means 13. Figure 5 shows an example of the calculation result of the epicenter azimuth θ(t).
Figures 6 and 7 are shown. When the seismic wave arrives, θ
(t) changes rapidly, indicating the epicenter direction. On the other hand, the above-mentioned vertical horizontal amplitude ratio calculation means 12
The amplitude ratio (VH) of vertical motion amplitude and horizontal motion amplitude is calculated by the following formula:
is calculated and the VH(t) for the past 3 seconds or so is stored. The vertical horizontal amplitude ratio calculating means 12 further calculates the past three
Calculate the average value (t) of about seconds and calculate the amplitude ratio
Seismic wave detection means 1 using VH (t) and average value (t)
Send to 3. In the case of Figure 5, when the seismic wave arrives, VH(t)
increases, and then when the S wave arrives, VH(t) decreases. In the case of Figure 7, the amplitude of the P wave is so small that it cannot be detected, and even if it is detected as an earthquake wave for the first time when the S wave arrives, it is difficult to identify that it is an S wave because VH(t) decreases. can. The seismic wave detection means 13 detects the value obtained by removing the offset level from the above-mentioned vertical sampling value.
The arrival of seismic waves is detected when x 1 (t)| exceeds seismic wave detection level A for n or more consecutive samples. The seismic wave detection level A is determined from the noise level NX 1 of x 1 (t). When a seismic wave is detected, the average value (t l ) of the amplitude ratio from the vertical and horizontal amplitude ratio calculation means 12 calculated and stored m samples before the detection and the average value (t l ) of the amplitude ratio obtained after m samples when the seismic wave is detected. Compare with the average value of the amplitude ratio (t l+2n ). As a result of the comparison, it is detected that (t l ) is larger than (t l +2n ), and the seismic wave is determined to be a P wave. After detecting the P waves, the average value of θ(t) is displayed on the display device 6 and sent to the communication control device 7 as the epicenter direction. (Effects of the Invention) According to the present invention, since the exponential smoothing value of the product of amplitudes is calculated, it is possible to accurately estimate the epicenter direction, and it is also possible to accurately determine whether it is a P wave or an S wave. .
第１図は初動振幅比から震央方位を推定する従
来の方法を説明する図、第２図は本発明の一実施
例を示すブロツク図、第３図は震央方位を説明す
る図、第４図は本発明の処理の概略フロー図、第
５図、第６図、第７図は実施例で得られた地震波
形の図である。
１……センサ、２……アンプ、１１２，１１４
……指数平滑値算出手段、１１６……震央方位算
出手段、１２……上下水平振幅比算出手段、１３
……地震波検出手段、１４……標本化手段。
Fig. 1 is a diagram explaining the conventional method of estimating the epicenter direction from the initial amplitude ratio, Fig. 2 is a block diagram showing an embodiment of the present invention, Fig. 3 is a diagram explaining the epicenter direction, and Fig. 4 is a schematic flowchart of the process of the present invention, and FIGS. 5, 6, and 7 are diagrams of seismic waveforms obtained in the example. 1...Sensor, 2...Amplifier, 112, 114
. . . Exponential smoothing value calculation means, 116 . . . Epicenter direction calculation means, 12 . . . Vertical horizontal amplitude ratio calculation means, 13
... Seismic wave detection means, 14 ... Sampling means.
Claims (1)
出するセンサと、前記センサから検出された上下
方向成分と南北方向成分との振幅値の積の指数平
滑値を算出する第１の指数平滑値算出手段と、前
記センサから検出された上下方向成分と東西方向
成分との振幅値の積の指数平滑値を算出する第２
の指数平滑値算出手段と、前記第１の指数平滑値
算出手段の出力と第２の指数平滑値算出手段の出
力とを比較して方位を算出する手段とからなる震
央方位推定装置。 ２ 上下方向、東西方向及び南北方向の地動を検
出するセンサと、前記センサから検出された上下
方向成分と南北方向成分との振幅値の積の指数平
滑値を算出する第１の指数平滑値算出手段と、前
記センサから検出された上下方向成分と東西方向
成分との振幅値の積の指数平滑値を算出する第２
の指数平滑値算出手段と、前記第１及び第２の指
数平滑値の出力を比較して方位算出する手段と、
上下方向の振幅の指数平滑値と水平方向の振幅指
数平滑値との振幅比を算出する上下水平振幅比算
出手段とを有することを特徴とする震央方位推定
装置。[Scope of Claims] 1. A sensor that detects ground motion in the vertical direction, east-west direction, and north-south direction, and a step that calculates an exponential smoothing value of the product of the amplitude values of the vertical component and the north-south component detected by the sensor. a second exponential smoothing value calculation means for calculating an exponential smoothing value of the product of the amplitude values of the vertical component and the east-west component detected by the sensor;
an epicenter azimuth estimating device comprising an exponential smoothing value calculating means, and a means for calculating an azimuth by comparing the output of the first exponential smoothing value calculating means and the output of the second exponential smoothing value calculating means. 2. A sensor that detects ground motion in the vertical direction, east-west direction, and north-south direction, and a first exponential smoothing value calculation that calculates an exponential smoothing value of the product of the amplitude values of the vertical component and the north-south component detected by the sensor. a second means for calculating an exponentially smoothed value of the product of the amplitude values of the vertical component and the east-west component detected by the sensor;
an exponentially smoothed value calculation means, and a means for comparing the outputs of the first and second exponentially smoothed values to calculate the direction;
An epicenter azimuth estimating device comprising: vertical horizontal amplitude ratio calculating means for calculating an amplitude ratio between an exponentially smoothed value of an amplitude in a vertical direction and an exponentially smoothed value of an amplitude in a horizontal direction.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP58064560A JPS59190683A (en) | 1983-04-14 | 1983-04-14 | Apparatus for estimating epicenter azimuth |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP58064560A JPS59190683A (en) | 1983-04-14 | 1983-04-14 | Apparatus for estimating epicenter azimuth |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS59190683A JPS59190683A (en) | 1984-10-29 |
| JPH0236196B2 true JPH0236196B2 (en) | 1990-08-15 |
Family
ID=13261734
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP58064560A Granted JPS59190683A (en) | 1983-04-14 | 1983-04-14 | Apparatus for estimating epicenter azimuth |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS59190683A (en) |
-
1983
- 1983-04-14 JP JP58064560A patent/JPS59190683A/en active Granted
Also Published As
| Publication number | Publication date |
|---|---|
| JPS59190683A (en) | 1984-10-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US4544892A (en) | Signal processing apparatus for frequency domain geophysical electromagnetic surveying system | |
| US6862541B2 (en) | Method and apparatus for concurrently estimating respective directions of a plurality of sound sources and for monitoring individual sound levels of respective moving sound sources | |
| KR101562904B1 (en) | Direction of Arrival Estimation Apparatus and Method therof | |
| US5341670A (en) | Method and apparatus for locating seepage from tanks | |
| CN105204065A (en) | Method and device for picking up preliminary wave | |
| CN114964259A (en) | Attitude calculation method and system based on infrared focal plane array time domain signal and infrared image | |
| US5177709A (en) | Method for determining velocity and confidence level of acoustic waves in penetrable ground | |
| US6601447B1 (en) | Acoustic anemometer for simultaneous measurement of three fluid flow vector components | |
| JP2009030990A (en) | Seismometer | |
| US4446541A (en) | Rotational geophone | |
| JP2008281462A (en) | Epicenter distance estimation apparatus, epicenter distance estimation system, and epicenter distance estimation method | |
| JPH0236196B2 (en) | ||
| CN115494546A (en) | Epicenter position detection method, device, terminal and storage medium | |
| JP2002333484A (en) | Moving target relative position detection method | |
| US20240027639A1 (en) | Seismometer with high sensitivity, broadband and all-dip | |
| JP3810323B2 (en) | Time difference direction detector | |
| JP2650935B2 (en) | Partial discharge location method | |
| JPS6319583A (en) | Device for detecting characteristics of surface layer ground | |
| JPH0263193B2 (en) | ||
| JPH0511871B2 (en) | ||
| JPH0750144B2 (en) | Partial discharge position location method | |
| JPH11183630A (en) | Method for evaluating amplification characteristic of ground at earthquake | |
| JP3106183B2 (en) | Ground exploration method and device | |
| CN206378503U (en) | A kind of device for detecting force-balanced accelerometer polarity | |
| JPH0236910B2 (en) |