Pemodelan Parameter Dinamis Tanah (VS, G) Berdasarkan Metode
Pemodelan Parameter Dinamis Tanah (VS, G) Berdasarkan Metode
Pemodelan Parameter Dinamis Tanah (VS, G) Berdasarkan Metode
Dosen Pembimbing I
Dr. Dwa Desa Warnana
NIP. 19760123 200003 1 001
Dosen Pembimbing II
Dr. Ir. Amien Widodo, M.Si.
NIP. 19591010 198803 1 001
Dosen Pembimbing I
Dr. Dwa Desa Warnana
NIP. 19760123 200003 1 001
Dosen Pembimbing II
Dr. Ir. Amien Widodo, M.Si.
NIP. 19591010 198803 1 002
i
HALAMAN INI SENGAJA DIKOSONGKAN
ii
UNDERGRADUATE THESIS - RF141501
Dosen Pembimbing I
Dr. Dwa Desa Warnana
NIP. 19760123 200003 1 001
Dosen Pembimbing II
Dr. Ir. Amien Widodo, M.Si.
NIP. 19591010 198803 1 002
iii
HALAMAN INI SENGAJA DIKOSONGKAN
iv
TUGAS AKHIR
v
HALAMAN INI SENGAJA DIKOSONGKAN
vi
PERNYATAAN KEASLIAN TUGAS AKHIR
vii
HALAMAN INI SENGAJA DIKOSONGKAN
viii
Pemodelan Parameter Dinamis Tanah (Vs, G0) Berdasarkan
Data Metode Multichannel Analysis Of Surface Wave (MASW)
Sebagai Evaluasi Tapak Lokal Surabaya
Nama Mahasiswa : Novien Ghoziana Indanartha
NRP : 03411440000065
Departemen : Teknik Geofisika FTSLK-ITS
Dosen Pembimbing : Dr. Dwa Desa Warnana
NIP. 19760123 200003 1 001
Dr. Ir. Amien Widodo, M.Si.
NIP. 19591010 198803 1 002
ABSTRAK
ix
HALAMAN INI SENGAJA DIKOSONGKAN
x
Dynamic Soil Parameter (Vs, G0) Modelling Based On
Multichannel Analysis Of Surface Wave (MASW) Method
For Evaluation of Site Local Surabaya
ABSTRACT
xi
HALAMAN INI SENGAJA DIKOSONGKAN
xii
KATA PENGANTAR
Puji syukur kehadirat Allah SWT karena dengan rahmat dan ridho-Nya
Sehingga laporan Tugas Akhir dengan judul “Pemodelan Parameter Dinamis
Tanah (Vs, G) Berdasarkan Metode Multichannel Analysis Of Surface Wave
(MASW) Sebagai Evaluasi Tapak Lokal Surabaya” telah terselesaikan.
Pelaksanaan penelitian dan penyusunan laporan Tugas Akhir ini tidak terlepas
dari bantuan,bimbingan dan dukungan dari berbagai pihak. Pada kesempatan ini
saya mengucapkan terimakasih kepada:
1. Allah SWT atas rahmat dan ridho-Nya sehingga Tugas Akhir ini dapat
berjalan dan terselesaikan dengan baik
2. Ayah dan Mama atas segala doa, pengorbanan, kesabaran, dan kasih
sayang yang telah diberikan tanpa rasa lelah.
3. Riza Ghozlina Indanartha, Nanang Shodiq, dan Nisrina Farah
Indanartha atas dukungan dan doa yang selalu ditujukan kepada penulis
selama hidupnya.
4. Semua keluarga atas motivasi dan dukungan yang sangat besar.
5. Dr. Dwa Desa Warnana dan Dr. Ir. Amien Widodo, M.Si. selaku
pembimbing yang telah memberikan bimbingan dan arahan dalam
pengerjaan Tugas Akhir sehingga laporan ini dapat terselesaikan.
6. Firman Syaifuddin, S.Si., M.T. selaku dosen pembimbing akademik
yang telah memberikan bimbingan dan motivasi selama penulis
menjalani proses perkuliahan
7. Seluruh Tim MASW (Ardianto, Moh. Fani Fatkhur Rohman, Husnia
Nur Annisa, Mochammad Novritsa Zulfikar) yang selalu menemani dan
mendukung penulis selama pengerjaan Tugas Akhir ini.
8. Seluruh dosen dan pegawai administrasi Departemen Teknik Geofisika
ITS yang telah banyak memberikan ilmu selama penulis melakukan
studi di Departemen Teknik Geofisika ITS.
9. Teman-teman Teknik Geofisika ITS angkatan 2014 atas dukungannya.
10. Semua pihak yang tidak dapat dituliskan satu per satu oleh penulis,
terima kasih banyak atas doa dan dukungannya.
Semoga Allah membalas semua kebaikan semua pihak. Penulis
menyadari tentunya penulisan Tugas Akhir ini masih banyak kekurangan. Oleh
karena itu, kritik dan saran yang membangun sangat diharapkan. Semoga Tugas
Akhir ini membawa manfaat bagi penulis pribadi maupun bagi pembaca.
Surabaya, 08 Juni 2018
xiv
DAFTAR ISI
xvi
DAFTAR GAMBAR
Gambar 1.1 Patahan Waru dan Patahan Surabaya (Tim Pusgen, 2017) ............. 1
Gambar 1.2 Peta Vs30 Kota Surabaya dengan Metode Mikrotremor
(Mufida,dkk. 2013). ............................................................................................ 4
Gambar 1.3 Peta Amplifikasi Surabaya Metode Mikrotremor (Rahmaningtyas
dkk,2017 ............................................................................................................. 5
Gambar 2.1 Kolom stratigrafi Sesar Kendeng (De genevraye, 1972)................. 9
Gambar 2.2 Peta geologi Kota Surabaya oleh Sukardi dkk, 1992 dimodifikasi11
Gambar 2.3 Sesar di Pulau Jawa (Tim Pusgen,2017) ....................................... 12
Gambar 2.4 Ilustrasi desain akuisisi MASW aktif (Park,2007) ........................ 16
Gambar 2.5 Pola gerakan partikel gelombang Rayleigh (Lowrie,2007) ........... 18
Gambar 2.6 Trace gelombang Rayleigh pada saat keadaan non-dispersif dan pada
keadaan dispersif (Foti et al., 2015) .................................................................. 18
Gambar 2.7 Contoh kurva dispersi Gelombang Rayleigh dengan normal dispersi
(Dal Moro et al., 2014) ..................................................................................... 19
Gambar 2.8 Model lapisan bumi dengan parameter kecepatan gelombang geser
(vs ), kecepatan gelombang kompresional (vp), densitas (ρ), dan ketebalan (h).
(Xia dkk., 1999) ................................................................................................ 23
Gambar 3.1 Titik lokasi pengukuran di Surabaya ............................................. 26
Gambar 3.2 Desain akuisisi MASW ................................................................. 27
Gambar 3.3 Diagram alir penelitian seacar umum ........................................... 28
Gambar 3.4 Diagram alir pengambilan data ..................................................... 29
Gambar 3.5 Diagram alir pengolahan data MASW .......................................... 30
Gambar 3. 6 Diagram alir pemodelan 3D ......................................................... 31
(b) Gambar 3.7 Hasil Inversi (a) Profil Vs 1D (b) Stratigrafi dan nilai error ... 35
Gambar 3.8 Hasil nilai vs dan ketebalan masing-masing lapisan ..................... 35
Gambar 4.1 Peta persebaran nilai Vs30 di Surabaya ....................................... 47
Gambar 4.2 Penampang 3D Vs dari arah Barat ke Utara (W-N) ...................... 48
Gambar 4.3 Penampang 3D Vs dari arah Selatan ke Barat (S-W) .................... 48
Gambar 4.4 Penampang 3D Vs dari arah Timur ke Selatan (E-S) .................... 49
Gambar 4.5 Penampang 3D Vs dari arah Utara ke Timur (N-E) ...................... 49
Gambar 4.6 Penampang 3D G0 dari arah Barat ke Utara (W-N) ..................... 50
Gambar 4.7 Penampang 3D G0 dari arah Timur ke Selatan ............................. 51
Gambar 4.8 Intersection penampang 2D (1-1’ = Inline 20 ; 2-2’ = Inline 26 ; 3-
3’ = Inline 65 ; 4-4’ = Xline 70 ; 5-5’ =Xline 30 ) ........................................... 51
Gambar 4.9 Inline 20 (intersection 1-1’) pada model Vs.................................. 52
xvii
Gambar 4.10 Inline 26 (Intersection 2-2’) pada model Vs ............................... 52
Gambar 4.11 Inline 65 (Intersection 3-3’) pada model Vs ............................... 53
Gambar 4.12 Xline 70 (Intersection 4-4’) pada model Vs ............................... 53
Gambar 4.13 Xline 30 (Intersection 5-5’) pada model Vs ............................... 54
Gambar 4.14 Crossection pada Xline 60 dan Inline 60 pada model Vs ........... 55
Gambar 4.15 Inline 20 (Intersection 1-1’) pada model G0 .............................. 55
Gambar 4.16 Inline 26 (Intersection 2-2’) pada model G0 .............................. 56
Gambar 4.17 Inline 65 (Intersection 3-3’) pada model G0 .............................. 56
Gambar 4.18 Xline 70 (Intersection 4-4’) pada model G0 ............................... 57
Gambar 4.19 Xline 30 (Intersection 5-5’) pada model G0 ............................... 57
Gambar 4.20 Crossection pada Xline 70 dan Inline 65 pada model G0 ........... 58
xviii
DAFTAR TABEL
Tabel 2.1 Klasifikasi nilai Vs30 menurut Eurocode 8 (Cen, 2004 ; SIST EN, 2005
; Carvalho,2008 ;Mufida dkk., 2013).............................................................. 144
Tabel 4.1 Hasil Vs30 Pengukuran MASW Surabaya ....................................... 37
xix
HALAMAN INI SENGAJA DIKOSONGKAN
xx
BAB I
PENDAHULUAN
1.1 Latar Belakang
Surabaya merupakan Ibukota Provinsi Jawa Timur, Indonesia. Kota ini
memiliki kepadatan penduduk yang cukup tinggi dengan jumlah penduduk
sebesar 2.848.583 jiwa pada tahun 2015 dengan luas wilayah 350,54 km2
(Badan Pusat Statistik, 2015). Kota ini dilewati oleh sesar aktif Baribis-Kendeng
fold-thrust zone yang memanjang mengarah Barat ke Timur di bagian Utara
Pulau Jawa yaitu segmen Surabaya sepanjang 25 Km dan segmen Waru
sepanjang 64 Km dengan pergerakan 0,05 mm/tahun dan Magnitudo maksimum
6,5 yang ditunjukkan pada gambar 1.1 (Pusgen, 2017). Oleh karena itu, Surabaya
merupakan daerah yang rawan akan terjadinya goncangan gempa bumi.
Gambar 1.1 Patahan Waru dan patahan Surabaya (Tim Pusgen, 2017)
Berdasarkan sejarah gempa bumi yang terjadi di sekitar Surabaya mulai
tahun 1800 hingga tahun 2016, rata-rata memiliki intensitas sebesar VII MMI.
Intensitas tersebut, setiap orang diharuskan untuk lari ke luar dan geataran
dirasakan di dalam kendaraan yang bergerak. Penelitian PGA oleh (pusgen,
2017) menyatakan bahwa Surabaya berada pada bahaya gempa bumi kelas 4
(dirasakan kuat oleh penduduk dan gedung mengalami guncangan).
1
Tingkat kerusakan dan bahaya gempabumi sangat dipengaruhi oleh kondisi
geologi lokal atau efek tapak lokal. Beberapa kasus gempabumi yang telah terjadi
menunjukkan bahwa kerusakan lebih parah terjadi pada dataran alluvial
dibandingkan dengan daerah perbukitan (Nakamura, dkk., 2000). Menurut peta
geologi regional lembar Surabaya oleh Soekardi (1992), 80% wilayah Surabaya
merupakan endapan aluvial, dan sisanya merupakan tanah hasil pelapukan
batuan tersier/tua. Batuan sedimen lunak mampu memperkuat getaran (Tuladhar,
2002). Oleh karena itu, rata-rata kerusakan gempa bumi lebih parah
dibandingkan dengan daerah dengan permukaan dekat yang terbentuk oleh
batuan keras (Bahri, 2016). Amplifikasi merupakan peristiwa perbesaran
gelombang akibat gelombang gempa bumi yang merambat melalui lapisan tanah
yang lebih lunak dibandingkan lapisan awal yang dilaluinya. Hal ini yang
meningkatkan potensi kerusakan akibat gempa bumi (Aisyah dkk., 2011). Oleh
sebab itu, diperlukan evaluasi kondisi lokal atau tapak lokal surabaya terhadap
efek gelombang seismik.
Evaluasi tapak lokal merupakan upaya pengurangan risiko bencana pada
suatu daerah dengan menganalisis seluruh potensi bahaya secara lengkap.
Pemahaman terhadap potensi bahaya dapat meningkatkan kapasitas persiapan
menghadapi risiko yang mungkin ditimbulkan. Oleh karena itu, kajian untuk
mengetahui potensi bahaya di suatu daerah sangat penting dilakukan. (Sunardi et
al., 2012).
Evaluasi dapat dilakukan dengan menanalisa perilaku dimanis tanah pada
bawah permukaan dangkal. Salah satu upaya yang dapat dilakukan adalah
melakukan pengukuran parameter dinamis tanah. Parameter ini berupa kecepatan
gelombang geser (Vs) an modulus geser (G), merupakan hal penting dalam
analisis respon dinamis tanah terhadap gelombang gempa bumi (Luna, 2000).
Penentuan parameter dinamis tanah dapat dilakukan dari pengujian
lapangan dengan metode cross-hole, down-hole, spectral analysis surface wave
(SASW), seismic cone penetration test, dan multichannel analysis surface wave
(MASW). Metode MASW merupakan metode seismik non-destructive yang
merekan perembatan gelombang permukaan (gelombang rayleigh) (Ariestianty
dkk, 2009).
Penelitian dengan menggunakan metode MASW (Multichannel Analysis
Surface Wave), yang merupakan salah satu metode seismik dengan ketepatan
dan resolusi yang tinggi pada struktur dekat permukaan untuk pemodelan
struktur geologi bawah permukaan bumi (Park, 2000; Kanli, 2006). Pemetaan
nilai kecepatan gelombang geser (Vs) konvensional dengan metode geofisika
borehole memiliki kekurangan pada segi biaya dan waktu. Sedangkan Penelitian
Vs menggunakan mikrotremor memiliki kekurangan karena data yang diperoleh
berupa gelombang alami bumi sehingga rentan terpengaruh oleh noise atau
gelombang yang tidak diinginkan baik akibat aktivitas alam maupun manusia
(Erik, 2009).
2
Aplikasi MASW telah diterapkan oleh banyak peneliti salah satunya
Tokeshi dkk. (2013) melakukan analisis MASW untuk karakterisasi dan
kepadatan tanah sebagai informasi penting dalam bidang geoteknik dalam
mendesain bangunan tahan guncangan. Penentuan modulus geser tanah
menggunakan metode MASW telah dilakukan oleh Arienstianty, dkk tahun
2009.
4
Gambar 1.3 Peta Amplifikasi Surabaya Metode Mikrotremor (Rahmaningtyas
dkk,2017
5
HALAMAN INI SENGAJA DIKOSONGKAN
6
BAB II
TINJAUAN PUSTAKA
2.1 Kondisi Regional Surabaya
Surabaya adalah kota terbesar kedua di Indonesia setelah Jakarta. Surabaya
terletak di tepi pantai Utara Pulau Jawa dan berhadapan dengan Selat
Madura serta Laut Jawa dengan koordinat 7°16′LU 112°43′BT. Luas wilayah
Kota Surabaya 33.048 Ha. Wilayah Surabaya secara umum terbagim menjadi 4
wilayah yaitu Surabaya Pusat, Surabaya Timur, Surabaya Barat, Surabaya Utara,
dan Surabaya Selatan. Secara administratif, Kota Surabaya terbagi dalam 31
Kecamatan, 163 Kelurahan, 1.363 Rukun Warga, dan 8.909 Rukun Tetangga.
Kota Surabaya memiliki ketinggian tanah antara 0 – 20 meter di atas permukaan
laut, sedangkan pada daerah pantai ketinggiannya berkisar antara 1–3 meter
diatas permukaan laut. Sebagian besar Kota Surabaya memiliki ketinggian tanah
antara 0 – 10 meter (80,72 % atau sekitar 26.345,19 Ha) yang menyebar di bagian
Timur, Utara, Selatan dan pusat kota. (Bahri, 2012).
2.1.1 Stratigrafi Surabaya
Daerah penelitian merupakan daerah yang dilewati oleh anticlinal Kendeng
atau Sesar Kendeng yang membentang sejauh 250 km dengan lebar 40 km.
Pringgoprawiro (1983) membagi morfologi Zona Kendeng menjadi tiga satuan
yang masing-masing membentang dari Barat ke Timur, yaitu:
1. Satuan morfologi perbukitan bergelombang, ditunjukan oleh jajaran bukit-
bukit rendah dengan ketinggian antara 50-100 meter di atas permukaan laut
yang mencerminkan lipatan batuan sedimen. Satuan ini nyaris secara
keseluruhan disusun oleh litologi napal abu-abu.
2. Satuan morfologi perbukitan terjal, yang merupakan inti Pegunungan
Kendeng dengan ketiggian rata-rata 350 meteer di atas permukaan laut, tipe
genetic sungainya adalah tipe konsekuen, subsekuen, dan insekuen. Litologi
yang menyusun satuan ini, sebagian besar adalah batugamping dan batupasir.
3. Satuan morfologi dataran rendah, yang disusun oleh endapan alluvial yang
terdapat di Ngawi (Bengawan Solo) dan dataran Sungai Brantas di Timur.
Zona Kendeng merupakan bagian tenah dari Cekungan Jawa Timur.
Sebagian besar litologinya menunjukkan pengaruh lingkungan laut dalam.
Menurut Pringgoprawiro (1983), stratigrafi Zona Kendeng dibagi kedalam unit-
unit sebagai berikut:
1. Formasi Pelang terdiri dari napal abu-abu yang massif sampai belapis yang
kaya fosil dan batulempung abu-abu dengan sisipan bagu gamping bioklastik.
Lapisan ini diendapkan pada lingkungan neritic dan berumur Oligosen akhir-
Miosen awal.
7
2. Formasi Kerek terdiri dari endapan turbidit dengan ketebalan 800 m,
sebagian besar terbentuk oleh lapisan yang menghalus dan menipis ke atas
dengan tipe struktur sedimen arus densitas. Litologinya terdiri atas batupasir
tufaan, batulempung, napal, dan batugamping. Formasi ini berumur Miosen
awal-Miosen akhir.
3. Formasi Kalibeng (Kalibeng Bawah) terdiri atas napal abu-abu kehijauan
kaya fosil dengan sisipan tuf berlapis tipis. Sedimen ini diendapkan pada
lingkungan batial. Bagian atas dari Formasi Kalibeng (anggota atasangin)
terdiri atas perlapisan batu pasir tufaan berukuran halus-kasar, tuff putih, dan
breksi volkanik. Sedimen ini diendapkan oleh mekanisme turbidit. Formasi
ini berumur Miosen akhir-Pliosen.
4. Formasi Sonde (Kalibeng atas) dibedakan atas Formasi Sonde bagian atas
dan Formasi Sonde bagian bawah. Formasi Sonde bagian bawah (Anggota
Klitik) didominasi oleh perlapisan napal pasiran, batupasir gampingan, dan
tuff. Formasi Sonde bagian atas terdiri atas batugamping mengandung
Balanus dan Grainstone. Formasi ini diendapkan di lingkungan laut dangkal
dan berumur Pliosen.
5. Formasi Pucangan terdiri atas batupasir-konglomeratan, batupasir, batupasir
tufaan, dan lempung hitam yang mengandung moluska air tawar. Di Zona
Kendeng bagian Barat dan tengah, Formasi Pucangan berkembang sebagai
fasies daratan. Sedangkan di bagian Timur Zona Kendeng, Formasi Pucangan
merupakan endapan laut dangkal. Formasi ini berumur Pliosen akir -
Pleistosen awal.
6. Formasi Kabuh terdiri dari perlapisan batupasir kasar dengan perlapisan
silang-siur, fosil vertebrata, lensa konglomerat, dan tuff. Di Zona Kendeng
bagian Barat dan tengah, Formasi Kabuh diendapkan pada lingkungan darat,
sedangkan di Zona Kendeng bagian Timur Formasi Kabuh mempunyai fasies
yang berbeda-beda, faies darat berangsur-angsur berubah menjadi fasies laut
yang makin ke atas berubah ke batuan volkanik yang diendapkan pada
lingkungan pantai. Formasi ini berumur Pleistosen.
7. Formasi Notopuro terdiri dari endapan lahar, tuff, dan batu pasir tuffaan
berumur Pleistosen yang diendapkan pada lingkungan darat.
Gambar 2.1 merupakan gambar kolom strategrafi pada zona kendeng yang
menjelaskan ururtan formasi beserta umur batuannya.
8
Gambar 2.1 Kolom stratigrafi Sesar Kendeng (De genevraye, 1972)
10
Gambar 2.2 Peta geologi Kota Surabaya oleh Sukardi dkk, 1992 dimodifikasi Novien, 2018
11
Berdasarkan pada gambar 2.3, Kota Surabaya dilewati oleh jalur
baribis-kendeng fold-thrust zone yang memanjang mengarah Barat ke Timur di
Pulau Jawa. Sesar ini merupakan sesar aktif yang memanjang dari Jawa Tengah
bagian Utara mengarah menuju Jawa Timur bagian Utara. Segment dari jalur
sesar tersebut yang melewati Surabaya adalah segment Surabya sepanjang 25
Km dan segment Waru sepanjang 64 Km dengan pergerakan 0,05 mm/tahun dan
Magnitudo maksimum 6,5 seperti yang ditunjukkan pada gambar 1.1 (Pusgen,
2017). Salah satu fenomena geologi yang menarik di Surabaya adalah
pembelokan sungai pada area antara dua antiklin yang tampak memiliki
pergeseran lokasi. Hal tersebut menjadi suatu indikasi adanya sesar yang
mengakibatkan sungai bergeser secara tidak wajar. (Riyantiyo N,D dkk.2017).
Antiklin ditunjukkan pada gambar 2.1 yang memanjang dari Barat ke Timur.
Tipe Vs30
Uraian Gambaran Stratigrafi
Tanah (m/s)
14
2.2.2 Modulus Geser (G0)
Modulus geser atau modulus rigidity adalah parameter yang dihitung
berdasarkan Vs menggunakan hubungan elastis sederhana, yaitu :
G0 = 𝜌×Vs2
Kepadatan massa sering diperkirakan atau diukur dengan sampel di bawah
permukaan atau menggunakan korelasi. Modulus geser digunakan untuk
melakukan pemodelan tanah yang lebih maju, dan respon dinamik interaksi
struktur tanah. Modulus geser pada tingkat lowstrain yang diukur dengan teknik
geofisika akan memberikan parameter elastis untuk analisis pondasi mesin atau
teknik gempa. Kegunaan penting dari parameter ini adalah untuk menghitung
frekuensi resonansi dan amplitudo getaran pada pondasi, interaksi tanah-struktur
akibat beban dinamis, menentukan koefisien kekakuan tanah dibawah pondasi
dan analisis perambatan gelombang. (Luna, 2000).
15
(Vs) untuk investigasi struktur dekat permukaan secara efektif, murah, efisien
dan mudah dalam pengolahannya. MASW memberikan informasi tentang
kecepatan gelombang geser (Vs) 1D. Biasanya untuk pengambilan data
digunakan frekuensi sumber 3-30 Hz dan menggunakan multichannel perekam
(geophone) yang disusun memanjang. Perekam yang digunakan bisa sejumlah
12 atau 24 dengan jarak antar perekam sama. (Park dkk., 2007).
Metode MASW sangat popular dalam identifikasi zona-zona rawan
guncangan gempa bumi, dan diklasifikasi berdasarkan site class yang mengacu
pada nilai kecepatan gelombang shear (Vs30) oleh Eurocode 8. MASW dapat
dilakukan dalam 3 tahap yaitu (1) akusisi; (2) ekstraksi kurva dispersi; (3) inversi
kurva dispersi. Metode ini menghasilkan profil kecepatan gelombang geser (Vs)
1D. Dalam akuisisi nya, penggunaan alat per titik dapat ditunjukkan pada gambar
2.4 yang menggambarkan tentang desain akuisisi MASW aktif. (Park dkk, 2007).
17
Gambar 2.5 Pola gerakan partikel gelombang Rayleigh (Lowrie,2007)
𝐹𝐹𝑇
U𝑗 (t) → Ū𝑗 (ω) j = 1,2,3, … . , n
Ūj (ω)
Ūj,norm (ω) = = 𝑝𝑗 (ω)
|Ūj(ω)|
21
𝑉𝑅,𝑡 ∶ 𝑃𝑒𝑛𝑔𝑢𝑗𝑖𝑎𝑛 𝑘𝑒𝑐𝑒𝑝𝑎𝑡𝑎𝑛 𝑔𝑒𝑙𝑜𝑚𝑏𝑎𝑛𝑔 𝑟𝑎𝑦𝑙𝑒𝑖𝑔ℎ
22
2.3.4 Analisa Inversi Gelombang Rayleigh
Tujuan utama dari kegiatan eksplorasi geofisika adalah untuk membuat
model bawah permukaan bumi dengan menggunakan data lapangan yang diukur
pada permukaan bumi. Dari model bawah permukaan yang diperoleh,
menggambarkan sifat fisis dari batuan yang diteliti. Sehingga pada akhirnya kita
dapat mengkarakterisasi kondisi geologi bawah permukaan. Model yang kita
dapat adalah representasi dari keadaan geologi bawah permukaan yang memiliki
besaran fisis dan geometri tertentu (Grandis, 2009). Pemodelan inversi adalah
pencarian parameter fisis bumi berdasarkan “data” yang kita peroleh dari akusisi.
Pemodelan inversi biasa juga disebut sebagai tahapan data fitting karena dalam
prosesnya kita mencari nilai parameter yang terbaik dan cocok dengan data
pengamatan (data observasi).
Karakteristik dispersi dari model tanah berlapis bergantung pada beberapa
kelompok sifat tanah, seperti jumlah dan ketebalan lapisan tanah dan sifat elastis
setiap lapisan. Gambar 2.8 menunjukkan contoh umum dari model tanah
berlapis. Sifat elastis setiap lapisan umumnya ditunjukkan oleh sifat fisik yang
mudah diukur, yaitu kecepatan gelombang kompresional (Vp), kecepatan
gelombang geser (Vs) dan densitas (𝜌) (Xia dkk., 1999 ). Sifat-sifat di dalam
setiap lapisan diambil sebagai konstan dan lapisan terakhir diasumsikan setengah
ruang (Xia dkk., 1999). Untuk model tanah berlapis, profil kecepatan gelombang
geser memiliki efek dominan pada kurva dispersi mode fundamental, diikuti oleh
ketebalan lapisan sehingga merupaka parameter yang signifikan (Xia dkk.,
1999).
Gambar 2.8 Model lapisan bumi dengan parameter kecepatan gelombang geser
(vs ), kecepatan gelombang kompresional (vp), densitas (ρ), dan ketebalan (h).
(Xia dkk., 1999)
23
Gambar 2.8 merupakan model matematis untuk material elastis berlapis
digunakan untuk menentukan kurva dispersi teoritis berdasarkan seperangkat
parameter masukan yang diberikan yaitu kecepatan gelombang geser (vs ),
kecepatan gelombang kompresional (vp), densitas (ρ), dan ketebalan (h). Profil
Vs dihitung dengan menggunakan inversi berbasis iterasi yang memerlukan data
inversi dan estimasi nilai poission ratio dan densitas. Set parameter yang berbeda
dimasukkan ke dalam model dengan cara yang iteratif untuk mencari kurva
dispersi teoretis yang paling sesuai dengan karakteristik dispersi yang diamati.
Dengan demikian, masalah memperoleh profil kecepatan gelombang geser yang
dapat diterima dapat diidentifikasi sebagai masalah optimasi multi parameter
dimana tujuannya adalah untuk meminimalkan perbedaan antara kurva dispersi
eksperimental dan teoritis. Kesalahan root mean square antara dua kumpulan
data biasanya digunakan sebagai indikator perbedaan antara kurva penyebaran
teoritis dan eksperimental. (Xia dkk., 1999).
Berbagai algoritma telah dikembangkan dan digunakan untuk inversi data
gelombang permukaan eksperimental. Prosedur inversi yang mungkin dapat
dibagi menjadi dua kategori; pencarian lokal dan prosedur pencarian global.
Pembagian ini didasarkan pada bagaimana parameter model diperbarui antara
iterasi selama pencarian parameter set yang paling mungkin. (Xia dkk., 1999).
24
BAB III
METODE PENELITIAN
3.1 Lokasi penelitian
Lokasi penelitian dengan metode multichannel analysis of surface wave
dilakukan di daerah Surabaya, Jawa Timur. Penelitian dilakukan di 45 titik yang
tersebar di daerah Surabaya ditunjukkan oleh gambar 3.1. Titik penelitian
tersebar berdasarkan hasil gridding dan pertimbangan geologinmya.
25
Gambar 3.1 Titik lokasi pengukuran di Surabaya
26
3.3 Desain Akuisisi
Berikut ini merupakan ilustrasi dari desain pengukuran multichannel
analysis of surface wave (MASW) di daerah surabaya :
Permukaan
tanah
= Shot
= Geophone
Gambar 3.2 Desain akuisisi MASW
27
Mulai
Raw Data
Pengolahan Data
Selesai
28
Gambar 3.4 Diagram alir pengambilan data
29
Mulai
Rekaman data
seismik
Filtering
FFT
Gambar Spektral
Dispersi
Pembentukan kurva dispersi Eksperimental
Proses inversi
Penentuan Model
Profil Vs
Lapisan Tanah
Pembentukan kurva Dispersi Teori
Proses iterasi
Selesai
30
Mulai
Data Vs dan G
Membuat Batas
Membuat Surface
Membuat Grid
Membuat Lapisan
Petrophysical Modelling
Model 3D
Selesai
31
3.5 Pengolahan Data
3.5.1 Filtering Data Menggunakan Seisee
Data seismik pada gambar 4.1 difilter frekuensi menggunakan bandpass
filter pada software pengolahan seismik dengan f1=4hz, f2=8hz,f3=40 hz dan
f4=60 hz.
Gambar 3.8 merupakan hasil trace gelombang line 25 pada near offset
setelah proses filtering. Selanjutnya data disimpan dalah format .sgy dan siap
dimasukkan ke software pengolahan MASW.
3.5.2 Pengolahan Data MASW
Setelah data difilter kemudian data dimasukkan ke software pengolahan
MASW dan dilakukan pengolahan untuk melihat profil Vs di Surabaya.
32
Gambar 3.9 Gambar data padaPengolahan MASW
Gambar 3.10 merupakan hasil picking surface wave line 25 shoot pertama
(near offset). Picking gelombang dapat diambil dari first break sampai
gelombang permukaannya, namun untuk hasil gambar spektral dispersi yang
lebih jelas cukup ambil gelombang permukaanya saja. Setelah dilakukan picking
surface wave kemudian masuk ke menu spectral analysis untuk mengatur
parameter X (frekuensi) dan Y (phase velocity) gambar spektral dispersi yang
ingin diperlihatkan pada sub menu analysis setting. Selanjutnya memasukkan
parameter berupa ketebalan lapisan dan Vs model pada sub menu synthetic
model. Tahap terakhir picking kurva dispersi pada fundamental mode (M0).
33
Gambar 3.11 Picking fundamental mode
(a)
34
(b)
Gambar 3.7 Hasil Inversi (a) Profil Vs 1D (b) Stratigrafi dan nilai error
35
HALAMAN INI SENGAJA DIKOSONGKAN
36
BAB IV
ANALISA DATA DAN PEMBAHASAN
4.1 Analisa Data
Tabel berikut ini merupakan hasil dari pengolahan data lapangan yang
menunjukkan nilai kecepatan gelombang geser (Vs) pada kedalaman 0-30 m,
nilai Vs per lapisan dan modulus geser (G0).
37
680678 9195845 11.22 195.04 68.47
680678 9195845 14.86 273.29 134.44
680678 9195845 18.69 393.77 279.1
680678 9195845 22.74 466.67 392
680678 9195845 37.17 531.54 508.57
686809 9192410 7.19 107.85 20.94
686809 9192410 12.44 326.03 191.33
6 686809 9192410 258.58 24.16 436.23 342.54 Qtp
686809 9192410 35.48 568.55 581.84
686809 9192410 46.94 676.59 823.99
680262 9190339 2.35 106.67 6.4
680262 9190339 6.02 159.98 10.26
7 680262 9190339 258.016 12.28 266.66 60.12 Qtp
680262 9190339 17.28 333.32 195.46
680262 9190339 34.44 533.26 456.43
682912 9194836 3.01 65.07 7.62
682912 9194836 5.51 160.00 46.08
8 682912 9194836 239.296 9.02 266.45 127.79 Qtp
682912 9194836 14.91 332.00 198.4
682912 9194836 43.59 399.98 287.97
683911 9189615 5.3 130.59 30.7
683911 9189615 8.12 134 32.31
9 683911 9189615 235.528 17 328.13 193.81 Qtp
683911 9189615 28.3 334.36 201.24
683911 9189615 44.33 525.48 497.03
688312 9192460 3 100 18
688312 9192460 15 200 72
10 688312 9192460 235.022 27 333.32 199.99 Qtp
688312 9192460 39 466.67 392
688312 9192460 42 600 648
11 677714 9199468 251.466 5.39 97.94 17.27 Qpk
38
677714 9199468 7.98 261.61 123.19
677714 9199468 11.48 333.26 199.91
677714 9199468 21.89 465.21 389.55
677714 9199468 37.15 533.25 511.85
685569 9196131 5.27 133.33 32
685569 9196131 14.08 266.67 128
12 242.062 Qpk
685569 9196131 28.08 333.33 200
685569 9196131 34.19 466.67 392
686889 9195469 4.28 102.77 19.01
686889 9195469 7.15 196.12 69.24
13 686889 9195469 238.606 12.07 249.11 111.7 Qpk
686889 9195469 27.15 363.79 238.22
686889 9195469 33.94 434.94 340.51
698681 9192413 2.66 92.81 15.51
698681 9192413 12.66 188.94 64.25
14 698681 9192413 218.714 18.67 260.46 122.11 Qa
698681 9192413 23.67 284.11 145.29
698681 9192413 31.17 308.27 171.05
683762 9197724 1.50 93.33 15.68
683762 9197724 9.01 137.57 34.06
15 683762 9197724 218.712 12.51 264.16 125.61 Qa
683762 9197724 17.51 331.16 197.41
683762 9197724 32.65 397.47 284.36
693880 9197348 3.7 104.57 19.68
693880 9197348 10.63 203.37 74.44
16 693880 9197348 215.472 26.28 272.95 134.11 Qa
693880 9197348 34.52 305.2 167.67
693880 9197348 38.08 417.7 314.05
682420 9191068 4.74 66.67 8
17 682420 9191068 215.03 14.51 333.33 200 Qa
682420 9191068 26.30 334.43 201.31
39
682420 9191068 33.24 429.39 331.87
682420 9191068 47.86 533.33 512
695794 9194355 6.29 131.18 30.97
695794 9194355 16.17 221.38 88.21
18 695794 9194355 214.872 35.8 302.69 164.92 Qa
695794 9194355 41.44 421.51 319.81
695794 9194355 45.56 440.2 348.8
700329 9193322 3.40 82.05 12.12
700329 9193322 5.9 160.00 46.08
19 700329 9193322 214.866 19.43 262.44 123.98 Qa
700329 9193322 39.17 329.32 195.21
700329 9193322 48.7 356.32 228.54
696495 9201873 3.86 107.46 20.79
696495 9201873 7.95 154.44 42.93
20 696495 9201873 214.388 27.31 278.51 139.62 Qa
696495 9201873 37.31 364.55 239.21
696495 9201873 47.31 372.39 249.62
696220 9190378 4 97.82 17.22
696220 9190378 6.5 119.7 25.79
21 696220 9190378 212.454 11.5 206.55 76.79 Qa
696220 9190378 27.5 330.76 196.92
696220 9190378 34.81 328.95 194.77
690628 9189550 9.97 117.6 24.89
690628 9189550 17.31 245.34 108.34
22 690628 9189550 208.33 34.27 318.63 182.75 Qa
690628 9189550 40.45 358.48 231.31
690628 9189550 46.99 411.77 305.2
699661 9196172 3.14 113.43 23.16
699661 9196172 9.14 142.68 36.64
23 208.13 Qa
699661 9196172 14.8 228.29 93.81
699661 9196172 22.48 276.43 137.55
40
699661 9196172 31.04 300.66 162.72
689948 9187996 2.66 92.81 15.51
689948 9187996 12.66 188.94 64.25
24 689948 9187996 205.042 18.67 260.46 122.11 Qa
689948 9187996 23.67 284.11 145.29
689948 9187996 31.17 308.27 171.05
696925 9189745 2.5 72.09 9.35
696925 9189745 10.81 167.17 50.3
25 696925 9189745 202.11 24.12 262.97 124.48 Qa
696925 9189745 27.53 272.97 134.12
696925 9189745 37.41 272.73 133.89
697132 9187419 6.87 111.22 22.27
697132 9187419 8.38 150.93 41
697132 9187419 16.04 237.15 101.23
26 201.326 Qa
697132 9187419 19.42 265.67 127.05
697132 9187419 26.35 332.56 199.08
697132 9187419 39.31 395.25 281.2
696949 9198850 2.5 83.95 12.69
696949 9198850 12.45 193.74 67.57
27 696949 9198850 201.322 26.78 270.62 131.82 Qa
696949 9198850 33.41 361.58 235.33
696949 9198850 42.96 466.16 391.15
687380 9197516 1.00 93.33 15.68
687380 9197516 6.42 114.5 31.95
28 687380 9197516 200.485 16.42 198.78 32.08 Qa
687380 9197516 23.6 333.33 200
687380 9197516 30.33 411.34 499.52
684735 9190587 6.93 101.87 18.68
684735 9190587 16.09 260.8 122.43
29 199.15 Qa
684735 9190587 26.92 344.17 213.22
684735 9190587 34.89 401.95 290.81
41
684735 9190587 44.4 499.94 449.89
696184 9195235 4.01 93.54 15.75
696184 9195235 11.95 133.52 32.09
696184 9195235 14.49 208.06 77.92
30 198.034 Qa
696184 9195235 16.49 278.94 140.06
696184 9195235 19.19 339.57 207.56
696184 9195235 27.42 373.36 250.91
691226 9193916 8.71 98.79 17.57
691226 9193916 17.38 175.07 55.17
691226 9193916 22.09 226.28 92.16
31 196.384 Qa
691226 9193916 30.1 248.49 111.14
691226 9193916 37.41 274.26 135.4
691226 9193916 45.07 326.76 192.19
694124 9189567 7.02 104.37 19.61
694124 9189567 9.58 141.92 36.25
694124 9189567 16.93 218.26 85.75
32 195.396 Qa
694124 9189567 22.48 261.04 122.65
694124 9189567 26.23 328.57 194.33
694124 9189567 32.45 383.02 264.07
692455 9190558 1.79 61.33 6.77
692455 9190558 5.16 102.26 18.82
33 692455 9190558 191.92 8.68 214.72 82.99 Qa
692455 9190558 27.57 299.63 161.6
692455 9190558 57.36 307.15 169.81
698909 9189042 4.06 62.38 7
698909 9189042 11.61 201.88 73.36
34 698909 9189042 190.488 31.15 327.31 192.83 Qa
698909 9189042 40.97 339.06 206.93
698909 9189042 43.81 402.65 291.83
699809 9190739 5.21 99.84 17.94
35 189.816 Qa
699809 9190739 17.21 147.56 39.19
42
699809 9190739 27.44 392.32 277.05
699809 9190739 39.41 466.19 391.19
699809 9190739 48.37 532.87 511.11
698505 9194835 5.65 115.74 24.11
698505 9194835 13.86 198.04 70.6
36 698505 9194835 189.484 23.88 219.75 86.92 Qa
698505 9194835 29.05 251.6 113.94
698505 9194835 34.73 410.8 303.77
687360 9199593 2.46 70.46 8.94
687360 9199593 6.97 101.50 18.54
37 687360 9199593 188.118 11.65 203.63 74.64 Qa
687360 9199593 16.67 251.73 114.07
687360 9199593 24.18 285.90 147.13
678657 9201357 4.86 93.3 14.94
678657 9201357 14.86 91.19 14.97
38 678657 9201357 186.834 27.98 333.3 199.86 Qa
678657 9201357 32.74 366.4 241.64
678657 9201357 49.83 439.97 348.44
680608 9199148 2.66 87.63 13.82
680608 9199148 12.55 171.46 52.92
39 185.522 Qa
680608 9199148 17.85 216.38 84.28
680608 9199148 32.17 317.5 181.45
693893 9199517 6.71 109.47 21.57
693893 9199517 9.26 159.69 45.9
693893 9199517 18.34 233.23 97.92
40 183.606 Qa
693893 9199517 24.95 292.96 154.49
693893 9199517 29.96 346.41 216.01
693893 9199517 43.13 399.26 286.93
681343 9198686 7.75 119.48 25.7
41 681343 9198686 180.426 13.19 119.46 25.69 Qa
681343 9198686 15.71 208.89 78.54
43
681343 9198686 17.9 249.31 111.88
681343 9198686 20.04 314.58 178.13
681343 9198686 37.48 358.04 230.75
686850 9200859 2.72 82.1 12.14
686850 9200859 12.15 176.5 56.07
42 179.574 Qa
686850 9200859 24.48 185.2 61.74
686850 9200859 39.82 294.26 155.86
682238 9202931 1.63 83.95 12.69
682238 9202931 10.22 199.84 71.89
43 682238 9202931 176.492 24.09 215.885 83.86 Qa
682238 9202931 30.14 215.5 83.93
682238 9202931 49.69 215.94 83.59
695412 9193225 5.02 149.91 40.45
695412 9193225 9.15 205.67 40.53
44 695412 9193225 284.25 14.16 354 288 Tpl
695412 9193225 22.12 466.67 392
695412 9193225 32.51 527.78 501.39
693820 9188848 2.95 93.33 15.68
693820 9188848 8.68 163.36 48.04
693820 9188848 18.02 293.68 155.25
45 243.04 Qpk
693820 9188848 25.44 453.46 370.13
693820 9188848 39.04 531.83 509.12
693820 9188848 50.6 599.94 647.87
Berdasarkan hasil dari tabel 4.1 menujukkan adanya perbedaan nilai pada
4 formasi di daerah Surabaya yaitu Formasi Lidah (Tpl), Formasi Pucangan
(Qtp), Formasi Kabuh (Qpk), dan Formasi Alluvium (Qa). Setiap formasi
geologi memiliki kisaran nilai Vs30 yang berbeda. Formasi lidah merupakan
formasi tertua di Surabaya, dimana formasi ini memliki kisaran nilai Vs30 yang
cukup tinggi dibandingkan formasi lainnya yaitu 252 m/s – 294 m/s. Formasi Qtp
(Pucangan) memiliki rentan nilai Vs30 antara 235-270. Formasi Qpk (Kabuh)
memiliki nilai Vs30 antara 238-243 m/s dan Formasi Aluvium memiliki rentan
Vs30 antara 176-218 m/s. Formasi aluvium memiliki nilai Vs30 yang rendah
karena komposisi penyusunnya merupakan kerikil, pasir, lempung juga fosil
44
kerang. Komposisi penyusunnya tersebut memiliki nilai Vs < 175 m/s dan pasir
kerikil memiliki Vs antara 175-350 m/s. Sementara untuk formasi Qpk dan Qtp
yang tersusun atas batu pasir dan formasi Tpl yang dominan batu lempung
sehingga nilai Vs30 nya cukup tinggi.
Hasil nilai Vs30 tersebut menunjukkan bahwa formasi yang paling muda di
Surabaya adalah formasi aluvium dengan nilai Vs paling rendah dibanding
formasi lain karena sedangkan, formasi lidah memiliki nilai Vs cukup tinggi
menggambarkan lapisannya cukup keras dan menunjukkan formasi paling tua.
Hal ini selaras dengan pernyataan dari sukardi tahun 1992 bahwa batuan dasar
untuk surabaya merupakan formasi Lidah yang berumur Pliosen (pre-tertiary).
Formasi yang paling tua seharusnya berada di paling bawah lapisan sehingga
akan lebih kompak dari lapisan diatasnya namun karena adanya gaya arah
vertikal, membuat lapisan tanah bergerak keatas dan membentuk antiklin. Maka
lapisan paling tua akan berada di tengah dan lapisan paling muda akan berada di
samping sehingga saat mengalami erosi lapisan paling tua akan berada di puncak
antiklin
Tabel 4.1 juga menunjukkan bahwa nilai kecepatan gelombang geser (Vs)
perlapisannya berkisar antara 93 m/s - 676 m/s. Nilai Vs lapisan mencapai 300
m/s pada kedalaman >30 m pada mayoritas titik pengukuran yang berupa
formasi alluvium. Hal ini menunjukkan bahwa mayoritas penyusun batuan di
surabaya merupakan pasir, lempung dan kerikil. Sementara untuk formasi Qpk
dan Qtp yang tersusun atas batu pasir dan formasi Tpl yang dominan batu
lempung, formasi ini memiliki nilai Vs diatas 350 m/s - 600 m/s. Sedangkan nilai
modulus gesesr berkisar antara 7 MPa – 823 MPa dengan kenaikan yang
berbanding lurus dengan nilai kecepatan gelombang geser (Vs) dan kedalaman.
4.2 Hasil
Data Vs30 dari tabel 4.1 dipetakan untuk mengetahui distribusi nilai Vs30
di kota Surabaya secara lokal. Hasil distribusi kecepatan gelombang geser pada
kedalaman 0-30 m (Vs30) di 45 lokasi di Kota Surabaya ditunjukkan oleh
gambar 4.1. Terdapat 5 zona berdasarkan klasifikasi nilai kecepatan gelombang
geser (Vs30). Zona pertama yaitu daerah dengan nilai Vs30 sebesar 170-200 m/s,
yang kedua 200-220 m/s, yang ketiga 230-247 m/s, yang keempat sebesar 247-
270 m/s, dan yang terakhir merupakan zona dengan nilai Vs30 sebesar 270-294
m/s.
Hasil pemetaan pada gambar 4.1 menunjukkan bahwa mayoritas daerah
surabaya berada pada zona kedua yaitu nilai Vs30 sebesar 200-220 m/s. Hal ini
menujukkan bahwa daerah Surabaya merupakan daerah bertanah lunak karena
mayoritas penyusunya merupakan formasi alluvium sehingga lapisannya belum
terkompaksi. Kearah Utara menujukkan bahwa nilai Vs semakin rendah. Hal ini
mengindikasikan adanya sedimentasi yang tebal pada daerah Utara dikarenakan
daerah pesisir dan menunjukkan bahwa lapisan pada daerah ini belum
45
terkompaksi. Sedangkan pada daerah Selatan sebelah Barat dan daerah Barat
nilai Vs tinggi dibandingkan daerah sekitar yang dikarenakan adanya dua
antiklin yaitu antiklin lidah dan guyangan. Pada daerah Surabaya bagian Timur
terdapat daerah dengan nilai Vs yang cukup tinggi. Hasil pemetaan pada gambar
4.1 menunjukkan bahwa pada daerah antiklin lidah terdapat daerah yang diapit
nilai tinggi, hal ini disebabkan oleh adanya patahan (Kotak warna hitam)
sehingga daerah tersebut termasuk zona lemah dibandingkan daerah pinggirnya.
Berdasarkan nilai Vs30, klasifikasi tanah di surabaya masuk dalam kategori
kelas C dan D beradsarkan Eurocode8. Kategori kelas C yang berada pada daerah
Antiklin merupakan endapan sand padat atau setengah padat yang tebal, gravel
atau clay padat dengan beberapa puluhan hingga ratusan meter. Endapan tanah
berbutir kasar kasar sedang atau tanah butiran halus menengah sedang.
Sedangkan kategori kelas D dengan niali kurang <180 m/s merupakan Endapan
tanah berbutir kasar dengan konsentris jelek atau tanah berbutir halus yang
mengandung konsentris yang buruk.
46
Gambar 4.1 Peta persebaran nilai Vs30 di Surabaya
47
Berikut merupakan pemodelan 3D hasil inversi kecepatan gelombang geser
(Vs) per lapisan di daerah Surabaya :
49
Pada sisi Utara ke arah Timur penampang Vs tidak terlalu menunjukkan
perbedaan, ditunjukkan oleh gambar 4.5. Nilai Vs pada kedalaman 30 m berkisar
antara 100-200 m/s.
50
Gambar 4.7 Penampang 3D G0 dari arah Timur ke Selatan
1’
2’
5’
5
2
3’ 2’ 2
4’
4 2
2
2
3 2 1
51
Gambar 4.8 merupakan gambar intersection antar titik pengukuran. Garis
biru menjukkan Xline yaitu perpotongan dari arah Timur ke Barat. Garis merah
menujukkan perpotongan dari arah Utara ke Selatan. Berikut merupakan
penampang 2D berdasarkan intersection antar ttitik pengukuran pada model Vs:
52
Inline 20, namun lebih keras karena pada kedalaman 30 m sudah mulai lapisan
keras dengan niali Vs >360 m/s.
53
lapisan lempung lunak, lapisan kedua dan ketiga merupakan endapan pasir atau
clay padat, sedangkan lapisan keempat berupa endapan pasir atau lempung yang
sangat padat atau sandstone. Pada penampang bagian sisi sebelah Barat terjadi
kenaikan struktur lapisan tanah, lapisan pertama berupa lapisan kedua pada
penampang sebelah Timur (kanan) yaitu nilai Vs sebesar 150-200 m/s dengan
kedalaman 0-12 m, selanjutnya lapisan kedua 12-26 m dengan nilai Vs 250-300
m/s , dan terakhir lapisan ketiga 26-50 m dengan Vs 350-550 m/s. Hal ini
dikarenakan adanya antiklin pada daerah penampang 4-4’ sehingga adanya
kenaikan struktur tanah. Sehingga dapat disimpulkan bahwa semakin ke arah
Barat (antiklin) lapisan tanah semakin keras ditunjukkan oleh nilai Vs yang
tinggi pada kedalaman 26 m.
54
Gambar 4.14 Crossection pada Xline 60 dan Inline 60 pada model Vs
55
Gambar 4.16 Inline 26 (Intersection 2-2’) pada model G0
56
Gambar 4.18 Xline 70 (Intersection 4-4’) pada model G0
57
100 MPa, lapisan kedua 100-200 MPa, dan lapisan terakhir 200-400 MPa.
Terlihat pada gambar 4.15 dan 4.19 yang merupakan perpotongan di daerah
antiklin bahwa lapisan terakhir mulai terlihat pada kedalaman 28 m. Namun pada
daerah aluvium seperti Xline 30, inline 20 dan inline 26 range lapisan terakhir
tidak mencapai nilai 200.
4.3 Pembahasan
Pengukuran data metode Multichannel Analysis Surface Wave (MASW)
pada tahun 2018 dilakukan sebanyak 45 titik di daerah Surabaya, Jawa Timur.
Pengukuran MASW bertujuan sebagai evaluasi kondisi lokal atau tapak lokal
Surabaya. Pengukuran dilakukan dengan konfigurasi pengukuran yaitu jarak
offset 5 m, spasi antar geophone 4 m, waktu perekaman 2 s dan sampling rate
0.001 s. Hasil akhir dari penelitian berupa pemetaan kecepatan gelombang geser
pada kedalaman 0-30 m (Vs30), serta pemodelan kecepatan gelombang geser
dan modulus geser di daerah Surabaya berdasarkan data lapangan.
Data mentah yang telah diperoleh dari lapangan diolah menggunakan
software pengolahan MASW untuk memperoleh nilai Vs,Vs30 dan G. Data
mentah yang akan diolah difilter terlebih dahulu menggunakan filter bandpass
pada software seismik dengan pengaturan f1= 4 Hz; f2 = 8Hz; f3 = 40 Hz; f4 =
60 Hz. Fiter bandpass digunakan untuk memilih frekuensi yang diinginkan dan
menghilangkan noise frekuensi tinggi yang berasal dari noise kabel dan
memperjelas gelombang Rayleigh sehingga memudahkan proses picking di
58
software pengolahan MASW. Proses pertama pengolahan MASW yaitu
mengambil gelombang permukaan dan mengatur geometri sesuai dengan
pengukuran dilapangan. Selanjutnya dilakukan proses pencitraan gambar
spektral dispersi dan dilakukan input inisial model dengan parameter input
berupa ketebalan dan Vs tiap kedalaman. Untuk input parameter ketebalan
formasi Qa dan 3 formasi yang ada di antiklin (QTp, Tpl, dan Qpk) diseragamkan
dengan lapisan pertama 3-5 meter, lapisan kedua 5-7 meter, lapisan ketiga 5-7
meter, lapisan keempat 7-9 meter, lapisan kelima 7-10 meter dan lapisan keenam
10 meter. Sementara untuk masukkan Vs pada formasi Qa dan formasi di antiklin
parameter modelnya berbeda. Untuk lapisan Qa Vs lapisan pertama adalah 70-
100 m/s, lapisan kedua 150 m/s, lapisan ketiga 200 m/s, lapisan keempat 250
m/s, lapisan kelima 300 m/s dan lapisan keenam 350 m/s. Sementara untuk
lapisan Qtp, Qpk dan Tpl Vs lapisan pertama 150-200 m/s, lapisan kedua 250
m/s, lapisan ketiga 300 m/s lapisan keempat 350 m/s, lapisan kelima 400 m/s dan
lapisan keenam 450 m/s. Perbedaan input Vs untuk model awal disesuaikan
dengan karakter formasi. Namun perbedaan sebenarnya hanya dilapisan pertama
dan kedua saja yang perbedaanya juga tidak terlalu signifikan karena lapisan Qa
dan formasi penyusun antiklin sama-sama terdiri dari batuan sedimen. Untuk
input data densitas digunakan 1800 Kg/m3 sesuai software. Hal ini dikarenakan
mayoritas tanah Surabaya berupa clay, sand, dan gravel. Densitas clay 1200-
2200 Kg/m3 dan densitas sand, gravel sebesar 1400-2300 Kg/m3 sehingga nilai
1800 Kg/m3 berada di pertengahan. Selajutnya picking pada fundamental mode
yaitu amplitudo atau energi terbesar dari gelombang rayleigh. Tahap terakhir
melakukan inversi untuk menndapatkan nilai Vs30, Vs dan nilai G.
Untuk mengetahui persebaran nilai Vs30, maka dilakukan pemetaan
menggunakan software arcgis. Distribusi hasil pemetaan kecepatan gelombang
geser pada kedalaman 0-30 m (Vs30) di 45 lokasi pada gambar 4.1 menunjukkan
5 zona berdasarkan klasifikasi kecepatan gelombang geser (Vs30), yaitu pada
zona 170-200 m/s ; 200-220 m/s ; 230-247 m/s; 270-294 m/s. Klasifikasi ini
sangat bersesuain dengan studi sebelumnya mengenai kondisi geologi Surabaya
yang diteliti oleh Mufida dkk pada tahun 2012. Mayoritas lokasi di surabaya
berada pada zona dengan Vs30 200-220 m/s. Berdasarkan nilai Vs30, klasifikasi
tanah di surabaya masuk dalam kategori kelas C dan D beradsarkan Eurocode8.
Kelas C merupakan Endapan sand padat atau setengah padat yang tebal, gravel
atau clay padat dengan beberapa puluhan hingga ratusan meter. dan kelas D
merupakan endapan tanah berbutir kasar dengan konsentris jelek atau tanah
berbutir halus yang mengandung konsentris yang buruk.
Pemodelan kecepatan gelombang geser (Vs) dan modulus geser (G0)
dilakukan dengan software pemodelan. Data Vs dan G0 disimpan dalam notepad
beserta keterangan nilai kedalaman dan koordinat menggunakan sistem UTM.
Selanjutnya input data kedalam software dan membuat batas titik pengukuran,
batas yang dibuat akan digunakan untuk pembuatan permukaan model. Proses
59
berikutnya adalah pembuatan permukaan atas dan bawah sebagai batas untuk
pembuatan grid sederhana. Setelah gridding, dilakukan pembuatan lapisan-
lapisan untuk model lapisan bumi, kemudian dilakukan well log up scalling pada
nilai yang ingin dijadikan model. Tujuan dari well log up scalling adalah agar
masukkan yang berupa titik-titik dapat dibaca sebagai data sumur. Tahap terakhir
adalah melakukan petrophysicall modelling.
Hasil dari pemodelan 3D Vs pada daerah Utara, daerah Timur, dan sebagian
daerah Selatan merupakan daerah lempung lunak yang cukup tebal. Hal ini sesuai
dengan hasil pemetaan distribusi nilai Vs30 pada gambar 4.1 bahwa daerah
tersebut bertanah lunak. Sebaran nilai Vs pemodelan 3D menujukkan bahwa sisi
bagian Barat memiliki nilai Vs yang tinggi dariapada sisi bagian Utara. Pada
kedalaman >20m nilai Vs mencapai 400-500 m/s dan mengalami penurunan nilai
Vs yang cukup signifikan kearah Utara, dikarenakan daerah uatara merupakan
formasi alluvium. Terdapat naikan yang ditunjukkan pada gambar 4.2 yang
diindikasikan sebagai antiklin lidah. Pada bagian Timur nilai Vs cukup rendah
antara 100-250 m/s sampai kedalaman 40 m. Pada sisi Selatan dari arah Timur
nilai Vs rendah berkisar antara 100-250 m/s sampai kedalaman 40 m, namun
semakin ke Barat nilai Vs semakin tinggi yaitu 300 m/s pada kedalaman 30 m.
Dilihat dari model 3D, terdapat 4 lapisan di daerah Surabaya, namun mayoritas
hanya terdiri dari 3 lapisan, karena lapisan pertama dengan nilai Vs <100 m/s
hanya ada di sisi Utara dan Timur dan lapisan ini sangat tipis pada kedalaman
maksimal 5 m. Hasil dari pemodelan modulus geser yang ditunjukkan gambar
4.6 dan gambar 4.7 memperlihatkan bahwa daerah Timur ke Selatan di sisi Timur
pada kedalaman 0-30 m nilai G0 sebesar <100-200 Mpa dan kearah Selatan pada
kedalaman 20-30 m nilainya sebesar 250 Mpa. Sedangkan pada daerah Utara
sama seperti daerah Timur yaitu <100-200 Mpa dan menipis kearah Barat pada
kedalaman 20-30 m sebesar 350 Mpa. Berdasarkan kedua gambar tersebut dapat
disimpulkan bahwa tren nilai pada modulus geser sama dengan nilai Vs.
Berdasarkan data hasil pengukuran parameter dinamis tanah yang berupa
kecepatan gelombang geser (Vs) dan modulus geser (G0) menggunakan metode
multichannel analysis of surface wave di Surabaya menunjukkan bahwa kondisi
lokal atau tapak lokal merpakan daerah bertanah lunak. Sehingga apabila terjadi
gempa bumi maka, daerah ini memungkin berdampak cukup besar terutama
daerah Utara dan Timur yang memiliki nilai Vs rendah. Hal ini selaras dengan
(Rasimeng, 2016) yang menyatakan jenis tanah yang lunak memiliki Vs yang
rendah, dimana jenis tanah ini biasanya adalah sedimen dengan umur masih
muda dan belum terkompaksi Jika gelombang seismik melewati jenis tanah ini,
gelombang tersebut akan mengalami amplifikasi yang cukup besar. Nilai
modulus geser semakin tinggi pada peningkatan kedalaman. Berdasarkan hasil
pemodelan menunjukkan bahwa daerah antiklin lebih padat daripada daerah
Utara dan Timur sehingga daerah ini tidak mudah untuk mengalami pergeseran
seperti yang dijelaskan oleh Braja pada buku mekanika tanah. Kegunaan kategori
60
tanah berdasarkan kecepatan gelombang geser dapat diilustrasikan oleh
kesamaan dalam pola distribusi kerusakan dan kecepatan. Namun, beberapa efek
lain (kualitas bangunan, resonansi frekuensi) juga harus dipertimbangkan. Selain
itu, kerusakan distribusi tidak dapat diprediksi sepenuhnya dengan hanya
menggunakan Vs dan G0. Peta kondisi tanah terperinci berdasarkan distribusi
kecepatan gelombang geser (Vs30) berada dalam EC8 standar dan parameter
geoteknik lainnya sangat penting untuk tidak hanya penyelesaian daerah
perkotaan tetapi juga untuk situs seleksi dan studi evaluasi keamanan situs.
61
HALAMAN INI SENGAJA DIKOSONGKAN
62
BAB V
PENUTUP
5.1 Kesimpulan
1. Nilai Vs30 Surabaya berkisar antara 176-294 m/s. Nilai kecepatan
gelombang geser (Vs) perlapisannya berkisar antara 93 m/s - 676 m/s. Nilai
modulus geser (G0) perlapisan berkisar antara 7 MPa – 823 MPa.
2. Hasil pemetaan Vs30 menunjukkan bahwa klasifikasi tanah di surabaya
masuk dalam kategori kelas C dan D berdasarkan Eurocode8. Sedangkan
sebaran nilai Vs dan G0 yang rendah yaitu 100-200 m/s dan <100-200 Mpa
pada hasil pemodelan mengalami penebalan dari arah Barat menuju ke
Utara dan dari arah Timur ke Selatan menujukkan nilai Vs dan G0 yang
relatif sama yaitu sekitar 100-200 m/s dan <100-200 Mpa pada kedalaman
0-30 m.
3. Nilai Vs menujukkan jenis tanah, semakin rendah nilai Vs maka tanah
semakin lunak, sedangkan nilai G0 ini menunjukkan bahwa semakin tinggi
nilainya maka lebih tidak mudah untuk mengalami pergeseran.
Berdasarkan hasil pemodelan ditunjukkan bahwa daerah antiklin memiliki
nilai Vs dan G0 yang lebih tinggi daripada daerah alluvium khususnya
daerah Utara, sehingga daerah antiklin lebih tidak berisiko daripada daerah
lainnya.
5.2 Saran
1. Menambahkan parameter dinamis lainnya untuk mendukung hasil penelitian
2. Menggunakan data SPT Surabaya sebagai data pendukung model sintetis.
Hal ini untuk menambahkan informasi densitas sesuai dengan litologi
batuan pada saat pembuatan model sintetis.
3. Pengukuran dapat dilakukan saat musim kemarau, karena apabila dilakukan
saat musim hujan tanah akan cenderung becek dan akan mempengaruhi nilai
Vs yang diukur.
63
HALAMAN INI SENGAJA DIKOSONGKAN
64
DAFTAR PUSTAKA
Aisyah, V., Suharna & S., A. (2011), “Groud Amplification Mapping Using
HVSR Method Horizontal to Vertical Spectral Ratio In Patuk,
Gunungkidul, Yogyakarta (Indonesia)”. Proceedings JCM Makassar, The
36th HAGI and 40th IAGI Annual Convention and Exhibition.
Ariestianty, Susy K. Raihan Taha, Mohd. Nayan, Khairul Anuar Mohd. Chik,
Zamri. (2009), “Penentuan Modulus Geser Tanah Menggunakan Metode
Analisis Multi-Channel Gelombang Permukaan”, Jurnal Ilmiah Semesta
Teknika Vol. 12, No. 2, 185-198, November 2009.
Badan Pusat Statistik. (2015), “Jawa Timur Dalam Angka. Badan Pusat
Statistik”, BPS Provinsi Jawa Timur.
Bahri, Syaeful., Madlazim. (2012), “Pemetaan Topografi, Geofisika, dan
Geologi Kota Surabaya”, Jurnal Penelitian Fisika dan Aplikasinya Vol 2
No 2, hal. 23-28.
Bessason, B. & Erlingsson, S. (2011), “Shear wave velocity in surface
sediments”, Jökull 61, 51-64.
Burger, Robert. (1992), “Exploration Geophysics of the Shallow Surface”,
Pentice-Hall, New Jersey
Carvalho, E. C. (2008), “Eurocode 8 General Rules and Seismic Actions”,
Dissemination of Information Workshop : Brussels.
Cen. (2004), “Eurocode 8—design of structures for earthquake resistance. Part
1: general rules, seismic actions and rules for buildings”, European standard
EN 1998-1, December 2004, European Committee for Standardization :
Brussels.
Cho, G. C. & Santamarina, J. C. (2001), “Unsaturated particulate materials:
particle-level studies”, J. Geotech. Geoenviron. Engng 127, No. 1, 84–96
Dal Moro, G. (2014), “Surface Wave Analysis for Near Surface Applications”,
Elsevier Inc.
Erik H. Saenger, Stefan M. Schmalholz, Marc-A. Lambert, dkk. (2009), “A
passive seismik survey over a gas field: Analysis of Low-Frequency
Anomalies”, Geophysic, No 74(2), hal. 29-40.
Evrett, M. E. (2013), “Near-Surface Applied Geophysics 1st ed”, Cambridge
University Press, New York.
Foti, S., Lai, C.G., Rix, G.J., Strobbia, C. (2015), “Surface Wave Methods for
Near-Surface Site Characterization”, CRC Press, Boca Raton.
65
Grandis, H. (2009), “Pengantar Pemodelan Inversi Geofisika”. Himpunan Ahli
Geofisika Indonesia (HAGI), Jakarta.
Kanlı, A. I., Tildy, P., Pr´onay, Z., Pınar, A., and Hermann, L. (2006), “Vs30
mapping and soil classification for seismic site effect evaluation in Dinar
region”, Geophys. J. Int., 165, 223-235, SW Turkey.
Lin, C.P., Chang, C.C. And Chang, T.S. (2004), “The use of Masw Method in
the Assessmentof Soil Liquefaction Potential”, Department of Civil
Engineering National Chiao Tung University, Taiwan.
Lowrie, W. (2007, “Fundamentals of Geophysics 2nd edition”, Cambridge
University Press, 381p. : London
Luna, R. & Jadi, H. (2000), “Determination of Dynamic Soil Properties Using
Geophysical Methods”, Preceedings of the First International Conference
on the Application of Geophysical and NDT Methodologies to
Transportation Facilities and Infrastructure, St. Louis, December 2000.
McMechan, G. & Yedlin, M. J. (1981), “Analysis of dispersive waves by wave
field transformation”, Geophysics, 46 (6), 869-874.
Mufida, A., Santosa, B.J., Warnana, D.D. (2013), “Profiling Kecepatan
Gelombang Geser (Vs) Surabaya Berdasarkan Pengolahan Data
Mikrotremor”,. Jurnal Sains Dan Seni Pomits Vol. 2, no.2., Surabaya.
Nakamura, Y., Sato, T., & Nishinaga, M. (2000), "Local Site Effect of Kobe
Based on Microtremor", 3–8.
Rahmaningtyas, A. P., Purwanto, M. S., & Widodo, A. (2017), "Identifikasi
Percepatan Tanah Maksimum (PGA) dan Kerentanan Tanah
Menggunakan Metode Mikrotremor Jalur Sesar Kendeng", Jurnal
Geosaintek, 3(2), 107.
Sunardi, B., Arifin, J., Susilanto, P., Ngadmanto, D., Nurdiyanto, B., & Shahzad,
S. (2012), "Kajian Potensi Bahaya Gempabumi Daerah Sumbawa
Berdasarkan Efek Tapak LokaL", Jurnal Meteorologi Dan Geofisika, 13,
131–137.
Ólafsdóttir, Elín Ásta. (2014), “Multichannel Analysis of Surface Waves
Methods for dispersion analysis of surface wave data”, University Of
Iceland.
Park, C. B., Miller, R. D. & Xia, J. (1997), “Summary report on surface-wave
project at Kansas Geological Survey (KGS)”, Lawrence, KS, Kansas
Geological Survey.
66
Park, C. B., Miller, R. D. & Xia, J. (1998). Imaging dispersion curves of surface
waves on multichannel record. 68th Annual International Meeting Society
of Exploration Geophysicists, Expanded Abstracts, 1377-1380.
Park, C. B., Miller, R. D. & Xia, J. (1999). Multichannel analysis of surface
waves. Geophysics, 64 (3), 800-808.
Park, C. B., Miller, R. D. & Miura, H. (2002). Optimum field parameters of an
MASW survey. Expanded Abstracts, SEG-J, Tokyo, May 22-23, 2002.
Park, C.B., Miller, D.R., Xia J., Ivanov J. (2007). Multichannel Analysis Of
Surface Waves (MASW)—active and passive methods. The Leading Edge,
Kansas Geological Survey, Lawrence : USA
R, A.P., Purwanto, M.S., Widodo, A. (2017). Identifikasi Percepatan Tanah
Maksimum (PGA) dan Kerentanan Tanah Menggunakan Metode
Mikrotremor Jalur Sesar Kendeng. J. Geosaintek 3, 107–114.
Rahmaningtyas, A.P., Widodo,A., Purwanto,M.S. (2017). Identifikasi
Percepatan Tanah Maksimum (PGA) Dan Kerentanan Tanah
Menggunakan Metode Mikrotremor I Jalur Sesar Kendeng.Surabaya.Jurnal
Sainsdan seni Vol 6 no.2 Pomits.
Rasimeng, Syamsurijal dkk. (2016). Analisis Kecepatan Gelombang Geser
(Vs30) Menggunakan Metode Seismik Multichannel Analysis Of Surface
Wave (Masw) Untuk Menentukan Resiko Bencana Gempabumi Di Kota
Bandar Lampung. Prosiding Seminar Nasional Kebencanaan Ikatan Ahli
Kebencanaan Indonesia Bandung.
Sheriff R.E.dan Geldart L.P. (1982), “Exploration seismology : History, theory
and data acquisition vol. 1”, Cambridge University press : London.
Riyantiyo,N.D,. Widodo,A,. Bahri,A.S. (2017), “Identifikasi Patahan Lokal
Menggunakan Metode Mikrotremor”. Jurnal Sainsdan seni Vol 6 no.2
Pomits, Surabaya.
Rošer, J. dan Gosar, A. (2010), “Determination of Vs30 for seismic ground
classifications in the Ljubljana area, Slovenia”, Acta Geotechnica Slovenia,
Slovenia
Ryden, N., Park, C. B., Ulriksen, P. & Miller, R. D. (2004), “Multimodal
approach to seismic pavement testing”, Journal of Geotechnical and
Geoenvironmental Engineering, 130, 636-645.
SIST EN 1998-1:2005/oA101. (2005). “Eurocode 8, design of structures for
earthquake resistance – part 1: general rules, seismic actions and rules for
67
buildings”, national Annex. Slovenian institute for standarization,
Ljubljana.
Sitharam TG, Anbazhagan P. (2006), “Measurements of dynamic properties and
soil profiling using multichannel analysis of surface waves”, Invited
keynote in the 4th Karl Terzaghi memorial workshops, Madgaon, Goa; 7
October 2006. p. 6.
Solihan, A., dan Santosa, B.J. (2009), “Analisis Dispersi Gelombang Rayleigh
Struktur Geologi Bawah Permukaan. Studi Kasus : Daerah pasir putih
dalegan Gresik”, Prosiding Seminar Nasional I\X Pasca Sarjana , FIS – 01,
Indonesia.
Sukardi, dkk. (1992), “Peta Geologi Lembar Surabaya dan Sapulu, Jawa”, Pusat
Penelitian dan Pengembangan Geologi (PPPG) , Bandung.
Sunardi, B., Daryono, Arifin, J., Susilanto, P., Ngadmanto, D., Nurdiyanto, B.,
Shahzad, S. (2012), “Kajian Potensi Bahaya Gempabumi Daerah Sumbawa
Berdasarkan Efek Tapak Lokal”, J. Meteorol. Dan Geofis. 13, 131–137.
Tim Pusat Gempa Nasional. (2017), “Peta Sumber dan Bahaya Gempa Indonesia
Tahun 2017”, Kementrian Pekerjaan Umum dan Perumahan Rakyat,
Indonesia.
Tokeshi, K., Harutoonian, p., Leo, C. J., Liyanapathirana, S. (2013), “Use of
surface waves for geotechnical engineering applications in Discussions
Western Sydney, Advances in Geosciences”, Copernicus Publisher.
Tuladhar, R. (2002), “Seismic Microzonation of Greater Bangkok Using
Microtremor Observations”, Thesis Asian Institute of Technology School of
Civil Engineering, Thailand
Vucetic, M. (1992), “Soil Properties and Seismic Response”, Earthquake
Engineering, Tenth World Conference, Balkerna, Rotterdam.
Xia, J., Miller, R. D. & Park, C. B. (1999), “Estimation of near-surface shear-
wave velocity by inversion of Rayleigh waves”, Geophysics, 64 (3), 691-
700.
Xia, J., Miller R. D., Park C. B. & Tian G. (2003), “Inversion of high frequency
surface waves with fundamental and higher modes”, Journal of Applied
Geophysics, 52, 45-57.
68
LAMPIRAN
69
HALAMAN INI SENGAJA DIKOSONGKAN
70
RIWAYAT PENULIS
71