�@位相速度の推定

図２−４８に示したフロ−に沿って位相速度の推定を行った。解析諸元は次のとおりである。なお、解析の過程では全データの解析区間を多数の小時間単位のブロックに分割し、各ブロック毎に所定の推定量（空間自己相関関数および係数など）を求め、それらを平均するいわゆる「ブロック平均法」を採用した。ブロック平均法は、時系列解析で各種推定値を求める際、それぞれの分散（ｖａｒｉａｎｃｅ）をできるだけ小さくするために広く用いられている。

・デ−タ編集 ：１０Ｈｚ（１００ｍｓ）のリサンプリング（高周波ノイズを取り除くために実施）

・観測デ−タ長：小、中アレー：６０分間

　　　　　　　　　　大アレー 　： １２０分間（６０分×２回）および

　　　　　　　　　　　　　　　　　　　２４０分間（１２０分間×４回、テスト１地点で実施）

　　　　　　　　　　　　　　　　　　　アレーサイズによって観測データ長を変えた

・ブロックデ−タ：小、中アレー：２０４．８秒

　　　　　　　　　　　大アレー 　 ：２０４．８秒および４０９．６秒の２通り

　※昨年度の解析結果をふまえ、長周期成分のデータを取得した

　※アレーサイズ（大アレー）は、ブロックデータを従来の２０４．８秒と４０９．６秒の２通りで主に解析を実施した。また、テスト１地点で１ブロック長８１９．２秒で解析を実施

・最大ブロック数：ブロックデータ２０４．８秒：３４個（４８ｻﾝﾌﾟﾙｵｰﾊﾞｰﾗｯﾌﾟ）

　　　　　　　　　　　ブロックデータ４０９．６秒：１９個（９６ｻﾝﾌﾟﾙｵｰﾊﾞｰﾗｯﾌﾟ）

・スペクトル分解バンド幅：０．０５Ｈｚ〜０．１Ｈｚ（パルゼンウィンドウ）

・座標入力： 各観測点の位置座標を１／２，５００より読みとり距離計算

・パワ−スペクトルおよびクロススペクトル計算：原デ−タの先頭、中間および最後部のパワ−スペクトル計算による時空間定常性の確認とノイズ区間除去後のデ−タの計算

・空間自己相関関数および空間自己相関係数の計算

・位相速度の推定：

アレ−サイズ別に空間自己相関係数の曲線形から、位相速度推定に採用する周波数範囲を決定した。すなわち、

（ａ）周波数の上限は、アレ−が見分け得る最短波長を最小のアレー半径の２倍と仮定し、それに対応する周波数（すなわち、空間自己相関係数＝Ｊ０（π） ≒−０．３０となる周波数）か、あるいはそれが不明瞭なデータの場合、空間自己相関係数が初めて極小（理論的には、空間自己相関係数≒−０．４０）となる周波数とした。

（ｂ）周波数の下限は、空間自己相関係数が推定されている最も低い周波数としたが、空間自己相関係数が周波数の減少とともに単調に減少する場合は、その減少し始める周波数とした。

平行層で近似することが適当でない複雑な地下構造の場合、同じアレ−サイズでもアレーの展開位置により空間自己相関係数に有意な差の生じることがある。今回の解析の場合は、アレーサイズ別にこの曲線から直接位相速度を求めた。また、同じアレーのうち、各地震計間の距離毎の空間自己相関係数のフィッティングを行い、距離別の位相速度を求める方法も採用した。ただし、空間自己相関係数に有意な差のなかった場合はそれらを平均化し位相速度を求めた。

各アレ−毎に求められた位相速度を統合し、最終位相速度を決定した。

�AＳ波速度構造を推定する逆解析

求められた位相速度の分散はレイリ−波の基本モ−ドであり、またアレ−直下の構造が平行層であると仮定して１次元のＳ波速度構造を求めた。

逆解析には、個体群探索分岐型遺伝的アルゴリズム（長ほか，１９９９、以下「ｆＧＡ」と称する）を用いた。図２−４９に逆解析のフロ−を示す。

Ｐ波速度および密度値についてはＬｕｄｗｉｇ，Ｎａｆｅ ａｎｄ Ｄｒａｋｅ（１９７０）の関係式を利用してＳ波速度値から換算した。