JP3017154B2 - Earthquake damage estimation system - Google Patents

Earthquake damage estimation system

Info

Publication number
JP3017154B2
JP3017154B2 JP10019774A JP1977498A JP3017154B2 JP 3017154 B2 JP3017154 B2 JP 3017154B2 JP 10019774 A JP10019774 A JP 10019774A JP 1977498 A JP1977498 A JP 1977498A JP 3017154 B2 JP3017154 B2 JP 3017154B2
Authority
JP
Japan
Prior art keywords
data
epicenter
temporary
earthquake
damage
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
JP10019774A
Other languages
Japanese (ja)
Other versions
JPH11211839A (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.)
NEC Corp
Original Assignee
NEC 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
Family has litigation
First worldwide family litigation filed litigation Critical https://patents.darts-ip.com/?family=12008692&utm_source=google_patent&utm_medium=platform_link&utm_campaign=public_patent_search&patent=JP3017154(B2) "Global patent litigation dataset” by Darts-ip is licensed under a Creative Commons Attribution 4.0 International License.
Application filed by NEC Corp filed Critical NEC Corp
Priority to JP10019774A priority Critical patent/JP3017154B2/en
Publication of JPH11211839A publication Critical patent/JPH11211839A/en
Application granted granted Critical
Publication of JP3017154B2 publication Critical patent/JP3017154B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Description

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

【0001】[0001]

【発明の属する技術分野】本発明は、地震による被害分
布を推定するための地震被害推定システムに関する。
The present invention relates to an earthquake damage estimating system for estimating an earthquake damage distribution.

【0002】[0002]

【従来の技術】地震被害推定システムは、震源データか
ら地域全体の被害分布を推定計算することのできるシス
テムであり、平常時において地域防災計画を立案する際
や地震発生直後において初動応急活動をする際などに、
被害状況を巨視的に把握する手段として用いられてい
る。
2. Description of the Related Art An earthquake damage estimation system is a system capable of estimating and calculating the distribution of damage in an entire area from epicenter data. At times,
It is used as a means to grasp the damage situation macroscopically.

【0003】ここで、主に平常時において地域防災計画
を立案する際に用いられる地震被害推定システムとして
は、従来、特開平8−329043号公報(以下、従来
例)に記載されているものが挙げられる。
[0003] Here, as an earthquake damage estimating system mainly used when drafting a regional disaster prevention plan in normal times, a system described in Japanese Patent Application Laid-Open No. Hei 8-329043 (hereinafter referred to as a conventional example) is conventionally known. No.

【0004】従来例の地震被害推定システムは、任意の
位置に任意の大きさで発生した地震による被害を過去の
データ等を用いてシミュレートするための地震被害の模
擬予測装置である。
[0004] The conventional earthquake damage estimation system is an earthquake damage simulation prediction device for simulating damage caused by an earthquake having an arbitrary size at an arbitrary position using past data and the like.

【0005】従来例のシステムは、データとして、地図
データ、建物データ、及び地質データを取り扱うもので
あり、これらのデータを記憶するための記憶手段と、地
図データの内容を表示するための表示手段と、震源位置
及び震源規模を設定するための地震位置設定手段とを備
えている。ここで、地図データとは、所定距離間隔で格
子状に区画されたメッシュ状地図の表示を可能とするメ
ッシュ地図データを備えており、その上、表示範囲の拡
大、縮小の選択を可能とするデータである。また、建物
データとは、建築構造物に係るデータであり、地質デー
タとは、メッシュ状地図における各区画毎の地盤の性質
を示すデータである。
The conventional system handles map data, building data, and geological data as data, and includes storage means for storing these data and display means for displaying the contents of the map data. And an earthquake position setting means for setting an epicenter and an epicenter scale. Here, the map data includes mesh map data that enables display of a mesh map partitioned into a grid at predetermined distance intervals, and further enables selection of enlargement or reduction of a display range. Data. The building data is data relating to a building structure, and the geological data is data indicating the properties of the ground for each section in the mesh map.

【0006】このような構成を備えた従来例のシステム
においては、建物データと地質データとを予め記憶手段
に登録しておき、地震位置設定手段により設定した任意
の震源データと記憶手段に格納されている地質データと
に従って、各区画ごとの揺れの強さを表す加速度を演算
すると共に、その加速度と建物データとから各建造物の
被害度を演算し、更にその演算結果をメッシュ状地図に
表示することにより、シミュレートされた地震の被害状
況を把握させることとしている。
In the conventional system having such a configuration, the building data and the geological data are registered in advance in the storage means, and are stored in the arbitrary epicenter data set by the earthquake position setting means and the storage means. According to the geological data that is being calculated, calculate the acceleration that indicates the strength of shaking for each section, calculate the degree of damage to each building from the acceleration and building data, and display the calculation result on a mesh map By doing so, it will be possible to grasp the damage situation of the simulated earthquake.

【0007】[0007]

【発明が解決しようとする課題】しかし、従来例のシス
テムは、前述の通り、平常時での利用を前提としている
システムであることから、基本的に地震発生直後の初動
応急活動に適用するには不向きな構成を備えているもの
であった。そのため、従来例のシステムを初動応急活動
に無理に適用しようとした場合、以下に示すような幾つ
かの問題が生じることとなる。
However, as described above, the conventional system is assumed to be used in normal times as described above, so that it is basically applied to the initial response immediately after an earthquake. Had an unsuitable configuration. For this reason, if the conventional system is forcibly applied to the first-response emergency activities, several problems as described below will occur.

【0008】初動応急活動が緊急を要することは、周知
の事実である。しかしながら、従来例のシステムにおい
ては、演算に先立って震源データを計算式に与えること
が演算の前提となっているため、震源データが確定する
まで演算を開始できない。即ち、従来例のシステムは、
演算開始までに、その構成から不可避な所定の時間を必
要とするという性質を有しており、結果として、初動応
急活動に迅速に応答することができないものであった。
It is a well-known fact that initial response is urgent. However, in the system of the related art, since it is premised on the calculation that the epicenter data is given to the calculation formula before the operation, the operation cannot be started until the epicenter data is determined. That is, the conventional system is
It has the property of requiring an inevitable predetermined time from the configuration before the start of the calculation, and as a result, it has not been possible to respond quickly to the initial response emergency action.

【0009】また、従来例において演算に用いられる揺
れの強さの関係式は、過去の地震データのみをもとにし
たものであり、地震発生直後に地域内で実際に観測され
た揺れの強さを反映したものではない。一方、過去の地
震データにはばらつきがあり、新たな地震における揺れ
の強さが、必ずしも過去の地震データの統計的な平均値
に一致するとは限らない。従って、従来例のシステムに
おいては、信頼性の高い演算結果を出力することが困難
であった。
Further, the relational expression of the magnitude of the shaking used in the calculation in the conventional example is based on only the past earthquake data, and the magnitude of the shaking actually observed in the area immediately after the earthquake occurs. It does not reflect that. On the other hand, past earthquake data varies, and the intensity of shaking in a new earthquake does not always match the statistical average value of past earthquake data. Therefore, in the conventional system, it was difficult to output a highly reliable operation result.

【0010】更に、従来例のシステムは、地表近くの比
較的浅い断層面で発生した地震について、演算できない
という問題点を有していた。一般に、此の種の地震にお
いては、距離的に幅のある断層面全体が震源となること
が知られている。しかしながら、従来例のシステムは、
1点の震源データを唯一の入力パラメータとしている。
従って、従来例のシステムでは、此の種の地震のよう
に、震源が距離的に幅を有する場合、所望とする精度の
出力結果を得ることができなかった。
[0010] Further, the conventional system has a problem that it is impossible to calculate an earthquake generated on a relatively shallow fault plane near the surface of the ground. Generally, it is known that in this kind of earthquake, the entire fault plane having a wide distance becomes the epicenter. However, the conventional system is
One source is the only input parameter.
Therefore, in the system of the conventional example, when the epicenter has a width in terms of distance as in this type of earthquake, it is not possible to obtain an output result with a desired accuracy.

【0011】そこで、本発明は、初動応急活動に適用で
きる地震被害推定システムであって、信頼性の高い演算
結果を得ることが出来ると共に、地表近くの比較的浅い
断層面で発生した地震についても所望とする精度の演算
結果を出力することのできる地震被害推定システムを提
供することを目的とする。
Therefore, the present invention is an earthquake damage estimating system applicable to initial response emergency activities, which can obtain a highly reliable calculation result and can also be applied to an earthquake generated on a relatively shallow fault plane near the ground surface. An object of the present invention is to provide an earthquake damage estimation system capable of outputting a calculation result with desired accuracy.

【0012】また、本発明は、上記システムに適用され
る地震被害推定手法、及び演算処理装置と共に当該地震
被害推定システムを構成し得る記録媒体を提供すること
を目的とする。
It is another object of the present invention to provide a seismic damage estimation method applied to the above-mentioned system, and a recording medium capable of constituting the earthquake damage estimation system together with an arithmetic processing unit.

【0013】[0013]

【課題を解決するための手段】本発明は、上述した課題
を解決すべく、緊急を要する地震発生直後の初動応急活
動において、地域内の複数の震度計で実際に観測された
揺れの強さ(計測震度、最大加速度、最大速度などの総
称)をもとにして、地域全体の被害分布を高速かつ高精
度に演算することにより、被害状況を巨視的に把握でき
る手段を提供する。
SUMMARY OF THE INVENTION In order to solve the above-mentioned problems, the present invention provides an initial response immediately after the occurrence of an urgent earthquake, and the magnitude of the shaking actually observed by a plurality of seismometers in the area. Based on (measured seismic intensity, maximum acceleration, maximum speed, etc.), it provides a means to macroscopically grasp the damage situation by calculating the damage distribution of the whole area at high speed and with high accuracy.

【0014】具体的には、本発明による解決手段は、以
下に示される。
Specifically, the solution according to the present invention is shown below.

【0015】即ち、本発明によれば、第1の地震被害推
定システムとして、地震を観測する複数の観測地点で、
夫々、得られた複数の実測データを入力とし、推定対象
地域における地震被害を推定するための地震被害推定シ
ステムにおいて、前記複数の実測データを入力するため
の入力手段と、震央位置からの地表上における距離であ
る震央距離と地震の揺れの強さとの関係を示す関係式デ
ータを含み、地震被害推定の演算に用いられる基礎情報
を予め記憶する記憶手段と、前記複数の実測データに基
づいて震源データを算出すると共に、当該算出結果を前
記関係式データに与えてメッシュ状の被害分布を演算す
るための演算処理装置と、当該演算処理装置からの演算
結果を受けて、地震被害推定結果として表示するための
表示手段とを備えることを特徴とする地震被害推定シス
テムが得られる。
That is, according to the present invention, as a first earthquake damage estimation system, at a plurality of observation points for observing an earthquake,
An earthquake damage estimation system for estimating seismic damage in the estimation target area, each of the plurality of obtained measured data being input, and an input means for inputting the plurality of measured data; Storage means for storing in advance the basic information used for the calculation of the earthquake damage estimation, including relational expression data indicating the relationship between the epicenter distance, which is the distance at, and the intensity of the earthquake tremor, and a hypocenter based on the plurality of measured data. An arithmetic processing device for calculating the data and applying the calculation result to the relational expression data to calculate a mesh-like damage distribution, and receiving the calculation result from the arithmetic processing device and displaying the result as an earthquake damage estimation result And a display means for performing the operation.

【0016】また、本発明によれば、第2の地震被害推
定システムとして、前記第1の地震被害推定システムに
おいて、前記演算処理装置は、前記震源データを算出す
る際に、複数の仮の震源を想定し、前記複数の仮の震源
の夫々に係る複数の仮震源データに従って、当該仮震源
データを前記関係式データに代入して求められた理論曲
線と揺れの強さに係る実測データとの誤差の総和を算出
する処理を、前記複数の仮震源データ複数の前記震源デ
ータの夫々について実行し、その算出結果から前記誤差
の総和を極小とする一の前記仮震源データを震源に係る
震源データとして推定決定するものであることを特徴と
する地震被害推定システムが得られる。
According to the present invention, as the second earthquake damage estimation system, in the first earthquake damage estimation system, the arithmetic processing unit may calculate a plurality of temporary earthquake sources when calculating the earthquake source data. Assuming that, according to a plurality of temporary hypocenter data according to each of the plurality of temporary hypocenters, the theoretical curve obtained by substituting the temporary hypocenter data into the relational expression data and the measured data related to the magnitude of shaking The process of calculating the sum of the errors is performed for each of the plurality of the hypocenter data and the plurality of the hypocenter data, and the one of the hypocenter data that minimizes the sum of the errors is calculated from the calculation result. The seismic damage estimation system is characterized in that it is determined by estimation.

【0017】また、本発明によれば、第3の地震被害推
定システムとして、前記いずれかの地震被害推定システ
ムにおいて、前記演算処理装置は、前記震央位置を地表
上における線状の位置座標として定義すると共に、前記
震央距離を当該線状の震央位置からの前記地表上におけ
る最短距離として算出し、当該算出結果を前記関係式デ
ータに代入することにより、前記被害分布を演算するも
のであることを特徴とする地震被害推定システムが得ら
れる。
According to the present invention, as a third earthquake damage estimation system, in any one of the above earthquake damage estimation systems, the arithmetic processing unit may define the epicenter position as a linear position coordinate on the ground surface. And calculating the epicenter distance as the shortest distance on the ground from the linear epicenter position, and substituting the calculation result into the relational expression data to calculate the damage distribution. The characteristic earthquake damage estimation system is obtained.

【0018】更に、本発明によれば、第4の地震被害推
定システムとして、前記第1乃至第3のいずれかの地震
被害推定システムにおいて、前記表示手段は、前記地震
被害推定結果を、視覚的に把握しやすくするために、グ
ラフ表示するものであることを特徴とする地震被害推定
システムが得られる。
Further, according to the present invention, as a fourth earthquake damage estimation system, in any one of the first to third earthquake damage estimation systems, the display means displays the earthquake damage estimation result visually. A seismic damage estimation system characterized in that it is displayed in a graph so that it can be easily grasped.

【0019】以上開示したように、本発明の地震被害推
定システムにおいては、まず、任意、且つ、複数の仮震
源を想定し、各観測地域内における実測データに基づい
て、各仮震源の夫々における仮震源データ(震央位置
と、震源深さと、地震規模を表すマグニチュード)を求
め、当該複数の仮震源データの中から、推定される震源
に係る震源データを決定し、次いで、その結果を震央距
離と揺れの強さの関係式に与えることにより、メッシュ
状の被害分布について単純計算にて高速に求めることを
可能とする。ここで、本発明の地震被害推定システムに
おいては、演算に用いられる関係式として、過去の地震
データをもとにして曲線の勾配に理論的根拠を持たせた
一般的な相関式を利用する。
As disclosed above, in the earthquake damage estimating system of the present invention, first, an arbitrary and plural temporary hypocenters are assumed, and each of the temporary hypocenters is determined based on measured data in each observation area. Obtain the hypocenter data (the epicenter location, epicenter depth, and magnitude indicating the magnitude of the earthquake), determine the epicenter data for the estimated epicenter from the multiple epicenter data, and then use the result as the epicenter distance And the magnitude of the swaying strength, it is possible to obtain the mesh-like damage distribution at high speed by simple calculation. Here, in the earthquake damage estimation system of the present invention, a general correlation equation in which the slope of a curve has a theoretical basis based on past earthquake data is used as a relational equation used for calculation.

【0020】より詳細には、この相関式に対して仮震源
データを与えることにより得られる理論曲線と、各観測
地域内で観測された実測データとの誤差の総和を各仮震
源データ毎に求め、当該誤差の総和が極小となるときの
仮震源データを、震源データ、即ち震源の推定値と判断
する。このような判断方法に従って得られる理論曲線
は、実測データのみを数学的に処理した単なる近似曲線
とは異なり、その勾配に理論的根拠を有するものである
一方、実測データで最適化されていることから実際に観
測された揺れの強さを反映しているものであると言え
る。従って、本発明の地震被害推定システムにおいて
は、実測データに近似しているにもかかわらず、測定誤
差の影響を受けにくい、高精度な被害推定が可能とな
る。
More specifically, the sum of the error between the theoretical curve obtained by giving the temporary hypocenter data to this correlation equation and the actual measurement data observed in each observation area is obtained for each temporary hypocenter data. The temporary epicenter data when the sum of the errors is minimized is determined as the epicenter data, that is, the estimated value of the epicenter. The theoretical curve obtained according to such a judgment method is different from a mere approximation curve obtained by mathematically processing only the actual measurement data, and has a theoretical basis for its gradient, while being optimized with the actual measurement data. Therefore, it can be said that this reflects the actually observed shaking strength. Therefore, in the earthquake damage estimation system of the present invention, it is possible to perform damage estimation with high accuracy, which is hardly affected by measurement errors, despite being close to actual measurement data.

【0021】更に、震源データ(及び仮震源データ)の
一要素である震央位置について、本発明では、1点では
なく線状に幅を持たせることとした。その結果、本発明
の地震被害推定システムは、点状の震源では想定が不可
能であった地表近くの断層面に沿って発生した地震にも
適応でき、理論曲線の精度を向上させることができる。
Further, in the present invention, the epicenter position, which is one element of the epicenter data (and the temporary epicenter data), is given a linear width instead of a single point. As a result, the earthquake damage estimation system of the present invention can be applied to an earthquake that occurred along a fault plane near the surface of the ground, which could not be estimated from a point-like source, and the accuracy of the theoretical curve can be improved. .

【0022】また、本発明においては、地震被害の推定
結果を直感的に把握しやすいようにグラフ化して表示す
ることとした。そのため、本発明によれば、求めた理論
曲線と実測データとの誤差の広がりを視覚的に確認し
て、推定結果の信頼性の度合いを直感的に把握すること
が可能となる。
Further, in the present invention, the estimation result of the earthquake damage is displayed in a graph so as to be easily intuitively grasped. Therefore, according to the present invention, it is possible to visually confirm the spread of the error between the obtained theoretical curve and the actually measured data, and to intuitively grasp the degree of reliability of the estimation result.

【0023】尚、本発明の概念は、当該推定方法を情報
処理装置に実行させるプログラムを格納した記録媒体と
して具現化され得るものである。
Note that the concept of the present invention can be embodied as a recording medium storing a program for causing an information processing apparatus to execute the estimation method.

【0024】即ち、本発明によれば、データを格納する
記憶手段を備えると共にプログラムに従う処理を実行す
る情報処理装置に対して、地震被害の推定対象地域内に
おける複数の観測地点において実際に観測された結果と
しての実測データに基づいて、前記推定対象地域におけ
る被害分布を推定させるために、複数の仮の震源を想定
し、前記複数の実測データの夫々に基づいて、前記複数
の震源に係る複数の仮震源データであって、夫々に震央
位置及び地震規模を示すマグニチュードを内包する仮震
源データを求める第1の処理と、震央距離と揺れの強さ
を理論的に示す関係式として、過去の地震データに基づ
いて、曲線の勾配に理論的な根拠を持たせてなる一般的
な相関式を用い、当該相関式に対して、前記第1の処理
において求めた前記仮震源データの夫々を代入すること
により、夫々の仮震源データに従った理論曲線を算出さ
せる第2の処理と、前記第2の処理の結果である各理論
曲線と前記複数の実測データとの誤差の総和を算出する
第3の処理と、前記複数の仮震源データの内、前記第3
の処理により算出された前記誤差の総和を極小とする理
論曲線に応じた前記仮震源データを前記震源に係る震源
データであるとして推定し、メッシュ上の被害分布につ
いて推定する第4の処理とを前記情報処理装置に実行さ
せる命令を含むプログラムを格納した記録媒体であっ
て、前記情報処理装置の読取可能な記録媒体が得られ
る。
That is, according to the present invention, an information processing apparatus having a storage means for storing data and executing a process according to a program is actually observed at a plurality of observation points in a target area for estimation of earthquake damage. Based on the actual measurement data as a result, in order to estimate the damage distribution in the estimation target area, a plurality of temporary hypocenters are assumed, and based on each of the plurality of actual measurement data, a plurality of hypocenters relating to the plurality of epicenters are estimated. The first processing for obtaining temporary hypocenter data including the magnitude indicating the epicenter position and the magnitude of the epicenter, respectively, and the relational expression theoretically indicating the epicenter distance and the magnitude of the quake, Based on the seismic data, a general correlation formula is used in which the slope of the curve has a theoretical basis. By substituting each of the temporary hypocenter data, a second process of calculating a theoretical curve according to each of the temporary hypocenter data, and a process of each of the theoretical curves, which are the result of the second process, and the plurality of measured data are performed. A third process for calculating a sum of errors; and a third process among the plurality of temporary hypocenter data.
And estimating the provisional hypocenter data according to the theoretical curve that minimizes the sum of the errors calculated by the above processing as the epicenter data relating to the epicenter, and estimating the damage distribution on the mesh. A recording medium storing a program including an instruction to be executed by the information processing device, the recording medium being readable by the information processing device is obtained.

【0025】ここで、前記記録媒体で取り扱われる前記
震源データは、その一要素として、震源の地表上におけ
る位置を示す震央位置を含み、且つ、当該震央位置は、
線状で示されるものであるものであっても良いことは言
うまでもない。
Here, the epicenter data handled by the recording medium includes, as one element thereof, an epicenter position indicating the position of the epicenter on the ground surface, and the epicenter position is:
Needless to say, it may be a linear one.

【0026】[0026]

【発明の実施の形態】本発明の実施の形態による地震被
害推定システムについて説明するにあたって、まず、当
該説明において用いられる用語について、図1及び図2
を用いて説明する。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS In describing an earthquake damage estimation system according to an embodiment of the present invention, terms used in the description will first be described with reference to FIGS.
This will be described with reference to FIG.

【0027】図1は、震央(震央位置)、震源深さ、震
央距離、観測地点、地表、表層地質、硬質基盤などの語
句の定義及びそれらの関係を示す。図1を参照すれば理
解されるように、地表に対して垂直に交わり且つ震源を
通る線と、地表との交点を震央(若しくは、震央位置)
と言い、震央位置から震源までの深さを震源深さと言
う。また、地表上において地震を観測する地点を観測地
点と言い、震央位置から観測地点までの距離を震央距離
と言う。更に、地表近傍の比較的浅い断層面における地
質を表層地質と言い、表層地質より深い位置に位置する
硬質な地質を硬質基盤と言う。
FIG. 1 shows definitions of terms such as epicenter (center of epicenter), epicenter depth, epicenter distance, observation point, ground surface, surface geology, and hard base, and their relationships. As can be understood with reference to FIG. 1, the intersection of the line perpendicular to the ground surface and passing through the epicenter and the ground surface is the epicenter (or epicenter position).
The depth from the epicenter to the epicenter is called the epicenter depth. The point on the surface of the earth where the earthquake is observed is called the observation point, and the distance from the epicenter to the observation point is called the epicenter distance. Furthermore, the geology at a relatively shallow fault plane near the surface is called surface geology, and the hard geology located deeper than the surface geology is called hard basement.

【0028】図2は、線状の震源を考慮した場合におけ
る震央位置及び震央距離を示す。前述のように、地表近
くの比較的浅い断層面で生じた地震の場合、距離的に幅
のある断層面全体が震源となることから、震央位置を点
状と仮定したのでは従来例と同様の問題点を有すること
となり、また、震央位置を点状でない形状に仮定した場
合であっても、単純に震央距離を震央位置から観測地点
までの距離と定義したのでは、震央距離の解釈に疑義の
生ずる恐れがある。そこで、一例として、線状の震源の
場合について、震央位置を「線状の震源を地表まで垂直
に移動した線分」と定義し、震央距離を「当該震央位置
(線分)までの最短距離」と定義することとする。尚、
他の形状の震源の場合についても、震源の形状に伴って
震央位置の形態を仮定することにより、震央距離につい
ては、同様の定義を適用し得ることは、言うまでもな
い。
FIG. 2 shows the epicenter position and epicenter distance when a linear epicenter is considered. As described above, in the case of an earthquake that occurred on a relatively shallow fault plane near the ground surface, the entire fault plane with a wide distance becomes the epicenter. Even if the epicenter position is assumed to be a non-point-like shape, simply defining the epicenter distance as the distance from the epicenter position to the observation point will not allow the interpretation of the epicenter distance. There is a risk of doubt. Therefore, as an example, in the case of a linear epicenter, the epicenter position is defined as "the line segment that has moved the linear epicenter vertically to the surface of the ground", and the epicenter distance is "the shortest distance to the epicenter position (line segment). ". still,
It is needless to say that the same definition can be applied to the epicenter distance for hypocenters of other shapes by assuming the form of the epicenter position according to the shape of the epicenter.

【0029】このような定義の下、本発明の実施の形態
による地震被害推定システムについて、以下に、図面を
参照して詳細に説明する。
Under such a definition, an earthquake damage estimation system according to an embodiment of the present invention will be described in detail below with reference to the drawings.

【0030】本発明の実施の形態による地震被害推定シ
ステムは、図3に示されるように、複数地点の揺れの強
さを観測する地震観測装置1と、演算に必要な基礎情報
などのデータを格納する記憶装置2と、プログラム制御
により動作する演算処理装置3と、演算結果を出力する
表示装置4とを備えている。
As shown in FIG. 3, the earthquake damage estimating system according to the embodiment of the present invention includes an earthquake observation device 1 for observing the intensity of shaking at a plurality of points and data such as basic information necessary for calculation. The apparatus includes a storage device 2 for storing, an arithmetic processing device 3 operated by program control, and a display device 4 for outputting an arithmetic result.

【0031】地震観測装置1は、地域内の複数地点に設
置された震度計11と、観測されたデータを集計する受
信処理装置12と、それらを結ぶ伝送路を備えており、
演算処理装置3に対して、揺れの強さに関する実測デー
タを伝達する。ここで、実測データとは、地震による揺
れの強さを表す物理量で、計測震度、最大加速度、最大
速度のいずれか一つを地表で観測した値である。
The earthquake observation device 1 includes seismometers 11 installed at a plurality of points in the area, a reception processing device 12 for summing up observed data, and a transmission line connecting them.
The measured data relating to the intensity of shaking is transmitted to the arithmetic processing unit 3. Here, the actual measurement data is a physical quantity indicating the strength of the shaking caused by the earthquake, and is a value obtained by observing one of the measured seismic intensity, the maximum acceleration, and the maximum speed on the ground surface.

【0032】記憶装置2は、演算に必要な基礎情報とし
て、関係式データ21と、地域特性データ22とを予め
記憶している。
The storage device 2 stores relational expression data 21 and regional characteristic data 22 in advance as basic information necessary for calculation.

