JP6518939B2 - Method of estimating moisture of concrete structure and system of estimating moisture of concrete structure - Google Patents

Method of estimating moisture of concrete structure and system of estimating moisture of concrete structure Download PDF

Info

Publication number
JP6518939B2
JP6518939B2 JP2015123727A JP2015123727A JP6518939B2 JP 6518939 B2 JP6518939 B2 JP 6518939B2 JP 2015123727 A JP2015123727 A JP 2015123727A JP 2015123727 A JP2015123727 A JP 2015123727A JP 6518939 B2 JP6518939 B2 JP 6518939B2
Authority
JP
Japan
Prior art keywords
concrete structure
water
moisture
response
concrete
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
JP2015123727A
Other languages
Japanese (ja)
Other versions
JP2017009371A (en
Inventor
光史 浅見
光史 浅見
和哉 西村
和哉 西村
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National Institute of Maritime Port and Aviation Technology
Original Assignee
National Institute of Maritime Port and Aviation Technology
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 National Institute of Maritime Port and Aviation Technology filed Critical National Institute of Maritime Port and Aviation Technology
Priority to JP2015123727A priority Critical patent/JP6518939B2/en
Publication of JP2017009371A publication Critical patent/JP2017009371A/en
Application granted granted Critical
Publication of JP6518939B2 publication Critical patent/JP6518939B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

本発明は、コンクリート構造物の水分推定方法及びコンクリート構造物の水分推定システムに関する。   The present invention relates to a method of estimating moisture of a concrete structure and a system of estimating moisture of a concrete structure.

コンクリート構造物にはジャンカやひび割れ等の欠損が生ずる可能性があり、そこに水が溜まってコンクリート構造物を構成する鉄筋や構造材を腐食したり脆弱化するおそれがある。したがって、コンクリート構造物内の水分の分布を推定できる技術の実現が望まれている。   In concrete structures, defects such as junk and cracks may occur, and water may be accumulated therein to corrode or weaken reinforcing bars and structural materials that constitute the concrete structure. Therefore, it is desirable to realize a technology that can estimate the distribution of water in a concrete structure.

例えば、コンクリート構造物を挟むように中性子線源とパネル型高感度中性子検出器を配置し、コンクリート構造物を透過させて中性子線源から放射された中性子を検出器によって検出することによりコンクリート構造物の欠損を検出する技術が開示されている(特許文献1)。   For example, a neutron beam source and a panel type high sensitivity neutron detector are disposed so as to sandwich the concrete structure, and the concrete structure is transmitted by the concrete structure and the neutron emitted from the neutron beam source is detected by the detector. Patent Document 1 discloses a technique for detecting a defect of

また、セメントを硬化材としてホウ素化合物等のトレーサを添加した注入剤をセメントの裏側に注入し、セメント表面側から中性子を照射してトレーサに救出された熱中性子の数を測定してセメントの注入度合を検出する技術が開示されている(特許文献2)。   Also, inject the cement into the back side of the cement with an injection agent added with a tracer such as a boron compound as a hardening material, irradiate neutrons from the cement surface side, measure the number of thermal neutrons rescued by the tracer, and inject cement. A technique for detecting the degree is disclosed (Patent Document 2).

特開2002−82073号公報JP 2002-82073 A 特開平8−94552号公報JP-A-8-94552

ところで、照射する中性子を発生させる中性子線源は、原子炉を使用する方式と加速器を使用する方式がある。透過型の中性子線検出では中性子線源とセンサとで測定対象物を挟み込む配置が必要であり、主として屋外に存在するコンクリート構造物を測定対象物とする場合には中性子線源やセンサの移動が不可能又は困難であったり、測定対象物の形状や厚みによっては背面側への配置が不可能であったりする。また、中性子線源と検出器との配置の関係において原理的にシステムが大型化してしまう。最近では屋外でも使用できる移動可能な加速器を開発している研究機関もあるが開発途上である。   By the way, as a neutron source for generating neutrons to be irradiated, there are a method using a nuclear reactor and a method using an accelerator. In the transmission type neutron beam detection, it is necessary to arrange the object to be measured between the neutron beam source and the sensor, and when using a concrete structure existing mainly outdoors as the object to be measured, movement of the neutron beam source or sensor is required. It is impossible or difficult, and depending on the shape and thickness of the object to be measured, placement on the back side may not be possible. In addition, in principle, the size of the system is increased due to the arrangement of the neutron source and the detector. Recently, some research institutes are developing mobile accelerators that can be used outdoors, but they are still developing.

また、ホウ素化合物等をトレーサとする技術は、トレーサが添加されていない既存のコンクリート構造物には適用することが原理的に不可能である。   Moreover, it is impossible in principle to apply the technique which makes a boron compound etc. a tracer to the existing concrete structure to which the tracer is not added.

本発明は、コンクリート構造物内のジャンカやひび割れ等に溜まった水分を検知するコンクリート構造物の水分推定方法及びコンクリート構造物の水分推定システムを提供する。   The present invention provides a method for estimating the moisture of a concrete structure and a system for estimating the moisture of the concrete structure, which detect moisture accumulated in junka, cracks, etc. in the concrete structure.

本発明の請求項1に係る中性子を用いたコンクリート構造物の水分推定方法は、コンクリートの元素密度組成及び結合水の含有量を条件として、中性子の線源の位置と計測点の位置との関係に基づいて前記線源から放射された中性子が前記コンクリートにより構成された調査対象であるコンクリート構造物を介して前記計測点において計測されるときの応答をモデル化してモンテカルロ法及びSN法を含む放射線輸送計算コードにより解析する第1のステップと、前記線源から放射された中性子が前記調査対象であるコンクリート構造物を介して前記計測点において計測されるときの応答を実測する第2のステップと、前記実測された応答と前記解析された応答とを用いて前記調査対象であるコンクリート構造物内の水分の存在を推定する第3のステップとを備えることを特徴とする。   In the method of estimating moisture of a concrete structure using neutrons according to claim 1 of the present invention, the relationship between the position of the radiation source of neutrons and the position of the measurement point, on the condition of the element density composition of concrete and the content of bound water Modeling the response when neutrons emitted from the radiation source are measured at the measurement point through the concrete structure composed of the concrete on the basis of the radiation including the Monte Carlo method and the SN method A first step of analyzing by a transport calculation code, and a second step of measuring a response when neutrons emitted from the radiation source are measured at the measurement point through the concrete structure to be investigated Estimating the presence of water in the concrete structure to be investigated using the measured response and the analyzed response Characterized in that it comprises of the steps.

ここで、前記第3のステップは、前記第1のステップにおいて複数の前記元素密度組成及び前記結合水の含有率の組み合わせであるコンクリートの種別毎に解析した前記解析された応答を記憶させたデータベースを参照して、前記実測された応答と前記調査対象であるコンクリート構造物を構成するコンクリートの種別に対して前記解析された応答とを用いて前記調査対象であるコンクリート構造物内の水分の存在を推定することが好適である。   Here, the third step is a database storing the analyzed response analyzed for each type of concrete, which is a combination of a plurality of the element density compositions and the content of the binding water in the first step. The presence of moisture in the concrete structure to be investigated using the measured response and the analyzed response to the type of concrete constituting the concrete structure to be investigated. It is preferable to estimate

また、前記第1のステップは、前記調査対象であるコンクリート構造物にそれぞれ検査対象領域及び周辺領域を設定してモデル化し、前記検査対象領域をm個の分割対象領域k(mは2以上の整数、k=1,…,m)にメッシュ分割して、複数の前記線源の位置i(iは2以上の整数)及び複数の前記計測点の位置j(jは2以上の整数)の組み合わせの各々についての前記分割対象領域kの各々からの寄与ri,j,kを求め、前記線源の位置iから放射された中性子が前記調査対象であるコンクリート構造物を介して前記計測点の位置jにおいて計測されるときの前記解析された応答Ri,j

Figure 0006518939
(ただし、Ari,jは、前記線源の位置iから放射された中性子が前記調査対象であるコンクリート構造物を介して前記計測点の位置jにおいて計測されるときの前記周辺領域からの応答の寄与)により求めることが好適である。 Further, in the first step, the inspection target area and the peripheral area are respectively set and modeled in the concrete structure to be investigated, and the inspection target area is divided into m division target areas k (m is 2 or more). The mesh is divided into integers, k = 1,..., M), and the positions i (i is an integer of 2 or more) of the plurality of radiation sources and the positions j (j The contribution ri , j, k from each of the division target regions k for each of the combinations is determined, and the neutrons emitted from the position i of the radiation source are measured via the concrete structure to be investigated The analyzed response R i, j when measured at position j of
Figure 0006518939
(Where Ar i, j is the response from the surrounding area when neutrons emitted from the position i of the radiation source are measured at the position j of the measurement point through the concrete structure to be investigated) It is preferable to obtain by

また、前記第1のステップは、前記線源の位置i及び前記計測点の位置jの組み合わせの各々について、前記分割対象領域kに前記コンクリートが体積率100%で存在するときの寄与rci,j,k、及び、前記分割対象領域kに前記水が体積率100%で存在するときの寄与rwi,j,kを求め、前記第2のステップは、前記線源の位置i及び前記計測点の位置jの組み合わせの各々について、前記線源から放射された中性子が前記調査対象であるコンクリート構造物を介して前記計測点において計測されるときの応答Di,jを実測し、前記第3のステップは、

Figure 0006518939
(ただし、Sは、前記線源の強度)
の関係を用いて前記分割対象領域kにおける水分の体積率ωを推定することが好適である。 In the first step, for each combination of the position i of the radiation source and the position j of the measurement point, contribution rc i when the concrete is present at a volume ratio of 100% in the division target area k j, k , and a contribution r w i, j, k when the water is present at a volume ratio of 100% in the division target area k , and the second step determines the position i of the radiation source and the measurement For each combination of point positions j, the response D i, j is measured when neutrons emitted from the radiation source are measured at the measurement point via the concrete structure to be investigated, Step 3 is
Figure 0006518939
(Where S is the intensity of the radiation source)
It is preferable to estimate the volume fraction ω k of water in the division target region k using the relationship of

また、前記数式(2)の誤差が最小となる前記分割対象領域kにおける水分の体積率
ωを推定することが好適である。
Further, it is preferable to estimate the volume fraction ω k of water in the division target region k which minimizes the error of the equation (2).

また、前記第3のステップにおいて推定された前記調査対象であるコンクリート構造物内の水分の存在を可視化する第4のステップを備えることが好適である。   In addition, it is preferable to include a fourth step of visualizing the presence of water in the concrete structure to be investigated, which is estimated in the third step.

また、前記第2のステップにおいて使用される前記線源は、前記調査対象であるコンクリート構造物に対面する領域以外を中性子遮蔽手段で覆った構造を有することが好適である。   Further, it is preferable that the radiation source used in the second step has a structure in which a region other than a region facing the concrete structure to be investigated is covered with a neutron shielding means.

また、前記第2のステップは、前記調査対象であるコンクリート構造物の表面上を前記線源及び前記計測点を順次移動させて前記実測された応答を実測することが好適である。   In the second step, it is preferable to move the radiation source and the measurement point sequentially on the surface of the concrete structure to be examined to measure the actually measured response.

本発明の請求項9に係る中性子を用いたコンクリート構造物の水分推定システムは、コンクリートの元素密度組成及び結合水の含有量を条件として、中性子の線源の位置と計測点の位置との関係に基づいて前記線源から放射された中性子が前記コンクリートにより構成された調査対象であるコンクリート構造物を介して前記計測点において計測されるときの応答をモデル化して解析した結果をデータベースとして記憶する記憶手段と、前記線源から放射された中性子が調査対象であるコンクリート構造物を介して前記計測点において計測されるときの応答を実測する実応答計測手段と、前記実測された応答と前記解析された応答とを用いて前記調査対象であるコンクリート構造物内の水分の存在を推定する水分推定手段とを備えることを特徴とする。   According to the ninth aspect of the present invention, there is provided a system for estimating the moisture content of a concrete structure using neutrons, wherein the relationship between the position of the radiation source of neutrons and the position of the measurement point is conditional on the element density composition of concrete and the content of bound water. And store the results of modeling and analyzing the response when neutrons emitted from the radiation source are measured at the measurement point through the concrete structure formed by the concrete based on the above. Storage means, actual response measurement means for measuring the response when neutrons emitted from the radiation source are measured at the measurement point through the concrete structure to be investigated, the measured response, and the analysis Water estimation means for estimating the presence of water in the concrete structure to be investigated using the determined response To.

ここで、前記水分推定手段において推定された前記調査対象であるコンクリート構造物内の水分の存在を可視化して表示する表示手段をさらに備えることが好適である。   Here, it is preferable to further include display means for visualizing and displaying the presence of the water in the concrete structure to be surveyed estimated by the water estimation means.

本発明の請求項1に係る中性子を用いたコンクリート構造物の水分推定方法によれば、コンクリートの元素密度組成及び結合水の含有量を条件として、中性子の線源の位置と計測点の位置との関係に基づいて前記線源から放射された中性子が前記コンクリートにより構成された調査対象であるコンクリート構造物を介して前記計測点において計測されるときの応答をモデル化してモンテカルロ法及びSN法を含む放射線輸送計算コードにより解析する第1のステップと、前記線源から放射された中性子が前記調査対象であるコンクリート構造物を介して前記計測点において計測されるときの応答を実測する第2のステップと、前記実測された応答と前記解析された応答とを用いて前記調査対象であるコンクリート構造物内の水分の存在を推定する第3のステップとを備えることによって、前記評価対象であるコンクリート構造物内の水分の分布を精度よく推定することができる。また、コンクリート構造物の一方向からの水分の存在の推定が可能となり、トレーサが不要で既存のコンクリート構造物への適用も可能となる。   According to the method of estimating the moisture content of a concrete structure using neutrons according to claim 1 of the present invention, the position of the radiation source of the neutron and the position of the measurement point, on the condition of the element density composition of the concrete and the content of the bound water Modeling the response when neutrons emitted from the radiation source are measured at the measurement point through the concrete structure to be surveyed composed of the concrete based on the relationship between the Monte Carlo method and the SN method And a second step of measuring a response when neutrons emitted from the radiation source are measured at the measurement point through the concrete structure to be investigated. Estimating the presence of moisture in the concrete structure to be investigated using the step, the measured response and the analyzed response Third by providing the steps of the distribution of moisture in the evaluation concrete structure which is the subject can be estimated accurately that. In addition, it becomes possible to estimate the presence of water from one direction of the concrete structure, and it becomes possible to apply the present invention to an existing concrete structure without requiring a tracer.

また、前記第3のステップは、前記第1のステップにおいて複数の前記元素密度組成及び前記結合水の含有率の組み合わせであるコンクリートの種別毎に解析した前記解析された応答を記憶させたデータベースを参照して、前記実測された応答と前記調査対象であるコンクリート構造物を構成するコンクリートの種別に対して前記解析された応答とを用いて前記調査対象であるコンクリート構造物内の水分の存在を推定することによって、予め前記データベースとして登録された前記解析された応答に基づいて前記実測された応答に対応する前記コンクリート構造物内の水分の存在を推定することができる。   In the third step, a database storing the analyzed response analyzed for each type of concrete, which is a combination of a plurality of the element density compositions and the content of the binding water in the first step Referring to the measured response and the analyzed response to the type of concrete constituting the concrete structure to be surveyed, the presence of water in the concrete structure to be surveyed is determined. By estimating, the presence of water in the concrete structure corresponding to the measured response can be estimated based on the analyzed response registered in advance as the database.

また、前記第1のステップは、前記調査対象であるコンクリート構造物にそれぞれ検査対象領域及び周辺領域を設定してモデル化し、前記検査対象領域をm個の分割対象領域k(mは2以上の整数、k=1,…,m)にメッシュ分割して、複数の前記線源の位置i(iは2以上の整数)及び複数の前記計測点の位置j(jは2以上の整数)の組み合わせの各々についての前記分割対象領域kの各々からの寄与ri,j,kを求め、前記線源の位置iから放射された中性子が前記調査対象であるコンクリート構造物を介して前記計測点の位置jにおいて計測されるときの前記解析された応答Ri,jを前記数式(1)により求めることにより、前記コンクリート構造物を分割した分割対象領域kにおける水分の体積比を求めることができる。 Further, in the first step, the inspection target area and the peripheral area are respectively set and modeled in the concrete structure to be investigated, and the inspection target area is divided into m division target areas k (m is 2 or more). The mesh is divided into integers, k = 1,..., M), and the positions i (i is an integer of 2 or more) of the plurality of radiation sources and the positions j (j The contribution ri , j, k from each of the division target regions k for each of the combinations is determined, and the neutrons emitted from the position i of the radiation source are measured via the concrete structure to be investigated The volume ratio of water in the division target area k obtained by dividing the concrete structure can be determined by determining the analyzed response R i, j when measured at the position j of .

また、前記第1のステップは、前記線源の位置i及び前記計測点の位置jの組み合わせの各々について、前記分割対象領域kに前記コンクリートが体積率100%で存在するときの寄与rci,j,k及び、前記分割対象領域kに前記水が体積率100%で存在するときの寄与rwi,j,kを求め、前記第2のステップは、前記線源の位置i及び前記計測点の位置jの組み合わせの各々について、前記線源から放射された中性子が前記調査対象であるコンクリート構造物を介して前記計測点において計測されるときの応答Di,jを実測し、前記第3のステップは、上記数式(2)の関係を用いて前記分割対象領域kにおける水分の体積率ωを推定することによって、前記寄与rci,j,k及び前記寄与rwi,j,kに基づいて前記分割対象領域kにおける水分の体積率ωを推定することができる。 In the first step, for each combination of the position i of the radiation source and the position j of the measurement point, contribution rc i when the concrete is present at a volume ratio of 100% in the division target area k j, k and contributions r w i, j, k when the water is present at a volume ratio of 100% in the division target area k are obtained, and in the second step, the position i of the radiation source and the measurement point The response D i, j when neutrons emitted from the radiation source are measured at the measurement point through the concrete structure to be examined is measured for each of the combinations of the position j of In the step of, by estimating the volume ratio ω k of the water in the division target area k using the relationship of the equation (2), the contribution rc i, j, k and the contribution r w i, j, k can be obtained . On the basis of It is possible to estimate the volume fraction omega k of the moisture in the serial division target region k.

また、前記数式(2)の誤差が最小となる前記分割対象領域kにおける水分の体積率ωを推定することによって、前記数式(2)から導かれる連立方程式の完全解を解くことなく、前記寄与rci,j,k及び前記寄与rwi,j,kに基づいて前記分割対象領域kにおける水分の体積率ωを推定することができる。 Further, by estimating the volume ratio ω k of water in the division target region k which minimizes the error of the equation (2), the solution of the simultaneous equations derived from the equation (2) is not solved. Based on the contribution rc i, j, k and the contribution r w i, j, k , it is possible to estimate the volume fraction ω k of water in the division target area k.

また、前記第3のステップにおいて推定された前記調査対象であるコンクリート構造物内の水分の存在を可視化する第4のステップを備えることによって、推定された水分の分布をユーザの視覚によって把握し易い態様で表示させることができる。   Further, by providing a fourth step of visualizing the presence of water in the concrete structure to be surveyed estimated in the third step, the distribution of the estimated water can be easily grasped by the user's vision It can be displayed in an aspect.

また、前記第2のステップにおいて使用される前記線源は、前記調査対象であるコンクリート構造物に対面する領域以外を中性子遮蔽手段で覆った構造を有することによって、中性子線源から前記コンクリート構造物内に入射・散乱せずに中性子検出器106に直接到達する中性子を低減することができる。   Further, the radiation source used in the second step has a structure in which the neutron shielding means covers a region other than the region facing the concrete structure to be investigated, so that the concrete structure from the neutron beam source Neutrons that directly reach the neutron detector 106 without being incident and scattered therein can be reduced.

また、前記第2のステップは、前記調査対象であるコンクリート構造物の表面上を前記線源及び前記計測点を順次移動させて前記実測された応答を実測することが好適である。   In the second step, it is preferable to move the radiation source and the measurement point sequentially on the surface of the concrete structure to be examined to measure the actually measured response.

本発明の請求項9に係る中性子を用いたコンクリート構造物の水分推定システムは、コンクリートの元素密度組成及び結合水の含有量を条件として、中性子の線源の位置と計測点の位置との関係に基づいて前記線源から放射された中性子が前記コンクリートにより構成された調査対象であるコンクリート構造物を介して前記計測点において計測されるときの応答をモデル化して解析した結果をデータベースとして記憶する記憶手段と、前記線源から放射された中性子が調査対象であるコンクリート構造物を介して前記計測点において計測されるときの応答を実測する実応答計測手段と、前記実測された応答と前記解析された応答とを用いて前記調査対象であるコンクリート構造物内の水分の存在を推定する水分推定手段とを備えることによって、前記評価対象であるコンクリート構造物内の水分の分布を精度よく推定することができる。また、コンクリート構造物の一方向からの水分の存在の推定が可能となり、トレーサが不要で既存のコンクリート構造物への適用も可能となる。   According to the ninth aspect of the present invention, there is provided a system for estimating the moisture content of a concrete structure using neutrons, wherein the relationship between the position of the radiation source of neutrons and the position of the measurement point is conditional on the element density composition of concrete and the content of bound water. And store the results of modeling and analyzing the response when neutrons emitted from the radiation source are measured at the measurement point through the concrete structure formed by the concrete based on the above. Storage means, actual response measurement means for measuring the response when neutrons emitted from the radiation source are measured at the measurement point through the concrete structure to be investigated, the measured response, and the analysis And means for estimating the presence of moisture in the concrete structure to be investigated using the determined response. , The distribution of moisture in the evaluation in the concrete structure is the object can be estimated accurately. In addition, it becomes possible to estimate the presence of water from one direction of the concrete structure, and it becomes possible to apply the present invention to an existing concrete structure without requiring a tracer.

ここで、前記水分推定手段において推定された前記調査対象であるコンクリート構造物内の水分の存在を可視化して表示する表示手段をさらに備えることによって、推定された水分の分布をユーザの視覚によって把握し易い態様で表示させることができる。   Here, by further comprising display means for visualizing and displaying the presence of the water in the concrete structure to be investigated estimated by the water estimation means, the distribution of the estimated water can be grasped by the user's vision It can be displayed in an easy manner.

本発明の実施の形態における水分推定システムの構成を示す図である。It is a figure which shows the structure of the water | moisture content estimation system in embodiment of this invention. 本発明の実施の形態における水分推定装置の構成を示す図である。It is a figure which shows the structure of the water | moisture content estimation apparatus in embodiment of this invention. 本発明の実施の形態における中性子線源の構成を示す模式図である。It is a schematic diagram which shows the structure of the neutron beam source in embodiment of this invention. 本発明の実施の形態における中性子の散乱モデルを示す図である。It is a figure which shows the scattering model of the neutron in embodiment of this invention. 本発明の実施の形態におけるコンクリート構造物の分割対象領域の設定例を示す図である。It is a figure which shows the example of a setting of the division | segmentation object area | region of the concrete structure in embodiment of this invention. 本発明の実施の形態における中性子の応答を解析するときの散乱モデルを示す図である。It is a figure which shows the scattering model when analyzing the response of the neutron in embodiment of this invention. 本発明の実施の形態における水分推定方法のフローチャートである。It is a flowchart of the water | moisture-content estimation method in embodiment of this invention.

本発明の実施の形態における中性子を用いたコンクリート構造物の水分推定システム100は、図1に示すように、水分推定装置102、中性子線源104及び中性子検出器106を含んで構成される。中性子線源104及び中性子検出器106は、調査対象であるコンクリート構造物200の同一表面側に配置されている。   The water estimation system 100 for concrete structures using neutrons according to the embodiment of the present invention includes a water estimation device 102, a neutron beam source 104, and a neutron detector 106, as shown in FIG. The neutron beam source 104 and the neutron detector 106 are disposed on the same surface side of the concrete structure 200 to be investigated.

水分推定装置102は、図2に示すように、一般的なコンピュータの基本構成である処理部10、記憶部12、入力部14、出力部16及びインターフェース部18を含んで構成される。   The water | moisture-content estimation apparatus 102 is comprised including the process part 10 which is the basic composition of a common computer, the memory | storage part 12, the input part 14, the output part 16, and the interface part 18 as shown in FIG.

処理部10は、CPU等の演算手段を含み、記憶部12に記憶された水分推定プログラムを実行することによって水分推定装置102において行われるコンクリート構造物200内の水分の推定処理を実現する。記憶部12は、処理部10からアクセス可能に接続された半導体メモリやハードディスク等の記憶装置を含んで構成される。記憶部12は、水分推定プログラム、各種パラメータ、測定データ等を記憶する。入力部14は、キーボードやマウス等の入力装置を含んで構成される。入力部14は、水分推定装置102によって使用される各種パラメータ等の入力に用いられる。出力部16は、ディスプレイやプリンタ等の出力装置を含んで構成される。出力部16は、水分推定装置102において得られた処理結果を可視化して表示する表示手段として機能する。出力部16は、水分推定装置102における処理に関する情報をユーザに提示する。インターフェース部18は、アナログ/デジタル変換器等を含んで構成される。インターフェース部18は、中性子検出器106で検出された中性子線の線量を示すデータを水分推定装置102に取り込むために用いられる。   The processing unit 10 includes computing means such as a CPU, and implements the moisture estimation process in the concrete structure 200 performed in the moisture estimation device 102 by executing the moisture estimation program stored in the storage unit 12. The storage unit 12 is configured to include storage devices such as a semiconductor memory and a hard disk that are connected so as to be accessible from the processing unit 10. The storage unit 12 stores a moisture estimation program, various parameters, measurement data, and the like. The input unit 14 is configured to include an input device such as a keyboard and a mouse. The input unit 14 is used to input various parameters and the like used by the moisture estimation device 102. The output unit 16 is configured to include an output device such as a display or a printer. The output unit 16 functions as a display unit that visualizes and displays the processing result obtained by the water estimation device 102. The output unit 16 presents, to the user, information related to the processing in the moisture estimation device 102. The interface unit 18 includes an analog / digital converter and the like. The interface unit 18 is used to load data indicating the dose of neutrons detected by the neutron detector 106 into the moisture estimation device 102.

中性子線源104は、中性子発生管や252Cf等のRI線源等の中性子発生源を含んで構成される。中性子線源104は、図3の断面図に示すように、中性子発生源20の周りを遮蔽することが好適である。遮蔽は、調査対象物であるコンクリート構造物200に向く方向以外、すなわちコンクリート構造物200に対面する領域以外を覆うように設けられる。中性子発生源20を遮蔽することによって、中性子線源104を調査対象であるコンクリート構造物200の表面に設置した場合に中性子やγ線がコンクリート構造物200内に入射及び散乱せず、中性子検出器106に中性子やγ線が直接到達することを抑制することができる。 The neutron source 104 includes a neutron source such as a neutron generator tube or an RI source such as 252 Cf. The neutron beam source 104 is preferably shielded around the neutron source 20, as shown in the cross-sectional view of FIG. The shielding is provided so as to cover the area other than the direction toward the concrete structure 200 which is the investigation target, that is, the area other than the area facing the concrete structure 200. By shielding the neutron source 20, when the neutron beam source 104 is installed on the surface of the concrete structure 200 to be investigated, neutrons and gamma rays do not enter and scatter in the concrete structure 200, and a neutron detector Direct arrival of neutrons or gamma rays at 106 can be suppressed.

例えば、複数種の遮蔽壁22,24,26を設けることが好適である。遮蔽壁22は、高速中性子を熱化させるために設けられる。遮蔽壁22は、例えばパラフィンを含む材料で構成される。遮蔽壁24は、熱化された中性子を吸収するために設けられる。遮蔽壁24は、カドミウムを含む材料で構成される。遮蔽壁26は、γ線を吸収するために設けられる。遮蔽壁26は、例えば鉛を含む材料で構成される。遮蔽壁22〜26は、例えば、それぞれ25cm、1mm、5cmの厚さとすればよい。   For example, it is preferable to provide multiple types of shielding walls 22, 24, 26. The shielding wall 22 is provided to heat fast neutrons. The shielding wall 22 is made of, for example, a material containing paraffin. The shielding wall 24 is provided to absorb the thermalized neutrons. The shielding wall 24 is made of a material containing cadmium. The shielding wall 26 is provided to absorb gamma rays. The shielding wall 26 is made of, for example, a material containing lead. The shielding walls 22 to 26 may have thicknesses of, for example, 25 cm, 1 mm, and 5 cm, respectively.

中性子検出器106は、中性子線源104から放出され、コンクリート構造物200によって反射又は散乱された中性子を検出する。中性子検出器106は、He−3計数管、BF3検出器等とすることができる。中性子検出器106は、中性子を検出すると、インターフェース部18を介して水分推定装置102へ出力する。水分推定装置102の処理部10は、中性子検出器106の検出信号を受けると単位時間当たりの中性子の検出数をカウントする。これにより、中性子線の線量率(計数率)が測定される。   The neutron detector 106 detects neutrons emitted from the neutron beam source 104 and reflected or scattered by the concrete structure 200. The neutron detector 106 can be a He-3 counter, a BF3 detector, or the like. When the neutron detector 106 detects neutrons, the neutron detector 106 outputs the neutrons to the water estimation device 102 through the interface unit 18. When receiving the detection signal of the neutron detector 106, the processing unit 10 of the moisture estimation device 102 counts the number of detected neutrons per unit time. Thereby, the dose rate (counting rate) of the neutron beam is measured.

本実施の形態における水分推定システム100では、図4に示すように、中性子線源104から放射される中性子をコンクリート構造物200の1つの表面側から入射させ、コンクリート構造物200から反射された中性子を同一の表面側に設けた中性子検出器106で検出することによってコンクリート構造物200内の水分の分布を評価する。このため、コンクリート構造物200の同一表面側、すなわち一方向からの水分の分布の評価が可能となる。中性子線源104から発生した高速中性子は、コンクリート構造物200やジャンカX等で反射又は散乱されて中性子検出器106に到達する。このとき、ジャンカX等に水が溜まった「水溜まり」が存在すると、中性子は水内の水素によって減速され、「水溜まり」が存在しない場合に比べて熱中性子化され易くなる。この核反応の違いに着目することによってコンクリート構造物200内の水分の分布を評価する。水分推定システム100によれば、コンクリート構造物200の表面から中性子の数平均自由行程程度の深さまでの水分の分布を推定することができる。   In the water estimation system 100 according to the present embodiment, as shown in FIG. 4, neutrons emitted from the neutron source 104 are made incident from one surface side of the concrete structure 200 and neutrons reflected from the concrete structure 200 are reflected. The distribution of water in the concrete structure 200 is evaluated by detecting the neutrons with a neutron detector 106 provided on the same surface side. For this reason, it becomes possible to evaluate the distribution of water from the same surface side of the concrete structure 200, that is, in one direction. The fast neutrons generated from the neutron source 104 are reflected or scattered by the concrete structure 200 or the junka X or the like to reach the neutron detector 106. At this time, if there is a "water pool" in which water is collected in the Janka X, etc., the neutrons are decelerated by hydrogen in the water, and it becomes easier to be converted into thermal neutrons than in the case where there is no "water pool". The distribution of water in the concrete structure 200 is evaluated by focusing attention on the difference in nuclear reaction. According to the moisture estimation system 100, it is possible to estimate the distribution of moisture from the surface of the concrete structure 200 to a depth about the number average free path of neutrons.

水分推定システム100では、図5に示すように、評価対象であるコンクリート構造物200を検査対象領域202とその周辺領域204に分けて、検査対象領域202内の水分の分布を評価する。このとき、検査対象領域202を複数(m個:mは2以上の整数)の分割対象領域k(kは領域番号:k=1,…,m)に分割して、それぞれの分割対象領域k内の水分を推定することによってコンクリート構造物200内の水分の分布を評価する。   In the water | moisture-content estimation system 100, as shown in FIG. 5, the concrete structure 200 which is evaluation object is divided into the test object area | region 202 and its periphery area | region 204, and distribution of the water | moisture content in the test object area | region 202 is evaluated. At this time, the inspection target area 202 is divided into a plurality of (m: m is an integer of 2 or more) division target areas k (k is an area number: k = 1,..., M). The distribution of water in the concrete structure 200 is evaluated by estimating the water in the container.

図5の例では、コンクリート構造物200の一つの表面(X−Y平面)を測定面として、測定面の検査対象領域202を6つの領域にメッシュ分割し、さらに深さ方向(Z方向)に2層に分割している。すなわち、検査対象領域202の測定面から遠い層は分割対象領域k(1≦k≦6)に分割され、測定面から近い層は分割対象領域k(7≦k≦12)に分割されている。もちろん、検査対象領域202の分割数は、これに限定されるものではなく、表面内の分割数は6以外でもよく、深さ方向の層数も2以外でもよい。   In the example of FIG. 5, with one surface (X-Y plane) of the concrete structure 200 as a measurement surface, the inspection target region 202 of the measurement surface is divided into six regions and meshed further in the depth direction (Z direction) It is divided into two layers. That is, the layer far from the measurement surface of the inspection target region 202 is divided into division target regions k (1 ≦ k ≦ 6), and the layer close to the measurement surface is divided into division target regions k (7 ≦ k ≦ 12) . Of course, the number of divisions of the inspection target area 202 is not limited to this, the number of divisions in the surface may be other than six, and the number of layers in the depth direction may also be other than two.

水分推定システム100では、コンクリート構造物200の測定面側の分割対象領域kに中性子線源104と中性子検出器106を配置して応答を測定し、実測された応答に対応する水分の分布を水分推定装置102により解析する。すなわち、中性子線源104を配置した分割対象領域kの位置iと中性子検出器106を配置した分割対象領域kの位置jとの組み合わせ毎に実際の応答Di,jと解析された応答Ri,jとが合致する水分の分布を求める。 In the water estimation system 100, the neutron beam source 104 and the neutron detector 106 are arranged in the division target region k on the measurement surface side of the concrete structure 200, the response is measured, and the distribution of water corresponding to the measured response is Analysis is performed by the estimation device 102. That is, for each combination of the position i of the division target region k where the neutron beam source 104 is arranged and the position j of the division target region k where the neutron detector 106 is arranged , the response R i analyzed as the actual response D i, j , J is determined to match the distribution of water.

例えば、図5の検査対象領域202において、領域番号7の分割対象領域に中性子線源104を配置すると共に、領域番号10,11,12の分割対象領域にそれぞれに中性子検出器106を順に配置した場合の実際の応答D7,10,D7,11,D7,12を測定する。また、領域番号8の分割対象領域に中性子線源104を配置すると共に、領域番号10,11,12の分割対象領域にそれぞれに中性子検出器106を順に配置した場合の実際の応答D8,10,D8,11,D8,12を測定する。また、領域番号9の分割対象領域に中性子線源104を配置すると共に、領域番号10,11,12の分割対象領域にそれぞれに中性子検出器106を順に配置した場合の実際の応答D9,10,D9,11,D9,12を測定する。さらに、中性子線源104の位置と中性子検出器106の位置を入れ替えて実際の応答を測定する。すなわち、領域番号10の分割対象領域に中性子線源104を配置すると共に、領域番号7,8,9の分割対象領域にそれぞれに中性子検出器106を順に配置した場合の実際の応答D10,7,D10,8,D10,9を測定する。また、領域番号11の分割対象領域に中性子線源104を配置すると共に、領域番号7,8,9の分割対象領域にそれぞれに中性子検出器106を順に配置した場合の実際の応答D11,7,D11,8,D11,9を測定する。また、領域番号12の分割対象領域に中性子線源104を配置すると共に、領域番号7,8,9の分割対象領域にそれぞれに中性子検出器106を順に配置した場合の実際の応答D12,7,D12,8,D12,9を測定する。 For example, in the inspection target area 202 of FIG. 5, the neutron beam source 104 is disposed in the division target area of area number 7, and the neutron detectors 106 are sequentially disposed in the division target areas of area numbers 10, 11 and 12, respectively. The actual response D 7,10 , D 7 , 11 , D 7 , 12 of the case is measured. The actual response D 8, 10 in the case where the neutron beam source 104 is disposed in the division target area of the area number 8 and the neutron detectors 106 are sequentially disposed in the division target areas of the area numbers 10, 11, 12 respectively. , D 8 , 11 , D 8 , 12 are measured. In addition, actual response D 9, 10 when neutron beam source 104 is arranged in the division target area of area number 9, and neutron detectors 106 are arranged in order in the division target areas of area numbers 10, 11 and 12, respectively. , D 9 , 11 and D 9 , 12 are measured. Furthermore, the position of the neutron beam source 104 and the position of the neutron detector 106 are interchanged to measure the actual response. That is, actual response D 10, 7 when neutron beam source 104 is arranged in the division target area of area number 10 and neutron detectors 106 are arranged in order in the division target areas of area numbers 7, 8, 9 respectively , D 10, 8 and D 10 , 9 are measured. The actual response D 11, 7 in the case where the neutron beam source 104 is disposed in the division target area of the area number 11 and the neutron detectors 106 are sequentially disposed in the division target areas of the area numbers 7, 8, and 9 respectively. , D 11 , 8 and D 11 , 9 are measured. The actual response D 12,7 when the neutron beam source 104 is arranged in the division target area of the area number 12 and the neutron detectors 106 are sequentially arranged in the division target areas of the area numbers 7, 8, 9 respectively , D 12 , 8 and D 12 , 9 are measured.

なお、測定時には、測定データの統計誤差が十分に小さくなるような測定時間で測定を行うことが好適である。   At the time of measurement, it is preferable to perform the measurement in such a measurement time that the statistical error of the measurement data is sufficiently reduced.

一方、水分推定装置102では、放射線輸送計算コードに基づく解析を行う。水分推定装置102は、コンクリート構造物200の分割対象領域kについて3次元形状を扱える放射線輸送計算コード、例えば汎用コードであるモンテカルロ法に基づくMCNPやSN法に基づくTORT等を使用して解析する。   On the other hand, the water estimation device 102 performs analysis based on the radiation transport calculation code. The moisture estimation device 102 analyzes the division target region k of the concrete structure 200 using a radiation transport calculation code that can handle a three-dimensional shape, for example, a general-purpose code such as MCNP based on Monte Carlo method or TERT based on SN method.

例えば、中性子線源104の位置iから毎秒1個の中性子が発生するものとし、中性子検出器106の位置jで得られる応答(計数率や線量率)は、分割対象領域kの各領域からの寄与を領域毎に評価する必要がある。また、周辺領域204からの影響も考慮する必要がある。   For example, it is assumed that one neutron is generated per second from the position i of the neutron beam source 104, and the response (counting rate and dose rate) obtained at the position j of the neutron detector 106 Contributions need to be assessed region by region. In addition, the influence from the peripheral area 204 also needs to be considered.

すなわち、中性子線源104の位置i及び中性子検出器106の位置jとした場合の応答Ri,jは、数式(3)により表わされる。

Figure 0006518939
ここで、ri,j,kは中性子線源104の位置i及び中性子検出器106の位置jとした場合の分割対象領域kからの寄与、Ari,jは周辺領域204からの寄与を示す。 That is, the response R i, j in the case where the position i of the neutron beam source 104 and the position j of the neutron detector 106 are represented by Formula (3).
Figure 0006518939
Here, r i, j, k represents the contribution from the division target area k when the position i of the neutron beam source 104 and the position j of the neutron detector 106 are taken, and Ar i, j represents the contribution from the peripheral area 204 .

コンクリート構造物200に入射した中性子の挙動は、図6に示すように、中性子線源104から反射点Aに至るまでの高速中性子の減衰と、反射点Aにおける散乱を非衝突の透過と衝突による散乱とによって表現される。   The behavior of neutrons incident on the concrete structure 200 is, as shown in FIG. 6, the attenuation of fast neutrons from the neutron source 104 to the reflection point A and the scattering at the reflection point A by noncollision transmission and collision. Expressed by scattering.

反射点Aにおける中性子線の非衝突線束Φ(E)は、数式(4)によって表現できる。ここで、Eは中性子の初期エネルギー、Sは中性子の初期強度、P1は中性子線源104から反射点Aまでの距離、σ(E)はエネルギーEの中性子の全断面積である。

Figure 0006518939
The non-collision flux Φ A (E) of the neutron beam at the reflection point A can be expressed by Equation (4). Here, E is the initial energy of neutrons, S is the initial intensity of neutrons, P1 is the distance from neutron source 104 to reflection point A, σ T (E) is the total cross section of neutrons of energy E.
Figure 0006518939

反射点Aにおける初回散乱線源FS(E’)は、数式(5)によって表現できる。ここで、σ(E,E’,θ1,θ2)は入射エネルギーEから反射エネルギーE’にエネルギー変化し、入射角θ1から反射角θ2に角度変化するときの微分散乱断面積である。

Figure 0006518939
The first-time scattered radiation source FS (E ′) at the reflection point A can be expressed by Equation (5). Here, σ S (E, E ′, θ1, θ2) is a differential scattering cross section when energy changes from incident energy E to reflected energy E ′ and angle changes from incident angle θ1 to reflection angle θ2.
Figure 0006518939

中性子検出器106における線束Φ(E’)は、数式(6)によって表現できる。ここで、B(E’)はエネルギーE’を有する中性子が反射点Aから中性子検出器106に至る経路における散乱成分の寄与を示すビルドアップ係数、P2は反射点Aから中性子検出器106までの距離、σ(E’)はエネルギーE’の中性子の全断面積である。

Figure 0006518939
The flux Φ B (E ′) in the neutron detector 106 can be expressed by equation (6). Here, B (E ′) is a buildup coefficient indicating the contribution of the scattered component in the path from the reflection point A to the neutron detector 106 by neutrons having the energy E ′, and P 2 is from the reflection point A to the neutron detector 106 The distance, σ T (E ′), is the total cross section of the neutron of energy E ′.
Figure 0006518939

中性子線源104の位置iからエネルギーEの中性子が放射され、各分割対象領域kを反射点Aとして散乱されて中性子検出器106の位置jにてエネルギーE’の中性子として検出される中性子の寄与ri,j,kは、数式(6)の線束Φ(E’)において中性子の初期強度Sを1にして算出することができる。すなわち、コンクリート構造物200の形状及び寸法と中性子線源104の位置i及び中性子検出器106の位置jとの関係から中性子線源104から反射点A(各分割対象領域k)までの距離P1、反射点Aから中性子検出器106までの距離P2、反射点Aに対する入射角θ1及び反射角θ2は幾何学的に定められるので、これらの値を数式(6)に代入することにより線束Φ(E’)を中性子の寄与ri,j,kとして求めることができる。 Neutrons of energy E are emitted from the position i of the neutron source 104, and each divided target area k is scattered as a reflection point A, and contribution of neutrons detected as neutrons of energy E 'at the position j of the neutron detector 106 r i, j, k can be calculated by setting the initial intensity S of neutrons to 1 in the flux Φ B (E ′) of equation (6). That is, from the relationship between the shape and size of the concrete structure 200 and the position i of the neutron beam source 104 and the position j of the neutron detector 106, the distance P1 from the neutron beam source 104 to the reflection point A (each divided target area k) Since the distance P2 from the reflection point A to the neutron detector 106 and the incident angle θ1 and the reflection angle θ2 with respect to the reflection point A are geometrically determined, the flux Φ B () can be obtained by substituting these values into equation (6). E ′) can be determined as the neutron contribution r i, j, k .

ここで、数式(4)及び数式(6)における全断面積σ(E),σ(E’)はコンクリート構造物200を構成するコンクリートの元素密度組成に基づいて定められる。また、反射点Aにおける微分散乱断面積σ(E,E’,θ1,θ2)は、コンクリート構造物200の元素密度組成に加えてジャンカ等に溜まった水分の原子数密度に基づいて定められる。近似的には、コンクリート構造物200を構成するコンクリートの微分散乱断面積と水の微分散乱断面積をそれぞれの反射点Aにおける体積比(体積率)から定めることができる。すなわち、コンクリート構造物200を構成するコンクリート以外に相当量の水分が含まれている対象領域の微分散乱断面積は、コンクリート構造物200を構成するコンクリートの微分散乱断面積と水の微分散乱断面積にそれぞれの体積比を掛けて加算すればよい。 Here, the total cross-sectional areas σ T (E) and σ T (E ′) in Equation (4) and Equation (6) are determined based on the element density composition of the concrete that constitutes the concrete structure 200. In addition to the element density composition of the concrete structure 200, the differential scattering cross section σ S (E, E ′, θ1, θ2) at the reflection point A is determined based on the atomic number density of water collected in a junka or the like. . Approximately, the differential scattering cross section of the concrete constituting the concrete structure 200 and the differential scattering cross section of water can be determined from the volume ratio (volume ratio) at each reflection point A. That is, the differential scattering cross-sectional area of the target region containing a considerable amount of water other than the concrete constituting the concrete structure 200 is the differential scattering cross-sectional area of the concrete constituting the concrete structure 200 and the differential scattering cross-sectional area of water Can be multiplied by each volume ratio and added.

そこで、反射点Aがコンクリート構造物200を構成するコンクリートである(コンクリートの体積比100%)として、コンクリート構造物200の元素密度組成に基づいて定められた反射点Aにおける微分散乱断面積σSC(E,E’,θ1,θ2)を適用したときの寄与ri,j,kを寄与rci,j,kとして求める。また、反射点Aが水である(水の体積比100%)として、水の反射点Aにおける微分散乱断面積σSW(E,E’,θ1,θ2)を適用したときの寄与ri,j,kを寄与rwi,j,kとして求める。そして、分割対象領域kに含まれる水の体積比ωを水分量とすると、中性子検出器106で計測される応答は、数式(7)で表わされる。

Figure 0006518939
Therefore, the differential scattering cross section σ SC at the reflection point A determined based on the element density composition of the concrete structure 200 as the reflection point A being concrete that constitutes the concrete structure 200 (100% by volume of concrete) The contribution r i, j, k when (E, E ′, θ 1, θ 2) is applied is obtained as the contribution r c i, j, k . In addition, contribution r i, when the differential scattering cross section σ SW (E, E ′, θ1, θ2) at the reflection point A of water is applied as the reflection point A is water (100% volume ratio of water) Find j, k as contributions r w i, j, k . Then, assuming that the volume ratio ω k of water contained in the division target region k is the water content, the response measured by the neutron detector 106 is expressed by Equation (7).
Figure 0006518939

数式(7)を展開すると数式(8)が得られる。

Figure 0006518939
Equation (8) is obtained by expanding Equation (7).
Figure 0006518939

評価対象であるコンクリート構造物200に中性子線源104及び中性子検出器106をそれぞれ位置i及び位置jに順次移動させて実測した応答Di,jを数式(8)に代入することにより、未知数を水の体積比ωとする連立方程式が得られる。したがって、この連立方程式を解くことにより分割対象領域kにおける水の体積比ω、すなわちコンクリート構造物200の検査対象領域202を分割した分割対象領域kの各々における水分の分布を求めることができる。 By substituting the measured response Di, j into the equation (8) by sequentially moving the neutron beam source 104 and the neutron detector 106 to the position i and the position j respectively in the concrete structure 200 to be evaluated, A simultaneous equation with the volume ratio of water ω k is obtained. Therefore, the volume ratio omega k of the water in the divided region of interest k by solving the simultaneous equations, that is, to determine the distribution of the water in each of the divided target area k obtained by dividing the inspection area 202 of the concrete structure 200.

しかしながら、寄与rci,j,k、寄与rwi,j,k及び周辺領域204からの寄与Ari,jを算出するミュレーションには計算誤差などが含まれ、実測された応答Di,jにも統計誤差等の不確かさが含まれる。したがって、数式(8)は厳密には成立しない。そこで、数式(8)の誤差が最小になるような水の体積比ωを求めれば、コンクリート構造物200の分割対象領域kの水分量を推定することができる。 However, the contribution rc i, j, k, contributing rw i, j, the contribution Ar i from k and the peripheral region 204, the simulation of calculating the j includes such calculation error, measured responses D i, j Also include uncertainty such as statistical error. Therefore, Equation (8) does not hold strictly. Then, if the volume ratio ω k of water such that the error of Expression (8) is minimized is obtained, the water content of the division target area k of the concrete structure 200 can be estimated.

例えば、最小二乗法を適用してコンクリート構造物200内の水分の分布を推定することができる。数式(8)について左辺の定数項Ci,jと右辺の計算値の残差Rを最小二乗法に基づく目的関数として数式(9)のように定義する。

Figure 0006518939
For example, the least squares method can be applied to estimate the distribution of water in the concrete structure 200. The residual R of the constant term C i, j on the left side and the calculated value on the right side of Expression (8) is defined as Expression (9) as an objective function based on the least squares method.
Figure 0006518939

この目的関数を最小にするために最急降下法を適用する。制約条件として水の体積比ωが0以上であるという条件を設定する。一般に数式(9)で与えられる問題を閉じた形の演算で解くことは不可能であり、何らかの繰り返し計算を行う必要がある。そこで、数式(10)となり、且つこの数列が目的関数の最小値に収束するようなベクトル列ω (0), ω (1),・・・,ω (n)を求めるアルゴリズムを用いる。

Figure 0006518939
The steepest descent method is applied to minimize this objective function. As a constraint condition, a condition that water volume ratio ω k is 0 or more is set. In general, it is impossible to solve the problem given by Equation (9) in a closed-form operation, and it is necessary to perform some repetitive calculations. Then, using an algorithm for obtaining a vector sequence ω k (0) , ω k (1) ,..., Ω k (n) such that equation (10) is obtained and this number sequence converges to the minimum value of the objective function .
Figure 0006518939

このアルゴリズムは一般に次の三つのステップで構成される。
(1)適当な水の体積比ωの初期点ω (0)を与え、繰り返し回数を表すパラメータnを0とする。
(2)現在のパラメータnに対して得られた水の体積比ω (n)が数式(9)の解であるか否かの判定を行う。水の体積比ω (n)が判定の条件を満たせば、計算を終了して水の体積比ω (n)を解とする。条件を満たさなければ次のステップに進む。
(3)パラメータnに対する計算結果を基にして、目的関数f(ω)を減少させる新しい水の体積比ω (n+1)を求める。パラメータnを1増加させてステップ(2)に戻る。
This algorithm generally consists of the following three steps:
(1) An initial point ω k (0) of a suitable volume ratio of water ω k is given, and a parameter n representing the number of repetitions is 0.
(2) It is determined whether the volume ratio ω k (n) of water obtained with respect to the current parameter n is a solution of equation (9). If the volume ratio of water ω k (n) satisfies the determination condition, the calculation is ended and the volume ratio of water ω k (n) is solved. If the condition is not satisfied, the process proceeds to the next step.
(3) Based on the calculation result for the parameter n, the volume ratio ω k (n + 1) of new water for reducing the objective function f (ω) is obtained . The parameter n is incremented by 1 and the process returns to step (2).

n回目の繰り返し計算で次の水の体積比ω (n)を定めるには、現在の水の体積比ω (n−1)からの探索方向とその方向に沿った歩み幅を決める必要がある。その結果、n+1回目の水の体積比ω (n+1)は数式(11)により得られる。

Figure 0006518939
In order to determine the next water volume ratio ω k (n) in the n-th repeated calculation, it is necessary to determine the search direction from the current water volume ratio ω k (n-1) and the step width along that direction There is. As a result, the volume ratio ω k (n + 1) of the (n + 1) th water is obtained by the equation (11).
Figure 0006518939

ここで、p (n)は方向を表す単位ベクトル、α (n)は歩み幅を表す正のスカラー値である。探索方向p (n)は、数式(12)を用いた最急降下法によって求められる。

Figure 0006518939
Here, p k (n) is a unit vector representing a direction, and α k (n) is a positive scalar value representing a step width. The search direction p k (n) is obtained by the steepest descent method using Expression (12).
Figure 0006518939

このように、最小二乗法を適用して数式(8)の連立方程式の左辺と右辺との差(誤差)を最小とする水の体積比ωを求める。これによって、コンクリート構造物200を分割した分割対象領域kの各々における水の体積比ωを推定することができる。 As described above, the least squares method is applied to determine the volume ratio ω k of water that minimizes the difference (error) between the left side and the right side of the simultaneous equations of Equation (8). By this, it is possible to estimate the volume ratio ω k of water in each of the division target regions k obtained by dividing the concrete structure 200.

水分推定装置102は、上記の方法により求められた水の体積比ωを出力部16に出力する。このとき、調査対象であるコンクリート構造物200の検査対象領域202を分割対象領域kに分割して表示し、分割対象領域kの各々の水の体積比ωを可視化して表示することが好適である。例えば、水の体積比ωを0%〜100%の範囲で色を変化させて表示させる等、ユーザが視覚から水の体積比ωを容易に把握できるような態様で表示させることが好適である。 The water estimation device 102 outputs the volume ratio ω k of water obtained by the above method to the output unit 16. At this time, it is preferable to divide and display the inspection target area 202 of the concrete structure 200 to be investigated into the division target area k and visualize and display the volume ratio ω k of each water of the division target area k. It is. For example, the volume ratio omega k of water or the like to be displayed by changing the color in the range of 0% to 100%, preferably by the user is displayed in a manner easily grasp the volume ratio omega k of water from the visual It is.

具体的には、水分推定装置102は、図7に示すフローチャートに沿ってコンクリート構造物の水分推定を行う。ステップS10では、コンクリート構造物200に検査対象領域202及び周辺領域204を設定すると共に、検査対象領域202を分割して分割対象領域kを設定する。ステップS12では、各分割対象領域kを反射点Aとして、分割対象領域kのコンクリートの体積比が100%である場合のパラメータを用いて中性子線源104の位置iと中性子検出器106の位置jの組み合わせ毎に寄与rci,j,kを算出する。このとき、計算に用いられるパラメータは評価対象であるコンクリート構造物200を構成するコンクリートの元素密度組成に基づいて定める。ステップS14では、各分割対象領域kを反射点Aとして、分割対象領域kの水の体積比が100%である場合のパラメータを用いて中性子線源104の位置iと中性子検出器106の位置jの組み合わせ毎に寄与rwi,j,kを算出する。ステップS16では、中性子線源104の位置iと中性子検出器106の位置jの組み合わせ毎に応答Di,jを実測する。ステップS18では、ステップS12及びS14で解析された寄与rci,j,k及び寄与rwi,j,k並びにステップS16において実測された応答Di,jにより表わされる数式(8)の連立方程式に基づいて最小二乗法により分割対象領域kの各々における水の体積比ωを推定する。ステップS20では、ステップS18において得られた水の体積比ωを可視化して表示する。 Specifically, the moisture estimation device 102 estimates the moisture of the concrete structure according to the flowchart shown in FIG. 7. In step S10, the inspection target area 202 and the peripheral area 204 are set in the concrete structure 200, and the inspection target area 202 is divided to set a division target area k. In step S12, the position i of the neutron beam source 104 and the position j of the neutron detector 106 are determined using parameters in the case where the volume ratio of concrete in the division target area k is 100%, with each division target area k as the reflection point A. The contribution rci , j, k is calculated for each combination of At this time, the parameters used for the calculation are determined based on the element density composition of the concrete constituting the concrete structure 200 to be evaluated. In step S14, the position i of the neutron beam source 104 and the position j of the neutron detector 106 are determined using parameters in the case where the volume ratio of water in the division target area k is 100%, with each division target area k as the reflection point A. The contribution r w i, j, k is calculated for each combination of In step S16, the response D i, j is measured for each combination of the position i of the neutron beam source 104 and the position j of the neutron detector 106. In step S18, the simultaneous equations of equation (8) represented by the contribution rc i, j, k and the contribution r w i, j, k analyzed in steps S12 and S14 and the response D i, j measured in step S16 are used. Based on the above, the volume ratio ω k of water in each of the division target regions k is estimated by the least squares method. In step S20, the volume ratio ω k of water obtained in step S18 is visualized and displayed.

なお、コンクリート構造物200のコンクリートの種別(コンクリートの元素密度組成)、形状及び大きさ、さらに予めコンクリート構造物200の検査対象領域202を分割した分割対象領域kの水の体積比ωを設定したときの応答Di,jを解析してデータベース化しておくことも好適である。データベースは、水分推定装置102の記憶部12に記憶させておけばよい。これにより、実際に中性子線源104及び中性子検出器106を設置して実測された応答Di,jとデータベースに登録されている解析された応答Di,jを比較することにより、コンクリート構造物200内の水分の分布を推定することができる。 In addition, the type (concrete elemental density composition), shape and size of concrete structure 200 and the volume ratio ω k of water of division target area k obtained by dividing inspection target area 202 of concrete structure 200 in advance are set. It is also preferable to analyze the response D i, j at the time of making a database. The database may be stored in the storage unit 12 of the water estimation device 102. Thus, the concrete structure can be obtained by comparing the actually measured response Di, j with the installed neutron beam source 104 and neutron detector 106 and the analyzed response Di, j registered in the database. The distribution of water within 200 can be estimated.

データベースの構築に当たっては、検査対象のコンクリート構造物200の元素密度組成や結合水の含有量を予め調査してデータベース化すること、複数のコンクリート構造物の元素密度組成や結合水の含有量を予め調査してデータベース化すること、またコンクリート構造物の設計書から元素密度組成や結合水の含有量を予測してデータベース化すること等が可能である。   In constructing the database, the element density composition of the concrete structure 200 to be inspected and the content of the bound water are previously investigated and made into a database, and the element density composition of the plurality of concrete structures and the content of the bound water are previously It is possible to investigate and create a database, and to predict and store elemental density composition and content of bound water from a design document of a concrete structure.

また、上記の記載から明らかなように本実施の形態における水分推定方法においては、特にトレーサを用いずにコンクリート構造物200内の水分の分布を推定することができる。このためトレーサが添加されていない既存のコンクリート構造物に適用することが可能である。また、本実施の形態においては、中性子線源104及び中性子検出器106は、調査対象であるコンクリート構造物200の同一表面側に配置されている例を示したが、コンクリート構造物の一方向からの水分の存在の推定が可能な範囲において、同一表面側でなくてもよい。なお、コンクリート構造物200の一方向とは、コンクリート構造物200の周囲の空間を2つ以上に分けた場合の1つの側からをいう。   Further, as is clear from the above description, in the water estimation method in the present embodiment, the distribution of water in the concrete structure 200 can be estimated without using a tracer. For this reason, it is possible to apply to the existing concrete structure to which the tracer is not added. In the present embodiment, the neutron beam source 104 and the neutron detector 106 are arranged on the same surface side of the concrete structure 200 to be investigated, but from one direction of the concrete structure It is not necessary to be on the same surface side as far as it is possible to estimate the presence of water. In addition, the one direction of the concrete structure 200 means from one side when the space around the concrete structure 200 is divided into two or more.

10 処理部、12 記憶部、14 入力部、16 出力部、18 インターフェース部、20 中性子発生源、22,24,26 遮蔽壁、100 水分推定システム、102 水分推定装置、104 中性子線源、106 中性子検出器、200 コンクリート構造物、202 検査対象領域、204 周辺領域。   Reference Signs List 10 processing unit, 12 storage unit, 14 input unit, 16 output unit, 18 interface unit, 20 neutron source, 22, 24, 26 shielding wall, 100 moisture estimation system, 102 moisture estimation device, 104 neutron beam source, 106 neutron Detector, 200 concrete structure, 202 inspection area, 204 peripheral area.

Claims (10)

コンクリートの元素密度組成及び結合水の含有量を条件として、中性子の線源の位置と計測点の位置との関係に基づいて前記線源から放射された中性子が前記コンクリートにより構成された調査対象であるコンクリート構造物を介して前記計測点において計測されるときの応答をモデル化してモンテカルロ法及びSN法を含む放射線輸送計算コードにより解析する第1のステップと、
前記線源から放射された中性子が前記調査対象であるコンクリート構造物を介して前記計測点において計測されるときの応答を実測する第2のステップと、
前記実測された応答と前記解析された応答とを用いて前記調査対象であるコンクリート構造物内の水分の存在を推定する第3のステップと
を備えることを特徴とする中性子を用いたコンクリート構造物の水分推定方法。
Under the condition of the element density composition of the concrete and the content of bound water, the neutron emitted from the radiation source based on the relationship between the position of the neutron radiation source and the position of the measurement point is the survey target composed of the concrete A first step of modeling a response measured at the measurement point through a concrete structure and analyzing it by a radiation transport calculation code including a Monte Carlo method and an SN method;
A second step of measuring a response when neutrons emitted from the radiation source are measured at the measurement point through the concrete structure to be investigated;
A concrete structure using neutrons, comprising: a third step of estimating the presence of water in the concrete structure to be investigated using the measured response and the analyzed response. Moisture estimation method.
請求項1に記載の中性子を用いたコンクリート構造物の水分推定方法であって、
前記第3のステップは、前記第1のステップにおいて複数の前記元素密度組成及び前記結合水の含有率の組み合わせであるコンクリートの種別毎に解析した前記解析された応答を記憶させたデータベースを参照して、前記実測された応答と前記調査対象であるコンクリート構造物を構成するコンクリートの種別に対して前記解析された応答とを用いて前記調査対象であるコンクリート構造物内の水分の存在を推定することを特徴とする中性子を用いたコンクリート構造物の水分推定方法。
A method of estimating water content of a concrete structure using neutrons according to claim 1, wherein
The third step refers to a database storing the analyzed response analyzed for each type of concrete, which is a combination of a plurality of the element density compositions and the content of the binding water in the first step. Estimate the presence of moisture in the concrete structure to be investigated using the measured response and the analyzed response to the type of concrete constituting the concrete structure to be investigated Method of estimating moisture of concrete structure using neutrons characterized by
請求項1又は2に記載の中性子を用いたコンクリート構造物の水分推定方法であって、
前記第1のステップは、
前記調査対象であるコンクリート構造物にそれぞれ検査対象領域及び周辺領域を設定してモデル化し、
前記検査対象領域をm個の分割対象領域k(mは2以上の整数、k=1,…,m)にメッシュ分割して、複数の前記線源の位置i(iは2以上の整数)及び複数の前記計測点の位置j(jは2以上の整数)の組み合わせの各々についての前記分割対象領域kの各々からの寄与ri,j,kを求め、
前記線源の位置iから放射された中性子が前記調査対象であるコンクリート構造物を介して前記計測点の位置jにおいて計測されるときの前記解析された応答Ri,j
Figure 0006518939
(ただし、Ari,jは、前記線源の位置iから放射された中性子が前記調査対象であるコンクリート構造物を介して前記計測点の位置jにおいて計測されるときの前記周辺領域からの応答の寄与)
により求めることを特徴とするコンクリート構造物の水分推定方法。
A method of estimating water content of a concrete structure using neutrons according to claim 1 or 2,
The first step is
Each of the concrete structures to be surveyed is modeled by setting the inspection target area and the surrounding area,
The inspection target area is divided into m division target areas k (m is an integer of 2 or more, k = 1,..., M), and the positions i of the plurality of radiation sources (i is an integer of 2 or more) And a contribution ri , j, k from each of the division target regions k for each combination of positions j (j is an integer of 2 or more) of the plurality of measurement points,
The analyzed response R i, j when neutrons emitted from the position i of the radiation source are measured at the position j of the measurement point through the concrete structure to be investigated
Figure 0006518939
(Where Ar i, j is the response from the surrounding area when neutrons emitted from the position i of the radiation source are measured at the position j of the measurement point through the concrete structure to be investigated) Contribution of
The moisture estimation method of the concrete structure characterized by calculating by
請求項3に記載のコンクリート構造物の水分推定方法であって、
前記第1のステップは、
前記線源の位置i及び前記計測点の位置jの組み合わせの各々について、前記分割対象領域kに前記コンクリートが体積率100%で存在するときの寄与rci,j,k及び、前記分割対象領域kに前記水が体積率100%で存在するときの寄与rwi,j,kを求め、
前記第2のステップは、
前記線源の位置i及び前記計測点の位置jの組み合わせの各々について、前記線源から放射された中性子が前記調査対象であるコンクリート構造物を介して前記計測点において計測されるときの応答Di,jを実測し、
前記第3のステップは、
Figure 0006518939
(ただし、Sは、前記線源の強度)
の関係を用いて前記分割対象領域kにおける水分の体積率ωを推定することを特徴とするコンクリート構造物の水分推定方法。
The method of estimating moisture of a concrete structure according to claim 3, wherein
The first step is
For each combination of position i of the radiation source and position j of the measurement point, contributions rc i, j, k when the concrete is present at a volume ratio of 100% in the division target area k, and the division target area The contribution r w i, j, k when the water exists at a volume ratio of 100% is obtained in k ,
The second step is
For each combination of position i of the source and position j of the measurement point, the response D when neutrons emitted from the source are measured at the measurement point via the concrete structure to be investigated Measure i, j ,
The third step is
Figure 0006518939
(Where S is the intensity of the radiation source)
A method of estimating water content in a concrete structure, comprising: estimating a volume ratio ω k of water in the division target area k using the relationship of
請求項4に記載のコンクリート構造物の水分推定方法であって、
前記数式(2)の誤差が最小となる前記分割対象領域kにおける水分の体積率ωを推定することを特徴とするコンクリート構造物の水分推定方法。
A method of estimating moisture of a concrete structure according to claim 4, wherein
A moisture estimation method for a concrete structure, comprising: estimating a volume fraction ω k of moisture in the division target region k which minimizes the error of the equation (2).
請求項1〜5のいずれか1項に記載のコンクリート構造物の水分推定方法であって、
前記第3のステップにおいて推定された前記調査対象であるコンクリート構造物内の水分の存在を可視化する第4のステップを備えることを特徴とするコンクリート構造物の水分推定方法。
It is the moisture estimation method of the concrete structure of any one of Claims 1-5, Comprising:
A method of estimating moisture of a concrete structure, comprising the fourth step of visualizing the presence of moisture in the concrete structure to be surveyed estimated in the third step.
請求項1〜6のいずれか1項に記載のコンクリート構造物の水分推定方法であって、
前記第2のステップにおいて使用される前記線源は、前記調査対象であるコンクリート構造物に対面する領域以外を中性子遮蔽手段で覆った構造を有することを特徴とするコンクリート構造物の水分推定方法。
It is a moisture estimation method of the concrete structure according to any one of claims 1 to 6,
The method of estimating moisture of a concrete structure according to claim 1, wherein the radiation source used in the second step has a structure in which a neutron shielding means covers a region other than a region facing the concrete structure to be investigated.
請求項1〜7のいずれか1項に記載のコンクリート構造物の水分推定方法であって、
前記第2のステップは、前記調査対象であるコンクリート構造物の表面上を前記線源及び前記計測点を順次移動させて前記実測された応答を実測することを特徴とするコンクリート構造物の水分推定方法。
It is a moisture estimation method of the concrete structure according to any one of claims 1 to 7,
In the second step, the water source and the measurement point are sequentially moved on the surface of the concrete structure to be surveyed to measure the measured response, and the water structure of the concrete structure is estimated. Method.
コンクリートの元素密度組成及び結合水の含有量を条件として、中性子の線源の位置と計測点の位置との関係に基づいて前記線源から放射された中性子が前記コンクリートにより構成された調査対象であるコンクリート構造物を介して前記計測点において計測されるときの応答をモデル化して解析した結果をデータベースとして記憶する記憶手段と、
前記線源から放射された中性子が調査対象であるコンクリート構造物を介して前記計測点において計測されるときの応答を実測する実応答計測手段と、
前記実測された応答と前記解析された応答とを用いて前記調査対象であるコンクリート構造物内の水分の存在を推定する水分推定手段と
を備えることを特徴とする中性子を用いたコンクリート構造物の水分推定システム。
Under the condition of the element density composition of the concrete and the content of bound water, the neutron emitted from the radiation source based on the relationship between the position of the neutron radiation source and the position of the measurement point is the survey target composed of the concrete Storage means for storing, as a database, the results of modeling and analyzing responses when measured at the measurement points via a concrete structure;
An actual response measurement means for measuring a response when neutrons emitted from the radiation source are measured at the measurement point via a concrete structure to be investigated;
A concrete structure using neutrons, comprising: water estimation means for estimating the presence of water in the concrete structure to be investigated using the measured response and the analyzed response. Moisture estimation system.
請求項9に記載のコンクリート構造物の水分推定システムであって、
前記水分推定手段において推定された前記調査対象であるコンクリート構造物内の水分の存在を可視化して表示する表示手段をさらに備えることを特徴とするコンクリート構造物の水分推定システム。
It is a moisture estimation system of the concrete structure according to claim 9, which is
A moisture estimation system for a concrete structure, further comprising display means for visualizing and displaying the presence of moisture in the concrete structure to be surveyed estimated by the moisture estimation means.
JP2015123727A 2015-06-19 2015-06-19 Method of estimating moisture of concrete structure and system of estimating moisture of concrete structure Active JP6518939B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2015123727A JP6518939B2 (en) 2015-06-19 2015-06-19 Method of estimating moisture of concrete structure and system of estimating moisture of concrete structure

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2015123727A JP6518939B2 (en) 2015-06-19 2015-06-19 Method of estimating moisture of concrete structure and system of estimating moisture of concrete structure

Publications (2)

Publication Number Publication Date
JP2017009371A JP2017009371A (en) 2017-01-12
JP6518939B2 true JP6518939B2 (en) 2019-05-29

Family

ID=57761348

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2015123727A Active JP6518939B2 (en) 2015-06-19 2015-06-19 Method of estimating moisture of concrete structure and system of estimating moisture of concrete structure

Country Status (1)

Country Link
JP (1) JP6518939B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112114349B (en) * 2019-06-21 2023-07-04 清华大学 Radiation analysis method, apparatus and computer readable storage medium

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5089094A (en) * 1973-12-07 1975-07-17
US5083029A (en) * 1991-02-04 1992-01-21 Halliburton Company Measuring water content by neutron thermalization
US5258622A (en) * 1992-03-13 1993-11-02 Troxler Electronic Laboratories, Inc. Apparatus and method for measuring physical characteristics of materials
JPH1048157A (en) * 1996-08-08 1998-02-20 Toray Ind Inc Apparatus for measuring and analyzing with simulation of molecule and method for analyzing chemical structure of substance
JPH1151799A (en) * 1997-08-05 1999-02-26 Haruo Chisaka Detector for leaky place
JP4951869B2 (en) * 2004-08-11 2012-06-13 富士通株式会社 Sample analysis method and analysis system
CN103336018B (en) * 2013-07-02 2016-05-11 李霞波 Single pole low-activity nuclear moisture Density Measuring Instrument

Also Published As

Publication number Publication date
JP2017009371A (en) 2017-01-12

Similar Documents

Publication Publication Date Title
Priyada et al. Intercomparison of gamma scattering, gammatography, and radiography techniques for mild steel nonuniform corrosion detection
JP2014512531A (en) Improved methods and systems for investigating radiation sources in specific areas
Haith et al. Radiographic modelling for NDE of subsea pipelines
US7022998B2 (en) Radiation dose rate measurement
Carasco Coupling gamma ray spectrometry and tomography in a Bayesian frame
Koshti X-ray ray tracing simulation and flaw parameters for crack detection
JP6518939B2 (en) Method of estimating moisture of concrete structure and system of estimating moisture of concrete structure
Scannavino Jr et al. A graphical tool for an analytical approach of scattering photons by the Compton effect
Hykes et al. Radiation source mapping with Bayesian inverse methods
JP2014228362A (en) Nuclear fuel subcriticality measuring and monitoring system and method
JP2018189516A (en) Device, method, and program for calculating surface contamination density distribution
KR102394944B1 (en) Apparatus for measuring and amending activity concentration density of radioactive material and operation calibration methode thereof
Szalóki et al. Reverse Monte Carlo iterative algorithm for quantification of X‐ray fluorescence analysis based on MCNP6 simulation code
Miller et al. Combining radiography and passive measurements for radiological threat localization in cargo
Rimlinger et al. Qualification of gamma spectrometry measurement for the radiological characterization of mixed VLLW cables in particle accelerators
JP6814049B2 (en) Subcriticality measurement method for nuclear fuel
Boin Developments towards the tomographic imaging of local crystallographic structures
Inanc Analysis of X-ray and gamma ray scattering through computational experiments
Bomsdorf et al. Quantitative simulation of coherent X-ray scatter measurements on bulk objects
Al-Bahi et al. Concept design of a PGNAA system for optimizing the performance of gravity separators
JP5926362B1 (en) Radioactivity concentration measuring apparatus and radioactivity concentration measuring method
WO2021230201A1 (en) Concentration detection device and concentration detection method
JP7451306B2 (en) Non-destructive structural analysis equipment, non-destructive structural inspection equipment and non-destructive structural analysis methods
JP7242390B2 (en) Surface contamination density distribution calculation device, method and program
Rebei et al. Quantitative and three-dimensional assessment of holdup material

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20180605

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20190206

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20190325

R150 Certificate of patent or registration of utility model

Ref document number: 6518939

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250