JP7246272B2 - Ground subsidence prediction system - Google Patents
Ground subsidence prediction system Download PDFInfo
- Publication number
- JP7246272B2 JP7246272B2 JP2019122228A JP2019122228A JP7246272B2 JP 7246272 B2 JP7246272 B2 JP 7246272B2 JP 2019122228 A JP2019122228 A JP 2019122228A JP 2019122228 A JP2019122228 A JP 2019122228A JP 7246272 B2 JP7246272 B2 JP 7246272B2
- Authority
- JP
- Japan
- Prior art keywords
- subsidence
- point
- parameter
- value
- learning
- 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
Links
Images
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Foundations (AREA)
Description
本発明は、地盤沈下を予測する地盤沈下予測システムに関する。 The present invention relates to a land subsidence prediction system for predicting land subsidence.
従来から、たとえば建設予定の未設建物の建設予定地の地盤の状態を予測することが成されている。特許文献1に示す技術では、まず、地盤推定箇所周辺の複数のデータに関して既に得られているボーリングデータから、各層の地質特性値を推定し、この地質特性値から地層特定値の等高線図を生成する。次に、生成した等高線図から地盤推定箇所の地質特性値を推定する。地盤推定箇所の地質特性値と、そこで構築される構造物の仕様値から、所定の算定式を適用して、構造物構築で地盤に生じる影響の値を算出している。
2. Description of the Related Art Conventionally, for example, prediction of the state of the ground of a planned construction site for an unconstructed building has been performed. In the technique shown in
しかしながら、特許文献1に示す地盤沈下予測システムで、地盤推定箇所(地盤沈下を推定する推定地点)の回りの等高線図を作成するには、推定地点周辺の多大な箇所のボーリング試験のデータを取得する必要がある。これに加えて、推定地点の地盤の状態がわかったとしても、推定地点の地盤沈下を正確に予測することは、容易ではない。
However, in the ground subsidence prediction system shown in
本発明は、このような点を鑑みてなされたものであり、本発明では、特定の地点を地盤沈下の推定地点として、推定地点の地盤沈下を正確に予測することができる地盤沈下予測システムを提供する。 The present invention has been made in view of these points, and the present invention provides a ground subsidence prediction system capable of accurately predicting ground subsidence at a specific point as an estimated point of ground subsidence. offer.
前記課題を鑑みて、本発明は、特定の地点を地盤沈下の推定地点として、前記推定地点の地盤沈下を予測する地盤沈下予測システムであって、複数の既設建物の建設地点のボーリング試験により得られた、前記建設地点の地表からの深度、前記深度に応じた土質を数値化して設定された土質パラメータの値、および前記深度に応じた標準貫入試験におけるN値、を含む学習用地盤情報と、前記学習用地盤情報に対応した前記既設建物の地盤沈下量に基づいて、地盤沈下のし易さに応じて設定された学習用沈下パラメータの値とを、教師データとして、任意の地点の地盤情報から前記任意の地点の沈下パラメータの値の算出を機械学習により学習した第1学習部と、所定の地域の地図データと、前記地図データの地図上の複数の既知地点において、前記既知地点ごとの位置情報、前記既知地点ごとのボーリング試験により得られた、地表からの深度、前記深度に応じた土質を数値化して設定された土質パラメータの値、および前記深度に応じた標準貫入試験におけるN値、を含む既知地点情報とが、格納された地域情報格納部と、前記地図上に前記推定地点を設定する推定地点設定部と、前記推定地点設定部で設定された前記推定地点の周辺に存在する前記既知地点ごとの前記既知地点情報を抽出する周辺情報抽出部と、前記周辺情報抽出部で抽出した前記既知地点ごとの前記既知地点情報を、前記第1学習部に入力し、前記既知地点ごとの前記沈下パラメータの値を前記第1学習部に算出させる周辺沈下パラメータ算出部と、前記周辺沈下パラメータ算出部に算出させた前記既知地点ごとの前記沈下パラメータに基づいて、前記推定地点の沈下パラメータの値を算出する沈下パラメータ推定部と、を備えることを特徴とする。 In view of the above problems, the present invention provides a ground subsidence prediction system for predicting ground subsidence at a specific point as an estimated ground subsidence point, which is obtained by boring tests at construction points of a plurality of existing buildings. Ground information for learning including the depth from the ground surface of the construction site, the value of the soil parameter set by quantifying the soil quality according to the depth, and the N value in the standard penetration test according to the depth. , based on the amount of land subsidence of the existing building corresponding to the ground information for learning, and the value of the learning subsidence parameter set according to the likelihood of ground subsidence, as teacher data, the ground at an arbitrary point a first learning unit that learns the calculation of the subsidence parameter value of the arbitrary point from the information by machine learning; map data of a predetermined area; position information, the depth from the ground surface obtained by the boring test for each known point, the value of the soil parameter set by quantifying the soil according to the depth, and the N in the standard penetration test according to the depth an estimated point setting unit for setting the estimated point on the map; and a location around the estimated point set by the estimated point setting unit. A surrounding information extraction unit for extracting the known point information for each existing known point, and the known point information for each known point extracted by the surrounding information extraction unit are inputted to the first learning unit, and the known a peripheral subsidence parameter calculation unit that causes the first learning unit to calculate the value of the subsidence parameter for each point; and a subsidence parameter estimator that calculates the value of the subsidence parameter.
本発明によれば、第1学習部で、ここまで地盤沈下した既設建物に対して、既設建物の建設地点の学習用地盤情報と、学習用沈下パラメータの値とを教師データとして、任意の地点の地盤情報から、その任意の地点の沈下パラメータの値の算出を機械学習により学習させる。この機械学習は、たとえば、後述する深層学習(ディープラーニング)の1つであるディープ・ニューラル・ネットワーク(DNN)により行ってもよく、その他のアルゴリズムを利用して機械学習を行ってもよい。 According to the present invention, in the first learning unit, for an existing building that has subsided so far, the ground information for learning of the construction site of the existing building and the value of the subsidence parameter for learning are used as teacher data, and an arbitrary point Calculation of the subsidence parameter value of the arbitrary point is learned by machine learning from the ground information. This machine learning may be performed by, for example, a deep neural network (DNN), which is one of the deep learnings described later, or may be performed using other algorithms.
ここで、この教師データの対象となる建設地点は、実際に地盤沈下が発生した地点であり、各建設地点において、地盤沈下量に対して既設建物の仕様(たとえば、建物の用途、回数、建物面積、建物の基礎の仕様を数値化した値)を数値化した値で正規化した学習用沈下パラメータの値を算出する。したがって、この沈下パラメータは、同じ仕様の建物に対して、地盤沈下のし易さを数値化したものであり、換言すると地盤沈下のリスクの高さを変数にしたものである。沈下パラメータの大きさが大きくなるに従って地盤沈下がし易くなるため、地盤沈下のリスクが高く、たとえば、地盤沈下し難い基礎構造が選定される。 Here, the construction sites targeted for this training data are the sites where land subsidence actually occurred. Area, values obtained by quantifying the specifications of the foundation of the building) are normalized by the quantified values of the learning subsidence parameters. Therefore, this subsidence parameter quantifies the likelihood of ground subsidence for buildings with the same specifications, in other words, it uses the high risk of ground subsidence as a variable. Since the larger the subsidence parameter, the easier it is for the ground to subside, the risk of ground subsidence is high, and for example, a foundation structure that is unlikely to subside is selected.
本発明では、特定の地点を地盤沈下の推定地点として、推定地点を含む所定の地域の地図データと、地図データの地図上の複数の既知地点において、既知地点ごとの既知地点情報を格納した地域情報格納部を設けている。 In the present invention, with a specific point as an estimated point of ground subsidence, map data of a predetermined area including the estimated point and known point information for each known point in a plurality of known points on the map of the map data are stored. An information storage unit is provided.
ここで、推定地点設定部で推定地点が設定されると、周辺情報抽出部で、推定地点の周辺に存在する既知地点ごとの既知地点情報が抽出される。この際、周辺情報抽出部は、推定地点から所定距離の範囲内に存在する既知地点を抽出してもよく、推定地点から近い特定の個数の既知地点を抽出してもよい。 Here, when the estimated point is set by the estimated point setting unit, the surrounding information extraction unit extracts known point information for each known point existing around the estimated point. At this time, the peripheral information extraction unit may extract known points existing within a predetermined distance from the estimated point, or may extract a specific number of known points close to the estimated point.
次に、周辺沈下パラメータ算出部で、これらの抽出した既知地点の既知地点情報を第1学習部に入力し、推定地点の周辺に存在する既知地点の沈下パラメータの値を算出する。この算出した沈下パラメータの値は、推定地点の沈下パラメータの値に影響のある値である。したがって、次に、沈下パラメータ推定部で、推定地点の周辺に存在する既知地点ごとの沈下パラメータに基づいて、推定地点の沈下パラメータの値を算出する。 Next, the surrounding subsidence parameter calculation unit inputs the known point information of these extracted known points to the first learning unit, and calculates subsidence parameter values of known points existing around the estimated point. This calculated subsidence parameter value is a value that affects the subsidence parameter value of the estimated point. Therefore, next, the subsidence parameter estimator calculates the value of the subsidence parameter of the estimated point based on the subsidence parameter of each known point existing around the estimated point.
このようにして、本発明では、実際に沈下の発生した建設地点における学習用地盤情報と、実際に沈下した既設建物の地盤沈下量に基づいて設定された学習用沈下パラメータの値と、を教師データとして学習するので、第1学習部は、推定地点を含む任意の地点での沈下パラメータの値の算出を精度良く学習することができる。そして、推定地点では、その周辺の既知地点の沈下パラメータの値を用いて、推定地点の沈下パラメータの値を算出するので、たとえば推定地点の地盤情報が無くても、推定地点の位置情報を入力するだけで、経験が豊富でない設計者であっても、合理的かつ客観的に、その推定地点の沈下のし易さの程度を精度良く予測することができる。これにより、地盤沈下による建物不具合を未然に防ぐことができ、建設費用を削減することができる。 In this way, in the present invention, the learning ground information at the construction site where subsidence actually occurred and the value of the learning subsidence parameter set based on the ground subsidence amount of the existing building that actually subsided are taught. Since learning is performed as data, the first learning unit can accurately learn the calculation of the subsidence parameter value at any point including the estimated point. Then, at the estimated point, the subsidence parameter value of the estimated point is calculated using the subsidence parameter values of the known points around it, so even if there is no ground information of the estimated point, for example, the position information of the estimated point is input. Only by doing so, even a designer with little experience can rationally and objectively predict the degree of easiness of subsidence at the estimated point with good accuracy. As a result, it is possible to prevent building failures due to land subsidence and reduce construction costs.
より好ましい態様としては、前記周辺情報抽出部は、前記推定地点から所定の距離内に存在する前記既知地点ごとの前記既知地点情報を抽出するものであり、前記地盤沈下予測システムは、前記周辺情報抽出部で抽出した前記既知地点の数と、前記沈下パラメータ推定部で算出した沈下パラメータの値とから、前記推定地点のボーリング試験および標準貫入試験の実施の要否を判定する試験実施判定部をさらに備える。 As a more preferable aspect, the surrounding information extraction unit extracts the known point information for each of the known points existing within a predetermined distance from the estimated point, and the ground subsidence prediction system extracts the surrounding information a test execution determination unit that determines whether or not to conduct a boring test and a standard penetration test at the estimated points based on the number of known points extracted by the extraction unit and the subsidence parameter values calculated by the subsidence parameter estimation unit; Prepare more.
ここで、推定地点から所定の距離内に存在する既知地点ごとの既知地点情報を抽出した際に、たとえば、抽出した既知地点の個数が、所定の個数よりも少ない場合、算出した推定地点の沈下パラメータの値が所定の値よりも大きい場合、または、これらの両方が成立する場合には、推定地点で地盤沈下する可能性が高い。したがって、この態様によれば、試験実施判定部で、このような場合を考慮して、推定地点のボーリング試験および標準貫入試験の実施の要否を判定することができるので、推定地点での地盤沈下の可能性を、より的確に予測することができる。 Here, when extracting known point information for each known point existing within a predetermined distance from the estimated point, for example, if the number of extracted known points is less than a predetermined number, the calculated subsidence of the estimated point If the value of the parameter is greater than the predetermined value, or if both of them hold, there is a high possibility that the ground will subside at the estimated point. Therefore, according to this aspect, in the test implementation determination unit, considering such a case, it is possible to determine whether or not to implement the boring test and the standard penetration test at the estimated point, so the ground at the estimated point The possibility of subsidence can be predicted more accurately.
さらに、好ましい態様としては、前記地盤沈下予測システムは、複数の既設建物の建設地点から選定した選定地点の学習用沈下パラメータの値と、前記選定地点から所定の範囲内の複数の周辺の建設地点の学習用沈下パラメータの値と、前記選定地点に対する前記周辺の建設地点ごとの相対的な位置情報を特定した学習用位置情報と、を1つのデータ群とし、異なる前記選定地点ごとに得られる複数のデータ群を教師データとし、任意の地点の周辺に存在する複数の周辺地点の沈下パラメータの値と、前記任意の地点に対する前記周辺地点ごとの相対的な位置情報とから、前記任意の地点の沈下パラメータの値の算出を機械学習により学習した第2学習部をさらに備えており、前記沈下パラメータ推定部は、前記周辺沈下パラメータ算出部で算出した前記既知地点ごとの前記沈下パラメータの値と、前記推定地点に対する前記既知地点ごとの相対的な位置情報とを、前記第2学習部に入力し、前記第2学習部に前記推定地点の沈下パラメータの値を算出させる。 Furthermore, as a preferred embodiment, the ground subsidence prediction system comprises a learning subsidence parameter value of a selected point selected from a plurality of existing building construction points, and a plurality of surrounding construction points within a predetermined range from the selected point. and the learning positional information specifying the relative positional information of each of the surrounding construction sites with respect to the selected site, as one data group, and a plurality of data obtained for each of the different selected sites Using the data group of as teacher data, the value of the subsidence parameter of a plurality of peripheral points existing around an arbitrary point, and the relative position information of each peripheral point with respect to the arbitrary point, the arbitrary point The subsidence parameter estimation unit calculates the subsidence parameter value for each of the known points calculated by the surrounding subsidence parameter calculation unit, and Relative position information of each of the known points with respect to the estimated point is input to the second learning unit, and the second learning unit is caused to calculate a subsidence parameter value of the estimated point.
この態様によれば、第2学習部は、所定の建設地点の沈下パラメータの値の算出を、その周辺の建設地点の相対的な位置情報とその沈下パラメータの値から学習したものである。したがって、沈下パラメータ推定部により、推定地点に対する周辺の既知地点の方位を加味した、推定地点の沈下パラメータの値を、第2学習部で算出させることができる。これにより、より正確な沈下パラメータの値を算出することができる。 According to this aspect, the second learning unit learns the calculation of the subsidence parameter value of the predetermined construction site from the relative positional information of the surrounding construction sites and the subsidence parameter value thereof. Therefore, the subsidence parameter estimating section can cause the second learning section to calculate the value of the subsidence parameter of the estimated point, taking into account the directions of known points around the estimated point. This makes it possible to calculate a more accurate subsidence parameter value.
より好ましい態様としては、前記第2学習部は、前記地域情報格納部で格納された前記所定の地域内の既知地点ごとの位置情報と、前記第1学習部に算出された前記所定の地域内の既知地点ごとの沈下パラメータの値と、を教師データとして機械学習したものである。この態様によれば、第2学習部は、推定地点を含む所定の地域において、沈下パラメータの値の算出を学習するので、推定地点の沈下パラメータの値の精度を高めることができる。 As a more preferable aspect, the second learning unit stores location information for each known point within the predetermined area stored in the area information storage unit, and the location information within the predetermined area calculated by the first learning unit. The value of the subsidence parameter for each known point of , was machine-learned as teacher data. According to this aspect, since the second learning unit learns the calculation of the subsidence parameter value in the predetermined area including the estimated point, it is possible to improve the accuracy of the subsidence parameter value of the estimated point.
より好ましい態様としては、前記沈下パラメータ推定部は、前記周辺沈下パラメータ算出部で算出した前記既知地点ごとの前記沈下パラメータの値と、前記推定地点から抽出した各既知地点までの距離に基づいて、前記推定地点の沈下パラメータの値を算出する。 In a more preferred embodiment, the subsidence parameter estimating unit, based on the value of the subsidence parameter for each of the known points calculated by the surrounding subsidence parameter calculating unit, and the distance from the estimated point to each known point extracted, A subsidence parameter value for the estimated point is calculated.
この態様によれば、推定地点の沈下パラメータの値は、既知地点ごとの前記沈下パラメータの値と、推定地点から抽出した各既知地点までの距離が大きく影響するので、これらの値を用いることで、推定地点の沈下パラメータの値の精度を高めることができる。 According to this aspect, the value of the subsidence parameter of the estimated point is greatly affected by the value of the subsidence parameter for each known point and the distance from the estimated point to each known point extracted. , the accuracy of the value of the subsidence parameter at the estimated point can be increased.
本発明によれば、特定の地点を地盤沈下の推定地点として、推定地点の地盤沈下を正確に予測することができる。 According to the present invention, it is possible to accurately predict ground subsidence at a specific point as an estimated ground subsidence point.
以下、本発明の第1および第2実施形態に係る地盤沈下予測システムを、図面を参照しながら説明する。 Hereinafter, ground subsidence prediction systems according to first and second embodiments of the present invention will be described with reference to the drawings.
〔第1実施形態〕
1.地盤沈下予測システム1の装置構成について
本実施形態に係る地盤沈下予測システム1は、たとえば、メインサーバ(ホストコンピュータ)10と、携帯端末(通信端末)20を用いて、特定の地点を地盤沈下の推定地点として、推定地点の地盤沈下を予測するものである。以下に、地盤沈下予測システム1を説明する。
[First embodiment]
1. About the device configuration of the ground
地盤沈下予測システム1のメインサーバ10は、地盤沈下の状態(具体的には沈下パラメータの値)を算出する(学習する)人工知能を有した演算装置10Aと、演算装置10Aにデータを入力するキーボードなどの入力装置10Bと、を備えている。演算装置10Aは、後述する所定の地域の地域情報などのデータ、人工知能により機械学習されたプログラムが記憶されたROMまたはRAM等からなる記憶部10aと、機械学習されたプログラムにより沈下パラメータの値を算出するCPU等からなる演算部10bと、を備えている。具体的には、記憶部10aには、図2に示す、沈下パラメータ学習部11のプログラム、地域情報格納部12で格納された地域情報などが記憶されており、演算部10bでは、図2に示す沈下パラメータ学習部11による学習(演算)が実行される。
The
地盤沈下予測システム1の携帯端末20は、ネットワークを介してメインサーバ10と通信可能であり、推定地点の沈下パラメータの値を算出する演算装置20Aと、データの入力および表示を行うタッチパネルなどの表示・入力部20Bと、を備えている。演算装置10Aは、メインサーバ10からのデータおよび表示・入力部20Bで入力されたデータを記憶するROMまたはRAM等からなる記憶部20aと、推定地点の沈下パラメータを算出するCPU等からなる演算部20bと、を備えている。具体的には、記憶部20aには、メインサーバ10からのデータおよび表示・入力部20Bで入力されたデータ、演算部20bで演算するためのアプリケーション(プログラム)などが記憶され、演算部20bでは、図2に示す周辺情報抽出部22、周辺沈下パラメータ算出部23、および沈下パラメータ推定部24などによる演算が実行される。
The
なお、本実施形態では、地盤沈下予測システム1は、メインサーバ10と、携帯端末20とを個別に備えており、携帯端末20は、上述するアプリケーションがインストールされていれば、特にその個数は限定されるものではない。
In this embodiment, the ground
また、地盤沈下予測システム1が、メインサーバ10と携帯端末20とで構成されるのではなく、これらが1つのシステムとして構成されていてもよい。具体的には、メインサーバ10の記憶部10aで、携帯端末20の記憶部20aで記憶されたデータを記憶し、メインサーバ10の演算部10bで、携帯端末20の演算部20bによる演算を行うことで、携帯端末20を省略してもよい。また、携帯端末20では、上述した記憶および演算を行わず、メインサーバ10ですべての記憶および演算を行って、携帯端末20は、推定地点の位置情報の入力と、メインサーバ10により演算された推定地点の沈下パラメータの値の出力とを実行するだけの役割を果たしてもよい。以下に、本実施形態では、地盤沈下予測システム1のソフトウエアの機能を、図2のブロック図を用いて説明する。
Further, the ground
本実施形態では、メインサーバ10の演算装置10Aは、沈下パラメータ学習部(第1学習部)11と、地域情報格納部12とを備えている。携帯端末20の演算装置20Aは、推定地点設定部21と、周辺情報抽出部22と、周辺沈下パラメータ算出部23と、沈下パラメータ推定部24と、を備えている。メインサーバ10の演算装置10Aのデータと、携帯端末20の演算装置20Aのデータとは、ネットワークを介して相互に送受信可能になっている。
In this embodiment, the
2.地盤沈下予測システム1のモデルについて
2-1.沈下パラメータ学習部11
沈下パラメータ学習部11(以下「第1学習部11」という)は、複数の既設建物A、B、C、…の建設地点P1、P2、P3、…における学習用地盤情報と、各建設地点P1、P2、P3、…に紐付けられた(関連付けられた)学習用沈下パラメータの値と、を教師データとして、任意の地点の地盤情報から、その任意の地点の沈下パラメータの値の算出を機械学習による学習したものである。ここで、第1学習部11により学習される教師データの対象となる「地点」は、既設建物が建設された建設地内に存在する建設地点であり、実際に地盤沈下し、地盤沈下量が測定された建設地内の地点、または、既設建物の仕様等から地盤沈下の沈下量が推定された建設地内の地点である。学習後の第1学習部11に入力されるデータの対象となる「地点」は、所定の地域においてボーリング試験および標準貫入試験が実施された既知地点であり、これらの試験結果から入力されたデータにより出力されるデータは、既知地点の沈下パラメータの値である。
2. About the model of the ground
The subsidence parameter learning unit 11 (hereinafter referred to as “
具体的には、「学習用地盤情報」は、既設建物A、B、C、…の建設地点P1、P2、P3、…のボーリング試験により得られた、建設地点P1、P2、P3、…の地表からの深度、深度に応じた土質を数値化して設定された土質パラメータの値、および深度に応じた標準貫入試験におけるN値、を含むものである。ここで、入力される「深度」は、たとえば、各土質を構成する層の最下部の深さなど、土質を構成する深度と、その土質の層厚みの情報が特定できるような「深度」を入力することができるのであれば、土質パラメータの値とN値に関連付けられた「深度」の値は特に限定されるものではない。 Specifically, the "ground information for learning" is obtained by boring tests at the construction points P1, P2, P3, . . . of the existing buildings A, B, C, . It includes the depth from the ground surface, the value of the soil parameter set by quantifying the soil quality according to the depth, and the N value in the standard penetration test according to the depth. Here, the "depth" to be input is, for example, the depth that constitutes the soil, such as the depth of the bottom of the layer that constitutes each soil, and the "depth" that can specify the information of the layer thickness of the soil. The "depth" value associated with the soil parameter value and the N value is not particularly limited as long as it can be input.
一般的に、ボーリング試験において、地表からの深度に応じて、土質の種類が調査される。ここで、「土質パラメータ」とは、土質ごとに数値化されて設定されたパラメータである。たとえば、土質には、粘土、シルト、細砂、粗砂、細礫、中礫、粗礫、コブル、またはボルタなどを挙げることができ、この土質ごとに、数値化されている。たとえば、土質の種類に応じた粒度分布に対応する数値が、設定されていてもよい。たとえば、土質の種類に応じた粒度分布が大きくなるに従って、建物からの荷重を支持する支持力が大きい(地盤沈下し難い)ことから、この順に、小さい数値となるように数値が割り当てられている。なお、たとえば、図3に示すように、土質A~土質Eに対して、土質パラメータの数値が割り当てられる。 Generally, in the boring test, the type of soil is investigated according to the depth from the ground surface. Here, the "soil parameter" is a parameter quantified and set for each soil property. For example, the soil quality can include clay, silt, fine sand, coarse sand, fine gravel, gravel, coarse gravel, cobble, or volta, and is quantified for each soil quality. For example, a numerical value corresponding to the particle size distribution according to the type of soil may be set. For example, as the particle size distribution according to the type of soil increases, the bearing capacity to support the load from the building increases (ground subsidence is less likely). . For example, as shown in FIG. 3, soil parameters A to E are assigned numeric values of soil parameters.
「標準貫入試験におけるN値」とは、標準貫入試験(JIS A 1219に準拠した試験)で測定されたN値、スウェーデン式サウンディング試験(JIS A 1212に準拠した試験)に基づく換算N値などを挙げることができ、N値が大きいほど、建物からの荷重を支持する支持力が大きい。 "N value in the standard penetration test" refers to the N value measured in the standard penetration test (test in accordance with JIS A 1219), the converted N value based on the Swedish sounding test (test in accordance with JIS A 1212), etc. The larger the N value, the greater the bearing capacity to support the load from the building.
図3に示すように、各建設地点P1、P2、P3、…において、その深度、その深度ごとの土質パラメータの値、およびその深度ごとのN値が教師データとして割り当てられる。なお、これらの学習用地盤情報のうち、各建設地点P1、P2、P3、…の土質パラメータの値を、たとえば、ボーリング柱状図の画像から土質を認識し、機械学習を利用した人工知能を用いて設定してもよい。 As shown in FIG. 3, at each construction point P1, P2, P3, . Among these ground information for learning, the values of the soil parameters of each construction site P1, P2, P3, . can be set.
各建設地点に関連付けられた「学習用沈下パラメータ」は、既設建物の地盤沈下量に基づいて、地盤沈下のし易さに応じた値に設定される変数であり、この変数の値が大きいほど地盤沈下がし易いような値に設定される。 The "learning subsidence parameter" associated with each construction site is a variable that is set to a value according to the likelihood of ground subsidence based on the amount of ground subsidence of the existing building. It is set to a value that facilitates land subsidence.
本実施形態では、「学習用沈下パラメータ」の値は、建設地点において、地盤沈下量に対して、既設建物の仕様(たとえば、建物の用途、回数、建物面積、建物の基礎の仕様を加味した仕様)を数値化した値で正規化された値である。たとえば、この沈下パラメータは、同じ仕様の建物に対して、地盤沈下のし易さ(地盤沈下のリスク)を数値化したものである。換言すると、沈下パラメータの大きさが大きくなるに従って、地盤沈下がし易くなるため、地盤沈下し難い基礎構造が選定される。このような学習用沈下パラメータの値は、学習の基礎となる教師データであることから、既設建物の沈下量と既設建物の仕様とから、熟練者がその経験に基づいて数値を設定してもよく、以下のように設定してもよい。 In this embodiment, the value of the "learning subsidence parameter" is the amount of ground subsidence at the construction site, considering the specifications of the existing building (for example, the use of the building, the number of times, the building area, specification) is normalized by a numerical value. For example, this subsidence parameter quantifies the easiness of ground subsidence (risk of ground subsidence) for buildings with the same specifications. In other words, as the magnitude of the subsidence parameter increases, the ground subsidence becomes more likely, so a foundation structure that is less likely to subside is selected. Since these learning subsidence parameter values are teaching data that serve as the basis for learning, the amount of subsidence of the existing building and the specifications of the existing building can be used by experts to set the values based on their experience. Well, you can set it like this:
具体的には、「学習用沈下パラメータ」の値は、既設建物の地盤沈下量を、以下に示す既設建物の建物仕様に対応した値で除算した値であってもよい。本実施形態では、「既設建物の建物仕様に対応した値」は、たとえば、(1)既設建物の用途を数値化して設定された用途パラメータの値、(2)既設建物の階数、(3)既設建物の建物面積、および(4)既設建物に選択された基礎形式を数値化して設定された基礎パラメータの値、を乗算した値である。 Specifically, the value of the "learning subsidence parameter" may be a value obtained by dividing the ground subsidence amount of the existing building by a value corresponding to the building specifications of the existing building shown below. In the present embodiment, the "value corresponding to the building specifications of the existing building" is, for example, (1) a usage parameter value set by quantifying the usage of the existing building, (2) the number of floors of the existing building, and (3) It is a value obtained by multiplying the building area of the existing building and (4) the value of the foundation parameter set by digitizing the foundation type selected for the existing building.
「用途パラメータ」は、既設建物の用途に応じたパラメータであり、その用途が、たとえば、病院、学校、工場、またはオフィスビル等に合わせて、既設建物の重量が変動することから、これらの用途に応じて数値化され、設定されている。たとえば、用途パラメータの値は、地盤沈下し易い用途の順に、大きい数値となるように数値が割り当てられている。 "Usage parameter" is a parameter according to the usage of the existing building. is digitized and set according to For example, the application parameter values are assigned numerical values in descending order of the degree of land subsidence.
「建物階数」は、既設建物または未設建物の階数であり、「建物面積」は、既設建物または未設建物の基礎部分に相当する面積である。 The "building floor number" is the number of floors of an existing building or a non-constructed building, and the "building area" is an area corresponding to the base portion of the existing building or non-constructing building.
「基礎パラメータ」は、既設建物または未設建物に選択された基礎形式を数値化して設定されたパラメータである。具体的には、(1)杭を設けた直接基礎構造、(2)杭を設けたパイル・ドラフト基礎構造、(3)摩擦杭による杭基礎構造、(4)支持層により支持される支持杭による杭基礎構造などを挙げることができ、これらの基礎形式(基礎構造の形式)を数値化して設定される。たとえば、地盤沈下し易い基礎構造の順に、大きい数値となるように数値が割り当てられている。 "Foundation parameters" are parameters set by digitizing the foundation type selected for an existing building or an unconstructed building. Specifically, (1) spread foundation structure with piles, (2) pile draft foundation structure with piles, (3) pile foundation structure with friction piles, and (4) support piles supported by support layers. The pile foundation structure etc. can be mentioned, and these foundation types (form of foundation structure) are set numerically. For example, numerical values are assigned in descending order of the subsidence of the foundation structure.
さらに、杭基礎構造の場合には、摩擦杭の長さ、摩擦杭の種類等により、さらに、地盤沈下し易い順に、これに応じた大きさの数値が設定されてもよい。さらに、基礎構造の他にも、例えば、セメントミルクにより地盤改良された基礎構造等を含んでもよく、この構造にも、上述した如き数値が設定される。 Furthermore, in the case of a pile foundation structure, depending on the length of the friction pile, the type of the friction pile, etc., the numerical value of the size may be set according to the order of susceptibility to ground subsidence. Furthermore, in addition to the foundation structure, for example, a foundation structure improved with cement milk may also be included, and the above-mentioned numerical values are set for this structure as well.
このようにして、教師データとして、図3に示すように、建設地点ごとの「学習用沈下パラメータ」の値が、その建設地点ごとの学習用地盤情報に紐付けて設定される。たとえば、建設地点P1では、学習用沈下パラメータの値が50であり、建設地点P2では、学習用沈下パラメータの値が32である。したがって、建設地点P1と建設地点P2に同じ仕様の建物を建設した場合には、建設地点P1の方が、地盤沈下がし易いことになる。図3に示す建設地点P1~P4では、建設地点P4、P1、P3、P2の順に、沈下し易い地盤を有することになる。 In this way, as shown in FIG. 3, the value of the "learning subsidence parameter" for each construction site is set as training data in association with the learning ground information for each construction site. For example, the value of the learning subsidence parameter is 50 at the construction point P1, and the value of 32 is the learning subsidence parameter at the construction point P2. Therefore, if buildings of the same specifications are constructed at the construction point P1 and the construction point P2, the construction point P1 is more likely to undergo land subsidence. The construction points P1 to P4 shown in FIG. 3 have ground that is likely to subside in the order of construction points P4, P1, P3, and P2.
図3に示す、各建設地点P1、P2、…ごとにおける学習用地盤情報と、学習用沈下パラメータの値を教師データとして、任意の地点の地盤情報から任意の地点の沈下パラメータの値の算出を機械学習により学習する。 Using the ground information for learning and the value of the settlement parameter for learning for each construction site P1, P2, . . . shown in FIG. Learn by machine learning.
本実施形態では、第1学習部11は、これまでに建設された既設建物A、B、C、…の入力値に対して算出される出力値が、各既設建物A、B、C、…で、実際の地盤沈下量に基づいて地盤沈下のし易さに応じて設定された沈下パラメータ(実際の沈下パラメータ)の値に収束するように、沈下パラメータの値の算出を機械学習により学習する。
In this embodiment, the
この学習は、例えば、学習用地盤情報の入力値の個数に応じた変数からなる所定の数式に対して、各変数に乗じられる補正係数を、繰り返し補正することにより行ってもよい。本実施形態では、図4に示すように、第1学習部11は、ディープニューラルネットワーク(DNN:以下「ニューラルネットワーク」という)11’を備えており、ニューラルネットワーク11’は、建設地点の地盤情報を入力値とし、沈下パラメータの値を出力値として算出するものである。
This learning may be performed, for example, by repeatedly correcting a correction coefficient by which each variable is multiplied to a predetermined mathematical formula consisting of variables corresponding to the number of input values of the ground information for learning. In this embodiment, as shown in FIG. 4, the
本実施形態では、ニューラルネットワーク11’は、入力層11Aを有している。入力層11Aは、各建設地点P1、P2、P3、P4、…に対して、地表からの深度N1、N2、…と、これらに応じた土質パラメータの値およびN値を入力値とする。入力層11Aは、この入力値の個数に合わせた複数の入力層ニューロン素子11aで構成される。入力層11Aの入力層ニューロン素子11aは、入力する条件データの個数に応じて、増減することができ、この増減による入力されるデータの個数に応じたニューラルネットワーク11’のモデルが構築される。
In this embodiment, the neural network 11' has an
ニューラルネットワーク11’は、出力層11Eを有している。出力層11Eは、沈下パラメータの値を出力する1つの出力層ニューロン素子11eで構成される。ニューラルネットワーク11’は、3つの中間層11B、11C、11Dを有している。3つの中間層11B、11C、11Dは、入力層11Aと出力層11Eとの間に設けられている。各中間層11B、11C、11Dは、これらの素子に直接的または間接的に結合された複数の中間層ニューロン素子11b、11c、11dを含む。
The neural network 11' has an
なお、本実施形態では、中間層が3つの層で構成されるが、たとえば、中間層が、1つ、2つの層、4つ以上の層で構成されていてもよく、中間層が3つの層に限定されるものではない。さらに、各中間層11B、11C、11Dを構成する中間層ニューロン素子11b、11c、11dは、入力層ニューロン素子11aの個数に応じた個数であるが、入力されるデータ数に応じて、入力層ニューロン素子11aの個数を変化させてもよく、入力されるデータ数と異なる個数の中間層のニューロン素子を、各中間層が備えてもよい。
In this embodiment, the intermediate layer is composed of three layers, but for example, the intermediate layer may be composed of one, two, four or more layers, It is not limited to layers. Further, the number of intermediate
中間層11B、11C、11Dは、入力層11A側の同じ層にあるニューロン素子のニューロンパラメータの値を用いて、所定の演算を行い、その演算結果を、出力層11E側のニューロン素子に出力するものである。具体的には、中間層ニューロン素子11b、11c、11dおよび出力層ニューロン素子11eは、入力層ニューロン素子11aまたは入力層11A側の中間層ニューロン素子11b、11c、11dから入力されるニューロンパラメータの値を用いて、所定の演算を行う。
The
具体的には、演算を実行する中間層ニューロン素子11b、11c、11dと、出力層ニューロン素子11eは、それぞれ所定の活性化関数を有しており、入力されたデータ(パラメータの値)をその活性化関数に代入することにより、ニューロンパラメータの値を演算する。たとえば、中間層11Bの各中間層ニューロン素子11bは、各入力層ニューロン素子11aの入力値がニューロンパラメータの値として入力され、活性化関数により、入力層ニューロン素子11aごとのニューロンパラメータの値が演算される。
Specifically, each of the intermediate
この際、各中間層11B、11C、11Dの中間層ニューロン素子11b、11c、11dで出力されるニューロンパラメータの値は、その素子内において、活性化関数により演算されたニューロンパラメータの値に対して、重み付け係数が乗算されることで算出され、出力層ニューロン素子11eまたは出力層11E側の中間層ニューロン素子11c、11dに出力される。
At this time, the values of the neuron parameters output from the intermediate
本実施形態では、出力層ニューロン素子11eで出力された既設建物A、B、C、…の沈下パラメータの値が、既設建物A、B、C、…の実際の沈下パラメータの値に対して予め設定した範囲内に収束するまで、各中間層ニューロン素子11b、11c、11dのニューロンパラメータの値に乗算される重み付け係数を繰り返し補正する。このようにして、本実施形態では、機械学習により、沈下パラメータの値の算出(方法)を予め学習することができる。このようにして、複数の建設地点P1、P2、P3、…における学習用地盤情報と学習用沈下パラメータの値とを教師データとして、第1学習部11の機械学習が実行され、任意の地点の地盤情報から、沈下パラメータの値を算出するモデルが生成される。
In this embodiment, the subsidence parameter values of the existing buildings A, B, C, . The weighting coefficients multiplied by the neuron parameter values of the intermediate
なお、本実施形態では、中間層ニューロン素子11b、11c、11dにおいて、活性化関数により算出されたニューロンパラメータの値に対して乗算される重み付け係数を、各中間層ニューロン素子11b、11c、11dに設けたが、これに加えて、出力層ニューロン素子11eにもさらに活性化関数により算出されたニューロンパラメータの値に対して乗算される重み付け係数を設けてもよい。
In the present embodiment, the weighting coefficients multiplied by the neuron parameter values calculated by the activation functions in the intermediate
2-2.地域情報格納部12
地域情報格納部12には、図5に示すように、所定の地域の地図データと、地図データの地図M1上の複数の既知地点T1、T2、…における既知地点情報が格納されている。この既知地点情報は、既知地点T1、T2、…ごとの位置情報、既知地点T1、T2、…ごとのボーリング試験により得られた、地表からの深度、深度に応じた土質を数値化して設定された土質パラメータの値、および深度に応じた標準貫入試験におけるN値、を含むものであり、これらのデータは既知地点ごとに紐付けられている。ここで示す、深度、土質パラメータの値、およびN値の地盤情報は、学習用地盤情報で説明したものと同じ種類の情報であるので、詳細な説明は省略する。なお、地図データの地図M1上の既知地点に、機械学習で用いた建設地点が含まれていてもよい。
2-2. Regional
As shown in FIG. 5, the area
具体的には、各既知地点は、ボーリング試験および標準貫入試験が実施された地点であり、これらの試験により得られた地盤情報が、地図M1上の各既知地点T1、T2、…の位置情報に紐付けられている。なお、この既知地点では、現在または過去に建物が建設されていなくてもよく、仮に、この既知地点において建物による地盤沈下が発生していた場合には、地域情報格納部12に、実際に測定されたまたは推定された地盤沈下量に関連付けられた沈下パラメータをさらに格納してもよい。なお、地図データの対象となる所定の地域とは、たとえば、区または市などの地域であってもよく、1つまたは複数の共通の土質の層が連続して含む領域であってもよい。このような地図データは、メインサーバ10から携帯端末20に送信され、携帯端末20の表示・入力部20Bに表示される。
Specifically, each known point is a point where a boring test and a standard penetration test were performed, and the ground information obtained by these tests is the position information of each known point T1, T2, ... on the map M1. tied to. It should be noted that, at this known point, it is not necessary that a building has been constructed at present or in the past. A subsidence parameter associated with the estimated or estimated ground subsidence amount may also be stored. Note that the predetermined area that is the target of the map data may be, for example, an area such as a ward or city, or an area that includes one or a plurality of common soil layers continuously. Such map data is transmitted from the
2-3.推定地点設定部21
図5に示すように、推定地点設定部21は、地図M1上の推定地点Gを設定する。本実施形態では、携帯端末20の表示・入力部20Bに表示された地図M1上に、操作者が推定地点Gの位置情報を入力する。これにより、携帯端末20の演算装置20Aにおいて、地図M1上の推定地点Gが設定される。
2-3. Estimated
As shown in FIG. 5, the estimated
2-4.周辺情報抽出部22
図2および図5に示すように、周辺情報抽出部22は、推定地点設定部21で設定された推定地点Gの周辺に存在する既知地点T1、T7、T8、…ごとの既知地点情報を抽出する(図6参照)。本実施形態では、図5、図7に示すように、推定地点Gを含む所定の範囲RG内に存在する既知地点T1、T7、T8、…ごとの既知地点情報を抽出する。
2-4. Peripheral
As shown in FIGS. 2 and 5, the surrounding
ここで、本実施形態では、周辺情報抽出部22は、推定地点Gから所定の距離r内に存在する既知地点ごとの既知地点情報を抽出する。具体的には、図5および図7に示すように、周辺情報抽出部22が抽出する領域は、推定地点Gから一定の距離r(たとえば3km)内にある円領域(半径3kmの円内)である。しかしながら、周辺情報抽出部22は、たとえば、推定地点Gからの距離が短い順に、特定(一定)の個数の既知地点の既知地点情報を抽出してもよい。
Here, in the present embodiment, the surrounding
2-5.周辺沈下パラメータ算出部23
周辺沈下パラメータ算出部23は、図6に示す周辺情報抽出部22で抽出した既知地点T1、T7、T8、…ごとの既知地点情報を、第1学習部11に入力し、図6に示す既知地点T1、T7、ごとの沈下パラメータの値を第1学習部11に算出させるものである。第1学習部11は、上述した如く、任意の地点の地盤情報から、その任意の地点の沈下パラメータの値の算出が学習されているので、既知地点T1、T7、T8、…ごとの沈下パラメータの値を精度良く算出することができる。
2-5. Surrounding
The surrounding
なお、ここで、地域情報格納部12において、抽出された既知地点に、実際の沈下パラメータの値が既に格納されている場合には、周辺沈下パラメータ算出部23は、この既知地点における沈下パラメータの値の算出を行わなくてもよい。また、別の態様としては、このように、演算によらない実際の沈下パラメータの値が既に格納されている場合であっても、周辺沈下パラメータ算出部23は、この既知地点における沈下パラメータの値の算出を行い、算出した沈下パラメータの値と実際の沈下パラメータの値との差分または比率から、他の既知地点において算出した沈下パラメータの値を補正してもよい。
Here, if the area
2-6.沈下パラメータ推定部24
沈下パラメータ推定部24は、図7に示すように、周辺沈下パラメータ算出部23で算出した既知地点T1、T7、T8、…ごとの沈下パラメータに基づいて、推定地点Gの沈下パラメータの値を算出する。ここで、推定地点Gの沈下パラメータの値を、たとえば、算出した既知地点T1、T7、T8、…ごとの沈下パラメータの平均値により算出してもよい。
2-6. Subsidence
As shown in FIG. 7, the
しかしながら、本実施形態では、周辺沈下パラメータ算出部23で算出した既知地点T1、T7、T8、…ごとの沈下パラメータの値と、推定地点Gからの抽出した既知地点ごとの距離L1、L7、L8…に基づいて、推定地点Gの沈下パラメータの値を算出する。これにより、推定地点Gの沈下パラメータの値は、既知地点T1、T7、T8、…ごとの沈下パラメータの値と、推定地点Gからの抽出した既知地点T1、T7、T8、…ごとの距離L1、L7、L8…が大きく影響するので、これらの値を用いることで、推定地点Gの沈下パラメータの値の精度を高めることができる。 However, in the present embodiment, the subsidence parameter value for each of the known points T1, T7, T8, . . . . . , the value of the subsidence parameter of the estimated point G is calculated. As a result, the subsidence parameter values of the estimated point G are the subsidence parameter values of the known points T1, T7, T8, . , L7, L8, .
より具体的には、本実施形態では、沈下パラメータ推定部24は、推定地点からの抽出した既知地点ごとの距離が大きくなるに従って、推定地点Gの沈下パラメータの値が小さくなるように、沈下パラメータの値を補正した後、これらの値の平均値を、推定地点Gの沈下パラメータの値と推定する。これにより、推定地点Gからより近い既知地点の沈下パラメータの値を、推定地点Gの沈下パラメータの値に反映することができる。このような算出を以下に示す数1に基づいて行ってもよい。
More specifically, in this embodiment, the subsidence
ここで、DGは、推定地点の沈下パラメータの値であり、Dkは、抽出した各既知地点の沈下パラメータの値であり、Lkは、推定地点からの抽出した既知地点ごとの距離であり、rは、推定地点からの既知地点を抽出する距離であり、Nは、抽出した既知地点の個数である。なお、αは、補正係数であり、この補正係数は、たとえば、数1で算出した推定地点の沈下パラメータの値が、推定地点における実際の沈下パラメータの値に、一致するように設定される定数である。但し、この式には、周辺情報抽出部22で抽出されていない既知地点のDkは含まれず、抽出されていない既知地点は、Nの個数にはカウントされない。
従って、この数1を、図7に示す場合に適用すると、以下の数2のようになる。
Here, DG is the value of the subsidence parameter of the estimated point, Dk is the value of the subsidence parameter of each extracted known point, Lk is the distance of each extracted known point from the estimated point, r is the distance for extracting known points from the estimated point, and N is the number of extracted known points. Note that α is a correction coefficient, and this correction coefficient is a constant set so that, for example, the value of the subsidence parameter at the estimated point calculated by
Therefore, applying this
なお、ここで、地域情報格納部12において、抽出された既知地点のいくつかに、実際の沈下パラメータの値が既に格納されている場合には、沈下パラメータ推定部24は、推定地点の沈下パラメータの値の算出に、実際の沈下パラメータを用いてもよい。
If actual subsidence parameter values are already stored for some of the extracted known points in the area
3.地盤沈下予測システム1を用いた地盤沈下予測方法について
以下に、図8を参照して、本実施形態に係る地盤沈下予測システム1を用いた地盤沈下予測方法を説明する。まず、ステップS1では、第1学習部11により、学習用地盤情報と学習用沈下パラメータの値を用いて機械学習を行う。これにより、地盤情報から沈下パラメータの値を出力するニューラルネットワーク11’のモデルが構築される。
3. Land Subsidence Prediction Method Using Land Subsidence Prediction System 1 A land subsidence prediction method using the land
次に、ステップS2で、地域情報格納部12に所定の地域の地図データと、地図データの地図上の複数の既知地点T1、T2、…における既知地点情報とを格納する。なお、ステップS1とステップS2の順序は問わず、メインサーバ10で実行される。次に、ステップS3で、携帯端末20のアプリケーションを起動し、メインサーバ10に格納された地図データの地図Mを、既知地点T1、T2、…とともに表示する。
Next, in step S2, map data of a predetermined area and known point information of a plurality of known points T1, T2, . . . Note that the steps S1 and S2 are executed by the
次に、ステップS4では、操作者が、携帯端末20に表示された地図M上に推定地点Gの位置情報を入力する。これにより、推定地点設定部21が、地図データに推定地点Gの位置情報を紐付けることで、地図M上の推定地点Gを設定する。次に、ステップS5では、周辺情報抽出部22が、推定地点Gの周辺に存在する既知地点T1、T7、T8、…ごとの既知地点情報(周辺情報)を抽出する。
Next, in step S<b>4 , the operator inputs the positional information of the estimated point G on the map M displayed on the
次に、ステップS6では、周辺沈下パラメータ算出部23は、携帯端末20からネットワークを介して各既知地点T1、T7、T8、…ごとの地盤情報をメインサーバ10の演算装置10Aに送信し、メインサーバ10内の第1学習部11で各既知地点T1、T7、T8、…ごとの沈下パラメータの値D1、D7、D8、…を算出する。第1学習部11は、算出した結果を、ネットワークを介して携帯端末20の演算装置20Aに送信し、この結果が沈下パラメータ推定部24に入力される。
Next, in step S6, the peripheral
次に、ステップS7では、沈下パラメータ推定部24が、推定地点Gからの抽出した既知地点ごとの距離L1、L7、L8…を算出する。ステップS8では、沈下パラメータ推定部24が、既知地点T1、T7、T8、…ごとの沈下パラメータの値D1、D7、D8、…と、推定地点Gから抽出した既知地点までの距離L1、L7、L8…に基づいて、推定地点Gの沈下パラメータの値を算出する。最後に、ステップS9で、推定地点Gの沈下パラメータの値を、携帯端末20の表示・入力部20Bに入力し、この結果が表示される。
Next, in step S7, the
本実施形態では、実際に沈下の発生した建設地点における学習用地盤情報と、実際に沈下した既設建物の沈下量に基づいて設定された学習用沈下パラメータの値と、を教師データとして学習するので、第1学習部11は、任意の地点での沈下パラメータの値の算出を精度良く学習することができる。推定地点Gでは、その周辺の既知地点の沈下パラメータの値を用いて、推定地点Gの沈下パラメータの値を算出するので、推定地点Gの地盤情報が無くても、推定地点Gの沈下のし易さの程度を精度良く予測することができる。
In this embodiment, the learning ground information at the construction site where subsidence actually occurred and the value of the learning subsidence parameter set based on the actual subsidence amount of the existing building are learned as teacher data. , the
〔第2実施形態〕
以下に、本発明の第2実施形態に係る地盤沈下予測システム1を、図9~図13を参照しながら説明する。第2実施形態に係る地盤沈下予測システム1が、第1実施形態のものと相違する点は、(1)図9に示すように、沈下推定学習部13を新たに設け、これを用いて沈下パラメータ推定部24で、沈下パラメータを算出した点と、(2)試験実施判定部25を新たに設けた点と、(3)第1学習部の学習に所定の地域の既知地点情報を教師データとして用いて補正した(再学習させた)点とである。したがって、第1実施形態と同じ構成については、その詳細な説明を省略し、相違する構成を以下に説明する。
[Second embodiment]
A ground
沈下推定学習部13(以下「第2学習部13」という)は、任意の地点の周辺に存在する複数の周辺地点の沈下パラメータの値と、任意の地点に対する前記周辺地点ごとの相対的な位置情報とから、任意の地点の沈下パラメータの値の算出を機械学習により学習したものである。
The subsidence estimation learning unit 13 (hereinafter referred to as the “
第2学習部13の教師データは、たとえば、図10に示すように、複数の既設建物の建設地点W1~W37から選定した選定地点の学習用沈下パラメータの値と、選定地点から所定の範囲内の複数の周辺の建設地点の学習用沈下パラメータの値と、選定地点に対する周辺の建設地点ごとの相対的な位置情報を特定した学習用位置情報とを、1組の教師データ(1つのデータ群)とし、異なる選定地点ごとに得られる複数のデータ群である。
For example, as shown in FIG. 10, the training data of the
たとえば、図10では、所定の地域の地図M2上の建設地点W1を選定地点に選定する。選定した建設地点W1に対して、建設地点W1から所定の範囲R1(たとえば距離r)内の周辺の建設地点W6、W7、W9、…を抽出する。そして、建設地点W1の沈下パラメータの値と、周辺の建設地点W6、W7、W9、…ごとの沈下パラメータの値とを、学習用沈下パラメータの値として用いる。建設地点(選定地点)W1に対する周辺の建設地点W6、W7、W9、…の相対位置を特定した学習用位置情報として用いる。なお、周辺の建設地点W6、W7、W9、…の学習用位置情報は、選定地点を原点とした直交座標系の座標であってもよく、極座標系の座標であってもよい。なお、第2学習部の学習で利用される建設地点に、第1学習部の学習で利用される建設地点が含まれていてもよい。 For example, in FIG. 10, the construction point W1 on the map M2 of the predetermined area is selected as the selected point. Peripheral construction points W6, W7, W9, . . . within a predetermined range R1 (for example, distance r) from the selected construction point W1 are extracted. Then, the subsidence parameter value of the construction site W1 and the subsidence parameter values of each of the surrounding construction sites W6, W7, W9, . . . are used as learning subsidence parameter values. It is used as learning position information specifying relative positions of surrounding construction sites W6, W7, W9, . . . with respect to the construction site (selected site) W1. The learning position information of the surrounding construction sites W6, W7, W9, . The construction sites used in the learning of the second learning section may include the construction sites used in the learning of the first learning section.
これらを1組の教師データとして、図11に示すように、周辺の建設地点W6、W7、W9、…ごとの学習用沈下パラメータの値と、周辺の建設地点W6、W7、W9、…の学習用位置情報とを、ニューラルネットワーク13’に入力し、この出力値が、建設地点W1の学習用沈下パラメータの値に収束するように、沈下パラメータの値の算出を機械学習により学習する。このような学習を、他の複数の建設地点(例えば図10のW4、W36等)に対しても選定地点に選定することにより、教師データを抽出し、この教師データを用いて、沈下パラメータの値の算出を機械学習により学習する。 Using these as a set of training data, as shown in FIG. The positional information for construction is input to the neural network 13', and the calculation of the subsidence parameter value is learned by machine learning so that the output value converges to the value of the subsidence parameter for learning of the construction site W1. By selecting such learning as a selection point for a plurality of other construction points (for example, W4, W36, etc. in FIG. 10), teacher data is extracted, and this teacher data is used to calculate settlement parameters. Learning to calculate values by machine learning.
なお、第2学習部13のニューラルネットワーク13’も、第1学習部11のニューラルネットワーク11’と同様に、入力層13A、中間層13B~13D、および出力層13Eを備えている。これらの層は、ニューロン素子13a~13eで構成されており、ニューロンパラメータの値に乗算される重み付け係数を繰り返し補正することにより、機械学習を行う。
As with the neural network 11' of the
ここで、たとえば、各建設地点W1~W34において教師データに用いる沈下パラメータの値は、第1学習部11により算出される各建設地点の沈下パラメータの値であってもよく、各地点の実際の沈下パラメータの値であってもよい。ここで、教師データのうち、実際の沈下パラメータの値が格納された建設地点を、選定地点とすることが好ましい。これにより、実際の沈下パラメータの値に収束するように、第2学習部13の機械学習が行われるので、より精度の高い沈下パラメータの値を算出することができる。
Here, for example, the subsidence parameter values used in the training data at each of the construction sites W1 to W34 may be the subsidence parameter values of each construction site calculated by the
さらに、本実施形態では、第2学習部13は、選定した建設地点W1に対して、建設地点W1から所定の範囲R1内に存在する7個の建設地点W6、W7、W9、…を抽出して学習したが、他の建設地点を選択したり、所定の範囲の大きさを変更したりすることで、この周辺の建設地点が、他の複数の個数(7個以外の個数)に対して、同様の学習を行ってもよい。これにより、周辺情報抽出部22で抽出された既知地点の個数に合わせて、後述する第2学習部13で推定地点Gの沈下パラメータの値をより正確に算出することができる。
Furthermore, in the present embodiment, the
ここで、第2学習部13は、教師データとして、地域情報格納部12で格納された所定の地域内の既知地点T1~T33ごとの位置情報と、第1学習部11に算出された既知地点T1~T33ごとの沈下パラメータの値とを用いて機械学習したものであってもよい。これにより、第2学習部13は、推定地点Gを含む所定の地域において、沈下パラメータの値の算出を学習するので、推定地点Gの沈下パラメータの値の精度を高めることができる。
Here, the
本実施形態では、図12に示すように、沈下パラメータ推定部24は、周辺沈下パラメータ算出部23で算出した既知地点T1、T7、T8、…ごとの沈下パラメータの値D1、D7、T8、…と、推定地点に対する既知地点ごとの相対的な位置情報(X1,Y1)、(X7,Y7)、(X8,Y8)、…とを、第2学習部13に入力し、第2学習部13に推定地点Gの沈下パラメータの値を算出させる。
In this embodiment, as shown in FIG. 12, the subsidence
ここで、図12に示すように、既知地点T1、T7、T8、…ごとの位置情報は、推定地点Gを原点(0,0)とした直交系座標における座標であってもよく、推定地点Gを原点とした極座標系における座標であってもよい。 Here, as shown in FIG. 12, the position information for each of the known points T1, T7, T8, . Coordinates in a polar coordinate system with G as the origin may be used.
第2学習部13では、選定した建設地点W1に対して、建設地点W1から所定の範囲R1内に存在する建設地点W6、W7、W9、…の個数に応じて生成されたニューラルネットワークのモデルを選択し、このモデルを用いて、推定地点の沈下パラメータの値を算出する。
In the
試験実施判定部25は、周辺情報抽出部22で抽出した既知地点の数と、沈下パラメータ推定部24で算出した沈下パラメータの値とから、推定地点のボーリング試験および標準貫入試験の実施の要否を判定する。具体的には、周辺情報抽出部22で抽出した既知地点の数が、所定の数よりも少ない(たとえば、6カ所より少ない)場合には、抽出した点が少ないため、試験実施判定部25は、推定地点Gのボーリング試験および標準貫入試験の実施を要と判定する。
Based on the number of known points extracted by the surrounding
一方、周辺情報抽出部22で抽出した既知地点の数が、所定の数以上(たとえば、6カ所以上)である場合には、抽出した点が十分であるため、試験実施判定部25は、推定地点Gのボーリング試験および標準貫入試験の実施を否と判定する。たとえば、本実施形態では、既知地点の数が7カ所であるので、試験の実施を否と判定する。
On the other hand, if the number of known points extracted by the surrounding
さらに、試験の実施を否と判定した場合であっても、沈下パラメータ推定部24で算出した沈下パラメータの値が、所定の値を超えている(たとえば30を超えている)場合には、地盤沈下の可能性が高いため、試験実施判定部25は、推定地点Gのボーリング試験および標準貫入試験の実施を要と判定する。一方、沈下パラメータ推定部24で算出した沈下パラメータの値が、所定の値以下(たとえば30以下)である場合には、地盤沈下の可能性が低いため、試験実施判定部25は、推定地点Gのボーリング試験および標準貫入試験の実施を否と判定する。
Furthermore, even if it is determined that the test should not be performed, if the subsidence parameter value calculated by the subsidence
このように、周辺情報抽出部22で抽出した既知地点の数が、所定の数よりも少ない場合、または、沈下パラメータ推定部24で算出した沈下パラメータの値が、所定の値を超えている場合に、試験実施判定部25は、推定地点Gのボーリング試験および標準貫入試験の実施を要と判定する。これにより、推定地点Gでの地盤沈下の可能性を、より的確に予測することができる。
In this way, when the number of known points extracted by the surrounding
以下に、図13を参照して、本実施形態に係る地盤沈下予測システム1を用いた地盤沈下予測方法を説明する。ステップS1、S2では、図8に示すステップS1、S2と同様のことを行う。本実施形態では、ステップS21において、S1で機械学習を行った第1学習部11に対して、ステップS2で格納した既知地点のデータのうち、実際に地盤沈下が発生した既知地点に対して、その地盤情報と沈下パラメータとを教師データとし、第1学習部11による学習を再度行う。これにより、ステップS1で生成したモデルが補正され、その所定の地域の地盤沈下の傾向を、第1学習部11に再学習させることができる。
A ground subsidence prediction method using the ground
ステップS3~ステップS6では、図8に示すステップS3~ステップS6と同様のことを行う。次に、ステップS61では、沈下パラメータ推定部24は、図12に示すように、既知地点の座標を算出する。
In steps S3 to S6, the same operations as steps S3 to S6 shown in FIG. 8 are performed. Next, in step S61, the
ステップS81では、沈下パラメータ推定部24は、携帯端末20からネットワークを介して、各既知地点T1、T7、T8、…ごとの位置情報(X1,Y1)、(X7,Y7)、(X8,Y8)、…と、これに紐付けられた沈下パラメータの値D1、D7、D8、…と、をメインサーバ10の演算装置10Aに送信する。メインサーバ10内の第2学習部13では、推定地点Gの沈下パラメータの値を算出する。第2学習部13は、算出した結果を、ネットワークを介して携帯端末20の演算装置20Aに送信し、この結果が試験実施判定部25に入力される。
In step S81, the
次に、ステップS82では、試験実施判定部25が、抽出した既知地点の数と、算出した推定地点Gの沈下パラメータの値とから、推定地点Gのボーリング試験および標準貫入試験の実施の要否を判定する。最後に、ステップS9では、推定地点Gの沈下パラメータの値と試験実施の要否の結果を、携帯端末20の表示・入力部20Bに入力し、この結果が表示される。
Next, in step S82, the test
以上、本発明の実施形態について詳述したが、本発明は、前記の実施形態に限定されるものではなく、特許請求の範囲に記載された本発明の精神を逸脱しない範囲で、種々の設計変更を行うことができるものである。 Although the embodiments of the present invention have been described in detail above, the present invention is not limited to the above-described embodiments, and various designs can be made without departing from the spirit of the invention described in the claims. Changes can be made.
第2実施形態では、試験実施判定部は、抽出した既知地点の数と、算出した推定地点の沈下パラメータの値とから、推定地点のボーリング試験および標準貫入試験の実施の要否を判定した。たとえば、試験実施判定部は、推定地点から抽出した既知地点までの距離の平均値、抽出した既知地点の数、または算出した推定地点の沈下パラメータの値から選択される1つ以上を用いて、推定地点のボーリング試験および標準貫入試験の実施の要否を判定してもよい。さらに、第1学習部を用いて、所定の地域における地図データに紐付けた沈下パラメータの値を示すマップを作成してもよい。 In the second embodiment, the test execution determination unit determines whether or not to conduct a boring test and a standard penetration test at the estimated point based on the number of extracted known points and the calculated subsidence parameter value at the estimated point. For example, the test implementation determination unit uses one or more selected from the average value of the distance from the estimated point to the known point extracted, the number of extracted known points, or the calculated subsidence parameter value of the estimated point, It may be determined whether or not a boring test and a standard penetration test should be conducted at the estimated point. Furthermore, the first learning unit may be used to create a map indicating values of subsidence parameters linked to map data in a predetermined area.
1:地盤沈下予測システム、10:メインサーバ、20:携帯端末、11:沈下パラメータ学習部(第1学習部)、12:地域情報格納部、13:沈下推定学習部(第2学習部)、21:推定地点設定部、22:周辺情報抽出部、23:周辺沈下パラメータ算出部、24:沈下パラメータ推定部、25:試験実施判定部、G:推定地点、T1~T33:既知地点 1: Ground subsidence prediction system 10: Main server 20: Mobile terminal 11: Subsidence parameter learning unit (first learning unit) 12: Regional information storage unit 13: Subsidence estimation learning unit (second learning unit) 21: Estimated point setting unit, 22: Surrounding information extractor, 23: Surrounding subsidence parameter calculator, 24: Subsidence parameter estimator, 25: Test implementation determination unit, G: Estimated point, T1 to T33: Known points
Claims (5)
複数の既設建物の建設地点のボーリング試験により得られた、前記建設地点の地表からの深度、前記深度に応じた土質を数値化して設定された土質パラメータの値、および前記深度に応じた標準貫入試験におけるN値、を含む学習用地盤情報と、前記学習用地盤情報に対応した前記既設建物の地盤沈下量に基づいて、地盤沈下のし易さに応じて設定された学習用沈下パラメータの値と、を教師データとして、任意の地点の地盤情報から前記任意の地点の沈下パラメータの値の算出を機械学習により学習した第1学習部と、
所定の地域の地図データと、前記地図データの地図上の複数の既知地点において、前記既知地点ごとの位置情報、前記既知地点ごとのボーリング試験により得られた、地表からの深度、前記深度に応じた土質を数値化して設定された土質パラメータの値、および前記深度に応じた標準貫入試験におけるN値、を含む既知地点情報とが、格納された地域情報格納部と、
前記地図上に前記推定地点を設定する推定地点設定部と、
前記推定地点設定部で設定された前記推定地点の周辺に存在する前記既知地点ごとの前記既知地点情報を抽出する周辺情報抽出部と、
前記周辺情報抽出部で抽出した前記既知地点ごとの前記既知地点情報を、前記第1学習部に入力し、前記既知地点ごとの前記沈下パラメータの値を前記第1学習部に算出させる周辺沈下パラメータ算出部と、
前記周辺沈下パラメータ算出部に算出させた前記既知地点ごとの前記沈下パラメータに基づいて、前記推定地点の沈下パラメータの値を算出する沈下パラメータ推定部と、を備えることを特徴とする地盤沈下予測システム。 A ground subsidence prediction system that predicts ground subsidence at a specific point as an estimated point of ground subsidence,
The depth from the ground surface of the construction site, the value of the soil parameter set by quantifying the soil quality according to the depth, and the standard intrusion according to the depth, obtained by boring tests at the construction site of a plurality of existing buildings. A learning subsidence parameter value set according to the likelihood of ground subsidence based on the ground subsidence information for learning including the N value in the test and the amount of land subsidence of the existing building corresponding to the ground information for learning. and a first learning unit that learns by machine learning the calculation of the subsidence parameter value of the arbitrary point from the ground information of the arbitrary point, using as teacher data,
Map data of a predetermined area, location information for each known point at a plurality of known points on the map of the map data, depth from the ground surface obtained by a boring test for each known point, and depending on the depth a region information storage unit that stores known point information including a soil parameter value set by quantifying the soil quality and the N value in the standard penetration test according to the depth;
an estimated point setting unit that sets the estimated point on the map;
a peripheral information extraction unit for extracting the known point information for each of the known points existing around the estimated point set by the estimated point setting unit;
The known point information for each of the known points extracted by the surrounding information extraction unit is input to the first learning unit, and the first learning unit calculates the value of the subsidence parameter for each known point. a calculation unit;
a subsidence parameter estimator that calculates the value of the subsidence parameter of the estimated point based on the subsidence parameter for each of the known points calculated by the surrounding subsidence parameter calculator. .
前記地盤沈下予測システムは、前記周辺情報抽出部で抽出した前記既知地点の数と、前記沈下パラメータ推定部で算出した沈下パラメータの値とから、前記推定地点のボーリング試験および標準貫入試験の実施の要否を判定する試験実施判定部をさらに備えることを特徴とする請求項1に記載の地盤沈下予測システム。 The peripheral information extraction unit extracts the known point information for each of the known points existing within a predetermined distance from the estimated point,
The ground subsidence prediction system uses the number of known points extracted by the surrounding information extraction unit and the subsidence parameter values calculated by the subsidence parameter estimation unit to perform a boring test and a standard penetration test on the estimated points. 2. The ground subsidence prediction system according to claim 1, further comprising a test execution determination unit that determines necessity.
前記沈下パラメータ推定部は、前記周辺沈下パラメータ算出部で算出した前記既知地点ごとの前記沈下パラメータの値と、前記推定地点に対する前記既知地点ごとの相対的な位置情報と、を前記第2学習部に入力し、前記第2学習部に前記推定地点の沈下パラメータの値を算出させることを特徴とする請求項1または2に記載の地盤沈下予測システム。 The ground subsidence prediction system provides a learning subsidence parameter value of a selected point selected from a plurality of existing building construction sites, and a learning subsidence parameter value of a plurality of surrounding construction points within a predetermined range from the selected point. and positional information for learning that specifies the relative positional information of each of the construction sites in the surrounding area with respect to the selected site as one data group, and a plurality of data groups obtained for each of the different selected sites as teacher data. Calculating the value of the subsidence parameter of the arbitrary point from the values of the subsidence parameters of a plurality of surrounding points existing around the arbitrary point and the relative position information of each of the surrounding points with respect to the arbitrary point further comprising a second learning unit learned by machine learning,
The subsidence parameter estimating unit calculates the value of the subsidence parameter for each of the known points calculated by the surrounding subsidence parameter calculating unit, and relative positional information of each known point with respect to the estimated point, to the second learning unit. 3. The land subsidence prediction system according to claim 1 or 2, wherein the subsidence parameter value of the estimated point is calculated by the second learning unit.
The subsidence parameter estimating unit calculates the subsidence parameter of the estimated point based on the value of the subsidence parameter for each of the known points calculated by the surrounding subsidence parameter calculating unit and the distance to each known point extracted from the estimated point. 3. The ground subsidence prediction system according to claim 1 or 2, wherein the value of is calculated.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019122228A JP7246272B2 (en) | 2019-06-28 | 2019-06-28 | Ground subsidence prediction system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019122228A JP7246272B2 (en) | 2019-06-28 | 2019-06-28 | Ground subsidence prediction system |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2021008729A JP2021008729A (en) | 2021-01-28 |
JP7246272B2 true JP7246272B2 (en) | 2023-03-27 |
Family
ID=74199480
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2019122228A Active JP7246272B2 (en) | 2019-06-28 | 2019-06-28 | Ground subsidence prediction system |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP7246272B2 (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112883646B (en) * | 2021-02-20 | 2023-11-24 | 首都师范大学 | Building settlement amount extraction method, system and device combining machine learning and soil mechanics model |
CN115419120A (en) * | 2022-06-08 | 2022-12-02 | 山东大学 | Highway subgrade settlement monitoring and predicting method |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006170739A (en) | 2004-12-15 | 2006-06-29 | Kajima Corp | Earthquake disaster prevention system using urgent earthquake prompt report |
JP2018194404A (en) | 2017-05-16 | 2018-12-06 | 株式会社パスコ | Machine learning method, and surface deformation determination method |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000002769A (en) * | 1998-06-17 | 2000-01-07 | Nippon Koei Co Ltd | Method and device for predicting the spatial distribution of geological structure to prepare geological map |
JP5516211B2 (en) * | 2010-08-09 | 2014-06-11 | 株式会社大林組 | Ground estimation method and ground estimation system |
JP2019065523A (en) * | 2017-09-29 | 2019-04-25 | 積水化学工業株式会社 | Ground investigation method |
-
2019
- 2019-06-28 JP JP2019122228A patent/JP7246272B2/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006170739A (en) | 2004-12-15 | 2006-06-29 | Kajima Corp | Earthquake disaster prevention system using urgent earthquake prompt report |
JP2018194404A (en) | 2017-05-16 | 2018-12-06 | 株式会社パスコ | Machine learning method, and surface deformation determination method |
Also Published As
Publication number | Publication date |
---|---|
JP2021008729A (en) | 2021-01-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shahin | State-of-the-art review of some artificial intelligence applications in pile foundations | |
Oh et al. | Neural network-based seismic response prediction model for building structures using artificial earthquakes | |
Shahin | Load–settlement modeling of axially loaded steel driven piles using CPT-based recurrent neural networks | |
Kuo et al. | ANN-based model for predicting the bearing capacity of strip footing on multi-layered cohesive soil | |
Hashash et al. | Comparison of two inverse analysis techniques for learning deep excavation response | |
Aydin | Fuzzy set approaches to classification of rock masses | |
Shahin | Use of evolutionary computing for modelling some complex problems in geotechnical engineering | |
JP7097274B2 (en) | Ground prediction system | |
JP7121891B2 (en) | Structural information display device, structural design support device, and structural design support model learning device | |
Khoshnevisan et al. | Robust design in geotechnical engineering–an update | |
JP7246272B2 (en) | Ground subsidence prediction system | |
Peng et al. | Reliability-based robust geotechnical design using Monte Carlo simulation | |
Fontan et al. | Soil–structure interaction: Parameters identification using particle swarm optimization | |
JP2019179538A (en) | Prediction method for predicting construction price of building and prediction method for predicting its construction period | |
Díaz et al. | Use of artificial neural networks to predict 3-D elastic settlement of foundations on soils with inclined bedrock | |
Hölter et al. | Optimal measurement design for parameter identification in mechanized tunneling | |
Harandizadeh | Developing a new hybrid soft computing technique in predicting ultimate pile bearing capacity using cone penetration test data | |
EP4123287A1 (en) | N-value prediction device and method using data augmentation-based artificial intelligence | |
Wang et al. | Robust optimization of the constructional time delay in the design of double-row stabilizing piles | |
CN109711063B (en) | Station underground excavation footage adjusting method and device suitable for soft upper and hard lower strata | |
Barari et al. | Determination of the ultimate limit states of shallow foundations using gene expression programming (GEP) approach | |
JP5516211B2 (en) | Ground estimation method and ground estimation system | |
Vargas et al. | Incremental dynamic analysis and pushover analysis of buildings. A probabilistic comparison | |
Monsalve-Giraldo et al. | Hybrid Parabolic Interpolation–Artificial Neural Network Method (HPI-ANNM) for long-term extreme response estimation of steel risers | |
JP4206980B2 (en) | Sinking amount prediction method and settlement amount prediction program |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20220602 |
|
TRDD | Decision of grant or rejection written | ||
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20230216 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20230221 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20230314 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7246272 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |