محاسبه عدم‌قطعیت اندازه‌گیری با روش‌های LPU و Monte Carlo؛ راهنمای آزمایشگاه‌های ISO/IEC 17025

راهنمای عملی انتخاب روش ارزیابی عدم‌قطعیت اندازه‌گیری با مقایسه کاربردها، محدودیت‌ها و معیارهای انتخاب LPU و شبیه‌سازی مونت‌کارلو

این مقاله مخصوص مسئول فنی، مدیر کیفیت، کارشناس فنی و هر آزمایشگاهی که می‌خواهد بودجه عدم‌قطعیت و نتیجه انطباقش در ارزیابی قابل دفاع باشد.

LPU یا Monte Carlo؛ مسئله واقعی انتخاب «روش» نیست، انتخاب «مدل قابل دفاع» است

آزمایشگاهی که بین قانون انتشار عدم‌قطعیت (LPU) و روش مونت‌کارلو (MCM) یکی را صرفاً چون ساده‌تر است انتخاب می‌کند، هنوز مسئله اصلی را حل نکرده است. روش درست، روشی است که با مدل اندازه‌گیری، توزیع ورودی‌ها، سطح غیرخطی‌بودن، هدف گزارش و ریسک تصمیم انطباق سازگار باشد.

در ISO/IEC 17025، ارزیابی عدم‌قطعیت یک خروجی تزئینی برای گواهی کالیبراسیون نیست؛ بخشی از شواهد اعتبار نتیجه است. برای کالیبراسیون، آزمایشگاه باید عدم‌قطعیت را ارزشیابی کند؛ برای آزمون نیز هرجا برای اعتبار نتایج، روش، گزارش یا تصمیم انطباق لازم باشد، باید آن را به‌صورت متناسب در نظر بگیرد. استاندارد یک نرم‌افزار یا یک فرمول واحد تجویز نمی‌کند؛ اما انتظار دارد مدل، داده‌ها، مفروضات، محاسبات و بازبینی فنی قابل ردیابی باشند.

قاعده طلایی: Monte Carlo ذاتاً «دقیق‌تر» از LPU نیست و LPU هم ذاتاً «کافی» نیست. هر دو فقط به اندازه مدل اندازه‌گیری، داده ورودی و کنترل مفروضات معتبرند.

LPU و Monte Carlo در یک نگاه

موضوعLPU؛ قانون انتشار عدم‌قطعیتMonte Carlo؛ انتشار توزیع‌ها
منطقاثر عدم‌قطعیت ورودی‌ها با تقریب موضعیِ مدل و ضرایب حساسیت به خروجی منتقل می‌شود.از هر توزیع ورودی نمونه‌گیری و مدل بارها اجرا می‌شود تا توزیع خروجی ساخته شود.
تناسب معمولمدل خطی یا نزدیک به خطی، عدم‌قطعیت‌های نسبتاً کوچک و خروجی تقریباً متقارن.مدل غیرخطی، ورودی‌های کراندار یا غیرنرمال، خروجی نامتقارن یا تصمیم نزدیک حد پذیرش.
مزیت عملیشفاف، سریع، مناسب روش‌های روتین و قابل ممیزی در کاربرگ.نمایش مستقیم شکل توزیع خروجی و فاصله پوشش نامتقارن در مدل‌های پیچیده.
ریسک اصلینادیده‌گرفتن غیرخطی‌بودن، هم‌بستگی یا شکل واقعی توزیع‌ها.شبیه‌سازیِ مفروضات اشتباه، تکرار ناکافی یا نبود کنترل عددی و بازتولیدپذیری.
آنچه باید ثبت شودمدل، ضرایب حساسیت، u هر ورودی، کوواریانس‌های مهم، u_c و ضریب پوشش.مدل، توزیع/پارامتر هر ورودی، وابستگی‌ها، تعداد تکرار، بذر/نسخه محاسبه و آزمون کفایت عددی.

منبع روش‌شناسی: JCGM 100 قواعد عمومی ارزشیابی و بیان عدم‌قطعیت را ارائه می‌کند و JCGM 101 انتشار توزیع‌ها از طریق مدل و پیاده‌سازی آن با Monte Carlo را شرح می‌دهد. JCGM GUM‑6 نیز توسعه و استفاده از مدل‌های اندازه‌گیری را محور قرار می‌دهد.

جدول انتخاب: چه زمانی LPU و چه زمانی Monte Carlo؟

در این جدول «ترجیح» به معنی حکم خودکار نیست. اگر چند علامت از ستون Monte Carlo هم‌زمان دیده شود، ابتدا مدل و داده‌ها را بازبینی کنید و سپس MCM را به‌عنوان روش اصلی یا کنترل مستقل اجرا کنید.

معیار تشخیصLPU معمولاً مناسب‌تر استMonte Carlo را جدی بررسی کنیداقدام مسئول فنی
رفتار مدلمدل جمعی/خطی است یا در بازه عدم‌قطعیت تقریب خطی پاسخ قابل قبولی می‌دهد.توان، نسبت، لگاریتم، آستانه، اشباع یا تبدیل غیرخطی اثر محسوس دارد.یک نمودار یا آزمون حساسیت در بازه ورودی بسازید؛ «خطی به نظر می‌رسد» کافی نیست.
ورودی‌هاتوزیع‌ها متقارن و پارامترهایشان قابل دفاع‌اند.توزیع‌ها کراندار، چوله، چندمدی یا حاصل داده کم‌اند.برای هر ورودی، منبع داده و دلیل انتخاب توزیع را ثبت کنید.
خروجیخروجی تقریباً نرمال/متقارن و فاصله پوشش دوطرفه منطقی است.فاصله پوشش خروجی نامتقارن است یا کران فیزیکی دارد.به‌جای اجبار U = k×u_c، صدک‌های پوشش را گزارش کنید.
نزدیکی به حدنتیجه فاصله معناداری از حد پذیرش دارد و تصمیم وابسته به دم توزیع نیست.نتیجه نزدیک حد است یا اظهار انطباق برای مشتری/مرجع صادر می‌شود.قاعده تصمیم و ریسک مصرف‌کننده/تولیدکننده را پیش از صدور نتیجه تعیین کنید.
وابستگی ورودی‌هاورودی‌ها مستقل‌اند یا هم‌بستگی محدود و قابل محاسبه است.وابستگی پیچیده یا مشترک میان تصحیحات، ضرایب و داده‌های برازش وجود دارد.منشأ وابستگی را در مدل نشان دهید؛ حذف آن بدون توجیه ممنوع است.
کنترل اجراکاربرگ LPU نسخه‌دار و بازبینی‌شده در دسترس است.کد/نرم‌افزار MCM نسخه‌دار، آزمون‌شده و قابل تکرار است.قبل از اتکا به خروجی، ورودی نمونه و نتیجه مرجع را مستقل بازبینی کنید.

درخت تصمیم عملی برای انتخاب روش

  1. کمیت مورد اندازه‌گیری (measurand) را با واحد، شرایط، روش و هدف استفاده از نتیجه تعریف کنید.
  2. مدل اندازه‌گیری را بنویسید: Y = f(X₁, X₂, …, Xₙ). هر تصحیح مؤثر، عامل رقیق‌سازی، مرجع، تکرارپذیری و اثر محیطی باید جای مشخصی در مدل داشته باشد.
  3. برای هر ورودی، داده، عدم‌قطعیت استاندارد، توزیع و هر وابستگی مهم را تعیین کنید.
  4. آزمون کنید: آیا مدل در بازه موردنظر تقریباً خطی است و خروجی انتظار می‌رود متقارن باشد؟ اگر بله، LPU را محاسبه کنید.
  5. اگر پاسخ منفی است، یا نتیجه به حد پذیرش نزدیک است، MCM را اجرا و پایداری صدک‌های خروجی را کنترل کنید.
  6. در مدل‌های حساس، LPU و MCM را کنار هم مقایسه کنید؛ اختلاف معنادار، نشانه بازبینی مدل یا ترجیح MCM است، نه خطای خودکار یکی از روش‌ها.
  7. روش انتخاب‌شده، نسخه محاسبه، نتیجه عدم‌قطعیت و در صورت لزوم قاعده تصمیم را در سوابق روش ثبت و در بازنگری‌های فنی به‌روز کنید.

خطاهای مدل‌سازی که نتیجه را غیرقابل‌دفاع می‌کنند

بخش بزرگی از اختلاف در عدم‌قطعیت از محاسبه به‌وجود نمی‌آید؛ از مدل ناقص یا مفروضه پنهان به‌وجود می‌آید. جدول زیر خطاهایی را نشان می‌دهد که در بازبینی فنی باید فعالانه جست‌وجو شوند.

خطانشانه در پروندهچرا خطرناک استکنترل پیشنهادی
تعریف مبهم measurandعنوان روش کلی است، اما کمیت، ماتریس، شرایط یا واحد دقیق نیست.ممکن است همان بودجه برای دو اندازه‌گیری متفاوت استفاده شود.یک جمله تعریف فنی بنویسید: چه چیزی، در چه شرایطی و با کدام روش اندازه‌گیری می‌شود؟
حذف تصحیح مؤثرگواهی کالیبراسیون، بلانک، بازیابی یا دما در گزارش هست اما در مدل نیست.سوگیری یا پراکندگی واقعی از خروجی حذف می‌شود.هر آیتم گزارش را با هر ترم مدل تطبیق دهید؛ «صفر بودن تصحیح» هم باید مستند باشد.
دوبارشماریتکرارپذیری روش و تکرارپذیری دستگاه بدون تفکیک منبع وارد شده‌اند.u_c مصنوعی بزرگ می‌شود و نتیجه گمراه‌کننده است.نقشه منشأ داده بسازید و استقلال هر مؤلفه را بررسی کنید.
اختلاط U و uعدم‌قطعیت گسترش‌یافته گواهی مستقیماً مانند u وارد محاسبه می‌شود.سطح پوشش دوباره اعمال می‌شود.پیش از ورود داده، U، k و سطح پوشش را به u تبدیل و ثبت کنید.
توزیع/تقسیم‌کننده بی‌دلیلبرای همه ورودی‌ها ضریب √3 یا k=2 به‌صورت پیش‌فرض انتخاب شده است.مدل آماری با اطلاعات موجود سازگار نیست.برای هر Type B، دلیل کران و توزیع را مستند کنید.
نادیده‌گرفتن هم‌بستگیضرایب منحنی کالیبراسیون، تصحیحات مشترک یا اندازه‌گیری‌های تکراری مستقل فرض شده‌اند.ممکن است u_c کم‌ یا بیش‌برآورد شود.کوواریانس یا مدل مشترک را وارد کنید؛ استقلال را فقط با دلیل فنی بپذیرید.
LPU برای مدل غیرخطیتقریب مرتبه اول در نسبت، توان یا آستانه بدون آزمون اعمال شده است.فاصله پوشش متقارن، دم واقعی خروجی را پنهان می‌کند.MCM را برای کنترل انجام دهید و تفاوت صدک‌ها را مقایسه کنید.
Monte Carlo صوریتعداد تکرار، نسخه فایل، توزیع‌ها و آزمون پایداری ثبت نشده‌اند.خروجی قابل بازتولید یا ممیزی نیست.کنترل نسخه، گزارش تنظیمات و بررسی ثبات صدک‌های پوشش الزامی باشد.
آمیختن عدم‌قطعیت و قاعده تصمیمU گزارش شده ولی نحوه اعلام انطباق با حد مشخص نیست.یک نتیجه اندازه‌گیری به ادعای انطباقِ پرریسک تبدیل می‌شود.قاعده تصمیم، guard band و سطح ریسک را جداگانه توافق و گزارش کنید.
هشدار فنی: قاعده «عدم‌قطعیت حداکثر یک‌سوم MPE باشد» یک الزام عمومی و جایگزین مدل عدم‌قطعیت یا قاعده تصمیم نیست. تناسب تجهیز با کاربرد باید با روش، تلرانس، ریسک، قابلیت اندازه‌گیری و سیاست تصمیم بررسی شود.

نمونه پروژه‌ای بازسازی‌شده: چرا MCM برای یک روش منحنی کالیبراسیون انتخاب شد؟

این مثال، بازسازی‌شده و بی‌نام است؛ برای آموزش طراحی شده و نماینده یک مشتری مشخص نیست. هدف آن نشان‌دادن منطق انتخاب روش است، نه ارائه دستورالعمل آماده برای همه روش‌های شیمیایی.

تعریف مسئله

یک آزمایشگاه کنترل کیفیت، غلظت یک آنالیت را با استفاده از سیگنال دستگاه و منحنی کالیبراسیون تعیین می‌کند. مدل ساده‌شده به‌صورت زیر است:

C = (S − b) / m

در این مدل C غلظت، S سیگنال نمونه، b عرض از مبدأ و m شیب منحنی است. در این پرونده، عدم‌قطعیت شیب در مقایسه با سایر مؤلفه‌ها قابل‌توجه بود؛ بنابراین قرار گرفتن m در مخرج، شکل خروجی را حساس و بالقوه نامتقارن می‌کرد.

