Voronoi-entropy characterisation of geometric disorder in two-dimensional porous networks
تُثبت هذه الورقة أن إنتروبيا شانون لتوزيع مساحة خلايا فورونوي تعمل كمقياس عددي مستقل عن المسامية لتكميم الاضطراب الهندسي في الشبكات المسامية ثنائية الأبعاد، مما يكمل مقاييس المسامية والنفاذية التقليدية من خلال رصد كيفية تأثير إعادة تنظيم الفراغات المسامية على انتقال السوائل.
المؤلفون الأصليون:Vladivostok Suxo, Alexsandro Kirch, Caetano Rodrigues Miranda
يجب على السوائل التي تتحرك عبر الأرض — سواء كانت مياهاً تتسرب إلى طبقة مياه جوفية، أو نفطاً ينساب عبر الصخور، أو حبراً يتشرب في الورق — أن تجتاز متاهة خفية وملتوية. وتعتمد سهولة مرور هذه السوائل على أمرين: مقدار المساحة الفارغة الموجودة داخل المادة، وكيفية ترتيب تلك المساحة. لطالما استخدم العلماء أداتين رئيسيتين لوصف هذا العالم تحت الأرض؛ الأولى هي المسامية، وهي مقياس بسيط لمدى نسبة الفراغ في الصخر مقابل الصخر الصلب. والثانية هي النفاذية، التي تصف مدى سهولة تدفق السائل عبر تلك المساحة عندما يُدفع بواسطة الضغط. وبينما تُعد هاتان القيمتان ضروريتين للمهندسين والجيولوجيين، إلا أنهما لا ترويان سوى جزء من القصة؛ إذ يمكن لعينتين من الصخور أن تمتلكا نفس مقدار المساحة الفارغة تماماً، ومع ذلك تسمحان للسوائل بالمرور بسرعات متفاوتة للغاية لأن شكل المسام واتصالها يختلفان. هذا التباين الخفي، المعروف باسم الاضطراب الهندسي، يصعب قياسه ولكنه أمر بالغ الأهمية لفهم كيفية عمل باطن الأرض.
لقد وضع فريق من الباحثين في جامعة ساو باولو هدفاً يتمثل في إيجاد طريقة لقياس هذا الاضطراب دون الضياع في تعقيدات الصخور الحقيقية. فقاموا بإنشاء نموذج مبسط ثنائي الأبعاد لوسط مسامي، متخيلين إياه كلوح مسطح مليء بعوائق صلبة دائرية، مثل عملات معدنية مبعثرة على طاولة. وتمثل المساحات بين هذه العملات المسام التي يتدفق من خلالها السائل. ولفهم هذا الترتيب، استخدموا تقنية رياضية تسمى "تجزئة فورونوي" (Voronoi tessellation). تقسم هذه الطريقة المساحة الفارغة إلى مناطق متميزة، حيث تنتمي كل منطقة إلى أقرب عائق صلب. ومن خلال قياس مساحة هذه المناطق، استطاع الباحثون حساب رقم واحد، يُعرف باسم "إنتروبيا شانون" (Shannon entropy)، والذي يعمل كدرجة لمدى اضطراب الترتيب. فالمخطط المنتظم تماماً للعوائق سيكون لديه درجة صفر، بينما سيكون للبعثرة العشوائية تماماً درجة أعلى. ثم استخدموا محاكاة حاسوبية قوية لمراقبة كيفية تحرك سائل، وتحديداً محلول مائي ملحي، عبر هذه الترتيبات المختلفة.
اختبر الباحثون نموذجهم عن طريق تغيير عامل واحد في كل مرة مع تثبيت العوامل الأخرى. أولاً، بحثوا في كيفية تأثير موقع العوائق على التدفق؛ حيث بدأوا بشبكة مرتبة ومنظمة، ثم قاموا تدريجياً ببعثرة مواقع العوائق لزيادة مستوى الاضطراب. ووجدوا أنه كلما أصبح الترتيب أكثر فوضوية، أصبح المسار الذي يجب أن يسلكه السائل أكثر التواءً، وهي خاصية تُسمى "التعرج" (tortuosity). ومع ذلك، فإن القدرة الإجمالية للسائل على المرور عبر المادة، أي النفاذية، تغيرت قليلاً. فحتى عندما نُثرت العوائق بشكل عشوائي، وجد السائل طريقاً للعبور بسهولة تقارب سهولة عبوره في الشبكة المرتبة؛ حدث هذا لأنه بينما يؤدي الاضطراب إلى تضييق بعض المسارات، فإنه يوسع مسارات أخرى في الوقت نفسه، مما يوازن التدفق فعلياً.
بعد ذلك، فحص الفريق ما يحدث عندما يغيرون حجم العوائق أو إجمالي المساحة الفارغة. فعندما حافظوا على مستوى الاضطراب ثابتاً ولكنهم جعلوا العوائق أكبر حجماً، زادت النفاذية بشكل كبير، حيث ارتفعت بمقدار سبعة أضعاف تقريباً. حدث هذا لأن العوائق الأكبر حجماً تعني الحاجة إلى عدد أقل منها لملء المساحة، مما ترك فجوات أوسع بينها ليندفع السائل من خلالها. وبالمثل، عندما زادوا كمية المساحة الفارغة مع الحفاظ على حجم العوائق ثابتاً، قفزت النفاذية بأكثر من سبعين ضعفاً. أكدت هذه النتائج أن حجم الفجوات، أو "الأعناق" بين الأجزاء الصلبة، هو المحرك الأساسي لسرعة حركة السائل. أما درجة الاضطراب، فبالرغم من أنها تغير شكل المتاهة، إلا أنها لا تتحكم في سرعة التدفق بنفس الطريقة.
تخلص الدراسة إلى أن درجة الإنتروبيا الجديدة هي أداة قيمة، ولكن ليس للتنبؤ بسرعة التدفق؛ بل تعمل بدلاً من ذلك كوصف دقيق للهندسة نفسها. فهي تلتقط إعادة التنظيم المحلي لمساحة المسام، وتلخص المزيج المعقد من القنوات الواسعة والضيقة في متغير واحد. وهذا يسمح للعلماء بالتمييز بين صخرتين تبدوان متطابقتين من حيث المسامية وسرعة التدفق، ولكنهما تمتلكان هياكل داخلية مختلفة تماماً. وبينما يعد النموذج المستخدم في هذه الدراسة تمثيلاً مبسطاً ثنائي الأبعاد ولا يأخذ في الاعية بعد التعقيد الكامل للصخور ثلاثية الأبعاد أو التفاعل بين سوائل متعددة، إلا أنه نجح في عزل دور الاضطراب. ويشير العمل إلى أنه لفهم نقل السوائل في الأرض بشكل كامل، نحتاج ليس فقط إلى قياس مقدار المساحة المتوفرة، بل أيضاً إلى قيضاء البصمة الإحصائية المحددة لكيفية ترتيب تلك المساحة.
ملخص تقني: توصيف الفوضى الهندسية في الشبكات المسامية ثنائية الأبعاد باستخدام إنتروبيا فورونوي
بيان المشكلة لا يزال قياس انتقال السوائل في الأوساط المسامية يشكل تحدياً لأن الواصفات القياسية —المسامية (كسر الفراغ) والنفاذية (سهولة التدفق)— غير كافية. فالعينات ذات المسامية المتطابقة يمكن أن تظهر نفاذيات متفاوتة تماماً بسبب الاختلافات في الترتيب الهندسي للفراغات المسامية. وبينما تقيس المسامية مقدار المساحة الموجودة، فإنها لا تقيس كيفية ترتيب تلك المساحة. غالباً ما تفشل الواصفات الحالية في عزل الفوضى الهندسية بشكل مستقل عن المسامية أو حجم المسام. تعالج هذه الدراسة الفجوة في توصيف هذه الفوضى الهندسية باستخدام مقياس قياسي مشتق من نظرية المعلومات، وتحديداً اختبار ما إذا كان إنتروبيا شانون لتوزيع مساحة خلايا فورونوي يمكن أن يعمل كمقياس لعدم النظام (الفوضى) مستقل عن المسامية، وكيف يرتبط بخصائص التدفق أحادي الطور.
المنهجية استخدم المؤلفون نموذجاً مثالياً ثنائي الأبعاد يتكون من عوائق دائرية تمثل الطور الصلب لوسط مسامي. استخدمت الدراسة نهجاً هرمياً يجمع بين خصائص السوائل المستمدة من الديناميكا الجزيئية ومحاكاة على مستوى المسام:
توليد الشبكة: تم ملء نطاقات ثنائية الأبعاد (بطول جانبي قدره 5.68μm) بعوائق دائرية. قامت الدراسة بتغيير أربعة عوامل بشكل منهجي: الفوضى الموضعية، حجم العوائق، المسامية، وإنتروبيا فورونوي الناتجة. ولعزل المتغيرات، قام المؤلفون بتوليد 1,000 تحقق مستقل لكل تكوين، تتراوح من الشبكات المربعة المنتظمة إلى التوزيعات العشوائية تماماً، مع فرض حد أدنى من الفصل بين الأسطح لضمان الاتصال.
تحليل فورونوي: تم قياس الفوضى الهندسية باستخدام إنتروبيا شانون (H) لتوزيع مساحات خلايا فورونوي. تم بناء تقسيم فوروي بناءً على مراكز العوائق. تم حساب الإنتروبيا عن طريق تقسيم توزيع مساحات الخلايا الموزونة (A/⟨A⟩) إلى فئات (bins). هذا المقياس متميز عن إنتروبيا فوروي الطوبولوجية (القائمة على جوانب الخلايا) كونه حساساً لمدى انتشار أحجام الخلايا.
محاكاة التدفق: تمت محاكاة تدفق أحادي الطور باستخدام طريقة "لاتيس بولتزمان" (Lattice-Boltzmann method - LBM) عبر حزمة "Taxila". تم استيراد خصائص السوائل (الكثافة، اللزوجة الكينماتيكية) ومعلمات التفاعل بين السائل والصلب (زاوية التلامس) من دراسات سابقة للديناميكا الجزيئية لمحلو API ملحي في مصفوفة مونتموريلونيت. أُجريت المحاكاة في نظام التدفق الزاحف (رقم رينولدز <1) لضمان الاستجابة الخطية.
مقاييس الانتقال: من حقول السرعة في الحالة المستقرة، استخرج المؤلفون الالتواء الهيدروليكي (hydraulic tortuosity) ونفاذية دارسي.
المساهمات الرئيسية
عزل الفوضى الهندسية: نجحت الدراسة في فصل تأثيرات الفوضى الموضعية، وحجم العوائق، والمسامية على خصائص الانتقال، وهو فصل يصعب تحقيقه في عينات الصخور الطبيعية.
إنتروبيا مساحة خلية فورونوي كواصف: تقدم الورقة وتتحقق من صحة إنتروبيا شانون لمساحة خلية فورونوي كواصف قياسي مدمج للفوضى الهندسية، يمكن حسابه فقط من مواقع العوائق، وبشكل مستقل عن محاكاة التدفق.
توضيح محركات النفاذية: يوضح العمل أنه بينما تعد المسامية والنفاذية معايير قياسية، إلا أنهما تتركان دون قياس إعادة التنظيم الهندسي المحدد للمساحة المسامية. يلتقط مقياس الإنتروبيا هذا "إعادة التنظيم المحلي" حيث تؤدي الفوضى إلى توسيع بعض الممرات وتضييق أخرى عند ثبات المسامية.
النتائج
اعتماد النفاذية: وُجد أن النفاذية الكلية محكومة أساساً بمقياس مسام العنق (pore-throat scale). وقد أدى قانون قوة واحد (k∝wmed2.9) إلى دمج البيانات عبر التباينات في المسامية وحجم العوائق.
تأثير الفوضى الموضعية: عند ثبات المسامية وحجم العوائق، أدى زيادة الفوضى الموضعية (من الشبكة المنتظمة إلى العشوائية) إلى زيادة طفيفة في النفاذية (من ≈4.75 إلى $5.01mD)وزيادةفيالالتواء(من1.04إلى1.16$). عُزيت زيادة النفاذية إلى اتساع توزيع عرض الممرات، حيث تعوض ظهور قنوات أوسع تفضيلية عن تضيق القنوات الأخرى.
تأثير الإنتروبيا: عند تغيير الإنتروبيا وحدها ضمن الشبكات العشوائية تماماً (مع ثبات المسامية وحجم العوائق)، لم يُلاحظ أي اتجاه منتظم في النفاذية أو الالتواء؛ حيث ظلت التباينات ضمن التشتت الإحصائي للتكوينات. يشير هذا إلى أنه بمجرد فقدان النظام البلوري، فإن أي زيادة إضافية في الإنتروبيا لا تغير خصائص الانتقال بشكل كبير.
حجم العائق والمسامية: أدت زيادة حجم العائق (عند ثبات المسامية) إلى زيادة كبيرة في النفاذية (ما يقرب من سبعة أضعاف) وفي الالتواء. كما تسببت زيادة المسامية (عند ثبات حجم العائق) في ارتفاع حاد في النفاذية (أكثر من سبعين ضعفاً) وانخفاض في الالتواء.
الإنتروبيا مقابل الالتواء: ارتبطت الإنتروبيا بالالتواء خلال الانتقال من الشبكة المنتظمة إلى الحالة المضطربة. ومع ذلك، فإن هذا الارتباط هو جزئياً نتيجة اصطفاف الشبكة المنتظمة مع اتجاه التدفق، مما يقلل من طول المسار. في الحالات المضطربة تماماً، يضعف هذا الارتباط.
الأهمية والادعاءات تدعي الورقة أن إنتروبيا شانون لمساحة خلية فورونوي تعمل كـ واصف مكمل وليس كمتنبئ مباشر للنقل.
الدور المكمل: تقيس الإنتروبيا الفوضى الهندسية التي تتركها المسامية والنفاذية دون قياس. وهي تختزل تعقيد إعادة تنظيم المساحة المسامية في متغير واحد.
القيود التنبؤية: تنص الدراسة صراحة على أن الإنتروبيا لا تتنبأ بالنفاذية بشكل أفضل من مقياس مسام العنق. ويتماشى الاعتماد الضعيف للنفاذية على الإنتروبيا مع نظرية انتقال المسار الحرج، والتي تشير إلى أن التدفق محكوم بمقياس الاختناق (أوسع المسام المتصلة) بدلاً من إحصائيات الترتيب العالمي.
التكرار: لاحظ المؤلفون أن الإنتروبيا مكررة رياضياً مع معامل الاختلاف (CV) لمساحات الخلايا بالنسبة لتوزيعات "غاما" الملاحظة في هذه التعبئات. وبالتالي، فهي تعمل كموجز رقمي واحد لمدى انتشار مساحة الخلية.
النطاق: تقتصر الاستنتاجات بدقة على النماذج المثالية ثنائية الأبعاد، والتدفق أحادي الطور في نظام التدفق الزاحف عبر عوائق دائرية. ولا يدعي المؤلفون أن هذه النتائج تنطبق مباشرة على صخور الخزانات المعقدة ثلاثية الأبعاد أو متعددة الأطوار دون مزيد من التحقق، رغم إمكانية حساب الإنتروبيا من مراكز الحبيبات في سير عمل الصخور الرقمية.
باختصار، يثبت هذا العمل أنه بينما تغير الفوضى الهندسية (المقاسة بإنتروبيا فورونوي) الطوبولوجيا المحلية للمساحة المسامية، فإن النفاذية الكلية في هذا النظام تهيمن عليها مقاييس مسام العنق. يعد مقياس الإنتروبيا مفيداً لتوصيف ترتيب الوسط، لكنه لا يحل محل المسامية أو النفاذية كواصف أساسي للنقل.