主な機器の性能と今回の処理での用途は、次のとおりである。

　種　類　　　　　　　名　称　　　　　　性能　　　　　用　途

コンピュータ　　　ＳＧＩ Ｐｏｗｅｒ−Ｃｈａｌｌｅｎｇｅ　６００Ｍｆｌｏｐｓ 　　データ処理

　　〃　　　　　　ＳＵＮ　Ｓｐａｒｃ−ＩＰＸ　　 　 ２０Ｍｉｐｓ　　　 静補正

　　〃　　　　　　Ｃｏｎｖｅｘ　Ｃ１−ＸＰ　　　 　４０Ｍｆｌｏｐｓ　　 データ転送

　プロッタ　　　　Ｖｅｒｓａｔｅｃ　７４２４　　　　白黒　　　　 プロット

　　〃　　　　　　ＨＰ　６５０Ｃ　　　　　　 カラー プロット

ソフトウェア　　　ＣＳＤ　Ｆｏｃｕｓ−２Ｄ　　　　　　　　　　　データ処理

今回のデータ処理で採用した手順は、図２−２−１２に示したとおりである。以下では、図２−２−１２に赤の星印で示した処理の詳細について述べる。

（１） 振幅調整

観測される地震波は時間と共に減衰してゆく。この原因は主として波頭面積の増大（エネルギーの幾何学的発散）にある。

振幅調整は、この減衰を補償するために行なわれるが、反射波の振幅を規格化された入射波に対する振幅（反射強度）に近づけることや、各種のフィルタ処理の効果を高めるために反射時系列の定常性を増大させることをも目的とする。

図２−２−１３、図２−２−１４に、振幅調整前後のショット・ギャザーの例を示す。

（２） デコンボリューション（重合前）

振源から発せられた地震パルスは、地層の非弾性や観測システムのフィルタ作用によって歪められ、分解能の悪い波形となる。また反射トレースには層内マルチプルが混入するが、この効果も一種のフィルタ作用である。

デコンボリューションは、これらの作用を補償するために用いる逆フィルタである。

デコンボリューションの理論は地震波の統計的な性質にいくつかの仮定をして組み立てられているが、実際の統計的な性質が未知であることから、適用にあたっては試行錯誤的なテストによって決定するパラメタの選定が重要である。

今回は通常の（ウィナー・レビンソン型）予測デコンボリューションを、オペレータ長３０ｍｓ、予測長２０ｍｓとして用いた。設計及び適用ゲートは、それぞれ５０〜６００ｍｓ、及び全長とした。図２−２−１５にデコンボリューション後のショット・ギャザーの例を示す。

（３） 周波数フィルタ（重合前）

反射トレースに含まれる各種雑音のうち、地震パルスと周波数帯域を異にする成分は周波数フィルタによって容易に除去できる。今回はテストによって、重合前フィルタの通過帯域を２０〜１００Ｈｚと決め適用した。図２−２−１６に周波数フィルタ適用後のショット・ギャザーの例を示す。

（４） 速度フィルタ

反射トレースに含まれる各種雑音のうち、地震パルスと周波数帯域が重なるものでも、反射波と軌跡の傾斜（見掛け速度）を異にするコヒーレントな雑音は、速度差を利用して除去できる。このために使用する２次元フィルタを総称して速度フィルタという。

今回は混入の認められたＳＶ波や音波を除去する目的で、時間−距離領域の速度フィルタを使用した。除去した傾斜は、トレース当たり２．５〜６．０ｍｓの範囲である。図２−２−１７に速度フィルタ適用後のショット・ギャザーの例を示す。

（５） 静補正計算（３Ｄ版）、静補正

観測点の標高の変化や、表層の速度や厚さの変化は、反射パルスに（観測における展開長＝スパンと比較して）長周期及び短周期の位相ずれを与える。これらはいずれも、ＣＤＰ重合法と重合後マイグレーションを基本とするイメージング法では補償できない。そこで浅部の速度構造を別途に解析し、それをもとに（一般に表層の下位に設定する）基準面（データム・レベル）での仮想的な観測走時に変換することを行なう。これを静補正という。

今回は読みとった初動走時をもとに、測線の屈曲を考慮して３次元のトモグラフィックなアルゴリズムを採用し、タイム・ターム法による解析を行なった。データム・レベルは、次のとおりとした（間は直線的に内挿）。

ＳＰ（距離程）　　　　データム・レベル

　 　０ｍ　　　　　　　　　ＧＬ ＋１０ｍ

　１４５５ｍ ＋２０ｍ

２４５５ｍ ＋４５ｍ

解析結果を図２−２−１８に示す。また静補正を施さない場合と施した場合の重合断面の一部を、図２−２−１９（ａ）、図２−２−１９（ｂ）に示す。

（６） ＤＭＯ

