�@　前処理

ａ）　データの転送と編集

磁気テープから処理装置（ＳＵＮ Ｓｐａｒｃ２０）にデータを転送する。次に不必要なショットレコード及びトレースを除く。フィルターテストの結果Ｒ−１，Ｒ−２，Ｒ−３測線のデータには２００Ｈｚ以上のシグナルに乏しかったので１ｍｓｅｃから２ｍｓｅｃにリサンプルした。ＲＳ−１測線のデータも同様の理由から、０．５ｍｓｅｃから１ｍｓｅｃにリサンプルした。

Ｂ）　ジオメトリーの定義

測量によって得られた受振点及び発振点の座標を入力する。各ショットの受振点パターンを定義する。処理ラインを設定しＣＤＰ ビンのサイズを定義する。ビンの測線方向の間隔は、Ｒ−１，Ｒ−２，Ｒ−３測線で５ｍ，ＲＳ−１測線で１ｍとした。

�A　重合前フィルター

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

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

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

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

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

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

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

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

本処理ではスパイキング・デコンボリューションを用いた。

Ｂ）　帯域通過フィルター

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

�B　静補正

図２−３−１−４に静補正の概念図を示す。

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

以下に図２−３−１−４を用いて概説する。

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

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

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

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

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

Ｔ（Ｉ，ｊ，ｋ） ＝ Ｔ（Ｉ，ｊ，ｋ−（Ｒｊ＋Ｓｉ）／Δｔ） （Δｔ はサンプリングレート）

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

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

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

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

本処理では、屈折法的解析の風化層補正を地形補正とともに実施した。

�C　速度解析

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

�D　ＮＭＯ補正及びＣＤＰ重合

図２−３−１−５にＮＭＯ補正及びＣＤＰ重合の概念図を示す。

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

測線長の長いＲ−１・Ｒ−２測線についてオフセット距離が長く深部からの反射波に乏しいので、ＣＤＰの近傍５０ｃｈの観測データを重合して処理を行った。解析精度均質化を図るためＲ−３測線もこれに習った。

図２−３−１−５　水平重合法（ＮＭＯ補正及びＣＤＰ重合）の概念図

�E　残差静補正

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

�F　重合後フィルター

帯域通過フィルター，Ｆ−Ｋフィルター，デコンボリューションフィルター，コヒーレンシーフィルターなどが使用される。

�G　マイグレーション

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

図２−３−１−６ マイグレーション概念図

�H　深度変換

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

図２−３−１−３−２　デ−タ処理フローチャート

図２−３−１−４ 静補正の概念図