A numerical framework for Newtonian-noise estimation at the Einstein Telescope: 2-D simulations beyond the plane-wave approximation

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

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

يهدف تلسكوب أينشتاين (ET)، وهو مرصد موجات جاذبية مستقبلي من الجيل الثالث يقع تحت الأرض، إلى اكتشاف الإشارات وصولاً إلى بضعة هرتزات. ومع ذلك، يُتوقع أن تكون حساسيته عند الترددات المنخفضة (تحديداً بين 1.7 و6 هرتز) محدودة بـ ضجيج نيوتن (NN). ينشأ ضجيج نيوتن من تقلبات قوة الجاذبية على كتل الاختبار الخاصة بالكاشف، والناتجة عن تقلبات الكثافة الناجمة عن النشاط الزلزالي في الصخور المحيطة.

تعتمد طرق تقدير ضجيج نيوتن الحالية بشكل كبير على نماذج تحليلية أو شبه تحليلية تفترض:

• وسطاً متجانساً أو طبقات أفقية.
• تقريبات الموجة المستوية ذات التماسك اللانهائي.
• إهمال عدم التجانس الجيولوجي، والتشتت، وتحول الأنماط، والتفاعلات الموجية المعقدة.

هذه التبسيطات تفشل في رصد التعقيدات الجيولوجية الواقعية، مثل التشتت من الواجهات الجيولوجية، والارتباطات المحددة بين الموجات الطولية (P-waves) والموجات المستعرضة (S-waves). علاوة على ذلك، فإن الحصول على بيانات مصفوفة زلزالية كثيفة مطلوبة لمخططات الإلغاء القائمة على التنبؤ بنقطة واحدة يعد أمراً صعباً من الناحية العملية في المنشآت الموجودة تحت الأرض. هناك حاجة ماسة إلى إطار عمل عددي كامل قادر على حل عدم التجانس الجيولوجي وحقول الموجات المعقدة لتوفير تقديرات خاصة بالموقع لضجيج نيوتن.

2. المنهجية

قام المؤلفون بتطوير إطار عمل عددي يعتمد على محاكاة العناصر الطيفية (spectral-element simulations) لنمذجة حقول الموجات الزلزالية وضجيج نيوتن اللاحق.

• المحلل العددي: تستخدم الدراسة برنامج Salvus، وهو محلل عناصر طيفية عالي التوازي ومعزز بوحدات معالجة الرسومات (GPU)، مما يجعله مناسباً للأوساط غير المتجانسة عبر الجمع بين القدرة على التكيف الهندسي للعناصر المحدودة والدقة العالية للطرق الطيفية.
• النموذج الفيزيائي:
  • معادلة الموجة: يتم حل حقل الموجة الزلزالية $\vec{\xi}(\vec{x}, t)$ باستخدام معادلة حفظ الزخم في وسط متصل.
  • تقلبات الكثافة: يتم اشتقاق تقلبات الكثافة $\delta\rho$ من حقل الإزاحة باستخدام معادلة الاستمرارية: $\delta\rho = -\vec{\nabla} \cdot (\rho \vec{\xi})$. وقد تم تفكيكها إلى حد مركزي (bulk term) (تغيرات الكثافة داخل الوسط) وحد سطحي (surface term) (قفزات الكثافة عند الواجهات، مثل جدران الكهوف).
  • حساب قوة نيوتن: يتم حساب تقلب قوة الجاذبية $\delta \vec{F}_M$ على كتلة الاختبار من خلال تكامل تقلبات الكثافة عبر الحجم (المعادلة 3). كما استخدم المؤلفون "معادلة ثنائية القطب" (Dipole equation) (المعادلة 8) التي توحد التكاملات المركزية والسطحية، مع مراعاة وجود كهف افتراضي حول كتلة الاختبار.
• سيناريوهات المحاكاة (إثبات مفهوم ثنائي الأبعاد):
  1. مرجع المصدر الواحد: مصدر قوة واحد يتجه للأعلى (Ricker wavelet، بتردد 2 هرتز) على سطح وسط متجانس للتحقق من النتائج العددية مقابل تنبؤات الموجة المستوية التحليلية.
  2. مصفوفة الضوضاء المحيطة: مصفوفة من 30 مصدراً سطحياً عشوائياً بمدخلات زمنية مستقلة (محدودة النطاق حتى 10 هرتز) لتقريب الضوضاء الزلزالية المحيطة المستقرة.
