🔬 مادة مكثفة

A Conservative Darcy-Brinkman-Stokes Formulation for PHREEQC-Coupled Reactive Transport with Cell-Based Surface Complexation

تقدم هذه الورقة صياغة "دارسي-برينكمان-ستوكس" حافظة للكتلة تم تنفيذها في برنامج GeoChemFoam وربطها ببرنامج PHREEQC، مما يتيح عمليات محاكاة دقيقة للنقل التفاعلي مع تعقيد سطحي قائم على الخلايا في كل من البنى المسامية المحلولة وغير المحلولة.

المؤلفون الأصليون: Hannah P. Menke, Florian Doster, Julien Maes

نُشر 2026-10-15

📖 4 دقيقة قراءة☕ قراءة في استراحة قهوة

CC BY 4.0

المؤلفون الأصليون: Hannah P. Menke, Florian Doster, Julien Maes

البحث الأصلي مرخَّص بموجب CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ هذا شرح مولَّده بالذكاء الاصطناعي للبحث أدناه. لم يكتبه المؤلفون ولم يصادقوا عليه. وللتحقق من الدقة التقنية، يرجى الرجوع إلى البحث الأصلي. اقرأ إخلاء المسؤولية الكامل

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

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

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

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

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

ملخص تقني: صياغة "دارسي-برينكمان-ستوكس" محافظة لنمذجة النقل التفاعلي المقترن ببرنامج PHREEQC

بيان المشكلة
يمثل النقل التفاعلي في الأوساط المسامية متعددة المقاييس، مثل الكربونات، تحدياً كبيراً في النمذجة. تحتوي هذه الأنظمة غالباً على مسام كبيرة مفتوحة إلى جانب مناطق مسامية دقيقة (microporous). وبينما يمكن للمحاكاة القائمة على الصور حل المسام الكبيرة، فإن مساحة السطح التفاعلية الداخلية للأطوار المسامية الدقيقة غالباً ما تكون أقل من دقة الصورة (على مقياس دون الخلية - sub-voxel). تعتمد نماذج "دارسي" الكلاسيكية على افتراضات رفع المقييس البنيوية التي لا تستطيع التقاط عدم التجانس المكاني، بينما تفشل نماذج "دارسي-برينكمان-ستوكس" (DBS) كاملة الدقة (مثل تلك التي تعتمد على الحدود الصريحة) عندما تقع المساحة التفاعلية المهيمنة داخل الخلايا الحسابية الخشنة بدلاً من حدود الشبكة. تعمل سير عمل PHREEQC الحالية عادةً على فرض الكيمياء كشروط تدفق على الواجهات المحلولة، مما يترك مساحات السطح الداخلية غير المحلولة إما متجاهلة أو ممثلة كمصارف حجمية عشوائية (ad hoc) تخل بحفظ الكتلة والتعريف الفيزيائي لمخزونات مقياس الخلية.

المنهجية
يقدم المؤلفون صياغة DBS محافظة على الكتلة تم تنفيذها في GeoDecoam (وهي مجموعة أدوات قائمة على OpenFOAM) ومقترنة بـ PHREEQC عبر واجهة PhreeqcRM. الابتكار الجوهري هو نهج تعقيد السطح القائم على الخلية حيث تُعامل المساحة السطحية التفاعلية كخاصية لحجم الخلية بدلاً من كونها حداً للشبكة.

المعادلات الحاكمة: تخضع عملية التدفق لمعادلات دارسي-برينكمان-ستوكس، حيث تتغير حقل المسامية ( $\epsilon$ ) الممركز في الخلية بسلاسة من الحالة الصلبة إلى المسام المفتوحة. يتم حل النقل لكسر الماء المتحرك الفيزيائي ( $\theta$ )، الذي يحكم التخزين المائي والانتشار.
تعريف المساحة السطحية: تُعرف المساحة السطحية التفاعلية ( $a_{surf}$ $a_{s u r f}$ ) ككثافة حجمية مركزة في الخلية ( $m^2 m^{-3}_{bulk}$ $m^{2} m_{b u l k}^{- 3}$ ) تتكون من حدين مضافين:
1. حد الواجهة ( $a_{interface}^{surf}$ ): يُعاد بناؤه من واجهة الصورة المحلولة بين طور المسام والطور الصلب باستخدام خوارزميات الواجهة الناعمة (مثل occupancyGaussian).
2. الحد الجوهري ( $a_{intrinsic}^{surf}$ ): مساهمة محددة مسبقاً لما دون مستوى الخلية (sub-voxel) مشتقة من كثافات الأطوار النوعية ( $s_k$ ) والموزونة بالكسر الهندسي للطور ( $f_k$ ) داخل الخلية. يسمح هذا بتمثيل المساحة السطحية المسامية الدقيقة الداخلية التي لا يمكن حلها بالصورة.
الاقتران المحافظ: يُعرف مخزون الكتلة المحلي كـ $M = \theta C + a_{surf} \Gamma$ $M = θ C + a_{s u r f} Γ$ ، حيث $C$ $C$ هو التركيز المائي، و $\Gamma$ $Γ$ هو التركيز السطحي. يستخدم الحل مخطط تقسيم "لي" (Lie splitting) من الدرجة الأولى:
1. خطوة النقل: تقوم بتقديم الأنواع المائية بشكل محافظ باستخدام نقل DBS ذي الحجم المحدود.
2. خطوة الكيمياء: تقوم برسم حالة الماء المحدثة، وكسر الماء المتحرك، والمخزون السطحي الممركز في الخلية إلى PHREEQC. يقوم PHREEQC بحل التوازن الديناميكي الحراري المحلي (تعقيد السطح) ويعيد الأنواع المائية والممتزة المحدثة، مما يضمن حفظ الكتلة بدقة بين الأطوار المائية والسطحية داخل كل خلية.
التنفيذ: تُعامل المساحة السطحية التفاعلية وكثافات المواقع كحقول ثابتة خلال المحاكاة. ويتم استبعاد الخلايا التي تقل فيها كسور الماء المتحرك عن عتبة ( $\theta_{mask}$ ) من العمليات الكيميائية.

المساهمات الرئيسية

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

النتائج والتحقق
تم التحقق من الصيغة من خلال مجموعة اختبار شاملة:

اختبار التأخير أحادي البعد: طابقت استجابة النقل-الانتشار-الامتزاز مع الحلول التحليلية والمراجع ذات النطاق المحدود، مما أظهر عوامل تأخير صحيحة وتقارباً من الدرجة الأولى في الزمان والمكان.
توازن النظام المغلق: أكدت اختبارات توازن الكتلة أن انتقال الكتلة بين الأطوار المائية والسطحية داخل PHREEQC يحافظ على المخزون الإجمالي مع بقاء البواقي أقل من $10^{-9}$ .
إعادة بناء المساحة السطحية: اختُبرت خوارزميات الواجهة الناعمة مقابل هندسات تحليلية (كرات وأعمدة دائرية). كانت الأخطاء محدودة ضمن مستويات أقل من 1% (على سبيل المثال، $<0.2\%$ للأعمدة ثنائية الأبعاد، و $<1\%$ للكرات ثلاثية الأبعاد عند الدقة الكافية) مقارنة بمرجع snappyHexMesh المتوافق مع الجسم.
مقارنة المقياس المجهري بـ DBS: أظهرت مقارنة نموذج "ليثيوم ميكروموديل" بين محاكاة كاملة الدقة لفتحات المسام وصيغة DBS المرفوعة المقياس، أن تمثيل المساحة الجوهرية الممركز في الخلية أعاد إنتاج كل من المخزونات المائية والمرتبطة بدقة تصل إلى 0.6% عبر أرقام بيكليه ( $Pe$ ) تتراوح من 0.1 إلى 10. وبينما كان إجمالي الامتصاص دقيقاً للغاية، أظهر ملف التركيز الطولي الممتز اختلافات في الشكل تصل إلى 4.75%، وهو ما عُزي إلى تنعيم تغيرات التركيز على مقياس المسام.
حالات تطبيقية: طُبق النموذج على ثلاثة سيناريوهات معقدة قائمة على الصور:
- شبكة كسور ثنائية الأبعاد: تأخير الرصاص ( $Pb^{2+}$ ) بواسطة طلاءات هيدروسيد الحديديك (HFO)، مما رصد عدم التجانس المكاني في إشغال المواقع القوية.
- خلية كيلفن ثلاثية الأبعاد: تبادل تنافسي بين المنجنيز ( $Mn^{2+}$ ) والمغنيسيوم ( $Mg^{2+}$ )، مما أظهر قدرة النموذج على التعامل مع التنافس متعدد الأنواع في شبكة مسامية دقيقة.
1. حجر جيري "إستاياليد" ثلاثي الأبعاد: تبادل السترونشيوم ( $Sr^{2+}$ ) والكالسيوم ( $Ca^{2+}$ ) لمحاكاة انتقال المواد المشعة، مما أظهر أنماط تبادل كروماتوغرافية حيث يحل السترونشيوم محل الكالسيوم في كل من الطور المائي والطور الممتز.

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

أبحاث مشابهة