〔反射法データ処理〕
（１）フォーマット変換及びデータ編集
データは全てＧＤＡＰＳ−３探鉱機を用いＳＥＧＹフォーマットで取得されているが、データ処理用の内部フォーマットへの変換を行なった。さらに不良トレース等を除去した。
（２）最小位相変換
スウィープ波形の自己相関であるクラウダーウェーブレットを最小位相変換するオペレータを求めデータに適用した。この処理により、以後のホワイトニングデコンボリューションによって震源波形はパルスに変換されることになる（図２−３−６）。
（３）インパクタデータに対する震源補正
呉羽山トンネル部分でのインパクタを震源とするデータとそれ以外のバイブレータデータを統合して処理するために、インパクタの震源波形をバイブレータの震源波形に変換するオペレータを求め、これをインパクタによる全データに適用した。オペレータは、ＶＰ．４８７の同一発震点で、震源以外のパラメータを同一として取得したインパクタとバイブレータ記録から求めた。図２−３−７−１には震源補正前後のショット記録を、図２−３−７−２には、重合記録におけるインパクタ記録とバイブレータ記録を示した。インパクタとバイブレータ記録間に顕著な差異は認められず、震源補正が適切であったことを示している。
（４）ＣＭＰ編集
図２−３−８−１に示すＣＭＰ分布図をもとに、受振測線に出来るだけ近づけるように重合測線を設定した（図２−３−８−２，図２−３−８−３）。この重合測線に沿って、ＣＭＰ間隔を５．０ｍとし、ＣＭＰギャザーはＣＭＰを中心として５．０ｍ（測線方向）×５００ｍ（測線と直交方向）の範囲のトレースデータから構成されるようなＣＭＰ編集を実施した（図２−３−１３−１−１，図２−３−１３−１−２，図２−３−１３−１−３，図２−３−１３−１−４。
（５）振幅回復、ミュート
幾何減衰を補償するために、図２−３−９に示す振幅回復テストを行なった。その結果、ウィンドウ長６００ｍｓの振幅調整（ＡＧＣ：Ａｕｔｏｍａｔｉｃ Ｇａｉｎ Ｃｏｎｔｒｏｌ） を実施した。この際、事前に屈折初動部分をミュートにより除去した。ミュートパラメータは以下の通りである（図２−３−１３−２−１，図２−３−１３−２−２）。
最小オフセット ： ０ｍ
スライディング開始時間 ： ０ ｓｅｃ
スライディング速度 ： １９００−２６００ｍ／ｓ（スペースバリアント）
テーパー ： ２００ ｍｓ
（６）デコンボリューション
図２−３−１０に示すテストの結果オペレータ長２４０ｍｓ、ウィンドウ長１０００ｍｓのホワイトニングデコンボリューションを実施した（Ｆｉｇ．１８（３））。
（７）Ｆｌｏａｔｉｎｇ ｄａｔｕｍ ｐｌａｎｅへの静補正
陸上部分のデータに対し改良タイムターム法を用いて表層速度構造を推定した。図２−３−１１にその結果を示す。各ＣＭＰ内の平均標高（ｆｌｏａｔｉｎｇ ｄａｔｕｍ ｐｌａｎｅ）を設定し、この点で受・発震が行なわれるような補正を実施した（図２−３−１４−４）。この際、表層は表層基底層速度で置き換えた。
（８）速度解析
定速度重合法により速度解析を行ない、速度プロファイルを作成した。本測線においては、水平方向の構造変化が激しいため、解析は標準５０ＣＭＰ（２５０ｍ）間隔で実施した。速度解析結果については最終記録断面図上に表示するとともに、巻末にまとめた。速度プロファイルを図２−３−１２に示す。
なお、明瞭な反射イベントの無い深部については区間速度として４５００ｍ／ｓを仮定した。
（９）ＮＭＯ補正及びミュート
速度解析により求まった速度関数を用いてＮＭＯ補正を実施した（図２−３−１４−５）。また、ＮＭＯ後のオフセット距離の大きいトレースで顕著な初動付近の波形歪を除去するために、ＮＭＯ補正により２．０倍以上に波形が伸びるデータを除去した。さらに、以下のパラメータのミュートも適用した。
最小オフセット ： ２００ｍ
スライディング開始時間 ： １００ ｍｓ
スライディング速度 ： １９００−２６００ ｍ／ｓ
テーパー ： ５０ ｍｓ
（１０）振幅調整
ウィンドウ長５０ｍｓのＡＧＣにより振幅を調整した。
（１１）残差静補正
トレース間の相互相関を求め、最大相関値を与えるラグのＣＭＰ内の平均値からのずれを誤差と定義し、その誤差を統計的に発震点、受振点の補正値に分離して補正を行なう残差静補正を適用した（図２−３−１４−６）。
（１２）ＣＭＰ重合
標準６０重合のＣＭＰ重合を行なった。
（１３）バンドパスフィルター
図２−３−１３に示す周波数成分解析により、以下の周波数レンジをもったバンドパスフィルターを、時間が大きくなるにつれて狭帯域になるようにして適用した。
０ − １．０ｓｅｃ １０−８０Ｈｚ
１．０ − １．３ｓｅｃ １０−７０Ｈｚ
１．３ − ５．０ｓｅｃ １０−６０Ｈｚ
（１４）Ｄａｔｕｍ Ｐｌａｎｅへの静補正
ここまでの処理において時刻ゼロは地表面でありＣＭＰ平均標高（ｆｌｏａｔｉｎｇｄａｔｕｍ）となっている。時刻ゼロをある一定標高（ｄａｔｕｍ ｐｌａｎｅ）とするような標高補正を実施した。この際ｄａｔｕｍ ｐｌａｎｅは７０ ｍ Ｍ．Ｓ．Ｌ．とし、補正速度としては表層基底層速度を用いた。
（１５）ＦＤマイグレーション
差分法マイグレーションを実施した。マイグレーション速度は速度解析結果を水平方向に平滑化したものを使用した。これはマイグレーション記録上に表示した。マイグレーション速度は上記の速度値を９０％にしたものを使用した。
（１６）深度変換
時間マイグレーション前後のそれぞれの記録に対し、マイグレーション速度を用いて深度に変換した。深度変換後のサンプル間隔は２ｍである。 〔反射法データ処理結果〕
重合記録、時間マイグレーション記録及びそれぞれに対する深度記録の４種類を、図２−３−１５、図２−３−１６、図２−３−１７、図２−３−１８に表示した。さらに、マイグレーション後の深度記録のカラー表示を図２−３−１９に示す。深度記録はいずれも縦横比１：１である。
図２−３−１５、図２−３−１６、図２−３−１７、図２−３−１８、図２−３−１９に示す最終記録においては、往復走時で約１秒、深度で約１２００ｍ付近まで多数の反射面が鮮明に捉えられている。記録上で特筆されるのは、呉羽山西縁から井田川の区間において、呉羽丘陵の東に位置する牛ｹ首用水付近を軸とする（ＣＭＰ６７５）背斜構造とこの東翼部に位置する断層が明瞭に認められることである。呉羽山東縁には断層は認められない。特に、断層は断層面自体が反射面として捉えられている（地表位置はＣＭＰ５４０〜５５５）。明瞭ではないが、この断層面は地下深部（２０００ｍ付近）まで直線的に延びているように見える。
図２−３−２０には、断層の両側においていくつかの反射面を選び、速度解析結果から求めた区間速度を深度記録上に表示した。図２−３−２１には断層両側での代表的な速度解析例を示す。これらによると、断層を境として両側で速度が劇的に変わっている。西側のﾎﾗｲｿﾞﾝＡは、呉羽丘陵東崖での地表地質から音川期上部（鮮新世）に対応しており、２１００〜２２００ｍ／ｓの区間速度は標準的な値を示している。また、ﾎﾗｲｿﾞﾝＣ−Ｄ間は八尾累層に対応すると考えれば、その区間速度２６００〜２８００ｍ／ｓは中新世の堆積層として一般的な値を示している。一方、断層の東側では浅部から区間速度が２７００ｍ／ｓを越えており、中新世の一般的な地層速度に対応する。このように、断層を境として両側で速度構造が大きく異なることは、図２−３−１１の表層速度構造や図２−３−１２の速度プロファイルにもはっきりと現れている。
当初呉羽丘陵東縁に推定されていた呉羽山断層は反射記録においては認められず、より東方の平野側に存在している。この断層は西に４５゜の傾斜を有し、断層両側での層序対比が不明なため正断層か逆断層か断定はできないが、背斜構造をともなうことから類推して逆断層の可能性が高い。この断層以外には、顕著な断層は認められないが、背斜構造の西側の変曲部において深度９００ｍ〜３５０ｍにわたり反射面の不連続が認められる。ただし、垂直変位はほとんどない。この不連続部は地表までは達しておらず、呉羽山西側に推定されていた、友坂断層の兆候は認められない。