各計算タイプにおけるハザード・リスクの計算方法
OpenQuake-engine で計算可能な以下の各計算タイプについて、簡単なサンプルを用いて説明を行う。
ハザードScenario CaseClassical PSHAEvent-based PSHADisaggregation

リスク


Scenario Damage

Scenario Risk

Classical Damage

Classical Risk

Classical BCR

ここでは、種々のパラメータと計算方法の説明にとどめ、地震の発生確率等のデータを用いて実際に各種ハザードの計算を実施する例は「第 4 章
よるハザード・リスクの計算例」（p.185~）で説明する。

105

実データを用いた OpenQuake に3.6.1

ハザード

Scenario Case
まず、一番簡単な例として、断層モデルを与えて地震動予測式に基づき地震動を計算する例を示す。OpenQuake-engine では Scenario Case（Scenario-based Hazard）という計算である。震源断層を特定した地震動と考えてもよい。この解析は、地震の発生確率を考慮せず、地震が発生した場合に、地震動の強さを計算するものである。

1)

必要なファイル

確率論的地震ハザード解析とは異なり、地震動予測式で使用する断層モデルとOpenQuake-engine 実行のための設定ファイルの以下の 2 ファイルが必要となる。
設定ファイル（ファイル 3.6.1）
断層モデル（ファイル 3.6.2）
震源モデルロジックツリーファイルや地震動予測式ロジックツリーファイルは必要ない。
設定ファイルをファイル 3.6.1 に、断層モデルファイルをファイル 3.6.2 に示した。Scenario
Case の計算では、hazard_outputs セクションには、何も記入しなくても、地震動強さの分布（Ground Motion Field）のみ出力される。
ファイル 3.6.1
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26

Scenario Case で使用する設定ファイル

ファイル内容[general]description = Scenario Calculation with Simple Fault Rupture
calculation_mode = scenario[sites]region_grid_spacing = 2

説明
計算概要
モードの選択region で指定した範囲内での計算地
点の刻み(km)
計算領域をカンマ区切りで指定(°)

region = 80 26, 88 26, 88 31, 80 31[erf]
rupture_model_file = earthquake_rupture_model.xml
rupture_mesh_spacing = 5

断層モデル
震源距離計算時に断層モデルを分割するサイズ(km)[site_parames]
reference_vs30_value = 760.0
reference_vs30_type = measured
reference_depth_to_2pt5km_per_sec = 5.0
reference_depth_to_1pt0km_per_sec = 100.0

Vs30(m/s)
Vs30 のタイプ
Vs=2.5km/s となる深さ(km)Vs=1.0km/s となる深さ(m)[correlation]
ground_motion_correlation_model =地震動補正モデル[hazard_calculation]
random_seed = 113
gsim = ChiouYoungs2008
intensity_measure_types = PGA
truncation_level = 3.0

統計処理のための乱数のシード値
地震動予測式

計算する地震動強さ
切断正規分布有界レベル

106

27
28

ファイル内容maximum_distance = 200
number_of_ground_motion_fields = 10

説明
計算する最大距離
計算回数*斜体は選択した GSIM によって、使用の有無や変数が異なる
ファイル 3.6.2 Scenario Case で使用する断層モデル
1
2
3
4
5
6

ファイル内容<?xml version='1.0' encoding='utf-8'?><nrml xmlns:gml="http://www.opengis.net/gml"
xmlns="http://openquake.org/xmlns/nrml/0.5"><simpleFaultRupture><magnitude>7.0</magnitude><rake>90.0</rake><hypocenter lat="84.4" lon="27.6" depth="30"/>

7<simpleFaultGeometry>

8
9
10
11
12
13
14
15

<gml:LineString><gml:posList>
85.7 27.3
83.8 27.8</gml:posList></gml:LineString><dip>90.0</dip><upperSeismoDepth>20.0</upperSeismoDepth>

16<lowerSeismoDepth>50.0</lowerSeismoDepth>

17
18
19

</simpleFaultGeometry></simpleFaultRupture></nrml>

この例では、「3.3.2

(2)

1) ③

説明
XML 宣言nrml: NRML 開始タグ
xmlns: XML 名前空間
simpleFaultRupture : Simple fault
magnitude: マグニチュードは 7
rake: すべり角
hypocenter: 震源位置
lat: 緯度 84.4 度、lon: 経度 27.7 度、
depth: 深さ 30km
simpleFaultGeometry: 断層面の形状
を simple fault タイプを使用する
gml:LineString: 線を使用
gml:posList: 位置のリスト
始点の経度 85.7 度 緯度 27.3 度
始点の経度 83.8 度 緯度 27.8 度
gml:posList 終了タグ
gml:LineString 終了タグ
dip: 傾斜角 90 度
upperSeismoDepth: 地震発生上面深さ 20.0 (km)lowerSeismoDepth: 地震発生下面深さ 50.0 (km)simpleFaultGeometry 終了タグ
simpleFaultRupture 終了タグ
nrml 終了タグ

Simple fault rupture」（p.45）で説明した Simple fault

rupture による断層モデルと Chiou and Youngs (2008)43)による地震動予測式を用いて、地表
最大加速度を求める。Chiou and Youngs (2008)の地震動予測式では、観測点から断層までの距離のほかに、すべり角、傾斜角、断層モデルの上端深さ、マグニチュード、地盤の S 波
速度が 100m/s となる深さ、Vs30（地表 30m の平均 S 波速度）が必要となる。ファイル 3.6.2に示したとおり、断層モデル中には地震の発生確率のデータは含まない。
計算範囲あるいは計算地点を変更するには、[sites]セクション（ファイル 3.6.1 中 5~7 行目）を変更することになる。観測点数が少ない場合には、ファイル 3.6.2(a)に示したように、
sites 変数を用いて、“経度(°) 緯度(°)”の組み合せをカンマ区切りで記入する。また、ファイル 3.6.2(b)に示したように、csv 形式ファイルを使用して、計算位置を指定することも可能である。csv ファイルは、1 行に 1 計算地点の“経度(°),緯度(°)”のデータが複数行記述されたファイルになる。また、特殊なケースとして、観測点地盤特性のファイルに含まれる
位置情報で指定する方法（site_mode_file 変数）や、exposure モデルファイルに含まれる各107

資産の位置情報をそのまま使用する（exposure_file 変数）ことも可能であるが、ここでは、説明を省略する。なお、その際には region 変数および region_grid_spacing 変数は不要である。
ファイル 3.6.2(a) 数が少ない場合の観測点指定方法
5
6

ファイル内容[Hazard sites]
sites = 85 26, 87 27, ...

説明
数が少ない場合には経度, 緯度を指定することもできる

7

ファイル 3.6.2(b)56csv 形式ファイルを使用した観測点指定方法

ファイル内容[Hazard sites]
sites_csv = site.csv

説明
csv 形式ファイルを使用した観測点指定方法

7

計算に使用する地震動予測式を変更する場合は、gsim 変数（ファイル 3.6.1 中 24 行目）
を変更する。使用する gsim 変数に合わせて intensity_measure_types 変数（ファイル 3.6.1 中
25 行目）を変更することで、計算する地震動指標を指定する。intensity_measure_types 変数
は、使用する地震動予測式に対応する地震動強さを指定する必要がある。地震動予測式が対応していない場合は、もちろん計算ができない。intensity_measure_types 変数で指定する
地震動指標は、複数指定することも可能である。Chiou and Youngs (2008)による地震動予測
式では、地表最大加速度（PGA）、地表最大速度（PGV）、加速度応答スペクトル（SA）が計算可能である。各種スペクトルの計算においては、対象とする周期の情報を、括弧付きで記述する必要がある（ファイル 3.6.2(c)参照）。
ファイル 3.6.2(c) スペクトル計算時の周期の記述方法
25

ファイル内容intensity_measure_types = PGA, SA(0.05), SA(0.1)

説明
計算する地震動強さ

OpenQuake-engine による地震動の計算では、地震動予測式に付属する標準偏差のデータを用いて、正規分布に従う乱数を使用した地震動のばらつきが適用され、偶然的不確定性
として空間補正がなされ、統計データとして使用される。random_seed 変数（ファイル 3.6.1
中 23 行目）はそのための乱数のシード値であり、number_of_ground_motion_fields 変数（ファイル 3.6.1 中 28 行目）はそのための計算回数である。number_of_ground_motion_fields 変数を 2 以上に設定すると、その回数の計算がなされる。偶然的不確定性を考慮しない場合
は、truncation_level 変数（ファイル 3.6.1 中 26 行目）を 0 にする必要がある。この場合、複数回処理しても地震動強さは同じであるため、number_of_ground_motion_fields 変数は 1108でよい。この truncation_level 変数そのものを記述しない場合は標準正規分布に従う乱数を使用して空間補正が適用され、0 でない場合は切断正規分布に従う乱数を使用して空間補正が適用される。ここで、truncation_level 変数を 0 にした場合、すなわち、偶然的不確定性を考慮しない場合は、ground_motion_correlation_model 変数（ファイル 3.6.1 中 20 行目）
で指定する地震動補正モデルを使用することもできないことに注意する。偶然的不確定性に関するフローチャートを図 3.6.1 に示した。

ハザードを計算する

偶然的不確定性を考慮するNo

Yes

切断正規分布
を使用する

No

Yes
truncation_level=3
有界レベルを設定する

truncation_level 変 数 を
Yes
記述しない

truncation_level=0
有界レベルを 0 にする

計算実行

図 3.6.1

地震動予測式に付随する標準偏差と乱数を使用した空間補正により偶然的不確定性を考慮するためのフローチャート

2)

Input Preparation Toolkit によるファイルの準備

Input Preparation Toolkit を使用して「1) 必要なファイル」（p.106~）で説明した 2 ファイ
ルを作成する方法を紹介する。Input Preparation Toolkit の使用方法は、「3.5

Input Preparation

Toolkit」（p.100~）で説明したとおりである。はじめに「Earthquake Rupture」タブで震源モ
デルファイルを作成し、次に「Configuration file」タブで設定ファイルを作成する。

①

震源ファイル

震源ファイルは以下の手順で作成する。
Input Preparation Toolkit の「Earthquake Rupture」タブを開く（図 3.6.2）。
必要な情報を入力する（図 3.6.3）。
ここでは、断層モデルを 2 点間のトレースと傾斜角で表現するため、Rupture type109

として「Simple Fault Rupture」を選択する。任意の断層面の場合は「Arbitrary Fault
Rupture」、複数の任意の断層面の場合は「Single Planar Rupture OR Multi-Planar
Rupture」、断層面の上端トレース・下端トレースを指定する場合には「Complex Fault
Rupture」を選択し、必要な情報を入力する。
「Convert NRML」ボタンをクリックして、nrml 形式に変換する（図 3.6.4）。
「Download」ボタンをクリックして、作成された震源モデルファイルをダウンロー
ドする（図 3.6.5）。

図 3.6.2

Input Preparation Toolkit（断層モデル作成）その 1
Earthquake Rupture タブを開いた画面

110

図 3.6.3

Input Preparation Toolkit（断層モデル作成）その 2
必要な情報の入力画面する

図 3.6.4

Input Preparation Toolkit（断層モデル作成）その 3
nrml 形式に変換した状態

111

図 3.6.5

Input Preparation Toolkit（断層モデル作成）その 4作成された震源モデルファイルのダウンロードする②

設定ファイル

設定ファイルは以下の手順で作成する。
Input Preparation Toolkit の「Configuration file」タブを開く（図 3.6.6）。ここでは、Scenario Case の計算なので、「Earthquake Scenarios」タブを開く。
ハザードの計算を行うので、Hazard にチェックを入れる（図 3.6.7）。
「Rupture model file」で参照ボタンをクリックして「①震源ファイル」で作成した earthquake_rupture_model.xml をアップロードする（図 3.6.8）。
earthquake_rupture_model.xml をリスクとから選択し、その他必要事項を記入する（図3.6.9、図 3.6.10）。
「Download」ボタンを押すと、必要なファイルが含まれた圧縮ファイルが作成される（図 3.6.11）。

112

図 3.6.6

Input Preparation Toolkit（ハザード設定ファイル）その 1「Configuration file」タブを開き、さらに「Earthquake Scenarios」タブを開く

図 3.6.7

Input Preparation Toolkit（ハザード設定ファイル）その 2
Hazard にチェックを入れる

113

図 3.6.8

Input Preparation Toolkit（ハザード設定ファイル）その 3

参照ボタンをクリックして「①震源ファイル」で作成した断層モデルファイルをアップロードする

図 3.6.9

Input Preparation Toolkit（ハザード設定ファイル）その 4

earthquake_rupture_model.xml を選択し、その他必要事項を記入する

114

図 3.6.10 Input Preparation Toolkit（ハザード設定ファイル）その 5
必要事項を選択・記入する

図 3.6.11 Input Preparation Toolkit（ハザード設定ファイル）その 6
ダウンロードボタンを押すと、必要なファイルが含まれた圧縮ファイルが作成できる

3)

計算方法

計算方法は「3.4.1

WebUI による実行」（p.81~）および「3.4.2 コマンドラインによる計算実行方法」（p.90~）に示したとおりであるが、確認のため、WebUI を使用した計算方法を説明する。
実行はウェブブラウザから「http:://localhost:8800/engine」にアクセスして行う（図3.6.12）。
「Run a Calculation」をクリックすると、ファイル選択画面が開かれるので、必要なファイルが含まれた圧縮ファイルを選択する（図 3.6.13）。

115

ファイルを選択すると計算が開始される（図 3.6.14）。
計算が終了すると、「Status」が「Complete」に変更され、「Outputs」ボタンと「RunRisk」ボタンが表示される（図 3.6.15）。
「Console」をクリックすると実行ログを表示できる（図 3.6.16）。

図 3.6.12 WebUI による計算実行その 1
ウェブブラウザから「http://localhost:8800/engine」を開く

116

図 3.6.13 WebUI による計算実行その 2
実行したい設定ファイルおよび震源モデルが含まれた圧縮ファイル(ZIP 形式)を選択してファイルをアップロードする

図 3.6.14 WebUI による計算実行その 3
ファイルをアップロードすると自動的に計算が始まる

117

図 3.6.15 WebUI による計算実行その 4計算が終了すると、「Status」が「complete」になり、Outputs、Run Risk ボタンが表示される

図 3.6.16 WebUI による計算実行その 5「Console」から開いた実行ログを表示すると、この計算では、108546 地点の地震動を計算し、計算にかかった時間は 71 秒であることがわかる

118

4)

計算結果

Scenario Case による計算結果は Ground Motion Field（地震動強さの分布）のデータのみである。計算結果は nrml 形式、csv 形式および npz 形式（NumPy 配列形式）の 3 パターンで出力可能である。また、各種データへ変換前の形式は hdf5 形式となっている。計算結果の出力形式一覧を表 3.6.1 に示した。
表 3.6.1

