・Ｐ波及びＳ波初動読み取り（図３−３−３�@）

選定した１４個の地震全てが観測された青葉区榎が丘小学校の観測点の波形記録を一例として初動読み取り結果と併せて図３−３−４に示す。

・初動走時読み取りデータのチェック（図３−３−３�A）

今回の解析では、次の３点につきチェックを行い、必要に応じて、初動走時を再度読み取り修正する。

ａ）初動走時の観測点別分布傾向：Ｐ波（ＵＤ成分）の走時とＳ波（ＮＳ及びＥＷ成分）走時の分布傾向の類似性とＮＳ成分とＥＷ成分の読み取り走時に大きな相違がないか（図３−３−５）。

ｂ）走時曲線のばらつき：各成分の走時曲線が、地震毎に一定のまとまった傾向となっているかどうか（図３−３−６）。

ｃ）Ｓ／Ｎ比（ここでは、読み取った初動走時の前後２秒間のパワーの比としたので、値が大きいほど良好な記録であることに相当する）のチェック：各成分について観測点別にＳ／Ｎ比が悪いものがないかチェック（図３−３−７）及び地震別にＳ／Ｎ比のヒストグラム（図３−３−８）を作成しデータセットとして解析で使用可能か検討した。

以上の結果から、１９９８／０６／１４（千葉県東方沖）の地震は、走時のばらつきも大きく、Ｓ／Ｎ比も他と比較して悪いため解析には使用していない。また、図３−３−４からわかるように１９９８／０３／２３（茨城県沖）の地震に対する記録はＳ／Ｎが悪く、初動を読み取ることが不可能であったたので、解析には使用していない。

・初動走時データセットの完成（図３−３−３�B）

Ｓ波の初動走時としては、水平２成分のデータを読みとり、最終的には、Ｓ／Ｎ比が良い方の成分のデータを採用。

・走時遅れ分布を計算（図３−３−３�C）

図３−３−９−１、図３−３−９−２、図３−３−９−３、図３−３−９−４における走時分布を見ると、観測走時には明らかに震源位置の影響によると推定される走時の差が見られる。震源位置による影響が含まれていると、地下構造による走時の変化がとらえにくい。そこで、観測走時に幾何学的に平面をフィッティングすることにより、震央位置の違いによる走時の差（トレンド成分）を求め、それを原記録から差し引いた。その結果が、図３−３−９−１、図３−３−９−２、図３−３−９−３、図３−３−９−４に示すＰ波、Ｓ波走時遅れの分布である。

同図からわかるように、地震の震央位置の違いによらず、横浜市の南東部と北西部に走時遅れが大きい領域が存在し、中央やや西よりの地域で走時遅れが小さいことがわかる。

・初期モデル作成（図３−３−３�D）

１０ｋｍ以深は科学技術庁防災技術センターが震源決定で使用している速度モデルで固定し、浅部は、既往文献調査結果を基に１次元モデルを作成（図３−３−１０）。

図３−３−１１は、各地震毎に、１０ｋｍ以浅についても科学技術庁のモデルとした場合の理論走時と、浅部を今回作成したモデルとした場合の理論走時を併せて示しているが、全般的に観測走時と合っているのは後者である。

なお、１９９８／０１／１６（千葉県南部）のＰ波のように、系統的にずれているものについては、震源情報に誤差が含まれているものと推定される。

・屈折を考慮した数値シミュレーション（図３−３−３�E）

３次元レイトレーシング（ホイヘンスの原理に基づく爆発法）により、震源までを含めた領域を１ｋｍの立方体のグリッドに分割して理論走時を計算する。（図３−３−１２参照）

・観測データとの対比・チェック及び検討と速度モデルの修正

（図３−３−３�F）

横浜市域内について初期モデルを修正し、再度理論走時を計算し観測データと対比・チェックを行う。観測データと修正後の速度モデルにおける走時遅れ分布の対比図を図３−３−１３に示す。

なお、全ての地震に関する波形記録・初動走時データは横浜市総務局災害対策室に保管している。（Ａ４キングファイル５冊分）