Feasibility demonstration of continuous signal-based neutron noise measurements by experiments and simulations

本論文は、シミュレーションおよび2基の研究用原子炉における実験を通じて、パルス波形デコンボリューションまたは検出器ペアを利用した連続信号中性子ノイズ解析が、従来のパルス計数法におけるデッドタイムおよびパイルアップの制限を効果的に克服し、高検出率においても偏りのない速度論的パラメータ推定を実現できることを実証するものである。

要約

大きな問題：「速すぎて数えられない」ジレンマ

屋根に降る雨粒を数えようとしている場面を想像してみてください。

• 従来の方法（パルス計数法）： バケツとカウンターを持ってそこに立っています。雨粒が当たるたびに、カウンターをクリックします。霧雨の程度であれば、これは非常にうまく機能します。
• 問題点： しかし、土砂降りになると、雨粒が激しく当たりすぎて、一つ一つの粒が重なってしまいます。どこで一つの粒が終わり、次の粒が始まるのか判別できなくなります。カウンターは「混乱」してしまいます（これはデッドタイムや**重なり（pile-up）**と呼ばれます）。その結果、カウント漏れが発生し、データは使い物にならなくなります。

原子炉において、科学者たちは原子炉がどのように動作しているかを理解するために、中性子（微小な粒子）を数える同様の手法を用いています。原子炉の出力が高くなったり、中性子が非常に高速で動いたりすると、中性の「雨」があまりに激しくなるため、従来のカウンターでは対応できなくなります。これでは、原子炉の安全性と効率を維持するために必要な、高速で重要な詳細を見逃してしまうことになります。

新しい解決策：粒を数えるのではなく、「ハミング」を聞く

この論文は、巧妙な回避策を提案しています。個々の雨粒を数えようとする代わりに、屋根を叩く雨の音を聞くことを想像してみてください。

• 連続信号： クリッカー（カウンター）を使う代わりに、屋根に伝わる連続的な「ハミング」や「振動」を記録するマイクロフォンを使用します。たとえ雨粒が重なり合っていたとしても、音波は雨がどれほど激しく、速く降っているかという情報を依然として保持しています。
• 目的： 科学者たちは、クリックを数えることで導き出していたものと同じ情報を、この「ハミング」（検出器からの連続的な電気信号）から導き出そうとしています。

検証方法：シミュレーションと実機実験

研究者たちは単に推測したわけではありません。彼らは2つの方法でこのアイデアをテストしました。

1. コンピュータ・シミュレーション（仮想ラボ）：
   コンピュータ上に仮想の原子炉を構築しました。中性の「嵐」をシミュレートし、従来の方法（クリックを数える）と新しい方法（ハミングを聞く）を比較しました。

   • 結果： 「嵐」が激しくなると、クリッカーは機能しなくなりました。しかし、「ハミング」による手法は、雨が非常に激しくなったときでも完璧に機能し続けました。また、クリッカーでは全く捉えられなかった「より速い」タイプの雨（高エネルギー中性子）も検出することができました。

2. 実機実験（現実世界）：
   このアイデアを、日本（KUCA）とハンガリー（BME TR）にある2つの実際の研究用原子炉に適用しました。

   • 特殊なマイクロフォン（核分裂計）を取り付け、連続的な電気信号を記録しました。
   • 原子炉を、非常に静かな状態からかなり騒がしい状態まで、異なる出力レベルで運転しました。
   • 結果： 静かな設定では、古いクリッカー法と新しいハミング法は一致しました。しかし、より騒がしい設定では、クリッカーは失敗（カウント漏れが発生）しましたが、ハミング法は正確な結果を出し続けました。

「ノイズ」の問題と「魔法のフィルター」

ただし、一つ問題がありました。マイクロフォンが風の音や電気的な静電気を拾ってしまうのと同様に、連続信号には、電子機器や信号自体の形状に起因する「ゴミ（不要な成分）」が含まれていました。これにより、信号が少し歪んでしまい、まるで電話の接続が悪い状態で話している声のように聞こえてしまいます。

これを解決するために、科学者たちは**デコンボリューション（逆畳み込み）**と呼ばれるデジタル技術を用いました。

