[go: up one dir, main page]

JP7451306B2 - Non-destructive structural analysis equipment, non-destructive structural inspection equipment and non-destructive structural analysis methods - Google Patents

Non-destructive structural analysis equipment, non-destructive structural inspection equipment and non-destructive structural analysis methods Download PDF

Info

Publication number
JP7451306B2
JP7451306B2 JP2020094955A JP2020094955A JP7451306B2 JP 7451306 B2 JP7451306 B2 JP 7451306B2 JP 2020094955 A JP2020094955 A JP 2020094955A JP 2020094955 A JP2020094955 A JP 2020094955A JP 7451306 B2 JP7451306 B2 JP 7451306B2
Authority
JP
Japan
Prior art keywords
measurement
destructive
predictor
simulation
structural analysis
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
JP2020094955A
Other languages
Japanese (ja)
Other versions
JP2021189050A (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.)
Toshiba Corp
Toshiba Energy Systems and Solutions Corp
Original Assignee
Toshiba Corp
Toshiba Energy Systems and Solutions Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Toshiba Corp, Toshiba Energy Systems and Solutions Corp filed Critical Toshiba Corp
Priority to JP2020094955A priority Critical patent/JP7451306B2/en
Publication of JP2021189050A publication Critical patent/JP2021189050A/en
Application granted granted Critical
Publication of JP7451306B2 publication Critical patent/JP7451306B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Analysing Materials By The Use Of Radiation (AREA)

Description

本発明の実施形態は、非破壊構造解析装置、これを用いた非破壊構造検査装置、および非破壊構造解析方法に関する。 Embodiments of the present invention relate to a non-destructive structural analysis device, a non-destructive structural inspection device using the same, and a non-destructive structural analysis method.

構造物内部の物質の分布や内部情報を測定する技術として、さまざまな非破壊検査手法が知られている。 Various non-destructive testing methods are known as techniques for measuring the distribution of substances and internal information inside structures.

代表的な技術としてエックス線を用いたラジオグラフィによる二次元透過画像作成、あるいはコンピュータ断層撮影(CT:Computed Tomography)による3次元画像作成による測定が挙げられる。その他に中性子線源を用いたラジオグラフィ、パルス中性子源を用いたパルス中性子イメージング、中性子散乱、回折法、宇宙線ミュオンの透過率を用いた二次元画像測定やミュオン散乱を用いたミュオントモグラフィ等の多数の手法が存在する。 Typical techniques include measurement using radiography using X-rays to create a two-dimensional transmission image, or computed tomography (CT) to create a three-dimensional image. Other methods include radiography using a neutron beam source, pulsed neutron imaging using a pulsed neutron source, neutron scattering, diffraction method, two-dimensional image measurement using cosmic ray muon transmittance, and muon tomography using muon scattering. A large number of methods exist.

これらの手法は物質に対して特定の放射線(光子や粒子)をプローブとして入射し、その放射線が物質中でエネルギーを失うことにより減衰する作用や、物質を構成する粒子との相互作用により散乱する現象を測定することで、物質の情報を間接的に取得する手法である。 These methods inject specific radiation (photons or particles) into a material as a probe, and the radiation attenuates as it loses energy within the material, or is scattered due to interaction with the particles that make up the material. This is a method of indirectly obtaining information about substances by measuring phenomena.

これらの手法は幅広い分野で用いられており、医療の分野ではレントゲンやCTによる診断、工業の分野では構造物の健全性の評価や材料の組成の評価等に用いられている。 These methods are used in a wide range of fields, including diagnosis using X-rays and CT in the medical field, and evaluation of the soundness of structures and composition of materials in the industrial field.

放射線を用いた構造解析手法では、物質内部で、プローブとする粒子が起こす物質との相互作用を解析することにより構造を予測することが可能であるため、測定値を再現するシミュレーション手法としてモンテカルロ法が幅広く用いられている。 In structural analysis methods using radiation, it is possible to predict the structure by analyzing the interaction between the probe particle and the substance inside the substance, so the Monte Carlo method is used as a simulation method to reproduce measured values. is widely used.

モンテカルロ法は、乱数を用いて多数の粒子の物質中での挙動を解析する手法であり、実際の測定における粒子挙動を高い精度で再現することができる。したがって、シミュレーションモデルとして、測定対象となる物質の形状や組成等を正確に再現し、プローブとなる粒子の種類やエネルギーなどの条件を再現することで、様々な手法における測定結果を非常に高い精度で再現することができる。 The Monte Carlo method is a method that uses random numbers to analyze the behavior of a large number of particles in a material, and can reproduce particle behavior in actual measurements with high accuracy. Therefore, as a simulation model, by accurately reproducing the shape and composition of the substance to be measured, and by reproducing conditions such as the type and energy of the probe particles, measurement results from various methods can be obtained with extremely high precision. It can be reproduced with.

これらのシミュレーション技術は、実物での測定を行う前の準備段階での測定条件の検討や、測定結果と予測値との整合性の評価等を行うために利用されている。特に、複雑な測定条件においては、測定結果の解釈が困難である。したがって、目的とする測定結果が得られているかを確認するために、測定値とシミュレーションの比較による詳細な測定結果の分析が必要となる。 These simulation techniques are used to examine measurement conditions in the preparation stage before measuring actual objects, and to evaluate consistency between measurement results and predicted values. Particularly under complex measurement conditions, it is difficult to interpret measurement results. Therefore, in order to confirm whether the desired measurement results are obtained, detailed analysis of the measurement results by comparing the measured values and simulations is required.

特開2013-231700JP2013-231700

前述したように、放射線をプローブとする非破壊検査手法は、モンテカルロ・シミュレーションにより高い精度で測定結果を再現することができる。しかしながら、モンテカルロ・シミュレーションを行うためには、測定体系を詳細に再現するモデルを作成する必要がある。このため、測定対象に対する情報が事前に得られない場合は、シミュレーシによる予測を行うことができないという課題があった。 As mentioned above, the non-destructive testing method using radiation as a probe can reproduce measurement results with high accuracy through Monte Carlo simulation. However, in order to perform Monte Carlo simulation, it is necessary to create a model that reproduces the measurement system in detail. For this reason, there has been a problem in that predictions cannot be made by simulation if information about the measurement target cannot be obtained in advance.

そこで、本発明の実施形態は、測定対象の形状や組成が未知の条件のもとでも、測定対象の情報を得ることを目的とする。 Therefore, an object of the embodiments of the present invention is to obtain information about a measurement target even under conditions where the shape or composition of the measurement target is unknown.

上述の目的を達成するため、本実施形態に係る非破壊構造解析装置は、放射線を用いた対象物の測定の結果から前記対象物の物質情報を得る非破壊構造解析装置であって、前記対象物が存在する測定領域を模擬したシミュレーション空間を複数の3次元セルに分割し、複数の前記3次元セルのそれぞれにランダムに前記物質情報を設定した複数のランダム物質モデルを作成するランダム物質モデル作成部と、数の前記ランダム物質モデルのそれぞれについて非破壊検査測定を再現するシミュレーションを実行するシミュレーション実行部と、前記シミュレーションのそれぞれの結果から前記非破壊検査測定の結果を模擬する実測再現値を算出する計測値演算部と、数の前記ランダム物質モデルとそれに対応する前記実測再現値の複数のセットを用いた機械学習により、前記実測再現値を説明変数とし、前記ランダム物質モデルの前記物質情報を目的変数とする予測器を作成する予測器作成部と、前記予測器を用いて、前記対象物の測定結果に基づいて前記対象物の前記物質情報を得る構造解析部と、を具備することを特徴とする。 In order to achieve the above-mentioned object, the non-destructive structural analysis device according to the present embodiment is a non-destructive structural analysis device that obtains material information of the target object from the results of measurement of the target object using radiation. A random material model that divides a simulation space that simulates a measurement region in which an object exists into a plurality of three-dimensional cells , and creates a plurality of random material models in which the material information is randomly set in each of the plurality of three -dimensional cells. a creation unit; a simulation execution unit that executes a simulation that reproduces the non-destructive inspection measurement for each of the plurality of random material models; and an actual measurement reproduction value that simulates the result of the non-destructive inspection measurement from the result of each of the simulations. A measured value calculation unit that calculates A predictor creation unit that creates a predictor with substance information as a target variable ; and a structure analysis unit that uses the predictor to obtain the substance information of the object based on the measurement results of the object. It is characterized by

第1の実施形態に係る非破壊構造検査装置およびその非破壊構造解析装置の構成を示すブロック図である。1 is a block diagram showing the configuration of a non-destructive structural inspection device and its non-destructive structural analysis device according to a first embodiment; FIG. 第1の実施形態に係る非破壊構造解析装置におけるランダム物質モデル作成のためのセル分割を示す概念的な分割図である。FIG. 3 is a conceptual division diagram showing cell division for creating a random material model in the non-destructive structural analysis apparatus according to the first embodiment. 第1の実施形態に係る非破壊構造解析装置におけるランダム物質モデルの一例を示す概念的なモデル図である。FIG. 3 is a conceptual model diagram showing an example of a random material model in the non-destructive structural analysis apparatus according to the first embodiment. 第1の実施形態に係る非破壊構造検査装置における測定系の構成を示すブロック図である。FIG. 2 is a block diagram showing the configuration of a measurement system in the non-destructive structural inspection apparatus according to the first embodiment. 第1の実施形態に係る非破壊構造解析装置による実施例の対象物を示す立体図である。FIG. 2 is a three-dimensional diagram showing an example object by the non-destructive structural analysis apparatus according to the first embodiment. 第1の実施形態に係る非破壊構造解析方法の手順を示すフロー図である。FIG. 2 is a flow diagram showing the procedure of a non-destructive structural analysis method according to the first embodiment. 第1の実施形態に係る非破壊構造解析装置による実施例のランダム物質モデルの一例を示す立体図である。FIG. 2 is a three-dimensional diagram showing an example of a random material model of an example using the non-destructive structural analysis apparatus according to the first embodiment. 第1の実施形態に係る非破壊構造解析装置による実施例の測定結果を示す立体図である。FIG. 3 is a three-dimensional diagram showing measurement results of an example using the non-destructive structural analysis apparatus according to the first embodiment. 第2の実施形態に係る非破壊構造検査装置における測定系の構成を示すブロック図である。FIG. 2 is a block diagram showing the configuration of a measurement system in a non-destructive structural inspection apparatus according to a second embodiment. 第3の実施形態に係る非破壊構造検査装置における測定系の構成を示すブロック図である。FIG. 3 is a block diagram showing the configuration of a measurement system in a non-destructive structural inspection apparatus according to a third embodiment. 第4の実施形態に係る非破壊構造検査装置およびその非破壊構造解析装置の構成を示すブロック図である。It is a block diagram showing the composition of the non-destructive structure inspection device and its non-destructive structure analysis device concerning a 4th embodiment. 第4の実施形態に係る非破壊構造検査装置における測定系の構成を示すブロック図である。FIG. 3 is a block diagram showing the configuration of a measurement system in a non-destructive structural inspection apparatus according to a fourth embodiment. 第4の実施形態に係る非破壊構造検査装置による物質重量予測結果の例を示すグラフである。非破壊構造検査装置It is a graph which shows an example of the material weight prediction result by the non-destructive structure inspection apparatus based on 4th Embodiment. Non-destructive structural inspection equipment

以下、図面を参照して、本発明の実施形態に係る非破壊構造解析装置および非破壊構造解析方法について説明する。ここで、互いに同一または類似の部分には、共通の符号を付して、重複する説明は省略する。 DESCRIPTION OF THE PREFERRED EMBODIMENTS A non-destructive structural analysis apparatus and a non-destructive structural analysis method according to embodiments of the present invention will be described below with reference to the drawings. Here, parts that are the same or similar to each other are given the same reference numerals, and redundant explanations will be omitted.

[第1の実施形態]
図1は、第1の実施形態に係る非破壊構造検査装置300およびその非破壊構造解析装置100の構成を示すブロック図である。非破壊構造検査装置300は、測定領域10内の対象物5の構造を解析する非破壊構造解析装置100および測定対象とする測定領域10の測定を行う測定系200を有する。なお、本実施形態では、後述するように、測定系200は、ミュオン散乱測定による測定装置210を用いた場合を例にとって示す。
[First embodiment]
FIG. 1 is a block diagram showing the configuration of a non-destructive structural inspection device 300 and its non-destructive structural analysis device 100 according to the first embodiment. The nondestructive structure inspection device 300 includes a nondestructive structure analysis device 100 that analyzes the structure of the object 5 within the measurement region 10 and a measurement system 200 that measures the measurement region 10 to be measured. In this embodiment, as will be described later, the measurement system 200 is exemplified using a measurement device 210 that performs muon scattering measurement.

非破壊構造解析装置100は、予測器関係演算部110、記憶部120、予測器130、構造解析部140、外部入力を受け入れる入力部150および外部へ出力する出力部160を有する。 The non-destructive structural analysis apparatus 100 includes a predictor relation calculation section 110, a storage section 120, a predictor 130, a structural analysis section 140, an input section 150 that receives external input, and an output section 160 that outputs to the outside.

予測器関係演算部110は、ランダム物質モデル作成部111、シミュレーション実行部112、計測値演算部113、予測器作成部114を有する。 The predictor relationship calculation unit 110 includes a random material model creation unit 111, a simulation execution unit 112, a measured value calculation unit 113, and a predictor creation unit 114.

まず、ランダム物質モデル作成部111について説明する。図2は、ランダム物質モデル作成のためのセル分割を示す概念的な分割図である。 First, the random material model creation section 111 will be explained. FIG. 2 is a conceptual division diagram showing cell division for creating a random material model.

ランダム物質モデル作成部111は、測定領域10を模擬するシミュレーション空間10aを、複数のセルまたはボクセル(以下3次元セルと略)11にセル分割し、各3次元セル11に指定した組成を配置する。ここで、セル分割は、各軸方向にすべて等間隔でもよく、各軸方向に等間隔であるが軸方向ごとの間隔は異なっていてもよい。さらに、等間隔ではなくともよい。図2では、7行×6列の分割の例を示している。図2では、図示の便宜上、2次元的な分割の例を示しているが、測定領域10が3次元の体系であれば、多くは、3次元的に複数の3次元セルに分割する。 The random material model creation unit 111 divides the simulation space 10a that simulates the measurement region 10 into a plurality of cells or voxels (hereinafter referred to as three-dimensional cells) 11, and arranges a specified composition in each three-dimensional cell 11. . Here, the cells may be divided at equal intervals in each axial direction, or evenly spaced in each axial direction, but the intervals in each axial direction may be different. Furthermore, the intervals may not be equal. FIG. 2 shows an example of division into 7 rows x 6 columns. Although FIG. 2 shows an example of two-dimensional division for convenience of illustration, if the measurement area 10 is a three-dimensional system, it is often divided three-dimensionally into a plurality of three-dimensional cells.

図3は、ランダム物質モデルの一例を示す概念的なモデル図である。図3では、測定領域10中の各3次元セル11に、白抜きの3次元セルが空気、黒塗りの3次元セルが鉄の2種類の場合を示している。なお、この際、鉄の密度として複数の組成を用いてもよい。また、図3では、空気と鉄の2通りの場合を示したが、測定領域10中に複数の種類の物質が含まれる場合は、体系に含まれる主要な物質を全て設定することで、複数の物質の分布を再現することができる。以下、各3次元セル11における物質種類あるいは同じ物質の密度区分などの情報を、3次元セル11内の物質情報と呼ぶものとする。 FIG. 3 is a conceptual model diagram showing an example of a random material model. FIG. 3 shows two types of three-dimensional cells 11 in the measurement area 10: white three-dimensional cells are air, and black three-dimensional cells are iron. Note that at this time, a plurality of compositions may be used as the iron density. In addition, although FIG. 3 shows two cases of air and iron, if multiple types of substances are included in the measurement area 10, multiple types of substances can be set by setting all the main substances included in the system. The distribution of substances can be reproduced. Hereinafter, information such as material type or density classification of the same material in each three-dimensional cell 11 will be referred to as material information in the three-dimensional cell 11.

ランダム物質モデル作成部111は、3次元セル中に鉄をランダムに分布させたシミュレーションモデルを大量に作成する。測定領域10中の3次元セルの数がn、物質および密度の種類の数がmの場合、組み合わせは、m通りとなる。すなわちランダム物質モデルの数Nは、最大m個となる。 The random material model creation unit 111 creates a large number of simulation models in which iron is randomly distributed in three-dimensional cells. When the number of three-dimensional cells in the measurement area 10 is n, and the number of types of substances and densities is m, there are m n combinations. In other words, the number N of random material models is m n at most.

しかしながら、後述するように、予測器130の作成を機械学習により行うことから、ランダム物質モデルの数Nは、前述の最大m個とする必要はない。ここで、作成するモデルの数は任意であるが、予測器130の作成を十分な精度で行うためには、千通り以上のモデルを作成することが好ましい。設定する3次元セルが十分に細かく、作成するモデルの数が十分に多い場合、ランダム物質モデル作成の手法により測定領域中で成立し得る鉄の分布条件のほとんどを網羅することができる。 However, as will be described later, since the predictor 130 is created by machine learning, the number N of random material models does not need to be the maximum m n pieces described above. Here, the number of models to be created is arbitrary, but in order to create the predictor 130 with sufficient accuracy, it is preferable to create 1,000 or more models. If the three-dimensional cells to be set are sufficiently fine and the number of models to be created is sufficiently large, most of the iron distribution conditions that can be established in the measurement region can be covered by the method of creating a random material model.

次に、シミュレーション実行部112について説明する。 Next, the simulation execution unit 112 will be explained.

シミュレーション実行部112は、ランダム物質モデル作成部111で作成された各ランダム物質モデルについて、モンテカルロ・シミュレーションにより、測定領域10におけるミュオンと対象物5との相互作用、および、それによるミュオンの挙動を再現するシミュレーションを実行する。 The simulation execution unit 112 reproduces the interaction between the muons and the object 5 in the measurement region 10 and the resulting behavior of the muons by Monte Carlo simulation for each random material model created by the random material model creation unit 111. Run a simulation.

