Tổng quan về luận án

Kiểm tra không phá hủy (Non-Destructive Testing – NDT) đóng vai trò then chốt trong việc bảo đảm an toàn cấu trúc của các công trình công nghiệp, từ đường ống dẫn dầu khí, thân tàu biển cho đến thành lò nung chịu lửa hoạt động trong điều kiện nhiệt độ và áp suất khắc nghiệt. Trong số các kỹ thuật NDT, phương pháp gamma tán xạ nổi bật nhờ ưu thế vượt trội: nguồn phát và đầu dò được bố trí cùng một phía, cho phép kiểm tra đơn diện và hoàn toàn không tiếp xúc. Tuy nhiên, rào cản kỹ thuật lớn nhất khi triển khai hệ đo di động tại hiện trường nằm ở sự đánh đổi giữa an toàn bức xạ và độ chính xác phân tích. Việc sử dụng nguồn phóng xạ hoạt độ thấp (như $^{137}\text{Cs}$ $5\text{ mCi}$) kết hợp đầu dò nhấp nháy $\text{NaI(Tl)}$ nhằm bảo đảm tính cơ động đòi hỏi phải mở rộng khẩu độ ống chuẩn trực, dẫn đến sự gia tăng mạnh mẽ của các photon tán xạ nhiều lần (multiple scattering) chồng chập phức tạp lên đỉnh tán xạ một lần (single scattering).

Khoảng trống nghiên cứu (research gap) cốt lõi được xác định là: các nghiên cứu trước đây (như Dumond 1930; Felsteiner và cộng sự 1975; Fernández 1991; Singh và cộng sự 2006; Tạ Duy Tám và cộng sự 2017) hoặc chỉ dừng lại ở các phép tính Monte Carlo lý thuyết thuần túy, hoặc đơn giản hóa phổ tán xạ nhiều lần bằng các hàm đa thức bậc bốn và hàm Gauss đơn lẻ chưa phản ánh đúng bản chất vật lý của các bậc tán xạ bậc cao, đồng thời thiếu vắng một công thức giải tích thống kê định lượng để xác định bề dày bão hòa ($T_s$).

Luận án tiến sĩ Vật lý nguyên tử và hạt nhân (Mã số: 62 44 05 01) của nghiên cứu sinh Võ Hoàng Nguyên, dưới sự hướng dẫn của GS. Châu Văn Tạo tại Trường Đại học Khoa học Tự nhiên – ĐHQG TP. Hồ Chí Minh (2023), đã giải quyết triệt để bài toán này thông qua hệ thống câu hỏi và giả thuyết nghiên cứu tường minh:

  • Câu hỏi nghiên cứu 1 (RQ1): Dạng phân bố năng lượng nội tại của từng thành phần tán xạ một lần, hai lần và trên hai lần biến thiên như thế nào theo hình học bia (phẳng và trụ rỗng)?
  • Câu hỏi nghiên cứu 2 (RQ2): Làm thế nào để xây dựng một quy trình bán thực nghiệm chuẩn hóa, phân tách chính xác các thành phần tán xạ từ phổ đo thực nghiệm bằng đầu dò có độ phân giải năng lượng thấp $\text{NaI(Tl)}$?
  • Câu hỏi nghiên cứu 3 (RQ3): Có thể thiết lập một mô hình toán học giải tích dựa trên bản chất thống kê phân bố Gauss để xác định chính xác bề dày bão hòa ($T_s$) và độ dày thành vật liệu ($T$) hay không?

Các giả thuyết nghiên cứu được kiểm chứng gồm:

  • Giả thuyết $H_1$: Phổ tán xạ hai lần có thể mô hình hóa chính xác bằng tổ hợp ba hàm phân bố Gauss $[G_1(x) + G_2(x) + G_3(x)]$, phản ánh các cơ chế tán xạ Compton-Compton, Compton-Rayleigh và Rayleigh-Compton.
  • Giả thuyết $H_2$: Tán xạ trên hai lần tuân theo một hàm mật độ xác suất liên tục $P(x)$ trải rộng về vùng năng lượng thấp.
  • Giả thuyết $H_3$: Việc loại trừ thành phần tán xạ nhiều lần bằng quy trình bán thực nghiệm sẽ cho phép xác định bề dày thép và nhôm với sai số tương đối dưới 4% trên toàn dải đo thực nghiệm.

Khung lý thuyết của luận án tích hợp sâu sắc giữa lý thuyết tán xạ lượng tử Klein–Nishina, phương trình vận chuyển photon Boltzmann và lý thuyết thăng giáng thống kê bức xạ. Luận án khảo sát thực nghiệm trên các bia nhôm 5052 (dày 2,0 mm đến 100,5 mm), thép C45 dạng tấm phẳng (dày 3,7 mm đến 32,2 mm) và ống thép C45 dạng trụ rỗng (đường kính ngoài 273 mm, dày 1,9 mm đến 12,0 mm; mô phỏng mở rộng đến 22,0 mm) với thời gian đo từ 14.400 s đến 18.000 s. Kết quả nghiên cứu mang lại đóng góp đột phá khi hạ sai số xác định bề dày thực tế xuống dưới 3,5%, cung cấp giải pháp NDT định lượng chính xác cao cho ngành công nghiệp kiểm định cơ khí và hóa dầu.

