JP7192150B2 - 加圧水型原子炉の管理方法および応分の管理システム - Google Patents

Description

本発明は、好ましくは原子力発電所内の加圧水型原子炉を管理する方法に関する。
本発明は同様に、応分の管理システムにも関する。

原子炉は、広い時間スケール内の高度の動力学を有する複雑な非線形系であり、その包括的制御は些細なものではない。
熱出力および冷却材温度という周知の制御以外に、原子炉の制御装置は、所与の限界内のみでの運転を可能にする運転上の安全性、燃焼の均一化、燃焼補償、毒作用の補償、出力密度分布の均一化、柔軟な発電の支援、運転経済性などの多くの他の側面を処理する。

世界中の原子炉の管理は現在、制御技術に基づいている。
しかしながら、このタスクのために使用される従来の制御技術は、その適用限界に達している。
このアプローチが有する主要な問題点としては、複合システムの逆問題を解く試み、原子炉毒物質の複雑で長時間スケールの動態およびそれらの空間分布、ならびに利用可能な動作主体の数よりも有意に多い制御目標数、がある。

したがって、本発明の目的は、加圧水型原子炉を制御または管理する方法を提供すること、および多数の制御目標を同時に考慮し平衡化することのできる応分の管理システムを提供することにある。
該方法および該応分のシステムは、異なる種類の既存のまたは新規の原子炉にとって実装および設定が容易なものである。
これらは、実時間制御または管理の能力を有するものである。

本発明によると、前記目的は、請求項１に係る方法によって達成される。
対応する管理システムは、請求項８で規定されている。

したがって、請求項１は、加圧水型原子炉を管理する方法において、
加圧水型原子炉が、炉心と、原子炉冷却材を保持する炉心用冷却回路とを含み、
原子炉の状態が、多数の測定可能な状態変数によって特徴付けされ、炉心の反応度が、多数の起動変数によって制御されかつ原子炉毒作用により影響され、
起動変数の所与の時間依存軌跡について、状態変数についての応分の軌跡が、状態変数の測定された現在の値、所与のまたは計算された毒作用値および反応度収支等式に基づいて予測される、
方法であって、
好ましくは１時間よりはるかに長く、典型的には２４時間前後の規模を有する将来の時間間隔について起動変数の多数の無作為に変動させられた起こりうる軌跡を、反復的に考慮するステップを含み、
起動変数の各軌跡には、起動変数、プロセス変数、および／またはこれらから導出された変数の、予め設定された条件または値によって特徴付けされる多数の事象または不利な炉心状態についての重み付けまたはペナルティ値を含む値表に基づいて、性能指数が割り当てられ、
起動変数軌跡は、性能指数が局所的極値を有するような形で選択され、
アクチュエータが相応して移動させられまたは設定される、
方法を規定している。

基本的概念のさらなる実施形態および進展は、従属請求項中および以下の説明の中で規定されている。

請求対象の方法および応分のデバイス／システム（「炉心ガバナ」の異名をもつ）は、制御技術に代って、コンピュータによる変分法（すなわち汎関数の最大値および／または最小値を求める）を使用する。
制御技術に基づく原子炉制御装置とは異なり、提案される炉心ガバナは、実時間でコマンドを制御するだけでなく、好ましくは２４時間という次の大きなタイムスパンにわたる制御動作についての包括的な計画をもコンパイルし、好ましくはそれを視覚化する。
多義の逆問題の解を求める代りに、本発明に係る炉心ガバナは、些細な順問題を何度も解き、各反復において制御動作についての計画を変動させ、所与の値表に基づいて性能指数を用いて最善の動作を探求する。
制御技術とは異なり、コンピュータによる変分法は、数および性質とは無関係に同時に全ての目標を組込む単一の性能指数を処理することから、制御目標数に制限が全く無い。

典型的な原子炉において、本発明が適切な形で平衡化し最適化することのできるこのような相反する目標の例としては、以下のものが含まれる。
－アキシャルオフセット（ＡＯ）の管理
－棒の運動の最小化
－ホウ酸消費量の最小化
－脱塩水消費量の最小化
－燃焼補償
－均一な燃焼

さらに、本発明は以下のことを提供し得る。
－全ての重要な炉心特性の長時間、好ましくは２４時間にわたる予測及び表示
－キセノンおよびサマリウム毒作用の監督および予測
－１００％までの急速出力上昇能力の監視および／または保証
－起こりうる全ての変化、例えば異なる負荷または炉心再構成のための些細な調整
－些細なパラメータ化
－起こりうる全ての加圧水型原子力発電所タイプに対する応用可能性

この根本的に新しいアプローチの対応する利点は、以下のことに関連するものである。
－安全性
－利用可能性
－操縦性（柔軟性ある運転）
－計画効率
－透明性
－操作性
－人間工学
－予測可能性

これにより、炉心の高度な負荷追従制御および他の制御体制を含めた、送電網サービスの完全自動化が可能となる。
しかしながら、原子炉運転の半自動または手動モードも同様に支援される。
詳細には、アクチュエータの運動は、手動式、半自動式または全自動式で行なうことができる。

提案される方法は、実時間よりはるかに高速のコンピュータ計算を必要とし提供する。
一連の好ましくは２４時間の軌跡を、瞬時に計算しなければならない。
提案されたアルゴリズムが、好ましくは２４時間の最適化タイムフレーム内で、好ましくは例えば各２００ｍｓについての反応度収支のみを処理することから、現代のコンピュータの出力はこれに十分なものである。
提案される方法は、このために単純な算術を必要とする。
炉心内の空間分布は、アキシャルオフセット（ＡＯ）だけが問題となるため、好ましくは、２点モデルにおいて考慮される。
２点モデルは、当該目的のために完全に十分なものであり、極めてわずかな計算能力しか必要としない。

提案される方法は、アルゴリズムがそれ自体未修正のままに留まることから、異なる加圧水型原子力発電所に対し容易に適応可能である。
原則として、アルゴリズムのパラメータのみ、更新する必要がある。
順問題のみが解かれることになるため、これらのパラメータは全て物理的で、各原子炉について周知である。
実際、これらのパラメータの大部分は、単純に反応度係数である。

提案される管理アルゴリズムの挙動は、値表中の定数値および値関数の修正によって、簡単に修正されることができる。
値表の補足によって、簡単に新しい特徴を埋め込むことができる。
そうすることで、性能指数アルゴリズムも同様に補足されるはずである。
性能指数は、異なる寄与の単なる合計であることから、この補足は些細なものである。

本発明のこれらのおよび他の態様は、以下の段落からより容易に明らかになるものであり、これらの段落では、本発明の例示的実施形態が添付図面を参照しながら論述されている。

加圧水型原子炉および本発明に係るその管理アーキテクチャの図式的概観を提供する。 本発明に係る「炉心ガバナ」機能のための管理アルゴリズムの図式的概観を提供する。 状態変数によってスパンされた、加圧水型原子炉のための状態空間の視覚化である。 起動変数によってスパンされた、加圧水型原子炉のための起動空間の視覚化である。 原子力発電所のための負荷計画の線図である。 図２に係る「炉心ガバナ」上で実行中のアルゴリズムにおいて使用するための例示的値表を示す。 図６に係る値表中で使用するための例示的価値関数を示す。 図６に係る値表中で使用するための別の例示的価値関数を示す。 図２に係る「炉心ガバナ」のディスプレイ上に示された例示的出力スクリーンである。 図２に係る「炉心ガバナ」のための好ましいアルゴリズムについてのフローチャートを示す。

図１は、加圧水型原子炉２および本発明に係る管理アーキテクチャを伴う原子力発電所４の簡略化された図式的概観を提供している。
本発明に係る炉心管理は、以下の説明から明らかであるように、従来の古典的制御理論を応用しないものの、管理アーキテクチャは、制御アーキテクチャと呼ぶことともできるものである。

例示的実施形態において、原子力発電所４は、加圧水型原子炉（ＰＷＲ）２を含む。
ＰＷＲは、原子炉圧力容器（ＲＰＶ）８の内部に核分裂性核物質を伴う炉心６を含む。
運転中、原子炉圧力容器８内部の核熱は、原子炉冷却材ポンプ（ＲＣＰ）１２の駆動力の下で一次冷却回路１０内を循環する一次冷却媒体（または冷却材）へと伝達される。
取付けられた蒸気発生器１４内で、一次冷却媒体の熱は、二次冷却回路１６内を循環する二次冷却媒体へと伝達され、こうして二次冷却媒体を蒸発させる。
冷却された一次冷却媒体は、その後、再び炉心６内へと導かれる。

こうして生産された生蒸気は、発電機２０に結合された蒸気タービン１８を駆動し、こうして、結び付けられた送電網２２内に補給すべき電気を発生させる。
蒸気タービン１８の下流側で、弛緩した蒸気は、凝縮器２４の中で凝縮され、その後再び給水として、給水ポンプ２６を用いて蒸気発生器１４まで導かれる。
二次冷却回路１６内部の給水タンク２８が、補償槽として作用し得る。

蒸気タービン１８内に入る蒸気の流量は、蒸気補給ライン３２内のタービンバルブ３０（ここでは簡略化のため１つだけが示されているが、通常は複数のものが並列で存在する）によって調整可能である。
いくつかの特殊なケース（プラントの操業開始、タービントリップ、所内負荷運転への切換え、など）においては、蒸気タービン１８に入らず、バイパス流量を設定するためにバイパスバルブ３６（ここでは簡略化のため１つだけが示されているが、通常は複数のものが並列で存在する）を含むバイパスライン３４を介して直接凝縮器２４に導かれる余剰の蒸気が存在する。
タービンバルブ３０およびバイパスバルブ３６の位置は、タービン制御装置４０およびバイパス制御装置４２を含む結び付けられた制御システムによって制御される。
制御ループ内に入る測定されたパラメータには、蒸気補給ライン３２内部の生蒸気圧力ｐ、蒸気タービン１８の回転速度ｎ、および／または発電機２０によって出力された電力Ｐが含まれ得る。

原子炉２の管理は、その炉心６内に挿入可能である多数の制御棒４４を介して可能である。
通常、加圧水型原子炉（ＰＷＲ、ＤＷＲ、ＷＷＥＲなど）内の制御棒は、制御アセンブリの形に統合される。
単一のアセンブリの複数の棒は、単一の棒駆動メカニズムによって駆動され、単一の燃料アセンブリの内部で共に移動する。
詳細には４つ（ＰＷＲ、ＤＷＲ）または６つ（ＷＷＥＲ）の対称的に位置設定されたアセンブリが、１つの制御グループを構成する。
これらのグループは、通常、さらに２つのバンクへと統合される。
小さい方のバンクＰ（ＰｏｗｅｒのＰ）が、原子炉出力を制御するために使用され、大きいバンクＨ（ＨｅａｖｙのＨ）は、原子炉の運転停止専用であり、正常運転中は活性ゾーンから完全に引き抜かれている。
ドイツ型およびフランス型などのいくつかの高度の制御概念においては、Ｈバンクは、同様に、活性ゾーン内の出力密度分布を調整するためにも使用される。
この目的で、Ｈバンクは、活性ゾーンの上部部分内にわずかに挿入され、いわゆる出力密度アキシャルオフセット（ＡＯ）制御装置によってわずかに移動させられ得る。
標準的吸収棒以外に、いくつかのフランス型原子炉は、特殊な「グレイ棒」を有する。
グレイ棒の吸収能力は、通常の（ブラック）棒よりも低い。
グレイ棒は、完全引抜き状態または完全挿入状態という２つの運転位置しか有さず、原子炉出力の大幅な段階で削減するために使用される。
図１に示されている例示的実施形態においては、ＰバンクおよびＨバンクが見られる。

