PENGOLAHAN DATA GRAVITY Pengolahan Data Complete Bouguer Anomali
1. Menghit Menghitung ung elevas elevasii sebenar sebenarnya nya (denga (dengan n GPS): GPS): Koreksi drift elevasi : • menit menit ke i−menit awal Drift = menit akhir−menit awal (elevasi akhir ! elevasi a"al)
III.1
#levasi terkoreksi: #levasi terkoreksi = elevasi $e%ba&aan alat ! drift III.2 Delta elevasi • Delta elevasi = elevasi di titik i ! elevasi di d i titik a"al III.3 #levasi sebenarnya • III.4 #levasi sebenarnya = delta elevasi ' elevasi referensi . Menghit Menghitung ung koreks koreksii $asa $asang ng suru surutt (Tide (Tide Correction) Correction) ilai $er&e$atan gravitasi yang diba&a oleh *lat gravi%eter gravi%eter *utograv *utograv S&ientre +G,- sudah •
terkoreksi tide. Sehingga nilai koreksi tidenya adalah . /. Menghit Menghitung ung $e%ba $e%ba&aan &aan g setel setelah ah koreks koreksii tide tide Pe%ba&aan koreksi tide = $e%ba&aan ' Tide Correction 0. Menghit Menghitung ung kore koreksi ksi drift drift $e%b $e%ba&aa a&aan n nilai nilai g (D+) (D+) menit kei −menit awal Drift = menit akhir−menit awal (GSakhir ! ! GS)
III.5
III.6
Di%ana2 GS = 3a&aan gravitasi terkoreksi $asut di 3S a"al GSakhir = = 3a&aan gravitasi terkoreksi $asut di 3S akhir -. Menghitung Menghitung nilai nilai $er&e$atan $er&e$atan gravit gravitasi asi (g) terkor terkoreksi eksi : Pe%b. erkoreksi erkoreksi = $e%ba&aan terkoreksi setelah Tide Correction , Correction , D+ 4. Meng Menghi hitu tung ng del delta ta g : Delta g = $e%ba&aan g setelah D+ ! $e%ba&aan g referensi 5. Menghit Menghitung ung g real real observ observasi asi : G real obs = g real referensi ' delta G 6. Mengh Menghit itun ung g korek koreksi si lint lintan ang g
III.7 III.8
III.9
2θ 2
1 + 0.0053024 0.0053024 sin θ− 0.0000058sin ❑¿
III.10
¿
g❑=978032.7 ¿ Di%ana2
θ
= latitude (radian)
7. Menghit Menghitung ung udara udara bebas bebas (8** (8**)) 8*+ 8*+ = 2/64 h
III.11
Di%ana2 h = elevasi sebenarnya. 1. Menghitung koreksi terrain (+): 9ika nilai dketahui dari data geologi : • 0.04191 ρ 2 2 2 2 R d + h – √ Rl + h √ + = (; ! ; ) ' < d N
III.12
Ket : ; < = ari luar ;d = ari dala% = u%lah sektor %ata angin h = ketinggian untuk setia$ arah %ata angin 9ika nilai tidak diketahui : • Digunakan %etode ha%%er &hart untuk esti%asi densitas batuan yang do%inan. 11. Menghitung koreksi 3ouguer (3+): 3+ = >Gh Karena2 >G = 2 0171 %aka 3+ = 20171h
III.13
ilai bisa di$eroleh dengan %enggunakan data geologi atau dengan %enggunakan %etode nettleton atau $arasnis. 1. Menghitung nilai +o%$lete 3ouguer *no%ali (+3*) *3< = gobs ! glintang '8** ! 3+ ' + Pengolahan Data Lebih Lanjut
Analisis Spektral 1. 3uka surfer . ?n$ut nilai $ada kolo% (latitude)2 y (longitude)2 @ (+3*).
/. Si%$an dala% bentuk +3*.bln 0. 3uat $lot baru dengan %e%ilih icon -. 3uat grid data +3* dengan tools Grid
Data
III.14
4. Pilih data +3*.bln ke%udian klik open 5. Pilih sheet di%ana data telah selesai diolah hingga mendapatkan Complete Bouguer Anomali. 6.
7. Setelah %elakukan gridding %aka dilakukanlah countouring 2 &aranya dengan %eng,klik tools Ma$s e"
+ountour Ma$ atau klik 2 lalu $ilih data +3*.grd sehingga ta%$il $eta
Complete Bouguer Anomali.
1. Mun&ulkan titik,titik $engukuran dengan $ilih icon
11. Setelah $roses countouring 2 selanutnya dilakukan $roses slicing . Perta%a buat garis $ada daerah yang akan di slice. Sebaiknya garis tersebut %engikuti lintasan agar tidak banyak data yang terinter$olasi. lalu klik Ma$
Digiti@e atau klik kanan ga%bar $ilih digiti@e.
1. Klik %asing,%asing uung garis sehingga %un&ul koordinat hasil digitasi.
1/. Save hasil digitasi terbut. 10.
yang akan di Sli&e dala% hal ini +3*.grd lalu $ilih file hasil digitasi ke%udian lakukan sli&ing. Save for%at dala% bentuk sli&e.bln dan sli&e.D*. 1-. 3uka folder te%$at $enyi%$an file hasil slice.
Gambar III.12 8or%at
hasil sli&e $ada #&el dengan s$asi baru dan +3* baru 15. +o$y nilai s$asi baru dan &ba baru ke note$ad. Save dala% bentuk sli&e.BC dan save data tersebut di folder yang sa%a dengan soft"are M#;?.#B# agar da$at terba&a $ada soft"are M#;?.#B#. 16. 3uka M#;?.#B# dan tekan 8-.
17. ekan - untuk %e%ilih ransfor%asi 8ourier Diskrit2 lalu tekan 12 1 dan sehingga %e%ilih Edata dari disketFhardisk
.
1. Kemudian klik Esc, Esc, Esc sehingga muncul pilihan sebagai berikut
. Pilih /2 -2 2 / dan ketikkan out$ut file yang diinginkan 2 lalu tekan 81 untuk %enyi%$an data.
/. Setelah %enda$atkan data yang telah ditransformasi fourier dala% bentuk .SPK %aka buka dengan ote$ad dan +o$y2 Paste ke dala% Mi&rosoft #&el. 0. Susun %enggunakan Data
et to +olu%n sehingga %enadi ra$i2 lalu Cut data $aling atas
dan Paste ke kolo% baru $ada baris ke / dan beri nilai nol $ada baris atasnya. -. 3eri label untuk tia$ kolo%. ilai yang ada beru$a bilangan real2 i%aginer2 dan frekuensi s$atial (8). 4. 3uat kolo% baru untuk s$ektru% a%$litudo (*)2 bilangan gelo%bang (k)2 log natural s$ketru% a%$litudo (ln *)
A = √ real + imaginer 2
2
k = >8
5. 3uat korelasi antara bilangan gelo%bang (k) dengan logarit%a natural a%$litude (ln *) sehingga ta%$il diagra% se$erti berikut :
6. 3agian sebelah kanan %eru$akan $en&er%inan dari bagian sebelah kiri2 sehingga untuk efisiensi $erkirakan data untuk regional dan residual.
RESIDUAL REGION
7. Plot data dan lakukan $enentuan titik $otong $ada grafik2 tentukan titik &ut of dan u%lah "indo".
+ut off = (&,&1)F(%1,%)
" = H$iF(H&ut off)
Pemisahan Anomali Regional dan Anomali residual (Filtering) 1. Data +3* yang telah selesai gridding $ada surfer %aka selanutnya akan dilakukan filterisasi untuk %enentukan ano%ali gravitasi regionalnya. Ial ini da$at dilakukan dengan tool filter2 &aranya klik Grid
8ilter
. Setelah itu %un&ul "indo" baru2 $ilih %etode filterasi yang akan dilakukan2 $ada hal ini digunakan Moving *verage dengan ukuran "indo" 17172 isi ;o"s 17 dan +ol 172 ganti out$ut file dengan na%a ;egional.grd ke%udian tekan AK.
/. 3uka file grid yang telah dilakukan filtering tadi2 dala% hal ini file ;egional.grd2 lalu klik A$en. Maka dida$atlah hasil kontur $eta *no%ali 3ouguer ;egional.
0. ntuk %e%$eroleh $eta *no%ali 3ouguer ;esidual2 dilakukan $engurangan $eta +o%$lete 3ouguer *no%ali dengan $eta *no%ali 3ouguer ;egional. +ara yang da$at dilakukan adalah dengan %engklik Grid
Math
-. Pilih file +3*.grd (*) dan ;egional.grd (3)2 lalu klik $ilih for%at fungsi *,3. Save out$ut file dengan na%a ;esidual.grd. Sehingga di$eroleh hasil $eta ano%ali bouguer residual.
4. 3uat sli&e $ada $eta ano%ali residual. 5. 3uka hasil sli&eres.D* $ada e&el. Data yang di&uat beru$a nilai 2 y2 arak2 dan ano%ali bouguer residual. 6. +o$y ano%ali bouguer residual dan arak ke note$ad.
7. Si%$an dala% bentuk sli&eres.dta dala% folder Gravd&. 1. 3uka soft"are G;*JD+.#B#. 11. ?n$ut file dala% bentuk .dta tadi.
1. Ke%udian $ilih file yang akan di%odelkan.
Sehingga %un&ul ta%$ilan se$erti berikut
1/. 8ile tersebut telah sia$ untuk di%odelkan.
10. 3uat sinkronisasi antara garis $utus,$utus (garis kalkulasi) dengan garis observasi dengan %e%buat body se%barang untuk %e"akili densitas batuan. Pilih add body untuk %ena%bahkan body. #dit densitasnya agar errornya berkurang sehingga %enghasilkan $e%odelan gravity dengan densitas yang beraga% untuk %engetahui la$isan $ada ba"ah $er%ukaan. 1-. ntuk lebih %e%$erelas densitasnya2 "arna body bisa diubah Dis$lay +hoose +olour S&he%e $ilih &olour skala.
Email auliakahar!"#gmail$c%m