JP7444811B2 - Building damage estimation method, building damage estimation system, building damage learning method, building damage learning system and program - Google Patents

Building damage estimation method, building damage estimation system, building damage learning method, building damage learning system and program Download PDF

Info

Publication number
JP7444811B2
JP7444811B2 JP2021050741A JP2021050741A JP7444811B2 JP 7444811 B2 JP7444811 B2 JP 7444811B2 JP 2021050741 A JP2021050741 A JP 2021050741A JP 2021050741 A JP2021050741 A JP 2021050741A JP 7444811 B2 JP7444811 B2 JP 7444811B2
Authority
JP
Japan
Prior art keywords
sec
building
information
degree
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.)
Active
Application number
JP2021050741A
Other languages
Japanese (ja)
Other versions
JP2022148889A (en
Inventor
貴浩 佐田
Original Assignee
パナソニックホームズ株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by パナソニックホームズ株式会社 filed Critical パナソニックホームズ株式会社
Priority to JP2021050741A priority Critical patent/JP7444811B2/en
Publication of JP2022148889A publication Critical patent/JP2022148889A/en
Application granted granted Critical
Publication of JP7444811B2 publication Critical patent/JP7444811B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

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

Description

本開示は、建物の被災度推定方法、建物の被災度推定システム、建物の被災度学習方法、建物の被災度学習システム及びプログラムに関する。 The present disclosure relates to a building disaster degree estimation method, a building disaster degree estimation system, a building disaster degree learning method, a building disaster degree learning system, and a program.

特許文献1は、建物被災推定システムおよび方法を開示する。この建物被災推定システムは、建物被災推定装置と、建物の基礎地盤面に設置された地震センサと、建物内の代表ポイントに設置された地震センサから構成される。建物被災推定装置は、建物の既知の建物情報に基づいて、建物の地震による動きを数式化した仮想建物モデルを導出する仮想建物モデル導出部と、地震発生中に地震センサで計測された地動加速度を仮想建物モデルに入力して、仮想建物モデルの地震応答解析を行う地震応答解析部と、地震応答解析の結果を用いて、地震情報と建物の地震による動きを示す動き情報を算出する地震・動き情報算出部と、地震終了後に地震情報と動き情報を用いて、建物の被災度を推定する被災度推定部とを備える。 Patent Document 1 discloses a building damage estimation system and method. This building damage estimation system is composed of a building damage estimation device, an earthquake sensor installed on the foundation ground surface of the building, and an earthquake sensor installed at representative points within the building. The building damage estimation device includes a virtual building model derivation unit that derives a virtual building model that mathematically represents the movement of a building due to an earthquake, based on known building information about the building, and a virtual building model derivation unit that derives a virtual building model that formulates the movement of the building due to an earthquake, and a ground motion acceleration measured by an earthquake sensor during an earthquake. an earthquake response analysis section that inputs information into a virtual building model and performs earthquake response analysis of the virtual building model; It includes a motion information calculation section and a disaster degree estimation section that estimates the degree of damage to a building using earthquake information and movement information after an earthquake.

特許文献1に記載の建物被災推定システムおよび方法にあっては、建物の全ての層に加速度センサを設置する必要がなく、建物の基礎地盤面と建物内の代表ポイントにのみ地震センサを設置すればよい。このため、地震発生後の建物の被災度を従来よりも安価に推定することができるものであった。 In the building damage estimation system and method described in Patent Document 1, there is no need to install acceleration sensors on all layers of the building, and earthquake sensors are installed only on the foundation ground surface of the building and at representative points within the building. Bye. Therefore, the degree of damage to a building after an earthquake can be estimated at a lower cost than before.

特開2016-197013号公報JP 2016-197013 Publication

しかしながら、特許文献1に記載の建物被災推定システムおよび方法にあっては、少なくとも建物内の代表ポイントに地震センサを設置する必要があるため、建物の全ての層に加速度センサを設置する場合と比較すれば安価であるものの、依然としてコストの低減をしにくいという問題があった。 However, in the building damage estimation system and method described in Patent Document 1, it is necessary to install earthquake sensors at least at representative points within the building, so compared to the case where acceleration sensors are installed on all layers of the building. Although it would be cheaper to do so, there was still the problem that it was difficult to reduce costs.

また、地震応答解析部が行う地震応答解析は、具体的には建物剛性の非線形性を考慮した時刻歴応答解析であるため、解析に時間を要するという問題があった。 Furthermore, the seismic response analysis performed by the seismic response analysis department is specifically a time history response analysis that takes into account the nonlinearity of building rigidity, so there is a problem in that the analysis takes time.

本開示は、上記従来の問題点に鑑みて発明したものであって、コストの低減をしやすく、被災度の推定に要する時間を短縮しやすい建物の被災度推定方法、建物の被災度推定システム、建物の被災度学習方法、建物の被災度学習システム及びプログラムを提供することを課題とする。 The present disclosure was invented in view of the above-mentioned conventional problems, and includes a method and system for estimating the degree of damage to a building, which facilitates cost reduction and the time required for estimating the degree of disaster. , our objective is to provide a method for learning the degree of damage to buildings, a system and program for learning the degree of damage to buildings.

本開示の一態様の建物の被災度推定方法は、取得ステップと、推定ステップと、を備える。前記取得ステップは、地震動の最大加速度および前記地震動の時系列波形の応答スペクトル解析におけるスペクトル値の情報を含む入力情報を取得するステップである。前記推定ステップは、予め構築された学習済みモデルを用いて、前記入力情報を基に前記出力情報を推定するステップである。前記学習済みモデルは、前記入力情報を入力データとするものであって、建物の前記地震動による被災度に関する出力情報を教師データとして、機械学習により構築される。 A disaster degree estimation method for a building according to one aspect of the present disclosure includes an acquisition step and an estimation step. The acquisition step is a step of acquiring input information including maximum acceleration of the earthquake motion and information on spectrum values in response spectrum analysis of the time-series waveform of the earthquake motion. The estimation step is a step of estimating the output information based on the input information using a learned model built in advance. The learned model uses the input information as input data, and is constructed by machine learning using output information regarding the damage level of the building due to the earthquake motion as training data.

本開示の一態様の建物の被災度推定システムは、取得部と、推定部と、を備える。前記取得部は、地震動の最大加速度および前記地震動の時系列波形の応答スペクトル解析におけるスペクトル値の情報を含む入力情報を取得する。前記推定部は、予め構築された学習済みモデルを用いて、前記入力情報を基に前記出力情報を推定する。前記学習済みモデルは、前記入力情報を入力データとするものであって、建物の前記地震動による被災度に関する出力情報を教師データとして、機械学習により構築される。 A disaster degree estimation system for a building according to one aspect of the present disclosure includes an acquisition unit and an estimation unit. The acquisition unit acquires input information including information on a maximum acceleration of an earthquake motion and a spectrum value in a response spectrum analysis of a time-series waveform of the earthquake motion. The estimation unit estimates the output information based on the input information using a learned model built in advance. The learned model uses the input information as input data, and is constructed by machine learning using output information regarding the damage level of the building due to the earthquake motion as training data.

本開示の一態様の建物の被災度学習方法は、取得ステップと、学習ステップと、を備える。前記取得ステップは、地震動の最大加速度および前記地震動の時系列波形の応答スペクトル解析におけるスペクトル値の情報を含む入力情報を取得するステップである。前記学習ステップは、前記入力情報を入力データとする学習済みモデルであって、建物の前記地震動による被災度に関する出力情報を教師データとして、前記学習済みモデルを機械学習により構築するステップである。 A building disaster degree learning method according to one aspect of the present disclosure includes an acquisition step and a learning step. The acquisition step is a step of acquiring input information including maximum acceleration of the earthquake motion and information on spectrum values in response spectrum analysis of the time-series waveform of the earthquake motion. The learning step is a step of constructing a trained model using the input information as input data, and using machine learning to construct the trained model using output information regarding the damage level of the building due to the earthquake motion as teacher data.

本開示の一態様の建物の被災度学習システムは、取得部と、学習部と、を備える。前記取得部は、地震動の最大加速度および前記地震動の時系列波形の応答スペクトル解析におけるスペクトル値の情報を含む入力情報を取得する。前記学習部は、前記入力情報を入力データとする学習済みモデルであって、建物の前記地震動による被災度に関する出力情報を教師データとして、前記学習済みモデルを機械学習により構築する。 A building disaster degree learning system according to one aspect of the present disclosure includes an acquisition unit and a learning unit. The acquisition unit acquires input information including information on a maximum acceleration of an earthquake motion and a spectrum value in a response spectrum analysis of a time-series waveform of the earthquake motion. The learning unit is a learned model that uses the input information as input data, and constructs the learned model by machine learning using output information regarding the damage level of the building due to the earthquake motion as teacher data.

本開示の一態様のプログラムは、1以上のプロセッサに、前記建物の被災度推定方法又は前記建物の被災度学習方法を実行させる。 A program according to an aspect of the present disclosure causes one or more processors to execute the method for estimating the degree of disaster for a building or the method for learning the degree of disaster for a building.

本開示の一態様の建物の被災度推定方法、建物の被災度推定システム、建物の被災度学習方法、建物の被災度学習システム及びプログラムによれば、コストの低減をしやすく、被災度の推定に要する時間を短縮しやすい。 According to the building damage estimation method, building damage estimation system, building damage learning method, building damage learning system, and program of one aspect of the present disclosure, costs can be easily reduced and disaster damage estimation It is easy to shorten the time required.

図1は、本開示の一実施形態に係る建物の被災度推定システム及び被災度学習システムのブロック図である。FIG. 1 is a block diagram of a building damage estimation system and a disaster damage learning system according to an embodiment of the present disclosure. 図2は、同上の実施形態の学習済みモデルにおけるニューラルネットワークの構成図である。FIG. 2 is a configuration diagram of a neural network in the trained model of the embodiment described above. 図3は、同上の実施形態における建物の解析モデルの構成図である。FIG. 3 is a configuration diagram of an analytical model of a building in the embodiment same as above. 図4は、同上の実施形態における建物の被災度学習方法のフロー図である。FIG. 4 is a flowchart of a method for learning the degree of damage to a building in the embodiment described above. 図5は、同上の実施形態のA工法の解析モデルのニューラルネットワークにおける重みW1の絶対値の平均値-ノード関係図である。FIG. 5 is a diagram showing the relationship between the average value of the absolute value of the weight W1 and the node in the neural network of the analytical model of method A according to the above embodiment. 図6は、同上の実施形態のB工法の解析モデルのニューラルネットワークにおける重みW1の絶対値の平均値-ノード関係図である。FIG. 6 is a diagram showing the relationship between the average absolute value of the weight W1 and the node in the neural network of the analytical model of method B according to the embodiment. 図7は、同上の実施形態における建物の被災度推定方法のフロー図である。FIG. 7 is a flowchart of a method for estimating the degree of damage to a building in the embodiment described above. 図8Aは、同上の実施形態のA工法の訓練データを用いた被災度の推定値-正解値の関係図である。図8Bは、同上の実施形態のA工法のテストデータを用いた被災度の推定値-正解値の関係図である。FIG. 8A is a relationship diagram between the estimated value of the degree of disaster and the correct value using the training data of method A according to the above embodiment. FIG. 8B is a relationship diagram between the estimated value of the degree of disaster and the correct value using the test data of construction method A according to the above embodiment. 図9Aは、同上の実施形態のB工法の訓練データを用いた被災度の推定値-正解値の関係図である。図9Bは、同上の実施形態のB工法のテストデータを用いた被災度の推定値-正解値の関係図である。FIG. 9A is a relationship diagram between the estimated value of the degree of disaster and the correct value using the training data of method B of the embodiment. FIG. 9B is a diagram showing the relationship between the estimated value of the degree of disaster and the correct value using the test data of method B of the embodiment. 図10Aは、同上の実施形態のB工法のテストデータを用いた被災度の推定値-正解値の関係図である。図10Bは、例1のB工法のテストデータを用いた被災度の推定値-正解値の関係図である。FIG. 10A is a relationship diagram between the estimated value of the degree of disaster and the correct value using the test data of method B of the embodiment same as above. FIG. 10B is a relationship diagram between the estimated value of the degree of disaster and the correct value using the test data of method B of Example 1. 図11Aは、同上の実施形態のB工法のテストデータを用いた被災度の推定値-正解値の関係図である。図11Bは、例2のB工法のテストデータを用いた被災度の推定値-正解値の関係図である。FIG. 11A is a relationship diagram between the estimated value of the degree of disaster and the correct value using the test data of construction method B of the embodiment same as above. FIG. 11B is a relationship diagram between the estimated value of the degree of disaster and the correct value using the test data of the B construction method of Example 2.

(概要)
本開示は、建物の被災度推定方法、建物の被災度推定システム、建物の被災度学習方法、建物の被災度学習システム及びプログラムについてであって、地震が発生した際に、迅速に建物の被災度を推定することを第一の目的としている。また、本開示は、迅速な建物の被災度の推定を低コストで実現することを第二の目的としている。
(overview)
The present disclosure relates to a building damage estimation method, a building damage estimation system, a building damage learning method, a building damage learning system, and a program. The primary purpose is to estimate the degree of A second objective of the present disclosure is to quickly estimate the degree of damage to a building at low cost.

地震による人的被害及び建物の被害は以前より大きな問題となっており、昨今、これらの被害を低減すべく様々な取り組みがなされている。人的被害としては、地震発生時及び地震発生直後に、人が転倒したり落下物と衝突して怪我をしたり、建物の損傷に起因する怪我をしたりする被害が挙げられる。また、地震発生時に人的被害が発生しなくても、建物の被害が発生すると、損傷した建物がさらに損傷したり崩壊したりすることによる人的被害が発生する可能性がある。いずれにしても、建物の損傷は、建物の被害であると共に人的被害の原因となるものであり、人的被害及び建物の被害を最小限に止めるためには、建物の損傷度合(被災度)を速やかに把握する必要がある。 Human casualties and damage to buildings caused by earthquakes have become bigger problems than before, and various efforts have recently been made to reduce these damages. Examples of human damage include injuries caused by people falling down or colliding with falling objects, and injuries caused by damage to buildings during and immediately after an earthquake. Furthermore, even if no human casualties occur during an earthquake, if buildings are damaged, human casualties may occur due to further damage or collapse of the damaged buildings. In any case, building damage is not only damage to the building, but also causes human damage, and in order to minimize human damage and damage to the building, it is necessary to determine the degree of damage to the building (disaster damage level). ) needs to be quickly understood.

このため、住宅メーカをはじめとする一部の企業においては、建物に地震センサを設置して、地震動による建物の振動を直接的に計測することにより、建物の被災度を把握している。しかしながら、このような方法では、建物の被災度を精度良く把握しやすいものの、多数の住戸に地震センサを設置すると共に地震センサからの計測データを送信する環境を整備する必要があり、コストがかかる。また、公的機関等が公表する震度分布から求まる震度のみでは、推定する被災度の精度がよくない。 For this reason, some companies, including housing manufacturers, install earthquake sensors in buildings to directly measure the vibrations of the building due to earthquake motion, thereby ascertaining the degree of damage to the building. However, although this method makes it easy to accurately assess the degree of damage to a building, it requires installing earthquake sensors in many residential units and creating an environment for transmitting measurement data from the earthquake sensors, which is costly. . Furthermore, the accuracy of estimating the degree of damage is not good if only the seismic intensity is determined from the seismic intensity distribution published by public institutions.

本開示の建物の被災度推定方法、建物の被災度推定システム、建物の被災度学習方法、建物の被災度学習システム及びプログラムにおいては、入力情報として、現実の住戸に設置した多数の地震センサからの地震動の時系列波形の情報を不要とし、代わりに、地震動の最大加速度および地震動の時系列波形の応答スペクトル解析におけるスペクトル値の情報を入力情報としている。これにより、建物に地震センサを全く又はほとんど設置する必要がなくなり、コストの低減をしやすくなると共に、建物剛性の非線形性を考慮した時刻歴応答解析を行わずに建物の被災度を推定するため、建物の被災度の推定に要する時間を短縮しやすくなる。 In the building damage estimation method, building damage estimation system, building damage learning method, building damage learning system and program disclosed herein, the input information is obtained from a large number of earthquake sensors installed in actual dwelling units. The information on the time-series waveform of the seismic motion is not required, and instead, the input information is the maximum acceleration of the seismic motion and information on the spectrum value in the response spectrum analysis of the time-series waveform of the seismic motion. This eliminates the need to install any or almost no seismic sensors in buildings, making it easier to reduce costs, and making it possible to estimate the degree of damage to buildings without performing time history response analysis that takes into account the nonlinearity of building stiffness. , it becomes easier to shorten the time required to estimate the degree of damage to a building.

(詳細)
以下、本開示の一実施形態に係る建物の被災度推定方法、建物の被災度推定システム、建物の被災度学習方法、建物の被災度学習システム及びプログラムについて、図11Bを参照して詳細に説明する。本実施形態の建物の被災度推定方法は、図1に示す建物の被災度推定システム1により実行され得る。被災度推定システム1は、通信ネットワーク8を介して地震動波形サーバ9と接続可能である。
(detail)
Hereinafter, a building damage estimation method, a building damage estimation system, a building damage learning method, a building damage learning system, and a program according to an embodiment of the present disclosure will be described in detail with reference to FIG. 11B. do. The building disaster degree estimation method of this embodiment can be executed by the building disaster degree estimation system 1 shown in FIG. The disaster damage estimation system 1 can be connected to the seismic motion waveform server 9 via the communication network 8 .

(1)被災度推定システム
被災度推定システム1は、上述したように、本実施形態の建物の被災度推定方法を実行するためのシステムである。被災度推定システム1は、通信部11と、記憶部12と、処理部13とを備える。被災度推定システム1は、1以上のサーバにより実現され得る。
(1) Disaster damage estimation system As described above, the disaster damage estimation system 1 is a system for executing the method for estimating the disaster damage of a building according to the present embodiment. The disaster damage estimation system 1 includes a communication section 11, a storage section 12, and a processing section 13. The disaster damage estimation system 1 can be realized by one or more servers.

通信部11は、通信インターフェースである。通信部11は、通信ネットワーク8に接続可能であり、通信ネットワーク8を通じた通信を行う機能を有する。通信部11は、例えば、トランスミッタとレシーバを含む。通信部11は、所定の通信プロトコルに準拠している。所定の通信プロトコルは、周知の様々な有線及び無線通信規格から選択され得る。 The communication unit 11 is a communication interface. The communication unit 11 is connectable to the communication network 8 and has a function of communicating through the communication network 8. The communication unit 11 includes, for example, a transmitter and a receiver. The communication unit 11 complies with a predetermined communication protocol. The predetermined communication protocol may be selected from a variety of well-known wired and wireless communication standards.

記憶部12は、処理部13が利用する情報を記憶するために用いられる。記憶部12は、1以上の記憶装置を含む。記憶装置は、例えば、RAM(Random Access Memory)、EEPROM(Electrically Erasable Programmable Read Only Memory)である。 The storage unit 12 is used to store information used by the processing unit 13. The storage unit 12 includes one or more storage devices. The storage device is, for example, a RAM (Random Access Memory) or an EEPROM (Electrically Erasable Programmable Read Only Memory).

処理部13は、被災度推定システム1の動作を制御する制御回路である。処理部13は、例えば、1以上のプロセッサ(マイクロプロセッサ)と1以上のメモリとを含むコンピュータシステムにより実現され得る。つまり、1以上のプロセッサが1以上のメモリに記憶された1以上の(コンピュータ)プログラム(アプリケーション)を実行することで、処理部13として機能する。プログラムは、ここでは処理部13のメモリに予め記録されているが、インターネット等の電気通信回線を通じて、又はメモリカード等の非一時的な記録媒体に記録されて提供されてもよい。 The processing unit 13 is a control circuit that controls the operation of the disaster damage estimation system 1. The processing unit 13 can be realized, for example, by a computer system including one or more processors (microprocessors) and one or more memories. That is, one or more processors function as the processing unit 13 by executing one or more (computer) programs (applications) stored in one or more memories. Although the program is pre-recorded in the memory of the processing unit 13 here, it may also be provided via a telecommunications line such as the Internet or by being recorded on a non-temporary recording medium such as a memory card.

処理部13は、前処理部131と、取得部132と、推定部133と、提示部134と、を備える。前処理部131と、取得部132と、推定部133と、提示部134とは実体のある構成を示しているわけではなく、処理部13によって実現される機能を示している。 The processing unit 13 includes a preprocessing unit 131, an acquisition unit 132, an estimation unit 133, and a presentation unit 134. The preprocessing unit 131, the acquisition unit 132, the estimation unit 133, and the presentation unit 134 do not represent actual configurations, but represent functions realized by the processing unit 13.

(1-1)前処理部
前処理部131は、被災度推定システム1が実行する被災度推定方法において必要となる入力情報を取得するため、地震動発生時の地表における変位の時系列データ、速度の時系列データ又は加速度の時系列データ(これらを総称して地震動の時系列波形というものとする)を基に、入力情報を求める処理を実行する。
(1-1) Preprocessing unit The preprocessing unit 131 acquires input information necessary for the disaster degree estimation method executed by the disaster degree estimation system 1. The process of obtaining input information is executed based on the time series data of or the time series data of acceleration (these are collectively referred to as the time series waveform of seismic motion).

地震動の時系列波形は、通信ネットワーク8を介して地震動波形サーバ9より取得可能である。地震動波形サーバ9としては、例えば、国立研究開発法人防災科学技術研究所が運営するインターネット上のサーバが挙げられ、このサーバから地震動の時系列波形として、「強震観測網(K-NET, KiK-net)」といったデータが取得可能である。「強震観測網(K-NET, KiK-net)」では、全国の多数の観測点におけるデータが取得可能である。 The time-series waveform of the seismic motion can be obtained from the seismic motion waveform server 9 via the communication network 8 . The seismic motion waveform server 9 may be, for example, a server on the Internet operated by the National Research Institute for Earth Science and Disaster Resilience. It is possible to obtain data such as "net)". The Strong Motion Observation Network (K-NET, KiK-net) allows you to obtain data from many observation points across the country.

前処理部131は、地震動波形サーバ9より取得した地震動の時系列波形より、被災度推定システム1が実行する被災度推定方法において必要となる入力情報を求める。入力情報は、地震動の最大加速度、所定加速度継続時間およびスペクトル値の情報を含む。 The preprocessing unit 131 obtains input information necessary for the disaster severity estimation method executed by the disaster severity estimation system 1 from the time-series waveform of the earthquake motion acquired from the earthquake motion waveform server 9 . The input information includes information on the maximum acceleration of the seismic motion, the predetermined acceleration duration, and the spectrum value.

