作成されたブーゲー異常図では負の異常が卓越しているが、これは本州の下に沈みこんでいる太平洋プレートとフィリピン海プレートおよびモホ面の形状等の影響によるものである（萩原、１９８６）。基盤構造解析のように表層付近の密度構造を解析する際にはこれらの広域的な成分は不要であるので、甲府盆地を中心に１００ｋｍ×１００ｋｍの範囲から求めた一次傾向面（図２−３−１０）を差し引いて広域成分を除去した図２−３−１１残差ブーゲー異常図を作成した。なお、図２−３−１２に示すフィリピン海プレ−トの重力値への影響を見ると、１次傾向面とは走行で約３５度の違いが認められるが、１次傾向面にはフィリピン海プレ−トの影響だけでなく、太平洋プレートやモホ面の形状の影響も含まれるものと判断し、１次傾向面による補正を実施した。断面解析用に設定した測線上のデータを残差ブーゲー異常値のグリッドデータから内挿で求め、二次元解析を行った。二次元解析に用いたプログラムは米国地質調査所が開発した「ＳＡＫＩ」である（Ｗｅｂｒｉｎｇ、１９８５）。

本プログラムでの重力値のフォワード計算はタルワニの方法が、またインバージョンアルゴリズムは修正マルカート法が用いられている。

解析では、堆積層に相当するポリゴンモデルを初期値として作成した。基盤に対する堆積層の密度差を仮定したポリゴンモデルによる重力値を計算し、計算された重力値と測定値との差を小さくするようにインバージョンアルゴリズムによってモデルの節点の座標の修正量を求め、修正したモデルから再度重力値を計算する。計算値と測定値の差が十分小さくなるまでこの過程を繰り返し、最終的な構造モデルを求める。

解析に際して、基盤に対する堆積層の密度差は國友・志知（１９８５）で使用された −０．３５〜−０．３７ｇ／ｃｍ^３ を採用した。また、地表で確認されている地質境界、断層等を考慮してモデルを作成した。Ａ−Ａ’断面では、反射法地震探査断面を参考に砂礫層を−０．６７ ｇ／ｃｍ^３と仮定して設定した。

図２−３−８にデータ処理及び解析のフローを示す。

図２−３−８　データ処理及び解析のフロー