JP6722952B2 - シンチレーション光検出装置および放射線検出装置 - Google Patents

Description

本装置は、複数の単位シンチレータが配列されて形成されたシンチレータ配列において、どの単位シンチレータに放射線が入射したかを特定するシンチレーション光検出装置およびこれを用いた放射線検出装置に関する。

癌の診断手法の一つとして、ＰＥＴ（Positron Emission Tomography）検査がある。このＰＥＴ検査では、癌細胞に集まりやすい糖に放射性物質を結合させた検査薬を患者に投与する。そして、癌細胞に集まった検査薬中の放射性物質が陽電子崩壊を起こし、生じた陽電子が消滅する際に発生する一対のγ線（以下、消滅γ線という）を検出する。特に、この一対の消滅γ線は、１８０度異なる反対方向に進む。従って、検体の周囲に配置した放射線検出器で同時に入射した一対の消滅γ線の位置を検出することで、検出位置を結んだ直線上に放射性物質が存在することから、これを繰り返し測定してＣＴ（Computed Tomography）のように再構成画像を作成することにより、癌細胞などの位置を特定することができる。

このようなＰＥＴ検査装置では、患者を取り巻く空間に多数の検出器を配置し、どの検出器ペアに消滅γ線が入射したかを特定する。検出器の数が多いほど放射線検出の効率が高くなるため、多数の検出器を配置することが望ましいが、検出器の数が増えるほど同時検出イベントの処理やコストも増大するため、シンチレータと位置敏感型光電子増倍管を用いた安価なシステムが採用されることが多い。

ここで、従来のＰＥＴ検出器では、多数の検出器を患者の周囲にリング状に配置しているが、患者に向かう方向（深さ方向）に対する分解能を有しておらず、位置の特定が困難である。従って、消滅γ線が検出器に対して斜め方向から入射する場合に位置分解能が低くなる。

これを解決する手段として、ＤＯＩ（Depth of Interaction）検出器が提案されている。このＤＯＩ検出器では、深さ方向に複数のシンチレータを配置し、これらのシンチレータからの発光信号を識別するための工夫が施されており、深さ方向の弁別が可能となっている。したがって、斜め方向から入射する消滅γ線についての位置分解能を高くできる。

ＷＯ２０１１／１２１７０７号公報

ＤＯＩ検出器では、深さ方向に配置された複数のシンチレータからの信号を１つの光電子増倍管に入力し、深さ方向のいずれにシンチレータからの信号かを弁別する。この信号の弁別には、特別の工夫が必要であり、必ずしも容易ではなかった。また、深さ方向におけるシンチレータの個数をあまり大きくできないという問題もあった。

本シンチレーション光検出装置は、放射線が入射するイベントによって生じる発光の減衰時定数が互いに異なる複数のシンチレータ素子が配列されたシンチレータ配列と、前記シンチレータ配列の複数のシンチレータ素子から出力される光を受光し、電気信号に変換する光電変換器と、前記イベントについて前記光電変換器からの電気信号の電圧波形におけるピーク値と積分電荷量を検出し、検出したピーク値と積分電荷量の比に応じて、当該電気信号が前記シンチレータ配列中のどのシンチレータ素子に放射線が入射したことに起因する電気信号かを特定する特演算処理装置と、を含む。なお、積分電荷量は、電気信号の積分値で代用することが可能である。

前記シンチレータ配列は、複数のシンチレータ素子が１列に配列され、当該１列の一端のシンチレータ素子から複数のシンチレータ素子からの光が出力されることができる。

前記シンチレータ配列を複数配列し、前記光電変換器は、個別のシンチレータ配列からの光を別々に受光して、別々の電気信号を出力し、前記演算処理装置は、各電気信号の電圧波形におけるピーク値と積分電荷量の比に応じて、各電気信号がどのシンチレータ配列のどのシンチレータ素子に放射線が入射したことに起因する電気信号かを特定することができる。

前記減衰時定数が互いに異なる複数のシンチレータ素子は、賦活材の濃度が互い異なるものを含むことができる。

前記減衰時定数が互いに異なる複数のシンチレータ素子は、シンチレータ素子を構成する材料が互いに異なるものを含むことができる。

前記減衰時定数が互いに異なる複数のシンチレータ素子は、賦活材の濃度が互いに異なるとともに、シンチレータ素子を構成する材料が異なるものを含むことができる。

消滅γ線のエネルギーは必ずしも１つのシンチレータ素子で吸収されるとは限らず、複数のシンチレータ素子にまたがってエネルギーが吸収され、誤った位置として検出される場合（以下、クロストークイベントという）がある。前記材料が異なる複数のシンチレータ素子についての出力エネルギーおよび各電気信号の電圧波形におけるピーク値と積分電荷量の比の相違を利用して、クロストークイベントを除去することができる。

