Tính toán và đo đạc suất liều chiếu tại một số vị trí của máy gia tốc uelr 10 15s2

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM TP HỒ CHÍ MINH Huỳnh Kiều Viết Lãm TÍNH TOÁN VÀ ĐO ĐẠC SUẤT LIỀU CHIẾU TẠI MỘT SỐ VỊ TRÍ CỦA MÁY GIA TỐC UELR-10-15S2 LUẬN VĂN THẠC SĨ KHOA HỌC VẬT CHẤT Thành phố Hồ Chí Minh – 2016 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM TP HỒ CHÍ MINH Huỳnh Kiều Viết Lãm TÍNH TỐN VÀ ĐO ĐẠC SUẤT LIỀU CHIẾU TẠI MỘT SỐ VỊ TRÍ CỦA MÁY GIA TỐC UELR-10-15S2 Chuyên ngành: Vật lí nguyên tử Mã số: 60 44 01 06 LUẬN VĂN THẠC SĨ KHOA HỌC VẬT CHẤT NGƯỜI HƯỚNG DẪN KHOA HỌC: TS TRẦN VĂN HÙNG Thành phố Hồ Chí Minh – 2016 LỜI CẢM ƠN Luận văn thực Trung tâm nghiên cứu triển khai công nghệ xạ (VINAGAMMA) số 202A đường 11, Phường Linh Xuân, Quận Thủ Đức Trong trình thực luận văn, nhận quan tâm, hướng dẫn, giúp đỡ động viên thầy cơ, đồng nghiệp, gia đình bạn Nay tơi chân thành bày tỏ lịng tri ân sâu sắc đến: TS.Trần Văn Hùng, trưởng phòng kế hoạch tài chính, nghiên cứu viên Trung tâm nghiên cứu triển khai công nghệ xạ - Viện Năng Lượng Nguyên tử Việt Nam Thầy định hình cho tơi lựa chọn đề tài đồng thời tận tình hướng dẫn tơi suốt thời gian nghiên cứu hồn thành tốt luận văn Đặc biệt, học Thầy phương pháp làm việc khoa học, chuyên nghiệp với trình độ chun mơn cao học thực tiễn, từ vốn sống, từ am hiểu thấu đáo riêng Thầy Quý thầy cô giáo khoa Vật Lý trường Đại học Sư phạm TP.Hồ Chí Minh, người đem lại cho tơi kiến thức bổ trợ, vơ có ích năm học vừa qua, giúp trưởng thành học tập nghiên cứu khoa học Ban giám hiệu, Phòng Đào tạo sau đại học trường Đại học Sư phạm TP.Hồ Chí Minh, Ban Chủ nhiệm Khoa Vật lí, tạo điều kiện tốt cho tơi q trình học tập, nghiên cứu Q thầy hội đồng bảo vệ luận văn Thạc sĩ đọc, đóng góp ý kiến, nhận xét đánh giá luận văn Sở Giáo dục-Đào tạo Thành phố Hồ Chí Minh, Ban giám hiệu Trường THPT Ernst Thalmann đồng nghiệp tạo nhiều điều kiện thuận lợi giúp đỡ tơi q trình học tập, thực luận văn Cuối xin khắc sâu công ơn Cha mẹ, anh em, người thân, bạn bè đặc biệt bạn lớp Vật Lý Nguyên tử K25, K24, người bên tôi, động viên khuyến khích tơi q trình thực đề tài nghiên cứu LỜI CAM ĐOAN Tơi tên là: Huỳnh Kiều Viết Lãm, học viên cao học lớp Vật lí Khóa 25, chun ngành Vật lí ngun tử, khố 2014-2016 Tơi xin cam đoan luận văn thạc sĩ ‘‘Tính toán đo đạc suất liều chiếu số vị trí máy gia tốc UELR-10-15S2’’là cơng trình nghiên cứu riêng tôi, số liệu nghiên cứu thu từ tính tốn, mơ thực nghiệm, khơng chép, trung thực chưa sử dụng để bảo vệ học vị Mọi giúp đỡ cho việc thực luận văn cảm ơn thơng tin trích dẫn luận văn rõ nguồn gốc, phép công bố Học viên Huỳnh Kiều Viết Lãm MỤC LỤC Lời cảm ơn Lời cam đoan Danh mục bảng Danh mục hình vẽ Danh mục chữ viết tắt MỞ ĐẦU Chương HIỆU ỨNG PHÁT BỨC XẠ VÀ CÁC KHÁI NIỆM VỀ LIỀU GIỚI HẠN HẰNG NĂM CHO CÁC ĐỐI TƯỢNG 1.1 Hiệu ứng phát xạ hãm 1.1.1 Tương tác electron với vật chất 1.1.2 Phổ lượng photon sinh chùm tia điện tử 1.2 Các giới hạn liều cho phép 1.2.1 Suất liều chiếu 1.2.2 Suất liều hấp thụ 1.2.3 Liều hiệu dụng 11 1.2.4 Các giới hạn an toàn cho nhân viên xạ dân cư 14 Chương GIỚI THIỆU VỀ MÁY GIA TỐC UELR-10-15S2 VÀ CHƯƠNG TRÌNH MCNP 16 2.1 Hệ thống chùm tia máy gia tốc UELR-10-15S2 17 2.2 Cấu trúc bê tông che chắn 20 2.3 Giới thiệu chương trình tính tốn MCNP 21 2.3.1 Giới thiệu chung 21 2.3.2 Dạng INPUT file MCNP 22 Chương TÍNH TỐN VÀ ĐO ĐẠC SUẤT LIỀU 23 3.1 Tính tốn an tồn 23 3.1.1 Tính tốn truyền qua xạ hãm 23 3.1.2 Tính tốn tán xạ chùm xạ hãm 25 3.2 Mơ tính tốn MCNP 28 3.2.1 Mô cấu trúc buồng chiếu 28 3.2.2 Tính tốn suất liều 29 3.3 Đo suất liều khu vực kiểm soát 32 3.3.2 Đo suất liều bên buồng chiếu 34 Chương KẾT QUẢ TÍNH TỐN, ĐO ĐẠC VÀ ĐÁNH G IÁ AN TOÀN 37 4.1 Các kết tính tốn đo đạc 37 4.1.1 Bản đồ liều tầng 37 4.1.2 Bản đồ liều tầng hai 39 4.1.3 Bản đồ liều buồng chiếu 40 4.2 Đánh giá an toàn xạ thiết bị vận hành 43 4.2.1 Tiêu chuẩn đánh giá an toàn xạ 43 4.2.2 Đánh giá an toàn xạ thiết bị vận hành 44 4.2.3 Khu vực cảnh báo xạ 45 KẾT LUẬN 46 TÀI LIỆU THAM KHẢO 48 PHỤ LỤC: File INPUT chương trình MCNP mơ tả tồn hệ máy gia tốc UELR-10-15S2 50 DANH MỤC CÁC BẢNG Bảng 1.1 Hệ số trọng số phóng xạ vài loại xạ (ICRP – 1990) 11 Bảng 1.2 Các trọng số mô đặc trưng cho mô thể W T (1990) 13 Bảng 2.1 Các thông số thiết bị UERL-10-15S2 19 Bảng 3.1 Thời gian đặt liều kế Fricke cho điểm buồng chiếu 34 Bảng 4.1 Kết tính tốn đo đạc tầng 38 Bảng 4.2 Kết tính tốn đo đạc suất liều tầng hai 40 Bảng 4.3 Kết tính đo suất liều buồng chiếu 41 Bảng 4.4 So sánh số liệu với báo 48 Bảng 4.5 So sánh suất liều thực tế với tiêu chuẩn ICRP 44 Bảng 4.6 Suất liều thời gian làm việc tối đa buồng chiếu 45 DANH MỤC CÁC HÌNH VẼ Hình 1.1 Sự xạ Synchrotron Hình 1.2 Phổ Photon từ bia Ta Hình 2.1 Đặc điểm, vị trí hai chùm tia 17 Hình 2.2 Phân bố cường độ chùm xạ hãm theo góc bay 18 Hình 2.3 Cấu trúc bê tơng che chắn phóng xạ máy gia tốc 20 Hình 3.1.Suất liều xạ hãm bia nặng theo hướng tia tới vng góc 24 Hình 3.2 Sự phụ thuộc TVL vào lượng vật liệu che chắn (a) bê tơng, (b) thép, (c) chì 25 Hình 3.3 Giản đồ chùm tia tán xạ nhiều lần 26 Hình 3.4 Hệ số tán xạ lên bê tông theo lượng 27 Hình 3.5 Mặt cắt vng góc trục Oy gốc tọa độ 29 Hình 3.6 Mặt cắt vng góc trục Ox gốc tọa độ 29 Hình 3.7.Tối ưu hóa khơng gian để tính liều phịng điều khiển 31 Hình 3.8 Mơ hình tính liều tán xạ điểm số 32 Hình 3.9 Phân bố điểm đặt liều kế Fricke 34 Hình 3.10 Phân bố điểm đặt liều kế TLD tầng 36 Hình 4.1 Phân bố điểm tính đo suất liều tầng 37 Hình 4.2 Phân bố điểm tính đo suất liều tầng hai 39 Hình 4.3 Phân bố điểm tính đo suất liều buồng chiếu 41 DANH MỤC CÁC TỪ VIẾT TẮT Từ Tiếng Anh Tiếng Việt TVL Tenth-Value-Layer Bề dày suy giảm mười lần LINAC Linear Accelerator Máy gia tốc tuyến tính Imp Important point Độ quan trọng MCNP Monte – Carlo N- particle Chương trình mơ Monte-carlo TLD Thermoluminescent Dosimeter Liều kế nhiệt phát quang IAEA International Atomic Enegy Agency Cơ quan Năng Lượng Nguyên Tử Quốc Tế ICRP International Commission on Ủy ban Quốc Tế Về An toàn Bức xạ Radiological Protection RF Radio Frequency Sóng vơ tuyến MỞ ĐẦU Trong năm gần đây, công nghệ xạ phát triển ngày nhanh chóng giới Việt Nam lĩnh vực Công nghiệp, Nông nghiệp, y tế mơi trường Vấn đề an tồn xạ quan tâm hàng đầu sở chiếu xạ có sử dụng đồng vị phóng xạ hay máy gia tốc Cùng với quy tắc nghiêm ngặt an ninh nguồn phóng xạ, an tồn xạ đánh giá nhiều phương diện khác nhau, nhằm đảm bảo tốt sức khỏe cho nhân viên vận hành, khách vãng lai dân cư xung quanh sở Đối với sở vận hành máy gia tốc bao gồm sử dụng biến đổi sang tia X dùng trực tiếp chùm điện tử, công việc đánh giá an toàn xạ đặt lên hàng đầu Năm 2013, Trung tâm Nghiên cứu Triển khai Công nghệ Bức xạ lắp đặt vào vận hành máy gia tốc chùm tia điện tử lượng trung bình 10 MeV Một nhiệm vụ cấp thiết đặt cần tính tốn đo đạc suất liều khu vực kiểm soát phương pháp mới, phù hợp khả thi nhằm cho số liệu đáng tin cậy, từ đánh giá chuẩn xác mức độ an toàn xạ máy gia tốc, để tiến hành thiết lập qui trình ứng phó trường hợp thiết bị có cố nhân viên bị chiếu liều Cũng năm 2013 nhóm cán trung tâm tính tốn đo đạc suất liều khu vực kiểm soát thiết bị [9] Tuy nhiên kết tính tốn kết đo đạc có chênh lệch lớn nên cần phải tiến hành tính tốn đo đạc lại Từ lý tơi chọn đề tài: ‘‘Tính tốn đo đạc suất liều chiếu số điểm máy gia tốc UELR-10-15S2’’, để đánh giá mức độ an toàn thiết bị cách xác hơn, đồng thời dựa số liệu đáng tin cậy mà thiết lập trung tâm sử dụng xây dựng kế hoạch ứng phó cố để cục an tồn phê duyệt Dựa theo cấu trúc bê tơng che chắn phóng xạ thiết kế xây dựng, phương pháp xấp xỉ triển khai áp dụng để tính tốn suất liều vị trí quan trọng vùng kiểm soát máy gia tốc UELR-10-15S2 Phương pháp xấp xỉ triển 42 Bộ số liệu phạm vi đề tài so sánh với số liệu mà nhóm cán trung tâm tính tốn đo đạc lúc vận hành chạy thử máy (khi mà máy hoạt động chưa thật ổn định, thời gian chạy máy giãn cách) bảng 4.4 Bảng 4.4 So sánh số liệu với báo Điểm Đề tài Bài báo [9] MCNP Đo thực nghiệm MCNP Đo thực nghiệm 4,2E+08 1,3E+08 (1,80±0,06).10 (1,50±0,07).10 6,9E+06 2,2E+06 (3,90±0,15).10 (3,20±0,16).10 3,2E+04 3,2E+04 (3,30±0,09).10 (3,40±0,15).10 4 1,4E-03 4,9E-01 (1,50±0,25).10 -3 (6,2±0,43).10 -1 4,4E+06 9,0E+05 (1,20±0,02).10 (0,99±0,03).10 6 1,3E-03 4,5E-01 (1,40±0,22).10 -3 (4,8±0,32).10 -1 2,5E-04 3,5E-01 (2,30±0,63).10 -4 (3,90±0,70).10 -2 4,6E+07 1,1E+07 (4,20±0,13).10 (4,10±0,12).10 6,4E+06 1,6E+06 (4,40±0,18).10 (4,6±0,14).10 10 9,4E+03 2,8E+4 (2,40±0,36).10 (2,5±0,15).10 11 2,6E-03 3,5E-01 (3,10±0,45).10-3 (3,3±0,25).10 -1 12 5,3E+00 1,5E+00 1,9±0,12 1,8±0,11 13 6,8E-01 6,2E-01 (6,5±0,59).10 -1 (6,4±0,44).10 -1 14 8,3E+01 1,8E+01 (3,31±0,33).10 (2,80±0,20).10 15 4,2E+02 2,47E+02 (4,00±0,24).10 (3,50±0,25).10 Ta thấy báo [9] có chênh lệch lớn mơ MCNP đo đạc thực nghiệm, đặc biệt điểm suất liều cao sai khác đến hai trăm đến ba trăm phần 43 trăm, điều gây khó khăn tính tốn an tồn lên kế hoạch khắc phục cố xạ, đánh giá mức độ an tồn cho cơng nhân viên Nhằm khắc phục điều đề tài sử dụng phương pháp làm giảm thời gian tính máy tăng độ xác như: Tăng độ quan trọng (Imp), tăng bán kính Detecto, tối ưu hóa khơng gian vật chất dựa phán đốn xác hướng hướng phản xạ chùm xạ hãm Nhờ tơi thu số liệu có độ xác cao hơn, đáng tin cậy xây dựng để tiến hành thiết lập quy trình, kế hoạch ứng phó cố xạ nhân viên bị chiếu liều 4.2 Đánh giá an toàn xạ thiết bị vận hành 4.2.1 Tiêu chuẩn đánh giá an toàn xạ Cơ sở đánh giá an toàn xạ dựa tiêu chuẩn ICRP-103 [14], theo liều lượng phép giới hạn tính cho chiếu xạ tồn thân phụ thuộc vào nhóm đối tượng Trong có hai nhóm đối tượng nhân viên xạ (những người làm việc thường xuyên với xạ) chịu mức liều tối đa khơng vượt q 20 mSv/năm (6,0 µSv/h), cịn lại dân chúng chịu mức liều không mSv/năm (0,12 µSv/h) Mức liều giới hạn cịn có đặc điểm cần lưu ý sau: - Liều giới hạn không bao gồm liều nhận từ dịch vụ y tế, - Liều giới hạn tính cho tổng liều hiệu dụng nhận từ chiếu xạ chiếu xạ trong thời gian - Với nhân viên xạ, liều hiệu dụng tính trung bình năm không vượt 50 mSv năm - Với dân chúng, trường hợp đặc biệt liều hiệu dụng cao mức cho phép với điều kiện trung bình năm khơng vượt q mSv/năm Ngồi ra, ICRP khuyến cáo mức liều cho trường hợp đặc biệt để khắc phục cố liên quan đến chiếu xạ Tuy nhiên, với máy gia tốc chùm tia điện tử tất tình xảy cố máy tự động tắt chùm tia nhân viên xạ khắc phục cố thiết bị phát chùm tia 44 4.2.2 Đánh giá an toàn xạ thiết bị vận hành Để đánh giá an toàn xạ cho nhóm đối tượng việc bên ngồi thiết bị chiếu xạ, cần so sánh đối chiếu suất liều nơi làm việc thời gian làm việc nhóm với liều lượng phép giới hạn theo tiêu chuẩn ICRP Trong bảng 4.4, suất liều thực đo tích lũy liều kế TLD kết đo bao gồm phơng phóng xạ mơi trường Đối với nhân viên xạ làm việc khu vực xác định bảng 4.4, suất liều thực tế (bằng phông, ngoại trừ điểm 12) thấp bậc so với tiêu chuẩn khuyến cáo ICRP, nên khu vực làm việc an toàn cho nhân viên xạ Tại khu vực bên buồng chiếu tính từ cửa vào băng tải (điểm 7) tới khu vực nạp dỡ hàng, suất liều đo TLD phơng mơi trường cịn kết từ MCNP khơng tính phơng cỡ 2,5x10 -4 (µSv/h) thấp bốn bậc so với phông Trong bảng 4.5 so sánh với suất liều thực tế với tiêu chuẩn ICRP cho thấy khu vực làm việc nhân viên xạ (các điểm 4, 6, 11, 12, 13), nhân viên tạp vụ khách vãng lai (điểm số 7) ngưỡng an toàn Bảng 4.5 So sánh suất liều thực tế với tiêu chuẩn ICRP Nhóm đối tượng Nhân viên xạ Nhân viên tạp vụ Khách vãng lai Điểm Suất liều thực, Tiêu chuẩn ICRP, (µSv/h) (µSv/h) (6,20±0,43).10 -1 6,0 (4,80±0,32).10 -1 6,0 11 (3,30±0,25).10-1 6,0 12 1,80±0,11 6,0 13 (6,40±0,44).10 -1 6,0 (3,90±0,7).10 -2 0,12 45 4.2.3 Khu vực cảnh báo xạ Ngoài khu vực làm việc thường xuyên nhân viên xạ, nhân viên nạp – dỡ hàng, khách vãng lai xác định an tồn cịn số khu vực suất liều cao cần có cảnh báo xạ Bảng 4.6 Suất liều thời gian làm việc tối đa buồng chiếu Điểm Suất liều (µSv/h) Thời gian làm việc tối đa (h) 14 (2,80±0,20).10 1,3 15 (3,50±0,25).10 0,2 Khu vực xác định suất liều cao buồng chiếu, đặc biệt phát chùm tia từ hướng thẳng lên Bảng 4.5 điểm có suất liều cao thời gian làm việc tối đa cần khắc phục cố Theo bảng 4.5, điểm 14 nhân viên xạ phép làm việc khoảng thời gian 1,3 giờ, điểm 15 phép làm việc 0,2 Tuy nhiên, với đặc trưng máy gia tốc cho phép tắt – mở chùm tia theo ý muốn nên trường hợp có cố độc lập với hệ điều khiển, người vận hành chủ động tắt chùm tia để khắc phục cố Trong thời gian chiếu xạ, cần có biển cảnh báo xạ cầu thang lên buồng chiếu nghiêm cấm việc leo lên buồng chiếu thiết bị vận hành 46 KẾT LUẬN Những kết thu được: Trong luận văn này, tính tốn xấp xỉ, thiết kế mơ tính tốn chương trình MCNP đo đạc thực nghiệm suất liều chiếu khu vực kiểm soát máy gia tốc chùm tia điện tử lượng 10 MeV, làm việc nghiêm túc thu kết sau: + Nắm bắt chế tương tác electron với vật chất, khái niệm liều giới hạn an toàn suất liều cho nhân viên vận hành xạ dân chúng + Hiểu cấu tạo, ngun tắc hoạt động, thơng số dịng máy gia tốc UELR-10-15S2 + Xác định suất liều khu vực kiểm soát máy gia tốc UELR-10-15S2 bao gồm: phòng điều khiển, phòng modulator, cửa vào buồng chiếu, khu nạp – dỡ hàng, đề tài thực phương pháp tính tốn xấp xỉ, mô MCNP đo đạc thực nghiệm, đưa bảng số liệu đáng tin cậy Mặc dù kết thu từ phương pháp tính, mơ phỏng, đo đạc khác nhau, thể chất phương pháp chúng mang ý nghĩa định: - Với phương pháp dùng công thức xấp xỉ, nguồn giả sử nguồn photon 10 MeV đẳng hướng gây suất liều bên buồng chiếu xạ hai đường truyền qua lớp bê tông che chắn tán xạ nhiều lần theo đường băng tải Nguồn photon lượng 10 MeV sử dụng để tính tốn đồng nghĩa với việc xem toàn lượng chùm electron chuyển thành photon chúng bị hãm vật liệu Kết thu thường cao suất liều thực tế sử dụng nhiều tính tốn, thiết kết che chắn theo tiêu chuẩn an tồn phải tính cho trường hợp cực đoan - Phương pháp mơ với cơng cụ chương trình MCNP mơ tả lại tồn cấu trúc bê tơng che chắn, vị trí phát chùm tia electron vật liệu hấp thụ electron sinh xạ hãm Toàn trình tương tác trình vật lý khác 47 MCNP mô tả dựa theo lý thuyết xác suất, từ chương trình cung cấp tiện ích cho phép tính suất liều bên buồng chiếu Tuy nhiên, với suất liều thấp bên buồng chiếu sai số điều khơng thể tránh khỏi, tính tốn kĩ thuật giảm sai số nghiên cứu, áp dụng triệt để Kết thu từ MCNP không bao gồm phông phóng xạ mơi trường, điểm bên ngồi buồng chiếu có suất liều thấp từ ba đến bốn bậc so với phông với sai số 20 % - Để có sở vững cho việc đánh giá an tồn xạ, phương pháp đo liều tích lũy liều kế TLD điểm bên buồng chiếu sử dụng Kết cho thấy suất liều khu vực làm việc nhân viên xạ, khu nạp – dỡ hàng ngang mức phơng mơi trường Điều chứng tỏ đóng góp xạ thiết bị vận hành khơng đáng kể hồn tồn nằm giới hạn cho phép ICRP Kiến nghị: Đề tài tính tốn, đo đạc suất liều buồng chiếu nơi xạ theo đường thơng gió Đây khu vực có suất liều cao, cần cảnh báo khuyến cáo thời gian làm việc khu vực đặc biệt nguy hiểm này, lên thiết kế che chắn an toàn 48 TÀI LIỆU THAM KHẢO Tiếng Việt Phan Quốc Hùng (2002), Phịng tránh phóng xạ an tồn hạt nhân, Nxb Đại học Quốc gia Hà Nội Ngô Quang Huy (2004), An tồn Bức xạ ion hóa, Nxb Đại học Quốc gia Hà Nội Trương Thị Hồng Loan (2008), Phương pháp Monte-Carlo, chuyên đề luận án, Trường Đại học Khoa học Tự nhiên TP.Hồ Chí Minh Trần Đại Nghiệp (2006), Giáo trình xử lý xạ sở công nghệ xạ, Nxb Đại học Quốc gia Hà Nội Trần Đức Nghiệp (2002), Máy gia tốc, Nxb Khoa học Kỹ thuật Mai Văn Nhơn (2001), Vật lý hạt nhân đại cương, Nxb Đại học Quốc gia TP.Hồ Chí Minh Châu Văn Tạo (2004), An tồn xạ ion hóa, Nxb Đại học Quốc gia TP.Hồ Chí Minh Châu Văn Tạo (2006), Liều lượng xạ ion hóa, Nxb Đại học Quốc gia TP.Hồ Chí Minh Nguyễn Anh Tuấn, Tính tốn đo đạc suất liều chiếu khu vực kiểm soát máy gia tốc chùm tia điện tử UELR-10-15S2 trung tâm nghiên cứu triển khai công nghệ xạ, báo cáo toàn văn kỷ yếu hội nghị khoa học lần thứ IX, Trường Đại học Khoa học Tự nhiên TP.Hồ Chí Minh Tiếng Anh 10 ASTM (2004), Standards practice for using the Fricke Reference-Standard dosimetry system, Standards on dosimetry for radiation processing, ASTM (2 nd edition), pp 261-268 11 Corad service (2009), Preliminatary Calculation of Radiation Shielding for Electron Beam System Deliverred Under Contract No 01/12-08-2, St Petertburg, Russia, 2009 49 12 IAEA (1979), Radiological Safety Aspects of the Operation of Electron Linear Accelerators, Technical Reports Series 188 13 ICRP, Recommendations of the International Commission on Radiological Protection, Annals of the ICRP Publication 103, 2007 14 Jean-Luc Biarrotte (2009), RF Cavities for Particle Acceleration, CNRS/IPN Orsay 15 Los Alamos National Laboratory (2000), Monte Carlo N-Partical Code System, Los Alamos, New Mexico, 2000 16 NCRP, Radiation Protection Design Guidelines for 0.1-100 MeV Particle Accelerator Facilities, National Council on Radiation and Measurements, 1997 17 Philip (1990) M.K.Leung, The Physical Basis of Radiotherapy, The Ontario Cancer Institute 18 www.jascoint.co.jp 50 PHỤ LỤC: File INPUT chương trình MCNP mơ tả tồn hệ máy gia tốc UELR-10-15S2 message: datapath=\user\trash testprob03 - Input tinh lieu may gia toc eb 1 -2.3 13 -18 -2 26 -27 imp:p=0 imp:e=0 -2.3 13 -18 -3 26 -27 imp:p=1 imp:e=1 -2.3 13 -14 -11 26 -27 imp:p=0 imp:e=0 -2.3 14 -15 -11 26 -27 imp:p=1 imp:e=1 -2.3 13 -18 11 -12 26 -27 imp:p=1 imp:e=1 -2.3 16 -17 -11 26 -27 29 30 imp:p=1 imp:e=1 -2.3 18 -19 -12 26 -27 imp:p=1 imp:e=1 -2.3 18 -19 -8 26 -27 imp:p=1 imp:e=1 -2.3 19 -21 -8 26 -27 imp:p=0 imp:e=0 10 -2.3 20 -21 -6 26 -27 imp:p=0 imp:e=0 11 -2.3 19 -23 -10 26 -27 imp:p=0 imp:e=0 12 -2.3 22 -23 -9 26 -27 imp:p=0 imp:e=0 13 -2.3 22 -23 -7 26 -27 imp:p=0 imp:e=0 14 -2.3 13 -24 -12 25 -26 imp:p=1 imp:e=1 $ nen dat 15 -2.3 13 -24 -12 27 -28 #39 #40 #52 imp:p=2 imp:e=2 $ tran nha tang 16 -0.00129 15 -18 -11 26 -27 #6 #23 #24 imp:p=1 imp:e=1 $ khong tang 17 -0.00129 18 -22 -9 26 -27 imp:p=0 imp:e=0 18 -0.00129 21 -22 -8 26 -27 imp:p=0 imp:e=0 19 -0.00129 19 -20 -6 26 -27 imp:p=0 imp:e=0 20 -0.00129 23 -24 -10 26 -27 imp:p=0 imp:e=0 51 21 -0.00129 19 -24 10 -12 26 -27 imp:p=0 imp:e=0 23 -0.4 -33 trcl=( 0 92) imp:p=1 imp:e=1 $hang chieu xa 24 -7.6 -34 trcl=(0 81.75) imp:p=1 imp:e=1 $conveyor 25 -2.3 13 -40 -43 28 -53 imp:p=0 imp:e=0 26 -2.3 13 -35 43 -44 28 -53 imp:p=0 imp:e=0 27 0.00129 35 -40 43 -44 28 -53 #44 #45 #46 & #47 #48 #49 #50 #51 imp:p=32 imp:e=32 28 -2.3 13 -40 44 -45 28 -53 #44 #45 #46 & #47 #48 #49 #50 #51 imp:p=16 imp:e=16 29 -2.3 13 -37 45 -46 28 -53 imp:p=2 imp:e=2 30 -2.3 13 -36 46 -12 28 -53 imp:p=2 imp:e=2 31 -2.3 36 -19 50 -12 28 -53 imp:p=0 imp:e=0 32 -2.3 36 51 49 -50 28 -53 imp:p=2 imp:e=2 33 -2.3 38 -39 47 -48 -52 28 -53 imp:p=2 imp:e=2 34 -2.3 39 -40 45 -48 28 -53 imp:p=2 imp:e=2 35 -0.00129 36 -40 45 -50 28 -53 #29 #32 #33 #34 #44 & #45 #46 #47 #48 #49 #50 #51 imp:p=8 imp:e=8 $ khong tang 36 -0.00129 40 -19 -50 28 -53 #37 imp:p=1 imp:e=1 37 -2.3 41 -19 44 -50 28 -53 imp:p=0 imp:e=0 38 -2.3 13 -19 -12 53 -54 #41 #42 #43 & imp:p=4 imp:e=2 $tran nha tang 39 -0.00129 65 -66 -67 -68 27 -61 imp:p=2 imp:e=2 $ chi ngan giua tang va 40 -0.00129 57 -58 59 -60 61 -28 69 imp:p=2 imp:e=2 $ khong giua tang va 41 -0.00129 13 -19 11 -12 56 -54 imp:p=1 imp:e=1 42 -0.00129 13 -19 -44 55 -54 imp:p=1 imp:e=1 52 43 -0.00129 40 -19 44 -48 56 -54 imp:p=1 imp:e=1 c kenh day cap nghieng 60 44 -1 -62 -63 64 trcl=(-132 -255 520 0.866 -0.5 0.5 0.866 0 1) & imp:p=16 imp:e=16 45 like 44 but trcl=(-152 -255 520 0.866 -0.5 0.5 0.866 0 1) 46 like 44 but trcl=(-172 -255 520 0.866 -0.5 0.5 0.866 0 1) 47 like 44 but trcl=(-192 -255 520 0.866 -0.5 0.5 0.866 0 1) 48 like 44 but trcl=(132 -255 520 0.866 0.5 -0.5 0.866 0 1) 49 like 44 but trcl=(152 -255 520 0.866 0.5 -0.5 0.866 0 1) 50 like 44 but trcl=(172 -255 520 0.866 0.5 -0.5 0.866 0 1) 51 like 44 but trcl=(192 -255 520 0.866 0.5 -0.5 0.866 0 1) c ong gia toc 52 -7.6 -69 61 -28 imp:p=2 imp:e=2 100 (-13:24:-1:12:-25:28) (-13:19:-1:12:-28:54) imp:p=0 imp:e=0 py -640 py -490 py -380 py -300 py py 170 py 70 $bo py 220 py 410 10 py 470 53 11 py 630 12 py 910 13 px -420 14 px -230 15 px -130 16 px 60 17 px 210 18 px 400 19 px 540 20 px 910 21 px 960 22 px 1150 23 px 1200 24 px 1320 $bo 25 pz -100 26 pz 27 pz 200 28 pz 280 29 box -75 -85 -52.5 150 0 170 0 105 30 box -75 -85 -52.5 150 0 170 0 105 31 px 85 32 px 85.1 33 box -20 -20 -10 40 0 40 0 20 34 box -25 -30 -0.15 50 0 60 0 0.3 35 px -385 54 36 px -230 37 px -160 38 px -70 39 px 160 40 px 330 41 px 490 43 py -625 44 py -370 45 py -140 46 py 200 47 py 410 48 py 630 49 py 670 50 py 790 51 p -110 670 220 60 790 220 60 790 300 52 p -70 460 220 160 630 220 160 630 300 53 pz 570 54 pz 790 55 pz 585 56 pz 670 57 px -32.5 58 px 32.5 59 py -32.5 60 py 32.5 61 pz 230 55 62 cy 5.5 63 py 135 64 py -135 65 py -7.5 66 py 7.5 67 p -44 -7.5 200 -44 7.5 200 -32.5 7.5 230 68 p 44 -7.5 200 44 7.5 200 32.5 7.5 230 69 cz 20 mode p e m1 1000 0.0847 8000 0.604 11000 0.0095 & 12000 0.0023 13000 0.0248 14000 0.24 & 19000 0.0068 20000 0.0205 26000 0.0046 $betong m2 6000 1.25e-4 7000 6.869e-1 8000 3.01248e-1 & 18000 1.17170e-2 $khong m3 6000 0.8 7000 1.2e-1 8000 0.04 1000 0.04 m4 26000 1.0 m5 82000 1.0 m6 1000 -0.1119 8000 -0.8881$ Nuoc sdef erg=10 pos 0 79 vec 0 dir par f5:p 0 500 30 fm5 9.4e20 f15:p -230 -400 520 30 fm15 9.4e20 f25:p -400 520 30 56 fm25 9.4e20 f35:p 230 -400 520 30 fm35 9.4e20 c f45:p 590 10 90 10 c fm45 9.4e20 de0 0.01 0.03 0.05 0.07 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 & 0.55 0.6 0.65 0.7 0.8 1.0 1.4 1.8 2.2 2.6 2.8 3.25 3.75 4.25 & 4.75 5.0 5.25 5.75 6.25 6.75 7.5 11 c photon flux to dose rate conversion factor c determined by (rem/hr)/(*photon/cm^2s)1 df0 3.96e-6 5.82e-7 2.9e-7 2.58e-7 2.83e-7 3.79e-7 5.01e-7 6.31e-7 & 7.59e-7 8.78e-7 9.85e-7 1.08e-6 1.17e-6 1.27e-6 1.36e-6 1.44e-6 & 1.52e-6 1.68e-6 1.98e-6 2.51e-6 2.99e-6 3.42e-6 3.82e-6 4.02e-6 & 4.41e-6 4.83e-6 5.23e-6 6.60e-6 5.80e-6 6.01e-6 6.37e-6 6.74e-6 & 7.11e-6 7.66e-6 8.77e-6 1.03e-5 tr1 135 -216 95 tr2 135 534 95 ctme 600 ... DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM TP HỒ CHÍ MINH Huỳnh Kiều Viết Lãm TÍNH TOÁN VÀ ĐO ĐẠC SUẤT LIỀU CHIẾU TẠI MỘT SỐ VỊ TRÍ CỦA MÁY GIA TỐC UELR- 10- 15S2 Chuyên ngành: Vật lí nguyên tử Mã số: ... Hệ thống chùm tia máy gia tốc UELR- 10- 15S2 Máy gia tốc UELR- 10- 15S2 thuộc hệ máy gia tốc tuyến tính (LINAC) [14] sử dụng sóng RF để gia tốc chùm electron đến lượng 10 MeV, công suất cực đại 15... gian chạy máy Từ kết đo liều tích lũy thời gian chạy máy tính suất liều (µSv/h) vị trí tương ứng 3.3.2 Đo suất liều bên ngồi buồng chiếu Bên buồng chiếu xạ máy gia tốc UELR- 10- 15S2 bao gồm: phòng