Embedded Files
Location in sites.google.com: https://sites.google.com/site/lulibistadi/
Lahirnya Mud Volcano LUSI (Sejarah dan Dampak dari Bencana Mud Volcano di Jawa Timur)
LUSI-Birth of a Mud Volcano: (History and impact of Mud Volcano disaster in East Java)
Dikaji kontekstual dan diterjemahkan ke Bahas Indonesia
oleh: Dr. Ir. Hardi Prasetyo
oleh: Dr. Ir. Hardi Prasetyo
Disusun kembali 24 Mei 2008
Disusun kembali 24 Mei 2008
Dalam Rangka Memperingati 2 Tahun Lumpur Sidoarjo
Dalam Rangka Memperingati 2 Tahun Lumpur Sidoarjo
Dipersembahkan kepada Para Pahlawan LUPSI
Dipersembahkan kepada Para Pahlawan LUPSI
Suber Cnet:
Studi Pemodelan Pertumbuhan dan Potensi Geohazard Mud Volcano Lusi: Jawa Timur, Indonesia
Studi Pemodelan Pertumbuhan dan Potensi Geohazard Mud Volcano Lusi: Jawa Timur, Indonesia
Bambang P. Istadi, Gatot H. Pramono,
Prihadi Sumintadireja, Syamsu Alam
Marine and Petroleum Geology
Volume 26 (2009) 17-24-1739
Ditinjau, Dianalisis dan Dialihbahasakan ke Indonesia Oleh: Dr. Hardi Prasetyo
Modeling study of growth and potential geohazard for LUSI mudvolcano: East Java, IndonesiaOriginal Research Article
Marine and Petroleum Geology, Volume 26, Issue 9, November 2009, Pages 1724-1739
Bambang P. Istadi, Gatot H. Pramono, Prihadi Sumintadireja, Syamsu Alam
Abstract
Abstract
The mud volcano known as LUSI first erupted in May 2006 in East Java, Indonesia. The eruption has continued for over two years, and potentially will continue for many years to come, impacting an ever larger area.
An obvious and significant question is how extensive the impacted area will become in the coming years. The answer is important for planning scenarios for the relocation of people and infrastructure and for managing the environment and economy. To make such a prediction, an understanding of the geological processes controlling the mud volcanic evolution is needed.
A three-dimensional simulation model was built to predict the area affected by the mudflow over a ten-year period, with a special focus on the period from December 2007 until June 2010. The primary model inputs are the mud debit rate, the rate of subsidence and the topography.
The model prediction was validated at the December 2007 time step by comparing the results with satellite images from the same period. The simulation was found to provide a good approximation for the mud overflow and growth.
The results indicate that the mudflow tends to spread to the west and particularly to the east and north from the currently inundated area. The model predicts that in June 2010 the peak of the mud volcano will have risen 26 m above the original ground level, and the maximum subsidence will have been −63 m below the original ground level. The dynamic subsurface condition in the area creates geohazard risks, and these are also discussed in this paper.
Beranda
TINJAUAN MAKALAH DALAM BAHASA INDONESIA DILENGKAPI ANALISIS KATA KUNCI
STUDI PEMODELAN PERTUMBUHAN DAN POTENSI GEOHAZARD MUD VOLCANO LUPSI: JAWA TIMUR, INDONESIA
Bambang P. Istadi a,*, Gatot H. Pramono b,
Prihadi Sumintadireja c, Syamsu Alam d
BASELINES DATA DAN KNOWLEDGE:
TERHADAP PILIHAN OPSI STRATEGIS PENANGGULANGAN SEMBURAN, PENANGANAN LUAPAN LUPSI DAN ANTISIPASI DAMPAK BERGANDA GEOHAZARD JANGKA PANJANG
Makalah ini dirancang untuk Lusi Library: Knowledge Management
Dengan Misi menyediakan dan menyampaikan informasi dan pemahaman terhadap fenomena Lusi mud volcano baik kepada publik umum maupun Pendidikan Kebumian di Indonesia pada khususnya
Apresiasi disampaikan kepada Penulis Bambang P. Istadi
Dikontribusikan Oleh: Hardi Prasetyo
29 Mei 2009 (3 Tahun Semburan Lupsi)
DOKUMEN BASELINE 01, 2009
POKOK BAHASAN DAN KATA KUNCI
Sari
• Mud volcano Lupsi:
• Durasi semburan dan kelanjutan dampak ke depan:
• Besarnya dampak ke depan?
• Menyediakan skenario untuk penanggulangan:
• Perlunya pemahaman terhadap proses geologi yang mengendalikan evolusi mud volcano:
• Membangun simulasi model 3-D untuk 10 tahun ke depan:
• Masukan model (kecepatan semburan, subsidence dan topografi):
• Validasi citra satelit Desember 2007:
• Model lumpur mengalir ke barat, utara dan timur:
Kesimpulan
• Mud volcano Lupsi terjadi pada bidang patahan:
• Dinamika tektonik telah dan akan merusak tatanan struktur permukaan:
• Kerusakan pada infrastruktur dan disebabkan oleh pergerakan tanah dan bubble gas:
• Tidak mungkin untuk menghentikan semburan dan kesulitan Relief Well baru bila akan diterapkan:
• Perkiraan kegagalan Relief Well ke depan: Implementasi
• Fenomena subsidence yang mempersulit pemboran, kemungkinan semburan muncul kembali sepanjang bidang patahan Watukosek:
• Parameter aliran lumpur akan berhenti:
• Model untuk memperkirakan waktu dan luasnya daerah terdampak:
• Model dikembangkan berdasarkan data set dan dikalibrasi dengan citra satelit:
• Model pertumbuhan mud volcano versus subsidence terhadap waktu:
• Aktualisasi model berdasarkan faktor kecepatan semburan, subsidence, runtuhnya tanggul:
Pendahuluan
• Awal dan lokasi Semburan Lupsi:
• Fluktuasi kecepatan semburan:
• Tabel intensitas semburan Juni 2006-Desember 2008
• Penyebab dan pemicu semburan masih diperdebatkan:
• Skenario underground blowout:
• Skenario reaktivasi patahan:
• Persepsi awal hubungan semburan dan BJP-1:
• Kesimpulan tidak ada hubungan antara semburan dan BJP-1 dan rasionalisasinya:
• Skenario reaktivasi P. Watukosek, sehingga semburan tidak dapat dihentikan:
• Lupsi tidak dapat lagi dihentikan, dan diprediksi akan mengalir berabad-abad lamanya.
• Mud volcano Lupsi menyediakan informasi saat lahir sampai tahap berkembang:
• Fokus makalah pada resiko geohazard ke depan
Geologi Umum
• Referensi Pemahaman Cekungan Jawa Timur
• Cekungan busur belakang:
• Tektonik Ekstensi Tersier Awal:
• Patahan membentuk sub-basin NE-SW:
• Konfigurasi basement dan basement ridge:
• Formasi Ngimbang:
• Batugamping Terumbu:
• Tektonik Transpresi:
• Batulumpur Formasi Kalibeng dan Pucangan:
• Volkanik di Jawa sebagai sumber batuan:
• Pengendapan yang cepat:
• Zona sedimen bersifat plastis:
• Formasi Kalibeng sebagai sumber dari mud volcano di Jawa:
• Jawa Timur ideal untuk pembentukan mud volcano:
• Peran P. Watukosek untuk pembentukan mud volcano Lupsi:
Resiko Bencana Geologi (Geohazard Risks)
• Dampak semburan Lupsi:
• Antisipasi skenario kasus terburuk:
• Studi kasus mud volcano di Jatim dan Bleduk Kuwu:
• Awal Lupsi terdiri dari 5 titik semburan:
• Bubble dengan gas terbentuk di zona lemah atau reaktifasi patahan:
• Bubble berjumlah 90 bersifat berfluktuatif:
• Gas bubble mudah dan sulit terbakar:
• Terjadinya kebocoran pada penyekat:
• Bubble dengan gas bersumber dari lapisan dalam:
• Pembentukan gas secara termogenik (thermogenic origin):
• Perhatian pada daerah terjadi reaktivasi patahan:
• Pergerakan horisontal dan vertikal:
• Dampak dari perekahan terhadap masyarakat dan infrastruktur:
• Indikasi pergerakan patahan pada rel kereta api:
• Patahan menyebabkan pipa gas dan PDAM patah atau pecah:
• Patahan mengendalikan timbulnya bubble gas:
• Berlanjutnya subsidence membentuk daerah depresi diisi lumpur:
• Viskositas Lupsi yang rendah, sulit akumulasi secara vertikal:
• Lumpur menyebar ke samping dan menekan tanggul-tanggul:
• Resiko daerah terdampak meluas makin besar:
• Kegagalan tanggul ke depan makin besar:
• Masyarakat harus dipersiapkan:
• Indikator resiko geohazard dan pengendali mekanismenya:
Menghentikan Semburan
• Upaya menghentikan semburan selama ini:
• Semburan akan keluar dari bidang Patahan Watukosek:
• Keberlanjutan semburan dan intensitas:
• Kondisi dan kejadian Lupsi berdasarkan data sumur eksplorasi:
• Skenario Lupsi keluar dari bidang patahan:
• RW-1 menghadapi daerah deformasi geologi atau zona Patahan:
• Diasumsikan sumur BJP-1 tidak menembus patahan:
• Aktivitas patahan aktif menimbulkan banyaknya bubble aktif:
• Masalah pada saat reentry BJP-1:
• Dampak subsidence dan pergerakan horisontal pada pengoperasian sumur:
• Pergerakan sepanjang bidang sesar dan dampaknya:
• Semburan sulit, aliran Lupsi lama, perhatian ditujukan pada geohazard:
• Perkiraan resiko geohazard untuk perencanaan wilayah:
Dimensi Lumpur dan Volume
• Pentingnya nilai atau dimensi aliran lumpur:
• Metoda yang diterapkan:
• Genangan Lupsi Juni dan Juli 2006:
• Genangan Lupsi Mei 2007:
• Hasil Gayaberat
• Nilai ‘gravity low’ di bawah semburan:
• Batas bidang sesar:
• Perubahan density mencerminkan kondisi struktur bawah permukaan:
• Bukti pengurangan masa pada pusat semburan:
• Meningkatnya nilai anomali negatif dipicu oleh subsidence sebagai konsekuensi meningkatnya volume yang hilang:
• Anomali gayaberat negatif karena adanya perbedaan densitas dikaitkan dengan peningkatan semburan:
• Anomali positif terkait dengan pergerakan masa sepanjang bidang sesar:
Model Simulasi Pertumbuhan Lumpur
• Pengembangan pemodelan tiga dimensi:
• Parameter model:
• Metoda/tekonologi yang digunakan:
• Estimasi Lama Semburan Lumpur
• Lama semburan berdasarkan besarnya volume sumber lumpur:
• Penentuan ketebalan lampiran serpih bertekanan tinggi:
• Kurva kompaksi (compaction curve):
• Model Sejarah Penguburan (burial history):
• Kompaksi Formasi Kalibeng Atas:
• Formasi Kalibeng Bawah:
• Perkiraan ketebalan Formasi Kalibeng over pressure:
• Luas Zona Seismik dengan kecepatan rendah (Low velocity zone): 2 dan 1,3 km2
• Perhitungan lama semburan dengan asumsi semburan 100.000m3/hari:
• Konstrain prediksi lama semburan dan mekanisme aliran:
• Lama semburan terkait aliran masuk:
• Aliran tidak beraturan (fluktuasi):
• Kecepatan Semburan
• Waktu pengukuran:
• Identifikasi sumber lumpur dari kurva kompaksi:
• Identifikasi dari data seismik refleksi:
• Perhitungan kecepatan antara volume dan semburan harian: Hasil perhitungan:
• Citra satelit untuk mengkoreksi model:
• Mendapatkan angka kecepatan semburan yang realistis:
• Subsidence
• Asumsi hubungan subsidence dan semburan:
• Konsepsi penyebab semburan:
• Perhitungan kecepatan dan cakupan/luas daerah:
• Waktu pengukuran GPS:
• Besaran (magnitude) Subsidence:
• Ekspresi permukaan subsidence seperti baskom: (basin-like shape).
• Analogi dengan struktur Runtuh Porong: Diapir dan piercement sedimen overpressure
• Hasil Penting
• Volume awal dan kecepatan:
• Simulasi Pertumbuhan Lusi mud volcano: Lumpur tidak akan berekspansi ke utara, dan ke selatan karena dibatasi oleh sungai.
• Bulan Desember 2008: Lumpur akan mulai membanjiri daerah barat dimana terdapat jalan raya dan rel kereta api.
• Juni 2009, beberapa kilometer jalan akan tergenang lumpur. Pengecekan pada Desember 2008 memperlihatkan bahwa lumpur tidak membanjiri jalan raya dan rel kereta api, sebagai hasil dari implementasi perkuatan tanggul-tanggul barat.
Referensi (Selected as baseline papers)
• Abidin, H.Z., Kusuma, M.A., Istadi, B., Darmoyo, A.B., Sumintadireja, P., Purwaman, I., Andreas, H., Gamal, M., 2007. GPS-based monitoring of subsidence phenomenon in the mud extrusion areas, Sidoarjo, East-Java. In: Proceedings of the International Union of Geodesy and Geophysics (IUGG) XXIV, Perugia, Italy.
• Davies, R., Swarbrick, R., Evans, R., Huuse, M., 2007. Birth of a mud volcano: East Java, 29 May 2006. GSA Today 17, 4–9.
• Davies, R.J., Brumm, M., Manga, M., Rubiandini, R., Swarbrick, R., Tingay, M., 2008. The East Java mud volcano (2006 to present): an earthquake or drilling trigger? Earth and
• Fukushima, Y., Mori, J., Hashimoto, M., Kano, Y. Subsidence associated with the LUSI mud eruption, East Java, investigated by SAR interferometry. Journal Marine and Petroleum Geology, Special issue: Mud volcanism, submitted for Q2 publication.
• Gafoer, S., Ratman, N., 1999. Peta Geologi Regional Jawa Bagian Timur, Skala 1:500.000. Pusat Penelitian dan Pengembangan Geologi, Bandung.
• Istadi, B.P., Kadar, A., Sawolo, N., 2008. Analysis & recent study results on East Java mud volcano. In: Proceedings of the Subsurface Sediment Remobilization and Fluid Flow in Sedimentary Basin Conference, October 2008. The Geological Society, Burlington House, Piccadilly, London.
• Kopf, A.J., 2002. Significance of mud volcanism. Reviews of Geophysics 40, 1–52.
• Mazzini, A., Svensen, H., Akhmanov, G.G., Aloisi, G., Planke, S., Malthe-Sørenssen, A., Istadi, B., 2007. Triggering and dynamic evolution of the LUSI mud volcano,Indonesia. Earth and Planetary Science Letters 261, 375–388.
• Pramono, G.H., Poniman, A., 2007. The validation of mud flow model using multitemporal satellite images. In: Proceedings of the 28th Asian Conference on Remote Sensing, Kuala Lumpur, Malaysia. Pramono, G.H., Sardjono, S.P., 2007. The application of geographic information systems to measure mud dimension in Porong, Sidoarjo. Forum Geografi 21 (2).
• Sawolo, N., Sutriono, E., Istadi, B.P., Darmoyo, A.B., 2008. East Java mud volcano (LUSI): drilling facts and analysis. In: Proceedings of the African Energy Global Impact, AAPG International Conference and Exhibition, October 2008, Cape Town, South Africa.
• Sawolo, N., Sutriono, E., Istadi, B.P., Darmoyo, A.B. Mud volcano triggering controversy: was it caused by drilling? Journal Marine and Petroleum Geology, Special issue: Mud volcanism, submitted for publication. Q3
• Sumintadireja, P., Purwaman, I., Istadi, B., Darmoyo, A.B., 2007. Geology and geophysics study in revealing subsurface condition of Banjarpanji mud extrusion, Sidoarjo, East Java, Indonesia. In: Proceedings of the International Union of Geodesy and Geophysics (IUGG) XXIV, Perugia, Italy.
• Tingay, M., Heidbach, O., Davies, R., Swarbrick, R., August 2008. Triggering of the Lusi mud eruption: earthquake versus drilling initiation. Geology 36 (8), 639– 642.
Sari
Mud volcano Lupsi:
Mud volcano yang dikenal sebagai Lupsi, pertamakali menyembur 29 Mei 2006 di Jawa Timur, Indonesia.
Durasi semburan dan kelanjutan dampak ke depan:
Semburan Lupsi telah berlanjut lebih dari dua tahun, dan dapat berlanjut beberapa tahun mendatang. Sehingga berpotensi menimbulkan dampak pada daerah yang lebih luas lagi.
Besarnya dampak ke depan?
Pertanyaan mendasar adalah sebesar besar dampak yang mungkin ditimbulkan pada daerah tersebut ke depan?
Menyediakan skenario untuk penanggulangan:
Jawaban dari pertayaan di atas menurut penulis buku yang penting adalah menyediakan suatu skenario yang terbaik (the best case scenario) untuk merelokasi warga, infrastruktur, mengelola kerusakan lingkungan dan ekonomi yang ditimbulkannya.
Perlunya pemahaman terhadap proses geologi yang mengendalikan evolusi mud volcano:
Untuk membuat suatu prediksi yang lebih akurat terhadap kemungkinan yang terjadi pada masa yang akan datang, diperlukan adanya pemahaman terhadap proses-proses geologi yang mengendalikan evolusi mud volcano.
Membangun simulasi model 3-D untuk 10 tahun ke depan:
Pada makalah ini, telah dibangun suatu simulasi model 3-dimensi untuk memperkirakan daerah yang dapat tekena dampak untuk periode lebih dari 10 tahun. Dengan fokus pada perioda dari Desember 2007 sampai Juni 2010.
Masukan model (kecepatan semburan, subsidence dan topografi):
Sebagai masukan dari model utama adalah kecepatan semburan lumpur, kecepatan subsidence dan topografi.
Validasi citra satelit Desember 2007:
Prediksi dari model tersebut telah divalidasi dengan membandingkan hasil citra satelit pada perioda yang sama, yaitu bulan Desember 2007.
Model lumpur mengalir ke barat, utara dan timur:
Hasil pemodalan mengindikasikan bahwa aliran lumpur cenderung menyebar ke arah barat dan khususnya ke timur dan utara dari daerah terdapak saat ini.
Dinamika kondisi bawah permukaan di daerah terdampak menimbulkan resiko geohazard, yang akan didiskusikan dalam makalah.
KESIMPULAN
Mud volcano Lupsi terjadi pada bidang patahan:
Telah terbentuk mud volcano secara alami pada lokasi Lupsi. Berdasarkan analisis diketahui bahwa pergerakan aliran lumpur berlangsung sepanjang bidang patahan.
Secara lebih regional, telah berkembang banyak mud volcano di daerah ini (Jawa Timur dan Jawa Tengah), dan karena itu sangat besar jumlah lumpur yang telah disemburkannya.
Dinamika tektonik telah dan akan merusak tatanan struktur permukaan:
Gerakan tektonik yang konstan telah dan akan terus merusak tatanan struktur permukaan, infrastruktur, rumah dan tanggul-tanggul serta menyebabkan luapan lumpur.
Kerusakan pada infrastruktur dan disebabkan oleh bubble gas:
Pergerakan horizontal dan vertikal akan lebih merusak infrastruktur, khususnya jalan raya, pipa gas dan pipa PDAM.
Berlanjutnya reaktivasi patahan dan pembentukan rekahan dan rongga-rongga menyebabkan munculnya gas bubble pada daerah dengan radius 2-3 km dari pusat semburan Lumpur, yang berada sepanjang kelurusan dari zona lemah.
Tidak mungkin untuk menghentikan semburan dan kesulitan Relief Well baru bila akan diterapkan:
Tidak ada peluang untuk dapat menghentikan semburan lumpur. Upaya menghentikan semburan lumpur tidak mungkin untuk dapat dilakukan dengan pemboran Relief Well yang baru (sebelumnya telah dilakukan 2 RW).
Karena beberapa masalah akan timbul pada sumur relief ke depan. Terutama adanya pergerakan tanah sehingga akan menyebabkan ketidakstabilan peralatan anjungan, lokasi, selubung pemboran akan runtuh dan memutuskan tail pemboran.
Perkiraan kegagalan Relief Well ke depan: Implementasi
Bila Relief Well baru tetap akan diimplementasikan, diperkirakan akan mengalami situasi loss dan kick sehingga dapat kehilangan tekanan di dalam selubung.
Selanjutnya akan menyebabkan lobang pemboran tersumbat, oleh adanya pergerakan.
Diperkirakan tidak ada lokasi yang jelas (ideal) untuk dapat melakukan pemboran sumur Relief Well tersebut.
Fenomena subsidence yang mempersulit pemboran, kemungkinan semburan muncul kembali sepanjang bidang patahan Watukosek:
Pemboran yang dilaksanakan dari sisi luar daerah subsidence tidak efektif, karena pergerakan horizontal juga akan sangat besar.
Bila semburan lumpur dapat dihentikan pada satu lokasi, diperkirakan akan muncul di lokasi lainnya sepanjang zona patahan Watukosek.
Parameter aliran lumpur akan berhenti:
Semburan lumpur secara alami akan berhenti ketika dicapai keseimbangan antara tekanan permukaan (surface) dan masa di bawah permukaan (subsurface).
Model untuk memperkirakan waktu dan luasnya daerah terdampak:
Dalam makalah ini, suatu modal simulasi tiga dimensi telah dibangun untuk memperkirakan waktu dan luasnya daerah genangan.
Model simulasi ini dikembangkan dengan menerapkan pendekatan GIS (Geographic Information System).
Model dikembangkan berdasarkan data set dan dikalibrasi dengan citra satelit:
Sebagai masukan model berasal dari beberapa set data. Model tersebut telah dikalibrasi menggunakan data lapangan luapan lumpur dan citra satelit bulan Desember 2007.
Model pertumbuhan mud volcano versus subsidence terhadap waktu:
Pada model prediksi, diperkirakan besaran, waktu pertumbuhan mud volcano dan penenggelaman tanah.
Pengamatan secara periodik di daerah yang ditentukan sebagai daerah bahaya pada resiko geohazard diperlukan untuk memutahirkan model, khususnya dalam mengidentifikasikan pengaktifan kembali patahan sepanjang sumbu kelurusan zona patahan Watukosek dan patahan-patahan lainnya.
Aktualisasi model berdasarkan faktor kecepatan semburan, subsidence, runtuhnya tanggul:
Kecepatan lumpur, subsidence, dan runtuhnya tanggul merupakan data yang dinamis, sehingga data baru yang tersedia harus diintegrasikan pada model untuk menyempurnakan suatu perkiraan secara lebih realistik.
KESIMPULAN
Mud volcano Lupsi terjadi pada bidang patahan:
Telah terbentuk mud volcano secara alami pada lokasi Lupsi. Berdasarkan analisis diketahui bahwa pergerakan aliran lumpur berlangsung sepanjang bidang patahan.
Secara lebih regional, telah berkembang banyak mud volcano di daerah ini (Jawa Timur dan Jawa Tengah), dan karena itu sangat besar jumlah lumpur yang telah disemburkannya.
Dinamika tektonik telah dan akan merusak tatanan struktur permukaan:
Gerakan tektonik yang konstan telah dan akan terus merusak tatanan struktur permukaan, infrastruktur, rumah dan tanggul-tanggul serta menyebabkan luapan lumpur.
Kerusakan pada infrastruktur dan disebabkan oleh bubble gas:
Pergerakan horizontal dan vertikal akan lebih merusak infrastruktur, khususnya jalan raya, pipa gas dan pipa PDAM.
Berlanjutnya reaktivasi patahan dan pembentukan rekahan dan rongga-rongga menyebabkan munculnya gas bubble pada daerah dengan radius 2-3 km dari pusat semburan Lumpur, yang berada sepanjang kelurusan dari zona lemah.
Tidak mungkin untuk menghentikan semburan dan kesulitan Relief Well baru bila akan diterapkan:
Tidak ada peluang untuk dapat menghentikan semburan lumpur. Upaya menghentikan semburan lumpur tidak mungkin untuk dapat dilakukan dengan pemboran Relief Well yang baru (sebelumnya telah dilakukan 2 RW).
Karena beberapa masalah akan timbul pada sumur relief ke depan. Terutama adanya pergerakan tanah sehingga akan menyebabkan ketidakstabilan peralatan anjungan, lokasi, selubung pemboran akan runtuh dan memutuskan tail pemboran.
Perkiraan kegagalan Relief Well ke depan: Implementasi
Bila Relief Well baru tetap akan diimplementasikan, diperkirakan akan mengalami situasi loss dan kick sehingga dapat kehilangan tekanan di dalam selubung.
Selanjutnya akan menyebabkan lobang pemboran tersumbat, oleh adanya pergerakan.
Diperkirakan tidak ada lokasi yang jelas (ideal) untuk dapat melakukan pemboran sumur Relief Well tersebut.
Fenomena subsidence yang mempersulit pemboran, kemungkinan semburan muncul kembali sepanjang bidang patahan Watukosek:
Pemboran yang dilaksanakan dari sisi luar daerah subsidence tidak efektif, karena pergerakan horizontal juga akan sangat besar.
Bila semburan lumpur dapat dihentikan pada satu lokasi, diperkirakan akan muncul di lokasi lainnya sepanjang zona patahan Watukosek.
Parameter aliran lumpur akan berhenti:
Semburan lumpur secara alami akan berhenti ketika dicapai keseimbangan antara tekanan permukaan (surface) dan masa di bawah permukaan (subsurface).
Model untuk memperkirakan waktu dan luasnya daerah terdampak:
Dalam makalah ini, suatu modal simulasi tiga dimensi telah dibangun untuk memperkirakan waktu dan luasnya daerah genangan.
Model simulasi ini dikembangkan dengan menerapkan pendekatan GIS (Geographic Information System).
Model dikembangkan berdasarkan data set dan dikalibrasi dengan citra satelit:
Sebagai masukan model berasal dari beberapa set data. Model tersebut telah dikalibrasi menggunakan data lapangan luapan lumpur dan citra satelit bulan Desember 2007.
Model pertumbuhan mud volcano versus subsidence terhadap waktu:
Pada model prediksi, diperkirakan besaran, waktu pertumbuhan mud volcano dan penenggelaman tanah.
Pengamatan secara periodik di daerah yang ditentukan sebagai daerah bahaya pada resiko geohazard diperlukan untuk memutahirkan model, khususnya dalam mengidentifikasikan pengaktifan kembali patahan sepanjang sumbu kelurusan zona patahan Watukosek dan patahan-patahan lainnya.
Aktualisasi model berdasarkan faktor kecepatan semburan, subsidence, runtuhnya tanggul:
Kecepatan lumpur, subsidence, dan runtuhnya tanggul merupakan data yang dinamis, sehingga data baru yang tersedia harus diintegrasikan pada model untuk menyempurnakan suatu perkiraan secara lebih realistik.
PENDAHULUAN
Gambar. 1. Lokasi Lupsi kira-kira 30m selatan dari Surabaya. Kerangka Tektonik Regional (Regional tectonic framework) Jawa Timur (atas kanan) memperlihatkan arah patahan regional NE–SW dan E–W. Peta gayaberat Bouguer (Bouguer gravity map) diproses dengan menggunakan nilai densitas sebesar 2.65 g/cc, daerah Jawa Timur memperlihatkan adanya pusat pengendapan dari Cekungan Jawa Timur (East Java Basin’s depositional centers) ditandai dengan warna biru yang dikontrol oleh patahan utama.
Awal dan lokasi Semburan Lupsi:
Semburan lumpur terjadi 29 Mei 2006 pada 5 lokasi tersebar/terpisah satu dengan lainnya, membentuk suatu kelurusan berarah NNE-SSW, kira-kira berjarak 200 m jauhnya dari sumur Eksplorasi Banjarpanji-1 (disebut BJP-1) di Sidoarjo.
Pusat Semburan Lupsi berlokasi kira-kira 30 km selatan dari Surabaya, Jawa Timur Indonesia. Lokasi Lupsi kira-kira berjarak 700 m dari kelurusan Zona Patahan Watukosek.
Fluktuasi kecepatan semburan:
Pada bulan-bulan pertama (Mei-Juni 2006) semburan lumpur panas terjadi kecepatan aliran (flow rate) sebesar 50.000 m3/hari. September 2006 kecepatan semburan meningkat menjadi 125.000 m3/hari.
Mencapai maksimum sebesar 156.000 m3/hari pada Desember 2006. Desember 2008 menurun kembali dengan kecepatan mendekati 90.000 m3/hari.
Tabel intensitas semburan Juni 2006-Desember 2008
Penyebab dan pemicu semburan masih diperdebatkan:
Sampai saat ini, penyebab semburan lumpur panas di Sidoarjo masih terus menjadi bahan perdebatan, antara lain Mazzini et al., 2007; Davies et al., 2007, 2008; Tingay et al., 2008; Sawolo et al., 2008; Istadi et al., 2008.
Sebagai konsekuensi itu untuk mendapatkan solusinya (menghentikan atau mengurangi debit semburan) masih menjadi masalah.
Skenario underground blowout:
Bila semburan lumpur bermula dari suatu ledakan bawah permukaan (underground blowout), sehingga secara teori ia akan dapat dihentikan dengan menggunakan pemboran sumur relief (drilling relief well) untuk memotong lubang asli BJP-1. Dilanjutkan dengan memompakan lumpur berat diikuti dengan penyemenan.
Skenario reaktivasi patahan:
Namun bila semburan ternyata dikendalikan oleh suatu pengaktifan kembali patahan-patahan yang sebelumnya telah ada, dan tidak ada hubungan dengan lubang sumur BJP-1.
Sehingga sangat tidak mungkin untuk menghentikan semburan.
Persepsi awal hubungan semburan dan BJP-1:
Pada awal kejadian, karena lokasi Lupsi dekat dari lokasi pemboran BJP-1, maka sangat rasional bila masyarakat umum telah mengasumsikan bahwa semburan disebabkan oleh UGBO dari sumur BJP-1
Kesimpulan tidak ada hubungan antara semburan dan BJP-1 dan rasionalisasinya:
Setelah seluruh data set secara komprehensif diintegrasikan dan dianalisis, didasarkan data lapangan dan analisis tekanan menjadi lebih jelas bahwa antara sumur BJP-1 dan fenomena semburan lumpur tidak mempunyai hubungan satu dengan lainnya (Sawolo et al., 2008). Dengan argumen yaitu:
Pertama, tekanan rendah:
Tekanan fluida di sumur bor sangat rendah untuk merekahkan lubang bor.
Kedua, tidak ada propagasi fluida:
Tidak ada fluida di lobang bor yang berpropagasi melalui rekahan ke permukaan, saat BOP pada posisi terbuka.
Ketiga, total mati:
Lobang sumur telah terbuka dan pada kondisi total mati saat sumur menyembur dengan kecepatan 50.000 m3/hari dan hanya berjarak 200m dari semburan.
Skenario reaktivasi P. Watukosek, sehingga semburan tidak dapat dihentikan:
Penulis makalah beranggapan bahwa yang paling mungkin sebagai pemicu semburan lumpur adalah reaktivasi dari sistem Patahan Watukosek (Mazzini et al., 2007).
Sehingga kasus semburan Lupsi tidak dapat lagi dihentikan, dan diprediksi akan mengalir berabad-abad lamanya.
Mud volcano Lupsi menyediakan informasi saat lahir sampai tahap berkembang:
Selama ini telah banyak dilaksanakan studi mud volcano yang telah ada. Namun sangat sedikit diketahui kondisi sebelum dan selama tahap awal dari semburan mud volcanic.
Lupsi merupakan suatu yang khusus, karena proses-proses geologi dari saat kelahirannya hingga kondisi saat ini (tahap perkembangan runtuh menjadi kaldera) dapat diamati secara langsung.
Fokus makalah pada resiko geohazard ke depan
Makalah ini akan memusatkan perhatian pada potensi dampak dan kaitannya dengan resiko geohazard terhadap kelanjutan semburan lumpur pada tahun-tahun mendatang (potential impact and related geohazard risks from the continued mud eruption in the years to come).
GEOLOGI UMUM
Gambar 2. Peta Geologi dan sebaran mud volcano di Jawa Tengah dan Jawa Timur. Titik merah mengidentifikasikan lokasi mud volcano, beberapa diantaranya digunakan sebagai analogi pengembangan mud volcano Lupsi. (Peta dimodifikasi dari Gafoer and Ratman, 1999.)
Referensi Pemahaman Cekungan Jawa Timur
Beberapa model telah diusulkan untuk menguraikan kompleksitas Cekungan Jawa Timur (East Jawa Basin), diantaranya van Bemmelen (1949), De Genevraye and Samuel (1972), Hamilton (1979), Hall (2002), Sribudiyani et al. (2003), Smyth et al. (2005), and Prasetyadi et al. (2006).
Cekungan busur belakang:
Cekungan Jawa Timur berkembang sebagai cekungan busur belakang (back-arc-basin), merupakan hasil dari subduksi Lempeng Samudera Australia ke arah baratlaut di bawah kontinen Sunda selama akhir Kapur.
Gambar. 3. Patahan Watukosek terdiri dari 2 patahan yang sejajar dimana Kali Porong mengalami pelurusan arah aliran searah bidang patahan, gawir patahan (escarpment) Watukosek menunjukkan adanya bagian blok terangkat pada bidang sesar. Kedudukan Lupsi berada sepanjang bidang sesar Watukosek.
Tektonik Ekstensi Tersier Awal:
Suatu sistem tektonik ekstensi (extensional tectonic system) berlangsung selama Tersier Awal disebabkan oleh interaksi antara lempeng-lempeng Pasifik, Eurasia, dan Australia yang komplek. Membentuk sistem graben ekstensi (extensional graben systems) dan selanjutnya berkembang cekungan peregangan (rift basins).
Berkembang seri setengah graben dan berarah sepanjang bidang lemah Pre-Tersier NE-SW, berubah menjadi arah Timur-Barat dan lebih ke selatan.
Dua arah yang berbeda ini dikenal sebagai konfigurasi struktur dari Cekungan Jawa Timur.
Patahan-patahan dengan arah timurlaut terbentuk sebagai struktur utama, selanjutnya ditafsirkan sebagai arah patahan-patahan batas lempeng yang berlangsung pada zona subduksi dan prisma akrasi pada zaman Kapur.
Patahan membentuk sub-basin NE-SW:
Arah kedua dari patahan berkembang di selatan dari cekungan dibatasi oleh jalur lipatan utara, dimana sub-cekungan ditentukan dengan patahan-patahan berarah Timur-Barat, dan deposenter berarah sejajar dengan pulau Madura dan bagian utara pantai Jawa Timur.
Orientasi NE-SW, membentuk struktur rendahan (structural low) dari beberapa sub cekungan yang memanjang dari lepas pantai di Laut Jawa sampai ke bagian daratan Jawa Timur.
Konfigurasi basement dan basement ridge:
Konfigurasi batuan dasar (basement) dengan orientasi struktur NE-SW terdiri dari deretan punggungan basement (basement ridge) yang berkembang baik.
Dengan selingan struktur graben yang membentuk deposenter, terdiri dari sedimen Tersier.
Formasi Ngimbang:
Pengendapan sedimen klastik dan bentukan karbonat (carbonate buildup) Formasi Ngimbang berlangsung selama Eosen dan Oligosen Awal.
Sekuen Oligosen Akhir dan Miosen dipisahkan dari sekuen di bawahnya oleh ketidakselarasan, yang menyediakan suatu bentukan karbonat dengan arah umum NNE-WSW.
Batugamping Terumbu:
Perkembangan platform, dikenal sebagai Batugamping Kujung, terjadi pada Oligosen. Sedangkan Terumbu Prupuh dan Tuban berkembang pada Miosen Awal dan Pertengahan
Tektonik Transpresi:
Perioda tektonisme berawal pada Miosen Akhir berlanjut sampai Pleistosen.
Tektonisme tersebut menghasilkan rezim transpresi disebabkan oleh pergerakan lateral yang berarah Timur-Barat.
Sebagai hasil adalah terbentuknya struktur antiklin dengan orientasi timur barat.
Batulumpur Formasi Kalibeng dan Pucangan:
Sedimentasi Pliosen dan Pleistosen terdiri dari batulumpur (mudstone) yang didominasi oleh klastika volkanik Formasi Kalibeng dan Pucangan yang berpropagasi ke arah timur, dengan ketebalan antara 2400-3000m.
Volkanik di Jawa sebagai sumber batuan:
Material volkanik tersebut berasal dari sumber klastik busur volkanik Jawa di selatan daerah Sidoarjo.
Batulumpur Formasi Kalibeng pada kondisi overpressure:
Batulumpur Formasi Kalibeng berada pada kondisi tekanan tinggi/berlebih (over pressure), yang menempati sebagian besar dari cekungan.
Pengendapan yang cepat:
Satuan batuan tersebut telah diendapkan dengan cepat. Kecepatan sedimentasi dan penimbunan batuan serpih yang tebal pada deposenter terjadi selama Pliosen dan Pleistosen, menghasilkan zona ‘under-compacted shales’.
Zona sedimen bersifat plastis:
Zona tersebut sangat plastis dan mengalami tekanan berlebih karena tersedianya pemerangkapan air yang berlebih dan pematangan (maturity) material yang kaya organik.
Formasi Kalibeng sebagai sumber dari mud volcano di Jawa:
Kebanyakan mud volcano yang diketemukan di Jawa Tengah dan Jawa Timur bersumber dari satuan sedimen Formasi Kalibeng.
14 mud volcano di Jawa:
Di sekitar Lupsi sekurang-kurangnya telah diketahui terdapat 14 mud volcano yang telah diidentifikasikan di Jawa Timur dan Jawa Tengah.
Jawa Timur ideal untuk pembentukan mud volcano:
Aktif tektonik berbarengan dengan tersedianya tekanan berlebih dari sedimen matang yang kaya dengan material organik, menempatkan Jawa Timur sebagai daerah yang sangat ideal untuk berkembangnya mud volcanism.
Peran P. Watukosek untuk pembentukan mud volcano Lupsi:
Keberadaan sistem patahan Watukosek di daerah Lupsi menyediakan suatu saluran (conduit) lumpur, fluida dan gas untuk menembus berasal dari suatu horizon dalam mengalir ke permukaan, selanjutnya membentuk mud volcano Lupsi.
RESIKO BENCANA GEOLOGI (GEOHAZARD RISKS)
Dampak semburan Lupsi:
Semburan lumpur panas di Sidoarjo telah mengubur rumah, desa, sekolah, perusahaan, dan mengungsikan ribuan warga.
Potensi resiko bencana geologi (geohazard risks) berlanjut pada daerah yang berpenduduk padat serta terdapat banyak aktivitas masyarakat dan infrastruktur.
Antisipasi skenario kasus terburuk:
Studi pemodelan terhadap pertumbuhan dan potensi geohazard merupakan bagian antisipasi resiko (risks anticipation) dan perencanaan untuk skenario kasus terburuk dari mud volcano Lupsi dan daerah sekitarnya.
Studi kasus mud volcano di Jatim dan Bleduk Kuwu:
Studi mud volcano lainnya di Jawa Timur telah digunakan sebagai pengkajian bencana geologi (geohazard assessment).
Khususnya mud volcano Bleduk Kuwu yang sampai sat ini masih aktif di Jawa Tengah, terdiri dari komplek mud volcano dengan lokasi erupsi lumpur yang berganda pada radius 3,5 km dari pusat semburan utama.
Serta menunjukkan beberapa penghubung atau saluran aktif (the presence of several active conduits).
Gambar 4. Subsidence diukur dengan GPS pada 33 lokasi. Dua pengukuran GPS berkelanjutan dicirikan dengan label A dan B, diplot kecepatan subsidence rata-rata (subsidence rate) 1,9 dan 3,9 cm/hari untuk lokasi RW-1 dan RW-2. Elip warna kuning adalah penafsiran daerah yang mengalami subsidence hasil dari analisis data satelit InSAR.
Awal Lupsi terdiri dari 5 titik semburan:
Pada awalnya Lupsi terdiri dari lima semburan lumpur, selanjutnya hanya satu yang terus menyembur atau aktif hingga sekarang ini.
Sangat mungkin bahwa saluran semburan lumpur yang tidak aktif tersebut suatu saat dapat diaktifkan kembali, atau terjadi suatu semburan baru pada lokasi yang lain.
Bubble dengan gas terbentuk di zona lemah atau reaktifasi patahan:
Studi ini berpendapat bahwa kemungkinan semburan lumpur atau yang lebih dikenal sebagai lokasi gas bubble atau lokasi baru semburan terjadi sepanjang zona lemah atau patahan-patahan yang diaktifkan kembali atau suatu zona sesar baru (new fault zone).
Bubble berjumlah 90 bersifat berfluktuatif:
Bubble-bubble yang menyemburkan gas metan telah diidentifikasi berjumlah lebih dari 90 lokasi, dan umumnya berasosiasi dengan rekahan-rekahan.
Kondisi bubble tersebut berfluktuatif, beberapa diantaranya lebih aktif daripada yang lainnya, sedangkan yang lainnya mengecil atau bahkan telah mati.
Gas bubble mudah dan sulit terbakar:
Pada kebanyakan kasus, kandungan gas metan tersebut tidak bersifat mudah terbakar (non-flammable), karena hanya mengandung konsentrasi yang rendah dari pemisahan gas yang cepat di udara.
Namun pada bubble dengan semburan yang kuat, maka konsentrasi gas metan cukup tinggi dan mudah terbakar.
Terjadinya kebocoran pada penyekat:
Rembesan gas dari rekahan-rekahan ini memberikan kepercayaan terjadinya kebocoran penyekat dan hilangnya kapasitas penyekatan dari patahan-patahan yang ada.
Di samping itu serpih pada struktur geologi yang menutupinya, mengandung akumulasi gas.
Bubble dengan gas bersumber dari lapisan dalam:
Terdapatnya bubble gas juga memberikan kepercayaan bahwa subsidence bukan merupakan fenomena yang berasal dari dekat permukaan.
Sebagai hasil dari pembebanan permukaan oleh berat dari lumpur atau kompaksi dari lapisan tanah lunak.
Tapi juga sebagai pengaruh dari lapisan yang lebih dalam di bawah permukaan bumi, dimana gas berasal dari sumber yang dalam.
Pembentukan gas secara termogenik (thermogenic origin):
Kromatografi dari contoh gas bubble berasal dari dekat pusat kawah semburan yang diambil Juli 2006 mengindikasikan bahwa komposisinya terutama terdiri dari metan.
Namun keberadaannya berasal dari gas berat dengan kuantitas yang kecil termasuk etana, propan, butana, dan pertan.
Ditafsirkan gas-gas tersebut berasal dari termogenik (thermogenic origin) dan mengalir melalui rekahan-rekahan yang panjang.
Perhatian pada daerah terjadi reaktivasi patahan:
Daerah yang secara khusus patut mendapatkan perhatian adalah di daerah yang mengalami reaktivasi patahan-patahan, dimana pergerakan diferensial telah menciptakan shear stress.
Pergerakan horisontal dan vertikal:
Daerah di dalam radius 2-3 km dari pusat semburan telah mengalami pergerakan horisontal dan vertikal yang sejajar dengan arah bidang patahan.
Secara umum pergerakan tersebut mempunyai variasi spasial dan temporal dalam besaran dan arahnya.
Dampak dari perekahan terhadap masyarakat dan infrastruktur:
Sebagai akibat dari terjadinya perekahan dari reaktivasi patahan-patahan yang sebelumnya telah eksis adalah rusaknya rumah warga, bangunan-bangunan dan infrastruktur yang berada di dekatnya.
Indikasi pergerakan patahan pada rel kereta api:
Pergerakan sesar tersebut adalah dextral (pergerakan menganan), yang dapat diamati dari pergerakan rel kereta api yang terjadi pada September 2006. Di mana pergerakan patahan berlangsung lebih besar dari pergerakan subsidence yang diukur dengan GPS.
Patahan menyebabkan pipa gas dan PDAM patah atau pecah:
Pecahnya pipa gas dan pipa PDAM pada lokasi yang sama mendukung bahwa pergerakan berlangsung sepanjang bidang patahan.
Patahan mengendalikan timbulnya bubble gas:
Berlanjutnya pergerakan sepanjang patahan-patahan juga tampaknya menyebabkan munculnya lebih banyak lagi bubble gas.
Berlanjutnya subsidence membentuk daerah depresi diisi lumpur:
Berlanjutnya subsidence telah membentuk suatu daerah depresi berbentuk baskom atau kerucut, yang terjadi sebagai akibat berlanjutnya semburan lumpur diikuti pembebanan.
Terjadinya subsidence telah membentuk suatu ruang, cekungan alami yang selanjutnya diisi lumpur.
Viskositas Lupsi yang rendah, sulit akumulasi secara vertikal:
Tingginya kandungan air dari lumpur yaitu dengan perbandingan antara 70-50% air terhadap lumpur, berarti bahwa lumpur mempunyai viskositas yang lebih rendah.
Sehingga secara vertikal tidak dapat mengendap untuk membentuk suatu struktur gunung yang tinggi (steep mountain-like structure).
Lumpur menyebar ke samping dan menekan tanggul-tanggul:
Lumpur cenderung menyebar ke samping sehingga bertambahnya tekanan pada lumpur akan mempengaruhi tanggul-tanggul, sampai pada kondisi dapat runtuh dan menimbulkan banjir.
Resiko daerah terdampak meluas makin besar:
Daerah di sekitar area terdampak dihadapkan pada suatu resiko banjir lumpur yang makin luas.
Tanggul-tanggul dan kolam penampungan lumpur sementara pada situasi darurat.
Tanggul tersebut dibangun hanya dari tanah dan material batuan sehingga tidak stabil.
Jebolnya tanggul-tanggul yang terjadi pada masa lalu telah menyebabkan banjir lumpur panas (temperatur 85-95oC) dan menyebabkan kecelakaan manusia atau kepunahan binatang.
Kegagalan tanggul ke depan makin besar:
Ke depan kegagalan tanggul berpotensi terjadi saat lumpur tumbuh dan meningkatkan tekanan hidrostatik pada dinding tanggul.
Masyarakat harus dipersiapkan:
Saat resiko kegagalan tanggul adalah besar, maka sebagai alternatif rute alternatif transportasi publik harus disiapkan untuk memelihara keberlanjutan akses.
Indikator resiko geohazard dan pengendali mekanismenya:
Dalam aspek resiko geohazard, bukti-bukti dan perhatian daerah mencakup:
• Retak dan pecahnya pipa gas dan pipa PDAM (shear and subsidence);
• Melengkungnya rel kereta api (shear/faulting and subsidence);
• Retak jalan (subsidence);
• Gangguan integritas dari selubung relief well (subsidence and shearing), 5) runtuh tanggul (subsidence);
• Bubble gas tampaknya terjadi di sepanjang rekahan atau sona lemah (shear/faulting and subsidence).
MENGHENTIKAN SEMBURAN
Upaya menghentikan semburan selama ini:
Beberapa upaya dengan menggunakan berbagai metoda dan teknik untuk menghentikan semburan telah dilaksanakan.
Sekurang-kurangnya pada masa Timnas PSLS empat metoda yang telah diimplementasikan untuk menghentikan semburan.
4 Metoda yang telah diterapkan yaitu:
(1) Snubbing Unit (re-entry of the Banjarpanji-1 well),
(2) Well Side Tracking,
(3) Relief Wells, dan
(4) Rangkaian bola-bola beton dengan densitas tinggi High Density Chained Balls (HDCB).
Upaya ini tanpa sukses dan telah berhadapan dengan isu keamanan dan tantangan bawah permukaan.
Keberhasilan menghentikan semburan sangat kecil:
Keberhasilan untuk menghentikan semburan sangat kecil, bahkan bila semburan dapat dihentikan dengan sukses.
Akan sangat mungkin bahwa semburan lumpur yang baru dapat terjadi dimana saja.
Semburan akan keluar dari bidang Patahan Watukosek:
Diperkirakan bahwa semburan lumpur lainnya (baru) akan bermula sepanjang bidang patahan Watukosek (kira-kira 1 km) ke selatan dari Lupsi.
Keberlanjutan semburan dan intensitas:
Aliran lumpur terus berlangsung kuat. Saat ini kecepatan semburan mendekati 90.000 m3/hari (Agustus 2008).
Baseline informasi dari RW-1 dan RW-2 untuk menghentikan semburan:
Informasi geologi penting tentang Lusi berasal dari relief wells drilled (RW-1 dan RW-2) untuk menghentikan semburan lumpur.
Kondisi dan kejadian Lupsi berdasarkan data sumur eksplorasi:
Operasi sumur eksplorasi dilaksanakan antara bulan Juni 2006 dan Januari 2007, dan menghasilkan kondisi dinamika geologi dari Lusi mud volcano yang didiskusikan di bawah ini.
Skenario Lupsi keluar dari bidang patahan:
Erupsi lumpur dapat dari bidang patahan, bukan dari lubang bor.
RW-1 menghadapi daerah deformasi geologi atau zona Patahan:
RW-1 menembus suatu bidang patahan dan saat pemboran diteruskan lumpur pemboran hilang dan terjadi kicks.
Seperti kondisi di bawah permukaan, jaringan patahan telah diaktifkan kembali dan dibentuk.
Diasumsikan sumur BJP-1 tidak menembus patahan:
Sedangkan sumur BJP-1 tidak menembus patahan dan fluida pemboran berhasil bersirkulasi setelah terjadinya semburan lumpur.
Aktivitas patahan aktif menimbulkan banyaknya bubble aktif:
Meningkatnya jaringan patahan aktif menghasilkan semakin banyak bubble gas di daerah tersebut.
Beberapa diantaranya berkembangnya di dekat dari sumur relief dimana memperlihatkan kondisi yang membahayakan terhadap keamanan pengoperasian operator dan peralatan untuk operasi relief well.
Masalah pada saat reentry BJP-1:
Pada awal dari re-entry sumur BJP-1, operasi tidak ada masalah, namun penurunan dan pergerakan horisontal bermula, sehingga permasalahan mulai timbul.
Dampak subsidence dan pergerakan horisontal pada pengoperasian sumur:
Subsidence dan pergerakan horisontal telah menyebabkan tight hole, key seating, equipment failure, casing collapse and buckling, moving upwards from deep to shallow zones in BJP-1 well as well as RW-1 and RW-2 wells. Even a wellhead drop was observed.
1) Integritas casing sulit untuk dicapai saat deformasi bertambah.
2) Pergerakan termasuk horizon dangkal dan dalam saat casing runtuh;
3) Keamanan pengoperasian dari sumur relief tergantung pada integritas dari tanggul-tanggul melindungi lumpur terhadap banjir di daerah operasi sumur pemboran.
Pergerakan sepanjang bidang sesar dan dampaknya:
Namun ketika tekanan pada tanggul-tanggul bertambah, ia dapat runtuh sehingga menyebabkan banjir lumpur panas.
Pergerakan yang berbeda sepanjang patahan telah menyebabkan rekahan dan pergerakan menghasilkan tanggul-tanggul runtuh pada beberapa waktu.
Semburan sulit, aliran Lupsi lama, perhatian ditujukan pada geohazard:
Ketika tidak ada indikasi bahwa semburan lumpur akan berhenti secepatnya dan upaya untuk menghentikan aliran lumpur tidak sukses, sebagai konsekuensi pemahaman terhadap geohazard risks menjadi sangat penting.
Perkiraan resiko geohazard untuk perencanaan wilayah:
Perkiraan terhadap resiko dan luas daerah yang dapat terkena dampak dari semburan yang berlangsung dan penurunan menjadi bagian penting didalam perencanaan untuk masyarakat dan prasarana.
Studi dan survei diperlukan untuk memahami perilaku dari mud volcano ini dan di dalam membangun model perkiraan.
DIMENSI LUMPUR DAN VOLUME
Pentingnya nilai atau dimensi aliran lumpur:
Memahami dimensi aliran lumpur menjadi bagian penting untuk menghitung kecepatan erupsi, hal ini merupakan suatu kunci masukan pada model simulasi.
Metoda yang diterapkan:
• Pelaksanaan survei dilakukan dengan melakukan 51 lokasi pengukuran pada daerah genangan lumpur tahun 2006. Ketinggian aliran lumpur diukur menggunakan pita dan membandingkan dengan daerah sekitarnya yang tidak tergenang. Luas genangan dari daerah banjir dihitung menggunakan citra satelit IKONOS.
• Luas dan ketebalan digunakan untuk menghitung volume dan merata-ratakan kecepatan semburan selama interval lebih dari satu bulan.
• Survei dilangsungkan pada Mei 2007, dengan 42 titik survei dalam upaya untuk mengkuantitatifkan perubahan dalam dimensi lumpur, volume dan kecepatan erupsi satu tahun setelah awal semburan.
• Volume lumpur telah dihitung dengan program GIS yang mampu untuk melakukan analisis 3-D. Volume dihitung menggunakan metoda cut and fill dimana volume dihitung ketika material telah dipindahkan atau dihasilkan.
HASIL
Genangan Lupsi Juni dan Juli 2006:
Untuk bulan Juni dan Juli 2006, daerah genangan lumpur mencapai 111 dan 179 hektar, perhitungan volume adalah 1,1 juta m3 dan 2,5 juta m3.
Genangan Lupsi Mei 2007:
Daerah genangan lumpur telah meningkat menjadi 628 ha walaupun telah dibangun tanggul-tanggul, atau lima kali lebih besar dari daerah genangan Juni 2006, perhitungan volume lumpur meningkat menjadi 37,3 juta m3.
Hasil Gayaberat
Nilai ‘gravity low’ di bawah semburan:
Sangat menarik adalah adanya lapisan yang mempunyai angka rendah di sekitar kawah semburan Lupsi dan dengan nilai tinggi di bagian timurlaut.
Rendahnya nilai anomali The low residual Bouguer anomaly values mencakup daerah seluas 1,3 km2 dan ditafsirkan sebagai dampak struktur bawah permukaan yang dangkal shallow subsurface structural effect.
Batas bidang sesar:
Batas antara nilai anomali Bouguer yang rendah dan tinggi ditafsirkan sebagai suatu zona patahan dengan komponen slip yang normal/strike
Perubahan density mencerminkan kondisi struktur bawah permukaan:
Perubahan bulk density ditafsirkan sebagai perubahan kondisi bawah permukaan karena tekanan dan perubahan densitas dipercaya penyebaran masa di bawah tanah disebabkan oleh pergerakan dari masa, fluida dan gas.
Bukti pengurangan masa pada pusat semburan:
Angka gayaberat mikro yang rendah mengindikasikan adanya pengurangan masa pada sekitar pusat semburan lumpur pada kedalaman 1700 m, sesuai dengan data set terdahulu.
Gambar 5. Peta gayaberat regional Anomali (Bouguer anomaly map) memperlihatkan terdapatnya deposenter sebagaimana dicirikan oleh nilai Bouger yang rendah (lower Bouguer values) di timurlaut Lupsi. Data residual gravity untuk Lupsi (kanan) dan daerah sekitarnya mencirikan suatu gayaberat rendah (gravity low), dimana ditafsirkan sebagai subsidence disebabkan oleh menghilangnya lumpur di daerah tersebut, dan adanya patahan-patahan pada daerah dengan kontras densitas yang tinggi (high density contrast areas) (perubahan signifikan pada kontur). Pergerakan horizontal dan vertikal mencirikan arah dan besaran pergerakan (surface movement direction and magnitude).
Meningkatnya nilai anomali negatif dipicu oleh subsidence sebagai konsekuensi meningkatnya volume yang hilang:
Pada Lupsi meningkatnya nilai anomali negatif sebagaimana dilihat dari peta microgravity time-lapse antara Agustus 2006-Agustus 2006, dipercaya adanya peningkatan subsidence disebabkan oleh peningkatan volume lumpur yang hilang.
Gambar 6. Peta anomali gayaberat-mikro selang waktu (Microgravity time-lapse anomaly map), Perioda bulan September 2006–Augustus 2006 (kiri) dan perioda Oktober 2006–Augustus 2006 (kanan).
Perubahan pada pola gayaberat mikro (microgravity patterns) terhadap waktu memberikan kepercayaan terhadap adanya penyebaran kembali masa dan mobilisasi lumpur (redistribution and mud remobilization), atau aliran fluida serpih yang bertekanan lebih (fluidized over-pressured shale flowing) diperlihatkan dengan panah merah ke dalam conduit yang plastis sepanjang kelurusan Patahan Watukosek.
Angka gayaberat mikro positif (Positive microgravity values) menunjukkan bahwa masa meningkat dan ada pemasukan (mass increase and recharge), dimana angka negatif memberikan dugaan subsidence disebabkan oleh hilangnya lumpur atau fluida (mud/fluid withdrawal).
Bila daerah pegunungan di baratdaya dan daerah baratlaut sebagai daerah imbuhan (recharged) dari dinamika hidrogeologi Lupsi, maka aliran lava akan berlangsung pada perioda yang lama.
Gambar 6. Peta anomali gayaberat-mikro selang waktu (Microgravity time-lapse anomaly map), Perioda bulan September 2006–Augustus 2006 (kiri) dan perioda Oktober 2006–Augustus 2006 (kanan).
Perubahan pada pola gayaberat mikro (microgravity patterns) terhadap waktu memberikan kepercayaan terhadap adanya penyebaran kembali masa dan mobilisasi lumpur (redistribution and mud remobilization), atau aliran fluida serpih yang bertekanan lebih (fluidized over-pressured shale flowing) diperlihatkan dengan panah merah ke dalam conduit yang plastis sepanjang kelurusan Patahan Watukosek.
Angka gayaberat mikro positif (Positive microgravity values) menunjukkan bahwa masa meningkat dan ada pemasukan (mass increase and recharge), dimana angka negatif memberikan dugaan subsidence disebabkan oleh hilangnya lumpur atau fluida (mud/fluid withdrawal).
Bila daerah pegunungan di baratdaya dan daerah baratlaut sebagai daerah imbuhan (recharged) dari dinamika hidrogeologi Lupsi, maka aliran lava akan berlangsung pada perioda yang lama.
Anomali gayaberat negatif karena adanya perbedaan densitas dikaitkan dengan peningkatan semburan:
Anomali negatif juga memperlihatkan adanya densitas kontras dengan daerah sekitarnya, dari 150-200 mGal sangat konsisten dengan adanya peningkatan kecepatan semburan lumpur dari 50.000 menjadi 100.000 m3/hari, September 2006.
Anomali positif terkait dengan pergerakan masa sepanjang bidang sesar:
Anomali positif dengan sebaran dengan arah NE-SW sangat selaras dengan kelurusan patahan di daerah ini menunjukkan bahwa anomali tersebut terkait dengan pergerakan aliran lumpur sepanjang bidang patahan.
Variasi anomali terutama searah dengan sistem Patahan Watukosek yang berarah NNE-SSW
MODEL SIMULASI PERTUMBUHAN LUMPUR
Pengembangan pemodelan tiga dimensi:
Suatu model simulasi tiga dimensi telah dibangun untuk memperkirakan daerah yang terdampak oleh aliran lumpur selama 10 tahun, dengan perhatian khusus pada bulan Desember 2007 sampai Juni 2010.
Parameter model:
Parameter model adalah 1) kecepatan semburan (rate of eruption), 2) kecepatan subsidence. Dan 3) topografi.
Metoda/tekonologi yang digunakan:
Pola subsidence berhubungan dengan sebaran patahan di daerah, sehingga beberapa data digunakan untuk menentukan termasuk 1) VLF, 2) resistivity, 3) microgravity anomaly changes, dan 4) daerah subsidence ditentukan dengan perkiraan subsidence dari InSAR data.
Estimasi Lama Semburan Lumpur
Lama semburan berdasarkan besarnya volume sumber lumpur:
Lamanya semburan lumpur diperkirakan dari volume formasi serpih bertekanan berlebih (over-pressure shale formation) yang dapat menyembur ke permukaan.
Penentuan ketebalan lampiran serpih bertekanan tinggi:
Ketebalan dari penampang serpih dengan tekanan tinggi ditentukan dari kurva kompaksi yang telah dikembangkan dari data sonic log, menggunakan formula the Reiga–Clemenceau equation (Issler, 1992):
4 ¼ 1 _ tma
tlog!1=x
Dimana: 1) 4 ¼ porosity, 2) tma ¼ transit time of the matrix (67.05 ms), 3) tlog ¼ transit time from the wire line log and 4) x ¼ acoustic formation factor (2.19).
Fig. 7. Hasil awal model simulasi lumpur untuk Desember 2007 (garis putus oranye) divalidasi dengan pengamatan dari citra satelit Crisp Januari 2008.
Model terhadap daerah genangan yang diperkirakan dan disini masukan model telah disesuaikan untuk mendapatkan keseusian dengan realitas lapangan.
Survey VLF (kanan atas) sepanjang lintasan garis merah, sedangkan survei ‘geolistrik’ (resistivity survey) pada kanan bawah sepanjang lintasan dengan warna kuning.
Kedua survei mengidentifikasikan terdapatnya patahan dangkal (shallow faults) yang tertutup oleh tanah atas. Patahan ini selaras dengan Sistem Patahan Watukosek.
Kurva kompaksi (compaction curve):
Kurva kompaksi di sumur BJP-1 (gambar 8) dapat dibagi menjadi tiga bagian. Bagian paling atas terdiri dari formasi Kabuh, Pucangan1 dan Pucangan 2 dari 0-980m.
Model Sejarah Penguburan (burial history):
Karena data tidak tersedia maka model sejarah penguburan menggunakan kurva kompaksi dari Sclater and Christie (1980).
Kompaksi Formasi Kalibeng Atas:
Bagian tengah terdiri dari Formasi Kalibeng Atas, dengan litologi didominasi dengan serpih, dan kurva kompaksi jatuh pada kategori formasi di bawah nilai kompaksi (under compacted) atau tekanan berlebih (over-pressured).
Gambar 8. Kurva kompasi (compaction curve) dari sumur Banjarpanji-1 memperlihatkan terdapatnya serpih pada kondisi under-compacted dan over-pressured shale pada Formasi Kalibeng 1 Atas.
Secara sismik, merupakan suatu sona dengan kecepatan rendah (low velocity zone) sebagaimana diperlihatkan pada profil interval velocity.
Zone ini merupakan sumber lumpur, dan ketebalannya telah digunakan untuk memperkirakan volume dari lumpur yang berpotensi menyembur.
Formasi Kalibeng Bawah:
Bagian lebih bawah terdiri dari Formasi Kalibeng 2, dimana litologi didominasi oleh batupasir volkanik, dan kurva kompaksi jatuh pada kisaran kompaksi normal (normal compaction).
Perkiraan ketebalan Formasi Kalibeng over pressure:
Perkiraan ketebalan dari Formasi Kalibeng 1 over pressured 640m. Zona tekanan berlebih ini didominasi oleh serpih antara interval 1860-1220m di sumur BJP berkorelasi dengan zona kecepatan rendah (low velocity) sebagaimana diperlihatkan oleh penampang interval velocity dari data BJP-1.
Gambar 10. Peta basis data seismik refleksi di overlay dengan citra satelit. Perbedaan garis lintasan menunjukkan perbedaan umur dan garis hitam adalah penampang seismik lintasan pada Gambar 9.
Daerah Lusi diidentifikasikan dengan peta residual gravity dari Gambar 5. Sebagai catatan daerah struktur Porong collapse di beri warna bayangan, dengan mencakup sebagaian daerah Lusi dibandingkan dengan Lusi ditafsirkan pada darah over-pressured shale area.
Porong berbatasan dengan mud volcano yang telah digunakan sebagai analogi untuk membangun model subsidence Lusi.
Luas Zona Seismik dengan kecepatan rendah (Low velocity zone):
Secara seismik zona kecepatan rendah diperlihatkan pada gambar 9. Perhitungan data diperkirakan 2 km2 dan 1,3 km.
Perhitungan lama semburan dengan asumsi semburan 100.000m3/hari:
Perhitungan lamanya semburan Lupsi adalah volume (ketebalan dikali dengan rata-rata daerah), dengan asumsi kecepatan semburan 100.000 m3/hari, menunjukkan durasi semburan dapat bervariasi antara 23 sampai 35 tahun.
Konstrain prediksi lama semburan dan mekanisme aliran:
Prediksi tersebut tergantung dari volume masa di bawah permukaan dan kondisi tekanan.
Dalam hal ini faktor-faktor seperti serpih overpressure, kegempaan (seismicity) dan tektonik mempengaruhi lamanya semburan dan mekanisme aliran.
Lama semburan terkait aliran masuk:
1) Skenario cepat berhenti: Aliran akan berhenti secepatnya bila tidak ada pemasukan atau imbuhan air dari luar, dan 2) Skenario lama berhenti: akan menjadi lama bila sistem berhubungan dengan suatu imbuhan air.
Aliran tidak beraturan (fluktuasi):
Kemungkinan lain bila aliran tidak beraturan, dimana aliran lumpur berhenti, tekanan yang terbangun dan meluasnya penutup, maka aliran akan mulai lagi.
KECEPATAN SEMBURAN
Waktu pengukuran:
Kecepatan semburan dihitung berdasarkan pada dimensi lumpur yang dilaksanakan 2006 dan 2007.
Identifikasi sumber lumpur dari kurva kompaksi:
Gambar 8 kurva kompaksi dari sumur BJP-1 memperlihatkan hadirnya zona serpih dari Formasi Kalibeng 1 Atas yang berada di bawah kompaksi (under-compacted) dan mengalami over-pressured.
Identifikasi dari data seismik refleksi:
Secara seismik, zona tersebut merupakan lapisan-lapisan dengan kecepatan rendah sebagaimana diperlihatkan pada penampang interval velocity.
Zona ini ditentukan sebagai sumber dari lumpur, dan besarnya ketebalan telah digunakan sebagai alat bantu untuk memperkirakan volume lumpur yang berpotensi menyembur.
Perhitungan kecepatan antara volume dan semburan harian:
Kecepatan semburan lumpur dihitung dengan membagi antara volume lumpur dengan jumlah hari semburan.
Hasil perhitungan:
Hasil perhitungan kecepatan semburan untuk bulan Juni dan Juli 2006 sebesar 50.785 m3/hari dan 44.671 m3/hari, sedangkan kecepatan semburan dari bulan Mei 2007 dengan dimensi sebesar 111.042m3/hari, lebih dua kali lebih besar dari tahun 2007.
Sehingga angka tersebut digunakan untuk melakukan Pemodelan.
Citra satelit untuk mengkoreksi model:
Data citra satelit IKONOS Januari 2008 digunakan untuk mengkoreksi model.
Mendapatkan angka kecepatan semburan yang realistis:
Disimpulan bahwa kecepatan 111.042m3/hari terlalu besar. Maka selanjutnya kecepatan 90.000 m3/hari memberikan angka yang lebih realistis.
Subsidence
Asumsi hubungan subsidence dan semburan:
Subsidence akan terus berlanjut bila semburan lumpur masih berlangsung.
Konsepsi penyebab semburan:
Subsidence kemungkinan merupakan hasil dari kombinasi dari releksasi tanah karena aliran lumpur, pembebanan karena berat lumpur menyebabkan daerah makin mengalami kompaksi, perpindahan tanah, transformasi geologi struktur dan tektonik (Abidin et al., 2007).
Perhitungan kecepatan dan cakupan/luas daerah:
Kecepatan subsidence dihasilkan terutama dari pengukuran GPS, sedangkan luas dari subsidence dihasilkan dari data InSAR. Pola subsidence juga dipengaruhi oleh sistem patahan di daerah (dalam hal ini Sistem Patahan Watukosek).
Waktu pengukuran GPS:
Pengukuran subsidence dari GPS antara 29 Juni 2006 sampai 12 April 2007 digunakan dalam model.
Besaran (magnitude) Subsidence:
1) Tertinggi: Besaran (magnitude) Subsidence tertinggi sebesar 5,53cm/hari terjadi berada di dekat dari pusat semburan. 2) Terendah: kecepatan terendah sebesar 0,04 cm/hari terkonsentrasi di luar daerah subsidence sebagaimana diidentifikasikan dari citra InSAR.
Ekspresi permukaan subsidence seperti baskom:
Ekspresi permukaan dari subsidence ditentukan dengan bentuk baskom atau bentuk seperti cekungan (basin-like shape).
Seperti halnya struktur runtuh Porong (Porong collapse structure) yang berjarak 7 km dari daerah Lusi.
Analogi dengan struktur Runtuh Porong:
Struktur Porong Collapse merupakan mud volcano purba (paleo mud volcano) yang memperlihatkan kawah, pipa penghubung, dan beberapa patahan menghasilkan slumping secara gayaberat (gravitational slumps) dan struktur intrusi.
Sehingga diperkirakan berkembang sebagai struktur piercement (pembubungan) pergerakan ke atas dari sedimen over-pressure di bawah diapir serpih (shale diaper).
HASIL PENTING
Volume awal dan kecepatan:
Volume awal lumpur adalah 37,3 juta m3 dan kecepatan aliran 90.000 m3/hari.
Simulasi model daerah lumpur dan volume sampai Juni tahun 2010 diperlihatkan pada Tabel 1.
Untuk perioda 6 bulan, ekspansi area ditutupi lumpur diperkirakan antara 25 sampai 292 hektar.
Berlanjutnya volume lumpur dengan masukan yang konstan. Sehingga Juni 2020 diprediksi volume sebesar 136 juta m3 dan maksimum pada daerah terdampak sebesar 1448 juta m3.
Gambar 11 memperlihatkan perluasan daerah kemungkinan terdampak oleh Lumpur.
Dari bulan Desember 2007 sampai Juni 2008, lumpur akan menyebar ke arah timur, karena daerah ini relatif datar dan tidak mempunyai infrastruktur yang signifikan.
Banjir dapat terjadi di daerah barat. Jalan raya dan rel kereta api berada di sisi ini, sehingga prioritas penanggulangan harus diberikan pada tanggul di bagian barat.
Gambar 11. Perkiraan pertumbuhan lumpur dari Desember 2007 sampai 2010. Catatan lebar penyebaran (lebar) lumpur adalah ke timur dan ke barat, karena viskositas lumpur rendah dan topografi datar (flat).
Sebaran ke utara dan selatan terhambat oleh adanya tanggul atau pematang Kali Porong dan Ketapang.
Simulasi memperlihatkan bahwa lumpur tidak akan berekspansi ke utara, dan ke selatan karena dibatasi oleh sungai.
Bulan Desember 2008, lumpur akan mulai membanjiri daerah barat dimana terdapat jalan raya dan rel kereta api. Tanggul bagian barat harus diperkuat untuk menjaga agar lumpur tidak membanjiri jalan dan rel.
Juni 2009, beberapa kilometer jalan akan tergenang lumpur. Karena itu direkomendasikan bahwa jalan raya dan rel alternatif harus disiapkan (antara lain direlokasi) untuk mengantisipasi kehilangan akses karena banjir lumpur.
Pengecekan pada Desember 2008 memperlihatkan bahwa lumpur tidak membanjiri jalan raya dan rel kereta api, sebagai hasil dari implementasi perkuatan tanggul-tanggul barat.
Gambar 12 memperlihatkan perkiraan evolusi pandangan penampang. Puncak relatif stabil pada ketinggian 21 m dan mencapai ketinggian maksimum 26 m dekat pusat semburan.
Tingginya kecepatan subsidence di pusat semburan dicerminkan oleh penenggelaman dari permukaan subsidence sebelumnya yang mencapai 63 m pada Juni 2010.
Kecepatan subsidence diperkirakan 9 m per 6 bulan sampai Desember 2008, lalu berkurang menjadi 8 m/6 bulan sampai 2010.
Gambar 13 memperlihatkan perkiraan pertumbuhan dari 4,1 dan 6,4 km pada Desember 2007 mendekati 6,4 menjadi 8,4 km Juni 2010.
Gambar 12. Pandangan penampang memperlihatkan prediksi pertumbuhan lumpur dan subsidence dari Desember 2007 sampai Juni 2010.
Subsidence terbesar didapatkan dekat dari semburan utama, dan membentuk suatu struktur depresi yang menyerupai kerucut. depression-cone-like structure.
Gambar 13. Subsidence yang berlangsung, ketika ia membentuk suatu ruang kosong secara alami diisi lumpur tapi juga menyebabkan kerusakan.
Perkembangan daerah subsidence diperlihatkan denga lingkaran merah, luas 6,4 km dari utara ke selatan pada simulasi Desember 2007. Pada bulan Juni 2010, diperkirakan akan memperluas 8,4 km.
Daerah yang dilindungi oleh tanggul-tanggul kemungkinan tidak tergenang oleh lumpur namun akan mengalami penurunan ke depan.
Banyak rumah-rumah pada daerah subsidence akan mengalami kerusakan seperti dinding retak karena danya pergerakan selama proses penurunan.
Referensi
Abidin, H.Z., Kusuma, M.A., Istadi, B., Darmoyo, A.B., Sumintadireja, P., Purwaman, I., Andreas, H., Gamal, M., 2007. GPS-based monitoring of subsidence phenomenon in the mud extrusion areas, Sidoarjo, East-Java. In: Proceedings of the International Union of Geodesy and Geophysics (IUGG) XXIV, Perugia, Italy.
Allis, R.G., Hunt, T.M., 1986. Analysis of exploration induced gravity changes at Wairakei geothermal field. Geophysics 51, 1647–1660.
Davies, R., Swarbrick, R., Evans, R., Huuse, M., 2007. Birth of a mud volcano: East Java, 29 May 2006. GSA Today 17, 4–9.
Davies, R.J., Brumm, M., Manga, M., Rubiandini, R., Swarbrick, R., Tingay, M., 2008. The East Java mud volcano (2006 to present): an earthquake or drilling trigger? Earth and Planetary Science Letters 272, 627–638.
De Genevraye, P., Samuel, L., 1972. Geology of the Kendeng zone (Central and East Java). In: Proceedings of the Indonesian Petroleum Association, 1st Annual Convention, pp. 17–30.
Deguchi, T., Maruyama, Y., Kato, M., Kobayashi, C., 2007. Surface displacement around mud volcano, East Java captured by Insar using Palsar data. In: Proceedings of the 28th Asian Conference on Remote Sensing.
Fukushima, Y., Mori, J., Hashimoto, M., Kano, Y. Subsidence associated with the LUSI mud eruption, East Java, investigated by SAR interferometry. Journal Marine and Petroleum Geology, Special issue: Mud volcanism, submitted for Q2 publication.
Gafoer, S., Ratman, N., 1999. Peta Geologi Regional Jawa Bagian Timur, Skala 1:500.000. Pusat Penelitian dan Pengembangan Geologi, Bandung.
Hall, R., 2002. Cenozoic geological and plate tectonic evolution of SE Asia and the SW Pacific: computer based reconstructions, model and animations. Journal of Asian Earth Sciences 20, 353–432.
Hamilton, W.R., 1979. Tectonics of the Indonesian Region. US Geological Survey Professional Paper 1078. US Government Printing Office.
Hogan, C.M., 1973. Analysis of highway Noise. Journal of Water, Air, & Soil Pollution 2 (3), 387–392. Biomedical and Life Sciences and Earth and Environmental Science Issue, Springer Verlag, Netherlands.
Issler, D.R., 1992. A new approach to shale compaction and stratigraphic restoration,Beaufort–Mackenzie Basin and Mackenzie Corridor, Northern Canada. AAPG Bulletin 76, 1170–1189.
Istadi, B.P., Kadar, A., Sawolo, N., 2008. Analysis & recent study results on East Java mud volcano. In: Proceedings of the Subsurface Sediment Remobilization and Fluid Flow in Sedimentary Basin Conference, October 2008. The Geological Society, Burlington House, Piccadilly, London.
Jentzsch, G., Calvache, M., Bermudez, A., Ordonez, M., 1999. Microgravity and GPS at Galeras volcano, Columbia: establishment of a new network and first measurements. In: Abstract for the International Symposium on Andean Geodynamics, Go¨ ttingen, October 1999.
Kadar, A.P., Kadar, D., Aziz, F., 2007. Pleistocene stratigraphy of Banjarpanji#1 well and the surrounding area. In: Proceedings of the International Geological Workshop on Sidoarjo Mud Volcano, Jakarta, IAGI-BPPT-LIPI, February 20–21, 2007. Indonesia Agency for the Assessment and Application of Technology, Jakarta.
Karous, M., Hjelt, S.E., 1985. Linear filtering of VLF dip angle measurements. Geophysical Prospecting 31 (5), 782–794.
Kholodov, V.N., 1983. Formation of gas–water fluids in sandy–clayey deposits of Elision Basins. In: Sedimentary Basins and their Oil–Gas Potential. Nauka, Moscow, pp. 28–44 (in Russian).
Kopf, A.J., 2002. Significance of mud volcanism. Reviews of Geophysics 40, 1–52.
Mazzini, A., Svensen, H., Akhmanov, G.G., Aloisi, G., Planke, S., Malthe-Sørenssen, A., Istadi, B., 2007. Triggering and dynamic evolution of the LUSI mud volcano,Indonesia. Earth and Planetary Science Letters 261, 375–388.
Pendrod, K., Ellsworth, A., Farrel, J., 2006. Application of GIS to Estimate the Volume of the Great Johnstown Flood. http://www2.nature.nps.gov/ParkScience/index. cfm?ArticleID¼74 (7 December 2006).
Pramono, G.H., Poniman, A., 2007. The validation of mud flow model using multitemporal satellite images. In: Proceedings of the 28th Asian Conference on Remote Sensing, Kuala Lumpur, Malaysia. Pramono, G.H., Sardjono, S.P., 2007. The application of geographic information systems to measure mud dimension in Porong, Sidoarjo. Forum Geografi 21 (2).
Prasetyadi, C., Suparka, E.R., Harsolumakso, A.H., Sapiie, B., 2006. The occurrence of a newly found Eocene tectonic Melange in Karangsambung Area, Central Java. In: Proceedings of the IAGI 35th Annual Meeting and Exhibition, Pekanbaru, October 2006.
Price, M., October–December 2002. Deriving Volume with ArcGIS Spatial Analyst. ArcUser, pp. 53–56.
Reynolds, J.M., 1997. An Introduction to Applied and Environmental Geophysics. John Wiley & Sons Ltd., England, p. 428.
Sawolo, N., Sutriono, E., Istadi, B.P., Darmoyo, A.B., 2008. East Java mud volcano (LUSI): drilling facts and analysis. In: Proceedings of the African Energy Global Impact, AAPG International Conference and Exhibition, October 2008, Cape Town, South Africa.
Sawolo, N., Sutriono, E., Istadi, B.P., Darmoyo, A.B. Mud volcano triggering controversy: was it caused by drilling? Journal Marine and Petroleum Geology, Special issue: Mud volcanism, submitted for publication. Q3
Schiller, D.M., Seubert, B.W., Musliki, S., Abdullah, M., 1994. The reservoir potential of globigerinid sands in Indonesia. In: Proceedings of the 23rd Annual Convention, Proceedings of the Indonesian Petroleum Association, vol. 1, p. 189ff.
Sclater, J.G., Christie, P.A.F., 1980. Continental stretching: an explanation of the post- Mid-Cretaceous subsidence of the Central North Sea Basin. Journal of Geophysical Research 85 (B7), 3711–3739.
Smyth, H., Hall, R., Hamilton, J., Kinny, P., 2005. East Java: Cenozoic Basins, volcano and S–N cross section line ancient basement. In: Proceedings of the 13th Annual Convention & Exhibition. Indonesian Petroleum Association.
Sribudiyani, Prasetya, I., Muchsin, N., Sapiie, B., Ryacudu, R., Asikin, S., Kunto, T., Harsolumakso, A.H., Astono, P., Yulianto, I., 2003. The collision of the East Java microplate and its implication for hydrocarbon occurrences in the East Java Basin. In: 29th Annual Convention Proceedings of the Indonesian PetroleumAssociation.
Sumintadireja, P., Purwaman, I., Istadi, B., Darmoyo, A.B., 2007. Geology and geophysics study in revealing subsurface condition of Banjarpanji mud extrusion, Sidoarjo, East Java, Indonesia. In: Proceedings of the International Union of Geodesy and Geophysics (IUGG) XXIV, Perugia, Italy.
Talwani, M., Worzel, J.L., Landisman, M., 1959. Rapid gravity computation for twodimensional bodies with application to the Mendocino submarine fracture zone. Journal of Geophysical Research 64 (1), 49–59.
Tingay, M., Heidbach, O., Davies, R., Swarbrick, R., August 2008. Triggering of the Lusi mud eruption: earthquake versus drilling initiation. Geology 36 (8), 639– 642.
van Bemmelen, R.W., 1949. General Geology of Indonesia and Adjacent Archipelagos. In: The Geology of Indonesia, vol. IA. Martinis Nijhoff, The Hague, 732 pp. Willumsen, P., Schiller, D.M., 1994. High quality volcaniclastic sandstone reservoirs in East Java, Indonesia. In: 23rd Annual Convention, vol. I. IPA, pp. 101–111.
Page updated
Google Sites
Report abuse