• 例え話： 響きが悪い部屋で音楽を聴いている場面を想像してください。あなたは、本来の曲が「どうあるべきか」を正確に知っています。コンピュータを使って、数学的に「部屋の悪い音響特性」を取り除き、元の曲を復元することができるのです。
• 結果： この「魔法のフィルター」（具体的にはウィーナー・フィルター）を使用することで、彼らは信号をクリーンアップしました。これにより、ノイズを打ち消すための別の検出器を併用することなく、単一の検出器からクリアな結果を得ることが可能になりました。

主なポイント

• 従来の方法： 個々の中性子を数える方法は、動きが遅いときは有効ですが、速すぎたり激しかったりすると機能しません。
• 新しい方法： 連続的な電気的「ハミング」を分析する方法は、速く、かつ激しい状況でも機能します。信号が重なり合っても混乱することはありません。
• 解決策： もし信号が電子機器や検出器自体の形状によって歪んでしまった場合は、数学を用いてクリーンアップすることができます（デコンボリューション）。
• 結論： この手法は、原子炉の挙動を「聞く」ための、信頼できる「デッドタイムのない」方法です。これにより、信号が速すぎたり混雑しすぎたりするために、これまで測定不可能だった事象を測定できるようになります。

この論文が主張していないこと：
この論文は、この手法ががん治療や都市への電力供給、あるいは地震の予測に使えるとは主張していません。この研究はあくまで、中性子が速すぎてカウントできないという制限を克服することで、研究用原子炉の挙動をどのように「測定」し「診断」するかを改善することに特化しています。

技術概要

技術要約：連続信号に基づく中性子ノイズ測定の実現可能性実証

問題提起
Rossi-$\alpha$法やFeynman-$\alpha$法などの従来の中性子ノイズ手法は、プロンプト中性子の減衰定数（$\alpha$）などの動力学的パラメータを決定するために、離散的な検出器パルスを記録することに依存している。しかし、これらのパルス計数法は、高検出率において重大な制限を受ける。高フラックス環境では、パルスの重なり（パイルアップ）やデッドタイム効果が発生し、カウントの損失や、ノイズ解析に必要な統計的モーメントの歪みを引き起こす。デッドタイム補正によって平均計数率は回復できるものの、ノイズ手法に不可欠な高次モーメント（分散、共分散）を正確に補正することはできない。さらに、デッドタイム効果は測定可能な$\alpha$の大きさを制限するため、高速炉システム、極めて低次臨界状態、あるいは熱中性子系よりも$\alpha$の値が著しく大きい高次動力学モードの正確な特性評価を妨げる要因となる。

手法
本研究では、離散的なパルス計数の代替案として、核分裂室からの生の連続電圧信号を使用することの実現可能性を調査する。このアプローチは、検出器電流を決定論的な平均パルス形状 $f(t)$ とランダムな検出イベントとの畳み込みとして扱う、Pál、Pázsit、およびElter（2014）による核分裂室信号の確率モデルに基づいている。

本手法には、主に以下の2つの分析経路が含まれる：

1. 理論的導出： 確率モデルを拡張し、Rossi-$\alpha$（自己共分散）およびFeynman-$\alpha$（分散対平均）の式の連続信号版を導出する。これらの導出では、検出器パルスの有限な減衰定数（$\alpha_e$）を考慮しており、これが相関関数に歪み項を導入することを考慮に入れている。
2. 信号処理技術： パルス形状の歪みやデッドタイム効果を軽減するために、本研究では以下を用いる：
   • 検出器ペア： 有限なパルス幅に起因する自己相関アーティファクトを抑制するために、2つの検出器間の相互共分散を使用する。
   • デコンボリューション（逆畳み込み）： 逆フーリエ変換およびウィーナーフィルタリングを適用して、平均パルス形状を連続信号からデコンボリューションする。これは、理想的なディラックデルタパルスのシーケンスを復元することを目的としており、実質的に解析から検出器の時間定数を取り除くものである。
