JP7477890B2 - γ線計測方法およびγ線計測装置 - Google Patents

Description

本発明は、γ線の線量とエネルギースペクトルとを、同時に計測して表示するγ線計測方法、および、前記γ線計測方法が組み込まれたγ線計測装置に関する。

原子炉や加速器施設などの核施設では、放射線の線量をモニターして、放射線作業従事者らの安全を確保することが法律で定められており、そのために、電離箱やＧＭ計数管などを用いたサーベイメータや個人線量計などの放射線計測装置を用いて、放射線の線量を計測することが行われている（非特許文献１）。

しかしながら、これらのサーベイメータや個人線量計は、一般的にかなりの重量がある。また、いずれも、表示は線量のみであり、しかも、その数値の評価は^１３７Ｃｓの一点校正によるものであり、γ線のような高エネルギーの放射線に対してまでの精度は保証されていなかった。

即ち、従来の一般的なサーベイメータや個人線量計では、^１３７Ｃｓγ線標準線源（６６２ｋｅＶ）での値に基づいて計測値の絶対値を規格化することにより、線量値の一点校正を行っている。このため、線量計にもよるが、１．５～２．０ＭｅＶ程度までであれば、数１０％以内の誤差に収めることができるが、前記した^１３７Ｃｓγ線標準線源から離れるに従って、計測の正確さが欠けてくることが避けられない。

実際のγ線エネルギーは、中性子がない状態（原子炉が停止している状態や、加速器が運転されていない状態など）では３ＭｅＶ以下であり、中性子がある状態（加速器が運転状態）では最大で１０ＭｅＶと大きい。このため、１０ＭｅＶはおろか、３ＭｅＶまでの高エネルギーにも対応できていない従来の一般的なサーベイメータや個人線量計においては、大きめの誤差を乗せた数値を保守的な数値として評価することで、過小評価となることにより実用上において発生する不具合を回避している。しかしながら、本発明者が確認したところ、従来の一般的なサーベイメータや個人線量計では、中性子が存在している状態において、理論値（正解値）の倍以上の数値がγ線の線量として表示されていることが分かった。

また、γ線エネルギーは、そのスペクトルに分布があり、形も一定していないため、線量をリアルタイムで知るには、そのエネルギースペクトルをリアルタイムに知る必要があるが、現在のところ、サーベイメータや個人線量計によりγ線のエネルギースペクトルをリアルタイムに知る技術は、まだ提案されていない。

このため、放射線作業従事者らは、作業を第一優先にせざるを得ないこととも相俟って、重くて身につけやすいとは言えない現状から、サーベイメータや個人線量計を使用することには抵抗感があり、放射線作業従事者らの放射線被曝への意識が低くなり、自身の被曝管理が疎かになりがちであった。

「放射線による健康影響等に関する統一的な基礎資料 平成２９年度版 第２章 放射線による被ばく（環境省）」、ｈｔｔｐｓ：／／ｗｗｗ．ｅｎｖ．ｇｏ．ｊｐ／ｃｈｅｍｉ／ｒｈｍ／ｋｉｓｏｓｈｉｒｙｏ／ｐｄｆ＿ｈ２９／２０１７ｔｋ１ｓ０２．ｐｄｆ

本発明は、上記の問題に鑑みて、線量だけでなくγ線のエネルギースペクトルをも、その場でリアルタイムに、精度高く計測、表示することができると共に、身の回りに装着し易い小型化により、放射線作業従事者らの放射線被曝への意識を高めさせることができるγ線計測技術を提供することを課題とする。

本発明者は、上記した課題の解決について鋭意検討を行い、以下に記載する発明により上記課題が解決できることを見出し、本発明を完成させるに至った。

請求項１に記載の発明は、
γ線の線量とエネルギースペクトルとを、同時に計測するγ線計測方法であって、
検知部に検知されたγ線の強度データに基づいて、ベイズ推定法を用いて、エネルギースペクトルを作成するエネルギースペクトル作成ステップと、
作成された前記エネルギースペクトルに、線量変換係数を乗じることにより、線量を算出する線量算出ステップとを備えており、
前記エネルギースペクトル作成ステップにおけるベイズ推定法が、逐次型ベイズ推定法であり、
修正の度合いを示すα値を、カウント数および減少指数に基づいて、下式で規定することを特徴とするγ線計測方法である。

請求項２に記載の発明は、
前記減少指数が、１０^－５～１０^－３であることを特徴とする請求項１に記載のγ線計測方法である。

請求項３に記載の発明は、
γ線の線量とエネルギースペクトルとを、同時に計測するγ線計測装置であって、
到来したγ線を検知する検知部と、
検知されたγ線の強度を電気信号に変換する変換部と、
変換された電気信号に基づいて、到来したγ線のエネルギースペクトルを作成するエネルギースペクトル作成部と、
作成されたエネルギースペクトルに基づいて、到来したγ線の線量を算出する線量算出部と、
作成されたエネルギースペクトルと算出された線量とを、計測データとして、同時に表示する表示部とを備えており、
前記エネルギースペクトル作成部に、請求項１または請求項２に記載のγ線計測方法が、予め、プログラムとして組み込まれていることを特徴とするγ線計測装置である。

請求項４に記載の発明は、
前記検知部が、到来したγ線のエネルギーに比例して発光するシンチレータであり、
前記変換部が、発光した光の光子をカウントして、γ線のエネルギーに比例した波高値に対応するパルス信号に変換するＭＰＰＣ（Ｍｕｌｔｉ－Ｐｉｘｅｌ Ｐｈｏｔｏｎ Ｃｏｕｎｔｅｒ）であることを特徴とする請求項３に記載のγ線計測装置である。

請求項５に記載の発明は、
前記シンチレータが、ヨウ化セシウム（ＣｓＩ）であることを特徴とする請求項４に記載のγ線計測装置である。

請求項６に記載の発明は、
前記シンチレータが、２～３ｃｍ角の立方体に加工されたヨウ化セシウムであることを特徴とする請求項５に記載のγ線計測装置である。

本発明によれば、線量だけでなくγ線のエネルギースペクトルをも、その場でリアルタイムに、精度高く計測、表示することができると共に、身の回りに装着し易い小型化により、放射線作業従事者らの放射線被曝への意識を高めさせることができるγ線計測技術を提供することができる。

本発明の一実施の形態に係るγ線計測装置における電子回路を説明するブロック図である。 γ線のエネルギースペクトルの不安定性を説明する図である。 減少指数ｋを説明する図である。 γ線の波高スペクトルからγ線のエネルギースペクトルへの展開を説明する図である。 Ｇｅ検出器によるγ線の波高スペクトルの計測結果を示す図である。 繰り返し数の検討を説明する図である。 本発明の一実施の形態において計測された線量の時間的変化を示す図である。 既知の複数のエネルギーのγ線場において得られた波高スペクトルおよび推定エネルギースペクトルを示す図である。 γ線フラックスの算出を説明する図である。 ^１３７Ｃｓのピークに基づいて算出されたγ線フラックスを説明する図である。 ^６０Ｃｏの低い方のピークに基づいて算出されたγ線フラックスを説明する図である。 ^６０Ｃｏの高い方のピークに基づいて算出されたγ線フラックスを説明する図である。 ^１３３Ｂａのピークに基づいて算出されたγ線フラックスを説明する図である。 図１０～図１３に基づいて得られた線量の推定値と、測定時間との関係を示す図である。 バックグランドでの６０秒計測によるエネルギースペクトルを示す図である。 バックグランドでの１時間計測によるエネルギースペクトルを示す図である。 バックグランドでの計測に適切な繰り返し数Ｉを調べた結果を示す図である。 バックグランドでの計測で得られた線量の例を示す図である。 本発明のベイズ推定法の適切な使用領域を示す図である。

［１］γ線のエネルギースペクトルのリアルタイム表示
前記したように、γ線の線量を精度高く計測するためには、γ線のエネルギースペクトルをリアルタイムに知る必要がある。

このため、最初に、このエネルギースペクトルのリアルタイム表示について説明する。

１．基本的な考え方
離散化した計測装置の出力（波高分布）ｙ_ｉは、放射線の情報（エネルギー分布）ψ_ｊ、および、あるエネルギーＥ_ｊにおける検出器応答Ｒ_ｉ，ｊと、以下の関係があることが知られており、この式（１）を解くことにより、エネルギー分布ψ_ｊを得ることができる。

しかし、実際には、波高分布ｙ_ｉに計測誤差、検出器応答Ｒ_ｉ，ｊに評価誤差があるため、工学的には式（１）が成立せず、一般的に式（２）のようになり、解くことが難しい。

なお、この逆問題の解法については、例えば、最尤法やモンテカルロ法などを用いた汎用コードを用いることにより近似的に解けることが知られているが、これらの方法は、原理的にポストプロセス（計測後処理）によるデータ処理であるため、その場でリアルタイムに解く本発明に適用することはできない。

２．ベイズ推定法
そこで、本発明者は、式（２）の逆問題の解法として、ベイズ推定法を用いることを考えた。

（１）ベイズ推定法
ベイズ推定法は、いわゆるベイズの定理に基づく推定法であり、ベイズ推定法を用いることにより、γ線計測装置で計測される波高分布からエネルギー分布が推定できる。なお、このことについては、東北大学の岩崎等により既に示されている（Ｓ．Ｉｗａｓａｋｉ，“Ａ ｎｅｗ ａｐｐｒｏａｃｈ ｆｏｒ ｕｎｆｏｌｄｉｎｇ ｐｒｏｂｌｅｍｓ ｂａｓｅｄ ｏｎｌｙ ｏｎ ｔｈｅ Ｂａｙｅｓ’Ｔｈｅｏｒｅｍ”，Ｐｒｏｃ．ｏｆ ｔｈｅ ９^ｔｈ Ｉｎｔ．Ｓｙｍｐ．ｏｎ Ｒｅａｃｔｏｒ Ｄｏｓｉｍｅｔｒｙ，ｐｐ．２４５－２５２（１９９６）参照）。

（２）スペクトル型ベイズ推定法
スペクトル型ベイズ推定法は、検出器による計測終了後にベイズ推定法によりエネルギー分布を推定する方法である。この方法は通常の放射線計測で使用されているものである。

しかしながら、この手法は、通常、測定が終了して測定データが取得された後（ポストプロセス）にエネルギー分布を求める場合に適用できる手法であり、その場でリアルタイムに解く本発明には適用できない。

（３）逐次型ベイズ推定法
そこで、本発明者は、リアルタイム計測を実現するため、逐次型ベイズ推定法を用いることを考えた。

逐次型ベイズ推定法は、検出器からの出力（波高値〉が１個（１カウント）得られるたびに、エネルギースペクトルの推定をベイズの定理により実施するもの（リアルタイム処理）であり、修正の度合いを示すα値を設定することにより、ある１カウント検出後の推定スペクトルを、以下の式で示すことができる。

上式は、例えば、α＝０．１の場合、新たな推定結果は、前回の推定結果に対して、１０％だけ上乗せされるということになることを示している。

しかしながら、この逐次型ベイズ推定法の採用については、従来、精度が悪く（測定が進んでカウント数が増えても結果が収束せず、統計精度が向上しない）、また、結果（カウント数に対する推定スペクトル）に不安定性があるという評価がなされていた。

本発明者は、この逐次型ベイズ推定法における精度の向上について、種々の実験と検討を行い、その結果、上記したα値を、計測場の計数率などの環境を勘案して設定、具体的には、計測を実施するたびにα値を良い結果が出る正確な数値に設定した場合には、実用上、十分に受け入れ可能な精度にまで推定誤差を抑え込んで、推定できることを見出した。

しかしながら、結果の不安定性については依然として認められ、未だ十分とは言えなかったため、さらに、α値の変化による収束の傾向などを、計測実験毎に詳しく調べ、その解決について検討を行った。

その結果、逐次型ベイズ推定法にα値の変化のパラメータとしての定数ｋ（減少指数）を導入し、α値をこのｋの関数として、カウント数に対して動的に変化させることにより、収束が著しく改善され、受け入れ可能な精度にまで十分に推定誤差を抑え込むことに加えて、不安定性が回避できることを見出した。

そして、このように、α値をｋの関数として、カウント数に対して動的に変化させて、逐次型ベイズ推定法に基づいて高い精度で不安定性が回避されたγ線エネルギースペクトルをリアルタイムに作成することにより、エネルギースペクトルに線量変換係数を乗じて算出される線量も、全てのエネルギーを厳密に考慮した状態で、リアルタイムに精度高く算出することができる。

そして、このような処理プロセスが組み込まれたγ線計測装置は、線量とエネルギースペクトルの両方を真の値に近い状態で、その場でリアルタイムに表示して評価することができる。また、γ線の検知素子としてＣｓＩのような小型化が可能な素子を用い、検知されたγ線の強度を電気信号に変換する変換部としてＭＰＰＣのようなＣｓＩと相性がよく、小型の変換器を用いることにより、身の回りに装着し易い程度にまで容易に小型化することができる。

この結果、このようなγ線計測装置は、放射線作業従事者らの放射線被曝への意識を高めさせることができ、これまでのサーベイメータや個人線量計を補完する、もしくは代替させることができる。

［２］実施の形態
以下、具体的な実施例を交えた実施の形態に基づき、図面を参照しつつ説明する。

１．本実施の形態に係るγ線計測装置
最初に、本実施の形態に係るγ線計測装置について説明する。

図１は、本実施の形態に係るγ線計測装置における電子回路を説明するブロック図である。図１に示すように、本実施の形態に係るγ線計測装置１は、以下の構成要素から形成されている。

（１）到来したγ線を検知する検知部として、γ線の強度（エネルギー）に比例して発光するシンチレータ１１
（２）検知されたγ線の強度を電気信号に変換する変換部として、シンチレータ１１に検知されることによって発光した光の光子をカウントして、γ線のエネルギーに比例した波高値に対応するパルス信号に変換するＭＰＰＣ１２
（３）変換された電気信号からＡＤＣ（Ａｎａｌｏｇ ｔｏ Ｄｉｇｉｔａｌ Ｃｏｎｖｅｒｔｅｒ）１６によりディジタル変換されたパルス信号に対して、後述する逐次型ベイズ推定法を適用することにより、到来したγ線のエネルギースペクトルを作成するエネルギースペクトル作成部１８
（４）作成されたエネルギースペクトルに基づいて、到来したγ線の線量を算出する線量算出部（エネルギースペクトル作成部１８の「ｄｏｓｅ」に対応）
（５）作成されたエネルギースペクトルと算出された線量とを、計測データとして、同時に表示する表示部（図示せず）

