Identifikasi Sistem Geothermal Menggunakan Metode Audio Magnetotelurik (AMT) Daerah Gedongsongo, Ungaran, Jawa Tengah
oleh: Dwiky Perdana Susanto, Dwi Noviyanto, Mochammad Husni Rizal, Adytia Laksamana Putra, Astya Brilliana, Setyarini Wijayanti Utami*
*)
Universitas Gadjah Mada, Yogyakarta Lab. Geofisika UGM Sekip Utara Unit III PO.BOX BLS 21 Yogyakarta 55281
Intisari Daerah Gedongsongo, Jawa Tengah merupakan salah satu area prospek panas bumi, dimana terdapat manifestasi yang berhubungan dengan aktivitas panas bumi. Manifestasi ini berupa fumarol, steaming fumarol, steaming ground dan hot spring . Berdasarkan penelitian geologi dan geokimia yang telah banyak dilakukan, diperkirakan daerah ini merupakan upflow. upflow. Telah dilakukan survei Audio Magnetotelurik di daerah prospek panas bumi Gedongsongo, yang bertujuan untuk mengidentifikasi beberapa komponen penyusun sistem panas bumi, yaitu kedalaman batuan penudung dan batuan reservoir. Pengukuran AMT telah dilakukan pada 32 titik pengukuran sepanjang 5 lintasan dengan spasi antar titik 500 meter. Berdasarkan hasil analisis data, didapat 3 zona resistivitas utama pada lokasi penelitian yang terdiri dari zona resistivitas tinggi, resistivitas sedang, dan resisitivitas rendah dengan nilai masing-masing zona 1000 – 1905,461 Ohmm , 275,42 – 724,436 Ohmm, 1.56 – 104.7129 Ohmm. Zona dengan nilai resistivitas 1000 – 1000 – 1905,461 1905,461 Ohmm pada Line A dan Line B diduga merupakan zona reservoir panas bumi pada lokasi penelitian. Zona dengan nilai resisitivitas di bawah 1000 diduga merupakan zona batuan pendung dan batuan ubahan pada sistem panas bumi lokasi penelitian. Kata kunci : Gedong Songo, Audio Magnetotelurik, Sistem Panas Bumi
1. Pendahuluan
reservoir dan batuan penudung pada suatu
Seperti sumber daya energi pada umumnya,
sistem panas bumi dari peta persebaran
sumber
bisa
resistivitas batuannya. Daya tembus dari
diperoleh melalui pengukuran karakteristik
metode ini adalah sekitar 1000 meter
fisika dan kimia dari suatu daerah yang
sehingga
dianggap memiliki manifestasi dan prospek
menjangkau kedalaman batuan penudung
panas bumi. Identifikasi terhadap komponen
dan batuan reservoir.
energi
panas
bumi
juga
penyusun sistem panas bumi (contohnya batuan penudung dan batuan reservoir) penting
untuk
menentukan
dianggap
cukup
untuk
dapat
2. Tinjauan Pustaka a. Tinjauan Geologi
prospeknya eksplorasi
Gunung Ungaran merupakan gunung api
lanjut dan eksploitasi. Jika prospek dari
kuarter yang terletak di sebelah utara deretan
suatu sistem panas bumi menjanjikan, maka
Pegunungan Serayu Utara, yaitu gunungapi
cadangan energi yang dimilikinya dapat
Ungaran – Telomoyo – Merbabu – Merapi.
sehingga
dapat
dilaksanakan
dikembangkan menjadi pembangkit tenaga
Geologi
permukaan
Gunung
listrik. Pemanfaatan sumber daya panas
Ungaran didominasi oleh batuan vulkanik
bumi
mendukung
berumur Kuarter. Menurut Budiardjo et. al.
terhadap
(1997), stratigrafi daerah Ungaran dari yang
juga
akan
berkurangnya
turut
ketergantungan
bahan bakar fosil dan mengurangi tingkat polusi. Di daerah Gunung Ungaran terdapat manifestasi yang tersebar di beberapa titik, berupa fumarol, mata air panas, batuan yang teralterasi, dan sisa manifestasi. Berdasarkan tinjauan
geologi,
mengindikasikan
manifestasi daerah
ini
tersebut memiliki
potensi panas bumi.
1. Batugamping volkanik 2. Breksi volkanik III 3. Batupasir volkanik 4. Batulempung volkanik 5. Lava andesitik 6. Andesit porfiritik 7. Breksi volkanik II 8. Breksi volkanik I
Akuisisi data dengan metode AMT (audio magnetotelluric) kali ini
tua ke yang muda adalah sebagai berikut:
bertujuan
untuk memperkirakan kedalaman batuan
9. Andesit porfiritik 10. Lava andesit 11. Aluvium
Gambar 1. Peta geologi regional daerah Ungaran (Budiardjo, et. al., 1997)
zona struktur berarah Utara-Selatan sebagai Sistem panas bumi yang berkembang
pembatas pemunculan dan struktur berarah
di Gunung Ungaran secara geologi berada di
Timurlaut – Baratdaya
zona depresi, berupa kerucut-kerucut muda.
Gedongsongo.
untuk
manifestasi
Prospek panas bumi daerah ini bersistem
Manifestasi suhu umumnya terjadi di
dominasi air panas, yang secara struktural
daerah bagian barat dari sungai kecil pada
dikontrol oleh struktur kaldera atau struktur
bagian selatan dari puncak Ungaran. Pada
runtuhan (collapse structure) Ungaran, yang
bagian timur dari sungai tersebut ditemukan
memanjang dari barat hingga tenggara dari
sisa manifestasi (relic manifestations), yang
arah Ungaran. Batuan vulkanik penyusun
diindikasikan dengan banyaknya pohon mati
pre-caldera dikontrol oleh sistem sesar yang
dan diselubungi oleh sulfur dan lumpur. Bau
berarah Baratlaut – Baratdaya dan Tenggara-
sulfur dapat tercium pada jarak dekat, dan
Barat. Pada batuan vulkanik penyusun post-
terlihat adanya struktur kecil di sepanjang
caldera hanya terdapat sedikit struktur yang
tanah yang juga diselubungi oleh lumpur
dikontrol
dan
(Budiardjo
oleh et.
sistem al.
sesar
1997).
regional
Pemunculan
manifestasi Gedongsongo dikontrol oleh
sulfur.
Kenampakan
ini
diinterpretasikan sebagai mulut mata air panas atau hangat di masa lampau.
Gambar 2. Cross-section arah barat-timur sepanjang manifestasi daerah Gedongsongo
Distribusi manifestasi menunjukkan
Dari pengukuran tersebut, diperoleh model
arah utara-selatan, sejajar dengan lereng
resistivitas (Gambar 3) yang menunjukkan
curam dan sungai kecil yang mungkin saja
fenomena geologi berupa kerucut intrusi
terbentuk akibat struktur runtuhan (collapse
dengan resistivitas sangat tinggi (>10000
structure) atau sesar normal seperti yang
ohm.m) yang diperkirakan merupakan salah
umumnya terjadi di daerah pegunungan.
satu kerucut parasitic di Ungaran. Kerucut
Struktur ini bersifat permeabel dan dianggap
ini diduga merupakan intrusi andesit yang
mengatur transfer panas secara konvektif.
masih menyimpan panas dan masih dalam
Patahan tersilifikasi yang ditemukan pada
tahap
daerah tanah beruap masa lampau (relic
penudung
steaming ground ) di bagian utara juga
resistivitas >1000 ohm.m berada sekitar 400
mengindikasikan adanya sesar normal.
meter
pendinginan. yang
dibawah
Perkiraan
dicirikan
Kawah
batuan
dengan
Item,
nilai
dengan
penyusun utama lava dan breksi lahar andesitik
b. Tinjauan geofisika
Pengukuran metode magnetotelurik telah dilakukan
di
lapangan
panas
bumi
Gedongsongo oleh S.Widarto et. al. (2007).
yang
terkompakkan.
Dibawah
lapisan yang terduga penudung tersebut ditemukan zona sangat konduktif (<10 ohm.m) yang diduga merupakan puncak
reservoir.
Dibawah
puncak
reservoir
andesit dengan resistivitas > 1000Ωm. Dan
tersebut, ditemukan persebaran resistivitas
batuan
antara 30 sampai 300 ohm.m yang terletak
gunung
pada kedalaman 600 hingga 1600 meter dari
sendimen
permukaan dan diduga merupakan zona
resistivitas 100-1000 Ωm dan 10-100 Ωm.
reservoir
Sementara
lapangan
panas
bumi
Gedongsongo.
reservoir api
utama
pre-caldera
laut
tersier
batuan
berupa dan dengan
penudungnya
batuan batuan harga
berupa
batuan post caldera yang mengalami alterasi
Gambar 3. Model tahanan jenis hasil inversi 2D ( S.Widarto et. al, 2007 )
Pengukuran dengan metode CS-AMT di
hidrotermal dengan nilai resistivitas sangat
area Gedongsongo juga telah dilakukan oleh
rendah <10 Ωm. Lapisan cap rock ini
Faulin, 2002 yakni di daerah manifestasi
diindikasikan menyebar kearah barat daya
fumarol dan sebelah selatan gedongsongo
berkaitan dengan struktur robohan yang
hingga banaran. Dari penelitian ini juga
mengontrol
diperoleh
gedongsongo.
ungaran
model yang
reservoir
diperkirakan
panas terdiri
bumi dari
batuan dasar berupa intrusi batuan beku
sistem
panas
bumi
dengan prinsip metode magnetotellurik.
3. Dasar Teori 1. Ko ns ep da sa r mag ne to tel lu ri k
pada frekuensi yang ditangkap. Secara
dan audio magnetotellurik
Dalam
metode
Perbedaan antara kedua metode ini terletak
magnetotellurik,
keberadaan medan elektromagnetik alam digunakan untuk mempelajari resistivitas listrik medium di bawah permukaan bmi,
sederhana
bumi. Konsep dasar metode magnetotellurik cukup
sederhana,
pengamatan diselidiki
di
suatu
titik
yang
akan
lapangan
nilai
permukaannya,
pada
resistivitas dilakukan
bawah
pengukuran
komponen tangensial medan listrik dan medan
magnetik
dari
medan
elektromagnetik yang berasal dari alam. Fluktuasi medan elektromagnetik dari alam terutama berasal dari aktivitas meteorologi dan aliran arus listrik di ionsfer. Medan elektromagnetik
ini
dalam
menghasilkan medan induksi sekunder yang dikontrol oleh sifat-sifat kelistrikan medium. Dalam survei magnetotellurik, receiver akan menerima medan sekunder dari medium. Secara magnetotellurik
umum
metode
menggunakan
sumber
gelombang elektromagnetik alam. Pada dasarnya
prinsip
Magnetotelluric
metode
(AMT)
hamper
Audio sama
AMT
elektromagnetik yang telah menginduksi bawah permukaan bumi. Perbandingan medan
listrik
dan
antara
intensitas
medan
magnetik
menunjukan sifat impedansi listrik. Pada kondisi
tertentu
impedansi
listrik
ini
merupakan fungsi dari sifat listrik medium atau batuan. Dengan menentukan impedansi pada
sederetan
diperoleh
frekuensi
informasi
konduktivitas
atau
maka
tentang
dapat
penampang
resistivitas
medium
sebagai fungsi kedalaman dibawah titik pengukuran
tersebut
dan
dapat
menggambarkan struktur geologi bawah permukaan. Pemakaiangelombang
interaksinya
dengan medium di bawah permukaan akan
metode
dilakukan dengan mengukur gelombang
dengan melakukan pengukuran medan listrik dan medan magnetik di atas permukaan
penggunaan
elektromagnetik dengan frekuensi tinggi hanya
memberikan
daya
tembus
yang
dangkal dikarenakan adanya penyerapan energy sebesar 37%. Kedalaman dimana amplitude
gelombang
elektromagnetik
terinduksi menjadi 1/e dari amplitude awal gelombang tersebut dipermukaan sisebut skin
depth.
Dengan
emikian
untuk
mendapatkan informasi mengenai resistivias atau konduktivitas pada lapisan yang lebih
dalam,
maka
digunakan
gelombang
Energy elektromagnetik pada frekuensi
elektromagnetik dengan frekueansi yang
rendah berosilasi atara permukaan bumi
rendah.
yang konduktif dan lapisan ionosfer di udara.
2. Me dan El ektro mag ne ti k
Medan
elektromagnetik
Hampir selalu terjadi klat setiap saat
mempunyai spectrum dengan frekuensi dari
di salah satu tempat di muka bumi ini. Ada
-4
10 Hz sampai dengan puluhan ribu Hz.
tiga tempat utama disaerah ekuator sebagai
Medan
dapat
pusat terjadinya kilat dengan rata-rata 200
ditimbulkan oleh berbagai sebab. Frekuensi
hari tiap tahun, yaitu Brasil, Afrika Tengah,
medan elektromagnetik <1Hz biasanya
dan Malysia.
elektromagnetik
alam
ditimbulkan oleh interaksi antar partikel yang dipancarkan oleh matahari dengan
3. Pe rs am aa n Ge lom ba ng EM
Metode
magnetotelluric
atmosfer bumi dan magnetosfer. Frekuensi
menggunakan teori elektromagnetik yang
medan elektromagnetik >1Hz dapat berasal
dijabarkan
dari kegiatan meteorology (guntur dan
dalam bentuk dfferensial , yaitu:
tujuan komunikasi) buatan manusia.
alam
memiliki
frekuensi
>1Hz
berasosiaso dengan guntur atau badai. Kilat
energi
yang
uniform.
Medan
elektromagnetik yang berasal dari kilat menjalar ke tempat yang lebih jauh. Dalam
(Hukum
Ampere)
→
(Hukum Fluk
Kontinyu)
yang terjadi pada tempat yang cukup jauh dari titik pengamatan merupakan sumber
(Hukum
→ →
yang berasal dari aktivitas meteorology (jaraknya cukup jauh), terutama kilat yang
→ → →
Dengan volt/meter),
(Hukum Coulomb)
adalah
adalah
medan
intensitas
magnetic (amper/meter), 2
→
medan
magnetic (weber/meter atau Tesla), 2
mengalami penguatan dengan cara osilasi.
listrik
adalh induksi
penjalarannya, frekuensi tingg mengalami pelemahan sedangkan frekuensi rendah
Maxwell
Faraday)
yang
terpenting adalah medan elektromagnetik
persamaan
→
kilat) serta dari sistem pemancar (untuk
Komponen medan elektromagnetik
dalam
adalah rapat arus (amper/meter ),
→ 2
→
adalah
perpindahan listrik (coulomb/meter ) dan q
adalah
rapat
muatan
listrik
3
(coulomb/meter ).
variasi terhadap waktu dinyatakan sebagai fungsi sinusoidal, maka:
Dalam keadaan homogeny isotrop berlaku hubungan :
dengan magnetic
→ → → → → → →
Dengan
4. Ab so rp si Gel om ba ng Bi da ng
Konduktivitas
adalah permeabilitas
permitivitas listrik (farad/meter), konduktivitas (mho/meter) dan
yang
batuan
merupakan
adalah
parameter
adalah
penentuan
adalah
Biasanya material bumi mempunyai harga
sunbstitusi
permitivitas
menentukan
struktur
konduktivitas
resistivitas medium (ohmmeter). melakukan
bawah
siemens
(3.2) kemudian delakukan operasi curl pada
sehingga
persamaan tersebut, maka akan didapatkan
Dengan
magnetic
menggunakan
secara
→ → →
vector
Selain bersifat konduktif, material bumi ada juga yang bersifat resistif. Untuk
maka
akan
didapatkan persamaan gelombang untuk
√
kasus ini maka harga
→ → → → → →
5. Im pe da ns i da n Tah ana n Je ni s Se mu
Solusi persamaan (3.8) dan (3.9) persamaan
gelombang
, sehingga
konstanta perambatan akan menjadi :
medan listrik dan medan magnetic:
merupakan
,
dan konstanta adalah:
terpisah.
identitas
dapat diabaikan. Pada kasus ini
persamaan gelombang untuk medan listrik medan
dan
farad/meter. Untuk
frekuensi kurang dari 100 Hz,
dalam
permukaan.
persamaan (3.5), (3.6) pada persamaan (3.1),
dan
adalah
frekuensi gelombang elektromagnetik (Hz).
(henry/meter),
Dengan
dan
yang
merupakan fungsi dari waktu dan jarak. Jika
Untuk matematis
merumuskan perambatan
elektromagnetik
di
dalam
secara gelombang medium
horizontal berlapis yang tiap lapisannya
homogeny isotrop berdasarkan pendekatan
menimbulkan medan listrik primer
Cagniard, maka dipakai beberapa asumsi
E0x yang tegak lurus terhadap medan
sebagai berikut:
magnetic H0y.
pusat koordinat berada di bawah
Di
dalam
permukaan bumi, dengan sumbu z
besarnya harga impedansi adalah modulus
tegak lurus bidang permukaan dan
dari Zxy :
Berdasarkan persamaan (3.46) dan jika
bumi terdiri dari n buah lapisan,
dimasukkan harga
masing-masing
ruang
ini
dianggap
konduksi yang mengalir kea rah sumbu x maka tidak ada aruus konduksi yang mengalir kea rah sumbu
y
dan
lapisan
ke-n
mempunyai tebal tak berhingga ke bawah,
dan ketebalan lapisannya
dengan
Sedangkan
magnetiknya
harga
adalah
dan
.
permeabilitas dasn
permeabilitas listriknya adalah dengan
hampa,
,
henry/m
farad/m,
sumber gelombang elektromagnetik
maka
permeabilitas
dapat
dituliskan
besarnya resistivitas dan fase, yaitu:
|| * + Jika
tanah
dan
sebagai
perambatan tidak homogeny maka menjadi
harga konduktivitas lapisan ke-m adalah
isotrop,
dengan bidang permukaan),
homogeny isotrop. Jika ada arus
homogeny
||
positif ke bawah (bidang x, y parallel
medium
resistivitas
medium
akan
(resistivitas semu). Definisi
semu
secara
umum
adalah
resistivitas yang terukur di atas permukaan medium berlapis lapis yang mempunyai perbedaan resistivitas dan ketebalan lapisan. Resistivitas
semu
pada
magnetotelluric
bergantung pada parameter sifat medium (ketebalan dan resistivitas batuan) dan frekuensi gelombang elektromagneti yang diukur.
berasal dari sebuah lempeng arus yang terletak diatas permukaan bumi
6. Sk in De pt h
dan menimbulkan medan magnetic
Hal ini menunjukan bahwa amplitude
uniform H0y. medan magnetic ini sebanding dengan waktu sehingga
pada kedalaman
besarnya medan akan
berkurang 1/e dari harganya dipermukaan
bumi, oleh karena itu fenomena penjalaran
3) Atmospheric
noise
disebabkan
oleh
gelombang elektromagnetik didalam bumi
hujan angina disertai petir dan Guntur
lebih mencerminkan kejadian difusi dari
serta efek static dari udara dan awan.
fenomena penjalaran gelombang itu sendiri.
Dihindari dengan melakukan low pass
Dalam SI, skin depth :
filtering serta stacking and averaging
( ) ( ) 7. Su mbe r De ra u (N oi se ) da lam
yang
meadai
ketika
dilakukan
pengolahan data 4) Wind noise disebabkan oleh angina yang dapat merubah posisi koil magnetic maka sebaiknya saat pengukuran koil dipendam didalam tanah yang juga
Me to de AM T
Noise dalam rekaman akuisisi AMT terbagi
berfungsi meminimalkan
efek panas
pada koil akibat pengaruh sinar matahari
empat kategori, yaitu: 1) Instrumentation noise ditimbulkan oleh alat seperti rendahnya input impedansi
langsung
dan
gerakan
koil
induksi
(gerakan kabel penghubung elektroda).
pada receiver , kemampuan filtering yang tidak memadai. Hal ini dapat dihindari dengan cara mengecek respon system secara periodic dan melakukan analisis terhadap
data
pengukuran
sehingga
dapat diketahui data yang baik dan yang
2) Listrik tegangan tinggi, jaringan listrik, kabel yang dibawah permukaan tanah merupakan sumber utama cultural noise. Cara meminimalisirnya adalah dengan
noise tersebut
Kedalaman yang besar berasosiasi dengan energy yang besar oada spektru dengan frekuensi rendah dari gelombang elektromagnetik. penggunaan
jelek
menjauhkan titik amat
8. Re so na ns i da n Ked al am an
Namun
gelombang
demikian
elektromagnetik
berfrekuensi rendah mempunyai resolusi yang rendah pula. Hal ini menyebabkan penggunaannya
tidak
mendeteksi lapisan tipis.
dari sumber Metode Penelitian Akuisisi Data
efektif
untuk
Akuisisi data dilakukan di area Gedongsongo
dan
desa
Candi
dengan
membuat titik survey dengan jumlah 32 titik dan spasi rata – rata 500 m (ditunjukan oleh Gambar 4). Proses akuisisi ini dilakukan pada 9 Juni – 18 Juni 2014.
Gambar 5. Diagram alir penelitian
Hasil dan Pembahasaan
Gambar 4. Peta titik en ukuran
Berdasarkan hasil pengolahan inversi Pengolahan Data
dari software IP2WINMT, dapat dilihat hasil
Pengolahan data dilakukan dengan
korelasi dari titik AMT 01 – AMT 08 (Line
menggunakan perangkat lunak IP2WIN(MT
A), AMT 09 – AMT 16 (Line B), dan AMT
)versi 2.0. untuk proses inversi 1D. Proses
17 – AMT 24 (Line C) yang menunjukkan
gridding data
adanya zona batuan yang memiliki nilai
dengan
di
kedalaman
menggunakan
dilakukan lunak
resistivitas berkisar dari 1,58 – 1905,461
Surfer versi 11. Proses gridding dilakukan
Ohmm (Gambar 6). Nilai ini dibagi lagi
dengan
natural
menjadi beberapa zona nilai log resistivitas
neighbour sehingga didapatkan gambaran
yang dianggap mewakili beberapa lapisan
sebarantahanan
yang berada di bawah permukaan. Secara
menggunakan
jenis
perangkat
metode
sebenarnya
(true
resistivity) sebagai fungsi kedalaman.
garis besar, kondisi bawah permukaan terdiri
Berikut diagram alir pengolahan data dan
dari 3 zona utama, zona pertama merupakan
interpretasi pada penelitian kali ini.
zona batuan resistivas yang digolongkan rendah dibandingkan sekitarnya dengan nilai
resisitivitas dari 1.56 – 104.7129 Ohmm.
yang juga dapat diduga menjadi reservoir
Zona ini berada pada zona selatan area
pada daerah penelitian. Sedangkan zona
penelitian pada Line A dan Line B, dan
dengan nilai resistivitas di bawahnya diduga
memanjang dari utara – selatan pada Line C.
merupakan
Zona kedua merupakan zona batuan yang
batuan ubahan pada area panas bumi lokasi
memiliki nilai resisitivitas 275,42 – 724,436
penelitian.
zona
batuan
penudung
dan
Ohmm. Zona ini terletak pada bagian utara di Line A dan Line B, serta pada zona
Kesimpulan dan Saran
sebelah selatan dari Line C. Zona ketiga
1. Terdapat 3 zona resistivitas utama pada
merupakan zona dengan nilai resistivitas
lokasi penelitian yang terdiri dari zona
1000 – 1905,461 Ohmm terletak pada
resistivitas tinggi, resistivitas sedang, dan
bagian utara di Line A dan Line B, serta
resisitivitas rendah dengan nilai masing-
pada bagian selatan di Line C.
masing zona 1000 – 1905,461 Ohmm ,
Berdasarkan
hasil
korelasi
dan
275,42 – 724,436 Ohmm, 1.56 – 104.7129
interpretasi, maka diduga bahwa batuan
Ohmm.
pada Line A dan Line Bsaling berkorelasi
2.
satu sama lain yang dapat dilihat pada
1905,461 Ohmm pada Line A dan Line B
penyebaranya yang cenderung sama-sama
diduga merupakan zona reservoir panas
mengarah ke arah utara untuk batuan dengan
bumi pada lokasi penelitian.
Zona dengan nilai resistivitas 1000 –
resistivtas sedang-tinggi dan mengarah ke
3. Zona dengan nilai resisitivitas di bawah
selatan
1000
untuk
zona
resistivitas
yang
diduga
merupakan
zona
batuan
dianggap rendah. Sementara pada Line C,
pendung dan batuan ubahan pada sistem
zona batuan dengan nilai resistivitas rendah
panas bumi lokasi penelitian.
terletak
memanjang
dari
utara-selatan
dengan posisi zona resistivitas tinggi berada di sebelah selatan (Gambar 7). Hasil interpretasi pada ketiga zona batuan ini menunjukkan bahwa zona batuan dengan nilai resistivitas 1000 – 1905,461 Ohmm pada Line A dan Line B dapat diinterpretasikan sebagai zona batuan beku
Daftar Pustaka
Budihardjo, B., Nugroho, Budiharti, M., 1997 , Resources Characteristics of the Ungaran Field, Central Java, Indonesia, Proceeding Seminar Nasional Sumber Daya Geologist Indonesia, Fakultas Teknik Geologi dan Mineral, UPN “Veteran”, Yogyakarta. Faulin, Taubah, 2002 , Interpretasi Resistivitas Batuan Daerah Prospek Panasbumi Gedongsongo, FMIPA UGM: Yogyakarta. Gaffar Z. Eddy, Dadan D. Wardhana, Djedi S. Widarto. 2007. Studi Geofisika Terpadu di Lereng Selatan Ungaran, Jawa Tengah, dan Implikasinya Terhadap Struktur PanasBumi. Jurnal Meteorologi dan Geofisika : 101-119. Nukman, M. 2014. Overview of Gedongsongo Manifestations of the Ungaran Geothermal Prospect, Central Java, Indonesia : a preliminary account. Geofisika, Universitas Gadjah Mada : Yogyakarta.
Lampiran
Gambar 6.Penampang Resistivitas Bawah Permukaan
Gambar 7.Korelasi Penampang Resistivitas
Gambar 8.Contoh Hasil Inversi IP2WINMT Titik AMT 01 Line A
Gambar 9.Contoh Hasil Inversi IP2WINMT Titik AMT 09 Line B
Gambar 10.Contoh Hasil Inversi IP2WINMT Titik AMT 17 Line C