Advanced Seismic Unix Processing Part-1

Radon Demultiple

Pada bagian ini, saya akan menunjukkan bagaimana caranya melakukan eliminasi multiple melalui Transformasi Radon.

Eliminasi multiple dengan teknik Radon atau dikenal juga dengan Radon Demultiple memiliki 4 tahapan utama:

1. Koreksi NMO (CMP gather)
2. Transformasi ke domain τ-slope index
   
3. Desain Mute
4. Inverse τ-slope index
   

Setelah koreksi NMO, karakteristik multiple akan berbeda dengan reflektor. Jika nilai kecepatan NMO yang dipilih tepat maka reflektor akan menjadi datar (flat) sedangkan multiple masih melengkung ke arah bawah (karena nilai kecepatannya lebih rendah). Jika ditransformasi ke dalam domain τ-slope index (baca: sumbu y adalah waktu dan x adalah kemiringan event), maka reflektor yang datar akan memiliki slope index bernilai sekitar nol dan multiple akan bernilai positif. Sedangkan untuk event yang melengkung ke atas (over NMO-ed, difraksi, side swipe) akan memiliki nilai slope index negatif.

Gambar di atas (kiri) menunjukkan CMP sintetik dengan reflektor pada waktu 0.2s dan reflektor yang bersamaan dengan multiple pada 0.4s, sedangkan gambar sebelah kanan adalah CMP dalam domain τ-slope index. Perhatikan reflektor yang flat pada waktu 0.2s dan 0.4s memiliki nilai slope index di sekitar nol, sedangkan multiple pada 0.4s memiliki nilai slope index sekitar 100.

Berikut adalah tahapan untuk menghasilkan gambar di atas:

1. Run kode Bourne Shell berikut:
   

#!/bin/sh
npl=3 # number of planes
nt=200 # number of time samples
ntr=100 # number of traces
taper=0 # no end-of-plane taper = 1 taper planes to zero at the end
offset=400 # offset
dt=0.004 # time sample interval in seconds
suplane npl=$npl nt=$nt ntr=$ntr taper=$taper offset=$offset \
dip1=0 len1=128 ct1=100 cx1=50 \
dip2=1 len2=128 ct2=115 cx2=50 \
dip3=0 len3=128 ct3=50 cx3=50 > planemul.su

2. Tampilkan:

sushw < planemul.su key=offset a=0 b=10 | sufilter f=0,5,50,60 > planemul_fil.su
suximage < planemul_fil.su title="Sebelum Radon Demultiple" key=offset label1="time (s)" label2="offset (meters)" wbox=500 hbox=500 &

3. Lakukan Tranformasi Radon dan tampilkan:

suradon < planemul_fil.su igopt=3 interp=4 choose=0 depthref=1000 interoff=0 offref=1190 pmin=-300 pmax=300 > radon.su
suximage < radon.su perc=99 label1="tau (s)" label2="slope index" title="suplane data Radon transformed" wbox=500 hbox=500 &

Gambar di atas adalah CMP dalam domain τ-slope index, sebelum (kiri) dan setelah mute.

Untuk melakukan Mute:

sumute < radon.su key=f2 \
xmute=-300,30,31,300 \
tmute=0.8,0.79,0.01,0 \
mode=1 ntaper=10 > radon_muted.su

suximage < radon_muted.su label1="tau (s)" label2="slope index" title="suplane data Radon transformed-muted" wbox=500 hbox=500 perc=99 &

Berikut adalah perintah untuk melakukan inversi τ-slope index setelah dimute serta menampilkan hasilnya:

suradon < radon_muted.su igopt=3 interp=4 choose=4 depthref=1000 interoff=0 offref=1190 pmin=-300 pmax=300 | suximage title="Setelah Radon Demultiple" key=offset label1="time (s)" label2="offset (meters)" wbox=500 hbox=500 &

Dari gambar di atas terlihat bahwa, multiple pada 0.4s telah dieliminasi dengan sempurna.

Gambar di bawah ini menunjukkan muting untuk data real pada domain τ-slope index . Kiri (sebelum mute) kanan (setelah mute).

Pemilihan batas pengambilan mute, harus dilakukan dengan hati-hati, jangan sampai melukai sinyal reflektor, mengingat berbagai hal seperti error pemilihan nilai kecepatan NMO, anisotropic, dll. Jadi, sebaiknya kita melakukan test dengan mengembalikannya ke domain waktu-offset, lalu lihatlah perilaku reflektornya. Gambar di bawah ini menunjukkan perbandingan CMP gather sebelum (kiri) dan setelah Radon Demultiple(kanan).