STUDI PRE-STACK DEPTH MIGRATION PADA STRUKTUR KOMPLEK
TUGAS AKHIR Disusun untuk memenuhi syarat kurikuler Program Sarjana Geofisika
Oleh
I KOMANG ANDIKA ARIS PERMANA NIM: 12403021
PROGRAM STUDI GEOFISIKA FAKULTAS TEKNIK PERTAMBANGAN DAN PERMINYAKAN INSTITUT TEKNOLOGI BANDUNG 2008
"Ya Devi Stuyate Nityam Vibhuhairvedaparagaih Same Vasatu Jihvagre Brahmarupa Saraswati" "Sarasvathi Namastubhyam, Varade Kaamaroopini Vidyaarambham Karishyaami, Siddhir Bhavatu Mey Sada" (saraswati sloka)
LEMBAR PENGESAHAN
TUGAS AKHIR Dengan Judul:
STUDI PRE-STACK DEPTH MIGRATION PADA STRUKTUR KOMPLEK
Oleh
I Komang Andika Aris Permana NIM: 12403021
Telah diperiksa dan disetujui oleh: Bandung, Juni 2008
Pembimbing,
Wahyu Trioyoso, Ph.D NIP. 131801350
i
KATA PENGANTAR
Om Swastyastu Puji syukur kehadapan Ida Sang Hyang Widi Wasa atas anugerah serta petunjuk yang diberikan, sehingga penulis mampu melaksanakan dan menyelesaikan penulisan tugas akhir dengan judul “Studi Pre-Stack Depth Migration pada Struktur Komplek ”. Penulis sangat menyadari bahwa penulisan dan penyajian tugas akhir ini jauh dari sempurna dan banyak memiliki kekurangan. Karenanya penulis sangat menerima masukan baik berupa saran ataupun kritik dari pembaca sekalian guna tercapainya tulisan yang lebih baik dan berguna bagi kita bersama. Selebihnya penulis meminta maaf apabila terdapat kata-kata yang kurang tepat pada laporan ini, akan tetapi penulis sangat berharap laporan ini akan berguna bagi penulis dan para pembaca. Pada kesempatan kali ini, penulis ingin mengucapkan terima kasih sebesar-besarnya kepada: 1.
Kedua Orang tua, Christna, Angga,Wenny, yang telah memberikan dorongan semangat dan doa.
2.
Wahyu Triyoso Ph.D, selaku dosen pembimbing
3.
Drs. Untoro M.S. dan Dr. Gunawan Ibrahim selaku dosen wali.
4.
Sonny Winardhie, Ph.D, Dr. Nanang T. Puspito, Dr. Hendra Grandis, Dr. rer. nat. Awali Priyono, Afnimar, Ph.D, Prof. Sri Widiantoro, Drs. Muhamad Ahmad, Teddy Yudistira, MSi, atas ilmu yang diberikan selama ini.
5.
Seluruh staf Tata Usaha ex-Departemen GM dan Teknik Geofisika Akhir kata penulis ingin mengucapkan selamat membaca laporan ini semoga dapat
bermanfaat. Om Santi Santi Santi Om Bandung, Juni 2008
ii
DAFTAR ISI LEMBAR PENGESAHAN .........................................................................................................
i
KATA PENGANTAR .................................................................................................................
ii
DAFTAR ISI................................................................................................................................
iii
ABSTRAK...................................................................................................................................
v
I. PENDAHULUAN .................................................................................................................
1
I.1 Latar Belakang.................................................................................................................
1
I.2 Tujuan..............................................................................................................................
2
I.3 Batasan Masalah..............................................................................................................
2
I.4 Sistematika Pembahasan .................................................................................................
2
II. TEORI DASAR......................................................................................................................
2
2.1 Migrasi Setelah Stack (Post Stack Migration)................................................................
3
2.2 Migrasi Sebelum Stack (Pre-Stack Migration)...............................................................
4
2.3 Pre-Stack Time Migration (PSTM) ................................................................................
4
2.4 Pre-Stack Depth Migration (PSDM)...............................................................................
4
2.5 Kirchhoff Migration .......................................................................................................
5
III. DATA DAN PENGOLAHAN DATA ..................................................................................
6
3.1 Data.................................................................................................................................
6
3.2 Pengolahan Data Sintetik................................................................................................
6
3.2.1 Pengolahan Tahap 1..............................................................................................
7
3.2.2 Pengolahan Tahap 2..............................................................................................
7
3.3 Pengolahan Data Riil ......................................................................................................
8
3.4 Post Time Kirchhoff Migration (PTM Kirchhoff) .........................................................
9
3.5 Post Depth Kirchhoff Migration (PDM Kirchhoff) .......................................................
10
3.6 Pre-Stack Time Kirchhoff Migration (PSTM Kirchhoff) ..............................................
10
3.7 Pre-Stack Depth Kirchhoff Migration (PSDM Kirchhoff) ............................................
11
IV. HASIL DAN ANALISIS.......................................................................................................
12
4.1 Data Sintetik ...................................................................................................................
12
4.1.1 Kecepatan HVA dan Well Adjusment..................................................................
12
4.1.2 Migrasi Domain Waktu ........................................................................................
13
4.1.3 Migrasi Domain Kedalaman.................................................................................
14
4.2 Data Riil..........................................................................................................................
15
4.2.1 Analisis Kecepatan dan Final Stack .....................................................................
15
iii
4.2.2 Migrasi Domain Waktu ........................................................................................
16
4.2.3 Migrasi Domain Kedalaman.................................................................................
16
V. KESIMPULAN DAN SARAN.............................................................................................
17
5.1 Kesimpulan .....................................................................................................................
17
5.2 Saran ..............................................................................................................................
18
DAFTAR PUSTAKA ..................................................................................................................
18
LAMPIRAN.................................................................................................................................
19
iv
STUDI PRE-STACK DEPTH MIGRATION PADA STRUKTUR KOMPLEK Oleh I Komang Andika Aris Permana NIM: 12403021 Pembimbing: Wahyu Triyoso, Ph.D
ABSTRAK Pencitraan kondisi bawah permukaan dapat dilakukan dengan melakukan pengolahan data seismik, salah satunya dengan melakukan tahapan migrasi. Pada tahapan pencitraan posisi reflektor tahap awal, data yang didapat masih berupa domain waktu. Pada kenyataannya, gambaran bawah permukaan dalam domain kedalaman sehingga perlu dilakukan pengolahan data hingga mendapatkan hasil akhir dalam domain kedalaman. Pada akhirnya, informasi kedalaman suatu kondisi geologi diperlukan guna mendapatkan ukuran volume suatu reservoir. Ketepatan dalam pencitraan kondisi geologi bawah permukaan akan menimbulkan kendala ketika daerah yang akan dicitrakan memiliki lapisan dengan variasi kecepatan lateral yang tinggi. Kendala diatas akan sangat berpengaruh ketika melakukan tahapan migrasi, terutama Pre-Stack Depth Migration pada data seismik. Pada tugas akhir ini akan dilakukan analisis pada tahapan migrasi terutama Pre-Stack Depth Migration. Analisis terhadap tahapan migrasi dilakukan pada data sintetik dan data riil yang memiliki struktur komplek dan variasi kecepatan lateral secara signifikan. Pembangunan model kecepatan dilakukan dengan menggunakan metoda HVA (Horizon Velocity Analysis) dan kurva Time-Depth dari data sumur. Ketepatan dalam mereposisi reflektor sangat dipengaruhi oleh kualitas kecepatan. Selain itu, migrasi pre stack dalam domain kedalaman lebih mereposisi reflektor ke posisi sebenarnya secara lateral dibandingkan dengan migrasi post atau pre stack lainnya. Kata Kunci : Pencitraan kondisi bawah permukaan, Pre-Stack Depth Migration, Model kecepatan
v
STUDY OF PRE-STACK DEPTH MIGRATION IN THE COMPLEX STRUCTURE By I Komang Andika Aris Permana NIM: 12403021 Supervisor: Wahyu Triyoso, Ph.D
ABSTRACT Migration is one of many methods in seismic processing for subsurface imaging. Initial seismic processing using data in time domain, in the fact geology structure is in depth domain. Processing data in depth domain important to reposition truth reflector in subsurface, and in some case geology image structure in depth domain used to occur volume of reservoir. Determine precise imaging of geology subsurface, will be difficult if the structure is complex with high lateral velocity variation. This problem will be influence in migration phase, especially in Pre Stack Depth Migration. Purpose of this paper is doing analysis in migration, especially in Pre Stack Depth Migration phase. The analysis of migration used synthetic and real data with complex structure and significant lateral velocity variation. Velocity model building, using HVA (Horizon Velocity Analysis) method and Time-Depth curve from well data. The precise of reflector reposition are depending by velocity quality. Reposition reflectors to the truth position using Pre Stack Depth Migration method, better than Post Migration or Pre Stack Time Migration. Keyword: Subsurface imaging, Pre Stack Depth Migration, Velocity model
vi
I.
dimana
PENDAHULUAN
berasal dari lapisan datar. Bila bidang
banyak metoda geofisika yang dapat
batas merupakan bidang miring, maka titik
digunakan untuk mencitrakan kondisi bawah
permukaan.
reflektor akan mengalami pergeseran dan
Metoda
akan memiliki kemiringan yang semu.
seismik merupakan salah satu metoda
Maka dari itu diperlukan proses migrasi
geofisika yang paling sering digunakan
yang
untuk memetakan kondisi geologi bawah
Keakuratan trace
migrasi dalam domain waktu ataupun
dipetakan
Terdapat 2 macam pengolahan migrasi dalam
merupakan
tepat
dan
pemilihan
domain
kedalaman
yakni
melakukan migrasi sebelum dan sesudah
struktur komplek. Penentuan kecepatan yang
seismik bisa
data migrasi dalam domain kedalaman.
tersebut sangat sering terjadi apabila yang
reflektor
didapatkan dengan melakukan pengolahan
mampu
memberikan hasil yang memuaskan. Hal struktur
titik
titik awalnya.
memberikan hasil yang kurang tepat. Hasil kurang
untuk
pada penampang seismik agar kembali ke
kondisi geologi bawah permukaan, masih
masih
bertujuan
mengembalikan/reposisi
permukaan. Hingga saat ini pemetaan
kedalaman
gelombang
seismik yang diterima di permukaan tidak
Sampai saat ini, telah terdapat
geologi
refleksi-refleksi
stack. Migrasi Kirchhoff merupakan salah
parameter
satu metoda perhitungan migrasi yang
migrasi merupakan 2 faktor penting dalam
memerlukan waku komputasi sedikit dan
menentukan kualitas hasil akhir dari
hasil akhir yang baik apabila diterapkan
pengolahan data dengan metoda seismik.
pada daerah dengan struktur komplek. Guna menghasilkan penampang seismik
I.1
Latar Belakang
yang menyerupai kondisi geologi daerah
Pengolahan data seismik secara
yang disurvey, perlu dilakukan pengaturan
tepat, akan mempengaruhi hasil akhir yang
parameter migrasi dengan tepat. Salah
didapat. Salah satu langkah penting dalam
satunya dengan mengatur nilai migration
pengolahan data seismik yakni migrasi.
aperture.
Data yang didapat dari suatu trace seismik kadang kala tidak mencerminkan kondisi geologi bawah permukaan daerah yang disurvey. Hal ini terjadi apabila bidang batas perlapisan daerah yang disurvey bukan
merupakan
bidang
horizontal,
1
I.2
Tujuan
Bab I. Pendahuluan
Adapun tujuan dari tugas akhir
Berisikan tentang latar belakang,
ini,yakni:
tujuan, batasan masalah, dan sistematika
Mengaplikasikan metoda Pre-Stack
pembahasan.
Depth Migration dan Pre-Stack Time
Migration untuk mengatasi variasi lateral
Bab II. Teori Dasar
kecepatan yang tinggi
Menjelaskan tentang teori yang
Mengaplikasikan metoda migrasi
mendasari pengerjaan tugas akhir ini.
pada struktur geologi yang komplek.
Teori dasar meliputi penjelasan mengenai
konsep migrasi, Post Stack Migration, Pre
Stack
Membandingkan hasil antara Post Migration
dengan
Pre
Stack
Stack Migration, dan migrasi Kirchhoff.
Migration
Bab III. Data dan Pengolahan Data
Mendapatkan citra yang paling Berisikan mengenai data yang
baik dari penampang geologi bawah
digunakan pada pengerjaan tugas akhir ini.
permukaan daerah yang disurvey
Pada bab ini juga dijelaskan proses I.3
Batasan Masalah
pengolahan data dari berbentuk data shot
Adapun hal-hal yang membatasi
gather, hingga menghasilkan penampang
dalam melakukan penyusunan tugas akhir
berupa hasil migrasi.
ini, diantaranya:
Bab IV. Hasil dan Analisa
Data yang digunakan adalah data
Berisikan hasil pengolahan data
Marmousi 2 yang merupakan data seismik
dan analisa dari hasil pengolahan data.
sintetik 2 dimensi, dan data riil di daerah
Bab V. Kesimpulan dan Saran
‘Lovina’ (bukan nama sebenarnya).
Berisikan
Pada tugas akhir ini menggunakan
kesimpulan
dari
beberapa metoda untuk menghasilkan
pengerjaan tugas akhir ini, serta saran
kecepatan yang baik, akan tetapi pada
untuk pengerjaan topik yang sama di masa
tugas akhir ini tidak membahas analisa
mendatang.
kecepatan secara keseluruhan. II. I.4
TEORI DASAR Tujuan dari eksplorasi seismik
Sistematika Pembahasan
yakni mendapatkan informasi tentang
Sistematika pembahasan tugas
kondisi geologi bawah permukaan dari
akhir ini terdiri atas 5 bab, yakni:
hasil
observasi
waktu
tiba,
variasi 2
amplitudo, frekuensi, fasa, dan bentuk gelombang (Telford, 1990). Data seismik yang terekam merupakan waktu tempuh gelombang dari source ke receiver yang terpantulkan ke permukaan oleh bidang reflektor. Data yang terekam, nantinya akan disorting kedalam Common Mid Point (CMP) gather sebelum dilakukan tahap processing. Migrasi
merupakan
salah
satu
proses pada processing data seismik yang
Gambar 1. Prinsip migrasi dengan penjumlahan
bertujuan menggambarkan kondisi struktur
difraksi. (a) penampang Zero-offset, (b) Hasil
geologi sebenarnya pada seismic section.
migrasi. Difraksi difokuskan di titik A. (Yilmaz, 1987)
Prinsip dasar migrasi yakni melakukan pemfokusan energi gelombang sebagai
2.1.
akibat propagasi gelombang yang melalui
Migrasi Setelah Stack (Post Stack Migration)
bidang batas dengan sifat yang berbeda-
Post stack migration merupakan
beda serta meletakkan reflektor ke posisi
proses migrasi terhadap data seismik
sebenarnya. Selain itu, migrasi bertujuan untuk
menguatkan
resolusi
setelah stack. Pengolahan data dengan Post
dengan
Stack Migration didasarkan pada asumsi
memfokuskan energi.
semua data merepresentasikan primary
Efek difraksi yang timbul sebagai akibat
hamburan
gelombang
reflection atau difraksi. Migrasi perlu
yang
berbentuk hiperbolik, dapat direkonstruksi mendekati Secara
titik
sederhana
reflektor
sebenarnya.
dapat
dirumuskan
kecepatan
dan
adanya
variabel
dipping
horizon
yang
terekam
data
dipermukaan memiliki posisi yang berbeda dengan subsurface. Post Stack Migration
(pers.1)
dapat
= jarak antara zero ofset dengan
menggambarkan
dengan
baik
penampang geologi suatu tempat yang tidak terlalu komplek. Waktu komputasi
geophone To
karena
mengakibatkan
sebagai berikut:
x
dilakukan
= two way time pada zero offset
yang tidak terlalu lama dan biaya yang
Vrms = kecepatan RMS di To
diperlukan
relatif
murah
merupakan
beberapa alasan Post Stack Migration
3
masih menjadi pilihan dalam melakukan
gelombang yang akan ditangkap oleh
pencitraan struktur geologi.
receiver. Perbedaan kecepatan lateral yang tinggi akan mengakibatkan perbedaan
2.2
Migrasi Sebelum Stack
respon
(Pre-Stack Migration)
dengan kurva hiperbola akan berbeda
Proses merupakan
migrasi proses
sebelum
non
zero
gelombang
antara satu dengan lainnya. Hal ini dapat
offset
mengakibatkan terjadi kesalahan dalam menentukan
stack pada data seismik untuk model bumi
sebenarnya.
kecepatan
konstan,
digambarkan
stack
migration. Konsep dasar migrasi sebelum dengan
yang
titik
reflektor
yang
dapat
Pada proses PSTM, pengolahan
dijelaskan dengan pendekatan metoda
data dilakukan pada domain waktu, dengan
difraksi hiperbola. Proses migrasi sebelum
menggunakan kecepatan RMS sebagai
stack dilakukan kerena pada kondisi
kecepatan migrasi. Pada pengolahan data
reflektor dengan kemiringan lebih besar
PSTM, masukan data awal pada umumnya
o
dari 15 kurang mendapatkan hasil yang
berupa common offset gather.
baik. Selain itu, dengan melakukan migrasi sebelum stack, akan mampu memperkuat
2.4
Pre-Stack Depth Migration
sinyal.
(PSDM) Proses migrasi sebelum stack lebih
sedikit
dilakukan
karena
Migrasi dalam domain kedalaman
memerlukan
menampilkan penampang geologi dalam
biaya yang besar dan waktu komputasi
domain kedalaman, selain itu algoritma
yang relatif lama apabila dibandingkan
pre stack akan menghidarkan berbagai
dengan melakukan migrasi setelah stack.
asumsi dan penyederhanaan yang akan mengakibatkan kesalahan dalam reposisi
2.3
domain
waktu.
Bidang
event
(PSTM)
reflektor tidak hanya dapat diposisikan
Migrasi proses
dalam
Pre-Stack Time Migration
dalam
merupakan metoda
salah
seismik
satu
dengan
untuk
tepat.
continuity
dan
lebih baik.
Penjalaran
Migrasi dalam domain kedalaman
gelombang seismik sangat dipengaruhi
menggunakan masukan berupa kecepatan
oleh sifat fisis medium yang dilaluinya.
interval domain kedalaman. Kecepatan
Kecepatan medium merupakan salah satu
awal yang didapat pada analisa kecepatan
faktor
merupakan kecepatan RMS. Menggunakan
yang
yang
dengan
discontinuities akan memetakan sesar/fault
mendapatkan penampang geologi bawah permukaan
baik,
mempengaruhi
respon
4
persamaan
DIX,
kecepatan
Secara umum, pengolahan data
RMS
dengan menggunakan metoda PSDM akan
dikonversikan menjadi kecepatan interval.
menghasilkan penampang geologi yang baik, akan tetapi perlu diingat bahwa (Blestein dan Gray,2001 ) (pers.2)
metoda PSDM sangat sensitif terhadap
x
= titik kedalaman
kualitas kecepatan yang dipakai. Kualitas
U
= trace seismik
model kecepatan yang akurat merupakan
F
= filter untuk merecover shape dari
kunci keberhasilan dalam metoda PSDM.
pulsa sumber. ξ
Apabila model kecepatan yang digunakan
= koordinat trace
tepat, metoda PSDM dapat digunakan
τ(x,ξ) = travel time dari source ke
untuk mengatasi segala masalah imaging
receiver yang melalui titik x β
= amplitudo migrasi pada titik x
W
= pembobotan migrasi
seismik.
2.5
Kirchhoff Migration
Persamaan 2 merupakan persamaan
Migrasi
Kirchhoff
yang akan menghasilkan amplitudo pada
memperhitungkan
titik kedalaman x.
spreading, dan penajaman wavelet dalam
Fungsi waktu tempuh TWT
directivity,
spherical
τ (x, ξ)
melakukan penjumlahan difraksi. Silva
diposisikan sebagai turunan kedua dari ray
(1992) mengeluarkan persamaan untuk
tracing. Trace seismik (U) berhubungan
migrasi Kirchhoff sebagai berikut:
dari trace yang dihitung dari ray tracing. Stacking
terhadap
pulsa,
dilakukan
sepanjang reflektor target sebagai turunan (pers. 3)
kedua dari zona Fresnel sepanjang titik refleksi yang diketahui. Teori zona Fresnel
Pout(x1,z1,t=0) = output sample p di titik
dapat diperlihatkan pada gambar 2.
x1,z1 pada saat t=0 P(x,z0,t=T)
= input sample p yang terdefinisi pada hiperbola difraksi pada pers. 1
To
= time pada sample yang termigrasi (x1,z1)
T
Gambar 2. Teori Fresnel Zone
= two way time pada jarak ke-x
5
V
= kec. RMS pada t=0 Migrasi
Kirchhoff
dapat
diaplikasikan pada data seismik yang memiliki kontras kecepatan yang tinggi, dan sudut dip yang besar. Selain itu, waktu komputasi
migrasi
efisien.
Pada
penjumlahan
Kirchhoff migrasi
difraksi
dengan
Kirchhoff,
dibatasi
migration
Jumlah Geophone
:120
Station Interval
:25 m
Source Interval
:50 m
Near Offset
:50 m
Frekuensi Dominan
:25 Hz
sangat dengan
Parameter
adanya pembatasan interval spasial yang disebut
data
daerah
Lovina
sebagai berikut:
aperture.
Kualitas migrasi yang baik ditentukan oleh nilai migration aperture yang optimal. Nilai migration aperture yang terlalu
Jumlah Geophone
:192
Station Interval
:12,5 m
Source Interval
:12,5 m
Near Offset
:30 m
3.2
Pengolahan Data Sintetik
tinggi akan menghasilkan data dengan rasio signal to noise (S/N) yang rendah disamping memerlukan waktu komputasi
Pengolahan data dilakukan dengan
yang lama tanpa meningkatkan kualitas
2
data.
tahapan.
Tahapan
pertama
yakni
melakukan processing data dengan tahap akhir Post Time Kirchhoff Migration (PTM
III.
DATA
DAN
PENGOLAHAN
Kirchhoff) untuk berikutnya dilakukan
DATA
proses Horizon Velocity Analysis (HVA).
Pengolahan data pada tugas akhir
Tahapan
ini menggunakan software ProMAX 2D
kedua
yakni
melakukan
processing data dengan menggunakan
dan MATLAB 7.0.
kecepatan yang didapat pada proses pertama dan membandingkan hasilnya
3.1
Data
dengan data yang diolah menggunakan
Tugas akhir ini menggunakan 2 data
yakni
data
Marmousi
2
data sumur. Pada tahapan ini, hasil akhir
yang
berupa penampang seismik dalam domain
merupakan data sintetik, dan data daerah Lovina
(bukan
merupakan
data
nama riil.
waktu dan kedalaman yang dihasilkan
sebenarnya)
Parameter
melalui
data
Migration
Marmousi 2 yang digunakan sebagai
proses dan
Pre-Stack Post
Stak
Kirchhoff Kirchhoff
Migration.
berikut:
6
pada stack section sebagai guide sehingga
3.2.1 Pengolahan Tahap 1 Data awal berupa rawdata, yang
pada saat melakukan picking kecepatan
kemudian dilakukan geometri dan sorting
terjadi
kedalam CDP gather. Proses berikutnya
kecepatan
dilakukan koreksi terhadap jarak dengan
(Kurniawan, 2007).
melakukan koreksi NMO terhadap data.
konsistensi
dalam menentukan
lapisan
secara
lateral
Hasil HVA berupa model kecepatan
Proses NMO juga diterapkan pada saat
RMS
melakukan
Hasil
kecepatan stacking data dan migrasi dalam
kecepatan yang dihasilkan merupakan
domain waktu dan kedalaman dengan
kecepatan
terlebih
analisa
kecepatan.
yang
RMS
nantinya
akan
kemudian
dahulu
digunakan
melakukan
sebagai
konversi
dismoothing sebelum digunakan untuk
menggunakan persamaan DIX. Sebelum
proses stacking. Kecepatan RMS dari hasil
digunakan sebagai kecepatan stacking dan
picking dikonversi menjadi kecepatan
migrasi,
interval
terlebih dahulu.
waktu
dengan
menggunakan
kecepatan
HVA
dismoothing
persamaan DIX. Kecepatan interval waktu yang
dihasilkan,
merupakan
masukan
3.2.2 Pengolahan Tahap 2
untuk melakukan proses HVA. Hasil
Alur pengolahan data tahap 2
stacking digunakan untuk melakukan PTM
secara umum dapat digambarkan sebagai
Kirchhoff dan kecepatan yang digunakan
berikut
merupakan
kecepatan
RMS
yang
menghasilkan hasil stack yang paling baik. Parameter
Kirchhoff
yang
digunakan
sebagai berikut:
First CDP
:140
Last CDP
:1160
CDP Interval
:12,5
CDP Spacing
:1
Frekuensi Maksimal :70
Migration Aperture
Gambar 3. Flowchart pengolahan data tahap 2
:1500
Hasil migrasi digunakan sebagai masukan Gambar
untuk HVA. Proses
HVA
secara
flowchart
teoritis
forward
diatas
merupakan
modelling.
Tujuan
utama dalam pengolahan data ini yakni
merupakan salah satu proses analisa
melakukan pemodelan data sintetik dengan
kecepatan dengan menggunakan horison
7
menggunakan
metoda
dan
Post
Faktor skala dari ketiga CDP,
dihasilkan
diekstrapolasi untuk keseluruhan CDP.
penampang yang mirip dengan kondisi
Matrik ini merupakan matrik adjusment
geologi awalnya.
dari kecepatan RMS HVA untuk tiap CDP.
migrasi
sehingga
Pre akan
Hasil kecepatan interval HVA pada
Matrik
adjusment
dikalikan
dengan
tahap pertama digunakan sebagai masukan
kecepatan RMS HVA akan menghasilkan
pada pengolahan data tahap 2. Pada tahap
kecepatan RMS baru. Hasil kecepatan RMS
ini dilakukan pengambilan data kecepatan
baru digunakan sebagai input velocity
interval dari model sintetik. Kecepatan
stack dan migrasi dalam domain waktu,
interval
diambil
sedangkan untuk migrasi dalam domain
merupakan kecepatan sebenarnya dari
kedalaman kecepatan RMS baru terlebih
model
dahulu dikonversikan menjadi kecepatan
kedalaman Marmousi
yang
2.
Data
kecepatan
interval model marmousi kita anggap
interval kedalaman.
sebagai data sonic well. Pengambilan data sumur dilakukan pada 3 CDP yakni CDP
3.3
Pengolahan Data Riil.
225, 715, dan 961 (lampiran, gambar 3).
Data
riil
yang
digunakan
Data sumur digunakan sebagai kontrol
merupakan data seismik 2D di daerah
terhadap data kecepatan yang didapat dari
‘Lovina’ (bukan nama sebenarnya). Bagan
hasil picking dan HVA.
pengolahan data riil sebagai berikut
Pada nomor CDP yang sama, nilai kecepatan interval kedalaman pada data sumur,
dibandingkan
dengan
nilai
kecepatan pada model kecepatan interval kedalaman kecepatan
Perbedaan
HVA. kedua
data
ini,
nilai dapat
dibandingkan dengan menggunakan kurva Time vs Depth (lampiran, gambar 4). Pada setiap kedalaman 25 meter pada tiap-tiap CDP, dapat dicari faktor skala dari perbandingan
waktu
kedua
data.
Perbandingan waktu ini akan digunakan sebagai
faktor
skala
yang
nantinya
digunakan sebagai matrik adjusment untuk mendapatkan kecepatan RMS yang baru.
Gambar 4. Flowchart pengolahan data riil
8
Data yang awal berupa rawgeom
kecepatan RMS dan kecepatan interval
dilakukan geometri untuk menyesuaikan
kedalaman yang akan digunakan untruk
dengan kondisi akuisisi data. Data yang
proses stack akhir dan migrasi.
telah digeometri disorting kedalam CDP gather untuk kemudian dilakukan tahapan
3.4
Post Time Kirchhoff Migration
preprocessing. Tahapan berikutnya yakni
(PTM Kirchhoff)
melakukan analisis kecepatan yang akan
Data awal berupa shot gather
digunakan sebagai masukan kecepatan
disorting kedalam CDP gather, kemudian
untuk melakukan brute stack. Proses brute
dilakukan pengolahan data awal dengan
stack diperlukan untuk melihat secara
melakukan filter. Hasil CDP gather yang
umum kondisi geologi bawah permukaan
telah di filter digunakan sebagai masukan
daerah survey.
analisis kecepatan. Hasil picking kecepatan
Proses berikutnya yakni melakukan
berupa kecepatan RMS digunakan sebagai
DMO guna mengoreksi posisi reflektor
kecepatan untuk melakukan stacking CDP
sebagai
yang
gather. Langkah berikutnya melakukan
reflektor miring. Hasil proses DMO
PTM kirchhoff dengan input berupa CDP
digunakan
sebagai
masukan
untuk
gather yang telah distack, dan kecepatan
melakukan
analisis
kecepatan
kedua.
yang digunakan berupa kecepatan RMS
Picking kecepatan dilakukan dari CDP 200
hasil picking. Parameter migrasi Kirchhoff
hingga CDP 3000 dengan interval 25.
yang digunakan untuk data sintetik sebagai
Model kecepatan yang dihasilkan menjadi
berikut:
masukan kecepatan untuk melakukan stack
First CDP.
:140
data. Hasil stack data digunakan untuk
Last CDP
:1160
melakukan migrasi dalam domain waktu
CDP Interval
:12.5
dengan menggunakan kecepatan yang
Frekuensi Max.
:70
sama dengan kecepatan stacking.
Migration Aperture
:1500
Max Dip to Migrate :180
akibat
Proses
adanya
HVA
bidang
dilakukan
pada
tahapan berikutnya dengan masukan data
Hasil akhir migrasi pada data riil
berupa hasil migrasi dalam domain waktu.
hanya menggunakan kecepatan dari proses
Hasil akhir dari proses HVA yakni model
HVA, parameter migrasi yang digunakan
kecepatan interval dalam domain waktu,
sebagai berikut:
lalu dilakukan proses smoothing pada
First CDP.
:200
model kecepatan tersebut. Hasil model
Last CDP
:3000
kecepatan interval dikonversikan menjadi
CDP Interval
:6.25
9
Frekuensi Max.
:80
migrasi
Migration Aperture
:3000
sebelum dikonversi, dilakukan smoothing
Max Dip to Migrate :180
terlebih
Hasil PTM pertama, digunakan
Kirchhoff yang digunakan untuk data
post
depth. Kecepatan
dahulu.
Parameter
HVA
migrasi
sebagai masukan untuk HVA sehingga
sintetik sebagai berikut:
akan menghasilkan kecepatan RMS yang
First CDP.
:140
diharapkan lebih baik jika dibandingkan
Last CDP
:1160
dengan hasil kecepatan dengan metoda
CDP Interval
:12.5
konvensional.
HVA
Frekuensi Maks.
:70
digunakan sebagai input kecepatan untuk
Migration Aperture
:1500
proses
dahulu
Max Dip to Migrate
:180
kecepatan
Max Depth to Migrate
:4000
PTM
melakukan
Hasil
kecepatan
dengan konversi
terlebih dari
interval domain waktu menjadi kecepatan Pada data riil, parameter migrasi
RMS.
yang digunakan sebagai berikut: Post Depth Kirchhoff Migration
First CDP.
:200
(PDM Kirchhoff)
Last CDP
:3000
Input PDM Kirchhoff sama dengan
CDP Interval
:6.25
input pada PTM. Masukan data awal
Frekuensi Maks.
:80
berupa CDP gather kemudian dilakukan
Migration Aperture
:3000
filter guna mendapatkan sinyal dengan
Max Dip to Migrate
:180
frekuensi yang diperlukan. Masukan untuk
Max Depth to Migrate
:3000
3.6
Pre-Stack
3.5
PDM berupa data CDP gather yang telah distack, dan kecepatan yang digunakan
Time
Kirchhoff
untuk migrasi adalah kecepatan RMS dari
Migration (PSTM Kirchhoff)
hasil picking dengan terlebih dahulu
Input data pada migrasi sebelum
dikonversikan menjadi kecepatan interval
stack berbeda dengan input data pada
domain kedalaman.
migrasi setelah stack. Data awal berupa
Pengolahan data pada PDM juga
cdp
gather
kemudian
difilter
guna
dilakukan dengan menggunakan kecepatan
mengurangi noise lalu disorting kedalam
hasil HVA. Kecepatan interval domain
common offset gather sebelum dimigrasi
waktu HVA dikonversi menjadi kecepatan
dalam domain waktu. Kecepatan yang
interval kedalaman, kemudian digunakan
digunakan untuk PSTM yakni kecepatan
sebagai
HVA yang sebelumnya dikonversi menjadi
kecepatan
untuk
melakukan
10
kecepatan
Parameter
RMS.
migrasi
Migration Aperture
3.7
Pre-Stack
:3000
Kirchhoff yang digunakan untuk data sintetik sebagai berikut:
Depth
Kirchhoff
First CDP.
:140
Migration (PSDM Kirchhoff)
Last CDP
:1160
Masukan data pada PSDM sama
CDP Interval
:12.5
dengan masukan data pada PSTM yakni
CDP Spacing
:1
sorting data dalam bentuk common offset
Frekuensi Max.
:70
gather. Parameter migrasi Kirchhoff yang
Migration Aperture
:1500
digunakan untuk data sintetik sebagai
Sebagai pembanding hasil migrasi
berikut:
menggunakan kecepatan HVA, maka pada
First CDP.
:140
proses migrasi kali ini digunakan juga
Last CDP
:1160
kecepatan yang didapat dari kombinasi
CDP Interval
:12.5
kecepatan sumur dan HVA. Kecepatan
CDP Spacing
:1
RMS baru ini didapat dari hasil perkalian
Frekuensi Maks.
:70
antara
dengan
Max Depth to Migrate
:4000
kecepatan HVA. Sebelum digunakan untuk
Migration Aperture
:1500
matriks
adjusment
Model kecepatan yang digunakan
migrasi, maka kedua jenis kecepatan ini
pada PSDM yakni kecepatan hasil HVA
dismoothing terlebih dahulu. migrasi
data
dengan
yang
dan
hasil
migrasi
interval kedalaman dengan menggunakan
kecepatan RMS baru kemudian dilakukan
persamaan DIX. Proses PSDM juga
koreksi NMO dan distack.
menggunakan kecepatan baru hasil dari
Hasil kecepatan
HVA
dikonversi
menjadi
kecepatan
riil,
kombinasi HVA dan ekstrapolasi sumur
kecepatan yang digunakan merupakan
yang terlebih dahulu dikonversi menjadi
hasil
Pada dari
pengolahan proses
menggunakan
HVA
kecepatan
data dan
tidak
kecepatan
interval
dari
hasil
kedalaman.
Kecepatan yang digunakan
ekstrapolasi sumur. Parameter migrasi data
untuk
migrasi,
riil, sebagai berikut:
dismoothing.
dalam terlebih
domain dahulu
First CDP.
:200
Last CDP
:3000
menggunakan parameter:
CDP Interval
:6.25
First CDP.
:200
CDP Spacing
:2
Last CDP
:3000
Frekuensi Max.
:80
CDP Interval
:6.25
Migrasi
pada
data
riil
11
CDP Spacing
:1
hanya digambarkan dengan model antiklin
Frekuensi Maks.
:80
pada
Max Depth to Migrate
:3000
(lampiran, gambar 2b). Pada kecepatan
Migration Aperture
:3000
interval domain waktu ataupun kedalaman,
model
model
kecepatan
kecepatan
RMS
mendekati
HVA
model
kecepatan yang sebenarnya. Perbedaan IV.
HASIL DAN ANALISIS
kecepatan pada struktur sesar terlihat
Pada bab ini akan membahas
cukup jelas pada model kecepatan interval
mengenai hasil migrasi dalam domain
(lampiran, gambar 2c dan 2d).
waktu dan kedalaman dengan metoda Pre
Hasil akhir pada proses HVA belum
Stack Migration dan Post Stack Migration
menghasilkan model kecepatan yang sama
menggunakan kecepatan yang didapatkan
dengan kecepatan model. Hal ini dapat
dari metoda HVA dan metoda Well
dilihat dari nilai error kecepatan yang
Adjusment.
dihasilkan (lampiran, gambar 5a), dimana nilai error kecepatan RMS melebihi 25%
4.1.
Data Sintetik
pada struktur komplek. Pada perhitungan
Proses pengolahan data pada data
nilai error kecepatan interval kedalaman,
sintetik, nantinya akan diterapkan pada data
riil
disesuaikan
dengan terhadap
parameter
yang
kondisi
data.
didapatkan nilai error maksimum yakni 62%. Sama halnya dengan kecepatan RMS, nilai error terjadi pada struktur komplek
Tahapan pengolahan data sintetik, secara
terutama pada kedalaman 2700 meter
garis besar dapat dilihat pada flowchart 1
(lampiran, gambar 5c). Pada daerah yang
(lampiran, gambar 1a dan 1b).
terisikan oleh fluida, error kecepatan yang terjadi cukup tinggi antara 40%-60%.
4.1.1. Kecepatan HVA dan Well
Kurang
Adjusment Penentuan
kecepatan
tepatnya
hasil
akhir
kecepatan pada metoda HVA kemungkinan
dengan
dikarenakan
menggunakan metoda HVA memiliki hasil
adanya
kesalahan
dalam
menentukan batas horison antar lapisan
yang cukup baik. Apabila dibandingkan
pada model. Metoda ini sangat dipengaruhi
dengan data kecepatan yang benar dari
oleh ketepatan dalam menentukan batas
model, kecepatan dengan metoda HVA
lapisan,
akan menghasilkan model kecepatan yang
serta
digunakan
sedikit menyerupai model sebenarnya.
kecepatan sebagai
dasar
guide
yang untuk
menghasilkan kecepatan yang baru.
Model kecepatan pada struktur komplek
12
Pada kurva T-D untuk masukan
kedalaman pada metoda sumur memiliki
pada penentuan kecepatan dengan metoda
nilai maksimum 60%. Kesalahan tertinggi
sumur, terlihat adanya perbedaan yang
terdapat
tidak terlalu jauh antara kurva T-D HVA
kedalaman 1200 (lampiran, gambar 5d)
dengan kurva T-D sumur untuk masing-
dimana error kecepatan hingga mencapai
masing CDP (lampiran, gambar 4). Hasil
60%.
ekstrapolasi dari ketiga data perbandingan
kesalahan
antara
struktur yang memiliki kandungan fluida.
kurva
T-D,
digunakan
untuk
pada
Daerah
CDP
800-900
yang
yang
besar
pada
memiliki
nilai
terdapat
pada
menghasilkan kecepatan RMS yang baru. Model kecepatan RMS baru (lampiran,
4.1.2 Migrasi Domain Waktu
gambar 2a) lebih menyerupai kecepatan model
sebenarnya
dalam
domain
waktu
dibandingkan
dilakukan dengan metoda post dan pre
dengan model kecepatan RMS dengan
stack. Hasil metoda post stack migrasi
metoda HVA.
untuk
Penentuan
jika
Migrasi
kecepatan
HVA
dan
dengan
kecepatan kombinasi HVA dan sumur
metoda sumur sangat dipengaruhi oleh
memiliki beberapa perbedaan. Perbedaan
kualitas data kecepatan yang dipakai
antara kedua penampang ini terlihat jelas
sebagai pembanding dengan sumur pada
pada struktur sesar, dimana tidak terjadi
kurva T-D. Penentuan posisi data sumur
kemenerusan
juga mempengaruhi hasil faktor skala dari
migrasi dengan kecepatan HVA pada
kurva T-D. Hasil stack menggunakan
waktu 2200-2400 ms (lampiran, gambar 7a
masukan kecepatan baru menghasilkan
dan 7b).
penampang
lebih
pada
penampang
Hasil migrasi pre stack domain
stack
waktu untuk kecepatan HVA menunjukkan
menggunakan kecepatan metoda picking
hasil yang kurang baik pada waktu 600-
dan HVA.
1600 ms. Penampang yang dihasilkan pada
dengan
baik
sinyal
jika
dibandingkan
yang
kecepatan
masukan
hasil
Besar kesalahan model kecepatan
satuan waktu tersebut memiliki amplitudo
ketika
kecepatan
yang kurang kuat, sehingga horison yang
kombinasi HVA dan sumur lebih kecil
ditampilkan kurang jelas. Hasil migrasi
dibandingkan dengan metoda HVA. Pada
pre stack dengan kecepatan kombinasi
struktur komplek yakni CDP 600-1000,
HVA dan sumur (lampiran, gambar 8a dan
besar nilai error yang terjadi berkisar dari
8b) menghasilkan penampang yang lebih
1-18 % (lampiran, gambar 5b). Kesalahan
baik dibandingkan dengan menggunakan
kecepatan
kecepatana
RMS
menggunakan
interval
dalam
domain
HVA.
Amplitudo
yang 13
dihasilkan lebih tinggi serta struktur sesar
post depth berkisar dari 1%-15% untuk
yang lebih jelas jika dibandingkan dengan
kecepatan
migrasi menggunakan kecepatan HVA.
kombinasi HVA dan sumur (lampiran,
HVA
ataupun
kecepatan
Perbedaan antara hasil migrasi post
gambar 6a). Pada model kecepatan interval
stack dan pre stack dalam domain waktu
kedalaman HVA CDP 700-800 untuk
untuk data sintetik terlihat jelas pada
kedalaman 500-800 meter, error yang
bagian struktur komplek yakni pada CDP
terjadi 25%. Jika dilihat pada kurva akan
620-860.
menimbulkan kesalahan reposisi reflektor
Hasil
pre
lebih
stack
menunjukkan kemiripan model dengan
hingga
model aslinya pada waktu 2200-2600 ms,
Sedangkan pada posisi yang sama, untuk
hasil stacking menunjukkan kelengkungan
kecepatan kombinasi HVA dan sumur yang
dari antiklin tersebut. Strukur sesar terlihat
memberikan nilai error 10%, terjadi
lebih jelas ketika dilakukan pengolahan
kesalahan reposisi reflektor sebesar 1%-
data menggunakan pre stack.
10%.
15%
Migrasi 4.1.3
pada
model
pre
migrasi.
stack
domain
Migrasi Domain Kedalaman
kedalaman menggunakan kecepatan HVA
Proses
memberikan hasil yang kurang baik jika
migrasi
dalam
domain
kedalaman diperlukan guna mengetahui
dibandingkan
posisi sebenarnya dari reflektor. Kecepatan
menggunakan kecepatan kombinasi HVA
yang
post
dan sumur. Perbedaan terlihat pada CDP
ataupun pre stack yakni kecepatan HVA
635-1075 (lampiran, gambar 10a dan 10b)
dan kecepatan kombinasi HVA dan sumur.
terutama
Hasil migrasi post dari kedua kecepatan ini
kedalaman
memiliki
adanya ketidakmenerusan reflektor ketika
digunakan
untuk
perbedaan
migrasi
yang
sedikit.
pada
dengan
daerah
1500-2000
migrasi
sesar. meter,
Pada terlihat
Perbedaan jelas terlihat pada CDP 195-275
menggunakan
kedalaman 1000-1400 meter. Pada zona
penampang dengan kecepatan kombinasi
gas
terlihat efek multipel yang dapat
HVA dan sumur, reflektor pada daerah
dikurangi ketika menggunakan kecepatan
sesar terlihat menerus dan memperlihatkan
kombinasi HVA dan sumur (lampiran,
amplitudo yang kuat. Pada CDP 855-1075
gambar 9a dan 9b).
di
Kesalahan
dalam
kedalaman
kecepatan
HVA.
1500-2000,
Pada
terjadi
memodelkan
kemenerusan reflektor pada penampang
kecepatan akan mengakibatkan terjadinya
seismik dengan menggunakan kecepatan
kesalahan dalm penentuan posisi reflektor.
kombinasi HVA dan sumur.
Besar kesalahan yang terjadi pada migrasi
14
Perbandingan antara hasil migrasi
signifikan antara layer, maka proses
dengan metoda post dan pre stack, terlihat
migrasi akan sangat sensitif terhadap
dengan adanya kemenerusan reflektor
model kecepatan yang digunakan.
yang tidak terlihat pada penampang post
Pada proses stack data dan migrasi
stack. Pada penampang hasil pre stack,
dalam domain waktu, model kecepatan
sesar dapat terlihat lebih baik dan antiklin
yang
pada daerah sesar lebih terlihat jelas.
kecepatan RMS (lampiran, gambar 11a).
Reposisi reflektor menggunakan
digunakan
Hasil
kecepatan
merupakan RMS
didapat
model dari
metoda pre-stack menghasilkan error yang
konversi model kecepatan interval domain
lebih
dengan
waktu (lampiran, gambar 11b) yang
stack.
dihasilkan dari proses HVA. CDP gather
yang
yang merupakan masukan dalam picking
kecil
dibandingkan
menggunakan Kesalahan
metoda reposisi
post reflektor
dihasilkan berkisar dari 1% - 9% untuk
kecepatan,
kecepatan
proses preconditioning guna meningkatkan
interval
kedalaman
menggunakan metoda HVA ataupun sumur
terlebih
dahulu
dilakukan
kualitas semblance.
(lampiran, gambar 6b). Hasil ini jauh lebih
Pada pemodelan kecepatan data riil
kecil dibandingkan dengan kesalahan yang
menggunakan metoda HVA, hasil yang
dihasilkan dengan menggunakan metoda
didapat cukup baik dalam menggambarkan
post stack.
struktur
bawah
permukaan.
Model
reflektor miring pada lapisan yang dalam 4.2
Data Riil
masih dapat termodelkan seperti yang
Alur pengolahan data riil, secara
terlihat pada model kecepatan interval
umum dapat digambarkan pada flowchart
kedalaman
(lampiran,
akhir
Keakuratan model kecepatan RMS yang
pengolahan data riil berupa penampang
dihasilkan, akan sangat mempengaruhi
seismik
hasil stack dari migrasi. Hasil kecepatan
gambar yang
1c).
Hasil
telah
dimigrasi
menggunakan kecepatan HVA.
(lampiran,
gambar
11c).
yang dihasilkan merupakan kecepatan dari data seimik, tingkat keakuratan dari model
4.2.1
Analisis Kecepatan dan Final
tersebut belum tentu tepat menggambarkan
Stack
kondisi kecepatan pada daerah yang
Kecepatan merupakan salah satu
disurvey.
Guna
mendapatkan
data
faktor penting yang menentukan kualitas
kecepatan yang mendekati nilai kecepatan
migrasi. Pada struktur komplek yang
kondisi asli bawah permukaan, perlu
memiliki
digunakan data sumur sebagai parameter.
variasi
kecepatan
yang
15
Pada analisis kecepatan kali ini, kecepatan
dilakukan pada tiap-tiap gathernya Pada
langsung diekstrak dari data seismik tanpa
proses pre stack, terjadi pemfokusan
menggunakan data sumur.
energi yang lebih baik pada posisi
Proses
masukan
reflektornya. Sesar pada CDP 2170-2600
kecepatan menggunakan kecepatan RMS
(lampiran, gambar 13b), tergambarkan
hasil dari proses HVA. Pada penampang
sangat jelas pada hasil stack proses pre
final stack, reflektor sudah terlihat cukup
stack. Sesar secara jelas tergambarkan
jelas serta amplitudo yang dihasilkan
hingga 1500 ms, sedangkan pada hasil post
sudah cukup kuat. Daerah sesar pada CDP
migrasi sesar hanya tegambarkan pada
2150-2600 (lampiran, gambar 12) serta
1400
reflektor berupa antiklin terlihat cukup
kemenerusan reflektor pada proses post
jelas. Kekurangan dari penampang ini
stack terlihat lebih jelas.
final
stack,
ms.
Pada
struktur
antiklin,
yakni sesar di atas struktur antiklin belum
Hasil migrasi dalam domain waktu
tergambarkan dengan baik dan masih
secara umum belum dapat memetakan
terdapat efek difraksi sebagai akibat
kondisi geologi sebenarnya pada bawah
reflektor yang tidak memiliki keselarasan.
permukaan. Perlu dilakukan pengolahan data dalam domain kedalaman guna
4.2.2
mengetahui posisi reflektor sebenarnya di
Migrasi Domain Waktu Proses
migrasi
dalam
domain
bawah
permukaan.
Berdasarkan
hal
waktu menggunakan kecepatan RMS hasil
tersebut, maka perlu dilakukannya proses
dari proses HVA. Pada penampang post
migrasi dalam domain kedalaman.
migrasi (lampiran, gambar 13a), sesar yeng terdapat pada CDP 2150-2600
4.2.3
mengalami kenaikan pada arah up-dip,
Migrasi Domain Kedalaman Migrasi
serta menghasilkan struktur yang lebih
menggunakan
kecil
kedalaman
apabila
dibandingkan
dengan
penampang hasil final stack.
kecepatan
interval
konversi
kecepatan
proses HVA menggunakan persamaan
menghasilkan
DIX. Kecepatan yang digunakan pada
penampang yang lebih baik dibandingkan
proses migrasi ini sebelum digunakan
dengan hasil post migrasi. Data masukan
dilakukan smoothing terlebih dahulu agar
dari proses pre stack yakni data offset
model kecepatan lebih menyerupai kondisi
stack
gather
yang
umum
kedalaman
interval domain waktu yang didapat dari
Pengolahan data menggunakan pre secara
hasil
domain
telah
dilakukan
reflektor bawah permukaan.
DMO,
sehingga pada proses pre stack migrasi
16
Pada proses migrasi ini, interval
V.
KESIMPULAN DAN SARAN
CDP yang digunakan yakni 6,25 m. Nilai
Pada bab kali ini akan membahas
ini merupakan setengah dari interval
tentang kesimpulan yang didapat dari
receiver. Apabila interval CDP yang
penelitian kali ini serta beberapa saran
digunakan terlalu besar akan terjadi efek
yang dapat digunakan sebagai masukan
stretching pada data seismik. Sedangkan
untuk melanjutkan penelitian ini.
apabila interval CDP yang digunakan lebih kecil dari semestinya akan mengakibatkan
5.1
Kesimpulan
terjadinya penumpukan antar CDP satu
Pada penelitian kali ini dapat
dengan lainnya. Hal ini akan terlihat jelas
disimpulkan beberapa hal, yakni: Pada struktur komplek, terbukti
pada reflektor miring (Kurniawan, 2007). Hasil post migrasi memperlihatkan reflektor
yang
terputus
pada
bahwa
daerah
migrasi
kedalaman
dalam
domain
menggambarkan
antiklin (lampiran, gambar 14a) jika
kondisi bawah permukan lebih
dibandingkan dengan hasil migrasi pre
tepat dibandingkan dengan migrasi
stack. Pada penampang hasil post stack,
dalam domain waktu
sesar pada CDP 2150-2600 terlihat cukup
Besar kesalahan dalam penentuan
jelas hingga kedalaman 1200 m. Apabila
kecepatan, akan memberikan nilai
dibandingkan dengan hasil pre stack, sesar
kesalahan yang hampir sama pada
pada CDP yang sama terlihat hingga
saat reposisi reflektor pada migrasi
kedalaman 1400 m (lampiran, gambar 14b).
Penentuan nilai migration aperture Migrasi post stack dilakukan pada
akan
data stack dengan asumsi zero offset,
mempengaruhi
tingkat
keakuratan hasil migrasi
proses stacking common mid point hanya Migrasi pre stack dalam domain
valid pada perlapisan horisontal. Energi pada struktur yang memiliki kemiringan
kedalaman
ataupun
reflektor
ketidakmenerusan,
akan
lebih ke
posisi
mereposisi sebenarnya
dijumlahkan tidak tepat pada daerah zero
secara lateral dibandingkan dengan
offset. Proses migrasi pre stack bertujuan
migrasi post atau pre stack dalam
untuk mengurangi hal tersebut dengan
domain waktu.
mengeleminasi conventional stacking dan menjumlahkan semua
energi pada titik
sebenarnya.
17
5.2
Program, Society of Exploration
Saran Pada
kekurangan
penelitian dalam
kali
ini,
Geophysics, New Orleans.
menggambarkan
Telford,
W.,
Geldart,
L.P.,
kondisi geologi bawah permukaan secara
Sheriff, R.E., 1990,
Applied
tepat dikarenakan oleh beberapa hal.
Geophysics
Edition,
Untuk mendapatkan citra yang lebih baik,
Cambridge University Press.
perlu
digunakan
data
sumur
lainnya
Second
Yilmaz, O., 1987, Seismic Data
misalnya data densitas batuan daerah yang
Processing, Society of Exploration
akan dicitrakan selain data kecepatan
Geophysicist, USA.
lapisan dari sumur. Selain itu, perlu dilakukan penentuan migration aperture yang lebih tepat salah satunya dengan menggunakan metoda yang berbeda seperti metoda Horizontal Displacement.
DAFTAR PUSTAKA
Kurniawan,
D.,
2007,
Studi
Kualitas Model Kecepatan untuk Melakukan
Pre-Stack
Depth
Migration,
Program
Studi
Geofisika,
Institut
Teknologi
Bandung.
Priyono, A., 2001, Buku Ajar Seismik Eksplorasi untuk Bidang Ilmu Kebumian, Institut Teknologi Bandung.
Robein, E., 2003, Velocities, Timeimaging and Depth-imaging in Reflection Seismics Principle and Methods, EAGE, Netherland.
Silva, R., 1992, Antialiasing and Aplication of Weighting Factor in Kirchhoff
Migration,
Technical
18