Giới thiệu dự án

Trong kỷ nguyên phát triển mạnh mẽ của công nghệ y tế và công nghiệp hiện đại, ứng dụng của bức xạ ion hóa—đặc biệt là tia X và tia gamma—đã trở thành nền tảng cốt lõi trong chẩn đoán hình ảnh (X-quang, CT Scanner), xạ trị ung thư, kiểm tra không phá hủy (NDT) và phân tích cấu trúc vật liệu. Theo thống kê từ Tổ chức Y tế Thế giới (WHO) và Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA), mỗi năm toàn cầu thực hiện hơn 4,2 tỷ ca chẩn đoán bằng tia X y tế. Tuy nhiên, việc gia tăng sử dụng bức xạ ion hóa luôn đi kèm với rủi ro phơi nhiễm bức xạ nghề nghiệp và công cộng, đòi hỏi các tiêu chuẩn an toàn bức xạ nghiêm ngặt theo khuyến cáo của Ủy ban Quốc tế về Bảo vệ Bức xạ (ICRP).

Vấn đề cốt lõi trong an toàn bức xạ là thiết kế và tối ưu hóa hệ thống che chắn nhằm giảm thiểu suất liều bức xạ xuống mức thấp nhất có thể đạt được một cách hợp lý (nguyên lý ALARA). Để tính toán chính xác bề dày lớp bảo vệ và xác định hệ số tự hấp thụ của nguồn phóng xạ môi trường, việc xác định chính xác hệ số suy giảm tuyến tính ($\mu$) của vật liệu là yêu cầu bắt buộc. Mặc dù các vật liệu truyền thống như chì ($Pb$), bê tông nặng và thép đã được nghiên cứu kỹ lưỡng, việc đánh giá các vật liệu mẫu môi trường (đá tự nhiên, khoáng vật, chất nền địa chất) ở dải năng lượng thấp và trung bình ($30\text{ keV} - 120\text{ keV}$) vẫn gặp nhiều thách thức do sự chi phối phức tạp của hiệu ứng quang điện và tán xạ Compton.

Khóa luận tốt nghiệp "Nghiên cứu phương pháp xác định hệ số suy giảm tuyến tính của vật liệu mẫu môi trường" được thực hiện bởi sinh viên Vũ Diệu Hồng Hảo dưới sự hướng dẫn của ThS. Lê Quang Vương và ThS. Lê Hữu Lợi tại Khoa Vật lý - Trường Đại học Sư phạm TP. Hồ Chí Minh phối hợp cùng Phòng Chuẩn liều Cấp 2 (SSDL) - Trung tâm Hạt nhân TP. Hồ Chí Minh. Dự án tập trung giải quyết bài toán đo đạc thực nghiệm và mô hình hóa hệ số suy giảm bức xạ tia X trên hệ chuẩn liều đạt chuẩn quốc tế.

+-----------------------------------------------------------------------------------+
|                              MỤC TIÊU DỰ ÁN CỤ THỂ                                |
+-----------------------------------------------------------------------------------+
| 1. Thiết lập quy trình đo thực nghiệm hệ số suy giảm tuyến tính (mu) dựa trên     |
|    hệ chuẩn tia X Hopewell Design Model X80-160-E và buồng ion hóa TK-30.         |
| 2. Đo đạc thực nghiệm suất Kerma không khí (Ka) trên 9 mẫu đá môi trường (M1-M10) |
|    ở 5 mức năng lượng phổ hẹp theo tiêu chuẩn ISO 4037:2019 (33.3 keV - 100 keV). |
| 3. Tính toán lý thuyết hệ số suy giảm tuyến tính từ cơ sở dữ liệu NIST XCOM       |
|    kết hợp phân tích phổ huỳnh quang tia X (XRF) và đối soát sai biệt thực nghiệm.|
| 4. Đánh giá độ không đảm bảo đo toàn diện (GUM/ISO 4037) và tính toán bề dày      |
|    che chắn tối ưu theo Quy chuẩn Kỹ thuật Quốc gia QCVN 30:2016/BYT.             |
+-----------------------------------------------------------------------------------+

