NUMOの確率論的評価手法の開発
原子力学会バックエンド部会
第30回「バックエンド」夏期セミナー
2014年 8月 7日（木）ビッグパレットふくしま
原子力発電環境整備機構
技術部 後藤 淳一
P.1

確率論的アプローチの検討の背景¾プレート運動の安定性を前提に，過去～現在の自然現象の変動傾向を将来に外挿し，地層の著しい変動を回避（決定論的アプローチ）⇒ 回避してもなお残る不確実性が存在¾プレート運動の安定性の前提が成り立たない超長期の将来予測には大きな不確実性が伴う⇒ 過度に保守的な評価となる可能性がある
9確率論的アプローチは，不確実性に関する定量的・系統的な情報を提供し，より合理的な評価を可能にする不確実性への対応策として，決定論的なアプローチを補強する位置づけで，確率論的アプローチの検討を開始

P.2

NUMOの確率論的手法開発プロジェクト¾ITM（International Tectonics Meeting）プロジェクト（2004～2008年度） 手法の構築 ⇒ ITM手法
9 プレート運動の安定性を前提とした議論
9 海外で採用されている確率論的モデル・解析手法の適用
9 広域的な火山活動，断層活動の確率論的な空間分布マップの作成¾TOPAZ（Tectonics Of Preliminary Assessment Zones）プロジェクト（2009～2013年度） ITM手法の拡張 ⇒ ITM-TOPAZ手法
9 プレート運動の安定性を前提にできず不確実性が大きくなる超長期9 火山活動，断層活動，隆起・侵食のシナリオの設定
9 ロジックツリーと専門家意見集約によるシナリオの発生確率の設定
9 サイト評価や安全評価に向けた情報提供

P.3

ITM-TOPAZ手法の概要
ステップ1： 長期変遷シナリオの設定

a: 広域変遷シナリオ（RES）
b: サイト変遷シナリオ（SES）
c: 影響シナリオ（IS）

ステップ2： ロジックツリーの構築
ステップ3： 専門家の意見集約
ステップ4： 確率論的評価
SES

専門家の意見集
約によるシナリオの確信度の設定（重み付け）

10 ky - 100 ky
0- 10 ky

処分施設周辺のシナリオの事象とその影響

サイト周辺のローカルな 0.4
IS-1.1.1
事象の変遷

RES

プレート運動に伴う広域的な 0.7
SES-1.1
事象の変遷

IS-1.1.2

0.6

IS-1.1.3

100 ky – 1 mySite

IS0.4

RES-1
RES-2

0.2
0.1

SES-1.2

起こりやすさP = 0.168
x
確率空間
分布マップ
（ITM手法）SES-1.3

ハザードマップ

処分システムへの影響とその発生確率

サイト評価

安全評価
P.4

ステップ１a 広域変遷シナリオ• 過去のプレート運動等の知見に基づき将来のプレート運動と広域的な現象の変遷を設定
10万年～
100万年

RES-1:
プレート運動の傾向が現在と同様

10万年～100万年

RES-2:
プレートの収束速度が2倍に増加

10万年～
100万年

RES-3:
プレート運動の方向が変化

10万年～
100万年

RES-4:
プレートの沈み込み角度が増大

P.5

ステップ1b サイト変遷シナリオ ／ ステップ1c 影響シナリオ• 広域変遷シナリオの下で生じうるサイト周辺の地殻変動とそれに伴う事象の変遷を設定
サイト変遷シナリオ（SES） 東北地方日本海側の検討例• サイト変遷シナリオの下で生じうる事象の影響の規模と発生確率を設定
影響シナリオ(IS)SES2
地殻変動の速度が2倍に変化
SES3
圧縮変形が局所化し西に移動

事象による影響の頻度

SES1
地殻変動の傾向が現在と同様事象の確率空間分布（ITM手法）

処分場深度＝300mと仮定

P.6

火山活動の空間分布モデル• 火山の活動時期，位置，マグマの生成・移動に関連する地質学的・地球物理学的な情報を確率モデルに適用し，ある期間のある領域に火山活動が発生する確率を推定

確率 （log)
確率▲