【0033】演算処理装置3は、震源推定手段31と、
被害推定手段33とを備え、これら手段による演算結果
としての震源データ32及び被害推定データ34を記憶
保持し得る機能を有する。震源推定手段31は、地震観
測装置1から入力された揺れの強さの実測データと、記
憶装置2に記憶されている関係式データ21及び地域特
性データ22を用いて、震源データ32を算出する。被
害推定手段33は、震源推定手段31により算出された
震源データ32と、記憶装置2からの関係式データ21
及び地域特性データ22とを用いて、メッシュ状の被害
推定データ34を算出する。
The arithmetic processing unit 3 includes a hypocenter estimation unit 31 and
It has a damage estimation means 33, and has a function of storing and holding the epicenter data 32 and the damage estimation data 34 as the calculation results by these means. The hypocenter estimation means 31 calculates hypocenter data 32 using the actually measured shaking intensity data input from the earthquake observation apparatus 1 and the relational expression data 21 and the regional characteristic data 22 stored in the storage device 2. . The damage estimating means 33 includes the epicenter data 32 calculated by the epicenter estimating means 31 and the relational expression data 21 from the storage device 2.
Then, the mesh-like damage estimation data 34 is calculated using the data and the regional characteristic data 22.

【0034】表示装置4は、震源推定手段31からの震
源データ32、及び被害推定手段33からの被害推定デ
ータ34を受けて、被害状況を直感的に把握しやすいグ
ラフとして、表示する。尚、表示装置4を単なるディス
プレイとし、震源データ32及び被害推定データ34を
受けて当該表示装置4に前述のグラフを表示させるよう
に制御するグラフィックコントローラを、別途、有して
いることとしても良い。
The display device 4 receives the epicenter data 32 from the epicenter estimation unit 31 and the damage estimation data 34 from the damage estimation unit 33, and displays the damage status as a graph that is easy to intuitively grasp. The display device 4 may be a mere display, and may further include a graphic controller that receives the epicenter data 32 and the damage estimation data 34 and controls the display device 4 to display the above-described graph. .

【0035】このような構成を備えた本実施の形態によ
る地震被害推定システムにおいては、概略、次のような
処理が行なわれ、地震の被害推定に係るデータが表示装
置4に表示されることとなる。即ち、地震が発生する
と、震度計11は、当該地震による地震波を観測し、揺
れの強さに関する実測データを出力する。受信処理装置
12は、震度計11からの実測データを収集して、演算
処理装置3に伝達する。これらのデータが受信処理装置
12から演算処理装置3に伝達されると、震源推定手段
31は、実測データ、並びに関係式データ21及び地域
特性データ22を用いて、直ちに、震源データ32を算
出する。続いて、被害推定手段33は、震源データ3
2、並びに関係式データ21及び地域特性データ22か
ら、メッシュ状の被害推定データ34を算出すると同時
に表示装置4に出力する。表示装置4は、震源データ3
2及び被害推定データ34から作成された被害状況を示
すグラフを表示する。
In the earthquake damage estimating system according to the present embodiment having such a configuration, the following processing is performed, and data relating to earthquake damage estimation is displayed on the display device 4. Become. That is, when an earthquake occurs, the seismic intensity meter 11 observes the seismic wave caused by the earthquake and outputs actual measurement data regarding the strength of the shaking. The reception processing device 12 collects measured data from the seismic intensity meter 11 and transmits it to the arithmetic processing device 3. When these data are transmitted from the reception processing device 12 to the arithmetic processing device 3, the epicenter estimating means 31 immediately calculates the epicenter data 32 using the measured data, the relational expression data 21 and the regional characteristic data 22. . Subsequently, the damage estimation means 33 calculates the epicenter data 3
2, and calculates the mesh-like damage estimation data 34 from the relational expression data 21 and the regional characteristic data 22 and outputs it to the display device 4 at the same time. The display device 4 displays the epicenter data 3
2 and a graph showing the damage situation created from the damage estimation data 34 is displayed.

【0036】次に、震源推定手段31における動作につ
いて、図4に示されるフローチャートを用いて、詳細に
説明する。
Next, the operation of the epicenter estimating means 31 will be described in detail with reference to the flowchart shown in FIG.

【0037】各観測地域内に夫々設けられた震度計11
において観測された実測データは、受信処理装置12を
介して、演算処理装置3に入力され、震源推定手段31
に伝達される(ステップS401)。
The seismic intensity meter 11 provided in each observation area
The measured data observed in the above is input to the arithmetic processing unit 3 via the reception processing unit 12, and the hypocenter estimation unit 31
(Step S401).

【0038】震源推定手段31に実測データが入力され
ると、震源推定手段31は、硬質基盤における揺れの強
さを算出するために、まず、記憶装置に格納されている
地域特性データ22を用いて、当該実測データを硬質基
盤における値に変換する(ステップS402)。
When the measured data is input to the epicenter estimation unit 31, the epicenter estimation unit 31 first uses the regional characteristic data 22 stored in the storage device to calculate the strength of the shaking on the hard base. Then, the actual measurement data is converted into a value on the hard base (step S402).

【0039】硬質基盤における揺れの強さは、理想的に
は地域特性を含まない一般化された物理量であり、以下
に、一例として、震央距離と基盤最大速度の関係式vb
=f(M,d,r)を適用して、算出する場合について
説明する。尚、dは震源深さを示し、Mは地震規模を表
すマグニチュードを示す。
The strength of shaking on a rigid base is ideally a generalized physical quantity that does not include regional characteristics. For example, the relational expression v b between the epicenter distance and the base maximum velocity will be described below.
= F (M, d, r) will be described. Note that d indicates the epicenter depth, and M indicates the magnitude indicating the magnitude of the earthquake.

【0040】実測データとして地表最大速度v(n)を
用いた場合、当該地表最大速度v(n)を表層地質の増
幅率w(n)で割ることにより、基盤最大速度v
b(n)[=v(n)÷w(n)]を算出することがで
きる。ここで、増幅率w(n)は、表層地質の種類によ
って異なるものであり、経緯度座標で一意に決まる地域
特性データ22の一つである。
When the maximum surface velocity v (n) is used as the measured data, the maximum surface velocity v (n) is divided by the amplification factor w (n) of the surface geology to obtain the maximum base velocity v (n).
b (n) [= v (n) ÷ w (n)] can be calculated. Here, the amplification factor w (n) differs depending on the type of surface geology, and is one of the regional characteristic data 22 uniquely determined by the latitude and longitude coordinates.

【0041】次いで、演算を行なうにあたって、震源デ
ータ32に仮定値を与えることにより、震源データ32
の初期設定を行なう(ステップS403)。例えば、ス
テップS402で硬質基盤における値に変換されたデー
タの中で、最大値を有する観測地点を震央と仮定する。
即ち、ここで設定された震源データ32は、仮の震源
(震央)に係るものであるから、便宜上、仮震源データ
32と呼ぶこととする。また、仮震源データ32は、震
源データと同様、震央位置を表す線分(x1,y1)−
(x2,y2)と、震源深さdと、地震規模を表すマグ
ニチュードMの各要素を含むデータである。
Next, in performing the calculation, an hypothetical value is given to the epicenter data 32 so that the epicenter data 32
(Step S403). For example, it is assumed that the observation point having the maximum value is the epicenter in the data converted into the value on the hard base in step S402.
That is, since the epicenter data 32 set here relates to a temporary epicenter (central epicenter), it is referred to as temporary epicenter data 32 for convenience. The temporary epicenter data 32 is, similarly to the epicenter data, a line segment (x1, y1)-
This is data including (x2, y2), an epicenter depth d, and each element of magnitude M representing the magnitude of the earthquake.

【0042】仮震源データ32の初期設定を行なった
ら、震央距離と揺れの強さの関係式データ21に対し
て、当該仮震源データ32の初期値を係数として代入
し、各観測地点の硬質基盤における揺れの強さを理論的
に求める(ステップS404)、即ち、仮震源データ3
2の初期値に従った理論曲線を求める。例えば、震央距
離と基盤最大速度の関係式vb=f(M,d,r)を用
いた場合、震央位置(x1,y1)−(x2,y2)
と、各観測地点の座標(x(n),y(n))から震央
距離r=f(x1,y1,x2,y2,x(n),y
(n))を求めて、同式に代入すれば、各観測地点での
硬質な基盤における最大速度vbが演算できる。尚、こ
こで言う最大速度vbは、前述の通り、理論値であり、
後述する繰り返しの処理を行ない得る関係から、初期理
論値と呼ぶこととしても良い。
After the initial setting of the temporary epicenter data 32, the initial value of the temporary epicenter data 32 is substituted into the relational expression data 21 of the epicenter distance and the shaking intensity as a coefficient, and the rigid base at each observation point is obtained. Is obtained theoretically (step S404), that is, the hypocenter data 3
A theoretical curve according to the initial value of 2 is obtained. For example, when the relational expression v b = f (M, d, r) between the epicenter distance and the base maximum velocity is used, the epicenter position (x1, y1) − (x2, y2)
And the epicenter distance r = f (x1, y1, x2, y2, x (n), y from the coordinates (x (n), y (n)) of each observation point
By calculating (n)) and substituting it into the same equation, the maximum velocity v b on the rigid base at each observation point can be calculated. Note that the maximum speed v b here is a theoretical value as described above,
It may be referred to as an initial theoretical value because of the relationship in which iterative processing described later can be performed.

