Giới thiệu dự án

Bối cảnh và thực trạng ngành quan trắc phóng xạ

Bức xạ ion hóa môi trường luôn hiện diện xung quanh đời sống con người, bao gồm bức xạ vũ trụ và các đồng vị phóng xạ tự nhiên (NORM - Naturally Occurring Radioactive Materials) tồn tại trong thạch quyển, thủy quyển và khí quyển. Theo báo cáo từ Ủy ban Khoa học của Liên Hợp Quốc về Tác động của Bức xạ Nguyên tử (UNSCEAR), liều bức xạ hiệu dụng trung bình toàn cầu mà một người tiếp nhận từ môi trường tự nhiên xấp xỉ 2,4 mSv/năm. Trong đó, đóng góp chính đến từ các chuỗi phân rã phóng xạ tự nhiên: chuỗi Uranium ($^{238}\text{U}$), chuỗi Thorium ($^{232}\text{Th}$), chuỗi Actinium ($^{235}\text{U}$) và đồng vị đơn lẻ Kali-40 ($^{40}\text{K}$). Chiếu xạ ngoài từ đất đá đóng góp khoảng 0,45 mSv/năm, trong khi nước và thực phẩm đóng góp khoảng 0,25 mSv/năm.

                  +-------------------------------------------------------+
                  |         BỨC XẠ MÔI TRƯỜNG TỰ NHIÊN (~2.4 mSv/năm)      |
                  +---------------------------+---------------------------+
                                              |
            +---------------------------------+---------------------------------+
            |                                 |                                 |
            v                                 v                                 v
+-----------------------+         +-----------------------+         +-----------------------+
|  Chuỗi Uranium-238    |         |  Chuỗi Thorium-232    |         |     Đồng vị NORM      |
|  - Ra-226 (1602 năm)  |         |  - Ra-228 (6.7 năm)   |         |  - K-40 (1.25x10^9 y) |
|  - Rn-222 (3.82 ngày) |         |  - Ac-228 (6.15 giờ)  |         |  - Bức xạ vũ trụ      |
|  - Pb-214, Bi-214     |         |  - Pb-212, Tl-208     |         |  (Chiếu trong/ngoài)  |
+-----------------------+         +-----------------------+         +-----------------------+

Vấn đề nghiên cứu và các điểm nghẽn kỹ thuật (Problem Statement)

Trong kiểm soát an toàn bức xạ và đánh giá phông phóng xạ môi trường, việc xác định chính xác suất liều chiếu (Exposure Dose Rate - $\dot{X}$) là yêu cầu bắt buộc. Tuy nhiên, các phương pháp hiện hành gặp phải những rào cản lớn:

  • Thiết bị đo liều tích phân truyền thống (Survey meters, buồng ion hóa, ống đếm GM): Chỉ ghi nhận tổng suất liều mà không cung cấp thông tin năng lượng photon, dẫn đến sai số đáp ứng phổ lớn khi phổ năng lượng môi trường thay đổi phức tạp.
  • Đầu dò bán dẫn siêu tinh khiết HPGe (High-Purity Germanium): Có độ phân giải năng lượng cực cao (FWHM $< 0,2%$ tại 1332 keV) giúp định danh chính xác từng đồng vị phóng xạ, nhưng đòi hỏi làm lạnh liên tục bằng Nitơ lỏng ($77\text{ K}$), kích thước cồng kềnh, chi phí bảo dưỡng cao và không tối ưu cho mạng lưới quan trắc dã ngoại cơ động.
  • Đầu dò nhấp nháy $\text{NaI(Tl)}$ truyền thống: Dù có hiệu suất ghi nhận cao và vận hành được ở nhiệt độ phòng, nhưng độ phân giải năng lượng thấp dẫn đến hiện tượng chồng chập đỉnh phổ gamma (overlapping peaks), gây khó khăn cho việc định danh và đánh giá suất liều trực tiếp từ phổ xung.

Mục tiêu của đồ án