第四紀単成火山●：第四紀火山¾ 活動時期と位置に基づく将来10万年間に ¾ 活動時期と位置に加えて重力と磁気データを用
いた将来10万年間に25km2の領域に１回以上の25km2 の領域に１回またはそれ以上の火山単成火山が発生する確率（コックスプロセス法：活動が生じる確率（カーネル法：東北地方へ中国地方への適用例）の適用例）

P.7ひずみ速度の空間分布モデル• 岩盤変形に関するデータに基づき将来のある期間までに蓄積されるひずみ速度を推定

ナノひずみ/年

ナノひずみ/年

ナノひずみ/年

地震

地表変位（断層変位 + 傾動）

GPS¾ 地震，地表変位，GPSのデータに基づく将来10万年間に25km2の領域に蓄積されるひずみ速度の空間分布（東北地方への適用例）
P.8

地表接近ハザードの空間分布モデル• 地表面の侵食により廃棄物が地表接近あるいは露出する可能性を相対的に評価
隆起速度（mm/年）

傾動速度（mm/km/年）

岩種分布
第四紀堆積岩新生代堆積岩と火山岩

0.0061 x 106 - 0.01 x 106

中生代と古生代堆積岩，火山岩，変成岩，新生代～古生代深成岩¾ 隆起速度，傾動速度，岩種分布（侵食されやすさの指標）に基づき，隆起速度＝侵食速度と仮定した場合の地下300mの廃棄物が地表に接近するハザードの空間分布（東北地方への適用例）
P.9

ステップ２ ロジックツリーの構築
火山活動の空間分布モデル火山活動の規模と頻度のモデル

経路（シナリオ）ごとの特定規模の火山の発生確率

RES-1

RES-2
広域変遷シナリオ（RES）が起こる相対的な確信度

火山活動の空間分布を推定するモデルの相対的なRESに対応したサイト
確信度
変遷シナリオ（SES）が起こる相対的な確信度

火山活動の規模と頻度を推定するモデルの相対的な確信度

このシナリオの確率と空間分布モデ
ルに基づきハザー
ドマップを作成

火山活動の影響に関するシナリオのロジックツリーの例
P.10

ステップ3 専門家の意見集約（Expert Elicitation)（RES-2: プレート収束速度2倍）
平均値
最適値

Expert 1

専門家 12
専門家 11
Expert 2

専門家 10
専門家 9
専門家 8

Expert 3

専門家 7
専門家 6

専門家の意見の統計学的な意見集約手法(Cook, 1991)

専門家 5
専門家 4
専門家 3
専門家 2
専門家 1
0

0.1

0.2

確信度（0~1）

0.3

0.4

0.5
P.11

ステップ４ サイト評価への適用： 統合ハザードマップ
火山活動ハザードと地表接近ハザードを統合した
相対的ハザード

火山ハザードと岩盤変形ハザードを統合した相対的ハザード（パーセンタイル）高中● 第四紀火山- - - - 活断層

低○ 第四紀火山

火山活動と岩盤変形の統合ハザードマップ（九州地方への適用例）

火山活動と地表接近の統合ハザードマップ（東北地方への適用例）
P.12

ステップ４ 安全評価への適用： 影響シナリオの発生確率

火山フロント海溝側
火山周辺
隆起高地
活断層周辺

第四紀火山
○ カルデラ火山

超過確率

確率(log)

内陸盆地●¾過去100万年間のカルデラ火山の活動履歴と第四紀火山の分布に基づき，将来10万年間に25km2の領域に半径7kmの一つまたはそれ以上のカルデラ火山が生じる確率（九州地方への適用例）

ひずみ速度（ナノひずみ/年）¾火山や断層との位置関係に基づき設定した複数の領域のひずみ速度の超過確率の累積頻度分布（東北地方への適用例）P.13

まとめ検討成果：• 一連の手法を開発し，以下の各段階で活用できる見通しを得た¾事前確認～文献調査
9ハザードマップを用いた評価結果の検証・補完¾概要調査
9現地調査データを用いた詳細検討による評価結果の検証・補完¾精密調査～安全審査
9超長期の影響シナリオの発生確率の提示
残された課題：• 地層処分システムへの影響シナリオの具体化• 国内専門家の理解促進と協力体制の構築
今後の予定：
 成果の取りまとめ・公表
P.14