特開2023-175110線源強度の評価方法、遮蔽解析方法、線源強度の評価装置及びプログラム
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2023175110
(43)【公開日】2023-12-12
(54)【発明の名称】線源強度の評価方法、遮蔽解析方法、線源強度の評価装置及びプログラム
(51)【国際特許分類】
G21C 17/06 20060101AFI20231205BHJP
G21C 17/00 20060101ALI20231205BHJP
【ＦＩ】
G21C17/06 080
G21C17/00 500
【審査請求】未請求
【請求項の数】12
【出願形態】ＯＬ
(21)【出願番号】P 2022087391
(22)【出願日】2022-05-30
(71)【出願人】
【識別番号】000006208
【氏名又は名称】三菱重工業株式会社
(74)【代理人】
【識別番号】100149548
【弁理士】
【氏名又は名称】松沼 泰史
(74)【代理人】
【識別番号】100162868
【弁理士】
【氏名又は名称】伊藤 英輔
(74)【代理人】
【識別番号】100161702
【弁理士】
【氏名又は名称】橋本 宏之
(74)【代理人】
【識別番号】100189348
【弁理士】
【氏名又は名称】古都 智
(74)【代理人】
【識別番号】100196689
【弁理士】
【氏名又は名称】鎌田 康一郎
(72)【発明者】
【氏名】上山 正彦
(72)【発明者】
【氏名】山路 和也
(72)【発明者】
【氏名】小池 啓基
(72)【発明者】
【氏名】浅野 耕司
【テーマコード（参考）】
2G075
【Ｆターム（参考）】
2G075CA38
2G075CA48
2G075DA08
2G075FB04
2G075FB05
2G075GA18
(57)【要約】
【課題】使用済み燃料の線源強度の不確かさを含めて評価する方法を提供する。
【解決手段】線源強度の評価方法は、核種の壊変や生成を計算する計算コードに入力する核データの値を変動させたときに計算コードが算出する核種生成量に基づいて、核データ起因の線源強度の不確かさを評価するステップと、核種生成量の実測値と計算コードの計算値の差に基づく線源強度の不確かさと、核データ起因の線源強度の不確かさとに基づいて、核データ起因の線源強度の不確かさを校正することにより、計算コードによる線源強度の不確かさを評価するステップと、炉心の燃料のパラメータが記録された炉心設計データベースに基づいて、燃料の条件に基づく線源強度の不確かさを評価するステップと、所定の基準となる線源強度を、計算コードによる線源強度の不確かさと、燃料の条件に基づく線源強度の不確かさとによって補正した補正後の線源強度を算出するステップと、を有する。
【選択図】図１
【解決手段】線源強度の評価方法は、核種の壊変や生成を計算する計算コードに入力する核データの値を変動させたときに計算コードが算出する核種生成量に基づいて、核データ起因の線源強度の不確かさを評価するステップと、核種生成量の実測値と計算コードの計算値の差に基づく線源強度の不確かさと、核データ起因の線源強度の不確かさとに基づいて、核データ起因の線源強度の不確かさを校正することにより、計算コードによる線源強度の不確かさを評価するステップと、炉心の燃料のパラメータが記録された炉心設計データベースに基づいて、燃料の条件に基づく線源強度の不確かさを評価するステップと、所定の基準となる線源強度を、計算コードによる線源強度の不確かさと、燃料の条件に基づく線源強度の不確かさとによって補正した補正後の線源強度を算出するステップと、を有する。
【選択図】図１
【特許請求の範囲】
【請求項1】
核種の壊変や生成を計算する計算コードに入力する核データの値を変動させたときに前記計算コードが算出する前記核種の生成量に基づいて、前記核データ起因の線源強度の不確かさを評価するステップと、
前記核種の生成量の実測値と前記計算コードの計算値の差に基づく線源強度の不確かさと、前記核データ起因の線源強度の不確かさとに基づいて、前記核データ起因の線源強度の不確かさを校正することにより、前記計算コードによる線源強度の不確かさを評価するステップと、
炉心の燃料のパラメータが記録された炉心設計データベースに基づいて、前記燃料の条件に基づく線源強度の不確かさを評価するステップと、
所定の基準となる線源強度を、前記計算コードによる線源強度の不確かさと、前記燃料の条件に基づく線源強度の不確かさとによって補正した補正後の線源強度を算出するステップと、
を有する線源強度の評価方法。
前記核種の生成量の実測値と前記計算コードの計算値の差に基づく線源強度の不確かさと、前記核データ起因の線源強度の不確かさとに基づいて、前記核データ起因の線源強度の不確かさを校正することにより、前記計算コードによる線源強度の不確かさを評価するステップと、
炉心の燃料のパラメータが記録された炉心設計データベースに基づいて、前記燃料の条件に基づく線源強度の不確かさを評価するステップと、
所定の基準となる線源強度を、前記計算コードによる線源強度の不確かさと、前記燃料の条件に基づく線源強度の不確かさとによって補正した補正後の線源強度を算出するステップと、
を有する線源強度の評価方法。
【請求項2】
前記計算コードによる線源強度の不確かさを評価するステップでは、前記核種の生成量の実測値に基づく線源強度のばらつきと前記核データ起因の線源強度の不確かさの関係式に基づいて、前記関係式の係数を算出し、前記関係式によって前記不確かさを校正する、
請求項１に記載の線源強度の評価方法。
請求項１に記載の線源強度の評価方法。
【請求項3】
前記計算コードに基づく線源強度の不確かさを評価するステップでは、生成量を測定することができない第１の核種と前記第１の核種を生成する又は前記第１の核種の壊変等によって生じる生成量を測定することができる第２の核種の関係に基づいて、前記第２の核種の数密度の変動に対する前記第１の核種の変動を表す感度と、前記第２の核種の生成量の実測値に基づく線源強度と前記第２の核種の生成量の計算値に基づく線源強度の差と、前記第１の核種の線源強度の不確かさの関係を示す計算式によって、前記第１の核種の線源強度の不確かさを推定し、
推定した前記第１の核種の線源強度の不確かさと前記第１の核種についての前記核データ起因の線源強度の不確かさの関係式に基づいて、前記関係式の係数を算出し、前記関係式によって前記第１の核種の線源強度の不確かさを校正する、
請求項１に記載の線源強度の評価方法。
推定した前記第１の核種の線源強度の不確かさと前記第１の核種についての前記核データ起因の線源強度の不確かさの関係式に基づいて、前記関係式の係数を算出し、前記関係式によって前記第１の核種の線源強度の不確かさを校正する、
請求項１に記載の線源強度の評価方法。
【請求項4】
前記燃料の条件に基づく線源強度の不確かさを評価するステップでは、線源強度に影響を与えるパラメータを変動させて前記計算コードに入力し、前記計算コードが出力する前記核種の生成量から線源強度を算出して、線源強度の不確かさを評価する、
請求項１から請求項３の何れか１項に記載の線源強度の評価方法。
請求項１から請求項３の何れか１項に記載の線源強度の評価方法。
【請求項5】
前記燃料の条件に基づく線源強度の不確かさを評価するステップでは、線源強度が所定値以上の燃料に関する前記パラメータに基づいて線源強度の不確かさを評価する、
請求項４に記載の線源強度の評価方法。
請求項４に記載の線源強度の評価方法。
【請求項6】
前記補正では、基準となる線源強度の最大値に、前記計算コードに基づく線源強度の不確かさの最大値と、前記燃料の条件に基づく線源強度の不確かさの最大値を加えて補正を行う、
請求項１から請求項３の何れか１項に記載の線源強度の評価方法。
請求項１から請求項３の何れか１項に記載の線源強度の評価方法。
【請求項7】
前記核種は、ガンマ線又は中性子線の発生源である、
請求項１から請求項３の何れか１項に記載の線源強度の評価方法。
請求項１から請求項３の何れか１項に記載の線源強度の評価方法。
【請求項8】
核種の壊変や生成を計算する計算コードに入力する核データの値を変動させたときに前記計算コードが算出する前記核種の生成量に基づいて、前記核データ起因の線源強度の不確かさを評価するステップと、
前記核種の生成量の実測値と前記計算コードの計算値の差に基づく線源強度の不確かさと、前記核データ起因の線源強度の不確かさとに基づいて、前記核データ起因の線源強度の不確かさを校正することにより、前記計算コードによる線源強度の不確かさを評価するステップと、
所定の基準となる線源強度を、前記計算コードによる線源強度の不確かさによって補正した補正後の線源強度を算出するステップと、
を有する線源強度の評価方法。
前記核種の生成量の実測値と前記計算コードの計算値の差に基づく線源強度の不確かさと、前記核データ起因の線源強度の不確かさとに基づいて、前記核データ起因の線源強度の不確かさを校正することにより、前記計算コードによる線源強度の不確かさを評価するステップと、
所定の基準となる線源強度を、前記計算コードによる線源強度の不確かさによって補正した補正後の線源強度を算出するステップと、
を有する線源強度の評価方法。
【請求項9】
炉心の燃料のパラメータが記録された炉心設計データベースに基づいて、燃料の照射履歴に基づく線源強度の不確かさを評価するステップと、
所定の基準となる線源強度を、前記照射履歴に基づく線源強度の不確かさによって補正した補正後の線源強度を算出するステップと、
を有する線源強度の評価方法。
所定の基準となる線源強度を、前記照射履歴に基づく線源強度の不確かさによって補正した補正後の線源強度を算出するステップと、
を有する線源強度の評価方法。
【請求項10】
請求項１から請求項３の何れか１項に記載の評価方法で算出した前記補正後の線源強度を入力として遮蔽解析を行う、
遮蔽解析方法。
遮蔽解析方法。
【請求項11】
核種の壊変や生成を計算する計算コードに入力する核データの値を変動させたときに前記計算コードが算出する前記核種の生成量に基づいて、前記核データ起因の線源強度の不確かさを評価する手段と、
前記核種の生成量の実測値と前記計算コードの計算値の差に基づく線源強度の不確かさと、前記核データ起因の線源強度の不確かさとに基づいて、前記核データ起因の線源強度の不確かさを校正することにより、前記計算コードによる線源強度の不確かさを評価する手段と、
炉心の燃料のパラメータが記録された炉心設計データベースに基づいて、前記燃料の条件に基づく線源強度の不確かさを評価する手段、
所定の基準となる線源強度を、前記計算コードによる線源強度の不確かさと、前記燃料の条件に基づく線源強度の不確かさとによって補正した補正後の線源強度を算出する手段と、
を有する線源強度の評価装置。
前記核種の生成量の実測値と前記計算コードの計算値の差に基づく線源強度の不確かさと、前記核データ起因の線源強度の不確かさとに基づいて、前記核データ起因の線源強度の不確かさを校正することにより、前記計算コードによる線源強度の不確かさを評価する手段と、
炉心の燃料のパラメータが記録された炉心設計データベースに基づいて、前記燃料の条件に基づく線源強度の不確かさを評価する手段、
所定の基準となる線源強度を、前記計算コードによる線源強度の不確かさと、前記燃料の条件に基づく線源強度の不確かさとによって補正した補正後の線源強度を算出する手段と、
を有する線源強度の評価装置。
【請求項12】
コンピュータに、
核種の壊変や生成を計算する計算コードに入力する核データの値を変動させたときに前記計算コードが算出する前記核種の生成量に基づいて、前記核データ起因の線源強度の不確かさを評価するステップと、
前記核種の生成量の実測値と前記計算コードの計算値の差に基づく線源強度の不確かさと、前記核データ起因の線源強度の不確かさとに基づいて、前記核データ起因の線源強度の不確かさを校正することにより、前記計算コードによる線源強度の不確かさを評価するステップと、
炉心の燃料のパラメータが記録された炉心設計データベースに基づいて、前記燃料の条件に基づく線源強度の不確かさを評価するステップと、
所定の基準となる線源強度を、前記計算コードによる線源強度の不確かさと、前記燃料の条件に基づく線源強度の不確かさとによって補正した補正後の線源強度を算出するステップと、
を実行させるプログラム。
核種の壊変や生成を計算する計算コードに入力する核データの値を変動させたときに前記計算コードが算出する前記核種の生成量に基づいて、前記核データ起因の線源強度の不確かさを評価するステップと、
前記核種の生成量の実測値と前記計算コードの計算値の差に基づく線源強度の不確かさと、前記核データ起因の線源強度の不確かさとに基づいて、前記核データ起因の線源強度の不確かさを校正することにより、前記計算コードによる線源強度の不確かさを評価するステップと、
炉心の燃料のパラメータが記録された炉心設計データベースに基づいて、前記燃料の条件に基づく線源強度の不確かさを評価するステップと、
所定の基準となる線源強度を、前記計算コードによる線源強度の不確かさと、前記燃料の条件に基づく線源強度の不確かさとによって補正した補正後の線源強度を算出するステップと、
を実行させるプログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本開示は、線源強度の評価方法、遮蔽解析方法、線源強度の評価装置及びプログラムに関する。
【背景技術】
【0002】
原子力プラントの使用済み燃料を貯蔵する施設やキャスク等の遮蔽設計を行う際には、使用済み燃料の線源強度が必要になる。線源強度は計算により評価するが、遮蔽設計時には、実際に格納する使用済み燃料が存在しないため線源強度の正確な値を評価することはできない。また、個々の使用済み燃料について、線源強度を逐一計算することは現実的ではない。遮蔽設計を精度よく行うためには、使用済み燃料のばらつきなど、様々な要因による不確かさを考慮して評価した線源強度が必要となる。
関連する技術として、特許文献１には、核燃料を構成する放射性物質の崩壊熱に関して、不確かさを含めて評価する方法が開示されている。
関連する技術として、特許文献１には、核燃料を構成する放射性物質の崩壊熱に関して、不確かさを含めて評価する方法が開示されている。
【先行技術文献】
【特許文献】
【0003】
【特許文献1】特開２０２０－７６６６９号公報
【発明の概要】
【発明が解決しようとする課題】
【0004】
使用済み燃料の線源強度に関して、不確かさを合理的に考慮して評価する方法が求められている。
【0005】
本開示は、上記課題を解決することができる線源強度の評価方法、遮蔽解析方法、線源強度の評価装置及びプログラムを提供する。
【課題を解決するための手段】
【0006】
本開示に係る線源強度の評価方法は、核種の壊変や生成を計算する計算コードに入力する核データの値を変動させたときに前記計算コードが算出する前記核種の生成量に基づいて、前記核データ起因の線源強度の不確かさを評価するステップと、前記核種の生成量の実測値と前記計算コードの計算値の差に基づく線源強度の不確かさと、前記核データ起因の線源強度の不確かさとに基づいて、前記核データ起因の線源強度の不確かさを校正することにより、前記計算コードによる線源強度の不確かさを評価するステップと、炉心の燃料のパラメータが記録された炉心設計データベースに基づいて、前記燃料の条件に基づく線源強度の不確かさを評価するステップと、所定の基準となる線源強度を、前記計算コードによる線源強度の不確かさと、前記燃料の条件に基づく線源強度の不確かさと、によって補正した補正後の線源強度を算出するステップと、を有する。
また、本開示に係る線源強度の評価方法は、核種の壊変や生成を計算する計算コードに入力する核データの値を変動させたときに前記計算コードが算出する前記核種の生成量に基づいて、前記核データ起因の線源強度の不確かさを評価するステップと、前記核種の生成量の実測値と前記計算コードの計算値の差に基づく線源強度の不確かさと、前記核データ起因の線源強度の不確かさとに基づいて、前記核データ起因の線源強度の不確かさを校正することにより、前記計算コードによる線源強度の不確かさを評価するステップと、所定の基準となる線源強度を、前記計算コードによる線源強度の不確かさによって補正した補正後の線源強度を算出するステップと、を有する。
また、本開示に係る線源強度の評価方法は、炉心の燃料のパラメータが記録された炉心設計データベースに基づいて、燃料の照射履歴に基づく線源強度の不確かさを評価するステップと、所定の基準となる線源強度を、前記照射履歴に基づく線源強度の不確かさによって補正した補正後の線源強度を算出するステップと、を有する。
また、本開示に係る線源強度の評価方法は、核種の壊変や生成を計算する計算コードに入力する核データの値を変動させたときに前記計算コードが算出する前記核種の生成量に基づいて、前記核データ起因の線源強度の不確かさを評価するステップと、前記核種の生成量の実測値と前記計算コードの計算値の差に基づく線源強度の不確かさと、前記核データ起因の線源強度の不確かさとに基づいて、前記核データ起因の線源強度の不確かさを校正することにより、前記計算コードによる線源強度の不確かさを評価するステップと、所定の基準となる線源強度を、前記計算コードによる線源強度の不確かさによって補正した補正後の線源強度を算出するステップと、を有する。
また、本開示に係る線源強度の評価方法は、炉心の燃料のパラメータが記録された炉心設計データベースに基づいて、燃料の照射履歴に基づく線源強度の不確かさを評価するステップと、所定の基準となる線源強度を、前記照射履歴に基づく線源強度の不確かさによって補正した補正後の線源強度を算出するステップと、を有する。
【0007】
本開示に係る遮蔽解析方法では、上記の評価方法で算出した前記補正後の線源強度を入力として遮蔽解析を行う。
【0008】
本開示に係る線源強度の評価装置は、核種の壊変や生成を計算する計算コードに入力する核データの値を変動させたときに前記計算コードが算出する前記核種の生成量に基づいて、前記核データ起因の線源強度の不確かさを評価する手段と、前記核種の生成量の実測値と前記計算コードの計算値の差に基づく線源強度の不確かさと、前記核データ起因の線源強度の不確かさとに基づいて、前記核データ起因の線源強度の不確かさを校正することにより、前記計算コードによる線源強度の不確かさを評価する手段と、炉心の燃料のパラメータが記録された炉心設計データベースに基づいて、前記燃料の条件に基づく線源強度の不確かさを評価する手段、所定の基準となる線源強度を、前記計算コードによる線源強度の不確かさと、前記燃料の条件に基づく線源強度の不確かさとによって補正した補正後の線源強度を算出する手段と、を有する。
【0009】
本開示に係るプログラムは、コンピュータに、核種の壊変や生成を計算する計算コードに入力する核データの値を変動させたときに前記計算コードが算出する前記核種の生成量に基づいて、前記核データ起因の線源強度の不確かさを評価するステップと、前記核種の生成量の実測値と前記計算コードの計算値の差に基づく線源強度の不確かさと、前記核データ起因の線源強度の不確かさと、に基づいて、前記核データ起因の線源強度の不確かさを校正することにより、前記計算コードによる線源強度の不確かさを評価するステップと、炉心の燃料のパラメータが記録された炉心設計データベースに基づいて、前記燃料の条件に基づく線源強度の不確かさを評価するステップと、所定の基準となる線源強度を、前記計算コードによる線源強度の不確かさと、前記燃料の条件に基づく線源強度の不確かさとによって補正した補正後の線源強度を算出するステップと、を実行させる。
【発明の効果】
【0010】
本開示の線源強度の評価方法、線源強度の評価装置及びプログラムによれば、不確かさを含めた線源強度を評価することができる。本開示の遮蔽設計方法によれば、不確かさを含めて評価した線源強度に基づいて、遮蔽設計を行うことができる。
【図面の簡単な説明】
【0011】
【図1】実施形態に係る評価装置の一例を示すブロック図である。
【図2A】実施形態に係る核データ起因の線源強度の不確かさを評価する処理について説明する第１の図である。
【図2B】実施形態に係る核データ起因の線源強度の不確かさを評価する処理について説明する第２の図である。
【図3】実施形態に係る燃焼チェーンの一例を示す図である。
【図4A】実施形態に係る炉心設計データベースに登録されたパラメータの一例を示す図である。
【図4B】実施形態に係る燃料の選定について説明する図である。
【図5】実施形態に係る不確かさを考慮した線量率について説明する図である。
【図6】実施形態に係る線源強度の評価処理の一例を示すフローチャートである。
【図7】実施形態の評価装置のハードウェア構成の一例を示す図である。
【発明を実施するための形態】
【0012】
＜実施形態＞
以下、本開示の線源強度の評価装置について、図１～図７を参照しながら説明する。
（構成）
図１は、実施形態に係る評価装置の一例を示すブロック図である。
評価装置１０は、使用済み燃料の線源強度や線量率を、不確かさを考慮して算出する。使用済み燃料の線源強度や線量率を評価する場合、計算コードを使用して、使用済み燃料から生じる核種（例えば、中性子やガンマ線等の評価対象となる放射線を放出する核種）の生成量を計算し、その生成量から線源強度等を計算する。計算コードには、様々なパラメータを入力するが、例えば、入力パラメータのうちの核データには、誤差・ばらつきが存在するため、どのような値を入力するかによって、計算コードが算出する核種の生成量や、核種の生成量に基づいて計算される線源強度や線量率にばらつきが生じる。また、計算コードが算出する核種の生成量と実際の使用済み燃料から発生する核種の生成量との間には差異がある。評価装置１０は、［１］核データのばらつきに起因する線源強度の不確かさや、［２］計算コードによって計算される核種生成量（核データのばらつきを含む）に基づく線源強度の不確かさ（計算コードによる線源強度の計算値と実際の線源強度のばらつきの差異）を評価する。
以下、本開示の線源強度の評価装置について、図１～図７を参照しながら説明する。
（構成）
図１は、実施形態に係る評価装置の一例を示すブロック図である。
評価装置１０は、使用済み燃料の線源強度や線量率を、不確かさを考慮して算出する。使用済み燃料の線源強度や線量率を評価する場合、計算コードを使用して、使用済み燃料から生じる核種（例えば、中性子やガンマ線等の評価対象となる放射線を放出する核種）の生成量を計算し、その生成量から線源強度等を計算する。計算コードには、様々なパラメータを入力するが、例えば、入力パラメータのうちの核データには、誤差・ばらつきが存在するため、どのような値を入力するかによって、計算コードが算出する核種の生成量や、核種の生成量に基づいて計算される線源強度や線量率にばらつきが生じる。また、計算コードが算出する核種の生成量と実際の使用済み燃料から発生する核種の生成量との間には差異がある。評価装置１０は、［１］核データのばらつきに起因する線源強度の不確かさや、［２］計算コードによって計算される核種生成量（核データのばらつきを含む）に基づく線源強度の不確かさ（計算コードによる線源強度の計算値と実際の線源強度のばらつきの差異）を評価する。
【0013】
また、線源強度を評価する場合、上述の核データ以外にも原子炉の出力や燃焼度などの入力パラメータを計算コードへ入力して核種の生成量を計算し、線源強度および線量率の評価を行うが、実際の使用済み燃料では、燃料の仕様や照射履歴、プラントの運転条件などの違いにより、線源強度に差が生じる。この照射履歴等に基づく差・ばらつきに関し、評価装置１０は、許認可炉心や取替炉心のパラメータが登録された炉心設計データベースから線源強度への感度が高い入力パラメータを抽出して、抽出した入力パラメータのばらつきによって生じる線源強度の不確かさを評価する。これにより、［３］照射履歴等の燃料の条件の違いによる線源強度の不確かさを評価する。このように、評価装置１０は、［１］核データ起因の線源強度の不確かさ、［２］計算コードによる線源強度の不確かさ、［３］照射履歴等の燃料の条件に関する入力パラメータに起因する線源強度の不確かさ（以下、入力起因の不確かさのように記載する場合がある。）、を評価し、これら［１］～［３］の不確かさを考慮して、使用済み燃料の線源強度および線量率を算出する。
【0014】
図示するように、評価装置１０は、入力受付部１１と、制御部１２と、記憶部１３と、を備える。
入力受付部１１は、キーボード、マウス、タッチパネル、ボタン等の入力装置を用いて入力された情報や指示などを受け付ける。入力受付部１１は、受け付けた情報を記憶部１３に記録したり、制御部１２へ出力したりする。
入力受付部１１は、キーボード、マウス、タッチパネル、ボタン等の入力装置を用いて入力された情報や指示などを受け付ける。入力受付部１１は、受け付けた情報を記憶部１３に記録したり、制御部１２へ出力したりする。
【0015】
制御部１２は、上記の［１］～［３］の不確かさを評価する処理、［１］～［３］の不確かさを考慮した線源強度や線量率の最大値を算出する処理などを実行する。また、制御部１２は、処理結果を表示装置に表示したり、電子データとして出力したりする。制御部１２は、第１評価部１２１と、第２評価部１２２と、第３評価部１２３と、遮蔽解析部１２４と、を備える。これらの機能については後述する。
記憶部１３は、計算コード、炉心設計データベース、処理中のデータなどを記憶する。
記憶部１３は、計算コード、炉心設計データベース、処理中のデータなどを記憶する。
【0016】
［１］核データ起因の不確かさの評価
第１評価部１２１は、［１］核データ起因の線源強度の不確かさを評価する。図２Ａに核データの一例を示す。図２Ａの幅ｅ１は、核データの誤差、ばらつきの幅を表している。核データおよびその誤差については、核データライブラリに登録されており、この値を使用することができる。第１評価部１２１は、幅ｅ１の範囲でランダムに複数の核データを選択する。例えば、第１評価部１２１は、ａ１～ｄ１の核データを選択する。核データｂ１は、基準状態の値である。第１評価部１２１は、選択した核データａ１～ｄ１を計算コードに入力する。計算コードは、例えば、ＳＣＡＬＥ(Standardized Computer Analyses for Licensing Evaluation)システムである。ＳＣＡＬＥシステムに、燃料集合体の諸元データ等（例えば、核データや照射履歴に関するパラメータ等）の入力パラメータを入力して実行すると、ＳＣＡＬＥシステムは、共鳴計算、輸送計算、燃焼計算等を行って、ガンマ線や中性子線を放出する核種の生成量を算出する。ＳＣＡＬＥシステムにて、照射期間を指定することで照射から冷却までの時系列の核種生成量を算出することができる。第１評価部１２１は、核データａ１および他の入力パラメータを計算コード（ＳＣＡＬＥシステム）に入力し、計算コードを実行する。第１評価部１２１は、計算コードが出力する核種の生成量に崩壊定数を乗じて線源強度を算出する。第１評価部１２１は、核データｂ１～ｄ１についても同様の計算を行い、核データａ１～ｄ１に対する線源強度を評価する。図２Ｂに、核データａ１～ｄ１について評価した線源強度の一例を示す。図２Ｂのグラフの縦軸は線源強度、横軸は時間である。線源強度ａ２，ｂ２，ｃ２，ｄ２は、それぞれ核データａ１，ｂ１，ｃ１，ｄ１に対応する線源強度の推移を示している。第１評価部１２１は、例えば、このグラフの時間Ｔ（燃料の冷却が完了した時間）における核データに起因する線源強度の不確かさを評価する。例えば、第１評価部１２１は、基準状態の核データｂ１を用いたときの時間Ｔにおける線源強度ｂ２を基準とする標準偏差を算出し、算出した標準偏差σを、［１］核データ起因の線源強度の不確かさとする。例えば、この例の場合、（（ａ２－ｂ２）２＋（ｃ２－ｂ２）２＋（ｄ２－ｂ２）２÷３）の平方根によって、［１］核データ起因の線源強度の不確かさを算出する。
第１評価部１２１は、［１］核データ起因の線源強度の不確かさを評価する。図２Ａに核データの一例を示す。図２Ａの幅ｅ１は、核データの誤差、ばらつきの幅を表している。核データおよびその誤差については、核データライブラリに登録されており、この値を使用することができる。第１評価部１２１は、幅ｅ１の範囲でランダムに複数の核データを選択する。例えば、第１評価部１２１は、ａ１～ｄ１の核データを選択する。核データｂ１は、基準状態の値である。第１評価部１２１は、選択した核データａ１～ｄ１を計算コードに入力する。計算コードは、例えば、ＳＣＡＬＥ(Standardized Computer Analyses for Licensing Evaluation)システムである。ＳＣＡＬＥシステムに、燃料集合体の諸元データ等（例えば、核データや照射履歴に関するパラメータ等）の入力パラメータを入力して実行すると、ＳＣＡＬＥシステムは、共鳴計算、輸送計算、燃焼計算等を行って、ガンマ線や中性子線を放出する核種の生成量を算出する。ＳＣＡＬＥシステムにて、照射期間を指定することで照射から冷却までの時系列の核種生成量を算出することができる。第１評価部１２１は、核データａ１および他の入力パラメータを計算コード（ＳＣＡＬＥシステム）に入力し、計算コードを実行する。第１評価部１２１は、計算コードが出力する核種の生成量に崩壊定数を乗じて線源強度を算出する。第１評価部１２１は、核データｂ１～ｄ１についても同様の計算を行い、核データａ１～ｄ１に対する線源強度を評価する。図２Ｂに、核データａ１～ｄ１について評価した線源強度の一例を示す。図２Ｂのグラフの縦軸は線源強度、横軸は時間である。線源強度ａ２，ｂ２，ｃ２，ｄ２は、それぞれ核データａ１，ｂ１，ｃ１，ｄ１に対応する線源強度の推移を示している。第１評価部１２１は、例えば、このグラフの時間Ｔ（燃料の冷却が完了した時間）における核データに起因する線源強度の不確かさを評価する。例えば、第１評価部１２１は、基準状態の核データｂ１を用いたときの時間Ｔにおける線源強度ｂ２を基準とする標準偏差を算出し、算出した標準偏差σを、［１］核データ起因の線源強度の不確かさとする。例えば、この例の場合、（（ａ２－ｂ２）２＋（ｃ２－ｂ２）２＋（ｄ２－ｂ２）２÷３）の平方根によって、［１］核データ起因の線源強度の不確かさを算出する。
【0017】
［２］計算コードによる不確かさの評価
第２評価部１２２は、［２］計算コードによる線源強度の不確かさを評価する。例えば、ある核種Ａについて、第１評価部１２１によって評価された［１］核データ起因の線源強度の不確かさがσＡであるとする。一方、計算コードで核種生成量を解析したものと同様の燃料集合体について実施した照射後試験（PIE：Post Irradiation Examination）にて測定された核種Ａの生成量の実測値に基づく線源強度のばらつき（例えば、実測値と計算値の差の標準偏差）がΔγ´Ａの場合、第２評価部１２２は、次式（１）を満たすようなαＡとβＡを推定する。
Δγ´Ａ＝αＡ・σＡ＋βＡ ・・・（１）
例えば、Δγ´ＡとσＡの組合せを複数用意して、αＡとβＡを推定してもよい。ガンマ線の線源強度を評価する場合、第２評価部１２２は、ガンマ線を放出する複数の核種Ａ１～Ａｎのうちの核種Ａ１について、αＡ１とβＡ１を推定し、核種Ａ２についてもαＡ２とβＡ２を推定し、・・・といった処理を核種Ａ１～Ａｎのそれぞれについて実行する。すると、評価対象のガンマ線の線源強度は、次式（２）によって算出することができる。
第２評価部１２２は、［２］計算コードによる線源強度の不確かさを評価する。例えば、ある核種Ａについて、第１評価部１２１によって評価された［１］核データ起因の線源強度の不確かさがσＡであるとする。一方、計算コードで核種生成量を解析したものと同様の燃料集合体について実施した照射後試験（PIE：Post Irradiation Examination）にて測定された核種Ａの生成量の実測値に基づく線源強度のばらつき（例えば、実測値と計算値の差の標準偏差）がΔγ´Ａの場合、第２評価部１２２は、次式（１）を満たすようなαＡとβＡを推定する。
Δγ´Ａ＝αＡ・σＡ＋βＡ ・・・（１）
例えば、Δγ´ＡとσＡの組合せを複数用意して、αＡとβＡを推定してもよい。ガンマ線の線源強度を評価する場合、第２評価部１２２は、ガンマ線を放出する複数の核種Ａ１～Ａｎのうちの核種Ａ１について、αＡ１とβＡ１を推定し、核種Ａ２についてもαＡ２とβＡ２を推定し、・・・といった処理を核種Ａ１～Ａｎのそれぞれについて実行する。すると、評価対象のガンマ線の線源強度は、次式（２）によって算出することができる。
【0018】
【数1】
【0019】
式（２）で算出できる線源強度は、［２］計算コードによる線源強度の不確かさを考慮した線源強度である。式（２）のγｃａｌｃは、従来の方法で計算した線源強度である。Δγは、［２］計算コードによる線源強度の不確かさである。上記のΔγの算出過程から理解できるように、Δγは、［１］核コード起因の線源強度の不確かさを含んでいる。式（２）のシグマの項は、上記の例の場合であれば、核種Ａ１～Ａｎについての核データ起因の線源強度のばらつきを式（１）によって校正した値の和である。核データの誤差によって生じる不確かさ、計算コードと実測値の差によって生じる不確かさ、という根拠に基づいた補正項Δγを加えることで、従来の方法で算出していたγｃａｌｃに比べて安全側に評価した、合理的根拠を有する線源強度を算出することができる。評価対象の線源強度に関する核種Ａ１～Ａｎについて、αＡ１、βＡ１等を算出すれば、それ以降は、式（２）によって、［２］計算コードによる線源強度の不確かさを算出することができる。
【0020】
（照射後試験にて実測値が得られない場合）
ところで、評価対象とする核種生成量の実測値が照射後試験にて測定できない場合がある。このような場合には燃焼チェーンを利用して、未測定の核種による線源強度を評価する。図３に燃焼チェーンの一例を示す。核種Ａは、照射後試験にて測定できない評価対象の核種（例えば、ガンマ線源となる核種）である。核種Ｂは、照射後試験による生成量の測定が可能であって、核分裂によって核種Ａを生成するアクニチド核種である。核種Ｃは、照射後試験による生成量の測定が可能であって、核種Ａの中性子捕獲、壊変などによって生じる核種である。ここで、核種Ａとの間に燃焼チェーンを用いて記述できる生成や壊変等の関係がある核種ｉについて、照射後試験にて測定された核種生成量の実測値に基づく線源強度と、計算コードで計算した核種生成量の計算値に基づく線源強度の差をΔγｉとする。生成量の実測値が得られる核種ｉの数密度が１％変動したときの核種Ａの数密度の変動（％）をＷｉ→Ａとする。Ｗｉ→Ａは計算コードで算出することができる。核種Ａの線源強度のばらつき、不確かさをΔγ´Ａとする。すると、核種Ｂから核種Ａが生成され、核種Ａから核種Ｃが生成されるという図３に例示する燃焼チェーンの関係から次式（３）が導出される。
ところで、評価対象とする核種生成量の実測値が照射後試験にて測定できない場合がある。このような場合には燃焼チェーンを利用して、未測定の核種による線源強度を評価する。図３に燃焼チェーンの一例を示す。核種Ａは、照射後試験にて測定できない評価対象の核種（例えば、ガンマ線源となる核種）である。核種Ｂは、照射後試験による生成量の測定が可能であって、核分裂によって核種Ａを生成するアクニチド核種である。核種Ｃは、照射後試験による生成量の測定が可能であって、核種Ａの中性子捕獲、壊変などによって生じる核種である。ここで、核種Ａとの間に燃焼チェーンを用いて記述できる生成や壊変等の関係がある核種ｉについて、照射後試験にて測定された核種生成量の実測値に基づく線源強度と、計算コードで計算した核種生成量の計算値に基づく線源強度の差をΔγｉとする。生成量の実測値が得られる核種ｉの数密度が１％変動したときの核種Ａの数密度の変動（％）をＷｉ→Ａとする。Ｗｉ→Ａは計算コードで算出することができる。核種Ａの線源強度のばらつき、不確かさをΔγ´Ａとする。すると、核種Ｂから核種Ａが生成され、核種Ａから核種Ｃが生成されるという図３に例示する燃焼チェーンの関係から次式（３）が導出される。
【0021】
【数2】
【0022】
図３の燃焼チェーンの場合、第２評価部１２２は、上式（３）にてｉ＝Ｂ，Ｃとして、生成量が測定できない核種Ａの燃焼強度の不確かさΔγ´Ａを推定する。
【0023】
なお、式（３）にて、ｉ＝Ａ，Ｂ，Ｃとしているが、これは、核種Ａの生成量が測定できる場合であっても、その値の信頼度が低いような場合であれば、上式（３）にてｉ＝Ａ，Ｂ，ＣとしてΔγ´Ａを計算することで、核種Ａの燃焼強度の不確かさをより精度よく推定することができるためである。
【0024】
第２評価部１２２は、式（３）を用いてΔγ´Ａを推定すると、照射後試験で核種生成量が測定できる場合と同様にして、上記の式（１）を用いて係数αＡ、βＡを算出する。そして、上記の式（２）によって、［２］計算コードによる線源強度の不確かさを考慮した線源強度を算出する。このように、評価対象の核種生成量が実測できない場合であっても、燃焼チェーンに基づく関係式（３）を利用することで、［２］計算コードによる線源強度の不確かさを考慮した線源強度の算出が可能である。
【0025】
［３］入力起因の不確かさの評価
第３評価部１２３は、［３］入力起因の線源強度の不確かさを評価する。これまでに、核データのばらつきや線源強度の計算値のばらつきと実測値のばらつきの差異に基づく不確かさを評価する方法について説明した。これに対し、第３評価部１２３は、燃料の仕様や照射履歴、運転環境などの差によって生じる線源強度の不確かさを評価する。炉心設計データベースには、許認可炉心や取替炉心など様々な炉心について、原子炉の出力、ほう素濃度、冷却期間、ＢＰ（バーナブルポイズン）の有無、燃料の組成などの情報が登録されている。第３評価部１２３は、炉心設計データベースに登録された、燃料の仕様や照射履歴、運転環境などに関する情報の中から線源強度の影響を与えるパラメータを抽出して解析し、線源強度の影響を与えるパラメータのばらつきを評価する。第３評価部１２３は、炉心設計データベースに登録されたパラメータを、線源強度への影響度に応じて「Ｈ」、「Ｍ」、「Ｌ」に分類する。「Ｈ」は線源強度への影響度が高いパラメータ、「Ｍ」は線源強度への影響度が中程度のパラメータ、「Ｌ」は線源強度への影響度が低いパラメータである。図４Ａにガンマ線に関して「Ｈ」、「Ｍ」、「Ｌ」に分類されたパラメータの一例を示す。「Ｈ」に分類されるパラメータは、例えば“出力”、「Ｍ」に分類されるパラメータは、例えば“燃焼度”、“組成”、“ＢＰの有無”、“照射期間”、「Ｌ」に分類されるパラメータは、例えば“ほう素濃度”、“隣接燃料”、“燃料温度”、“減速材温度”などとなっている。「Ｈ」、「Ｍ」、「Ｌ」の分類は、知見ある技術者が設定を行い、この設定に従って、第３評価部１２３が、線源強度への影響度に応じたパラメータを炉心設計データベースから抽出する。
第３評価部１２３は、［３］入力起因の線源強度の不確かさを評価する。これまでに、核データのばらつきや線源強度の計算値のばらつきと実測値のばらつきの差異に基づく不確かさを評価する方法について説明した。これに対し、第３評価部１２３は、燃料の仕様や照射履歴、運転環境などの差によって生じる線源強度の不確かさを評価する。炉心設計データベースには、許認可炉心や取替炉心など様々な炉心について、原子炉の出力、ほう素濃度、冷却期間、ＢＰ（バーナブルポイズン）の有無、燃料の組成などの情報が登録されている。第３評価部１２３は、炉心設計データベースに登録された、燃料の仕様や照射履歴、運転環境などに関する情報の中から線源強度の影響を与えるパラメータを抽出して解析し、線源強度の影響を与えるパラメータのばらつきを評価する。第３評価部１２３は、炉心設計データベースに登録されたパラメータを、線源強度への影響度に応じて「Ｈ」、「Ｍ」、「Ｌ」に分類する。「Ｈ」は線源強度への影響度が高いパラメータ、「Ｍ」は線源強度への影響度が中程度のパラメータ、「Ｌ」は線源強度への影響度が低いパラメータである。図４Ａにガンマ線に関して「Ｈ」、「Ｍ」、「Ｌ」に分類されたパラメータの一例を示す。「Ｈ」に分類されるパラメータは、例えば“出力”、「Ｍ」に分類されるパラメータは、例えば“燃焼度”、“組成”、“ＢＰの有無”、“照射期間”、「Ｌ」に分類されるパラメータは、例えば“ほう素濃度”、“隣接燃料”、“燃料温度”、“減速材温度”などとなっている。「Ｈ」、「Ｍ」、「Ｌ」の分類は、知見ある技術者が設定を行い、この設定に従って、第３評価部１２３が、線源強度への影響度に応じたパラメータを炉心設計データベースから抽出する。
【0026】
線源強度への影響度に応じたパラメータを炉心設計データベースから抽出すると、第３評価部１２３は、各パラメータについて値を設定する。例えば、第３評価部１２３は、「Ｈ」や「Ｍ」に分類されるパラメータについては、炉心設計データベースからランダムに抽出した値を設定する。第３評価部１２３は、例えば、各パラメータの最大値や最小値を含むようにして炉心設計データベースから値を抽出してもよい。また、第３評価部１２３は、「Ｌ」に分類されるパラメータについては、代表値を設定する。例えば、“燃料温度”の代表値として９００Ｋを設定し、“減速材温度”の代表値には、炉心特性説明書に記載された値を設定する。これらの値は、知見を有する技術者が設定を行い、この設定に従って、第３評価部１２３が、代表値を各パラメータへ設定する。
【0027】
第３評価部１２３は、「Ｈ」、「Ｍ」のパラメータにはランダムに抽出した値を設定し、「Ｌ」のパラメータには代表値を設定し、これらのパラメータを計算コードに入力して、例えば、ガンマ線を放出する核種の生成量を計算し、各種計算量に崩壊定数を乗じることにより、核種ごとに線源強度を算出する。第３評価部１２３は、「Ｈ」、「Ｍ」のパラメータに設定する値を変化させて（「Ｌ」は代表値、つまり固定値）、繰り返し、同様にして線源強度を計算する。これにより、複数の線源強度が算出される。第３評価部１２３は、算出した複数の線源強度から、燃料の仕様や照射履歴のばらつきに基づく線源強度のばらつき（つまり、［３］入力起因の不確かさ）を評価する。安全性を担保するために線源強度が高くなる場合のばらつき具合を評価したいので、第３評価部１２３は、例えば、次に説明するように線源強度が高い計算結果を抽出して、抽出した計算結果に係る燃料を対象としてばらつきを評価する。なお、線源強度は、遮蔽解析により線量率に変換されるため、以下では、線源強度を線量率に変換して、ばらつきを評価する。
【0028】
例えば、第３評価部１２３は、基準となる許認可炉心の最大線量率に対する不確かさを評価するため、線量率が高い燃料を評価対象として選定する。図４Ｂに様々な値を「Ｈ」、「Ｍ」のパラメータに設定して計算した線源強度を線量率に変換したデータを示す。図４Ｂのグラフの横軸は計算結果の線量率と基準となる許認可炉心の最大線量率との差、縦軸は各差に対応する計算結果の頻度を示している。横軸の値が「０」のデータが、許認可炉心の最大線量率と同じ線量率となったときの計算結果の頻度である。第３評価部１２３は、この「０」を基準として線量率が高い燃料を選定する。範囲Ｐ１は、選定された燃料である。第３評価部１２３は、選定した燃料の線量率を線源強度に変換し、変換後の線源強度を標本として、ブートストラップ法により、線源強度の不確かさの信頼係数を評価する。例えば、第３評価部は、信頼係数９５％（信頼度９５％）の信頼区間を算出し、この区間幅を（３）入力起因による線源強度の不確かさとする。また、第３評価部は、線源強度の変化に対する線量率の感度、線源強度のばらつき及び信頼係数から線量率の不確かさを評価する。
【0029】
遮蔽解析部１２４は、第２評価部１２２によって評価された［２］計算コードによる線源強度の不確かさと、第３評価部１２３によって評価された［３］入力起因の線源強度の不確かさとによって、所定の基準となる線源強度（例えば、許認可炉心などの代表的な炉心の線源強度）を補正して補正後の線源強度を算出する。例えば、遮蔽解析部１２４は、基準となる線源強度に、［２］計算コードによる線源強度の不確かさの最大値と［３］入力起因の線源強度の不確かさの最大値を加算して補正後の線源強度を算出する。そして、遮蔽解析部１２４は、補正後の線源強度を入力として、使用済み燃料の貯蔵施設やキャスクの遮蔽解析を行って遮蔽部材の厚さなどを設計し、線量率の評価を行う。
【0030】
図５に線量率の評価結果を示す。図５の横軸は線量率を示している。ｒ０は、不確かさを考慮しない許認可炉心の最大線量率である。ｒ１は、[２]計算コードによる不確かさ（[１]核コード起因の不確かさを含む）を加算した線源強度を入力として遮蔽解析を行って得られる線量率である。ｒ１は、合理的な根拠に基づく許認可炉心の最大線量率である。ｒ２は、許認可炉心の最大線量率に相当する線源強度に[２]計算コードによる不確かさと[３]入力起因の不確かさを加算した線源強度を入力として遮蔽解析を行って得られた線量率である。ｒ２は、不確かさを考慮した線量率評価値である。実際の使用済み燃料は、許認可炉心と仕様や照射履歴などが異なるが、その差は、[３]入力起因の不確かさによって吸収される。ｒ０を用いてキャスク等の遮蔽解析を行うと、安全性に影響が出る可能性が高いため、従来は、過度に保守的に線源強度を見積もって遮蔽解析を行っていた。しかし、本実施形態によれば、核データ起因や計算コードによる線源強度の不確かさ、照射履歴等のばらつきによる線源強度の不確かさ、を考慮することによって、線源強度の不確かさを合理的に評価する。これにより、過度に保守的にならず、且つ適切に安全側に評価した線源強度を算出することができる。実際、従来の線源強度を用いた場合と線量率の比較を行うと、本実施形態の方法で評価した線源強度を入力として評価した線量率は、従来の保守的な線源強度を入力として評価した線量率よりも下回ることが確認できた。
【0031】
（動作）
次に図６を用いて、本実施形態の線源強度の評価処理流れを説明する。
図６は、実施形態に係る線源強度の評価処理の一例を示すフローチャートである。
前提として、記憶部１３には、計算コード、核データのライブラリ、炉心設計データベース、計算コードに入力する入力パラメータに対する「Ｈ」、「Ｍ」、「Ｌ」の設定、「Ｌ」のパラメータに設定する代表値が登録されている。
まず、オペレータが評価対象の放射線（ガンマ線、中性子線など）を放出する核種を指定し、線源強度の評価を指示する。入力受付部１１は、この指示を受け付ける。制御部１２は、線源強度を評価する処理を実行する。まず、第１評価部１２１が、［１］核データ起因の線源強度の不確かさを評価する（ステップＳ１）。例えば、第１評価部１２１は、評価対象の核種の反応断面積について、ランダムサンプリングにより、複数の値を核データライブラリの誤差データの中から選択し、選択した値を入力パラメータとして、選択した値の数だけ計算コードを実行し、核種生成量を算出する。第１評価部１２１は、選択した値の数だけ得られた核種生成量を線源強度に変換し、そのデータライブラリの誤差データに起因する標準偏差（核データ起因の線源強度の不確かさ）を算出する。第１評価部１２１は、核種ごとに核データ起因の不確かさを算出して、記憶部１３に記録する。
次に図６を用いて、本実施形態の線源強度の評価処理流れを説明する。
図６は、実施形態に係る線源強度の評価処理の一例を示すフローチャートである。
前提として、記憶部１３には、計算コード、核データのライブラリ、炉心設計データベース、計算コードに入力する入力パラメータに対する「Ｈ」、「Ｍ」、「Ｌ」の設定、「Ｌ」のパラメータに設定する代表値が登録されている。
まず、オペレータが評価対象の放射線（ガンマ線、中性子線など）を放出する核種を指定し、線源強度の評価を指示する。入力受付部１１は、この指示を受け付ける。制御部１２は、線源強度を評価する処理を実行する。まず、第１評価部１２１が、［１］核データ起因の線源強度の不確かさを評価する（ステップＳ１）。例えば、第１評価部１２１は、評価対象の核種の反応断面積について、ランダムサンプリングにより、複数の値を核データライブラリの誤差データの中から選択し、選択した値を入力パラメータとして、選択した値の数だけ計算コードを実行し、核種生成量を算出する。第１評価部１２１は、選択した値の数だけ得られた核種生成量を線源強度に変換し、そのデータライブラリの誤差データに起因する標準偏差（核データ起因の線源強度の不確かさ）を算出する。第１評価部１２１は、核種ごとに核データ起因の不確かさを算出して、記憶部１３に記録する。
【0032】
次に第２評価部１２２が、［２］計算コードによる線源強度の不確かさを評価する（ステップＳ２）。照射後試験で評価対象の核種の生成量が測定できる場合、第２評価部１２２は、核種の生成量の実測値に基づく線源強度のばらつき（Δγ´Ａ）と、核データ起因の線源強度の不確かさ（σＡ）の関係式（上記の式（１））に基づいて、関係式の係数αＡ、βＡを算出する。第２評価部１２２は、核種ごとに式（１）の係数αＡ、βＡを算出
し、式（２）によって、［２］計算コードによる線源強度の不確かさを算出する。
し、式（２）によって、［２］計算コードによる線源強度の不確かさを算出する。
【0033】
照射後試験で評価対象の核種の生成量が測定できない場合、第２評価部１２２は、生成量が測定できない核種Ａと、核種Ａを生成する測定可能な核種Ｂ又は核種Ａの壊変等によって生じる測定可能な核種Ｃの関係に基づいて、核種ｉ（ｉ＝Ｂ，Ｃ）の数密度の変動に対する核種Ａの変動を表す感度（Ｗｉ→Ａ）と、核種ｉの生成量の実測値に基づく線源強度と計算値に基づく線源強度の差（Δγｉ）と、核種Ａの線源強度の不確かさ（Δγ´Ａ）の関係を示す計算式（上記の式（３））によって、核種Ａの線源強度の不確かさを推定する。そして、第２評価部１２２は、上記の式（１）によって係数αＡ、βＡを算出し、式（２）によって、［２］計算コードによる線源強度の不確かさを算出する。第２評価部１２２は、算出した計算コードによる線源強度の不確かさを記憶部１３に記録する。
【0034】
次に第３評価部１２３が、［３］入力起因の線源強度の不確かさを評価する（ステップＳ３）。まず、第３評価部１２３は、炉心設計データベースに登録されたパラメータを「Ｈ」、「Ｍ」、「Ｌ」に分類する。第３評価部１２３は、「Ｈ」と「Ｍ」に分類されたパラメータについては、炉心設計データベースからランダムに抽出した値を設定し、「Ｌ」に分類されたパラメータについては、所定の代表値を設定する。第３評価部１２３は、各パラメータに設定した値を組み合わせて、計算コードに入力し核種生成量を算出し、核種生成量から線源強度を算出する。第３評価部１２３は、各パラメータに設定した値の組み合わせを変えて複数回、核種生成量および線源強度を算出する。これにより、燃料の仕様や照射履歴のばらつきに応じた線源強度が算出される。次に第３評価部１２３は、算出した複数の線源強度を線量率に変換し、許認可炉心など代表的な炉心の線量率の最大値に基づいて、線量率が所定値以上となるときの各パラメータの値の組み合わせが示す燃料を選定する。図４Ｂの例では、対象となる燃料の線量率から許認可炉心燃焼の線量率を引いた値が－０．０５μＳｖ／ｈ以上となる燃料を選定する。（この値は一例であって、これに限定されない。）。第３評価部１２３は、選定した燃料の線量率を線源強度に代えて、選定した燃料の線源強度を標本としてブートストラップ法などにより、［３］入力起因による線源強度の不確かさ（９５％信頼区間など）を算出する。第３評価部１２３は、算出した入力起因による線源強度の不確かさを記憶部１３に記録する。
【0035】
次に遮蔽解析部１２４が、不確かさを考慮した線量率を評価する（ステップＳ４）。遮蔽解析部１２４は、許認可炉心の線源強度の最大値に、記憶部１３に記録された計算コードによる線源強度の不確かさの最大値と、入力起因による線源強度の不確かさの最大値とを加算した値を入力として使用済み燃料を格納する施設やキャスクの遮蔽解析を行い、許認可炉心の最大線量率（図５のｒ０）に基づく、不確かさを考慮した線量率（図５のｒ２）を算出する。
【0036】
（効果）
以上、説明したように、本実施形態によれば、計算コード及び入力起因の不確かさを合理的に考慮したガンマ線や中性子線の線源強度を算出することができる。これにより、合理的且つ安全側に評価された線源強度に基づいて、ＳＦＰ（Spent fuel pools）やキャスクの遮蔽設計を行うことができる。
以上、説明したように、本実施形態によれば、計算コード及び入力起因の不確かさを合理的に考慮したガンマ線や中性子線の線源強度を算出することができる。これにより、合理的且つ安全側に評価された線源強度に基づいて、ＳＦＰ（Spent fuel pools）やキャスクの遮蔽設計を行うことができる。
【0037】
なお、上記の実施形態では、許認可炉心の線源強度に［１］～［３］の不確かさを積み上げて線源強度を評価することとしたが、線源強度を評価する目的に応じて、例えば、［１］と［２］の不確かさだけを積み上げて線源強度を評価してもよいし、［３］の不確かさだけを積み上げて線源強度を評価してもよい。また、［１］と［２］の不確かさと、［３］以外の方法で評価した不確かさを積み上げて線源強度を評価してもよいし、［３］の不確かさと、［１］、［２］以外の方法で評価した不確かさを積み上げて線源強度を評価してもよい。
【0038】
図７は、評価装置のハードウェア構成の一例を示す図である。
コンピュータ９００は、ＣＰＵ９０１、主記憶装置９０２、補助記憶装置９０３、入出力インタフェース９０４、通信インタフェース９０５を備える。
上述の評価装置１０は、コンピュータ９００に実装される。そして、上述した各機能は、プログラムの形式で補助記憶装置９０３に記憶されている。ＣＰＵ９０１は、プログラムを補助記憶装置９０３から読み出して主記憶装置９０２に展開し、当該プログラムに従って上記処理を実行する。また、ＣＰＵ９０１は、プログラムに従って、記憶領域を主記憶装置９０２に確保する。また、ＣＰＵ９０１は、プログラムに従って、処理中のデータを記憶する記憶領域を補助記憶装置９０３に確保する。
コンピュータ９００は、ＣＰＵ９０１、主記憶装置９０２、補助記憶装置９０３、入出力インタフェース９０４、通信インタフェース９０５を備える。
上述の評価装置１０は、コンピュータ９００に実装される。そして、上述した各機能は、プログラムの形式で補助記憶装置９０３に記憶されている。ＣＰＵ９０１は、プログラムを補助記憶装置９０３から読み出して主記憶装置９０２に展開し、当該プログラムに従って上記処理を実行する。また、ＣＰＵ９０１は、プログラムに従って、記憶領域を主記憶装置９０２に確保する。また、ＣＰＵ９０１は、プログラムに従って、処理中のデータを記憶する記憶領域を補助記憶装置９０３に確保する。
【0039】
なお、評価装置１０の全部または一部の機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより各機能部による処理を行ってもよい。ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。また、「コンピュータシステム」は、ＷＷＷシステムを利用している場合であれば、ホームページ提供環境（あるいは表示環境）も含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、ＣＤ、ＤＶＤ、ＵＳＢ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。また、このプログラムが通信回線によってコンピュータ９００に配信される場合、配信を受けたコンピュータ９００が当該プログラムを主記憶装置９０２に展開し、上記処理を実行しても良い。また、上記プログラムは、前述した機能の一部を実現するためのものであっても良く、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであってもよい。
【0040】
以上のとおり、本開示に係るいくつかの実施形態を説明したが、これら全ての実施形態は、例として提示したものであり、発明の範囲を限定することを意図していない。これらの実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で種々の省略、置き換え、変更を行うことができる。これらの実施形態及びその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれる。
【0041】
＜付記＞
実施形態に記載の線源強度の評価方法、遮蔽解析方法、線源強度の評価装置及びプログラムは、例えば以下のように把握される。
実施形態に記載の線源強度の評価方法、遮蔽解析方法、線源強度の評価装置及びプログラムは、例えば以下のように把握される。
【0042】
（１）第１の態様に係る線源強度の評価方法は、核種の壊変や生成を計算する計算コードに入力する核データの値を変動させたときに前記計算コードが算出する前記核種の生成量に基づいて、前記核データ起因の線源強度の不確かさを評価するステップと、前記核種の生成量の実測値と前記計算コードの計算値の差に基づく線源強度の不確かさと、前記核データ起因の線源強度の不確かさとに基づいて、前記核データ起因の線源強度の不確かさを校正することにより、前記計算コードによる線源強度の不確かさを評価するステップと、炉心の燃料のパラメータが記録された炉心設計データベースに基づいて前記燃料の条件に基づく線源強度の不確かさを評価するステップと、所定の基準となる線源強度を、前記計算コードによる線源強度の不確かさと、前記燃料の条件に基づく線源強度の不確かさとによって補正した補正後の線源強度を算出するステップと、を有する。
計算コード及び入力起因の不確かさを考慮した線源強度を算出することにより、遮蔽設計に用いる不確かさを合理的に考慮したガンマ線・中性子線等の線源強度を算出することができる。
計算コード及び入力起因の不確かさを考慮した線源強度を算出することにより、遮蔽設計に用いる不確かさを合理的に考慮したガンマ線・中性子線等の線源強度を算出することができる。
【0043】
（２）第２の態様に係る線源強度の評価方法は、（１）の評価方法であって、前記計算コードによる線源強度の不確かさを評価するステップでは、前記核種の生成量の実測値に基づく線源強度のばらつきと前記核データ起因の線源強度の不確かさの関係式（式（１））に基づいて、前記関係式の係数を算出し、前記関係式によって前記不確かさを校正する。
これにより、計算コードによる線源強度の不確かさを算出することができる。
これにより、計算コードによる線源強度の不確かさを算出することができる。
【0044】
（３）第３の態様に係る線源強度の評価方法は、（１）の評価方法であって、前記計算コードに基づく線源強度の不確かさを評価するステップでは、生成量を測定できない第１の核種と前記第１の核種を生成する又は前記第１の核種の壊変等によって生じる生成量を測定することができる第２の核種の関係に基づいて、前記第２の核種の数密度の変動に対する前記第１の核種の変動を表す感度と、前記第２の核種の生成量の実測値に基づく線源強度と前記第２の核種の生成量の計算値に基づく線源強度の差と、前記第１の核種の線源強度の不確かさの関係を示す計算式（式（３））によって、前記第１の核種の線源強度の不確かさを推定し、推定した前記第１の核種の線源強度の不確かさと前記第１の核種についての前記核データ起因の線源強度の不確かさの関係式（式（１））に基づいて、前記関係式の係数を算出し、前記関係式によって前記第１の核種の線源強度の不確かさを校正する。
これにより、線源強度の元となる核種が測定できない場合でも、当該核種に関する計算コードによる線源強度の不確かさを算出することができる。
これにより、線源強度の元となる核種が測定できない場合でも、当該核種に関する計算コードによる線源強度の不確かさを算出することができる。
【0045】
（４）第４の態様に係る線源強度の評価方法は、（１）～（３）の評価方法であって、前記燃料の条件に基づく線源強度の不確かさを評価するステップでは、前記線源強度に影響を与えるパラメータを変動させて前記計算コードに入力し、前記計算コードが出力する前記核種の生成量から線源強度を算出して、前記線源強度の不確かさを評価する。
これにより、様々な燃料の仕様や照射履歴に基づく線源強度の不確かさを算出することができる。
これにより、様々な燃料の仕様や照射履歴に基づく線源強度の不確かさを算出することができる。
【0046】
（５）第５の態様に係る線源強度の評価方法は、（１）～（４）の評価方法であって、前記燃料の条件に基づく線源強度の不確かさを評価するステップでは、線源強度が所定値以上の燃料に関する前記パラメータに基づいて線源強度の不確かさを評価する。
これにより、線源強度を安全側に評価することができる。
これにより、線源強度を安全側に評価することができる。
【0047】
（６）第６の態様に係る線源強度の評価方法は、（１）～（５）の評価方法であって、前記補正では、基準となる線源強度の最大値に、前記計算コードに基づく線源強度の不確かさの最大値と、前記燃料の条件に基づく線源強度の不確かさの最大値を加えて補正を行う。
これにより、不確かさを合理的且つ安全側に考慮した線源強度を算出することができる。
これにより、不確かさを合理的且つ安全側に考慮した線源強度を算出することができる。
【0048】
（７）第７の態様に係る線源強度の評価方法は、（１）～（６）の評価方法であって、前記核種は、ガンマ線又は中性子線の発生源である。
これにより、ＳＦＰの遮蔽解析で評価対象となるガンマ線、キャスクの遮蔽解析で評価対象となるガンマ線および中性子線の線源強度を評価することができる。
これにより、ＳＦＰの遮蔽解析で評価対象となるガンマ線、キャスクの遮蔽解析で評価対象となるガンマ線および中性子線の線源強度を評価することができる。
【0049】
（８）第８の態様に係る線源強度の評価方法は、核種の壊変や生成を計算する計算コードに入力する核データの値を変動させたときに前記計算コードが算出する前記核種の生成量に基づいて、前記核データ起因の線源強度の不確かさを評価するステップと、前記核種の生成量の実測値と前記計算コードの計算値の差に基づく線源強度の不確かさと、前記核データ起因の線源強度の不確かさとに基づいて、前記核データ起因の線源強度の不確かさを校正することにより、前記計算コードによる線源強度の不確かさを評価するステップと、所定の基準となる線源強度を、前記計算コードによる線源強度の不確かさによって補正した補正後の線源強度を算出するステップと、を有する。
これにより、計算コードの不確かさを考慮した線源強度を提供することができる。
これにより、計算コードの不確かさを考慮した線源強度を提供することができる。
【0050】
（９）第９の態様に係る線源強度の評価方法は、炉心の燃料のパラメータが記録された炉心設計データベースに基づいて、燃料の照射履歴に基づく線源強度の不確かさを評価するステップと、所定の基準となる線源強度を、前記照射履歴に基づく線源強度の不確かさによって補正した補正後の線源強度を算出するステップと、を有する。
これにより、入力起因の不確かさを考慮した線源強度を提供することができる。
これにより、入力起因の不確かさを考慮した線源強度を提供することができる。
【0051】
（１０）第１０の態様に係る遮蔽解析方法は、（１）～（９）の何れかに記載の線源強度の評価方法で算出した前記補正後の線源強度を入力として遮蔽解析を行う。
計算コード及び入力起因の不確かさを考慮した線源強度を入力とした遮蔽解析を行うことにより、安全側に評価した線源強度に基づいて、使用済み燃料の貯蔵施設や貯蔵容器の遮蔽設計を精度よく行うことができる。
計算コード及び入力起因の不確かさを考慮した線源強度を入力とした遮蔽解析を行うことにより、安全側に評価した線源強度に基づいて、使用済み燃料の貯蔵施設や貯蔵容器の遮蔽設計を精度よく行うことができる。
【0052】
（１１）第１１の態様に係る評価装置は、核種の壊変や生成を計算する計算コードに入力する核データの値を変動させたときに前記計算コードが算出する前記核種の生成量に基づいて、前記核データ起因の線源強度の不確かさを評価する手段、前記核種の生成量の実測値と前記計算コードの計算値の差に基づく線源強度の不確かさと、前記核データ起因の線源強度の不確かさとに基づいて、前記核データ起因の線源強度の不確かさを校正することにより、前記計算コードによる線源強度の不確かさを評価する手段、炉心の燃料のパラメータが記録された炉心設計データベースに基づいて、前記燃料の条件に基づく線源強度の不確かさを評価する手段、所定の基準となる線源強度を、前記計算コードによる線源強度の不確かさと、前記燃料の条件に基づく線源強度の不確かさとによって補正した補正後の線源強度を算出する手段、を有する。
【0053】
（１２）第１２の態様に係るプログラムは、コンピュータに、核種の壊変や生成を計算する計算コードに入力する核データの値を変動させたときに前記計算コードが算出する前記核種の生成量に基づいて、前記核データ起因の線源強度の不確かさを評価するステップと、前記核種の生成量の実測値と前記計算コードの計算値の差に基づく線源強度の不確かさと、前記核データ起因の線源強度の不確かさと、に基づいて、前記核データ起因の線源強度の不確かさを校正することにより、前記計算コードによる線源強度の不確かさを評価するステップと、炉心の燃料のパラメータが記録された炉心設計データベースに基づいて、前記燃料の条件に基づく線源強度の不確かさを評価するステップと、所定の基準となる線源強度を、前記計算コードによる線源強度の不確かさと、前記燃料の条件に基づく線源強度の不確かさとによって補正した補正後の線源強度を算出するステップと、を実行させる。
【符号の説明】
【0054】
１０・・・評価装置
１１・・・入力受付部
１２・・・制御部
１２１・・・第１評価部
１２２・・・第２評価部
１２３・・・第３評価部
１２４・・・遮蔽解析部
１３・・・記憶部
９００・・・コンピュータ
９０１・・・ＣＰＵ
９０２・・・主記憶装置
９０３・・・補助記憶装置
９０４・・・入出力インタフェース
９０５・・・通信インタフェース
１１・・・入力受付部
１２・・・制御部
１２１・・・第１評価部
１２２・・・第２評価部
１２３・・・第３評価部
１２４・・・遮蔽解析部
１３・・・記憶部
９００・・・コンピュータ
９０１・・・ＣＰＵ
９０２・・・主記憶装置
９０３・・・補助記憶装置
９０４・・・入出力インタフェース
９０５・・・通信インタフェース
【図1】
【図2A】
【図2B】
【図3】
【図4A】
【図4B】
【図5】
【図6】
【図7】