地震動の最大加速度は、いわゆる地表面最大加速度(Peak Ground Acceleration、略してPGAとする)であり、地震動発生時の地表における加速度の最大値の測定値である。PGAは、地震動の時系列波形より求めることができる。 The maximum acceleration of seismic motion is the so-called peak ground acceleration (PGA), which is a measured value of the maximum acceleration on the earth's surface at the time of seismic motion. PGA can be obtained from the time-series waveform of seismic motion.

所定加速度継続時間は、地震動の所定の継続時間に占める、PGAに1未満の正数値を乗じてなる所定加速度を超える地震動の加速度が継続する時間である。PGAに乗じる正数値は、本実施形態では0.8である。所定加速度継続時間は、地震動の時系列波形より求めることができる。 The predetermined acceleration duration is the time during which the acceleration of the seismic motion exceeds the predetermined acceleration obtained by multiplying PGA by a positive value less than 1, which occupies the predetermined duration of the seismic motion. The positive value by which PGA is multiplied is 0.8 in this embodiment. The predetermined acceleration duration can be determined from the time-series waveform of the seismic motion.

スペクトル値は、地震動の加速度の時系列波形の応答スペクトル解析(本実施形態では以下に示すように加速度応答スペクトル解析)により得られる。具体的には、減衰定数を一定とした(本実施形態では5%とした)線形の1質点系せん断モデルに対象の地震動波形を入力して時刻歴応答解析を行うことで最大応答加速度が得られ、1質点系せん断モデルの固有周期を変化させて本解析を繰り返すことで、固有周期と最大応答加速度の関係を示す加速度応答スペクトルが表される。 The spectral value is obtained by a response spectrum analysis (in this embodiment, an acceleration response spectrum analysis as shown below) of a time-series waveform of acceleration of the seismic motion. Specifically, the maximum response acceleration is obtained by inputting the target seismic motion waveform into a linear one-mass system shear model with a constant damping constant (5% in this embodiment) and performing time history response analysis. By repeating this analysis by changing the natural period of the one-mass shear model, an acceleration response spectrum showing the relationship between the natural period and the maximum response acceleration is expressed.

(1-2)取得部
取得部132は、被災度推定システム1が実行する被災度推定方法において必要となる入力情報を取得する。入力情報は、上述したPGA、所定加速度継続時間および地震動の時系列波形の加速度応答スペクトル解析におけるスペクトル値である。
(1-2) Acquisition Unit The acquisition unit 132 acquires input information required in the disaster severity estimation method executed by the disaster severity estimation system 1. The input information is the above-mentioned PGA, the predetermined acceleration duration time, and the spectrum value in the acceleration response spectrum analysis of the time-series waveform of the seismic motion.

(1-3)推定部
推定部133は、予め構築された学習済みモデルを用いて、取得部132で取得した入力情報を基に、建物の地震動による被災度に関する出力情報を推定する。本実施形態においては、推定部133は出力情報として、建物の最大変形量の情報を出力する。
(1-3) Estimating Unit The estimating unit 133 estimates output information regarding the degree of damage caused by earthquake motion of the building based on the input information acquired by the acquiring unit 132 using a trained model built in advance. In this embodiment, the estimation unit 133 outputs information on the maximum amount of deformation of the building as output information.

学習済みモデルは、入力情報を入力データとするものであって、出力情報を教師データとして、機械学習により構築される。本実施形態においては、学習済みモデルは、ニューラルネットワークを用いたものである。学習済みモデルは、建物の被災度学習システム2により構築される。被災度学習システム2については後述する。 A trained model is constructed by machine learning using input information as input data and output information as training data. In this embodiment, the learned model uses a neural network. The trained model is constructed by the building damage level learning system 2. The disaster severity learning system 2 will be described later.

(1-4)提示部
提示部134は、推定部133により推定された出力情報を提示する。本実施形態では、推定部133により推定された出力情報を、通信ネットワーク8を通じて端末装置7に送ることで、端末装置7の入出力部71にて出力情報を提示可能である。
(1-4) Presentation Unit The presentation unit 134 presents the output information estimated by the estimation unit 133. In this embodiment, by sending the output information estimated by the estimation unit 133 to the terminal device 7 through the communication network 8, the output information can be presented at the input/output unit 71 of the terminal device 7.

(2-1)通信ネットワーク
通信ネットワーク8は、インターネットを含み得る。通信ネットワーク8は、単一の通信プロトコルに準拠したネットワークだけではなく、異なる通信プロトコルに準拠した複数のネットワークで構成され得る。通信プロトコルは、周知の様々な有線及び無線通信規格から選択され得る。通信ネットワーク8は、リピータハブ、スイッチングハブ、ブリッジ、ゲートウェイ、ルータ等のデータ通信機器を含み得る。
(2-1) Communication Network The communication network 8 may include the Internet. The communication network 8 may include not only a network based on a single communication protocol but also a plurality of networks based on different communication protocols. The communication protocol may be selected from a variety of well-known wired and wireless communication standards. Communication network 8 may include data communication equipment such as repeater hubs, switching hubs, bridges, gateways, routers, and the like.

(2-2)端末装置
端末装置7は、建物の被災度推定システム1からの情報の表示や、例えば被災度推定システム1への安否情報の送信に利用可能である。端末装置7は、入出力部71と、通信部72と、処理部73とを備える。端末装置7は、デスクトップコンピュータ、ラップトップコンピュータ、又は、携帯端末(スマートフォン、タブレット端末、ウェアラブル端末等)により実現され得る。
(2-2) Terminal Device The terminal device 7 can be used to display information from the building disaster degree estimation system 1 and to transmit safety information to the disaster degree estimation system 1, for example. The terminal device 7 includes an input/output section 71, a communication section 72, and a processing section 73. The terminal device 7 can be realized by a desktop computer, a laptop computer, or a mobile terminal (smartphone, tablet terminal, wearable terminal, etc.).

入出力部71は、端末装置7を操作するための入力装置を備える。入力装置は、例えば、キーボード、マウス、タッチパッド等を含み得る。また、入出力部71は、情報を表示するための画像表示装置を備える。画像表示装置は、液晶ディスプレイや有機ELディスプレイ等の薄型のディスプレイ装置を含み得る。 The input/output unit 71 includes an input device for operating the terminal device 7. Input devices may include, for example, a keyboard, mouse, touch pad, and the like. The input/output unit 71 also includes an image display device for displaying information. The image display device may include a thin display device such as a liquid crystal display or an organic EL display.

通信部72は、通信インターフェースであり、上述した通信部11と同様のものであるため、詳細な説明を省略する。 The communication unit 72 is a communication interface and is similar to the communication unit 11 described above, so a detailed description thereof will be omitted.

処理部73は、例えば、1以上のプロセッサ(マイクロプロセッサ)と1以上のメモリとを含むコンピュータシステムにより実現され得る。 The processing unit 73 can be realized, for example, by a computer system including one or more processors (microprocessors) and one or more memories.

(3)被災度学習システム
被災度学習システム2は、通信部21と、記憶部22と、処理部23とを備える。被災度学習システム2は、1以上のサーバにより実現され得る。
(3) Disaster damage learning system The disaster damage learning system 2 includes a communication section 21, a storage section 22, and a processing section 23. The disaster severity learning system 2 can be realized by one or more servers.

通信部21は、通信インターフェースであり、上述した通信部11と同様のものであるため、詳細な説明を省略する。記憶部22は、処理部23が利用する情報を記憶するために用いられるものであり、上述した記憶部12と同様のものであるため、詳細な説明を省略する。処理部23は、被災度学習システム2の動作を制御する制御回路であり、上述した記憶部12と同様のハードウェア構成を有している。 The communication unit 21 is a communication interface and is similar to the communication unit 11 described above, so a detailed description thereof will be omitted. The storage unit 22 is used to store information used by the processing unit 23, and is similar to the storage unit 12 described above, so a detailed explanation will be omitted. The processing unit 23 is a control circuit that controls the operation of the disaster severity learning system 2, and has the same hardware configuration as the storage unit 12 described above.

処理部23は、前処理部231と、取得部232と、学習部233と、提示部234と、を備える。前処理部231と、取得部232と、学習部233と、提示部234とは実体のある構成を示しているわけではなく、処理部23によって実現される機能を示している。 The processing unit 23 includes a preprocessing unit 231, an acquisition unit 232, a learning unit 233, and a presentation unit 234. The preprocessing unit 231, the acquisition unit 232, the learning unit 233, and the presentation unit 234 do not represent actual configurations, but represent functions realized by the processing unit 23.

(3-1)前処理部
前処理部231は、被災度学習システム2が実行する被災度学習方法において必要となる教師データとなる入力情報を取得するため、地震動の時系列波形を基に、入力情報を求める処理を実行する。教師データの基となる地震動の時系列波形は、被災度推定システム1と同様に、通信ネットワーク8を介して地震動波形サーバ9より取得可能である。
(3-1) Pre-processing unit The pre-processing unit 231 uses the time-series waveform of the seismic motion to acquire input information that serves as teacher data required in the disaster severity learning method executed by the disaster severity learning system 2. Executes a process that requests input information. The time-series waveform of the seismic motion, which is the basis of the training data, can be obtained from the seismic motion waveform server 9 via the communication network 8, similarly to the disaster degree estimation system 1.

前処理部231は、地震動波形サーバ9より取得した地震動の時系列波形より、被災度学習システム2が実行する被災度学習方法において必要となる入力情報を求める。入力情報として、被災度推定システム1と同様のPGA、所定加速度継続時間およびスペクトル値の情報を含む。 The preprocessing unit 231 obtains input information necessary for the disaster severity learning method executed by the disaster severity learning system 2 from the time-series waveform of the earthquake motion acquired from the earthquake motion waveform server 9. The input information includes information on the PGA, predetermined acceleration duration, and spectrum value similar to the disaster damage estimation system 1.

(3-2)取得部
取得部232は、被災度学習システム2が実行する被災度学習方法において必要となる入力情報を取得する。入力情報は、上述したPGA、所定加速度継続時間およびスペクトル値の情報である。
(3-2) Acquisition Unit The acquisition unit 232 acquires input information necessary for the disaster severity learning method executed by the disaster severity learning system 2. The input information is the above-mentioned PGA, predetermined acceleration duration time, and spectrum value information.

(3-3)学習部
学習部233は、入力情報を入力データとする学習済みモデルであって、建物の地震動による被災度に関する出力情報を教師データとする学習済みモデルを機械学習により構築する。学習済みモデルは、図2に示すニューラルネットワーク4を用いたものである。本実施形態では、入力データとして、入力層41のノード41000に計測震度のデータを入力し、ノード41001にPGAのデータを入力し、ノード41002に所定加速度継続時間のデータ(正数値は0.8)を入力する。また、入力層41のノード41003~ノード41102に、地震動の加速度応答スペクトル値を入力する。スペクトル値は、ノード41003には周期帯0.0(sec)以上かつ0.1(sec)未満の最大値、ノード41004には周期帯0.1(sec)以上かつ0.2(sec)未満の最大値、・・・ノード41102には周期帯9.9(sec)以上かつ10.0(sec)未満の最大値をそれぞれ入力する。入力層41は、ノード41000~ノード41102の103個のノードを有している。
(3-3) Learning Unit The learning unit 233 uses machine learning to construct a trained model that uses input information as input data, and uses output information regarding the degree of damage caused by earthquake motion of buildings as training data. The trained model uses the neural network 4 shown in FIG. 2. In this embodiment, as input data, measured seismic intensity data is input to the node 41000 of the input layer 41, PGA data is input to the node 41001, and data of the predetermined acceleration duration is input to the node 41002 (a positive value is 0.8 ). Furthermore, the acceleration response spectrum values of earthquake motion are input to nodes 41003 to 41102 of the input layer 41. The spectrum value is the maximum value of the period band 0.0 (sec) or more and less than 0.1 (sec) for node 41003, and the maximum value of the period band 0.1 (sec) or more and less than 0.2 (sec) for node 41004. The maximum value of . The input layer 41 has 103 nodes, node 41000 to node 41102.

ニューラルネットワーク4は、中間層42及び中間層43の2層の中間層を有している。中間層42は、ノード4200~ノード4263の64個のノードを有している。中間層43は、ノード4300~ノード4363の64個のノードを有している。出力層44は、1個のノード440を有している。ノード440から出力される出力情報は、建物の最大変形量の情報である。教師データとしての建物の最大変形量は、建物の解析モデルについて、前処理部231で取得した地震動の時系列波形を時刻歴応答解析して得たデータを用いる。建物の解析モデルについては後述する。 The neural network 4 has two intermediate layers, an intermediate layer 42 and an intermediate layer 43. The middle layer 42 has 64 nodes, node 4200 to node 4263. The intermediate layer 43 has 64 nodes, node 4300 to node 4363. Output layer 44 has one node 440. The output information output from the node 440 is information on the maximum amount of deformation of the building. The maximum amount of deformation of the building as training data uses data obtained by time-history response analysis of the time-series waveform of seismic motion acquired by the preprocessing unit 231 for the analytical model of the building. The analytical model of the building will be described later.

(3-4)提示部
図1に示すように、提示部234は、処理部23からの各種の出力結果を提示する。本実施形態では、処理部23からの各種の出力結果を、通信ネットワーク8を通じて端末装置7に送ることで、端末装置7の入出力部71にて出力情報を提示可能である。
(3-4) Presentation Unit As shown in FIG. 1, the presentation unit 234 presents various output results from the processing unit 23. In this embodiment, by sending various output results from the processing section 23 to the terminal device 7 through the communication network 8, the output information can be presented at the input/output section 71 of the terminal device 7.

(4)建物の解析モデル
図3に示すように、本実施形態では、建物3の解析モデルは、非線形せん断剛性と線形粘性減衰を有するn質点系せん断モデルとする。建物3への外乱入力は地盤30の振動のみとし、地盤30のロッキングバネの影響はないものとする。
(4) Analytical model of building As shown in FIG. 3, in this embodiment, the analytical model of building 3 is an n-mass point system shear model having nonlinear shear stiffness and linear viscous damping. It is assumed that the disturbance input to the building 3 is only the vibration of the ground 30, and there is no influence of the rocking spring of the ground 30.

図1に示すように、建物3のi(i=1~n)層(階)の質量をmで表わし、i層の層剛性をk、i層の層減衰係数をc、時刻tにおけるi層の絶対加速度をu″(t)で表わすと、i層における時刻tの層せん断力Q(t)は、運動方程式に基づいて下記の[数1]に示す式で表現される。式(1)のu′(t)は時刻tの絶対速度であり、u(t)は時刻tの絶対変位である。 As shown in FIG. 1, the mass of the i (i=1 to n) layer (story) of the building 3 is expressed as m i , the layer stiffness of the i layer is k i , the layer damping coefficient of the i layer is c i , and the time When the absolute acceleration of the i layer at time t is expressed as u'' i (t), the layer shear force Q i (t) at time t in the i layer is expressed by the formula shown in [Equation 1] below based on the equation of motion. In equation (1), u' i (t) is the absolute velocity at time t, and u i (t) is the absolute displacement at time t.

Figure 0007444811000001
Figure 0007444811000001

時刻tにおけるi層の層間変位D(t)=u(t)-u-1(t)、層間速度D′(t)=u′(t)-u-1′(t)、層間加速度D″(t)=u″(t)-u-1″(t)とおくと、式(1)は下記の[数2]に示す式で表現される。 Interlayer displacement D i (t)=u i (t)−u i −1(t) of i layer at time t, interlayer velocity D i ′(t)=u i ′(t)−u i −1′( t), interlayer acceleration D i ″(t)=u i ″(t)−u i −1″(t), then equation (1) is expressed by the equation shown in [Equation 2] below.

Figure 0007444811000002
Figure 0007444811000002

1(すなわちi=1)層の絶対加速度u″(t)、絶対速度u′(t)及び絶対変位u(t)は、地震動の時系列波形より求められる。そして、建物3の解析モデルの時刻歴応答解析により、全ての絶対加速度u″(t)、絶対速度u′(t)及び絶対変位u(t)、全ての層間変位D(t)、層間速度D′(t)及び層間加速度D″、各層の最大変形量が求められる。 The absolute acceleration u 1 ″(t), absolute velocity u 1 ′(t), and absolute displacement u 1 (t) of the 1st (i=1) layer are obtained from the time-series waveform of the earthquake motion. Through time history response analysis of the analytical model, all absolute accelerations u i ″(t), absolute velocities u i ′(t), and absolute displacements u i (t), all interstory displacements D i (t), and interstory speeds D i ′(t), interlayer acceleration D i ″, and the maximum deformation amount of each layer.

本実施形態では、いわゆる異なる工法であるA工法とB工法の2パターンについて、それぞれ建物3の解析モデルを設定している。 In this embodiment, analytical models of the building 3 are set for two patterns of construction methods A and B, which are so-called different construction methods.

(5)動作
(5-1)建物の被災度学習方法
図4に示すように、被災度学習システム2において、まず、前処理ステップが実行される(S1)。ステップS1において、前処理部231により、入力情報を求める処理が実行される。具体的には、ステップS1において、計測震度のデータが取得される。更に、前処理ステップにおいて、地震動の時系列波形が取得されて、PGAのデータ及び所定加速度継続時間のデータが抽出される。更に、ステップS1において、地震動の時系列波形が加速度応答スペクトル解析されて、各周期帯におけるスペクトル値が取得される。
(5) Operation (5-1) Building disaster level learning method As shown in FIG. 4, in the disaster level learning system 2, a preprocessing step is first performed (S1). In step S1, the preprocessing unit 231 executes processing to obtain input information. Specifically, in step S1, measured seismic intensity data is acquired. Furthermore, in the preprocessing step, a time-series waveform of the seismic motion is acquired, and PGA data and data of a predetermined acceleration duration are extracted. Furthermore, in step S1, the time-series waveform of the seismic motion is subjected to acceleration response spectrum analysis to obtain spectrum values in each period band.

本実施形態では、国立研究開発法人防災科学技術研究所が運営する地震動波形サーバ9より、近年日本で発生した1250の地震動の時系列波形を取得して、被災度の学習を行っている。 In this embodiment, the time-series waveforms of 1250 earthquakes that have occurred in Japan in recent years are acquired from the seismic motion waveform server 9 operated by the National Research and Development Agency, National Research Institute for Earth Science and Disaster Prevention, to learn the degree of disaster.

次に、取得ステップが実行される(S2)。ステップS2において、取得部232により、入力情報が取得される。 Next, an acquisition step is executed (S2). In step S2, the acquisition unit 232 acquires input information.

次に、学習ステップが実行される(S3)。ステップS3において、学習部233により、学習済みモデルが構築される。本実施形態では、1250の地震動の時系列波形を基にしたデータのうち、75%にあたる937の地震動の時系列波形を基にしたデータを訓練データとし、25%にあたる313の地震動の時系列波形を基にしたデータを汎化能力の確認に用いるテストデータとした。 Next, a learning step is executed (S3). In step S3, the learning unit 233 constructs a trained model. In this embodiment, among the data based on the time-series waveforms of 1250 earthquake motions, data based on the time-series waveforms of 937 earthquake motions, which is 75%, is used as training data, and the time-series waveforms of 313 earthquake motions, which are 25%, are used as training data. The data based on this was used as test data to confirm generalization ability.

次に、提示ステップが実行される(S4)。ステップS4において、提示部234により、ニューラルネットワーク4の重みW1~W3が提示される。 Next, a presentation step is executed (S4). In step S4, the presenting unit 234 presents the weights W1 to W3 of the neural network 4.

図5に、A工法の解析モデルのニューラルネットワーク4の学習済みモデルにおける重みW1の平均値-ノード関係図を示す。この結果より、ノード41000(計測震度)、ノード41001(PGA)、ノード41002(所定加速度継続時間)、ノード41003(周期帯0.0(sec)以上かつ0.1(sec)未満のスペクトル値)~ノード41007(周期帯0.4(sec)以上かつ0.5(sec)未満のスペクトル値)の値に乗じる重みW1の値が、ノード41008~ノード41102の値に乗じる重みW1の値よりも大きく、ノード41000~ノード41107の値が出力情報に与える影響が大きいことが分かる。 FIG. 5 shows an average value of the weight W1 in the trained model of the neural network 4 of the analytical model of method A--a node relationship diagram. From this result, node 41000 (measured seismic intensity), node 41001 (PGA), node 41002 (predetermined acceleration duration), node 41003 (spectral value in period band 0.0 (sec) or more and less than 0.1 (sec)) ~The value of weight W1 multiplied by the value of node 41007 (spectral value of period band 0.4 (sec) or more and less than 0.5 (sec)) is greater than the value of weight W1 multiplied by the value of nodes 41008 to 41102. It can be seen that the values of nodes 41000 to 41107 have a large influence on the output information.

同様に、図6に、B工法の解析モデルのニューラルネットワーク4の学習済みモデルにおける重みW1の平均値-ノード関係図を示す。この結果より、図5と同様に、B工法においても、ノード41008~ノード41102の値に乗じる重みW1の値よりも大きく、ノード41000~ノード41107の値が出力情報に与える影響が大きいことが分かる。 Similarly, FIG. 6 shows an average value of the weight W1 in the trained model of the neural network 4 of the analytical model of method B--a node relationship diagram. From this result, it can be seen that in Method B, the values of nodes 41000 to 41107 have a large influence on the output information, which is larger than the value of weight W1 multiplied by the values of nodes 41008 to 41102, as in FIG. 5. .

(5-2)建物の被災度推定方法
図7に示すように、被災度推定システム1において、まず、前処理ステップが実行される(S11)。ステップS11において、前処理部131により、入力情報を求める処理が実行される。具体的には、ステップS11において、計測震度のデータが取得される。更に、ステップS11において、地震動の時系列波形が取得されて、PGAのデータ及び所定加速度継続時間のデータが抽出される。更に、ステップS11において、地震動の時系列波形が加速度応答スペクトル解析されて、各周期帯におけるスペクトル値と各周期帯におけるスペクトル値の最大値が取得される。
(5-2) Method for estimating disaster degree of building As shown in FIG. 7, in the disaster degree estimation system 1, first, a preprocessing step is executed (S11). In step S11, the preprocessing unit 131 executes processing to obtain input information. Specifically, in step S11, measured seismic intensity data is acquired. Furthermore, in step S11, a time-series waveform of the seismic motion is acquired, and PGA data and data of a predetermined acceleration duration are extracted. Further, in step S11, the time-series waveform of the seismic motion is subjected to acceleration response spectrum analysis, and the spectral value in each periodic band and the maximum value of the spectral value in each periodic band are obtained.

次に、取得ステップが実行される(S12)。ステップS12において、取得部132により、入力情報が取得される。 Next, an acquisition step is executed (S12). In step S12, the acquisition unit 132 acquires input information.

次に、推定ステップが実行される(S13)。ステップS13において、推定部133により、学習済みモデルを用いて、被災度の推定が実行される。本実施形態では、937の地震動の時系列波形を基にした訓練データと、313の地震動の時系列波形を基にしたテストデータとについて、推定した。 Next, an estimation step is executed (S13). In step S13, the estimation unit 133 estimates the degree of disaster using the trained model. In this embodiment, estimation was made for training data based on time-series waveforms of 937 earthquake motions and test data based on time-series waveforms of 313 earthquake motions.

次に、提示ステップが実行される(S14)。ステップS14において、提示部234により、出力情報として、建物の最大変形量が提示される。 Next, a presentation step is executed (S14). In step S14, the presentation unit 234 presents the maximum amount of deformation of the building as output information.

図8Aに、A工法の訓練データを用いた被災度の推定値-正解値(教師データ)の関係図を示し、図8Bに、A工法のテストデータを用いた被災度の推定値-正解値の関係図を示す。なお、単位は任意単位(a.u.)とする。この結果より、テストデータを用いた被災度の推定値の正解度は、訓練データを用いた被災度の推定値の正解度に比べて遜色なく、被災度推定方法における汎化能力が確認できる。 Figure 8A shows a relationship diagram between the estimated value of the degree of disaster using the training data of method A - the correct value (teacher data), and Figure 8B shows the estimated value of the degree of disaster using the test data of method A - the correct value. A relationship diagram is shown. Note that the unit is an arbitrary unit (au). From this result, the accuracy of the disaster damage estimation value using the test data is comparable to the accuracy of the disaster damage estimation value using the training data, confirming the generalization ability of the disaster damage estimation method.

図9Aに、B工法の訓練データを用いた被災度の推定値-正解値の関係図を示し、図9Bに、B工法のテストデータを用いた被災度の推定値-正解値の関係図を示す。この結果からも同様に、テストデータを用いた被災度の推定値の正解度は、訓練データを用いた被災度の推定値の正解度に比べて遜色なく、被災度推定方法における汎化能力が確認できる。 Figure 9A shows the relationship diagram between the estimated value of the disaster degree using the training data of method B and the correct value, and Figure 9B shows the relationship diagram between the estimated value of the degree of disaster and the correct value using the test data of method B. show. This result also shows that the accuracy of disaster damage estimation using test data is comparable to the accuracy of disaster damage estimation using training data, and the generalization ability of the disaster damage estimation method is You can check it.

(6)本実施形態のまとめ
上述した被災度推定システム1、被災度推定方法、被災度学習システム2、被災度学習方法にあっては、入力情報として、地震動の時系列波形の情報を不要とし、学習済みモデルを用いることにより、地震動の最大加速度および地震動の時系列波形の加速度応答スペクトル解析におけるスペクトル値の情報を要するだけですむようにしている。これにより、従来のような時刻歴応答解析を行う必要がなくなり、被災度の推定に要する時間を短縮しやすい。具体的には、時刻歴応答解析を行う場合、解析に要する時間が膨大となっていたところ、推定部133における推定に要する時間は数分ですむ。なお、加速度応答スペクトル解析を行う前処理部131における処理に幾分かの時間を要するが、これを合わせても被災度推定方法に要する時間は時刻歴応答解析を行う場合と比較して10分の1以下ですみ、従来よりも大幅な時間短縮が可能である。
(6) Summary of the present embodiment In the above-described disaster degree estimation system 1, disaster degree estimation method, disaster degree learning system 2, and disaster degree learning method, information on time-series waveforms of earthquake motion is unnecessary as input information. By using a trained model, we only need information on the maximum acceleration of earthquake motion and the spectrum value in the acceleration response spectrum analysis of the time-series waveform of earthquake motion. This eliminates the need to perform time history response analysis as in the past, making it easier to shorten the time required to estimate the degree of disaster. Specifically, when performing time history response analysis, the time required for the analysis would be enormous, but the time required for estimation in the estimator 133 is now only a few minutes. Note that although some time is required for the processing in the preprocessing unit 131 that performs acceleration response spectrum analysis, the time required for the disaster degree estimation method is 10 minutes in total, compared to the time history response analysis. 1 or less, and it is possible to significantly shorten the time compared to the conventional method.

また、建物に地震センサを全く又はほとんど設置する必要がなくなり、コストの低減が図られる。 Furthermore, there is no need to install any or almost no earthquake sensors in the building, resulting in cost reduction.

(6-1)例1
図10Aに、同上の実施形態のB工法のテストデータを用いた被災度の推定値-正解値の関係図(図9Bと同)を示し、図10Bに、例1のB工法のテストデータを用いた被災度の推定値-正解値の関係図を示す。
(6-1) Example 1
FIG. 10A shows a relationship diagram (same as FIG. 9B) of the estimated degree of disaster using the test data of Method B of the above embodiment (same as FIG. 9B), and FIG. 10B shows the test data of Method B of Example 1. A relationship diagram between the estimated value of the degree of disaster used and the correct value is shown.

例1は、入力データとして、所定加速度継続時間のデータを含めない場合の例である。図10Bに示す例1の推定値の正解値に対する一致度は、図10Aに示す本実施形態の推定値の正解値に対する一致度よりも低く、図10Bに示す例1の被災度の推定値の正解度は、図10Aに示す本実施形態の被災度の推定値の正解度に比べて劣っている。このことから、入力データとして所定加速度継続時間のデータを含めることの有効性が確認できる。 Example 1 is an example in which data on the predetermined acceleration duration time is not included as input data. The degree of agreement between the estimated value of Example 1 shown in FIG. 10B and the correct value is lower than the degree of agreement between the estimated value of this embodiment and the correct value shown in FIG. 10A, and the degree of agreement between the estimated value of Example 1 shown in FIG. The degree of accuracy is inferior to the degree of accuracy of the estimated value of the degree of disaster according to the present embodiment shown in FIG. 10A. From this, the effectiveness of including data on the predetermined acceleration duration time as input data can be confirmed.

(6-2)例2
図11Aに、同上の実施形態のB工法のテストデータを用いた被災度の推定値-正解値の関係図(図9Bと同)を示し、図11Bに、例2のB工法のテストデータを用いた被災度の推定値-正解値の関係図を示す。
(6-2) Example 2
FIG. 11A shows a relationship diagram (same as FIG. 9B) of the estimated degree of disaster using the test data of Method B of the above embodiment and the correct value, and FIG. 11B shows the test data of Method B of Example 2. A relationship diagram between the estimated value of the degree of disaster used and the correct value is shown.

例2は、入力データとして、周期帯3.0(sec)以上かつ10.0(sec)未満におけるスペクトル値の最大値のデータを含めない場合の例である。すなわち、本実施形態では、入力層41がノード41000~ノード41102の103個のノードを有するのに対して、例2では、入力層41がノード41000~ノード41032の33個のノードのみを有する。 Example 2 is an example in which the input data does not include data of the maximum value of the spectrum value in a period band of 3.0 (sec) or more and less than 10.0 (sec). That is, in this embodiment, the input layer 41 has 103 nodes from node 41000 to node 41102, whereas in Example 2, the input layer 41 has only 33 nodes from node 41000 to node 41032.

図11Bに示す例2の推定値の正解値に対する一致度は、図11Aに示す本実施形態の推定値の正解値に対する一致度よりも低く、図11Bに示す例2の被災度の推定値の正解度は、図11Aに示す本実施形態の被災度の推定値の正解度に比べて劣っている。このことから、入力データのスペクトル値の周期帯0.0(sec)以上かつ10.0(sec)未満の範囲とすることの有効性が確認できる。 The degree of agreement between the estimated value of Example 2 shown in FIG. 11B and the correct value is lower than the degree of agreement between the estimated value of this embodiment and the correct value shown in FIG. 11A, and the degree of agreement between the estimated value of Example 2 shown in FIG. The degree of accuracy is inferior to the degree of accuracy of the estimated value of the degree of disaster according to the present embodiment shown in FIG. 11A. From this, it can be confirmed that it is effective to set the periodic band of the input data spectrum value to a range of 0.0 (sec) or more and less than 10.0 (sec).

(7)変形例
本開示の実施形態は、上記実施形態に限定されない。上記実施形態は、本開示の課題を達成できれば、設計等に応じて種々の変更が可能である。以下に、上記実施形態の変形例を列挙する。以下に説明する変形例は、適宜組み合わせて適用可能である。
(7) Modifications Embodiments of the present disclosure are not limited to the above embodiments. The embodiments described above can be modified in various ways depending on the design, etc., as long as the objects of the present disclosure can be achieved. Modifications of the above embodiment are listed below. The modified examples described below can be applied in combination as appropriate.

本開示における被災度推定システム1、被災度学習システム2は、例えば、コンピュータシステムを含んでいる。コンピュータシステムは、ハードウェアとしてのプロセッサ及びメモリを主構成とする。コンピュータシステムのメモリに記録されたプログラムをプロセッサが実行することによって、本開示における建物の被災度推定方法や被災度学習方法としての機能が実現される。プログラムは、コンピュータシステムのメモリに予め記録されてもよく、電気通信回線を通じて提供されてもよく、コンピュータシステムで読み取り可能なメモリカード、光学ディスク、ハードディスクドライブ等の非一時的記録媒体に記録されて提供されてもよい。コンピュータシステムのプロセッサは、半導体集積回路(IC)又は大規模集積回路(LSI)を含む1ないし複数の電子回路で構成される。ここでいうIC又はLSI等の集積回路は、集積の度合いによって呼び方が異なっており、システムLSI、VLSI(Very Large Scale Integration)、又はULSI(Ultra Large Scale Integration)と呼ばれる集積回路を含む。さらに、LSIの製造後にプログラムされる、FPGA(Field- Programmable Gate Array)、又はLSI内部の接合関係の再構成若しくはLSI内部の回路区画の再構成が可能な論理デバイスについても、プロセッサとして採用することができる。複数の電子回路は、1つのチップに集約されていてもよいし、複数のチップに分散して設けられていてもよい。複数のチップは、1つの装置に集約されていてもよいし、複数の装置に分散して設けられていてもよい。ここでいうコンピュータシステムは、1以上のプロセッサ及び1以上のメモリを有するマイクロコントローラを含む。したがって、マイクロコントローラについても、半導体集積回路又は大規模集積回路を含む1ないし複数の電子回路で構成される。 The disaster damage estimation system 1 and the disaster damage learning system 2 in the present disclosure include, for example, a computer system. A computer system mainly consists of a processor and a memory as hardware. When a processor executes a program recorded in the memory of the computer system, the functions of the building damage estimation method and disaster damage learning method in the present disclosure are realized. The program may be pre-recorded in the memory of the computer system, may be provided through a telecommunications line, or may be recorded on a non-transitory storage medium readable by the computer system, such as a memory card, optical disc, hard disk drive, etc. may be provided. A processor in a computer system is comprised of one or more electronic circuits including semiconductor integrated circuits (ICs) or large scale integrated circuits (LSIs). The integrated circuits such as IC or LSI referred to herein have different names depending on the degree of integration, and include integrated circuits called system LSI, VLSI (Very Large Scale Integration), or ULSI (Ultra Large Scale Integration). Furthermore, FPGAs (Field-Programmable Gate Arrays), which are programmed after the LSI is manufactured, or logic devices that can reconfigure the connections inside the LSI or reconfigure the circuit sections inside the LSI, may also be used as processors. I can do it. The plurality of electronic circuits may be integrated into one chip, or may be provided in a distributed manner over a plurality of chips. A plurality of chips may be integrated into one device, or may be distributed and provided in a plurality of devices. The computer system herein includes a microcontroller having one or more processors and one or more memories. Therefore, the microcontroller is also composed of one or more electronic circuits including semiconductor integrated circuits or large-scale integrated circuits.

