SSS25-01

会場:A04

時間:5 月 25 日 09:00-09:15

M6.5 中国雲南省昭通地震における建物被害と強震動特性の関係―その 1：現地調査
報告および建物被害分析
Relationships between building damage and characteristics of strong ground motions during the M6.5 Zhaotong Earthquake-I
王 欣 1∗ ; 馬 強 2 ; 司 宏俊 3 ; 党 紀 4 ; 呉 浩 5 ; 倉橋 奨 6
WANG, Xin1∗ ; MA, Qiang2 ; SI, Hongjun3 ; DANG, Ji4 ; WU, Hao5 ; KURAHASHI, Susumu6
1

東北大学災害科学国際研究所, 2 中国地震局工程力学研究所, 3 東京大学地震研究所, 4 埼玉大学, 5 愛知工業大学地域防災
研究センター, 6 愛知工業大学
1
IRIDeS, Tohoku University, 2 Institute of Engineering Mechanics, China Earthquake Administration, 3 Earthquake Research Institute, The University of Tokyo, 4 Saitama University, 5 DPREC, Aichi Institute of Technology, 6 Aichi Institute of Technology

2014 年 8 月 3 日 16 時 30 分（北京時間）に中国雲南省昭通市魯甸 (ルーデン) 県龍頭山鎮付近（27.110N、103.330E）でマグニチュード 6.5 の地震が発生した。震央付近で最大震度 IX（JMA 震度の 5 強と相当）で、最大加速度は震央から約8km を離れた龍頭山鎮で 949gal で記録された。死者数が 617 人（2014 年 8 月 8 日 15 時現在）で、このうち約 85 ％（526人）は龍頭山鎮に集中した。建物倒壊による人的被害が圧倒的に多く、建物の耐震性能不足と強震動特性の両面から原因が考えられる。
被害現場の建物被害状況および地盤の増幅と周期特性を究明するため、本研究では 2014 年 11 月 9 日〜12 日に現地調査を実施した。現地調査の内容として、断層近傍の強震動観測点および昭通市内において、建物被害の調査および地盤常時微動観測を行っている。現地調査を実施した強震動観測点を図 1 に示す。本震時に図 1 に示す強震動観測点の記録の最大加速度（PGA）が昭通市 15gal、茨院 45gal、龍頭山 949gal、鉛場 137gal、馬樹 135gal である。PGA の距離減衰式によって、断層最短距離の増大と共に PGA が急速に減少することが分かった。
被災地ではレンガ造建物がよく使用され、龍頭山鎮では建物の倒壊が多く見られる。王らが提案した 2008 年四川地
震の建物の被害関数に基づいて予測した龍頭山鎮の倒壊率は 79 ％で、この予測値と建物被害現地調査の結果がよく一致する。鉛場と馬樹にレンガ造建物の壁にひび割れが散見できるが、建物倒壊が殆どない。昭通市内では建物被害が殆ど見えない。
地盤常時微動観測は高感度速度計を用いて 100Hz のサンプリングレートで各地点に 30 分で実施した。H/V スペク
トルにより龍頭山鎮の強震動観測点の地盤の卓越振動数は約 4Hz であり、他の強震動観測点とほぼ同じであることが分かった。しかい、龍頭山の建物被害が地盤状況によって異なることが見られている。

キーワード: 中国雲南省昭通地震, 建物被害調査, 地盤常時微動, レンガ造建物
Keywords: Zhaotong Earthquake, Onsite Investigation, Ground Microtremor Measurement, Masonry Buildings

1/1

SSS25-02

会場:A04

時間:5 月 25 日 09:15-09:30

2014 年中国雲南省昭通地震 (M6.5) における建物被害と強震動特性の関係―その２：大加速度の生成メカニズム
Relationships between building damage and characteristics of strong ground motions during the M6.5 Zhaotong Earthquake-2
倉橋 奨 1∗ ; 入倉 孝次郎 1 ; 王 欣 2 ; 司 宏俊 3 ; 馬 強 4 ; 呉 浩 1
KURAHASHI, Susumu1∗ ; IRIKURA, Kojiro1 ; WANG, Xin2 ; SI, Hongjun3 ; MA, Qiang4 ; WU, Hao1
1

愛知工業大学, 2 東北大学災害科学国際研究所, 3 東京大学地震研究所, 4 中国地震局工程力学研究所
Aichi Institute of Technology, 2 IRIDeS, Tohoku University, 3 Earthquake Research Institute, The University of Tokyo, 4 Instituteof Engineering Mechanics, China Earthquake Administration

1

2014 年 8 月 3 日 16 時 30 分（北京時間）に中国雲南省昭通市を震源としたマグニチュード 6.5 の地震が発生した。この地震により震央付近の龍頭山鎮では、617 人の死者が出るなどの被害が発生した。この地震では、震央の極近傍の龍頭山鎮観測点（LLT）での最大加速度は 949gal となっており、特に、観測波形には、特徴的なパルス波形が 2 つ観測されている（図 (b)）。
地震動と被害との関係を解明するためには、この大加速度の地震動の生成メカニズムの解明が必要である。そこで、本研究では、観測記録を用いて本地震の震源モデルを構築し、大加速度の地震動の生成メカニズムの解明することを目的としている。
本地震の余震分布（中国地震局、2014）を図 (a) に示す。図のように、南東-北西方向に沿った面（ここでは断面 A と呼ぶ）と南西-北東方向に沿った面（ここでは断面 B と呼ぶ）が見られ、すなわち、くの字になっていることがわかる。
中国地震局による被害分布は、断面 A に沿っていることから、断面 A が主要な断面である可能性がある。
一方で、断面 A および断面 B において、遠地実体波による波形インバージョンが実施されている（中国地震局）。この結果では、どちらの断面が動いたかは明瞭に示されていない。
そこで、ここでは、まずは経験的グリーン関数法により、短周期地震動の大局的な震源モデルの構築を試みた。ただし、断層面等が明確でないため、震源極近傍の LLT 観測点については、除外して行った。利用した観測点は４観測点であるが、震源の南側に偏っている。要素地震の地震モーメントやメカニズムなど、震源パラメータの情報は非常に乏し
いが、地震モーメントはマグニチュードから、メカニズムは本震と近似していると仮定して実施している。小断層の足し合わせ数と応力降下量比は、本震と要素地震のスペクトル比から計算している。
以上、情報が乏しいながらも解析した結果、断面 A でも断面 B でも、遠くの観測点については、概ね表現できることがわかった。各断面は 10km × 10km、応力降下量は 10MPa 程度で表現できた。しかしどちらの断面が動いたか、両方とも動いたかの判断はできない。
以上のことから、震源極近傍の観測記録から生成メカニズムを解明していく必要があることがわかった。そこで次に、震源極近傍の観測点 LLT の波形に注目し、波形のパーティクルモーションから、S 波の到来方向を推定する。
P 波初動の水平動のパーティクルモーションを見ると、やや北西-南東方向に振動している。最初の地震動は観測点の北西方向から来ていることが推測される（図 (c) 左図）。次に、S 波の第 1 パルスおよび第 2 パルスのパーティクルモー
ションを検討する。その結果、第 1 パルスの初動は北西-南東方向に振動し（図 (c) 中図）、第 2 パルスは北東-南西方向に振動している（図 (c) 右図）ことが見て取れる。これらのパルスが、SH 波が卓越していると考えると、振動方向と直交する方向から波が到来していることなる。すなわち、第 1 パルスと第 2 パルスでは、到来方向が異なることがわかる。
余震分布と観測点との幾何学的な関係から考えると、第 1 パルスは断面 A から、第 2 パルスは断面 B から到来していることが考えられる。
今後、断層の位置を確定させるとともに、各断層面を仮定して震源極近傍の観測点 LLT 観測点も含めた、強震動生成域の推定を試みる。

キーワード: 中国雲南昭通地震, 大加速度, 強震動生成域
Keywords: Zhaotong Earthquake, large acceleration, strong-motion generation area

1/2

SSS25-02

会場:A04

時間:5 月 25 日 09:15-09:30

2/2

SSS25-03

会場:A04

時間:5 月 25 日 09:30-09:45

2014 年中国雲南省昭通地震 (Ms6.5) における建物被害と強震動特性の関係―その３：
龍頭山強震点周辺建物被害要因及び地盤震動特性
Relationships between building damage and strong ground motions characteristics during
the 2014 Zhaotong Earthquake-III
呉 浩 1∗ ; 王 欣 2 ; 司 宏俊 3 ; 党 紀 4 ; 馬 強 5 ; 林 国良 6 ; 入倉 孝次郎 1 ; 倉橋 奨 1
WU, Hao1∗ ; WANG, Xin2 ; SI, Hongjun3 ; DANG, Ji4 ; MA, Qiang5 ; LIN, Guoliang6 ; IRIKURA, Kojiro1 ;KURAHASHI, Susumu1
1

愛知工業大学, 2 東北大学, 3 東京大学, 4 埼玉大学, 5 中国地震局工程力学研究所, 6 雲南省地震局
Aichi Institute of Technology, 2 Tohoku University, 3 The University of Tokyo, 4 Saitama University, 5 China Earthquake Administration, IEM, 6 Yunnan Province Earthquake Administration

1

Many of the buildings near the Longtoushan (LTS) strong-motion station were heavily damaged or collapsed during the 2014
Zhaotong earthquake (Ms 6.5) which occurred on Aug. 3, 2014. The largest PGA with approximately 1 g was recorded at LTS
during the mainshock. In contrast, the damage of buildings was minor near other stations, and the strong-motion records there
were lower than 150 gal.
We conducted an almost three-day field survey near five strong-motion stations around the epicenter. We found that the causes
for heavy damage of buildings near the LTS station were complicated, such as geological failure (e.g. landslide), insufficiently
seismic resistant structural design, inappropriate construction, and site effect (transfer function and nonlinear effect), and so on.
We also found that the collapse of many buildings in the EW direction was obviously heavier than that in the NS direction near
the LTS station. It can be explained by the higher pseudo-velocity response spectrum in the range of 0.5 s to 1.0 s in the EW
direction.
We also conducted microtremor measurement at six sites around the LTS station to examine the variation of site effects in thisarea. We found that the predominant periods at heavily damaged sites were similar, i.e., 0.25 s, while the predominant periods atlight damaged sites were about 0.15 s.
In addition, we obtained several aftershock records as well as mainshock record, and the boring data near the LTS station.
These data can be used to clarify the site effect characteristic at the LTS station.

キーワード: 昭通地震, 強震動, 常時微動, サイト特性
Keywords: Zhaotong Earthquake, ground motion, microtremor, site effect

1/1

SSS25-04

会場:A04

時間:5 月 25 日 09:45-10:00

STRONG GROUND MOTION SIMULATION of THE 24 MAY 2014 NORTH AEGEANSEA EARTHQUAKE (Mw 6.9) in TEKIRDAG and SURROUNDING ARE
STRONG GROUND MOTION SIMULATION of THE 24 MAY 2014 NORTH AEGEANSEA EARTHQUAKE (Mw 6.9) in TEKIRDAG and SURROUNDING AREA
KARAGOZ Ozlem1∗ ; 地元 孝輔 1 ; 山中 浩明 1 ; Ozel Oguz3 ; citak seckin ozgur4
KARAGOZ, Ozlem1∗ ; CHIMOTO, Kosuke1 ; YAMANAKA, Hiroaki1 ; OZEL, Oguz3 ; CITAK, Seckin ozgur4
1

Department of Environmental Science and Technology, Tokyo Institute of Technology, Tokyo, Japan, 2 Department of Geophysical Engineering, Canakkale Onsekiz Mart University, Canakkale, Turkey, 3 Department of Geophysical Engineering, Istanbul
University, Istanbul, Turkey, 4 JAMSTEC
1
Department of Environmental Science and Technology, Tokyo Institute of Technology, Tokyo, Japan, 2 Department of Geophysical Engineering, Canakkale Onsekiz Mart University, Canakkale, Turkey, 3 Department of Geophysical Engineering, Istanbul
University, Istanbul, Turkey, 4 JAMSTEC
The Marmara Region (NW Turkey) was affected from destructive earthquakes since historical times. The North Anatolian
Fault with 1,200 km length is the main source of the earthquakes in the region. The most recent 1999 Kocaeli Earthquake (Mw
7.4) damaged residential and industrial areas of the large cities in Marmara. The future earthquake is expected on the northwestern segment of the fault close to the city of Tekirdag.
In this study, we simulated strong ground motion records of the 24 May 2014 North Aegean Sea Earthquake (Mw 6.9) inTekirdag and surrounding area. We used one-dimensional homogeneous horizontal layer model at each AFAD (Republic of
Turkey Prime Ministry Disaster & Emergency Management Presidency Earthquake Department) strong motion station site inTekirdag and surroding area (i.e. 5904, 5902, 5907, 5906), Canakkale (i.e. 1701, 1710), Gokceada (1711), Edirne-Enez (2201).
We combined our shallow model (0-250 m) that obtained from our previous microtremor explorations (Karagoz et al., 2014) andthe deeper parts were taken from pervious crustal studies in the region. The outer fault parameters of the mainshock (seismic
moment, strike, dip, and rake) were determined by previous focal mechanism solutions studies while the inner fault parameters
were estimated by following the recipe of Irikura and Miyake (2011).
The fault plane (background) and asperities were divided into several subfaults that were assumed as single double-couple
point source. We defined appropriate Kostrov-like slip-velocity function (modified by Nakamura and Miyatake, 2000) for the asperity areas to simulate high frequency ground motions. The synthetic seismograms were obtained using a deterministic discrete
wave number method for each sub-fault and were summed to get full waveform at the station around the epicentre in a broadband
frequency range (0.1-10Hz).
The simulated peak ground velocities at the ground surface were estimated by multiplying the simulated ground motion at the
top layer of Vs=780 m/s from the discrete wave number method with 1D amplification factors of S-waves in the shallow soil
layers derived from the microtremor explorations. For validation, the results converted accelerations and were compared with
S-wave portions of the recorded acceleration waveforms at the strong motion stations.

キーワード: Discreete wave form method, earthquake waveform simulation, Gokceada, site effect, Tekirdag
Keywords: Discreete wave form method, earthquake waveform simulation, Gokceada, site effect, Tekirdag

1/1

SSS25-05

会場:A04

時間:5 月 25 日 10:00-10:15

強震波形から推定した 2014 年長野県北部地震の震源過程
Source Rupture Process of the 2014 Northern Nagano Earthquake Estimated by Strong
Motion Data
浅野 公之 1∗ ; 岩田 知孝 1 ; 久保 久彦 1ASANO, Kimiyuki1∗ ; IWATA, Tomotaka1 ; KUBO, Hisahiko1
1
1

京都大学防災研究所
Disaster Prevention Research Institute, Kyoto University

2014 年 11 月 22 日 22 時 8 分に長野県北部の北安曇郡白馬村付近で発生した MJM A 6.7 の地震（以下、2014 年長野県北部の地震）では、最大震度６弱の地震動が観測された。Global CMT Project や F-net から公開されているモーメントテン
ソル解によれば、北北東−南南西の走向をもつ東側隆起の逆断層で発生した地震と考えられる。また、この地震は糸魚川−静岡構造線断層帯の神城断層の活動に関連していると報告されている（例えば、地震調査研究推進本部地震調査委
員会，2014）。本研究では、この地震の震源像の実態を把握するため、強震波形記録を用いたインバージョン解析によって、2014 年長野県北部の地震の震源過程を推定した。
解析には震源域周辺の K-NET、KiK-net、気象庁、長野県の強震観測点から 12 点の波形記録を用いた。加速度記録を
速度に積分し、0.05 − 1Hz のバンドパスフィルターを適用した波形記録から、S 波到着 1 秒前より 15 秒間を抽出した波形をインバージョンのためのデータとした。この地域の地下構造は複雑であり、単一の速度構造モデルを用いてグリーン関数を計算することは適切ではない。本研究では、地震調査研究推進本部の全国１次地下構造モデル（Koketsu et al.,2012）から観測点ごとに直下の一次元構造を抽出した速度構造モデルを使用し、離散化波数法（Bouchon, 1981）及び透過・反射係数行列法（Kennett and Kerry, 1979）によりグリーン関数を計算した。
余震分布（例えば、防災科学技術研究所, 2014）及び地表地震断層の分布から断層面が複雑であることが推測される。
本研究では、震源断層の北部と南部で傾斜の異なる 2 枚の平面断層を組み合わせたモデルで断層モデルを表現した。南部は走向 25°、傾斜 65°、北部は走向 25°、傾斜 50°と仮定した。南部の断層面の上端は神城断層の地表トレースに対応するようにしている。全体の断層サイズは長さ 22 km、幅 14 km とし、これを 2 km×2 km の小断層で分割した。各小断
層には時間幅 1 秒の smoothed ramp 関数を 0.5 秒ずらしで 6 個並べて、その小断層でのモーメント時間関数を表現した。
インバージョン手法はマルチタイムウィンドウ線形波形インバージョン（Hartzell and Heaton, 1983）を用い、時空間のすべりの平滑化（Sekiguchi et al., 2000）の強さ及び第一タイムウィンドウの破壊伝播速度は ABIC 最小規準（Akaike,1980）により決定した。すべり方向は 45° ±45°の範囲に拘束している。
破壊開始点から約 5km 北のやや深いところにすべり量の大きな領域が推定され、最大すべり量は 1.8m と推定された。
この領域では周辺に比べ余震数が少ない。各機関によるモーメントテンソル解のセントロイドが北側深部に推定されていることとも整合している。また、地表に近いところでも 0.3〜0.5m 程度のすべりがあり、これらのすべりは地表地震
断層とも関係していると考えられる。全体の地震モーメントは 3.85×1018 Nm、モーメントマグニチュードは 6.3 である。
断層全体の平均すべり量は約 0.4m である。これらの断層パラメータを既往の内陸地殻内地震と比較したところ、これまでの日本の内陸地殻内地震の断層パラメータのばらつきの範囲内であった。地表地震断層に関しては、地震規模に比べ、地表変位量が大きいことが指摘されているものの、強震動に関わる震源過程としては、平均的な地震であったといってもよい。
今後も詳細な余震分布や地表地震断層調査の結果が報告されると考えられる。これらの情報を踏まえ、震源モデルの見直しを図っていく必要がある。
謝辞：独立行政法人防災科学技術研究所強震観測網（K-NET 及び KiK-net）、気象庁震度計、東京大学地震研究所首都圏強震動総合ネットワーク（SK-net）より公開されている長野県震度情報ネットワークシステムの強震波形、気象庁一元化震源カタログを使用した。データ公開に関わられる皆様に記して感謝する。

キーワード: 2014 年長野県北部地震, 震源過程, 強震記録
Keywords: the 2014 northern Nagano earthquake, source process, strong motion data

1/1

SSS25-06

会場:A04

時間:5 月 25 日 10:15-10:30

2014 年長野県北部の地震の余震と微動の観測による白馬村の地盤増幅特性の評価
Site amplification in Hakuba from microtremor and aftershock observation of the 2014
Northern Nagano Earthquake
地元 孝輔 1∗ ; 山中 浩明 1 ; 佐口 浩一郎 1 ; 津野 靖士 2 ; 盛川 仁 1 ; 飯山 かほり 1 ; 後藤 浩之 3
CHIMOTO, Kosuke1∗ ; YAMANAKA, Hiroaki1 ; SAGUCHI, Koichiro1 ; TSUNO, Seiji2 ; MORIKAWA, Hitoshi1 ;IIYAMA, Kahori1 ; GOTO, Hiroyuki3
1
1

東京工業大学大学院総合理工学研究科, 2 鉄道総合技術研究所, 3 京都大学防災研究所
Tokyo Institute of Technology, 2 Railway Technical Research Institute, 3 DPRI

2014 年 11 月 22 日に発生した長野県北部の地震では、最大震度 6 弱が観測され、震度 5 強が観測された K-NET 白馬から南におよそ 5km の白馬村堀之内や三日市場では建物被害が顕著であった。そこで、11 月 24 日から 12 月 4 日にかけて K-NET 白馬から白馬村堀之内、三日市場の周辺で余震観測と微動探査を実施し、地下構造による地盤増幅特性について調べた。
臨時余震観測は、加速度計 JEP6A3 とデータロガー LS7000XT または LS8800 を用いた 13 観測点で 10 日ほど行った。
震度 1 以上の余震は 30 ほど観測され、最大加速度にして 30gal 程度の余震記録が得られた。建物被害が顕著であった堀之内と三日市場からその西側の沖積低地にかけて大きな加速度が得られ、長い後続位相も確認される。また、速度応答スペクトルは K-NET 白馬に比べて周期 0.5 から 1 秒程度で大きくなっている。三日市場の山裾の余震観測点に対するス
ペクトル比を求めたところ、堀之内では周期 0.5 秒から 1 秒程度にピークを持つことがわかった。
各余震観測点では 20m 程度以下の微動アレイと単点微動観測を行った。さらに、堀之内近くでは、神城断層の東側と西側で最大 680m ほどの微動アレイ観測も実施し、SPAC 法を適用して位相速度の推定を行った。20m 程度のアレイ観測
記録からは、3Hz から 30Hz 程度までの位相速度分散曲線が得られ、K-NET 白馬で最も速い。堀之内や三日市場では、
100 から 200m/s 程度の位相速度が得られた。大きなアレイでは周期 1 秒以上で 1000m/s 以上の位相速度が観測された。
位相速度の逆解析により S 波速度構造を推定し、K-NET 白馬では S 波速度 300m/s 程度の表層が 10m 未満であったものの、堀之内や三日市場では 200m/s 以下の低速度層が 10m 程度存在し、さらに 400m/s 以下の層が 50m 程度と深いことがわかった。大きなアレイ観測により得られた位相速度の逆解析では、断層の西側では東側に比べて深い構造であることがわかった。S 波速度 1000m/s 以上の層までの深さは、西側で 700m ほど、東側で、400m ほどであった。K-NET 観測点の 1km 西にある KiK-net 観測点の検層結果では、それが 100m 以下であることから、被害が大きかった地域付近では深部地盤も深いことがわかる。
推定された S 波速度構造によって S 波地盤増幅特性を求めた。表層地盤による地盤増幅は、K-NET 白馬では周期 0.1
秒程度で卓越するものの、堀之内や三日市場では周期 1 秒程度で卓越する。余震観測記録のスペクトル比と同様にして地盤増幅特性の比を求めて比較したところ、卓越周期はおおよそ類似したものとなった。しかし、概して観測されたス
ペクトル比のほうが大きく、深部地盤による増幅効果の影響なども考慮する必要がある。

キーワード: 2014 年長野県北部の地震, 地盤増幅特性, 余震観測, 微動, 神城断層
Keywords: 2014 Northern Nagano Earthquake, Site Amplification, Aftershock observation, Microtremors, Kamishiro Fault

1/1

SSS25-07

会場:A04

時間:5 月 25 日 10:30-10:45

2014 年 11 月 22 日長野県北部の地震による長野善光寺の石灯籠の被害と強震動
Damages of Stone Lanterns at Zenkoji Temple, Nagano, Caused by Northeastern Nagano
Earthquake, November 22, 2014加藤 護 1∗
KATO, Mamoru1∗
1
1

京大院人間・環境学
GSHE, Kyoto University

2014 年 11 月 22 日に長野県北部を震源とする気象庁マグニチュード 6.7 の地震が発生した。震源地付近で震度 6 弱の
強い揺れとなり、白馬村や小谷村などで家屋等に多大な被害が出た。長野市中心部に存在する長野善光寺でも鐘楼の土台の破損や多数の石灯籠の倒壊など被害が出た。この地域の歴史地震のひとつに 1714 年の地震があるが、この際にも善光寺の石灯籠の被害が生じたとの記録が残っている。白馬盆地の被害分布を考え合わせると今回の地震と揺れの特徴が類似していた可能性が示唆される。
糸魚川・静岡構造線北部で発生する地震による地震動予測を精緻化するうえで今回の地震で生じた長野善光寺に代表される長野市街の被害を整理しその原因となった強震動を検討することは重要であると考えられる。以上を踏まえわれ
われは 11 月 22 日の地震によって生じた長野善光寺の石灯籠の現地調査を行った。本発表では調査の結果を整理し、この被害を引き起こした強震動の特徴に関する予察を示す。
被害の調査は地震翌日の 11 月 23 日に行った。11 月 30 日及び 12 月 1 日に追加調査を行ったが、既に一部で復旧が始まっており両日の結果に差異がある場合は 23 日に確認した被害を採用した。われわれは 1847 年善光寺地震による石灯籠の被害の把握を目的として 2013 年夏に善光寺石灯籠の悉皆調査を行っている（日岡，2014; 加藤・日岡，2014）。今回の調査は 2013 年の調査記録をベースとしている。日岡 (2014) などでは石灯籠のみを調査対象としていたが、今回は石塔の被害も確認できたので合わせて報告する。
なお、長野市中心部では長野地方気象台で震度 5 強、防災科研 K-net NGN004 観測点で震度４が観測されている。
被害発生後の改変の可能性がほぼない石灯籠 65 基、石塔 12 基について部位の落下や転倒を確認した。石灯籠の被災率は約 3 割 5 分程度である。この他に部位の回転などの被害も確認された。被害を受けた石灯籠などの分布は境内内で偏りがないようにみえる。落下や倒壊の方位は偏って分布しており、被害を引き起こした強震動の特徴を反映していると考えられる。
石灯籠などの小型構造物は一般に 0.3〜0.5Hz 程度の固有周期を持つ。善光寺に近い長野地方気象台の波形はこの周期帯の成分に富んでいたことから長野善光寺の石灯籠の被害はこの短周期の揺れを反映している可能性が考えられる。
謝辞：気象庁と防災科研で公開されている強震波形記録を用いました。

キーワード: 強震動, 2014 年長野県北部地震
Keywords: String Motion, Northern Nagano Earthquake, 2014

1/1

SSS25-08

会場:A04

時間:5 月 25 日 11:00-11:15

2011 年東日本大震災に伴う死者の年令等依存性 (5) 死亡率算定法の再改訂
Age-dependent Mortality in the 2011 East Japan Giant Earthquake (5) Additional Revision of the Current Equation
太田 裕 1∗ ; 志垣 智子 2 ; 小山 真紀 3
OHTA, Yutaka1∗ ; SHIGAKI, Tomoko2 ; KOYAMA, Maki311

東濃地震科学研究所, 2 高齢者住宅研究所, 3 京大院医学研究科
Tono Res. Inst. Earthq. Science, 2 Inst. Elderly Housing Sci., 3 Medical Inst.,Kyoto Univ.

1. 前承け前報までは地震・津波等に伴って発生する急性期の死者について，年令区分別の発生状況を調べ, 算定式の構成を行ってきた１）．その当初段階では伝統に従って，被災域における年令区分別［死亡者数/母数人口］を基本式としてきた．その結果，この関係式は高令者の死亡実態はよく反映するものの，乳幼児・若年者については過少評価となっている場合が間々あり，年令区分でみた死者発生の実態を正しく評価し得ていないことを例示してきた（周知の片田 2）が活用する算定式も上記の伝統手法と同等であり，同様の問題を内蔵している）。そこで， 2 つの算定法を新たに導入することで事態の改善を試みた.1 つは「平常時，1 年間」の年令区分別死者数と対比する尾崎法 3) を小改訂したものであり，他の１つは経済学者 Sen が主唱する『乳幼児. 若年者の死亡は，相応の生存ができたならばあったであろう可能性の全てを剥奪するものである４）』との主唱に添う形で 年令の逆数 を重みとする損失量を導入した新しい構成式を提示した．そして，これらの式を活用することで，伝統評価法がもつ問題点のうち，特に乳幼児・若年層に対する評価の是正を試みたものである．
2. 年令区分別死亡率の算定法の再改定
本論はこの延長線上にあり，先ず上記の改訂尾崎法がもつ欠点について再度の改訂を行った．すなわち，伝統の算定法では死亡率を最終出力としているのに，前報までの改訂版尾崎式では，縦軸が死者の絶対数を与える形に止まっており，両者の直接対比が出来ないままとなっていた．本報では先ず，この点に注目し，縦軸を死亡率に変換・表示できるよう，さらなる改訂を試みた．その結果，この再改訂版算定法を主体とすること，そして，余命特性損失評価法を
適宜介在させることで地震・津波等の災害に伴う急性期の死亡率を年令区分別にみる算定法は一応の到達点に達し得た
ものと，筆者らは考えている次第である．
以下では, 今回得た再改訂版算定式を活用して，県別および市町村別にみた年令依存性についていくつかの考察を行っている．
3. 県および市町村単位でみた年令区分別死亡率1) 岩手県/釜石の奇跡
今回の地震津波に対して，東北３県のうち特に岩手県においては乳幼児〜若年層の高い生存率がよくいわれている．このことは在来の伝統法によった場合でも，年令区分別平滑化曲線とのズレ（実測値＜平滑曲線）を見ることで容易に想像できるところである．つまり，岩手県（詳しくは釜石市と周辺地域）の奇蹟といわれる根拠をあたえている 2）．しかし，
今回得た再改訂版を使って算定し直してみることで「釜石等周辺地域も乳幼児・若年者の，津波による死亡率は平常時・１年間のそれ以下とはとてもいえない」結果となっていることが確認出来る．つまり大方の努力にも拘わらず，この年代の死亡率はなお相当に高いことを示している．同様の結果は，期待特性損失量の算定結果からもハッキリと確認できる．
2）その他の特異変化
同様の視点で，県・市町村別に実際の年令依存性曲線と平滑曲線とのズレに注目してみると，特異事例のいくつかが浮上する．年令を昇順にみていくと，県別のみならず，市町村別にみても多くの場合，20〜30 代で実測値＞平滑値（乳幼児等とは逆のズレ）となるゾーンが目立つ．これは若者〜壮年なるが故の対津波公務活動（沿岸域水門の管理とか避難誘導等）に従事した結果の Negative な反映であろうか．一方，さらに高令者側では 3 県共に 65 才前後で死亡率が局所的に急増していることを示すズレ（実測値＞平滑曲線）が目立っている．これについては東北沿岸域の家族構成（複数
世代の多少等）とか相互支援を考慮することで実態解明への手掛かりが得られることを期待している．
4. 終りに
従来のように，伝統の死者発生率を単純にみた場合，年令に対して右上がりの形，Ｊ字型分布の関係が得られ，乳幼児の死亡率が低いようにみえたとしても，視点を変えて平常時のそれと比べてみることで彼らの死亡率が未だ相当に高いことが，今回までの試算で初めて明瞭になった．今回までの成果は，先に得た U 字型の分布 1) についても，今一度立ち入った考察が必要なことを示唆している．
文献等
1) 太田・小山,2011 年東日本大震災に伴う人間被害の激甚性 (2〜4)2013,14 年春 J ｐ GU 大会.
2) 片田，人が死なない防災，集英社新書，2012．
3) 尾崎，地震災害時および災害後の健康被害，厚生の指標，59，2012．
1/2

SSS25-08

会場:A04

時間:5 月 25 日 11:00-11:15

4)Amatya Sen, Wikpedia で Sen の項を検索.

キーワード: 東日本大震災, 年令区分別死亡率, 算定法改訂, 災害弱者
Keywords: East Japan Earthquake, Age-band specific Mortality, Revision of Evaluation Equation, Disaster-Vulnerable Persons

2/2

SSS25-09

会場:A04

時間:5 月 25 日 11:15-11:30

S 波ピーク領域と距離別の速度を使って解析した２０１１東北沖地震の多重震源の解析
Multiple Seismic Origins of the 2011 Tohoku Earthquake Analyzed by S-Wave Peak andRegions
沼倉 政幸 1∗
NUMAKURA, Masayuki1∗
1
1

元 千葉県 佐倉東高校 教諭
Last position Sakura Higashi High School

1, Introduction
Within the field of research on the 2011 Tohoku earthquakes, analyses of the vicinity of the ocean trench axis have maderemarkable progress. In contrast, opinion is divided regarding regions closer to land. One example concerns the question ofwhether tremors that reached the Kanto region originated from the ocean trench axis. Many studies emphasize land-based
seismic origins for these tremors, but the precise locations from which these tremors originated remains to be determined indetail. Considering that this earthquake was a ones-in-several-hundred-years event, it can hardly be said that our understandingof the earthquake as a whole has progressed to a fully satisfactory state.
The S-wave and P-wave data available for this earthquake are mixed together and hence difficult to use. Multiple largeamplitude contributions following S-wave initial shocks, S-wave peaks, were present. These data may be used for analyses,
despite their mixed nature, if they can be well separated. Due to the paucity of previous research on this subject, we have tested
our methods on other large earthquakes in Japan. Method is indicated to 3-1-2 Results of validation.
2, Method
2-1 We determined epicenters by following the Omori method of drawing three circles on map.
2-2 We use Omori’s formula. (However, we use distant-dependent velocities.)
r=k×t k1=(Vp×Vs)/(Vp-Vs)
2-3 We used S-wave peak travel times and velocities. The method is indicated to 3-1-2.
3, Results
3.1 Validation.
3-1-2 2004 Niigata Chuetsu earthquake
Epicenter 37.29N 138.87E Depth 12km Start ; 17:56:00
Ojiya Distance13.9km Traveltime S-wave peak 8.73s
Velocity peak3.19km/s Acceleration0.365 Constant1.16Time arrival 56:08
Tokamachi Distance24.1km Travel time S-wave peak 12.60s
Velocity peak 3.83km/s Acceleration 0.304 Constant1.16 Time arrival 56:12
Kashiwazaki Distance31.0km Travel time S-wave peak 14.90s
Velocity peak4.16km/s Acceleration 0.280 Constant1.16 Time arrival 56:14
Travel time is peak time at many stations. (2)r=1/2×a×tˆ2 (3)r=1/2 √ (ktˆ3 )(4) r=1/2×v×t (5)k=v×a k=4rˆ2/tˆ3 (6)v=2r/t v=dr/dt (7) a=2r/tˆ2 a=dv/dt
It is sufficiently possible to calculate in this equation (2), but coefficient ″ a ″ is not constant. Therefore, I prepare theequation (3). The equations (2) and (3) are function that the distance and time of the hypocenter and observation point. In thoseequations, the function (5), which is in inverse proportion to ″ v ″ and ″ a ″, is contained.
3-2 Seismic origins of 2011 Tohoku Earthquake
P-1-2 (Epicenter announced by United States Geological Survey)
Epicenter 37.291N 138.867N Depth 30km Start: 46minites, 24 seconds
Oshika Distance 81.5km Travel time S-wave peak 43.3s
Velocity peak 3.76km/s Acceleration 0.0868 Constant 0.327
Time arrival 47:07
Utatsu Distance 94.1km Travel time S-wave peak 47.6s
Velocity peak 3.95km/s Acceleration 0.0829 Constant 0.327
Time arrival 47:11
Tsukidate Distance 130.9km Travel time S-wave peak 59.3s

1/2

SSS25-09

会場:A04

時間:5 月 25 日 11:15-11:30

Velocity peak 4.41km/s Acceleration 0.0744 Constant 0.327
Time arrival 47:23
P-3, Epicenter 38.045N141.47E Depth30km Start: 47minites, 37 seconds
Oshika Distance 41.8km Travel time peak 20.0s
Velocity peak 4.18km/s Acceleration 0.209 Constant 0.873
Time arrival 47:57
Kamaishi Distance 143.6km Travel time peak 45.4s
Velocity peak 6.33km/s Acceleration 0.139 Constant 0.873
Time arrival 48:22
Hitachi Distance 180.2km Travel time peak 53.0s
Velocity peak 6.80km/s Acceleration 0.128 Constant 0.873
Time arrival 48:30
P5 (May be divided into 5 epicenters in all)In the first phase, P1-2 and P1-3 gave a strong motion to Miyagi and Iwate. In the third phase, Max acceleration of Tsukidate
and Oshika was observed by P3. It is both epicenter third phase and trench axis that gave a strong motion to Kanto.

キーワード: S 波ピーク, ２０１１東北沖地震
Keywords: S-Wave Peak, 2011 Tohoku Warthquake

2/2

SSS25-10

会場:A04

時間:5 月 25 日 11:30-11:45

多段階不均質震源モデルに基づく 2011 年 Mw9.0 東北地震の強震動の再現
Simulation of strong ground motions from the Evaluation of the 2011 Mw 9.0 Tohoku
earthquake
入倉 孝次郎 1∗ ; 倉橋 奨 1
IRIKURA, Kojiro1∗ ; KURAHASHI, Susumu1
1
1

