JP3307840B2 - Measurement method of seismic motion with seismograph - Google Patents

Measurement method of seismic motion with seismograph

Info

Publication number
JP3307840B2
JP3307840B2 JP27668996A JP27668996A JP3307840B2 JP 3307840 B2 JP3307840 B2 JP 3307840B2 JP 27668996 A JP27668996 A JP 27668996A JP 27668996 A JP27668996 A JP 27668996A JP 3307840 B2 JP3307840 B2 JP 3307840B2
Authority
JP
Japan
Prior art keywords
value
seismic motion
acceleration
component
components
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
Application number
JP27668996A
Other languages
Japanese (ja)
Other versions
JPH10123257A (en
Inventor
善久 清水
健一 小金丸
貴 簗田
光 田久保
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.)
Tokyo Gas Co Ltd
Azbil Corp
Original Assignee
Tokyo Gas Co Ltd
Azbil Corp
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 Tokyo Gas Co Ltd, Azbil Corp filed Critical Tokyo Gas Co Ltd
Priority to JP27668996A priority Critical patent/JP3307840B2/en
Publication of JPH10123257A publication Critical patent/JPH10123257A/en
Application granted granted Critical
Publication of JP3307840B2 publication Critical patent/JP3307840B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Geophysics And Detection Of Objects (AREA)

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は制御用地震計等の地
震計における地震動の計測方法に関するものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a method for measuring seismic motion in a seismometer such as a control seismometer.

【0002】[0002]

【従来の技術】地震が発生した場合に、その強度に応じ
て各種のシステムを制御して、被害の拡大や二次災害の
発生を防止するための装置として制御用地震計があり、
制御用地震計は交通機関、都市ガス、電力、水道等の各
種施設等において、大地震時における自動緊急停止装置
に組み入れられて利用されている。
2. Description of the Related Art In the event of an earthquake, there is a control seismometer as a device for controlling various systems according to the intensity of the earthquake to prevent the spread of damage and the occurrence of secondary disasters.
2. Description of the Related Art A control seismometer is used in various facilities such as transportation, city gas, electric power, and water supply by being incorporated in an automatic emergency stop device at the time of a large earthquake.

【0003】このような制御用地震計においては構造物
の被害の程度と相関の高い制御を行うためにSI値を地
震動の強度の尺度として計測する方法が提案され、そし
て実施されている。(例えば、特開昭62−12884
号公報、特開昭62−12884号公報、特開昭62−
12884号公報、特開昭62−12885号公報、特
開昭62−12886号公報、特開平6−214040
号公報、特開平8−36062号公報等を参照のこ
と。)
In such a control seismometer, a method of measuring the SI value as a measure of the intensity of earthquake motion has been proposed and implemented in order to perform control having a high correlation with the degree of damage to the structure. (For example, Japanese Patent Application Laid-Open No. 62-12884
JP, JP-A-62-12884, JP-A-62-12884
No. 12884, Japanese Patent Application Laid-Open No. Sho 62-12885, Japanese Patent Application Laid-Open No. Sho 62-12886, Japanese Patent Application Laid-Open No. 6-214040
See JP-A-8-36062 and JP-A-8-36062. )

【0004】このような方法において、SI値は、地震
計により計測される地震動の加速度波形を1自由度振動
系の運動方程式を満たす演算部に入力して速度応答を時
々刻々と求め、速度応答の最大値のスペクトル、即ち、
速度応答スペクトルSvから所定の演算を行って求めて
いる。
In such a method, the SI value is obtained by inputting an acceleration waveform of a seismic motion measured by a seismograph into an arithmetic unit that satisfies the equation of motion of a one-degree-of-freedom vibration system, and obtaining a speed response every moment. Spectrum of the maximum value of
It is obtained by performing a predetermined calculation from the speed response spectrum Sv.

【0005】通常、地震計は直交する水平方向2成分、
例えば、南北方向と東西方向の成分の加速度を計測する
ように設置され、計測されるこれらの方向の成分の加速
度波形から夫々SI値を演算するように構成している。
Usually, a seismometer has two orthogonal horizontal components,
For example, it is installed so as to measure the acceleration of the components in the north-south direction and the east-west direction, and is configured to calculate the SI value from the measured acceleration waveforms of the components in these directions.

【0006】[0006]

