研究概要 |
本研究の目的は高精度かつ高解像度の3次元地球内部構造(地震波速度構造)を推定することである.このような研究は理論地震波形及び偏微分係数の計算が膨大になり,これらを効率的かつ高精度に計算することが必要不可欠である.このために,高精度理論波形計算のための新しい離散化手法を開発した(Geller & Takeuchi 1995,GJI).本年度は実際のデータ解析により内部構造推定を行なう予定であったが,より高精度のインバージョン実現のための手法開発を先に行なうことにした. 本年度は第一に新しい離散手法を任意3次元不均質媒質の場合に適用できるように,定式化を拡張した.また,この定式化を用いて2次元不均質媒質に対する理論波形計算を行なうソフトウェアを開発した(Cummins et al.1997,GJI).本ソフトウェアでは,従来よく用いられてきた,不均質構造を球対称構造からの無限小の摂動であるとする仮定を用いていないため,強度の不均質構造に対しても正確な計算が可能である.理論波形計算により計算された波面と,波形計算と独立に波線理論を用いて計算された波面と良く一致することが確認され,精度が十分あることが確認された. 第二に,新しい離散化手法を時間領域の計算の場合に拡張した(Geller & Takeuchi 1997,IASPEI meeting;Geller,Takeuchi & Mizutani 1997,OHP Symposium).これまで,本手法の応用は周波数領域の計算に限られていたが,時間領域の場合にも同様に拡張できることが確認された.
|