デ−タ処理は、以下にまとめるように、基本処理と解釈のための処理に大別される。

●基本処理

ａ．初期編集‥‥‥デ−タ整理・ＣＭＰ編集．
ｂ．静補正 ‥‥‥表層静補正・標高補正・残留静補正．
ｃ．波形処理‥‥‥プレフィルタ−・位相振幅補償・デコンボリュ−ション．
ｄ．速度解析‥‥‥重合のための速度を求める．
ｅ．重　合

これらは時間断面の作成を目的とする処理である。この過程では記録の質が良好な限り、あまり人為的な処理操作は含まれない。

●解釈のための処理

ａ．マイグレ−ション‥‥地層の傾斜・回折による見かけの断面を補正する．
ｂ．深度変換‥‥‥‥‥‥時間断面を深度断面に変換する．

これらは時間軸上で与えられた見かけ上の記録から、真の地質構造を推定するための処理である。この処理においては、与えた速度分布により結果が大きく異なることがあり、できるだけ正確な速度分布を与える必要がある。
各探査測線におけるデ−タ処理の手順およびパラメ−タ−を表３−３−１（１）〜（３）に、またデ−タ処理のブロックダイアグラムを図３−３−１に示す。

表３−３−１（１）　デ−タ処理手順および処理諸元一覧表[岡部Ｓ波Ａ測線]

表３−３−１（２）　デ−タ処理手順および処理諸元一覧表[岡部Ｓ波Ｂ測線]

表３−３−１（３）　デ−タ処理手順および処理諸元一覧表[本庄Ｓ波測線]

（１）初期編集

デ−タ整理…ノイズ除去、不良デ−タの除去、測量デ−タの整理等を行う。

ＣＭＰ編集…測定は通常道路上で行い、障害物を避けて実施するため、ＣＭＰ（反射中点）は一般に平面上に分布する。重合用測線はこれらのＣＭＰの分布を最適に通る折れ線として設定し、この測線沿いにＣＭＰ編集を行った。なお、重合測線はオフセット距離（震源−受震器間距離）が比較的短いものを考慮して決定した。

（２）静補正

一般的には屈折法により表層をはぎ取る方法が用いられるが、特に、ミラ−ジ的な速度変化を示すような速度構造地盤では、必ずしも精度の高い補正値を得られるとは限らない。ここでは、「屈折波を用いたトモグラフィ−」により静補正値を算出し、表層に起因する乱れを補正した。この処理の手順は次のとおりである。

ａ．観測波形よりＰ波の初動走時を読み取る。
ｂ．差分格子点に適当な初期速度分布値を与える。
ｃ．アイコナ−ル法によりある震源点で起震した場合の各格子点の初動走時を計算する。
ｄ．初動走時分布をもとに波線を求める。
ｅ．各波線の観測走時と計算走時の比を修正係数とし、波線周辺の格子点に記憶する。
ｆ．ｃ〜ｅの操作を全震源点について行う。
ｇ．格子に配られた修正係数をもとに新たな速度分布を算出する。
ｈ．ｃ〜ｇの操作を収束するまで繰り返す。

ＣＭＰアンサンブル内での標高静補正・・・ＮＭＯ補正前に、各アンサンブルごとにその平均標高までの標高差補正を行った。なお、補正速度は各探査測線によって異なるが、１４０〜２００ｍ／ｓｅｃを用いた。

重合後標高補正・・・マイグレ−ション、深度変換後に各ＣＭＰの平均から基準標高までの標高補正を行った（基準標高はT．P．５０ｍ）。また、時間断面図のプロットの際も、地表平均標高（ｆｌｏａｔ ｄａｔｕｍ）から基準標高までの速度を各探査測線に応じて１４０〜２００ｍ／ｓｅｃに仮定して、標高補正を実施した。

図３−３−２（１）〜（３）に各探査測線のトモグラフィー解析結果を示す。

図３−３−２（１）　屈折波トモグラフィー解析による表層速度構造[岡部Ｓ波Ａ測線]

図３−３−２（２）　屈折波トモグラフィー解析による表層速度構造[岡部Ｓ波Ｂ測線]

図３−３−２（３）　屈折波トモグラフィー解析による表層速度構造[本庄Ｓ波測線]

（３）波形処理

初期処理

波形処理のうち最も重要な処理は、パルスの短縮、短周期の多重反射の除去、スペクトルの平滑化等を目的に実施するデコンボリュ−ションである。この処理を良好にするため、次の前処理を実施した。
　 ・プレフィルタリング…再帰型のフィルタ−を用い、位相特性は次の位相補償で併せてミニマムフェ−ズに直した。
　 ・位相補償…デコンボリュ−ションが有効に働くためには、トレ−スがミニマムフェ−ズ特性であることが条件の一つである。測定系でもっともこの条件を満たさないものは、探鉱機のフィルタ−の位相特性である。これをミニマム位相特性に戻すフィルタ−を設計し、補償を行った。
　 ・振幅補償…次の２段階に分けて実施した。
　　 ａ．全トレースよりオフセット距離（震源−受震器間距離）別に振幅の時間減衰特性を統計的に求め、この特性の逆数で振幅補償を行った。
　　 ｂ．次に各トレース別に、所望のゲ−ト幅で平均振幅の時間変化を求め、振幅補償を行った（ＡＡＣ）。

デコンボリュ−ション

自己相関演算のゲ−ト長 ７００ｍｓｅｃ､フィルタ−長 ８０ｍｓｅｃ､ホワイトノイズ３％のタイムバリアント型プレディクティブ・デコンボリュ−ションを用いた。

（４）ＮＭＯ補正と速度解析

ＮＭＯ補正

ＮＭＯ補正はオフセットの違いによる走時のずれをゼロオフセットの記録に補正するものである。通常、直線の波線を図３−３−３に示すように仮定した下記の方法で行う。

図３−３−３　ＮＭＯ補正における波線の考え方

地下構造が水平２層でオフセットがＸのとき、第１層下端から反射して受震点に到達する波の走時は次式で表される。

　　　　　　　　　　　 Ｘ２
　ｔ２　＝ ４ｔ０２ ＋───
　　　　　　　　　　　 Ｖ２

　ただし、　　　　　　　　　　　　Ｖ；波の伝播速度
　　　　　　　　　 ２Ｚ 　　　　　Ｘ；震源点と受信点の距離
　　　　　　ｔ０＝───
　　　　　　　　　　Ｖ　　　　　ｔ０；ゼロオフセットの走時

　多層構造の場合も、オフセットＸに比べて反射面深度が十分に大きければ、

　　　　　　　　　　　　　　　　　　 Ｘ２
　　　　Ｔｎ（Ｘ）２ = Tn（０）２ ＋───
　　　　　　　　　　　　　　　　　　ＶＲ２
　と近似できる。

　ここで、ＶＲはＲＭＳ速度と呼ばれるものであり、次式のように定義される。

　　　　　　　ｎ
　　　　　　　 Σ Ｖi２△tｉ ここで Ｖｉ ；第ｉ層の速度
　　　　　　　i=1
　　　ＶＲ２ =─────── △tｉ；第ｉ層の鉛直往復走時
　　　　　　　　Ｔｎ（０）

速度解析

速度解析で得られる速度値は重合速度と呼ばれ、地層がほぼ水平の場合には近似的にＲＭＳ速度に等しいと見なされている。速度解析は次の手順で行った。

ａ．推定される重合速度の範囲のなかで、１２０種の速度を等分に仮定する。
ｂ．各速度でＮＭＯ補正を行い、オフセットによらず反射波が同じ時刻に並ぶ速度と時間を求める。
ｃ．図から直接速度を読み取ることもできるが、労力および精度に問題があるため、ＣＭＰアンサンブルのデ−タを定速度でＮＭＯ補正した後に、
イ．狭い時間ゲ−ト内でのトレ−ス間の相関（センブランス）を求める（速度スペクトル法）
ロ．全てのトレ−スを加算して１本のトレ−スとする（Ｃ．Ｖ．Ｓ．法、ｃｏｎｓｔａｎｔ ｖｅｌｏｃｉｔｙ ｓｃａｎ）などの方法で整理する。
ｄ．実際には全デ−タを用いた場合には、イ．，ロ．の方法でもＣＭＰの数だけ図が得られるため、読み取りを自動的に行い、これを整理して重合速度を決定する。

以上のデータ処理を行った結果は、デコンボリューションテスト結果、速度解析結果、および１００％断面図（１重合断面図）として各測線についてまとめて巻末資料に添付する。

（５）重合・マイグレ−ション・深度変換

上記した各種の解析結果をもとに、重合・マイグレーション・深度変換などのデータ処理をへて時間断面図と深度断面図が作成されるが、それぞれの結果は次の３．４探査結果に示す。以下に各処理過程の概要をまとめて示す。

ミュ−ト・重合
速度解析で得た速度分布を用いてＮＭＯ補正し、ミュ−ト後、重合を行った。

マイグレ−ション
重合後時間断面を水平方向の速度分布が均一なストレッチ断面（深度方向には速度変化を認める）に変換した後、波動場補外法（位相移動によるイメ−ジングと下方接続）によるＦＫマイグレ−ションを行った。

深度変換
マイグレ−ション後の時間断面を、マイグレ−ションに用いたＲＭＳ速度より求めた平均区間速度を用いて深度変換した。

零位相表示断面
位相デ−タとは、各トレ−スを複素変換したのち、振幅情報を無視し位相情報のみを取り出したものである。この位相デ−タのうち零位相に近い部分を黒線で示す。成層構造をなさない場所では、位相が乱れるため破砕部の把握に役立つと思われる。また、振幅が相対的に小さい場所の情報も得られる。

相対振幅強度＋零位相表示断面
相対振幅強度表示とは、深度断面の振幅をその大きさに応じた色で表示したものである。通常の表示では、大きい振幅は隣のトレ−スの上に描かれるため細部が不明瞭となるが、この表示では細部が明瞭となり、また負の振幅も情報として得られる。この断面では、上記の零位相表示と重ね合わせて図化している。