Sunday, December 26, 2010

Advanced Seismic Unix Processing Part-4

Waveform Modelling

Seismic Unix memiliki fungsi sufdmod2 yang berguna untuk melakukan Waveform Modelling.

sufdmod2 memerlukan model kecepatan dalam domain kedalaman dengan format "C-style binary floating point numbers".

Berikut ini adalah script untuk melakukan Waveform Modelling untuk model kecepatan Marmousi.

Oleh karena sufdmod2 memerlukan memory yang cukup excessive, maka langkah pertama yang saya lakukan adalah mereduksi jumlah sampel (suresamp) velocity.su menjadi 25% (rf=0.25). Anda dapat menggunakan surange untuk mengevaluasi jumlah sampel sebelum dan setelah suresamp:
suresamp < velocity.su > velocity1.su rf=0.25

Selanjutnya melakukan suwind dengan j=4 untuk mereduksi jumlah trace menjadi 25%:
suwind < velocity1.su key=tracl j=4 > velocity2.su
mv velocity2.su velocity.su


Mengubah format dari su ke binary:
sustrip < velocity.su > velocity.h@ ftn=0


Konversi dari time ke depth domain:
suttoz < velocity.su > velocitydepth.su vfile=velocity.h@

Membatasi depth dengan maksimum depth=2400m:
suwind < velocitydepth.su > velocitydepth1.su tmax=2400
mv velocitydepth1.su velocitydepth.su


Mengubah format dari su ke binary untuk model kecepatan dalam domain kedalaman:
sustrip < velocitydepth.su > velocitydepth.h@ ftn=0

Ketika anda melakukan sustrip anda akan memperoleh informasi parameter berikut:
n1=348 n2=575 nt=348 ntr=575 ns=348


Gunakan text editor untuk membuat kode Waveform Modelling:

#! /bin/sh
##gunakan surange velocitydepth.su untuk memperoleh informasi d1 dan d2
n1=348
n2=575
d1=6.912000
nt=348
ntr=575
d2=9.216
ns=348
velfile=velocitydepth.h@

ximage < $velfile wbox=400 hbox=200 n1=$n1 n2=$n2 d1=$d1 d2=$d2 \
wbox=800 hbox=400 legend=1 title="Marmousi model" label1="Depth (m)" cmap=hsv2 \
label2="distance m " units="m/s" &

# input sufdmod2
xs=2000 zs=10 hsz=10 vsx=1000 verbose=2
vsfile="vsfile.su" ssfile="ssfile.su" hsfile="hsfile.su"
tmax=2.5 mt=18
label1="Depth m"
label2="Distance m"

sufdmod2 < $velfile nz=$n1 dz=$d1 nx=$n2 dx=$d2 verbose=1 \
xs=$xs zs=$zs hsz=$hsz vsx=$vsx dt=0.002 hsfile=$hsfile fmax=70 mono=0 \
vsfile=$vsfile ssfile=$ssfile verbose=$verbose \
tmax=$tmax abs=1,1,1,1 mt=$mt |
suxmovie clip=1.0 fframe=70 dframe=10 title=" %g ms" \
label1=$label1 label2=$label2 \
n1=$n1 d1=$d1 f1=$f1 n2=$n2 d2=$d2 f2=$f2 \
loop=1 sleep=20000 -geometry 800x400 &
exit 0


Gambar di bawah ini menunjukkan model Marmousi yang telah di-downsampling dengan parameter di atas:

Gambar di bawah ini menunjukkan snapshot untuk wavefield pada Model Marmousi dengan posisi source pada X=1000m dan depth=10m.


Gambar di bawah ini menunjukkan animasi wavefield pada Model Marmousi:

marmousi_fd

No comments: