観測時に得られた記録中には、通過する車両や構造物の中の大きな機械などより発する振動等のノイズや、電磁波によるノイズ、また屈折波や表面波といった発震に伴うノイズ（コヒーレントノイズ）等さまざまなノイズが含まれている。反射法の解析処理の主な目的は、これらのノイズを多く含んだデータから必要な反射波だけを抽出することである。

現場でハードディスクに収録した各発震毎のデジタル記録は、ワークステーション

（ＳＵＮ　Ｕｌｔｒａ）に転送し、反射法探査解析システムｐｒｅｍｉｘ（ＬａｎｄＭａｒｋ社製）を用いて処理を行った。図２−５−１−２に解析処理のフローチャートを示す。

図２−５−１−２　反射法弾性波探査解析フローチャート

以下に主要な処理の概要を述べる。

（２）ジオメトリ編集

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

（３）静補正

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

（４）バンドパスフィルター（Ｂａｎｄ Ｐａｓｓ Ｆｉｌｔｅｒ）　　　　　　　　

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

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

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

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

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

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

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

＊ＣＤＰギャザー

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

図２−５−１−３ ＣＤＰギャザーの概念

（８）速度解析　　　　　　　　　　　　　　　　　　　　　　　

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

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

　　　　　　ここに、　　ｉ：ｉ番目の反射面を表す添字　　　　　　　　　　　　　　　　　　　　　　　　　Ｘ ：オフセット距離　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　Ｔ ：垂直（Ｘ＝０での）往復走時

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

（９）ＮＭＯ補正、ミュート、ＣＤＰスタック　　　　　　　　　　　

ＣＤＰスタックの目的は、ＣＤＰギャザー内の記録を加算（重合）し、ＣＤＰ位置における地下情報を表す１個の波形記録を作成することである。ＣＤＰスタックに先立ち、ＣＤＰギャザー内の各オフセット距離の波形記録をオフセットがゼロの場合の記録に変換する必要がある。この処理をＮＭＯ（Ｎｏｒｍａｌ Ｍｏｖｅｏｕｔ）補正と呼ぶ。次にＮＭＯ補正によって波形が大きく歪んだ部分や初動付近の屈折波等の不要な部分を消去する。この処理をミュートという。最後に、ＣＤＰスタックして各オフセット距離の波形記録を重合する。ＣＤＰスタックを行うことによって、速度Ｖｓｔｋを持つ反射位相だけが重ね合わされ強調され、一方、多重反射波や表面波など、このＶｓｔｋと異なる速度を持つ波の振幅は相対的に抑制される。重合後は、各ＣＤＰ地点につき１本の波形記録となり、各々ＣＤＰ点の記録として断面表示される（時間断面）。

（１０）Ｆ−Ｘデコンボリューション

スタック断面中の反射波を強調するために、コヒーレントな信号を強調し、ランダムなノイズを低減するのに効果があるＦ−Ｘデコンボリューションを使用した。Ｆ−Ｘデコンボリューション処理は、周波数領域に変換したスタック波形に予測フィルターの手法を適用し、波形の類似性を強調する処理である。

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

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

（１２）深度変換

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