典型的には、Ｐバンク位置用に結び付けられた制御システムは、原子炉の操業開始に使用されるφ－制御機構としても知られる中性子束制御装置４６、および出力運転を担当する平均冷却材温度（ＡＣＴ）制御装置４８を含む。
φ－制御機構は、典型的に炉外計装５０によって測定される炉心６内部の中性子束φについての測定値に依存する。
ＡＣＴ制御機構は、一次冷却媒体の温度についての測定値、詳細には、炉心６との関係における流入物温度Ｔ１および流出物温度Ｔ２から導出される平均冷却材温度（ＡＣＴ）に依存する。
ＰＷＲおよびＤＷＲとは異なり、ＷＷＥＲプラントは、この時点で、平均原子炉冷却材温度ＡＣＴではなく生蒸気圧力Ｐを制御する。
しかしながら、本発明に係る管理概念は、この事実によるさらなる影響を全く受けることがない。

さらに、炉心６内部の反応度は、一次冷却回路１０内のホウ素濃度の調整によって影響される。
このために、簡略化のためここでは明示的に示されていない理由から、いわゆる体積制御システム（ＶＣＳ）を介して、一次冷却回路１０に取付けられた、一方ではホウ酸用（ホウ素濃度を増大させる）そして他方では脱塩水用（ホウ素濃度を減少させる）のための補給システムが存在する。

合わせてＢＯＤＥとしても公知であるホウ酸および脱塩水の補給量、ひいては一次冷却回路１０内のホウ素濃度は、ＢＯＤＥ補給ライン６０、６２内の応分の補給ポンプ５２、５４および制御バルブ５６、５８を介して調整可能である。
補給ポンプ５２、５４および制御バルブ５６、５８は、ＢＯＤＥ制御装置６４によって設定される。

図１の例示的実施形態において、「負荷ガバナ」６６と呼ばれる上級管理モジュールは、スクリーン８０内および図５中で視覚化されているように送電網オペレータから負荷計画を受信し、そこには将来の電力需要計画が、通常は来る２４時間について、通常は１５分毎の時間の関数として示されている。
原子力発電所オペレータは、この負荷計画を編集し発出することができる。
発出された負荷計画を用いて、負荷ガバナ６６は、電力についての設定点を、実時間でタービン制御装置４０に提供する。
さらに、負荷ガバナ６６は、ＡＣＴ（またはｐ）制御装置４８を介してＰバンクのフィードフォワード制御を行なう。

炉心ガバナ６８と呼ばれる関連する上級管理モジュールが、負荷ガバナ６６から有効な負荷計画を、そして炉内計装７０から現在の出力密度を受信する。
この情報および出力履歴から導出された現在の原子炉の毒作用およびその空間分布についてのその独自の情報を用いて、炉心ガバナ６８は、以下でさらに詳細に説明するように、将来、好ましくは来る２４時間についての全ての重要な炉心プロセス変数の予測を行なう。
これらの軌跡は、炉心ガバナ６８のスクリーン８２上に示されることになる。
これらの軌跡の原点（ｔ＝０）は、以下で説明した通りのＢＯＤＥ補給システムに作用するＢＯＤＥ制御装置６４について、およびＨバンクについての現在の設定を示している。

包括的な視点から見ると、所与の時間における図１の実施形態における原子炉２の状態は、本質的に以下の３つの測定可能なおよび独立したプロセス変数または状態変数によって特徴付けされ得る。
１．核分裂出力（ＰＲ）
２．原子炉冷却材温度（ＡＣＴ）、またはＷＷＥＲの場合には、生蒸気圧力（ｐ）
３．炉心出力密度のアキシャルオフセット（ＡＯ）

詳細には、核分裂出力は、核分裂反応によって生成される熱出力である。
それは、一方では測定された冷却材温度上昇（原子炉冷却材の出口および入口温度間の差）および公知の冷却材流量から、そして他方では炉内中性子検出器によって測定される中性子束から導出可能である。

平均原子炉冷却材温度ＡＣＴは、原子炉冷却材の入口および出口温度の平均を表わす。
代替的には、ＷＷＥＲの場合、一次冷却材温度の代りに二次冷却回路１６内の生蒸気圧力ｐを独立プロセス変数として取り上げることができる。

アキシャルオフセットＡＯは、上部炉心半分と下部炉心半分の間の核分裂出力の正規化された差異を表わす。
それは、通常、炉内中性子束計装によって測定される。
炉内計装の無い原子力発電所は、この目的で炉外束測定を使用する。

したがって、所与の時間における原子炉２の状態を、図３で視覚化されているような前記状態変数によってスパンされた３次元状態空間内の点または状態ベクトルによって表わすことができる。
このとき、状態変数の時間的推移は、前記状態空間内の軌跡Ｔｓによって表わされる。
換言すると、軌跡Ｔｓは、時間の関数として状態空間内の原子炉状態の座標を含む。

これら３つの独立した状態変数以外に、原子炉毒作用および炉心内のその空間分布が、反応度に対して有意な寄与を提供する炉心の状態の重要な特性である。
核分裂出力、冷却材温度およびアキシャルオフセットとは異なり、原子炉毒作用は、アクチュエータを用いて直接制御できない。
これはむしろ、主として過去２４時間由来の炉心運転履歴の結果としてもたらされ、明らかに原子炉の制御中考慮に入れるべきものである。

同様にして、図１の実施形態において原子炉２の状態を独立して管理する３つの主要な起動変数が存在する。
１．Ｐバンクの位置
２．一次冷却回路１０、そして場合によってはグレイ棒内のホウ素濃度
３．原子炉が可動Ｈバンクを有する場合、Ｈバンクの位置

一部の特別なフランス型原子炉は、さらに、通常の（ブラック）棒よりも吸収度の低いいわゆるグレイ棒を有する。
グレイ棒には、可能な位置が２つしかない。
すなわち、完全に引抜かれた位置または完全に挿入された位置である。
グレイ棒の利点は、軸方向出力密度分布を変形せず、したがってボレート化／希釈などの炉心内のアキシャルオフセット（ＯＡ）を変更しないことにある。
このような理由から、グレイ棒は、ボレート化／希釈と共に単一の起動変数とみなすことができる。

ホウ素濃度は、ＢＯＤＥ補給ポンプ５２、５４および制御バルブ５６、５８を上述のようにＢＯＤＥ制御装置６４を介して起動させることによって調整可能である。

したがって、起動変数の現在値は、図４で視覚化されているように、起動変数の３次元空間内の点またはベクトルによって表わされ得る。
状態空間についての以上の説明と同様に、このいわゆる起動空間内の軌跡Ｔａは、起動変数の時間的推移を表わしている。

いずれの所与の時点においても、起動変数は原則として、物理的に実現可能な限度内で任意の値に設定可能である。
したがって、起動空間内の軌跡Ｔａは、自由に選択可能であり、起動変数の実用的な最大値および最小値によって、ならびに原子炉の保護および制限システムが与える限界によってのみ制限される。

しかしながら、状態空間内の軌跡Ｔｓは、現在の原子炉毒作用および二次冷却回路による現在の熱除去と共に起動変数の軌跡Ｔａによって左右される。
したがって、軌跡Ｔｓは、内在する原子炉の物理特性に起因して、「自動的に自己調整する」。
管理物理特性の固有の非線形性のため、依存性は複雑なものであり得、これには主に複合毒作用効果に起因して遅延応答が関与し得る。
この意味合いにおいて、起動空間を「一次」空間とみなし、状態空間を従属する「二次」空間とみなしてよい。

それでもなお、単純で周知の反応度収支等式（恒常な原子炉出力Σρ＝０である場合）により、例えば来る２４時間といった一定のタイムスパンについて起動変数の所与の軌跡Ｔａに基づいて原子炉の状態変数の軌跡Ｔｓを予測することが可能である。
換言すると、起動空間内の所与の軌跡Ｔａについて、状態空間内の応分の軌跡Ｔｓは、原則として、反応度収支等式に基づく単純な予測装置モジュールを用いて予測可能である。

平均冷却材温度（ＡＣＴ）の周知の制御以外に、原子炉制御タスクには、所与の限界内のみでの運転を可能にする運転安全性、燃焼の均一化、燃焼補償、毒作用の補償、出力密度分布の均一化、柔軟な電力生産の支援、運転経済性などの他の多くの側面が含まれている。
したがって、制御目標の数は、通常、起動変数の数よりもはるかに多い。
この側面のため、すでに言及した制御すべきシステムの非線形性および複雑な動態、ならびに、即発中性子についての１０μｓから始まって遅発中性子および熱伝達プロセスについての数秒超、毒作用効果についての数十時間超、そして、燃料減損についての数年に至るまでの異なる作用時間スケールの莫大な範囲の広さと合わさって、従来の制御理論は、不適切なものになっている。
多重入出力（ＭＩＭＯ）などの特殊な方法でさえ、不適切に思われる。

上述の問題を克服するため、本発明は、数値変分法を使用する、以上で説明した軌跡表現に基づく全く異なるアプローチを提案するものである。

要約すると、本発明に係る管理スキームは、好ましくは所与の原子炉出力計画の下で、予め設定された将来の大きな時間間隔（２４時間の規模）について多数の起動変数の多数の無作為に変動させられた起こりうる軌跡Ｔａを反復的に考慮するステップを含んでおり、ここで、起動変数の各軌跡Ｔａには、起動変数、プロセス変数および／またはこれらから導出された変数の、予め設定された条件または値によって特徴付けされる多数の事象または不利な炉心状態についての重み付けまたはペナルティ値を含む値表に基づいて、性能指数Σが割当てられる。
その後、応分のアクチュエータを設定するために使用される起動変数の実際の軌跡Ｔａは、性能指数Σが局所極値を有するような形で選択される。

本発明に係る炉心ガバナ６８のこの一般的作動原理は、図１に示されている原子炉２の当該実施形態について図２で概略的に視覚化されており、図１０は、炉心ガバナ６８についての好ましいアルゴリズムを示している。

炉心ガバナ６８は、［状態変数の所与の現在値セット、（負荷ガバナ６６から受信した）来る２４時間についての電力の所与の負荷計画について、および（炉心ナビゲータによって提案される）Ｈバンクの軌跡について］、例えば来る２４時間といった所与の将来のタイムフレームについて状態変数、そして場合によってはそれらから導出される他の変数の軌跡Ｔｓを予測する炉心予測装置と呼ばれるモジュールを含む。
詳細には、すでに言及された状態変数、熱原子炉出力ＰＲ、アキシャルオフセットＡＯ、および原子炉冷却材温度ＡＣＴ以外に、他の変数の時間的推移を導出し予測することができる。
反応度収支等式を用いて、Ｐバンクの軌跡を予測することができる。
さらに、炉心予測装置は、炉心の反応度に有意な影響を及ぼす炉心６内部のキセノン（Ｘｅ）およびサマリウム（Ｓｍ）などの原子炉毒物質の積算濃度および空間分布を予測する。
要約すると、炉心予測装置は、所与の負荷計画、計算された毒作用、そして可動Ｈバンク（存在する場合）およびグレイ棒（存在する場合）の位置およびホウ素濃度などの緩徐動作主体についての提案された軌跡Ｔａからの状態変数およびＰバンク位置の対応する軌跡Ｔｓを導出する。

さらに、炉心ガバナ６８は、炉心ナビゲータと呼ばれるモジュールを含む。
炉心ナビゲータは、一定程度の自由度を有するホウ素濃度、Ｈバンク（存在する場合）およびグレイ棒（存在する場合）という緩徐動作主体についての提案およびその最適化を担当する。

炉心ナビゲータは、原子炉毒作用を考慮に入れて、起動変数の現在の軌跡Ｔａおよび状態変数の応分の軌跡Ｔｓを考慮して、単一の性能指数Σを割当てることによって、現在の軌跡を査定することができる。
この割当ては、図６に一例が示されている値表に基づいている。
値表は、起動変数、プロセス変数、および／またはこれらから導出された変数の予め設定された条件または値により特徴付けされる多数の事象または不利な炉心状態についての重み付けまたはペナルティ値を含んでいる。
これに関連して、「条件」または「値」なる用語は、広義で理解されるべきものである。

明解さおよび描写性を理由として、１つのスカラー量であるそれぞれのペナルティ値は、貨幣価値として、例えばユーロ単位の価格として示され得る（図６）。
詳細には、バンクの運動の数は最小化されるものとし、こうして各バンクステップに一定の値を割当てることができる。
同様にして、ＢＯＤＥシステムによるボレート化および希釈作用には、ペナルティ値が割当てられる。
より一層重要なこととして、状態変数および／または起動変数から導出可能な一定のパラメータまたは変数は、運転安定性の理由から所与の帯域（原子炉保護および原子炉制限システムによって与えられた最小値および／または最大値）を超えてはならない。
したがって、所与の境界または（予備）限界（動的に調整可能なもの）の超過および類似の事象には、比較的高いペナルティ値が割当てられる。

図７および８は、このような予め設定された条件のいくつかの例を示す。
図７は、大きいＡＯについてのペナルティを示す。
すなわち、ペナルティ値は、ＡＯの予め設定された関数である。
頂点は、オペレータステーションのキーボードおよび／またはマウスを用いて予め設定され得る。
図８は、大きい出力密度についてのペナルティを示す。
ここでは頂点は手動で予め設定され得、または原子炉制限システムから連続的に受信され得る。性能指数についての単一の寄与は同様に、原子炉制限システムによって生成される動的限界にも左右され得る。このような動的限界についての軌跡は、炉心予測装置によって計算され得る。

当然のことながら、図６に示されているペナルティ値の所与のリストおよび図７および８に示されているペナルティ関数は、単なる一例にすぎない。概して、リストは拡張および／または補正可能である。この場合、性能指数アルゴリズムも同様に適応させられるべきである。性能指数Σ（全体的価格）は好ましくは、個別のペナルティ値（個別価格）の単純な合計であることから、この適応は些細なことである。

炉心ナビゲータは、ランダム発生器を用いて周期的に、ボレート化／希釈についての軌跡、およびＨバンク（存在する場合）およびグレイ棒（存在する場合）の動きを変動させる。
ボレート化および希釈ならびに棒の運動は、離散的事象である（図９参照）ことから、それらの修正は、これらの動作ならびに対応する振幅の無作為変更の「オン」および「オフ」フロントの時間シフトのように見える可能性がある。
炉心ナビゲータは、同様に、一部の動作を消去するかまたは新たな動作を創出し、既存の動作を一定のタイムスパンで分離されたより小さい振幅または持続時間の２つの動作に分割するか、または２つの隣接する動作を単一の動作に合併する試みを行なうこともできる。
各々のこのような反復サイクルにおいて、炉心予測装置は、Ｐバンクの軌跡を相応して補正し、原子炉状態変数についての軌跡Ｔｓを計算する。
ここでもまた、各々の反復ステップにおいて、炉心ナビゲータは、このように修正された軌跡を計算し、この軌跡に対して単一の性能指数Σを割当てる。
炉心ナビゲータは、修正された軌跡が現在のものよりも優れた性能指数Σを提供する場合に、修正を受諾する。
この場合、修正された軌跡は、現在のものに置換する。
そうでなければ、修正は破棄され、現在の軌跡の新しい無作為変動が試行される。

一般的に言うと、ボレート化／希釈についての実際の軌跡、およびＨバンク（存在する場合）およびグレイ棒（存在する場合）の動きは、性能指数Σが局所極値を有するような形で選択される。
この最適化アルゴリズムは、常時バックグラウンドで周期的に実行され毎秒数千回もの反復を行ない、オペレータには、炉心ガバナ６８のスクリーン８２上で、連続的にゆっくりと変化する画像として現われる（図９）。

バックグラウンドで可能なかぎり最高の速度で実行する上述の高速最適化プロセスと並列で、いわゆる実時間プロセスが実行される。
実時間プロセスは、好ましくは２００ｍｓの時間ステップで実行され、図９に係るオペレータスクリーン８２上の全ての軌跡を左へシフトさせる。
２００ｍｓ毎に全ての軌跡は左へ２００ｍｓだけシフトされ、水平軸上の時間スケールは不動の状態にとどまる。
実行時間は、オペレータには左への全ての軌跡の低速連続クリープとして現われる。
この低速クリープを通して、炉心ナビゲータによって提案されるボレート化／希釈などの動作ならびにＨバンク（存在する場合）およびグレイ棒（存在する場合）の動きは、時として、「現時点」ラインに達することになる。
この時点で、対応する動作が実行のために提案されることになる。
実行は、手動式、半自動式（許可ボタン）または全自動式に行なわれることができる。
同時に、全ての軌跡の出発点（ｔ＝０）は、原子炉出力ＰＲ、冷却材温度ＡＣＴおよびアキシャルオフセットＡＯの測定値を用いて更新されることになる。
「現時点」ラインから左側のスクリーンセグメントは、過去からの測定値を示す。
手動または半自動モードで、提案された動作がオペレータによって無視される場合、炉心予測装置は迅速に全ての軌跡を相応して更新することになり、上述のその標準アルゴリズムはそのために十分なものであり、特別なアルゴリズムは全く不要である。
同様にして、負荷計画は、送電網オペレータまたはプラントオペレータによっていつでも更新可能である。
上述のアルゴリズムは、このような再送出に対処するのに十分なものであり、付加的なアルゴリズムは全く不要である。
運転中の値表の起こりうる手動式変更についても同様であり、上述のアルゴリズムはそれ自体でこのような変更に対処するのに十分なものである。

