فیزیک امواج لرزهای و اختلاف زمان رسیدن S−P
مبنای فیزیکی تعیین موقعیت مرکز سطحی زلزله، تفاوت در سرعت انتشار امواج لرزهای مختلف درون زمین است. هنگام وقوع زلزله، امواج متراکمکننده یا اولیه (P-waves) و امواج برشی یا ثانویه (S-waves) به طور همزمان از کانون زلزله آزاد میشوند. از آنجا که امواج P با سرعت بیشتری نسبت به امواج S حرکت میکنند، فاصله زمانی بین ثبت این دو موج در ایستگاههای لرزهنگاری با افزایش فاصله از منبع بزرگتر میشود. این اختلاف زمان ثبت که با عنوان بازه S−P شناخته میشود، امکان محاسبه فاصله مستقیم ایستگاه تا منبع لرزهای را فراهم میکند.
در یک مدل پوسته یکنواخت با سرعتهای ثابت، رابطه بین زمان سیر و فاصله به صورت ریاضی فرمولبندی میشود. با فرض اینکه سرعت موج P برابر با Vₚ و سرعت موج S برابر با Vₛ باشد، اختلاف زمان رسیدن (Δ t) با فاصله مستقیم (D) از طریق رابطه Δ t = D(1/Vs - 1/Vp) مرتبط است. با بازآرایی این فرمول، تقریب فاصله به صورت زیر به دست میآید:
D = Δ t ÷ (1/Vs - 1/Vp)
این ابزار آموزشی با فرض یک منبع کمعمق و سرعتهای موج ثابت، این فاصله مستقیم را به عنوان فاصله روی سطح زمین در نظر میگیرد.
هندسه مثلثسازی و بهینهسازی کمترین مربعات
تعیین موقعیت دقیق مرکز سطحی زلزله به توزیع هندسی ایستگاههای ثبتکننده بستگی دارد. ثبت بازه S−P در یک ایستگاه منفرد، تنها یک دایره به شعاع مشخص به دور آن ایستگاه ایجاد میکند و جهت منبع را مشخص نمیسازد. اضافه شدن ایستگاه دوم، نقاط تقاطع را به دو نقطه محدود میکند. برای حل این ابهام و یافتن یک موقعیت منحصربهفرد، به حداقل سه ایستگاه نیاز است.
در دنیای واقعی، به دلیل خطاهای اندازهگیری و ناهمگنیهای زمین، دایرههای ترسیمشده از ایستگاههای مختلف دقیقاً در یک نقطه متقاطع نمیشوند. برای حل این مشکل، از روش بهینهسازی کمترین مربعات استفاده میشود تا بهترین نقطه انطباق ریاضی پیدا شود. این روش با کمینهسازی مجموع مجذور تفاوتهای بین فاصلههای محاسباتی و فاصلههای رصدشده کار میکند:
برازش مرکز سطحی: کمینهسازی Σ(dᵢ − Dᵢ)² در تمام ایستگاهها
در این رابطه، Dᵢ فاصله استخراجشده از بازه S−P برای ایستگاه i و dᵢ فاصله سطحی روی مدل کروی زمین است که با استفاده از فرمول هاورسین محاسبه میشود:
فاصله سطحی: d = R × زاویه مرکزی هاورسین، با R = 6371 km
راهنمای وارد کردن دادهها و پارامترهای ورودی
برای استفاده از محاسبهگر، کاربر باید اطلاعات مربوط به حداقل ۳ و حداکثر ۸ ایستگاه لرزهنگاری را وارد کند. مقادیر ورودی و محدودیتهای مجاز آنها به شرح زیر است:
* **عرض جغرافیایی**: باید در محدوده -90^° تا 90^° باشد.
* **طول جغرافیایی**: باید در محدوده -180^° تا 180^° باشد.
* **بازه S−P**: باید بزرگتر از 0 و حداکثر 2000 ثانیه باشد.
* **سرعت موج P**: باید بزرگتر از 0 و حداکثر 20 km/s بوده و از سرعت موج S بیشتر باشد.
* **سرعت موج S**: باید بزرگتر از 0 و حداکثر 20 km/s بوده و از سرعت موج P کمتر باشد.
* **عدم قطعیت زمانی (1σ)**: باید بزرگتر از 0 و حداکثر 60 ثانیه باشد.
کاربران میتوانند برای بررسی یک سناریوی واقعی، روی دکمه بارگذاری نمونه توکیو کلیک کنند تا دادههای نمونه در جدول قرار گیرند، یا با کلیک روی دکمه پاک کردن تمامی فیلدها را بازنشانی نمایند.
تحلیل خروجیها و معیارهای کیفیت برازش
پس از وارد کردن دادهها، ابزار نتایج محاسباتی را در قالب بخشهای مختلف نمایش میدهد:
بخش پاسخ مرکز سطحی
این بخش شامل مختصات جغرافیایی نقطه بهینه تحت عنوان مرکز سطحی تخمینی است. کیفیت این محاسبات با معیارهای زیر سنجیده میشود:
- با کیفیت برازش فاصله: به صورت یکی از سه حالت برازش فاصله نزدیک، برازش فاصله مخلوط یا برازش فاصله ضعیف گزارش میشود.
- کیفیت هندسه ایستگاه: به صورت هندسه ایستگاه قوی، هندسه ایستگاه متوسط یا هندسه ایستگاه ضعیف نمایش داده میشود که نشاندهنده میزان پراکندگی مناسب ایستگاهها در اطراف مرکز سطحی است.
- باقیمانده فاصله RMS و بزرگترین باقیمانده فاصله.
- موقعیت تقریبی 1σ و عدم قطعیت فاصله 1σ به ازای هر ایستگاه.
طرح هندسه ایستگاهها
یک تصویر شماتیک از موقعیت نسبی ایستگاهها و مرکز سطحی تخمینزدهشده نمایش داده میشود که زیر آن عبارت زیر درج شده است:
فقط هندسه نسبی — بدون نقشه، مرزها یا دادههای زنده زلزله.
جدول برازش ایستگاه
این جدول شامل ستونهای ایستگاه، فاصله S−P، فاصله کروی و باقیمانده برای تکتک ایستگاهها است. باقیمانده از تفاضل فاصله کروی برازششده و فاصله استخراجشده از زمان S−P به دست میآید.
قوانین، خطاهای ورودی و محدودیتهای مدل
این ابزار مجهز به سیستم اعتبارسنجی دقیق ورودیها است و در صورت بروز خطا، پیامهای زیر را نمایش میدهد:
- تعداد نامناسب ایستگاهها:
بین سه تا هشت ایستگاه استفاده کنید. - ترتیب نادرست سرعتها:
سرعت موج P باید بیشتر از سرعت موج S باشد. - خارج از محدوده بودن عدم قطعیت زمانی:
عدم قطعیت زمانی باید بزرگتر از 0 و حداکثر 60 ثانیه باشد. - هندسه ناپایدار (مانند قرارگیری ایستگاهها روی یک خط مستقیم یا تجمع در یک نقطه):
این ایستگاهها موقعیت دو بعدی پایداری را محدود نمیکنند. ایستگاهی دور از خط یا خوشه فعلی اضافه کنید. - سرریز محاسباتی:
یک مقدار یا نتیجه میانی از محدوده عددی پشتیبانیشده فراتر رفته است.
محدودیتهای علمی مدل
این ابزار یک مدل آموزشی سادهشده است. زمین به عنوان یک کره کامل با شعاع ثابت 6371 km فرض شده است. این مدل از عمق کانون زلزله، مسیرهای پرتو منحنی درون زمین، ساختارهای سرعتی لایهای و سهبعدی پوسته، اشتباهات در انتخاب فاز و انحراف زمانی ساعت ایستگاهها چشمپوشی میکند. برای کارهای حرفهای و ثبت رسمی در کاتالوگهای لرزهنگاری، باید از مدلهای زمان سیر منطقهای استفاده شود.
همچنین، تمام فرآیندهای محاسباتی به طور کامل درون مرورگر کاربر انجام میشود و هیچ دادهای به سرورهای خارجی ارسال نمیگردد.
پرسشهای متداول (FAQ)
چرا به سه ایستگاه نیاز است؟
یک اندازهگیری S−P یک دایره فاصله را مشخص میکند، نه یک جهت را. دو دایره میتوانند در دو نقطه با هم تلاقی داشته باشند. یک ایستگاه سوم معمولاً این انتخاب را حل میکند، در حالی که ایستگاه چهارم یا ایستگاههای بعدی، عدم تطابق را از طریق باقیماندهها نمایان میسازند. ایستگاههایی که در اطراف منبع پراکنده شدهاند، عرض و طول جغرافیایی را بهتر از ایستگاههای واقع در یک خط مستقیم محدود میکنند.
آیا فاصله S−P همان فاصله دقیق روی سطح زمین است؟
خیر. معادله سرعت ثابت، یک مسیر حرکت مشترک را بر اساس اختلاف زمان رسیدن تخمین میزند. در نظر گرفتن این شعاع به عنوان فاصله روی سطح کروی، فرض را بر یک منبع کمعمق و ساختار سرعت ساده میگذارد. منحنیهای زمان سیر واقعی به عمق، فاصله و موادی که هر پرتو از آنها عبور میکند بستگی دارند.
باقیمانده مثبت یا منفی چه معنایی دارد؟
باقیمانده عبارت است از فاصله کروی برازششده منهای فاصله استخراجشده از S−P. مقدار مثبت به این معنی است که نقطه برازششده نسبت به شعاع S−P آن، از آن ایستگاه دورتر است؛ مقدار منفی یعنی نزدیکتر است. باقیماندههای بزرگ و مخلوط نشان میدهند که دایرهها تحت سرعتهای انتخابشده، تقاطع تمیز و واحدی ندارند.
چه مواردی در عدم قطعیت موقعیت گنجانده شده است؟
این مقدار، تک عدم قطعیت زمانی 1σ را از طریق ضریب فاصله S−P و هندسه ایستگاه در نزدیکی نقطه برازششده منتقل میکند. این محاسبه شامل مدل سرعت نادرست، عمق کانون، فازهای اشتباه برچسبگذاریشده، انحراف ساعت یا ساختار سه بعدی زمین نمیشود، بنابراین یک ادعای دقت کامل نیست.
منابع علمی و مراجع
- سازمان زمینشناسی ایالات متحده (USGS)، «علم زلزله»: USGS: علم زلزله
- مواد آموزشی USGS، «امواج لرزهای و موقعیتیابی زلزله»: USGS: امواج لرزهای و موقعیتیابی زلزله
- سازمان زمینشناسی ایالات متحده (USGS)، «زمانهای سیر زلزله»: USGS: زمانهای سیر زلزله