本放射線検出装置は、放射線が入射するイベントによって生じる発光の減衰時定数が互いに異なる複数のシンチレータ素子が配列されたシンチレータ配列を複数空間的に異なる位置に配置したシンチレータ配列群と、前記シンチレータ群の前記シンチレータ配列の複数のシンチレータ素子からそれぞれ出力される光を受光し、シンチレータ素子毎に電気信号に変換する光電変換器と、前記イベントについて前記光電変換器からの電気信号の電圧波形におけるピーク値とを検出し、検出したピーク値との比に応じて、当該電気信号が前記シンチレータ配列中のどのシンチレータ素子に放射線が入射したことに起因する電気信号かを特定するとともに、放射線の入射した複数のシンチレータ素子を特定することで、放射線源の空間位置を検出するように構成されている演算処理装置と、を含む。

本装置によれば、複数のシンチレータ素子から入力する光を弁別することができる。従って、複数のシンチレータ素子を並べて、その一端からの光信号を前記光電変換器で電気信号に変換し、電圧波形におけるピーク値と積分電荷量の比に応じて、どのシンチレータ素子に放射線が入射したことに起因する電気信号かを正確に特定することが可能になる。例えば、ＤＯＩ検出器などにおいて効果的な信号の弁別を行え、深さ方向に配置するシンチレータ素子の数を大きくすることもできる。

シンチレーション光検出装置の構成の例を示す図である。 シンチレータ配列の構成の例を示す図である。 複数のシンチレータ配列を有する例を示す図である。 電圧波形の例を示す図である。 エネルギー別のイベント数の例を示す図である。 減衰時定数別のＶｐ／Ｑｔｏｔａｌに対するイベント数の例を示す図である。 減衰時定数別のエネルギー毎にＶｐ／Ｑｔｏｔａｌの例を示す図である。 エネルギー５００ｋｅＶ付近における減衰時定数別のＶｐ／Ｑｔｏｔａｌの例を示す図である。 材料が異なるシンチレータ素子を含むシンチレータ配列の構成の一例を示す図である。 ＧＳＯおよびＬｕＡＧについての減衰時定数別のＶｐ／Ｑｔｏｔａｌに対するイベント数の例を示す図である。 ＧＳＯおよびＬｕＡＧの光子数に対するＶｐ／Ｑｔｏｔａｌの例を示す図である。 エネルギーによるクロストークイベントの除去の例を示す図である。 クロストーク除去後の光子数に対するＶｐ／Ｑｔｏｔａｌの例を示す図である。 エネルギーによるクロストークイベント除去の例を示す図である。 生データのクロストークイベント数の例を示す図である。 エネルギーによるクロストークイベント除去後のイベント数の例を示す図である。 ＤＯＩ検出器の構成の例を示す図である。 閾値決定処理の他の例を説明するための図である。 決定された閾値を利用して弁別したイベントを示す図である。

以下、本シンチレーション光検出装置および放射線検出装置の実施形態について、図面に基づいて説明する。なお、本装置は、ここに記載される実施形態に限定されるものではない。

＜本シンチレーション光検出装置の概要＞
本シンチレーション光検出装置の概要について説明する。まず、本出願人は、先に特願２０１５−３５７８８号（先願）を出願している。この先願においては、ＬａＢｒ_３シンチレータを使用して、γ線を検出するに当たり、自己放射能によるα崩壊イベントを分離することを課題としている。そして、放射線入射による発光によって得られる信号の電圧波形（パルス波形）の相違に着目し、各イベントにおける電圧波形の電圧ピーク値Ｖｐと、積分電荷量（全電荷量）の比（Ｖｐ／Ｑｔｏｔａｌ）を算出する。そして、Ｖｐ／Ｑｔｏｔａｌを指標とすることで、α崩壊イベントを除外して、検出対象のγ線のイベントをカウントする。なお、積分電荷量は、電気信号である電圧Ｖを１イベント（１つのピーク）について時間積分した結果であり、電気信号の積分値として代用することができ、本実施形態では、積分電荷量として電気信号である電圧Ｖの積分値を用いる。

ここで、Ｖｐ／Ｑｔｏｔａｌによるイベントの弁別について、理論的考察を進めたところ、この弁別性能は、シンチレータによる発光の減衰時間に依存していることがわかった。すなわち、検出対象のγ線のイベントにおけるシンチレータの発光の減衰時間と、α崩壊のイベントにおける発光の減衰時間が異なり、これによってＶｐ／Ｑｔｏｔａｌが異なっている。

