Nghiên Cứu Phương Pháp Monte Carlo Trong Luận Văn Tốt Nghiệp

Khóa luận tốt nghiệp nghiên cứu tốt nghiệp vật lý đánh giá kỹ thuật xử lý phổ cải tiến áp dụng cho việc xác định bề dày nhiều loại, vận dụng lý thuyết vào thực tế, đề xuất giải

• Tác giả: Nguyen Thi My Le
• Người hướng dẫn: TS. Hoang Duc Tam
• Trường học: Trường Đại Học Sư Phạm Thành Phố Hồ Chí Minh
• Chuyên ngành: Sư Phạm Vật Lý
• Đề tài: Nghiên Cứu Phương Pháp Monte Carlo Trong Vật Lý Hạt Nhân
• Loại tài liệu: Luận Văn Tốt Nghiệp
• Năm xuất bản: 2016
• Địa điểm: Thành Phố Hồ Chí Minh

Mục lục chi tiết

• LỜI CÁM ƠN
• MỞ ĐẦU
• 1. CHƯƠNG 1: TỔNG QUAN
• 1.1. Sơ lược về tình hình nghiên cứu phương pháp gamma tán xạ ngược
• 1.2. Tình hình nghiên cứu trên thế giới
• 1.3. Tình hình nghiên cứu trong nước
• 1.4. Giới thiệu phương pháp Monte Carlo và chương trình MCNP
• 1.4.1. Phương pháp Monte Carlo
• 1.4.2. Chương trình MCNP
• 1.4.3. Mô hình tương tác của photon với vật chất trong MCNP
• 1.4.3.1. Tán xạ Compton
• 1.4.3.2. Hiệu ứng quang điện
• 2. CHƯƠNG 2: XÂY DỰNG MÔ HÌNH MÔ PHỎNG HỆ ĐO GAMMA TÁN XẠ NGƯỢC
• 2.1. Đầu dò và ống chuẩn
• 2.2. Đường chuẩn năng lượng của hệ đo
• 2.3. Tổng kết chương 2
• 3. CHƯƠNG 3: KẾT QUẢ VÀ THẢO LUẬN
• 3.1. Khảo sát các đặc trưng của phô gamma tán xạ ngược
• 3.2. Dạng phô gamma tán xạ ngược
• 3.3. Khảo sát theo bậc số nguyên tử Z
• 3.4. Khảo sát theo bề dày bia
• 3.5. Khảo sát theo năng lượng
• 3.6. Tính toán bề dày bia
• 3.7. Kết luận chương 3
• DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT
• DANH MỤC CÁC BẢNG BIỂU
• DANH MỤC CÁC HÌNH VẼ VÀ ĐỒ THỊ
• DANH MỤC CÁC CÔNG TRÌNH CỦA TÁC GIẢ
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng Quan Về Nghiên Cứu Phương Pháp Monte Carlo Trong Vật Lý Hạt Nhân

Phương pháp Monte Carlo đã trở thành một công cụ quan trọng trong nghiên cứu vật lý hạt nhân. Phương pháp này cho phép mô phỏng các quá trình phức tạp mà không thể giải quyết bằng các phương pháp phân tích truyền thống. Đặc biệt, trong lĩnh vực vật lý hạt nhân, Monte Carlo giúp nghiên cứu các tương tác giữa photon và vật chất, từ đó cung cấp những hiểu biết sâu sắc về cấu trúc và tính chất của vật liệu. Việc áp dụng phương pháp này đã mở ra nhiều hướng nghiên cứu mới, giúp cải thiện độ chính xác trong các phép đo và mô phỏng.

1.1. Tình Hình Nghiên Cứu Phương Pháp Monte Carlo Trên Thế Giới

Trên thế giới, nhiều nghiên cứu đã chỉ ra rằng phương pháp Monte Carlo có thể áp dụng hiệu quả trong việc mô phỏng các quá trình hạt nhân. Các nghiên cứu như của Paramesh và cộng sự (1983) đã xác định sự phụ thuộc của bề dày bão hòa vào bậc số nguyên tử Z, cho thấy độ chính xác cao của phương pháp này.

1.2. Tình Hình Nghiên Cứu Phương Pháp Monte Carlo Tại Việt Nam

Tại Việt Nam, nghiên cứu về phương pháp Monte Carlo cũng đã được thực hiện với nhiều kết quả đáng chú ý. Hoàng Đức Tâm và các cộng sự đã sử dụng chương trình MCNP để mô phỏng bề dày bão hòa của vật liệu thép, cho thấy sự phù hợp giữa kết quả mô phỏng và thực nghiệm.

II. Vấn Đề Và Thách Thức Trong Nghiên Cứu Phương Pháp Monte Carlo

Mặc dù phương pháp Monte Carlo mang lại nhiều lợi ích, nhưng cũng tồn tại một số thách thức trong quá trình nghiên cứu. Một trong những vấn đề chính là độ chính xác của các số liệu đầu vào. Nếu số liệu không chính xác, kết quả mô phỏng sẽ bị sai lệch. Hơn nữa, việc tính toán với số lượng lớn sự kiện cũng đòi hỏi tài nguyên tính toán cao, điều này có thể gây khó khăn cho các nghiên cứu quy mô lớn.

2.1. Độ Chính Xác Của Dữ Liệu Đầu Vào

Độ chính xác của dữ liệu đầu vào là yếu tố quyết định đến kết quả mô phỏng. Việc sử dụng các số liệu hạt nhân không chính xác có thể dẫn đến sai số lớn trong các phép đo và phân tích.

