レシーバ関数の計算に使用した地震や検証に用いた地震観測点は、表５−５−１の通りである。ただし、データを調べた結果、観測点２２２０５Ａについては明瞭なピークが見られなかったため、また、観測点Ｇ４１２１ＡおよびＫＯＧＡに関しては１つの地震しか観測記録が得られていないため検証には使用していない。

レシーバ関数は、Ｐ波初動から３秒間のデータを抜き出し、ラディアル成分と上下成分の波形のフーリエスペクトルを計算し、ラディアル成分と上下成分のクロススペクトルを上下成分のフーリエスペクトルで除し、さらに時間領域に変換することにより求めた。時間領域に変換したあとに０．５Ｈｚのローカットフィルターおよび５Ｈｚのハイカットフィルターを施している。データ解析長は３秒間を基本としているが、Ｓ波主要動が含まれてしまう場合には適宜データ解析長を調整した（表５−５−１）。

一方、観測点それぞれについて、３次元地下構造モデルから抜き出したＰ波、Ｓ波の速度構造を用いて理論レシーバ関数を計算した。計算は水平成層構造を仮定し、離散化波数法（Ｂｏｕｃｈｏｎ Ｍ．，１９７９）により求めた。観測点毎に表５−５−２に示したＰ波の入射角の平均値を求め、その入射角で中心周期０．５秒のリッカーウエーブレットを入射させ、地表応答として現れる水平成分と上下成分の波形を観測波形に対して同様の処理を施し理論レシーバ関数を求めた。

図５−５−１には観測レシーバ関数と理論レシーバ関数の比較図を、図５−５−２には観測レシーバ関数より得られるＰＳ−Ｐ時間と理論レシーバ関数より得られるＰＳ−Ｐ時間の差の分布を示す。各観測地点で観測によるＰＳ−Ｐ時間と地盤モデルによるＰＳ−Ｐ時間は概ねよく対応しているが、濃尾平野南部の観測点１２１５１Ａ、１２１５２Ａ、１２１５３Ａでは、観測によるＰＳ−Ｐ時間に比べて地盤モデルによるＰＳ−Ｐ時間が０．２秒程度短くなっている。図５−５−３には観測ＰＳ−Ｐ時間と理論ＰＳ−Ｐ時間の比較図を示す。同図に示した理論ＰＳ−Ｐ時間は、観測点毎に表５−５−２に示したＰ波の入射角の平均値を求め、その入射角の場合の理論ＰＳ−Ｐ時間を示している。図５−５−２より、１２１５１Ａ、１２１５２Ａ、１２１５３Ａの３地点でのずれが大きいものの、観測地点間の相対的なＰＳ−Ｐ時間は観測と理論で概ね対応しており、修正した地下構造モデルの基盤深度の概略的な傾向は妥当なものと考えられる。

１２１５１Ａ、１２１５２Ａ、１２１５３Ａの３地点は主に基盤の深度とＳ波速度構造に着目したＨ／Ｖスペクトルによるピーク周波数は観測値とモデルによる理論値の対応がよいことから、ＰＳ−Ｐ時間でモデルによる理論ＰＳ−Ｐ時間が観測ＰＳ−Ｐ時間より短いのは、Ｓ波速度に対する相対的なＰ波速度が現実より遅いことが考えられる。図５−５−４に各検証地点で観測ＰＳ−Ｐ時間を満足するようにＰ波速度を調整したときのモデルのＰ波速度とＳ波速度の関係を、既往のＰＳ検層データのＰ波速度とＳ波速度の関係と重ねて示す。観測ＰＳ−Ｐ時間を満足させることでモデルのＰ波速度とＳ波速度の関係はばらつきが見られるが、既往のＰＳ検層データのばらつきの範囲内にあることがわかる。すなわち観測と理論でＰＳ−Ｐ時間に若干のずれが生じているのは、Ｐ波速度とＳ波速度の関係において本来は地点や地域による違いがあるが、濃尾平野全域に対して共通のＰ波速度とＳ波速度の関係式を用いていることが要因の１つと考えられる。

本来は、Ｐ波速度とＳ波速度の関係に地域による違いがあることを念頭に置かなければならないが、現状ではＰ波速度とＳ波速度を一部の地域ごとに単独で設定できるほどデータは十分ではなく、それをカバーするために速度境界面を適切に設定することが重要である。

表５−５−１　解析に用いた地震の最大振幅一覧

表５−５−２　Ｐ波入射角の一覧（鉛直下向を０度とする）

図５−５−１　観測地震波によるレシーバ関数と作成した地盤モデルより計算したＰＳ−Ｐ時間との比較

図５−５−２　観測ＰＳ−Ｐ時間と理論ＰＳ−Ｐ時間の差の分布

図５−５−３　観測ＰＳ−Ｐ時間と理論ＰＳ−Ｐ時間の比較

図５−５−４　Ｐ波速度とＳ波速度の関係