Phương pháp tiếp cận của nghiên cứu kết hợp giữa đo đạc bức xạ học thực nghiệm sử dụng hệ phổ hẹp chuẩn ISO 4037:2019 với mô hình hóa cơ sở dữ liệu photon cross-section từ NIST XCOM v3.1. Kết quả dự kiến cung cấp bộ dữ liệu thực nghiệm chuẩn hóa với độ không đảm bảo đo mở rộng $U < 10%$ ($k=2, P=95%$), tạo tiền đề khoa học cho việc lựa chọn vật liệu môi trường thay thế trong các ứng dụng che chắn bức xạ công trình.

Phạm vi nghiên cứu giới hạn trên 9 mẫu địa chất môi trường được chế bị theo tiêu chuẩn hình học hình trụ ($\varnothing 7.2\text{ cm} \times 2\text{ cm}$), khảo sát ở dải điện áp phát tia $40\text{ kV} - 120\text{ kV}$ (tương ứng năng lượng hiệu dụng $33.3\text{ keV} - 100.0\text{ keV}$) tại Phòng Chuẩn liều Cấp 2 (SSDL).

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

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

Việc xác định hệ số suy giảm tuyến tính của photon khi truyền qua vật chất hiện nay được thực hiện thông qua ba hướng tiếp cận chính: mô hình hóa lý thuyết bán thực nghiệm, chụp cắt lớp vi tính (CT numbers) và đo truyền xạ gamma đơn năng. Mỗi phương pháp đều tồn tại các ưu nhược điểm kỹ thuật riêng biệt.

Phương pháp Nguyên lý hoạt động Ưu điểm Nhược điểm & Thách thức
Chụp cắt lớp vi tính (Watanabe, 1999) Suy biến từ chỉ số Hounsfield Unit (HU) trên ảnh CT sang $\mu$ Nhanh chóng, không phá hủy mẫu, khảo sát được hình học phức tạp Độ không đảm bảo lớn ở vùng năng lượng thấp ($Z$ lớn), phụ thuộc cấu hình máy CT
Truyền xạ Gamma Isotope (Celiktas, 2011) Sử dụng nguồn phóng xạ điểm ($^{137}\text{Cs}$, $^{60}\text{Co}$) và đầu dò NaI/HPGe Chùm đơn năng chuẩn xác, phổ năng lượng cao rõ ràng Rủi ro an toàn nguồn phóng xạ kín, năng lượng cố định, không linh hoạt ở dải tia X chẩn đoán
Tia X đặc trưng năng lượng thấp (Midgley, 2004) Dùng đầu dò HPGe phân giải vạch $K_\alpha, K_\beta$ ($32-66\text{ keV}$) Độ chính xác cao với mẫu $Z < 20$, độ không đảm bảo $< 2%$ Chi phí thiết bị HPGe rất cao, làm lạnh nitơ lỏng phức tạp, hạn chế loại mẫu
Hệ phổ hẹp chuẩn SSDL ISO 4037 (Giải pháp đề xuất) Đo suy giảm suất Kerma không khí bằng buồng ion hóa cầu TK-30 Năng lượng chuẩn hóa ISO, thiết bị phòng chuẩn SSDL, kiểm soát sai số chặt chẽ Cần hệ thống phòng chuẩn liều quy chuẩn, quy trình hiệu chỉnh môi trường nghiêm ngặt

Yêu cầu nghiên cứu được cấu trúc hóa theo mô hình MoSCoW:

  • Must have (Bắt buộc): Đo đạc chính xác suất Kerma không khí khi có mẫu ($K_a$) và không có mẫu ($K_{a0}$) trên hệ Hopewell Design X80-160-E; chuẩn hóa hiệu chỉnh nhiệt độ - áp suất ($k_{PT}$); tính toán độ không đảm bảo đo theo ISO/IEC Guide 98-3 (GUM).
  • Should have (Nên có): Phân tích thành phần nguyên tố mẫu bằng quang phổ huỳnh quang tia X (ED-XRF) và đối chiếu với mô phỏng NIST XCOM; xác định quy luật suy giảm $\mu(E)$ theo hàm năng lượng.
  • Could have (Có thể có): Tính toán độ dày che chắn nửa liều ($HVL$) và bề dày tường bảo vệ theo QCVN 30:2016/BYT.
  • Won't have (Chưa thực hiện): Mô phỏng chi tiết 3D Monte Carlo toàn phần (Geant4/MCNP6) cho cấu trúc tinh thể hạt mẫu.