Scenario Case の計算結果の出力形式一覧

計算結果

名前

出力可能な形式

地震動強さの分布

gmf_data

nrml 形式
csv 形式
npz 形式*太字はデフォルトで出力される
nrml 形式のファイルは、OpenQuake Platform や GEM の開発者が作成したその他のツールで使用可能であるが、汎用のソフトウェアでは単純に使用することはできない。csv 形式ファイルは、各種 GIS ソフトウェアや表計算ソフトで簡単に開くことができる。npz 形
式は、Python を使用してデータを処理する場合以外には使用しない。
WebUI を使用しデータの出力方法は、図 3.6.17~図 3.6.19 に示したとおりである。

図 3.6.17 WebUI による計算実行その 6
計算が終了した際に「Outputs」ボタンが表示される

119

図 3.6.18 WebUI による計算実行その 7「Outputs」ボタンをクリックすると、この計算で出力されたファイルに一覧が表示される

図 3.6.19 WebUI による計算実行その 8
データをダウンロードする
複数の計算結果がある場合には、それぞれの計算結果が 1 つに圧縮されたファイルとなっている。図 3.6.19 に示したように出力ファイル名は、計算 ID、地震動予測式、対応する

120

地震動強さおよび出力 ID などを使用して自動的に命名されている。
nrml 形式の出力サンプルをファイル 3.6.3 に、csv 形式の出力サンプルをファイル 3.6.4に示した。参考までに、偶然的不確定性を考慮しない場合(truncation_level = 0)の計算結果
をファイル 3.6.5 に示した。偶然的不確定性を考慮しない場合は、10 回の計算がすべて同じ値になっていることがわかる。
前述したとおり、csv 形式ファイルはほかのソフトウェアを用いた図化が可能であるが、nrml 形式のファイルについては、OpenQuake 開発者が作成したツール群を使用して図化を行う必要がある。ただし、現時点（2017 年 3 月時点）では gmf_data が記述された nrml 形
式をそのまま図化する方法は開発者からは提供されていない。
ファイル 3.6.3 output-1-gmf-rlz-000_1.xml サンプル
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19

ファイル内容<?xml version="1.0" encoding="utf-8"?><nrml
xmlns="http://openquake.org/xmlns/nrml/0.5"
xmlns:gml="http://www.opengis.net/gml"><gmfCollection
gsimTreePath="b1"
sourceModelTreePath="b_1"

説明
XML 宣言nrml 開始タグ
xmlns: XML 名前空間
xmlns: XML 名前空間
gmfCollection: gmf コレクション
gsimTreePath: gsim ツリーパス
sourceModelTreePath: 震源モデルツ
リーパス><gmfSetstochasticEventSetId="1"

gmfSet: gmf セット
stochasticEventSetId: 統計イベント
セット 1

>
<gmf
IMT="PGA"
ruptureId="scenario-0000000000~ses=1"><node
gmv="4.6253479E-03"
lat="2.6024265E+01" lon="8.0050040E+01"/><node
gmv="4.6253479E-03"
lat="2.6024265E+01" lon="8.0050040E+01"/><node
gmv="4.6253479E-03"
lat="2.6024265E+01" lon="8.0050040E+01"/>

gmf: ground motion field
IMT : PGA データ
ruptureId: 断層モデル ID
node: ノード
gmv: 地震動強さ、lat: 緯度、lon: 経度ファイル 3.6.4 output-1-ChiouYoungs2008()-PGA_1.csv
1
23

ファイル内容lon,lat,000,001,002,003,004,005,006,007,008,009
80.02001,26.00628,0.15265,0.15293,0.30426,0.18915,0.20476,0.16088,0.13338,0.15385,0.16175,0.15885
80.02002,26.02427,0.03401,0.01250,0.01378,0.01210,0.02689,0.01106,0.02466,0.02735,0.04120,0.02169

ファイル 3.6.5 truncation_level = 0 のサンプル
1
2
3

ファイル内容lon,lat,000,001,002,003,004,005,006,007,008,009
80.02001,26.00628,0.21676,0.21676,0.21676,0.21676,0.21676,0.21676,0.21676,0.21676,0.21676,0.21676
80.02002,26.02427,0.02351,0.02351,0.02351,0.02351,0.02351,0.02351,0.02351,0.02351,0.02351,0.02351

121Classical PSHA（Point source）
次に、Sources for modelling distributed seismicity の震源モデルを使用した Classical PSHAの計算例を説明する。Sources for modelling distributed seismicity の震源モデルのうち、もっとも単純な Point Source を使用した計算例を示す。Classical PSHA は、一般的に確率論的地
震ハザード解析に使用されるものであり、ハザードカーブを直接的に計算する。

1)

必要なファイル

Scenario Case とは異なり、OpenQuake-engine 実行のための設定ファイルと震源モデルロ
ジックツリーおよび地震動モデルロジックツリーが必要となる。
設定ファイル（ファイル 3.6.6）
震源モデルロジックツリーファイル（ファイル 3.6.7）
地震動モデルロジックツリーファイル（ファイル 3.6.8）
震源モデルロジックツリー内で指定される震源モデル（ファイル 3.6.9）もまた必要となる。
設定ファイルをファイル 3.6.6、震源モデルロジックツリーファイルをファイル 3.6.7、地震動モデルロジックツリーファイルをファイル 3.6.8、震源モデルロジックツリー内で指定される震源モデルをファイル 3.6.9 にそれぞれ示す。
ファイル 3.6.6 Classical PSHA（Point source）で使用する設定ファイル
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27

ファイル内容[general]description = Classical PSHA Calculation with Point Source
calculation_mode = classical

説明
計算概要
モードの選択[sites]region = 80 26, 88 26, 88 31, 80 31
region_grid_spacing = 5

計算領域をカンマ区切りで指定(度)region 範囲内の計算地点の刻み(km)[logic_tree]source_model_logic_tree_file = source_model_logic_tree.xml
gsim_logic_tree_file = gmpe_logic_tree.xml

震源モデルロジックツリー
地震動予測式ロジックツリー[erf]width_of_mfd_bin = 0.1
rupture_mesh_spacing = 2

マグニチュード別度数分布の刻み
震源距離計算時に断層モデルを分割するサイズ(km)[site_params]
reference_vs30_value = 760.0
reference_vs30_type = measured
reference_depth_to_2pt5km_per_sec = 5.0
reference_depth_to_1pt0km_per_sec = 100.0

Vs30(m/s)
Vs30 のタイプ
Vs=2.5km/s となる深さ(km)Vs=1.0km/s となる深さ(m)[correlation]
ground_motion_correlation_model =地震動補正モデル[hazard_calculation]
random_seed = 113

乱数のシード値

122

28

29
30
31
32
33
34
35
36

ファイル内容intensity_measure_types_and_levels = {"PGA": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(0.025)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13], "SA(0.05)": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(0.1)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13], "SA(0.2)": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(0.5)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13], "SA(1.0)": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(2.0)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13]}
truncation_level = 3.0
maximum_distance = 200
investigation_time = 50.0

説明

切断正規分布の有界レベル計算する最大距離
発生確率の対象期間[hazard_outputs]
hazard_maps = true
uniform_hazard_spectra = true
poes = 0.1 0.02

ハザードマップの出力
uhs の出力
超過確率 10%、2%

計算する地震動指標とその強さ
複数選択時には、{ }、[]、””、カンマなどの区切り文字を使用して適切に記述する必要がある*斜体は選択した GSIM によって、使用の有無や変数が異なる
ファイル 3.6.7 Classical PSHA（Point source）で使用する震源モデルロジックツリー
ファイル
1
2

ファイル内容<?xml version="1.0" encoding="UTF-8"?><nrml xmlns="http://openquake.org/xmlns/nrml/0.5">3<logicTree logicTreeID="lt1">

4<logicTreeBranchingLevel branchingLevelID="bl1">

5<logicTreeBranchSet uncertaintyType="sourceModel"
branchSetID="bs1">

6<logicTreeBranch branchID="b1">

7
8
9
10
11
12
13

<uncertaintyModel>source_model.xml</uncertaintyModel><uncertaintyWeight>1.0</uncertaintyWeight></logicTreeBranch></logicTreeBranchSet></logicTreeBranchingLevel></logicTree></nrml>123

説明XML 宣言nrml: NRML 開始タグ
xmlns: XML 名前空間
logicTree: ロジックツリー
logicTreeID: ID"lt1"
logicTreeBranchingLevel: ブランチ
ングレベル
branchingLevelID: ID"bl1"
logicTreeBranchSet: ブランチセット
branchSetID: ID"bs1、uncertaintyType:sourceModel
logicTreeBranch: ブランチ
branchID: ID b1"
uncertaintyModel: 不確定性は
source_model.xml
uncertaintyWeight: 重みは 1.0
logicTreeBranch 終了タグ
logicTreeBranchSet 終了タグ
logicTreeBranchingLevel 終了タグ
logicTree 終了タグ
nrml 終了タグ

ファイル 3.6.8 Classical PSHA（Point source）で使用する地震動予測式ロジック
ツリーファイル
1
2

ファイル内容<?xml version="1.0" encoding="UTF-8"?><nrml xmlns="http://openquake.org/xmlns/nrml/0.5">3<logicTree logicTreeID='lt1'>

4<logicTreeBranchingLevel branchingLevelID="bl1">

5<logicTreeBranchSet uncertaintyType="gmpeModel"
branchSetID="bs1" applyToTectonicRegionType="Active
Shallow Crust">

6<logicTreeBranch branchID="b1">

7
8
9
10
11
12
13

<uncertaintyModel>AkkarBommer2010</uncertaintyModel><uncertaintyWeight>1.0</uncertaintyWeight></logicTreeBranch></logicTreeBranchSet></logicTreeBranchingLevel></logicTree></nrml>

説明XML 宣言nrml: NRML 開始タグ
xmlns: XML 名前空間
logicTree: ロジックツリー
logicTreeID: ID"lt1"
logicTreeBranchingLevel: ブランチ
ングレベル
branchingLevelID: ID"bl1"
logicTreeBranchSet: ブランチセット
branchSetID: ID"bs1、uncertaintyType:
gmpeModel、
applyToTectonicRegionType: 適用するテクトニックタイプは Active
Shallow Crust
logicTreeBranch: ブランチ
branchID: ID b1"
uncertaintyModel: 不確定性は
AkkarBommer2010
uncertaintyWeight: 重みは 1.0
logicTreeBranch 終了タグ
logicTreeBranchSet 終了タグ
logicTreeBranchingLevel 終了タグ
logicTree 終了タグ
nrml 終了タグ

ファイル 3.6.9 Classical PSHA（Point source）で使用する震源モデルファイル
1
2
3
4
5

ファイル内容<?xml version='1.0' encoding='utf-8'?><nrml xmlns:gml="http://www.opengis.net/gml"
xmlns="http://openquake.org/xmlns/nrml/0.5"><sourceModel name="Example Source Model Containing a
Point Source"><sourceGroup name="g1" tectonicRegion="Active Shallow
Crust"><pointSource id="2" name="Point Source"
tectonicRegion="Active Shallow Crust">

6
7
8
9
10

<pointGeometry><gml:Point><gml:pos>84.4 27.6</gml:pos></gml:Point><upperSeismoDepth>20.0</upperSeismoDepth>11<lowerSeismoDepth>50.0</lowerSeismoDepth>

1213</pointGeometry><magScaleRel>WC1994</magScaleRel>

14<ruptAspectRatio>1.5</ruptAspectRatio>124

説明
xml 宣言nrml: NRML 開始タグ
xmlns: XML 名前空間
sourceModel: 震源モデルファイル
name: 簡単な説明
sourceGroup: 震源グループ
tectonicRegion: テクトニックタイプ
pointSource : Point Source を使用
name: 簡単な説明、tectonicRegion:
テクトニックタイプは Active
Shallow Crust
pointGeometry: Point ジオメトリ
gml:Point: 点
gml:pos:経度 84.4 度、緯度 27.6 度
gml:Point 終了タグ
upperSeismoDepth: 地震発生上面深さ 20(km)lowerSeismoDepth: 地震発生下面深さ 50(km)pointGeometry 終了タグ
magScaleRel: スケーリング則はWC1994
ruptAspectRatio: 断層面のアスペク
ト比は 1.5

15

16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33

ファイル内容<truncGutenbergRichterMFD aValue="3.0" bValue="1.0"minMag="5.0" maxMag="7.0" /><nodalPlaneDist><nodalPlane probability="0.2" strike="265.0" dip="90.0"
rake="90.0" /><nodalPlane probability="0.2" strike="275.0" dip="90.0"
rake="90.0" /><nodalPlane probability="0.2" strike="285.0" dip="90.0"
rake="90.0" /><nodalPlane probability="0.2" strike="295.0" dip="90.0"
rake="90.0" /><nodalPlane probability="0.2" strike="305.0" dip="90.0"
rake="90.0" /></nodalPlaneDist><hypoDepthDist><hypoDepth probability="0.2" depth="25.0" /><hypoDepth probability="0.2" depth="30.0" /><hypoDepth probability="0.2" depth="35.0" /><hypoDepth probability="0.2" depth="40.0" /><hypoDepth probability="0.2" depth="45.0" /></hypoDepthDist></pointSource></sourceGroup></sourceModel></nrml>説明
truncGutenbergRichterMFD: G-R 則に基づくマグニチュード別度数分布
aValue: 3.0、bValue: 1.0、minMag: 5.0、
maxMag: 7.0
nodalPlaneDist: 断層面のばらつき
nodalPlane: 断層面
probability: 発生確率、strike: 走向、
dip: 傾斜、rake: すべり角
傾斜 90 度、すべり角 90 度は固定で、走向が 265 度、275 度、285 度、295
度、305 度の断層面の発生確率が各20%

nodalPlaneDist 終了タグ
hypoDepthDist: 震源深さのばらつき
hypoDepth: 震源深さ
probability: 発生確率、depth: 深さ
深さ 25km、30km、35km、40km、45kmの発生確率が各 20%
hypoDepthDist 終了タグ
pointSource 終了タグ
sourceGroup 終了タグ
sourceModel 終了タグ
nrml 終了タグ

この例では、
ファイル 3.6.6 で示した設定ファイルと、ファイル 3.6.7 およびファイル 3.6.8
で示した各ロジックツリーファイル、ファイル 3.6.9 で示した震源モデルに従い、確率論的
地震ハザード解析を行う。
震源モデルロジックツリーファイル（ファイル 3.6.7）中に記述されるように震源モデル
はファイル 3.6.9 の 1 つのみであるが、複数の震源モデルでロジックツリー解析を実施する
場合は、ファイル 3.6.7(a)に示したように、ファイル 3.6.7 の 6~9 行目のようなブランチをもう一つ作成し、重みを足して 1.0 となるように変更する。複数の震源モデルを組み込む
計算は、「4.5 確率論的地震ハザード解析その 3」（p.210~）で紹介するファイル 3.6.7(a) 複数の震源モデルを使用する震源モデルロジックツリーファイル
5

6

ファイル内容<logicTreeBranchSet uncertaintyType="sourceModel"
branchSetID="bs1"><logicTreeBranch branchID="b1">

7
8
9
10

<uncertaintyModel>source_model_1.xml</uncertaintyModel><uncertaintyWeight>0.5</uncertaintyWeight></logicTreeBranch><logicTreeBranch branchID="b2">125

説明
logicTreeBranchSet: ブランチセット
branchSetID: ID"bs1、uncertaintyType:sourceModel
logicTreeBranch: ブランチ
branchID: ID b1"
uncertaintyModel: 不確定性は
source_model_1.xml
uncertaintyWeight: 重みは 0.5
logicTreeBranch 終了タグ
logicTreeBranch: ブランチ
branchID: ID b1"

ファイル内容
11
12
13
14

<uncertaintyModel>source_model_2.xml</uncertaintyModel><uncertaintyWeight>0.5</uncertaintyWeight></logicTreeBranch></logicTreeBranchSet>

説明
uncertaintyModel: 不確定性は
source_model_2.xml
uncertaintyWeight: 重みは 0.5
logicTreeBranch 終了タグ
logicTreeBranchSet 終了タグ

簡単化のため、今回の計算では、震源モデル中の断層モデルも 1 つのみとしているが、必要に応じて断層モデルを追加することができる。注意点として、ファイル 3.6.9 の震源モ
デル中に指定される tectonicRegion に対応する地震動予測式がファイル 3.6.8 の地震動モデ
ルのロジックツリー中に存在する必要がある。
Classical PSHA の計算では、設定ファイルにおいて、mean_hazard_curves 変数およびquantile_hazard_curves 変数の設定もできるが、ロジックツリーによる経路が 1 通りである
今回の計算では意味がないため、この設定ファイル中では設定していない。
多くの設定は Scenario Case の計算時と同様であるが、source_model_logic_tree_file 変数（ファイル 3.6.6 中 10 行目）、gsim_logic_tree_file 変数（同 11 行目）、計算する地震動指標とその強さを記述する intensity_measure_types_and_levels 変数（同 28 行目）、超過確率を計算する期間である investigation_time 変数（同 31 行目）および計算したい超過確率を指定す
る poes 変数（同 36 行目）を適切に記述する必要がある。
この例では、断層面やマグニチュードとは独立の震源深さのばらつきを複数設定しているが（ファイル 3.6.9 中 16~29 行）、一意に決まる場合は、probability=1.0 にして、いずれか一つだけ設定すればよい。
Classical PSHA の計算では、ハザードカーブはデフォルトで出力される。ハザードマッ
プ や 、 一 様 ハ ザ ー ド ス ペ ク ト ル を 出 力 し た い 場 合 に は 、 hazard_maps 変 数 や
uniform_hazard_spectra 変数（ファイル 3.6.6 中 34、35 行目）を true にする必要がある。

2)

計算方法

計算方法は「3.4.1

WebUI による実行」（p.81~）および「3.4.2 コマンドラインによる計算実行方法」（p.90~）に示したとおりである。

3)

計算結果

Classical PSHA による計算では、設定ファイル中のパラメータの指定およびロジックツ
リーでの経路によって複数の計算結果が出力される。今回の計算では、ロジックツリーは1 経路であることから、ハザードカーブ、ハザードマップおよび一様ハザードスペクトルのみが計算可能である。計算結果は出力の種類ごとに、csv 形式、nrml 形式、geojson 形式
および npz 形式の 4 パターンで出力可能である（表 3.6.2）。ロジックツリーパスと震源グ
ループについては、csv 形式で出力される。

126

表 3.6.2

Classical PSHA（Point source）の計算結果の出力形式一覧

計算結果

名前

出力形式

ハザードカーブhcurves

ハザードマップ

hmap

一様ハザードスペクトルuhs

nrml 形式
csv 形式
geojson 形式
npz 形式

ロジックツリーパス

realizations

震源グループsourcegroups

csv 形式*太字はデフォルトで出力される
csv 形式ファイルは汎用ソフトウェアを使用して図化を行う。nrml 形式ファイルはOpenQuake 開発者が作成したツールを使用して図化することが可能である。OpenQuake 開
発者が作成したツールを用いて図化した計算結果の例を図 3.6.20 に示す。

図 3.6.20 出力結果の例（Classical PSHA（Point source））
左: ハザードカーブ、中: 一様ハザードスペクトル、右: ハザードマップ

127Classical PSHA（Simple fault source）
次に、Fault sources with floating ruptures の震源モデルによる Classical PSHA の計算例を説明する。Fault sources with floating ruptures のうち、もっとも単純な Simple fault source を使
用する。

1)

必要なファイルClassical PSHA（Point source）と同様に、OpenQuake-engine 実行のための設定ファイルと震源モデルロジックツリーおよび地震動モデルロジックツリーが必要となる。
設定ファイル（ファイル 3.6.6）
震源モデルロジックツリーファイル（ファイル 3.6.7）
地震動モデルロジックツリーファイル（ファイル 3.6.8）
震源モデルロジックツリー内で指定される震源モデル（ファイル 3.6.10）もまた必要となる。ここでは、設定ファイル、震源モデルロジックツリーファイルおよび地震動モデルロ
ジックツリーファイルは震源モデルロジックツリーおよび地震動モデルロジックツリーは、「(2) Classical PSHA（Point source）」（p.122~）と同じファイル（ファイル 3.6.6、ファイル 3.6.7 およびファイル 3.6.8）を使用する。震源モデルロジックツリー内で指定される震源モデルをファイル 3.6.10 に示す。
ファイル 3.6.10

Classical PSHA（Simple fault source）で使用する震源モデル

ファイル
1
2
3
4
5

ファイル内容<?xml version='1.0' encoding='utf-8'?><nrml xmlns:gml="http://www.opengis.net/gml"
xmlns="http://openquake.org/xmlns/nrml/0.5"><sourceModel name="Example Source Model Containing a
Simple Fault Source"><sourceGroup name="g1" tectonicRegion="Active Shallow
Crust"><simpleFaultSource id="3" name="Simple Fault Source"
tectonicRegion="Active Shallow Crust">

6
7
8
9
10
11
12
13
14

<simpleFaultGeometry><gml:LineString><gml:posList>
85.7 27.3
83.8 27.8</gml:posList></gml:LineString><dip>90.0</dip><upperSeismoDepth>20.0</upperSeismoDepth>

15<lowerSeismoDepth>50.0</lowerSeismoDepth>

16</simpleFaultGeometry>128

説明
xml 宣言nrml: NRML 開始タグ
xmlns: XML 名前空間
sourceModel: 震源モデルファイル
name: 簡単な説明
sourceGroup: 震源グループ
tectonicRegion: テクトニックタイプ
simpleFaultSource : Simple fault
name: 簡単な説明、tectonicRegion:
テクトニックタイプは Active
Shallow CrustsimpleFaultGeometry: Simple fault
gml:Point: 線
gml:posList: 位置情報のリスト
始点の経度 85.7 度、緯度 27.3 度
終点の経度 83.8 度、緯度 27.8 度
gml:posList 終了タグ
gml:LineString 終了タグ
dip: 傾斜角 90 度
upperSeismoDepth: 地震発生上面深さ 20kmlowerSeismoDepth: 地震発生下面深さ 50kmsimpleFaultGeometry 終了タグ

17

ファイル内容<magScaleRel>WC1994</magScaleRel>18<ruptAspectRatio>1.5</ruptAspectRatio>

19<truncGutenbergRichterMFD aValue="3.0" bValue="1.0"minMag="5.0" maxMag="7.0" />

20
21
22
23
24

<rake>90.0</rake></simpleFaultSource></sourceGroup></sourceModel></nrml>説明
magScaleRel: スケーリング則はWC1994
ruptAspectRatio: 断層面のアスペク
ト比は 1.5
truncGutenbergRichterMFD: G-R 則に基づくマグニチュード別度数分布
aValue: 3.0、bValue: 1.0、minMag: 5.0、
maxMag: 7.0
rake: すべり角 90 度
simpleFaultSource 終了タグ
sourceGroup 終了タグ
sourceModel 終了タグ
nrml 終了タグ

この例では、
ファイル 3.6.6 に示した設定ファイルと、ファイル 3.6.7 およびファイル 3.6.8に示した各ロジックツリーファイルに従い、確率論的地震ハザード解析を行う。多くの設定は Classical PSHA（Point source）と同様であり、震源モデルファイルのみ異なる。

2)

計算方法

計算方法は「3.4.1

WebUI による実行」（p.81~）および「3.4.2 コマンドラインによる計算実行方法」（p.90~）に示したとおりである。

3)

計算結果

Classical PSHA による計算結果は、構成ファイル中のパラメータの指定およびロジック
ツリーでの経路によって複数の計算結果が出力される。今回の計算でも、ロジックツリー
は 1 経路であることから、ハザードカーブ、ハザードマップおよび一様ハザードスペクト
ルのみが計算可能である。計算結果は出力の種類ごとに、csv 形式、nrml 形式、geojson 形式および npz 形式の 4 パターンで出力可能である（表 3.6.3）。
表 3.6.3 出力結果一覧
計算結果

名前

出力形式

ハザードカーブhcurves

ハザードマップ

hmap

一様ハザードスペクトルuhs

nrml 形式
csv 形式
geojson 形式
npz 形式

ロジックツリーパス

realizations

震源グループsourcegroups

csv 形式*太字はデフォルトで出力される
csv 形式ファイルは汎用ソフトウェアを使用して図化を行う。nrml 形式ファイルはOpenQuake 開発者が作成したツールを使用して図化することが可能である。OpenQuake 開
発者が作成したツールを用いて図化した計算結果の例を図 3.6.21 に示す。
129

図 3.6.21 出力結果の例（Classical PSHA（Simple fault source））
左: ハザードカーブ、中: 一様ハザードスペクトル、右: ハザードマップ
参考のため、「(2) Classical PSHA（Point source）」（p.122~）において、Point Source を使用
して、断層面を 1 面（ファイル 3.6.9 中 15~28 行を変更）とした場合（ファイル 3.6.11）の計算結果を図 3.6.22 に示した。設定ファイルは同じものを使用し、断層面の設定はほぼ同
じような設定としているが、震源モデルの設定方法の違いで、計算結果が大きく異なることが確認できる。
ファイル 3.6.11(2) Classical PSHA（Point source）の震源モデルのうち、断層面を1 面とした震源モデル（ファイル 3.6.9 改）
1
2
3
4
5

ファイル内容<?xml version='1.0' encoding='utf-8'?><nrml xmlns:gml="http://www.opengis.net/gml"
xmlns="http://openquake.org/xmlns/nrml/0.5"><sourceModel name="Example Source Model Containing a
Point Source"><sourceGroup name="g1" tectonicRegion="Active Shallow
Crust"><pointSource id="2" name="Point Source"
tectonicRegion="Active Shallow Crust">

6
7
8<pointGeometry><gml:Point><gml:pos>84.4 27.6</gml:pos>

910</gml:Point><upperSeismoDepth>20.0</upperSeismoDepth>11<lowerSeismoDepth>50.0</lowerSeismoDepth>

1213</pointGeometry><magScaleRel>WC1994</magScaleRel>

14<ruptAspectRatio>1.5</ruptAspectRatio>130

説明
xml 宣言nrml: NRML 開始タグ
xmlns: XML 名前空間
sourceModel: 震源モデルファイル
name: 簡単な説明
sourceGroup: 震源グループ
tectonicRegion: テクトニックタイプ
pointSource : Point source
name: 簡単な説明、tectonicRegion:
テクトニックタイプは Active
Shallow Crust
pointGeometry: Point
gml:Point: 点
gml:pos: 震央経度 84.4 度、緯度 27.6
度
gml:Point 終了タグ
upperSeismoDepth: 地震発生上面深さ 20．0(km)lowerSeismoDepth: 地震発生下面深さ 50.0(km)pointGeometry 終了タグ
magScaleRel: スケーリング則はWC1994
ruptAspectRatio: 断層面のアスペク
ト比は 1.5

15

16
17

18
19
20

21
22
23
24
25

ファイル内容<truncGutenbergRichterMFD aValue="3.0" bValue="1.0"minMag="5.0" maxMag="7.0" /><nodalPlaneDist><nodalPlane probability="1.0" strike="285.0" dip="90.0"
rake="90.0" /></nodalPlaneDist><hypoDepthDist><hypoDepth probability="1.0" depth="25.0" /></hypoDepthDist></pointSource></sourceGroup></sourceModel></nrml>説明
truncGutenbergRichterMFD: G-R 則に基づくマグニチュード別度数分布
aValue: 3.0、bValue: 1.0、minMag: 5.0、
maxMag: 7.0
nodalPlaneDist: 断層面のばらつき
nodalPlane: 断層面
probability: 発生確率、strike: 走向、
dip: 傾斜、rake: すべり角
傾斜 90 度、すべり角 90 度、走向 285
度の断層面の発生確率が 100%
nodalPlaneDist 終了タグ
hypoDepthDist: 震源深さのばらつき
hypoDepth: 震源深さ
probability: 発生確率、depth: 深さ
深さ 35km の発生確率が各 100%
hypoDepthDist 終了タグ
pointSource 終了タグ
sourceGroup 終了タグ
sourceModel 終了タグ
nrml 終了タグ

図 3.6.22 出力結果の比較
上：Classical PSHA（Simple fault source）の計算結果
下：Point Source を用いて断層面を一枚とした場合
左: ハザードカーブ、中: 一様ハザードスペクトル、右: ハザードマップ

131Event-based PSHA
次に、「(3) Classical PSHA（Simple fault source）」（p.128~）で使用した震源モデルによる Event-based PSHA の計算例を説明する。Event-based PSHA による解析は、Classical PSHAの解析とは異なり、最初に多数の地震動場を計算する。計算された多数の地震動場からモンテカルロ法を用いることで、確率論的な統計データとして処理することになる。

1)

必要なファイルClassical PSHA と同様に、OpenQuake-engine 実行のための設定ファイルと震源モデルロ
ジックツリーおよび地震動モデルロジックツリーが必要となる。
設定ファイル（ファイル 3.6.12）
震源モデルロジックツリーファイル（ファイル 3.6.7）
地震動モデルロジックツリーファイル（ファイル 3.6.8）
震源モデルロジックツリー内で指定される震源モデル（ファイル 3.6.10）もまた必要となる。震源モデルロジックツリーおよび地震動モデルロジックツリーは「(2)Classical PSHA

