Bối cảnh và vấn đề nghiên cứu

Kỹ thuật kiểm tra không phá hủy (Non-Destructive Testing - NDT) đóng vai trò quan trọng trong việc đánh giá chất lượng máy móc, cấu trúc vật liệu và sản phẩm công nghiệp mà không làm suy giảm hay biến đổi khả năng sử dụng của đối tượng khảo sát. Các phương pháp NDT truyền thống đều tồn tại những ưu điểm và giới hạn kỹ thuật riêng:

  • Phương pháp siêu âm có chi phí thấp và dễ triển khai nhưng chỉ kiểm tra được bề mặt hoặc độ sâu nông, không xuyên sâu để xác định chính xác các khuyết tật ẩn bên trong vật liệu.
  • Phương pháp chụp ảnh phóng xạ có độ chính xác cao nhưng mang lại rủi ro an toàn bức xạ lớn và hình ảnh trên phim không thể hiện trực tiếp thông tin về chiều sâu khuyết tật.
  • Phương pháp thẩm thấu chất lỏng phát hiện tốt vết nứt trên bề mặt kim loại phi từ tính nhưng đòi hỏi bề mặt phải sạch và khô, không phù hợp với các mẫu bị bám bẩn hoặc có độ nhám cao.
  • Phương pháp gamma tán xạ tuy có tiềm năng khảo sát nhưng phức tạp về cơ sở lý thuyết và cách bố trí hệ đo.

Trong bối cảnh đó, phương pháp gamma truyền qua nổi lên như một công cụ hiệu quả để xác định các thông số tương tác của photon với vật chất, bao gồm hệ số suy giảm, tiết diện tương tác, mật độ electron và nguyên tử số hiệu dụng ($Z_{eff}$). Đối với các đơn chất (nhôm, đồng, chì), nguyên tử số $Z$ là một đại lượng cố định. Tuy nhiên, đối với các hợp chất hoặc hỗn hợp phức tạp (dung dịch lỏng, polyme, hợp kim, mô sinh học), một nguyên tử số đơn lẻ không thể đại diện cho toàn bộ vật liệu. Theo nghiên cứu của Hine, $Z_{eff}$ không phải là hằng số cố định mà là một tham số phụ thuộc vào năng lượng photon và các quá trình tương tác tương ứng.

Bên cạnh $Z_{eff}$, mật độ khối $\rho$ là một tham số vật lý then chốt của chất lỏng. Phương pháp đo mật độ truyền thống dựa trên tỉ số giữa khối lượng và thể tích tuy chính xác nhưng tiềm ẩn nguy cơ mất an toàn nghiêm trọng đối với các hóa chất độc hại, dễ bay hơi như axit mạnh hay benzene. Việc áp dụng kỹ thuật gamma truyền qua cho phép xác định mật độ chất lỏng trực tiếp bên trong các bình chứa đóng kín mà không cần tiếp xúc vật lý.

Đề tài đặt ra 3 nhiệm vụ nghiên cứu trọng tâm:

  1. Tính toán các giá trị nguyên tử số hiệu dụng $Z_{eff}$ cho các loại axit bằng phương pháp trực tiếp và phương pháp nội suy.
  2. Thực hiện đo đạc thực nghiệm và mô phỏng Monte Carlo bằng phần mềm MCNP6, sau đó so sánh và đánh giá kết quả thu được với cơ sở dữ liệu chuẩn của Viện Tiêu chuẩn và Kỹ thuật quốc gia Hoa Kỳ (NIST).
  3. Xây dựng đường chuẩn mật độ đối với chất lỏng axit dựa trên tỉ số số đếm thực nghiệm và mô phỏng.

Đối tượng nghiên cứu của khóa luận là các dung dịch axit hữu cơ và vô cơ lỏng (Acetic, Formic, Lactic, Nitric, Phosphoric, Sunfuric, Perchloric, v.v.) cùng mẫu chuẩn nước cất. Đề tài được thực hiện bởi sinh viên Nguyễn Thu Hằng dưới sự hướng dẫn của TS. Hoàng Đức Tâm tại Bộ môn Vật lý hạt nhân, Khoa Vật lý, Trường Đại học Sư phạm TP.HCM (hoàn thành vào tháng 7 năm 2020).

Cơ sở lý thuyết và phương pháp

Khung lý thuyết tương tác bức xạ gamma với vật chất

Bức xạ gamma khi truyền qua vật chất sẽ tương tác với lớp vỏ electron hoặc hạt nhân nguyên tử thông qua ba hiệu ứng chính:

  • Hiệu ứng quang điện: Photon truyền toàn bộ năng lượng cho electron liên kết ở lớp trong, làm electron bứt ra khỏi nguyên tử và phát ra tia X đặc trưng. Tiết diện vi phân quang điện tỉ lệ với: $$\sigma_{photo} \sim \frac{Z^5}{E^{7/2}}$$
  • Tán xạ Compton: Photon truyền một phần năng lượng cho electron lớp ngoài, làm electron bật ra mang động năng và photon bị tán xạ lệch hướng với năng lượng giảm. Tiết diện tán xạ Compton tỉ lệ với: $$\sigma_{Compton} \sim \frac{Z}{E}$$
  • Hiệu ứng tạo cặp: Photon có năng lượng lớn hơn $1022\text{ keV}$ tương tác với trường hạt nhân để tạo thành một cặp electron - positron, sau đó quá trình hủy hạt tạo ra hai photon $511\text{ keV}$. Tiết diện tạo cặp tỉ lệ với: $$\sigma_{pair} \sim Z^2 \ln(E)$$

Tiết diện tương tác tổng cộng được tính bằng: $$\sigma_{total} = \sigma_{photo} + \sigma_{Compt} + \sigma_{pair}$$

Định luật suy giảm chùm tia gamma

Khi chùm tia gamma hẹp đã được chuẩn trực truyền qua lớp vật chất có bề dày $x\text{ (cm)}$ và mật độ khối $\rho\text{ (g/cm}^3)$, cường độ chùm photon suy giảm theo định luật Beer - Lambert: $$I = I_0 e^{-\mu x} = I_0 e^{-\mu_m \rho x}$$ Trong đó:

  • $I_0$ và $I$ lần lượt là cường độ chùm tia trước và sau khi đi qua mẫu.
  • $\mu$ là hệ số suy giảm tuyến tính ($\text{cm}^{-1}$), được xác định bởi: $\mu = \frac{1}{x} \ln\left(\frac{I_0}{I}\right)$.
  • $\mu_m$ là hệ số suy giảm khối ($\text{cm}^2\text{/g}$), liên hệ qua công thức: $\mu_m = \frac{\mu}{\rho}$.
  • Đối với hợp chất, $\mu_m$ tuân theo quy tắc cộng tính: $\mu_m = \sum w_i (\mu_m)_i$ với $w_i$ là phần trăm khối lượng của nguyên tố thứ $i$.

Phương pháp nghiên cứu và công cụ thực nghiệm

  1. Phương pháp thực nghiệm:
    • Hệ đo sử dụng đầu dò nhấp nháy NaI(Tl) kết nối bộ xử lý tín hiệu và phần mềm thu nhận phổ ADMCA.
    • Nguồn phóng xạ gamma đơn năng $^{137}\text{Cs}$ phát bức xạ năng lượng $662\text{ keV}$.
    • Mẫu chất lỏng được chứa trong các ống trụ có kích thước hình học xác định.
    • Phổ thực nghiệm được phân tích vùng đỉnh hấp thụ toàn phần (FEA) bằng phần mềm Colegram (khớp đỉnh dạng Gauss kết hợp đường nền One Step) và Origin.
  2. Phương pháp mô phỏng Monte Carlo:
    • Sử dụng mã tính MCNP6 (Los Alamos National Laboratory) để mô phỏng sự vận chuyển của photon qua hình học hệ đo và các loại vật liệu axit.
  3. Phương pháp xử lý sai số:
    • Sử dụng giải thuật truyền sai số cho các hàm nhiều biến. Phương sai của tỉ số số đếm $R = \frac{N_{Acid}}{N_{H_2O}}$ được tính với hệ số $2,58$, tương ứng độ tin cậy $99,99%$.