• المعلمات: تستخدم المحاكاة وسطاً مرجعياً علوياً متجانساً ($c_P = 3000$ م/ث، $c_S \approx 1764$ م/ث، $\rho \approx 2224$ كجم/م³) في نطاق مساحته 20 كم $\times$ 8 كم مع حدود امتصاص.

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

• تطوير إطار عمل عددي: إنشاء سير عمل لمحاكاة الحقول الموجية الزلزالية وحساب ضجيج نيوتن مباشرة من حقل الإزاحة، متجاوزين بذلك التقريبات التحليلية.
• التحقق مقابل النماذج التحليلية: نجحت عملية مطابقة النتائج العددية مع تنبؤات الموجة المستوية التحليلية في وسط متجانس، مما أظهر توافقاً ممتازاً في كل من المساهمات المركزية ومساهمات الكهف.
• تحديد نسبة الموجات الطولية (P-wave fraction): تقديم طريقة لتقدير نسبة الموجات الطولية ($p$) في حقل موجي متعدد المصادر عبر تحليل نسبة كثافة الطيف القدرتي لتباعد (divergence) وحمل (curl) حقل الإزاحة.
• إثبات إمكانية تعزيز التخفيف: تقديم دليل على أن محتوى الموجات الطولية في الضوضاء الزلزالية المحيطة قد يكون أقل بكثير مما كان يُفترض سابقاً، مما يشير إلى آفاق أفضل لاستراتيجيات تخفيف ضجيج نيوتن.

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

• التحقق (المصدر الواحد):
  • أعادت المحاكاة العددية لمصدر واحد إنتاج التنبؤات التحليلية بدقة.
  • طابق التكامل المركزي (المنحنى الأزرق) تنبؤ الموجة الطولية التحليلي، بينما نجح تكامل ثنائي القطب (المنحنى البرتقالي)، الذي يتضمن تأثير الكهف الافتراضي، في رصد مساهمات الموجات الطولية والمستعرضة مجتمعة.
  • أظهرت النت أن مساهمة الموجات المستعرضة (S-waves) تنشأ فقط من الحد السطحي (جدار الكهف)، حيث لا تسبب الموجات المستعرضة تقلبات في الكثافة المركزية في وسط متجانس.
• محاكاة الضوضاء المحيطة (30 مصدراً):
  • وقعت كثافة الطيف السعة (ASD) لضجيج نيوتن المحسوب ضمن الحدود النظرية العليا (الموجات الطولية النقية) والدنيا (الموجات المستعرضة النقية)، لكنها تذبذبت بشكل كبير، مما يشير إلى ارتباطات معقدة بين أنواع الموجات لا تغطيها النماذج التحليلية البسيطة.
  • نسبة الموجات الطولية ($p$): قدرت الدراسة نسبة الموجات الطولية في الحقل الموجي بالقرب من كتلة الاختبار. ووجد أن متوسط الوسيط لمتوسط التردد هو $p \approx 0.137 \pm 0.02$.
  • هذه القيمة أقل بكثير من القيمة المفترضة شائعاً وهي $p = 1/3$ في التقديرات التحليلية.

5. الأهمية والآفاق المستقبلية

• الآثار على التخفيف: بما أن الموجات الطولية (P-waves) عادة ما تكون أصعب في التخفيف من الموجات المستعرضة (S-waves) (بسبب قدرتها على إحداث تقلبات في الكثافة المركزية)، فإن انخفاض نسبة الموجات الطولية ($p \approx 0.14$ مقابل $0.33$) يعني أن مستوى ضجيج نيوتن قد يكون أقل مما تم تقديره سابقاً، وأن إمكانية التخفيف الناجح (مثل استخدام مصفوفات المستشعرات الزلزالية) أعلى.
• الاتجاهات المستقبلية:
  • تعمل هذه الدراسة ثنائية الأبعاد المتجانسة كإثبات للمفهوم. الإطار مصمم للامتداد إلى محاكاة ثلاثية الأبعاد.
  • سيعمل العمل المستقبلي على دمج نماذج جيولوجية محلية مفصلة، وتضاريس، وعدم تجانس لحساب التشتت وتحولات الأنماط.
  • الهدف النهائي هو توليد تقديرات خاصة بالموقع لضوة نيوتن في تلسكوب أينشتاين لتحسين تصميم المصفوفات الزلزالية المخصصة لإلغاء الضوضاء.

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