【発明が解決しようとする課題】しかしながら、実際の
地震では、上述した特定の方向が必ずしも最もSI値が
大きくなるような振動をしているとは限らない。即ち、
SI値は、その性質上角度依存性が大きく、地震計によ
り計測する加速度の成分の方向により大きく値が異なる
ので、特定の方向の成分の加速度波形を単に用いて、そ
の方向のSI値を算出しても、必ずしも最大のSI値を
算出することができるとはいえず、地震による被害を評
価する際に誤差が大きくなる可能性がある。
However, in an actual earthquake, the above-mentioned specific direction does not always vibrate so as to maximize the SI value. That is,
Since the SI value has a large angle dependency due to its nature, and the value differs greatly depending on the direction of the acceleration component measured by the seismometer, the SI value in that direction is calculated simply by using the acceleration waveform of the component in a specific direction. However, it cannot always be said that the maximum SI value can be calculated, and there is a possibility that an error may increase when evaluating damage due to an earthquake.

【0007】例えば、図5は実際に起きた地震(日本海
中部地震(津軽大橋))における地震計のデータから、
全方位について算出したSI値の例を示すものである。
図において角度0゜は南北方向、90゜は東西方向に対
応するもので、これらの方向で計測される加速度波形に
よるSI値出力は43.2[kine]となる。しかしなが
ら図に示されるように、実際に最もSI値が大きくなる
のは、角度が約45゜の方向で、この方向のSI値は6
7.1[kine]であり、上記のSI値出力より約1.5
5倍も大きい。
[0007] For example, FIG. 5 shows seismograph data from an actual earthquake (Nipponkai Chubu Earthquake (Tsugaru Ohashi)).
It shows an example of SI values calculated for all directions.
In the figure, the angle 0 ° corresponds to the north-south direction and 90 ° corresponds to the east-west direction, and the SI value output by the acceleration waveform measured in these directions is 43.2 [kine]. However, as shown in the figure, the SI value actually becomes the largest in the direction at an angle of about 45 °, and the SI value in this direction is 6 °.
7.1 [kine], which is about 1.5 from the above SI value output.
Five times larger.

【0008】地震により施設に被害が生じるのは、振動
の方向には依存しないので、上述した南北方向及び東西
方向の振動のSI値出力では、被害を過小評価する恐れ
がある。
[0008] Damage to the facility due to the earthquake does not depend on the direction of the vibration. Therefore, the SI value output of the north-south direction and the east-west direction vibration described above may underestimate the damage.

【0009】ところで、地震発生後の各時点における最
大加速度の瞬時値の大きさ及びその方向は、図6に示す
ように上記南北方向及び東西方向の加速度の単純なベク
トル合成、即ち、ベクトル和を演算することにより求め
ることができるが、この最大加速度の方向は、図中x軸
方向(東西方向)とy軸方向(南北方向)の夫々の計測
値の比に依存して刻々と変化し、常時一方向に定まるも
のではない。
By the way, the magnitude and the direction of the instantaneous value of the maximum acceleration at each time point after the occurrence of the earthquake are determined by a simple vector synthesis of the north-south direction and the east-west direction acceleration as shown in FIG. The direction of the maximum acceleration changes momentarily depending on the ratio of the measured values in the x-axis direction (east-west direction) and the y-axis direction (north-south direction) in the figure, It is not always determined in one direction.

【0010】従って、このような単純なベクトル合成に
より順次求めた最大加速度から速度応答を求め、そして
これからSI値を求めても、不連続な面についての積分
を行っていることになって物理的に意味がなく、結局、
この方法でも、最大のSI値を求めることはできない。
Therefore, even if the velocity response is obtained from the maximum acceleration sequentially obtained by such a simple vector synthesis and the SI value is obtained therefrom, integration for a discontinuous surface is performed. Makes no sense, after all,
Even with this method, the maximum SI value cannot be obtained.

【0011】本発明は、以上の点に鑑みて創案されたも
のであり、その目的は、地震計により計測した地震動の
複数方向成分から、予め設定した所望の方向の地震動成
分を計測可能とすることにある。
The present invention has been made in view of the above points, and an object of the present invention is to make it possible to measure a ground motion component in a predetermined desired direction from a plurality of directional components of a ground motion measured by a seismometer. It is in.

【0012】[0012]