Nội dung chính theo từng chương

Chương 1. Tổng quan về gamma truyền qua

Chương này hệ thống hóa lịch sử phát hiện tia gamma từ công trình của Paul Villard năm 1900 đến các định danh và mô tả bản chất của Ernest Rutherford năm 1903. Tác giả phân tích chi tiết cơ chế vi mô của ba quá trình tương tác (quang điện, tán xạ Compton, tạo cặp) và phân bố năng lượng trên phổ gamma thực tế với các thành phần: đỉnh hấp thụ toàn phần (FEA), tán xạ đơn (SS), tán xạ nhiều lần (MS), tán xạ ngược (BS), đỉnh thoát đơn (SPE) và đỉnh thoát đôi (DPE). Chương 1 cũng trình bày định luật suy giảm bức xạ cho chùm tia hẹp, chùm tia rộng (kèm hệ số tích lũy $B$) và mối liên hệ xác định mật độ khối $\rho$.

Chương 2. Phương pháp xác định nguyên tử số hiệu dụng

Chương 2 trình bày tổng quan 4 phương pháp xác định $Z_{eff}$ trong các nghiên cứu thực nghiệm và lý thuyết:

  1. Phương pháp tính trực tiếp: Xác định $Z_{eff}$ từ công thức giải tích dựa trên thành phần khối lượng nguyên tố và hệ số suy giảm khối.
  2. Phương pháp tính theo tiết diện nguyên tử: Tính toán $Z_{eff}$ thông qua mối tương quan giữa tiết diện nguyên tử tổng và tiết diện điện tử.
  3. Phương pháp nội suy: Xây dựng đường cong chuẩn tiết diện tương tác của các nguyên tố tinh khiết theo nguyên tử số $Z$, sau đó nội suy giá trị $Z_{eff}$ của mẫu hỗn hợp.
  4. Phương pháp sử dụng công thức XMuDat: Áp dụng phần mềm và cơ sở dữ liệu chuyên dụng để tính toán nhanh các tham số tương tác.

Chương 3. Khảo sát nguyên tử số hiệu dụng và mật độ của axit bằng mô phỏng Monte Carlo và thực nghiệm

Nội dung chương 3 tập trung mô tả chi tiết quy trình triển khai:

  • Thiết lập hệ đo thực nghiệm gamma truyền qua gồm nguồn $^{137}\text{Cs}$, hệ chuẩn trực chì, giá đỡ ống chứa mẫu axit và detector NaI(Tl).
  • Xây dựng mô hình MCNP6 mô phỏng cấu trúc đầu dò NaI(Tl), kích thước ống nghiệm và các thẻ vật liệu chứa thành phần hóa học của axit.
  • Quy trình xử lý phổ thực nghiệm qua phần mềm Colegram nhằm bóc tách diện tích đỉnh thực nghiệm $N$ và trừ phông nền môi trường trên Origin.
  • Thiết lập hệ thống công thức đánh giá sai số thực nghiệm và thuật toán truyền sai số cho tỉ số số đếm $R$.

Chương 4. Xác định nguyên tử số hiệu dụng và mật độ của một số axit

Chương 4 công bố toàn bộ dữ liệu đo đạc thực nghiệm trên 6 loại axit và dữ liệu mô phỏng trên 24 loại axit. Tác giả đối chiếu hệ số suy giảm khối $\mu_m$, tính toán giá trị $Z_{eff}$, kiểm tra tiêu chuẩn hợp lệ $Z_{min} < Z_{eff} < Z_{max}$ và phân tích độ sai biệt tương đối (RD%) so với NIST. Đồng thời, chương 4 thiết lập đường chuẩn tương quan tuyến tính giữa mật độ $\rho$ và $\ln(R)$, chứng minh khả năng đo mật độ chất lỏng không tiếp xúc.

Kết quả và đóng góp

Kết quả đo hệ số suy giảm khối và nguyên tử số hiệu dụng thực nghiệm

