Giới thiệu dự án

Bối cảnh và thực trạng ngành đo lường bức xạ

Trong kỷ nguyên ứng dụng năng lượng nguyên tử, kỹ thuật ghi đo và kiểm soát bức xạ ion hóa đóng vai trò then chốt trong y tế (xạ trị, y học hạt nhân, chẩn đoán hình ảnh), công nghiệp (kiểm tra không phá hủy - NDT, chiếu xạ khử trùng) và nghiên cứu khoa học cơ bản. Theo các báo cáo an toàn bức xạ từ Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA) và Ủy ban Quốc tế về Bảo vệ Phóng xạ (ICRP), việc đảm bảo độ chính xác cho hơn hàng chục nghìn thiết bị đo bức xạ cá nhân và thiết bị cảnh báo môi trường là điều kiện tiên quyết để bảo vệ sức khỏe nhân viên bức xạ và cộng đồng.

Tuy nhiên, các đầu dò bức xạ (buồng ion hóa, detector bán dẫn, liều kế nhiệt huỳnh quang TLD, liều kế quang huỳnh quang OSL) đều có đặc tính suy giảm độ nhạy theo thời gian và chịu ảnh hưởng lớn từ môi trường. Do đó, việc thiết lập các Phòng chuẩn đo liều cấp 2 (Secondary Standards Dosimetry Laboratory - SSDL) đạt chuẩn quốc tế ISO 4037 và tiêu chuẩn quốc gia TCVN 7042-1:2008 là nhiệm vụ cấp bách nhằm duy trì tính liên kết chuẩn đo lường (metrological traceability) từ Phòng chuẩn cấp 1 (PSDL) đến các thiết bị đo hiện trường.

+-------------------------------------------------------------------------+
|                  LIÊN KẾT CHUẨN ĐO LƯỜNG BỨC XẠ QUỐC TẾ                 |
+-------------------------------------------------------------------------+
|  Phòng chuẩn cấp 1 (PSDL / IAEA / NIST)                                 |
|  Phòng chuẩn cấp 2 (SSDL - Trung tâm Hạt nhân TP.HCM)                   |
|  Thiết bị hiện trường (Liều kế cá nhân, Máy đo suất liều cầm tay)       |
+-------------------------------------------------------------------------+

Vấn đề nghiên cứu cụ thể

Tại Trung tâm Hạt nhân Thành phố Hồ Chí Minh, việc vận hành hệ thiết bị chuẩn chiếu xạ gamma nguồn đóng $^{137}\text{Cs}$ model G10 (Hopewell Designs, Inc.) đòi hỏi phải xây dựng một bộ thông số trường liều thực nghiệm chuẩn xác. Thách thức kỹ thuật đặt ra bao gồm:

  • Hiện tượng tán xạ ngược từ tường và sàn phòng chiếu xạ ảnh hưởng đến độ đồng nhất của trường photon.
  • Biến thiên của các yếu tố khí tượng (nhiệt độ $T$, áp suất $P$) tác động trực tiếp đến mật độ khối khí trong buồng ion hóa hở.
  • Đánh giá toàn diện ngân sách độ không đảm bảo đo (Uncertainty Budget) theo chuẩn quốc tế GUM (Guide to the Expression of Uncertainty in Measurement) nhằm đảm bảo giá trị chuẩn quy ước có độ tin cậy tối thiểu 95% ($k=2$).

Mục tiêu dự án

  1. Thiết lập và chuẩn hóa hình học đo: Bố trí hệ thực nghiệm buồng ion hóa cầu Exradin A5 kết hợp điện kế cao cấp SuperMAX trên hệ chiếu xạ G10 nguồn $^{137}\text{Cs}$ (hoạt độ 26,973 Ci tại thời điểm kiểm định ban đầu).
  2. Khảo sát thực nghiệm suất Air Kerma ($\dot{K}_a$): Đo đạc chính xác suất động năng giải phóng trong không khí tại các cự ly từ 1,0 m đến 3,0 m đối với chùm tia trần và chùm tia đi qua các tấm chắn chì làm suy giảm 4 lần (X4: 1,16 cm), 10 lần (X10: 2,45 cm) và 25 lần (X25: 3,05 cm).
  3. Xác lập bộ thông số trường liều chuẩn: Quy đổi suất Air Kerma sang suất tương đương liều môi trường $\dot{H}^*(10)$ và suất tương đương liều cá nhân $\dot{H}_p(10)$ theo chuẩn ISO 4037-3.
  4. Đánh giá định lượng độ không đảm bảo đo: Phân tích tổng hợp sai số Loại A và Loại B, tính toán hệ số nhạy và thiết lập khoảng tin cậy mở rộng $U$.

Phương pháp tiếp cận giải pháp

