Fizik Gelombang Seismik dan Selang S−P
Δt = D(1/Vs − 1/Vp)
D = Δt ÷ (1/Vs − 1/Vp)
d = R × θ, R = 6371 km
min Σ(dᵢ − Dᵢ)²
Penentuan lokasi episentrum gempa bumi bergantung pada sifat fizikal gelombang seismik yang merambat melalui jasad Bumi. Apabila sesuatu gempa bumi berlaku, ia membebaskan tenaga dalam bentuk gelombang jasad, terutamanya gelombang P (primer atau mampatan) dan gelombang S (sekunder atau ricih). Gelombang P bergerak lebih laju berbanding gelombang S dalam mana-mana medium batuan. Perbezaan kelajuan ini menyebabkan gelombang P sentiasa tiba dahulu di stesen pencatat, diikuti oleh gelombang S beberapa saat kemudian.
Apabila jarak dari punca gempa bumi ke stesen pencatat bertambah, selang masa antara ketibaan gelombang P dan gelombang S (selang S−P) juga akan meningkat secara berkadar terus. Dengan mengukur selang masa ini (Δt) pada seismogram, ahli seismologi boleh mengira anggaran jarak dari stesen tersebut ke episentrum gempa bumi. Di bawah andaian kelajuan gelombang yang malar, hubungan antara jarak (D) dan selang masa ketibaan memenuhi persamaan fizikal berikut:
Persamaan ini boleh disusun semula untuk mencari jarak anggaran dari stesen ke episentrum:
Geometri Penyegitigaan pada Model Bumi Sfera
Satu stesen seismik tunggal yang merekodkan selang S−P hanya dapat menentukan jarak radius sahaja, yang bermaksud episentrum boleh berada di mana-mana titik pada bulatan dengan jejari tersebut di sekeliling stesen. Apabila terdapat dua stesen, persilangan kedua-dua bulatan jarak mengecilkan kemungkinan lokasi kepada dua titik unik. Kehadiran stesen ketiga diperlukan untuk menyelesaikan kekaburan ini dan menentukan satu titik episentrum yang unik.
Kalkulator Penyegitigaan Episentrum Gempa Bumi menggunakan model Bumi sfera dengan jejari malar R = 6,371 km untuk mengira jarak permukaan antara koordinat stesen dan anggaran episentrum. Jarak permukaan (d) antara dua titik pada permukaan sfera dikira menggunakan formula sudut tengah haversine:
d = R × sudut tengah haversine, dengan R = 6,371 km
Dalam keadaan dunia sebenar, data pemerhatian daripada stesen seismik sering kali mengandungi ralat kecil akibat ketidakpastian masa ketibaan atau variasi struktur geologi. Oleh itu, bulatan-bulatan jarak daripada tiga atau lebih stesen jarang sekali bersilang tepat pada satu titik. Untuk mengatasi masalah ini, kaedah matematik kuasa dua terkecil digunakan untuk meminimumkan jumlah kuasa dua perbezaan antara jarak sfera yang disesuaikan (dᵢ) dan jarak yang dikira daripada selang S−P (Dᵢ) bagi setiap stesen:
Kesesuaian episentrum: minimumkan Σ(dᵢ − Dᵢ)² merentasi semua stesen
Parameter Input dan Batasan Model
Untuk menggunakan alat pengajaran ini, pengguna perlu memasukkan parameter pemerhatian stesen dan model perambatan yang khusus. Alat ini menyokong analisis menggunakan minimum 3 stesen sehingga maksimum 8 stesen.
Pemerhatian Stesen
Bagi setiap stesen (antara 3 hingga 8 stesen), pengguna perlu memasukkan data berikut:
- Latitud: Nilai mestilah dalam julat −90° hingga 90°.
- Longitud: Nilai mestilah dalam julat −180° hingga 180°.
- Selang S−P: Nilai dalam saat, mestilah lebih besar daripada 0 dan tidak melebihi 2,000 saat.
Pengguna boleh menekan butang Muat contoh Tokyo untuk mengisi alat dengan data sampel, atau menekan Kosongkan untuk menetapkan semula semua input.
Model Perambatan
- Kelajuan gelombang P: Mestilah lebih besar daripada 0 dan tidak melebihi 20 km/s, serta mestilah lebih besar daripada kelajuan gelombang S.
- Kelajuan gelombang S: Mestilah lebih besar daripada 0 dan tidak melebihi 20 km/s, serta mestilah kurang daripada kelajuan gelombang P.
- Ketidakpastian masa (1σ): Nilai dalam saat, mestilah lebih besar daripada 0 dan tidak melebihi 60 saat.
Batasan Model Pengajaran
Alat ini direka sebagai model pengajaran kerak seragam yang dipermudahkan. Pengiraan jarak S−P mengandaikan gempa bumi punca cetek dan kelajuan gelombang yang malar. Model ini mengabaikan kedalaman fokus (fokal), laluan sinar melengkung, struktur kelajuan berlapis atau tiga dimensi, kesilapan pemilihan fasa, dan pincang masa stesen. Untuk kerja seismologi profesional, ahli sains menggunakan model masa perjalanan serantau yang lebih kompleks untuk menghasilkan lokasi katalog seismik yang berautoriti.
Analisis Sisaan dan Geometri Stesen
Kualiti penyelesaian episentrum dinilai melalui analisis sisaan (residuals) dan kekuatan geometri susunan stesen. Sisaan bagi setiap stesen dikira sebagai jarak sfera yang disesuaikan minus jarak yang disimpulkan daripada selang S−P.
Sisaan positif menunjukkan bahawa episentrum yang disesuaikan berada lebih jauh dari stesen berbanding jejari S−P yang dikira. Sebaliknya, sisaan negatif menunjukkan episentrum yang disesuaikan adalah lebih dekat dengan stesen tersebut. Sisaan RMS (Root Mean Square) memberikan gambaran keseluruhan tentang kejituan kesesuaian model.
Kekuatan geometri stesen dipengaruhi oleh taburan spatial stesen di sekeliling episentrum. Jika stesen-stesen terletak dalam satu garisan lurus atau berkelompok rapat, penyelesaian matematik menjadi tidak stabil kerana ia tidak dapat mengekang lokasi dua dimensi dengan baik. Taburan stesen yang mengelilingi episentrum dari pelbagai arah akan menghasilkan geometri stesen yang kuat.
Alat ini mengklasifikasikan kualiti penyelesaian berdasarkan kriteria berikut:
- Kesesuaian jarak: Dikategorikan sebagai kesesuaian jarak rapat, kesesuaian jarak bercampur, atau kesesuaian jarak lemah.
- Geometri stesen: Dikategorikan sebagai geometri stesen kuat, geometri stesen sederhana, atau geometri stesen lemah.
Perambatan Ralat dan Ketidakpastian Kedudukan
Ketidakpastian dalam menentukan masa ketibaan gelombang pada seismogram (ketidakpastian masa 1σ) akan merambat menjadi ketidakpastian kedudukan spatial episentrum. Ralat masa ini ditukarkan kepada ralat jarak menggunakan faktor jarak S−P yang diperoleh daripada kelajuan gelombang yang ditetapkan.
Seterusnya, geometri susunan stesen akan menentukan bagaimana ralat jarak ini merambat ke arah mendatar. Alat ini mengira anggaran ketidakpastian kedudukan 1σ sepanjang paksi mengufuk utama. Jika geometri stesen adalah lemah, ralat kedudukan pada paksi tertentu akan menjadi sangat besar walaupun ralat masa ketibaan adalah kecil.
Aturan Pengesahan dan Mesej Ralat
Alat ini melaksanakan aturan pengesahan input yang ketat untuk memastikan kestabilan pengiraan. Jika terdapat ralat pada input, mesej ralat berikut akan dipaparkan:
- Had Bilangan Stesen: Jika stesen kurang daripada 3 atau lebih daripada 8, alat memaparkan:
Gunakan antara tiga hingga lapan stesen. - Turutan Kelajuan Gelombang Salah:
Kelajuan gelombang P mestilah lebih besar daripada kelajuan gelombang S. - Ketidakpastian Masa Luar Julat:
Ketidakpastian masa mestilah lebih besar daripada 0 dan tidak melebihi 60 saat. - Geometri Tidak Stabil:
Stesen-stesen ini tidak mengekang lokasi dua dimensi yang stabil. Tambah stesen yang jauh daripada garisan atau kelompok semasa. - Limpahan Nombor:
Nilai atau keputusan perantara melebihi julat nombor yang disokong.
Apabila pengiraan berjaya dijalankan, status kaedah akan memaparkan salah satu daripada mesej berikut:
Edit pemerhatian stesen untuk menyesuaikan semula episentrum.
Privasi Data dan Pemprosesan
Semua data yang dimasukkan ke dalam kalkulator ini, termasuk koordinat stesen, selang masa S−P, parameter kelajuan gelombang, dan keputusan pengiraan yang dihasilkan, diproses secara tempatan di dalam pelayar web peranti pengguna. Tiada data yang dihantar atau dimuat naik ke pelayan luaran.
Soalan Lazim (FAQ)
Mengapakah tiga stesen diperlukan? Satu pengukuran S−P memberikan bulatan jarak, bukan arah. Dua bulatan boleh bertemu di dua tempat. Stesen ketiga biasanya menyelesaikan pilihan tersebut, manakala stesen keempat atau seterusnya menjadikan ketidakpadanan kelihatan melalui sisaan. Stesen yang tersebar di sekeliling punca mengekang kedua-dua latitud dan longitud dengan lebih baik berbanding stesen di sepanjang satu garisan.
Adakah jarak S−P merupakan jarak permukaan yang tepat? Tidak. Persamaan kelajuan malar menganggarkan laluan perjalanan biasa daripada perbezaan masa ketibaan. Menganggap jejari tersebut sebagai jarak permukaan sfera mengandaikan punca cetek dan struktur kelajuan yang mudah. Keluk masa perjalanan sebenar bergantung pada kedalaman, jarak dan bahan yang dilalui oleh setiap sinar.
Apakah maksud sisaan positif atau negatif? Sisaan ialah jarak sfera yang disesuaikan tolak jarak yang disimpulkan daripada S−P. Positif bermakna titik yang disesuaikan adalah lebih jauh dari stesen tersebut berbanding jejari S−P miliknya; negatif bermakna ia adalah lebih dekat. Sisaan bercampur yang besar menunjukkan bahawa bulatan-bulatan tersebut tidak berkongsi satu persilangan yang bersih di bawah kelajuan yang dipilih.
Apakah yang termasuk dalam ketidakpastian kedudukan? Ia merambatkan nilai tunggal ketidakpastian masa 1σ melalui faktor jarak S−P dan geometri stesen berhampiran titik yang disesuaikan. Ia tidak mengambil kira model kelajuan yang salah, kedalaman fokus, fasa yang salah label, pincang jam atau struktur Bumi tiga dimensi, jadi ia bukan tuntutan ketepatan yang lengkap.