地震

Size: px
Start display at page:

Download "地震"

Transcription

1 静岡県立大学市民勉強会 2015 平成 27 年 3 月 14 日 地震の物理学 - 地震予知はなぜ難しいか? 大阪大学理学研究科 川村光

2 自己紹介 所属 : 大阪大学理学研究科宇宙地球科学専攻 兼務学振 学術システム研究センター主任研究員 ( 数物系 ) 専門分野 : 物性物理 ( 磁性 相転移 ) 統計物理 ( 多体系 相互作用系 ) 地震研究統計物理 物性物理的観点から ( モデル系の数値シミュレーション )

3 Part I 今日のお話 地震とは地震の2 面性地震観測と最新テクノロジー地震先行現象固着すべり現象としての地震 Part II 地震の統計物理モデルと数値シミュレーションバネ-ブロックモデル OFCモデル ( セルラーオートマトンモデル )

4 地震の古文書 東北沖地震と津波 関東大震災 : 被災後の銀座 兵庫県南部地震 : 高速道の倒壊

5 地震の世界地図

6 プレートテクトニクス 地球表面はプレートに覆われ 各プレートはお互いに動いている

7 横ずれ型プレート境界 ( 北米サンアンドレアス断層 )

8 日本付近のプレート収束型プレート境界

9 地震はプレート運動に伴う断層の 固着ー滑り不安定性 断層 断層 プレート プレート

10 地震の物理プロセスの探求 地震は プレートの運動によって駆動された断層の固着すべり ( 破砕 ) 現象 弾性反発理論 (H.F. Reid, 1910) 摩擦 破砕の物性物理 非平衡開放系の統計物理

11 地震現象の 2 面性 ー地震のジキルとハイド 決定論的 規則的 固有地震 アスペリティ 非決定論的 不規則的 臨界的 ( 自己組織化臨界現象 ;SOC)

12 地震現象の 2 面性 I 決定論的 規則的側面 固有地震 アスペリティ 特徴的スケールを持つ

13 海岸段丘地形 ( ニュージーランド )

14 固有地震 アスペリティ - 特徴的スケールの存在 釜石沖のアスぺリティ アスぺリティ概念

15 固有地震 アスペリティ - 特徴的スケールの存在 釜石沖のアスぺリティ 東北沖以前 東北沖以後

16 東海 南海トラフにおける大地震の歴史記録 過去 白鳳地震 仁和地震 永長地震 明応地震 宝永地震 安政南海地震 安政東海地震 現在 昭和南海地震 昭和東南海地震

17 地震の再帰 歪み率プレート境界 プレート内部 歪み e = 変化した長さ 元の長さ 地震時の歪み変化 地震再帰時間 プレート境界 プレート内部

18 地震現象の 2 面性 II 非決定論的 不規則的側面 臨界的 ( 自己組織化臨界現象 ;SOC) 特徴的スケールを持たない

19 地震の統計的性質 * 地震は複雑な現象であるが 多数のイベントを平均して得られる統計的性質には しばしば明瞭な規則性 再現性が現れる. * そのような統計的性質には しばしばべき乗則が現れる 例 : グーテンベルグ - リヒター則 大森則 現象に特徴的なスケールが存在しない 自己組織化臨界現象 (SOC) [P. Bak et al, 89] vs. 固有地震

20 グーテンベルグ - リヒター則 地震モーメント : M w 地震のエネルギー 大きさ M w の地震の起きる頻度 N (M w ) ~ a M w - B, B ~ 2/3 べき乗側の出現 特徴的なエネルギースケールの欠如 マグニチュード : m w イベント数 m w =(2/3)logM w 6 マグニチュード

21 ベキ乗則とスケール不変性 ベキ乗則 y = k x a x x = cx y y = c a y と変換する このとき と変換することにより y = k x a と 最初と同型の関係式が得られる 尺度 ( 単位 ) の変更により法則が不変

22 大森則 * 本震に引き続いて起きる余震の頻度は, 時間に関するべき乗則に従う n(t) = c / t a a ~1 t: 本震後の経過時間 * 本震に先立って起きる前震の頻度も しばしば時間に関するべき乗則に従う ( 逆大森則 ). t : 本震前の時間

23 濃尾地震の余震頻度の時間変化 (1891; M8.0) 1 日あたりの余震の数 濃尾地震 1891 で地表に現れた地震断層 本震後の日数 N(t) ~ 1/t

24 地震は 臨界的 か 固有的 か? ー地震のジキルとハイド 地震現象に特徴的スケールは 存在するか?

25 地震観測と最新テクノロジー GNSS(GPS) 深海掘削 ( ちきゅう ) 津波堆積物とボーリング

26 GNSS(GPS) Global Navigation Satellite System (GPS は米国による最初の GNSS) GEONET: 国土地理院による GPS 連続観測システムリアルタイムで地殻変動をモニター

27 GNSS により観測された地表の動き 東北沖地震前 東北沖地震時

28 2003 十勝沖地震の前後の地表の動き 余効滑り 本震 [by K. Heki]

29 GNSS によって観測された スロースリップイベント (SSE) スロー地震 [by K. Heki]

30 GNSS で観測した 琉球諸島近傍での地表の動き 33 個のスロースリップイベントと 3 個の通常地震 [K. Heki et al., 2008]

31 深海掘削船 : ちきゅう

32 南海トラフや東北沖での深海掘削

33 長期抗内観測装置によるリアルタイム観測 コア資料の回収と分析

34 津波堆積物とボーリング 浜松東部 太田川流域での津波堆積物調査 [by O. Fujiwara]

35 浜松東部地域の過去の大津波 1 2 白鳳地震仁和地震永長地震 3 明応地震 4 4 回の大地震に対応する津波堆積物 [by O. Fujiwara]

36 地震予知 地震が固有的 ( 固有的な地震 ) なら GNSS 等を利用して ある程度の予知は可能になりつつある地震が臨界的 ( 臨界的な地震 ) なら 実際的な予知は困難 現状では 個々の地震に関し その時間 場所 規模を信頼性を持って予知するのは 困難 地震先行現象の解明

37 地震先行現象 [ex.] 2014 April 1 チリ北部 Iquique ( イキケ ) 地震 (M W =8.1) [B.Schurr et al., Nature (2014) 512, 219]

38 チリ北部 2014 Iquique 地震 (M W =8.1) における地震先行現象 1877 Mw=8.8 以来の地震空白域 2007 国際的集中観測体制スタート 地震空白域 2014 April 1 本震 Mw=8.1 発生 2014 April, 3 最大余震 Mw= 個の前震 1400 個の余震 ナスカプレート南米レート

39 前震 余震 ロッキング 前震 余震分布 前震 本震 余震 最大余震 深さ ロッキング 本震エリアは中間的なロッキングのエリアに対応

40 本震までの前震頻度 地表変位の時間変化 マグニチュード イベント積算数 GPS 信号 緯度 本震前の地表変位の GPS データ 本震までの日数 経度 経度 本震 3 年前くらいより b 値の減少と合わせて前震活動が静穏化 本震 9 か月前くらいに元のレベルに復し 本震 18 日前からは急激に前震活動が増大

41 地震前のマグニチュード分布と b 値の変化 積算頻度 マグニチュード 本震エリアで本震 3 年前くらいから b 値が減少 b 値 本震までの日数

42 固着すべり現象 としての地震 速度状態依存摩擦則 地震のミクロとマクロ

43 地震の物理プロセスの探求 地震は プレートの運動によって駆動された 断層の固着すべり ( 破砕 ) 現象 摩擦 破砕の物性物理学 非平衡開放系の統計物理学

44 摩擦と構成則 静摩擦力 > 動摩擦力 固着すべり不安定性 動的不安定性の起源 地震の物理学にとって最も重要であるが 最も不確定性が大きい部分

45 摩擦構成則 --- 動的不安定性の起源 地震の物理学にとって最も重要であるが 最も不確定性が大きい部分 [ 例 ] 速度弱化摩擦則 [Carlson and Langer, 91] f 1 早くすべるほど 摩擦が弱くなる friction 1-s slope=2a velocity a: 速度弱化性の程度を表す ~ 他にも様々な構成則が提案されている ~

46 岩石摩擦の室内実験 静摩擦の時間変化 (Dieterich,1972) 静止摩擦は静止時間に対し対数的に増大 速度変化に伴う動摩擦の時間変化 (Dieterich,1981) 動摩擦には 速度変化に対し 直ちに反応する速度強化成分と 時間をかけて反応する速度弱化成分が有る 摩擦は速度だけではなく状態による複雑な力

47 速度状態依存摩擦則 Dieterich(1979)-Ruina(1983) 新たに状態変数 q i (t ) ( 摩擦面の状態を記述する現象論的変数 ) を導入し 摩擦面の状態に依存する摩擦構成則を定式化 F = c + a logv i + b logq i 直接的速度強化項状態依存項 摩擦不安定性を導く V i : 速度 ; c, a, b : 摩擦則を特徴付ける物質定数 q i (t ) の発展方程式 : スローネス則 dq i /dt = 1 - v i q i 熱的過程低速領域で良い記述 a-b > 0 : 速度強化 安定すべり ( クリープ ) a-b < 0 : 速度弱化 不安定固着すべり ( 地震 )

48 地震現象のミクロとマクロ 断層のジオメトリー地殻の弾性的性質断層での摩擦構成則プレート運動による外部駆動 ミクロ 多自由度の協力現象 巨視現象としての地震 マクロ

49 Part II 地震の統計物理的モデルと 数値シミュレーション

50 地震現象のミクロとマクロ 断層のジオメトリー地殻の弾性的性質断層での摩擦構成則 プレート運動による外部駆動 地震の統計物理モデルミクロ 多自由度の協力現象 数値計算機シミュレーション 巨視現象としての地震 マクロ

51 地震の統計モデルと計算機シミュレーション 様々な統計モデルが提案されている e.g., バネーブロックモデル (Burridge-Knopoff model) coupled-lattice モデル (OFC model) モデル系の数値シミュレーションの利点 [Burridge, Knopoff, 67] [Olami,Feder,Christensen, 92] - 極めて多数のイベントを生成し統計的に有意な解析が出来る - 摩擦構成則等の物性パラメータを制御出来る

52 地震の統計モデル バネ - ブロックモデル (Burridge-Knopoff モデル ) 断層の摩擦則 ( 構成則 ) パラメータや地殻の弾性パラメータと地震の統計的性質との関係

53 バネ - ブロックモデルの 数値シミュレーション 地震のサイズ ( マグニチュード ) 地震の破壊核形成過程 スロースリップ

54 1D モデルでの典型的なブロックの変位曲線 マグニチュードの定義 : m

55 破壊伝播の様子 I 2D モデル

56 破壊伝播の様子 II 2D モデル 複雑 異方的かつ段階的な破壊伝播 結果的にアスぺリティが連動して動いたように見える

57 バネーブロックモデルのマグニチュード分布 ー速度弱化摩擦則ー 1 次元 2 次元 臨界的 a: 摩擦弱化パラメータ 超臨界的 臨界的 亜臨界的

58 バネーブロックモデルのマグニチュード分布 ー速度状態依存擦則ー 2 次元 強い摩擦不安定性 (b 大 ) 弱い摩擦不安定性 (b 小 ) [S.Kakui & H.Kawamura, 11] [S.Kakui & H.Kawamura, 11] 臨界的 超臨界的 小地震 : 臨界的大地震 : 固有的 より固有的性格が顕著! 速度状態依存摩擦則に従う一様均一な断層での大地震は固有的な振る舞いを示す