Literature Review và Positioning

Nghiên cứu về bức xạ gamma tán xạ đã trải qua gần một thế kỷ phát triển với nhiều trường phái tiếp cận lý thuyết, mô phỏng và thực nghiệm. Ngay từ năm 1930, Dumond đã đặt nền móng giải tích khi ứng dụng lý thuyết tán xạ cổ điển Thompson kết hợp hiệu chỉnh Breit cho bức xạ phân cực để tính toán cường độ tán xạ hai lần trên mẫu khối cầu. Đến năm 1937, Kirkpatrick mở rộng mô hình cho góc tán xạ $90^\circ$ trên mẫu hình trụ. Tuy nhiên, các giải pháp giải tích cổ điển này nhanh chóng bộc lộ hạn chế khi xử lý các môi trường tán xạ thực tế đa nguyên tố và hình học bất đối xứng.

Sự phát triển của kỹ thuật tính toán đã mở ra kỷ nguyên mô phỏng ngẫu nhiên. Felsteiner và cộng sự (1975) đã công bố các tính toán Monte Carlo xác định đóng góp của tán xạ nhiều lần trong phép đo Compton, chỉ ra vai trò quyết định của quãng đường tự do trung bình photon. Halonen và cộng sự (1979) so sánh tiết diện tán xạ phi tương đối Thompson, Klein–Nishina và tiết diện Ribberfors (hiệu chỉnh phân bố động lượng electron), chứng minh rằng tại mức năng lượng gamma trung bình ($>200\text{ keV}$), tiết diện tương đối tính Klein–Nishina là hoàn toàn phù hợp để tính toán tán xạ bậc cao.

Bước đột phá về lý thuyết vận chuyển giải tích thuộc về Fernández (1991, 1999), khi giải phương trình vector Boltzmann để phân tích chi tiết bốn kênh tương tác của tán xạ hai lần: Rayleigh-Rayleigh, Compton-Rayleigh, Rayleigh-Compton và Compton-Compton. Fernández chứng minh rằng thành phần Compton-Compton tạo ra một phổ năng lượng liên tục trải rộng từ $E'/(1 + 2E'/m_e c^2)$ đến $E$, trong khi Compton-Rayleigh và Rayleigh-Compton đạt đỉnh cực đại ngay tại năng lượng tán xạ một lần $E'$, gây ra hiện tượng chồng lấn phổ nghiêm trọng. Gần đây hơn, Tarim và cộng sự (2013) sử dụng Monte Carlo để tách phổ từ bậc hai đến bậc chín, minh chứng sự phụ thuộc phi tuyến của tán xạ ngược nhiều lần vào bề dày và năng lượng chùm tia.

Trong thực nghiệm, các khảo sát của Paramesh và cộng sự (1983), Singh và cộng sự (2006, 2007) và Kiran và cộng sự (2016) trên các bia C, Al, Fe, Cu, Pb ở các góc tán xạ $40^\circ - 180^\circ$ đã khẳng định: cường độ tán xạ nhiều lần tăng dần theo bề dày bia và đạt tới giá trị bão hòa ($T_s$). Tuy nhiên, tồn tại một tranh luận học thuật và hạn chế cố hữu trong các công trình này:

  • Quan điểm thứ nhất (Nhóm Singh et al. 2006; Kiran et al. 2016): Ghi nhận hiện tượng bão hòa thực nghiệm nhưng không cung cấp được công thức toán học tường minh để định lượng giá trị $T_s$, dẫn đến việc xác định $T_s$ mang tính trực quan, thiếu cơ sở thống kê.
  • Quan điểm thứ hai (Nhóm ứng dụng NDT: Mullin et al. 1994; Priyada et al. 2011, 2012): Khắc phục nhiễu tán xạ nhiều lần bằng cách sử dụng đầu dò bán dẫn siêu tinh khiết $\text{HPGe}$ có độ phân giải năng lượng cực cao. Dù đạt độ chính xác tốt, phương pháp này đòi hỏi hệ thống làm lạnh bằng nitơ lỏng cồng kềnh, triệt tiêu tính cơ động tại hiện trường công nghiệp.
  • Quan điểm thứ ba (Nhóm Tạ Duy Tám et al. 2017; Huỳnh Đình Chương et al. 2017): Sử dụng đầu dò $\text{NaI(Tl)}$ và làm khớp phổ Compton bằng một hàm Gauss cho tán xạ một lần, một Gauss cho tán xạ hai lần và đa thức bậc bốn cho tán xạ trên hai lần. Tuy nhiên, cách tiếp cận này chưa phản ánh đúng bản chất vật lý đa đỉnh của tán xạ hai lần và phân bố phi đối xứng của tán xạ bậc cao.

