Retrieving the Seismic Moment Release Rate of Small Seismic Events Using the Adjoint Method - Theory and Numerical Validation

عندما تهتز الأرض، فإن حتى الهزة الطفيفة تحمل في طياتها قصة خفية حول كيفية تكسر التربة. لقد تمكن علماء الزلازل منذ زمن طويل من تحديد مكان بدء الزلزال وكمية الطاقة التي أطلقها، لكن فهم كيفية إطلاق تلك الطاقة بالضبط بمرور الوقت ظل لغزاً صعباً. تُروى هذه القصة من خلال ما يسمى "دالة زمن المصدر"، وهي سجل لنبض التصدع من أول شرخ له حتى توقفه النهائي. معرفة هذا الإيقاع أمر بالغ الأهمية لأنه يكشف عن سرعة التصدع، واتجاه حركته، وتغيرات الإجهاد التي حدثت في أعماق الأرض. ومع ذلك، فإن استخراج هذا الإيقاع من الموجات الفوضوية المسجلة على السطح أمر صعب للغاية. تعتمد الطرق التقليدية إما على إيجاد زلزال أصغر ومماثل لاستخدامه كمرجع — وهي عملية يمكن أن تكون ذاتية وغير دقيقة — أو تحاول تشغيل عمليات محاكاة ضخمة تستهلك قدرات حاسوبية هائلة وتكون بطيئة جداً للاستخدام الروتيني.

في دراسة جديدة، اقترح الباحثان فويتشيك ديبسكي وكاميل واسكيفيتش من معهد الجيوفيزياء في الأكاديمية البولندية للعلوم طريقة مختلفة للاستماع إلى هذه الإيقاعات الخفية. فقد طورا تقنية تعتمد على الفكرة الفيزيائية البسيطة المتمثلة في أن الصوت والموجات الزلزالية يمكن أن تسافر إلى الوراء في الزمن تماماً كما تسافر إلى الأمام. ومن خلال أخذ التسجيلات من السطح، وقلبها زمنياً، وإرسالها مرة أخرى إلى باطن الأرض باستخدام نموذج حاسوبي، تتجمع الموجات بشكل طبيعي في النقطة التي بدأ منها الزلزال. ومع تقارب هذه الموجات المعكوسة، فإنها تعيد بناء النمط الأصلي لإطلاق الطاقة. اختبر الفريق هذه الفكرة باستخدام محاكاة حاسوبية لزلازل صغيرة ناتجة عن نشاط بشري، مثل تلك الناجمة عن التعدين، باستخدام بيانات اصطناعية خالية من الضجيج. وتظهر نتائجهم أن هذه الطريقة يمكنها بنجاح استعادة التوقيت التفصيلي للتصدع دون الحاجة إلى عمليات عكس رياضية معقدة أو البحث عن زلزال مرجعي. إنها تقدم طريقة أسرع وأكثر قوة لفهم ميكانيكا الأحداث الزلزالية الصغيرة، محولةً مشكلة رياضية صعبة إلى عملية فيزيائية مباشرة.

يعتمد جوهر هذا النهج الجديد على خاصية أساسية في فيزياء الموجات: عكس الزمن. في وسط مرن تماماً، حيث لا تضيع الطاقة على شكل حرارة أو احتكاك، تعمل المعادلات التي تصف حركة الموجات بنفس الطوال سواء كان الزمن يتحرك للأمام أو للخلف. تخيل فيديو لحجر يسقط في بركة؛ إذا قمت بتشغيل الفيديو بالعكس، ستعود التموجات للداخل ويقفز الحجر خارج الماء. وبينما لا يمكنك عكس الزمن فعلياً في العالم الحقيقي، إلا أن الرياضيات تسمح للعلماء بمحاكاة هذه العملية. استخدم الباحثون هذا المفهوم لإنشاء خوارزمية أطلقوا عليها اسم TRSTF. تبدأ العملية بأخذ الموجات الزلزالية المسجلة في محطات مختلفة على الس surfaced. ثم يتم عكس هذه الموجات زمنياً وتغذيتها مرة أخرى في نموذج حاسوبي للقشرة الأرضية. وبينما تتحرك الموجات المعكوسة عبر النموذج، فإنها ترتد عن الحدود وتتفاعل مع الجيولوجيا المعقدة، لتتقارب في النهاية عائدة إلى موقع المصدر الأصلي.

وعندما تلتقي هذه الموجات المعكوسة عند المصدر، فإنها تتداخل مع بعضها البعض. في الموقع الدقيق للزلزال، تصطف الموجات بشكل مثالي، مما يخلق إشارة مركزة وقوية. أما في أي موقع آخر، فتصل الموجات في أوقات غير متزامنة وتلغي بعضها البعض، مما يترك فقط ضجيجاً خافتاً. هذا التأثير التركيزي هو مفتاح الطريقة. فالإشارة التي تظهر عند نقطة التركيز هي إعادة بناء مباشرة لدالة زمن المصدر، أي إيقاع إطلاق الطاقة للزلزال. وجد الباحثون أن هذه التقنية تعمل بشكل جيد للغاية في تجاربهم الاصطناعية الخالية من الضجيج، حتى عندما لا يكون نموذج باطن الأرض معروفاً بدقة. في تجاربهم الحاسوبية، قاموا بمحاكاة ثلاثة سيناريوهات مختلفة لمحاكاة إعدادات المراقبة في العالم الحقيقي: سيناريو حيث وضعت المستشعرات بالقرب من مستوى نشاط التعدين، وسيناريو حيث كانت المستشعرات على السطح فوق الحدث، وثالث بتوزيع مختلف على السطح. وفي كل حالة، نجحت الطريقة في استعادة شكل نمط التصدع الأصلي.

كما استكشفت الدراسة كيف تؤثر أخطاء نموذج سرعة الصوت في الأرض على النتائج. في العالم الحقيقي، لا يعرف الجيولوجيون السرعة الدقيقة للموجات الزلزالية في كل نقطة تحت الأرض؛ بل يعتمدون على التقديرات. اختبر الباحثون طريقتهم من خلال إدخال أخطاء عشوائية وتغيرات كتلية في نماذجهم الحاسوبية، لمحاكاة حالة يكون فيها نموذج السرعة دقيقاً بنسبة 90 بالمائة فقط. وأشارت النتائج إلى أن استخدام نماذج سرعة غير دقيقة يؤدي إلى أخطاء منهجية في دالة زمن المصدر المعاد بناؤها. وقد ثبت أن الطريقة أكثر حساسية للأخطاء المنظمة (الكتلية) منها للأخطاء العشوائية المشتتة. ومع ذلك، رغم هذه الأخطاء، ظلت الإشارة المستعادة قابلة للتمييز ووفرت معلومات مفيدة، مما يثبت أن الطريقة تحتفظ بجدواها حتى مع النماذج غير المثالية، رغم أن دقتها ليست محصنة ضد عدم دقة النماذج.

كانت واحدة من أبرز النتائج هي كيف أثر ترتيب المستشعرات على جودة النتيجة. اكتشف الباحثون أن وجود مستشعرات أقرب إلى حدود المنطقة التي تتم دراستها يساعد في الواقع عملية إعادة البناء. وذلك لأن الانعكاسات من هذه الحدود تعمل كمستشعرات افتراضية إضافية، مما يوفر مزيداً من المعلومات للنظام. في إحدى حالات الاختبار، حيث وضعت المستشعرات بالقرب من الحافة العلوية للنموذج، كانت الإشارة المعاد بناؤها أكثر حدة ودقة من التجهيزات الأخرى، على الرغم من أن البيانات الخام بدت أكثر تشوهاً بسبب التداخل. هذه النتيجة غير المتوقعة تشير إلى أن الطريقة تستفيد من تعقيد مجال الموجة، مستخدمة الانعكاسات لصالحها بدلاً من معاملتها كضجيج.