愛知工業大学地域防災研究センター
Aichi Institute of Technology, Disaster Prevention Research Center１．はじめに2011 年東北地震の震源近傍域で得られた強震加速度波形に、顕著な複数のパルス的波群が見られる。これらの波群は、
強震動観測記録を用いたセンブランス解析や逆伝播法解析により、震源断層の中で西端に近いプレート沈み込みのダウ
ン・ディップに沿って、南北方向に並ぶ 4〜5 つの位置から生成された、と考えられる。
我々は、過去の起こった中小地震の強震動記録を経験的グリーン関数として特性化震源モデルを仮定して、本震の地
震動のシミュレーションを行い、観測記録との比較により、5 つの強震動生成域からなる短周期生成モデルを求めた。この解析で、震源近傍域の強震動記録で見られる顕著な波群は、沈みこむプレートのダウン・ディップに位置する強震動生成域から生成された、ことを明らかにしなった。さらに、ここで得られた強震動生成域モデルと数値的グリーン関数
を用いて、より長周期の地震動のシミュレーションを行い、この短周期震源モデルの有効な周期範囲は 0.1˜10 秒である
ことを確認した。さらに、各波群に見られるパルス的波形をより精度良く再現するには、強震動生成域は一様な実効応力ではなく、強震動生成域ごとに立ち上がりに大きな実効応力を有する不均質なモデルの設定が必要なことを確認した。
切迫性の指摘されている南海トラフ地震など、すべり込み型巨大地震に対する防災・減災対策を効果的に進めるためには、信頼性の高い強震動や津波の予測技術の確立が必要とされている。本研究は、2011 年東北地震の強震動観測デー
タから明らかになった多段階不均質モデルに基づき、信頼性ある強震動予測手法について 検討することを目的としている。
２．2011 年東北地震の強震動再現のための多段階不均質震源モデル
この地震の短周期震源モデルは、多くの著者により試みられ、強震動生成域の数は 4〜5 個、それらの位置はダウン・
ディップの近くに位置するなど、共通性のある結果が得られている (Asano and Iwata, 2012; 佐藤, 2012; 川辺・釜江, 2013;
Kurahashi and Irikura, 2013)。しかしながら、強震動観測点は、震源断層の西側のみで、震源域の広さ（400 x 200 km2）に比べると、観測点の数は十分でないため、強震動生成域の位置は必ずしもユニークには決まらない。 Kurahashi and
Irikura(2013) は、強震動生成域内で応力パラメータの不均質を設定しているが、他の著者が示しているように、強震動生
成域は一様と設定しても、観測波形に見られる 5 つの波群とほぼ同様のスペクトル特性をもつ合成波形のシミュレーショ
ンは可能である。また、これらの SMGA（強震動生成域）モデルを用いると、工学的に重要な 2 〜 10 秒の長周期地震動の再現に有効であることがわかってきた（川辺, 2012; 倉橋・入倉, 2014）。
野津 (2014) は、これまでに発表されたこの地震の短周期震源モデルを用いて強震動予測の有効性について検討した。
野津らによる SPGA（強震動パルス生成域）モデル、疑似点震源モデルと 3 つの SMGA モデルを用いたシミュレーショ
ン結果を比較した結果、構造物に影響の大きい周期 1〜2 秒の卓越するパルス波形の再現には、SPGA モデルや疑似点震源モデルは有効だが、SMGA モデルは過小評価になると結論している。一方で、野津らのモデルは、強震動生成域を極めて狭くしているため、構造物への影響の大きい破壊的強震動の特徴である破壊の指向性効果や 2〜10 秒の長周期地震動の再現が可能かなど、巨大地震に対する強震動予測としての有効性に問題がある。
もう１つの問題として、この地震の強震動波形の特徴の１つであるパルス的波形の再現には、破壊開始点付近に大きな実効応力をもつ不均質モデルが必要と考えられる。そのためには、適切な経験的グリーン関数の選定が必要とされる。
Kurahashi and Irikura (2013) では、強震動のシミュレーションのために、経験的グリーン関数として 2005 年宮城沖地震 (Mw 7.2) の後半部分を用いている。この地震の強震動再現のための震源モデルは、震源域で発生した小地震（M 4.1）
記録を経験的グリーン関数として用いると、東西に並ぶ２つの SMGA からなっていることがわかる。それらの SMGAの大きさは約 8 x 8 km2 である。それぞれ SMGA 内で、破壊は東から西（海側から陸側）方向に進んでおり、震源に近い宮城県内の観測点では破壊伝播の指向性パルスが顕著に見られる。これは先行研究の Suzuki and Iwata (2008) とほぼ同じ結果である。
このことから、東北地震の強震動は、Mw7.2 の宮城県沖地震の記録を用いなくても、Mw ４クラスの地震記録を用いて、約 8 x 8 km2 の要素地震を合成し、それを経験的グリーン関数として、Kurahashi and Irikura (2013) の提唱する多段階不均質モデルの手順に従って強震動シミュレーションを行えば、観測に見られるパルシブな形状をもつ強震動でかつ 0.1 ˜10 秒の広帯域強震動の再現が可能となる。

1/2

SSS25-10

会場:A04

時間:5 月 25 日 11:30-11:45

キーワード: 2011 年東北地方太平洋沖地震, 沈み込み地震, 強震動, 特性化震源モデル, 強震動生成域, 経験的グリーン関数法
Keywords: the 2011 Tohoku earthquake, subduction earthquake, strong ground motion, characterized source model, strong
motion generation area, the empirical Green’s function method

2/2

SSS25-11

会場:A04

時間:5 月 25 日 11:45-12:00

神奈川県における関東大震災の慰霊碑・記念碑・遺構調査 (2) 県西部＋熱海・伊東
Field Survey for the Memorial Matters from the 1923 Great Kanto Earthquake in Western
Kanagawa Prefecture, Japan
武村 雅之 1∗
TAKEMURA, Masayuki1∗
11名古屋大学減災連携研究センター
Disaster Mitigation Research Center, Nogoya-Univ.

関東大震災の研究は震災直後に被害調査を含む多くの研究がなされたあと、1990 年代から 2000 年にかけて武村らを中心として、地震計による地震記録の発掘と解析による本震・余震の解明、被害データの整理分析による詳細震度分布の評価などが行われた。また東京両国の被服廠跡の火災による大量死や小田原市根府川・米神における土砂災害など関東
大震災を構成する個々の大惨事についても実態解明にあたってきた [武村 (2009)]。一方、市民がより身近なものとして関東大震災を理解し防災意識の向上に役立てるべく、当時の人々が記した日記などに科学的な解説を加えることによって、体験談のもつ臨場感と科学的客観性を融合させ時代を超えて震災を理解できるようにする取り組みも行ってきた [武村（2005、2008）]。
さらに現代人が震災を日常的に意識するために有効なものとして、様々なところに人知れず建つ慰霊碑や記念碑や震災の爪痕を残す遺構などがあることに目を付けた。関東大震災の慰霊碑や記念碑については力武（1994-1997）、神奈川県立歴史博物館（2003）、上西 (2012) や Web サイト「関東大震災の跡と痕を訪ねて」などでも紹介されている。各自治体の教育委員会などによる石造物調査もある。しかしながらこれらの調査結果はそれぞれ限定的であり誤りも多く、震災を生き抜いた当時の人々の意思を知り、現代人に震災を伝えるには不十分であると判断した。このため、従来の調査結果をすべて踏まえさらに調査対象を広げ、正確を期すため全ての対象物の現地調査を行うことにした。
武村（2012）による『関東大震災を歩く』吉川弘文館では、東京都 23 区内を対象として、震災の慰霊碑、記念碑、多くの犠牲者が出た場所や逆に多くの避難者の命を救った場所の現在の様子、震災で破壊された跡またはその再生を伝えるもの、復興過程で生まれた建物や施設、震災後の帝都復興事業による土地区画整理で郊外へ移転を余儀なくされた多くの寺院、さらには関東大震災以前から江戸・東京をたびたび襲った自然災害についての記念碑や遺構も調査しまとめた。調査地点は 180 カ所、対象物は 260 にのぼった。
今回その調査を神奈川県下（静岡県伊豆地方も含む）に広げた。神奈川県を西部、中部、東部の 3 地域に分け、2014 年
3 月に中部地域、2015 年 3 月には西部地域ならびに熱海・伊東地域の調査結果をまとめた。今回まとまった西部地域では、慰霊碑 31、復興碑など記念碑 58、その他エピソードを伝えるもの 38、消滅したもの 1、他の災害 42 の合計 170 を対象とした。他の災害として主なものは元禄地震や富士山の宝永噴火や関東地震の影響で起った水害である。現状の調査結果はすでに防災講演会や市民向け現地ツアー、文書館や博物館、防災団体へのデータ提供などの形で活用を進めている。
史上最悪の自然災害である関東大震災の慰霊碑・記念碑や被害の実態を今に伝える遺構、ならびに復興の記念物などを通じて、当時の人々の意志を受け止める。そのことを通じて現代社会に生きる人々の震災への理解を深め、防災意識の向上に繋げること。それこそが本調査研究の究極の目的である。本研究調査は JSPS KAKENHI 25350496 の助成を受けたものである。

キーワード: 関東大震災, 慰霊碑, 神奈川県
Keywords: Great Kanto Earthquake, Memorial tower, Kanagawa Prefecture

1/1

SSS25-12

会場:A04

時間:5 月 25 日 12:00-12:15

地震ザブトンを用いた地震動体験の防災効果測定に向けた検討
Effects on disaster prevention by The Portable Earthquake Simulator ”Jishin-The-Vuton”
東 宏樹 1∗ ; 内藤 昌平 1 ; 藤原 広行 1
AZUMA, Hiroki1∗ ; NAITO, Shohei1 ; FUJIWARA, Hiroyuki1
1
1

防災科学技術研究所
NIED

１． はじめに地震ザブトンは地震動の忠実な再現を目指して作られた可搬型シミュレータである。兵庫県南部地震の JR 鷹取波（震度 7）や東北地方太平洋沖地震の K-NET 仙台波（震度 6 強）のような観測波形をはじめ、新宿の高層階で実測されたデー
タ、想定東海地震の新宿高層 30F 建てのような想定の長周期地震動にいたるまで 14 の波形を、VUTON クローラという機構を持ったロボットによって水平 2 次元で再現できる。体験時には同一の揺れによって一般的な部屋がどうなるかを振動台で撮影された映像を、音とともに同時再生し、よりリアルな地震時の状況再現のための工夫がなされている。
２． 調査概要
防災科研では 2014 年から地震ザブトンの運用を正式に開始し、所内見学や近隣で行われる各種防災訓練、ワークショッ
プなどでの利用に供している。調査はこれらのイベント等における体験者に対し体験終了時のタイミングで簡単な質問用紙を記入してもらうことで行なった。質問用紙は一人あたりの地震動体験の時間を考慮し、記入待ち行列が発生しな
い概ね 1 分程度で記入を終える５?６項目程度の無記名のチェックシートで構成した。
３． 調査結果
以下に 2014 年 8 月 2 日から 12 月 7 日まで所内見学を含む断続的な 14 日間のイベント時に行われた調査結果の概要を示す。身長等の安全上の理由で地震ザブトンを体験できない未就学児を除く 12 歳以下の子供から 66 歳以上の高齢者まで幅広いデータを得たが、割合としては 56 歳?65 歳が 32%、66 歳以上が 42%と高齢者側に偏りが見られ、これは実際に防災訓練等に参加する年齢構成を反映しているものと思われる (N=137)。体験後の率直な感想としてはこわかった (40%)・
おもしろかった (31%)・ひとにすすめたい (29%) がほぼ同割合の結果となっている (N=220, 複数選択可)。体験した地震
動は震度 7 が体験できる兵庫県南部地震の JR 鷹取波 (29%) と新潟県中越地震の K-NET 小千谷波 (21%) が多いが、直下型以外にも東北地方太平洋沖地震で観測された各地の波形（あわせて 15%）や長周期地震動（あわせて 17%）など適度な分散が見られた。これら地震動の体験後に地震の対策を「する」と答えた人は 91%にのぼり、地震ザブトンの地震動
体験は耐震対策の意識を向上させる効果が見られた（図 1）。
４． おわりに
実際に地震の揺れを体験することで地震対策を行う意識が向上することが明らかになった。これらは他の災害と比べても比較的低頻度である地震において、平時にわがこと化を促進し対策行動を引き起こすことが困難であるという既往の問題の解決に有効なアプローチであることを示唆している。今後はこうした検証を高度化しつつ他機関とも連携して母数を増やし、上記課題の解決に向けた取り組みを加速していきたい。さいごに実験にご協力いただいた防災科研アウトリーチグループと阿見町、つくば市の皆様に、御礼を申し上げます。ありがとうございました。
参考文献）
松平 昌之 (白山工業 (株)) 山口 龍介 (東工大) 平山 義治 (白山工業 (株)) 王猛 (東工大) 吉田 稔 (白山工業 (株)) 翠川三郎
(東工大) 広瀬茂男 (東工大)：全方向移動型可搬式地震動シミュレータ 地震ザブトン の開発、第 27 回日本ロボット学
会学術講演会
キーワード: 地震ザブトン, 地震動体験, 防災効果測定
Keywords: ”Jishin-The-Vuton”, Strong motion experience, measurement on countermeasure

1/2

SSS25-12

会場:A04

時間:5 月 25 日 12:00-12:15

2/2

SSS25-13

会場:A04

時間:5 月 25 日 12:15-12:30

波相論と液状化
Wave features Theory and Liquefaction
西澤 勝 1∗NISHIZAWA, Masaru1∗
1
1

なし
none

1．まえがき
図-1 に従って液状化した場合の波形について、詳しく述べる。日本では昨年地震学会が、新潟地震 50 年ということで
行われた。
2．要旨
NZ 地震で倒壊した CTV ビルに近い、クライスチャーチにおける地震波相（Wave Features Theory）を図-1 に示してある。右側に主として、CTV ビルにとって大切なことを、二つ示した。つまり、地盤が液状化すると、(1) 地震波の位相がずれるため、建物に回転力というか、ねじれが生ずる。そのねじれは、東西南北はむろん、上下方向にも生じて建物に影響を与える。(2) これは CTV ビルだけでなく、クライスチャーチの平野 (?) 一帯が液状化（日本の浦安と同様）しているため、ここ全体の液状化対策が必要となる。以下図-1 に従って説明する。
図-1 の東西、南北の水平波相（Wave features）が軟弱地盤（Soft Ground）の波形。振動数も主要動以後のそれも異なっている。一番下の鉛直方向の波相が実際の地震波（Seismic Wave）に最も近い。このことから、この軟弱地盤は、そんなに厚くないことが判る。主要動（7〜8 秒）以後の波相（Wave features）も波打っていることから、他の地震計（この地区だけで 10 個程の地震計が設置されている）これら一帯が全体として波打っている。つまり、東京浦安も同様であるが、これら一帯が広く液状化し、波打ったと考えられる。全体を波動理論（Wave theory）の考えで扱う必要がある。一番下の鉛直方向の波相（Wave features）の中心線の上下で、下の揺れの振幅が上側の振幅より、振幅が小さいのが目立つ。他の地震計でもかなり認められる。これは軟弱地盤の下によりかたい地層があるためと考えられる。
つまり、軟弱な地層の下に硬い地層があるため、波が反射を受け、中央線上部より下部の振幅が小さくなる。これは、地震波だけでなく、日常の生活でもよく見られる現象である。トランポリンのようなものも一種のこの考えで良いと考える。
この、クライスチャーチの波形は日本では、1964 年の新潟地震の川岸町県営アパート 3 号棟の地下室の強震計で見事にとらえれている。図-2 にそれを示した。図-3 は 30 秒間のランニング・スペクトルで、最初と最後では、ずいぶん異なっていることが判る。
なお、新潟地震の波形というランニング・スペクトルの図は参考文献にも示したが、大崎順彦先生の著書「新・地震動のスペクトル解析入門」より転用させていただいた。新潟平野は大きな地震が来ると、NZ や浦安のように広域な液状化をおこし易いと考える。
参考文献(1) 大崎順彦 著：新・地震動のスペクトル解析入門(2)MASARU NISHIZAWA：The Great Kanto Earthquake (of 1923) and TOKOHAMA(1). Japan Geoscience Union Meeting
2014

Keywords: Wave Features Theory, Liquefaction, Phase, Rolation

1/2

SSS25-13

会場:A04

時間:5 月 25 日 12:15-12:30

2/2

SSS25-14

会場:A04

時間:5 月 25 日 12:30-12:45

海溝型プレート内地震のための地震動予測式の補正項に関する検討
An additional correction term of ground-motion prediction equation for intra-plate earthquakes
森川 信之 1∗ ; 藤原 広行 1
MORIKAWA, Nobuyuki1∗ ; FUJIWARA, Hiroyuki1
1
1

防災科学技術研究所
NIED

沈み込んだプレート内で発生するスラブ内地震は、同規模のプレート間地震と比べて短周期地震動が強く励起され、その強さには震源深さ依存性があること（例えば、Asano et al., 2003）やプレートによる違いがあること（例えば、池田・他、2004）が指摘されている。一方で、海溝軸より沖側で発生する、いわゆる「アウターライズ地震」は沈み込むプレー
ト内の浅い地震であるものの、その短周期地震動の強さは震源が深いスラブ内と同程度であるとする検討結果もある（例
えば、川端・他、2010）。
海溝型プレート内地震を対象とした強震動予測において、短周期地震動の励起特性を震源の深さやプレートの違いを踏まえて適切にモデル化することが重要である。我々は新しい地震動予測式を提案している（Morikawa and Fujiwara, 2013；
以下、MF13 式）が、海溝型プレート内地震に対する地震動強さの評価に関して、伝播経路特性としての異常震域の補正項を既に求めているものの、上述のような地震動強さの震源深さ依存性やプレートの違いは考慮されていない。そこで、
MF13 式に付与する補正項を求めるという観点から強震動記録の分析を行う。
はじめに、地震動強さの最大振幅について MF13 式による予測値との残差を求め、続いて地震ごとの残差の平均値（震源特性値）と観測点ごとの残差の平均値（観測点特性値）を求める。このとき、予測値は異常震域の補正も含めて求める。このようにして得られた各地震の「震源特性値」について震源の深さやプレートの違いに関する比較を行う。検討には、MF13 式を求めるために用いた強震動データベースの中から、5 地点以上地表で観測されている Mw5.0 以上の地震の強震動記録を用いる。
フィリピン海プレートの地震と比べて、太平洋プレートの地震の「震源特性値」が大きい傾向が見られる一方で、震源が深くなるほど「震源特性値」が大きくなるという傾向も見られる。ただし、現状のデータにおいてばらつきが大きいことに加え、太平洋プレートは震源が深い（50km 以深）地震が多く、フィリピン海プレートは震源が浅い（50km 以浅）地震が多いため、両者の影響を分離することは困難である。従って、予測式の補正項としては当面どちらかの影響のみをモデル化し、不確実さ（ばらつき）を別途考慮することで対応せざるをえないと考えられる。
参考文献：
Asano, K., T. Iwata and K. Irikura (2003): Earth Planets and Space, 55, pp. e5-e8.
池田孝・武村雅之・加藤研一（2004）：日本建築学会構造系論文集，586，pp. 53-61．
川端渉・笹谷努・高井伸雄・前田宜浩（2010）：第 13 回日本地震工学シンポジウム論文集， pp. 275-280．
Morikawa, N. and H. Fujiwara (2013): Journal of Disaster Research, 8, pp. 878-888.

キーワード: スラブ内地震, アウターライズ地震, 地震動予測式
Keywords: intra-slab earthquakes, outer-rise earthquakes, ground-motion prediction equation

1/1

SSS25-15

会場:A04

時間:5 月 25 日 14:45-15:00

中国・四国・九州地方におけ地震動の距離減衰分布
Attenuation characteristics of strong ground motions in Chugoku, Shikoku and Kyushu
distrints
池浦 友則 1∗
IKEURA, Tomonori1∗
1
1

鹿島技術研究所
Kajima Tech. Res. Inst.

■はじめに
西南日本の中国・四国・九州地方における強震動の距離減衰特性を検討する．この地域はフィリピン海プレートが屈曲して沈み込む領域にあたり，周辺には大規模な火山帯が存在するなど，構造的な条件がやや複雑である．また，東北
日本に比べて地震活動度が低いため，強震動の距離減衰特性を正しく把握するのがやや難しい地域でもある．とはいえ，1997 年 3 月 26 日と 5 月 13 日の鹿児島県北西部地震，1997 年山口県北部地震，2000 年鳥取県西部地震，2005 年福岡県西方沖地震と内陸の被害地震がいくつかあり，2001 年芸予地震のようなスラブ内の被害地震も発生している．これらの地震による強震動の距離減衰特性を正しく把握しておくことは有意義であろう．
■西南日本地域における K-NET・KiK-net のサイトファクター評価
中国・四国・九州地方における K-NET・KiK-net 地点で隣接観測点のネットワークを構成し，最小二乗法を用いてす
べての隣接観測点間の平均スペクトル比を同時に満足するように，参照観測点に対する相対サイトファクターを評価した．隣接観測点としては地点間距離 25km 以内を目安とし，各隣接 2 地点間の平均スペクトル比を求める記録は，➀震源から 2 観測地点に伸ばした直線の角が 5 度以内，➁両地点の PGA が 100Gal 以下，➂両地点の震源距離が 300km 以内，の3 条件を満足することとしている．
最終的に中国・四国・九州地方の K-NET と KiK-net で隣接観測点 2173 ペアを用い，630 点（地表と地中は別地点扱い）の相対サイトファクターを評価した．最短の隣接 2 地点間距離は 0.1km，隣接 2 地点間で用いた地震数はペア毎に大きく異なり 4〜108 地震，最大規模の地震は 2004 年 9 月 5 日 23 時 57 分の紀伊半島沖地震 (M7.4) である．参照観測点
としては YMGH12 美東の地中 (GL-102m,Vp=6000m/s, Vs=3250m/s) を採用した．なお，九州地方北西部では，隣接する観測点ペアで上記の条件を満足する地震が乏しく，隣接観測点間の相対サイトファクターが安定的に評価できないため除外した地点も多い．
■内陸地震とスラブ内地震の距離減衰特性の検討
各 K-NET・KiK-net 地点の観測地震動スペクトルを，上記で得られた相対サイトファクターで割算して距離減衰特性を検討した．2000 年 10 月 6 日鳥取県西部地震 (M7.3,h9km) による高周波数領域の距離減衰特性は地点毎の揺れ易さの違いが排除されて非常によくまとまり，両対数図上でみると震源距離 10〜500km の範囲でほぼ線状の分布となった．ただし，震源付近から震源距離 200km 前後までの区間とそれ以遠では明らかに傾斜が異なり，遠方の傾斜が急である．震
源距離 200km は概ね九州地方にさしかかる距離であり，この折れ曲がりは九州地方に入って以降の高減衰を示唆する．
2001 年 3 月 24 日芸予地震 (M6.7,h46km) による距離減衰分布も，周波数 10Hz 程度以下の範囲では 2000 年鳥取県西部地震ほどではないが，ある程度よくまとまり，震源距離 50〜400km の範囲で比較的単純な距離減衰分布の特徴を示した．ただし，20Hz 程度の高周波数になると，距離減衰効果に地点によるばらつきが認められ，早めに減衰する地点とそ
うでない地点との違いがやや明瞭になった．
さらに，この特徴は 2014 年 3 月 14 日伊予灘のスラブ内地震 (M6.2,h78km) でより顕著に観察された．この地震で震源
距離 80〜350km の範囲の距離減衰分布を検討したところ，観測地震動の場合は高周波数になるほどダンゴ状の分布であるが，相対サイトファクターで割算することにより，いくつか枝状の距離減衰分布が浮かび上がってくることが分かった．枝状の距離減衰分布の様子は周波数帯域によって変化し，2Hz 付近の振幅は中国＞四国＞九州，一方，10Hz 以上の振幅は四国＞中国＞九州であった．また，20Hz 付近では forearc ＞ backarc という単純な構図も見える．
最近の伊予灘のスラブ内地震としては他に 1998 年 5 月 23 日 (M5.4,h86km)，2006 年 9 月 26 日 (M5.3,h70km) があり，
これらでも類似の特徴が見いだされた．そのため，このような複雑な距離減衰特性は，この地域特有の現象と考えることができ，筧 (2014) が指摘しているとおり，沈み込んでいるフィリピン海プレートの形状や高減衰の火山帯の分布といった構造的な要因を踏まえて定量的に解釈される必要がある．
■参考文献：池浦・加藤 (2011) 日本地震工学会論文集，筧 (2014) 日本地震学会 2014 年秋季大会■謝辞：防災科研 K-NET・KiK-net の記録を使用させていただきました．記して感謝します．

キーワード: 強震動, 距離減衰特性, サイトファクター, 西南日本
Keywords: strong motions, attenuation characteristics, site factor, southwest Japan

1/1

SSS25-16

会場:A04

時間:5 月 25 日 15:00-15:15

周波数・震源距離依存型のラディエーション係数のモデル化に関する検討
Study on frequency and hypocentral distance dependent radiation coefficient
長坂 陽介 1∗ ; 野津 厚 1
NAGASAKA, Yosuke1∗ ; NOZU, Atsushi1
1
1

独立行政法人 港湾空港技術研究所
Port and Airport Research Institute

強震動シミュレーションにおいて，揺れの向きは震源メカニズムに応じたラディエーション係数により考慮することができる．しかし，理論的なラディエーション係数 (Aki and Richards, 2002) を用いると，観測よりも極端な方向性が出てしまう．実際には地中での散乱などによりある程度平均的（＝等方的）な値に近づくと考えられるため，適切にラディ
エーション係数をモデル化することができれば強震動シミュレーションの精度の向上が期待される．
ラディエーション係数のモデル化については多くの既往研究がある．Liu and Helmberger (1985) はラディエーション係
数には周波数依存性が存在し，高周波数ほど平均的な傾向になることを指摘している．強震動シミュレーションのためのモデルの提案としては，例えば釜江ら (1990), Pitarka et al. (2000), Pulido (2002) などがある．これらはいずれも Boore
and Boatwright(1984) の方法により理論的ラディエーション係数を空間的に平均する方法を採用しているが，高周波数になるほど平均をとる範囲を広げることで周波数依存性を考慮している．また，佐藤 (2002) は 2000 年鳥取県西部地震とその余震記録の分析を行い、経験的モデルを提案している。このとき、観測フーリエスペクトルの水平 2 成分の比をラディ
エーション係数の水平 2 成分の比とみなしており、また、周波数だけでなく震源距離にも依存したモデルとなっている．
いずれの既往研究においても，ラディエーション係数が理論的から平均的に遷移する周波数帯は，たとえば Pitarka et al.
(2000) では 1˜3Hz など，構造物に対する影響が大きい周波数帯となっている．
著者らは観測記録に見られるラディエーション係数は理論的な値と平均的な値の重みづけ平均として表されると仮定
し，その重みづけ係数（α とする）が周波数や震源距離によってどのように変化するかを観測記録に基づいて調べてきた． 2000 年鳥取県西部地震の観測記録の分析により α は周波数・震源距離双方に依存することが分かり，最終的には
α=exp(-πfr/QR VS ) としてモデル化することを提案している．この式は震源から対象地点までの波数が多いほど平均的なラディエーション係数に近づくことを表しており，QR は周波数・震源距離への依存度を決定する係数となっている．
2000 年鳥取県西部地震は観測記録が豊富であり，また震源過程も詳細に研究されている地震であるという利点がある
が，震源過程の複雑さから理論的ラディエーション係数を正確に知ることが困難である．そのため，今後は点震源とみなせるような地震についても提案モデルの適用性と係数 QR の値について検討を行い，報告する予定である．
キーワード: ラディエーション係数, 強震動シミュレーション
Keywords: radiation coefficient, strong ground motion simulation

1/1

SSS25-17

会場:A04

時間:5 月 25 日 15:15-15:30

長周期地震動に対する構造および震源起因の付加体効果
Effects of Accretionary Prisms on Long-Period Ground Motions Associated with Velocity
Structure Models and Sources
郭 雨佳 1∗ ; 纐纈 一起 1 ; 三宅 弘恵 1
GUO, Yujia1∗ ; KOKETSU, Kazuki1 ; MIYAKE, Hiroe1
11

東京大学地震研究所
Earthquake Research Institute, University of Tokyo

南海トラフで発生する海溝型地震が、大阪・濃尾・関東平野などに長周期地震動をもたらすことが知られている。そのような長周期地震動の生成・伝播過程に影響を及ぼす要素の一つとして、南海トラフに沿って発達している付加体の存在が挙げられる。Yamada and Iwata (2005), Yoshimura et al. (2008), 渡辺・他 (2014) は、付加体は地震波の直達部の振幅を減衰させる一方で、波全体の継続時間を延ばす効果があることを指摘した。後藤・永野 (2013), 渡辺・加藤 (2013) は付加体効果における震源位置依存性を示した。一般に、陸側に比べて海域の速度構造モデルは高い不確定性を伴っている。特に、地震波に支配的な影響を及ぼす S 波速度構造に対する検討は一部の地域を除き、まだ十分になされていない。
このように、付加体効果の評価において、付加体モデルを含む海域の速度構造モデルの不確定性だけでなく、海溝型地震の震源の広域性や多様性を考慮に入れることが必要不可欠である。また、大阪・濃尾・関東の各平野の長周期地震動の卓越周期が異なるため、付加体効果を周波数ごとに議論することも重要である。
本研究ではまず、構造モデルの違いから生じる付加体効果の変化を明らかにするため、 2004 年 9 月 5 日 19 時 7 分頃
(JST) に発生した紀伊半島南東沖地震の前震 (Mw 7.2) を対象に三次元波動場シミュレーションを実施した。シミュレー
ションでは、以下に示す付加体の層構造が異なる 3 つの速度構造モデルを使用した。 (A) 全国１次地下構造モデル (Koketsu
et al. , 2008, 2012)、(B) 付加体の S 波速度に関する既往研究 (Takahashi et al. , 2002; Tsuji et al. , 2011, 2014) や J-SHIS 深部
地盤モデル (藤原・他, 2009, 2012) を参考に、モデル A の付加体領域内の層構造を変化させたモデル、(C) モデル A の付加体層の S 波速度を 3.2 km/s に置換し、付加体なしと仮定したモデル。これらのモデルによる計算結果に対して周波数ごとに比較検討を行った。波形計算にはボクセル有限要素法 (Ikegami et al. , 2008) を使用し、地形や海水および広帯域 Q
値一定の減衰を計算に取り入れた。計算の有効周波帯域は 0.05-0.3 Hz である。震源は Yamada and Iwata (2005) の点震源を使用し、震源時間関数は八木 (2004) を使用した。
次に、先行研究で示唆されている付加体効果の震源依存性に着目し、上記と同じ計算手法と速度構造モデルを用いて、震源の位置や破壊伝播効果を考慮した場合の付加体効果を調べた。震源として、内閣府 (2012) で想定されている四国沖
から東海沖にかけての強震動生成域の位置を参考に、複数の面震源を仮定した。
シミュレーションの結果、以下の知見が得られた。付加体ありの場合の直達波の振幅は、付加体なしの場合と比べて必ずしも小さくなるとは限らない。大阪・濃尾平野では、堆積層の厚い地点における最大振幅が付加体モデルに依存する傾向にあるが、地震波の伝播経路の大部分を付加体が占める関東平野ではこうした特徴は見られず、付加体ありの場
合の最大振幅がおおむね小さい。付加体東縁における長周期地震動の伝播特性はモデル A と B の間で異なり、モデル Bの場合、伊豆半島や関東平野において周期 9-12 秒の後続波がより発達する可能性がある。また、付加体の後続波に対する増幅効果（継続時間伸長効果）は、トラフ寄りの浅い震源の破壊伝播方向側で顕著である。
キーワード: 長周期地震動, 付加体, 南海トラフ, 海溝型地震, 地震動シミュレーション, 速度構造モデル
Keywords: Long-period ground motion, Accretionary prism, Nankai Trough, Subduction earthquake, Ground motion simulation,
Velocity structure model

1/1

SSS25-18

会場:A04

時間:5 月 25 日 15:30-15:45

南海トラフにおけるやや長周期地震動の観測とシミュレーションによる再現
Long-period ground-motion observations and simulations in the Nankai Trough, southwest Japan中村 武史 1∗ ; 竹中 博士 2 ; 岡元 太郎 3 ; 大堀 道広 4 ; 坪井 誠司 1
NAKAMURA, Takeshi1∗ ; TAKENAKA, Hiroshi2 ; OKAMOTO, Taro3 ; OHORI, Michihiro4 ; TSUBOI, Seiji11

海洋研究開発機構, 2 岡山大学, 3 東京工業大学, 4 福井大学
Japan Agency for Marine-Earth Science and Technology, 2 Okayama University, 3 Tokyo Institute of Technology, 4 University of
Fukui

1

海洋研究開発機構では、紀伊半島沖の水深 1,900-4,400 m の海底に 20 点から成る地震観測網（DONET）を 2010 年に設置した。各海底観測点には強震計を備えており、リアルタイムで気象庁や防災科学技術研究所などの機関にデータ配信を行っている。南海トラフをはじめとする観測点周辺の地域で発生する大地震に対し、海底における強震観測データ
は、早期地震警報のためのリアルタイム解析や断層震源解析への活用が期待されている。
本研究では、海底強震観測データを用いて、2013 年 4 月に淡路島で発生した地震（Mw 5.8）時の海底におけるやや長周期成分の震動の解析を行った。陸上の強震観測データ（K-NET）と比較を行ったところ、4 秒以上の帯域において、海底観測点では最大振幅値が数倍大きい値を示すことが分かった。一方、1 秒以下の短周期成分では振幅値に大きな違いを見出すことはできなかった。震動の継続時間については、陸上では数 10 秒程度の継続時間であった。一方、海底では減衰が小さく、100 秒以上にわたって震動が継続していることが分かった。震源及び海底観測点周辺の領域における 3 次元
波動場のシミュレーションを試みたところ、10-20 秒の帯域において、振幅の増幅や震動長大化の特徴について観測デー
タと整合的な結果を得ることができた。シミュレーションから得られたスナップショットでは、表面波が付加体へ入射
後、振幅の増幅や波群の伝播速度が低下する様子が確認できた。また、表面波の波線追跡でも同様に、付加体で表面波群の伝播速度が急激に低下することを確認した。これらは、陸上の堆積平野部でこれまで多く観測された堆積層表面波と同様の特徴を示しており、海底観測点で見られたやや長周期成分の特徴的な波群は、低速度構造かつ盆地形状を成す
付加体内で発達した表面波の伝播によるものと考えられる。
やや長周期地震動の帯域（2-20 秒）は、早期地震警報における地震の規模推定や波形インバージョンなどでの断層震源解析時に重要な帯域である。本研究のシミュレーションでは、2 秒までの計算精度があるものの、やや長周期地震動の帯域の内、10 秒より短周期側の観測波形の再現はできていない。これは、ミュレーションで用いた S 波速度構造モデルと実際の構造との間に乖離があるためと考えている。今後は、コアロギングデータや変換波データなどを活用して、付加体における S 波構造の詳細な解析を行い、構造モデルを更新する必要があると考えている。

キーワード: 長周期地震動, やや長周期地震動, 強震動, 南海トラフ, 海底観測
Keywords: Long-period ground motion, strong motion, Nankai Trough, seafloor observation, DONET

1/1

SSS25-19

会場:A04

時間:5 月 25 日 15:45-16:00

理論的手法による長周期地震動評価の広帯域化に向けた検討
Approach to broaden the period-range of long-period ground motion evaluation based ontheoretical method
前田 宜浩 1∗ ; 岩城 麻子 1 ; 森川 信之 1 ; 今井 隆太 2 ; 青井 真 1 ; 藤原 広行 1
MAEDA, Takahiro1∗ ; IWAKI, Asako1 ; MORIKAWA, Nobuyuki1 ; IMAI, Ryuta2 ; AOI, Shin1 ;FUJIWARA, Hiroyuki1
1
1

防災科学技術研究所, 2 みずほ情報総研
NIED, 2 Mizuho Information & Research Institute, Inc.

はじめに
我々は、3 次元差分法による大規模地震動シミュレーションに基づいた海溝型巨大地震による長周期地震動のハザード評価に関する研究を行っている。これまで、不確実性を考慮して設定した多数の震源モデルや、堆積層やプレートなどが考慮された 3 次元地下構造モデルを用いた地震動シミュレーションにより、超高層建物や石油タンクなどに影響を及
ぼす周期 3〜20 秒程度の長周期地震動を対象とした検討を行ってきている。一方、多数存在するより小規模な構造物への影響を評価するためにはより短周期帯までの地震動評価が必要であった。そのためには、理論的手法による地震動計
算を短周期側へ拡張するアプローチと、統計的または経験的手法を取り入れた広帯域地震動計算手法の導入によるアプ
ローチ（岩城・他、本大会）が考えられる。本研究では前者のアプローチに関して、相模トラフの地震を対象として、これまでよりも精緻にモデル化された地下構造、震源モデルを用いて、長周期地震動評価の広帯域化に向けた検討を行う。
地下構造モデル
本研究では、防災科学技術研究所において南関東地域を対象として作成中の浅部・深部統合地盤モデル［例えば、Senna
et al.(2013、JDR)］を用いる。浅部・深部統合地盤モデルは、これまで別々にモデル化が行われてきた工学的基盤よりも浅い地盤モデル（浅部地盤構造）と，それよりも深い地盤モデル（深部地盤構造）を統合し、両方の地盤構造の影響を受ける周期帯（周期 0.5 秒から 2 秒）における増幅特性などをより適切に評価することで、強震動の予測精度の向上を目指
すものである。浅部・深部統合地盤モデルを用いて Vs=350m/s 以下の速度層（ここでは Vs=250m/s 層と仮定）までが考慮された計算を行い、これまでの Vs=500m/s の工学的基盤上での計算結果と比較したところ、低速度層がモデル化された地域では地震動の最大振幅や継続時間が増大する傾向がみられた。また、フーリエスペクトルの比較からは、両者の差異は主として周期 2 秒よりも短周期帯域において認められた。
震源モデル
これまでの長周期地震動シミュレーションでは特性化震源モデルを用いているため、アスペリティサイズよりも小さなスケールの不均質性が十分に考慮されていない。そこで、マルチスケールの不均質性（関口・吉見，2006，月間地球）を考慮した震源モデルを作成し、その効果について検討した。ここでは、相模トラフで想定される M8 級の地震に対して、破壊開始点やアスペリティ配置を変えた 274 個のシナリオについて、破壊伝播に対してそれぞれ異なる不均質性を付与した震源モデルを作成した。不均質性を付与した影響は、長周期帯域に比べ短周期帯域において大きくなる傾向がみられた。また、破壊開始点やアスペリティ配置に依存している傾向もみられた。
まとめ
これまで行ってきた長周期地震動シミュレーションを周期 1 秒程度まで広帯域化する上で、より精緻にモデル化された地下構造モデルや震源モデルを用いる重要性が示唆された。一方で、本研究では、低速度層までを考慮した周期 1 秒程度までの長周期地震動の計算を差分法により行っているが、これまでよりも格子サイズや時間刻みが小さくなっており、計算に要する時間が増大するとともに、多タイムステップの計算を安定して行うことが難しくなりつつある。長周期地震動評価の広帯域化を図るためには、地下構造モデルや震源モデルのみならず、適切な計算手法についての検討も重要な課題である。
謝辞 本研究は文部科学省の「平成 26 年度長周期地震動ハザードマップ等支援事業」による。

キーワード: 長周期地震動, 差分法, 地下構造モデル, マルチスケール不均質
Keywords: long-perido ground motion, FDM, velocity structure model, multi-scale heterogeneity

1/1

SSS25-20

会場:A04

時間:5 月 25 日 16:15-16:30

相模トラフの海溝型巨大地震を対象とした広帯域地震動計算手法の検討
Broadband ground motion simulation techniques applied to megathrust earthquakes in the
Sagami trough岩城 麻子 1∗ ; 前田 宜浩 1 ; 森川 信之 1 ; 藤原 広行 1
IWAKI, Asako1∗ ; MAEDA, Takahiro1 ; MORIKAWA, Nobuyuki1 ; FUJIWARA, Hiroyuki1
1
1

防災科学技術研究所
NIED

