Untuk mendemonstrasikan pembuatan Super CMP Gather dengan Seismic Unix, saya akan mengunakan CMP gather Marmousi yang terkontaminasi oleh noise.
Berikut adalah beberapa CMP gathers model Marmousi sebelum terkontaminasi noise:
Lakukan suaddnoise untuk mengkontaminasi model Marmuosi dengan noise:
suaddnoise < data_cmp.su > data_cmp_noise.su sn=10
suwind < data_cmp_noise.su key=cdp min=100 max=103 | suxwigb perc=90 &
Jalankan Kode Bourne Shell berikut untuk menjumlahkan 2 CMP yang berdekatan:
#!/bin/sh
##set -x
cdpmin=1
cdpmax=574
filein=data_cmp_noise.su
# Memecah data menjadi beberapa CMP gathers:
n=$cdpmin
while [ $n -le $cdpmax ]
do
suwind < $filein key=cdp min=$n max=$n > cedepe_1$n.su
n=$(( n+1 ))
done
# Menjumlahkan 2 CMP yang berdekatan dengan suop2 op=sum:
n=$cdpmin
while [ $n -le $((cdpmax-1)) ]
do
suop2 cedepe_1$n.su cedepe_1$((n+1)).su op=sum > cedepe2_$n.su
n=$(( n+1 ))
done
# Untuk menjumlahkan 4 CMP yang berdekatan, lakukan proses
# di atas untuk Super CMP gather yang baru terbentuk.
# Banyaknya CMP yang dijumlahkan tentunya harus mempertimbangkan
# karakter lapisan dan geologi. Untuk lapisan yang datar menjumlahkan 6 CMP
# masih bisa dilakukan. Sedangkan untuk lapisan miring serta perubahan geologi
# yang tiba-tiba menjumlahkan 2 atau 4 CMP mungkin sudah cukup.
# Saya sarankan untuk selalu membandingkannya dengan Single Gathernya.
rm *cedepe_1*
fileout=data_cmp_noise_supergather.su
# Membuat Bourne Shell yang bernama gabung untuk menggabungkan Super gather:
echo "#!/bin/sh" >> gabung
echo "##set -x" >> gabung
echo "cat \\" >> gabung
n=$cdpmin
while [ $n -le $((cdpmax-2)) ]
do
echo "cedepe2_$n.su \\" >> gabung
n=$(( n+1 ))
done
echo "| sushw key=tracl a=1 > $fileout" >> gabung
echo "rm *cedepe*" >> gabung
Setelah kode di atas dijalankan, run shell gabung, dengan sh gabung.
Tampilkan beberapa Single CMP gather:
suwind < data_cmp_noise.su key=cdp min=100 max=103 | suxwigb perc=90 title="Single Gather" &
Tampilkan beberapa Super CMP gather:
suwind < data_cmp_noise_supergather.su key=cdp min=100 max=103 | suxwigb perc=90 title="Super Gather (sum 2 adjacent CMPs)" &
Perbandingannya:
Berikut adalah semblance dari Single Gather:
suwind < data_cmp_noise.su key=cdp min=100 max=102 | suvelan | suximage cmap=hsv1 legend=1 title="Single Gather" &
Semblance Super Gather:
suwind < data_cmp_noise_supergather.su key=cdp min=100 max=102 | suvelan | suximage cmap=hsv1 legend=1 title="Super Gather (sum 2 adjacent CMPs)"&
Perbandingannya:
No comments:
Post a Comment