ここで、一般的に、対象物5と放射線との相互作用としては、エックス線、ガンマ線等の電磁波および中性子線、ミュオン等の粒子線の減衰、エックス線、中性子線、ミュオン等の散乱および回折などがある。 Here, in general, interactions between the object 5 and radiation include attenuation of electromagnetic waves such as X-rays and gamma rays and particle beams such as neutron beams and muons, and scattering and diffraction of X-rays, neutron beams, and muons. be.

この際、測定系200の構成、特性および配置、測定領域10へのミュオンの照射条件などのシミュレーション条件は、あらかじめ、入力部150に外部入力として受け入れられ、記憶部120に収納、記憶される。 At this time, simulation conditions such as the configuration, characteristics, and arrangement of the measurement system 200 and conditions for muon irradiation onto the measurement region 10 are received in advance as external inputs by the input unit 150 and stored and stored in the storage unit 120.

シミュレーション実行部112は、このシミュレーション条件の下で、準備されたすべてのランダム物質モデルについて、測定系200による測定を再現するシミュレーションを実行する。 The simulation execution unit 112 executes a simulation that reproduces the measurement by the measurement system 200 for all the prepared random material models under this simulation condition.

次に、計測値演算部113について説明する。計測値演算部113は、シミュレーション実行部112でシミュレートされたミュオンと対象物5との相互作用によるミュオンの挙動に基づいて、ミュオン散乱測定の体系の下でミュオン散乱測定の再現により計測値を算出する。すなわち、実際の測定と同様に、検出器位置でのミュオン軌跡のデータが出力される。このデータは実測値と同等の形式を持つことから、実測再現値と呼称する。ランダム物質モデルをN通り作成した場合は、測定を再現するシミュレーションもN通り実行され、実測再現値もN通りの値が出力される。なお、ミュオン散乱測定のみの場合はN通りとなるが、複数の異なるM種類の測定手法を適用する場合は、1つのランダム物質モデルに対してM個の実測再現値を持つデータセットがN組作成される。 Next, the measured value calculation section 113 will be explained. The measured value calculation unit 113 calculates the measured value by reproducing the muon scattering measurement under the muon scattering measurement system, based on the behavior of the muon due to the interaction between the muon and the object 5 simulated by the simulation execution unit 112. calculate. That is, data of the muon trajectory at the detector position is output, similar to the actual measurement. Since this data has the same format as the actual measurement value, it is called the actual measurement reproduction value. When N different random material models are created, N different simulations are executed to reproduce the measurement, and N different actual measurement reproduction values are output. Note that in the case of only muon scattering measurement, there are N types, but when applying multiple M different measurement methods, there are N sets of data sets with M actual measurement reproduction values for one random material model. Created.

次に、予測器作成部114について説明する。予測器作成部114は、測定値を説明変数Xとし、目的とする値を目的変数Yとした関数(Y=f(X))を作成する。ここで、目的とする値は、たとえば、後述する回帰分析の場合には、説明変数である測定値Xに対応するランダム物質モデルのモデル情報である。 Next, the predictor creation unit 114 will be explained. The predictor creation unit 114 creates a function (Y=f(X)) with the measured value as the explanatory variable X and the target value as the objective variable Y. Here, the target value is, for example, model information of a random material model corresponding to the measured value X, which is an explanatory variable, in the case of regression analysis to be described later.

予測器作成部114は、関数f(X)を、機械学習に基づいて作成する。以下に、機械学習における教師あり学習の場合として、回帰分析およびクラス分類の例を示す。 The predictor creation unit 114 creates a function f(X) based on machine learning. Below, examples of regression analysis and class classification are shown as cases of supervised learning in machine learning.

まず、回帰分析の場合を説明する。回帰分析では、目的変数Yを、ランダム物質モデルが持つ数値情報とし、説明変数Xを実測再現値とする。 First, the case of regression analysis will be explained. In regression analysis, the objective variable Y is numerical information possessed by a random material model, and the explanatory variable X is an actual measurement reproduction value.

説明変数Xが1種類である単回帰分析においては、ある変数xから、目的となる値yを予測する。また、説明変数が複数でm個ある重回帰分析では、x1、x2~xmまでのm個の変数から、目的となるyの値を予測する関数を作成する。 In simple regression analysis in which there is one type of explanatory variable X, a target value y is predicted from a certain variable x. In addition, in multiple regression analysis with a plurality of m explanatory variables, a function is created to predict the target value of y from m variables x1, x2 to xm.

回帰分析を用いる場合は、目的変数Yを、測定領域中の鉄の重量や平均密度などの連続量として設定することによって、説明変数Xであるミュオン散乱データから、直接的にモデルが持つ情報を推測することができる。 When using regression analysis, by setting the objective variable Y as a continuous quantity such as the weight or average density of iron in the measurement area, the information held by the model can be directly extracted from the muon scattering data, which is the explanatory variable X. You can guess.

次に、クラス分類の場合を説明する。クラス分類においては、たとえば、ある3次元セルに対して3次元セルの組成が空気か鉄かの2値の分類問題を扱う。あるいは、設定した物質が複数ある場合、同じ物質でも密度のレベルが複数ある場合等においては、多クラス分類問題となる。クラス分類の代表的な手法として、SVM(サポートベクターマシン)が挙げられる。 Next, the case of class classification will be explained. In class classification, for example, we deal with a binary classification problem of whether the composition of a certain three-dimensional cell is air or iron. Alternatively, if there are multiple substances set, or if the same substance has multiple density levels, it becomes a multi-class classification problem. SVM (support vector machine) is a typical method for class classification.

クラス分類を行うことによって、各3次元セルが鉄であるか空気であるかが判定される。この分類を行うことにより、測定値から3次元セル中の鉄の分布を予測することができる。あるいは、他クラス分類の場合は、そのクラスに応じた分類に基づいて3次元セル中の物質の分布を予測することができる。 By performing classification, it is determined whether each three-dimensional cell is iron or air. By performing this classification, the distribution of iron in the three-dimensional cell can be predicted from the measured values. Alternatively, in the case of other class classification, the distribution of substances in the three-dimensional cell can be predicted based on the classification according to the class.

上記の手法以外にも機械学習、ディープラーニングの各手法を用いて、3次元セルの物質の分類、または測定体系の密度や物質の含有率などの特徴量を推定する予測手法を適用することで、測定値を入力としてランダム物質モデルで作成されたモデル情報を出力する予測器を作成することができる。 In addition to the above methods, machine learning and deep learning methods can be used to classify materials in three-dimensional cells, or to apply prediction methods to estimate feature values such as density and material content in measurement systems. , it is possible to create a predictor that inputs measured values and outputs model information created using a random material model.

次に、記憶部120について説明する。記憶部120は、モデル記憶部121、シミュレーション条件記憶部122、および計測値演算結果記憶部123を有する。 Next, the storage unit 120 will be explained. The storage unit 120 includes a model storage unit 121 , a simulation condition storage unit 122 , and a measurement value calculation result storage unit 123 .

まず、モデル記憶部121は、ランダム物質モデル作成部111が作成したすべてのランダム物質モデルを収納、記憶する。それぞれのランダム物質モデルは、図3に概念的に示したような3次元セル11の体系における各3次元セル11内の物質情報の集合として収納される。 First, the model storage section 121 stores and stores all the random material models created by the random material model creation section 111. Each random material model is stored as a set of material information within each three-dimensional cell 11 in a system of three-dimensional cells 11 as conceptually shown in FIG.

次に、シミュレーション条件記憶部122は、外部入力として入力部150が受け入れたシミュレーション条件を収納、記憶する。シミュレーション条件としては、前述のような測定系200の構成、特性および配置、測定領域10へのミュオンの照射条件の外に、3次元セル11の寸法等の情報がある。 Next, the simulation condition storage unit 122 stores and stores the simulation conditions accepted by the input unit 150 as external input. The simulation conditions include information such as the dimensions of the three-dimensional cell 11 in addition to the configuration, characteristics, and arrangement of the measurement system 200 as described above, and muon irradiation conditions to the measurement region 10.

次に、計測値演算結果記憶部123は、計測値演算部113による検出器位置での測定値の演算結果である実測再現値を、それぞれのランダム物質モデルに対応付けた形で収納する。なお、計測値演算結果記憶部123は、シミュレーション実行部112による演算結果等をさらに収納、記憶してもよい。 Next, the measured value calculation result storage unit 123 stores actual measurement reproduction values, which are the calculation results of the measured values at the detector positions by the measured value calculation unit 113, in a form that is associated with each random material model. Note that the measurement value calculation result storage unit 123 may further store and store calculation results etc. by the simulation execution unit 112.

次に、予測器130について説明する。予測器130は、前述のように、予測器作成部114により作成された、測定値を説明変数Xとし、目的とする値を目的変数Yとした関数(Y=f(X))を有している。予測器130が有するこの関数は、前述のように、機械学習により、使用される測定系200の条件の下での測定値(説明変数X)と目的変数Yとの対応付けが可能となっている関数である。 Next, the predictor 130 will be explained. As described above, the predictor 130 has a function (Y=f(X)) created by the predictor creation unit 114 with the measured value as the explanatory variable X and the target value as the objective variable Y. ing. As described above, this function possessed by the predictor 130 is made possible by machine learning to associate the measured value (explanatory variable X) with the objective variable Y under the conditions of the measurement system 200 used. This is a function that has

