[論文レビュー] Real-time Tomography-based Bayesian Inference from TCV Bolometry Data

D. Hamm, C. Theiler+4 more|arXiv (Cornell University)|Mar 12, 2026

Magnetic confinement fusion research被引用数 0

ひとこと要約

この論文は、事前に計算された係数を用い、不確実性定量を提供する、TCV tokamak のボロメトリデータから放射電力を推定するリアルタイム断層 tomography に基づくベイズ法を提案します。

ABSTRACT

Radiated power information is crucial to diagnose and optimize the performance of fusion plasmas. Traditionally, at the TCV tokamak, radiated power analysis has only ever been possible following plasma discharge termination. However, recently, TCV bolometer data have become available in real-time. This offers the opportunity of integrating the radiated power information into the TCV plasma control system. In this work, we propose a novel real-time tomography-based Bayesian technique allowing estimation of the power radiated from user-defined regions of interest in the plasma. The real-time estimates are obtained as computationally cheap linear combinations of bolometer measurements, using pre-computed coefficients that are optimized for the specific discharge planned. This method is not, thus, trained on a set of synthetic or tomographically reconstructed emissivity profiles. We detail the derivation of the technique and show its equivalence to traditional tomographic estimates under suitable conditions. We then demonstrate that this technique enables accurate real-time estimation of the total, core, divertor and main chamber radiated power, by its application to a representative and heterogeneous set of TCV discharges. Finally, we discuss the robustness of the technique to faulty detectors, showing that simple precautions allow safe handling of many common issues. The computational routines implementing the described technique are provided as open-source code.

研究の動機と目的

• TCV ボロメトリデータから放射電力（総量、コア、ダイバタ、メインチャンバー）のリアルタイム推定を提供する。
• 磁気平衡情報を拡散事前分布として組み込み、空間的現実性を高める。
• ベイズ事後分散に基づく放射電力の不確実性定量を可能にする。
• 故障したディテクタへの頑健性を示し、事後放電後のトモグラフィー検証と比較して妥当性を検証する。
• 再現性と再利用性のためのオープンソース実装を公開する。

提案手法

• ボロメトリ inversion を拡散事前分布を用いるガウス後方としてモデル化し、磁気平衡情報をエンコードする。
• 放射電力の関心量が発光度の线形汎関数であることを導出し、閉形式の後方モーメントを得る。
• 計画放電平衡（FBT）からリアルタイム係数 beta_j を事前計算し、放射電力をボロメータ測定の線形結合として表現する。
• 後方平均 mu_post = A y を線形写像として計算し、放射電力推定を P_rad_tot = b^T x として表現することで、不確実性を伴うリアルタイム線形推定器を得る。
• 関心領域の任意の領域について、ピクセルを選択し対応する beta_j 係数を形成することでリアルタイム推定を提供する。
• チャネル故障などの実務的課題に対して、複数のチャネル選択戦略を評価し、頑健性を実証する。

実験結果

リサーチクエスチョン

• RQ1 オンラインのトモグラフィック逆解なしで、リアルタイム放射電力推定（総量、コア、ダイバタ、メインチャンバー）をボロメトリデータから取得できるか。
• RQ2 平衡情報を事前に取り込んだトモグラフィーに基づくベイズアプローチは、放射電力の信頼できる不確実性定量を提供するか。
• RQ3 故障した/挙動不良のボロメトリチャネルに対してリアルタイム推定はどれくらい頑健か。
• RQ4 異なる磁気構成にわたって、リアルタイム推定はポスト放電のトモグラフィックベンチマークとどれくらい近いか。
• RQ5 事前計算された係数を用いて、実際のリアルタイム平衡再構成ではなく、計画された放電平衡（FBT）に基づくアプローチを展開できるか。

主な発見

Strategy	E_tot	Δ_tot	r_tot	n_sigma_tot	E_core	Δ_core	r_core	n_sigma_core	E_div	Δ_div	r_div	n_sigma_div	E_main	Δ_main	r_main	n_sigma_main
S1	254.5 kW	37.9%	0.776	3.2	126.1 kW	5.5%	0.741	2.7	67.9 kW	10.1%	0.655	2.1	257.4 kW	46.1%	0.762	3.1
S2	10.5 kW	-0.1%	0.992	0.7	12.6 kW	6.0%	0.974	1.2	10.4 kW	7.1%	0.975	1.5	10.8 kW	-1.3%	0.994	0.7
S3	10.6 kW	-0.2%	0.992	0.7	12.6 kW	5.9%	0.974	1.2	10.6 kW	6.8%	0.970	1.5	10.9 kW	-1.3%	0.994	0.8
S4	10.5 kW	0.0%	0.992	0.7	12.5 kW	6.0%	0.974	1.2	10.6 kW	7.0%	0.970	1.5	10.8 kW	-1.0%	0.994	0.8

• 総量、コア、ダイバタ、メインチャンバーの放射電力のリアルタイム推定は、放電後のトモグラフィック推定を概ね追従する。
• ガウス後方を介した意味のある不確実性定量を提供し、放射電力の分散を計算可能にする。
• FBT から事前計算された係数は、複数の磁気構成に対してリアルタイム推定の正確性を満たす。
• 故障したチャネルに対しても、複数のチャネル選択戦略を試すことで頑健性と精度を維持する。
• オープンソースのコードと実装を公開し、コミュニティの利用と再現性を促進する。
• 50件の TCV 放電を横断して、リアルタイム推定はトモグラフィック参照と整合する信頼区間を提供した。