なお、１３はＭＰＰＣへ電力を供給する電源部であり、１４はＭＰＰＣ１２により変換されたパルス信号を増幅する増幅器である。

そして、本実施の形態に係るγ線計測装置１においては、後述するように、作成されたエネルギースペクトルが、高い精度を有しており、また、安定性にも優れているため、線量も精度高い値で算出される。

なお、本実施の形態において、シンチレータ１１としては、潮解性が無く、容易に加工して小型化が可能なＣｓＩが好ましく、通常は、タリウム（Ｔｌ）を少量添加することにより活性化させたＣｓＩ（Ｔｌ）（タリウム活性化ヨウ化セシウム）が使用される。なお、シンチレータ１１として使用する際のサイズとしては、２～３ｃｍ角の立方体に加工したものが好ましく、具体的な一例として、ＣｓＩ（Ｔｌ）の２．６ｃｍ角の立方体（重量：１００ｇ弱）が使用される。

このＣｓＩは、光電子増倍管に比べて格段に軽いＭＰＰＣとの相性がよく、ＣｓＩとＭＰＰＣを組み合わせることにより、優れたエネルギー分解能と検出効率の確保を図ると共に、γ線計測装置の小型化、軽量化を図ることができる。

２．γ線のエネルギースペクトルの作成
次に、本実施の形態におけるγ線のエネルギースペクトルの作成について説明する。

ベイズの定理によると、チャンネルｉにおける波高分布ｙ_ｉから、ｊにおけるエネルギースペクトルの事後確率を、以下の式のように、検出器の応答関数Ｒで導出することができる。

そして、逐次型ベイズ推定法においては、検出されたチャンネルｉでの信号（波高信号）１カウントごとに寄与を加算していくことにより、下式のように、エネルギースペクトルψ_ｊを改定して、推定していく。なお、αは、事前確率として、どの程度、事後確率が制御されるのかを調整する係数（修正の度合いを示す係数）である。

しかし、上記の式を適用した場合、γ線のエネルギースペクトルのピークの位置がカウント数によってずれたり、強度が変動する場合があり、得られるγ線のエネルギースペクトルに不安定性がある。

具体的な一例を図２に示す。図２は、天然ウランの核燃料庫の内部における計測例であり、線量は、～２μＳｖ／ｈｒ程度である。なお、図２においては、α＝０．０１として、５時間（下図）および３０秒間（上図）、γ線の計測を行って、図２左に示すγ線の波高スペクトルを得た後、逐次型ベイズ推定法を用いて、図２右に示すエネルギースペクトルを得ている。

図２より、γ線の波高スペクトルでは両者の波高分布の形状変化が小さいが、統計精度は５時間計測の方が格段に優れていることが分かる。しかし、エネルギースペクトルでは、両者の間に、ピーク位置のずれが生じており、しかも、５時間計測の方がむしろ精度が悪くなっており、エネルギースペクトルが安定していないことが分かる。これに対して、３０秒計測の方では、１ＭｅＶ付近にはっきりとしたピークが認められ、天然ウランからのピークであることが確認されている。

本発明者が、さらに、α値を変化させて、その不安定性を見た場合、α値が大き過ぎるとこの不安定性が大きくなり、一方、α値が小さ過ぎると改定スピードが遅くなることが分かり、結果として、計測の開始時にはα値を大きくし、その後は、徐々に小さくしていく必要性があることが分かった。

本発明者は、このようなα値を適切に変化させる手法について検討し、図３に示すようにα値の減少指数ｋを規定して、この減少指数ｋおよびカウント数Ｎに基づいて、下式のようにα値を規定すればよいことを見出した。

本発明者は、次に、上記したｋの好ましい値について検討を行った。具体的には、図４に示すように、１秒から５時間までのγ線の計測時間におけるγ線の波高スペクトルに対し、ｋを１０^－２０から１０^－１まで変化させて、逐次型ベイズ法を用いて、リアルタイムのエネルギースペクトルに展開（ｕｎｆｏｌｄｉｎｇ）した。その結果、ｋは１０^－５～１０^－３であることが好ましく、特に、１０^－４の場合、シンチレータによるγ線の強度の計測時間が６０秒あれば、図４左に示す波高スペクトルの３つのピーク、即ち、０．１ＭｅＶ付近における^２３４Ｔｈのピーク、０．８ＭｅＶ付近、および１ＭｅＶ付近における^２３４Ｐａのピークが、図４右に示すエネルギースペクトルにおいても再現されており、特に好ましいことが分かった。なお、図４は、天然ウランの核燃料庫の内部における計測例（線量～２μＳｖ／ｈｒ）であり、波高スペクトルからエネルギースペクトルへの展開を説明する図である。

そして、上記で得られたγ線のエネルギースペクトルを、Ｇｅ検出器によるγ線の波高スペクトルの計測結果（図５）と比較すると、概ね似た形となっており、また、同じ位置にピークが現れていることが確認でき、正しく計測できていることが確認できた。

そして、本実施の形態に係るγ線計測装置の場合、上記したように、線量が２μＳｖ／ｈｒ程度の強度の場では、少なくとも６０秒の計測を行うことにより、リアルタイムで精度の高いエネルギースペクトルが得られることが確認できた。

一方、結果の不安定性は、逐次型ベイズ推定における繰り返し数Ｉに関係しているため、次に、好ましい繰り返し数Ｉについて検討した。なお、この繰り返し数Ｉは、逐次型ベイズ推定中に、１／ｋカウントの計測が何回繰り返されるかを示す指標であり（図６（ａ）参照）、下式によって定義されている。

その結果、図６（ｂ）～（ｅ）に示すように、繰り返し数Ｉは約５０とすることが特に好ましいことが分かった。

即ち、Ｉが４５である図６（ｂ）や、Ｉが４７である図６（ｃ）では、ピークが明確に現れており、位置もずれていない。一方、Ｉが１００を超える図６（ｄ）では、ピークが歪んだり、その位置がずれるなど、αが大きいことによる不安定性が現れている。そして、Ｉが１未満の図６（ｅ）では、十分な改定が行われないため、カウント数が増えても、スペクトルの形に変化が現れず、ピークを把握することが難しい。なお、図６（ｄ）における具体的なＩは、（１）４５０、（２）２７００、（３）５４００、（４）２７０００である。また、図６（ｅ）における具体的なＩは、（１）１．５Ｅ－７、（２）１．５Ｅ－６、（３）４．５Ｅ－６、（４）９Ｅ－６、（５）１．３５Ｅ－５、（６）４．５Ｅ－５、（７）２．７Ｅ－４、（８）５．４Ｅ－４、（９）２．７Ｅ－３である。

３．γ線のエネルギースペクトルから線量への変換
次に、得られたγ線のエネルギースペクトルから線量への変換について説明する。

前記したように、エネルギースペクトルに線量変換係数を乗じることにより、線量を求めることができる。

図７は、図４と同様にして、ｋを１０^－２（０．０１）から１０^－５（０．００００１）まで変化させて得られたγ線のエネルギースペクトルから求められた線量の時間的変化を示す図である。なお、図７には、比較のために、従来のサーベイメータで計測された線量も記載している。

図７より、安定した線量値は、ｋ≦１０^－４の場合に得られており、従来のサーベイメータで計測された線量よりも低い値となっていることが分かる。これは、従来のサーベイメータで計測された線量が、^１３７Ｃｓで一点校正されているためである。

次に、本実施の形態のγ線計測装置において得られた線量値が正しい値であるか否かを確認するために、３つの標準線源、即ち、^１３７Ｃｓ（６６２ＫｅＶ）、^６０Ｃｏ（１．３ＭｅＶ）、^１３３Ｂａ（３５６ＫｅＶ）について、本実施の形態のγ線計測装置、および、従来のサーベイメータを用いて線量の計測を行い、得られた各値を理論値と比較した（ｋ＝１０^－４）。結果を表１に示す。

表１から分かるように、本実施の形態のγ線計測装置では、理論値とほぼ一致した値が得られている。これに対して、従来のサーベイメータでは、^１３７Ｃｓ以外は理論値と合致していない。これは、従来のサーベイメータが^１３７Ｃｓのγ線一点で校正されているためである。

そして、核燃料庫における平均エネルギーは、^１３７Ｃｓよりも低いことを考慮すると、たとえ、本実施の形態のγ線計測装置で、^１３７Ｃｓで一点校正された従来のサーベイメータよりも低い値が表示されたとしても、図７に示すように、正しい値を示していることが分かる。

そして、本実施の形態に係るγ線計測装置は、スマートフォン程度のコンパクトさで、重量も２００ｇ程度とすることができるため、放射線作業従事者らが、身の回りに装着し易い。

４．標準線源場での測定
上記においては、天然ウランの核燃料庫の内部における計測例を示したが、以下においては、線源強度が理論的に分かっている標準線源場での計測例を挙げて、本発明においては、短時間にγ線のフラックスのピークのエネルギーを推定して、３０秒以内に、線量を正しく評価できることを説明する。

具体的には、まず、^１３７Ｃｓ、^６０Ｃｏ、^１３３Ｂａ、３種類の標準線源を用いて、既知の複数のエネルギーのγ線場を製作し、６００秒間の計測を行い、図８上に示す波高スペクトルを得た後、得られた波高スペクトルに対して、ｋの値を１０^－３、１０^－４、１０^－５と変化させて、逐次型ベイズ法を用いて展開し、図８下に示す２０秒後の推定エネルギースペクトルを得た。

得られた推定エネルギースペクトルでは、図８下に示すように、^１３７Ｃｓに基づく１つのピーク、^１３３Ｂａに基づく１つのピーク、^６０Ｃｏに基づく２つのピークがあることが分かる。

次に、各ピークのγ線フラックスについて、逐次型ベイズ法に基づく推定値Ｇｆ_ｅｓｔの、図９に示す式から算出される理論値Ｇｆ_ｒｅｆに対する比（Ｇｆ_ｅｓｔ／Ｇｆ_ｒｅｆ）と、測定時間との関係について、ｋの値を１０^－３、１０^－４、１０^－５と変化させて、算出した。結果を図１０～図１３に示す。（Ｇｆ_ｅｓｔ／Ｇｆ_ｒｅｆ）が１に近いほど、推定値と理論値が一致していることを示す。

図１０は、^１３７Ｃｓに基づくピークについての結果であり、ｋ＝１０^－４の場合には、６０秒で（Ｇｆ_ｅｓｔ／Ｇｆ_ｒｅｆ）が、１．１に安定的に収束しているが、ｋ＝１０^－３の場合には、不安定性が生じており、ｋ＝１０^－５の場合には、不十分な収束となっている。この結果より、^１３７Ｃｓの場合、ｋを１０^－４と設定することにより、６０秒ほどで（Ｇｆ_ｅｓｔ／Ｇｆ_ｒｅｆ）が１近くに収束して、６６２ｋｅＶのピークにおけるγ線フラックスの推定が可能であることが分かる。

図１１は、^６０Ｃｏに基づく低い方のピークについての結果であり、ｋ＝１０^－４の場合には、６０秒で（Ｇｆ_ｅｓｔ／Ｇｆ_ｒｅｆ）が、０．９に安定的に収束しているが、ｋ＝１０^－３の場合には、不安定性が生じており、ｋ＝１０^－５の場合には、不十分な収束となっている。この結果より、^６０Ｃｏの場合、ｋを１０^－４と設定することにより、６０秒ほどで（Ｇｆ_ｅｓｔ／Ｇｆ_ｒｅｆ）が１近くに収束して、１．１７ＭｅＶのピークにおけるγ線フラックスの推定が可能であることが分かる。

図１２は、^６０Ｃｏに基づく高い方のピークについての結果であり、ｋ＝１０^－４の場合には、６０秒で（Ｇｆ_ｅｓｔ／Ｇｆ_ｒｅｆ）が、０．９に収束し、ｋ＝１０^－５の場合には、６０秒で（Ｇｆ_ｅｓｔ／Ｇｆ_ｒｅｆ）が、０．９５に収束しているが、ｋ＝１０^－３の場合には、不安定性が生じている。この結果より、^６０Ｃｏの場合、ｋを１０^－４あるいは１０^－５と設定することにより、６０秒ほどで（Ｇｆ_ｅｓｔ／Ｇｆ_ｒｅｆ）が１近くに収束して、１．３３ＭｅＶのピークにおけるγ線フラックスの推定が可能であることが分かる。

図１３は、^１３３Ｂａに基づくピークについての結果であり、ｋ＝１０^－３、１０^－４、１０^－５のいずれの場合においても、３０秒で（Ｇｆ_ｅｓｔ／Ｇｆ_ｒｅｆ）が、１．４近くに収束している（ｋ＝１０^－３の場合、若干、不安定な収束）。ピークが０．３５６ＭｅＶ、０．３８４ＭｅＶと６００ｋｅＶ以下であり、後方散乱の影響があるため、（Ｇｆ_ｅｓｔ／Ｇｆ_ｒｅｆ）は１からはズレているが、ｋを１０^－４あるいは１０^－５と設定することにより、３０秒ほどで（Ｇｆ_ｅｓｔ／Ｇｆ_ｒｅｆ）が素早く収束することが確認できた。

なお、ｋは前記した通り、１０^-５～１０^－３が好ましく、図１０～図１３の説明において、１０^-５、１０^－３の場合に、不十分ないし不安定という説明をしたが、これは最も好ましい１０^－４との比較において説明したものであり、ｋが１０^-５～１０^－３であれば、実用的に適用可能な収束はできている。

次に、図１０～図１３の結果に基づいて、線量の評価を行う。線量は、（測定時の放射能）×（空気衝突カーマ率定数）／（シンチレータと線源の距離）^２の式より求めることができ、理論値は、^６０Ｃｏ、^１３３Ｂａ、^１３７Ｃｓ、それぞれの線量の合計として示すことができる。一方、推定値としては、推定された各エネルギースペクトルに線量変換係数をかけて合計した値で示すことができる。

図１４は、図１０～図１３に基づいて得られた線量の推定値と、測定時間との関係を示す図であり、ｋ＝１０^－3やｋ＝１０^－5でも推定可能であるが、ｋ＝１０^－3の場合には、若干不安定性が見られ、ｋ＝１０^－5の場合には、少し収束が遅いことが分かる。ｋ＝１０^－４の場合、最も安定しており、３０秒以内で正しい線量を評価できることが分かる。具体的な結果として、図１４において、ｋ＝１０^－４、３０秒の場合、推定値は６．００６μＳｖ／ｈｒとなり、理論値６．０１３μＳｖ／ｈｒとほぼ一致している。

