JP2018159669A - 中性子検出器等の信号のみに基づいて核***性物質の組成、未臨界度、遅発中性子割合、中性子世代時間、即発中性子寿命を測定する方法。 - Google Patents

Abstract

【課題】例えば事故後の原子炉のように、核燃料体系の内部の構造や全体の形状の情報が得られない状況であっても、反応度を適切に測定することができる方法及びその方法を用いた装置を提供すること。【解決手段】中性子検出器によって検出される、核燃料内の時刻tにおける中性子数をn(t)、反応度等の変動後に安定した状態となった時の中性子数をn∞、未臨界度の関数をαy、中性子数n(t)の重み付き時間微分n・(t)をq(t)としたとき、n(t) = αq(t) + n∞の式が成り立つこと、及びαyqとn - n∞（中性子検出信号からノイズ成分を除いた計数率であるn とその安定時の値n∞の差）の比が時間によらず1になることを利用して、遅発中性子比率と遅発中性子崩壊定数を定める。【選択図】図１

Description

本発明は、例えば事故後の原子炉における燃料デブリなど、あらかじめ内部構造や形状が不明である未臨界状態の核燃料体系の動特性を、中性子検出器の信号のみに基づいて測定する方法に関する。

ウラン235、プルトニウム239など、中性子を吸収して核***を生じる核***性物質か、核***性物質と構造材等の核***性ではない物質の混在する物質（例えば原子炉では炉心全体）を以下では核燃料体系と呼ぶ。核燃料体系について、そのうちの核***性物質の組成（ウラン235が100%とか、ウラン235が90%で残りがプルトニウム239といった割合）を調べる主な方法は、少量のサンプルを採取して行う化学分析が知られている。

また、核燃料体系の未臨界度を測定する手法としては、従来、例えば稼働中の原子炉においては、未臨界度測定法の代表的なものとして、例えば（１）負のペリオド法、（２）制御棒落下法（特許文献１）、（３）補償法、（４）中性子源増倍法（特許文献２）、（５）逆動特性法（特許文献３）、（６）炉雑音解析法、（７）パルス中性子源法が知られている。

上述の７つの従来手法の内、最初の３つの手法は、初期状態として臨界状態にしてやる必要や（負のペリオド法、制御棒落下法）、予め校正された制御棒がすでに存在する必要が（補償法）ある。また、中性子源増倍法では未臨界度が既知である基準体系を必要とする。また、逆動特性法は、動特性パラメータの精度良い推定が必要であり、臨界近傍の浅い未臨界（中性子束分布の形状が臨界時の形状と等しいとみなせる程度の未臨界）でのみ用いることができる。また、炉雑音法は、工程異常が起きる場合では溶液条件や 即発中性子寿命が時間とともに変化するために、異常と判定されるアルファ値が変動する。さらに測定システムの時間追随性が問題となる。高次モードの影響の除去が困難である。また、最後に掲げたパルス中性子源法は、一定周期のパルス中性子源を必要とする上、入射中性子に励起される中性子束分布の高次モードの影響を補正するためには、あらかじめ数値解析が必要である。

実効遅発中性子割合の測定手法としては、²⁵²Cf中性子源法、Rossi-アルファ、Nelson 数法、改良Bennett 法、炉雑音法が知られている。これらの手法では、あらかじめ数値解析を行うために核燃料体系の内部構造に関する詳細な寸法を必要とする。²⁵²Cf中性子源法とNelson 数法では、核燃料体系が未臨界であるが臨界に近いことを前提としている。

さらにまた、中性子世代時間、即発中性子寿命については、通常、核***性物質やそれを含む（原子炉で言えば）炉心全体などの内部構造に関する詳細な寸法を用いて数値解析により求められる。

特開平８−３１３６６９号公報 特開平０９−１７８８８７号公報 特開２０１６−１４２６２３号公報

核燃料体系内部の核***性物質の組成を調べるための化学分析では、対象のサンプルを採取する必要がある。また、従来の未臨界度測定手法のうち、上述のパルス中性子源法を除くいずれの方法も、動特性パラメータの解析値の精度に依存する。また、パルス中性子源法は補正用の数値解析の精度に依存する。またこれらの数値解析には、核燃料体系の内部の構造や全体の形状の情報もしくは制御棒等の反応度を必要とする。

総じて、従来の手法では、数値解析のために核燃料体系の内部構造に関する詳細な寸法を必要とする。それ故に、未臨界度、実効遅発中性子割合、中性子世代時間、即発中性子寿命の評価結果はそれらの寸法の精度や解析に用いる計算コードの計算精度に依存する。

本発明の目的は、広範囲の未臨界（中性子束分布の形状を限定しない未臨界）において、通常の原子炉にあっては、炉内の寸法の精度や数値計算コードの精度の影響を受けないで、核***性物質の組成、未臨界度、遅発中性子割合、中性子世代時間、即発中性子寿命を測定する方法を提供することにある。

例えば事故後の原子炉のように、核燃料体系の内部の構造や全体の形状の情報が得られない状況であっても、核***性物質の組成、未臨界度、遅発中性子割合、中性子世代時間、即発中性子寿命を適切に測定することができる方法を提供することにある。

本発明は、未臨界状態の核燃料体系について、内部構造や全体の形状についてあらかじめ分かっていないために予備解析等によって正確なパラメータを事前に求めることができない状況又は、予備解析由来の誤差を排除したい状況、予め構成された制御棒等が使用できない状況、事前に臨界にできない状況、初めて到達した未臨界度において測定したい状況で、正又は負の反応度の添加により臨界に近づくもしくは遠ざかる場合かつ/又は外部中性子源の出力変動があった場合に、中性子検出器等からの核***率に比例する信号のみに基づいて、核***性物質の組成、未臨界度、遅発中性子割合、中性子世代時間、即発中性子寿命を測定する方法である。

具体的には、中性子検出器によって検出される、核燃料体系内の時刻tにおける中性子数をn(t)、未臨界度をρ(t)、中性子源強度をS(t)とすると、未臨界度と中性子源強度のどちらかもしくは両方が変動した後でその変動がほとんど停止した時（ρ(t) ≒ S(t) ≒ 0）から最終的に安定した状態になってn(t)が一定の値n_∞となるまでの間に、以下の式が成り立つことを用いる。

ここでα_yは、実効遅発中性子割合をβとすると、ドル単位の未臨界度ρ_$=ρ/βにより
α_y = 1/ρ_$−1と表され、q(t)はn(t)の時間的変化率n^・(t)と以下の式で表される時間の関数μ(t)の比（n(t)の重み付き時間微分）である。

Ci^・(t)は、n(t)が中性子数である場合には、第i群の遅発中性子先行核の量Ci(t)の時刻tにおける時間的変化率を表す。通常は6群が用いられるが、以下の測定では群数に制限はない。λ_iは第i群の遅発中性子先行核の崩壊定数である。

式１において、α_yはドル単位の未臨界度ρ_$のみ、n_∞は未臨界度ρ及び中性子源強度Sの関数で、βは未臨界度ρの関数であるから、ρ及びSの時間的変化が十分に小さければ、α_yとn_∞は定数とみなせる。このとき、式１はα_yを傾き、n_∞を切片とする直線の方程式とみなすことができる。以下の測定では、このことを利用する。

中性子源強度としては、中性子源から放出される中性子のうち、対象となる核燃料体系に到達するもののみを考慮する（実効的な中性子源強度）。

n(t)については、中性子検出器信号（計数率）から校正定数が求められている場合にはそれを用いれば容易に得ることができる。あらかじめ校正できない場合には、中性子検出器信号を用いる。不感時間の補正など、必要かつ実施可能な補正や変換を行う。

具体的には、本発明は、式１が成り立つことを用いた、以下の方法にある。

（１）α_yqとn-n_∞の比の時系列データと1との差が最も小さくなるように、遅発中性子比率と遅発中性子減衰定数を定める方法。

（２）時系列データ（q, n）が直線n(t)=α_yq(t)+n_∞上にあることを利用して、未臨界度ρ_$とn(t)の安定時の値n_∞を求める方法。

（３）上記項目（１）または（２）において、n = exp(f (t))を用いて（n^・, n）の組を求める方法。