要約すると、さまざまな制約の下で状態Ａから状態Ｂに原子炉を移す制御動作を発見するという多義的でかつ数学的に過剰決定された制御問題（逆問題）を解こうとする代りに、本発明に係る方法は、無作為に修正された制御動作を試行し、最高／最低（定義および符号に応じた）の性能指数が割り当てられたものを選択することによって、応分の順問題を何回も反復する。

換言すると、炉心ガバナ６８は、例えば来る２４時間といった所与の最適化タイムフレームについて、起動変数（例えば棒の運動、ボレート化および希釈動作）の起こりうる軌跡Ｔａを考慮し、各々の起こりうる軌跡Ｔａについて性能指数Σを計算する。
軌跡Ｔａを修正し反復して、炉心ガバナ６８は、最高／最低の性能指数Σを与える最良の軌跡Ｔａを探し、この最良の発見された軌跡Ｔａを、好ましくは対応する状態軌跡Ｔｓと共に表示し、実時間で制御動作を提案する。
したがって、性能指数Σは、組合わされた欠点を自動的に重み付けまたは平衡化する、組合わされた欠点の定量化された全体値としてみなされ得る。
作動原理は、人工頭脳の分野に属し、ストリートナビゲータまたはチェスコンピュータの基礎にあるものと類似している。

システムは、好ましくは、好適な入力および出力デバイスで補完される。
図９に示されている炉心ガバナ６８の見本の出力スクリーンが、炉心予測装置により来る２４時間について予測された状態変数の軌跡Ｔｓ、および炉心ナビゲータによって提案されたアクチュエータについての軌跡Ｔａを表示する。

以下の段落では、システムの初期化も同様に網羅するいくつかの実装の詳細が、図１および２に基づく特定の実施形態について開示されている。

好ましいアルゴリズムは、図１０に示されている。
１．
予測は、好ましくは来る２４時間について時間の関数としての電力を示す図５の負荷計画で始まる。
負荷計画は、送電網オペレータ２２に由来するものであり、負荷ガバナ６６内に記憶される。
負荷計画は、いつでも送電網オペレータによって再送出され得、プラントオペレータによる起こりうる変更の発出後に、運転のために有効となる。
送電網オペレータがプラントオペレータに負荷計画を提供しない場合、プラントオペレータは自ら独自の「最良の推定計画」をコンパイルすることができる。
電気的出力が来る２４時間以内に低下するか増大するかまたは一定にとどまるかのいずれであれ、非常に大まかな計画でも有用である。
２．
負荷計画を用いて、プラントの効率ηおよび熱出力と電力の間の遅延τを考慮に入れて、炉心予測装置モジュール９０は、来る２４時間についての原子炉の熱出力ＰＲ_ｔｈの軌跡を計算する。
予測装置モジュール９２は、出力反応度（ドップラ）ρｐ＝Γ_ＰＲ・ＰＲ_ｔｈについての軌跡を予測する。
３．
予測装置モジュール９６は、核分裂出力および核分裂率についての軌跡を計算する。
ＰＲ_ｆｉｓｓ＝ＰＲ_ｔｈ－Ｐ_{ｄｅｃａｙ}
Ｒ_ｆｉｓｓ＝Ｐ_ｆｉｓｓ／Ｅ_ｆｉｓｓ
４．
核分裂率の軌跡を用いて、モジュール９８は、ＸｅおよびＳｍ濃度およびそれらの反応度ρ_{ｐｏｉｓｏｎ}についての軌跡を計算する。
５．
出力反応度ρｐおよび毒物質反応度ρ_{ｐｏｉｓｏｎ}についての軌跡を用いて、公知の特性を有するＡＣＴ制御装置４８によって制御されるＡＣＴを考慮に入れて、予測装置モジュール９４は、Ｐバンクについての軌跡を計算する。
この計算は、反応度収支等式Σρ＝０に基づいている。
６．
Ｐバンクのこのように計算された軌跡のいくつかの部分が、棒挿入制限システムにより許可されて（将来の時間周期については限界についての予測された軌跡が使用されなければならない）、帯域を離れた場合、モジュール１０４は、反応度収支等式Σρ＝０を用いて対応する時間間隔についてグレイ棒の挿入／引抜きによってかまたはボレート化／希釈によって棒の運動を置換する。
１００％までの高速出力上昇能力が望まれる場合、それは、引抜かれるＰバンクの制限によって保証され得る。
この場合、特殊なモジュールが、現在の出力から１００％までの高速出力上昇に必要とされるＰバンクの反応度についての軌跡を計算し、相応してＰバンクの引抜きを制限する。
１００％の出力上昇保証は、予測間隔全体にわたって、またはその一部について活動化され得る。
７．
原子炉熱出力ＰＲ_ｔｈ、原子炉核分裂出力ＰＲ_ｆｉｓｓならびにＰおよびＨバンクの位置についての軌跡に基づいて、炉心の上半分および下半分における熱出力Ｐ_Ｕ、Ｐ_Ｄ、アキシャルオフセットＡＯおよび原子炉毒物質Ｘｅ_Ｕ、Ｘｅ_Ｄ、Ｓｍ_Ｕ、Ｓｍ_Ｄの２点分布についての第一近似軌跡は、モジュール１００により計算されることになる。
この第一近似は、軌跡の初期化のためにだけ必要とされる。
さらなる周期的運転中、この中間計算は省略されることになる。
８．
モジュール１００によって計算されたアキシャルオフセットＡＯの軌跡についての第一近似を用いて、炉心ナビゲータのモジュール１０２は、Ｈバンクの軌跡についての第１の提案を行なう。
この第一近似は、軌跡の初期化のためにだけ必要とされる。
さらなる周期的運転中、この中間計算は省略されることになる。
９．
モジュール１０４は、モジュール１０２により計算されたＨバンク軌跡の第一近似を考慮に入れて、ボレート化、希釈およびグレイ棒の挿入についての軌跡を更新する。
この更新は、軌跡の初期化のためにだけ必要とされる。
さらなる周期的運転中、この中間計算は省略されることになる。
１０．
原子炉熱出力ＰＲ_ｔｈ、原子炉核分裂出力ＰＲ_ｆｉｓｓならびにＰおよびＨバンクの位置についての軌跡に基づいて、炉心の上半分および下半分における熱出力Ｐ_Ｕ、Ｐ_Ｄ、アキシャルオフセットＡＯおよび原子炉毒物質Ｘｅ_Ｕ、Ｘｅ_Ｄ、Ｓｍ_Ｕ、Ｓｍ_Ｄの２点分布についての軌跡は、モジュール１０６により計算されることになる。

