Data Gathers Marine 2D: Radon Demultiple

Sebagaimana yang telah saya jelaskan di sini, data ini memiliki multiple yang cukup dominan.

Salah satunya adalah water bottom multiple yang terlihat pada semblance yang rendah (~1500m/s), perulangan event dan frekuensi yang tinggi pada data gathers maupun brute stack. Salah satu cara untuk mengeliminasi multiple adalah Radon Demultiple.

Berikut adalah cara melakukan Radon Demultiple pada data marine_cmp_nmo.su (Ingat Radon Demultiple diterapkan setelah koreksi NMO)

1. Dengan menggunakan text editor (vi, gedit, kedit, dll.), copy kode Bourne Shell berikut lalu berinama radon. Jalankan dengan perintah sh radon

#!/bin/sh
##set -x
cdpmin=171
cdpmax=2142
filein=marine_cmp_nmo.su
n=$cdpmin
while [ $n -le $cdpmax ]
do
suwind <$filein key=cdp min=$n max=$n | suradon choose=3 igopt=4 pmin=-600 pmax=150 nmax=400 offref=-6000 > marine_radon_cmp$n.su
n=$((n+1))
done

Kode tersebut memecah data marine ke masing-masing cdp lalu menjalankan perintah radon. Saya menggunakan choose=3 untuk menghitung forward dan reverse tranform sekaligus dan igopt=4 (abs linear), pmin dan pmax mengontrol kemiringan event yang akan dieliminasi. pmin mengontrol event yang miring ke atas dan pmax event yang miring ke bawah.

Karena multiple adalah event yang miring ke bawah maka saya melakukan pembatasan yang cukup agresif yakni 150. nmax=400 adalah jumlah trace maksimum per cdp dan offref=-6000 adalah offset reference yang jauh lebih besar dari offset data tersebut. Untuk mempelajari lebih jauh anda bisa melakukan test untuk merubah parameter tersebut dan mengetikkan suradon pada terminal linux.

Gambar di bawah ini adalah data input dan hasil transformasi radon dengan option di atas. Perhatikan bahwa setiap sample data dengan amplitudo nol akan menghasilkan artifact pada Radon. Untuk menghilangkan artifact tersebut saya harus melakukan mute.

Data Input

Hasil Radon Demultiple

2. Ekstrak mute dari NMO data yang memiliki stretch mute.

Kode berikut memecah data menjadi masing-masing cdp.

#!/bin/sh
##set -x
cdpmin=171
cdpmax=2142
filein=marine_cmp_nmo.su
n=$cdpmin
while [ $n -le $cdpmax ]
do
suwind <$filein key=cdp min=$n max=$n > marine_nmo_cmp$n.su
n=$((n+1))
done


3. Ekstrak mute

Perintah berikut menghasilkan dummy gathers dimana sample data dengan nilai amplitudo nol akan nol dan yang lainnya akan bernilai 1.

#!/bin/sh
##set -x
cdpmin=171
cdpmax=2142
n=$cdpmin
while [ $n -le $cdpmax ]
do
suop2 marine_nmo_cmp$n.su marine_nmo_cmp$n.su op=quo > marinecmp_nmo_mute_cmp$n.su
n=$((n+1))
done

Contoh dummy gather dari perintah di atas.

4. Menerapkan mute:
Intinya meng-kalikan hasil Radon dengan dummy gather.

#!/bin/sh
##set -x
cdpmin=171
cdpmax=2142
filein=angolacmp_nmo.su
n=$cdpmin
while [ $n -le $cdpmax ]
do
suop2 marine_radon_cmp$n.su marinecmp_nmo_mute_cmp$n.su op=prod > marine_radon_muted_cmp$n.su
n=$((n+1))
done

Hasil Radon yang telah dimute.

Perbandingannya:

5. Menggabungkan semua cdp radon

Copy kode berikut lalu beri nama gabung dan jalankan dengan sh gabung

#!/bin/sh

cdpmin=171
cdpmax=2142

echo "#!/bin/sh" >> gabungradon
echo "##set -x" >> gabungradon
echo "cat \\" >> gabungradon
n=$cdpmin
while [ $n -le $((cdpmax)) ]
do
echo "marine_radon_muted_cmp$n.su \\" >>gabungradon
n=$(( n+1 ))
done
echo " > marine_radon_all_muted.su" >>gabungradon


6. Jalankan sh gabungradon dan stack.
sustack < marine_radon_all_muted.su > stack_radon.su
suximage < stack_radon.su perc=90 d1=0.004 f1=0