鉛直方向任意不均質弾性媒質における平面波入射問題の時間領域差分解法

  • 田中 宏樹
    九州大学理学部地球惑星科学教室 現在: 気象庁仙台管区気象台
  • 竹中 博士
    九州大学理学部地球惑星科学教室

書誌事項

タイトル別名
  • An Efficient FDTD Solution for Seismic Plane-wave Responses of Vertically Heterogeneous Media
  • エンチョク ホウコウ ニンイ フキンシツ ダンセイ バイシツ ニ オケル ヘイメンハ ニュウシャ モンダイ ノ ジカン リョウイキサ ブンカイホウ
  • An efficient FDTD solution for plane-wave response of vertically heterogeneous media

この論文をさがす

抄録

Plane-wave responses of vertically heterogeneous structure models (1-D media) are often computed in seismology. For horizontally layered media, they can be calculated by semi-analytical methods such as the propagator matrix method. However, for gradient velocity or randomly heterogeneous structures, we have to use numerical methods such as the finite-difference method. Conventional codes for the 2-D or 3-D finite-difference method require huge computer memory and long computation time even for calculating plane-wave responses of 1-D media. In this study we propose an efficient procedure to calculate plane-wave responses of arbitrary 1-D media using the finite-difference method in the time domain (FDTD). We first derive an elastodynamic equation of plane-wave incidence problem for vertically heterogeneous media by applying the Snell's law to 3-D elastodynamic equation. We then discretize the velocity-stress formulation of the derived elastodynamic equation using a staggered-grid finite-difference scheme of fourth-order accurate in space and second-order accurate in time. We also investigate the implementation of the stress-free surface condition for the scheme, and perform a stability check of the total scheme through actual computations. We computed plane-wave responses of a structure model with a velocity gradient using the derived finite-difference method. We focused on the PS-converted phase and found a “offset” phase appearing between the PS-converted phases generated at the top and bottom boundaries of the velocity-gradient layer on the surface responses of the structure model with the velocity gradient due to a P-wave incidence. This phase can be emphasized by calculating the receiver function from the radial and vertical waveforms. In this study we also investigate the “offset” phase attributed to the velocity gradient by numerical computations using the derived finite-difference method.

収録刊行物

  • 地震 第2輯

    地震 第2輯 57 (3), 343-354, 2005

    公益社団法人 日本地震学会

被引用文献 (4)*注記

もっと見る

参考文献 (23)*注記

もっと見る

詳細情報 詳細情報について

問題の指摘

ページトップへ