なお、本処理には反射法探査データ処理ソフトＰｒｏＭＡＸ を使用した。

（１）　前処理

（ａ）　データの転送と編集

８ｍｍ磁気テープから処理装置（ワークステーションＳＵＮ Ｕｌｔｒａ６０）にデータを転送し、不必要なデータを取り除く。

（ｂ）　共通反射点編集

発震点・受振点の組み合わせと、測量によって得られた各座標を用いて、各データのＣＭＰ（発震点−受振点の中点）の位置を計算し、ショット記録を共通反射点記録に編集した。屈曲した測線においては反射点が２次元的に分布するので、反射点の平均的な位置に仮想の測線（重合測線）を設定し、ばらついた反射点をあたかもその測線上で反射しているように集めた。Ｐ１、Ｐ２測線の重合測線の位置を、図２−１−２、図２−１−４に示す。

（２）　重合前フィルタ

重合前には、できるだけ反射波を強振幅かつインパルス（ｉｍｐｕｌｓｅ：瞬間的に振幅を有する波形）に近い波形に変換し、反射波の周波数帯域以外のノイズを弱めたり除去したりするため、デコンボリューション・フィルタ（下記）や帯域通過フィルタなどが適用される。

（ａ）　デコンボリューション・フィルタ＊

このフィルタは震源波形、地層特性など反射地震記録に、コンボリューションの関係で含まれている基本波形をインパルスあるいはそれに近い波形に変換するフィルタの一種である。デコンボリューション・フィルタの具体的な効果は以下のとおりである。

・　様々な周波数成分をもつ間延びした反射波を，インパルスに近い（高周波かつ分解能の高い）波に変換する。

・　主として浅部の影響による重複反射波を除去または弱め、独立した反射波に変換する。

＊（注）線形フィルターをかけることをコンボリューション（ｃｏｎｖｏｌｕｔｉｏｎ）、その逆操作をデコンボリューション（ｄｅｃｏｎｖｏｌｕｔｉｏｎ）と言う。

（ｂ）　帯域通過フィルタ

信号である反射波とノイズである他の波との周波数帯域が異なっている場合には、反射波の帯域のみを通す帯域通過フィルタをかけることにより、Ｓ／Ｎ比（ｓｉｇｎａｌ − ｎｏｉｓｅ ｒａｔｉｏ）の向上が期待できる。そのためには周波数領域でフィルタを設計し、それをフーリエ変換して時間領域のフィルタ・オペレータを求め、反射記録にコンボリューションする。具体的には何種類かの帯域のフィルタ・テストを実施して、最適なものを適用する。

（３）　静補正

発震点および受振点の標高差、表層付近での弾性波の速度差による反射波の遅速、表層の厚さの変化による反射波の遅速などを補正する（前一者は標高（地形）補正、後二者は表層（風化層）補正と呼ばれる）。これを実施しない場合には構造解釈を誤ったり、水平重合法において反射波そのものを消去することがある。

具体的には、各発震点・受振点の標高の移動平均を仮想の基準面とし、あたかもその基準面上で測定が行われたかのように各発震点、各受振点の記録を上下させる。

なお、本処理での最終の基準面は、Ｐ１、Ｐ２測線とも標高１００ｍとした。

表層（風化層）補正は、ショット記録の初動から屈折法の解析を行って求められる。

本処理では、全ショット記録を用いた屈折トモグラフィーによる表層解析を行い、補正量を求めている。屈折トモグラフィー解析により得られたＰ１、Ｐ２測線の表層付近の速度分布を図２−１−６、図２−１−７に示す。

（４）　速度解析

ＣＭＰ重合に用いる重合速度分布を求める処理である。具体的には、様々な重合速度でＮＭＯ補正（後述）、水平重合を施し、そのうち反射波の振幅強度が最大となったときの速度を最適な重合速度として採用する。速度解析はこのような最適な重合速度を時間の関数として選んでいくものである。速度解析は、測線に沿ったある間隔で行われ、その間の速度は内挿により求められる。ここで求められる速度は重合速度と呼ばれ、地層が水平な場合は ＲＭＳ速度に近似的に一致するが、地層が傾斜している場合にはＲＭＳ速度より速くなり、また傾斜が急であるほど速い値となる。速度解析の方法としては、定速度走査法と定速度重合法があるが、ここでは両方を用いた。

（５）　ＮＭＯ補正およびＣＭＰ重合

ＣＭＰアンサンブル（各ＣＭＰ毎に対応するデータを集めたもの）に対して発震点・受振点間距離の違いによる反射波の到達時間（走時）の遅れを補正し、発震点と受振点が同一位置にある場合（ゼロオフセット）の走時に合わせる操作が ＮＭＯ補正である（ＮＭＯ：ｎｏｒｍａｌ ｍｏｖｅ−ｏｕｔ ｃｏｒｒｅｃｔｉｏｎ）。ＮＭＯ補正を行うにはあらかじめ地下の速度分布を設定する必要があり、通常は速度解析により得られた速度分布を用いる。このＮＭＯ補正後のＣＭＰアンサンブルを足し合わせて、反射波を強調する操作が水平重合（ＣＭＰ重合）である。共通反射点Ｐでは、Ａ−Ｐ−ａ，Ｂ−Ｐ−ｂ，Ｃ−Ｐ−ｃ，Ｄ−Ｐ−ｄ，Ｅ−Ｐ−ｅ およびＦ−Ｐ−ｆ のそれぞれ異なった経路をもつ反射記録（�T） が得られる。Ｐ点からの反射波の走時は、水平距離Ｘの増加とともに遅くなっている。次に速度解析で得られた速度を用いてｐ−Ｐ−ｐ 経路の仮想の記録に合わせる（ＮＭＯ補正）と、Ｐ点からの反射波はいずれの経路のものもほぼ同じ走時の記録（�U） が得られる。さらにこれらを足し合わせる（水平重合）と、ランダムなノイズや反射波以外の波は打ち消しあって相対的に弱くなり、逆にＰ点からの反射波は加算され、強調された記録（�V）が得られる。

（６）　残差静補正

高度補正や表層補正を施した後でも，初動屈折波と反射波の経路の違いによる時間の不規則性や、表層補正で行う２層構造の仮定を採用したことによる局所的な速度の異常に起因するものは完全には補正されず、ＣＭＰアンサンブル内での同一反射の到達時間は一定ではないのが普通である。水平重合反射法探査においては、最適なＣＭＰアンサンブル群が得られるように統計的処理を施してこの時間差を補正し、各発震点および受振点における２次補正値を求める。

（７）　重合後フィルタ

重合後のフィルタとしては、重合断面をより地下構造を反映した断面にすることを目的に用いられ、通常、帯域通過フィルタ、Ｆ−Ｋフィルタ、デコンボリューション・フィルタなどが使用される。

（８）　マイグレーション

傾斜した地層から反射した波でも各データの直下に表示されるため、重合記録上の反射面の傾斜や位置は地下の実際の反射面とは異なってくる。記録上での見かけの傾斜や位置を移動させて、真の傾斜および位置に復元することをマイグレーション処理という。マイグレーションの方法としては波動方程式マイグレーションが最も一般的で、計算手法により差分法、Ｆ−Ｋ法などがある。

ここでは、地下の速度分布が不明確な場合やノイズの多い記録に有効な差分法を用いた。

マイグレーションに使用する速度としては、通常、水平多層構造の場合には重合速度をそのまま用いるが、反射波が傾斜している場合に重合速度をそのまま用いるとオーバーマイグレーションの可能性が出てくる。

（９）　深度変換

（８）までの処理によって得られる断面図の縦軸は時間である。地下構造を解釈するためには時間軸を深度軸に変換する必要がある。地下の速度分布を仮定して縦軸を深度に変換する。深度変換に使用する速度には、通常、重合速度から得られる区間速度が用いられるが、ここでは地層の傾斜等を考慮した。

（１０） 断面図作成

以上の処理を施した両測線のマイグレーション前の時間断面図、深度断面図、マイグレーション後の時間断面図、深度断面図およびカラー表示深度断面図を作成した（巻末に掲載）。

本データ処理に用いられたパラメータを、表２−１−５にまとめて示す。

表２−１−５ 浅層反射法探査のデータ処理