【0043】各観測地点の実測データと、ステップS4
04で求めた各観測地点における理論値を比較し、誤差
の総和を求める(ステップS405)。例えば、実測デ
ータから変換された基盤最大速度vb(n)と、関係式
データ21から演算された基盤最大速度vbにより、誤
差の総和q=Σ(vb−vb(n))2を求める。
Actual measurement data at each observation point,
The theoretical values at each observation point obtained in step 04 are compared, and the total error is obtained (step S405). For example, a converted from measured data were based maximum speed v b (n), the foundation maximum velocity v b which is calculated from the relationship data 21, the error of the sum q = Σ (v b -v b (n)) 2 Ask for.

【0044】ステップS405で求めた誤差の総和が極
小値であるか否かについて、判断する(ステップS40
6)。
It is determined whether or not the sum of the errors obtained in step S405 is a minimum value (step S40).
6).

【0045】ステップS406において誤差の総和が極
小値であった場合、対応する仮震源データ32を推定さ
れた震源に係る震源データ32として決定し、演算処理
装置3の内部に記録すると共に、理論曲線と各実測デー
タとの誤差の関係を、図6に示した例の様に直感的に把
握できるグラフとしてグラフ化し、当該グラフを表示装
置4に出力する(ステップS407)。
If the sum of the errors is a minimum value in step S406, the corresponding provisional hypocenter data 32 is determined as the hypocenter data 32 relating to the estimated hypocenter, and is recorded in the arithmetic processing unit 3 and the theoretical curve is calculated. The relationship between the error and the actual measurement data is graphed as a graph that can be intuitively grasped, as in the example shown in FIG. 6, and the graph is output to the display device 4 (step S407).

【0046】一方、ステップS406における判断の結
果、誤差の総和が極小値でなかった場合、仮震源データ
32の各要素の組み合わせを変更(ステップS408)
し、ステップS403からステップS406までを繰り
返す。具体的には、先ず、仮震源データ32の各要素の
変化量を大きくして広範囲で演算を繰り返し、その中で
誤差の総和が一番小さくなる各要素の組み合わせを概算
値として求め、次に、その付近を微少量変化させて詳細
な組み合わせを求める。
On the other hand, if the result of determination in step S 406 is that the sum of the errors is not the minimum value, the combination of the elements of the temporary hypocenter data 32 is changed (step S 408).
Then, steps S403 to S406 are repeated. Specifically, first, the amount of change of each element of the temporary hypocenter data 32 is increased, and the calculation is repeated over a wide range. Among them, a combination of each element with the smallest sum of errors is obtained as an approximate value. , A small amount is changed in the vicinity thereof to obtain a detailed combination.

【0047】次に、被害推定手段33の詳細な動作を、
図5に示されるフローチャートを参照して説明する。
Next, the detailed operation of the damage estimation means 33 will be described.
This will be described with reference to the flowchart shown in FIG.

【0048】震源推定手段31で求めた震源データ32
を、ステップS404で用いた式と同一の式を示す関係
式データ21に係数として代入し、また地域内のメッシ
ュ状の各地点について震央距離を計算して同式の変数と
することで、地域全体の硬質基盤の揺れの強さを単純計
算にて演算する(ステップS501)。例えばメッシュ
状の各地点(x(m),y(m))における基盤最大速
度vbは、震央距離と基盤最大速度の関係式vb=f
(M,d,r)と震央距離r=f(x1,y1,x2,
y2,x(m),y(m))から演算できる。
Epicenter data 32 obtained by epicenter estimation means 31
Is substituted into the relational expression data 21 indicating the same expression as the expression used in step S404 as a coefficient, and the epicenter distance is calculated for each mesh-like point in the region, and is used as a variable of the same expression. The strength of the sway of the entire hard base is calculated by a simple calculation (step S501). For example, the maximum base velocity v b at each mesh-like point (x (m), y (m)) is expressed by the relational expression v b = f between the epicenter distance and the maximum base velocity.
(M, d, r) and epicenter distance r = f (x1, y1, x2,
y2, x (m), y (m)).

【0049】メッシュ状の各地点ごとに、ステップS5
01で求めた基盤における揺れの強さと、ステップS4
02で用いた地域特性データ22とを用いて、地表にお
ける揺れの強さを求める。更に、地表における揺れの強
さと、別の地域特性データ22とを用いて、メッシュ状
の各地点における被害推定データ34を演算する(ステ
ップS502)。例えば、地域特性データ22の一つで
ある座標位置に依存する人口データa(m)と、地表に
おける揺れの強さを示す一つの尺度である地表最大速度
v(m)とを用いて、被害推定データの一要素である人
的被害データA(m)=f(v(m),a(m))を演
算する。ここで、被害推定データ34は、上記人的被害
の他に、一次被害(揺れの強さ)、建物被害、火災被害
などの各要素を含むものであり、一次被害の他の要素
は、夫々、地域特性に応じて一次被害との分布傾向が異
なるものである。
Step S5 is performed for each of the mesh points.
01 and the strength of the shaking on the base determined in step S4
The shaking strength on the ground surface is obtained using the regional characteristic data 22 used in 02. Further, damage estimation data 34 at each point in a mesh is calculated using the intensity of shaking on the ground surface and another regional characteristic data 22 (step S502). For example, using population data a (m) depending on the coordinate position, which is one of the regional characteristic data 22, and the maximum ground surface speed v (m), which is one measure of the intensity of shaking on the ground, The human damage data A (m) = f (v (m), a (m)), which is one element of the estimated data, is calculated. Here, the damage estimation data 34 includes each element such as primary damage (strength of shaking), building damage, and fire damage in addition to the above-mentioned human damage. The distribution tendency with primary damage differs depending on the regional characteristics.

【0050】ステップS502で求めた被害推定データ
34の各要素を演算処理装置3の内部に記録すると同時
に、演算結果を表示装置4に出力する(ステップS50
3)。
Each element of the damage estimation data 34 obtained in step S502 is recorded in the arithmetic processing device 3 and at the same time, the calculation result is output to the display device 4 (step S50).
3).

【0051】以上説明したような処理を実行することに
より、本実施の形態による地震被害推定システムにおい
ては、地震被害に関する推定を高速に行なうことが可能
となる。
By executing the processing as described above, the earthquake damage estimation system according to the present embodiment can estimate earthquake damage at high speed.

【0052】また、本実施の形態による地震被害推定シ
ステムは、過去の地震データをもとにして曲線の勾配に
理論的根拠を持たせた関係式を実測データで最適化して
いることから、実測データに近似していながらも、測定
誤差の影響を受けにくい演算結果を出力することができ
るという特徴を有する。
Further, the earthquake damage estimation system according to the present embodiment optimizes the relational expression in which the slope of the curve has a theoretical basis based on the past seismic data with the actually measured data. Although it is close to data, it has a feature that it can output a calculation result that is hardly affected by a measurement error.

【0053】更に、本実施の形態による地震被害推定シ
ステムによれば、例えば、震央付近の揺れが激しい地点
における計測震度計の故障などにより、値の大きい実測
データが入手できなくても、震源地を推定することがで
きる。同様に、震央が遠方にあり、地域内の実測データ
が震央付近の値より小さい場合も推定できる。本発明に
よるシステムがこのような推定を行なうことができるの
は、上述した高精度の被害推定を行なうことができる理
由に起因する、即ち、実測データが欠測しても、他の実
測データでこれを補うことができるためである。
Further, according to the earthquake damage estimating system according to the present embodiment, even if actual measurement data having a large value cannot be obtained due to, for example, a failure of the measuring seismometer at a location where the shaking near the epicenter is severe, the epicenter of the epicenter can be obtained. Can be estimated. Similarly, it can be estimated when the epicenter is far away and the measured data in the area is smaller than the value near the epicenter. The reason that the system according to the present invention can perform such an estimation is due to the reason that the above-described high-precision damage estimation can be performed, that is, even if actual measurement data is missing, other actual measurement data can be used. This is because this can be compensated.

【0054】また、本発明の地震被害推定システムは、
震源データの要素に線状の震央位置を定義して、震源推
定を行っていることから、従来例では対応できなかった
震源が幅を有するような地震、例えば、地表近くの断層
面で発生した地震にも適応することができる。
Also, the earthquake damage estimation system of the present invention
Since the hypocenter estimation was performed by defining the linear epicenter position as an element of the epicenter data, an earthquake that had a width of the epicenter that could not be handled in the conventional example, such as a fault near the ground surface, occurred It can adapt to earthquakes.

