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

Xác định mật độ chất lỏng là yêu cầu kỹ thuật nền tảng trong nhiều ngành công nghiệp như kiểm soát chất lượng xăng dầu, nhiên liệu sinh học theo tiêu chuẩn ASTM, giám sát nồng độ cồn và hàm lượng đường trong sản xuất thực phẩm - đồ uống, cũng như định lượng dung tích hóa chất và dược phẩm. Các phương pháp đo tỷ trọng truyền thống (tỷ trọng kế, bình đo tỷ trọng pycnometer, cân thủy tĩnh dựa trên lực đẩy Archimedes) tồn tại nhiều điểm hạn chế: đòi hỏi thao tác thủ công, yêu cầu lượng mẫu lớn, phụ thuộc vào kỹ năng người đo, dễ phát sinh sai số do hiện tượng dính ướt, bọt khí hoặc độ nhớt cao. Quan trọng hơn, việc phải tiếp xúc trực tiếp với mẫu thử gây nguy cơ mất an toàn nghiêm trọng khi đối tượng đo là các acid mạnh, chất ăn mòn, chất độc hại hoặc chất dễ bay hơi.

Để khắc phục nhược điểm trên, kỹ thuật kiểm tra không phá hủy (NDT) ứng dụng bức xạ hạt nhân được xem là giải pháp thay thế tối ưu nhờ khả năng đo gián tiếp xuyên qua thành ống chứa mà không làm biến đổi mẫu. Tuy nhiên, các kỹ thuật đo bức xạ hiện hành vẫn tồn tại những rào cản kỹ thuật nhất định:

  • Phương pháp gamma tán xạ ngược Compton đòi hỏi cấu hình hình học phức tạp, hiệu ứng tán xạ nhiều lần gây nhiễu phổ lớn và cần thời gian đo kéo dài để đảm bảo độ tin cậy thống kê nếu dùng nguồn hoạt độ thấp.
  • Phương pháp gamma truyền qua có thời gian đo nhanh hơn và bố trí đơn giản hơn, nhưng việc lập đường chuẩn thực nghiệm trực tiếp thường bị ảnh hưởng lớn bởi điều kiện môi trường và độ trôi thiết bị (như nghiên cứu của Priyada và cộng sự cho thấy độ lệch ở mẫu hexane lên tới 21,54%).
  • Phương pháp bán thực nghiệm của Chuong và cộng sự đã cải thiện độ chính xác bằng cách xây dựng đường chuẩn từ mô phỏng Monte Carlo dựa trên tỷ số cường độ tương đối $R$, nhưng lại gặp hạn chế thực tiễn: mỗi khi thay đổi đường kính ống chứa mẫu $d$, người phân tích phải thiết lập lại một đường chuẩn mô phỏng mới, gây tốn kém thời gian tính toán và làm mất tính linh hoạt khi triển khai công nghiệp.

Khóa luận tốt nghiệp của tác giả Nguyễn Thu Hằng đặt ra mục tiêu và các nhiệm vụ nghiên cứu sau:

  1. Nghiên cứu cơ sở lý thuyết về sự suy giảm chùm tia gamma truyền qua vật chất và phương pháp bán thực nghiệm xác định mật độ chất lỏng.
  2. Xây dựng mô hình mô phỏng Monte Carlo bằng chương trình MCNP6 để thiết lập hệ thống 20 đường chuẩn tỷ số cường độ bức xạ theo mật độ tương ứng với 20 kích thước đường kính ống chứa (từ 1,0 cm đến 10,0 cm).
  3. Thiết kế, huấn luyện và tối ưu hóa mô hình mạng nơ-ron nhân tạo (ANN) dạng Perceptron đa lớp (MLP) nhằm dự đoán tự động mật độ chất lỏng với hai thông số đầu vào là tỷ số cường độ $R$ và đường kính ống $d$.
  4. Bố trí hệ đo thực nghiệm sử dụng đầu dò nhấp nháy NaI(Tl) và nguồn phóng xạ $^{137}\text{Cs}$ cường độ thấp để đo phổ truyền qua của 7 mẫu chất lỏng, từ đó kiểm chứng độ chính xác và tính khả thi của mô hình mạng nơ-ron đã phát triển.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Sự suy giảm bức xạ gamma năng lượng 661,7 keV truyền qua các chất lỏng có dải mật độ từ 0,78 $\text{g/cm}^3$ đến 1,84 $\text{g/cm}^3$; mô hình mạng nơ-ron nhân tạo MLP.
  • Phạm vi thực nghiệm: Hệ phổ kế gamma tại Phòng thí nghiệm Bộ môn Vật lý hạt nhân, Khoa Vật lý, Trường Đại học Sư phạm TP.HCM. Cấu hình hình học thực nghiệm cố định trên ống nghiệm thủy tinh đường kính trong 2,68 cm (đường kính ngoài 3,05 cm); nguồn $^{137}\text{Cs}$ có hoạt độ tại thời điểm đo là $(0,76 \pm 0,03)\ \mu\text{Ci}$.
  • Phạm vi mô phỏng và tính toán: Dữ liệu mô phỏng Monte Carlo trên 16 hợp chất lỏng với 20 kích thước đường kính ống hình trụ từ 1,0 cm đến 10,0 cm (bước nhảy 0,5 cm). Thời gian hoàn thành khóa luận: tháng 5 năm 2021.

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

Khóa luận xây dựng trên nền tảng lý thuyết tương tác giữa bức xạ gamma với vật chất, bao gồm ba hiệu ứng chính:

  1. Hiệu ứng quang điện: Photon truyền toàn bộ năng lượng cho electron liên kết ở lớp vỏ trong nguyên tử (lớp K, L). Tiết diện tương tác $\sigma_{\text{photo}}$ tỷ lệ thuận với $Z^5$ và tỷ lệ nghịch với $E^{7/2}$.
  2. Hiệu ứng tán xạ Compton: Va chạm giữa photon và electron liên kết yếu ở lớp vỏ ngoài, làm tán xạ photon và truyền một phần động năng cho electron. Tiết diện tán xạ $\sigma_{\text{Compt}}$ tỷ lệ thuận với số nguyên tử $Z$ và tỷ lệ nghịch với năng lượng photon $E$.
  3. Hiệu ứng tạo cặp electron - positron: Chỉ xuất hiện khi năng lượng photon vượt ngưỡng $2m_0c^2 \approx 1,022\text{ MeV}$. Với nguồn $^{137}\text{Cs}$ ($E = 661,7\text{ keV}$), hiệu ứng này không xảy ra.

Sự suy giảm cường độ của chùm tia gamma đơn năng, hẹp truyền qua lớp vật chất có bề dày $x$ (cm) và mật độ $\rho$ ($\text{g/cm}^3$) tuân theo định luật Beer - Lambert: $$I = I_0 e^{-\mu x} = I_0 e^{-\mu_m \rho x}$$ Trong đó $I_0$ là cường độ ban đầu, $I$ là cường độ sau khi truyền qua, $\mu$ ($\text{cm}^{-1}$) là hệ số suy giảm tuyến tính, $\mu_m$ ($\text{cm}^2/\text{g}$) là hệ số suy giảm khối. Khi chùm tia không đạt điều kiện hẹp lý tưởng, phương trình được hiệu chỉnh bằng hệ số tích lũy chùm tia $B(E, Z)$: $$I = I_0 B(E, Z) e^{-\mu_m \rho x}$$