将来起こりうる海溝型巨大地震に対する地震動評価において、周期 1∼ 数秒程度以上の長周期地震動については、震源
破壊過程と 3 次元波動伝播過程を適切にモデル化した上で理論的手法によって評価することができると考えられる。岩
城・他 (2014) では、相模トラフで発生しうるマグニチュード 8 程度以上の海溝型巨大地震を対象として、震源パラメー
タの不確実性を考慮した多数の震源モデル（シナリオ）について三次元差分法による地震動シミュレーションに基づく
長周期地震動ハザード評価を行った。その際、震源モデルや地下構造モデルの精度を考慮して周期 3 秒以上の帯域を対象とした。
一方で、相模トラフの地震は震源域が首都圏直下まで広がることや、マグニチュード 7 程度の地震の発生も切迫しているとされていることから、これまでの長周期地震動ハザード評価を踏まえた上で、より短周期の地震動も含めた広帯域地震動ハザード評価を行うことが重要な課題となっている。そのためには、理論的手法による地震動計算を短周期側へ拡張するアプローチ（前田・他、本大会）と、統計的または経験的手法を取り入れた広帯域地震動計算手法の導入によるアプローチが考えられる。
本検討では、大正関東地震に相当する震源域を持つ数ケースの地震を対象として、広帯域地震動計算手法として従来広く用いられているハイブリッド法（先名・他, 2004）による計算を行う。また、岩城・藤原 (2013) などで提案された広帯域地震動計算手法の適用を試み、結果について計算手法同士の比較や理論手法による地震動との比較を行う。後者の手法は周波数帯域間の加速度エンベロープの比の経験的特性を利用して、長周期地震動波形から短周期へ広げるものである。ここでは便宜上 ERF (envelope ratio function) 法と呼ぶことにする。
どちらの手法についても、長周期側は全国一次地下構造モデル（暫定版；地震調査委員会, 2012）を基本としたモデル（以下 mod-JIVSM）を用いて三次元差分法（GMS；青井・他, 2004）で計算する。その際、数値計算上は周期 1 秒まで有効となるように計算を行うが、mod-JIVSM が周期 2 秒以上を対象としていることから、波形合成における接続周期は 2秒とする。短周期側は、ハイブリッド法では統計的グリーン関数法によって計算し、ERF 法では経験的な ERF と長周期側の波形を利用して短周期地震動波形を作成する。統計的グリーン関数法に用いる地下構造モデルは、周期 0.5-2 秒程度の地震動をより適切に評価することを目的として作成中の浅部深部統合地盤モデル（Senna et al., 2013）とする。
関東平野内の代表地点において、ハイブリッド法、ERF 法、および差分法による波形やフーリエスペクトルをそれぞれ比較した。ハイブリッド法と ERF 法は接続周期である 2 秒から 1 秒の間で差分法に比べてスペクトル振幅が有意に大きい。ハイブリッド法と ERF 法の速度波形における最初に到達する主要動部分の振幅レベルは互いにおおむね一致するが、後続動ではハイブリッド法の短周期成分が小さくなる傾向が見られ、特に振幅の大きい後続波群を持つ場合にハイブリッド法と ERF 法の違いが顕著であった。現行のハイブリッド法では、長周期と短周期の波形を重ね合わせるタイミ
ングに任意性があり、特に巨大地震では両者の波形の走時のずれが顕著に表れる可能性がある。
ここで検討した手法は、今後面的評価も含めた広帯域地震動ハザード評価に現実的な計算量で寄与できることが期待される。地震動予測式等との比較によってこれらの手法の適用性検証や手法同士の比較検討を行うことも重要である。

キーワード: 海溝型巨大地震, 長周期地震動, 広帯域地震動, 相模トラフ, 関東平野
Keywords: megathrust earthquake, long-period ground motion, broadband ground motion, Sagami trough, Kanto plain

1/1

SSS25-21

会場:A04

時間:5 月 25 日 16:30-16:45

弾塑性有限変形解析による断層的破壊に伴う不規則波生成シミュレーション
Simulation of irregular wave generation due to fault formation by an elasto-plastic finite
deformation analysis
山田 正太郎 1∗ ; 野田 利弘 1 ; 浅岡 顕 2
YAMADA, Shotaro1∗ ; NODA, Toshihiro1 ; ASAOKA, Akira2
1
1名古屋大学, 2 (公財) 地震予知総合研究振興会
Nagoya University, 2 Association for the Development of Earthquake Prediction

著者らはこれまでに，慣性力を考慮した水〜土連成弾塑性有限変形解析コードGEOASIA 1) を用いて横ずれ断層に伴う上部地盤のせん断帯の生成シミュレーションを行っている 2) ．本研究では，同解析コードを用いて，高い脆性を有する地盤を対象に，正断層および逆断層の生成と，それに伴う地震動の発生シミュレーションを実施した．この解析コードが搭載する構成式 SYS Cam-clay model3) は広範な土材料を同じ理論的枠組で記述する弾塑性構成式である．また，使用した解析コードは速度型の運動方程式を忠実に時間積分して運動を求める動的問題にも対応可能なコードであり，せん
断帯発生に伴う波動生成や波動伝播も再現可能 4),5) である．この特徴を活かし，今回は特に，断層破壊に伴う振動の発生に着目して解析を行った．平面ひずみ条件下で，解析断面側方から変位制御により地盤を圧縮したところ，ひずみの局所化を伴う進行性破壊現象として逆断層的な破壊が発生した（図 1）．この破壊は加速度的に進行し，それに伴って地盤内に振動が発生した．この際，非破壊領域に蓄えられた弾性エネルギーが一気に解放されることを確認した．水平成層地盤を対象にした解析では，破壊が急速に進行するため，初動で最大加速度に達し，その後指数関数的に減衰する人工地震と同様な特徴を持った振動が発生した（図 2）．一方で，ランダムな初期不整を与えた地盤では，大きな破壊イベ
ントの中に小さな破壊イベントが幾つも存在することに起因して，自然の地震波にも似た不規則な振動が発生した（図3）．なお，解析断面側方から変位制御により地盤を伸長した場合には，正断層的な破壊が生成され，やはりランダムな
初期不整を有する地盤では，不規則な振動が発生した．

1) Noda, T., Asaoka, A. and Nakano, M. (2008): Soil-water coupled finite deformation analysis based on a rate-type equationof motion incorporating the SYS Cam-clay model, Soils and Foundations , 48(6), 771-790.
2) 野田利弘，山田正太郎，浅岡顕，川合裕太 (2014): 横ずれ断層に伴う上部地盤のせん断帯発生の数値シミュレーショ
ン, 日本地球惑星科学連合大会 2014 年大会, SSS31-08.
3) Asaoka, A., Noda, T., Yamada, E., Kaneda, K. and Nakano, M. (2002): An elasto-plastic description of two distinct volume
change mechanisms of soils, Soils and Foundations , 42(5), 47-57.
4) Noda, T., Xu, B. and Asaoka, A. (2013): Acceleration generation due to strain localization of saturated clay specimen basedon dynamic soil-water coupled finite deformation analysis, Soils and Foundations , 53(5), 653-670.
5) 浅岡顕，山田正太郎，野田利弘 (2013): 地表面載荷による地盤の破壊現象と破壊に伴う振動の発生に関する数値解析,
日本地球惑星科学連合大会 2013 年大会, SSS28-18.
キーワード: 正断層, 逆断層, 地震波, ひずみの局所化, 慣性力, 弾塑性体
Keywords: natural fault, reverse fault, seismic wave, strain localization, inertial force, elasto-plastic body

1/1

SSS25-22

会場:A04

時間:5 月 25 日 16:45-17:00

Estimation of shallow shear wave velocity in Bandung basin, Indonesia using horizontal
to vertical (H/V) spectral ratio
Estimation of shallow shear wave velocity in Bandung basin, Indonesia using horizontal
to vertical (H/V) spectral ratio
Pramatadie Andi Muhamad1∗ ; 地元 孝輔 1 ; AFNIMAR Afnimar2 ; 山中 浩明 1
PRAMATADIE, Andi muhamad1∗ ; CHIMOTO, Kosuke1 ; AFNIMAR, Afnimar2 ; YAMANAKA, Hiroaki1
1
1

Tokyo Institute of Technology, 2 Bandung Institute of Technology
Tokyo Institute of Technology, 2 Bandung Institute of Technology

Bandung, the third populous city in Indonesia with population of around 2.7 million, is located in the western part of Javaisland. The city lies on a basin structure called Bandung basin. The concern of seismic risk in this area becomes importantas existing of 24 km-long Lembang fault in the northern part. As part of potential seismic hazard estimation, Microtremors
measurement is one of effective tool to estimate shear wave velocity profile and site amplification factor, especially in the urban
area. A 3 component microtremors measurement is conducted in 76 sites to observe horizontal-to-vertical (H/V) spectral ratio thatreflected to the ground characteristic. In the previous work, the obtained shear wave velocity models from Spatial Autocorrelation(SPAC) are used as reference model. Estimation of the ground structures are using the tuning factor to fitting the observed H/Vspectrum with the theoretical ellipticity of fundamental mode of Rayleigh wave from the velocity model. From the obtained
velocity model, we observed the thicknesses of soft layer (<500 m/s) from north to south of basin are changing from around 8
to 12 m in edge of basin area and around 40 m at the central of basin. The obtain profiles provide detail structure informationin Bandung basin. Also the calculation of AVS30s and site amplification factor were conducted in each site, to understand thepotential seismic hazard in the area. We also discuss a spatial variation of the amplification of earthquake ground motion using
the obtained profiles.
キーワード: H/V spectral ratio, S-wave velocity, shallow soil, Bandung basin, amplification
Keywords: H/V spectral ratio, S-wave velocity, shallow soil, Bandung basin, amplification

1/1

SSS25-23

会場:A04

時間:5 月 25 日 17:00-17:15

地震波干渉法結果のデコンボリューションに基づく位相情報の抽出―数値計算に基づく手法の検証と実記録への適用―
An evaluation to trace phase of surface wave using the seismic interferometry
元木 健太郎 1∗ ; 加藤 研一 1
MOTOKI, Kentaro1∗ ; KATO, Kenichi1
1
1

小堀鐸二研究所
Kobori Research Complex

１．はじめに
地震波干渉法の適用例として、長時間の微動記録から 2 地点間の群速度が評価されている。地震波干渉法によってグリーン関数が再現できるとすれば、表面波の波群の情報だけでなく、位相情報も有効に活用することができると考えられる。Behm et al.(2014) は地震波干渉法の結果による近傍観測点の位相の時間差は、群速度から求められる時間差よりもばらつきが小さいことを示しているが、実記録に基づくのみで検証はなされていない。本検討では、地震波干渉法によって得られる波の位相伝播を活用するための評価法を示し、解析的な検討から適用性を確認した後、実観測記録に適用した事例を示す。
２．評価法について
グリーン関数の位相を用いて地下構造モデルとの関連付けを行うためには、伝播する波を追跡できるように、波長に対して観測点間隔を充分に短くすることが望ましい。しかし、地震波干渉法によって評価できる波長範囲は、空間エリアジングが起きる短波長のみに限られる。そこで本検討では、対象地域の外の観測点に対して、対象地域内の観測点の相互相関関数を求める。対象地域の端の評価波形と、対象地域内の評価波形のデコンボリューションにより、対象地域
内の波の伝播を評価する。本検討では、上下成分を対象とした。
解析的検討では擬似的に equipartition を表現するために、周辺からの波を相反定理により有限差分法を用いて評価する。周辺の発振点は 40km × 90km のエリアに 0.2km 間隔に 90651 地点を配置する。4km の測線を対象地域として 0.4km
間隔で評価点（以降、評価点）を配置し、地震波干渉法を適用するための点は、測線の延長上に 18km 離れた点（以降、基準点）を測線の両側に配置した。モデルは単純化のために 2 層モデルとし、水平成層構造モデル、測線に対して直交する方向に不整形なモデル（2 次元不整形モデル）、堆積層の形状が盆地となるような 3 次元不整形モデルと、徐々に不整形性が強くなるモデルを設定した。基準点を加振した時の各評価点の差分法の結果をグリーン関数として、地震波干渉法による結果をデコンボリューションして波の伝播を評価する手法の適用性について検討した。
３．解析的検討による適用性の確認
水平成層構造モデルの場合、相互相関関数のデコンボリューション波形は、グリーン関数のデコンボリューション波
形に近く再現性は高い。2 次元不整形の場合、デコンボリューションをとる方向、すなわち波の伝播方向によって結果が異なった。硬い地盤上の評価波形に対するデコンボリューション波形においては、地震波干渉法の結果はグリーン関数とよく対応する。軟らかい層上の評価波形に対するデコンボリューション波形は、層境界に達するまでの波の伝播は地震波干渉法とグリーン関数でよく対応するものの、層境界で反射する位相については地震波干渉法とグリーン関数で若干の位相ずれが生じた。3 次元不整形の場合は、デコンボリューションをとる方向に加え、地震波干渉法を適用する際の波の伝播方向によって結果が異なる。硬い地盤から軟らかい地盤へ波が伝播することを評価する際は、地震波干渉法を適用する際の伝播方向と、デコンボリューションをとる方向を合わせると、地震波干渉法の結果はグリーン関数の結果とよく対応する結果が得られた。しかし、やわらかい層から硬い層へ伝播する方向に対しては、解析の方向を合わせても地震波干渉法の結果とグリーン関数の結果は対応せず、グリーン関数を再現することはできなかった。硬い地盤から軟らかい地盤への波の伝播であれば、解析する処理の方向に注意を払えば、3 次元不整形地盤においても地震波干渉法によって波の伝播を評価できることを示した。
４．実データへの適用例
実記録の検討では、Hi-net 一関西 (N.ICWH) と一関東 (N.ICEH) 間を対象地域として、その 2 点の延長線上近くに位置している藤沢（N.FSWH）と東和（N.TOWH）を、地震波干渉法を適用するための基準点として選定した。1 か月間の記録を用いて、地震波干渉法を適用した。記録は計器特性を補正し、データ処理方法は地元・山中（2012）の方法を参照
した。記録から求めた波形と J-SHIS の 3 次元地下構造モデルを用いた解析解の比較を行った。
藤沢を基準として地震波干渉法の結果をデコンボリューションした波形と、東和を基準として処理した結果の波形はよく対応しており、基準点によらず安定して波形が求められることを確認した。3 次元地下構造モデルによる解析解と、観測記録から得られた波形は概ね対応する結果が得られ、実観測記録でも適用できることを示した。
今後の課題として、例えば関東平野の盆地端部付近の観測点に適用して表面波の生成が本手法で追跡できるかどうかの検討を行いたい。

1/2

SSS25-23

会場:A04

時間:5 月 25 日 17:00-17:15

キーワード: 地震波干渉法, デコンボリューション
Keywords: seismic interferometry, deconvolution

2/2

SSS25-P01

会場:コンベンションホール

時間:5 月 25 日 18:15-19:30

震度分布からみた明治三陸地震津波の強震動の特徴
Characteristic of Strong Motion in Sanriku Tsunami of 1896 by the point of view ofseismic intensity distribution
相澤 幸治 1∗AIZAWA, Koji1∗
11

気象庁地震火山部管理課
Administration Division, Seismology and Volcanology Department

明治 29 年（1896 年）6 月 15 日に発生した明治三陸地震津波は、表面波などの地震動の大きさから推定される地震の規模に比べ津波の規模が大きい「津波地震」であると考えられている（Kanamori，1972；阿部，1989）。
明治三陸地震津波の震度分布は２つの報告がある。中央気象台（1896）の等震度分布図では、東北北部を中心に弱震、関東南部から北海道にかけての広い範囲を微震としている。大森（1901）は、等震度線と観測した震度とを示し、等
震度線は弱震部と微震部に分けて中央気象台（1896）と同様の分布を示しているが、弱震部に強震が報告された地点を示している。今回、気象庁本庁に保管されている地震調査原簿に記載されている地点の震度と大森（1901）の図中の地
点の震度を対照すると、数点の相違を除き同じことが分かった。このため大森（1901）を基本として、観測された震度階の分布について検討を行った。東北地方北部では、北上山地では微震が、北上河谷帯、小川原低地および脊梁山地西側に弱震、強震の地点が多い。震度は震源からの距離による減衰が明瞭ではない。また、震度分布と深部地盤の一次固有周期（中央防災会議、2008；横田ほか、2010）を重ね合わせると、強または弱を観測した地点は深部地盤構造の一次固有周期が２〜３秒以上の地域に分布し、微の地点は１〜２秒以下の地域に多い。
地震調査原簿や大森（1901）は、明治三陸地震津波の揺れの性質について「緩慢」「緩」と表現している。野田ほか（2000）は、人間は、振動の様子について、周期 1.5 秒以上の揺れについては「ゆっくり」と表現するとしている。このため、明治三陸地震津波は少なくとも周期 1.5 秒以上の周期成分が卓越していたと考えられる。東北地方北部で震度階が大きな地点は平成 23 年（2011 年）東北地方太平洋沖地震における長周期地震動階級（相澤ほか、2014）が相対に大きい
地点と類似する。このため、長周期地震動階級の計算に用いる絶対速度応答スペクトルについて、渡辺（1998）の震央位置を用い、Dhakal（2013）、気象庁（2014）の手法を用いて、マグニチュードを変化させ、周期毎に距離減衰特性を調べた。その結果、マグニチュードが 7.6 以上、周期５秒以上場合に、絶対速度応答スペクトルと震度階の距離減衰の特徴が近いことが分かった。平成 23 年（2011 年）東北地方太平洋沖地震の強震波形を周期５秒でカットオフして計測震度を計算すると、東北地方北部では概ね震度２から３となり、北上山地に比べて脊梁山地西側で震度が大きくなる。
以上の結果から、明治三陸地震津波における震度は、観測された地震動の性質と合わせ、周期５秒程度以上の長周期地震動による体感を示していると考えられる。

キーワード: 震度, 長周期地震動, 津波地震, 強震動
Keywords: Seismic Intensity, Long-term Period Ground Motion, Tsunami Earthquake, Strong Motion

1/1

SSS25-P02

会場:コンベンションホール

時間:5 月 25 日 18:15-19:30

大規模シミュレーションによる沈み込み帯地震からの短周期表面波の再現性検討
Reproducing the short-period seismic surface waves from subduction zone earthquakes
based on large-scale simulation
岡元 太郎 1∗ ; 竹中 博士 2 ; 原 辰彦 3 ; 中村 武史 4 ; 青木 尊之 1
OKAMOTO, Taro1∗ ; TAKENAKA, Hiroshi2 ; HARA, Tatsuhiko3 ; NAKAMURA, Takeshi4 ; AOKI, Takayuki1
1
1

東京工業大学, 2 岡山大学, 3 建築研究所, 4 海洋研究開発機構
Tokyo Institute of Technology, 2 Okayama University, 3 Building Research Institute, 4 JAMSTEC

