JPH04252998A - 原子炉の安定性予測方法及びその装置 - Google Patents

Abstract

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

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は、沸騰水型原子炉におけ
る燃料バンドルの熱水力学的振動（チャンネル安定性）
並びに原子炉全体の反応度帰還効果による炉心出力振動
（炉心安定性）の現状監視及び将来予測を行なう原子炉
の安定性予測方法及びその装置に関する。

【０００２】

【従来の技術】一般に、沸騰水型原子炉の炉心内に配置
された燃料バンドルを通る冷却材は水と蒸気の二相流を
なして、炉心内に発生する熱を除去するが、この二相流
には冷却材中の蒸気泡の量（ボイド量）、圧力損失、流
量間のフィードバックに基づく熱水力学的振動の可能性
があることが知られている。炉心内には多数の燃料チャ
ンネルがあり、これらは並行流路を形成するので、各燃
料チャンネルにとってはその熱水力学的振動は入口、出
口間の差圧一定の条件の下で他の燃料バンドルと一応独
立に生ずるものと考えてよい。このような振動モードを
“チャンネル安定性”と呼ぶ。

【０００３】一方、燃料チャンネルの二相流のフィード
バック現象に炉心の核特性及び再循環流路の動特性が関
係する振動のモードを“炉心安定性”と呼んでいる。通
常、炉心安定性の振動は炉心全体での中性子束（或いは
出力）の振動として現われ、チャンネル安定性の振動は
、振動した燃料チャンネル近傍の中性子束の振動として
現われる。

【０００４】このような振動状態が持続する状態を発振
状態といい、燃料の破損の原因にもなるため、原子炉の
安全上からは、発振状態を避ける必要がある。

【０００５】

【発明が解決しようとする課題】ところで、実際の原子
力発電所の運転時の炉心安定性に関する発振状態は、炉
心平均出力検出器信号、再循環流量等の全体的信号によ
って検知することは可能であるが、チャンネル安定性に
ついては、それが局所的な現象であるため、炉心内に多
数の局所出力検出器が設置されているものの、炉心平均
出力検出器信号のように常時監視に用いられているわけ
ではないので、発振を検知することは困難である。

【０００６】このように、原子炉の発振状態は、計装系
によりある程度検出することは可能であるが、発振状態
なる以前に、安定性余裕を適確に把握することは困難で
ある。

【０００７】一方、現在の運転計画では、実際の運転中
に再循環ポンプがトリップして不安定な状態になったと
きにも、安定性余裕を確保することが要求されているが
、現実の運転状態以外の状態の安定性余裕を知ることは
できないという問題がある。

【０００８】さらに、今後の運転計画を立案する場合に
も、安定性余裕を確保する必要があり、将来の炉心状態
の安定性余裕を正確に予測することが重要となっている
。

【０００９】本発明は、このような点を考慮してなされ
たもので、原子炉運転中に現状の安定性余裕を評価する
ことができるとともに、現状の運転状態から、外乱が発
生したときの炉心状態での安定性余裕を予測でき、また
今後の運転計画における炉心状態での安定性を予測でき
る原子炉の安定性予測方法及びその装置を提供すること
を目的とする。

【００１０】

【課題を解決するための手段】本発明に係る原子炉の安
定性予測方法は、前記目的を達成する手段として、原子
炉の炉心内出力分布を、３次元核熱水力計算モデルを用
いて計算するとともに、原子炉内の局所出力検出器また
は移動式炉心内計装系で検出される中性子束分布を用い
て前記炉心内出力分布の補正因子を決定し、この補正因
子により前記炉心内出力分布を補正して現状及び将来の
炉心内出力分布を求め、求められた炉心内出力分布を用
いて安定性解析を行ない、原子炉の現状及び将来の安定
度を予測するようにしたことを特徴とする。

【００１１】また、本発明に係る原子炉の安定性予測装
置は、前記目的を達成する手段として、原子炉の炉心内
出力分布を３次元核熱水力計算モデルを用いて計算する
炉心現状監視計算部と；原子炉内の局所出力検出器また
は移動式炉心内計装系で検出される中性子束分布を用い
て炉心内出力分布の補正因子を決定し、この補正因子を
前記炉心現状監視計算部に与えて炉心内出力分布を補正
させる補正因子計算部と；補正後の炉心内出力分布の入
力により、現状の炉心状態の安定性余裕を計算する安定
性現状監視計算部と；３次元核熱水力計算モデルを用い
て将来の炉心内出力分布を計算するとともに、前記補正
因子計算部からの補正因子により、将来の炉心内出力分
布を補正する炉心状態予測計算部と；補正後の将来の炉
心内出力分布の入力により、将来の炉心状態の安定性余
裕を計算する安定性予測計算部と；をそれぞれ設けるよ
うにしたことを特徴とする。

【００１２】

【作用】本発明に係る原子炉の安定性予測方法において
、３次元核熱水力計算モデルを用いて計算された原子炉
の炉心内出力分布は、原子炉内の局所出力検出器または
移動式炉心内計装系で検出した中性子束分布を用いて決
定される補正因子により補正され、現状および将来の炉
心内出力分布が求められる。そして、この炉心内出力分
布を用いて安定性解析が行なわれ、原子炉の現状および
将来の安定度が予測される。このため、精度の高い予測
が可能となる。

【００１３】また、本発明に係る原子炉の安定性予測装
置においては、炉心現状監視計算部により、３次元核熱
水力計算モデルを用いて原子炉の炉心内出力分布が計算
されるとともに、補正因子計算部により、原子炉内の局
所出力検出器または移動式炉心内計装系で検出した中性
子束分布を用いて炉心内出力分布の補正因子が決定され
、この補正因子により前記炉心内出力分布が補正される
。そして、この補正後の炉心内出力分布の入力により、
安定性現状監視計算部において、現状の炉心状態の安定
性余裕が計算される。このため、高精度の安定性余裕が
計算される。

【００１４】一方、将来の炉内出力分布は、炉心状態予
測計算部において、３次元核熱水力計算モデルを用いて
計算されるとともに、前記補正因子により将来の炉心内
出力分布が補正される。そして、この補正後の炉心内出
力分布の入力により、安定性予測計算部において、将来
の炉心状態の安定性余裕が計算される。このため、高精
度の安定性余裕が計算される。

【００１５】

【実施例】以下、本発明の一実施例を図面を参照して説
明する。

【００１６】図１は、本発明に係る原子炉の安定性予測
装置の一例を示すもので、図中、符号１は炉心２が格納
された原子炉であり、前記炉心２に設置されている局所
出力検出器または移動式炉心内計装系等の各検出器から
は、制御棒パターン、炉心流量、出入口温度、炉心圧力
等のプロセス量が出力されるようになっている。そして
、このプロセス量は、炉心性能計算装置３に入力される
ようになっている。

【００１７】炉心性能計算装置３は、炉心現状監視計算
部４と、補正因子計算部５と、炉心状態予測計算部６と
を備えており、前記炉心現状監視計算部４は、前記プロ
セス量の入力により、ヒートバランスから炉心熱出力を
計算するとともに、３次元核熱水力計算モデルを用いて
炉心内出力分布を計算するようになっている。

【００１８】また、前記補正因子計算部５は、炉心内の
局所出力検出器または移動式炉心内計装系で検出された
中性子束分布を使って、計算モデルで得られる出力分布
計算結果に対する補正因子を求めるようになっている。 具体的には、検出された中性子束分布と計算モデルで得
られる中性子束分布との比として補正因子を求め、これ
を記憶するようになっている。そして、この補正因子計
算部５で求められた補正因子は、前記炉心現状監視計算
部４に送られ、この炉心現状監視計算部４において計算
モデルで求めた炉心内出力分布計算値に補正因子を掛け
ることにより、現状の炉心内出力分布が求められるよう
になっている。

【００１９】また、前記炉心状態予測計算部６は、将来
の時点、例えば再循環ポンプトリップの外乱が発生した
時点や今後の運転計画点に対し、炉心状態の予測計算を
実行するようになっている。すなわち、炉心状態予測計
算部６は、予測しようとする運転点に対する炉心データ
（制御棒パターン、炉心熱出力、炉心流量等）を与える
ことにより、前記炉心現状監視計算部４と同一の３次元
核熱水力計算モデルを用いて炉心内出力分布を計算する
とともに、補正因子計算部５に記憶されている補正因子
を掛けることにより、将来の炉心内出力分布を精度よく
求めることができるようになっている。そして、この将
来の炉心内出力分布および前記炉心現状監視計算部４で
求められた現状炉心内出力分布は、安定性計算装置７に
与えられるようになっている。

【００２０】この安定性計算装置７は、安定性現状監視
計算部８と、安定性予測計算部９とを備えており、安定
性現状監視計算部８は、現状の炉心状態の安定性余裕を
計算するとともに、安定性予測計算部９は、将来の炉心
状態の安定性余裕を計算するようになっている。

【００２１】すなわち、前記安定性現状監視計算部８は
、炉心現状監視計算部４で求めた炉心内出力分布の入力
により、また前記安定性予測計算部９は、炉心状態予測
計算部６で求めた炉心内出力分布の入力により、炉心安
定性とチャンネル安定性とを計算するようになっている
。

【００２２】次に、本実施例の作用について説明する。 チャンネル安定性は、二相流体のチャンネル内のスイー
ブ時間遅れ等によって生じるチャンネル流量の不安定を
取扱うものであり、局部の安定性を調べることを目的と
するものである。

【００２３】このチャンネル安定性は、入力された炉心
内出力分布のうち、最も厳しくなる単一チャンネルにつ
いて、熱水力特性モデルを解いてラプラス変換し、チャ
ンネル全体の開ループ伝達関数を求め、これから周波数
応答法により安定余裕を求める。

【００２４】一方、炉心安定性は、熱水力学的に特性の
異なる多数のチャンネルが炉心内に存在していて、炉心
出力が各部で変動し、流量の再配分が起こって炉心の特
性が各部で変化するような場合に、炉心出力の安定性を
調べることを目的とするものである。

【００２５】この炉心安定性は、入力された炉心内出力
分布を用い、熱水力学的にみて類似したいくつかのチャ
ンネルを１つのグループにまとめて、炉心全体をいくつ
かの領域に分け、各領域の熱水力特性モデルを解いてラ
プラス変換し、周波数応答法により安定余裕を求める。

【００２６】ところで、チャンネル安定性及び炉心安定
性の安定度を表わす指標としては、次の減幅比を用いる
。減幅比は、出力振動が起こった場合の第１波と第２波
との振幅の比として定義する。減幅比が１以上になると
、振動は持続または増幅することを示し、いわゆる発振
状態を示す。したがって、減幅比が１より小さいことが
安定であることを示し、減幅比が小さいほど安定性余裕
が大きいことを示す。

【００２７】そこで、プラント運転中に、予期されるあ
らゆる運転状態に対し、減幅比が１より小さくなるよう
に設計する必要がある。また、実際の運転に際しては、
充分な安定性余裕を確保するため、通常運転状態で、チ
ャンネル安定性の減幅比は０．５以下、炉心安定性の減
幅比は０．２５以下になることを目標基準とする。

【００２８】したがって、通常運転時及び運転計画作成
時には、安定性計算装置７で計算された減幅比が、前記
目標基準を満足していることを確認する。運転計画作成
時に、目標基準が満足できない場合には計画を変更する
。また、現状運転状態に対し、目標基準を満足していな
い場合には、再循環ポンプトリップ時の炉心状態を炉心
性能計算装置３で予測し、安定性計算装置７で安定性を
予測する。予測された減幅比が１よりも小さい場合には
、そのまま運転状態を継続するが、減幅比が１以上と予
測される場合には、運転状態を変更する処置をとる。 通常、制御棒を挿入して炉心流量を増加させると、安定
性余裕は増加する。

【００２９】なお、本発明は、前記実施例に限定される
ものではなく、例えば以下のような変形例も考えられる
。

【００３０】まず、炉心性能計算装置３における３次元
核熱水力計算モデルを用い計算された炉心内出力分布の
計算誤差補正方法の他の例としては、以下の方法が考え
られる。

【００３１】すなわち、３次元核熱水力計算モデルの拡
散方程式に含まれるバックリング項や無限増倍率を修正
し、計算される中性子束分布が、炉心内の局所出力検出
器または移動式炉心内計装系で測定された中性子束分布
に一致するようにするものであり、得られたバックリン
グ項や無限増倍率の補正量を記憶しておき、予測計算に
も使用するものである。

【００３２】また、安定性計算装置７の安定性計算モデ
ルとしては、次のものが考えられる。すなわち、安定性
計算モデルは、一般に動特性モデルを解く必要がある。 前記実施例では、動特性方程式を解く際に、炉心内出力
分布等を簡略化しているが、そのモデル化の精度を上げ
る必要がある場合には、例えば動特性方程式を解く際に
、炉心を１次元化、２次元化、３次元化と上げることが
考えられ、次元を上げることにより、計算精度はよくな
る。ただし、動特性方程式を解くためには、計算時間が
かかるため、計算機の能力と許容できる計算時間との関
係で、モデル化の次元数を決定する必要がある。

