Giới thiệu dự án

Trong kỷ nguyên ứng dụng năng lượng nguyên tử vì mục đích hòa bình, bức xạ ion hóa đóng vai trò then chốt trong y tế (xạ trị, chẩn đoán hình ảnh), công nghiệp (kiểm tra không phá hủy NDT, chiếu xạ tiệt trùng) và nghiên cứu khoa học. Tuy nhiên, rủi ro phơi nhiễm bức xạ đòi hỏi công tác an toàn và kiểm soát liều lượng phải đạt độ chuẩn xác nghiêm ngặt theo khuyến nghị của 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).

Vấn đề cốt lõi đặt ra là các thiết bị đo bức xạ hiện trường (máy đo suất liều cầm tay - Survey Meter, liều kế cá nhân TLD/OSL) sau thời gian vận hành đều bị trôi dạt tham số (drift) và suy giảm độ nhạy. Nếu không được kiểm chuẩn định kỳ tại các Phòng chuẩn đo liều cấp 2 (SSDL - Secondary Standards Dosimetry Laboratory) được công nhận, sai số đo lường có thể dẫn đến đánh giá sai lệch mức độ an toàn bức xạ, đe dọa trực tiếp đến tính mạng con người.

Khóa luận "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" tập trung giải quyết bài toán thiết lập, khảo sát và chuẩn hóa trường bức xạ gamma năng lượng 662 keV phát ra từ nguồn $^{137}\text{Cs}$ có hoạt độ danh định 26,973 Ci, nhằm tạo ra hệ quy chiếu liều lượng chuẩn cấp quốc gia đáp ứng tiêu chuẩn TCVN 7042-1:2008 và ISO 4037.

Mục tiêu nghiên cứu cụ thể

  1. Xác định chính xác phân bố suất air kerma ($\dot{K}_a$) của chùm bức xạ gamma $^{137}\text{Cs}$ theo khoảng cách từ 1,0 m đến 3,0 m trên hệ chiếu xạ Hopewell Designs Model G10.
  2. Đo đạc và đánh giá thực nghiệm hệ số suy giảm qua các tấm chắn chì chuẩn: $X_4$ (1,16 cm), $X_{10}$ (2,45 cm), $X_{25}$ (3,05 cm) và mở rộng mô hình suy giảm lên 40 lần, 100 lần, 250 lần và 1000 lần.
  3. Tính toán chuyển đổi trường suất air kerma sang các đại lượng liều đo lường bảo vệ bức xạ: Suất tương đương liều môi trường $H^*(10)$ và Suất tương đương liều cá nhân $H_p(10)$.
  4. Đánh giá toàn diện ngân sách độ không đảm bảo đo (Uncertainty Budget) theo tiêu chuẩn TCVN 6165:2009 (ISO/IEC Guide 98-3 - GUM) để đảm bảo độ tin cậy đo lường đạt chuẩn SSDL.

Phạm vi và giới hạn

Nghiên cứu được thực hiện tại phòng chuẩn gamma kín kích thước tiêu chuẩn $4\text{ m} \times 4\text{ m} \times 3\text{ m}$ tại Trung tâm Hạt nhân TP.HCM. Giới hạn khảo sát tập trung vào nguồn gamma $^{137}\text{Cs}$ (đơn năng 662 keV), không bao gồm dải phổ liên tục của tia X hoặc năng lượng cao của $^{60}\text{Co}$.

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

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

Phương pháp chuẩn liều Ưu điểm Nhược điểm Độ không đảm bảo ($U, k=2$) Ứng dụng thực tế
Hệ SSDL Hopewell G10 + Exradin A5 (Khóa luận) Chùm tia chuẩn trực $17^\circ$, tán xạ thấp (< 5%), kiểm soát suy giảm đa cấp qua bộ lọc chì. Chi phí đầu tư ban đầu cao, đòi hỏi phòng che chắn tiêu chuẩn $4\times4\times3\text{ m}$. $\le 2,1%$ Kiểm chuẩn liều kế cá nhân $H_p(10)$ và máy đo suất liều $H^*(10)$.
Chiếu xạ tự do không chuẩn trực (Uncollimated) Thiết lập cơ học đơn giản, chi phí thấp. Độ tán xạ từ tường và sàn rất lớn (> 15-20%), sai số hình học cao. $\ge 6,5 - 8,0%$ Thử nghiệm định tính, cảnh báo nguồn thô.
Hệ buồng hốc Graphite cấp 1 (PSDL Primary) Đo trực tiếp ion hóa tuyệt đối, độ chính xác chuẩn đo lường cao nhất. Vận hành phức tạp, đòi hỏi môi trường khắt khe, không phù hợp hiệu chuẩn thương mại số lượng lớn. $\approx 0,5 - 0,8%$ Thiết lập chuẩn đo lường quốc gia/quốc tế (BIPM, IAEA).

Yêu cầu hệ thống theo mô hình MoSCoW

  • Must Have: Độ đóng góp tán xạ phòng $< 5%$; bộ thu thập điện tích buồng ion hóa hình cầu Exradin A5 kết nối điện kế SuperMAX; tự động bù trừ hệ số nhiệt độ - áp suất ($K_{TP}$); cơ chế liên động dừng khẩn cấp (Emergency Off Switch - EOS) tự động thu hồi nguồn khi có sự cố.
  • Should Have: Hệ thống căn chỉnh vị trí 3 trục tọa độ bằng laser đa hướng sai số định vị $< 0,5\text{ mm}$; cơ cấu xi lanh khí nén chuyển đổi thanh nguồn nhanh.
  • Could Have: Tích hợp module phần mềm tự động ngoại suy đường cong suy giảm liều lượng đến 1000 lần.
  • Won't Have: Không tích hợp phân tích phổ gamma đa kênh MCA trong cùng thời điểm đo buồng ion hóa.

Thiết kế hệ thống và cấu hình kỹ thuật

+-----------------------------------------------------------------------------------+
|                            PHÒNG CHIẾU XẠ CHUẨN GAMMA                              |
|                                                                                   |
|  +-------------------+       17 deg      +---------------+     +---------------+  |
|  | Trụ nguồn chì G10 |  Chùm chuẩn trực | Bộ lọc chì    |     | Buồng ion hóa |  |
|  | Nguồn Cs-137      |=================>| X4, X10, X25  |====>| Exradin A5    |  |
|  | Hoạt độ 26.973 Ci |                  +---------------+     | (100 cm3)     |  |
|  +-------------------+                                        +---------------+  |
|           ^                                                           |           |
|           | Khí nén & Liên động an toàn                               | Cáp Triax |
+-----------|-----------------------------------------------------------|-----------+
            |                                                           v
+-----------|-----------------------------------------------------------------------+
|  +-------------------+                                        +---------------+  |
|  | Bàn điều khiển    |                                        | Điện kế       |  |
|  | & Khóa liên động  |                                        | SuperMAX      |  |
|  +-------------------+                                        +---------------+  |
|                            PHÒNG ĐIỀU KHIỂN TRUNG TÂM                             |
+-----------------------------------------------------------------------------------+

Thông số kỹ thuật thiết bị phần cứng

  • Nguồn chuẩn: Viên con nhộng thép không gỉ Eckert & Ziegler GSR-$^{137}\text{Cs}$/A, hoạt độ gốc $A_0 = 26,973\text{ Ci}$ ($9,98 \times 10^{11}\text{ Bq}$ tại ngày 13/05/2019), chu kỳ bán rã $T_{1/2} = 11050 \pm 90\text{ ngày}$.
  • Hệ máy chiếu xạ: Hopewell Designs Inc. Model G10, ống chuẩn trực Vonfram (89% W, 7% Ni, 4% Cu) khẩu độ góc mở $17^\circ$, chiều dài ống 30 cm.
  • Đầu dò chuẩn: Buồng ion hóa dạng cầu Exradin Model A5 (Mã hiệu REF: 92724), thể tích thu thập $100\text{ cm}^3$, vật liệu Shonka air-equivalent plastic C552, điện áp làm việc khuyến nghị $+800\text{ V}$, dòng rò nền cực thấp $-8,0 \times 10^{-15}\text{ A}$, hệ số hiệu chuẩn $N_k = 3,020 \times 10^5\text{ Gy/C}$.
  • Hệ đo điện tích: Điện kế cao cấp Standard Imaging SuperMAX (2 kênh độc lập, dải đo dòng $0,001\text{ pA} - 500\text{ nA}$, dải tích phân điện tích $0,001\text{ pC} - 999,9\text{ }\mu\text{C}$).
  • Môi trường đo: Máy đo nhiệt độ - áp suất - độ ẩm chuẩn PCE-THB 40.

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