Nhằm loại bỏ sự phụ thuộc vào hiệu suất tuyệt đối của đầu dò cũng như các yếu tố nhiễu môi trường giữa thực nghiệm và mô phỏng, tác giả kế thừa cách tiếp cận sử dụng tỷ số cường độ tương đối $R$: $$R_{\text{exp}} = \frac{N_{\text{exp}}(\text{chất lỏng})}{N_{\text{exp}}(\text{nước})}, \quad R_{\text{sim}} = \frac{N_{\text{sim}}(\text{chất lỏng})}{N_{\text{sim}}(\text{nước})}$$ Với $N$ là diện tích đỉnh hấp thụ toàn phần năng lượng 661,7 keV. Do cùng một hình học bố trí, hiệu suất ghi của đầu dò triệt tiêu trong phép chia, dẫn đến $R_{\text{exp}} \approx R_{\text{sim}}$. Mối quan hệ giữa tỷ số cường độ $R$ và mật độ $\rho$ được tuyến tính hóa qua phương trình: $$\ln(R) = a + b\rho$$

Phương pháp nghiên cứu kết hợp giữa mô phỏng tính toán và thực nghiệm kiểm chứng:

  • Mô phỏng Monte Carlo bằng MCNP6: Mô hình hóa cấu hình hình học thực nghiệm với 10 tỷ hạt phát xạ ($10^{10}$ histories) để đảm bảo sai số thống kê nhỏ cho phổ phân bố độ cao xung (PHS) trên 8192 kênh. Tham số làm rộng phổ FWHM được đưa vào mã nguồn thông qua hàm thực nghiệm: $$\text{FWHM}(\text{MeV}) = a + b\sqrt{E} + cE^2$$ với các hệ số $a = -0,00526\text{ MeV}$, $b = 0,06096\text{ MeV}^{1/2}$, $c = 0,01401\text{ MeV}^{-1}$.
  • Xử lý phổ thực nghiệm và mô phỏng bằng Colegram: Sử dụng thuật toán Levenberg - Marquardt kết hợp phương pháp bình phương tối thiểu để khớp đỉnh hấp thụ toàn phần bằng hàm Gaussian và hàm phông bậc thang One Step. Vùng quan tâm (ROI) được chọn nghiêm ngặt trong khoảng $[\text{vị trí đỉnh} - 2,5\text{FWHM}; \text{vị trí đỉnh} + 2,5\text{FWHM}]$ để khắc phục hiện tượng dịch đỉnh (peak drift) do nhiệt độ.
  • Xây dựng mô hình trí tuệ nhân tạo (ANN): Sử dụng cấu trúc mạng Perceptron đa lớp (MLP) được lập trình và huấn luyện trên nền tảng phần mềm Wolfram Mathematica.

Nguồn dữ liệu của nghiên cứu:

  • Dữ liệu mô phỏng: 304 điểm dữ liệu sinh ra từ MCNP6, gồm 16 chất lỏng hữu cơ và acid có hệ số suy giảm khối biến thiên đồng nhất (dao động trong khoảng 11%) trên 20 giá trị đường kính ống.
  • Dữ liệu thực nghiệm: Đo đạc trực tiếp phổ bức xạ của 7 mẫu chất lỏng chuẩn (nước cất, acetic, formic, lactic, nitric, phosphoric, sulfuric) với thời gian đo liên tục 8 giờ cho mỗi mẫu.

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

Chương 1: Tổng quan về phương pháp gamma truyền qua và phương pháp xác định mật độ chất lỏng

Chương 1 hệ thống hóa lý thuyết vật lý hạt nhân về bức xạ gamma và các phương pháp đo mật độ chất lỏng. Tác giả phân tích chi tiết cơ chế tương tác vi mô của photon với lớp vỏ điện tử và hạt nhân nguyên tử, chỉ ra đóng góp của từng thành phần trên phổ năng lượng gamma của đồng vị $^{137}\text{Cs}$ (gồm đỉnh hấp thụ toàn phần 661,7 keV, nền tán xạ Compton, đỉnh tán xạ ngược và gờ Compton).

Nội dung chương so sánh ưu và nhược điểm của các dụng cụ đo tỷ trọng cơ học và điện tử:

  • Bình tỷ trọng pycnometer và tỷ trọng kế thủy tinh: phụ thuộc vào kỹ năng đọc, bọt khí, sự dính ướt thành bình và thay đổi nhiệt độ.
  • Máy đo tỷ trọng hiện số (dao động ống chữ U): độ chính xác cao nhưng chi phí thiết bị đắt đỏ, quy trình bảo dưỡng phức tạp.
  • Phương pháp bức xạ gamma không tiếp xúc: đảm bảo an toàn tuyệt đối khi phân tích các hóa chất ăn mòn mạnh và độc hại.