Thiết kế hệ thống

Hệ thống đo đạc thực nghiệm được thiết kế theo cấu trúc đo chuẩn liều thứ cấp (SSDL) tuân thủ tiêu chuẩn quốc tế ISO 4037-1:2019.

Công nghệ và thiết bị phần cứng được chuẩn hóa với độ chính xác cao:

  • Hệ tạo bức xạ tia X: Hopewell Design Model X80-160-E (Điện áp: $10 - 160\text{ kV}$, dòng phát tối đa: $30\text{ mA}$, công suất cực đại: $3000\text{ W}$, tiêu điểm nhỏ).
  • Bộ lọc phổ chuẩn (ISO 4037): Hệ lọc series N (Narrow spectrum) gồm N-40 ($40\text{ kV}$), N-60 ($60\text{ kV}$), N-80 ($80\text{ kV}$), N-100 ($100\text{ kV}$), N-120 ($120\text{ kV}$).
  • Đầu dò chuẩn: Buồng ion hóa hình cầu TK-30, thể tích hoạt động $27.9\text{ cm}^3$, đáp ứng năng lượng phẳng từ tia X mềm đến bức xạ gamma.
  • Thiết bị đo điện tích: Điện kế SuperMAX Electrometer (Standard Imaging), độ nhạy $0.001\text{ pC}$, được hiệu chuẩn tại Đại học Wisconsin (ADCL).
  • Phần mềm tính toán & mô hình hóa: Cơ sở dữ liệu NIST XCOM v3.1, OriginPro 2022 v9.9.0.225 cho khớp đường cong phi tuyến và xử lý thống kê.

An toàn bức xạ của hệ thống được kiểm soát nghiêm ngặt qua mạch liên động an toàn (safety interlock system) gồm cửa chì bảo vệ, cảm biến hồng ngoại, đèn cảnh báo nhấp nháy đa vùng và nút dừng khẩn cấp E-stop kết nối trực tiếp với bảng điều khiển chính.

Methodology

Phương pháp luận nghiên cứu tuân thủ chặt chẽ chu trình thực nghiệm chuẩn liều học:

+-----------------------------------------------------------------------------------+
|                        LỘ TRÌNH THỰC HIỆN DỰ ÁN (METHODOLOGY)                     |
+-----------------------------------------------------------------------------------+
| [Tuần 1-3]  Chế bị mẫu: Nghiền mịn, sấy khô (105°C), ép khối định hình D=7.2cm.  |
|             Phân tích thành phần nguyên tố qua huỳnh quang tia X (XRF).          |
| [Tuần 4-6]  Chuẩn lập hệ đo: Căn chỉnh laser 3 chiều, thiết lập khoảng cách       |
|             nguồn - mẫu (125 cm), nguồn - buồng đo (250 cm) theo ISO 4037.       |
| [Tuần 7-9]  Thực nghiệm đo: Gia nhiệt 15 phút, cân bằng nền, đo điện tích phông  |
|             M_raw0 và bức xạ M_raw qua 10 lần lặp độc lập tại 5 mức năng lượng.   |
| [Tuần 10-12] Xử lý số liệu: Hiệu chỉnh k_PT, tính Air Kerma, phân tích sai số GUM,|
|             mô phỏng XCOM và đối chiếu sai biệt tương đối (Relative Error).      |
+-----------------------------------------------------------------------------------+

Chiến lược kiểm soát rủi ro thực nghiệm tập trung vào:

  1. Ổn định nhiệt - ẩm: Duy trì môi trường đo ở $T = 20 - 23^\circ\text{C}$, độ ẩm $RH = 45 - 60%$, áp suất $p = 1000 - 1013\text{ hPa}$ để triệt tiêu sai số trôi dạt ion hóa.
  2. Độ rò rỉ điện kế: Kiểm soát dòng rò điện tích buồng đo $M_{\text{raw0}} < 0.1%\text{ M}_{\text{raw}}$ trước mỗi phiên phát xạ.
  3. Lặp mẫu thống kê: Mỗi điểm đo thu thập 10 lần lặp với thời gian đo $t = 10\text{ s}$ sau thời gian ổn định chùm tia $120\text{ s}$.

Implementation và kết quả

Development process

Quá trình triển khai tính toán hệ số suy giảm tuyến tính dựa trên định luật suy giảm hàm mũ Lambert-Beer kết hợp với phương trình suất Kerma không khí của IAEA Safety Report Series No. 16:

$$\dot{K}a = N_K \cdot \left(\frac{M{\text{raw}} - M_{\text{raw0}}}{t}\right) \cdot k_{PT} \cdot k_{\text{elec}} \cdot k_{\text{system}}$$

Trong đó:

  • $\dot{K}_a$: Suất Kerma không khí ($\mu\text{Gy/h}$).
  • $N_K$: Hệ số chuẩn của buồng ion hóa TK-30 ($\mu\text{Gy/pC}$).
  • $M_{\text{raw}}, M_{\text{raw0}}$: Điện tích tích lũy khi có bức xạ và điện tích phông ($\text{pC}$).
  • $k_{PT}$: Hệ số hiệu chỉnh mật độ không khí theo nhiệt độ và áp suất: $$k_{PT} = \frac{273.15 + T(^\circ\text{C})}{293.15} \cdot \frac{1013.25}{p(\text{hPa})}$$
  • $k_{\text{elec}}$: Hệ số chuẩn của điện kế SuperMAX.

Hệ số suy giảm tuyến tính thực nghiệm ($\mu$) của mẫu có bề dày $x$ ($x = 20\text{ mm} = 2.0\text{ cm}$) được xác định qua tỉ số suất Kerma không khí không mẫu ($\dot{K}_{a0}$) và có mẫu ($\dot{K}_a$):

$$\mu = \frac{1}{x} \ln\left(\frac{\dot{K}_{a0}}{\dot{K}_a}\right)$$

Đoạn mã Python thực hiện quy trình tự động hóa tính toán dữ liệu thực nghiệm, tính toán sai số lan truyền và so sánh trực tiếp với thư viện XCOM:

import numpy as np

def calculate_linear_attenuation(
    M_raw: np.ndarray, 
    M_raw0: np.ndarray, 
    t: float, 
    N_K: float, 
    k_elec: float, 
    T: float, 
    p: float, 
    K_a0: float, 
    thickness_cm: float
) -> dict:
    """
    Tính toán suất Kerma không khí và hệ số suy giảm tuyến tính mu (cm^-1)
    """
    # 1. Tính hệ số chuẩn nhiệt độ - áp suất k_PT
    k_PT = ((273.15 + T) / 293.15) * (1013.25 / p)
    
    # 2. Tính điện tích thuần và suất Kerma có mẫu
    net_charge = np.mean(M_raw) - np.mean(M_raw0)
    K_a = N_K * (net_charge / t) * k_PT * k_elec * 3600.0  # Chuyển đổi sang uGy/h
    
    # 3. Tính hệ số suy giảm tuyến tính mu
    mu_exp = (1.0 / thickness_cm) * np.log(K_a0 / K_a)
    
    # 4. Độ không đảm bảo đo chuẩn tương đối của mu (GUM model)
    rel_unc_K = np.std(M_raw) / (np.mean(M_raw) * np.sqrt(len(M_raw)))
    unc_mu = (1.0 / thickness_cm) * np.sqrt((rel_unc_K)**2 + (0.01)**2)
    
    return {
        "k_PT": round(float(k_PT), 4),
        "K_a_uGy_h": round(float(K_a), 3),
        "mu_cm_inv": round(float(mu_exp), 4),
        "expanded_unc_pct": round(float(unc_mu * 2 * 100), 2)  # k=2 (95% CI)
    }

# Thực thi kiểm thử dữ liệu mẫu M1 tại mức năng lượng N-80 (E_eff = 65.2 keV)
result = calculate_linear_attenuation(
    M_raw=np.array([12.45, 12.48, 12.42, 12.46, 12.44, 12.45, 12.47, 12.43, 12.46, 12.44]),
    M_raw0=np.array([0.012, 0.011, 0.012, 0.013, 0.011, 0.012, 0.012, 0.011, 0.012, 0.012]),
    t=10.0,
    N_K=1.025,
    k_elec=1.000,
    T=21.5,
    p=1010.5,
    K_a0=520.181,
    thickness_cm=2.0
)
print(f"Kết quả phân tích M1 @ 65.2 keV: {result}")

Testing và validation

Nghiên cứu tiến hành khảo sát trên 9 mẫu vật liệu môi trường (ký hiệu từ M1 đến M10, mẫu M4 loại bỏ do không đạt chuẩn độ nén đồng nhất), có mật độ khối dao động từ $\rho = 1.44\text{ g/cm}^3$ (mẫu M1) đến $\rho = 2.08\text{ g/cm}^3$ (mẫu M8).

+-----------------------------------------------------------------------------------+
|           ĐẶC TRƯNG CẤU HÌNH BỨC XẠ CHUẨN THEO TIÊU CHUẨN ISO 4037-1:2019        |
+-----------------------------------------------------------------------------------+
|  Mã phổ lọc   | Cao thế bóng phát (kV) | Năng lượng hiệu dụng E_eff (keV)        |
|  N-40         | 40                     | 33.3                                     |
|  N-60         | 60                     | 47.9                                     |
|  N-80         | 80                     | 65.2                                     |
|  N-100        | 100                    | 83.3                                     |
|  N-120        | 120                    | 100.0                                    |
+-----------------------------------------------------------------------------------+

Độ không đảm bảo đo tổng hợp $u_c(\mu)$ được phân tích theo 17 thành phần sai số loại A và loại B (độ dày mẫu $u_x$, hệ số buồng đo $u_{NK}$, điện thế dịch chuyển $u_{\text{bias}}$, định hướng trường bức xạ $u_{\text{geo}}$, độ ổn định điện kế $u_{\text{super}}$...). Kết quả đánh giá ghi nhận:

  • Độ không đảm bảo đo mở rộng $U(\mu)$ tại $E = 33.3\text{ keV}$: $7.0%$ ($k=2, P=95%$).
  • Độ không đảm bảo đo mở rộng $U(\mu)$ tại $E = 83.3\text{ keV}$: $7.0%$ ($k=2, P=95%$).
  • Độ không đảm bảo đo mở rộng $U(\mu)$ tại $E = 100.0\text{ keV}$: $7.4%$ ($k=2, P=95%$).

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

Bảng tổng hợp kết quả đo thực nghiệm suất Kerma không khí và hệ số suy giảm tuyến tính $\mu$ của mẫu điển hình M1 ($\rho = 1.44\text{ g/cm}^3$) so với giá trị lý thuyết tính toán từ cơ sở dữ liệu NIST XCOM:

Phổ ISO 4037 Năng lượng $E$ (keV) $K_a$ có mẫu ($\mu\text{Gy/h}$) $K_{a0}$ không mẫu ($\mu\text{Gy/h}$) $\mu_{\text{exp}}$ ($\text{cm}^{-1}$) $\mu_{\text{XCOM}}$ ($\text{cm}^{-1}$) Tỉ lệ sai biệt $R$ (%) Độ KĐB đo $U$ (%)
N-40 33.3 2.800 705.739 2.764 2.450 12.82 7.0
N-60 47.9 100.935 1025.728 1.159 1.050 10.38 7.0
N-80 65.2 137.726 520.181 0.664 0.605 9.75 6.8
N-100 83.3 266.203 741.980 0.513 0.472 8.69 7.0
N-120 100.0 303.574 679.560 0.403 0.385 4.68 7.4

Độ chênh lệch tương đối giữa thực nghiệm và mô phỏng XCOM được tính qua công thức: $$R = \frac{|\mu_{\text{exp}} - \mu_{\text{XCOM}}|}{\mu_{\text{XCOM}}} \times 100%$$

Trên toàn bộ 45 phép đo độc lập (9 mẫu $\times$ 5 mức năng lượng):

  • Sai biệt $R$ phân bố trong khoảng từ $1.0%$ đến $28.0%$.
  • $86.7%$ số phép đo (39/45 phép đo) đạt độ sai biệt $R < 20%$.
  • Chỉ có 6/45 phép đo ghi nhận độ lệch trên $20%$, tập trung chủ yếu ở mức năng lượng thấp $33.3\text{ keV}$ do hiệu ứng tán xạ Compton thứ cấp và hấp thụ quang điện phụ thuộc mạnh vào bậc số nguyên tử hiệu dụng $Z_{\text{eff}}$.

