なお，本処理には反射法地震探査データ処理システムＰｒｏＭＡＸ を用いている．
（１）前処理
（ａ） データの転送と編集
８ｍｍ磁気テープから処理装置（ＳＵＮ Ｓｐａｒｃ２０）にデータを転送し，不必要なショット記録及びトレースを取り除く．
（ｂ） 共通反射点編集
発震点・受振点の組み合わせと，測量によって得られた各座標を用いて，各トレースのＣＭＰの位置を計算し，ショット記録を共通反射点記録に編集する．北広島測線のように屈曲した測線においては反射点が２次元的に分布するので，分布する反射点の平均的な位置に仮想の測線（重合測線）を設定し，ばらついた反射点をあたかもその測線上で反射しているように集める．
（２） 重合前フィルタ
重合前には，できるだけ反射波を強振幅かつインパルスに近い波形に変換し，反射波の周波数帯域以外のノイズを弱めたり除去したりするために，デコンボリューション・フィルタや帯域通過フィルタなどが適用される．
（ａ）デコンボリューション・フィルタ
このフィルタは震源波形，地層特性など反射地震記録に，コンボリューションの関係で含まれている基本波形をインパルスに短縮するフィルタの一種である．入力波形は以下の条件を有すると仮定する．
・反射係数時系列は無規則（ランダム）で定常である．
・基本波形（基本ウェーブレット）は最小位相型である．
一般に，反射記録は反射係数時系列に基本波形がコンボリューションされたものと考えられる．基本波形が既知であれば，これをインパルスに変換するフィルタ（これを逆フィルタと呼ぶ）は正確に設計できる．しかし，反射記録の性質上基本波形を正確に知ることは非常に困難であるため，デコンボリューション・フィルタでは，上記最小位相型の仮定を設けることにより統計的に処理されている．
デコンボリューション・フィルタの具体的な効果は以下のとおりである．
Ａ．様々な周波数成分をもつ間延びした反射波を，インパルスに近い（高周波かつ分解能の高い）波に変換する．
Ｂ．主として浅部の影響による重複反射波を除去または弱め，独立した反射波に変換する．
（ｂ）帯域通過フィルタ
信号である反射波とノイズである他の震動との周波数帯域が異なっている場合には，反射波の帯域のみを通す帯域通過フィルタをかけることにより，Ｓ／Ｎ比の向上が期待できる．そのためには周波数領域でフィルタを設計し，それをフーリエ変換して時間領域のフィルタ・オペレータを求め，反射記録にコンボリューションする．
（３）静 補 正
発震点および受振点の標高差，表層付近での弾性波の速度差による反射波の遅速，表層の厚さの変化による反射波の遅速などを補正する（前一者は標高（地形）補正，後二者は表層（風化層）補正と呼ばれる．）．これを実施しない場合には構造解釈に誤りを犯したり，水平重合法において反射波そのものを破壊することがある．
具体的には，各発震点・受振点の標高の移動平均したものを仮想の基準面とし，あたかもその基準面上で測定が行われたかのように各発震点，各受振点の記録を上下する．
なお，本処理での最終の基準面は，新篠津測線では標高５０ｍ，北広島測線では標高１００ｍ，とした．
表層（風化層）補正は，ショット記録の初動から屈折法的解析を行って求められる．
本処理では，全データを用いた屈折トモグラフィーによる表層解析を行った．屈折トモグラフィー解析により得られた表層付近の速度分布を図３−６−１４（新篠津測線），図３−６−１８（北広島測線）に示す．
（４）速度解析
ＣＭＰ重合に用いる重合速度分布を求める処理である．具体的には，様々な重合速度でＮＭＯ補正，水平重合を施し，そのうち反射波の振幅強度が最大となったときの速度が最適な重合速度として採用し，そしてこのような最適な重合速度を時間の関数として選んでいくものである．速度解析は，測線に沿ったある間隔で行われ，その間の速度は内挿により求められる．ここで求める速度は重合速度と呼ばれ，水平多層構造を仮定した場合は ＲＭＳ速度に近似的に一致する．速度解析の方法としては，定速度走査法と定速度重合法があり，本処理では両者を用いた．
速度解析を行った位置及び速度情報は時間断面図の上欄に記載している．
（５）ＮＭＯ補正及びＣＭＰ重合
ＣＭＰアンサンブル（各ＣＭＰ［発震点−受振点の中点］毎に対応するトレースを集めたもの）に対して発震点・受振点間距離の違いによる反射波の到達時間（走時）の遅れを補正し，発震点と受振点が同一位置にある場合（ゼロオフセット）の走時に合わせる操作が ＮＭＯ補正である．ＮＭＯ補正を行うにはあらかじめ地下の速度分布を設定する必要があり，通常は速度解析により得られた速度分布を用いる．このＮＭＯ補正後のＣＭＰアンサンブルを足し合わせて，反射シグナルを強調する操作が水平重合（ＣＭＰ重合）である．図３−６−１１１は，６重合の場合のＮＭＯ補正と，ＣＭＰ重合（水平重合）処理の概念図である．共通反射点Ｐでは，Ａ−Ｐ−ａ，Ｂ−Ｐ−ｂ，Ｃ−Ｐ−ｃ，Ｄ−Ｐ−ｄ，Ｅ−Ｐ−ｅ 及びＦ−Ｐ−ｆ の各々異なった経路をもつ反射記録（T） が得られる．Ｐ点からの反射波の走時は，水平距離Ｘの増加とともに遅くなっている．次に速度解析で得られた速度を用いてｐ−Ｐ−ｐ 経路の仮想の記録に合わせる（ＮＭＯ補正）と，Ｐ点からの反射波はいずれの経路のものもほぼ同じ走時の記録（U） が得られる．さらにこれらを足し合わせる（水平重合）と，ランダムなノイズや反射波以外の波は打ち消しあって相対的に弱くなり，逆にＰ点からの反射波は加算されて反射波が強調された記録（V）が得られる．
（６）残差静補正
高度補正や表層補正を施した後でも，初動屈折波と反射波の経路の違いによる時間の不規則性や，表層補正で行う２層構造仮定を採用したことによる局所的な速度の異常に起因するものは完全には補正されず，ＣＭＰアンサンブル内での同一反射の到達時間は一定ではないのが普通である．水平重合反射法地震探査においては，最適なＣＭＰアンサンブル群が得られるように統計的処理を施してこの時間差を補正し，各発震点及び受振点における２次補正値を求める．
（７）重合後フィルタ
重合後のフィルタとしては，重合断面をより地下構造を反映した断面にすることを目的に用いられ，通常，帯域通過フィルタ，Ｆ−Ｋフィルタ，デコンボリューション・フィルタなどが使用される．
本処理では，反射イベントをより際だたせるためにＦ−Ｘデコンボリューション・フィルタを適用した．
（８）マイグレーション
傾斜した地層から反射した波は各トレースの直下に表示されるため，重合記録上の反射面の傾斜や位置は地下の実際の反射面とは異なってくる．記録上での見かけの傾斜や位置を移動させて，真の傾斜および位置に復元することをマイグレーション処理という（図３−６−１２２参照）．マイグレーションの方法としては，現在波動方程式マイグレーションが最も一般的で，計算手法により差分法，Ｆ−Ｋ法などがある．
本処理では，地下の速度分布が不明確な場合に有効な差分法を採用した．マイグレーションに使用する速度としては，通常水平多層構造の場合には重合速度をそのまま用いるが，反射波が傾斜している場合に重合速度をそのまま用いるとオーバーマイグレーションの可能性が出てくる．
（９）深度変換
（８）までの処理によって得られる断面は縦軸を時間とするものである．地下の速度分布を仮定し、これを用いて縦軸を深度に変換する．深度変換に使用する速度には通常重合速度から得られる区間速度が用いられる．
（１０）断面図作成
以上の処理を施した両測線の表層速度分布，マイグレーション前の時間断面図，深度断面図，マイグレーション後の時間断面図，深度断面図及び解釈深度断面図を作成した（新篠津測線：図３−６−１３、図３−６−１４、図３−６−１５、図３−６−１６、図３−６−１７，北広島測線：図３−６−１８、図３−６−１９、図３−６−２０、図３−６−２１、図３−６−２２）．