## Effect of Conductive Minerals in the Determination of Water Saturation in “Y” Field

## (Pengaruh Mineral Konduktif Dalam Penentuan Saturasi Air Di Lapangan “Y”)

## Panduwinata Arifin, Ratnayu Sitaresmi Hendri, Benyamin
Prodi Magister Teknik Perminyakan, Fakultas Teknologi Kebumian & Energi, Universitas Trisakti, Jakarta

## Sari
Evaluasi  formasi  adalah  salah  satu  cabang  ilmu  dari  teknik  perminyakan  untuk  mempelajari  karakteristik  reservoir  serta permasalahan  yang  berhubungan  dengan  keberhasilan  dalam  penemuan  cadangan  hidrokarbon.  Cadangan  hidrokarbon dipengaruhi oleh karakteristik batuan reservoir yang terdiri dari porositas (  ), saturasi air (S w ) dan permeabilitas (k)Pada wilayah  resistivitas  rendah  sering  ditemukan  ketidaknormalan  dalam  evaluasi  formasi.  Ketidaknormalan  itu  berupa  nilai resistivitas  yang  dihasilkan  rendah  sehingga  pada  saat  perhitungan  saturasi  air,  didapat  nilai  yang  tidak  akurat.  Pada penelitian  ini,  dilakukan  proses  evaluasi  terhadap  fenomena  ketidaknormalan  dalam  evaluasi  daerah  resistivitas  rendah. Evaluasi dilakukan dengan menganalisa data log dan data core. Evaluasi pada daerah resistivitas rendah bertujuan untuk memperoleh  faktor  koreksi  dari  ketidaknormalan  yang  terjadi.  Ketidaknormalan  dalam  perhitungan  nilai  resistivitas diakibatkan oleh adanya mineral konduktif yang mengakibatkan penurunan dari nilai resistivitas formasi. Faktor koreksi ini dapat digunakan dalam perhitungan Archie yang menghasilkan harga saturasi air, sehingga menghasikan nilai yang lebih akurat dalam perhitungan harga saturasi.
Kata-kata kunci: porositas, resistivitas, saturasi, permeabilitas, core, log .

## Abstract
Formation evaluation is a part of petroleum engineering to study reservoir characteristics as well as issues related to success in the discovery of hydrocarbon reserves. Hydrocarbon reserves are influenced by reservoir rock characteristics consisting of porosity (  ), water saturation (S w ), and permeability (k.) In low resistivity regions are often found abnormalities in formation evaluation. The abnormality is a low resistivity value produced so that an inaccurate value obtained when water saturation is calculated. In this study, an evaluation process of the abnormality phenomena was performed on low resistivity areas. The evaluation is done by analyzing both log and core data. Evaluation in low resistivity areas aims to obtain a correction factor of the  abnormality  that  occurs.  Abnormalities  in  the  calculation  of  resistivity  values  caused  by  the  existence  of  conductive minerals which resulting in a decrease of the resistivity value of the formation. This correction factor can be used in Archie to calculate water saturations, therefore more accurate values of saturation were found.
Keywords: porosity, resistivity, saturation, permeability, core, log
*Ratnayu Sitaresmi ( corresponding author ): E-mail: rsitaresmi@trisakti.ac.id

## I. INTRODUCTION
Evaluasi formasi adalah salah satu cabang ilmu dari  teknik  perminyakan  untuk  mempelajari karakteristik  reservoir  serta  permasalahan  yang berhubungan dengan keberhasilan dalam penemuan cadangan  hidrokarbon.  Cadangan  hidrokarbon dipengaruhi  oleh  karakteristik  batuan  reservoir yang terdiri dari porositas (  ), saturasi air (S w ) dan permeabilitas (k). Untuk mengetahui permasalahan diatas  diperlukan  beberapa  jenis  pekerjaan,  antara lain: pengambilan contoh batuan (coring), logging, well testing, dan lain-lain.
Pada  penelitian  ini,  dilakukan  proses  evaluasi terhadapat  kondisi  ketidak  normalan  dalam evaluasi  daerah  resistivitas  rendah.  Evaluasi  pada daerah resistivitas rendah bertujuan untuk
memperoleh  faktor  koreksi  dari  ketidaknormalan yang terjadi, sehingga menghasikannilai yang lebih akurat dalam perhitungan harga saturasi.

