14
BAB III TEORI DASAR
3.1.
Metoda GPR (Ground Penetrating Radar)
3.1.1. Pengertian Ground Penetrating Radar (GPR) adalah salah satu metode survey untuk soil, bangunan dan kondisi bawah permukaan (dalam interval beberapa centimeter hingga kedalaman 60 meter). Metode GPR ini menggunakan analisa refleksi/pantulan dari gelombang elektromagnetik yang dihasilkan akibat dari perbedaan sifat /konstanta dielektrik benda-benda di bawah permukaan.
Secara umum peralatan GPR terdiri dari dua komponen utama yaitu peralatan pemancar gelombang radar (transmitter) dan peralatan penerima pantulan/ refleksi gelombang radar (tranceiver). Sistem yang digunakan adalah merupakan sistem aktif dimana dilakukan ‘penembakan’ pulsa-pulsa gelombang elektromagnetik (pada interval gelombang radar) untuk kemudian dilakukan perekaman intensitas gelombang radar yang berhasil dipantulkan kembali ke permukaan (Quan dan Haris, 1997).
15
3.1.2. Prinsip dasar ground penetrating radar (GPR) Ground Penetraling Radar (GPR) terdiri dari control unit, transmitter, receiver, note book, kabel serat optik dan tambahan alat lain untuk trigger (Gambar 3.1 dan Gambar 3.2). Dengan frekuensi antenna yang bermacam-macam, seperti ; 25 MHz, 50 MHz, 100 Mhz, 200 MHz, 500 Mhz dan 1000 MHz (Lane, dkk., 1996).
Gambar 3.1. Komponen/peralatan GPR/RAMAC (Lane, dkk., 1996).
Gambar 3.2. Diagram kerja GPR (Arisona, 2009)
16
Control unit berfungsi sebagai pengatur pengumpulan data. Komputer memberikan informasi lengkap bagaimana prosedur yang harus dilakukan, dan saat sistem diaktifkan, control unit mengatur transmitter dan receiver. Control unit menyimpan data mentah dalam sebuah buffer sementara dan saat dibutuhkan, dapat diambil dan ditransfer ke komputer.
Transsmitter menghasilkan energi elektromagnetik dan mengirimnya pada daerah sekitar, khususnya ke dalam medium yang diobservasi. Energi dalam bentuk pulsa pada amplitudo tinggi (370 V) yang dipindahkan ke bagian antena.
Receiver mengkonversi sinyal yang diterima oleh antena menjadi nilai integer. Dalam unit receiver terdapat dua konektor optik, pertama digunakan untuk mentransfer sinyal terkontrol dari control unit (bertanda R) dan lainnya mengirim data yang diperoleh ke control unit (bertanda D). Antenna receiver menerima pulsa yang tidak terabsorbsi oleh bumi tetapi dipantulkan dalam domain waktu tertentu.
Mode konfigurasi antenna transmitter dan receiver pada GPR terdiri dari mode monostatik dan bistatik. Mode monostatik yaitu bila transmitter dan receiver digabung dalam satu antena sedangkan mode bistatik bila kedua antena memiliki jarak pemisah.
Prinsip kerja GPR adalah Transmitter membangkitkan pulsa gelombang elektromagnetik pada frekuensi tertentu sesuai dengan karakteristik antenna
17
tersebut (10 Mhz-4Ghz). Receiver diset untuk melakukan scan yang secara normal mencapai 32 – 512 scan per detik. Setiap hasil scan ditampilkan pada layar monitor sebagai fungsi waktu two-way time travel time, yaitu waktu tempuh gelombang elektromagnetik menjalar dari tranmitter – target – receiver. Tampilan ini disebut dengan radargram (Lane, dkk., 1996).
3.1.3. Persamaan yang mendasari pengukuran GPR Persamaan yang menghubungkan sifat fisik medium dengan medan yang timbul pada medium tersebut dapat dinyatakan dengan (Vasco, 1997):
⃑=
⃑=
⃑= dimana;
⃑
(1)
⃑ =
(3)
⃑
(2)
⃑ adalah pergeseran vector dielektrik, ⃑ adalah rapat fluks magnet, ⃑
merupakan kepadatan saat vektor, c yaitu kecepatan cahaya, dan t adalah waktu masing-masing vektor.
Untuk menyederhanakan masalah, sifat fisik medium diasumsikan tidak bervariasi terhadap waktu dan posisi (homogen isotropi). Untuk pendekatan tersebut, persamaan Maxwell dapat ditulis sebagai berikut (Yee, 1966): ∇⃑ × ⃑ = −
⃑/
(4)
18
∇⃑ × ⃑ =
∇⃑. ⃑ = 0
⃑+
⃑/
∇⃑. ⃑ = 0 dimana
(5) (6) (7)
⃑ merupakan intensitas medan magnet dengan satuan ampere/m,
⃑ adalah perpindahan listrik yang memiliki satuan coulomb/m2, ε merupakan
permitivitas listrik dengan satuan farad/m, dan σ yaitu konduktivitas dengan satuan 1/ohm-m.
Persamaan Maxwel ini adalah landasan teori dari perambatan gelombang elektromagnet. Pada material dielektrik murni, suseptibilitas magnetik (μ), dan permitivitas listrik (ε) adalah konstan dan tidak terdapat atenuasi dalam perambatan gelombang. Tidak sama halnya jika berhadapan dengan material dielektrik (Reynolds, 1996).
Sifat-sifat dari material bumi bergantung dari komposisi dan kandungan air material tersebut. Kedua hal ini mempengaruhi cepat rambat perambatan gelombang dan atenuasi gelombang elektromagnet.
Keberhasilan dari metoda GPR bergantung pada variasi bawah permukaan yang dapat menyebabkan gelombang tertransmisikan. Perbandingan energi yang direfleksikan disebut koefisien refleksi (R) yang ditentukan oleh perbedaan cepat rambat gelombang elektromagnet dan lebih mendasar lagi adalah perbedaan dari konstanta dielektrik relatif dari media yang berdekatan. Hal ini dapat terlihat pada persamaan berikut (Kunz dan Luebbers, 1993):
19
= =
(
(
)
(8)
ℰ
(9)
) ℰ ℰ
ℰ
dimana V1 merupakan cepat rambat geombang elektromagnet pada lapisan 1, V2 adalah cepat rambat geombang elektromagnet pada lapisan 2 , dan V1 < V2, sedangkan ε1 dan ε2 yaitu konstanta dielektrik relatif lapisan 1 dan lapisan 2.
Dalam semua kasus, besarnya R terletak antara -1 dan 1. Bagian dari energi yang ditransmisikan sama dengan 1-R. Persamaan di atas daplikasikan untuk keadaan normal pada permukaan bidang datar. Dengan asumsi tidak ada sinyal yang hilang sehubungan dengan amplitudo sinyal.
Jejak yang terdapat pada rekaman GPR merupakan konvolusi dari koefisien refleksi dan impulse GPR ditunjukkan oleh persamaan (Kunz dan Luebbers, 1993):
A (t) = r(t)*F(t)*n(t)
(10)
Dimana r(t) adalah koefisien refleksi, A(t) yaitu amplitudo rekaman GPR, F(t) adalah impulse radar, dan n(t) merupakan noise radar.
Besar amplitudo rekaman GPR r(t) akan tampak pada penampang rekaman GPR berupa variasi warna. Refleksi atau transmisi di sekitar batas lapisan menyebabkan energi hilang. Jika kemudian ditemukan benda yang memiliki dimensi yang sama dengan panjang gelombang dari sinyal gelombang
20
elektromagnet maka benda ini menyebabkan penyebaran energi secara acak. Absorbsi ( mengubah energi elektromagnet menjadi energi panas ) dapat menyebabkan energi hilang. Penyebab yang paling utama hilangnya energi karena atenuasi fungsi kompleks dari sifat lstrik dan dielektrika media yang dilalui sinyal radar. Atenuasi (α) tergantung dari konduktivitas (σ), permeabilitas magnetik (μ), dan permitivitas (ε) dari media yang dilalui oleh sinyal dan frekuensi dari sinyal itu sendir (2πf). Sifat bulk dari material ditentukan oleh sifat fisik dari unsur pokok yang ada dan komposisinya (Reynolds, 1996).
3.1.4. Akuisisi data GPR Ada tiga cara penggunaan sistem radar yaitu: reflection profiling (antena monostatik ataupun bistatik), wide-angle reflection and refraction (WARR) atau common-mid point (CMP) sounding, dan transillumination atau radar tomography. Pemilihan cara tersebut di atas tergantung kepada tujuan survei. 3.1.4.1. Radar reflection profiling
Cara ini dilakukan dengan membawa antenna radar bergerak bersamaan diatas permukaan tanah di mana nantinya hasil tampilan pada radargram merupakan kumpulan tiap titik pengamatan (Gambar 3.3).
R
T
Trace 1
T
R
Trace n
reflector
Gambar 3.3. Radar Reflection Profiling (Reynolds, 1997)
21
3.1.4.2. Wide agle reflection and refraction (WARR) atau common mid point
Cara Wide Angle Reflection and Refraction (WARR) Sounding ini dilakukan dengan menaruh transmitter pada posisi yang tetap dan receiver dibawa pada area penyelidikan (Gambar 3.4). WARR sounding diterapkan pada kasus dimana bidang reflector relatif datar atau memiliki kemiringan yang rendah, karena asumsi ini tidak selalu benar pada kebanyakan kasus maka digunakan CMP sounding untuk mengatasi kelemahan tersebut. Pada CMP sounding kedua antenna bergerak menjauhi satu sama lainnya dengan titik tengah pada posisi yang tetap (Gambar 3.5).
Tx
Rx
Rx
Rx
Rx
Gambar 3.4. Wide Angle Reflection and Refraction (Reynolds, 1997)
T3
T2
T1
Rx
R2
Gambar 3.5. Common Mid Point (Reynolds, 1997)
R3
22
3.1.4.3. Transillumination atau radar tomography
Metoda ini dilakukan dengan cara menempatkan transmitter dan receiver pada posisi yang berlawanan. Sebagai contoh jika transmitter diletakan pada satu sisi, maka receiver diletakan pada sisi yang lain dan saling berhadapan. Umumya metoda ini digunakan pada kasus non-destructive testing (NDT) dengan menggunakan frekuensi antenna yang tinggi sekitar 900 MHz (Gambar 3.6).
T
R1
R2
R3
R4
Gambar 3.6. Radar Tomografi (Reynolds, 1997)
Pemilihan frekuensi antena ditentukan oleh tujuan survei yaitu tergantung pada kedalaman atau resolusi yang dibutuhkan dan kondisi material setempat. Semakin
23
tinggi resolusi yang dibutuhkan maka frekuensi antena yang digunakan semakin tinggi (Tabel 3.1 dan Tabel 3.2). Tabel 3.1.
Penentuan frekuensi antena yang digunakan (Lehmann dan Green, 1999).
Frekuensi Antena (MHz)
Ukuran Target (m)
Batas Kedalaman (m)
Kedalaman Penetrasi Maksimum (m)
25 50 100 200 400 1000
≥1 ≥ 0.5 0.1 - 1.0 0.05 - 0.50 ≈ 0.05 Cm
5 – 30 5 – 20 2 – 15 1 – 10 1–5 0.05 – 2
35 – 60 20 – 30 15 – 25 5 – 15 3 – 10 0.5 – 4
Tabel 3.2. Rekomendasi penggunaan sampling frekuensi (Lehmann dan Green, 1999). Frekuensi Antena (MHz) 25 50 100 200 400 1000
Rekomendasi Sampling Frekuensi (MHz) 150 – 600 400 – 800 800 – 1800 1600 – 3500 3200 – 5000 25000 – 110000
Rekomendasi Trace Interval (m) 0.30 - 0.75 0.20 - 0.50 0.10 - 0.30 0.03 - 0.10 0.02 - 0.10 0.01 - 0.05
Proses trigging pada RAMAC/GPR dapat dilakukan dengan cara yang berbedabeda. Pemilihan trigger bergantung pada beberapa faktor, yaitu: 1. Tujuan survei 2. Jenis antena 3. Keperluan untuk proses data lebih lanjut. 4. Peralatan yang mendukung.
24
RAMAC/GPR memberikan 4 macam cara untuk proses trigging ini, yaitu: 1. Dengan menggunakan alat pengukur panjang (Hip Chain atau Measuring Wheel) 2. Melalui keyboard PC eksternal (dengan menekan tombol enter) 3. Dengan menggunakan kotak trig 4. Dengan pembacaan pada interval waktu tertentu. Untuk pengalihan data dari control unit ke dalam external PC, bisa digunakan penghubung secara paralel atau secara serial. Mode paralel dapat mentransfer data lebih cepat dibanding mode serial (www.malags.com).
3.1.5. Processing data GPR Terkadang terdapat beberapa komponen yang dapat merusak sinyal yang di perlukan, sinyal-sinyal ini menyebar secara acak. oleh Karena itu, sinyal tersebut harus di hilangkan untuk meningkatkan mutu hasil pencitraan (Daniels, 1996). Sehingga digunakan beberapa sekuen filtering dalam penerapannya. Sekuen filtering yang akan digunakan untuk pengolahan data GPR ini antara lain:
1. Static correction Proses filtering pertama adalah static correction, filter ini digunakan untuk setiap trace, tidak bergantung satu sama lain. Filter ini digunakan untuk mengoreksi data terhadap elevasi dan waktu tempuh gelombang akibat pengurangan kecepatan.
25
2. subtract-mean (dewow) Dewow merupakan salah satu noise frekuensi rendah yang terekam oleh system. Hal ini terjadi karena instrumen elektronik tersaturasi oleh nilai amplitude besar dari gelombang langsung dan gelombang udara (Van overmeeren, 1997). 3. Gain Filter ini digunakan karena pada lapisan tanah, frekuensi tinggi diserap lebih cepat dibandingkan dengan frekuensi rendah dan terjadi juga spherical divergensi, yaitu energi gelombang yang menjalar berkurang berbanding terbalik dengan kuadrat dari sumber dan hal ini sejalan dengan jarak dan waktu, maka untuk menghilangkannya dilakukan penguatan kembali amplitude yang hilang sehingga seolah-olah di setiap titik energinya sama. Tampilan setelah static correction merupakan Manual Gain yang menggunakkan proses energy decay (Pasasa, 1999).
4. Background Removal Proses pengolahan data dilanjutkan dengan filter Background Removal yang bertindak atas angka terpilih dari lintasan. Filter satu mengurangi lintasan ratarata (tracerange) yaitu memberi jarak jangkauan secara aktual pada suatu bagian . Filter melaksanakan pembersihan latar belakang. Untuk itu, hal ini didefinisikan sebagai waktu / pemberian jarak jangkauan rata-rata.
5. Bandpass butterworth Filter ini berfungsi untuk menghilangkan frekuensi-frekuensi yang tidak diinginkan. Hal ini di maksudkan untuk menghilangkan sinyal horizontal (Hugenschmidt, dkk., 1997).
26
6. F-K filter Dan proses terakhir yaitu F-K filter. Filtering ini berfungsi untuk membatasi area yang akan di filter, dimana Amplitudo spectrum F-K yang terpilih akan memperlihatkan profil asli (Sandmeier, 2012).
Sekuen tersebut merupakan sekuen filter standar untuk data GPR yang merupakan data mentah (raw data). Seluruh proses pengolahan data ini dilakukan berdasarkan tujuan yang ingin dicapai.
3.1.6. Interpretasi GPR Pekerjaan akhir dalam penyelidikan geofisika adalah menerjemahkan data-data sinyal yang telah diperoleh dari akuisisi untuk kemudian diplot ke dalam suatu bentuk konfigurasi agar dapat dibaca dan diambil kesimpulan, pekerjaan ini adalah
interpretasi.
Beberapa
hal
yang
lazim
diperhatikan
dalam
penginterpretasian adalah :
a)
Interpretasi grafik
Kecepatan gelombang dapat diketahui dengan berasumsi pada suatu konstanta dielektrik relative yang mendekati atau sesuai dengan nilai material yang diselidiki, dengan cara demikian two-way travel time (TWT) dapat diterjemahkan menjadi kedalaman, dan jika ditambahkan dengan pengidentifikasian sinyal pantulan dari target (refleksi), maka peta TWT dapat dihasilkan guna menunjukkan kedalaman, ketebalan, perlapisan, dll. Dari sini dapat diketahui nilai sebenarnya dari kecepatan gelombang.
27
b)
Analisa kuantitatif
Dengan
menggunakan
beberapa
analisa,
kedalaman
interpretasi
sinyal
juga kedalaman target atau reflektor dapat dideterminasi tergantung kepada cukup tidaknya nilai yang diketahui dari analisa kecepatan juga variasi konstanta dielektrik relatif material yang dilewati, juga kepada analisa amplitude dan koefisian refleksi.
3.2.
Metoda Geolistrik
3.2.1. Pengertian Metode geolistrik merupakan salah satu metode geofisika yang dapat memberikan gambaran tentang susunan litologi atau struktur bawah permukaan suatu daerah serta kedalaman lapisan batuan berdasarkan sifat kelistrikan batuan (Telford dkk, 1990). Tujuan survey geolistrik tahanan jenis adalah mengetahui perbedaan tahanan jenis (resistivitas) bawah permukaan bumi dengan melakukan pengukuran di permukaan bumi.
3.2.2. Prinsip dasar geolistrik 3.2.2.1. Hukum Coulomb Dua buah medan listrik yang berbeda muatan listriknya akan tarik menarik dan sebaliknya jika bermuatan sama akan tolak menolak. Besarnya gaya interaksi antara dua muatan listrik telah diselidiki oleh Charles Agustin de Coulomb yang menghasilkan hokum suatu muatan uji q diletakan dekat suatu muatan sumber Q yang berjarak r dari muatan uji, maka besarnya gaya yang dialami muatan uji adalah
28
⃑ = ̂
(11)
dimana ⃑ merupakan gaya Coulomb (N), Q adalah muatan sumber (Coulomb), q merupakan muatan uji (Coulomb), r yaitu jarak antar kedua muatan (meter), dan ε0 adalah konstanta permitivitas ruang hampa dengan nilai 8.854 *10-12C2/Nm2. 3.2.2.2. Medan listrik Besarnya medan listrik ⃑ yang ditimbulkan oleh muatan sumber Q adalah:
⃑ =
⃑
̂
(12)
Besar medan listrik dapat dihitung dari persamaan di atas. Sedangkan arahnya yaitu jika muatan sumber Q bermuatan positif, maka arah medan listrik meninggalkan sumber, sebaliknya jika muatan sumber Q bermuatan negative, maka arah medan listrik menuju sumber.
3.2.2.3. Potensial listrik Energi potensial suatu benda adalah kemampuan benda tersebut melakukan kerja. Energi potensial listrik adalah usaha yang diperlukan untuk memindahkan muatan q dari titik tak berhingga ke titik tersebut.
= ∫∝ ⃑ . ⃑ =
(13)
Sedangkan potensial listrik (V) adalah energy potensial U per satuan muatan uji.
29
(14) Energi potensial dan potensial listrik merupakan besaran skalar, sedangkan gaya Coulomb dan medan listrik merupakan besaran vector.
3.2.2.4. Hukum Ohm Kemampuan benda untuk menahan arus listrik berbeda-beda. Hukum Ohm memberikan gambaran hubungan antara besarnya beda potensial listrik (V), kuat arus (I), dan besarnya tahanan listrik kawat penghantar (R). Hal ini diperlihatkan pada Gambar 3.7.
Gambar 3.7. Resistansi (sumber: http://www.predict.its.ac.id). Hubungan tersebut adalah
(15) Sedangkan resistivitas (ρ) berhubungan dengan luas (A) dan panjang (L) mediumnya (Gambar 3.8.).
Gambar 3.8. Resistivitas (sumber:http://www.predict.its.ac.id.).
30
Resistivitas merupakan turunan dari resistansi, yaitu:
=
(16)
3.2.2.5. Dua sumber titik arus pada permukaan Sumber arus ganda pada permukaan memungkinkan bagi kita untuk menganalisa sifat kelistrikan batuan yang ada diantara dua sumber arus tersebut sebagaimana diperlihatkan pada Gambar 3.9.
Gambar 3.9. Arus listrik yang di timbulkan oleh dua sumber (Telford, dkk., 1990).
Potensial di setiap titik terdekat permukaan akan dipengaruhi oleh kedua elektroda arus, sehingga equipotensial yang dihasilkan dari dua titik sumber dekat permukaan akan dipengaruhi oleh kedua elektroda tersebut. Secara jelas dapat dilihat pada Gambar 3.10.
Gambar 3.10. Dua elektroda arus dan potensial pada permukaan homogen isotropis (Telford, dkk., 1990).
31
Dari Gambar 3.10, maka di dapatkan persamaan berikut ini:
V 2 1 1 I 1 1 r r R R 1 2 1 2 2 K 1 1 1 1 r r R R 1 2 1 2
(17)
(18)
dimana V merupakan beda potensial antara P1 dan P2 (volt), K adalah faktor geometri yang merupakan besaran yang berubah terhadap jarak spasi elektroda dan bergantung pada susunan elektroda, I
adalah kuat arus yang dialirkan
melalui elektroda arus C1 dan C2 (ampere), r1 yaitu jarak antara C1 dan P1 (meter), r2 adalah jarak antara C2 dan P1 (meter), R1 merupakan jarak antara C1 dan P2 (meter), dan R2 adalah jarak antara C2 dan P2 (meter).
K
V I
(19)
3.2.3. Akuisisi data geolistrik 3.2.3.1. Teknik pengukuran Teknik pengukuran dalam geolistrik adalah sebagai berikut (Barker, 1981):
1. Sounding : untuk informasi bawah permukaan secara vertikal (model bumi berlapis)
2. Profilling : untuk informasi bawah permukaan secara mendatar (variasi lateral)
32
3. Offset Sounding : untuk informasi bawah permukaan profil sounding yang kontinyu secara lateral
3.2.3.2. Tahapan akusisi Tahapan dalam akuisisi geolistrik adalah:
1. Tentukan konfigurasi elektroda yang ingin dipakai. Dalam pelaksanaan survey dikenal beberapa metoda pengambilan data sesuai dengan peletakan eloktroda yang dilakukan (Drecun, 1976). Hal ini berpengaruh terhadap faktor geometri peneletian resistivity yang kita lakukan. Adapun aturan/metoda tersebut antara lain :
a. Metoda Wenner b. Metoda Gradien c. Metoda Schlumberger d. Metoda Dipole-dipole e. Metoda Pole-dipole
2.
Pasang elektroda sesuai dengan konfigurasi yang di pilih. Pada hal ini,
konfigurasi yang digunakan adalah konfigurasi wenner-schlumberger. Dimana konfigurasi Wenner-Schlumberger adalah konfigurasi dengan sistem aturan spasi yang konstan dengan catatan faktor “n” untuk konfigurasi ini adalah perbandingan jarak antara elektroda C1-P1 (atau C2-P2) dengan spasi antara P1-P2 seperti pada Gambar 3.11. Jika jarak antar elektroda potensial (P1 dan P2) adalah a maka jarak antar elektroda arus (C1 dan C2) adalah 2na + a. Proses penentuan
33
resistivitas menggunakan 4 buah elektroda yang diletakkan dalam sebuah garis lurus.
Gambar 3.11. Pengaturan elektroda konfigurasi Wenner-Schlumberger
3.
Ukur besar resistivity semunya. Dalam hal ini Resistivitymeter memberikan nilai resistansi R = V/I sehingga nilai resistivitas dapat dihitung dengan:
Gambar 3.12. Bentuk konfigurasi Wenner-Schlumberger
Bumi tersusun atas lapisan-lapisan tanah yang nilai resistivitas suatu lapisan tanah atau batuan tertentu berbeda dengan nilai resistivitas lapisan tanah atau batuan lainnya (Gambar 3.12). Nilai resistivitas ini dapat diketahui dengan menghubungkan battery dengan sebuah Ammeter dan elektroda arus untuk mengukur sejumlah arus yang mengalir ke dalam tanah, selanjutnya ditempatkan dua elektroda potensial dengan jarak a untuk mengukur perbedaan potensial antara dua lokasi.
34
3.2.4. Processing data geolistrik Dalam bagian ini perangkat lunak yang akan digunakan untuk pengolahan data geolistrik adalah Res2dinv, yang meliputi langkah – langkah penginputan data, read data, inversi, sampai pada pencetakan hasil inverse, dan tampilan topografi (Grifiths, 1993). Namun, sebelum itu dilakukan pengolahan data geolistrik hasil pengukuran di lapangan terlebih dahulu dengan menggunakan Ms. Excel.
Selanjutnya buat worksheet baru kemudian copy dp (datum point), a (spasi terkecil), dan Rho ke worksheet baru sebagai sarana dalam pembacaan di RES2DINV ( Loke, 1997). Sebelum proses inversi dilakukan maka terlebih dahulu dilakukan proses pembacaan data dengan langkah – langkah: file-read data file. Untuk melakukan inverse data maka pada menu inversion, pilih "use logarithm of apparent resistivity" sehingga muncul kotak dialog “use logarithm of apparent resistivity”, pilih “Use apparent resistivity” kemudian “OK” untuk memulai proses inversi.
Pada jendela yang baru ini, pada kolom “file name” input data yang telah di read data file tadi (latihan.inv), kemudian “save” maka akan muncul 3 model, dimana model pertama adalah data awal, model kedua adalah data inverse, dan model ketiga adalah hasil inversi. Jika tampilan yang didapat tidak sesuai dengan kondisi geologi yang ada, dapat melakukan editing data dengan cara klik pada Edit / terminate bad datum points / OK. Penghilangan data ini dimaksudkan untuk menghilangkan data yang kurang baik sehingga dapat mengganggu model dan menurunkan RMS Errors yang diperoleh.
35
Jika model yang diperoleh sudah memuaskan kita dapat menyimpannya dalam file gambar (*.bmp) secara permanen.
3.2.5. Interpretasi hasil pengolahan data geolistrik
Langkah dalam menentukkan jenis batuan yang terlihat pada tampilan hasil processing data geolistrik yaitu dengan mengkorelasikan nilai resistivitas yang didapat dengan tabel resistivitas yang ada (Tabel 3.3).
Tabel 3.3. Resistivitas material-material bumi (Telford, dkk., 1990) Material Resistivity Pyrite (Pirit) Quartz (Kwarsa)
(Ohm-meter) 0,01 – 100 500 – 800.000
Calcite (Kalsit)
1 x 1012 – 1 x 1013
Rock Salt (Garam Batu) Granite (Granit)
30 – 1 x 1013 200 – 100.000
Andesite (Andesit) Basalt (Basal) Limestones (Gamping) Sandstones (Batu Pasir) Shales (Batu Tulis) Sand (Pasir) Clay (Lempung) Ground Water (Air Tanah) Sea Water (Air Asin) Magnetite (Magnetit) Dry Gravel (Kerikil Kering) Alluvium (Aluvium) Gravel (Kerikil)
1,7 x 102 – 45 x 104 200- 100.000 500 – 10.000 200 – 8.000 20 – 2.000 1 – 1.000 1 – 100 0.5 – 300 0.2 0.01 – 1.000 600 – 10.000 10 – 800 100 – 600
36
3.3.
Metode Gravity
3.3.1. Pengertian
Metoda Gravity merupakan suatu metoda eksplorasi geofisika yang didasarkan atas adanya anomali medan gravitasi bumi, akibat adanya variasi densitas batuan ke arah lateral maupun vertikal dibawah titik ukur. Pemodelan pada metoda ini dilakukan berdasarkan atas fungsi variasi densitas (density) dan kedalaman z (Sarkowi, 2010).
3.3.2. Pengukuran gravity
1. Pengukuran absolut biasanya dilakukan di laboratorium-laboratorium, sukar untuk mendapatkan harga gravity absolut yang akurat, karena banyak kendalakendala yang sangat mempengaruhi hasil pengukuran. Oleh karena itu, pengukuran secara absolut jarang sekali dilakukan karena terlalu sukar dan melibatkan banyak faktor maupun alat (Grant dan West, 1969). 2. Pengukuran relatif lebih umum dan mudah dilakukan, pada penelitian gravity (Gambar 3.13). Pengukuran relatif dilakukan dengan membandingkan hasil pengukuran titik yang tidak diketahui nilai gayaberatnya dengan titik yang sudah diketahui yang telah diikat kepada titik-titik referensi (Postdam, IGSN dsb).
37
A
B
Pembacaan alat s1 Pada Titik A yang telah di ketahui harga gayaberatnya (ga)
Pembacaan alat s2 Pada titik-titik yang di ketahui harga gb
Gambar 3. 13. Pengukuran gravity secara relatif
3.3.3. Koreksi gravity
a. Koreksi tide Koreksi tidal (koreksi pasang surut) koreksi ini timbul disebabkan efek tarikan massa yang disebabkan oleh benda-benda langit, terutama bulan dan matahari. Harga koreksi ini berubah-ubah setiap waktu secara periodik tergantung dari kedudukan benda-benda langit tersebut. Jika koreksi ini merupakan gaya tarik bulan dan matahari pada permukaan bumi maka harga baca tersebut ditambahkan pada harga baca dan pengamatan. Dan jika koreksi tersebut merupakan lawan dari gaya tarik maka perlu dikurangkan (Goodkind, 1986). Grav+Tide = Grav.Observasi + Koreksi Tide
(20)
38
b. Koreksi Drift Nilai bacaan gayaberat pada titik yang sama setelah dikoreksi pasang surutgravity, kadang memberikan nilai yang berbeda baik lebih kecil maupun lebih besar (Longman, 1959).
Perbedaan nilai bacaan setelah dikoreksi
pasang surut ini biasanya disebut dengan koreksi drift. Perbedaan ini disebabkan oleh sensor pegas pada gravimeter mengalami kelelahan akibat goncangan, pemakaian terlalu lama dan lain-lain. Untuk menghindari atau mengoreksi perbedaan ini :
c.
Pengukuran gravity dilakukan pada lintasan tertutup
Koreksi Drift diasumsikan linier terhadap waktu.
Koreksi gravity observasi Untuk mendapatkan gravity observasi maka nilai gayaberat hasil pengukuran setelah dikoreksi tide dan drift harus diikatkan pada suatu titik yang sudah diketahui nilai gravitynya. Gobs = g obs base st + Δ g mGal
d.
(21)
Koreksi lintang Koreksi ini dilakukan karena bentuk bumi yang tidak bulat sempurna, terdapat perbedaan antara jari-jari bumi di kutub dengan di katulistiwa. Nilai gravity dikutub akan lebih besar dibandingkan nilai gravity di katulistiwa. Nilai gravity teoritik pada lintang menurut Moritz (1980) :
g 978.032,7(1 5.3024 10 -3 Sin 2 5,8 10 -6 Sin 2 2 ) mGal
(22)
39
e.
Koreksi udara bebas (Free Air Correction) Koreksi udara bebas ini dilakukan terhadap titik-titik pengukuran yang terletak pada ketinggian h dari permukaan air laut. Koreksi gravity yang dihitung dari persamaan gravity normal bumi dengan bentuk ellipsoid, Sehingga koreksi udara bebas di permukaan bumi memiliki nilai :
FAC = -0.308765 h
f.
(23)
Koreksi Bouger
Pada koreksi udara bebas perhitungan dilakukan dengan menghitung selisih ketinggian titik amat dengan datum tanpa memperhitungkan adanya massa di antara titik tersebut. Untuk itu diperlukan koreksi untuk menghilangkan pengaruh adanya massa dari datum sampai ketinggian titik tersebut . Nilai koreksi Bouguer untuk model slab tak hingga diberikan oleh persamaan :
g B 41.93h
(24)
dimana: gB adalah koreksi bouguer dengan satuan microGall, merupakan densitas permukaan dengan satuan gr/cc, adalah porositas, dan
h
merupakan ketinggian dengan satuan cm.
g.
Koreksi terrain
Koreksi ini merupakan koreksi akibat adanya efek massa disekitar titik observasi misalnya gunung, gedung, limbah dan lain-lain. Adanya gunung atau massa menonjol keluar akan menarik titik di permukaan g(R) ke atas
40
sehingga akan mengurangi g(R), demikian juga lembah akan mengurangi tarikan ke bawah pada g(R) yang berarti mengurangi g(R).
CBA = SBA + TC
(25)
Dimana: CBA adalah complete bouguer anomaly, SBA yaitu simple bouguer anomaly, dan TC adalah terrain correction.