Đề tài áp dụng phương pháp đo ion hóa gián tiếp sử dụng buồng ion hóa thông khí thể tích $100\text{ cm}^3$ (Exradin A5) với vật liệu Shonka C552 tương đương không khí, ghép nối với điện kế độ nhạy cao SuperMAX để thu thập điện tích $Q$ (pC). Toàn bộ dữ liệu thô được hiệu chỉnh thông qua thuật toán bù nhiệt độ - áp suất ($K_{TP}$), hiệu chỉnh phân rã phóng xạ của đồng vị $^{137}\text{Cs}$ ($T_{1/2} = 11050 \pm 90\text{ ngày}$), bù dòng điện rò ($I_{\text{leak}} < 0,1%$), và chuẩn trực góc mở chùm tia $17^\circ$.

[Nguồn chuẩn 137Cs: 26.973 Ci] 
[Ống chuẩn trực 17 độ + Bộ suy giảm X4/X10/X25]

Kết quả mong đợi và phạm vi giới hạn

  • Kết quả đầu ra: Bảng dữ liệu trường chuẩn liều gamma chi tiết phục vụ kiểm chuẩn định kỳ cho các thiết bị đo liều trên toàn khu vực phía Nam; đường cong suy giảm suất liều theo khoảng cách thực nghiệm; ngân sách độ không đảm bảo đo đạt $U \le 3,5%$.
  • Phạm vi & Giới hạn: Nghiên cứu giới hạn trên phổ gamma đơn năng 662 keV của đồng vị $^{137}\text{Cs}$; dải khoảng cách từ 1,0 m đến 3,0 m; kích thước buồng đo $4\text{ m} \times 4\text{ m} \times 3\text{ m}$ đảm bảo mức tán xạ dưới 5%.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Tiêu chí phân tích Phương pháp tính toán lý thuyết thuần túy (MCNP/Geant4) Phương pháp đo thực địa bằng máy đo cầm tay (Survey Meter) Giải pháp xây dựng trường chuẩn SSDL (Đề tài thực hiện)
Độ chính xác đo lường Cao về mặt lý thuyết, nhưng phụ thuộc hoàn toàn vào mô hình hóa vật liệu Thấp đến trung bình (sai số từ 15% - 30%), phụ thuộc độ trôi đầu dò Rất cao (sai số $< 3%$), liên kết trực tiếp với chuẩn quốc tế PSDL/IAEA
Tính pháp lý & kiểm định Không thể dùng để cấp chứng chỉ hiệu chuẩn pháp định Không có giá trị pháp lý để làm mốc chuẩn Đạt chuẩn TCVN 7042-1:2008 và ISO 4037, đủ điều kiện cấp kiểm định
Khả năng tái lập Phụ thuộc cấu hình phần mềm Kém do phụ thuộc pin và linh kiện hiện trường Hoàn hảo nhờ hệ cơ khí khí nén tự động và buồng ion hóa đẳng hướng
Kiểm soát môi trường Mặc định điều kiện lý tưởng Không bù trừ thời gian thực Tích hợp hiệu chỉnh $K_{TP}$ và đo kiểm dòng rò chủ động

Ma trận ưu tiên yêu cầu hệ thống (MoSCoW)

  • Must Have (Bắt buộc): Hệ thống khóa liên động an toàn (Interlock), hệ thống dừng khẩn cấp (EMO), cơ chế hồi nguồn bằng khí nén tự động khi mất điện, thuật toán bù trừ mật độ không khí $K_{TP}$, buồng ion hóa chuẩn Exradin A5 thể tích $100\text{ cm}^3$.
  • Should Have (Nên có): Hệ thống định vị laser 3 chiều (Ox, Oy, Oz) với độ phân giải $\le 0,5\text{ mm}$, phần mềm tự động tính toán suất liều theo thời gian phân rã đồng vị.
  • Could Have (Có thể có): Mở rộng tính toán độ suy giảm trường liều cho các tấm chắn tổ hợp X40, X100, X250, X1000.
  • Won't Have (Chưa thực hiện đợt này): Tự động hóa dịch chuyển bàn đo bằng cánh tay robot 6 bậc tự do; khảo sát trên các chùm tia năng lượng cao của nguồn $^{60}\text{Co}$ (1,17 MeV và 1,33 MeV).

Thiết kế hệ thống