一方、ＰＥＴ検査装置において使用されるＧＳＯ（Ｇｄ_２ＳｉＯ_５）：Ｃｅシンチレータ（セシウム添加ＧＳＯシンチレータ）では、Ｃｅ濃度によって、発光の減衰時定数が変化する。例えば、文献（C. L. Melcher and J. S. Schweitzer "SCINTILLATION PROPERTIES OF GSO" IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 37, NO. 2, APRIL 1990）には、Ｃｅ濃度と減衰時定数の関係が示されている。この例では、Ｃｅ濃度：0.1, 0.5, 0.95, 1.5, 2.0, 5.0ｍｏｌ％としたときの減衰時定数が190, 56, 32, 27, 19nsであったとされている。

そこで、本シンチレーション光検出装置は、Ｃｅ濃度を変更したＧＳＯシンチレータ素子を複数用意し、これらから得られた電圧波形のＶｐ／Ｑｔｏｔａｌを調べることで、いずれのシンチレータ素子からの信号かを弁別する。

さらに、ＰＥＴ検査において使用されるシンチレータとして、ＬｕＡＧ（ルテニウム・アルミニウム・ガーネット）：Ｐｒシンチレータ（プラセオジム添加ＬｕＡＧシンチレータ）も知られており、Ｐｒの添加量により発光の減衰時定数が変化する。従って、これを利用することもできる。

さらに、ＧＳＯ（Ｇｄ_２ＳｉＯ_５）：Ｃｅシンチレータと、ＬｕＡＧ（ルテニウム・アルミニウム・ガーネット）：Ｐｒシンチレータとでは、発光の減衰時定数が異なる。そこで、両者を組み合わせて信号の弁別を行うことで、より多くのシンチレータ素子からの信号を弁別することが可能となる。なお、他のシンチレータ材料を用いることもできる。

ここで、ＰＥＴ検査において、その位置解像度を上昇するものとして、ＤＯＩ検出器が知られているが、ＤＯＩ検出器では深さ方向に配置された複数のシンチレータ素子からの信号を分離する。本シンチレーション光検出装置では、上述したように複数のシンチレータ素子からの信号を分離できる。そこで、本シンチレーション光検出装置によって、深さ方向に配置された複数のシンチレータ素子からの信号の分離を効果的に行うことで、より簡単な構成で位置分解能の高いＤＯＩ検出器を構成することができる。

なお、本シンチレーション光検出装置は、１つの光電子増倍管（光電変換器）に入射してくる、複数のシンチレータ素子からの光に基づく信号を弁別できるので、その適用はＤＯＩ検出器に限定されない。

＜全体構成＞
図１は、本実施形態におけるシンチレーション光検出装置の全体構成を示す図である。シンチレータ配列１０は、放射線の入射によって発光する複数のシンチレータ素子１２（１２ａ〜１２ｅ）が１列に配列して形成されている。シンチレータ配列１０の一端は、光電変換器である光電子増倍管１４に接続されている。光電子増倍管１４のシンチレータ配列１０が接続される部分には、受光面が形成され、シンチレータ素子１２ａ〜１２ｅで発光した光がここより入射する。

光電子増倍管１４によって、入射光が光電変換され、入射光の強度に応じた電気信号が出力される。光電子増倍管１４からの電気信号は電気波形取得装置１６に入力され、ここで時間軸に沿った電圧波形信号が得られる。なお、電気波形取得装置１６としては、オシロスコープなどが採用できる。

電気波形取得装置１６の出力は、演算処理装置１８を介し、記憶装置２０に記憶される。すなわち、電圧波形信号の生データ（波形整形などを施さないデータ（Ｒａｗ））がそのまま記憶装置２０に記憶される。

