A Spectral-Domain Pseudo-Inverse Method for True 3D Gravity Inversion

ملخص تقني: طريقة المعكوس الزائف في المجال الطيفي لعكس الجاذبية ثلاثي الأبعاد الحقيقي

بيان المشكلة
تتناول الورقة التحديات الجوهرية المتأصلة في عكس الجاذبية ثلاثي الأبعاد: عدم التحديد، وعدم الاستقرار، وضعف الدقة الرأسية. تفتقر الملاحظات السطية التقليدية، التي توفر بيانات فقط على مستوى ثنائي الأبعاد ($z=0$)، إلى معلومات الموجة العددية الرأسية اللازمة لحل توزيع الكثافة تحت السطحي ثلاثي الأبعاد $\rho(x,y,z)$. يؤدي هذا إلى عدد غير كافٍ من القيود المستقلة. علاوة على ذلك، يتلاشى مجال الجاذبية بسرعة مع العمق، مما يؤدي إلى إضعاف الإشارات القادمة من المصادر العميقة بشدة، مما يجعل مشكلة العكس سيئة التحديد وحساسة للغاية للضوضاء. إن النهج الحالي، مثل التنظيم التكراري في المجال المكاني (مثل تيكهونوف)، مكلف حاسوبياً وحساس للنماذج الأولية، بينما تقتصر طرق المجال الترددي عادةً على إعادة بناء الواجهات بدلاً من العكس الحجمي الكامل.

المنهجية
يقترح المؤلف طريقة عكس "ثلاثية الأبعاد حقيقية" تعتمد على نظرية المعكوس الزائف في المجال الطيفي للأنظمة القابلة للقطرية الوحدوية. تسير المنهجية عبر مستويين نظريين متميزين:

1. إكمال فضاء البيانات عبر الاستمرار التحليلي:
   لتحقيق متطلبات المعكوس الزائف في المجال الطيفي للبيانات ثلاثية الأبعاد، تستغل الورقة حقيقة أن شذوذ الجاذبية $g_z$ يحقق معادلة لابلاس ($\nabla^2 g_z = 0$) في طبقة الهواء الخالية من المصدر ($z < 0$). ومن خلال تطبيق تحويل فورييه ثنائي الأبعاد على بيانات السطح $g_z(x,y,0)$، يشتق المؤلف صيغة استمرار صاعد تحليلي. يسمح هذا ببناء حجم بيانات ثلاثي الأبعاد كامل ومتسق فيزيائياً فوق مستوى المراقبة. والأهم من ذلك، تشتق الورقة تعبيراً تحليلياً عاماً لطيف نصف الفضاء:
   $$\tilde{g}^{obs}_z(k_x, k_y, k_z) = \tilde{g}_z(k_x, k_y, 0) \frac{1}{a - i k_z}$$
   حيث $a = \sqrt{k_x^2 + k_y^2}$. يوفر هذا التعبير فصل الطيف الأفقي عن نواة انتشار رأسية عالمية، مما يوفر أخذ عينات الموجة العددية الرأسية اللازمة لتحويل فورييه ثلاثي الأبعاد دون أخطاء الاستيفاء العددي.

2. العكس بالمعكوس الزائف في المجال الطيفي:
   يُظهر أن عامل التمرير ثلاثي الأبعاد للجاذبية قابل للقطرية الوحدوية في مجال فورييه. وتُعبر العلاقة بين طيف الكثافة $\tilde{\rho}$ وطيف الجاذبية $\tilde{g}_z$ كما يلي:
   $$\tilde{g}_z(\mathbf{k}) = \lambda(\mathbf{k}) \tilde{\rho}(\mathbf{k}), \quad \text{حيث } \lambda(\mathbf{k}) = \frac{-i 4\pi G k_z}{K^2}$$
   حيث $K^2 = k_x^2 + k_y^2 + k_z^2$. وبناءً على نظرية المؤلف السابقة للأنظمة القابلة للقطرية الوحدوية، يتم إنشاء مرشح معكوس معمّم مستقر:
   $$q_\alpha(\lambda) = \frac{\bar{\lambda}}{|\lambda|^2 + \alpha}$$
   حيث $\alpha > 0$ هو معامل التنظيم. يُطبق هذا المرشح على طيف البيانات ثلاثي الأبعاد الناتج للحصول على طيف الكثافة. ويتم استعادة توزيع الكثافة المكاني النهائي عبر تحويل فورييه ثلاثي الأبعاد العكسي. ونظراً لبنية القطب في طيف البيانات (تحديداً القطب عند $k_z = -ia$)، يتم تنفيذ تحويل فورييه $k_z$ العكسي عبر تكامل مسار في المستوى المركب، مغلق في النصف السفلي من أجل $z > 0$.

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

• إكمال فضاء البيانات: تثبت الورقة أن الاستمرار الصاعد عبر معادلة لابلاس يمكنه بناء حجم بيانات ثلاثي الأبعاد بشكل صارم من بيانات سطحية ثنائية الأبعاد، مما يخفف من عدم التحديد في المشكلة ويوفر أخذ عينات الموجة العددية الرأسية المطلوبة للتحليل الطيفي ثلاثي الأبعاد.
• العكس ذو الصيغة المغلقة: تثبت الورقة أن عامل التمرير ثلاثي الأبعاد للجاذبية قابل للقطرية الوحدوية تماماً. وبناءً على ذلك، يتقلص العكس إلى عملية حسابية حتمية ذات صيغة مغلقة تتكون من تحويل فورييه ثلاثي الأبعاد أمامي واحد، وخطوة قياس طيفي باستخدام المرشح المشتق، وتحويل عكسي واحد، مما يلغي الحاجة إلى التحسين التكراري.
• الاستقرار والاتساق: ثبت أن المرشح العكسي المنشأ يتمتع باستقرار محدود ($\|q_\alpha\|_\infty \leq 1/(2\sqrt{\alpha})$) واتساق (يتقارب إلى معكوس مور-بينروز المعمم عندما $\alpha \to 0^+$).
• الإطار التحليلي: يوفر اشتقاق طيف نصف الفضاء (المعادلة 7) إطار مدخلات موحد يفصل المكونات الأفقية والرأسية، وهو قابل للتطبيق على أي بيانات جاذبية سطحية بغض النظر عن تعقيد النموذج تحت السطحي.

النتائج
تم التحقق من صحة الطريقة باستخدام نموذج كرة متجانسة. وتشير النتائج إلى أن:

• عملية العكس تستعيد الموقع ثلاثي الأبعاد للشذوذ بشكل صحيح.
• مع تقارب معامل التنظيم $\alpha \to 0^+$، تنجح الطريقة في إزالة تأثير انتشار دالة غرين (Green's function spreading effect)، مستعادةً السلوك المنفرد (دالة ديراك دلتا) عند نقطة المصدر.
• بالنسبة لقيمة $\alpha$ المحدودة، يكون الحل توزيعاً منعمًا، ويحافظ التكامل عبر الفضاء بأكمله بدقة على الكتلة الإجمالية للمصدر.
• تستعيد الطريقة كثافة الاضطراب (نسبة إلى خلفية غير معروفة) بمتوسط صفري؛ تتطلب الكثافات المطلقة معلومات سابقة مستقلة. يعتمد استعادة توزيع الكثافة الكمي على اختيار $\alpha$: حيث تؤدي $\alpha$ الأصغر إلى نتائج أكثر حدة وقيم داخلية أكثر دقة، مقتربة من الحد المتباعد عند المصدر.

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