JP6037968B2 - Radiation measurement apparatus and radiation measurement method - Google Patents
Radiation measurement apparatus and radiation measurement method Download PDFInfo
- Publication number
- JP6037968B2 JP6037968B2 JP2013163859A JP2013163859A JP6037968B2 JP 6037968 B2 JP6037968 B2 JP 6037968B2 JP 2013163859 A JP2013163859 A JP 2013163859A JP 2013163859 A JP2013163859 A JP 2013163859A JP 6037968 B2 JP6037968 B2 JP 6037968B2
- Authority
- JP
- Japan
- Prior art keywords
- response function
- measurement object
- measured
- radiation
- correction
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000005855 radiation Effects 0.000 title claims description 90
- 238000005259 measurement Methods 0.000 title claims description 79
- 238000000691 measurement method Methods 0.000 title claims description 9
- 238000005316 response function Methods 0.000 claims description 137
- 238000009826 distribution Methods 0.000 claims description 71
- 238000012937 correction Methods 0.000 claims description 70
- 238000004364 calculation method Methods 0.000 claims description 36
- 238000007689 inspection Methods 0.000 claims description 32
- 238000001514 detection method Methods 0.000 claims description 30
- 239000000941 radioactive substance Substances 0.000 claims description 20
- 238000001228 spectrum Methods 0.000 claims description 16
- 238000000034 method Methods 0.000 claims description 11
- 239000000126 substance Substances 0.000 claims 2
- 238000003860 storage Methods 0.000 description 18
- 235000013305 food Nutrition 0.000 description 8
- 238000012545 processing Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000010354 integration Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 239000012857 radioactive material Substances 0.000 description 4
- 230000035945 sensitivity Effects 0.000 description 4
- 238000004980 dosimetry Methods 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 2
- 238000010298 pulverizing process Methods 0.000 description 2
- 238000011002 quantification Methods 0.000 description 2
- 235000002566 Capsicum Nutrition 0.000 description 1
- 241000758706 Piperaceae Species 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 238000000084 gamma-ray spectrum Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Landscapes
- Measurement Of Radiation (AREA)
Description
この発明は、測定対象の試料である被測定体を破壊することなく、被測定体の放射能を測定する放射能分析装置および放射能分析方法に関する。 The present invention relates to a radioactivity analyzer and a radioactivity analysis method for measuring the radioactivity of a measurement object without destroying the measurement object which is a sample to be measured.
食品、水、土壌等に含まれる放射性物質の分析や定量を行う放射能分析装置は、測定対象の試料である被測定体からの放射線を検出し、そのエネルギーおよび強度を測定する。すなわち、測定対象核種が発するエネルギーを持つ放射線の強度を測定することにより、分析や定量が行われる。 A radioactivity analyzer that analyzes and quantifies radioactive substances contained in food, water, soil, and the like detects radiation from a measurement object that is a sample to be measured, and measures its energy and intensity. That is, analysis and quantification are performed by measuring the intensity of radiation having energy generated by the measurement target nuclide.
ここで、放射線のエネルギーおよび強度を高精度に測定する方法として、検出器の応答関数を用いた逆問題演算を適用することができる。具体的には、測定されたパルス波高分布から光子エネルギースペクトルを求めるために、応答分布行列による逆問題演算(アンフォールディング)を行い、光子エネルギーの強度を測定して放射性物質を定量する分光光子線量測定法が提案されている(例えば、特許文献1参照)。 Here, as a method for measuring the energy and intensity of radiation with high accuracy, inverse problem calculation using the response function of the detector can be applied. Specifically, in order to obtain the photon energy spectrum from the measured pulse wave height distribution, inverse problem calculation (unfolding) using the response distribution matrix is performed, and the intensity of the photon energy is measured to quantify the radioactive substance. A measurement method has been proposed (see, for example, Patent Document 1).
この分光光子線量測定法では、まず、ある決まった形状や密度の被測定体を定められた測定場所に配置し、ある場所に固定された検出器がどのような応答を示すかを、あらかじめシミュレーション等により求めておく。 In this spectrophoton dosimetry method, first, a measured object of a certain shape and density is placed at a predetermined measurement location, and the response of a detector fixed at a certain location is simulated in advance. Etc.
また、この分光光子線量測定法では、放射線のエネルギーを様々な値に変化させたときの検出器の応答関数を、応答分布行列としてデータベースに蓄えておき、実際に被測定体を測定したときの結果に対して、この応答分布行列を用いた逆問題演算を行うことにより、入射する放射線のエネルギーおよび強度を求めている。 In this spectrophoton dosimetry method, the response function of the detector when the radiation energy is changed to various values is stored in a database as a response distribution matrix, and the measured object is actually measured. The energy and intensity of the incident radiation are obtained by performing an inverse problem calculation using this response distribution matrix on the result.
なお、食品に含まれる放射性物質による検出器からの出力は、食品の形状に依存し、その形状は不定形なので、測定毎にばらつきを生じる。そこで、測定毎のばらつきを低減するために、食品を粉砕してマリネリ容器等の特定の形状を有する容器に充填し、GeやNaIといったシンチレータで、γ線スペクトルを測定して含有する放射能を定量していた。そのため、食品については主に抜き取りの破壊検査となり、全数検査が困難であった。 In addition, since the output from the detector by the radioactive substance contained in the food depends on the shape of the food, and the shape is indefinite, variation occurs in each measurement. Therefore, in order to reduce the variation in each measurement, the food is crushed and filled into a container having a specific shape such as a marinelli container, and the radioactivity contained by measuring the γ-ray spectrum with a scintillator such as Ge or NaI. It was quantified. For this reason, food was mainly subjected to sampling destructive inspection, and 100% inspection was difficult.
これに対して、形状の不定な被測定体に対して粉砕等破壊することなく放射性物質を定量する方法として、被測定体が設置される測定範囲を格子状に分割し、その個々の格子点に点状の放射線源をおいた場合に得られる検出器の感度を計測し、被測定体の3次元形状内に含まれる格子点での各々の検出器の感度を積分して、その体積積分空間を体積線源とした検出効率を得る放射能の測定方法が提案されている(例えば、特許文献2参照)。 On the other hand, as a method of quantifying radioactive substances without pulverizing or destroying the object to be measured with an indefinite shape, the measurement range in which the object to be measured is installed is divided into grids, and the individual lattice points are divided. Measure the sensitivity of the detector obtained when a point-like radiation source is placed at the point, integrate the sensitivity of each detector at the lattice points included in the three-dimensional shape of the object to be measured, and integrate the volume There has been proposed a radioactivity measurement method that obtains detection efficiency using space as a volume radiation source (see, for example, Patent Document 2).
しかしながら、従来技術には、以下のような課題がある。
特許文献1に示された分光光子線量測定法では、被測定体の形状や密度、幾何学的配置を仮定した上で応答分布行列を作成しているので、被測定体の形状や密度が仮定していたものと異なる場合には、正確に放射線のエネルギーおよび強度を測定することができないという問題がある。
However, the prior art has the following problems.
In the spectrophoton dosimetry shown in
すなわち、被測定体の形状や密度が様々に変化する場合、特に生体のように放射性物質の蓄積に部位依存性があるような場合には、逆問題演算(アンフォールディング)により得られたエネルギースペクトルの精度が低下し、測定対象核種の放射性物質の定量精度が著しく低下するという問題がある。 In other words, when the shape and density of the object to be measured change variously, especially when the accumulation of radioactive material is site-dependent, such as a living body, the energy spectrum obtained by inverse problem calculation (unfolding) There is a problem that the accuracy of quantification of the radioactive material of the target nuclide is significantly reduced.
また、特許文献2に示された放射能の測定方法では、形状による感度分布依存性を補正することはできるものの、均一の線源分布を仮定し、検出効率のみの体積積分を行って一点の代表点に転換しているので、表面汚染や線源が生体の一部に蓄積されるような、放射線源に密度分布を持つような被測定体に適用することができないという問題がある。
In addition, although the radioactivity measurement method disclosed in
また、この放射線の測定方法は、逆問題演算に適用することができず、エネルギー分析精度に課題があるとともに、線源に対するエネルギースペクトル分布の線源位置依存性に対応することができず、逆問題演算によるスペクトル分解能の向上、感度の向上に供することができないという問題もある。 In addition, this radiation measurement method cannot be applied to inverse problem calculation, has a problem in energy analysis accuracy, cannot cope with the dependence of the energy spectrum distribution on the source on the source position, and There is also a problem that spectral resolution cannot be improved by problem calculation and sensitivity cannot be improved.
この発明は、上記のような課題を解決するためになされたものであり、試料毎に形状が不定形の被測定体に対しても、最適な検出器応答関数を算出し、この応答関数を用いて逆問題演算を行うことにより、被測定体を破壊することなく全数検査を行うとともに、被測定体から放射される放射線のエネルギースペクトルを高精度に測定することができる放射線測定装置および放射線測定方法を得ることを目的とする。 The present invention has been made to solve the above-described problems. An optimum detector response function is calculated even for an object having an indefinite shape for each sample, and this response function is calculated. used by performing arithmetic inverse problem, performs total inspection without destroying the object to be measured, radiation can be measured energy spectrum of the radiation emitted from the object to be measured with high accuracy measuring device and the radiation The purpose is to obtain a line measurement method.
この発明に係る放射線測定装置は、被測定体から放射される放射線を検出し、出力パルスに応じたパルス波高分布を出力する放射線検出部と、被測定体の形状情報を取得する被測定体情報取得部と、被測定体情報取得部で取得された被測定体の形状情報に基づいて、被測定体の形状情報に応じた応答関数の補正値を、補正応答関数として演算する補正応答関数生成部と、を備え、放射線検出部から出力されたパルス波高分布に対して、補正応答関数生成部で演算された補正応答関数を用いた逆問題演算を行うことにより、被測定体から放射される放射線のエネルギースペクトルを演算する放射線測定装置であって、被測定体情報取得部は、検査領域内の少なくとも1つの被測定体の3次元形状を実時間で計測する3次元形状計測器を有し、補正応答関数生成部は、被測定体の存在しうる空間における所定間隔の格子点の個々について放射線検出部の応答関数を有し、被測定体が検査領域内に占有する空間に包含される個々の格子点についての応答関数を積算することにより、補正応答関数を生成するものである。 Radiation measuring device according to this invention detects radiation emitted from the object to be measured, and the radiation detecting section for outputting a pulse-height distribution in accordance with the output pulse, the object to be measured for obtaining the shape information of the object to be measured A correction response function that calculates a correction value of a response function corresponding to the shape information of the measurement object as a correction response function based on the information acquisition unit and the shape information of the measurement object acquired by the measurement object information acquisition unit comprising a generator, a, the output pulse height distribution from the radiation detecting unit, by performing the inverse problem calculation using the corrected response function calculated by the corrected response function generating unit, is emitted from the object to be measured The measurement object information acquisition unit includes a three-dimensional shape measuring instrument that measures in real time the three-dimensional shape of at least one measurement object in the examination region. And complement The response function generation unit has a response function of the radiation detection unit for each of the lattice points at a predetermined interval in the space where the measured object can exist, and each of the individual included in the space occupied by the measured object in the examination region The corrected response function is generated by integrating the response functions for the lattice points .
また、この発明に係る放射線測定方法は、被測定体から放射される放射線を検出し、出力パルスに応じたパルス波高分布を出力する放射線検出ステップと、被測定体の形状情報を取得する被測定体情報取得ステップと、被測定体情報取得ステップで取得された被測定体の形状情報に基づいて、被測定体の形状情報に応じた応答関数の補正値を、補正応答関数として演算する補正応答関数生成ステップと、を備え、放射線検出ステップから出力されたパルス波高分布に対して、補正応答関数生成ステップで演算された補正応答関数を用いた逆問題演算を行うことにより、被測定体から放射される放射線のエネルギースペクトルを演算する放射線測定方法であって、被測定体情報取得ステップは、検査領域内の少なくとも1つの被測定体の3次元形状を実時間で計測する3次元形状計測器を用い、補正応答関数生成ステップは、被測定体の存在しうる空間における所定間隔の格子点の個々について放射線検出ステップで用いられる応答関数に基づいて、被測定体が検査領域内に占有する空間に包含される個々の格子点についての応答関数を積算することにより、補正応答関数を生成するものである。 Further, radiation measurement method according to this invention detects radiation emitted from the object to be measured, and the radiation detecting step of outputting a pulse-height distribution corresponding to the output pulse, obtains the shape information of the object to be measured to be Based on the measured object information acquisition step and the measured object information acquisition step, the correction for calculating the correction value of the response function corresponding to the measured object shape information as a corrected response function includes a response function generating step, the relative pulse height distribution outputted from the radiation detecting step, by performing the inverse problem calculation using the corrected response function calculated by the corrected response function generating step, the object to be measured a radiation measuring method for calculating an energy spectrum of the radiation emitted, measured body information acquiring step, three-dimensional least one object to be measured in the inspection area The correction response function generation step uses a three-dimensional shape measuring instrument that measures the shape in real time, and the correction response function generation step is based on the response function used in the radiation detection step for each of the lattice points at predetermined intervals in the space where the measured object can exist. The correction response function is generated by integrating the response functions for the individual lattice points included in the space occupied by the measured object in the inspection region .
この発明に係る放射線測定装置および放射線測定方法によれば、補正応答関数生成部(ステップ)は、被測定体の形状情報を取得する被測定体情報取得部(ステップ)で取得された被測定体の形状情報に基づいて、被測定体の形状情報に応じた応答関数の補正値を、補正応答関数として演算し、被測定体から放射される放射線を検出し、出力パルスに応じたパルス波高分布を出力する放射線検出部(ステップ)から出力されたパルス波高分布に対して、補正応答関数生成部(ステップ)で演算された補正応答関数を用いた逆問題演算を行うことにより、被測定体から放射される放射線のエネルギースペクトルを演算する。
このとき、被測定体情報取得部は、検査領域内の少なくとも1つの被測定体の3次元形状を実時間で計測する3次元形状計測器を有し、補正応答関数生成部は、被測定体の存在しうる空間における所定間隔の格子点の個々について放射線検出部の応答関数を有し、被測定体が検査領域内に占有する空間に包含される個々の格子点についての応答関数を積算することにより、補正応答関数を生成する。
また、被測定体情報取得ステップは、検査領域内の少なくとも1つの被測定体の3次元形状を実時間で計測する3次元形状計測器を用い、補正応答関数生成ステップは、被測定体の存在しうる空間における所定間隔の格子点の個々について放射線検出ステップで用いられる応答関数に基づいて、被測定体が検査領域内に占有する空間に包含される個々の格子点についての応答関数を積算することにより、補正応答関数を生成する。
そのため、試料毎に形状が不定形の被測定体に対しても、最適な検出器応答関数を算出し、この応答関数を用いて逆問題演算を行うことにより、被測定体を破壊することなく全数検査を行うとともに、被測定体から放射される放射線のエネルギースペクトルを高精度に測定することができる。
According to radiation measurement apparatus and the radiation beam measuring method according to the present invention, the correction response function generating unit (step) is to be obtained by the measured object information acquisition unit that acquires shape information of the object to be measured (step) based on the shape information of the measured sample, the correction value of the response function corresponding to the shape information of the object to be measured, calculated as a correction response function, to detect radiation emitted from the object to be measured, corresponding to the output pulse pulse By performing inverse problem calculation using the corrected response function calculated by the correction response function generation unit (step) on the pulse wave height distribution output from the radiation detection unit (step) that outputs the pulse height distribution, the measurement target Calculate the energy spectrum of radiation emitted from the body.
At this time, the measured object information acquisition unit includes a three-dimensional shape measuring instrument that measures the three-dimensional shape of at least one measured object in the inspection region in real time, and the correction response function generation unit includes the measured object A response function of the radiation detection unit is provided for each of the lattice points at a predetermined interval in the space where the object can exist, and the response functions for the individual lattice points included in the space occupied by the measured object in the inspection region are integrated. Thus, a corrected response function is generated.
The measured object information acquisition step uses a three-dimensional shape measuring instrument that measures the three-dimensional shape of at least one measured object in the inspection region in real time, and the correction response function generation step includes the presence of the measured object. Based on the response function used in the radiation detection step for each of the lattice points at predetermined intervals in the possible space, the response functions for the individual lattice points included in the space occupied by the measured object in the examination region are integrated. Thus, a corrected response function is generated.
For this reason, the optimal detector response function is calculated even for an object with an indefinite shape for each sample, and the inverse function is calculated using this response function without destroying the object. While performing 100% inspection, it is possible to measure the energy spectrum of the radiation emitted from the measurement object with high accuracy.
以下、この発明に係る放射能分析装置および放射能分析方法の好適な実施の形態につき図面を用いて説明するが、各図において同一、または相当する部分については、同一符号を付して説明する。 Hereinafter, preferred embodiments of a radioactivity analyzer and a radioactivity analysis method according to the present invention will be described with reference to the drawings. In the drawings, the same or corresponding parts will be described with the same reference numerals. .
実施の形態1.
図1は、この発明の実施の形態1に係る放射能分析装置を示すブロック図である。図1において、この放射能分析装置は、被測定体1から放射される放射線を検出する放射線検出部10、被測定体情報取得部20、補正応答関数生成部30、放射能値演算部40および表示部50を備えている。
FIG. 1 is a block diagram showing a radioactivity analyzer according to
放射線検出部10は、例えば検出器11、増幅器12および波高分析器13を有し、被測定体1から放射されるγ線等の放射線を検出し、その出力パルスの波高に対してパルス数をカウントする。
The
検出器11は、入射放射線のエネルギーに対し、それに応じた確率密度分布に従ってパルス出力を発生する。例えば、NaIシンチレータを用いたγ線計測の場合であれば、全吸収ピークやコンプトン散乱による連続波高分布等が発生する。この分布関数を用いて逆演算を行い、高精度にエネルギーを分析する方法が逆問題演算(アンフォールディング)であるが、この応答関数は、検出器11に対する放射線源の位置により変化する。
The
図2は、この発明の実施の形態1に係る放射能分析装置における応答関数分布を示す説明図である。図2において、検出器11に対し、被測定体1の設置される領域を格子状にL、M、N分割した検出空間分割格子2を定義し、応答関数をRで表す。
FIG. 2 is an explanatory diagram showing a response function distribution in the radioactivity analyzer according to
このとき、検出空間分割格子2の格子点(i、j、k)に、放射線のエネルギーEの点状の放射線源を置いて得られる放射線検出部10の出力応答関数R(E、i、j、k)と、格子点(l、m、n)の出力応答関数R(E、l、m、n)とは、互いに異なる分布となる。
At this time, the output response function R (E, i, j) of the
図3は、この発明の実施の形態1に係る放射能分析装置の補正応答関数生成部30を示すブロック図である。図3において、補正応答関数生成部30は、検出器応答関数記憶部31、被測定体分布演算部32、積分演算処理部33、補正応答関数出力部34、密度分布情報記憶部35(後述する)および放射性物質含有率分布情報記憶部36(後述する)を有している。
FIG. 3 is a block diagram showing the correction response
図2に示された放射線検出部10で得られた出力応答関数は、検出器応答関数記憶部31のVoxelデータとして記憶される。一方、被測定体情報取得部20では、図4に示されるように、3次元形状計測器21により被測定体1の形状情報が取得され、被測定体分布演算部32に出力される。
The output response function obtained by the
被測定体分布演算部32は、入力された被測定体1の形状情報に基づいて、格子点(i、j、k)に重み付けを行い、重み付けデータW(i、j、k)を算出する。最も簡単には、格子点が被測定体情報取得部20で取得された被測定体1の内部にある場合には1を、外部にある場合には0を入力することが考えられるが、被測定体1の境界において、隣接格子点との境界内に存在する被測定体1の占有率により、荷重するようにしてもよい。
The measured object
次に、積分演算処理部33は、被測定体分布演算部32で算出された重み付けデータW(i、j、k)と、検出器応答関数記憶部31に記憶されたVoxelデータとを合わせて次式(1)で表される演算を行い、補正応答関数出力部34は、式(1)で演算された被測定体1の応答関数R(E)を出力する。
Next, the integral
放射能値演算部40は、放射線検出部10から出力された被測定体1からの放射線によるパルス波高分布に対して、補正応答関数出力部34から出力された被測定体1の応答関数R(E)を用いて、例えば逐次近似法等を用いた逆問題演算を行い、エネルギースペクトルを算出して、表示部50にその結果を出力する。
The radioactivity value calculation unit 40 responds to the pulse wave height distribution due to the radiation from the
なお、表示部50は、エネルギースペクトルそのものを表示するものに限定されず、ベクレル値に変換して表示してもよい。また、スクリーニング等を行うような場合には、一定量の放射能の有無を判定表示するだけであってもよいし、食品の選別等、被測定体1の処理結果として表示されるものであってもよい。
In addition, the
また、被測定体情報取得部20として高速のコンピュータを用いることにより、3次元形状計測器21により得られる被測定体1の分布を実時間で(リアルタイムに)演算することができ、例えばコンベア等で食品を連続で検査する全数検査システムとして適用することができる。
Further, by using a high-speed computer as the measured object
また、複数の被測定体1が検出領域に存在する場合であっても、3次元形状計測器21により得られるデータから補正応答関数を正確に算出することができ、正確な検査を行うことができる。
Even when a plurality of objects to be measured 1 exist in the detection region, the correction response function can be accurately calculated from the data obtained by the three-dimensional
以上のように、実施の形態1によれば、補正応答関数生成部は、被測定体の形状情報を取得する被測定体情報取得部で取得された被測定体の形状情報に基づいて、被測定体の形状情報に応じた応答関数の補正値を、補正応答関数として演算し、放射能値演算部は、被測定体から放射される放射線を検出し、出力パルスに応じたパルス波高分布を出力する放射線検出部から出力されたパルス波高分布に対して、補正応答関数生成部で演算された補正応答関数を用いた逆問題演算を行うことにより、被測定体に含まれる放射能量を演算する。
すなわち、点状線源に対するエネルギースペクトルの応答関数を持つとともに、被測定体の3次元形状情報および線源の分布情報を入力する事により、応答関数から逆問題演算に適用する応答関数の補正値が演算される。
また、被測定体情報取得部は、検査領域内の少なくとも1つの被測定体の3次元形状を実時間で計測する3次元形状計測器を有し、補正応答関数生成部は、被測定体の存在しうる空間における所定間隔の格子点の個々について放射線検出部の応答関数を有し、被測定体が検査領域内に占有する空間に包含される個々の格子点についての応答関数を積算することにより、補正応答関数を生成する。
そのため、試料毎に形状が不定形の被測定体に対しても、最適な検出器応答関数を算出し、この応答関数を用いて逆問題演算を行うことにより、被測定体を破壊することなく全数検査を行うとともに、被測定体に含まれる放射能量を高精度に測定することができる。
また、不定形の連続して流れる被測定体に対して、被測定体を粉砕することなく全数検査行うことができる。
As described above, according to the first embodiment, the correction response function generation unit is based on the shape information of the measurement target acquired by the measurement target information acquisition unit that acquires the shape information of the measurement target. The correction value of the response function corresponding to the shape information of the measurement object is calculated as a correction response function, and the radioactivity value calculation unit detects the radiation emitted from the measurement object, and calculates the pulse height distribution corresponding to the output pulse. By calculating the inverse problem using the correction response function calculated by the correction response function generation unit for the pulse wave height distribution output from the output radiation detection unit, the amount of radioactivity contained in the measured object is calculated. .
That is, it has a response function of the energy spectrum for the point source, and by inputting the three-dimensional shape information of the object to be measured and the distribution information of the source, the correction value of the response function applied to the inverse problem calculation from the response function Is calculated.
The measured object information acquisition unit includes a three-dimensional shape measuring instrument that measures a three-dimensional shape of at least one measured object in the inspection region in real time, and the correction response function generation unit A response function of the radiation detection unit is provided for each lattice point at a predetermined interval in a space that can exist, and the response function for each lattice point included in the space occupied by the measured object in the examination region is integrated. To generate a corrected response function.
For this reason, the optimal detector response function is calculated even for an object with an indefinite shape for each sample, and the inverse function is calculated using this response function without destroying the object. While performing 100% inspection, the amount of radioactivity contained in the measurement object can be measured with high accuracy.
Further, 100% inspection can be performed on an object to be measured that flows continuously indefinitely without pulverizing the object to be measured.
実施の形態2.
例えばピーマンのように、食品には、外形に対して空洞等が存在する等、密度が均一でない場合が多い。そこで、図5に示されるように、被測定体1に空洞が存在する場合には、それを除外した領域での積分演算により、応答関数R(E)を算出する必要がある。
For example, like peppers, foods often have non-uniform densities, such as cavities with respect to the outer shape. Therefore, as shown in FIG. 5, when a cavity exists in the device under
そこで、図3に示されるように、この発明の実施の形態2に係る放射能分析装置は、補正応答関数生成部30において、密度分布情報記憶部35をさらに有している。密度分布情報記憶部35は、被測定体情報として、マンマシンインターフェースや画像認識等により、被測定体の種別を密度分布情報とともに記憶している。
Therefore, as shown in FIG. 3, the radioactivity analyzer according to
被測定体分布演算部32は、密度分布情報記憶部35に記憶された被測定体の種別に対応する密度分布情報を用いて、被測定体情報取得部20により得られる外形形状から、内部状態を推定して被測定体分布データW(i、j、k)を演算する。また、積分演算処理部33は、この被測定体分布データW(i、j、k)と検出器応答関数記憶部31に記憶されたVoxelデータとを合わせて、補正応答関数を演算する。
The measured object
以上のように、実施の形態2によれば、被測定体情報取得部は、検査領域内の少なくとも1つの被測定体の3次元形状を実時間で計測する3次元形状計測器を有し、補正応答関数生成部は、被測定体の密度分布情報を有するとともに、被測定体が検査領域内に占有する空間に包含される格子点について、格子点における密度を乗算した応答関数を積算することにより、補正応答関数を生成する。
そのため、密度分布を有する被測定体に対して、被測定体を粉砕することなく全数検査行うことができる。
As described above, according to the second embodiment, the measured object information acquisition unit includes a three-dimensional shape measuring instrument that measures the three-dimensional shape of at least one measured object in the inspection region in real time. The correction response function generation unit accumulates the response function obtained by multiplying the density at the lattice points with respect to the lattice points included in the space occupied by the measured object in the inspection region, as well as having the density distribution information of the object to be measured. To generate a corrected response function.
Therefore, 100% inspection can be performed on a measurement object having a density distribution without crushing the measurement object.
実施の形態3.
例えば生体では、放射性物質を取り込み蓄積する部位が局在する場合があり、たとえ密度が一様であったとしても、図6に示されるように、第1放射性物質濃度領域1aおよび第2放射性物質濃度領域1b等、複数の放射性物質濃度の領域に分割される。ここで、放射性物質を蓄積する部位は、生物種で一定であり、種別および外形が分かれば、領域を特定することができる。
For example, in a living body, there are cases where a site for taking in and accumulating radioactive substances is localized, and even if the density is uniform, as shown in FIG. 6, the first radioactive substance concentration region 1a and the second radioactive substance It is divided into a plurality of radioactive substance concentration regions such as the concentration region 1b. Here, the site for accumulating radioactive substances is constant for each species, and the region can be specified if the type and outline are known.
そこで、図3に示されるように、この発明の実施の形態3に係る放射能分析装置は、補正応答関数生成部30において、放射性物質含有率分布情報記憶部36をさらに有している。放射性物質含有率分布情報記憶部36は、被測定体情報として、外形形状に対する相対的な領域分割データと含有率としての重み値を記憶している。
Therefore, as shown in FIG. 3, the radioactivity analyzer according to
被測定体分布演算部32は、放射性物質含有率分布情報記憶部36に記憶された外形形状に対する相対的な領域分割データと含有率としての重み値を用いて、被測定体情報取得部20により得られる外形形状から、被測定体分布データW(i、j、k)を演算する。また、積分演算処理部33は、この被測定体分布データW(i、j、k)と検出器応答関数記憶部31に記憶されたVoxelデータとを合わせて、補正応答関数を演算する。
The to-be-measured object
なお、この発明の実施の形態3に係る放射能分析装置において、放射性物質含有率分布情報記憶部36は、上述した実施の形態2における密度分布情報記憶部35と併用してもよく、生体に密度分布および蓄積部位双方が存在する場合に、より正確な補正応答関数を生成することができる。
In the radioactivity analyzer according to
以上のように、実施の形態3によれば、被測定体情報取得部は、検査領域内の少なくとも1つの被測定体の3次元形状を実時間で計測する3次元形状計測器を有し、補正応答関数生成部は、被測定体の放射性物質の含有率分布情報を有するとともに、被測定体が検査領域内に占有する空間に包含される格子点について、格子点における放射性物質の含有率を乗算した応答関数を積算することにより、補正応答関数を生成する。
そのため、放射性物質の含有量に分布を有する被測定体に対して、被測定体を粉砕することなく全数検査行うことができる。
As described above, according to the third embodiment, the measured object information acquisition unit includes the three-dimensional shape measuring instrument that measures the three-dimensional shape of at least one measured object in the inspection region in real time. The correction response function generation unit has the content distribution information of the radioactive material of the measured object, and the content of the radioactive material at the lattice point for the lattice points included in the space occupied by the measured object in the inspection region. A corrected response function is generated by integrating the multiplied response functions.
Therefore, 100% inspection can be performed on the measurement object having a distribution in the content of the radioactive substance without crushing the measurement object.
1 被測定体、1a 第1放射性物質濃度領域、1b 第2放射性物質濃度領域、2 検出空間分割格子、10 放射線検出部、11 検出器、12 増幅器、13 波高分析器、20 被測定体情報取得部、21 3次元形状計測器、30 補正応答関数生成部、31 検出器応答関数記憶部、32 被測定体分布演算部、33 積分演算処理部、34 補正応答関数出力部、35 密度分布情報記憶部、36 放射性物質含有率分布情報記憶部、40 放射能値演算部、50 表示部。
DESCRIPTION OF
Claims (6)
前記被測定体の形状情報を取得する被測定体情報取得部と、
前記被測定体情報取得部で取得された前記被測定体の形状情報に基づいて、前記被測定体の形状情報に応じた応答関数の補正値を、補正応答関数として演算する補正応答関数生成部と、を備え、
前記放射線検出部から出力されたパルス波高分布に対して、前記補正応答関数生成部で演算された前記補正応答関数を用いた逆問題演算を行うことにより、前記被測定体から放射される放射線のエネルギースペクトルを演算する放射線測定装置であって、
前記被測定体情報取得部は、検査領域内の少なくとも1つの前記被測定体の3次元形状を実時間で計測する3次元形状計測器を有し、
前記補正応答関数生成部は、前記被測定体の存在しうる空間における所定間隔の格子点の個々について前記放射線検出部の応答関数を有し、前記被測定体が前記検査領域内に占有する空間に包含される個々の前記格子点についての応答関数を積算することにより、前記補正応答関数を生成する
放射線測定装置。 A radiation detector that detects radiation emitted from the measurement object and outputs a pulse height distribution according to the output pulse;
A measured object information acquisition unit that acquires shape information of the measured object;
A correction response function generation unit that calculates a correction value of a response function according to the shape information of the measurement object as a correction response function based on the shape information of the measurement object acquired by the measurement object information acquisition unit and, with a,
By performing inverse problem calculation using the correction response function calculated by the correction response function generation unit on the pulse wave height distribution output from the radiation detection unit, the radiation emitted from the measurement object is measured . a radiation measuring device for calculating the energy spectrum,
The measured object information acquisition unit includes a three-dimensional shape measuring instrument that measures a three-dimensional shape of at least one of the measured objects in an inspection region in real time,
The correction response function generation unit has a response function of the radiation detection unit for each of lattice points at a predetermined interval in a space where the measured object can exist, and a space occupied by the measured object in the inspection region Generating the corrected response function by integrating the response functions for each of the grid points included in
Radiation measurement device .
前記被測定体の形状情報を取得する被測定体情報取得部と、
前記被測定体情報取得部で取得された前記被測定体の形状情報に基づいて、前記被測定体の形状情報に応じた応答関数の補正値を、補正応答関数として演算する補正応答関数生成部と、を備え、
前記放射線検出部から出力されたパルス波高分布に対して、前記補正応答関数生成部で演算された前記補正応答関数を用いた逆問題演算を行うことにより、前記被測定体から放射される放射線のエネルギースペクトルを演算する放射線測定装置であって、
前記被測定体情報取得部は、検査領域内の少なくとも1つの前記被測定体の3次元形状を実時間で計測する3次元形状計測器を有し、
前記補正応答関数生成部は、前記被測定体の密度分布情報を有するとともに、前記被測定体が前記検査領域内に占有する空間に包含される格子点について、前記格子点における密度を乗算した応答関数を積算することにより、前記補正応答関数を生成する
放射線測定装置。 A radiation detector that detects radiation emitted from the measurement object and outputs a pulse height distribution according to the output pulse;
A measured object information acquisition unit that acquires shape information of the measured object;
A correction response function generation unit that calculates a correction value of a response function according to the shape information of the measurement object as a correction response function based on the shape information of the measurement object acquired by the measurement object information acquisition unit and, with a,
By performing inverse problem calculation using the correction response function calculated by the correction response function generation unit on the pulse wave height distribution output from the radiation detection unit, the radiation emitted from the measurement object is measured . a radiation measuring device for calculating the energy spectrum,
The measured object information acquisition unit includes a three-dimensional shape measuring instrument that measures a three-dimensional shape of at least one of the measured objects in an inspection region in real time,
The correction response function generation unit has a density distribution information of the object to be measured, and a response obtained by multiplying the lattice points included in the space occupied by the object to be measured in the inspection region by the density at the lattice points. Generate the corrected response function by integrating the function
Radiation measurement device .
前記被測定体の形状情報を取得する被測定体情報取得部と、
前記被測定体情報取得部で取得された前記被測定体の形状情報に基づいて、前記被測定体の形状情報に応じた応答関数の補正値を、補正応答関数として演算する補正応答関数生成部と、を備え、
前記放射線検出部から出力されたパルス波高分布に対して、前記補正応答関数生成部で演算された前記補正応答関数を用いた逆問題演算を行うことにより、前記被測定体から放射される放射線のエネルギースペクトルを演算する放射線測定装置であって、
前記被測定体情報取得部は、検査領域内の少なくとも1つの前記被測定体の3次元形状を実時間で計測する3次元形状計測器を有し、
前記補正応答関数生成部は、前記被測定体の放射性物質の含有率分布情報を有するとともに、前記被測定体が前記検査領域内に占有する空間に包含される格子点について、前記格子点における放射性物質の含有率を乗算した応答関数を積算することにより、前記補正応答関数を生成する
放射線測定装置。 A radiation detector that detects radiation emitted from the measurement object and outputs a pulse height distribution according to the output pulse;
A measured object information acquisition unit that acquires shape information of the measured object;
A correction response function generation unit that calculates a correction value of a response function according to the shape information of the measurement object as a correction response function based on the shape information of the measurement object acquired by the measurement object information acquisition unit and, with a,
By performing inverse problem calculation using the correction response function calculated by the correction response function generation unit on the pulse wave height distribution output from the radiation detection unit, the radiation emitted from the measurement object is measured . a radiation measuring device for calculating the energy spectrum,
The measured object information acquisition unit includes a three-dimensional shape measuring instrument that measures a three-dimensional shape of at least one of the measured objects in an inspection region in real time,
The correction response function generation unit has content distribution information of the radioactive substance of the object to be measured, and for the lattice points included in the space occupied by the object to be measured in the inspection region, the radioactivity at the lattice points The corrected response function is generated by integrating the response function multiplied by the substance content rate.
Radiation measurement device .
前記被測定体の形状情報を取得する被測定体情報取得ステップと、
前記被測定体情報取得ステップで取得された前記被測定体の形状情報に基づいて、前記被測定体の形状情報に応じた応答関数の補正値を、補正応答関数として演算する補正応答関数生成ステップと、を備え、
前記放射線検出ステップから出力されたパルス波高分布に対して、前記補正応答関数生成ステップで演算された前記補正応答関数を用いた逆問題演算を行うことにより、前記被測定体から放射される放射線のエネルギースペクトルを演算する放射線測定方法であって、
前記被測定体情報取得ステップは、検査領域内の少なくとも1つの前記被測定体の3次元形状を実時間で計測する3次元形状計測器を用い、
前記補正応答関数生成ステップは、前記被測定体の存在しうる空間における所定間隔の格子点の個々について前記放射線検出ステップで用いられる応答関数に基づいて、前記被測定体が前記検査領域内に占有する空間に包含される個々の前記格子点についての応答関数を積算することにより、前記補正応答関数を生成する
放射線測定方法。 A radiation detection step of detecting radiation emitted from the measurement object and outputting a pulse height distribution according to the output pulse;
Measurement object information acquisition step for acquiring shape information of the measurement object;
A correction response function generating step for calculating a correction value of a response function according to the shape information of the measurement object as a correction response function based on the shape information of the measurement object acquired in the measurement object information acquisition step and, with a,
By performing inverse problem calculation using the correction response function calculated in the correction response function generation step on the pulse wave height distribution output from the radiation detection step, the radiation emitted from the measured object a radiation measuring method for calculating the energy spectrum,
The measurement object information acquisition step uses a three-dimensional shape measuring instrument that measures a three-dimensional shape of at least one of the measurement objects in an inspection region in real time,
In the correction response function generation step, the measurement object occupies the inspection area based on the response function used in the radiation detection step for each of the lattice points at predetermined intervals in the space where the measurement object can exist. The corrected response function is generated by integrating response functions for the individual lattice points included in the space to be
Radiation measurement method .
前記被測定体の形状情報を取得する被測定体情報取得ステップと、
前記被測定体情報取得ステップで取得された前記被測定体の形状情報に基づいて、前記被測定体の形状情報に応じた応答関数の補正値を、補正応答関数として演算する補正応答関数生成ステップと、を備え、
前記放射線検出ステップから出力されたパルス波高分布に対して、前記補正応答関数生成ステップで演算された前記補正応答関数を用いた逆問題演算を行うことにより、前記被測定体から放射される放射線のエネルギースペクトルを演算する放射線測定方法であって、
前記被測定体情報取得ステップは、検査領域内の少なくとも1つの前記被測定体の3次元形状を実時間で計測する3次元形状計測器を用い、
前記補正応答関数生成ステップは、前記被測定体の密度分布情報に基づいて、前記被測定体が前記検査領域内に占有する空間に包含される格子点について、前記格子点における密度を乗算した応答関数を積算することにより、前記補正応答関数を生成する
放射線測定方法。 A radiation detection step of detecting radiation emitted from the measurement object and outputting a pulse height distribution according to the output pulse;
Measurement object information acquisition step for acquiring shape information of the measurement object;
A correction response function generating step for calculating a correction value of a response function according to the shape information of the measurement object as a correction response function based on the shape information of the measurement object acquired in the measurement object information acquisition step and, with a,
By performing inverse problem calculation using the correction response function calculated in the correction response function generation step on the pulse wave height distribution output from the radiation detection step, the radiation emitted from the measured object a radiation measuring method for calculating the energy spectrum,
The measurement object information acquisition step uses a three-dimensional shape measuring instrument that measures a three-dimensional shape of at least one of the measurement objects in an inspection region in real time,
In the correction response function generation step, a response obtained by multiplying a density at the grid point with respect to a grid point included in a space occupied by the measured object in the inspection region based on the density distribution information of the measured object. Generate the corrected response function by integrating the function
Radiation measurement method .
前記被測定体の形状情報を取得する被測定体情報取得ステップと、
前記被測定体情報取得ステップで取得された前記被測定体の形状情報に基づいて、前記被測定体の形状情報に応じた応答関数の補正値を、補正応答関数として演算する補正応答関数生成ステップと、を備え、
前記放射線検出ステップから出力されたパルス波高分布に対して、前記補正応答関数生成ステップで演算された前記補正応答関数を用いた逆問題演算を行うことにより、前記被測定体から放射される放射線のエネルギースペクトルを演算する放射線測定方法であって、
前記被測定体情報取得ステップは、検査領域内の少なくとも1つの前記被測定体の3次元形状を実時間で計測する3次元形状計測器を用い、
前記補正応答関数生成ステップは、前記被測定体の放射性物質の含有率分布情報に基づいて、前記被測定体が前記検査領域内に占有する空間に包含される格子点について、前記格子点における放射性物質の含有率を乗算した応答関数を積算することにより、前記補正応答関数を生成する
放射線測定装置。 A radiation detection step of detecting radiation emitted from the measurement object and outputting a pulse height distribution according to the output pulse;
Measurement object information acquisition step for acquiring shape information of the measurement object;
A correction response function generating step for calculating a correction value of a response function according to the shape information of the measurement object as a correction response function based on the shape information of the measurement object acquired in the measurement object information acquisition step and, with a,
By performing inverse problem calculation using the correction response function calculated in the correction response function generation step on the pulse wave height distribution output from the radiation detection step, the radiation emitted from the measured object a radiation measuring method for calculating the energy spectrum,
The measurement object information acquisition step uses a three-dimensional shape measuring instrument that measures a three-dimensional shape of at least one of the measurement objects in an inspection region in real time,
In the correction response function generation step, based on the content distribution information of the radioactive substance of the object to be measured, for the lattice points included in the space that the object to be measured occupies in the inspection area, The corrected response function is generated by integrating the response function multiplied by the substance content rate.
Radiation measurement device .
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013163859A JP6037968B2 (en) | 2013-08-07 | 2013-08-07 | Radiation measurement apparatus and radiation measurement method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013163859A JP6037968B2 (en) | 2013-08-07 | 2013-08-07 | Radiation measurement apparatus and radiation measurement method |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2015031675A JP2015031675A (en) | 2015-02-16 |
JP6037968B2 true JP6037968B2 (en) | 2016-12-07 |
Family
ID=52517094
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2013163859A Active JP6037968B2 (en) | 2013-08-07 | 2013-08-07 | Radiation measurement apparatus and radiation measurement method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6037968B2 (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6709490B2 (en) * | 2016-03-07 | 2020-06-17 | 株式会社北川鉄工所 | Radioactivity concentration measuring device and radioactivity concentration measuring method |
JP6925842B2 (en) * | 2017-04-04 | 2021-08-25 | 株式会社東芝 | Radioactivity measuring device and radioactivity measuring method |
JP2020027079A (en) * | 2018-08-17 | 2020-02-20 | 三菱電機株式会社 | Radioactivity analysis device |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3616036B2 (en) * | 2001-06-20 | 2005-02-02 | 三菱重工業株式会社 | Article radioactivity detection apparatus and article radioactivity detection system |
JP4576108B2 (en) * | 2003-11-07 | 2010-11-04 | 三菱重工業株式会社 | Radioactivity measurement method, radioactivity measurement program, and radioactivity measurement apparatus |
WO2006131172A1 (en) * | 2005-06-10 | 2006-12-14 | Gesellschaft für Schwerionenforschung mbH | Method for the spectrometric photon dosimetry for x-ray and gamma radiation |
JP2007024795A (en) * | 2005-07-20 | 2007-02-01 | Horiba Ltd | Configuration measuring device |
JP2008139094A (en) * | 2006-11-30 | 2008-06-19 | Toshiba Corp | Radioactivity measuring method and instrument |
WO2013031897A1 (en) * | 2011-08-31 | 2013-03-07 | 富士電機株式会社 | Food product inspection system |
JP2013061254A (en) * | 2011-09-13 | 2013-04-04 | Kimoto Denshi Kogyo Kk | Radioactive suspended particulate matter measuring instrument and radioactive suspended particulate matter measuring method |
JP3179617U (en) * | 2012-08-29 | 2012-11-08 | 株式会社島津製作所 | Food radioactivity screening inspection equipment |
-
2013
- 2013-08-07 JP JP2013163859A patent/JP6037968B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
JP2015031675A (en) | 2015-02-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5340524B2 (en) | Radiation detector and radiation detection method | |
CN106353828B (en) | The method and apparatus that checked property body weight is estimated in safe examination system | |
JP4903654B2 (en) | Radiation detection apparatus and efficiency calibration method thereof | |
KR100930681B1 (en) | Monitoring Method of Environmental Radioactivity Monitoring System for Dose Rate Measurement by Radionuclide | |
JP5988890B2 (en) | Radioactivity analyzer and radioactivity analysis method | |
RU2019143155A (en) | METHOD AND DEVICE FOR MULTI-ELEMENT ANALYSIS BASED ON NEUTRON ACTIVATION, AS WELL AS APPLICATION | |
De Mesquita et al. | Industrial tomography using three different gamma ray | |
Kim et al. | Industrial gamma-ray tomographic scan method for large scale industrial plants | |
JP6037968B2 (en) | Radiation measurement apparatus and radiation measurement method | |
JP6833605B2 (en) | Surface contamination density distribution calculation device, method and program | |
Nazemi et al. | Proposing an intelligent approach for measuring the thickness of metal sheets independent of alloy type | |
KR102394944B1 (en) | Apparatus for measuring and amending activity concentration density of radioactive material and operation calibration methode thereof | |
Johansen | Gamma-ray tomography | |
JP7000264B2 (en) | Radioactivity concentration evaluation system and radioactivity concentration evaluation method | |
JP6674836B2 (en) | Waste measuring device and waste measuring method | |
JP6925842B2 (en) | Radioactivity measuring device and radioactivity measuring method | |
KR102249120B1 (en) | Analysis program and analysis system for radioactive contamination in and outside the body of the subject | |
JP4994225B2 (en) | Two-parameter spectrum processing method and apparatus | |
JP5926362B1 (en) | Radioactivity concentration measuring apparatus and radioactivity concentration measuring method | |
Miller et al. | Combining radiography and passive measurements for radiological threat localization in cargo | |
Vella et al. | An artificial X-ray wire test emitter and calculations on the resolution and field of view of X-ray pinhole optics by simulation | |
Conlin et al. | On using code emulators and Monte Carlo estimation to predict assembly attributes of spent fuel assemblies for safeguards applications | |
Stepanov et al. | Application of gamma-ray imager for non-destructive testing | |
JP7183206B2 (en) | Radioactivity inspection device | |
RU2730107C2 (en) | Device for in situ analysis of radioactive wastes containing isotope chlorine-36 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20150918 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20160705 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20160902 |
|
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: 20161004 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20161101 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6037968 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |