Optimization of a cosmic muon tomography scanner for cargo border control inspection

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

بيان المشكلة
يعد التصوير المقطعي بالميونات الناتجة عن الأشعة الكونية، وتحديداً التصوير المقطعي بتشتت الميونات (MST)، تقنية واعدة للأمن الحدودي، قادرة على اكتشاف المواد ذات العدد الذري المرتفع (Z) مثل المتفجرات والمواد النو|^{المهربة داخل حاويات الشحن. ومع ذلك، يعتمد أداء أنظمة MST بشكل كبير على تكوين الكاشف؛ إذ يمكن للتصاميم غير المثالية أن تؤدي إلى انخفاض في دقة التتبع، وانخفاض معدلات الكشف، وعدم كفاءة أوقات المسح. وبينما اعتمدت جهود التحسين السابقة على المحاكاة مقترنة بالبحث الشبكي المنفصل، لم يتم استكشاف استراتيجية تحسين شاملة ومنتظمة تربط مباشرة بين معايير تصميم الكاشف وأداء التمييز بين المواد بشكل كامل. يعالج هذا العمل الحاجة إلى تحسين تصميم ماسح "SilentBorder"، وهو نظام MST نمطي، لتعظيم فعاليته في تحديد التهديدات المخفية.

المنهجية
تستخدم الدراسة نهجاً مزدوج الطريقة يجمع بين محاكاة عالية الدقة وتقنيات التحسين التفاضلية:

1. إطار عمل محاكاة GEANT4:

• تم بناء بيئة محاكاة تفصيلية باستخدام GEANT4 لنمذجة هندسة ماسح SilentBorder (48 هودوسكوب مرتبة في شبكة 3×4 على جوانب الحاوية الطويلة) وسيناريوهات شحن متنوعة (بما في ذلك المعادن، والمتفجرات، والمخدرات، والبضائع القانونية).
• قام مولد MuSiBo بمحاكاة تدفق الميونات الناتجة عن الأشعة الكونية، مع مراعاة التحلل وانحناء الأرض.
• طبقت محاكاة الكاشف الزائف قيوداً واقعية، بما في ذلك التلطيخ المكاني (1 ملم) وعتبات الطاقة (0.4 ميجا إلكترون فولت)، لنمذجة ضربات الكاشف.
• استخدمت عملية إعادة البناء خوارزمية "نقطة التقارب الأقصى" (PoCA) لتحديد رؤوس التشتت.
• الدراسات التي أجريت:
  • تموضع الهودوسكوب: تم تغيير المسافات بين ألواح الكاشف في الأبعاد $x$ و $y$ و $z$ لتقييم التأثيرات على كفاءة جمع الميونات والدقة الزاوية.
  • ضربات الجسيمات الثانوية: تم تحليل تأثير الجسيمات الثانوية (بشكل أساس قبل أشعة دلتا) على إعادة بناء الأحداث، وتقييم ما إذا كانت تضعف التتبع أو توفر تمييزاً إضافياً للمواد.

2. TomOpt والتحسين البايزي (Bayesian Optimization):

• TomOpt: تم استخدام إطار عمل محاكاة تفاضلي قائم على لغة بايثون لتحسين هندسة الكاشف عبر نزول التدرج. وهو يدمس توليد الميونات، والانتشار، ونمذجة التشتت (بناءً على صيغ مجموعة بيانات الجسيمات PDG)، وإعادة بناء الحجم.
• خوارزميات إعادة البناء: تم اختبار طريقتين للاستدلال:
  • تقدير طول الإشعاع ($X_0$): يستخدم جذر متوسط مربعات زوايا التشتت لاستنتاج كثافة المادة.
  • خوارزمية التجميع المبوب (BCA): تقوم بتجميع رؤوس التشتت بناءً على التقارب المكاني وزوايا التشتت لاستنتاج ملفات الكثافة.
• التحسين البايزي (BO): لمعالجة الضوضاء المتأصلة في استدلال الحجم (خاصة القيم المتطرفة في تقدير $X_0$)، تم دمج وحدة BO باستخدام مكتبة Ax. وعلى عكس نزول التدرج، يستخدم BO نموذج بديل لعملية غاوس (Gaussian Process) ودالة استعلام "التحسين المتوقع" للتنقل في مساحات الأهداف الصاخبة. استهدف التحسين متوسط مربع الخطأ (MSE) بين توزيعات المواد المستنتجة والحقيقية.

النتائج الرئيسية

• تموضع الهودوسكوب (GEANT4):

  • الفجوات داخل المستوى ($x, y$): يؤدي تقليل الفجوات بين الهودوسكوبات إلى تحسين كفاءة النظام من خلال زيادة حجم الميونات التي يتم اعتراضها. ومع ذلك، فإن زيادة الفجوات يقلل من نقاط PoCA القابلة لإعادة البناء داخل الشحنة، مما يستلزم أوقات تعرض أطول لنفس مستوى الحساسية.
  • المباعدة الرأسية ($z$): تؤدي زيادة الفصل الرأسي بين الألواح داخل الهودوسكوب إلى تحسين الدقة الزاوية بشكل كبير بسبب تأثير "ذراع الرافعة". تشير الدراسة إلى أن زيادة المسافة حتى ~20 سم توفر مقايضة مجدية، حيث تحسن الدقة مع تقليل الكفاءة بشكل طفيف فقط. بعد ذلك، تستقر مكاسب الدقة.

• الجسيمات الثانوية (GEANT4):

  • نُسبت حوالي 8.3% من الضربات المسجلة في الهودوسكوبات إلى نشاط ثانوي غير الميونات.
  • وجدت الدراسة أن الضربات الثانوية لا توفر وسيلة للتمييز بين أنواع الشحنات المختلفة (مثل الموز مقابل الكوكايين).
  • والأهم من ذلك، أن وجود الضربات الثانوية يسبب تدهوراً محدوداً جداً في تتبع الميونات. فحتى مع عتبات الطاقة المنخفضة، تظل "نسبة التلوث" (الأحداث التي تعيق فيها الثنائيات مسارات الميون) أقل من 6%. وبالتالي، فإن اتخاذ تدابير منفصلة لاستبعاد الثنائيات ليس ضرورياً تماماً لتكوين الماسح الحالي.

• أداء التحسين (TomOpt و BO):

  • BO مقابل نزول التدرج: نجح BO في تحسين تكوينات الكاشف في البيئات الصاخبة حيث قد يواجه نزول التدرج صعوبة. بالنسبة لطريقة BCA، قلل BO من دالة الهدف بنسبة 13.3% في حوالي 26 دقيقة (أسرع من وقت تشغيل نزول التدرج الذي كان أطول بنسبة ~38% بسبب التكلفة الحسابية لـ BCA). وبالنسبة لطريقة $X_0$، حقق BO انخفاضاً بنسبة 22.8% في حوالي 11 دقيقة.
  • متانة الاستدلال: أظهرت طريقة BCA قدرة أكبر على تحمل الضوضاء وأداءً متفوقاً (MSE أقل، ومقياس تشابه بنيوي أعلى) مقارنة بطريقة تقدير $X_0$، خاصة في التمييز بين المواد منخفضة Z.
  • التكوينات المحسنة: نجحت عملية BO في توجيه وضع الهودوسكوبات لتعظيم تغطية الحجم الخامل. ومع ذلك، لم تتفوق التكوينات المحسنة بشكل كبير على التصميم المرجعي للكاشف المستمد من الحدس البشري.

الأهمية والادعاءات
يزعم البحث تقديم مسار واضح لتحسين ماسحات التصوير المقطعي بالميونات المستقبلية لرقابة الحدود. وتتمثل مساهماته الرئيسية في:

1. إرشادات التصميم: إثبات أن تقليل الفجوات داخل المستوى وتحسين المباعدة الرأسية (حتى ~20 سم) أمران حاسمان لموازنة الكفاءة والدقة الزاوية.
2. تقييم الجسيمات الثانوية: إثبات أن الجسيمات الثانوية لا تضعف جودة إعادة البناء بشكل كبير في إعداد SilentBorder، مما يبسط متطلبات تصميم الكاشف.
3. إطار عمل التحسين: تقديم وحدة تحسين بايزي ضمن إطار عمل TomOpt للتعامل مع دالات الأهداف الصاخبة، مما يوفر بديلاً قوياً للطرق القائمة على التدرج للمهام المعقدة لإعادة البناء.
4. المقارنة المنهجية: تسليط الضوء على تفوق خوارزمية التجميع المبوب (BCA) على تقدير طول الإشعاع البسيط من حيث الاستقرار ضد الضوضاء واستخدامها لأغراض التحسين.

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