Luận văn thạc sĩ: Khảo sát và mô phỏng tính toán che chắn bức xạ gamma

Khảo sát và mô phỏng tính toán cấu hình che chắn nguồn bức xạ gamma bằng phương pháp Monte Carlo trong luận văn thạc sĩ vật lý kỹ thuật.

• Tác giả: Huỳnh Thụy Vi
• Người hướng dẫn: TS. Huỳnh Quang Linh
• Trường học: Đại học Bách Khoa - ĐHQG TP.HCM
• Chuyên ngành: Vật lý kỹ thuật
• Đề tài: Khảo Sát Và Mô Phỏng Tính Toán Các Cấu Hình Che Chắn Của Một Số Nguồn Bức Xạ Gamma Bằng Phương Pháp Monte Carlo
• Loại tài liệu: luận văn thạc sĩ
• Năm xuất bản: 2013
• Địa điểm: TP. Hồ Chí Minh

Mục lục chi tiết

• DANH MỤC CÁC BẢNG
• DANH MỤC CÁC HÌNH VẼ
• LỜI MỞ ĐẦU – MỤC TIÊU VÀ NHIỆM VỤ LUẬN VĂN
• 1. CHƯƠNG I: TƯƠNG TÁC CỦA BỨC XẠ GAMMA VỚI VẬT CHẤT
• 1.1. Bản chất của tia gamma
• 1.2. Tính chất của tia gamma
• 1.3. Tương tác của tia gamma với vật chất
• 1.3.1. Hiệu ứng quang điện
• 1.3.2. Hiệu ứng Compton
• 1.3.3. Hiệu ứng tạo cặp
• 1.4. Sự suy giảm của chùm bức xạ gamma khi đi qua vật chất
• 1.4.1. Sự suy giảm của chùm gamma hẹp
• 1.4.2. Sự suy giảm của chùm gamma rộng. Hệ số tích luỹ
• 1.4.3. Hệ số hấp thụ tuyến tính
• 1.4.4. Hệ số suy giảm khối
• 1.5. Các đơn vị, khái niệm trong an toàn phóng xạ
• 1.5.1. Liều hấp thụ
• 1.5.2. Suất liều hấp thụ
• 1.5.3. Suất liều chiếu
• 1.5.4. Liều tương đương
• 1.5.5. Liều hiệu dụng
• 1.5.6. Mối tương quan giữa các đơn vị đo bức xạ
• 1.5.7. Mức chiếu xạ được phép giới hạn
• 2. CHƯƠNG II: SỰ CHE CHẮN BỨC XẠ GAMMA
• 2.1. Tổng quan những nghiên cứu che chắn bức xạ từ trước đến nay
• 2.2. Kỹ thuật bảo vệ chống bức xạ
• 2.2.1. Thời gian chiếu xạ
• 2.2.2. Khoảng cách tới nguồn
• 2.2.3. Nguồn có kích thước
• 2.2.4. Che chắn bức xạ
• 2.2.4.1. Bảo vệ sơ cấp
• 2.2.4.2. Bảo vệ thứ cấp
• 2.3. Các tính toán bảo vệ chống bức xạ gamma (theo độ bội quy giảm) có chú ý tới hệ số tích lũy
• 2.3.1. Sự suy giảm cường độ của chùm bức xạ. Hệ số tích lũy
• 2.3.2. Các loại hệ số tích lũy liên quan đến đặc trưng của bức xạ tới
• 2.3.3. Các loại hệ số tích lũy liên quan đến đặc trưng của môi trường
• 2.3.4. Đối với môi trường che chắn có nhiều lớp
• 2.3.5. Độ bội suy giảm
• 2.3.6. Tính toán bề dày lớp bảo vệ theo các lớp suy giảm một nửa
• 2.3.6.1. Đối với nguồn gamma đơn năng
• 2.3.6.2. Đối với nguồn gamma đa năng
• 2.4. Che chắn nguồn gamma có dạng hình học khác nhau
• 2.4.1. Che chắn nguồn bức xạ gamma dạng hình học điểm
• 2.4.2. Che chắn nguồn gamma có dạng hình học bất kì
• 2.4.3. Che chắn nguồn gamma có dạng hình học dây thẳng
• 2.4.4. Che chắn nguồn bức xạ gamma dạng hình tròn
• 2.4.5. Che chắn nguồn gamma dạng hình học cầu rỗng
• 2.4.6. Che chắn nguồn gamma dạng hình học cầu đặc
• 2.4.7. Suất liều chiếu khi nguồn mặt trụ
• 3. CHƯƠNG III: THỰC HÀNH TÍNH TOÁN THIẾT KẾ CHE CHẮN BỨC XẠ GAMMA
• 3.1. Phương pháp giải tích
• 3.1.1. Nguồn bức xạ gamma dạng hình học điểm
• 3.1.2. Nguồn bức xạ gamma dạng trụ đặc
• 3.2. Phương pháp Monte Carlo - Chương trình Penelope
• 3.2.1. Chương trình mô phỏng PENELOPE
• 3.2.2. Tính suất liều của nguồn bức xạ dạng hình học điểm và nguồn hình trụ đặc
• 3.3. Sự phụ thuộc của liều hấp thụ vào góc tán xạ
• TÀI LIỆU THAM KHẢO
• Phụ lục 1: Các hằng số của các nguồn phóng xạ
• Phụ lục 2: Các hằng số A1, α1, α2, D, μ của 1 số vật liệu che chắn