The seismic-wave propagation from shallow subduction zone earthquakes is strongly affected by the heterogeneous structure:
there are the thick layers with low-seismic wave velocities, such as the oceanic water and the sedimentary layers, and these layers
are separated with three-dimensionally irregular interfaces. These heterogeneities affect both the excitation and propagation of
the seismic waves. Therefore they must be considered in generating the synthetic waveforms for the waveform analysis of theearthquake sources and structural heterogeneities, otherwise the analysis could result in erroneous solutions [1,2]. In this paper westudy the effect of the structure of the Japan trench on the strong-motion seismic waveforms, especially the surface waves, froma shallow subduction zone earthquake in connection with the analysis of the 2011 Tohoku-Oki earthquake (Mw9.1). We selectedan small event (2003/11/1, Mw5.8) because this event occurred in the source area of the 2011 Tohoku-Oki earthquake, and has
a similar mechanism (inter-plate thrusting) as that of the 2011 Tohoku-Oki earthquake. The strong-motion records (K-NET,
KiK-net) were integrated to obtain velocity components. We assumed a three-dimensional structure model for the Japan trenchby compiling models for topography [3], sediments [4], crust, and subducting plates [5,6]. We applied a GPU-accelerated finitedifference program developed by ourselves [7,8]. We used the TSUBAME-2.5 supercomputer in Tokyo Institute of Technology
for the finite-difference computation. By using the 3D structure model, the strong-motion seismograms were well reproduced fora period band of 12-40 s. For periods shorter than around 10 s, however, the misfit of surface waves was large: the amplitudes ofthe synthetic surface waves were smaller than that of the observations. That is, the excitations of the short period surface waves
from shallow earthquakes were not sufficient in the assumed 3D structure model. As an experiment we reduced the S-wave
velocities of the oceanic sediments by 30 % and re-computed the synthetics. However, the fit between the observed and synthetic
waveforms were not improved. We will discuss these results and present findings based on further experiments (e.g., results bymodifying the thickness of the sediments).
References[1] Okamoto, T. and Takenaka, H., Advances in Geosciences , 13 (Solid Earth), p.215-p.229, World Scientific Pubishing, 2009.
[2] Okamoto, T. and Takenaka, H., Earth Planets Space , 61, e17-e20, 2009.
[3] Kisimoto, K., Geological Survey of Japan, Open-file Report , No. 353, 2000.
[4] Fujiwara, H. et al., Technical Note of the National Research Institute for Earth Science and Disaster Prevention , No.337,
2009.
[5] Baba, T. et al., Abstr. Japan Geoscience Union Meet. , S111-006, Makuhari, Japan, 2006.
[6] Nakamura, T. et al., Abstr. Seism. Soc. Japan , P1-06, Hiroshima, Japan, 2006.
[7] Okamoto, T. et al., Earth Planets Space , 62, 939-942, 2010.
[8] Okamoto, T. et al., GPU Solutions to Multi-scale Problems in Science and Engineering , 375-389, Springer-Verlag, 2013.
Keywords: subduction zone earthquake, waveform modeling, short period wave, GPU computing

1/1

SSS25-P03

会場:コンベンションホール

時間:5 月 25 日 18:15-19:30「全国地震動予測地図 2014 年版〜全国の地震動ハザードを概観して〜」 の公表についてOn National Seismic Hazard Maps of Japan 2014 edition
地震調査研究推進本部 事務局 1∗
JISHINHONBU, Jimukyoku1∗
1
1

文部科学省
Ministry of Education, Culture, Sports, Science and Technology

地震調査研究推進本部地震調査委員会は、活断層で発生する地震と海溝型地震の長期的な発生確率を評価するとともに、いくつかの震源断層を対象に強震動を予測し、公表してきた。それらに基づき、2005 年 3 月に「全国を概観した地
震動予測地図」を公表し、以来、毎年評価の改訂を行いその結果を公表してきた。2009 年 7 月には全面的な改訂を行い、名称を「全国地震動予測地図」に変更した。
2011 年も全国地震動予測地図 2011 年版を公表する予定であったが、平成 23 年（2011 年）東北地方太平洋沖地震が発生し、確率論的地震動予測地図 磯 について解決すべき多くの課題が指摘されたことなどにより、公表を見送り、その作成手法の基本的な枠組みの有効性を確かめるとともに、指摘された課題の検討を開始した。それらの検討を踏まえて作成したモデルに基づいた確率論的地震動予測地図は、「今後の地震動ハザード評価に関する検討〜2011 年・2012 年における検討結果〜」、「今後の地震動ハザード評価に関する検討〜2013 年における検討結果〜」として、それぞれ 2012 年、2013 年に公表した。
2014 年も引き続き課題の検討を行うとともに、モデルの作成を行ってきた。このたび、これまでに行った検討結果を踏まえた確率論的地震動予測地図の作成が完了したため、「全国地震動予測地図 2014 年版〜全国の地震動ハザードを概観して〜」（以降、2014 年版）を公表する。

※日本とその周辺で発生する地震の位置・規模・発生間隔などに基づいて、各地点がどの程度の確率でどの程度揺れ
るのかなどを計算し、地図に示したもの。
キーワード: 全国地震動予測値図２０１４年版, 地震動ハザード, 地震調査委員会
Keywords: National Seismic Hazard Maps of Japan 2014 edition, seismic hazard, Earthquake research committee

1/1

SSS25-P04

会場:コンベンションホール

時間:5 月 25 日 18:15-19:30

関東堆積盆地の地震波速度構造モデルと長周期地震動シミュレーション
Simulation of long-period ground motion using 3D structure model of the Kanto Basin
吉本 和生 1∗ ; 武村 俊介 1 ; 鍛治川 謙吾 1 ; 増田 啓 1 ; 玉置 大志 1 ; 北澤 楽奈 1 ; 小林 憩加 1
YOSHIMOTO, Kazuo1∗ ; TAKEMURA, Shunsuke1 ; KAJIKAWA, Kengo1 ; MASUDA, Kei1 ; TAMAKI, Taishi1 ;KITAZAWA, Rana1 ; KOBAYASHI, Keika1
11

横浜市立大学
Yokohama City University

はじめに
首都圏における長周期地震動の正確な評価のためには，関東堆積盆地の精緻な地震波速度構造モデルが求められる．
本研究では，近年構築された関東堆積盆地の地震波速度構造モデル（例えば，SBVSM: 増田・他, 2014; Takemura et al.,2015）を用いて中規模地震についての長周期地震動シミュレーションを実施し，主に関東堆積盆地の北部及び西部で観測された長周期地震動（周期 4-8 秒）の再現性について検討した．
差分法による長周期地震動シミュレーション
2004 年 10 月 27 日に新潟県中越地方で発生した地震（深さ 11 km，Mw 5.8）と 2013 年 2 月 25 日に栃木県北部で発生した地震（深さ 8 km，Mw 5.8）について，3 次元差分法による長周期地震動シミュレーションを行った．計算領域は，震源域及び関東平野北部を含む 201.6×127.5×60 km3 の領域で，水平方向に 0.15 km，鉛直方向に 0.075 km で離散化した．
発震機構解には F-net による CMT 解を採用し，震源には点震源モデルを仮定した．地震モーメントの値は岩盤点における観測地震動の振幅を利用して調整した．長周期地震動シミュレーションは，SBVSM と JIVSM（Koketsu et al., 2008）に加えて，SBVSM に堆積盆地端部における局所的な S 波速度構造の推定結果（武村・他, 2014）を取り入れたモデル（以下，SBVSM2）についても実施し，観測された長周期地震動の再現性を比較した．SBVSM 及び SBVSM2 の堆積層以深の構造については，JIVSM による地震基盤形状，Moho 面形状，海洋プレートの上面形状を仮定した．
長周期地震動シミュレーションの結果
SBVSM を用いた計算結果は，JIVSM による計算結果と比較して，堆積盆地端部の一部の観測点を除いて観測された長周期地震動の振幅や発現時刻などの特徴をより良く再現した．さらに，堆積盆地端部の地震波速度構造を改善した SBVSM2では，関東堆積盆地の北端部で観測された長周期地震動の再現性に一段の向上が見られた．この結果は，武村・他（2014）の手法による堆積盆地端部の地震波速度構造モデルの精緻化により，堆積盆地端部での表面波の励起・伝播特性を実用的に正確に評価できる可能性を示している．その一方で，堆積盆地端部の地震波速度構造モデルが精緻化されていない
関東堆積盆地西縁（埼玉県西部）における長周期地震動については，観測記録の再現性があまり高くないことが確認された．学会発表では，SBVSM 及び SBVSM2 の特徴と現状の問題点について詳しく紹介するとともに，堆積盆地端部における地震波速度構造と表面波の励起・伝播特性の関係についても説明する．
謝辞
本研究では防災科学研究所の K-NET/KiK-net 強震記録ならびに F-net の CMT 解を利用させていただきました．また，
首都圏強震動総合ネットワーク SK-net の強震記録を使用しました．地震動シミュレーションには東京大学地震研究所の計算機システムを利用させていただきました．記して感謝いたします．
キーワード: 長周期地震動, 関東堆積盆地, 地震動シミュレーション
Keywords: long-period ground motion, Kanto Basin, simulation of seismic ground motion

1/1

SSS25-P05

会場:コンベンションホール

時間:5 月 25 日 18:15-19:30

首都圏の地震動応答継続時間の面的分布に関する検討
Study on Spatial Distribution of Response Duration Time of Earthquake Motions in Tokyo
Metropolitan Area
新井 健介 1∗ ; 石井 透 1 ; 平田 直 2
ARAI, Kensuke1∗ ; ISHII, Toru1 ; HIRATA, Naoshi2
1
1

清水建設, 2 東京大学
Shimizu Corporation, 2 University of Tokyo

１．背景と目的
石井（2012）は、地震動の継続時間の周期特性に注目し、これを定量評価するための新たな地震動評価指標として応答継続時間スペクトルを定義・提案した上で、東北地方太平洋沖地震とその前震・余震・誘発地震の際に東京都江東区
越中島で得られた複数の観測記録のスペクトルを求めてその特性を定量評価し、検討した。特に周期 1 秒以上の長周期
地震動の継続時間は、短周期地震動や今まで考えられていた値よりも著しく長くなったり狭帯域でピークを持ったりしていることが明らかになった。また、震央方向の深い地下構造の違いが長周期表面波の伝播特性に違いをもたらし、主に観測記録の長周期成分に反映されたと推定した。特に長周期地震動の振幅・継続時間の特性には深い地下構造を反映した方位依存性が含まれると考え、多くの良質の観測記録に基づく検討の必要性を指摘した。
この課題に取り組むためには一地点での記録による検討では限界があり、面的な検討が必要となる。首都圏地震観測網（MeSO-net）では、首都圏の 300 以上の観測点において広帯域・長継続時間の記録が面的に得られている。そこで、新井・ほか（2014）は、MeSO-net による東北地方太平洋沖地震とその前震・最大余震を含む計 14 地震の広域・広帯域・
長継続時間地震観測記録を用い、応答継続時間スペクトルを算出し、その特性を検討した。引き続き本報では、特に長周期地震動の時刻歴・応答スペクトル・応答継続時間スペクトルの特性とその空間変動について、深い地下構造の変化との関係、震央の方位依存性や震源深さとの関係などについて検討する。
２．データと解析手法
本報では、新井・ほか（2014）が検討した首都圏地震観測網（MeSO-net）の記録の中から 2011 年 3 月 11 日から 2011年 4 月 12 日までの計 6 地震の記録の速度応答スペクトル SV [cm/s] および速度応答継続時間スペクトル TSV [s] を検討対象とする。この中には、2011 年東北地方太平洋沖地震とその最大余震の記録も含まれる。TSV のパラメータは、石井（2012）や新井・ほか（2014）と同じく、減衰定数 5%、パラメータ p1 = 0.03、p2 = 0.95 である。MeSO-net の観測点の中から複数の測線を選び、それに沿ったスペクトル特性の変化と地下構造の変化を対比させ、考察する。考察には、防災科学技術研究所の地震ハザードステーション J-SHIS の深い地下構造モデルを参照する。
３．結果
首都圏の地震動の応答継続時間スペクトルは、深い地下構造の影響により、空間的に変化している。応答継続時間ス
ペクトルの性状は、地震基盤以浅の構造が比較的単純な地点では特にやや長周期の狭帯域で比較的単純なピークを持つ場合が多いが、地震基盤以浅の構造が多層で複雑化するにつれて複数のピークを持つように見える。また、同一地点でも、地震規模や震源距離だけでなく、震央の方位や震源深さ等に応じて応答継続時間スペクトルの性状が異なり、伝播する表面波の性状を反映していると推定される。
謝辞：
本研究では、文部科学省研究委託事業「首都直下地震防災・減災特別プロジェクト」および「都市の脆弱性が引き起こ
す激甚災害の軽減化プロジェクト」による首都圏地震観測網（MeSO-net）の記録を用いた。
参考文献：
石井 透（2012）：日本建築学会構造系論文集，第 77 巻，第 676 号，pp.843〜850，2012.6
新井健介・ほか（2014）：日本地震学会講演予稿集，2014 年度秋季大会，S16-P19，2014.11
防災科学技術研究所： 地震ハザードステーション J-SHIS； http://www.j-shis.bosai.go.jp/

キーワード: 東京, 地震動, 継続時間, 表面波, 地下構造
Keywords: Tokyo, Earthquake Motion, Duration, Surface Wave, Structure

1/1

SSS25-P06

会場:コンベンションホール

時間:5 月 25 日 18:15-19:30

強震動の継続時間と変位振幅を用いた地震規模決定の可能性
Possible quantification of earthquake size using strong motion duration and displacement
amplitude
原 辰彦 1∗
HARA, Tatsuhiko1∗
1
1

建築研究所 国際地震工学センター
IISEE, BRI

Hara (2014) showed the correlation between strong motion duration and duration of high band-pass filtered tele-seismic P
waves. In this study, we investigated the possibility to quantify earthquake sizes using strong motion durations and displacement
amplitudes referring to Hara (2007), who showed that it is possible to determine magnitude for huge earthquakes using highfrequency energy radiation durations and displacement amplitudes from tele-seismic P waves. We used horizontal componentsof strong motion seismograms from the KiK-net borehole stations for nine large earthquakes that occurred in and around Japan
since 2003, which were retrieved and processed by Hara (2014). We obtained displacement records by integration of strong
motion seismograms and measured amplitudes after applying high-pass filter with the corner frequency of 0.01 Hz. We used thestrong motion durations obtained by Hara (2014), who used high band-pass filter with the corner frequencies of 5 and 10 Hz.
He suggested the duration dependence on epicentral distance. We calculated the logarithms of the products of the strong motion
durations and displacement amplitudes divided by the epicentral distance twice and compared them to the moment magnitudes.
Although the scatter is large, they correlate well, which suggests the possibility to quantify the earthquake size using durations
and displacement amplitudes measured from strong motion seismograms.
Acknowledgements. We used strong motion seismograms recorded by the KiK-net of the National Research Institute for Earth
Science and Disaster Prevention. We used the hypocenters of the unified hypocenter catalog of the Japan Meteorological Agency.

キーワード: 強震動, 継続時間
Keywords: strong motion duration

1/1

SSS25-P07

会場:コンベンションホール

時間:5 月 25 日 18:15-19:30

地盤の経験的伝達関数に基づく線形時と非線形時の計測震度の関係式
Relations between seismic intensity by linear and nonlinear site responses using empiricaltransfer functions
佐藤 智美 1∗SATOH, Toshimi1∗
1
1

清水建設
Shimizu Corporation

1923 年関東地震や南海トラフ沿いの巨大地震などの歴史地震の震度分布に合う広帯域震源モデルの推定や強震動予測
では、地盤の非線形性を考慮して震度を推定する必要がある。しかし、経験的地盤増幅率を用いた統計的グリーン関数法や、経験的グリーン関数法での計算波は、地盤が線形時のものである。また、地盤の非線形性を考慮した解析手法を用いる場合には、地盤モデルと動的変形特性等を設定する必要があるが、強震観測点でも S 波検層やボーリング調査がなされているとは限らず、動的変形特性試験はほとんどされない。さらに、動的変形特性の設定方法は様々であり、地盤応答解析手法の違いでも結果は異なる。そこで、本研究では、多数の KiK-net 記録に基づき、地盤が線形時と非線形時の経験的伝達関数を求め、線形時と非線形時の計測震度の関係式を作成した。
はじめに、以下のようにして、地盤の線形解析と等価線形解析で用いる理論伝達関数の代わりに、経験的伝達関数を用いて、地盤の非線形性の影響を受けている地表の観測波から、地盤が線形時の波形を推定する。用いたデータは、全国の KiK-net 観測点の強震動と、その観測点での弱震動である。強震動は、68 地震の PGA が 300 cm/s2 以上の記録、弱震動は 520 地震の PGA が 20〜100cm/s2 の記録を選択した。この 124 観測点における記録から、さらに、観測波形に、
ロッキング振動や浮き上がり振動の影響のあるデータを除いた。最も影響が大きかった記録は、東北地方太平洋沖地震
で震度７相当となった KiK-net 芳賀 (TCGH16) であった。水平成分の経験的伝達関数は、強震動の地中の上下成分に対する地表の水平成分のスペクトル比に対する、弱震動の地中の上下成分に対する地表の水平成分のスペクトル比の比で定義した。上下成分の経験的伝達関数は、地表の水平成分の代わりに、地表の上下成分を用いることにより、同様に定義した。そして、複素フーリエスペクトルである経験的伝達関数を強震記録の地表のスペクトルに乗じ、フーリエ逆変
換して、線形時の地表での波形を推定した。地中の上下成分には、地盤の非線形性の影響がほとんどないと考えられる
ため、地表の記録から非線形地盤応答で地中の上下成分のスペクトルを推定し、これに線形地盤応答を考慮して地表でのスペクトルを算出することに相当する。この手法では、地盤の非線形性を考慮した解析手法で必要となる地下構造モ
デルや動的変形特性の設定が必要ない。
そして、地盤が線形時の計算波形から算出した計測震度 ILIN と等価卓越周波数 (fe) をパラメータとした地盤が非線形時の計測震度 IN ON の経験式を、ILIN ≧ IN ON の条件付き最小二乗法により作成した。得られた式は、IN ON -ILIN =6.155-1.669 ILIN + 0.110 ILIN 2 -0.688log10 feである。ここで、fe は、計算波の最大加速度 PGALIN と最大速度 PGVLIN から fe= PGALIN /(2πPGVLIN ) で定義した。fe の代わりに Vs30 を用いた検討も行ったが、fe の方がデータの再現性が良かった。これは、fe には、サイト特性のみではなく、震源特性と伝播特性の影響も含まれているためと考えられる。図には、IN ON と ILIN の関係と、得られた経験式を示す。最大の観測計測震度 IN ON は、2000 年鳥取県西部地震の KiK-net 日野 (TTRH02) での 6.6 である。ILINの範囲は 4.5〜7.0 であり、用いた経験的伝達関数は 192 個である。得られた経験式は、震度 5 強以上で、INON と ILINに差が生じはじめる。ひずみが 10−3 程度以上に大きくなると動的変形特性の変化率が小さくなることを考慮するため、ILIN 2 項を導入して、差にはゆるやかな頭打ちを持たせている。fe が大きいほど ILIN と IN ON の差が大きくなる特徴があるが、これは、高周波数ほど、地盤の非線形性の影響により地盤増幅率が低下するためと考えられる。作成した経験
式では、ILIN が 7.0 の場合、fe =5Hz では IN ON は約 6.4、fe=2Hz では IN ON は約 6.7 となる。
謝辞：本研究は、科学研究費補助金基盤研究 (A) 26242034 及び (A) 23241054 による成果である。本研究では、防災科学技術研究所の KiK-net 強震記録、PS 検層結果、気象庁の震源情報を用いました。

1/2

SSS25-P07

会場:コンベンションホール

2/2

時間:5 月 25 日 18:15-19:30

SSS25-P08

会場:コンベンションホール

時間:5 月 25 日 18:15-19:30