2.2. Tài Nguyên Tính Toán Cần Thiết

Phương pháp Monte Carlo yêu cầu một lượng lớn tài nguyên tính toán để thực hiện các mô phỏng. Điều này có thể gây khó khăn cho các nghiên cứu không có đủ nguồn lực.

III. Phương Pháp Mô Phỏng Monte Carlo Trong Nghiên Cứu Vật Lý Hạt Nhân

Phương pháp mô phỏng Monte Carlo sử dụng các số ngẫu nhiên để mô phỏng các quá trình vật lý phức tạp. Trong nghiên cứu vật lý hạt nhân, phương pháp này cho phép mô phỏng các tương tác giữa photon và vật chất, từ đó giúp xác định các đặc trưng của vật liệu. Chương trình MCNP là một trong những công cụ phổ biến nhất được sử dụng để thực hiện các mô phỏng này.

3.1. Nguyên Tắc Hoạt Động Của Phương Pháp Monte Carlo

Phương pháp Monte Carlo hoạt động dựa trên việc tạo ra các số ngẫu nhiên và sử dụng chúng để mô phỏng các sự kiện vật lý. Điều này cho phép nghiên cứu các quá trình phức tạp mà không cần phải giải quyết các phương trình vật lý phức tạp.

3.2. Chương Trình MCNP Và Ứng Dụng Của Nó

Chương trình MCNP (Monte Carlo N-Particles) được phát triển để mô phỏng quá trình vận chuyển của neutron, photon và electron. Chương trình này cho phép nghiên cứu các tương tác giữa các hạt và vật chất, từ đó cung cấp những thông tin quý giá cho các nghiên cứu vật lý hạt nhân.

IV. Ứng Dụng Thực Tiễn Của Phương Pháp Monte Carlo Trong Vật Lý Hạt Nhân

Phương pháp Monte Carlo đã được áp dụng rộng rãi trong nhiều lĩnh vực của vật lý hạt nhân, từ nghiên cứu cấu trúc vật liệu đến phân tích các quá trình tương tác hạt. Các ứng dụng này không chỉ giúp cải thiện độ chính xác của các phép đo mà còn mở ra nhiều hướng nghiên cứu mới. Việc sử dụng phương pháp này trong các thí nghiệm thực tế đã cho thấy những kết quả khả quan.

4.1. Nghiên Cứu Cấu Trúc Vật Liệu

Phương pháp Monte Carlo cho phép nghiên cứu cấu trúc vật liệu một cách chi tiết, từ đó giúp xác định các đặc trưng vật lý của chúng. Các nghiên cứu này có thể ứng dụng trong việc phát triển vật liệu mới.

4.2. Phân Tích Các Quá Trình Tương Tác Hạt

Phương pháp này cũng được sử dụng để phân tích các quá trình tương tác giữa các hạt, giúp hiểu rõ hơn về các hiện tượng vật lý phức tạp trong vật lý hạt nhân.

V. Kết Luận Về Nghiên Cứu Phương Pháp Monte Carlo Trong Vật Lý Hạt Nhân

Nghiên cứu phương pháp Monte Carlo trong vật lý hạt nhân đã chứng minh được giá trị của nó trong việc mô phỏng và phân tích các quá trình phức tạp. Mặc dù còn tồn tại một số thách thức, nhưng những lợi ích mà phương pháp này mang lại là không thể phủ nhận. Tương lai của nghiên cứu này hứa hẹn sẽ tiếp tục phát triển với sự hỗ trợ của công nghệ tính toán ngày càng tiên tiến.

5.1. Triển Vọng Phát Triển Của Phương Pháp Monte Carlo

Với sự phát triển của công nghệ tính toán, phương pháp Monte Carlo sẽ ngày càng trở nên mạnh mẽ và chính xác hơn. Điều này mở ra nhiều cơ hội cho các nghiên cứu mới trong lĩnh vực vật lý hạt nhân.

5.2. Tương Lai Của Nghiên Cứu Vật Lý Hạt Nhân

Nghiên cứu vật lý hạt nhân sẽ tiếp tục đóng vai trò quan trọng trong việc phát triển các công nghệ mới, từ năng lượng hạt nhân đến y học hạt nhân. Phương pháp Monte Carlo sẽ là một công cụ không thể thiếu trong các nghiên cứu này.

Nội dung chính

Tổng quan nghiên cứu

Trong kỹ thuật kiểm tra không phá hủy hiện đại (Non-Destructive Testing - NDT), việc xác định độ dày và phát hiện ăn mòn thành bình, đường ống kim loại đóng vai trò sống còn đối với an toàn công nghiệp dầu khí, hóa chất và năng lượng hạt nhân. Theo ước tính từ các báo cáo kỹ thuật ngành, hơn 60% sự cố rò rỉ cơ học bắt nguồn từ hiện tượng ăn mòn cục bộ mà các phương pháp đo truyền qua thông thường khó tiếp cận được do hạn chế về không gian hình học. Phương pháp gamma tán xạ ngược nổi lên như một giải pháp vượt trội với ưu thế không cần tiếp xúc trực tiếp và chỉ yêu cầu tiếp cận đối tượng từ một phía duy nhất, cho phép thực hiện đo đạc ngay cả khi hệ thống đang vận hành ở nhiệt độ và áp suất cao.