Tác giả thiết lập công thức toán học dẫn dắt từ định luật Beer - Lambert đến phương trình hàm khớp tuyến tính $\ln(R) = a + b\rho$, làm sáng tỏ tính độc lập của tỷ số $R$ đối với hiệu suất ghi nhận của đầu dò bức xạ.

Chương 2: Khảo sát mật độ chất lỏng bằng phương pháp bán thực nghiệm

Chương 2 trình bày chi tiết cách thiết lập thực nghiệm, thông số kỹ thuật thiết bị và quy trình mô phỏng hạt nhân:

Thông số hình học và thiết bị Giá trị cấu hình thực nghiệm
Nguồn bức xạ $^{137}\text{Cs}$ ($E = 661,7\text{ keV}$, hoạt độ $0,76\ \mu\text{Ci}$)
Khối chì che chắn nguồn $12,39\text{ cm} \times 12,39\text{ cm} \times 5,48\text{ cm}$
Đường kính lỗ chuẩn trực nguồn $0,96\text{ cm}$
Đầu dò bức xạ Đầu dò nhấp nháy $\text{NaI(Tl)}$ hãng AMPTEK
Kích thước tinh thể $\text{NaI(Tl)}$ Đường kính $7,62\text{ cm} \times$ Chiều dài $7,62\text{ cm}$
Mật độ tinh thể $\text{NaI(Tl)}$ $3,670\text{ g/cm}^3$
Đường kính lỗ chuẩn trực đầu dò $1,00\text{ cm}$
Khoảng cách nguồn - tâm ống mẫu $12,4\text{ cm}$
Khoảng cách tâm ống mẫu - đầu dò $12,4\text{ cm}$
Kích thước ống mẫu thủy tinh Đường kính trong $2,68\text{ cm}$; đường kính ngoài $3,05\text{ cm}$
Chế độ phân tích đa kênh MCA 8192 kênh

Danh mục 7 mẫu chất lỏng đo thực nghiệm:

  • Nước cất ($\text{H}_2\text{O}$, $1,000\text{ g/cm}^3$, 100%, Merck)
  • Acetic acid ($\text{C}_2\text{H}_4\text{O}_2$, $1,045\text{ g/cm}^3$, 99%, Xilong)
  • Formic acid ($\text{CH}_2\text{O}_2$, $1,220\text{ g/cm}^3$, 88%, Xilong)
  • Lactic acid ($\text{C}_3\text{H}_6\text{O}_3$, $1,206\text{ g/cm}^3$, 88%, Xilong)
  • Nitric acid ($\text{HNO}_3$, $1,513\text{ g/cm}^3$, 65%, Merck)
  • Phosphoric acid ($\text{H}_3\text{PO}_4$, $1,390\text{ g/cm}^3$, 85%, VWR Chemicals)
  • Sulfuric acid ($\text{H}_2\text{SO}_4$, $1,839\text{ g/cm}^3$, 98%, Merck)

Tác giả phân tích hiện tượng dịch đỉnh phổ do thời gian tích lũy xung kéo dài 8 giờ/mẫu dưới tác động của nhiệt độ môi trường và nhiễu điện tử. Quy trình xử lý bằng phần mềm Colegram được cố định trong dải $\pm 2,5\text{FWHM}$ quanh tâm đỉnh Gauss, kết hợp hàm nền One Step giúp chuẩn hóa diện tích đỉnh hấp thụ toàn phần giữa các lần đo. Trong phần mô phỏng MCNP6, tác giả thiết lập thẻ định nghĩa vật liệu cho 16 hợp chất lỏng chuẩn (mật độ từ 0,79 đến 1,83 $\text{g/cm}^3$) và chạy bài toán truyền hạt với 10 tỷ photon để thu thập phổ truyền qua cho 20 kích thước ống.