測定系200で得られた測定結果は、入力部150により受け入れられ、予測器130に入力される。予測器130は、この測定結果を説明変数Xとして、関数f(X)を用いて、目的変数Y、すなわち、測定結果に対応する3次元セル11内の物質情報を出力する。 The measurement results obtained by the measurement system 200 are accepted by the input unit 150 and input to the predictor 130. The predictor 130 uses the measurement result as the explanatory variable X and uses the function f(X) to output the objective variable Y, that is, the material information in the three-dimensional cell 11 corresponding to the measurement result.

次に、構造解析部140は、予測器作成部114により作成された予測器130を用いて、入力部150を介して受け入れた測定系200で測定した測定結果を、関数f(X)の入力として、目的変数Yである物質構造を得る。 Next, the structure analysis unit 140 uses the predictor 130 created by the predictor creation unit 114 to input the measurement results measured by the measurement system 200 received via the input unit 150 into the function f(X). , we obtain the material structure with the objective variable Y.

次に、測定系200について説明する。図4は、非破壊構造検査装置300における測定系200の構成を示すブロック図である。 Next, the measurement system 200 will be explained. FIG. 4 is a block diagram showing the configuration of the measurement system 200 in the non-destructive structural inspection apparatus 300.

本第1の実施形態においては、測定系200として、ミュオン散乱法による測定装置210が用いられている。 In the first embodiment, a measurement device 210 using a muon scattering method is used as the measurement system 200.

ミュオン散乱法は、宇宙線ミュオンが物質を透過する際のミュオン散乱を測定することにより、物質の構造を解析する手法である。ミュオン散乱は、次の式(1)に示されるように、散乱角θが、物質の原子番号に固有の放射長Xに依存して変化する。 The muon scattering method is a method for analyzing the structure of materials by measuring muon scattering when cosmic ray muons pass through materials. In muon scattering, the scattering angle θ 0 changes depending on the radiation length X 0 specific to the atomic number of the substance, as shown in the following equation (1).

θ=[13.6/(βp)]・z・√(x/X)[1+0.038ln(x/X)]
…(1)
なお、β、p、zは、それぞれミュオンの速度、運動量、電荷であり、xはミュオンが物質中を通過する距離である。
θ 0 = [13.6/(β c p)]・z・√(x/X 0 ) [1+0.038ln(x/X 0 )]
...(1)
Note that β c , p, and z are the velocity, momentum, and charge of the muon, respectively, and x is the distance that the muon passes through the substance.

この特徴を利用して、ミュオンの散乱角θの大きさの統計量と散乱位置とを解析することにより、測定領域内の物質の3次元分布を測定することができる。 Utilizing this characteristic, by analyzing the statistical value of the magnitude of the muon scattering angle θ 0 and the scattering position, it is possible to measure the three-dimensional distribution of the substance within the measurement region.

ミュオン散乱法による測定のためには、2つのミュオン検出器211、212が必要となる。すなわち、ミュオン検出器211は、ミュオンが対象物5に入射する前のミュオン軌跡を、また、ミュオン検出器212は、ミュオンが対象物5を通過する際に散乱された後のミュオン軌跡を、それぞれ測定する。 Two muon detectors 211 and 212 are required for measurement using the muon scattering method. That is, the muon detector 211 detects the muon trajectory before the muon enters the object 5, and the muon detector 212 detects the muon trajectory after the muon is scattered when passing through the object 5. Measure.

ここで、測定領域10については、測定に用いる制御された放射線以外には、外部から測定領域10に放射線が入射することを防止するために、測定領域10の周囲は、図示しない遮蔽構造により遮蔽されている。なお、これに代えて、宇宙からのミュオンを測定用の放射線として用いてもよい。この場合は、ミュオンの照射条件を宇宙からの照射の条件とする。 Here, regarding the measurement area 10, in order to prevent radiation from entering the measurement area 10 from the outside other than the controlled radiation used for measurement, the periphery of the measurement area 10 is shielded by a shielding structure (not shown). has been done. Note that instead of this, muons from space may be used as radiation for measurement. In this case, the conditions for muon irradiation are the conditions for irradiation from space.

信号処理部213は、ミュオン検出器211およびミュオン検出器212からの信号を受け入れて、たとえば、増幅等を行った上で、対象物5への入射前後のミュオン軌跡の方向の差分を、ミュオンの散乱角θとして算出する。算出されたミュオンの散乱角θは、入力部150に出力される。 The signal processing unit 213 receives the signals from the muon detector 211 and the muon detector 212, performs amplification, etc., and calculates the difference in the direction of the muon trajectory before and after the muon is incident on the object 5. Calculate as scattering angle θ 0 . The calculated muon scattering angle θ 0 is output to the input section 150.

最後に出力部160は、構造解析部140で解析された結果に基づいて、3次元画像化のためのデータを生成し、表示する。 Finally, the output unit 160 generates and displays data for three-dimensional imaging based on the results analyzed by the structure analysis unit 140.

次に、以上のように構成された非破壊構造検査装置300を用いた非破壊構造解析方法について説明する。図5は、本実施形態に係る非破壊構造解析装置100による実施例の対象物5を示す立体図である。実施例の対象物5は、外形が立方体のアルミニウム製の容器内に球体が収納されたものである。対象物5が設置されている測定領域を模したシミュレーション空間10aは、3次元的に、X軸、Y軸、Z軸の各方向に200分割されており、8×10個の3次元セルに分割されている。以下、本実施形態による非破壊構造解析方法のステップに従って、実施例を説明する。 Next, a nondestructive structural analysis method using the nondestructive structural inspection apparatus 300 configured as described above will be described. FIG. 5 is a three-dimensional diagram showing an example object 5 by the non-destructive structural analysis apparatus 100 according to the present embodiment. The object 5 of the example is a sphere housed in an aluminum container having a cubic outer shape. The simulation space 10a, which simulates the measurement area in which the object 5 is installed, is three-dimensionally divided into 200 parts in each direction of the X-axis, Y-axis, and Z-axis, and is divided into 8×10 6 three-dimensional cells. It is divided into. Examples will be described below according to the steps of the non-destructive structural analysis method according to this embodiment.

図6は、第1の実施形態に係る非破壊構造解析方法の手順を示すフロー図である。図6では、予測器作成関係ステップS10、対象物5の測定ステップS20、および構造解析ステップS30を示している。 FIG. 6 is a flow diagram showing the procedure of the non-destructive structural analysis method according to the first embodiment. FIG. 6 shows a predictor creation step S10, a measurement step S20 of the object 5, and a structure analysis step S30.

非破壊構造解析方法は、この中の、予測器作成関係ステップS10および構造解析ステップS30を有する。 The non-destructive structural analysis method includes a predictor creation related step S10 and a structural analysis step S30.

まず、予測器作成関係ステップS10について説明する。予測器作成関係ステップS10は、ランダム物質モデル作成ステップ(S11)、シミュレーション実行ステップ(S12)、実測値再現ステップ(S13)、および予測器作成ステップ(S14)を有する。 First, the predictor creation related step S10 will be explained. The predictor creation related step S10 includes a random material model creation step (S11), a simulation execution step (S12), an actual measurement value reproduction step (S13), and a predictor creation step (S14).

ランダム物質モデル作成ステップ(S11)においては、ランダム物質モデル作成部111が、複数の3次元セル11に分割されたシミュレーション空間10a上に、測定領域10中で成立し得る条件に基づきランダムに設定されたシミュレーションモデル、すなわちランダム物質モデルを、大量に作成する。作成されたこれらのランダム物質モデルは、モデル記憶部121に収納、記憶される。 In the random material model creation step (S11), the random material model creation section 111 randomly sets a model on the simulation space 10a divided into a plurality of three-dimensional cells 11 based on conditions that can be established in the measurement region 10. Create a large number of simulation models, that is, random material models. These created random material models are stored and stored in the model storage unit 121.

図7は、非破壊構造解析装置100による実施例のランダム物質モデルの一例を示す立体図である。各軸方向にそれぞれ200に分割されたそれぞれの3次元セルについて、物質情報としては、空気、および鉄の区分とした。図7に示したランダム物質モデルの例において、色塗りされたセルは鉄に区分されたセル、色塗りされていないセルは空気に区分されたセルを示す。 FIG. 7 is a three-dimensional diagram showing an example of a random material model according to an embodiment of the non-destructive structural analysis apparatus 100. For each three-dimensional cell divided into 200 cells in each axial direction, the material information was classified into air and iron. In the example of the random material model shown in FIG. 7, colored cells indicate cells classified as iron, and uncolored cells indicate cells classified as air.

シミュレーション実行ステップ(S12)においては、シミュレーション実行部112が、すべてのランダム物質モデルについて、測定系200による測定を再現するシミュレーションを実行する。すなわち、作成された各ランダム物質モデルについて、モンテカルロ・シミュレーションにより、測定領域10におけるミュオンと対象物5との相互作用、および、それによるミュオンの挙動を再現するシミュレーションを実行する。この際、入力部150で受け入れてシミュレーション条件記憶部122で記憶されたシミュレーション条件に基づいて、シミュレーションを実行する。 In the simulation execution step (S12), the simulation execution unit 112 executes a simulation that reproduces the measurement by the measurement system 200 for all random material models. That is, for each created random material model, a Monte Carlo simulation is performed to reproduce the interaction between the muons and the object 5 in the measurement region 10 and the resulting behavior of the muons. At this time, the simulation is executed based on the simulation conditions accepted by the input unit 150 and stored in the simulation condition storage unit 122.

実測値再現ステップ(S13)においては、計測値演算部113が、シミュレーション実行部112でシミュレートされたそれぞれの結果に基づいて、ミュオン散乱測定の体系の下でのミュオン散乱測定の再現による計測値である実測再現値を算出する。この結果、N通りのランダム物質モデルについてN通りの実測再現値が得られる。なお、測定手法がM種類ある場合、このM種類について実測再現値を算出する場合には、M個の実測再現値を持つデータセットがN組作成される。これらの実測再現値は、計測値演算結果記憶部123に収納、記憶される。 In the actual measurement value reproduction step (S13), the measurement value calculation unit 113 calculates the measurement value by reproducing the muon scattering measurement under the muon scattering measurement system, based on each result simulated by the simulation execution unit 112. Calculate the actual measured reproduction value. As a result, N types of actually measured reproduction values are obtained for N types of random material models. Note that in the case where there are M types of measurement methods, when calculating actual measurement reproduction values for these M types, N sets of data sets having M actual measurement reproduction values are created. These actually measured reproduction values are stored and stored in the measured value calculation result storage section 123.

予測器作成ステップ(S14)においては、予測器作成部114が、機械学習に基づいて、測定値を説明変数Xとし、目的とする値を目的変数Yとした関数(Y=f(X))を作成する。ここで、説明変数Xとなる測定値として、計測値演算部113により算出され計測値演算結果記憶部123に記憶されている実測再現値が用いられる。 In the predictor creation step (S14), the predictor creation unit 114 uses a function (Y=f(X)) with the measured value as the explanatory variable X and the target value as the objective variable Y, based on machine learning. Create. Here, as the measured value serving as the explanatory variable X, an actual measured reproduction value calculated by the measured value calculation section 113 and stored in the measured value calculation result storage section 123 is used.

機械学習が、教師あり学習であり回帰分析が用いられる場合は、目的変数Yは、測定領域中の鉄の重量や平均密度などの連続量として設定することができる。 When machine learning is supervised learning and regression analysis is used, the objective variable Y can be set as a continuous quantity such as the weight or average density of iron in the measurement area.

また、機械学習が、教師あり学習でありクラス分類が用いられる場合は、クラスに応じた分類に基づいて3次元セル中の物質の分布を、目的変数Yとすることができる。 Further, when the machine learning is supervised learning and class classification is used, the distribution of substances in the three-dimensional cell can be set as the objective variable Y based on the classification according to the class.

次に、対象物5の測定ステップS20について説明する。 Next, the step S20 of measuring the object 5 will be explained.

測定系200としてのミュオン散乱法による測定装置210により対象物5におけるミュオン散乱が測定される。すなわち、ミュオン検出器211によりミュオンが対象物5に入射する前のミュオン軌跡を、また、ミュオン検出器212によりミュオンが対象物5を通過する際に散乱された後のミュオン軌跡がそれぞれ測定される。これらの結果は信号処理部213に受け入れられ、対象物5への入射前後のミュオン軌跡の方向の差分がミュオンの散乱角θとして算出され、算出されたミュオンの散乱角θが入力部150に出力される。 Muon scattering in the object 5 is measured by a measuring device 210 that uses a muon scattering method as the measuring system 200 . That is, the muon detector 211 measures the muon trajectory before the muon enters the object 5, and the muon detector 212 measures the muon trajectory after the muon is scattered as it passes through the object 5. . These results are accepted by the signal processing unit 213, the difference in the direction of the muon trajectory before and after the incidence on the object 5 is calculated as the muon scattering angle θ 0 , and the calculated muon scattering angle θ 0 is input to the input unit 150. is output to.

次に、構造解析ステップS30について説明する。 Next, the structural analysis step S30 will be explained.

まず、予測器130が、測定系200による対象物5の測定結果に基づいて、物質情報を予測する(ステップS31)。すなわち、予測器130は、測定系200による測定で得られた実測値を入力として受け入れる。測定系200による実測値と、計測値演算部113で算出した実測再現値は同じ形式であるため、予測器130は、実測再現値を入力した場合と同様に、実測値から予測される最も適した物質モデルの情報を出力する。 First, the predictor 130 predicts material information based on the measurement results of the object 5 by the measurement system 200 (step S31). That is, the predictor 130 receives as input the actual measured value obtained by measurement by the measurement system 200. Since the actual measurement value by the measurement system 200 and the actual measurement reproduction value calculated by the measurement value calculation unit 113 are in the same format, the predictor 130 calculates the most suitable value predicted from the actual measurement value in the same way as when inputting the actual measurement reproduction value. Outputs material model information.

次に、構造解析部140が、予測器130による予測の結果から、物質構造を解析する(ステップS32)。 Next, the structure analysis unit 140 analyzes the material structure based on the prediction result by the predictor 130 (step S32).

以上の手順により、対象物5の組成や構造が未知の条件で合っても、測定系200による測定値から、その組成や構造等の内部情報を、予測値として出力することができる。 Through the above procedure, even if the composition and structure of the object 5 meet unknown conditions, internal information such as the composition and structure can be output as a predicted value from the measured values by the measurement system 200.

図8は、本実施形態に係る非破壊構造解析装置による実施例の測定結果を示す立体図である。予測器130の作成時点で、鉄の分布を学習しているため、アルミ容器は3次元画像には現れずに、鉄球のみが画像化されている。 FIG. 8 is a three-dimensional diagram showing measurement results of an example using the non-destructive structural analysis apparatus according to the present embodiment. Since the distribution of iron is learned at the time of creation of the predictor 130, the aluminum container does not appear in the three-dimensional image, and only the iron ball is imaged.

すなわち、学習した鉄の分布に関しては、図5の実際の対象物5と比較して、鉄球の位置、形状、大きさ等が、再現されている。 That is, regarding the learned distribution of iron, the position, shape, size, etc. of the iron ball are reproduced compared to the actual object 5 in FIG.

以上のように、本実施形態によれば、測定対象である対象物5の形状や組成が未知の条件のもとでも、その情報を得ることができる。 As described above, according to the present embodiment, information on the shape and composition of the object 5 to be measured can be obtained even under conditions where the shape and composition are unknown.

[第2の実施形態]
図9は、第2の実施形態に係る非破壊構造検査装置における測定系の構成を示すブロック図である。本第2の実施形態は、第1の実施形態の変形である。
[Second embodiment]
FIG. 9 is a block diagram showing the configuration of a measurement system in a non-destructive structural inspection apparatus according to the second embodiment. The second embodiment is a modification of the first embodiment.

第1の実施形態ではミュオン散乱法に場合を例にとって、非破壊構造検査装置および非破壊構造解析方法を示したが、本第2の実施形態においては、測定系200として中性子ラジオグラフィを用いた測定装置220の場合を示す。 In the first embodiment, a non-destructive structural inspection device and a non-destructive structural analysis method were shown using the muon scattering method as an example, but in the second embodiment, neutron radiography was used as the measurement system 200. The case of measuring device 220 is shown.

測定装置220は、中性子発生装置221、中性子検出器222、および信号処理部223を有する。 The measurement device 220 includes a neutron generator 221, a neutron detector 222, and a signal processing section 223.

中性子発生装置221は、たとえば、研究用原子炉などの原子炉、あるいは加速器である。 The neutron generator 221 is, for example, a nuclear reactor such as a research nuclear reactor, or an accelerator.

中性子発生装置221から発生する中性子を、測定領域10中のサンプルである対象物5に照射し、対象物5に入射し、対象物5の中で減衰または散乱した中性子を中性子検出器222で測定する。 Neutrons generated from the neutron generator 221 are irradiated onto the target object 5, which is a sample in the measurement area 10, and the neutrons that are incident on the target object 5 and attenuated or scattered within the target object 5 are measured by the neutron detector 222. do.

中性子検出器222の出力は、信号処理部223に受け入れられ、増幅され、AD変換され非破壊構造解析装置100の入力部150(図1)に入力される。 The output of the neutron detector 222 is received by the signal processing unit 223, amplified, AD converted, and input to the input unit 150 (FIG. 1) of the non-destructive structural analysis apparatus 100.

一方、非破壊構造解析装置100においては、第1の実施形態と同様に、シミュレーション実行部112が、対象物5に中性子を照射し対象物5を透過することにより減衰または散乱する中性子の挙動のシミュレーションを実施し、計測値演算部113が、実測再現値を算出する。予測器作成部114は、ランダム物質モデルと実測再現値との組み合わせに基づいて機械学習を行い、予測器130を作成する。 On the other hand, in the nondestructive structural analysis apparatus 100, similarly to the first embodiment, the simulation execution unit 112 irradiates the object 5 with neutrons and calculates the behavior of the neutrons that are attenuated or scattered by passing through the object 5. The simulation is performed, and the measured value calculation unit 113 calculates the actual measurement reproduction value. The predictor creation unit 114 performs machine learning based on the combination of the random material model and the actual measured reproduction value, and creates the predictor 130.