【0055】尚、本発明は、上述した本発明の実施の形
態による地震被害推定システムのみに制限される訳でな
く、以下に示すような種々の応用が可能である。
The present invention is not limited to only the earthquake damage estimation system according to the above-described embodiment of the present invention, but can be applied to various applications as described below.

【0056】例えば、上述した実施の形態に対し、次の
ような応用例が考えられる。
For example, the following application is conceivable with respect to the above-described embodiment.

【0057】即ち、上述した実施の形態においては(図
3参照)、地震が発生すると、震度計11で観測された
揺れの強さの実測データが受信処理装置12から演算処
理装置3に渡される。ここで、全ての実測データが揃わ
ない場合において、演算処理装置3は、図4及び図5に
示したフローチャートによる処理と同一の処理を少数の
実測データにて実行するものとする。その後、受信処理
装置12から残りの実測データを受け取ると、図4及び
図5に示したフローチャートによる処理と同一の処理を
再度実行する。このように入手した実測データを最大限
に利用するための段階的な処理を行なうことにより、実
測データの入手に伴い迅速に演算結果を得ることが出来
ると共に、更なる実測データの入手により精度の向上し
た演算結果を得ることができる。
That is, in the above-described embodiment (see FIG. 3), when an earthquake occurs, actual measurement data of the shaking intensity observed by the seismometer 11 is passed from the reception processing device 12 to the arithmetic processing device 3. . Here, when all the measured data are not collected, the arithmetic processing unit 3 executes the same processing as the processing according to the flowcharts shown in FIGS. 4 and 5 with a small number of measured data. Thereafter, when the remaining actual measurement data is received from the reception processing device 12, the same processing as the processing according to the flowcharts shown in FIGS. 4 and 5 is executed again. By performing the stepwise processing for maximizing the use of the actually measured data obtained in this way, it is possible to obtain a calculation result quickly with the acquisition of the actually measured data, and to obtain the accuracy by further obtaining the actually measured data. An improved calculation result can be obtained.

【0058】また、上述した実施の形態による地震被害
推定システムは、図7に示されるようなシステムとして
変形応用することが可能である。
The earthquake damage estimation system according to the above-described embodiment can be modified and applied as a system as shown in FIG.

【0059】即ち、図7に示される地震被害推定システ
ムは、図3における地震観測装置1に変えて、データ入
力装置5を備えている。データ入力装置5は、キーボー
ド、マウス、タッチパネルその他の入力装置であってよ
い。
That is, the earthquake damage estimation system shown in FIG. 7 includes a data input device 5 instead of the earthquake observation device 1 in FIG. The data input device 5 may be a keyboard, a mouse, a touch panel, or another input device.

【0060】このような構成を備えたシステムにおいて
は、地震が発生して地域内外の揺れの強さの実測データ
が観測されると、その実測データがデータ入力装置5か
ら演算処理装置3に入力され、演算処理装置3では、図
2と図3に示したフローチャートによる処理と同一の処
理が実行される。結果として、上述した実施の形態によ
るシステムと同等の効果を得ることができる。
In a system having such a configuration, when an actual measurement data of shaking intensity inside and outside the area is observed due to the occurrence of an earthquake, the actual measurement data is input from the data input device 5 to the arithmetic processing device 3. Then, in the arithmetic processing unit 3, the same processing as the processing according to the flowcharts shown in FIGS. 2 and 3 is executed. As a result, an effect equivalent to that of the system according to the above-described embodiment can be obtained.

【0061】更に、図8に示されるように、本発明の概
念による地震被害推定処理を演算処理装置に実行させ得
るプログラムを記録した記録媒体へ適用することも可能
である。
Further, as shown in FIG. 8, it is also possible to apply a seismic damage estimation process according to the concept of the present invention to a recording medium that stores a program that allows an arithmetic processing unit to execute the process.

【0062】図8に示される記録媒体6は、情報処理装
置7で読み取り可能なものであれば、磁気ディスク、光
ディスク、半導体メモリその他の記録媒体であってよ
い。この記録媒体には、図4及び図5を用いて説明した
処理を情報処理装置7に実行させ得るプログラムが格納
されている。このプログラムは、記録媒体6から情報処
理装置7に読み込まれ、情報処理装置7における動作を
制御する。情報処理装置7は、地震被害推定プログラム
に従って、上述した実施の形態における演算処理装置3
による処理と同一の処理を実行する。
The recording medium 6 shown in FIG. 8 may be a magnetic disk, an optical disk, a semiconductor memory or any other recording medium as long as it can be read by the information processing device 7. This recording medium stores a program that allows the information processing device 7 to execute the processing described with reference to FIGS. 4 and 5. This program is read from the recording medium 6 into the information processing device 7 and controls the operation in the information processing device 7. The information processing device 7 executes the arithmetic processing device 3 according to the above-described embodiment in accordance with the earthquake damage estimation program.
Execute the same processing as the processing by.

【0063】[0063]

【発明の効果】本発明によれば、地震被害に関する推定
を高速に行なうことが可能となる。詳しくは、本発明に
おいては、地域内の揺れの強さの実測データのみの入力
により、これら複数の実測データから震源データが推定
でき、その後は単純計算にて地域全体の被害分布が演算
できる。結果として、本発明の地震被害推定システム
は、地震が発生して気象庁の震源情報が発令される前
に、震源情報と地域全体の被害分布を高速に算出でき
る。
According to the present invention, it is possible to quickly estimate earthquake damage. More specifically, in the present invention, by inputting only the actual measurement data of the shaking strength in the area, the epicenter data can be estimated from the plurality of actual measurement data, and thereafter, the damage distribution of the entire area can be calculated by simple calculation. As a result, the earthquake damage estimation system of the present invention can quickly calculate the epicenter information and the damage distribution of the entire area before the occurrence of an earthquake and the issuance of epicenter information by the Japan Meteorological Agency.

【0064】また、本発明の地震被害推定システムは、
高精度の被害推定を行なうことができる。即ち、本発明
の地震被害推定システムでは、過去の地震データをもと
にして曲線の勾配に理論的根拠を持たせた関係式を実測
データで最適化していることから、本発明のシステム
は、実測データに近似していながらも、測定誤差の影響
を受けにくい演算結果を出力することができるという特
徴を有する。
Further, the earthquake damage estimation system of the present invention
Highly accurate damage estimation can be performed. That is, in the earthquake damage estimation system of the present invention, since the relational expression giving the theoretical basis to the slope of the curve based on the past earthquake data is optimized by the actually measured data, the system of the present invention is: Although it is close to the actual measurement data, it has a feature that it can output a calculation result which is hardly affected by the measurement error.

【0065】また、本発明によれば、震央付近の揺れが
激しい地点における計測震度計の故障などにより、値の
大きい実測データが入手できなくても、震源地を推定す
ることができる。同様に、震央が遠方にあり、地域内の
実測データが震央付近の値より小さい場合も推定でき
る。本発明によるシステムがこのような推定を行なうこ
とができるのは、上述した高精度の被害推定を行なうこ
とができる理由に起因する、即ち、実測データが欠測し
ても、他の実測データでこれを補うことができるためで
ある。
Further, according to the present invention, the epicenter can be estimated even if actual measurement data with a large value cannot be obtained due to a failure of the measuring seismometer at a location where the shaking near the epicenter is severe. Similarly, it can be estimated when the epicenter is far away and the measured data in the area is smaller than the value near the epicenter. The reason that the system according to the present invention can perform such an estimation is due to the reason that the above-described high-precision damage estimation can be performed, that is, even if actual measurement data is missing, other actual measurement data can be used. This is because this can be compensated.

【0066】更に、本発明の地震被害推定システムは、
従来例では対応できなかった震源が幅を有するような地
震、例えば、地表近くの断層面で発生した地震にも適応
することができる。その理由は、震源データの要素に線
状の震央位置を定義して、震源推定を行っているためで
ある。
Further, the earthquake damage estimation system of the present invention
It can also be applied to earthquakes whose epicenters have a width that could not be handled in the conventional example, for example, earthquakes that occurred on a fault plane near the ground surface. The reason is that the hypocenter is defined by defining the linear epicenter position as an element of the epicenter data.

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

【図1】本発明の実施の形態を説明に用いられる語句の
定義を示す図である。
FIG. 1 is a diagram showing definitions of terms used for describing an embodiment of the present invention.

【図2】図1と同様に語句の定義を示す図である。FIG. 2 is a diagram showing definitions of words and phrases as in FIG. 1;

【図3】本発明の実施の形態による地震被害推定システ
ムの構成を示すブロック図である。
FIG. 3 is a block diagram illustrating a configuration of an earthquake damage estimation system according to an embodiment of the present invention.

【図4】図3に示される震源推定手段31における震源
推定処理の手順を示すフローチャートである。
FIG. 4 is a flowchart showing a procedure of a hypocenter estimation process in a hypocenter estimation unit 31 shown in FIG. 3;

【図5】図3に示される被害推定手段33における被害
推定処理の手順を示すフローチャートである。
FIG. 5 is a flowchart showing a procedure of damage estimation processing in damage estimation means 33 shown in FIG. 3;

