Saturday, November 20, 2010

Seismic Processing with Seismic Unix - Part1

Sebagaimana yang kita ketahui, Seismic Unix memiliki ratusan fungsi yang berguna untuk mengolah data seismic.

Untuk melakukan latihan pengolahan data seismik, berikut adalah tahapan yang bisa diikuti:

Tahap1:
Download data seismik 2 dimensi dengan format segy (388.9MB) di sini

Tahap2:
Simpanlah file tersebut dalam folder tertentu, katakanlah processing, lalu ekstrak sehingga kita akan memperoleh beberapa file berikut: Line_001.sgy, Line_001.SPS, Line001.RPS, Line_001.XPS, dan Line_001.TXT. File dengan ektensi SPS dan RPS berisikan informasi sumber-penerima seperti indeks (nomor), koordinat, elevasi, statik, dll. File XPS berisikan informasi hubungan sumber-penerima (relational) dan TXT berisikan informasi tentang parameter survey. Berikut adalah isi dari Line_001.TXT:


Tahap 3:

Untuk melakukan processing dengan Seismic Unix, terlebih dahulu kita melakukan konversi format data seismik dari segy ke su.

segyread tape=Line_001.sgy verbose=1 endian=0 > Line_001.su


Pada perintah di atas, saya mendefiniskan endian=0 karena saya menggunakan mesin little endian (Laptop/PC). Jika anda menggunakan workstation anda harus medefinisikan endian=1.

Tahap 4:
Analisa header data su dengan perintah

surange < Line_001.su

Sehingga diperoleh:

Dari informasi di atas terlihat bahwa data seismik ini memiliki jumlah trace 71284, ep adalah shot number dari 32 sampai 282 = 251 shots, serta indeks trace dalam setiap shot tracf dari -1 sampai 282. Dengan kata lain jumlah trace dalam setiap shot adalah 284.

Tahap 5:
Untuk mengevaluasi data yang kita miliki lakukan peritah berikut:

suwind < Line_001.su key=ep min=32 max=32 | suxwigb perc=80

Perintah di atas, kita memilih data dengan suwind untuk shot ke 32, lalu ditampilkan sebagai wiggle dengan suxwigb dengan persentase amplitudo 80.

Gambar di bawah ini (kiri) adalah hasil dari perintah di atas, serta zoom in (kanan) dari trace-trace awal. Terlihat bahwa 2 trace pertama adalah data source signature yang didefinisikan dengan tracf=-1 dan tracf=0.


Lakukanlah proses ini untuk beberapa tempat dengan ep berbeda. Kesimpulan yang saya peroleh adalah semua shot memiliki 284 trace dengan 2 trace pertama sebagai source signature yang ditanamkan pada setiap shot record.

Tahapan berikutnya

6 comments:

eko_pcs18 said...

halo mas agus.... saya mau tanya ni....
setelah saya input di terminal command seperti yang anda tulis (suwind < Line_001.su key=ep min=32 max=32 suxwigb perc=80 ), muncul (suwind: fputtr.c: segy output can't be tty)... pada terminal dan gambar evaluasi data seperti mas agus katakan tidak muncul.... mohon pencerahan... terimakasih

Eko, jogja

Agus Abdullah, PhD said...

Tanda pipe '|' nya terlewat. Coba test lagi dengan kode ini:

suwind < Line_001.su key=ep min=32 max=32 | suxwigb perc=80

Nanda said...

mas agus, script di part 1 ini tu untuk shotgathers ya?

Agus Abdullah, PhD said...

Ya, data input nya shot gathers.

Paulala said...

Halo mas agus, saya lagi coba SU, dan ada problem "surange: fgettr.c: bad first header"

kira-kira bagaimana ya mas? terimakasih.

Agus Abdullah, PhD said...

Mungkin headernya/sebagian header masih null.
Perlu dilakukan sushw. Lihat contoh penerapannya di:
http://ensiklopediseismik.blogspot.co.id/2010/11/seismic-processing-with-seismic-unix_25.html?m=1