Giới thiệu dự án

Bức xạ ion hóa môi trường là thành phần hiện diện thường trực trong sinh quyển Trái Đất, bắt nguồn từ các nguồn phóng xạ tự nhiên (tia vũ trụ, các chuỗi phân rã phóng xạ nguyên thủy trong vỏ Trái Đất) và các nguồn nhân tạo (cơ sở hạt nhân, y tế bức xạ, ứng dụng chiếu xạ công nghiệp). Theo các báo cáo từ UNSCEAR (Ủy ban Khoa học của Liên Hợp Quốc về Ảnh hưởng của Bức xạ Nguyên tử), 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ừ bức xạ tự nhiên đạt khoảng $2.4\text{ mSv/năm}$, trong đó chiếu xạ ngoài từ đất đá đóng góp khoảng $0.45\text{ mSv/năm}$, nước biển và thực phẩm đóng góp $0.25\text{ mSv/năm}$, phần còn lại chủ yếu đến từ khí Radon ($^{222}\text{Rn}$) tích tụ trong không gian kín (khi nồng độ đạt mức $150\text{ Bq/m}^3$, nguy cơ tử vong do ung thư phổi tăng từ 1% đến 3%).

Thách thức cốt lõi trong quan trắc an toàn bức xạ hiện đại là nhu cầu giám sát liên tục, nhanh chóng và chính xác suất liều bức xạ môi trường. Các đầu dò bán dẫn siêu tinh khiết High-Purity Germanium (HPGe) sở hữu độ phân giải năng lượng vượt trội nhưng cồng kềnh, chi phí đắt đỏ và bắt buộc phải làm lạnh liên tục bằng nitơ lỏng ($77\text{ K}$). Ngược lại, detector nhấp nháy Thallium-doped Sodium Iodide ($\text{NaI(Tl)}$) có hiệu suất ghi nhận cao, hoạt động ở nhiệt độ phòng và tính linh động tuyệt vời nhưng lại có độ phân giải năng lượng thấp, dẫn đến hiện tượng chồng chập đỉnh phổ gamma phức tạp.

Đề tài "Xác định suất liều phóng xạ môi trường sử dụng detector $\text{NaI(Tl)}\ 7,6\text{cm} \times 7,6\text{cm}$" (thực hiện bởi sinh viên Ngô Vũ Thiên Quang dưới sự hướng dẫn của TS. Võ Hồng Hải tại Bộ môn Vật lý Hạt nhân, Trường Đại học Sư phạm TP.HCM phối hợp cùng 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) giải quyết trực tiếp bài toán này.

[Trường bức xạ Gamma Môi trường]

Mục tiêu cụ thể của dự án:

  1. Ghi nhận và phân tích phổ năng lượng gamma phông môi trường liên tục trong 24 giờ ($86.400\text{ s}$) bằng cả hai hệ đo: detector bán dẫn HPGe Canberra GC2518 và detector nhấp nháy $\text{NaI(Tl)}\ 7,6\text{cm} \times 7,6\text{cm}$ (InSpector™ 1000).
  2. Ứng dụng phần mềm chuyên dụng Colegram kết hợp cơ sở dữ liệu Laraweb và Bảng đồng vị hạt nhân để phân tách các vùng phổ chồng chập, định danh các hạt nhân phóng xạ thuộc chuỗi Uranium ($^{238}\text{U}$), Thorium ($^{232}\text{Th}$), Actinium ($^{235}\text{U}$) và $^{40}\text{K}$.
  3. Thiết lập và tối ưu hóa hàm chuyển đổi suất liều $G(E)$ theo hai hướng photon tới (song song và vuông góc với trục detector) thông qua mô hình đa thức bậc 9 với hệ số tương quan $R^2 > 0.9999$.
  4. Xác định chính xác suất liều chiếu tổng ($\mu\text{R/h}$) của môi trường phòng thí nghiệm từ phổ năng lượng thực nghiệm $\text{NaI(Tl)}$.

Phạm vi nghiên cứu tập trung vào dải năng lượng gamma tự nhiên từ $50\text{ keV}$ đến $3000\text{ keV}$ tại môi trường phòng thí nghiệm kỹ thuật hạt nhân, thiết lập cơ sở phương pháp luận vững chắc cho các trạm quan trắc bức xạ tự động ngoài hiện trường.

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

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

Việc lựa chọn thiết bị đo bức xạ môi trường luôn đòi hỏi sự cân bằng giữa độ chính xác phổ, hiệu suất ghi và tính khả thi trong triển khai:

Tiêu chí kỹ thuật Detector bán dẫn HPGe Detector nhấp nháy $\text{NaI(Tl)}$ Đầu dò ống đếm Geiger-Müller (GM)
Độ phân giải năng lượng (FWHM) Cực cao ($< 0.2%$ tại $1332\text{ keV}$) Trung bình/Thấp ($6.5% - 7.5%$) Không có khả năng phân giải năng lượng
Hiệu suất ghi nhận photon Trung bình Rất cao (thể tích tinh thể $3" \times 3"$) Rất thấp đối với tia gamma ($< 1% - 2%$)
Khả năng định danh đồng vị Xuất sắc, tách biệt hoàn toàn các đỉnh Cần thuật toán phân rã phổ (deconvolution) Không thể thực hiện
Điều kiện vận hành Nitơ lỏng ($LN_2$), bảo trì phức tạp Nhiệt độ phòng, cơ động cao Đơn giản, nguồn pin gọn nhẹ
Xác định suất liều theo năng lượng Phức tạp, hiệu chỉnh hình học lớn Trực tiếp qua hàm chuyển đổi $G(E)$ Sai số lớn khi năng lượng thay đổi

Yêu cầu kỹ thuật theo mô hình MoSCoW:

  • Must-have: Khả năng ghi nhận phổ liên tục không gián đoạn trong 24 giờ; tích hợp thuật toán hàm chuyển đổi $G(E)$ để tính suất liều tổng; tách được các đỉnh gamma chủ đạo của $^{214}\text{Bi}$, $^{208}\text{Tl}$, $^{40}\text{K}$, $^{212}\text{Pb}$.
  • Should-have: Mô hình hóa toán học hàm $G(E)$ cho các góc tới khác nhau (song song và vuông góc); đối chứng song song với hệ phổ kế HPGe độ chính xác cao.
  • Could-have: Tách phổ tự động các vùng đa đỉnh phức tạp (multiplets) bằng hàm phân bố Gaussian kết hợp đa thức nền trên Colegram.
  • Won't-have: Đo đạc phổ neutron hoặc phân tích bức xạ alpha/beta trực tiếp trong cùng chu kỳ đo.

Thiết kế hệ thống

Kiến trúc hệ thống đo lường và xử lý số liệu bức xạ bao gồm 2 khối chính: Khối thu nhận phần cứng (Hardware Acquisition Layer) và Khối phân tích số liệu phổ hạt nhân (Nuclear Spectroscopy Processing Layer).

|                         HỆ THỐNG PHẦN CỨNG THỰC NGHIỆM                      |
|                                                                             |
|                     KHỐI XỬ LÝ & TÍNH TOÁN LIỀU BỨC XẠ                      |
|                                                                             |
|  1. Chuẩn hóa năng lượng & Hiệu chỉnh thời gian chết (Dead time < 0.5%)      |
|  2. Phần mềm Colegram: Khớp đỉnh Gaussian + Nền phi tuyến                   |
|  3. Laraweb Database: Định danh đồng vị chuỗi U-238, Th-232, K-40           |
|  4. Tích phân phổ kết hợp hàm chuyển đổi liều G(E) theo mô hình đa thức     |

Công nghệ và công cụ sử dụng:

  • Hệ đo bán dẫn HPGe: Đầu dò đồng trục Canberra model GC2518, đường kính tinh thể $61.8\text{ mm}$, bề dày $35.8\text{ mm}$, tích hợp bình chứa nitơ lỏng và hệ giá đỡ cơ động; kết nối bộ phân tích đa kênh MCA 8k kênh ($8192\text{ channels}$).
  • Hệ đo nhấp nháy $\text{NaI(Tl)}$: Canberra InSpector™ 1000 với đầu dò tinh thể hình trụ đồng trục kích thước $7,6\text{cm} \times 7,6\text{cm}$ ($3" \times 3"$), tích hợp mạch tiền khuếch đại, bộ xử lý tín hiệu số DSP và hiển thị suất liều trực tiếp.
  • Phần mềm xử lý phổ: Colegram (Nuclear Data Processing Software), bộ thư viện dữ liệu phân rã hạt nhân trực tuyến Laraweb (CEA/LNHB) và Bảng đồng vị phóng xạ (Table of Isotopes 8th Edition).

