JPS6232392A - 原子炉のオンライン出力分布監視装置 - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 ［発明の技術分野］ 本発明は原子炉のオンライン出力分布監視装置に関する
。

［発明の技術的背景とその問題点］ 一般に、原子炉において炉心がその健全性を保ちつつ、
必要な性能を発揮しているか否かを監視することは、原
子炉の安全およびその効率的運用にとって不可欠のこと
である。このため、通常原子炉にはオンライン出力分布
監視装置が設置されており、この監視装置によって随時
原子炉の健全性を監視するための主要なパラメータの１
つである炉内出力分布を監視することができるように構
成されている。

ところで、従来のオンライン出力分布監視装ａにおける
炉内出力分布の算出方法は以下に説明するような手順で
なされていた。すなわち、まず炉心内中性子束測定器の
測定値と、これに隣接する燃料集合体の出力との間の関
係を予めオフライン計算により評価し、その結果をフィ
ッティング式として出力分布監視装置に内蔵しておく。

次に、炉心内中性子束測定器の測定値をオンラインで入
力し、前記フィッティング式によって燃料集合体の出力
を算出してオンラインでの炉内出力分布を推定していた
。しかしながら上述した従来の炉内出力分布推定方法は
、炉心内中性子束測定器が中性子束を測定する際の測定
誤差が無視しくワる程小さい場合にはよい精度で炉内出
力分布を推定することができる。しかし実際にはこの中
性子束の測定誤差は無視することができない。それは測
定器を案内する導管ないしは測定器自体の偏位が原子炉
の運用によって不可避的に発生ずることであり、これに
よって、従来の炉内出力分布推定方法によるときには、
不可避な誤差または無視できない誤差が発生するのでこ
れらの誤謬を運転マージンとして考慮する必要が生じ、
このことは原子炉の効率的な運用の妨げとなっていた。

また、近年上述した従来の炉内出力分ｔ５推定方法のよ
うに、炉心内中性子束測定器の測定値とこれに隣接する
燃料集合体の出力との間の関係を予め内蔵しておいたフ
ィッティング式により評価するのではなく、オンライン
出力分布監視装置に炉心の物理モデルを内蔵した出力分
布演算装置を備え、この計算装置により算出された炉内
出力分有ｉが炉心内中性子束測定器の測定値に一致Ｊ゛
るように物理モデルを適合することによって炉内出力分
布を推定する方法が提案されている。しかしながら、こ
の炉内出力分布推定方法によっ−Ｃも測定値を全面的に
真と仮定して炉内出力分布を推定するときには、上述し
たように正確な炉内出力分布を推定することができない
という不具合があった。

［発明の目的１ 本発明は、上記事情に鑑みてなされたちので、その目的
は、炉心径方向の変異に対してもその測定精度がほとん
ど変化しないγ線束測定器を使用して、最良の炉心内の
出力分布を出力することのできる原子炉のオンライン出
力分布監視装置を促供づることにある。

「発明の概要」 本発明は、上記目的を達成するために、炉心現状データ
測定器と、この炉心現状データ測定器からの測定値を入
力し予め内蔵された物理モデルに基づき炉心内の出力分
布を算出する出力分布演算装置と、炉心内γ線束測定器
と、この炉心内γ線束測定器の測定値を入力してその炉
心軸方向成分のみを抽出するとともに前記出力分布演算
装置で算出された出力分布を入力して得られた炉心軸方
向成分を前記炉心内γ線束測定器の測定値の炉心軸方向
成分により補正することにより最良出力分布を推定する
最良出力分布推定装置とを備えた原子炉のオンライン出
力分布監視装置に関するものである。

［発明の実施例］ 本発明の一実施例を図面について説明する。