Đồ án tập trung giải quyết các mục tiêu khoa học và ứng dụng cụ thể:

  1. Thiết lập hệ đo thực nghiệm phổ kế gamma: Tiến hành đo phổ phông phóng xạ môi trường liên tục trong 24 giờ (86.400 giây) bằng hai hệ đầu dò song song: detector bán dẫn HPGe GC2518 (Canberra Inc.) và detector nhấp nháy $\text{NaI(Tl)}$ kích thước $7,6\text{ cm} \times 7,6\text{ cm}$ ($3'' \times 3''$) tích hợp trên hệ InSpector™ 1000.
  2. Định danh phổ và bóc tách đỉnh phổ: Ứng dụng phần mềm chuyên dụng Colegram kết hợp cơ sở dữ liệu phân rã hạt nhân Laraweb và Table of Isotopes để phân tích, bóc tách các đỉnh phổ chồng chập trên phổ $\text{NaI(Tl)}$ dựa trên đối chiếu phổ chuẩn độ phân giải cao HPGe.
  3. Mô hình hóa hàm chuyển đổi suất liều $G(E)$: Khớp hàm toán học đa thức bậc 9 cho hệ số chuyển đổi phổ năng lượng sang suất liều chiếu $G(E)$ đối với đầu dò $\text{NaI(Tl)}$ $7,6\text{ cm} \times 7,6\text{ cm}$ theo hai cấu hình hình học photon tới (song song và vuông góc với trục đối xứng của tinh thể).
  4. Xác định suất liều chiếu môi trường tổng: Tính toán tích phân suất liều chiếu môi trường ($\mu\text{R/h}$) và đánh giá mức độ tin cậy của phương pháp đo phổ liều lượng (Dose Rate Spectroscopy).

Giải pháp kỹ thuật và kết quả kỳ vọng

  • Giải pháp: Ứng dụng phương pháp chuyển đổi ma trận đáp ứng nghịch đảo phổ năng lượng thông qua hàm $G(E)$ (đơn vị: $\mu\text{R}\cdot\text{h}^{-1}/\text{cps}$), cho phép chuyển đổi trực tiếp tốc độ đếm $n(E)$ tại từng kênh năng lượng thành suất liều chiếu theo công thức giải tích: $$\dot{X} = \int_{E_{\min}}^{E_{\max}} n(E) G(E) , dE$$
  • Chỉ số kỳ vọng: Xác định suất liều chiếu môi trường phòng thí nghiệm đạt độ tin cậy cao, hệ số tương quan làm khớp hàm $R^2 \ge 0,9999$, đồng thời nhận diện đầy đủ các nhóm đồng vị phát gamma tự nhiên thuộc chuỗi $^{238}\text{U}$, $^{232}\text{Th}$ và $^{40}\text{K}$.

Phạm vi và giới hạn nghiên cứu

  • Phạm vi không gian: Thực hiện tại Phòng thí nghiệm Kỹ thuật Hạt nhân, Trường Đại học Khoa học Tự nhiên - ĐHQG TP.HCM.
  • Phạm vi năng lượng: Khảo sát dải phổ gamma từ $50\text{ keV}$ đến $3000\text{ keV}$.
  • Giới hạn: Khóa luận tập trung tính toán suất liều chiếu tổng dựa trên phân bố phổ năng lượng, chưa đi sâu vào bóc tách suất liều riêng phần cho từng đồng vị cá lẻ trong điều kiện phông biến thiên nhanh.

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

Phân tích hiện trạng công nghệ đo bức xạ

Tiêu chí Buồng Ion hóa / Survey Meter GM Hệ phổ kế HPGe (GC2518) Hệ phổ kế NaI(Tl) 7,6cm x 7,6cm + G(E)
Bản chất đo Tích phân dòng ion hóa / đếm xung thô Phổ kế bán dẫn siêu tinh khiết Phổ kế nhấp nháy + Chuyển đổi $G(E)$
Độ phân giải năng lượng Không có thông tin năng lượng Rất cao ($< 2\text{ keV}$ tại $1332\text{ keV}$) Trung bình (~$7%$ tại $662\text{ keV}$)
Hiệu suất ghi nhận Thấp đến trung bình Trung bình (tinh thể nhỏ) Rất cao (tinh thể tích tích $344\text{ cm}^3$)
Điều kiện vận hành Đơn giản, pin nguồn Phức tạp, bắt buộc làm lạnh $\text{LN}_2$ ($77\text{ K}$) Linh động, hoạt động ở nhiệt độ phòng
Khả năng định danh NORM Không khả thi Rất chính xác (tách đỉnh $< 1\text{ keV}$) Tốt (kết hợp thuật toán deconvolution)
Chi phí đầu tư & bảo trì Thấp ($< 2.000\text{ USD}$) Rất cao ($> 60.000\text{ USD}$ + chi phí $\text{LN}_2$) Vừa phải ($8.000 - 15.000\text{ USD}$)

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

  • Must have (Bắt buộc có):
    • Hệ thu nhận xung MCA đa kênh (tối thiểu 1024 kênh cho NaI và 8192 kênh cho HPGe).
    • Thuật toán làm khớp hàm đa thức $G(E)$ bậc cao đảm bảo sai số phần dư tối thiểu ($R^2 > 0,999$).
    • Module phân tách đỉnh phổ chồng chập cho các năng lượng đặc trưng ($^{214}\text{Pb}$, $^{214}\text{Bi}$, $^{208}\text{Tl}$, $^{228}\text{Ac}$).
  • Should have (Nên có):
    • Cơ chế bù trừ hướng tới của chùm bức xạ (hàm $G(E)$ theo phương song song và vuông góc).
    • Giao diện trực quan hóa phổ năng lượng và phân bố đóng góp liều lượng.
  • Could have (Có thể mở rộng):
    • Khả năng tự động tính toán suất liều riêng phần (Individual Dose Rate Spectroscopy).
  • Won't have (Chưa thực hiện):
    • Tự động bù trôi đỉnh phổ theo nhiệt độ thời gian thực bằng phần cứng analog.