Methodology

Phương pháp nghiên cứu kết hợp giữa đo đạc thực nghiệm vật lý hạt nhân và mô phỏng giải tích hàm đáp ứng detector:

  1. Thiết lập hình học đo: Đặt hai đầu dò HPGe và $\text{NaI(Tl)}$ tại cùng vị trí trung tâm trong Phòng thí nghiệm Kỹ thuật Hạt nhân (Đại học Khoa học Tự nhiên TP.HCM).
  2. Thu nhận dữ liệu 24 giờ liên tục: Đo phông môi trường trong chu kỳ $86.400\text{ s}$ nhằm đạt độ thống kê số đếm cao, giảm sai số thống kê Poisson ($\sigma = \sqrt{N}$) xuống mức tối thiểu.
  3. Phân rã phổ và khớp hàm: Sử dụng Colegram để thực hiện khớp phi tuyến Levenberg-Marquardt trên từng cửa sổ năng lượng hẹp ($190 - 430\text{ keV}$, $430 - 550\text{ keV}$, ..., $2400 - 2900\text{ keV}$).
  4. Tính toán suất liều thông qua tích chập ma trận đáp ứng: Áp dụng hệ số chuyển đổi phổ $G(E)$ thu được từ nghịch đảo ma trận đáp ứng của đầu dò kết hợp hệ số hấp thụ năng lượng khối của không khí.

Implementation và kết quả

Development process & Cơ sở toán học

Mối liên hệ giữa trường thông lượng photon $\psi(E)$ và suất liều chiếu $\dot{X}$ được xây dựng từ lý thuyết Kerma (Kinetic Energy Released per unit Mass). Khi photon tương tác với vật chất, năng lượng truyền cho các hạt mang điện thứ cấp tạo nên Kerma va chạm ($K_c$). Trong không khí khô, suất liều chiếu $\dot{X}$ liên hệ với Kerma va chạm qua biểu thức:

$$\dot{X} = K_c \cdot \frac{e}{W_{\text{air}}} = \int_{0}^{E_{\max}} \psi(E) \cdot \left(\frac{\mu_{\text{en}}}{\rho}\right){\text{air}} \cdot \frac{e}{W{\text{air}}} , dE$$

Trong đó $W_{\text{air}}/e = 33.97\text{ J/C}$ là năng lượng trung bình tạo một cặp ion trong không khí khô.

Đối với đầu dò $\text{NaI(Tl)}$, tốc độ đếm $n(E)$ tại kênh năng lượng đo được liên hệ với thông lượng photon tới thông qua hàm đáp ứng $R(E, E_0)$. Bằng phương pháp ma trận đáp ứng nghịch đảo, suất liều chiếu tổng được tính trực tiếp từ phổ thực nghiệm $n(E)$ không cần giải phổ ngược phức tạp:

$$\dot{X} = \int_{E_{\min}}^{E_{\max}} n(E) \cdot G(E) , dE = \sum_{i} n(E_i) \cdot G(E_i) \cdot \Delta E_i$$

Hệ số chuyển đổi $G(E)$ ($\mu\text{R}\cdot\text{h}^{-1}/\text{cps}$) được làm khớp chính xác bằng hàm đa thức bậc 9 theo $\log E$ trên hai cấu hình hình học:

# Mô hình giải tích tính toán hệ số chuyển đổi suất liều G(E)
import numpy as np

def calculate_g_factor(energy_kev, orientation='parallel'):
    """
    Tính hệ số chuyển đổi G(E) [uR/h/cps] cho detector NaI(Tl) 7.6cm x 7.6cm
    Energy dải từ 40 keV đến 3000 keV.
    """
    log_e = np.log10(energy_kev)
    
    if orientation == 'parallel':
        # Bức xạ tới song song trục đối xứng detector (R^2 = 0.99998)
        coeffs = [
            -1232.80706,  5288.20962, -9584.54726,  9752.82964,
            -6193.51923,  2559.84742,  -691.24776,   117.93858,
              -11.56348,     0.49737
        ]
    elif orientation == 'perpendicular':
        # Bức xạ tới vuông góc trục đối xứng detector (R^2 = 0.99996)
        coeffs = [
            -1603.43163,  6788.24172, -12176.31261, 12217.50175,
            -7758.24172,  3176.31261,  -849.01750,   142.69653,
              -13.75718,     0.57981
        ]
    else:
        raise ValueError("Hướng đo không hợp lệ.")
        
    log_g = sum(c * (log_e ** i) for i, c in enumerate(coeffs))
    return 10 ** log_g

Testing và validation

Quá trình đo đạc thực nghiệm liên tục trong 24 giờ ($1\text{ ngày đêm}$) mang lại bộ thông số vận hành ổn định:

Thông số đo thực nghiệm Detector HPGe Canberra GC2518 Detector $\text{NaI(Tl)}\ 7,6\text{cm} \times 7,6\text{cm}$
Thời gian đo thực tế (Real time) $86.400\text{ s}$ $86.774\text{ s}$
Thời gian sống của hệ (Live time) $86.088\text{ s}$ $86.358\text{ s}$
Thời gian chết của hệ (Dead time) $0.36%$ $0.48%$
Tổng số đếm tích lũy (Counts/day) $15.728.205\text{ counts}$ $47.480.549\text{ counts}$
Tốc độ đếm trung bình (Count rate) $183\text{ cps}$ $549\text{ cps}$
Độ rộng kênh MCA $0.36\text{ keV/kênh}$ ($8192\text{ kênh}$) $6.7\text{ keV/kênh}$ ($1024\text{ kênh}$)

Tốc độ đếm của đầu dò $\text{NaI(Tl)}$ cao gấp 3.0 lần so với HPGe ($549\text{ cps}$ so với $183\text{ cps}$), chứng minh ưu thế vượt trội về hiệu suất thu nhận thể tích tinh thể lớn trong trường bức xạ môi trường có hoạt độ thấp.

Kết quả đạt được

Hệ phổ kế HPGe phân giải được hơn 50 đỉnh gamma đơn lẻ, định danh chính xác sự có mặt của các chuỗi phóng xạ tự nhiên và nguồn chuẩn nhân tạo tồn tại trong phòng thí nghiệm. Phổ $\text{NaI(Tl)}$ sau khi qua xử lý Colegram đã bóc tách thành công các đỉnh năng lượng then chốt:

Năng lượng ($E_\gamma$) Đồng vị định danh 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 $\text{NaI(Tl)}$
$74.8 - 87.2\text{ keV}$ $^{212}\text{Pb}$ $10.64\text{ giờ}$ Thorium ($^{232}\text{Th}$) Có (tách 4 đỉnh) Chập thành dải X-ray
$186.2\text{ keV}$ $^{226}\text{Ra}$ / $^{235}\text{U}$ $1602\text{ năm}$ Uranium ($^{238}\text{U}$) Có Không phân tách
$238.6\text{ keV}$ $^{212}\text{Pb}$ $10.64\text{ giờ}$ Thorium ($^{232}\text{Th}$) Có Có
$295.2\text{ keV}$ $^{214}\text{Pb}$ $26.8\text{ phút}$ Uranium ($^{238}\text{U}$) Có Có
$352.5\text{ keV}$ $^{214}\text{Pb}$ $26.8\text{ phút}$ Uranium ($^{238}\text{U}$) Có Có
$583.1\text{ keV}$ $^{208}\text{Tl}$ $3.083\text{ phút}$ Thorium ($^{232}\text{Th}$) Có Có
$609.3\text{ keV}$ $^{214}\text{Bi}$ $19.9\text{ phút}$ Uranium ($^{238}\text{U}$) Có Có
$661.6\text{ keV}$ $^{137}\text{Cs}$ $30.17\text{ năm}$ Nhân tạo (Nguồn chuẩn Lab) Có Bị chập đỉnh
$727.2\text{ keV}$ $^{212}\text{Bi}$ $60.55\text{ phút}$ Thorium ($^{232}\text{Th}$) Có Có
$834.8\text{ keV}$ $^{54}\text{Mn}$ $312.3\text{ ngày}$ Nhân tạo (Nguồn chuẩn Lab) Có Bị chập đỉnh
$911.1 - 969.1\text{ keV}$ $^{228}\text{Ac}$ $6.15\text{ giờ}$ Thorium ($^{232}\text{Th}$) Có Có
$1120.3\text{ keV}$ $^{214}\text{Bi}$ $19.9\text{ phút}$ Uranium ($^{238}\text{U}$) Có Có
$1460.8\text{ keV}$ $^{40}\text{K}$ $1.25 \times 10^9\text{ năm}$ Phóng xạ tự nhiên đơn lẻ Có Có
$1764.5\text{ keV}$ $^{214}\text{Bi}$ $19.9\text{ phút}$ Uranium ($^{238}\text{U}$) Có Có
$2614.7\text{ keV}$ $^{208}\text{Tl}$ $3.083\text{ phút}$ Thorium ($^{232}\text{Th}$) Có Có

Tích phân phổ năng lượng $\text{NaI(Tl)}$ với hàm $G(E)$ theo hai hướng tới cho kết quả suất liều chiếu môi trường tại phòng thí nghiệm:

  • Hướng photon song song trục detector: $\dot{X}_{\parallel} = 12.43 \pm 0.53\ \mu\text{R/h}$ ($0.109 \pm 0.005\ \mu\text{Sv/h}$).
  • Hướng photon vuông góc trục detector: $\dot{X}_{\perp} = 10.66 \pm 0.45\ \mu\text{R/h}$ ($0.093 \pm 0.004\ \mu\text{Sv/h}$).
  • Suất liều thực tế trong không gian đẳng hướng phòng thí nghiệm dao động trong dải tin cậy $10.66 - 12.43\ \mu\text{R/h}$, hoàn toàn phù hợp với giới hạn an toàn phông tự nhiên theo TCVN ($< 20\ \mu\text{R/h}$).

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

  1. Chuẩn hóa quy trình chuyển đổi phổ sang suất liều bằng hàm giải tích $G(E)$: Khóa luận áp dụng thành công mô hình ma trận đáp ứng nghịch đảo chuyển đổi trực tiếp toàn bộ phổ biên độ đa kênh thành suất liều vật lý thực tế mà không cần giả định trước phân bố trường bức xạ.
  2. Khớp giải tích đa thức bậc cao với độ chính xác tuyệt đối: Xây dựng 2 hàm giải tích bậc 9 cho hàm $G(E)$ trên dải năng lượng rộng $40 - 3000\text{ keV}$ đạt hệ số tương quan $R^2 = 0.99998$ (song song) và $R^2 = 0.99996$ (vuông góc), loại bỏ sai số nội suy dạng bảng thông thường.
  3. Quy trình kết hợp Colegram giải mã phổ nhấp nháy $\text{NaI(Tl)}$: Đưa ra phương pháp luận chuẩn sử dụng dữ liệu đối chứng HPGe làm chuẩn vị trí đỉnh, hỗ trợ thuật toán phân tách phổ Colegram bóc tách chính xác các cụm đỉnh chập phổ biến của $^{214}\text{Pb}$/ $^{212}\text{Pb}$ và $^{228}\text{Ac}$/ $^{214}\text{Bi}$.
  4. Chứng minh tính khả thi thay thế HPGe trong giám sát an toàn bức xạ công nghiệp: Giảm thiểu chi phí đầu tư thiết bị hơn 70% và loại bỏ 100% chi phí vận hành bảo trì nitơ lỏng định kỳ.

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

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

  • Mạng lưới trạm quan trắc cảnh báo phóng xạ quốc gia: Triển khai các trạm đo bức xạ tự động liên tục tại các khu vực biên giới, vùng lân cận nhà máy điện hạt nhân hoặc cơ sở lò phản ứng nghiên cứu.
  • Giám sát phế liệu kim loại và vật liệu xây dựng: Ứng dụng đầu dò $\text{NaI(Tl)}$ tích hợp hàm $G(E)$ tại các cổng kiểm soát phóng xạ cảng biển, nhà máy xi măng, chế biến đá granite để phát hiện dị thường phóng xạ ($^{226}\text{Ra}$, $^{232}\text{Th}$, $^{137}\text{Cs}$).
  • Ứng phó sự cố bức xạ và chống khủng bố hạt nhân: Hệ thống thiết bị InSpector™ 1000 gọn nhẹ cho phép lực lượng an ninh/cứu hộ triển khai cơ động, đo quét nhanh suất liều và nhận diện nguồn phát tán trong vòng vài giây.
[Môi trường / Hiện trường]

Phân tích hiệu quả kinh tế & Lộ trình triển khai

  • Chi phí thiết bị: Đầu tư hệ HPGe hoàn chỉnh dao động từ $$80.000 - $120.000$ kèm chi phí nitơ lỏng định kỳ khoảng $$3.000 - $5.000/\text{năm}$. Hệ đo $\text{NaI(Tl)}\ 3" \times 3"$ hoàn chỉnh có chi phí chỉ từ $$15.000 - $25.000$, không tiêu hao phụ phẩm, tuổi thọ tinh thể $> 10\text{ năm}$.
  • Lộ trình triển khai:
    • Giai đoạn 1 (1 - 3 tháng): Hiệu chuẩn phổ năng lượng với các nguồn chuẩn đơn năng ($^{137}\text{Cs}$, $^{60}\text{Co}$, $^{241}\text{Am}$).
    • Giai đoạn 2 (4 - 6 tháng): Nhúng thuật toán tính toán suất liều $G(E)$ và tách phổ vào bộ vi điều khiển (MCU/DSP) trên bo mạch thiết bị.
    • Giai đoạn 3 (7 - 12 tháng): Triển khai thử nghiệm 5 trạm quan trắc ngoài trời đạt chuẩn IP67, truyền dữ liệu IoT thời gian thực về máy chủ giám sát.

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

Hạn chế kỹ thuật

  • Độ phân giải năng lượng giới hạn: Không thể bóc tách các đỉnh năng lượng quá gần nhau có độ chênh $< 50\text{ keV}$ ở vùng năng lượng thấp (ví dụ: đỉnh $^{235}\text{U}\ 185.7\text{ keV}$ và $^{226}\text{Ra}\ 186.2\text{ keV}$; các đỉnh thoát đơn, thoát kép và tán xạ ngược).
  • Hiện tượng trôi đỉnh do nhiệt độ: Tinh thể $\text{NaI(Tl)}$ và ống nhân quang điện (PMT) nhạy cảm với biến thiên nhiệt độ môi trường, đòi hỏi phải có mạch ổn định đỉnh số (Digital Spectrum Stabilizer) khi đặt ngoài trời.
  • Nguồn nhân tạo nồng độ vết: Các đỉnh nhân tạo cường độ thấp như $^{137}\text{Cs}$ ($661.6\text{ keV}$) và $^{54}\text{Mn}$ ($834.8\text{ keV}$) bị chìm vào nền phổ liên tục của $\text{NaI(Tl)}$, chỉ phát hiện được khi dùng HPGe.

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

  • Xác định suất liều riêng phần (Individual Dose Rate Spectrometry): Ứng dụng kỹ thuật $G(E)$ tách biệt để định lượng chính xác tỷ lệ đóng góp suất liều của từng đồng vị riêng lẻ ($^{40}\text{K}$, $^{232}\text{Th}$, $^{238}\text{U}$) thay vì chỉ tính suất liều tổng.
  • Nghiên cứu mở rộng trường bức xạ ngoài trời (In-situ Environmental Monitoring): Tiến hành các đợt đo dã ngoại tại các vùng có phông bức xạ cao (mỏ sa khoáng titan, vùng đất đỏ bazan Tây Nguyên) để hoàn thiện hàm đáp ứng $G(E)$ trong điều kiện tán xạ tự do.

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

  • Sinh viên & Học viên chuyên ngành Vật lý/Kỹ thuật Hạt nhân: Cung cấp tài liệu mẫu mực về phân tích phổ gamma, kỹ thuật phân rã phổ bằng Colegram và lý thuyết định lượng liều bức xạ hiện đại.
  • Kỹ sư An toàn Bức xạ & Cơ quan Quản lý Môi trường: Sở hữu phương pháp luận được chuẩn hóa để xây dựng các trạm giám sát bức xạ trực tuyến, giá thành thấp, độ tin cậy cao.
  • Doanh nghiệp Xử lý Chất thải & Công nghiệp Cảng biển: Nắm bắt công nghệ kiểm soát liều tự động và phát hiện nguồn phóng xạ mồ côi nhanh chóng, tối ưu hóa chi phí tuân thủ quy chuẩn an toàn bức xạ hạt nhân.
  • Cộng đồng Nghiên cứu Khoa học: Bộ dữ liệu thực nghiệm so sánh đồng thời giữa hai hệ đo HPGe và $\text{NaI(Tl)}$ trong cùng điều kiện hình học 24 giờ, làm cơ sở kiểm chứng các mô hình mô phỏng Monte Carlo (MCNP, Geant4).

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

1. Yêu cầu kỹ thuật phần cứng tối thiểu để triển khai phương pháp đo này là gì?