Chương 3: Phương pháp xây dựng mô hình mạng nơ-ron nhân tạo để dự đoán mật độ của một số chất lỏng

Chương 3 tập trung phát triển mô hình học máy ANN giải quyết bài toán hồi quy mật độ chất lỏng đa biến. Cấu trúc mạng Multilayer Perceptron (MLP) được thiết lập với các thông số:

  • Lớp vào (Input layer): Gồm 2 nơ-ron nhận giá trị tỷ số cường độ $R$ và đường kính trong của ống chứa $d$ (cm).
  • Lớp ẩn (Hidden layer): Gồm 1 lớp ẩn. Số lượng nơ-ron trong lớp ẩn được khảo sát và lựa chọn thông qua chỉ số lỗi căn bậc hai trung bình (RMSE) và hệ số xác định ($R^2$).
  • Lớp ra (Output layer): Gồm 1 nơ-ron xuất giá trị mật độ dự đoán $\rho_{\text{pre}}$ ($\text{g/cm}^3$).
  • Hàm kích hoạt (Activation function): Sử dụng hàm phi tuyến Softsign: $$\text{Softsign}(x) = \frac{x}{1 + |x|}$$
Đặc trưng cấu trúc mạng Thông số thiết lập
Kiểu kiến trúc Multilayer Perceptron (MLP) truyền thẳng
Số nơ-ron lớp vào 2 (Tỷ số $R$, Đường kính ống $d$)
Số lớp ẩn 1
Số nơ-ron lớp ra 1 (Mật độ chất lỏng $\rho$)
Hàm kích hoạt lớp ẩn Softsign
Công cụ huấn luyện Wolfram Mathematica
Thuật toán tối ưu Lan truyền ngược giảm thiểu sai số bình phương trung bình
Quy mô tập dữ liệu 304 mẫu từ mô phỏng MCNP6
Phân chia dữ liệu 85% tập huấn luyện (training), 15% tập kiểm tra (testing)
Tập dữ liệu đánh giá độc lập Dữ liệu đo đạc thực nghiệm trên hệ phổ kế $\text{NaI(Tl)}$

Sai số của mật độ dự đoán được lan truyền từ sai số của tỷ số thực nghiệm $\Delta R$ theo công thức thống kê với hệ số độ tin cậy 99,99% ($k = 2,58$).

Chương 4: Xác định mật độ của một số chất lỏng

Chương 4 công bố và thảo luận các kết quả tính toán thực nghiệm và mô phỏng:

  1. Xây dựng 20 đường chuẩn mô phỏng: Tác giả xây dựng thành công 20 hàm tuyến tính dạng $\ln(R_{\text{sim}}) = a + b\rho$ tương ứng với 20 đường kính từ 1,0 cm đến 10,0 cm. Khi đường kính ống tăng dần, độ dốc $b$ của đường chuẩn tăng theo, phản ánh mức độ suy giảm bức xạ mạnh hơn khi chiều dày chất lỏng tăng.
  2. Tính toán mật độ theo phương pháp bán thực nghiệm truyền thống (ống 2,68 cm): Thay giá trị $R_{\text{exp}}$ của 7 chất lỏng vào phương trình đường chuẩn mô phỏng tương ứng ($d = 2,68\text{ cm}$), mật độ tính toán cho thấy độ tương hợp cao với mật độ danh định của nhà sản xuất.
  3. Đánh giá mô hình ANN: Quá trình huấn luyện trên phần mềm Mathematica với 85% dữ liệu mô phỏng đạt độ hội tụ cao. Kiểm tra mô hình trên 15% dữ liệu thử nghiệm mô phỏng còn lại cho thấy các chỉ số $R^2$ đều tiệm cận 1 và RMSE ở mức thấp trên toàn bộ 20 kích thước đường kính ống.
  4. Đối sánh kết quả giữa đường chuẩn và mạng ANN: Mật độ dự đoán từ mô hình ANN cho 7 chất lỏng thực nghiệm được đặt cạnh kết quả tính từ đường chuẩn mô phỏng và mật độ chuẩn để phân tích độ lệch tương đối $\text{RD}\ (%)$. Cả hai phương pháp đều kiểm soát sai số tốt, chứng minh mạng ANN hoàn toàn có khả năng thay thế việc lập từng đường chuẩn riêng lẻ.

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