【課題を解決するための手段】上述した課題を解決する
ために、本発明では、まず、地震動の複数方向成分を検
出するように設置した地震計の、これらの複数方向成分
の夫々の検出信号を、予め設定した所望方向にベクトル
合成して所望方向の地震動成分を算出し、この所望方向
の地震動成分を連続的に計測することを提案する。
In order to solve the above-mentioned problems, according to the present invention, first, a detection signal of each of these multi-directional components of a seismometer installed to detect the multi-directional components of the seismic motion is provided. It is proposed to calculate a seismic motion component in a desired direction by performing vector synthesis in a desired direction set in advance , and to continuously measure the seismic motion component in the desired direction.

【0013】本発明では、上記の構成において、検出す
る複数方向成分は水平方向2成分として、所望の水平方
向の地震動成分を連続的に計測することができる。
According to the present invention, in the above-described configuration, a desired horizontal seismic motion component can be continuously measured as two components in the horizontal direction to be detected.

【0014】また本発明では、上記の構成において、検
出する複数方向成分は水平方向2成分と上下方向成分と
し、3次元の所望方向の地震動成分を連続的に計測する
ことができる。
According to the present invention, in the above configuration, the detected multi-directional components include two components in the horizontal direction and a component in the vertical direction, so that a three-dimensional earthquake motion component in a desired direction can be continuously measured.

【0015】また本発明では、上記の構成において、ベ
クトル合成により算出して連続的に計測する地震動成分
の方向は、一方向又は複数方向とすることができる。
According to the present invention, in the above configuration, the direction of the seismic motion component calculated by vector synthesis and continuously measured can be one direction or a plurality of directions.

【0016】以上の本発明によれば、地震動の複数方向
成分の夫々の検出信号をベクトル合成することにより、
地震計が検出する方向以外の所望方向の地震動成分を算
出して連続的に計測することができる。
According to the present invention described above, the detection signals of the plurality of components of the seismic motion are vector-synthesized,
Seismic motion components in desired directions other than the direction detected by the seismometer can be calculated and continuously measured.

【0017】従ってこのように連続的に計測した複数の
方向の地震動成分を比較することにより、最大の地震動
強度を有する成分を選択することができる。従って、地
震動強度の尺度をSI値とすれば、地震計の設置方向に
かかわらず、SI値の最大値を求めることができる。
Therefore, by comparing the seismic ground motion components in a plurality of directions continuously measured as described above, the component having the maximum ground motion intensity can be selected. Therefore, if the magnitude of the seismic intensity is the SI value, the maximum value of the SI value can be obtained regardless of the installation direction of the seismometer.

【0018】[0018]

【発明の実施の形態】次に本発明の実施の形態を図を参
照して説明する。図1は本発明に係る地震動の計測方法
を実施する構成の一例を概念的に示す系統図である。ま
ず図中の実線の構成を説明すると、符号1は地震計であ
り、2a,2bは直交する水平方向2成分を検出する検
出部を示すものである。これらの検出部2a,2bは、
例えば従来と同様に夫々南北方向,東西方向の成分の加
速度を計測するように設置している。
Next, an embodiment of the present invention will be described with reference to the drawings. FIG. 1 is a system diagram conceptually showing an example of a configuration for implementing a method of measuring an earthquake motion according to the present invention. First, the configuration of the solid line in the figure will be described. Reference numeral 1 denotes a seismometer, and reference numerals 2a and 2b denote detection units for detecting two orthogonal horizontal components. These detectors 2a, 2b
For example, they are installed so as to measure the acceleration of the components in the north-south direction and the east-west direction, respectively, as in the conventional case.