また、被災度推定システム1、被災度学習システム2における複数の機能が、1つの筐体内に集約されてもよい。さらに、被災度推定システム1、被災度学習システム2の少なくとも一部の機能、例えば、処理部23の一部の機能がいわゆるクラウド(クラウドコンピューティング)等によって実現されてもよい。 Further, a plurality of functions in the disaster damage estimation system 1 and the disaster damage learning system 2 may be integrated into one housing. Furthermore, at least some functions of the disaster damage estimation system 1 and the disaster damage learning system 2, for example, some functions of the processing unit 23, may be realized by a so-called cloud (cloud computing) or the like.

地震動の時系列波形として、公的機関が公表する地震動の時系列波形を利用しているが、自ら計測して地震動の時系列波形を取得してもよい。 As the time-series waveform of seismic motion, the time-series waveform of seismic motion published by a public institution is used, but the time-series waveform of seismic motion may be obtained by measuring it by yourself.

地震動の時系列波形は、必ずしも通信ネットワーク8を介して取得しなくてもよい。例えば、通信ネットワーク8を介さず、被災度推定システム1が地震動波形サーバ9より直接的に地震動の時系列波形を取得してもよい。また、被災度推定システム1が地震計を直接備えてもよい。 The time-series waveform of seismic motion does not necessarily need to be acquired via the communication network 8. For example, the disaster degree estimation system 1 may directly acquire the time-series waveform of the earthquake motion from the earthquake motion waveform server 9 without going through the communication network 8 . Moreover, the disaster degree estimation system 1 may be directly equipped with a seismometer.

PGAに乗じる正数値は、0.8でなくてもよく、0.9、0.7、0.6、0.5等、適宜に設定可能である。 The positive value by which PGA is multiplied does not have to be 0.8, and can be set as appropriate, such as 0.9, 0.7, 0.6, 0.5, etc.

推定部133からの出力情報は、建物の最大変形量の情報に限定されない。 The output information from the estimation unit 133 is not limited to information on the maximum amount of deformation of the building.

機械学習のアルゴリズムとしては、教師あり学習のアルゴリズムであるサポートベクトルマシン(Support Vector Machine)等、ニューラルネットワークを用いるもの以外のアルゴリズムを用いてもよい。 As a machine learning algorithm, an algorithm other than one using a neural network may be used, such as a supervised learning algorithm such as a support vector machine.

建物3の解析モデルは、非線形せん断剛性と線形粘性減衰を有するn質点系せん断モデルに限定されない。 The analytical model of the building 3 is not limited to an n-mass point system shear model having nonlinear shear stiffness and linear viscous damping.

推定部133における推定に用いられる入力情報は、計測震度の情報を含むことが好ましい。これにより、建物の被災度の推定の精度を高めやすくなる。 It is preferable that the input information used for estimation in the estimation unit 133 includes information on measured seismic intensity. This makes it easier to improve the accuracy of estimating the degree of damage to a building.

推定部133において推定される出力情報は、複数の層についての最大変形量の情報を含むことが好ましい。これにより、建物の被災度の推定の精度を高めやすくなる。 The output information estimated by the estimation unit 133 preferably includes information on the maximum deformation amount for a plurality of layers. This makes it easier to improve the accuracy of estimating the degree of damage to a building.

推定部133において推定される出力情報は、最大層間変形量の情報を含むことが好ましい。これにより、建物の被災度の推定の精度を高めやすくなる。 It is preferable that the output information estimated by the estimation unit 133 includes information on the maximum amount of interlayer deformation. This makes it easier to improve the accuracy of estimating the degree of damage to a building.

推定部133において推定される出力情報は、建物の累積塑性変形量の情報を含むことが好ましい。これにより、建物の被災度の推定の精度を高めやすくなる。 The output information estimated by the estimation unit 133 preferably includes information on the amount of cumulative plastic deformation of the building. This makes it easier to improve the accuracy of estimating the degree of damage to a building.

推定部133において推定される出力情報は、被災ランクの情報を含むことが好ましい。これにより、建物の被災度の推定の精度を高めやすくなる。 The output information estimated by the estimation unit 133 preferably includes information on the disaster rank. This makes it easier to improve the accuracy of estimating the degree of damage to a building.

応答スペクトル解析としては、加速度応答スペクトル解析の他に、速度応答スペクトル解析、変位応答スペクトル解析等が挙げられるが、応答スペクトル解析としてはこれらのうち任意のものが採用可能であって、応答スペクトル解析は特に限定されない。 In addition to acceleration response spectrum analysis, response spectrum analysis includes velocity response spectrum analysis, displacement response spectrum analysis, etc. Any of these can be adopted as response spectrum analysis, and response spectrum analysis is not particularly limited.

(8)態様
上記実施形態及び変形例から明らかなように、本開示は、下記の態様を含む。
(8) Aspects As is clear from the above embodiments and modifications, the present disclosure includes the following aspects.

第1の態様は、建物の被災度推定方法であって、取得ステップと、推定ステップと、を備える。取得ステップは、地震動の最大加速度および地震動の時系列波形の応答スペクトル解析におけるスペクトル値の情報を含む入力情報を取得するステップである。推定ステップは、予め構築された学習済みモデルを用いて、入力情報を基に出力情報を推定するステップである。学習済みモデルは、入力情報を入力データとするものであって、建物の地震動による被災度に関する出力情報を教師データとして、機械学習により構築される。この態様によれば、入力情報として、地震動の時系列波形の情報を不要とし、学習済みモデルを用いることにより、地震動の最大加速度および地震動の時系列波形の応答スペクトル解析におけるスペクトル値の情報を要するだけである。このため、従来のような時刻歴応答解析を行う必要がなくなり、被災度の推定に要する時間を短縮しやすい。また、建物に地震センサを全く又はほとんど設置する必要がなくなり、コストの低減をしやすい。 The first aspect is a method for estimating the degree of damage to a building, and includes an acquisition step and an estimation step. The acquisition step is a step of acquiring input information including maximum acceleration of earthquake motion and information on spectrum values in response spectrum analysis of time-series waveforms of earthquake motion. The estimation step is a step of estimating output information based on input information using a learned model built in advance. The learned model uses input information as input data, and is constructed by machine learning using output information regarding the degree of damage caused by earthquake motion of a building as training data. According to this aspect, information on the time-series waveform of the earthquake motion is not required as input information, and information on the maximum acceleration of the earthquake motion and the spectrum value in the response spectrum analysis of the time-series waveform of the earthquake motion is required by using the trained model. Only. Therefore, there is no need to perform time history response analysis as in the past, and the time required for estimating the degree of disaster can be easily shortened. Furthermore, there is no need to install any or almost no earthquake sensors in the building, making it easier to reduce costs.

第2の態様は、第1の態様に基づく建物の被災度推定方法である。第2の態様では、入力情報は、所定加速度継続時間の情報を含む。所定加速度継続時間は、地震動の所定の継続時間に占める、地震動の最大加速度に1未満の正数値を乗じてなる所定加速度を超える地震動の加速度が継続する時間である。この態様によれば、建物の被災度の推定の精度を高めやすくなる。 The second aspect is a method for estimating the degree of damage to a building based on the first aspect. In the second aspect, the input information includes information on the predetermined acceleration duration. The predetermined acceleration duration is the time during which the acceleration of the seismic motion exceeds the predetermined acceleration obtained by multiplying the maximum acceleration of the seismic motion by a positive value less than 1, which occupies the predetermined duration of the seismic motion. According to this aspect, it becomes easier to improve the accuracy of estimating the degree of damage to a building.

第3の態様は、第1又は第2の態様に基づく建物の被災度推定方法である。第3の態様では、入力情報は、計測震度の情報を含む。この態様によれば、建物の被災度の推定の精度を高めやすくなる。 The third aspect is a method for estimating the degree of damage to a building based on the first or second aspect. In the third aspect, the input information includes information on measured seismic intensity. According to this aspect, it becomes easier to improve the accuracy of estimating the degree of damage to a building.

第4の態様は、第1~第3のいずれかの態様に基づく建物の被災度推定方法である。第4の態様では、出力情報は、建物の最大変形量の情報を含む。この態様によれば、建物の被災度の推定の精度を高めやすくなる。 A fourth aspect is a method for estimating the degree of damage to a building based on any one of the first to third aspects. In the fourth aspect, the output information includes information on the maximum amount of deformation of the building. According to this aspect, it becomes easier to improve the accuracy of estimating the degree of damage to a building.

第5の態様は、第4の態様に基づく建物の被災度推定方法である。第5の態様では、建物は、複数の層を備える。出力情報は、複数の層についての最大変形量の情報を含む。この態様によれば、建物の被災度の推定の精度を高めやすくなる。 A fifth aspect is a method for estimating the degree of damage to a building based on the fourth aspect. In a fifth aspect, the building comprises multiple tiers. The output information includes information on maximum deformation amounts for multiple layers. According to this aspect, it becomes easier to improve the accuracy of estimating the degree of damage to a building.

第6の態様は、第1~第5のいずれかの態様に基づく建物の被災度推定方法である。第6の態様では、建物は、複数の層を備える。出力情報は、最大層間変形量の情報を含む。この態様によれば、建物の被災度の推定の精度を高めやすくなる。 A sixth aspect is a method for estimating the degree of damage to a building based on any one of the first to fifth aspects. In a sixth aspect, the building comprises multiple tiers. The output information includes information on the maximum amount of interlayer deformation. According to this aspect, it becomes easier to improve the accuracy of estimating the degree of damage to a building.

第7の態様は、第1~第6のいずれかの態様に基づく建物の被災度推定方法である。第7の態様では、出力情報は、建物の累積塑性変形量の情報を含む。この態様によれば、建物の被災度の推定の精度を高めやすくなる。 A seventh aspect is a method for estimating the degree of damage to a building based on any one of the first to sixth aspects. In the seventh aspect, the output information includes information on the amount of cumulative plastic deformation of the building. According to this aspect, it becomes easier to improve the accuracy of estimating the degree of damage to a building.

第8の態様は、第1~第7のいずれかの態様に基づく建物の被災度推定方法である。第8の態様では、出力情報は、被災ランクの情報を含む。この態様によれば、建物の被災度の推定の精度を高めやすくなる。 An eighth aspect is a method for estimating the degree of damage to a building based on any one of the first to seventh aspects. In the eighth aspect, the output information includes information on disaster rank. According to this aspect, it becomes easier to improve the accuracy of estimating the degree of damage to a building.