以上より、本実施の形態によれば、複数の標準ガンマ線源を用いた線量の計測が可能であることが確認できた。

５．バックグランドγ線場での測定
上記においては、天然ウランの核燃料庫の内部における計測例を示したが、以下においては、本発明を実際に使用する場合、最も線量が低い場と考えられるバックグランドγ線場での計測例を示し、バックグランドγ線場においても、即ち、ほぼ全ての場所で、本発明が適用できることを説明する。

図１５は、放射線源が存在しない室内、即ち、バックグランドγ線場でγ線を６０秒、計測した例であり、図１５上に示すγ線の波高スペクトルを得た後、本発明による手法を用いて、図１５下に示すエネルギースペクトルを推定している。

図１５下のエネルギースペクトルにおいては、^４０Ｋ（１４６１ＫｅＶ）と、^２０８Ｔｌ（２６１４．５ＫｅＶ）のピークが現れている一方、ピークがないことが予め分かっている？の位置にピークが見られ、安定した測定ができていないことが分かる。

図１６は、図１５の場合と同様の条件で、γ線を１時間計測した結果である。図１６の場合は、図１５下に現れていた？の位置のピークが、消えていることが分かる。このことは、バックグランドγ線場でも、測定は可能であるが、放射線が弱いため、安定した測定をするには、時間がかかることを示している。

なお、図１５、図１６の場合、測定時間、１週間で、総カウント数は２８００万カウント（４６ＣＰＳ）である。

次に、図１７は、バックグランドγ線場での計測に適切な繰り返し数Ｉを調べた結果であり、繰り返し数Ｉを、Ｉ＝１、Ｉ＝２０、Ｉ＝５０にした時のエネルギースペクトルを示している。（ａ）のＩ＝１の場合、ピークがないことが予め分かっている丸囲みした部分にピークが表れており、安定した測定ができていないことが分かる。

一方、（ｂ）のＩ＝２０の場合、Ｉ＝１の場合に丸囲みした部分に表われていたピークが消えていることが分かる。また、得られたピークは、Ｉ＝１の場合より明瞭に表れており、精度高く測定できていることが分かる。

（ｃ）のＩ＝５０の場合においても、Ｉ＝１の場合に丸囲みした部分に表われていたピークが消えていることが分かり、また、各ピークが、Ｉ＝２０の場合より明瞭に表れており、精度高く測定できていることが分かる。一方、Ｉ＝５０を超えると、安定した測定ができないことが実験で確認できている。このため、Ｉ＝２０～５０の間で測定することが好ましい。

そして、繰り返し数Ｉを変化させて得られた図１７の３つのエネルギースペクトルのいずれにおいても、ｋの値に拘わらず、ほぼ同じエネルギースペクトルを示しており、推定の精度は、ｋの値に関係なく、繰り返し数Ｉによって制御されていることが分かる。

次に、バックグランドでの計測に必要な測定時間を検討した。

前記した通り、Ｉ＝ｋＮであるため、繰り返し数Ｉが決まっている場合、ｋが小さいほど、カウント数Ｎが多くなり、必要な測定時間が長くかかる。

前記の通り、推定の精度は、ｋに関係なく、繰り返し数Ｉによって制御されるため、測定時間の短縮のためには、ｋは支障のない範囲で大きいことが好ましく、図１７の場合、ｋが０．００１の場合でも、精度高い測定ができているため、Ｉ＝２０で、ｋ＝０．００１とすることが好ましい。この場合、カウント数Ｎは２００００となり、安定した測定に必要な時間は７．５分となる。

なお、Ｉ＝２０において、ｋ＝０．０００１とすると、カウント数Ｎは２０００００となり、安定した測定に必要な時間は１．２時間となる。また、Ｉ＝２０において、ｋ＝０．００００１とすると、カウント数Ｎは２００００００となり、安定した測定に必要な時間は１２時間となる。

以上より、天然ウランの核燃料庫の内部における計測の場合、上記した通り、６０秒で、精度の高いエネルギースペクトルが得られていることと比較すると、バックグランドでの計測の場合、長い計測時間を要することが、実験の結果から分かる。

図１８は線量の結果を示すグラフであり、線量変換係数で推定したエネルギースペクトルに線量変換係数を乗じることで求めている。図１８より、ｋが０．００１より小さいと、線量は、２０秒程度で不安定性もなく収束することが分かる。

６．まとめ
上記において、本発明によるベイズ推定法に基づいた天然ウランの核燃料庫の内部における計測例および、バックグランドγ線場での計測例を示したが、上記以外の実験結果も含めて、種々の実験により得られた結果を表２に示す。そして、表２より以下のことが分かる。

空間線量が、核燃料庫の内部や、加速器運転後、適当な時間経過後の線量である２μＳｖ／ｈｒ程度であれば、スペクトルも、線量も十分に推定が可能である。空間線量が０．０８μＳｖ／ｈｒ程度のバックグランドの場合でも、スペクトル、線量共に、十分に推定が可能であるが、スペクトルについては、７．５分程度の時間を要す。

