水平成層構造モデルから理論位相速度を数値計算するのに最低限必要なモデルパラメータは，層数ｎ，各層毎の層厚ｈ_ｉ，弾性波速度（Ｐ波速度α_ｉ及びＳ波速度β_ｉ）及び密度

ρ_ｉである（ｉ＝１，２，…，ｎ）。実際の逆解析作業では，数値計算の効率化及び安定化を図るため，未知のモデルパラメータを層厚ｈ_ｉ及びＳ波速度β_ｉの２種類とし，Ｐ波速度α_ｉ及び密度ρ_ｉは松岡・白石・梅沢（２０００）が本調査地域の北部に位置する埼玉県南部地域における微動アレイ解析の際に提案した統計回帰式

β_ｉ＝０．６４α_ｉ−０．５１，

ρ_ｉ＝０．２α_ｉ＋１．５１，

単位はα_ｉ［ｋｍ／ｓｅｃ］，β_ｉ［ｋｍ／ｓｅｃ］，ρ_ｉ［ｇ／ｃｍ^３］，

を用いてＳ波速度β_ｉから換算した。

理論位相速度の計算においては，計算作業の効率向上を図るため，従来の解析例と同じく「基本モードのパワーが最も優勢」という仮定を採用し，基本モード分散のみを考慮した。そのため，最初にフォワードモデリングによるテスト計算を行い，その結果基本モードのみでは再現困難と判断した観測位相速度を，逆解析計算の際に除外した。

遺伝的アルゴリズムに基づく逆解析計算（筒井・藤本，１９９４；長ほか，１９９９）において，未知のモデルパラメータに係る広域探索範囲は，前出の表２−２−７−１、表２−２−７−２に示すとおりに設定した。局所探索範囲幅は広域探索範囲幅の１０％とし，遺伝子ビット長は広域探索，局所探索にいずれについても７ビットとした。試行５０００回×実験２０回の計算により求めた計２０個の候補解を平均し，必要に応じてさらにパラメータの微調整を加え，最終解とした。

調査地点毎に求めた候補解及び最終解を，図２−２−３５、図２−２−３６、図２−２−３７、図２−２−３８、図２−２−３９、図２−２−４０に示す。各図の上段には候補解，同じく下段には最終解を示す。各図の右上図及び右下図はＳ波速度深度分布（ただし，第７層までを表示），同じく左上図及び左下図は観測位相速度と理論位相速度との比較である。観測位相速度の表示において，細破線付き○印で水色表示したデータは，逆解析計算の際に除外したデータを意味する。

図２−２−４１，及び図２−２−４２には，各調査地点のＳ波速度深度分布のみ（ただし，第７層までを表示）を一括表示してある。また，最終解に対するモデルパラメータ（層厚，各層上面深度，Ｐ波速度，Ｓ波速度及び密度）を，表２−２−９に一括して示す。