ভূকম্পন কেন্দ্র ট্রায়াঙ্গুলেশন ক্যালকুলেটর

ভূমিকম্পের কেন্দ্র নির্ণয় করা সিসমোলজির একটি মৌলিক প্রক্রিয়া। "Earthquake Epicenter Triangulation Calculator" একটি নিখরচায় উপলব্ধ অনলাইন শিক্ষামূলক টুল, যা ব্যবহারকারীদের একটি গোলাকার পৃথিবী মডেলে তিন থেকে আটটি রেকর্ডিং স্টেশনের পর্যবেক্ষণ ব্যবহার করে ভূমিকম্পের কেন্দ্র বা এপিসেন্টার অনুমান করার সুযোগ দেয়। প্রতিটি স্টেশনের স্থানাঙ্ক এবং S−P তরঙ্গের আগমন সময়ের ব্যবধান প্রবেশ করিয়ে, এবং সেই সাথে তরঙ্গের গতিবেগ ও টাইমিং অনিশ্চয়তা নির্ধারণ করে ব্যবহারকারীরা একটি সর্বোত্তম ফিট করা ভূকম্পন কেন্দ্র অবস্থান হিসাব করতে পারেন।

এই টুলটি সমস্ত স্টেশন জুড়ে দূরত্বের অমিলগুলো সর্বনিম্ন করতে একটি ন্যূনতম-বর্গের (least-squares) হিসাব ব্যবহার করে। এটি হিসাবের পর আনুমানিক ভূকম্পন কেন্দ্র, স্টেশন জ্যামিতির শক্তি এবং প্রতিটি স্টেশনের অবশিষ্ট মান (residuals) প্রদর্শন করে দেখায় যে পর্যবেক্ষণগুলো কতটা সামঞ্জস্যপূর্ণ।

সিসমিক তরঙ্গের পদার্থবিজ্ঞান এবং S−P ব্যবধান

ভূমিকম্পের সময় বিভিন্ন ধরনের সিসমিক তরঙ্গ উৎপন্ন হয়। এর মধ্যে P-তরঙ্গ (প্রাথমিক বা সংকোচনশীল তরঙ্গ) এবং S-তরঙ্গ (দ্বিতীয় বা শিয়ার তরঙ্গ) অন্যতম। P-তরঙ্গ সবসময় S-তরঙ্গের চেয়ে দ্রুত ভ্রমণ করে। এর ফলে, ভূমিকম্পের উৎস থেকে দূরত্ব যত বাড়ে, এই দুটি তরঙ্গের স্টেশনে পৌঁছানোর সময়ের ব্যবধান (S−P interval) তত বৃদ্ধি পায়।

টুলটির ধ্রুব-গতিবেগ অনুমানের অধীনে, আগমন-সময়ের পার্থক্য নিচের সমীকরণটি মেনে চলে:

Δt = D(1/Vs − 1/Vp)

এই সমীকরণটিকে পুনর্বিন্যাস করলে দূরত্বের আনুমানিক হিসাব পাওয়া যায়:

D = Δt ÷ (1/Vs − 1/Vp)

এখানে D হলো স্টেশন থেকে দূরত্ব, Δt হলো S−P ব্যবধান, Vp হলো P-তরঙ্গের গতিবেগ এবং Vs হলো S-তরঙ্গের গতিবেগ।

ট্রায়াঙ্গুলেশনের জ্যামিতি এবং ন্যূনতম-বর্গের অপ্টিমাইজেশন

একটি একক স্টেশন থেকে প্রাপ্ত S−P ব্যবধান কেবল একটি দূরত্বের ব্যাসার্ধ প্রদান করে, যা স্টেশনের চারপাশে একটি বৃত্ত তৈরি করে। এই বৃত্তের যেকোনো বিন্দুতে ভূমিকম্পের কেন্দ্রটি অবস্থিত হতে পারে। দুটি স্টেশন ব্যবহার করলে তাদের বৃত্ত দুটি দুটি বিন্দুতে পরস্পরকে ছেদ করে। একটি অনন্য ভূকম্পন কেন্দ্র নির্ধারণের জন্য কমপক্ষে তিনটি স্টেশনের প্রয়োজন হয়।

বাস্তব সিসমিক ডেটাতে পরিমাপের ত্রুটি এবং ভূত্বকের অসমতার কারণে তিনটি বৃত্ত সবসময় একটি নিখুঁত বিন্দুতে মিলিত হয় না। এই অমিল সমাধানের জন্য ক্যালকুলেটরটি ন্যূনতম-বর্গের অপ্টিমাইজেশন পদ্ধতি ব্যবহার করে। এটি এমন একটি বিন্দু খুঁজে বের করে যেখানে সমস্ত স্টেশন জুড়ে গোলাকার দূরত্ব এবং S−P থেকে প্রাপ্ত দূরত্বের পার্থক্যের বর্গের সমষ্টি সর্বনিম্ন হয়।

ইনপুট এবং আউটপুট প্যারামিটার

ব্যবহারকারীকে টুলে নিম্নলিখিত ইনপুটগুলো প্রদান করতে হয়:

• স্টেশন পর্যবেক্ষণ (৩ থেকে ৮টি স্টেশন):
  • স্টেশনের নাম
  • অক্ষাংশ: অবশ্যই −90° থেকে 90° এর মধ্যে হতে হবে।
  • দ্রাঘিমাংশ: অবশ্যই −180° থেকে 180° এর মধ্যে হতে হবে।
  • S−P ব্যবধান: অবশ্যই 0-এর বেশি এবং 2,000 সেকেন্ডের বেশি নয় এমন হতে হবে।
• প্রসারণ মডেল প্যারামিটার:
  • P-তরঙ্গের গতিবেগ: অবশ্যই 0-এর বেশি এবং 20 km/s-এর বেশি হওয়া যাবে না, এবং S-তরঙ্গের গতিবেগের চেয়ে বেশি হতে হবে।
  • S-তরঙ্গের গতিবেগ: অবশ্যই 0-এর বেশি এবং 20 km/s-এর বেশি হওয়া যাবে না, এবং P-তরঙ্গের গতিবেগের চেয়ে কম হতে হবে।
  • টাইমিং অনিশ্চয়তা (1σ): অবশ্যই 0-এর বেশি এবং 60 সেকেন্ডের বেশি হওয়া যাবে না।

ব্যবহারকারীরা চাইলে টোকিওর উদাহরণ লোড করুন-এ ক্লিক করে নমুনা ডেটা দেখতে পারেন, অথবা ইনপুটগুলো রিসেট করতে মুছে ফেলুন বোতাম ব্যবহার করতে পারেন।

হিসাব সম্পন্ন হওয়ার পর টুলটি নিম্নলিখিত ফলাফলগুলো প্রদর্শন করে:

• স্টেশন পর্যবেক্ষণ বিভাগে প্রতিটি স্টেশনের জন্য হিসাবকৃত ব্যাসার্ধ দেখানো হয়।
• ভূকম্পন কেন্দ্রের সমাধান বিভাগে আনুমানিক ভূকম্পন কেন্দ্র (অক্ষাংশ ও দ্রাঘিমাংশ), RMS দূরত্বের অবশিষ্ট মান, বৃহত্তম দূরত্বের অবশিষ্ট মান, আনুমানিক অবস্থানের 1σ, এবং প্রতি স্টেশনে দূরত্বের 1σ প্রদর্শিত হয়।
• ফিটের গুণমান অনুযায়ী নিকটবর্তী দূরত্বের ফিট, মিশ্র দূরত্বের ফিট, অথবা দুর্বল দূরত্বের ফিট লেবেল প্রদর্শিত হয়।
• স্টেশন বিন্যাসের ওপর ভিত্তি করে শক্তিশালী স্টেশন জ্যামিতি, মাঝারি স্টেশন জ্যামিতি, অথবা দুর্বল স্টেশন জ্যামিতি লেবেল দেখানো হয়।
• স্টেশন জ্যামিতির রেখাচিত্র অংশে একটি ভিজ্যুয়াল উপস্থাপনা থাকে যার নিচে লেখা থাকে: "কেবল আপেক্ষিক জ্যামিতি — কোনো মানচিত্র, সীমানা বা লাইভ ভূমিকম্পের ডেটা নেই।"
• স্টেশন ফিট টেবিলে স্টেশন, S−P দূরত্ব, গোলাকার দূরত্ব, এবং অবশিষ্ট মান (Residual) কলামগুলো থাকে।