Nonlinear Site Response Analysis in Tokyo Meteropolitan Area Using Equivalent Linear
ApproachNonlinear Site Response Analysis in Tokyo Meteropolitan Area Using Equivalent Linear
Approach
IBRAHIM, Rami1∗ ; MASUDA, Tetsu1 ; KOKETSU, Kazuki1 ; HIROSE, Takeshi2
IBRAHIM, Rami1∗ ; MASUDA, Tetsu1 ; KOKETSU, Kazuki1 ; HIROSE, Takeshi2
1
1Earthq. Res. Inst., Univ. Tokyo, 2 East Nippon Expressway Company Limited
Earthq. Res. Inst., Univ. Tokyo, 2 East Nippon Expressway Company Limited

Nonlinear site effect analysis was carried out for 210 simulated strong ground motion records in the Tokyo Metropolitan area.
The simulated strong motion records synthesized for two scenario earthquakes, i.e. plate-boundary and intra-slab ones, were theoutputs from the Special Project for Earthquake Disaster Mitigation in Tokyo Metropolitan Area. Fault plane of plate boundary
earthquake extends along the surface of subducting Philippine Sea Plate (PHS), located beneath the northern part of Tokyo bay
area. Its presumed magnitude of Mw 7.3. The fault plane of intra-slab earthquake is expected within the PHS at a focal depth of
50 km, extended beneath Chiba and Ibaraki prefectures. Its presumed magnitude of Mw 7.1. All waveforms were synthesized onengineering base of shear-wave velocity 500 m/s.
Site response analyses were achieved by an equivalent linear method using DYNEQ program developed by Yoshida and
Suetomi (1996). G-γ and h-γ relationships of Central Disaster Mitigation Council are recognized to express the shear deformationof soil (clay, sand, and gravel). Resultant waveforms on free surface shows a systematic dependence on thickness of the soft
structure above the engineering base. Large amplification is dominant at corresponding short periods above shallow soft structure,
whereas the peak amplitudes shifted to longer periods for sites located above deep soft structure. Nonlinear site effects weretypically obvious at short periods of approximately 0.2 s and shorter. The deformations mostly concentrated in the shallow 20 mof the soil inferred from maximum shear strains analysis. Our analysis showed the considerable effects of the non-linear responseof surface layers to large seismic inputs on the engineering base.

キーワード: Nonlinear site effect, Plate boundary earthquake, Intra-slab earthquake
Keywords: Nonlinear site effect, Plate boundary earthquake, Intra-slab earthquake

1/1

SSS25-P09

会場:コンベンションホール

時間:5 月 25 日 18:15-19:30

東濃地震科学研究所の地震観測点戸狩および TRIES における雑微動最小振幅と地震
動最大振幅との関係
Relation between smallest microtremor amplitudes and largest seismic oscillations observed by TRIES
田中 寅夫 1∗ ; 大久保 慎人 1
TANAKA, Torao1∗ ; OKUBO, Makoto1
1
1

無職
NOTHING

東濃地震科学研究所（TRIES) の地震観測点は 1999 年 5 月に 02TRIES 点、8 月には 00 戸狩点が設置されて、同年内に10 点の観測網が完成した。この新設時期に観測された雑微動と地震動の関係を調べるため、02TRIES 観測点を「基準点」とし、他の観測点における雑微動と地震動の観測値は、それぞれ 00TRIES の値で除した相対値で与えることとした。設
置された地震計は固有振動数 3.0Hz の過減衰型加速度計であり、1.0〜10.0Hz の帯域でほぼ水平な特性である。トリガーが作動して、100Hz サンプリング記録が始まり、遅延がかけられた初動までの雑微動の最小振幅、それに続く地震動から最大振幅を求めている。雑微動については数〜10 秒間、地震動振幅については P 波、S 波、表面波部分を数秒〜10 秒
間、それぞれ離散フーリエ変換した。東濃地域における地震動観測から明らかにされてきた site effect に基づいて、2.0〜
4.0Hz の振動数帯域が地震動災害に大きく関わると考え、上記帯域を幅 0.1Hz、21 区間に分割して、各区間での地震動振幅最大値と雑微動振幅最小値を計算した。02TRIES 観測点では 70 個を超える雑微動記録から最小振幅はかなり高い信頼性で決定できたが、今回報告する 00 戸狩観測点は 14 個の雑微動データしか解析できておらず、まだその信頼性は高くない。結果は全て雑微動振幅の最小値、地震動振幅の最大値ともに、02TRIES 点での値で除した相対値で与えて地震のマグニチュードや震央距離などの影響を少なくして site effect を強調することとした。2 観測点だけの予備的解析結果であるが、雑微動最小振幅値と地震動最大振幅値の組み合わせから siteeffect との有意な関連性を導くことが出来そうである。雑微動最小振幅から震源、マグニチュード、メカニズム等を想定すれば 02TRIES 観測点に対する最大地震動振幅の予想が出来ると考えている。
キーワード: 雑微動, 地震動, 離散フーリエ変換, 最小振幅, 最大振幅, サイト・エフェクト
Keywords: microtremor, seismic waves, discrete Fourier transform, ground soil, maximal amplitude, site effect

1/1

SSS25-P10

会場:コンベンションホール

時間:5 月 25 日 18:15-19:30

沖縄島中部から奄美群島南部におけるＳ波速度構造探査
Microtremor surveys for S-wave velocity structure in Okinawa and Amami islands, Japan山田 伸之 1∗ ; 竹中 博士 2
YAMADA, Nobuyuki1∗ ; TAKENAKA, Hiroshi2
1
1

福岡教育大学 教育学部, 2 岡山大学 自然科学研究科
Fukuoka University of Education, 2 Okayama University

琉球海溝に沿って点在する南西諸島では，島ごとに頻度の程度にばらつきはあるものの，過去に主にプレート境界で発生する地震による被害に見舞われている．また，一部の島々では，近い将来に震度 6 弱以上の強い揺れに見舞われる可能性が高いと評価されているとともに，活断層の存在も指摘される島もあり [例えば，地震調査研究推進本部]，将来的に地震被害を受ける危険性は低くないといえる．
これまでに，著者らは，奄美大島から八重山諸島の南西諸島の各島々での地震観測点周辺を中心に微動アレイ探査を実施し，暫定版として報告してきた（例えば，山田・竹中 (2013, 2014 各秋季大会)）．その中でも，Yamada & Takenaka(2014)
では，沖縄島南東部の南城市から那覇市，座間味島，久米島にかけての東西を縦断する地点での微動アレイ探査の結果を報告し，地点毎のＳ波速度構造を解明してきた．こうして，具体的な物性値であるＳ波速度構造情報を蓄積させることは，地震動評価の精度向上のためのモデル化に必要なことであると考えられる．今回は，これまでにまだ報告してきていなかった地点の深部地盤のＳ波速度構造を明らかにしたものの報告を行う．なお，この地域では，防災科学技術研究所の J-SHIS の深部地盤モデルが公表されているものの，Ｓ波速度の探査はこれまでになされてはいない．
本発表では，沖縄島中北部から奄美群島南部にかけての強震観測点 (K-NET) 周辺で実施した微動アレイ探査結果について報告する．対象地点は，沖縄島内のうるま (URM)，名護 (NGO)，国頭 (KGS) および奄美群島南部の与論島茶花 (YRN)，
沖永良部島知名 (OKE)，徳之島伊仙 (TKS) の各地点である．アレイ観測の実施は，URM と NGO は 2014 年 2 月，KGS,
YRN, OKE は 2009 年 3 月，TKS は 2010 年 2 月および 2011 年 9 月に行われたもので，観測には，強震計とロガーの組み合わせを 7 セット用いた．なお，探査サイトにより使用した観測機器が異なるが，データ解析への支障はなかった．観測については，これまでの報告と同様に，常時微動の収録を行った．与論島・沖永良部島・徳之島の各島でのスペクトル形
状は，周期 4 秒前後に非常に明瞭なピークを呈していたが，沖縄島内の地点ではあまり明瞭でなかった．こうして得られた観測データから，位相速度を求めると，いずれの地点も，周期 1 秒前後で 2.0 km/s を超える似た傾向を示していた．
今後，詳細な速度構造の解明を行うとともに，これまでの探査結果との対応関係や J-SHIS の速度構造モデルとの比較を行い，強震記録を参考にしながら，逆解析で得られた速度構造の妥当性についての検討も視野に入れたい．
なお，この研究の一部は，科学研究費基盤研究 (B) (課題番号：26282105) の補助を受けています．記して，感謝いた
します．
キーワード: 沖縄島・与論島・沖永良部島・徳之島, S 波速度構造, 微動アレイ探査
Keywords: Okinawa, Yoron, Okinoerabu and Tokunoshima island, S-wave velocity structure, Microtremor array measurement

1/1

SSS25-P11

会場:コンベンションホール

時間:5 月 25 日 18:15-19:30

地震波干渉法による別府湾および周辺地域の表面波群速度の推定ー広帯域稠密地震
観測記録に基づく検討
Estimation of surface wave group velocities beneath the Beppu bay, western Japan, usinga dense broadband array
林田 拓己 1∗ ; 吉見 雅行 2 ; 徳丸 哲義 3 ; 杉山 長志 4
HAYASHIDA, Takumi1∗ ; YOSHIMI, Masayuki2 ; TOKUMARU, Tetsuyoshi3 ; SUGIYAMA, Takeshi41
1

建築研究所 国際地震工学センター, 2 産業技術総合研究所 活断層・火山研究部門, 3 徳丸技術士事務所, 4 なし
IISEE, Building Research Institute, 2 Geological Survey of Japan, AIST, 3 Tokumaru Professioal Engineer’s Office, 4 Free

１．はじめに別府湾および周辺地域における深部地盤構造モデルの検証を目的として，別府湾周辺の連続微動記録を用いた地震波
干渉法解析を実施している．林田・吉見（2014，地震学会）では，大分県中部地域および愛媛県佐田岬に設置されている（独）防災科学技術研究所 Hi-net 連続観測記録を用いて表面波群速度を推定し，既存の地震波速度構造モデル（J-SHIS,
2014）による理論値との比較を行った．しかしながら，湾の空間スケールに対して観測点間距離が長大であるため，詳細な地盤構造の空間変化に関する検討を十分に行うことが出来ない．また，短周期地震計記録に基づく解析では 0.1Hz 前後の波群の同定が困難であることから，地盤構造の空間変化スケールを考慮した，より高密度かつ広帯域での連続観測を行う必要がある．
２．別府湾周辺地域における地震動の連続観測2014 年 8 月以降，別府湾を中心とした領域の 12 か所（南北 30km，東西 65km の範囲）において地震動の連続観測を実施している．使用機器は広帯域速度計（Nanometrics 社製 Trillium Compact：固有周期 120 秒）およびデータ収録装置（白山工業製 LS-8800：100Hz サンプリング）であり，地震計の埋設深度は 30cm である．2014 年 9 月-11 月の観測記録からパワースペクトル密度を確認したところ，1Hz 以下の低周波帯域では全地点の微動レベルに大差がなく，0.08Hz-0.1Hz，
0.15Hz-0.25Hz，0.5Hz 前後においてスペクトルのピークが見られた．
３．地震波干渉法による表面波群速度の推定
本研究では，取得した連続波形データ（３成分）を１時間毎のセグメントに分割し，振幅値の正規化処理（ Bensen
et al., 2007）を行った後，0.05Hz-2.0 Hz の帯域において全 66 観測点ペア（観測点間距離：6km-65km）の相互相関関数を導出した．本検討では上下動成分（Z-Z），水平動成分（R-R, T-T）および上下動?水平動成分（Z-R, Z-T）の 5 成分を対象とし，観測期間全体のスタッキング処理により得られたものを最終的な相互相関関数とした．多くの相互相関関数には明瞭な波群が見られ，堆積層が厚く分布すると考えられる領域を跨ぐペアにおいては，概ね 0.1Hz-1.0Hz の周波数帯域において波群の分散性が見られた．但し，他地域での検討事例［林田・他（2014, 地震 2）］と同様，相互相関関数のS/N 比が 2 を下回る場合には群速度の推定が困難となり，多くの観測点ペアでは未だ十分な S/N 比が得られていないことから，今後も引き続き観測記録を蓄積しながら，地震波干渉法の適用条件について検討を重ねる必要がある．
謝辞：本研究は，文部科学省研究委託事業「別府−万年山断層帯（大分平野−由布院断層帯東部）における重点的な調査観測」の一環として実施しています．

キーワード: 微動, 地震波干渉法, 表面波, 群速度, 深部地盤構造
Keywords: ambient noise, seismic interferometry, surface wave, group velocity, beppu sedimentary basin

1/1

SSS25-P12

会場:コンベンションホール

時間:5 月 25 日 18:15-19:30

別府湾周辺域（大分平野および別府扇状地）における微動アレイ探査
Microtremor array survay in Oita and Beppu area
吉見 雅行 1∗ ; 林田 拓己 2 ; 徳丸 哲義 3 ; 杉山 長志 4
YOSHIMI, Masayuki1∗ ; HAYASHIDA, Takumi2 ; TOKUMARU, Tetsuyoshi3 ; SUGIYAMA, Takeshi41
1

産総研 活断層・火山研究部門, 2 建築研 国際地震工学センター, 3 徳丸技術士事務所, 4 なし
Geological survey of Japan, AIST, 2 IISEE, Building Research Institute, 3 Tokumaru Professioal Engineer’s office, 4 Free

別府湾周辺域の S 波速度構造を明らかにするため，大分平野内 12 カ所，別府扇状地５カ所にて微動アレイ探査を行った．アレイは正三角形頂点と重心からなる 4 点アレイとし，1 カ所あたり半径 20m?1000m 程度の６アレイの観測を実施した．このうち半径約 300m，500m，1000m の３アレイは交通雑音を避けるため夜間の同時観測（12 点）とした．使用センサーは東京測振製速度計 SE-321（固有周期 10 秒，5V/kine），白山工業製 LS-8800 データロガーにて 24bit，100Hz
サンプリングで 3 成分の微動記録を収録した．
取得データの上下成分のうち明らかなノイズ等を除去して SPAC 法，V 法（Tada et al,2007），CCA 法（Cho et al. 2006）
解析を実施し，周波数ごとの位相速度を推定した．概ね 0.2Hz から 5Hz 程度までの周波数帯域にて位相速度が求められた．
観測位相速度を既存地下構造モデルの理論位相速度と比較した．既往の微動アレイ探査が実施されている大分平野では，観測とモデルからの理論値が概ね良い対応を示した．別府扇状地では既往地下構造モデルの理論位相速度は１ Hz 以
下の低周波数領域で観測値より大きな値を示し，速度構造に修正が必要なことがわかった．
観測位相速度曲線から１次元速度構造を推定した．逆解析には山中・石田（1995）の遺伝的アルゴリズム（GA）を用
い，観測位相速度の残差が小さくなるような S 波速度構造モデルを探索した．
本研究は文部科学省研究委託事業「別府ー万年山断層帯（大分平野ー由布院断層帯東部）における重点的な調査観測」の一環として実施した。微動アレイデータの解析には BIDO2.0（Tada et al, 2010，http://staff.aist.go.jp/ikuo-chou）を使用
しました．
キーワード: 速度構造, 別府堆積盆地, 微動
Keywords: velocity structure, beppu sedimentary basin, ambient noise

1/1

SSS25-P13

会場:コンベンションホール

時間:5 月 25 日 18:15-19:30

微動・重力探査に基づく鳥取県鹿野地域における地盤震動特性の評価と地盤構造の推定
Estimation of Subsurface Structure based on Microtremor and Gravity Survey in the
Shikano Area, Tottori Prefecture
野口 竜也 1∗ ; 香川 敬生 1 ; 津吉 祐典 1 ; 小倉 康平 1 ; 吉田 昌平 1
NOGUCHI, Tatsuya1∗ ; KAGAWA, Takao1 ; TSUYOSHI, Yusuke1 ; OGURA, Kohei1 ; YOSHDA, Shohei1
1
1

鳥取大学
Tottori University

1943 年鳥取地震の際に，地表地震断層として鹿野断層および吉岡断層が生じた。本研究では，鳥取県の鹿野地域において，微動および重力探査により地下構造を推定し，地震被害が生じた地域の地盤震動特性を調べた。微動探査では 156
地点の単点 3 成分観測および 12 地点のアレイ観測を行った。微動探査の解析では，アレイ観測記録より得られる位相速度より S 波速度構造，単点 3 成分記録から得られる H/V スペクトルのピーク周期より堆積層の層厚分布を推定した。重力探査ではラコスト・ロンバーグ重力計を用いて 38 点で重力観測を実施し，解析では重力異常データベースおよび鳥取東部地域の既存重力データに今回得られた観測記録を含めて，仮定密度 2.4t/m3 として重力異常を求め，2 次元および 3次元の密度構造を推定した。
キーワード: 微動探査, 重力探査, 地盤構造, 鹿野断層
Keywords: microtremor survey, gravity survey, subsurface structure, hikano area in Tottori Prefecture

1/1

SSS25-P14

会場:コンベンションホール

時間:5 月 25 日 18:15-19:30

1943 年鳥取地震の学校被害地点における常時微動による地盤構造推定と地震動評価
Estimation of Subsurface Structure from Microtremor and Strong Ground Motion Evaluation of the 1943 Tottori Earthquake
吉田 昌平 1∗ ; 貝崎 達憲 1 ; 野口 竜也 1 ; 香川 敬生 1
YOSHDA, Shohei1∗ ; KAIZAKI, Tatsunori1 ; NOGUCHI, Tatsuya1 ; KAGAWA, Takao1
1
1

鳥取大学
Tottori University

1943 年 9 月 10 日，鳥取市東部を震源に鳥取地震が発生し，鳥取市街地を中心に甚大な被害を出した．地震発生から70 年が経過したが，被害状況を記す資料はほとんど残っておらず当時の被害状況は十分に把握できていない．しかし，学校建物被害に関する資料や鳥取市街地の被害写真など，わずかに残った資料から被害が明確にされている場所も存在
する．本研究では，残存資料から把握された建物被害地点を対象に微動 3 成分単点観測と微動アレイ観測を実施し，詳
細な地下構造の推定と推定モデルを用いた地震動評価により被害の原因について検討した．3 成分単点観測記録を用いた
水平動と上下動のスペクトル比 (H/V) により対象地点の卓越周期を推定したところ，地盤の卓越周期と建物の固有周期が一致する地点は少なく，H/V だけでは十分に説明できなかった． 微動アレイ観測記録を用いたフォワード解析により
S 波速度 700m/s までの地下構造モデルを得ることができ，鳥取市街地では S 波速度 100˜200m/s 程度の層が 30m 程度堆積していることがわかった． 地震動評価では統計的グリーン関数法を用いて工学基盤波形を推定し，等価線形応答解析
により地表面地震動の推定を行った．周期 1.0 秒程度の応答値が大きい地点で被害が大きい傾向があることがわかった．
キーワード: 微動観測, 鳥取地震, S 波速度構造, H/V
Keywords: Microtremor observation, Tottori Earthquake, S-wave velocity structure, H/V

1/1

SSS25-P15

会場:コンベンションホール

時間:5 月 25 日 18:15-19:30

地震波干渉法を用いた立川断層帯周辺におけるレイリー波群速度の評価
Evaluation of Rayleigh-wave group velocities using seismic interferometry in the vicinityof Tachikawa fault zone
石毛 宏和 1∗ ; 地元 孝輔 1 ; 佐口 浩一郎 1 ; 山中 浩明 1 ; 酒井 慎一 2 ; 蔵下 英司 2 ; 平田 直 2
ISHIGE, Hirokazu1∗ ; CHIMOTO, Kosuke1 ; SAGUCHI, Koichiro1 ; YAMANAKA, Hiroaki1 ; SAKAI, Shin’ichi2 ;KURASHIMO, Eiji2 ; HIRATA, Naoshi2
1

