Continuous Wavelet Transform (CWT) sangat bermanfaat untuk menganalisa lapisan tipis, analis anomali 'sub-seismic', ataupun untuk mempelajari karakteristik kandungan frekuensi dari seismik.
Pada Seismic Unix, CWT dapat dilakukan dengan menggunakan fungsi sugabor. Input untuk sugabor adalah trace-trace seismik pada domain waktu, sedangkan outputnya adalah data seismik pada domain waktu-frekuensi yang dihasilkan melalui transformasi Gabor.
Pada bagian ini, saya akan menunjukkan CWT pada data seismik yang bisa diperoleh di sini (Line A).
Pada Seismic Unix, CWT dapat dilakukan dengan menggunakan fungsi sugabor. Input untuk sugabor adalah trace-trace seismik pada domain waktu, sedangkan outputnya adalah data seismik pada domain waktu-frekuensi yang dihasilkan melalui transformasi Gabor.
Pada bagian ini, saya akan menunjukkan CWT pada data seismik yang bisa diperoleh di sini (Line A).
Lakukan konversi dari format sgy menjadi su:
segyread tape=lineA.sgy conv=0 endian=0 > lineA.su
Tampilkan:
suximage < lineA.su
Lakukan surange untuk mempelajari headers dari line A:surange < lineA.su
Lakukan suwind, untuk mengetest sugabor:
suwind < lineA.su key=cdp min=100 max=100 > lineA_cdp100.su
Lakukan sugabor:
sugabor < lineA_cdp100.su > lineA_cdp100_cwt.su
suximage < lineA_cdp100_cwt.su
Run kode Bourne Shell berikut untuk melakukakn CWT pada keseluruhan line A:#!bin/sh
n3=390 #no of seismic traces
n1=1001 #no of time samples
d1=0.004 #time sampling interval
d3=1 #trace increment
n2=20 #no of frequency traces
ntr=7800 #no of seismic traces * freq traces
d2=6.25 #frequency sample interval
filein=lineA.su
sugabor < $filein | suxmovie n1=$n1 n2=$n2 ntr=$ntr n3=$n3 d1=$d1 d2=$d2 d3=$d3 title="CWT" label1="Time (sec)"
label2="Freq (hz)" loop=1 bclip=4 wclip=0 title="trace no %g"
loop=1 sleep=20000 -geometry 300x500 exit 0 &
Berikut adalah hasil CWT untuk beberapa trace dari line A:
Dari hasil di atas, saya melakukan time-variant bandpass filter dengan sutvband:
sutvband < lineA.su tf=.5,1.5,2.5,3.5 f=15,17.5,50,55 f=15,17.5,40,45 f=15,17.5,35,40 f=15,17.5,30,35 > lineA_tvf.su
Gambar di bawah ini menunjukkan penampang sebelum sutvband:
suwind < lineA.su tmax=1.5 | suximage cmap=rgb1 perc=90 title="Before TVF" &
Sesudah sutvband:
suwind < lineA_tvf.su tmax=1.5 | suximage cmap=rgb1 perc=90 title="After TVF" &
Perbandingannya:
2 comments:
Wah, datanya harus pesan dulu ya mas.
Atau anda bisa menggunakan data marine http://ensiklopediseismik.blogspot.com/2011/08/data-gathers-marine-2d.html?m=1
Post a Comment