Tổng quan luận án

Kiểm tra không phá hủy mẫu (Non-Destructive Testing – NDT) đóng vai trò then chốt trong việc đánh giá khuyết tật, độ ăn mòn và xác định bề dày kết cấu công nghiệp mà không làm ảnh hưởng đến tính toàn vẹn của đối tượng. Trong số các kỹ thuật NDT, phương pháp gamma tán xạ sở hữu đặc điểm kỹ thuật: nguồn phát bức xạ gamma và đầu dò ghi nhận được bố trí về cùng một phía so với đối tượng đo, đồng thời thiết bị không yêu cầu tiếp xúc trực tiếp với bề mặt mẫu. Đặc tính này mang lại ưu thế khi kiểm tra các đối tượng vận hành trong môi trường khắc nghiệt (nhiệt độ cao, áp suất lớn) hoặc các kết cấu chỉ có thể tiếp cận từ một hướng duy nhất như thân vỏ tàu thủy, thành lò nung, đường ống dẫn dầu và khí đốt. Đề tài bắt nguồn từ việc tiếp nối và mở rộng kết quả nghiên cứu đề tài cấp Bộ Công Thương (2012–2014) về hệ đo độ ăn mòn thành lò chịu lửa bằng phương pháp gamma tán xạ không tiếp xúc do Bộ môn Vật lý Hạt nhân – Kỹ thuật Hạt nhân, Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia TP. Hồ Chí Minh thực hiện.

Trong các ứng dụng hiện trường đòi hỏi tính cơ động cao, đầu dò nhấp nháy NaI(Tl) thường được ưu tiên lựa chọn nhờ ưu điểm không cần làm lạnh bằng nitơ lỏng và có hiệu suất ghi nhận bức xạ cao. Tuy nhiên, đầu dò NaI(Tl) có độ phân giải năng lượng thấp hơn đáng kể so với đầu dò bán dẫn siêu tinh khiết (HPGe). Mặt khác, để bảo đảm an toàn bức xạ và thuận tiện khi vận chuyển tại hiện trường, hệ đo thường sử dụng nguồn phóng xạ có hoạt độ thấp; điều này đòi hỏi không thể thu hẹp kích thước ống chuẩn trực quá mức nhằm duy trì tốc độ đếm thống kê trong thời gian đo cho phép. Việc mở rộng trường nhìn dẫn đến sự gia tăng mạnh mẽ của thành phần tán xạ nhiều lần (thành phần gây nhiễu). Sự chồng chập phân bố năng lượng giữa thành phần tán xạ một lần và thành phần tán xạ nhiều lần trong phổ đo NaI(Tl) làm suy giảm độ chính xác của phép xác định bề dày vật liệu. Luận án lựa chọn giải quyết khoảng trống này thông qua các mục tiêu nghiên cứu cụ thể:

  1. Nghiên cứu đặc trưng phân bố năng lượng của các thành phần tán xạ (tán xạ một lần, hai lần và trên hai lần) trong phổ gamma tán xạ toàn phần.
  2. Xây dựng các hàm toán học đại diện và phát triển phương pháp phân tách phổ gamma tán xạ thành các thành phần riêng biệt.
  3. Thiết lập quy trình bán thực nghiệm hoàn chỉnh để xử lý phổ đo từ đầu dò NaI(Tl), phục vụ tính toán bề dày vật liệu và bề dày bão hòa với các tiêu chuẩn: đơn giản, chính xác và ổn định.

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

  • Đối tượng nghiên cứu: Phổ bức xạ gamma tán xạ ghi nhận từ nguồn phóng xạ Cs-137 trên các bia tán xạ kim loại có hình dạng tấm phẳng (hợp kim nhôm 5052, thép C45) và bia dạng trụ rỗng (ống thép C45).
  • Phạm vi không gian và thời gian: Nghiên cứu được thực hiệ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, hoàn thành và bảo vệ năm 2023. Phạm vi kỹ thuật bao gồm hệ đo thực nghiệm gamma tán xạ góc $120^\circ \pm 0{,}1^\circ$, kết hợp mô phỏng Monte Carlo bằng hai chương trình chuyên dụng GEANT4 và MCNP6.

Tổng quan tài liệu và vị trí của luận án

Tổng quan tài liệu trong luận án hệ thống hóa các công trình nghiên cứu trong và ngoài nước theo ba hướng tiếp cận chính:

1. Các nghiên cứu giải tích và mô hình hóa tán xạ nhiều lần trên thế giới

  • Dumond (1930) sử dụng lý thuyết tán xạ cổ điển Thompson kết hợp hiệu chỉnh Breit để tính toán giải tích cường độ và phân bố phổ tán xạ hai lần theo góc tán xạ trên mẫu khối cầu.
  • Kirkpatrick (1937) công bố các tính toán giải tích cường độ tán xạ hai lần tại góc $90^\circ$ cho mẫu dạng cầu và trụ với photon tới có trạng thái phân cực bất kỳ.
  • Fernández (1991) giải phương trình vận chuyển Boltzmann để suy ra cường độ tán xạ hai lần của các hiệu ứng Compton và Rayleigh (gồm 4 cơ chế: Compton-Compton, Compton-Rayleigh, Rayleigh-Compton, Rayleigh-Rayleigh) cho mẫu dày vô hạn bị chiếu bởi chùm tia gamma đơn năng. Sau đó, Fernández và cộng sự (1999) mở rộng phương trình vectơ Boltzmann cho chùm photon phân cực.
  • Kakutani và cộng sự (2004) mô phỏng Monte Carlo phổ năng lượng tán xạ Compton hai lần từ tia X phân cực tuyến tính hoặc elip trên các bia có số nguyên tử $Z$ khác nhau.

2. Nghiên cứu thực nghiệm và mô phỏng Monte Carlo khảo sát đặc trưng tán xạ

  • Felsteiner và cộng sự (1975) dùng Monte Carlo khảo sát ảnh hưởng của bề dày mẫu, vật liệu và quãng đường tự do trung bình của photon lên số lượng tán xạ nhiều lần.
  • Halonen và cộng sự (1979) so sánh tiết diện Thompson phi tương đối, tiết diện tương đối tính Klein-Nishina và hiệu chỉnh động lượng electron Ribberfors trong tính toán tán xạ hai lần trên mẫu nhôm và niken hình trụ.
  • Tarim và cộng sự (2013) xây dựng chương trình Monte Carlo phân tách riêng biệt phổ năng lượng từ tán xạ một lần đến chín lần, làm rõ mức độ chồng chập của các bậc tán xạ trong phổ tổng.
  • Paramesh và cộng sự (1983) khảo sát sự phụ thuộc của bề dày bão hòa vào số nguyên tử $Z$ với nguồn Cs-137 (662 keV) ở góc $120^\circ$, thu được bề dày bão hòa của nhôm, sắt, đồng và chì lần lượt là 8,3 cm; 2,65 cm; 2,2 cm và 0,4 cm.
  • Singh và cộng sự (2006, 2007) khảo sát sự biến thiên của tán xạ nhiều lần và bề dày bão hòa theo năng lượng tới (279, 662, 1120 keV), đường kính ống chuẩn trực đầu dò và góc tán xạ ($50^\circ$ đến $130^\circ$), ghi nhận cường độ tán xạ nhiều lần đạt giá trị tối thiểu tại góc $90^\circ$.
  • Kiran và cộng sự (2016) kết hợp mô phỏng MCNP và thực nghiệm trên các bia C, Al, Fe, Cu ở các góc từ $40^\circ$ đến $180^\circ$.

3. Ứng dụng phương pháp gamma tán xạ trong kiểm tra không phá hủy

  • Hussein và cộng sự (1989) dùng chương trình SIMPHO và hệ đo Cs-137/NaI(Tl) để kiểm tra thanh sắt và lỗ rỗng trong khối bê tông.
  • Anjos và cộng sự (1989) phát hiện vết nứt bề mặt khối nhôm có kích thước trên 1,6 mm.
  • Mullin và cộng sự (1994) kết hợp nguồn Co-60, đầu dò HPGe và mô phỏng Monte Carlo để phát hiện khuyết tật bên trong nhôm có kích thước từ 0,5 mm.
  • Silva và cộng sự (1999) xây dựng đường chuẩn năng lượng đỉnh tán xạ để đo bề dày tấm nhôm (sai số dưới 4%) và phát hiện khuyết tật $\ge 2\text{ mm}$ trên ống sắt dày 6,3 mm.
  • Priyada và cộng sự (2011, 2012) ứng dụng đầu dò HPGe đánh giá độ ăn mòn thép mềm (sai số tương đương chụp ảnh cắt lớp và phóng xạ), phát hiện mặt phân cách lỏng–lỏng, lỏng–khí và lỗ rỗng bê tông.
  • Salgado và cộng sự (2021) kết hợp mô phỏng MCNP6 và mạng nơ-ron nhân tạo với hệ 3 đầu dò NaI(Tl) để dự đoán bề dày lớp cặn $\text{BaSO}_4$ trong đường ống (hơn 90% kết quả có sai số $< 10%$).

4. Tình hình nghiên cứu trong nước và khoảng trống luận án giải quyết

Tại Việt Nam, các tác giả T. Nghiệp (2005), Loan và cộng sự (2011), Phương và cộng sự (2012), Tuấn và cộng sự (2014), Nghiệp và cộng sự (2014), Kiên và cộng sự (2016), Chương và cộng sự (2017, 2020), Sang và cộng sự (2021) đã thực hiện nhiều khảo sát thực nghiệm và mô phỏng MCNP về tán xạ gamma nhằm xác định mật độ chất lỏng, bề dày tấm phẳng và khuyết tật vật liệu. Đặc biệt, Tam và cộng sự (2016) từng đề xuất mô hình làm khớp phổ tán xạ trên thép C45 bằng ba hàm: 1 hàm Gauss cho tán xạ một lần, 1 hàm Gauss cho tán xạ hai lần và 1 đa thức bậc bốn cho tán xạ trên hai lần.

Tuy nhiên, luận án chỉ ra hai hạn chế lớn trong các công trình đi trước:

  1. Việc mô tả thành phần tán xạ hai lần và trên hai lần bằng 1 hàm Gauss và 1 đa thức bậc bốn chưa phản ánh đúng bản chất phân bố năng lượng phi đối xứng và phức tạp của các bậc tán xạ bậc cao.
  2. Khái niệm "bề dày bão hòa" ($T_s$) được sử dụng phổ biến nhưng chưa từng có công thức giải tích mang tính định lượng xác thực dựa trên lý thuyết thống kê sai số đo. Luận án đã tập trung giải quyết trọn vẹn hai khoảng trống cốt lõi này.

Cơ sở lý thuyết và phương pháp nghiên cứu

1. Cơ sở lý thuyết tán xạ gamma và dẫn xuất công thức giải tích

Khi photon gamma năng lượng $E$ va chạm với electron trong vật liệu bia, tán xạ Compton làm photon đổi hướng theo góc $\theta$ và giảm năng lượng xuống $E'$, tuân theo hệ thức:

$$E' = \frac{E}{1 + \frac{E}{m_e c^2}(1 - \cos\theta)}$$

Tiết diện tán xạ vi phân Compton được xác định qua công thức Klein–Nishina:

$$\frac{d\sigma}{d\Omega} = \frac{r_0^2}{2} \left[\frac{E'}{E}\right]^2 \left[\frac{E'}{E} + \frac{E}{E'} - \sin^2\theta\right] = \frac{r_0^2}{2 \left[1 + \frac{E}{m_e c^2}(1 - \cos\theta)\right]^2} \left[1 + \cos^2\theta + \frac{\frac{E^2}{(m_e c^2)^2}(1 - \cos\theta)^2}{1 + \frac{E}{m_e c^2}(1 - \cos\theta)}\right]$$

Quá trình tán xạ một lần được mô hình hóa qua ba giai đoạn: (1) photon truyền từ nguồn đến điểm tương tác $P$ theo quãng đường $x$; (2) xảy ra tán xạ Compton vi phân tại $P$ trong thể tích $dV$; (3) photon tán xạ năng lượng $E'$ di chuyển quãng đường $x'$ thoát khỏi bia đến đầu dò. Tích phân trên toàn bộ bề dày $T$ của bia dẫn đến phương trình đường cong bão hòa của cường độ tán xạ một lần $I$:

$$I = I_s \left(1 - e^{-\mu_{eff} T}\right)$$

Trong đó:

  • $\mu_{eff} = \mu(E)\sec\theta_1 + \mu(E')\sec\theta_2$ là hệ số suy giảm hiệu dụng tuyến tính.
  • $I_s = \frac{k I_0}{\mu_{eff}}$ là cường độ tán xạ một lần tại bề dày vô hạn.

Từ phương trình trên, bề dày $T$ của bia tán xạ và sai số lan truyền $u_T$ được xác định theo công thức:

$$T = -\frac{1}{\mu_{eff}} \ln\left(1 - \frac{I}{I_s}\right)$$

$$u_T = \sqrt{\left(\frac{\ln(1 - I/I_s)}{\mu_{eff}^2} u_{\mu_{eff}}\right)^2 + \left(\frac{I}{\mu_{eff} I_s^2 (1 - I/I_s)} u_{I_s}\right)^2 + \left(\frac{1}{\mu_{eff} I_s (1 - I/I_s)} u_I\right)^2}$$

Về bề dày bão hòa, tác giả định nghĩa $T_s$ là bề dày nhỏ nhất mà kể từ đó sự gia tăng cường độ tán xạ không thể phân biệt được với thăng giáng thống kê của phép đo ở mức tin cậy 95% ($n = 1{,}96$):

$$I(T_s) + 1{,}96\sqrt{I(T_s)} \ge I_s - 1{,}96\sqrt{I_s}$$

Giải bất phương trình trên, tác giả thu được công thức tính bề dày bão hòa $T_s$ và sai số tương ứng $u_{T_s}$:

$$T_s = -\frac{1}{\mu_{eff}} \ln\left[1 - \left(1 - \frac{1{,}96}{\sqrt{I_s}}\right)^2\right]$$

$$u_{T_s} = \sqrt{\left(\frac{\ln\left[1 - \left(1 - \frac{1{,}96}{\sqrt{I_s}}\right)^2\right]}{\mu_{eff}^2} u_{\mu_{eff}}\right)^2 + \left(\frac{1{,}96 \left(1 - \frac{1{,}96}{\sqrt{I_s}}\right)}{\mu_{eff} I_s^{3/2} \left[1 - \left(1 - \frac{1{,}96}{\sqrt{I_s}}\right)^2\right]} u_{I_s}\right)^2}$$

2. Phương pháp thực nghiệm và công cụ mô phỏng Monte Carlo

Tham số / Thiết bị Chi tiết kỹ thuật
Nguồn phóng xạ Nguồn kín Cs-137 hãng Eckert & Ziegler sản xuất, hoạt độ 5 mCi (tính tại ngày 18/04/2013), phát tia gamma chính 661,657 keV (hiệu suất 84,99%) và tia X (31,8–37,3 keV). Vỏ bọc thép không gỉ hai lớp hình trụ đường kính 6,0 mm, chiều dài 8,0 mm.
Đầu dò bức xạ Đầu dò nhấp nháy $\text{NaI(Tl)}$ model 802 (Canberra), kích thước tinh thể hình trụ $76{,}2\text{ mm} \times 76{,}2\text{ mm}$.
Khối điện tử & Phần mềm Khối xử lý tích hợp Osprey (chứa cao thế, tiền khuếch đại, khuếch đại số và bộ MCA); phần mềm điều khiển và thu thập phổ Genie 2000 phiên bản 3.3.
Hệ chuẩn trực & Che chắn Khối chì che chắn nguồn (đường kính ngoài 130 mm, lỗ trong 10/20 mm); ống chuẩn trực đầu dò bằng chì kích thước $143 \times 112 \times 112\text{ mm}$, đường kính lỗ trong 84 mm, nắp mở lỗ chuẩn trực 10 mm (đo bia phẳng) và 20 mm (đo ống thép).
Bàn tán xạ cơ khí Khung thép không gỉ dày 8 mm, hai cánh tay đòn điều khiển bằng động cơ bước, thiết lập góc tán xạ cố định $120^\circ \pm 0{,}1^\circ$.
Mẫu bia thực nghiệm - Bia phẳng: Nhôm 5052 (dày 2,0–100,5 mm) và Thép C45 (dày 3,7–32,2 mm), kích thước $100 \times 300\text{ mm}$, thời gian đo 14.400 giây/mẫu.
- Bia trụ rỗng: Ống thép C45 dài $250 \pm 0{,}5\text{ mm}$, đường kính ngoài $273 \pm 0{,}1\text{ mm}$, bề dày thành ống 1,9–12,0 mm, thời gian đo 18.000 giây/mẫu.
Chương trình MCNP6 Mô phỏng phổ năng lượng F8:P (10 tỷ lịch sử hạt/phổ), sử dụng thư viện hạt nhân ENDF/B-VI, áp dụng thẻ mở rộng phổ FT8 GEB theo hàm $FWHM = a + b\sqrt{E + cE^2}$ với $a = -0{,}0036\text{ MeV}$, $b = 0{,}0621\text{ MeV}^{1/2}$, $c = 0{,}3146\text{ MeV}^{-1}$.
Chương trình GEANT4 Mô hình hóa chi tiết hình học hệ đo; kích hoạt thể tích nhạy (Sensitive Detector) tại vùng bia để ghi nhận chính xác lịch sử va chạm từng bậc (1 lần, 2 lần, trên 2 lần).

Thành phần hóa học của vật liệu bia sử dụng trong thực nghiệm và mô phỏng được quy chuẩn theo tiêu chuẩn kỹ thuật:

Vật liệu Thành phần nguyên tố định lượng
Hợp kim nhôm 5052 (TCVN 12513-7: 2018) $\text{Mg}: 2{,}200 - 2{,}500%;\ \text{Cr}: 0{,}150 - 0{,}350%;\ \text{Fe} < 0{,}400%;\ \text{Si} < 0{,}250%;\ \text{Cu} < 0{,}100%;\ \text{Mn} < 0{,}100%;\ \text{Zn} < 0{,}100%;\ \text{Al: phần còn lại.}$
Thép hợp kim C45 (TCVN 12142-1: 2017) $\text{C}: 0{,}420 - 0{,}500%;\ \text{Mn}: 0{,}500 - 0{,}800%;\ \text{Si}: 0{,}100 - 0{,}400%;\ \text{Cr} < 0{,}400%;\ \text{Ni} < 0{,}400%;\ \text{Cu} < 0{,}300%;\ \text{Mo} < 0{,}100%;\ \text{P} < 0{,}045%;\ \text{S} < 0{,}045%;\ \text{Fe: phần còn lại.}$

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

Chương 1: Nghiên cứu tổng quan

Chương 1 trình bày tổng quan tình hình nghiên cứu lý thuyết, thực nghiệm và mô phỏng về bức xạ gamma tán xạ trên thế giới và tại Việt Nam. Tác giả phân tích cơ chế tương tác Compton, Rayleigh, biểu thức vi phân Klein–Nishina và hiện tượng chồng chập phổ. Chương này thiết lập cơ sở toán học: thiết lập phương trình đường cong bão hòa $I = I_s (1 - e^{-\mu_{eff} T})$, dẫn xuất công thức giải tích xác định bề dày $T$ của bia tán xạ kèm công thức lan truyền sai số $u_T$, đồng thời đề xuất công thức xác định bề dày bão hòa $T_s$ và sai số $u_{T_s}$ dựa trên mức tin cậy thống kê Gauss 95%.

Chương 2: Bố trí thí nghiệm và các chương trình mô phỏng Monte Carlo

Chương 2 mô tả chi tiết thiết kế phần cứng của hệ đo gamma tán xạ:

  • Cấu trúc thanh gắn nguồn bằng chì, hộp chứa nguồn chì hai tầng trụ rỗng (đường kính ngoài 130 mm), khối chuẩn trực đầu dò NaI(Tl) kích thước $143 \times 112 \times 112\text{ mm}$ có nắp chuẩn trực tháo lắp linh hoạt (lỗ tròn 10 mm cho bia phẳng và 20 mm cho ống trụ).
  • Thiết lập bố trí hình học đo bia phẳng (khoảng cách nguồn–tâm bia 200,0 mm; tâm bia–đầu dò 180,0 mm; góc $120^\circ$) và bia trụ rỗng ống thép $\Phi 273\text{ mm}$ (khoảng cách nguồn–mặt ống 284,3 mm; mặt ống–đầu dò 203,3 mm).
  • Trình bày phương pháp mô phỏng Monte Carlo: ứng dụng MCNP6 với thẻ xử lý độ phân giải đỉnh FT8 GEB cùng bộ tham số thực nghiệm ($a = -0{,}0036$, $b = 0{,}0621$, $c = 0{,}3146$) để tái tạo phổ tổng; ứng dụng GEANT4 xây dựng thể tích nhạy (Sensitive Detector) tại bia để phân tách lịch sử photon.

Chương 3: Phát triển các phương pháp xử lý phổ gamma tán xạ

Chương 3 là trọng tâm học thuật của luận án, phát triển hai phương pháp xử lý phổ:

  1. Phương pháp 1 (Phân tách hai thành phần trên bia phẳng):

    • Khảo sát dạng phổ tán xạ một lần $s(x)$ và phổ tán xạ nhiều lần $m(x)$ từ mô phỏng trên các bia nhôm và thép phẳng.
    • Xây dựng thuật toán làm khớp phổ hai thành phần để bóc tách thành phần tán xạ một lần trực tiếp từ phổ đo thực nghiệm.

  2. Phương pháp 2 (Phân tách ba thành phần – Quy trình bán thực nghiệm cho ống thép):

    • Khảo sát chi tiết bằng GEANT4 trên các ống thép C45 có bề dày từ 2,0 mm đến 22,0 mm, chia phổ tán xạ toàn phần thành ba thành phần:
      • Thành phần tán xạ một lần: Được mô tả bằng hàm $S(x)$ xác định từ hàm mật độ xác suất chuẩn hóa.
      • Thành phần tán xạ hai lần: Phát hiện phân bố có dạng đa đỉnh bất đối xứng, được làm khớp chính xác bằng tổng của ba hàm Gauss: $D(x) = G_1(x) + G_2(x) + G_3(x)$.
      • Thành phần tán xạ trên hai lần: Được làm khớp bằng hàm phân bố $M(x) = P(x)$.
    • Xây dựng Quy trình bán thực nghiệm 3 bước:
      • Bước 1: Tính toán các hàm mật độ xác suất trung bình $S(x)$, $D(x)$, $M(x)$ từ mô phỏng Monte Carlo.
      • Bước 2: Làm khớp phổ thực nghiệm bằng tổ hợp tuyến tính $F(x) = h_1 S(x) + h_2 D(x) + h_3 M(x)$ để trích xuất các trọng số $h_1, h_2, h_3$.
      • Bước 3: Thu nhận diện tích đỉnh tán xạ một lần $I_1$ phục vụ xác định bề dày vật liệu.

Chương 4: Kết quả áp dụng các phương pháp xử lý phổ gamma tán xạ

Chương 4 kiểm chứng tính chính xác của các phương pháp đề xuất:

  • Bia nhôm và thép dạng tấm phẳng: Làm khớp đường cong bão hòa thực nghiệm với nguồn Cs-137, xác định chính xác các hệ số $I_s$ và $\mu_{eff}$. Kết quả tính toán bề dày cho các tấm nhôm (2,0–100,5 mm) và tấm thép (3,7–32,2 mm) cho thấy sự tương quan chặt chẽ với số liệu đo thực tế bằng thước kẹp điện tử.
  • Ống thép C45 dạng trụ rỗng:
    • Áp dụng quy trình bán thực nghiệm trên tập phổ mô phỏng MCNP6 (bề dày ống từ 2,0 mm đến 22,0 mm) và phổ đo thực nghiệm (bề dày thành ống từ 1,9 mm đến 12,0 mm).
    • Trích xuất thành công các hệ số $h_1, h_2, h_3$. Bề dày thành ống tính toán từ diện tích tán xạ một lần hoàn toàn phù hợp với bề dày thực tế trong phạm vi sai số cho phép.

Chương 5: Kết luận và kiến nghị

Tổng kết toàn bộ các phát hiện lý thuyết và thực nghiệm, khẳng định các điểm mới của luận án và đưa ra 3 kiến nghị định hướng cho các nghiên cứu tiếp nối.

Kết quả và những đóng góp mới

1. Đóng góp mới về mặt lý luận và khoa học

  • Xây dựng công thức giải tích tính bề dày bão hòa $T_s$: Luận án lần đầu tiên thiết lập công thức toán học tường minh xác định bề dày bão hòa trong phép đo gamma tán xạ dựa trên lý thuyết thăng giáng thống kê Gauss ở độ tin cậy 95%. Công thức này cung cấp cơ sở định lượng để xác định giới hạn đo bề dày cực đại của một cấu hình hệ đo cụ thể.
  • Làm rõ bản chất phân bố năng lượng tán xạ bậc cao: Sử dụng GEANT4 chứng minh phổ tán xạ hai lần có cấu trúc phân bố phức hợp và đề xuất mô hình toán học tổ hợp 3 hàm Gauss $G_1(x) + G_2(x) + G_3(x)$ để mô tả thành phần này, khắc phục nhược điểm mô hình hóa bằng đa thức bậc bốn đơn thuần trong các nghiên cứu trước đây.

2. Đóng góp mới về mặt thực tiễn

  • Quy trình bán thực nghiệm xử lý phổ ba thành phần: Xây dựng thành công quy trình chuẩn hóa kết hợp giữa mô phỏng Monte Carlo và đo đạc thực nghiệm, cho phép phân tách tự động các thành phần tán xạ trên phổ ghi nhận bằng đầu dò nhấp nháy NaI(Tl).
  • Giải pháp NDT hiện trường tin cậy: Cung cấp giải pháp đo đạc không tiếp xúc, cơ động cao, sử dụng nguồn phóng xạ Cs-137 hoạt độ thấp (5 mCi) kết hợp đầu dò NaI(Tl) hoạt động ở nhiệt độ phòng, đo chính xác bề dày tấm phẳng kim loại và thành ống thép công nghiệp.

3. Kiến nghị và giải pháp đề xuất

  • Tác giả kiến nghị ứng dụng quy trình bán thực nghiệm để phát triển các phần mềm tự động hóa quá trình phân tích số liệu trên các thiết bị đo NDT hiện trường.
  • Đề xuất áp dụng kỹ thuật cho việc kiểm tra chất lượng định kỳ các đường ống dẫn xăng dầu, khí đốt và thành bồn bể công nghiệp đang vận hành mà không cần dừng hệ thống.

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

Theo nội dung trình bày trong luận án, các hạn chế và hướng nghiên cứu tiếp theo bao gồm:

Hạn chế của nghiên cứu:

  • Hệ đo thực nghiệm mới khảo sát cố định tại góc tán xạ $120^\circ \pm 0{,}1^\circ$ với một loại nguồn phóng xạ Cs-137 năng lượng 662 keV và đầu dò NaI(Tl) kích thước $76{,}2 \times 76{,}2\text{ mm}$.
  • Đối tượng kiểm chứng thực nghiệm tập trung vào các kim loại đồng nhất (hợp kim nhôm 5052 và thép C45) dạng tấm phẳng và ống trụ đơn lớp có đường kính ngoài cố định 273 mm.
  • Chưa khảo sát thực nghiệm trong điều kiện bên trong ống có chứa lưu chất chuyển động đa pha (dầu–khí–nước) hoặc môi trường nhiệt độ thực tế của lò nung.

Hướng nghiên cứu tiếp theo:

  1. Mở rộng khảo sát đặc trưng phân bố của các thành phần tán xạ ở các góc đo khác nhau và với các loại nguồn bức xạ gamma đơn năng khác.
  2. Ứng dụng và kiểm chứng quy trình bán thực nghiệm trên các đối tượng hình học phức tạp hơn (như ống cong, mối hàn, cấu trúc nhiều lớp kim loại - lớp cách nhiệt).
  3. Nghiên cứu tích hợp các thuật toán học máy hoặc mạng nơ-ron nhân tạo vào quy trình xử lý phổ nhằm nâng cao khả năng phân tách phổ tự động trong thời gian thực.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu trong các lĩnh vực:

  • Nghiên cứu sinh và Giảng viên ngành Vật lý hạt nhân, Kỹ thuật hạt nhân: Có thể tham khảo phương pháp luận kết hợp song song giữa hai phần mềm Monte Carlo (MCNP6 và GEANT4) để phân tách các cơ chế tương tác vi mô của bức xạ với vật chất, cũng như kỹ thuật xây dựng hàm toán học mô hình hóa phổ nhấp nháy.
  • Kỹ sư và Chuyên gia Kiểm tra không phá hủy (NDT): Tham khảo bản vẽ thiết kế cơ khí hệ chuẩn trực chì, bảng thông số hình học tối ưu cho hệ đo tán xạ góc $120^\circ$, và quy trình toán học bán thực nghiệm để phát triển các hệ đo bề dày ống và bồn bể công nghiệp bằng nguồn phóng xạ hoạt độ thấp.
  • Cơ quan Quản lý an toàn bức xạ và Doanh nghiệp công nghiệp: Tham khảo giải pháp kỹ thuật thay thế việc sử dụng nguồn phóng xạ hoạt độ cao bằng nguồn hoạt độ thấp (5 mCi) kết hợp thuật toán xử lý phổ, nâng cao mức độ an toàn bức xạ tại hiện trường.

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

1. Tại sao luận án lựa chọn đầu dò NaI(Tl) thay vì đầu dò HPGe cho hệ đo hiện trường?

Đầu dò HPGe có độ phân giải năng lượng cao nhưng đòi hỏi phải làm lạnh liên tục bằng nitơ lỏng hoặc hệ thống làm lạnh điện tử cồng kềnh, gây khó khăn lớn khi di chuyển tại hiện trường. Đầu dò nhấp nháy NaI(Tl) hoạt động ổn định ở nhiệt độ môi trường, có hiệu suất ghi nhận bức xạ cao và kích thước nhỏ gọn, đáp ứng yêu cầu về tính cơ động trong kiểm tra công nghiệp.

2. Nguyên lý xây dựng công thức tính bề dày bão hòa $T_s$ trong luận án là gì?

Bề dày bão hòa $T_s$ được xác định dựa trên tiêu chuẩn thống kê: là bề dày nhỏ nhất mà tại đó tốc độ gia tăng cường độ tán xạ một lần không còn phân biệt được với thăng giáng thống kê của phép đo (tuân theo phân bố Gauss) ở mức độ tin cậy 95% ($n = 1{,}96$). Từ bất phương trình $I(T_s) + 1{,}96\sqrt{I(T_s)} \ge I_s - 1{,}96\sqrt{I_s}$, tác giả giải ra công thức giải tích $T_s = -\frac{1}{\mu_{eff}} \ln\left[1 - \left(1 - \frac{1{,}96}{\sqrt{I_s}}\right)^2\right]$.

3. Thành phần tán xạ hai lần được mô hình hóa bằng dạng hàm toán học nào trong Phương pháp 2?

Qua mô phỏng chi tiết bằng GEANT4, tác giả phát hiện phổ tán xạ hai lần có cấu trúc đa đỉnh bất đối xứng và chứng minh thành phần này được làm khớp chính xác nhất bằng tổ hợp tuyến tính của ba hàm phân bố Gauss: $D(x) = G_1(x) + G_2(x) + G_3(x)$, thay vì chỉ dùng một hàm Gauss hoặc đa thức bậc bốn như các nghiên cứu trước đây.

4. Cấu hình nguồn phóng xạ và hình học đo thực nghiệm trong luận án được thiết lập như thế nào?

Luận án sử dụng nguồn phóng xạ kín Cs-137 hoạt độ 5 mCi (ngày 18/04/2013) do Eckert & Ziegler sản xuất. Góc tán xạ giữa trục nguồn và trục đầu dò được cố định ở $120^\circ \pm 0{,}1^\circ$. Đối với bia phẳng, khoảng cách nguồn–bia là 200,0 mm và bia–đầu dò là 180,0 mm. Đối với ống thép đường kính 273 mm, khoảng cách nguồn–mặt ngoài ống là 284,3 mm và mặt ngoài ống–đầu dò là 203,3 mm.

5. Thẻ xử lý nào được sử dụng trong chương trình MCNP6 để tái tạo độ phân giải năng lượng của đầu dò NaI(Tl)?

Luận án sử dụng thẻ xử lý FT8 GEB (Gaussian Energy Broadening) trong MCNP6 với hàm phi tuyến phụ thuộc năng lượng $FWHM = a + b\sqrt{E + cE^2}$. Các tham số thực nghiệm được xác định cho đầu dò NaI(Tl) trong luận án là $a = -0{,}0036\text{ MeV}$, $b = 0{,}0621\text{ MeV}^{1/2}$, $c = 0{,}3146\text{ MeV}^{-1}$.

Kết luận

Luận án của nghiên cứu sinh Võ Hoàng Nguyên đã giải quyết bài toán phân tách thành phần tán xạ nhiều lần trong phổ gamma tán xạ ghi nhận bằng đầu dò nhấp nháy NaI(Tl). Tác giả đã thiết lập thành công công thức giải tích xác định bề dày bão hòa dựa trên cơ sở thống kê, đồng thời làm rõ cấu trúc phân bố của các bậc tán xạ bằng chương trình GEANT4. Quy trình bán thực nghiệm ba thành phần được đề xuất đã chứng minh độ chính xác và tính ổn định cao trong việc xác định bề dày tấm phẳng nhôm, thép và ống thép C45, đóng góp giải pháp kỹ thuật có giá trị ứng dụng thực tiễn cho lĩnh vực kiểm tra không phá hủy tại hiện trường.