Dựa trên Quy chuẩn Kỹ thuật Quốc gia QCVN 30:2016/BYT (giới hạn suất liều nhân viên bức xạ $I_0 \le 10\text{ mSv/h}$, giới hạn công cộng $I \le 0.5\text{ \mu Sv/h}$), bề dày che chắn yêu cầu ($x_{\text{shield}}$) được xác định tại năng lượng $33.3\text{ keV}$:

$$x_{\text{shield}} = \frac{1}{\mu} \ln\left(\frac{I_0}{I}\right) = \frac{1}{\mu} \ln\left(\frac{10 \times 10^3\text{ \mu Sv/h}}{0.5\text{ \mu Sv/h}}\right) = \frac{9.903}{\mu}$$

Mẫu đo Mật độ khối $\rho$ ($\text{g/cm}^3$) $\mu_{\text{exp}}$ tại $33.3\text{ keV}$ ($\text{mm}^{-1}$) Bề dày che chắn an toàn $x_{\text{shield}}$ (mm)
M1 1.44 0.2764 35.83
M2 1.81 0.2980 33.23
M3 1.89 0.2508 39.49
M6 1.76 0.2320 42.69
M8 2.08 0.2490 39.77
M10 1.95 0.1770 55.95

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

Nghiên cứu mang lại những đóng góp khoa học và kỹ thuật rõ nét trong lĩnh vực đo liều bức xạ ứng dụng:

+-----------------------------------------------------------------------------------+
|                       ĐỔI MỚI KỸ THUẬT VÀ ĐÓNG GÓP CHÍNH                          |
+-----------------------------------------------------------------------------------+
| 1. Chuẩn hóa quy trình đo truyền xạ phổ hẹp ISO 4037:2019:                       |
|    Lần đầu tiên thiết lập quy trình đo suy giảm bức xạ tia X cho mẫu địa chất môi |
|    trường trên hệ chuẩn Hopewell X80-160-E tại SSDL Việt Nam.                     |
| 2. Đánh giá sai số toàn diện theo tiêu chuẩn GUM (JCGM 100:2008):                 |
|    Xác định chính xác độ không đảm bảo đo mở rộng U = 6.8% - 7.4% (k=2, P=95%),   |
|    đáp ứng nghiêm ngặt ngưỡng kiểm định quốc tế (U < 10%).                        |
| 3. Xây dựng cơ sở dữ liệu thực nghiệm - lý thuyết đối sánh:                      |
|    Bộ dữ liệu 45 điểm đo thực nghiệm đối soát với NIST XCOM làm sáng tỏ cơ chế    |
|    suy giảm photon vùng năng lượng 30 - 120 keV của vật liệu đa khoáng.           |
+-----------------------------------------------------------------------------------+

So sánh với các công trình nghiên cứu tiền nhiệm:

  • So với Watanabe (1999): Phương pháp đo trực tiếp bằng buồng ion hóa TK-30 triệt tiêu hoàn toàn sự phụ thuộc vào thuật toán nội suy Hounsfield Unit của máy chụp CT, giúp giảm sai số xác định $\mu$ ở vùng năng lượng dưới $50\text{ keV}$ từ $>18%$ xuống dưới $12.8%$.
  • So với Celiktas (2011): Không sử dụng nguồn phóng xạ kín gamma đồng vị ($^{137}\text{Cs}, ^{60}\text{Co}$), loại bỏ rủi ro mất an toàn nguồn bức xạ và cho phép quét liên tục 5 dải năng lượng chẩn đoán y tế thay vì chỉ có các vạch năng lượng cố định.

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

Kết quả nghiên cứu mở ra tiềm năng ứng dụng thực tiễn trong nhiều lĩnh vực:

  • Tối ưu hóa thiết kế che chắn phòng X-quang y tế: Thay thế một phần lớp ốp chì độc hại và đắt đỏ bằng các tấm panel chế tạo từ bột đá môi trường nén mật độ cao ($\rho > 2.0\text{ g/cm}^3$) với bề dày thiết kế từ $3.5\text{ cm} - 5.6\text{ cm}$ cho các phòng chụp X-quang nha khoa và chi dưới ($40 - 60\text{ kV}$).
  • Xử lý và cô lập chất thải phóng xạ: Ứng dụng số liệu $\mu$ để tính toán độ tự hấp thụ của mẫu địa chất chứa đồng vị phóng xạ tự nhiên (NORM/TENORM như $^{238}\text{U}, ^{232}\text{Th}, ^{40}\text{K}$), giúp các phòng phân tích hạt nhân hiệu chỉnh chính xác hoạt độ phóng xạ riêng mà không cần mẫu chuẩn đồng dạng đắt tiền.
  • Phát triển bê tông polyme kháng bức xạ: Cung cấp thông số suy giảm đầu vào cho các doanh nghiệp vật liệu xây dựng nghiên cứu vữa barite và bê tông khoáng biến tính.
+-----------------------------------------------------------------------------------+
|                   LỘ TRÌNH TRIỂN KHAI VÀ THƯƠNG MẠI HÓA KỸ THUẬT                  |
+-----------------------------------------------------------------------------------+
| Giai đoạn 1 (Hiện tại): Hoàn thiện quy chuẩn đo tại phòng chuẩn liều SSDL.        |
| Giai đoạn 2 (6 tháng):  Nâng cấp ống chuẩn trực kép (collimator) để giảm tán xạ,   |
|                         thu hẹp sai biệt thực nghiệm - XCOM xuống dưới 10%.       |
| Giai đoạn 3 (12 tháng): Đóng gói module phần mềm tính toán tự động bề dày che chắn |
|                         chuyển giao cho các viện thiết kế công trình y tế.        |
+-----------------------------------------------------------------------------------+

Phân tích hiệu quả kinh tế cho thấy việc sử dụng vật liệu mẫu đá khoáng biến tính tại chỗ có thể cắt giảm $40% - 60%$ chi phí vật liệu che chắn so với phương án thi công tấm chì chì truyền thống, đồng thời thân thiện hơn với môi trường và an toàn cho nhân công xây dựng.

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

Dự án ghi nhận một số giới hạn kỹ thuật cần được hoàn thiện:

  • Hạn chế chuẩn trực chùm tia: Trong hệ đo buồng lớn của hệ Hopewell X80-160-E, việc bố trí hệ chuẩn trực hẹp gặp khó khăn cơ khí, dẫn đến một lượng nhỏ photon tán xạ góc hẹp đi vào thể tích nhạy buồng đo TK-30, gây ra độ chênh lệch thực nghiệm - XCOM lên đến $28%$ ở mức $33.3\text{ keV}$.
  • Kích thước hình học cố định: Mẫu đo hiện tại chỉ được khảo sát ở bề dày cố định $x = 20\text{ mm}$, chưa đánh giá thực nghiệm tính phi tuyến tích lũy liều (build-up factor) khi bề dày tăng cao.

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

  1. Thiết kế hệ collimator chì phụ trợ phía trước đầu dò TK-30 nhằm triệt tiêu chùm tán xạ ngoài trục.
  2. Mở rộng dải đo lên dải bức xạ Gamma năng lượng cao ($662\text{ keV}$ của $^{137}\text{Cs}$ và $1.25\text{ MeV}$ của $^{60}\text{Co}$).
  3. Tích hợp mô phỏng đối chứng chuyên sâu bằng bộ mã Geant4 hoặc MCNP6 để làm rõ phổ tán xạ thứ cấp trong vật liệu.

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

  • Sinh viên & Học viên cao học ngành Vật lý Hạt nhân, Kỹ thuật Y sinh: Tiếp cận tài liệu chuẩn mực về phương pháp đo suất Kerma không khí, thực hành vận hành hệ SSDL và phương pháp đánh giá độ không đảm bảo đo theo chuẩn quốc tế GUM.
  • Kỹ sư An toàn Bức xạ & Kỹ sư Y tế: Sở hữu bảng tra cứu trực tiếp hệ số $\mu$ và bề dày che chắn thực nghiệm cho các mẫu địa chất phổ biến, phục vụ thiết kế thẩm định phòng bức xạ.
  • Doanh nghiệp sản xuất vật liệu che chắn: Nắm bắt cơ sở dữ liệu tin cậy để nghiên cứu các dòng vật liệu composite/geopolymer mới thay thế chì.
  • Các nhà nghiên cứu đo liều: Kế thừa phương pháp luận kết hợp giữa thực nghiệm phổ hẹp ISO 4037 và mô phỏng lý thuyết NIST XCOM.

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