第１図は本発明の一実施例のブロック構成図であり、同
図に示すように、原子炉１内には炉心２が収納されてお
り、この炉心２内には固定型中性子測定器３おび可変型
γ線束測定器４が設けられている。また、原子炉１内に
は炉心現状データすなわち炉心冷却水流量、炉心圧力、
炉心出入口温度、４１１制棒位置などを測定する測定器
５（以下代表して符＠５で示す）が設けられており、こ
れら各測定器３〜５で測定された出力信号は原子炉のオ
ン、ライン出力分布監視装置６に入力される。この監視
装置６はデータサンシラ７、出力分布計ロ装置８、最良
出力分布推定装置９および入出力装置１０とから構成さ
れている。

ところで、炉心内γ線束測定器４は通常第２図に示され
るように、燃料集合体１１を構成するコーナーギャップ
１２に導管１３（以下ストリングという）を設置し、そ
のストリング１３の内部に可動型γ線束測定器４が挿入
され、炉心軸方向に沿ってスキャンすることにより、炉
心軸方向に沿う連続したγ線束分布を測定するものであ
る。このス１−リング１３内部には約４個の固定型中性
子束測定器３も設置されている。また、制御棒１４はス
トリング１３の設置されるギャップの燃料集合体を挾む
反対側に挿入される。

しかして、データサンプラー７は原子炉１の炉心２内に
設置された固定型中性子束測定器３（以下ＬＰＰＭとい
う）および可動型γ線束測定器４（以下ＴＩＰという）
の出力信号を入力し、これらの出力信号を最良出力分布
推定装置９に伝送する。ここで、ＬＰＰＭの出力信号は
常時発生しているが、データサンプラ−７は出力分布ｉ
ｔ　Ｒ装置ε３からの伝送要求があった場合にのみ前記
伝送を実行する。また、ＴＩＰによる炉心内γ線束測定
は通常は１か月に一回程度の傾度で行われ、この場合に
はデータサンプラー７は常にＴＩＰの出力信号を最良出
力分布推定装置９に伝送する。

データサンプラー７は、炉心現状データずなわら炉心冷
却水流量、炉心圧力、炉心出入口温度および制御棒位置
等を測定する測定器５からの出力信号を入力し、これら
の信号を出力分布計算装置８に伝送Ｊ゛る。ところで、
測定器５からの出力信号は常時発生しているが、これら
の出力信号の伝送は前記Ｌ　Ｐ　Ｐ　Ｍ出力信号ないし
ＴＩＰ出力信号の最良出力分布推定装置９への伝送と同
時に実行される。出力分布計ｐ装置８は、予め定められ
た周期に従って自動的に起動するかまたは原子炉運転員
の入出力装置１０を介しての要求により起動リ−るかま
たは前記ＴＩＰによる炉心内中性子束測定の終了を示す
信号をトリが−として自動的に起動づるか等してデータ
サンプラ−７に対して前記ＴＩ　Ｐの出力信号ないしＬ
ＰＰＭの出力信号の最良出力分布推定装置９への伝送要
求を行ない、データリンブラー７から炉心現状データを
入力し、その内蔵する物理モデルによって炉内出力分布
を算出する。このような演算に用いられる物理モデルと
しては、装置全体の即応性を高めるため多少計輝精度が
悪くても計算時間が短くて済む、例えば、いわゆる１群
・粗格子点拡散モデルを使用するのが望ましい。このよ
うな１群・粗格子点拡散モデルを用いて炉内の出力分布
計算を行う場合には、第２図に示されたバンドル１１を
空間的に多数の接点（ノード）の集合で表し、各ノード
における中性子束が次の方程式（Ｉ）を満足するものと
する。

Ｌ・ΦトＢ　・Φ＝Ｏ・・・（Ｉ） ここで、 Ｌニラプラス演算子を有限階差近似した、ノードの総数
側に等しい行・列数の正方行列Φ：各ノードにおける中
性子束の全体を要素とする、ノードの総数側に等しい次
元のベクトル Ｂ２：各ノードの物質バックリングを対角要素とする、
ノードの総数側に等しい行・列数の正方対角行列 ざらに、原子炉炉心を三次元直角座標系で記述し、炉心
軸方向を７軸、炉心径方向をＸ軸、ｙ軸で夫々表し、１
方向でに番目、×方向で１番目、ｙ方向で、３番目に存
在するノードを添字に、Ｉ。

Ｊで示すことにし、Φの要素たるノードの中性子束をΦ
（Ｋ、ｌ、Ｊ）と書くこととする。すると、炉内出力分
布は各ノードにおける出力密度Ｐ　（Ｋ。

１、Ｊ）は次の方程式（［）の全体から成る集合で与え
られる。

Ｐ（Ｋ、１．Ｊ） ＝Σす（Ｋ、Ｉ、Ｊ）　　・　Φ　（Ｋ、１．Ｊ）　　
・・・　（Ｉｔ）ここで、 Ｐ（Ｋ、１．Ｊ）：ノード（Ｋ、１．Ｊ）における出力
密度、ただし単位 は単位時間・単位体積当り の核分裂数 殉（Ｋ、Ｉ、Ｊ）：ノード（Ｋ、ｒ、Ｊ）における巨視
的核分裂断面積。

上記方程式（Ｉ）および方程式（ｎ）に基づいて、出力
分布計算装置８は周知の方法で炉内出力分布を算出する
。この計算法は、例えば「築城諒外著“］ンバージエン
ス　アンド　アクセレーション　オブ　ボイド　イｔテ
レ−シコン　インボイリング　ウォータ　リアクター　
コア　力リキュレーション”ニュクリア　サイエンス　
エンジニア　６４，７２４−７３２（１９７７）　　［
Ｍ、　Ｔｓｕｉｋｉ　ｅｔａｌ：“Ｃ０ｎＶｅｒ＜Ｉｅ
ｎＣｅ　　ａｎｄ　　ＡＣＣｅｌｅｒａｔｉＯｎＯｆＶ
ｏｉｄ　　Ｉ　ｔｅｒａｔｉｏｎｓ　　ｉｎ　　Ｂｏｉ
ｌｉｎａ　　ＷａｔｅｒＲｅａｃｔｏｒ　　Ｃｏｒｅ　
　Ｃａ１ｃｕｌａｔｉｏｎｓ　ｕ　、　ＮＬＩＣＩ　。

Ｓｃｉ、　Ｅ　ｎａ、、　６４，７２４−７３２　（１
９７７）　］　Ｊに記載されている。出力分布計算装置
８はこのようにして棹出された炉内出力分布を最良出力
分布推定装置９に伝送する。

最良出力分布推定装ぎ９は、前記データサンプラー７の
伝送するＴＩＰおよびＬＰＲＭの出力信号および出力分
布計算装置８の口出した炉内出力分布を取込み、以下の
ようにしてまず最良ＴＩＰ測定値を推定する。すなわち
、 まず、前記出力分布計算装置８から伝送された炉内出力
分布から、以下の方程式（Ｉ［［）によりこれに対応す
るＴＩＰ測定値に対する３１ｒ１値を算出する。

ＥＲＣ（Ｋ、Ｌ） 一８ｕｍＰ　（Ｋ、１．Ｊ）ｘＣ（Ｋ、ｌ、、Ｊ＞−・
・・・・・・・（Ｉ［Ｉ） ここで、 ＥＲＣ（Ｋ、Ｌ）：炉心軸方向ノード位置Ｋ、ストリン
グしくこれを測定 位δ（Ｋ、Ｌ）と略称する） におけるＴ　［Ｐｉｌｌｌｌ定値に対 する計算値 ３　ｕＬＩ：測定位置（Ｋ、Ｌ）に隣接する全ノードに
関する和 Ｃ（Ｋ、Ｉ、Ｊ）：　１．Ｊ、に位置燃料の無限格子系
の中央位置にあるＴ 線束検出器の定格出力時の 測定値に対する計算値 ところで、上記計ＤＩＧ　（Ｋ、１．Ｊ）は例えば次の
ようにして厳密に求めることができる。寸なわら、通常
行われる無限格子系の中−燃料集合体核定数計鋒から得
られる各燃料棒又はチトンネルボックス等の各領域に於
ける各核種の核分裂反応数又は吸収反応数にその核種の
単位核分裂又は吸収及応当りに発生するγ線をエネルギ
ー依存ぐ求める。そのγ線源を使用して単−燃料集合体
無限格子系の中央位置にγ線束検出器を置いた体系のγ
線輸送計算を行い、γ線束測定器が単一種類の燃料集合
体で囲まれた場合の測定値を篩出しτおき、その結果を
その燃料集合体タイプの定数として燃焼度及びボイド率
でフィッティングしで記憶しておく。一方、各ノードの
燃焼度及びボイド率と前記のフィッティング係数を用い
て計算した各ノードの定格出力に対するγ線束測定器の
測定値に対する計算値を、出力分布計ｎ装置により得ら
れたそのストリングを囲む各燃料集合体の出力を重みに
して平均する事により、そのストリングを囲む燃料集合
体の中央位置に存在するγ線束測定器の測定値に対する
計算値が求められる。