Số liệu đo diện tích đỉnh thực nghiệm $N_{TN}$, hệ số suy giảm khối $\mu_m$ và $Z_{eff}$ của 6 loại axit và nước cất khi chiếu xạ bởi nguồn $^{137}\text{Cs}$ ($662\text{ keV}$) được tổng hợp dưới đây:

STT Loại mẫu Diện tích đỉnh $N_{TN}$ (số đếm) $\mu_m$ Thực nghiệm ($\text{cm}^2\text{/g}$) $\mu_m$ NIST ($\text{cm}^2\text{/g}$) Độ sai biệt $\mu_m$ (%) $Z_{eff}$ Thực nghiệm $Z_{eff}$ NIST Độ sai biệt $Z_{eff}$ (%)
1 Axit Lactic $70669,83 \pm 685,86$ $0,079 \pm 0,005$ 0,082 4,14 $4,13 \pm 0,09$ 4,00 3,24
2 Axit Acetic $71883,80 \pm 691,73$ $0,085 \pm 0,006$ 0,082 3,23 $3,83 \pm 0,07$ 4,00 4,28
3 Axit Formic $70488,80 \pm 684,98$ $0,079 \pm 0,005$ 0,080 2,18 $4,70 \pm 0,07$ 4,80 2,02
4 Axit Nitric $67690,20 \pm 671,25$ $0,080 \pm 0,005$ 0,078 2,05 $6,42 \pm 0,07$ 6,43 0,07
5 Axit Photphoric $61533,98 \pm 640,00$ $0,087 \pm 0,004$ 0,079 10,09 $6,04 \pm 0,05$ 6,19 2,44
6 Axit Sunfuric $62405,97 \pm 644,51$ $0,077 \pm 0,004$ 0,079 2,20 $7,87 \pm 0,06$ 7,14 10,17
7 Nước cất $72571,39 \pm 695,03$ $0,085 \pm 0,006$ 0,086 0,58 $3,31 \pm 0,09$ 3,33 0,57

Tất cả giá trị $Z_{eff}$ thực nghiệm đều thỏa mãn điều kiện $Z_{min} < Z_{eff} < Z_{max}$. Ví dụ đối với $\text{H}2\text{SO}4$ cấu thành từ H ($Z=1$), O ($Z=8$), S ($Z=16$), giá trị $Z{eff} = 7,87$ nằm hoàn toàn trong khoảng từ 1 đến 16. Mẫu Axit Nitric có độ sai biệt $Z{eff}$ nhỏ nhất ($0,07%$), trong khi Axit Sunfuric có độ sai biệt lớn nhất ($10,17%$).

So sánh nguyên tử số hiệu dụng giữa mô phỏng MCNP6, thực nghiệm và NIST

Hệ số suy giảm khối lượng của các nguyên tố đơn chất dùng trong tính toán gồm: H ($0,1531$), C ($0,0771$), N ($0,0771$), O ($0,0771$), S ($0,0770$), Cl ($0,0740$) và P ($0,0747\text{ cm}^2\text{/g}$).

Bảng so sánh kết quả $Z_{eff}$ giữa ba phương pháp:

Loại axit $Z_{eff}$ NIST $Z_{eff}$ Thực nghiệm $Z_{eff}$ Mô phỏng Độ sai biệt RD (%) giữa Mô phỏng và Thực nghiệm
Axit Lactic 4,00 4,13 3,74 9,44
Axit Acetic 4,00 3,83 3,75 2,09
Axit Formic 4,80 4,70 4,49 4,47
Axit Nitric 6,43 6,42 5,87 8,57
Axit Photphoric 6,19 6,04 5,76 4,64
Axit Sunfuric 7,14 7,87 6,65 15,50

Kết quả mô phỏng trên 24 loại axit cho thấy độ lệch tương đối của hệ số suy giảm khối $\mu_m$ so với NIST đều dưới $8%$. Dữ liệu $Z_{eff}$ mô phỏng cũng phù hợp tốt với các công bố trước đây tại năng lượng $1\text{ MeV}$ (như Arachidic: 2,84; Behenic: 2,82; Cerotic: 2,80; Lauric: 2,95; Oleic: 2,93; Palmitic: 2,88; Stearic: 2,86).