（Point source）」（p.122~）、震源モデルファイルは、「(3) Classical PSHA（Simple fault source）」（p.128~）と同じファイル（ファイル 3.6.7、ファイル 3.6.8、ファイル 3.6.10）を使用する。
今回の計算で使用する設定ファイルを（ファイル 3.6.12）に示す。
ファイル 3.6.12
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26

Event-based PSHA で使用する設定ファイル

ファイル内容[general]description = Event Based PSHA calculation with Simple fault
calculation_mode = event_based

説明
計算概要
モードの選択[Hazard site]
region = 80 26, 88 26, 88 31, 80 31
region_grid_spacing = 2

計算領域をカンマ区切りで指定(度)region 範囲内の計算地点の刻み(km)[logic_tree]source_model_logic_tree_file = source_model_logic_tree.xml
gsim_logic_tree_file = gmpe_logic_tree.xml

震源モデルロジックツリー
地震動予測式ロジックツリー[erf]width_of_mfd_bin = 0.1
rupture_mesh_spacing = 5

マグニチュード別度数分布の刻み
断層モデル作成時のメッシュ分割サイズ (km)[site_params]
reference_vs30_value = 760.0
reference_vs30_type = measured
reference_depth_to_2pt5km_per_sec = 5.0
reference_depth_to_1pt0km_per_sec = 100.0

Vs30(m/s)
Vs30 のタイプ
Vs=2.5km/s となる深さ(km)Vs=1.0km/s となる深さ(m)[correlation]
ground_motion_correlation_model =[Hazard calculation]

132

27
28

29
30
31
32
33
34
35
36
37
38
39
40

ファイル内容random_seed = 113
intensity_measure_types_and_levels = {"PGA": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(0.025)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13], "SA(0.05)": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(0.1)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13], "SA(0.2)": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(0.5)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13], "SA(1.0)": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(2.0)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13]}
truncation_level = 3.0
maximum_distance = 500
investigation_time = 50.0

説明
乱数のシード値

計算する地震動強さ
地震動強さの数値も必要となる

切断正規分布の有界レベル計算する最大距離
発生確率の対象期間[event_based_params]
ses_per_logic_tree_path=100

統計処理の回数[hazard_outputs]
ground_motion_fields = true
hazard_curves_from_gmfs = true
hazard_maps = true
poes = 0.1 0.02

GMF の出力
GMF からのハザードカーブを出力
ハザードマップの出力超過確率 10%、2%

この例では、ファイル 3.6.12 の設定ファイルとファイル 3.6.7 およびファイル 3.6.8 の各
ロジックツリーファイル、ファイル 3.6.10 の震源モデルに従い、確率論的地震ハザード解析を行う。多くの設定は Classical PSHA と同様であるが、ses_per_logic_tree_path 変数（ファイル 3.6.12 中 34 行目）を設定し、統計処理の回数を指定する。モンテカルロ法により確率を近似的に求めるため、回数が少なければ近似は荒く、回数が多ければよい近似となる。
多ければ、計算に時間がかかることは言うまでもない。

2)

計算方法

計算方法は「3.4.1

WebUI による実行」（p.81~）および「3.4.2 コマンドラインによる計算実行方法」（p.90~）に示したとおりである。

3)

計算結果

Event Based PSHA による計算結果は、設定ファイル中のパラメータの指定およびロジッ
クツリーでの経路によって複数の計算結果が出力される。今回の計算では、ロジックツリ

133

ることから、ハザードカ
カーブ、ハザザードマップ
プおよび一様様ハザードスペクーは 1 経路である
トルののみが計算可可能である。計算結果は
は出力の種類
類ごとに、csv 形式と nrm
ml 形式、geo
ojson
形式、、npz 形式のの 4 パターンンで出力可能能である（表表 3.6.4）。
表 3.6
6.4

Event--based PSHA
A の出力結果
果一覧

計算結果
果

名前

出力
力形式

地震動強
強さの分布

gmf_ddatanrml 形式
csv 形式形
npz 形式形

ハザードカーブhcurvees

ハザードマップ

hmap

一様ハザ
ザードスペクトル

uhs

断層モデ
デルデータ

rup_daata

ロジックツリーパス

realizaations

震源グル
ループ

sourceegroups

nrml 形式
csv 形式形
geojsson 形式
npz 形式形
csv 形式形

イルは汎用ソフトウェ アを使用しして図化を行行う。nrml 形
形式ファイルは
csvv 形式ファイ
OpennQuake 開発者
者が作成したたツールを使使用して図化化することがが可能である
る。図化した計算
結果のの例を図 3.66.23 に示す。

3
Eve nt-based PS
SHA の計算結結果
図 3.6.23
左: ハザードカ ーブ、右: ハザードマッ
ハップ
考のため、sses_per_logicc_tree_path 変
変数を 10、100、1000 とした場合の 計算結果と「(3)参考Classsical PSHA（SSimple fault source）
s」（p.1128~）によるる計算結果をを図 3.6.24 ににそれぞれ示す。
震源
源モデルフ ァイルや設定ファイル
ルは、ほとんんど「(3)

Cllassical PSH
HA（Simple fault

sourcce）」（p.128~~）と同様のもものを使用ししているが、ses_per_log
gic_tree_path 変数の値で結果が異ななることが確確認できる。計算に要すする時間につついては、Cllassical PSHA
A の場合は、626
秒、E
Event -based PSHA
Pの場合
合は、ses_perr_logic_tree__path の値が 10 で 32 秒、、100 で 132 秒、1000 で 312 秒となっているが、ses_per__logic_tree_ppath の値が 1000 の場合
合でも収束しているとはは言えないたため、より多多い数での計計算が必要でであると考ええる。

134

ses_per_logic_tree_path = 10

ses_per_logic_tree_path = 100

ses_per_logic_tree_path = 1000

Classical PSHA

図 3.6.24 Event-based PSHA と Classical PSHA の計算結果の比較
左: ハザードカーブ、右: ハザードマップ

135

Disaggregation
ハザード解析の最後に Disaggregation の計算を解説する。Disaggregation はハザードの再分解のことである。

1)

必要なファイル

OpenQuake-engine 実行のための設定ファイルと震源モデルロジックツリーおよび地震動
モデルロジックツリーが必要となる。
設定ファイル（ファイル 3.6.13）
震源モデルロジックツリーファイル（ファイル 3.6.7）
地震動モデルロジックツリーファイル（ファイル 3.6.14）
震源モデルロジックツリー内で指定される震源モデル（ファイル 3.6.15）もまた必要となる。地震動モデルロジックツリーは「(2) Classical PSHA（Point source）」（p.122~）と同じファイル（ファイル 3.6.7）を使用する。今回の計算で使用する設定ファイル、地震動モデ
ルロジックツリーファイルおよび震源モデルファイルをファイル 3.6.13、ファイル 3.6.14
およびファイル 3.6.15 に示す。
ファイル 3.6.13
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26

Disaggregation で使用する設定ファイル

ファイル内容[general]description = Disaggregation calculation with Source Modelcontaining an Area Source and a Simple Fault Source
belonging to different tectonic region types
calculation_mode = disaggregation

説明
計算概要
モードの選択[sites]sites = 0.5 -0.5

サイトを点で指定[logic_tree]source_model_logic_tree_file = source_model_logic_tree.xml
gsim_logic_tree_file = gmpe_logic_tree.xml
number_of_logic_tree_samples = 0

震源モデルロジックツリー
地震動予測式ロジックツリー
ロジックツリーはすべて使用[erf]
rupture_mesh_spacing = 2
width_of_mfd_bin = 0.1
area_source_discretization = 5.0

震源距離計算時に断層モデルを分割するサイズ(km)マグニチュード別度数分布の刻み
area source の離散化サイズ (km)[site_params]
reference_vs30_type = measured
reference_vs30_value = 600.0
reference_depth_to_2pt5km_per_sec = 5.0
reference_depth_to_1pt0km_per_sec = 100.0

Vs30(m/s)
Vs30 のタイプ
Vs=2.5km/s となる深さ(km)Vs=1.0km/s となる深さ(m)[hazard_calculation]
random_seed = 23
investigation_time = 50.0

乱数のシード値
発生確率の対象期間

136

27
28
29
30
31
32
33
34
35
36

ファイル内容intensity_measure_types_and_levels = {"PGA": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13]}
truncation_level = 3
maximum_distance = 200.0

説明

切断正規分布の有界レベル計算する最大距離[disaggregation]
poes_disagg = 0.1
mag_bin_width = 1.0
distance_bin_width = 10.0
coordinate_bin_width = 0.2
num_epsilon_bins = 3

超過確率の分解
マグニチュードの分解の刻み
距離の分解の刻み
座標の分解の刻み
epsilon の数計算する地震動強さ
地震動強さの数値も必要となる*斜体は選択した GSIM によって、使用の有無や変数が異なる
ファイル 3.6.14

Disaggregation で使用する地震動予測式ロジックツリーファイル

1
2

ファイル内容<?xml version="1.0" encoding="UTF-8"?><nrml xmlns="http://openquake.org/xmlns/nrml/0.5">3<logicTree logicTreeID='lt1'>

4<logicTreeBranchingLevel branchingLevelID="bl1">

5<logicTreeBranchSet uncertaintyType="gmpeModel"
branchSetID="bs1" applyToTectonicRegionType="Active
Shallow Crust">

6<logicTreeBranch branchID="b11">

7
8
9
10
11
12

13

<uncertaintyModel>ChiouYoungs2008</uncertaintyModel><uncertaintyWeight>1.0</uncertaintyWeight></logicTreeBranch></logicTreeBranchSet></logicTreeBranchingLevel><logicTreeBranchingLevel branchingLevelID="bl2"><logicTreeBranchSet uncertaintyType="gmpeModel"
branchSetID="bs2" applyToTectonicRegionType="Stable
Continental Crust">

14<logicTreeBranch branchID="b21">

15<uncertaintyModel>ToroEtAl2002</uncertaintyModel>

16
17
18
19

<uncertaintyWeight>1.0</uncertaintyWeight></logicTreeBranch></logicTreeBranchSet></logicTreeBranchingLevel>137

説明XML 宣言nrml: NRML 開始タグ
xmlns: XML 名前空間
logicTree: ロジックツリー
logicTreeID: ID"lt1"
logicTreeBranchingLevel: ブランチ
ングレベル
branchingLevelID: ID"bl1"
logicTreeBranchSet: ブランチセット
branchSetID: ID"bs1”、
uncertaintyType: gmpeModel、
applyToTectonicRegionType: 適用するテクトニックタイプは Active
Shallow Crust
logicTreeBranch: ブランチ
branchID: ID b1"
uncertaintyModel: 不確定性は
ChuouYoungs2008
uncertaintyWeight: 重みは 1.0
logicTreeBranch 終了タグ
logicTreeBranchSet 終了タグ
logicTreeBranchingLevel 終了タグ
logicTreeBranchingLevel: ブランチ
ングレベル
branchingLevelID: ID"bl1"
logicTreeBranchSet: ブランチセット
branchSetID: ID"bs2”、
uncertaintyType: gmpeModel、
applyToTectonicRegionType: 適用するテクトニックタイプは Stable
Continental Crust
logicTreeBranch: ブランチ
branchID: ID b1"
uncertaintyModel:
不確定性は ToroEtAl2002
uncertaintyWeight: 重みは 1.0
logicTreeBranch 終了タグ
logicTreeBranchSet 終了タグ
logicTreeBranchingLevel 終了タグ

20
21

ファイル内容</logicTree></nrml>ファイル 3.6.15
1
2
3
4
5

説明
logicTree 終了タグ
nrml 終了タグ

Disaggregation で使用する震源モデルファイル

ファイル内容<?xml version='1.0' encoding='utf-8'?><nrml xmlns:gml="http://www.opengis.net/gml"
xmlns="http://openquake.org/xmlns/nrml/0.5"><sourceModel name="Example Source Model Containing an
Area Source and a Simple Fault Source belonging to different
tectonic region types"><sourceGroup name="g1" tectonicRegion="Active Shallow
Crust"><areaSource id="1" name="Area Source"
tectonicRegion="Stable Continental Crust">

6<areaGeometry>

7
8
9
10
11
12
13
14
15
16
17
18
19

<gml:Polygon><gml:exterior><gml:LinearRing><gml:posList>-0.5 -0.5-0.3 -0.1
0.1 0.2
0.3 -0.8</gml:posList></gml:LinearRing></gml:exterior></gml:Polygon><upperSeismoDepth>0.0</upperSeismoDepth>

20<lowerSeismoDepth>10.0</lowerSeismoDepth>

21
22</areaGeometry><magScaleRel>WC1994</magScaleRel>23<ruptAspectRatio>1.0</ruptAspectRatio>

24<truncGutenbergRichterMFD aValue="4.5" bValue="1.0"minMag="5.0" maxMag="7.0" />25
26<nodalPlaneDist><nodalPlane probability="1.0" strike="0.0" dip="90.0"
rake="0.0" />

27
28
29

</nodalPlaneDist><hypoDepthDist><hypoDepth probability="1.0" depth="5.0" />138

説明
xml 宣言nrml: NRML 開始タグ
xmlns: XML 名前空間
sourceModel: 震源モデルファイル
name: 簡単な説明
sourceGroup: 震源グループ
tectonicRegion: テクトニックタイプ
areaSource : Area Source
name: 簡単な説明、tectonicRegion:
テクトニックタイプは Stable
Continental Crust
areaGeometry: Area ジオメトリを使用
gml:Polygon: ポリゴン
gml:exterior: 通常の領域
gml:LinearRing: ポリゴンの線
gml:posList: ポジションリスト

gml:posList 終了タグ
gml:LinearRing 終了タグ
gml:exterior 終了タグ
gml:Polygon 終了タグ
upperSeismoDepth: 地震発生上面深さ 0．0(km)lowerSeismoDepth: 地震発生下面深さ 10.0(km)areaGeometry 終了タグ
magScaleRel: スケーリング則はWC1994
ruptAspectRatio: 断層面のアスペク
ト比は 1.5
truncGutenbergRichterMFD: G-R 則に基づくマグニチュード別度数分布
aValue: 4.5、bValue: 1.0、minMag: 5.0、
maxMag: 7.0
nodalPlaneDist: 断層面のばらつき
nodalPlane: 断層面
probability: 発生確率、strike: 走向、
dip: 傾斜、rake: すべり角
傾斜 90 度、すべり角 0 度、走向 0 度の断層面の発生確率が 100%
nodalPlaneDist 終了タグ
hypoDepthDist: 震源深さのばらつき
hypoDepth: 震源深さ
probability: 発生確率、depth: 深さ、深さ 5km の発生確率が 100%

30
31
32

ファイル内容</hypoDepthDist></areaSource><simpleFaultSource id="2" name="Simple Fault Source"
tectonicRegion="Active Shallow Crust">

33
34
35
36
37
38
39
40
41
42

<simpleFaultGeometry><gml:LineString><gml:posList>
1.0 -0.5
1.4 0.0
1.4 0.3</gml:posList></gml:LineString><dip>30.0</dip><upperSeismoDepth>8.0</upperSeismoDepth>

43<lowerSeismoDepth>20.0</lowerSeismoDepth>

44
45</simpleFaultGeometry><magScaleRel>WC1994</magScaleRel>46<ruptAspectRatio>2.0</ruptAspectRatio>

47<truncGutenbergRichterMFD aValue="3.2" bValue="0.9"minMag="6.5" maxMag="7.5" />

48
49
50
51
52

<rake>90.0</rake></simpleFaultSource></sourceGroup></sourceModel></nrml>説明
hypoDepthDist 終了タグ
areaSource 終了タグ
simpleFaultSource : Simple fault
name: 簡単な説明、tectonicRegion:
テクトニックタイプは Active
Shallow CrustsimpleFaultGeometry: Simple fault
gml:Point: 線
gml:posList: 位置情報のリスト
始点の経度 1.0 度、緯度-0.5 度
始点の経度 1.4 度、緯度 0.0 度
始点の経度 1.4 度、緯度 0.3 度
gml:posList 終了タグ
gml:LineString 終了タグ
dip: 傾斜角 30 度
upperSeismoDepth: 地震発生上面深さ 8．0kmlowerSeismoDepth: 地震発生下面深さ 20.0km
simpleFaultGeometry 終了タグ
magScaleRel: スケーリング則はWC1994
ruptAspectRatio: 断層面のアスペク
ト比は 2.0
truncGutenbergRichterMFD: G-R 則に基づくマグニチュード別度数分布
aValue: 4.5、bValue: 1.0、minMag: 5.0、
maxMag: 7.0
rake: すべり角 90 度
simpleFaultSource 終了タグ
sourceGroup 終了タグ
sourceModel 終了タグ
nrml 終了タグ

この例では、ファイル 3.6.13 の設定ファイルと、ファイル 3.6.7 およびファイル 3.6.14の各ロジックツリーファイル、ファイル 3.6.15 の震源モデルに従い、ある観測点（ファイル 3.6.13 の 6 行目）のハザードの再分解を行う。

2)

計算方法

計算方法は「3.4.1

WebUI による実行」（p.81~）および「3.4.2 コマンドラインによる計算実行方法」（p.90~）に示したとおりである。

3)

計算結果

計算結果については、ここでは省略する。

139

3.6.2

リスク

Scenario Damage
まず、リスク計算で一番簡単な例として、Scenario Damage の計算例を示す。Scenario
Damage は、Scenario Case（Scenario-based Hazard）による地震動強さの分布を使用して、exposure モデルと fragility モデルを用いて、建物被害などを求めるものである。

1)

必要なファイル

OpenQuake-engine による計算では、まずハザードの計算を行い、その計算結果をもとにリスクの計算を行うことになる。ここでは、Scenario Case タイプのハザードの計算結果を用いて損害分布を求めるため、以下の計 4 ファイルが必要となる。オプションとして、
consequence モデルファイルを指定することも可能である。
設定ファイル（ファイル 3.6.16）
震源モデルファイル（ファイル 3.6.2）
exposure モデルファイル（ファイル 3.6.17）
fragility モデルファイル（ファイル 3.6.18）（オプションとして）consequence モデルファイル（ファイル 3.6.19）
設定ファイル、exposure モデルファイル、fragility モデルファイルおよび consequence モデ
ルファイルをファイル 3.6.17、ファイル 3.6.17、ファイル 3.6.18 およびファイル 3.6.19 にそれぞれ示す。震源モデルファイルは「3.6.1(1)

Scenario Case」（p.106~）で使用したファ

イ ル 3.6.2 を そ の ま ま 使 用 す る 。 ハ ザ ー ド と リ ス ク の 設 定 フ ァ イ ル は 別 フ ァ イ ル（job_hazard.ini、job_risk.ini など）とすることも可能であるが、今回の計算では同じファ
イルに両方の設定を記述した。
ファイル 3.6.16
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18

Scenario Damage で使用する設定ファイル

ファイル内容[general]description = Scenario Damage calculation with Simple fault
rupturecalculation_mode = scenario_damage

説明
計算概要
モードの説明[sites]region_grid_spacing = 5
region = 80 26, 88 26, 88 31, 80 31

計算間隔 (km)範囲（経度 緯度）

[erf]
rupture_model_file = earthquake_rupture_model.xml
rupture_mesh_spacing = 2

断層モデル
断層モデルファイル
断層分割サイズ (km)[site_params]
reference_vs30_value = 760
reference_vs30_type = measured
reference_depth_to_2pt5km_per_sec = 5.0
reference_depth_to_1pt0km_per_sec = 100.0

観測点パラメータ
Vs30 (m/s)
Vs30 のタイプ
Vs=2.5km/s となる深さ(km)Vs=1.0km/s となる深さ(m)

140

19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36

ファイル内容[correlation]
ground_motion_correlation_model =説明
観測点パラメータ
地震動補正モデル[hazard_calculation]
random_seed = 113
gsim = ChiouYoungs2008
truncation_level = 3.0
maximum_distance = 500.0
number_of_ground_motion_fields = 10

乱数のシード値
距離減衰式
切断正規分布の有界レベル
最大距離 (km)
地震動場の計算回数[exposure]
exposure_file = exposure_model.xml

exposure モデルファイル[fragility]
structural_fragility_file = structural_fragility_model.xml

fragility モデルファイル[consequence]
structural_consequence_file = consequence_model.xml

consequence モデルファイル

ファイル 3.6.17
1
2
3

Scenario Damage で使用する exposure モデルファイルの一部

ファイル内容<?xml version="1.0" encoding="UTF-8"?><nrml xmlns="http://openquake.org/xmlns/nrml/0.5"><exposureModel id="ep" category="buildings">

45<description>Exposure Model for Nepal </description><conversions>

6
7<costTypes><costType name="structural" type="per_asset"
unit="USD"/><costType name="nonstructural" type="per_asset"
unit="USD"/><costType name="contents" type="per_asset"
unit="USD"/></costTypes><insuranceLimit isAbsolute="false"/>

8
9
10
11
12

<deductible isAbsolute="false"/>

13
14
15

</conversions><assets><asset id="a1" number="1" taxonomy="Wood">

16<location lon="83.313823" lat="29.461172"/>

17
18<costs><cost type="structural" value="11340"
deductible="0.1" insuranceLimit="0.8"/><cost type="nonstructural" value="19845"
deductible="0.1" insuranceLimit="0.8"/><cost type="contents" value="5670" deductible="0.1"
insuranceLimit="0.8"/>

19
20

141

説明
XML 宣言nrml 開始タグ
exposureModel: exposure モデル
id: ID "ep"、category: カテゴリーはbuildings
description: 説明
conversions: コンバージョンセクシ
ョン
costTypes: コストタイプの集合
costType: コストタイプ
name: 分類は構造物、非構造物、コンテンツ、type: コストタイプの単位は資産ごと、unit: 単価は USD
costTypes 終了タグ
insuranceLimit: 保険限度額
isAbsolute: 絶対値か否かは false
deductible: 控除免責金額
isAbsolute: 絶対値か否かは false
conversions 終了タグ
assets: asset の集合
asset: 資産
id: ユニークな文字列、taxonomy: 分
類は木造、number: 数は 1
location: 位置情報
lon: 経度 83.313823 度、lat: 緯度
29.461172 度
costs: cost の集合
cost: コスト
type: タイプは構造物、非構造物、コンテンツ、value: 価値は USD、
deductible: 控除免責金額は 10%、insuranceLimit: 保険限度額は 80%

21
22
23

ファイル内容</costs><occupancies><occupancy occupants="5" period="night"/>

24
25
26
27
28

</occupancies></asset></assets></exposureModel></nrml>ファイル 3.6.18
1
2

Scenario Damage で使用する fragility モデルファイル

ファイル内容<?xml version="1.0" encoding="utf-8"?><nrml xmlns="http://openquake.org/xmlns/nrml/0.5">>

3<fragilityModel assetCategory="building" id="fm"lossCategory="structural">4<description>Fragility model for Nepal(discrete)</description><limitStates>slight moderate extreme
complete</limitStates><fragilityFunction format="discrete" id="Concrete"shape="logncdf">

5
6

7<imls imt="PGA" noDamageLimit="0.05 ">0.0 0.2 0.4 0.6
0.8 1.0 1.2 1.4</imls>

8<poes ls="slight">0.0 0.788 0.97 0.994 0.998 1.0 1.0
1.0</poes><poes ls="moderate">0.0 0.6 0.909 0.976 0.992 0.997
0.999 1.0</poes><poes ls="extreme">0.0 0.341 0.75 0.904 0.96 0.982 0.992
0.996</poes><poes ls="complete">0.0 0.109 0.441 0.686 0.825 0.9
0.942 0.965</poes></fragilityFunction><fragilityFunction format="discrete"
id="Unreinforced-Brick-Masonry" shape="logncdf">

9
10
11
12
13

14

<imls imt="PGA" noDamageLimit="5.0000000E-02">0.0
0.2 0.4 0.6 0.8 1.0 1.2 1.4</imls>

15<poes ls="slight">0.0 0.771 0.95 0.989 0.997 0.999 1.0
1.0</poes><poes ls="moderate">0.0 0.5 0.861 0.957 0.985 0.994
0.997 0.999</poes><poes ls="extreme">0.0 0.231 0.636 0.837 0.924 0.962
0.981 0.989</poes>

16
17

説明
costs 終了タグ
occupancies: occupancy の集合
occupancy: 居住者
occupants: 居住者 5、period: 期間は夜
occupancies 終了タグ
asset 終了タグ
assets 終了タグ
exposureModel 終了タグ
nrml 終了タグ

142

説明
XML 宣言nrml 開始タグ
xmlns: XML 名前空間
fragilityModel:
id: ID "fm"、assetCategory: 資産はbuilding、lossCategory: 損失はstructuraldescription: 簡単な説明
limitStates: 被害状態は、lightmoderate extreme complete の 4 状態
fragilityFunction: fragility 関数
id: "Concrete"、format: フォーマット
は discrete、shape: 形状は logncdf
対数正規累積分布関数による fragility
関数
imls: intensity measure levels
imt: intensity measure type は PGA、
noDamageLimit: 被害なしの制限値
poes: probabilities of exceedances
ls: limit state ごとに設定

fragilityFunction 終了タグ
fragilityFunction: fragility 関数
id: "Unreinforced-Brick-Masonry"、format: フォーマットは discrete、
shape: 形状は logncdf
対数正規累積分布関数による fragility
関数
imls: intensity measure levels
imt: intensity measure type は PGA、
noDamageLimit: 被害なしの制限値
poes: probabilities of exceedances
ls: limit state ごとに設定

18
19
20

ファイル内容<poes ls="complete">0.0 0.097 0.414 0.661 0.806 0.887
0.933 0.959</poes></fragilityFunction><fragilityFunction format="discrete" id="Wood"
shape="logncdf">

21<imls imt="PGA" noDamageLimit="5.0000000E-02">0.0
0.2 0.4 0.6 0.8 1.0 1.2 1.4</imls>

22<poes ls="slight">0.0 0.5 0.861 0.957 0.985 0.994 0.997
0.999</poes><poes ls="moderate">0.0 0.204 0.6 0.813 0.909 0.954
0.976 0.986</poes><poes ls="extreme">0.0 0.041 0.255 0.49 0.664 0.78 0.855
0.903</poes><poes ls="complete">0.0 0.007 0.088 0.236 0.394 0.532
0.642 0.728 </poes></fragilityFunction><fragilityFunction format="discrete" id="Adobe"
shape="logncdf">

23
24
25
26
27

28

<imls imt="PGA" noDamageLimit="5.0000000E-02">0.0
0.2 0.4 0.6 0.8 1.0 1.2 1.4</imls>

29<poes ls="slight">0.0 0.75 0.96 0.992 0.998 0.999 1.0
1.0</poes><poes ls="moderate">0.0 0.6 0.909 0.976 0.992 0.997
0.999 1.0</poes><poes ls="extreme">0.0 0.341 0.75 0.904 0.96 0.982 0.992
0.996</poes><poes ls="complete">0.0 0.168 0.548 0.775 0.886 0.94
0.967 0.981</poes></fragilityFunction><fragilityFunction format="discrete" id="Stone-Masonry"
shape="logncdf">

30
31
32
33
34

35

<imls imt="PGA" noDamageLimit="5.0000000E-02">0.0
0.2 0.4 0.6 0.8 1.0 1.2 1.4</imls>

36<poes ls="slight">0.0 0.673 0.937 0.985 0.996 0.998 0.999
1.0</poes><poes ls="moderate">0.0 0.532 0.878 0.964 0.988 0.995
0.998 0.999</poes><poes ls="extreme">0.0 0.3 0.711 0.883 0.95 0.977 0.989
0.994</poes><poes ls="complete">0.0 0.148 0.516 0.75 0.869 0.929
0.96 0.977</poes></fragilityFunction></fragilityModel></nrml>

37
38
39
40
41
42

143

説明

fragilityFunction 終了タグ
fragilityFunction: fragility 関数
id: "Wood"、format: フォーマットはdiscrete、shape: 形状は logncdf
対数正規累積分布関数によるフラジ
リティ関数
imls: intensity measure levels
imt: intensity measure type は PGA、
noDamageLimit: 被害なしの制限値
poes: probabilities of exceedances
ls: limit state ごとに設定

fragilityFunction 終了タグ
fragilityFunction: fragility 関数
id: "Adobe"、format: フォーマットはdiscrete、shape: 形状は logncdf
対数正規累積分布関数による fragility
関数
imls: intensity measure levels
imt: intensity measure type は PGA
noDamageLimit: 被害なしの制限値
poes: probabilities of exceedances
ls: limit state ごとに設定

fragilityFunction 終了タグ
fragilityFunction: fragility 関数
id: "Stone-Masonry"、format: フォー
マットは discrete、shape: 形状はlogncdf
対数正規累積分布関数による fragility
関数
imls: intensity measure levels
imt: intensity measure type は PGA、
noDamageLimit: 被害なしの制限値
poes: probabilities of exceedances
ls: limit state ごとに設定

fragilityFunction 終了タグ
fragilityModel 終了タグ
nrml 終了タグ

ファイル 3.6.19
1
2

Scenario Damage で使用する consequence モデルファイル

ファイル内容<?xml version="1.0" encoding="utf-8"?><nrml xmlns="http://openquake.org/xmlns/nrml/0.5">

3<consequenceModel id="cm-italy" assetCategory="buildings"
lossCategory="structural">

4
5<description>Consequence Model - Italy</description><limitStates>slight moderate extreme collapse</limitStates>

6<consequenceFunction id="Concrete" dist="LN">

7
8
9
10
11
12

<params ls="slight" mean="0.04" stddev="0.00"/><params ls="moderate" mean="0.31" stddev="0.00"/><params ls="extreme" mean="0.60" stddev="0.00"/><params ls="collapse" mean="1.00" stddev="0.00"/></consequenceFunction><consequenceFunction id="Unreinforced-Brick-Masonry"
dist="LN">

13
14
15
16
17
18

<params ls="slight" mean="0.04" stddev="0.00"/><params ls="moderate" mean="0.31" stddev="0.00"/><params ls="extreme" mean="0.60" stddev="0.00"/><params ls="collapse" mean="1.00" stddev="0.00"/></consequenceFunction><consequenceFunction id="Wood" dist="LN">

19
20
21
22
23
24

<params ls="slight" mean="0.04" stddev="0.00"/><params ls="moderate" mean="0.31" stddev="0.00"/><params ls="extreme" mean="0.60" stddev="0.00"/><params ls="collapse" mean="1.00" stddev="0.00"/></consequenceFunction><consequenceFunction id="Adobe" dist="LN">

25
26
27
28
29
30

<params ls="slight" mean="0.04" stddev="0.00"/><params ls="moderate" mean="0.31" stddev="0.00"/><params ls="extreme" mean="0.60" stddev="0.00"/><params ls="collapse" mean="1.00" stddev="0.00"/></consequenceFunction><consequenceFunction id="Stone-Masonry" dist="LN">

31
32
33
34
35
36
37

<params ls="slight" mean="0.04" stddev="0.00"/><params ls="moderate" mean="0.31" stddev="0.00"/><params ls="extreme" mean="0.60" stddev="0.00"/><params ls="collapse" mean="1.00" stddev="0.00"/></consequenceFunction></consequenceModel></nrml>144

説明XML 宣言nrml 開始タグ
xmlns: XML 名前空間
consequenceModel: consequence モデ
ル
id: "cm-italy"、assetCategory: 資産はbuildings、lossCategory: 損失はstructuraldescription: 簡単な説明
limitStates: 被害の状態は slightmoderate extreme collapse の 4 状態
consequenceFunction: 被害関数
id: "Concrete"、dist: LN は対数正規分
布による consequence 関数
params: 対数正規分布のパラメータls: limit state は slight、mean: 平均値、
stddev: 標準偏差
consequenceFunction 終了タグ
consequenceFunction: 被害関数
id: "Unreinforced-Brick-Masonry "、dist: LN は対数正規分布による
consequence 関数
params: 対数正規分布のパラメータls: limit state は slight、mean: 平均値、
stddev: 標準偏差
consequenceFunction 終了タグ
consequenceFunction: 被害関数
id: "Wood"、dist: LN は対数正規分布
による consequence 関数
params: 対数正規分布のパラメータls: limit state は slight、mean: 平均値、
stddev: 標準偏差
consequenceFunction 終了タグ
consequenceFunction: 被害関数
id: "Adobe "、dist: LN は対数正規分
布による consequence 関数
params: 対数正規分布のパラメータls: limit state は slight、mean: 平均値、
stddev: 標準偏差
consequenceFunction 終了タグ
consequenceFunction: 被害関数
id: "Adobe"、dist: LN は対数正規分布
による consequence 関数
params: 対数正規分布のパラメータls: limit state は slight、mean: 平均値、
stddev: 標準偏差
consequenceFunction 終了タグ
consequenceModel 終了タグ
nrml 終了タグ

fragility モデル、consequence モデル中には、exposure モデルに含まれる建物分類に対応する fragility 関数、consequence 関数がそれぞれ必要である（表 3.6.5 参照）。
表 3.6.5

exposure、fragility および consequence モデルに含まれる建物分類
exposure モデル（ファイル 3.6.17）

fragility モデル（ファイル 3.6.18）

consequence モデル（ファイル 3.6.19）

Adobe（アドベ）

省略

27 行目

24 行目

Concrete（コンクリート）

省略

6 行目

6 行目

Stone-Masonry（石造建築）

省略

34 行目

30 行目

Unreinforced-Brick-Masonry（無補強煉瓦造建築）

省略

13 行目

12 行目

Wood（木造）

15 行目

29 行目

18 行目

建物分類

ファイル 3.6.16 では、region 変数と region_grid_spacing 変数を用いてハザードの計算地点を指定した（ファイル 3.6.16 の 5~7 行目を参照）。しかしながら、リスクの計算を実行する場合には、exposure モデルファイルに含まれる位置のみの計算でもよい。その場合は、ファイル 3.6.16 の 5~7 行目を削除すれば、自動的に exposure_model 変数で指定した exposure
モデルファイル内の資産の位置のみを計算するようになる。

2)

Input Preparation Toolkit によるファイルの準備「1) 必要なファイル」（p.140~）で説明したファイルのうち震源ファイルと設定ファイ
ルの 2 ファイルを作成するため、Input Preparation Toolkit を使用する。Input Preparation
Toolkit の使用方法は、「3.5

Input Preparation Toolkit」で説明したとおりである。ここでは、

まず「Earthquake Rupture」タブで震源モデルファイルを作成し、「Configuration file」タブ
で設定ファイルを作成する。

①

震源ファイル

震源ファイルは Scenario Case と同じものを使用する。

②

設定ファイル

設定ファイルは以下の手順で作成する。
Input Preparation Toolkit の「Configuration file」タブを開く（図 3.6.25）。ここでは、シナリオケースの計算なので、「Earthquake Scenarios」タブを開く。
シナリオダメージの計算を行うので、Hazard、Risk にチェックを入れ、Damage を選択する（図 3.6.26）。
「Rupture model file」で参照ボタンをクリックして「3.6.1

(1)

2) ①

震源ファ

イル」で作成した earthquake_rupture_model.xml をアップロードする（図 3.6.27）。
「exposure model file」で参照ボタンをクリックして exposure_model.xml ほかをアップロードする（図 3.6.29）。

145

その他必要事項を記入する（図 3.6.30）。
ダウンロードボタンを押すと、必要なファイルが含まれた圧縮ファイルが作成される（図 3.6.31）。

図 3.6.25 Input Preparation Toolkit（リスク設定ファイル作成）その 1「Configuration file」タブを開き、さらに「Earthquake Scenarios」タブを開く

図 3.6.26 Input Preparation Toolkit（リスク設定ファイル作成）その 2

146

Hazard と Risk にチェックを入れ、Damage を選択する

図 3.6.27 Input Preparation Toolkit（リスク設定ファイル作成）その 3「Rupture model file」で参照ボタンをクリックして「3.6.1(1) 2) ①

震源ファイル」（p.109）で作成した earthquake_rupture_model.xml をアップロードする

図 3.6.28 Input Preparation Toolkit（リスク設定ファイル作成）その 4
earthquake_rupture_model.xml をリスクとから選択し、その他必要事項を記入する

147

図 3.6.29 Input Preparation Toolkit（リスク設定ファイル作成）その 5
exposure_model.xml などをアップロード

図 3.6.30 Input Preparation Toolkit（リスク設定ファイル作成）その 6
必要事項を記入する

148

図 3.6.31 Input Preparation Toolkit（リスク設定ファイル作成）その 7
ダウンロードボタンを押すと、必要なファイルが含まれた圧縮ファイルが作成できる

3)

計算方法

計算方法は「3.4.1

WebUI による実行」（p.81~）および「3.4.2 コマンドラインによる計算実行方法」（p.90~）に示したとおりである。

4)

計算結果

Scenario Damage による計算結果は、資産ごと、建物分類ごとの被害分布とトータルの被害分布および建物倒壊分布図である。consequence モデルを指定した場合にも、資産ごと、
建物分類ごとの被害分布とトータルの consequence データもまた出力される。計算結果は出力の種類ごとに、csv 形式と nrml 形式、npz 形式の 3 パターンで出力可能である（表 3.6.6）。
表 3.6.6

Scenario Damage の出力結果一覧

計算結果

名前

出力形式

地震動強さの分布

gmf_data

nrml 形式
csv 形式
npz 形式

資産ごとの被害分布図と倒壊分布図

dmg_by_asset_and_collapse_map

建物分類ごとの被害分布
図

dmg_by_taxon

トータルの被害

dmg_total

資産ごとの consequence

csq_by_asset

資産ごとの consequence

csq_by_taxon

トータルの consequence

csv_total

nrml 形式
csv 形式

csv 形式WebUI を使用したデータの出力方法は、図 3.6.17~図 3.6.19 に示したとおりである。csv
形式ファイルは汎用ソフトウェアを使用して図化を行う。nrml 形式ファイルは OpenQuake
開発者が作成したツールを使用して図化することが可能である。OpenQuake 開発者が作成したツールを用いて図化した計算結果の例を図 3.6.32 および図 3.6.33 に示した。

149

図 3.6.32

建物分類ごとの建物倒壊分布図および統計

図 3.6.33 建物分類ごとの損害率

150

Scenario Risk
次に、地震によるリスク計算の簡単な例として、Scenario Risk の計算例を示す。Scenario
Damage との違いは、Scenario Damage が建物等の被害数・損害数を fragility モデルから算出するのに対して、Scenario Risk は vulnerability モデルを用いて経済被害あるいは人的被害を算出することである。

1)

必要なファイル

Scenario Damage の計算と同様に、まずハザードの計算を行い、その計算結果をもとに地震によるリスクの計算を行うことになる。ここでは、Scenario Case タイプのハザードの結果を用いて、損失を求めるため、以下の 4 ファイルが必要となる。vulnerability はコンテン
ツごとに用意する。
設定ファイル（ファイル 3.6.20）
震源モデルファイル（ファイル 3.6.2）
exposure モデルファイル（ファイル 3.6.17）
vulnerability モデルファイル（ファイル 3.6.21、ファイル 3.6.22、ファイル 3.6.23）
vulnerability モデルファイルは、exposure モデルで指定したコストタイプに対応する
vulnerability モデルが必要であり、コストタイプの分だけ、vulnerability モデルファイルが必要となる。ここでは、構造物（structural）、非構造物（nonstructural）、居住者（occupants）の各 vulnerability モデルファイルを用いて計算する。今回のサンプルでは、3 種のvulnerability モデルを用意したが、必要でなければ、いずれか一つの vulnerability モデルによる計算でもかまわない。
ファイル 3.6.20
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19

Scenario Risk で使用する設定ファイル

ファイル内容[general]description = Scenario Risk calculation with Simple fault
rupturecalculation_mode = scenario_risk

説明
計算概要
モードの説明[sites]region_grid_spacing = 5
region = 80 26, 88 26, 88 31, 80 31

計算間隔 (km)範囲（経度 緯度）

[erf]
rupture_model_file = earthquake_rupture_model.xml
rupture_mesh_spacing = 2

断層モデル
断層モデルファイル
断層分割サイズ (km)[site_params]
reference_vs30_value = 760
reference_vs30_type = measured
reference_depth_to_2pt5km_per_sec = 5.0
reference_depth_to_1pt0km_per_sec = 100.0

観測点パラメータ
Vs30(m/s)
Vs30 のタイプ
Vs=2.5km/s となる深さ(km)Vs=1.0km/s となる深さ(m)[correlation]

観測点パラメータ

151

20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38

ファイル内容ground_motion_correlation_model =説明
地震動補正モデル[hazard_calculation]
random_seed = 113
gsim = ChiouYoungs2008
truncation_level = 3.0
maximum_distance = 500.0
number_of_ground_motion_fields = 10

乱数のシード値
距離減衰式
切断正規分布の有界レベル
最大距離(km)
地震動場の計算回数[exposure]
exposure_file = exposure_model.xml

exposure モデルファイル[vulnerability]
structural_vulnerability_file =structural_vulnerability_model.xml
nonstructural_vulnerability_file =nonstructural_vulnerability_model.xml
occupants_vulnerability_file =occupants_vulnerability_model.xml

vulnerability
構造物の vulnerability モデルファイ
ル
非構造物の vulnerability モデルファ
イル
居住者の vulnerability モデルファイ
ル[risk_outputs]
insured_losses = True

保険損失額の出力

ファイル 3.6.21
1
2
3

4

5
6
7
8
9

10
11
12
13

Scenario Risk で使用する構造物の vulnerability モデルファイル

ファイル内容<?xml version="1.0" encoding="utf-8"?><nrml xmlns="http://openquake.org/xmlns/nrml/0.5"><vulnerabilityModel id="Nepal" assetCategory="buildings"
lossCategory="structural"><vulnerabilityFunction dist="LN" id="Wood"><imls imt="PGA">0.0001 0.2 0.4 0.6 0.8 1.0 1.2
1.4</imls><meanLRs>0.0001 0.05 0.21 0.4 0.56 0.67 0.76
0.82</meanLRs><covLRs>0.3 0.3 0.3 0.1 0.1 0.1 0.1 0.1 </covLRs></vulnerabilityFunction><vulnerabilityFunction dist="LN" id="Stone-Masonry"><imls imt="PGA">0.0001 0.2 0.4 0.6 0.8 1.0 1.2
1.4</imls><meanLRs>0.0001 0.25 0.63 0.82 0.91 0.95 0.98
0.99</meanLRs><covLRs>0.3 0.3 0.1 0.1 0.1 0.1 0.1 0.1</covLRs></vulnerabilityFunction>

152

説明XML 宣言nrml 開始タグ
xmlns: XML 名前空間
vulnerabilityModel: vulnerability モデ
ル
id: "Nepal"、assetCategory: 資産はbuildings、lossCategory: 損失はstructuralvulnerabilityFunction: vulnerability 関
数
id: "Wood"、dist: LN は対数正規分布関数による vulnerability 関数
imls: intensity measure levels
imt: intensity measure type は PGA
meanLRs: 平均損害率
covLRs: 損害率の変動係数
vulnerabilityFunction 終了タグ
vulnerabilityFunction: vulnerability 関
数
id: " Stone-Masonry"、dist: LN は対数
正規分布関数による vulnerability 関
数
imls: intensity measure levels
imt: intensity measure type は PGA
meanLRs: 平均損害率
covLRs: 損害率の変動係数
vulnerabilityFunction 終了タグ

14

15
16
17
18
19

20
21
22
23
24

25
26
27
28
29
30

