5.4. すべり速度時間関数の設定

　　選定されたケース3 の断層パラメータを用いて、三次元有限差分法に用いるすべり速度時間関数を中村・宮武（2000）の方法により算定した。中村・宮武（2000）のすべり速度時間関数の近似式を計算するためには以下のパラメータが必要となる。それぞれのパラメータを求めるために与えるべき定数は、震源スペクトルから見積もられる（source-controlled ）、すべり時に働く動的な応力降下量、断層幅、ローパスフィルターのコーナー周波数（sourcecontrolled と同等）、破壊速度、剛性率、最終すべり量の計6 つである。

• （最大すべり速度到達時間）
  　震源スペクトルから見積もられる（source-controlled ）との関係（中村・宮武（2000）の(9)式）から求める。
  　　　（ に反比例）
• （最大すべり速度）
  　すべり時に働く動的な応力降下量、断層幅、ローパスフィルターのコーナー周波数（source-controlled fmaxaと同等）、破壊速度、剛性率との関係（中村・宮武（2000）の(10)式の誤植を下式のように修正）から求める。
  　　　（ に比例、、の平方根に比例、に反比例）
  　　　　　　　　　　　　　　　　 　　　　　(5.4−1)
• （ライズタイム）
  　断層幅と破壊伝播速度を用いて近似式から求める。
  　　　（ に比例、に反比例）
• （すべり速度振幅が に比例するKosttov 型関数に移行する時間）
  　すべり速度時間関数の近似式（中村・宮武（2000）の(7)式）で最終すべり量を与えることにより自動的に求めることができる。

　以上のような考え方で設定した、A1 断層、A2 断層に対するすべり速度時間関数のパラメータを、それぞれ、表5.4−1、表5.4−2 に示す。また、A1 断層、A2 断層に対するすべり速度時間関数を、それぞれ、図5.4−1、図5.4−2 に示す。実際の計算では、すべり速度時間関数を三角パルスの重ね合わせで近似して与えるため、図5.4−1、図5.4−2 には中村・宮武（2000）に基づくすべり速度時間関数の他、計算用のすべり速度時間関数を表示した。