JP2017009371A - Moisture estimation method of concrete structure and moisture estimation system of concrete structure - Google Patents

Moisture estimation method of concrete structure and moisture estimation system of concrete structure Download PDF

Info

Publication number
JP2017009371A
JP2017009371A JP2015123727A JP2015123727A JP2017009371A JP 2017009371 A JP2017009371 A JP 2017009371A JP 2015123727 A JP2015123727 A JP 2015123727A JP 2015123727 A JP2015123727 A JP 2015123727A JP 2017009371 A JP2017009371 A JP 2017009371A
Authority
JP
Japan
Prior art keywords
concrete structure
moisture
response
neutron
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.)
Granted
Application number
JP2015123727A
Other languages
Japanese (ja)
Other versions
JP6518939B2 (en
Inventor
光史 浅見
Mitsufumi Asami
光史 浅見
和哉 西村
Kazuya Nishimura
和哉 西村
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

Landscapes

  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

PROBLEM TO BE SOLVED: To estimate a distribution of moisture in a concrete structure which is an estimation object.SOLUTION: When a neutron emitted from a neutron beam source 104 with a content of the chemical element density composition of the concrete and bound water as a condition based on a relation between the position of a neutron beam source 104 and the position of a neutron beam detector 106 is measured in the neutron detector 106 via a concrete structure 200, the response is modeled and analyzed by a radiation transportation calculation code including the Monte Carlo method and the SN method, and when the neutron emitted from the neutron beam source 104 is measured in the neutron detector 106 via the concrete structure 200, the response is actually measured, and the existence of the moisture in the concrete structure 200 is estimated using the actually measured response and the analyzed response.SELECTED DRAWING: Figure 1

Description

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

コンクリート構造物にはジャンカやひび割れ等の欠損が生ずる可能性があり、そこに水が溜まってコンクリート構造物を構成する鉄筋や構造材を腐食したり脆弱化するおそれがある。したがって、コンクリート構造物内の水分の分布を推定できる技術の実現が望まれている。   There is a possibility that a defect such as a jumper or a crack may occur in the concrete structure, and there is a possibility that water accumulates there and corrodes or weakens the reinforcing bars and the structural material constituting the concrete structure. Therefore, realization of a technique capable of estimating the moisture distribution in the concrete structure is desired.

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

また、セメントを硬化材としてホウ素化合物等のトレーサを添加した注入剤をセメントの裏側に注入し、セメント表面側から中性子を照射してトレーサに救出された熱中性子の数を測定してセメントの注入度合を検出する技術が開示されている(特許文献2)。   Also, injecting cement containing a hardener as a hardener and adding tracers such as boron compounds into the back of the cement, irradiating the neutron from the cement surface, measuring the number of thermal neutrons rescued by the tracer, and injecting the 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, there are two types of neutron sources that generate neutrons to be irradiated: a method using a nuclear reactor and a method using an accelerator. In transmission-type neutron beam detection, it is necessary to place an object to be measured between a neutron source and a sensor. When a concrete structure existing mainly outdoors is used as a measurement object, movement of the neutron source or sensor is not possible. It may be impossible or difficult, and depending on the shape and thickness of the measurement object, it may be impossible to place the object on the back side. In addition, the system becomes larger in principle in relation to the arrangement of the neutron source and the detector. Recently, some research institutions have developed mobile accelerators that can be used outdoors.

また、ホウ素化合物等をトレーサとする技術は、トレーサが添加されていない既存のコンクリート構造物には適用することが原理的に不可能である。   In addition, a technique using a boron compound or the like as a tracer cannot be applied in principle to an existing concrete structure to which no tracer is added.

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

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

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

Figure 2017009371
(ただし、Ari,jは、前記線源の位置iから放射された中性子が前記調査対象であるコンクリート構造物を介して前記計測点の位置jにおいて計測されるときの前記周辺領域からの応答の寄与)により求めることが好適である。 In the first step, an inspection target region and a peripheral region are respectively set and modeled on the concrete structure to be investigated, and the inspection target region is divided into m divided target regions k (m is 2 or more). The mesh is divided into integers, k = 1,..., M), and a plurality of the source positions i (i is an integer of 2 or more) and a plurality of measurement point positions j (j is an integer of 2 or more). The contributions r i, j, k from each of the division target regions k for each combination are obtained, and the measurement points are obtained through the concrete structure that is the object of investigation by neutrons emitted from the position i of the radiation source. The analyzed response R i, j when measured at position j of
Figure 2017009371
(However, Ar i, j is a response from the surrounding region when neutrons radiated 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 (contribution).

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

Figure 2017009371
(ただし、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 exists in the division target region k at a volume ratio of 100% , j, k and a contribution rw i, j, k when the water is present in the division target region k at a volume ratio of 100% are obtained, and the second step includes the position i of the radiation source and the measurement For each combination of point positions j, a response D i, j when neutrons radiated from the radiation source are measured at the measurement point via the concrete structure to be investigated is measured, The third step is
Figure 2017009371
(Where S is the intensity of the radiation source)
It is preferable to estimate the water volume ratio ω k in the division target region k using the relationship.

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

また、前記第3のステップにおいて推定された前記調査対象であるコンクリート構造物内の水分の存在を可視化する第4のステップを備えることが好適である。   Moreover, it is preferable to provide the 4th step which visualizes the presence of the water | moisture content in the concrete structure which is the said investigation object estimated in the said 3rd step.

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

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

本発明の請求項9に係る中性子を用いたコンクリート構造物の水分推定システムは、コンクリートの元素密度組成及び結合水の含有量を条件として、中性子の線源の位置と計測点の位置との関係に基づいて前記線源から放射された中性子が前記コンクリートにより構成された調査対象であるコンクリート構造物を介して前記計測点において計測されるときの応答をモデル化して解析した結果をデータベースとして記憶する記憶手段と、前記線源から放射された中性子が調査対象であるコンクリート構造物を介して前記計測点において計測されるときの応答を実測する実応答計測手段と、前記実測された応答と前記解析された応答とを用いて前記調査対象であるコンクリート構造物内の水分の存在を推定する水分推定手段とを備えることを特徴とする。   A concrete structure moisture estimation system using neutrons according to claim 9 of the present invention is a relationship between the position of a neutron source and the position of a measurement point on condition that the element density composition of concrete and the content of bound water are the conditions. The result of modeling and analyzing the response when neutrons radiated from the radiation source are measured at the measurement point through the concrete structure to be investigated constituted by the concrete is stored as a database. Storage means; actual response measurement means for actually measuring a 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 a moisture estimation means for estimating the presence of moisture in the concrete structure that is the object of investigation using the measured response. To.

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

本発明の請求項1に係る中性子を用いたコンクリート構造物の水分推定方法によれば、コンクリートの元素密度組成及び結合水の含有量を条件として、中性子の線源の位置と計測点の位置との関係に基づいて前記線源から放射された中性子が前記コンクリートにより構成された調査対象であるコンクリート構造物を介して前記計測点において計測されるときの応答をモデル化してモンテカルロ法及びSN法を含む放射線輸送計算コードにより解析する第1のステップと、前記線源から放射された中性子が前記調査対象であるコンクリート構造物を介して前記計測点において計測されるときの応答を実測する第2のステップと、前記実測された応答と前記解析された応答とを用いて前記調査対象であるコンクリート構造物内の水分の存在を推定する第3のステップとを備えることによって、前記評価対象であるコンクリート構造物内の水分の分布を精度よく推定することができる。また、コンクリート構造物の一方向からの水分の存在の推定が可能となり、トレーサが不要で既存のコンクリート構造物への適用も可能となる。   According to the method for estimating the moisture content of a concrete structure using neutrons according to claim 1 of the present invention, the position of the neutron source and the position of the measurement point are obtained on condition that the element density composition of concrete and the content of bound water. Based on the relationship, the response when the neutron radiated from the radiation source is measured at the measurement point through the concrete structure to be investigated constituted by the concrete is modeled, and the Monte Carlo method and the SN method are performed. A first step of analyzing by a radiation transport calculation code including a second step of measuring a response when neutrons radiated from the radiation source are measured at the measurement point via 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 is possible to estimate the presence of moisture from one direction of the concrete structure, and no tracer is required, and application to an existing concrete structure is possible.

また、前記第3のステップは、前記第1のステップにおいて複数の前記元素密度組成及び前記結合水の含有率の組み合わせであるコンクリートの種別毎に解析した前記解析された応答を記憶させたデータベースを参照して、前記実測された応答と前記調査対象であるコンクリート構造物を構成するコンクリートの種別に対して前記解析された応答とを用いて前記調査対象であるコンクリート構造物内の水分の存在を推定することによって、予め前記データベースとして登録された前記解析された応答に基づいて前記実測された応答に対応する前記コンクリート構造物内の水分の存在を推定することができる。   In the third step, the database storing the analyzed response analyzed for each type of concrete, which is a combination of the plurality of element density compositions and the combined water content in the first step. Referring to the measured response and the analyzed response to the type of concrete that constitutes the concrete structure that is the investigation target, the presence of moisture in the concrete structure that is the investigation target. By estimating, it is possible to estimate the presence of moisture in the concrete structure corresponding to the actually measured response 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における水分の体積比を求めることができる。 In the first step, an inspection target region and a peripheral region are respectively set and modeled on the concrete structure to be investigated, and the inspection target region is divided into m divided target regions k (m is 2 or more). The mesh is divided into integers, k = 1,..., M), and a plurality of the source positions i (i is an integer of 2 or more) and a plurality of measurement point positions j (j is an integer of 2 or more). The contributions r i, j, k from each of the division target regions k for each combination are obtained, and the measurement points are obtained through the concrete structure that is the object of investigation by neutrons emitted from the position i of the radiation source. By determining the analyzed response R i, j when measured at the position j of the following equation (1), the volume ratio of moisture in the division target region k obtained by dividing the concrete structure can be obtained. .

また、前記第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 exists in the division target region k at a volume ratio of 100% , j, k and a contribution rw i, j, k when the water is present in the division target region k at a volume ratio of 100% are obtained, and the second step includes the position i of the radiation source and the measurement point For each of the combinations of positions j, a response D i, j when neutrons radiated from the radiation source are measured at the measurement point via the concrete structure to be investigated is measured, and the third The step of estimating the volume ratio ω k of the moisture in the division target region k using the relationship of the formula (2), to the contribution rc i, j, k and the contribution rw i, j, k 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 the moisture in the division target region k that minimizes the error of the mathematical formula (2), the above-mentioned simultaneous equations derived from the mathematical formula (2) are not solved, and Based on the contribution rc i, j, k and the contribution rw i, j, k , the volume ratio ω k of water in the division target region k can be estimated.

また、前記第3のステップにおいて推定された前記調査対象であるコンクリート構造物内の水分の存在を可視化する第4のステップを備えることによって、推定された水分の分布をユーザの視覚によって把握し易い態様で表示させることができる。   Moreover, it is easy to grasp | ascertain the estimated moisture distribution by a user's vision by providing the 4th step which visualizes the presence of the water | moisture content in the concrete structure which is the investigation object estimated in the said 3rd step. It can be displayed in a manner.

また、前記第2のステップにおいて使用される前記線源は、前記調査対象であるコンクリート構造物に対面する領域以外を中性子遮蔽手段で覆った構造を有することによって、中性子線源から前記コンクリート構造物内に入射・散乱せずに中性子検出器106に直接到達する中性子を低減することができる。   Moreover, the said radiation source used in the said 2nd step has a structure which covered the area | region other than the area | region which faces the said concrete structure to be investigated with a neutron shielding means, From the neutron radiation source to the said concrete structure It is possible to reduce neutrons that reach the neutron detector 106 without being incident / scattered.

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

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

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

本発明の実施の形態における水分推定システムの構成を示す図である。It is a figure which shows the structure of the moisture estimation system in embodiment of this invention. 本発明の実施の形態における水分推定装置の構成を示す図である。It is a figure which shows the structure of the moisture 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 neutron scattering model 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の同一表面側に配置されている。   A concrete structure moisture estimation system 100 using neutrons according to an embodiment of the present invention includes a moisture estimation device 102, a neutron source 104, and a neutron detector 106 as shown in FIG. 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.

水分推定装置102は、図2に示すように、一般的なコンピュータの基本構成である処理部10、記憶部12、入力部14、出力部16及びインターフェース部18を含んで構成される。   As shown in FIG. 2, the moisture estimation apparatus 102 includes a processing unit 10, a storage unit 12, an input unit 14, an output unit 16, and an interface unit 18 that are basic configurations of a general computer.

処理部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 a calculation unit such as a CPU, and realizes a moisture estimation process in the concrete structure 200 performed in the moisture estimation device 102 by executing a moisture estimation program stored in the storage unit 12. The storage unit 12 includes a storage device such as a semiconductor memory or a hard disk connected 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 includes an input device such as a keyboard and a mouse. The input unit 14 is used to input various parameters used by the moisture estimation device 102. The output unit 16 includes 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 in the moisture estimation device 102. The output unit 16 presents information related to processing in the moisture estimation apparatus 102 to the user. The interface unit 18 includes an analog / digital converter and the like. The interface unit 18 is used to import data indicating the dose of the neutron beam 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 beam source 104 is configured to include a neutron generation source such as a neutron generation tube or an RI beam source such as 252 Cf. The neutron 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 direction other than the direction facing the concrete structure 200 that is the object to be investigated, that is, the area other than the area facing the concrete structure 200. By shielding the neutron generation source 20, when the neutron beam source 104 is installed on the surface of the concrete structure 200 to be investigated, neutrons and γ rays do not enter and scatter into the concrete structure 200, and a neutron detector Direct arrival of neutrons and γ 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 a plurality of types of shielding walls 22, 24, 26. The shielding wall 22 is provided to heat fast neutrons. The shielding wall 22 is made of a material containing paraffin, for example. The shielding wall 24 is provided to absorb the heated neutrons. The shielding wall 24 is made of a material containing cadmium. The shielding wall 26 is provided to absorb γ rays. The shielding wall 26 is made of a material containing lead, for example. The shielding walls 22 to 26 may have a thickness 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 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 detecting the neutron, the neutron detector 106 outputs the detected neutron to the moisture estimation apparatus 102 via the interface unit 18. When receiving the detection signal from the neutron detector 106, the processing unit 10 of the moisture estimation apparatus 102 counts the number of detected neutrons per unit time. Thereby, the dose rate (count rate) of a neutron beam is measured.

本実施の形態における水分推定システム100では、図4に示すように、中性子線源104から放射される中性子をコンクリート構造物200の1つの表面側から入射させ、コンクリート構造物200から反射された中性子を同一の表面側に設けた中性子検出器106で検出することによってコンクリート構造物200内の水分の分布を評価する。このため、コンクリート構造物200の同一表面側、すなわち一方向からの水分の分布の評価が可能となる。中性子線源104から発生した高速中性子は、コンクリート構造物200やジャンカX等で反射又は散乱されて中性子検出器106に到達する。このとき、ジャンカX等に水が溜まった「水溜まり」が存在すると、中性子は水内の水素によって減速され、「水溜まり」が存在しない場合に比べて熱中性子化され易くなる。この核反応の違いに着目することによってコンクリート構造物200内の水分の分布を評価する。水分推定システム100によれば、コンクリート構造物200の表面から中性子の数平均自由行程程度の深さまでの水分の分布を推定することができる。   In the moisture estimation system 100 according to the present embodiment, as shown in FIG. 4, neutrons emitted from the neutron source 104 are incident from one surface side of the concrete structure 200 and reflected from the concrete structure 200. Is detected by a neutron detector 106 provided on the same surface side to evaluate the moisture distribution in the concrete structure 200. For this reason, it becomes possible to evaluate the moisture distribution from the same surface side of the concrete structure 200, that is, from one direction. Fast neutrons generated from the neutron source 104 are reflected or scattered by the concrete structure 200, the jumper X or the like, and reach the neutron detector 106. At this time, if there is a “puddle” in which water is accumulated in the jumper X or the like, neutrons are decelerated by hydrogen in the water, and thermal neutrons are more easily formed than in the case where there is no “puddle”. The distribution of moisture in the concrete structure 200 is evaluated by paying attention to the difference in the nuclear reaction. According to the moisture estimation system 100, it is possible to estimate the moisture distribution from the surface of the concrete structure 200 to a depth of 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 moisture estimation system 100, as shown in FIG. 5, the concrete structure 200 to be evaluated is divided into an inspection target region 202 and its peripheral region 204, and the moisture distribution in the inspection target region 202 is evaluated. At this time, the inspection target region 202 is divided into a plurality of (m: m is an integer of 2 or more) division target regions k (k is a region number: k = 1,..., M). The distribution of moisture in the concrete structure 200 is evaluated by estimating the moisture in the interior.

図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, one surface (XY plane) of the concrete structure 200 is used as a measurement surface, and the inspection target region 202 of the measurement surface is divided into six regions, and 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 near the measurement surface is divided into division target regions k (7 ≦ k ≦ 12). . Of course, the number of divisions of the inspection target region 202 is not limited to this, and the number of divisions in the surface may be other than 6, and the number of layers in the depth direction may be other than 2.

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

例えば、図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 region 202 of FIG. 5, the neutron source 104 is arranged in the division target region of region number 7 and the neutron detectors 106 are sequentially arranged in the division target regions of region numbers 10, 11, and 12, respectively. The actual responses D 7,10 , D 7,11 , D 7,12 of the case are measured. Further, the actual response D 8,10 when the neutron radiation source 104 is arranged in the division target area of the area number 8 and the neutron detectors 106 are arranged in order in the division target areas of the area numbers 10, 11, 12 respectively. , D 8,11 , D 8,12 are measured. In addition, the neutron radiation source 104 is arranged in the division target area of the area number 9, and the actual responses D 9,10 when the neutron detectors 106 are arranged in order in the division target areas of the area numbers 10, 11, 12 respectively. , D 9,11 , D 9,12 are measured. Further, the actual response is measured by switching the position of the neutron source 104 and the position of the neutron detector 106. That is, the actual response D 10,7 when the neutron source 104 is arranged in the division target area of the area number 10 and the neutron detectors 106 are arranged in order in the division target areas of the area numbers 7,8,9. , D 10,8, to measure the D 10,9. Further, the actual response D 11,7 when the neutron source 104 is arranged in the division target area of the area number 11 and the neutron detectors 106 are sequentially arranged in the division target areas of the area numbers 7, 8, 9 respectively. , D 11,8 and D 11,9 are measured. Further, the actual response D 12,7 when the neutron radiation source 104 is arranged in the division target region of the region number 12 and the neutron detectors 106 are sequentially arranged in the division target regions of the region numbers 7, 8, 9 respectively. , D 12,8, to measure the D 12,9.

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

一方、水分推定装置102では、放射線輸送計算コードに基づく解析を行う。水分推定装置102は、コンクリート構造物200の分割対象領域kについて3次元形状を扱える放射線輸送計算コード、例えば汎用コードであるモンテカルロ法に基づくMCNPやSN法に基づくTORT等を使用して解析する。   On the other hand, the moisture estimation apparatus 102 performs analysis based on the radiation transport calculation code. The moisture estimation apparatus 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 the Monte Carlo method or TORT based on the SN method.

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

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

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

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

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

Figure 2017009371
The non-collisional 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 the neutron, S is the initial intensity of the neutron, P1 is the distance from the neutron source 104 to the reflection point A, and σ T (E) is the total cross-sectional area of the neutron of energy E.
Figure 2017009371

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

Figure 2017009371
The initial 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 the energy changes from the incident energy E to the reflected energy E ′ and changes from the incident angle θ1 to the reflected angle θ2.
Figure 2017009371

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

Figure 2017009371
The flux Φ B (E ′) in the neutron detector 106 can be expressed by Equation (6). Here, B (E ′) is a build-up coefficient indicating the contribution of scattering components in a path from a neutron having energy E ′ to the neutron detector 106 from the reflection point A, and P2 is a distance from the reflection point A to the neutron detector 106. The distance, σ T (E ′), is the total cross section of the neutron with energy E ′.
Figure 2017009371

中性子線源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 with energy E are radiated from position i of neutron source 104, scattered from each division target region k as reflection point A, and contributed by neutrons detected as neutrons with energy E ′ at position j of neutron detector 106 r i, j, k can be calculated by setting the initial neutron intensity S to 1 in the line bundle Φ B (E ′) of the equation (6). That is, from the relationship between the shape and size of the concrete structure 200 and the position i of the neutron source 104 and the position j of the neutron detector 106, the distance P1 from the neutron source 104 to the reflection point A (each division target region k), Since the distance P2 from the reflection point A to the neutron detector 106, the incident angle θ1 and the reflection angle θ2 with respect to the reflection point A are geometrically determined, by substituting these values into the equation (6), the bundle Φ B ( E ′) can be determined as 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 the formulas (4) and (6) are determined based on the element density composition of the concrete constituting the concrete structure 200. Further, the differential scattering cross section σ S (E, E ′, θ1, θ2) at the reflection point A is determined based on the atomic density of water accumulated in the junkers and the like in addition to the element density composition of the concrete structure 200. . Approximately, the differential scattering cross section and the water differential scattering cross section of the concrete constituting the concrete structure 200 can be determined from the volume ratio (volume ratio) at each reflection point A. That is, the differential scattering cross section of the target region containing a considerable amount of water other than the concrete constituting the concrete structure 200 is the differential scattering cross section of the concrete constituting the concrete structure 200 and the differential scattering cross section of water. Is multiplied by the respective volume ratios 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 2017009371
Therefore, assuming that the reflection point A is the concrete constituting the concrete structure 200 (concrete volume ratio 100%), the differential scattering cross section σ SC at the reflection point A determined based on the element density composition of the concrete structure 200 The contribution r i, j, k when (E, E ′, θ1, θ2) is applied is determined as the contribution rc i, j, k . Further, assuming that the reflection point A is water (water volume ratio 100%), the contribution r i, when applying the differential scattering cross section σ SW (E, E ′, θ1, θ2) at the reflection point A of water . Find j, k as contributions rw i, j, k . Then, when the volume ratio ω k of the water contained in the division target region k is the amount of water, the response measured by the neutron detector 106 is expressed by Expression (7).
Figure 2017009371

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

Figure 2017009371
When Formula (7) is expanded, Formula (8) is obtained.
Figure 2017009371

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

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

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

Figure 2017009371
For example, the distribution of moisture in the concrete structure 200 can be estimated by applying the least square method. With respect to Equation (8), the left side constant term C i, j and the residual R of the calculated value on the right side are defined as Equation (9) as an objective function based on the least square method.
Figure 2017009371

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

Figure 2017009371
The steepest descent method is applied to minimize this objective function. As a constraint condition, a condition that the volume ratio ω k of water is 0 or more is set. In general, it is impossible to solve the problem given by Equation (9) by a closed-type operation, and it is necessary to perform some kind of repeated calculation. Therefore, an algorithm is used to obtain a vector sequence ω k (0) , ω k (1) ,..., Ω k (n) such that the numerical formula (10) is converged to the minimum value of the objective function. .
Figure 2017009371

このアルゴリズムは一般に次の三つのステップで構成される。
(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 appropriate initial point ω k (0) of a volume ratio ω k of water is given, and a parameter n representing the number of repetitions is set to 0.
(2) It is determined whether or not the volume ratio ω k (n) of water obtained for the current parameter n is a solution of Equation (9). If the water volume ratio ω k (n) satisfies the determination condition, the calculation is terminated and the water volume ratio ω k (n) is taken as a solution. If the conditions are not met, proceed to the next step.
(3) A new water volume ratio ω k (n + 1) that reduces the objective function f (ω) is obtained based on the calculation result for the parameter n. The parameter n is incremented by 1, and the process returns to step (2).

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

Figure 2017009371
In order to determine the next water volume ratio ω k (n) in the n-th iteration, 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 can be obtained by Expression (11).
Figure 2017009371

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

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

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

水分推定装置102は、上記の方法により求められた水の体積比ωを出力部16に出力する。このとき、調査対象であるコンクリート構造物200の検査対象領域202を分割対象領域kに分割して表示し、分割対象領域kの各々の水の体積比ωを可視化して表示することが好適である。例えば、水の体積比ωを0%〜100%の範囲で色を変化させて表示させる等、ユーザが視覚から水の体積比ωを容易に把握できるような態様で表示させることが好適である。 The moisture estimation device 102 outputs the water volume ratio ω k obtained by the above method to the output unit 16. At this time, it is preferable that the inspection target area 202 of the concrete structure 200 to be investigated is divided into the division target areas k and displayed, and the volume ratio ω k of each water in the division target area k is visualized and displayed. 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 apparatus 102 estimates the moisture of the concrete structure along the flowchart shown in FIG. 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 source 104 and the position j of the neutron detector 106 are set using the parameters when each division target region k is the reflection point A and the volume ratio of the concrete in the division target region k is 100%. The contribution rc i, j, k is calculated for each combination. At this time, 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 source 104 and the position j of the neutron detector 106 are set using parameters when each division target region k is a reflection point A and the volume ratio of water in the division target region k is 100%. The contribution rw i, j, k is calculated for each combination. In step S16, the response D i, j is measured for each combination of the position i of the neutron 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 rw i, j, k analyzed in steps S12 and S14 and the response D i, j measured in step S16 are changed. Based on this, the water volume ratio ω k in each of the division target regions k is estimated by the least square method. In step S20, the water volume ratio ω k obtained in step S18 is visualized and displayed.

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

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

また、上記の記載から明らかなように本実施の形態における水分推定方法においては、特にトレーサを用いずにコンクリート構造物200内の水分の分布を推定することができる。このためトレーサが添加されていない既存のコンクリート構造物に適用することが可能である。また、本実施の形態においては、中性子線源104及び中性子検出器106は、調査対象であるコンクリート構造物200の同一表面側に配置されている例を示したが、コンクリート構造物の一方向からの水分の存在の推定が可能な範囲において、同一表面側でなくてもよい。なお、コンクリート構造物200の一方向とは、コンクリート構造物200の周囲の空間を2つ以上に分けた場合の1つの側からをいう。   Further, as is clear from the above description, in the moisture estimation method in the present embodiment, the moisture distribution in the concrete structure 200 can be estimated without using a tracer. Therefore, it can be applied to an existing concrete structure to which no tracer is added. Moreover, in this Embodiment, although the neutron beam source 104 and the neutron detector 106 showed the example arrange | positioned on the same surface side of the concrete structure 200 which is investigation object, from one direction of a concrete structure As long as it is possible to estimate the presence of water, they may not be on the same surface side. 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 周辺領域。   DESCRIPTION OF SYMBOLS 10 Processing part, 12 Storage part, 14 Input part, 16 Output part, 18 Interface part, 20 Neutron generation source, 22, 24, 26 Shielding wall, 100 Moisture estimation system, 102 Moisture estimation apparatus, 104 Neutron source, 106 Neutron Detector, 200 concrete structure, 202 inspection object area, 204 peripheral area.

Claims (10)

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

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112114349A (en) * 2019-06-21 2020-12-22 清华大学 Radiation analysis method, apparatus and computer-readable storage medium
CN113358674A (en) * 2021-04-01 2021-09-07 西安交通大学 Neutron resonance CT imaging system and method designed for reinforced concrete member

Citations (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
JPH07504748A (en) * 1992-03-13 1995-05-25 トロクスラー、エレクトロニック、ラボラトリーズ、インコーポレーテッド Instruments and methods for measuring physical properties of substances
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
JP2006078464A (en) * 2004-08-11 2006-03-23 Fujitsu Ltd Apparatus and method for measuring small-angle scattering and sample analysis method
CN103336018A (en) * 2013-07-02 2013-10-02 李霞波 Single-pole low-activity nucleon moisture density instrument

Patent Citations (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
JPH07504748A (en) * 1992-03-13 1995-05-25 トロクスラー、エレクトロニック、ラボラトリーズ、インコーポレーテッド Instruments and methods for measuring physical properties of substances
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
JP2006078464A (en) * 2004-08-11 2006-03-23 Fujitsu Ltd Apparatus and method for measuring small-angle scattering and sample analysis method
CN103336018A (en) * 2013-07-02 2013-10-02 李霞波 Single-pole low-activity nucleon moisture density instrument

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
兼松学他: "熱中性子ラジオグラフィによるコンクリート中の給水過程の解明", 日本建築学会構造系論文集, vol. 78, no. 690, JPN6019006366, August 2013 (2013-08-01), pages 1339 - 1347, ISSN: 0003985714 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112114349A (en) * 2019-06-21 2020-12-22 清华大学 Radiation analysis method, apparatus and computer-readable storage medium
CN113358674A (en) * 2021-04-01 2021-09-07 西安交通大学 Neutron resonance CT imaging system and method designed for reinforced concrete member
CN113358674B (en) * 2021-04-01 2024-05-24 西安交通大学 Neutron resonance CT imaging system and method designed for reinforced concrete member

Also Published As

Publication number Publication date
JP6518939B2 (en) 2019-05-29

Similar Documents

Publication Publication Date Title
Vidmar EFFTRAN—A Monte Carlo efficiency transfer code for gamma-ray spectrometry
Santisteban et al. ENGIN-X: a third-generation neutron strain scanner
JP5289542B2 (en) Air dose evaluation apparatus and method
Priyada et al. Intercomparison of gamma scattering, gammatography, and radiography techniques for mild steel nonuniform corrosion detection
Haith et al. Radiographic modelling for NDE of subsea pipelines
JP6518939B2 (en) Method of estimating moisture of concrete structure and system of estimating moisture of concrete structure
JP2018189516A (en) Device, method, and program for calculating surface contamination density distribution
JP5931494B2 (en) Activation radioactivity evaluation method and activation radioactivity evaluation system
Scannavino Jr et al. A graphical tool for an analytical approach of scattering photons by the Compton effect
JP2014228362A (en) Nuclear fuel subcriticality measuring and monitoring system and method
Hykes et al. Radiation source mapping with Bayesian inverse methods
JP5461413B2 (en) Method for quantification of radiation dose and related isodose curve
Miller et al. Combining radiography and passive measurements for radiological threat localization in cargo
Sharma et al. Experimental response function of NaI (Tl) scintillation detector for gamma photons and tomographic measurements for defect detection
Khorsandi et al. Development of image reconstruction for Gamma-ray CT of large-dimension industrial plants using Monte Carlo simulation
Rimlinger et al. Qualification of gamma spectrometry measurement for the radiological characterization of mixed VLLW cables in particle accelerators
Cheng et al. Feasibility study for wax deposition imaging in oil pipelines by PGNAA technique
JP6814049B2 (en) Subcriticality measurement method for nuclear fuel
JP7489094B2 (en) Concentration detection device and concentration detection method
JP2017122656A (en) Radioactive waste measurement system and method
Bomsdorf et al. Quantitative simulation of coherent X-ray scatter measurements on bulk objects
JP5926362B1 (en) Radioactivity concentration measuring apparatus and radioactivity concentration measuring method
ElHoushy et al. Use of gamma ray back scattering method for corrosion assessment in insulated pipes
Wang et al. High-efficiency Monte Carlo simulation based on CADIS method for Gamma Density Measurement
Harara Deposit thickness measurement in pipes by tangential radiography using gamma ray sources

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