+-----------------------------------------------------------------------------------+
|                     KIẾN TRÚC HỆ THỐNG ĐO VÀ PHÂN TÍCH SUẤT LIỀU                  |
+-----------------------------------------------------------------------------------+
  [ Photon Gamma Môi Trường (NORM: U-238, Th-232, K-40) ]
                            |
                            v
  +------------------------------------------------------+
  | Đầu dò NaI(Tl) 7.6cm x 7.6cm + Ống nhân quang (PMT)  |
  +------------------------------------------------------+
                            | (Tín hiệu xung điện tích)
                            v
  +------------------------------------------------------+
  | Bộ tiền khuếch đại (Preamplifier) & Định hình xung   |
  +------------------------------------------------------+
                            | (Tín hiệu điện thế tương tự)
                            v
  +------------------------------------------------------+
  | Khối phân tích đa kênh (Canberra InSpector 1000 MCA)  |
  +------------------------------------------------------+
                            | (Tệp dữ liệu phổ xung n(E))
                            v
  +------------------------------------------------------+
  | Phần mềm Colegram: Phân tách đỉnh & Khử nhiễu phông |
  +------------------------------------------------------+
                            | (Phổ năng lượng đã xử lý)
                            v
  +------------------------------------------------------+
  | Module Tích phân Số: n(E) x G(E) -> Suất liều X_dot  |
  +------------------------------------------------------+

Công nghệ và thiết bị sử dụng

  1. Đầu dò NaI(Tl) InSpector™ 1000 (Canberra Inc.):
    • Tinh thể nhấp nháy Thallium-doped Sodium Iodide hình trụ kích thước $\varnothing 76\text{ mm} \times 76\text{ mm}$.
    • Tích hợp ống nhân quang điện (Photomultiplier Tube - PMT) và mạch tiền khuếch đại điện tích.
  2. Đầu dò HPGe GC2518 (Canberra Inc.):
    • Tinh thể Germanium siêu tinh khiết đồng trục, đường kính $61,8\text{ mm}$, bề dày $35,8\text{ mm}$.
    • Bình chứa cryostat làm lạnh bằng Nitơ lỏng, kết nối MCA 8192 kênh (8k).
  3. Phần mềm phân tích phổ Colegram:
    • Công cụ chuyên dụng giải chập phổ gamma, khớp dạng đỉnh Gauss biến tính (Gaussian with tailing) và nền Compton phi tuyến.
  4. Cơ sở dữ liệu hạt nhân: Thư viện đồng vị phóng xạ Laraweb (CEA/LNHB) và bảng số liệu phân rã Table of Isotopes (8th Edition).

Implementation và kết quả thực nghiệm

Cơ sở lý thuyết và mô hình toán học

Tương tác của chùm bức xạ photon với môi trường được mô tả thông qua hệ số Kerma ($K$), liều hấp thụ ($D$), và liều chiếu ($X$). Kerma va chạm trong không khí $(K_c){\text{air}}$ liên hệ với thông lượng năng lượng $\psi$ thông qua hệ số hấp thụ năng lượng khối $(\mu{\text{en}}/\rho){\text{air}}$: $$(K_c){\text{air}} = \psi \cdot \left(\frac{\mu_{\text{en}}}{\rho}\right)_{E, \text{air}}$$

Mối liên hệ giữa suất liều chiếu $\dot{X}$ và thông lượng photon $\psi_j$ tại mức năng lượng $E_j$: $$\dot{X} = 65,8 \cdot \left(\frac{\mu_{\text{en}}}{\rho}\right)_{\text{air}} \cdot E_j \cdot \psi_j = \left(\frac{\dot{X}}{\psi_j}\right) \cdot \psi_j$$