59 バネーブロックモデルの連続極限 マグニチュード分布 ブロックサイズ d 0 ; l ~ 1/d バネ - ブロックモデルでは空間を離散化している 連続極限ではどうなっているだろうか? 速度状態依存摩擦則に従うバネ - ブロックモデルの連続極限は固有的

60 2D バネーブロックモデルにおける 本震前の b 値の時間変化 日本の地震カタログ (JUNEC) 本震前 通常時 本震前 通常時 Iquique 地震前震 B 値は本震前に減少する!

61 破壊域サイズの分布速度状態依存摩擦則 強い摩擦不安定性 連続極限 弱い摩擦不安定性 L 0 100~200 km 特徴的破壊域長さ L 0 一様均一な無限断層においても地震には固有の長さスケールがある?

62 平均変位と応力降下 平均変位 平均応力降下 平均変位 平均応力降下 摩擦不安的性大 摩擦不安的性小 平均変位と応力降下は破壊域サイズによらない

63 地震の破壊核形成過程 [ Ohnaka, 2000 ] Lc 高速破壊フェーズ L < Lsc 安定かつ準静的に破壊が進行 加速フェーズ Lsc < L < Lc 不安定かつ加速的に破壊が進行 初期フェーズ Lsc 震源 Lc < L 高速度で定常的に破壊が進行 破壊核の臨界サイズ : Lsc と Lc の 2 つ存在 バネ - ブロックモデルで破壊核形成過程はどのように再現され理解できるだろうか?

64 バネ ブロックモデルの核形成過程 b>b c or b<b c (L sc >1 or L sc < 1) で異なる振る舞い b c =2l 2 +1 摩擦不安定性弱 b< b c (L sc >1) 摩擦不安定性強 b> b c (L sc <1) block number 初期破壊フェーズを伴う block number block number 初期破壊フェーズを伴わない 連続極限 :l 常に初期破壊フェーズを伴う

65 駆動を停止させたときの振舞い L=10 で駆動を停止 L=11 で駆動を停止 step block number block number Lsc = 10 破壊が止まるか 止まらないか Lsc は破壊の安定 不安定に対応しているとみなせる

66 臨界破壊核長さ Lsc スティフネス 摩擦弱化率 動いているブロック数 Lsc 不安定性! スティフネス : 理論 摩擦弱化率理論 ~ b 物質パラメータ b と l で決まる Lsc : 理論 Lsc > 1 : b < bc = 2l 2 +1 弱い摩擦不安定性 Lsc < 1 : b > bc = 2l 2 +1 強い摩擦不安定性

67 核形成ダイナミクス L c 震源ブロックの運動 ( すべり速度 v) 初期フェーズ すべり速度 v 状態変数 L c L sc L sc q 加速フェーズ Vq L c L sc Lc : 震源ブロックのすべり速度が最大になる点

68 臨界破壊核長さ L sc & L c L c L sc 連続極限における臨界破壊核長さ L sc ~ (GL v s /s n wb) 1/2 c.f. 連続弾性体 Lsc ~ (GL / s n b) ~ a few km L c と地震のサイズは無相関 ~ 10 km

69 各フェーズの継続時間 L c 高速破壊フェーズ継続時間 :T γ 地震の大きさと各継続時間 v=10-4 から本震までの時間 :T β ' 加速フェーズ継続時間 :T β 破壊核形成過程の継続時間と地震の大きさには相関なし v = 10-4 L sc 初期破壊フェーズ継続時間 :T α 各継続時間のパラメータ依存性 T β T β : v* 依存性有りただし T β 収束傾向 ~ 数時間程度 顕著な相関なし

70 地震予知に役立つだろうか? 本震のマグニチュードを予知できるか? Lsc は本震のマグニチュードに依存しない Lc も依存性は弱い 本震の開始時間を予知できるか? T b, T b, T g は本震のマグニチュードにほとんど依存しない Lsc Lc Tg Tb Ta

71 地震概念の拡大ースロースリップ * 前駆すべり (pre-slip) 破壊核形成過程 * 余効すべり (after-slip) * スロー地震 * 非地震性クリープ * 深部低周波微動

72 余効すべり スロースリップイベント スロー地震 ( 本震 ) [Heki et al, 1997] 変位 ( 本震 ) ( 本震 ) 変位 displacement 1 year 1 year ( 年 ) 三陸はるか沖地震 (1994) 豊後水道スロースリップ ( )

73 バネ-ブロックモデルでの余効すべり スロー地震 a b 速度強化項の係数 a と速度弱化の係数 b の大小関係により 地震を以下の 3 つの種類に分類できる a < b 通常の高速破壊 a b 余効すべりを伴う高速破壊 a >> b スロー地震 領域 Ⅰ: a < b 領域 Ⅱ: a b 領域 Ⅲ: a >> b 高速破壊 余効すべりを伴う高速破壊 スロー地震 a < b b b < a a

74 t 領域 Ⅰ: a < b ~ 高速破壊と前駆すべり ~ 三陸はるか沖地震前後の変位 破壊伝播速度 [V r =3.6]~ [l =4] [Heki et al, 1997] 高速破壊フェーズ 前駆すべり ( 破壊核形成過程 ) 断層ブロックのすべり速度 すべり継続時間 ~1 程度 最大すべり速度 = 数十 ~ 数百 t 時間 t 領域 Ⅱ: block a number b ~ 余効すべりを伴う高速破壊 ~ すべり変位 高速破壊 余効すべり t 破壊伝播速度 [V r =3.0] ~ [l =4] 余効すべり すべり速度 最大すべり速度 v ~ 60 前駆すべり ( 破壊核形成過程 ) 余効すべり ~ 数百程度 t 時間

75 領域 Ⅲ: a >> b ~ スロー地震 ~ t すべり速度 最大すべり速度 ~0.027 すべり継続時間 ~ 数万程度 破壊伝播速度 [V r =0.0043] << [l =4] t 時間 破壊伝播速度 すべり速度が遅く 継続時間が長い スロー地震の特徴 block number すべり変位 豊後水道スロースリップ ( ) t 時間

76 部分破壊か全破壊か? N= パラメータ 領域 Ⅰ[a < b] : 高速破壊 a=9 a=3 変位 領域 u Ⅰ[ a < b ] : 部分破壊領域 Ⅱ[ a b ] : 全破壊領域 Ⅲ[a >> b] : 全破壊 a=9 block number 領域 Ⅱ[a b] : 余効すべりを伴う高速破壊 N=1600 変位 u a=20 部分破壊 block number 領域 Ⅲ[a >> b] : スロー地震 変位 u 全破壊 変位 u 全破壊 block number block number

77 破壊伝搬速度 すべり速度 すべり継続時間 弾性波速度 破壊伝播速度 V r 高速破壊地震 余効すべりを伴う地震 Ⅰ Ⅱ Ⅲ スロー地震 最大すべり速度 v 高速破壊地震 Ⅰ 余効すべりを伴う地震 Ⅱ スロー地震 Ⅲ すべり継続時間 T Ⅰ Ⅱ Ⅲ 高速破壊地震 余効すべりを伴う地震 スロー地震

78 OFC モデル ( セルラーオートマトンモデル ) の数値シミュレーション モデルのダイナミクス地震のサイズ 再帰固有性と臨界性アスペリティ長時間変動

79 OFC モデル サイズ (s)=2 の地震イベント 2 次元正方格子上の各格子点 ( サイト ) ストレス変数 f 0 (Ⅰ) 閾値 (=1) に達するまで一定の割合で全サイトの f を一様に増加 (Ⅱ) ストレスの移送 *4 ヶ所の最近接サイト f n,n f n,n + α f i,j * 閾値サイト f i, j 0 0 α 0.25 a < 0.25 非保存 a = 0.25 保存

80 初期状態 ( ランダム ) ストレス分布の時間発展 シミュレーション スタート α=0.20, L=256

81 10event ストレス分布の時間発展

82 10 2 event ストレス分布の時間発展

83 10 3 event ストレス分布の時間発展

84 10 4 event ストレス分布の時間発展

85 10 5 event ストレス分布の時間発展

86 10 6 event ストレス分布の時間発展

87 ストレス分布の時間発展 event ストレス値が揃ったパッチ状のドメインが境界に出現

88 event ストレス分布の時間発展

89 ストレス分布の時間発展 event 系の内部に周期 1-4α の s=1 イベントが出現

90 event ストレス分布の時間発展

91 ストレス分布の時間発展 event パッチ状のドメインが境界から系の内部に次第に侵入

92 event ストレス分布の時間発展

93 event ストレス分布の時間発展

94 event ストレス分布の時間発展

95 ストレス分布の時間発展 10 7 event パッチ状のドメインが内部まで浸透 以後は定常状態

96 10 8 event ストレス分布の時間発展

97 10 9 event ストレス分布の時間発展

98 event ストレス分布の時間発展

99 OFC モデルは SOC 的な振る舞いを示す [Z.Olami et al, 1992] [S. Hergarten et al, 2002] 前震の頻度 グーテンベルグリヒター則 余震の頻度 大森則 サイズ分布 [H. Kawamura et al, 2010]

100 局所的再帰時間分布 T =1-4α 周期的ピークに対応した地震のサイズ分布 P(T) T T 分布は周期的なピークとベキ乗則的な部分から成る 時間に関して周期的だが サイズに関しては臨界的 臨界性と固有性 ( 周期性 ) が共存

101 a=0.2 T = 1-4a = 0.2 OFC モデルにおけるアスペリティ的挙動 震源 地震 破壊域 ストレスプロット loading 震源 地震 アスペリティ 破壊域 同じ領域 ( アスペリティ ) が何度も繰り返し破壊

102 震源位置 震源破壊域 アスペリティ的イベント 非アスペリティ的イベント OFC モデルでは 震源は破壊域の端に位置する

103 1 回目 アスペリティ的な連続イベント列における破壊時のストレス値の時間発展 P(f) f アスペリティ的イベントと非アスペリティ的イベントにおけるストレス分布 P(f) f 14 回目 アスペリティ的イベントでは自己組織化的なストレス集中が起きる

104 アスペリティ サイクルの終焉 モード I: 分裂 モード II: 合体 アスぺリティは永遠には続かない - モデルは並進対称

105 アスペリティ的イベントと非アスペリティ的イベントの時間発展 (a) (b) 非アスペリティ的イベント アスペリティ的イベント

106 まとめ 地震の研究の現状を 地震の決定論的 固有的性格と非決定論的 臨界的性格の対比に着目しつつ 概観した * 固有的な性格の強い地震については GNSS 等の観測手段の発展により ある程度の予知が可能になりつつある一方 臨界的な性格の強い地震については現実的には 予知は未だ困難 代表的な地震の統計モデルである OFCモデルとバネ ブロックモデルの性質を 数値シミュレーションにより探査した * 固有的な振る舞いと臨界的な振る舞いの共存 前駆すべり ( 破壊核形成過程 ) や余効すべり スロー地震のスロースリップ等の多様なすべり現象が観測された

107 さいごに 現状では 予知はあくまで補助的に位置付け * 街や住居を地震に強い構造に作る * リアルタイムの警報システムを整備する ことが重要 一見迂遠なようでも 地震の本性を理解するための地道な分野横断的な基礎研究が大変重要 * 最新技術を駆使し 地震データを蓄積していく * 地震の本性を 物理学の視点から理解する

<4D F736F F F696E74202D AD482C682E882DC82C682DF90E096BE8E9197BF C C C816A2E B93C782DD8EE682E890EA97705D>