【００３３】さらに、安定性計算装置７の安定性計算方
法を、前記実施例よりもさらに簡略化し、計算時間を短
縮する方法としては、次の方法が考えられる。

【００３４】すなわち、安定性計算装置７において、炉
心性能計算装置３で得られる現状及び将来の炉心内出力
分布を、炉心内の燃料バンドルの炉心平均出力に対する
ピーキング係数の二乗平均値と炉心平均軸方向ピーキン
グ係数とその位置との３つのパラメータに整理し、その
パラメータと原子炉熱出力及び炉心流量の関数として炉
心安定性を予測し、チャンネル安定性については、炉心
性能計算装置３で得られる現状及び将来の炉心内出力分
布を、炉心内の燃料バンドルの軸方向ピーキング係数の
最大値と最大出力バンドルの軸方向ピーキング係数とそ
の位置との３つのパラメータに整理し、そのパラメータ
と原子炉熱出力および炉心流量の関数として予測するよ
うにしたものである。なお、この場合も、炉心安定性及
びチャンネル安定性の安定度を表わす指標としては、減
幅比を用いる。

【００３５】次に、チャンネル安定性減幅比及び炉心安
定性減幅比の評価方法について説明する。一般に、原子
炉の安定性に影響を及ぼすパラメータとしては、数多く
のものが考えられる。炉心安定性は、熱水力学的な要因
が反応度フィードバックすることにより、炉心内で空間
的な出力振幅を起こすもので、炉心出力が各部で変動し
、各バンドルの流量の再配分が起こって、炉心の特性が
各部で変化することにより起こるため、炉心出力、流量
及び径方向出力分布が炉心安定性に大きな影響を与える
。

【００３６】すなわち、径方向の出力分布のバラ付きが
大きい程、炉心安定性の安定度は悪化する。このことは
、径方向の出力分布の大小の差が大きい程、各バンドル
の流量の再配分が起こり易くなり、それが核的なフィー
ドバックとなって炉心内の出力の振動が起こり易くなる
結果、炉心安定性の安定度が悪化するためと考えること
ができる。

【００３７】一方、チャンネル安定性については軸方向
のボイド分布が特に影響を与える。すなわち、炉心内に
単相で流れこんできた冷却材がやがて加熱されて気液二
相となるが、気液二相となっている領域が長い程チャン
ネル安定性は悪化する。チャンネル安定性に影響を与え
るパラメータとしてはまず出力と流量が考えられる。出
力が高い程、流量が低い程チャンネル安定性の安定度は
低下する。また、これら以外でチャンネル安定性に大き
な影響を及ぼすものとして、軸方向の出力分布形があげ
られる。

【００３８】すなわち、この軸方向の出力分布形はチャ
ンネル安定性に大きな影響を及ぼす上に、実際の運転状
態において、制御棒の出入れや燃焼に伴なう核種ゼノン
の効果等によってかなり変化し、かつ、燃料チャンネル
毎に異なった形になっているので、現実の運転状態に即
した安定度評価を行なうためには、軸方向の出力分布形
の効果を考慮する必要がある。

【００３９】ところで、軸方向の出力分布形とチャンネ
ル安定性の安定度との関係は、一般的には、軸方向の出
力分布形のピークの高さが高ければ高い程、かつ、ピー
ク位置が炉心（燃料バンドル）の下部であればある程、
安定度が小さくなるということができる。

【００４０】以上に示したような性質を反映して、チャ
ンネル安定性および炉心安定性には、図２ないし図６に
示すような傾向があることが知られている。

【００４１】炉心安定性は炉心全体での現象であること
から図２及び図３に示すように、炉心内の燃料バンドル
の炉心平均出力に対するピーキング係数の二乗平均値（
Ｒ値）、炉心平均軸方向平均出力（Ｐａ　）及びその発
生軸方向ノード（Ｎ）について傾向があり、Ｒが小さい
程、またＰａ　が大きくＮが小さい程安定である。言い
換えれば、炉心内の径方向出力分布がより平坦で、炉心
平均の軸方向出力ピーキングが大きく、その発生位置が
より上部にある程安定であることを示している。

