Giới thiệu dự án
Trong lĩnh vực kiểm tra không phá hủy (NDT - Non-Destructive Testing) và vật lý hạt nhân ứng dụng, việc xác định các đặc trưng tương tác photon của vật liệu đóng vai trò then chốt trong thiết kế che chắn bức xạ, tính toán liều xạ trị y khoa và kiểm chuẩn quy trình công nghiệp. Đối với các đơn chất (như $Pb, Al, Cu$), nguyên tử số $Z$ là một đại lượng cố định. Tuy nhiên, trong thực tế công nghiệp và hóa chất, phần lớn môi trường tương tác là các hợp chất hoặc dung dịch phức tạp. Một nguyên tử số đơn lẻ không thể đại diện cho toàn bộ hợp chất; do đó, khái niệm nguyên tử số hiệu dụng ($Z_{eff}$) được thiết lập để mô tả tương đương khả năng tương tác bức xạ của hỗn hợp đa nguyên tố.
Đặc biệt, các dung dịch axit vô cơ và hữu cơ có mức độ độc hại, tính ăn mòn cao và dễ bay hơi, gây nguy hiểm nghiêm trọng cho kỹ thuật viên khi sử dụng các phương pháp đo tỉ trọng cổ điển (bình đo tỉ trọng pycnometer hoặc tỷ trọng kế hở). Đồ án tập trung giải quyết bài toán: Xác định nguyên tử số hiệu dụng $Z_{eff}$ và mật độ khối $\rho$ của dung dịch chất lỏng (trọng tâm là các hệ axit) thông qua phương pháp gamma truyền qua kết hợp mô phỏng Monte Carlo MCNP6 và thực nghiệm phổ kế gamma.
+------------------+ Chùm tia gamma chuẩn trực +--------------------+ Tín hiệu xung +--------------------+
| Nguồn Cs-137 | ----------------------------------> | Ống mẫu chất lỏng | ---------------------> | Đầu dò NaI(Tl) |
| (E = 661.7 keV) | | (Axit vô cơ/hữu cơ)| | (76.2 x 76.2 mm) |
+------------------+ +--------------------+ +--------------------+
|
v
+----------------------------+ +--------------------+
| Đánh giá Zeff, mật độ rho | <------------------ | Phân tích đa kênh |
| (Colegram / NIST / Origin) | Phổ PHS (8192 ch) | MCA & ADMCA |
+----------------------------+ +--------------------+
Mục tiêu nghiên cứu cụ thể
- Xác định $Z_{eff}$ lý thuyết: Tính toán giải tích $Z_{eff}$ của các dung dịch axit thông qua phương pháp trực tiếp và phương pháp tỉ số tiết diện nguyên tử/electron dựa trên cơ sở dữ liệu tiêu chuẩn NIST.
- Khảo sát mô phỏng Monte Carlo: Xây dựng mô hình hình học chi tiết của hệ đo gamma truyền qua trên mã nguồn MCNP6 phiên bản 1.0 với $6 \times 10^9$ đến $6 \times 10^{12}$ lịch sử hạt, khảo sát 24 loại axit (vô cơ và hữu cơ) có dải mật độ từ $0,43\text{ g/cm}^3$ đến $1,84\text{ g/cm}^3$.
- Đo đạc và kiểm chứng thực nghiệm: Thiết lập hệ phổ kế gamma sử dụng đầu dò nhấp nháy $NaI(Tl)$ ($76,2\text{ mm} \times 76,2\text{ mm}$) và nguồn phóng xạ $^{137}Cs$ ($661,7\text{ keV}$, hoạt độ $0,76 \pm 0,03\ \mu Ci$), đo phổ truyền qua của 6 loại axit thực tế và nước cất.
- Xây dựng đường chuẩn mật độ không phá hủy: Thiết lập hàm tương quan bán thực nghiệm $\ln(R) = f(\rho)$ cho phép xác định mật độ chất lỏng bên trong bình chứa kín mà không cần tiếp xúc trực tiếp.
Phân tích và thiết kế giải pháp
Phân tích hiện trạng
Trong phân tích chất lỏng và kiểm định cấu trúc vật liệu, các phương pháp truyền thống bộc lộ nhiều hạn chế về an toàn và độ xuyên sâu:
| Phương pháp kiểm tra |
Cơ chế hoạt động |
Ưu điểm |
Nhược điểm / Hạn chế |
| Đo tỉ trọng cổ điển (Pycnometer / Hydrometer) |
Đo thể tích và cân khối lượng |
Độ chính xác cao, chi phí thiết bị thấp |
Tiếp xúc trực tiếp; nguy cơ rò rỉ hóa chất độc hại, ăn mòn ($H_2SO_4, HNO_3$); sai số do bay hơi |
| Siêu âm (UT) |
Phân tích sóng phản xạ bề mặt/biên |
Chi phí thấp, không dùng nguồn phóng xạ |
Khó xuyên sâu, chịu ảnh hưởng lớn bởi bọt khí và độ nhớt chất lỏng |
| Chất lỏng thẩm thấu (PT) |
Ngấm mao dẫn và hiện màu |
Phát hiện nứt bề mặt nhanh |
Chỉ áp dụng cho bề mặt hở, vật liệu không xốp, không phân tích tính chất khối |
| Gamma truyền qua (Phương pháp đề xuất) |
Độ suy giảm photon theo luật Beer-Lambert |
Đo xuyên thấu, không tiếp xúc, kiểm tra bình kín thời gian thực |
Đòi hỏi che chắn bức xạ an toàn và hiệu chuẩn hình học chuẩn trực |
Phân tích yêu cầu hệ thống theo mô hình MoSCoW:
- Must have: Hệ chuẩn trực chì đôi (Collimator) cho nguồn và detector để đảm bảo chùm tia hẹp ($B \approx 1$); thuật toán khớp đỉnh Gauss và trừ nền One Step trên Colegram; bộ dữ liệu tiết diện photon chuẩn từ NIST.
- Should have: Mô phỏng MCNP6 với hàm mở rộng phổ FWHM thực nghiệm; hàm tương quan tuyến tính $R_{MP}$ để nội suy mật độ.
- Could have: Mở rộng dải năng lượng đa nguồn ($^{241}Am, ^{133}Ba, ^{60}Co$) để khảo sát $Z_{eff}$ phụ thuộc năng lượng.
- Won't have (trong phạm vi khóa luận): Tái tạo ảnh cắt lớp 3D (Computed Tomography) của dòng chảy động.
Thiết kế hệ thống
[ Khối Chì Che Chắn Nguồn: 24 x 22 x 5 cm ]
+----------------------------------------+
| Nguồn Cs-137 (661.7 keV, 0.76 uCi) |
+-------------------+--------------------+
| Chuẩn trực nguồn (r = 0.251 cm)
v
[ Tia Gamma Hẹp (I_0) ]
|
v (Khoảng cách d_1 = 12.4 cm)
+-------------------+--------------------+
| Ống Đong Thủy Tinh (Mẫu Axit Khảo Sát) |
| r_in = 1.340 cm, d_out = 1.520 cm |
+-------------------+--------------------+
|
v (Khoảng cách d_2 = 12.4 cm)
[ Tia Gamma Truyền Qua (I) ]
|
v Chuẩn trực đầu dò (r = 0.500 cm)
+-------------------+--------------------+
| Đầu dò NaI(Tl) (Amptek 76.2 x 76.2 mm) |
+----------------------------------------+
|
v
+----------------------------------------+
| Bộ Phân Tích Đa Kênh MCA (8192 Kênh) |
+----------------------------------------+
Hệ thống được thiết kế đồng trục tuyệt đối bằng phương pháp căn chỉnh laser quang học bước sóng đỏ (vùng khả kiến). Các thông số kỹ thuật cốt lõi:
- Nguồn bức xạ: $^{137}Cs$, phát photon gamma đơn năng $661,7\text{ keV}$ với xác suất phát $84,99%$.
- Hệ chuẩn trực (Collimation): Khối chì bảo vệ kích thước $24\text{ cm} \times 22\text{ cm} \times 5\text{ cm}$. Bán kính lỗ chuẩn trực nguồn $r_1 = 0,251\text{ cm}$, bán kính lỗ chuẩn trực đầu dò $r_2 = 0,500\text{ cm}$.
- Bộ ghi nhận tín hiệu: Đầu dò nhấp nháy $NaI(Tl)$ kích thước $76,2\text{ mm} \times 76,2\text{ mm}$ (Amptek), kết nối khối xử lý MCA cấu hình 8192 kênh năng lượng thông qua phần mềm ADMCA.
- Ống chứa mẫu: Ống thủy tinh chuyên dụng dung tích $100\text{ ml}$, bán kính trong $r = 1,340 \pm 0,009\text{ cm}$, độ dày thành ống được xác định chính xác để trừ suy giảm nền.
Phương pháp luận (Methodology)
Quy trình nghiên cứu áp dụng phương pháp nghiên cứu kết hợp song song giữa tính toán giải tích, mô phỏng ngẫu nhiên Monte Carlo và kiểm chứng thực nghiệm:
[Phổ Thực Nghiệm (ADMCA 8192 ch)] [Mô Phỏng MCNP6 (6x10^12 hạt)]
| |
v v
[Khớp phổ Colegram (Gauss+Step)] [Phổ PHS & FWHM thực nghiệm]
| |
+------------------+-------------------+
|
v
[Tính toán diện tích đỉnh toàn phần S]
|
v
[Hệ số suy giảm khối mu_m = (1/rho*x)*ln(I_0/I)]
|
v
[Tính tiết diện sigma_a, sigma_e -> Z_eff,II = sigma_a/sigma_e]
|
v
[So sánh đối chuẩn NIST & Xây dựng đường chuẩn R_MP]
- Kiểm soát rủi ro bức xạ: Tuân thủ nguyên tắc ALARA (As Low As Reasonably Achievable), đo phông môi trường độc lập, vận hành hệ thống chì che chắn dày $5\text{ cm}$ hấp thụ bức xạ tán xạ bên ngoài.
- Khử sai số tán xạ nhiều lần: Duy trì chùm tia hẹp song song, cho phép xấp xỉ hệ số tích lũy $B(E, Z) \approx 1$ trong định luật Beer-Lambert.
Implementation và kết quả
Quá trình triển khai và thuật toán xử lý
1. Định luật thực nghiệm Beer-Lambert và hệ số suy giảm khối
Khi chùm photon hẹp đơn năng đi qua môi trường vật chất có mật độ $\rho$ và bề dày $x$:
$$I = I_0 e^{-\mu x} = I_0 e^{-\mu_m \rho x} \implies \mu_m = \frac{1}{\rho x} \ln\left(\frac{I_0}{I}\right)$$
Trong đó:
- $I_0, I$: Cường độ chùm photon trước và sau khi qua mẫu (diện tích đỉnh hấp thụ năng lượng toàn phần FEA).
- $\mu$: Hệ số suy giảm tuyến tính ($\text{cm}^{-1}$).
- $\mu_m = \frac{\mu}{\rho}$: Hệ số suy giảm khối ($\text{cm}^2/\text{g}$).
2. Thuật toán xác định nguyên tử số hiệu dụng $Z_{eff}$
Phương pháp tính $Z_{eff}$ theo tiết diện tổng hiệu dụng ($\sigma_a$) và tiết diện electron hiệu dụng ($\sigma_e$):
$$\sigma_a = \frac{\mu_{m, Compound}}{N_A \sum_{i} \frac{w_i}{A_i}} = \frac{1}{N_A} \sum_{i} n_i A_i (\mu_m)_i$$
$$\sigma_e = \frac{1}{N_A} \sum_{i} \frac{f_i A_i}{Z_i} (\mu_m)i = \frac{\sigma_a}{Z{eff}}$$
$$\implies Z_{eff, II} = \frac{\sigma_a}{\sigma_e}$$
Trong đó $w_i, f_i, A_i, Z_i, n_i$ lần lượt là phần trăm khối lượng, tỉ lệ số nguyên tử, khối lượng nguyên tử, nguyên tử số và chỉ số nguyên tử của nguyên tố thứ $i$; $N_A = 6,022 \times 10^{23}\text{ mol}^{-1}$.
3. Mô hình tham số độ phân giải năng lượng trên MCNP6
Độ rộng nửa chiều cao cực đại (FWHM) được nạp vào thẻ hàm Gauss của MCNP6 để làm khớp đáp ứng phổ của detector $NaI(Tl)$:
$$FWHM(E) = a + b\sqrt{E + cE^2}$$
Với các hệ số thực nghiệm tối ưu hóa:
$$a = -0,00526\text{ MeV}, \quad b = 0,06096\text{ MeV}^{1/2}, \quad c = 0,01401\text{ MeV}^{-1}$$
4. Thuật toán phân tích phổ trên Colegram
Colegram sử dụng giải thuật phi tuyến Marquardt-Levenberg để khớp hàm dạng đỉnh Gaussian $G(x)$ chồng trên hàm nền One-Step $STEP(x)$:
$$G(x) = \frac{A}{\sigma \sqrt{2\pi}} \exp\left(-\frac{(x - x_0)^2}{2\sigma^2}\right)$$
$$STEP(x) = B + \frac{S}{2} \text{erfc}\left(\frac{x - x_0}{\sigma \sqrt{2}}\right)$$
// Logic thuật toán tính toán Z_eff và truyền sai số thực nghiệm
double calculate_Zeff(double mu_m_compound, const std::vector<Element>& elements) {
double sum_w_over_A = 0.0;
double sum_f_A_mu_over_Z = 0.0;
const double N_A = 6.02214076e23; // Avogadro constant
for (const auto& el : elements) {
sum_w_over_A += el.w / el.A;
sum_f_A_mu_over_Z += (el.f * el.A * el.mu_m) / el.Z;
}
double sigma_a = mu_m_compound / (N_A * sum_w_over_A);
double sigma_e = (1.0 / N_A) * sum_f_A_mu_over_Z;
return sigma_a / sigma_e; // Tra ve gia tri Z_eff
}
Kiểm thử và đánh giá độ chính xác
Thời gian đo mỗi mẫu thực nghiệm được thiết lập cố định ở $t = 28.800\text{ giây}$ (8 giờ) ở nhiệt độ phòng ổn định $24^\circ\text{C}$ nhằm đảm bảo số đếm thống kê đỉnh $N > 60.000\text{ xung}$, giảm sai số thống kê Poisson ($\sigma_N = \sqrt{N}$) xuống dưới $0,4%$.
Sai số bán kính trong của ống đo được đo lặp 6 lần bằng thước kẹp panme có độ chia $0,02\text{ mm}$:
$$r = 1,340 \pm 0,00898\text{ cm} \quad (\Delta r / r \approx 0,67%)$$
Sai số tương đối ($RD%$) được đánh giá theo công thức chuẩn:
$$RD(%) = \frac{|X_{Calculated} - X_{NIST}|}{X_{NIST}} \times 100%$$
Kết quả đạt được
1. Hệ số suy giảm khối lượng thực nghiệm và độ sai biệt với NIST ($E = 661,7\text{ keV}$)
| STT |
Loại mẫu Axit |
Mật độ $\rho$ ($\text{g/cm}^3$) |
Số đếm thực nghiệm $N$ |
$\mu_{m, NIST}$ ($\text{cm}^2/\text{g}$) |
$\mu_{m, TN}$ ($\text{cm}^2/\text{g}$) |
Độ sai biệt $RD$ (%) |
| 1 |
Nước cất ($H_2O$) |
1,000 |
$72571 \pm 695$ |
0,086 |
$0,085 \pm 0,006$ |
0,58% |
| 2 |
Nitric ($HNO_3$) |
1,513 |
$67690 \pm 671$ |
0,078 |
$0,080 \pm 0,005$ |
2,05% |
| 3 |
Formic ($CH_2O_2$) |
1,220 |
$70488 \pm 685$ |
0,080 |
$0,079 \pm 0,005$ |
2,18% |
| 4 |
Sulfuric ($H_2SO_4$) |
1,839 |
$62406 \pm 645$ |
0,079 |
$0,077 \pm 0,004$ |
2,20% |
| 5 |
Acetic ($C_2H_4O_2$) |
1,045 |
$71884 \pm 692$ |
0,082 |
$0,085 \pm 0,006$ |
3,23% |
| 6 |
Lactic ($C_3H_6O_3$) |
1,206 |
$70670 \pm 686$ |
0,082 |
$0,079 \pm 0,005$ |
4,14% |
| 7 |
Phosphoric ($H_3PO_4$) |
1,390 |
$61534 \pm 640$ |
0,079 |
$0,087 \pm 0,004$ |
10,09% |
2. Kết quả đối sánh nguyên tử số hiệu dụng $Z_{eff}$
| STT |
Mẫu chất lỏng |
Công thức |
$Z_{eff}$ (NIST) |
$Z_{eff}$ Thực nghiệm |
$RD_{TN}$ (%) |
$Z_{eff}$ Mô phỏng MCNP6 |
$RD_{MP}$ (%) |
| 1 |
Nước cất |
$H_2O$ |
3,33 |
$3,31 \pm 0,09$ |
0,57% |
3,12 |
6,30% |
| 2 |
Axit Nitric |
$HNO_3$ |
6,43 |
$6,42 \pm 0,07$ |
0,07% |
5,87 |
6,69% |
| 3 |
Axit Formic |
$CH_2O_2$ |
4,80 |
$4,70 \pm 0,07$ |
2,02% |
4,49 |
6,34% |
| 4 |
Axit Phosphoric |
$H_3PO_4$ |
6,19 |
$6,04 \pm 0,05$ |
2,44% |
5,76 |
6,80% |
| 5 |
Axit Lactic |
$C_3H_6O_3$ |
4,00 |
$4,13 \pm 0,09$ |
3,24% |
3,74 |
6,42% |
| 6 |
Axit Acetic |
$C_2H_4O_2$ |
4,00 |
$3,83 \pm 0,07$ |
4,28% |
3,75 |
6,14% |
| 7 |
Axit Sulfuric |
$H_2SO_4$ |
7,14 |
$7,87 \pm 0,06$ |
10,17% |
6,65 |
6,84% |
So sánh giá trị Zeff thực nghiệm vs NIST:
Axit Nitric : [TN: 6.42] vs [NIST: 6.43] ---> Lệch 0.07% (Trùng khớp xuất sắc)
Nước cất : [TN: 3.31] vs [NIST: 3.33] ---> Lệch 0.57%
Axit Formic : [TN: 4.70] vs [NIST: 4.80] ---> Lệch 2.02%
Axit Phosphoric: [TN: 6.04] vs [NIST: 6.19] ---> Lệch 2.44%
Axit Lactic : [TN: 4.13] vs [NIST: 4.00] ---> Lệch 3.24%
Axit Acetic : [TN: 3.83] vs [NIST: 4.00] ---> Lệch 4.28%
Axit Sulfuric : [TN: 7.87] vs [NIST: 7.14] ---> Lệch 10.17%
3. Đường chuẩn mật độ chất lỏng theo tỉ số đếm $R$
Khảo sát tỉ số số đếm chuẩn hóa theo nước:
$$R = \frac{N_{Acid}}{N_{H_2O}} \implies \ln(R) = \ln\left(\frac{N_{Acid}}{N_{H_2O}}\right)$$
Dữ liệu mô phỏng trên 22 hợp chất axit mở rộng xác lập quan hệ tuyến tính chặt chẽ giữa $\ln(R)$ và mật độ $\rho$ ($\text{g/cm}^3$):
- Dải chất lỏng nhẹ ($\rho = 0,781 - 0,894\text{ g/cm}^3$: Oleic, Palmitic, Montanic...): $R \approx 1,02 - 1,05; \ln(R) > 0$.
- Dải chất lỏng trung bình ($\rho = 0,990 - 1,051\text{ g/cm}^3$: Myristic, Acetic, Acrylic...): $R \approx 1,00; \ln(R) \approx 0,00$.
- Dải axit đậm đặc ($\rho = 1,170 - 1,830\text{ g/cm}^3$: Nitric, Phosphoric, Perchloric, Sulfuric): $R \approx 0,97 - 0,87; \ln(R) = -0,03 \text{ đến } -0,14$.
Phương trình đường chuẩn cho phép nội suy trực tiếp mật độ chất lỏng thực nghiệm với độ tin cậy $99,99%$ ($k = 2,58$).
Đổi mới và đóng góp
- Phương pháp luận đo mật độ không tiếp xúc cho chất ăn mòn mạnh: Giải quyết triệt để vấn đề rủi ro hóa chất độc hại bằng kỹ thuật gamma truyền qua; không cần mở nắp ống chứa, cho phép đo inline trên đường ống công nghiệp.
- Cơ sở dữ liệu $Z_{eff}$ mở rộng cho 24 loại axit: Cung cấp bộ thông số thực nghiệm và mô phỏng chi tiết cho các axit hữu cơ chuỗi dài (Arachidic, Behenic, Cerotic, Lignoceric...) và axit vô cơ đậm đặc vốn rất hiếm dữ liệu trong các tài liệu quốc tế trước đây.
- Mô hình hóa chính xác hệ phổ kế NaI(Tl) trên MCNP6: Tích hợp thành công hàm làm rộng năng lượng $FWHM(E)$ thực nghiệm vào mô phỏng hạt Monte Carlo, giúp độ lệch giữa mô phỏng và NIST duy trì ổn định dưới $8%$ trên toàn dải mật độ.
- Kiểm chứng tính hợp lệ vật lý của $Z_{eff}$: Mọi giá trị $Z_{eff}$ xác định được đều thỏa mãn điều kiện biên vật lý nghiêm ngặt $Z_{min} < Z_{eff} < Z_{max}$ của các nguyên tố cấu thành.
Ứng dụng thực tế và triển khai
Tình huống ứng dụng thực tế (Use Cases)
[Nhà máy hóa chất / Lọc dầu] [Hệ thống kiểm chuẩn y sinh / Che chắn xạ]
| |
v v
+------------------------------------+ +------------------------------------+
| Kiểm chuẩn nồng độ Axit tự động | | Thiết kế vật liệu cản xạ lỏng |
| Giám sát inline bồn chứa kín | | Đánh giá liều hấp thụ xạ trị |
| Phân loại phân đoạn dầu khí | | Mô phỏng tương đương mô người |
+------------------------------------+ +------------------------------------+
- Giám sát nồng độ axit trong công nghiệp hóa chất: Kiểm tra tự động độ tinh khiết và nồng độ của các bồn chứa axit vô cơ ($H_2SO_4, HNO_3$) trong dây chuyền sản xuất phân bón, bán dẫn mà không làm gián đoạn dòng chảy.
- Phân loại nhiên liệu và sản phẩm dầu mỏ: Ứng dụng đường chuẩn $R - \rho$ để nhận diện các phân đoạn hydrocacbon, dầu bôi trơn, phát hiện lẫn nước hoặc tạp chất tỉ trọng thấp.
- Tính toán liều xạ trị và thiết kế tấm chắn phóng xạ y tế: Cung cấp giá trị $Z_{eff}$ chính xác cho các dung dịch tương đương mô (tissue-equivalent liquids), hỗ trợ lập kế hoạch xạ trị gamma/X-ray.
Lộ trình và yêu cầu triển khai kỹ thuật
- Phần cứng yêu cầu: Nguồn chuẩn $^{137}Cs$ hoạt độ thấp ($< 1\ \mu Ci$, đảm bảo an toàn bức xạ cơ bản); đầu dò nhấp nháy $NaI(Tl)$ hoặc $LaBr_3(Ce)$; khối chuẩn trực chì định hình chùm tia hẹp ($d \le 5\text{ mm}$).
- Thời gian phân tích: Tối ưu hóa từ $28.800\text{ s}$ xuống còn $300 - 600\text{ s}$ trong điều kiện công nghiệp bằng cách sử dụng nguồn hoạt độ $10 - 20\ \mu Ci$ đi kèm che chắn chuẩn.
Hạn chế và hướng phát triển
Hạn chế kỹ thuật
- Sai số hình học ống trụ: Việc sử dụng ống đong hình trụ khiến bề dày xuyên quang thực tế phụ thuộc nghiêm ngặt vào độ chính xác căn trục laser; độ lệch tâm nhỏ có thể gây biến thiên chiều dày lớp chất lỏng $x$, dẫn đến độ lệch $10,17%$ ở axit sulfuric.
- Độ phân giải năng lượng của đầu dò: Đầu dò $NaI(Tl)$ có độ phân giải năng lượng trung bình ($\sim 7%$ tại $662\text{ keV}$), thấp hơn đáng kể so với đầu dò bán dẫn siêu tinh khiết $HPGe$.
- Mô hình năng lượng đơn lẻ: Khóa luận mới khảo sát tại mốc năng lượng $661,7\text{ keV}$ (vùng tán xạ Compton chiếm ưu thế tuyệt đối), chưa mở rộng sang vùng hiệu ứng quang điện ($E < 100\text{ keV}$).
Hướng nâng cấp đề xuất
- Tích hợp nguồn bức xạ đa năng lượng: Kết hợp nguồn $^{241}Am$ ($59,5\text{ keV}$), $^{133}Ba$ ($81, 356\text{ keV}$) và $^{60}Co$ ($1173, 1332\text{ keV}$) để vẽ đồ thị biến thiên $Z_{eff}(E)$ liên tục theo năng lượng.
- Chuyển đổi hình học buồng đo: Sử dụng cuvet phẳng (planar container) có bề dày chuẩn cố định thay cho ống trụ tròn để loại bỏ triệt để sai số đường truyền tia.
- Tự động hóa xử lý tín hiệu: Phát triển module Python nhúng giải thuật Levenberg-Marquardt kết nối trực tiếp cổng USB của MCA để xuất mật độ $\rho$ và $Z_{eff}$ theo thời gian thực (Real-time dashboard).
Đối tượng hưởng lợi
- Sinh viên ngành Vật lý hạt nhân / Kỹ thuật hạt nhân: Nắm vững phương pháp dựng file input MCNP6, xử lý phổ thực nghiệm qua Colegram, Origin và quy tắc truyền sai số thống kê.
- Kỹ sư NDT và tự động hóa công nghiệp: Mô hình tham chiếu để thiết kế hệ đo mật độ chất lỏng, kiểm chuẩn vật liệu che chắn bức xạ không tiếp xúc.
- Doanh nghiệp sản xuất hóa chất & Lọc hóa dầu: Giải pháp kiểm soát chất lượng an toàn, tiết kiệm chi phí, loại bỏ nguy cơ tai nạn lao động do rò rỉ hóa chất ăn mòn.
- Nhóm nghiên cứu chuyên sâu về tương tác bức xạ: Bộ dữ liệu tham chiếu thực nghiệm và mô phỏng đáng tin cậy về $\mu_m$ và $Z_{eff}$ cho 24 loại dung dịch axit.
Câu hỏi thường gặp
1. Tại sao giá trị $Z_{eff}$ của cùng một hợp chất lại thay đổi theo năng lượng photon?
Khác với nguyên tử số $Z$ cố định của đơn chất, $Z_{eff}$ phụ thuộc vào tỉ trọng đóng góp của 3 cơ chế tương tác chính: Hiệu ứng quang điện ($\sigma_{photo} \sim Z^5/E^{7/2}$), Tán xạ Compton ($\sigma_{Compton} \sim Z/E$) và Tạo cặp ($\sigma_{pair} \sim Z^2 \ln E$). Ở mỗi dải năng lượng, cơ chế chiếm ưu thế khác nhau dẫn đến trọng số tương tác của từng nguyên tố thành phần thay đổi, làm cho $Z_{eff}$ là một hàm theo năng lượng $E$.
2. Làm thế nào để đảm bảo hệ số tích lũy $B \approx 1$ trong thực nghiệm?
Để loại bỏ bức xạ tán xạ nhiều lần lọt vào đầu dò, hệ đo bắt buộc phải thiết lập hình học chùm tia hẹp thông qua hệ chuẩn trực chì đôi: đường kính lỗ chuẩn trực nguồn chỉ $0,502\text{ cm}$ và đầu dò $1,00\text{ cm}$, kết hợp căn chỉnh thẳng trục tuyệt đối bằng tia laser.
3. Tại sao kết quả $Z_{eff}$ thực nghiệm của Axit Sulfuric ($H_2SO_4$) lại có độ lệch lớn nhất (10,17%)?
Axit Sulfuric có mật độ rất cao ($\rho = 1,839\text{ g/cm}^3$) và chứa nguyên tố lưu huỳnh ($Z = 16$). Sai số hình học khi chùm tia đi lệch tâm ống trụ thủy tinh sẽ bị khuếch đại mạnh nhất ở các môi trường có độ suy giảm lớn, dẫn đến sai số diện tích đỉnh cao hơn so với các axit nhẹ.
4. Phần mềm MCNP6 đóng vai trò gì so với đo đạc thực nghiệm?
MCNP6 cho phép mô phỏng sự vận chuyển hạt trong môi trường lý tưởng không bị nhiễu bởi phông bức xạ môi trường không kiểm soát, đồng thời mở rộng khảo sát được các chất nguy hiểm, đắt tiền hoặc khó tổng hợp trong phòng thí nghiệm (khảo sát thành công 24 loại axit).
5. Yêu cầu an toàn bức xạ khi triển khai hệ đo này trong thực tế là gì?
Nguồn $^{137}Cs$ sử dụng có hoạt độ rất thấp ($0,76\ \mu Ci$), được bọc trong khối chì bảo vệ $5\text{ cm}$, đảm bảo suất liều bức xạ tại vị trí thao tác tương đương phông tự nhiên ($< 0,2\ \mu Sv/h$), an toàn tuyệt đối cho người vận hành theo tiêu chuẩn TCVN.
Kết luận
Đồ án tốt nghiệp đã giải quyết thành công bài toán xác định nguyên tử số hiệu dụng $Z_{eff}$ và mật độ $\rho$ của các dung dịch axit bằng phương pháp gamma truyền qua. Sự kết hợp giữa thực nghiệm phổ kế gamma $NaI(Tl)$, mô phỏng Monte Carlo MCNP6 và đối sánh chuẩn NIST đã chứng minh tính khả thi, độ chính xác cao (độ lệch thực nghiệm $HNO_3$ chỉ $0,07%$, $H_2O$ chỉ $0,57%$) và độ an toàn vượt trội của kỹ thuật kiểm tra không phá hủy. Kết quả nghiên cứu mở ra tiềm năng ứng dụng to lớn trong tự động hóa công nghiệp hóa chất, dầu khí và thiết kế hệ thống che chắn an toàn bức xạ hiện đại.