演算処理装置１８は、通常のコンピュータで形成され、電圧波形（電圧波形信号）についてのデータを解析し、解析結果を出力する。すなわち、演算処理装置１８は、所定時間に亘る１回の検査の終了時において、次のような解析処理を行う。
（ａ）電圧波形信号に対するフィルタ処理
（ｂ）積分電荷量Ｑｔｏｔａｌおよび電圧ピーク値Ｖｐを算出して、Ｖｐ／Ｑｔｏｔａｌを算出する処理
（ｃ）積分電荷量ＱｔｏｔａｌおよびＶｐ／Ｑｔｏｔａｌに対して設定された閾値を用いて、各電圧波形信号(電気信号）がいずれのシンチレータ素子からの信号かを特定する処理

ここで、図２には、シンチレータ配列１０と、光電子増倍管１４について示してある。シンチレータ配列１０は、５つのシンチレータ素子１２ａ〜１２ｅを１列に積層したものであって、シンチレータ素子１２ｅの底面側が光電子増倍管１４の受光面に接続されている。そこで、５つのシンチレータ素子１２ａ〜１２ｅのいずれに放射線が入射し、発光した場合にもその光は光電子増倍管１４に入力される。従って、光電子増倍管１４からの出力（電気信号）においては、いずれのシンチレータ素子１２ａ〜１２ｅへの放射線の入射に応じた発光に基づくものかはわからない。演算処理装置１８の解析によって、いずれのシンチレータ素子１２ａ〜１２ｅの発光かを特定する。

図３には、シンチレータ配列１０を複数個（４個）マトリクス状に配置してシンチレータブロックを形成した例を示す。すなわち、シンチレータ配列１０−１〜１０−４を光電子増倍管１４の受光面上にまとめて配置する。各シンチレータ配列１０−１〜１０−４は、一番上にシンチレータ素子１２−１ａ〜１２−４ａが位置し、一番下のシンチレータ素子１２−１ｅ〜１２−４ｅが光電子増倍管１４の受光面上に位置する。

図３の右下には、光電子増倍管１４の受光面を示してある。このように、シンチレータ配列１０−１〜１０−４のそれぞれに対応して、受光部１４−１〜１４−４が設けられており、対応するシンチレータ配列１０−１〜１０−４からの光が受光部１４−１〜１４−４からそれぞれ入射する。光電子増倍管１４には、４つのアノードに接続される信号線１４−１ｂ〜１４−４ｂを有しており、シンチレータ配列１０−１〜１０−４からの光に応じた信号が個別に電気波形取得装置１６に供給される。

なお、図３の例では、シンチレータブロックを４つのシンチレータ配列１０−１〜１０−４をマトリクス状に配置してブロックを構成したが、例えば５×５＝２５個や、それ以上でも構わない。

＜シンチレータ素子の構成＞
本実施形態においては、シンチレータ素子１２ａ〜１２ｅにおいて、発光の減衰時定数が異なるものを用いる。発光の時定数の変化は、電気波形取得装置１６で得られる電圧波形の減衰時定数に対応する。

例えば、ＰＥＴ用のシンチレータとして広く用いられている、ＧＳＯ（Ｇｄ_２ＳｉＯ_５）：Ｃｅ（セシウム添加ＧＳＯ）シンチレータでは、ＧＳＯの単結晶に賦活材としてＣｅを含み、そのＣｅ濃度によって、発光の減衰時定数が変化することが知られている。また、ＬｕＡＧ（ルテニウム・アルミニウム・ガーネット）：Ｐｒ（プラセオジム添加ＬｕＡＧ）シンチレータにおいても、ＬｕＡＧ結晶に対する賦活材Ｐｒの添加量によって発光の減衰時定数が変化する。なお、その他の材料を用いるシンチレータにおいても、セシウム等の賦活材の添加量により発光の減衰時定数を制御することが可能であり、減衰時定数が制御できればどのような材料を用いてもよい。

このように、賦活材の添加量を変更することで、減衰時定数が変化する。本実施形態では、このようにして減衰時定数が異なるシンチレータ素子１２ａ〜１２ｅを１列に配列し、いずれからの光も光電子増倍管１４に入射する。

＜演算部における処理＞
（ａ）フィルタ処理
電気波形取得装置１６から出力される電圧波形信号は、ノイズを含んでいるため、その後の処理（Ｖｐ／Ｑｔｏｔａｌの算出）のために、ノイズを除去することが好ましい。そこで、電圧波形信号について、例えば数ｎｓ（２ｎｓ）の時間幅で移動平均したり、５０ＭＨｚ以上の高周波ノイズをローパスフィルタで除去することで、ノイズを除去することが好ましい。

図４には、ＧＳＯ：Ｃｅシンチレータ素子において観測した１イベントについての電気波形取得装置１６から出力される電圧波形が示してある。縦軸は、規格化された強度（Normalized Intensity）であり、横軸が時間(Time[ns])である。減衰時定数３０，６０，９０ｎｓの３種類のシンチレータ素子についての波形を示してある。

（ｂ）Ｖｐ／Ｑｔｏｔａｌの算出
図４に示された電圧波形から、ピーク電圧Ｖｐおよび全電荷Ｑｔｏｔａｌが算出できる。すなわち、ピーク電圧は、電圧値の絶対値が最も大きな電圧であり、減衰時定数３０ｎｓ，６０ｎｓ，９０ｎｓで、それぞれ−０．００８，−０．０１，−０．０１６程度である。一方、全電荷Ｑｔｏｔａｌは電圧Ｖを１イベントについて時間積分した結果（積分電荷量）であり、図４の各曲線のピークの面積を積分することで得られる。従って、全電荷Ｑｔｏｔａｌは、規格化されたＶ＊ｎｓ（時間）のディメンジョンを有する。なお、ここでの電圧の時間積分値は、シンチレータへの放射線入射時（１イベント）の光量に応じたものであり、光電子増倍管１４において受け取った電荷量に対応するため、ここでは全電荷Ｑｔｏｔａｌとしている。すなわち、本実施形態では、光電子増倍管１４からの電流を電気波形取得装置１６で電流／電圧変換している。このため、この電流／電圧変換後の電圧の値は、入力インピーダンスＲ［Ω］によって変化し、Ｖ＝ＩＲ(例えばＲ＝５０Ω）の関係があり、この電圧は電流、その積分は電荷量に対応している。

本実施形態では、ＶｐとＱｔｏｔａｌの比（Ｖｐ／Ｑｔｏｔａｌ）を算出してこれをパラメータにすることが目的であり、電圧について規格化してあってもよく、またＱｔｏｔａｌについての単位もどのようなものでもよい。すなわち、Ｖｐ／Ｑｔｏｔａｌは、１イベントにおける（発光強度のピーク）／（トータル発光量）に対応した値であり、電圧波形でいえば、（電圧のピーク値）／（ピークの面積：積分値）が、Ｖｐ／Ｑｔｏｔａｌに対応する。

図５は、各種の減衰時定数のＧＳＯ（Ｇｄ_２ＳｉＯ_５）：Ｃｅシンチレータについて、観測イベント数をシミュレーションし、結果をエネルギーに対しプロットしたものである。Ｃｅ濃度に応じて減衰時定数が３０，４５，６０，８０，１００，１４０，１８０，２２０ｎｓの８種類に設定されたＧＳＯ：Ｃｅシンチレータについての結果であり、各シンチレータにおいて有意の差はなかった（図５において、シンチレータ別の結果は判別できない）。

図６は、減衰時定数３０ｎｓ，３６．８ｎｓ，４４．７ｎｓ，５６．５ｎｓ，７３．４ｎｓ，１００ｎｓ，１４５ｎｓのシンチレータ素子による観測結果をＶｐ／Ｑｔｏｔａｌごとのイベント数としてプロットしたものである。この結果から、７つの減衰時定数のシンチレータ素子からの信号を分離することが可能であることがわかる。

＜閾値の決定処理＞
図７は、ＧＳＯ（Ｇｄ_２ＳｉＯ_５）：Ｃｅシンチレータにおいて観測した各イベントについて、Ｖｐ／Ｑｔｏｔａｌを計算し、横軸をエネルギーとして、プロットしたものである。消滅γ線は５１１ｋｅＶのエネルギーを持っているため、これを含む領域に対して閾値を決定する。また、各減衰時定数のシンチレータ素子に応じたＶｐ／Ｑｔｏｔａｌ境界値を閾値として設定することにより、複数のシンチレータ素子からの信号を弁別できる。
このように、減衰時定数の相違によって、観測したイベントについてのＶｐ／Ｑｔｏｔａｌが変化し、減衰時定数に基づいて、イベントを分離可能であることがわかる。図７より理論的には７段階の分離が可能であり、測定器に起因する誤差を考慮しても、５段階の分離は確実に行えることがわかる。

図８には、エネルギー５１１ｋｅＶのγ線イベントについてのＶｐ／Ｑｔｏｔａｌの分布を示してある。図中「＜」で示した部分は、クロストークイベントが発生している可能性がある領域であり、同じ種類のシンチレータ材料で異なる減衰時定数をもつシンチレータ素子を用いたとしても、完全な弁別は困難である。

なお、閾値の設定は、減衰時定数の相違によるγ線イベントを識別できれば、どのようなものでもよく、多項式や、指数関数などの係数を適切なものに設定することで得ることもできる。

＜他の構成例＞
図９には、シンチレータ配列１０の他の構成例を示してある。この例では、１０個のシンチレータ素子１２１ａ，１２２ｂ，１２１ｃ，１２２ｄ，１２１ｅ，１２２ｆ，１２１ｇ，１２２ｈ，１２１ｉ，１２２ｊがこの順で積層されて構成されている。そして、シンチレータ素子１２１ａ，１２１ｃ，１２１ｅ，１２１ｇ，１２１ｉには、ＧＳＯ：ＣｅシンチレータでＣｅ濃度を変えて減衰時定数の相違するものを採用し、シンチレータ素子１２２ｂ，１２２ｄ，１２２ｆ，１２２ｈ，１２２ｊにはＬｕＡＧ：ＰｒシンチレータでＰｒ濃度を変えて減衰時定数の相違するものを採用している。なお、図においては、シンチレータ素子１２２の正面領域に斜線のハッチングを付し、２種類のシンチレータ素子１２１，１２２を区別して示してある。また、１つのシンチレータ素子１２１，１２２の大きさは一辺２．５ｍｍの立方体とした。

この例において、１０個のシンチレータ素子１２１ａ，１２２ｂ，１２１ｃ，１２２ｄ，１２１ｅ，１２２ｆ，１２１ｇ，１２２ｈ，１２１ｉ，１２２ｊ，１２１ａ，１２１ｃ，１２１ｅ，１２１ｇ，１２１ｉの減衰時定数は、それぞれ３０ｎｓ，１５ｎｓ，３６．８ｎｓ，１６．６ｎｓ，４４．７ｎｓ，１８．６ｎｓ，５６．５ｎｓ，２１．１ｎｓ，７３．４ｎｓ，２４．３ｎｓに設定されている。

このようなシンチレータ配列１０について、５１１ｋｅＶのγ線を入力して電圧波形を得るシミュレーションを行い、Ｖｐ／Ｑｔｏｔａｌを算出してプロットしたところ、図１０に示すように、ＧＳＯシンチレータと、ＬｕＡＧシンチレータとでは、Ｖｐ／Ｑｔｏｔａｌの値がずれている。従って、図９のシンチレータ配列１０において、１０個のシンチレータ素子からの信号弁別が可能であることがわかった。なお、図１０において縦軸はイベント数を規格化した強度（Intensity）としている。

図１１には、横軸に光子数（Number of photon）、縦軸にＶｐ／Ｑｔｏｔａｌをとったイベントのプロットを示してある。光子数は入射した消滅γ線のエネルギーに比例するが、シンチレータの特性によるため、５１１ｋｅＶのγ線を入射した場合において、ＧＳＯシンチレータと、ＬｕＡＧシンチレータとでは、光子数が異なっている。一方、観測データからは、ＧＳＯシンチレータと、ＬｕＡＧシンチレータの中間の光子数のイベントも出現しており、これらはクロストークイベントによるノイズと考えられる。そこで、図１２に四角で示すよう中間の光子数のイベントを削除する。これによって、図１３に示すように、クロストークが抑制された観測結果が得られる。

すなわち、ＧＳＯ：Ｃｅのシンチレータ素子と、ＬｕＡＧ：Ｐｒのシンチレータ素子とを交互に配置することによって、それらのクロストークによるイベントをそれが持つ発光量に応じて除去することができる。

図１４には、１０個のシンチレータ素子１２ａ〜１２ｊ（図においてａ〜ｊで示す）についてのイベントカウント数のプロットが示してある。黒丸が生データ、四角がエネルギー（光子数）閾値を利用してクロストークを排除した場合のデータ、白丸がシミュレーションによって求めたクロストークを含まないデータ（真のデータ）である。このように、エネルギー閾値を利用してクロストークを排除することによって、カウント数がほぼ真のデータに近づくことがわかる。なお、エネルギーの閾値についても、１種類の材料のシンチレータ素子の観測結果により決定するとよい。

図１５は、シミュレーションによって得られた生データにおけるイベントカウント数の内のクロストークイベントの割合を示し、図１６はエネルギー閾値を用いて弁別を行った場合のイベントカウント数の内のクロストークイベントの割合を示している。このように、生データではクロストークイベント数の割合が８．１６％あるが、エネルギー閾値を用いて弁別した場合には、クロストークイベント数の割合が０．８５％となり、クロストークイベントの大部分が除去されていることがわかる。

このように、２種類の異なる物質のシンチレータを利用することで、Ｖｐ／Ｑｔｏｔａｌが広い範囲に広がるため、より多くのシンチレータ素子からの信号の弁別が可能になる。さらに、２種類のシンチレータからなるシンチレータ素子を交互に配置することで、その出力エネルギー（光子数）の相違に基づき、その間のエネルギーとなるイベントを除去することでクロストークイベントを排除することができ、正しいカウント数に補正することが可能である。

＜ＰＥＴ検査装置＞
図１７には、ＰＥＴ検査装置（放射線検出装置）の模式図が示してある。患者が置かれる空間を取り囲むように、多数のシンチレータ配列１０が配置されたシンチレータ配列群が設けられる。患者の癌細胞に集積した陽電子放出核種に起因する５１１ｋｅＶの消滅γ線が放射線源となる。各シンチレータ配列１０にからの信号は光電子増倍管１４でそれぞれ別々に光電変換される。１または複数のシンチレータ配列１０に対応して１つの光電子増倍管１４を設け、各シンチレータ配列１０の各シンチレータ素子１２から電気信号が別々に得られる。得られた各シンチレータ素子からの電気信号が電気波形取得装置１６にてそれぞれ電圧波形信号に変換され、得られた各シンチレータ素子１２からの電圧波形信号が演算処理装置１８に供給される。演算処理装置１８は、供給される各電圧波形信号からいずれのシンチレータ素子１２からの信号であるかを特定するとともに、同時に計測された一対のシンチレータ素子１２を特定する。さらに、経時的に得られる複数の一対のシンチレータ素子１２による一対のシンチレータ素子１２の特定結果に基づき、消滅γ線が発生された癌細胞の３次元位置（空間位置）を特定する。すなわち、一対のシンチレータ素子１２によって形成される複数の直線の交点を特定することで、空間位置を特定することができる。なお、ＰＥＴ検査装置の各構成部材について、上述したような構成を採ることができる。

＜その他の構成＞
Ｖｐ／Ｑｔｏｔａｌの代わりにＱｔｏｔａｌ／Ｖｐを用いても同様のことが行える。

シンチレータ配列は光電子増倍管から一方向に伸びる構成としたが、これに限らず光電子増倍管の１つの受光面に複数のシンチレータ素子を配列して、ここからの光を光電子増倍管に入射するようにしてもよい。これによって、１つの光電子増倍管からの出力を処理することで、複数のシンチレータ素子のいずれに放射線が入射したのかを検出することができる。

シンチレータ素子の発光の減衰時定数は、シンチレータ素子における賦活材（例えば、ＧＳＯ：ＣｅのＣｅ）の濃度を変更することで、制御することができる。また、シンチレータ素子の材料（ＧＳＯ，ＬｕＡＧ）を変更することでも制御できる。この両方を利用することで、弁別できるシンチレータ素子の数を大きくすることができる。シンチレータ素子の材料については、３種類以上利用してもよい。

また、シンチレータ素子の材料（ＧＳＯ，ＬｕＡＧ）が異なる場合、１イベントにおける発生光子数（エネルギー）が異なる。そこで、このエネルギー差により、いずれのシンチレータ素子に放射線が入射したのかを弁別することも可能である。

＜閾値の決定処理の他の例＞
図７においては、エネルギーに対し、Ｖｐ／Ｑｔｏｔａｌをプロットし、５１１ｋｅＶ近辺において、遅延時間が異なるＶｐ／Ｑｔｏｔａｌを分離できるように閾値を決定した。シンチレータを積層する場合、光電子増倍管から遠いシンチレータほど光量が低下し、同一エネルギーのガンマ線に対する信号でも、シンチレータの違いにより光量が変化するため、エネルギーの代わりにＱｔｏｔａｌを用いて説明する。なお、Ｑｔｏｔａｌを単にＱと表記する。

図１８には、Ｃｅ濃度の異なる４つのＧＳＯシンチレータ素子を用いたイベントの計測結果を示してある。Ｃｅ濃度の異なる４つのＧＳＯシンチレータ素子として、ＧＳＯ：Ｃｅ Ａ ｍｏｌ％、ＧＳＯ：Ｃｅ Ｂ ｍｏｌ％、ＧＳＯ：Ｃｅ Ｃ ｍｏｌ％、ＧＳＯ：Ｃｅ Ｄ ｍｏｌ％を用意した。ここで、Ｃｅ濃度を示す、Ａ，Ｂ，Ｃ，Ｄは、０．２５ｍｏｌ％〜１．７５ｍｏｌ％程度であり、Ａ＞Ｂ＞Ｃ＞Ｄである。

図１８（Ａ）に示すように、全電荷Ｑ（Ｉｎｔｅｇｒａｔｅｄ ｃｈａｒｇｅ）に対し、Ｖｐ／Ｑをプロットしてある。全電荷Ｑｔｏｔａｌ（Ｑ）の３−４［Ｖ＊ｎｓ］の範囲では、比較的４つのＧＳＯシンチレータの検出結果が離れている。Ｖｐ／Ｑ値は、基本的にガンマ線のエネルギーに依らず一定となることから、この範囲のイベントをピックアップし、図１８（Ｂ）に示すようなＶｐ／Ｑ値に対するヒストグラムを作成することにより、詳細なシンチレータごとのＶｐ／Ｑ閾値を設定することができる。この図から隣り合うピークの中心の中央を各シンチレータを弁別するための閾値を設定することができる。

このようにして得られた閾値によって分離されたＱに対するＶｐ／Ｑ値の分布を図１８（Ｃ）に示す。このように、４つのＧＳＯシンチレータ素子の検出結果を分離できることがわかる。

そして、上述のようにして決定された閾値毎に測定データを分割した上で、エネルギースペクトルを作成した結果を図１９に示す。点線は、５１１ｋｅＶのイベントをフィッティングした結果である。図１９（Ａ）はＧＳＯ：Ｃｅ Ａ ｍｏｌ％、図１９（Ｂ）はＧＳＯ：Ｃｅ Ｂ ｍｏｌ％、図１９（Ｃ）はＧＳＯ：Ｃｅ Ｃ ｍｏｌ％、図１９（Ｄ）はＧＳＯ：Ｃｅ Ｄ ｍｏｌ％における全電荷Ｑに対するイベント数を示す。

なお、図７における５１１ｋｅＶのイベントは０．５ＭｅＶ付近にある円状の集団であり、図１８（Ｃ）においても、シンチレータごとに分離されたデータの中央付近にある楕円状の集団である。このように、４つのＧＳＯシンチレータ素子のイベントを明確に分離することができる。

＜実施形態の効果＞
本実施形態によれば、発光の減衰時定数が相違するシンチレータ素子を用いることで、光電変換器で得られた複数のシンチレータ素子からの信号を分離できる。そこで、ＤＯＩ検出器における信号の弁別（放射線入射シンチレータの弁別）を効果的に行うことができる。

シンチレータ素子を１列に配置したシンチレータ配列を用い、各シンチレータからの光を１つの光電変換器に入射することで、ＤＯＩ検出器を容易に構成することが可能となる。

さらに、シンチレータ配列を複数まとめてブロックを形成し、光電変換器に個別の受光面から入射させることができる。

シンチレータ素子において、賦活材の濃度を変更することで、シンチレータ素子の減衰時定数を容易に変更することができる。

シンチレータ素子において、材料を変更することで、シンチレータ素子の減衰時定数を変更することもできる。

シンチレータ素子において、賦活材の濃度および材料の両方を変更することで、より多くのシンチレータ素子からの信号の弁別が可能になる。

実施形態に係るシンチレーション光検出装置を利用することで、ＤＯＩ検出器における信号の弁別を効果的に行え、ＰＥＴ検査を効果的に行うことができる。

１０ シンチレータ配列、１２，１２１，１２２ シンチレータ素子、１４ 光電子増倍管、１６ 電気波形取得装置、１８ 演算処理装置、２０ 記憶装置。

Claims (9)

1. 放射線が入射するイベントによって生じる発光の減衰時定数が互いに異なる複数のシンチレータ素子が配列されたシンチレータ配列と、
   前記シンチレータ配列の複数のシンチレータ素子から出力される光を受光し、電気信号に変換する光電変換器と、
   前記イベントについて前記光電変換器からの電気信号の電圧波形におけるピーク値と積分電荷量を検出し、検出したピーク値と積分電荷量の比に応じて、当該電気信号が前記シンチレータ配列中のどのシンチレータ素子に放射線が入射したことに起因する電気信号かを特定するように構成されている演算処理装置と、
   を含む、
   シンチレーション光検出装置。
2. 請求項１に記載のシンチレーション光検出装置であって、
   前記シンチレータ配列は、複数のシンチレータ素子が１列に配列され、当該１列の一端のシンチレータ素子から複数のシンチレータ素子からの光が出力される、
   シンチレーション光検出装置。
3. 請求項２に記載のシンチレーション光検出装置であって、
   前記シンチレータ配列を複数配列し、
   前記光電変換器は、個別のシンチレータ配列からの光を別々に受光して、別々の電気信号を出力し、
   前記演算処理装置は、各電気信号の電圧波形におけるピーク値と積分電荷量の比に応じて、各電気信号がどのシンチレータ配列のどのシンチレータ素子に放射線が入射したことに起因する電気信号かを特定する、
   シンチレーション光検出装置。
4. 請求項１〜３のいずれか１つに記載のシンチレーション光検出装置であって、
   前記減衰時定数が互いに異なる複数のシンチレータ素子は、賦活材の濃度が互い異なるものを含む、
   シンチレーション光検出装置。
5. 請求項１〜３のいずれか１つに記載のシンチレーション光検出装置であって、
   前記減衰時定数が互いに異なる複数のシンチレータ素子は、シンチレータ素子を構成する材料が互いに異なるものを含む、
   シンチレーション光検出装置。
6. 請求項１〜３のいずれか１つに記載のシンチレーション光検出装置であって、
   前記減衰時定数が互いに異なる複数のシンチレータ素子は、賦活材の濃度が互いに異なるとともに、シンチレータ素子を構成する材料が異なるものを含む、
   シンチレーション光検出装置。
7. 請求項５に記載のシンチレーション光検出装置であって、
   前記演算処理装置は、前記材料が異なる複数のシンチレータ素子についての出力エネルギーおよび各電気信号の電圧波形におけるピーク値と積分電荷量の比の相違を利用して、クロストークイベントを除去する、
   シンチレーション光検出装置。
8. 請求項６に記載のシンチレーション光検出装置であって、
   前記演算処理装置は、前記材料が異なる複数のシンチレータ素子についての出力エネルギーの相違を利用して、クロストークイベントを除去する、
   シンチレーション光検出装置。
9. 放射線が入射するイベントによって生じる発光の減衰時定数が互いに異なる複数のシンチレータ素子が配列されたシンチレータ配列を複数空間的に異なる位置に配置したシンチレータ配列群と、
   前記シンチレータ群の前記シンチレータ配列の複数のシンチレータ素子からそれぞれ出力される光を受光し、シンチレータ素子毎に電気信号に変換する光電変換器と、
   前記イベントについて前記光電変換器からの電気信号の電圧波形におけるピーク値と積分電荷量の比に応じて、当該電気信号が前記シンチレータ配列中のどのシンチレータ素子に放射線が入射したことに起因する電気信号かを特定するとともに、放射線の入射した複数のシンチレータ素子を特定することで、放射線源の空間位置を検出するように構成されている演算処理装置と、
   を含む、
   放射線検出装置。
   

Images