付録B．詳細法計算手法

　ここでは、統計的グリーン関数法の計算手法について示す。

　B.1　地震基盤における統計的グリーン関数の作成

　地震基盤における統計的グリーン関数は、佐藤（1994a, 1994b）が仙台地域で観測された主に海溝型地震の記録から推定したパラメータを用いたスペクトルモデルと経時特性モデルを用いて、Boore（1983）と同様の手法で作成した。
　地震基盤における S 波主要動のスペクトルは、次式で表される地震動の加速度フーリエスペクトルのモデルによって表現した。

ここに、 は要素断層に関する添字で、 は地震動の加速度フーリエスペクトル、 は地震波の放射特性、 および は要素断層における地殻の密度および剪断波速度、 は地震モーメント、 は臨界振動数、 は高周波遮断振動数、 は定数、 は震源距離、 は地殻の 値、 および は地震基盤の密度および剪断波速度である。 および は震源における密度3.03 、S波速度3.93 である。最終項は、自由表面の影響および要素断層における地殻のインピーダンスと地震基盤のインピーダンスとの相違 (佐藤、 1978) を考慮した項である（壇ほか、2000）。
　地殻の 値である は、

(B-1)

で表される値とした。（B-1）式は、佐藤ほか(1994b)の推定値であるが、最近の研究では、0.8Hz前後以下で 値が頭うちするあるいは減少することが指摘されていることから、0.8Hz以下で一定値とする。
　なお、 は、以下の式から算出できる。

　さらに、強震動を評価する領域内で、ある１つの地点を考えると、その地点は各要素断層からみて常に同じ方位に位置しているわけではないこと、および短周期の地震記録には明確な放射特性が見られないことから、 は平均的な値を用いることとした。計算地点と断層面との幾何学的関係、および断層のすべりのタイプを考慮して、 Boore and Boatwright (1984) に示された の値のうち、SH 波 (S 波のうち重力の方向に直交する成分) とSV 波 (S 波のうち SH 波に直交する成分) の値の大きいほうを採用する。
　一方、時刻歴波形の作成には経時特性あるいは位相特性が必要であるが、地震基盤におけるS波主要動の経時特性に関しては、現在までに研究成果が得られていない。しかし、本検討では、最終的に統計的グリーン関数を定義する位置を「詳細法工学的基盤」上とするため、佐藤ほか（1994a)の仙台地域の工学的基盤における地震記録から求められた経験的経時特性モデルを準用した。

ここに、

　気象庁マグニチュード は、佐藤(1989) による次式に基づき、地震モーメント から次式に基づき算出した。

　B.2　工学的基盤上面での統計的グリーン関数の作成

　各計算ポイント直下の三次元地盤モデルから、各計算ポイントでの1次元地盤モデルを作成し、S波の一次元重複反射理論により、B.1 において作成された地震基盤における統計的グリーン関数を入射波とした工学的基盤上面での統計的グリーン関数を計算する。

　B.3　工学的基盤上面における統計的グリーン関数を用いた波形合成

　B.2 で作成された工学的基盤上面における統計的グリーン関数を用いて、壇・佐藤（1998）の断層の非一様すべり破壊を考慮した半経験的波形合成法に従い、波形合成を行う。この際、大地震の要素断層のすべり量と応力降下量が再現されるように、震源スペクトルの補正を行った。このようにして算定された大地震の要素断層の波形を断層面全体の破壊過程を考慮して合成を行う。