こ う えい フ ォ ー ラ ム第 24 号 / 2016.3

ミャンマー国エーヤワディ地域における地震 ・ 津波のリスク評価
EARTHQUAKES AND TSUNAMI RISK ESTIMATION IN AYEYAWADY DELTA OF MYANMAR
櫻庭 雅明 * ・ 佐藤 誠一 ** ・ 石見 和久 *** ・ 木村 健太郎 ***Masaaki SAKURABA, Seiichi SATO, Kazuhisa IWAMI and Kentaro KIMURA

Myanmar is at a high risk from natural disasters such as ﬂoods, cyclones, earthquakes andtsunami. Thus, it is necessary to estimate the damage risks of natural disaster in advance. Thispaper presents earthquake and tsunami risk estimation in Ayeyawady Delta of Myanmar.
In this study, recent epicenter and magnitude of earthquake data was collected. Analysis ofearthquake frequency in Myanmar was carried out as well as tsunami damage estimation bynumerical simulation. The results present the frequency and the tsunami scale in the typicalarea of Ayeyawady Delta in Myanmar.

Keywords ： Ayeyawady Delta, Earthquake, Tsunami, Numerical Simulation

ルギーを保存しており、 今後巨大地震が発生する危険性が1． はじめにある。 さらに、 ミャンマーのデルタ地域 （エーヤワディ管区）

ミャンマーは、 人口 6,242 万人、 国土は日本の 1.8 倍で主要産業が農業の国である。 ミャンマーは 2011 年 3 月の軍政は標高の低い地域が広がり津波の来襲による甚大な被害が懸念される。

から民政への移行で世界の注目を集め、 その経済発展の潜本検討では、 ミャンマーのエーヤワディ管区周辺における在能力の高さから今後大規模かつ急速な開発が見込まれてい地震および津波の被害想定を行ったものである。 地震の検る。 一方、 ミャンマーはその地形的特性から、 国内に活断層

討は、過去に発生した地震の震源および規模 （マグニチューをいくつもかかえ、これまで多くの地震による被害が生じている。

ド）から沿岸部で発生する地震の発生頻度と確率を算定した。

また、 さらにエーヤワディ管区の沿岸にあるアンダマン海はス

また、 津波に寄与する断層を設定して、 発生頻度の高い津マトラ沖から南北に断層が貫いている。 図－ 1 に検討対象位波 （レベル 1 津波） および最大クラスの津波 （レベル 2 津置図を示す。 2004 年 12 月に来襲したスマトラ沖地震によるイ波）を想定した。これらの津波条件に対して津波シミュレーショ

ンド洋津波により被害が生じた。 さらに、 2008 年のサイクロンンを実施して、 エーヤワディ管区全域およびモデル地域であNargis では強風および高波 ・ 高潮により死者行方不明者約る Labutta の浸水想定を行った。 この結果からレベル 1、 214 万人という甚大な被害が生じたことが知られている 。

津波に対する浸水被害規模について考察を行った。

1）

ミャンマーは、 被害規模の記録が残っている最近 60 年程
2）度では、 洪水およびサイクロンの頻度が高い 。 また、 頻度

なお、 本研究は JICA プロジェクト 「ヤンゴン港 ・ 内陸水運施設改修プロジェクト」 の成果の一部をまとめたものである。

は少ないが発生の際には甚大な被害が予想されるのは地震および津波である。 ミャンマーに発生した地震は、 1912 年にシャン州で発生したマグニチュード 8 が記録上最大であるが、1839 年、 1930 年にもマンダレー管区とヤンゴン管区で大地震があり、 それぞれ被害が生じた。 また、 2004 年のスマトラ沖地震によるインド洋津波ではエーヤワディ管区で死者 61 人、負傷者 42 人、 被災者 2,592 人、 倒壊家屋 601 棟の被害が生じている 2）。
ミャンマーでは、 前述のとおり、 沿岸部に断層があり、 西側を抜けているが、 近傍で地震は発生しておらず地震エネ*
技術本部 中央研究所 総合技術開発部** コンサルタント国内事業本部 流域 ・ 都市事業部 地盤環境部*** コンサルタント海外事業本部 アジア推進事業部 アジア事業推進室

図－ 1 検討対象位置図
出典 ： OpenStreetMap
21

ミャンマー国エーヤワディ地域における地震 ・ 津波のリスク評価

2． 現地調査（1） 津波被害の状況
2004 年 12 月 26 日に発生したスマトラ沖地震による津波はの浸水が見られているのに対して、 津波による被害は殆ど見られていない。 Ywe 川および Labutta はエーヤワディのやや西側の地域であり、 当時の津波の規模が比較的小さかったためであると考えられる。

ミャンマーのデルタ地域まで来襲しており、 一部で浸水被害が生じている。 現地の浸水深や浸水範囲について詳細に報告された事例は少ないが、 エーヤワディ東部の一部集落 （AingHlaing village） で浸水が生じた事例がある。 図－ 2 にその浸水被害の写真を示す。 家屋の一部が損壊しており、 1m 程度の浸水が生じていることが分かる。

3． エーヤワディ地域における地震発生頻度の評価（1） ミャンマー国周辺の地震の発生状況
次にミャンマー全域における地震の発生状況を検討した。
図－ 4 にミャンマー近傍における地震の発生状況 （震源分布）を示す。 これまでミャンマーでは多くの地震が発生しており、（2） 現地踏査

1900 年からマグニチュード 4.8 以上の地震が合計で 760 回

地震 ・ 津波の検討に先立ってエーヤワディの西側にある記録されている。 この図では、 スマトラ沖周辺で数多くの記録

Ywe 川およびその上流の集落である Labutta の現地踏査をが見られているが、 ミャンマー国内でも高い頻度で発生してい行った。 なお、 踏査の際に現地ヒアリングを行って、 被害状況ることが分かる。

を推定するための基礎資料とした。 なお、 被害については津波被害のみでなく、 Nargis 来襲時の高潮による浸水状況も合（2） 地震の確率的評価わせて行っている。 ヒアリング実施の状況は図－ 3 に示すとおUSGS の 地 震 デ ー タ ベ ー ス よ り、 イ ン ド － オ ー ス ト ラ リ アりである。 ほぼすべての地域でサイクロン Nargis で 1m 前後

プレートで発生した過去の地震記録 （1973 ～ 2012 年） を抽出した。 このデータベースから抽出した地震の震源分布を図－ 5 に示す。 また、 マグニチュード 5.0 以上の発生状況を示したもの （1973 年以降を抜粋） が図－ 6 である。 この結果から分かるように、 2004 年のスマトラ沖地震以降にマグニチュード 8.0 以上の地震が多数発生していることが分かる。
本検討では、 ミャンマーのエーヤワディ沿岸における津波の来襲 ・ 浸水規模について検討することを目的としている。 地震および津波の規模は、 日本の津波レベルの考え方と同様に位置づけるものとした。 具体的には以下の 2 種類とした。

図－ 2 インド洋津波によるデルタ地域の浸水被害

図－ 3 現地踏査 ・ ヒアリング調査の状況
22

図－ 4 ミャンマー周辺の既往地震の震源分布こ う えい フ ォ ー ラ ム第 24 号 / 2016.3●●

発生頻度の高い津波 （レベル 1 津波、 以下 L1 津波）

リヒター式に当てはめる。 また、 図－ 8 に示すような断層を L2に対しては、 可能な限り構造物で人命 ・ 財産を守り切津波の条件として想定した場合に、 断層規模とモーメントマグる 「防災」 を目指すもの。

ニチュードの関係から Mw を求めると 9.4 と推定できる。 なお、発生頻度は低いが沿岸域への影響が甚大な最大クラス