東京工業大学大学院総合理工学研究科, 2 東京大学地震研究所
Interdisciplinary Graduate School of Science and Engineering, Tokyo Institute of Technology, 2 Earthquake Research Institute,
the University of Tokyo1

近年注目を集めている地盤探査手法の一つに地震波干渉法がある。本研究ではこの手法を用いて立川断層帯周辺の表
面波群速度の推定と評価を行った。解析はまず各観測点間において相互相関関数に基づくグリーン関数を推定し、表面波群速度の分散曲線を求めた。次に、地元・山中 (2013) によって提案された地震波干渉法の適用範囲を用いて抽出したデータをもとに、表面波のスローネストモグラフィ解析を行った。解析されたトモグラフィマップから、立川断層帯の南西側では長周期の表面波ほど群速度の大きい範囲を南側に広げ、北東側とは明らかに構造が異なることが分かった。最後に、解析結果と地震調査研究推進本部の「長周期地震動予測地図」による既往の地下構造モデルを用いて算出した群速度を比較した。その結果、断層の南西側で推定された群速度の大きい領域が、既往のモデルに比べ短周期帯で広く現れており、既往のモデルよりも表層の層厚が小さいことを示唆する結果となった。
キーワード: 立川断層帯, 地震波干渉法, レイリー波, 表面波群速度
Keywords: Tachikawa fault zone, Seismic interferometry, Rayleigh wave, surface-wave group velocities

1/1

SSS25-P16

会場:コンベンションホール

時間:5 月 25 日 18:15-19:30

立川断層近傍における短周期地震動の空間変動特性
Characteristics of Spatial Variation of Short-Period Ground Motion in the Vicinity ofTachikawa-Fault
津野 靖士 1∗ ; 佐口 浩一郎 2 ; 地元 孝輔 2 ; 佐藤 浩章 3 ; 松島 信一 4 ; 重藤 迪子 5 ; 高井 伸雄 5 ;神野 達夫 6 ; 山中 浩明 2 ; 川瀬 博 4
TSUNO, Seiji1∗ ; SAGUCHI, Koichiro2 ; CHIMOTO, Kosuke2 ; SATO, Hiroaki3 ; MATSUSHIMA, Shinichi4 ;MICHIKO, Shigefuji5 ; TAKAI, Nobuo5 ; KANNO, Tatsuo6 ; YAMANAKA, Hiroaki2 ; KAWASE, Hiroshi4
1

鉄道総合技術研究所, 2 東京工業大学, 3 電力中央研究所, 4 京都大学, 5 北海道大学, 6 九州大学
Railway Technical Research Institute, 2 Tokyo Insitute of Technology, 3 Central Research Institute of Electric Power Industry,
4
Kyoto University, 5 Hokkaido University, 6 Kyusyu University
1

表層地盤内には、不規則な形状の地層境界や S 波速度の横方向の変化など様々な不均質性が存在している。表層地盤の不均質性は、局所的な地震被害集中や大加速度生成などの現象（山中ら，2011）の主要な原因となることが多い。本研
究では、こうした特異な地震動の分布に関した現象を理解することを目的として、表層地盤内に段差状の構造がある立川断層近傍に於いて、断層直交方向にアレー地震観測を実施し、収録された観測地震動の空間変動特性について検討した。
東京都武蔵村山市に位置する真如苑グラウンド内で地震観測を実施し、加速度計 JEP-6A3 と収録器 Datamark LS7000XT（200Hz サンプリング）を使用した。観測期間は 2014 年 6 月 17 日〜20 日の計 4 日間である。立川断層に直交する形で、
約 650m の測線上に地震計を 8 点設置した。6 月 20 日に発生した千葉県北西部の地震（Mj 4.6）と茨城県沖の地震（Mj3.9）が収録され、それらデータをもとに短周期地震動の空間変動特性を評価した。測線のある観測点を境に観測波形が顕著に異なっており、波形全体のフーリエ・スペクトルはその観測点の東側と西側で、4Hz 以上の周波数帯で振幅が異なり、10Hz 付近でその差が最も顕著である。このことは、表層地盤の速度構造が異なるためであると考えられる。一方
で、測線の最も東側に位置する観測点を基準点としたコヒーレンスは、水平距離が離れる観測点ほど低下し、この場合
0.5-1Hz の周波数でその変化が最も顕著である。このことは、水平方向に伝播する表面波の影響が現れているためと考えられる。
本研究では、立川断層に直交する測線上の観測点に収録された地震データをスペクトル解析することより、表層地盤の不均質性による短周期地震動の空間変動特性を検討した。今後は、Goodness-of-fit の手法を用いて、より定量的に地震動の空間変動特性を評価する予定である。
キーワード: 短周期地震動, 空間変動, 表層地盤構造, 地震観測, 立川断層
Keywords: short-period ground motion, spatial variation, shallow soil structure, earthquake observation, Tachikawa-Fault

1/1

SSS25-P17

会場:コンベンションホール

時間:5 月 25 日 18:15-19:30

微動アレイ観測による大崎市古川地区の浅層地盤構造の推定
Estimation for velocity structure of shallow sediments using microtremor array observation in Furukawa, Japan
満永 仁志 1∗ ; 飯山 かほり 1 ; 盛川 仁 1 ; 後藤 浩之 2 ; 稲谷 昌之 2 ; 羽田 浩二 3 ; 池田 貴昭 2 ;高屋 俊康 2 ; 木村 紗也佳 2 ; 秋山 良平 2 ; 澤田 純男 2
MITSUNAGA, Hitoshi1∗ ; IIYAMA, Kahori1 ; MORIKAWA, Hitoshi1 ; GOTO, Hiroyuki2 ; INATANI, Masayuki2 ;HADA, Koji3 ; IKEDA, Takaaki2 ; TAKAYA, Toshiyasu2 ; KIMURA, Sayaka2 ; AKIYAMA, Ryohei2 ;SAWADA, Sumio2
1
1

東京工業大学, 2 京都大学, 3 株式会社ニュージェック
Tokyo Institute of Technology, 2 Kyoto University, 3 Newjec INC.

On March 11, 2011, the Tohoku earthquake (Mw 9.0) brought the large ground motion and tsunami in the eastern part ofmainland Japan. As a result, many buildings were seriously damaged and many people were killed. Furukawa of Miyagi prefecture is one of the most seriously damaged areas by the earthquake ground motion. The anomaly of the damage distribution
was observed in this area. It is considered that the anomaly caused by difference in ground motion characteristics, and especiallyshallow ground structure affects ground motion characteristics. To identify the characteristics, some observations have been carried out in Furukawa, for example, very dense seismic array observations (Goto ea al., 2012) and so on. The observation covers
2?3km2 area by 23 sensors and provides the relative information on depth to ground basement. However, it is still difficult todetermine the absolute values of the depth to the basement.
In order to clarify the shallow ground structure in the area, the array observation of microtremors was carried out from December 13 to 16, 2014. The three components seismometer was adopted, that is, moving-coil type velocity sensor with the natural
frequency of 2 Hz. We set three arrays in the north-western area of Furukawa station, which were the equilateral triangle arrayswith radii of 2.9 m to 45 m.
The microtremor data were analyzed by the spatial auto-correlation (SPAC) method (aki, 1957). In the surface layer of ground,the shear-wave velocity is approximately 100 m/s at all the sites and the depth is 10 or 15 m. Also, the depth to the groundbasement is around 35 m at all the sites. On the other hand, the differences of the ground properties are found in middle layers.
It, however, is difficult to recognize the differences on the horizontal-to-vertical spectral ratios (H/V), because of the similarityof the predominant frequencies.
The estimated S-wave velocity structures suggest that the properties of middle layers affect to the anomaly of damage distribution. This means that the impedance ratio between middle layer and basement plays the important roles to the amplification of
the ground motions in this area.
Keywords: microtremor, velocity structure, array observation, Furukawa, ground motion, spatial auto-correlation method

1/1

SSS25-P18

会場:コンベンションホール

時間:5 月 25 日 18:15-19:30

大崎市古川地区の浅層地盤構造をターゲットとした超高密度常時微動観測
Super dense microtremor observation evaluating shallow ground structure in Furukawa
area
後藤 浩之 1∗ ; 稲谷 昌之 2 ; 羽田 浩二 3 ; 池田 貴昭 4 ; 高屋 俊康 4 ; 木村 紗也佳 2 ; 秋山 良平 2 ;澤田 純男 1 ; 飯山 かほり 5 ; 盛川 仁 5
GOTO, Hiroyuki1∗ ; INATANI, Masayuki2 ; HADA, Koji3 ; IKEDA, Takaaki4 ; TAKAYA, Toshiyasu4 ;KIMURA, Sayaka2 ; AKIYAMA, Ryohei2 ; SAWADA, Sumio1 ; IIYAMA, Kahori5 ; MORIKAWA, Hitoshi51
1

京大防災研, 2 京大院工, 3 ニュージェック, 4 京大工, 5 東工大院総合理工
DPRI, Kyoto Univ., 2 Eng. Dept., Kyoto Univ., 3 Newjec Inc., 4 Eng., Kyoto Univ., 5 Science Eng., Tokyo Inst. Tech.

2011 年東北地方太平洋沖地震では，東北地方の内陸部において地震動による被害も局所的に発生している．宮城県大崎市古川地区は，地震動による被害が顕著であった地区の１つである．古川地区内でも被害の発生したエリアは限定的であったため，地盤震動特性の違いなどがその原因として予想された．そこで，2011 年 9 月より，地震観測点を高密度に配置して地盤震動特性の違いを明らかにする研究が進められている（Goto et al., 2012）．地震記録に基づいて浅層地盤構造を推定した結果，被害の発生したエリアでは相対的に表層地盤が厚く堆積していることが示唆されている（稲谷他，
2013）．
地震観測点は，市街地中心部（1km x 1.5km）に 25 台が展開されており，およそ 350m 四方に 1 台程度の密度である．
地震観測網としては世界最高クラスの密度ではあるが，実際に地盤震動特性を詳細に議論しようとするには十分でなかった．そこで，地震観測点間を埋めた高精細な浅層地盤構造を把握することを目的として，より高密な単点常時微動観測を実施した．
単点常時微動観測は 2014 年 9 月 30 日から 10 月 2 日，および 12 月 14 日から 17 日の計７日間実施した．対象とする
古川地区市街地の 527 点で計測した．なお，測定点の間隔は約 50m となるように設定している．各測定点で 11 分間の常時微動を計測し，計測記録からそれぞれ H/V スペクトル比を求める．H/V スペクトル比のピーク値を与える周波数を各測定点で抽出し，地図上にプロットした．
H/V スペクトル比のピーク周波数は，古川地区市街地で概ね 1-2Hz の範囲に求められる．測定エリア内でピーク周波数が相対的に高い（1.5-2.0Hz）エリアは，地震記録により表層厚さが浅めに推定されたエリアと概ね整合する．また，地震動被害の生じたエリアはピーク周波数が相対的に低く（1.0-1.5Hz）求められている．ただし，ピーク周波数は空間的に相当ばらつくため，常時微動により浅層地盤を分析するためには，そのばらつきを如何に取り扱うかが重要である．
参考文献
Goto et al.,: Very dense seismic array observations in Furukawa district, Japan, Seism. Res. Lett., 83(5), pp.765-774, 2012.
稲谷他：大崎市古川高密度地震観測を利用した表層地盤構造の推定，土木学会論文集 A1，69(4)，pp.I 758-766，2013．

キーワード: 常時微動観測, H/V スペクトル比, 地盤震動, 古川地区
Keywords: micrtoremor observation, H/V spectral ratio, site factor, Furukawa area

1/1

SSS25-P19

会場:コンベンションホール

時間:5 月 25 日 18:15-19:30

断層近傍の詳細地盤モデル作成の検討−深谷断層帯・綾瀬川断層帯を例として−
Construction of underground structure model in the fault neighborhood(An example of
Fukaya fault and Ayasegawa fault).
先名 重樹 1∗ ; 松山 尚典 2 ; 神 薫 2 ; 藤原 広行 1
SENNA, Shigeki1∗ ; MATSUYAMA, Hisanori2 ; JIN, Kaoru2 ; FUJIWARA, Hiroyuki1
1
1

防災科学技術研究所, 2 応用地質
NIED, 2 OYO Corp

１．はじめに防災科研では、南関東全域においてボーリングデータ等の地質・地盤資料に基づく初期地質モデルの構築と、そのモデルを基に地震記録と微動観測記録により、物性値（主に S 波速度）を調整し、周期・増幅特性を考慮した浅部・深部
統合地盤構造モデルの試作を行ってきている。これらは広域地盤モデル構築手法の標準化の取組として、地震本部にて「地下構造モデル作成のレシピ」を構築している。さらに、地盤モデル高度化の検討として、詳細な地盤モデル（詳細地盤モデル）の構築時に必要となる、活断層近傍における地震動評価やモデル作成時の地盤の不整形性の検討および強震動時の非線形特性を評価できる地盤構造モデルの構築等の構築を検討している。本検討では、南関東地域に存在する深谷断層帯・綾瀬川断層帯を例として詳細地盤モデル構築の検討結果について報告する。
２．断層帯周辺の既往調査と地質構造の概要
深谷断層帯・綾瀬川断層は、地震本部の長期評価見直しで関東平野北西縁断層帯から変更されている。この断層自体の活動度は低いが、連動すると M8 クラスの地震が発生すると予測される。深谷断層は、明瞭な重力異常分布境界となっている。なお、深谷断層については杉山・他 (2009) により、反射法地震波探査、ボーリング、トレンチ等の調査を行われており、断層の構造や活動性の検討がなされている。反射法探査や既往の反射断面の解釈により、深谷断層から北東側では、基盤岩上面が深度３ｋｍ付近まで落ちていること、その上に中新世以降の地層が厚く堆積していることが確認されている。深谷断層は、南側の平井断層、櫛引断層と合わせて地質構造の形成過程が検討されている。
３．調査概要上記の断層近傍の地盤構造モデルを構築するために、地震観測および微動観測を実施した。観測は断層を挟むように5 測線を展開し、単点微動約 200m 間隔、極小アレイ 1km 間隔、大アレイが 2km 間隔で実施している。大アレイの位置には地震観測も実施した。５測線のうち１測線は、杉山他 (2009) における反射法地震波探査断面の測線近傍で実施している。なお、断層の北東側は沖積低地、南東側がローム台地で構成されている。
（１） 地震観測
断層近傍の１４か所において臨時地震観測を実施した。観測した記録について、フーリエスペクトルを地点ごとにまとめたところ、地震によって異なる特性が確認できるが、地点ごとにスペクトルに一定の傾向が確認出来た。そこで、各地点とのスペクトル比を取ったところ、各地点ともスペクトル比は地震に寄らず安定しており、サイト増幅特性として利用できる。
（２）常時微動観測
微動観測の結果、単点の H/V スペクトルのピーク周期は、断層の両側で異なる傾向を示す。断層の落ち側（北東側）では、上がり側（南西側）に比べて、H/V スペクトルのピーク周期が周期の長い方へシフトする（周波数が小さくなる）傾
向が確認出来た。微動探査 (アレイ) 地点に近い、地震動の R/V スペクトルと比較したところ、低周波数域のスペクトル
特性も調和的であることを確認している。微動アレイ解析で得られた S 波速度構造と J-SHIS の深部地盤モデルの S 波速度構造を比較すると、地震基盤相当の Vs=3.2km 層の上面は、今回の結果の方が全体的に約 500m 程度深くなっている。
一方で、Vs=1.5Km 層の上面は、断層落ち側で 300〜500 ｍ程度浅くなっている。また、Vs=0.9km 層の上面は、断層落ち側で最大で 900m 程度深くなることが分かった。工学的基盤以浅の地盤構造では、断層落ち側の断層直近で Vs=0.3km な
いし 0.2km 以下の層が厚い（十数ｍ）。同様の傾向は、現在作成中の防災科研の南関東地盤モデルでも確認されている。
断層周辺のボーリングデータでみると、本数が少ないので、詳細は不明であるが、断層落ち側の断層直近で N 値の小さい完新統（粘性土層など）が厚くなっている傾向があるようにみえ、前述の S 波速度構造分布と調和的である。
４．まとめと今後の検討
本研究の結果、断層の落ち側では、上り側に比べて微動の H/V スペクトルの周波数が小さくなる傾向が明瞭で双方の構造の変化が比較的良く確認でき、観測された地震動の R/V スペクトルも同様の特性を示すことがわかる。また、深谷
断層を挟む両側の地盤では、スペクトル特性、S 波速度構造にも相違が明瞭にみられ、ボーリングデータからも判別できることから、この断層においては地震観測記録・微動記録による詳細地盤モデルの構築・検証は十分に可能であること

1/2

SSS25-P19

会場:コンベンションホール

時間:5 月 25 日 18:15-19:30が分かった。
今後、地震観測データを用いた、スペクトルインバージョン等による、断層両側の地盤での地震動の増幅特性の定量的な検討および、ボーリングデータから読み取れる地下浅部の地質構成、地質構造や重力探査結果をふまえた浅部・深部地盤の速度構造モデルの作成を試み、観測された地震動の作成した地盤モデルを用いた検証や工学的基盤の不整形性の影響評価を実施する予定である。

キーワード: 常時微動, 活断層, 強震動予測, 地盤モデル
Keywords: microtrmor observation, active fault, strong motion prediction, underground structure model

2/2

SSS25-P20

会場:コンベンションホール

時間:5 月 25 日 18:15-19:30

単点常時微動観測を用いた工学的基盤の傾斜箇所の判定に関する基礎的検討
A study for detection of the site on inclined engineering bedrock based on microtremor
observation
坂井 公俊 1∗ ; 田中 浩平 1 ; 宇佐美 敦浩 1 ; 小林 薫 2 ; 平林 雅也 2
SAKAI, Kimitoshi1∗ ; TANAKA, Kohei1 ; USAMI, Atsuhiro1 ; KOBAYASHI, Kaoru2 ; HIRABAYASHI, Masaya2
1
1

鉄道総合技術研究所, 2 東日本旅客鉄道
Railway Technical Research Institute, 2 East Japan Railway Company

工学的基盤が傾斜しているような地点では，地震動が局所的に増幅する可能性があり，土木構造物の耐震設計の観点からも注意が必要な箇所に該当する．ここで，2 点の常時微動同時観測記録の鉛直スペクトル比がグリーン関数の比になるという理論的な報告がされており，この方法を用いることで，工学的基盤の傾斜箇所を簡易に抽出できる可能性がある．そこで，工学的基盤が傾斜している地点において常時微動観測を実施することで，地盤傾斜抽出の可能性について検討を行った．その結果，基盤が概ね水平な地点間の鉛直スペクトル比は概ね 1.0 であるのに対し，基盤が傾斜している
地点間では振幅比の大きさ，山と谷の振動数も大きく変化しており，この方法によって基盤傾斜が抽出できる可能性があることを確認した．今後は本地点を対象とした数値解析を実施することで結果の妥当性を確認するとともに，基盤傾
斜箇所の簡易抽出法を構築する予定である．
キーワード: 常時微動, 表面波, 傾斜地盤
Keywords: microtremor, surface wave, irregular ground

1/1