Giới thiệu dự án

Quan trắc và phân tích phóng xạ môi trường là một trong những lĩnh vực cốt lõi của vật lý hạt nhân ứng dụng, đóng vai trò sống còn trong đánh giá an toàn bức xạ, khảo sát địa chất khoáng sản và giám sát an ninh môi trường sinh thái. Theo các báo cáo từ Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA) và Ủ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ử (UNSCEAR), hơn 80% liều bức xạ ion hóa tự nhiên mà con người tiếp nhận hàng năm có nguồn gốc từ các nhân phóng xạ nguyên thủy thuộc chuỗi phân rã Uranium ($^{238}\text{U}$), Thorium ($^{232}\text{Th}$) và đồng vị Kali ($^{40}\text{K}$). Trong đó, việc định lượng hoạt độ riêng của $^{226}\text{Ra}$ và $^{238}\text{U}$ trong đất đá, trầm tích và mẫu sinh học là tiêu chí hàng đầu để xác định nguy cơ phơi nhiễm khí phóng xạ Radon ($^{222}\text{Rn}$).

Chuỗi phân rã Uranium-238 tự nhiên:
238U (4.468 Ga) ---> 234Th (24.1 d) ---> 234Pa (1.17 m) ---> 234U (245.5 ka)
      ---> 226Ra (1600 a) ---> 222Rn (3.82 d) ---> 214Pb (26.8 m) ---> 214Bi (19.9 m)

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

Phương pháp phân tích phổ kế Gamma truyền thống thường dựa vào nguyên lý cân bằng thế kỷ (Secular Radioactive Equilibrium). Để xác định hoạt độ của $^{226}\text{Ra}$ và $^{238}\text{U}$, mẫu môi trường sau khi xử lý buộc phải được đóng gói kín trong các hộp tiêu chuẩn và lưu kho từ 20 đến 30 ngày (tương đương 5 đến 7 lần chu kỳ bán rã của $^{222}\text{Rn}$, $T_{1/2} = 3.82$ ngày) để các sản phẩm phân rã chuỗi sau gồm $^{214}\text{Pb}$ và $^{214}\text{Bi}$ đạt trạng thái cân bằng phóng xạ với $^{226}\text{Ra}$.

Quy trình cân bằng truyền thống tồn tại ba nhược điểm nghiêm trọng:

  1. Độ trễ thời gian cực lớn: Thời gian chờ 30 ngày làm tê liệt khả năng cảnh báo sớm khi xảy ra sự cố bức xạ và làm tăng chi phí lưu kho bảo quản mẫu.
  2. Nguy cơ sai số do rò rỉ khí Radon ($^{222}\text{Rn}$): Việc đóng kín mẫu đá hoặc bột đất nếu không đạt độ kín tuyệt đối ($100%$) sẽ gây thất thoát khí Radon, làm giảm số đếm ghi nhận tại các đỉnh gamma của $^{214}\text{Pb}$ ($351.9\text{ keV}$) và $^{214}\text{Bi}$ ($609.3\text{ keV}$, $1764.3\text{ keV}$), dẫn đến kết quả phân tích bị đánh giá thấp từ $15% - 40%$.
  3. Sự suy giảm chất lượng mẫu sinh học: Các mẫu rau củ, nước, thủy hải sản tươi sống bị phân hủy, thối rữa hoặc biến đổi cấu trúc hình học thể tích trong thời gian 30 ngày chờ cân bằng.

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

Dự án tập trung nghiên cứu, xây dựng và chuẩn hóa phương pháp xác định trực tiếp hoạt độ mẫu môi trường bằng hệ phổ kế gamma độ phân giải cao sử dụng đầu dò bán dẫn Germanium siêu tinh khiết (HPGe), nhằm đạt được các mục tiêu cụ thể:

  1. Xác định trực tiếp hoạt độ của $^{226}\text{Ra}$ ngay sau khi chuẩn bị mẫu thông qua đỉnh năng lượng $186.21\text{ keV}$, loại bỏ hoàn toàn thời gian chờ cân bằng phóng xạ 30 ngày.
  2. Thiết lập thuật toán phân tách sự đóng góp tỷ lệ số đếm tại vùng năng lượng $186\text{ keV}$ giữa $^{226}\text{Ra}$ ($186.21\text{ keV}$) và $^{235}\text{U}$ ($185.71\text{ keV}$) dựa trên dữ liệu hạt nhân mới nhất.
  3. Xây dựng đường cong hiệu suất đỉnh năng lượng toàn phần (Full-Energy Peak Efficiency - FEPE) của đầu dò HPGe sử dụng mẫu chuẩn đất $IAEA-RGU-1$, khống chế sai số thực nghiệm $< 8%$.
  4. Tính toán và hiệu chỉnh hệ số tự hấp thụ bức xạ gamma ($F_{abs}$) của mẫu thể tích dạng trụ dựa trên thành phần nguyên tố từ huỳnh quang tia X tán sắc bước sóng (WDXRF) và thư viện cơ sở dữ liệu NIST XCOM.
  5. Kiểm chứng thực nghiệm trên 5 mẫu đá thực tế, đánh giá độ sai biệt giữa phương pháp trực tiếp và phương pháp cân bằng truyền thống đạt mức tin cậy $< 20%$.

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

Phân tích hiện trạng và giải pháp kỹ thuật

Tiêu chí kỹ thuật Phương pháp đo cân bằng truyền thống Mô phỏng Monte Carlo thuần túy Phương pháp đo trực tiếp kết hợp hiệu chỉnh bán thực nghiệm (Đề xuất)
Thời gian phân tích $20 - 30\text{ ngày}$ (chờ cân bằng Radon) Vài giờ mô phỏng nhưng cần cấu hình chi tiết đầu dò $2 - 4\text{ giờ}$ đo trực tiếp ngay sau khi đóng mẫu
Ảnh hưởng rò rỉ $^{222}\text{Rn}$ Rất cao; sai số lớn nếu hộp mẫu hở Không phụ thuộc Radon nhưng nhạy với sai số kích thước Hoàn toàn không bị ảnh hưởng do đo trực tiếp $^{226}\text{Ra}$
Bảo toàn mẫu sinh học Kém; mẫu tươi dễ hỏng Không áp dụng trực tiếp Rất tốt; phân tích tức thì mẫu tươi
Độ phức tạp tính toán Thấp (chỉ tính diện tích đỉnh) Rất cao; đòi hỏi phần mềm MCNP/GEANT4 có bản quyền Vừa phải; kết hợp giải tích thực nghiệm và tra cứu NIST XCOM
Sai số tương đối ($%$) $5% - 15%$ (nếu kín khí) $5% - 10%$ (phụ thuộc dead layer) $< 10% - 15%$ sau khi hiệu chỉnh tự hấp thụ 
Ma trận yêu cầu kỹ thuật MoSCoW:

Thiết kế hệ thống đo và Technology Stack

   +-----------------------------------------------------------+
   |                BUỒNG CHÌ CHẮN PHÔNG THẤP                   |
   |  +-----------------------------------------------------+  |
   |  |           Mẫu môi trường (Hộp nhựa SG50)            |  |
   |  +-----------------------------------------------------+  |
   |                            | gamma ray                    |
   |  +-----------------------------------------------------+  |
   |  |  Đầu dò HPGe p-type (Bán dẫn Germanium siêu tinh khiết)  | |
   |  +-----------------------------------------------------+  |
   |                            | Tín hiệu điện tích (Preamplifier)
   +----------------------------|------------------------------+
                                v
     +----------------------------------------------------+
     |         Khối xử lý phổ kỹ thuật số (MCA)           |
     |         (Bộ khuếch đại tạo xung + ADC)             |
     +----------------------------------------------------+
                                | Dữ liệu phổ (SPE/CNF)
                                v
     +----------------------------------------------------+
     |             Phần mềm và Công cụ Tính toán          |
     |  * Thu nhận phổ: Genie 2000 v3.4 / Maestro v7.01    |
     |  * Tra cứu tiết diện tương tác: NIST XCOM v3.1     |
     |  * Phân tích & Tách đỉnh: Python 3.10 (SciPy/NumPy)|
     +----------------------------------------------------+

Ngăn chặn phông bức xạ môi trường

Hệ phổ kế gamma phông thấp sử dụng đầu dò HPGe đồng trục (Canberra GEM P4-83) được bố trí trong buồng chì che chắn có độ dày thành $100\text{ mm}$ chì nguyên chất, lót trong bằng $1.5\text{ mm}$ Cadmium và $1.0\text{ mm}$ Đồng nhằm hấp thụ tia X đặc trưng của chì ($K_\alpha \approx 74\text{ keV}$). Buồng đo được thổi khí Nitơ ($N_2$) lỏng bay hơi liên tục để đẩy toàn bộ khí Radon trong không khí ra khỏi không gian đo, duy trì mức phông môi trường cực thấp.

Cơ sở phương pháp luận và thuật toán toán học

  1. Phương pháp tách đỉnh chập $186\text{ keV}$: Tại năng lượng $186\text{ keV}$, phổ gamma ghi nhận sự trùng lặp giữa tia gamma $186.21\text{ keV}$ (xác suất phát $I_\gamma = 3.59%$) của $^{226}\text{Ra}$ và tia gamma $185.71\text{ keV}$ ($I_\gamma = 57.2%$) của $^{235}\text{U}$. Trong điều kiện mẫu tự nhiên đạt cân bằng đồng vị với tỷ lệ tự nhiên $^{235}\text{U}/^{238}\text{U} = 0.00720$ và $^{238}\text{U}$ cân bằng với $^{226}\text{Ra}$, tỷ lệ số đếm thực nghiệm đóng góp của $^{226}\text{Ra}$ tại đỉnh $186\text{ keV}$ được xác định bởi: $$k_{\text{Ra}} = \frac{I_{\gamma,\text{Ra}}}{I_{\gamma,\text{Ra}} + \left(\frac{A_{\text{U-235}}}{A_{\text{Ra-226}}}\right) \cdot I_{\gamma,\text{U-235}}} = \frac{0.0359}{0.0359 + 0.0461 \times 0.572} \approx 57.2%$$ Do đó, diện tích đỉnh thực của $^{226}\text{Ra}$ là: $$S_{\text{Ra-226}}(186.21) = 0.572 \times S_{\text{net}}(186\text{ keV})$$

  2. Công thức tính hoạt độ riêng ($A$): $$A = \frac{S_{\text{net}}}{T_{\text{live}} \cdot \varepsilon(E) \cdot I_\gamma \cdot m \cdot F_{\text{abs}}}$$ Trong đó:

  • $S_{\text{net}}$: Diện tích đỉnh năng lượng toàn phần thuần (đã trừ phông nền compton).
  • $T_{\text{live}}$: Thời gian đo sống thực tế (giây).
  • $\varepsilon(E)$: Hiệu suất đỉnh năng lượng toàn phần của đầu dò tại năng lượng $E$.
  • $I_\gamma$: Xác suất phát photon gamma trên một phân rã.
  • $m$: Khối lượng mẫu đo ($\text{kg}$).
  • $F_{\text{abs}}$: Hệ số hiệu chỉnh tự hấp thụ gamma của mẫu so với mẫu chuẩn.
  1. Hiệu chỉnh hệ số tự hấp thụ ($F_{abs}$): Đối với mẫu hình học trụ có bề dày $d$, hệ số suy giảm gamma tự hấp thụ được tính toán qua hệ số suy giảm tuyến tính $\mu(E)$: $$F_{\text{abs}} = \frac{1 - e^{-\mu(E) \cdot d}}{\mu(E) \cdot d}$$ Hệ số suy giảm khối lượng $(\mu/\rho)$ được tính toán chính xác bằng việc tổng xạ các nguyên tố thành phần hóa học phân tích từ WDXRF qua NIST XCOM: $$\left(\frac{\mu}{\rho}\right){\text{mẫu}} = \sum{i} w_i \cdot \left(\frac{\mu}{\rho}\right)_i$$

Implementation và kết quả

Quy trình phát triển và Code Snippet phân tích

Quy trình tự động hóa giải phổ, làm khớp đường cong hiệu suất và hiệu chỉnh tự hấp thụ được thực hiện thông qua module giải tích xây dựng trên nền tảng Python:

import numpy as np
from scipy.optimize import curve_fit

def efficiency_curve(energy, a, b, c, d):
    """
    Hàm chuẩn hóa hiệu suất đỉnh năng lượng toàn phần (FEPE)
    theo hàm log-polynomial bậc 3: ln(eff) = a + b*ln(E) + c*(ln(E))^2 + d*(ln(E))^3
    """
    ln_e = np.log(energy)
    return np.exp(a + b*ln_e + c*(ln_e**2) + d*(ln_e**3))

def calculate_self_absorption(mu_linear, thickness_cm):
    """
    Tính hệ số hiệu chỉnh tự hấp thụ F_abs cho mẫu hình học trụ
    mu_linear: Hệ số suy giảm tuyến tính (cm^-1) = (mu/rho) * rho
    thickness_cm: Chiều cao lớp mẫu trong hộp SG50 (cm)
    """
    if mu_linear * thickness_cm == 0:
        return 1.0
    return (1.0 - np.exp(-mu_linear * thickness_cm)) / (mu_linear * thickness_cm)

def direct_radium_activity(s_net_186, live_time, eff_186, mass_kg, f_abs, ratio=0.572):
    """
    Tính hoạt độ riêng của Ra-226 bằng phương pháp trực tiếp từ đỉnh 186 keV
    I_gamma của Ra-226 tại 186.21 keV là 0.0359
    """
    i_gamma_ra = 0.0359
    s_net_ra = s_net_186 * ratio
    activity = s_net_ra / (live_time * eff_186 * i_gamma_ra * mass_kg * f_abs)
    return activity

Dữ liệu thực nghiệm chuẩn hóa hiệu suất đỉnh (Mẫu chuẩn IAEA-RGU-1)

Mẫu chuẩn $IAEA-RGU-1$ ($m = 100.0\text{ g}$, $A(^{238}\text{U}) = 4940\text{ Bq/kg}$) được đo trong thời gian $54000\text{ s}$ để xây dựng đường chuẩn hiệu suất.

Đồng vị Năng lượng $E$ (keV) Xác suất phát $I_\gamma$ (%) Hiệu suất thực nghiệm $\varepsilon_{\text{exp}}$ Hiệu suất làm khớp $\varepsilon_{\text{fit}}$ Độ sai biệt $\Delta$ (%)
$^{234}\text{Th}$ $63.30$ $3.75$ $0.0312 \pm 0.0015$ $0.0321$ $+2.88$
$^{226}\text{Ra}$ $186.21$ $3.59$ $0.0524 \pm 0.0011$ $0.0518$ $-1.15$
$^{214}\text{Pb}$ $295.22$ $18.42$ $0.0398 \pm 0.0006$ $0.0401$ $+0.75$
$^{214}\text{Pb}$ $351.93$ $35.60$ $0.0335 \pm 0.0004$ $0.0332$ $-0.90$
$^{214}\text{Bi}$ $609.31$ $45.49$ $0.0195 \pm 0.0002$ $0.0197$ $+1.02$
$^{214}\text{Bi}$ $1120.29$ $14.92$ $0.0108 \pm 0.0002$ $0.0107$ $-0.93$
$^{214}\text{Bi}$ $1764.49$ $15.31$ $0.0069 \pm 0.0001$ $0.0071$ $+2.89$

Nhận xét: Độ sai biệt giữa giá trị hiệu suất thực nghiệm và hàm làm khớp log-polynomial đều nằm dưới $3%$, chứng minh độ tin cậy tuyệt đối của hàm chuẩn hiệu suất trên toàn dải năng lượng từ $60\text{ keV}$ đến $1800\text{ keV}$.

Kết quả đo đạc và hiệu chỉnh trên các mẫu đá phân tích

Dưới đây là kết quả phân tích hoạt độ riêng của $^{226}\text{Ra}$ trên 5 mẫu đá thực nghiệm ($M1$ đến $M5$) so sánh giữa phương pháp đo trực tiếp (đo ngay) và phương pháp đo cân bằng truyền thống (sau 30 ngày lưu mẫu):

Mẫu phân tích Hoạt độ trực tiếp chưa hiệu chỉnh $F_{abs}$ (Bq/kg) Hệ số tự hấp thụ $F_{abs}$ ($186\text{ keV}$) Hoạt độ trực tiếp sau hiệu chỉnh $F_{abs}$ (Bq/kg) Hoạt độ đo cân bằng 30 ngày (Bq/kg) Sai biệt sau hiệu chỉnh $\Delta$ (%)
Mẫu 1 (M1) $65.29 \pm 1.37$ $0.912$ $71.59 \pm 1.50$ $80.64 \pm 1.56$ $11.22$
Mẫu 2 (M2) $52.36 \pm 1.16$ $0.925$ $56.61 \pm 1.25$ $64.60 \pm 1.32$ $12.37$
Mẫu 3 (M3) $92.67 \pm 1.70$ $0.908$ $102.06 \pm 1.87$ $105.60 \pm 1.86$ $3.35$
Mẫu 4 (M4) $62.52 \pm 1.31$ $0.931$ $67.15 \pm 1.41$ $68.46 \pm 1.39$ $1.91$
Mẫu 5 (M5) $61.63 \pm 1.21$ $0.915$ $67.36 \pm 1.32$ $76.35 \pm 1.39$ $11.77$ 
Biểu đồ so sánh độ sai biệt hoạt độ Ra-226 (Trực tiếp vs Cân bằng):
Mẫu M1: [============>       ] 11.22%
Mẫu M2: [=============>      ] 12.37%
Mẫu M3: [===>                ]  3.35%  <-- Độ tương đồng cực cao
Mẫu M4: [==>                 ]  1.91%  <-- Độ tương đồng cực cao
Mẫu M5: [============>       ] 11.77%
Tất cả các mẫu đều đạt tiêu chuẩn chấp nhận sai số phóng xạ môi trường (< 20%).

Đổi mới và đóng góp khoa học

  1. Tối ưu hóa thời gian phân tích vượt bậc ($99.5%$): Giảm chu kỳ phân tích mẫu từ $720\text{ giờ}$ ($30\text{ ngày}$) xuống còn $2 - 4\text{ giờ}$ thao tác và đo trên hệ phổ kế, cho phép phản ứng nhanh trong các kịch bản kiểm tra an toàn phóng xạ khẩn cấp.
  2. Loại trừ sai số do thoát khí Radon: Khác với phương pháp đo gián tiếp qua $^{214}\text{Pb}/^{214}\text{Bi}$ vốn chịu ảnh hưởng nặng bởi độ kín của bao bì mẫu, phương pháp trực tiếp đo chính photon phát ra từ hạt nhân mẹ $^{226}\text{Ra}$, loại bỏ hoàn toàn nguy cơ sai số do rò rỉ khí gas.
  3. Mô hình hóa chính xác hiệu ứng tự hấp thụ thể tích: Ứng dụng tích hợp thành công dữ liệu phân tích hóa học nguyên tố WDXRF với cơ sở dữ liệu NIST XCOM, giúp bù trừ sai số suy giảm photon ở vùng năng lượng thấp ($63.3\text{ keV}$ của $^{234}\text{Th}$ và $186.21\text{ keV}$ của $^{226}\text{Ra}$), đưa độ sai lệch tổng thể về dưới $15%$.

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

Lộ trình triển khai ứng dụng (Implementation Roadmap):
Quý 1: Chuẩn hóa hệ đầu dò HPGe & Thư viện ma trận hấp thụ XCOM
Quý 2: Thử nghiệm phân tích mẫu địa chất khai khoáng & Đất hiếm
Quý 3: Tích hợp module tự động hóa xử lý phổ và trích xuất báo cáo
Quý 4: Chuyển giao quy trình cho các Trung tâm Quan trắc Môi trường

Các kịch bản ứng dụng thực tế

  • Giám sát an toàn mỏ khai khoáng và vật liệu xây dựng: Kiểm tra nhanh hàm lượng phóng xạ tự nhiên trong đá granit, quặng bauxite, xỉ than và cát xây dựng trước khi đưa vào sản xuất thương mại.
  • Quan trắc thực phẩm và nông sản xuất khẩu: Kiểm tra chỉ tiêu an toàn bức xạ trên các lô hàng nông sản tươi sống, thủy hải sản mà không làm thối rữa, biến chất mẫu trong quá trình chờ đợi phân tích.
  • Đánh giá địa tầng và thủy văn: Định lượng tỷ lệ mất cân bằng phóng xạ $^{238}\text{U}/^{226}\text{Ra}$ phục vụ nghiên cứu vận chuyển địa chất và xói mòn đất.

Phân tích Chi phí - Lợi ích (Cost-Benefit Analysis)

  • Tăng thông lượng phòng thí nghiệm: Nâng công suất phân tích từ $10\text{ mẫu/tháng/đầu dò}$ lên hơn $150\text{ mẫu/tháng/đầu dò}$ (tăng gấp 15 lần).
  • Tiết kiệm chi phí vận hành: Tiết kiệm hàng ngàn hộp kín chuyên dụng chống thoát radon và giảm $80%$ diện tích kho lưu mẫu chờ cân bằng, mang lại tỷ suất hoàn vốn đầu tư (ROI) thiết bị chỉ trong vòng 12 - 18 tháng.

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

Hạn chế kỹ thuật

  • Độ bất định ở vùng năng lượng rất thấp ($< 100\text{ keV}$): Tại đỉnh $63.3\text{ keV}$ của $^{234}\text{Th}$, hệ số hấp thụ photon thay đổi rất dốc theo thành phần nguyên tố nặng ($Z > 20$), dẫn đến sai số phép đo có thể lên tới $20% - 30%$ nếu mẫu có chứa nhiều hợp chất sắt hoặc kim loại nặng chưa được định lượng chính xác.
  • Giới hạn tỷ lệ làm giàu Uranium: Thuật toán phân tách tỷ lệ $57.2%$ tại đỉnh $186\text{ keV}$ chỉ áp dụng cho các mẫu có tỷ lệ đồng vị Uranium tự nhiên ($^{235}\text{U}/^{238}\text{U} = 0.72%$). Đối với các mẫu bị làm giàu hoặc cạn kiệt nhân tạo, cần phải bổ sung kênh đo đồng vị phụ trợ.

Hướng nghiên cứu tiếp theo

  • Nghiên cứu ứng dụng mã mô phỏng Monte Carlo MCNP6 hoặc GEANT4 để mô hình hóa chi tiết hiệu ứng trùng phùng thực (True Coincidence Summing - TCS) của các đồng vị phát cascade gamma ($^{214}\text{Bi}$, $^{214}\text{Pb}$).
  • Tự động hóa quá trình xác định thành phần nền thông qua kỹ thuật đo tán xạ Compton/Rayleigh trực tiếp trên phổ gamma nhằm loại bỏ bước phân tích WDXRF trung gian.

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

Ma trận phân phối giá trị đến các nhóm thụ hưởng:
  • Sinh viên ngành Vật lý hạt nhân/Kỹ thuật hạt nhân: Tiếp cận hệ thống tài liệu chuẩn mực về lý thuyết tương tác bức xạ, quy trình hiệu chuẩn thiết bị HPGe thực tế và phương pháp xử lý dữ liệu thực nghiệm nâng cao.
  • Kỹ sư vận hành hệ phổ kế tại các viện nghiên cứu: Ứng dụng quy trình đo trực tiếp để nâng cao năng suất xử lý mẫu khẩn cấp, tiết kiệm thời gian và tài nguyên phòng thí nghiệm.
  • Cơ quan quản lý an toàn bức xạ và Doanh nghiệp: Có công cụ kỹ thuật tin cậy để thẩm định nhanh các chỉ số an toàn bức xạ môi trường, đẩy nhanh tiến độ dự án khai khoáng và xử lý môi trường.

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 này là gì?

Hệ thống đòi hỏi đầu dò bán dẫn HPGe có độ phân giải năng lượng (FWHM) $\le 2.0\text{ keV}$ tại đỉnh $1332\text{ keV}$ của $^{60}\text{Co}$ và $\le 1.2\text{ keV}$ tại đỉnh $186\text{ keV}$. Buồng chắn chì phải đạt độ dày tối thiểu $10\text{ cm}$ với lớp lót $Cu/Cd$ để triệt tiêu phông bức xạ tự nhiên.

2. Khi nào phương pháp đo trực tiếp tại đỉnh 186 keV không còn chính xác?

Phương pháp sẽ xuất hiện sai số lớn nếu mẫu môi trường bị nhiễm bẩn bởi nguồn Uranium nhân tạo (đã qua làm giàu hoặc cạn kiệt $^{235}\text{U}$) vì tỷ lệ đóng góp $57.2%$ dựa trên giả định cân bằng đồng vị tự nhiên ($^{235}\text{U}/^{238}\text{U} = 0.0072$).

3. Làm thế nào để hiệu chỉnh tự hấp thụ nếu phòng thí nghiệm không có máy phân tích WDXRF?

Có thể sử dụng phương pháp đo gamma truyền qua (Transmission Method) bằng cách đặt một nguồn chuẩn dạng điểm ($^{152}\text{Eu}$ hoặc $^{133}\text{Ba}$) chiếu qua mẫu đo để xác định trực tiếp hệ số suy giảm tuyến tính thực nghiệm $\mu(E)$ mà không cần biết thành phần nguyên tố chi tiết.

4. Chi phí bảo trì hệ thống phổ kế HPGe hàng năm là bao nhiêu?

Hệ thống HPGe làm mát bằng Nitơ lỏng tiêu tốn khoảng 15 - 25 lít $N_2$ lỏng mỗi tuần (khoảng 20 - 30 triệu VNĐ/năm). Nếu sử dụng hệ thống làm mát điện cơ học tử (Cryo-electric cooler), chi phí vận hành sẽ chủ yếu là điện năng tiêu thụ liên tục.

5. Độ sai biệt dưới 20% giữa phương pháp trực tiếp và phương pháp cân bằng có đạt tiêu chuẩn quốc tế không?

Hoàn toàn đạt tiêu chuẩn. Theo các quy chuẩn so sánh liên phòng thí nghiệm của IAEA (IAEA Interlaboratory Comparisons), đối với phân tích mẫu phóng xạ môi trường mức phông thấp, độ sai biệt tương đối giữa các phương pháp nằm trong khoảng $\pm 15% - 20%$ được xếp loại chấp nhận hoàn toàn (Acceptable Performance).

Kết luận

Khóa luận tốt nghiệp đã chứng minh tính khả thi, độ tin cậy và hiệu quả vượt trội của phương pháp xác định trực tiếp hoạt độ mẫu môi trường bằng phổ kế gamma đầu dò HPGe. Bằng việc làm chủ kỹ thuật phân tách đỉnh năng lượng chập $186.21\text{ keV}$ của $^{226}\text{Ra}$ và kết hợp mô hình hiệu chỉnh tự hấp thụ thể tích bán thực nghiệm thông qua NIST XCOM, nghiên cứu đã giải quyết triệt để bài toán thắt nút cổ chai về thời gian chờ cân bằng 30 ngày của phương pháp truyền thống.

Kết quả thực nghiệm trên các mẫu chuẩn $IAEA-RGU-1$ và 5 mẫu đá thực tế cho thấy độ sai biệt giữa phương pháp đo trực tiếp và phương pháp đo cân bằng luôn duy trì ở mức $< 12.5%$ đối với $^{226}\text{Ra}$, khẳng định độ chính xác cao và tiềm năng ứng dụng thực tiễn rộng rãi trong công tác quan trắc môi trường và an toàn bức xạ hạt nhân.