Cần một đầu dò nhấp nháy $\text{NaI(Tl)}$ kích thước $3" \times 3"$ (ví dụ Canberra Model 802 hoặc tương đương), nguồn nuôi cao thế phân cực từ $+800\text{V}$ đến $+1200\text{V}$, mạch tiền khuếch đại tích hợp và một bộ phân tích đa kênh (MCA) tối thiểu 1024 kênh kết nối máy tính qua chuẩn USB/RS-232.

2. Giới hạn dải năng lượng đo lường và tính toán suất liều là bao nhiêu?

Phương pháp hàm $G(E)$ cho đầu dò $\text{NaI(Tl)}\ 7,6\text{cm} \times 7,6\text{cm}$ hoạt động tối ưu trong dải năng lượng photon gamma từ $40\text{ keV}$ đến $3000\text{ keV}$. Dưới $40\text{ keV}$, hiệu ứng hấp thụ vỏ bọc đầu dò làm sai lệch hàm đáp ứng; trên $3000\text{ keV}$, xác suất tương tác tạo cặp chiếm ưu thế và thông lượng gamma tự nhiên gần như triệt tiêu.

3. Làm thế nào để tích hợp thuật toán này vào hệ thống giám sát SCADA hiện hữu?

Số liệu phổ từ MCA được truy xuất định kỳ qua tập lệnh API/Modbus. Chương trình điều khiển trung tâm (C/C++ hoặc Python) thực hiện phép nhân ma trận giữa vector tốc độ đếm kênh $n_i$ với vector hệ số $G_i$ đã được lập trình sẵn để xuất kết quả suất liều trực tiếp theo thời gian thực mỗi chu kỳ $60\text{ s} - 300\text{ s}$.

4. Thiết bị cần bảo dưỡng và hiệu chuẩn định kỳ như thế nào?

Cần kiểm tra chuẩn hóa năng lượng định kỳ mỗi 3 đến 6 tháng bằng nguồn chuẩn kép gắn kín hoạt độ thấp ($^{137}\text{Cs}\ 661.6\text{ keV}$ và $^{60}\text{Co}\ 1173.2/1332.5\text{ keV}$). Do tinh thể $\text{NaI(Tl)}$ bị hút ẩm, đầu dò phải luôn được duy trì bọc kín nguyên khối; không yêu cầu bổ sung chất làm lạnh.

5. Chi phí đầu tư và thời gian hoàn vốn (ROI) ước tính ra sao?

Chi phí ban đầu khoảng $$15.000$ (tiết kiệm hơn $75%$ so với hệ thống bán dẫn HPGe). Thời gian hoàn vốn thể hiện qua việc giảm thiểu rủi ro vi phạm an toàn bức xạ, cắt giảm $100%$ chi phí mua nitơ lỏng định kỳ hàng năm và tiết kiệm nhân lực trực tiếp đo đạc tại hiện trường nhờ cơ chế tự động hóa hoàn toàn.

Kết luận

Khóa luận tốt nghiệp của tác giả Ngô Vũ Thiên Quang đã giải quyết trọn vẹn và chuẩn xác bài toán xác định suất liều bức xạ gamma môi trường thông qua việc kết hợp đầu dò nhấp nháy $\text{NaI(Tl)}\ 7,6\text{cm} \times 7,6\text{cm}$ và kỹ thuật hàm chuyển đổi $G(E)$. Nghiên cứu không chỉ định danh thành công các đồng vị phóng xạ tự nhiên chủ đạo ($^{228}\text{Ac}$, $^{212}\text{Bi}$, $^{214}\text{Bi}$, $^{208}\text{Tl}$, $^{40}\text{K}$, $^{212}\text{Pb}$, $^{214}\text{Pb}$) với sự đối chứng tin cậy từ phổ kế bán dẫn HPGe Canberra GC2518, mà còn đưa ra giá trị suất liều chiếu tổng chuẩn xác tại phòng thí nghiệm ($10.66 - 12.43\ \mu\text{R/h}$).

Thành công của đề tài khẳng định khả năng làm chủ các công cụ xử lý phổ hạt nhân tiên tiến (Colegram, Laraweb) và mở ra tiềm năng ứng dụng to lớn trong việc xây dựng các trạm quan trắc bức xạ tự động hóa diện rộng, đáp ứng yêu cầu cấp thiết về an toàn hạt nhân và bảo vệ môi trường trong kỷ nguyên phát triển công nghệ cao.