【００４２】また炉心安定性は、原子炉の出力及び流量
に対して図６に示すような傾向があることが知られてい
る。すなわち減幅比の出力及び流量に対する変化は、１
次関数または２次関数で表現することができる。

【００４３】次にチャンネル安定性については、バンド
ル単位の局所的な現象であることから、図４及び図５に
示すように、最大出力バンドル（ホットチャンネル）の
径方向出力ピーキング（Ｐｒ　）、ホットチャンネル軸
方向出力ピーキング（Ｐａｈ）及びその発生ノード（Ｎ
）について傾向があり、Ｐｒ　が小さい程、またＰａｈ
が小さくＮが大きい程安定である。言い換えれば、ホッ
トチャンネルの出力が炉心平均出力により近く、またホ
ットチャンネルの軸方向分布がより平坦で、そのピーク
位置がより上部にある程安定であることを示している。 チャンネル安定性と原子炉出力、流量との関係について
は、炉心安定性の場合と同様に考えてよい。

【００４４】本実施例では以上の点を考慮に入れて、チ
ャンネル安定性減幅比（Ｄｃｈａｎ）についてはＰｒ　
、Ｐａｈ、Ｎを用い、炉心安定性減幅比（ＤＣＯＲＥ）
についてはＲ、Ｐａ　、Ｎを用いて次式で評価する。

【００４５】すなわち、炉心安定性については式１で評
価し、またチャンネル安定性については式２で評価する
。 〔式１〕 Ｄ′ＣＯＲＥ＝ａ０　＋ａ１　＊Ｒ−ａ２　＊Ｐａ　−
ａ３　＊Ｎ　　ＤＣＯＲＥ＝Ｄ′ＣＯＲＥ＋ｂ１　＊Δ
Ｐ＋ｂ２　＊ΔＰ２　＋ｂ３　＊ΔＷ＋ｂ４　＊ΔＷ２
　ただし、ΔＰ：基準出力からの熱出力変化量ΔＷ：基
準流量からの流量変化量 〔式２〕 Ｄ′ｃｈａｎ＝ｃ０　＋ｃ１　＊Ｐｒ　＋ｃ２　＊Ｐａ
ｈ−ｃ３　＊Ｎ　　Ｄｃｈａｎ＝Ｄ′ｃｈａｎ＋ｄ１　
＊ΔＰ＋ｄ２　＊ΔＰ２　＋ｄ３　＊ΔＷ＋ｄ４　＊Δ
Ｗ２　上式においてａ０　〜ａ３　、ｂ１　〜ｂ４　、
ｃ１　〜ｃ３　、ｄ１　〜ｄ４　を決定することにより
ＤＣＯＲＥ及びＤｃｈａｎを決めることができるが、こ
れらの係数は事前に予測計算を数ケース行い、その計算
結果を用いて容易にしかも精度よく決めることができる
。

【００４６】これらの式を予め決定しておくことにより
、現在並びに将来における安定性の度合はＰａ　やＲな
ど、炉心性能計算３により得られる値のみで検知できる
こととなり、従来のようにその都度安定性解析を行う必
要がなくなるため計算時間を大幅に短縮できる。また炉
心状態予測計算部６で求まった値を入力として用いれば
将来における安定性の度合も極めて短時間のうちに知る
ことができる。

【００４７】従って、通常運転時および運転計画作成時
に現在及び将来の減幅比を計算する場合、炉心性能計算
装置３及び炉心状態予測計算部６で得られた結果を用い
て極めて短時間のうちに精度よく行うことが可能となり
、計算量を大幅に削減できることができる。

【００４８】なお、本実施例では、安定性計算装置７の
安定性計算方法を簡略化し、計算時間を短縮する場合に
ついて示したが、このことについて他の実施例も考えら
れる。まず、炉心安定性減幅比とチャンネル安定性減幅
比を表現する式について、Ｐａ　、Ｒなどを含む項を２
次まで拡張するものである。例えば前述の実施例におけ
るＤ′ＣＯＲＥ及びＤ′ｃｈａｎは次のようになる。