Khi photon tương tác với detector, phổ năng lượng thực nghiệm $n(E)$ liên hệ với thông lượng thông qua ma trận đáp ứng $R_{ij}$. Hệ số chuyển đổi suất liều chiếu $G_j$ biểu diễn đáp ứng nghịch đảo: $$G_j = R_{ij}^{-1} \cdot \left(\frac{\dot{X}}{\psi_j}\right)$$

Khớp hàm đa thức hệ số chuyển đổi $G(E)$

Hệ số $G(E)$ của đầu dò $\text{NaI(Tl)} \ 7,6\text{ cm} \times 7,6\text{ cm}$ được xây dựng từ mô phỏng Monte Carlo (MCNP/Geant4) và được khớp bởi phương trình đa thức bậc 9 theo hàm $\log$:

$$\log[G(E)] = \sum_{k=0}^{9} a_k \cdot (\log E)^k$$

  1. Trường hợp photon tới song song với trục detector:

    • Hệ số hồi quy:
      • $a_0 = -1232,80706$
      • $a_1 = +5288,20962$
      • $a_2 = -9584,54726$
      • $a_3 = +9752,82964$
      • $a_4 = -6193,51923$
      • $a_5 = +2559,84742$
      • $a_6 = -691,24776$
      • $a_7 = +117,93858$
      • $a_8 = -11,56348$
      • $a_9 = +0,49737$
    • Hệ số tương quan: $R^2 = 0,99998$.

  2. Trường hợp photon tới vuông góc với trục detector:

    • Hệ số hồi quy:
      • $a_0 = -1603,32194$
      • $a_1 = +6794,83543$
      • $a_2 = -12203,37818$
      • $a_3 = +12318,43163$
      • $a_4 = -7758,24172$
      • $a_5 = +3176,31261$
      • $a_6 = -848,01750$
      • $a_7 = +142,69653$
      • $a_8 = -13,75718$
      • $a_9 = +0,57981$
    • Hệ số tương quan: $R^2 = 0,99996$.

Thuật toán tính toán suất liều chiếu tích phân (Python Implementation)

Đoạn mã xử lý tích phân phổ năng lượng và tính toán suất liều môi trường:

import numpy as np

def calculate_ge_factor(energy_kev: np.ndarray, orientation: str = 'parallel') -> np.ndarray:
    """
    Tính hệ số chuyển đổi suất liều G(E) cho đầu dò NaI(Tl) 7.6cm x 7.6cm.
    Năng lượng E tính bằng keV, G(E) tính bằng (uR/h)/cps.
    """
    log_e = np.log10(energy_kev)
    
    if orientation == 'parallel':
        coeffs = [
            -1232.80706, 5288.20962, -9584.54726, 9752.82964, -6193.51923,
            2559.84742, -691.24776, 117.93858, -11.56348, 0.49737
        ]
    elif orientation == 'perpendicular':
        coeffs = [
            -1603.32194, 6794.83543, -12203.37818, 12318.43163, -7758.24172,
            3176.31261, -848.01750, 142.69653, -13.75718, 0.57981
        ]
    else:
        raise ValueError("Hướng photon phải là 'parallel' hoặc 'perpendicular'")

    log_ge = np.zeros_like(log_e)
    for power, coeff in enumerate(coeffs):
        log_ge += coeff * (log_e ** power)
        
    return 10 ** log_ge

def compute_dose_rate(spectrum_cps: np.ndarray, energy_bins_kev: np.ndarray, orientation: str = 'parallel'):
    """
    Tính toán suất liều chiếu tổng X_dot (uR/h) từ phổ năng lượng thực nghiệm.
    """
    ge_values = calculate_ge_factor(energy_bins_kev, orientation)
    # Tích phân số dạng hình thang: Sum(n(E) * G(E) * dE)
    dose_rate = np.sum(spectrum_cps * ge_values)
    return dose_rate

Thông số thực nghiệm và so sánh phổ đo 24 giờ

Thông số đo kiểm Đầu dò bán dẫn HPGe (GC2518) Đầu dò nhấp nháy NaI(Tl) (7,6cm x 7,6cm)
Thời gian thực (Real time) $86.400\text{ s}$ (24 giờ) $86.774\text{ s}$
Thời gian sống (Live time) $86.088\text{ s}$ $86.400\text{ s}$
Thời gian chết (Dead time) $0,36%$ $0,48%$
Tổng số xung ghi nhận $15.205.000\text{ counts}$ $47.480.000\text{ counts}$
Tốc độ đếm trung bình $183\text{ cps}$ $549\text{ cps}$
Độ rộng kênh năng lượng $0,3\text{ keV/kênh}$ $6,7\text{ keV/kênh}$ 
Tốc độ đếm ghi nhận (cps):
HPGe GC2518       [█████████                           ] 183 cps
NaI(Tl) 7.6x7.6cm [███████████████████████████         ] 549 cps (Gấp 3.0 lần)

Kết quả định danh đồng vị phóng xạ môi trường

Phân tích phổ đối chứng HPGe và giải chập phổ $\text{NaI(Tl)}$ bằng Colegram đã định danh chính xác các đồng vị:

Năng lượng ($E_\gamma$, keV) Đồng vị nhận diện Chu kỳ bán rã ($T_{1/2}$) Chuỗi nguồn gốc Ghi nhận trên HPGe Ghi nhận trên NaI(Tl)
74,8 - 89,8 $^{212}\text{Pb}$ $10,64\text{ giờ}$ Thorium ($^{232}\text{Th}$) Có (tách rõ 4 đỉnh) Chồng chập nền X-ray
186,2 $^{226}\text{Ra}$ $1602\text{ năm}$ Uranium ($^{238}\text{U}$) Có Không tách được
238,6 $^{212}\text{Pb}$ $10,64\text{ giờ}$ Thorium ($^{232}\text{Th}$) Có Có (Colegram)
295,2 & 352,5 $^{214}\text{Pb}$ $26,8\text{ phút}$ Uranium ($^{238}\text{U}$) Có Có (Colegram)
583,1 $^{208}\text{Tl}$ $3,083\text{ phút}$ Thorium ($^{232}\text{Th}$) Có Có
609,3 $^{214}\text{Bi}$ $19,9\text{ giây}$ Uranium ($^{238}\text{U}$) Có Có
661,6 $^{137}\text{Cs}$ $30,17\text{ năm}$ Nguồn chuẩn nhân tạo Có Ẩn trong nền
834,8 $^{54}\text{Mn}$ $312,3\text{ ngày}$ Nguồn chuẩn nhân tạo Có Ẩn trong nền
911,1 & 969,1 $^{228}\text{Ac}$ $6,15\text{ giờ}$ Thorium ($^{232}\text{Th}$) Có Có (Colegram)
1460,8 $^{40}\text{K}$ $1,25 \times 10^9\text{ năm}$ NORM tự nhiên Có (đỉnh đơn cực mạnh) Có (đỉnh nổi bật)
1764,5 $^{214}\text{Bi}$ $19,9\text{ giây}$ Uranium ($^{238}\text{U}$) Có Có
2614,7 $^{208}\text{Tl}$ $3,083\text{ phút}$ Thorium ($^{232}\text{Th}$) Có (năng lượng cao) Có (đỉnh đơn rõ)

Kết quả tính toán suất liều chiếu tổng môi trường

Áp dụng hàm chuyển đổi $G(E)$ cho toàn bộ dải phổ đo được trên detector $\text{NaI(Tl)} \ 7,6\text{ cm} \times 7,6\text{ cm}$:

  • Cấu hình hướng photon song song trục detector: $$\dot{X}_{\parallel} = 12,43 \pm 0,53 \ \mu\text{R/h}$$
  • Cấu hình hướng photon vuông góc trục detector: $$\dot{X}_{\perp} = 10,66 \pm 0,45 \ \mu\text{R/h}$$
  • Kết luận mức phông: Suất liều chiếu phông phóng xạ tự nhiên tại phòng thí nghiệm dao động ổn định trong khoảng $10,66 - 12,43 \ \mu\text{R/h}$ (tương đương $0,093 - 0,109 \ \mu\text{Sv/h}$ hay $\approx 0,81 - 0,95 \text{ mSv/năm}$), hoàn toàn nằm trong giới hạn an toàn bức xạ theo TCVN và tiêu chuẩn ICRP đối với người dân ($< 1 \text{ mSv/năm}$ chiếu xạ ngoài gia tăng).

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