測定時間に関しては、実験において、１インチ角のＣｓｌのシンチレータを用いたが、シンチレータのサイズを大きくすることにより、精度が向上し、測定時間を短くすることができる。

繰り返し数Ｉについては、数十程度が適切である。

また、αの減少指数ｋについては、バックグランドの場合や、計数率が小さい場合には、１０^-３程度、加速器運転後などの場合には１０^－４程度、ＢＮＣＴなどのＯｎ-ｂｅａｍの場合には１０^－５程度が適当である。

次に、図１９に、本発明のベイズ推定法の適切な使用領域を示す。Ｉ＝１０の場合は、スペクトルの評価の安定度が悪く、Ｉ＝４０付近の方が好ましい。しかし、バックグランドでの計測まで考慮すると、Ｉ＝４０では、スペクトルの評価に少し時間がかかるので、計数値を上げる別の方策も考慮することが好ましい。

ただし、バックグランドでの計測の場合、スペクトルの評価まで必要な場合が、どの程度あるかを考えておく必要がある。また、医療現場での計測では、バックグランドよりずっと高い線量である。測定時間については、シンチレータの結晶を大きくすればよいと言えるが、実用上は、あまり大きくしたくないものである。このように種々の考慮すべき事項があり、これらを総合的に考慮して、適切な使用領域を決めることになる。

次に、実用面からの考察結果について説明する。実用上、効率的な計測には、ｋの値を短時間に適切に決めることが重要となる。本発明者が行ったこれまでの結果を考慮するとｋは、計数率ＣＰＳ、計測時間、そしてＩが決まれば決めることができ、具体的には、ｋは以下の式で求められる。

ここで、Ｉは上記した通り、数十程度が適当である。計測時間ｔは、実用上、数十秒以下にしたいものである。このため、Ｉ、ｔは適宜決めることができる。そして、１秒当たりのカウント数ＣＰＳは、計測場所で決まる。このため、適切なｋを自動的に決めることができる。ｋを自動的に決めることができることは、実用上、大きな意義を有すると言える。また、ｋを１０^－３、１０^－４、１０^－５などからスイッチで選択できるようにしておく方法も、実用上、好ましい方法と言える。

以上、本発明を実施の形態に基づいて説明したが、本発明は、上記の実施の形態に限定されるものではない。なお、本発明と同一および均等の範囲内において、上記の実施の形態に対して種々の変更を加えることが可能である。

１ γ線計測装置
１１ シンチレータ
１２ ＭＰＰＣ
１３ 電源部
１４ 増幅器
１６ ＡＤＣ
１８ エネルギースペクトル作成部

Claims (6)

1. γ線の線量とエネルギースペクトルとを、同時に計測するγ線計測方法であって、
   検知部に検知されたγ線の強度データに基づいて、ベイズ推定法を用いて、エネルギースペクトルを作成するエネルギースペクトル作成ステップと、
   作成された前記エネルギースペクトルに、線量変換係数を乗じることにより、線量を算出する線量算出ステップとを備えており、
   前記エネルギースペクトル作成ステップにおけるベイズ推定法が、逐次型ベイズ推定法であり、
   修正の度合いを示すα値を、カウント数および減少指数に基づいて、下式で規定することを特徴とするγ線計測方法。
   

   
2. 前記減少指数が、１０^－５～１０^－３であることを特徴とする請求項１に記載のγ線計測方法。
3. γ線の線量とエネルギースペクトルとを、同時に計測するγ線計測装置であって、
   到来したγ線を検知する検知部と、
   検知されたγ線の強度を電気信号に変換する変換部と、
   変換された電気信号に基づいて、到来したγ線のエネルギースペクトルを作成するエネルギースペクトル作成部と、
   作成されたエネルギースペクトルに基づいて、到来したγ線の線量を算出する線量算出部と、
   作成されたエネルギースペクトルと算出された線量とを、計測データとして、同時に表示する表示部とを備えており、
   前記エネルギースペクトル作成部に、請求項１または請求項２に記載のγ線計測方法が、予め、プログラムとして組み込まれていることを特徴とするγ線計測装置。
4. 前記検知部が、到来したγ線のエネルギーに比例して発光するシンチレータであり、
   前記変換部が、発光した光の光子をカウントして、γ線のエネルギーに比例した波高値に対応するパルス信号に変換するＭＰＰＣ（Ｍｕｌｔｉ－Ｐｉｘｅｌ Ｐｈｏｔｏｎ Ｃｏｕｎｔｅｒ）であることを特徴とする請求項３に記載のγ線計測装置。
5. 前記シンチレータが、ヨウ化セシウム（ＣｓＩ）であることを特徴とする請求項４に記載のγ線計測装置。
6. 前記シンチレータが、２～３ｃｍ角の立方体に加工されたヨウ化セシウムであることを特徴とする請求項５に記載のγ線計測装置。

Images