Tóm tắt

I. Tổng quan về bức xạ gamma

Bức xạ gamma là một dạng bức xạ ion hóa có năng lượng cao, phát ra từ sự phân rã hạt nhân. Tia gamma có khả năng xuyên thấu mạnh mẽ, chỉ có thể bị chặn lại bởi các vật liệu có nguyên tử lượng lớn như chì hoặc bê tông. Các quá trình tương tác chính của bức xạ gamma với vật chất bao gồm hiệu ứng quang điện, hiệu ứng Compton và sự sinh cặp. Những quá trình này không chỉ ảnh hưởng đến sự suy giảm cường độ bức xạ mà còn có thể gây ra các tác động nguy hiểm đến sức khỏe con người. Việc hiểu rõ về bức xạ gamma và các quá trình tương tác của nó với vật chất là rất quan trọng trong việc thiết kế các cấu hình che chắn hiệu quả.

1.1. Tính chất của tia gamma

Tia gamma không thể nhìn thấy bằng mắt thường và không cảm nhận được bằng các giác quan. Chúng có khả năng ion hóa vật chất, đặc biệt là chất khí, và có thể gây hại cho tế bào sống. Độ xuyên thấu của tia gamma phụ thuộc vào năng lượng của photon và tính chất của vật chất mà chúng đi qua. Các định luật cơ bản của ánh sáng cũng áp dụng cho tia gamma, bao gồm phản xạ và khúc xạ. Việc nắm vững các tính chất này là cần thiết để phát triển các phương pháp che chắn bức xạ hiệu quả.

II. Kỹ thuật che chắn bức xạ gamma

Che chắn bức xạ gamma là một kỹ thuật quan trọng nhằm bảo vệ con người và môi trường khỏi tác hại của bức xạ ion hóa. Các phương pháp che chắn thường được áp dụng bao gồm sử dụng vật liệu có mật độ cao và thiết kế cấu hình che chắn hợp lý. Việc tính toán độ dày lớp che chắn cần thiết dựa trên các thông số như cường độ bức xạ, khoảng cách từ nguồn đến người tiếp xúc và thời gian chiếu xạ. Kỹ thuật che chắn không chỉ giúp giảm thiểu liều lượng bức xạ mà còn đảm bảo an toàn trong các hoạt động liên quan đến bức xạ.

2.1. Các phương pháp tính toán che chắn

Các phương pháp tính toán che chắn bức xạ gamma thường được chia thành hai loại: phương pháp giải tích và phương pháp mô phỏng. Phương pháp giải tích thường được sử dụng cho các cấu hình hình học đơn giản, trong khi phương pháp mô phỏng, đặc biệt là phương pháp Monte Carlo, cho phép xử lý các cấu hình phức tạp hơn. Chương trình PENELOPE 2006 là một trong những công cụ mạnh mẽ để mô phỏng quá trình lan truyền bức xạ trong vật chất, giúp kiểm tra độ tin cậy của các phương pháp tính toán giải tích.

III. Ứng dụng của phương pháp Monte Carlo trong che chắn bức xạ

Phương pháp Monte Carlo đã trở thành một công cụ quan trọng trong nghiên cứu và thiết kế che chắn bức xạ gamma. Bằng cách mô phỏng quá trình tương tác của bức xạ với vật chất, phương pháp này cho phép đánh giá chính xác hơn về hiệu quả của các cấu hình che chắn. Việc sử dụng chương trình PENELOPE 2006 để mô phỏng các cấu hình che chắn giúp các nhà nghiên cứu có được những kết quả đáng tin cậy, từ đó đưa ra các giải pháp tối ưu cho việc bảo vệ an toàn bức xạ.

3.1. Kết quả mô phỏng và ứng dụng thực tiễn

Kết quả từ các mô phỏng cho thấy rằng việc áp dụng phương pháp Monte Carlo có thể cải thiện đáng kể độ chính xác trong việc tính toán liều lượng bức xạ. Các ứng dụng thực tiễn của phương pháp này không chỉ giới hạn trong lĩnh vực y tế mà còn mở rộng ra các ngành công nghiệp khác như năng lượng hạt nhân, nghiên cứu môi trường và an toàn lao động. Việc phát triển kỹ năng sử dụng phương pháp Monte Carlo trong nghiên cứu che chắn bức xạ là cần thiết để nâng cao hiệu quả và độ tin cậy trong các ứng dụng thực tế.

Nội dung chính

Tổng quan nghiên cứu

Kỹ thuật bức xạ hạt nhân, đặc biệt là bức xạ photon năng lượng cao như tia gamma từ 50 keV đến trên 10 MeV, đóng vai trò then chốt trong y học xạ trị, kiểm tra không phá hủy và nông nghiệp công nghệ cao. Tuy nhiên, khả năng ion hóa mạnh mẽ của tia gamma tiềm ẩn nhiều hiểm họa sinh học, đòi hỏi các tiêu chuẩn an toàn nghiêm ngặt theo khuyến cáo của Ủy ban Bảo vệ Phóng xạ Quốc tế với mức giới hạn liều hiệu dụng cho nhân viên bức xạ là 20 mSv/năm và công chúng là 1 mSv/năm. Vấn đề trọng tâm trong thực tế là việc tính toán thiết kế tường che chắn bức xạ thường dựa vào phương pháp giải tích truyền thống. Phương pháp này chỉ phù hợp với các dạng hình học đơn giản và thường dẫn đến sai số ước lượng đáng kể hoặc lãng phí chi phí vật liệu do tăng bề dày bảo vệ quá mức cần thiết.

Nghiên cứu này được thực hiện nhằm khảo sát toàn diện các phương pháp tính toán giải tích kỹ thuật che chắn bức xạ gamma, đồng thời ứng dụng kỹ thuật mô phỏng Monte Carlo thông qua mã chương trình chuyên dụng để kiểm chứng và tối ưu hóa các cấu hình che chắn. Phạm vi nghiên cứu tập trung vào các nguồn phóng xạ phổ biến như Coban-60 và Xesi-137 với các dạng hình học cơ bản gồm nguồn điểm, nguồn dây, nguồn cầu và nguồn hình trụ đặc. Công trình được hoàn thành tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh vào tháng 7 năm 2013. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn trong việc xây dựng chuẩn dữ liệu tính toán, hỗ trợ giảm thiểu từ 15% đến 25% chi phí vật tư che chắn bằng chì, sắt hoặc bê tông nhưng vẫn đảm bảo duy trì suất liều khu vực kiểm soát dưới ngưỡng an toàn 8 microSievert/giờ.

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 xây dựng trên nền tảng lý thuyết tương tác của photon năng lượng cao với môi trường vật chất. Quá trình truyền bức xạ gamma qua vật chất chịu sự chi phối của ba hiệu ứng vật lý cơ bản: hiệu ứng quang điện chiếm ưu thế ở vùng năng lượng thấp và phụ thuộc bậc ba vào năng lượng photon; hiệu ứng tán xạ Compton chiếm ưu thế trong khoảng năng lượng từ 0.1 MeV đến 5.0 MeV; và hiệu ứng tạo cặp electron-positron diễn ra khi năng lượng photon vượt ngưỡng 1.022 MeV. Sự suy giảm cường độ chùm photon hẹp tuân theo định luật hàm mũ Lambert-Beer thông qua hệ số hấp thụ tuyến tính tính bằng cm^-1 và hệ số suy giảm khối tính bằng cm^2/g.

Đối với chùm bức xạ rộng thường gặp trong thực tế, hiện tượng tán xạ nhiều lần tạo ra bức xạ thứ cấp làm tăng cường độ đo được tại đầu dò. Do đó, lý thuyết hệ số tích lũy liều và hệ số tích lũy năng lượng được tích hợp vào phương trình suy giảm chùm tia. Khung liều lượng học áp dụng các đại lượng chuẩn hóa quốc tế: liều hấp thụ đo bằng Gray (1 Gy = 100 rad), liều tương đương đo bằng Sievert với trọng số bức xạ bằng 1 đối với photon gamma, và liều hiệu dụng tích hợp các trọng số mô như 0.20 cho tuyến sinh dục và 0.12 cho phổi, dạ dày và tủy xương. Khái niệm bề dày suy giảm một nửa và bề dày suy giảm một phần mười được sử dụng làm thước đo định lượng cho khả năng cản xạ của vật liệu.

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

Nghiên cứu sử dụng hai cách tiếp cận song song là phương pháp giải tích và phương pháp mô phỏng ngẫu nhiên Monte Carlo. Về mặt giải tích, tác giả thiết lập các tích phân giải tích kết hợp công thức thực nghiệm Taylor và Berger để xác định suất liều cho các cấu hình nguồn điểm, nguồn dây thẳng, nguồn vòng tròn và nguồn trụ đặc được che chắn bởi các tấm vật liệu phẳng hoặc đồng tâm.

Về phương pháp mô phỏng, nghiên cứu ứng dụng chương trình PENELOPE 2006, một công cụ tính toán vận chuyển photon và electron tiên tiến. Lý do lựa chọn PENELOPE là khả năng xử lý chi tiết từng tương tác vật lý ngẫu nhiên của photon trong dải năng lượng từ 50 eV đến 1 GeV với độ chính xác cao. Cỡ mẫu mô phỏng được thiết lập từ 10^6 đến 10^7 lịch sử hạt cho mỗi cấu hình hình học, áp dụng phương pháp chọn mẫu phân tầng và kỹ thuật giảm phương sai để đảm bảo sai số thống kê của suất liều luôn được kiểm soát dưới 1.5%. Các dữ liệu đầu vào bao gồm phổ năng lượng đặc trưng của Coban-60 (1.173 MeV và 1.332 MeV) và Xesi-137 (661.66 keV), cùng các thông số mật độ vật liệu chuẩn như chì (11.34 g/cm3), sắt (7.87 g/cm3) và bê tông (2.35 g/cm3). Quá trình khảo sát và thu thập dữ liệu được thực hiện liên tục trong mười hai tháng nghiên cứu.

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

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

Thứ nhất, đối với nguồn bức xạ Coban-60 dạng điểm che chắn bằng tấm chì, giá trị bề dày suy giảm một nửa tính toán được là 1.25 cm và bề dày suy giảm một phần mười là 4.12 cm. Khi so sánh giữa phương pháp giải tích chùm hẹp và mô phỏng Monte Carlo, độ lệch tương đối chỉ ở mức 3.8%, xác nhận độ tin cậy của mô hình mô phỏng khi số lượng tán xạ thứ cấp ở mức tối thiểu.

Thứ hai, trong điều kiện hình học chùm rộng với bề dày lớp chắn chì vượt quá 5.0 cm, hiện tượng tích lũy bức xạ tán xạ Compton làm tăng suất liều thực tế từ 25% đến 65% so với mô hình giải tích thuần túy không xét hệ số tích lũy. Việc áp dụng hệ số tích lũy thực nghiệm dạng Berger giúp kéo giảm sai số này xuống dưới 6.2% so với kết quả mô phỏng PENELOPE.

Thứ ba, khảo sát trên nguồn phóng xạ dạng hình trụ đặc cho thấy sự phân bố suất liều bất đối xứng rõ rệt giữa phương che chắn dọc trục và phương che chắn ngang trục. Do hiệu ứng tự hấp thụ bức xạ bên trong lòng nguồn, suất liều cực đại theo phương ngang giảm 18.4% so với tính toán lý thuyết xem nguồn là tập hợp các điểm không tự hấp thụ.

Thứ tư, khi đánh giá lớp bảo vệ thứ cấp cho tia tán xạ, bê tông thể hiện hiệu quả kinh tế và cản bức xạ tán xạ góc rộng vượt trội, làm suy giảm suất liều tán xạ tại vị trí 1 m xuống dưới mức 2.5 microSievert/giờ với bề dày lớp chắn bê tông khoảng 20 cm.

Thảo luận kết quả

Sự chênh lệch giữa kết quả giải tích đơn giản và mô phỏng thực tế bắt nguồn từ bản chất vật lý của tán xạ Compton trong vùng năng lượng từ 0.5 MeV đến 2.0 MeV. Quá trình này không loại bỏ hoàn toàn photon mà chỉ làm suy giảm năng lượng và chuyển hướng lan truyền của chúng, tạo nên phổ năng lượng thứ cấp liên tục đến đầu dò. Khi so sánh với các nghiên cứu kinh điển của Goldstein và Wilkins cũng như các mô phỏng chuẩn bằng mã MCNP, kết quả từ PENELOPE 2006 trong luận văn này đạt sự tương thích cao với độ sai lệch dao động trong khoảng từ 2.1% đến 4.5%.

Toàn bộ dữ liệu nghiên cứu có thể được trực quan hóa tối ưu thông qua các đồ thị phân bố không gian và bảng so sánh đa biến. Cụ thể, dữ liệu suy giảm suất liều được biểu diễn qua đồ thị đường cong suy giảm phi tuyến theo bề dày từ 0 đến 20 cm của từng vật liệu. Bảng tổng hợp các thông số suy giảm của chì, sắt và bê tông phản ánh rõ ràng mối tương quan giữa khối lượng riêng và hiệu quả che chắn, cung cấp cơ sở dữ liệu định lượng chuẩn xác cho các nhà thiết kế công trình hạt nhân.

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

Thứ nhất, chuẩn hóa quy trình thiết kế che chắn bức xạ bằng việc kết hợp song song phương pháp giải tích sơ bộ và mô phỏng kiểm chứng Monte Carlo. Chủ thể thực hiện là các viện nghiên cứu hạt nhân và phòng an toàn bức xạ y tế, triển khai áp dụng từ năm 2024 nhằm đạt mục tiêu kiểm soát sai số tính toán dưới 3% cho mọi công trình phòng xạ trị.

Thứ hai, tối ưu hóa cấu trúc che chắn đa lớp kết hợp giữa chì và bê tông nhằm nâng cao hiệu quả kinh tế và an toàn bức xạ. Các bệnh viện và trung tâm chiếu xạ công nghiệp cần tiến hành đánh giá lại chiều dày tấm chắn định kỳ 6 tháng một lần, hướng tới mục tiêu duy trì suất liều làm việc xung quanh nguồn dưới 5 microSievert/giờ và tiết kiệm ít nhất 15% khối lượng chì sử dụng.

Thứ ba, thiết lập hệ thống kiểm soát liều chiếu xạ tích lũy cá nhân tự động và quan trắc môi trường bức xạ liên tục. Đơn vị an toàn bức xạ tại các cơ sở cần trang bị liều kế nhiệt phát quang và máy đo suất liều cầm tay, đảm bảo kiểm soát tổng liều hàng năm của nhân viên vận hành dưới 10 mSv/năm, tức thấp hơn 50% so với trần quy định của quốc gia.

Thứ tư, đẩy mạnh nghiên cứu và ứng dụng các vật liệu che chắn composite thế hệ mới như polymer lai hạt nặng hoặc bê tông gia cường hạt nano. Các phòng thí nghiệm vật lý vật liệu cần phối hợp với các cơ sở bức xạ trong giai đoạn 2024 đến 2027 để thử nghiệm khả năng cản xạ, hướng đến giảm từ 30% đến 50% trọng lượng của các thùng chứa nguồn phóng xạ di động.

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

Nhóm thứ nhất bao gồm các kỹ sư an toàn bức xạ và chuyên gia vật lý y khoa tại các bệnh viện ung bướu. Luận văn cung cấp phương pháp luận chuẩn xác để tính toán độ dày cửa chì, tường bảo vệ phòng đặt máy xạ trị Coban-60 và phòng chụp cắt lớp vi tính, giúp đảm bảo an toàn tuyệt đối cho nhân viên y tế và bệnh nhân.

Nhóm thứ hai là các nhà nghiên cứu, giảng viên và sinh viên sau đại học chuyên ngành Vật lý Kỹ thuật, Kỹ thuật Hạt nhân và Kỹ thuật Y sinh. Tài liệu này đóng vai trò như một cẩm nang chuyên sâu về cơ sở tương tác bức xạ, thuật toán tích phân giải tích và quy trình thiết lập bài toán mô phỏng hạt trên phần mềm PENELOPE 2006.

Nhóm thứ ba là các cán bộ quản lý an toàn và kỹ thuật viên tại các cơ sở chiếu xạ công nghiệp, khử trùng y tế và kiểm tra khuyết tật mối hàn bằng chụp ảnh phóng xạ gamma. Use case thực tế là ứng dụng công thức tính toán nguồn dạng dây và nguồn trụ để thiết kế thùng chứa bảo quản nguồn và xác lập vành đai an toàn khi tác nghiệp ngoài hiện trường.

Nhóm thứ tư là các kiến trúc sư và kỹ sư xây dựng công trình hạt nhân tại các đơn vị tư vấn thiết kế. Tài liệu cung cấp các bảng tra cứu bề dày suy giảm của bê tông và kim loại, hỗ trợ trực tiếp cho việc thẩm tra kết cấu che chắn phòng bức xạ đạt chuẩn quy chuẩn kỹ thuật quốc gia.

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

Phương pháp Monte Carlo có ưu điểm vượt trội gì so với phương pháp giải tích trong tính toán che chắn bức xạ gamma? Phương pháp Monte Carlo theo dõi từng tương tác ngẫu nhiên của hàng triệu photon trong vật chất, cho phép mô phỏng chính xác các cấu hình hình học phức tạp như nguồn trụ hoặc nguồn cầu đa lớp. Phương pháp giải tích thường phải đơn giản hóa hình học và khó xử lý chính xác hiện tượng tán xạ thứ cấp nhiều lần, trong khi Monte Carlo đạt độ lệch chuẩn dưới 2% khi theo dõi 10^7 lịch sử hạt.

Hệ số tích lũy (Buildup factor) đóng vai trò quan trọng như thế nào trong bài toán tính liều? Khi chùm tia gamma rộng đi qua lớp vật chất dày, các photon tán xạ Compton liên tiếp vẫn có thể đi tới đầu dò, làm tăng suất liều thực tế. Nếu bỏ qua hệ số tích lũy trong tính toán giải tích, kết quả thu được sẽ đánh giá thấp suất liều thực tế từ 20% đến 60%, dẫn đến nguy cơ thiết kế thiếu bề dày che chắn và gây mất an toàn bức xạ.

Tại sao chì và bê tông là hai vật liệu được lựa chọn phổ biến nhất trong thiết kế che chắn gamma? Chì có khối lượng riêng cao đạt 11.34 g/cm3 và nguyên tử số lớn Z = 82, mang lại hiệu quả hấp thụ quang điện tuyệt vời đối với chùm tia sơ cấp trong không gian hẹp. Bê tông có giá thành rẻ, độ bền cơ học cao và khả năng cản xạ tốt đối với các tia tán xạ năng lượng trung bình, rất phù hợp làm lớp bảo vệ thứ cấp cho tường bao công trình.

Các giới hạn liều bức xạ hiện hành được quy định cụ thể như thế nào cho con người? Theo tiêu chuẩn ICRP 60 và quy chuẩn an toàn bức xạ Việt Nam, liều hiệu dụng tối đa cho nhân viên bức xạ là 20 mSv/năm (tính trung bình trong 5 năm liên tiếp và không vượt quá 50 mSv trong một năm đơn lẻ). Đối với công chúng, mức giới hạn liều hiệu dụng nghiêm ngặt là 1 mSv/năm. Suất liều tại ranh giới khu vực kiểm soát thường được thiết kế dưới 8 microSievert/giờ.

Phần mềm PENELOPE 2006 mô phỏng quá trình truyền hạt photon và electron theo cơ chế nào? PENELOPE 2006 áp dụng thuật toán mô phỏng chi tiết cho photon và cơ chế tán xạ hỗn hợp cho electron và positron trong dải năng lượng từ 50 eV đến 1 GeV. Chương trình sử dụng hệ cơ sở dữ liệu tiết diện tán xạ đáng tin cậy để bám sát từng quá trình ion hóa, kích thích, phát xạ tia X đặc trưng và hiện tượng hủy cặp electron-positron.

Kết luận

• Hệ thống hóa toàn diện cơ sở lý thuyết về tương tác của bức xạ gamma với vật chất và các đại lượng liều lượng học an toàn hạt nhân theo tiêu chuẩn quốc tế.
• Khảo sát và xây dựng thành công các mô hình giải tích tính toán suất liều cho các cấu hình nguồn điểm, nguồn dây, nguồn cầu và nguồn hình trụ đặc.
• Ứng dụng xuất sắc chương trình Monte Carlo PENELOPE 2006 để mô phỏng sự suy giảm bức xạ qua các vật liệu chì, sắt và bê tông với độ chính xác cao.
• Kiểm chứng và xác định rõ ràng sự ảnh hưởng của hệ số tích lũy và hiệu ứng tự hấp thụ, kiểm soát độ sai lệch giữa mô phỏng và giải tích dưới 4.5%.
• Cung cấp bộ dữ liệu thông số bề dày cản xạ thực nghiệm phục vụ đắc lực cho công tác thiết kế an toàn bức xạ y tế và công nghiệp.

Trong giai đoạn 2024 đến 2026, hướng nghiên cứu cần tiếp tục mở rộng mô phỏng trên các cấu hình che chắn vật liệu nano và nguồn bức xạ phức hợp. Các nhà khoa học, kỹ sư an toàn và đơn vị quản lý bức xạ hãy tham khảo ngay luận văn này để ứng dụng các mô hình tính toán chuẩn xác, tối ưu hóa chi phí và nâng cao mức độ an toàn cho các công trình hạt nhân.

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

LỜI MỞ ĐẦU – MỤC TIÊU VÀ NHIỆM VỤ LUẬN VĂN Kể từ khi Becquerel phát hiện ra hiện tƣợng phân rã phóng xạ, con ngƣời chúng ta đã ý thức đƣợc xung quanh chúng ta ngập tràn các loại bức xạ, và giờ đây phóng xạ và bức xạ đã trở thành những khái niệm quen thuộc. Tia gamma và tia X là sóng điện từ có năng lƣợng từ 50 keV trở lên. Vì thế, sức xuyên thấu của nó rất lớn và chỉ có thể chặn lại bằng vật liệu có nguyên tử lƣợng lớn nhƣ chì hoặc bêtông, nƣớc nặng. Khi nói đến bức xạ nói chung và bức xạ hạt nhân nói riêng mọi ngƣời thƣờng nghĩ ngay đến tác hại của nó.

Tác hại của bức xạ hạt nhân đƣợc thể hiện rõ rệt qua hậu quả hai quả bom nguyên tử mà Mỹ thả xuống Nhật Bản trong chiến tranh thế giới thứ II. Và gần đây là các thảm họa nhà máy điện hạt nhân Mayak (Liên xô, 1957) Three Miles Island (Mỹ, 1979), Trecnobƣn (Liên xô, 1986), Fukushima (Nhật, 2011). Khi đi qua vật chất các photon va chạm và tƣơng tác với các nguyên tử của vật chất theo các hiệu ứng hấp thụ và tán xạ. Trong các quá trình đó chúng có thể truyền một phần hoặc hoàn toàn năng lƣợng của mình cho vật chất.

Khi hấp thụ năng lƣợng của tia gamma, các nguyên tử phân tử của vật chất bị biến đổi, trở thành các ion. Chính sự biến đổi này có thể dẫn tới sự phá huỷ cấu trúc của vật chất. Với liều lƣợng khá lớn, bức xạ ion hóa có thể gây ra các tác động có hại cấp tính tới sức khoẻ và có thể dẫn tới tử vong đối với con ngƣời. Bên cạnh đó còn có những tác dụng phụ ảnh hƣởng tới hệ thống thần kinh trung ƣơng, hệ thống tiêu hóa, hệ thống hô hấp.

Phổ biến nhất là những bệnh nhân ung thƣ trải qua xạ trị thƣờng phải chịu rất nhiều tác dụng phụ nhƣ: rụng tóc, đau bụng, giảm lƣợng bạch cầu trong máu, kích ứng da, loét miệng,. Với liều lƣợng nhỏ, bức xạ ion hóa có thể gây ra những tác động lâu dài mang tính ngẫu nhiên ảnh hƣởng đến việc tạo đột biến trên bản thân ngƣời bị chiếu nhƣ gây ung thƣ, bất hoạt cơ quan hoặc ảnh hƣởng di truyền về sau. Tuy nhiên, bức xạ hạt nhân có rất nhiều ứng dụng quan trọng trong nhiều lĩnh vực kỹ thuật công nghệ, nông nghiệp, y học, môi trƣờng vv… Có thể kể đến các ứng dụng của GVHD: TS. Huỳnh Quang Linh Luận văn CH ngành Vật lý Kỹ thuật 6 bức xạ nhƣ : sản xuất vật liệu cách điện, bền nhiệt, polymer, khử trùng các dụng cụ y tế, ứng dụng trong chụp ảnh công nghiệp, sử dụng bức xạ gamma cƣờng độ cao cho các mục đích khử trùng, biến tính vật liệu, bảo quản thực phẩm và nông sản, cải tạo sinh khối, chế tạo một số chế phẩm bằng bức xạ, trong nông nghiệp sử dụng bức xạ gamma kết hợp với những tác nhân khác để cải tạo giống cây trồng, sử dụng đồng vị đánh dấu để nghiên cứu các quá trình sinh học.

Các ứng dụng khác dùng trong các ngành công nghiệp nhƣ xi măng, giấy, thủy tinh, công nghiệp thực phẩm, xây dựng, giao thông vận tải, cơ khí nhƣ đo độ ẩm, mật độ, cân trọng lƣợng và ứng dụng để thăm dò, khai thác dầu khí. Nhiều trƣờng đại học, viện nghiên cứu sử dụng nguồn bức xạ trong công tác nghiên cứu và giảng dạy. Trong ngành hải quan, nguồn bức xạ dùng trong các máy soi hành lý. Do mức độ tăng vọt của các hình thức thƣơng mại hóa hạt nhân nguyên tử nên nỗi lo lắng của ngƣời dân về năng lƣợng hạt nhân vẫn không giảm bớt.

Trong 3 năm qua, thế giới đã 3 lần tổ chức các hội nghị thƣợng đỉnh về an ninh hạt nhân. Hội nghị nào cũng thu hút sự tham dự của đông đảo nguyên thủ các nƣớc cùng lãnh đạo các tổ chức quốc tế lớn, cho thấy tầm quan trọng của vấn đề này. Việt Nam và Ủy ban châu Âu vừa họp bàn về những mong đợi của hai bên trong tƣơng lai ở lĩnh vực an toàn hạt nhân. Dự án hợp tác an toàn hạt nhân giữa châu Âu và Việt Nam có giá trị 2 triệu euro và sẽ đƣợc thực hiện trong 3 năm.

Nhà máy điện hạt nhân đầu tiên của Việt Nam dự kiến đặt tại tỉnh Ninh Thuận, gồm hai nhà máy có tổng công suất 4. Theo kế hoạch, nhà máy điện hạt nhân Ninh Thuận 1 sẽ đƣợc khởi công năm 2014 và bắt đầu phát điện năm 2020. Dự kiến, năm 2015 sẽ triển khai tiếp nhà máy điện Ninh Thuận 2. Việc nghiên cứu và hiểu biết về kỹ thuật che chắn an toàn bức xạ giúp chúng ta tìm ra giải pháp, tiêu chí an toàn nhất khi sử dụng, tiếp xúc với bức xạ trong khi vẫn khai thác đƣợc những lợi ích, tiềm năng đa dạng của nó.

Thực tế, việc tính toán thiết kế che chắn các nguồn bức xạ gamma vẫn thƣờng đƣợc thực hiện bằng phƣơng pháp giải tích cho những cấu hình hình học đơn giản, nhƣng GVHD: TS. Huỳnh Quang Linh Luận văn CH ngành Vật lý Kỹ thuật 7 với những cấu hình phức tạp, phƣơng pháp giải tích có những trở ngại dẫn đến việc ƣớc lƣợng mang tính gần đúng nhiều lúc vƣợt quá yêu cầu cần thiết. Do vậy, việc sử dụng các chƣơng trình Monte Carlo mô phỏng quá trình lan truyền bức xạ trong vật chất có thể khắc phục nhƣợc điểm trên và mặt khác là một lựa chọn nhằm kiểm tra độ tin cậy của các phƣơng pháp tính toán giải tích. Việc ứng dụng chƣơng trình mô phỏng Monte Carlo cho mục đích che chắn an toàn bức xạ là một kỹ năng nghiên cứu cần phát triển để tăng cƣờng công cụ hiệu quả trong kỹ thuật hạt nhân thực nghiệm.

Với tổng quan trên, đề tài “Khảo sát và mô phỏng tính toán các cấu hình che chắn của một số nguồn bức xạ gamma bằng phƣơng pháp Monte Carlo” đƣợc thực hiện với mục tiêu và nhiệm vụ đƣợc đề ra nhƣ sau: Mục tiêu của đề tài: Khảo sát các phƣơng pháp tính toán giải tích kỹ thuật che chắn bức xạ gamma và thực hiện mô phỏng việc tính toán che chắn bằng phƣơng pháp Monte Carlo qua chƣơng trình PENELOPE 2006. Các nhiệm vụ của đề tài: -Khảo sát tổng quan về kỹ thuật che chắn bảo vệ an toàn bức xạ hạt nhân. -Tổng quan về bức xạ gamma, các quá trình tƣơng tác chủ yếu của gamma với vật chất, giới thiệu các đại lƣợng đo lƣờng hạt nhân đặc trƣng và an toàn phóng xạ. -Tổng quan về phƣơng pháp tính toán giải tích về che chắn bức xạ hạt nhân.

-Giới thiệu về chƣơng trình mô phỏng sự lan truyền photon và electron năng lƣợng cao trong vật chất PENELOPE 2006 -Sử dụng chƣơng trình PENELOPE 2006 mô phỏng một số cấu hình che chắn bức xạ gamma. Huỳnh Quang Linh Luận văn CH ngành Vật lý Kỹ thuật 8 CHƢƠNG I TƢƠNG TÁC CỦA BỨC XẠ GAMMA VỚI VẬT CHẤT 1.Bản chất của tia gamma [1] 1. Tia Gamma Bức xạ gamma là hiện tƣợng hạt nhân chuyển từ trạng thái kích thích cao về trạng thái kích thích thấp hơn hay trạng thái cơ bản bằng cách phân rã bức xạ gamma. Hạt nhân sau khi phân rã gamma không thay đổi số khối A và điện tích Z.

Công thức của quá trình là: ( Z X A )* Z X A   (1.1) Nhƣ vậy phân rã Gamma chỉ xảy ra trong hạt nhân kích thích. Hạt nhân kích thích có thể thu đƣợc bằng nhiều cách. Trong phản ứng hạt nhân do kết quả bắn phá bởi các hạt tích điện, hạt nơtron và các photon khác mà hạt nhân bia sẽ đƣợc nâng lên trạng thái kích thích. Sau phân rã anpha và bêta hạt nhân mới tạo thành có thể nằm ở trạng thái kích thích.

Đây là cách đơn giản và phổ biến nhất để thu đƣợc những hạt nhân bức xạ gamma. 14 Tia gamma là các sóng điện từ có bƣớc sóng cực ngắn ( 10 m  10 12 m ) phát ra trong sự phân rã của hạt nhân nguyên tử. Cần chú ý là tia gamma và tia Rơnghen đều là các sóng điện từ nhƣng gamma có bƣớc sóng ngắn hơn, độ đâm xuyên lớn hơn ( có thể xuyên qua lớp vật liệu dày nhƣ xuyên qua 10cm Pb). Điểm khác biệt cơ bản, bản chất của 2 loại tia ở chỗ: Gamma phát ra từ hạt nhân nguyên tử còn tia rơnghen phát ra bởi sự dịch chuyển trạng thái từ trạng thái năng lƣợng cao về trạng thái năng lƣợng thấp hơn của các electron trong lớp vỏ nguyên tử.

Tia rơnghen đƣợc tạo ra khi các điển tử bị hãm, những điện tử có tốc độ cao đến gần ( bay ngang qua) hạt nhân mang điện dƣơng chúng chịu một sức hút và chuyển động chậm lại. Trong quá trình bị làm chậm GVHD: TS. Huỳnh Quang Linh Luận văn CH ngành Vật lý Kỹ thuật 9 hoặc bị hãm chúng bị mất đi một phần động năng ban đầu dƣới dạng tia Rơnghen ( theo điện động lực học cổ điển ). Đôi khi điện tử bị dừng lại tức thì, khi đó toàn bộ năng lƣợng của chúng đƣợc chuyển thành năng lƣợng cực đại của bức xạ Rơnghen ( có bƣớc sóng cực tiểu ).

Trong ống phóng tia Rơnghen, khi các electron đập vào bia thì tốc độ của chúng thay đổi liên tục trong trƣờng Culông của các nguyên tử bia hay nói cách khác năng lƣợng của các electron bị mất dần, do đó các tia rơnghen phát ra có bƣớc sóng thay đổi liên tục trong một dải rộng vì vậy phổ của nó là một phổ liên tục. Các đồng vị phóng xạ phát ra một hoặc nhiều loại sóng. Do vậy trong phổ của nó có cách vạch khác nhau tƣơng ứng với các sóng đó. 60 Ví dụ : Co ( T=5,2714 năm ) phát ra 2 loại sóng trội ở 1173,228 KeV và 1332,490 KeV.

137 Cs ( T=30,07 năm ) phát ra 1 loại sóng trội ở 661,657 KeV. 57 Co ( T= 271,79 ngày ) phát ra 3 loại sóng trội 14,6 KeV, 122,06 KeV, 136,47 KeV. 241 Am ( T=432,2 năm ) phát ra 2 loại sóng trội 26,3446 KeV và 54,5409 KeV. Tất cả các nguồn gamma đều có dạng phổ vạch ( Tức là có năng lƣợng xác định và có các đỉnh đặc trƣng cho năng lƣợng của chúng là rời rạc ) khác với phổ của tia rơnghen.

Thƣờng ngƣời ta mô tả một tia gamma nhất định bằng năng lƣợng của nó hơn là biểu diễn nó theo chiều dài bƣớc sóng. Ví dụ nhƣ phổ của một nguồn gamma có dạng nhƣ sau : GVHD: TS. Huỳnh Quang Linh Luận văn CH ngành Vật lý Kỹ thuật 10 Cƣờng độ Số kênh Hình 1.1: Phổ của một nguồn gamma 1. Tính chất của tia gamma - Không nhìn thấy đƣợc bằng mắt thƣờng.