Đồ án kỹ thuật hạt nhânĐồ án kỹ thuật hạt nhânĐồ án kỹ thuật hạt nhânĐồ án kỹ thuật hạt nhânĐồ án kỹ thuật hạt nhânĐồ án kỹ thuật hạt nhânĐồ án kỹ thuật hạt nhânĐồ án kỹ thuật hạt nhânĐồ án kỹ thuật hạt nhânĐồ án kỹ thuật hạt nhânĐồ án kỹ thuật hạt nhân
LỜI CẢM ƠN Lời đầu tiên, để hoàn thành đồ án này, em xin gửi lời cảm ơn sâu sắc tới thầy KS Trần Văn Thống anh chị kỹ sư vật lý công tác Trung tâm y học hạt nhân ung bướu Bệnh viện Bạch Mai quan tâm giúp đỡ em suốt trình thực tập làm đồ án tốt nghiệp Em xin gửi lời cảm ơn chân thành tới thầy Mai Đình Thủy thầy Bùi Ngọc Hà ln giúp em giải đáp thắc mắc kịp thời, động viên em suốt trình làm đồ án Em gửi lời cảm ơn tới cô Trần Thùy Dương hướng dẫn chi tiết cho em phương pháp Monte Carlo để em ứng dụng mơ phân bố liều thiết bị gamma ART-6000 Do lượng kiến thức em nhiều hạn chế nên khơng tránh khỏi thiếu sót, em mong nhận ý kiến đóng góp q báu từ phía thầy bạn để em hồn thiện Cuối cùng, em xin chân thành cảm ơn gia đình, bạn bè- người ln kịp thời động viên giúp đỡ em vượt qua khó khăn sống LỜI CAM ĐOAN Tôi xin cam đoan kết đồ án tự mơ chương trình MCNP5 hướng dẫn thầy cô viện Kỹ thuật hat nhân Vật lý mơi trường Để hồn thành đề tài này, sử dụng tài liệu ghi mục tài liệu tham khảo, ngồi khơng sử dụng tài liệu mà khơng ghi Nếu sai tơi xin chịu hình thức kỷ luật theo qui định Hà Nội, tháng năm 2017 Sinh viên thực Ngô Đắc Lê MỤC LỤC DANH MỤC CÁC BẢNG DANH MỤC HÌNH VẼ VÀ ĐỒ THỊ MỞ ĐẦU……………………… CHƯƠNG 1: CƠ SỞ LÝ THUYẾT VỀ XẠ PHẪU 1.1 Phương pháp xạ trị .9 1.2 Phương pháp xạ phẫu 16 1.3 Các khái niệm liên quan liều xạ phẫu 18 1.4 Cơ sở phép tính phân bố liều 21 CHƯƠNG 2: TỔNG QUAN VỀ THIẾT BỊ XẠ PHẪU GAMMA ART-6000 .24 2.1 Giới thiệu dao gamma 24 2.2 Cấu tạo thiết bị dao gamma ART-6000 [2] 25 2.3 Nguyên lý hoạt động dao gamma ART-6000 31 2.4 Đặc điểm xạ phẫu gamma ART so với xạ trị thông thường 33 CHƯƠNG 3: MÔ PHỎNG MONTE CARLO TRONG THIẾT BỊ XẠ PHẪU GAMMA ART-6000 SỬ DỤNG CHƯƠNG TRÌNH MCNP5 34 3.1 Giới thiệu phương pháp Monte Carlo [8] 34 3.2 Chương trình MCNP[10] 35 3.3 Mô nguồn đơn kênh thiết bị xạ phẫu Gamma ART-6000 .44 3.4 Mô MCNP5 cho 30 nguồn thiết bị xạ phẫu Gamma ART-6000 46 CHƯƠNG 4: KẾT QUẢ VÀ BÀN LUẬN 49 4.1 Kết phân bố liều nguồn đơn kênh 49 4.2 Kết phân bố liều 30 nguồn 53 KẾT LUẬN…………………… .72 TÀI LIỆU THAM KHẢO 73 PHỤ LỤC…………………… 74 DANH MỤC TỪ VIẾT TẮT Từ viết tắt ACTL Tên tiếng Anh The Activation Library Tên tiếng Việt Thư viện hoạt độ chất phóng xạ ART Auto Rotating Treatment Hệ thống điều trị cách quay tự động CT Computed Tomography Ghi hình cắt lớp vi tính CTV Clinical Target Volume Thể tích bia lâm sàng ENDF The Evaluated Nuclear Data File Dữ liệu hạt nhân đánh giá ENDL The Evaluated Nuclear Data Library Thư viện liệu hạt nhân FWHM Full Width Half Maximum Bề rộng nửa đỉnh GTV Gross Tumor Volume Thể tích khối u IV Irradiated Volume Thể tích chiếu xạ LINAC Linear Particle Accelerator Máy gia tốc thẳng MCNP Monte Carlo-N Particle Phương pháp Monte Carlo mô hạt OR Organ at Risk Các tổ chức nhạy cảm PTV Planning Target Volume Thể tích lập kế hoạch RGS Rotating Gamma Knife System Hệ thống dao gamma quay SAD Source to Axis Distance Khoảng cách từ nguồn đến tâm điều trị SSD Source to Surface Distance Khoảng cách từ nguồn đến bề mặt da TV Treatment Volume Thể tích điều trị DANH MỤC CÁC BẢNG Bảng Chỉ tiêu kỹ thuật hệ thống thiết bị Gamma ART-6000 .25 Bảng 2 So sánh đặc điểm xạ phẫu với xạ trị thông thường 33 Bảng Các loại tally 39 Bảng Một số vật liệu dùng mô .48 Bảng So sánh kết FWHM[12] .64 Bảng So sánh đường đồng liều 80% 65 Bảng So sánh đường đồng liều 90% 66 Bảng 4 So sánh khoảng cách đường đồng liều 90% FWHM 67 Bảng So sánh FWHM dọc theo trục Ox shot chiếu .68 Bảng So sánh FWHM dọc theo trục Oy shot chiếu 69 Bảng So sánh FWHM dọc theo trục Oz shot chiếu 70 DANH MỤC HÌNH VẼ VÀ ĐỒ THỊ Hình 1 Hiệu ứng quang điện Hình Hiệu ứng Compton 10 Hình Quá trình sinh cặp hủy cặp 12 Hình Đường cong sống sót tế bào động vật có vú theo liều hấp thụ 14 Hình Nguyên lý phương pháp xạ phẫu 16 Hình Sơ đồ biểu diễn thể tích bia kỹ thuật xạ trị .19 Hình Cách tính liều sâu phần trăm .22 Hình Liều sâu phần trăm mơi trường nước kích cỡ khác ống chuẩn trực [9] 23 Hình Dao Gamma ART-6000TM 25 Hình 2 Hệ thống dao gamma quay .27 Hình Viên nang nguồn Co-60 27 Hình Khối nguồn 28 Hình Khối ống chuẩn trực 29 Hình Hệ thống ống chuẩn trực 29 Hình Cửa che chắn 30 Hình Giường điều trị 30 Hình Khung định vị lập thể 31 Hình 10 Giá đỡ khung đầu 31 Hình 12 Nguyên lý hoạt động dao gamma ART-6000 32 Hình 13 Chùm tia hội tụ 32 Hình Giao diện sử dụng chương trình MCNP5 36 Hình Minh họa tally FMESH .40 Hình 3 Hình học 3D 43 Hình Hình học đơn giản nguồn đơn kênh mơ MCNP 44 Hình Mơ hình chi tiết nguồn đơn kênh thiết bị xạ phẫu Gamma ART-6000 44 Hình Mặt cắt nguồn đơn kênh dọc theo mặt phẳng Oxy .45 Hình Phantom plastic cầu 45 Hình Phân bố nguồn Co-60 thiết bị gamma ART-6000 46 Hình Biểu diễn góc phương vị nguồn so với mặt phẳng Oxy 47 Hình 10 Hình ảnh mặt cắt 10 nguồn chương trình MCNP5 .47 Hình Phân bố liều dọc theo trục x collimator 4mm .49 Hình Phân bố liều tương đối mặt phẳng Oxz collimator mm .49 Hình Phân bố liều tương đối mặt phẳng Oxz collimator mm .50 Hình 4 Phân bố liều dọc theo trục x collimator 8mm .50 Hình Phân bố liều dọc theo trục x collimator 8mm .51 Hình Phân bố liều mặt phẳng Oxz collimator 8mm .51 Hình Phân bố liều dọc theo trục Ox collimator 14mm 52 Hình Phân bố liều dọc theo trục y collimator 18mm .52 Hình Phân bố liều dọc theo trục Ox collimator 4mm 53 Hình 10 Phân bố liều dọc theo trục Oy collimator 4mm 54 Hình 11 Phân bố liều dọc theo trục Oz collimator 4mm 54 Hình 12 Phân bố liều mặt phẳng Oxy collimator 4mm 55 Hình 13 Phân bố liều theo mặt phẳng Oxz collimator 4mm 55 Hình 14 Phân bố liều dọc theo trục Ox collimator 8mm 56 Hình 15 Phân bố liều dọc theo trục Oy collimator 8mm 56 Hình 16 Phân bố liều dọc theo trục Oz collimator 8mm 57 Hình 24 Phân bố liều theo mặt phẳng Oxy collimator 8mm 57 Hình 18 Phân bố liều theo mặt phẳng Oxz collimator 8mm 58 Hình 19 Phân bố liều dọc theo trục Ox collimator 14mm 58 Hình 20 Phân bố liều dọc theo trục Oy collimator 14mm 59 Hình 21 Phân bố liều dọc theo trục Oz collimator 14mm 59 Hình 22 Phân bố liều theo mặt phẳng Oxy collimator 14mm 60 Hình 23 Phân bố liều theo mặt phẳng Oxz collimator 14mm 60 Hình 24 Phân bố liều dọc theo trục Ox collimator 18mm 61 Hình 25 Phân bố liều dọc theo trục Oy collimator 18mm 61 Hình 26 Phân bố liều dọc theo trục Oz collimator 18mm 62 Hình 27 Phân bố liều theo mặt phẳng Oxy collimator 18mm 62 Hình 28 Phân bố liều theo mặt phẳng Oxz collimator 18mm 63 Hình 29 So sánh liều tương đối trục Ox 68 Hình 30 So sánh liều tương đối trục Oy 69 Hình 31 So sánh liều tương đối trục Oz 70 MỞ ĐẦU Theo thống kê tổ chức y tế giới, tỷ lệ tử vong giới bệnh ung thư cao Hàng năm có khoảng gần 10 triệu trường hợp mắc ung thư triệu người chết bệnh Ở Việt Nam, năm ước tính có khoảng 150000 ca ung thư có 50000 ca tử vong.[1] Với phát triển khoa học kĩ thuật thập kỷ gần đây, thiết bị chuẩn đoán điều trị ung thư có tiến vượt bậc Vì mà chuẩn đốn điều trị xác khối u Hiện Việt Nam thiết bị chuẩn đoán điều trị tia xạ đưa vào sử dụng phổ biến bệnh viện thiết bị chuẩn đốn đồng vị phóng xạ: PET, SPECT, CT, Gamma Camera thiết bị điều trị xạ ion hóa đại máy gia tốc tuyến tính Đặc biệt gần bệnh viện Bạch Mai đưa vào máy xạ phẫu gamma ART-6000, coi thiết bị tiên tiến để chữa trị u não Cho đến có 500000 bệnh nhân lựa chọn điều trị Gamma Knife toàn giới.[2] Đồ án nhằm mục đích tìm hiểu sâu thiết bị xạ phẫu gamma ART6000 cấu tạo, nguyên lý hoạt động, ưu điểm xạ phẫu gamma knife so với phương pháp xạ trị thơng thường Qua việc tìm hiểu cấu tạo cách xếp phân bố nguồn Co-60 thiết bị xạ phẫu, chương trình mơ xây dựng để tính tốn phân bố liều kết so sánh với kết nghiệm thu thiết bị nhằm kiểm tra tính đắn q trình mơ Chương trình dùng để mô đồ án MCNP5, phương pháp Monte Carlo, xem xác đại việc tính tốn phân bố liều Đề tài "Mơ phân bố liều thiết bị xạ phẫu Gamma ART-6000 sử dụng phương pháp Monte Carlo" mở hướng nghiên cứu việc ứng dụng phương pháp MCNP5 kỹ thuật tính liều thiết bị xạ phẫu Gamma ART-6000 Quá trình nghiên cứu việc tìm hiểu tổng quan cấu tạo nguyên tắc hoạt động thiết bị xạ phẫu Gamma ART-6000, tìm hiểu đặc điểm, sở vật lý code, cấu trúc thành phần file input chương trình MCNP5 Sau vận dụng phương pháp để mơ q trình tính liều cho thiết bị xạ phẫu Gamma ART-6000 Vì thời gian thực đề tài có hạn nên tơi mơ q trình tính liều phân bố liều theo độ sâu phantom vật liệu plastic-A150 cho nguồn đơn kênh 30 nguồn, shot chiếu khác tâm điều trị Từ mục đích nội dụng cơng việc đồ án có bố cục sau : Chương 1: Các kiến thức liên quan tới xạ trị xạ phẫu Giới thiệu kĩ thuật tính liều xạ trị Các kiến thức tảng sở để xây dựng q trình mơ tương tác xạ với môi trường vật chất, định nghĩa liều chiếu, kích thước trường chiếu, hướng chiếu liều hấp thụ Chương 2: Tổng quan cấu tạo, nguyên lý hoạt động ưu điểm xạ phẫu gamma ART so với xạ trị thông thường Chương 3: Giới thiệu phương pháp Monte Carlo Nội dung chương bao gồm: điều kiện để mơ chương trình MCNP, chương trình MCNP5, thư viện liệu sử dụng, cấu trúc chương trình MCNP5, tally MCNP5, đánh giá sai số chương trình, ưu điểm MCNP5 Mô MCNP5 cho thiết bị xạ phẫu Gamma ART-6000 nguồn đơn kênh 30 nguồn thiết bị xạ phẫu gamma ART-6000 Chương 4: Đưa kết phân bố liều nguồn đơn kênh thiết bị xạ phẫu gamma ART-6000 Biện luận kết thu so sánh kết thu kết nghiệm thu thiết bị năm 2003 Bên cạnh tốn mơ phân bố liều sử dụng shot chiếu kích thước khác vị trí điều trị Bảng So sánh đường đồng liều 80% Đường đồng liều 80% (mm) Đường kính collimator Trục Kết tính Giá trị trung Kết nghiệm tốn bình thu(2003) 3.73 ± 0.72 4.4 6.73 ± 0.87 6.6 14.73 ± 1.08 13.70 18.13 ± 1.10 17.00 4mm Ox 3.80 ± 0.74 4mm Oy 3.60 ± 0.74 4mm Oz 3.80 ± 0.67 8mm Ox 7.80 ± 0.73 8mm Oy 6.80 ± 0.74 8mm Oz 5.60 ± 1.02 14mm Ox 15.60 ±0.67 14mm Oy 15.40 ±1.17 14mm Oz 13.20 ± 1.25 18mm Ox 18.80 ±0.70 18mm Oy 18.20 ±1.32 18mm Oz 17.40 ±1.28 65 Bảng So sánh đường đồng liều 90% Đường đồng liều 90% (mm) Đường kính collimator Trục Kết tính Giá trị trung Kết nghiệm tốn bình thu(2003) 2.47 ± 0.72 3.10 4.5 ± 0.87 4.55 11.33 ± 1.08 11.30 13.27 ± 1.10 13.95 4mm Ox 2.60 ± 0.74 4mm Oy 2.20 ± 0.74 4mm Oz 2.60 ± 0.67 8mm Ox 6.80 ± 0.73 8mm Oy 4.80 ± 0.74 8mm Oz 1.80 ± 1.02 14mm Ox 13.80 ± 0.67 14mm Oy 13.80 ± 1.17 14mm Oz 6.40 ± 1.25 18mm Ox 17.40 ± 0.70 18mm Oy 16.00 ± 1.32 18mm Oz 6.40 ± 1.28 66 Bảng 4 So sánh khoảng cách đường đồng liều 90% FWHM Đường kính Khoảng cách đường đồng liều 90% tới FWHM (mm) collimator Giá trị trung bình Kết nghiệm thu(2003) 4mm 2.03 ± 0.72 1.9 mm 3.15 ± 0.87 2.65 14mm 3.1 ± 1.08 3.35 18mm 4.50 ± 1.10 4.05 Căn vào kết bảng 4.2, 4.3, 4.4 ta thấy kết FWHM trung bình luận văn thu được, kết nghiệm thu (2003) có sai khác, đường đồng liều 80% đường đồng liều 90% tương tự Vấn đề sai lệch chưa gieo đủ số hạt photon cần thiết, khoảng chia mơ chưa đủ nhỏ Bên cạnh đó, cấu hình mơ gần tối ưu đơn giản hóa để tiện cho việc mơ tả cấu trúc hình học nên dẫn tới cấu trúc mô chưa sát so với thực tế, kết mơ bị sai lệch Đường đồng liều 80% thu từ bảng 4.3 phù hợp với tiêu chuẩn nhà sản xuất đưa Sai số lớn 1,27 mm collimator mm, sai số nhỏ 0,13 mm collimator 18mm (sai số không vượt 1,5 mm so với kích thước collimator danh định) Khoảng cách đường đồng liều 90% 50% từ bảng 4.5 nhỏ tiêu chuẩn nhà sản xuất (6mm) Khoảng cách nhỏ collimator mm (2,03 mm) khoảng cách lớn collimator 18 mm (4,05 mm) Khi tăng kích thước collimator vùng điều trị tổn thương mở rộng Do kết hợp nhiều shot chiếu ta thu hiệu điều trị cao So sánh shot chiếu Bài toán đặt ra: So sánh phân bố liều chiếu theo cách: Cách 1: 5Gy collimator 8mm Cách 2: 5Gy collimator 14mm Cách 3: Gy collimator 8mm +2 Gy collimator 14mm 67 Giả thiết toán đặt : Phantom sử dụng phantom cầu plastic A-150 có bán kính 16 cm Số hạt photon gieo 2x108 hạt Một số kết thu nhận : Bảng So sánh FWHM dọc theo trục Ox shot chiếu Đường đồng liều Shot chiếu FWHM(mm) 8mm 11.4 7.80 14mm 20.00 15.60 8mm+14mm 15.40 8.20 80% Phân bố liều theo trục x shot chiếu khác 1.00 0.90 0.80 Liều tương đối 0.70 0.60 0.50 0.40 0.30 0.20 0.10 -2.5 -2 -1.5 -1 -0.5 0.5 1.5 2.5 Trục Ox(cm) Collimator 8mm Collimator 14mm Collimator 8mm+ collimator 14mm Hình 29 So sánh liều tương đối trục Ox 68 Bảng So sánh FWHM dọc theo trục Oy shot chiếu Đường đồng liều Shot chiếu FWHM(mm) 8mm 10.00 6.80 14mm 18.60 15.40 8mm+14mm 14.00 7.40 80% Phân bố liều dọc theo trục y shot chiếu khác 1.00 0.90 0.80 Liều tương đối 0.70 0.60 0.50 0.40 0.30 0.20 0.10 - -2.5 -2 -1.5 -1 -0.5 0.5 1.5 2.5 Trục Oy(cm) collimator 8mm collimator 14mm collimator8mm+collimator 14mm Hình 30 So sánh liều tương đối trục Oy 69 Bảng So sánh FWHM dọc theo trục Oz shot chiếu Đường đồng liều Shot chiếu FWHM(mm) 8mm 10.20 5.60 14mm 14.20 13.20 8mm+14mm 12.60 6.40 80% (mm) Phân bố liều theo trục z shot chiếu khác 1.00 0.90 0.80 Liều tương đối 0.70 0.60 0.50 0.40 0.30 0.20 0.10 - -2.5 -2 -1.5 -1 -0.5 0.5 1.5 2.5 Trục Oz(cm) collimator 8mm collimator 14mm collimator 8mm+collimator 14mm Hình 31 So sánh liều tương đối trục Oz Nhận xét : - Phân bố liều hấp thụ đối xứng theo trục Ox, Oy, Oz, tập trung vùng trung tâm - Dựa vào độ nhọn đỉnh cho ta thấy ứng với collimator 8mm, phân bố liều tập trung xung quanh tâm nên điều trị cho tổn thương có kích 70 thước nhỏ Trong đó, với trường hợp sử dụng collimator 14mm, điều trị tổn thương có kích thước lớn Đối với trường hợp sử dụng kết hợp collimator 8mm collimator 14mm ta điều trị kích thước khối u lớn trường hợp sử dụng collimator 8mm nhỏ trường hợp sử dụng collimator 14mm để điều trị - Trong thực tế, tùy theo phức tạp cấu trúc hình dạng tổn thương mà bác sỹ kỹ sư vật lý cần phải sử dụng linh hoạt loại collimator khác nhau, kết hợp nhuần nhuyễn để lập kế hoạch điều trị tốt 71 KẾT LUẬN Việc tính tốn liều xạ phẫu toán phức tạp tính chất phức tạp tương tác xạ với vật chất tính chất phức tạp cấu hình thực tế Trong năm gần việc nghiên cứu sử dụng code MCNP trình mơ tính tốn liều kĩ thuật xạ phẫu trở nên phổ biến ưu điểm chương trình mang lại Đồ án với mục đích tìm hiểu giới thiệu phương pháp tính chương trình MCNP5 cung cấp cho kỹ sư vật lý tài liệu tham khảo bổ ích Đồng thời xây dựng tốn tính tốn phân bố liều thiết bị xạ phẫu gamma ART-6000 chương trình MCNP5 để kiểm tra tính đắn q trình tính tốn so với thực tế.[8] Qua việc mơ phân bố liều nguồn đơn kênh, 30 nguồn thiết bị gamma ART-6000, cho thấy kết phân bố liều thu từ chương trình MCNP5 phù hợp so với thực tế Điều mở hướng nghiên cứu ứng dụng chương trình MCNP5 việc lập kế hoạch xạ phẫu cho thiết bị gamma ART-6000 Do hạn chế thời gian tính chất phức tạp tốn tính liều, cấu trúc hình học tốn mơ chưa xác tuyệt đối so với cấu trúc hình học thực tế Do vậy, kết tốn phân bố liều nhiều sai lệch so với kết công bố Đồ án áp dụng mô phân bố liều cho phantom plastic dạng hình cầu, shot chiếu với thời gian chiếu vào vị trí tâm phantom Để phát huy vai trò chương trình MCNP5 tính tốn phân bố liều cho kỹ thuật xạ phẫu thiết bị gamma ART-6000, tơi hy vọng tương lai tính tốn phân bố liều cụ thể cho vị trí khối u thực tế, mơ tính liều cho khối u bất kì, với nhiều shot chiếu khác khơng trùng vị trí 72 TÀI LIỆU THAM KHẢO [1] N Đ T Nguyễn Thái Hà, “Y học hạt nhân kỹ thuật xạ trị.” [2] Phạm Thị Tuyết, “QA tính tốn liều máy Gamma Knife quay phương pháp mô Monte Carlo.” 2016 [3] Ngô Quang Huy, “Cơ sở vật lý hạt nhân.” [4] Phùng Văn Duân, “An toàn xạ bảo vệ môi trường.” 2006 [5] Nguyễn Thị Phương An, “Áp dụng chương trình penelope để mơ phân bố liều xạ trị dao gamma leksell,” 2011 [6] Quản Tùng Lâm, “Quy trình kỹ thuật xạ phẫu dao gamma Gyro knife bệnh viện C Thái Nguyên.” [7] Ngơ Quang Huy, “An tồn xạ ion hóa.” 2004 [8] Đặng Trương Ka My, “Mơ thiết bị xạ phẫu Leksell Gamma Knife chương trình MCNP5,” 2009 [9] Arndt Jürgen, “Gamma Knife Dosimetry & Treatment Planning,” 1999 [10] M A, M Carlo, and X.-M C Team, “MCNP-A General Monte Carlo NParticle Transport Code Version Volume I Overview and Theory,” vol 836, 2003 [11] Đặng Nguyên Phương, “Hướng dẫn sử dụng MCNP cho hệ điều hành Windows,” 2015 [12] I American Radiosurgery, “Acceptance test procedure Gamma ART-6000.” 2003 73 PHỤ LỤC Phụ lục A : Chương trình mơ MCNP cho nguồn đơn kênh c cell cards 10 -0.001205 -6 -2 imp:p=1 $ khong 20 -19.35 -9 -2 imp:p=1 $ tungsten 10 -0.001205 -6 -3 imp:p=1 $ khong 30 -11.35 -9 -3 imp:p=1 $ chi 10 -0.001205 -6 -4 imp:p=1 $ khong 20 -19.35 -9 -4 imp:p=1 $ tungsten 40 -8.9 -8 -5 imp:p=1 $ Co-60 50 -7.87 -9 -5 imp:p=1 $ Fe 10 -0.001205 -7 -5 imp:p=1 $ khong 10 60 -1.127 12 -13 14 -15 16 -17 imp:p=1 $ phantom 90 10 -0.001205 #10 -90 #(-9 -5) imp:p=1 91 90 imp:p=0 c surface cards pz 16.5 pz 22.5 pz 33 pz 39 pz 41 TRC 0 39 0 -22.5 0.025025 0.19 $ collimator 8mm cz 0.05005 cz 0.05 cz 12 so 90 so 50 c data cards 74 mode p sdef pos=0 40 cell=7 par=2 erg=1.25 fmesh4:p geom=xyz origin=-4 -0.2 -0.2 & imesh=4 iints=100 & jmesh=0.2 jints=1 & kmesh=0.2 kints=1 & out=ik m10 7000 -0.755268 8000 -0.231781 18000 -0.012827 6000 -0.000124 $ khong m20 74000 $ tungsten m30 82000 $ chi m40 29000 $ cobalt m50 26000 $Fe m60 1000 -0.10133 6000 -0.775498 7000 -0.035057 & 8000 -0.052315 9000 -0.167415 20000 -0.176585 nps 1000000000 75 Phụ lục B Chương trình mơ MCNP5 cho 30 nguồn c cell card 10 -0.001205 -6 -4 u=7 imp:p=1 $khong 20 -19.35 -2 u=7 imp:p=1 $tungsten 30 -11.35 -3 u=7 imp:p=1 $chi 20 -19.35 -4 u=7 imp:p=1 $tungsten 40 -8.9 -8 u=7 imp:p=1 $Co60 50 -7.87 u=7 imp:p=1 $Fe 10 -0.001205 -7 u=7 imp:p=1 $khong -5 -9 fill=7 imp:p=1 11 like but trcl=(0 0 0.97437 -0.22495 0.22495 0.97437 0 1) 12 like but trcl=(0 0 0.93667 -0.35021 0.35021 0.93667 0 1) 13 like but trcl=(0 0 0.88295 -0.46947 0.46947 0.88295 0 1) 14 like but trcl=(0 0 0.81412 -0.5807 0.5807 0.81412 0 1) 15 like but trcl=(0 0 0.73135 -0.682 0.682 0.73135 0 1) 16 like but trcl=(0 0 -0.97437 -0.22495 0.22495 -0.97437 0 1) 17 like but trcl=(0 0 -0.93667 -0.35021 0.35021 -0.93667 0 1) 18 like but trcl=(0 0 -0.88295 -0.46947 0.46947 -0.88295 0 1) 19 like but trcl=(0 0 -0.81412 -0.5807 0.5807 -0.81412 0 1) 20 like but trcl=(0 0 -0.73135 -0.682 0.682 -0.73135 0 1) 21 like but trcl=(0 0 0.48719 -0.22495 0.84383 0.11248 & 0.97437 0.19481 -0.86603 0.5) 22 like but trcl=(0 0 0.46834 -0.35021 0.81118 0.1751 & 0.93667 0.30329 -0.86603 0.5) 23 like but trcl=(0 0 0.44147 -0.46947 0.76466 0.23474 & 0.88295 0.40657 -0.86603 0.5) 24 like but trcl=(0 0 0.40706 -0.5807 0.70504 0.29035 & 0.81412 0.5029 -0.86603 0.5) 25 like but trcl=(0 0 0.36568 -0.682 0.63337 0.341 0.73135 & 76 0.59063 -0.86603 0.5) 26 like but trcl=(0 0 -0.48719 -0.22495 -0.84383 0.11248 & -0.97437 0.19481 -0.86603 0.5) 27 like but trcl=(0 0 -0.46834 -0.35021 -0.81118 0.1751 & -0.93667 0.30329 -0.86603 0.5) 28 like but trcl=(0 0 -0.44147 -0.46947 -0.76466 0.23474 & -0.88295 0.40657 -0.86603 0.5) 29 like but trcl=(0 0 -0.40706 -0.5807 -0.70504 0.29035 & -0.81412 0.5029 -0.86603 0.5) 30 like but trcl=(0 0 -0.36568 -0.682 -0.63337 0.341 & -0.73135 0.59063 -0.86603 0.5) 31 like but trcl=(0 0 -0.48719 -0.22495 0.84383 -0.11248 & 0.97437 0.19481 -0.86603 -0.5) 32 like but trcl=(0 0 -0.46843 -0.35021 0.81118 -0.1751 & 0.93667 0.30329 -0.86603 -0.5) 33 like but trcl=(0 0 -0.44147 -0.46947 0.76466 -0.23474 & 0.88295 0.40657 -0.86603 -0.5) 34 like but trcl=(0 0 -0.40706 -0.5807 0.70504 -0.29035 & 0.81412 0.5029 -0.86603 -0.5) 35 like but trcl=(0 0 -0.36568 -0.682 0.6337 -0.341 & 0.73135 0.59063 -0.86603 -0.5) 36 like but trcl=(0 0 0.48719 -0.22495 -0.84383 -0.11248 & -0.97437 0.19481 -0.86603 -0.5) 37 like but trcl=(0 0 0.46843 -0.35021 -0.81118 -0.1751 & -0.93667 0.30329 -0.86603 -0.5) 38 like but trcl=(0 0 0.44147 -0.46947 -0.76466 -0.23474 & -0.88295 0.40657 -0.86603 -0.5) 39 like but trcl=(0 0 0.40706 -0.5807 -0.70504 -0.29035 & -0.81412 0.5029 -0.86603 -0.5) 77 40 like but trcl=(0 0 0.36568 -0.682 -0.6337 -0.341 & -0.73135 0.59063 -0.86603 -0.5) 60 -1.127 -12 imp:p=1 $ phantom 500 10 -0.001205 #9 #8 -500 #11 #12 #13 #14 #15 & #16 #17 #18 #19 #20 #21 #22 #23 #24 #25 #26 #27 #28 #29 #30 #31 #32 #33 " #35 #36 #37 #38 #39 #40 imp:p=1 501 500 imp:p=0 c surface cards px 16.5 px 22.5 px 33 px 39 px 41 kx 41.055 5.93E-4 -1 $ collimator 18mm cx 0.05005 cx 0.05 cx 0.6 12 so 500 so 50 c data cards mode p sdef pos=d1 par=2 erg=1.25 si1 L 38.9748 -8.998 37.4669 -14.0083 35.3179 -18.7789 & 32.5646 -23.2281 29.2542 -27.2799 -38.9748 -8.9980 & -37.4669 -14.0083 -35.3179 -18.7789 -32.56462 -23.2281 & -29.2542 -27.2799 19.6305 -9.0746 34 18.5052 -13.8058 32 & 17.3201 -18.423 30 16.1637 -23.0575 28 14.4351 -26.919 25 & 78 -19.6305 -9.0746 -34 -18.5052 -13.8058 -32 -17.3201 -18.423 -30 & -16.1637 -23.0575 -28 -14.4351 -26.919 -25 -19.6305 -9.0746 34 & -18.5052 -13.8058 32 -17.3201 -18.423 30 -16.1637 -23.0575 28 & -14.4351 -26.919 25 19.6305 -9.0746 -34 18.5052 -13.8058 -32 & 17.3201 -18.423 -30 16.1637 -23.0575 -28 14.4351 -26.919 -25 sp1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 fmesh4:p geom=xyz origin=-2.5 -0.5 -0.5 & imesh=2.5 iints=1000 & jmesh=0.5 jints=1 & kmesh=0.5 kints=1 & out=ij m10 7000 -0.755268 8000 -0.231781 18000 -0.012827 6000 -0.000124 $ khong m20 74000 $ tungsten m30 82000 $ chi m40 29000 $ cobalt m50 26000 $Fe m60 1000 -0.10133 6000 -0.775498 7000 -0.035057 & 8000 -0.052315 9000 -0.167415 20000 -0.176585 nps 1000000000 79 ... 65 Bảng So sánh đường đồng liều 90% 66 Bảng 4 So sánh khoảng cách đường đồng liều 90% FWHM 67 Bảng So sánh FWHM dọc theo trục Ox shot chiếu .68 Bảng So sánh FWHM dọc theo trục... chịu hình thức kỷ luật theo qui định Hà Nội, tháng năm 2017 Sinh viên thực Ngô Đắc Lê MỤC LỤC DANH MỤC CÁC BẢNG DANH MỤC HÌNH VẼ VÀ ĐỒ THỊ MỞ ĐẦU……………………… CHƯƠNG... 2 So sánh đặc điểm xạ phẫu với xạ trị thông thường 33 Bảng Các loại tally 39 Bảng Một số vật liệu dùng mô .48 Bảng So sánh kết FWHM[12] .64 Bảng So sánh đường đồng