【図6】震源データを直感的に把握しやすくするための
グラフである。
FIG. 6 is a graph for making it easier to intuitively grasp epicenter data.

【図7】図3に示される地震被害推定システムの変形例
を示す図である。
FIG. 7 is a diagram showing a modification of the earthquake damage estimation system shown in FIG.

【図8】図3に示される地震被害推定システムの応用例
を示す図である。
FIG. 8 is a diagram showing an application example of the earthquake damage estimation system shown in FIG.

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

1 地震観測装置 11 震度計 12 受信処理装置 2 記憶装置 21 関係式データ 22 地域特性データ 3 演算処理装置 31 震源推定手段 32 震源データ 33 被害推定手段 34 被害推定データ 4 表示装置 5 データ入力装置 6 記憶媒体 7 情報処理装置 DESCRIPTION OF SYMBOLS 1 Earthquake observation device 11 Seismic intensity meter 12 Reception processing device 2 Storage device 21 Relational expression data 22 Regional characteristic data 3 Arithmetic processing device 31 Epicenter estimation means 32 Epicenter data 33 Damage estimation means 34 Damage estimation data 4 Display device 5 Data input device 6 Storage Medium 7 Information processing device

───────────────────────────────────────────────────── フロントページの続き (72)発明者 原 正仁 東京都港区芝五丁目7番1号 日本電気 株式会社内 (72)発明者 清水 洋 東京都港区芝五丁目7番1号 日本電気 株式会社内 (72)発明者 葛野 浩史 東京都港区芝五丁目7番1号 日本電気 株式会社内 (72)発明者 芝野 洋 東京都港区芝五丁目7番1号 日本電気 株式会社内 (72)発明者 坂本 朗一 東京都中央区日本橋久松町4番7号 日 本橋エビスビル 株式会社防災・情報研 究所内 (72)発明者 吉井 博明 神奈川県茅ヶ崎市行谷1100 学校法人 文教大学学園内 (56)参考文献 特開 平6−324160(JP,A) 特開 平8−329043(JP,A) 特開 平8−75865(JP,A) (58)調査した分野(Int.Cl.7,DB名) G01V 1/28 G01V 1/00 ────────────────────────────────────────────────── ─── Continuing on the front page (72) Inventor Masahito Hara 5-7-1 Shiba, Minato-ku, Tokyo NEC Corporation (72) Inventor Hiroshi Shimizu 5-7-1 Shiba, Minato-ku, Tokyo NEC (72) Inventor Hiroshi Kuzuno 5-7-1 Shiba, Minato-ku, Tokyo NEC Corporation (72) Inventor Hiroshi Shibano 5-7-1 Shiba, Minato-ku, Tokyo NEC Corporation ( 72) Inventor Ryuichi Sakamoto 4-7 Nihonbashi Hisamatsu-cho, Chuo-ku, Tokyo Nihonbashi Ebisu Building Disaster Prevention and Information Research Institute Co., Ltd. 56) References JP-A-6-324160 (JP, A) JP-A-8-329043 (JP, A) JP-A-8-75865 (JP, A) (58) Fields investigated (Int. Cl. 7 , DB name) G01V 1/2 8 G01V 1/00

Claims (10)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】 地震を観測する複数の観測地点で、夫
々、得られた複数の実測データを入力とし、推定対象地
域における地震被害を推定するための地震被害推定シス
テムにおいて、 前記複数の実測データを入力するための入力手段と、 震央位置からの地表上における距離である震央距離と地
震の揺れの強さとの関係を示す関係式データを含み、地
震被害推定の演算に用いられる基礎情報を予め記憶する
記憶手段と、 前記複数の実測データに基づいて震源データを算出する
と共に、当該算出結果を前記関係式データに与えてメッ
シュ状の被害分布を演算するための演算処理装置であっ
て、前記震源データを算出する際に、複数の仮の震源を
想定し、前記複数の仮の震源の夫々に係る複数の仮震源
データに従って、当該仮震源データを前記関係式データ
に代入して求められた理論曲線と揺れの強さに係る実測
データとの誤差の総和を算出する処理を、前記複数の仮
震源データ複数の前記震源データの夫々について実行
し、その算出結果から前記誤差の総和を極小とする一の
前記仮震源データを震源に係る震源データとして推定決
定する演算処理装置と、 当該演算処理装置からの演算結果を受けて、地震被害推
定結果として表示するための表示手段とを備えることを
特徴とする地震被害推定システム。
An earthquake damage estimating system for estimating earthquake damage in an estimation target area by inputting a plurality of obtained actual measurement data at a plurality of observation points for observing an earthquake, respectively, Input means for inputting the information, and relational expression data indicating the relationship between the epicenter distance, which is the distance on the ground from the epicenter position, and the intensity of the shaking of the earthquake. Storage means for storing, and an arithmetic processing device for calculating hypocenter data based on the plurality of actually measured data, and applying a calculation result to the relational expression data to calculate a mesh-like damage distribution, When calculating the epicenter data, assuming a plurality of temporary epicenters, and according to the plurality of temporary epicenter data pertaining to each of the plurality of temporary epicenters, convert the temporary epicenter data into the relational expression. The process of calculating the sum of the error between the theoretical curve obtained by substituting the data and the measured data related to the magnitude of the sway is executed for each of the plurality of temporary source data and the plurality of source data, and the calculation result An arithmetic processing unit for estimating and determining one of the temporary epicenter data as the epicenter data relating to the epicenter, which minimizes the sum of the errors, and receiving the arithmetic result from the arithmetic processing unit and displaying the result as an earthquake damage estimation result. An earthquake damage estimation system, comprising:
【請求項2】 請求項1に記載の地震被害推定システム
において、 前記演算処理装置は、前記震央位置を地表上における線
状の位置座標として定義すると共に、前記震央距離を当
該線状の震央位置からの前記地表上における最短距離と
して算出し、当該算出結果を前記関係式データに代入す
ることにより、前記被害分布を演算するものであること
を特徴とする地震被害推定システム。
2. The earthquake damage estimation system according to claim 1, wherein the arithmetic processing unit defines the epicenter position as a linear position coordinate on the surface of the ground, and defines the epicenter distance as the linear epicenter position. And calculating the damage distribution by calculating the distance as the shortest distance from the ground surface on the ground surface and substituting the calculation result into the relational expression data.
【請求項3】 請求項1又は請求項2のいずれかに記載
の地震被害推定システムにおいて、 前記表示手段は、前記地震被害推定結果を、視覚的に把
握しやすくするために、グラフ表示するものである こ
とを特徴とする地震被害推定システム。
3. The earthquake damage estimation system according to claim 1, wherein the display means displays the earthquake damage estimation result in a graph so that the result can be easily grasped visually. An earthquake damage estimation system, characterized in that:
【請求項4】 地震被害の推定対象地域内における複数
の観測地点において実際に観測された結果としての実測
データに基づいて、前記推定対象地域における被害分布
を推定するための推定方法であって、 複数の仮の震源を想定し、前記複数の実測データに基づ
いて、前記複数の仮の震源に係る複数の仮震源データで
あって、夫々に震央位置及び地震規模を示すマグニチュ
ードを内包する仮震源データを求める第1のステップ
と、 震央距離と揺れの強さを理論的に示す関係式に対して、
前記仮震源データの夫々を代入することにより、震源に
係る震源データを決定し、当該震源データに基づいて、
メッシュ上の被害分布について推定する第2のステップ
とを有することを特徴とする地震被害の推定方法。
4. An estimation method for estimating a damage distribution in the estimation target area based on actual measurement data as a result of actually being observed at a plurality of observation points in the estimation target area of the earthquake damage, Assuming a plurality of temporary hypocenters, based on the plurality of actual measurement data, a plurality of temporary epicenter data relating to the plurality of temporary epicenters, each including a magnitude indicating an epicenter position and a magnitude of an earthquake. The first step of obtaining the data, and the relational expression that theoretically shows the epicenter distance and the strength of the sway,
By substituting each of the temporary epicenter data, the epicenter data relating to the epicenter is determined, and based on the epicenter data,
And a second step of estimating a damage distribution on the mesh.
【請求項5】 請求項4に記載の地震被害の推定方法で
あって、 前記第2のステップは、前記関係式として、過去の地震
データに基づいて、曲線の勾配に理論的な根拠を持たせ
てなる一般的な相関式を利用するものであることを特徴
とする地震被害の推定方法。
5. The method for estimating earthquake damage according to claim 4, wherein the second step has a theoretical basis for the slope of the curve based on past earthquake data as the relational expression. A method for estimating earthquake damage, characterized by utilizing a general correlation formula.
【請求項6】 請求項5に記載の地震被害の推定方法で
あって、 前記第2のステップは、前記相関式に対して前記仮震源
データを与えることにより得られる理論曲線と、前記各
観測地点における実測データとの誤差の総和を前記仮震
源データの夫々について算出し、当該算出結果に従い、
前記誤差の総和を極小とする一の仮震源データを震源デ
ータとして決定し、当該震源データを震源に係る推定結
果とするものであることを特徴とする地震被害の推定方
法。
6. The method for estimating earthquake damage according to claim 5, wherein in the second step, a theoretical curve obtained by giving the temporary hypocenter data to the correlation equation, and each of the observations Calculate the sum of the errors with the actual measurement data at each point for each of the temporary hypocenter data, and according to the calculation results,
A method for estimating earthquake damage, wherein one temporary epicenter data that minimizes the sum of the errors is determined as epicenter data, and the epicenter data is used as an estimation result of the epicenter.
【請求項7】 請求項4乃至請求項6のいずれかに記載
の地震被害の推定方法において、 前記震源データ及び仮震源データは、その一要素とし
て、震源の地表上における位置を示す震央位置を含み、 当該震央位置は、線状で示されるものであることを特徴
とする地震被害の推定方法。
7. The method for estimating earthquake damage according to claim 4, wherein the epicenter data and the temporary epicenter data include, as one element thereof, an epicenter position indicating a position of the epicenter on the ground surface. A method for estimating earthquake damage, wherein the epicenter position is indicated by a line.
【請求項8】 請求項4乃至請求項7のいずれかに記載
の地震被害の推定方法において、 前記第2のステップにおける推定結果をグラフとして示
す第3のステップを更に備えることを特徴とする地震被
害の推定方法。
8. The method for estimating earthquake damage according to claim 4, further comprising: a third step of displaying a result of the estimation in the second step as a graph. Damage estimation method.
【請求項9】 データを格納する記憶手段を備えると共
にプログラムに従う処理を実行する情報処理装置に対し
て、地震被害の推定対象地域内における複数の観測地点
において実際に観測された結果としての実測データに基
づいて、前記推定対象地域における被害分布を推定させ
るために、 複数の仮の震源を想定し、前記複数の実測データの夫々
に基づいて、前記複数の震源に係る複数の仮震源データ
であって、夫々に震央位置及び地震規模を示すマグニチ
ュードを内包する仮震源データを求める第1の処理と、 震央距離と揺れの強さを理論的に示す関係式として、過
去の地震データに基づいて、曲線の勾配に理論的な根拠
を持たせてなる一般的な相関式を用い、当該相関式に対
して、前記第1の処理において求めた前記仮震源データ
の夫々を代入することにより、夫々の仮震源データに従
った理論曲線を算出させる第2の処理と、 前記第2の処理の結果である各理論曲線と前記複数の実
測データとの誤差の総和を算出する第3の処理と、 前記複数の仮震源データの内、前記第3の処理により算
出された前記誤差の総和を極小とする理論曲線に応じた
前記仮震源データを前記震源に係る震源データであると
して推定し、メッシュ上の被害分布について推定する第
4の処理とを前記情報処理装置に実行させる命令を含む
プログラムを格納した記録媒体であって、前記情報処理
装置の読取可能な記録媒体。
9. An information processing apparatus comprising a storage unit for storing data and executing a process according to a program is provided with actual measurement data as a result of actual observation at a plurality of observation points in a target area for estimation of earthquake damage. In order to estimate the damage distribution in the estimation target area based on the above, a plurality of temporary hypocenters are assumed, and based on each of the plurality of actually measured data, a plurality of temporary hypocenter data relating to the plurality of hypocenters are obtained. Based on past earthquake data, a first process for obtaining temporary hypocenter data containing magnitudes indicating the epicenter position and the magnitude of the earthquake, respectively, and a relational expression theoretically indicating the epicenter distance and the intensity of shaking, Using a general correlation equation having a theoretical basis for the slope of the curve, each of the temporary hypocenter data obtained in the first processing is used for the correlation equation. By substituting, a second process for calculating a theoretical curve according to the respective temporary hypocenter data, and calculating a sum of errors between each theoretical curve as a result of the second process and the plurality of measured data. A third process, wherein, among the plurality of temporary source data, the temporary source data corresponding to a theoretical curve that minimizes the sum of the errors calculated by the third process is source data relating to the source; And a fourth process for estimating a damage distribution on a mesh, the program including a command for causing the information processing apparatus to execute the processing, the storage medium being readable by the information processing apparatus.
【請求項10】 請求項9に記載の記録媒体であって前
記震源データは、その一要素として、震源の地表上にお
ける位置を示す震央位置を含み、且つ、当該震央位置
は、線状で示されるものであることを特徴とする記録媒
体。
10. The recording medium according to claim 9, wherein the epicenter data includes, as one element thereof, an epicenter position indicating a position of the epicenter on the ground surface, and the epicenter position is indicated by a line. A recording medium characterized by being recorded.
JP10019774A 1998-01-30 1998-01-30 Earthquake damage estimation system Expired - Lifetime JP3017154B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP10019774A JP3017154B2 (en) 1998-01-30 1998-01-30 Earthquake damage estimation system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP10019774A JP3017154B2 (en) 1998-01-30 1998-01-30 Earthquake damage estimation system

