Độ Không Đảm Bảo Đo (Uncertainty of Measurement): Tính Toán & Báo Cáo Theo GUM/ISO 17025
Tìm hiểu toàn diện về Measurement Uncertainty (MU) – khái niệm bắt buộc theo ISO 17025:2017. Bài viết hướng dẫn từ nguyên lý GUM, phân loại đánh giá Type A & Type B, xây dựng Uncertainty Budget, đến cách trình bày kết quả trong báo cáo thử nghiệm/hiệu chuẩn.
Độ không đảm bảo đo (Measurement Uncertainty – MU) là thông số định lượng đặc trưng cho sự phân tán của các giá trị có thể gán cho đại lượng đo. Theo ISO 17025:2017 (Điều 7.6), mọi phòng thí nghiệm được công nhận đều phải ước lượng và báo cáo MU khi kết quả ảnh hưởng đến việc đánh giá sự phù hợp.
Bài viết này cung cấp hướng dẫn thực hành từ nền tảng lý thuyết GUM đến xây dựng Uncertainty Budget hoàn chỉnh, giúp phòng thí nghiệm đáp ứng yêu cầu kiểm định và công nhận.
1. Khái Niệm Độ Không Đảm Bảo Đo & Tầm Quan Trọng
Định nghĩa (VIM 2.26): "Tham số không âm đặc trưng cho sự phân tán của các giá trị được gán một cách hợp lý cho đại lượng đo, dựa trên thông tin sử dụng."
| Thuật ngữ | Ký hiệu | Ý nghĩa |
|---|---|---|
| Độ không đảm bảo chuẩn | u(x) | Biểu diễn ở mức 1 sigma (≈68% confidence) |
| Độ không đảm bảo chuẩn tổng hợp | uc(y) | Tổng hợp tất cả nguồn u(xi) |
| Độ không đảm bảo mở rộng | U | U = k × uc, thường k=2 (≈95%) |
| Hệ số phủ | k | Thường k=2; k=3 cho 99.7% |
Tại sao quan trọng?
- ISO 17025:2017 (7.6): Bắt buộc ước lượng MU cho kết quả định lượng
- Quyết định phù hợp (Conformity Assessment): So sánh kết quả ± U với giới hạn quy chuẩn (ILAC-G8:09/2019)
- Độ tin cậy kết quả: Kết quả không có MU là kết quả không hoàn chỉnh
2. Tài Liệu Tham Chiếu & Tiêu Chuẩn Áp Dụng
| Tài liệu | Phạm vi | Ghi chú |
|---|---|---|
| GUM (JCGM 100:2008) | Hướng dẫn chung ước lượng MU | Tài liệu nền tảng, áp dụng mọi lĩnh vực |
| ISO 17025:2017 | Yêu cầu năng lực PTN | Điều 7.6 – Ước lượng MU |
| EURACHEM/CITAC Guide | MU trong phân tích hóa học | Ví dụ thực tế chi tiết (3rd ed. 2012) |
| EA-4/02 M:2022 | Biểu diễn MU trong hiệu chuẩn | Yêu cầu của tổ chức công nhận châu Âu |
| ILAC-G8:09/2019 | Quy tắc quyết định phù hợp | Guard band, shared risk |
| ĐLVN 263:2020 | Hướng dẫn tính MU tại Việt Nam | Dựa trên GUM, áp dụng cho cơ sở hiệu chuẩn |
3. Đánh Giá Type A – Phương Pháp Thống Kê
Type A sử dụng phân tích thống kê từ dữ liệu đo lặp để ước lượng độ không đảm bảo.
3.1. Công thức cơ bản
Từ n phép đo lặp, tính:
- Giá trị trung bình: x̄ = (1/n) × Σxi
- Độ lệch chuẩn mẫu: s = √[Σ(xi - x̄)² / (n-1)]
- Độ không đảm bảo chuẩn Type A: uA = s / √n
3.2. Ví dụ thực tế – Cân phân tích
| Lần đo | Khối lượng (g) |
|---|---|
| 1 | 10.0023 |
| 2 | 10.0025 |
| 3 | 10.0021 |
| 4 | 10.0024 |
| 5 | 10.0022 |
| 6 | 10.0026 |
| 7 | 10.0023 |
| 8 | 10.0024 |
| 9 | 10.0025 |
| 10 | 10.0022 |
x̄ = 10.00235 g; s = 0.000158 g; uA = 0.000158/√10 = 0.000050 g
💡 Lưu ý: Cần tối thiểu n ≥ 6 phép đo (khuyến nghị n ≥ 10) để Type A có ý nghĩa thống kê. Với n nhỏ, phải dùng hệ số Student t thay cho k=2.
4. Đánh Giá Type B – Phương Pháp Phi Thống Kê
Type B ước lượng MU từ các nguồn thông tin khác ngoài dữ liệu đo lặp: chứng chỉ hiệu chuẩn, thông số kỹ thuật thiết bị, kinh nghiệm, dữ liệu tham khảo.
4.1. Các phân bố phổ biến
| Phân bố | Khi nào dùng | Công thức uB | Ví dụ |
|---|---|---|---|
| Chuẩn (Normal) | Chứng chỉ hiệu chuẩn cho U và k | uB = U/k | Chứng chỉ ghi U = 0.02 mg, k=2 → u = 0.01 mg |
| Chữ nhật (Rectangular) | Chỉ biết giới hạn ±a, không có thêm thông tin | uB = a/√3 | Độ phân giải ±0.1 mg → u = 0.058 mg |
| Tam giác (Triangular) | Biết giới hạn ±a, giá trị trung tâm có xác suất cao hơn | uB = a/√6 | Nhiệt độ phòng 23 ± 2°C, thường ở giữa |
4.2. Nguồn Type B thường gặp
- Chứng chỉ hiệu chuẩn: MU của chuẩn tham chiếu (CRM, quả cân, nhiệt kế chuẩn)
- Độ phân giải thiết bị: ±0.5 × digit nhỏ nhất → phân bố chữ nhật
- Độ trôi (Drift): Sai lệch giữa 2 kỳ hiệu chuẩn → phân bố chữ nhật
- Ảnh hưởng nhiệt độ: Hệ số giãn nở nhiệt × biên độ nhiệt
- Độ tinh khiết chất chuẩn: Từ CoA (Certificate of Analysis)
- Hiệu ứng matrix: Recovery từ spike experiment
5. Mô Hình Đo & Tổng Hợp Độ Không Đảm Bảo (Law of Propagation)
5.1. Thiết lập mô hình đo
Bước quan trọng nhất: viết phương trình toán học liên hệ đại lượng đo Y với các đại lượng đầu vào Xi:
Y = f(X1, X2, ..., XN)
Ví dụ – Nồng độ dung dịch pha loãng:
Cs = C0 × V0 / Vs
Trong đó: C0 = nồng độ gốc, V0 = thể tích lấy, Vs = thể tích định mức.
5.2. Luật truyền độ không đảm bảo
Với các đại lượng không tương quan:
uc²(y) = Σ [ci × u(xi)]²
Trong đó ci = ∂f/∂xi là hệ số nhạy (sensitivity coefficient).
5.3. Trường hợp đặc biệt – Công thức tương đối
Khi mô hình chỉ gồm phép nhân/chia:
[uc(y)/y]² = Σ [u(xi)/xi]²
Đây là cách nhanh nhất để tổng hợp khi biết relative uncertainty (%) của từng thành phần.
5.4. Xử lý tương quan (Correlation)
Nếu hai đại lượng có tương quan (ví dụ: cùng hiệu chuẩn bằng một chuẩn), phải thêm hạng tử hiệp phương sai: 2 × ci × cj × u(xi) × u(xj) × r(xi, xj)
6. Xây Dựng Uncertainty Budget – Bảng Tổng Hợp
Uncertainty Budget là bảng tổng hợp tất cả các nguồn đóng góp, là yêu cầu cốt lõi khi trình bày MU.
6.1. Ví dụ: Uncertainty Budget cho phép cân mẫu trên cân phân tích
| Nguồn | Giá trị | Phân bố | Ước số | u(xi) (mg) | ci | ui(y) (mg) | % đóng góp |
|---|---|---|---|---|---|---|---|
| Lặp lại (Type A, n=10) | s = 0.05 mg | Normal | √10 | 0.0158 | 1 | 0.0158 | 3.6% |
| Chứng chỉ hiệu chuẩn | U = 0.10 mg, k=2 | Normal | 2 | 0.0500 | 1 | 0.0500 | 36.1% |
| Độ phân giải (2 lần đọc) | d = 0.1 mg | Rectangular | √3 | 0.0408 | √2 | 0.0577 | 48.1% |
| Độ lệch tâm (Eccentricity) | ±0.05 mg | Rectangular | √3 | 0.0289 | 1 | 0.0289 | 12.1% |
uc = √(0.0158² + 0.0500² + 0.0577² + 0.0289²) = 0.0832 mg
U = 2 × 0.0832 = 0.17 mg (k=2, p ≈ 95%)
📝 Nhận xét: Độ phân giải (48.1%) và chứng chỉ hiệu chuẩn (36.1%) là hai nguồn đóng góp lớn nhất → ưu tiên cải thiện bằng cách dùng cân có độ phân giải cao hơn hoặc chuẩn có MU nhỏ hơn.
7. Báo Cáo Kết Quả & Quy Tắc Làm Tròn
7.1. Cách trình bày theo GUM & ISO 17025
Kết quả đo phải bao gồm:
- Giá trị đo y
- Độ không đảm bảo mở rộng U
- Hệ số phủ k và mức tin cậy p
Ví dụ chuẩn:
"Khối lượng mẫu: m = 10.0024 ± 0.00017 g
Độ không đảm bảo mở rộng U được tính với hệ số phủ k = 2, tương ứng mức tin cậy xấp xỉ 95%."
7.2. Quy tắc làm tròn (EA-4/02)
- U làm tròn lên tối đa 2 chữ số có nghĩa
- Giá trị đo làm tròn phù hợp với vị trí thập phân cuối cùng của U
- Không bao giờ làm tròn U xuống (để không đánh giá thấp MU)
7.3. Quyết định phù hợp (ILAC-G8)
| Phương pháp | Mô tả | Rủi ro |
|---|---|---|
| Simple Acceptance | Kết quả nằm trong giới hạn → Đạt (không xét MU) | Rủi ro chấp nhận sai tới 50% tại biên |
| Guard Band (W = U) | Thu hẹp giới hạn chấp nhận: LSL+U đến USL−U | Rủi ro chấp nhận sai < 2.5% |
| Shared Risk | Hai bên thỏa thuận mức rủi ro chấp nhận | Theo thỏa thuận |
8. Ứng Dụng Tính MU Theo Lĩnh Vực
8.1. Phân tích hóa học (EURACHEM/CITAC)
Các nguồn MU đặc thù:
- Chuẩn bị mẫu: cân, pha loãng, chiết xuất (recovery)
- Thiết bị phân tích: đường chuẩn (calibration curve), LOD/LOQ
- Matrix effect: spike recovery ± %RSD
- Phương pháp top-down: sử dụng dữ liệu QC, proficiency testing
8.2. Hiệu chuẩn thiết bị (EA-4/02)
- MU của chuẩn tham chiếu (traceable to SI)
- Độ phân giải thiết bị cần hiệu chuẩn
- Độ trôi giữa các kỳ hiệu chuẩn
- Ảnh hưởng môi trường (nhiệt độ, độ ẩm, áp suất)
8.3. Vi sinh (ISO 19036)
MU trong vi sinh chủ yếu dựa trên Type A (độ tái lặp) do khó kiểm soát nguồn Type B. Thường biểu diễn bằng Standard Deviation of Reproducibility (sR) tính trên log10 CFU.
8.4. Quy tắc tỷ lệ TUR (Test Uncertainty Ratio)
TUR = Tolerance / (2 × Ucal)
Yêu cầu: TUR ≥ 4:1 (lý tưởng). Khi TUR < 4:1, phải áp dụng guard band hoặc nêu rõ trong báo cáo.
9. Câu Hỏi Thường Gặp & Tài Nguyên
FAQ
- MU có bắt buộc ghi trên báo cáo thử nghiệm không?
Theo ISO 17025:2017 (7.8.3): Phải báo cáo MU khi: (1) ảnh hưởng đến đánh giá phù hợp, (2) khách hàng yêu cầu, (3) MU ảnh hưởng đến giới hạn quy chuẩn. - Phân biệt Error và Uncertainty?
Error = giá trị đo − giá trị thật (có thể hiệu chính). Uncertainty = vùng phân tán quanh kết quả (không thể loại bỏ, chỉ ước lượng). - Khi nào dùng Type A, khi nào Type B?
Type A khi có dữ liệu đo lặp. Type B khi thông tin đến từ chứng chỉ, catalog, kinh nghiệm. Thực tế luôn dùng cả hai. - Top-down vs Bottom-up approach?
Bottom-up (GUM): liệt kê từng nguồn → phù hợp hiệu chuẩn. Top-down: dùng dữ liệu validation, QC, PT → phù hợp phân tích hóa/vi sinh. - HTVLAB hỗ trợ gì?
HTVLAB cung cấp thiết bị đo lường chính xác cao kèm chứng chỉ hiệu chuẩn có truy xuất nguồn gốc, hỗ trợ tư vấn IQ/OQ/PQ và tính toán MU. Liên hệ 0909.860.489 hoặc Sales@htvsci.com.
Tài nguyên tham khảo
- JCGM 100:2008 (GUM)
- EURACHEM/CITAC Guide QUAM
- EA-4/02 M:2022 – Evaluation of the Uncertainty of Measurement in Calibration
- ILAC-G8:09/2019 – Guidelines on Decision Rules