A Parallel Imaging Simulation Framework for Lunar-based
Radio Telescope Arrays
تقترح هذه الورقة إطار عمل للمعالجة المتوازية باستخدام مجموعات من وحدات المعالجة المركزية (CPU) ووحدات معالجة الرسومات (GPU) لتعزيز كفاءة ودقة محاكاة مصفوفات التلسكوبات الراديوية الضخمة القائمة على القمر، وتحديداً لدعم مهمة "خطة هونغ مينغ" من خلال تحسين محاكة الخط القاعدي، ونمذجة الرصد، وإعادة بناء الصور.
المؤلفون الأصليون:Xiang Zhao, Ran Duan, Li Deng, You Song
لعقود من الزمن، ظل علماء الفلك عميانًا عن ركن شاسع وصامت من الكون. فأسفل تردد 30 مليون دورة في الثانية، يعمل غلاف الأرض الجوي كستارة ثقيلة، تحني وتمتص الهمسات اللاسلكية الخافتة القادمة من الكون. يحمل هذا النطاق المخفي مفاتيح فهم "العصور المظلمة الكونية"، وهي الحقبة الغامضة التي سبقت اشتعال النجوم الأولى، وما تلاها من "فجر كوني" عندما بدأت تلك النجوم في إضاءة الكون. ولرؤية ذلك، يجب على العلماء النظر من الفضاء، بعيدًا عن التداخل الناتج عن كوكبنا. وتتخيل إحدى الخطط الطموحة، المعروفة باسم "خطة هونغمينغ" (Hongmeng Plan)، أسطولًا من الأقمار الصناعية التي تدور حول القمر لإنشاء تلسكوب راديوي ضخم عائم. ومع ذلك، فإن الحجم الهائل للبيانات التي ستولدها مثل هذه المهمة مذهل، مما يتطلب طريقة جديدة لمعالجة المعلومات لا تستطيع الحواسيب الحالية التعامل معها بمفردها.
ولحل هذه المشكلة، طور الباحثون شيانغ تشاو، وران دوان، ولي دينج، ويو سونغ، إطار عمل برمجياً قوياً مصمماً لمحاكاة كيفية عمل هذا التلسكوب القمري. ويركز عملهم على التحدي الحسابي الهائل المتمثل في تحويل الإشارات الراديوية الخام إلى صور واضحة للسماء بأكملها. في التلسكوب الراديوي التقليدي، تُجمع الإشارات بواسطة أزواج من الهوائيات، وتحدد المسافة بينهما — والتي تسمى "خط الأساس" — تفاصيل الصورة. تستخدم خطة "هونغمينغ" ثمانية أقمار صناعية صغيرة تتحرك في رقصة منسقة، وتغير مسافاتها من بعضها البعض باستمرار لبناء خريطة ثلاثية الأبعاد للسماء. ولأن هذه الأقمار الصناعية تتحرك باستمرار، ولأن القمر يحجب الضوضاء الراديوية للأرض فقط في جزء من مدارها، فإن البيانات التي تجمعها معقدة للغاية وغير موزعة بشكل منتظم.
لقد أنشأ الفريق محاكاة تحاكي العملية بأكملها لهذه المهمة، بدءاً من حساب مواقع الأقمار الصناعية وصولاً إلى إعادة بناء الصورة النهائية. ووجدوا أن محاولة القيام بذلك على حاسوب واحد ستكون بطيئة بشكل مستحيل. وبدلاً من ذلك، بنوا نظاماً يقسم العمل عبر معالجات متعددة، باستخدام كل من الرقائق الحاسوبية القياسية وبطاقات الرسوميات المتخصصة للتعامل مع أجزاء مختلفة من المهمة. وللقيام بالعمل الشاق المتمثل في حساب مواقع الأقمار الصناعية وعدّ كم عدد نقاط البيانات التي تقع في مناطق محددة من السماء، استخدموا معالج حاسوب متعدد النواة. أما بالنسبة للعمل الرياضي المكثف المتمثل في تحويل تلك البيانات إلى صورة، فقد سخروا قوة بطاقات الرسوميات، المصممة لأداء ملايين العمليات الحسابية في وقت واحد. وقد سمح هذا النهج الهجين بمعالجة الكميات الهائلة من البيانات الناتجة عن المهمة في جزء بسيط من الوقت الذي كانت ستستغرقه لولا ذلك.
يعد عدم الانتظام في توزيع البيانات عقبة رئيسية في هذه العملية، حيث تدور الأقمار الصناعية بزاوية محددة، مما يعني أنها لا تستطيع رؤية كل أجزاء السماء بنفس الكفاءة. ويخلق هذا فجوات في البيانات، مما يؤدي إلى صور تبدو أكثر سطوعاً في بعض المناطق وأقل سطوعاً في مناطق أخرى، بغض النظر عما هو موجود فعلياً هناك. واكتشف الباحثون أن هذا التفاوت لم يكن عيباً في تصميم التلسكوب، بل كان نتيجة طبيعية لمسار الأقزاء الصناعية. ولإصلاح ذلك، طوروا طريقة تصحيحية تحلل توزيع نقاط البيانات وتعدل سطوع الصورة النهائية وفقاً لذلك. ومن خلال تقسيم السماء إلى مناطق مختلفة وتطبيق تعديلات رياضية محددة لكل منها، تمكنوا من تنعيم هذه البقع الساطعة الاصطناعية وإظهار صورة أكثر دقة للكون.
كانت نتائج محاكاتهم مذهلة. فعندما اختبروا النظام عند تردد 10 ملايين دورة في الثانية، نجح إطار العمل في إعادة بناء صورة كاملة للسماء في غض_ن ساعات، وهي مهمة كانت ستستغرق أياماً أو أسابيع باستخدام الطرق القديمة. وأظهرت الصور المصححة مستوى عالياً من التفاصيل، حيث حافظت على أشكال وهياكل المصادر الراديوية الكونية مع إزالة التشوهات الناتجة عن تشكيل الأقمار الصناعية. وقاس الفريق جودة هذه الصور باستخدام اختبارات الدقة القياسية ووجدوا أن النسخ المصححة كانت أقرب بكثير إلى السماء الحقيقية من النسخ غير المصححة. ويشير هذا النجاح إلى أن إطار العمل المقترح ليس مجرد تمرين نظري، بل هو أداة عملية يمكن استخدامها لتوجيه مهمة "هونغمينغ" الفعلية عند انطلاقها.
بينما تعد هذه المحاكاة خطوة كبيرة للأمام، يقر الباحثون بأن النظام ليس مثالياً بعد. فالنظام الحالي لم يتم اختباره بعد عند أعلى الترددات التي قد تصل إليها المهمة، كما أنه لا يأخذ في الاعتبار بعد حجب القمر لأجزاء من الرؤية خلال أوقات معينة. ومع ذلك، فإن العمل يثبت أنه من خلال الجمع الصحيح بين القدرة الحاسوبية والخوارزميات الذكية، يمكن التغلب على تحديات البيانات الضخمة لعلم الفلك الراديوي القمري. ومن خلال توفير مسار واضح لمعالجة هذه الإشارات المعقدة، يضع إطار العمل هذا الأساس التقني لمستقبل يمكننا فيه أخيراً رفع الستار عن اللحظات الأولى للكون.
ملخص تقني: إطار محاكاة التصوير المتوازي لمصفوفات التلسكوبات الراديوية القائمة على القمر
1. بيان المشكلة
تعتبر عمليات الرصد الراديوي الأرضي تحت مستوى 30 ميجاهرتز مقيدة بشدة بسبب الانكسار الأيوني والتداخل الراديوي (RFI)، مما يخلق عنق زجاجة لدراسة العصور المظلمة الكونية وعصر إعادة التأين عبر إشعاع الهيدروجين المحايد (21 سم) منزاح نحو الأحمر. وللتغلب على ذلك، تقترح "خطة هونغ مينغ" (Hongmeng Plan) نشر مصفوفة أقمار صناعية خطية في مدار قمري لإجراء تصوير كامل السماء.
ومع ذلك، تواجه هذه المهمة تحديات حوسبية كبيرة:
حجم بيانات هائل: يتضمن تصميم المهمة 8 أقمار صناعية فرعية تشكل خطوط أساس (baselines) تخضع لتحولات تلسكوبية "تشبه الزنبرك" على مدار مهمة مدتها 3 سنوات. عند تردد 30 ميجاهرتز، يُتوقع أن يتجاوز العدد الإجمالي لخطوط الأساس 3 مليارات.
هندسة تصوير معقدة: تستخدم المصفوفة هوائيات شاملة الاتجاهات مع خطوط أساس غير مستوية في مستوى واحد ومجال رؤية واسع. وتعتبر تقريبات تحويل فورييه ثنائية الأبعاد التقليدية غير صالحة؛ حيث يجب مراعاة حد الـ w (تأثيرات خطوط الأساس غير المستوية في المستوى) باستخدام تحويلات فورييه ثلاثية الأبعاد.
عيوب التصوير: بسبب ميل المدار القمري بمقدار 30 درجة، يكون توزيع خطوط الأساس في نطاق التردد ثلاثي الأبعاد غير مكتمل، مما يترك فجوات مخروطية (بزاوية نصف مخروط قدرها 60 درجة) على طول المحور القطبي. يؤدي هذا إلى عدم انتظام يعتمد على خط العرض في صور درجة حرارة السطوع المعاد بناؤها، وتحديداً في درجات الحرارة المرتفعة ضمن نطاق ±30 درجة عرض.
عنق الزجاجة الحوسبي: تتطلب معالجة البيانات بمقياس البيتابايت لمثل هذه المهام قدرات حوسبة فائقة (exascale). والأطر الحالية المصممة للمصفوفات الأرضية ليست قابلة للتطبيق مباشرة على الديناميكيات المدارية وخصائص البيانات الفريدة لخطة "هونغ مينغ".
2. المنهجية
يقترح المؤلفون إطار عمل للحوسبة المتوازية غير المتجانسة يدمج بنيات متعددة من المعالجات المركزية (CPU) ومعالجات الرسوميات (GPU) للتعامل مع خط معالجة التصوير بأكمله: محاكاة خط الأساس، وإحصاءات أعداد خطوط الأساس، ومحاكاة الرصد، وإعادة بناء الصورة.
2.1 بنية سير العمل المتوازي
ينقسم إطار العمل إلى أربع مراحل أساسية، محسنة لقدرات عتادية محددة:
محاكاة خط الأساس (متعدد الـ GPU):
يتم توليد مواقع الأقمار الصناعية لكل دورة تحول تلسكوبي (14 يوماً).
تُحسب المواقع النسبية بين 8 أقمار صناعية فرعية لتشكيل متجهات خط الأساس.
التوازي: تعالج وحدات معالجة رسوميات متعددة دورات مختلفة بالتوازي. داخل وحدة GPU واحدة، تشمل المهام توليد مواقع الأقمار الصناعية، وحساب خطوط الأساس الصالحة، وإجراء عمليات الفرز وإزالة التكرار باستخدام مكتبة Thrust (مجموعة أدوات CUDA) لإدارة مصفوفات خط الأساس العالمية بكفاءة.
إحصاءات أعداد خطوط الأساس (متعدد الـ CPU):
تحسب هذه المرحلة معاملات الوزن (Δck) بناءً على كثافة خطوط الأساس في شبكات التردد ثلاثية الأبعاد (المقسمة حسب 0.5 طول موجي).
اختيار العتاد: رغم التوازي العالي لوحدات GPU، اختار المؤلفون نهج متعدد المعالجات المركزية (Multi-CPU) باستخدام لغة بايثون (Pandas) وتعدد الخيوط (multi-threading).
التعليل: نظراً لقيود ذاكرة GPU، فإن تقسيم مجموعة البيانات الضخمة (مليارات خطوط الأساس) إلى كتل GPU (M) سيتطلب وقتاً طويلاً جداً لنقل البيانات مقارنة بنهج الـ CPU (N كتل). أثبت إطار عمل الـ CPU أنه أكثر كفاءة لهذه المهمة المحددة المتعلقة بالإدخال/الإخراج والتجميع، حيث يعالج البيانات في أجزاء أكبر مع عدد أقل من التكرارات.
محاكاة الرصد (متعدد الـ GPU):
تحاكي دالات الرؤية (V) بناءً على خريطة درجة حرارة السطوع كاملة السماء (I) والإحداثيات السماوية (l,m,n).
التحسين: يستخدم MPI و OpenMP للتحكم المتوازي عبر وحدات GPU.
تنفيذ النواة (Kernel): تُخزن الإحداثيات السماوية وبيانات سطوع السماء في ذاكرة المشاركة (shared memory) الخاصة بـ GPU لتقليل زمن التأخير. تنفذ نواة محاكاة خط الأساس تحويل فورييه (المعادلة 2) بالتوازي، مع استخدام تقنية "فك الحلقات" (loop unrolling) لتسريع الحسابات.
إعادة بناء الصورة (متعدد الـ GPU):
تعيد بناء صورة السماء الكاملة عبر تحويل فورييه العكسي (المعادلة 5).
التحسين: يتم توزيع بيانات خطوط الأساس، والرؤى، والأعداد عبر وحدات GPU. يتم تحميل البيانات القابلة لإعادة الاستخدام (خطوط الأساس، الرؤى، الإحداثيات) مسبقاً في الذاكرة المشتركة. تُنفذ نواة التحويل العكسي بالتوازي مع استخدام تقنية فك الحلقات.
2.2 استراتيجية تصحيح الصورة
لمعالجة عيوب درجة حرارة السطوع المعتمدة على خط العرض والناتجة عن عدم اكتمال تغطية خطوط الأساس:
تقسيم المناطق: يتم تقسيم السماء إلى منطقة مركزية (C، ≈±30∘)، ومنطقة انتقالية (T)، ومنطقة خطوط العرض العالية (H) بناءً على كثافة خطوط الأساس.
معاملات التصحيح:
المنطقة المركزية: يتم حساب معامل تصحيح كنسبة بين متوسط خطوط الأساس في المناطق غير المركزية إلى تلك الموجودة في المنطقة المركزية.
المنطقة الانتقالية: تُشتق المعاملات من نسبة متوسط درجات حرارة السطوع بين مناطق خطوط العرض العالية والمناطق الانتقالية.
المعايرة التكيفية: يطبق خط المعالجة مرشح وسيط (median filter) واسع النطاق لعزل المكونات الخلفية، ويحسب تسلسل كسب عرضي (latitudinal gain sequence)، ويطبق تنعيماً غاوسياً وعامل تخميد (α=0.5) لإنشاء مجال كسب لكل بكسل. يُستخدم هذا المجال للتعويض التفاضلي لمكونات الخلفية والهيكل، يليه إزاحة شدة عالمية وتقليم للقيم الشاذة.
3. المساهمات الرئيسية
إطار عمل غير متجانس: تطوير إطار محاكاة قابل للتوسع يخصص المهام استراتيجياً للمعالجات المركزية (لتجميع البيانات الضخمة/الإحصاء) ومعالجات الرسوميات (للتحويلات الفورية وخطوط الأساس المكثفة حوسبياً).
تصوير كامل السماء ثلاثي الأبعاد: تكييف طرق تحويل فورييه ثلاثية الأبعاد خصيصاً للهندسة غير المستوية واسعة المجال للمرصد المداري القمري، متجاوزاً تقريبات 2D الأرضية.
التصحيح المدفوع بخطوط الأساس: اقتراح طريقة مبتكرة لتصحيح درجة حرارة السطوع تستخدم التوزيع المحدد لخطوط الأساس للتخفيف من عيوب التصوير الناتجة عن الميل المداري، مما يقلل بفعالية من عدم الانتظام المعتمد على خط العرض.
تحسين الأداء: تنفيذ تحسينات محددة لـ CUDA، بما في ذلك استخدام الذاكرة المشتركة للبيانات القابلة لإعادة الاستخدام، وفك الحلقات، واستخدام مكتبة Thrust للفرز وإزالة التكرار بكفاءة.
4. النتائج
أُجريت التجارب باستخدام عنقود حوسبي يحتوي على 2 من معالجات Intel Xeon Gold 8336C و3 من معالجات NVIDIA GeForce RTX 4090.
الكفاءة:
عند 1 ميجاهرتز، أكمل خط المعالجة الكامل (محاكاة الرصد + إعادة بناء الصورة) في حوالي 52 دقيقة (25 دقيقة + 27 دقيقة).
عند 10 ميجاهرتز، نجح إطار العمل في إعادة بناء الصور في 344 ساعة (167 ساعة + 177 ساعة) باستخدام 3 وحدات GPU. وهذا يثبت جدوى معالجة البيانات عالية الدقة (Nside=4096، ~200 مليون بكسل) والتي ستكون مستعصية باستخدام الطرق التسلسلية.
جودة الصورة:
المقاييس الكمية: حسنت طريقة تصحيح الصورة الدقة بشكل كبير.
SSIM: ارتفع من 0.7847 إلى 0.9369 عند 1 ميجاهرتز ومن 0.8195 إلى 0.9290 عند 10 ميجاهرتز.
PSNR: تحسن من ~20.2 ديسيبل إلى 29.5 ديسيبل (1 ميجاهرتز) ومن ~20.9 ديسيبل إلى 28.1 ديسيبل (10 ميجاهرتز).
MSE: انخفض بمقدار يقارب عشرة أضعاف.
التحليل البصري: نجح التصحيح في كبح درجات حرارة السطوع المرتفعة في المنطقة المركزية (±30∘). وبينما تظل بعض العيوب الطفيفة موجودة بالقرب من حدود المناطق المركزية والانتقالية، فقد تم استعادة الدقة المورفولوجية العامة.
5. الأهمية والادعاءات
يزعم البحث أن هذا العمل يوفر أساساً تقنياً حاسماً لـ خطة هونغ مينغ. ومن خلال إثبات وجود إطار عمل متوازٍ قابل للتوسع قادر على التعامل مع أحجام البيانات الهائلة والقيود الهندسية المعقدة للمصفوفة المدارية القمرية، يضع المؤلفون مساراً لـ:
المعالجة في الوقت الفعلي وعالية الدقة: تمكين معالجة البيانات بمقياس البيتابايت المطلوبة لمشاريع علم الفلك الراديوي المستقبلية ذات الحوسبة الفائقة.
الجدوى العلمية: التحقق من أن التصوير عالي الدقة لكامل السماء للعصور المظلمة الكونية هو أمر ممكن حوسبياً باستخدام الأجهزة غير المتجانسة الحالية، بشرما تم تطبيق تصحيحات محددة للهندسة المدارية.
مرجع منهجي: تقديم بنية مرجعية لمهام التداخل القمرية المستقبلية، لا سيما في التعامل مع خطوط الأساس غير المستوية وإحصاءات البيانات واسعة النطاق.
يقر المؤلفون بالقيود، مشيرين إلى أن المحاكاة الحالية لا تتضمن بعد تأثيرات الاحتجاب القمري، وأن المحاكاة كاملة الدقة عند 30 ميجاهرتز مقيدة حالياً بذاكرة الـ GPU، مما يتطلب عملاً مستقبلياً في الحوسبة الموزعة وتخصيص النواة (kernel customization).