ａ）フォーマット変換及びトレース編集

現場で記録された原データ（ＧＤＡＰＳ−４ ＳＥＧＹ Ｆｏｒｍａｔ）からＳＳＶ８ Ｆｏｒｍａｔ（ＪＧＩ Ｉｎｔｅｒｎａｌ Ｆｏｒｍａｔ）への変換を行った。発震記録上で不良トレースを選択し、処理対象から除外した。

ｂ）トレースヘッダーへの測線情報の入力

各トレースのヘッダーに対して、発震点・受振点の座標・標高を入力した。

ｃ）表層補正解析（Ｒｅｆｒａｃｔｉｏｎ Ａｎａｌｙｓｉｓ）

全ての現場原記録から屈折初動走時を読みとり、その値から発震点・受振点におけるタイムターム、及び表層基底層速度を未知数とするインバージョン（改良タイムターム法）を行い、表層の構造を求めた。

インバージョンの結果得られた各発震点・受振点のタイムターム値を図３−２−４上段に示す。赤色の点が発震点、青色の点が受振点のタイムターム値を表している。横軸は測線の受信点番号である。同図中段には表層速度及び表層基底層速度、下段にはタイムターム値、速度情報から計算した深度構造を示す。

これらの結果は、表層第１層の厚さの変化及び標高変化に対する走時変化の補正（表層静補正）を行う際に使用される。

ｄ）最小位相変換（Ｍｉｎｉｍｕｍ Ｐｈａｓｅ Ｃｏｎｖｅｒｓｉｏｎ）

相互相関処理後の震源波形は、ゼロ位相型となっているが、後述のデコンボリューション処理によりゼロ位相波形では波形の歪が生じてしまう。これを避けるために、震源波形をゼロ位相から最小位相変換するオペレータを設計し、それを原記録に適用する最小位相変換を行った。

ｅ）共通反射点編集（ＣＤＰ編集）

全記録の発震点・受振点座標を用いて、反射点（発震点と受振点の中点座標）の分布図を作成し、反射点の分布が密な位置を選び、重合測線（ＣＤＰ測線）を設定した。図３−２−５に反射点分布図及び重合測線を示す。

・ＣＤＰ間隔 ： １２．５ｍ（受振点間隔２５ｍの１／２）

・総ＣＤＰ数 ： １０９６

ｆ）振幅補償（Ｇａｉｎ Ｒｅｃｏｖｅｒｙ）

反射波の減衰効果を補正するため、ｔ１．２（時間の１．２乗）型の振幅補償に加え、ゲート長 ６００ ｍｓｅｃの自動振幅調整（ＡＡＣ）を行った。

ｇ）デコンボリューション（Ｄｅｃｏｎｖｏｌｕｔｉｏｎ）

発震点・受振点の特性の相違を補正し、分解能の高いデータを得るためのデコンボリューション（ホワイトニング デコンボリューション）を行った。テストの結果、以下のパラメータを使用した。

・ゲート長　　　　　： １０２４ ｍｓｅｃ

・オペレータ長　　： ９６ ｍｓｅｃ

・ホワイトノイズ　　： ０．５ ％

・予測距離　　　　　： ８ ｍｓｅｃ

ｈ）浮動基準面に対する静補正

表層補正解析ｃ）の結果を用いて、浮動基準面に対する静補正を実施した。浮動基準面として各ＣＤＰ毎に計算される平均標高を用いた。平均標高の計算には各ＣＤＰ内のオフセット範囲０〜５００ｍのトレースを使用した。浮動基準面に不連続を生じないように平滑化を行っている。

ｉ）トレースバランシング

以下の自動振幅調整を重合前に適用した。

・ゲート長 ： １００ ｍｓｅｃ

ｊ）速度解析（Ｖｅｌｏｃｉｔｙ Ａｎａｌｙｓｉｓ）

定速度重合法（Ｃｏｎｓｔａｎｔ Ｖｅｌｏｃｉｔｙ Ｓｔａｃｋ）を用いて速度解析を行った。速度解析例を図３−２−６−１、図３−２−６−２、図３−２−６−３、図３−２−６−４、図３−２−６−５、図３−２−６−６、図３−２−６−７、図３−２−６−８、図３−２−６−９、図３−２−６−１０に示す。

ｋ）ＮＭＯ補正（ＮＭＯ Ｃｏｒｒｅｃｔｉｏｎ）

各速度解析点で決定した速度関数（Ｔｏ，Ｖ）を測線方向に内外挿することによって得られた重合速度プロファイル（図３−２−７）を用いてＮＭＯ補正を行った。補正に伴う波形のひずみを抑制するために最大の波形の伸びを２．２倍とした。

ｌ）ミュート

ＮＭＯ補正後に残存する屈折波やＮＭＯによる過伸長部を除去することを目的にミュートを実施した。

ｍ）残差静補正（Ｒｅｓｉｄｕａｌ Ｓｔａｔｉｃｓ）

既述したｈ）で実施した浮動基準面に対する静補正では、比較的長周期の補正値は精度良く補正されるが、短周期の受振点・発震点固有の補正は不十分である。そこで、ＮＭＯ補正後のデータに見られる連続性の良い反射波を利用して残差静補正を行った。補正の際に許される相関関数のラグ値は ８ ｍｓｅｃとした。

ｎ）ＣＤＰ重合（ＣＤＰ Ｓｔａｃｋ）

ＮＭＯ補正、残差静補正終了後、各ＣＤＰ内の反射波走時は、同一時間に並び、屈折波・表面波・ノイズ等は同一走時とならない。そこで、これらを足し合わせる（重合する）ことで、Ｓ／Ｎの良い反射記録が得られる。重合には、オフセット範囲を０〜３０００ ｍのトレースを使用した。

ｏ）Ｆ−ＸＹ予測フィルター（Ｆ−ＸＹ Ｐｒｅｄｉｃｔｉｏｎ Ｆｉｌｔｅｒ）

重合断面上で連続性のあるイベントを強調し、ランダムノイズを抑制するために、Ｆ−ＸＹ予測フィルターを適用した。

ｐ）最終断面図（Ｆｉｌｔｅｒｅｄ Ｓｔａｃｋ）

重合断面に以下のバンドパスフィルターを施し、フィルター後重合断面とした。

・時間 ０ 〜 ２５００ ｍｓ ４／６ 〜 ４５／５０ Ｈｚ

ｑ）時間マイグレーション（Ｔｉｍｅ Ｍｉｇｒａｔｉｏｎ）

重合断面図上では、反射波は、各ＣＤＰ位置から反射面までの往復垂直走時がそのＣＤＰ位置に表現されている。したがって、傾斜した反射面に対しては、重合断面図上の反射面の傾斜及び位置が、真の構造から若干ずれてくる。これを補正し、各ＣＤＰ直下の構造形態を表す様にする処理がマイグレーション処理である。

本測線では差分法マイグレーション処理を行い、マイグレーション速度は平滑化した重合速度（１００％）を用いた。

ｒ）深度変換（Ｄｅｐｔｈ Ｃｏｎｖｅｒｓｉｏｎ）

マイグレーション後の記録に対し、平滑化した速度関数を用いて、時間軸から深度軸へ変換を行った。深度断面の基準面は平均海水面とした。

以上の処理内容のうち、次の７項目についてはデータ処理の品質管理のために、中間段階の記録として、重合数１の断面（１００％セクション）を作成している（巻末資料７）。

（ａ）最小位相変換

（ｂ）振幅補償

（ｃ）デコンボリューション

（ｄ）浮動基準面に対する静補正

（ｅ）ＮＭＯ補正

（ｆ）ＮＭＯ補正後のミュート

（ｇ）残差静補正