地殻構造探査学

データ解析演習（その4）

NMO補正

この課題ではサンプルデータ demo_nmo.su を使用する。データがなければここを参考にしてコピーする。

まず、データ(CMPギャザー)を表示しよう。

$ suximage < demo_nmo.su perc=98 &

このCMPギャザーを一定速度 v=1500 m/s でNMO補正し、結果を表示してみる。

$ sunmo < demo_nmo.su vnmo=1500 | suximage perc=98 &

sunmo はNMO補正をするプログラム。
vnmo, tnmo の組合せで1次元速度構造を与えることができる。
この例では vnmo の値が1つしかないので（tmnoは省略している）、一定の速度で補正する。

図：NMO補正、（左）オリジナル、（右）速度1500 m/sでNMOした結果。

vnmo の値を 1600, 1700, 1800, 1900, 2000, 2100, 2200 と変化させて、NMO補正の結果がどのようになるか調べよう。

$ sunmo < demo_nmo.su vnmo=1600 | suximage perc=98 &
$ sunmo < demo_nmo.su vnmo=1700 | suximage perc=98 &
.....

NMO補正に使用する速度が過小 → 補正しすぎ (Over-correction) → 反射波が右上がり
NMO補正に使用する速度が適正 → 適正な補正 → 反射波が水平
NMO補正に使用する速度が過大 → 補正不足 (Under-correction) → 反射波が右下がり

速度を変えながら一定速度でNMO補正を行う定速度スキャン(CVS)を表示するシェルスクリプト create_CVS_panel.sh を動かしてみる。スクリプトの最初を見ると、

------------------------------------
# Parameters
vmin=1200
vmax=2500
dv=50
------------------------------------

ここでは速度 1200 m/s から 2500 m/s まで 50 m/s 間隔の速度でNMO補正をして並べて表示する。

$ sh create_CVS_panel.sh

図：CVSパネルのイメージ表示、（左）オリジナル、（右）速度1200 m/sから 2500 m/sまで50 m/s間隔の速度でNMOした結果。

図の横軸が重合速度を表していないのでわかりにくいが、それぞれ反射波が上から 1650, 1850, 2000, 2200 m/s あたりで適切にNMO補正されていることがわかる。これにより重合速度を求めることができる。

速度を変えながらNMO補正を行った結果を動画で表示してみる。

$ suxmovie < cvspanel_demo.su n1=1601 n2=69 n3=28 loop=1 sleep=1 perc=98 &

suxmovie は動画表示するプログラム。
時間軸方向のサンプル数は1601個。 'n1='は時間軸方向のサンプル数であるが、指定しなくてもヘッダーから読み込まれるので省略可。
一枚につき (64トレース) + (5空トレース) = 69 トレースある。速度は 1200 m/s〜2500 m/s まで27通りあって最初のCMPギャザーと合わせて画像は28枚あるので、'n2=69 n3=28' と指定する。
loop=0,1,2 は繰り返しの方式を与える。0: 繰り返さない、1: 繰り返し、2: 行き帰り。
'sleep=1' はコマ送り時間(sec)。'=1' は1秒。
その他の操作は suximage と同じ。ズームもできる。

NMO補正をした後にそれぞれのCMPの波形を重合（足し合わせる）。その効果を以下で確認する。

$ susort d2 < cvspanel_demo.su | sustack key=d2 | suxwigb &
$ susort d2 < cvspanel_demo.su | sustack key=d2 | suxwigb key=d2 &

susort は波形の順番をソート（並び替え）するプログラム。ここでは、ヘッダ d2 の値でソートする。ヘッダ d2 にはNMO補正に使用した速度値が納めてある（シェルスクリプト create_CVS_panel.sh を参照）ので、補正速度が等しいトレースが集められる。
sustack は波形の重合（スタック：足し合わせ）するプログラム。'key= ' で指定したヘッダ（ここでは d2）の値が等しいトレースを重合する。

重合の結果を表示すると、適切な速度でNMO補正されている波形の振幅が大きくなることがわかる。

図：NMO補正後に重合をした結果の波形表示。横軸は重合結果のトレース番号であるが、NMO補正で使用した速度の順番に対応している。

作成した cvspanel_demo.su は以後不要なので消去してよい。

速度解析

上記のようにNMO補正を利用して速度構造を求める方法が速度解析である。

定速度重合(CVS)したパネルを表示する。

$ suvelan fv=1200 dv=50 nv=30 < demo_nmo.su | suximage f2=1200 d2=50 &

suvelan は定速度重合した速度解析のパネルを作成するプログラム。与えられた速度で重合（振幅の足し合わせ）をする。
'fv=1200 dv=50 nv=30' で速度 1200 m/s から 50 m/s ごとに30個の速度値で計算。
suximage で表示する際に、 'f2'（横軸の開始値）に suvelan の 'fv', 'd2'（横軸の間隔）に 'dv' の値を用いる。横軸の値を速度値にするため。
suximage ではマウスの中ボタンでマウスの位置を画面に表示する。
's' キーを押すと、その時点のマウスの位置を出力（ターミナルに表示）する。ボタンをクリックする必要はない。

重合速度が正しいと足しあわされた振幅が強く（この図の表示では赤く）なる。そこで、マウスを振幅の大きい位置に移動させ、中ボタンを押すか、's' キーを押すと、時間と速度の値を読み取ることができる。

図：定速度重合(CVS)センブランスパネル。縦軸は時間 (s)、横軸は重合に用いた速度 (m/s)。

次は、's' で読んだ値をファイルに書くようにする。'mpicks= ' でデータを書き込むファイル名を与える。

$ suvelan fv=1200 dv=50 nv=30 < demo_nmo.su | suximage f2=1200 d2=50 perc=98 mpicks=pick.txt &

【注意】timeが増える方向（上から下へ）に順番に読んでいくこと（次のNMO補正に使うため）。もし、順番が混乱したら（自信がなかったら）、あとで順番を並び替えることもできる。
```
$ sort -n pick.txt > pick_sorted.txt
```
【注意】読み取りが終わったら必ずキーボードで 'q' を押して suximage を終了させること。マウスでウィンドウを閉じてはいけない。プログラムを正しく終了させないと読み取り値がファイルに書き込まれない。

読んだ時間と速度の値のファイルを次のNMO補正のプログラムの入力に使えるように形式を整える。

$ mkparfile < pick.txt > pick.par string1=tnmo string2=vnmo

pick.txt と pick.par の内容を確認しておく。

$ cat pick.txt
$ cat pick.par

'tnmo=1.0,1.5,2.0,2.5,..... vnmo=1400,1500,1600,1700,.....' で速度構造を指定する。tnmo で指定した時間における速度を vnmo で指定する。両者のデータ数は同じでなければならない。

NMO補正とCMP重合

NMO補正をする。速度構造は tnmo と vnmo の組で下のように与えるのであるが、

$ sunmo tnmo=1.0,1.5,2.0,2.5,..... vnmo=1400,1500,1600,1700,.....  < demo_nmo.su | .....

tnmo と vnmo の内容を pick.par に書き込んであるので、これを読み込む。

$ sunmo par=pick.par smute=5 < demo_nmo.su | suximage perc=98 &
$ sunmo par=pick.par smute=1.5 < demo_nmo.su | suximage perc=98 &

NMO補正をすると、オフセットの大きいところで波形が伸びてしまったり、並びがずれてしまうことがわかる (smute=5の場合)。そのような部分を切る操作をミュート(mute)といい、'smute=' で指定する。ここでは、元の波形の1.5倍以上延びた部分をミュートする ('smute=1.5')。 sunmo では 'smute=1.5' がデフォルトの設定である。

図：NMO補正の結果。反射波が水平に並んでいることがわかる。

CMP重合（スタック）してみる。

$ sunmo par=pick.par < demo_nmo.su | sustack > trace.su
$ sunmo par=pick.par < demo_nmo.su | sustack repeat=1 nrepeat=10 > trace10.su
$ suxwigb perc=98 < trace.su &
$ suxwigb perc=98 < trace10.su &

上の図のNMO補正された複数のトレースを重合する（足す）ので上の一本のトレースになる。2番目のように、'repeat=1 nrepeat=10' とすると重合結果を10回出力する。どのくらいS/N比が向上したか、元の図と比較してみる。

図：NMO補正と重合、（左）が波形、（右）がイメージ表示、それぞれの図内で左が重合結果（5本並べてある）、比較しやすいように間隔を開けて（振幅が 0 のトレースを並べてある）、右が元の波形。S/Nが向上していることがわかる。

CMP重合後のトレースは、オフセットがゼロ、すなわち、震源と受振点とが同一の点にある場合に得られる記録になる。これをゼロオフセット記録という。

課題

サンプルデータ demo_nmo.su にさらにノイズを加えてS/N比を下げたデータを作る。

$ suaddnoise sn=5 < demo_nmo.su > demo_nmo_more_noise.su

suaddnoise はノイズを加えるプログラム。出力ファイル名は何でもよい。
'sn= 'でS/N比を与える。小さいほどS/Nが低い。'sn= 'の値はいくつでもよいが、5〜10くらいが適当か。'sn=10' はたぶん楽勝。sn<5では厳しいかもしれない。

このデータを速度解析してNMO補正したのち、CMP重合しなさい。

ノイズを加えて作成したデータ（とSN比の値）と
最終的に重合して得られたトレースの波形を表示させて出力し（やり方は下記参照）、
速度解析で読み取った値 (pick.txtまたはpick.par)

とともに提出しなさい。ノイズに反射波が埋もれてしまったデータに対しても速度解析や重合はかなり有効であることがわかるだろう。

※重合結果のトレースと元のデータを同じ振幅スケールで表示させる方法

重合結果のトレース（ここでは trace.su とする）をわかりやすいように5本並べ、その横に元のデータを書きたい場合、cat を利用してデータを追記すればよい。

$ .... | sustack repeat=1 nrepeat=5 > trace5.su
$ cat trace5.su demo_nmo_more_noise.su > stacked_result_and_original.su

重合結果 (trace.su) が5本並んだ後に、元のデータ (demo_nmo_more_noise.su) が出力される。

※参考

演習の例では、見やすいように重合結果の波形（5本）と元のデータとの間に、値がゼロのトレースを5本並べてある。

$ sunull ntr=5 nt=1601 > null5.su
$ cat trace5.su null5.su demo_nmo_more_noise.su > stacked_result_and_original.su

sunull は値がゼロのトレースを出力するプログラム。
'ntr=5' はトレースの数、'nt=1601' は各トレースのサンプル数。

トレースのサンプル数はヘッダの値の範囲を出力する surange で調べることができる。

課題その3へ < 課題その4 > 課題その5へ

戻る

Last modified: Fri Sep 27 13:47:04 JST 2024