今や全ての軌跡が初期化され、修正プロセスを開始することができる。
第１の修正ルーチンは、実行時間に対する更新である。
１．
好ましくは２００ｍｓである各更新サイクルにおいて、全ての軌跡が２００ｍｓだけ左へシフトされることになる。
測定可能なプロセス変数および起動変数の原点は、現在の測定値にしたがって更新されることになる。
２．
モジュール９０、９２、９６、９８、９４および１０６は、現在測定されているプロセス変数およびＰバンク位置に応じて対応する軌跡を更新するために呼出されることになる。

第２の修正ルーチンは、Ｈバンク、ボレート化、希釈およびグレイ棒についての軌跡の最適化を目標としており、恒久的にバックグラウンドで実行される。
それは、可能なかぎり高速で実行し、好ましくは毎秒１００回の最適化サイクルを行なう高速非実時間周期的ルーチンである。
１．
炉心ナビゲータのモジュール１０８は、Ｈバンク、ボレート化、希釈およびグレイ棒の軌跡の小さい無作為修正を行なう。
ボレート化および希釈ならびに棒の運動は離散的事象である（図９参照）ことから、それらの修正は、これらの動作ならびに対応する振幅の無作為変更の「オン」および「オフ」フロントの無作為時間シフトのように見える可能性がある。
炉心ナビゲータは、同様に、一部の動作を消去するかまたは新たな動作を創出し、既存の動作を一定のタイムスパンで分離されたより小さい振幅または持続時間の２つの動作に分割するか、または２つの隣接する動作を単一の動作に合併する試みを行なうこともできる。
２．
モジュール９４は、反応度収支等式Σρ＝０を用いてＰバンクの軌跡を更新し、モジュール１０６は、Ｐ_Ｕ、Ｐ_Ｄ、ＡＯおよびＸｅ_Ｕ、Ｘｅ_Ｄ、Ｓｍ_Ｕ、Ｓｍ_Ｄについての軌跡を更新する。
３．
モジュール１１０は、図６の値表を用いて現在のおよび修正された軌跡セットについての性能指数を計算する。
修正された軌跡が現在のものより優れた性能指数を生み出す場合、現在の軌跡は、修正された軌跡により置換されることになり、そうでない場合には、修正は破棄される。
その後に新しい修正の試みが開始されることになる。

再送出動作は、負荷計画の離散的修正を導く。
モジュール９０は、電力についての現在の有効な軌跡をゆっくりと連続的に新しいものに再整形することができる。
このような変換には数分間かかる可能性があり、現在の軌跡から開始してそれを負荷ガバナの新規要求に合わせて連続的に再整形する。
この低速変換により、追加の特別なアルゴリズムを導入することなく、両方の周期的実行ルーチンが全ての軌跡を再形成することが可能になる。

同様に、軌跡を初期化するための代替的可能性も存在する。
炉心ガバナ６８の始動またはリセットの直後に、来る２４時間にわたる全ての軌跡の定常状態が仮定されることになり（現在測定されている状態変数および動作変数に応じた恒常な軌跡）、両方の周期的ルーチンが開始されることになる。
周期的ルーチンを実行する間のモジュール９０による電力についての有効な軌跡の低速再整形は、全ての軌跡を再整形することになる。
この場合、Ｈバンク、ボレート化、希釈およびグレイ棒についての軌跡は明示的に初期化されることはなく、モジュール１０８の億兆回の無作為動作の中で創出されることになる。
この場合、モジュール１００、１０２および１０４は省略され得る。

全ての軌跡の実行時間および後続する左へのクリープ（図９）に起因して、炉心ナビゲータによって提案されるボレート化／希釈ならびにＨバンク（存在する場合）およびグレイ棒（存在する場合）の運動などの動作は、時として「現時点」ラインに達することになる。
その時点で、対応する動作が実行のために提案されることになる。
実行は手動式、半自動式（許可ボタン）または全自動式で行なうことができる。

手動または半自動モードにおいて、提案された動作がオペレータにより破棄される場合、炉心予測装置は迅速に全ての軌跡を相応して更新することになり、上述の標準アルゴリズムはそのために十分なものであり、特別なアルゴリズムは全く不要である。

状態変数、動作変数および全体的性能指数の全ての軌跡を含めた、来る２４時間についてのこのように計算された全体的計画は、炉心ガバナ６８のスクリーン８２上に絶えず表示される。
画像は、ゆっくりと左へクリープし、連続的にそれ自体を改善する。性能指数は同様に表示されることになり、この連続的改善を例示する。

炉心ガバナは、その主要な機能以外に、来る２４時間についての柔軟な運転に関する発電所の可能性を調査することを可能にする。
この目的で、マウスによりプロットの頂点を移動させて、負荷ガバナのスクリーン上の負荷計画の変更を試みることができる。
コンピュータの性能に応じて数秒または数分以内に、全ての状態変数を描写する全ての軌跡に対するこの試験的変更の影響を観察することが可能である。

