札幌市街地の道路に沿って測線を設定したため，車両通行に伴うノイズによるＳ／Ｎの低下が懸念されたが，交通量が多いにもかかわらず，弾性波の伝播状況は良く，全体的に品質の良いデータが取得できた。図３−４−１−１、図３−４−１−２、図３−４−１−３、図３−４−１−４、図３−４−１−５、図３−４−１−６に示す現場モニター記録を例にとって，記録に現れている反射波，屈折波の特徴について，以下に述べる。

（ｱ）　豊平川の東側での発震記録〜図３−４−１−１

全体的に低周波成分が卓越しており，屈折初動の繰り返しが見られる。屈折初動の伝播距離は，他区間の記録と比較してあまり良くない。発震点近傍での屈折初動の見かけ速度は１，８００〜１，９００ｍ／ｓである。

（ｲ）　道道丘珠空港東線での発震記録〜図３−４−１−２

交通量のわりには記録の品質は良好であり，発震点から受振区間の端まで屈折初動は伝播している。豊平川の東側の区間と比較して高周波成分が含まれ，波形が鮮明となっている。西側に展開長が長く，第２層以下の屈折波（見かけ速度約３，５００ｍ／ｓ）が見られる。往復走時で１秒程度までの区間には，明瞭な反射波が見られる。

（ｳ）　丘珠空港線・篠路道での発震記録〜図３−４−１−３、図３−４−１−４

記録の品質は良好である。発震点近傍での屈折波の見かけ速度は１，６００〜１，８００ｍ／ｓ程度である。地表から往復走時で１．５秒までの区間には明瞭な反射波が見られる。特に，１〜１．５秒の区間に存在する反射波の頂点が発震点から西側にずれているので，反射面は東方向に傾斜していると考えられる

（ｴ）　主要道札幌北広島環状線での発震記録〜図３−４−１−５、図３−４−１−６

記録の品質は良好である。発震点近傍での屈折波の見かけ速度は１，６００〜１，８００ｍ／ｓ程度である。地表から往復走時で１．２秒までの区間には明瞭な反射波が見られる。反射波の頂点はほぼ発震点の直下に存在するので，反射面は概ね平らであると考えられる。

イ　データ処理方法

データの処理は，共通反射点重合法に基づき，図３−４−２に示す標準データ処理フローにしたがって実施した。

（ｱ）　データ編集

磁気テープに記録された測定データの編集を行い，データの品質を確認した。

（ｲ）　ジオメトリ作成

発震点と受振点の中点を反射点と定義し，各測線において反射点分布（図３−４−３）に基づき，重合測線（図３−４−４）を設定した。共通反射点間隔は１２．５ｍとして，重合測線に沿って共通反射点位置（ＣＤＰ位置）を設定した。

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

スイープ信号との相互相関処理後の起震車の震源波形は，ゼロ位相型となっているが，後述するデコンボリューション処理を行うと波形の歪が起こる。これを避けるために，震源波形をゼロ位相から最小位相変換するフィルターを原記録に適用した。

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

測線の屈曲の影響で，反射点は重合測線の周辺に帯状に分布するので，各共通反射点位置（ＣＤＰ）に対し，原則として５００ｍ以内に分布する反射点を集めて，ＣＤＰギャザーに編集した。

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

反射波の減衰を補正するために，ゲート長 １２００ｍｓｅｃの自動振幅調整（ＡＧＣ）を適用した。

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

発震点・受振点での応答特性の相違を補正し，多重反射波の除去など分解能の高いデータを得るための波形変換処理（デコンボリューション）を行った。

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

・オペレータ長　： ２４０ｍｓｅｃ

・予測距離　　　： ４８ｍｓｅｃ

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

（ｷ）　表層補正（静補正）

地表付近には，弾性波速度が極端に遅い（一般に，Ｐ波速度で１０００ｍ／ｓ以下）表層部（未固結表土）が存在する。この表層部の形状の不規則性は，地下からの反射波の時間のずれとなって現れ，重合後の記録において著しい品質低下を招く。この対策として地表付近の表層を除去し，その下の基底層の速度層に置換する表層補正（時間補正）を行った。表層補正値を求めるためには。タイムターム法が実用的な方法として用いられる。

タイムターム法とは，基底層からの屈折波の観測走時を，表層を波が伝わる時間の項（タイムターム）と基底層を波が伝わる時間の項とに分解し，多数の観測値からこれらの未知数を連立方程式として解く方法である。タイムターム法により求めた表層構造を図３−４−５に示す。

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

速度解析とは，ＮＭＯ補正を行うために必要な速度情報を得るための解析である。反射イベントの走時はオフセット距離に対して，近似的に（３．２）式（式３−２）のような双曲線として与えられ，その曲率は速度の関数となる。

（式３−２）

Ｖ_ｒｍｓ は，水平多層構造の場合は，地表と反射面までの地層の平均的な速度として（３．３）式（式３−３）で与えられる。反射法における速度解析とは，このＲＭＳ速度を求めることに他ならない。

（式３−３） 　　　　　　　　　　　　

実際の速度解析では，個々の反射イベントの曲率を解析してＲＭＳ速度を求めるのではなく，様々な速度に対してＮＭＯ補正と重合を繰り返し行って，最も反射イベントが鮮明で振幅が大きくなるような速度をＲＭＳ速度として抽出する。通常は，一定の速度でＮＭＯ補正を行って重合したパネル上から最適な速度を選択する定速度重合法が一般的に用いられる。

ＲＭＳ速度が求まると，（３．４）式（式３−４）から区間速度や平均速度を計算し，深度と往復走時との関係式を得ることができる。

（式３−４）

ＮＭＯ補正値（−）は，速度解析から得られたを（３．２）式に代入して式３−５

（式３−５）

で与えられる。

定速度重合法（Ｃｏｎｓｔａｎｔ Ｖｅｌｏｃｉｔｙ Ｓｔａｃｋ）により，測線に沿って１ｋｍ間隔に速度解析を実施した。速度解析例を図３−４−６−１、図３−４−６−２、図３−４−６−３、図３−４−６−４に，速度解析結果に基づく速度プロファイルを図３−４−７に示す。

また，後述するマイグレーション処理のために，反射面の傾斜の影響を除去した速度解析（ＤＭＯ速度解析）も合わせて実施した（図３−４−８）。

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

共通の反射点をもつ反射波（ＣＤＰギャザー）の受振点距離による到達時間の遅れを受振点距離０ｍの時間に補正するため，速度プロファイル（図３−４−７）に基づき，ＮＭＯ補正を行った。

（ｺ）　残差静補正

各トレースの反射イベントは，表層補正とＮＭＯ補正により，理論的には同一の走時となるが，実際には表層補正での誤差の影響などでわずかな相違が見られる。その相違を解消するための補正値をトレース間の相関から求める残差静補正を実施した。残差静補正後，再び速度解析を行い，その精度を向上させた。

（ｻ）　ＣＤＰ重合（共通反射点重合）

ＮＭＯ補正のＣＤＰギャザーを重合した。

（ｼ）　ノイズ抑制フィルタ

重合記録断面上に卓越するランダムノイズの抑制と擬似イベントの除去を目的として，Ｆ−Ｘ予測フィルタを適用し，連続性のある反射イベントを強調した。フィルタ適用後の重合記録断面図を図３−４−９に示す。

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

差分法に基づいたマイグレーション処理を行った。速度については，反射面の傾斜の影響を除去したＤＭＯ速度プロファイル（図３−４−８）の１００％値を適用した。

（ｾ）　帯域通過フィルタ

マイグレーション後の記録に対し，反射イベントを強調するために，帯域通過フィルタを適用した。浅部と深部では反射波の卓越する周波数が異なるので，往復走時に応じて，次の帯域幅のフィルタを適用し，マイグレーション時間断面図とした。

・０ 　　 〜 １．５ｓｅｃ ６／８ 〜 ３５／４０ Ｈｚ

・１．５ｓｅｃ 〜 ８．０ｓｅｃ ６／８ 〜 ２５／３０ Ｈｚ

マイグレーション時間断面図を図３−４−１０に示す。

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

マイグレーション後の記録に対し，平滑化した速度関数を用いて，時間軸から深度軸へ変換を行い，図３−４−１１及び図３−４−１２（カラー表示）の深度断面図を作成した。加えて，平成１４年度の反射断面と合わせてパネルダイヤグラム（立体視）も作成した（図３−４−１３−１、図３−４−１３−２）。なお，深度断面図の深度軸の表示は，標高０ｍを基準としている。