（４）上記項目（２）において、a = 1/ρ₀-1/ρ_k、b = 1-r、r = n_0∞/n_k∞とすると、ρ_k$→ρ_0$ のとき、a/bが未臨界度ρ₀におけるβに収束することを用いて実効遅発中性子割合βを求める方法。

（５）上記項目（２）において、Λ／β=l(-ρ)+l／β が成り立つことを利用してl[s]、Λ[s]を求める方法。

未臨界状態の核燃料について、内部構造や全体の形状についてあらかじめ分かっていないために予備解析等によって正確なパラメータを事前に求めることができない状況又は、予備解析由来の誤差を排除したい状況、予め構成された制御棒等が使用できない状況、事前に臨界にできない状況、初めて到達した未臨界度において測定したい状況で、正又は負の反応度の添加により臨界に近づくもしくは遠ざかる場合かつ/又は外部中性子源の出力変動があった場合に、中性子検出器等からの信号のみに基づいて、遅発中性子比率と遅発中性子減衰定数、ドル単位の反応度を求めることが可能となった。

本発明に係る未臨界状態の核燃料体系の反応度測定装置の概略説明図。 検出器信号（青色マーカー）に直線をフィッティングして傾きα_yと切片n_∞の値を求めた図（フィッティングの結果はα_y = −1.094、n_∞ = 6.902×10^-5）。 ウラン235が100%であるような場合に、時々刻々のα_yq/(n−n_∞) の値を、遅発中性子比率と遅発中性子先行核崩壊定数の値としてウラン235の値を用いた場合（青）とプルトニウム239の値を用いた場合（赤）の比較を示した図。 反応度が添加され、未臨界度が浅くなったあとの中性子検出器信号の変化の様子を表した図。 図４の時系列データを加工して(n_∞−n)の時系列データを作成した図。 図５の時系列データを加工して(n_∞−n)の自然対数、ln(n_∞−n)、 の時系列データを作成しこれに直線y=at+bをフィッティングして、係数a、bを求めた図。 図６で求めた係数a, bの値を用いてn = n_∞ + e^at+bの時系列データを作成し、元の図４のデータと比較した図。 異なった反応度ρ_0$, ρ_1$, ρ_2$，・・・に対して(1-r)をx軸、(1/ρ_0$-r/ρ_k$)をy軸(r = n_0∞/n_k∞)としてプロットし、原点と最も原点に近い点を用いて傾きβを求めた図。

本発明に係る未臨界状態の核燃料体系の反応度測定装置の基本構成を図１に示す。反応度測定装置１０は、中性子検出器等から成る検出部１、数値計算を行うためのデータを入力するための入力部２、数値計算を行う演算部３、演算部３に接続された、数値計算に必要なデータを保存する記憶部４、測定結果を外部に伝える出力部５からなる。検出部１の信号は核***で放出される出力の値に相当するように予め校正しておく。出力部５は、数値データを外部機器に伝送したり、結果を分かり易く画面に表示したりすることができる構成になっている。

本発明に係る反応度測定装置１０の内、検出部１を除く部分は、基本的にマイクロコンピューターで構成できる。したがって、臨界安全管理されていない容器に本発明に係る装置を据え付けることで、核燃料物質の誤移送による臨界近接を検知することもできる。

次に、上述の反応度測定装置１０の演算部３で行う、未臨界状態の核燃料体系の反応度測定方法の演算手順について説明する。

本発明で使用する方程式は、次の考え方に基づいて一点炉動特性方程式から導出したものである。核燃料体系に対して未臨界度ρの変化や中性子源強度Sの変化といった外的擾乱が与えられたとき、それに対する応答としての中性子数nの変動は未臨界度ρと密接に結びついているはずであり、特にρとSの時間的変化率（dρ/dt及びdS/dt）が非常に小さい場合には、中性子数nの変動（nの値及びその時間変化率dn/dt）は未臨界度ρの影響を受けつつ、遅発中性子先行核の崩壊（遅発中性子先行核濃度Ci、遅発中性子先行核崩壊定数λi、遅発中性子比率βiで特徴づけられる）、により供給される中性子に支配されると考えられるので、一点炉動特性方程式にρとSの時間的変化率（dρ/dt及びdS/dt）が非常に小さい条件を適用して当該方程式を導出した。なお、本明細書では、便宜上、例えば、nの時間微分をn^・のように示す。また、Ciについても同様にCi^・のように示す。