ここで、ストリングを囲む燃料の出力を重みとづる理由
は各燃料集合体の出力の大小でγ線源のレベルが線形に
上下する事と各γ線源がそれによって生じる各領域のγ
線束に対して独立であるというγ線の特徴によるもので
ある。また多数の単一燃料集合体無限格子のγ線輸送計
算の結果を利用すれば、γ線輸送計算を行わなくても各
燃料棒又は領域のγ線源とガドリニア棒のガドリニア濃
度や位置及び燃料集合体の平均濃縮度やボイド率等のパ
ラメータから簡易計算によりγ線束検出器の測定値が筒
易的に算出できるので、通常の場合は単一燃料集合体無
限格子体系のγ線輸送計算の必要はなくなる事となる。

以上説明したようにして、炉心内の全測定位置に対する
γ線束測定器の測定値に対する計算値ＥＲＣが得られる
が、これらは以下で説明する炉心内γ線束測定器の出力
信号である測定値とは全く独立に得られるものである。

次に、最良出力分布推定装置９はデータサンプラー７か
ら伝送されたＴＩＰないしＬＰＰＭの出力信号からその
炉心軸方向（７方向）成分を抽出する。まず、ＶＩＰの
場合について説明する。

測定位置（Ｋ、Ｌ）におけるＴＩＰの出力信号をＥＲＭ
　（Ｋ、Ｌ）とすると、その炉心軸方向成分の抽出は、
以下のＮＶ）式及び（Ｖ）式を実行することにより得ら
れる。

ＥＲＭＺ　（Ｋ、Ｌ） ＝ＥＲＭ　（Ｋ、Ｌ）／ＥＲＭＲ（Ｌ）　・　（ＩＶ）
ＥＲＭＲ（Ｌ） ・・・・・・（Ｖ） ここで、 ３　ｕｌ：全炉心軸方向側定位ｉＫに関する和に 上記の（ｒＶ）式及び（Ｖ）式から次の（Ｖｌ）式が得
られる。

ＥＲＭ　（Ｋ、Ｌ） ＝ＥＲＭＲ（Ｌ）・ＥＲＭＺ　（Ｋ−、Ｌ）・・・（Ｖ
ｌ）この（Ｖｌ　）式をみればＴＩＰの出力信号は２個
の成分ＥＲＭＲとＥＲＭＺに分割できることが分る。し
かしてこの（Ｖｌ　）式においてＥＲＭＲ（Ｌ）は、着
目するストリングＬの仝軸方向測定位置に４３ける出力
信号を相加平均したものであるので、そのスｌ−リング
における出力信号の平均的な大きさを表すものである。

また、ストリングが炉心径り向にほぼ一様に設置されて
いることから、上記ＥＲＭＲ（Ｌ）はＴＩＰ測定値の炉
心径方向分布を表す成分であると解釈することができる
。一方、１ヨＲＭＺ　（Ｋ、Ｌ）はもとの出力信号ＥＲ
Ｍ　（Ｋ。

し）をそのストリングの平均値ＥＲＭＲ（Ｌ）で割った
ものであるから、上記の意味でのＴ　Ｉ　Ｐ測定値の炉
心径方向成分に関する情報を失っており、そのス１−リ
ングにおけるＴＩＰ測定値の炉心軸方向成分のみを表す
隋となっている。ここでは、この操作によりＥＲＭＺ　
（Ｋ、Ｌ）をｎ出することによりＴＩＰの出力信号たる
ＴＩＰ測定値の炉心軸方向成分を抽出するものである。

次に、最良出力分布推定装置９は、式（ＩＩＩ＞より算
出されるＴＩＰ測定値に対する計算値ＥＲＣ（Ｋ−、Ｌ
）から、次の（Ｖｌ）式及び（■）式を実行することに
より炉心径方向成分を抽出することができる。すなわら
、 ＥＲＣＺ（Ｋ、Ｌ） ＝ＥＲＣ（Ｋ、Ｌ）／ＥＲＣＲ（Ｌ）・・・（Ｖｌ　）
ＦＲＣＲ（Ｌ） ＝Ｓｕｎ＋Ｅ　ＲＣ（Ｋ、　ｌ−）　／Ｓｕｍ［１、Ｏ
］　−に　　　　　　　　　　Ｋ ・・・・・・（■） これらの（ＶＩ）式及び（■）式から、次の（ＩＸ）式
が得られる。

ＥＲＣ（Ｋ、Ｌ、） ＝ＥＲＣＲ（１）・ＥＲＣＺ（Ｋ、ｌ−）・・・（ＩＸ
）この（ＩＸ）式をみると王ＩＰ測定値に対する計算値
を２つの成分ＥＲＣＲ（Ｌ）　とＥＲＣＺ　（Ｋ。

Ｌ）とに分割できることが分る。ここでＥＲＣＲ（Ｌ）
は、上記ＥＲＭＲ（１−）がＴＩＰ３１１１定値の炉心
径方向分布を表す成分であると解釈できるのと同様にＴ
ＩＰ測定値に対する計算値の炉心径方向分布を表す成分
であると解釈することができる。

したがって、ここでは上記（■）式によりＥＲＣＲ（し
）を算出することをもって、Ｔ［Ｐ３１１定値に対する
計算値の炉心径方向成分を抽出するものである。

次に、最良出力分布推定装置９は次の（Ｘ）式により最
良ＴＩＰ測定値ＥＲＢ（Ｋ、Ｌ）を推定する。すなわち
、 ＦＲＢ　（Ｋ、Ｌ） ＝ＥＲＣＲ（Ｌ）・ＥＲＭＺ　（Ｋ、Ｌ）・・・（Ｘ）
この（Ｘ）式から分るように、最良ＴＩＰ測定値ＥＲＢ
　（Ｋ、Ｌ）は、ＴＩＰ測定値に対する計算値からその
炉心径方向成分ＥＲＣＲ（Ｌ）をとり、一方ＴＩＰの出
力信号であるＴＩＰ測定値からその炉心軸方向成分ＥＲ
ＭＺ　（Ｋ、Ｌ）をとり、両成分の積により算出される
。ここで、最良ＴＩＰ測定値とは、測定に関わる測定誤
差が全く存在しないと仮定した場合に得られるであろう
ところの理想的なＴＩＰ出力信号にできるだけ近いもの
という意味である。

ところで、上記の最良ＴＩＰ測定値の推定に当り、ＴＩ
Ｐの出力信号からその炉心径方向成分のみを使用する理
由は、ＴＩＰによる測定の誤差がＴＩＰないしその導管
の炉心径方向各側定設定点での偏位ではほとんど生じな
いことによる。このことを第３図を用いて説明する。す
なわち、ＴＩＰによるγ線束測定の際には、燃料集合体
１１間のギャップ１２の正確な中心位置１５を通過する
ようにＴＩＰを案内しなければならないが、このことは
ＴＩＰの導管が原子炉の運用によって変異・変形を起し
たり、またＴＩＰ自体の導管内での位置決めの技術的困
難さなどにより実際上困難であり、中心からずれた位置
１６を直線的にまたは中心に対して湾曲した位置１７を
案内されるＴＩＰが不可避的に存在することになる。し
かしながら、γ線束はコーナーギャップ内で径方向に関
してほとんど平坦な形であるため、この誤差は上記の位
＠１６または１７をＴＩＰが案内されてもそのストリン
グにおけるＴＩＰ測定値の炉心軸方向成分ＥＲＭＺ（Ｋ
、Ｌ）は正しいものにほぼ一致していることになる。