Publications (2)

Publication Number Publication Date
JPH11211839A JPH11211839A (en) 1999-08-06
JP3017154B2 true JP3017154B2 (en) 2000-03-06

Family

ID=12008692

Family Applications (1)

Application Number Title Priority Date Filing Date
JP10019774A Expired - Lifetime JP3017154B2 (en) 1998-01-30 1998-01-30 Earthquake damage estimation system

Country Status (1)

Country Link
JP (1) JP3017154B2 (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002116263A (en) * 2000-10-05 2002-04-19 Shimizu Corp System and method for processing/transmitting earthquake information
JP2002148348A (en) * 2000-11-06 2002-05-22 Taisei Corp Earthquake risk evaluation system
JP2004233053A (en) * 2000-12-12 2004-08-19 Hiroyuki Inubushi Apparatus and method for solving problem occurring immediately before or after earthquake
JP5281284B2 (en) * 2004-03-11 2013-09-04 リスク マネジメント ソリューションズ, インコーポレイテッド System and method for manufacturing a flexible geographic grid
US7707050B2 (en) 2004-03-11 2010-04-27 Risk Management Solutions, Inc. Systems and methods for determining concentrations of exposure
US8229766B2 (en) 2004-07-30 2012-07-24 Risk Management Solutions, Inc. System and method for producing a flexible geographical grid
EP1894159A1 (en) 2005-06-10 2008-03-05 Swiss Reinsurance Company Computer system and method for determining an earthquake damage index
JP4805993B2 (en) * 2008-09-30 2011-11-02 株式会社竹中工務店 Seismic amplitude level prediction device, control device, and program
CN103376465B (en) * 2012-04-25 2016-01-27 林沛旸 The earthquake instant analysis system of buildings floor and method thereof

Also Published As

Publication number Publication date
JPH11211839A (en) 1999-08-06

Similar Documents

Publication Publication Date Title
JP6027797B2 (en) Building earthquake resistance evaluation system and building earthquake resistance evaluation method
Lund et al. Calculating horizontal stress orientations with full or partial knowledge of the tectonic stress tensor
JP4253435B2 (en) Method and apparatus for estimating seismic intensity
Gallovič et al. On strong ground motion synthesis with k− 2 slip distributions
JP6552928B2 (en) Weather forecasting device, weather forecasting method and program
CN105445756A (en) Safety diagnosis system for structure
JP3017154B2 (en) Earthquake damage estimation system
D’Amico et al. Fling effects from near‐source strong‐motion records: Insights from the 2016 Mw 6.5 Norcia, central Italy, earthquake
JP5126143B2 (en) Earthquake motion prediction system
Price et al. Interactions between the Landers and Hector Mine, California, earthquakes from space geodesy, boundary element modeling, and time-dependent friction
CA2744089A1 (en) Method of processing measured data
WO2023121612A1 (en) A digital twin system for predicting response of structures to natural disasters
CN114207479A (en) Earthquake observation device, earthquake observation method, and recording medium for recording earthquake observation program
JP6441869B2 (en) Building earthquake resistance evaluation system and building earthquake resistance evaluation method
Neighbors et al. Solving for source parameters using nested array data: A case study from the Canterbury, New Zealand earthquake sequence
JP6103516B2 (en) Tsunami warning system, tsunami warning method, and program for tsunami warning system
Ghayamghamian et al. The effect of torsional ground motion on structural response: code recommendation for accidental eccentricity
Touhami et al. A holistic approach of numerical analysis of the geology effects on ground motion prediction: Argostoli site test
JP7402336B2 (en) Estimation system, estimation method, and estimation program
JP4846779B2 (en) Disaster emergency ground fluctuation analysis system and disaster emergency ground fluctuation analysis method
Kohler et al. Detection of building damage using Helmholtz tomography
JP6642232B2 (en) Earthquake damage estimation system, structure with earthquake damage estimation system, and earthquake damage estimation program
JP2004362311A (en) Earthquake damage evaluation system, earthquake damage evaluation method, earthquake damage evaluation program and storage medium recorded with earthquake damage evaluation program
JP7258682B2 (en) Plant evaluation system and plant evaluation method
Geist Rapid tsunami models and earthquake source parameters: far-field and local applications

Legal Events

Date Code Title Description
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 19991201