<4D F736F F F696E74202D AD482C682E882DC82C682DF90E096BE8E9197BF C C C816A2E B93C782DD8EE682E890EA97705D> 南海トラフの巨大地震モデル検討会中間とりまとめポイント はじめに Ⅰ 章 中間とりまとめの位置づけ 南海トラフの巨大地震モデルの想定震源域 想定津波波源域の設定の考え方や最終とりまとめに向けた検討内容等をとりまとめたもの 南海トラフの最大クラスの巨大な地震 津波に関する検討スタンス Ⅱ 章 これまでの対象地震 津波の考え方 過去数百年間に発生した地震の記録 (1707 年宝永地震以降の 5 地震 )

More information

別添資料 3 南海トラフ沿いの大規模地震の 予測可能性に関する調査部会 ( 報告 ) 南海トラフ沿いの大規模地震の予測可能性について 平成 25 年 5 月 -0-

別添資料 3 南海トラフ沿いの大規模地震の 予測可能性に関する調査部会 ( 報告 ) 南海トラフ沿いの大規模地震の予測可能性について 平成 25 年 5 月 -0- 別添資料 3 南海トラフ沿いの大規模地震の 予測可能性に関する調査部会 ( 報告 ) 南海トラフ沿いの大規模地震の予測可能性について 平成 25 年 5 月 -0- 1. はじめに東南海 南海地震は今世紀前半にも発生するおそれがあり 想定東海地震が発生していない現状を考慮すると 東海地震と東南海 南海地震が連動して発生する可能性が生じてきており 南海トラフ巨大地震対策検討ワーキンググループ( 主査

More information

本ワーキンググループにおけるこれまでの検討事項

本ワーキンググループにおけるこれまでの検討事項 本ワーキンググループにおけるこれまでの検討事項 資料 1 本ワーキンググループの検討内容 第 3 回資料 1 を一部修正 地震発生予測について ( 予測可能性調査部会において検討 ) 〇現状の地震発生予測の可能性 確度 予測可能性に関する科学的知見を整理 社会が混乱するおそれがある 4 つのケースについて検討 〇南海トラフにおけるリアルタイムモニタリング 南海トラフで発生している現象を分析 評価し

More information

Hazard_ pptx

Hazard_ pptx 1 南海トラフの巨大地震 : 新想定 予測可能性 長期評価 京都大学防災研究所橋本学 2 この 2 年間の主な所外での活動 日本地震学会東北地方太平洋沖地震対応臨時委員会委員 地震調査研究推進本部地震調査委員会長期評価部会海溝型分科会 ( その 2) 委員 内閣府南海トラフの巨大地震モデル検討会委員 内閣府南海トラフ沿いの大規模地震の予測可能性に関する調査部会委員 総合科学技術会議評価専門調査会 日本海溝海底地震津波観測網の整備及び緊急津波速報

More information

日本海溝海底地震津波観測網の整備と緊急津波速報 ( 仮称 ) システムの現状と将来像 < 日本海溝海底地震津波観測網の整備 > 地震情報 津波情報 その他 ( 研究活動に必要な情報等 ) 海底観測網の整備及び活用の現状 陸域と比べ海域の観測点 ( 地震計 ) は少ない ( 陸上 : 1378 点海域

日本海溝海底地震津波観測網の整備と緊急津波速報 ( 仮称 ) システムの現状と将来像 < 日本海溝海底地震津波観測網の整備 > 地震情報 津波情報 その他 ( 研究活動に必要な情報等 ) 海底観測網の整備及び活用の現状 陸域と比べ海域の観測点 ( 地震計 ) は少ない ( 陸上 : 1378 点海域 資料 2 総合科学技術会議評価専門調査会 日本海溝海底地震津波観測網の整備及び緊急津波速報 ( 仮称 ) に係るシステム開発 評価検討会 ( 第 2 回 ) 資料 平成 23 年 11 月 10 日 文部科学省 研究開発局地震 防災研究課 日本海溝海底地震津波観測網の整備と緊急津波速報 ( 仮称 ) システムの現状と将来像 < 日本海溝海底地震津波観測網の整備 > 地震情報 津波情報 その他 ( 研究活動に必要な情報等

More information

<4D F736F F F696E74202D208E9197BF A957A8E9197BF816A205B8CDD8AB B83685D>

<4D F736F F F696E74202D208E9197BF A957A8E9197BF816A205B8CDD8AB B83685D> 資料 2 内閣府における 長周期地震動の検討 ( 内閣府検討結果の概要 ) 1 平成 27 年 12 月 17 日 内閣府の公表資料一覧 (1) 南海トラフ沿いの巨大地震による長周期地震動に関する報告 (2) 南海トラフ沿いの巨大地震による長周期地震動に関する報告図表集 (3) 別冊 1-1 南海トラフ沿いの過去地震の強震断層モデル及び津波断層モデル (4) 別冊 1-2 南海トラフ沿いの過去地震の強震断層モデル

More information

Microsoft Word - 223回予知連表紙.doc

Microsoft Word - 223回予知連表紙.doc 資料 8 第 2 2 3 回 地震予知連絡会資料 2019 年 5 月 22 日 京都大学防災研究所 第 223 回地震予知連絡会提出資料 目 次 Ⅰ. 日向灘の地殻活動等 (3c) 1. 日向灘の地震 (2019 年 5 月 10 日,Mj = 6.3) について 1 2.GNSS 観測による日向灘の地震 (M6.3) の地殻変動 4 3.ETAS モデルを用いた南海トラフ 日向灘における地震活動の異常性の評価

More information

目 次 1. 想定する巨大地震 強震断層モデルと震度分布... 2 (1) 推計の考え方... 2 (2) 震度分布の推計結果 津波断層モデルと津波高 浸水域等... 8 (1) 推計の考え方... 8 (2) 津波高等の推計結果 時間差を持って地震が

目 次 1. 想定する巨大地震 強震断層モデルと震度分布... 2 (1) 推計の考え方... 2 (2) 震度分布の推計結果 津波断層モデルと津波高 浸水域等... 8 (1) 推計の考え方... 8 (2) 津波高等の推計結果 時間差を持って地震が 別添資料 1 南海トラフ巨大地震対策について ( 最終報告 ) ~ 南海トラフ巨大地震の地震像 ~ 平成 25 年 5 月 中央防災会議 防災対策推進検討会議 南海トラフ巨大地震対策検討ワーキンググループ 目 次 1. 想定する巨大地震... 1 2. 強震断層モデルと震度分布... 2 (1) 推計の考え方... 2 (2) 震度分布の推計結果... 2 3. 津波断層モデルと津波高 浸水域等...

More information

12-9 地震活動の再現と固有地震

12-9 地震活動の再現と固有地震 12-9 地震活動の再現と固有地震 Reappearance of Similar Seismicity Pattern preceding Characteristic Earthquakes 防災科学技術研究所松村正三 National Research Institute for Earth Science and Disaster Prevention 1. 茨城県沖の地震活動プレート境界固有地震の発生前には,

More information

日本海地震・津波調査プロジェクト

日本海地震・津波調査プロジェクト 17/9/29 (2-5-3) 構成岩 モデルの構築 1. 構成岩 モデルを推定 2. 地震発 層下限の推定 震源断層モデル構築へ貢献 28-1-2-5-3 29-1-2-5-3 横浜国 学 川正弘 l 然地震データ解析 地殻構成岩 と断層下限の推定 u 捕獲岩の弾性波速度と地震波速度構造を 較することで地殻構成岩 を推定 秋 県 ノ 潟産捕獲岩の弾性波速度と東北地 の地震波速度構造を 較 (Nishimoto

More information

図 東北地方太平洋沖地震以降の震源分布図 ( 福島第一 第二原子力発電所周辺 ) 図 3 東北地方太平洋沖地震前後の主ひずみ分布図 ( 福島第一 第二原子力発電所周辺 )

図 東北地方太平洋沖地震以降の震源分布図 ( 福島第一 第二原子力発電所周辺 ) 図 3 東北地方太平洋沖地震前後の主ひずみ分布図 ( 福島第一 第二原子力発電所周辺 ) 平成 3 年 8 月 30 日東京電力株式会社 平成 3 年東北地方太平洋沖地震を踏まえた新耐震指針に照らした既設発電用原子炉施設等の耐震安全性の評価結果の報告に係る 原子力安全 保安院における検討に際しての意見の追加への対応について ( 追加指示 ) に基づく報告 概要版 当社は 平成 3 年 3 月 日に発生した東北地方太平洋沖地震 (M9.0) 以降の地震の発生状況及び地殻変動 ( 地盤の動き

More information

新潟県中越沖地震を踏まえた地下構造特性調査結果および駿河湾の地震で敷地内の揺れに違いが生じた要因の分析状況について

新潟県中越沖地震を踏まえた地下構造特性調査結果および駿河湾の地震で敷地内の揺れに違いが生じた要因の分析状況について < 別紙 > 新潟県中越沖地震を踏まえた地下構造特性調査結果 および 駿河湾の地震で敷地内の揺れに違いが生じた要因の分析状況について 新潟県中越沖地震を踏まえた地下構造特性調査 地下構造特性にかかわる既往の調査結果の信頼性を確認するとともに 知見をより一層充実させるため 敷地および敷地周辺の地下構造特性の調査を実施しました 調査項目 1 微動アレイ観測 調査箇所 調査内容 敷地内および敷地周辺 :147

More information

Microsoft PowerPoint - システム創成学基礎2.ppt [互換モード]

Microsoft PowerPoint - システム創成学基礎2.ppt [互換モード] システム創成学基礎 - 観測と状態 - 古田一雄 システムの状態 個別の構成要素の状態の集合としてシステムの状態は記述できる 太陽系の状態 太陽の状態 s 0 = {x 0,y 0,z 0,u 0,v 0,w 0 } 水星の状態 s 1 = {x 1,y 1,z 1,u 1,v 1,w 1 } 金星の状態 s 2 = {x 2,y 2,z 2,u 2,v 2,w 2 } 太陽系の状態 S={s 0,s

More information

Microsoft Word - 概要版(案)_ docx

Microsoft Word - 概要版(案)_ docx 第 2 編地震による自然現象の予測 1 調査の条件 1.1 想定地震 1.1.1 想定地震の設定方針本調査は 沖縄県の陸地部および周辺海域で想定される大規模地震により予想される物的 人的被害の状況を総合的に把握し 災害対策の基礎資料とするものであり 解析のための想定地震は 以下の点を考慮して設定した 過去の調査と整合を保つため 過去の調査 ( 平成 21 年度沖縄県地震被害想定調査 平成 23 24

More information

Microsoft PowerPoint - 西村先生【30枚】 公開講座.pptx

Microsoft PowerPoint - 西村先生【30枚】 公開講座.pptx 西南日本のひずみ集中帯と内陸地震 京都大学 防災研究所 地震予知研究センター 西村 卓也 京都大学防災研究所公開講座 2017年10月11日 地殻変動とその観測方法 地殻変動の例(熊本地震 測位衛星 みちびき の軌道 JAXA GPS衛星の軌道 地殻に生ずる種々の変形 変位 隆起 沈 降 断層や褶曲などの地殻変動は通常と ても小さいため 測定するためには精密な 測量が必要 GNSS 衛星測位システム

More information

陦ィ邏・3

陦ィ邏・3 研 究 ニ ュ ー ス 地震波で覗いた マントル最下部まで沈んだ 表面地殻の岩石質 ロバート ゲラー 地球惑星科学専攻 教授 私たちの立っている地殻のもとには D" 層はマントル対流における熱境界層 行った 図 1 その結果 他の地域で 地球の全体積の 8 割を超える 岩石で であり そこでは温度の不均質や組成の の D 領域構造と異なる S 波速度の 構成されているマントル そしてさらに 分化の可能性が示唆されており

More information

Microsoft PowerPoint - 資料4-1.ppt [互換モード]

Microsoft PowerPoint - 資料4-1.ppt [互換モード] 資料 4-1 南海トラフの巨大地震モデル検討会第 1 回会合 東海地震に関する専門調査会 及び 東南海 南海地震等に関する専門調査会 における検討 東海地震に関する専門調査会における検討 想定震源域 基本的な考え : 昭和東南海地震の未破壊領域を震源域とする 1 北側の境界当該地域の最近の震源分布からみて両プレートの存在が明瞭である領域までとする 2 北西側の境界深さ約 30km より浅い領域とする

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 平成 29 年度防災セミナー開催日時 : 平成 30 年 2 月 6 日開催場所 : 発明会館ホール 南海トラフ地震に関連する情報 について 気象庁地震火山部地震予知情報課 池田雅也 内容 1. 大規模地震対策特別措置法と東海地震の予知 2. 南海トラフ沿いの地震観測 評価に基づく防災対応検討ワーキンググループ における検討結果 3. 南海トラフ地震に関連する情報 について 4. 最近の南海トラフ周辺の地殻活動

More information

スライド 1

スライド 1 国立大学附置研究所 センター長会議第 1 部会シンポジウム 激変する地球環境の現状と未来像 スロー地震の発見と 巨大地震発生予測に対する貢献 地震研究所小原一成 2014 年 10 月 9 日 15:00~15:45 @ 京王プラザホテル札幌プラザホール 激変 地震活動やプレートの動きが激変? 地震に対する理解が激変 東北地方太平洋沖地震 スロー地震 講演内容 東北地方太平洋沖地震 それまでの地震観

More information

スライド 1

スライド 1 非線形数理秋の学校 パターン形成の数理とその周辺 - 反応拡散方程式理論による時 空間パターンの解析を中心に - 2007 年 9 月 25 日 -27 日 モデル方程式を通してみるパターン解析ー進行波からヘリカル波の分岐を例としてー 池田勉 ( 龍谷大学理工学部 ) 講義概要, 講義資料, 講義中に使用する C 言語プログラムと初期値データ, ヘリカル波のアニメーションをウェブで公開しています :

More information

2011 年度第 41 回天文 天体物理若手夏の学校 2011/8/1( 月 )-4( 木 ) 星間現象 18b 初代星形成における水素分子冷却モデルの影響 平野信吾 ( 東京大学 M2) 1. Introduction 初代星と水素分子冷却ファーストスター ( 初代星, PopIII) は重元素を

2011 年度第 41 回天文 天体物理若手夏の学校 2011/8/1( 月 )-4( 木 ) 星間現象 18b 初代星形成における水素分子冷却モデルの影響 平野信吾 ( 東京大学 M2) 1. Introduction 初代星と水素分子冷却ファーストスター ( 初代星, PopIII) は重元素を 2011 年度第 41 回天文 天体物理若手夏の学校 2011/8/1( 月 )-4( 木 ) 星間現象 18b 初代星形成における水素分子冷却モデルの影響 平野信吾 ( 東京大学 M2) 1. Introduction 初代星と水素分子冷却ファーストスター ( 初代星, PopIII) は重元素を含まない原始ガスから形成される 宇宙で最初に誕生する星である 初代星はその後の星形成や再電離など宇宙初期の天文現象に強く関係し

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション PID 制御の基礎 ON/OFF 制御 PID 制御 P 制御 過渡特性を改善しよう PD 制御と P-D 制御 定常特性を改善しよう PI-D 制御 4.2 節 I-PD 制御 角度制御実験装置 0 [deg] 30 [deg] 角度制御実験装置 目標値 コントローラ ( マイコン ) アクチュエータ (DC モータ ) 制御対象 ( アーム ) 角度 センサ ( ロータリエンコーダ ) ON/OFF

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 復習 ) 時系列のモデリング ~a. 離散時間モデル ~ y k + a 1 z 1 y k + + a na z n ay k = b 0 u k + b 1 z 1 u k + + b nb z n bu k y k = G z 1 u k = B(z 1 ) A(z 1 u k ) ARMA モデル A z 1 B z 1 = 1 + a 1 z 1 + + a na z n a = b 0

More information

以下 変数の上のドットは時間に関する微分を表わしている (ex. 2 dx d x x, x 2 dt dt ) 付録 E 非線形微分方程式の平衡点の安定性解析 E-1) 非線形方程式の線形近似特に言及してこなかったが これまでは線形微分方程式 ( x や x, x などがすべて 1 次で なおかつ

以下 変数の上のドットは時間に関する微分を表わしている (ex. 2 dx d x x, x 2 dt dt ) 付録 E 非線形微分方程式の平衡点の安定性解析 E-1) 非線形方程式の線形近似特に言及してこなかったが これまでは線形微分方程式 ( x や x, x などがすべて 1 次で なおかつ 以下 変数の上のドットは時間に関する微分を表わしている (e. d d, dt dt ) 付録 E 非線形微分方程式の平衡点の安定性解析 E-) 非線形方程式の線形近似特に言及してこなかったが これまでは線形微分方程式 ( や, などがすべて 次で なおかつそれらの係数が定数であるような微分方程式 ) に対して安定性の解析を行ってきた しかしながら 実際には非線形の微分方程式で記述される現象も多く存在する

More information

紀伊半島 ~ 愛知県における歪 傾斜 地下水観測結果 (2013 年 9 月 ~11 月 ) 産業技術総合研究所 産業技術総合研究所 ( 産総研 ) の紀伊半島 ~ 愛知県の10 観測点 ( 図 1) における2013 年 9 月 1 日から11 月 14 日までの歪 傾斜 地下水および気象観測結果を図 2~19に示す. また, 産総研がエンベロープ相関法を用いて決定した, 深部低周波微動の震源時空間分布および個数を図

More information

GMT v3.4.6 Document from psxy

GMT v3.4.6 Document from psxy 平成 年 月 日気象庁地震火山部 最近の東海地域とその周辺の地殻活動 現在のところ 東海地震に直ちに結びつくとみられる変化は観測していません. 地震の観測状況 月 日から 日 及び 日から 日にかけて 愛知県のプレート境界付近を震源とする深部低周波地震 ( 微動 ) を観測しています また 月 日以降 長野県のプレート境界付近を震源とする深部低周波地震 ( 微動 ) を観測しています. 地殻変動の観測状況

More information

2009 年 11 月 16 日版 ( 久家 ) 遠地 P 波の変位波形の作成 遠地 P 波の変位波形 ( 変位の時間関数 ) は 波線理論をもとに P U () t = S()* t E()* t P() t で近似的に計算できる * は畳み込み積分 (convolution) を表す ( 付録

2009 年 11 月 16 日版 ( 久家 ) 遠地 P 波の変位波形の作成 遠地 P 波の変位波形 ( 変位の時間関数 ) は 波線理論をもとに P U () t = S()* t E()* t P() t で近似的に計算できる * は畳み込み積分 (convolution) を表す ( 付録 遠地 波の変位波形の作成 遠地 波の変位波形 ( 変位の時間関数 ) は 波線理論をもとに U () t S() t E() t () t で近似的に計算できる は畳み込み積分 (convolution) を表す ( 付録 参照 ) ここで St () は地震の断層運動によって決まる時間関数 1 E() t は地下構造によって生じる種々の波の到着を与える時間関数 ( ここでは 直達 波とともに 震源そばの地表での反射波や変換波を与える時間関数

More information

紀伊半島 ~ 愛知県における歪 傾斜 地下水観測結果 (2012 年 9 月 ~11 月 ) 産業技術総合研究所 産業技術総合研究所 ( 産総研 ) の紀伊半島 ~ 愛知県の10 観測点 ( 図 1) における2012 年 9 月 1 日から11 月 15 日までの歪 傾斜 地下水および気象観測結果を図 2~19に示す. また, 産総研がエンベロープ相関法を用いて決定した, 深部低周波微動の震源時空間分布および個数を図

More information

Microsoft PowerPoint - LectureB1handout.ppt [互換モード]

Microsoft PowerPoint - LectureB1handout.ppt [互換モード] 本講義のスコープ 都市防災工学 後半第 回 : イントロダクション 千葉大学大学院工学研究科建築 都市科学専攻都市環境システムコース岡野創 耐震工学の専門家として知っていた方が良いが 敷居が高く 入り口で挫折しがちな分野をいくつか取り上げて説明 ランダム振動論 地震波形に対する構造物応答の理論的把握 減衰と地震応答 エネルギーバランス 地震動の各種スペクトルの相互関係 震源モデル 近年では震源モデルによる地震動予測が良く行われている

More information

プレス発表資料 平成 27 年 3 月 10 日独立行政法人防災科学技術研究所 インドネシア フィリピン チリにおけるリアルタイム 津波予測システムを公開 独立行政法人防災科学技術研究所 ( 理事長 : 岡田義光 以下 防災科研 ) は インドネシア フィリピン チリにおけるリアルタイム地震パラメー

プレス発表資料 平成 27 年 3 月 10 日独立行政法人防災科学技術研究所 インドネシア フィリピン チリにおけるリアルタイム 津波予測システムを公開 独立行政法人防災科学技術研究所 ( 理事長 : 岡田義光 以下 防災科研 ) は インドネシア フィリピン チリにおけるリアルタイム地震パラメー プレス発表資料 平成 27 年 3 月 10 日独立行政法人防災科学技術研究所 インドネシア フィリピン チリにおけるリアルタイム 津波予測システムを公開 独立行政法人防災科学技術研究所 ( 理事長 : 岡田義光 以下 防災科研 ) は インドネシア フィリピン チリにおけるリアルタイム地震パラメータ推定システム (SWIFT) と連動した津波予測システムを公開します 今回公開するのは SWIFT

More information

多次元レーザー分光で探る凝縮分子系の超高速動力学

多次元レーザー分光で探る凝縮分子系の超高速動力学 波動方程式と量子力学 谷村吉隆 京都大学理学研究科化学専攻 http:theochem.kuchem.kyoto-u.ac.jp TA: 岩元佑樹 iwamoto.y@kuchem.kyoto-u.ac.jp ベクトルと行列の作法 A 列ベクトル c = c c 行ベクトル A = [ c c c ] 転置ベクトル T A = [ c c c ] AA 内積 c AA = [ c c c ] c =

More information

             論文の内容の要旨

             論文の内容の要旨 論文の内容の要旨 論文題目 Superposition of macroscopically distinct states in quantum many-body systems ( 量子多体系におけるマクロに異なる状態の重ね合わせ ) 氏名森前智行 本論文では 量子多体系におけるマクロに異なる状態の重ねあわせを研究する 状態の重ね合わせ というのは古典論には無い量子論独特の概念であり 数学的には

More information

線形粘弾性 a.応力緩和とマクスウェル模型

線形粘弾性 a.応力緩和とマクスウェル模型 地殻 上部マントルのレオロジー 高温 長時間では, 岩石は流動する. 線形粘弾性体 ( 歪速度はせん断応力に比例 ) 流体力学 ( ナビエ - ストークスの式 ) べき乗型流動則 ( 歪速度はせん断応力のべき乗に比例 ) 比例係数の温度依存性が大きい. 参考書岩波講座地球科学 2 地球の物質科学 Ⅰ レオロジーと地球科学唐戸俊一郎グローバルテクトニクス杉村新 色々な深度の断層岩 ( 産業技術総合研究所,

More information

国土技術政策総合研究所 研究資料

国土技術政策総合研究所 研究資料 3. 解析モデルの作成汎用ソフトFEMAP(Ver.9.0) を用いて, ダムおよび基礎岩盤の有限要素メッシュを8 節点要素により作成した また, 貯水池の基本寸法および分割数を規定し,UNIVERSE 2) により差分メッシュを作成した 3.1 メッシュサイズと時間刻みの設定基準解析結果の精度を確保するために, 堤体 基礎岩盤 貯水池を有限要素でモデル化する際に, 要素メッシュの最大サイズならびに解析時間刻みは,

More information

スライド 1

スライド 1 相対論的プラズマにおける PIC シミュレーションに伴う数値チェレンコフ不安定の特性ついて 宇宙物理学研究室 4 年池谷直樹 研究背景と目的 0 年 Ie Cube 国際共同実験において超高エネルギーニュートリノを検出 780Tev-5.6PeV 890TeV-8.5PeV 相互作用が殆んど起こらないため銀河磁場による軌道の湾曲が無く 正確な到来方向の情報 を得られる可能性がある ニュートリノから高エネルギー宇宙線の起源を追う

More information

Chap2.key

Chap2.key . f( ) V (V V ) V e + V e V V V V ( ) V V ( ) E. - () V (0 ) () V (0 ) () V (0 ) (4) V ( ) E. - () V (0 ) () V (0 ) O r θ ( ) ( ) : (r θ) : { r cos θ r sn θ { r + () V (0 ) (4) V ( ) θ θ arg( ) : π π

More information

Probit , Mixed logit

Probit , Mixed logit Probit, Mixed logit 2016/5/16 スタートアップゼミ #5 B4 後藤祥孝 1 0. 目次 Probit モデルについて 1. モデル概要 2. 定式化と理解 3. 推定 Mixed logit モデルについて 4. モデル概要 5. 定式化と理解 6. 推定 2 1.Probit 概要 プロビットモデルとは. 効用関数の誤差項に多変量正規分布を仮定したもの. 誤差項には様々な要因が存在するため,

More information

スライド 1

スライド 1 P.1 NUMO の確率論的評価手法の開発 原子力学会バックエンド部会第 30 回 バックエンド 夏期セミナー 2014 年 8 月 7 日 ( 木 ) ビッグパレットふくしま 原子力発電環境整備機構技術部後藤淳一 確率論的アプローチの検討の背景 P.2 プレート運動の安定性を前提に, 過去 ~ 現在の自然現象の変動傾向を将来に外挿し, 地層の著しい変動を回避 ( 決定論的アプローチ ) 回避してもなお残る不確実性が存在

More information

東日本大震災 鳴らされていた警鐘

東日本大震災 鳴らされていた警鐘 .5m 9. 311 11 11869 15 3 1131116 13kmkm 9. 7 6 5 311 M7.7 M7.5M7. 7 M7.1 J A X A 3 km M8. 5 1 1 1319 17 7 6689 15853 855 1936 8 87km 8 16 5 11 6 5 311 13kmkm M9. 5km 1m 1896 1933 31m 1 km8m 63mm M7.3 M9.

More information

物理学 II( 熱力学 ) 期末試験問題 (2) 問 (2) : 以下のカルノーサイクルの p V 線図に関して以下の問題に答えなさい. (a) "! (a) p V 線図の各過程 ( ) の名称とそのと (& きの仕事 W の面積を図示せよ. # " %&! (' $! #! " $ %'!!!

物理学 II( 熱力学 ) 期末試験問題 (2) 問 (2) : 以下のカルノーサイクルの p V 線図に関して以下の問題に答えなさい. (a) ! (a) p V 線図の各過程 ( ) の名称とそのと (& きの仕事 W の面積を図示せよ. #  %&! (' $! #!  $ %'!!! 物理学 II( 熱力学 ) 期末試験問題 & 解答 (1) 問 (1): 以下の文章の空欄に相応しい用語あるいは文字式を記入しなさい. 温度とは物体の熱さ冷たさを表す概念である. 物体は外部の影響を受けなければ, 十分な時間が経過すると全体が一様な温度の定常的な熱平衡状態となる. 物体 と物体 が熱平衡にあり, 物体 と物体 が熱平衡にあるならば, 物体 と物体 も熱平衡にある. これを熱力学第 0

More information

資料 1 南海トラフの巨大地震モデル検討会 第 6 回会合 深部地盤モデルの作成の考え方 平成 23 年 12 月 12 日 1. 震度分布の推計方法 中央防災会議 (2003) 1 は 強震波形計算によって求められた地表の震度と経験的手法によって求められた地表の震度を比較検討し 強震波形計算による結果を主に それにより表現できていないところについては 経験的手法による結果も加えて 最終的な震度分布を求めている

More information

風力発電インデックスの算出方法について 1. 風力発電インデックスについて風力発電インデックスは 気象庁 GPV(RSM) 1 局地気象モデル 2 (ANEMOS:LAWEPS-1 次領域モデル ) マスコンモデル 3 により 1km メッシュの地上高 70m における 24 時間の毎時風速を予測し

風力発電インデックスの算出方法について 1. 風力発電インデックスについて風力発電インデックスは 気象庁 GPV(RSM) 1 局地気象モデル 2 (ANEMOS:LAWEPS-1 次領域モデル ) マスコンモデル 3 により 1km メッシュの地上高 70m における 24 時間の毎時風速を予測し 風力発電インデックスの算出方法について 1. 風力発電インデックスについて風力発電インデックスは 気象庁 GPV(RSM) 1 局地気象モデル 2 (ANEMOS:LAWEPS-1 次領域モデル ) マスコンモデル 3 により 1km メッシュの地上高 70m における 24 時間の毎時風速を予測し 2000kW 定格風車の設備利用率として表示させたものです 数値は風車の定格出力 (2000kW)

More information

詳細な説明 2016 年 4 月 16 日に発生した熊本地震 ( マグニチュード (M) 7.3)( 図 1) は 熊本県 大分県を中心に甚大な被害をもたらしました 九州地方は 北東 - 南西方向に縦走する 別府 - 島原地溝帯 と呼ばれる顕著な地殻の裂け目によって特徴づけられます 別府 - 島原地

詳細な説明 2016 年 4 月 16 日に発生した熊本地震 ( マグニチュード (M) 7.3)( 図 1) は 熊本県 大分県を中心に甚大な被害をもたらしました 九州地方は 北東 - 南西方向に縦走する 別府 - 島原地溝帯 と呼ばれる顕著な地殻の裂け目によって特徴づけられます 別府 - 島原地 平成 30 年 10 月 23 日 報道機関各位 東北大学大学院理学研究科 九州を南北に分裂させる地溝帯の構造を解明 -2016 年熊本地震の発生とも関連 - 発表のポイント 別府 - 島原地溝帯 周辺の地下構造を 初めて高分解能で解明した この地溝帯は 活火山下の熱いマントル上昇流 沖縄トラフ 及び中央 構造線の影響が複合して形成されたと考えられる 2016 年熊本地震の発生には この地溝帯の不均質構造と震源直下の水の挙動が影響した

More information

第 109 回 火山噴火予知連絡会資料 2008 年 2 月 15 日 東北大学大学院理学研究科

第 109 回 火山噴火予知連絡会資料 2008 年 2 月 15 日 東北大学大学院理学研究科 第 19 回 火山噴火予知連絡会資料 28 年 2 月 15 日 学大学院理学研究科 (a) 5 4 微動活動と震源分布 昨年の夏以降, 岩手火山における火山性微動 低周波地震の活動がやや活発になっている. Tremor & LF at Iwate Volcano DURATION DURATION (MIN) 3 2 AMPLITUDE ( 1microm/s) 1 (b) 1995 1996 1997

More information

数値計算で学ぶ物理学 4 放物運動と惑星運動 地上のように下向きに重力がはたらいているような場においては 物体を投げると放物運動をする 一方 中心星のまわりの重力場中では 惑星は 円 だ円 放物線または双曲線を描きながら運動する ここでは 放物運動と惑星運動を 運動方程式を導出したうえで 数値シミュ

数値計算で学ぶ物理学 4 放物運動と惑星運動 地上のように下向きに重力がはたらいているような場においては 物体を投げると放物運動をする 一方 中心星のまわりの重力場中では 惑星は 円 だ円 放物線または双曲線を描きながら運動する ここでは 放物運動と惑星運動を 運動方程式を導出したうえで 数値シミュ 数値計算で学ぶ物理学 4 放物運動と惑星運動 地上のように下向きに重力がはたらいているような場においては 物体を投げると放物運動をする 一方 中心星のまわりの重力場中では 惑星は 円 だ円 放物線または双曲線を描きながら運動する ここでは 放物運動と惑星運動を 運動方程式を導出したうえで 数値シミュレーションによって計算してみる 4.1 放物運動一様な重力場における放物運動を考える 一般に質量の物体に作用する力をとすると運動方程式は

More information

Microsoft PowerPoint - suta.ppt [互換モード]

Microsoft PowerPoint - suta.ppt [互換モード] 弾塑性不飽和土構成モデルの一般化と土 / 水連成解析への適用 研究の背景 不飽和状態にある土構造物の弾塑性挙動 ロックフィルダム 道路盛土 長期的に正確な予測 不飽和土弾塑性構成モデル 水頭変動 雨水の浸潤 乾湿の繰り返し 土構造物の品質変化 不飽和土の特徴的な力学特性 不飽和土の特性 サクション サクション s w C 飽和度が低い状態 飽和度が高い状態 サクションの効果 空気侵入値 B. サクション増加

More information

<4D F736F F D208C46967B926E906B82CC96C6906B8C9A95A8899E939A89F090CD>

<4D F736F F D208C46967B926E906B82CC96C6906B8C9A95A8899E939A89F090CD> 平成 29 年 9 月 1 日 観測記録に基づく免震住宅の地震応答解析 - 216 年熊本地震 - 1. はじめに 216 年 4 月 16 日 1 時 25 分に発生した熊本地震は マグニチュード 7.3 最大震度 7 と発表されています 防災科学技術研究所では 強震観測網 (K-NET KiK-net) により観測されたデータを公開データしています この観測地震動を用いて 免震住宅の地震応答解析を実施しました

More information

Microsoft Word - note02.doc

Microsoft Word - note02.doc 年度 物理化学 Ⅱ 講義ノート. 二原子分子の振動. 調和振動子近似 モデル 分子 = 理想的なバネでつながった原子 r : 核間距離, r e : 平衡核間距離, : 変位 ( = r r e ), k f : 力の定数ポテンシャルエネルギー ( ) k V = f (.) 古典運動方程式 [ 振動数 ] 3.3 d kf (.) dt μ : 換算質量 (m, m : 原子, の質量 ) mm

More information

1 想定地震の概要南海トラフで発生する地震は 多様な地震発生のパターンが考えられることから 次の地震の震源域の広がりを正確に予測することは 現時点の科学的知見では困難です そのため 本市では 南海トラフで発生する地震として 次の2つの地震を想定して被害予測調査を行いました (1) 過去の地震を考慮し

1 想定地震の概要南海トラフで発生する地震は 多様な地震発生のパターンが考えられることから 次の地震の震源域の広がりを正確に予測することは 現時点の科学的知見では困難です そのため 本市では 南海トラフで発生する地震として 次の2つの地震を想定して被害予測調査を行いました (1) 過去の地震を考慮し 1 想定地震の概要南海トラフで発生する地震は 多様な地震発生のパターンが考えられることから 次の地震の震源域の広がりを正確に予測することは 現時点の科学的知見では困難です そのため 本市では 南海トラフで発生する地震として 次の2つの地震を想定して被害予測調査を行いました 南海トラフ沿いでは 宝永地震 (1707) 安政東海 安政南海地震(1854) 昭 和東南海 (1944) 昭和南海地震(1946)

More information

Microsoft Word - thesis.doc

Microsoft Word - thesis.doc 剛体の基礎理論 -. 剛体の基礎理論初めに本論文で大域的に使用する記号を定義する. 使用する記号トルク撃力力角運動量角速度姿勢対角化された慣性テンソル慣性テンソル運動量速度位置質量時間 J W f F P p .. 質点の並進運動 質点は位置 と速度 P を用いる. ニュートンの運動方程式 という状態を持つ. 但し ここでは速度ではなく運動量 F P F.... より質点の運動は既に明らかであり 質点の状態ベクトル

More information

Microsoft PowerPoint _T3 [互換モード]

Microsoft PowerPoint _T3 [互換モード] 課題研究 T3( 固体圏 ) この課題では さまざまな観測 調査 実験的 理論的な研究 各種の数値シミュレーションなどにより 固体地球の構造 物性の解明 固体地球でのさまざまな時間 空間スケールでの変動やそのメカニズムの解明を目指す T3 がカバーする分野地震学 測地学 活構造学 地球熱学 火山学 担当教員 固体地球物理学講座平原和朗 ( 地震 ) 中西一郎 ( 地震 ) 久家慶子 ( 地震 ) 福田洋一

More information

昭和東南海地震以降初めての中規模プレート境界地震発生過程をとらえた - リアルタイム観測網の統合的データ解析 - 1. 概要国立研究開発法人海洋研究開発機構 ( 理事長平朝彦 以下 JAMSTEC という ) 地震津波海域観測研究開発センターの荒木英一郎主任技術研究員 鹿児島大学大学院理工学研究科の

昭和東南海地震以降初めての中規模プレート境界地震発生過程をとらえた - リアルタイム観測網の統合的データ解析 - 1. 概要国立研究開発法人海洋研究開発機構 ( 理事長平朝彦 以下 JAMSTEC という ) 地震津波海域観測研究開発センターの荒木英一郎主任技術研究員 鹿児島大学大学院理工学研究科の 昭和東南海地震以降初めての中規模プレート境界地震発生過程をとらえた - リアルタイム観測網の統合的データ解析 - 1. 概要国立研究開発法人海洋研究開発機構 ( 理事長平朝彦 以下 JAMSTEC という ) 地震津波海域観測研究開発センターの荒木英一郎主任技術研究員 鹿児島大学大学院理工学研究科の小林励司准教授 京都大学防災研究所地震予知研究センターの木下千裕博士課程大学院生らは テキサス大学 ペンシルバニア州立大学等と共同で

More information

<4D F736F F D E9197BF31817A975C91AA907D C4816A82C982C282A282C491CE8FDB926E906B82CC90E096BE2E646F63>

<4D F736F F D E9197BF31817A975C91AA907D C4816A82C982C282A282C491CE8FDB926E906B82CC90E096BE2E646F63> 資料 1 < 新たな津波浸水予測図 ( 素案 ) について > 今後の津波対策を構築するにあたっては 二つのレベルの津波を想定する 最大クラスの津波 : 住民避難を柱とした総合的防災対策を構築する上で設定する津波であり 発生頻度は極めて低いものの 発生すれば甚大な被害をもたらす最大クラスの津波 頻度の高い津波 : 防潮堤など構造物によって津波の内陸への侵入を防ぐ海岸保全施設等の整備を行う上で想定する津波

More information

第6章 実験モード解析

第6章 実験モード解析 第 6 章実験モード解析 6. 実験モード解析とは 6. 有限自由度系の実験モード解析 6.3 連続体の実験モード解析 6. 実験モード解析とは 実験モード解析とは加振実験によって測定された外力と応答を用いてモードパラメータ ( 固有振動数, モード減衰比, 正規固有モードなど ) を求める ( 同定する ) 方法である. 力計 試験体 変位計 / 加速度計 実験モード解析の概念 時間領域データを利用する方法

More information

パソコンシミュレータの現状

パソコンシミュレータの現状 第 2 章微分 偏微分, 写像 豊橋技術科学大学森謙一郎 2. 連続関数と微分 工学において物理現象を支配する方程式は微分方程式で表されていることが多く, 有限要素法も微分方程式を解く数値解析法であり, 定式化においては微分 積分が一般的に用いられており. 数学の基礎知識が必要になる. 図 2. に示すように, 微分は連続な関数 f() の傾きを求めることであり, 微小な に対して傾きを表し, を無限に

More information

相対性理論入門 1 Lorentz 変換 光がどのような座標系に対しても同一の速さ c で進むことから導かれる座標の一次変換である. (x, y, z, t ) の座標系が (x, y, z, t) の座標系に対して x 軸方向に w の速度で進んでいる場合, 座標系が一次変換で関係づけられるとする

相対性理論入門 1 Lorentz 変換 光がどのような座標系に対しても同一の速さ c で進むことから導かれる座標の一次変換である. (x, y, z, t ) の座標系が (x, y, z, t) の座標系に対して x 軸方向に w の速度で進んでいる場合, 座標系が一次変換で関係づけられるとする 相対性理論入門 Lorentz 変換 光がどのような座標系に対しても同一の速さ で進むことから導かれる座標の一次変換である. x, y, z, t ) の座標系が x, y, z, t) の座標系に対して x 軸方向に w の速度で進んでいる場合, 座標系が一次変換で関係づけられるとすると, x A x wt) y y z z t Bx + Dt 弨弱弩弨弲弩弨弳弩弨弴弩 が成立する. 図 : 相対速度

More information

Earthquakes 101

Earthquakes 101 東日本大震災で地震学者が学んだこと 青木陽介 東京大学地震研究所 Email: yaoki@eri.u-tokyo.ac.jp 2011 年 11 月 12 日中央大学白門物理同窓会 地震とは? 地下の断層の摩擦すべりによって地震が発生する. 岩石の温度が約 350 度以下でないと摩擦すべりをおこせない. 地震発生は地下のごく浅い ( 地球半径 6371km に比べれば ) 部分に限定される. なぜ日本には地震が多いのか?

More information

PowerPoint Presentation

PowerPoint Presentation Non-linea factue mechanics き裂先端付近の塑性変形 塑性域 R 破壊進行領域応カ特異場 Ω R R Hutchinson, Rice and Rosengen 全ひずみ塑性理論に基づいた解析 現段階のひずみは 除荷がないとすると現段階の応力で一義的に決まる 単純引張り時の応カーひずみ関係 ( 構成方程式 ): ( ) ( ) n () y y y ここで α,n 定数, /

More information

CAEシミュレーションツールを用いた統計の基礎教育 | (株)日科技研

CAEシミュレーションツールを用いた統計の基礎教育 | (株)日科技研 CAE シミュレーションツール を用いた統計の基礎教育 ( 株 ) 日本科学技術研修所数理事業部 1 現在の統計教育の課題 2009 年から統計教育が中等 高等教育の必須科目となり, 大学でも問題解決ができるような人材 ( 学生 ) を育てたい. 大学ではコンピューター ( 統計ソフトの利用 ) を重視した教育をより積極的におこなうのと同時に, 理論面もきちんと教育すべきである. ( 報告 数理科学分野における統計科学教育

More information

Microsoft PowerPoint - 発表II-3原稿r02.ppt [互換モード]

Microsoft PowerPoint - 発表II-3原稿r02.ppt [互換モード] 地震時の原子力発電所燃料プールからの溢水量解析プログラム 地球工学研究所田中伸和豊田幸宏 Central Research Institute of Electric Power Industry 1 1. はじめに ( その 1) 2003 年十勝沖地震では 震源から離れた苫小牧地区の石油タンクに スロッシング ( 液面揺動 ) による火災被害が生じた 2007 年中越沖地震では 原子力発電所内の燃料プールからの溢水があり

More information

Microsoft PowerPoint - 04_01_text_UML_03-Sequence-Com.ppt

Microsoft PowerPoint - 04_01_text_UML_03-Sequence-Com.ppt システム設計 (1) シーケンス図 コミュニケーション図等 1 今日の演習のねらい 2 今日の演習のねらい 情報システムを構成するオブジェクトの考え方を理解す る 業務プロセスでのオブジェクトの相互作用を考える シーケンス図 コミュニケーション図を作成する 前回までの講義システム開発の上流工程として 要求仕様を確定パソコンを注文するまでのユースケースユースケースから画面の検討イベントフロー アクティビティ図

More information

Autodesk Inventor Skill Builders Autodesk Inventor 2010 構造解析の精度改良 メッシュリファインメントによる収束計算 予想作業時間:15 分 対象のバージョン:Inventor 2010 もしくはそれ以降のバージョン シミュレーションを設定する際

Autodesk Inventor Skill Builders Autodesk Inventor 2010 構造解析の精度改良 メッシュリファインメントによる収束計算 予想作業時間:15 分 対象のバージョン:Inventor 2010 もしくはそれ以降のバージョン シミュレーションを設定する際 Autodesk Inventor Skill Builders Autodesk Inventor 2010 構造解析の精度改良 メッシュリファインメントによる収束計算 予想作業時間:15 分 対象のバージョン:Inventor 2010 もしくはそれ以降のバージョン シミュレーションを設定する際に 収束判定に関するデフォルトの設定をそのまま使うか 修正をします 応力解析ソルバーでは計算の終了を判断するときにこの設定を使います

More information

2019 年1月3日熊本県熊本地方の地震の評価(平成31年2月12日公表)

2019 年1月3日熊本県熊本地方の地震の評価(平成31年2月12日公表) 平成 3 年 月 日地震調査研究推進本部地震調査委員会 9 年 月 3 日熊本県熊本地方の地震の評価 月 3 日 8 時 分に熊本県熊本地方の深さ約 km でマグニチュード (M)5. の地震が発生した この地震により熊本地方の震央近傍で最大震度 6 弱を観測した その後 北西 - 南東方向に延びる約 5 kmの領域で地震活動が減衰しつつも継続している 月 日までに発生した最大の地震は 月 6 日に深さ約

More information

WTENK5-6_26265.pdf

WTENK5-6_26265.pdf 466 2014年秋季 極域 寒冷域研究連絡会 の報告 海 カラ海 北大西洋 北米大陸の北部 東アジアで が多重に見られることが多い 南極昭和基地 69.0 S, 寒気質量の減少傾向が 中央シベリアの内陸部とベー 39.6 E における PANSY レーダー Sato et al.2014 リング海で寒気質量の増加傾向が5つの再解析データ のデータは このような小さな に共通して見られた 中央シベリアの内陸部の寒気質

More information

概論 : 人工の爆発と自然地震の違い ~ 波形の違いを調べる前に ~ 人為起源の爆発が起こり得ない場所がある 震源決定の結果から 人為起源の爆発ではない事象が ある程度ふるい分けられる 1 深い場所 ( 深さ約 2km 以上での爆発は困難 ) 2 海底下 ( 海底下での爆発は技術的に困難 ) 海中や

概論 : 人工の爆発と自然地震の違い ~ 波形の違いを調べる前に ~ 人為起源の爆発が起こり得ない場所がある 震源決定の結果から 人為起源の爆発ではない事象が ある程度ふるい分けられる 1 深い場所 ( 深さ約 2km 以上での爆発は困難 ) 2 海底下 ( 海底下での爆発は技術的に困難 ) 海中や 地震波からみた自然地震と爆発の 識別について 平成 22 年 9 月 9 日 ( 財 ) 日本気象協会 NDC-1 概論 : 人工の爆発と自然地震の違い ~ 波形の違いを調べる前に ~ 人為起源の爆発が起こり得ない場所がある 震源決定の結果から 人為起源の爆発ではない事象が ある程度ふるい分けられる 1 深い場所 ( 深さ約 2km 以上での爆発は困難 ) 2 海底下 ( 海底下での爆発は技術的に困難

More information

地震の将来予測への取組 -地震調査研究の成果を防災に活かすために-

地震の将来予測への取組 -地震調査研究の成果を防災に活かすために- 地震調査研究推進本部は 地震調査研究を一元的に推進する政府の特別の機関です 地震調査研究推進本部は 平成7年1月に 発生した阪神 淡路大震災の教訓 地震調査 基本的な目標 分に伝達 活用される体制になっていなかっ たこと を踏まえ 同年7月 地震防災対策 役 割 特別措置法 に基づき設置された政府の特別 1 総合的かつ基本的な施策の立案 の機関です 行政施策に直結すべき地震調査研究の責任 体制を明らかにし

More information

沈降 隆起 南海トラフで発生した過去の地震には明らかに大きさの違いがある. 図 7 は室戸岬の近くにある室津港での地震毎の隆起 沈降を示したもので, 横軸は西暦である. ハッキリわかっているのは濃い黒で示した宝永, 安政, 昭和の 3 回の隆起量であり, あとは沈降の勾配が同じと仮定して描いてある.

沈降 隆起 南海トラフで発生した過去の地震には明らかに大きさの違いがある. 図 7 は室戸岬の近くにある室津港での地震毎の隆起 沈降を示したもので, 横軸は西暦である. ハッキリわかっているのは濃い黒で示した宝永, 安政, 昭和の 3 回の隆起量であり, あとは沈降の勾配が同じと仮定して描いてある. 地震の発生確率 (Ⅱ) 2012 年 9 月 30 日 産総研梅田康弘 はじめに前回の基礎講座 地震の発生確率 (Ⅰ) (4 月 24 日の ) の続きです. 前回は地震調査研究推進本部が行っている地震発生確率についてご説明し, 発生確率の求め方をお話ししました. 今回は発生確率を求める際の条件などについてご説明し, 確率予測の問題点についてお話しします. 8. 酔っ払いのモデル BPT 分布と更新過程

More information

<4D F736F F D2091E F195AC89CE975C926D D89EF5F97B089A99387>

<4D F736F F D2091E F195AC89CE975C926D D89EF5F97B089A99387> 141 15' 141 18' 141 21' 141 24' 24 51' 24 51' 24 48' IJ T 24 48' IJ M IJ E IJ S 24 45' 24 45' 141 15' 141 18' 141 21' 141 24' 硫黄島の地震活動 (28/8/1~29/1/28) N 24.8 IJM IJT IJE 24.8 24.75 IJS NIED MOD 24.75

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 内陸地震の発生過程について 京都大学防災研究所飯尾能久 内陸地震の発生過程はよくわかって いない. 何が問題なのか? なぜプレート境界から遠くで起こる? なぜ発生間隔が長い? 何が発生間隔を決めるか? 内陸地震とは? 本来は Intraplate earthquake( プレート内部で発生する地震 ) のこと プレート内部とは? プレート境界でないところ プレート境界とは??? Diffused plate

More information

資料2 森田科学官による12月30日提供資料

資料2 森田科学官による12月30日提供資料 地震 火山噴火予知研究の歴史 地震予知研究は昭和 40 年度から, 火山噴火予知研究は昭和 49 年度から, 社会の要請を踏まえて開始された 研究計画の推進のため, 各大学に観測センター, 観測所, 研究部門が設置され, 研究者数が増加し, 研究が進歩し, 地震 火山噴火の知見が得られ, 防災 減災に活かされている 地震予知研究においては,1~3 次の研究体制と観測網の整備,4~6 次の研究体制の充実,7

More information

Seismological Society of Japan - NAIFURU No.99 October, 2014 Report 1 02

Seismological Society of Japan - NAIFURU No.99 October, 2014 Report 1 02 2014.10 No. 99 Contents 2 4 6 8 Seismic Activity in 3 months Seismological Society of Japan - NAIFURU No.99 October, 2014 Report 1 02 03 Seismological Society of Japan - NAIFURU No.99 October, 2014 ひずみ集中帯

More information

8km M km M M8.4 1M M M 東北地方太平洋沖で想定されていた地震 Fig % 8 9% M8. 6 3m M % Fig.1 Distribution of

8km M km M M8.4 1M M M 東北地方太平洋沖で想定されていた地震 Fig % 8 9% M8. 6 3m M % Fig.1 Distribution of 東日本大震災 A Catastrophic Earthquake in Tohoku, Japan 1) 東北地方太平洋沖地震の地震 地震動について 1) Earthquake and Strong Ground Motion of the 211 off the Pacific Coast of Tohoku Earthquake 小林喜久二 Kikuji Kobayashi *1 211 3 11

More information

目次 第 Ⅰ 編本編 第 1 章調査の目的 Ⅰ-1 第 2 章検討体制 Ⅰ-2 第 3 章自然 社会状況 Ⅰ-3 第 4 章想定地震 津波の選定条件等 Ⅰ-26 第 5 章被害想定の実施概要 Ⅰ-37 第 6 章被害想定結果の概要 Ⅰ-48 第 7 章防災 減災効果の評価 Ⅰ-151 第 8 章留意

目次 第 Ⅰ 編本編 第 1 章調査の目的 Ⅰ-1 第 2 章検討体制 Ⅰ-2 第 3 章自然 社会状況 Ⅰ-3 第 4 章想定地震 津波の選定条件等 Ⅰ-26 第 5 章被害想定の実施概要 Ⅰ-37 第 6 章被害想定結果の概要 Ⅰ-48 第 7 章防災 減災効果の評価 Ⅰ-151 第 8 章留意 広島県地震被害想定調査 報告書 平成 25 年 10 月 広島県 目次 第 Ⅰ 編本編 第 1 章調査の目的 Ⅰ-1 第 2 章検討体制 Ⅰ-2 第 3 章自然 社会状況 Ⅰ-3 第 4 章想定地震 津波の選定条件等 Ⅰ-26 第 5 章被害想定の実施概要 Ⅰ-37 第 6 章被害想定結果の概要 Ⅰ-48 第 7 章防災 減災効果の評価 Ⅰ-151 第 8 章留意事項 Ⅰ-158 第 Ⅱ 編結果編

More information

Microsoft PowerPoint - H21生物計算化学2.ppt

Microsoft PowerPoint - H21生物計算化学2.ppt 演算子の行列表現 > L いま 次元ベクトル空間の基底をケットと書くことにする この基底は完全系を成すとすると 空間内の任意のケットベクトルは > > > これより 一度基底を与えてしまえば 任意のベクトルはその基底についての成分で完全に記述することができる これらの成分を列行列の形に書くと M これをベクトル の基底 { >} による行列表現という ところで 行列 A の共役 dont 行列は A

More information

素粒子物理学2 素粒子物理学序論B 2010年度講義第4回

素粒子物理学2 素粒子物理学序論B 2010年度講義第4回 素粒子物理学 素粒子物理学序論B 010年度講義第4回 レプトン数の保存 崩壊モード 寿命(sec) n e ν 890 崩壊比 100% Λ π.6 x 10-10 64% π + µ+ νµ.6 x 10-8 100% π + e+ νe 同上 1. x 10-4 Le +1 for νe, elμ +1 for νμ, μlτ +1 for ντ, τレプトン数はそれぞれの香りで独立に保存

More information

<4D F736F F D D082B882DD90AC89CA95F18D908F F312D385F895E896388CF2E646F63>

<4D F736F F D D082B882DD90AC89CA95F18D908F F312D385F895E896388CF2E646F63> 3.1.8 伸張場におけるひずみ集中メカニズムに関する研究 業務の目的日本海拡大による伸張とその後の圧縮によって形成されたひずみ集中帯の機構を解明するためには 圧縮場のみならず伸張場における応力の境界条件とこれに対する地殻の応答を比較することが重要であるため 現在日本列島の中でも顕著な伸張場を形成している別府島原地溝帯を含む九州地域において 地殻の挙動を正確に把握する検知システム 及び局所的な構造を調べるための観測システムを整備し

More information

自然地理学概説

自然地理学概説 世界と日本の大地形 プレートテクトニクスと世界の大地形 (8.1) 世界の火山と日本の火山 (8.4) 日本列島の成立 日本の山地形成 (8.3) 世界の地震の分布 世界的な火山の分布 世界的な火山の分布を見ると, 太平洋の周りに集中 = 環太平洋火山帯 それ以外の地域も帯状に分布するところがある プレート (p76 図 8.1) 地球の表面はプレートと呼ばれる薄い ( 厚さ約 100~ 150km)

More information

s ss s ss = ε = = s ss s (3) と表される s の要素における s s = κ = κ, =,, (4) jωε jω s は複素比誘電率に相当する物理量であり ここで PML 媒質定数を次のように定義する すなわち κξ をPML 媒質の等価比誘電率 ξ をPML 媒質の

s ss s ss = ε = = s ss s (3) と表される s の要素における s s = κ = κ, =,, (4) jωε jω s は複素比誘電率に相当する物理量であり ここで PML 媒質定数を次のように定義する すなわち κξ をPML 媒質の等価比誘電率 ξ をPML 媒質の FDTD 解析法 (Matlab 版 2 次元 PML) プログラム解説 v2.11 1. 概要 FDTD 解析における吸収境界である完全整合層 (Perfectl Matched Laer, PML) の定式化とプログラミングを2 次元 TE 波について解説する PMLは異方性の損失をもつ仮想的な物質であり 侵入して来る電磁波を逃さず吸収する 通常の物質と接する界面でインピーダンスが整合しており

More information

Microsoft PowerPoint _FPCJ_hirata?v3.pptx

Microsoft PowerPoint _FPCJ_hirata?v3.pptx FPCJ プレス ブリーフィング 日本の地震 火山活動の現状と今後 東京大学地震研究所 地震予知研究センター長 教授平田直 日時 : 2017 年 2 月 15 日 ( 木 )12:00~13:30 会場 : 記者会見室 ( 千代田区内幸町 2-2-1 日本プレスセンタービル 6 2017/2/15 FPCJ プレス ブリーフィング 1 内容 1. 2016 年熊本地震 2. 首都直下地震 3. 地震と火山活動

More information

令和元年6月 地震・火山月報(防災編)

令和元年6月 地震・火山月報(防災編) (2) 地震活動ア. 地震の発生場所の詳細及び地震の発生状況 2019 年 6 月 18 日 22 時 22 分に山形県沖の深さ 14km で M6.7 の地震 ( 最大震度 6 強 ) が発生した この地震発生以降 北東 - 南西方向に長さ約 20km の領域で 地震活動が本震 - 余震型で推移している 最大規模の余震は 6 月 19 日 00 時 57 分に発生した M4.2 の地震 ( 最大震度

More information

佐賀県の地震活動概況 (2018 年 12 月 ) ( 1 / 10) 平成 31 年 1 月 15 日佐賀地方気象台 12 月の地震活動概況 12 月に佐賀県内で震度 1 以上を観測した地震は1 回でした (11 月はなし ) 福岡県 佐賀県 長崎県 熊本県 図 1 震央分布図 (2018 年 1

佐賀県の地震活動概況 (2018 年 12 月 ) ( 1 / 10) 平成 31 年 1 月 15 日佐賀地方気象台 12 月の地震活動概況 12 月に佐賀県内で震度 1 以上を観測した地震は1 回でした (11 月はなし ) 福岡県 佐賀県 長崎県 熊本県 図 1 震央分布図 (2018 年 1 佐賀県の地震活動概況 (2018 年 12 月 ) ( 1 / 10) 平成 31 年 1 月 15 日佐賀地方気象台 12 月の地震活動概況 12 月に佐賀県内で震度 1 以上を観測した地震は1 回でした (11 月はなし ) 福岡県 佐賀県 長崎県 熊本県 図 1 震央分布図 (2018 年 12 月 1 日 ~31 日 深さ 30km 以浅 ) 灰色の線は地震調査研究推進本部の長期評価による活断層を示しています

More information

<4D F736F F F696E74202D208CB48E7197CD8A7789EF8AE989E6835A C478CB3816A2E B8CDD8AB B83685D>

<4D F736F F F696E74202D208CB48E7197CD8A7789EF8AE989E6835A C478CB3816A2E B8CDD8AB B83685D> SSHAC 手法に基づく原子力発電所サイトの確率論的地震ハザード評価 (BC0504) 内陸地殻内地震を対象とした確率論的地震動ハザード評価課題とその解決に向けて 隈元崇 @ 岡山大学 Lettis, 2013 http://www.pge.com/includes/docs/pdfs/shared/edusafety/ systemworks/dcpp/sshac/sugmworkshops/ 特性化震源モデル

More information

PowerPoint Presentation

PowerPoint Presentation 資料 11 第 191 回 地震予知連絡会資料 2011 年 6 月 13 日 統計数理研究所 確率的除群法 (Stochastic declustering) Original data 1926 2011 Mar. M 5.0, h 100km De-clustered data 1 Latitude Time Time Time Latitude Time Longitude Longitude

More information

様々なミクロ計量モデル†

様々なミクロ計量モデル† 担当 : 長倉大輔 ( ながくらだいすけ ) この資料は私の講義において使用するために作成した資料です WEB ページ上で公開しており 自由に参照して頂いて構いません ただし 内容について 一応検証してありますが もし間違いがあった場合でもそれによって生じるいかなる損害 不利益について責任を負いかねますのでご了承ください 間違いは発見次第 継続的に直していますが まだ存在する可能性があります 1 カウントデータモデル

More information

4 段階推定法とは 予測に使うモデルの紹介 4 段階推定法の課題 2

4 段階推定法とは 予測に使うモデルの紹介 4 段階推定法の課題 2 4 段階推定法 羽藤研 4 芝原貴史 1 4 段階推定法とは 予測に使うモデルの紹介 4 段階推定法の課題 2 4 段階推定法とは 交通需要予測の実用的な予測手法 1950 年代のアメリカで開発 シカゴで高速道路の需要予測に利用 日本では 1967 年の広島都市圏での適用が初 その後 1968 年の東京都市圏など 人口 30 万人以上の 56 都市圏に適用 3 ゾーニング ゾーニングとネットワークゾーン間のトリップはゾーン内の中心点

More information

構造力学Ⅰ第12回

構造力学Ⅰ第12回 第 回材の座屈 (0 章 ) p.5~ ( 復習 ) モールの定理 ( 手順 ) 座屈とは 荷重により梁に生じた曲げモーメントをで除して仮想荷重と考える 座屈荷重 偏心荷重 ( 曲げと軸力 ) 断面の核 この仮想荷重に対するある点でのせん断力 たわみ角に相当する曲げモーメント たわみに相当する ( 例 ) 単純梁の支点のたわみ角 : は 図 を仮想荷重と考えたときの 点の支点反力 B は 図 を仮想荷重と考えたときのB

More information

線積分.indd

線積分.indd 線積分 線積分 ( n, n, n ) (ξ n, η n, ζ n ) ( n-, n-, n- ) (ξ k, η k, ζ k ) ( k, k, k ) ( k-, k-, k- ) 物体に力 を作用させて位置ベクトル A の点 A から位置ベクトル の点 まで曲線 に沿って物体を移動させたときの仕事 W は 次式で計算された A, A, W : d 6 d+ d+ d@,,, d+ d+

More information

日向灘 佐伯市で震度 2 を観測 8 日 08 時 33 分に日向灘で発生した M3.9 の地震 ( 深さ 31km) により 佐伯市 愛媛県西予市 高知県宿毛市などで震度 2 を観測したほか 大分県 宮崎県 愛媛県および高知県で震度 1 を観測しました ( 図 1) 今回の地震の震源付近 ( 図

日向灘 佐伯市で震度 2 を観測 8 日 08 時 33 分に日向灘で発生した M3.9 の地震 ( 深さ 31km) により 佐伯市 愛媛県西予市 高知県宿毛市などで震度 2 を観測したほか 大分県 宮崎県 愛媛県および高知県で震度 1 を観測しました ( 図 1) 今回の地震の震源付近 ( 図 (1/10) 大分県の地震活動概況 (2018 年 12 月 ) 平成 31 年 1 月 15 日大分地方気象台 地震活動概要 (2018 年 12 月 ) 12 月に大分県内で震度 1 以上を観測した地震は 1 回 ( 下図範囲外 1 回 ) でした (11 月 :4 回 ) 震央分布図内の地震の総数 :592 マク ニチュート 大分県 深さ 震央分布図 (2018 年 12 月 1 日 ~31

More information

切片 ( 定数項 ) ダミー 以下の単回帰モデルを考えよう これは賃金と就業年数の関係を分析している : ( 賃金関数 ) ここで Y i = α + β X i + u i, i =1,, n, u i ~ i.i.d. N(0, σ 2 ) Y i : 賃金の対数値, X i : 就業年数. (

切片 ( 定数項 ) ダミー 以下の単回帰モデルを考えよう これは賃金と就業年数の関係を分析している : ( 賃金関数 ) ここで Y i = α + β X i + u i, i =1,, n, u i ~ i.i.d. N(0, σ 2 ) Y i : 賃金の対数値, X i : 就業年数. ( 統計学ダミー変数による分析 担当 : 長倉大輔 ( ながくらだいすけ ) 1 切片 ( 定数項 ) ダミー 以下の単回帰モデルを考えよう これは賃金と就業年数の関係を分析している : ( 賃金関数 ) ここで Y i = α + β X i + u i, i =1,, n, u i ~ i.i.d. N(0, σ 2 ) Y i : 賃金の対数値, X i : 就業年数. ( 実際は賃金を就業年数だけで説明するのは現実的はない

More information

H19年度

H19年度 Ⅲ-6. 津波高及び津波浸水の分布 1. 元禄型関東地震の津波数値シミュレーション ( 東京湾 ) 1.1 津波数値シミュレーションの計算条件 津波の影響も考慮すべき地震として採用した元禄型関東地震行谷ほか (2011) モデルについて 以下の計算条件で津波遡上の数値シミュレーションを実施した 使用した断層モデル ( 波源モデル ) 元禄型関東地震行谷ほか (2011) モデル計算条件 メッシュサイズ

More information

indd

indd 45 Crustal Deformation and a Fault Model of the West - off Fukuoka Prefecture Earthquake in 2005 測地観測センター小清水寛 畑中雄樹 根本盛行 Geodetic Observation Center Hiroshi KOSHIMIZU, Yuki HATANAKA and Moriyuki NEMOTO

More information

Information Theory

Information Theory 前回の復習 情報をコンパクトに表現するための符号化方式を考える 情報源符号化における基礎的な性質 一意復号可能性 瞬時復号可能性 クラフトの不等式 2 l 1 + + 2 l M 1 ハフマン符号の構成法 (2 元符号の場合 ) D. Huffman 1 前回の練習問題 : ハフマン符号 符号木を再帰的に構成し, 符号を作る A B C D E F 確率 0.3 0.2 0.2 0.1 0.1 0.1

More information

Microsoft PowerPoint - aep_1.ppt [互換モード]

Microsoft PowerPoint - aep_1.ppt [互換モード] 物理計測法特論 No.1 第 1 章 : 信号と雑音 本講義の主題 雑音の性質を理解することで 信号と雑音の大きさが非常に近い状態での信号の測定技術 : 微小信号計測 について学ぶ 講義の Web http://www.g-munu.t.u-tokyo.ac.jp/mio/note/sig_mes/tokuron.html 物理学の基本は実験事実の積み重ねである そして それは何かを測定することから始まる

More information

Microsoft PowerPoint - シミュレーション工学-2010-第1回.ppt

Microsoft PowerPoint - シミュレーション工学-2010-第1回.ppt シミュレーション工学 ( 後半 ) 東京大学人工物工学研究センター 鈴木克幸 CA( Compter Aded geerg ) r. Jaso Lemo (SC, 98) 設計者が解析ツールを使いこなすことにより 設計の評価 設計の質の向上を図る geerg の本質の 計算機による支援 (CA CAM などより広い名前 ) 様々な汎用ソフトの登場 工業製品の設計に不可欠のツール 構造解析 流体解析

More information

(Microsoft PowerPoint - \221\34613\211\361)

(Microsoft PowerPoint - \221\34613\211\361) 計算力学 ~ 第 回弾性問題の有限要素解析 (Ⅱ)~ 修士 年後期 ( 選択科目 ) 担当 : 岩佐貴史 講義の概要 全 5 講義. 計算力学概論, ガイダンス. 自然現象の数理モデル化. 行列 場とその演算. 数値計算法 (Ⅰ) 5. 数値計算法 (Ⅱ) 6. 初期値 境界値問題 (Ⅰ) 7. 初期値 境界値問題 (Ⅱ) 8. マトリックス変位法による構造解析 9. トラス構造の有限要素解析. 重み付き残差法と古典的近似解法.

More information

Microsoft Word doc

Microsoft Word doc . 正規線形モデルのベイズ推定翠川 大竹距離減衰式 (PGA(Midorikawa, S., and Ohtake, Y. (, Attenuation relationships of peak ground acceleration and velocity considering attenuation characteristics for shallow and deeper earthquakes,

More information

鉱物と類似の構造を持つ白雲母の鉱物表面に挟まれた塩化ナトリウム (NaCl) 水溶液が 厚さ 1 ナノメートル ( 水分子約 3 個分の厚み ) 以下まで圧縮されても著しい潤滑性を示すことを実験的に明らかにしてきました しかし そのメカニズムについては解明されておらず 世界的にも存在が珍しいクリープ

鉱物と類似の構造を持つ白雲母の鉱物表面に挟まれた塩化ナトリウム (NaCl) 水溶液が 厚さ 1 ナノメートル ( 水分子約 3 個分の厚み ) 以下まで圧縮されても著しい潤滑性を示すことを実験的に明らかにしてきました しかし そのメカニズムについては解明されておらず 世界的にも存在が珍しいクリープ 報道関係者各位 地震を起こさない断層すべりのメカニズムを世界で初めて解明 平成 23 年 8 月 5 日 東京工業大学お茶の水女子大学高エネルギー加速器研究機構岡山大学東京工業大学大学院理工学研究科の佐久間博特任助教 お茶の水女子大学大学院人間文化創成科学研究科の近藤敏啓教授 高エネルギー加速器研究機構 (KEK) 物質構造科学研究所の中尾裕則准教授 東京工業大学大学院理工学研究科の河村雄行教授 (

More information