Monday, January 9, 2012

Advanced SU: Median Filter

Seismic Unix menyediakan fungsi sumedian yang berguna untuk mengeliminasi event seismik dengan moveout (baca: kemiringan) tertentu.

Fungsi ini sangat berguna baik dalam pengolahan data seismic refleksi maupun Vertical Seismic Profiling (VSP).

Untuk mendemonstrasikan kegunaan sumedian pada data VSP silakan download data VSP di sini.

Link di atas direlease oleh Pemerintah Negara Bagian Victoria, dimana data tersebut hanya boleh digunakan untuk keperluan Studi dan Penelitian.

Selain data VSP terdapat juga link yang berisi data seismik (link1-link2), well-log, dll.

Untuk data VSP, saya akan menggunakan data VSI_007_A_gac_wavefield_z.sgy

Convert sgy ke su:

segyread tape=VSI_007_A_gac_wavefield_z.sgy verbose=1 endian=0 > VSI_007_A_gac_wavefield_z.su

Gunakan shell berikut untuk melakukan windowing, menampilkan dan median filter.

#!/bin/sh
suwind < VSI_007_A_gac_wavefield_z.su key=tracl
| suwind key=tracl min=20 max=120 tmin=0.1 tmax=0.6 | suximage perc=90 d1=0.002 title='ORIGINAL'&
suwind < VSI_007_A_gac_wavefield_z.su key=tracl
| sumedian tshift=1.25,0.6 xshift=1.11,149.436 median=1 | sumedian tshift=0.396025,0.972893 xshift=3.0093, 47.5488 median=1 | suwind key=tracl min=20 max=120 tmin=0.1 tmax=0.6 | suximage perc=90 d1=0.002 title='DOWNGOING WAVES'&

suwind < VSI_007_A_gac_wavefield_z.su key=tracl
| sumedian tshift=0.11,0.76 xshift=1.65,215.53 median=1 | sumedian tshift=0.396025,0.972893 xshift=3.0093, 47.5488 median=1 | suwind key=tracl min=20 max=120 tmin=0.1 tmax=0.6 | suximage perc=90 d1=0.002 title='UPGOING WAVES'&

Parameter penting sumedian adalah tshift (nilai sumbu Y) dan xshift (nilai sumbu X). tshift-xshift merupakan pasangan beberapa point untuk event-event dengan kemiringan tertentu. tshift dapat berupa waktu, depth, dll. dan xshift dapat berupa offset, tracl, dll.

Di bawah ini adalah rekaman VSP yang masih memiliki upgoing waves dan downgoing waves.
Downgoing waves:

Upgoing waves:
Perbandingannya:
Photobucket

No comments: