９－５

近地強震記録によって推定された 2016 年 10 月 21 日に鳥取県中部で発生した地震の震源過程
Rupture process of the 2016 Central Tottori prefecture earthquake derived fromnear-source strong-motion records
防災科学技術研究所
National Research Institute for Earth Science and Disaster Resilience

2016 年 10 月 21 日 14 時 7 分に鳥取県中部で発生した地震（M 6.6; 気象庁）について，防災科学技術研究所の K-NET，KiK-net の強震波形記録を用いた震源インバージョン解析を行った．
◎断層面モデルと震源過程のパラメータ化
F-net のモーメントテンソル逆解析 1)で得られたメカニズム解や震源位置の情報から同地震は地殻内で発生した浅い横ずれ断層型地震と考えられる．余震活動が北北西―南南東方向に広がっている
ことから，断層面として F-net モーメントテンソル解析結果で得られた二つの節面のうち余震活動に調和的な走向を持つ解（走向 162 度，傾斜 88 度）を採用した（第 1 図）．断層面の大きさは，予備的な解析の結果から走向方向の長さを 20 km，傾斜方向の幅を 16 km とした（第 1 図）．破壊開始点は，Hi-net 震源に基づき，北緯 35.3806 度，東経 133.8545 度，深さ 11.58 km においた．
本解析ではマルチタイムウィンドウ線型波形インバージョン法 2,3)に基づき，断層破壊過程を時空間的に離散化した．空間的には，断層面を 2 km 四方の小断層で，走向方向 10 個，傾斜方向 8 個に分割した．時間的には，各小断層でのすべり時間関数を，破壊開始点から一定速度 Vftw で広がる同心円が到達した時刻から，0.8 秒幅のスムーズドランプ関数を 0.4 秒ずらして 5 個並べることにより表現した．これにより，各小断層からの要素波形（グリーン関数）を通じて，断層破壊過程と各観測点での波形は線型の方程式で結び付けられる．
◎波形データとグリーン関数
防災科学技術研究所の K-NET 地表観測点 7 点，KiK-net 地中観測点 5 点，KiK-net 地中観測点 3点の計 15 観測点での強震加速度波形記録を積分して得た速度波形を用いた（第 2 図）．これらの速度波形に 0.1－1.0 Hz のバンドパスフィルタをかけ，5 Hz にリサンプリングし，S 波到達 1 秒前から10 秒間を切り出し解析データとした．
各小断層からの要素波形（グリーン関数）は，一次元地下構造モデルを仮定し，離散化波数積分法 4)と反射・透過係数行列法 5)により点震源の波形を計算して求めた．地下構造モデルは，全国深部地盤構造モデル 6)の各観測点直下の情報を用いて観測点ごとに構築した．KiK-net 観測点については速度検層の情報も利用した．また小断層内部の破壊伝播の効果を 25 個の点震源（走向方向，傾斜方向それぞれ 5 列）を小断層内に分布させることにより表現した．
◎波形インバージョン
各小断層の各タイムウィンドウでのすべり量を，観測波形と合成波形の差を最小とするように，最小二乗法を用いて求めた．不等式拘束条件をつけた最小二乗法 7)を用いて，各小断層でのすべり方向の変化を，F-net メカニズム解のすべり角である–11 度の±45 度以内に収めた．また時空間的に－ 360 －近接するすべりを平滑化する拘束条件 8)を付加した．平滑化の強さは ABIC9)を参考とし，Vftw は残差を最小とするものを選んだ．
◎結果
第 2 図に推定された断層面上の最終すべり分布を示す．第 3 図に断層破壊の時間進展過程を示す．
第 4 図に観測波形と理論波形の比較を示す．Vftw は 3.1 km/s，最大すべり量は 0.7 m，断層面全体での地震モーメントは 2.2×1018 Nm（Mw 6.2）である．大きなすべりが破壊開始点から浅い側にかけた領域（最大すべり量 0.7m）と破壊開始点から見て北北西の領域（最大すべり量 0.6m）において見られる． 各領域での破壊はそれぞれ地震発生から 0–3 秒後と 3–5 秒後に生じている．第 5 図に前震および余震の空間分布と最終すべり分布の比較を示す．本震発生後一時間の間に余震活動が活発であった領域と本震時に大きくすべった領域は重ならない．また最大前震（2016 年 10 月 21 日 12時 12 分，M 4.2）およびそれに続いた地震活動は本震の震源付近で発生していた．
注：なお，本解析は暫定的なものであり，今後修正される可能性がある．
（文責：久保久彦・鈴木亘・青井真（防災科研）・関口春子（京大防災研・防災科研））