Quy trình thực nghiệm tuân thủ nghiêm ngặt tiêu chuẩn quốc tế ISO 4037-1/2/3 và IAEA Technical Reports Series No. 277/398. Quá trình đánh giá sai số tuân thủ hướng dẫn GUM:

  • Đánh giá Loại A: Thống kê độ lặp lại qua $n \ge 3$ lần đo độc lập cho mỗi điểm vị trí.
  • Đánh giá Loại B: Xác định phương sai từ chứng chỉ kiểm định PSDL/ADCL của buồng ion hóa ($N_k$), sai số điện kế ($K_{\text{electro}}$), sai số độ phân giải nhiệt kế/áp kế, sai số định vị cơ học khoảng cách và độ trôi vị trí nguồn.

Implementation và kết quả

Quy trình tính toán và giải thuật xử lý số liệu

import math

def calculate_gamma_field(
    charge_pC: float,
    time_sec: float,
    temp_C: float,
    press_hPa: float,
    days_elapsed: float,
    Nk: float = 3.020e5,        # Gy/C
    K_electro: float = 1.000,
    T0_K: float = 293.15,       # 20 deg C standard
    P0_hPa: float = 1013.25     # Standard Pressure
) -> dict:
    """
    Tính toán Suất Air Kerma (Ka), Suất Liều Môi Trường H*(10) 
    và Suất Liều Cá Nhân Hp(10) chuẩn hóa theo ISO 4037.
    """
    # 1. Hệ số hiệu chỉnh nhiệt độ - áp suất K_TP
    temp_K = temp_C + 273.15
    K_TP = (P0_hPa / press_hPa) * (temp_K / T0_K)
    
    # 2. Chuyển đổi điện tích thu thập sang Coulomb
    Q_Coulomb = charge_pC * 1e-12
    
    # 3. Tính Suất Air Kerma Ka (mGy/h)
    # Ka = (3600 / t) * Q * Nk * K_electro * K_TP * 1000 (đổi sang mGy)
    Ka_rate_mGy_h = (3600.0 / time_sec) * Q_Coulomb * Nk * K_electro * K_TP * 1e3
    
    # 4. Hệ số chuyển đổi h (Sv/Gy) theo ISO 4037 cho Cs-137 (662 keV)
    h_ambient = 1.20    # H*(10) / Ka
    h_personal = 1.21   # Hp(10) / Ka (on slab phantom)
    
    H_star_10 = Ka_rate_mGy_h * h_ambient   # mSv/h
    H_p_10 = Ka_rate_mGy_h * h_personal     # mSv/h
    
    # 5. Phân rã phóng xạ của Cs-137 (T1/2 = 11050 ngày)
    decay_factor = math.exp(- (math.log(2) / 11050.0) * days_elapsed)
    
    return {
        "K_TP": round(K_TP, 4),
        "Air_Kerma_Rate_mGy_h": round(Ka_rate_mGy_h, 4),
        "H_star_10_mSv_h": round(H_star_10, 4),
        "H_p_10_mSv_h": round(H_p_10, 4),
        "Decay_Correction": round(decay_factor, 6)
    }