以上のように、中性子検出器222による測定結果に基づいて、予測器130が測定結果に対応するランダム物質モデルを出力する点では、第1の実施形態と同様である。 As described above, this embodiment is similar to the first embodiment in that the predictor 130 outputs a random material model corresponding to the measurement results based on the measurement results by the neutron detector 222.

以上のように、中性子の照射によっても測定対象である対象物5の形状や組成が未知の条件のもとでも、その情報を得ることができる。また、例として中性子を用いたラジオグラフィの場合を示したが、プローブをX線やガンマ線に変更した場合でも同等の構成で測定することができる。 As described above, even under conditions where the shape and composition of the object 5 to be measured are unknown, information can be obtained by irradiation with neutrons. Furthermore, although radiography using neutrons has been shown as an example, measurements can be made with the same configuration even when the probe is changed to X-rays or gamma rays.

[第3の実施形態]
図10は、第3の実施形態に係る非破壊構造検査装置における測定系の構成を示すブロック図である。本第3の実施形態は、第2の実施形態の変形である。
[Third embodiment]
FIG. 10 is a block diagram showing the configuration of a measurement system in a non-destructive structural inspection apparatus according to the third embodiment. The third embodiment is a modification of the second embodiment.

第2の実施形態では、サンプル中で減衰または散乱した中性子の測定による中性子ラジオグラフィを用いる場合を例にとって示したが、本第3の実施形態においては、測定系200として中性子放射化分析を用いた測定装置230の場合を示す。 In the second embodiment, neutron radiography is used to measure neutrons attenuated or scattered in a sample. However, in the third embodiment, neutron activation analysis is used as the measurement system 200. The case of the measuring device 230 shown in FIG.

測定装置230は、中性子発生装置221、中性子検出器222、ガンマ線検出器232、および信号処理部233を有する。中性子発生装置221および中性子検出器222は、第2の実施形態と同様である。 The measurement device 230 includes a neutron generator 221, a neutron detector 222, a gamma ray detector 232, and a signal processing section 233. The neutron generator 221 and the neutron detector 222 are the same as those in the second embodiment.

ガンマ線検出器232は、入射中性子に対する対象物5の放射化により発生するガンマ線を測定する。発生するガンマ線スペクトルは物質内の組成により異なるため、ガンマ線スペクトルから測定領域中に含まれる物質の種類を推定することができる。 The gamma ray detector 232 measures gamma rays generated by activation of the object 5 with respect to incident neutrons. Since the generated gamma ray spectrum differs depending on the composition within the substance, the type of substance contained in the measurement region can be estimated from the gamma ray spectrum.

非破壊構造解析装置100においては、第1の実施形態と同様に、シミュレーション実行部112が、対象物5に中性子を照射し対象物5が放射化する挙動のシミュレーションを実施し、計測値演算部113が、実測再現値を算出する。予測器作成部114は、予測器作成部114は、ランダム物質モデルと実測再現値との組み合わせに基づいて機械学習を行い、予測器130を作成する。 In the nondestructive structural analysis apparatus 100, similarly to the first embodiment, the simulation execution unit 112 executes a simulation of the behavior in which the target object 5 is irradiated with neutrons and becomes radioactive, and the measured value calculation unit 113 calculates the actual measurement reproduction value. The predictor creating unit 114 performs machine learning based on the combination of the random material model and the actual measurement reproduction value, and creates the predictor 130.

なお、同様の体系を用いて、中性子を入射し、対象物5中で発生する核分裂中性子を測定する手法により非破壊構造検査を実施することができる。対象物5中で発生する中性子量は、対象物5に含まれる核分裂性物質の量に対応するため、対象物5中の核物質量を評価する場合は、この方法が有効である。 Note that, using a similar system, a non-destructive structural inspection can be performed by a method of injecting neutrons and measuring fission neutrons generated in the object 5. Since the amount of neutrons generated in the object 5 corresponds to the amount of fissile material contained in the object 5, this method is effective when evaluating the amount of nuclear material in the object 5.

[第4の実施形態]
図11は、第4の実施形態に係る非破壊構造検査装置300aおよびその非破壊構造解析装置100aの構成を示すブロック図である。
[Fourth embodiment]
FIG. 11 is a block diagram showing the configuration of a non-destructive structural inspection device 300a and its non-destructive structural analysis device 100a according to the fourth embodiment.

本第4の実施形態における非破壊構造検査装置300aの非破壊構造解析装置100aは、第1の実施形態の変形である。本第4の実施形態における非破壊構造解析装置100aは、2つの予測器130、すなわち、第1の予測器131および第2の予測器132を有する。 The non-destructive structural analysis device 100a of the non-destructive structural inspection device 300a in the fourth embodiment is a modification of the first embodiment. The nondestructive structural analysis apparatus 100a in the fourth embodiment has two predictors 130, that is, a first predictor 131 and a second predictor 132.

本第4の実施形態における測定系200は、ミュオン散乱法測定と中性子測定を組み合わせた測定装置240である。 The measurement system 200 in the fourth embodiment is a measurement device 240 that combines muon scattering measurement and neutron measurement.

図12は、第4の実施形態に係る非破壊構造検査装置における測定系の構成を示すブロック図である。本第4の実施形態の測定系200としての測定装置240は、第2の実施形態における測定装置220および第3の実施形態における測定装置230の組み合わせである。すなわち、測定装置240は、ミュオン検出器211、212、中性子発生装置221、中性子検出器222、および信号処理部243を有するミュオン散乱法測定と中性子測定を組み合わせた測定体系である。なお、ミュオン散乱法測定と中性子測定の組み合わせ以外にも、X線やガンマ線の測定を組み合わせる、あるいはこれらの組み合わせを加えることでもよい。 FIG. 12 is a block diagram showing the configuration of a measurement system in a non-destructive structural inspection apparatus according to the fourth embodiment. The measurement device 240 as the measurement system 200 of the fourth embodiment is a combination of the measurement device 220 in the second embodiment and the measurement device 230 in the third embodiment. That is, the measurement device 240 is a measurement system that combines muon scattering measurement and neutron measurement, and includes muon detectors 211 and 212, a neutron generator 221, a neutron detector 222, and a signal processing section 243. In addition to the combination of muon scattering measurement and neutron measurement, it is also possible to combine X-ray or gamma ray measurement, or add a combination of these.

非破壊構造解析装置100aにおいては、まず、シミュレーション実行部112は、ミュオン散乱法測定の挙動のシミュレーションを実施し、計測値演算部113が、実測再現値を算出する。予測器作成部114は、ランダム物質モデルと実測再現値との組み合わせに基づいて機械学習を行い、第1の予測器131を作成する。次に、シミュレーション実行部112は、中性子測定の挙動のシミュレーションを実施し、計測値演算部113が、実測再現値を算出する。予測器作成部114は、ランダム物質モデルと実測再現値との組み合わせに基づいて機械学習を行い、第2の予測器132を作成する。なお、この順番は、逆でもよいし、並行して行われてもよい。 In the non-destructive structural analysis apparatus 100a, the simulation execution unit 112 first performs a simulation of the behavior of muon scattering measurement, and the measured value calculation unit 113 calculates the actual measurement reproduction value. The predictor creation unit 114 performs machine learning based on the combination of the random material model and the measured reproduction value, and creates the first predictor 131. Next, the simulation execution unit 112 executes a simulation of the behavior of neutron measurement, and the measurement value calculation unit 113 calculates an actual measurement reproduction value. The predictor creation unit 114 performs machine learning based on the combination of the random material model and the actual measured reproduction value, and creates the second predictor 132. Note that this order may be reversed or may be performed in parallel.

第1の予測器131により得られた予測結果と、第2の予測器132により得られた予測結果とが照合され、たとえば、両者の平均化処理等により、非破壊構造解析装置100aとしての予測がなされる。 The prediction result obtained by the first predictor 131 and the prediction result obtained by the second predictor 132 are collated, and, for example, by averaging the two, the prediction result as the non-destructive structural analysis device 100a is calculated. will be done.

図11では、2通りの測定手法に基づいて2つの予測器130を作成する場合を示したが、3通り以上の測定手法の場合も、同様である。 Although FIG. 11 shows a case where two predictors 130 are created based on two measurement techniques, the same applies to the case of three or more measurement techniques.

なお、以上は、複数通りの測定手法のそれぞれについて予測器130を互いに独立に作成して、それぞれの予測器130によって予測を行う場合について示したが、これに限定されない。 In addition, although the case has been described above in which the predictors 130 are created independently for each of a plurality of measurement methods and prediction is performed by each predictor 130, the present invention is not limited to this.

すなわち、複数手法を統合した予測器130を作成することができる。複数手法を統合した予測器130を作成する場合は、予測器作成部114において、単回帰分析を行う代わりに、重回帰分析をすることによって、複数の説明変数X(X1,…,Xm)から一つの目的変数Yを出力する。 That is, it is possible to create a predictor 130 that integrates multiple methods. When creating a predictor 130 that integrates multiple methods, the predictor creation unit 114 performs multiple regression analysis instead of simple regression analysis to calculate the results from multiple explanatory variables X (X1,...,Xm). Output one objective variable Y.

