2.5-D Electrical Resistivity Forward Modelling with Undulating Topography using a Modified Half-Space Analytical Solution

إليك ملخص تقني مفصل للورقة البحثية بعنوان: "النمذجة الأمامية للمقاومة الكهربائية ثنائية ونصف الأبعاد (2.5-D) مع تضاريس متموجة باستخدام حل تحليلي معدل لنصف الفضاء."

1. بيان المشكلة

تُستخدم تصوير المقاومة الكهربائية للتيار المستمر (DC) على نطاق واسع لتوصيف باطن الأرض، ولكنها تُجرى غالباً فوق تضاريس متموجة. وتعد النمذجة الأمامية الدقيقة أمراً ضرورياً لعمليات الانعكاس (inversion) الموثوقة، ومع ذلك فإن دمج التضاريس المعقدة يفرض تحديات عددية كبيرة:

• إزالة التفرد (Singularity Removal): تواجه الطرق العددية القياسية صعوبة في التعامل مع تفرد الجهد الكهربائي عند أقطاب المصدر النقطي. ويُستخدم تقنية تقسيم "الجهد الأولي-الثانوي" (primary-secondary potential split) الشائعة، حيث يتم حل الجهد الأولي المفرد تحليلياً، بينما يُحل الجهد الثانوي السلس عددياً.
• محدودية الطرق الحالية: تعتمد النهج الحالية (مثل Penz et al., 2013) على جهد أولي تحليلي مشتق لـ نصف فضاء متجانس مسطح. وهذا الافتراض يكون صالحاً فقط عندما تكون التضاريس بالقرب من المصدر سلسة ومستوية.
• المشكلة الجوهرية: في التضاريس الطبيعية ذات الزوايا الحادة، أو الانحناء العالي، أو عدم استمرارية المنحدرات، أو حتى عند تمثيل التضاريس السلسة باستخدام عناصر محدودة خطية (linear finite elements)، يفشل افتراض السطح المسطح.
  • تقوم العناصر الخطية بتقريب الأسطح المنحنية إلى قطع مستقيمة خطية مجزأة، مما يخلق زوايا اصطناعية على شكل حرف "V" عند المصدر.
  • الحل التحليلي لنصف الفضاء المسطح لا يحقق شرط نيومان الحدّي (تدفق صفري) على هذه الهندسات التي تتخذ شكل حرف "V".
  • يؤدي هذا التباين الهندسي إلى أخطاء نمذجة جسيمة (تتجاوز غالباً 1%)، وهي أخطاء تستمر حتى مع تحسين الشبكة (mesh refinement) أو تنعيم التضاريس هندسياً.

2. المنهجية

يقترح المؤلفون استراتيجية معدلة لإزالة التفرد تشتق جهداً أولياً تحليلياً جديداً مصمماً خصيصاً للهندسة المحلية للمصدر.

• المعادلات الحاكمة: تستخدم الدراسة صياغة 2.5-D للمقاومة الكهربائية للتيار المستمر، حيث يتم تطبيق تحويل جيب التمام فوريه (Fourier cosine transform) على معادلة بواسون ثلاثية الأبعاد لتقليلها إلى مسألة ثنائية الأبعاد في المستوى $x-z$.
• تقسيم الأولي-الثانوي: يتم تقسيم الجهد الكلي $u$ إلى جهد أولي ($u_p$) وجهد ثانوي ($u_s$).
• اشتقاق الجهد الأولي المعدل:
  • بدلاً من افتراض نصف فضاء مسطح (الزاوية المجسمة $S = 2\pi$)، يشتق المؤلفون حلاً تحليلياً لمصدر نقطي يقع عند رأس وتد على شكل حرف V (رُقْعَة كروية/spherical lune).
  • يأخذ الحل في الاعتبار الزاوية المجسمة ($S$) التي تشغلها التضاريس المحلية عند المصدر.
  • التعبير التحليلي الجديد للجهد الأولي هو:
    $$u_p(r) = \frac{I}{\sigma_0 S r}$$
    وفي مجال فوريه لـ 2.5-D:
    $$\tilde{u}_p(x, k_y, z) = \frac{I}{\sigma_0 S} K_0(k_y r_{xz})$$
    حيث $S = 4\pi - 2\beta$ (مع كون $\beta$ هي زاوية الوتد)، وبالنسبة للسطح المسطح، تعود $S$ إلى $2\pi$.