＜（１）反応度測定方法＞
検出部からの信号をnとし、その時間微分をn^・で表す。nより次式でμを求める。

ここで、βは実効遅発中性子割合、β_iは第i群の遅発中性子比率、Λは中性子世代時間、λ_iは第i群の遅発中性子先行核崩壊定数、C_iは第i群の遅発中性子先行核量である。 Ci^・はC_iの時間微分を表す。
n^・とμより次式でqを求める。

データの組（q, n）を用いて自己回帰分析によりα_yとn_∞ を求める。α_yよりドル単位の反応度ρ_$を次式で求める。

＜（２）遅発中性子比率測定方法＞
図３を参照する。（１）において、α_yq とn - n_∞の比を時系列データとして求める。この比が計算対象となる時間の範囲においてできる限り1に近くなるように、（β_i、λ_i）の値を調整する。多変量回帰分析の手法を用いても良い。あらかじめ想定される核***性物質がある場合は、その値の組（β_i、λ_i）を物質の種類ごとに用意し、（例えばU-235, Pu-239など）それらの混合割合に応じて（β_i、λ_i）の内挿値をもとめ、得られた値の組（β_i、λ_i）を適用してα_yq とn - n_∞との比を計算する。この比が最も１に近くなったときの混合割合が、核***性物質の混合割合となる。

＜（３）時刻ｔにおける（n^・, n）の値の組を求める方法＞
n（図４）もしくはn - n_∞またはn_∞ - n（図５）の時系列データの対数をとって得られたデータをn₁とする（図６）。得られたn₁データを時間の関数で近似し（図７）、その関数をもとに、時刻tにおける（n^・, n）値の組を求める。

＜（４）あらたな未臨界度においてα_y とn_∞の値を推定する方法＞
（３）において、（１）を一度以上実施した後であれば、α_y とn_∞の最初の推定値を以下のように得ることができる。先に実施して得られた未臨界度ρ₀ ( = ρ_$0×β )と安定時中性子数n_0∞を用いて得られる次の値p₀を用いると、

未臨界度がそれほど大きく変わらなければ、前述の式１は常に点（p₀, −p₀）を通るので、あらたな未臨界度ρ₁におけるα_1y とn_1∞の最初の推定値を次のようにして得ることができる。

ここで、t₀は未臨界度と中性子源強度の変動が停止した直後（未臨界度はρ₁）以降の任意の時刻である。

＜（４）実行遅発中性子割合βを求める方法＞
上記（１）により求めた、互いに異なった未臨界度の二組以上のデータ(ρ_k$, n_k∞)（うち実効遅発中性子割合を求めたい反応度におけるデータの組を（ρ_0$, n_0∞）とする）を用いて、a= 1/ρ_0$ -r/ρ_k$ 及びb = 1 - r_kを求める（ここで r_{k =} n_0∞/ n_k∞ とする）。データの組(a, b)が原点に近づくにつれて比a/bが実効遅発中性子割合βに収束することを用いて、βを求める。

以下、核燃料物質に設置された中性子検出器からの信号（中性子計数率n[数/s]）を用いたデータ処理手順をさらに詳細に説明する。

中性子検出器の信号もしくはそれを校正、加工した信号は、一定もしくは任意の時間間隔（Δt）で得られるものとし、その得られたデータ列n₁、n₂、n₃、・・・、n_j、を以下では時系列データと呼ぶ。以下ではこれらをまとめてn(t)と書く。