異なる測定手法では、異なる情報が得られるため、組み合わせる測定方法を増やすほど、より精度の高い構造解析が可能となる。 Different measurement methods yield different information, so the more measurement methods that are combined, the more accurate structural analysis becomes possible.

図13は、第4の実施形態に係る非破壊構造検査装置300aによる物質重量予測結果の例を示すグラフである。図13は、鉄とウランが混合した物質中のウランの重量についてのグラフであり、横軸は、真のUO重量[kg]、縦軸は、真のUO重量[kg]の予測値を示す。 FIG. 13 is a graph showing an example of the material weight prediction result by the non-destructive structure inspection apparatus 300a according to the fourth embodiment. FIG. 13 is a graph of the weight of uranium in a material containing a mixture of iron and uranium, where the horizontal axis is the true UO2 weight [kg], and the vertical axis is the predicted value of the true UO2 weight [kg]. shows.

図13は、ミュオン散乱法、中性子計数率測定、ガンマ線測定の3通りの手法を組み合わせることにより予測した結果である。 FIG. 13 shows the results predicted by combining three methods: muon scattering, neutron count rate measurement, and gamma ray measurement.

このような評価の結果から複数の物質が混合した状態において、その中の一部の物質の量を抽出する目的でも本実施形態を適用することができる。 Based on the results of such evaluation, the present embodiment can also be applied to the purpose of extracting the amount of some of the substances in a state where a plurality of substances are mixed.

[その他の実施形態]
以上、本発明の実施形態を説明したが、実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。たとえば、実施形態では、放射線を用いた計測方法の場合を例にとって示しているがこれに限定されず、たとえば電磁気的方法、あるいは超音波を用いた方法など、その他の非破壊方式による計測方法であってもよい。
[Other embodiments]
Although the embodiments of the present invention have been described above, the embodiments are presented as examples and are not intended to limit the scope of the invention. For example, in the embodiment, a measurement method using radiation is shown as an example, but the measurement method is not limited to this, and other non-destructive measurement methods such as an electromagnetic method or a method using ultrasonic waves may be used. There may be.

また、各実施形態の特徴を組み合わせてもよい。また、実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。 Moreover, the features of each embodiment may be combined. Further, the embodiments can be implemented in various other forms, and various omissions, substitutions, and changes can be made without departing from the gist of the invention.

実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。 The embodiments and their modifications are included within the scope and gist of the invention as well as within the scope of the invention described in the claims and its equivalents.

5…対象物、10…測定領域、10a…シミュレーション空間、11…3次元セル、100…非破壊構造解析装置、110…予測器関係演算部、111…ランダム物質モデル作成部、112…シミュレーション実行部、113…計測値演算部、114…予測器作成部、120…記憶部、121…モデル記憶部、122…シミュレーション条件記憶部、123…計測値演算結果記憶部、130…予測器、140…構造解析部、150…入力部、160…出力部、200…測定系、210…測定装置、211、212…ミュオン検出器、213…信号処理部、220…測定装置、221…中性子発生装置、222…中性子検出器、223…信号処理部、230…測定装置、232…ガンマ線検出器、233…信号処理部、240…測定装置、243…信号処理部、300…非破壊構造検査装置 5... Object, 10... Measurement region, 10a... Simulation space, 11... Three-dimensional cell, 100... Non-destructive structural analysis device, 110... Predictor relation calculation section, 111... Random material model creation section, 112... Simulation execution section , 113... Measured value calculation section, 114... Predictor creation section, 120... Storage section, 121... Model storage section, 122... Simulation condition storage section, 123... Measured value calculation result storage section, 130... Predictor, 140... Structure Analysis section, 150... Input section, 160... Output section, 200... Measurement system, 210... Measurement device, 211, 212... Muon detector, 213... Signal processing section, 220... Measurement device, 221... Neutron generator, 222... Neutron detector, 223...Signal processing section, 230...Measuring device, 232...Gamma ray detector, 233...Signal processing section, 240...Measuring device, 243...Signal processing section, 300...Non-destructive structure inspection device

Claims (11)

放射線を用いた対象物の測定の結果から前記対象物の物質情報を得る非破壊構造解析装置であって、
前記対象物が存在する測定領域を模擬したシミュレーション空間を複数の3次元セルに分割し、複数の前記3次元セルのそれぞれにランダムに前記物質情報を設定した複数のランダム物質モデルを作成するランダム物質モデル作成部と、
数の前記ランダム物質モデルのそれぞれについて非破壊検査測定を再現するシミュレーションを実行するシミュレーション実行部と、
前記シミュレーションのそれぞれの結果から前記非破壊検査測定の結果を模擬する実測再現値を算出する計測値演算部と、
数の前記ランダム物質モデルとそれに対応する前記実測再現値の複数のセットを用いた機械学習により、前記実測再現値を説明変数とし、前記ランダム物質モデルの前記物質情報を目的変数とする予測器を作成する予測器作成部と、
前記予測器を用いて、前記対象物の測定結果に基づいて前記対象物の前記物質情報を得る構造解析部と、
を具備することを特徴とする非破壊構造解析装置。
A non-destructive structural analysis device that obtains material information of an object from the results of measurement of the object using radiation,
A random method in which a simulation space simulating a measurement region in which the object exists is divided into a plurality of three-dimensional cells , and a plurality of random material models are created in which the material information is randomly set in each of the plurality of three -dimensional cells. Material model creation department,
a simulation execution unit that executes a simulation that reproduces non-destructive inspection measurements for each of the plurality of random material models;
a measurement value calculation unit that calculates an actual measurement reproduction value that simulates the result of the non-destructive inspection measurement from each result of the simulation;
A predictor that uses the actual measured reproduction value as an explanatory variable and the substance information of the random material model as an objective variable, using machine learning using a plurality of sets of the plurality of random material models and the corresponding actual measurement reproduction values. a predictor creation unit that creates
a structural analysis unit that uses the predictor to obtain the material information of the object based on the measurement results of the object;
A non-destructive structural analysis device comprising:
前記シミュレーション実行部は、前記ランダム物質モデル作成部で作成された前記ランダム物質モデルを用いて、非破壊検査手法をモンテカルロ・シミュレーションにより再現することを特徴とする請求項1に記載の非破壊構造解析装置。 The non-destructive structural analysis according to claim 1, wherein the simulation execution unit reproduces a non-destructive inspection method by Monte Carlo simulation using the random material model created by the random material model creation unit. Device. 前記予測器作成部は、
前記ランダム物質モデルとこれに対応する前記実測再現値を一組とするデータセットを複数作成し、
前記ランダム物質モデルを前記目的変数とし、
前記実測再現値を前記説明変数として、
回帰分析またはクラス分類を用いた機械学習を行う、
ことを特徴とする請求項1または請求項2記載の非破壊構造解析装置。
The predictor creation unit includes:
creating a plurality of data sets including the random material model and the corresponding actual measurement reproduction value;
The random material model is the objective variable,
With the actual measured reproduction value as the explanatory variable,
Perform machine learning using regression analysis or class classification,
The non-destructive structural analysis device according to claim 1 or 2, characterized in that:
前記放射線は、エックス線、ガンマ線等の電磁波または中性子線、電子線、ミュオン、陽子線の少なくともいずれか一つの粒子線であることを特徴とする請求項1ないし請求項3のいずれか一項に記載の非破壊構造解析装置。 According to any one of claims 1 to 3, the radiation is an electromagnetic wave such as an X-ray or a gamma ray, or a particle beam of at least one of a neutron beam, an electron beam, a muon beam, and a proton beam. non-destructive structural analysis equipment. 前記請求項1ないし請求項4のいずれか一項に記載の非破壊構造解析装置と、
前記対象物の前記測定を行う測定装置と、
を備えることを特徴とする非破壊構造検査装置。
The non-destructive structural analysis device according to any one of claims 1 to 4,
a measuring device that performs the measurement of the object;
A non-destructive structural inspection device comprising :
前記測定装置は、測定用の前記放射線以外に外部から前記測定領域に入射する放射線を遮蔽する構造を有することを特徴とする請求項5に記載の非破壊構造検査装置。 6. The non-destructive structural inspection device according to claim 5 , wherein the measuring device has a structure that shields radiation that enters the measurement region from the outside in addition to the radiation for measurement . 前記測定装置は、測定用の前記放射線の種類に応じて、放射性同位体、加速器、または核融合装置の少なくとも一つを有することを特徴とする請求項5または請求項6に記載の非破壊構造検査装置。 The nondestructive method according to claim 5 or 6, wherein the measurement device includes at least one of a radioisotope, an accelerator, or a nuclear fusion device, depending on the type of the radiation for measurement. Structural inspection equipment. 前記測定装置は、ミュオンとして大気中で自然に生成される宇宙線ミュオンを測定することを特徴とする請求項または請求項に記載の非破壊構造検査装置。 7. The nondestructive structural inspection device according to claim 5 , wherein the measurement device measures cosmic ray muons that are naturally generated in the atmosphere as muons . 前記測定装置は、前記測定領域中を通過するエックス線、ガンマ線等の電磁波および、中性子線、ミュオン等の粒子線のいずれかの減衰を測定することを特徴とする請求項5ないし請求項8のいずれか一項に記載の非破壊構造検査装置。 9. The measuring device measures attenuation of any of electromagnetic waves such as X-rays and gamma rays, and particle beams such as neutron beams and muons, which pass through the measurement region. The non-destructive structural inspection device according to item (1) . 前記測定装置は、前記測定領域中を通過するエックス線、中性子線、ミュオン等のいずれかの散乱および回折を測定することを特徴とする請求項ないし請求項のいずれか一項に記載の非破壊構造検査装置。 9. The measuring device measures scattering and diffraction of any of X-rays , neutron beams, muons, etc. passing through the measurement region. non-destructive structural inspection equipment. 放射線を用いた対象物の測定の結果から前記対象物の物質情報を得る非破壊構造解析方法であって、
ランダム物質モデル作成部が、前記対象物が存在する測定領域を模擬したシミュレーション空間を複数の3次元セルに分割し、複数の前記3次元セルのそれぞれにランダムに前記物質情報を設定した複数のランダム物質モデルを作成するするランダム物質モデル作成ステップと、
シミュレーション実行部が、複数の前記ランダム物質モデルのそれぞれについて非破壊検査測定を再現するシミュレーションを実行するシミュレーション実行ステップと、
計測値演算部が、前記シミュレーションのそれぞれの結果から前記非破壊検査測定の結果を模擬する実測再現値を算出する計測値演算ステップと、
予測器作成部が、前記ランダム物質モデルとそれに対応する前記実測再現値の複数のセットを用いた機械学習により、前記実測再現値を説明変数とし、前記ランダム物質モデルの前記物質情報を目的変数とする予測器を作成する予測器作成ステップと、
前記予測器作成ステップの後に、構造解析部が、前記予測器をもちいて前記対象物の測定結果に基づいて前記物質情報を予測する予測ステップと、
を有することを特徴とする非破壊構造解析方法
A non-destructive structural analysis method for obtaining material information of an object from the results of measurement of the object using radiation, the method comprising:
A random material model creation unit divides a simulation space simulating the measurement region in which the target object exists into a plurality of three-dimensional cells, and randomly sets the material information in each of the plurality of three-dimensional cells. a random material model creation step for creating a material model;
a simulation execution step in which the simulation execution unit executes a simulation to reproduce non-destructive inspection measurements for each of the plurality of random material models;
a measurement value calculation step in which the measurement value calculation unit calculates an actual measurement reproduction value that simulates the result of the non-destructive inspection measurement from each result of the simulation;
The predictor creation unit uses machine learning using the plurality of sets of the random material model and the corresponding actual measurement reproduction values to set the actual measurement reproduction value as an explanatory variable and the substance information of the random material model as an objective variable. a predictor creation step of creating a predictor for
After the predictor creation step, a prediction step in which the structure analysis unit predicts the material information based on the measurement results of the object using the predictor;
A non-destructive structural analysis method characterized by having the following .
JP2020094955A 2020-05-29 2020-05-29 Non-destructive structural analysis equipment, non-destructive structural inspection equipment and non-destructive structural analysis methods Active JP7451306B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2020094955A JP7451306B2 (en) 2020-05-29 2020-05-29 Non-destructive structural analysis equipment, non-destructive structural inspection equipment and non-destructive structural analysis methods

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2020094955A JP7451306B2 (en) 2020-05-29 2020-05-29 Non-destructive structural analysis equipment, non-destructive structural inspection equipment and non-destructive structural analysis methods