• التنفيذ:
  • يُستخدم الجهد الأولي المعدل لتحديد الشروط الحدية للجهد الثانوي.
  • تم تنفيذ الطريقة ضمن إطار طريقة العناصر المحدودة (FEM) باستخدام عناصر خطية (Q4) وعناصر من رتب أعلى (Q9, Q16).
  • يتم حساب الزاوية المجسمة $S$ صراحةً بناءً على هندسة العناصر المتصلة بعقدة المصدر، مما يستوعب كل من التجزئات الخطية والمنحنية.

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

1. تحديد التباين الهندسي: تثبت الورقة البحثية بصرامة أن أخطاء النمذجة في التضاريس المتموجة تنبع من عدم التطابق بين المجال الأولي للسطح المسطح والزاوية المجسمة الحقيقية التي تشغلها التضاريس المجزأة (خاصة مع العناصر الخطية).
2. الصياغة التحليلية المبتكرة: اشتقاق حل تحليلي مغلق للجهد الأولي لمصدر نقطي في وتد على شكل حرف V، مما يعمم حل نصف الفضاء القياسي.
3. المتانة مع العناصر الخطية: تسمح هذه الطريقة باستخدام شبكات عناصر محدودة خطية خشنة دون الحاجة إلى تنعيم اصطناعي للتضاريس أو تكثيف مفرط للشبكة بالقرب من المصدر، وهي الحلول البديلة المكلفة حاسوبياً في الطرق الموجودة.
4. الاستقلال عن الخطأ: يظهر التحليل النظري أن الخطأ في النهج التقليدي للسطح المسطح يعتمد فقط على نسبة الزاوية المجسمة الحقيقية إلى $2\pi$ ($|1 - S/2\pi|$)، مما يجعل الخطأ ثابتاً بغض النظر عن المسافة من المصدر.

4. النتائج

تم التحقق من صحة الطريقة المقترحة من خلال تجارب عددية على نماذج مسطحة، ونموذج خندق على شكل حرف V، ونموذج تلة ووادي جيبية:

• النماذج المتجانسة والطبقية المسطحة: حققت الطريقة أخطاء نسبية أقل من 0.1% مقابل الحلول التحليلية، مما يؤكد صحة التنفيذ.
• نموذج الخندق على شكل حرف V:
  • الطريقة التقليدية (افتراض السطح المسطح): أنتجت أخطاء بنسبة ~16.5% بسبب عدم استمرارية المنحدر الحاد عند المصدر.
  • الطريقة المقترحة: حققت أخطاء أقل من 0.1% حتى مع الشبكات الخشنة (حجم عنصر 1 متر) والعناصر الخطية.
  • ظل الخطأ ثابتاً عبر النطاق، مما يؤكد صحة الاشتقاق النظري.
• نموذج التلة والوادي الجيبي:
  • الطريقة التقليدية: تراوحت الأخطاء بين 1% إلى 2% مع العناصر الخطية، ووصلت إلى 10% عندما كان المصدر عند التقاء منحدر، حتى مع استخدام عناصر من رتب أعلى.
  • الطريقة المقترحة: حققت باستمرار أخطاء ≤ 0.1% عبر جميع أنواع الشبكات (الخطية والمنحنية) ومواقع المصدر (وادي، سطح مستوٍ، أو التقاء منحدر).
  • ومن الجدير بالذكر أن الطريقة المقترحة باستخدام العناصر الخطية تفوقت على الطريقة التقليدية باستخدام العناصر المنحنية من رتب أعلى في بعض السيناريوهات.

5. الأهمية والخاتمة

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

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