## II. TEORI DASAR
Evaluasi  logging  akan  menghasilkan  data-data yang  menunjukan  sifat-sifat  dari  batuan  reservoir. Sifat-sifat  yang  dihasilkan  dari  evaluasi  logging berupa  nilai  porositas  (  ),  resistivitas  batuan,  dan saturasi air (S w ).

## Porositas
Porositas  merupakan  total  volume  dari  ruang pori  batuan  berbanding  dengan  volume  total batuan, poroitas pada formasi yang tidak solid akan bergantung  terhadap  distribusi  ukuran  butiran,
tidak  padaukuran  butiran  mutlak.  Porositas  akan memiliki  nilai  yang  tinggi  pada  saat  butiran berukuran  hampir  sama,  porositas  akan  bernilai rendah pada kondisi butiran batuan yang bervariasi. Pada  lubang  pori  butiran  yang  kecil  akan  mengisi ruang diantara butiran yang besar. (Dewan. 1983).

## Permeabilitas
Permeabilitas  adalah  suatu  sifat  dari  media berpori  yang  mengukur  kapsitas  dan  kemampuan dari  batuan  untuk  mengalirkan  fluida.  Darcy mendefinisikan  permeabilitas  untuk  mengetahui kemampuan batuan mengalirkan fluida, yaitu "Laju aliran  fluida  homogen  melalui  media  berpori berbanding lurus dengan gradient tekanan dan arah normal  terhadap  luas  penampang  lintang  dari  arah aliran serta berbanding terbalik dengan viskositas.

## Saturasi
Menurut  Harsono  (1997)  saturasi  air  atau kejenuhan  air  merupakan  bagian  dari  ruang  pori yang  terisi  oleh  air.  Apabila  terisi  hidrokarbon maka  disebut  saturasi  hidrokarbon.  Hidrokarbon mengisi  ruang  pori  batuan  berawal  dari  batuan yang  mula-mula  terisi  air  dan  selang  masa perubahan geologi, minyak dan gas yang terbentuk di  tempat  lain  berpindah  ke  formasi  batuan  yang berpori,  menggantikan  air  pada  ruang  yang  lebih besar.
Saturasi dapat dihitung menggunakan persamaan (NN, 1996: 8).
√
Archie Equation :
(  )
atau Indonesian Equation :


## Mineral Konduktif
Mineral  konduktif  merupakan  mineral  yang dapat  menghantarkan  arus  listrik  dengan  baik, sehingga dengan adanya kehadiran mineral tersebut pada  formasi  akan  menghasilkan  nilai  resistivitas yang rendah.
Ketika clean reservoir dipengaruhi oleh mineral konduktif  dan  clay,  resistivitas  akan  mengalami penurunan, hal ini diakibatkan karena  mineral clay meningkatkan konduktivitas batuan dengan
meningkatkan  konduktivitas  air  yanng  ada  pada pori batuan.
Clavier,  Heim  dan  Scala  menemukan  bahwa dengan  mengasumsikan  interpretasi  dari  efek mineral  konduktif  merupakan  sebuah  model rangkaian  paralalel.  Hubungan  antara  pengaruh kunduktivitas  mineral  dengan  volume  mineral didapat  hubungan  yang  linier.  Dengan  adanya hubungan  yang  linier  ini  maka  pengaruh  mineral konduktif  dapat  dievaluasi  dengan  menggunakan nilai konduktivitasnya.(Clavier, 1990).
Model resistivitas dari zona air dan hidrokarbon tanpa pengaruh mineral konduktif

(  )

## Konversi kedalam konduktifitas
Model konduktifitas dari mineral konduktif dan hubungan  resistivitas  mineral  konduktif  dengan konduktifitas  air  dan  konduktifitas  hydrokarbon. Air + konduktif dan  hidrokarbon  +  konduktif, menggunakan teori sirkuit paralalel (Clavier,1990):
Konversi  konduktifitas ke dalam nilai resistivitas:
Nilai resistivitas yang dipengaruhi oleh mineral konduktif  dibandingkan  dengan  nilai  resistivitas awal,  dari  hasilkan  perbandingan  tersebut  kita  kan mendapat pengaruh dari mineral konduktif terhadap nilai dari resistivitas pada formasi.
Dengan menghubungkan antara nilai konduktivitas dari mineral konduktif tersebut, pada daerah  clean  sand  kita  dapat  menggunakan  rumus archie  untuk  menghitung  nilai  saturasi  air  (S w ). Dengan  menghilangkan  pengaruh  konduktivitas dari mineral konduktifnya.
(  )

## Identifikasi Mineral Konduktif
Mineral  konduktif  pada  batuan  dapat  diidentifikasi dengan  menggunakan  log  resistivitas,  log  densitas, dan  log  sonic  dan  data  core.  Kehadiran  mineral konduktif pada batuan reservoir akan mengakibatkan  pembacaan  resistivitas  menjadi rendah,  karena  pengaruh  dari  sifatnya  yang konduktif.

## III. PEMBAHASAN Pengumpulan data
Untuk mendukung hasil penelitian agar pembahasan  dapat  dilakukan  dengan  lebih  efektif, maka dikumpulkan data-data sebagai penunjang.Hal  pertama  yang  dilakukan  adalah analisa kualitatif dengan tujaun menemukan lapisan batuan  yang  diperkirakan  merupana  lapisan prospek minyak bumi dan gas.
Dalam analisa ini digunaka data dari pembacaan  log  dan  data  dari  batuan  core,  agar didapat data yang lebih lengkap.

## Data dari kurva log
Data  –  data  yang  diperoleh  dari  kurva  log merupakan  data-data  terpenting  dalam  penelitian ini  karena  data-data  tersebut  menekankan  seluruh pengukuran dalam lubang yang meliputi parameter-parameter  sifat  fisik  batuan  yang diidentifikasi  sebagi  reservoir.  Data  yang  dapat diperoleh  dari  kurva  log  ada  bermacam-macam tergantung  alat  yang  digunakan  dalam  pemboran. Umumnya data yang banyak digunakan dan sangat penting  adalah  log  GR,  log  SP,  log  densitas,  log Neutron, log resistivitas, log Caliper.

## Analisa petrofisik
Berdasarkan  data  log  dari  lapangan  X  sumur “X-1”,  dan  sumur  “X3”  maka  dapat  dilakukan analisa  kualitatif  yang  meliputi  penentuan  litologi, lapisan  permeable  sebagai  batuan  reservoir  dan penetuan jenis fluidanya.

## Zonasi
Penentuan  litologi  diawali  dengan  zonasi  pada interval  kedalaman  sumur  menurut  gabungan pembacaan  log  yang  mendeskripsikan  kesamaan jenis  batuan.  Gabungan  pembacaan  log  tersebut adalah  log  gamma  ray,  log  resistivitas,  log densitas-neutron.  Berdasarkan  hasil  zonasi  didapat
beberapa  daerah  yang  merupakan  daerah  prospek hidrokarbon.  Pada  sumur  X1  didapat  daerah prospek pada kedalaman 1272 –  1273, 1273.5-1275, 1276.5-1280, 1652-1658 mku. Kandungan Lempung (Vsh) Kandungan lempung (Vsh) dihitung
menggunakan pendekatan deterministic dari metoda log gamma ray (GR-log) dengan pemilihan GR-min  dan  GR-max masing-masing yang mencerminkan kandungan yang bersih dan kandungan  yang  kotor.  Untuk  mengatisipasi adanya  mineral  penganggu  (batuan  mengandung radioaktif  selain  shale/lempung)  dalam  pembacaan gamma  ray  log  pemilihan  GR-min  dan  GR-max dipilih  berlainan  untuk  setiap  interval/zona  dan setiap sumur.
Pada  sumur  X1  zona  A  terdapat  beberapa interval  yang  merupaka  zona  prospek  dari reservoir. Interval
kedalaman
1272-1273,1273.5-1275,  dan  1276.5-1280  pada interval  tersebut  terlihat  ada  defleksi  dari  kurva Gamma ray kearah kiri.pada kurva log densitas dan neutron  terjadi  separasi  pada  interval  tersebut terdapat  kandungan  lumping  yang  rendah.  Pada kedalam 1272 terdapat kandungan lempung sebesar 0.0553.pada kedalaman 1273.8 terdapat kandungan lempung  sebanyak  0.1531,  dan  pada  kedalaman 1276.5 terdapat kandungan lempung sebesar 0.19.

## Porositas
Ada  beberapa  metode  penentuan  porositas  dari data log yaitu porositas dari densitas, neutron, sonic dan  kombinasi  diantara  ketiganya.  Pada  penelitian ini  digunakan  perhitungan  porositas  metoda kombinasi dari Neutron-Density. Untuk mendapatkan nilai porositas efektifnya (Tabel 1).
Table 1. Nilai pororsitas pada sumur X1 Depth Rhob  D  N  DN 1274 2.2228 0.2387 0.2854 0.262 1274.3 2.1895 0.252 0.2772 0.2646 1276.6 2.2323 0.2407 0.248 0.2491 1278.3 2.3272 0.0735 0.4127 0.2447 1279.6 2.5149 0.1646 0.2596 0.1646 1600.35 2.3852 0.09849 0.2838 0.1882 1653 2.3959 0.12386 0.2738 0.17266 1655 2.3695 0.15815 0.284 0.1668 1655.6 2.4106 0.09904 0.22607 0.17052 1657 2.2392 0.2845 0.2938 0.11089

## Saturasi Fluida
Identifikasi  Mineral  Konduktif.  Pada  penelitian ini dilakukan identifikasi mineral konduktif, hal ini di  akibat  karena  adanya  kesalahan  dalam pengukuran  nilai  resistivitas  yang  terbaca  pada pembacaan log.
Identifikasi mineral konduktif dilakukan dengan menidentifiaksi  nilai  resistivitas  pada  daerah
prospek  dan  melihat  dari  hasil  uji  core  pada laboratorium.Dan berdasarkan hasil evaluasi batuan core.Dari hasil evaluasi batuan core diketahui nilai persentasi  volume  dari  mineral  konduktif  yang terdapat dibatuan reservoir.
Pada  zona  prospek  sumur-sumur  yang  di  uji memperlihatkan  ada  ketidak  normalan  dalam  nilai log resistivitasnya. Pada zona prospek terlihat nilai resistivitas  rendah.  Reservoir  resistivitas  rendah merupakan  kondisi  dimana  resistivitas  yang dihasilkan  dari  perhitungan  berada  pada  0,5  -  5 Ohm meter, sehingga evaluasi secara konvensional tidak  memberikan  nilai  dari  saturasi  air  yang akurat(Pierre.1997,).  Pada  sumur  X1  zona  A  pada kedalaman  1279  mku  nilai  dari  kurva  resistivitas menunjukan  angka  1.2871  ohm  .m.  Dari  hail  uji core  pada  kedalaman  tersebut  di  dapat  bahwa adanya  kandungan  mineral  konduktif  berupa mineral  pirit  sebesar  3  %  dan  mineral  clay  yang berupa  mineral  kaolinite  dan  chlorite.  Kandungan mineral-mineral konduktif dan clay ini mengakibatkan  terjadinya  ketidak  akuratan  dalam pembacaaan nilai resistivitas. Pada sumur X1 zona B  pada  kedalaman  1653  juga  terjadi  ketidak akuratan  dalam  pembacaan  nilai  resistivitas,  pada kedalaman  tersebut  pembacaan  nilai  resistivitas menujukan  angka  1.3112.  Dari  hasil  uji  core  pada kedalaman  tersebut  terdapat  adanya  mineral konduktif  ,  mineral  pirit  sebesar  1%  dan  mineral clay kaolinte serta chlorite.

## Evaluasi  pengaruh    mineral    konduktif  pada log resistivitas
Nilai resistivitas yang dipengaruhi oleh mineral konduktif  dibandingkan  dengan  nilai  resistivitas awal,  dari  hasilkan  perbandingan  tersebut  kita  kan mendapat pengaruh dari mineral konduktif terhadap nilai dari resistivitas pada formasi.
Dengan menghubungkan
antara
nilai konduktivitas dari mineral konduktif tersebut, pada daerah  clean  sand  kita  dapat  menggunakan  rumus Archie  untuk  menghitung  nilai  saturasi  air  (Sw). Dengan  menghilangkan  pengaruh  konduktivitas dari mineral konduktifnya.
Sebagai  contoh  perhitungan  adalah  sumur  “X1 pada  zona  2”.  Pada  kedalaman  1655.00  mku  dari hasil core di alboratorium didapat data bahwa pada kedalaman  tersebut  terdapat  mineral  konduktif yang mempengaruhi pembacaan dari log resisitivitas,  yaitu  mineral  pirit.  Dengan  adanya mineral  pirit  yang  ada  pada  batuan  reservoir mengakibatkan terjadinya penurunan nilai resistivitas pada pembacaan log. Dari hasil data uji lab  core  di  dapat  nilai  persen  volumenya  sebesar 1%.  Hubungan  antara  pengaruh  kunduktivitas mineral  dengan  volume  mineral  didapat  hubungan yang  linier.  Dengan  adanya  hubungan  yang  linier ini  maka  pengaruh  mineral  konduktif  dapat dievaluasi dengan menggunakan nilai konduktivitasnya. (Clavier.1990).
Dengan  menggunakan  hukum  kirchof  1  nilai total  dari  resistivitas  yang  terbaca  merupakan  nilai total  dari  penjumlahan  resistivitas  hidrokarbon  + konduktif. Sehingga dapat diasumsikan :
atau dengan persamaan
maka.
Diketahui :
R t : 1.2596 ohm Volume mineral konduktif: 1 %
Resistivitas mineral konduktif : 1.5 ohm
Porositas  : 0.2738 R w
: 0.0299
Nilai
Clavier,  Heim  dan  Scala  menemukan  bahwa dengan  mengasumsikan  interpretasi  dari  efek mineral  konduktif  merupakan  sebuah  model rangkaian paralalel. Sehingga
maka dimana sehingga
Dari evaluasi yang dilakukan pada nilai Rt yang dipengaruhi oleh mineral konduktif terjadi peningkatan  nilai  dari  R t .  Nilai  R t yang  terbaca pada  log  pada  awalnya  sebesar  1.2596  pada  saat dilakukan  koreksi  terhadap  nilai  resistivitas tersebut  yang di pengaruhi oleh  mineral konduktif. Nilai  resistivitas  mengalami  peningkatan  menjadi 5.2579  ohm.  Dengan  menggunakan  nilai  dari koreksi  resisitivitas  tersebut.  Nilai  resistivitas tersebut  dapat  digubakan  untuk  menghitung  S w dengan menggunakan rumus Archie:

## ( )
Dari hasil koreksi nilai resistivitas dapat terlihat adanya  perubahan  dari  nilai  S w .  Pada  perhitungan menggunakan  nilai  S w dari  hasil  pengamatan  data log  didapat  nilai  S w sebesar  =  0.8786,  dengan adanya  koreksi  terhadapat  nilai  resistivitas  yang dipengaruhi  oleh  mineral  konduktif  didapat penurunan nilai Sw dari kondisi awalnya. Nilai Sw yang  dihasilkan  dari  koreksi  nilai  mineral konduktif didapat sebesar = 0.33168.
Perhitungan  nilai  Sw  juga  dilakukan  dengan menggunakan  metode  Sw  Indonesia  Equation,  hal ini  dilakukan  untuk  perbandingan  hasil  Sw  yang dikoreksi  mineral  konduktifnya  dengan  hasil perhitungan  S w yang  dipengaruhi  oleh  mineral lempungnya. Hal ini untuk
melihat karakteristiknya.  Mineral  apa  yang  mempengaruhi dengan signifikan nilai Sw dari formasi tersebut.
Dalam menghitung  nilai  S w Indonesia  Equation perlu diketahui terlebih dahulu nilai dari Resistivitas  lempungnya  (Rsh).  Untuk  menghitung nilai resistivitas lempung digunakan nilai resistivitas  dari  pembacaan  log  pada  daerah  yang memiliki  nilai  Gamma  ray  paling  tinggi,  Rsh  = 0.7023.
Sw Indonesia Equation :

## Tekanan Kapiler dan Nilai Saturasi
Berdasarkan  hasil  perhitugan  evaluasi  logging yang  dilakukan  dengan  metode  archie  dan  metode Indonesia  Equation,  menghasilkan  nilai  Saturasi yang  beragam.  Sehingga  perlu  adanya  dasar perbandingan yang dapat menjadi tolak ukur dalam mengambil kesimpulan terhadap nilai saturasi yang menjadi acuan dalam penilai formasi tersebut.
Dari  hasil  evaluasi  data  core.  Didapat  data tekanan  kapiler  dari  formasi.  Menggunakan  2  core sebagai  sample  dalam  menguji  nilai  tekanan kapilernya.
Berdasarkan Hasil analisa batuan core terhadap hubungan  antara  nilai  tekanan  Kapiler  dengan Saturasi  air  didapat  bahwa.  Nilai  maksimal  dari tekanan  kapiler  pada  formasi  berada  pada  tekanan 35  psi.dengan  nilai  saturasi  air  yang  di  dapat sebesar  34.29%  pada  sampel  1  Hasil  uji menunjukan  batasan  dari  nilai  satursi  air  yang memungkinkan  terjadi  pada  sumur  reservoir tersebut.  Nilai  S wirr yang  didapat  34.29%  pada wilayah  tersebut  saturasi  air  tidak  dapat  mengalir karena  sudah  di  desak  oleh  minyak,  sehingga  nilai dari saturasi air menjadi minimum.
Gambar 2. Hubungan antara tekanan kapiler dan nilai saturasi sampel 1.

## IV. KESIMPULAN DAN SARAN
Berdasarkan analisis dalam pembahasan di atas, beberapa penyataan dibuat sebagai berikut:
1.  Selain  mineral  konduktif  yang  mempengaruhi nilai  resistivitas  juga  adanya  mineral  clay  pada formasi.  Untuk  mengkoreksi  adanya  mineral clay  pada  formasi  di  gunakan  Indonesia equation dan rumus Archie.
2.  Bedasarkan  hasil  evaluasi  terhadap  mineral konduktif  yang  mempengaruhi  resistivitas formasi,  didapat  penurunan  nilai  dari  saturasi air pada formasi tersebut. Pada kedalaman 1655 mku  wilayah tersebut dipengaruhi oleh  mineral konduktif  pirit  dan  mineral  clay.  Setelah dilakukan  evaluasi  terhadap  kedua  mineral tersebut  didapatkan  hasil  nilai  saturasi  airnya sebesar  0.35391,  dimana  pada  perhitungan sebelumnya  dengan  persamaan  Archie  didapat nilai sebesar 0.8558
3.  Pada  kedalaman  1655  evaluasi  menggunakan persamaan Indonesia menghasilkan nilai saturasi sebesar 0.34528
4.  Berdasarkam  hasil  evaluasi  core  didapat perbandingan  antara  nilai  tekanan  kapiler  pada sumur  dan  nilai  saturasi  sebesar  35  psi  dimana nilai  saturasinya  sebesar  0.34  dan  0.4.  Hal  ini berarti  nilai  S wirr pada  formasi  tersebut  sebesar 0.34  dan  0.43,  merupakan  nilai  minimum dimana  air  tidak  dapat  bergerak,  sehingga  itu merupakan  batas  minimum  dari  nilai  Sw  pada formasi tersebut
5.  Perhitungan  S w dengan  menggunakan  metode archie  menghasilkan  hasil  yang  lebih  akurat. Hal ini di buktikan dengan nilai yang terletak di tas nilai S wirr -nya

## DAFTAR PUSTAKA
1.   Archie.G.E., "The Electrical Resistivitas Log as an  Aid  in  Determining  Some  Reservoir
Characteristics", Petroleum Technology, Januari 1942.
2.  Anonim,  Indonesia  Basin  Summaries.  Jakarta. Patra Nusa Data, 2006.
3.  Berger.P.”Detecting  Hydrocarbon  in  Low Resistivity Environments”, Schlumberger,
1997.
4.  Bishop, Michele G., Petroleum Systems of The Northwest  Java  Province,  Java  and  Offshore Southeast  Sumatra  Indonesia,  USGS,  USA, 2000.
5.  Clavier  C.,Heim  A.,  Scala  C.”effect  of  mineral konduktif  on  resistivity  and  other  logging
measurement”.  SPWLA  17th  Annual  Logging Symposium. June 1976.
6.  Dewan,  John  T.,  “Modern  Open  Hole  Log Interpretation”, Tulsa, Oklahoma, 1983.
7.  Fickry.  2004.  Panduan  Kuliah  Gamma  Ray Log. Universitas Indonesia. Jakarta
8.  Grandis,  H., Diktat Kuliah Geo-Elektromagnet, Departemen Geofisika, FIKTM, Institut Teknologi Bandung, Bandung, 2006.
9.  Hamilton,  W.,  1979,  Tectonics  of  the Indonesian  Region.  USGS  Professional  Paper, 1078
10. Harsono  Adi,  1997,  Evaluasi  Formasi  dan Aplikasi  Log,  Edisi  revisi,  Shlumberger  Oil Services, 8 mei 1997.