ここで取り扱う物理量について説明する。
ウランなどの核***により生じる核***片のうち、一定の時定数で崩壊して中性子を放出する核種を遅発中性子先行核という。有限の大きさを持つ現実の体系で遅発中性子の中性子の全数（即発中性子と遅発中性子の和）に対する割合が実効遅発中性子割合（βであらわす）である。遅発中性子先行核は多数存在するが、通常、崩壊定数の近いものをまとめて一つの群として扱い、その濃度（単位体積当たりの数）を遅発中性子先行核濃度と呼び、Ci(r,t)であらわす。iは群の番号で、通常６個の群が用いられる。Ci(r,t)を体系全体で空間積分したものをCi(t)であらわす。これは第i群の遅発中性子先行核の時刻tにおける総量を表す。遅発中性子のうち、Ci(t)から生じたものの比率（第i群の遅発中性子比率）をbiであらわす。biは式９を満たす。第i群の遅発中性子先行核の崩壊定数をλiとする。６個の(bi、λi)の組は、核***性物質（ウラン２３５、プルトニウム２３９など）ごとに異なった値が提案されている（文献James J. Duderstadt and Louis J. Hamilton, “Nuclear Reactor Analysis”, John Wiley & Sons, ISBN 0-471-22363-8を参照）。

評価する物理量は順に、中性子検出器信号からノイズ成分を除いた係数率n [数/s]、この中性子計数率の時間変化率n^・（=dn/dt）[数/s²]、遅発中性子先行核相当量Ci(t)、Ci(t)の時間変化率Ci(t)、変数μ(t)[1/s]、変数q(t)[1/s]、未臨界度ρの関数α(ρ)[-]および安定時のnの値n_∞(ρ)[数/s]、ドル単位の未臨界度ρ_$(=ρ／β)[-]、第i群の遅発中性子先行核崩壊定数λi[1/s]、第i群の遅発中性子比率bi[-]、実効遅発中性子割合β[-]、即発中性子寿命l[s]、中性子世代時間Λ[s]である。nと n^・を含む時間の関数となっている変数は時系列データである。手順は、未臨界度の変化もしくは中性子源強度の変化により中性子検出器信号が上昇する場合、手順A、下降する場合、手順Bを用いる。

＜手順1-A（ノイズ成分を除いたnおよび n^・の評価）＞
検出器の信号が安定した後ならば、n_∞の測定値を、変動中であれば適当な予測値を用いてx = ln(n_∞-n )を計算し、この時系列データx(t)を適当な時間の関数f(t)で近似する。たとえばx(t)のデータに対してy=at+bをフィッティングして定数a,bを求めることで、f(t)=at+bでx(t)を近似する。このf(t)を用いてnと n^・はそれぞれ下式１０、１１とあらわされる。

＜手順1-B（ノイズ成分を除いたnおよびn^・の評価）＞
検出器の信号が安定した後ならば、n_∞の測定値を、変動中であれば適当な予測値を用いてx = ln(n-n_∞)を計算し、この時系列データx(t)を適当な時間の関数f(t)で近似する。たとえばx(t)のデータに対してy=at+bをフィッティングして定数a,bを求めることで、f(t)=at+bでx(t)を近似する。このf(t)を用いてnと n^・はそれぞれ下式１２、１３とあらわされる。

上述の手順1-A、1-Bにおいて、n_∞に予測値を使った場合は、得られたnの時系列データと中性子検出器信号の時系列データを比較して、すべてのデータの二乗和を評価し、これが最も小さくなるようにn_∞の値を調節して、手順1-A,1-Bを再度実施する。この手順を繰り返して、差が最も小さくなるようにする。

手順1において、手順5を一度以上実施した後であれば、n_∞の最初の推定値を以下のように得ることができる。先に実施して得られた未臨界度ρ₀ ( = ρ_$0×β )と安定時中性子数n_0∞を用いて得られる次の値p₀を用いると、

未臨界度がそれほど大きく変わらなければ、式１は常に点（p₀, -p₀）を通るので、あらたな未臨界度ρ₁におけるn_1∞の最初の推定値を

ここで、t₀は未臨界度と中性子源強度の変動が停止した直後（未臨界度はρ₁）以降の任意の時刻である。

これまでのやり方では、一点炉動特性方程式に基づくため、その解の形から、例えば手順1-Aの場合、n = n_∞ - ΣAie^f(ti) のような関数形を考える必要があった。しかし、n = α_yq - n_∞ が成立つことにもとづく本手法では、n = α_yq - n_∞ より解の形はn =n_∞ - e^f(t) となることが示されるため、より簡単にノイズ成分のないnと n^・の値を求めることができる。

＜手順2（Ci(t)、Ci^・(t)の評価）＞
bi、λi、Λ、について、想定される核燃料物質の値もしくは、もっとも近いと想定できる仮の値を用いて、以下の式により、Ci^・(t)を求める。ここでβi=β×biである。