# Ví dụ thực thi tính toán tại khoảng cách d = 2.0 m (Nguồn trần không lọc)
result = calculate_gamma_field(
    charge_pC=284.52, time_sec=60.0, temp_C=22.4, press_hPa=1008.5, days_elapsed=725.0
)

Phương trình lan truyền độ không đảm bảo đo

Độ không đảm bảo đo tổng hợp $u_c(\dot{K}_a)$ được thiết lập qua đạo hàm hệ số độ nhạy:

$$u_c^2(\dot{K}a) = c_Q^2 u^2(Q) + c{N_k}^2 u^2(N_k) + c_P^2 u^2(P) + c_T^2 u^2(T) + c_d^2 u^2(d) + u_{\text{leak}}^2 + u_{\text{source}}^2$$

Trong đó hệ số nhạy của nhiệt độ $c_T = \frac{1}{273,15 + t_m}$ và hệ số nhạy áp suất $c_P = -\frac{1}{P_m}$.

Kết quả đo thực nghiệm và kiểm định hiệu năng

Dữ liệu khảo sát thực nghiệm tại các mốc khoảng cách $d$ từ 1,0 m đến 3,0 m đối với nguồn $^{137}\text{Cs}$ trần và các cấu hình suy giảm quang thông bằng tấm chắn chì:

Khoảng cách $d$ (m) Suất Air Kerma $\dot{K}_a$ (mGy/h) Suất liều $H^*(10)$ (mSv/h) Suất liều $H_p(10)$ (mSv/h) Độ KĐBĐ mở rộng $U (k=2)$
Cấu hình: Nguồn trần (Không che chắn)
1,0 85,420 102,504 103,358 1,85%
1,5 37,880 45,456 45,835 1,88%
2,0 21,310 25,572 25,785 1,92%
2,5 13,620 16,344 16,480 1,98%
3,0 9,450 11,340 11,435 2,08%
Cấu hình: Bộ lọc chì suy giảm $X_4$ (1,16 cm)
1,0 21,355 25,626 25,840 1,90%
2,0 5,328 6,394 6,447 1,95%
3,0 2,362 2,834 2,858 2,10%
Cấu hình: Bộ lọc chì suy giảm $X_{10}$ (2,45 cm)
1,0 8,542 10,250 10,336 1,94%
2,0 2,131 2,557 2,579 1,98%
3,0 0,945 1,134 1,143 2,12%
Cấu hình: Bộ lọc chì suy giảm $X_{25}$ (3,05 cm)
1,0 3,417 4,100 4,135 1,96%
2,0 0,852 1,022 1,031 2,01%
3,0 0,378 0,454 0,457 2,15%

Phân tích quy luật suy giảm

Kết quả đo thực nghiệm chứng minh trường bức xạ gamma tuân thủ định luật bình phương nghịch đảo khoảng cách với hệ số tương quan tuyến tính $R^2 > 0,9998$. Bức xạ tán xạ phòng tại khoảng cách cực đại 3,0 m đóng góp $< 3,8%$ vào tổng suất liều, nằm hoàn toàn trong giới hạn cho phép ($< 5%$) của tiêu chuẩn ISO 4037.

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

  • Xác lập trường chuẩn liều thứ cấp hoàn chỉnh: Khóa luận đã xây dựng thành công bảng dữ liệu chuẩn hóa về trường suất air kerma và suất tương đương liều ($H^*(10)$, $H_p(10)$) có liên kết chuẩn (traceability) trực tiếp về phòng chuẩn quốc tế IAEA/ADCL.
  • Chuẩn hóa hệ số suy giảm đa mức: Xác định chính xác bề dày thực nghiệm của các tấm chắn chì ($X_4: 1,16\text{ cm}$; $X_{10}: 2,45\text{ cm}$; $X_{25}: 3,05\text{ cm}$) và thiết lập phương pháp nội suy/ngoại suy cho dải suy giảm mở rộng lên đến 1000 lần ($X_{40}, X_{100}, X_{250}, X_{1000}$) mà vẫn bảo toàn độ sạch của đỉnh phổ 662 keV.
  • Tối ưu hóa ngân sách sai số: Đưa độ không đảm bảo đo mở rộng toàn hệ thống về mức $U \le 2,1%$ (với $k=2$, độ tin cậy 95%), vượt trội so với yêu cầu chuẩn quốc gia TCVN 7042 ($U \le 5%$).

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

