A Finite Volume Method for Elastic Waves in Heterogeneous, Anisotropic and Fractured Media

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

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

ركز الباحثون على نهج عددي محدد يسمى "طريقة الحجم المحدود" (finite volume method). تخيل النطاق الموجود تحت الأرض ليس ككتلة ملساء ومستمرة، بل كمجموعة من القطع الصغيرة غير منتظمة الشكل، مثل فسيفساء مكونة من أحجار ذات أحجام وأشكال مختلفة. تتيح هذه المرونة للنموذج أن يتناسب تماماً حول صدع متعرج أو طبقة صخرية معقدة دون إجبار البيانات على الدخول في شبكة صلبة. وضمن هذا الإطار، طبقوا تقنية تُعرف باسم "تقريب الإجهاد متعدد النقاط مع التماثل الضعيف" (Multi-Point Stress Approximation with weak symmetry). وببساطة، تقوم هذه الطريقة بحساب كيفية ضغط الصخور وسحبها لجيرانها من خلال النظر في القوى المؤثرة على قطع متعددة محيطة في آن واحد، بدلاً من الاكتفاء بالقطع المجاورة مباشرة فقط. وهذا يسمح للنموذটি بالتعامل مع الصخور "تباينية الخواص" (anisotropic)، مما يعني أن صلابتها تتغير اعتماداً على الاتجاه الذي تنظر إليه، تماماً كما هو حال الخشب الذي يكون أقوى على طول الألياف منه عبرها.

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

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

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

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