Publications (2)

Publication Number Publication Date
JP2021189050A JP2021189050A (en) 2021-12-13
JP7451306B2 true JP7451306B2 (en) 2024-03-18

Family

ID=78849294

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020094955A Active JP7451306B2 (en) 2020-05-29 2020-05-29 Non-destructive structural analysis equipment, non-destructive structural inspection equipment and non-destructive structural analysis methods

Country Status (1)

Country Link
JP (1) JP7451306B2 (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005114475A (en) 2003-10-06 2005-04-28 Rigaku Corp Method and apparatus for analyzing film structure
JP2013231700A (en) 2012-05-01 2013-11-14 Tokyo Electron Ltd X-ray inspection method and x-ray inspection device
JP2014044079A (en) 2012-08-24 2014-03-13 Sumitomo Bakelite Co Ltd Evaluation method for transparent composite sheet, transparent composite sheet product, and transparent composite sheet
JP2014523522A (en) 2011-06-07 2014-09-11 アトミック エナジー オブ カナダ リミテッド Detection of high atomic number materials using cosmic ray muon tomography
JP2017044642A (en) 2015-08-28 2017-03-02 株式会社東芝 Structure inspection device and inspection method therefor
JP2018036156A (en) 2016-08-31 2018-03-08 株式会社東芝 Heavy element inclusion measuring apparatus and measuring method thereof
WO2019030449A1 (en) 2017-08-10 2019-02-14 Safran Method for non-destructive automatic inspection of mechanical parts
WO2020002704A1 (en) 2018-06-29 2020-01-02 Universiteit Antwerpen Item inspection by dynamic selection of projection angle
JP2021081843A (en) 2019-11-15 2021-05-27 東レ株式会社 Design support device, design support method, and design support program

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005114475A (en) 2003-10-06 2005-04-28 Rigaku Corp Method and apparatus for analyzing film structure
JP2014523522A (en) 2011-06-07 2014-09-11 アトミック エナジー オブ カナダ リミテッド Detection of high atomic number materials using cosmic ray muon tomography
JP2013231700A (en) 2012-05-01 2013-11-14 Tokyo Electron Ltd X-ray inspection method and x-ray inspection device
JP2014044079A (en) 2012-08-24 2014-03-13 Sumitomo Bakelite Co Ltd Evaluation method for transparent composite sheet, transparent composite sheet product, and transparent composite sheet
JP2017044642A (en) 2015-08-28 2017-03-02 株式会社東芝 Structure inspection device and inspection method therefor
JP2018036156A (en) 2016-08-31 2018-03-08 株式会社東芝 Heavy element inclusion measuring apparatus and measuring method thereof
WO2019030449A1 (en) 2017-08-10 2019-02-14 Safran Method for non-destructive automatic inspection of mechanical parts
WO2020002704A1 (en) 2018-06-29 2020-01-02 Universiteit Antwerpen Item inspection by dynamic selection of projection angle
JP2021081843A (en) 2019-11-15 2021-05-27 東レ株式会社 Design support device, design support method, and design support program

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
北村 尚斗,中性子・放射光X線全散乱測定による結晶性酸化物の欠陥分布シミュレーション,粉体および粉末冶金,2017年09月,第64巻、第9号,pp. 489-494,https://www.jstage.jst.go.jp/article/jjspm/64/9/64_489/_pdf/-char/ja,DOI:https://doi.org/10.2497/jjspm.64.489

Also Published As

Publication number Publication date
JP2021189050A (en) 2021-12-13

Similar Documents

Publication Publication Date Title
Parker et al. The use of ionising radiation to image nuclear fuel: A review
Chatzidakis et al. Analysis of spent nuclear fuel imaging using multiple coulomb scattering of cosmic muons
Priyada et al. Intercomparison of gamma scattering, gammatography, and radiography techniques for mild steel nonuniform corrosion detection
Craft et al. Characterization of a neutron beam following reconfiguration of the Neutron Radiography Reactor (NRAD) core and addition of new fuel elements
JP2018036156A (en) Heavy element inclusion measuring apparatus and measuring method thereof
Jafari et al. A Monte Carlo evaluation of neutron images quality in a research reactor based neutron radiography facility
JP2020094906A (en) Nuclear material quantity measuring apparatus and nuclear material quantity measuring method
US20060193421A1 (en) Monitoring
JP7451306B2 (en) Non-destructive structural analysis equipment, non-destructive structural inspection equipment and non-destructive structural analysis methods
Herald et al. Improving the accuracy of PEPT algorithms through dynamic parameter optimisation
Morichi et al. Developments and experiences of the CHANCE, MICADO and PREDIS projects in radioactive waste characterization
Kilby et al. Multi-modal tomographic imaging system for poolside characterization of nuclear test fuels: Design considerations and studies
Liu et al. Simulated imaging of spent nuclear fuel using associated-particle-neutron-induced gamma rays
Shaver et al. Radiation detection computational benchmark scenarios
Ocampo Giraldo et al. Design and validation of a gamma-ray scanning system for measuring irradiated nuclear fuel
Amann et al. High-energy test of proton radiography concepts
Hartman et al. Combined photon–neutron radiography for nondestructive analysis of materials
Núñez-Chongo et al. Advancing Neutron Safety and Dosimetry in Nuclear Facilities: Applications and Current Status of the Development of Nereida
Herald Development of Monte Carlo models for the optimisation of positron emission particle tracking experiments
Chung et al. Preliminary results of a spent nuclear fuel internal tomography (SFIT) system
Jin Characterization of the acquisition parameters of a submersible gamma-ray computed tomography system
Laloy et al. First steps towards productionizing probabilistic radwaste characterization
Hansen Application guide to tomographic gamma scanning of uranium and plutonium
Micklich et al. Key research issues in the pulsed fast-neutron analysis technique for cargo inspection
Paaren Development of a switchable radioisotope generator

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230206

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20231113

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20231121

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20240118

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20240306

R150 Certificate of patent or registration of utility model

Ref document number: 7451306

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150