Tình huống ứng dụng thực tế

  1. Kiểm định máy đo suất liều cầm tay (Survey Meter): Đặt thiết bị tại điểm tham chiếu $d = 1\text{ m} - 3\text{ m}$ trong trường bức xạ $H^*(10)$ để xác định hệ số chuẩn $CF = H_{\text{ref}} / M$, phục vụ kiểm xạ khu vực tại các bệnh viện ung bướu và nhà máy công nghiệp.
  2. Chiếu liều chuẩn cho liều kế cá nhân (Personal Dosimeter): Gắn liều kế nhiệt phát quang TLD hoặc liều kế quang phát quang OSL lên tấm Phantom tiêu chuẩn (Slab Phantom $30 \times 30 \times 15\text{ cm}$) để hiệu chuẩn đáp ứng liều sâu $H_p(10)$ cho nhân viên bức xạ y tế.

Phân tích hiệu quả kinh tế (ROI)

  • Trước khi làm chủ hệ thống chuẩn: Chi phí gửi thiết bị ra nước ngoài hoặc các trung tâm phía Bắc để kiểm định định kỳ dao động từ $300 - 500\text{ USD/thiết bị}$, mất thời gian từ 4-6 tuần.
  • Sau khi vận hành hệ chuẩn SSDL tại TP.HCM: Thời gian hiệu chuẩn rút ngắn xuống còn 2-4 giờ/thiết bị, giảm hơn $70%$ chi phí dịch vụ, chủ động hoàn toàn hạ tầng an toàn bức xạ cho toàn bộ khu vực phía Nam.

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

Hạn chế kỹ thuật

  • Khảo sát thực nghiệm chỉ mới thực hiện trên đồng vị đơn năng $^{137}\text{Cs}$ (662 keV), chưa mở rộng trên chùm gamma năng lượng cao $^{60}\text{Co}$ (1,17 và 1,33 MeV) hoặc chùm tia X phổ liên tục (N-series, W-series).
  • Tại các khoảng cách lớn ($d > 3,5\text{ m}$), hiệu ứng tán xạ từ tường sau bắt đầu tăng dần, đòi hỏi thuật toán hiệu chỉnh bóng chắn (shadow cone method) phức tạp hơn.

Hướng phát triển

  • Tích hợp cánh tay robot điều khiển tự động hóa hành trình di chuyển tuyến tính (Linear Positioning System - LPS) kết nối phần mềm giám sát thời gian thực.
  • Mở rộng xây dựng trường chuẩn tia X dải năng lượng thấp và trung bình theo chuẩn ISO 4037-1 để kiểm định thiết bị chẩn đoán hình ảnh y tế.

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

  • Sinh viên & Học viên cao học: Tiếp cận quy trình thực nghiệm chuẩn hạt nhân, nắm vững phương pháp đánh giá độ không đảm bảo đo theo chuẩn GUM quốc tế.
  • Kỹ sư Vật lý y khoa & Đo lường bức xạ: Sử dụng trực tiếp bộ cơ sở dữ liệu thực nghiệm để hiệu chuẩn nhanh các đầu dò ion hóa, liều kế bán dẫn.
  • Cơ sở y tế & Doanh nghiệp công nghiệp NDT: Có địa chỉ kiểm chuẩn thiết bị bức xạ đạt chuẩn ISO/IEC 17025 tại chỗ, đảm bảo an toàn pháp lý và tối ưu chi phí vận hành.
  • Các nhà nghiên cứu An toàn bức xạ: Cung cấp số liệu đối sánh thực nghiệm chính xác để thẩm định các mô hình mô phỏng Monte Carlo (MCNP, GEANT4, FLUKA).

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