第9の態様は、第1~第8のいずれかの態様に基づく建物の被災度推定方法である。第9の態様では、学習済みモデルは、ニューラルネットワーク4を用いたものである。この態様によれば、建物の被災度の推定の精度を高めやすくなる。 A ninth aspect is a method for estimating the degree of damage to a building based on any one of the first to eighth aspects. In the ninth aspect, the learned model uses the neural network 4. According to this aspect, it becomes easier to improve the accuracy of estimating the degree of damage to a building.

第10の態様は、建物の被災度推定システム1であって、取得部132と、推定部133と、を備える。取得部132は、地震動の最大加速度および地震動の時系列波形の応答スペクトル解析におけるスペクトル値の情報を含む入力情報を取得する。推定部133は、予め構築された学習済みモデルを用いて、入力情報を基に出力情報を推定する。学習済みモデルは、入力情報を入力データとするものであって、建物の地震動による被災度に関する出力情報を教師データとして、機械学習により構築される。この態様によれば、入力情報として、地震動の時系列波形の情報を不要とし、学習済みモデルを用いることにより、地震動の最大加速度および地震動の時系列波形の応答スペクトル解析におけるスペクトル値の情報を要するだけである。このため、従来のような時刻歴応答解析を行う必要がなくなり、被災度の推定に要する時間を短縮しやすい。また、建物に地震センサを全く又はほとんど設置する必要がなくなり、コストの低減をしやすい。 A tenth aspect is a building damage degree estimation system 1, which includes an acquisition section 132 and an estimation section 133. The acquisition unit 132 acquires input information including information on maximum acceleration of earthquake motion and spectrum values in response spectrum analysis of time-series waveform of earthquake motion. The estimation unit 133 estimates output information based on input information using a learned model built in advance. The learned model uses input information as input data, and is constructed by machine learning using output information regarding the degree of damage caused by earthquake motion of a building as training data. According to this aspect, information on the time-series waveform of the earthquake motion is not required as input information, and information on the maximum acceleration of the earthquake motion and the spectrum value in the response spectrum analysis of the time-series waveform of the earthquake motion is required by using the trained model. Only. Therefore, there is no need to perform time history response analysis as in the past, and the time required for estimating the degree of disaster can be easily shortened. Furthermore, there is no need to install any or almost no earthquake sensors in the building, making it easier to reduce costs.

第11の態様は、建物の被災度学習方法であって、取得ステップと、学習ステップと、を備える。取得ステップは、地震動の最大加速度および地震動の時系列波形の応答スペクトル解析におけるスペクトル値の情報を含む入力情報を取得するステップである。学習ステップは、入力情報を入力データとする学習済みモデルであって、建物の地震動による被災度に関する出力情報を教師データとする学習済みモデルを機械学習により構築するステップである。この態様によれば、入力情報として、地震動の時系列波形の情報を不要とし、学習済みモデルを用いることにより、地震動の最大加速度および地震動の時系列波形の応答スペクトル解析におけるスペクトル値の情報を要するだけである。このため、従来のような時刻歴応答解析を行う必要がなくなり、被災度の推定に要する時間を短縮しやすい。また、建物に地震センサを全く又はほとんど設置する必要がなくなり、コストの低減をしやすい。 An eleventh aspect is a method for learning the degree of damage to a building, which includes an acquisition step and a learning step. The acquisition step is a step of acquiring input information including maximum acceleration of earthquake motion and information on spectrum values in response spectrum analysis of time-series waveforms of earthquake motion. The learning step is a step of constructing, by machine learning, a trained model that uses input information as input data and uses output information regarding the degree of damage caused by earthquake motion of buildings as training data. According to this aspect, information on the time-series waveform of the earthquake motion is not required as input information, and information on the maximum acceleration of the earthquake motion and the spectrum value in the response spectrum analysis of the time-series waveform of the earthquake motion is required by using the trained model. Only. Therefore, there is no need to perform time history response analysis as in the past, and the time required for estimating the degree of disaster can be easily shortened. Furthermore, there is no need to install any or almost no earthquake sensors in the building, making it easier to reduce costs.

第12の態様は、建物の被災度学習システム2であって、取得部232と、学習部233と、を備える。取得部232は、地震動の最大加速度および地震動の時系列波形の応答スペクトル解析におけるスペクトル値の情報を含む入力情報を取得する。学習部233は、入力情報を入力データとする学習済みモデルであって、建物の地震動による被災度に関する出力情報を教師データとする学習済みモデルを機械学習により構築する。この態様によれば、入力情報として、地震動の時系列波形の情報を不要とし、学習済みモデルを用いることにより、地震動の最大加速度および地震動の時系列波形の応答スペクトル解析におけるスペクトル値の情報を要するだけである。このため、従来のような時刻歴応答解析を行う必要がなくなり、被災度の推定に要する時間を短縮しやすい。また、建物に地震センサを全く又はほとんど設置する必要がなくなり、コストの低減をしやすい。 A twelfth aspect is a building disaster degree learning system 2, which includes an acquisition section 232 and a learning section 233. The acquisition unit 232 acquires input information including information on maximum acceleration of earthquake motion and spectrum values in response spectrum analysis of time-series waveform of earthquake motion. The learning unit 233 constructs, by machine learning, a trained model that uses input information as input data and uses output information regarding the degree of damage caused by earthquake motion of buildings as training data. According to this aspect, information on the time-series waveform of seismic motion is not required as input information, and information on the maximum acceleration of seismic motion and the spectrum value in response spectrum analysis of the time-series waveform of seismic motion is required by using a trained model. Only. Therefore, there is no need to perform time history response analysis as in the past, and the time required to estimate the degree of disaster can be easily shortened. Furthermore, there is no need to install any or almost no earthquake sensors in the building, making it easier to reduce costs.

第13の態様は、プログラムであって、第1~第9の態様のいずれかの態様の建物の被災度推定方法又は第11の態様の建物の被災度学習方法を、1以上のプロセッサに実行させる。この態様によれば、入力情報として、地震動の時系列波形の情報を不要とし、学習済みモデルを用いることにより、地震動の最大加速度および地震動の時系列波形の応答スペクトル解析におけるスペクトル値の情報を要するだけである。このため、従来のような時刻歴応答解析を行う必要がなくなり、被災度の推定に要する時間を短縮しやすい。また、建物に地震センサを全く又はほとんど設置する必要がなくなり、コストの低減をしやすい。 A thirteenth aspect is a program that executes the building damage estimation method according to any one of the first to ninth aspects or the building disaster learning method according to the eleventh aspect on one or more processors. let According to this aspect, information on the time-series waveform of the earthquake motion is not required as input information, and information on the maximum acceleration of the earthquake motion and the spectrum value in the response spectrum analysis of the time-series waveform of the earthquake motion is required by using the trained model. Only. Therefore, there is no need to perform time history response analysis as in the past, and the time required for estimating the degree of disaster can be easily shortened. Furthermore, there is no need to install any or almost no earthquake sensors in the building, making it easier to reduce costs.

1 被災度推定システム
132 取得部
133 推定部
2 被災度学習システム
232 取得部
233 学習部
3 建物
4 ニューラルネットワーク
1 Disaster damage estimation system 132 Acquisition unit 133 Estimation unit 2 Disaster damage learning system 232 Acquisition unit 233 Learning unit 3 Building 4 Neural network

Claims (13)

地震動の最大加速度および前記地震動の時系列波形の応答スペクトル解析におけるスペクトル値の周期帯0.0(sec)以上かつ10.0(sec)未満の範囲の情報を含む入力情報を取得する取得ステップと、
前記入力情報を入力データとする学習済みモデルであって、建物の前記地震動による被災度に関する出力情報を教師データとして機械学習により予め構築された前記学習済みモデルを用いて、前記入力情報を基に前記出力情報を推定する推定ステップと、を備え
前記スペクトル値の周期帯0.0(sec)以上かつ10.0(sec)未満の範囲の情報は、周期帯0.0(sec)以上かつ0.1(sec)未満の最大値、周期帯0.1(sec)以上かつ0.2(sec)未満の最大値、・・・・周期帯9.9(sec)以上かつ10.0(sec)未満の最大値の各情報を含む、
建物の被災度推定方法。
an acquisition step of acquiring input information including maximum acceleration of earthquake motion and information on a periodic band of spectral values in a range of 0.0 (sec) or more and less than 10.0 (sec) in response spectrum analysis of the time-series waveform of the earthquake motion; ,
A trained model that uses the input information as input data, and that is constructed in advance by machine learning using output information regarding the degree of damage caused by the earthquake motion of a building as training data, based on the input information. an estimation step of estimating the output information ,
The information on the range of the periodic band of 0.0 (sec) or more and less than 10.0 (sec) of the spectrum value is the maximum value of the periodic band of 0.0 (sec) or more and less than 0.1 (sec), the periodic band A maximum value of 0.1 (sec) or more and less than 0.2 (sec), ... Contains each information of a maximum value of a period band of 9.9 (sec) or more and less than 10.0 (sec),
A method for estimating the degree of damage to buildings.
前記入力情報は、所定加速度継続時間の情報を含み、
前記所定加速度継続時間は、前記地震動の所定の継続時間に占める、前記地震動の最大加速度に1未満の正数値を乗じてなる所定加速度を超える前記地震動の加速度が継続する時間である、
請求項1記載の建物の被災度推定方法。
The input information includes information on a predetermined acceleration duration,
The predetermined acceleration duration is the time during which the acceleration of the seismic motion exceeds a predetermined acceleration obtained by multiplying the maximum acceleration of the seismic motion by a positive value less than 1, which occupies the predetermined duration of the seismic motion.
The method for estimating the degree of damage to a building according to claim 1.
前記入力情報は、計測震度の情報を含む、
請求項1又は2記載の建物の被災度推定方法。
The input information includes information on measured seismic intensity;
The method for estimating the degree of damage to a building according to claim 1 or 2.
前記出力情報は、前記建物の最大変形量の情報を含む、
請求項1~3のいずれか一項に記載の建物の被災度推定方法。
The output information includes information on the maximum amount of deformation of the building.
The method for estimating the degree of damage to a building according to any one of claims 1 to 3.
前記建物は、複数の層を備え、
前記出力情報は、複数の前記層についての前記最大変形量の情報を含む、
請求項4に記載の建物の被災度推定方法。
The building comprises a plurality of layers,
the output information includes information on the maximum deformation amount for the plurality of layers;
The method for estimating the degree of damage to a building according to claim 4.
前記建物は、複数の層を備え、
前記出力情報は、最大層間変形量の情報を含む、
請求項1~5のいずれか一項に記載の建物の被災度推定方法。
The building comprises a plurality of layers,
The output information includes information on the maximum interlayer deformation amount.
The method for estimating the degree of damage to a building according to any one of claims 1 to 5.
前記出力情報は、前記建物の累積塑性変形量の情報を含む、
請求項1~6のいずれか一項に記載の建物の被災度推定方法。
The output information includes information on the amount of cumulative plastic deformation of the building.
The method for estimating the degree of damage to a building according to any one of claims 1 to 6.
前記出力情報は、被災ランクの情報を含む、
請求項1~7のいずれか一項に記載の建物の被災度推定方法。
The output information includes information on disaster rank.
The method for estimating the degree of damage to a building according to any one of claims 1 to 7.
前記学習済みモデルは、ニューラルネットワークを用いたものである、
請求項1~8のいずれか一項に記載の建物の被災度推定方法。
The trained model uses a neural network,
The method for estimating the degree of damage to a building according to any one of claims 1 to 8.
地震動の最大加速度および前記地震動の時系列波形の応答スペクトル解析におけるスペクトル値の周期帯0.0(sec)以上かつ10.0(sec)未満の範囲の情報を含む入力情報を取得する取得部と、
前記入力情報を入力データとする学習済みモデルであって、建物の前記地震動による被災度に関する出力情報を教師データとして機械学習により予め構築された前記学習済みモデルを用いて、前記入力情報を基に前記出力情報を推定する推定部と、を備え
前記スペクトル値の周期帯0.0(sec)以上かつ10.0(sec)未満の範囲の情報は、周期帯0.0(sec)以上かつ0.1(sec)未満の最大値、周期帯0.1(sec)以上かつ0.2(sec)未満の最大値、・・・・周期帯9.9(sec)以上かつ10.0(sec)未満の最大値の各情報を含む、
建物の被災度推定システム。
an acquisition unit that acquires input information including maximum acceleration of earthquake motion and information on a periodic band of spectrum values in a range of 0.0 (sec) or more and less than 10.0 (sec) in response spectrum analysis of the time-series waveform of the earthquake motion; ,
A trained model that uses the input information as input data, and that is constructed in advance by machine learning using output information regarding the degree of damage caused by the earthquake motion of the building as training data, based on the input information. an estimation unit that estimates the output information ,
The information on the range of the periodic band of 0.0 (sec) or more and less than 10.0 (sec) of the spectrum value is the maximum value of the periodic band of 0.0 (sec) or more and less than 0.1 (sec), the periodic band A maximum value of 0.1 (sec) or more and less than 0.2 (sec), ... Contains each information of a maximum value of a periodic band of 9.9 (sec) or more and less than 10.0 (sec),
Building damage estimation system.
地震動の最大加速度および前記地震動の時系列波形の応答スペクトル解析におけるスペクトル値の周期帯0.0(sec)以上かつ10.0(sec)未満の範囲の情報を含む入力情報を取得する取得ステップと、
前記入力情報を入力データとする学習済みモデルであって、建物の前記地震動による被災度に関する出力情報を教師データとする前記学習済みモデルを機械学習により構築する学習ステップと、を備え
前記スペクトル値の周期帯0.0(sec)以上かつ10.0(sec)未満の範囲の情報は、周期帯0.0(sec)以上かつ0.1(sec)未満の最大値、周期帯0.1(sec)以上かつ0.2(sec)未満の最大値、・・・・周期帯9.9(sec)以上かつ10.0(sec)未満の最大値の各情報を含む、
建物の被災度学習方法。
an acquisition step of acquiring input information including maximum acceleration of earthquake motion and information on a periodic band of spectral values in a range of 0.0 (sec) or more and less than 10.0 (sec) in response spectrum analysis of the time-series waveform of the earthquake motion; ,
a learning step of constructing, by machine learning, a trained model that uses the input information as input data, and that uses output information regarding the degree of damage caused by the earthquake motion of the building as training data ;
The information on the range of the periodic band of 0.0 (sec) or more and less than 10.0 (sec) of the spectrum value is the maximum value of the periodic band of 0.0 (sec) or more and less than 0.1 (sec), the periodic band A maximum value of 0.1 (sec) or more and less than 0.2 (sec), ... Contains each information of a maximum value of a period band of 9.9 (sec) or more and less than 10.0 (sec),
A method for learning the degree of damage to buildings.
地震動の最大加速度および前記地震動の時系列波形の応答スペクトル解析におけるスペクトル値の周期帯0.0(sec)以上かつ10.0(sec)未満の範囲の情報を含む入力情報を取得する取得部と、
前記入力情報を入力データとする学習済みモデルであって、建物の前記地震動による被災度に関する出力情報を教師データとする前記学習済みモデルを機械学習により構築する学習部と、を備え
前記スペクトル値の周期帯0.0(sec)以上かつ10.0(sec)未満の範囲の情報は、周期帯0.0(sec)以上かつ0.1(sec)未満の最大値、周期帯0.1(sec)以上かつ0.2(sec)未満の最大値、・・・・周期帯9.9(sec)以上かつ10.0(sec)未満の最大値の各情報を含む、
建物の被災度学習システム。
an acquisition unit that acquires input information including maximum acceleration of earthquake motion and information on a periodic band of spectrum values in a range of 0.0 (sec) or more and less than 10.0 (sec) in response spectrum analysis of the time-series waveform of the earthquake motion; ,
a learning unit that constructs, by machine learning, a trained model that uses the input information as input data, and that uses output information regarding the degree of damage caused by the earthquake motion of a building as training data ;
The information on the range of the periodic band of 0.0 (sec) or more and less than 10.0 (sec) of the spectrum value is the maximum value of the periodic band of 0.0 (sec) or more and less than 0.1 (sec), the periodic band A maximum value of 0.1 (sec) or more and less than 0.2 (sec), ... Contains each information of a maximum value of a periodic band of 9.9 (sec) or more and less than 10.0 (sec),
Building damage level learning system.
1以上のプロセッサに、
請求項1~9のいずれか一項に記載の建物の被災度推定方法又は請求項11に記載の建物の被災度学習方法を実行させる、
プログラム。
one or more processors,
Executing the building disaster degree estimation method according to any one of claims 1 to 9 or the building disaster degree learning method according to claim 11,
program.
JP2021050741A 2021-03-24 2021-03-24 Building damage estimation method, building damage estimation system, building damage learning method, building damage learning system and program Active JP7444811B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2021050741A JP7444811B2 (en) 2021-03-24 2021-03-24 Building damage estimation method, building damage estimation system, building damage learning method, building damage learning system and program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2021050741A JP7444811B2 (en) 2021-03-24 2021-03-24 Building damage estimation method, building damage estimation system, building damage learning method, building damage learning system and program

Publications (2)

Publication Number Publication Date
JP2022148889A JP2022148889A (en) 2022-10-06
JP7444811B2 true JP7444811B2 (en) 2024-03-06

Family

ID=83462342

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021050741A Active JP7444811B2 (en) 2021-03-24 2021-03-24 Building damage estimation method, building damage estimation system, building damage learning method, building damage learning system and program

Country Status (1)

Country Link
JP (1) JP7444811B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116433991B (en) * 2023-06-14 2023-08-22 中国地质大学(武汉) Post-earthquake building damage assessment method for emergency rescue

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008298704A (en) 2007-06-04 2008-12-11 Sekisui Chem Co Ltd Damage diagnosis system of building
JP2009156743A (en) 2007-12-27 2009-07-16 Osaka Gas Co Ltd Seismoscope, emergency shutoff system, and using method of emergency shutoff system
JP2011057420A (en) 2009-09-14 2011-03-24 Mitsubishi Electric Corp Elevator earthquake damage prediction device
JP2013152197A (en) 2012-01-26 2013-08-08 Daiwa House Industry Co Ltd Evaluation device, evaluation method and evaluation program for earthquake damage loss of building
JP2014153191A (en) 2013-02-08 2014-08-25 Ohbayashi Corp Calculation system of response value of building caused by earthquake, and method for creating research table of damage to building caused by earthquake
JP2016197014A (en) 2015-04-02 2016-11-24 アズビル株式会社 Building damage intensity estimating system, and method
US20180336652A1 (en) 2017-05-16 2018-11-22 One Concern, Inc. Estimation of damage prevention with building retrofit
JP2020101941A (en) 2018-12-20 2020-07-02 旭化成ホームズ株式会社 Damage level deriving system, damage level deriving device, damage level deriving method, and program

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008298704A (en) 2007-06-04 2008-12-11 Sekisui Chem Co Ltd Damage diagnosis system of building
JP2009156743A (en) 2007-12-27 2009-07-16 Osaka Gas Co Ltd Seismoscope, emergency shutoff system, and using method of emergency shutoff system
JP2011057420A (en) 2009-09-14 2011-03-24 Mitsubishi Electric Corp Elevator earthquake damage prediction device
JP2013152197A (en) 2012-01-26 2013-08-08 Daiwa House Industry Co Ltd Evaluation device, evaluation method and evaluation program for earthquake damage loss of building
JP2014153191A (en) 2013-02-08 2014-08-25 Ohbayashi Corp Calculation system of response value of building caused by earthquake, and method for creating research table of damage to building caused by earthquake
JP2016197014A (en) 2015-04-02 2016-11-24 アズビル株式会社 Building damage intensity estimating system, and method
US20180336652A1 (en) 2017-05-16 2018-11-22 One Concern, Inc. Estimation of damage prevention with building retrofit
JP2020101941A (en) 2018-12-20 2020-07-02 旭化成ホームズ株式会社 Damage level deriving system, damage level deriving device, damage level deriving method, and program

Also Published As

Publication number Publication date
JP2022148889A (en) 2022-10-06

Similar Documents

Publication Publication Date Title
Chiou et al. NGA project strong-motion database
Lai et al. Semi‐supervised structural linear/nonlinear damage detection and characterization using sparse identification
Zhang et al. Using smartphones to detect and identify construction workers’ near-miss falls based on ANN
Rupakhety et al. Can simple pulses adequately represent near-fault ground motions?
CN105975749A (en) Bearing health assessment and prediction method and system
Kuwabara et al. A new approach to system identification and damage detection of high‐rise buildings
WO2017136234A1 (en) Motion sensing wi-fi sensor networks for continuous 3d modeling and prediction of facility responses to disturbances
JP2017044614A (en) System for highly accurate evaluation of transfer function of structure, earthquake response estimation and deterioration diagnosis thereof, and its method
US20160189176A1 (en) Systems and methods for response rate determination and offer selection
Torbol Real‐time frequency‐domain decomposition for structural health monitoring using general‐purpose graphic processing unit
Hu et al. Conditional simulation of spatially variable seismic ground motions based on evolutionary spectra
Elefante et al. Disaggregation-based response weighting scheme for seismic risk assessment of structures
JP7444811B2 (en) Building damage estimation method, building damage estimation system, building damage learning method, building damage learning system and program
MacPherson et al. A stochastic extreme sea level model for the German Baltic Sea coast
Du et al. An algorithm for selecting spatially correlated ground motions at multiple sites under scenario earthquakes
Niedzwecki et al. System identification of distributed-parameter marine riser models
Hattab et al. An adaptive modeling approach to structural health monitoring of multistory buildings
Erazo et al. High‐resolution seismic monitoring of instrumented buildings using a model‐based state observer
Aydin et al. Damage detection in Timoshenko beam structures by multilayer perceptron and radial basis function networks
Poskus et al. Damage identification for hysteretic structures using a mode decomposition method
Vahedian et al. A new method for generating aftershock records using artificial neural network
Peng et al. Data driven structural damage assessment using phase space embedding and Koopman operator under stochastic excitations
Rasouli et al. A new method for damage prognosis based on incomplete modal data via an evolutionary algorithm
Domizio et al. Experimental and numerical analysis to collapse of a framed structure subjected to seismic loading
Zhang et al. Damage detection of nonlinear structures using probability density ratio estimation

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210929

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20221025

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20221226

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20230322

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230515

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20230725

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20231023

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20231109

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: 20240130

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20240222

R150 Certificate of patent or registration of utility model

Ref document number: 7444811

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150