A Computational Framework for Pulmonary Assessing Wave Intensity Following Simulated Lung Resection

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

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

لا يزال سرطان الرئة سبباً رئيسياً للوفاة الناجمة عن السرطان، ويعد استئصال الرئة وسيلة علاجية أساسية. وبينما يعد الاستئصال فعالاً لإزالة الورم، فمن المعروف أنه يضعف وظيفة البطين الأيمن (RV). والآلية الكامنة وراء هذا الضعف غير واضحة؛ فرغم الاشتباه في زيادة الحمل البعدي (afterload)، إلا أن المقاييس التقليدية مثل ضغط الشريان الرئوي (PA) والمقاومة الوعائية غالباً ما تفشل في إظهار زيادات ملحوية بعد الجراحة.

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

2. المنهجية

تستخدم الدراسة مساراً حوسبياً متعدد المراحل يجمع بين معالجة الصور، ونظرية المخططات (graph theory)، والنمذجة الهيدروديناميكية أحادية الأبعاد (1D)، وتحليل شدة الموجة.

أ. الحصول على البيانات والمعالجة المسبقة

• مصدر البيانات: 48 شبكة سطحية للشريان الرئوي (بصيغة STL) مستمدة من صور الأشعة المقطعية لمرضى لا يعانون من أمراض رئوية (من دراسة Goubergrits وآخرون).
• بناء النطاق أحادي الأبعاد (1D):
  • استخراج الخط المركزي (Centerline Extraction): استُخدمت مجموعة أدوات النمذجة الوعائية (VMTK) لاستخراج الخطوط المركزية من الأسطح ثلاثية الأبعاد.
  • التمثيل بالمخطط (Graph Representation): تم تحويل الشبكة الوعائية إلى مخطط شجري موجه وذو جذر $G=(V, E)$.
  • التجزئة (Segmentation): جُمعت الحواف في "قطاعات" (مجموعات من الحواف بين نقاط الالتقاء أو النهايات) لتقليل التعقيد مع الحفاظ على بنية التفرع.
  • تصحيح الأخطاء: تم تنفيذ خوارزميات محددة لتصحيح العيوب الشائعة في التجزئة الآلية:
    • الأوعية التوأم (Twin Vessels): تم اكتشاف وإزالة الحلقات (القطاعات التوأم) عبر دمجها بناءً على عدد الحواف أو نصف القطر.
    • توسيع الخط المركزي: تم استكمال القطاعات الخارجية لتتطابق مع حدود السطح الحقيقية باستخدام استكمال نصف القطر.
    • تنعيم نقاط الالتقاء: تم إزاحة نقاط التفرع نحو الطرف البعيد لإزالة الانثناءات الاصطناعية وضمان أن تكون القطاعات أطول من أنصاف أقطارها.
    • إزالة القطاعات القصيرة: تمت إزالة القطاعات القصيرة بشكل مريب للحفاظ على طوبولوجيا الشجرة.
  • البيانات النهائية: تم الاحتفاظ بـ 44 شبكة صالحة قبل العملية بعد استبعاد 4 شبكات بسبب أخطاء التجزئة.

ب. محاكاة استئصال الرئة

• الاستئصال الافتراضي: لمحاكاة حالات ما بعد الجراحة، قام المؤلفون بإزالة فروع الأوعية الدموية من الـ 44 شبكة بشكل منهجي.
  • إزالة الأوعية غير الجذعية (أي شيء ليس الشريان الرئوي الرئيسي MPA، أو الشريان الرئوي الأيسر LPA، أو الشريان الرئوي الأيمن RPA).
  • إزالة مجموعات من الأوعية الطرفية (الأشقاء والعمات) لمحاكاة عمليات الاستئصال الأكبر.
• النطاق: أنتجت هذه العملية 1,602 شبكة معدلة، مما نتج عنه إجمالي 1,646 عملية محاكاة (44 قبل العملية و1,602 بعد العملية).

ج. النمذجة الهيدروديناميكية

• المعادلات الحاكمة: استُخدم نموذج غير خطي أحادي البعد (1D) يعتمد على معادلات نافييه-ستوكس لتدفق السائل غير القابل للانضغاط في أوعية مرنة متناقصة القطر.
  • تم الحل لإيجاد التدفق ($q$)، ومساحة التجويف ($A$)، والضغط ($p$).
  • العلاقة التكوينية: قانون الجدار المرن الخطي.
• الشروط الحدية:
  • المدخل (Inlet): ملف تعريف معد معدل التدفق (بناءً على بيانات MRI) عند الشريان الرئوي الرئيسي (MPA).
  • المخرج (Outlet): شروط حدية لـ "الشجرة المهيكلة" تمثل الشعيرات الدموية الدقيقة (الأوعية الصغيرة) كشبكة مقاومة. هذا يلتقط المقاومة الطرفية دون الحاجة لنمذجة ملايين الشعيرات الدموية صراحةً.
• المحلل (Solver): طريقة فرق نهائي "لاكس-ويندروف" (Lax-Wendroff) ذات الخطوتين. استمرت عمليات المحاكاة لمدة 20 دورة قلبية حتى الوصول للتقارب (مجموع مربعات الأخطاء < 1).

د. تحليل شدة الموجة (WIA)

• الحساب: تم حساب شدة الموجة ($dI$) باستخدام معادلة "مطرقة الماء"، مع فصل الموجات إلى مكونات أمامية ($+$) وخلفية ($-$).
• التفكيك: تم تفكيك الموجات إلى أربعة مكونات بناءً على تغيرات الضغط ($dp$):
  1. موجة الضغط الأمامية (FCW)
  2. موجة تخفيف الضغط الأمامية (FDCW)
  3. موجة الضغط الخلفية (BCW)
  4. موجة تخفيف الضغط الخلفية (BDCW)
• المعايث (Parameters): تم استخراج ثلاثة معايث لكل مكون:
  • ذروة الشدة ($dI_{peak}$)
  • الوقت للوصول للذروة ($t_{peak}$)
  • شدة الموجة المتكاملة ($I_W$)
• التحقق من الصحة: تم تطبيق مرشح "سافيتزكي-غولاي" (Savitzky-Golay) والتقريب (interpolation) لمطابقة الدقة الزمنية (1 مللي ثانية) وخطوات المعالجة المستخدمة في الدراسة السريرية لـ Glass وآخرون.

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

1. مسار حوسبي مبتكر: تطوير إطار عمل قوي لتحويل الأسطح الرئوية ثلاثية الأبعاد المستمدة من الأشعة المقطعية إلى نطاقات حوسبية أحادية البعد، بما في ذلك تصحيح الأخطاء الآلي لطوبولوجيا المخطط (الأوعية التوأم، وتنعيم نقاط الالتقاء).
2. محاكاة واسعة النطاق: نجحت الدراسة في محاكاة تدفق الدم في 1,646 شبكة شريانية رئوية متميزة (44 قبل العملية و1,602 بعد العملية)، وهو مقياس لم يكن ممكناً الوصول إليه سابقاً باستخدام النماذج ثلاثية الأبعاد الخاصة بكل مريض.
3. العزل الهندسي: أثبتت الدراسة أن التغيرات في شدة الموجة بعد الاستئصال يمكن إرجاعها حصرياً إلى التغيرات الهندسية في الشجرة الشريوية، دون الحاجة لاستدعاء إعادة التشكيل البيولوجي أو التغيرات في خصائص الدم.
4. التحقق مقابل البيانات السريرية: أظهرت الدراسة توافقاً نوعياً قوياً بين تنبؤات النموذج الحوسبي وبيانات الرنين المغناطيسي السريرية لـ 27 مريضاً، مما يؤكد صحة الفرضية القائلة بأن تغيرات الهندسة هي التي تقود تغيرات انعكاس الموجة.

4. النتائج

• الخط المرجعي قبل العملية: أظهر التحليل الإحصائي (اختبارات t المزدوجة) أنه بينما أظهرت الشرايين الرئوية اليسرى (LPA) واليمنى (RPA) سلوكيات تدفق متشابهة قبل العملية، إلا أنها اختلفت بشكل ملحوظ عن الشريان الرئوي الرئيسي (MPA).
• التغيرات بعد العملية:
  • إعادة توزيع التدفق: نجح النموذج في رصد إعادة توزيع التدفق بين LPA وRPA عقب الاستئصال المحاكي.
    // توجات شدة الموجة: اتفقت النتائج الحوسبية مع النتائج السريرية في 29 من أصل 36 تغيراً في المعايث المرصودة (FCW, FDCW, BCW).
    • مثال: في الجانب الذي خضع للعملية، توقع النموذج انخفاضاً في ذروة شدة موجة الضغط الأمامية (FCW) وانزياحاً في توقيتها، مما طابق الملاحظات السريرية.
  • الدلالة الإحصائية: بعد العملية، أصبحت الاختلافات في معايث شدة الموجة بين الجانب الذي خضع للعملية والجانب الآخر ذات دلالة إحصائية ($p < 0.001$)، في حين لم تكن ملحوحة إحصائياً قبل العملية.
• الفرق في المقدار: بينما طابق اتجاه التغيير البيانات السريرية، كان مقدار التغيير في المحاكاة أصغر مما هو عليه في الدراسة السريرية. ويعزو المؤلفون ذلك إلى استخدام ملف تعريف تدفق موحد لجميع المرضى بدلاً من ملفات تعريف خاصة بكل مريض.

5. الأهمية والقيود

الأهمية:

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

القيود:

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

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