تقدم هذه المقالة أداة بروتوكول ية وبرامج لتحديد كمية عدم اليقين في معايرة وتحليل البيانات لمحلل كربون عضوي/عنصري حراري شبه مستمر.
وكثيراً ما يسعى باحثون من عدد لا يحصى من المجالات إلى تحديد تركيزات الهباء الجوي الكربوني وتصنيفها على أنها كربون عضوي أو كربون عنصري (EC). ويتم ذلك عادة باستخدام محللات OC/EC الحرارية البصرية (TOAs)، التي تمكن من القياس عن طريق الانحلال الحراري الخاضع للرقابة والأكسدة بموجب بروتوكولات درجة حرارة محددة وضمن الأجواء المقيدة. وتوجد عدة وكالات اتصال تجارية، بما في ذلك أداة شبه مستمرة تتيح إجراء تحليلات على الإنترنت في الميدان. يستخدم هذا الجهاز إجراء معايرة في الاختبار يتطلب معايرة متكررة نسبياً. تفاصيل هذه المقالة بروتوكول معايرة لهذا TOA شبه المستمر ويقدم أداة البرمجيات مفتوحة المصدر لتحليل البيانات والتحديد الكمي الصارم ة مونتي كارلو من عدم اليقين. وتجدر الإشارة إلى أن أداة البرمجيات تتضمن وسائل جديدة لتصحيح الانجراف في الأجهزة وتحديد وقياس عدم اليقين في نقطة الانقسام بين OC/EC. وهذا تحسن كبير في تقدير عدم اليقين في برنامج الشركة المصنعة، الذي يتجاهل عدم اليقين في نقطة الانقسام ويستخدم المعادلات الثابتة للأخطاء النسبية والمطلقة (مما يؤدي عموما إلى عدم اليقين الناقص ة و غالباً ما تسفر عن نتائج غير مادية كما هو موضح في عدة مجموعات بيانات على سبيل المثال). يتم تقاسم بروتوكول المعايرة المثبت وأداة البرمجيات الجديدة التي تمكن من التحديد الكمي الدقيق للشكوك المجمعة من المعايرة والتكرار ونقطة تقسيم OC/EC بهدف مساعدة الباحثين الآخرين في تحقيق أفضل قياسات OC، EC، وكتلة الكربون الكلي في عينات الهباء الجوي.
والقدرة على قياس التركيزات الجوية للأنواع الكربونية بدقة أمر بالغ الأهمية بالنسبة للعديد من الباحثين. وقد اقتُرح أن تكون الأنواع الكربونية في الجسيمات المحيطة (PM، أكبر عامل خطر بيئي للوفاة المبكرة1)العنصر الرئيسي في PM المسؤول عن الآثار الصحية الضارة والنتائج2و3 ،4. والكربون الجسيمي في الغلاف الجوي هو ملوث مناخي حرج، حيث من المعروف أن لمختلف الأنواع الكربونية تأثيرات متغيرة، بل ومعاكسة. الكربون الأسود هو يحتمل أن يكون ثاني أقوى قوة مشعة مباشرة في الغلاف الجوي للأرض5،6،7،8. عند إيداعها على الثلج والجليد، والكربون الأسود يقلل أيضا من انعكاسية المناظر الطبيعية في القطبالشمالي، وتعزيز امتصاص أشعة الشمس، وزيادة معدل ذوبان9،10،11،12 . وعلى النقيض من ذلك، تعمل جزيئات الكربون العضوية التنظيرية كنوى لتكثيف السحب، مما يزيد من انعكاسية الأرض المتوسطة، ويسبب تأثير تبريد13. ومن ثم فإن التصنيف الدقيق للمواد الكربونية المأخوذة عينات والتحديد الكمي المتزامن لأوجه عدم اليقين في القياس هما جانبان أساسيان في قياسات الجسيمات.
ويمكن تحقيق التمييز بين الكربون العضوي والكربون الأولي في عينة محملة بالجسيمات باستخدام تحليل حراري بصري14. تم إنشاء أنظمة تجارية مختبرية لتحليل الكربون الحراري البصري15و16و17 بما في ذلك محلل على الخط وشبه المستمر18 التي تمكن من تنفيذ التحليلات الحرارية البصرية في هذا المجال. ويصف هذا العمل إجراء ً مفصلاً لمعايرة هذا الجهاز الأخير من أجهزة OCEC (انظر جدول المواد)ويشترك في أداة برمجيات مفتوحة المصدر للتقدير الكمي الدقيق لأوجه عدم اليقين في المعايرة والتحليل. وعلى الرغم من أن الإصدار الأولي للبرمجيات المفتوحة المصدر مصمم لشكل ملف إخراج الأداة شبه المستمرة، فإن أداة البرمجيات يمكن أن توسعها الجهات الأخرى بسهولة في المستقبل للعمل مع النواتج التي تولدها أدوات أخرى.
محلل الكربون العضوي/الأولي الحراري البصري شبه المستمر (OCEC) يحدد كمية الكربون العضوي (OC) والكربون الأولي (EC) في حجم العينة. ويتضمن إجراء التحليل أربع مراحل مبينة في الشكل 1. أولا، يتم سحب حجم عينة من خلال الصك، حيث يتم إيداع الجسيمات على، ويتم امتصاص المواد العضوية في مرحلة الغاز من قبل، مرشح الكوارتز. عند إنهاء أخذ العينات ، يتم تسخين فلتر الكوارتز من خلال بروتوكول درجة الحرارة المنصوص عليهفي في الغلاف الجوي الخامل والهيليوم (He). خلال هذا الإجراء، يتم البيروزيلمنت جزء من المواد الكربونية حراريا من مرشح الكوارتز. يتم توجيه العادم الغازي إلى درجة حرارة ثابتة، أكسيد المنغنيز (MnO2)الفرن الذي يحول الأنواع الكربونية البيرالى إلى ثاني أكسيد الكربون (CO2). ويقاس التركيز المحدد زمنياً لثاني أكسيد الكربون المفرج عنه فيما بعد بكاشف الأشعة تحت الحمراء غير المشتتة(NDIR). بعد التدفئة الأولية في بيئة He، يتم تسخين العينة من خلال بروتوكول مماثل في بيئة مؤكسدة (Ox). في وجود الأكسجين، يتم أكسدة الأنواع الكربونية الحرارية المتبقية على مرشح الكوارتز ومن ثم توجيهها من خلال فرن MnO2 وكاشف NDIR بنفس الطريقة. وبمجرد أن يتم تطوير الأنواع الكربونية التي تم أخذ عينات منها بشكل كامل من مرشح الكوارتز، يتم إجراء المعايرة النهائية في الاختبار. يتم إدخال كمية ثابتة (اسميا0.8 مل) من خليط الميثان 5٪ (CH4)-الهليوم في الصك، المؤكسدة في فرن MnO2، وتحويلها إلى CO2،وتقاس في وقت لاحق من قبل NDIR. تتوافق إشارة NDIR المدمجة خلال مرحلة المعايرةداخل الاختبار (التي تسمى CH 4-loop) مع كتلة الكربون المعروفة (التي تم تقديمها باسم CH4)وبالتالي تحدد حساسية NDIR، والتي يمكن أن تنجرف في الوقت المناسب. ثم يستخدم هذا المقياس من حساسية NDIR لاستنتاج كتل الكربون من إشارة NDIR خلال مراحل He- وOx السابقة من التحليل.
الشكل 1: إجراء التحليل الحراري البصري. تحليل إجراء من ال [ثرمب-بصري] [أسيك] جهاز. بعد الحصول على عينة على عامل تصفية الكوارتز (الخطوة 0)، يتم تنفيذ ثلاث خطوات التحليل الرئيسية. يتم تنفيذ بروتوكولين حراريين أولاً في الغلاف الجوي للهيليوم (المرحلة هي، الخطوة 1) ثم في جو مؤكسد (مرحلة الثور، الخطوة 2)، حيث يتم استخدام المكونات الكربونية/أكسدة من فلتر الكوارتز، وأكسدة المزيد فيثاني أكسيد الكربون في غضون فرن تحفيز الثانوية، وتقاس في وقت لاحق من قبل كاشف NDIR. يتم تنفيذ إجراء المعايرةالنهائية في الاختبار (CH 4-حلقة، الخطوة 3)، حيث الأكسدة من كتلة معروفة من الميثان يوفر قدرا من الحساسية للكشف NDIR. وتقترن حساسية الكاشفبثاني أكسيد الكربون المقاس بـ NDIR أثناء مرحلتي “هي” و”أوكس” لتحديد كمية التحميل الكتلي الكربوني على فلتر الكوارتز. الرجاء النقر هنا لعرض نسخة أكبر من هذا الرقم.
يمكن أن تكون كتلة الكربون المحقون أثناء حلقة CH4حساسة لظروف التشغيل بحيث يلزم إجراء معايرة متقطعة. تستخدم هذه المعايرة محلول السكروز المائي من التركيز المعروف (حوالي 0.99٪م / م)كمعيار خارجي. يتم إجراء الاختبارات المتكررة حيث يتم إدخال وحدات تخزين مختلفة معروفة من محلول السكروز بترتيب عشوائي في الأداة ويتم إجراء التحليل الحراري البصري. وينتج عن كل اختبار متكرر (أي نتائج كل حقن وتحليل لاحق) إشارة NDIR متكاملة أثناء حلقة CH4(“منطقة المعايرة”) وإشارة NDIR متكاملة لإجمالي الكربون (أي الإشارة أثناء مرحلتي الهي والثور المشار إليه؛ إلى ك “منطقة إجماليّة”), أيّ يماثل إلى الكتلة معروفة كربون في سكروز. يوفر الانحدار الخطي لكتلة الكربون المعروفة مع “المساحة الإجمالية” المبلغ عنها بالأداة مقياسًا لحساسية NDIR المتوسطة. ثم تقترن هذه الحساسية مع متوسط “منطقة المعايرة” لإنتاج المعرفةمعايرة من كتلة الكربون حقن خلال CH 4-حلقة.
وإلى جانب المعايرة، يتمثل أحد التحديات الرئيسية في تفسير النتائج من أداة OCEC في تحديد الكسور النسبية لـ OC وEC في العينة المقاسة. كما OC pyrolyzes خلال مرحلة He من بروتوكول درجة الحرارة، والأحرف كسر على مرشح بدلا من أن يطلق سراحه، أكسدة في فرن MnO2، وكشف من قبل NDIR. من الناحية النظرية، يبقى هذا OC الحرارية المتفحمة (تسمى الكربون البيرولز، PC) على مرشح حتى مرحلة الثور، عندما يتم أكسدة جنبا إلى جنب مع الجماعة الأوروبية. وبالتالي، فإن وضع علامات ساذجة على جميع الكربون المتطورة خلال المرحلة هي كـ OC والكربون المتطور في مرحلة الثور حيث أن الجماعة الأوروبية تؤدي إلى تقديرات متحيزة للجزء الحقيقي من OC وEC. وهناك وسيلة شائعة لتحديد نقطة الانقسام (أي اللحظة التي يعتبر فيها كل تطور الكربون السابق OC وكل تطور الكربون اللاحق ة كللجماعة الأوروبية) هو الأسلوب19للنفاذ الحراري/البصري ( TOT). هنا، يتم توجيه الليزر من خلال مرشح الكوارتز أثناء التحليل الحراري ويتم الكشف عن قوتها (المصب بصريا من مرشح الكوارتز) من قبل كاشف ضوئي. بافتراض أن OC ليست نشطة بصريا في الطول الموجي الليزر (أي، OC يمتص الضوء negligibly) وPC يشارك الخصائص البصرية مع EC، يمكن تقدير نقطة الانقسام. الفرضية هي قياس قوة الليزر المخففة أولاً في بداية التحليل. كما OC يتطور (جزئيا إلى جهاز الكمبيوتر امتصاص الضوء)، يتم تعزيز التوهين من الليزر بحيث يسقط إشارة الكشف الضوئي. كما يتم إدخال مرحلة الثور وEC / PC تتطور بشكل مشترك، يتم تقليل التوهين، وإشارة كاشف ضوئي يبدأ في الزيادة. يتم تعريف نقطة الانقسام على أنها المثيل في الوقت الذي تعود فيه قوة الليزر المقاسة إلى قيمتها الأولية. وفي حين أن منطق هذا النهج قوي، فإن النتيجة تعتمد على الافتراضات المشار إليها أعلاه. وعلى هذا النحو، من الشائع الإعلان عن أن نتائج OC وEC المبلغ عنها “محددة من الناحية التشغيلية” – أي أنها خاصة بالأسلوب المستخدم لتقييم نقطة الانقسام14و20و21.
على الرغم من أنه من الصحيح من الناحية النظرية أن يتم تطوير OC في مرحلة He وPC / EC تطورت في مرحلة الثور، وقد لوحظ أن تطور PC / EC يمكن أن يحدث في الواقع خلال المرحلة He-بسبب آليات مختلفة22،23،24 ،25، بحيث يمكن أن تحدث نقطة الانقسام الحقيقية قبل إدخال الأكسجين. ويشير هذا الغموض في التنبؤ بالمكان الذي ينبغي أن تكمن فيه نقطة الانقسام، إلى جانب عدم اليقين في افتراضات الخصائص البصرية للأوسي، والكمبيوتر الشخصي، والجماعة الأوروبية، إلى أن عدم اليقين في نقطة الانقسام يمكن أن يكون مصدراً مهيمناً لعدم اليقين في الكربون المقاس الجماهير. ولحسن الحظ، فإن التقدير المنهجي لنقطة الانقسام عن طريق طريقة TOT يمكّن من تقدير موضوعي لعدم اليقين في نقطة الانقسام. غير أنه لا يوجد، حسب علم المؤلفين، تقدير مباشر (وما يترتب على ذلك من انتشار) لعدم اليقين في نقطة الانقسام في برمجيات الشركة المصنعة؛ بدلا ً من ذلك يتم حساب إجمالي حالات عدم اليقين المبلغ عنها مع مكونات نسبية ومطلقة ثابتة26. كجزء من هذا العمل، يتم عرض تقنية جديدة لتقدير عدم اليقين نقطة الانقسام – تقنية “انخفاض التوهين”. وفي هذه التقنية، يُحدد مقدار عدم اليقين في نقطة الانقسام كفرق بين كتلة الكربون المتطورة عند التقسيم الاسمي (عن طريق طريقة TOT) وكتلة الكربون المتطورة في نقطة لاحقة حيث انخفض التوهين بالليزر إلى ما يتجاوز بعض الكمية الحرجة، جزء محدد من التوهين الأولي. ويُقدَّر انخفاض الإنتوهين الحرج استناداً إلى عدم اليقين في التوهين بالليزر بالنسبة لقيمته الأولية؛ من الناحية النظرية، هذا النهج يلتقط عدم اليقين في مبدأ مطابقة التوهين الرئيسي لطريقة TOT. وعلاوة على ذلك، وللنظر (جزئياً على الأقل) في عدم اليقين في نقطة الانقسام بسبب الخصائص البصرية المفترضة للحاسوب الشخصي والجماعة الأوروبية، يتم توسيع الانخفاض المقترح في التوهين الحرج بمعامل اثنين.
تقدم هذه المقالة بروتوكول مفصل لمعايرة أداة OCEC جنبا إلى جنب مع أداة برمجية لتحديد بدقة المعايرة وتحليل الشكوك. أولاً، الأقسام من 1 إلى 3 من البروتوكول تحدد تعليمات إنشاء محلول السكروز المائي، وإعداد الأداة للمعايرة، والحصول على بيانات المعايرة. يستخدم القسم 4 أداة البرنامج الجديدة المفتوحة المصدر (انظر جدول المواد)لتحليل بيانات المعايرة عبر واجهة المستخدم الرسومية للبرنامج. ويحدد القسم 5 اعتبارات الحصول على عينة باستخدام أداة OCEC، ويصف الفرع 6 استخدام البرمجيات المشار إليها أعلاه لحساب كتل الكربون وما يرتبط بها من أوجه عدم اليقين، بما في ذلك المساهمات المستمدة من تقدير التقسيم نقطه. ويرد في وثائق البرنامج على الإنترنت وصف للتقنيات الجديدة لتحسين معالجة بيانات OCEC – بما في ذلك تقنية “انخفاض التوهين” الواردة أعلاه.
ضمن أداة البرامج المعروضة، يتم حساب ثوابت المعايرة وكتل الكربون المقاسة والشكوك المرتبطة بها باستخدام طريقة مونتي كارلو (MC). ينشر هذا الإجراء الأخطاء التي لا، على حد علم المؤلفين، التي تعتبر حاليافي البرامج الخاصة بالشركة المصنعة. بالنسبة للمعايرة، تشمل مصادر الخطأ هذه عدم اليقين في تركيز السكروز في المحلول المائي، والدقة في حجم محلول السكروز التطبيقي (دقة الأداة، وإمكانية التكرار بين المستخدمين، والتكرار داخل المستخدم)، و عدم اليقين في الانحدار الخطي. وفيما يتعلق بتحليل البيانات، تشمل مصادر الخطأ المدروسة عدم اليقين في المعايرة وإمكانية التكرار، والأهم من ذلك، عدم اليقين المقدر لنقطة الانقسام. وفي نهاية المطاف، يمكّن البرنامج المستخدم من تحديد كمية عدم اليقين بدقة في معايرة الأداة ونشر هذا الغموض إلى جانب تقدير النقاط المنقسمة في حساب كتل الكربون. ويمثل ذلك تحسناً ملحوظاً في بروتوكول الشركة المصنعة، وذلك بالنظر مباشرة في مصادر الخطأ الرئيسية في القياس، بدلاً من تقدير أوجه عدم اليقين الجزئية باستخدام معادلة ثابتة.
ويبين الجدول 3 مساهمة المصادر المحددة لعدم اليقين في ثابت المعايرة الكتلية للحالة المثال الموصوفة في الجدول 1 والشكل 2. يتم سرد عدم اليقين المعايرة التراكمية الناجمة عن خطأ التحيز في كاشف NDIR، وخطأ التحيز في تركيز السكروز، وخطأ الدقة والتحيز في وحدات التخزين pipetted. ويميل الخطأ في التحيز في كاشف NDIR (أي التباين في “منطقة المعايرة”) إلى الهيمنة، مع التحيز في إجراء الأنابيب الذي يكون ثاني أهم (على الرغم من صغره في المثال التمثيلي). وبالتالي فإن التقدير السليم لخطأ الأنابيب أمر بالغ الأهمية لضمان التحديد الكمي الدقيق لعدم اليقين في المعايرة عموماً؛ وبالإشارة إلى الخطوة 4-3-2 من البروتوكول، يُقترح بالتالي تقييم قابلية التكرار داخل المستخدم وإمكانية استنساخه فيما بين المستخدمين لكل مجموعة من المستعملين والماصة. وعلى النقيض من ذلك، فإن عدم اليقين الناجم عن تركيز السكروز في المعيار الخارجي صغير للغاية. وعلاوة على ذلك، يبدو أن هناك مساهمة ضئيلة من عدم اليقين الانحدار، نتيجة محتملة للخطية جيدة من الصك – معامل التحديد (R2)لاحتواء خطي من بيانات المعايرة يتجاوز عموما 99.95٪ . إذا لم تكن بيانات المعايرة خطية بما فيه الكفاية، يقوم البرنامج تلقائياً بتحذير المستخدم، الذي يكون قادراً بعد ذلك على استكشاف أخطاء مجموعة البيانات وإصلاحها عبر عنصر التحكم المنطقي المشار إليه في الخطوة 4.2 من البروتوكول؛ ثم يمكن للمستخدم تعديل مجموعة بيانات المعايرة الخاصة بهم عن طريق الحصول على بيانات الاستبدال حسب الحاجة.
أوجه عدم اليقين التي تم النظر فيها | معايرة كاملة (6 نقاط) |
انحياز NDIR | ± 2.61٪ |
+ محلول السكروز | ± 2.61٪ |
+ ماصة | ± 2.78٪ |
النتيجة الاسمية [ميكروغرامC] | 18.49 |
الجدول 3: أوجه عدم اليقين في ثبات المعايرة الجماعية. مساهمة تمثيلية من عدم اليقين في معايرة أداة OCEC للمعايرة من ست نقاط على سبيل المثال (انظر الجدول 1). ويهيمن على عدم اليقين المعايرة عموما من قبل التحيز في كاشف NDIR مع خطأ بسبب الدقة في الأنابيب من محلول السكروز، بما في ذلك الأخطاء البشرية (استنساخ بين المستخدمين والتكرار داخل المستخدم)، ويجري ثاني أهم، يليها عدم اليقين في الانحدار الخطي وتركيز السكروز (وكلاهما لا يكاد يذكر).
معايرة جهاز OCEC هو إجراء يستغرق وقتاطويلا، وعادة ما يتطلب 2-3 ساعات لإكمال اعتمادا على طول البروتوكول الحراري المستخدمة. من المستحسن إجراء معايرة أسرع. ولهذه الغاية، تم تحليل فعالية بروتوكول معايرة معدل واقتطاع باستخدام أداة البرامج المعروضة. تم تنفيذ إجراء MC المطور باستخدام جميع المجموعات الفرعية الممكنة من بيانات معايرة المثال المدرجة في الجدول 1 — تقتصر على الحالات التي تحتوي على ثلاثة بيانات أو أكثر وقياس فارغ واحد على الأقل. يتم رسم كافة ثوابت المعايرة الكتلية الناتجة من هذا التحليل في الشكل 4 كدالة لعدد بيانات المعايرة المستخدمة، حيث تم تطبيع ثوابت المعايرة بواسطة نتيجة المعايرة الكاملة (6 نقاط). ومن غير المستغرب أن يزداد عدم اليقين في ثبات المعايرة مع تقليل بيانات المعايرة المتوفرة. ومع ذلك، فإن وسائل جميع المعايرة المقتطعة تقع ضمن CI 2י من نتيجة المعايرة الكاملة، والتي هي نتيجة للخطية المذكورة أعلاه من الصك. يشير هذا الاتساق في المتوسط MC إلى أنه يمكن استخدام معايرة معدلة وأسرع تتكون من بيانات معايرة قليلة كفحص “اختبار عثرة” لمعايرة أداة OCEC. وهذا يعني أنه إذا كان متوسط MC لمجموعة بيانات معايرة من 3 نقاط ضمن CI 2י للمعايرة الحالية، فمن المحتمل أن أداة OCEC لا تتطلب إعادة المعايرة. ومن الواضح أيضاً في الشكل 4 أن عدم اليقين في المعايرة ينخفض مع المزيد من بيانات المعايرة، ولكن الانخفاض في عدم اليقين يعاني من تناقص العائدات. بالإشارة إلى الجدول 3 ومناقشته أعلاه، حيث أن عدم اليقين في المعايرة يهيمن عليه انحياز NDIR (كمياً مع الخطأ القياسي في “مناطق المعايرة”)، فإن الانخفاض الهامشي في عدم اليقين في المعايرة عن طريق إدراج الـ nth يمكن تقدير نقطة البيانات مع العامل √(1−1/n). وبالتالي، فإن الانخفاض الهامشي في عدم اليقين، في المثال التمثيلي، أكبر عند الانتقال من معايرة من ثلاث إلى أربع نقاط منه عند الانتقال من معايرة من خمس إلى ست نقاط. وتتيح أداة البرمجيات المطورة، التي يمكن تنفيذها بعد الحصول على كل نقطة بيانات معايرة (أي بعد كل تكرار للخطوة 3-3 من البروتوكول)، للمستخدم تحديد مدى عدم اليقين في المعايرة طوال عملية الحصول على البيانات. ومن الأهمية بمكان أن هذه القدرة تمكن المستخدم ليس فقط من اتخاذ قرار بشأن التقارب الكافي للمعايرة في سياق عدم اليقين، ولكن أيضا الكشف عن وجود بيانات زائفة – وهذا هو الانخفاض المتزايد في عدم اليقين المعايرة التي لا سيما مختلفة من المتوقع يسلط الضوء على المستخدم أن أحدث نقطة بيانات المعايرة المكتسبة يمكن أن تكون معيبة.
الشكل 4: تقييم حجم العينة على عدم يقين المعايرة. ثابت المعايرة الكتلية المحسوبة لجميع المجموعات الممكنة من بيانات المعايرة المدرجة في الجدول 1 (تتطلب ما لا يقل عن ثلاثة بيانات، بما في ذلك واحد على الأقل فارغ) تم تطبيعها بالنتيجة باستخدام جميع البيانات الستة. ينخفض عدم اليقين النسبي في ثابت المعايرة مع زيادة في عدد بيانات المعايرة. وتتطابق المنطقة المظللة الزرقاء في الشكل مع CI 2י من ثابت المعايرة المحسوب باستخدام جميع بيانات المعايرة. ومن الواضح أن جميع النتائج الاسمية تقع ضمن هذا CI مما يشير إلى أنه – على الرغم من عدم اليقين – يمكن استخدام إجراء معايرة مقتطع من ثلاث نقاط بيانات معايرة فقط كفحص “اختبار عثرة” لمعايرة الأداة. الرجاء النقر هنا لعرض نسخة أكبر من هذا الرقم.
وترد في الجدول 4تفاصيل كتل الكربون المحسوبة وأوجه عدم اليقين فيما بين مجموعات البيانات الأربعة على سبيل المثال؛ ويتم الحصول على هذه البيانات من قياسات الانبعاثات الكربونية من اللهب الهيدروكربوني27و29وتوربينات الغاز29والشكل الدقيق (< 2 μm) الجسيمات الكربونية التي تم الحصول عليها من عينات الرواسب تحت المحيطات30. وترد كتل OC وEC وTC المحسوبة في الجدول 4 إلى جانب انخفاض التوهين الحرج المحسوب ونسبة EC/TC لكل مجموعة بيانات، مما يبين اتساع نطاق البيانات النموذجية في سياق تكوين الكربون المعون. ويرد موجز في الجدول 4 للمصادر وأوجه عدم اليقين العامة في كتلة الكربون باستخدام أداة البرمجيات المعروضة مقارنة بالمصادر التي أبلغ عنها الصك. ضمن أداة البرنامج، ينتج عن الجمع (توزيع المنتج) من عدم اليقين والمعايرة عدم اليقين الكلي في كتلة TC المكلورة MC (الكمية كـ −8.32/+8.40% في العمل الحالي)، وهو أمر مستقل عن عدم اليقين في انقسام نقطة وبالتالي بمثابة الحد الأدنى في عدم اليقين من OC والجماهير EC. وتطبق هذه الشكوك في المعايرة التمثيلية على كل مجموعة بيانات على سبيل المثال مع الحفاظ على ثابت المعايرة الكتلية الاسمية المستخدم في التحليلات الأصلية.
الجدول 4: أوجه عدم اليقين في تحليل البيانات. المساهمة في عدم اليقين في كتل الكربون التي تقيسها OCEC لأربعة مجموعات بيانات على سبيل المثال من مجموعة واسعة من المصادر وتؤديهامختبرات مختلفة27و29و30. (أ)النتائج العددية الرئيسية من مجموعات البيانات المثال: OC، EC، وTC الجماهير، وانخفاض التوهين الحرجة للقياس الكمي من عدم اليقين نقطة الانقسام، ونسبة الكربون من عنصري إلى إجمالي. (ب)موجز للشكوك في كتل الكربون المحسوبة. وتشمل مصادر عدم اليقين المساهمة ثابت المعايرة الكتلية، وتكرار إجراء المعايرة، وعدم اليقين في نقطة الانقسام (بالنسبة لكتلة TC) المقابلة لانخفاض التوهين الحرج المدرجة في (أ). وترد أيضاً في (ب) أوجه عدم اليقين في كتل الكربون باستخدام المعادلات الثابتة (مكافئ (1)) التي تستخدمها الأداة. ويتوافق تسليط الضوء الأحمر والأصفر على البيانات مع التقدير الناقص والإفراطي لعدم اليقين، على التوالي، عند استخدام المعادلات الثابتة بالنسبة للطريقة الحالية. وفي معظم الحالات، فإن الصك لا يقدر قدرة قدرة الكربون على الكتلة، إلا أنه إذا كانت كتلة الوسواس القهري أو الجماعة الأوروبية المقاسة صغيرة، فإن الصك قد يبالغ في تقدير عدم اليقين مقارنة بالبرنامج الحالي. الرجاء النقر هنا لتحميل هذا الملف.
وبالنسبة لهذه الأمثلة، استخدمت تقنية انخفاض التوهين لتحديد مقدار عدم اليقين في نقطة الانقسام. وتراوحتالقيم الحرجة المحسوبة لتراجع التوهين من 1.342% إلى 2.059% مما أدى إلى حالات عدم يقين من نقطة الانقسام من 0.10% إلى 4.50% من كتلة TC. وفي حين أن القيم المستخدمة لانخفاض التوهين هي في الواقع ذاتية إلى حد ما – ولا سيما العامل المستخدم لتقدير عدم اليقين في نقطة الانقسام بسبب نقطة الخصائص البصرية – فإن هذه الأمثلة تبرز اعتماد عدم اليقين في نقطة الانقسام على بيانات التحليل المحددة. فعلى سبيل المثال، يكون عدم اليقين عند نقطة الانقسام حساساً لميل قطعة الرسم التي تحتوي عليها أفيا بالقرب من نقطة التقسيم الاسمية. النظر في مجموعة البيانات “A” المقابلة للبيانات المثال في الشكل 3 ومجموعة البيانات “D”؛ على الرغم من وجود انخفاضات مماثلة في التوهين الحرج، فإن المنحدرات الضحلة والحادة نسبياً لأراضي AVEC الخاصة بكل منها (على سبيل المثال، انظر الشكل 3ب لمجموعة البيانات “A”) تسفر عن أكبر وأصغر نقطة انقسام الشكوك بنسبة 4.50٪ و 0.10٪ من TC الجماعي. وبالإضافة إلى ذلك، تبين البيانات النموذجية أن تأثير عدم اليقين في نقطة الانقسام يعتمد إلى حد كبير على نطاقه بالنسبة لكتل OC وEC الاسمية. النظر في مجموعات البيانات على سبيل المثال “باء” و “جيم”، التي لديها شبه متطابقة نقطة الانقسام الشكوك (1.22٪ من كتلة TC)؛ تحتوي مجموعة البيانات “C” على 43% OC بينما تحتوي مجموعة البيانات “B” على 8%؛ ويؤدي انخفاض الكمية النسبية لـ OC في هذه الأخيرة إلى مضاعفة عدم اليقين الشامل تقريباً. وتبرز هذه النتائج بشكل حاسم شرط النظر مباشرة في عدم اليقين في نقاط الانقسام في سياق بيانات الـ AVEC للتحليل وكتل الكربون الإجمالية.
ويبين الجدول 4أيضاً أوجه عدم اليقين التي أبلغ عنها الصانعون في كتل الكربون. ولا تراعي هذه التقديرات مباشرة أوجه عدم اليقين في المعايرة ونقطة الانقسام، بل تُحسب باستخدام العلاقات الثابتة المبينة في التاق (1)26،حيث تمثل مط الكتلة الاسمية لمكون الكربون المحدد.
(1)
وتسمح هذه العلاقات الثابتة بأن تكون أوجه عدم اليقين المقدرة في كتلة OC و/أو EC أقل من كتلة التعاون التقني بشكل مصطنع – ويحدث هذا الشرط عندما تكون كتلة OC أو EC أقل من ثلث كتلة التعاون التقني، كما هو الحال بالنسبة لمجموعات البيانات “ألف” و”باء” و”دال”. هذه النتيجة غير مادية، لأن عدم اليقين النسبي في كتل OC وEC يجب أن يحدها أدناه من قبل كتلة TC، نتيجة لنشر عدم اليقين نقطة الانقسام في كتل OC وEC المحسوبة. تتوافق الخلايا التي تم تمييزها باللونين الأحمر والأصفر في الجدول مع التقديرات الناقصة والمفرطة لعدم اليقين في كتلة الكربون عند استخدام Eq. (1) من الشركة المصنعة. المعادلات الثابتة أقل من تقدير عدم اليقين الشامل TC لجميع الأمثلة الأربعة، نتيجة لكتل TC المحسوبة كبيرة بما فيه الكفاية. وفي معظم الحالات، كانت المعادلات الثابتة أيضاً أقل من تقدير عدم اليقين الكتلي للجماعة الأوروبية وOC، باستثناء الحالات التي كانت فيها OC (مجموعة البيانات “باء”) والجماعة الأوروبية (مجموعة البيانات “دال”) صغيرة بما يكفي للتسبب في الإفراط في التقدير عن طريق مكافئ (1). هذه الزيادة في عدم اليقين المتقاربة عن طريق Eq. (1) يتفق مع البرنامج الحالي منذ عدم اليقين في كتلة OC وEC بسبب أن من نقطة الانقسام تعتمد على حجمها المطلق; ومع ذلك، فإن أوجه عدم اليقين عبر المعادلات الثابتة لا تتبع مع تلك الخاصة بالبرمجيات الحالية، التي تنظر مباشرة في حالة عدم اليقين في نقطة الانقسام وتنشرها في سياق بيانات تحليل محددة.
إن استخدام إطار MC في أداة البرمجيات المعروضة ضروري لنشر أوجه عدم اليقين في المكونات بدقة من خلال الخوارزميات غير الخطية لتحليلات OC/EC الحرارية البصرية. ومع ذلك، من المهم ملاحظة أن أساليب MC، من خلال طبيعتها العشوائية بطبيعتها، ليست حتمية وتميل إلى أن تسفر عن نتائج غير متسقة إذا كان العدد المستخدم من السحوبات/التكرارات MC (انظر الخطوتين 4-3-3-3 و6-5-2) غير كاف – على غرار ما يلي: عينة صغيرة الحجم إحصائيا. ولذلك، هناك اتساق متأصل مقابل حساب الوقت المفاضلة التي ينبغي النظر فيها عند معالجة البيانات باستخدام OCECgo. وبالتالي فمن المستحسن للمستخدمين لإجراء المعالجة الأولية واستكشاف الأخطاء وإصلاحها من البيانات باستخدام عدد صغير من السحوبات MC (على سبيل المثال، 104). بمجرد أن تسفر الحسابات عن نتائج مرضية، يجب علىالمستخدم ثم زيادة عدد السحوبات MC (إلى 10 6-108)لتحقيق نتيجة أقل تأثراً بالطبيعة المنفصلة والعشوائية لطريقة MC. وإلى جانب “التعريف التشغيلي” اللازم لتحليلات OC/EC الحرارية البصرية، هناك قيود أخرى في معالجة هذه البيانات يجب الاعتراف بها أيضاً عند استخدام OCECgo لحساب بيانات OC/EC والإبلاغ عنها. أولاً، الصكوك المستندة إلى NDIR (مثل تلك المدرجة في جدول المواد)تعاني من الانجراف في إشارة NDIR التي يجب تصحيحها. وفي هذا البروتوكول (انظر الخطوة 6-3-2 ووثائق OCECgo)، يمكن للمستخدم أن يستخدم نهجاً جديداً لتصحيح الانجراف في كاشف NDIR اختيارياً. في حين أن تجربة المؤلفين هذا يؤدي إلى نتائج محسنة على تصحيح NDIR الخطي القياسي للشركة المصنعة، يجب أن يلاحظ أن عدم اليقين في هذا التصحيح NDIR هو التحدي إن لم يكن من المستحيل تحديد كمي، وبالتالي لا يزال غير معروف عنصر من عدم اليقين في حساب كتل الكربون. وعلى نفس المنوال، من الصعب أيضا تحديد كمية عدم اليقين في الافتراض اللازم بأن الخصائص البصرية للحاسوب وEC تشتركان في خصائص بصرية. إذا تم اختيارها (انظر الخطوة البروتوكولية 6-4-1)، تحاول تقنية انخفاض التوهين الحرج أن تحد بتحفظ من تأثير هذا الافتراض من خلال عامل توسع ذاتي. بيد أن الأهم من ذلك هو أن هذا ليس بالضرورة سوى تقدير، ويقترح على المستخدم تقييم أثر عامل التوسع هذا (أي انخفاض التوهين الحرج) على بياناته المحددة. وتقدم OCECgo كأداة مفتوحة المصدر بحيث يمكن بسهولة توسيعنطاقها من قبل المؤلفين وغيرهم من المتعاونين المهتمين ليس فقط للتواصل مع الصكوك الأخرى ولكن أيضا لتشمل وظائف أخرى مفيدة، ومحددة المجال. وعموماً، فإن أداة البرمجيات المفتوحة المصدر المطورة إلى جانب إجراء المعايرة التفصيلية تهدف إلى المساعدة في تحقيق قياسات أكثر دقة لكتلة OC وEC وTC في عينات الهباء الجوي مع تبسيط الحساب القوي للقياس الشكوك.
The authors have nothing to disclose.
وقد تم دعم هذا العمل من قبل مجلس العلوم الطبيعية والبحوث الهندسية في كندا (NSERC) الشبكة الاستراتيجية FlareNet (منحة # 479641)، منح البحوث اكتشاف NSERC (منحة #06632 و 522658)، والموارد الطبيعية كندا (مدير المشروع، مايكل طبقة). ويعرب المؤلفون عن امتنانهم لأولئك الذين تقاسموا ملفات البيانات الأولية لاستخدامها كأمثلة تمثيلية في هذا العمل.
10% oxygen gas in helium | Local gas supplier | – – – | Primary or certified standard preferred |
5% methane gas in helium | Local gas supplier | – – – | Primary or certified standard preferred |
Distilled, de-ionized water | Harleco | 6442-85 | ASTM D1193-91 Type II or Type I (preferred) |
Filter punch tool | Sunset Laboratories Inc. | – – – | Included with carbon analyzer |
Filter removal tool | Sunset Laboratories Inc. | – – – | Included with carbon analyzer |
Glass jar (4 oz.) | ULINE | S-17982P-BL | Or suitable equivalent; borosillicate glass preferred |
Helium gas | Local gas supplier | – – – | Ultra-high purity (> 99.999%) or better preferred |
High-accuracy thermal gas mass flow meter | Bronkhorst | EL-FLOW Prestige | For accurate measurement of sample volume (see Protocol step 5) |
High-purity sucrose | Sigma Aldrich | S9378 | Purity ≥ 99%m/m or higher preferred |
Lint-free tissues | Kimtech | 34155 | Or suitable equivalent |
MatLab Runtime (R2016a or newer) | MathWorks Inc. | mathworks.com | Search "runtime compiler" and install appropriate version for the operating system |
Non-powdered, disposable, plastic gloves | VWR | 89428-752 | Or suitable, properly-sized equivalent |
OCECgo software | Carleton University, Energy and Emissions Research Lab. | GitHub Repository | Source and build distributions of the software are available on GitHub |
Oxygen trap | Supelco | 22449 | Or suitable GC-quality equivalent |
Pipette | Eppendorf | 3120000020 | Model: Research® Plus 0.5 – 10 μL – Or any single-channel, adjustable volume, manual pipette |
Pipette tips | Eppendorf | 022492012 | Model: epT.I.P.S.® Standard, 0.1 – 20 μL |
Precision scale / balance | AND | FX-3000IWP | Precision balance with capacity > 1 kg |
Quartz filters | Pall | 7202 | Model: Tissuquartz 2500 QAT-UP – 47 mm |
Semi-continuous thermal-optical organic/elemental carbon analyzer | Sunset Laboratories Inc. | – – – | Model 4 semi-continuous analyzer |