強震アレイを用いたアスペリティのイメージングに関する研究

2003年9月26日04時50分ころに発生した十勝沖地震は、 断層面積が100km四方に及ぶM8クラスの巨大地震でした。この地震は、文部科学省によって全国におよそ20km間隔で強震計（K-NET：防災科学技術研究所）が設置・運用を開始して以来はじめての巨大地震であり、多くの研究者がその実像を確かめるべくさまざまなデータと解析手法を駆使してアスペリティ分布などを分析しました。

一般に断層面上のすべり量を推定する場合、断層面や地下構造などを仮定して理論波形を計算し、どのようなすべりが起これば実際の観測波形を説明できるか、理論波形と観測波形の一致の具合で判断します。そのような解析の結果（例えば、Honda et al. 2004, Koketsu 2004, Yagi 2004など）、襟裳岬の東方沖およそ80km付近の断層面上で、もっともすべり量が大きかったことが分かりました（図1のB）。
この点については複数の研究者の結果が一致しています。しかし、多くの観測点、特に北海道東部（道東）の観測点では、この最大のアスペリティ以外から来たと考えられる波群が観測されており、その発生源についてははっきりとしたことが分かっていませんでした。Honda et al. (2004)の解析（図１）では、Cと書かれたアスペリティがその発生源としてやや不明瞭ながら見えています。

道東は厚い堆積層が発達し、観測点ごとに地下構造が大きくことなります。また振幅の大きな表面波も発生します。このような場所では理論波形を正確に計算することが難しくなります。そのため道東で観測された波形をなるべく使わないようにして解析を行ったため、道東ではっきりと観測される波群の震源（アスペリティ）が結果にうまく現れなかった可能性があります。
　そこで本研究では、従来の解析手法とは違う手法を採用し、道東で観測される波群を発生させたアスペリティの位置を正確に推定しました。

本研究では、物理探査や散乱体の研究でよく使われているセンブランス解析の手法をアスペリティ位置の推定に応用します。

解析手法を図２に示します。断層面を仮定し2kmの小断層に区切ります。各小断層から各観測点までの走時（到達するまでにかかる時間）を計算し、　その時間差に応じて各観測点での観測波形の時間をずらしてセンブランス値Sを計算します。センブランス値は、おおまかには波形がどのくらい似ているか、という指標で、値が大きいほど波形が似ているといえます。ある小断層から地震波が放射されたと仮定してセンブランス値を計算し、そのセンブランス値を小断層上に投影すると、どこからその地震波が来たと考えるのが妥当か、ということが分かります。

このような解析をするためには、解析したい地震波の波長に対して十分狭い範囲になるべくたくさんの観測点（観測点の集まりをアレイといいます）が必要となります。そこで、今回は道東に設置されている防災科研のK-NETの観測点を仮想的にアレイとみなして解析を行いました（図３）。

さらに、センブランス値を計算する際には、なるべく精度よく断層面から各観測点までの走時が わかっている必要があります。そこで、本震の断層面周辺で発生した余震の波形を使い、その波形が再現できるような速度構造モデルを作成し、走時計算に使います。さらに、観測点ごとのローカルな速度構造の差を考慮するため、理論走時と観測走時を比較し、その差を補正してセンブランス値を計算しました。

本研究では、従来の波形逆解析では位置を推定するのが難しかった、道東で観測された波群の震源（アスペリティ）を精度よく求めることができました。図4に解析の結果得られたアスペリティの位置を示します。道東で観測されたパルス的な波形の震源（アスペリティ）は、釧路沖約65kmの深さ52km付近に存在することが分かりました。このアスペリティの位置は、M4以上の余震（Ito et al., 2004）や地震後のすべりの大きかった領域（Mityazaki et al., 2004）と相補的な関係にあり、この領域に蓄積されていた応力は本震時にほとんどが開放されたことを示唆しています。

また、そのアスペリティが破壊した時刻は、本震の震源で破壊がはじまってから約29秒後で、破壊開始点からこのアスペリティまでの平均的な破壊伝播速度は2.5km/sとなります。これは、最大アスペリティ（アスペリティB）の平均破壊伝播速度（3.6km/s）よりも有意に遅く、断層面上の破壊伝播速度が方位によって異なっていたと考えられます。

本研究では、従来の波形逆解析とは異なる新しい手法によって、大地震のアスペリティ分布を推定できることを示しました。アスペリティの位置やすべり量分布の精度の良い推定は、応力の蓄積・開放過程や地震の再来間隔の推定につながります。また、地震が発生したときの地震動を推定するための、震源モデルの高精度化にも貢献できます。 従来の解析でアスペリティ分布（すべり量分布）を推定するためには、断層の破壊伝播速度や震源メカニズムなど多くの制約条件を課して計算を安定させる必要がありますが、本研究で用いた手法では原理的にはそのような制約をかけなくても比較的安定して解が得られます。
また、理論波形の計算など時間のかかる重たい処理をする必要がないので、迅速にすべり領域の広がりなどを把握したい場合にも対応できると考えられます。
神奈川県西部地域では、以前からM7クラスの地震が発生する可能性が指摘されています。今回の解析で使用したK-NETの観測網ではM7クラスの地震を解析するには粗すぎますが、足柄平野などに展開されているSK-netなどのデータを利用すれば、そのような地震についても詳細な震源像を推定できると考えられます。