So sánh với hai nghiên cứu quốc tế điển hình gần đây:

  1. Priyada và cộng sự (2011): Sử dụng hệ đo gamma tán xạ với nguồn $^{137}\text{Cs}$ chuẩn trực hẹp và đầu dò $\text{HPGe}$ để đánh giá ăn mòn thép mềm, đạt độ nhạy cao nhưng chi phí thiết bị đắt đỏ và không thể áp dụng cho các cấu hình đo nhanh tại hiện trường.
  2. Salgado và cộng sự (2021): Sử dụng mô phỏng MCNP6 kết hợp mạng nơ-ron nhân tạo (ANN) với ba đầu dò $\text{NaI(Tl)}$ để dự đoán độ dày cáu cặn $\text{BaSO}_4$ trong đường ống dẫn dầu. Phương pháp của Salgado mang tính chất hộp đen (black-box), phụ thuộc lớn vào tập dữ liệu huấn luyện và thiếu tính giải thích cơ chế vật lý tương tác.

Luận án của Võ Hoàng Nguyên định vị chính xác vào khoảng trống khoa học này: kết hợp sức mạnh mô phỏng vi mô của GEANT4 để giải mã bản chất các phổ thành phần, kiểm chứng bằng MCNP6 và chuẩn hóa thành một quy trình bán thực nghiệm tường minh, vận hành trên phần cứng gọn nhẹ $\text{NaI(Tl)}$ với nguồn $^{137}\text{Cs}$ $5\text{ mCi}$.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã mở rộng và hoàn thiện khung lý thuyết tương tác bức xạ gamma trong môi trường vật chất giới hạn thông qua ba đóng góp then chốt:

  1. Mô hình hóa giải tích quá trình suy giảm và tán xạ một lần: Xuất phát từ thông lượng photon tới tại điểm tán xạ $P$ bên trong vật liệu: $$I_1(P) = I_0 \exp\left(-\left[\frac{\mu}{\rho}\right] \rho x\right) = I_0 \exp(-\mu(E) x)$$ và cường độ tán xạ Compton vi phân theo góc $\theta$ xác định bởi tiết diện vi phân Klein–Nishina: $$\frac{d\sigma(E, \theta)}{d\Omega} = \frac{r_0^2}{2} \left[\frac{E'}{E}\right]^2 \left[\frac{E}{E'} + \frac{E'}{E} - \sin^2\theta\right]$$ Luận án thiết lập phương trình tích phân cường độ chùm gamma tán xạ một lần ghi nhận tại đầu dò từ bản phẳng có bề dày $T$: $$I = I_s \left[1 - \exp(-\mu_{\text{eff}} T)\right]$$ Trong đó: $$\mu_{\text{eff}} = \mu(E)\sec\theta_1 + \mu(E')\sec\theta_2$$ $$I_s = \frac{k I_0}{\mu(E)\sec\theta_1 + \mu(E')\sec\theta_2} = \frac{k I_0}{\mu_{\text{eff}}}$$ Hệ số $k = \rho_e \left[\frac{d\sigma(E, \theta)}{d\Omega}\right] S(E, \theta, Z) d\Omega$. Từ đó, công thức xác định bề dày vật liệu được dẫn xuất tường minh: $$T = -\frac{1}{\mu_{\text{eff}}} \ln\left(1 - \frac{I}{I_s}\right)$$

  2. Dẫn xuất công thức thống kê định lượng bề dày bão hòa ($T_s$): Đây là đóng góp lý thuyết breakthrough của luận án. Bề dày bão hòa $T_s$ được định nghĩa chặt chẽ là bề dày nhỏ nhất mà tại đó sự gia tăng cường độ tán xạ không còn phân biệt được với thăng giáng thống kê của phép đo ở mức độ tin cậy xác định ($1 - \alpha = 95%$, tương ứng $n = 1{,}96$ theo chuẩn phân bố Gauss): $$I(T_s) + n \sqrt{I(T_s)} \ge I_s$$ Giải bất phương trình phi tuyến này dẫn đến công thức giải tích chuẩn xác: $$T_s = -\frac{1}{\mu_{\text{eff}}} \ln\left[1 - \left(1 - \frac{n}{\sqrt{I_s}}\right)^2\right]$$ Kèm theo công thức lan truyền sai số toàn phần: $$u(T_s) = \sqrt{\left(\frac{\partial T_s}{\partial \mu_{\text{eff}}}\right)^2 u^2(\mu_{\text{eff}}) + \left(\frac{\partial T_s}{\partial I_s}\right)^2 u^2(I_s)}$$

  3. Chuyển dịch Paradigm trong biểu diễn phổ tán xạ bậc cao: Thay vì xem tán xạ nhiều lần là một hàm nền trơn (smooth background), luận án chứng minh sự tồn tại của cấu trúc đa đỉnh vi mô trong phổ tán xạ hai lần và phân bố bất đối xứng của phổ tán xạ bậc cao ($>2$ lần), xác lập cơ chế phân rã phổ chính xác về mặt vật lý hạt nhân.

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp đồng thời ba trụ cột lý thuyết và phương pháp:

+-----------------------------------------------------------------------------------+
|                            KHUNG PHÂN TÍCH ĐỘC ĐÁO                                |
+-----------------------------------------------------------------------------------+
|  1. Lý thuyết Tương tác Bức xạ         2. Mô phỏng Đa mã       3. Phân tích Dữ liệu|
|  - Tiết diện Klein-Nishina             - GEANT4 (Theo vết,     - Quy trình Bán     |
|  - Phương trình Vận chuyển Boltzmann     phân rã bậc tán xạ)     thực nghiệm       |
|  - Thăng giáng thống kê Gauss (95% CI) - MCNP6 (F8:P, FT8 GEB) - Làm khớp phi tuyến|
+-----------------------------------------------------------------------------------+
+-----------------------------------------------------------------------------------+
|                        MÔ HÌNH HÀM TOÁN HỌC ĐẶC TRƯNG                             |
+-----------------------------------------------------------------------------------+
|  * Tán xạ 1 lần:    S(x) -> Hàm Gauss đơn đỉnh chuẩn tại E' (288 keV)             |
|  * Tán xạ 2 lần:    D(x) = G1(x) + G2(x) + G3(x) (Triple-Gaussian Model)          |
|  * Tán xạ >2 lần:   M(x) = h3 * P(x) (Hàm phân bố xác suất bất đối xứng)          |
+-----------------------------------------------------------------------------------+
+-----------------------------------------------------------------------------------+
|                        KẾT QUẢ ĐẦU RA ĐỊNH LƯỢNG                                  |
+-----------------------------------------------------------------------------------+
|  - Tách triệt để cường độ tán xạ 1 lần (I)                                       |
|  - Xác định chính xác Bề dày bão hòa (Ts) và Bề dày thành ống (T)                |
+-----------------------------------------------------------------------------------+

Điều kiện biên (boundary conditions) của mô hình được xác định rõ: áp dụng cho chùm gamma đơn năng ($^{137}\text{Cs}$, $E_0 = 661{,}657\text{ keV}$), góc tán xạ hình học cố định $120^\circ \pm 0{,}1^\circ$, vật liệu kim loại đồng nhất (Al 5052, thép C45) trong môi trường không khí tiêu chuẩn.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ triết lý thực chứng thực nghiệm (positivism) kết hợp mô hình hóa tính toán đa cấp (multi-level computational modeling). Thiết kế nghiên cứu gồm ba tầng bổ trợ lẫn nhau:

  • Tầng 1 (Vi mô - Cơ chế tương tác): Sử dụng bộ công cụ GEANT4 để theo vết từng photon, phân lập chính xác các sự kiện tán xạ theo số lần va chạm (1 lần, 2 lần, trên 2 lần).
  • Tầng 2 (Mô phỏng Hệ đo hoàn chỉnh): Sử dụng MCNP6 để mô hình hóa toàn diện hình học thực tế của nguồn, ống chuẩn trực chì, bia tán xạ và tinh thể đầu dò $\text{NaI(Tl)}$.
  • Tầng 3 (Thực nghiệm kiểm chứng): Đo đạc thực tế trên hệ thiết bị chuẩn tại Phòng thí nghiệm Kỹ thuật Hạt nhân, Trường Đại học Khoa học Tự nhiên – ĐHQG TP. HCM.

Quy trình nghiên cứu rigorous

  1. Thiết lập phần cứng và chuẩn trực thực nghiệm:
  • Nguồn phóng xạ: $^{137}\text{Cs}$ (mã hiệu chuẩn Eckert & Ziegler, hoạt độ $5\text{ mCi}$ tại thời điểm 18/04/2013), phát tia gamma đơn năng $661{,}657\text{ keV}$ (xác suất phát $84{,}99%$) và tia X đặc trưng $31{,}8 - 37{,}3\text{ keV}$. Nguồn được bọc trong vỏ thép trụ ($\varnothing 6{,}0\text{ mm} \times 8{,}0\text{ mm}$) và đặt trong hộp bảo vệ bằng chì ($\varnothing_{\text{ngoài}} = 130{,}0\text{ mm}$, $\varnothing_{\text{trong}} = 10{,}0\text{ mm}$ và $20{,}0\text{ mm}$, chiều dài $152\text{ mm}$).
  • Đầu dò và điện tử hạt nhân: Đầu dò nhấp nháy $\text{NaI(Tl)}$ Model 802 (Canberra Inc.), tinh thể kích thước lớn $76{,}2\text{ mm} \times 76{,}2\text{ mm}$ ($3'' \times 3''$) tích hợp ống nhân quang điện (PMT). Khối xử lý số tích hợp Canberra Osprey kết nối phần mềm Genie 2000 phiên bản 3.3. Đầu dò được che chắn bằng vỏ chì dày ($143\text{ mm} \times 112\text{ mm} \times 112\text{ mm}$), nắp chuẩn trực có lỗ tròn $\varnothing 10\text{ mm}$ (đối với bia phẳng) hoặc $\varnothing 20\text{ mm}$ (đối với ống trụ).
  • Cấu hình hình học: Bàn tán xạ tự động điều khiển góc bằng động cơ bước, thiết lập góc tán xạ cố định $120^\circ \pm 0{,}1^\circ$. Khoảng cách từ nguồn đến tâm bia: $200{,}0 \pm 0{,}5\text{ mm}$ (bia phẳng) và $284{,}3 \pm 0{,}5\text{ mm}$ (ống trụ); khoảng cách từ tâm bia đến bề mặt đầu dò: $180{,}0 \pm 0{,}5\text{ mm}$ (bia phẳng) và $203{,}3 \pm 0{,}5\text{ mm}$ (ống trụ).
                        [ Nguồn 137Cs trong hộp chì ]
                                      \
                                       \ Chùm tia tới (E = 662 keV)
                                        \
                                         v
   [ Ống thép / Bia phẳng ] <---- (Tâm tán xạ P)
                                         /
                                        / Chùm tán xạ (E' = 288 keV)
                                       /  (Góc tán xạ theta = 120 độ)
                                      v
                        [ Đầu dò NaI(Tl) + Ống chuẩn trực ]
                        [ Khối MCA số Canberra Osprey ]
                        [ Phần mềm phân tích Genie 2000 ]
  1. Quy trình mô phỏng Monte Carlo:
  • Trong GEANT4 (phiên bản 10.x): Xây dựng các lớp G4VSolid, G4LogicalVolume và gán bia tán xạ là Sensitive Detector. Thông qua đối tượng Hit và Track, chương trình trích xuất riêng biệt phân bố năng lượng của các photon sau đúng 1 lần va chạm, 2 lần va chạm và $>2$ lần va chạm trong bia.
  • Trong MCNP6: Thiết lập tệp vào chi tiết với số lịch sử hạt khổng lồ: $10\text{ tỷ phép thử}$ ($10^{10}$ NPS) cho mỗi phổ mô phỏng nhằm triệt tiêu sai số thống kê Monte Carlo ($<0{,}1%$). Sử dụng thẻ đánh giá độ cao xung F8:P kết hợp thẻ mở rộng năng lượng Gauss FT8 GEB để tái tạo hiệu ứng thăng giáng điện tích và nhiễu điện tử của đầu dò thực tế: $$\text{FWHM}(E) = a + b\sqrt{E} + c E^2$$ với các hệ số hiệu chuẩn thực nghiệm xác định chính xác: $a = -0{,}0036\text{ MeV}$, $b = 0{,}0621\text{ MeV}^{1/2}$, $c = 0{,}3146\text{ MeV}^{-1}$. Thư viện tiết diện tương tác photon ENDF/B-VI được sử dụng nhất quán.

Data và phân tích

Mẫu khảo sát gồm:

  • 13 mẫu nhôm 5052 (thành phần chuẩn TCVN 12513-7:2018: $2{,}2-2{,}5%\text{ Mg}$, $<0{,}4%\text{ Fe}$, $<0{,}25%\text{ Si}$, còn lại là $\text{Al}$) với độ dày từ $2{,}0\text{ mm}$ đến $100{,}5\text{ mm}$.
  • 10 mẫu thép C45 dạng tấm phẳng (TCVN 12142-1:2017: $0{,}42-0{,}50%\text{ C}$, $0{,}50-0{,}80%\text{ Mn}$, $0{,}10-0{,}40%\text{ Si}$, còn lại là $\text{Fe}$) với độ dày từ $3{,}7\text{ mm}$ đến $32{,}2\text{ mm}$.
  • Các đoạn ống thép C45 công nghiệp đường kính ngoài chuẩn $273 \pm 0{,}1\text{ mm}$, dài $250 \pm 0{,}5\text{ mm}$, bề dày thành ống từ $1{,}9\text{ mm}$ đến $12{,}0\text{ mm}$ (và mô phỏng dải dày $2{,}0 - 22{,}0\text{ mm}$). Bề dày được kiểm chuẩn bằng thước kẹp điện tử Mitutoyo độ chính xác $\pm 0{,}01\text{ mm}$.

Dữ liệu phổ thực nghiệm và mô phỏng được phân tích bằng thuật toán tối ưu hóa phi tuyến Levenberg-Marquardt, kiểm tra độ tương thích bằng chỉ số chi bình phương rút gọn ($\chi^2/\text{ndf} \approx 1{,}0$) và hệ số tương quan $R^2 > 0{,}998$.

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Phát hiện cấu trúc vi mô của phổ tán xạ hai lần: Nghiên cứu GEANT4 chỉ ra rằng phổ tán xạ hai lần không thể mô tả bằng một hàm Gauss đơn giản. Do sự kết hợp của các hiệu ứng Compton-Compton, Compton-Rayleigh và Rayleigh-Compton, phổ tán xạ hai lần phân bố trải dài từ $180\text{ keV}$ đến trên $400\text{ keV}$ và được làm khớp hoàn hảo bởi tổ hợp ba hàm Gauss $D(x) = G_1(x) + G_2(x) + G_3(x)$ với độ lệch chuẩn cực tiểu.

  2. Quy luật bất biến của hàm mật độ xác suất tán xạ trên các ống thép: Một phát hiện mang tính quy luật là: đối với các ống thép C45 có bề dày từ $2{,}0\text{ mm}$ đến $22{,}0\text{ mm}$, phổ mật độ xác suất của từng thành phần tán xạ (chuẩn hóa diện tích bằng 1) gần như đồng nhất và bất biến theo bề dày. Cụ thể:

  • Phổ mật độ xác suất trung bình của tán xạ một lần $S_0(x)$ có đỉnh cực đại tại kênh tương ứng năng lượng lý thuyết $E' = 288\text{ keV}$.
  • Phổ mật độ xác suất trung bình của tán xạ hai lần $D_0(x)$ và tán xạ trên hai lần $P_0(x)$ giữ nguyên dạng hình học, chỉ có biên độ tỷ lệ thay đổi theo độ dày thành ống.

  1. Xác lập Quy trình Bán thực nghiệm 3 thành phần hoàn chỉnh: Phổ tán xạ tổng cộng $Y(x)$ trên ống thép được phân tách chính xác theo phương trình: $$Y(x) = h_1 S_0(x) + h_2 D_0(x) + h_3 P_0(x)$$ Trong đó $h_1, h_2, h_3$ là các trọng số diện tích tương ứng của thành phần tán xạ một lần, hai lần và trên hai lần. Quy trình cho phép trích xuất số đếm tán xạ một lần thực sự ($I = h_1$) mà không bị méo dạng bởi phông tán xạ bậc cao.

  2. Độ chính xác vượt trội trong xác định bề dày:

  • Đối với bia nhôm phẳng: Bề dày bão hòa thực nghiệm xác định được là $T_s = 84{,}5 \pm 2{,}1\text{ mm}$; sai số xác định bề dày thực tế $<3{,}8%$.
  • Đối với bia thép C45 phẳng: Bề dày bão hòa $T_s = 28{,}2 \pm 0{,}9\text{ mm}$; sai số xác định bề dày $<3{,}2%$.
  • Đối với ống thép C45 ($\varnothing 273\text{ mm}$): Xử lý phổ mô phỏng MCNP6 cho độ lệch bề dày $<2{,}0%$; xử lý phổ đo thực nghiệm cho độ lệch bề dày $<3{,}5%$ trên toàn dải từ $1{,}9\text{ mm}$ đến $12{,}0\text{ mm}$.
+---------------------------------------------------------------------------------------------+
|               BẢNG SO SÁNH KẾT QUẢ ĐO BỀ DÀY VÀ ĐẶC TRƯNG HỆ THỐNG                          |
+----------------------+--------------------------+--------------------+----------------------+
| Loại mẫu / Cấu hình  | Bề dày thực tế / Dải đo  | Bề dày bão hòa Ts  | Độ sai lệch tối đa   |
+----------------------+--------------------------+--------------------+----------------------+
| Tấm nhôm phẳng 5052  | 2,0 mm - 100,5 mm        | 84,5 ± 2,1 mm      | < 3,8 %              |
| Tấm thép phẳng C45   | 3,7 mm - 32,2 mm         | 28,2 ± 0,9 mm      | < 3,2 %              |
| Ống thép C45 (MCNP6) | 2,0 mm - 22,0 mm         | 26,8 ± 0,8 mm      | < 2,0 %              |
| Ống thép C45 (Thực tế)| 1,9 mm - 12,0 mm        | 26,5 ± 0,9 mm      | < 3,5 %              |
+----------------------+--------------------------+--------------------+----------------------+

Implications đa chiều

  • Về mặt Lý thuyết & Phương pháp luận: Luận án cung cấp một phương pháp luận chuẩn mực kết hợp giữa mô phỏng vi mô (GEANT4) và vĩ mô (MCNP6) để giải mã các bài toán chồng phổ phức tạp trong vật lý bức xạ ứng dụng.
  • Về mặt Kỹ thuật Thực tiễn: Cho phép chế tạo các thiết bị NDT xách tay sử dụng đầu dò $\text{NaI(Tl)}$ giá thành thấp, hoạt động bền bỉ ở nhiệt độ thường, thay thế hoàn toàn nhu cầu sử dụng đầu dò $\text{HPGe}$ trong kiểm tra bề dày và ăn mòn thành ống.
  • Về mặt An toàn Bức xạ: Việc thuật toán xử lý phổ tối ưu hóa cho phép giảm hoạt độ nguồn phóng xạ $^{137}\text{Cs}$ xuống mức thấp ($5\text{ mCi}$), giảm thiểu trọng lượng vỏ chì che chắn, bảo đảm an toàn bức xạ tối đa cho nhân viên vận hành tại hiện trường theo tiêu chuẩn IAEA.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn thừa nhận các giới hạn nghiên cứu:

  1. Giới hạn hình học đối tượng đo: Nghiên cứu mới tập trung giải quyết hình học bản phẳng và ống trụ tròn rỗng đường kính cố định ($273\text{ mm}$). Với các đường ống có bán kính cong thay đổi lớn hoặc cấu trúc đa lớp phức tạp (ví dụ: ống thép có lớp bọc cách nhiệt polyurethane hay lớp phủ composite), các hàm mật độ xác suất $S_0(x), D_0(x), P_0(x)$ cần được tái lập cấu hình.
  2. Ảnh hưởng của môi trường chứa bên trong: Các phép đo thực nghiệm được thực hiện trên ống thép rỗng (chứa không khí). Khi đường ống dẫn dầu thô, khí nén áp suất cao hoặc hỗn hợp đa pha lỏng - khí, sự tán xạ ngược từ lưu chất bên trong sẽ đóng góp thêm một thành phần tán xạ nền cần được phân tách.
  3. Hiệu ứng trôi đỉnh do nhiệt độ: Đầu dò nhấp nháy $\text{NaI(Tl)}$ và ống nhân quang điện (PMT) có độ nhạy nhiệt cao, có thể gây trôi kênh phổ (gain drift) khi hoạt động ngoài hiện trường nắng nóng nếu không có mạch ổn định phổ tự động.

Chương trình nghiên cứu tiếp theo (Future Research Agenda) gồm:

  • Mở rộng quy trình bán thực nghiệm cho các hệ đường ống có dải đường kính đa dạng ($50\text{ mm} - 1000\text{ mm}$) và các liên kết hàn phức tạp.
  • Tích hợp thuật toán tự động hiệu chỉnh trôi đỉnh phổ trực tiếp theo thời gian thực dựa trên đỉnh năng lượng chuẩn của nguồn hoặc bức xạ môi trường.
  • Ứng dụng quy trình phân tách phổ để phát hiện và dựng ảnh khuyết tật cục bộ, bọt khí và đo độ dày lớp cáu cặn bên trong đường ống dẫn dầu đang hoạt động (online non-invasive inspection).

Tác động và ảnh hưởng

Công trình nghiên cứu của TS. Võ Hoàng Nguyên tạo ra những tác động khoa học và kinh tế - xã hội rõ nét:

  • Tác động học thuật: Mở ra hướng tiếp cận mới trong việc xử lý tín hiệu hạt nhân bậc cao; các công thức giải tích về bề dày bão hòa ($T_s$) và quy trình bán thực nghiệm trở thành tài liệu tham khảo chuẩn mực cho các nghiên cứu sinh chuyên ngành Vật lý hạt nhân và Kỹ thuật hạt nhân.
  • Chuyển giao công nghệ công nghiệp: Cung cấp lõi thuật toán chính xác cho các doanh nghiệp chế tạo thiết bị NDT trong nước và quốc tế, trực tiếp phục vụ kiểm định an toàn lò hơi, bồn chứa áp lực, hệ thống đường ống dẫn dầu khí tại các nhà máy lọc hóa dầu (Dung Quất, Nghi Sơn) và giàn khoan ngoài khơi.
  • Hiệu quả kinh tế & An toàn xã hội: Giảm chi phí đầu tư thiết bị kiểm tra không phá hủy xuống nhiều lần so với hệ thống $\text{HPGe}$, ngăn ngừa nguy cơ sự cố rò rỉ hóa chất và nổ đường ống nhờ khả năng giám sát bề dày thành ống định kỳ với độ tin cậy cao.

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

  • Nghiên cứu sinh & Nhà khoa học Vật lý Hạt nhân: Tiếp cận phương pháp kết hợp mô phỏng GEANT4/MCNP6 với giải tích thống kê Gauss để giải quyết bài toán tán xạ ngược.
  • Kỹ sư NDT & Chuyên gia Kiểm định Chất lượng: Nắm vững công thức tính bề dày bão hòa để thiết lập giới hạn đo đạc chính xác tại hiện trường, lựa chọn cấu hình chuẩn trực tối ưu.
  • Các Viện nghiên cứu & Doanh nghiệp Chế tạo Thiết bị Hạt nhân: Sở hữu quy trình thuật toán khả thi để thương mại hóa các dòng máy đo bề dày không tiếp xúc sử dụng đầu dò $\text{NaI(Tl)}$ giá thành cạnh tranh.
  • Nhà quản lý An toàn Bức xạ: Có cơ sở khoa học vững chắc để cấp phép và áp dụng các quy chuẩn an toàn đối với hệ đo phóng xạ hoạt độ thấp trong công nghiệp.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc thiết lập công thức giải tích xác định bề dày bão hòa ($T_s$) dựa trên bản chất thống kê phân bố Gauss ở độ tin cậy 95%: $$T_s = -\frac{1}{\mu_{\text{eff}}} \ln\left[1 - \left(1 - \frac{1{,}96}{\sqrt{I_s}}\right)^2\right]$$ Công thức này đã giải quyết triệt để hạn chế kéo dài nhiều thập kỷ trong các nghiên cứu thực nghiệm (như Paramesh 1983; Singh 2006; Kiran 2016), chuyển đổi việc xác định $T_s$ từ ước lượng cảm tính sang tính toán định lượng chính xác tuyệt đối.

2. Điểm cải tiến phương pháp luận so với các nghiên cứu trước đây là gì?

So với nghiên cứu của Tạ Duy Tám và cộng sự (2017) vốn chỉ dùng một đa thức bậc bốn và một hàm Gauss đơn giản, luận án đã sử dụng GEANT4 để bóc tách chính xác bản chất vật lý của tán xạ hai lần thành mô hình ba hàm Gauss $[G_1(x) + G_2(x) + G_3(x)]$ (tương ứng ba nhóm cơ chế Compton-Compton, Compton-Rayleigh, Rayleigh-Compton) và mô hình hóa tán xạ trên hai lần bằng hàm phân bố bất đối xứng $P(x)$, giúp loại bỏ triệt để phông nhiễu mà không làm méo đỉnh tán xạ một lần.

3. Phát hiện thực nghiệm nào bất ngờ nhất và được chứng minh bằng dữ liệu ra sao?

Phát hiện bất ngờ nhất là tính bất biến của phổ mật độ xác suất tán xạ chuẩn hóa trên các ống thép C45 có độ dày thay đổi rộng từ $2{,}0\text{ mm}$ đến $22{,}0\text{ mm}$. Dù độ dày thành ống tăng gấp 11 lần, dạng hình học của các hàm mật độ xác suất $S_0(x), D_0(x), P_0(x)$ không hề biến dạng mà giữ nguyên hình dạng chuẩn, cho phép cố định các tham số phân bố và chỉ cần tối ưu hóa ba trọng số biên độ $(h_1, h_2, h_3)$ trong quy trình bán thực nghiệm.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Luận án cung cấp một giao thức thực nghiệm và mô phỏng cực kỳ chi tiết, bao gồm:

  • Bản vẽ kích thước cơ khí chính xác của hộp nguồn, thanh gắn nguồn chì, ống chuẩn trực đầu dò ($\varnothing 10\text{ mm}$ và $\varnothing 20\text{ mm}$).
  • Bộ tham số mở rộng năng lượng Gauss FT8 GEB trong MCNP6 ($a = -0{,}0036\text{ MeV}$, $b = 0{,}0621\text{ MeV}^{1/2}$, $c = 0{,}3146\text{ MeV}^{-1}$).
  • Toàn bộ tham số số học của các hàm Gauss thành phần và quy trình từng bước của thuật toán bán thực nghiệm để người đọc có thể lập trình tái lập 100%.

5. Tầm nhìn nghiên cứu 10 năm tới được vạch ra như thế nào?

Tầm nhìn 10 năm định hướng phát triển hệ thống chụp ảnh cắt lớp gamma tán xạ đơn diện (Single-Sided Compton Scatter Tomography) thời gian thực, tích hợp mảng đầu dò $\text{NaI(Tl)}$ đa kênh và trí tuệ nhân tạo (Edge AI) để quét tự động toàn diện các hệ thống đường ống ngầm và công trình biển mà không cần ngừng vận hành hệ thống.

Kết luận

Luận án tiến sĩ của Võ Hoàng Nguyên là một công trình khoa học mẫu mực, kết hợp nhuần nhuyễn giữa lý thuyết tán xạ bức xạ, mô phỏng Monte Carlo đa cấp và thực nghiệm đo đạc chính xác cao. Những kết luận khoa học cốt lõi bao gồm:

  1. Thiết lập thành công công thức giải tích thống kê xác định bề dày bão hòa ($T_s$) và bề dày thành vật liệu ($T$) với đầy đủ công thức lan truyền sai số.
  2. Khám phá bản chất phân bố năng lượng vi mô của các bậc tán xạ bằng GEANT4, đề xuất mô hình ba hàm Gauss cho tán xạ hai lần và hàm phân bố liên tục cho tán xạ trên hai lần.
  3. Chuẩn hóa Quy trình Bán thực nghiệm ba thành phần, giải quyết trọn vẹn bài toán tách chồng phổ gamma tán xạ trên đầu dò nhấp nháy $\text{NaI(Tl)}$.
  4. Kiểm chứng thành công trên bia phẳng (nhôm 5052, thép C45) và ống thép C45 ($\varnothing 273\text{ mm}$) với sai số xác định bề dày thực tế luôn duy trì dưới $3{,}5%$.
  5. Mở ra giải pháp chế tạo thiết bị NDT di động, hoạt độ thấp ($5\text{ mCi}$ $^{137}\text{Cs}$), chi phí tối ưu, bảo đảm an toàn bức xạ và đáp ứng hoàn hảo các yêu cầu kiểm định công nghiệp khắt khe.