গাণিতিক সূত্র এবং প্রতিস্থাপন

ক্যালকুলেটরটি তার হিসাবের ধাপগুলো সূত্র এবং মান প্রতিস্থাপন বিভাগে বিস্তারিতভাবে প্রদর্শন করে:

• দূরত্বের আনুমানিক হিসাব: D = Δt ÷ (1/Vs − 1/Vp)
• পৃষ্ঠের দূরত্ব: d = R × হ্যাভারসাইন কেন্দ্রীয় কোণ, যেখানে R = 6371 km
• ভূকম্পন কেন্দ্র ফিট: সমস্ত স্টেশন জুড়ে Σ(dᵢ − Dᵢ)² সর্বনিম্ন করা

এছাড়াও এটি দূরত্ব গুণক হিসাব দেখায়: দূরত্ব গুণক = vp × vs ÷ (vp − vs) = প্রতি সেকেন্ড S−P-এর জন্য factor কিমি

প্রতিটি স্টেশনের জন্য দূরত্ব নির্ধারণ: station: interval সেকেন্ড × factor কিমি/সেকেন্ড = distance কিমি

এবং টাইমিং অনিশ্চয়তা প্রসারণ: timing সেকেন্ডের টাইমিং 1σ দেয় distance কিমি দূরত্বের 1σ; স্টেশন জ্যামিতি এটিকে প্রধান অনুভূমিক অক্ষে প্রায় position কিমি পর্যন্ত প্রসারিত করে।

ত্রুটি এবং বৈধতা নিয়ম

টুলটি ব্যবহারের সময় ইনপুট সংক্রান্ত কোনো অসঙ্গতি থাকলে নির্দিষ্ট ত্রুটি বার্তা প্রদর্শিত হয়:

• স্টেশন সংখ্যা ৩ থেকে ৮ এর বাইরে হলে: তিনটি থেকে আটটি স্টেশন ব্যবহার করুন।
• সংখ্যা ব্যতীত অন্য কোনো মান প্রবেশ করালে: 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 স্টেশনের স্থানাঙ্ক একই। ভিন্ন ভিন্ন স্টেশনের অবস্থান ব্যবহার করুন।
• স্টেশনগুলো সরলরেখায় বা খুব কাছাকাছি থাকলে: এই স্টেশনগুলো একটি স্থিতিশীল দ্বিমাত্রিক অবস্থান নির্ধারণ করতে পারছে না। বর্তমান রেখা বা ক্লাস্টার থেকে দূরে একটি স্টেশন যোগ করুন।
• গাণিতিক ওভারফ্লো ঘটলে: একটি মান বা মধ্যবর্তী ফলাফল সমর্থিত সংখ্যার সীমা অতিক্রম করেছে।

সরলীকৃত মডেল বনাম বাস্তব সিসমোলজি

এই ক্যালকুলেটরটি শিক্ষামূলক উদ্দেশ্যে ডিজাইন করা হয়েছে এবং এটি পৃথিবীকে ৬,৩৭১ কিমি ধ্রুবক ব্যাসার্ধের একটি নিখুঁত গোলক হিসেবে বিবেচনা করে। এটি ধরে নেয় যে ভূমিকম্পের উৎসটি অগভীর এবং তরঙ্গের গতিবেগ সর্বত্র অভিন্ন।

বাস্তব ক্ষেত্রে, সিসমিক তরঙ্গের গতিপথ পৃথিবীর অভ্যন্তরীণ স্তরগুলোর ঘনত্ব এবং গঠনের পার্থক্যের কারণে বাঁকা হয়ে থাকে। পেশাদার সিসমোলজিস্টরা ভূমিকম্পের কেন্দ্র এবং ফোকাল গভীরতা (focal depth) নিখুঁতভাবে নির্ধারণ করতে ত্রিমাত্রিক গতিবেগ মডেল এবং জটিল ভ্রমণ-সময় সারণী ব্যবহার করেন। এই টুলটি কোনো সিসমিক-ক্যাটালগ অবস্থান বা ভূমিকম্পের পূর্বাভাস প্রদান করে না।