L1 は後述のとおり Mw=9.0 としている。 これより当該プレートの津波 （レベル 2 津波、 以下 L2 津波） に対しては、境界部において L1 津波、 L2 津波を生じさせる規模の地震発

最低限人命を守り、 被害をできるだけ小さくする 「減災」

生確率をグーテンベルグ ・ リヒター式から算定すると以下のとを目指す。

おりとなる。

なお、 日本では L1 津波は設計津波、 L2 は想定最大クラ1.0E+03スの位置づけになるが、 ミャンマーにおいて日本と同等の対策2.5
近似直線
地震数施設の整備は現実的でない。 ここでいう L1、 L2 津波はそのN

1.0E+021.5
b値

地震数た。 地震の発生頻度の設定は、USGS の地震データを用いて、地震規模と発生確率の関係を求めた。 地震規模と発生確率の関係を求めるにあたっては、 グーテンベルグ ・ リヒター式を用2.0

b値

発生頻度に対する津波浸水想定の把握を目的として位置づけ

1.0
1.0E+01

いた。 グーテンベルグ ・ リヒター式は、 データの対象期間にお0.5ける地震規模別の地震発生回数から年発生回数を求め、 この逆数をとることにより、 地震規模ごとの発生確率を統計的見地から算出するものである。 図－ 7 に地震規模 （マグニチュード）と地震発生回数の関係を示す。 この結果をグーテンベルグ ・

1.0E+00

0.0
5.0

6.0

7.0

8.0

9.0

10.0

マグニチュード：M

図－ 7 地震規模と地震発生回数の関係

図－ 5 震源分布 （USGS のデータベースから記載）
10.0

8.0
7.0

6.0
5.0
4.0

西暦

図－ 6 地震発生状況の時刻歴 （M>5.0）

2011年

2009年

2007年

2005年

2003年

2001年

1999年

1997年

1995年

1993年

1991年

1989年

1987年

1985年

1983年

1981年

1979年

1977年

1975年

3.0
1973年

地震規模（マグニチュード）

9.0

図－ 8 ミャンマー沿岸を想定した断層モデル（緑色：2004 年
スマトラ沖地震、 赤色 ： 本検討で追加）
出典 ： USGS
23

ミャンマー国エーヤワディ地域における地震 ・ 津波のリスク評価●

L1 津波（Mw=9.0）: 年間想定発生回数 :0.0217 回（約45 年に一度発生する確率）●

L2 津波（Mw=9.4）: 年間想定発生回数 :0.0113 回（約90 年に一度発生する確率）

4． 津波リスクの検討（1） 計算条件
前述の地震 ・ 断層モデルの検討結果から、 L1 津波および上記の結果より、 L1 と L2 の発生頻度が 2 倍程度と小さい。

L2 津波の条件を設定した。 断層パラメータは表－ 1 に示すと日本周辺におけるグーテンベルグ ・ リヒター式の算定結果によおりである。 なお、 この断層パラメータは、 2004 年インド洋津れば M9.4 クラスで 300 年に一回程度と推定できたが、 イン波における断層パラメータを東北大 4） が設定したものであり、ド－オーストラリアプレートを対象とした場合、 発生確率が非常

L2 津波の条件はそれを拡張したものである。 津波の来襲 ・ 浸に高くなる。 これは 2004 年スマトラ沖地震以降、 当該プレー水状況の検証 ・ 予測にあたっては、 津波シミュレーションを採トの活動が非常に活発であることによる。 日本で想定される L1用した。 津波シミュレーションのモデルおよび計算方法は、 「津および L2 の発生頻度とは違いが生じているが、 震源の位置

波浸水想定設定の手引き」 5） に示されたものと同等のものをにより津波の規模が異なり、 次に示す断層モデルにより最終的用いているため、 説明は割愛する。 計算条件は、 表－ 1 に示な条件を設定した。

すとおりであり、 計算格子は対象とする計算範囲の格子間隔を徐々に細かくするネスティング手法を採用した。 なお、 今回の（3） 断層モデルの検討

検 討 で は、 Ywe 川 上 流 30km 程 度 に 位 置 す る Labutta に津波の来襲 ・ 浸水予測を検討するにあたり、 断層モデルを着目している。 計算範囲は図－ 9 に示すとおりであり、 広域検討した。 断層モデルの基本として、 2004 年スマトラ地震津はベンガル湾およびアンダマン海を包括できるようにしている。

波の断層モデルを用いた。 なお、 ミャンマーでの津波規模が最小メッシュ幅は 10m であり、 Labutta 周辺を包括している。

大きくなることを想定して、 北部空白域も連動することを仮定

津波シミュレーションに必要となるメッシュごとの粗度係数は、 小した （図－ 8 緑色のセグメント参照）。 北部空白域では、 スマ

谷らの手法を基本とするが、 それに必要となる土地利用は 「ミャトラ地震の余震はほとんど発生していないが （図－ 8 の黄丸

ンマー国 国家復興開発計画のための地理情報データベースが余震震源）、 地震そのものが発生しない区域ではないため、将来的に連動した地震が発生する可能性がある。 なお、 当該地域では過去約 3000 年に 4 回の活動が 900 年間隔で発生していることが地形や堆積物年代測定調査で明らかとされている （西海岸 ManAung 島）。 既往最大津波に加え、 科学的表－ 1 津波シミュレーションの計算条件
項

目

ブダクション構造も津波断層として考慮し、 断層モデルを設定した。 なお、 宍倉ら 3） によれば、 ミャンマーの西方沖には、2004 年スマトラ－アンダマン地震の震源域から連続する沈み込みプレート境界が横たわっている。 ミャンマー西方沖の沈み込み境界でインドプレートがビルマプレートに対して年間約 101800m ⇒ 600m ⇒ 200m ⇒50m ⇒ 10m の格子を接続して、水位 ・流量を境界条件として与えた。

計算時間間隔

1.0 秒

計算時間

地震発生から 16 時間

計算格子の測地系

WGS 1984

計算格子の座標系

UTM 54N

基本方程式と解法

非線形長波理論式
Leap-Frog 差分法

沖側境界と接続境界

沖側境界は自由透過とし、 接続境界は水位と流量を接続

遡上境界条件

陸域では小谷ほか (1998) の遡上境界条件とし、 他は完全反射。

地震断層モデル

2004
IndianOceanEarthquake （東北大モデル）
北部空白域も上記断層パラメータを参考に設定

潮位条件

1.44m （L2 津波、 満潮位）

地盤変位条件

沈降は考慮し、 隆起は考慮しない

底面摩擦モデルManning の粗度係数を土地利用に応じて設定

構造物モデル

200m、 50m、 10m 領域において、YWE 川右岸のみ設定

mm ／年の速度で、 北北東に沈み込みつつあると推定され、数百年間隔で海溝型プレート間地震が発生していると考えられているとされている。 地震の発生間隔は数百年と考えられているが、 前述のとおり、 スマトラ沖地震以降の地震活動の活発化により発生頻度は短くなる。 また、 2004 年スマトラ沖地震津波の断層とインド洋津波の断層 + ミャンマー西部サブダクション構造を含めた断層のそれぞれモーメントマグニチュードは 9.0
および 9.4 と推定されるが、 対象とする地域の到達する津波の規模は断層の存在する位置に依存する。 このことから、 2004
年スマトラ沖地震津波 （当該地域における既往最大津波） をL1 津波、 インド洋津波の断層 + ミャンマー西部サブダクション構造を含めた断層を L2 津波にすることで概ね津波規模を説明できると仮定して採用している。

24

容計算格子間隔知見を持って考えうる可能性を考慮した最大クラスの津波 （L2津波） を設定する必要があることから、 ミャンマー国西部のサ内こ う えい フ ォ ー ラ ム第 24 号 / 2016.3

構築調査 （2004 年 8 月）」 の報告に従った。 また、 この報告

タは Ywe 川周辺以外は 200m メッシュとなっているため、 解で得られない範囲の土地利用は、 衛星画像 （Landsat-7） で像度はそれ以下として評価されないため、 上記方法は概略で補足した。 図－ 10 に土地利用の分布を示す。 また、 Ywe 川はあるが、 水深データの得られない地域としては十分であるとの堤防高を設定にあたっては、 現地調査の結果を反映した。

判断した。 最終的に作成したメッシュ （河道内標高分布） を図－ 12 に示す。

（2） 河川の地形モデルの設定
エーヤワディのデルタ地帯は、 図－ 10 の土地利用から判（3） 2004 年インド洋津波による検証 （L1 津波の検討）

断できるように、 大小の河川が非常に多くあることが分かる。

津波シミュレーションの妥当性を検証するために、 2004 年津波の来襲の状況を推定するにあたっては、 これらの河川をインド洋津波の再現計算を行った。 エーヤワディ全域の最大考慮する必要がある。 しかしながら、 ほぼ全ての河川は縦横津波高さの分布を図－ 13 に示す。 エーヤワディデルタ周辺は断測量の成果が存在しないことから、 地形データとして作成す複雑な海岸線を有しており、 水深は非常に浅い地域となってるためには工夫が必要となる。 本検討では、 すべての河川のいる。 この地形特性の影響により、 最大津波高は各地で局所測量を行っていないため、 YWE 川で行われた縦横断測量、的に異なるが、 概ね 3m 程度となる。 なお、 検討対象とした

Pathein 川および Pya Ma Law 川で行われた縦断測量の成Labutta では 2m 程度の高さとなった。 なお、 2004 年のイン果を用い、 河道の類型設定を行い、 その他の河川に対して川ド洋津波は、 当該地域では観測記録・痕跡記録は存在しない。

幅から設定できるようにした。 具体的には、 図－ 11 に示すよ計算の整合性を確認するためには、 現地でのヒアリング結果にうな河道に対して左右岸、 中央に流路があるケースに分類し、基づき、 各地点の河川水位が堤防や船着き場の天端近くまで川幅から水深を設定した。 なお、 本検討に用いたメッシュデー上がったという証言に基づき比較した。 図－ 14 に水位観測地点 （目視観測の BM があるのみ） の標高と本計算の最大水位の比較を示す。当該地区は浸水しない結果であり、シミュレーションも同等の結果となった。

図－ 9 計算範囲

図－ 12 河川の水深を考慮した水深分布図

図－ 10 土地利用状況(ⅲ)(ⅳ)-40-200x0x

20204040(ⅲ)(ⅳ)-0.5

y-1.0-1.5-20

y-1.5 -1.0 -0.5

y-1.0-1.5-3-4-4-40(ⅰ)(ⅱ)

0.00.0-0.5-1-2y-2-3yy-1.5 -1.0 -0.5

-10.0(ⅰ)(ⅱ)

(ⅲ)

00(ⅱ)

0.0

(ⅰ)

-40-40-20-200x0x20204040-40-40-20-20

図－ 11 河道横断の類型区分の例

0x0x20204040

図－ 13 L1 津波の最大水位分布25

ミャンマー国エーヤワディ地域における地震 ・ 津波のリスク評価

水位（m）

87

水位観測地点のBM標高

6最大水位

54
3
210
㻜

㻡㻜㻜㻜

㻝㻜㻜㻜㻜

㻝㻡㻜㻜㻜
㻞㻜㻜㻜㻜
㻞㻡㻜㻜㻜
河口からの距離（m）

㻟㻜㻜㻜㻜

㻟㻡㻜㻜㻜

㻠㻜㻜㻜㻜

図－ 14 Ywe 川における津波縦断水位（4） 最大クラスの津波における想定 （L2 津波の検討）

