Elastic wave propagation in fractured media with spring-type and frictional contact deformation laws

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

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

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

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

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

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

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