3. 検証： 実現可能性の評価は以下を通じて行われた：
   • シミュレーション： BMEトレーニングリアクターに基づく1速度点モンテカルロモデルを用い、様々なソース強度および$\alpha$値（最大$10^5$ s$^{-1}$）の下での検出器信号を生成した。
   • 実験： 京都大学臨界実験装置（KUCA）およびBMEトレーニングリアクター（BME TR）において測定を実施した。データ収集には、高速FPGAデジタイザ（125 MHzサンプリング）を使用し、連続信号を従来のパルスタイミングチェーンと共に記録した。

主な貢献と結果

• シミュレーション結果：

  • 高レート性能： シミュレーションにより、パイルアップのために従来のパルス計数が失敗する計数率（検出器あたり数$10^5$ events/s）においても、連続信号解析が正確性を維持できることが示された。
  • 高$\alpha$の測定： 連続的な手法は、パルスベースの手法が急速に劣化する高$\alpha$値（最大$10^5$ s$^{-1}$）を成功裏に推定した。
  • デコンボリューションの有効性： パルス形状のデコンボリューションは単一検出器による評価を大幅に改善し、検出器のパルス減衰定数がシステムの$\alpha$と同程度である場合でも、正確な$\alpha$推定を可能にした。ウィーナーフィルタリングは、単純な逆フーリエデコンボリューションよりも電子ノイズ（ノイズ対信号比2.5%まで）に対して頑健であることが示された。

• 実験結果：

  • KUCA測定： 亜臨界および低出力臨界構成において、連続信号解析は従来のパルス計数と一致する$\alpha$値を与えた。しかし、高出力レベル（CR-1, CR-2）では、非圧縮形式で記録された連続信号は、データ収集（DAQ）チェーンの非線形周波数伝達特性によるアーティファクトを示した。圧縮記録（信号の振れ幅のみを保存）を用いることで、低周波ノイズとDAQアーティファクトが軽減され、解析が可能となった。
  • デコンボリューションの適用： KUCAにおいて、高出力であるCR-2構成は、パイルアップのため標準的なパルス計数による解析が不可能であった。しかし、連続信号にウィーナーデコンボリューションを適用することで、デッドタイムの減少したパルスベースの信号が再構成され、低出力時の結果と一致するRossi-$\alpha$およびFeynman-$\alpha$の評価に成功した。
  • BME TR測定： 実験は本手法の実行可能性を裏付けたが、検出効率の低さ（検出器がコアの外側に配置されていること）による制限も浮き彫りにした。データには二次的な減衰定数（$\sim 1500$ s$^{-1}$）が観察されたが、これはOpenMCシミュレーションにより、コアの動力学ではなく、周囲の減速材や反射材からの中性子の散乱に起因するものと特定された。
  • アーティファクトの軽減： 本研究では、DAQチェーンの周波数伝達特性が人工的な短ラグ相関を導入する場合があることを特定した。圧縮記録および特定の平滑化技術が、これらのアーティファクトを抑制するのに有効であった。

意義と主張
本論文は、連続信号による中性子ノイズ解析が、高レートの原子炉診断において、従来のパルス計数に代わる堅牢で偏りのない選択肢であることを結論付けている。主な意義は、デッドタイムの制限によりパルスベースの手法では到達できない、より高い計数率およびより高い$\alpha$値への動作範囲の拡大にある。

著者らは以下のことを主張している：

1. デッドタイムからの独立性： 信号が適切に処理される限り、連続的な手法はパルス計数を悩ませるデッドタイムの制限を本質的に回避する。
2. 高次モードへのアクセス： より高い$\alpha$値を測定できる能力により、本技術は高速炉システム、深く低次臨界の状態、および高次動力学モードの特性評価に適している。
3. 実用的な実現可能性： 電子ノイズやDAQの周波数伝達特性が課題となるものの、これらは圧縮記録、検出器ペアリング、およびウィーナーデコンボリューションを通じて管理可能である。

著者らは現在の限界について謙虚な姿勢を保っており、実験結果の精度は現在、理論的な手法そのものよりも、データ収集ハードウェアや検出器の配置効率によって制約されていると述べている。技術の可能性を完全に実現するための次のステップとして、電子機器の改良と高効率検出器の使用が挙げられている。