ＣＤＰ重合法（より正確にはＣＭＰ重合法）においては、速度構造の設定やイメージング操作が地下座標（ＣＤＰ）系ではなく、地表座標（ＣＭＰ）系で行なわれる。その結果、傾斜した反射面上の各反射点には、ＣＤＰ毎に別々の座標が割り振られることになる。そしてこの不具合は、イメージング速度の見掛けの増大となって現れる。従って例えば不整合部のような、同一の地下座標に異なる傾斜の反射面が混在するような場所では、各々に対して適切な見掛け速度を与えることができないため、各々を同時にイメージングすることができないという不都合が生ずる。

ＤＭＯ補正は、ＣＭＰ重合法と地下座標系でのイメージング法（重合前マイグレーション）の中間的な手法であり、各反射点と地表座標を１対１に対応させることによって、地表座標系を用いる通常の手続きの範囲内で不具合を除去する手法である。今回の調査地には顕著な不整合構造が認められたので、コモンオフセット・ギャザーでのキルヒホッフ積分法によるＤＭＯを適用した。静補正とＤＭＯを施した重合断面の一部を、図２−２−１９（ｃ）に示す。

（７） 速度解析

ＣＤＰギャザーに含まれている（地表座標系での）共通反射点からの反射波は、重合速度を与えてＮＭＯ補正することによって同位相となる。この性質を逆に利用して、さまざまな速度を与えてＮＭＯ補正することによって、反射波が同位相となる場合を探し出せば重合速度が求まる。これが速度解析である。

今回は１４ヶ所で速度解析を行なった。解析用の図面例を図２−２−２０に示す。また解析結果を、重合断面に重ねる様式で図２−２−２１に示す。

（８） ミュート

入射角の大きい反射パルスは、鉛直方向の見掛け速度（ＤＭＯに関連して述べた見掛け速度とは意味が異なる）が速いため分解能が悪い。そこでそのような部分を切取る操作を施す。これをミュートという。

今回は次のスキームでのミュートを行なった。

オフセット距離　　　　ミュート走時

　　　０ｍ　　　　　　　　　 ０ｍｓ

１０９ｍ ０ｍｓ

１１０ｍ ７０ｍｓ

３５０ｍ １７０ｍｓ

４５０ｍ ４１０ｍｓ

（９） 残留静補正

静補正を行なっても、ごく高周波の位相ずれは補正しきれないのが普通である。そこで統計的な手法によって、このようなずれをさらに補正することを行なう。これを残留静補正という。静補正とＤＭＯに加えて残留静補正をも施した重合断面の一部を、図２−２−１９（ｄ）に示す。

（１０）デコンボリューション（重合後）

重合前のデコンボリューションが理想的には働かないことによって残留する位相を、重合後に再びデコンボリューションを適用することによって取り除くことを、重合後デコンボリューション（ＤＡＳ）という。

今回のデータには長周期と短周期の繰り返し位相が認められたので、次の要領でＤＡＳを２回適用した。

オペレータ長　　予測長　　 設計ゲート　　　適用ゲート

　　６０ｍｓ　　　　 　８ｍｓ ３０ 〜 ４００ｍｓ 　　 全長

４００ｍｓ ２００ｍｓ　 １００ 〜 １０００ｍｓ 全長

ＤＡＳ適用前後の重合断面の一部を、図２−２−２２（ｂ）に示す。

（１１）周波数フィルタ（重合後）

ＤＡＳによって増幅された低Ｓ／Ｎ帯域の周波数成分を除去するために、周波数フィルタを再び適用した。通過帯域は３０〜１００Ｈｚとした。周波数フィルタ後の重合断面の一部を、図２−２−２２（ｃ）に示す。

（１２）速度構造モデリング

速度解析の結果得られる重合速度は、簡単な計算によって区間速度に変換できる。しかしこうして得られる区間速度は一般に精度が良くないので、地質学的なセンスをも加味しつつスムージングによって速度構造をモデリングし、マイグレーションと深度変換に利用する。

マイグレーションと深度変換に用いた速度構造モデルを、それぞれ図２−２−２３、図２−２−２４に示す。両図はそれぞれマイグレーション断面と深度断面に重ねて表示してある。

（１３）マイグレーション

ＣＤＰ重合が地表座標（ＣＭＰ）系で行なわれるイメージングであることから、重合断面は真の座標系（地下座標系）では歪んで見える。この歪みは、例えば傾斜した反射面の見掛けの傾斜角や不連続部からの回折波として現れる。

そこで地表座標系と地下座標系を対応づける方程式（波動方程式）を用いて座標変換を行ない、歪みを取り除くことを行なう。これをマイグレーションという。今回は時間領域の差分法によってこれを行なった。マイグレーション後の断面の一部を、図２−２−２２（ｄ）に示す。

（１４）深度変換

一連の処理の最後に時間軸で観測、処理された断面を、速度モデルを用いて伸縮させ深度軸で表現し直すことを行なう。これを深度変換という。