Khóa luận đã đạt được các kết quả định lượng cụ thể:

  • Thiết lập thành công cơ sở dữ liệu gồm 304 điểm dữ liệu mô phỏng Monte Carlo MCNP6, bao phủ dải mật độ từ 0,78 $\text{g/cm}^3$ đến 1,84 $\text{g/cm}^3$ trên 20 cỡ ống chứa từ 1,0 cm đến 10,0 cm.
  • Hệ thống hóa 20 phương trình đường chuẩn thực nghiệm mô phỏng $\ln(R) = a + b\rho$ với hệ số tương quan tuyến tính cao.
  • Phát triển mô hình mạng nơ-ron nhân tạo MLP sử dụng hàm kích hoạt Softsign huấn luyện trên Mathematica, có khả năng nhận dạng và ngoại suy mật độ chất lỏng tức thời chỉ từ hai tham số hình học và bức xạ ($R, d$).
  • Thực nghiệm đo đạc thành công 7 hóa chất acid và dung môi thông dụng trên hệ đầu dò $\text{NaI(Tl)}$ với nguồn $^{137}\text{Cs}$ hoạt độ thấp ($0,76\ \mu\text{Ci}$).

Đóng góp khoa học và kỹ thuật so với các công trình đi trước:

  • Khắc phục hạn chế của Priyada và cộng sự (2018): Thay vì dùng đường chuẩn từ thực nghiệm dễ phát sinh sai số lớn (từng bị lệch tới 21,54% đối với hexane), khóa luận dùng dữ liệu mô phỏng Monte Carlo chuẩn hóa theo tỷ số $R$, giảm thiểu triệt để sai số do môi trường.
  • Khắc phục hạn chế của Chuong và cộng sự (2020, 2021): Loại bỏ hoàn toàn sự phụ thuộc vào việc phải xây dựng lại đường chuẩn cho từng đường kính ống chứa cụ thể. Mô hình ANN đóng vai trò là một hàm nội suy tổng quát hóa, cho phép xác định mật độ ở bất kỳ kích thước ống nào nằm trong dải khảo sát mà không cần chạy lại mô phỏng MCNP6 tốn nhiều thời gian.
  • Cải thiện an toàn bức xạ: Sử dụng nguồn phóng xạ $^{137}\text{Cs}$ hoạt độ rất thấp ($0,76\ \mu\text{Ci}$) kết hợp đầu dò nhấp nháy phổ thông $\text{NaI(Tl)}$ thay vì đầu dò bán dẫn siêu tinh khiết HPGe đắt tiền hoặc các nguồn hoạt độ cao, nâng cao tính khả thi khi ứng dụng vào phòng thí nghiệm học đường và dây chuyền công nghiệp.

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

Khóa luận ghi nhận một số hạn chế kỹ thuật:

  • Hiện tượng dịch chuyển đỉnh phổ (peak drift) do nhiệt độ phòng biến thiên trong các phép đo kéo dài 8 giờ vẫn gây khó khăn cho việc xử lý số liệu, đòi hỏi phải can thiệp thủ công bằng phần mềm Colegram để điều chỉnh vùng ROI theo từng phổ.
  • Dữ liệu thực nghiệm mới chỉ được tiến hành kiểm chứng trên một kích thước hình học duy nhất là ống thủy tinh có đường kính trong 2,68 cm; chưa có điều kiện đo thực nghiệm trực tiếp trên 19 kích cỡ đường kính còn lại do hạn chế về trang thiết bị ống đo tại phòng thí nghiệm.
  • Dải chất lỏng mô phỏng chỉ mới giới hạn trong khoảng mật độ từ 0,78 $\text{g/cm}^3$ đến 1,84 $\text{g/cm}^3$, chưa bao quát các chất lỏng có mật độ rất cao (như thủy ngân, bromoform) hoặc các hệ chất lỏng đa pha, huyền phù và nhũ tương.

Hướng phát triển tiếp theo được đề xuất:

  • Mở rộng hệ thống đo thực nghiệm với nhiều loại vật liệu ống (thép không gỉ, nhựa PVC, thủy tinh) và nhiều kích thước đường kính khác nhau để kiểm định toàn diện khả năng thích ứng của mạng nơ-ron.
  • Nâng cấp hệ thống đo bằng các mạch ổn định nhiệt hoặc phần mềm tự động hiệu chỉnh dịch đỉnh phổ trực tiếp theo thời gian thực nhằm rút ngắn thời gian tích lũy xung.
  • Mở rộng tập dữ liệu huấn luyện cho mạng ANN với các chùm bức xạ đa năng lượng hoặc kết hợp đồng thời phương pháp truyền qua và tán xạ gamma để giải quyết bài toán nhận diện chất lỏng có bọt khí hoặc dòng chảy đa pha trong đường ống dẫn dầu.

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 và học viên cao học chuyên ngành Sư phạm Vật lý, Vật lý hạt nhân, Kỹ thuật hạt nhân: Cung cấp phương pháp luận chi tiết về việc kết hợp kỹ thuật mô phỏng hạt nhân MCNP6 với các công cụ xử lý phổ chuyên dụng (Colegram) và hệ đo thực nghiệm phổ kế gamma $\text{NaI(Tl)}$.
  • Các kỹ sư trong lĩnh vực kiểm tra không phá hủy (NDT) và tự động hóa công nghiệp: Cung cấp quy trình hoàn chỉnh để phát triển hệ đo tỷ trọng trực tuyến trên đường ống công nghiệp, ứng dụng thuật toán học máy để giải bài toán định lượng không tiếp xúc đối với chất lỏng độc hại hoặc dễ cháy nổ.
  • Nội dung đáng tham khảo nhất: Phương pháp chuẩn hóa tỷ số cường độ $R$ để loại trừ hiệu suất đầu dò; kỹ thuật xử lý vùng ROI $[\text{đỉnh} \pm 2,5\text{FWHM}]$ để chống trôi đỉnh; và quy trình thiết kế cấu trúc mạng MLP 2 đầu vào trên Mathematica.

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

1. Tại sao nghiên cứu lại sử dụng tỷ số cường độ tương đối $R$ thay vì sử dụng trực tiếp cường độ chùm tia truyền qua $I$ để lập đường chuẩn?
Việc sử dụng cường độ tuyệt đối $I$ làm đường chuẩn khiến kết quả bị phụ thuộc chặt chẽ vào hiệu suất ghi nhận của đầu dò bức xạ, vốn thay đổi liên tục theo nhiệt độ môi trường, độ trôi điện tử và hình học đo. Bằng cách lấy tỷ số giữa cường độ truyền qua mẫu chất lỏng và cường độ truyền qua nước cất trong cùng điều kiện hình học ($R = N_{\text{chất lỏng}} / N_{\text{nước}}$), hiệu suất tuyệt đối của đầu dò được triệt tiêu hoàn toàn. Nhờ đó, giá trị $R$ tính từ mô phỏng Monte Carlo ($R_{\text{sim}}$) tương đương với giá trị đo được trên thực nghiệm ($R_{\text{exp}}$).

2. Cấu hình hệ đo thực nghiệm sử dụng loại nguồn phóng xạ và đầu dò nào?
Nghiên cứu sử dụng nguồn phóng xạ đồng vị $^{137}\text{Cs}$ phát tia gamma đơn năng 661,7 keV với xác suất phát 84,99%, hoạt độ tại thời điểm đo là $(0,76 \pm 0,03)\ \mu\text{Ci}$, đặt trong khối chì bảo vệ có kích thước $12,39\text{ cm} \times 12,39\text{ cm} \times 5,48\text{ cm}$ cùng ống chuẩn trực chì đường kính 0,96 cm. Đầu dò sử dụng là loại nhấp nháy $\text{NaI(Tl)}$ của hãng AMPTEK với tinh thể kích thước $7,62\text{ cm} \times 7,62\text{ cm}$, chuẩn trực đầu dò đường kính 1,00 cm và bộ phân tích đa kênh MCA cài đặt 8192 kênh.