一方、前記の最良ＴＩＰ測定値の推定に当り、ＴＩＰ測
定値に対する計算値からその炉心径方向成分ＥＲＣＲ（
Ｌ）をとり、その炉心軸方向成分ＥＲＣＺ　（Ｋ、Ｌ）
をとらない。その理由は、出力分布演ｆ１装首８の内蔵
する物理モデルにより算出される出力分布は、その炉心
軸方向成分の８１ｆ３精度が不十分であっても、その炉
心径方向成分が精度良く計算されることが実験的にｉ認
されたからである。すなわち第４図は、横軸を炉心の中
心から燃料集合体への距離とし、縦軸を各燃料集合体に
おける平均のＬａ１４０濃度の誤差とした時に、照射後
のガンマ・スキャン実験による実験値に対する物理モデ
ルにより針筒された出力分布から変換される等価なＬａ
１４０１度の誤差を点によって示したものである。図か
ら分るようにこれらの誤差は、実用上はとんど問題が無
く、十分に無視できるものである。このことは、物理モ
デルによる出力分布露１算値の炉心径方向成分がよい精
度を持つことと同時に、出力分布に対する一次変換の（
［［）式によって締出されるＴＩＰ測定値に対する計算
値の炉心径方向成分ＥＲＣＲ（Ｌ）がよい精度を持つこ
とを裏づけでいる。

以−Ｌの理由により、最良Ｔ　Ｉ　Ｐ３１１１定値の推
定に当り、ＴＩＰ測定値に対するｇＩｎ値からその炉心
（￥方向成分ＥＲＣＲ（＋−）をとり、その炉心軸方向
成分ＥＲＣＺ　（Ｋ、Ｌ）をとらｒＺいこととする。

次に、最良出力分布推定装置９は、ノードの出力密度が
上記で推定された最良ＴＩＰ測定値に適合するように出
力分布計算値を下記（ＸＩ）式及び（Ｘｎ）式により修
正する。すなわら、Ｇ（Ｋ、Ｉ、Ｊ） ＝ＥＲＢ　（Ｋ、Ｌ）／ＥＲＣ（Ｋ、Ｌ）・・・・・・
・・・（ＸＩ） ＰＣ（Ｋ、Ｉ、Ｊ） ＝Ｇ　（Ｋ、１．Ｊ）・Ｐ（Ｋ、Ｉ、Ｊ）・・・・・・
・・・（ＸＩＩ） ここで、ノード（Ｋ、１．、））は測定位置（Ｋ。

Ｌ）に隣接するもの全てを指す。また（ＸＩＩ）式によ
り針環されるＰＣ（Ｋ、Ｉ、Ｊ）を（ＩＩＩ）式のＰ（
Ｋ、１．Ｊ）に代入することにより、補正因子Ｇ（Ｋ、
１．Ｊ）の定義と（ＩＩＩ）式とから明らかなように、
ＥＲＣ（Ｋ、Ｌ）として最良ＴｌＰ測定値ＥＲ［３（Ｋ
、Ｌ）を得る。この意味でＰＣ（Ｋ、Ｉ、Ｊ）は最良Ｔ
ＩＰ測定値に適合すると言う。