【0019】符号3は演算処理手段であり、この演算処
理手段3は上記検出部2a,2bから出力される南北方
向と東西方向の振動の加速度信号を後述するようにベク
トル合成して予め設定した所望方向の加速度成分を算出
する構成である。そして演算処理手段3によりベクトル
合成された所望方向の加速度成分は記憶手段4に記憶さ
れて以降の演算処理等に供される。図示の例では、ベク
トル合成により加速度成分を算出する方向は、全方位を
8等分した方向、即ち、図に示すようにθ=0゜,2
2.5゜,45゜,67.5゜,90゜,112.5
゜,135゜,157.5゜の8方向であり、符号5a
〜5hは夫々の記憶部を示している。尚、θ=0゜と9
0゜の方向は夫々検出部2a,2bによる検出方向であ
るからベクトル合成は不要である。
Reference numeral 3 denotes an arithmetic processing means. The arithmetic processing means 3 sets the acceleration signals of the north-south and east-west vibrations output from the detection sections 2a and 2b by vector synthesis as described later . This is a configuration for calculating an acceleration component in a desired direction. Then, the acceleration component in the desired direction obtained by the vector synthesis by the arithmetic processing means 3 is stored in the storage means 4 and used for the subsequent arithmetic processing and the like. In the illustrated example, the direction in which the acceleration component is calculated by vector synthesis is a direction obtained by dividing all directions into eight, that is, as shown in FIG.
2.5 ゜, 45 ゜, 67.5 ゜, 90 ゜, 112.5
{135, 157.5} in eight directions,
-5h indicate respective storage units. Note that θ = 0 ° and 9
Since the direction of 0 ° is the direction detected by the detection units 2a and 2b, vector synthesis is unnecessary.

【0020】図2は本発明に係るベクトル合成の原理
を、θ(=22.5゜)とθ′(=135゜)の方向に
ついて説明するもので、x軸は東西方向、y軸は南北方
向に対応している。図中X,Yは、夫々、ある時点にお
いて検出部2b,2aで検出した加速度を示すもので、
これらの加速度ベクトルをθ方向に射影して加えること
により、θ方向の成分の加速度を算出することができ
る。即ち、θ方向の成分の加速度をαθとすると、この
αθは、αθ=Xcosθ+Ysinθの演算式により算出さ
れる。同様に、θ′方向の成分の加速度αθは、αθ′
=Xcosθ′+Ysinθ′の演算式により算出される。こ
のように、ある時点において検出部2b,2aで検出し
た加速度に対して、予め設定した方向の夫々について以
上の演算を行うことにより、夫々の方向の成分の加速度
を算出することができ、これを適宜時間間隔で連続的に
行うことにより、所望方向、この例では水平な8方向
(演算の不要な2方向を含む)の加速度波形を連続的に
計測することができる。
FIG. 2 illustrates the principle of vector synthesis according to the present invention with respect to the directions of θ (= 22.5 °) and θ ′ (= 135 °). The x-axis is the east-west direction, and the y-axis is the north-south direction. It corresponds to the direction. In the figure, X and Y indicate accelerations detected by the detection units 2b and 2a at a certain time, respectively.
By projecting and adding these acceleration vectors in the θ direction, the acceleration of the component in the θ direction can be calculated. That is, assuming that the acceleration of the component in the θ direction is αθ, this αθ is calculated by an arithmetic expression of αθ = Xcosθ + Ysinθ. Similarly, the acceleration αθ of the component in the θ ′ direction is αθ ′
= Xcosθ '+ Ysinθ'. In this way, by performing the above-described calculation on the acceleration detected by the detection units 2b and 2a at a certain point in time in each of the preset directions, the acceleration of the component in each direction can be calculated. Is performed continuously at appropriate time intervals, thereby enabling continuous measurement of acceleration waveforms in desired directions, in this example, eight horizontal directions (including two directions requiring no calculation).

【0021】以上のようにして連続的に計測した8方向
の加速度成分は、上述したとおり、記憶手段4の各記憶
部5a〜5hを介して以降の演算処理等に供される。こ
の演算処理としては、最大加速度を求める処理や、後述
するように各方向のSI値を算出して、最大のSI値を
求める処理等が含まれる。
As described above, the acceleration components in the eight directions continuously measured as described above are supplied to the subsequent arithmetic processing and the like via the storage units 5a to 5h of the storage unit 4. This calculation processing includes processing for obtaining the maximum acceleration, processing for calculating the SI value in each direction to obtain the maximum SI value as described later, and the like.

【0022】そこで次に、後者の処理、即ち、各方向の
SI値を算出して、最大のSI値を求める処理を含めた
本発明の全体の処理の流れの一例を図3の流れ図を参照
して説明する。まずステップS1では、加速度射影角度
θを初期化、即ちθ=0とする。次いでステップS2で
は、地震計の検出部によって経時的に検出される加速度
信号の時刻を初期化、即ち、t=0として計測を開始す
る。次いでステップS3では、時刻tにおいて、各軸、
即ち、x軸,y軸の加速度を入力する。これらの加速度
信号を夫々αxt,αytとする。次いでステップS4で
は、θ方向につき上記演算式の演算を行って、時刻t,
θ方向の加速度αθtを算出する。次いでステップS5
では、ステップS4で算出した加速度αθtを、固有周
期ωを有する1自由度振動系の運動方程式を満たす演算
を行って、時刻t,θ方向、固有周期ωの速度応答Sv
θtωを算出する。次いでステップS6では、速度応答
スペクトルを計算すると共に、ステップS7において、
それまでの処理のループで算出され、最大値として記憶
されている速度応答と比較して、今回のものの方が大き
い場合には、ステップS7において最大値として記憶す
る、即ち最大値を更新すると共に、前回のものの方が大
きい場合には最大値の更新を行わずに次のステップS8
に移行する。ステップS8では、時刻を所定の時間間隔
Δtだけ進ませて、時刻t=t+ΔtとしてステップS
9に移行し、このステップS9では地震の終了の有無を
判定する。ステップS9において地震が終了していない
と判定した場合には、ステップS3に戻り、以降、時刻
t(=t+Δt)の加速度信号について、ステップS3
〜S7の処理が繰り返して行われて、速度応答スペクト
ルの最大値の更新が行われる。そしてステップS9にお
いて地震の終了の判定が成された場合には、次いで、ス
テップS10に移行して、図中の演算式により、θ方向
のSI値を算出する。このようにしてステップS10に
おいてSI値を算出した後、ステップS11に移行し、
このステップS11では、射影角度θを所定の角度幅Δ
θだけ加えて、θ=θ+Δθとして次のステップS12
に移行する。ステップS12では、最大SI値の算出を
行う。即ち、ステップS13において、今回のSI値、
即ちSIθと、それまでの処理のループで算出され、最
大値として記憶されているSImaxと比較して、今回の
ものの方が大きい場合には、ステップS13において最
大値として記憶する、即ち最大値を更新すると共に、前
回のものの方が大きい場合には最大値の更新を行わずに
次のステップS14に移行する。ステップS14では、
全ての方向についての以上の処理が完了したか否かを判
定し、θ<180゜、即ち全ての方向についての処理が
完了していないと判定した場合には、ステップS2に戻
り、ステップS11において設定された射影角度θにつ
き、以上の処理が繰り返し行われる。ステップS11で
射影角度がθ=180゜と設定されると、ステップS1
4では、全ての方向についての以上の処理が完了したと
判定して一連の処理が終了する。そして、この終了時点
においては、ステップS13の処理において記憶されて
いる最大SI値、即ち、SImaxを得ることができる。
Next, an example of the latter process, that is, the process of calculating the SI value in each direction and obtaining the maximum SI value, will be described with reference to the flowchart of FIG. I will explain. First, in step S1, the acceleration projection angle θ is initialized, that is, θ = 0. Next, in step S2, the time of the acceleration signal detected over time by the detector of the seismometer is initialized, that is, measurement is started with t = 0. Next, at step S3, at time t, each axis,
That is, the x-axis and y-axis accelerations are input. These acceleration signals are referred to as α xt and α yt , respectively. Next, in step S4, the above equation is calculated in the θ direction, and the time t,
to calculate the θ direction of the acceleration αθ t. Next, step S5
In the acceleration Arufashita t calculated in step S4, performs an operation that satisfies the equation of motion of one degree of freedom vibration system having a natural period omega, time t, theta direction, velocity response Sv natural period omega
Calculate θ t ω. Next, in step S6, a velocity response spectrum is calculated, and in step S7,
If the current response is larger than the speed response calculated in the processing loop up to that time and stored as the maximum value, the current response is stored as the maximum value in step S7, that is, the maximum value is updated. If the previous one is larger, the next step S8 is performed without updating the maximum value.
Move to In step S8, the time is advanced by a predetermined time interval Δt, and as time t = t + Δt, step S8 is performed.
Then, in step S9, it is determined whether or not the earthquake has ended. If it is determined in step S9 that the earthquake has not ended, the process returns to step S3, and thereafter, the acceleration signal at time t (= t + Δt) is returned to step S3.
Steps S7 to S7 are repeatedly performed, and the maximum value of the speed response spectrum is updated. If it is determined in step S9 that the earthquake has ended, the process proceeds to step S10, where the SI value in the θ direction is calculated by the operation formula in the figure. After calculating the SI value in step S10 in this manner, the process proceeds to step S11,
In this step S11, the projection angle θ is set to a predetermined angle width Δ
The following step S12 is performed by adding θ only and setting θ = θ + Δθ.
Move to In step S12, the maximum SI value is calculated. That is, in step S13, the current SI value,
That is, if the current value is larger than SIθ and SImax calculated in the processing loop up to that time and stored as the maximum value, the maximum value is stored in step S13, that is, the maximum value is stored. In addition to the updating, if the previous one is larger, the process moves to the next step S14 without updating the maximum value. In step S14,
It is determined whether or not the above process has been completed for all directions. If it is determined that θ <180 °, that is, if it has been determined that the process has not been completed for all directions, the process returns to step S2 and returns to step S11. The above processing is repeatedly performed for the set projection angle θ. When the projection angle is set to θ = 180 ° in step S11, step S1
In step 4, it is determined that the above processes in all directions have been completed, and the series of processes ends. Then, at this end point, the maximum SI value stored in the processing of step S13, that is, SImax, can be obtained.

【0023】以上に説明した処理においては、図1の構
成に関する上述した処理、即ち、演算処理手段3におい
て、各方向の加速度成分を予め算出した後、これらの加
速度成分につき、以降の処理を行うものとは異なり、各
方向の夫々について、加速度成分の算出に引き続いて、
以降のSI値算出等の処理を行うことが異なっており、
これらの処理の具体的手順は、この他、適宜に設定する
ことができるものである。
In the processing described above, the processing described above with reference to the configuration of FIG. 1, ie, after the acceleration components in each direction are calculated in advance by the arithmetic processing means 3, the subsequent processing is performed on these acceleration components. Unlike the ones, for each of the directions, following the calculation of the acceleration component,
The difference is that the subsequent processing such as SI value calculation is performed.
Specific procedures of these processes can be appropriately set in addition to the above.

【0024】次に本発明は、以上に説明した実施の形態
の他、例えば次のような形態をとることもできるもので
ある。 a.地震動を検出する複数方向成分は水平方向2成分に
加えて、上下方向の成分を加えて3成分とすることがで
き、この場合には3次元の所望方向の地震動成分を連続
的に計測することができる。尚、この上下方向の成分を
検出する検出部は、図1中に2点鎖線で示す検出部2c
である。 b.ベクトル合成により算出して連続的に計測する地震
動成分の方向は、上述した水平8方向とする他、その数
を適宜増減することができ、増加させれば演算処理に要
する資源及び時間が大となるが、この場合には算出した
SI値の最大値は、定義通りに求めたSI値との誤差が
小さくなって行くという利点があり、また逆に減少させ
れば、ある程度のSI値の誤差は許容した上で、SI値
の算出に要する資源及び時間を小さくすることができる
という利点がある。尚、図4は、計測する方向の数を変
化させた場合の、定義通りに求めたSI値に対する95
%信頼区間を示すものである。 c.各方向の地震動成分の算出は、単一の演算処理手段
により行うほか、複数の演算処理手段を並列処理して行
うことができる。
Next, in addition to the above-described embodiments, the present invention can take the following forms, for example. a. The multi-directional component for detecting ground motion can be made into three components by adding a vertical component in addition to two horizontal components, and in this case, it is necessary to continuously measure three-dimensional ground motion components in a desired direction. Can be. The detection unit for detecting the component in the vertical direction is a detection unit 2c indicated by a two-dot chain line in FIG.
It is. b. The direction of the seismic motion component that is calculated and continuously measured by the vector synthesis is the above-described eight horizontal directions, and the number thereof can be appropriately increased or decreased. However, in this case, the calculated maximum value of the SI value has an advantage that the error from the SI value obtained as defined decreases, and conversely, if the maximum value of the SI value decreases, the error of the SI value becomes a certain value. Has the advantage that the resources and the time required for calculating the SI value can be reduced, while permitting. FIG. 4 shows 95% of the SI value obtained as defined when the number of directions to be measured is changed.
It shows the% confidence interval. c. The calculation of the seismic motion component in each direction can be performed by a single arithmetic processing unit or by performing a plurality of arithmetic processing units in parallel.

【0025】[0025]

【発明の効果】本発明は以上のとおり、地震計により計
測した地震動の複数方向成分から、予め設定した所望の
方向の地震動成分を計測可能であるので、次のような効
果がある。 a.最小限の数の地震計の検出部により、多数の方向の
地震動を連続的に計測することができ。 b.従って方向性がある地震動強度、例えばSI値の最
大値を容易に確実に求めることができる。したがって地
震の被害を過小評価することを防止することができる。
As described above, according to the present invention, a seismic wave component in a desired direction set in advance can be measured from a plurality of directional components of a seismic wave measured by a seismometer. a. The minimum number of seismometer detectors can continuously measure ground motion in many directions. b. Therefore, it is possible to easily and reliably obtain the direction of the seismic intensity, for example, the maximum value of the SI value. Therefore, underestimation of the damage caused by the earthquake can be prevented.

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

【図1】 本発明に係る地震動の計測方法を実施する構
成の一例を概念的に示す系統図である。
FIG. 1 is a system diagram conceptually showing an example of a configuration for implementing a method for measuring a seismic motion according to the present invention.

【図2】 本発明に係るベクトル合成の原理を、θ(=
22.5゜)とθ′(=135゜)の方向について説明
するベクトル図である。
FIG. 2 shows the principle of vector synthesis according to the present invention as θ (=
FIG. 2 is a vector diagram illustrating directions of 22.5 °) and θ ′ (= 135 °).

【図3】 各方向のSI値を算出して、最大のSI値を
求める処理を含めた本発明の全体の処理の流れの一例を
示す流れ図である。
FIG. 3 is a flowchart showing an example of an overall processing flow of the present invention including processing of calculating an SI value in each direction and obtaining a maximum SI value.

【図4】 計測する方向の数を変化させた場合の、定義
通りに求めたSI値に対する95%信頼区間を示す説明
図である。
FIG. 4 is an explanatory diagram showing 95% confidence intervals for SI values obtained as defined when the number of directions to be measured is changed.

【図5】 実際に起きた地震(日本海中部地震(津軽大
橋))における地震計のデータから、全方位について算
出したSI値の例を示す説明図である。
FIG. 5 is an explanatory diagram showing an example of SI values calculated in all directions from seismograph data of an actually occurred earthquake (Chubu-Japan Sea Earthquake (Tsugaru Ohashi)).

【図6】 単純なベクトル合成を示すベクトル図であ
る。
FIG. 6 is a vector diagram showing simple vector synthesis.

【符号の説明】[Explanation of symbols]

1 地震計 2a,2b 検出部 3 演算処理手段 4 記憶手段 5a〜5h 記憶部 DESCRIPTION OF SYMBOLS 1 Seismometer 2a, 2b detection part 3 Operation processing means 4 Storage means 5a-5h storage part

フロントページの続き (72)発明者 田久保 光 東京都渋谷区渋谷2丁目12番19号 山武 ハネウエル株式会社 本社内 (56)参考文献 特開 昭62−12884(JP,A) 特開 昭62−12886(JP,A) 特開 平6−214040(JP,A) 特開 平8−36062(JP,A) 特開 昭63−121718(JP,A) 特開 平6−313729(JP,A) 特開 平2−302687(JP,A) (58)調査した分野(Int.Cl.7,DB名) G01V 1/28 G01H 1/00 G01V 1/16 Continuation of the front page (72) Inventor Hikaru Takubo 2-12-19 Shibuya, Shibuya-ku, Tokyo Yamatake Honeywell Co., Ltd. Head Office (56) References JP-A-62-128484 (JP, A) JP-A-62-1886 (JP, A) JP-A-6-214040 (JP, A) JP-A-8-36062 (JP, A) JP-A-63-121718 (JP, A) JP-A-6-313729 (JP, A) Kaihei 2-302687 (JP, A) (58) Fields investigated (Int. Cl. 7 , DB name) G01V 1/28 G01H 1/00 G01V 1/16

Claims (4)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】 地震動の複数方向成分を検出するように
設置した地震計の、これらの複数方向成分の夫々の検出
信号を、予め設定した所望方向にベクトル合成して所望
方向の地震動成分を算出し、この所望方向の地震動成分
を連続的に計測することを特徴とする地震計における地
震動の計測方法。
1. A seismometer installed so as to detect a plurality of directional components of a seismic motion, and a detection signal of each of the plurality of directional components is vector-synthesized in a predetermined desired direction to calculate a seismic motion component in a desired direction. And measuring the seismic motion component in the desired direction continuously.
【請求項2】 検出する複数方向成分は水平方向2成分
として、所望の水平方向の地震動成分を連続的に計測す
ることを特徴とする請求項1記載の地震計における地震
動の計測方法。
2. The method for measuring seismic motion in a seismometer according to claim 1, wherein the detected multi-directional components are two horizontal components, and a desired horizontal seismic motion component is continuously measured.
【請求項3】 検出する複数方向成分は水平方向2成分
と上下方向成分とし、3次元の所望方向の地震動成分を
連続的に計測することを特徴とする請求項1記載の地震
計における地震動の計測方法。
3. The seismometer according to claim 1, wherein the detected multi-directional components are two horizontal components and a vertical component, and a three-dimensional seismic motion component in a desired direction is continuously measured. Measurement method.
【請求項4】 ベクトル合成により算出して連続的に計
測する地震動成分の方向は、複数方向であることを特徴
とする請求項1〜3までのいずれか1項に記載の地震計
における地震動の計測方法。
4. The seismometer according to claim 1, wherein the direction of the seismic motion component calculated by vector synthesis and continuously measured is a plurality of directions. Measurement method.
JP27668996A 1996-10-18 1996-10-18 Measurement method of seismic motion with seismograph Expired - Lifetime JP3307840B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP27668996A JP3307840B2 (en) 1996-10-18 1996-10-18 Measurement method of seismic motion with seismograph

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP27668996A JP3307840B2 (en) 1996-10-18 1996-10-18 Measurement method of seismic motion with seismograph

Publications (2)

Publication Number Publication Date
JPH10123257A JPH10123257A (en) 1998-05-15
JP3307840B2 true JP3307840B2 (en) 2002-07-24

Family

ID=17572960

Family Applications (1)

Application Number Title Priority Date Filing Date
JP27668996A Expired - Lifetime JP3307840B2 (en) 1996-10-18 1996-10-18 Measurement method of seismic motion with seismograph

Country Status (1)

Country Link
JP (1) JP3307840B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016205920A (en) * 2015-04-20 2016-12-08 株式会社ミツトヨ Earthquake intensity value calculation system and earthquake intensity calculation method

Also Published As

Publication number Publication date
JPH10123257A (en) 1998-05-15

Similar Documents

Publication Publication Date Title
CN106382946B (en) Parameter calibrating method and device
JP5015970B2 (en) Earthquake motion prediction system
JP3307840B2 (en) Measurement method of seismic motion with seismograph
RU99100109A (en) SATELLITE RECEIVER WITH COMPUTED SPEED RELIABILITY CONTROL
JP3314912B2 (en) How to measure seismic intensity
Tichavsky et al. Quasi-fluid-mechanics-based quasi-Bayesian Crame/spl acute/r-Rao bounds for deformed towed-array direction finding
US3803916A (en) Direct measurement of earth{40 s vertical deflection using ship{40 s inertial navigation system
JP3513057B2 (en) Wide-area control support system
CN115183785A (en) Vehicle state detection method, device, terminal, program product and vehicle
JP2886862B2 (en) Loose parts monitoring device in the fluid channel of the reactor
JP2611686B2 (en) Target detection device
JPH0682067B2 (en) Earthquake detector
JPH0750144B2 (en) Partial discharge position location method
RU2107897C1 (en) Method of inertia navigation
JP2001042049A (en) Method for measuring vibration intensity
RU2154810C2 (en) Gimballess attitude control system
JPS62277573A (en) Gps position measuring instrument
JPH05272964A (en) Method and apparatus for displaying moving body posture angle
JP2001147272A (en) Early earthquake motion detecting/warning method and its device
JP2009098127A (en) Navigator
RU2098765C1 (en) Method for development of place navigational parameters and vertical
JP2001059761A (en) Gas meter with isolation valve
RU1789019C (en) Method for detection of location of miners trapped by rock fall
Gelb Geodetic and geophysical uncertainties-Fundamental limitations on terrestrial inertial navigation
RU2020417C1 (en) Inertial navigation system

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090517

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090517

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100517

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110517

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110517

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120517

Year of fee payment: 10

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130517

Year of fee payment: 11

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140517

Year of fee payment: 12

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

EXPY Cancellation because of completion of term