地震波のサイト増幅特性評価の高度化とそれを用いた地下構造モデルの検証に関する研究

過年度までに開発した，スペクトルインバージョン法に用いる観測方程式の震源特性に破壊伝播による影響を含む項を加えた定式化を行い，2016年熊本地震の地震系列を対象として適用した．今年度は，先行研究を参照してこの破壊伝播項を明示的にバイラテラル破壊伝播モデルを仮定した定式化を進め，破壊伝播方向，破壊伝播速度，破壊伝播比（ユニラテラルからバイラテラルの間）のパラメータを推定する定式化を行った．この方法により，これまでと同じく2016年熊本地震系列のデータに適用した．震源を取り囲む観測点データが準備できた場合には，これらのパラメータは信頼性高く求められることがわかった一方，十分にまわりを取り囲んでいない場合は，パラメータ推定精度が低くなることがわかった．また，この方法によって推定された震源項を元に応力降下量を推定した．応力降下量の推定値分布は，コンベンショナルな方法に比べてバラツキが小さくなったものの，イベント毎の特徴（規模，深さ，発生位置）は残っており，イベント毎の応力降下量のバラツキは大きいことがわかった．加えて，スペクトルインバージョンで得られた各観測点の地震波サイト増幅特性を用いて，そのサイト増幅が観測点直下の地盤構造によってもたらされるものと仮定して，サイト増幅特性による地盤構造モデルの推定を行った．解析周波数帯（0.2-10Hz）では，深い地盤構造のみならず，浅い地盤構造による増幅の影響が大きく，このサイト増幅特性を用いた地盤構造（S波速度構造）のモデル化や，既往のモデルがある場合の検証に使うことができる．

過年度までに開発した，スペクトルインバージョン法に用いる観測方程式の震源特性に破壊伝播による影響を含む項を加えた定式化を行い，2016年熊本地震の地震系列を対象として適用した．今年度は，先行研究を参照してこの破壊伝播項を明示的にバイラテラル破壊伝播モデルを仮定した定式化を進め，破壊伝播方向，破壊伝播速度，破壊伝播比（ユニラテラルからバイラテラルの間）のパラメータを推定する定式化を行った．この方法により，これまでと同じく2016年熊本地震系列のデータに適用した．震源を取り囲む観測点データが準備できた場合には，これらのパラメータは信頼性高く求められることがわかった一方，十分にまわりを取り囲んでいない場合は，パラメータ推定精度が低くなることがわかった．また，この方法によって推定された震源項を元に応力降下量を推定した．応力降下量の推定値分布は，コンベンショナルな方法に比べてバラツキが小さくなったものの，イベント毎の特徴（規模，深さ，発生位置）は残っており，イベント毎の応力降下量のバラツキは大きいことがわかった．加えて，スペクトルインバージョンで得られた各観測点の地震波サイト増幅特性を用いて，そのサイト増幅が観測点直下の地盤構造によってもたらされるものと仮定して，サイト増幅特性による地盤構造モデルの推定を行った．解析周波数帯（0.2-10Hz）では，深い地盤構造のみならず，浅い地盤構造による増幅の影響が大きく，このサイト増幅特性を用いた地盤構造（S波速度構造）のモデル化や，既往のモデルがある場合の検証に使うことができる．