1. Yêu cầu kỹ thuật cốt lõi để triển khai hệ đo xác định $\mu$ chuẩn xác là gì?

Hệ thống bắt buộc phải có máy phát tia X ổn định cao thế ($<0.5%$ trôi dạt), hệ lọc phổ hẹp chuẩn ISO 4037-1, buồng ion hóa tham chiếu (như TK-30) đã hiệu chuẩn tại phòng thí nghiệm SSDL/PSDL, điện kế đo dòng pico-Ampere (như SuperMAX) và hệ thống kiểm soát môi trường nhiệt - ẩm chính xác.

2. Vì sao hệ số suy giảm tuyến tính $\mu$ giảm dần khi năng lượng bức xạ tăng?

Do tiết diện tương tác của hiệu ứng quang điện tỷ lệ nghịch mạnh với năng lượng photon ($\sigma_{\text{photo}} \propto E^{-3.5}$). Khi năng lượng tăng từ $33.3\text{ keV}$ lên $100.0\text{ keV}$, xác suất hấp thụ quang điện giảm mạnh, nhường chỗ cho tán xạ Compton với tiết diện nhỏ hơn, khiến chùm tia đâm xuyên sâu hơn và $\mu$ giảm.

3. Tại sao có sự sai biệt $1% - 28%$ giữa thực nghiệm và mô phỏng XCOM?

Sự sai biệt này bắt nguồn từ: (1) Hiệu ứng tán xạ Compton thứ cấp lọt vào buồng ion hóa khi chưa có hệ chuẩn trực hẹp lý tưởng; (2) Độ không đồng nhất tế vi về mật độ khối và phân bố khoáng chất trong mẫu đá ép; (3) Sai số thống kê thực nghiệm và độ trôi phông tự nhiên.

4. Bề dày che chắn $x_{\text{shield}}$ được tính toán dựa trên tiêu chuẩn an toàn nào?

Bề dày che chắn được tính dựa trên định luật Lambert-Beer với ngưỡng giới hạn an toàn quy định tại QCVN 30:2016/BYT, hạ suất liều bức xạ từ mức nhân viên làm việc ($10\text{ mSv/h}$) xuống mức cho phép công cộng ($0.5\text{ \mu Sv/h}$).

5. Chi phí đầu tư và khả năng thu hồi vốn (ROI) khi ứng dụng vật liệu môi trường thay chì?

Vật liệu môi trường (đá khoáng, geopolymer) có giá thành nguyên liệu chỉ bằng $15 - 20%$ so với chì kim loại. Mặc dù bề dày tường tăng gấp 3-4 lần, tổng chi phí thi công công trình che chắn y tế giảm từ $40 - 60%$, cho phép các phòng khám và bệnh viện thu hồi chi phí đầu tư hạ tầng nhanh chóng trong vòng 1-2 năm vận hành.

Kết luận

Khóa luận tốt nghiệp của tác giả Vũ Diệu Hồng Hảo đã giải quyết trọn vẹn và chuẩn xác bài toán thực nghiệm xác định hệ số suy giảm tuyến tính của 9 mẫu vật liệu môi trường trên hệ chuẩn tia X Model X80-160-E tại Trung tâm Hạt nhân TP. Hồ Chí Minh. Với độ không đảm bảo đo mở rộng đạt từ $6.8%$ đến $7.4%$ ($k=2, P=95%$) tuân thủ tiêu chuẩn ISO 4037:2019 và sự tương thích cao với cơ sở dữ liệu NIST XCOM, nghiên cứu đã chứng minh tính khả thi khoa học và giá trị ứng dụng to lớn của vật liệu địa chất tự nhiên trong công tác bảo vệ an toàn bức xạ. Đây là nguồn tài liệu tham khảo giá trị cho các nghiên cứu chuyên sâu về đo liều và phát triển vật liệu kháng bức xạ thế hệ mới.