反射法の解析処理の主な目的は，これらのノイズを多く含んだデータから必要な反射波だけを抽出することである。

現場でハードディスクに収録した各発震毎のデジタル記録は，ワークステーション（ＳＵＮ　Ｕｌｔｒａ）に転送し，反射法探査解析システムＰｒｏＭＡＸ（Ｌａｎｄ Ｍａｒｋ社製）を用いて処理を行った。

図３−２１に解析処理の流れ図を示し，次頁に主要な処理の概要を述べる。

ａ．ジオメトリ編集

各観測波形データと，それが得られた震源位置，受振点位置などを関連づけるための処理である。

ｂ．静補正

測線が起伏に富んでいたり，弾性波速度が非常に小さい表層の層厚変化があった場合，仮に地下深部の反射面が水平であっても，反射波の到達時間にばらつきが生じる。このような地表付近の不均質に起因する時間ずれを補正する処理が静補正である。静補正にはいくつかの方法があるが，ここでは，屈折波の初動走時を読み取り，この走時より各発震点・受振点におけるＤｅｌａｙ　Ｔｉｍｅを求め，このＤｅｌａｙ Ｔｉｍｅ（表層部の伝播時間）で補正を施す処理（Ｒｅｆｒａｃｔｉｏｎ Ｓｔａｔｉｃｓ）を用いた。

ｃ．バンドパスフィルター（ｂａｎｄ ｐａｓｓ ｆｉｌｔｅｒ）　　　　　　　　　

周波数フィルターの一種。観測された記録には，表面波のような反射波以外の波やバックグラウンドノイズが含まれている。これらのいわゆるノイズと反射波の周波数帯域の違いに着目して，反射波の信号と異なる周波数を持つノイズを減少させる処理である。

ｄ．ＡＧＣ（Ａｕｔｏｍｔｉｃ Ｇａｉｎ Ｃｏｎｔｒｏｌ）　　　　　　　　　　　　

観測された記録は，屈折波や表面波の振幅が大きく，反射波の振幅はこれらの波に比べて小さいのが普通である。このような振幅の小さい反射波を初動付近の波の振幅と同程度の大きさになるように強制的に増幅する処理をＡＧＣと呼ぶ。

ｅ．デコンボリューションフィルター（ｄｅｃｏｎｖｏｌｕｔｉｏｎ ｆｉｌｔｅｒ）

観測された反射波形は，地層の音響インピーダンス変化にともなう反射係数列と地下を伝わる波の基本波形のコンボリューションであると考えられる。したがって，基本波形の逆特性を持つフィルタを設計し，これに観測波形を入力すると，地下の反射係数列を得ることができる。このような処理をデコンボリューションフィルタと呼ぶ。この処理により，多重反射波が除去され（弱められ），反射波はインパルスに近い（周波数が高く独立している）波に変換される。

ｆ．ＣＤＰソーティング（ＣＤＰ Ｓｏｒｔｉｎｇ）

観測に際しては，１回の発震で９６受振点の波形記録が得られ，１発震点毎の記録として収録される。以下の処理を行うためには，すべての記録がＣＤＰギャザーごとに並んでいる方が扱い易いため，発震点毎の記録をＣＤＰギャザーごとに並び変える作業を行う。この並び変えを，ＣＤＰソーティングと呼ぶ。

＜ＣＤＰギャザー＞

図３−２２（ａ）に示したような発震点受振点配置の観測データを並び替え，図３−２２（ｂ）に示すように反射点が共通な記録，すなわち発震点と受振点の中点が同じ位置となる記録を集める。このような記録群をＣＤＰギャザーと呼ぶ。最終的には微弱な反射波を強調させる目的でこの記録群内の記録を加算する。このような手法は，ＣＤＰ重合法（ＣＤＰスタック）と呼ばれ，反射法探査の標準的な解析法として用いられている。

図３−２２　ＣＤＰギャザー模式図

ｇ．速度解析　　　　　　　　　　　　　　　　　　　　　　　

速度解析は，ＣＤＰスタックを実行する際に必要な速度を知るために，ＣＤＰギャザー内の反射走時Ｔｉ（Ｘ）が，オフセット距離Ｘ（発震点と受振点の距離），２−ｗａｙ ｔｉｍｅ Ｔ（Ｘ＝０での反射面までの往復走時），ＣＤＰギャザー内での反射位相のみかけの平均速度Ｖｓｔｋによって次式のように表されることを利用し，ＶｓｔｋとＴを決定する作業である。

Ｔｉ（Ｘ）＝｛（Ｔ，Ｉ）２ ＋（Ｘ／Ｖｓｔｋ，Ｉ）２ ｝１／２

ここに，　ｉ ：ｉ番目の反射面を表す添字

　　　　　Ｘ ：オフセット距離

　　　　　Ｔ ：垂直（Ｘ＝０での）往復走時

　　　　　Ｖｓｔｋ：ＣＤＰギャザー内の反射位相のみかけの平均速度

ｈ．ＮＭＯ補正，ミュート，ＣＤＰスタック　　　　　　　　　　　

図３−２３にＮＭＯ補正からＣＤＰスタック処理までの過程を示す。

ＣＤＰスタックの目的は，ＣＤＰギャザー内の記録を加算（重合）し，ＣＤＰ位置における地下情報を表す１個の波形記録を作成することである。

ＣＤＰスタックに先立ち，ＣＤＰギャザー内の各オフセット距離の波形記録をオフセットがゼロの場合の記録に変換する必要がある。この処理をＮＭＯ（Ｎｏｒｍａｌ Ｍｏｖｅｏｕｔ）補正と呼ぶ（図３−２３ｂ）。

次にＮＭＯ補正によって波形が大きく歪んだ部分や初動付近のの屈折波等の不要な部分を消去する。この処理をミュートという。

最後に，ＣＤＰスタックして各オフセット距離の波形記録を重合する（図３−２３ｃ）。ＣＤＰスタックを行うことによって，速度Ｖｓｔｋを持つ反射位相だけが重ね合わされ強調され，一方，多重反射波や表面波など，このＶｓｔｋと異なる速度を持つ波の振幅は相対的に抑制される。

重合後は，各ＣＤＰ地点につき１本の波形記録となり，各々ＣＤＰ点の記録として断面表示される（時間断面）。

図３−２３　ＣＤＰ重合の概念図

ｉ．残差静補正

静補正では補正できなかった発震点や受振点近傍地表条件の違いによる反射波の走時のばらつきを補正する処理。ここで用いたＭａｘｉｍｕｍ Ｐｏｗｅｒ ａｕｔｏｓｔａｔｉｃｓは ，ＮＭＯ補正後のＣＤＰギャザーにおいて，あらかじめ指定したウインドウ内でのスタック結果のパワーが最も大きくなるような補正値を各受振点，発振点別に求め，これをすべてのＣＤＰに対して行い，測線全体を通して最もスタック結果のパワーが大きくなる各受振点，発振点の補正量を自動反復計算によって求めるものである。

ｊ．マイグレーション

ＣＤＰスタックにより得られる時間断面は，反射面が傾斜していたり，凹凸があった場合には見かけの構造しか示さない。このような時間断面を真の構造に近い断面に変換する処理をマイグレーションと呼ぶ。ここでは，周波数―空間領域で真の傾斜へ変換するＦ−Ｋマイグレーションを用いた。

ｋ．深度変換

ここまで述べたような処理を行って得られる時間断面において，その縦軸は時間を表している。縦軸を深度で表す深度断面を得るためには，速度解析で求めた速度値あるいは，ＶＳＰ探査やボーリング資料から推定された速度値を用いて，時間を深度に変換する必要がある。この処理を深度変換という。