図－ 15 L2 津波における最大水位分布

次に、 前述の断層モデルの検討で設定した L2 津波について検討した。断層モデルは図－ 8 に示したとおり、ミャンマー6.0の西岸を抜ける条件となっており、 2004 年インド洋津波に比水位（m）

4.0

べてミャンマー全域に厳しい条件となっている。 この条件で

2.0
時間

0.0

エーヤワディデルタ周辺および Labutta 沿岸における津波最0123456789101112

-2.0

大水位分布および津波の時間変化を図－ 15 ～図－ 16 に示

6.0す。 Ywe 川河口での最大津波高は 4.3m、 Labutta で 3.0m

Labutta

4.0
水位（m）の津波が来ることが予測される。 また、 Pathein 川では約 30分で 6m 近い津波が来襲することが予想され、 エーヤワディ

2.0
時間

0.00123456789101112-2.0

西側で津波が高くなることがこの結果より分かる。
（5） 津波アニメーションによる防災啓発

Pathein River

図－ 16 各 地点における津波高時刻歴（上図：Pathein 川河口、下図：Labutta 周辺）

本検討で実施した津波シミュレーションをできるだけ関係者にわかりやすく見せることを目的として、 津波動画 （アニメーシ ョン） を作成している。 津波動画の作成にあたっては、 レンダリングソフトの Pov-Ray 6） を用いて、 津波計算結果のコマ撮りからアニメーションを作成している。 また、 現地の状況
は オ ル ソ 画 像 を 重 ね て 表 現 し て い る。 広 域 の 計 算 か ら 津 波の伝播状況が分かるように広域、 エーヤワディデルタ全域、Labutta 周辺のそれぞれにおける鳥瞰図をベースとしてアニメーションを作成した。 アニ メ ー シ ョ ン の ス ナ ッ プ シ ョ ッ トを図－ 17 津 波アニメーションのスナップショット （左図：第1波Labutta周辺来襲時、右図：第2波エーヤワディ来襲時）

図－ 17 に示す。 この結果より予想される L2 津波の条件ではエーヤワディデルタの広い範囲で浸水が生じることが見て取

参考文献

れ、 Labutta 周辺でも広い範囲で浸水することが分かる。

1）

DMH、 DRR et al.: Hazard Proﬁle of Myanmar, 20092）

田平由希子、 川崎昭如、 市原裕之 : 民政移行後のミャンマー中
央政府の防災体制と今後の課題、 地域安全学会論文集、 No.21、5． おわりに本研究では、 既存資料からミャンマーにおける地震の発生

pp.189-198、 20133）

報告、 活断層研究センターニュース、 No.76、 pp.1-3、 2007

頻度および確率を算定した。 また、 現地ヒアリング ・ 測量の結果を用いて津波リスクに関する検討を数値シミュレーションを用宍倉正展 ・ 岡村行信 ・ 藤野滋弘 : ミャンマー西海岸古地震調査

4）

木村健太郎 ・ 石見和久 ・ 佐藤誠一 ・ 福田忠弘 ・ 三上勉 ・ 菊池茂 :いて行った。 その結果、 最大クラスとして想定した津波に対しスマトラ沖地震津波被害実態調査に基づく波源モデルの検証と建物てはエーヤワディデルタのほぼ全域で浸水し、 場所によっては被害関数の提案、 海岸工学論文集、 Vol.53、 pp.301-305、 2006

地震発生後 1 時間以内で津波が来襲することが確認された。

5）

政策総合研究所河川研究部海岸研究室 : 津波浸水想定設定の手

今後、 ミャンマーにおける津波規模の定義、 浸水想定および引き ver2.0、 2012

防災啓発のための基礎資料およびツールとして本成果は活用できるものと考えている。
26

国土交通省水管理 ・ 国土保全局海岸室 ・ 国土交通省国土技術

6）

Pov-Ray ホームページ :http://www.povray.org/