Advanced SU Part-7: Velocity Model

Seismic Unix memiliki kapabilitas untuk membuat model kecepatan. Pekerjaan ini dilakukan dengan menggunakan perintah trimodel.

Gambar di bawah ini adalah contoh model kecepatan untuk sebuah reservoir karbonat dimana windward di sebelah kanan (ditandai dengan kemiringan yang curam) serta leeward di sebelah kiri (kemiringan landai). Di atasnya terdapat lapisan klastik dan di bawahnya terdapat base karbonat serta basement dan/atau source rocks.

Berikut adalah kode Bourne Shell yang digunakan untuk menghasilkan gambar di atas:

1 #! /bin/sh
2 set -x
3
4 trimodel xmin=0 xmax=30 zmin=0 zmax=4 \
5 1 xedge=0,30 \
6 zedge=0,0 \
7 sedge=0,0 \
8 2 xedge=0,12,24,30 \
9 zedge=1,0.6,0.3,0 \
10 sedge=0,0,0,0 \
11 3 xedge=0,3,6,9,12,15,18,21,24,27,30 \
12 zedge=2.5,2.3,2.1,1.5,1.5,1.3,1.2,1.1,2.1,2.3,2.5 \
13 sedge=0,0,0,0,0,0,0,0,0,0,0 \
14 4 xedge=0,12,24,30 \
15 zedge=3,3,3,3 \
16 sedge=0,0,0,0 \
17 5 xedge=0,30 \
18 zedge=4,4 \
19 sedge=0,0 \
20 kedge=1,2,3,4,5 \
21 sfill=0,0.5,0,0,0.44,0,0 \
22 sfill=0,2.0,0,0,0.16,0,0 \
23 sfill=0,2.7,0,0,0.06,0,0 \
24 sfill=0,3.5,0,0,0.11,0,0 > modelfile.bin
25 spsplot < modelfile.bin > modelfile.eps \
26 gedge=0.5 gtri=2.0 gmin=0 gmax=1 \
27 labelz="Depth (km)" labelx="Distance (km)" \
28 dxnum=1.0 dznum=0.5 wbox=10 hbox=4
29
30 gv modelfile.eps
31 exit

Baris ke-4 memberikan batasan model dengan sumbu X antara 0 s/d 30km dan kedalaman 0 s/d 4km.
Baris 5-7: membuat permukaan bumi yang datar antara 0 s/d 30 dengan kedalaman 0 dan gradient sloth(1/kecepatan²)
Baris 8-10: membuat horizon kedua, dst. sampai baris ke 17-19 untuk horizon ke-5.
Baris 20: menunjukkan jumlah interface (batas).
Baris 21-24: mempopulasi masing-masing lapisan dengan nilai sloth.
Dimana masing-masing angka mencerminkan sfill=x,z,x0,z0,sloth,ds/dx,ds/dz
Baris 23 menunjukkan x=0 dan z=2.7 yang mempopulasi nilai karbonat dengan sloth=0.06 (4000m/s).
Baris 25-28: memplot model binary menjadi eps, baris 30 menampilkan file eps dengan perintah gv.
Jika anda belum memiliki gv, maka install dengan perintah sudo apt-get install gv

Untuk menampilkan model binary, gunakan perintah berikut:
sxplot < modelfile.bin tricolor=none cmap=hsv1 label1=depth edgecolor=red