A Nonmatching-Grid Discontinuous Galerkin Method for Nonlinear Reactive Transport in Hybrid-Dimensional Fractured Media

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

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

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

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

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

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

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

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