【００４９】〔式３〕 　　Ｄ′ＣＯＲＥ＝ａ０　＋ａ１　＊Ｒ＋ａ１　′＊Ｒ
２　−ａ２　＊Ｐａ　−ａ２　′＊Ｐａ　２　　　　　
　　　　　　　　−ａ３　＊Ｎ−ａ３　′＊Ｎ２　　　
Ｄ′ｃｈａｎ＝ｃ０　＋ｃ１　＊Ｐｒ　＋ｃ１　′＊Ｐ
ｒ　２　＋ｃ２　＊Ｐａｈ＋ｃ２　′＊Ｐａｈ２　　　
　　　　　　　　　　−ｃ３　＊Ｎ−ｃ３　′＊Ｎ２　
このように２次まで拡張することによって精度を上げる
ことができるが、事前に求めるべき係数が増えるので、
２次まで拡張する項を別途検討する必要がある。またそ
の他の実施例としては、炉心安定性減幅比とチャンネル
安定性減幅比を表現する式に、安定性評価パラメーター
の積の項を導入することも考えられる。例えばＤ′ＣＯ
ＲＥを評価する時に、Ｒ＊Ｐａ　の項も含めて評価する
ことである。このことにより安定性評価パラメータ同士
が互いに相関をもって減幅比に与えていた効果も評価す
ることができ、炉心安定性減幅比とチャンネル安定性減
幅比を表現する式の信頼性をより高めることができる可
能性がある。

【００５０】

【発明の効果】以上説明したように、本発明に係る原子
炉の安定性予測方法及びその装置においては、３次元核
熱水力計算モデルを用いて計算した原子炉の炉心内出力
分布を、原子炉内の局所出力検出器または移動式炉心内
計装系で検出した中性子束分布から求められる補正因子
で補正して現状及び将来の炉心内出力分布を求め、この
炉心内出力分布を用いて安定性解析を行ない、原子炉の
現状及び将来の安定度を予測するようにしているので、
常に安定性を確保した運転が実現でき、原子炉の安定性
を向上させることができる。

【図面の簡単な説明】

【図１】本発明の一実施例に係る原子炉の安定性予測装
置を示すブロック図。

【図２】径方向分布指標と炉心安定性減幅比との関係を
示すグラフ。

【図３】炉心平均軸方向ピーキング及びその発生ノード
と炉心安定性減幅比との関係を示すグラフ。

【図４】ホットチャンネル径方向出力ピーキングとチャ
ンネル安定性減幅比との関係を示すグラフ。

【図５】ホットチャンネル軸方向出力ピーキング及びそ
の発生ノードとチャンネル安定性減幅比との関係を示す
グラフ。

【図６】原子炉出力及びポンプ速度と炉心安定性減幅比
との関係を示すグラフ。

【符号の説明】

１　　原子炉 ２　　炉心 ３　　炉心性能計算装置 ４　　炉心現状監視計算部 ５　　補正因子計算部 ６　　炉心状態予測計算部 ７　　安定性計算装置 ８　　安定性現状監視計算部 ９　　安定性予測計算部

Claims (2)

【特許請求の範囲】
1. 【請求項１】原子炉の炉心内出力分布を、３次元核熱水
   力計算モデルを用いて計算するとともに、原子炉内の局
   所出力検出器または移動式炉心内計装系で検出される中
   性子束分布を用いて前記炉心内出力分布の補正因子を決
   定し、この補正因子により前記炉心内出力分布を補正し
   て現状及び将来の炉心内出力分布を求め、求められた炉
   心内出力分布を用いて安定性解析を行ない、原子炉の現
   状及び将来の安定度を予測することを特徴とする原子炉
   の安定性予測方法。
2. 【請求項２】原子炉の炉心内出力分布を３次元核熱水力
   計算モデルを用いて計算する炉心現状監視計算部と；原
   子炉内の局所出力検出器または移動式炉心内計装系で検
   出される中性子束分布を用いて炉内出力分布の補正因子
   を決定し、この補正因子を前記炉心現状監視計算部に与
   えて炉心内出力分布を補正させる補正因子計算部と；補
   正後の炉心内出力分布の入力により、現状の炉心状態の
   安定性余裕を計算する安定性現状監視計算部と；３次元
   核熱水力計算モデルを用いて将来の炉心内出力分布を計
   算するとともに、前記補正因子計算部からの補正因子に
   より、将来の炉心内出力分布を補正する炉心状態予測計
   算部と；補正後の将来の炉心内出力分布の入力により、
   将来の炉心状態の安定性余裕を計算する安定性予測計算
   部と；を具備することを特徴とする原子炉の安定性予測
   装置。