燃料集合体の置き違い事故の分析方法
【課題】少ない分析対象から分析する燃料集合体の置き違い事故の分析方法を提供する。
【解決手段】ステップ１１で単一の燃料集合体の置き違いの状況を計算し、ステップ１２で各個別の状況に対して炉心全体の燃料消耗を計算し、ステップ１３で置き違いの組み合わせに対応する二つの個別の置き違いの限界出力比のマトリックスを検出し、二組のＣＰＲマトリックスを重ね合わせると、一つの置き違いの組み合わせのＣＰＲマトリックスを得ることができ、ステップ１４で置き違いの各組み合わせのＣＰＲマトリックスから最小限界出力比の最大変化量を求め、ステップ１５で上述の置き違いの組み合わせをＭａｘＤＣＰＲに基づいて数値の大きいものから小さい順に１０番目までの置き違いの組み合わせを選び出し、炉心全体の燃料消耗を再度計算して、１０個の組み合わせから最大のＭａｘＤＣＰＲ値を検出し、最終的な分析結果を得る。
【解決手段】ステップ１１で単一の燃料集合体の置き違いの状況を計算し、ステップ１２で各個別の状況に対して炉心全体の燃料消耗を計算し、ステップ１３で置き違いの組み合わせに対応する二つの個別の置き違いの限界出力比のマトリックスを検出し、二組のＣＰＲマトリックスを重ね合わせると、一つの置き違いの組み合わせのＣＰＲマトリックスを得ることができ、ステップ１４で置き違いの各組み合わせのＣＰＲマトリックスから最小限界出力比の最大変化量を求め、ステップ１５で上述の置き違いの組み合わせをＭａｘＤＣＰＲに基づいて数値の大きいものから小さい順に１０番目までの置き違いの組み合わせを選び出し、炉心全体の燃料消耗を再度計算して、１０個の組み合わせから最大のＭａｘＤＣＰＲ値を検出し、最終的な分析結果を得る。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、燃料集合体の置き違い事故の分析方法に関し、特に沸騰水型軽水炉における燃料集合体の置き違いの事故を分析する燃料集合体の置き違い事故の分析方法に関する。
【背景技術】
【０００２】
原子力発電所で使用する核燃料は、一定の運転期間を経過すると、次の運転に備えて新しい燃料に交換しなければならない。炉心には数百体もの燃料集合体が設けられている。よって、新しい燃料に交換する場合、複数の燃料集合体を出し入れしなければならない。
【０００３】
燃料集合体を置き違えることによって事故が発生する確率は低い。しかしながら、一度事故が発生すると、原子力発電の安全に重大な影響を与える。よって、原子力発電の安全から言うと、それぞれの燃料レイアウトについていずれも実際に運転する前に燃料集合体の置き違い事故を想定した安全のための分析を行い、原子炉のための分析を行い、原子炉の安全を確実に保ち、燃料集合体の置き違いによる沸騰水型軽水炉の事故発生を防がなければならない。
【０００４】
従来の燃料集合体の置き違い事故の分析は、エンジニアの経験に基づき「トライ・アンド・エラー（ｔｒｙ−ａｎｄ−ｅｒｒｏｒ）」という方法で燃料集合体の最も深刻な置き違いの状況を検出していた。しかしながら、該「トライ・アンド・エラー」の分析方法は、分析に多くの人手を必要とするとともに、最も深刻な燃料集合体の置き違いの状況を必ず導き出せるわけではなかった。よって、その後、「簡易選別法（ｓｉｍｐｌｅ ｆｉｌｔｅｒｉｎｇ ｍｅｔｈｏｄ）」、「パラメーター分析法（ｐａｒａｍｅｔｅｒ ａｎａｌｙｓｉｓ ｍｅｔｈｏｄ）」及びフォークス式検索法（Ｆａｗｋ’ｓ ｓｅａｒｃｈｉｎｇ ｍｅｔｈｏｄ）（アメリカ合衆国特許公報５，７９０，６１８号）などを含むその他の分析方法が開発された。
【０００５】
該「簡易選別法」は、二次元の炉心分析方法か、もしくはハリング燃料消耗法（Ｈａｌｉｎｇ ｐｒｉｎｃｉｐｌｅ ｔｏ ａｎａｌｙｚｅ）など簡易な炉心の計算方法で、最も深刻な燃料集合体の置き違いの状況を分析し、三次元の炉心分析ツールで比較的正確に分析を行い、最も深刻な燃料集合体の置き違いの状況を求める。
【０００６】
「パラメーター分析法」は、最も深刻な燃料集合体の置き違いの状況に関連する物理的なパラメーターを定義した後、それぞれの燃料集合体のパラメーターの値の排列に基づいて最も深刻な燃料集合体の置き違いの状況を決定する。
【０００７】
フォークス式検索法は１９９８年、特許登録した。フォークス式検索法は、スワップバリュー（ＳｗａｐＶａｒｉａｂｌｅ）によって最も深刻な燃料集合体の置き違いの状況を検索するための範囲を定義する。フォークス式検索法は、初期段階（Ｉｎｉｔｉａｌｉｚａｔｉｏｎ Ｐｈａｓｅ）と、検索段階（Ｓｅａｒｃｈ Ｐｈａｓｅ）の２つの段階に分けることができる。
【０００８】
初期段階では、燃料集合体が炉心の効率における貢献度がスワップバリュー（Ｓｗａｐ Ｖａｌｕｅ）を決定する。該スワップバリューはそれぞれ、炉心の燃料集合体１つ１つに対応している。検索段階では、任意で２つのスワップバリューを選択し、該スワップバリューに対応する燃料集合体の位置を交換して燃料集合体の置き違えた状況を発生させる。その後、炉心の分析ツールで燃料集合体の配置を誤った状況が与える影響を分析する。さらに、次の一対のスワップバリュー値を交換する。変換のルールは軽水炉により決定され、異なる軽水炉はそれぞれ異なる変換ルールを具えており、且つエンジニアの経験によりこられのルールが決定される。
【０００９】
「簡易選別法」もしくは「パラメーター分析法」に関わらず、いずれも正確性に欠ける。よって、最も深刻な燃料集合体の置き違いの状況を含め、三次元の炉心分析のための計算を大量に行わなければならないため、作業時間が延長され、実用的でない。また、フォークス式検索法はスワップバリューを調整する場合に共通のルールがなく、異なる軽水炉がそれぞれ独自の調整ルールを具え、且つエンジニアの経験によりルールが決定される。基本的にフォークス式は検索法に属し、燃料集合体の配置の誤りによる最も深刻な状況を確実に求めることができない。フォークス式の分析方法で最も深刻な燃料集合体の置き違いの状況を検索する確率を高めようとする場合、検索過程におけるスワップバリューの数、及び三次元の炉心分析の時間を延長しなければならないため、実用的ではない。
【００１０】
上述のそれぞれの分析方法は燃料集合体の置き違いの最も深刻な状況を求めることができるが、非常に長い作業時間を必要とする。沸騰水型軽水炉を例にとると、炉心には数百体もの燃料集合体が納められている。よって、すべての燃料集合体の置き違いの状況の全周期の最小限界出力比（Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭＣＰＲ）の最大変化量を分析すると、分析対象件数が膨大な数になり、非常に多くの時間を必要とする。よって、従来の方法では、実用性に欠ける。
【特許文献１】アメリカ合衆国特許公報５，７９０，６１８号
【発明の開示】
【発明が解決しようとする課題】
【００１１】
この発明は、分析対象となる燃料集合体の置き違いの数が少ない燃料集合体の置き違いの事故の分析方法を提供することを課題とする。
【課題を解決するための手段】
【００１２】
そこで本発明者は、従来の技術に鑑み鋭意研究を重ねた結果、分析対象となる燃料集合体の置き違いの数が少ない燃料集合体の置き違いの事故の分析方法であって、燃料集合体の置き違いの個別の計算結果を相互に重ね合わせてそれぞれの置き違いの結果を求める沸騰水型軽水炉の燃料集合体の置き違いの事故の分析方法であって、
該分析方法は少なくとも次ぎのステップにより実施し：ステップａではプログラムを利用して少なくとも単一の燃料集合体の置き違いの状況を計算し、ステップｂでは炉心の中性子の分析方法を利用して上述の燃料の置き違いのそれぞれの個別の状況に対して炉心全体の燃料消耗を計算し、ステップｃでは別のプログラムを介して燃料集合体の置き違いの組み合わせを発生させて、置き違いの組み合わせに対応する二つの個別の置き違いの限界出力比（Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＣＰＲ）のマトリックスを検出し、該二組のＣＰＲマトリックスを重ね合わせると、少なくとも一つの置き違いの組み合わせのＣＰＲマトリックスを得ることができ、
ステップｄでは燃料集合体の置き違いのそれぞれの組み合わせのＣＰＲマトリックスから最小限界出力比（Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭＣＰＲ）の最大変化量（Ｍａｘｉｍｕｍ ｏｆ Ｄｉｆｆｅｒｅｎｃｅ ｏｆ Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭａｘＤＣＰＲ）を求め、ステップｅでは上述の置き違いの組み合わせをＭａｘＤＣＰＲに基づいて数値の大きいものから小さい順に並べ、大きいものから１０番目までの置き違いの組み合わせを選び出し、炉心全体の燃料消耗を再度計算して、該１０個の置き違いの組み合わせの中から最大のＭａｘＤＣＰＲ値を検出し、該数値が燃料集合体の置き違いの事故の最終的な分析結果になる。
以下、この発明について具体的に説明する。
【００１３】
請求項１に記載の燃料集合体の置き違いの事故の分析方法は、燃料集合体の置き違いの個別の計算結果を相互に重ね合わせてそれぞれの置き違いの結果を求める沸騰水型軽水炉の燃料集合体の置き違いの事故の分析方法であって、該分析方法は少なくとも次ぎのステップにより実施し：ステップａではプログラムを利用して少なくとも単一の燃料集合体の置き違いの状況を計算し、ステップｂでは炉心の中性子の分析方法を利用して上述の燃料の置き違いのそれぞれの個別の状況に対して炉心全体の燃料消耗を計算し、ステップｃでは別のプログラムを介して燃料集合体の置き違いの組み合わせを発生させて、置き違いの組み合わせに対応する二つの個別の置き違いの限界出力比（Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＣＰＲ）のマトリックスを検出し、該二組のＣＰＲマトリックスを重ね合わせると、少なくとも一つの置き違いの組み合わせのＣＰＲマトリックスを得ることができ、ステップｄでは燃料集合体の置き違いのそれぞれの組み合わせのＣＰＲマトリックスから最小限界出力比（Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭＣＰＲ）の最大変化量（Ｍａｘｉｍｕｍ ｏｆ Ｄｉｆｆｅｒｅｎｃｅ ｏｆ Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭａｘＤＣＰＲ）を求め、ステップｅでは上述の置き違いの組み合わせをＭａｘＤＣＰＲに基づいて数値の大きいものから小さい順に並べ、大きいものから１０番目までの置き違いの組み合わせを選び出し、炉心全体の燃料消耗を再度計算して、該１０個の置き違いの組み合わせの中から最大のＭａｘＤＣＰＲ値を検出し、該数値が燃料集合体の置き違いの事故の最終的な分析結果になる。
【００１４】
請求項２に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法は、請求項１におけるステップａの個別の置き違いの状況の数値が、燃料集合体の個数に燃料形態の個数から１を引いた数をかけた数値になる。
【００１５】
請求項３に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法は、請求項１におけるステップａの個別の置き違いの状況の計算には次の段階を実施し：ａ１の段階では燃料集合体を正確に配置した状況の燃料形態の配置データを読み取り、ａ２の段階では炉心の四分の一の領域において、炉心のレイアウトの順序に基づいてそれぞれもともとの燃料形態とその他の燃料形態に分けてから置き換えて、それぞれの個別の置き違いの状況を計算する。
【００１６】
請求項４に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法は、請求項１におけるステップｃにおける燃料集合体の置き違いの組み合わせのＣＰＲマトリックスの形成には次の段階を実施し：ｃ１の段階では燃料集合体を正確に配置した状況の燃料集合体のレイアウトのデータを読み取り、ｃ２の段階では炉心の四分の一の領域において炉心のレイアウトに基づいて燃料形態を得て、再度、炉心全体のレイアウトに基づいて別の燃料形態を得て、置き違いの組み合わせを発生させて、ｃ３の段階では燃料集合体の置き違いの組み合わせに相対する第一及び第二の置き違いの組み合わせの炉心全体の燃料消耗を計算した結果を読み取り、第一及び第二のＣＰＲマトリックスをそれぞれ求め、正確に装荷した状況の計算結果を読み取り、第三のＣＰＲマトリックスを求め、ｃ４の段階では第二のＣＰＲマトリックスを第四、第五及び第六のＣＰＲマトリックスと置き換えて、第一のＣＰＲマトリックスに重ね合わせてから第三のＣＰＲマトリックスを引くと、少なくとも一つの燃料集合体の置き違いの状況の組み合わせのＣＰＲマトリックスを得る。
【００１７】
請求項５に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法は、請求項１におけるステップｄでは少なくとも次の段階を実施し：ｄ１の段階ではステップｃで得た燃料集合体の置き違いの状況のそれぞれの組み合わせのＣＰＲマトリックスからそれぞれの周期の燃料消耗点のＭＣＰＲ値を読み取り、ｄ２の段階では燃料集合体を正確に装荷した状況のＣＰＲマトリックスからそれぞれの周期の燃料消耗点のＭＣＰＲ値を読み取り、ｄ３の段階では燃料集合体の置き違いのそれぞれの状況の組み合わせのＣＰＲマトリックスにおけるそれぞれの周期の燃料消耗点のＭＣＰＲ値から、燃料集合体を正確に装荷した状況のＣＰＲマトリックスから相対する周期の燃料消耗点のＭＣＰＲ値を引くと、それぞれの周期の燃料消耗点のＭＣＰＲの変化量（Ｄｉｆｆｅｒｅｎｃｅ ｏｆ Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＤＣＰＲ）を求め、ｄ４の段階ではそれぞれの周期の燃料消耗点のＤＣＰＲ値からＭａｘＤＣＰＲを求める。
【００１８】
請求項６に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法は、請求項１における炉心の中性子の分析で使用する解析コードがＳＩＭＵＬＡＴＥ−３コードである。
【００１９】
請求項７に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法は、請求項１における二つのプログラムはＦｏｒｔｒａｎ言語から形成する。
【００２０】
請求項８に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法は、請求項１における沸騰水型軽水炉は従来の沸騰水型軽水炉、及び最新の沸騰水型軽水炉である。
【００２１】
請求項９に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法は、請求項１におけるステップｅではＭａｘＤＣＰＲの数値の大きいものから１００番目までの置き違いの状況の組み合わせを選び出し、炉心全体の燃料消耗を再度計算する。
【発明の効果】
【００２２】
この発明の燃料集合体の置き違いの分析方法は、分析する対象数が少ないため、分析する時間を短縮することができ、実用的で炉心の安全性を確保することができるという効果を有する。
【発明を実施するための最良の形態】
【００２３】
この発明は、分析対象となる燃料集合体の置き違いの数が少ない燃料集合体の置き違いの事故の分析方法であって、燃料集合体の置き違いの個別の計算結果を相互に重ね合わせてそれぞれの置き違いの結果を求める沸騰水型軽水炉の燃料集合体の置き違いの事故の分析方法であって、該分析方法は少なくとも次ぎのステップにより実施し：ステップａではプログラムを利用して少なくとも単一の燃料集合体の置き違いの状況を計算し、ステップｂでは炉心の中性子の分析方法を利用して上述の燃料の置き違いのそれぞれの個別の状況に対して炉心全体の燃料消耗を計算し、ステップｃでは別のプログラムを介して燃料集合体の置き違いの組み合わせを発生させて、置き違いの組み合わせに対応する二つの個別の置き違いの限界出力比（Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＣＰＲ）のマトリックスを検出し、該二組のＣＰＲマトリックスを重ね合わせると、少なくとも一つの置き違いの組み合わせのＣＰＲマトリックスを得ることができ、ステップｄでは燃料集合体の置き違いのそれぞれの組み合わせのＣＰＲマトリックスから最小限界出力比（Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭＣＰＲ）の最大変化量（Ｍａｘｉｍｕｍ ｏｆ Ｄｉｆｆｅｒｅｎｃｅ ｏｆ Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭａｘＤＣＰＲ）を求め、ステップｅでは上述の置き違いの組み合わせをＭａｘＤＣＰＲに基づいて数値の大きいものから小さい順に並べ、大きいものから１０番目までの置き違いの組み合わせを選び出し、炉心全体の燃料消耗を再度計算して、該１０個の置き違いの組み合わせの中から最大のＭａｘＤＣＰＲ値を検出し、該数値が燃料集合体の置き違いの事故の最終的な分析結果になる。
係る構成の燃料集合体の置き違いの事故の分析方法について、その構造と特徴を詳述するために具体的な実施例を挙げ、以下に説明する。
【実施例１】
【００２４】
図１は、この発明のフローチャートを示したものである。図面によれば、この発明は燃料集合体の置き違いの事故の分析方法であって、二つの個別の燃料集合体の置き違いの状況の限界出力比（Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＣＰＲ）のマトリックスと、燃料集合体を正しく装荷した状況のＣＰＲマトリックスとを交互に重ね合わせてそれぞれの配置を誤った組み合わせのＣＰＲマトリックスを求めた。また、この発明は対象となる炉心の四分の一の領域にのみ分析を行い、分析は次のステップで行った。
【００２５】
ステップａ（１１）では燃料集合体の置き違いの個別の状況に関する数値を計算した。プログラムを利用して少なくとも一つの個別の配置を誤った状況に関する数値を計算した。該プログラムはＦｏｒｔｒａｎ言語により作成した。個別の置き違いの状況の数は対象になる炉心の燃料集合体の個数の４分の１に燃料集合体の形態の個数から１を引いた数値をかけた数になった。
【００２６】
ステップｂ（１２）ではそれぞれの個別の燃料集合体の置き違いに対して炉心の燃料の消耗を計算した。炉心全体の中性子の分析方法を利用して、上述のそれぞれの個別の燃料集合体の置き違いの状況対する炉心全体の燃料の消耗を計算した。該炉心の中性子の分析方法に使用する解析コードはＳＩＭＵＬＡＴＥ−３コードとした。
【００２７】
ステップｃ（１３）では複数のＣＰＲマトリックスを重ね合わせると燃料集合体の置き違いの状況のＣＰＲマトリックスを得ることができた。別のＦｏｒｔｒａｎプログラムにより燃料集合体の置き違いの状況の組み合わせを作成し、燃料集合体の置き違いの状況の組み合わせに対応する二つの個別の燃料集合体の置き違いのＣＰＲマトリックスを求め、該二組のＣＰＲマトリックスと燃料集合体を正しく装荷した状況のＣＰＲマトリックスとを重ね合わせると、少なくとも一つの燃料集合体の置き違いの組み合わせのＣＰＲマトリックスを得ることができた。
【００２８】
ステップｄ（１４）ではそれぞれの燃料集合体の置き違いの状況の組み合わせのＭａｘＤＣＰＲを求めた。それぞれの燃料集合体の置き違いの状況の組み合わせの最小限界出力比（Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭＣＰＲ）の最大変化量（Ｍａｘｉｍｕｍ ｏｆ Ｄｉｆｆｅｒｅｎｃｅ ｏｆ Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭａｘＤＣＰＲ）を求めた。
【００２９】
ステップｅ（１５）では燃料集合体の置き違いの事故の分析結果を得た。該燃料集合体の置き違いの組み合わせを排列し、最終的な計算を行うと、燃料集合体の置き違いの事故の分析の最終結果を得た。
【実施例２】
【００３０】
この発明の分析方法を初装荷炉心に応用した場合を説明した。
【００３１】
図１はこの発明のフローチャートを示した。図１Ａはステップａのフローチャートを示した。図１Ｂはステップｃのフローチャートを開示した。図１Ｃはステップｄのフローチャートを示した。図１Ｄは実例１の初装荷炉心における燃料集合体のレイアウトを示した。
【００３２】
図面によれば、実例１では初装荷炉心１６に対して燃料集合体の置き違いの事故における分析を行った。この実例で分析する初装荷炉心１６は、対象となる炉心の縦、横軸が直交して四等分された領域の四分の一とした。該初装荷炉心１６は６２４体の燃料集合体を具えてなり、第一燃料形態１６１、第二燃料形態１６２及び第三燃料形態１６３など三種類の燃料形態を具えたものとした。この発明の分析方法は炉心の四分の一の領域１６４に対してのみ分析を行い、分析には次のステップを実施した。
【００３３】
ステップａではＦｏｒｔｒａｎ言語で作成されたプログラムで３１２体の燃料集合体の置き違いの状況の数値を計算した。計算は次の段階で行った。
【００３４】
（ａ１）の段階において燃料集合体を正確に配置した燃料形態の配置データを読み取った１１１のステップを行い、（ａ２）の段階においてもともとの燃料形態と他の二種類の燃料形態とを置き換えて、それぞれの炉心のに二つの個別の誤った配置を形成した１１２のステップを行った。炉心の四分の一の領域１６４に、炉心のレイアウトに基づいてもともとの燃料形態とその他の二種類の燃料形態とを置き換えると、それぞれの炉心のレイアウトに二つの個別の置き違いの状況が発生した。前述の方法で計３１２体の個別の燃料集合体の置き違いの状況の数値を計算した。
【００３５】
上述の（ａ２）の段階で個別の置き違いの数が炉心の四分の一の領域１６４の燃料集合体の個数に燃料形態の個数から１を引いた数値をかけた数値であることが得られた。よって、この実施例では３１２個（６２４／４×（３−１）＝３１２）の個別の燃料集合体の置き違いを具えたものとなった。
【００３６】
ステップｂではＳＩＭＵＬＡＴＥ−３コードの炉心の中性子を分析する方法で、上述の３１２体の置き違いのそれぞれの状況に対して炉心全体の燃料の消耗を計算した。
【００３７】
ステップｃではＦｏｒｔｒａｎ言語で作成したプログラムで燃料集合体の置き違いの状況の組み合わせを発生させて、３１２体の燃料集合体の置き違いの状況から計算した炉心全体の燃料の消耗の結果から相対する二つの個別の置き違いのＣＰＲマトリックスを求めた。また、該二組のＣＰＲマトリックスと正確に装荷した状況のＣＰＲマトリックスとを重ね合わせると、少なくとも一つの燃料集合体の置き違いの状況の組み合わせのＣＰＲマトリックスを得ることができ、マトリックスを得るために次の段階を実施した。
【００３８】
（ｃ１）の段階では燃料集合体を正確に配置した燃料形態の配置データを読み取った１３１のステップを実施した。
【００３９】
（ｃ２）の段階では燃料形態の情報を得た後、さらにレイアウトの順序に従って別の燃料形態を得て、置き違いの組み合わせを発生させ、炉心の四分の一の領域１６４における炉心のレイアウトの順序に基づいて燃料形態の情報を読み取り、炉心全体のレイアウトの順序に基づいて別の燃料形態の情報を読み取り置き違いの組み合わせを発生させた１３２のステップを実施した。
【００４０】
（ｃ３）の段階では第一及び第二ＣＰＲマトリックス及び正確に配置した状況の計算結果から第三のＣＰＲマトリックスを得た１３３のステップを実施した。誤った配置の組み合わせに対応する二つの個別の燃料集合体の置き違いの状況の計算結果を読み取り、第一及び第二のＣＰＲマトリックスを求めた。また、正確に配置した計算結果を読み取り第三のＣＰＲマトリックスを得た。
【００４１】
（ｃ４）では第二のＣＰＲマトリックスを置き換えて第四、第五及び第六のＣＰＲマトリックスを得たステップ１３４を実施した。また、第一のＣＰＲマトリックスを重ね合わせてから第三のＣＰＲマトリックスを引くと誤った配置の組み合わせのＣＰＲマトリックスを得ることができた。燃料集合体を置き違えた炉心の第二のＣＰＲマトリックスに対応するマトリックスを置き換え、第四、第五及び第六のＣＰＲマトリックス（追って説明する）を取得し、それぞれ第一のＣＰＲマトリックスに重ね合わせて、第三のＣＰＲマトリックスを引くと、少なくとも一つの誤った配置の組み合わせのＣＰＲマトリックスを得られた。
【００４２】
ステップｄではそれぞれの誤った配置の組み合わせのＭａｘＤＣＰＲ値を求めた。検出には次のステップを実施した。
【００４３】
（ｄ１）では、燃料集合体の置き違いのそれぞれの組み合わせのＭＣＰＲを読み取った１４１のステップを実施した。ステップｃで得た燃料集合体の置き違いのそれぞれの組み合わせのＣＰＲマトリックスから、各周期の燃料の消耗点のＭＣＰＲを読み取った（この実施例では計２７３３２個の置き違いの組み合わせを具えていたため、２７３３２個のＣＰＲマトリックスを摘出する必要があった。また全周期では４５個の燃料消耗点を具えていた）。
【００４４】
（ｄ２）では燃料集合体を正確に配置した状況のＭＣＰＲを得た１４２のステップを実施した。燃料集合体を正確に配置した状況のＣＰＲマトリックスからそれぞれの周期の燃料の消耗点におけるＭＣＰＲを読み取った。
【００４５】
（ｄ３）では燃料集合体の置き違いの状況の組み合わせのそれぞれのＭＣＰＲ値から燃料集合体を正確に装荷した状況のＭＣＰＲ値を引くと、燃料集合体の置き違いの状況の組み合わせのＤＣＰＲを得られた１４３のステップを実施した。（ｄ２）の段階のＭＣＰＲ値から（ｄ１）の段階のＭＣＰＲ値を引くと、燃料集合体の置き違いのそれぞれの状況の組み合わせの各周期の燃料の消耗点のＭＣＰＲの変化量を得ることができた（ＤＣＰＲ；Ｄｉｆｆｅｒｅｎｃｅ ｏｆ Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ）。
【００４６】
（ｄ４）の段階ではＤＣＰＲ値からＭａｘＤＣＰＲを求めた１４４のステップを実施した。（ｄ３）の段階のＤＣＰＲ値から最大値を求めると、燃料集合体の置き違いのそれぞれの状況の組み合わせの各周期の燃料の消耗点におけるＭａｘＤＣＰＲを得ることができた。
【００４７】
ステップｅでは上述の燃料集合体の置き違いの組み合わせの状況であって、そのＭａｘＤＣＰＲを大きな値から小さな値に並べ、最大値から１０番目までの数値か、もしくは１００番目までの誤った配置の組み合わせを選び出して、再度、炉心全体の燃料の消耗を計算した。また、該１０個、もしくは１００個の燃料集合体の置き違いの状況の組み合わせにおいて最大のＭａｘＤＣＰＲ値を求めた。該数値は燃料集合体の置き違いの事故の最終的な分析結果になった。
【００４８】
上述の方法は、燃料集合体の置き違いの状況を検出し、求められた燃料集合体の置き違いの状況に対して分析を行うため、軽水炉の運転前に実施する燃料集合体の置き違いの事故の安全性を分析する時間を縮小することができた。
【実施例３】
【００４９】
この発明の分析方法を再装荷した炉心に応用する場合を説明した。
【００５０】
図１Ｅは、再装荷した炉心の燃料集合体に応用した場合を示した。図面によれば、この発明の実例２は再装荷した炉心１７で燃料集合体を誤って配置した場合の事故を分析する方法とした。この実例では再装荷した炉心１７の四分の一の領域に対して分析を実施した。該再装荷した炉心１７は６２４体の燃料集合体を具えてなり、第一の燃料形態８１、第二の燃料形態８２及び第四の燃料形態８３など三種類の燃料形態を具えてなるようにした。
【００５１】
該第一の燃料形態８１及び第二の燃料形態８２は、第一次に消耗した燃料とした。よって、第一の燃料形態８１及び第二の燃料形態８２は燃料消耗の分布を有し、同一形態の燃料集合体はその燃料の消耗に差があるため、２体の同一形態の燃料集合体でも置き違いの状況の組み合わせが発生する可能性があった。これは初装荷炉心ではみられない状況であった。この実例の燃料形態は三種類とした。しかしながら、第一の燃料形態に見られる燃料消耗の分布は７．２５７〜１３．１９５ＧＷＤ／ＭＴ（ｇｉｇａｗａｔｔ−ｄａｙｓ ｐｅｒ ｍｅｔｒｉｃ ｔｏｎ）であって、燃料消耗の分布の範囲が比較的大きかった。よって第一の燃料形態は更に二種類の形態に分けることができた。一種類目の形態は燃料消耗の分布が１０ＧＷＤ／ＭＴよりも大きく、もう一種類の形態は燃料消耗の分布が１０ＧＷＤ／ＭＴよりも小さくなった。また、第二種の燃料形態の燃料消耗の分布は１０．８３１〜１１．３８７ＧＷＤ／ＭＴになり、分布の範囲が大きくなかった。このため更に細かく燃料形態を分ける必要がなく、この実施例では計４種類の燃料形態を具えていた。また、この発明の分析方法は、炉心の四分の一の領域１７４に対してのみ分析を行い、分析には下記のステップを実施した。
【００５２】
ステップａではＦｏｒｔｒａｎ言語で作成したプログラムで４６８体の個別の燃料集合体の配置を誤った状況を計算した。計算には次ぎの段階を実施した。
【００５３】
（ａ１）の段階では、燃料集合体を正確に配置した燃料形態のデータを読み取った。
【００５４】
（ａ２）の段階では炉心の四分の一の領域１７４において炉心のレイアウトの順序に基づいてもともとの燃料形態に分けて、その他の三種類の燃料形態に置き換え、それぞれの炉心のレイアウトに対して燃料集合体の置き違いの三つの状況を形成した。この方法により計４６８体の個別の置き違いの状況を計算した。
【００５５】
上述の（ａ２）の段階は個別の置き違いの状況の数値が炉心の四分の一の領域１７４の燃料集合体の数に燃料形態の数から１を引いた数をかけたものになった。よって、この実施例では４６８個（６２４／４×（４−１）＝４６８）の置き違いの状況が得られた。
【００５６】
ステップｂではＳＩＭＵＬＡＴＥ−３コードを炉心の中性子を分析する方法として、上述の４６８体の個別の置き違いのレイアウトにおける炉心全体の燃料の消耗を計算した。
【００５７】
ステップｃでは別のＦｏｒｔｒａｎ言語から作成したプログラムで燃料集合体の置き違いの組み合わせを発生させた。４６８体の燃料集合体の個別の置き違いの炉心全体の燃料消耗を計算した結果から相対する二つの個別の置き違いのＣＰＲマトリックスを求める。該二組のＣＰＲマトリックスと燃料集合体を正確に装填した状況のＣＰＲマトリックスを重ね合わせると、少なくとも一つの燃料集合体の置き違いの組み合わせのＣＰＲマトリックスが得られた。マトリックスを得るために次ぎの段階を実施した。
【００５８】
（ｃ１）の段階では燃料集合体を正確に配置した燃料形態の配置データを読み取った。
【００５９】
（ｃ２）の段階では炉心の四分の一の領域１７４において炉心のレイアウトの順序に基づいて一種類の燃料形態の情報を得た。その後、炉心全体の配置順序に基づいてもう一種類の燃料形態を得て燃料集合体を置き違えた状況の組み合わせを形成した。
【００６０】
（ｃ３）の段階では燃料集合体を置き違えた状況の組み合わせに対応する二つの個別の燃料集合体の置き違えた炉心全体の燃料の消耗を計算した結果を読み取り、第一及び第二のＣＰＲマトリックスをそれぞれ求めた。また、燃料集合体を正確に装荷した状況の炉心全体の燃料の消耗を計算した結果を読み取ると、第三のＣＰＲマトリックスを得られた。
【００６１】
（ｃ４）の段階では燃料集合体を置き違えた状況の炉心の位置により第二のＣＰＲマトリックスと対応するマトリックスのレイアウトの交換を行った。レイアウトの交換は第四、第五及び第六のＣＰＲマトリックス（追って説明する）とし、それぞれ第一のＣＰＲマトリックスと重ね合わせた後、第三のＣＰＲマトリックスを引くと、少なくとも一つの燃料集合体を置き違えた状況の組み合わせのＣＰＲマトリックスを得られた。
【００６２】
ステップｄではそれぞれの誤った配置の組み合わせのＭａｘＤＣＰＲ値を求めた。検出には、次ぎのステップを実施した。
【００６３】
（ｄ１）の段階ではステップｃで得たそれぞれの誤った配置の組み合わせのＣＰＲマトリックスから、各周期の燃料消耗点のＭＣＰＲ値を求めた（この実施例では３２１１６個の誤った配置の組み合わせを具えていたため、３２１１６個のＣＰＲマトリックスを求めた。また、全周期では計３５個の燃料消耗点を具えていた）。
【００６４】
（ｄ２）の段階では燃料集合体を正確に装荷した状況のＣＰＲマトリックスから、それぞれの周期の燃料消耗点のＭＣＰＲ値を読み取った。
【００６５】
（ｄ３）の段階では（ｄ２）の段階のＭＣＰＲ値からステップｄ１のＭＣＰＲ値を引くと、それぞれの置き違いの組み合わせの各周期の燃料消耗点のＤＣＰＲ値を得た。
【００６６】
（ｄ４）の段階では（ｄ３）の段階のＤＣＰＲ値から最大値を選び出すと、それぞれの誤った配置の組み合わせの各周期の燃料消耗点のＭａｘＤＣＰＲ値を得られた。
【００６７】
ステップｅでは上述の誤った配置の組み合わせのＭａｘＤＣＰＲを大きい数値から小さくなるように並べた。最も大きい数値から１０番目か、もしくは１００番目までの誤った組み合わせを選び出し、全ての炉心の燃料消耗を再度計算し、１０個、もしくは１００個の誤った組み合わせの中から最大のＭａｘＤＣＰＲ値を求めた。その数値は燃料集合体の設置を誤った事故の最終的な分析結果になった。
【００６８】
発生する可能性のある燃料集合体の置き違いの状況を検出し、燃料集合体を置き違えた状況に対して上述の方法で分析を行うと、軽水炉の運転前に実施する燃料集合体の配置の誤りによる事故の安全性に関する分析の時間を短縮することができた。
【００６９】
図２Ａは燃料集合体を正確に装荷した状況の燃料の配置を示した。図２Ｂは燃料集合体を置き違えた第一の状況を示した。図２Ｃは燃料集合体を置き違えた第一の個別の状況を示した。図２Ｄは燃料集合体を置き違えた第二の個別の状況を示した。図２Ｅは燃料集合体を置き違えた第三の個別の状況を示した。図２Ｆは燃料集合体を置き違えた第二の組み合わせを示した。図２Ｇは燃料集合体を置き違えた第四の個別の状況を示した。図２Ｈは燃料集合体を置き違えた第三の状況を示した。図２Ｉは燃料集合体を置き違えた第五の個別の状況を示した。図２Ｊは燃料集合体を置き違えた第四の状況を示した。
【００７０】
図面によれば、ステップｃにおいて複数のＣＰＲマトリックスを重ね合わせた方法を下記に示した。図２Ａでは、炉心の四分の一の対称となる領域２１に燃料集合体を正確に装荷したと仮定した。また、そのＣＰＲマトリックスは第三のＣＰＲマトリックス（ｆｏｌ）であって、次のように表示された。
【数１】
【００７１】
図２Ｂでは燃料集合体の置き違いの第一の状況１０を仮定した。また、図２Ｃは上述の燃料集合体の置き違いの状況の第一の組み合わせ２２１の第一の個別の誤ったレイアウト２３１であって、そのＣＰＲマトリックスは第一のＣＰＲマトリックス（Ｐ）だった。該ＣＰＲマトリックスを下記に表示した。
【数２】
【００７２】
図２Ｄでは上述の燃料集合体の第一の置き違いの状況の組み合わせ２２１の第二の個別の置き違いの状況２３２を示した。該ＣＰＲマトリックスは第二のＣＰＲマトリックス（Ｑ）であって、下記に表示した。
【数３】
【００７３】
よって、上述の第一のＣＰＲマトリックス（Ｐ）と第二のＣＰＲマトリックス（Ｑ）を加えて、第三のＣＰＲマトリックス（ｆｏ１）を引くと、第一の誤った配置の組み合わせのＣＰＲマトリックスを得ることができた。該ＣＰＲマトリックスを下記に表示した。
【数４】
【００７４】
第二のＣＰＲマトリックスはＹ軸を通して反射すると、第２Ｅに開示した第三の個別の置き違いの状況２３３のＣＰＲマトリックスを得ることができた。該ＣＰＲマトリックスは第四のＣＰＲマトリックスであって、第一のＣＰＲマトリックスに第四のＣＰＲマトリックスを加えた後、第三のＣＰＲマトリックスを引くと、第二の置き違いの状況の組み合わせ２２２（図２Ｆ参照）のＣＰＲマトリックスを得ることができた。該ＣＰＲマトリックスを下記に表示した。
【数５】
【００７５】
第二のＣＰＲマトリックスがＸ軸を通して反射すると、図２Ｇに開示する第四の個別の置き違いの状況２３４のＣＰＲマトリックスを得られた。該ＣＰＲマトリックスは第五のＣＰＲマトリックスであって、第一のＣＰＲマトリックスに第五のＣＰＲマトリックスを加えた後、第三のＣＰＲマトリックスを引くと、第三の置き違いの状況の組み合わせ２２３（図２Ｈ参照）のＣＰＲマトリックスを得られた。該マトリックスを下記に表示した。
【数６】
【００７６】
第二のＣＰＲマトリックスがＸ軸及びＹ軸を通して２度反射されると、図２Ｉに開示する第五の個別の置き違いの状況２３５のＣＰＲマトリックスを得ることができた。該ＣＰＲマトリックスは第六のＣＰＲマトリックスであって、第一のＣＰＲマトリックスに第六のＣＰＲマトリックスを加えて、第三のＣＰＲマトリックスを引くと、第四の置き違いの状況の組み合わせ２２４（図２Ｊ参照）のＣＰＲマトリックスを得られた。該マトリックスを下記に表示した。
【数７】
【００７７】
上述の方法を通して初装荷炉心を分析する場合、３１２個の単一置き違いのＣＰＲマトリックスと燃料集合体を正確に装荷したＣＰＲマトリックスとを重ね合わせた後、２７３３２個の置き違いの状況のＣＰＲマトリックスを得た。再装荷した炉心を分析する場合には、４６８個の個別の誤った配置のＣＰＲマトリックスを正確に装荷したＣＰＲマトリックスを重ね合わせて３２１１６個の置き違いの状況の組み合わせのＣＰＲマトリックスを得た。
【００７８】
上述のように、この発明の沸騰水型軽水炉における燃料集合体の置き違いの事故の分析方法は、従来の分析方法の欠点を有効に改善し、分析に必要な誤った組み合わせの数を減らし、運転前に実施する燃料集合体の置き違いの事故の安全性を分析する時間を短縮して、炉心の安全を確保することができた。
【００７９】
以上は、この発明の好ましい実施例であって、この発明の実施の範囲を限定するものではない。よって、当業者のなし得る修正、もしくは変更であって、この発明の精神の下においてなされ、この発明に対して均等の効果を有するものは、いずれもこの発明の特許請求の範囲に属するものとする。
【図面の簡単な説明】
【００８０】
【図１】この発明の沸騰水型軽水炉の燃料集合体を誤って配置した事故の分析方法を示したフローチャートである。
【図１Ａ】図１に開示する沸騰水型軽水炉の燃料集合体を誤って配置した事故の分析方法のステップａを示したフローチャートである。
【図１Ｂ】図１に開示する沸騰水型軽水炉の燃料集合体を誤って配置した事故の分析方法のステップｃを示したフローチャートである。
【図１Ｃ】図１に開示する沸騰水型軽水炉の燃料集合体を誤って配置した事故の分析方法のステップｄを示したフローチャートである。
【図１Ｄ】この発明の沸騰水型軽水炉の燃料集合体の実例１の初装荷炉心を示した断面図である。
【図１Ｅ】この発明の沸騰水型軽水炉の燃料集合体の実例２の再装荷した炉心を示した断面図である。
【図２Ａ】この発明の沸騰水型軽水炉に燃料集合体を正確に装荷した状況を示した断面図である。
【図２Ｂ】この発明の燃料集合体を誤って配置した第一の状況を示した断面図である。
【図２Ｃ】この発明の燃料集合体を誤って配置した第一の個別の状況を示した断面図である。
【図２Ｄ】この発明の燃料集合体を誤って配置した第二の個別の状況を示した断面図である。
【図２Ｅ】この発明の燃料集合体を誤って配置した第三の個別の状況を示した断面図である。
【図２Ｆ】この発明の燃料集合体を誤って配置した第二の状況を示した断面図である。
【図２Ｇ】この発明の燃料集合体を誤って配置した第四の個別の状況を示した断面図である。
【図２Ｈ】この発明の燃料集合体を誤って配置した第三の状況を示した断面図である。
【図２Ｉ】この発明の燃料集合体を誤って配置した第五の個別の状況を示した断面図である。
【図２Ｊ】この発明の燃料集合体を誤って配置した第四の状況を示した断面図である。
【符号の説明】
【００８１】
１６ 初装荷炉心
１７ 再装荷した炉心
１６１ 第一燃料形態
１６２ 第二燃料形態
１６３ 第三燃料形態
１６４ 炉心の四分の一の領域
１７１ 第一燃料形態
１７２ 第二燃料形態
１７３ 第四燃料形態
１７４ 炉心の四分の一の領域
２１ 炉心の四分の一の対称となる領域
２２１ 第一の置き違いの状況の組み合わせ
２２２ 第二の置き違いの状況の組み合わせ
２２３ 第三の置き違いの状況の組み合わせ
２２４ 第四の置き違いの状況の組み合わせ
２３２ 第二の個別の置き違いの状況
２３３ 第三の個別の置き違いの状況
２３４ 第四の個別の置き違いの状況
２３５ 第五の個別の置き違いの状況
【技術分野】
【０００１】
この発明は、燃料集合体の置き違い事故の分析方法に関し、特に沸騰水型軽水炉における燃料集合体の置き違いの事故を分析する燃料集合体の置き違い事故の分析方法に関する。
【背景技術】
【０００２】
原子力発電所で使用する核燃料は、一定の運転期間を経過すると、次の運転に備えて新しい燃料に交換しなければならない。炉心には数百体もの燃料集合体が設けられている。よって、新しい燃料に交換する場合、複数の燃料集合体を出し入れしなければならない。
【０００３】
燃料集合体を置き違えることによって事故が発生する確率は低い。しかしながら、一度事故が発生すると、原子力発電の安全に重大な影響を与える。よって、原子力発電の安全から言うと、それぞれの燃料レイアウトについていずれも実際に運転する前に燃料集合体の置き違い事故を想定した安全のための分析を行い、原子炉のための分析を行い、原子炉の安全を確実に保ち、燃料集合体の置き違いによる沸騰水型軽水炉の事故発生を防がなければならない。
【０００４】
従来の燃料集合体の置き違い事故の分析は、エンジニアの経験に基づき「トライ・アンド・エラー（ｔｒｙ−ａｎｄ−ｅｒｒｏｒ）」という方法で燃料集合体の最も深刻な置き違いの状況を検出していた。しかしながら、該「トライ・アンド・エラー」の分析方法は、分析に多くの人手を必要とするとともに、最も深刻な燃料集合体の置き違いの状況を必ず導き出せるわけではなかった。よって、その後、「簡易選別法（ｓｉｍｐｌｅ ｆｉｌｔｅｒｉｎｇ ｍｅｔｈｏｄ）」、「パラメーター分析法（ｐａｒａｍｅｔｅｒ ａｎａｌｙｓｉｓ ｍｅｔｈｏｄ）」及びフォークス式検索法（Ｆａｗｋ’ｓ ｓｅａｒｃｈｉｎｇ ｍｅｔｈｏｄ）（アメリカ合衆国特許公報５，７９０，６１８号）などを含むその他の分析方法が開発された。
【０００５】
該「簡易選別法」は、二次元の炉心分析方法か、もしくはハリング燃料消耗法（Ｈａｌｉｎｇ ｐｒｉｎｃｉｐｌｅ ｔｏ ａｎａｌｙｚｅ）など簡易な炉心の計算方法で、最も深刻な燃料集合体の置き違いの状況を分析し、三次元の炉心分析ツールで比較的正確に分析を行い、最も深刻な燃料集合体の置き違いの状況を求める。
【０００６】
「パラメーター分析法」は、最も深刻な燃料集合体の置き違いの状況に関連する物理的なパラメーターを定義した後、それぞれの燃料集合体のパラメーターの値の排列に基づいて最も深刻な燃料集合体の置き違いの状況を決定する。
【０００７】
フォークス式検索法は１９９８年、特許登録した。フォークス式検索法は、スワップバリュー（ＳｗａｐＶａｒｉａｂｌｅ）によって最も深刻な燃料集合体の置き違いの状況を検索するための範囲を定義する。フォークス式検索法は、初期段階（Ｉｎｉｔｉａｌｉｚａｔｉｏｎ Ｐｈａｓｅ）と、検索段階（Ｓｅａｒｃｈ Ｐｈａｓｅ）の２つの段階に分けることができる。
【０００８】
初期段階では、燃料集合体が炉心の効率における貢献度がスワップバリュー（Ｓｗａｐ Ｖａｌｕｅ）を決定する。該スワップバリューはそれぞれ、炉心の燃料集合体１つ１つに対応している。検索段階では、任意で２つのスワップバリューを選択し、該スワップバリューに対応する燃料集合体の位置を交換して燃料集合体の置き違えた状況を発生させる。その後、炉心の分析ツールで燃料集合体の配置を誤った状況が与える影響を分析する。さらに、次の一対のスワップバリュー値を交換する。変換のルールは軽水炉により決定され、異なる軽水炉はそれぞれ異なる変換ルールを具えており、且つエンジニアの経験によりこられのルールが決定される。
【０００９】
「簡易選別法」もしくは「パラメーター分析法」に関わらず、いずれも正確性に欠ける。よって、最も深刻な燃料集合体の置き違いの状況を含め、三次元の炉心分析のための計算を大量に行わなければならないため、作業時間が延長され、実用的でない。また、フォークス式検索法はスワップバリューを調整する場合に共通のルールがなく、異なる軽水炉がそれぞれ独自の調整ルールを具え、且つエンジニアの経験によりルールが決定される。基本的にフォークス式は検索法に属し、燃料集合体の配置の誤りによる最も深刻な状況を確実に求めることができない。フォークス式の分析方法で最も深刻な燃料集合体の置き違いの状況を検索する確率を高めようとする場合、検索過程におけるスワップバリューの数、及び三次元の炉心分析の時間を延長しなければならないため、実用的ではない。
【００１０】
上述のそれぞれの分析方法は燃料集合体の置き違いの最も深刻な状況を求めることができるが、非常に長い作業時間を必要とする。沸騰水型軽水炉を例にとると、炉心には数百体もの燃料集合体が納められている。よって、すべての燃料集合体の置き違いの状況の全周期の最小限界出力比（Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭＣＰＲ）の最大変化量を分析すると、分析対象件数が膨大な数になり、非常に多くの時間を必要とする。よって、従来の方法では、実用性に欠ける。
【特許文献１】アメリカ合衆国特許公報５，７９０，６１８号
【発明の開示】
【発明が解決しようとする課題】
【００１１】
この発明は、分析対象となる燃料集合体の置き違いの数が少ない燃料集合体の置き違いの事故の分析方法を提供することを課題とする。
【課題を解決するための手段】
【００１２】
そこで本発明者は、従来の技術に鑑み鋭意研究を重ねた結果、分析対象となる燃料集合体の置き違いの数が少ない燃料集合体の置き違いの事故の分析方法であって、燃料集合体の置き違いの個別の計算結果を相互に重ね合わせてそれぞれの置き違いの結果を求める沸騰水型軽水炉の燃料集合体の置き違いの事故の分析方法であって、
該分析方法は少なくとも次ぎのステップにより実施し：ステップａではプログラムを利用して少なくとも単一の燃料集合体の置き違いの状況を計算し、ステップｂでは炉心の中性子の分析方法を利用して上述の燃料の置き違いのそれぞれの個別の状況に対して炉心全体の燃料消耗を計算し、ステップｃでは別のプログラムを介して燃料集合体の置き違いの組み合わせを発生させて、置き違いの組み合わせに対応する二つの個別の置き違いの限界出力比（Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＣＰＲ）のマトリックスを検出し、該二組のＣＰＲマトリックスを重ね合わせると、少なくとも一つの置き違いの組み合わせのＣＰＲマトリックスを得ることができ、
ステップｄでは燃料集合体の置き違いのそれぞれの組み合わせのＣＰＲマトリックスから最小限界出力比（Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭＣＰＲ）の最大変化量（Ｍａｘｉｍｕｍ ｏｆ Ｄｉｆｆｅｒｅｎｃｅ ｏｆ Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭａｘＤＣＰＲ）を求め、ステップｅでは上述の置き違いの組み合わせをＭａｘＤＣＰＲに基づいて数値の大きいものから小さい順に並べ、大きいものから１０番目までの置き違いの組み合わせを選び出し、炉心全体の燃料消耗を再度計算して、該１０個の置き違いの組み合わせの中から最大のＭａｘＤＣＰＲ値を検出し、該数値が燃料集合体の置き違いの事故の最終的な分析結果になる。
以下、この発明について具体的に説明する。
【００１３】
請求項１に記載の燃料集合体の置き違いの事故の分析方法は、燃料集合体の置き違いの個別の計算結果を相互に重ね合わせてそれぞれの置き違いの結果を求める沸騰水型軽水炉の燃料集合体の置き違いの事故の分析方法であって、該分析方法は少なくとも次ぎのステップにより実施し：ステップａではプログラムを利用して少なくとも単一の燃料集合体の置き違いの状況を計算し、ステップｂでは炉心の中性子の分析方法を利用して上述の燃料の置き違いのそれぞれの個別の状況に対して炉心全体の燃料消耗を計算し、ステップｃでは別のプログラムを介して燃料集合体の置き違いの組み合わせを発生させて、置き違いの組み合わせに対応する二つの個別の置き違いの限界出力比（Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＣＰＲ）のマトリックスを検出し、該二組のＣＰＲマトリックスを重ね合わせると、少なくとも一つの置き違いの組み合わせのＣＰＲマトリックスを得ることができ、ステップｄでは燃料集合体の置き違いのそれぞれの組み合わせのＣＰＲマトリックスから最小限界出力比（Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭＣＰＲ）の最大変化量（Ｍａｘｉｍｕｍ ｏｆ Ｄｉｆｆｅｒｅｎｃｅ ｏｆ Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭａｘＤＣＰＲ）を求め、ステップｅでは上述の置き違いの組み合わせをＭａｘＤＣＰＲに基づいて数値の大きいものから小さい順に並べ、大きいものから１０番目までの置き違いの組み合わせを選び出し、炉心全体の燃料消耗を再度計算して、該１０個の置き違いの組み合わせの中から最大のＭａｘＤＣＰＲ値を検出し、該数値が燃料集合体の置き違いの事故の最終的な分析結果になる。
【００１４】
請求項２に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法は、請求項１におけるステップａの個別の置き違いの状況の数値が、燃料集合体の個数に燃料形態の個数から１を引いた数をかけた数値になる。
【００１５】
請求項３に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法は、請求項１におけるステップａの個別の置き違いの状況の計算には次の段階を実施し：ａ１の段階では燃料集合体を正確に配置した状況の燃料形態の配置データを読み取り、ａ２の段階では炉心の四分の一の領域において、炉心のレイアウトの順序に基づいてそれぞれもともとの燃料形態とその他の燃料形態に分けてから置き換えて、それぞれの個別の置き違いの状況を計算する。
【００１６】
請求項４に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法は、請求項１におけるステップｃにおける燃料集合体の置き違いの組み合わせのＣＰＲマトリックスの形成には次の段階を実施し：ｃ１の段階では燃料集合体を正確に配置した状況の燃料集合体のレイアウトのデータを読み取り、ｃ２の段階では炉心の四分の一の領域において炉心のレイアウトに基づいて燃料形態を得て、再度、炉心全体のレイアウトに基づいて別の燃料形態を得て、置き違いの組み合わせを発生させて、ｃ３の段階では燃料集合体の置き違いの組み合わせに相対する第一及び第二の置き違いの組み合わせの炉心全体の燃料消耗を計算した結果を読み取り、第一及び第二のＣＰＲマトリックスをそれぞれ求め、正確に装荷した状況の計算結果を読み取り、第三のＣＰＲマトリックスを求め、ｃ４の段階では第二のＣＰＲマトリックスを第四、第五及び第六のＣＰＲマトリックスと置き換えて、第一のＣＰＲマトリックスに重ね合わせてから第三のＣＰＲマトリックスを引くと、少なくとも一つの燃料集合体の置き違いの状況の組み合わせのＣＰＲマトリックスを得る。
【００１７】
請求項５に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法は、請求項１におけるステップｄでは少なくとも次の段階を実施し：ｄ１の段階ではステップｃで得た燃料集合体の置き違いの状況のそれぞれの組み合わせのＣＰＲマトリックスからそれぞれの周期の燃料消耗点のＭＣＰＲ値を読み取り、ｄ２の段階では燃料集合体を正確に装荷した状況のＣＰＲマトリックスからそれぞれの周期の燃料消耗点のＭＣＰＲ値を読み取り、ｄ３の段階では燃料集合体の置き違いのそれぞれの状況の組み合わせのＣＰＲマトリックスにおけるそれぞれの周期の燃料消耗点のＭＣＰＲ値から、燃料集合体を正確に装荷した状況のＣＰＲマトリックスから相対する周期の燃料消耗点のＭＣＰＲ値を引くと、それぞれの周期の燃料消耗点のＭＣＰＲの変化量（Ｄｉｆｆｅｒｅｎｃｅ ｏｆ Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＤＣＰＲ）を求め、ｄ４の段階ではそれぞれの周期の燃料消耗点のＤＣＰＲ値からＭａｘＤＣＰＲを求める。
【００１８】
請求項６に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法は、請求項１における炉心の中性子の分析で使用する解析コードがＳＩＭＵＬＡＴＥ−３コードである。
【００１９】
請求項７に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法は、請求項１における二つのプログラムはＦｏｒｔｒａｎ言語から形成する。
【００２０】
請求項８に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法は、請求項１における沸騰水型軽水炉は従来の沸騰水型軽水炉、及び最新の沸騰水型軽水炉である。
【００２１】
請求項９に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法は、請求項１におけるステップｅではＭａｘＤＣＰＲの数値の大きいものから１００番目までの置き違いの状況の組み合わせを選び出し、炉心全体の燃料消耗を再度計算する。
【発明の効果】
【００２２】
この発明の燃料集合体の置き違いの分析方法は、分析する対象数が少ないため、分析する時間を短縮することができ、実用的で炉心の安全性を確保することができるという効果を有する。
【発明を実施するための最良の形態】
【００２３】
この発明は、分析対象となる燃料集合体の置き違いの数が少ない燃料集合体の置き違いの事故の分析方法であって、燃料集合体の置き違いの個別の計算結果を相互に重ね合わせてそれぞれの置き違いの結果を求める沸騰水型軽水炉の燃料集合体の置き違いの事故の分析方法であって、該分析方法は少なくとも次ぎのステップにより実施し：ステップａではプログラムを利用して少なくとも単一の燃料集合体の置き違いの状況を計算し、ステップｂでは炉心の中性子の分析方法を利用して上述の燃料の置き違いのそれぞれの個別の状況に対して炉心全体の燃料消耗を計算し、ステップｃでは別のプログラムを介して燃料集合体の置き違いの組み合わせを発生させて、置き違いの組み合わせに対応する二つの個別の置き違いの限界出力比（Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＣＰＲ）のマトリックスを検出し、該二組のＣＰＲマトリックスを重ね合わせると、少なくとも一つの置き違いの組み合わせのＣＰＲマトリックスを得ることができ、ステップｄでは燃料集合体の置き違いのそれぞれの組み合わせのＣＰＲマトリックスから最小限界出力比（Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭＣＰＲ）の最大変化量（Ｍａｘｉｍｕｍ ｏｆ Ｄｉｆｆｅｒｅｎｃｅ ｏｆ Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭａｘＤＣＰＲ）を求め、ステップｅでは上述の置き違いの組み合わせをＭａｘＤＣＰＲに基づいて数値の大きいものから小さい順に並べ、大きいものから１０番目までの置き違いの組み合わせを選び出し、炉心全体の燃料消耗を再度計算して、該１０個の置き違いの組み合わせの中から最大のＭａｘＤＣＰＲ値を検出し、該数値が燃料集合体の置き違いの事故の最終的な分析結果になる。
係る構成の燃料集合体の置き違いの事故の分析方法について、その構造と特徴を詳述するために具体的な実施例を挙げ、以下に説明する。
【実施例１】
【００２４】
図１は、この発明のフローチャートを示したものである。図面によれば、この発明は燃料集合体の置き違いの事故の分析方法であって、二つの個別の燃料集合体の置き違いの状況の限界出力比（Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＣＰＲ）のマトリックスと、燃料集合体を正しく装荷した状況のＣＰＲマトリックスとを交互に重ね合わせてそれぞれの配置を誤った組み合わせのＣＰＲマトリックスを求めた。また、この発明は対象となる炉心の四分の一の領域にのみ分析を行い、分析は次のステップで行った。
【００２５】
ステップａ（１１）では燃料集合体の置き違いの個別の状況に関する数値を計算した。プログラムを利用して少なくとも一つの個別の配置を誤った状況に関する数値を計算した。該プログラムはＦｏｒｔｒａｎ言語により作成した。個別の置き違いの状況の数は対象になる炉心の燃料集合体の個数の４分の１に燃料集合体の形態の個数から１を引いた数値をかけた数になった。
【００２６】
ステップｂ（１２）ではそれぞれの個別の燃料集合体の置き違いに対して炉心の燃料の消耗を計算した。炉心全体の中性子の分析方法を利用して、上述のそれぞれの個別の燃料集合体の置き違いの状況対する炉心全体の燃料の消耗を計算した。該炉心の中性子の分析方法に使用する解析コードはＳＩＭＵＬＡＴＥ−３コードとした。
【００２７】
ステップｃ（１３）では複数のＣＰＲマトリックスを重ね合わせると燃料集合体の置き違いの状況のＣＰＲマトリックスを得ることができた。別のＦｏｒｔｒａｎプログラムにより燃料集合体の置き違いの状況の組み合わせを作成し、燃料集合体の置き違いの状況の組み合わせに対応する二つの個別の燃料集合体の置き違いのＣＰＲマトリックスを求め、該二組のＣＰＲマトリックスと燃料集合体を正しく装荷した状況のＣＰＲマトリックスとを重ね合わせると、少なくとも一つの燃料集合体の置き違いの組み合わせのＣＰＲマトリックスを得ることができた。
【００２８】
ステップｄ（１４）ではそれぞれの燃料集合体の置き違いの状況の組み合わせのＭａｘＤＣＰＲを求めた。それぞれの燃料集合体の置き違いの状況の組み合わせの最小限界出力比（Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭＣＰＲ）の最大変化量（Ｍａｘｉｍｕｍ ｏｆ Ｄｉｆｆｅｒｅｎｃｅ ｏｆ Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭａｘＤＣＰＲ）を求めた。
【００２９】
ステップｅ（１５）では燃料集合体の置き違いの事故の分析結果を得た。該燃料集合体の置き違いの組み合わせを排列し、最終的な計算を行うと、燃料集合体の置き違いの事故の分析の最終結果を得た。
【実施例２】
【００３０】
この発明の分析方法を初装荷炉心に応用した場合を説明した。
【００３１】
図１はこの発明のフローチャートを示した。図１Ａはステップａのフローチャートを示した。図１Ｂはステップｃのフローチャートを開示した。図１Ｃはステップｄのフローチャートを示した。図１Ｄは実例１の初装荷炉心における燃料集合体のレイアウトを示した。
【００３２】
図面によれば、実例１では初装荷炉心１６に対して燃料集合体の置き違いの事故における分析を行った。この実例で分析する初装荷炉心１６は、対象となる炉心の縦、横軸が直交して四等分された領域の四分の一とした。該初装荷炉心１６は６２４体の燃料集合体を具えてなり、第一燃料形態１６１、第二燃料形態１６２及び第三燃料形態１６３など三種類の燃料形態を具えたものとした。この発明の分析方法は炉心の四分の一の領域１６４に対してのみ分析を行い、分析には次のステップを実施した。
【００３３】
ステップａではＦｏｒｔｒａｎ言語で作成されたプログラムで３１２体の燃料集合体の置き違いの状況の数値を計算した。計算は次の段階で行った。
【００３４】
（ａ１）の段階において燃料集合体を正確に配置した燃料形態の配置データを読み取った１１１のステップを行い、（ａ２）の段階においてもともとの燃料形態と他の二種類の燃料形態とを置き換えて、それぞれの炉心のに二つの個別の誤った配置を形成した１１２のステップを行った。炉心の四分の一の領域１６４に、炉心のレイアウトに基づいてもともとの燃料形態とその他の二種類の燃料形態とを置き換えると、それぞれの炉心のレイアウトに二つの個別の置き違いの状況が発生した。前述の方法で計３１２体の個別の燃料集合体の置き違いの状況の数値を計算した。
【００３５】
上述の（ａ２）の段階で個別の置き違いの数が炉心の四分の一の領域１６４の燃料集合体の個数に燃料形態の個数から１を引いた数値をかけた数値であることが得られた。よって、この実施例では３１２個（６２４／４×（３−１）＝３１２）の個別の燃料集合体の置き違いを具えたものとなった。
【００３６】
ステップｂではＳＩＭＵＬＡＴＥ−３コードの炉心の中性子を分析する方法で、上述の３１２体の置き違いのそれぞれの状況に対して炉心全体の燃料の消耗を計算した。
【００３７】
ステップｃではＦｏｒｔｒａｎ言語で作成したプログラムで燃料集合体の置き違いの状況の組み合わせを発生させて、３１２体の燃料集合体の置き違いの状況から計算した炉心全体の燃料の消耗の結果から相対する二つの個別の置き違いのＣＰＲマトリックスを求めた。また、該二組のＣＰＲマトリックスと正確に装荷した状況のＣＰＲマトリックスとを重ね合わせると、少なくとも一つの燃料集合体の置き違いの状況の組み合わせのＣＰＲマトリックスを得ることができ、マトリックスを得るために次の段階を実施した。
【００３８】
（ｃ１）の段階では燃料集合体を正確に配置した燃料形態の配置データを読み取った１３１のステップを実施した。
【００３９】
（ｃ２）の段階では燃料形態の情報を得た後、さらにレイアウトの順序に従って別の燃料形態を得て、置き違いの組み合わせを発生させ、炉心の四分の一の領域１６４における炉心のレイアウトの順序に基づいて燃料形態の情報を読み取り、炉心全体のレイアウトの順序に基づいて別の燃料形態の情報を読み取り置き違いの組み合わせを発生させた１３２のステップを実施した。
【００４０】
（ｃ３）の段階では第一及び第二ＣＰＲマトリックス及び正確に配置した状況の計算結果から第三のＣＰＲマトリックスを得た１３３のステップを実施した。誤った配置の組み合わせに対応する二つの個別の燃料集合体の置き違いの状況の計算結果を読み取り、第一及び第二のＣＰＲマトリックスを求めた。また、正確に配置した計算結果を読み取り第三のＣＰＲマトリックスを得た。
【００４１】
（ｃ４）では第二のＣＰＲマトリックスを置き換えて第四、第五及び第六のＣＰＲマトリックスを得たステップ１３４を実施した。また、第一のＣＰＲマトリックスを重ね合わせてから第三のＣＰＲマトリックスを引くと誤った配置の組み合わせのＣＰＲマトリックスを得ることができた。燃料集合体を置き違えた炉心の第二のＣＰＲマトリックスに対応するマトリックスを置き換え、第四、第五及び第六のＣＰＲマトリックス（追って説明する）を取得し、それぞれ第一のＣＰＲマトリックスに重ね合わせて、第三のＣＰＲマトリックスを引くと、少なくとも一つの誤った配置の組み合わせのＣＰＲマトリックスを得られた。
【００４２】
ステップｄではそれぞれの誤った配置の組み合わせのＭａｘＤＣＰＲ値を求めた。検出には次のステップを実施した。
【００４３】
（ｄ１）では、燃料集合体の置き違いのそれぞれの組み合わせのＭＣＰＲを読み取った１４１のステップを実施した。ステップｃで得た燃料集合体の置き違いのそれぞれの組み合わせのＣＰＲマトリックスから、各周期の燃料の消耗点のＭＣＰＲを読み取った（この実施例では計２７３３２個の置き違いの組み合わせを具えていたため、２７３３２個のＣＰＲマトリックスを摘出する必要があった。また全周期では４５個の燃料消耗点を具えていた）。
【００４４】
（ｄ２）では燃料集合体を正確に配置した状況のＭＣＰＲを得た１４２のステップを実施した。燃料集合体を正確に配置した状況のＣＰＲマトリックスからそれぞれの周期の燃料の消耗点におけるＭＣＰＲを読み取った。
【００４５】
（ｄ３）では燃料集合体の置き違いの状況の組み合わせのそれぞれのＭＣＰＲ値から燃料集合体を正確に装荷した状況のＭＣＰＲ値を引くと、燃料集合体の置き違いの状況の組み合わせのＤＣＰＲを得られた１４３のステップを実施した。（ｄ２）の段階のＭＣＰＲ値から（ｄ１）の段階のＭＣＰＲ値を引くと、燃料集合体の置き違いのそれぞれの状況の組み合わせの各周期の燃料の消耗点のＭＣＰＲの変化量を得ることができた（ＤＣＰＲ；Ｄｉｆｆｅｒｅｎｃｅ ｏｆ Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ）。
【００４６】
（ｄ４）の段階ではＤＣＰＲ値からＭａｘＤＣＰＲを求めた１４４のステップを実施した。（ｄ３）の段階のＤＣＰＲ値から最大値を求めると、燃料集合体の置き違いのそれぞれの状況の組み合わせの各周期の燃料の消耗点におけるＭａｘＤＣＰＲを得ることができた。
【００４７】
ステップｅでは上述の燃料集合体の置き違いの組み合わせの状況であって、そのＭａｘＤＣＰＲを大きな値から小さな値に並べ、最大値から１０番目までの数値か、もしくは１００番目までの誤った配置の組み合わせを選び出して、再度、炉心全体の燃料の消耗を計算した。また、該１０個、もしくは１００個の燃料集合体の置き違いの状況の組み合わせにおいて最大のＭａｘＤＣＰＲ値を求めた。該数値は燃料集合体の置き違いの事故の最終的な分析結果になった。
【００４８】
上述の方法は、燃料集合体の置き違いの状況を検出し、求められた燃料集合体の置き違いの状況に対して分析を行うため、軽水炉の運転前に実施する燃料集合体の置き違いの事故の安全性を分析する時間を縮小することができた。
【実施例３】
【００４９】
この発明の分析方法を再装荷した炉心に応用する場合を説明した。
【００５０】
図１Ｅは、再装荷した炉心の燃料集合体に応用した場合を示した。図面によれば、この発明の実例２は再装荷した炉心１７で燃料集合体を誤って配置した場合の事故を分析する方法とした。この実例では再装荷した炉心１７の四分の一の領域に対して分析を実施した。該再装荷した炉心１７は６２４体の燃料集合体を具えてなり、第一の燃料形態８１、第二の燃料形態８２及び第四の燃料形態８３など三種類の燃料形態を具えてなるようにした。
【００５１】
該第一の燃料形態８１及び第二の燃料形態８２は、第一次に消耗した燃料とした。よって、第一の燃料形態８１及び第二の燃料形態８２は燃料消耗の分布を有し、同一形態の燃料集合体はその燃料の消耗に差があるため、２体の同一形態の燃料集合体でも置き違いの状況の組み合わせが発生する可能性があった。これは初装荷炉心ではみられない状況であった。この実例の燃料形態は三種類とした。しかしながら、第一の燃料形態に見られる燃料消耗の分布は７．２５７〜１３．１９５ＧＷＤ／ＭＴ（ｇｉｇａｗａｔｔ−ｄａｙｓ ｐｅｒ ｍｅｔｒｉｃ ｔｏｎ）であって、燃料消耗の分布の範囲が比較的大きかった。よって第一の燃料形態は更に二種類の形態に分けることができた。一種類目の形態は燃料消耗の分布が１０ＧＷＤ／ＭＴよりも大きく、もう一種類の形態は燃料消耗の分布が１０ＧＷＤ／ＭＴよりも小さくなった。また、第二種の燃料形態の燃料消耗の分布は１０．８３１〜１１．３８７ＧＷＤ／ＭＴになり、分布の範囲が大きくなかった。このため更に細かく燃料形態を分ける必要がなく、この実施例では計４種類の燃料形態を具えていた。また、この発明の分析方法は、炉心の四分の一の領域１７４に対してのみ分析を行い、分析には下記のステップを実施した。
【００５２】
ステップａではＦｏｒｔｒａｎ言語で作成したプログラムで４６８体の個別の燃料集合体の配置を誤った状況を計算した。計算には次ぎの段階を実施した。
【００５３】
（ａ１）の段階では、燃料集合体を正確に配置した燃料形態のデータを読み取った。
【００５４】
（ａ２）の段階では炉心の四分の一の領域１７４において炉心のレイアウトの順序に基づいてもともとの燃料形態に分けて、その他の三種類の燃料形態に置き換え、それぞれの炉心のレイアウトに対して燃料集合体の置き違いの三つの状況を形成した。この方法により計４６８体の個別の置き違いの状況を計算した。
【００５５】
上述の（ａ２）の段階は個別の置き違いの状況の数値が炉心の四分の一の領域１７４の燃料集合体の数に燃料形態の数から１を引いた数をかけたものになった。よって、この実施例では４６８個（６２４／４×（４−１）＝４６８）の置き違いの状況が得られた。
【００５６】
ステップｂではＳＩＭＵＬＡＴＥ−３コードを炉心の中性子を分析する方法として、上述の４６８体の個別の置き違いのレイアウトにおける炉心全体の燃料の消耗を計算した。
【００５７】
ステップｃでは別のＦｏｒｔｒａｎ言語から作成したプログラムで燃料集合体の置き違いの組み合わせを発生させた。４６８体の燃料集合体の個別の置き違いの炉心全体の燃料消耗を計算した結果から相対する二つの個別の置き違いのＣＰＲマトリックスを求める。該二組のＣＰＲマトリックスと燃料集合体を正確に装填した状況のＣＰＲマトリックスを重ね合わせると、少なくとも一つの燃料集合体の置き違いの組み合わせのＣＰＲマトリックスが得られた。マトリックスを得るために次ぎの段階を実施した。
【００５８】
（ｃ１）の段階では燃料集合体を正確に配置した燃料形態の配置データを読み取った。
【００５９】
（ｃ２）の段階では炉心の四分の一の領域１７４において炉心のレイアウトの順序に基づいて一種類の燃料形態の情報を得た。その後、炉心全体の配置順序に基づいてもう一種類の燃料形態を得て燃料集合体を置き違えた状況の組み合わせを形成した。
【００６０】
（ｃ３）の段階では燃料集合体を置き違えた状況の組み合わせに対応する二つの個別の燃料集合体の置き違えた炉心全体の燃料の消耗を計算した結果を読み取り、第一及び第二のＣＰＲマトリックスをそれぞれ求めた。また、燃料集合体を正確に装荷した状況の炉心全体の燃料の消耗を計算した結果を読み取ると、第三のＣＰＲマトリックスを得られた。
【００６１】
（ｃ４）の段階では燃料集合体を置き違えた状況の炉心の位置により第二のＣＰＲマトリックスと対応するマトリックスのレイアウトの交換を行った。レイアウトの交換は第四、第五及び第六のＣＰＲマトリックス（追って説明する）とし、それぞれ第一のＣＰＲマトリックスと重ね合わせた後、第三のＣＰＲマトリックスを引くと、少なくとも一つの燃料集合体を置き違えた状況の組み合わせのＣＰＲマトリックスを得られた。
【００６２】
ステップｄではそれぞれの誤った配置の組み合わせのＭａｘＤＣＰＲ値を求めた。検出には、次ぎのステップを実施した。
【００６３】
（ｄ１）の段階ではステップｃで得たそれぞれの誤った配置の組み合わせのＣＰＲマトリックスから、各周期の燃料消耗点のＭＣＰＲ値を求めた（この実施例では３２１１６個の誤った配置の組み合わせを具えていたため、３２１１６個のＣＰＲマトリックスを求めた。また、全周期では計３５個の燃料消耗点を具えていた）。
【００６４】
（ｄ２）の段階では燃料集合体を正確に装荷した状況のＣＰＲマトリックスから、それぞれの周期の燃料消耗点のＭＣＰＲ値を読み取った。
【００６５】
（ｄ３）の段階では（ｄ２）の段階のＭＣＰＲ値からステップｄ１のＭＣＰＲ値を引くと、それぞれの置き違いの組み合わせの各周期の燃料消耗点のＤＣＰＲ値を得た。
【００６６】
（ｄ４）の段階では（ｄ３）の段階のＤＣＰＲ値から最大値を選び出すと、それぞれの誤った配置の組み合わせの各周期の燃料消耗点のＭａｘＤＣＰＲ値を得られた。
【００６７】
ステップｅでは上述の誤った配置の組み合わせのＭａｘＤＣＰＲを大きい数値から小さくなるように並べた。最も大きい数値から１０番目か、もしくは１００番目までの誤った組み合わせを選び出し、全ての炉心の燃料消耗を再度計算し、１０個、もしくは１００個の誤った組み合わせの中から最大のＭａｘＤＣＰＲ値を求めた。その数値は燃料集合体の設置を誤った事故の最終的な分析結果になった。
【００６８】
発生する可能性のある燃料集合体の置き違いの状況を検出し、燃料集合体を置き違えた状況に対して上述の方法で分析を行うと、軽水炉の運転前に実施する燃料集合体の配置の誤りによる事故の安全性に関する分析の時間を短縮することができた。
【００６９】
図２Ａは燃料集合体を正確に装荷した状況の燃料の配置を示した。図２Ｂは燃料集合体を置き違えた第一の状況を示した。図２Ｃは燃料集合体を置き違えた第一の個別の状況を示した。図２Ｄは燃料集合体を置き違えた第二の個別の状況を示した。図２Ｅは燃料集合体を置き違えた第三の個別の状況を示した。図２Ｆは燃料集合体を置き違えた第二の組み合わせを示した。図２Ｇは燃料集合体を置き違えた第四の個別の状況を示した。図２Ｈは燃料集合体を置き違えた第三の状況を示した。図２Ｉは燃料集合体を置き違えた第五の個別の状況を示した。図２Ｊは燃料集合体を置き違えた第四の状況を示した。
【００７０】
図面によれば、ステップｃにおいて複数のＣＰＲマトリックスを重ね合わせた方法を下記に示した。図２Ａでは、炉心の四分の一の対称となる領域２１に燃料集合体を正確に装荷したと仮定した。また、そのＣＰＲマトリックスは第三のＣＰＲマトリックス（ｆｏｌ）であって、次のように表示された。
【数１】
【００７１】
図２Ｂでは燃料集合体の置き違いの第一の状況１０を仮定した。また、図２Ｃは上述の燃料集合体の置き違いの状況の第一の組み合わせ２２１の第一の個別の誤ったレイアウト２３１であって、そのＣＰＲマトリックスは第一のＣＰＲマトリックス（Ｐ）だった。該ＣＰＲマトリックスを下記に表示した。
【数２】
【００７２】
図２Ｄでは上述の燃料集合体の第一の置き違いの状況の組み合わせ２２１の第二の個別の置き違いの状況２３２を示した。該ＣＰＲマトリックスは第二のＣＰＲマトリックス（Ｑ）であって、下記に表示した。
【数３】
【００７３】
よって、上述の第一のＣＰＲマトリックス（Ｐ）と第二のＣＰＲマトリックス（Ｑ）を加えて、第三のＣＰＲマトリックス（ｆｏ１）を引くと、第一の誤った配置の組み合わせのＣＰＲマトリックスを得ることができた。該ＣＰＲマトリックスを下記に表示した。
【数４】
【００７４】
第二のＣＰＲマトリックスはＹ軸を通して反射すると、第２Ｅに開示した第三の個別の置き違いの状況２３３のＣＰＲマトリックスを得ることができた。該ＣＰＲマトリックスは第四のＣＰＲマトリックスであって、第一のＣＰＲマトリックスに第四のＣＰＲマトリックスを加えた後、第三のＣＰＲマトリックスを引くと、第二の置き違いの状況の組み合わせ２２２（図２Ｆ参照）のＣＰＲマトリックスを得ることができた。該ＣＰＲマトリックスを下記に表示した。
【数５】
【００７５】
第二のＣＰＲマトリックスがＸ軸を通して反射すると、図２Ｇに開示する第四の個別の置き違いの状況２３４のＣＰＲマトリックスを得られた。該ＣＰＲマトリックスは第五のＣＰＲマトリックスであって、第一のＣＰＲマトリックスに第五のＣＰＲマトリックスを加えた後、第三のＣＰＲマトリックスを引くと、第三の置き違いの状況の組み合わせ２２３（図２Ｈ参照）のＣＰＲマトリックスを得られた。該マトリックスを下記に表示した。
【数６】
【００７６】
第二のＣＰＲマトリックスがＸ軸及びＹ軸を通して２度反射されると、図２Ｉに開示する第五の個別の置き違いの状況２３５のＣＰＲマトリックスを得ることができた。該ＣＰＲマトリックスは第六のＣＰＲマトリックスであって、第一のＣＰＲマトリックスに第六のＣＰＲマトリックスを加えて、第三のＣＰＲマトリックスを引くと、第四の置き違いの状況の組み合わせ２２４（図２Ｊ参照）のＣＰＲマトリックスを得られた。該マトリックスを下記に表示した。
【数７】
【００７７】
上述の方法を通して初装荷炉心を分析する場合、３１２個の単一置き違いのＣＰＲマトリックスと燃料集合体を正確に装荷したＣＰＲマトリックスとを重ね合わせた後、２７３３２個の置き違いの状況のＣＰＲマトリックスを得た。再装荷した炉心を分析する場合には、４６８個の個別の誤った配置のＣＰＲマトリックスを正確に装荷したＣＰＲマトリックスを重ね合わせて３２１１６個の置き違いの状況の組み合わせのＣＰＲマトリックスを得た。
【００７８】
上述のように、この発明の沸騰水型軽水炉における燃料集合体の置き違いの事故の分析方法は、従来の分析方法の欠点を有効に改善し、分析に必要な誤った組み合わせの数を減らし、運転前に実施する燃料集合体の置き違いの事故の安全性を分析する時間を短縮して、炉心の安全を確保することができた。
【００７９】
以上は、この発明の好ましい実施例であって、この発明の実施の範囲を限定するものではない。よって、当業者のなし得る修正、もしくは変更であって、この発明の精神の下においてなされ、この発明に対して均等の効果を有するものは、いずれもこの発明の特許請求の範囲に属するものとする。
【図面の簡単な説明】
【００８０】
【図１】この発明の沸騰水型軽水炉の燃料集合体を誤って配置した事故の分析方法を示したフローチャートである。
【図１Ａ】図１に開示する沸騰水型軽水炉の燃料集合体を誤って配置した事故の分析方法のステップａを示したフローチャートである。
【図１Ｂ】図１に開示する沸騰水型軽水炉の燃料集合体を誤って配置した事故の分析方法のステップｃを示したフローチャートである。
【図１Ｃ】図１に開示する沸騰水型軽水炉の燃料集合体を誤って配置した事故の分析方法のステップｄを示したフローチャートである。
【図１Ｄ】この発明の沸騰水型軽水炉の燃料集合体の実例１の初装荷炉心を示した断面図である。
【図１Ｅ】この発明の沸騰水型軽水炉の燃料集合体の実例２の再装荷した炉心を示した断面図である。
【図２Ａ】この発明の沸騰水型軽水炉に燃料集合体を正確に装荷した状況を示した断面図である。
【図２Ｂ】この発明の燃料集合体を誤って配置した第一の状況を示した断面図である。
【図２Ｃ】この発明の燃料集合体を誤って配置した第一の個別の状況を示した断面図である。
【図２Ｄ】この発明の燃料集合体を誤って配置した第二の個別の状況を示した断面図である。
【図２Ｅ】この発明の燃料集合体を誤って配置した第三の個別の状況を示した断面図である。
【図２Ｆ】この発明の燃料集合体を誤って配置した第二の状況を示した断面図である。
【図２Ｇ】この発明の燃料集合体を誤って配置した第四の個別の状況を示した断面図である。
【図２Ｈ】この発明の燃料集合体を誤って配置した第三の状況を示した断面図である。
【図２Ｉ】この発明の燃料集合体を誤って配置した第五の個別の状況を示した断面図である。
【図２Ｊ】この発明の燃料集合体を誤って配置した第四の状況を示した断面図である。
【符号の説明】
【００８１】
１６ 初装荷炉心
１７ 再装荷した炉心
１６１ 第一燃料形態
１６２ 第二燃料形態
１６３ 第三燃料形態
１６４ 炉心の四分の一の領域
１７１ 第一燃料形態
１７２ 第二燃料形態
１７３ 第四燃料形態
１７４ 炉心の四分の一の領域
２１ 炉心の四分の一の対称となる領域
２２１ 第一の置き違いの状況の組み合わせ
２２２ 第二の置き違いの状況の組み合わせ
２２３ 第三の置き違いの状況の組み合わせ
２２４ 第四の置き違いの状況の組み合わせ
２３２ 第二の個別の置き違いの状況
２３３ 第三の個別の置き違いの状況
２３４ 第四の個別の置き違いの状況
２３５ 第五の個別の置き違いの状況
【特許請求の範囲】
【請求項１】
燃料集合体の置き違いの個別の計算結果を相互に重ね合わせてそれぞれの置き違いの結果を求める沸騰水型軽水炉の燃料集合体の置き違いの事故の分析方法であって、
該分析方法は少なくとも次ぎのステップにより実施し：
ステップａではプログラムを利用して少なくとも単一の燃料集合体の置き違いの状況を計算し、
ステップｂでは炉心の中性子の分析方法を利用して上述の燃料の置き違いのそれぞれの個別の状況に対して炉心全体の燃料消耗を計算し、
ステップｃでは別のプログラムを介して燃料集合体の置き違いの組み合わせを発生させて、置き違いの組み合わせに対応する二つの個別の置き違いの限界出力比（Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＣＰＲ）のマトリックスを検出し、該二組のＣＰＲマトリックスを重ね合わせると、少なくとも一つの置き違いの組み合わせのＣＰＲマトリックスを得ることができ、
ステップｄでは燃料集合体の置き違いのそれぞれの組み合わせのＣＰＲマトリックスから最小限界出力比（Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭＣＰＲ）の最大変化量（Ｍａｘｉｍｕｍ ｏｆ Ｄｉｆｆｅｒｅｎｃｅ ｏｆ Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭａｘＤＣＰＲ）を求め、
ステップｅでは上述の置き違いの組み合わせをＭａｘＤＣＰＲに基づいて数値の大きいものから小さい順に並べ、大きいものから１０番目までの置き違いの組み合わせを選び出し、炉心全体の燃料消耗を再度計算して、該１０個の置き違いの組み合わせの中から最大のＭａｘＤＣＰＲ値を検出し、該数値が燃料集合体の置き違いの事故の最終的な分析結果になることを特徴とする沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法。
【請求項２】
前記ステップａの個別の置き違いの状況の数値が、燃料集合体の個数に燃料形態の個数から１を引いた数をかけた数値になることを特徴とする請求項１に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法。
【請求項３】
前記ステップａの個別の置き違いの状況の計算には次の段階を実施し：
ａ１の段階では燃料集合体を正確に配置した状況の燃料形態の配置データを読み取り、
ａ２の段階では炉心の四分の一の領域において、炉心のレイアウトの順序に基づいてそれぞれもともとの燃料形態とその他の燃料形態に分けてから置き換えて、それぞれの個別の置き違いの状況を計算することを特徴とする請求項１に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法。
【請求項４】
前記ステップｃにおける燃料集合体の置き違いの組み合わせのＣＰＲマトリックスの形成には次の段階を実施し：
ｃ１の段階では燃料集合体を正確に配置した状況の燃料集合体のレイアウトのデータを読み取り、
ｃ２の段階では炉心の四分の一の領域において炉心のレイアウトに基づいて燃料形態を得て、再度、炉心全体のレイアウトに基づいて別の燃料形態を得て、置き違いの組み合わせを発生させて、
ｃ３の段階では燃料集合体の置き違いの組み合わせに相対する第一及び第二の置き違いの組み合わせの炉心全体の燃料消耗を計算した結果を読み取り、第一及び第二のＣＰＲマトリックスをそれぞれ求め、正確に装荷した状況の計算結果を読み取り、第三のＣＰＲマトリックスを求め、
ｃ４の段階では第二のＣＰＲマトリックスを第四、第五及び第六のＣＰＲマトリックスと置き換えて、第一のＣＰＲマトリックスに重ね合わせてから第三のＣＰＲマトリックスを引くと、少なくとも一つの燃料集合体の置き違いの状況の組み合わせのＣＰＲマトリックスを得ることができることを特徴とする請求項１に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法。
【請求項５】
前記ステップｄでは少なくとも次の段階を実施し：
ｄ１の段階ではステップｃで得た燃料集合体の置き違いの状況のそれぞれの組み合わせのＣＰＲマトリックスからそれぞれの周期の燃料消耗点のＭＣＰＲ値を読み取り、
ｄ２の段階では燃料集合体を正確に装荷した状況のＣＰＲマトリックスからそれぞれの周期の燃料消耗点のＭＣＰＲ値を読み取り、
ｄ３の段階では燃料集合体の置き違いのそれぞれの状況の組み合わせのＣＰＲマトリックスにおけるそれぞれの周期の燃料消耗点のＭＣＰＲ値から、燃料集合体を正確に装荷した状況のＣＰＲマトリックスから相対する周期の燃料消耗点のＭＣＰＲ値を引くと、それぞれの周期の燃料消耗点のＭＣＰＲの変化量（Ｄｉｆｆｅｒｅｎｃｅ ｏｆ Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＤＣＰＲ）を求め、
ｄ４の段階ではそれぞれの周期の燃料消耗点のＤＣＰＲ値からＭａｘＤＣＰＲを求めることを特徴とする請求項１に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法。
【請求項６】
前記炉心の中性子の分析で使用する解析コードはＳＩＭＵＬＡＴＥ−３コードであることを特徴とする請求項１に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法。
【請求項７】
前記二つのプログラムはＦｏｒｔｒａｎ言語から形成することを特徴とする請求項１に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法。
【請求項８】
前記沸騰水型軽水炉は従来の沸騰水型軽水炉、及び最新の沸騰水型軽水炉であることを特徴とする請求項１に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法。
【請求項９】
前記ステップｅではＭａｘＤＣＰＲの数値の大きいものから１００番目までの置き違いの状況の組み合わせを選び出し、炉心全体の燃料消耗を再度計算することを特徴とする請求項１に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法。
【請求項１】
燃料集合体の置き違いの個別の計算結果を相互に重ね合わせてそれぞれの置き違いの結果を求める沸騰水型軽水炉の燃料集合体の置き違いの事故の分析方法であって、
該分析方法は少なくとも次ぎのステップにより実施し：
ステップａではプログラムを利用して少なくとも単一の燃料集合体の置き違いの状況を計算し、
ステップｂでは炉心の中性子の分析方法を利用して上述の燃料の置き違いのそれぞれの個別の状況に対して炉心全体の燃料消耗を計算し、
ステップｃでは別のプログラムを介して燃料集合体の置き違いの組み合わせを発生させて、置き違いの組み合わせに対応する二つの個別の置き違いの限界出力比（Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＣＰＲ）のマトリックスを検出し、該二組のＣＰＲマトリックスを重ね合わせると、少なくとも一つの置き違いの組み合わせのＣＰＲマトリックスを得ることができ、
ステップｄでは燃料集合体の置き違いのそれぞれの組み合わせのＣＰＲマトリックスから最小限界出力比（Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭＣＰＲ）の最大変化量（Ｍａｘｉｍｕｍ ｏｆ Ｄｉｆｆｅｒｅｎｃｅ ｏｆ Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＭａｘＤＣＰＲ）を求め、
ステップｅでは上述の置き違いの組み合わせをＭａｘＤＣＰＲに基づいて数値の大きいものから小さい順に並べ、大きいものから１０番目までの置き違いの組み合わせを選び出し、炉心全体の燃料消耗を再度計算して、該１０個の置き違いの組み合わせの中から最大のＭａｘＤＣＰＲ値を検出し、該数値が燃料集合体の置き違いの事故の最終的な分析結果になることを特徴とする沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法。
【請求項２】
前記ステップａの個別の置き違いの状況の数値が、燃料集合体の個数に燃料形態の個数から１を引いた数をかけた数値になることを特徴とする請求項１に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法。
【請求項３】
前記ステップａの個別の置き違いの状況の計算には次の段階を実施し：
ａ１の段階では燃料集合体を正確に配置した状況の燃料形態の配置データを読み取り、
ａ２の段階では炉心の四分の一の領域において、炉心のレイアウトの順序に基づいてそれぞれもともとの燃料形態とその他の燃料形態に分けてから置き換えて、それぞれの個別の置き違いの状況を計算することを特徴とする請求項１に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法。
【請求項４】
前記ステップｃにおける燃料集合体の置き違いの組み合わせのＣＰＲマトリックスの形成には次の段階を実施し：
ｃ１の段階では燃料集合体を正確に配置した状況の燃料集合体のレイアウトのデータを読み取り、
ｃ２の段階では炉心の四分の一の領域において炉心のレイアウトに基づいて燃料形態を得て、再度、炉心全体のレイアウトに基づいて別の燃料形態を得て、置き違いの組み合わせを発生させて、
ｃ３の段階では燃料集合体の置き違いの組み合わせに相対する第一及び第二の置き違いの組み合わせの炉心全体の燃料消耗を計算した結果を読み取り、第一及び第二のＣＰＲマトリックスをそれぞれ求め、正確に装荷した状況の計算結果を読み取り、第三のＣＰＲマトリックスを求め、
ｃ４の段階では第二のＣＰＲマトリックスを第四、第五及び第六のＣＰＲマトリックスと置き換えて、第一のＣＰＲマトリックスに重ね合わせてから第三のＣＰＲマトリックスを引くと、少なくとも一つの燃料集合体の置き違いの状況の組み合わせのＣＰＲマトリックスを得ることができることを特徴とする請求項１に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法。
【請求項５】
前記ステップｄでは少なくとも次の段階を実施し：
ｄ１の段階ではステップｃで得た燃料集合体の置き違いの状況のそれぞれの組み合わせのＣＰＲマトリックスからそれぞれの周期の燃料消耗点のＭＣＰＲ値を読み取り、
ｄ２の段階では燃料集合体を正確に装荷した状況のＣＰＲマトリックスからそれぞれの周期の燃料消耗点のＭＣＰＲ値を読み取り、
ｄ３の段階では燃料集合体の置き違いのそれぞれの状況の組み合わせのＣＰＲマトリックスにおけるそれぞれの周期の燃料消耗点のＭＣＰＲ値から、燃料集合体を正確に装荷した状況のＣＰＲマトリックスから相対する周期の燃料消耗点のＭＣＰＲ値を引くと、それぞれの周期の燃料消耗点のＭＣＰＲの変化量（Ｄｉｆｆｅｒｅｎｃｅ ｏｆ Ｍｉｎｉｍｕｍ Ｃｒｉｔｉｃａｌ Ｐｏｗｅｒ Ｒａｔｉｏ、ＤＣＰＲ）を求め、
ｄ４の段階ではそれぞれの周期の燃料消耗点のＤＣＰＲ値からＭａｘＤＣＰＲを求めることを特徴とする請求項１に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法。
【請求項６】
前記炉心の中性子の分析で使用する解析コードはＳＩＭＵＬＡＴＥ−３コードであることを特徴とする請求項１に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法。
【請求項７】
前記二つのプログラムはＦｏｒｔｒａｎ言語から形成することを特徴とする請求項１に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法。
【請求項８】
前記沸騰水型軽水炉は従来の沸騰水型軽水炉、及び最新の沸騰水型軽水炉であることを特徴とする請求項１に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法。
【請求項９】
前記ステップｅではＭａｘＤＣＰＲの数値の大きいものから１００番目までの置き違いの状況の組み合わせを選び出し、炉心全体の燃料消耗を再度計算することを特徴とする請求項１に記載の沸騰水型軽水炉の燃料集合体の置き違い事故の分析方法。
【図１】
【図１Ａ】
【図１Ｂ】
【図１Ｃ】
【図１Ｄ】
【図１Ｅ】
【図２Ａ】
【図２Ｂ】
【図２Ｃ】
【図２Ｄ】
【図２Ｅ】
【図２Ｆ】
【図２Ｇ】
【図２Ｈ】
【図２Ｉ】
【図２Ｊ】
【図１Ａ】
【図１Ｂ】
【図１Ｃ】
【図１Ｄ】
【図１Ｅ】
【図２Ａ】
【図２Ｂ】
【図２Ｃ】
【図２Ｄ】
【図２Ｅ】
【図２Ｆ】
【図２Ｇ】
【図２Ｈ】
【図２Ｉ】
【図２Ｊ】
【公開番号】特開２００７−２５６２５２（Ｐ２００７−２５６２５２Ａ）
【公開日】平成１９年１０月４日（２００７．１０．４）
【国際特許分類】
【出願番号】特願２００６−２０１０８６（Ｐ２００６−２０１０８６）
【出願日】平成１８年７月２４日（２００６．７．２４）
【出願人】（５９５１６５６５６）行政院原子能委員会核能研究所 (51)
【Ｆターム（参考）】
【公開日】平成１９年１０月４日（２００７．１０．４）
【国際特許分類】
【出願日】平成１８年７月２４日（２００６．７．２４）
【出願人】（５９５１６５６５６）行政院原子能委員会核能研究所 (51)
【Ｆターム（参考）】
[ Back to top ]