１ｄ模に、最良出力分布推定装置９は上記で得られたＰ
Ｃ（Ｋ、Ｔ、Ｊ）を、その径方向成分が物理モデルによ
る出力分布計算値のそれに一致するように補正する。す
なわち、 ＰＲ（１，Ｊ） ・・・・・・（Ｘｌ［Ｉ） ＰＺ　（Ｋ、Ｉ、Ｊ） ＝Ｐ　（Ｋ、Ｉ、Ｊ）／ＰＲ（ｒ、Ｊ）・・・（ＸＩＶ
）によりＰ（Ｋ、Ｉ、Ｊ）を炉心軸方向成分ＰＺ（Ｋ、
１．Ｊ）と炉心径方向成分ＰＲ（１，Ｊ）に分割する。

ここで、前記のように物理モデルの計算する出力分布の
炉心径方向成分の精度が良いことを考慮し、上記ＰＣ（
Ｋ、Ｉ、Ｊ）の径方向成分がＰＲ（Ｔ、Ｊ）に一致する
ように下記（ＸＶ）式で補正する。

ここで、上記（ＸＶ）式で計算されるＰＢ　（Ｋ。

１、Ｊ）が本発明の最良出力分布推定装置９の最終的に
求めるものであり、これをＴＩＰの出力信号および物理
モデルの計算値各々からの精度の良い成分のみを抽出し
て推定されたものという意味で最良出力分布と呼ぶ。

なお、上記実施例では炉心内移動型γ線束測定器である
ＴＩＰの測定が行われた場合について述べているが、Ｌ
ＰＰＭの場合には以下の様にしてＴＩＰの測定値が求め
られ、上記で述べたことが全くそのまま行われる。

すなわち、第３図において、同一測定点におけるＬＰＰ
Ｍ３とＴＩＰ４の両側定値は、測定対象の違いまたは検
出効率の違い等により一致しない。

しかしながら、両者を一致さけるような係数をＬｐ　Ｒ
Ｍ　３に対して移動型γ線束測定器の測定の際に測定点
毎に求める事ができる。従って各測定点での各測定設定
点におけるγ線束の測定値は、そのコーナーギャップ内
での現時点に４３ける熱中性子束測定器の測定値にこの
係数を掛けたものと、熱中性子束測定器の配置された位
置における現時点に一番近い時点でのγ線束測定器の測
定値との差を内挿し、γ線束測定器の測定値に加える事
により求める事ができる。

［梵明の効果］ 以上説明したように、本発明の原子炉のオンライン出力
分布監８！装置によれば、原子炉内の出力分布をオンラ
インで高精度で監視りることができるので、運転マージ
ンを減少させることかできる。

したがって、原子炉運用の効率を従来装置に比べてより
一層向上させることができる。

【図面の簡単な説明】

第１図は、本発明の一実施例のブロック構成図、第２図
は本発明が適用される炉心内γ線束測定器の配置を示す
τ１視図、第３図は第２図の炉心内γ線束測定器の偏位
を説明するための図、第４図は本発明による物理モデル
により計口される出力分布の炉心径方向成分の計口粘度
を説明するだめの図である。 ３・・・固定型中性子測定器 ４・・・可変型α線束測定器 ５・・・炉心現状データを測定Ｊる測定器６・・・出力
分布監視装置 ７・・・データサンプラー ８・・・出力分布監視装置 ９・・・最良出力分ｔｆｉ推定装置 １０・・・入出力装置 （８７３３）代理人　弁理士　猪　股　祥　晃くほか１
名） 第２図

Claims (1)

【特許請求の範囲】
1. 炉心現状データ測定器と、この炉心現状データ測定器か
   らの測定値を入力し予め内蔵された物理モデルに基づき
   炉心内の出力分布を算出する出力分布演算装置と、炉心
   内γ線束測定器と、この炉心内γ線束測定器の測定値を
   入力してその炉心軸方向成分のみを抽出するとともに前
   記出力分布演算装置で算出された出力分布を入力して得
   られた炉心軸方向成分を前記炉心内γ線束測定器の測定
   値の炉心軸方向成分により補正することにより最良出力
   分布を推定する最良出力分布推定装置と、を備えたこと
   を特徴とする原子炉のオンライン出力分布監視装置。