（１）前処理

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

　磁気テープから処理装置（ＳＵＮ Ｓｐａｒｃ２０）にデータを転送する。次に不必要なショットレコード及びトレースを除く。

（ｂ）ジオメトリーの定義

　測量によって得られた受振点及び発振点の座標を入力する。各ショットの受振点パターンを定義する。処理ラインを設定しＣＤＰ ビンのサイズを定義する。

（２）重合前フィルター

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

　このフィルターは振源波形、地層特性等反射地震記録に、コンボリュー ションの関係で含まれている基本波形をインパルスに短縮するフィルタ ーの一種である。入力波形は以下の条件を有すると仮定する。

反射係数時系列は無規則（ランダム）で定常である。

基本波形（基本ウェーブレット）は最小位相型である。

　一般に、反射記録は反射係数時系列に基本波形がコンボリューションされたものと考えられる。基本波形が既知であれば、これをインパルスに変換するフィルター（これを逆フィルターと呼ぶ）は正確に設計できる。しかし、反射記録の性質上基本波形を正確に知ることは非常に困難であるため、デコンボリューション・フィルターでは、上記最小位相型の仮定を設けることにより統計的に処理されている。

デコンボリューション・フィルターの具体的な効果は以下のとおりである。

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

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

本処理ではＳｐｉｋｉｎｇ Ｄｅｃｏｎｖｏｌｕｔｉｏｎ Ｆｉｌｔｅｒを用いた。

（ｂ）帯域通過フィルター

　信号である反射波とノイズである他の振動との周波数帯域が異なっている場合には、反射波の帯域のみを通す帯域通過フィルターをかけることにより、Ｓ／Ｎ比の向上が期待できる。そのためには周波数領域でフィルターを設計し、それをフーリエ変換して時間領域のフィルタ・オペレータを求め、地震記録にコンボリューションする。

（３）静補正

　表層付近での弾性波の速度差による反射波の遅速、表層の厚さの変化による反射波の遅速、発振点及び受振点の標高差等を補正する。これを実施しない場合には構造解釈に誤りをおかしたり、水平重合法において反射波そのものを破壊することがある。 具体的には、ある基準面（Ｄａｔｕｍ Ｐｌａｎｅ）を設けて、あたかもその基準面上で測定が行われたかのように各発振点及び各受振点の記録を上下する。

　以下にモデルを用いて概説する。図２−３−６ （ａ） のモデルに対する反射波の走時は（ｂ） の○印のようになる。（ａ） において実際の標高（○印）の移動平均をとりその標高値を×印で表した。受振点ｊ の実際の標高をＨｊ、移動平均標高をＡｊ、風化層以下の伝幡速度をＶ とすると、受振点ｊ の補正値Ｒｊ は、

　Ｒｊ ＝ （Ｈｊ − Ａｊ）／Ｖ

となり、発振点Ｉ の補正値Ｓｉ は、

　Ｓｉ ＝ （Ｈｉ − Ａｉ）／Ｖ

となる。各トレースについて、

　Ｔ（Ｉ，ｊ，ｋ） ＝ Ｔ（Ｉ，ｊ，ｋ−（Ｒｊ＋Ｓｉ）／Δｔ） （Δｔ はｻﾝﾌﾟﾘﾝｸﾞﾚｰﾄ）

となるように時間をずらす。具体的には（ｂ） の×印のようになる。（ｃ） は（ｂ）の○印（標高補正行わず）の反射波を重合したものであり、（ｄ） は×印（標高補正実施後）の反射波を重合したものである。単純な水平多層構造であるにもかかわらず、（ｃ） では複雑な地下構造であるかのような結果となっているが、（ｄ）ではほぼモデルの地下構造を再現している。

　ここでは地形の補正のみについて説明したが風化層補正の場合も同様に処理できる。具体的には風化層の厚さによる各点の走時差を Ｗｓｉ，Ｗｒｊ とすると上式は、

　Ｔ（Ｉ，ｊ，ｋ） ＝ Ｔ（Ｉ，ｊ，ｋ−（Ｒｊ＋Ｓｉ＋Ｗｓｉ＋Ｗｒｊ）／Δｔ）

となる。補正値Ｗｓｉ、Ｗｒｊ は、各トレースの初動から屈折法的解析を行って求められる。

　本処理では，取得されたデータの初動の質は良くないが、比較的読みとり可能な振源近くのデータ（オフセット距離：３００ｍ以下）を用いて、屈折法的解析の風化層補正を行った。

（４）速度解析

　ＣＤＰ重合に用いる重合速度分布を求める処理である。速度解析の方法としては、定速度走査法と定速度重合法があり、本処理では両者を用いた。ここで求める速度は、重合速度と呼ばれ、水平多層構造を仮定した場合は ＲＭＳ 速度に一致する。

（５）ＮＭＯ補正及びＣＤＰ重合

　ＣＤＰアンサンブル（各 ＣＤＰ［発振点−受振点の中点］毎に対応するトレースを集めたもの）を発振点・受振点間距離の違いによる反射波の到達時間の遅れを補正し、発振点と受振点が同じ場合（ゼロオフセット）の時間に合わせる操作が ＮＭＯ補正である。ＮＭＯ補正を行うにはあらかじめ地下の速度分布を設定する必要があり，通常は速度解析により得られた速度分布を用いる。このＮＭＯ補正後のＣＤＰアンサンブルを足し合わせて、反射シグナルを強調する操作が水平重合（ＣＤＰ重合）である。図２−３−７は、６重合の場合のＮＭＯ補正と、ＣＤＰ重合（水平重合）処理の概念図である。共通反射点Ｐでは、Ａ−Ｐ−ａ，　Ｂ−Ｐ−ｂ，Ｃ−Ｐ−ｃ，Ｄ−Ｐ−ｄ，Ｅ−Ｐ−ｅ及びＦ−Ｐ−ｆの各々異なった経路をもつ反射記録（Ｉ） が得られる。Ｐ点からの反射波の走時は、水平距離Ｘの増加とともに遅くなっている。次に速度解析で得られた速度を用いてｐ−Ｐ−ｐ の経路の仮想の記録に合わせる（ＮＭＯ補正）と，Ｐ点からの反射波はいずれの経路のものもほぼ同じ走時の記録（ｉｉ） が得られる。さらにこれらを足し合わせる（水平重合）と、ランダムなノイズは打ち消しあい、反射波以外のシグナルは相対的に弱くなり、逆にＰ点からの反射波は加算されて反射波が強調された記録（ｉｉｉ） が得られる。

（６）残差静補正

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

（７）重合後フィルター

　帯域通過フィルター，Ｆ−Ｋフィルター，デコンボリューション・フィルター，コヒーレンシー・フィルターなどが使用される。本処理では、予測デコンボリューション・フィルター及び帯域通過フィルターを用いた。

（８）マイグレーション

　傾斜している反射面を空間的に正しい位置に戻す操作である（図２−３−８参照）。マイグレーションの具体的な方法としては、ディフラクション・マイグレーション、反射事象のピッキング・マイグレーション、波動方程式マイグレーション、Ｆ−Ｋ領域マイグレーションがあるが、本処理では波動方程式マイグレーションを用いた。

（９）深度変換

　ＮＭＯ補正で用いた重合速度を用いて得られた、反射断面の縦軸を時間から深度に変換する。