２ 加圧水型原子炉（ＰＷＲ）
４ 原子力発電所
６ 炉心
８ 原子炉圧力容器
１０ 一次冷却回路
１２ 原子炉冷却材ポンプ
１４ 蒸気発生器
１６ 二次冷却回路
１８ 蒸気タービン
２０ 発電機
２２ 送電網
２４ 凝縮器
２６ 給水ポンプ
２８ 給水タンク
３０ タービンバルブ
３２ 蒸気補給ライン
３４ バイパスライン
３６ バイパスバルブ
４０ タービン制御装置
４２ バイパス制御装置
４４ 制御棒
４６ 中性子束制御装置
４８ 平均冷却温度（ＰＷＲ、ＤＷＲ）または生蒸気圧力（ＷＷＥＲ）制御装置
５０ 炉外計装
５２ 補給ポンプ
５４ 補給ポンプ
５６ 制御バルブ
５８ 制御バルブ
６４ ＢＯＤＥ制御装置
６６ 負荷ガバナ
６８ 炉心ガバナ
７０ 炉内計装
８０ 負荷ガバナのスクリーン
８２ 炉心ガバナのスクリーン
９０ 原子炉熱出力用予測装置モジュール
９２ 出力（ドップラ）反応度用予測装置モジュール
９４ 平均冷却材温度およびＰバンクの位置用の予測装置モジュール
９６ 核分裂出力用予測装置モジュール
９８ 原子炉毒作用（Ｘｅ、Ｓｍ）用予測装置モジュール
１００ 炉心内の出力密度および毒物質密度用の第１の予測装置モジュール
１０２ Ｈバンクの位置の第１の提案用のナビゲータモジュール
１０４ ＢＯＤＥおよびグレイ棒用のナビゲータモジュール
１０６ 炉心内の出力密度および毒物質密度用の第２の予測装置モジュール
１０８ Ｈバンク、ＢＯＤＥおよびグレイ棒についての軌跡をランダム化するナビゲータモジュール
１１０ 性能指数を計算するナビゲータモジュール

ＡＣＴ 平均冷却材温度
ＡＯ アキシャルオフセット
ＢＯＤＥ ホウ酸／脱塩水（ＢＯｒｉｃ ａｃｉｄ／ＤＥｍｉｎｅｒａｌｉｚｅｄ ｗａｔｅｒ）
ＤＷＲ Ｄｒｕｃｋ Ｗａｓｓｅｒ Ｒｅａｋｔｏｒ（加圧水型原子炉）（独語）
ＰＲ 原子炉出力
ＰＷＲ 加圧水型原子炉
Ｔａ 起動変数の軌跡
Ｔｓ 状態変数の軌跡
ＷＷＥＲ 水－水型原子炉（旧ソ連型ＰＷＲ）
Σ 性能指数

Claims (8)

1. 加圧水型原子炉（２）を管理する方法において、
   前記加圧水型原子炉（２）が、炉心（６）と、原子炉冷却材を保持する前記炉心（６）用冷却回路（１０）とを含み、
   前記原子炉の状態が、多数の測定可能な状態変数によって特徴付けされ、前記炉心の反応度が、多数の起動変数によって制御されかつ原子炉毒作用により影響され、
   起動変数の所与の時間依存軌跡（Ｔａ）について、前記状態変数についての応分の軌跡（Ｔｓ）が、前記状態変数の測定された現在の値、計算された毒作用値および反応度収支等式に基づいて予測される、
   方法であって、
   時間の関数として電力を示す負荷計画を受信するステップと、
   将来の時間間隔について前記起動変数の多数の無作為に変動させられた起こりうる軌跡（Ｔａ）を反復的に考慮するステップであって、起動変数軌跡（Ｔａ）の前記変動が前記原子炉出力についての所与の負荷計画の制約の下で実現されるステップと、
   を含み、
   起動変数の各軌跡（Ｔａ）には、前記起動変数、プロセス変数、および／またはこれらから導出された変数の予め設定された条件または値によって特徴付けされる多数の事象または不利な炉心状態についての重み付けまたはペナルティ値を含む値表に基づいて、性能指数（Σ）が割り当てられ、
   前記起動変数軌跡（Ｔａ）が、性能指数（Σ）が局所的極値を有するような形で選択され、
   対応するアクチュエータが相応して移動させられる、
   方法。
2. 前記プロセス変数が、熱原子炉出力（ＰＲ）、原子炉冷却材温度（ＡＣＴ）、生蒸気圧力（Ｐ）および／または炉心出力密度のアキシャルオフセット（ＡＯ）、または等価の数量のうちの単数または複数のものを含む、請求項１に記載の方法。
3. 前記加圧水型原子炉（２）がさらに、好ましくは制御バンクとして統合された、前記炉心（６）内に挿入可能な多数の制御棒（４４）と、前記原子炉冷却材のホウ素濃度を設定するためのボレート化および希釈システムとを含み、前記起動変数が制御棒の位置または運動および／またはボレート化動作および／または希釈動作に特徴的なものである、請求項１または２に記載の方法。
4. 前記起動変数軌跡（Ｔａ）の前記変動が、前記炉心の所与の毒作用の制約の下でさらに実現される、請求項１から３のいずれか一つに記載の方法。
5. 前記状態変数の前記測定された現在値が実時間で更新される、請求項１から４のいずれか一つに記載の方法。
6. 前記計算アルゴリズムが乱数発生器を使用する、請求項１から５のいずれか一つに記載の方法。
7. 前記考慮対象の将来の時間間隔が、１時間よりもはるかに長く、典型的には２４時間前後の規模を有する、請求項１から６のいずれか一つに記載の方法。
8. 請求項１から７のいずれか一つに記載の方法が実装されるコンピュータ化された管理モジュールを伴う、原子炉用の管理システムにおいて、前記計算された最適な起動変数軌跡（Ｔａ）が、対応するアクチュエータを制御するために使用される、管理システム。

Images