ഭൂകമ്പ എപിസെന്റർ ട്രയാംഗുലേഷൻ കാൽക്കുലേറ്റർ
ഭൂകമ്പ തരംഗങ്ങളുടെ പ്രവേഗവും വിവിധ നിരീക്ഷണ സ്റ്റേഷനുകളിൽ അവ രേഖപ്പെടുത്തിയ സമയവും അടിസ്ഥാനമാക്കി ഭൂകമ്പത്തിന്റെ പ്രഭവകേന്ദ്രം (Epicenter) കണ്ടെത്താൻ സഹായിക്കുന്ന ഒരു വിദ്യാഭ്യാസ ഉപകരണമാണ് ഭൂകമ്പ എപിസെന്റർ ട്രയാംഗുലേഷൻ കാൽക്കുലേറ്റർ. ഗോളീയ മാതൃകയിലുള്ള ഭൂമിയിൽ മൂന്ന് മുതൽ എട്ട് വരെ സ്റ്റേഷനുകളിൽ നിന്നുള്ള വിവരങ്ങൾ ഉപയോഗിച്ച് ലീസ്റ്റ്-സ്ക്വയേഴ്സ് (Least-squares) രീതിയിലൂടെ ഏറ്റവും അനുയോജ്യമായ എപിസെന്റർ ലൊക്കേഷൻ കണ്ടെത്താൻ ഈ ഉപകരണം സഹായിക്കുന്നു.
ഉപയോക്താക്കൾ നൽകുന്ന സ്റ്റേഷൻ കോർഡിനേറ്റുകൾ, S−P തരംഗങ്ങളുടെ ആഗമന സമയ വ്യത്യാസം (S−P interval), തരംഗ വേഗതകൾ, ടൈമിംഗ് അനിശ്ചിതത്വം എന്നിവ വിശകലനം ചെയ്താണ് കാൽക്കുലേറ്റർ ഫലം കണക്കാക്കുന്നത്. ഈ പ്രക്രിയ പൂർണ്ണമായും ഉപയോക്താവിന്റെ ബ്രൗസറിൽ മാത്രമാണ് നടക്കുന്നത്; വിവരങ്ങൾ യാതൊരുവിധ ബാഹ്യ സെർവറുകളിലേക്കും അപ്ലോഡ് ചെയ്യപ്പെടുന്നില്ല.
സീസ്മിക് തരംഗങ്ങളുടെ ഭൗതികശാസ്ത്രവും ട്രയാംഗുലേഷനും
ഭൂകമ്പം ഉണ്ടാകുമ്പോൾ വിവിധ തരത്തിലുള്ള സീസ്മിക് തരംഗങ്ങൾ പുറപ്പെടുന്നു. ഇതിൽ പ്രൈമറി തരംഗങ്ങളും (P-waves) സെക്കൻഡറി തരംഗങ്ങളും (S-waves) പ്രധാനപ്പെട്ടവയാണ്. P-തരംഗങ്ങൾ S-തരംഗങ്ങളേക്കാൾ വേഗത്തിൽ സഞ്ചരിക്കുന്നു. അതിനാൽ, പ്രഭവകേന്ദ്രത്തിൽ നിന്ന് അകലം കൂടുംതോറും ഈ രണ്ട് തരംഗങ്ങളും സ്റ്റേഷനിൽ എത്തുന്ന സമയങ്ങൾ തമ്മിലുള്ള വ്യത്യാസം (S−P interval) വർദ്ധിച്ചുകൊണ്ടിരിക്കുന്നു.
ഒരു നിരീക്ഷണ സ്റ്റേഷനിലെ S−P സമയ വ്യത്യാസം ഉപയോഗിച്ച് ആ സ്റ്റേഷനിൽ നിന്ന് ഭൂകമ്പ കേന്ദ്രത്തിലേക്കുള്ള ദൂരം (ആരം) കണക്കാക്കാം. എന്നാൽ ഒരു സ്റ്റേഷനിൽ നിന്നുള്ള ദൂരം കൊണ്ട് മാത്രം കൃത്യമായ ദിശ കണ്ടെത്താൻ കഴിയില്ല. രണ്ട് സ്റ്റേഷനുകളിലെ വിവരങ്ങൾ ലഭ്യമാകുമ്പോൾ, രണ്ട് വൃത്തങ്ങൾ പരസ്പരം കൂട്ടിമുട്ടുന്ന രണ്ട് സാധ്യതാ ബിന്ദുക്കൾ ലഭിക്കുന്നു. മൂന്നാമതൊരു സ്റ്റേഷന്റെ വിവരങ്ങൾ കൂടി ലഭിക്കുന്നതോടെ, ഈ വൃത്തങ്ങൾ പൊതുവായി സന്ധിക്കുന്ന കൃത്യമായ ഒരു പ്രഭവകേന്ദ്രം (എപിസെന്റർ) കണ്ടെത്താൻ സാധിക്കുന്നു. പ്രൊഫഷണൽ വിശകലനങ്ങളിൽ കൂടുതൽ സ്റ്റേഷനുകൾ ഉപയോഗിക്കുന്നത് കൃത്യത വർദ്ധിപ്പിക്കാൻ സഹായിക്കുന്നു.
ഇൻപുട്ട് പാരാമീറ്ററുകളും പരിധികളും
കൃത്യമായ കണക്കുകൂട്ടലിനായി ഉപയോക്താക്കൾ താഴെ പറയുന്ന വിവരങ്ങൾ നൽകേണ്ടതുണ്ട്:
- സ്റ്റേഷൻ നിരീക്ഷണങ്ങൾ (കുറഞ്ഞത് 3, പരമാവധി 8 സ്റ്റേഷനുകൾ):
- സ്റ്റേഷന്റെ പേര്
- അക്ഷാംശം: −90° മുതൽ 90° വരെ ആയിരിക്കണം.
- രേഖാംശം: −180° മുതൽ 180° വരെ ആയിരിക്കണം.
- S−P ഇന്റർവെൽ: 0-ൽ കൂടുതലും 2,000 സെക്കൻഡിൽ കൂടാത്തതുമായിരിക്കണം.
- പ്രൊപഗേഷൻ മോഡൽ:
- P-തരംഗ വേഗത: 0-ൽ കൂടുതലും 20 km/s-ൽ കൂടാത്തതുമായിരിക്കണം. ഇത് എപ്പോഴും S-തരംഗ വേഗതയേക്കാൾ കൂടുതലായിരിക്കണം.
- S-തരംഗ വേഗത: 0-ൽ കൂടുതലും 20 km/s-ൽ കൂടാത്തതുമായിരിക്കണം. ഇത് എപ്പോഴും P-തരംഗ വേഗതയേക്കാൾ കുറവായിരിക്കണം.
- ടൈമിംഗ് അനിശ്ചിതത്വം (1σ): 0-ൽ കൂടുതലും 60 സെക്കൻഡിൽ കൂടാത്തതുമായിരിക്കണം.
ഉപയോക്താക്കൾക്ക് മാതൃകാ വിവരങ്ങൾ പരിശോധിക്കുന്നതിനായി ടോക്കിയോ ഉദാഹരണം ലോഡ് ചെയ്യുക എന്ന ബട്ടണും, ഇൻപുട്ടുകൾ പുനഃക്രമീകരിക്കുന്നതിനായി മായ്ക്കുക എന്ന ബട്ടണും ഉപയോഗിക്കാവുന്നതാണ്.
കണക്കുകൂട്ടൽ രീതികളും സമവാക്യങ്ങളും
ഈ ഉപകരണം ഭൂമിയെ 6,371 കി.മീറ്ററോളം സ്ഥിരമായ ആരമുള്ള (R = 6,371 km) ഒരു ഗോളമായി കണക്കാക്കുന്നു. ആഴം കുറഞ്ഞ ഉറവിടമുള്ള ഭൂകമ്പങ്ങളും സ്ഥിരമായ തരംഗ വേഗതയും അനുമാനിച്ചാണ് കണക്കുകൂട്ടലുകൾ നടത്തുന്നത്.
ഉപയോഗിക്കുന്ന പ്രധാന സമവാക്യങ്ങൾ താഴെ പറയുന്നവയാണ്:
- ദൂര ഏകദേശ കണക്ക്: D = Δt ÷ (1/Vs − 1/Vp)
ഇവിടെ
Δtഎന്നത് S−P സമയ വ്യത്യാസവും,Vp,Vsഎന്നിവ യഥാക്രമം P, S തരംഗങ്ങളുടെ വേഗതയുമാണ്. - ഉപരിതല ദൂരം: d = R × ഹാവർസൈൻ സെൻട്രൽ ആംഗിൾ, R = 6371 km.
- എപിസെന്റർ ഫിറ്റ്: എല്ലാ സ്റ്റേഷനുകളിലുമായി Σ(dᵢ − Dᵢ)² കുറയ്ക്കുക
യഥാർത്ഥ നിരീക്ഷണങ്ങളിലെ ചെറിയ വ്യത്യാസങ്ങൾ പരിഹരിക്കാൻ ലീസ്റ്റ്-സ്ക്വയേഴ്സ് രീതി ഉപയോഗിച്ച് സ്റ്റേഷനുകളിലെ ദൂര വ്യത്യാസങ്ങളുടെ വർഗ്ഗങ്ങളുടെ തുക ഏറ്റവും കുറഞ്ഞ അളവിലേക്ക് പരിമിതപ്പെടുത്തി മികച്ച ഫിറ്റ് കണ്ടെത്തുന്നു.
ഔട്ട്പുട്ടുകളും ഫല വിശകലനവും
കണക്കുകൂട്ടൽ പൂർത്തിയാകുമ്പോൾ താഴെ പറയുന്ന ഫലങ്ങൾ ലഭ്യമാകും:
- കണക്കാക്കിയ എപിസെന്റർ: കണക്കാക്കിയ അക്ഷാംശ, രേഖാംശ കോർഡിനേറ്റുകൾ.
- ദൂര ഫിറ്റ് ഗുണനിലവാരം: ഇതിനെ അടുത്ത ദൂര ഫിറ്റ്, മിശ്രിത ദൂര ഫിറ്റ്, അല്ലെങ്കിൽ മോശം ദൂര ഫിറ്റ് എന്ന് തരംതിരിക്കുന്നു.
- സ്റ്റേഷൻ ജ്യാമിതി ഗുണനിലവാരം: ഇതിനെ ശക്തമായ സ്റ്റേഷൻ ജ്യാമിതി, മിതമായ സ്റ്റേഷൻ ജ്യാമിതി, അല്ലെങ്കിൽ ദുർബലമായ സ്റ്റേഷൻ ജ്യാമിതി എന്ന് രേഖപ്പെടുത്തുന്നു.
- റെസിഡ്യുവലുകൾ: ഫിറ്റ് ചെയ്ത ഗോളീയ ദൂരത്തിൽ നിന്ന് S−P-യിൽ നിന്ന് അനുമാനിച്ച ദൂരം കുറയ്ക്കുന്നതാണ് റെസിഡ്യുവൽ. പോസിറ്റീവ് റെസിഡ്യുവൽ എന്നാൽ ഫിറ്റ് ചെയ്ത എപിസെന്റർ സ്റ്റേഷനിൽ നിന്ന് അതിന്റെ S−P ആരത്തേക്കാൾ അകലെയാണെന്നും, നെഗറ്റീവ് എന്നാൽ അത് കൂടുതൽ അടുത്താണെന്നും അർത്ഥമാക്കുന്നു.
കൂടാതെ, സ്റ്റേഷനുകളുടെ ആപേക്ഷിക സ്ഥാനം കാണിക്കുന്ന ഒരു സ്റ്റേഷൻ ജ്യാമിതി സ്കീമാറ്റിക് ചിത്രവും ഇതിലൊപ്പം ഉണ്ടാകും. ഇതിന് താഴെയായി "ആപേക്ഷിക ജ്യാമിതി മാത്രം — മാപ്പോ അതിരുകളോ തത്സമയ ഭൂകമ്പ വിവരങ്ങളോ ഇല്ല." എന്ന അടിക്കുറിപ്പ് നൽകിയിരിക്കുന്നു.
പിശക് സന്ദേശങ്ങളും പരിഹാരങ്ങളും
ഇൻപുട്ടുകളിൽ തെറ്റായ വിവരങ്ങൾ നൽകുമ്പോൾ താഴെ പറയുന്ന പിശക് സന്ദേശങ്ങൾ കാണിക്കുന്നു:
- സ്റ്റേഷനുകളുടെ എണ്ണം പരിധി ലംഘിച്ചാൽ:
മൂന്ന് മുതൽ എട്ട് വരെ സ്റ്റേഷനുകൾ ഉപയോഗിക്കുക. - സംഖ്യയല്ലാത്ത വിവരങ്ങൾ നൽകിയാൽ:
input: “abc” എന്നത് ഒരു സംഖ്യയല്ല. - അക്ഷാംശം തെറ്റായാൽ:
സ്റ്റേഷൻ 1: അക്ഷാംശം −90° മുതൽ 90° വരെ ആയിരിക്കണം. - രേഖാംശം തെറ്റായാൽ:
സ്റ്റേഷൻ 1: രേഖാംശം −180° മുതൽ 180° വരെ ആയിരിക്കണം. - S−P ഇന്റർവെൽ തെറ്റായാൽ:
സ്റ്റേഷൻ 1: 0-ൽ കൂടുതലും 2000 സെക്കൻഡിൽ കൂടാത്തതുമായ ഒരു S−P ഇന്റർവെൽ നൽകുക. - വേഗത പരിധി ലംഘിച്ചാൽ:
Vₚ/Vₛ 0-ൽ കൂടുതലും 20 km/s-ൽ കൂടാത്തതുമായിരിക്കണം. - വേഗതയുടെ ക്രമം തെറ്റായാൽ:
P-തരംഗ വേഗത S-തരംഗ വേഗതയേക്കാൾ കൂടുതലായിരിക്കണം. - ടൈമിംഗ് അനിശ്ചിതത്വം തെറ്റായാൽ:
ടൈമിംഗ് അനിശ്ചിതത്വം 0-ൽ കൂടുതലും 60 സെക്കൻഡിൽ കൂടാത്തതുമായിരിക്കണം. - ഒരേ കോർഡിനേറ്റുകൾ നൽകിയാൽ:
1, 2 എന്നീ സ്റ്റേഷനുകൾക്ക് ഒരേ കോർഡിനേറ്റുകളാണുള്ളത്. വ്യത്യസ്ത സ്റ്റേഷൻ ലൊക്കേഷനുകൾ ഉപയോഗിക്കുക. - സ്റ്റേഷനുകൾ ഒരേ വരിയിലോ കൂട്ടമായോ വന്നാൽ:
ഈ സ്റ്റേഷനുകൾ സ്ഥിരതയുള്ള ഒരു ദ്വിമാന ലൊക്കേഷൻ നിർണ്ണയിക്കുന്നില്ല. നിലവിലെ ലൈനിൽ നിന്നോ ക്ലസ്റ്ററിൽ നിന്നോ മാറ്റി ഒരു സ്റ്റേഷൻ ചേർക്കുക. - സംഖ്യാ പരിധി കവിഞ്ഞാൽ:
ഒരു മൂല്യമോ ഇടനില ഫലമോ പിന്തുണയ്ക്കുന്ന സംഖ്യാ പരിധിയേക്കാൾ കൂടുതലാണ്.
രീതിയും വ്യാപ്തിയും
ഈ ഉപകരണം യുഎസ്ജിഎസ് (USGS) വിദ്യാഭ്യാസ റഫറൻസുകൾ പിന്തുടർന്നാണ് രൂപകൽപ്പന ചെയ്തിരിക്കുന്നത്. ഇതിന്റെ രീതിയും വിവരങ്ങളും പരിശോധിച്ച തീയതി 6 ഓഗസ്റ്റ് 2026 ആണ്.
ഇതൊരു വിദ്യാഭ്യാസ ആവശ്യത്തിനായുള്ള ലളിതമായ മാതൃക മാത്രമാണ്. യഥാർത്ഥ ഭൂകമ്പ തരംഗങ്ങളുടെ സഞ്ചാരം ഭൂമിയുടെ ആഴം, പാളികളുള്ള ഘടന, വളഞ്ഞ തരംഗ പാതകൾ എന്നിവയെ ആശ്രയിച്ചിരിക്കുന്നു. അതിനാൽ പ്രൊഫഷണൽ ആവശ്യങ്ങൾക്കായി ഔദ്യോഗിക സീസ്മിക് കാറ്റലോഗ് വിവരങ്ങൾ ഉപയോഗിക്കേണ്ടതാണ്.
പതിവായി ചോദിക്കുന്ന ചോദ്യങ്ങൾ (FAQ)
ചോദ്യം: S−P ദൂരം എന്നത് കൃത്യമായ ഉപരിതല ദൂരമാണോ?
ഉത്തരം: അല്ല. സ്ഥിരമായ വേഗതയുള്ള സമവാക്യം, ആഗമന സമയ വ്യത്യാസത്തിൽ നിന്ന് ഒരു പൊതു യാത്രാ പാത കണക്കാക്കുകയാണ് ചെയ്യുന്നത്. ആ ആരത്തെ ഗോളീയ ഉപരിതല ദൂരമായി കണക്കാക്കുന്നത് ആഴം കുറഞ്ഞ ഉറവിടവും ലളിതമായ വേഗത ഘടനയും അനുമാനിക്കുന്നു. യഥാർത്ഥ യാത്രാ സമയ വക്രങ്ങൾ ആഴം, ദൂരം, ഓരോ തരംഗവും കടന്നുപോകുന്ന മാധ്യമം എന്നിവയെ ആശ്രയിച്ചിരിക്കുന്നു.
ചോദ്യം: എന്തുകൊണ്ടാണ് മൂന്ന് സ്റ്റേഷനുകൾ ആവശ്യമായി വരുന്നത്?
ഉത്തരം: ഒരു S−P അളവ് ഒരു ദൂര വൃത്തം മാത്രമാണ് നൽകുന്നത്, ദിശയല്ല. രണ്ട് വൃത്തങ്ങൾക്ക് രണ്ട് സ്ഥലങ്ങളിൽ പരസ്പരം കൂട്ടിമുട്ടാൻ കഴിഞ്ഞ് രണ്ട് സാധ്യതാ ബിന്ദുക്കൾ നൽകുന്നു. ഒരു മൂന്നാമത്തെ സ്റ്റേഷൻ സാധാരണയായി ശരിയായ ഓപ്ഷൻ കണ്ടെത്താൻ സഹായിക്കുന്നു, അതേസമയം നാലാമതോ അതിനുശേഷമുള്ളതോ ആയ സ്റ്റേഷൻ റെസിഡ്യുവലുകളിലൂടെ പൊരുത്തക്കേടുകൾ ദൃശ്യമാക്കുന്നു. ഉറവിടത്തിന് ചുറ്റും വ്യാപിച്ചുകിടക്കുന്ന സ്റ്റേഷനുകൾ, ഒരു വരിയിലുള്ള സ്റ്റേഷനുകളേക്കാൾ നന്നായി അക്ഷാംശവും രേഖാംശവും നിർണ്ണയിക്കാൻ സഹായിക്കുന്നു.
ചോദ്യം: പോസിറ്റീവ് അല്ലെങ്കിൽ നെഗറ്റീവ് റെസിഡ്യുവൽ എന്നാൽ എന്താണ് അർത്ഥമാക്കുന്നത്?
ഉത്തരം: ഫിറ്റ് ചെയ്ത ഗോളീയ ദൂരത്തിൽ നിന്ന് S−P-യിൽ നിന്ന് അനുമാനിച്ച ദൂരം കുറയ്ക്കുന്നതാണ് റെസിഡ്യുവൽ. പോസിറ്റീവ് എന്നാൽ ഫിറ്റ് ചെയ്ത ബിന്ദു ആ സ്റ്റേഷനിൽ നിന്ന് അതിന്റെ S−P ആരത്തേക്കാൾ അകലെയാണെന്നാണ്; നെഗറ്റീവ് എന്നാൽ അത് കൂടുതൽ അടുത്താണെന്നാണ്. വലിയ മിശ്രിത റെസിഡ്യുവലുകൾ കാണിക്കുന്നത് തിരഞ്ഞെടുത്ത വേഗതകളിൽ വൃത്തങ്ങൾക്ക് കൃത്യമായ ഒരു പൊതുവായ വിഭജനം പങ്കിടാൻ കഴിയില്ലെന്നാണ്.
ചോദ്യം: സ്ഥാനത്തിന്റെ അനിശ്ചിതത്വത്തിൽ എന്താണ് ഉൾപ്പെടുത്തിയിരിക്കുന്നത്?
ഉത്തരം: ഇത് ഒരൊറ്റ ടൈമിംഗ് 1σ മൂല്യത്തെ S−P ദൂര ഘടകത്തിലൂടെയും ഫിറ്റ് ചെയ്ത ബിന്ദുവിന് സമീപമുള്ള സ്റ്റേഷൻ ജ്യാമിതിയിലൂടെയും വ്യാപിപ്പിക്കുന്നു. ഇതിൽ തെറ്റായ വേഗത മോഡൽ, ഫോക്കൽ ആഴം, തെറ്റായി അടയാളപ്പെടുത്തിയ ഫേസുകൾ, ക്ലോക്ക് ബയാസ് അല്ലെങ്കിൽ ത്രിമാന ഭൂഘടന എന്നിവ ഉൾപ്പെടുന്നില്ല, അതിനാൽ ഇത് ഒരു സമ്പൂർണ്ണ കൃത്യത അവകാശവാദമല്ല.