قام الباحثون أيضاً بقياس خصائص محددة للإشارات المعاد بناؤها، مثل سرعة بدء التصدع ومدته. ووجدوا أنه بينما كانت الإشارات المعاد بناؤها أعرض قليلاً من الإشارات الأصلية، وهو ما يعود على الأرجح إلى العدد المحدود من المستشعرات المستخدمة في المحاكاة، إلا أن الشكل العام والتوقيت قد تم الحفاظ عليهما. إن القدرة على تقدير مدة التصدع أمر حيوي لأنها تساعد العلماء على حساب انخفاض الإجهاد الديناميكي، وهو مقياس لمقدار الإجهاد الذي تم إطلاقه أثناء الحدث. هذه المعلومات ضرورية لفهم ميكانيكا النشاط الزلزالي المستحث، مثل الزلازل الناجمة عن التعدين أو حقن السوائل. أظهرت الدراسة أن الطريقة الجديدة يمكنها تقدير هذه المعايير بدقة تضاهي أو تفوق التقنيات الموجودة، ولكن بجهد حاسوبي أقل بكثير.

على عكس طريقة "قلب الموجة الكامل" (full-waveform inversion)، التي تتطلب تشغيل آلاف عمليات المحاكاة وحلقات تحسين معقدة للوصول إلى أفضل مطابقة، تتطلب هذه الطريقة الجديدة تمريرة واحدة فقط للموجات المعكوسة عبر النموذج. وهذا يجعلها فعالة حاسوبياً، حيث يمكن تشغيلها على أجهزة مكتبية عادية للنماذج الصغيرة أو على أجهزة حاسوب فائقة للنماذج الإقليمية. كما أنها تتجنب الحاجة إلى اختيار زلزال مرجعي، وهي خطوة قد تدخل التحيز البشري وعدم اليقين في الطرق التقليدية. ويؤكد الباحثون أن نهجهم مناسب بشكل خاص للأحداث الزلزالية الصغيرة، حيث يكون المصدر صغيراً بما يكفي ليعامل كنقطة واحدة. وهم يقرون بأن الطريقة تفترض أن الأرض تتصرف كمادة مرنة مثالية، مما يعني أنها لا تمتص الطاقة، وأن المصدر صغير مقارنة بطول موجة الموجات. وهذه افتراضات معقولة للعديد من أحداث التعدين، لكن الباحثين يشيرون إلى أن هناك حاجة لمزيد من الدراسة لفهم كيفية أداء الطريقة في البيئات الأكثر تعقيداً في العالم الحقيقي والتي تشهد فقداً كبيراً للطاقة.

تختتم الورقة البحثية بالتأكيد على إمكانية أن تصبح هذه التقنية أداة قياسية في علم الزلازل. فمن خلال توفير طريقة سريعة وموثوقة وفطرية فيزيائياً لاستعادة دالة زمن المصدر، فإنها تفتح الباب أمام التحليل الروتيني للأحداث الزلزالية الصغيرة التي كان من الصعب سابقاً توصيفها بالتفصيل. ويقترح الباحثون إمكانية دمج هذه الطريقة في أنظمة المراقبة الحالية لتوفير رؤى فورية حول طبيعة النشاط الزلزالي المستحث. وبينما اقتصرت الدراسة الحالية على المحاكاة الحاسوبية، فإن الأساس النظري متين والنتائج واعدة. والخطوة التالية، وفقاً للمؤلفين، ستكون اختبار الطريقة على بيانات زلزالية حقيقية وتدقيق فهم القيود وعدم اليقين المرتبط بها. وفي الوقت الحالي، تقف الدراسة كإثبات لمفهوم أن عكس الزمن يمكن أن يكون عدسة قوية للنظر إلى ماضي الأرض، وكشف القصة الخفية لكيفية تكسر التربة.