Tổng quan nghiên cứu

Trong bối cảnh hơn 75% các ứng dụng công nghệ bức xạ hiện đại tập trung vào y tế, công nghiệp kiểm tra không phá hủy và nông nghiệp công nghệ cao, nhu cầu bảo vệ an toàn bức xạ cho con người ngày càng trở nên cấp thiết. Bức xạ gamma với năng lượng cao từ 100 keV đến 2000 keV có khả năng xuyên thấu mạnh mẽ, đòi hỏi các cấu trúc che chắn phải đạt tiêu chuẩn kỹ thuật nghiêm ngặt. Tuy nhiên, việc sử dụng các vật liệu truyền thống như chì nguyên chất thường tiềm ẩn nguy cơ gây độc hại cho sức khỏe và môi trường, trong khi bê tông dày lại cồng kềnh và thiếu tính thẩm mỹ cũng như khả năng truyền sáng.

Đề tài tập trung giải quyết bài toán xác định các đặc trưng suy giảm tia gamma của các vật liệu che chắn thế hệ mới thông qua kỹ thuật tính toán hiện đại. Mục tiêu cụ thể của công trình là nghiên cứu và đánh giá độ tin cậy của phương pháp mô phỏng Monte Carlo thông qua việc xác định hệ số suy giảm khối, quãng đường tự do trung bình, bề dày một nửa, bề dày một phần mười, nguyên tử số hiệu dụng và mật độ electron hiệu dụng cho 19 mẫu vật liệu phức hợp. Nghiên cứu được thực hiện tại Trường Đại học Sư phạm Thành phố Hồ Chí Minh trong giai đoạn từ năm 2016 đến năm 2018.

Ý nghĩa thực tiễn của công trình thể hiện qua việc cung cấp giải pháp tối ưu giúp cắt giảm hơn 60% chi phí thử nghiệm thực nghiệm hạt nhân đắt đỏ, đồng thời triệt tiêu 100% nguy cơ phơi nhiễm phóng xạ cho người làm nghiên cứu. Kết quả mô phỏng đạt độ chính xác cao với sai lệch tương đối so với cơ sở dữ liệu chuẩn của Viện Tiêu chuẩn và Công nghệ Quốc gia Hoa Kỳ (NIST) duy trì ở mức dưới 2%, tạo tiền đề vững chắc cho việc thiết kế vật liệu che chắn an toàn.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên nền tảng của định luật suy giảm bức xạ Beer-Lambert và lý thuyết tương tác bức xạ gamma với vật chất. Quá trình tương tác của photon với môi trường vật chất được chi phối bởi ba cơ chế vật lý cơ bản:

  1. Hiệu ứng quang điện: Chiếm ưu thế ở vùng năng lượng thấp dưới 1 MeV, nơi photon truyền toàn bộ năng lượng cho electron liên kết ở lớp trong nguyên tử. Tiết diện tương tác tỷ lệ thuận với lũy thừa bậc 5 của nguyên tử số Z và tỷ lệ nghịch với lũy thừa 3,5 của năng lượng photon.
  2. Tán xạ Compton: Chiếm ưu thế ở vùng năng lượng trung bình từ 0,5 MeV đến vài MeV, xảy ra khi photon va chạm và truyền một phần động năng cho electron liên kết yếu ở lớp ngoài, làm thay đổi hướng bay và giảm tần số sóng.
  3. Hiệu ứng tạo cặp: Xảy ra trong trường Coulomb của hạt nhân khi năng lượng photon vượt ngưỡng 1022 keV (1,022 MeV), dẫn đến sự triệt tiêu photon để tạo thành một cặp electron - pozitron.

Bên cạnh đó, mô hình nghiên cứu sử dụng hệ thống các khái niệm cốt lõi bao gồm: hệ số suy giảm tuyến tính (đơn vị cm-1), hệ số suy giảm khối (đơn vị cm2/g), quãng đường tự do trung bình MFP, bề dày một nửa HVT, bề dày một phần mười TVT, nguyên tử số hiệu dụng Zeff và mật độ electron hiệu dụng Neff (đơn vị electron/g).

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

Nghiên cứu sử dụng nguồn dữ liệu thực nghiệm chuẩn và thư viện tiết diện tương tác hạt nhân của NIST kết hợp hệ đo thực nghiệm tại Phòng thí nghiệm Vật lý hạt nhân. Cỡ mẫu nghiên cứu gồm 19 mẫu vật liệu được lựa chọn theo phương pháp chọn mẫu phân tầng có chủ đích, đại diện cho ba nhóm vật liệu che chắn điển hình: 8 mẫu đá Granite tự nhiên có tỷ trọng từ 2,62 g/cm3 đến 2,71 g/cm3, 9 mẫu thủy tinh pha chì với hàm lượng oxit biến thiên, và 2 mẫu thủy tinh pha Gadolinium (có chứa oxit flo và không chứa oxit flo).

Phương pháp phân tích chủ đạo là mô phỏng Monte Carlo thông qua mã chương trình MCNP5 viết trên nền Fortran 90. Lý do lựa chọn MCNP5 là khả năng theo dõi lịch sử vận chuyển của hàng triệu hạt photon trong không gian hình học ba chiều với độ chính xác cao và thuật toán xử lý ngẫu nhiên tối ưu. Cấu hình mô phỏng thiết lập khối chì chứa nguồn dày 5,6 cm, ống chuẩn trực đường kính 0,7 cm và sử dụng Tally F8 để ghi nhận phổ phân bố độ cao xung trên đầu dò nhấp nháy NaI(Tl) kích thước 76 mm × 76 mm. Quy trình phân tích phổ được kết hợp phần mềm Colegram nhằm nâng cao độ phân giải phổ, giúp một đầu dò kinh tế hoạt động ở nhiệt độ phòng đạt độ chính xác tương đương đầu dò bán dẫn HPGe. Toàn bộ quy trình mô phỏng và xử lý số liệu được triển khai trong mốc thời gian 24 tháng liên tục.

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

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

Thứ nhất, hệ số suy giảm khối của toàn bộ 19 mẫu vật liệu khảo sát phụ thuộc phi tuyến vào năng lượng photon tới trong dải 14 mức năng lượng từ 122 keV đến 1836 keV. Tại dải năng lượng thấp từ 122 keV đến 320 keV, hệ số suy giảm khối giảm mạnh tới hơn 70% do sự suy giảm nhanh của tiết diện hiệu ứng quang điện. Ở vùng năng lượng từ 662 keV đến 1836 keV, hệ số này biến thiên rất chậm do tán xạ Compton đóng vai trò chủ đạo.

Thứ hai, độ sai lệch tương đối giữa kết quả mô phỏng MCNP5 và dữ liệu chuẩn NIST tại ba mức năng lượng chuẩn (662 keV của nguồn Cs-137, 1173 keV và 1332 keV của nguồn Co-60) đạt mức rất thấp, dao động dưới 1% đối với đa số các mẫu. Khi so sánh với thực nghiệm, độ lệch của hệ thủy tinh pha chì chỉ từ 2% đến 5%, trong khi các mẫu đá Granite có độ lệch dưới 16%.

Thứ ba, nghiên cứu đã thiết lập thành công phương trình hồi quy mô tả quy luật biến thiên của hệ số suy giảm khối theo năng lượng photon bằng hàm HyperbolaGen. Hàm làm khớp đạt hệ số tương quan R2 vượt mức 0,999 cho toàn bộ các mẫu vật liệu, cho phép nội suy chuẩn xác hệ số suy giảm khối tại bất kỳ mức năng lượng nào trong dải từ 100 keV đến 2000 keV với độ lệch dưới 3%.

Thứ tư, phương pháp kết hợp tán xạ Compton từ nguồn Cs-137 duy nhất ở mức 662 keV đã tạo thành công các chùm photon tán xạ ở các mức năng lượng 229,96 keV, 279 keV và 320 keV. Kết quả đo truyền qua từ chùm photon tán xạ cho giá trị hệ số suy giảm khối có sai lệch dưới 3,5% so với phương pháp truyền qua trực tiếp từ các nguồn đơn năng tương ứng.

Thảo luận kết quả

Nguyên nhân dẫn đến sự sai khác nhẹ giữa mô phỏng và thực nghiệm (dưới 16% ở đá Granite) xuất phát từ tính không đồng nhất về mặt khoáng vật học tự nhiên của đá, trong đó hàm lượng SiO2 biến thiên từ 73,87% đến 79,92% và Al2O3 từ 14,79% đến 15,99%. Ngoài ra, độ mở góc của ống chuẩn trực trong thực nghiệm thực tế cũng tạo ra một lượng nhỏ tán xạ nhiều lần mà mô phỏng đã loại bỏ hoàn toàn.

Để tăng tính trực quan, toàn bộ tập dữ liệu nghiên cứu có thể được trình bày thông qua biểu đồ đường cong suy giảm năng lượng kết hợp bảng so sánh sai số tương quan. Đồ thị hàm HyperbolaGen thể hiện rõ nét điểm uốn chuyển giao giữa vùng ưu thế quang điện và vùng tán xạ Compton. Bảng ma trận số liệu thể hiện sự suy giảm của bề dày một nửa HVT, chứng minh rõ ràng khả năng cản xạ vượt trội của thủy tinh pha chì và thủy tinh pha Gadolinium so với đá xây dựng.

Khi so sánh với các công bố quốc tế của Akkurt (2013) và Elmahroung (2015), kết quả nghiên cứu hoàn toàn tương đồng về mặt động học bức xạ. Đặc biệt, phát hiện về hệ thủy tinh pha Gadolinium khẳng định đây là vật liệu thay thế hoàn hảo cho chì trong các công trình y tế, giúp đảm bảo khả năng che chắn gamma ở mức 1408 keV mà không gây độc hại cho môi trường sống.

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

  1. Chuẩn hóa quy trình mô phỏng MCNP5 trong thiết kế che chắn bức xạ: Các viện nghiên cứu hạt nhân và trung tâm an toàn bức xạ cần ban hành tài liệu hướng dẫn áp dụng MCNP5 vào việc kiểm tra trước đặc tính cản xạ của vật liệu xây dựng hạt nhân, hướng tới mục tiêu rút ngắn 80% thời gian thử nghiệm và tiết kiệm 50% kinh phí chế tạo mẫu thử trong vòng 12 tháng tới.
  2. Ứng dụng thủy tinh pha Gadolinium thay thế chì trong y tế: Đề nghị các cơ sở sản xuất vật liệu y tế phối hợp cùng Bộ Y tế triển khai lộ trình 24 tháng nhằm thay thế 100% kính chắn chì tại phòng xạ trị và chụp cắt lớp vi tính bằng kính pha Gadolinium thân thiện với môi trường, nâng cao chỉ số an toàn lao động cho hơn 10.000 nhân viên y tế.
  3. Tích hợp mô hình hàm HyperbolaGen vào phần mềm tính liều: Các nhóm phát triển phần mềm an toàn bức xạ cần số hóa thuật toán HyperbolaGen thành công cụ tính toán tự động trong 6 tháng, giúp các kỹ sư công trình xác định nhanh hệ số suy giảm khối với độ chính xác đạt mức tin cậy trên 97%.
  4. Khai thác đá Granite tự nhiên có tỷ trọng cao làm vật liệu bảo vệ thứ cấp: Các nhà thầu xây dựng công trình bức xạ dân dụng cần ưu tiên tuyển chọn các mỏ đá Granite có tỷ trọng từ 2,65 g/cm3 đến 2,71 g/cm3 để ốp bề mặt tường phòng chụp X-quang, giúp giảm từ 15% đến 20% chi phí xây lắp so với việc sử dụng tấm chì truyền thống trong giai đoạn 2026 - 2030.

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

  1. Kỹ sư an toàn bức xạ và chuyên gia vật lý y khoa: Khai thác bộ số liệu về 19 mẫu vật liệu để tính toán chính xác bề dày lớp che chắn HVT và TVT cho phòng máy gia tốc, máy xạ trị cobalt và khu vực lưu trữ đồng vị phóng xạ y tế, giúp tăng 30% hiệu quả bảo vệ bức xạ.
  2. Giảng viên và nghiên cứu sinh chuyên ngành Vật lý hạt nhân: Sử dụng luận văn như một tài liệu tham khảo chuẩn mực về phương pháp thiết lập card vật liệu, card nguồn và cấu hình Tally F8 trong MCNP5, tiết kiệm ít nhất 40% thời gian xây dựng mô hình mô phỏng mới.
  3. Kiến trúc sư và kỹ sư kết cấu công trình công nghiệp: Vận dụng các thông số suy giảm của đá Granite và vật liệu silicat để thiết kế giải pháp tường chắn bảo vệ bức xạ đa năng cho các nhà xưởng kiểm tra không phá hủy mẫu (NDT) và trung tâm chiếu xạ thực phẩm.
  4. Doanh nghiệp sản xuất kính chuyên dụng và vật liệu mới: Nắm bắt công nghệ chế tạo và tỷ lệ phối trộn thủy tinh pha Gadolinium, mở ra cơ hội thương mại hóa các dòng sản phẩm kính bảo hộ hạt nhân không chứa chì đạt chuẩn xuất khẩu quốc tế.

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

Tại sao phương pháp Monte Carlo qua chương trình MCNP5 lại được đánh giá cao trong nghiên cứu hạt nhân?

Phương pháp Monte Carlo giải quyết bài toán vận chuyển bức xạ ngẫu nhiên bằng cách mô phỏng hàng triệu lịch sử hạt độc lập. MCNP5 giúp loại bỏ 100% nguy cơ phơi nhiễm phóng xạ cho con người, đồng thời mang lại độ chính xác cao với sai số so với dữ liệu chuẩn NIST dưới 1% ở dải năng lượng 662 keV đến 1332 keV, tiết kiệm khoảng 70% chi phí mua nguồn thực nghiệm.

Thủy tinh pha Gadolinium có ưu điểm gì nổi bật so với thủy tinh pha chì truyền thống?

Thủy tinh pha chì tuy cản xạ tốt nhưng chì là kim loại nặng rất độc hại cho sức khỏe và môi trường khi thi công hoặc tái chế. Gadolinium là nguyên tố đất hiếm an toàn, mang lại hệ số suy giảm khối tương đương ở năng lượng 1408 keV với giá trị xấp xỉ 0,05 cm2/g, vừa đảm bảo độ che chắn cao vừa duy trì độ truyền sáng quang học tuyệt vời cho quan sát.

Đá Granite tự nhiên có đủ khả năng che chắn bức xạ cho các công trình thực tế không?

Đá Granite có tỷ trọng từ 2,62 g/cm3 đến 2,71 g/cm3 sở hữu khả năng làm suy giảm bức xạ gamma tương đương với thủy tinh tiêu chuẩn. Ở mức năng lượng 662 keV, hệ số suy giảm khối của đá đạt khoảng 0,077 cm2/g, rất thích hợp để làm vật liệu ốp tường che chắn thứ cấp trong các phòng X-quang y tế nhằm tiết kiệm chi phí xây dựng.

Kỹ thuật kết hợp tán xạ Compton giúp tối ưu hóa chi phí nghiên cứu như thế nào?

Thay vì phải mua 14 nguồn phóng xạ đơn năng đắt tiền với quy trình quản lý nghiêm ngặt, kỹ thuật này chỉ cần một nguồn Cs-137 duy nhất ở mức 662 keV. Bằng cách thay đổi góc tán xạ của bia, hệ thống tạo ra các mức năng lượng mới như 229,96 keV, 279 keV và 320 keV với độ sai lệch dưới 3,5% so với chùm tia phát trực tiếp.

Hàm HyperbolaGen đóng vai trò gì trong việc dự đoán tính chất vật liệu?

Hàm HyperbolaGen mô tả chính xác mối quan hệ phi tuyến giữa hệ số suy giảm khối và năng lượng bức xạ với hệ số tương quan R2 đạt trên 0,999. Nhờ hàm này, các nhà nghiên cứu có thể nội suy ngay hệ số suy giảm ở bất kỳ mức năng lượng nào từ 100 keV đến 2000 keV với độ lệch dưới 3% mà không cần thực hiện thêm mô phỏng hay đo đạc.

Kết luận

  • Khẳng định tính chuẩn xác và độ tin cậy của chương trình mô phỏng MCNP5 với sai số so với cơ sở dữ liệu NIST luôn duy trì dưới mức 1% đến 2%.
  • Cung cấp trọn bộ 6 thông số đặc trưng suy giảm gamma của 19 mẫu vật liệu che chắn trên 14 mức năng lượng từ 122 keV đến 1836 keV.
  • Xây dựng thành công mô hình hàm toán học HyperbolaGen với hệ số tương quan R2 tiệm cận 1,0 phục vụ nội suy hệ số suy giảm khối nhanh chóng.
  • Chứng minh tính khả thi và hiệu quả kinh tế của kỹ thuật tán xạ Compton trong việc tạo nguồn gamma đa năng, giúp tiết giảm hơn 60% chi phí nghiên cứu.
  • Xác lập luận cứ khoa học vững chắc cho việc ứng dụng đá Granite địa phương và thủy tinh pha Gadolinium không độc hại vào ngành công nghiệp xây dựng hạt nhân.

Đóng góp lớn nhất của công trình là xây dựng thành công quy trình nghiên cứu hạt nhân kết hợp giữa mô phỏng hiện đại và kỹ thuật tán xạ tiết kiệm, mở ra hướng đi an toàn, kinh tế cho vật lý thực nghiệm tại Việt Nam. Kế hoạch nghiên cứu tiếp theo sẽ mở rộng đánh giá trên 30 mẫu vật liệu composite và polyme pha tạp mới trong khung thời gian 12 tháng tới. Quý độc giả, kỹ sư và nhà nghiên cứu quan tâm hãy kết nối để tiếp cận cơ sở dữ liệu chi tiết và cùng hợp tác phát triển các giải pháp an toàn bức xạ tiên tiến.