Đổi mới kỹ thuật nổi bật

  1. Tối ưu hóa mô hình toán học giải tích $G(E)$: Khóa luận xây dựng thành công hàm hồi quy đa thức bậc 9 với độ chính xác xấp xỉ gần như tuyệt đối ($R^2 = 0,99998$), loại bỏ sự phụ thuộc vào các bảng tra cứu rời rạc truyền thống.
  2. Kỹ thuật kết hợp đối chuẩn chéo phổ (Cross-Spectrometry Deconvolution): Sử dụng phổ HPGe làm "bản đồ chuẩn" (ground-truth mapping) để định hướng cho thuật toán phân rã peak trên phần mềm Colegram đối với phổ $\text{NaI(Tl)}$. Phương pháp này cho phép nhận dạng thành công 7 nhóm đồng vị NORM phức tạp từ phổ có độ phân giải thấp.
  3. Hiệu suất thu nhận thống kê vượt trội: Tinh thể $\text{NaI(Tl)}$ $7,6\text{ cm} \times 7,6\text{ cm}$ cho tốc độ đếm cao gấp 300% (549 cps so với 183 cps của HPGe), giúp giảm thời gian đo cần thiết để đạt cùng mức độ tin cậy thống kê (sai số $< 5%$) từ 24 giờ xuống còn dưới 2 giờ trong các ứng dụng thực địa.
So sánh hiệu năng tổng thể:
Hiệu suất ghi nhận (Counts):     NaI(Tl) [==============================] +212% vs HPGe
Tính cơ động hiện trường:         NaI(Tl) [==============================] Vượt trội (Không LN2)
Khả năng tách đỉnh đơn lẻ:       HPGe    [==============================] Vượt trội (<2 keV)
Tối ưu chi phí đầu tư (TCO):     NaI(Tl) [==============================] Tiết kiệm 75-80%

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

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

  • Trạm quan trắc cảnh báo sớm phóng xạ quốc gia: Triển khai các trạm giám sát tự động cố định dùng đầu dò $\text{NaI(Tl)}$ $3'' \times 3''$ kết hợp thuật toán $G(E)$ để truyền dữ liệu suất liều và phổ thời gian thực về trung tâm điều hành.
  • Kiểm soát phóng xạ cửa khẩu, cảng biển và hải quan: Lắp đặt trên hệ thống cổng kiểm xạ container (Radiation Portal Monitors) giúp phân biệt tức thời giữa bức xạ NORM từ nông sản, gốm sứ với các nguồn phóng xạ nguy hiểm hoặc vật liệu hạt nhân trái phép.
  • Khảo sát bức xạ môi trường dã ngoại và lập bản đồ phông: Tích hợp đầu dò $\text{NaI(Tl)}$ gọn nhẹ lên xe chuyên dụng hoặc thiết bị bay không người lái (UAV) để quét phổ diện rộng với chi phí tối ưu.
[ Môi trường thực địa / Cảng biển / Trạm NORM ]
                       |
                       v
[ Trạm đầu dò cơ động NaI(Tl) 7.6cm x 7.6cm ]
                       | (RS-485 / 4G LTE Telemetry)
                       v
[ Máy chủ quan trắc Trung tâm / Cloud Server ]

Phân tích chi phí - lợi ích (TCO & ROI)

Khoản mục chi phí Phương án truyền thống (HPGe Station) Phương án tối ưu (NaI(Tl) + Hàm G(E)) Tỷ lệ tiết kiệm
Chi phí đầu tư phần cứng ban đầu $65.000\text{ USD}$ $12.000\text{ USD}$ $81,5%$
Chi phí tiêu hao hàng năm ($\text{LN}_2$, bảo dưỡng) $4.500\text{ USD/năm}$ $300\text{ USD/năm}$ $93,3%$
Thời gian thiết lập và đo kiểm $4 - 6\text{ giờ}$ (chờ cân bằng nhiệt $\text{LN}_2$) Tức thời ($< 5\text{ phút}$) $95,0%$
Độ bền thiết bị trong môi trường khắc nghiệt Nhạy cảm với sốc nhiệt/cơ khí Rất bền bỉ, chống rung xóc Vượt trội

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

Hạn chế kỹ thuật

  • Độ phân giải năng lượng hữu hạn: Do cơ chế nhấp nháy phát quang, đầu dò $\text{NaI(Tl)}$ có độ phân giải FWHM $\approx 7%$ tại $661,6\text{ keV}$, khiến các đỉnh năng lượng lân cận hoặc các nguồn có hoạt độ rất thấp ($^{137}\text{Cs}$, $^{54}\text{Mn}$, $^{226}\text{Ra}$) dễ bị chìm trong nền tán xạ Compton.
  • Tính toán suất liều riêng phần (Individual Dose Rate): Nghiên cứu hiện tại giải quyết bài toán suất liều chiếu tổng; việc bóc tách chính xác tỷ lệ đóng góp liều của từng đồng vị riêng lẻ ($^{214}\text{Bi}$ đóng góp bao nhiêu $\mu\text{R/h}$, $^{40}\text{K}$ bao nhiêu $\mu\text{R/h}$) đòi hỏi giải thuật khử ma trận đáp ứng phức tạp hơn.