Δt後のCi(t+Δt)を以下により求める。

このような数値積分はルンゲクッタ法や陰解法など多数存在するので、どれを用いてもよい。Ciの初期値は、中性子検出器信号が安定している状態のときに

として求めておく。

＜手順3（μの評価）＞
時刻tにおけるμ(t)の値は以下の計算により評価する。

＜手順4（qの評価）＞
時刻tにおけるq(t)の値は以下の計算により評価する。

＜手順5（α_y、n_∞、ρ_$評価）＞
上述の手順4で求めたq(t)の時系列データと手順（１）で求めたn(t)の時系列データを用いる。横軸にq(t)、縦軸にnをとって、時刻ごとに点（q(t),n(t)）をプロットする。このとき、の点の軌跡に直線の方程式をフィッティングして得られる傾きの値がα_y 、切片の値がn_∞である。次式でρ_$を計算する。

＜手順6（λi[1/s]、bi[-]の評価）＞
上述の手順5で得られたα_y 、n_∞の値を用いて以下の時系列データy(t)を計算する。

縦軸をy、横軸をtとして点(t,y(t))をプロットする。このとき、y(t)が直線y=1に最も近くなるように、λi[1/s]、bi[-]を調整する。このとき、biの全群の総和は1となるようにする。

より簡単には、核***性物質の種類に応じて、ウラン235、プルトニウム239など、それぞれに（λi[1/s]、bi[-]）の組（i＝1,・・・,6）が文献で与えられている。以下では、ウラン２３５のそれをλ5i[1/s]、b5i[-]、プルトニウム２３９のそれをλ9i[1/s]、b9i[-]とする。たとえば対象としている核燃料物質中にウラン２３５とプルトニウム２３９が存在することが想定されているが、その割合が不明のときにλi[1/s]、bi[-]を調整する場合、核***性物質全体のうちのウラン２３５の個数割合をRとすると、λi[1/s]、bi[-]は以下の式で表される。

この計算をi=1〜6のすべてで同時に行う。λi[1/s]、bi[-]を個別に変更しない。Rを変えて上述の手順2から手順6を繰り返し、(y(t)-1)²の総和が最も小さくなったときのRが、対象とする核燃料物質におけるウラン２３５の個数割合である。

＜手順7（βの評価）＞
上述の手順5で求めたn_∞、ρ_$の値が一組あるとする。それをn_∞0、ρ_$0であらわす。未臨界度が変化して新たな未臨界度ρ_$になったとき、上述の手順1から手順6を適用してn_∞を求める。二つの安定出力の比をr(=n_∞／ n_∞0)とすると、βは次式で求められる。

別の例としては、複数の値の組（n_∞i、ρ_$i）を用いて、点（n_∞0、ρ_$0）に対してyi = 1/ρ_$0-r/ρ_$i とxi = 1 - ri (ri=n_∞0／ n_∞i)を求め、これらの点（xi, yi）をプロットする。これらの点に原点を通る直線をフィッティングしたとき、この直線の傾きからβを得る。ここで得られたβの値を手順２で用いることにより、仮の値であったβの値を正確な値とすることができるため、手順２を安定的に実施することができる。

＜手順8（l[s]、Λ[s]の評価）＞
未臨界度がρ_$で一定の状態が続き、中性子検出器信号が安定した値ｎ_∞となっている状態を初期状態とする。外部中性子源強度はSとする。外部中性子源（パルサトロンなど）の電源をOFFにしたり遮蔽したりして、瞬時にS=0の状態にする。上述の手順1-Bで得られたｎとn^・からn^・/nの最大値ωを求める。Λ/β=ρ_$/ωによりΛ/βを求める。未臨界度が異なる２組以上のデータの組（−ρ_$k、Λk/β）をプロットする。直線y=ax+bをこれらの点にフィッティングして得られるa, bから、l = aもしくは l＝b×βによりlを評価する。Λ=βρ_$/ωもしくは、Λ= l/keff （keff=1/(1−ρ)）によりΛを求める。ここで得たΛの値を上述の手順2で用いることで、仮の値であったΛの値を正確な値とすることができるため、上述の手順2を安定的に実施することができる。

未臨界状態の核燃料について、内部構造や全体の形状についてあらかじめ分かっていないために予備解析等によって正確なパラメータを事前に求めることができない状況又は、予備解析由来の誤差を排除したい状況、予め構成された制御棒等が使用できない状況で、正又は負の反応度の添加により臨界に近づくもしくは遠ざかる場合かつ/又は外部中性子源の出力変動があった場合に、上述の手順1により、中性子検出器等からの核***率に比例する信号からノイズのない信号を抽出することが可能となった。上述の手順2から4により中性子数の重み付き時間微分を計算できる。手順5により、中性子検出器等からの核***率に比例する信号のみに基づいてドル単位の未臨界度を求めることが可能となった。手順6により遅発中性子比率と遅発中性子先行核崩壊定数を求めることが可能となり、さらにここで求めた遅発中性子比率と遅発中性子先行核崩壊定数を用いて上述の手順1から5の計算を行うことにより、初期値として与えた遅発中性子比率等の値に依存しない、中性子検出器等からの核***率に比例する信号のみに基づいた、ドル単位の未臨界度を求めることが可能となった。さらに、上述の手順7により実効遅発中性子割合を求めることが可能となり、これをドル単位の反応度に乗ずることで、ドル単位ではない未臨界度を求めることが可能となった。手順8により中性子世代時間及び即発中性子寿命を測定することが可能となった。

本発明は、再処理施設においても、使用済み燃料を受け入れる際に中性子を照射することで核***性物質の比率（U-235, Pu-239, 241 など）を測定できるので、本発明は燃焼度測定に利用できる。また、臨界安全管理されていない容器に本発明に係る装置を据え付けることで、核燃料物質の誤移送による臨界近接を検知することができる。その際、流入する核燃料物質の種類が変化しても、自動的に補正して中性子実効増倍率を計算し続けることができる。中身が不明な物質に対して中性子を照射することで、核***性物質の同定ができる。原子炉の臨界近接の自動化が可能となる。また、本発明は、ADS（加速器駆動未臨界炉）での運転制御に必要な未臨界度測定方法の有力な候補となる。

１・・・検出部
２・・・入力部
３・・・演算部
４・・・記憶装置
５・・・出力部

Claims (6)

1. 中性子検出器の検出信号もしくはこれを校正することによって得られる、核燃料内の時刻tにおける中性子数をn(t)、反応度等の変動後に安定した状態となった時の検出信号もしくはこれを校正することによって得られる中性子数をn_∞、未臨界度の関数をα_y、検出信号もしくはこれを校正することによって得られる中性子数n(t)の重み付き時間微分n^・(t)をq(t)としたとき、
   n(t) = αq(t) + n_∞ の式が成り立つことを利用して、α_yqとn - n_∞（中性子検出信号からノイズ成分を除いた計数率であるn とその安定時の値n_∞ の差）の比が測定期間中のいずれの時刻においても最も１に近くなるように、遅発中性子比率と遅発中性子崩壊定数を定める方法。
2. 請求項１において、中性子検出器の検出信号を基に求められる、時系列データ（q, n）が、直線n = α_yq + n_∞ 上にあることを利用して、未臨界度ρ_$とn（ｔ）の安定値n_∞ を求める方法。
3. 請求項１又は２において、n = exp(f (t))を用いて（n^・, n）の組を求める方法。
   ここで、f(t) はn(t)が検出器信号を再現するように定めたtの関数である。
4. 請求項１又は３において、前回の測定値（未臨界度ρ_0$、中性子検出器の安定値n_0∞ ）を用いて新たな未臨界度におけるn_∞ を推定する方法。
5. 請求項２において、a = 1/ρ₀ - ｒ/ρ_k、b = 1 - ｒ, ｒ = n_0∞/n_k∞とすると、ρ_k→ρ₀のとき、a/bが未臨界度ρ₀におけるβに収束することを用いて実効遅発中性子割合βを求める方法。
6. 請求項2において、Λ/β= l（-ρ）+ l/β が成り立つことを利用して、即発中性子寿命l [s]、中性子世代時間Λ[s]を求める方法。

Images