Nghiên cứu tập trung khảo sát sự biến thiên đặc trưng phổ gamma tán xạ ngược theo bậc số nguyên tử $Z$ bằng phương pháp mô phỏng Monte Carlo trên nền tảng chương trình MCNP5. Mục tiêu trọng tâm là thiết lập mô hình hệ đo tán xạ tại góc $120^\circ$, đánh giá kỹ thuật xử lý phổ cải tiến và tính toán chính xác bề dày cũng như bề dày bão hòa trên 9 loại vật liệu tấm phẳng ($50\text{ cm} \times 50\text{ cm}$) gồm Cacbon ($Z=6$), Nhôm ($Z=13$), Sắt ($Z=26$), Đồng ($Z=29$), Kẽm ($Z=30$), Thiếc ($Z=50$), Bạc ($Z=47$), Vàng ($Z=79$) và Chì ($Z=82$). Nghiên cứu được thực hiện tại phòng thí nghiệm Bộ môn Vật lý Hạt nhân, Trường Đại học Sư phạm Thành phố Hồ Chí Minh năm 2016. Kết quả nghiên cứu kiểm soát sai số tính toán bề dày dưới 7,6%, cung cấp bộ thông số chuẩn xác cho việc chế tạo và tối ưu hóa hệ đo NDT thực nghiệm ngoài hiện trường.

---

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

Khung lý thuyết áp dụng

Quá trình tương tác của bức xạ gamma với vật chất trong nghiên cứu được xây dựng trên nền tảng 3 hiệu ứng vật lý lượng tử cốt lõi: tán xạ Compton, hiệu ứng quang điện và hiệu ứng tạo cặp:

1. Tán xạ Compton: Đóng vai trò chủ đạo ở vùng năng lượng trung bình. Khi photon tới có năng lượng $E_0 = 661,657\text{ keV}$ va chạm với electron liên kết yếu, photon tán xạ với góc $\theta = 120^\circ$ mang năng lượng $E$ tuân theo hệ thức bảo toàn năng lượng và xung lượng:

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

Tiết diện tán xạ vi phân được mô tả chính xác qua công thức Klein-Nishina kết hợp hàm tán xạ không kết hợp $S(q, Z)$.

2. Hiệu ứng quang điện: Chiếm ưu thế tuyệt đối ở vùng năng lượng thấp và vật liệu có bậc số nguyên tử $Z$ cao, với tiết diện tỉ lệ thuận xấp xỉ $Z^4 / (h\nu)^3$. Hiệu ứng này tạo ra các lỗ trống điện tử, dẫn đến bức xạ tia X đặc trưng (điển hình là đỉnh tia X $75\text{ keV}$ của Chì) và phát xạ electron Auger.

3. Hiệu ứng tạo cặp: Xảy ra khi năng lượng photon vượt ngưỡng $1,022\text{ MeV}$, tạo cặp electron - positron trong điện trường hạt nhân.

Sự suy giảm cường độ chùm photon hẹp truyền qua vật chất tuân theo định luật Beer-Lambert:

$$I = I_0 \exp(-\mu x)$$

Trong đó $\mu$ là hệ số suy giảm tuyến tính. Khi xét photon tán xạ một lần từ lớp vật liệu có bề dày $T$, cường độ chùm photon $I(T)$ và bề dày bão hòa $T_0$ được mô hình hóa qua hàm mũ:

$$I(T) = I_0 [1 - \exp(-\mu_t T)]$$

Với $\mu_t$ là tổng hệ số suy giảm toàn phần hiệu dụng.

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

Nghiên cứu ứng dụng chương trình mô phỏng hạt nhân MCNP5 (Monte Carlo N-Particles) để giải bài toán vận chuyển photon trong môi trường không gian 3 chiều:

• Nguồn dữ liệu và tiết diện: Sử dụng thư viện dữ liệu hạt nhân chuẩn ENDF và bảng tiết diện tương tác photon từ Viện Tiêu chuẩn và Công nghệ Quốc gia Hoa Kỳ (NIST XCOM) cho toàn bộ các nguyên tố từ $Z=1$ đến $Z=94$.
• Cỡ mẫu và phương pháp chọn mẫu: Phương pháp Monte Carlo thực hiện lấy mẫu thống kê ngẫu nhiên từ phân bố đều trên đoạn $[0, 1]$. Để đảm bảo phổ phân bố độ cao xung (PHS) có sai số tương đối $R < 0,05$ (đạt độ tin cậy 95% theo định lý giới hạn trung tâm), cỡ mẫu mô phỏng được thiết lập với $6 \times 10^7$ (60 triệu) lịch sử hạt nguồn.
• Cấu hình hình học hệ đo: Nguồn phóng xạ đồng vị Cs-137 ($661,657\text{ keV}$, chu kỳ bán rã 30,05 năm) đặt trong hộp chì có ống chuẩn trực dài $20\text{ cm}$, đường kính trong $1\text{ cm}$. Đầu dò nhấp nháy NaI(Tl) kích thước $76\text{ mm} \times 76\text{ mm}$ gắn ống nhân quang điện nhôm ($\varnothing 83,2\text{ mm} \times 30\text{ mm}$). Khoảng cách nguồn - bia là $34\text{ cm}$, bia - đầu dò là $31\text{ cm}$, góc tán xạ thiết lập ở $120^\circ$.
• Xử lý phổ và phân tích dữ liệu: Hệ phổ kế phân tích đa kênh (MCA) chia thành 8192 kênh với phương trình chuẩn năng lượng tuyến tính thực nghiệm: $E\text{ (keV)} = 0,20407 \times \text{Ch} - 13,5363$ ($R^2 = 1,000$). Chương trình Colegram được sử dụng để tách đỉnh tán xạ đơn khỏi nền tán xạ nhiều lần trong dải kênh 604 đến 1772 bằng cách khớp đồng thời hàm phân bố Gauss (thẻ GEB) và hàm đa thức bậc 4:

$$G(x) = \frac{A}{\sigma \sqrt{2\pi}} \exp\left[-\frac{(x - x_0)^2}{2\sigma^2}\right]$$

$$\text{poly}(x) = a_0 + a_1(x - x_0) + a_2(x - x_0)^2 + a_3(x - x_0)^3 + a_4(x - x_0)^4$$

Phương án 1 (cho $Z$ nhỏ và trung bình) sử dụng 2 hàm Gauss kết hợp đa thức bậc 4. Phương án 2 (cho $Z$ lớn) sử dụng 1 hàm Gauss kết hợp đa thức bậc 4. Dữ liệu phi tuyến sau đó được làm khớp trên phần mềm Origin 8.

---

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

1. Độ chuẩn xác về vị trí năng lượng đỉnh tán xạ đơn: Trên phổ mô phỏng của bia Sắt (Fe) dày $2,334\text{ cm}$, đỉnh tán xạ một lần xuất hiện tại kênh 1165, tương ứng năng lượng $224,205\text{ keV}$. So với giá trị lý thuyết $224,882\text{ keV}$, độ sai lệch chỉ đạt $0,33%$. Đối với 4 nguồn phóng xạ khác nhau gồm Cs-137 ($662\text{ keV}$), Mn-54 ($835\text{ keV}$), Zn-65 ($1116\text{ keV}$) và Co-60 ($1173/1332\text{ keV}$), độ lệch năng lượng tán xạ giữa mô phỏng và lý thuyết đều được duy trì dưới mức $1,0%$.

2. Quy luật biến thiên của bề dày bão hòa theo số hiệu nguyên tử $Z$: Kết quả mô phỏng làm khớp hàm mũ $I(T)$ cho thấy bề dày bão hòa $T_0$ giảm mạnh khi $Z$ tăng: Cacbon ($Z=6$) đạt $8,00\text{ cm}$, Nhôm ($Z=13$) đạt $7,06\text{ cm}$, Sắt ($Z=26$) đạt $2,03\text{ cm}$, Đồng ($Z=29$) đạt $2,33\text{ cm}$, Kẽm ($Z=30$) đạt $1,83\text{ cm}$, Vàng ($Z=79$) đạt $0,34\text{ cm}$, và Chì ($Z=82$) đạt $0,40\text{ cm}$.

Vật liệu	Số hiệu nguyên tử ($Z$)	Bề dày bão hòa mô phỏng (cm)	Hệ số suy giảm $\mu_t$ lý thuyết ($\text{cm}^{-1}$)	Hệ số suy giảm $\mu_t$ mô phỏng ($\text{cm}^{-1}$)	Độ sai biệt $\mu_t$ (%)
Cacbon (C)	6	8,00	0,382	0,387	1,3%
Nhôm (Al)	13	7,06	0,587	0,600	2,2%
Sắt (Fe)	26	2,03	1,793	1,909	6,1%
Đồng (Cu)	29	2,33	2,081	2,196	5,2%
Kẽm (Zn)	30	1,83	1,692	1,783	5,1%
Vàng (Au)	79	0,34	16,798	17,983	6,7%
Chì (Pb)	82	0,40	10,635	11,359	6,4%

3. Độ tương thích cao của hệ số suy giảm tuyến tính: Hệ số suy giảm tuyến tính $\mu_t$ trích xuất từ đường cong làm khớp mô phỏng so với dữ liệu chuẩn NIST XCOM đạt độ sai biệt dao động từ $1,3%$ (đối với C) đến $6,7%$ (đối với Au). Toàn bộ các phép làm khớp phi tuyến đều đạt hệ số tương quan xuất sắc $R^2 \ge 0,999$.

4. Độ chính xác trong xác định độ dày vật liệu: Sai số tương đối trong tính toán bề dày thực tế từ diện tích đỉnh tán xạ đơn cho tất cả các nhóm vật liệu $Z$ thấp, $Z$ trung bình và $Z$ cao đều không vượt quá $7,6%$.

Thảo luận kết quả

Cơ chế vật lý chi phối hình dạng phổ gamma tán xạ ngược là sự cạnh tranh gay gắt giữa tán xạ Compton và hấp thụ quang điện:

• Đối với các vật liệu có $Z$ trung bình ($Z = 20 - 30$ như Fe, Cu, Zn), tiết diện tán xạ Compton chiếm ưu thế áp đảo, tạo ra diện tích đỉnh tán xạ một lần có cường độ cao nhất trên phổ PHS. Điều này mang lại độ tin cậy thống kê tối ưu khi sử dụng nguồn Cs-137 để đo kiểm độ dày nhóm kim loại công nghiệp này.
• Ngược lại, ở vật liệu có $Z$ thấp ($Z \le 13$ như C, Al), tán xạ nhiều lần (chủ yếu là tán xạ đôi) đóng góp mạnh làm dâng cao chân trái của đỉnh phổ, đòi hỏi bắt buộc phải tách phổ theo Phương án 1 (khớp 2 hàm Gauss). Ở vật liệu có $Z$ cao ($Z \ge 79$ như Au, Pb), hiệu ứng quang điện triệt tiêu phần lớn photon tán xạ, làm bề dày bão hòa thu hẹp chỉ còn dưới $0,4\text{ cm}$.

Đồ thị đường cong bão hòa biểu diễn quan hệ $I(T)$ theo bề dày bia minh họa rõ nét hiện tượng tiệm cận: khi bề dày $T$ vượt qua ngưỡng bão hòa $T_0$, xác suất photon bị hấp thụ trên đường thoát ra ngoài cân bằng hoàn toàn với xác suất tán xạ mới tạo thành, khiến số đếm ghi nhận tại đầu dò không còn thay đổi. So sánh với các nghiên cứu thực nghiệm trước đây, giá trị bề dày bão hòa của Sắt ($2,03\text{ cm}$) hoàn toàn tương thích với các công bố thực nghiệm của TS. Hoàng Đức Tâm ($1,7 - 2,0\text{ cm}$), và độ chênh lệch bề dày bão hòa với nghiên cứu của Paramesh (1983) đối với Nhôm là $1,24\text{ cm}$ (tương đương $14,9%$). Sự chênh lệch độ rộng FWHM mô phỏng ($22,6\text{ keV}$ ở Fe dày $0,34\text{ cm}$) so với hàm làm khớp Gauss ($20,7\text{ keV}$) chứng minh rằng độ rộng đỉnh không chỉ phụ thuộc vào hiệu ứng mở rộng năng lượng của tinh thể NaI(Tl) mà còn chịu tác động đáng kể từ tán xạ bội trong thể tích hình học của bia.

---

Đề xuất và khuyến nghị

1. Tối ưu hóa thiết kế ống chuẩn trực và che chắn đầu dò: Nhóm kỹ sư phần cứng hạt nhân cần gia công nâng cấp ống chuẩn trực chì với góc mở hẹp $15^\circ$ và bổ sung lớp đệm đồng - nhôm lót trong để triệt tiêu đỉnh huỳnh quang tia X $75\text{ keV}$ của chì, mục tiêu giảm trên $20%$ nhiễu nền tán xạ trong vòng 3 tháng đầu triển khai.
2. Tự động hóa quy trình phân tích phổ thời gian thực: Nhóm phát triển phần mềm NDT cần tích hợp thuật toán làm khớp 2 Gauss kết hợp đa thức bậc 4 vào trực tiếp bộ phân tích đa kênh MCA, nhằm rút ngắn thời gian xử lý diện tích đỉnh tán xạ từ 15 phút xuống dưới 30 giây cho mỗi chu kỳ đo trong lộ trình 6 tháng.
3. Mở rộng dải nguồn bức xạ đa năng lượng: Các viện nghiên cứu hạt nhân chuyên ngành nên tích hợp cấu hình nguồn kép Co-60 ($1,17/1,33\text{ MeV}$) và Cs-137 ($0,662\text{ MeV}$) để nâng giới hạn đo bề dày bão hòa cho kết cấu thép dày trên $50\text{ mm}$, khống chế sai số đo dưới $5,0%$ trong vòng 12 tháng.
4. Ban hành quy chuẩn đo kiểm NDT một phía cho đường ống công nghiệp: Cơ quan kiểm định kỹ thuật an toàn công nghiệp cần áp dụng mô hình toán học giải tích từ nghiên cứu này để xây dựng bộ tiêu chuẩn đánh giá độ mỏng thành ống dẫn khí và bình chịu áp lực, đảm bảo độ lặp lại phép đo đạt $98%$ trong thời hạn 9 tháng.

---

Đối tượng nên tham khảo luận văn

• Kỹ sư kiểm tra không phá hủy (NDT): Ứng dụng mô hình hình học tán xạ $120^\circ$ để chế tạo thiết bị đo bề dày thành ống dẫn dầu, bồn chứa hóa chất bị ăn mòn mà không cần dừng dây chuyền sản xuất.
• Nghiên cứu sinh và Giảng viên Vật lý Hạt nhân: Khai thác phương pháp viết code thẻ MCNP5 (khai báo cell card, surface card, bộ đánh giá xung F8, card mở rộng GEB) và quy trình chuẩn năng lượng 8192 kênh làm học liệu giảng dạy chuyên sâu.
• Chuyên gia an toàn bức xạ và thiết kế che chắn: Sử dụng hệ số suy giảm tuyến tính $\mu_t$ và bề dày bão hòa của 9 loại vật liệu từ Cacbon đến Chì để tính toán chiều dày che chắn bức xạ gamma tối ưu cho các cơ sở hạt nhân.
• Nhà phát triển thiết bị đo lường công nghiệp: Sử dụng thuật toán tách phổ Colegram cải tiến để lập trình vi điều khiển cho các máy đo mật độ và đo độ dày trực tuyến di động.

---

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

1. Tại sao phương pháp gamma tán xạ ngược chỉ đo được độ dày nhỏ hơn bề dày bão hòa?

Khi độ dày vật liệu tăng đến giá trị bão hòa ($T_0$), số lượng photon tán xạ ngược tăng lên nhưng đồng thời quãng đường photon di chuyển trong vật chất cũng dài hơn, khiến xác suất hấp thụ quang điện tăng theo. Hai quá trình này triệt tiêu lẫn nhau làm cường độ chùm tia đến đầu dò không đổi, khiến mọi biến thiên bề dày lớn hơn $T_0$ không thể ghi nhận được.

2. Vì sao cần áp dụng hai phương án xử lý phổ khác nhau cho các nhóm vật liệu $Z$ khác nhau?

Với vật liệu $Z$ thấp và trung bình ($Z \le 30$), xác suất tán xạ bậc hai (tán xạ đôi) khá lớn, làm méo mó chân đỉnh tán xạ đơn, do đó cần 2 hàm Gauss để tách riêng. Với vật liệu $Z$ cao ($Z \ge 79$), hiệu ứng quang điện chiếm ưu thế làm triệt tiêu tán xạ đôi, do đó chỉ cần 1 hàm Gauss kết hợp đa thức bậc 4 để khử nền.

3. Ý nghĩa của việc chọn $6 \times 10^7$ lịch sử hạt trong mô phỏng MCNP5 là gì?

Trong mô phỏng Monte Carlo, sai số tương đối $R$ tỉ lệ nghịch với căn bậc hai của số lịch sử hạt ($\sqrt{N}$). Việc sử dụng 60 triệu hạt bảo đảm sai số thống kê phổ xung F8 đạt mức dưới $1,0%$, đưa giá trị $R < 0,05$, đáp ứng tiêu chuẩn nghiêm ngặt về độ tin cậy của đầu dò điểm trong mô phỏng hạt nhân.

4. Nguồn phóng xạ Cs-137 có ưu thế gì so với Co-60 hay Mn-54 trong bài toán này?

Cs-137 phát photon đơn năng $661,657\text{ keV}$ rõ ràng và có chu kỳ bán rã dài 30,05 năm, giúp chùm tán xạ một lần không bị chồng chập bởi nhiều đỉnh năng lượng phức tạp như Co-60 ($1173\text{ keV}$ và $1332\text{ keV}$), đặc biệt tối ưu cho dải kim loại $Z$ trung bình từ 20 đến 30.

5. Nguyên nhân nào làm tăng độ rộng FWHM của đỉnh tán xạ trong mô phỏng so với tính toán lý thuyết?

Ngoài sự giãn nở năng lượng dạng Gauss do đặc tính tinh thể NaI(Tl) và hệ điện tử (khai báo qua card GEB), độ rộng FWHM thực tế trong mô phỏng còn bị mở rộng bởi yếu tố hình học của bia tán xạ dày $50\text{ cm} \times 50\text{ cm}$ và sự đóng góp của các photon tán xạ nhiều lần có góc tán xạ lệch nhẹ quanh $120^\circ$.

---

Kết luận

• Xây dựng thành công mô hình mô phỏng MCNP5 hệ đo gamma tán xạ ngược góc $120^\circ$ sử dụng nguồn Cs-137 ($661,657\text{ keV}$) và đầu dò NaI(Tl) $76\text{ mm} \times 76\text{ mm}$.
• Xác lập chính xác quy luật suy giảm bề dày bão hòa từ Cacbon ($8,00\text{ cm}$) đến Chì ($0,40\text{ cm}$), khẳng định độ nhạy vượt trội trên dải vật liệu $Z$ trung bình.
• Chứng minh tính hiệu quả của kỹ thuật xử lý phổ cải tiến 2 phương án (Gauss + đa thức bậc 4), tách triệt để đỉnh tán xạ đơn với sai số năng lượng dưới $0,33%$.
• Kiểm soát sai số tính toán bề dày vật liệu dưới $7,6%$ và sai số hệ số suy giảm tuyến tính $\mu_t$ so với chuẩn NIST dưới $6,7%$.
• Cung cấp cơ sở dữ liệu số học và giải tích chuẩn xác để chế tạo hệ thiết bị NDT đo độ dày một phía ứng dụng trong công nghiệp kiểm định chất lượng.

Trong giai đoạn tiếp theo (Quý 1 - Quý 2), nhóm nghiên cứu sẽ tiến hành thử nghiệm hệ đo thực nghiệm tại hiện trường và tích hợp thuật toán phân tích phổ tự động. Quý độc giả, kỹ sư kiểm định và nhà nghiên cứu quan tâm vui lòng kết nối để trao đổi học thuật và chuyển giao công nghệ mô phỏng hạt nhân ứng dụng.

Trích đoạn nội dung tài liệu

CHƯƠNG 1 TONG QUAN 1. Sơ lược về tinh hình nghiên cứu phương pháp gamma tán xa ngược 1. Tình hình nghiên cứu trên thế giới Năm 1983, Paramesh và cộng sự [6] đã xác định sự phụ thuộc của bẻ dày bão hòa vào bậc số nguyên tử Z sử dụng nguồn Os, Kết quả mà họ đã thu được cho thay rang, độ dày bão hòa của Al, Fe, Cu, Pb lần lượt là 8,3, 2,65, 2. Năm 2011, nhóm nghiên cứu Priyada và cộng sự [7] đã so sánh độ chính xác giữa phương pháp gamma tắn xạ ngược.

chụp anh gamma và chụp ảnh tia X trong việc xác định độ ăn mòn của thép mém. Kết quả mô phỏng Monte Carlo được nhóm tác giả sử dụng chương trình MCNP4C thực hiện. sau đó so sánh với các dữ liệu thực nghiệm và cho thay sự phù hợp khá tốt. Nhóm tác giả cũng chi ra rằng độ chính xác của các kỹ thuật này là tương đương nhau, đặc biệt kỹ thuật gamma tin xạ được sử dụng phô biến vì có nhiều ưu điểm như: có thể đo trong các điều kiện khắc nghiệt (nhiệt độ và áp suất cao), chỉ cần tiếp cận đối tượng từ một phía thậm chí là không cần dừng hoạt động các thiết bị.

Singh và cộng sự [8] đã tiến hành thực nghiệm các phép do gamma tán xạ ngược trên các bia vật liệu AI, Fe, Zn, Sn, Pb đề nghiên cứu sự phụ thuộc của bê day bão hòa vào bậc số nguyên tử Z ứng với năng lượng 662keV và đầu dò HPGe. Theo đó, nhóm đã đánh giá sự phụ thuộc của bé day bão hòa vào bậc số nguyên tử Z của vật liệu làm bia, cụ thé là bé day bão hòa tăng cùng với sự tăng của số hiệu nguyên tử. Các kết quá nghiên cứu được sử dụng đề làm cơ sở cho các phép đo thực nghiệm. Tình hình nghiên cứu trong nước Hoàng Đức Tâm và các cộng sự [1] đã mô phỏng Monte Carlo bằng chương trình MCNP§ dé tìm bề dày bão hòa của vật liệu thép chịu nhiệt đối với chùm tia nm gamma tới có năng lượng 662keV của nguồn '“’Cs, kết quả cho thay bề day bão hòa đối với góc tán xạ 135” là 17mm.

Hoàng Đức Tâm và các cộng sự [5] đã tiến hành thực nghiệm và mô phỏng bằng chương trình MCNPS cho các phép đo gamma tan xạ ngược trên vật liệu thép C45, sử dụng nguồn ' “Cs và đầu đò Nal(TD). Trong công trình này, các tác giả đã đưa ra biêu thức giải tích đẻ tính toán bẻ dày vật liệu và kỹ thuật xử lý pho cải tiền dé phân tích phô tán xạ từ dit liệu mô phỏng và thực nghiệm. Kết quả cho thay sự phù hợp khá tốt giữa thực nghiệm và mô phỏng, độ sai biệt giữa độ dày thực tế và độ day tính toán được la đưới 4%. Nhằm thực hiện dé tài nghiên cứu với tên dé tài: "Đánh giá kỹ thuật xử lý pho cải tiến áp dụng cho việc xác định be day nhiều loại vật liệu Z khác nhau bằng phương pháp Monte Carlo”, trong luận văn nảy tác giả đã thiết lập một mô hình mô phỏng đề xác định cường độ chùm photon tại góc tán xạ 120” trên các bia vật liệu C, AI, Fe, Cu, Zn, Sn, Ag, Au, Pb dang tắm phăng.

Bên cạnh đó, tác giả áp dụng kỹ thuật xử lý phô cải tiến trong nghiên cứu [5] dé phân tích phô tán xạ. Dựa trên kết qua mô phóng dé tìm bề day bão hòa của vật liệu, từ đó xác định bé day của các loại vật liệu khác nhau. Giới thiệu phương pháp Monte Carlo và chương trình MCNP 1. Phương pháp Monte Carlo Mô phỏng bằng phương pháp Monte Carlo là phương pháp mô phỏng trên máy tính dựa vào sự phát sinh các số ngẫu nhiên.

Do đó, phương pháp Monte Carlo cung cấp những lời giải gần đúng cho các bài toán bằng cách thực hiện các thí nghiệm lay mẫu thống kê sử dụng số ngẫu nhiên [9]. Dé giải bài toán Monte Carlo ta phải: ~ Tạo các sô ngẫu nhiên phân bố đều trên đoạn [0,1]. — Lấy mau các đại lượng ngẫu nhiên từ các phân bỗ cho trước của chúng dựa trên các sô ngẫu nhiên phân bo đều trên đoạn [0,1]. — Tính các đặc trưng trung bình được quan tâm dựa trên các giá trị của các đại lượng ngẫu nhiên được lựa chọn và xử lý thông kê.

Phương pháp Monte Carlo sử dụng các số ngẫu nhiên nên phương pháp này có thành công hay không còn tùy thuộc vào tập hợp các số ngẫu nhiên có đáng tin cậy hay không. Da phần các số ngẫu nhiên được phát sinh từ các ngôn ngữ lập trình trên máy tính có giá trị trong khoảng từ 0 đến 1, gọi là các giả ngẫu nhiên và chúng có tính tuần hoàn lặp lại sau một chu kỳ nhất định tương ứng với mỗi thuật toán. Tuy nhiên, rất khó có khả năng để những mỗi liên hệ tinh vi giữa các giá trị giả ngau nhiên có ảnh hưởng đáng kẻ lên kết quả bài toán. Cho nên trong phạm vi giới han của chu kỳ tuần hoàn, ta có thé ứng dụng các số giả ngẫu nhiên cho bài toán Monte Carlo.

Phương pháp Monte Carlo có hai đặc điểm chính: đặc điểm thứ nhất của phương pháp là thuật toán đơn giản. Khi mô phỏng, ta chỉ cần xây dựng thuật toán 31 (thuật toán tạo số ngẫu nhiên phân bố đều trên đoạn [0,1]) cho sự kiện va sau đó lặp cho tất cá các sự kiện còn lại. Vì vậy phương pháp Monte Carlo còn được gọi là phương pháp thử thông kê. Đặc điểm thứ hai của phương pháp là sai số của kết quả nhận được tỉ lệ với đại lượng ÍN gọi là độ lệch chuẩn.

trong đó D là phương sai còn N là số sự kiện dùng để mô phỏng. Theo công thức này ta thấy nêu lấy mẫu thông kê càng lớn thì kết quả tính toán Monte Carlo càng chính xác. Chương trình MCNP MCNP (Monte Carlo N-Partuicles) là một chương trình máy tính ứng dụng phương pháp Monte Carlo dé mô phóng quá trình vận chuyên của các hat neutron, photon và electron riêng biệt. hoặc kết hợp sự vận chuyển của neutron/photon (trong đó các photon được tạo ra bởi tương tác của neutron với vật liệu), neutron/photon/ electron, photon/electron, electron/photon trong những môi trường vật liệu khác nhau.

Chương trình MCNP được phát triển bởi Trung tâm thí nghiệm quốc gia Los Alamos = My (Los Alamos National Laboratory = USA). Ké tir khi các phiên ban đầu tiên của MCNP được đưa vào ứng dung trong thập niên 1980, các nhà lập trình MCNP đã không ngừng nâng cấp và cho ra đời những phiên bản mới hơn trong một thời gian ngắn cùng với sự hoàn thiện hơn về các quá trình vật lý của các hạt, các thư viện ứng dụng và những tính năng tiện ích của chúng. Chương trình MCNP sử dụng các thư viện số liệu hạt nhân và nguyên tử năng lượng liên tục. Các nguồn số liệu hạt nhân chủ yếu là các đánh giá từ các số liệu hạt nhân ENDE, thư viện các số liệu hạt nhân ENDL, các thư viện kích hoạt AC TL từ Livemore và các đánh giá từ nhóm khoa học hạt nhân ứng dụng ở Los Alamos [9].

Các bảng tương tác gamma được xây dựng đối với tat cả các nguyên tổ từ Z = 1 đến Z = 94. Các số liệu trong các bảng tương tác photon cho phép MCNP tính toán tán xạ kết hợp, không kết hợp, hấp thụ quang điện với khả năng phát huỳnh quang và tạo cặp. Trong luận văn này, phiên bản MCNPS được sử dụng đẻ mô phỏng bài toán vận chuyên của photon trong mô hình hệ đo tan xạ ngược nhằm thu được các phé tan xa. Mô hình tương tác của photon với vật chat trong MCNP Chương trình MCNP giải bài toán va chạm của các photon (xem như là hạt) với vật chất qua hai mô hình: đơn giản và chi tiết, dựa trên lý thuyết của ba loại tương tác là tán xạ Compton, hiệu ứng quang điện và hiệu ứng tạo cặp [2].

Sau đây, cơ sở lý thuyết của sự mô phỏng ba loại tương tác theo mô hình chỉ tiết sẽ được trình bày: 1. Tấn xạ Compton Tán xạ Compton lả quả trình tương tác của photon với electron liên kết yếu của nguyên tử, trong đó photon truyền một phan năng lượng cho electron và bị lệch đi so với hướng ban đầu (hình 1. Do năng lượng của photon lớn hơn rất nhiều so với năng lượng liên kết của electron trong nguyên tử nên electron được xem là electron tự đo. Từ định luật bảo toàn năng lượng va xung lượng ta có mỗi liên hệ giữa năng lượng sau tan xạ E, góc tan xạ 6 và năng lượng ban dau Ea, như sau: Eq 1+—*, (1-cos@) m.c* Năng lượng truyền cho electron: E.

hv, ; với Ky= mc’ E,= hy, và góc tán xạ của electron được tinh theo công thức: tang = cot! (1.3) ltkẹ 2 , * Đối với góc tán xạ rất nhỏ (0=0"), năng lượng electron Egnin = 0, khi đó photon tán xạ có năng lượng gần bằng với năng lượng của photon tới. * Dối với góc tán xạ lớn (8 =180”), electron giật lùi có năng lượng cực đại và bằng: E. Hiệu ứng quang điện Hiệu ứng quang điện là quá trình tương tác giữa photon và các electron liên kết trong nguyên tử (hình 1. Trong quá trình tương tác này, photon truyền toàn bộ năng lượng của minh cho electron liên kết của nguyên tử, năng lượng này đủ đẻ birt các electron ra khỏi lớp vỏ nguyên tử, các electron này được gọi là electron quang điện (photoelectron).

Một phần nhỏ năng lượng này giúp electron thắng năng lượng liên kết giữa electron với nguyên tử dé đảm bảo định luật bảo toàn năng lượng và xung lượng. Hiệu ứng quang điện chỉ xảy ra khi năng lượng photon tới lớn hơn năng lượng liên kết của clectron trong nguyên tử. Trong hiệu ứng quang điện. năng lượng giật lùi của nguyên tử xem như không đáng kê, do đó động năng của clectron được xác định: E.5) với E, là năng lượng liên kết của electron thứ i (i = K, L, M,.

bức xạ đặc trưng Hình 1. Minh họa hiệu ứng quang điện Hiệu ứng quang điện là một kênh trội của tương tác giữa photon với vật chất ở vùng năng lượng tương đối thấp. Tiết diện của hiệu ứng quang điện phụ thuộc chủ yếu vào năng lượng của photon tới và điện tích Z của hạt nhân nguyên tử. Đối với những vật liệu nặng (Z lớn) thì xác suất xảy ra hiệu ứng quang điện lớn ngay cả với những photon có năng lượng cao.

đối với những vật liệu nhẹ thì hiệu ứng quang điện chi có ý nghĩa với những photon năng lượng thấp. Biểu thức thé hiện tiết điện quang điện ty lệ theo hệ thức sau: G9 z* (1.