Hướng phát triển tiếp theo

  • Ứng dụng học máy (Machine Learning/Deep Learning): Xây dựng mạng nơ-ron tích chập (1D-CNN) hoặc Autoencoder để tự động phân rã phổ chồng chập trên $\text{NaI(Tl)}$ thời gian thực mà không cần đối chiếu thủ công qua phần mềm Colegram.
  • Mở rộng khảo sát dã ngoại In-situ: Triển khai các đợt đo thực địa ngoài trời tại các khu vực có dị thường phóng xạ tự nhiên (mỏ đất hiếm, mỏ than, khu vực đá hoa cương) nhằm hoàn thiện bộ tham số hàm $G(E)$ dưới điều kiện tán xạ môi trường bán không gian vô hạn ($2\pi$ geometry).

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

                                +-----------------------------------+
                                |    ĐỐI TƯỢNG HƯỞNG LỢI DỰ ÁN      |
                                +-----------------+-----------------+
                                                  |
            +--------------------+----------------+--------------------+--------------------+
            |                    |                                     |                    |
            v                    v                                     v                    v
+-----------------------+ +-----------------------+   +-----------------------+ +-----------------------+
|  Sinh viên / Học viên | |   Kỹ sư An toàn BX    |   |  Cơ quan Quản lý MT   | |  Nhà nghiên cứu HN    |
| - Tài liệu thực nghiệm| | - Quy trình đo chuẩn  |   | - Thiết lập trạm NORM | | - Dữ liệu khớp G(E)   |
| - Phương pháp Colegram| | - Giảm 80% TCO đầu tư |   | - Giám sát tự động 24h| | - Mở rộng phổ kế liều |
+-----------------------+ +-----------------------+   +-----------------------+ +-----------------------+
  • Sinh viên, học viên chuyên ngành Vật lý Hạt nhân & Kỹ thuật Hạt nhân: Cung cấp tài liệu thực chứng toàn diện về kỹ thuật đo phổ gamma, phương pháp hiệu chuẩn năng lượng và quy trình xử lý phổ thực tế.
  • Kỹ sư an toàn bức xạ và cán bộ quan trắc: Nắm bắt giải pháp kỹ thuật tối ưu chi phí để chế tạo hoặc cấu hình thiết bị đo liều phổ kế mà không cần đầu tư các hệ thống bán dẫn đắt tiền.
  • Cơ quan quản lý nhà nước và Môi trường: Có cơ sở khoa học tin cậy để triển khai các mạng lưới quan trắc phông phóng xạ tự nhiên diện rộng, đảm bảo an ninh nguồn phóng xạ và an toàn sức khỏe cộng đồng.

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

1. Yêu cầu phần cứng và phần mềm tối thiểu để triển khai hệ đo $\text{NaI(Tl)} + G(E)$ là gì?

Hệ thống yêu cầu đầu dò nhấp nháy $\text{NaI(Tl)}$ chuẩn $7,6\text{ cm} \times 7,6\text{ cm}$ kết nối bộ tiền khuếch đại, hệ thu nhận MCA tối thiểu 1024 kênh (khuyến nghị Canberra InSpector 1000 hoặc tương đương). Về phần mềm, cần công cụ phân tích đỉnh phổ (như Colegram) và module tính toán giải tích hỗ trợ tích phân số (Python/MATLAB/C++).

2. Sự khác biệt về hướng photon tới (song song vs vuông góc) ảnh hưởng thế nào đến suất liều?

Hình dạng hình học trụ của tinh thể $\text{NaI(Tl)}$ dẫn đến chiều dày hấp thụ hiệu dụng khác nhau theo các góc tới. Khi photon tới vuông góc, diện tích tiết diện hình học lớn hơn nhưng bề dày tương tác trung bình khác so với hướng song song dọc trục. Kết quả thực nghiệm cho thấy suất liều theo hướng song song đạt $12,43 \pm 0,53 \ \mu\text{R/h}$, trong khi hướng vuông góc đạt $10,66 \pm 0,45 \ \mu\text{R/h}$ (chênh lệch khoảng $14%$). Trong thực địa bức xạ môi trường phân bố đẳng hướng ($4\pi$), giá trị thực tế sẽ là giá trị trung bình góc của hai hàm đáp ứng này.

