Youngster Physics Journal ISSN : 2302 - 7371

Vol. 07, No. 1, Januari 2018, Hal. 11-18

Identifikasi struktur sesar daerah manifestasi panas bumi “X” di Kabupaten

Manggarai Nusa Tenggara Timur berdasarkan analisis Horizontal Gradient
Anik Fitrochaton Chasanah1), Agus Setyawan1), Dadan Dani Wardhana2)
1)
Departemen Fisika, Fakultas Sains dan Matematika, Universitas Diponegoro, Semarang
2)
Pusat Penelitian Geoteknologi LIPI Bandung, Indonesia
E-mail:[email protected]

ABSTRACT
Area “X” located in Manggarai of East Nusa Tenggara is volcanic area in Indonesia which has the potential
for geothermal to be indicated by the appearance of a geothermal manifestation of hot springs and fumarole or sulfatar.
Processing data by analysis gradient horizontal has been done to determine the structure of fault. From the map
complete Bouguer anomali (ABL) shows anomalies in the middle of the study area which is a crater form of Mount
Poco Ranakah, Poco Manggung, Poco Mandasawu and Poco Kasteno. In the southwest of the study area there is a low
anomaly associated with alteration zone of rocks and hot springs. The high anomaly to the North of the study area lies
in the crumpled hills, so that the area is dominated by high anomaly values. Gradient analysis were performed
determine the limits of lithogy and the structure of fault associated with the appearance of hot springs and fumarol or
sulfate in the study area. Geothermal manifestastios MAP 6, MAP 7, MAP 8, F/S 1 and F/S 2 controlled by strike slip
fault of horizontal gradient. Geothermal manifestations of MAP 1, MAP 4 and MAP 11 are located on geological fault
trending Northeast-Southwest and Northwest-Southeast.
Keywords: gravity data, gradient analysis

ABSTRAK
Daerah “X” yang terletak di Kabupaten Manggarai Nusa Tenggara Timur merupakan daerah vulkanik di
Indonesia yang memiliki potensi panas bumi diindikasikan oleh kemunculan manifestasi panas bumi berupa mata air
panas dan fumarol atau sulfatar. Pengolahan data gayaberat meliputi analisis horizontal gradien dilakukan struktur
berupa sesar. Dari peta anomali Bouguer lengkap (ABL) memperlihatkan adanya anomali rendah di sebelah tengah
daerah penelitian yang berasosiasi dengan bentuk kawah dari Gunung Poco Ranakah, Poco Manggung, Poco
Mandasawu dan Poco Kasteno. Disebelah Barat Daya daerah penelitian terdapat anomali rendah yang berasosiasi
dengan zona alterasi batuan dan kemunculan beberapa mata air panas. Anomali tinggi di sebelah Utara daerah
penelitian terletak pada daerah perbukitan yang mengalami pengkikisan, sehingga daerah tersebut didominasi oleh
nilai anomali tinggi. Analisis gradien dilakukan untuk mengetahui batas litologi dan jenis sesar yang berhubungan
dengan munculnya mata air panas dan fumarol atau sulfatar di daerah penelitian. Kemunculan manifestasi.
Kemunculan manifestasi MAP 6, MAP 7, MAP 8, F/S 1 dan F/S 2 terkontrol oleh sesar mendatar pada horizontal
gradient. Kemunculan MAP 1, MAP 4 dan MAP 11 berada pada sesar dari informasi geologi berarah Timur Laut –
Barat Daya dan Barat Laut – Tenggara.
Kata kunci: Data gaya berat, analisis gradien

PENDAHULUAN Energi panas bumi adalah energi yang

dapat diperbaharui ramah lingkungan, dan
Pemenuhan kebutuhan energi nasional bersih karena emisi CO2 yang dihasilkan
saat ini masih mengandalkan energi yang berasal rendah. Energi panas bumi ini juga cenderung
dari energi fosil seperti bahan bakar minyak dan tidak akan habis karena proses
gas. Ketersediaan akan sumber energi fosil pembentukannya yang terus menerus selama
semakin menipis mengingat sumber energi ini kondisi lingkungannya (geologi dan hidrologi)
tidak dapat diperbaharui. Sehingga dibutuhkan dapat terjaga keseimbangannya [1]. Pemilihan
suatu energi yang mampu menjadi sumber daya energi panas bumi sebagai energi prioritas ke
prioritas untuk mencukupi ketersedian energi depan mempunyai peranan penting terutama
tersebut. Salah satunya adalah energi panas bumi. sebagai energi domestik mengingat sifat energi

11
Anik Fitrochaton Chasanah, dkk Identifikasi Struktur Sesar dan Litologi.....

ini yang tidak dapat diekspor, dengan demikian DASAR TEORI

energi panas bumi akan menjadi energi prioritas
andalan dan vital karena dapat mengurangi Geologi Daerah Penelitian
ketergantungan Indonesia terhadap sumber energi
Kepulauan Nusa Tenggara terbentuk
fosil dan dapat memberikan nilai tambah dalam
sebagai akibat dari subduksi lempeng Indo-
rangka optimalisasi pemanfaatan aneka ragam
Australia di bawah busur Sunda - Banda
sumber energi di Indonesia [2].
selama Tersier atas. Subduksi ini membentuk
Daerah “X” yang terletak di kabupaten
Volcanic Arc atau Busur Vulkanik di
Manggarai kepulauan Flores Nusa Tenggara
kepulauan Nusa Tenggara [9]. Terdapat
Timur merupakan salah satu daerah vulkanik di
beberapa gunung api yang masih aktif yaitu
Indonesia yang memiliki potensi panas bumi
gunung api Poco Mandasawu dan gunung
diindikasikan oleh kemunculan manifestasi panas
Ranakah yang merupakan anak gunung api
bumi berupa mata air panas dan fumarol atau
yang dihasilkan oleh letusan gunung api Poco
sulfatar [3]. Beberapa penelitian gaya berat telah
Ranakah yang berada di sayap bagian timur
dilakukan di daerah manifestasi panas bumi “X”
laut [10].
kabupaten Manggarai oleh Dadan D.Wardhana
Struktur geologi yang terdapat di
dan Fuad yang menyatakan terdapat sesar berarah
daerah penelitian berupa sesar, lipatan dan
barat laut-tenggara [4]. Penelitian lain juga telah
kelurusan. Sesar yang terdapat pada batuan
dilakukan oleh T.Padmawidjaja (2010) yang
Miosen Tengah dan Pliosen berarah Barat Laut
menyatakan bahwa terdapat sesar berarah barat
– Tenggara dan Timur Laut – Barat Daya.
laut-tenggara dan dicirikan oleh adanya kelurusan
Kelurusan yang terdapat pada batuan Miosen
belokan kontur, kelurusan pantai, dan adanya
Tengah sampai termuda yaitu
offset kontur. Diduga sesar ini sesar mendatar
batuangunungapi Holosen berarah Barat Laut
menganan dengan kemiringan ke arah barat daya
– Tenggara dan Barat Daya – Timur Laut [3].
[5].
Berdasarkan Peta Geologi Lembar
Dalam eksplorasi sumber energi panas
Ruteng skala 1:125.000 oleh Koesoemadinata
bumi, metode geofisika memiliki peranan yang
dkk (1994) dalam Nasution (2014) terlihat
penting yaitu untuk menentukan struktur geologi
pada gambar 1 bahwa daerah penelitian
dan komposisi bawah permukaan dengan
tersusun oleh beberapa satuan batuan dari
parameter fisika [6]. Metode geofisika yang
muda ke tua yang terdiri dari: satuan vulkanik
digunakan adalah metode gaya berat. Metode
Poco Ranakah komplek, satuan vulkanik Poco
gaya berat adalah metode yang didasarkan pada
Mandasawu-Kasteno komplek, satuan
pengukuran variasi medan gravitasi dipermukaan
vulkanik Poco Leok 3, satuan vulkanik Poco
bumi akibat perbedaan rapat massa batuan di
Leok 2, satuan vulkanik Poco Rii, satuan
bawah permukaan bumi [7]. Metode gaya berat
vulkanik Poco Leok 1, satuan vulkanik Poco
pada sistem panas bumi digunakan untuk
Cuing komplek, satuan sedimen tersier[11].
mengetahui keberadaan struktur yang mengontrol
keberadaan manifestasi panas bumi yang muncul Teori gayaberat Newton
di permukaan.
Analisis gradien merupakan salah satu Teori yang mendasari metode gaya
jenis filtering untuk menjelaskan kontras densitas berat adalah hukum Newton tentang gaya
dari hasil anomali data gaya berat. Analisis gravitasi. Hukum gravitasi Newton tersebut
gradien dengan horizontal gradient mampu menyatakan bahwa gaya tarik menarik antara
memetakan keberadaan struktur berupa sesar dan dua buah benda sebanding dengan massa
batas kontak litologi batuan pada nilai kedua benda tersebut dan berbanding terbalik
maksimumnya [8].

12
Youngster Physics Journal ISSN : 2302 - 7371
Vol. 07, No. 1, Januari 2018, Hal. 11-18

dengan jarak kuadrat antara pusat massa kedua minus pada persamaan (1) penting karena arah
benda [12]. r̂ berlawanan dengan arah gaya tarik menarik
[13].

Percepatan Gaya berat

Besaran yang diukur dalam metode

gaya berat adalah percepatan gaya berat ⃗g.
Hukum Newton II menyatakan bahwa sebuah
gaya adalah hasil perkalian antara massa dan
percepatan. Hukum Newton II tersebut
dijelaskan oleh persamaan (2).

⃗ (𝑟) = 𝑚g
F ⃗ (2)
dengan F⃗ adalah gaya pada sebuah massa
(Newton), m adalah massa benda (kg), dan 𝑔
adalah percepatan gayaberat (m/s2).
Interaksi antara bumi (bermassa M)
dengan benda di permukaan bumi (bermassa
m) sejauh jarak r dari pusat keduanya juga
memenuhi hukum tersebut, maka dari
persamaan (1) dan (2) didapatkan persamaan
(3).
𝑀
⃗g(𝑟) = −𝐺 𝑟 2 r̂ (3)
2
dengan satuan ⃗g dalam m/s (SI), atau Gal
Gambar 1. Peta geologi daerah penelitian [3]. (Galileo) yaitu 1 cm/s2. Karena pengukuran
dilakukan dalam variasi percepatan gaya berat
Hukum gravitasi Newton diilustrasikan yang begitu kecil, maka satuan yang sering
pada Gambar 2 dan persamaan (1). digunakan adalah miliGal (mGal) [14].
Persamaan (3) menunjukkan bahwa besarnya
percepatan yang disebabkan oleh gaya berat di
bumi (g ⃗ ) adalah berbanding lurus dengan
massa bumi (M) dan berbanding terbalik
dengan kuadrat jari-jari bumi (r).
Gambar 2. Gaya tarik menarik antara dua benda Analisis Gradien
bermassa m dan M [6].
Analisis gradient yang digunakan yaitu
⃗ (𝑟) = −𝐺 𝑀𝑚
F r̂ (1) Horizontal gradient berupa First Horizontal
𝑟2
Gradient dan Second Horizontal Gradient.
dengan 𝐹⃑ adalah gaya tarik menarik (Newton), G Horizontal gradient digunakan untuk
adalah konstanta universal gaya berat menekankan anomali tinggi yang terkandung
(6,67×10−11 m3kg−1s−2), m adalah massa benda di dalam data gaya berat, karena nilai
(kg), M adalah massa bumi (kg), r adalah jarak maksimum yang menunjukkan kepadatan
antar massa m dan massa M (m), dan r̂ adalah lateral pada kontras yang diindikasikan
satuan vektor yang mengarah dari M ke m. Tanda

13
Anik Fitrochaton Chasanah, dkk Identifikasi Struktur Sesar dan Litologi.....

sebagai sesar. Horizontal gradient dibagi menjadi HASIL DAN PEMBAHASAN

