JP2006105862A - Method for predicting real-time earthquake risk - Google Patents
Method for predicting real-time earthquake risk Download PDFInfo
- Publication number
- JP2006105862A JP2006105862A JP2004294960A JP2004294960A JP2006105862A JP 2006105862 A JP2006105862 A JP 2006105862A JP 2004294960 A JP2004294960 A JP 2004294960A JP 2004294960 A JP2004294960 A JP 2004294960A JP 2006105862 A JP2006105862 A JP 2006105862A
- Authority
- JP
- Japan
- Prior art keywords
- time
- earthquake
- seismic
- epicenter
- point
- 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 abstract description 84
- 238000003795 desorption Methods 0.000 claims description 2
- 230000008034 disappearance Effects 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 5
- 238000007689 inspection Methods 0.000 description 13
- 230000033001 locomotion Effects 0.000 description 12
- 238000005259 measurement Methods 0.000 description 11
- 238000002360 preparation method Methods 0.000 description 9
- 238000009826 distribution Methods 0.000 description 7
- 238000011156 evaluation Methods 0.000 description 6
- 230000005484 gravity Effects 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 4
- 239000002245 particle Substances 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000002265 prevention Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000013598 vector Substances 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000007500 overflow downdraw method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000000528 statistical test Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
地震観測網が一定の水準にあって、データが即時的（リアルタイム）あるいは準リアルタイムで収集されている環境で、地震の発生後かつ、地震波主要動が到達するまでにその情報をユーザーに伝え、地震の早期警報あるいは自動防災対応措置を行うためのリアルタイム地震危険度予測の方法に関する。 In an environment where the seismic network is at a certain level and data is collected immediately (real time) or near real time, that information is communicated to the user after the earthquake occurs and before the main seismic wave arrives, The present invention relates to a method for real-time earthquake risk prediction for early earthquake warning or automatic disaster prevention measures.
全国的または地域的に整備された、一定規模の地震観測網からのデータに基づく、リアルタイム震源要素決定の方法については、代表的にユレダス法（鉄道総研報告１９巻３号、ｐｐ１〜１２）および、着未着法（特願２００１−３０８２１８号公報参照）が提案され、実証試験も行われている。前者の特徴は、観測点が一点から震源の情報が求まることであり、後者では、２ケ所以上の観測点データから相当に高い精度で震源パラメータが求まることである。 For the real-time seismic element determination method based on data from a certain-scale seismic network established nationally or regionally, the Uredas method (Railway Research Institute Vol. 19, No. 3, pp 1-12) and A non-arrival method (see Japanese Patent Application No. 2001-308218) has been proposed and a verification test has been conducted. The feature of the former is that the information of the epicenter is obtained from one observation point, and the latter is that the epicenter parameter is obtained with considerably high accuracy from the data of two or more observation points.
しかし、このようなリアルタイム地震危険度予測の方法において、特定のユーザーがその地震危険度を十分に早くかつ精確に推定し、当該の対象物の保全に供するためには、従来の方法では十分とは言えない面があり、その改善の方法が求められている。
地震の発生と共に、一定の密度をもった地震観測網からリアルタイム地震情報が発信され始める。これ等の情報を主に用いて、高信頼度が要求される緊急対応自動的制御などに供するため、特定のユーザーがその地震危険度を、より早く、十分な精度で、かつ、精度の範囲を明示して推定できるシステムが必要とされている。
However, in such a real-time earthquake risk prediction method, the conventional method is not sufficient for a specific user to estimate the earthquake risk sufficiently quickly and accurately and to provide maintenance of the target object. There are aspects that cannot be said, and there is a need for an improvement method.
As earthquakes occur, real-time earthquake information begins to be transmitted from seismic observation networks with a certain density. Because this information is mainly used for emergency response automatic control that requires high reliability, a specific user can quickly and sufficiently measure the earthquake risk within a range of accuracy. There is a need for a system that can be estimated explicitly.
ユレダス法の特徴は、検測点が一点の地震波形から震源の情報が求める。この方法では、精度がやや低いこと、震源の深さが求まらないという課題がある。また、着未着法では、迅速性に優れており、２観測点以上のデータを用いて相当に高い精度で震源パラメータが求められるが、２個以上のデータの到達を待つ必要があること、推定解の精度が不明であることが課題である。 The feature of the Yuredas method is that the information on the epicenter is obtained from the seismic waveform at one inspection point. This method has a problem that the accuracy is slightly low and the depth of the epicenter cannot be obtained. In addition, the non-arrival method is excellent in quickness, and the epicenter parameter is calculated with considerably high accuracy using data of two or more observation points, but it is necessary to wait for arrival of two or more data, The problem is that the accuracy of the estimated solution is unknown.
また、着未着法により当初それほど多くない複数の観測データを用いて震源地決定を行う場合、同時に異なる地点で発生した地震があった時には、同一の地震に属する地震波と見なされ、推定パラメータに大きな誤差が入る可能性がある。
本発明では、これらの問題を解決する方法を提示する。すなわち、リアルタイム地震危険度予測の方法の課題を解決し、実用性を向上させるため、次の２項を提供する。
１）震源情報の高精度・高速度決定の方法
２）ユーザー側でのリアルタイム高精度地震危険度推定の方法（以下「高度化着未着法」という）
In addition, when determining the epicenter using multiple observation data that are not so large at the beginning by the arrival and non-arrival method, if there are earthquakes occurring at different points at the same time, they are regarded as seismic waves belonging to the same earthquake, and the estimated parameters There may be a big error.
The present invention presents a method for solving these problems. That is, the following two items are provided to solve the problem of the real-time earthquake risk prediction method and improve the practicality.
1) Method of high-precision / high-speed determination of epicenter information 2) Real-time high-accuracy earthquake risk estimation method on the user side (hereinafter referred to as “advanced landing / non-arrival method”)
本発明では、新たな方法の提示により、以上に説明した課題に応えるものである。
この目的を達成するため、本発明においては、特許請求の範囲に記載するように構成している。すなわち、本発明は、請求項１に記載のように、
地震波観測データによる震源の時刻および位置の決定方法であって、１個以上の地震波観測データを用い、高度化着未着法によって震源要素および、その精度を推定することを特徴とするリアルタイム地震危険度予測の方法を構成する。
また、本発明は、請求項２に記載のように、
地震波観測データによる震源の時刻および位置の決定方法であって、前記１個以上の地震波観測データを用い、前記高度化着未着法とユレダス法を融合して震源を決定することを特徴とするリアルタイム地震危険度予測の方法を構成する。
また、本発明は、請求項３に記載のように、
地震波観測データによる震源の時刻および位置の決定方法であって、ほぼ同時に２カ所以上で発生した可能性のある地震を、前記高度化着未着法によって識別することを特徴とするリアルタイム地震危険度予測の方法を構成する。
また、本発明は、請求項４に記載のように、
地震波観測データによる震源の時刻および位置の決定方法を使い、特定地区での利用に供するために高精度で、特定サイトの地震危険度を推定することを特徴とするリアルタイム地震危険度予測の方を構成する。
The present invention meets the above-described problems by presenting a new method.
In order to achieve this object, the present invention is configured as described in the claims. That is, the present invention as described in
A method for determining the time and position of a hypocenter from seismic observation data, characterized by using one or more seismic observation data and estimating the source element and its accuracy by the advanced arrival and desorption method. Configure the degree prediction method.
Further, the present invention provides the following, as described in
A method for determining the time and position of an epicenter based on seismic observation data, wherein the seismic source is determined by using the one or more seismic observation data and merging the advanced arrival / disappearance method and the Uredas method. Configure real-time earthquake risk prediction methods.
Further, the present invention provides a method as claimed in
A method for determining the time and position of a hypocenter based on seismic wave observation data, characterized by identifying earthquakes that may have occurred at two or more locations at almost the same time by the advanced arrival and non-arrival method, Configure the prediction method.
Further, the present invention provides the following, as described in
Using real-time earthquake risk prediction, which uses the method of determining the time and position of the epicenter from seismic observation data, and estimates the earthquake risk at a specific site with high accuracy for use in a specific area. Constitute.
本発明に係るリアルタイム地震危険度予測の方法によって得られる効果を請求項ごとに説明する。
まず、請求項１の発明においては、高度化着未着法によって一観測点データを用いて、迅速かつ、かなりの高精度で震源時刻および位置が、その信頼限界と共に求まる。
また、請求項２の発明においては、ユレダス法と高度化着未着法による融合法によって、一つの観測データを用いて、高精度に震源の位置が、その信頼限界と共に決定できる。
また、請求項３の発明においては、許容の値を基準として、ほぼ同時に２ケ所以上で観測された地震波が、同一の地震によるものであるか否かを判定することが可能である。
また、請求項４の発明においては、特定サイトの地震危険度を、高精度に推定することが可能であり、緊急地震速報の利活用に供される。
The effect obtained by the method of real-time earthquake risk prediction according to the present invention will be described for each claim.
First, in the first aspect of the invention, the time and position of the epicenter is obtained together with its reliability limit quickly and with considerably high accuracy using one observation point data by the advanced arrival and non-arrival method.
Moreover, in the invention of
In the invention of
Further, in the invention of
以下に実施例1〜４を示す。 Examples 1-4 are shown below.
図１は、本発明に係る地震観測網の例を示す図である。図１に例示するように、日本全国に観測点数Nの観測網が一定のレベルに整備されているとする。例として高感度地震観測網（Hi-net）や気象庁の津波地震早期検知システム、多機能型地震観測網を想定する。Hi-netでは、現在約800点の観測点がある。間隔は約20kmであるが、北海道と本州の間に津軽海峡があり、そこでは測点の密度は小さいなど観測点密度は一様ではない。又、海の測点数は、現在、若干であるに過ぎない。 FIG. 1 is a diagram showing an example of an earthquake observation network according to the present invention. As illustrated in FIG. 1, it is assumed that an observation network with N observation points has been established at a certain level throughout Japan. As an example, a high-sensitivity seismic network (Hi-net), a tsunami earthquake early detection system of the Japan Meteorological Agency, and a multi-function seismic network are assumed. Hi-net currently has about 800 observation points. The distance is about 20km, but there is a Tsugaru Strait between Hokkaido and Honshu, where the station density is small and the station density is not uniform. Also, the number of sea stations is currently only small.
地震観測データは連続的に一定時間幅（Hi-netでは１秒）のパケットで伝送されてくるので、図２に例示するように、パケット単位で地震の到達した測点の到着時間が読み取られ、震源パラメータ（位置、震源時）が推定される。又、振幅から地震の規模を示すマグニチュードが推定される。 Since the seismic observation data is continuously transmitted in packets of a certain time width (1 second in Hi-net), the arrival time of the station where the earthquake arrived is read in packet units as illustrated in Fig. 2. The epicenter parameters (location, epicenter) are estimated. Also, the magnitude indicating the magnitude of the earthquake is estimated from the amplitude.
この発明では、小数点のデータから、震源パラメータを推定する新たな方法を提示する。なお、パケット毎の伝送データの場合の取り扱いは、連続伝送の一つの特殊ケースとして扱えるので、ここでは連続的に伝送されているとする。 In the present invention, a new method for estimating the epicenter parameter from the decimal point data is presented. In addition, since the handling in the case of the transmission data for every packet can be handled as one special case of continuous transmission, it is assumed here that it is continuously transmitted.
地震波が検知された時刻をｔ1、ｔ2、ｔ3、・・・ｔNとし、その測点をS1、S2、S3・・・SNとする（図２）。Nは測点の数である。S1、S2、S3・・・の位置ベクトルをｒ1、ｒ2、ｒ3・・・とする。簡単の為に観測点が水平地表面上にあると仮定し、座標系は３次元直角座標とする。水平面上にｘ、ｙ軸をとり、鉛直上方にｚ軸の正方向をとる。 Let t1, t2, t3,... TN be the times when the seismic waves were detected, and let S1, S2, S3. N is the number of stations. The position vectors of S1, S2, S3... Are r1, r2, r3. For the sake of simplicity, it is assumed that the observation point is on the horizontal ground surface, and the coordinate system is a three-dimensional rectangular coordinate. The x and y axes are taken on the horizontal plane, and the positive direction of the z axis is taken vertically upward.
震源の位置をR、震源時をｔ0とする。地震波の速度をP波に対してはVp、S波に対してはVsとする。ここでは簡単のため一様媒質と仮定し、概ねVp を６km/s、Vsを3.6km/sとする。
（１）１点検測（ｔ1≦ｔ＜ｔ2）の場合
一定の基準で当該イベントが地震によるものと判定されたとする。この判定方法そのものが一つの重要な課題であるが、本稿では取り扱わない。時間がｔ1から第２の観測点に検知される時刻ｔ2までの区間で震源のパラメータが、どのように推定されるかを調べる。
任意の時刻ｔ∈（ｔ1、ｔ2）では、震源をRとして以下の条件が成り立つ。
The location of the epicenter is R, and the time of the epicenter is t0. The velocity of the seismic wave is Vp for P waves and Vs for S waves. Here, for simplicity, a uniform medium is assumed, and Vp is approximately 6 km / s and Vs is 3.6 km / s.
(1) In the case of 1 inspection measurement (t1 ≦ t <t2) Assume that the event is determined to be due to an earthquake based on a certain standard. This judgment method itself is an important issue, but it is not dealt with in this paper. It is examined how the parameters of the epicenter are estimated in the section from time t1 to time t2 detected at the second observation point.
At an arbitrary time t∈ (t1, t2), the following conditions are satisfied with the epicenter as R.
｜R −r1｜／Vp ＝ (ｔ1 − ｔ0) (1-1)
｜R −r2｜／Vp ＞ (ｔ − ｔ0) (1-2)
・・・・・・・・・・・・・・・・
｜R −ri｜／Vp ＞ (ｔ − ｔ0) (1-i)
・・・・・・・・・・・・・・・・
式(1-1)、(1-i)から、
｜R ー rj｜ ー ｜R − r1｜ ＞ Δt・Vp (2)
２ ≦ ｉ ≦ N
となる。すなわち、測点Siと測点S1を焦点とする双曲面Ci1の、測点S1を内部に含む双曲面
| R −r1 | / Vp = (t1 − t0) (1-1)
| R −r2 | / Vp> (t − t0) (1-2)
・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・
| R −ri | / Vp> (t − t0) (1-i)
・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・
From formulas (1-1) and (1-i),
｜ R ー rj ｜ ー ｜ R − r1 ｜ ＞ Δt ・ Vp (2)
2 ≤ i ≤ N
It becomes. That is, the hyperboloid Ci1 of the hyperboloid Ci1 with the focus on the station Si and the station S1 and including the station S1 inside
の内側に震源があることになる（図３参照）。ここで座標軸を、点SiとS1を結ぶ線をXi軸とし、それに直交して鉛直線上にZ軸、さらに右手示になるようにYi軸をとる。原点Oiは、点S1と点Siを結ぶ線の中点とする。 There is an epicenter inside (see Figure 3). Here, the coordinate axis is the line connecting the points Si and S1 as the Xi axis, the Z axis is perpendicular to the vertical line, and the Yi axis is on the right hand side. The origin Oi is the midpoint of the line connecting the points S1 and Si.
この場合、式(3)のa、bは、 ｔ − ｔ1 ＝ Δtとして
｜ri − r1｜ ＝ ２ｌi1
a ＝ Vp ・ Δt (4)
In this case, a and b in the equation (3) are expressed as follows: t−t1 = Δt
｜ ri − r1 ｜ = 2li1
a = Vp · Δt (4)
となる。双曲面Ci1の当該片業がX軸と交わる点の座標は、（−ａ、０、０）である。ｔ＝ｔ1の時にはａ＝０で、Ci1は平面（X＝０）となり、震源が点S1、Siの対称面（X＝０）の点S1側の領域（Di1）にあることが言える。任意の観測点Siに関しても、同じような制限が課されるので、震源の存在可能領域はDi1（i=1,2,3,,,,）の共通領域である角柱の内部にあることになる。この場合、深さに関しては、地震活動から知られる制限を使うことになる。
Δｔが大きくなると、Ci1はより鋭い曲面となり、点S1に向かって移動すると共に収縮し、震源の存在できる領域Dilは小さくなる。
It becomes. The coordinates of the point where the one side of the hyperboloid Ci1 intersects the X axis are (−a, 0, 0). When t = t1, a = 0 and Ci1 is a plane (X = 0), and it can be said that the epicenter is in the region (Di1) on the point S1 side of the point S1 and the symmetry plane (X = 0) of Si. The same restriction is imposed on any observation point Si, so that the region where the epicenter can exist is inside the prism that is the common area of Di1 (i = 1,2,3 ,,,). Become. In this case, the limit known from seismic activity is used for depth.
As Δt increases, Ci1 becomes a sharper curved surface, moves toward the point S1 and contracts, and the region Dil where the epicenter can exist becomes smaller.
すなわち、震源位置Rの推定精度が高くなる。
Di1(t)（ｉ＝２、・・・、N）の共通領域をD1 (t)とする
That is, the estimation accuracy of the epicenter R is increased.
The common area of Di1 (t) (i = 2,..., N) is D1 (t).
地震発生域をＤと考える。Ｄは鉛直方向には、地表面Z＝０と、最大震源深度ｄmとで挟まれる領域で、水平方向は、地震発生帯で区切られるとして良い。従って全測点（i＝２、・・・、N）の検測結果を使って、震源の存在領域D1(t)は、 The earthquake occurrence area is considered as D. D is an area between the ground surface Z = 0 and the maximum epicenter depth dm in the vertical direction, and the horizontal direction may be divided by an earthquake occurrence zone. Therefore, using the inspection results of all stations (i = 2, ..., N), the existence area D1 (t) of the epicenter is
と表される、 Expressed as
D1(t)は、３次元の立体領域となる。D1を十分な精度で推定するには、点S1を囲み、最も小さな多角形をなす測点群を考えれば十分である。
震源位置の推定値を、立方領域D1(t)の重心と定める。
D1 (t) is a three-dimensional solid region. To estimate D1 with sufficient accuracy, it is sufficient to consider a station group that encloses the point S1 and forms the smallest polygon.
The estimated value of the epicenter is determined as the center of gravity of the cubic region D1 (t).
信頼限界 σ ＝ （σx、σy、σz）は、 The confidence limit σ = (σx, σy, σz) is
と定めることができる。 Can be determined.
地震活動度の空間分布S(ｒ)がわかっているとして、それに比例する重みを使って、
より合理的な推定値を定義することもできる。
Assuming that the spatial distribution S (r) of seismic activity is known, using a proportional weight,
A more reasonable estimate can also be defined.
信頼限界は、 The confidence limit is
と定義できる。
一方、読み取り誤差Δt1による震源位置推定誤差は以下のように見積られる。
Can be defined.
On the other hand, the hypocenter position estimation error due to the reading error Δt1 is estimated as follows.
いま、D1+(t) ＝ D（ｔ｜ｔ1＋Δｔ1）とすると、Δｔ1＞０の時には、 Assuming that D1 + (t) = D (t | t1 + Δt1), when Δt1> 0,
となる。このときの震源位置の推定値をR+(t)とすると、 It becomes. If the estimated value of the epicenter location at this time is R + (t),
となる。RとR＋の差の絶対値で読み取り誤差によって推定誤差ΔRを定める、 It becomes. The estimation error ΔR is determined by the reading error with the absolute value of the difference between R and R +.
となる。
（２）２点検測（ｔ2≦ｔ＜ｔ3）の場合
ｔが大きくなり２番目の観測点S2にも地震動が到達した後、３番目の点S3に到達する前までの時間帯を対象にする。
It becomes.
(2) In the case of two inspection measurements (t2 ≤ t <t3) The time period from when t becomes large and the earthquake motion reaches the second observation point S2 and before reaching the third point S3 is taken into consideration. .
この時間帯では、震源は、 During this time, the epicenter is
なる曲面C12上にある。
推定値を求める前に、ｔ＝ｔ2で測点S2にＰ波（初微動）が検出されたことの妥当性の検証を行う。点S2がないとして、点S1での検測、S3、S4・・・では未到着であるとすると、震源の存在できる領域は
Is on the curved surface C12.
Before obtaining the estimated value, the validity of the P wave (initial fine movement) detected at the measurement point S2 at t = t2 is verified. Assuming that there is no point S2, if the inspection at point S1, S3, S4 ...
であり、点S1がないとして、点S2のみの検測とした場合の震源の存在できる領域は、 Assuming that there is no point S1, the region where the epicenter can exist when only the point S2 is inspected is
となる。 It becomes.
この時には、S1またはS2のＰ波観測に誤りがあったことになる。どれが誤りであるかは不明である。 At this time, there was an error in the P-wave observation of S1 or S2. It is unclear which is wrong.
判定する方法として、地震動による粒子運動の方向と許容領域と測点との相対関係から予想される粒子運動面との整合性を調べることが考えられる。また、Ｐ波か否かについても、粒子運動の極性を使うことが出来る。なお、この粒子運動対照法の確実度を定量化しておくことも必要である。 As a determination method, it is conceivable to check the consistency of the particle motion surface predicted from the relative relationship between the direction of particle motion due to the earthquake motion and the allowable area and the measurement point. The polarity of particle motion can also be used for whether or not it is a P wave. It is also necessary to quantify the certainty of this particle motion control method.
この場合には、震源の存在可能な領域が、 In this case, the region where the epicenter can exist is
と推測される。震源位置の推定値として、領域D12の重心と定めることができる。信頼限界として、領域D12の広がりで定義できる。
（３）３点検測の場合
前項２で述べた方法を順次拡張して行くことで、３点S1、S2、S3検測の条件での、解の推定ができる。３点以上になると条検測の妥当性の判定に多数決法も導入できる。
（４）データ収集が断続的に行われる場合
地震のデータが一定の時間Δtごとのフレーム伝送によって収集される場合には、任意の時間での決定ではなく、t＝Δt、２Δt、３Δt・・・・・という時間ごとに決定することになるのみで、これまでの導き方に変更はない。
（５）簡便な計算法
実用に供するために簡易な計算の方法を呈示する。点S0を囲む４つの測点からなる四辺形S′1S′2S′3S′4の近似長方形をS1S2S3S4とする（図４参照）。図４は、平面図である。ここで近似長方形Ｓ0とＳ1のデータから、曲面C01の左葉のｘ軸負側に震源が制限される。C01を表現する式は、
It is guessed. As the estimated value of the epicenter, it can be determined as the center of gravity of the region D12. The confidence limit can be defined by the extent of the region D12.
(3) In the case of 3-inspection measurement By sequentially expanding the method described in 2 above, the solution can be estimated under the conditions of 3-point S1, S2, S3 inspection. If the score is 3 or more, the majority method can be introduced to determine the validity of the inspection.
(4) When data collection is performed intermittently When earthquake data is collected by frame transmission every fixed time Δt, t = Δt, 2Δt, 3Δt,. It will only be decided every time, and there is no change in the way it has been guided so far.
(5) Simple calculation method A simple calculation method is presented for practical use. Let S1S2S3S4 be an approximate rectangle of the quadrilateral S'1S'2S'3S'4 consisting of four stations surrounding the point S0 (see FIG. 4). FIG. 4 is a plan view. Here, from the data of the approximate rectangles S0 and S1, the epicenter is limited to the x-axis negative side of the left leaf of the curved surface C01. The expression expressing C01 is
また、C02を表現する式は、 The expression expressing C02 is
となる。ここで、 It becomes. here,
二つの双曲面C01、C02のZ＝０（地表）上の交点P＋012、P−012の座標を求める（図5参照）。 The coordinates of the intersections P + 012, P-012 on the Z = 0 (ground surface) of the two hyperboloids C01, C02 are obtained (see FIG. 5).
但し、二つの解の内
−ｌ02 ＜ X ＜ ｌ01
の方を選ぶ。式（18a）、（18b）から
However, of the two solutions
-L02 <X <l01
Choose the one. From equations (18a) and (18b)
ここで、ｌ01 ＝ ｌ02 の時には、
X ＝ ０
となる。又b＝b′であるので、双曲面C01、C02は合同となる。この時は、C01、C02で囲まれる領域の重心のX方向はS0のX座標と一致する。
Here, when l01 = l02,
X = 0
It becomes. Since b = b ′, the hyperboloids C01 and C02 are congruent. At this time, the X direction of the center of gravity of the area surrounded by C01 and C02 coincides with the X coordinate of S0.
となる。他の解は、X＜−ｌ02でこれは使わない。式(18a)にX02を入れて求まるｙの値を±Y02とする。XZ面で対応する点は（X012、−Y02）である。一般の場合には、式(19)の二つの解X1、X2（X1＜X2）の内、X2を使う。このX2を式(18a)に代入して求められるｙをY2とする。この領域C01C02重心のX座標をXとする。 It becomes. The other solutions are not used because X <-102. The value of y obtained by adding X02 to the equation (18a) is defined as ± Y02. The corresponding point on the XZ plane is (X012, -Y02). In the general case, X2 is used out of the two solutions X1 and X2 (X1 <X2) of Equation (19). Let y2 be y obtained by substituting X2 into equation (18a). Let X be the X coordinate of the center of gravity of this region C01C02.
積分は台形公式（台形数３ー４個）で求める。範囲は、 Integral is calculated by the trapezoidal formula (number of trapezoids 3-4). The range is
同じようにS3、S4に対する曲面C03C04から決まる領域C03C04の重心のY座標をInv.(Y)とする。 Similarly, the Y coordinate of the center of gravity of the region C03C04 determined from the curved surface C03C04 for S3 and S4 is Inv. (Y).
（Inv.(X)、Inv.(Y)）を領域C01C02C03C04の重心の水平方向の座標と近似できる。
（ΔZ）については、
(Inv. (X), Inv. (Y)) can be approximated to the horizontal coordinate of the center of gravity of the region C01C02C03C04.
For (ΔZ),
として、絶対値の小さい方をとる。 As a result, the smaller absolute value is taken.
小さい値を与える領域をC0iC0i＋１とすると、この間の最深点座標がY0i（i＋1）であり、Zの範囲を If the area giving a small value is C0iC0i + 1, the deepest point coordinate between them is Y0i (i + 1), and the range of Z is
とすることが出来る。Y0i(i＋1)が地震活動度深度ｈＳを越えるときは、不定扱いとする。
境界点の場合
点S0の周囲点４つの内１つまたは２つが欠けている場合である。この時には欠けた点の方向での震源の存在領域の制限に当たっては、地震活動から知られる範囲の知識を使う。
以上のように、震源位置の最大値および、分散を求めることで、一観測点データによる、震源時刻および位置が求まる（高度化着未着法）。その結果、リアルタイム地震危険度の予測の高速度・高精度化が可能となる。
It can be. Y0i (i + 1) is when more than seismicity depth h S, and indefinite treatment.
In the case of a boundary point, one or two of the four surrounding points of the point S0 are missing. At this time, in order to limit the existence area of the epicenter in the direction of the missing point, knowledge within the range known from the seismic activity is used.
As described above, by obtaining the maximum value and variance of the epicenter position, the time and position of the epicenter from one observation point data can be obtained (advanced landing method). As a result, real-time earthquake risk can be predicted at high speed and high accuracy.
ユレダス法と高度化着未着法のハイブリッド（融合）法による一つの観測データを用いた震源パラメータの推定方法を説明する。ユレダス法では、地震波の震央到達方向θと、震央距離Δが求まる。これ等の推定誤差をｄθ，ｄΔとする。これから計算される震源の位置をRｕ、確率分布をProb（Rｕ）とする。一方、高度化着未着法による震源の推定位置をR aとし、確率分布をProb（R a）とする。なお、この分布は、ユレダス法では、測定におけるデータのＳ／Ｎ（信号対雑音）比などによって決まり、少数着未着法では、主として観測点数および、場所の関数となる。
ハイブリット法では、ユレダス法と高度化着未着法との結果を重みつき平均で推定する。すなわち、
R h = λR a＋（1−λ）R u （２６）
ここで、正数λは、両方法に対する重みであり、確率分布から求められる値である。
ユレダス法と高度化着未着法との確率分布の分散を、それぞれΔｕ，ΔＡとすると、重みλは、例えば、以下のように定める。
λ＝Δｕ −１／（Δｕ −１＋ΔＡ −１） （２７）
震源時間についても、同様な方法で推定する。
以上説明したように、ユレダス法と高度化着未着法によるハイブリッド法によって、一つの観測データを用いて、高精度に震源の位置が決定できる。
The hypocenter parameter estimation method using one observation data by the hybrid (fusion) method of the Yuredas method and the advanced non-arrival method will be explained. In the Yuredas method, the epicenter arrival direction θ of the seismic wave and the epicenter distance Δ are obtained. These estimation errors are defined as dθ and dΔ. The location of the hypocenter calculated from now is Ru, and the probability distribution is Prob (Ru). On the other hand, the estimated location of the epicenter by the advanced arrival and non-arrival method is Ra, and the probability distribution is Prob (Ra). This distribution is determined by the S / N (signal-to-noise) ratio of the data in the measurement in the Uredas method, and is mainly a function of the number of observation points and the location in the non-arrival method.
In the hybrid method, the results of the Uredas method and the advanced arrival / absence method are estimated by a weighted average. That is,
R h = λR a + (1−λ) R u (26)
Here, the positive number λ is a weight for both methods and is a value obtained from the probability distribution.
If the variances of the probability distributions of the Uredas method and the advanced arrival / absence method are Δ u and Δ A , respectively, the weight λ is determined as follows, for example.
λ = Δ u -1 / (Δ u -1 + Δ A -1) (27)
The hypocenter time is estimated in the same way.
As described above, the location of the epicenter can be determined with high accuracy using a single observation data by the hybrid method based on the Uredas method and the advanced arrival and non-arrival method.
次に、同一の地震によるものであるか否かの判定方法について説明する。
高度化着未着法では、二つ以上の観測データまたは、通常の地震観測では、３点以上の地震観測であるから、その際、データが同一の地震に対応するか否かを、確かめておくことが必要となる。ここでは、その方法を提示する。
同一の地震によるものであるか否かの判定は、第１の実施形態、あるいは第２の実施形態の方法を、地震波を新たな観測点で検知する度に行う。すなわち、第一の測点データによる震源位置、震源時間を（R １，ｔ01）、第二の測点による震源位置、時間を（R ２，
ｔ0２）とする。二つのベクトルの差を
｜ΔR ｜＝｜R １−R ２｜＋Ｖｐ｜ｔ01−ｔ0２｜ （２８）
と定め、許容の値｜ΔR ｜Ｓと比較して、
｜ΔR ｜≦｜ΔR｜Ｓ；同一地震 （２９ａ）
｜ΔR｜＞｜ΔR ｜Ｓ；異なる地震 （２９ｂ）
と判定する。｜ΔR ｜Ｓは、それぞれの方法による誤差を基準として、式（２８）になら
って、決めるものとする。統計的検定（帰無仮説）を行って判定してもよい。
３点以上になったときは、多数法を占めるグループに入るか否かを判定する。
以上説明したように、許容の値｜ΔR ｜Ｓを基準として、ほぼ同時に２ケ所で検測された、同一の地震によるものであるか否かを判定することが可能である。
Next, a method for determining whether or not they are caused by the same earthquake will be described.
In the advanced arrival and non-arrival method, there are two or more observation data, or in ordinary earthquake observations, there are three or more earthquake observations. At that time, check whether the data corresponds to the same earthquake. It is necessary to keep it. Here, the method is presented.
The determination as to whether or not they are caused by the same earthquake is performed each time the seismic wave is detected at a new observation point by the method of the first embodiment or the second embodiment. That is, the location and time of the epicenter by the first station data are (R 1 , t 01 ), and the location and time of the seismic center by the second station (R 2 ,
t02 ). The difference between the two vectors is expressed as | ΔR | = | R 1 −R 2 | + Vp | t 01 −t 02 | (28)
Compared with the allowable value | ΔR | S ,
| ΔR | ≦ | ΔR | S ; Same earthquake (29a)
| ΔR |> | ΔR | S : Different earthquake (29b)
Is determined. | ΔR | S is determined according to the equation (28) based on the error of each method. The determination may be made by performing a statistical test (null hypothesis).
When the score is 3 or more, it is determined whether or not the group occupies the majority method.
As described above, it is possible to determine whether or not they are caused by the same earthquake, which was measured at two locations almost simultaneously with the allowable value | ΔR | S as a reference.
次に、リアルタイムで特定サイトの地震危険度を推定する方法を説明する。
４．１ 危険度推定方法
図６は、危険度推定方法の説明図である。図６に示す例では、地震の震央Ｓ（符号１）とユーザーＵ（符号６（Ｕ））が、ある程度離れている場合（観測網の平均的間隔をＬとして、３Ｌ〜４Ｌ）について説明する。
Next, a method for estimating the earthquake risk of a specific site in real time will be described.
4.1 Risk Estimation Method FIG. 6 is an explanatory diagram of the risk estimation method. In the example shown in FIG. 6, the case where the epicenter S of the earthquake (reference numeral 1) and the user U (reference numeral 6 (U)) are separated to some extent (3L to 4L, where L is the average interval of the observation network) will be described. .
ある防災緊急対応を行うのに必要な最小時間をＴa（allowed time）とし、準備に要する時間をＴｐ（preparatory time）とする。Ｔａはほぼ１〜数秒、Ｔｐは１０秒前後と考えてよい。というよりこのような特性を有する機器がリアルタイム地震情報の活用に適すると思われる。ユーザーは、主要動の到達前Tａに緊急制御を行い、Ｔｐより前に予備操作を開始する。予備操作は何段階かに分かれていても良い。但し、直近で地震が発生する場合には、到達前情報が、最小時間Ｔａより小さいことがあることに注意するべきである。このような場合には、ユーザーサイド（特定サイト）の地震計データを参照し緊急処置を取ることとなる。
以下、処理方法を説明する。但し全ての観測データが融合して扱える場合を例としてあげるが、ユーザーサイドの地震計データが独立に扱う必要がある時も、ほぼ同様な考え方で処理できる。
Let Ta (allowed time) be the minimum time required to perform a certain emergency disaster response, and Tp (preparatory time) be the time required for preparation. It can be considered that Ta is approximately 1 to several seconds and Tp is approximately 10 seconds. Rather, it seems that devices with such characteristics are suitable for using real-time earthquake information. The user performs emergency control before Ta reaches the main motion, and starts the preliminary operation before Tp. The preliminary operation may be divided into several stages. However, it should be noted that the pre-arrival information may be smaller than the minimum time Ta when an earthquake occurs most recently. In such a case, emergency measures are taken with reference to seismometer data on the user side (specific site).
The processing method will be described below. However, a case where all the observation data can be handled in combination is given as an example, but when seismometer data on the user side needs to be handled independently, it can be processed with almost the same concept.
図７は、危険度（震源値、時間）の算定アルゴリズムを示す図である。図７に示すように、観測データ、震源情報が逐次送達されてくるものとする。ユーザーゾーン（注意ゾーンＰＺ、動作ゾーンＲＺ）にあるデータをユーザー領域データＤ（ｕ）として、以下、図７の危険度（震源値、時間）の算定アルゴリズムを４種類（Ａ〜Ｄ）について説明する。なお、それがネットワークのものであることもあり、又ユーザーサイドのものであることもあるが、マージされているものとする。こうすることによって公用のデータと、ユーザーの固有のデータを統一的に活用できる。
Ａ．センターからｉ番目の震源情報Ｐｉが出されたとする。その確度がＰｒｏｂ（Ｐｉ）である。その時の危険度（地震強度、余裕時間）を算定する。それらをまとめてＲｉとし、確立分布をProb（Ｒｉ）とする。
Ｂ．Ｐ波がＰＺに到達する予定の時刻まで：
発信される情報のチェックを送信側のアルゴリズムと独立な算出方法で行う。すなわち、広域データＤ（ｒ）とユーザー領域データＤ（ｕ）を使って、推定震源パラメーターの妥当性を検定する。
FIG. 7 is a diagram showing an algorithm for calculating the degree of risk (seismic source value, time). As shown in FIG. 7, it is assumed that observation data and epicenter information are delivered sequentially. The data in the user zone (attention zone P Z , operation zone R Z ) is assumed to be user area data D (u), and the calculation algorithms for risk (seismic value, time) in FIG. Will be described. It may be a network or user side, but it is assumed that it has been merged. In this way, public data and user-specific data can be used in a unified manner.
A. Assume that the i-th epicenter information Pi is issued from the center. The accuracy is Prob (Pi). The risk level (earthquake intensity, margin time) at that time is calculated. These are collectively referred to as R i , and the probability distribution is referred to as Prob (R i ).
B. P-wave until the time that you plan to reach the P Z:
The transmitted information is checked by a calculation method independent of the algorithm on the transmission side. That is, the validity of the estimated hypocenter parameter is verified using the wide area data D (r) and the user area data D (u).
（１）変更が必要でないとき：
このときは、確かにまだP波がＰＺに到達していないとして良いので、準備的行動を取るか否か待機する体制となる。但し、安全サイドのリスクマネージメントでは、準備的行動を取る選択もある。
(1) When no change is required:
At this time, certainly still a P-wave is so good as not to reach the P Z, a system to wait whether or not to take the preparatory action. However, in risk management on the safety side, there is an option to take preparatory actions.
（２）変更が必要となったとき：
Ｐ波が既にＰＺに到達しているか、または当該地震以外の原因で当該測点にトリガーがかかったことになる。当該トリガーが妥当であるか否かを次の方法で判定する。
(1)近傍の測点に既に地震P波が到達しているか否か
(2)予想される地震強度か否か
(3)見かけ速度に矛盾がないか
(4)波形データが有る場合に質点の運動方向が震源方向と矛盾がないか
Ｃ．Ｐ波がPＺに入った時：
パラメータをチェック、補正する（例、(特願２００１−２５７７６５号公報参照)、発明者（藤縄等））。なお、特願２００１−２５７７６５号では、ユーザー側のデータを用いて補正するとしているが、本発明ではユーザー側と広域観測データとを統一的に扱ってより正しい補正を期す。前提として、測定条件を把握しておく必要がある。
（ａ）観測が正しく行われているか否か。
（ｂ）バックグランドノイズレベルの特性を常時調べる。
この上で、Ｐ波がＰＺに入ったことの真偽を検定し、その確度を求める。予め定めた確度の、しきい値以上で真であるとすれば、Ｐ波がＰＺに入ったことと判定する。
(2) When changes are necessary:
Or P-wave has already reached the P Z, or so that the trigger is applied to the measuring point for reasons other than the earthquake. Whether or not the trigger is valid is determined by the following method.
(1) Whether or not an earthquake P wave has already reached a nearby station
(2) Whether the earthquake intensity is expected
(3) Is there any contradiction in apparent speed?
(4) If there is waveform data, is the motion direction of the mass point consistent with the epicenter direction? When the P-wave has entered the P Z:
The parameter is checked and corrected (eg, the inventor (Fujinawa et al.) (See Japanese Patent Application No. 2001-257765). In Japanese Patent Application No. 2001-257765, correction is performed using data on the user side. However, in the present invention, the user side and the wide-area observation data are handled in a unified manner and more correct correction is expected. As a premise, it is necessary to grasp the measurement conditions.
(A) Whether the observation is performed correctly.
(B) Always examine the characteristics of the background noise level.
On this, we assayed the authenticity of the P-wave enters the P Z, determine its accuracy. Of a predetermined accuracy, if it is true at or above the threshold, it is determined that the P-wave enters the P Z.
（ｃ−１）ＰＺ内一点検測のとき：
推定P波振幅との比較：
推定値との整合性が真と判定されれば、この検測を真とする。偽となった場合には、
PＺに検測されたイベントと、当該地震は、同一でないとする。
(C-1) P Z in one inspection measurement of time:
Comparison with estimated P-wave amplitude:
If the consistency with the estimated value is determined to be true, this inspection is determined to be true. If it becomes false,
And events gage to P Z, the earthquake, and not identical.
（ｃ−２）ＰＺ内2点検測のとき：
(1)推定P波振幅との比較
(2)前述の整合性の判定
これらの判定を双方とも使って総合的に行う。（ｃ−１）の場合と同様に、イベントの同一性がないと判定された場合には、別扱いとなる。ＰＺ内の地震データから求められる危険度Ｒiが一定の閾値以上なら準備操作信号を出す。この操作は、補正震源パラメータを用いて、Ｐ波が動作ゾーンＲＺに入ると思われる時間まで行う。
(C-2) P Z in the 2 Inspection measurement of time:
(1) Comparison with estimated P-wave amplitude
(2) Consistency determination The above-mentioned determination is performed comprehensively using both of these determinations. As in the case of (c-1), when it is determined that there is no event identity, it is handled differently. Risk Ri obtained from seismic data in the P Z issues a preparation operation signal if more than a certain threshold. This operation using the correction source parameters is performed until the time you think that the P-wave enters the operating zone R Z.
（ｃ−３）別イベントの時：
同時に２つ以上の地震が発生したことになる。ＰＺ内の地震データのみを用い、先に提示した少数データによる震源パラメータの推定を行なって、危険度Ｒiが一定の閾値以上なら余裕時間を勘案して、準備操作信号単独または緊急制御も殆ど同時に出す。
(C-3) At another event:
Two or more earthquakes occurred at the same time. Using only seismic data in P Z, perform the estimation of source parameters with a small number data presented above, the risk Ri is considering the Save margin time above a certain threshold, preparation operation signal alone or the emergency control almost At the same time.
Ｄ．Ｐ波が危険領域Ｒｚ（図６の７）に入った後：
Ｒｚ内の二つ以上測点で検知され、地震危険度及び確度が一定レベル以上の時には、制御開始信号を出す。
以上、方法Ａ〜ＤをＳ波がＲｚに到達するまで繰り返す。
ユーザーサイドの地震観測網
以上の結果を用いて、ユーザー地点Ｕに当該施設があり、この施設を地震から守るための地震観測網を設計する。ユーザー地点からＲＺまでの距離をρｒ，ＰＺまでの距離を
ρｐとする。警戒時間をＴｒ，準備時間をＴｐとすると、
ρｒ＝ＶｓｘＴｒ （３０a）
ρp＝ＶｓｘＴｒ （３０b）
と定める。
大体の値は、Vs＝３．6ｋｍ／ｓ、Ｔｒ＝１ｓとして、ρｒは３．１ｋｍとなる。
D. After the P wave enters the danger zone Rz (7 in FIG. 6):
When detected at two or more stations in Rz and the seismic risk and accuracy are above a certain level, a control start signal is issued.
The methods A to D are repeated until the S wave reaches Rz.
User-side seismic observation network Using the above results, the facility is located at user point U, and an earthquake observation network is designed to protect this facility from earthquakes. The distance from the user point to R Z is ρ r , and the distance to P Z is ρ p . If the warning time is Tr and the preparation time is Tp,
ρ r = VsxTr (30a)
ρ p = VsxTr (30b)
It is determined.
Approximate value, Vs = 3.6km / s, as Tr = 1s, ρ r is the 3.1km.
そして、ＰＺに一定の地震計を配置して、公式発表データの検証を行うこととする。地震計Ｕを中心とする円形にかつ、方位角θで均等に並べるとする。正ｎ角形の頂点に配置すると、
ρｒ／ρｒ'＝ｃｏｓθ （３１）
θ＝２π／ｎ
となる。
たとえば、θ＝６０°とすると、ｎ＝６、ρｒ／ρp＝１／２となり、ρp≒７．２ｋｍとなり、ほぼ合理的な距離に設置できる。ちなみに、ｎ＝4では、ρp＝∞と、ほとんど非実用的な距離になる。
このようにすると、任意の到来方向の地震到来に対して、少なくとも２点でのデータを使用して準備操作ができる。この場合の１点検測による予備時間Ｔｐ（１）は、
√3／２ｓｅｃ≦Ｔｐ（１）≦１ｓｅｃ
また、２点検測による準備時間Ｔｐ（２）は、
０≦Ｔｐ（２）≦√3／２ｓｅｃ
の範囲となる。
さらに、Ｔｐ（２）をＴｐ'としたい場合には、
Δρｒ＝Ｔｐ'ｘＶＳ
だけ広げる。これに応じて、ρp はΔ２ρｒだけ大きくする。なお、対象をシナリオ地震に限る場合には、当該地震の方向をカバーすればよいので、その分、必要な地震計の数は少なくなる。
４．３ 地震危険度推定値の評価
推定値が妥当か否かの判定が必要になることが多い。それの基本となる考え方を説明する。図８は、地震の強度（通常は計測震度が使われるが、対象に応じて、パラメータが選択される）の推定の評価グラフを示す。図８に示す例では、計測震度を使っている。先の広いロート状の領域（Oabcd）が合格の範囲としている。被害発生最小震度Iｔ ｍ以下の震度では、被害が無いという結果になればよいので、推定値がOABCのなかに入れば合格と言うことになるが、最大無被害震度などは、対象によって異なること、被害が震度に依存することを考慮すれば、先の広いロート状の領域を合格範囲とする合理性はあるであろう。
一方、図９は、地震危険度の予測値の評価グラフのうち、地震危険度Ｒの今ひとつの指標である余裕時間に関するものである。緊急防災対応を行うのに必要な最小時間（Ｔｒ）、準備行動を取るのに必要な時間（Ｔｐ）を使って、余裕時間の合否範囲の例を示している。
Then, by placing certain of seismometers to P Z, and to perform the verification of the official announcement data. It is assumed that they are arranged in a circle centered on the seismometer U and evenly at an azimuth angle θ. When placed at the apex of a regular n-gon,
ρ r / ρ r ′ = cos θ (31)
θ = 2π / n
It becomes.
For example, if θ = 60 °, n = 6, ρ r / ρ p = 1/2, and ρ p ≈7.2 km, so that it can be installed at a substantially reasonable distance. Incidentally, when n = 4, ρ p = ∞, which is an almost impractical distance.
In this way, a preparatory operation can be performed using data at at least two points for an earthquake arrival in any direction of arrival. In this case, the preliminary time Tp (1) by one inspection measurement is:
√3 / 2sec ≦ Tp (1) ≦ 1sec
In addition, the preparation time Tp (2) by two-inspection measurement is
0 ≦ Tp (2) ≦ √3 / 2sec
It becomes the range.
Furthermore, when it is desired to set Tp (2) to Tp ′,
Δρ r = Tp′xV S
Just spread. Accordingly, ρ p is increased by Δ2ρ r . Note that when the target is limited to scenario earthquakes, it is only necessary to cover the direction of the earthquake, so the number of necessary seismometers is reduced accordingly.
4.3 Evaluation of estimated earthquake risk level It is often necessary to determine whether the estimated value is valid. Explain the basic idea. FIG. 8 shows an evaluation graph for estimating the intensity of an earthquake (usually measured seismic intensity is used, but parameters are selected depending on the object). In the example shown in FIG. 8, the measured seismic intensity is used. The wide funnel-shaped area (Oabcd) is the acceptable range. If the seismic intensity is less than the damage occurrence minimum seismic intensity I t m , it should be a result that there is no damage, so if the estimated value enters the OABC, it will be said to pass, but the maximum undamaged seismic intensity varies depending on the target Considering that the damage depends on the seismic intensity, there is a reasonable reason for accepting the wide funnel area.
On the other hand, FIG. 9 relates to a margin time which is another index of the earthquake risk R in the evaluation graph of the predicted value of the earthquake risk. The example of the pass / fail range of the surplus time is shown by using the minimum time (T r ) necessary for performing emergency disaster prevention and the time (T p ) necessary for taking a preparatory action.
正しい到達時間をＴ（ｔ）とし、推定値をＴ（e）、誤差をΔＴとすると、
Ｔ（e）＝Ｔ（ｔ）＋ΔＴ （３２）
となり、場合として3つが考えられる。
１）Ｔｒ≧Ｔ（ｔ）
この場合には、主要動が到達するまでに終えたい操作を時間内に完了させることが出来ない。この場合には、対処として、対象システムの特性にあわせて、何も行わない、あるいは、完了しないまでも、制御を開始することが選択される。
２）Ｔｐ≧Ｔ（ｔ）≧Ｔｒ：
このケースでは、準備操作には間に合わないが、緊急操作を行うには、時間の余裕が有ることになる。推定の結果が閾値以上で有れば、緊急対応操作を行う。
３）Ｔ（ｔ）≧Ｔp:
このケースでは、準備操作も行える。推定の結果がそれぞれ閾値以上で有れば、準備操作および緊急対応操作を行う。
以上説明したように、入力された震源情報をリアルタイムに処理して、特定サイトの地震危険度を推定することが可能である。
If the correct arrival time is T (t) , the estimated value is T (e) and the error is ΔT,
T (e) = T (t) + ΔT (32)
There are three possible cases.
1) T r ≧ T (t)
In this case, the operation that the user wants to complete before the main motion arrives cannot be completed in time. In this case, as a countermeasure, in accordance with the characteristics of the target system, it is selected to start control even if nothing is done or it is not completed.
2) T p ≧ T (t) ≧ T r :
In this case, it is not in time for the preparation operation, but there is time to perform the emergency operation. If the estimation result is equal to or greater than the threshold value, an emergency response operation is performed.
3) T (t) ≧ Tp:
In this case, preparation operations can also be performed. If the estimation results are each equal to or greater than the threshold value, the preparation operation and the emergency response operation are performed.
As described above, it is possible to estimate the earthquake risk of a specific site by processing the input source information in real time.
１ 震央
２〜５ 震央近傍の地震波観測点
６（ｕ） ユーザーの場所
７ 緊急指令を出す時点の主要動の波面の包絡線（Ｒｚ）
８ 準備指令を出す時点の主要動の波面の包絡線（Ｏｚ）
１０〜１５ ユーザー近傍の地震波観測点
２１ 地震波観測点（最初の地震波到達点）
２２ ある時間の地震波が最初の観測点２１に到達した時の波面の位置
２３〜３２ 地震波が未到達の観測点
1 Epicenter 2-5
8 Waveform envelope (Oz) of the main motion at the time of issuing the preparation command
10-15 Seismic observation point near the user 21 Seismic observation point (first seismic wave arrival point)
22 Position of wave front when seismic wave for a certain time reaches first observation point 21 to 32 Observation point where seismic wave has not yet reached
Claims (4)
A method for real-time earthquake risk prediction, which uses the method of determining the time and position of the epicenter based on seismic observation data and estimates the earthquake risk of a specific site with high accuracy for use in a specific area. .
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2004294960A JP4518551B2 (en) | 2004-10-07 | 2004-10-07 | Real-time earthquake risk prediction method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2004294960A JP4518551B2 (en) | 2004-10-07 | 2004-10-07 | Real-time earthquake risk prediction method |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2009293818A Division JP4868471B2 (en) | 2009-12-25 | 2009-12-25 | Real-time earthquake risk prediction method |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| JP2006105862A true JP2006105862A (en) | 2006-04-20 |
| JP2006105862A5 JP2006105862A5 (en) | 2009-09-10 |
| JP4518551B2 JP4518551B2 (en) | 2010-08-04 |
Family
ID=36375771
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2004294960A Expired - Fee Related JP4518551B2 (en) | 2004-10-07 | 2004-10-07 | Real-time earthquake risk prediction method |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4518551B2 (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009236497A (en) * | 2008-03-25 | 2009-10-15 | Railway Technical Res Inst | Quantitative evaluation method of earthquake warning information for railway |
| JP2010101903A (en) * | 2009-12-25 | 2010-05-06 | Real Time Jishin Joho Riyo Kyogikai | Method for predicting the degree of risk of an earthquake in real time |
| JP4510128B1 (en) * | 2009-05-29 | 2010-07-21 | 株式会社エイツー | Earthquake judgment system and earthquake analysis method |
| JP2010230407A (en) * | 2009-03-26 | 2010-10-14 | Taisei Corp | Earthquake motion predicting system |
| JP2012078196A (en) * | 2010-10-01 | 2012-04-19 | Taisei Corp | Seismic ground motion forecasting system |
| JPWO2020196292A1 (en) * | 2019-03-26 | 2020-10-01 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0980163A (en) * | 1995-09-18 | 1997-03-28 | Koichi Iwata | Method and apparatus for predicting earthquake |
| JPH09105781A (en) * | 1995-08-29 | 1997-04-22 | Kagaku Gijutsucho Bosai Kagaku Gijutsu Kenkyusho | Method for observing electromagnetic field related to premonitory symptoms of an earthquake and device therefor |
| JP2003066152A (en) * | 2001-08-28 | 2003-03-05 | National Research Institute For Earth Science & Disaster Provention | System for predicting and instantaneously reporting earthquake |
| JP2003114281A (en) * | 2001-10-04 | 2003-04-18 | National Research Institute For Earth Science & Disaster Provention | Determining method for focus of earthquake |
-
2004
- 2004-10-07 JP JP2004294960A patent/JP4518551B2/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH09105781A (en) * | 1995-08-29 | 1997-04-22 | Kagaku Gijutsucho Bosai Kagaku Gijutsu Kenkyusho | Method for observing electromagnetic field related to premonitory symptoms of an earthquake and device therefor |
| JPH0980163A (en) * | 1995-09-18 | 1997-03-28 | Koichi Iwata | Method and apparatus for predicting earthquake |
| JP2003066152A (en) * | 2001-08-28 | 2003-03-05 | National Research Institute For Earth Science & Disaster Provention | System for predicting and instantaneously reporting earthquake |
| JP2003114281A (en) * | 2001-10-04 | 2003-04-18 | National Research Institute For Earth Science & Disaster Provention | Determining method for focus of earthquake |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009236497A (en) * | 2008-03-25 | 2009-10-15 | Railway Technical Res Inst | Quantitative evaluation method of earthquake warning information for railway |
| JP2010230407A (en) * | 2009-03-26 | 2010-10-14 | Taisei Corp | Earthquake motion predicting system |
| JP4510128B1 (en) * | 2009-05-29 | 2010-07-21 | 株式会社エイツー | Earthquake judgment system and earthquake analysis method |
| WO2010137246A1 (en) * | 2009-05-29 | 2010-12-02 | 株式会社エイツー | Earthquake determination system and seismic analysis method |
| JP2010276536A (en) * | 2009-05-29 | 2010-12-09 | A2 Corp | Earthquake determination system and method of analyzing earthquake |
| US9366770B2 (en) | 2009-05-29 | 2016-06-14 | A2 Co., Ltd | Earthquake determination system and seismic analysis method |
| JP2010101903A (en) * | 2009-12-25 | 2010-05-06 | Real Time Jishin Joho Riyo Kyogikai | Method for predicting the degree of risk of an earthquake in real time |
| JP2012078196A (en) * | 2010-10-01 | 2012-04-19 | Taisei Corp | Seismic ground motion forecasting system |
| JPWO2020196292A1 (en) * | 2019-03-26 | 2020-10-01 | ||
| WO2020196292A1 (en) * | 2019-03-26 | 2020-10-01 | 国立研究開発法人海洋研究開発機構 | Earthquake estimation method, earthquake estimation program, and earthquake estimation device |
| US11686868B2 (en) | 2019-03-26 | 2023-06-27 | Japan Agency For Marine-Earth Science And Technology | Earthquake estimation method, non-transitory computer readable medium, and earthquake estimation device |
| JP7406828B2 (en) | 2019-03-26 | 2023-12-28 | 国立研究開発法人海洋研究開発機構 | Earthquake estimation method, earthquake estimation program, and earthquake estimation device |
Also Published As
| Publication number | Publication date |
|---|---|
| JP4518551B2 (en) | 2010-08-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Poston et al. | Indoor footstep localization from structural dynamics instrumentation | |
| CN103460073A (en) | Earthquake warning system | |
| US10598756B2 (en) | System and method for determining the source location of a firearm discharge | |
| Li et al. | An automated operational modal analysis algorithm and its application to concrete dams | |
| US20130046475A1 (en) | System and method for on-site instant seismic analysis | |
| Husen et al. | Probabilistic earthquake relocation in three-dimensional velocity models for the Yellowstone National Park region, Wyoming | |
| CN105759311A (en) | Near-real time earthquake source position positioning method | |
| US11835670B2 (en) | Seismic observation device, seismic observation method, and recording medium in which seismic observation program is recorded | |
| CN107290787A (en) | A kind of monitoring signals correlating method of earthquake infrasonic sound with the location station | |
| JP4868471B2 (en) | Real-time earthquake risk prediction method | |
| JP2006105862A (en) | Method for predicting real-time earthquake risk | |
| RU2464594C2 (en) | Method of estimating main characteristics of anticipated strong tsunamigenic earthquake and system for realising said method | |
| JP2006105862A5 (en) | ||
| CN102410525A (en) | Acoustic monitoring and location method for furnace tube leakage and device thereof | |
| KR101697227B1 (en) | Method for determining epicenter location in order to prevent false alarms of an earthquake early warning system using forced association of adjacent observation station information | |
| CN103376465B (en) | The earthquake instant analysis system of buildings floor and method thereof | |
| CN202253571U (en) | Acoustic monitoring location device for boiler tube leakage | |
| JP7015523B2 (en) | Earthquake warning system | |
| JP5770051B2 (en) | Immediate estimation method of the rupture area of a large earthquake | |
| JP2013007728A (en) | Real time estimation method of hypocentral region of giant earthquake | |
| CN109709598B (en) | False earthquake alarm identification method based on pseudo source attenuation method | |
| RU2509890C1 (en) | Control and determination method of coordinates of dangerous state of mine rock massif at underground mining operations | |
| JP2005114380A (en) | Method for real-time prediction of risk of earthquake | |
| Liu et al. | Particle swarm optimization for integrity monitoring in BDS/DR based railway train positioning | |
| Xie et al. | Integrating distributed acoustic sensing and computer vision for real-time seismic location of landslides and rockfalls along linear infrastructure |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20070612 |
|
| RD03 | Notification of appointment of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7423 Effective date: 20070612 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070622 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090723 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20091028 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20091225 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100210 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100304 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20100512 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20100517 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130528 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130528 Year of fee payment: 3 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313113 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130528 Year of fee payment: 3 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130528 Year of fee payment: 3 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| 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 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |