Giới thiệu dự án
Trong bối cảnh công nghiệp hóa và hiện đại hóa, nhu cầu kiểm tra không phá hủy (Non-Destructive Testing - NDT) đối với các kết cấu kim loại, đường ống dẫn áp lực, bồn chứa hóa chất và lò hơi công nghiệp tăng trưởng hơn 8.5% mỗi năm trên toàn cầu. Các khuyết tật ngầm bên trong vật liệu như rỗ khí, nứt ngầm, hoặc sự ăn mòn làm mỏng thành kim loại là nguyên nhân hàng đầu dẫn đến các sự cố thảm họa nếu không được phát hiện kịp thời.
Phương pháp truyền thống như chụp ảnh phóng xạ truyền qua (Transmission Radiography) hay siêu âm kiểm tra đòi hỏi phải tiếp cận được từ hai phía của vật thể hoặc yêu cầu dừng toàn bộ dây chuyền sản xuất để tháo dỡ lớp bảo ôn, gây tổn thất kinh tế rất lớn. Để giải quyết triệt để rào cản này, phương pháp gamma tán xạ ngược (Gamma Backscattering) nổi lên như một giải pháp tối ưu cho phép bố trí cả nguồn phóng xạ và đầu dò ở cùng một phía, kiểm tra trực tiếp mà không cần dừng thiết bị. Đề tài "Xác định vị trí khuyết tật của vật liệu bằng phương pháp Monte-Carlo" tập trung nghiên cứu mô hình hóa và khảo sát định lượng khả năng phát hiện khuyết tật trên tấm thép chịu nhiệt C45 thông qua tương tác tán xạ Compton.
+-----------------------------+
| Nguồn bức xạ gamma Co-60 |
| (1173 & 1332 keV) |
+--------------+--------------+
| Chùm tia hẹp (Góc tới 90 độ)
v
+-----------------------------------------------------------+
| Tấm thép C45 (Chiều dày 25 mm) |
| [ Vùng khuyết tật: Lỗ rỗng d=2.2 - 3.0 cm ] |
+-----------------------------+-----------------------------+
| Tán xạ ngược (Góc 120 độ)
v
+-----------------------------+
| Đầu dò nhấp nháy NaI(Tl) |
| (3" x 3" Scintillator) |
+--------------+--------------+
| Xung quang phổ
v
+-----------------------------+
| Thuật toán khớp phổ |
| (1-Gauss vs 2-Gauss) |
+-----------------------------+
Mục tiêu dự án
- Thiết lập mô hình hình học và vật lý chi tiết của hệ đo gamma tán xạ ngược gồm nguồn phóng xạ $^{60}\text{Co}$, buồng chuẩn trực chì (Collimator), đầu dò nhấp nháy $\text{NaI(Tl)}$ ($76 \times 76\text{ mm}$) và mẫu thép C45 dày 25 mm bằng chương trình mô phỏng MCNP5 (Monte Carlo N-Particle v5).
- Mô phỏng quá trình truyền và tán xạ của $3 \times 10^9$ hạt photon gamma nhằm xác định quang phổ tán xạ tại góc tán xạ $120^\circ$.
- Xây dựng đường chuẩn năng lượng và khảo sát mối tương quan giữa số đếm đỉnh tán xạ đơn (Single Scattering) với bề dày còn lại của vật liệu ($5\text{ mm} \le d \le 23\text{ mm}$) và đường kính khuyết tật ($2.2\text{ cm} \le \Phi \le 3.0\text{ cm}$).
- Đánh giá hiệu quả của hai phương pháp khớp phổ: 1 đỉnh Gauss + nền đa thức so với 2 đỉnh Gauss (phân tách tán xạ đơn và tán xạ đôi) + nền đa thức bằng phần mềm COLEGRAM.
- Xác định giới hạn bề dày bão hòa và ngưỡng nhận biết vị trí khuyết tật hiệu dụng theo hệ số khuyết tật $C (%)$.
Phạm vi và giới hạn nghiên cứu
- Đối tượng: Hiện tượng gamma tán xạ ngược do hiệu ứng Compton chi phối.
- Mẫu khảo sát: Tấm thép C45 đồng nhất dày 25 mm, chứa các khuyết tật hình trụ dạng lỗ tròn có đường kính $\Phi \in {2.2\text{ cm}, 2.6\text{ cm}, 3.0\text{ cm}}$, độ sâu đáy lỗ biến thiên tương ứng bề dày kim loại còn lại từ 5 mm đến 23 mm (bước nhảy 2 mm).
- Hệ đo: Nguồn $^{60}\text{Co}$ bức xạ đa năng lượng (1.173 MeV và 1.332 MeV), góc chiếu vuông góc $90^\circ$, góc tán xạ thu nhận $120^\circ$.
Phân tích và thiết kế giải pháp
Phân tích hiện trạng
| Tiêu chí |
Phương pháp truyền qua (Transmission NDT) |
Phương pháp siêu âm (Ultrasonic Testing - UT) |
Phương pháp Gamma tán xạ ngược (Monte-Carlo Simulation) |
| Khả năng tiếp cận |
Bắt buộc tiếp cận từ 2 phía (Nguồn & Phim/Đầu dò) |
Tiếp cận 1 phía (cần tiếp xúc bề mặt trực tiếp & chất tiếp âm) |
Tiếp cận 1 phía, không cần tiếp xúc trực tiếp |
| Trạng thái vận hành |
Thường phải dừng máy, tháo lớp cách nhiệt |
Cần làm sạch bề mặt, khó đo ở nhiệt độ siêu cao |
Đo online khi hệ thống đang vận hành, xuyên lớp bảo ôn |
| Chi phí thử nghiệm |
Cao (chi phí an toàn bức xạ, thiết bị chụp) |
Trung bình (đòi hỏi thợ bậc cao) |
Tối ưu chi phí nhờ mô phỏng MCNP5 trước khi chế tạo hệ đo |
| Độ chính xác vị trí khuyết tật ngầm |
Cao đối với chiều sâu tổng thể |
Tốt với vết nứt phẳng, kém với rỗ khí đa hướng |
Định lượng rõ qua suy giảm số đếm đỉnh Compton |
MoSCoW PRIORITIZATION MATRIX
+-----------------------------------------------------------------------------------------+
| MUST HAVE (Bắt buộc) |
| - Mô phỏng chuẩn hình học đầu dò NaI(Tl) 76x76mm + vỏ nhôm + lớp phản xạ Al2O3. |
| - Sử dụng nguồn Co-60 với năng lượng chuẩn 1.173 MeV và 1.332 MeV. |
| - Thuật toán Tally 8 (Pulse Height Distribution) thu nhận phổ năng lượng thực tế. |
+-----------------------------------------------------------------------------------------+
| SHOULD HAVE (Nên có) |
| - Khảo sát đa kích thước khuyết tật (d = 2.2 cm, 2.6 cm, 3.0 cm). |
| - So sánh trực tiếp 2 mô hình làm khớp hàm (1-Gaussian vs 2-Gaussian) trên COLEGRAM. |
+-----------------------------------------------------------------------------------------+
| COULD HAVE (Có thể mở rộng) |
| - Mở rộng sang các nguồn đơn năng như Cs-137 (662 keV) hoặc Ir-192. |
| - Tự động hóa quá trình trích xuất Tally Output sang Origin qua Python script. |
+-----------------------------------------------------------------------------------------+
| WON'T HAVE THIS TIME (Chưa thực hiện) |
| - Thử nghiệm phá hủy thực tế với nguồn phóng xạ hoạt độ cao trên mẫu thép thực nghiệm. |
+-----------------------------------------------------------------------------------------+
Thiết kế hệ thống
Cơ sở vật lý tương tác bức xạ
Quá trình tán xạ Compton giữa photon gamma và electron tự do/liên kết yếu tuân theo phương trình bảo toàn năng lượng và xung lượng:
$$E' = \frac{E}{1 + \frac{E}{m_e c^2}(1 - \cos\theta)}$$
Trong đó:
- $E$: Năng lượng photon tới ban đầu ($E_{\text{avg}} = \frac{1173 + 1332}{2} = 1253\text{ keV}$ đối với $^{60}\text{Co}$).
- $m_e c^2 = 511\text{ keV}$: Năng lượng nghỉ của electron.
- $\theta = 120^\circ$: Góc tán xạ hình học của hệ đo.
- $E'$: Năng lượng tán xạ lý thuyết tính được:
$$E' = \frac{1253}{1 + \frac{1253}{511}(1 - \cos 120^\circ)} = \frac{1253}{1 + 2.452 \times 1.5} \approx 268\text{ keV}$$
Mật độ electron của vật liệu tán xạ $n_e$ được xác định thông qua khối lượng riêng $\rho$, số hiệu nguyên tử $Z$ và nguyên tử khối $A$:
$$n_e = \rho \cdot N_A \cdot \frac{Z}{A}$$
Khi xuất hiện khuyết tật rỗng bên trong thép C45, mật độ electron cục bộ giảm mạnh về $\approx 0$, làm suy giảm đột ngột cường độ chùm tia tán xạ ngược. Để định lượng mức độ khuyết tật, hệ số khuyết tật $C (%)$ được xác lập:
$$C = \frac{C_0 - C_i}{C_0} \times 100%$$
Với $C_0$ là diện tích đỉnh tán xạ của tấm thép nguyên vẹn (dày 25 mm) và $C_i$ là diện tích đỉnh tán xạ tại vị trí có khuyết tật.
MÔ HÌNH HÌNH HỌC HỆ ĐO TRÊN MCNP5
(Hệ tọa độ Descartes)
| y
|
| Collimator Chì
+-------+======+
| Nguồn | Co-60| ---> Chùm tia gamma tới
+-------+======+ (Trục Y = 0)
|
Góc tán xạ 120 độ /
/
/
v
+-------------------+
| Đầu dò NaI(Tl) | [Tinh thể NaI: Phi 76 x 76 mm]
| + Cửa sổ Si (2mm) | [Vỏ nhôm ngoài: Phi 83 mm]
| + Lớp Al2O3 (2mm) | [Mô phỏng PMT: Al đặc 30mm]
+-------------------+
|
+-----------------------------+
|
x = 0 v
+---------------------------------------------------------+
| Tấm thép C45: 25 mm |
| +-------------------+ |
| | Khuyết tật (Lỗ) | |
| +-------------------+ |
| |<-- Bề dày còn lại: 5 - 23 mm -->| |
+---------------------------------------------------------+
Technology Stack & Môi trường tính toán
- Mã nguồn mô phỏng Monte Carlo: MCNP5 Version 1.40 (Los Alamos National Laboratory).
- Thư viện dữ liệu hạt nhân: ENDF/B-VI Photoatomic Data Library (xử lý hiệu ứng quang điện, tán xạ kết hợp Coherent/Rayleigh, tán xạ không kết hợp Incoherent/Compton, và hiệu ứng tạo cặp).
- Công cụ xử lý và khớp phổ: COLEGRAM (Phân tích khớp đỉnh phi tuyến sử dụng hàm Gauss và đa thức nền), ORIGIN Lab 8.5.
- Ngôn ngữ xử lý số liệu tự động: Bash Script / Python 3.x trên nền tảng tính toán hiệu năng cao.
Methodology
+-------------------------------------------------------------------------------+
| GIAI ĐOẠN 1: MÔ HÌNH HÓA VẬT LÝ & HÌNH HỌC (MCNP5) |
| - Định nghĩa Surface Cards: Mặt phẳng (P, PX, PY, PZ), Mặt trụ (C/X, C/Y, CZ).|
| - Định nghĩa Cell Cards: Gán vật liệu NaI, Al, Si, Pb, Thép C45, Không khí. |
| - Thiết lập nguồn SDEF: Vị trí nguồn Co-60, phân bố hướng qua Collimator chì. |
+---------------------------------------+---------------------------------------+
|
v
+-------------------------------------------------------------------------------+
| GIAI ĐOẠN 2: THỰC THI MÔ PHỎNG MONTE CARLO N-PARTICLE |
| - Khởi tạo NPS = 3,000,000,000 hạt photon. |
| - Ghi nhận phổ Tally F8:P trên Cell tinh thể NaI kết hợp thẻ mở rộng Broadening|
| (FT8 GEB - Gaussian Energy Broadening) để tái tạo độ phân giải thực tế. |
+---------------------------------------+---------------------------------------+
|
v
+-------------------------------------------------------------------------------+
| GIAI ĐOẠN 3: XÂY DỰNG ĐƯỜNG CHUẨN NĂNG LƯỢNG |
| - Hiệu chuẩn vị trí kênh với các đỉnh năng lượng chuẩn đã biết. |
| - Thiết lập phương trình chuẩn: E(keV) = 0.20407 * Channel - 13.53648 |
+---------------------------------------+---------------------------------------+
|
v
+-------------------------------------------------------------------------------+
| GIAI ĐOẠN 4: PHÂN TÍCH QUANG PHỔ & ĐỊNH LƯỢNG KHUYẾT TẬT |
| - Trích xuất dữ liệu từ Tally 8 Output. |
| - Khớp phổ trên COLEGRAM: Đánh giá phương án 1-Gauss vs 2-Gauss. |
| - Tính diện tích đỉnh tán xạ Si, sai số delta S và hệ số khuyết tật C (%). |
+-------------------------------------------------------------------------------+
Implementation và kết quả
Development process
Cấu trúc file đầu vào MCNP5 (Input File Deck) được xây dựng chặt chẽ theo cấu trúc 3 khối: Cell Cards, Surface Cards, và Data Cards.
c =================================================================
c FILE INPUT MCNP5: MO PHONG GAMMA TAN XA NGUOC XAC DINH KHUYET TAT
c =================================================================
c --- KHOI 1: DINH NGHIA CELL CARDS ---
1 1 -3.67 -10 11 -12 imp:p=1 $ Tinh the NaI (Mat do 3.67 g/cm3)
2 2 -2.70 (-13 14 -15):#1 imp:p=1 $ Vo boc Nhom bao quanh dau do
3 3 -7.85 -20 21 -22 #30 imp:p=1 $ Tam thep C45 co chua lo khuyet tat
30 0 -31 32 -33 imp:p=1 $ Vung khuyet tat (Chankhong / Khongkhi)
4 4 -11.34 -40 41 -42 imp:p=1 $ Buong chi chuan truc (Collimator)
99 0 #1 #2 #3 #4 imp:p=1 $ Khong khi moi truong xung quanh
999 0 50 imp:p=0 $ Vung bien ngoai (Graveyard)
c --- KHOI 2: DINH NGHIA SURFACE CARDS ---
10 CZ 3.8 $ Mat tru tinh the NaI (Ban kinh R = 38 mm)
11 PZ 0.0 $ Mat phang day tinh the
12 PZ 7.6 $ Mat phang dinh tinh the (Dai L = 76 mm)
20 RPP -15.0 15.0 -15.0 15.0 0 2.5 $ Hop thep C45 kich thuoc 30x30x2.5 cm
31 C/Z 0.0 0.0 1.1 $ Lo khuyet tat hinh tru ban kinh R = 1.1 cm
50 SO 100.0 $ Mat cau gioi han toan he tinh toan
c --- KHOI 3: DINH NGHIA DATA CARDS ---
mode p
c Dinh nghia vat lieu thep C45 (Mat do 7.85 g/cm3)
m3 26000 0.98 6000 0.0045 25055 0.007 14000 0.003 15031 0.0004
c Dinh nghia nguon Co-60 (Hai vach nang luong 1.173 MeV va 1.332 MeV)
sdef pos=0 0 10.0 erg=d1 dir=1 vec=0 0 -1 par=p
si1 L 1.1732 1.3325
sp1 0.5 0.5
c Ghi nhan pho nang luong Tally 8 trong Cell 1 (Tinh the NaI)
f8:p 1
e8 0.01 2000i 1.5
ft8 GEB 0.005 0.05 0.001 $ Hieu chinh mo rong pho thuc te (GEB Card)
nps 3000000000
Testing và validation
1. Xây dựng đường chuẩn năng lượng
Để chuyển đổi từ số thứ tự kênh ghi nhận $Ch$ trên máy phân tích đa kênh (MCA) sang giá trị năng lượng thực $E (\text{keV})$, một tập hợp 6 điểm đo chuẩn từ các nguồn phóng xạ chuẩn đã biết trước năng lượng được khớp tuyến tính:
$$\text{Phương trình chuẩn hóa: } E (\text{keV}) = 0.20407 \times Ch - 13.53648$$
Năng lượng (keV)
^
1400 | * Co-60 (1332 keV)
1200 | * Co-60 (1173 keV)
1000 |
800 |
600 | * Cs-137 (662 keV)
400 | * Ba-133 (356 keV)
200 | * Am-241 (59.5 keV)
0 +---|---|---|---|---|---|---|---|---|---|---|---|---> Kênh (Channel)
0 1000 2000 3000 4000 5000 6000 7000
- Kiểm chứng đỉnh tán xạ Compton:
- Vị trí kênh đỉnh tán xạ một lần ghi nhận trên phổ: $Ch = 1355$.
- Năng lượng thực nghiệm suy ra: $E_{\text{sim}} = 0.20407 \times 1355 - 13.53648 = 263.02\text{ keV}$.
- Sai số tương đối so với tính toán lý thuyết tương đối tính ($E_{\text{theory}} = 268\text{ keV}$):
$$\text{Error} = \frac{|268 - 263.02|}{268} \times 100% = 1.85% \quad (\text{Độ tin cậy > 98.1%})$$
Kết quả đạt được
1. So sánh phương pháp làm khớp 1-Gauss và 2-Gauss
Toàn bộ quang phổ tán xạ sau khi thu thập được đưa vào phần mềm COLEGRAM để phân tích diện tích đỉnh tán xạ $S$ và sai số tương đối $\delta S (%)$.
| Bề dày còn lại $d$ (mm) |
Diện tích đỉnh (1-Gauss) $S_1$ |
Sai số $\delta S_1$ (%) |
Diện tích đỉnh (2-Gauss) $S_2$ |
Sai số $\delta S_2$ (%) |
Hệ số $C_{1\text{-Gauss}}$ (%) |
Hệ số $C_{2\text{-Gauss}}$ (%) |
| 5 |
114,832 |
0.224 |
112,795 |
0.231 |
47.96 |
48.28 |
| 7 |
139,120 |
0.228 |
137,014 |
0.230 |
36.96 |
37.17 |
| 9 |
165,490 |
0.228 |
163,737 |
0.229 |
25.01 |
24.92 |
| 11 |
185,210 |
0.230 |
183,110 |
0.232 |
16.07 |
16.03 |
| 13 |
199,091 |
0.231 |
197,450 |
0.233 |
9.78 |
9.46 |
| 15 |
208,410 |
0.232 |
206,120 |
0.235 |
5.56 |
5.48 |
| 17 |
214,500 |
0.233 |
212,300 |
0.236 |
2.80 |
2.65 |
| 21 |
219,080 |
0.235 |
216,900 |
0.238 |
0.72 |
0.54 |
| 25 (Không khuyết tật) |
220,678 |
0.235 |
218,085 |
0.238 |
0.00 |
0.00 |
Số đếm đỉnh tán xạ (Counts)
^
220k | *---*---* (Vùng 2: Bão hòa > 13mm)
200k | *---*
180k | *---*
160k | *---*
140k | *---*
120k | *---* (Vùng 1: Tuyến tính 5 - 9mm)
100k |
+-----+-----+-----+-----+-----+-----+-----+-----+-----> Bề dày còn lại (mm)
5 7 9 11 13 15 17 21 25
2. Quy luật biến thiên và các vùng đặc trưng
- Vùng nhạy tuyến tính ($5\text{ mm} \le d \le 9\text{ mm}$): Số đếm đỉnh tán xạ tăng mạnh theo độ dày còn lại của thép. Hệ số khuyết tật $C$ giảm nhanh từ $48.28%$ xuống $24.92%$, tạo tương phản cực lớn để xác định chính xác độ sâu đáy khuyết tật.
- Vùng chuyển tiếp ($9\text{ mm} < d \le 13\text{ mm}$): Tốc độ tăng diện tích đỉnh giảm dần, hệ số $C$ giảm từ $24.92%$ về $\approx 9.46%$.
- Vùng bão hòa ($d > 13\text{ mm}$): Chiều dày bão hòa của thép C45 đối với nguồn bức xạ $^{60}\text{Co}$ được xác định chính xác tại khoảng 12 - 13 mm. Với các khuyết tật nằm sâu hơn 13 mm (bề dày còn lại $> 13\text{ mm}$), số đếm đỉnh đạt trạng thái tiệm cận bão hòa, hiện tượng chồng chập phổ xảy ra khiến việc phân biệt kích thước và độ sâu khuyết tật không còn đạt độ tin cậy cao.
Đổi mới và đóng góp
- Làm rõ bản chất phổ tán xạ nhiều lần: Phân lập thành công thành phần tán xạ đơn (Single Scattering) và tán xạ đôi (Double Scattering). Mô hình 2-Gauss chứng minh rằng tán xạ đôi chiếm tỉ trọng lớn nhất trong thành phần tán xạ nhiều lần và việc tách đỉnh giúp làm mịn sai số thống kê ở mức $\delta S \le 0.238%$.
- Xác lập giới hạn vật lý định lượng: Khẳng định bằng thực nghiệm số rằng bề dày bão hòa của vật liệu thép C45 với chùm gamma $^{60}\text{Co}$ góc $120^\circ$ là $12.5 \pm 0.5\text{ mm}$.
- Hiệu quả tiết kiệm và an toàn bức xạ: Giảm thiểu 100% rủi ro phơi nhiễm bức xạ trong giai đoạn thiết kế ban đầu, đồng thời cắt giảm hơn $90%$ chi phí gia công mẫu cơ khí thử nghiệm nhờ tối ưu hóa kích thước Collimator và góc thu nhận trước khi lắp đặt thực tế.
Ứng dụng thực tế và triển khai
Tình huống ứng dụng công nghiệp
- Kiểm tra hiện tượng ăn mòn cục bộ đáy bồn bể chứa hóa chất/dầu khí: Cho phép rà quét từ mặt ngoài bồn mà không cần rút cạn dung dịch hay làm sạch bên trong.
- Đo độ mỏng thành ống dẫn hơi áp lực cao trong nhà máy nhiệt điện: Phát hiện hiện tượng xói mòn/ăn mòn do dòng chảy (Flow-Accelerated Corrosion) ở các vị trí cút cong.
- Giám sát chất lượng mối hàn kết cấu thép chịu lực: Phát hiện các lỗ rỗ khí ngầm có đường kính từ 2.2 cm trở lên với chiều sâu vỏ kim loại bảo vệ dưới 9 mm.
+-----------------------------------------------------------------------------------+
| LỘ TRÌNH TRIỂN KHAI HỆ THỐNG TỪNG GIAI ĐOẠN (ROADMAP) |
+-----------------------------------------------------------------------------------+
| Phase 1: Benchmark mô phỏng & Tối ưu Collimator (Hoàn thành) |
| - Hoàn thiện mô hình MCNP5, sai số lý thuyết < 1.9%. |
+-----------------------------------------------------------------------------------+
| Phase 2: Chế tạo Module đầu dò NaI(Tl) + Collimator cơ khí (3 - 6 tháng) |
| - Gia công cụm chuẩn trực chì góc phát 90 độ, góc thu 120 độ. |
| - Tích hợp khối xử lý tín hiệu số MCA nhúng. |
+-----------------------------------------------------------------------------------+
| Phase 3: Thử nghiệm Pilot tại nhà máy cơ khí / lọc dầu (6 - 12 tháng) |
| - Kiểm chuẩn trên mẫu ống thép C45 chuẩn hóa. |
| - Tích hợp thuật toán tự động khớp phổ cảnh báo khuyết tật thời gian thực.|
+-----------------------------------------------------------------------------------+
Hạn chế và hướng phát triển
- Hạn chế kỹ thuật:
- Giới hạn chiều sâu phát hiện khuyết tật: Chỉ đạt độ nhạy tối ưu khi bề dày kim loại còn lại nhỏ hơn 9 mm; mất khả năng phân biệt khi vượt ngưỡng bão hòa 13 mm.
- Đầu dò nhấp nháy $\text{NaI(Tl)}$ có độ phân giải năng lượng ở mức trung bình ($\approx 7 - 8%$ tại đỉnh $662\text{ keV}$), dễ bị chồng lấn giữa đỉnh tán xạ đơn và phổ liên tục khi góc mở collimator rộng.
- Hướng phát triển:
- Ứng dụng đầu dò bán dẫn siêu tinh khiết HPGe (High-Purity Germanium) hoặc CZT (Cadmium Zinc Telluride) nhằm tăng độ phân giải phổ năng lượng mà không cần làm lạnh bằng Nitơ lỏng phức tạp.
- Sử dụng mạng nơ-ron tích chập (CNN / Deep Learning) để tự động hóa quá trình "Unfolding Spectrum" (giải chập quang phổ), nâng cao độ nhạy phát hiện khuyết tật ở vùng sau bão hòa ($> 13\text{ mm}$).
- Khảo sát mở rộng trên các nguồn bức xạ khác như $^{137}\text{Cs}$ ($662\text{ keV}$) để tăng diện tích tán xạ vi phân Compton ở vùng năng lượng trung bình.
Đối tượng hưởng lợi
- Sinh viên & Học viên ngành Kỹ thuật Hạt nhân / Vật lý tính toán: Nguồn tài liệu tham khảo chuẩn mực về phương pháp xây dựng Input File cho chương trình MCNP5, cách khai báo thẻ F8 kết hợp xử lý phổ bằng COLEGRAM.
- Kỹ sư NDT & Đảm bảo chất lượng (QA/QC): Nắm bắt các thông số vật lý then chốt (bề dày bão hòa, hệ số khuyết tật $C$) để thiết kế cấu hình đo tán xạ một phía tại hiện trường.
- Doanh nghiệp vận hành nhà máy công nghiệp nặng: Giải pháp công nghệ kiểm tra không phá hủy online giúp giảm thiểu thời gian dừng máy (downtime), tiết kiệm hàng trăm triệu đồng chi phí bảo trì định kỳ.
Câu hỏi thường gặp
1. Yêu cầu cấu hình phần cứng để chạy mô phỏng MCNP5 với $3 \times 10^9$ hạt?
Để mô phỏng $3 \times 10^9$ photon với Tally 8 đạt sai số thống kê dưới $0.2%$, hệ thống khuyến nghị sử dụng CPU đa nhân (từ 8 - 16 Cores, xung nhịp $\ge 3.0\text{ GHz}$), RAM tối thiểu 16GB. Mô phỏng hỗ trợ chạy song song qua giao thức OpenMP hoặc MPI trên nền tảng Linux/Windows.
2. Tại sao năng lượng đỉnh tán xạ thực tế (263 keV) lại thấp hơn tính toán lý thuyết (268 keV)?
Sự chênh lệch 5 keV ($1.85%$) xuất phát từ hiệu ứng mở rộng Doppler (Doppler Broadening) do electron trong nguyên tử sắt của thép C45 có xung lượng nội tại ban đầu (không phải electron tự do đứng yên tuyệt đối) và góc tán xạ thực tế thu vào đầu dò là một nón góc khối tích phân quanh góc $120^\circ$, không phải một tia hình học đơn lẻ.
3. Phương pháp này có thể áp dụng cho các vật liệu khác ngoài thép C45 không?
Hoàn toàn có thể. Tuy nhiên, chiều dày bão hòa $d_{\text{sat}}$ sẽ thay đổi tỉ lệ nghịch với mật độ electron của vật liệu. Với nhôm ($Z=13, \rho = 2.7\text{ g/cm}^3$), bề dày bão hòa sẽ lớn hơn ($\approx 35 - 40\text{ mm}$); với chì ($Z=82, \rho = 11.34\text{ g/cm}^3$), bề dày bão hòa chỉ vào khoảng vài milimet.
4. Tại sao lại chọn góc tán xạ 120 độ thay vì 180 độ (tán xạ ngược hoàn toàn)?
Tại góc $120^\circ$, tiết diện tán xạ vi phân Compton Klein-Nishina lớn hơn so với tại góc $180^\circ$, giúp tăng hiệu suất ghi nhận xung đếm của đầu dò. Đồng thời, cấu hình $120^\circ$ cho phép bố trí buồng chì chuẩn trực của nguồn và vỏ bọc đầu dò cạnh nhau mà không bị cản trở không gian vật lý.
5. So sánh sai số giữa phương pháp làm khớp 1 đỉnh Gauss và 2 đỉnh Gauss?
Cả hai phương pháp đều cho kết quả hệ số khuyết tật $C (%)$ gần như tương đương nhau (chênh lệch $< 0.3%$). Tuy nhiên, phương pháp 2-Gauss phản ánh chính xác bản chất vật lý của sự đóng góp từ tán xạ đôi, giúp ổn định đường nền khi xử lý các phổ có độ sâu khuyết tật lớn.
Kết luận
Đề tài đã ứng dụng thành công phương pháp Monte Carlo thông qua mã MCNP5 để mô phỏng và khảo sát toàn diện hệ đo gamma tán xạ ngược dùng nguồn $^{60}\text{Co}$ và đầu dò $\text{NaI(Tl)}$ nhằm xác định vị trí khuyết tật trên thép C45. Kết quả chứng minh tính khả thi vượt trội của phương pháp trong việc phát hiện khuyết tật ngầm có bề dày kim loại còn lại dưới 9 mm với độ tin cậy quang phổ đạt trên $98.1%$. Đây là tiền đề khoa học và kỹ thuật vững chắc để tiến tới chế tạo các thiết bị NDT gamma tán xạ ngược "Made in Vietnam", phục vụ đắc lực cho công tác kiểm định an toàn công nghiệp và dầu khí.