Xây dựng đường chuẩn xác định mật độ chất lỏng

Tỉ số số đếm $R = \frac{N_{Acid}}{N_{H_2O}}$ được khảo sát trên 21 axit mô phỏng và mẫu nước cất để thiết lập tương quan với mật độ $\rho$:

STT Loại axit Mật độ $\rho\text{ (g/cm}^3)$ Tỉ số $R$ $\ln(R)$
1 Oleic 0,781 1,05 0,04
2 Montanic 0,819 1,04 0,04
3 Cerotic 0,820 1,04 0,04
4 Lignoceric 0,821 1,04 0,04
5 Behenic 0,822 1,04 0,03
6 Arachidic 0,824 1,04 0,03
7 Palmitic 0,852 1,03 0,03
8 Palmitoleic 0,894 1,02 0,02
9 Myristic 0,990 1,00 0,00
10 Lauric 1,007 1,00 0,00
11 Acetic 1,045 1,00 0,00
12 Acrylic 1,051 1,00 0,00
13 Leucine 1,170 0,97 -0,03
14 Lactic 1,206 0,97 -0,03
15 Formic 1,220 0,97 -0,03
16 Aminobutanoic 1,230 0,96 -0,04
17 Adipic 1,360 0,94 -0,06
18 Nitric 1,513 0,92 -0,08
19 Photphoric 1,685 0,89 -0,12
20 Perchloric 1,770 0,88 -0,13
21 Sunfuric 1,830 0,87 -0,14
22 Nước cất 1,000 1,00 0,00

Vì tỉ số mô phỏng $R_{MP}$ và thực nghiệm $R_{TN}$ có độ tương thích cao, đường chuẩn $\ln(R)$ theo mật độ $\rho$ cho phép nội suy trực tiếp mật độ của các mẫu axit thực nghiệm với độ tin cậy cao.

Hạn chế và hướng nghiên cứu tiếp

  • Hạn chế thực nghiệm: Việc sử dụng các ống chứa mẫu hình trụ tròn làm tăng nguy cơ sai lệch hình học. Rất khó để đảm bảo chùm tia gamma đã chuẩn trực truyền qua chính xác tâm đối xứng của ống trụ và đi thẳng vào tâm bề mặt đầu dò mà không bị lệch góc. Đây là nguyên nhân chính dẫn đến độ sai biệt lớn ở một số mẫu (Axit Sunfuric chênh lệch $10,17%$ ở thực nghiệm và $15,50%$ so với mô phỏng).
  • Hạn chế mô phỏng: Mô hình hình học trong MCNP6 cần tiếp tục được tinh chỉnh chi tiết hơn để giảm thiểu sự chênh lệch có tính hệ thống đối với thực nghiệm.
  • Hướng mở rộng: Cải tiến hệ chuẩn trực và cấu trúc giá đỡ ống chứa mẫu, tối ưu hóa các thông số hình học trong tệp đầu vào MCNP6 và mở rộng dải khảo sát sang các dung dịch hóa chất công nghiệp khác.

Giá trị tham khảo

Khóa luận là tài liệu tham khảo chuyên khảo hữu ích cho sinh viên, học viên cao học và nghiên cứu sinh thuộc các ngành Vật lý học, Vật lý hạt nhân, Kỹ thuật hạt nhân và Kỹ thuật an toàn bức xạ.

Các nội dung có giá trị ứng dụng cao gồm:

  • Quy trình chuẩn hóa việc phân tích phổ kế gamma thực nghiệm bằng phần mềm Colegram và Origin.
  • Phương pháp xây dựng tệp dữ liệu đầu vào mô phỏng tương tác bức xạ gamma trên MCNP6.
  • Bộ số liệu thực nghiệm và mô phỏng về $\mu_m$, $Z_{eff}$ và tỉ số $R$ của 24 loại axit.
  • Phương pháp luận đo mật độ chất lỏng nguy hại trong bình kín bằng kỹ thuật NDT hạt nhân.

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