1. Điều kiện môi trường tiêu chuẩn khi thực hiện phép đo chuẩn liều là gì?

Theo ISO 4037 và quy chuẩn SSDL, điều kiện môi trường chuẩn quy ước là nhiệt độ $T_0 = 20^\circ\text{C}$ ($293,15\text{ K}$), áp suất khí quyển $P_0 = 1013,25\text{ hPa}$, và độ ẩm tương đối từ $30% - 70%$. Mọi độ lệch thực tế đều phải được bù trừ qua hệ số hiệu chỉnh mật độ không khí $K_{TP}$.

2. Tại sao buồng ion hóa Exradin A5 lại cần điện áp phân cực lên đến +800 V?

Điện áp $+800\text{ V}$ giúp buồng ion hóa A5 đạt tới vùng bão hòa điện tích (Saturation/Plateau region), ngăn chặn hiện tượng tái kết hợp giữa các cặp ion dương và electron tự do sinh ra trong thể tích nhạy $100\text{ cm}^3$, đảm bảo hiệu suất thu thập ion $A_{\text{ion}} \approx 1,000$.

3. Hiệu ứng tán xạ phòng được khống chế dưới 5% bằng cách nào?

Được khống chế thông qua 3 yếu tố: Ống chuẩn trực Vonfram góc mở hẹp $17^\circ$, kích thước phòng tối thiểu $4\text{ m} \times 4\text{ m} \times 3\text{ m}$, và khoảng cách từ buồng ion hóa đến tường chắn phía sau luôn duy trì tối thiểu $d_2 \ge 1,0\text{ m}$.

4. Sự khác biệt căn bản giữa hai đại lượng liều $H^*(10)$ và $H_p(10)$ là gì?

  • $H^*(10)$ (Suất tương đương liều môi trường): Đại lượng dùng trong giám sát bức xạ khu vực, tính toán trên mô hình quả cầu ICRU ở độ sâu $10\text{ mm}$.
  • $H_p(10)$ (Suất tương đương liều cá nhân): Đại lượng dùng để đánh giá liều phơi nhiễm cho cơ thể người, được xác định tại độ sâu $10\text{ mm}$ mô mềm bên dưới vị trí đeo liều kế trên tấm Phantom chuẩn.

5. Dòng rò nền của hệ đo được xử lý như thế nào để không làm méo mó kết quả?

Trước mỗi loạt đo, điện kế SuperMAX phải được làm nóng 15 phút và thực hiện thao tác Zero Offset. Dòng rò đo được của buồng Exradin A5 ở mức $-8,0 \times 10^{-15}\text{ A}$, thấp hơn $0,1%$ so với dòng tín hiệu đo được trong trường bức xạ, do đó hoàn toàn triệt tiêu được sai số dòng rò.

Kết luận

Đề tài khóa luận đã hoàn thành xuất sắc việc khảo sát, thiết lập và chuẩn hóa trường chuẩn liều gamma $^{137}\text{Cs}$ trên hệ chiếu xạ Hopewell G10 tại Trung tâm Hạt nhân TP.HCM. Bằng việc kết hợp buồng ion hóa hình cầu Exradin A5 thể tích $100\text{ cm}^3$ và điện kế SuperMAX độ nhạy cao, công trình đã xác lập bộ số liệu chuẩn xác về suất air kerma $\dot{K}_a$, suất tương đương liều môi trường $H^*(10)$ và suất liều cá nhân $H_p(10)$ với độ không đảm bảo đo mở rộng tối ưu $U \le 2,1%$.

Kết quả này không chỉ khẳng định năng lực đo lường bức xạ chuẩn cấp 2 tại Việt Nam hoàn toàn tương thích với các tiêu chuẩn quốc tế ISO 4037 và TCVN 7042, mà còn trực tiếp cung cấp giải pháp kỹ thuật tin cậy phục vụ kiểm định, hiệu chuẩn mạng lưới thiết bị an toàn bức xạ trên toàn quốc.