dua tahap yaitu FHG (First Horizontal Gradient)
dan SHG (Second Horizontal Gradient). Ketinggian dan Medan Gayaberat
Besarnya horizontal gradien (x,y) didefinisikan Obseravasi
oleh persamaan (4) dan (5).
Penelitian ini dilakukan di daerah “X”
𝜕∆𝑔 𝜕∆𝑔 2 Kabupaten Manggarai, Nusatenggara Timur
𝐻𝐺𝑥 = = √( 𝜕𝑥 ) (4)
𝜕𝑥 dan merupakan daerah yang didominasi oleh
𝜕∆𝑔 𝜕∆𝑔 2 perbukitan. Gambar 4 menunjukkan rentang
𝐻𝐺𝑦 = = √( 𝜕𝑦 ) (5)
𝜕𝑦 ketinggian daerah penelitian antara 355,8
maka besarnya horizontal gradien (x,y) pada hingga 2247,1 meter. Ketinggian rendah
persamaan (6) [15]. berada pada rentang nilai 355,8 hingga 766,0
meter terletak di sebelah Tenggara dan Barat
𝜕∆𝑔 2 𝜕∆𝑔 2 Daya daerah penelitian yang merupakan
𝐻𝐺 (𝑥, 𝑦) = √( 𝜕𝑥 ) + ( 𝜕𝑦 ) (6)
daerah dataran rendah. Ketinggian tinggi
dengan ∆g adalah nilai Bouguer gaya berat total, berada pada rentang nilai 1384,5 hingga
𝜕x adalah nilai gradien sumbu x dan 𝜕𝑦 adalah 2247,1 meter terletak di sebelah tengah daerah
nilai gradien sumbu y. penelitian yang merupakan kenampakan dari
puncak gunung api. Medan gaya berat
METODOLOGI PENELITIAN observasi merupakan nilai gaya berat
sebenarnya pada daerah penelitian.
Data yang digunakan adalah data Gambar 5 menunjukkan rentang nilai
sekunder dari Pusat Penelitian Geoteknologi LIPI medan gaya berat observasi antara 977731,8
Bandung sebanyak 104 titik. Data yang telah hingga 978156,4 mGal. Nilai medan gaya
didapat dalam proses pengambilan data kemudian berat observasi rendah berada di sebelah
dilakukan pengolahan data. Pengolahan data tengah daerah penelitian yang merupakan
dilakukan dengan mengoreksi pasang surut, dan kenampakan dari puncak gunung api dan nilai
apungan pada data gayaberat sehingga didapat medan gaya berat observasi tinggi terletak di
nilai medangayaberat observasi. Nilai medan sebelah tenggara dan barat daerah penelitian
gayaberat observasi kemudian dikoreksi lintang, yang merupakan daerah dengan dataran rendah
koreksi udara bebas, koreksi Bouguer, dan yang berbatasan dengan sungai. Dari gambar 4
koreksi medan sehingga didapat nilai anomali dan gambar 5 menunjukkan bahwa pada
Bouguer lengkap (ABL). Anomali Bouguer daerah yang memiliki ketinggian yang tinggi,
lengkap kemudian dipisahkan menggunakan nilai medan gaya berat observasinya rendah
upward continuation sehingga didapat kontur begitupun sebaliknya. Hal ini sesuai dengan
anomali regional dan residual. Anomali Bouguer teori Hukum Newton tentang gaya tarik-
juga digunakan untuk analisis gradien secara menarik bahwa besar nilai medan gaya berat
horizontal (first horizontal dan second horizontal observasi berbanding terbalik dengan kuadrat
gradient) sehingga didapatkan kontur hasil jarak antar kedua massa. Dalam hal ini jarak
analisis gradient, kemudian diinterpretasi bersama antara kedua massa adalah ketinggian. Dari
dengan informasi geologi untuk mengetahui Hukum Newton tersebut bisa disimpulkan
struktur sesar daerah panasbumi “X” di bahwa semakin tinggi ketinggian suatu daerah
Kabupaten Manggarai, Nusa Tenggara Timur. maka nilai medan gaya berat observasinya
Diagram alir penelitian ditunjukkan pada gambar semakin rendah dan begitupun sebaliknya.
3.

14
Youngster Physics Journal ISSN : 2302 - 7371
Vol. 07, No. 1, Januari 2018, Hal. 11-18

Gambar 5. Medan gaya berat observasi

Anomali Bouguer Lengkap

Interpretasi anomali Bouguer

dilakukan berdasarkan kontur anomali
Bouguer lengkap pada Gambar 6. Gambar 6
menunjukkan rentang nilai anomali Bouguer
antara 50,8 hingga 122,5 mGal. Nilai anomali
rendah ditunjukkan pada rentang nilai 50,8
hingga 91,0 mGal terletak di sebelah barat
Gambar 3. Diagram alir pengolahan data
daya dan tengah daerah penelitian. Anomali
rendah di sebelah tengah berada pada topografi
tinggi yang merupakan daerah perbukitan yang
merupakan puncak dari gunung Poco Ranakah,
Poco Manggung, Poco Mandasawu dan Poco
Kasteno dan anomali rendah sebelah barat
daya terletak pada daerah dataran rendah yang
memperlihatkan adanya kemunculan beberapa
manifestasi panasbumi berupa mata air panas.
Anomali sedang ditunjukkan pada rentang
nilai 91,0 hingga 103,4 mGal yang tersebar di
sebelah Tenggara dan Selatan serta
memanjang di bagian tengah daerah penelitian.
Anomali tinggi ditunjukkan pada rentang nilai
103,4 hingga 122,5 mGal berada di sebelah
Utara dan beberapa titik di sebelah tengah
daerah penelitian yang terletak pada topografi
Gambar 4. Peta ketinggian daerah penelitian sedang. Terdapat zona lemah yang
memisahkan antara nilai anomali Bouguer
tinggi dan nilai anomali Bouguer rendah
disebelah tengah daerah penelitian yang
diinterpretasikan sebagai struktur sesar.

15
Anik Fitrochaton Chasanah, dkk Identifikasi Struktur Sesar dan Litologi.....

Kemunculan manifestasi pada daerah sumbu y. Proses perhitungan menggunakan

panasbumi “X” di kabupaten Manggarai berupa software Geosoft Oasis Montaj menghasilkan
mata air panas dan fumarol atau sulfatar. kontur first horizontal gradient pada gambar 7
Keberadaan manifestasi MAP 1, MAP 2, MAP 3, dan second horizontal gradient pada gambar 8.
MAP 4, MAP 5 dan MAP 6 berada pada nilai Gambar 7 menunjukkan sebaran nilai first
anomali Bouguer rendah. Keberadaan MAP 7, horizontal gradient daerah penelitian. Daerah
MAP 11, F/S 1 dan F/S 2 berada pada nilai penelitian memiliki sebaran nilai first
anomali Bouguer sedang yang di indikasikan horizontal gradient antara 0,001 hingga
sebagai zona lemah dan diinterpretasikan sebagai 0,0374 mGal/m. Gambar 8 menunjukkan
sesar. Keberadaan manifestasi MAP 8, MAP 9, sebaran nilai second horizontal gradient
MAP 10 berada pada nilai anomali Bouguer daerah penelitian. Daerah penelitian memiliki
cenderung tinggi, dimungkinkan kemunculan sebaran nilai second horizontal gradient antara
MAP 10 dikontrol oleh heat source dari sistem 0,0000001 hingga 0,0002408 mGal/m2.
geothermal didaerah tersebut. Hal tersebut Dari peta kontur horizontal gradient
menunjukkan bahwa manifestasi panasbumi “X” daerah penelitian menunjukkan beberapa
di kabupaten Manggarai muncul pada anomali struktur sesar yang menyebar dari tengah
rendah hingga sedang atau dikontrol oleh hingga selatan daerah penelitian. Keberadaan
keberadaan struktur. manifestasi MAP 6, MAP 7, MAP 8, F/S 1 dan
F/S 2 terkorelasi dengan nilai maksimum
horizontal gradien yang di interpretasikan
sebagai sesar mendatar berarah barat laut –
tenggara dan timur laut – barat daya.

Gambar 6. Peta anomali Bouguer lengkap

Analisis Gradien

Horizontal gradient digunakan untuk Gambar 7. Peta first horizontal gradient

menekankan informasi yang terkandung dalam
data gaya berat, karena nilai horizontal gradient
maksimal didentifikasikan sebagai patahan dan
kontak batuan [15]. Interpretasi data gaya berat
menggunakan horizontal gradient didasarkan
pada penurunan persamaan Laplace pada nilai
anomali Bouguer lengkap terhadap sumbu x dan

16
Youngster Physics Journal ISSN : 2302 - 7371
Vol. 07, No. 1, Januari 2018, Hal. 11-18

Panasbumi Flores Nusa Tenggara Timur,

LIPI, Bandung.

[5]. T. Patmawidjaja, 2010, Kondisi Geologi

Daerah Ruteng Ditafsir Pada Data
Gayaberat , JSDG Vol. 20 No. 5 Oktober
2010.

[6]. Rosid, S., 2005, Lecture Notes : Gravity

Methodin Exploration Geophysics: Depok,
University of Indonesia.

[7]. Telford, W. M., Geldart, L. P., & Sheriff,

R. E., 1990, Applied Geophysics 2nd
Edition, Cambridge University Press,
Cambridge.
Gambar 8. Peta second horizontal gradient
[8]. Setyawan, A., Yudianto, H., Nishijima, J.,
Hakim, S., 2015, Horizontal Gradient
KESIMPULAN Analysis for Gravity and Magnetic Data
Beneath Gedongsongo Geothermal
Telah dilakukan analisis gradien yang Manifestations, Ungaran, Indonesia,
meliputi horizontal gradient. Hasil analisis Proceedings World Geothermal Congress
horizontal gradient menginterpretasikan adanya 2015, Melbourne, Australia, 19-25 April
sesar mendatar yang mengontrol munculnya 2015.
manifestasi MAP 6, MAP 7, MAP 8, F/S 1 dan
F/S 2 didaerah penelitian. Keberadaan [9]. Darman, H. dan Sidi, F.H., 2000, An
manifestasi MAP 1, MAP 4 dan MAP 11 Outline of The Geology of Indonesia,
terkontrol oleh sesar geologi berarah timur laut - IAGI, Jakarta.
barat daya dan barat laut - tenggara.
[10]. Gurniwa, Agus., 2005, Penyelidikan
Geokimia Bersistem Lebar Ruteng Barat
DAFTAR PUSTAKA Kabupaten Manggarai dan Manggarai
Barat Provinsi Nusa Tenggara Timur,
[1]. Daud, Yunus, 2010, Diktat Kuliah: Kolokium Hasil Lapangan-DIM.
Introduction to Geothermal System and
Technology, Laboratorium Geofisika, FMIPA [11]. Nasution., 2014, Laporan Kemajuan
Universitas Indonesia. Kegiatan Lapangan Panasbumi Flores
Nusa Tenggara Timur, LIPI, Bandung.
[2]. Kasbani, 2009, Tipe Sistem Panasbumi di
Indonesia dan Estimasi Potensi Energinya, [12]. Jacobs, J. A., Russel, R. D., & Wilson, J.,
Kolokium Sumber Daya Geologi Vol. 4 No. 3. 1974, Physics and Geology,Mc Graw-Hill
Book Company, New York.
[3]. Koesoemadinata, S. Noya, Y. dan
Kadarisman, D., 1994, Peta Geologi Lembar [13]. Wijayanti R,D., Yulianto, Tony., Irham,
Ruteng, Nusa Tenggara, Skala 1:250.000, N., 2009, Prediksi Batas Cekungan Jawa
Pusat Penelitan dan Pengembangan Geologi, Timur Utara Dengan Pegunungan Selatan
Bandung. Di Daerah Magetan-Plaosan Dengan
Metode Gravitasi, Jurnal Sains dan
[4]. Wardhana, Dadan D., Firmansyah, Fuad., Matematika, Nomor 1, Volome 17.
2015, Laporan Kemajuan Kegiatan Lapangan

17
Anik Fitrochaton Chasanah, dkk Identifikasi Struktur Sesar dan Litologi.....

[14]. Nettleton, L.L., 1976, Gravity and Magnetics

in Oil Prospecting. New York: Mc Graw-Hill.

[15]. Graunch VSJ, Cordell L, 1987, Limitations of

determining density or magnetic bounderies
from the horizontal gradient of gravity or
pseudogravity data. Short note, Geophysics 52
(1), 118-121.

18