3. Tại sao đầu dò $\text{NaI(Tl)}$ không nhận diện được trực tiếp $^{137}\text{Cs}$ ($661,6\text{ keV}$) và $^{54}\text{Mn}$ ($834,8\text{ keV}$) trong khi HPGe lại thấy rõ?

Do độ phân giải năng lượng của $\text{NaI(Tl)}$ ở mức trung bình (~$7%$), đỉnh $661,6\text{ keV}$ của $^{137}\text{Cs}$ bị chìm và chồng lấn hoàn toàn với đỉnh năng lượng rất mạnh $609,3\text{ keV}$ của $^{214}\text{Bi}$ thuộc chuỗi Uranium tự nhiên. Tương tự, đỉnh $834,8\text{ keV}$ của $^{54}\text{Mn}$ bị che khuất bởi các đỉnh $794,9\text{ keV}$ và $911,1\text{ keV}$ của $^{228}\text{Ac}$. Đầu dò HPGe với FWHM cực hẹp ($< 2\text{ keV}$) tách biệt hoàn toàn các đỉnh này.

4. Hệ số Kerma và liều hấp thụ, liều chiếu khác nhau như thế nào trong bài toán này?

  • Kerma ($K$): Tổng động năng ban đầu của tất cả các hạt mang điện thứ cấp sinh ra trên một đơn vị khối lượng vật chất bởi bức xạ gián tiếp (photon), đơn vị là Gray ($\text{Gy} = \text{J/kg}$).
  • Liều hấp thụ ($D$): Phần năng lượng bức xạ thực tế truyền lại và hấp thụ trong môi trường thông qua ion hóa và kích thích nguyên tử ($\text{Gy}$).
  • Liều chiếu ($X$): Thước đo khả năng ion hóa không khí khô của chùm tia X hoặc gamma, tính bằng điện tích sinh ra trên một đơn vị khối lượng không khí ($\text{C/kg}$ hoặc Roentgen, $1\text{ R} = 2,58 \times 10^{-4}\text{ C/kg}$).

5. Phương pháp này có thể áp dụng để đo phóng xạ liên tục ngoài trời (outdoor) không?

Hoàn toàn khả thi. Tuy nhiên, khi chuyển sang môi trường ngoài trời, cần bổ sung: (1) Hộp bảo vệ đầu dò chống ẩm và ổn nhiệt (nhiệt độ làm trôi đỉnh quang điện của NaI); (2) Hiệu chuẩn lại hàm $G(E)$ với điều kiện nguồn bức xạ phân bố dạng nửa không gian ($2\pi$ soil geometry) thay vì hình học phòng thí nghiệm khép kín.

Kết luận

Khóa luận tốt nghiệp đã giải quyết thành công bài toán xác định suất liều phóng xạ môi trường bằng phương pháp phổ kế gamma trên đầu dò nhấp nháy $\text{NaI(Tl)} \ 7,6\text{ cm} \times 7,6\text{ cm}$, thông qua các đóng góp then chốt:

  1. Làm chủ kỹ thuật phân tích phổ: Kết hợp hiệu quả giữa độ phân giải cao của hệ chuẩn HPGe và phần mềm Colegram để bóc tách, định danh chính xác các chuỗi phóng xạ tự nhiên $^{238}\text{U}$, $^{232}\text{Th}$ và $^{40}\text{K}$ trên phổ $\text{NaI(Tl)}$.
  2. Chuẩn hóa hàm chuyển đổi $G(E)$: Xây dựng bộ tham số đa thức bậc 9 hoàn chỉnh với độ tin cậy vượt bậc ($R^2 > 0,9999$), thiết lập nền tảng tính toán suất liều chiếu tức thời từ phổ đếm thô.
  3. Đánh giá định lượng an toàn môi trường: Xác định chính xác suất liều phông tự nhiên tại phòng thí nghiệm đạt $10,66 - 12,43 \ \mu\text{R/h}$, khẳng định tính an toàn bức xạ của cơ sở nghiên cứu.

Nghiên cứu mở ra hướng đi thực tiễn cho việc phát triển các thiết bị đo bức xạ cơ động, trạm quan trắc cảnh báo sớm thông minh với chi phí thấp, đóng góp thiết thực cho công tác kiểm soát an toàn bức xạ và an ninh hạt nhân tại Việt Nam.