গোপনীয়তা এবং প্রক্রিয়াকরণ

ব্যবহারকারীর গোপনীয়তা রক্ষার্থে এই টুলের সমস্ত হিসাব সরাসরি ব্যবহারকারীর ব্রাউজারে সম্পন্ন হয়। স্টেশনের স্থানাঙ্ক, S−P সময় এবং হিসাবের কোনো ফলাফল কোনো বাহ্যিক সার্ভারে আপলোড করা হয় না।

---

প্রায়শই জিজ্ঞাসিত প্রশ্ন (FAQ)

প্রশ্ন: কেন তিনটি স্টেশনের প্রয়োজন হয়? উত্তর: একটি S−P পরিমাপ কেবল একটি দূরত্বের বৃত্ত দেয়, কোনো দিক নির্দেশ করে না। দুটি বৃত্ত দুটি স্থানে মিলিত হতে পারে। সাধারণত একটি তৃতীয় স্টেশন এই বিকল্পটি সমাধান করে, যেখানে চতুর্থ বা পরবর্তী স্টেশনগুলো অবশিষ্ট মানের মাধ্যমে অমিলগুলো দৃশ্যমান করে তোলে। একটি লাইনে থাকা স্টেশনের চেয়ে উৎসের চারপাশে ছড়িয়ে থাকা স্টেশনগুলো অক্ষাংশ এবং দ্রাঘিমাংশ উভয়ই আরও ভালোভাবে নির্ধারণ করতে পারে।

প্রশ্ন: S−P দূরত্ব কি একদম সঠিক পৃষ্ঠের দূরত্ব? উত্তর: না। ধ্রুব-গতিবেগের সমীকরণটি আগমন-সময়ের পার্থক্য থেকে একটি সাধারণ ভ্রমণ পথ অনুমান করে। সেই ব্যাসার্ধকে গোলাকার পৃষ্ঠের দূরত্ব হিসেবে বিবেচনা করার ক্ষেত্রে একটি অগভীর উৎস এবং একটি সাধারণ গতিবেগ কাঠামো ধরে নেওয়া হয়। প্রকৃত ভ্রমণ-সময়ের বক্ররেখা গভীরতা, দূরত্ব এবং প্রতিটি তরঙ্গ দ্বারা অতিক্রান্ত উপাদানের ওপর নির্ভর করে।

প্রশ্ন: ধনাত্মক বা ঋণাত্মক অবশিষ্ট মান বলতে কী বোঝায়? উত্তর: অবশিষ্ট মান (Residual) হলো ফিট করা গোলাকার দূরত্ব থেকে S−P থেকে প্রাপ্ত দূরত্বের বিয়োগফল। ধনাত্মক মানের অর্থ হলো ফিট করা বিন্দুটি সেই স্টেশন থেকে তার S−P ব্যাসার্ধের চেয়ে দূরে অবস্থিত; ঋণাত্মক মানের অর্থ হলো এটি আরও কাছে অবস্থিত। বড় ও মিশ্র অবশিষ্ট মানগুলো নির্দেশ করে যে বৃত্তগুলোর নির্বাচিত গতিবেগের অধীনে কোনো একটি নিখুঁত সাধারণ ছেদবিন্দু নেই।

প্রশ্ন: অবস্থানের অনিশ্চয়তার মধ্যে কী কী অন্তর্ভুক্ত রয়েছে? উত্তর: এটি একক টাইমিং 1σ মানকে S−P দূরত্বের গুণক এবং ফিট করা বিন্দুর কাছাকাছি স্টেশন জ্যামিতির মাধ্যমে প্রসারিত করে। এর মধ্যে ভুল গতিবেগ মডেল, ফোকাল গভীরতা, ভুলভাবে চিহ্নিত ফেজ, ঘড়ির পক্ষপাত বা পৃথিবীর ত্রিমাত্রিক কাঠামো অন্তর্ভুক্ত নয়, তাই এটি কোনো সম্পূর্ণ নির্ভুলতার দাবি নয়।