ALTYN Inversion: A MATLAB-Based 2D/2.5D Finite-Strike Gravity Inversion Framework with Box-Constrained IRLS Regularization
تقدم هذه الورقة ALTYN، وهو إطار عمل شفاف لقلب الجاذبية ثنائي الأبعاد/2.5D يعتمد على برنامج MATLAB ويستخدم تنظيم خوارزمية المربعات الصغرى التكرارية (IRLS) المقيدة بصندوق وحلال نيوتن المسقط لتوليد نماذج كثافة مقبولة جيوفيزيائياً، مما يثبت من خلال الاختبارات الاصطناعية والميدانية أنه بينما تنجح الطريقة بفعالية في مطابقة البيانات وإبراز السمات الهيكلية المشتركة، إلا أنها لا تستطيع إعادة بناء الهندسات المعقدة لباطن الأرض بشكل فريد بسبب عدم التفرد المتأصل.
إن الأرض ليست كرة متجانسة من الصخور؛ بل هي كعكة معقدة متعددة الطبقات من مواد ذات أوزان مختلفة. فبعض الأجزاء كثيفة وثقيلة، مثل الجرانيت الصلب أو الخامات الغنية بالمعادن، بينما أجزاء أخرى خفيفة ومسامية، مثل الرماد البركاني المتصدع أو التربة المشبعة بالمياه. ولأن الجاذبية قوة تجذب الكتلة، فإن هذه الاختلافات الخفية في الكثافة تخلق تباينات طفيفة يمكن قياسها في قوة الجاذبية عند السطح. يستخدم علماء الجيوفيزياء هذه الحقيقة للنظر فيما تحت أقدامهم دون الحاجة إلى الحفر. فهم يقيسون قوة الجاذبية في نقاط عديدة على طول خط معين، ثم يحاولون العمل بشكل عكسي لمعرفة ما يوجد هناك. وتُعرف هذه العملية، المعروفة باسم "القلب الجاذبي" (gravity inversion)، بأنها تشبه محاولة تخمين شكل جسم مخفي داخل صندوق من خلال الشعور فقط بمدى ثقل الصندوق من الخارج. وهي أداة قوية للعثور على النفط والمعادن والطاقة الحرارية الأرضية، لكنها أيضاً عملية معقدة للغاية. إذ يمكن غالباً لنفس النمط من الجاذبية أن ينتج عن العديد من الأشكال والأعماق المختلفة للصخور، مما يعني عدم وجود إجابة واحدة واضحة لما يكمن في الأسفل.
ولمواجهة هذا التحدي، طور باحث يدعى أحمد علي موخاروم إطاراً رقمياً جديداً يسمى "ألثين" (ALTYN)، صُمم لرسم خرائط لهذه التباينات الخفية في الكثافة بوضوح وأمانة أكبر. تم بناء هذه الأداة للتعامل مع نوع محدد من المشكلات الجيولوجية: عندما تمتد التشكيلات الصخرية جانبياً ولكن ليس إلى الأبد. ففي كثير من السيناريوهات الواقعية، لا يمتد طبقة من الصخور أو جسم بركاني إلى ما لا نهاية في كلا الاتجاهين كما افترضت النماذج الأقدم والأبسط؛ بل يكون لها طول محدد، مثل لوح حجري طويل وسميك ينتهي بعد مسافة معينة. ويأخذ نظام "ألثين" هذا في الاعتبار من خلال التعامل مع باطن الأرض كسلسلة من الكتل المستطيلة التي لها عرض محدد، مما يسمح للكمبيوتر بحساب قوة الجاذبية بدقة أكبر مما لو افترض أن الصخر يستمر إلى ما لا نهاية. كما يحترم النظام الشكل الفعلي للأرض فوقها، حيث يعدل حساباته لتتطابق مع التلال والوديان بدلاً من افتراض أن السطح مستوٍ.
اختبر الباحثون هذا الإطار الجديد باستخدام ثلاثة نماذج محاكاة حاسوبية تم إنشاؤها بعناية، حيث عرفوا الإجابة الدقيقة مسبقاً. في الاختبار الأبسط، أنشأوا كتلة واحدة مدمجة من الصخور الثقيلة تحت الأرض. نجح نظام "ألثين" في تحديد مكان الكتلة وحدد بشكل صحيح أنها ثقيلة، لكنه لم يستطع إعادة إنشاء حوافها الحادة بدقة؛ وبدلاً من ذلك، أظهرت الصورة المستعادة أن الصخر الثقيل كان على شكل ضبابي وممتد قليلاً. وهذا قيد شائع في دراسات الجاذبية: يمكن للنظام أن يخبرك بمكان وجود جسم ثقيل، لكنه يواجه صعوبة في تحديد مدى عمقه بدقة أو مدى حدة حدوده. وفي اختبار ثانٍ يتضمن كتلتين منفصلتين، إحداهما ثقيلة والأخرى خفيفة، وجد النظام كلتيهما وميز طبيعتهما المتضادة، ومع ذلك ظهرت الكتلة الخفيفة أكثر عمقاً وانتشاراً مما كانت عليه في الواقع. أما الاختبار الأكثر صعوبة فقد تضمن ترتيباً معقداً من الصخور الخفيفة والثقيلة المتراكمة فوق بعضها البعض تحت سطح تلي. وفي هذا الاختبار، تمكن النظام من مطابقة قياسات الجاذبية تماماً، ومع ذلك فشل في إعادة بناء الصخر الثقيل الموجود في مكان أعمق في موقعه الحقيقي؛ وبدلاً من ذلك، قام بتوزيع الوزن بطريقة أخرى أنتجت نفس القراءة السطحية. يسلط هذا النتيجة الضوء على حقيقة جوهرية: وهي أن المطابقة المثالية لبيانات الجاذبية لا تضمن بالضرورة أن الصورة المستخرجة لباطن الأرض صحيحة.
ولمعرفة كيف تؤثر الاختيارات الرياضية المختلفة على الصورة النهائية، طبق الفريق الأداة على مسح جاذبية حقيقي من حقل "ديينج" للطاقة الحرارية الأرضية في إندونيسيا. هذه المنطقة عبارة عن مشهد بركاني يمتد فيه خط طويل من نقاط القياس لمسافة تزيد عن 9 كيلومترات. أجرى الباحثون البيانات عبر النظام مرتين، باستخدام نهجين رياضيين مختلفين لتنعيم النتائج. النهج الأول، المعروف باسم "التنظيم من النوع L2"، أنتج نموذجاً كانت فيه تغيرات الكثة تدريجية ومنتشرة، مثل لوحة مائية ناعمة. أما النهج الثاني، المسمى "IRLS"، فقد سمح بوجود حدود أكثر حدة وتميزاً، مما خلق نموذجاً يبدو أكثر شبهاً بالهيكل المقطعي أو الكتل المتراصة. وقد ناسب كلا الأسلوبين قياسات الجواذبية الفعلية بنفس القدر، محققين المتطلبات الإحصائية للمطابقة الجيدة. ومع ذلك، بدت الخريطتان الناتجتان لباطن الأرض مختلفتين تماماً في تفاصيلهما. واعتبرت المناطق التي اتفقت فيها الخريطتان — وتحديداً موقع منطقة ثقيلة كبيرة تحت جبل ومنطقة واسعة خفيفة بالقرب منها — هي النتائج الأكثر موثوقية. أما الاختلافات بين الخريطتين، مثل مدى عمق المنطقة الثقيلة أو مدى حدة حوافها، فقد تم تحديدها على أنها أقل يقيناً.
تخلص الدراسة إلى أنه بينما يعد إطار "ألثين" طريقة شفافة وفعالة لبناء نماذج جيولوجية معقولة، إلا أنه لا يمكنه الكشف سحرياً عن صورة فريدة ومثالية لما تحت السطح. يُنظر إلى النظام كأداة لبناء ومقارنة نسخ مختلفة ممكنة لما تحت الأرض، وليس كآلة تخرج حقيقة واحدة وحتمية. إن المعلومات الأكثر موثوقية تأتي من السمات العامة التي تظل ثابتة عبر الأساليب المختلفة، مثل الموقع العام والوزن للأجسام الصخرية الرئيسية. أما التفاصيل الدقيقة، مثل العمق الدقيق لصدع ما أو الشكل الدقيق لتداخل بركاني، فتظل مفتوحة للتفسير وتتطلب أنواعاً أخرى من البيانات للتأكيد. وفي النهاية، يوفر هذا العمل طريقة أكثر وضوحاً وأمانة للنظر في كثافة الأرض الخفية، مذكراً إيانا بأنه حتى مع الرياضيات المتقدمة، فإن أعماق الأرض غالباً ما تحتفظ ببعض أسرارها.
ملخص تقني: إطار عمل ALTYN للقلب (Inversion)
بيان المشكلة يعتبر قلب الجاذبية (Gravity inversion) بطبيعته عملية غير محددة بشكل جيد وغير فريدة؛ حيث يمكن لتوزيعات كثافة تحت سطحية متميزة أن تولد مجالات جاذبية متطابقة أو متطابقة تقريبًا عند السطح. ويتفاقم هذا الغموض بسبب المقايضة بين مطابقة البيانات وتعقيد النموذج، وتلاشي الحساسية مع العمق، وصعوبة تمثيل الهياكل الجيولوجية ذات أطوال الضرب (strike lengths) المحدودة بدلاً من الامتداد اللانهائي. غالبًا ما تفترض نماذج 2D التقليدية ضربًا لانهائيًا، بينما تنطوي عمليات القلب ثلاثية الأبعاد الكاملة (3D) على تكاليف حوسبة عالية. علاوة على ذلك، تستخدم الأدوات الحالية مثل عائلة برمجيات GROWTH نهج "نمو النموذج" حيث يتم تعيين تباينات كثافة محددة مسبقًا لخلايا مختارة لبناء الهندسة. هناك حاجة إلى إطار عمل شفاف ومقيد لتقدير تباينات الكثافة المستمرة على شبكة ثابتة مع مراعاة فيزياء الضرب المحدود، والتضاريس، واستراتيجيات التنظيم (regularization) المختلفة صراحةً، دون ادعاء إعادة بناء جيولوجية فريدة.
المنهجية يقدم البحث إطار عمل ALTYN Inversion، وهو إطار عمل يعتمد على MATLAB ثنائي الأبعاد/2.5D مصمم لقلب بيانات الجاذبية للمقاطع (profile gravity inversion). تدمج المنهجية عدة مكونات رئيسية:
النمذجة الأمامية (2.5D ذات الضرب المحدود): يتم تقسيم ما تحت السطح إلى شبكة X–Z ثنائية الأبعاد حيث تُعامل كل خلية كمنشور مستطيل ثلاثي الأبعاد بـ طول ضرب محدد (Ly). يستخدم العامل الأمامي تعبيرات تحليلية للجذب الجاذبي الرأسي لمنشور مستطيل (Nagy et al., 2000)، مع مراعاة الامتداد المحدود العمودي على المقطع. يتجنب هذا افتراض الضرب اللانهائي للنماذج القياسية 2D دون الحاجة إلى بارامترية كثافة ثلاثية الأبعاد كاملة.
التضاريس والتقسيم (Discretization): يستخدم النموذج شبكة X–Z متوافقة مع التضاريس. يتم استبعاد الخلايا فوق السطح، وتُقطع الخلايا النشطة العلوية لتتطابق مع ارتفاعات السطح المحلية. يتم تمرير هذه الهندسة مباشرة إلى نواة العمليات الأمامية (forward kernel)، مما يضمن الاتساق بين التضاريس، ومصفوفة الحساسية، وعناصر التنظيم.
مطابقة البيانات وعدم اليقين: يستخدم إطار العمل مطابقة بيانات مرجحة بعدم اليقين. إذا لم تتوفر انحرافات معيارية محددة، فإنه يبني نموذج عدم يقين مشترك (مطلق-نسبي). تعمل دالة الهدف على تقليل حد عدم مطابقة البيانات المرجح (ϕd) جنبًا إلى جنب مع حد تنظيم النموذج (ϕm).
استراتيجيات التنظيم: تم تنفيذ مخططين رئيسيين للتنظيم:
التنظيم من النوع L2: نهج البنية الدنيا التربيعي (P=2) الذي يفضل التباينات الكثافية الموزعة والناعمة.
خوارزمية IRLS (المربعات الصغرى الموزونة تكراريًا): نهج غير تربيعي (P<2) لتعزيز النماذج المركزة أو الثابتة جزئيًا. يسمح التنفيذ بالتحكم الاتجاهي؛ فعلى سبيل المثال، يستخدم تطبيق المجال Pz=0.5 للتدرجات الرأسية مع الحفاظ على Px=Ps=2، مما يشجع على الانتقالات الرأسية الحادة مع الحفاظ على النعومة الأفقية.
القيود والمحلل (Solver): يفرض القلب قيود الصندوق الصريحة (حدود الكثافة) على كل خلية نشطة أثناء عملية التحسين، بدلاً من قص النتائج بعد عملية القلب. يستخدم محلل Newton/active-set preconditioned conjugate-gradient (PCG) لعملية التحسين. يتم تقدير مكون الإزعاج الإقليمي (اتجاه ثابت أو خطي) في وقت واحد مع نموذج الكثافة لمنع التوجهات ذات الطول الموجي الطويل من أن تُفسر خطأً كبنية تحت سطحية.
التقارب: تستخدم الخوارزمية جدول تبريد لمعامل مقايضة التنظيم (β) وتنتهي بناءً على معيار التفاوت الطبيعي (الهدف WRMS ≈ 1) أو تغير النموذج النسبي.
المساهمات الرئيسية المساهمة الأساسية لـ ALTYN ليست تقديم منظم جديد، بل هي التكامل الشفاف للمفاهيم الراسخة في إطار عمل واحد قابل لإعادة الإنتاج لقلب المقاطع. تشمل الميزات الرئيسية ما يلي:
الجمع بين فيزياء المنشور ذي الضرب المحدود وتقسيم الخلايا النشطة المتوافق مع التضاريس.
التقدير المتزامن لتباينات الكثافة ومعاملات الإزعاج الإقليمية ضمن حلقة تحسين مقيدة بالصندوق.
تحليل مقارن لتنظيم L2 و IRLS على نفس مجموعة البيانات لإظراك كيف تختار المقدمات (priors) المختلفة نماذجًا مقبولة جيوفيزيائيًا.
التمييز الصارم بين المطابقة في فضاء البيانات (البواقي) والاسترداد في فضاء النموذج (الدقة الهندسية)، مع التأكيد على أن المطابقة الجيدة للبيانات لا تضمن نموذجًا جيولوجيًا صحيحًا.
النتائج تم التحقق من صحة إطار العمل من خلال ثلاثة اختبارات معيارية اصطناعية تزداد تعقيدًا، وطُبق على مقطع جاذبية مكون من 171 محطة (9.23 كم) في منطقة ديينج (Dieng) الحرارية الأرضية، إندونيسيا.
الاختبارات الاصطناعية:
المعيار 1 (شذوذ بسيط): نجح في استعادة الموقع الجانبي والقطبية لمصدر مدمج، لكنه أظهر تلطيخًا رأسيًا كبيرًا، وهي سمة من سمات عمليات الحد الأدنى للبنية.
المعيار 2 (قطبية مزدوجة): حدد بشكل صحيح الأجسام الموجبة والسالبة المنفصلة والاتجاه الإقليمي، رغم أن السعات تعرضت للتوهن واتسعت الامتدادات الرأسية.
المعيار 3 (هيكل معقد): أظهر حدود عدم اليقيد. وبينما حقق القلب توافقًا في WRMS يتوافق مع مستوى الضجيج، إلا أنه فشل في إعادة بناء التدخل الإيجابي العميق بشكل متميز في موقعه الحقيقي، حيث قام بدلاً من ذلك بتوزيع الكثافة على حجم أكبر. أكد هذا أن المطابقة المقبولة إحصائيًا يمكن أن تتعايش مع خطأ كبير في فضاء النموذج.
التطبيق الميداني (ديينج):
أنتج كل من حلي L2 و IRLS نتائج حققت معيار التفاوت (WRMS ≈ 0.95–0.97) وحددا نفس المعالم الأولية: منطقة كثافة موجبة قوية تحت 2–3 كم (مرتبطة بجبل براو)، ومنطقة موجبة ثانوية عند 6–7 كم، مفصولة بنطاقات سالبة واسعة.
حل L2: أنتج تباينات كثافة أكثر نعومة وتوزيعًا مع انتقالات تدريجية.
حل IRLS: أنتج انتقالات رأسية أكثر حدة وهياكل كثافة عمودية أكثر بسبب حد التدرج الرأسي غير التربيعي.
البواقي (Residuals): أظهر كلا النموذجين عدم مطابقة موضعي (البواقي المعيارية > 4σ) بالقرب من أعلى تدرج للجاذبية (2.6–3.0 كم)، مما يشير إلى أن بارامتريّة 2.5D أو النموذج الإقليمي قد يكون غير كافٍ لهذا الشذوذ عالي السعة تحديدًا.
الأهمية والادعاءات يضع البحث ALTYN بتواضع ليس كأداة لإعادة بناء جيولوجية فريدة، بل كإطار عمل شفاف لبناء ومقارنة نماذج كثافة مقبولة جيوفيزيائيًا.
عدم اليقين (Non-Uniqueness): تؤكد الدراسة أن قلب الجاذبية غير فريد بطبيعته. تظهر النتائج الاصطناعية صراحةً أن تلبية معيار تفاوت البيانات لا يثبت صحة الهندسة المستعادة.
الاستقرار مقابل التفاصيل: يزعم المؤلف أن المعالم الأولية (قطبية الشذوذ والموقع الجانبي التقريبي) أكثر دفاعًا واستقرارًا عبر مخططات التنظيم المختلفة من الحدود التفصيلية، أو قيم الكثافة القصوى، أو الامتدادات العميقة.
الاعتماد على المقدمات (Prior Dependence): تسلط النتائج الضوء على أن النموذج المستعاد يعتمد بشدة على الافتراضات المسبقة، بما في ذلك نوع التنظيم، وقيود الضرب، ودقة الشبكة. لا يدعي البحث المتانة الكمية تجاه اضطرابات هذه المعلمات، حيث لم يتم إجراء تجارب اضطراب مخصصة.
التفسير الجيولوجي: بالنسبة لتطبيق ديينج، يخلص المؤلف إلى أن النماذج المستعادة تمثل تباينات كثافة نسبية تتوافق مع البيانات والافتراضات. وبينما تعد هذه التباينات معقولة جيولوجيًا (مثل: الشذوذات الموجبة التي قد تشير إلى أجسام اقتحامية، والشذوذات السالبة التي قد تشير إلى تغيرات أو مواد بركانية فتاتية)، يذكر البحث صراحةً أن هذه التفسيرات تظل فرضيات تتطلب قيودًا مستقلة (مثل بيانات الآبار، أو الموجات الكهرومغناطيسية، أو الزلازل) للتحقق منها. يعمل إطار العمل على جعل الافتراضات الكامنة وراء عملية القلب صريحة بدلاً من إنتاج مقطع عرضي جيولوجي "صحيح" واحد.