3. Tác giả đã xử lý hiện tượng dịch chuyển đỉnh phổ (peak drift) trong quá trình đo thực nghiệm như thế nào?
Do mỗi mẫu chất lỏng phải đo tích lũy xung liên tục trong 8 giờ để đảm bảo thống kê, sự thay đổi nhiệt độ môi trường và nhiễu vi mạch làm đỉnh phổ bị trượt sang trái hoặc phải. Tác giả không dùng một vùng kênh ROI cố định mà sử dụng phần mềm Colegram để khớp đỉnh Gauss và phông One Step cho từng phổ, sau đó cắt vùng ROI động theo công thức $[\text{vị trí đỉnh} - 2,5\text{FWHM}; \text{vị trí đỉnh} + 2,5\text{FWHM}]$. Cách xử lý này đảm bảo lấy trọn vẹn diện tích đỉnh hấp thụ toàn phần mà không bị ảnh hưởng bởi hiện tượng trôi kênh.

4. Mô hình mạng nơ-ron nhân tạo (ANN) trong khóa luận có cấu trúc và hàm kích hoạt như thế nào?
Mô hình ANN thuộc dạng Perceptron đa lớp (MLP) gồm 3 lớp: lớp đầu vào có 2 nơ-ron (nhận 2 biến số là tỷ số cường độ $R$ và đường kính ống $d$), một lớp ẩn duy nhất với số nơ-ron được tối ưu hóa dựa trên chỉ số RMSE và $R^2$, và lớp đầu ra có 1 nơ-ron (xuất giá trị mật độ dự đoán $\rho$). Hàm kích hoạt phi tuyến được lựa chọn cho mạng là hàm Softsign với công thức $f(x) = x / (1 + |x|)$, được lập trình và huấn luyện trên phần mềm Wolfram Mathematica.

5. Dữ liệu huấn luyện mạng ANN được tạo ra từ đâu và có quy mô bao nhiêu?
Dữ liệu huấn luyện mạng ANN hoàn toàn không lấy từ thực nghiệm mà được tạo ra từ chương trình mô phỏng Monte Carlo MCNP6 nhằm tiết kiệm thời gian đo đạc và dễ dàng thay đổi thông số hình học. Tập dữ liệu mô phỏng gồm tổng cộng 304 điểm dữ liệu tương ứng với 16 hợp chất lỏng chuẩn được khảo sát trên 20 kích thước đường kính ống từ 1,0 cm đến 10,0 cm. Tập dữ liệu này được phân chia ngẫu nhiên thành 85% cho tập huấn luyện (training) và 15% cho tập kiểm tra (testing).

Kết luận

Khóa luận tốt nghiệp của tác giả Nguyễn Thu Hằng đã giải quyết thành công bài toán ước lượng mật độ chất lỏng không tiếp xúc bằng cách tích hợp phương pháp bức xạ gamma truyền qua với mô phỏng Monte Carlo MCNP6 và mạng nơ-ron nhân tạo MLP. Việc ứng dụng mạng nơ-ron huấn luyện trên 304 điểm dữ liệu mô phỏng đã khắc phục triệt để hạn chế phải lập lại đường chuẩn cho từng kích thước ống chứa trong các nghiên cứu trước đây. Kết quả đo kiểm chứng thực nghiệm trên 7 mẫu chất lỏng bằng đầu dò $\text{NaI(Tl)}$ và nguồn $^{137}\text{Cs}$ hoạt độ thấp khẳng định tính đúng đắn, độ an toàn và tiềm năng ứng dụng thực tế của mô hình trong kiểm soát chất lượng công nghiệp.