فيزياء الموجات الزلزالية وفواصل الوصول الزمنية
تنتقل الطاقة الناتجة عن كسر الصدع الزلزالي عبر الأرض في شكل موجات زلزالية مختلفة. وتعد الموجات الأولية (P-waves) والموجات الثانوية (S-waves) الأطوار الأساسية المستخدمة لتحديد موقع مصدر الزلزال. تتميز الموجات P بأنها موجات انضغاطية سريعة تنتقل كأول وصول يتم تسجيله على أجهزة قياس الزلازل، تليها الموجات S التي تمثل موجات قص مستعرضة وتتحرك بسرعة أبطأ.
بسبب هذا التفاوت في السرعة، يتسع الفارق الزمني بين وصول الموجتين P و S كلما زادت المسافة التي تقطعها الموجات من مركز الزلزال إلى محطة الرصد. ويُعرف هذا الفارق باسم "فاصل S−P الزمني". ومن خلال معرفة سرعة كل من الموجتين P و S في الوسط الصخري، يمكن استخدام هذا الفارق الزمني لحساب المسافة المباشرة من المحطة إلى مصدر الزلزال.
تحت افتراض ثبات السرعة في نموذج القشرة الموحدة، ترتبط المسافة D بالفاصل الزمني Δ t عبر العلاقة الفيزيائية التالية: Δ t = D × (1 / Vₛ - 1 / Vₚ)
وعند إعادة ترتيب هذه المعادلة لعزل المسافة، نحصل على صيغة تقريب المسافة المستخدمة في الحساب: D = Δ t ÷ (1/Vs - 1/Vp)
حيث يمثل Vₚ سرعة الموجة P، ويمثل Vₛ سرعة الموجة S.
هندسة التثليث على نموذج الأرض الدائرية
عند تسجيل فاصل S−P زمني في محطة رصد واحدة، فإن ذلك يحدد مسافة نصف قطرية دائرية تحيط بالمحطة، مما يعني أن المركز السطحي للزلزال يمكن أن يقع في أي نقطة على محيط هذه الدائرة.
- محطة واحدة: توفر دائرة مسافة واحدة دون تحديد الاتجاه.
- محطتان: تتقاطع دائرتا المسافة في نقطتين محتملتين كحد أقصى.
- ثلاث محطات أو أكثر: تتقاطع الدوائر الثلاث في نقطة فريدة تمثل المركز السطحي للزلزال، وهو ما يُعرف بعملية التثليث الزلزالي.
تستخدم هذه الحاسبة نموذجًا كرويًا للأرض بنصف قطر ثابت يبلغ R = 6371 km. ولحساب المسافات على هذا السطح المنحني، تعتمد الأداة على صيغة الهافرسين لحساب الزاوية المركزية بين إحداثيات المحطات والمركز السطحي المقدر: d = R × زاوية الهافرسين المركزية، مع R = 6371 km
تتطلب الأداة إدخال بيانات ما بين 3 محطات كحد أدنى و 8 محطات كحد أقصى لضمان استقرار الحل الهندسي وتجنب الأخطاء الناتجة عن قلة البيانات.
معالجة المربعات الصغرى والملاءمة الرياضية
في التطبيقات العملية، نادراً ما تتقاطع دوائر المسافات المحسوبة من المحطات المختلفة في نقطة واحدة بدقة متناهية، ويعود ذلك إلى وجود أخطاء في قراءة أوقات الوصول أو عدم تجانس سرعات الموجات في باطن الأرض. لتجاوز هذه المشكلة، تطبق الحاسبة منهجية معالجة المربعات الصغرى (Least-Squares Treatment) للوصول إلى أفضل ملاءمة هندسية.
تقوم هذه الطريقة بتقليل مجموع مربعات الفوارق بين المسافة السطحية الفعلية (dᵢ) والمسافة المستنتجة من فاصل S−P الزمني (Dᵢ) عبر جميع المحطات المستخدمة: ملاءمة المركز السطحي: تقليل Σ(dᵢ - Dᵢ)² عبر جميع المحطات
ويُعرف الفرق بين المسافة الدائرية الملاءمة والمسافة المستنتجة من زمن الوصول بالفارق المتبقي (Residual):
- الفارق المتبقي الموجب: يعني أن المركز السطحي المقدر يقع على مسافة أبعد من المحطة مقارنة بنصف قطر S−P المحسوب لها.
- الفارق المتبقي السالب: يعني أن المركز السطحي المقدر يقع في نقطة أقرب إلى المحطة مقارنة بنصف قطر S−P المحسوب لها.
مدخلات ومحددات الأداة
لإجراء عملية التثليث وحساب المركز السطحي، يتعين على المستخدم إدخال البيانات التالية بدقة ضمن النطاقات المحددة:
أرصاد المحطة (بين 3 و 8 محطات)
- اسم المحطة: معرف نصي لكل محطة رصد.
- خط العرض: يجب أن تكون القيمة بين -90^° و 90^°. وإذا تم إدخال قيمة خارج هذا النطاق، يظهر الخطأ:
المحطة 1: يجب أن يكون خط العرض من −90° إلى 90°. - خط الطول: يجب أن تكون القيمة بين -180^° و 180^°. وإذا تم إدخال قيمة خارج هذا النطاق، يظهر الخطأ:
المحطة 1: يجب أن يكون خط الطول من −180° إلى 180°. - فاصل S−P الزمني: يجب أن يكون أكبر من 0 ولا يزيد عن 2000 ثانية. وفي حال المخالفة يظهر الخطأ:
المحطة 1: أدخل فاصل S−P زمني أكبر من 0 ولا يزيد عن 2000 ثانية.
نموذج الانتشار
- سرعة الموجة P: يجب أن تكون أكبر من 0 ولا تزيد عن 20 km/s، ويجب أن تكون أكبر من سرعة الموجة S.
- سرعة الموجة S: يجب أن تكون أكبر من 0 ولا تزيد عن 20 km/s، ويجب أن تكون أقل من سرعة الموجة P.
- عدم اليقين في التوقيت (1σ): يجب أن يكون أكبر من 0 ولا يزيد عن 60 ثانية.
يمكن للمستخدم النقر على زر تحميل مثال طوكيو لتعبئة الحقول تلقائيًا ببيانات تجريبية متناسقة، أو النقر على زر مسح لإعادة تعيين كافة الحقول والبدء من جديد.
مخرجات الحساب وتحليل النتائج
تعرض الحاسبة النتائج مقسمة إلى عدة أقسام تفصيلية لتسهيل قراءتها وتحليلها:
1. حل المركز السطحي
- المركز السطحي المقدر: يعرض إحداثيات خطوط الطول والعرض الجغرافية للمركز السطحي المحسوب.
- ملاءمة مسافة متقاربة أو ملاءمة مسافة مختلطة أو ملاءمة مسافة ضعيفة: لتقييم مدى تطابق دوائر المسافات عند الحل المشترك.
- هندسة محطات قوية أو هندسة محطات متوسطة أو هندسة محطات ضعيفة: لتقييم التوزيع الجغرافي للمحطات حول المركز السطحي.
- جذر متوسط مربعات (RMS) الفوارق المتبقية للمسافة و أكبر فارق متبقي للمسافة.
- عدم اليقين التقريبي في الموقع 1σ و عدم اليقين في المسافة 1σ لكل محطة.
2. مخطط هندسة المحطات
رسم تخطيطي يوضح التوزيع الهندسي النسبي للمحطات والمركز السطحي المقدر، ويحمل التوضيح التالي: "الهندسة النسبية فقط — لا توجد خريطة أو حدود أو بيانات زلازل حية."
3. ملاءمة المحطة (جدول النتائج)
جدول تفصيلي يحتوي على الأعمدة التالية لكل محطة:
- المحطة
- مسافة S−P (المسافة المستنتجة من الفارق الزمني)
- المسافة الدائرية (المسافة الفعلية على النموذج الكروي من المحطة إلى المركز السطحي المقدر)
- الفارق المتبقي
4. المعادلة والتعويض
يعرض هذا القسم الخطوات الرياضية المتبعة في الحساب بالتفصيل:
- صيغة تقريب المسافة وصيغة المسافة السطحية الكروية.
- حساب عامل المسافة:
عامل المسافة = vp × vs ÷ (vp − vs) = factor كم لكل ثانية من S−P - اشتقاق مسافة كل محطة بشكل منفرد:
station: interval ثانية × factor كم/ثانية = distance كم - إحداثيات الملاءمة النهائية:
أفضل ملاءمة دائرية = latitude°، longitude°؛ جذر متوسط مربعات الفوارق المتبقية (RMS) = rms كم - حساب انتشار الخطأ:
قيمة عدم اليقين في التوقيت 1σ البالغة timing ثانية تعطي عدم يقين في المسافة 1σ بمقدار distance كم؛ وتنقله هندسة المحطات إلى حوالي position كم على المحور الأفقي الرئيسي.
قواعد التحقق من البيانات ومعالجة الأخطاء
تطبق الحاسبة قواعد صارمة للتحقق من صحة المدخلات لضمان استقرار العمليات الحسابية وتجنب الانهيار البرمجي:
- عدد المحطات: إذا حاول المستخدم الحساب باستخدام أقل من 3 محطات أو أكثر من 8 محطات، يظهر تنبيه:
يرجى استخدام ما بين ثلاث وثماني محطات. - القيم غير الرقمية: عند إدخال نص أو رمز غير رياضي في حقول الأرقام، يظهر الخطأ:
where: «token» ليس رقمًا. - ترتيب السرعات: إذا تم إدخال سرعة موجة S مساوية أو أكبر من سرعة موجة P، يظهر الخطأ:
يجب أن تكون سرعة الموجة P أكبر من سرعة الموجة S. - نطاق السرعات: إذا خرجت قيم السرعة عن الحدود المسموحة، يظهر الخطأ:
يجب أن يكون field أكبر من 0 ولا يزيد عن 20 km/s. - عدم اليقين في التوقيت: إذا كانت القيمة خارج النطاق المسموح، يظهر الخطأ:
يجب أن يكون عدم اليقين في التوقيت أكبر من 0 ولا يزيد عن 60 ثانية. - تطابق الإحداثيات: لمنع تداخل الدوائر الحسابية بشكل متطابق، يمنع استخدام نفس الإحداثيات لمحطتين مختلفتين، ويظهر الخطأ:
المحطتان station وother تملكان الإحداثيات نفسها. يرجى استخدام مواقع محطات متميزة. - الهندسة غير المستقرة: إذا كانت المحطات موزعة على خط مستقيم واحد أو متجمعة في نقطة ضيقة للغاية، فلن تتمكن الخوارزمية من تحديد موقع ثنائي الأبعاد مستقر، ويظهر الخطأ:
هذه المحطات لا تحدد موقعًا ثنائي الأبعاد مستقرًا. أضف محطة بعيدة عن الخط أو التجمع الحالي. - الفيض العددي: في حال تجاوزت الحسابات الوسيطة الحدود الرياضية للحاسوب، يظهر الخطأ:
تجاوزت قيمة أو نتيجة وسيطة نطاق الأرقام المدعوم.
حدود النموذج التعليمي مقارنة بالواقع
تعد هذه الحاسبة أداة تعليمية مبسطة تعتمد على افتراضات مثالية لتسهيل فهم المبادئ الأساسية لعلم الزلازل. وتختلف هذه الافتراضات عن الطرق الاحترافية المتبعة في المراصد الزلزالية لإنتاج الفهارس الرسمية للزلازل:
| ميزة النموذج التعليمي في الحاسبة | الواقع في المراصد والكتالوجات الاحترافية |
|---|---|
| يفترض أن مصدر الزلزال ضحل للغاية (على السطح مباشرة). | يأخذ في الاعتبار العمق البؤري (Focal Depth) للزلزال في الفضاء ثلاثي الأبعاد. |
| يعتمد على سرعات موجات ثابتة وموحدة في جميع الاتجاهات. | يستخدم نماذج سرعة طبقية وثلاثية الأبعاد تعكس تغير كثافة طبقات الأرض. |
| يفترض انتقال الموجات في مسارات مستقيمة على طول السطح الكروي. | يحسب مسارات الأشعة المنحنية الناتجة عن انكسار وانعكاس الموجات في باطن الأرض. |
| يقتصر حساب عدم اليقين على أخطاء التوقيت المدخلة فقط. | يشمل مصادر أخطاء متعددة مثل انحياز ساعات المحطات، وأخطاء تحديد الأطوار، والتباين الإقليمي للسرعات. |
خصوصية البيانات ومعالجتها
تتميز هذه الأداة بالخصوصية التامة والأمان المحلي؛ حيث تتم جميع العمليات الحسابية ومعالجة البيانات مباشرة داخل متصفح الويب الخاص بالمستخدم وعلى جهازه المحلي. لا يتم رفع إحداثيات المحطات، أو أوقات S−P، أو نتائج تحديد الموقع إلى أي خوادم خارجية على الإطلاق، مما يضمن بقاء بياناتك محليًا طوال فترة الاستخدام.
الأسئلة الشائعة
لماذا نحتاج إلى ثلاث محطات؟
يعطي قياس S−P واحد دائرة مسافة، وليس اتجاهًا. ويمكن لدائرتين أن تلتقيا في موقعين. وعادةً ما تحسم المحطة الثالثة هذا الخيار، بينما تجعل المحطة الرابعة أو ما يليها عدم التطابق مرئيًا من خلال الفوارق المتبقية. وتعمل المحطات الموزعة حول المصدر على تقييد خطوط الطول والعرض بشكل أفضل من المحطات الواقعة على خط واحد.
هل مسافة S−P هي المسافة السطحية الدقيقة؟
لا. تقدّر معادلة السرعة الثابتة مسار انتقال مشترك من فرق وقت الوصول. وافتراض أن نصف القطر ذلك يمثل مسافة سطحية كروية يفترض وجود مصدر ضحل وبنية سرعة بسيطة. تعتمد منحنيات وقت انتقال الزلازل الحقيقية على العمق والمسافة والمادة التي تعبرها كل موجة.
ماذا يعني الفارق المتبقي الموجب أو السالب؟
الفارق المتبقي هو المسافة الدائرية الملاءمة مطروحًا منها المسافة المستنتجة من S−P. القيمة الموجبة تعني أن النقطة الملاءمة أبعد عن تلك المحطة من نصف قطر S−P الخاص بها؛ والقيمة السالبة تعني أنها أقرب. تشير الفوارق المتبقية المختلطة الكبيرة إلى أن الدوائر لا تشترك في تقاطع واحد واضح تحت السرعات المحددة.
ما الذي يشتمل عليه عدم اليقين في تحديد الموقع؟
يقوم بحساب تأثير قيمة عدم اليقين الفردية في التوقيت 1σ عبر عامل مسافة S−P وهندسة المحطات بالقرب من النقطة الملاءمة. ولا يشمل ذلك نموذج السرعة غير الصحيح، أو عمق البؤرة، أو الأطوار الخاطئة، أو انحياز ساعة المحطة، أو بنية الأرض ثلاثية الأبعاد، لذا فهو لا يمثل ادعاءً بالدقة الكاملة.