1. Tại sao nguyên tử số hiệu dụng $Z_{eff}$ không phải là một hằng số cố định của vật liệu?

Theo Hine (1952), đối với các hợp chất hoặc hỗn hợp, bức xạ gamma tương tác đồng thời với nhiều loại nguyên tử khác nhau. Do xác suất xảy ra của từng hiệu ứng (quang điện tỉ lệ $Z^5/E^{7/2}$, Compton tỉ lệ $Z/E$, tạo cặp tỉ lệ $Z^2\ln(E)$) phụ thuộc vào mức năng lượng của photon tới, giá trị $Z_{eff}$ sẽ thay đổi tương ứng theo năng lượng của chùm bức xạ.

2. Mẫu axit nào có độ sai biệt $Z_{eff}$ thực nghiệm thấp nhất và cao nhất so với dữ liệu chuẩn NIST?

Trong các phép đo thực nghiệm với nguồn $^{137}\text{Cs}$, Axit Nitric đạt độ sai biệt thấp nhất so với NIST là $0,07%$ ($Z_{eff}$ thực nghiệm đạt $6,42 \pm 0,07$ so với NIST là $6,43$). Axit Sunfuric có độ sai biệt lớn nhất là $10,17%$ ($Z_{eff}$ thực nghiệm đạt $7,87 \pm 0,06$ so với NIST là $7,14$).

3. Phương pháp gamma truyền qua khắc phục nhược điểm gì của phương pháp đo mật độ truyền thống?

Phương pháp truyền thống (cân khối lượng và đo thể tích) gây nguy hiểm trực tiếp cho người thực hiện khi tiếp xúc với các hóa chất độc hại, ăn mòn mạnh hoặc dễ bay hơi như axit đặc, benzene. Phương pháp gamma truyền qua cho phép xác định mật độ gián tiếp khi dung dịch vẫn được niêm phong an toàn trong bình chứa kín.

4. Tỉ số $R$ dùng để dựng đường chuẩn mật độ được tính toán như thế nào?

Tỉ số $R$ là tỉ số giữa diện tích đỉnh hấp thụ toàn phần của mẫu axit ($N_{Acid}$) và diện tích đỉnh của mẫu nước cất ($N_{H_2O}$) đo trong cùng điều kiện hình học: $$R = \frac{N_{Acid}}{N_{H_2O}}$$ Khi mật độ chất lỏng tăng từ $0,781\text{ g/cm}^3$ (Oleic) đến $1,830\text{ g/cm}^3$ (Sunfuric), tỉ số $R$ giảm từ $1,05$ xuống $0,87$, tạo nên đường chuẩn tuyến tính theo $\ln(R)$.

5. Tiêu chuẩn nào được sử dụng để đánh giá tính hợp lệ của giá trị $Z_{eff}$ tính toán được?

Giá trị nguyên tử số hiệu dụng của một hợp chất được xem là hợp lệ khi nằm trong khoảng giới hạn giữa nguyên tử số nhỏ nhất và lớn nhất của các nguyên tố cấu thành: $Z_{min} < Z_{eff} < Z_{max}$. Toàn bộ các giá trị $Z_{eff}$ thực nghiệm và mô phỏng trong nghiên cứu đều thỏa mãn tiêu chuẩn này.

Kết luận

Khóa luận tốt nghiệp của tác giả Nguyễn Thu Hằng đã nghiên cứu và ứng dụng thành công phương pháp gamma truyền qua kết hợp mô phỏng Monte Carlo (MCNP6) để 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. Kết quả thực nghiệm và mô phỏng trên nguồn $^{137}\text{Cs}$ ($662\text{ keV}$) cho thấy sự phù hợp cao với cơ sở dữ liệu chuẩn NIST và các công bố quốc tế. Nghiên cứu đã chứng minh tính khả thi của kỹ thuật kiểm tra không phá hủy bằng bức xạ hạt nhân trong việc phân tích tham số vật lý của các chất lỏng độc hại trong bình chứa kín.