(57)【要約】      （修正有）【課題】セル摩擦を、原子炉サイクル炉心の設計、最適化、ライセンス付与、および監視の一部として管理する方法を提供する。【解決手段】チャネル面高速フルエンス・チャネル面被制御運転パラメータを、すべてのチャネルについて判定する５１０。総バウ値を、チャネルごとに、チャネル面高速フルエンスおよび／またはチャネル面制御パラメータに基づいて計算する５２０。チャネル壁圧力降下パラメータを判定し、総バルジ値を、チャネルごとに、チャネル面高速フルエンスパラメータおよびチャネル壁圧力降下パラメータを使用して計算する５３０。セルの指定されたチャネル軸方向高度での総変形を、総バウ値および総バルジ値に基づいて判定する５４０。制御ブレード軸方向摩擦力値を、各軸方向高度で、チャネル剛性値およびチャネル−制御ブレード摩擦係数値と共に総変形に基づいて計算する５５０。最大摩擦値を、セルのセル摩擦測定基準として選択する。【選択図】図５
2008203162.tif 000002
【特許請求の範囲】
【請求項１】
原子炉の制御セルのセル摩擦測定基準を判定する方法であって、
前記制御セル内のチャネルごとに、チャネル面高速フルエンスパラメータおよびチャネル面被制御運転パラメータのうちの一方または両方を判定すること（５１０）と、
前記制御セル内のチャネル面ごとに、複数のチャネル軸方向高度のそれぞれで、総バウ値を計算すること（５２０）と、
前記制御セル内のチャネル面ごとに、各チャネル軸方向高度で、総バルジ値を計算すること（５３０）と、
前記制御セルの各チャネル軸方向高度で、前記総バウ値および前記総バルジ値に基づいて、総変形を判定すること（５４０）と、
前記軸方向高度のそれぞれで、前記総変形に基づいて、セル軸方向摩擦力値を計算すること（５５０）と、
前記計算されたセル摩擦力値の最大値を、前記制御セルの前記セル摩擦測定基準として選択することと
を含む方法。
【請求項２】
チャネルゆがみと前記セルの制御ブレードの運転に対する結果の制御ブレード軸方向摩擦負荷との重大度を査定するために、複数の重大度閾値に対して前記セル摩擦測定基準を比較することをさらに含むことを特徴とする請求項１記載の方法。
【請求項３】
総バウ値を計算することが、炉心内の各セルの計算された高速フルエンス勾配によって誘導されたバウ値（Ｓ２１２）および計算されたシャドウ腐食によって誘導されたバウ値（Ｓ２２２）のうちの一方または両方に基づいて炉心平均セル平均バウ値を判定すること（Ｓ２３３）を含むことを特徴とする請求項１記載の方法。
【請求項４】
前記セル平均バウ値を判定することが、
前記セル内のチャネルごとに高速フルエンス勾配によって誘導されたバウ値を計算すること（Ｓ２１０）と、
総バウ値を得るために、前記高速フルエンス勾配によって誘導されたバウを最初の製造されたバウ値に加算すること（Ｓ２１１）と、
高速フルエンス勾配バウ不確実性値を計算すること（Ｓ２１２）と、
前記セル平均バウ値とセル平均バウの不確実性とを判定するために、前記総バウ値をその計算された不確実性と組み合わせること（Ｓ２１３）と
を含むことを特徴とする請求項３記載の方法。
【請求項５】
セル平均バウ値を判定することが、
前記セル内のチャネルごとにシャドウ腐食によって誘導されたバウ値を計算すること（Ｓ２２０）と、
総バウ値を得るために、前記シャドウ腐食によって誘導されたバウを前記セルの最初の製造されたバウ値に加算すること（Ｓ２２１）と、
シャドウ腐食バウ不確実性値を計算すること（Ｓ２２２）と、
前記セル平均バウ値とセル平均バウの不確実性とを判定するために、前記総バウ値をその計算された不確実性と組み合わせること（Ｓ２２３）とを含むことを特徴とする請求項３記載の方法。
【請求項６】
セル平均バウ値を判定することが、
前記セル内のチャネルごとに高速フルエンス勾配によって誘導されたバウ値およびシャドウ腐食によって誘導されたバウ値を計算すること（Ｓ２３０）と、
総バウ値を得るために、前記高速フルエンス勾配によって誘導されたバウ値および前記シャドウ腐食によって誘導されたバウを前記セルの最初の製造されたバウ値に加算すること（Ｓ２３１）と、
高速フルエンス勾配バウ不確実性およびシャドウ腐食バウ不確実性値を計算すること（Ｓ２３２）と、
前記セル平均バウ値とセル平均バウの不確実性とを判定するために、前記総バウ値をその計算された不確実性と組み合わせること（Ｓ２３３）とを含むことを特徴とする請求項３記載の方法。
【請求項７】
前記制御セルの前記総バルジ値を計算することが、
前記制御セルの制御ブレードのブレードウィングに面しているチャネルの面ごとに、各軸方向高度で弾性バルジ値およびクリープ値を計算すること（Ｓ３００）と、
ブレードウィングに面するチャネル面ごとの各軸方向位置での総変形値を有するために、ブレードウィングに面しているチャネルの各面での軸方向高度ごとに前記総バウ値および前記総バルジ値を加算することを含めるために総変形を判定すること（Ｓ３０１）とを含むことを特徴とする請求項１記載の方法。
【請求項８】
前記軸方向高度のそれぞれでセル摩擦力値を計算することが、
ブレードウィングに面しているチャネルの各面での各軸方向高度で、公称摩擦力値（Ｓ３０３）および摩擦力値の不確実性（Ｓ３０５ｂ）を判定することと、
前記所与の軸方向高度での前記セルの公称上界摩擦力値（Ｓ３０５ｃ）および統計的上界摩擦力値（Ｓ３０６）を判定するために、すべての面の前記公称摩擦力値および前記不確実性を組み合わせることを含むことを特徴とする請求項７記載の方法。と
【請求項９】
公称摩擦力を判定することが、
各軸方向高度で、その軸方向高度での前記総変形に基づいて、所与のチャネル面とそれに面する制御ブレードウィングとの間の公称干渉値を計算すること（Ｓ３０２）と、
各面の前記計算された干渉およびチャネル−制御ブレード摩擦係数に基づいて、チャネル剛性値を使用して、前記計算された公称干渉値を公称摩擦力値に変換すること（Ｓ３０３）とを含むことを特徴とする請求項８記載の方法。
【請求項１０】
上界摩擦力を判定することが、
各軸方向高度で、その高度での前記総変形に基づいて、所与のチャネル面とそれに面する制御ブレードウィングとの間の上界干渉値を計算することと、
各面の前記計算された干渉およびチャネル−制御ブレード摩擦係数に基づいて、チャネル剛性値を使用して、前記計算された上界干渉値を上界力値に変換することとを含むことを特徴とする請求項８記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、全般的には、原子炉の制御セル（ｃｏｎｔｒｏｌ  ｃｅｌｌ）内の制御ブレードの軸方向移動に関するセル摩擦測定基準（ｃｅｌｌ  ｆｒｉｃｔｉｏｎ  ｍｅｔｒｉｃ）の判定に関する。
【背景技術】
【０００２】
図１は、燃料チャネルを示すための例の燃料バンドルの上面図である。沸騰水型原子炉（ＢＷＲ）では、燃料バンドル１１０が、通常、比較的薄く長方形の燃料チャネル１２０に入れられる。グリッド内で、図１に、燃料バンドル１１０の２次元（２Ｄ）レイアウトが示されており、この２Ｄレイアウトは、この例では、燃料チャネル１２０に囲まれた棒グリッド位置１１５の１０×１０マトリックスからなる。１０×１０棒グリッド位置の一部は、燃料バンドル１１０内の異なる構成要素を示す、より大きい円を形成するために組み合わされている。しかし、そのような詳細は、この議論には関連しない。というのは、図１が、単に、燃料バンドル１１０を囲む燃料チャネル１２０を示すために提供されているからである。燃料チャネル１２０は、ＢＷＲ炉心内の燃料支持板の上に約１６５インチ延び、約０．１０インチの厚さを有する。
【０００３】
図２に、ＢＷＲの制御ブレードセル２００を見下ろす２Ｄ上面図を示す。燃料バンドル１１０を、燃料チャネル１２０と一緒に、通常は燃料アセンブリ２１０と称する。図２に示されているように、ＢＷＲ炉心内では、４つの燃料アセンブリ２１０が、１つの十字形の制御ブレード２３０によって制御される形で位置決めされる。４つのブレードウィング２３５を有する制御ブレード２３０は、図２に示されているように、４つの燃料アセンブリ２１０の中心を通過する。この形で、各チャネル２２０の２つの面２２５は、個々の制御セル２００内の合計８つの面について、必ずブレードウィング２３５に面する。ＢＷＲ炉心は、多数のそのような制御セルすなわち、制御ブレード２３０の周囲の４つの燃料アセンブリ２１０のグループの繰り返しとして表すことができ、これらの制御セルは、炉心にまたがって多数の位置に配置される。たとえば、ＢＷＲ炉心は、通常、数百個の制御セル２００および数百個の燃料チャネル２２０からなる。
【０００４】
図３に、制御セル内のチャネルゆがみの２つの機構を示す。図示のために、制御セル２００は、図３では間隔をおかれた燃料チャネル２２０を有し、運転中のＢＷＲでの真の寸法を示さない。燃料チャネル２２０は、運転中のＢＷＲ内のさまざまな放射線応答および機械的応答に起因するチャネル変形を受ける。図３には、２つのそのような種類のチャネル変形すなわち、「バウ（ｂｏｗ）」および「バルジ（ｂｕｌｇｅ）」が示されている。図（ａ１）および（ｂ１）は、炉心内の所与の軸方向高度での制御セル２００の２Ｄ上面図であり、図（ａ２）は、大写しの正面図であり、制御ブレードウィング２３５に関するチャネル面位置に対するチャネルバウ機構の影響を示し、図（ｂ２）は、大写しの正面図であり、制御ブレードウィング２３５に関するチャネル面位置に対するチャネルバルジ機構の影響を示す。
【０００５】
バウ機構では、チャネル面２２５が、（ａ２）に示されているように、ブレードウィングに向かってまたはブレードウィング２３５から離れてのいずれかで変形し、軸方向には全体的に正弦波の形状を有するが、この全体的な形状の多少の変形が発生し得る。図３では、図（ａ１）のオフセットした線が、この変形を示し、制御セル２００内のバウに起因するチャネル２２０のオフセットを示し、このオフセットは、単にチャネル２２０の変位である。バルジ機構では、チャネル２２０は、図３の図（ｂ２）に示されているように膨らむ。このバルジは、図（ｂ１）に示されているように、チャネル２２０の４つのすべての面２２５について外向きである。
【発明の開示】
【発明が解決しようとする課題】
【０００６】
チャネル変形は、ＢＷＲの多数の運転パラメータおよび安全パラメータに影響し、したがって、原子炉サイクル炉心の設計、最適化、ライセンス付与、および監視の一部として対処されなければならない。チャネル変形は、チャネル−制御ブレード干渉をもたらす可能性があり、このチャネル−制御ブレード干渉は、ブレード移動中のブレードに対する軸方向摩擦負荷（セル摩擦とも称する）をもたらし、この軸方向摩擦負荷は、制御セル２００内の制御ブレード２３０の動作を妨げる場合がある。
【０００７】
図４に、制御セル内のセル摩擦を示す。チャネル−制御ブレード干渉は、図４では、２つの制御セルについて示されており、図４には、炉心を見下ろしているかのように、制御セルの上面図が示されている。図４では、４つのチャネル面が、別の図で後に参照するために、１つのチャネルについて１〜４の番号を付けられている。制御ブレードおよびチャネルの寸法は、制御ブレードウィングと隣接するチャネル面との間にギャップを作るように指定される。左側のセル２００は、変形が全くない４つのチャネル２２０を示し、したがって、中央の制御ブレード２３０は、製造状態のクリアランスを維持し、したがって、炉心に出入りして（すなわち、ページの平面を通って）自由に移動するのに十分なクリアランスを維持する。その一方で、右側の制御セル２００’には、実質的な変形を有する４つのチャネル２２０’が示されており、この変形は、制御ブレードウィングとチャネルとの間の製造状態のギャップを減らす。チャネルに隣接する２つの面のオフセットした線は、４１０で実際に互いに触れ、制御ブレード２３０のためのギャップが残されていない。そのような状態は、チャネル−制御ブレード干渉をもたらす可能性があり、このチャネル−制御ブレード干渉は、制御ブレード２３０’が制御セル２００’に出入りして移動される時に制御ブレード軸方向摩擦力をもたらす。この状態を、「セル摩擦」と称する。ＢＷＲ内の制御ブレードは、原子炉を安全かつ効率的に運転するのに能動的に使用される。制御ブレード動作に対するすべての妨害が、望ましくない緩和処置につながる可能性がある。そのような処置は、たとえば、最大出力運転戦略にペナルティを与え、失われた発電所利用率および／またはこうむられる置換電力コストをもたらす。したがって、セル摩擦を、原子炉サイクル炉心の設計、最適化、ライセンス付与、および監視の一部として管理し（査定し、緩和し）なければならない。
【課題を解決するための手段】
【０００８】
本発明の例の実施形態は、原子炉の制御セルのセル摩擦測定基準を判定する方法を対象とする。この方法では、チャネル面高速フルエンス（ｆａｓｔ  ｆｌｕｅｎｃｅ）および／またはチャネル面被制御運転パラメータが、制御セルのチャネルについて判定される。総バウ値が、チャネル面高速フルエンスパラメータおよび／またはチャネル面制御パラメータに基づいてチャネルごとに計算される。チャネルごとに、チャネル壁圧力降下パラメータが、判定され、総バルジ値が、チャネル面高速フルエンスパラメータおよびチャネル壁圧力降下パラメータを使用してチャネルごとに計算される。制御セルの指定されたチャネル軸方向高度での総変形が、総バウ値および総バルジ値に基づいて判定される。制御ブレード軸方向摩擦力値が、チャネル剛性値およびチャネル制御ブレード摩擦係数値と共に総変形に基づいて、各軸方向高度で計算される。最大の摩擦力値が、制御セルのセル摩擦測定基準として選択される。
【０００９】
本発明のもう１つの実施形態は、複数のセルを有する原子炉炉心の炉心平均セル平均バウ値を判定する方法を対象とする。この方法では、炉心内のセルごとに、セル平均バウ値が、計算された高速フルエンス勾配によって誘導されたバウ値と計算されたシャドウ腐食（ｓｈａｄｏｗ  ｃｏｒｒｏｓｉｏｎ）によって誘導されたバウ値とのうちの一方または両方に基づいて判定される。判定された値をセルごとに統計的に組み合わせて、炉心平均セル平均バウ値と、炉心の炉心平均セル平均バウの不確実性とを得る。
【００１０】
本発明の例の実施形態は、本明細書で下で示す詳細な説明および添付図面からより十分に理解されるようになるが、添付図面では、類似する要素が、類似する符号によって表され、添付図面は、例示のみのために与えられ、したがって、本発明の例の実施形態について制限的ではない。
【発明を実施するための最良の形態】
【００１１】
本明細書で使用される制御セルは、「ブレード−センタードセル（ｂｌａｄｅ−ｃｅｎｔｅｒｅｄ  ｃｅｌｌ）」とも称するが、燃料バンドルのグループの間で説明できる制御ブレードとして表すことができる。もう１つの例では、セルは、インスツルメント−センタードセル（ｉｎｓｔｒｕｍｅｎｔ−ｃｅｎｔｅｒｅｄ  ｃｅｌｌ）として理解することができ、このインスツルメント−センタードセルは、燃料バンドルのグループの間で説明できる機器チューブとして表すことができる。したがって、いくつかの実例のセルは、「インスツルメント−センタード」セルまたはブレード−センタードセルと見なすことができる。というのは、ＢＷＲにおいて、いくつかの位置が、４つのバンドルによって囲まれたプラント計装チューブを有するからである。
【００１２】
図５は、ＢＷＲ制御セルのセル摩擦測定基準（ＣＦＭ）を判定する全般的な方法を示す処理流れ図である。ＣＦＭは、ＢＷＲの運転サイクルごとの仕様の炉心の設計フェーズ、最適化フェーズ、ライセンス付与フェーズ、および／または監視フェーズのうちの１つまたは複数の間に、制御セルの１つまたは複数のチャネル内のチャネルゆがみおよび結果のセル摩擦の影響を査定し、緩和する手段を炉心設計者に提供する。この例の方法は、炉心の設計、最適化、ライセンス付与、および監視に使用される既存の方法またはソフトウェア実装にリンクし、かつ／またはそれにコーディングすることができる。１例として、この方法を、ＢＷＲ炉心の設計、最適化、ライセンス付与、および監視に使用されるコンピュータプログラムの組の一部であるコンピュータソフトウェアモジュール内で実装することができる。
【００１３】
一般に、制御セルのセル摩擦測定基準を判定する方法には、総変形とも称する総チャネル面変位値を入手するために組み合わせなければならない複数のチャネル変位を計算することが含まれる。次に、計算された総チャネル変位を使用して、制御ブレードウィングと隣接するチャネルとの間の干渉の量を計算する。摩擦負荷は、計算された干渉値と、干渉依存チャネル剛性値と、嵌合するチャネルの材料および制御ブレードの材料の既知のまたは測定された摩擦係数とを使用して計算される。総チャネル変位（変形）値、チャネル−制御ブレード干渉値、および結果の摩擦力は、セルの軸方向高度ごとに計算することができ、その後、セル摩擦力値またはセル摩擦測定基準を判定することができる。例のセル摩擦方法は、たとえばＢＷＲ炉心の設計、最適化、ライセンス付与、および監視に使用されるプログラムの一部として、対話型最適化プロセスに使用されるコンピュータプログラムのモジュールの一部として実装することができる。
【００１４】
方法５００では、設計されまたは評価される炉心の１つまたは複数の制御セルごとに、チャネル壁圧力降下（すなわち、チャネルの内側と外側の圧力の差）、チャネル壁高速フルエンス、および被制御運転パラメータを含むがこれらに限定されない、チャネル変形に影響する運転要因を、ＰＡＮＡＣＥＡ（登録商標）などの炉心シミュレータを使用してそのような放射線応答および機械的応答を定量化することによって判定５１０することができる。チャネル面の変形または変位（すなわち、５２０ではバウ値、５３０ではバルジ値）を、理論的考慮からまたは経験に基づく関係から分析的に導出された既知の数学的表現と共にこれらの計算された運転パラメータを使用して計算することができる。
【００１５】
制御セルの複数の軸方向高度のそれぞれでの総（チャネル面）変形値を、総バウ値および総バルジ値に基づいて判定することができる（５４０）。セル摩擦力値を、総変形および結果のチャネル−制御ブレード干渉に基づいて、軸方向高度のそれぞれで判定することができる。軸方向高度のそれぞれの計算されたセル摩擦力値の最大値を、その制御セルのセル摩擦測定基準として採用する（５５０）。セル摩擦力の計算された軸方向分布は、実際に達成される正味セル摩擦に寄与するので、この軸方向分布に対処するために、代替の処理を適用することができる。
【００１６】
図６〜８に、図５に示された方法のプロセスフロー機能をより詳細に示す。図６では、機能Ｓ１００からＳ２００は、ユーザまたは炉心設計者の計算オプションに応じて、高速フルエンスパラメータおよび／またはチャネル被制御運転パラメータを計算するプロセスを示す。図６に示された高速フルエンスパラメータおよび被制御運転パラメータの計算は、標準炉心シミュレータアプリケーションの通常の部分として実行することができ、あるいは、必要な情報を、標準炉心シミュレータアプリケーションの通常の部分として実行される計算からたやすく抽出可能とすることができる。図６には示されていないが、ユーザまたは炉心設計者は、熱膨張、製造状態での応力緩和によって誘導されたチャネルバウ、チャネルねじれ／回転、または製造状態での冷間加工によって誘導されたチャネルバウなど、他のチャネル変形機構に関する計算を選択することもできる。
【００１７】
例示として、一般化された高速フルエンス蓄積は、式（１）によって例示される。
【００１８】
【数１】
2008203162.tif 000003
【００１９】
式（１）において、ＦＬＵＮＣＥ（ｋ，ｉ，ｊ，ｎ）は、チャネル（ｉ，ｊ）のチャネル面ｎの軸方向高度ｋでのチャネル面中性子フルエンス［中性子／ｃｍ２］である。炉心座標（ｉ，ｊ）は、炉心内のチャネル位置を一意に識別する。ＦＬＵＮＣＥｔ−１（ｋ，ｉ，ｊ，ｎ）は、現在の時間ステップＤＴの始めでのフルエンスである。時間ステップＤＴは、出力を作るためのコアバーン（ｃｏｒｅ  ｂｕｒｎ）を追跡する際の時間増分を表す。ＦＬＵＸＳ（ｋ，ｉ，ｊ，ｎ）は、「高速」中性子の特徴を表すために指定されたエネルギレベルを超えるチャネル面中性子束［中性子／ｃｍ２−ｓｅｃ］である。
【００２０】
ＦＬＵＸＳ（ｋ，ｉ，ｊ，ｎ）計算は、通常の炉心シミュレーションで実行される標準的な計算の結果の単純な産物である。チャネル照射成長は、特定のチャネル材料に関する蓄積された高速フルエンスの既知の関数であるが、集合組織（ｔｅｘｔｕｒｅ）、残留冷間加工、およびチャネル水素含量などであるがこれらに限定されないチャネル材料特性の関数でもある。チャネル材料に関する照射成長の関係を計算されたチャネル高速フルエンスと共に用いて、対向するチャネル面の総照射成長を計算することができる。高速フルエンス勾配によって誘導されたチャネルバウは、チャネル幾何形状パラメータと共に、対向するチャネル面の差成長からたやすく計算することができる。代替案では、高速フルエンス勾配によって誘導されたチャネルバウを、チャネルバウ測定値と、個々のチャネルにまたがる計算された照射線量勾配の使用によるなどの蓄積された高速フルエンス勾配の近似とから導出される経験的関係から計算することができる。
【００２１】
第２の例示として、一般化されたチャネル面被制御運転パラメータは、式（２）によって例示することができる。
【００２２】
【数２】
2008203162.tif 000004
【００２３】
式（２）において、ＥＣＢＥ（ｉ，ｊ）は、被制御運転パラメータであり、ＥＣＢＥｔ−１（ｉ，ｊ）は、現在の時間ステップの始めまでに蓄積された被制御パラメータである。ＬＥＮＧＴＨ（ｉ，ｊ）は、係数ｆによって重みを付けられた、現在の時間増分ＤＴに制御されるチャネル長である。係数ｆは、チャネルの総滞留時間に依存する有効被制御照射線量重み付け係数であり、０．０から１．０までの範囲にわたることができる。ｆの定義は、測定されたチャネル変形に対する予測されたチャネル変形の比較から決定される。この係数は、制御ブレードに対するチャネル照射線量の相対的な重要性に対処することができる。というのは、この重要性が、たとえば制御ブレードハンドルによるより大きい寄与などの軸方向感度依存性を反映すると同時にチャネル運転寿命中に変化する可能性があり、あるいは、ピーク高速フルエンスの軸方向位置での制御に伴って変化する可能性があるからである。この一般化されたチャネル面被制御運転パラメータを用いると、制御ブレードシャドウ腐食によって誘導されたチャネルバウを、被制御運転パラメータおよびたとえば総蓄積チャネル平均照射線量などの他の重要な性能パラメータの関数としてのチャネルバウの経験に基づく関係から計算することができる。
【００２４】
図７では、機能Ｓ２００からＳ３３０は、全体的に、ユーザまたは炉心設計者によって選択された計算オプション（１、２、または３）に応じた、炉心平均セル平均バウ（ＢＯＷＡＶＧ）値を判定するための高速フルエンス勾配によって誘導されたチャネルバウおよび／または制御ブレードシャドウ腐食によって誘導されたチャネルバウの計算を記述する。これらの計算は、バウ値を判定するための理論に基づくまたは経験に基づく関係と共に、適当な運転中性子パラメータの炉心シミュレータ計算に基づいて実行される。不確実性値または炉心シミュレータ運転パラメータの関数としての不確実性関係は、測定されたチャネル変形に対する、理論に基づくまたは分析に基づくモデルを用いて予測されたチャネル変形の比較から決定される。
【００２５】
機能Ｓ２００は、チャネルの最初の製造状態での（製造された）バウを考慮に入れる。ＰＡＮＡＣＥＡ（登録商標）炉心シミュレータなどの既知の炉心シミュレータでは、このパラメータは、１つの包括的な値を割り当てるためのチャネルタイプおよびプラントタイプの依存性を反映する複数の包括的な値に基づく。機能Ｓ２１０、Ｓ２２０、およびＳ２３０では、高速フルエンス勾配によって誘導されたバウおよび／またはシャドウ腐食によって誘導されたバウの計算が、炉心内のチャネルごとに、各軸方向高度で各チャネル面で実行される。高速フルエンス勾配によって誘導されたバウは、対向するチャネル面の差成長歪みから計算される。図４の番号を付けられたチャネル面からわかるように、面の２つのそのような組すなわち、チャネル面１−３および面２−４がある。チャネル成長歪みは、各チャネル面への高速フルエンスによって誘導される。高速フルエンスの関数としてのジルカロイ照射成長の経験的相関を、下で、理論に基づくまたは経験に基づく関係と共に、中性子パラメータの使用を例示するために式の組（３）に示す。
GROW(k,i,j,n)=C1*G*Max(RL,RH)   (3)
ただし、
G=C2+C3*(FLUNCE(k,i,j,n))C4、FLUNCE(k,i,j,n)≦C5 n/cm２の場合
=C6+C7*(FLUNCE(k,i,j,n)+C15*(FLUNCE(k,i,j,n))C16
+C8*(FLUNCE(k,i,j,n))C9、
それ以外の場合
RL=C10+C11T+C12T2
RH=C13+C14T
上の式の組（３）において、ＧＲＯＷ（ｋ，ｉ，ｊ，ｎ）は、無次元成長歪みであり、Ｔは、照射温度であり、係数Ｃ１からＣ１６は、理論的基礎および経験的基礎を有する。無次元成長歪を使用することによって、高速フルエンス勾配によって誘導されたバウを、単純な数学的関係を使用して計算することができる。これらの関係は、原子炉技術の当業者に既知であり、したがって、説明を簡潔にするために省略する。代替案では、高速フルエンス勾配によって誘導されたチャネルバウを、チャネルバウ測定値と、個々のチャネルにまたがる計算された照射線量勾配の使用によるなどの蓄積された高速フルエンス勾配の近似とから導出される経験的関係から計算することができる。
【００２６】
シャドウ腐食バウは、ＰＡＮＡＣＥＡ（登録商標）炉心シミュレータでは、上で示したチャネル面被制御運転パラメータＥＣＢＥ（ｉ，ｊ）に依存する一般化された非線形モデルとして特徴を表される。このモデルの例示が、次の多項式関係である。
【００２７】
【数３】
2008203162.tif 000005
【００２８】
式（４）では、μｓは、シャドウ腐食バウの量であり、蓄積されたチャネル照射線量にも依存する。チャネル照射線量は、そのチャネルが運転中の炉心にどれほど長く滞留したかを表す。この関係の係数Ａ３、Ａ４、Ａ５、．．．およびチャネル照射線量依存性は、予測されたチャネル変形と測定されたチャネル変形の比較から決定され、異なる原子炉クラス、セル幾何形状、および水の化学的性質の環境について異なる場合がある。
【００２９】
炉心平均セル平均バウ（ＢＯＷＡＶＧ）値は、通常、ブレードウィングに面するすべてのチャネル面のすべての軸方向高度からの最大バウ値を使用して計算される。ＢＯＷＡＶＧ計算は、当技術分野で周知であり、プラントのライセンス付与に必要である。
【００３０】
図８では、チャネル壁圧力降下および総バルジ（弾性およびクリープ）を、ＰＡＮＡＣＥＡ（登録商標）などの炉心シミュレータで経験的基礎および／または数学的／物理学的基礎と共に諸式を使用して、チャネルバルジ変形に影響する計算された核運転条件および機械的運転条件を直接に反映することによって計算する［Ｓ３００］ことができる。たとえば、チャネル高速フルエンス、チャネル壁圧力降下、およびチャネル照射線量は、チャネルバルジ変形に影響する計算された核運転条件および機械的運転条件を表す。計算されたチャネル高速フルエンス値、チャネル壁圧力降下値、およびチャネル照射線量値は、バルジの弾性成分およびクリープ成分を計算するのに使用される。総バルジは、弾性成分およびクリープ成分の合計である。
【００３１】
Ｓ３００ａに示されているように、各制御セル内で、その制御セルの各軸方向高度で、機能Ｓ３０１からＳ３０６を実行する。各軸方向高度で、機能Ｓ３０１からＳ３０５ｂを、その軸方向位置でブレードウィングに面する８つの面（４つのチャネルについて１チャネルあたり２つ）のそれぞれについて実行する。
【００３２】
制御セル内の所与の軸方向高度の面ごとに、チャネルゆがみのすべての源（図６のＳ１１０、Ｓ１２０、またはＳ１３０のうちの１つからの結果と、図７のＳ２１３、Ｓ２２３、またはＳ２３３のうちの１つからの結果）を、機能Ｓ３１０で組み合わせて、総変形を得る。一般に、図８の機能Ｓ３０２からＳ３０７は、この総変形を使用して、変形の不確実性を計算し（図８のＳ３０３ａ、Ｓ３０５ｂ、およびＳ３０５ｄ）、公称のおよび統計的な上界チャネル−ブレード干渉ならびに結果の軸方向摩擦力をも計算する（図８のＳ３０３、Ｓ３０４、Ｓ３０５、Ｓ３０５ｂ、Ｓ３０５ｃ、およびＳ３０６）。セル摩擦力を、すべてのチャネル面について各軸方向高度で計算し（Ｓ３０６）、１つのセル内のすべての軸方向高度からの最大の力を、そのセルのセル摩擦測定基準として採用する（Ｓ３０７）。
【００３３】
チャネル高速フルエンスパラメータおよびチャネル被制御運転パラメータの計算（機能Ｓ１１０、Ｓ１２０、およびＳ１３０）と、高速フルエンスチャネルバウおよびシャドウ腐食によって誘導されたバウの計算（Ｓ２１０、Ｓ２２０、Ｓ２３０）と、バウの不確実性の計算（Ｓ２１２、Ｓ２２２、Ｓ２３２）と、弾性バルジおよびクリープバルジの計算（Ｓ３００）と、セル摩擦測定基準の計算（機能Ｓ３０１、Ｓ３０２、Ｓ３０３、Ｓ３０４、Ｓ３０５、Ｓ３０５ｃ、Ｓ３０６、およびＳ３０７）と、ＣＦＭの不確実性の計算（機能Ｓ３０３ａ、Ｓ３０５ａ、Ｓ３０５ｂ、Ｓ３０５ｄ）は、理論的基礎または経験的基礎を有する一般化された式に基づく。そのような式は、原子炉技術の当業者に既知であるか、原子炉技術の当業者によってたやすく導出され、したがって、説明を簡潔にするために省略する。
【００３４】
図９は、特定のチャネル応用例に関する支持する技術的基礎の例を示すための、測定されたチャネルバウ対予測された高速フルエンス勾配によって誘導されたバウのグラフである。図９は、予測された核運転条件およびチャネル材料の測定された物理的特性に基づいてＰＡＮＡＣＥＡ（登録商標）炉心シミュレータによって予測されたチャネルバウに対する測定されたチャネルバウの比較を表す。図９の実線は、予測値と測定値との間の完全な一致を表す。記号は、チャネル面ごとの予測されたバウと測定されたバウの実際の比較を表し、実線の周囲の記号の散乱は、バウ予測の不確実性を表す。
【００３５】
したがって、上で説明したように、使用可能な技術的現状のモデルが、チャネルゆがみの各成分を計算する。この例の方法が、高精度炉心シミュレーションコード（ＰＡＮＡＣＥＡ（登録商標）炉心シミュレータなど）に結合され、チャネルゆがみ計算への入力と強いて炉心シミュレーション結果を使用するように構成される場合に、改善された頑健性を達成することができる。たとえば、チャネル壁高速フルエンスおよびチャネル壁圧力降下を、チャネルが運転される時に、明確に確立された経験的基礎および／または数学的／物理学的基礎を有する式からなる一般的方法の組を使用して計算することができる。そのような式は、原子炉技術の当業者に周知であり、したがって、説明を簡潔にするために省略する。この例の方法での炉心シミュレーション結果のそのような使用は、計算されたチャネル変形がそのチャネルの実際の運転または計画された運転を反映することを保証する。
【００３６】
次に、チャネルゆがみの異なる成分を一緒に加算して、チャネルの各面での総変形の最良の推定値を与えることができる（図８のＳ３０１を参照されたい）。この最良の推定値を、それに関連する不確実性と組み合わせて、安全性パラメータの評価用のチャネルゆがみ入力をもたらすことができる。たとえば、炉心平均セル平均バウ（ＢＯＷＡＶＧ）およびその不確実性が、安全限界最小限界出力比（ＭＣＰＲ）計算への入力のために計算される。もう１つの例では、制御セル内の８つのチャネルのゆがみすなわち、ブレードウィングに面する８つのチャネル面（４つのチャネルについて１チャネルあたり２面）のゆがみを組み合わせて、チャネル−制御ブレード干渉とそのブレードへの結果の軸方向摩擦力を入手することができる。
【００３７】
最良の推定チャネル面ゆがみの計算を用いて、公称のまたは期待されるセル摩擦力を計算することができる（図８のＳ３０５ｃ）。計算されたパラメータに不確実性を含めることによって（図８のＳ３０５ｂ）、セル摩擦力の統計に基づく上界値（図８のＳ３０６）を、下の例の式の組（５）に示されているように計算することができる。
ＦＵｐｐｅｒ  ＝ＦＮｏｍｉｎａｌ＋ＴσＦ      （５）
ただし
ＦＵｐｐｅｒ  ＝統計に基づく上界制御ブレード摩擦力
ＦＮｏｍｉｎａｌ＝公称制御ブレード摩擦力
Ｔ        ＝統計的係数
σＦ      ＝セル摩擦の不確実性
式の組（５）では、セル摩擦力の不確実性（σＦ）を、モンテカルロシミュレーションまたは標準的な誤差の伝搬などの従来の統計的方法を使用して、既知の入力パラメータ不確実性に基づいて判定することができる。統計的係数Ｔは、所望のレベルの統計的信頼を提供するように選択され、高い摩擦を有する制御ブレードに関する原子炉内の経験の特性記述を基礎として決定することができる。すべての軸方向高度からのセル摩擦力（ＦＵｐｐｅｒ）の最大値が、セル摩擦測定基準（ＣＦＭ）として採用される（図８のＳ３０７を参照されたい）。
【００３８】
統計的係数Ｔは、機能Ｓ３０５ｂおよびＳ３０６での計算に含めることができる。統計的係数Ｔは、この例の方法の結果の表現の統計的信頼を高めるために、高い干渉および高い摩擦を有するブレードに関する使用可能な産業界の経験を利用する。したがって、この例の方法は、干渉および摩擦力の予測について使用可能な最良の態様を利用し、それと同時に、問題制御セルに関する実際の産業界の経験を反映することによって高いレベルの統計的な確かさを提供する。
【００３９】
図１０Ａおよび１０Ｂに、エネルギサイクルを介してＢＷＲ炉心が運転される時の炉心平均セル平均バウとその標準偏差とを示す。セル平均バウは、ある制御セル内の８つのチャネル面の最大バウの平均値として計算される。いくつかの例で、セルを、ブレード−センタードセルではなく「インスツルメント−センタード」セルと見なすことができる。というのは、ＢＷＲにおいて、いくつかの位置が、４つのバンドルによって囲まれたプラント計装チューブを有するからである。しかし、同一の概念が、インスツルメント−センタードセルにあてはまる。
【００４０】
炉心平均セル平均バウは、セル平均バウの統計的サンプリングを行うことによって計算され、ここで、このサンプリングには、制御セル内の８つのチャネル面の最大バウに関する不確実性が含まれる。さらに、このプロセスには、寿命のどの時にも指定された最小値を超えるバンドル出力を有するセルだけが統計的サンプリングに含まれるという要件を含めることができる。図１０Ａを参照すると、炉心平均セル平均バウＢＯＷＡＶＧ（１）および図１０Ｂのその標準偏差ＢＯＷＡＶＧ（２）は、炉心が燃焼する時のサイクル照射線量の関数として、ｍｉｌ（１インチの１／１０００）単位で示されている。サイクルの始め（ＢＯＣ）での負の約−１５ｍｉｌからサイクルの終り（ＥＯＣ）での正の２４ｍｉｌまでのＢＯＷＡＶＧ（１）の漸進的な遷移が観察される。それに対応して、ＢＯＷＡＶＧ（２）は、ＢＯＣでの約５ｍｉｌからＥＯＣでの約３８ｍｉｌまで増加する。
【００４１】
図１０Ａでは、チャネルバウの「最良推定値」を、関連する不確実性およびその不確実性の中での統計的サンプリングと共に、ＢＯＷＡＶＧを得るのに使用することができ、これらの値は、高い度合の統計的な確かさを有する。高速フルエンス勾配によって誘導されたバウとシャドウ腐食によって誘導されたバウの両方を、この分析に含めることができる。これらの値は、原子炉炉心の設計、最適化、ライセンス付与、および監視の他の態様に関連する入力である。たとえば、バンドル出力計算を含む、中性子計算に使用される核物質断面は、チャネルバウの関数である。もう１つの例で、安全限界最小限界出力比（ＭＣＰＲ）計算は、高速フルエンス勾配によって誘導されたバウおよびシャドウ腐食によって誘導されたバウの不確実性を利用する。
【００４２】
図１１に、ＢＷＲ炉心のセル摩擦を査定する、例の出力表示を示す。この例では、図１１は、ＢＷＲ炉心内のすべての制御セルのセル摩擦測定基準マップを示す。セルは、この図ではレベル０〜３として示されたＣＦＭ重大度の異なるレベルを示すために、カラーコーディング（または別の判断基準によって区別）することができる。同一の情報を、制御セルごとに重大度キーを含めることによって表示することができる。
【００４３】
１つの例の応用例で、結果データを、ユーザコマンドを介して、査定用の所望の表示結果に出力することができる。もう１つの例で、結果データを、炉心内の全制御セルのあるレベルを超えるＣＦＭを計算し、緩和するための自動化されたプロセスでこの例の方法を実施するように設計されたコンピュータプログラムの組の一部として保管することができる。図１１には、水平の棒または垂直の棒によって、運転中に遭遇した実際の問題セルも示されている。この視覚的な例では、水平の棒は、セル摩擦が垂直の棒を有する制御セルと比較してより重大である制御セルを表す。
【００４４】
図１１を参照すると、各番号は、制限する値に対して正規化された、その制御セルの計算されたＣＦＭ値を表す。ＣＦＭレベル３は、レベル１より重大である。たとえば、レベル３は、動作不能ブレードの設計限度または重大度閾値とすることができ、レベル１は、安定しないブレード（ｎｏ−ｓｅｔｔｌｅ  ｂｌａｄｅ）の可能性に関する重大度閾値とすることができる。レベル０は、非常に高いレベルの統計的信頼で、その制御セルに高いセル摩擦の危険性がないことを示す、最も重大でないレベルである。
【００４５】
図１１に示されているように、ユーザまたは炉心設計者は、計算されたＣＦＭ値に基づいて、異なるレベルの重大度を有するチャネル（バンドル）をグループ化することができる。そのようなグループ化は、各グループの緩和処置を定式化する際にユーザまたは炉心設計者を援助することができる。たとえば、理論に基づく閾値を使用して、異なるレベルをセットすることができる。ユーザは、たとえばプラント固有運転経験の入手可能性に基づいて、経験に基づく閾値を使用するオプションも有する。理論に基づく閾値は、より高い経験に基づく閾値を使用することに信頼がある時に、ユーザまたは炉心設計者によって調整できる保守的な推定値として採用されるはずである。
【００４６】
問題セルの個数を減らすことは、チャネルゆがみの懸念およびセル摩擦の懸念を管理することの経済的側面に重大な影響を有する。これに関して、この例の方法は、ユーザが望む信頼の度合に基づいてそれらの懸念を緩和するための柔軟性を提供する。ユーザは、図１１の例では各セルのＣＦＭだけを見るが、ＣＦＭに寄与する異なるチャネルゆがみ成分に基づいて１つまたは複数の特定のセルの問題の大きさを査定する、炉心シミュレータ出力ファイル内でユーザが使用でき、アクセスできる詳細な情報がある。
【００４７】
たとえば、ＣＦＭに関する詳細な軸方向情報が、セルの異なる軸方向高度での問題の重大度を査定するのに利用可能であり、ここで、軸方向高度は、炉心シミュレーションのために事前にセットされている。さらに、ＣＦＭに寄与する変形および不確実性に関する詳細な情報が、制御セル内の４つのチャネルのそれぞれについて炉心シミュレータ出力ファイル内で利用可能である。この情報は、原子炉サイクル炉心の設計、最適化、ライセンス付与、および監視のさまざまなステージ中に緩和処置を定式化し、実施する際にユーザを援助することができる。１例として、炉心設計ステージ中および／または炉心最適化ステージ中に（ここで、炉心構成が設計されまたはモデル化されつつある）、チャネル（バンドル）を、問題チャネルを単一のセル内に一緒に置くか再チャネル化のために「犠牲セル」を作成するようにセルのグループに置くように配置することができる。代替案では、すべてのセルの干渉および摩擦を最小にする努力において、問題チャネル（バンドル）を炉心全体に分散させることができる。所望の制御ブレード監視戦略を確立するために、この例の方法は、炉心が燃焼する時に、確立された判断基準に基づく単純で柔軟な形で、疑わしいセルを識別する能力を提供することができる（たとえば、図１１に示されているように）。
【００４８】
したがって、本明細書で説明したように、原子炉の制御セルのセル摩擦測定基準を判定するこの例の方法は、１セル内の１つまたは複数のチャネルからのゆがみの影響を緩和することができる。というのは、このゆがみが、セル制御ブレードとの干渉を生じ、おそらくは制御セル軸方向摩擦と制御ブレード移動を妨げることとにつながるからである。
【００４９】
この方法では、特定のチャネル軸方向高度でのチャネル面変位を、チャネルバウおよびチャネルバルジを含むがこれらに限定されない、チャネル材料の既知の物理的特性およびチャネル運転条件から計算することができる。制御セルに関する例では、チャネル面変位を、個々の制御セル内の制御ブレードウィングに隣接する８つのチャネル面のそれぞれについて計算することができる。個々の制御ブレードウィングとその２つの隣接するチャネルとの間の、各軸方向高度での干渉は、その軸方向高度での計算されたチャネル面変位から計算される。チャネル−制御ブレード干渉を、同様に、セル内の制御ブレードウィングごとに計算して、各軸方向高度での総干渉を判定することができる。
【００５０】
特定の軸方向高度でのセル摩擦を、その軸方向高度での計算された総チャネル−制御ブレード干渉から、ならびに既知のチャネル剛性と嵌合するチャネル材料および制御ブレード材料の既知の摩擦係数とから、計算することができる。次に、制御セル内の複数の軸方向高度のそれぞれでのそのような総変形、干渉、およびセル摩擦を、判定することができる。すべての軸方向高度での計算されたセル摩擦力の最大値が、その制御セルのセル摩擦測定基準として選択される。
【００５１】
したがって、この例の方法は、運転中の原子炉サイクル炉心の設計、最適化、ライセンス付与、および監視のさまざまなステージ中にチャネルゆがみの懸念およびセル摩擦の懸念に焦点を合わせる手段を提供することができる。チャネルゆがみおよびセル摩擦の重大度を定量化する（セル摩擦測定基準によって）能力は、炉心設計判断を行うための、ならびに／または計算されたおよび計画されたチャネル運転に基づいて緩和処置を講じるための基礎を提供する。この例の方法を、炉心の設計、最適化、ライセンス付与、および監視に使用される既存の方法にリンクし、これにコーディングすることができる。たとえば、この方法を、炉心の設計、最適化、ライセンス付与、および監視に使用される既存コンピュータプログラムの組の一部であるコンピュータソフトウェアモジュール内で実装することができる。
【００５２】
したがって、本発明の例の実施形態を説明したので、この実施形態を多数の形で変更できることは明白である。そのような変形形態は、本発明のこの例の実施形態の趣旨および範囲からの逸脱と見なされてはならず、当業者に明白な修正形態のすべてが、添付の特許請求の範囲に含まれることが意図されている。
【図面の簡単な説明】
【００５３】
【図１】燃料チャネルを示すための例の燃料バンドルを示す上面図である。
【図２】ＢＷＲの制御セルを見下ろして示す２Ｄ上面図である。
【図３】制御セル内のチャネルゆがみの２つの機構を示す図である。
【図４】制御セル内のセル摩擦を示す図である。
【図５】ＢＷＲの制御セルのセル摩擦測定基準を判定する全般的な方法を示す処理流れ図である。
【図６】図５に示された方法のプロセスフロー機能をより詳細に示す図である。
【図７】図５に示された方法のプロセスフロー機能をより詳細に示す図である。
【図８】図５に示された方法のプロセスフロー機能をより詳細に示す図である。
【図９】特定のチャネル応用例に関する支持する技術的基礎の例を示すための、測定されたチャネルバウ対予測された高速フルエンス勾配によって誘導されたバウを示すグラフである。
【図１０Ａ】運転サイクル中に炉心が進行する時の炉心平均セル平均バウとその標準偏差とを示す図である。
【図１０Ｂ】運転サイクル中に炉心が進行する時の炉心平均セル平均バウとその標準偏差とを示す図である。
【図１１】原子炉サイクル炉心の設計、最適化、ライセンス付与、および監視中にセル摩擦を査定し、緩和する際に使用される、ＢＷＲ炉心内のすべての制御セルのセル摩擦測定基準の例の出力表示を示す図である。
【符号の説明】
【００５４】
１１０  燃料バンドル
１１５  棒グリッド位置
１２０  燃料チャネル
２００  制御セル
（ａ１）  制御セルの２Ｄ上面図
（ｂ１）  制御セルの２Ｄ上面図
２００’  制御セル
２１０  燃料アセンブリ
２２０  燃料チャネル
２２０’  燃料チャネル
２２５  チャネル面
（ａ２）  大写しの正面図
（ｂ２）  大写しの正面図
２３０  制御ブレード
２３０’  制御ブレード
２３５  ブレードウィング
４１０  ２つの面に隣接するチャネルが接する点
【図１】
2008203162.tif 000006
【図２】
2008203162.tif 000007
【図３】
2008203162.tif 000008
【図４】
2008203162.tif 000009
【図５】
2008203162.tif 000010
【図６】
2008203162.tif 000011
【図７】
2008203162.tif 000012
【図８】
2008203162.tif 000013
【図９】
2008203162.tif 000014
【図１０Ａ】
2008203162.tif 000015
【図１０Ｂ】
2008203162.tif 000016
【図１１】
2008203162.tif 000017