(57)【要約】【課題】  従来の技術では中性子の入射方向を特定するために複数の検出器を配置するが、入射方向を特定する原理上、該検出器を配置する位置が決定あるいは制限されており、配置に対する自由度が無かった。これは検出器を設計する上で外観や形状に対する制限が大きく、特に外観形状の変化に十分な冗長性が求められるウェアラブル型の検出器には採用が困難であった。【解決手段】  複数の中性子検出部を、中性子を減速する作用を有する人体、含水体、ポリエチレン等（以下、減速体という）に装着して用い、該検出部の計数を比較することによって、中性子線源の方向を特定できる。また複数の中性子検出部は減速体に装着する際に互いに重ならなければ該検出部の配置を制限されないことから布等のフレキシブルな素材に取り付けられることで、外観形状の変化に十分な冗長性を持たせることが可能となる。【選択図】  図４
2018013378.tif 000002
【特許請求の範囲】
【請求項１】
装着して用いる中性子検出器であって、装着時には異なる位置に位置することになる複数の中性子検出部と、該検出部で検出された信号の計数を比較する手段と、該比較結果から中性子線源の方向を特定する手段を具備する装着型中性子検出器。
【請求項２】
中性子検出部がフレキシブルな布状素材上に取り付けられてなる請求項１記載の装着型中性子検出器。
【請求項３】
フレキシブルな布状素材が、上衣形状をしている請求項２記載の装着型中性子検出器。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は中性子シンチレーターを用いる中性子検出器に関する。詳しくは中性子線源の方向を特定する中性子検出器に関する。
【背景技術】
【０００２】
放射線の入射方向を検出する技術として特許文献１〜３に記載にされる従来技術がある。特許文献１に記載される技術は、特にガンマ線ではあるが放射線の入射方向を検出するため、同じ材質の独立した柱状の複数本のシンチレーターを一本の柱を形成するように束ねて配置している。複数本のシンチレーターが入射する放射線に対して互いに影を形成するように並べて配置することで、各シンチレーターに対して直接入射する放射線と他のシンチレーターの影になって間接的に入射する放射線の割合の組み合わせが放射線の入射方向によって変化するようにされている。
【０００３】
特許文献２に記載される技術は、中性子の入射方向を検出するため４つの無指向性中性子検出器は前後と左右で対象性を持つように密に配置され、前後の該検出器の中性子の計数差と左右の該検出器の中性子の計数差から合成されるベクトルから中性子の入射方向を求めている。
【０００４】
特許文献３に記載される技術は、６つの中性子検出器を規則正しくH型に並べることで、入射する中性子に向いた検出器とその検出器に遮蔽された検出器の各々の計数率の差から角度分解能を得ている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００８−１３４２００号公報
【特許文献２】米国特許番号５３４５０８４号
【特許文献３】米国特許番号５６５９１７７号
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、特許文献１〜３の技術では中性子の入射方向を特定するために複数の検出器を配置するが、入射方向を特定する原理上、該検出器を配置する位置が決定（制限）されており、配置に対する自由度が無かった。これは検出器を設計する上で外観や形状に対する制限が大きい。特に外観形状の変化に十分な冗長性が求められるウェアラブル型の検出器には採用が困難である。
【０００７】
該ウェアラブル型検出器は特にテロリストが所持、あるいは秘匿、隠蔽した核物質を探知する場合に有効性を発揮する。従来の検出器は検出器本体をバックパック等に入れて背負い、核物質を探しながら巡回して使用している。この場合、検出器を所持していることが第３者から見て分かり易く、目立ち易いため、ひいてはテロリストの警戒を招く可能性がある。ウェアラブル型の検出器であれば装着して使用することで外見から検出器を所持していることが第３者に分かり難くすることが出来、ひいてはテロリストを警戒させずに核物質を探索することが可能になる。
【０００８】
本発明は上記のような問題点を解決することを目的になされた。中性子の入射方向を特定でき検出器の配置に制限の少ない中性子検出器を提供する。
【課題を解決するための手段】
【０００９】
本発明者は中性子の入射方向を特定できる検出器について鋭意検討を重ねた。その結果、複数の中性子検出部を、中性子を減速する作用を有する人体、含水体、ポリエチレン等（以下、減速体という）に装着して用い、該検出部の計数を比較することによって、中性子線源の方向を特定できることを見出した。さらに、複数の中性子検出部は減速体に装着する際に互い重ならなければ配置を制限されないことを見出した。
【００１０】
即ち、本発明は、装着して用いる中性子検出器であって、装着時には異なる位置に位置することになる複数の中性子検出部と、該検出部で検出された信号の計数を比較する手段と、該比較結果から中性子線源の方向を特定する手段を具備する装着型中性子検出器である。
【００１１】
また本発明は、中性子検出部がフレキシブルな布状素材上に取り付けられてなる装着型中性子検出器であって、フレキシブルな布状素材は上衣形状をしていることが好ましい。
【発明の効果】
【００１２】
本発明によると中性子検出器により中性子線源の方向を特定でき、中性子検出部の配置に制限が少ないことから検出器の外観形状の変化に十分な冗長性を持たせることが可能な中性子検出器を提供出来る。
【図面の簡単な説明】
【００１３】
【図１】検出部が中性子を検出する様子を表した図である。
【図２】２個の検出部の減速体への配置と、減速体と線源の位置関係を表した図である。
【図３】３個の検出部の減速体への配置と、減速体と線源の位置関係を表した図である。
【図４】４個の検出部の減速体への配置と、減速体と線源の位置関係を表した図である。
【図５】３個の検出部を３等分した区域に１個ずつ配置した図である。
【図６】３個の検出部を等間隔に配置した図である。
【図７】３個の検出部で方位角の強度比を比較したスペクトルである。
【図８】４個の検出部で仰俯角の強度比を比較したスペクトルである。
【図９】実施例１~３の４個の検出部の減速体への配置を表した図である。
【図１０】実施例１~３の減速体と線源の位置関係を表した図である。
【図１１】実施例１で得られた波高分布スペクトルである。
【図１２】実施例１の線源の方向を特定するため方位角に対する最小二乗法で求めた誤差を表した図である。
【図１３】図１２の方位角０°付近を拡大表示した図である。
【図１４】実施例２で得られた波高分布スペクトルである。
【図１５】実施例２の線源の方向を特定するため方位角に対する最小二乗法で求めた誤差を表した図である。
【図１６】図１５の方位角９０°付近を拡大表示した図である。
【図１７】実施例３で得られた波高分布スペクトルである。
【図１８】実施例３の線源の方向を特定するため方位角に対する最小二乗法で求めた誤差を表した図である。
【図１９】図１８の方位角４５°付近を拡大表示した図である。
【発明を実施するための形態】
【００１４】
以下、本発明に係る中性子検出器の実施形態を説明する。ただし、本発明はここで取り挙げた実施形態に限定されるものではなく、その発明の技術的思想を逸脱しない範囲で適宜組み合わせや改良が可能である。
【００１５】
本発明の検出部は、減速されてエネルギーが低下した中性子、所謂熱中性子を検出する。そのため線源から飛来する中性子のうち直接検出部に飛来した中性子は減速する前でエネルギーが低下していないため検出されない。図１に検出部が中性子を検出する様子を示す。
【００１６】
線源から飛来した中性子のうち減速体に入射した中性子は減速体を構成する物質と衝突しエネルギーを衝突相手に与えながら減速する。中性子は減速されながら減速体内部を通過・散乱して、減速体に装着された検出部に到達し検出される。ここで、減速体は、中性子を減速すると同時に遮蔽もする。したがって、検出部に到達するまでに減速体内部を通過する距離が長いほど、中性子の計数は低下する。本発明は、この現象を応用し、複数の中性子検出部を減速体に装着し、該複数の検出部の計数を比較することによって、中性子線源の方向を特定する。
【００１７】
以下、本発明の中性子検出器の作用機序を、図２~４を用いて具体的に説明する。なお、便宜的に水平面上に検出器を設置し、方位角と仰俯角を特定する態様について説明するが、本発明の中性子検出器の使用方法はかかる態様に限定されない。
【００１８】
例えば図２に示すように、減速作用を有する人体の胸部に２個の検出部を装着する。正面の位置（Ａ）に線源がある場合には、線源と２個の検出部の位置が対称的であるため２個の検出部の計数は同じである。一方、側面（Ｃ）の位置に線源がある場合には、中性子が検出部aに到達するまでに減速体を通過する距離は、検出部ｂに到達するまでに減速体を通過する距離に比較して短い。
【００１９】
したがって、検出部ａの計数は検出部ｂの計数に比較して大きくなる。また、斜め前方の位置(Ｂ)に線源がある場合にも、検出部ａの計数は検出部ｂの計数に比較して大きくなるが、検出部ｂの計数に対する検出部ａの計数の比（ａ／ｂ）は、Ｃの位置の場合に比較して小さくなる。したがって、予め各方位角に線源がある場合の検出部ａの計数と検出部ｂの計数の比較結果を求めておき、実際の計測の際に得られた比較結果と照合することによって、線源の方向を特定することができる。
【００２０】
なお、図２に示した態様では、斜め前方の位置(Ｂ)に線源がある場合と斜め後方の位置(Ｂ’)に線源がある場合とでは検出部ａの計数と検出部ｂの計数の比較結果が類似するため、線源の方向を特定することが困難となる。したがって、図３に示すように検出部を３個以上用いることが好ましい。
【００２１】
図３に示す態様では、胸部に加えて背部に検出部ｃを装着している。斜め前方の位置(Ｂ)に線源がある場合に中性子が検出部ｃに到達するまでに減速体を通過する距離は、斜め後方の位置(Ｂ’)に線源がある場合と比較して長い。したがって、検出部aおよびｂの計数に対する検出部ｃの計数の比（それぞれｃ／ａおよびｃ／ｂ）は、線源の位置が斜め前方の位置(Ｂ)にある場合に、斜め後方の位置(Ｂ’)にある場合と比較して小さくなる。
【００２２】
したがって、予め各方位角に線源がある場合の検出部ａ，ｂ及びｃの計数の相互の比較結果を求めておき、実際の計測の際に得られた比較結果と照合することによって、線源の方向を特定することができる。よって、中性子検出部が３個以上あれば中性子の入射方向の方位角を３６０°の方向で特定することができる。
【００２３】
さらに図４に示す態様では、胸部に２個の検出部を地上から同じ高さで装着し、背部に２個の検出部を地上から異なる高さで装着している。これによって中性子線源の方位角を特定したのと同様の原理で仰俯角を特定することが出来る。水平の位置（Ｅ）に線源がある場合には、線源と検出部ｆと検出部ｇの位置が対称的であるため2個の検出部の計数は同じである。
【００２４】
一方、頭上（Ｇ）の位置に線源がある場合には、中性子が検出部ｆに到達するまでに減速体を通過する距離は、検出部ｇに到達するまでに減速体を通過する距離に比較して短い。したがって、検出部ｆの計数は検出部ｇの計数に比較して大きくなる。
【００２５】
また、斜め上方の位置(Ｆ)に線源がある場合にも、検出部ｆの計数は検出部ｇの計数に比較して大きくなるが、検出部ｇの計数に対する検出部ｆの計数の比（ｆ／ｇ）は、Ｇの位置の場合に比較して小さくなる。したがって、予め各仰俯角に線源がある場合の検出部ｆの計数と検出部ｇの計数の比較結果を求めておき、実際の計測の際に得られた比較結果と照合することによって、線源の仰俯角を特定することができる。
【００２６】
なお、斜め上方の位置(Ｆ)に線源がある場合と斜め上方の位置(Ｆ’)に線源がある場合とでは検出部ｆの計数と検出部ｇの計数の比較結果が類似するため、該比較結果から線源の方向を特定することが困難となるが、該比較結果に加えて検出部ｄならびに検出部ｅとの比較結果を組み合わせることで特定できる。即ち、斜め上方の位置(Ｆ)に線源がある場合に中性子が検出部ｄならびに検出部ｅに到達するまでに減速体を通過する距離は、斜め上方の位置(Ｆ’)に線源がある場合と比較して短い。
【００２７】
したがって、検出部fならびに検出部gの計数に対する検出部ｄならびに検出部ｅの計数の比は、線源の位置が斜め上方の位置(Ｆ)にある場合、斜め上方の位置(Ｆ’)にある場合と比較して大きい。したがって、予め各方向に線源がある場合の検出部d、e、f及びgの計数の相互の比較結果を求めておき、実際の計測の際に得られた比較結果と照合することによって、線源の方向を特定することができる。よって、中性子検出部が４個以上あれば、方位角に加えて仰俯角を特定することができ、全球型の検出器とすることができる。
【００２８】
本発明による中性子検出器の中性子検出部には無機蛍光体と樹脂とを含有する樹脂組成物からなる中性子シンチレーターや、ＬｉＦ／ＺｎＳ：Ａｇ，Ｌｉガラス（サンゴバン社製、ＧＳ−２０等)、Ｃｅ：Ｃｓ２ＬｉＹＣｌ６等の中性子シンチレーターを用いたシンチレーション検出器、Ｈｅ−３比例計数管、ＢＦ３比例計数管、ホウ素コート比例計数管等の比例計数管等を特に制限なく用いることができる。
【００２９】
シンチレーション検出器とは、シンチレーターと光検出器を接続したもので、入射した中性子の数に応じてシンチレーターが蛍光を発し、蛍光を検知した光検出器が光を電子に変換し増幅する。そしてパルス状の信号として出力し、該パルス状信号の強度（波高値）によって中性子の数を計測している。
【００３０】
比例計数管とは、ガス中での電子のなだれ増殖を利用した検出器である。円筒形の容器の中に不活性ガスをつめ、中に張った芯線に高電圧をかけることで計数管として動作させる。中性子の入射によってガスが電離され、中性子のエネルギーに比例した数の一次電子雲が作られる。電子は電場に引かれ芯線に向かって加速と電離を繰り返しながら増幅する。電子は芯線に到達すると電気パルスを発生し、該電気パルスの強度（波高値）と頻度によって中性子の数を計測している。
【００３１】
中でも前記樹脂組成物からなる中性子シンチレーターを中性子検出部に用いることによって、特に形状自由度の高い中性子検出器を得ることができるため好ましい。該中性子シンチレーターは例えば特開ＷＯ２０１４／０９２２０２に開示されており、リチウム６及びホウ素１０から選ばれる少なくとも１種の中性子捕獲同位体を含む無機蛍光体と樹脂とを含有する樹脂組成物からなる中性子シンチレーターである。樹脂組成物からなるため柔軟性があり外力による変形が容易であるため減速体の形状に合わせて中性子検出部を装着することができる。
【００３２】
中性子検出器の中性子検出部は減速体に装着して用いる。減速体は中性子を減速する作用を有していれば構成成分に制限はない。また、減速体の形状に制限はなく、柱状、球状、錐体状、や凹凸のある形状でも使用できる。なお、線源の方向を特定する精度を高めるため、計測対象とする各方向に線源を置いた場合に、全ての検出部で中性子を計数可能なように減速体の形状を調整することが好ましい。
【００３３】
中性子検出部は互いに重ならなければ配置は制限されないが、減速体を中心とする３６０°方向を、検出部の数で等角度分割した区域を設定し、各区域に１個ずつ中性子検出部を配置することが好ましい。即ち、例えば図５に示すように中性子検出部の数が３個の場合であれば、減速体を中心とする３６０°方向を、3個の区域に等角度分割（１２０°毎）し、各区域に１個ずつ中性子検出部を配置すれば良い。このように配置することで全方位に偏りなく検出感度を持たせることが出来る。より好ましくは図６に示すように減速体を中心として中性子検出部を等角度間隔に配置することで、全方位に特に偏りなく検出感度を持たせることが出来る。
【００３４】
本発明において、中性子検出器を装着するとは中性子検出部を減速体の表面近傍に保持することをいう。中性子検出部と減速体表面の距離は特に制限されないが、中性子の入射方向を精度良く特定するため、中性子検出部と減速体表面の距離を短くすることが好ましい。当該中性子検出部と減速体表面の距離は、一般には０〜２０ｃｍの範囲内であり、好ましくは０〜１０ｃｍ、より好ましくは０〜５ｃｍである。
【００３５】
計数を比較する手段は特に制限されず、個々の検出部について従来公知の計数法を採用し、個々の検出部の計数を求めてお互いに比較すれば良い。具体的な計数法として、前記シンチレーション検出器或いは比例計数管等の検出部から出力される出力信号を、増幅器等を介して増幅した後、所定の閾値を設けたコンパレータ等に入力して該閾値を超える信号の頻度を計数する方法が好適に採用される。また、各検出部に前記計数法を行う機器をそれぞれ接続し、各検出部の計数を得て、該計数をお互いに比較する手段が好適に採用される。
【００３６】
以下、計数の比較結果から中性子線源の方向を特定する手段について具体的に説明するが、以下に説明する手段に限定されない。
【００３７】
例えば図３の態様において、各方向に中性子線源を置いた場合に得られる各検出部の計数の比較結果を図７に示す。該比較結果と実際の計測の際に得られた比較結果とを照合することによって、線源の方向を特定することができる。
【００３８】
照合を簡便にするため、予め各方向に中性子線源を置いた場合に得られる各検出部の計数を応答関数で近似して表すことが好ましい。前記図７で得られた各検出部の計数は、例えば下記式（１）の応答関数で近似して表すことができる。
【００３９】
Ｎｉ＝ａ・ｓｉｎ(θ+αi)＋ｂ      （１）
（ただし、Ｎｉはｉ番目の検出部の計数、θは線源の方位角を表し、αiはi番目の検出部の装着位置に基づく定数であり、ａおよびｂは定数である。）
【００４０】
前記式（１）から得られる各θにおける計算値（Ｎｉ）と、実際の計測で得られたｉ番目の検出部の計数の実測値（Ｎ’ｉ）を照合し、最小二乗法によって最も誤差が小さくなるθを求め、線源の方向として特定することができる。
なお、前記応答関数は、用いる検出部の種類や減速体の形状等によって異なるため、所望の態様に応じて適宜求めることが好ましい。応答関数としては、前記三角関数の他、例えば下記の多項式を好適に用いることができる。
【００４１】
Ｎｉ＝ａｎ・(θ+αi)ｎ＋ａｎ−１・(θ+αi)ｎ−１＋・・・＋ａ１・(θ+αi)＋ａ０
（ただし、Ｎｉはｉ番目の検出部の計数、θは線源の方位角を表し、αiはi番目の検出部の装着位置に基づく定数であり、ａ０~ａｎは定数である。）
次に計数の比較結果から中性子線源の仰俯角を特定する手段について具体的に説明するが、以下に説明する手段に限定されない。
【００４２】
例えば図４の態様において、各仰俯角に中性子線源を置いた場合に得られる各検出部の計数の比較結果を図８に示す。該比較結果と実際の計測の際に得られた比較結果とを照合することによって、線源の仰俯角を特定することができる。
【００４３】
照合を簡便にするため、予め各方向に中性子線源を置いた場合に得られる各検出部の計数を応答関数で近似して表すことが好ましい。前記図７で得られた各検出部の計数は、例えば下記式（１）の応答関数で近似して表すことができる。
【００４４】
Ｎｉ＝ａ・ｓｉｎ(θ+αi)＋ｂ      （１）
（ただし、Ｎｉはｉ番目の検出部の計数、θは線源の方位角を表し、αiはi番目の検出部の装着位置に基づく定数であり、ａおよびｂは定数である。）
【００４５】
前記式（１）から得られる各θにおける計算値（Ｎｉ）と、実際の計測で得られたｉ番目の検出部の計数の実測値（Ｎ’ｉ）を照合し、最小二乗法によって最も誤差が小さくなるθを求め、線源の方向として特定することができる。
【００４６】
なお、前記応答関数は、用いる検出部の種類や減速体の形状等によって異なるため、所望の態様に応じて適宜求めることが好ましい。応答関数としては、前記三角関数の他、例えば下記の多項式を好適に用いることができる。
【００４７】
Ｎｉ＝ａｎ・(θ+αi)ｎ＋ａｎ−１・(θ+αi)ｎ−１＋・・・＋ａ１・(θ+αi)＋ａ０
（ただし、Ｎｉはｉ番目の検出部の計数、θは線源の方位角を表し、αiはi番目の検出部の装着位置に基づく定数であり、ａ０~ａｎは定数である。）
【００４８】
次に計数の比較結果から中性子線源の仰俯角を特定する手段について具体的に説明するが、以下に説明する手段に限定されない。
【００４９】
例えば図４の態様において、各仰俯角に中性子線源を置いた場合に得られる各検出部の計数の比較結果を図８に示す。該比較結果と実際の計測の際に得られた比較結果とを照合することによって、線源の仰俯角を特定することができる。
【００５０】
照合を簡便にするため、予め各仰俯角に中性子線源を置いた場合に得られる各検出部の計数を応答関数で近似して表すことが好ましい。前記図８で得られた各検出部の計数は、例えば下式の応答関数で近似して表すことができる。
【００５１】
Ｎｉ＝ａ・ｃｏｓ(θ+αi)＋ｂ
（ただし、Ｎｉはｉ番目の検出部の計数、θは線源の仰俯角を表し、αiはi番目の検出部の装着位置に基づく定数であり、ａおよびｂは定数である。）
【００５２】
前記式から得られる各θにおける計算値（Ｎｉ）と、実際の計測で得られたｉ番目の検出部の計数の実測値（Ｎ’ｉ）を照合し、最小二乗法によって最も誤差が小さくなるθを求め、線源の仰俯角として特定することができる。
【００５３】
なお、前記応答関数は、用いる検出部の種類や減速体の形状等によって異なるため、所望の態様に応じて適宜求めることが好ましい。応答関数としては、前記三角関数の他、例えば下記の多項式を好適に用いることができる。
【００５４】
Ｎｉ＝ａｎ・(θ+αi)ｎ＋ａｎ−１・(θ+αi)ｎ−１＋・・・＋ａ１・(θ+αi)＋ａ０
（ただし、Ｎｉはｉ番目の検出部の計数、θは線源の仰俯角を表し、αiはi番目の検出部の装着位置に基づく定数であり、ａ０~ａｎは定数である。）
【００５５】
本発明の装着型中性子検出器の形状は、装着対象の形状に従って適宜設定すればよいが、装着を簡便に行うため、或いは中性子検出部の柔軟性及び可撓性を維持しながら機械的な強度を高めるため、中性子検出部が布等のフレキシブルな素材に取り付けられていることが好ましい。また、人体に装着する際の操作性を鑑みて、該布がシャツやベスト等の上衣の形状を有していることが特に好ましい。
【００５６】
なお、中性子検出部の取り付け位置は上述した通りであるが、その他の装置の取り付け位置は特に限定されない。例えば特定された中性子線源の位置を装着者が直接確認できるよう、表示機器を具備することが好ましく、また、操作性を鑑みて手の平に収まる程度の大きさの無線の表示機器を具備することが特に好ましい。
【実施例１】
【００５７】
上述した本発明の中性子入射方向検出器の実施例について説明する。  減速体として高密度ポリエチレンを用意した。中性子シンチレーターには前記樹脂組成物からなる中性子シンチレーターを使用し、光検出器と接続して中性子検出部とした。中性子検出部を４個用意し、減速体に装着した。
【００５８】
なお、減速体を中心とする３６０°方向を、4個の区域に等角度分割（９０°毎）し、各区域に１個ずつ中性子検出部を配置した。具体的には図９のように、検出部ａと検出部ｂ、検出部ｃと検出部ｄの間の角が１２０°となるように、検出部ａと検出部ｄ、検出部ｂと検出部ｃの間の角が６０°となるように装着した。ここで減速体を検出部ａと検出部ｂを同時に見る側を基準の０°とし、検出部ａと検出部ｄを同時に見る側を９０°、またその間を二等分する位置を４５°とした。
【００５９】
中性子線源を図１０のように減速体からみて０°の方向となる線源位置（Ａ）に配置し中性子を検出したところ、各中性子検出部から図１１のような波高分布スペクトルを得た。得られた各スペクトルから閾値を超えた信号の頻度を計数し、各中性子検出部の計数を比較した。
【００６０】
該比較結果と予め各方位角に線源がある場合の各検出部の計数の比較結果とを照合し線源の方向を求めた。なお、本実施例では、予め各方位角に線源がある場合の各検出部の計数を、前記式（１）の応答関数で表し、該応答関数から得られる各θにおける計算値（Ｎｉ）と、実際の計測で得られたｉ番目の検出部の計数の実測値（Ｎ’ｉ）を照合し、最小二乗法によって最も誤差が小さくなるθを求めた。
【００６１】
即ち、方位角１°毎に計算値と実測値の残差平方和｛Σ（Ｎｉ−Ｎ’ｉ）２｝を求め、該残差平方和が最小となる方位角を線源の方向とした。各方位角における残差平方和を図１２、図１３に示す。この結果より、方位角０°と特定することができた。
【実施例２】
【００６２】
実施例１と同様のセッティングを行い、中性子線源の位置を図１０のように減速体からみて９０°の方向となる線源位置（Ｃ）に配置し中性子を検出した。各中性子検出部から図１４のような波高分布スペクトルを得た。得られた各スペクトルから閾値を超えた信号の頻度を計数し、各中性子検出部の計数を比較した。
【００６３】
実施例１と同様にして、該比較結果と予め各方位角に線源がある場合の各検出部の計数の比較結果とを照合し線源の方向を求めたところ、図１５、図１６で示すように方位角８７°と特定することができた。
【実施例３】
【００６４】
実施例１と同様のセッティングを行い、中性子線源の位置を図１０のように減速体からみて４５°の方向となる線源位置（Ｂ）に配置し中性子を検出した。各中性子検出部から図１７のような波高分布スペクトルを得た。得られた各スペクトルから閾値を超えた信号の頻度を計数し、各中性子検出部の計数を比較した。
【００６５】
実施例１と同様にして、該比較結果と予め各方位角に線源がある場合の各検出部の計数の比較結果とを照合し線源の方向を求めたところ、図１８、図１９で示すように方位角４３°と特定することができた。
【図１】
2018013378.tif 000003
【図２】
2018013378.tif 000004
【図３】
2018013378.tif 000005
【図４】
2018013378.tif 000006
【図５】
2018013378.tif 000007
【図６】
2018013378.tif 000008
【図７】
2018013378.tif 000009
【図８】
2018013378.tif 000010
【図９】
2018013378.tif 000011
【図１０】
2018013378.tif 000012
【図１１】
2018013378.tif 000013
【図１２】
2018013378.tif 000014
【図１３】
2018013378.tif 000015
【図１４】
2018013378.tif 000016
【図１５】
2018013378.tif 000017
【図１６】
2018013378.tif 000018
【図１７】
2018013378.tif 000019
【図１８】
2018013378.tif 000020
【図１９】
2018013378.tif 000021