المنهجية
تم بناء نموذج نافذة انزلاقية مراقبة للتنبؤ بقيمة PM2.5 المستقبلية عند الزمن t + 1 من خلال الملاحظات السابقة كل ساعة. تم تقييم أطوال نوافذ الإدخال المرشحة التي تبلغ 12 و24 و48 خطوة زمنية كل ساعة باستخدام أداء التحقق، واستخدم التكوين النهائي ل CORTA-Net تسلسل إدخال لمدة 24 ساعة. تم استخدام الاستيفاء الخطي لملء القيم المفقودة؛ تمت إزالة القيم الشاذة باستخدام طريقة IQR؛ وتم تكبير جميع المتغيرات باستخدام التطبيع الأقصى-الأعلى قبل توليد ميزات PM2.5 المتأخرة. ثم تم تطبيق إجراء CorrXGBoost-Rank لاختيار الميزات. أولا، مرشح بيرسون للارتباط يحتفظ بالمتغيرات مع |r| ≥ 0.30. ثانيا، أزواج الميزات المتكررة للغاية ذات ارتباط زوجي |corr(xi, xj)| تم ترشيح ≥ 0.85 لتقليل التعدد التعاوني. ثالثا، تم تدريب مسجل XGBoost على المتغيرات المتبقية، وتم الاحتفاظ بالمتنبئات التي تحمل درجات الأهمية القائمة على كسب XGBoost ≥ 0.015 لإدخال النموذج النهائي. تكونت بنية كورتا-نت النهائية من طبقات LSTM مكدسة للترميز الزمني، وطبقة انتباه متعددة الرؤوس لوزن مستوى الميزات والخطوات الزمنية، وطبقات انحدار كثيفة لتقدير PM2.5 . تم تدريب النموذج باستخدام محسن آدم مع متوسط تربيع فقدان الخطأ والتوقف المبكر. تم إجراء تقييم النموذج باستخدام RMSE و R2 عبر أقسام التدريب والتحقق من الصحة والاختبار والتحقق المتقاطع بعشرة أضعاف. يمثل الجدول 2 ملخص معالجة البيانات المسبقة وهندسة الميزات.
| الخطوة | الطريقة المستخدمة | المعامل / العتبة | الغرض |
| حساب القيمة المفقودة | نسبة الملاحظات المفقودة | التقرير حسب النسبة المئوية المتغيرة | يقييس اكتمال البيانات |
| النسبة القصيرة الفجوة | الاستيفاء الخطي | طول الفجوة ≤ 6 ساعات | تملأ الفترات القصيرة المفقودة |
| اكتشاف القيم الشذرية | طريقة IQR | Q1 − 1.5 × IQR، Q3 + 1.5 × IQR | يزيل القيم القصوى غير الصالحة |
| التطبيع | الحد الأدنى لأقصى القياس | الحد الأدنى والأقصى لمجموعة التدريب | توحيد نطاق الميزات |
| PM₂.₅ lags | المتغيرات المتراجعة | lag₁، lag₂، lag₃ | تلتقط الاستمرارية الزمنية |
| الإحصائيات المتحركة | المتوسطات المتحركة | 3 ساعات، 6 ساعات، 12 ساعة، 24 ساعة | يلتقط التراكم قصير المدى |
| ميزة عدد النيران | موديس فايركونت | عد اليوم نفسه / اليوم السابق | يمثل تأثير الحرائق الإقليمي |
الجدول 2: ملخص معالجة البيانات وهندسة الميزات. في الجدول 2، تتم معالجة البيانات بطريقتين: أولا، من خلال تنظيف وتحويل البيانات الخام (المعالجة المسبقة)، وثانيا، عبر هندسة الميزات لإنشاء أو اختيار أو تعديل الخصائص (الميزات). كوحدة واحدة، تساعد هاتان العمليتان في القضاء على أو تقليل الضوضاء؛ يتعامل مع القيم المفقودة؛ تحسين اتساق البيانات؛ وإنشاء نماذج أكثر تنبؤية.
تم تقدير معلمات التطبيع فقط من مجموعة التدريب ثم تم تطبيقها دون تغيير على مجموعات الاختبار والتحقق لتجنب تسرب المعلومات.
الصياغة الرياضية لرتبة CorrXGBoost
ليكن X = {x1,x 2, ...,x n} يرمز إلى مجموعة المتغيرات المرشحة للمدخلات، وy يرمز إلى تركيز PM₂.₅ المستهدف. لكل ميزة xi، تم حساب معامل ارتباط بيرسون مع المتغير المستهدف كالتالي:
(1)
حيث أن cov(xi, y) هو التغاير بين الميزةx i والهدف y، وσxi و σy هما انحرافات معيارية. الميزات المحفوظة مرضية: |ri| ≥ τr حيث τr = 0.30 في هذه الدراسة.
تم حساب الترابطات الزوجية بين جميع features_xi وx j، وإذا كان |corr(xi, xj)| >= τأحمر، حيث τالأحمر = 0.85، يتم حذف الميزة ذات الارتباط المطلق الأدنى مع هدف PM2.5 من مجموعة الميزات. تقلل هذه العملية من المتغيرات من مجموعة الميزات التي لها تعدد التعاون. بعد ذلك، تم ملاءمة نموذج انحدار XGBoost للميزات المتبقية، وتم حساب درجات أهمية الميزة لكل ميزة باستخدام أهمية XGBoost، وتم الاحتفاظ بميزات تحقق i_i ≥ τxgb حيث τxgb = 0.015 في مجموعة الميزات النهائية (S") المعرفة ك Sالنهائية =S corr ∪ Sxgb، أي مجموعة الميزات المعرفة من خلال تصفية الميزات للتكرار بالنسبة للارتباط وإزالة أي ميزات بناء على أهمية متغير XGBoost. النهج العام المتبع في اختيار الميزات هو كالتالي: حساب ارتباط بيرسون المستهدف على الميزة، التخلص من الميزات التي لها |r_i|< 0.30، التخلص من الميزات التي لها أزواج لديها ارتباطات زوجية ≥ 0.85، ملاءمة XGBoost للميزات المتبقية، الحفاظ على الميزات التي لديها درجة أهمية XGBoost ≥ 0.015، تحديد الميزات النهائية Sالنهائية =S corr ∪ Sxgb، حيث Scorr هي الميزات التي تحتفظ بها بعد تصفية الارتباط المتكررة، وSxgb هي الميزات المختارة باستخدام درجة أهمية ميزات XGBoost. لذا، فإن سير عمل CorrXGBoost-Rank هو: حساب ارتباط بيرسون المستهدف بين ميزات، وإزالة الميزات التي تحتوي على |r i| < 0.30، إزالة الميزات الزائدة للغاية ذات ارتباط زوجي ≥ 0.85، تدريب XGBoost على المتغيرات المتبقية، الاحتفاظ بالمتغيرات ذات درجة الأهمية XGBoost ≥ 0.015، واستخدام اتحاد المتغيرات المختارة بالارتباط والمتغيرات المختارة من XGBoost كمجموعة الميزات النهائية. المعلمات المستخدمة في هذه الدراسة هي τr = 0.30، τالأحمر = 0.85، و τxgb = 0.015.
مصادر البيانات
دمج المؤلفون ثلاث مجموعات بيانات كبيرة للدراسة؛ هذه بيانات عن ملوثات الهواء (PM2.5، PM10، NO2، CO وSO₂) حصل عليها المؤلفون من خلال CPCB (مجلس مكافحة التلوث المركزي) / DPCC (لجنة مكافحة التلوث في دلهي) لمراقبة جودة الهواء، والبيانات الجوية (درجة الحرارة، الرطوبة، سرعة/اتجاه الرياح والضغط) من إدارة الأرصاد الجوية الهندية (IMD)، ومعلومات الأقمار الصناعية عن الحرائق النشطة من منتجات MODIS النشطة للحرائق النشطة من ناسا. تغطي هذه المجموعات الثلاث الفترة من يناير 2012 إلى ديسمبر 2023 والتي تمثل بذلك أنواعا مختلفة من الانبعاثات. تساهم حوادث الحرائق في تلوث الهواء، كما هو موضح في الشكل التكميلي 2، الذي يوضح اتجاهات FIRECOUNT من 2012 إلى 2024.
الاتجاهات طويلة الأمد للملوثات خلال فترة الدراسة
يوضح الشكل 1 الاتجاهات طويلة الأمد للملوثات خلال سنوات الدراسة (2012–2024). تعد أعداد الحرائق السنوية ومتوسط تركيزات PM₂.₅ بين عامي 2012 و2024 علاقة إيجابية قوية إلى حد ما (r = 0.688)، مما يشير إلى أن زيادة نشاط الحرائق عادة ما ترتبط بتركيزات أعلى من PM₂.₅. تتبع بيانات PM₂ والمتوقعة اتجاها مشابها، مما يدعم فكرة أن حرق الكتلة الحيوية يساهم في تلوث الجسيمات. على الرغم من وجود تباين كبير بين السنين، فإن الانبعاثات المرتبطة بالحرائق تعد عاملا مهما في تحديد تباين PM₂.₅، مما يؤكد الحاجة إلى إدارة حرائق إقليمية لتحسين جودة الهواء. الارتباط الذاتي لPM2.5 المتوقع عند ثلاثين تأخيرا، كما هو موضح في الشكل التكميلي 3، له ارتباط إيجابي كبير في جميع التأخيرات الثلاثين، مما يؤكد أن PM2.5 في دلهي له اعتماد زمني قوي واستمرارية متعددة الأيام.

الشكل 1: اتجاهات PM₣ طويلة الأمد وعدد الحرائق من 2012 إلى 2024. يقارن الشكل 1 التغير السنوي في عدد الحرائق مع تركيزات PM₂₅ المرصودة والمتوقعة. يتم الإبلاغ عن PM₂.₅ بوحدة μg/m3. يمثل FIRECOUNT نشاط الحريق المستخرج من الأقمار الصناعية من منتجات MODIS. يرجى الضغط هنا لعرض نسخة أكبر من هذا الشكل.
تحلل STL (التحليل الموسمي والاتجاه باستخدام اللويس)
تقنية STL أو تحليل الموسمية والاتجاهات باستخدام لوس هي خوارزمية تكرارية تسمح بتقسيم بيانات السلاسل الزمنية إلى ثلاثة مكونات إضافية كما هو موضح في الشكل 2: الاتجاه (الاتجاه طويل الأجل) (الشكل 2A)، (الشكل 2D)، (الشكل 2G)، (الشكل 2J)، الموسمي (الفترات الدورية) (الشكل 2B)، (الشكل 2E)، (الشكل 2H)، (الشكل 2K)، والباقي (الضوضاء/البقايا) (الشكل 2C)، (الشكل 2F)، (الشكل 2I)، (الشكل 2L). لقد أثبتت STL نجاحها في التعامل مع الخصائص غير الخطية المعقدة للبيانات البيئية (مثل PM2.5) والتي فشلت العديد من التقنيات التحليلية الأخرى بسبب تغير السعة في الإشارة الموسمية ووجود القيم الشاذة. يسمح تنعيم LOESS بفصل دقيق لهذه المكونات، مما ينتج تفسيرا أوضح للاتجاهات أو الأنماط. من خلال معالجة STL، من الممكن فصل الاتجاهات الصاعدة الأوسع ل PM2.5 عن الدورات اليومية/الموسمية الأقصر وعن البواقي المتقلبة داخل تلك الدورات. يسمح هذا الفصل بتحديد الانبعاثات التي أثرت على PM2.5، على عكس العوامل الجوية التي أثرت على PM2.5. تساعد نتائج هذا الفصل في التنبؤ الدقيق بانبعاثات PM2.5 المستقبلية؛ توفير بيانات للقرارات التنظيمية؛ والامتثال للمعايير المحددة لتفكيك بيانات السلاسل الزمنية الدقيقة لدراسات جودة الهواء34. يشير متوسط قياسات تركيز PM2.5 اليومية التي أجريت في محطات المراقبة في دلهي من 2012 إلى 2024 إلى تغييرات طفيفة في مستوى PM2.5 (أي عدم وجود تغيير كبير من 2022-2024) بشكل عام، وسلوك ثابت محدود للغاية (أو ثبات في مستويات PM2.5 بين المواقع الأربعة للمراقبة). مستويات PM2.5 في قطاعي دواركا 8 وموندكا-DPCC تظهر اتجاها هبوطا مستمرا (أي استمرار في الانخفاض)، بينما لدى أناند فيهار اتجاه تصاعدي (أي في ارتفاع ملحوظ)، وسونيا فيهار لديه اتجاه صعودي طفيف (بعد زيادة طفيفة جدا منذ 2022). بالإضافة إلى ذلك، أشارت بقايا تركيزات PM2.5 اليومية لكل من محطات المراقبة الأربع إلى أن تركيزات PM2.5 تأثرت بشدة بالتغيرات الموسمية (الأرصادية)، مما أدى إلى تقلبات كبيرة في متوسط تركيزات PM2.5 اليومية للمحطات الفردية. يتم عرض تركيزات المراقبة اليومية2.5 مساء لمحطات المراقبة الأربع في دلهي (2012–2024) في الشكل 2.

الشكل 2: تحليل التصنيف الموسمي والاتجاهات باستخدام اللويس (STL) لتركيزات PM₂₅ اليومية عبر أربع محطات مراقبة في دلهي من 2022 إلى 2024. يتم تقسيم السلسلة الزمنية لكل محطة مراقبة إلى ثلاثة مكونات إضافية: الاتجاه طويل الأمد، التغير الموسمي، والباقي (الباقي). (أ) مكون الاتجاه لقطاع دواركا 8. (ب) المكون الموسمي لقطاع دواركا 8. (ج) المكون المتبقي لقطاع دواركا 8. (د) مكون الاتجاه لأناند فيهار. (ه) المكون الموسمي لأناند فيهار. (و) المكون المتبقي لأناند فيهار. (G) مكون الاتجاه ل Mundka-DPCC. (ح) المكون الموسمي لموندكا-DPCC. (I) المكون المتبقي ل Mundka-DPCC. (J) مكون الاتجاه لسونيا فيهار. (K) المكون الموسمي لسونيا فيهار. (L) المكون المتبقي لسونيا فيهار. STL، تحليل المواسم والاتجاهات باستخدام اللويس. يرجى الضغط هنا لعرض نسخة أكبر من هذا الشكل.
يعرض الشكل 2 التحللات الزمنية لتركيزات PM2.5 اليومية في أربعة مواقع مراقبة حضرية خلال 2022-2024، كاشفا عن أنماط من الانخفاضات طويلة الأمد، والزيادات التدريجية، والتغيرات اليومية الكبيرة. يمكن عزو هذه الاختلافات إلى حد كبير إلى التأثيرات الفيزيائية للطقس، أي التغيرات في ارتفاع طبقة حدود الكوكب (PBL)، مثل التمدد النهار الذي يزيد من التشتت الرأسي ويقلل من التركيزات، وانكماش الليل، مما يترك الملوثات قريبة من الأرض ويخلق ذروة أعلى ليلا. تأتي التأثيرات الجوية الأخرى من انبعاثات النشاط البشري المحلي المرتبطة بساعات الحجم الأكبر (ذروات الصباح/المساء)، والصناعات، وتفاعلاتها مع: العوامل الخاصة بالموقع مثل التضاريس المحلية، سرعة الرياح، الرطوبة النسبية، والموسم (أي الارتباط الأكبر في الشتاء) التي قد تكون تسببت في التغيرات الخاصة بالمحطة والانخفاضات طويلة الأمد بشكل عام قد تأثرت بالضوابط التنظيمية على الانبعاثات. مما أدى إلى اتجاهات هبوطية بشكل عام.
الارتباط بين السمات
يعرض الشكل 3 توزيع جميع ميزات الإدخال المستخدمة في نموذج CORTA-Net. في الألواح، تظهر متغيرات الملوثات (PM10 (الشكل 3A)، NO₂ (الشكل 3B)، CO (الشكل 3C)، SO₂ (الشكل 3D)) توزيعات مائلة إلى اليمين نموذجية لبيانات جودة الهواء الحضرية، بينما تظهر O₃ (الشكل 3E) والضغط (الشكل 3I) أنماطا شبه طبيعية35. تظهر درجة الحرارة (الشكل 3F) بنية موسمية ثنائية الشكل واضحة، والرطوبة (الشكل 3G) تتبع انتشارا موحدا واسعا، وسرعة الرياح (الشكل 3H) تظهر توزيعا خفيفا. تسلط هذه الأنماط الضوء على السلوك الإحصائي المتباين للمتنبعات وتبرر الحاجة إلى هندسة الميزات وتطبيره قبل تدريب النماذج.

الشكل 3: توزيع ميزات المدخلات المستخدمة في نموذج CORTA-Net، بما في ذلك تركيزات ملوثات الهواء والمتغيرات الجوية الجوية. توضح التوزيعات الخصائص الإحصائية للمتنبعات قبل المعالجة المسبقة وتدريب النماذج. (أ) تركيز PM₁₀. (ب) تركيز لا يذكر. (ج) تركيز أول أكسيد الكربون. (د) تركيز SO₂. (E) تركيز O₃. (F) درجة حرارة الهواء. (G) الرطوبة النسبية. (H) سرعة الرياح. (1) الضغط الجوي. تظهر المتغيرات الملوثة، وخاصة PM₁₀، NO₂، CO، وSO₂، توزيعات مائلة نحو اليمين نموذجية لبيانات جودة الهواء الحضري، بينما يظهر O₃ وضغط الغلاف الجوي توزيعات طبيعية تقريبا. تظهر درجة الحرارة نمطا موسميا ثنائي الوضع، والرطوبة لها توزيع واسع، وسرعة الرياح تتركز عند قيم أقل مع توزيع ذيل خفيف. تدعم هذه التوزيعات غير المتجانسة استخدام هندسة الميزات وتطبيعه قبل تطوير النموذج. يرجى الضغط هنا لعرض نسخة أكبر من هذا الشكل.
يعرض الشكل 4 تصنيف أهمية الميزات القائم على الكسب XGBoost بعد معالجة مسبق في رتبة CorrXGBoost لتوقع PM₂.₅. يظهر مخطط الأهمية الحالي للميزة أن PM₂.₅ في اليوم السابق هو الأعلى تصنيفا بدرجة أهمية 0.280، يليه PM10 = 0.180، FIRECOUNT = 0.150، NO₂ = 0.120، CO = 0.080، سرعة الرياح = 0.070، درجة الحرارة = 0.040، الرطوبة = 0.030، ويوم السنة = 0.020. يمثل الخط الرأسي المتقطع عتبة اختيار الميزات الفعلية ل XGBoost وهي 0.015. تم الاحتفاظ بالمتنبئات التي كانت درجات الأهمية أكبر أو مساوية من 0.015 لمجموعة المدخلات النهائية ل CORTA-Net، بينما تم استبعاد المتنبئات التي تقل عن هذا العتبة، بما في ذلك SO₂ = 0.010، O₃ = 0.010، الضغط = 0.005، هطول الأمطار = 0.005، عطلة نهاية الأسبوع = 0.002، والعطلة = 0.001. تمثل هذه القيم درجات أهمية الميزات القائمة على مكاسب XGBoost ولا ينبغي تفسيرها كمعاملات ارتباط بيرسون أو تأثيرات سببية.
لذلك، تم تضمين الميزات التي حصلت على مساهمة تجاوزت هذا الحد فقط في النموذج. تستخدم نماذج XGBoost مجموعة من أشجار القرار، التي تبني شجرة بشكل تكراري لتقليل دالة الفقد، وتعرف باسم تعزيز التدرج. يحسب XGBoost أهمية الميزات باستخدام أحد ثلاثة مقاييس: الكسب (مدى تحسن تقسيم الميزة لأداء النموذج)، الوزن (عدد مرات اختيار الميزة كتقسيم لشجرة)، أو التغطية (عدد الملاحظات المتأثرة بالتقسيم). تم بناء نموذج XGBoost على بيانات تتعلق بجودة الهواء (الملوثات، المتغيرات الجوية الجوية) ويركز على PM2.5 المتأخر للتنبؤ التلقائي ب PM2.5، وهو أمر شائع في نماذج PM2.5 في دلهي، ويساعد في التقاط الاستمرارية الزمنية ل PM2.5. Delhi_PM2.5 مساهم مهم بسبب خصائص الارتباط الذاتي العالية للجسيمات الدقيقة، والتي تشير أيضا إلى وجود قصور ذاتي للتلوث من مصادر الانبعاث المتسقة. تعد سرعة الرياح مساهما مهما لأنها توفر آلية لتشتت الملوثات؛ بينما تسمح الرياح المنخفضة بتراكم الملوثات خلال الشتاء عندما تكون هناك انعكاسات. يرتبط NO2 ب PM2.5 بسبب حركة المرور/الانبعاثات، بينما يوم السنة هو مؤشر على دورة الانبعاثات السنوية للملوثات (أي انخفاض الانبعاثات في عطلة نهاية الأسبوع). في دلهي، يتمتع PM2.5 بثبات ذاتي عالي بسبب ركود الأرصاد الجوية الشتوي والانبعاثات المستمرة من المركبات والصناعة والكتلة الحيوية28. تساعد الرياح في تشتت الهباء الجوي، لكن ضعف الرياح < 2 م/ث ساهم في زيادة تراكم PM2.5 بسبب الانعكاسات. العلاقة بين NO2 و PM2.5 هي نتيجة لأصلهما المشترك (أي الاحتراق) وتتعرض لتأثيرات زمنية (مثل اليومية مقابل الأسبوعية). العوامل الأربعة الرئيسية المسؤولة عن تفسير 93٪ أو أكثر من التغير في مؤشر جودة الهواء هي التغيرات في تركيزات الملوثات الناتجة عن تراكم PM2.5 المتأخر (93٪+) وانخفاض سرعات الرياح التي تحبس الانبعاثات، وانبعاثات NO2/PM من المركبات والصناعة، والتغيرات اليومية في كمية حركة المرور37. بالإضافة إلى ذلك، فإن جغرافية دلهي عامل مساهم في استمرار الانعكاسات في المنطقة، وبالتالي تزيد من تركيزات PM2.5 . تعديل المعاملات الفائقة، والتنظيم، وإضافة متغيرات إضافية (مثل درجة الحرارة) سيساعد في تقليل فرص التركيب الزائد وسيحسن الأداء العام للنموذج. يجب أن يشمل تنفيذ استراتيجيات تشغيلية لتقليل انبعاثات PM تنفيذ حدود يومية لانبعاثات PM، واستخدام أجهزة توفر مراقبة فورية لوجود PM، واستخدام الرياح لتوزيع PM عبر المساحات الخضراء الحضرية.

الشكل 4: تصنيف أهمية الميزات القائم على المكاسب في XGBoost بعد فحص ميزات CorrXGBoost-Rank لتوقع PM₂.₅. أشرطة التنبؤ مرتبة بترتيب تنازلي حسب الأهمية. يمثل الخط الرأسي المتقطع عتبة اختيار الميزات الفعلية ل XGBoost وهي 0.015. تم الاحتفاظ بالمتنبئين الذين حصلت على درجات ≥ 0.015 لمجموعة المدخلات النهائية ل CORTA-Net، بينما تم استبعاد المتنبئات التي تقل عن هذا العتبة. يستخدم هذا التصنيف لفحص الميزات وتفسير سلوك النماذج ولا ينبغي تفسيره على أنه نسب سببي. يرجى الضغط هنا لعرض نسخة أكبر من هذا الشكل.
بنية نموذج CORTA-net
يقدم الشكل 5 إطار عمل مقترح للتنبؤ بنظام PM2.5. توضح الأشكال التكميلية 4 و5 الترتيب الداخلي لوحدة LSTM وتمثيل مخطط كتلة الانتباه متعدد الرؤوس. جميع مصادر البيانات (أي الظروف المحيطة، الملاحظات الجوية الجوية، نشاط الحرائق، والجوانب الزمنية والسياقية للبيانات) خضعت لمعالجة مسبقة لضمان توافقها زمنيا، وأن هناك بيانات مفقودة (إذا لزم الأمر)، وأنها خضعت لإزالة الشاذة وتطبيعها قبل استخدامها في أي من عمليات بناء النماذج. يتم عرض بنية النموذج، هندسة الميزات، تكوين التدريب، تكوين البيانات، تعلم النقل، مقاييس التقييم، التنبؤات المستقبلية، وتفاصيل التنفيذ في الجدول التكميلي 1. باستخدام وحدة CorrXGBoost-Rank، يتم تحديد المجموعة المثلى من الميزات قبل إدخال جزء منها في بنية LSTM، والتي تم تعزيزها بدمج طبقات الانتباه متعددة الرؤوس لأخذ سلوك السلاسل الزمنية في الاعتبار خلال مرحلة النمذجة. تم توفير كل من توقعات PM2.5 وميزات PM2.5 لأهمية ميزات الإدخال، بالإضافة إلى أوزان الانتباه متعددة الرؤوس كخرائط انتباه، كمخرجات. تم تضمين هيكل وحدة LSTM الداخلية ومخطط كتلة الانتباه متعدد الرؤوس كشكلين تكميليين 4 و5 على التوالي، بينما يوفر الجدول التكميلي 1 معلمات البنية لكل نوع من الطبقات ضمن هذه البنية.

الشكل 5: البنية العامة لنموذج التنبؤ PM₂.₅ من CORTA-Net. تشمل عملية التنبؤ المدخلات من مؤشرات جودة الهواء، ومؤشرات الطقس، ومؤشرات الوقت، ومؤشرات عداد الحرائق MODIS (مقياس طيفي التصوير متوسط الدقة). لكل من هذه الخطوات، تقوم الخوارزمية بمحاذاة الطوابع الزمنية، وتتعامل مع القيم المفقودة، وتصفية القيم الشاذة، وتطبيع البيانات، وهندسة الميزات. ثم تم اختيار المتنبئين النهائيين باستخدام عملية CorrXGBoost-Rank. بعد ذلك، تم وضع الميزات المختارة في نوافذ منزلقة بتسلسل زمني وتمريرها عبر مشفر LSTM (الذاكرة طويلة المدى القصيرة) باستخدام التعلم التحويلي. لكل ميزة وخطوة زمنية، تطبق آلية الانتباه متعددة الرؤوس أوزانا قبل إرسال المخرج المجمع إلى طبقة الانحدار النهائية الكثيفة لإنتاج توقع PM₂.₅. يرجى الضغط هنا لعرض نسخة أكبر من هذا الشكل.
التدريب والتقييم
تم استخدام نهج نافذة منزلقة لتشكيل تسلسلات تعلم تحت الإشراف. استخدم التدريب محسن آدم ومتوسط تربيع فقدان الخطأ (الجدول التكميلي 2). تم تقييم أداء النموذج باستخدام RMSE وR2 عبر أقسام التدريب والتحقق من الصحة والاختبار والتحقق المتبادل. يتم تجميع بنية النموذج، وهندسة الميزات، وتكوين التدريب، وتكوين البيانات، وتعلم النقل، ومقاييس التقييم، والتنبؤات المستقبلية، وتفاصيل التنفيذ في الجدول التكميلي 1. يمثل الشكل 6 تشخيصات التدريب لنموذج التنبؤ PM₂.₅ من CORTA-Net. في اللوحات، يمثل الشكل 6A منحنيات فقدان التدريب والتحقق، موضحا تقليل الخطأ التدريجي، مع التوقف الأمثل عند العصر 106. يتضح تطور فجوة التعميم أثناء تقييم تحليل الملاءمة الزائدة من خلال التباعد بين بيانات التحقق والتدريب في الشكل 6B. على وجه الخصوص، يظهر فترات يتجاوز فيها التباعد الحدود المحددة مسبقا؛ توفر هذه المعلومات دليلا على أن سلوك معدل التعلم وأنماط تقليل الفقدان توضح فوائد تدهور معدل التعلم المجدول (LR) في الفترات الحرجة كوسيلة لتعزيز استقرار التقارب كما هو موضح في الشكل 6C. تقدم هذه الأشكال مجتمعة ملخصا لديناميكيات التعلم، وقدرة التعميم، وتكوين التدريب الموصى به للنموذج.

الشكل 6: تشخيصات التدريب لنموذج التنبؤ PM₂.₅ من CORTA-Net. (أ) منحنيات فقدان التدريب والتحقق تظهر انخفاضا تدريجيا في الخطأ أثناء تدريب النموذج، مع التوقف المبكر عند العصر 106 بناء على أداء التحقق. (ب) فجوة التعميم بين فقدان التدريب والتحقق عبر العصور، مما يوضح تطور تعميم النماذج وفترات زيادة التباعد التي تشير إلى إمكانية التكيف الزائد. (ج) جدول معدل التعلم الذي يوضح تأثير تدهور معدل التعلم المجدول على التحسين طوال فترة التدريب. (د) ملخص سلوك تقارب النموذج، يعرض التحسين المستقر والتكوين النهائي للتدريب المختار لنموذج كورتا-نت. معا، توضح هذه التشخيصات ديناميكيات النموذج، وخصائص التقارب، وأداء التعميم أثناء التدريب. يرجى الضغط هنا لعرض نسخة أكبر من هذا الشكل.