(6)逆解析

水平成層構造モデルから理論位相速度を数値計算するのに最低限必要なモデルパラメータは,層数n,各層毎の層厚h,弾性波速度(P波速度α及びS波速度β)及び密度

ρである(i=1,2,…,n)。実際の逆解析作業では,数値計算の効率化及び安定化を図るため,未知のモデルパラメータを層厚h及びS波速度βの2種類とし,P波速度α及び密度ρは松岡・白石・梅沢(2000)が本調査地域の北部に位置する埼玉県南部地域における微動アレイ解析の際に提案した統計回帰式

β=0.64α−0.51,

ρ=0.22α+1.51,

単位はαi[km/sec],β[km/sec],ρ[g/cm3],

を用いてS波速度βから換算した。

遺伝的アルゴリズムに基づく逆解析計算(筒井・藤本,1994;長ほか,1999)において,未知のモデルパラメータに係る広域探索範囲は,前出の表2−2−7に示すとおりに設定した。局所探索範囲幅は広域探索範囲幅の10%とし,遺伝子ビット長は広域探索,局所探索にいずれについても7ビットとした。理論位相速度の計算においては,従来から採用されている「基本モードパワーが最も優勢」の仮定に則り,基本モード分散のみを考慮した。試行5000回×実験20回の計算により求めた計20個の候補解を平均して,最終解とした。

各調査地点毎に求めた候補解及び最終解を,図2−2−35図2−2−36図2−2−37に示す。各図の上2図は候補解,同じく下2図は最終解である。各図の右上図及び右下図はS波速度深度分布(ただし,第5層までを表示),同じく左上図及び左下図は観測位相速度と理論位相速度との比較である。

図2−2−38には,各調査地点のS波速度深度分布のみ(ただし,第5層までを表示)を一括表示してある。また,最終解に対するモデルパラメータ(層厚,各層上面深度,P波速度,S波速度及び密度)を,表2−2−9に一括して示す。

なお,図2−2−35図2−2−36図2−2−37の左下図には参考として,最終解を利用して計算した高次モード(1次モード)の理論位相速度(赤色破線)と,各モード毎のmedium response(基本モード=水色実線,1次モード=水色破線)とを併せて表示した。高次モードの位相速度,あるいはmedium responseについては,考察(後節2.2.4)で詳しく述べる。