アレイ解析から推定した新潟県中越沖地震の震源過程

地震の震源過程の解析には従来、波形バージョン法がよく用いられます。この手法では、コンピュータを用いて計算した理論波形と、実際の観測波形の一致の具合から、地震の際に断層面上のどの部分が何メートル滑った、というように物理的な値を推定できます。しかし、今回の地震のように観測点が厚い堆積層の上や地下の構造が複雑な場所にあるときは、理論波形を正確に計算することが困難であり、正確な震源象の推定することが難しい場合があります。

 震央距離約16kmにある柏崎刈羽原子力発電所敷地内では、97台の地震計のうち33台で記録された本震の近地強震動の波形が公開されています。 本報告では、震源断層に近い位置で、しかも空間的に高密度な観測点配置で得られた観測記録の利点を生かすために、アレイ解析という手法を用いて断層の破壊過程を推定することにしました。
この手法では、すべり量などの絶対値を推定することはできませんが、反射波や表面波など震源由来以外の波が含まれる場合でも、比較的安定して震源を特定することができ、また破壊伝播速度も自然に求めることができます。

震源過程の解析には、観測された水平動2成分の加速度波形に1Hzから20Hzのバンドパスフィルターをかけたものを使用しました。
断層モデルはAoi et al. （2008）のモデルと同じ30km×24kmの南東傾斜の断層面を仮定し、2km四方の小断層に分割しています。小断層と観測点のペアごとに理論走時を計算し、ある小断層について成分ごとに全ての加速度波形をスラントスタック（例えばYilmaz、1987）することによって、その小断層からのコヒーレントな波だけが強調され、残りの位相はキャンセルされた重合波形ができます。これをその小断層からのエネルギーの放射履歴とみなしました。

スタッキングを行う際、振幅が大きいだけでなくよりコヒーレントな波だけを強調するために、サンプリング時間の前後0.5秒間の時間窓でセンブランス値を計算し、重みとしてスタック波形に掛け合わせるsemblance enhanced stacking法（松本ほか、1999）を用いました。

本報告で用いた手法を応用すると、地震が発生し断層破壊が進展していく様子をほぼリアルタイムに近い形でイメージできる可能性があります。
その結果、例えば緊急地震速報で流される震度の予測（揺れはじめる時間や震度の大きさ）をより精度よく行うことができるようになるかもしれません。 神奈川県内や周辺は、東京大学のSK-netや防災科学技術研究所のK-NETなど、高密度な観測網があります。これらの観測網のデータを利用して、上にあげたような地震防災に直結するような研究が行えるかもしれません。

 各小断層から放出されたエネルギーの相対量を比較すると、破壊の開始点付近と断層の南西部分にピークが存在し、これらが観測波形に見られる2つの大きなパルスの震源、すなわちアスペリティと推定されます。アスペリティの分布は、従来の波形インバージョンによって得られた結果（例えば、Aoi et al., 2008）と調和的です。断層破壊は、震源から南西方向に向かってunilateral（一方向）に進み、柏崎市沖で停止したことが分かりました。また破壊伝播速度は約2450m/sと推定されました。

東京電力株式会社によって記録・配布された柏崎刈羽原発敷地内の強震動波形を使用させていただきました。記して感謝いたします