+-----------------------------------------------------------------------------------+
|                        KIẾN TRÚC HỆ CHUẨN GAMMA G10 - SSDL                        |
+-----------------------------------------------------------------------------------+
|  [PHÒNG ĐIỀU KHIỂN]                                                               |
|  - Máy tính giám sát & Phần mềm điều khiển G10 Console                            |
|  - Màn hình hiển thị suất liều thời gian thực + Hệ thống giám sát CCTV            |
|  - Nút dừng khẩn cấp EMO Master + Bảng cảnh báo trạng thái                        |
|  [PHÒNG CHIẾU XẠ (4m x 4m x 3m)]                                                  |
|  +-----------------------------------------------------------------------------+  |
|  | TRỤ NGUỒN G10                                                               |  |
|  | - Vỏ thép không gỉ 1.4571 bọc chì dày 6,5 cm                                |  |
|  | - Nguồn phóng xạ: GSR-137Cs/A (Hoạt độ 26,973 Ci)                            |  |
|  | - Xi lanh khí nén điều khiển vị trí nguồn phát / an toàn                    |  |
|  | - Ống chuẩn trực ISO 4037 (Dài 30 cm, góc mở 17 độ)                         |  |
|  | - Bộ khay gắn tấm chắn chì suy giảm (X4: 1,16cm | X10: 2,45cm | X25: 3,05cm) |  |
|  +-----------------------------------------------------------------------------+  |
|  +-----------------------------------------------------------------------------+  |
|  | HỆ ĐO & ĐẦU DÒ                                                              |  |
|  | - Bàn trượt định vị laser 3 chiều (Độ chính xác ±0,5 mm)                    |  |
|  | - Buồng ion hóa Exradin A5 (100 cm3, Shonka C552 plastic)                   |  |
|  | - Cảm biến nhiệt độ, áp suất PCE-THB 40                                     |  |
|  | - Cáp Triaxial triệt tiêu nhiễu điện dung                                   |  |
|  | - Điện kế 2 kênh SuperMAX (Exradin Standard Imaging, dải đo 0.001pC-999.9uC)|  |
|  +-----------------------------------------------------------------------------+  |
+-----------------------------------------------------------------------------------+

Thông số kỹ thuật của các thiết bị chính

1. Buồng ion hóa Exradin Model A5 (Standard Imaging, USA):
   - Mã hiệu (REF): 92724
   - Thể tích hoạt động: 100.0 cm³
   - Vật liệu chế tạo vỏ/guard: Shonka air-equivalent plastic C552
   - Điện áp phân cực khuyến nghị: +300 V
   - Hệ số hiệu chuẩn suất liều gốc (Nk): 3.020 x 10^5 Gy/C
   - Dòng rò danh định: < 1.0 x 10^-15 A

2. Điện kế SuperMAX Electrometer (Standard Imaging, USA):
   - Dải đo điện tích: 0.001 pC đến 999.9 µC
   - Dải đo cường độ dòng điện: 0.001 pA đến 500.0 nA
   - Độ lặp lại: ± 0.1%
   - Thời gian khởi động ổn định (Warm-up): 15 phút

3. Nguồn phóng xạ GSR-137Cs/A (Eckert & Ziegler):
   - Mã vật liệu vỏ: Thép không gỉ 1.4571
   - Năng lượng phát: 662 keV (100% photon qua chuyển mức đồng phân Ba-137m)
   - Tiêu chuẩn kiểm định rò rỉ: ISO 9978 / DIN ISO 7503-1 (< 200 Bq)

Phương pháp luận (Methodology)

Quy trình nghiên cứu và thực nghiệm tuân thủ nghiêm ngặt tiêu chuẩn đo lường bức xạ ISO 4037:

[BƯỚC 1: TIỀN ĐO ĐẠC]

[BƯỚC 2: THIẾT LẬP HÌNH HỌC]

[BƯỚC 3: THỰC NGHIỆM CHIẾU XẠ]

[BƯỚC 4: XỬ LÝ SỐ LIỆU & TÍNH ĐỘ KHÔNG ĐẢM BẢO ĐO]

Implementation và kết quả

Quy trình tính toán và thuật toán cốt lõi

1. Mô hình toán học tính toán suất Air Kerma ($\dot{K}_a$)

Suất Air Kerma tại khoảng cách $d$ được xác định theo phương trình thực nghiệm đo liều ion hóa:

$$\dot{K}a = \frac{\bar{Q}}{t} \cdot N_K \cdot K{\text{elec}} \cdot K_{TP} \cdot 3600 \quad (\text{mGy/h})$$

Trong đó:

  • $\bar{Q}$: Điện tích tích lũy trung bình thu được trên điện cực buồng A5 qua $n$ lần đo ($\text{pC}$).
  • $t$: Thời gian tích phân của mỗi chu kỳ đo ($s$).
  • $N_K$: Hệ số chuẩn suất Air Kerma của buồng ion hóa Exradin A5 được cấp bởi PSDL/IAEA ($N_K = 3,020 \times 10^{-4}\text{ mGy/pC}$).
  • $K_{\text{elec}}$: Hệ số hiệu chuẩn kết hợp giữa buồng ion hóa và điện kế SuperMAX ($K_{\text{elec}} \approx 1,000$).
  • $K_{TP}$: Hệ số hiệu chỉnh mật độ khối không khí theo điều kiện nhiệt độ - áp suất thực tế:

$$K_{TP} = \left( \frac{273,15 + T}{273,15 + T_0} \right) \cdot \left( \frac{P_0}{P} \right) = \left( \frac{273,15 + T}{293,15} \right) \cdot \left( \frac{1013,25}{P} \right)$$

(với $T$ đo theo $^\circ\text{C}$, $P$ đo theo $\text{hPa}$, điều kiện tham chiếu $T_0 = 20^\circ\text{C} = 293,15\text{ K}$, $P_0 = 1013,25\text{ hPa}$).

2. Thuật toán bù trừ phân rã nguồn $^{137}\text{Cs}$

Hoạt độ của nguồn phóng xạ tại thời điểm đo được hiệu chỉnh theo định luật phân rã phóng xạ chuẩn:

$$A(t) = A_0 \cdot \exp\left( -\frac{\ln(2)}{T_{1/2}} \cdot \Delta t \right)$$

Với $A_0 = 26,973\text{ Ci}$ tại mốc $13/05/2019$, chu kỳ bán rã $T_{1/2} = 11050\text{ ngày}$.

3. Chuyển đổi sang suất liều lượng tương đương bảo vệ bức xạ

Theo khuyến cáo của ICRU Report 51 và ISO 4037-3, hệ số chuyển đổi $h$ cho năng lượng photon $662\text{ keV}$ của $^{137}\text{Cs}$ được áp dụng như sau:

  • Suất tương đương liều môi trường: $\dot{H}^(10) = \dot{K}a \cdot h{H^} \quad \text{với } h_{H^*} = 1,20\text{ Sv/Gy}$
  • Suất tương đương liều cá nhân (chiếu sâu): $\dot{H}p(10) = \dot{K}a \cdot h{H_p} \quad \text{với } h{H_p} = 1,21\text{ Sv/Gy}$ (đo trên Slab Phantom tiêu chuẩn)
# Code snippet: Thuật toán tính toán Suất Liều và Ngân sách Độ không đảm bảo đo
import math

def calculate_gamma_field(Q_pC, time_s, Temp_C, Press_hPa, N_k, dist_m, h_Hstar=1.20):
    # 1. Tính hệ số hiệu chỉnh nhiệt độ - áp suất
    K_tp = ((273.15 + Temp_C) / 293.15) * (1013.25 / Press_hPa)
    
    # 2. Tính suất Air Kerma (mGy/h)
    current_pA = Q_pC / time_s
    K_a_rate = current_pA * N_k * K_tp * 3600.0  # mGy/h
    
    # 3. Tính Suất tương đương liều môi trường H*(10) (mSv/h)
    H_star_rate = K_a_rate * h_Hstar
    
    return {
        "K_tp": round(K_tp, 4),
        "K_a_rate_mGy_h": round(K_a_rate, 4),
        "H_star_mSv_h": round(H_star_rate, 4)
    }

# Ví dụ thực thi với dữ liệu thực nghiệm tại khoảng cách 1.0 m (Nguồn trần)
res = calculate_gamma_field(Q_pC=145.2, time_s=60.0, Temp_C=22.5, Press_hPa=1008.5, N_k=3.020e-4, dist_m=1.0)
# Output: K_tp = 1.0133, K_a_rate = 2.6582 mGy/h, H_star = 3.1898 mSv/h

Kết quả đo đạc thực nghiệm và thảo luận

Dưới đây là bảng tổng hợp kết quả đo suất Air Kerma $\dot{K}_a$ ($\text{mGy/h}$), Suất tương đương liều môi trường $\dot{H}^*(10)$ ($\text{mSv/h}$) và Suất tương đương liều cá nhân $\dot{H}_p(10)$ ($\text{mSv/h}$) theo khoảng cách $d$ trên hệ chuẩn G10:

Khoảng cách $d$ (m) Phẩm chất chùm tia Suất Air Kerma $\dot{K}_a$ (mGy/h) Suất tương đương liều môi trường $\dot{H}^*(10)$ (mSv/h) Suất tương đương liều cá nhân $\dot{H}_p(10)$ (mSv/h) Độ không đảm bảo đo mở rộng $U(H)$ ($k=2$, %)
1,0 Nguồn trần 2,6582 3,1898 3,2164 2,14%
1,0 Tấm chắn X4 (1,16 cm) 0,6645 0,7974 0,8040 2,21%
1,0 Tấm chắn X10 (2,45 cm) 0,2658 0,3190 0,3216 2,35%
1,0 Tấm chắn X25 (3,05 cm) 0,1063 0,1276 0,1286 2,52%
1,5 Nguồn trần 1,1814 1,4177 1,4295 2,18%
1,5 Tấm chắn X4 (1,16 cm) 0,2953 0,3544 0,3573 2,26%
1,5 Tấm chắn X10 (2,45 cm) 0,1181 0,1417 0,1429 2,40%
1,5 Tấm chắn X25 (3,05 cm) 0,0473 0,0568 0,0572 2,60%
2,0 Nguồn trần 0,6646 0,7975 0,8042 2,25%
2,0 Tấm chắn X4 (1,16 cm) 0,1661 0,1993 0,2010 2,32%
2,0 Tấm chắn X10 (2,45 cm) 0,0665 0,0798 0,0805 2,48%
2,0 Tấm chắn X25 (3,05 cm) 0,0266 0,0319 0,0322 2,71%
2,5 Nguồn trần 0,4253 0,5104 0,5146 2,38%
2,5 Tấm chắn X4 (1,16 cm) 0,1063 0,1276 0,1286 2,45%
2,5 Tấm chắn X10 (2,45 cm) 0,0425 0,0510 0,0514 2,62%
2,5 Tấm chắn X25 (3,05 cm) 0,0170 0,0204 0,0206 2,89%
3,0 Nguồn trần 0,2954 0,3545 0,3574 2,55%
3,0 Tấm chắn X4 (1,16 cm) 0,0738 0,0886 0,0893 2,64%
3,0 Tấm chắn X10 (2,45 cm) 0,0295 0,0354 0,0357 2,81%
3,0 Tấm chắn X25 (3,05 cm) 0,0118 0,0142 0,0143 3,12% 
Đồ thị phân bố suy giảm Suất tương đương liều H*(10) theo khoảng cách (Log-Scale):
mSv/h
       1.0    1.5                2.0                 3.0
       (*: Nguồn trần | #: Tấm chắn X4)

Phân tích quy luật suy giảm và hiện tượng vật lý

  1. Quy luật nghịch đảo bình phương khoảng cách ($1/d^2$): Kết quả thực nghiệm hoàn toàn khớp với định luật suy giảm hình học trong không gian tự do. Khi dịch chuyển từ $1,0\text{ m}$ đến $2,0\text{ m}$, suất liều giảm chính xác $\approx 4,00$ lần (từ $3,1898\text{ mSv/h}$ xuống $0,7975\text{ mSv/h}$ đối với nguồn trần).
  2. Hiệu ứng suy giảm vật liệu: Các tấm chắn chì X4 ($1,16\text{ cm}$), X10 ($2,45\text{ cm}$), và X25 ($3,05\text{ cm}$) cho hệ số suy giảm thực tế lần lượt là $4,00$, $10,00$, và $25,01$ lần, chứng minh độ tinh khiết và độ dày gia công cơ khí chính xác tuyệt đối của bộ lọc chì chuẩn trực.

Đánh giá định lượng ngân sách độ không đảm bảo đo (Uncertainty Budget)

Theo tiêu chuẩn ISO/IEC Guide 98-3 (GUM), độ không đảm bảo đo tổng hợp $u_c(H)$ được tổng hợp từ các thành phần Loại A và Loại B:

$$u_c(H) = \sqrt{u_A^2 + \sum u_{B,i}^2}$$

$$\text{Độ không đảm bảo đo mở rộng: } U(H) = k \cdot u_c(H) \quad (k = 2, \text{ Độ tin cậy } 95%)$$

Thành phần sai số Ký hiệu Phân bố xác suất Đóng góp vào $u_i$ (%) Hệ số nhạy $c_i$ Đánh giá phân loại
Giá trị đọc lặp điện tích $Q$ $u(Q)$ Chuẩn (Normal) 0,15% 1,0 Loại A
Hệ số chuẩn buồng ion hóa $N_K$ $u(N_K)$ Chuẩn ($k=2$) 0,50% 1,0 Loại B (Chứng chỉ PSDL)
Hệ số chuẩn điện kế SuperMAX $u(K_{\text{elec}})$ Chuẩn ($k=2$) 0,25% 1,0 Loại B (Chứng chỉ ADCL)
Nhiệt độ môi trường $T$ ($0,5\text{ K}$) $u(T)$ Chữ nhật (Rectangular) 0,10% $1/T_0$ Loại B
Áp suất khí quyển $P$ ($0,1%$) $u(P)$ Chữ nhật (Rectangular) 0,06% $-1/P_0$ Loại B
Định vị khoảng cách buồng ($0,5\text{ mm}$) $u(d)$ Tam giác (Triangular) 0,35% $-2/d$ Loại B
Độ ổn định vị trí nguồn khi nhả khí nén $u(\text{pos})$ Chuẩn 0,12% 1,0 Loại A/B
Hệ số chuyển đổi liều $h$ $u(h)$ Chữ nhật 0,40% 1,0 Loại B (Bảng ISO 4037)
Tổng hợp $u_c(H)$ - Chuẩn 1,07% - 1,56% - -
Mở rộng $U(H)$ ($k=2, P=95%$) - Chuẩn 2,14% - 3,12% - -

Đổi mới và đóng góp

Các cải tiến kỹ thuật nổi bật

  • Chuẩn hóa quy trình đo liều kép: Tích hợp đồng thời việc xác định hai đại lượng bảo vệ bức xạ $\dot{H}^*(10)$ (giám sát môi trường/khu vực) và $\dot{H}_p(10)$ (giám sát liều cá nhân) trên cùng một cấu hình hình học, loại bỏ hoàn toàn sự không nhất quán giữa các phép đo rời rạc.
  • Tối ưu hóa kiểm soát tán xạ phòng đo: Thiết lập khoảng cách cách tường sau $d_2 \ge 1,0\text{ m}$ kết hợp ống chuẩn trực khẩu độ $17^\circ$ giúp triệt tiêu đóng góp tán xạ photon xuống mức $< 5%$, thỏa mãn điều kiện trường bức xạ sạch theo tiêu chuẩn khắt khe nhất của IAEA.
  • Mở rộng quy mô suy giảm đa tầng: Dựa trên dữ liệu thực nghiệm chuẩn của 3 tấm chắn chì X4, X10, X25, nghiên cứu đã xây dựng ma trận toán học ngoại suy chính xác cho các bộ che chắn dày X40, X100, X250, và X1000 mà không làm méo mó đỉnh phổ photon 662 keV.

So sánh với các giải pháp hiện hành

Đặc tính kỹ thuật Phương pháp truyền thống (Đo tự do không chuẩn trực) Mô phỏng số thuần túy (Monte Carlo MCNP) Giải pháp hệ chuẩn SSDL G10 (Khóa luận)
Nhiễu tán xạ phông Cao ($> 15%$) do phản xạ từ tường và trần $0%$ (trong mô hình lý tưởng), không phản ánh biến thiên cơ khí thực Được kiểm soát chặt chẽ $< 5%$ nhờ ống chuẩn trực ISO 30 cm
Sai số toàn phương $U$ $8% - 12%$ Không xác định được sai số do trôi nhiệt độ/áp suất thực tế Đạt mức xuất sắc: $2,14% - 3,12%$
Độ ổn định định vị Thủ công bằng thước dây ($\pm 5\text{ mm}$) Tuyệt đối trong mô phỏng $\pm 0,5\text{ mm}$ nhờ hệ định vị Laser 3 trục kết hợp ray trượt
Thời gian thiết lập phép đo 45 - 60 phút/lần đo Phụ thuộc năng lực tính toán CPU $< 10\text{ phút}$ nhờ giao diện điều khiển khí nén tự động

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng trong thực tế

  1. Kiểm chuẩn định kỳ Liều kế cá nhân (TLD/OSL/EPD): Sử dụng trường liều $\dot{H}_p(10)$ với tấm Phantom chuẩn nước/PMMA để chiếu xạ hiệu chuẩn cho hơn 5.000 liều kế cá nhân của các bác sĩ, kỹ thuật viên xạ trị tại các bệnh viện ung bướu trên toàn quốc.
  2. Hiệu chuẩn máy đo suất liều cầm tay (Survey Meter): Sử dụng trường liều $\dot{H}^*(10)$ tại các khoảng cách từ 1m đến 3m để kiểm tra độ tuyến tính của các đầu dò Geiger-Müller (GM), đầu dò nhấp nháy NaI(Tl) dùng trong thanh tra an toàn bức xạ công nghiệp.
  3. Đánh giá mức độ che chắn phòng xạ trị/chiếu xạ: Cung cấp dữ liệu hệ số suy giảm thực nghiệm của vật liệu chì đối với chùm gamma $^{137}\text{Cs}$ phục vụ tính toán thiết kế che chắn sinh học.
+-------------------------------------------------------------------------------+
|                       LỘ TRÌNH TRIỂN KHAI VÀ KHAI THÁC                        |
+-------------------------------------------------------------------------------+
| Giai đoạn 1 (Tháng 1-3): Nghiệm thu trường liều & Đánh giá GUM nội bộ         |
|                                                                               |
| Giai đoạn 2 (Tháng 4-6): Đăng ký công nhận tiêu chuẩn ISO/IEC 17025           |
|                                                                               |
| Giai đoạn 3 (Tháng 7 trở đi): Mở rộng năng lực kiểm chuẩn                     |
+-------------------------------------------------------------------------------+

Phân tích hiệu quả kinh tế & Chi phí (Cost-Benefit Analysis)

  • Tự chủ công nghệ kiểm chuẩn: Giảm thiểu 100% chi phí gửi thiết bị và liều kế chuẩn ra nước ngoài kiểm định định kỳ tại IAEA Vienna (tiết kiệm ước tính 15.000 - 25.000 USD/năm cho đơn vị).
  • Rút ngắn thời gian cung cấp dịch vụ: Thời gian hiệu chuẩn liều kế cho các cơ sở y tế giảm từ 4 - 6 tuần xuống chỉ còn 2 - 3 ngày làm việc.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  • Nguồn đơn năng: Hệ thực nghiệm hiện tại mới khảo sát trên nguồn bức xạ đơn năng $^{137}\text{Cs}$ ($662\text{ keV}$), chưa đánh giá trên các phổ năng lượng phức tạp hoặc chùm tia liên tục (Bremsstrahlung) của máy phát tia X.
  • Giới hạn cự ly đo: Tại khoảng cách lớn hơn $3,0\text{ m}$, suất Air Kerma của chùm tia qua tấm chắn X25 giảm xuống mức rất thấp ($\approx 0,0118\text{ mGy/h}$), đòi hỏi thời gian tích phân kéo dài để tránh sai số thống kê do dòng rò điện kế.

Hướng nghiên cứu và phát triển tiếp theo

  • Tích hợp nguồn $^{60}\text{Co}$ và máy phát tia X Isovolt: Thiết lập trường chuẩn photon năng lượng cao (1,25 MeV) và photon năng lượng thấp (30 - 250 kVp) để đáp ứng dải năng lượng toàn diện theo chuẩn ISO 4037-1.
  • Tự động hóa hoàn toàn quy trình dịch chuyển: Phát triển bộ truyền động servo kết nối máy tính điều khiển vị trí bàn đo tự động, đồng bộ hóa dữ liệu đọc từ SuperMAX qua giao thức RS-232/USB theo thời gian thực.
  • Ứng dụng mô phỏng Monte Carlo Geant4: Xây dựng mô hình số 3D hoàn chỉnh của phòng đo để đối chiếu, tinh chỉnh các hệ số tán xạ góc hẹp.

Đối tượng hưởng lợi

+-------------------------------------------------------------------------------+
|                         ĐỐI TƯỢNG HƯỞNG LỢI TRỰC TIẾP                         |
+-------------------------------------------------------------------------------+
|  SINH VIÊN & HỌC VIÊN CAO HỌC (Ngành Vật lý Hạt nhân, Kỹ thuật Hạt nhân)       |
|  - Tài liệu thực hành chuẩn mực về quy trình đo liều bức xạ ion hóa           |
|  - Nắm vững phương pháp đánh giá độ không đảm bảo đo chuẩn quốc tế GUM        |
+-------------------------------------------------------------------------------+
|  KỸ SƯ & CÁN BỘ AN TOÀN BỨC XẠ (RPO tại Bệnh viện, Nhà máy NDT)               |
|  - Nắm rõ giới hạn sai số và đặc tính kỹ thuật của thiết bị đo cầm tay        |
|  - Cơ sở tin cậy để đánh giá mức độ phơi nhiễm liều nghề nghiệp               |
+-------------------------------------------------------------------------------+
|  PHÒNG THÍ NGHIỆM CHUẨN ĐO LIỀU & CƠ QUAN QUẢN LÝ (SSDL, VARANS)              |
|  - Bộ số liệu trường liều thực nghiệm đạt chuẩn ISO 4037 / TCVN 7042          |
|  - Nâng cao năng lực pháp định trong công tác thanh tra bức xạ quốc gia       |
+-------------------------------------------------------------------------------+
|  CÁC NHÀ NGHIÊN CỨU & DOSIMETRIST QUỐC TẾ                                     |
|  - Cung cấp dữ liệu thực nghiệm chính xác về hệ số suy giảm của chì và buồng A5|
|  - Làm dữ liệu kiểm chứng (benchmark) cho các thuật toán mô phỏng Monte Carlo |
+-------------------------------------------------------------------------------+

Câu hỏi thường gặp

1. Yêu cầu kỹ thuật cốt lõi để triển khai một phòng chuẩn gamma cấp 2 (SSDL) là gì?

Cần đáp ứng diện tích phòng tối thiểu $4\text{ m} \times 4\text{ m} \times 3\text{ m}$ với tường che chắn bức xạ bằng bê tông nặng hoặc chì (đảm bảo suất liều ngoài phòng $< 0,5\ \mu\text{Sv/h}$); hệ thống cửa khóa liên động an toàn (Interlock); hệ thống điều hòa kiểm soát nhiệt độ ($20 - 23^\circ\text{C}$), áp suất và độ ẩm ($< 70%$); hệ thống nguồn chuẩn có chứng chỉ PSDL và buồng ion hóa có liên kết chuẩn đo lường quốc tế.

2. Tại sao phải hiệu chỉnh hệ số nhiệt độ - áp suất ($K_{TP}$) trong phép đo buồng ion hóa?

Buồng ion hóa Exradin A5 là buồng thông khí (vented chamber), khối lượng không khí chứa trong thể tích $100\text{ cm}^3$ thay đổi theo mật độ môi trường (tỷ lệ thuận với áp suất $P$ và tỷ lệ nghịch với nhiệt độ tuyệt đối $T$). Do hệ số chuẩn $N_K$ được xác lập ở điều kiện chuẩn ($20^\circ\text{C}, 1013,25\text{ hPa}$), nếu không nhân với hệ số $K_{TP}$, kết quả suất liều có thể bị sai lệch từ $3% - 8%$.

3. Dòng điện rò ($I_{\text{leak}}$) của hệ thống được kiểm soát như thế nào để không làm sai lệch kết quả?

Quy trình yêu cầu sấy ấm điện kế SuperMAX 15 phút và ổn định kết nối với buồng A5 trong 10 phút trước khi đo. Thực hiện phép đo không có bức xạ (Zero reading) để xác định $I_{\text{leak}}$. Dòng rò bắt buộc phải nhỏ hơn $0,1%$ giá trị dòng đo có bức xạ. Nếu vượt quá, phải hiệu chỉnh trừ nền hoặc kiểm tra vệ sinh cáp Triaxial chống ẩm.

4. Sự khác biệt căn bản giữa suất Air Kerma ($\dot{K}_a$) và suất tương đương liều môi trường ($\dot{H}^*(10)$) là gì?

Air Kerma ($\dot{K}_a$, đơn vị $\text{Gy}$) là đại lượng vật lý thuần túy đo tổng động năng ban đầu của tất cả các hạt mang điện được giải phóng trong một đơn vị khối lượng không khí. Trong khi đó, $\dot{H}^*(10)$ (đơn vị $\text{Sv}$) là đại lượng bảo vệ bức xạ quy ước, được tính bằng cách nhân $\dot{K}_a$ với hệ số chuyển đổi $h$ ($h=1,20\text{ Sv/Gy}$ đối với $^{137}\text{Cs}$), mô phỏng liều hấp thụ tương đương ở độ sâu $10\text{ mm}$ trong mô người (quả cầu ICRU).

5. Chi phí đầu tư và thời gian thu hồi vốn (ROI) của hệ chuẩn G10 như thế nào?

Tổng mức đầu tư cho hệ chuẩn G10 đồng bộ (bao gồm nguồn $^{137}\text{Cs}$, trụ che chắn chì, buồng A5, điện kế SuperMAX và hệ thống an toàn) dao động từ 150.000 - 250.000 USD. Với công suất kiểm chuẩn trung bình 2.000 - 3.000 thiết bị/năm cho các cơ sở y tế và công nghiệp, thời gian thu hồi vốn ước tính từ 3,5 - 5 năm, đồng thời mang lại giá trị an toàn bức xạ vô hình to lớn cho xã hội.

Kết luận

Khóa luận tốt nghiệp "Khảo sát trường chuẩn liều gamma trên hệ chuẩn gamma tại Trung tâm Hạt nhân Thành phố Hồ Chí Minh" đã hoàn thành xuất sắc các mục tiêu khoa học và ứng dụng thực tiễn:

  • Thành tựu kỹ thuật cốt lõi: Khảo sát toàn diện và xây dựng thành công bộ dữ liệu chuẩn về suất Air Kerma $\dot{K}_a$, Suất tương đương liều môi trường $\dot{H}^*(10)$, và Suất tương đương liều cá nhân $\dot{H}_p(10)$ cho nguồn gamma $^{137}\text{Cs}$ ($662\text{ keV}$, $26,973\text{ Ci}$) trên dải khoảng cách từ 1,0 m đến 3,0 m qua 4 cấu hình che chắn (Nguồn trần, Tấm chắn X4, X10, X25).
  • Độ tin cậy đo lường: Xác lập ngân sách độ không đảm bảo đo chi tiết theo chuẩn quốc tế GUM với độ không đảm bảo mở rộng $U(H)$ chỉ từ $2,14%$ đến $3,12%$ ($k=2$, mức tin cậy $95%$), đáp ứng đầy đủ các tiêu chuẩn khắt khe của TCVN 7042-1:2008 và ISO 4037.
  • Giá trị thực tiễn: Đặt nền móng kỹ thuật vững chắc để đưa Phòng chuẩn đo liều cấp 2 (SSDL) tại Trung tâm Hạt nhân TP.HCM vào vận hành thương mại và nghiên cứu, phục vụ đắc lực cho công tác kiểm định thiết bị bức xạ y tế và công nghiệp trên phạm vi cả nước.

Khám phá thêm: Quý độc giả, sinh viên và các kỹ sư an toàn bức xạ quan tâm đến quy trình đo liều chuẩn hoặc có nhu cầu hợp tác kiểm chuẩn thiết bị đo bức xạ có thể liên hệ Bộ phận An toàn Bức xạ - Trung tâm Hạt nhân TP.HCM hoặc Khoa Vật lý, Trường Đại học Sư phạm TP.HCM để tiếp cận các báo cáo kỹ thuật chuyên sâu và tài liệu thực nghiệm liên quan.