引用文献
1) 福山英一・石田瑞穂・D. S. Dreger・川井啓廉 (1998), オンライン広帯域地震データを用いた完全自動メカニズム決定, 地震 2, 51, 149-156.
2) Olson, A. H. and R. J. Apsel (1982), Finite faults and inverse theory with applications to the 1979 Imperial
Valley earthquake, Bull. Seismol. Soc. Am., 72, 1969-2001.
3) Hartzell, S. H. and T. H. Heaton (1983), Inversion of strong ground motion and teleseismic waveform datafor the fault rupture history of the 1979 Imperial Valley, California, earthquake, Bull. Seism. Soc. Am., 73,1553-1583.
4) Bouchon, M. (1981), A simple method to calculate Green's function for elastic layered media, Bull. Seism.
Soc. Am., 71, 959-971.
5) Kennett, B. L. and N. J. Kerry (1979), Seismic waves in a stratified half space, Geophys. J. R. astr. Soc.,
57, 557-583.
6) 藤原広行・河合伸一・青井 真・森川信之・ 先名重樹・工藤暢章・大井昌弘・はお憲生・早川 讓・遠山信彦・松山尚典・岩本鋼司・鈴木晴彦・劉瑛 (2009), 強震動評価のための全国深部地盤構造モ
デル作成手法の検討, 防災科学技術研究所研究資料, No.337.
7) Lawson, C. L., and R. J. Hanson (1974). Solving Least Squares Problems, Prentice-Hall, Inc., New Jersey,
340 pp.
8) Sekiguchi, H., K. Irikura, and T. Iwata (2000), Fault geometry at the rupture termination of the 1995
Hyogo-ken Nanbu earthquake, Bull. Seism. Soc. Am., 90, 117–133.
9) Akaike, H. (1980), Likelihood and the Bayes procedure, in Bayesian statics, edited by J. M. Bernardo, M.
H. DeGroot, D. V. Lindlely, and A. F. M. Smith, University Press, Valencia, Spain.

－ 361 －第 1 図 解析に用いた観測点の分布および断層面の地表投影．星印は破壊開始点を示す．
Fig.1 Station distribution on the map. The assumed fault plane is also projected on the map. A star denotes
the rupture starting point.

第2図
Fig.2

断層面上の最終すべり分布．ベクトルは上盤のすべり方向とすべり量を示す．星印は破壊開始点を示す．
Total slip distribution on the fault. The vectors denote the direction and amount of the slip of thehanging wall side. A star denotes the rupture starting point.

－ 362 －

第 3 図 破壊の時間進展過程．1 秒ごとのすべり分布を示す．
Fig.3 Rupture progression in terms of slip amount for each 1 s time window.

第 4 図 観測波形（黒）と理論波形（赤）の比較．右上の数字はそれぞれの最大値（cm/s）を示す．
Fig.4 Comparison between the observed and the synthetic waveforms. The maximum values are shown on
the upper right of each waveform.

－ 363 －

第5図
Fig.5

前震および余震の空間分布と最終すべり分布の比較．震源情報には Hi-net 震源カタログを用いた．下図は上図の緑枠内の領域に対応する断面図である．
Comparison of spatial distribution of foreshocks and aftershocks (M ≥ 2) determined by Hi-net withthe final slip distribution of the mainshock. The lower figure is the cross section for the region of thegreen rectangle in the upper figure.

－ 364 －