ファイル内容<vulnerabilityFunction dist="LN" id="Adobe"><imls imt="PGA">0.0001 0.2 0.4 0.6 0.8 1.0 1.2
1.4</imls><meanLRs>0.0001 0.28 0.67 0.85 0.93 0.96 0.98
0.99</meanLRs><covLRs>0.3 0.3 0.1 0.1 0.1 0.1 0.1 0.1</covLRs></vulnerabilityFunction><vulnerabilityFunction dist="LN" id="Concrete"><imls imt="PGA">0.0001 0.2 0.4 0.6 0.8 1.0 1.2
1.4</imls><meanLRs>0.0001 0.25 0.61 0.8 0.9 0.94 0.97
0.98</meanLRs><covLRs>0.3 0.3 0.1 0.1 0.1 0.1 0.1 0.1</covLRs></vulnerabilityFunction><vulnerabilityFunction dist="LN"
id="Unreinforced-Brick-Masonry"><imls imt="PGA">0.0001 0.2 0.4 0.6 0.8 1.0 1.2
1.4</imls><meanLRs>0.0001 0.2 0.55 0.76 0.87 0.93 0.96
0.98</meanLRs><covLRs>0.3 0.3 0.1 0.2 0.1 0.1 0.1 0.1</covLRs></vulnerabilityFunction></vulnerabilityModel></nrml>ファイル 3.6.22
1
2
3<vulnerabilityModel id="Nepal" assetCategory="buildings"
lossCategory="nonstructural"><vulnerabilityFunction dist="LN" id="Wood">5<imls imt="PGA">0.0001 0.2 0.4 0.6 0.8 1.0 1.2
1.4</imls><meanLRs>0.0001 0.08 0.23 0.38 0.51 0.62 0.7
0.76</meanLRs><covLRs>0.3 0.3 0.3 0.3 0.1 0.1 0.1 0.1</covLRs></vulnerabilityFunction>

7
8

covLRs: 損害率の変動係数
vulnerabilityFunction 終了タグ
vulnerabilityFunction: バルネラビリ
ティ関数
id: " Concrete"、dist: LN は対数正規分布関数による vulnerability 関数
imls: intensity measure levels
imt: intensity measure type は PGA
meanLRs: 平均損害率
covLRs: 損害率の変動係数
vulnerabilityFunction 終了タグ
vulnerabilityFunction: vulnerability 関
数
id: " Unreinforced-Brick-Masonry"、
dist: LN は対数正規分布関数による
vulnerability 関数
imls: intensity measure levels
imt: intensity measure type は PGA
meanLRs: 平均損害率
covLRs: 損害率の変動係数
vulnerabilityFunction 終了タグ
vulnerabilityModel 終了タグ
nrml 終了タグ

Scenario Risk で使用する非構造物の vulnerability モデルファイル

ファイル内容<?xml version="1.0" encoding="utf-8"?><nrml xmlns="http://openquake.org/xmlns/nrml/0.5">

4

6

説明
vulnerabilityFunction: vulnerability 関
数
id: " Adobe"、dist: LN は対数正規分
布関数による vulnerability 関数
imls: intensity measure levels
imt: intensity measure type は PGA
meanLRs: 平均損害率153

説明XML 宣言nrml 開始タグ
xmlns: XML 名前空間
vulnerabilityModel: vulnerability モデ
ル
id: "Nepal"、assetCategory: 資産はbuildings、lossCategory: 損失はnonstructuralvulnerabilityFunction: vulnerability 関
数
id: "Wood"、dist: LN は対数正規分布関数による vulnerability 関数
imls: intensity measure levels
imt: intensity measure type は PGA
meanLRs: 平均損害率
covLRs: 損害率の変動係数
vulnerabilityFunction 終了タグ

9

10
11
12
13
14

15
16
17
18
19

20
21
22
23
24

25
26
27
28
29
30

ファイル内容<vulnerabilityFunction dist="LN" id="Stone-Masonry"><imls imt="PGA">0.0001 0.2 0.4 0.6 0.8 1.0 1.2
1.4</imls><meanLRs>0.0001 0.07 0.22 0.38 0.51 0.62 0.7
0.76</meanLRs><covLRs>0.3 0.3 0.3 0.3 0.1 0.1 0.1 0.1</covLRs></vulnerabilityFunction><vulnerabilityFunction dist="LN" id="Adobe"><imls imt="PGA">0.0001 0.2 0.4 0.6 0.8 1.0 1.2
1.4</imls><meanLRs>0.0001 0.07 0.22 0.38 0.51 0.62 0.7
0.76</meanLRs><covLRs>0.3 0.3 0.3 0.3 0.1 0.1 0.1 0.1</covLRs></vulnerabilityFunction><vulnerabilityFunction dist="LN" id="Concrete"><imls imt="PGA">0.0001 0.2 0.4 0.6 0.8 1.0 1.2
1.4</imls><meanLRs>0.0001 0.08 0.23 0.38 0.51 0.62 0.7
0.76</meanLRs><covLRs>0.3 0.3 0.3 0.3 0.1 0.1 0.1 0.1</covLRs></vulnerabilityFunction><vulnerabilityFunction dist="LN"
id="Unreinforced-Brick-Masonry"><imls imt="PGA">0.0001 0.2 0.4 0.6 0.8 1.0 1.2
1.4</imls><meanLRs>0.0001 0.07 0.22 0.38 0.51 0.62 0.7
0.76</meanLRs><covLRs>0.3 0.3 0.3 0.3 0.1 0.1 0.1 0.1</covLRs></vulnerabilityFunction></vulnerabilityModel></nrml>ファイル 3.6.23
1
2
3

説明
vulnerabilityFunction: vulnerability 関
数
id: " Stone-Masonry"、dist: LN は対数
正規分布関数による vulnerability 関
数
imls: intensity measure levels
imt: intensity measure type は PGA
meanLRs: 平均損害率
covLRs: 損害率の変動係数
vulnerabilityFunction 終了タグ
vulnerabilityFunction: vulnerability 関
数
id: "Adobe"、dist: LN は対数正規分布関数による vulnerability 関数
imls: intensity measure levels
imt: intensity measure type は PGA
meanLRs: 平均損害率
covLRs: 損害率の変動係数
vulnerabilityFunction 終了タグ
vulnerabilityFunction: vulnerability 関
数
id: "Concrete"、dist: LN は対数正規分
布関数による vulnerability 関数
imls: intensity measure levels
imt: intensity measure type は PGA
meanLRs: 平均損害率
covLRs: 損害率の変動係数
vulnerabilityFunction 終了タグ
vulnerabilityFunction: Vulnerability 関
数
id: "Unreinforced-Brick-Masonry"、
dist: LN は対数正規分布関数による
vulnerability 関数
imls: intensity measure levels
imt: intensity measure type は PGA
meanLRs: 平均損害率
covLRs: 損害率の変動係数
vulnerabilityFunction 終了タグ
vulnerabilityModel 終了タグ
nrml 終了タグ

Scenario Risk で使用する居住者の vulnerability モデルファイル

ファイル内容<?xml version="1.0" encoding="utf-8"?><nrml xmlns="http://openquake.org/xmlns/nrml/0.5"><vulnerabilityModel id="Nepal" assetCategory="buildings"
lossCategory="occupants">154

説明XML 宣言nrml 開始タグ
xmlns: XML 名前空間
vulnerabilityModel: vulnerability モデ
ル
id: "Nepal"、assetCategory: 資産はbuildings、lossCategory: 損失はnonstructural

4

5
6
7
8
9

10
11
12
13
14

15
16
17
18
19

20
21
22
23
24

25
26
27
28
29
30

ファイル内容<vulnerabilityFunction dist="LN" id="Wood"><imls imt="PGA">0.0001 0.2 0.4 0.6 0.8 1.0 1.2
1.4</imls><meanLRs>0.0001 0.001 0.009 0.024 0.039 0.053 0.064
0.073 </meanLRs><covLRs>0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 </covLRs></vulnerabilityFunction><vulnerabilityFunction dist="LN" id="Stone-Masonry"><imls imt="PGA">0.0001 0.2 0.4 0.6 0.8 1.0 1.2
1.4</imls><meanLRs>0.0001 0.022 0.077 0.112 0.13 0.139 0.144
0.147 </meanLRs><covLRs>0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0</covLRs></vulnerabilityFunction><vulnerabilityFunction dist="LN" id="Adobe"><imls imt="PGA">0.0001 0.2 0.4 0.6 0.8 1.0 1.2
1.4</imls><meanLRs>0.0001 0.025 0.082 0.116 0.133 0.141 0.145
0.147</meanLRs><covLRs>0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0</covLRs></vulnerabilityFunction><vulnerabilityFunction dist="LN" id="Concrete"><imls imt="PGA">0.0001 0.2 0.4 0.6 0.8 1.0 1.2
1.4</imls><meanLRs>0.0001 0.022 0.088 0.137 0.165 0.18 0.188
0.193</meanLRs><covLRs>0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0</covLRs></vulnerabilityFunction><vulnerabilityFunction
dist="LN"
id="Unreinforced-Brick-Masonry"><imls imt="PGA"> 0.0001 0.2 0.4 0.6 0.8 1.0 1.2
1.4</imls><meanLRs>0.0001 0.014 0.062 0.099 0.121 0.133 0.14
0.144</meanLRs><covLRs>0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0</covLRs></vulnerabilityFunction></vulnerabilityModel></nrml>155

説明vulnerabilityFunction: vulnerability 関
数
id: "Wood"、dist: LN は対数正規分布関数による vulnerability 関数
imls: intensity measure levels
imt: intensity measure type は PGA
meanLRs: 平均損害率
covLRs: 損害率の変動係数
vulnerabilityFunction 終了タグ
vulnerabilityFunction: vulnerability 関
数
id: " Stone-Masonry"、dist: LN は対数
正規分布関数による vulnerability 関
数
imls: intensity measure levels
imt: intensity measure type は PGA
meanLRs: 平均損害率
covLRs: 損害率の変動係数
vulnerabilityFunction 終了タグ
vulnerabilityFunction: vulnerability 関
数
id: " Adobe"、dist: LN は対数正規分
布関数による vulnerability 関数
imls: intensity measure levels
imt: intensity measure type は PGA
meanLRs: 平均損害率
covLRs: 損害率の変動係数
vulnerabilityFunction 終了タグ
vulnerabilityFunction: vulnerability 関
数
id: "Concrete"、dist: LN は対数正規分布関数による vulnerability 関数
imls: intensity measure levels
imt: intensity measure type は PGA
meanLRs: 平均損害率
covLRs: 損害率の変動係数
vulnerabilityFunction 終了タグ
vulnerabilityFunction: vulnerability 関
数
id: "Unreinforced-Brick-Masonry"、
dist: LN は対数正規分布関数による
vulnerability 関数
imls: intensity measure levels
imt: intensity measure type は PGA
meanLRs: 平均損害率
covLRs: 損害率の変動係数
vulnerabilityFunction 終了タグ
vulnerabilityModel 終了タグ
nrml 終了タグ

各 vulnerability モデルファイル中には、exposure モデルで指定した建物分類に対応する
vulnerability 関数が必要である（表 3.6.7 参照）。
表 3.6.7

exposure モデル、各 vulnerability モデルに含まれる建物分類
structural vulnerability モデル（ファイル 3.6.21）
exposure モデル
nonstructural vulnerability モデル（ファイル 3.6.22）（ファイル 3.6.17）
occupants vulnerability モデル（ファイル 3.6.23）

建物分類
Adobe（アドベ）

省略

14 行目

Concrete（コンクリート）

省略

19 行目

Stone-Masonry（石造建築）

省略

9 行目

Unreinforced-Brick-Masonry（無補強煉瓦造建築）

省略

24 行目

Wood（木造）

15 行目

4 行目

2)

計算方法

計算方法は「3.4.1

WebUI による実行」（p.81~）および「3.4.2 コマンドラインによる計算実行方法」（p.90~）に示したとおりである。

3)

計算結果

Scenario Risk による計算結果は、建物分類ごとの経済被害あるいは人的被害の分布が出力される。計算結果は csv 形式と nrml 形式、geojson 形式の 3 パターンで出力可能である（表 3.6.8）。
表 3.6.8

Scenario Risk の出力結果一覧

計算結果

名前

出力形式

地震動強さの分布

gmf_data

nrml 形式
csv 形式
npz 形式

総合損失額

agglosses-rlzs

csv 形式

資産ごとの損失データ

losses_by_asset

nrml 形式
csv 形式
geojson 形式

WebUI を使用したデータの出力方法は、図 3.6.17~図 3.6.19 に示したとおりである。csv
形式ファイルは汎用ソフトウェアを使用して図化を行う。nrml 形式ファイルは OpenQuake
開発者が作成したツールを使用して図化することが可能である。OpenQuake 開発者が作成したツールを用いて図化した計算結果の例を図 3.6.34 に示した。

156

図 3.6.34 Scenario Risk の計算による建物分類ごとのロスマップ

157

Classical Damage
次に、Classical Damage の計算例を示す。Scenario Damage との違いは、Scenario Damage
が Scenario Case タイプのハザードの結果を用いて損害分布を求めるのに対し、ClassicalDamage は、Classical PSHA の解析結果から、確率論的な建物被害を算出することである。

1)

必要なファイル

Scenario Damage の計算と同様に、まずハザードの計算を行い、その計算結果をもとにリスクの計算を行うことになる。Classical Damage では、Classical PSHA タイプのハザードの結果を用いて損害分布を求めるため、以下の計 6 ファイルが必要となる。オプションとして、consequence モデルファイルを指定することも可能である。
設定ファイル（ファイル 3.6.24）
震源モデルロジックツリーファイル（ファイル 3.6.7）
震源モデルロジックツリー内で指定される震源モデル（ファイル 3.6.10）
地震動モデルロジックツリーファイル（ファイル 3.6.8）
exposure モデルファイル（ファイル 3.6.17）
fragility モデルファイル（ファイル 3.6.18）（オプションとして）consequence モデルファイル（ファイル 3.6.19）
設定ファイルをファイル 3.6.24 に示す。設定ファイル以外は既出のもの（ファイル 3.6.7、
ファイル 3.6.8、ファイル 3.6.10、ファイル 3.6.17、ファイル 3.6.18 およびファイル 3.6.19）を使用する。
ファイル 3.6.24
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21

Classical Damage で使用する設定ファイル

ファイル内容[general]description = Classical Damage calculation with Simple fault
calculation_mode = classical_damage[logic_tree]source_model_logic_tree_file = source_model_logic_tree.xml
gsim_logic_tree_file = gmpe_logic_tree.xml[erf]
width_of_mfd_bin = 0.1
rupture_mesh_spacing = 2

説明
計算概要
モードの選択震源モデルロジックツリー
地震動予測式ロジックツリー

マグニチュード別度数分布の刻み
震源距離計算時に断層モデルを分割するサイズ(km)[site_params]
reference_vs30_value = 760.0
reference_vs30_type = measured
reference_depth_to_2pt5km_per_sec = 5.0
reference_depth_to_1pt0km_per_sec = 100.0

観測点パラメータ
Vs30 (m/s)
Vs30 のタイプ
Vs=2.5km/s となる深さ(km)Vs=1.0km/s となる深さ(m)[correlation]
ground_motion_correlation_model =地震動補正モデル

158

22
23
24

25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43

2)

ファイル内容[hazard_calculation]
random_seed = 113
intensity_measure_types_and_levels = {"PGA": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(0.025)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13], "SA(0.05)": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(0.1)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13], "SA(0.2)": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(0.5)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13], "SA(1.0)": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(2.0)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13]}
truncation_level = 3.0
maximum_distance = 200
investigation_time = 50.0

説明
乱数のシード値

計算する地震動指標とその強さ
複数選択時には、{ }、[ ]、” ”、およびカンマを使用して適切に記述する必要がある

切断正規分布の有界レベル計算する最大距離(km)発生確率の対象期間(年)[risk_calculation]
risk_investigation_time = 1
steps_per_interval = 4

リスクの超過確率の対象期間(年)
fragility 関数の補間間隔[exposure]
exposure_file = exposure_model.xml

exposure モデルファイル[fragility]
structural_fragility_file = structural_fragility_model.xml
structural_consequence_file = consequence_model.xml

構造物の fragility モデルファイル
構造物の consequence モデルファイル[hazard_outputs]
hazard_maps = true
uniform_hazard_spectra = true
poes = 0.1 0.02

ハザードマップを出力する
一様ハザードマップを出力する
超過確率 10%、2%

計算方法

計算方法は「3.4.1

WebUI による実行」（p.81~）および「3.4.2 コマンドラインによる計算実行方法」（p.90~）に示したとおりである。

3)

計算結果

Classical damage による計算結果は、資産ごとの確率論的な被害データが出力される。計
算結果は csv 形式と nrml 形式、geojson 形式および npz 形式の 3 パターンで出力可能である（表 3.6.9）。

159

表 3.6.9Classical Damage の出力結果一覧

計算結果

名前

出力形式

ハザードカーブhcurves

ハザードマップ

hmap

一様ハザードスペクトルuhs

nrml 形式
csv 形式
geojson 形式
npz 形式

ロジックツリーパス

realizations

震源グループsourcegroups

資産ごとの確率論的な被害データ

damages-rlzs

csv 形式
csv 形式160

Classical Risk
1)

必要なファイル

Scenario Risk の計算と同様に、まずハザードの計算を行い、その計算結果をもとにリスクの計算を行うことになる。ここでは、Classical PSHA タイプのハザードの結果を用いて、損害分布を求めるため、以下の計 6 ファイルが必要となる。
設定ファイル（ファイル 3.6.25）
震源モデルロジックツリーファイル（ファイル 3.6.7）
震源モデルロジックツリー内で指定される震源モデル（ファイル 3.6.10）
地震動モデルロジックツリーファイル（ファイル 3.6.8）
Exposure モデルファイル（ファイル 3.6.17）
Vulnerability モデルファイル（ファイル 3.6.21、ファイル 3.6.22、ファイル 3.6.23）
設定ファイルをファイル 3.6.25 に示す。設定ファイル以外は既出のもの（ファイル 3.6.7、
ファイル 3.6.8、ファイル 3.6.10、ファイル 3.6.21～ファイル 3.6.23）を使用する。
ファイル 3.6.25
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23

Classical Risk で使用する設定ファイル

ファイル内容[general]description = Classical Risk calculation with Simple fault
calculation_mode = classical_risk[logic_tree]source_model_logic_tree_file = source_model_logic_tree.xml
gsim_logic_tree_file = gmpe_logic_tree.xml[erf]
width_of_mfd_bin = 0.1
rupture_mesh_spacing = 2

説明
計算概要
モードの選択震源モデルロジックツリー
地震動予測式ロジックツリー

マグニチュード別度数分布の刻み
震源距離計算時に断層モデルを分割するサイズ(km)[site_params]
reference_vs30_value = 760.0
reference_vs30_type = measured
reference_depth_to_2pt5km_per_sec = 5.0
reference_depth_to_1pt0km_per_sec = 100.0

観測点パラメータ
Vs30(m/s)
Vs30 のタイプ
Vs=2.5km/s となる深さ(km)Vs=1.0km/s となる深さ(m)[correlation]
ground_motion_correlation_model =地震動補正モデル[hazard_calculation]
random_seed = 113

乱数のシード値

161

24

25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43

2)

ファイル内容intensity_measure_types_and_levels = {"PGA": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(0.025)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13], "SA(0.05)": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(0.1)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13], "SA(0.2)": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(0.5)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13], "SA(1.0)": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(2.0)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13]}
truncation_level = 3.0
maximum_distance = 200
investigation_time = 50.0

切断正規分布の有界レベル計算する最大距離
発生確率の対象期間(年)[risk_calculation]
lrem_steps_per_interval= 1
conditional_loss_poes= 0.10, 0.02

vulnerability 関数の補間間隔
条件付き超過確率 10%、2%[exposure]
exposure_file = exposure_model.xml

exposure モデルファイル[Vulnerability model]
structural_vulnerability_file
structural_vulnerability_model.xml
nonstructural_vulnerability_file
nonstructural_vulnerability_model.xml
occupants_vulnerability_file
occupants_vulnerability_model.xml===[hazard_outputs]
hazard_maps = true
uniform_hazard_spectra = true
poes = 0.1 0.02

説明

計算する地震動指標とその強さ
複数選択時には、{ }、[ ]、” ”、およびカンマを使用して適切に記述する必要がある構造物の Vulnerability モデルファイ
ル
非構造物の Vulnerability モデルファ
イル
居住者の Vulnerability モデルファイ
ル

ハザードマップを出力する
一様ハザードマップを出力する
超過確率 10%、2%

計算方法

計算方法は「3.4.1

WebUI による実行」（p.81~）および「3.4.2 コマンドラインによる計算実行方法」（p.90~）に示したとおりである。

3)

計算結果

Classical Risk による計算結果は、算出された建物分類ごとの確率論的な経済被害あるいは人的被害の分布および資産ごとのロスカーブが出力される。計算結果は nrml 形式、csv
形式、geojson 形式、npz 形式の 4 パターンで出力可能である（表 3.6.10）。WebUI を使用し162たデータの出力方法は、図 3.6.17~図 3.6.19 に示したとおりである。
表 3.6.10 Classical Risk の出力結果一覧
計算結果

名前

出力形式

ハザードカーブhcurves

ハザードマップ

hmap

一様ハザードスペクトルuhs

nrml 形式
csv 形式
geojson 形式
npz 形式

ロジックツリーパス

realizations

震源グループsourcegroups

ロスカーブ

ロスマップ

csv 形式

loss_curves

nrml 形式
csv 形式
geojson 形式
npz 形式

loss_maps

nrml 形式
csv 形式
geojson 形式
npz 形式

csv 形式ファイルは汎用ソフトウェアを使用して図化を行う。nrml 形式ファイルはOpenQuake 開発者が作成したツールを使用して図化することが可能である。OpenQuake 開
発者が作成したツールを用いて図化した計算結果の例を図 3.6.35 および図 3.6.36 に示した。

163

図 3.6.35 建物分類ごとのロスマップ

図 3.6.36 資産ごとのロスカーブ

164Event-based Risk
1)

必要なファイルClassical Risk の計算と同様に、まずハザードの計算を行い、その計算結果をもとにリスクの計算を行うことになる。ここでは、Event-based PSHA タイプのハザードの結果を用いて、損害分布を求めるため、以下の計 6 ファイルが必要となる。
設定ファイル（ファイル 3.6.26）
震源モデルロジックツリーファイル（ファイル 3.6.7）
震源モデルロジックツリー内で指定される震源モデル（ファイル 3.6.10）
地震動モデルロジックツリーファイル（ファイル 3.6.8）
exposure モデルファイル（ファイル 3.6.17）
Vulnerability モデルファイル（ファイル 3.6.21、ファイル 3.6.22、ファイル 3.6.23）
設定ファイルをファイル 3.6.26 に示した。設定ファイル以外は既出のもの（ファイル 3.6.7、
ファイル 3.6.8、ファイル 3.6.10、ファイル 3.6.21～ファイル 3.6.23）を使用する。
ファイル 3.6.26
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24

Event-based Risk で使用する設定ファイル

ファイル内容[general]description =Event-based Risk calculation with Simple fault
calculation_mode = event_based_risk[logic_tree]source_model_logic_tree_file = source_model_logic_tree.xml
gsim_logic_tree_file = gmpe_logic_tree.xml
number_of_logic_tree_samples = 0[erf]
width_of_mfd_bin = 0.1
rupture_mesh_spacing = 2

説明
計算概要
モードの選択震源モデルロジックツリー
地震動予測式ロジックツリー

マグニチュード別度数分布の刻み
震源距離計算時に断層モデルを分割するサイズ(km)[site_params]
reference_vs30_value = 760.0
reference_vs30_type = measured
reference_depth_to_2pt5km_per_sec = 5.0
reference_depth_to_1pt0km_per_sec = 100.0

観測点パラメータ
Vs30(m/s)
Vs30 のタイプ
Vs=2.5km/s となる深さ(km)Vs=1.0km/s となる深さ(m)[correlation]
ground_motion_correlation_model =地震動補正モデル[hazard_calculation]
random_seed = 113

乱数のシード値

165

25

26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49

2)

ファイル内容intensity_measure_types_and_levels = {"PGA": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(0.025)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13], "SA(0.05)": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(0.1)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13], "SA(0.2)": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(0.5)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13], "SA(1.0)": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(2.0)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13]}
truncation_level = 3.0
maximum_distance = 200
investigation_time = 50.0

切断正規分布の有界レベル
最大距離
発生確率の対象期間(年)[event_based_params]
ses_per_logic_tree_path=100

統計処理の回数[risk_calculation]
master_seed = 42
lrem_steps_per_interval= 1
conditional_loss_poes= 0.10, 0.02

リスク計算のための乱数のシード値
vulnerability 関数の補間間隔
条件付き超過確率 10%、2%[exposure]
exposure_file = exposure_model.xml

exposure モデルファイル[Vulnerability model]
structural_vulnerability_file
structural_vulnerability_model.xml
nonstructural_vulnerability_file
nonstructural_vulnerability_model.xml
occupants_vulnerability_file
occupants_vulnerability_model.xml===[hazard_outputs]
hazard_maps = true
uniform_hazard_spectra = true
poes = 0.1 0.02

計算する地震動指標とその強さ
複数選択時には、{ }、[ ]、” ”、およびカンマを使用して適切に記述する必要がある構造物の Vulnerability モデルファイ
ル
非構造物の Vulnerability モデルファ
イル
居住者の Vulnerability モデルファイ
ル

ハザードマップを出力する
一様ハザードマップを出力する
超過確率 10%、2%

計算方法

計算方法は「3.4.1

WebUI による実行」（p.81~）および「3.4.2 コマンドラインによる計算実行方法」（p.90~）に示したとおりである。

3)

説明

計算結果

計算結果については省略する。
166

Classical BCR
1)

必要なファイルClassical Risk の計算と同様に、まずハザードの計算を行い、その計算結果をもとに費用便益分析の計算を行うことになる。ここでは、Classical PSHA タイプのハザードの結果を用いて、地震による損失および改修費を求めるため、以下の計 7 ファイルが必要となる。
設定ファイル（ファイル 3.6.27）
震源モデルロジックツリーファイル（ファイル 3.6.7）
震源モデルロジックツリー内で指定される震源モデル（ファイル 3.6.10）
地震動モデルロジックツリーファイル（ファイル 3.6.8）
exposure モデルファイル（ファイル 3.6.17）
vulnerability モデルファイル（ファイル 3.6.21）
改修費の vulnerability モデルファイル
設定ファイルをファイル 3.6.27 に示す。設定ファイル以外は既出のもの（ファイル 3.6.7、
ファイル 3.6.8、ファイル 3.6.10、ファイル 3.6.21～ファイル 3.6.23）を使用する。改修費のvulnerability モデルファイルと vulnerability モデルファイルの違いは、指定する数値が異なるのみであるので、ここでの説明は省略する。
ファイル 3.6.27
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23

Classical BCR で使用する設定ファイル

ファイル内容[general]description = Classical BCR calculation with Simple fault
calculation_mode = classical_bcr[logic_tree]source_model_logic_tree_file = source_model_logic_tree.xml
gsim_logic_tree_file = gmpe_logic_tree.xml[erf]
width_of_mfd_bin = 0.1
rupture_mesh_spacing = 2

説明
計算概要
モードの選択震源モデルロジックツリー
地震動予測式ロジックツリー

マグニチュード別度数分布の刻み
震源距離計算時に断層モデルを分割するサイズ(km)[site_params]
reference_vs30_value = 760.0
reference_vs30_type = measured
reference_depth_to_2pt5km_per_sec = 5.0
reference_depth_to_1pt0km_per_sec = 100.0

観測点パラメータ
Vs30(m/s)
Vs30 のタイプ
Vs=2.5km/s となる深さ(km)Vs=1.0km/s となる深さ(m)[correlation]
ground_motion_correlation_model =地震動補正モデル[hazard_calculation]
random_seed = 113

乱数のシード値

167

24

25
26
27
28
29
30

31
32
33
34
35
36
37
38
39
40
41
42
43
44

2)

ファイル内容intensity_measure_types_and_levels = {"PGA": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(0.025)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13], "SA(0.05)": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(0.1)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13], "SA(0.2)": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(0.5)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13], "SA(1.0)": [0.005, 0.007,
0.0098, 0.0137, 0.0192, 0.0269, 0.0376, 0.0527, 0.0738, 0.103,
0.145, 0.203, 0.284, 0.397, 0.556, 0.778, 1.09, 1.52, 2.13],"SA(2.0)": [0.005, 0.007, 0.0098, 0.0137, 0.0192, 0.0269,
0.0376, 0.0527, 0.0738, 0.103, 0.145, 0.203, 0.284, 0.397,
0.556, 0.778, 1.09, 1.52, 2.13]}
truncation_level = 3.0
maximum_distance = 200
investigation_time = 50.0
[risk_calculation]interest_rate = 0.05

計算する地震動指標とその強さ
複数選択時には、{ }、[ ]、” ”、およびカンマを使用して適切に記述する必要がある

切断正規分布の有界レベル
最大距離
発生確率の対象期間(年)将来キャッシュフローを割り引くことで潜在的な将来利益の現在価値の計算に使われる資産の耐用年数
vulnerability 関数の補間間隔

asset_life_expectancy = 50
lrem_steps_per_interval = 1[exposure]
exposure_file = exposure_model.xml[Vulnerability model]
structural_vulnerability_file
structural_vulnerability_model.xml
structural_vulnerability_retrofitted_file
vulnerability_model_retrofitted.xml[hazard_outputs]
hazard_maps = true
uniform_hazard_spectra = true
poes = 0.1 0.02=

構造物の vulnerability モデル=

構造物の改修費の vulnerability モデル

ハザードマップを出力する
一様ハザードマップを出力する
超過確率 10%、2%

計算方法

計算方法は「3.4.1

WebUI による実行」（p.81~）および「3.4.2 コマンドラインによる計算実行方法」（p.90~）に示したとおりである。

3)

説明

計算結果

計算結果については省略する。

168