ورودیمقدار نمونه بازسازی‌شدهu استاندارد / مدل ورودیتوضیح
S؛ سیگنال نمونه0.2450.003؛ نرمالاز تکرار خوانش در شرایط کنترل‌شده.
b؛ عرض از مبدأ0.0100.002؛ نرمالاز خروجی برازش، با کنترل سازگاری داده‌ها.
m؛ شیب منحنی0.004700.00050؛ نرمال در این مثالعدم‌قطعیت شیب نسبت به مقدار آن مؤثر است و در مخرج مدل قرار دارد.
وابستگیدر مدل کامل باید بررسی شوددر مثال آموزشی برای تمرکز بر اثر غیرخطی لحاظ نشده است.در پروژه واقعی، کوواریانس b و m نباید بدون آزمون حذف شود.

مقایسه محاسبات و تصمیم فنی

خروجیLPU؛ تقریب مرتبه اولMCM؛ ۵۰۰٬۰۰۰ تکرار آموزشیبرداشت
برآورد مرکزی C50.00میانه 50.00؛ میانگین 50.58در مدل نسبت‌دار، میانگین و میانه لزوماً برابر نیستند.
عدم‌قطعیت/پوشش ۹۵٪u_c ≈ 5.37 و U(k=2) ≈ 10.75صدک 2.5٪ تا 97.5٪: 41.27 تا 63.23بازه Monte Carlo نامتقارن است و دم بالایی بلندتر دیده می‌شود.
تصمیم روشبرای کنترل مستقل مفید است.برای گزارش پوشش و تصمیم نزدیک حد، انتخاب شد.علت انتخاب MCM غیرخطی‌بودن ناشی از m در مخرج و اثر آن بر دم خروجی بود.
حد و تصمیم: اگر حد بالای مشخصات مشتری 60 باشد، صرفِ مقدار مرکزی 50 برای صدور «منطبق» کافی نیست؛ باید قاعده تصمیم، ریسک و guard band پیش از اظهار انطباق مشخص شود. در این مثال، صدک بالای 97.5٪ از 60 عبور می‌کند.

کنترل‌های لازم قبل از نهایی‌کردن چنین پرونده‌ای عبارت‌اند از: بازبینی اعتبار منحنی، بررسی کوواریانس b و m، آزمون با داده یا نمونه مرجع، کنترل تعداد تکرار تا پایداری صدک‌ها و ثبت نسخه ابزار محاسباتی. MCM مدل اشتباه را اصلاح نمی‌کند؛ فقط پیامد آماری همان مدل را با شفافیت بیشتری نشان می‌دهد.

کاربرگ کوتاه انتخاب و بازبینی روش

این کاربرگ را می‌توان پیش از تأیید بودجه عدم‌قطعیت یا بازنگری روش تکمیل کرد. اگر در هر یک از سه ردیف اول پاسخ «بله» است، بررسی MCM باید در پرونده تصمیم ثبت شود.

پرسش کنترلپاسخ/شاهدی که باید ثبت شود
آیا مدل در بازه عدم‌قطعیت غیرخطی است یا ورودی در مخرج/توان/آستانه قرار دارد؟بله/خیر + نمودار حساسیت یا محاسبه کنترل.
آیا ورودی یا خروجی توزیع کراندار، چوله یا نامتقارن دارد؟نوع توزیع، حدود و منبع توجیه.
آیا نتیجه نزدیک حد پذیرش است یا اظهار انطباق صادر می‌شود؟حد، قاعده تصمیم، مسئول تأیید و ریسک مورد توافق.
آیا هم‌بستگی میان ورودی‌ها بررسی شده است؟منشأ وابستگی و نحوه اعمال/حذف با توجیه.
آیا داده‌های Type A و Type B از نظر هم‌پوشانی و دوبارشماری کنترل شده‌اند؟نقشه منشأ داده و تأیید بازبین.
اگر MCM اجرا شده، آیا تکرار، نسخه و پایداری صدک‌ها ثبت شده‌اند؟شناسه فایل/کد، تعداد تکرار، بذر یا تنظیمات و نتیجه آزمون عددی.
آیا گزارش عدم‌قطعیت از قاعده تصمیم جدا شده است؟عبارت گزارش، سطح پوشش و متن اظهار انطباق.

گزارش نتیجه: عدم‌قطعیت، فاصله پوشش و قاعده تصمیم را مخلوط نکنید

گزارش فنی باید روشن کند چه کمیتی اندازه‌گیری شده، نتیجه مرکزی چیست، عدم‌قطعیت یا فاصله پوشش با چه سطحی بیان شده و کدام مدل/روش به کار رفته است. اگر نتیجه انطباق با یک مشخصه اعلام می‌شود، قاعده تصمیم باید جداگانه قابل تشخیص باشد؛ یعنی مخاطب بداند عدم‌قطعیت چگونه در تصمیم لحاظ شده و ریسک چه کسی پذیرفته شده است.

راهنمای ILAC G8 برای استفاده از قواعد تصمیم و اظهار انطباق در چارچوب ISO/IEC 17025 تهیه شده است. JCGM 106 نیز نقش عدم‌قطعیت را در ارزیابی انطباق توضیح می‌دهد. هیچ‌کدام مجوز نمی‌دهند که آزمایشگاه، پس از دیدن نتیجه، قاعده تصمیم دلخواه انتخاب کند.

فرمول گزارش پیشنهادی: «نتیجه … است. فاصله پوشش ۹۵٪ بر پایه [LPU/MCM] و مدل مستند روش محاسبه شده است. هر اظهار انطباق، مطابق قاعده تصمیم [شناسه/نام] صادر شده است.»

پرسش‌های پرتکرار

آیا برای همه روش‌ها باید Monte Carlo اجرا شود؟

خیر. در مدل‌های خطی یا نزدیک به خطی با مفروضات روشن و خروجی متقارن، LPU می‌تواند روش کارآمد و کاملاً قابل دفاع باشد. تصمیم باید مستند و متناسب با ریسک باشد.

اگر LPU و MCM پاسخ متفاوت دادند، کدام درست است؟

ابتدا مدل، توزیع ورودی‌ها، هم‌بستگی‌ها و پیاده‌سازی هر دو روش را بازبینی کنید. اختلاف ممکن است نشانه غیرخطی‌بودن، خروجی نامتقارن یا یک خطای داده/مدل باشد. «بزرگ‌تر بودن» به‌تنهایی معیار انتخاب نیست.

آیا تعداد ثابت تکرار برای MCM وجود دارد؟

خیر. تعداد تکرار باید با دقت عددی موردنیاز و پایداری کمیت‌های گزارش‌شده، به‌ویژه صدک‌های پوشش، توجیه شود. ثبت تعداد تکرار بدون نشان‌دادن کنترل پایداری، کفایت فنی را ثابت نمی‌کند.

آیا عدم‌قطعیت گواهی کالیبراسیون همان عدم‌قطعیت روش است؟

خیر. عدم‌قطعیت گواهی یکی از ورودی‌های بالقوه مدل است. عدم‌قطعیت روش ممکن است شامل تکرارپذیری، بازیابی، نمونه‌برداری، اثر ماتریس، تصحیحات و عوامل دیگری باشد که با هدف اندازه‌گیری مرتبط‌اند.

آیا می‌توان برای هر تجهیز از یک نسبت ثابت به MPE استفاده کرد؟

خیر. MPE، کاربرد تجهیز، روش، تلرانس محصول، ریسک تصمیم و قابلیت اندازه‌گیری باید با هم بررسی شوند. یک نسبت ثابت، جایگزین تحلیل فنی و قاعده تصمیم نمی‌شود.

جمع‌بندی مدیریتی و دعوت به اقدام

انتخاب میان LPU و Monte Carlo یک انتخاب نرم‌افزاری نیست؛ انتخاب یک مسیر دفاع‌پذیر برای تبدیل داده به نتیجه است. LPU را وقتی به‌کار بگیرید که مدل و مفروضاتش واقعاً اجازه می‌دهند؛ Monte Carlo را وقتی وارد کنید که غیرخطی‌بودن، توزیع‌ها یا ریسک تصمیم، تقریب ساده را ناکافی می‌کند.

پیشنهاد هیوا: برای بازبینی بودجه عدم‌قطعیت، انتخاب روش مناسب، تعریف معیار پذیرش و تدوین قاعده تصمیم، ارزیابی فنی روش خود را درخواست کنید. تجربه هیوا در بیش از ۱۵۰ پروژه استقرار ISO/IEC 17025 منجر به اخذ تأیید صلاحیت NACI، مبنای تمرکز ما بر نتیجه اجرایی و قابل دفاع است.