3 0 3 MB
STUDI AWAL PERCEPATAN GETARAN TANAH GEMPA LOMBOK TAHUN 2018 BERDASARKAN MODEL ATENUASI DAN DATA AKSELEROGRAF
SKRIPSI
Oleh
HASTI WAHYU CAHYANINGRUM K1C015025
KEMENTERIAN RISET TEKNOLOGI DAN PENDIDIKAN TINGGI UNIVERSITAS JENDERAL SOEDIRMAN FAKULTAS MATEMATIKA DAN ILMU PENGETAHUAN ALAM PURWOKERTO 2019
STUDI AWAL PERCEPATAN GETARAN TANAH GEMPA LOMBOK TAHUN 2018 BERDASARKAN MODEL ATENUASI DAN DATA AKSELEROGRAF
SKRIPSI
Oleh
HASTI WAHYU CAHYANINGRUM K1C015025
Sebagai Salah Satu Persyaratan untuk Memperoleh Gelar Sarjana Strata Satu (S1) pada Jurusan Fisika Fakultas Matematika dan Ilmu Pengetahuan Alam Universitas Jenderal Soedirman
KEMENTERIAN RISET TEKNOLOGI DAN PENDIDIKAN TINGGI UNIVERSITAS JENDERAL SOEDIRMAN FAKULTAS MATEMATIKA DAN ILMU PENGETAHUAN ALAM PURWOKERTO 2019
i
LEMBAR PENGESAHAN
STUDI AWAL PERCEPATAN GETARAN TANAH GEMPA LOMBOK TAHUN 2018 BERDASARKAN MODEL ATENUASI DAN DATA AKSELEROGRAF
Oleh HASTI WAHYU CAHYANINGRUM K1C015025
Diterima dan disahkan Pada tanggal : ……………………………
Pembimbing I
Pembimbing II
ii
Sukmaji Anom Raharjo, M.Si
Sigit Pramono, M.Si.
NIP. 19780417 200812 1 002
NIP. 19820608 200312 1 002
Mengetahui Dekan Fakultas Matematika dan Ilmu Pengetahuan Alam
Drs. Sunardi, M.Si. NIP.19590715 199002 1 001
PERNYATAAN Dengan ini saya menyatakan bahwa skripsi yang berjudul:
STUDI AWAL PERCEPATAN GETARAN TANAH GEMPA LOMBOK TAHUN 2018 BERDASARKAN MODEL ATENUASI DAN DATA AKSELEROGRAF adalah benar merupakan hasil karya saya sendiri dan semua sumber data serta informasi yang digunakan telah dinyatakan secara jelas dan dapat diperiksa kebenarannya. iii
Bila pernyataan ini tidak benar, maka saya bersedia menerima sanksi pencabutan gelar kesarjanaan yang telah saya peroleh.
Purwokerto,
November 2019
Hasti Wahyu Cahyaningrum K1C015025
iv
PEDOMAN PENGGUNAAN SKRIPSI
Skripsi ini terdaftar dan tersedia di Pusat Informasi Ilmiah Fakultas Matematika dan Ilmu Pengetahuan Alam dan terbuka untuk umum dengan ketentuan bahwa hak cipta ada pada penulis dengan mengikuti aturan HaKI yang berlaku di Universitas Jenderal Soedirman. Pengutipan dan atau peringkasan hanya dapat dilakukan dengan mengikuti kebiasaan ilmiah untuk menyebutkan sumbernya.
v
KATA PENGANTAR Puji syukur kehadirat Allah SWT yang telah memberikan rahmat serta karuniaNya sehingga penulis dapat menyelesaikan tugas akhir dengan judul “Studi Awal Percepatan Getaran Tanah Gempa Lombok Tahun 2108 Berdasarkan Model Atenuasi dan Data Akselerograf” dengan baik dan lancar, bukan karena kuat dan gagah manusia, namun karena penyertaan-Nya semua ini dapat dilalui dengan penuh rasa syukur. Pada kesempatan ini penulis mengucapkan terima kasih kepada: 1. Allah Yang Maha Esa yang telah memberikan rahmat dan karunia-Nya kepada penulis. 2. Orangtua dan keluarga yang senantiasa memberikan do’a, kasih sayang, nasihat, dan dukungan baik secara moril maupun materil. 3. Bapak Sukmaji Anom Raharjo, M.Si selaku pembimbing I yang selalu memberikan
ilmu,
saran,
dan
bimbinganya
kepada
penulis
dalam
menyelesaikan skripsi. 4. Bapak Sigit Pramono, M.Si selaku pembimbing II yang selalu memberikan ilmu, saran, dan bimbinganya kepada penulis dalam pengolahan data skripsi. 5. Bapak R. Farzand Abdullatif, Ph.D. selaku penelaah yang memberikan masukan dalam penulisan skripsi ini. 6. Bapak Jamrud Aminuddin, Ph.D. dan bapak Sehah, M.Si selaku penguji I dan penguji II yang telah mengevaluasi penyusunan skripsi. 7. Seluruh dosen dan staf Jurusan Fisika yang tidak dapat saya sebutkan satu per satu. 8. Seluruh staf bidang Seismologi Teknik Badan Meteorlogi Klimatologi, dan Geofisika Jakarta Pusat atas dukungan dan sharing ilmu selama penelitian 9. Teman-teman satu penelitian dari Universitas Indonesia dan Universitas Brawijaya yang telah menjadi partner sharing ilmu selama penelitian di Bidang Seismologi Teknik, Badan Meteologi Klimatologi, dan Geofisika. vi
10. Teman-teman Gravitasi 2015 yang telah membersamai selama empat tahun dalam proses belajar di fisika Unsoed. Penulis menyadari bahwa keterbatasaan kemampuan dan pengetahuan dalam menyusun skripsi ini masih jauh dari sempurna. Oleh karena itu, kritik dan saran yang bersifat membangun sangat diharapkan penulis guna perbaikan dimasa yang akan datang. Semoga skripsi ini dapat bermanfaat bagi semua pihak yang membutuhkannya. Terima kasih.
vii
DAFTAR ISI
HALAMAN JUDUL................................................................................................i HALAMAN PENGESAHAN.................................................................................ii PERNYATAAN.....................................................................................................iii PEDOMAN PENGGUNAAN SKRIPSI................................................................iv KATA PENGANTAR.............................................................................................v DAFTAR ISI.........................................................................................................vii DAFTAR SINGKATAN DAN LAMBANG..........................................................x DAFTAR GAMBAR.............................................................................................xii DAFTAR TABEL...................................................................................................xi DAFTAR LAMPIRAN.........................................................................................xii ABSTRAK............................................................................................................xiii ABSTRACT............................................................................................................xiv BAB I PENDAHULUAN.......................................................................................1 1.1. Latar Belakang....................................................................................1 1.2. Perumusan Masalah............................................................................2 1.3. Cakupan dan Batasan Masalah...........................................................2 1.4. Tujuan Penelitian................................................................................3 1.5. Manfaat Penelitihan............................................................................3 BAB II TINJAUAN PUSTAKA............................................................................4 2.1. Gelombang Seismik............................................................................4 2.2. Gempa Bumi.......................................................................................7 2.3. Tinjauan Seismistas Lombok..............................................................8 2.4. Akselerograf......................................................................................10 2.5. Percepatan Tanah..............................................................................12 2.6. Faktor yang mempengaruhi persamaan atenuasi gerakan tanah.......13 2.4.1. Magnitudo gempa (earthquake magnitude)...........................14 2.4.2. Jarak ke site............................................................................16 viii
2.4.3. Pengaruh mekanisme sumber gempa (source mechanism)....16 2.4.4. Pengaruh kondisi situs (lokal site condition).........................16 2.7. Prediction Equation (GMPE)...........................................................17 2.8. Analisis Regresi................................................................................19 2.9. Analisis Korelasi...............................................................................20 2.10. Intensitas Seismik.............................................................................21 BAB III METODE PENELITIAN.......................................................................22 3.1
Tempat dan Waktu............................................................................22
3.2
Alat dan Bahan..................................................................................22
3.3
Tahapan.............................................................................................22 3.3.1 Tahapan menghitung jarak stasiun ke sumber gempa ( Rhiposenter ¿..........................................................................22 3.3.2 Analisi regresi untuk menghasilkan persamaan atenuasi.......23 3.3.3 Validasi...................................................................................24 3.3.4 Penerapan persamaan atenuasi dalam pemetaan wilayah terdampak akibat gempa bumi................................................25
3.4
Diagram Alir Penelitian....................................................................26
BAB IV HASIL DAN PEMBAHASAN..............................................................28 4.1
Menentukan persamaan empiris percepatan getaran tanah...............28 4.1.1. Pemilihan data katalog gempa bumi di wilayah Lombok dan sekitarnya................................................................................28 4.1.2. Pemilihan data percepatan getaran tanah hasil analisis rekaman..................................................................................29 4.1.3. Hasil penentuan persamaan empiris percepatan tanah...........30
4.2
Validasi hasil percepatan tanah model dengan percepatan tanah observasi...........................................................................................31 4.2.1. Hasil analisis model................................................................31 4.2.2. Perbandingan nilai percepatan tanah model dengan metode percepatan tanah Mc Guire (1997).........................................33
4.3
Perbandingan antara Peta Kontur PGA Observasi, PGA Model Atenuasi, PGA Mc Guire (1997)......................................................34 ix
BAB V KESIMPULAN DAN SARAN................................................................37 5.1
Kesimpulan.......................................................................................37
5.2
Saran.................................................................................................37
DAFTAR PUSTAKA............................................................................................38 LAMPIRAN...........................................................................................................41
x
DAFTAR SINGKATAN DAN LAMBANG
SINGKATA N PGA
Peak Ground Acceleration
USGS
United States Geological Survey
BMKG
Badan Meteorlogi, Klimatologi, dan Geofisika
PSHA
Probabilistik Seismic Hazard Analysis
DSHA
Deterministik Seismic Hazard Analysis
TWSI
Sensor akselerograf di daerah Taliwang
MASE
Sensor akselerograf di daerah Mataram
PHA
Peak Horizontal Acceleration
GMPE
Ground Motion Prediction Equation
MMI
Modified mercalli intensity
SIGMA
Seismic Ground Motion Assessment
LAMBANG b1
Konstanta model
b2
Koefisien magnitudo
b3 Y ML Mw mb Ms Rh Repic h
Koefisien jarak hiposenter Nilai percepatan tanah Magnitudo Lokal Magnitudo momen Magnitudo body Magnitudo permukaan Jarak hiposenter Jarak episenter Kedalaman centimeter
cm xi
s
sekon
km
kilometer
g
Satuan percepatan gravitasi
gal
Satuan percepatan tanah
d
Data
m
Model
G
Matrix kernel
R2
R square
xii
DAFTAR GAMBAR Gambar 2.1 Refleksi Refraksi Gelombang Seismik (Gadallah, 2009)...................5 Gambar 2.2 Prinsip Huygens (Gadallah, 2009).....................................................6 Gambar 2.3 Prinsip Fermat (Gadallah, 2009)........................................................6 Gambar 2.4 Deformasi batuan akibat berbagai bentuk stress (Sunarjo, 2012)......8 Gambar 2.5 Jaringan akselerograf di wilayah Lombok.......................................12 Gambar 2.6 Korelasi antara ML dan Mw..............................................................15 Gambar 2.7 Pengaruh kondisi tanah terhadap PHA........................................16YY Gambar 3.1 Ilustrasi jarak stasiun ke sumber gempa...........................................23 Gambar 3.2 Titik grid daerah penelitian..............................................................25 Gambar 3.3 Diagram alir penelitian.....................................................................27 Gambar 4.1 Sebaran Gempa bumi Lombok tahun 2018......................................29 Gambar 4.2 Perbandingan nilai PGA hasil observasi, model, dan Mc Guire (1997) untuk magnitudo 4,5..............................................................32 Gambar 4.3 Perbandingan nilai PGA hasil observasi dan model.........................33 Gambar 4.4 Perbandingan nilai PGA antara rekaman akselerograf dan model...34 Gambar 4.5 Peta pecepatan tanah dari data rekaman akselerograf (BMKG)......36 Gambar 4.6 Peta pecepatan tanah dari hasil model atenuasi................................36
xiii
DAFTAR TABEL YTabel 2.1 Jaringan Akselerograf di Wilayah Lombok............................................
Tabel 2.2 Jenis-jenis Magnitudo...........................................................................14 Tabel 2.3 Korelasi konversi magnitudo untuk wilayah Indonesia........................15 Tabel 2.4 Klasifikasi jenis tanah berdasarkan nilai Vs30 menurut NEHRP.........17 Tabel 2.5 Interpretasi Koefisien Korelasi..............................................................21 Tabel 2.6 Intensitas Seismik21 Y
Tabel 4.1 Nilai PGA hasil rekaman akselerograf..................................................30 Tabel 4.2 Koefisien hasil model............................................................................32
xi
DAFTAR LAMPIRAN YLampiran A Peta Geologi Lombok........................................................................
Lampiran B Data parameter gempa rekaman akselerograf..................................43 Lampiran C Hasil Koefisien tiap Residual..........................................................56 Lampiran D Nilai PGA hasil rekaman akselerograf dari Shakemap untuk.........56 Lampiran E Nilai PGA hasil rekaman akselerograf dari SIGMA untuk............58 Lampiran F Riwayat Hidup..................................................................................59
xii
ABSTRAK Percepatan getaran tanah merupakan salah satu parameter yang sangat berperan dalam menentukan tingkat kerusakan tanah dan bangunan yang terjadi di permukaan bumi akibat guncangan gempa bumi. Belum adanya persamaan empiris yang dikembangkan untuk wilayah Lombok serta kekhawatiran akan dampak gempa bumi besar dan merusak, mendorong penelitian untuk mengembangkan suatu persamaanan fungsi atenuasi percepatan getaran tanah maksimum. Pada penelitian ini, dilakukan penentuan model persamaan empiris PGA berdasarkan persamaanan empiris umum Mc Guire (1997) dengan menggunakan 544 data dari sensor akselerograf yang bersesuaian dengan event gempa bumi sekitar wilayah Lombok dengan periode tahun 20018. Berdasarkan analisis regresi, maka diperoleh koefisien geometrical spreading b 2 = 0,712, koefisien magnitudo b 3 = -0,395, serta konstanta model persamaan empiris b 1 = -2,525. Penerapan (implementasi) nilai percepatan tanah model atenuasi dan rekaman akselerograf dalam bentuk pemetaan daerah terdampak akibat gempa bumi di wilayah Lombok. Kata kunci: koefisien geometrical spreading, koefisien magnitudo, konstanta, PGA
xiii
ABSTRACT Peak ground acceleration (PGA) is one of the parameters that were instrumental in determining the level of damage to the land and buildings that occur in the earth’s surface due to earthquake shocks. Absence of empirical formula developed to the Lombok as well as concerns about the impact of large and destructive earthquakes encourage researchers to develop a formulation of the attenuation function for peak ground acceleration. In this research, made the determination based on general model empirical formula from empirical formula Mc Guire (1997) using 522 data from akselerograph that corresponds to event of earthquakes around Lombok with the period 2018. Based on regression analysis, the obtained of geometrical spreading coefficients b 2= 0,712, magnitudo coefficient b 3= -0,395 , and also constants models PGA empirical b 1= -2,525. The implementation of the attenuation equation of the model and akselerograph in the mapping of affected areas due to earthquake in the Lombok. Keywords: geometrical spreading coefficient, magnitudo coefficient, constants, PGA
xiv
2
BAB I PENDAHULUAN 1.1. Latar Belakang Secara tektonik Lombok merupakan kawasan seismik aktif. Hal ini karena wilayah Lombok terletak di antara dua pembangkit gempa dari selatan dan utara. Dari selatan terdapat zona subduksi lempeng Indo Australia yang menunjam ke bawah Pulau Lombok, sedangkan dari utara terdapat struktur geologi sesar naik Flores (Flores back arc thusting). Sesar naik ini jalurnya memanjang dari laut Bali ke timur hingga Laut Flores. Hal ini yang menyebabkan Lombok rawan terjadi gempa bumi. Gempa bumi dengan kekuatan yang cukup besar akan sangat berdampak pada kerusakan bangunan. Salah satu faktor yang dapat menentukan besar kecilnya kerusakan tersebut adalah percepatan tanah maksimum atau Peak Ground Acceleration (PGA) (Edwiza, 2018). Setiap gempa akan menghasilkan satu nilai percepatan tanah pada suatu tempat. Nilai percepatan tanah yang dibutuhkan untuk menyesuaikan kekuatan bangunan yang akan dibangun adalah percepatan tanah maksimum. Pemetaan nilai percepatan tanah maksimum akan menjadi informasi penting dalam menunjang tata ruang dan wilayah di daerahdaerah rawan gempa (Kapojos dkk, 2015). Percepatan tanah merupakan salah satu parameter yang mempengaruhi tingkat bahaya suatu wilayah saat terjadi gempa bumi, yaitu tingkat kerusakan yang dapat ditimbulkan. Dengan mengetahui nilai pecepatan tanah pada suatu wilayah, maka dapat diketahui tempat mana yang berpotensi tinggi yang bisa terjadi kerusakan dan manakah yang berpotensi rendah terjadi kerusakan saat terjadi gempa bumi. Nilai percepatan tanah umumnya dapat di ukur langsung menggunakan instrumen yang disebut akselerograf. Selain menggunakan hasil pembacaan instrumen, nilai percepatan tanah dapat diketahui dengan menggunakan pendekatan persamaan empiris.
3
Pada wilayah Lombok telah terpasang beberapa akselerograf yang tersebar di beberapa wilayah seperti di Mataram dan Taliwang. Ketersediaan jaringan akselerograf yang terbatas dan belum adanya persamaan empiris yang dikembangkan untuk wilayah Lombok, serta kekhawatiran akan dampak gempa bumi besar dan merusak yang mendorong penelitian ini dilakukan dengan mengembangkan suatu persamaan untuk menentukan persamaan atenuasi percepatan tanah maksimum. Persamaan ini yang dapat menggambarkan intensitas gerakan tanah yang terjadi karena gempa dan dapat dijadikan pertimbangan dalam pengembangan yang lebih akurat di masa depan dalam hal analisis bahaya gempa dan kode bangunan, dan pemetaan wilayah rawan gempa bumi di wilayah penelitian. Penelitian ini akan menyempurnakan rumusan umum persamaan atenuasi Mc Guire (1997). Sehingga berdasarkan rumusan tersebut diperoleh nilai setiap koefisien dan konstanta yang sesuai dengan wilayah Lombok. Data yang digunakan adalah data hasil pembacaan akselerograf wilayah Lombok. 1.2. Perumusan Masalah Penentuan daerah rawan gempa bumi wilayah Lombok dilakukan berdasarkan nilai percepatan tanah maksimum. Nilai percepatan tanah tersebut diperoleh berdasarkan hasil pembacaan data accelerograph. Data hasil pembacaan ini akan diolah untuk memperoleh rumusan empiris yang sesuai dengan kondisi geologi Lombok. Berikut ini beberapa perumusan masalah terkait dengan penelitian ini, antara lain : 1. Bagaimana persamaan atenuasi percepatan tanah di wilayah Lombok ? 2. Bagaimana perbandingan nilai percepatan tanah hasil persamaan atenuasi dengan model dan hasil pembacaan rekaman akselerograf (Observation data)? 3. Bagaimana pemetaan percepatan getaran tanah dengan persamaan model dan data rekaman akselerograf? 1.3. Cakupan dan Batasan Masalah Cakupan masalah dalam penelitian ini, adalah:
4
1. Data rekaman akselerograf, data parameter gempa, dan site class
yang
digunakan bersumber dari Badan Meteorlogi Klimatologi dan Geofisika (BMKG). 2. Data yang digunakan dalam penelitian ini adalah data gempa bumi periode 1 januari sampai 31 desember 2018. 3. Pendekatan persamaan yang digunakan dalam pengolahan data, yaitu menggunakan persamaan Mc Guire 1997. 4. Magnitudo yang digunakan adalah magnitudo lokal ( M L ¿ . 5. Kedalaman yang digunakan kurang dari 50 km 6. Jarak hiposenter yang digunakan kurang dari 150 km. 7. Menggunakan metode analisis regresi non linier untuk mendapatkan konstanta dan koefisien dari persamaan percepatan tanah. Batasan masalah dalam penelitian ini adalah: 1. Komponen Z pada rekaman akselerograf tidak digunakan dalam penelitian ini. 2. Data gempa bumi yang digunakan belum terrelokasi. 3. Formula empiris diasumsikan hanya dipengaruhi oleh magnitudo dan jarak hiposenter. 1.4. Tujuan Penelitian Tujuan dalam penelitian ini adalah: 1. Menentukan persamaan atenuasi percepatan tanah di wilayah Lombok. 2. Menentukan perbandingan nilai percepatan tanah maksimum hasil perhitungan model dan hasil pembacaan sensor akselerograf (Observation Data). 3. Memetakan nilai percepatan tanah persamaan atenuasi model dan rekaman akselerograf. 1.5. Manfaat Penelitihan Hasil dari penelitian ini diharapkan dapat memberikan manfaat dan pertimbangan untuk pembaca. Beberapa manfaat yang diharapkan dari penelitian ini yaitu, mengetahui nilai percepatan gerakan tanah maksimum di wilayah Lombok sebagai upaya mitigasi bencana, mengetahui daerah yang rawan
5
mengalami kerusakan akibat terjadinya gempa bumi di wilayah Lombok, dan sebagai data dukung analisis resiko gempa bumi Probabilistik Seismic Hazard Analysis (PSHA), dan Deterministik Seismic Hazard Analysis (DSHA) dalam bidang teknik.
6
BAB II TINJAUAN PUSTAKA 2.1. Gelombang Seismik Gelombang seismik merupakan gelombang yang menjalar di dalam bumi yang disebabkan karena adanya deformasi struktur di bawah bumi akibat adanya tekanan ataupun tarikan karena sifat keelastisan bumi. Gelombang ini membawa energi kemudian menjalarkan ke segala arah di seluruh bagian bumi dan mampu dicatat oleh seismograf maupun akselerograf [CITATION Bor02 \l 1033 ]. Pada gelombang seismik terjadi osilasi partikel-partikel medium. Hal ini terjadi sebagai akibat interaksi antara gaya gangguan dan gaya elastik. Oleh karena itu, gelombang seismik disebut juga dengan gelombang elastik. Gelombang seismik pada keadaan tidak teredam dinyatakan dengan persamaan berikut 2
∇ ψ=
1 ∂2 ψ v2 ∂ t 2
(2.1)
dengan ∇=î
∂ ∂ ∂ +ĵ +¿ k ∂x ∂y ∂z
(2.2)
Dimana ψ merupakan fungsi gelombang, v merupakan kecepatan gelombang (m/s), t merupakan waktu (s). Pelepasan energi tegangan pada rekahan pada tepian lempeng tektonik merupakan penyebab utama dari aktivitas gempa bumi. Hal ini yang menyebabkan menjalarnya getaran pada bagian bumi dalam bentuk gelombang. Gelombang gempa bumi terdiri dari gelombang badan (body wave) dan gelombang permukaan (surface wave) [ CITATION Kra96 \l 1033 ]. Gelombang badan (body wave) yang merupakan gelombang menjalar melalui bagian dalam bumi dan biasa disebut dengan free wave karena dapat menjalar ke segala arah di dalam bumi. Gelombang badan terdiri atas gelombang
7
longitudinal (comperessional wave) dan gelombang transversal (shear wave). Gelombang P atau gelombang primer, merupakan gelombang longitudinal dimana gerakan partikel medium yang dilewati searah dengan arah penjalarannya. Karena bergerak lebih cepat dari pada gelombang S, gelombang P merupakan yang pertama tiba pada detektor gempa bumi (sehingga disebut “primer”). Gelombang S atau gelombang sekunder, merupakan gelombang transversal dimana gerakan partikelnya tegak lurus dengan arah penjalaran gelombangnya [ CITATION Pet13 \l 1033 ].
Pada penjalarannya, gelombang seismik dianggap sebagai sinar yang memenuhi hukum Snellius, prinsip Huygens, dan prinsip Fermat. Bumi sebagai medium penjalaran gelombang seismik dianggap berlapis-lapis sehinggga memiliki nilai kecepatan yang berbeda-beda. Pada bidang batas antar lapisan, gelombang seismik akan menjalar dengan kecepatan lapisan bawah bidang tersebut. Pertambahan nilai kecepatan sebanding dengan pertambahan kedalaman lapisan batuan. Hukum snellius menyatakan apabila gelombang jatuh pada bidang batas dengan nilai densitas yang berbeda, maka gelombang tersebut akan dibiaskan. Hal ini terjadi apabila sudut datang gelombang lebih kecil atau bernilai sama dengan sudut kritisnya. Jika nilai sudut datang lebih besar daripada sudut kritis maka gelombang tersebut akan dipantulkan. Secara umum, Hukum Snellius dinyatakan dengan persamaan berikut sin θ1 v 1 n1 = = sin θ2 v 2 n2
(2.3)
dimana θ1 merupakan sudut datang, θ2 merupakan sudut bias, v1 merupakan kecepatan cahaya sinar datang, v 2 merupakan kecepatan cahaya sinar bias, n1 merupakan indeks bias medium yang dilalui sinar datang, dan n2 merupakan indeks bias medium yang dilalui sinar bias.
8
Gambar 2.1 Refleksi Refraksi Gelombang Seismik [ CITATION Gad09 \l 1033 ]
Prinsip Huygens menyatakan bahwa setiap titik-titik penganggu yang berada di depan muka gelombang utama akan menjadi sumber bagi terbentuknya deretan gelombang yang baru. Titik-titik penganggu tersebut dapat berupa patahan, rekahan, pembajian, dan juga antiklin. Deretan gelombang baru yang terbentuk merupakan gelombang difraksi.
Gambar 2.2 Prinsip Huygens [ CITATION Gad09 \l 1033 ]
Prinsip fermat menyatakan bahwa jika sebuah gelombang merambat dari satu titik ke titik yang lain, maka gelombang tersebut akan memilih jejak yang tercepat.
9
Gam bar 2.3 Prinsip Fermat[ CITATION Gad09 \l 1033 ]
Gelombang seismik dalam penjalarannya akan mengalami refleksi, refraksi dan transmisi. Selain itu, gelombang seismik juga akan mengalami atenuasi atau peredaman oleh medium batuan bawah permukaan. Hal ini disebabkan karena adanya disipasi penjalarannya. Atenuasi gelombang seismik adalah proses pengurangan energi gelombang akibat penyerapan dan penyebaran. Pengaruh atenuasi terhadap sinyal seismik akan menyebabkan adanya penurunan amplitudo dan pelebaran sinyal, karena medium yang dilewati gelombang seismik berbeda-beda maka penyerapan frekuensi oleh medium tidak sama rata [CITATION Sut \l 1033 ].
2.2. Gempa Bumi Waktu dan lokasi terjadinya gempa bumi tidak dapat ditentukan secara spesifik, karena gempa bumi cenderung terjadi di tempat-tempat tertentu di lapisan kulit bumi. Lapisan kulit bumi terluar atau litosfer terdiri atas lempeng-lempeng tektonik yang kaku. Dan lempeng ini terapung di atas batuan yang relatif tidak kaku serta bergerak satu sama lain. Daerah pertemuan dua lempeng disebut sebagai plate margin atau batas lempeng, yang bisa berupa zona subduksi, pemekaran dasar samudra, atau pengangkatan, pelipatan di zona tumbukan. Umumnya gempa bumi terjadi di sekitar batas lempeng, yang membentuk jalur gempa bumi dunia, dan sekitar sesar[ CITATION Sun12 \l 1033 ].
10
Pergerakan antar lempeng akan terjadi diantara dua lempeng yang berada di sekitar sesar. Pergerakan yang akan terjadi yaitu keduanya bergerak saling menjauhi, saling mendekati, atau saling bergeser. Pergerakan yang terjadi berlangsung lambat dan tidak dapat dirasakan oleh manusia. Namun pergerakan ini dapat diukur, terjadi pergerakan sebesar 0-15 cm per tahun. Pergerakan yang terjadi juga tidak selalu lancar, terdapat keadaan gerakan lempeng yang terjadi tertahan. Hal ini menyebabkan antar lempeng saling mengunci sehingga terjadi pengumpulan energi yang berlangsung terus. Akumulasi dari energi ini akan terlepas pada saat antar lempeng tidak mampu lagi menahan stress tersebut. Akibatnya akan terjadi patahan secara spontan dan terlepasnya energi dalam bentuk getaran yang dikenal sebagai gempa bumi [ CITATION Sun12 \l 1033 ].
2.3. Tinjauan Seismistas Lombok Lombok merupakan salah satu pulau di gugusan kepulauan Nusa Tenggara. Secara geologi Pulau Lombok mempunyai batuan yang tergolong relatif muda, didominasi oleh batuan gunung api: breksi volkanik, lava dan batupasir [CITATION Agu052 \l 1033 ]. Batuan yang tertua di Pulau Lombok adalah batuan dari Formasi
Pengulung dan Kawangan, berumur Oligosen yang terbentuk dari kegiatan gunungapi bawah laut akibat adanya gejala tektonik. Gejala tektonik ini menyebabkan sesar normal dan sesar geser lurus yang berarah barat laut-tenggara [CITATION SAn94 \l 1033 ] . Selain sesar atau kelurusan, kekar juga banyak dijumpai
di Pulau Lombok (Agustawijaya, 2006). Formasi Pengulung dan Formasi Kawangan berada di bagian selatan Pulau Lombok, yang secara tektonik berada di bagian depan (fore arc). Di beberapa tempat batuan breksi volkanik dan batupasir dari Formasi Pengulung dan Formasi Kawangan ini diterobos oleh batuan beku basal. Sedangkan di bagian utara Pulau Lombok, yaitu di sekitar Gunungapi Rinjani, batuan terdiri dari batuan gunungapi
11
bersifat lepas dan berumur Kuarter. Lapisan batuan ini termasuk cukup tebal, dan menutupi hampir dua pertiga bagian Pulau Lombok. Bagian utara ini secara tektonik termasuk ke dalam bagian tengah (volcanic arc). Kaitannya dengan kegempaan, batuan-batuan ini memberi pengaruh kepada potensi bencana, terutama dalam merambatkan getaran gempa. USGS menyebutkan bahwa kerusakan akibat gempa adalah yang paling utama karena getaran. Intensitas getaran yang dialami oleh sebuah struktur bangunan adalah persamaan dari tiga faktor: kekuatan gempa, jarak terhadap sumber gempa, dan jenis material penyusun bumi. Makin kuat gempa, makin besar getaran, maka makin besar kerusakan. Makin dekat dengan patahan akibat gempa, makin besar getaran. Makin lunak jenis material, makin besar getaran. Posisi tektonik suatu wilayah menentukan jarak gempa. Makin dekat dengan pusat gempa, maka risiko gempa makin besar. Pulau Lombok berada persis di depan jalur tumbukan yang merupakan pusat gempa (Agustawijaya D, 2006). Jarak jalur tumbukan dengan titik terdekat Pulau Lombok adalah sekitar 45 km. Ketebalan lapisan batuan dan tanah di Pulau Lombok bervariasi dari sekitar 2 m hingga lebih dari 25 m. Jenis tanah yang dijumpai adalah didominasi oleh jenis tanah pasir lanauan, terutama untuk ketebalan 2 m. Sedangkan untuk ketebalan lapisan tanah mulai dari 5 m hingga lebih dari 25 m, jenis tanah yang dijumpai adalah tanah pasir lepas. Lapisan tanah pasir merupakan akuifer yang baik, sehingga mempunyai potensi likuifaksi jika terjadi gempa. Likuifaksi merupakan peristiwa mengubah tanah menjadi lumpur, seperti cairan. Sehingga tanah akan kehilangan kekuatannya dan menjadi lembek atau cair. Gempa bumi sebagai pemicu terjadinya likuifaksi bersumber dari aktifitas sesar naik yang terletak di sebelah selatan akibatnya akan ada gedung yang amblas, jempatan patah, jalanan terputus, dan lain sebagainya. Khusus di sekitar Gunungapi Rinjani, batuan yang berada di sini adalah batuan sedimen lepas hasil letusan gunungapi. Material batuan berupa kerikil,
12
lapili, dan bom, bersifat mudah lepas. Tanah penutup yang berupa pasir batuapung dengan ketebalan sekitar 2 m, bersifat lepas. Berdasarkan kondisi geologi di atas, tampak bahwa Pulau Lombok mempunyai potensi bencana yang cukup tinggi. Hal ini didasarkan pada kondisi tanah yang bersifat lepas, tebal, mempunyai akuifer, kondisi geologi struktur kekar dan sesar, serta dibagian utara pulau adalah merupakan kompleks gunungapi yang mempunyai batuan sedimen tidak kompak dan mudah lepas. Untuk peta geologi dapat dilihat pada Lampiran A.
2.4. Akselerograf Nilai percepatan tanah didapatkan dari analisis sinyal akselerograf, yang merekam tiga komponen getaran tanah, yaitu dua komponen horizontal (utara-selatan dan timur-barat) dan satu komponen vertikal. Untuk kepentingan persamaan atenuasi, data yang digunakan yaitu hasil rekaman dari kedua komponen horizontal dengan mengaplikasikan resultan percepatan tanah sedangkan komponen vertikal jarang digunakan karena struktur bangunan sudah didesain mampu menahan gaya axial tekan [ CITATION Oct16 \l 1033 ] . Pengukuran akselerograf dilakukan dengan menghitung amplitudo maksimum komponen sinyal (zero to peak) dari rekaman akselerograf. Berkembangnya teknologi semikonduktor memungkinkan untuk membuat suatu sistem sensor getaran. Dengan ditemukannya sistem piezoresistive maka sensor akselerasi menjadi semakin praktis. Efek piezoresistif adalah perubahan resistivitas listrik semikonduktor atau logam saat regangan mekanis diterapkan. Piezoresistivitas adalah suatu kemampuan yang dimiliki sebagian kristal maupun bahan-bahan tertentu lainya yang dapat menghasilkan arus listrik jika mendapatkan perlakuan tekanan. Sistem kerja ini berdasarkan gaya yang bekerja pada salah satu permukaan semikonduktor. Berdasarkan hukum Newton, gaya sebanding dengan massa dikali percepatan[ CITATION Dad14 \l 1033 ].
13
F=m∙ a
(2.4)
dimana F adalah Gaya dalam (N), m adalah massa (kg), dan a adalah percepatan ( m/s 2 ¿ ¿ Piezoelectric tranducer adalah perangkat yang menggunakan efek piezoelectric
untuk
mengukur
tekanan,
percepatan,
regangan
dengan
mengubahnya ke sinyal istrik. Efek piezoelectric adalah muatan yang terakumulasi pada bahan padat tertentu seperti kristal dan keramik akibat dikenai tekanan
(pressure). Pada
prinsipnya,
efek
piezoelectric
diperoleh
dari
ketidakseimbangan distribusi arus listrik pada bahan piezoelectric yang disebabkan oleh terjadinya tegangan (stress) pada bahan tersebut dan selanjutnya mengakibatkan terjadinya regangan. Bila kedua permukaan bahan tersebut dilapisi dengan bahan logam dan lempengan kecil tembaga ditempelkan padanya, maka perubahan arus listrik ini akan mengakibatkan terjadinya sinyal tegangan listrik pada lempeng tembaga tersebut. Efek ini dapat dibalik yaitu bila tegangan listrik dikenakan pada bahan tersebut maka akan terjadi regangan. Hal ini dapat di ilustrasikan sebagai plat paralel kapasitor. Muatan total yang di hasilkan oleh tranduser sebanding dengan gaya yang bekerja, yang ditunjukan oleh persamaan dibawah ini [ CITATION Boo04 \l 1033 ]. q=k ∙ F
(2.5)
dimana q adalah muatan listrik (coulomb), k adalah konstanta Dari persamaan 2.4 dan persamaan 2.5 didapat q =m ∙ a k
(2.6)
Nilai perubahan muatan listrik yang terjadi pada piezoelektric akan sebanding dengan nilai percepatan (q sebanding a). Oleh karena itu maka keluaran dari piezoelectric sebanding dengan nilai percepatan yang diukur (dalam hal ini adalah percepatan permukaan tanah).
14
Cara kerja accelerograph yaitu getaran seismik yang merupakan besaran fisis akan ditangkap oleh sensor dalam hal ini akselerograf yang berfungsi merubah besaran fisis menjadi besaran elektik. Keluaran dari akselerograf sudah berupa tegangan analog. Tegangan analog ini kemudian dirubah oleh ADC menjadi count yang sudah berbentuk nilai digital [ CITATION Sun18 \l 1033 ]. Untuk mendukung monitoring gempa bumi di wilayah Lombok telah dipasang seismograf dan akselerograf. Alat ini dapat menerima penjalaran gelombang yang terjadi akibat lempeng tektonik. Dari data yang diperoleh seismograf dapat digunakan unuk menentukan parameter gempa bumi seperti waktu, episenter, magnitudo, dan kedalaman. Sedangkan akselerograf mempunyai tingkat kemampuan yang lebih baik karena sudah mengukur percepatan gerakan tanah (strong motion).
Gambar 2.4 Jaringan akselerograf di wilayah Lombok
Jaringan seismograf dan akselerograf
yang dipasang di wilayah Indonesia
memiliki dua jenis jaringan, yaitu Libra (Indonesia) dan Jaringan GFZ (Jerman).
15
Untuk jaringan akselerograf di wilayah Lombok sendiri telah terpasang beberapa akselerograf yaitu: Tabel 2.1 Jaringan Akselerograf di Wilayah Lombok Koordinat Kode
Stasiun Lintang
Elevasi
Jaringan
Bujur
MASE
STA MET Selampang Mataram
-8,562
116.167
192
Libra
TWSI
Taliwang
-8,738
116,882
124
Libra
2.5. Percepatan Tanah Percepatan merupakan parameter yang menyatakan perubahan kecepatan mulai dari keadaan diam sampai pada kecepatan tertentu. Nilai percepatan terbagi menjadi dua bagian yaitu percepatan tanah maksimum dan percepatan tanah sesaat. Percepatan tanah merupakan parameter yang perlu dikaji pada setiap terjadinya gempa bumi untuk dipetakan agar bisa memberikan informasi tentang efek paling parah yang pernah dialami suatu lokasi. Percepatan tanah maksimum atau Peak Ground Acceleration (PGA) adalah nilai terbesar percepatan tanah pada suatu tempat akibat getaran gempa bumi dalam periode waktu tertentu. Nilai percepatan tanah dipengaruhi oleh nilai kekuatan gempa bumi (magnitudo), jarak, kedalaman sumber gempa bumi, dan geologi permukaan suatu daerah. Semakin besar nilai magnitudo, semakin besar nilai percepatan tanah. Semakin dekat jarak sumber gempa bumi dengan daerah penelitian, semakin besar juga nilai percepatan tanah yang diperoleh. Pengukuran percepatan tanah dengan cara empiris dapat dilakukan dengan pendekatan dari beberapa rumus yang diturunkan dari magnitudo gempa bumi atau data intensitas. Untuk melakukan perhitungan nilai percepatan tanah maksimum, terdapat hubungan antara intensitas getaran tanah, magnitudo dan jarak suatu daerah dari sumber gempa bumi yang dapat digambarkan dengan menggunakan fungsi atenuasi [ CITATION Paw12 \l 1033 ].
16
2.6. Faktor-faktor yang mempengaruhi persamaan atenuasi gerakan tanah Atenuasi gerakan tanah (ground motion attenuation) adalah proses atau persamaan suatu gerakan tanah akibat gempa akan mengecil pada jarak yang semakin jauh dengan sumber gempa. Secara matematis dapat dijelaskan bahwa atenuasi gerakan tanah adalah suatu hubungan antara parameter gempa (percepatan, kecepatan, simpangan, intensitas gempa, atau magnitudo gempa) dengan jarak ke lokasi pencatat gempa (jarak episenter, jarak hiposenter, atau jarak terdekat). Misalnya hubungan antara percepatan tanah dengan jarak episenter untuk setiap magnitudo gempa yang berbeda. Dengan persamaan atenuasi yang sudah diketahui maka gerakan tanah ataupun intensitas gempa di suatu tempat relatif terhadap sumber gempa dapat diprediksi.Terdapat beberapa faktor atau parameter yang secara dominan maupun kurang dominan akan mempengaruhi atenuasi gerakan tanah. Parameter yang dimaksud adalah sebagai berikut: 2.6.1. Magnitudo gempa (earthquake magnitude) Magnitudo gempa merupakan suatu ukuran yang merepresentasikan kekuatan gempa yang terjadi. Magnitudo gempa umumnya merupakan skala logaritmik kekuatan gempa berdasarkan pada pengukuran amplitudo maksimum gelombang gempa. Amplitudo tersebut diperoleh dari hasil pencatatan peralatan rekam gelombang gempa yang disebut seismograf. Magnitudo gempa merupakan skala ukuran kekuatan gempa yang dicatat pada pusat gempa, sehingga sebuah gempa hanya mempunyai satu ukuran atau besaran magnitudo. Kekuatan guncangan tanah yang diakibatkan oleh suatu gempa sangat bergantung pada besarnya magnitudo gempa yang terjadi. Semakin besar gempa, semakin besar pula guncangan tanah pada suatu lokasi demikian pula sebaliknya. Magnitudo gempa bumi terdiri dari 4 jenis, yaitu magnitudo lokal, magnitudo permukaan, magnitudo badan, dan magnitudo momen. Dalam persamaan percepatan tanah, magnitudo yang dibutuhkan adalah magnitudo momen karena magnitudo momen merupakan jenis magnitudo yang tidak
17
mengalami saturasi atau kejenuhan saat dilakukan perhitungan. Perbedaan dan definisi dari setiap magnitudo, yaitu:
Tabel 2.2 Jenis-jenis Magnitudo
No Jenis Magnitudo 1 Magnitudo lokal ( M L )
Definisi
Penggunaan
Magnitudo gempa lokal T s± 1 Untuk gempa bumi detik. Saturasi pada magnitudo dengan M L = 3-7 sekitar 7.
2
Magnitudo permukaan Ms ¿
Magnitudo gempa bumi untuk ( gelombang permukaan teleseismik T −wave ± 20 detik
Untuk gempa bumi dengan M s= 5-7,5
3
Magnitudo badan (M b ¿
Hubungan magnitudo untuk Untuk gempa bumi gelombang badan teleseismik dengan M b = 5-7 (P,PP,S) T −wave1-3 detik
4
Magnitudo momen (M w ¿
Dihitung berdasarkan yang dilepaskan sumber
energi Untuk gempa bumi dengan M w ˃ 7
Tabel 2.3 Korelasi konversi magnitudo untuk wilayah Indonesia Korelasi Konversi
Jumlah
Range Data
Data 2
M w =0,143 M s −1,051 M s +7,285 M w =0,114 m b2−0,556 m b +5,560 M w =0,787 M E +1,537 M s m b=0,125 M L2−0,389 M L +3,513 M L =0,717 M D +1,003
Kesesuaian (R2)
(Events) 3173
4,5≤ M s ≤8,6
93,9%
978
4,9≤m b ≤8,2
72,0%
154 722
5,2≤ M E ≤ 7,3 3,0≤ M L ≤ 6,2
71,2% 56,1%
384
3,0 ≤ M D
29,1%
(magnitudo durasi ≤ 5,8) [ CITATION Asr10 \l 1033 ] memberikan hubungan antara jenis-jenis magnitudo
gempa bumi dengan menghimpun banyak data gempa bumi di Indonesia. Setelah melakukan regresi, didapat persamaan korelasi konversi magnitudo untuk wilayah Indonesia disajikan pada Tabel 4.4. Berdasarkan penilitian yang dilakukan oleh
18
Tim Pusat Studi Gempa Nasional (2017), magnitudo lokal ( M L ) tidak perlu dikonversi menjadi M W karena hasil plot pada Gambar 2.7 menunjukkan hasil yang hampir sebanding di antara keduanya dan dianggap magnitudo lokal ( M L ) dapat merepresentasikan magnitudo momen
9 8 7 6 5 y = 0,7473x + 1,0651
4 3 2 2
3
4
5
6
7
8
Gambar 2.5 Korelasi antara M L dan M w dari katalog ISC, BMKG, dan USGS untuk gempa dari tahu 1907-2013 sebanyak 2.282 gempa dengan magnitudo M L 2,0 hingga 7,5 dan Mw dari magnitudo 3,0 hingga 7,8 [ CITATION Irs10 \l 1033 ]
2.6.2. Jarak ke site Tingkat guncangan tanah akibat gempa di suatu lokasi (site) sangat bergantung pada jarak lokasi tersebut terhadap sumber gempa. Semakin jauh jarak site dengan sumber gempa, semakin kecil site terhadap sumber gempa, semakin besar goncangan yang dirasakan pada lokasi tersebut. 2.6.3. Pengaruh mekanisme sumber gempa (source mechanism) Mekanisme sumber atau solusi bidang sesar adalah pernyataan tentang orientasi bidang sesar dan gerakan relatif (slip) yang terjadi pada sistem koordinat geografi [ CITATION Pur17 \l 1033 ] . Mekanisme sumber digunakan untuk memahami
orientasi geometri sesar secara aktual, mengidentifikasi jenis sesar, dan memperkirakan kondisi tegangan (stress) pada sumber gempa bumi. Solusi mekanisme sumber ditentukan berdasarkan pola radiasi gelombang seismik yang
19
dipancarkan oleh sumber gempa bumi dengan melihat arah gerakan sesar dan arah bidang sesar. 2.6.4. Pengaruh kondisi situs (lokal site condition) Suatu energi gempa yang datang dari tempat dan jarak yang sama, kemudian direkam diatas tanah rock dan tanah endapan akan mempunyai karakter rekaman (percepatan
tanah,
durasi
gempa,
dan
kandungan
frekuensi)
yang
berbeda[ CITATION Idr90 \l 1033 ]. Dengan demikian kondisi tanah situs (site condition) di bawah seimograf merupakan parameter penting yang diperhitungkan dalam menentukan persamaan atenuasi. Contoh pengaruh site condition terhadap atenuasi misalnya adalah pada Gambar 2.9. Pada gambar tampak bahwa pada jarak dekat (kurang dari 5 km) PHA untuk rock site lebih besar daripada soil site, sementara itu berbanding terbalik untuk pada jarak lebih dari 5 km. Hal ini berarti bahwa pada soil site, atenuasi Gerakan tanah akan berlangsung lebih lambat dibandingkan di rock site . Daya redam energi soil site lebih kecil daripada rock
Jarak terdekat (km)
Jarak terdekat (km)
site merupakan akibat dari hal tersebut. Batas tersebut sedikit bergeser atau membesar pada magnitudo gempa yang semakin besar [ CITATION Abr97 \l 1033 ]. Gambar 2.6 Pengaruh kondisi tanah terhadap PHA
Dalam Rudyanto (2013), Aki dan Richards menyatakan bahwa amplitudo ground motion pada permukaan tanah sangat dipengaruhi oleh kerapatan dan
20
kecepatan shear wave pada medium dekat permukaan. Borcherdt dan Glassmoyer, 1994 mengatakan bahwa pendekatan praktis dalam memperkirakan profil kecepatan shear wave ini dapat diperoleh dengan menggunakan kecepatan pada kedalaman kurang dari 30 m di stasiun atau disebut dengan Vs30 [ CITATION Rud13 \l 1033 ]. Selanjutnya tipe tanah dapat dikategorikan berdasarkan klasifikasi dari
National Hazard Reduction Program (NEHRP) seperti pada tabel dibawah ini. Tabel 2.4 Klasifikasi jenis tanah berdasarkan nilai Vs30 menurut NEHRP
Jenis Tanah A B C D E 2.7.
Deskripsi Hard Rock Rock Very dense soil and soft rock Stiff soil Soft clay soil
Vs30 ≥1500 m/s 760 to1500 m/s 360 to 760 m/s 180 to 360 m/s ≤ 180 m/s
Ground Motion Prediction Equation (GMPE)
GMPE disebut juga persamaan atenuasi, yaitu persamaan matematis yang berkaitan dengan parameter getaran tanah kuat terhadap karakteristik gempa bumi, media penjalaran, dan kondisi lokal [ CITATION Sri12 \l 1033 ]. GMPE berguna untuk analisis bahaya seismik terkait dengan rancangan bangunan tahan gempa bumi dan sebagai pendekatan jika suatu wilayah tidak memiliki sensor yang memadai. GMPE yang digunakan dalam penelitian ini adalah GMPE Mc Guire (1997).
Mc Guire (1997) melakukan penelitian resiko seismik dari ruangan
bawah tanah suatu bangunan atau dari lapangan bebas yang terjadi di California Selatan. Menggunakan kurang lebih tujuh catatan dari gempa bumi yang sama dan tidak lebih dari sembilan catatan dari satu site untuk meminimalisasi terlalu rendahnya perkiraan variasi yang dihitung. Rekaman yang digunakan dengan jarak yang besar, kedalaman yang dangkal, dan magnitudo besar. Penelitian yang dilakukan juga menggunakan 11 rekaman dengan kategori site rock untuk batuan sedimen atau tanah dengan ketebalan kurang dari 10 meter, dan 59 rekaman dengan kategori site soil untuk batuan alluvium atau tanah lunak dengan ketebalan
21
lebih dari 10 meter. Persamaan atenuasi yang digunakan oleh [ CITATION McG77 \l 1033 ] sebagai berikut:
log ( Y )=b1 +b2 ∙ M −b3 log( Rh +25)
(2.7)
Persamaan 2.7 tersebut dapat disederhanakan menjadi, Y =b1 ∙ 10b ∙ M (Rh +25)−b 2
(2.8)
3
Dimana Y adalah nilai percepatan tanah maksimum (gal), M adalah magnitudo, Rh adalah jarak hiposenter (km), b 1 , b2 , dan b 3 adalah suatu koefisien.
Pengembangan persamaan atenuasi dilakukan berdasarkan analisis regresi secara statistik pada rekaman gerakan tanah atau akselerogram. Hal ini juga berlaku di Indonesia, sampai saat ini belum tersedia persamaan atenuasi yang spesifik untuk wilayah Indonesia karena keterbatasan rekaman gerakan tanah untuk wilayah Indonesia [ CITATION Irs07 \l 1033 ]. Oleh karena itu, studi terkait bahaya gempa bumi di Indonesia menggunakan persamaan atenuasi yang dikembangkan untuk wilayah lain di dunia, karena mempunyai kemiripan karakter gempa bumi dan kondisi tektonik maka pemakaian persamaan ini bisa diaplikasikan untuk wilayah Indonesia. Dalam[ CITATION Dou11 \l 1033 ] beberapa persamaan empiris yang digunakan dalam perhitungan percepatan tanah maksimum di antaranya adalah persamaan Fukushima dan Tanaka (1990), Esteva, dan Setiawan. Bentuk persamaan Fukushima dan Tanaka (1990) diterapkan di wilayah Jepang, dengan bentuk persamaan sebagai berikut: log(Y )=b1 ∙ M S −log ( Rh +b 2 10b ∙ M )−b3 ∙ Rh +1,30 1
s
(2.9)
Dimana Y adalah nilai percepatan tanah maksimum (gal), M S adalah magnitudo permukaan, Rh adalah jarak hiposenter (km), b 1= 0,41, b 2=0,032 ,dan b 3= 0,0034.
22
Rekaman yang digunakan data Jepang dengan kedalaman < 30 km, M JMA ˃ 5,0 dan PGA ≥ 0,1 m/s 2. Bentuk persamaan Esteva (1970) adalah sebagai berikut b 1( expb ∙ M ) 2
Y=
( Rh +40 )
b3
(2.10)
Dimana Y adalah nilai percepatan tanah maksimum (gal), M adalah magnitudo, Rh adalah jarak hiposenter (km), b 1= 5600, b 2=0,8 , b 3= 25, dan b 4=2. Rekaman yang digunakan dari tanah yang sebanding dengan tanah liat atau batuan konglomerat yang padat. Bentuk persamaan Lin dan Wu yang diterapkan di wilayah India. Persamaan formulanya adalah sebagai berikut: log 10 (Y )=b 1 · log 10 ( R h ) +b2 · M b + b3
(2.11)
Dimana Y adalah nilai percepatan tanah maksimum (gal), M badalah magnitudo body, Rh adalah jarak hiposenter (km), b 1= -2,0663, b 2=0,9019 , dan b 3= 0,1091. 2.8.
Analisis Regresi Metode regresi dapat dikelompokkan menjadi dua macam, yaitu regresi linier
dan regresi non linier. Metode regresi linier digunakan dalam penelitian apabila pola hubungan antara variabel dependen dan variabel independen adalah linier. Sedangkan metode regresi non linier digunakan apabila pola hubungan antara variabel dependen dan variabel independen adalah tidak linier dari satu atau lebih parameter. Hubungan antar variabel dalam regresi non linier membentuk fungsi tertentu [ CITATION Sar06 \l 1033 ]. Permasalahan regresi secara umum dapat diformulasikan dengan melibatkan variabel atau parameter yang ada untuk dinyatakan sebagai notasi vektor atau matriks yang mempresentasikan variabel dengan banyak komponen atau elemen. Jika data (d) dan model (m) masing-masing dinyatakan dengan vektor: d=[d 1 , d 2, d 3 , … , dN ]T
(2.12)
m= [ m 1, m 2, m3 ,… , mM ] T
(2.13)
Hubungan antara data dan model menjadi:
23
G ( m )=d
(2.14)
Dimana G merupakan fungsi umum pemodelan ke depan (forward modeling) yang memetakan model menjadi besaran dalam “domain” data. Dengan kata lain, fungsi G memungkinkan kita memprediksi data suatu model m. Elemen data yang bervariasi direpresentasikan sebagai d. Penyusunan persamaan tersebut untuk mendapatkan perkiraan model (m) sebagai solusi persamaan inversi menjadi sebagai berikut :
[
G11 G12 … … G1 M m1 d1 G21 G22 … …G2 M m2 d2 ∗¿ . . = . . . . GN 1 GN 2 … … G NM mM d3
(2.15)
m=[GT G]−1 GT d
(2.16)
] [ ][]
Matriks G T G adalah matriks bujur sangkar (M x M ) sesuai dengan parameter model yang dicari. Jika bukan merupakan matriks singular maka inversi matriks tersebut dapat dihitung menggunakan teknik inversi matriks yang umum seperti eliminasi Gauss-Jordan, dekomposisi LU, dan sebagainya. Teknik inversi yang lebih stabil untuk matriks yang mendekati singular adalah teknik dekomposisi nilai singular (Singular Value Decomposition). Elemen matriks parameter tersebut adalah nilai percepatan tanah maksimum sebagai variabel terikat, kemudian, magnitudo dan jarak hiposenter sebagai variabel bebas.
2.9. Analisis Korelasi Koefisien korelasi merupakan pengukuran statistik kovarian atau asosiasi antara dua variabel. Besarnya koefisien korelasi berkisar +1 s/d -1. Koefisien korelasi menunjukan kekuatan (strength) hubungan linier dan arah hubungan dua variabel acak. Jika koefisien korelasi positif, maka kedua variabel mempunyai hubungan searah. Jika nilai variabel X tinggi, maka nilai variabel Y akan tinggi pula. Sebaliknya, jika koefisien korelasi negatif, maka kedua variabel mempunyai
24
hubungan terbalik. Jika nilai variabel X tinggi, maka nilai variabel Y akan menjadi rendah dan berlaku sebaliknya [ CITATION Sar06 \l 1033 ]. Untuk memudahkan melakukan interpretasi tentang kekuatan hubungan antara kedua variabel, maka dilakukan sebagai berikut: R square=R2 Tabel 2.5 Interpretasi Koefisien Korelasi
Nilai Korelasi
Tingkat Hubungan
0,00 – 0,29
Sangat Rendah
0,30 – 0,59
Rendah
0,60 – 0,79
Kuat
0,80 – 1,00
Sangat Kuat
2.10.
Intensitas Seismik
Intensitas seismik adalah skala yang dihitung berdasarkan kerusakan wilayah atau bangunan di dekat terjadinya gempa bumi. Skala intensitas yang digunakan di Indonesia adalah skala MMI (Modified Mercalli Intensity). Berikut hubungan antara nilai percepatan tanah dan intensitas seismik gempa bumi.
Tabel 2.6 Intensitas Seismik
Skala BMKG
Percepatan Tanah (gal)
Skala MMI
I
< 2,9
I-II
Dirasakan oleh beberapa orang tetapi terekam oleh alat.
II
2,9 – 88
III-V
Dirasakan oleh orang banyak tetapi tidak menimbulkan kerusakan.
III
88 – 167
VI
Mengalami kerusakan ringan seperti
Deskripsi
25
genteng IV
168 – 564
VII-VIII
V
564
IX-XII
bergeser
kebawah
dan
Terjadi kerusakan sedang seperti retakan yang terjadi pada dinding bangunan dan sebagainya. Struktur bangunan kerusakan berat.
mengalami
[ CITATION BMK17 \l 1033 ]
Intensitas berbeda dengan magnitudo, karena intensitas merupakan hasil dari pengamatan visual pada suatu tempat tertentu, sedangkan magnitudo merupakan hasil pengamatan alat instrumental menggunakan seismograf. Pada suatu kejadian gempa bumi besarnya intensitas pada suatu tempat yang berbeda dapat sama atau berlainan, sedangkan besarnya magnitudo selalu sama walaupun dicatat atau dirasakan di tempat yang berbeda.
26
BAB III METODE PENELITIAN 3.1 Tempat dan Waktu Penelitian dilakukan di Badan Meteorologi dan Klimatologi Geofisika (BMKG) Pusat, Jakarta tanggal 24 Juni 2019 sampai dengan 10 Oktober2019. 3.2 Alat dan Bahan Alat yang digunakan pada penelitian ini, yaitu perangkat lunak SPSS yang digunakan untuk mencari konstanta dan koefisien yang di butuhkan dalam persamaan empiris, perangkat lunak SIGMA dan Shakemap yang digunakan untuk mendapatkan nilai percepatan tanah hasil rekaman akselerograf, perangkat lunak ArcMap 10.4 digunakan untuk pemetaan nilai percepatan tanah di daerah penelitian. Bahan-bahan yang digunakan dalam penelitian ini adalah data historis gempa bumi periode 2018 dengan magnitudo kurang dari 3,5 M S . Parameter meliputi lokasi gempa bumi (koordinat episenter), origin time, magnitudo dan jenis magnitudo serta kedalaman pusat gempa bumi, dan data PGA hasil pembacaan akselerograf MASE dan TWSI yang berada di Lombok dengan menggunakan data dari BMKG tahun 2018. Nilai PGA dalam penelitian ini dihitung dari resultan komponen gerak horizontal (horizontal motion). 3.3 Tahapan Penelitian tentang persamaan atenuasi yang paling sesuai untuk wilayah Lombok dilakukan melalui beberapa tahapan. Tahapan- tahapan tersebut akan dijelaskan sebagai berikut: 3.3.1
Tahapan menghitung jarak stasiun ke sumber gempa ( Rhiposenter ¿
Dalam perhitungan persamaan atenuasi, jarak dari stasiun ke sumber gempa menggunakan jarak hiposenter. Penggunaan jarak hiposenter dalam menghasilkan persamaan atenuasi di wilayah Lombok karena hiposenter merupakan model sederhana tentang jarak terdekat antara sumber dan stasiun pencatat serta hiposeter merupakan sumber dari energi radiasi gempa bumi. 22
23
Gambar 3.1 Ilustrasi jarak stasiun ke sumber gempa
Kemudian jarak hypocenter bisa dihitung dengan rumus phytagoras:
√
2
2
2
(3.1)
Rh = ( X h− X s ) + ( Y h −Y s ) + ( Z h−Z s )
Rh =
2
2
√[ (( R−h) ∙ cos ϕ ∙ cos θ )−( R ∙ cos ϕ −cos θ ) ] + ¿[ (( R−h ) ∙ cos ϕ ∙ sin θ )−( R ∙cos ϕ −sin θ ) ] + ¿[ (( h
h
s
s
h
h
s
(3.2) X s=R ∙ cos ϕh ∙ cos θh
(3.3)
Y s =R ∙cos ϕs ∙ cos θs
(3.4)
Z s=R ∙ sin ϕs
(3.5)
Dimana Rh dalah jarak hiposenter (km), h merupakan kedalaman (km), ϕ h adalah lintang hiposenter¸ϕ s adalah lintang stasiun, θh adalah bujur hiposenter, θ s adalah bujur stasiun, dan R adalah jari-jari bumi (6371 km). 3.3.2
Analisi regresi untuk menghasilkan persamaan atenuasi
Persamaan atenuasi menghubungkan parameter gerakan tanah (ground motion) ke bermacam variabel bebas (independen) seperti magnitudo gempa bumi, jarak dari sumber ke site, dan kondisi lokal. Variabel terikat (dependen) merupakan variabel yang diprediksi yaitu nilai puncak percepatan tanah (PGA). Hubungan antara puncak percepatan tanah horizontal (PGA) dengan parameter prediksi dalam
s
24
bentuk matematis bukanlah hal yang mudah untuk diselesaikan dengan hubungan non linier. Untuk membangun persamaan atenuasi secara konversional digunakan analisis regresi. Melalui metode regresi ini dapat ditentukan seberapa baik kurva regresi dari persamaan atenuasi sesuai dengan data pengamatan yang terlihat dari plot residualnya (Boore dan Atkinson, 1982). Kebanyakan bentuk persamaan menggunakan logaritma dari nilai percepatan tanah yang diprediksi, karena penggunaan logaritma akan menghasilkan residual yang sesuai distribusi lognormal dengan variasi yang bebas dari jarak pengukuran (Campbell,1981). Suatu model persamaan regresi secara fisis dapat menghitung saturasi dari amplitudo percepatan tanah di wilayah dekat sumber yang terkait. Model juga setidaknya memiliki karakteristik simple point source termasuk geometric spreading dalam jarak menengah sampai jauh. Melakukan regresi non linier dengan masukan data parameter yang ada dari hasil sortir dengan nilai residual untuk mendapatkan nilai konstanta dan koefisien regresi dengan metode regresi pada persamaan atenuasi Mc Guire (1997) yaitu: log ( Y )=b1 +b2 ∙ M −b3 log( Rh +25)
(3.6)
Dimana Y adalah nilai percepatan tanah maksimum (gal), M L adalah magnitudo lokal, Rh adalah jarak hiposenter (km), b 1 , b2 , dan b 3 adalah suatu koefisien. 3.3.3
Validasi
Pada tahap validasi ini akan dilakukan uji coba terhadap persamaan atenuasi yang telah diperoleh dari penelitian. Validasi dilakukan dengan membandingkan hasil perhitungan persamaan atenuasi tersebut dengan data pengukuran percepatan tanah dari gempa bumi pada tahun 2018 yang tercatat oleh sensor akselerograf di wilayah Lombok. Perbandingan dilakukan dengan cara membuat kurva atenuasi dari model persamaan atenuasi yang terpilih terhadap jarak hiposenter dengan nilai magnitudo tertentu. Nilai percepatan tanah hasil pengukuran kemudian diplot pada kurva atenuasi. Jika plot percepatan tanah hasil pengukuran rekaman akselerograf berada di antara kurva model dan ambang batas standar deviasi,
25
maka model dianggap sudah sesuai untuk wilayah Lombok pada nilai parameter gempa tersebut. 3.3.4 Penerapan persamaan atenuasi dalam pemetaan wilayah terdampak akibat gempa bumi Pemetaan menggunakan rekaman akselerograf dengan cara menentukan jarak hiposenter dengan setiap kejadian gempa. Kemudian nilai percepatan tanah maksimum pada tiap kejadian gempa yang direkam oleh akselerograf dipilih nilai tertinggi pada tiap sensor akselerograf, selanjutnya diplot menggunakan program ArcGis 10.4 sehingga didapatkan peta percepatan tanah maksimum akibat gempa bumi untuk wilayah Lombok. Setelah diperoleh modifikasi persamaan empiris percepatan tanah, maka digunakan untuk mengidentifikasi nilai percepatan tanah untuk kasus event gempa. Pemetaan dilakukan dengan membagi daerah penelitian menjadi grid dengan ukuran 0,5°×0,5° dan menentukan jarak hiposenter dengan setiap titik grid. Kemudian nilai percepatan tanah maksimum pada tiap grid dihitung dengan menggunakan fungsi atenuasi model atenuasi pada setiap kejadian gempa bumi yang berpusat di daerah penelitian. Dari hasil perhitungan kemudian dipilih nilai tertinggi pada tiap grid dan selanjutnya diplot menggunakan program ArcGis 10.4 sehingga
26
didapatkan peta percepatan tanah maksimum akibat gempa bumi untuk wilayah Lombok. Gambar 3.2 Titik grid daerah penelitian
3.4 Diagram Alir Penelitian Adapun diagram alir penelitian yaitu:
27
Mulai
Persiapan alat dan bahan
Sortir pemilihan data PGA rekaman akselerograf yang sesuai dengan katalog parameter gempabumi Data PGA dan data gempabumi yang bersesuaian
Percepatan tanah observasi, Magnitudo
Jarak episenter, kedalaman Menghitung jarak hiposenter terhadap sensor akselerograf Jarak hiposenter
Analisis regresi untuk mendapatkan b1, b2, b3 untuk fungsi atenuasi PGA model Persamaan empiris PGA model berdasarkan data akselerograf
28
Menentukan nilai PGA menggunakan formula empiris model
PGA
Validasi nilai PGA menggunakan persamaan empiris model dan rekaman akselerograf Grafik Perbandingan PGA
Pemetaan nilai PGA daerah terdampak akibat gempabumi dengan persamaan empiris model dan rekaman akselerograf Peta PGA model dan PGA rekaman akselerograf Selesai Gambar 3.3 Diagram alir penelitian
BAB IV HASIL DAN PEMBAHASAN
4.1 Menentukan persamaan empiris percepatan getaran tanah Untuk menentukan persamaan empiris percepatan tanah, dalam penelitian ini digunakan data observasi akselerograf BMKG di wilayah Lombok selama tahun 2018 beserta dengan parameter gempa yang diperoleh dari katalog gempa BMKG tahun 2018. Tahapan dalam menentukan formula atau rumus empiris percepatan getaran tanah adalah melalui pemilihan data katalog gempa bumi, perhitungan nilai percepatan getaran tanah (PGA) hasil observasi peralatan akselerograf, menghitung nilai jarak hiposenter terhadap senso akselerograf, dan penentuan rumus empiris percepatan getaran tanah menggunakan analisis regresi. 4.1.1. Pemilihan data katalog gempa bumi di wilayah Lombok dan sekitarnya Langkah pertama dalam menentukan rumus empiris percepatan tanah adalah dengan mengunduh database kejadian gempa bumi di repository gempa bumi berdasarkan klasifikasi wilayah penelitian 8,03 - 9,27 °LS dan 115.75 – 117.21 °BT dengan magnitudo lebih dari 3,5 M L selama kurun waktu tahun 2018. Parameter yang diperlukan antara lain magnitudo, dan jarak hiposenter terhadap sensor akselerograf dari katalog Badan Meteorlogi Klimatologi dan Geofisika (BMKG) dari tahun 2018 (data lampiran B). Setelah itu sebaran gempa bumi dari database tersebut dipetakan menggunakan program ArcGis 10.4. dalam pemetaan tersebut ditampilkan koordinat stasiun akselerograf di dalam lingkup wilayah penelitiam berdasarkan dari database bidang Seismologi Teknik BMKG.
28
29
Gambar 4.1 Sebaran Gempa bumi Lombok tahun 2018
4.1.2. Pemilihan data percepatan getaran tanah hasil analisis rekaman Akselerograf Pemilihan data rekaman akselerograf untuk mendapatkan nilai percepatan tanah, dengan mengumpulkan katalog gempa bumi yang berisikan parameter gempa bumi. Data akselerograf yang digunakan dalam penelitian adalah data katalog rekaman akselerograf BMKG yang sudah diolah berdasarkan data katalog gempa bumi di wilayah Lombok dengan bantuan algoritma phyton yang terhubung oleh software SIGMA dan ShankeMap dengan memasukkan informasi seperti kejadian, lokasi, kedalaman, dan magnitudo tiap event. Hasil pembacaan amplitudo maksimum dari respon sensor akselerograf untuk mendapatkan nilai percepatan dari masing-masing komponen yaitu dua komponen horizontal dan satu komponen vertikal dalam satuan g ¿. Kemudia nilai percepatan tanah tersebut dikonversi satuannya menjadi gal ¿. Data rekaman akselerograf juga memberikan informasi mengenai letak sensor akselerograf.
30
Letak sensor akselerograf mewakili variasi keadaan lokal tanah suatu wilayah yang mempengaruhi nilai percepatan tanah yang dihasilkan. Sensor akselerograf di wilayah Lombok yang terpasang terdiri dari dua sensor akselerograf yaitu MASE dan TWSI dimana kedua sensor tersebut berada pada site class B, dengan tanah lokal yang didominasi oleh batuan dasar (bedrock) cenderung memiliki nilai percepatan tanah yang lebih kecil dibandingkan yang didominasi oleh sedimen atau tanah lunak (soil). Hal ini disebabkan karena gelombang seismik yang melewati lapisan bedrock akan terjadi perendaman amplitudo gelombang yang berakibat pada kecilnya nilai percepatan tanah yang dirasakan di permukaan. Sebaliknya pada lapisan soil
akan terjadi amplifikasi gelombang seismik
sehingga lebih rentan terhadap kerusaka akibat percepatan tanah yang tinggi di permukaan. Salah satu nilai percepatan tanah hasil rekaman akselerograf dapat dilihat pada tabel di bawah ini. Tabel 4.1 Nilai PGA hasil rekaman akselerograf Gempa bumi 5 Agustus 2018, jam 18:46:37 WIB, mag:7,0, Lat 8,35 LS, Long : 116,47 BT, Kedalaman : 15 Km, Lombok, Indonesia No
KODE
STASIUN
LINTANG
BUJUR
JARAK HIPOSENTER
1
TWSI TALIWANG
-8,738
116,8821
58,48
2
MASE STA MET BIL MATARAM
-8,766
116,279
48
3
SUBE STA MET SUMBAWA BESAR
-8,488
117,4135
DEBI
-8,738 -8,740
4 5
BALAI BESAR WIL III DENPASAR
DEMO STA MET NGURAH RAI DENPASAR
Z (gal)
N (gal)
E (gal)
INTENSITAS
SITE CLASS
11,1
14,4
18,89
III-IV
B
17,8762 43,4444 29,8665
IV
B
92,02
5,2783
IV
D
115,179
136,8
7,5891
NaN
16,7325
IV
D
115,18
137,1
6,0339
13,0301
20,334
IV
10,8682 13,5573
6
NEKI
STA KLIM NEGARA BALI
-8,340
114,62
192,69
2,0629
3,9876
2,7509
III
7
JAGI
JAJAG BANYUWANGI
-8,470
114,15
244,1
3,0831
1,6121
5,9408
III II-III
8
JBJN
D
-8,160
113,72
292,78
0,6566
1,6121
1,878
9
LABA STA MET LABUHAN BAJO
-8,480
119,88
362,04
0,1411
0,3518
0,3489
I-II
D
10
WGNI STA MET WANGAPU
-9,670
120,3
431,59
0,5576
1,0378
1,0388
II-III
E
11
BASO BARU SULAWESI
-4,530
119,62
538,14
0,2999
0,2744
0,2607
II
12
BJMI
-3,440
114,75
568,68
0,1656
0,4106
0,3254
II
13
MMRI MAUMERE
-8,640
122,24
621,99
0,148
3,9876
0,3381
I-II
14
UGM
WANAGAMA
-7,910
110,52
645,89
0,4381
NaN
1,0388
II-III
15
SOET
SOEI
-9,760
124,27
837,06
0,0715
NaN
13,5573
I-II
STA MET BANJARMASIN
D
4.1.3. Hasil penentuan persamaan empiris percepatan getaran tanah Dalam penelitian ini, diharapkan memperoleh suatu model rumus empiris percepatan tanah maksimum untuk wilayah Lombok. Untuk mendapatkan rumus empiris percepatan tanah, dalam penelitian ini menggunakan data observasi akselerograf BMKG di wilayah Lombok selama tahun 2018 beserta dengan parameter gempa yang diperoleh dari katalog gempa BMKG tahun 2018. Tercatat sebanyak 532 data rekaman akselerograf dari event gempa bumi di wilayah
31
Lombok dan sekitarnya. Melalui perhitungan analisis regresi non linier, maka diperoleh koefisien b 1 , b2, dan b 3 untuk rumus empiris percepatan getaran tanah di wilayah Lombok b 1=−2,525, b 2=0,712, dan b 3=0,395 Sehingga rumusan empiris model berdasarkan data penelitian tahun 2018 adalah sebagai berikut: log 10 ( Y )=−2,525+ 0,712∙ M −0,395 log 10 ( Rh +25)
(4.1)
Dimana Y adalah nilai percepatan tanah maksimum (gal), M L adalah magnitudo lokal, dan Rh adalah jarak hiposenter (km). Data survey dan pengelolah data terlampir (lampiran B) 4.2 Validasi hasil percepatan tanah model dengan percepatan tanah observasi Setelah memperoleh persamaan model yang baru, langkah selanjutnya adalah mencoba melakukan validasi dengan cara membandingkan hasil perhitungan dari persamaan yang baru dengan data observasi sensor accelerograf BMKG beserta persamaan awal dari Mc Guire (1997) sudah ada. Data yang dipergunakan untuk validasi adalah berjumlah 544 data observasi PGA pada kategori rock. Untuk memudahkan validasi, maka parameter PGA model yang baru dinormalisasikan pada magnitudo lokal ( M L ) 4,5 dan 5,1 dan site class rock. Hasil perbandingan nilai PGA model atenuasi dari data yang telah dinormalisasi tersebut selanjutnya akan dibandingkan berdasarkan fungsi jarak dari sumber ke stasiun dari 0-1000 km. 4.2.1. Hasil analisis model Fungsi atenuasi adalah hubungan antara intensitas getaran tanah dengan magnitudo dan jarak (parameter gempa bumi) terhadap suatu titik pengamatan. Model persamaan akhir ini secara fisis telah memperhitungkan saturasi dari amplitudo percepatan di daerah dekat sumber gempa. Dan model juga memiliki karakteristik simple point source
termasuk geometric spreading dalam jarak
32
menengah sampai jauh. Bentuk asumsi fungsi dasar dapat mewakili dua karakteristik tersebut adalah (Mc Guire, 1997). log ( Y )=b1 +b2 ∙ M −b3 log( Rh +25)
(4.2)
Pada hasil model ini menggunakan magnitudo lokal ( M L ) dan jarak hiposenter dari pusat gempa bumi. Dengan perhitungan menggunakan analisi regresi non linier, maka didapatkan koefisien dan konstanta sebagai berikut: Tabel 4.2 Koefisien hasil model
Parameter
Koefisien model
Nilai eror
Batas bawah
Batas atas
b1
-2.525
0.452
-3.413
-1.636
b2
0.712
0.034
0.645
0.779
b3
-0.395
0.226
-0.838
0.049
Bentuk persaman dari model setelah dimasukan semua koefisien yaitu: log ( Y )=−2,525+0,712∙ M −0,395 log(R h +25)
(4.3)
Model ini mempunyai nilai koefisien korelasi ( R2 ¿ yang rendah yaitu 0,45 dengan standar deviasi sebesar 0,544. persamaan model ini kemudian dibandingkan hasilnya untuk nilai magnitudo tertentu terhadap data observasi rekaman akselerograf.
Peak Ground Acceleration (gal)
00,010
00,001
00,000 00,001
Data Observasi Model
00,010 00,100 Jarak Hiposenter (km)
01,000
Gambar 4.2 Perbandingan nilai PGA hasil observasi dan model untuk magnitudo 4,6
33
Gambar 4.3 Perbandingan nilai PGA hasil observasi dan model untuk magnitudo 5,1
Perbandingan model hasil PGA dari persamaan baru hasil regresi pada beberapa tingkat magnitudo lokal ( M L ) yaitu 4,5 dan 5,1 Pada model tersebut terlihat pengaruh dari magnitudo terhadap nilai PGA. Semakin besar magnitudo gempa bumi yang terjadi, maka energi seismik yang menjalar dari sumber gempa bumi juga semakin besar dengan menghasilkan nilai percepatan tanah yang terekam oleh akselerograf. Semakin besar kecil perbedaan PGA terhadap perubahan magnitudo, maka pengaruh magnitudo terhadap perubahan nilai PGA seiring perubahan jarak hiposenter dari sumber gempa bumi akan semakin kecil. Garis model dapat mengikuti pola data observasi dari mulai jarak hiposenter 21,77 km hingga 100 km. Kemudian garis model pada jarak 100 km terus menurun signifikan hingga jarak 1000 km dan berada diantara data observasi. 4.2.2. Perbandingan nilai percepatan tanah model dengan metode percepatan tanah Mc Guire (1997) Validasi yang pertama dilakukan dengan membandingkan rumus model atenuasi yang dalam hal ini akan penulis sebut sebagai model dengan rekaman akselerograf. Hasil perbandingan ditampilkan dalam bentuk grafik. Berdasarkan Gambar 4.4 terlihat bahwa trend hasil perhitungan menggunakan rumus empiris
34
model memiliki pola yang paling mendekati dengan nilai PGA observasi. Nilai PGA dari persamaan empiris model sudah diduga sudah mendekati nilai PGA hasil rekaman alat akselerograf karena model dibentuk berdasarkan kondisi geologi yang sesuai dengan wilayah penelitian.
Percepatan tanah (gal)
10
1
0.1
0.01 1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
Jarak hiposenter (km) Rekaman akselerograf
Model
Gambar 4.4 Perbandingan Nilai PGA antara Observasi sensor, Model, dan rumus Mc Guire (1997)
4.3 Perbandingan antara Peta Kontur PGA Observasi, PGA Model Atenuasi, PGA Mc Guire (1997) Implementasi persamaan atenuasi model dalam bentuk peta digunakan untuk melihat besaran nilai percepatan tanah model atenuasi terhadap wilayah penelitian. Kegunaan lainnya untuk pemetaan suatu area yang terdampak dari kejadian gempa bumi dengan cepat dan membantu dalam memperkirakan nilai percepatan tanah di lokasi yang tidak terdapat stasiun akselerograf. Selain itu juga berguna untuk studi penilaian bahaya kegempaan. Pemetaan hasil model atenuasi akan dibandingkan dengan pemetaan data observasi dari rekaman akselerograf. Implementasi pemetaan dari nilai percepatan masing-masing stasiun menggunakan software Arcgis versi 10.4, data nilai percepatan tanah kemudiam diinterpolasi dengan teknik IDW, sehingga menghasilkan peta percepatan tanah
35
dari data asli dan data model yang di dapatkan. Skala nilai percepatan tanah pada pemetaan mengikuti rentang sebaran nilai percepatan tanah hasil model. Dalam percepatan tanah, apabila suatu gelombang melalui lapisan sedimen maka akan timbul suatu resonansi. Oleh karena itu, bangunan-bangunan yang berada diatasnya akan menerima getaran-getaran tersebut, dimana arahnya dapat diuraikan menjadi dua komponen yaitu komponen horizontal dan komponen vertikal. Untuk getaran vertikal, pada umumnya kurang membahayakan sebab searah
dengan
gaya
gravitasi.
Sedangkan
untuk
komponen
horizontal
menyebabkan keadaan bangunan seperti diayun. Bila bangunan itu tinggi maka dapat diumpamakan seperti bandul yang mengalami getaran paksaan (force vibration), dan ini sangat membahayakan. Perbandingan sebaran nilai percepatan tanah hasil rekaman akselerograf dan model pada Gambar 4.5 dan Gambar 4.6 menunjukkan nilai percepatan tanah terbesar tersebar di sekitar wilayah Lombok Utara dan Lombok Timur, dimana pada wilayah tersebut struktur geologinya tersusun oleh material tanah lunak berupa material pasir dan aluvium. Karakteristik tanah lunak semacam ini dapat menimbulkan resonansi gelombang gempa yang menyebabkan terjadinya amplifikasi guncangan gempa. Nilai percepatan hasil rekaman akselerograf dan nilai percepatan dari model masing–masing sebesar 53,44 gal dan 65,22 gal. Nilai percepatan tanah model sudah mendekati nilai percepatan tanah hasil rekaman akselerograf walaupun ada sedikit perbedaan nilai percepatan tersebut bisa diakibatkan karena nilai percepatan hasil rekaman akselerograf memiliki jarak dari sumber gempa buminya sebesar 48 km, sedangkan nilai percepatan tanah hasil model diasumsikan jarak sumber gempa dengan titik penelitian sebesar 17,83 km. Hal tersebut yang membuat nilai percepatan tanah hasil rekaman akselerograf berbeda dengan model, karena besarnya nilai percepatan tanah dipengaruh oleh besarnya magnitudo dan jarak pusat gempa dengan titik penelitian. Semakin dekat suatu wilayah dengan dari pusat gempa bumi maka nilai percepatan tanah maksimum akan semakin besar dan sebaliknya.
36
Gambar 4.5 Peta pecepatan tanah dari data rekaman akselerograf (BMKG)
37
Gambar 4.6 Peta pecepatan tanah dari hasil model
BAB V KESIMPULAN DAN SARAN 5.1 Kesimpulan Berdasarkan pembahasan yang telah dipaparkan, maka kesimpulan yang dapat diberikan yaitu: 1. Persamaan atenuasi yang sesuai dengan data rekaman akselerograf untuk gempa bumi di wilayah Lombok pada jarak menengah hingga jauh, yaitu:
log ( A ) =−2,525+0,712 ∙ M −0,395 log ( R+25) 2. Nilai percepatan tanah model atenuasi berdasarkan persamaan yang diperoleh sesuai dengan percepatan tanah observasi yang memiliki nilai korelasi sebesar 0,45 berarti berkorelasi rendah. 3. Penerapan (implementasi) persamaan atenuasi hasil model dalam penelitian daerah terdampak akibat gempa bumi di wilayah Lombok, menunjukan hasil yang lebih mendekati data observasi.
5.2 Saran Saran pada penelitian selanjutnya, yaitu 1. Menggunakan model persamaan atenuasi yang melibatkan parameter seperti kondisi geologis setempat dengan masukan variabel bebas (independent) dalam perhitungan atenuasi. 2. Menggunakan data parameter gempa bumi yang sudah direlokasi dalam penelitian 3. Untuk meningkatkan akurasi model rumus empiris percepatan tanah ini diperlukan data observasi yang banyak mengingat distribusi sensor akselerograf yang masih belum merata di wilayah Lombok. Oleh karena itu perlu adanya peningkatan jaringan akselerograf, dengan dukungan akurasi alat yang baik.
37
DAFTAR PUSTAKA Abrahamson, N. a. (1997). Empirical Response Spectral Attenuation Relations For Shallow Crustal Earthquake. Seismological Research Letters, 68(1), pp.94-127. Agustawijaya D, S. S. (2006). Pengkajian Bahaya Longsor Tipe Jatuhan Batuan (Rockfall). Laporan Tahap II Penelitian Hibah Bersaing, Direktorat Jenderal Pendidikan Tinggi, Departemen Pendidikan Nasional. Agustawijaya, S Sulistyowati, T Suroso, dan Hadi. (2006). Laporan Tahap II Penelitian Hibah Bersaing, Direktorat Jenderal Pendidikan Tinggi, Departemen Pendidikan Nasional. Agustawijaya, S. S. (2005). Laporan Tahap I Penelitian Hibah Bersaing, Direktorat Jenderal Pendidikan Tinggi, Departemen Pendidikan Nasional. Andi, M. S. (1994). Geological Map Sheet Lombok, NTB. Geological Research and Development Center. Asrurifak, M. (2010). Peta Respon Spektra Indonesia Untuk Perencanaan Struktur Bangunan Tahan Gempa Berdasarkan Model Sumber Gempa Tiga Dimensi Dalam Analisis Probabilitas, Disertasi. Bandung: Institut Teknologi Bandung. BMKG. (2017). Indeks Seismisitas. Jakarta Pusat: BMKG. Boore, D. B. (2004). Processing of strong-motion accelerograms: needs,options and consequences. Soil Dynamics and Earthquake Engineering, 25(2), 93115. Bormann, P. (2002). IASPEI New Manual of Seismological Observatory Practice (NMSOP). Vol. 1 dan 2. Dadafshar, M. (2014). Accelerometer and Gyroscopes Sensors: Operation, Sensing, and Applications. Aplication Note 5830 . Douglass. (2011). Ground Motion Prediction Equations. Pasific Earthquake Engineering Research Center College, Engineering University of California, Barkeley. doi:2011/102
38
Fukushima, Y. d. (1990). A New Attenuation Relation For Peak Horizontal Acceleration Of Strong Earthquake Ground Motion In Japan. Bulletin of the Seismological Society of America, 80(4), 757-783. Gadallah, M. F. (2009). Exploration Geophysics. Houston: Springer-Verlag Berlin Heidelberg. Grandis, D. H. (2009). Pengantar Pemodelan Inversi Geofisika. Bandung: Himpunan Ahli Geofisika Indonesia (HAGI). Idris, L. (1990). Response of Soft Soil Sites during Earthquakes. Proc.,H. Boltom Seed Memorial Symp, 2, 273-290. Irsyam M., H. D. (2007, Agustus 21-22). Usulan Revisi Peta Hazard Kegempaan Wilayah Indonesia . Paper Seminar HAKI, Konstruksi Tahan Gempa Di Indonesia. Irsyam, M. S. (2010). Ringkasan Hasil Studi Tim Revisi Peta Gempa Bumi Indonesia 2010. Kementrian Pekerja Umum. Kramer, S. L. (1996). Geotechnical Earthquake Engineering. United States of America: Prentice Hall, Inc. Mc Guire, R. K. (1977). Seismic Design Spectra And Mapping Procedures Using Hazard Analysis Based Directly On Oscillator Response. Earthquake Engineering and Structural Dynamics, 5, 211-234. Octantyo, A. Y. (2016). Model Atenuasi Percepatan Tanah Wilayah Bali dan Nusa Tenggara Barat Untuk Gempabumi Interface dan Intralsab,Skripsi. Tanggerang Selatan: Sekolah Tinggi Meteorologi Klimatologi dan Geofisika. Pawirodikromo,
W.
(2012).
Seismologi
Teknik
Rekayasa
Kegempaan.
Yogyakarta: Pusaka Pelajar. Petrus, D. S. (2013). Penentuan seismisitas dan tingkat risiko gempa bumi. Malang: UB Press. Purwana, I. (2017). Materi Kuliah Seismologi Lanjut Seri 1. Tangerang Selatan: Sekolah Tinggi Meteorlogi Klimatologi dan Geofisika.
39
Rudyanto, A. (2013). Development of Strong Motion Database For The SumatraJava Region, Tesis, Program Master Filosofi. Australia: The Australian National Universirty. Sarwono. (2006). Metode Penelitian Kuantitatif dan Kualitatif. Yogyakarta: Graha Ilmu. Setiawan, Y. (2012). Kajian Rumus Empiris Percepatan Tanah di Daerah Bali dan Sekitarnya, Skripsi. Bandung: Program Studi Meteorlogi Institut Teknologi Bandung. Srinivasan, C. W. (2012). , Strong Ground Motion Prediction Equation for Low Magnitude and Near-field Earthquake Data for Shield Region In India. Bengaluru: National Institute of Rock Mechanics. Sunarjo, G. M. (2012). Gempa bumi Edisi Populer. Sungkowo, A. (2018). Perhitungan Nilai Percepatan Tanah Maksimum Berdasar Rekaman Sinyal Accelerograph di Stasiun Pengukuran UNSO Surakarta. Indonesian Journal of Applied Physics, 8 No 1, 43. Sutopo. (2007). Studi Faktor Kualitas (Q-Faktor) Gelombang Kompresi (P-wave) Pada Batuan Sedimen Dangkal. 10(1).
40
LAMPIRAN
41
Lampiran A Peta Geologi Lombok
42
Lampiran B Data parameter gempa rekaman akselerograf Origin Time 30/07/201 8 29/07/201 8 29/07/201 8 20/08/201 8 20/08/201 8 19/08/201 8 19/08/201 8 19/08/201 8 19/08/201 8 18/08/201 8 18/08/201 8 15/08/201 8 13/08/201 8 12/08/201 8 10/08/201 8 11/08/201 8 10/08/201 8 07/08/201 8 06/08/201 8 06/08/201 8 06/08/201 8
Lintang
Bujur
Kedalama n
Magnitudo
Resultan PGA
Rhypo
18:07:13
-8.37
116.43
10
3.5
0.1071
37.267
11:31:55
-8.85
116.27
10
3.5
0.0642
35.344
1:04:51
-8.63
116.34
12
3.5
0.1083
23.754
19:28:41
-8.36
116.33
17
3.5
0.0219
33.357
17:40:01
-8.61
116.58
27
3.5
0.0180
53.080
16:54:07
-8.28
116.6
10
3.5
0.0376
57.871
7:26:37
-8.32
116.61
10
3.5
0.0458
56.492
6:26:09
-8.43
116.56
10
3.5
0.0629
46.684
4:48:46
-8.35
116.62
10
3.5
0.0347
56.000
23:47:36
-8.29
116.22
15
3.5
0.3760
34.258
16:06:45
-8.25
116.23
14
3.5
0.2333
38.068
3:10:57
-8.25
116.23
10
3.5
0.0890
36.785
11:06:53
-8.71
116.29
10
3.5
0.0908
23.531
7:15:32
-8.27
116.16
10
3.5
0.2025
34.004
23:14:03
-8.26
116.08
11
3.5
0.1894
36.592
19:16:38
-8.25
116.26
10
3.5
0.0368
37.459
10:33:57
-8.41
116.58
10
3.5
0.0587
49.422
8:40:01
-8.39
116.14
25
3.5
0.0800
31.583
17:51:51
-8.23
116.36
27
3.5
0.0424
50.436
17:37:32
-8.3
116.24
17
3.5
0.0665
34.656
14:24:29
-8.37
116.12
11
3.5
0.1318
24.597
43
06/08/201 8 06/08/201 8 06/08/201 8 05/08/201 8 01/08/201 8 01/08/201 8 03/09/201 8 26/11/201 8 29/07/201 8 30/07/18 29/07/201 8 29/07/18 29/07/18 20/08/18 20/08/18 19/08/18 19/08/18 19/08/18 19/08/18 18/08/18 15/08/201 8 13/08/201 8 12/08/18 12/08/201 8 10/08/18 11/08/18 10/08/18 10/08/201 8 10/08/18 07/08/201 8 06/08/201 8
9:22:07
-8.27
116.18
15
3.5
0.0644
35.795
4:43:31
-8.37
116.19
10
3.5
0.8671
23.712
4:03:05
-8.33
116.01
14
3.5
0.1554
34.015
20:32:56
-8.28
116.19
15
3.5
0.0815
34.799
10:06:34
-8.36
116.63
11
3.5
0.0416
56.678
0:03:15
-8.43
116.51
13
3.5
0.1510
42.535
21:13:08
-8.24
116.12
11
3.5
0.0210
37.834
10:46:14
-8.22
116.5
10
3.5
0.0234
53.739
21:10:17
-8.81
116.29
10
3.5
0.0887
66.259
18:07:13
-8.37
116.43
10
3.5
0.3152
65.073
11:31:55
-8.85
116.27
10
3.5
0.1543
69.127
7:24:47 1:04:51 19:28:41 17:40:01 16:54:07 7:26:37 6:26:09 4:48:46 23:47:36
-8.25 -8.63 -8.36 -8.61 -8.28 -8.32 -8.43 -8.35 -8.29
116.49 116.34 116.33 116.58 116.6 116.61 116.56 116.62 116.22
10 12 17 27 10 10 10 10 15
3.5 3.5 3.5 3.5 3.5 3.5 3.5 3.5 3.5
0.1855 0.2952 0.8976 1.6882 0.9823 5.6562 8.4972 1.5200 0.1437
70.018 61.875 75.733 45.080 60.433 56.197 50.206 52.855 89.466
3:10:57
-8.25
116.23
10
3.5
0.0519
90.454
11:06:53
-8.71
116.29
10
3.5
0.2940
65.864
16:38:35
-8.37
116.35
28
3.5
0.2060
76.694
7:15:32
-8.27
116.16
10
3.5
0.2064
95.425
23:14:03 19:16:38 16:08:52
-8.26 -8.25 -8.58
116.08 116.26 116.14
11 10 10
3.5 3.5 3.5
0.3348 0.0346 0.0953
103.486 87.871 83.997
10:33:57
-8.41
116.58
10
3.5
2.4115
50.304
21:11:54
-8.28
116.03
11
3.5
0.1898
107.166
8:40:01
-8.39
116.14
25
3.5
0.0352
93.697
17:51:51
-8.23
116.36
27
3.5
0.1424
84.907
44
06/08/201 8 06/08/201 8 06/08/201 8 06/08/201 8 06/08/201 8 01/08/201 8 01/08/201 8 27/12/18 29/07/201 8 29/07/201 8 29/07/201 8 29/07/201 8 19/08/201 8 18/08/201 8 09/08/201 8 09/08/201 8 07/08/201 8 07/08/201 8 06/08/201 8 06/08/201 8 05/08/201 8 05/08/201 8 05/08/201 8 05/08/201 8 05/08/201 8
17:37:32
-8.3
116.24
17
3.5
0.1834
87.370
14:24:29
-8.37
116.12
11
3.5
0.1709
93.847
9:22:07
-8.27
116.18
15
3.5
0.1281
94.301
4:43:31
-8.37
116.19
10
3.5
1.4490
86.977
4:03:05
-8.33
116.01
14
3.5
0.1968
106.921
10:06:34
-8.36
116.63
11
3.5
0.3417
51.489
0:03:15
-8.43
116.51
13
3.5
0.9113
54.862
11:48:54
-8.36
116.19
10
3.5
0.1317
87.473
22:34:07
-8.44
116.45
12
3.6
0.1123
35.961
14:10:46
-8.88
116.27
10
3.6
0.0583
38.424
6:52:09
-8.42
116.44
12
3.6
0.0696
35.961
6:25:18
-8.86
116.28
32
3.6
0.1026
47.720
20:11:59
-8.45
116.68
23
3.6
0.0394
62.115
12:39:07
-8.36
116.19
10
3.6
0.6883
24.714
9:06:55
-8.27
116.5
10
3.6
0.0714
49.912
3:27:13
-8.53
116.14
20
3.6
0.1534
20.522
12:02:45
-8.5
116.13
11
3.6
0.1333
13.601
5:47:15
-8.25
116.18
10
3.6
0.1741
36.112
20:13:01
-8.24
116.3
10
3.6
0.1706
39.971
15:57:11
-8.37
116.44
10
3.6
0.1097
38.134
15:49:04
-8.37
116.4
10
3.6
0.0748
34.769
15:18:31
-8.31
116.45
10
3.6
0.1033
43.077
13:32:07
-8.33
116
10
3.6
0.2757
33.240
13:30:21
-8.22
116.27
10
3.6
0.3943
40.940
12:48:07
-8.21
116.16
10
3.6
0.0883
40.359
45
05/08/201 8 29/08/201 8 26/08/201 8 03/09/201 8 02/10/201 8 18/12/201 8 29/07/201 8 29/07/201 8 29/07/201 8 29/07/201 8 29/07/201 8 29/07/18 19/08/201 8 19/08/18 18/08/201 8 14/08/18 09/08/201 8 07/08/201 8 07/08/201 8 06/08/201 8 06/08/201 8 06/08/18 29/08/201 8 26/08/201 8 23/08/18 18/12/201 8 30/07/201 8
12:40:09
-8.35
116.2
10
3.6
0.2061
25.815
16:46:15
-8.2
116.55
15
3.6
0.0258
60.102
19:52:38
-8.28
116.02
10
3.6
0.0458
36.689
3:35:45
-8.29
116.2
10
3.6
0.0962
32.002
2:07:55
-8.74
116.47
13
3.6
0.0884
40.825
7:41:42
-8.33
116.55
10
3.6
0.3168
50.402
22:34:07
-8.44
116.45
12
3.6
0.2494
59.130
14:10:46
-8.88
116.27
10
3.6
0.3630
66.295
6:52:09
-8.42
116.44
12
3.6
0.1923
69.721
6:25:18
-8.86
116.28
32
3.6
0.3836
61.286
6:25:18
-8.86
116.28
32
3.6
0.3836
74.701
2:33:45
-8.51
116.36
18
3.6
0.1725
65.233
20:11:59
-8.45
116.68
23
3.6
2.5554
45.277
4:27:49
-8.36
116.64
10
3.6
0.4027
50.696
12:39:07
-8.36
116.19
10
3.6
3.1205
87.473
15:41:25
-8.35
116.69
10
3.6
0.1381
49.031
3:27:13
-8.53
116.14
20
3.6
0.1830
67.643
12:02:45
-8.5
116.13
11
3.6
0.1745
87.029
5:47:15
-8.25
116.18
10
3.6
0.1267
87.494
20:13:01
-8.24
116.3
10
3.6
0.2017
94.829
15:57:11
-8.37
116.44
10
3.6
0.4801
85.189
15:57:11
-8.37
116.44
10
3.6
0.4801
64.283
16:46:15
-8.2
116.55
15
3.6
0.6918
71.687
19:52:38
-8.28
116.02
10
3.6
0.0722
108.064
20:06:59
-8.41
115.93
15
3.6
0.7716
111.811
7:41:42
-8.33
116.55
10
3.6
0.6427
59.053
16:59:26
-8.33
116.48
10
3.7
0.0795
44.147
46
30/07/201 8 30/07/201 8 30/07/201 8 29/07/201 8 19/08/201 8 11/08/201 8 07/08/201 8 07/08/201 8 06/08/201 8 06/08/201 8 05/08/201 8 05/08/201 8 05/08/201 8 05/08/201 8 05/08/201 8 05/08/201 8 04/08/201 8 25/08/201 8 03/09/201 8 21/09/201 8 11/11/201 8 29/07/201 8 5/8/18 25/08/18 21/09/18 21/09/18 11/12/201
1:57:44
-8.43
116.55
10
3.7
0.0546
45.707
0:42:47
-8.34
116.47
10
3.7
0.0450
42.591
23:35:03
-8.4
116.56
10
3.7
0.1979
47.856
17:55:43
-8.27
116.61
10
3.7
0.0929
59.644
23:38:36
-8.26
116.59
10
3.7
0.0790
58.166
1:50:59
-8.2
116.21
10
3.7
0.2367
41.716
16:26:56
-8.22
116.23
12
3.7
0.0888
40.422
3:18:07
-8.25
116.25
10
3.7
0.1544
37.267
20:04:38
-8.32
116.1
11
3.7
0.2313
29.990
4:21:08
-8.21
116.5
10
3.7
0.0647
54.525
19:20:04
-8.49
116.17
17
3.7
1.9068
18.788
17:59:42
-8.23
116.61
10
3.7
0.1678
61.913
16:40:14
-8.24
116.2
10
3.7
0.2901
37.363
15:46:56
-8.24
116.03
10
3.7
0.2893
40.068
13:24:36
-8.34
116.3
10
3.7
0.1576
30.392
13:02:23
-8.28
116.02
10
3.7
0.5165
36.689
7:00:10
-8.28
116.64
10
3.7
0.0710
61.518
4:15:07
-8.28
116.16
10
3.7
0.6230
32.859
23:18:37
-8.25
116.11
11
3.7
0.0776
36.879
7:54:35
-8.32
116.14
13
3.7
0.4189
30.057
17:30:12
-8.3
116.06
11
3.7
0.2698
33.271
18:32:39
-8.2
116.5
10
3.7
0.2386
36.522
23:40:13 11:15:07 14:54:35 21:43:28 00:30:11
-8.21 -8.28 -8.32 -8.06 -8.27
116.21 116.16 116.12 116.24 116.06
19 10 16 10 10
3.7 3.7 3.7 3.7 3.7
0.2976 0.6307 0.4435 0.0248 0.2210
36.222 56.492 54.980 79.234 61.025
47
8 30/07/201 8 30/07/201 8 30/07/201 8 30/07/201 8 29/07/201 8 29/07/18 19/08/201 8 13/08/18 11/08/201 8 07/08/201 8 07/08/201 8 06/08/201 8 06/08/201 8 04/08/201 8 02/08/18 25/08/201 8 18/12/201 8 29/07/201 8 29/07/201 8 30/07/201 8 25/08/201 8 29/07/201 8 29/07/201 8 18/08/201 8 12/08/201 8 08/08/201
16:59:26
-8.33
116.48
10
3.7
0.2245
64.085
1:57:44
-8.43
116.55
10
3.7
1.4355
50.990
0:42:47
-8.34
116.47
10
3.7
0.2154
64.085
23:35:03
-8.4
116.56
10
3.7
1.1671
52.560
17:55:43
-8.27
116.61
10
3.7
2.6853
60.828
1:55:31
-8.86
116.23
15
3.7
0.5150
74.427
23:38:36
-8.26
116.59
10
3.7
2.0053
62.900
14:33:22
-8.37
116.1
32
3.7
0.0155
100.434
1:50:59
-8.2
116.21
10
3.7
0.1451
95.525
16:26:56
-8.22
116.23
12
3.7
0.0754
92.680
3:18:07
-8.25
116.25
10
3.7
0.1216
88.666
20:04:38
-8.32
116.1
11
3.7
0.3983
98.317
4:21:08
-8.21
116.5
10
3.7
0.2625
72.889
7:00:10
-8.28
116.64
10
3.7
0.8403
58.265
11:17:00
-8.26
116.69
10
3.7
0.2566
58.068
4:15:07
-8.28
116.16
10
3.7
0.4100
94.829
15:53:32
-8.26
116.18
10
3.7
0.1198
94.232
8:54:29
-8.86
116.23
15
3.7
0.7148
62.623
18:32:39
-8.2
116.5
10
3.7
1.5501
73.781
8:57:44
-8.43
116.55
10
3.7
1.4160
51.088
11:15:07
-8.28
116.16
10
3.7
0.4171
94.829
18:07:57
-8.42
116.47
10
3.8
0.1420
38.231
0:26:10
-8.28
116.61
10
3.8
0.0761
58.757
15:12:37
-8.27
116.03
10
3.8
0.2992
37.170
0:59:16
-8.32
116.65
15
3.8
0.0332
61.362
11:41:32
-8.29
116.14
10
3.8
0.4445
32.002
48
8 06/08/201 8 05/08/201 8 05/08/201 8 02/08/201 8 01/08/201 8 23/08/201 8 21/08/201 8 12/09/201 8 11/09/201 8 29/07/201 8 29/07/201 8 13/09/201 8 29/07/201 8 29/07/18 29/07/18 29/07/201 8 18/08/201 8 12/08/201 8 11/08/18 08/08/201 8 06/08/201 8 02/08/18 01/08/18 23/08/18 21/08/18 13/06/201 8 29/07/201 8
3:15:23
-8.22
116.01
18
3.8
0.2677
45.419
15:45:05
-8.23
116.5
10
3.8
0.1002
52.953
13:56:25
-8.22
116.48
10
3.8
0.1448
52.266
15:49:25
-8.26
116.48
10
3.8
0.0953
49.129
10:12:21
-8.26
116.67
10
3.8
0.0830
65.468
18:53:50
-8.22
116.14
10
3.8
0.2155
39.390
18:51:42
-8.28
116.64
10
3.8
0.0408
61.518
21:49:36
-8.61
116.07
11
3.8
0.1214
16.205
17:59:34
-8.24
116.65
10
3.8
0.0859
64.777
7:26:10
-8.28
116.61
10
3.8
0.0821
61.806
9:21:50
-8.4
116.49
10
3.8
0.9305
68.138
04:49:35
-8.31
116.18
10
3.8
0.1201
52.756
18:07:57
-8.42
116.47
10
3.8
0.4218
59.328
1:54:29 1:43:15
-8.63 -8.75
116.34 116.32
17 18
3.8 3.8
0.5151 1.2025
59.865 64.272
0:26:10
-8.28
116.61
10
3.8
1.0725
59.841
15:12:37
-8.27
116.03
10
3.8
0.3207
107.566
0:59:16
-8.32
116.65
15
3.8
0.2660
55.082
4:59:59
-8.29
116.69
10
3.8
0.3077
55.016
11:41:32
-8.29
116.14
10
3.8
0.2920
96.122
3:15:23
-8.22
116.01
18
3.8
0.2566
113.240
15:49:25 10:12:21 18:53:50 18:51:42
-8.26 -8.26 -8.22 -8.28
116.48 116.67 116.14 116.64
10 10 10 10
3.8 3.8 3.8 3.8
0.4700 0.4641 0.2219 0.3363
69.820 58.856 100.300 58.265
18:23:41
-8.49
116.07
10
3.8
0.5109
93.934
7:26:10
-8.28
116.61
10
3.8
1.0649
59.940
49
29/07/201 8 29/07/201 8 17/11/201 8 15/11/201 8 29/07/201 8 29/07/201 8 29/07/201 8 29/07/201 8 29/07/201 8 19/08/201 8 07/08/201 8 06/08/201 8 06/08/201 8 06/08/201 8 06/08/201 8 05/08/201 8 05/08/201 8 05/08/201 8 29/07/201 8 5/8/18 8/8/18 10/05/201 8 15/11/201 8 17/11/201 8 17/11/201 8 29/07/18
9:21:50
-8.4
116.49
10
3.8
3.8743
58.068
9:29:54
-8.1
116.64
10
3.8
0.5780
76.457
5:41:17
-8.35
116.09
11
3.9
1.0825
27.313
12:46:31
-8.21
116.63
5
3.9
0.2303
64.974
19:48:27
-8.39
116.57
10
3.9
0.2450
49.226
7:54:11
-8.37
116.55
10
3.9
0.4867
48.248
3:27:35
-8.2
116.45
10
3.9
0.0891
51.873
2:56:08
-8.51
116.45
13
3.9
0.1340
33.802
0:49:12
-8.27
116.5
14
3.9
0.1074
50.865
16:59:31
-8.34
116.42
21
3.9
0.1071
42.718
0:13:27
-8.29
116.18
10
3.9
0.1721
31.908
21:24:46
-8.29
116.64
10
3.9
0.0387
60.926
14:52:07
-8.24
116.2
10
3.9
0.1959
37.363
3:48:46
-8.25
116.7
10
3.9
0.0560
68.830
2:34:32
-8.23
116.3
15
3.9
0.2353
42.439
15:34:03
-8.21
116.23
10
3.9
0.4461
40.940
15:20:43
-8.22
116.17
10
3.9
0.3385
39.294
12:33:31
-8.3
116.04
29
3.9
0.4523
43.408
7:49:12
-8.27
116.5
14
3.9
0.0969
47.856
22:20:42 18:30:03
-8.21 -7.83
116.17 116.79
10 10
3.9 3.9
0.3643 0.0692
49.070 25.854
9:52:11
-8.42
116.87
11
3.9
0.3780
76.297
19:46:30
-8.21
116.63
5
3.9
0.2171
73.071
2:41:17
-8.35
116.09
11
3.9
1.0739
51.880
14:16:48
-8.01
116.49
10
3.9
0.0410
87.772
19:48:27
-8.39
116.57
10
3.9
1.6769
52.658
50
29/07/18 29/07/18 29/07/201 8 29/07/201 8 19/08/201 8 15/08/18 10/08/18 09/08/18 07/08/201 8 06/08/201 8 06/08/201 8 06/08/201 8 06/08/201 8 29/07/201 8 29/07/201 8 08/08/18 05/10/18 28/07/201 8 19/08/201 8 18/08/201 8 15/08/201 8 15/08/201 8 10/08/201 8 08/08/201 8 06/08/201 8 05/08/201 8 05/08/201 8 23/08/201
7:54:11 3:27:35
-8.37 -8.2
116.55 116.45
10 10
3.9 3.9
3.8493 0.4350
55.705 77.052
2:56:08
-8.51
116.45
13
3.9
0.5038
55.348
0:49:12
-8.27
116.5
14
3.9
0.5195
68.349
16:59:31
-8.34
116.42
21
3.9
1.0002
70.500
9:23:35 20:59:40 8:29:20
-8.32 -8.41 -8.28
116.6 116.1 116.07
18 16 10
3.9 3.9 3.9
1.2906 0.6989 0.2353
58.727 94.761 103.186
0:13:27
-8.29
116.18
10
3.9
0.0868
92.343
21:24:46
-8.29
116.64
10
3.9
0.1852
57.378
14:52:07
-8.24
116.2
10
3.9
0.3468
93.735
3:48:46
-8.25
116.7
10
3.9
0.4240
58.659
2:34:32
-8.23
116.3
15
3.9
0.4025
86.609
7:49:12
-8.27
116.5
14
3.9
0.5368
68.349
9:56:08
-8.56
116.43
12
3.9
0.5236
54.829
18:30:03 9:52:11
-7.83 -8.42
116.79 116.87
10 11
3.9 3.9
1.7335 30.9720
101.991 37.070
23:03:41
-8.34
116.55
12
4
0.2164
50.254
22:28:19
-8.2
116.62
10
4
0.0940
64.777
20:12:57
-8.3
116.06
11
4
0.5610
33.271
9:17:17
-8.26
116.59
10
4
0.0338
58.166
9:08:41
-8.23
116.58
10
4
0.0540
59.349
10:40:01
-8.42
116.57
10
4
0.0814
48.052
0:50:05
-8.21
116.09
12
4
0.8366
41.761
14:18:14
-8.3
116.15
11
4
0.6234
31.570
22:44:03
-8.29
116.02
19
4
0.9126
39.211
15:25:57
-8.45
116.08
10
4
0.1588
18.614
21:26:30
-8.31
116.68
15
4
0.0786
64.761
51
8 08/12/201 8 7/8/18 11/8/18 25/08/201 8 09/04/201 8 09/11/201 8 11/08/201 8 28/07/201 8 15/08/18 15/08/201 8 15/08/201 8 10/08/201 8 08/08/201 8 06/08/201 8 23/08/201 8 29/07/201 8 07/08/18 11/08/18 25/08/201 8 29/07/201 8 29/07/201 8 28/07/201 8 19/08/201 8 11/08/201 8 06/08/201 8 06/08/201 8
14:00:41
-8.21
116.61
10
4
0.1003
63.295
21:04:38 22:42:08 12:56:4 6
-8.23 -8.19
116.23 116.32
14 10
4.0 4.0
0.3069 0.8842
106.869 45.805
-8.73
116.17
14
4
0.7928
18.835
2:41:41
-8.28
116.12
12
4
0.3060
58.054
09:40:05
-8.32
116.19
10
4
0.1180
51.481
21:58:10
-8.26
116.79
20
4
0.0640
81.977
23:03:41
-8.34
116.55
12
4
1.5769
58.543
19:48:23
-8.25
116.63
10
4
4.7980
61.716
9:17:17
-8.26
116.59
10
4
1.5540
62.900
9:08:41
-8.23
116.58
10
4
0.5730
66.259
10:40:01
-8.42
116.57
10
4
4.9776
50.206
0:50:05
-8.21
116.09
12
4
0.6176
105.683
14:18:14
-8.3
116.15
11
4
0.4017
94.763
21:26:30
-8.31
116.68
15
4
4.8511
54.601
6:03:41
-8.34
116.55
12
4
3.4855
58.641
21:04:38 22:42:08
-8.23 -8.19
116.23 116.32
14 10
4 4
0.1825 1.0369
92.367 87.374
12:56:46
-8.73
116.17
14
4
0.9156
79.542
1:14:21
-8.65
116.37
13
4.1
0.1637
27.647
0:32:02
-8.2
116.6
10
4.1
0.1601
63.097
23:38:04
-8.3
116.57
10
4.1
0.2427
53.935
4:45:35
-8.38
116.58
10
4.1
0.2629
50.696
22:21:05
-8.36
116.59
10
4.1
0.1360
52.560
23:16:03
-8.21
116.38
10
4.1
0.3489
46.684
19:27:31
-8.2
116.33
10
4.1
0.7182
45.122
52
06/08/201 8 06/08/201 8 05/08/201 8 05/08/201 8 05/08/201 8 05/08/201 8 05/08/201 8 27/08/201 8 11/09/201 8 03/09/201 8 02/09/201 8 29/07/201 8 7/8/18 27/08/201 8 28/08/201 8 29/07/201 8 29/07/201 8 28/07/201 8 19/08/201 8 11/08/201 8 06/08/201 8 06/08/201 8 06/08/201 8 27/08/201 8 02/09/201 8
15:22:11
-8.23
116.43
12
4.1
0.1187
48.411
0:34:59
-8.29
116.25
12
4.1
0.4869
33.802
22:28:10
-8.24
116.06
11
4.1
0.1360
39.272
19:44:56
-8.23
116.25
10
4.1
0.9115
39.294
17:11:09
-8.3
116.41
10
4.1
0.1524
40.746
17:09:57
-8.36
116.57
10
4.1
0.6358
50.696
13:45:21
-8.36
116.13
11
4.1
0.8547
25.315
14:12:28
-8.29
116.59
10
4.1
0.1082
56.394
2:40:05
-8.28
116.19
10
4.1
0.0577
32.954
19:41:42
-8.28
116.12
12
4.1
0.4321
35.440
21:48:02
-8.25
116.65
10
4.1
0.1044
64.184
6:38:04
-8.3
116.57
10
4.1
0.2448
55.610
7:13:25
-8.17
116.23
10
4.1
0.1798
40.122
21:12:29
-8.29
116.59
10
4.1
0.1139
63.789
10:13:29
-8.42
116.64
10
4.1
0.2544
56.197
1:14:21
-8.65
116.37
13
4.1
0.6351
58.561
0:32:02
-8.2
116.6
10
4.1
1.2683
68.138
23:38:04
-8.3
116.57
10
4.1
0.6949
60.433
4:45:35
-8.38
116.58
10
4.1
8.0334
52.756
22:21:05
-8.36
116.59
10
4.1
3.5572
53.837
23:16:03
-8.21
116.38
10
4.1
0.6683
81.218
19:27:31
-8.2
116.33
10
4.1
0.8852
85.785
15:22:11
-8.23
116.43
12
4.1
0.5876
76.151
14:12:28
-8.29
116.59
10
4.1
1.8319
60.039
21:48:02
-8.25
116.65
10
4.1
0.6871
60.729
53
22/12/18 29/07/201 8 29/07/201 8 07/08/18 27/08/201 8 28/08/201 8 29/07/201 8 29/07/201 8 07/08/201 8 06/08/201 8 28/08/201 8 29/07/201 8 29/07/201 8 6/8/18 19/08/201 8 20/08/201 8 09/03/201 8 29/07/201 8 29/07/201 8 07/08/201 8 06/08/201 8 28/08/201 8 29/07/201 8 29/07/201 8 19/08/201 8 20/08/201 8
23:59:50
-8.31
116.08
10
4.1
0.7470
100.698
6:38:04
-8.3
116.57
10
4.1
0.9169
60.433
7:32:02
-8.2
116.6
10
4.1
1.2910
68.138
7:13:25
-8.17
116.23
10
4.1
0.1384
96.122
21:12:29
-8.29
116.59
10
4.1
1.8470
60.137
10:13:29
-8.42
116.64
10
4.1
10.8954
45.415
4:22:51
-8.31
116.49
13
4.2
0.1597
47.032
1:53:17
-8.61
116.36
15
4.2
0.2892
26.544
14:04:40
-8.27
116.21
15
4.2
0.3416
36.067
17:49:28
-8.28
116.25
12
4.2
0.5255
34.738
3:13:29
-8.43
116.6
10
4.2
0.2335
55.213
10:05:01
-8.33
116.44
10
4.2
0.1968
42.591
0:18:34
-8.19
116.58
10
4.2
0.1447
51.186
6:46:40
-8.18
116.12
24
4.2
0.6232
94.741
16:30:31
-8.29
116.71
10
4.2
0.1301
67.643
15:12:06
-8.27
116.46
10
4.2
0.2691
59.447
08:32:43
-8.17
116.84
10
4.2
0.0960
91.051
4:22:51
-8.31
116.49
13
4.2
1.5550
65.503
1:53:17
-8.61
116.36
15
4.2
0.7044
60.974
14:04:40
-8.27
116.21
15
4.2
0.1848
91.537
17:49:28
-8.28
116.25
12
4.2
0.5169
86.932
3:13:29
-8.43
116.6
10
4.2
10.3465
45.317
10:05:01
-8.33
116.44
10
4.2
1.2814
64.249
0:18:34
-8.19
116.58
10
4.2
2.7330
57.376
16:30:31
-8.29
116.71
10
4.2
1.6323
54.230
15:12:06
-8.27
116.46
10
4.2
1.9839
70.414
54
03/09/18 19/08/201 8 09/08/201 8 05/08/201 8 05/08/201 8 02/08/201 8 01/08/201 8 01/08/201 8 28/07/201 8 1/8/18 14/09/201 8 19/09/201 8 19/08/201 8 16/08/18 09/08/18 02/08/201 8 01/08/201 8 01/08/18 28/07/201 8 01/08/18 13/08/201 8 29/07/201 8 19/08/201 8 19/08/201 8 05/08/201 8 05/08/201 8 29/07/201 8
08:32:43
-8.17
116.84
10
4.2
0.4518
64.085
9:30:31
-8.31
116.68
10
4.3
0.1282
67.643
20:55:24
-8.21
116.2
10
4.3
0.3644
40.552
23:46:32
-8.54
116.66
10
4.3
0.6154
55.115
20:35:14
-8.28
116.22
10
4.3
0.4650
33.431
19:17:43
-8.25
116.64
10
4.3
0.1763
63.295
10:44:02
-8.24
116.68
10
4.3
0.2480
67.544
10:04:01
-8.25
116.68
10
4.3
0.1145
66.951
23:19:37
-8.34
116.58
16
4.3
0.2506
52.857
17:44:02
-8.24
116.68
10
4.3
0.2556
26.038
11:57:08
-8.16
116.37
12
4.3
0.4321
69.149
01:45:06
-8.31
116.08
10
4.3
1.1747
56.098
9:30:31
-8.31
116.68
10
4.3
1.5978
54.230
22:12:13 20:55:24
-8.31 -8.21
116.64 116.2
10 10
4.3 4.3
12.2612 0.2077
55.410 95.724
19:17:43
-8.25
116.64
10
4.3
3.2626
61.222
10:44:02
-8.24
116.68
10
4.3
0.9274
60.433
10:04:01
-8.25
116.68
10
4.3
1.0457
59.447
23:19:37
-8.34
116.58
16
4.3
2.6646
57.568
17:44:02
-8.24
116.68
10
4.3
0.9549
60.532
8:42:23
-8.27
116.14
12
4.3
1.6200
97.541
0:05:30
-8.32
116.48
12
4.4
0.3120
44.830
16:23:58
-8.3
116.66
10
4.4
0.2525
62.308
4:19:36
-8.42
116.55
10
4.4
0.3439
46.098
12:25:02
-8.34
116.38
12
4.4
0.3333
36.056
12:07:05
-8.31
116.16
24
4.4
0.6186
36.878
7:05:30
-8.32
116.48
12
4.4
0.3180
56.776
55
17/08/201 8 29/07/201 8 19/08/201 8 19/08/201 8 29/07/201 8 17/08/201 8 23/08/201 8 29/07/201 8 29/07/201 8 29/07/201 8 12/08/201 8 05/08/201 8 14/09/201 8 05/08/201 8 29/07/201 8 07/08/201 8 29/07/201 8 29/07/201 8 5/8/18 21/08/201 8 24/08/201 8 23/08/201 8 29/07/201 8 29/07/201 8 29/07/201 8
22:21:10
-8.38
116.53
11
4.4
0.3113
46.131
0:05:30
-8.32
116.48
12
4.4
0.6550
65.214
16:23:58
-8.3
116.66
10
4.4
4.1364
55.410
4:19:36
-8.42
116.55
10
4.4
3.5275
51.775
7:05:30
-8.32
116.48
12
4.4
0.6561
65.312
22:21:10
-8.38
116.53
11
4.4
1.0892
56.580
18:08:36
-8.23
116.16
10
4.5
1.1690
35.152
18:03:16
-8.32
116.5
10
4.5
0.1844
46.488
3:03:24
-8.31
116.5
13
4.5
0.2360
47.172
0:43:54
-8.3
116.62
16
4.5
0.1275
58.560
18:53:54
-8.28
116.45
12
4.5
0.0371
45.800
14:53:01
-8.24
116.22
10
4.5
3.0175
37.652
20:55:05
-8.21
116.49
10
4.5
0.2251
53.739
17:42:35
-8.31
116.03
10
4.5
1.2788
33.335
3:04:29
-8.14
116.45
10
4.5
0.3719
45.220
6:01:11
-8.4
116.03
11
4.5
1.3025
25.947
2:21:50
-8.4
116.49
10
4.5
0.5458
41.292
10:03:24
-8.31
116.5
13
4.5
0.3681
73.880
21:53:01
-8.24
116.22
10
4.5
3.0449
43.742
05:28:36
-8.3
116.49
10
4.5
0.3971
57.674
01:08:36
-8.26
116.14
10
4.5
1.2835
59.151
18:08:36
-8.23
116.16
10
4.5
0.6681
97.813
18:03:16
-8.32
116.5
10
4.5
0.6290
63.394
3:03:24
-8.31
116.5
13
4.5
0.9755
64.817
0:43:54
-8.3
116.62
16
4.5
1.5119
58.818
56
12/08/201 8 05/08/18 22/12/18 29/07/201 8 07/08/201 8 29/07/201 8 29/07/201 8 07/08/18 21/08/201 8 24/08/201 8 22/12/201 8 29/08/201 8 08/08/201 8 05/08/201 8 05/08/201 8 05/08/201 8 25/08/201 8 28/07/201 8 28/07/201 8 29/07/201 8 20/08/201 8 20/08/201 8 30/08/201 8 09/10/201 8 10/10/201 8 29/08/201 8
18:53:54
-8.28
116.45
12
4.5
0.1707
70.627
1:53:45 5:58:27
-8.22 -8.63
116.3 116.3
10 10
4.5 4.5
0.0194 8.0267
86.679 65.864
3:04:29
-8.14
116.45
10
4.5
0.9755
61.984
6:01:11
-8.4
116.03
11
4.5
1.3762
101.498
2:21:50
-8.4
116.49
10
4.5
2.2261
58.068
10:03:24
-8.31
116.5
13
4.5
0.8204
64.817
13:01:11
-8.34
116.03
11
4.5
1.3645
104.182
05:28:36
-8.3
116.49
10
4.5
1.1472
65.765
01:08:36
-8.26
116.14
10
4.5
0.6899
97.912
12:56:26
-8.56
116.31
10
4.5
11.051904
66.654
23:10:34
-8.3
116.08
11
4.6
0.6180
32.193
2:10:55
-8.31
116.16
13
4.6
1.2462
30.871
15:26:52
-8.26
116.47
10
4.6
0.3148
48.346
12:35:19
-8.26
116.04
10
4.6
1.0413
37.749
12:08:41
-8.37
116.01
10
4.6
1.1608
29.168
16:27:14
-8.22
116.67
10
4.6
0.1822
67.841
22:58:47
-8.35
116.54
10
4.6
1.0314
49.716
23:15:56
-8.3
116.08
11
4.6
0.6180
41.425
7:37:35
-8.31
116.49
11
4.6
0.2596
50.010
04:44:53
-8.29
116.56
10
4.6
0.7370
62.110
05:21:31
-8.21
116.62
10
4.6
0.2822
72.988
06:10:33
-8.3
116.08
11
4.6
0.6484
57.266
15:49:17
-8.12
116.46
10
4.6
0.2700
75.168
09:16:02
-8.18
116.82
10
4.6
0.0930
88.765
23:10:34
-8.3
116.08
11
4.6
1.4864
101.195
57
29/07/18 08/08/201 8 25/08/201 8 28/07/201 8 28/07/201 8 29/07/201 8 30/08/201 8 28/07/201 8 19/08/201 8 15/08/201 8 05/08/201 8 18/10/201 8 29/07/201 8 15/08/201 8 20/08/201 8 18/09/201 8 28/07/201 8 15/08/201 8 29/07/201 8 29/07/201 8 29/07/201 8 15/08/201 8 20/08/201 8 28/07/201 8 19/08/201 8
0:37:35
-8.31
116.49
11
4.6
0.6256
79.564
2:10:55
-8.31
116.16
13
4.6
1.3984
93.409
16:27:14
-8.22
116.67
10
4.6
1.4052
62.900
22:58:47
-8.35
116.54
10
4.6
3.7031
55.803
23:15:56
-8.3
116.08
11
4.6
1.4864
46.684
7:37:35
-8.31
116.49
11
4.6
0.7736
65.136
06:10:33
-8.3
116.08
11
4.6
1.4658
101.398
22:55:39
-8.3
116.58
13
4.7
0.4843
55.446
21:44:52
-8.27
116.58
10
4.7
0.3038
56.689
9:04:10
-8.21
116.6
10
4.7
0.1682
61.025
13:12:13
-8.28
116.2
11
4.7
5.1322
33.365
1:13:00
-8.31
116.25
10
4.7
0.4417
31.149
5:55:39
-8.3
116.58
13
4.7
0.6367
59.974
16:04:10
-8.23
116.6
10
4.7
0.2666
61.124
20:25:13
-8.31
116.81
10
4.7
0.3219
77.944
14:27:16
-8.14
116.27
10
4.7
0.1482
70.315
22:55:39
-8.3
116.58
13
4.7
5.2376
60.318
9:04:10
-8.21
116.6
10
4.7
2.7381
65.172
5:55:39
-8.3
116.58
13
4.7
5.2162
60.415
5:58:47
-8.36
116.56
10
4.7
3.7188
55.902
7:18:34
-8.35
116.57
16
4.7
2.8990
57.376
16:04:10
-8.23
116.6
10
4.7
2.7612
65.271
20:25:13
-8.31
116.81
10
4.7
12.7397
49.324
23:39:21
-8.26
116.48
11
4.8
0.2682
49.129
9:50:11
-8.29
116.69
10
4.8
0.2798
65.765
58
13/08/201 8 04/08/201 8 01/11/201 8 29/07/201 8 20/08/201 8 09/09/201 8 1810/2018 28/07/201 8 19/08/201 8 13/08/201 8 09/08/18 04/08/201 8 29/07/201 8 10/08/18 20/08/201 8 29/07/201 8 19/08/201 8 19/08/201 8 29/07/201 8 20/08/201 8 29/07/201 8 19/08/201 8 19/08/201 8 29/07/201 8 08/08/18 20/08/201 8 05/08/201
1:42:24
-8.29
116.15
13
4.8
1.8366
32.971
5:44:55
-8.28
116.64
10
4.8
0.6285
61.518
1:06:23
-7.92
117.41
10
4.8
0.1174
48.547
6:39:21
-8.26
116.48
11
4.8
0.2814
65.864
06:29:55
-8.17
116.51
10
4.8
0.6211
71.701
11:53:35
-8.13
116.58
5
4.8
0.3560
78.260
08:13:01
-8.2
116.24
10
4.8
0.4677
63.887
23:39:21
-8.26
116.48
11
4.8
0.9448
69.970
9:50:11
-8.29
116.69
10
4.8
3.8602
55.016
1:42:24
-8.29
116.15
13
4.8
1.4016
95.489
20:33:04
-8.33
116.58
10
4.8
17.2923
57.083
5:44:55
-8.28
116.64
10
4.8
10.6473
58.265
6:39:21
-8.26
116.48
11
4.8
0.9169
69.970
3:33:06
-8.44
116.5
10
4.8
22.5875
54.427
06:29:55
-8.17
116.51
10
4.8
1.7919
75.961
1:44:24
-8.28
116.55
10
4.9
1.0102
53.444
23:29:56
-8.26
116.53
14
4.9
0.6267
53.948
18:18:36
-8.23
116.6
10
4.9
0.6070
61.025
6:19:37
-8.34
116.58
16
4.9
0.2590
61.393
01:18:36
-8.23
116.6
10
4.9
0.8398
70.018
1:44:24
-8.28
116.55
10
4.9
1.6170
63.740
23:29:56
-8.26
116.53
14
4.9
1.7835
67.175
18:18:36
-8.23
116.6
10
4.9
4.7587
65.172
6:19:37
-8.34
116.58
16
4.9
2.7342
57.568
9:10:54
-8.31
116.17
12
4.9
1.4179
92.383
01:18:36
-8.23
116.6
10
4.9
100.7176
65.271
20:41:27
-8.24
116.26
10
5
0.5370
38.521
59
8 05/08/201 8 05/08/201 8 02/08/201 8 27/06/201 8 29/07/201 8 6/8/18 10/8/18 09/08/201 8 09/10/201 8 10/08/201 8 02/08/201 8 29/07/201 8 06/08/18 10/08/18 08/10/18 21/12/201 8 29/07/201 8 28/07/201 8 28/07/201 8 05/08/201 8 07/09/201 8 19/08/201 8 19/08/201 8 20/08/201 8 22/08/201 8 26/08/201 8
16:49:42
-8.32
116.27
10
5
0.8339
30.865
13:07:35
-8.21
116.47
10
5
0.6886
52.364
19:07:21
-8.29
116.63
11
5
1.3706
60.213
12:37:40
-8.12
117.82
15
5
0.1809
62.859
6:15:58
-8.16
116.54
10
5
1.3992
61.716
22:50:55 22:57:38
-8.42 -8.05
116.03 116.49
23 10
5.0 5.0
7.9068 0.2108
45.512 67.149
00:14:36
-8.26
116.66
10
5
0.3740
70.117
10:50:15
-8.16
116.6
10
5
0.4040
75.961
00:57:22
-8.43
116.76
10
5
0.5350
70.414
19:07:21
-8.29
116.63
11
5
8.9570
58.052
6:15:58
-8.16
116.54
10
5
7.1785
74.375
22:50:55 22:57:38 00:57:22
-8.42 -8.05 -8.43
116.03 116.49 116.76
23 10 10
5 5 5
0.3326 0.5962 26.3263
102.026 87.672 37.845
17:43:56
-8.35
116.74
10
5
21.2242
46.293
2:20:09
-8.24
116.5
11
5.1
0.8822
50.122
23:05:40
-8.34
116.51
13
5.1
0.5120
46.195
22:56:46
-8.32
116.52
11
5.1
0.6225
48.248
12:42:32
-8.26
116.21
10
5.1
1.7009
35.344
17:14:36
-8.28
116.67
11
5.1
0.4749
64.544
21:50:01
-8.66
116.41
10
5.1
0.7269
30.487
3:21:08
-8.33
116.92
10
5.1
1.0906
85.288
7:31:34
-8.23
116.47
10
5.1
0.8996
63.097
21:50:22
-8.15
116.82
10
5.1
0.2206
90.454
10:54:05
-8.36
116.95
10
5.1
0.5570
86.182
60
31/08/201 8 09/06/201 8 10/08/201 8 29/07/201 8 28/07/201 8 19/08/201 8 20/08/201 8 22/08/201 8 26/08/201 8 31/08/201 8 06/09/18 08/10/18 29/07/201 8 06/08/201 8 10/09/201 8 20/08/201 8 21/08/201 8 29/07/201 8 06/08/201 8 20/08/201 8 21/08/201 8 06/08/201 8 6/8/18 20/08/18 09/02/201 8 09/11/201 8 12/06/201
09:37:15
-8.37
116.06
10
5.1
5.2721
50.304
13:43:09
-8.29
116.99
10
5.1
0.3480
94.033
00:57:21
-8.41
116.81
10
5.1
0.5350
70.414
2:20:09
-8.24
116.5
11
5.1
3.7907
69.127
23:05:40
-8.34
116.51
13
5.1
2.8178
61.588
3:21:08
-8.33
116.92
10
5.1
23.9990
45.805
7:31:34
-8.23
116.47
10
5.1
2.0290
72.295
21:50:22
-8.15
116.82
10
5.1
0.6291
65.666
10:54:05
-8.36
116.95
10
5.1
19.8418
42.980
09:37:15
-8.37
116.06
10
5.1
2.6440
98.907
13:43:09 00:57:21
-8.29 -8.41
116.99 116.81
10 10
5.1 5.1
5.8393 26.3263
51.382 37.845
1:50:32
-8.22
116.46
10
5.2
1.8436
53.301
18:21:19
-8.36
116.39
10
5.2
3.6087
34.673
20:22:28
-8.32
116.47
10
5.2
0.8793
43.953
01:23:25
-8.42
116.93
13
5.2
0.9376
81.345
08:09:24
-8.44
116.94
10
5.2
0.7650
80.821
1:50:32
-8.22
116.46
10
5.2
2.9250
74.574
18:21:19
-8.36
116.39
10
5.2
4.2878
69.226
01:23:25
-8.42
116.93
13
5.2
100.6088
37.149
08:09:24
-8.44
116.94
10
5.2
74.9808
34.291
0:28:20
-8.36
116.26
27
5.3
3.6110
31.782
7:28:19 08:30:32
-8.37 -8.35
116.26 116.73
20 13
5.3 5.3
3.7095 0.6231
31.212 69.330
07:15:47
-8.15
116.37
16
5.3
0.7520
71.805
03:22:26
-8.41
116.52
10
5.3
0.6740
48.737
08:02:47
-8.5
116.06
10
5.3
10.2160
39.874
61
8 20/08/18 02/09/18 06/08/201 8 19/08/201 8 19/08/201 8 09/11/201 8 06/08/201 8 19/08/201 8 19/08/18 06/12/201 8 7/8/18 19/8/18 26/08/18 07/08/18 19/08/18 26/08/201 8 26/08/18 19/08/18 26/08/18 26/08/18 05/08/201 8 29/07/201 8 5/8/18 29/07/201 8 09/08/18 08/09/201 8 08/09/18 19/08/201 8 19/08/18 09/08/18 19/08/201 8
08:30:32 07:15:47
-8.35 -8.15
116.73 116.37
13 16
5.3 5.3
5.2383 2.2886
48.958 88.262
15:50:55
-8.37
116.03
10
5.4
7.2345
27.391
4:06:13
-8.29
116.62
10
5.4
1.4001
59.151
11:06:13
-8.29
116.62
10
5.4
4.7870
58.856
03:22:27
-8.32
116.48
10
5.4
0.8110
55.016
15:50:55
-8.37
116.03
10
5.4
10.8368
101.693
4:06:13
-8.29
116.62
10
5.4
17.0858
58.363
11:06:13
-8.29
116.62
10
5.4
56.3890
58.166
1:02:48
-8.5
116.06
10
5.5
13.6672
21.771
1:21:19 23:37:38 1:33:17 1:21:19 23:37:38
-8.18 -8.25 -8.55 -8.18 -8.25
116.29 116.84 116.93 116.29 116.84
10 10 11 10 10
5.5 5.5 5.5 5.5 5.5
3.5494 1.0973 3.0526 4.2834 54.5765
45.220 83.997 75.703 89.659 54.918
01:33:17
-8.55
116.93
11
5.5
69.8036
25.676
01:33:20 22:16:37 01:33:18 01:33:18
-8.53 -8.35 -8.47 -8.47
116.93 116.53 116.93 116.93
30 10 10 10
5.5 5.6 5.6 5.6
69.8022 10.2468 3.0526 69.8036
47.634 54.132 78.440 32.098
19:49:52
-8.36
116.16
10
5.7
4.9380
60.213
8:50:32
-8.22
116.46
10
5.7
1.9098
57.684
19:49:52
-8.36
116.15
11
5.7
7.0250
40.359
8:50:32
-8.22
116.46
10
5.7
3.1954
73.385
12:25:32
-8.44
116.21
14
5.9
16.7159
81.117
12:25:33
-8.4
116.21
12
6
33.5880
22.177
12:25:33
-8.4
116.21
12
6
12.4090
83.863
4:10:23
-8.44
116.59
18
6.2
5.1380
49.387
4:10:23 12:25:32
-8.44 -8.36
116.59 116.22
18 12
6.2 6.2
53.2410 16.7159
50.114 81.686
11:10:23
-8.45
116.57
14
6.3
5.1380
48.074
62
19/08/18 19/08/201 8 19/08/18 28/07/201 8 29/07/201 8 28/07/18 29/07/201 8 19/08/18 19/08/18 5/8/18 05/08/18
11:10:23
-8.45
116.57
14
6.3
4.7870
45.793
11:10:23
-8.45
116.57
14
6.3
53.2420
48.820
11:10:23
-8.45
116.57
14
6.3
72.5724
46.650
10:47:39
-8.4
116.5
24
6.4
7.9980
47.387
5:47:39
-8.4
116.5
24
6.4
8.4042
60.729
10:47:39
-8.4
116.5
24
6.4
40.7540
61.322
5:47:39
-8.4
116.5
24
6.4
51.4665
60.834
11:10:22 11:10:22 18:46:35 18:46:35
-8.24 -8.24 -8.37 -8.37
116.66 116.66 116.48 116.48
10 10 15 15
6.5 6.5 7.0 7
7.2381 72.5724 53.7961 23.7752
60.433 56.098 37.652 52.872
Lampiran C Hasil Koefisien tiap Residual Sebelum proses residual Std. Lower Upper Parameter Estimate Error Bound Bound a -2.525 0.452 -3.413 -1.636 b 0.712 0.034 0.645 0.779 -0.838000 c -0.395 0.226 0.049 Dependent variable: Percepatan Tanah a. R squared = 1 - (Residual Sum of Squares) /(Corrected Sum of Squares) =.450.
Lampiran D Nilai PGA hasil rekaman akselerograf dari Shakemap untuk magnitudo lebih dari 4,5
63
64
Lampiran E Nilai PGA hasil rekaman akselerograf dari SIGMA untuk magnitudo kurang dari 4,5
65
Lampiran F Riwayat Hidup
RIWAYAT HIDUP Nama
: Hasti Wahyu Cahyaningrum
Tempat Tanggal Lahir
: Jakarta, 23 Juni 1997
Jenis Kelamin
: Perempuan
Agama
: Islam
No Telp
: 08157696749
Email
: [email protected]
Alamat
: Jalan Belimbing Gang Kramat RT 011 RW 007 No.80 Jagakarsa, Jakarta selatan
RIWAYAT PENDIDIKAN Tahun 2003 – 2009 2009 – 2012 2012 – 2015 2015 – 2019
Pendidikan SDN Jagakarsa 02 Pagi SMP Negeri 175 Jakarta SMA Suluh Jakarta Jurusan Fisika FMIPA Universitas Jenderal Soedirman
RIWAYAT ORGANISASI Tahun 2016 2017 2017 2017 2018
Organisasi Himpunan Mahasiswa Fisika FMIPA Unsoed Himpunan Mahasiswa Fisika FMIPA Unsoed Komunitas Geofisika Unsoed Ikatan Himpunan Mahasiswa Fisika Indonesia Wilayah 3 Himpunan Mahasiswa Fisika FMIPA Unsoed
66
Jabatan Staf Humas dan Advokasi Kepala Departemen Hubungan Masyarakat Staff Divisi Akademik Kegiatan Staff Divisi Kekeluargaan Dewan Pertimbangan Organisasi