ứng dụng phương pháp monte carlo trong xạ trị ung thư phổi bằng máy gia tốc tuyến tính

ĐẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN KHOA VẬT LÝ-VẬT LÝ KỸ THUẬT BỘ MÔN VẬT LÝ HẠT NHÂN ……… o0o………… KHÓA LUẬN TỐT NGHIỆP ĐẠI HỌC Đề tài: ỨNG DỤNG PHƯƠNG PHÁP MONTE CARLO TÍNH LIỀU HẤP THỤ TRONG XẠ TRỊ UNG THƯ PHỔI BẰNG MÁY GIA TỐC TUYẾN TÍNH SVTH: TRẦN VĂN NGHĨA CBHD: ThS NGUYỄN THỊ CẨM THU ThS LÊ THANH XUÂN GVPB: ThS NGUYỄN TẤN CHÂU TP HỒ CHÍ MINH - 2014 MỤC LỤC DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT a DANH MỤC BẢNG VÀ HÌNH c LỜI MỞ ĐẦU CHƯƠNG 1: TỒNG QUAN VỀ UNG THƯ VÀ PHƯƠNG PHÁP XẠ TRỊ 1.1 Tổng quan ung thư xạ trị 1.1.1 Ung thư 1.1.2 Các phương pháp chẩn đoán điều trị ung thư 1.2 Phương pháp xạ trị 1.2.1 Xạ trị (Brachytherapy) 1.2.2 Xạ trị (Teletherapy) CHƯƠNG 2: MÁY GIA TỐC TUYẾN TÍNH VÀ CHƯƠNG TRÌNH MCNP5 10 2.1 Máy gia tốc 10 2.1.1 Lịch sử hình thành phát triển máy gia tốc 10 2.1.2 Cấu hình máy gia tốc dùng xạ trị 11 2.1.3 Nguyên lý hoạt động 12 2.2 Chương trình MCNP5 14 2.2.1 Khái quát 14 2.2.2 Các đặc trưng MCNP 15 2.2.2.1 Số liệu hạt nhân 15 2.2.2.2 Các đặc trưng nguồn 15 2.2.2.3 Đánh giá sai số 16 2.2.2.4 Giảm sai số 17 2.2.2.5 Kết toán 17 CHƯƠNG 3: LẬP KẾ HOẠCH XẠ TRỊ UNG THƯ PHỔI DỰA TRÊN PHẦN MỀM PROWESS PANTHER 18 3.1 Quy trình xạ trị 18 3.2 Lập kế hoạch điều trị ca ung thư phổi 19 3.2.1 Yêu cầu 19 3.2.2 Quy trình lập kế hoạch 19 CHƯƠNG 4: MÔ PHỎNG MÁY GIA TỐC TUYẾN TÍNH BẰNG CHƯƠNG TRÌNH MCNP5 27 4.1 Mục đích 27 4.2 Cấu hình đầu máy gia tốc PRIMUS mơ MCNP5 27 4.2.1 Cấu hình đầu máy gia tốc PRIMUS HPD 27 4.2.2 Mơ hình đầu máy gia tốc MCNP5 29 4.2.3 Khai báo nguồn 29 4.2.4 Cấu hình Phantom 30 4.2.5 Các thông số mô cho MCNP5 31 4.3 Đánh giá kết tính tốn phân bố liều phantom CT 33 KẾT LUẬN VÀ KIẾN NGHỊ 42 TÀI LIỆU THAM KHẢO 44 DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT Các ký hiệu FOM: số chất lượng tính tốn mơ N: số lần tính R: tỉ số thăng giáng chuẩn giá trị trung bình T: thời gian tính (s) : thăng giáng chuẩn : giá trị trung bình phép đo Danh mục chữ viết tắt 2D: Two – Dimensional 3D: Three – Dimensional ADN: Acid Deoxyribo Nucleic BEV: Beam Eye View CTV: Clinical Target Volume CT: Computed Tomography DVH: Dose Volume Histogram DICOM: Digital Imaging Communications Medicine GTV: Gross Tumor Volume MRI: Magentic Resonance System MCNP: Monte Carlo N – Particle a MCNG: Monte Carlo Neutron - Gamma MCP: Monte Carlo Photon PET: Positron Emission Tomography PP: Prowess Panther PTV: Planning Target Volume SIM: Simulator SPECT: Single Photon Emission Computed Tomography b DANH MỤC BẢNG VÀ HÌNH Danh mục bảng Bảng 2.1: Ý nghĩa giá trị R .………………………………………………… 18 Bảng 3.1: Thống kê liều hấp thụ quan theo bảng DVH ……………… 28 Bảng 4.1: Khoảng lượng xác suất phát chùm electron để tạo photon lượng 15 MV ……………………………………………………………………… 32 Bảng 4.2: Tọa độ liều tương đối tính MCNP ………………………… 41 Danh mục hình Hình 1.1: Phổi bị ung thư …………………………………………………………….4 Hình 1.2: Ung thư dày…… …………………………………………………… Hình 1.3: Xạ trị ung thư vú giai đoạn sớm xạ trị áp xát ………………………8 Hình 1.4: Xạ trị áp sát trị ung thư vòm họng… Hình 1.5: Xạ trị ngồi máy Colbalt 60 ……………………………………… 10 Hình 2.1: Các phận máy gia tốc xạ trị ……………………………… 12 Hình 3.1: Biểu tượng nút chọn cài dặt tâm chì thiết kế phần mềm Prowess Panther .…………………………………… 22 Hình 3.2: Hình ảnh Prowess Panther hiển thị sau chọn Cross section of marker ……………………………22 Hình 3.3: Đường contour bao quanh khối u quan lành: đường màu đỏ bao quanh khối u, đường màu vàng bao quanh tủy sống, đường màu cam bao quanh phổi phải, c đường màu xanh dương bao quanh tim, đường màu hồng nhạt bao quanh phổi trái …………………… 23 Hình 3.4: Hướng trường chiếu ứng với góc quay: 200 hướng từ suống chếch sang phải, 1950 hướng từ lên chếch sang trái, 2700 hướng từ trái qua ……………… 24 Hình 3.5: Hình dạng trường chiếu nêm tương ứng ………………… … 24 Hình 3.6: Khối che (block) ứng với hai trường chiếu, (a) hình dạng khối che chắn ứng với trường chiếu 200, (b) hình dạng khối che chắn ứng với trường chiếu 1950 .……………… ……………23 Hình 3.7: Khối che chắn ứng với trường chiếu 2700 24 Hình 3.8: Dạng hình học 3D khối u quan lành ………………………… 26 Hình 3.9: Bảng DVH ……………………………………………………………….27 Hình 4.1: Mơ hình máy gia tốc 2D vẽ Visual Editor MCNP5 phát chùm photon 15 MV …………………………………………………………………… 31 Hình 4.2: Kích thước voxel phantom thu từ ảnh CT ……………… ….33 Hình 4.3: Giao diện chuyển đổi chương trình CODIM (CONVERT VIEW) 36 Hình 4.4: Kết tính liều theo hướng chiếu MCNP5 vẽ chương trình CODIM (ISODOSE MESHTAL) cho mặt CT trung tâm 38 Hình 4.5: Kết tính liều theo hướng chiếu MCNP5 vẽ chương trình CODIM (ISODOSE MESHTAL) cho mặt CT trung tâm +1 (cách slice trung tâm mm) .39 Hình 4.6: Phân bố đường đồng liều mô MCNP5 .40 d Hình 4.7: Phân bố đường đồng liều mô Prowess Panther .40 e LỜI MỞ ĐẦU Bệnh ung thư sức khỏe cộng đồng vấn đề ngày quan tâm hầu hết quốc gia giới Theo số thống kê Cơ quan quốc tế ung thư tổ chức Y tế Thế giới (WHO) dự báo, số bệnh nhân mắc bệnh ung thư phát toàn giới vào năm 2035 lên tới xấp xỉ 24 triệu người, tăng lên 70% so với năm 2012 14,1 triệu người Trong đó, số ca tử vong ung thư tăng gần gấp đôi từ 8,2 triệu người năm 2012 lên tới 14,6 triệu người vào năm 2035 Riêng với nước nghèo, chậm phát triển, số ca mắc bệnh ung thư tăng 40% năm, nước giàu có 20% năm [12] Ở Việt Nam, dẫn theo nguồn nghiên cứu Viện Nghiên cứu phòng chống ung thư Việt Nam đăng báo Lao Động số 27, đăng ngày tháng năm 2014 nước có từ 240 nghìn đến 250 nghìn người mắc bệnh ung thư Mỗi năm, có 150 nghìn ca mắc bệnh có từ 75 nghìn tới 85 nghìn người tử vong bệnh Một số đáng báo động có xu hướng tăng lên theo năm Ung thư chữa trị nhiều phương pháp khác tập trung chủ yếu vào ba phương pháp là: xạ trị, hóa trị phẫu thuật Trong đó, xạ trị phương pháp đem lại hiệu cao sử dụng phổ biến nước ta Việc ứng dụng phần mềm mơ để tính tốn liều hấp thụ khối u quan lành thiếu bệnh viện xạ trị Vì thời gian lập kế hoạch có giới hạn nên phương pháp Monte Carlo thường không áp dụng bệnh viện Tuy nhiên, phương pháp Monte Carlo cho kết có độ tin cậy cao… Với ý nghĩa quan trọng đó, khóa luận thực việc tính liều hấp thụ lên khối u phổi xạ trị máy gia tốc tuyến tính chương trình mơ Monte Carlo MCNP5 so sánh kết mô với kết có từ chương trình lập kế hoạch bệnh viện Với nội dung trên, khóa luận trình bày gồm chương sau: Chương 1: Tổng quan ung thư phương pháp xạ trị Chương 2: Giới thiệu máy gia tốc tuyến tính chương trình MCNP5 Chương 3: Lập kế hoạch xạ trị ung thư phổi dựa phần mềm Prowess Panther Chương 4: Mô máy gia tốc tuyến tính chương trình MCNP5 voxel, voxel phantom CT có kích thước kích thước pixel nhân với độ dày lát cắt CT (thông tin thu từ phần header file hình ảnh CT) Hình 4.2 Kích thước voxel phantom thu từ ảnh CT Hình ảnh CT đem mơ hình ảnh thu nhận từ bệnh viện Đa khoa Kiên Giang thực máy chụp cắt lớp CT SIEMENS Hình ảnh cấu thành từ 512 x 512 điểm ảnh (pixel) Kích thước pixel 0,8945 x 0,8945 mm2 độ dày lát cắt CT mm (xem Hình 4.2) 4.2.5 Các thơng số mô cho MCNP5  Tally Tally phương thức tính tốn MCNP cung cấp để thu kết ghi nhận hạt truyền mơi trường vật chất Tùy vào mục đích tính tốn khác mà ta lựa chọn sử dụng loại tally tương ứng Các phiên 31 MCNP cũ cung cấp tất loại tally gồm có F1, F2, F4, F5, F6, F7, F8 Bên cạnh phiên MCNP5 cịn cung cấp thêm loại tally tally FMESH Trong phạm vi luận văn này, cần sử dụng hai loại tally F8 FMESH để xác định liều hấp thụ phantom sau [6]:  Tally F8: tally cung cấp phân bố lượng xạ phát thu nhận cell gọi tally độ cao xung Ở luận văn này, tally F8 dùng để tính liều hấp thụ voxel trục trung tâm (phân bố liều hấp thụ theo độ sâu) phantom nước thơng qua việc tính tốn lượng để lại voxel  Tally FMESH: tally tính tốn ước lượng chiều dài vết thơng lượng hạt tạo cell, giá trị tính trung bình mạng lưới cell, đơn vị tally FMESH hạt/cm2 Ở luận văn này, tally FMESH dùng để chia nhỏ phantom nước phantom CT thành mạng lưới voxel nhỏ tính phân bố liều lên voxel cách chuyển đổi từ thông lượng sang liều thông qua hệ số chuyển đổi (flux-to-dose rate conversion factors)  Số hạt (NPS) Quy định số hạt phát cửa thoát máy gia tốc hay nói cách khác số biến cố sử dụng mô Monte Carlo người dùng thiết lập  Thời gian (CTME) Quy định thời gian tính tốn mơ phỏng, thời gian nhiều số lịch sử hạt phát lớn Trong q trình mơ MCNP, chương trình dừng lại hai thông số NPS CTME đạt tới ngưỡng quy định trước  Kĩ thuật giảm phương sai Nhằm mục đích giảm thời gian tính tốn chương trình mà khơng ảnh hưởng đến độ xác toán Trong luận văn áp dụng hai kỹ thuật lượng cắt (energy cut-off) phân chia bremsstrahlung (bremsstrahlung splitting) 32 Do việc mô vận chuyển electron vật chất tốn nhiều thời gian chạy chương trình nên để rút ngắn thời gian, chọn mức lượng cắt electron vào khoảng 0,7 MeV mức lượng cắt photon vào khoảng 0,01 MeV Nếu q trình mơ phỏng, lượng electron photon rơi xuống ngưỡng cắt chương trình MCNP ngưng khảo sát tiếp hạt Đồng thời, áp dụng phương pháp phân chia bremsstrahlung có chương trình MCNP Trong phương pháp này, tương tác phát xạ bremsstrahlung, electron phát N photon thay photon phương pháp mơ thơng thường, photon phát có trọng số 1/N, chọn số photon phát N = 80 [6] 4.3 Đánh giá kết tính tốn phân bố liều phantom CT Phantom CT xây dựng dựa số liệu chụp ảnh CT ca ung thư phổi bệnh nhân cụ thể bệnh viện Đa khoa Kiên Giang Cùng với tiến hành sử dụng chương trình CODIM để đọc số liệu hình ảnh, chuyển đổi liệu thành input file tiến hành tính liều MCNP5 Hình bên mơ tả giao diện làm việc chương trình CODIM 33 Hình 4.3: Giao diện chuyển đổi chương trình CODIM ( CONVERT VIEW )  Kích thước phantom CT Mỗi voxel xây dựng có kích thước 0,08945×0,08945×0,6 cm3 Kích thước phantom xây dựng MCNP5 sau: x: -22,8992 cm đến 22,8992 cm y: 22,8992 cm đến 22,8992 cm z: -70,425 cm đến -116,198 cm  Các thông số thiết lập hướng chiếu Q trình mơ thực theo ba trường chiếu với chùm photon lượng 15 MV ứng với trường, tọa độ vị trí trung tâm X = 0,8 cm, Y = -2,5 cm Z = -90,3 cm với trường chiếu sau: 34 - Trường ứng với góc quay 600 so với phương thẳng đứng từ xuống, với kích thước trường chiếu 11,3 cm x 10,8 cm (tương ứng với độ mở rộng theo hai hướng trục x y) với chất lượng chùm tia thiết lập 100% - Trường hai ứng với góc quay 1540 so với phương thẳng đứng từ xuống, với kích thước trường chiếu 10,3 cm x 10,8 cm (tương ứng với độ mở rộng theo hai hướng trục x y) với chất lượng chùm tia 100% - Trường ba ứng với góc quay 3300 so với phương thẳng đứng từ xuống, với kích thước trường chiếu 10,4 cm x 10,6 cm (tương ứng với độ mở rộng theo hai hướng trục x y) với chất lượng chùm tia 100%  Tính liều so sánh với chương trình Prowess Panther(PP) Q trình tính liều mơ MCNP5 thực hỗ trợ chương trình CODIM dựa thông số cung cấp cách đầy đủ Sau đó, kết mơ đem so sánh với hình ảnh chương trình lập kế hoạch điều trị PP Hình 4.4; 4.5; 4.6 4.7 trình bày kết tính liều MCNP5 cho slice trung tâm, slice cách slice trung tâm mm so sánh đường đồng liều thu sau mô MCNP5 với đường thu từ PP 35 Hình 4.4: Kết tính liều theo hướng chiếu MCNP5 vẽ chương trình CODIM (ISODOSE MESHTAL) cho mặt CT trung tâm 36 Hình 4.5: Kết tính liều theo hướng chiếu MCNP5 vẽ chương trình CODIM (ISODOSE MESHTAL) cho mặt CT trung tâm +1 ( cách slice trung tâm mm) 37 Hình 4.6: Phân bố đường đồng liều mơ MCNP5 Hình 4.7: Phân bố đường đồng liều mô Prowess Panther 38 Hình 4.6 4.7 biểu diễn vùng nhận mức liều khác hai chương trình MCNP5 PP Từ vùng liều cao ứng với khu vực có khối u tới vùng liều thấp xa khối u thể màu sắc khác Ở hình 4.6, để tiện cho việc phân biệt vùng liều, gọi vùng với liều khác với tên: màu đỏ thẩm, màu đỏ tươi, màu đỏ màu cam Vùng màu đỏ đậm (màu đỏ thẩm), có liều tác động cao vào khối u ứng với đường 100% liều đánh dấu điểm O, vùng màu đỏ đậm (màu đỏ tươi) đánh dấu điểm A ứng với đường 93,6% liều, vùng ký hiệu điểm B ứng với vùng màu đỏ ứng với đường 92,5% liều vùng màu cam ký hiệu điểm C ứng với đường 90,1% liều Bên cạnh từ hình 4.6, cho thấy có phần liều gây tác động lên lớp lưới giữ cố định bên từ hai hướng chiếu 600 3300, nguyên nhân MCNP5 tự loại bỏ yếu tố lưới giữ cố định với hướng chiếu 1540 chương trình khơng loại bỏ yếu tố giường nằm Vì mơ khơng tách mật độ vật chất giường lưới với thể bệnh nhân, dẫn tới tính tốn liều phân bố phần lên giường nằm lưới giữ trước gây tác động lên lớp da Bảng 4.2: Tọa độ liều tương đối tính MCNP Phần trăm Tọa độ Liều tương đối Liều tương đối điểm ảnh MCNP (%) Prowess Panther (%) O (297;251) 100 100 A (275;261) 93,6 95,0 2,4 B (264;269) 92,5 92,0 0,5 C (251;278) 90,1 88,0 2,1 Vị trí 39 chênh lệch (%) Bảng 4.2 trình bày suất liều tính chương trình MCNP chương trình Prowess Panther cho bốn điểm mà ta đánh dấu hình 4.6 Vì khơng có đủ thơng tin thông số kỹ thuật máy gia tốc nên giá trị so sánh với chương trình Prowess Panther giá trị suất liều tương đối Ở đây, thấy việc so sánh suất liều tương đối mơ MCNP chương trình Prowess Panther có giá trị phù hợp với nhau, so sánh vùng suất liều độ chênh lệch chương trình nhỏ 2% Bảng 4.2 cho thấy chương trình mơ phân biệt vùng liều 98%, 95%, 92% 88% chương trình Prowess, điều chứng tỏ việc ứng dụng mơ MCNP5 cho trường hợp có độ tin cậy cao Như vậy, phân bố đường đồng liều thu từ hai chương trình MCNP5 Prowess Panther trùng khớp với Vùng nhận liều cao gần trùng khớp hai chương trình Có phân biệt vùng liều khác thông qua màu khác Điểm liều cực đại (max isodose) đảm bảo nằm vào khối u giống tính Prowess Panther Mặc dù cịn hạn chế định, chẳng hạn phân biệt vùng liều khác nhau, phân biệt chưa rõ nét dứt khoát Các vùng liều nhiều chỗ chồng lên nhau, gây liều nhiễu vùng liều, dẫn tới việc chênh lệch suất liều hai chương trình Nguyên nhân dẫn tới hạn chế nêu số nguyên nhân sau: Thăng giáng MCNP5 cao, việc tính tốn liều cịn sai lệch, dẫn tới sai số việc phân tách vùng liều Thuật tốn Monter Carlo địi hỏi thời gian tính tốn lâu thống kê lớn Trong đó, chương trình Prowess Panther (PP) sử dụng liệu bán thực nghiệm phương pháp tất định nên có thăng giáng thống kê tốt Trong q trình tính tốn, chương trình Prowess Panther tự động loại bỏ vật mốc CT bệnh nhân lưới giữ cố định bệnh nhân, mặt nạ, dây chì giường nằm chương trình mơ MCNP5 lại khơng thể loại bỏ Một nguyên nhân kể đến hạn chế chương trình CODIM 40 việc sử dụng số liệu chuyển đổi liều kerma sang liều hấp thụ mô thể khác so với hệ số mà phần mềm Prowess Panther sử dụng Đồng thời, yếu tố cường độ, thời gian phát tia, trường chiếu, kích thước thành phần cấu tạo xác thành phần đầu máy gia tốc có sai khác định nên dẫn đến sai lệch Một vấn đề làm ảnh hưởng không nhỏ đến độ sai lêch thành phần vật chất, mật độ cấu tạo mô thể người hai chương trình khác 41 KẾT LUẬN VÀ KIẾN NGHỊ Trong năm qua, phương pháp mô Monter Carlo nói chung phần mềm mơ MCNP nói riêng nghiên cứu ứng dụng rộng rãi nhiều lĩnh vực khác Trong lĩnh vực Y học Hạt Nhân, việc ứng dụng thành công chương trình mơ MCNP mơ tính tốn liều xạ trị điều trị ung thư đem lại nhiều lợi ích đáng kể ưu điểm mà phương pháp đem đến Thơng qua tìm hiểu cấu hình đầu máy gia tốc tuyến tính PRIMUS sử dụng xạ trị dựa số liệu ca xạ trị ung thư cụ thể Trong khóa luận này, kế thừa chương trình mơ máy gia tốc PRIMUS chương trình MCNP5 chương trình CODIM tác giả Lê Thanh Xuân [6] để mô cho trường hợp lập kế hoạch mô cho bệnh nhân ung thư phổi Bệnh viện Kiên Giang kết thu phù hợp tốt so sánh với chương trình mơ lập kế hoạch Prowess Panther có bệnh viện Những kết đạt được: - Xây dựng hình ảnh chụp CT bệnh nhân thực tế lấy từ bệnh viện Đa khoa Kiên Giang để đưa vào MCNP5 thông qua phần mềm CODIM - Thiết lập xác trường chiếu máy gia tốc lên phantom người phù hợp với trường hợp lập kế hoạch cho bệnh nhân thực tế Bệnh viện - Kết thu q trình mơ sau so sánh với kết thực nghiệm từ bệnh viện Đa khoa Kiên Giang cho thấy kết thu từ mô tốt, có độ tin cậy cao Các đường đồng liều liều nhận vào quan khối u trùng khớp, có sai số bảo đảm tính đắn kết thu 42 Tuy nhiên, bên cạnh kết đạt cịn có số hạn chế định - Chưa linh hoạt việc thiết lập trường chiếu khác mô - Chưa xử lý phần nhiễu thăng giáng mơ MCNP5 chương trình CODIM - Mặc dù áp dụng phương pháp giúp rút ngắn thời gian tính tốn MCNP5 thời gian tính tốn MCNP5 hạn chế, thời gian chạy hồn thành mơ vào khoảng từ tới 10 ngày Vì vậy, sai số thống kê thăng giáng mơ cịn lớn dẫn tới khó phân tách vùng có suất liều cao làm ảnh hưởng tới phân bố suất liều vùng Để cải thiện khó khăn, hạn chế cần phải xây dựng phương án làm giảm thời gian tính tốn Sử dụng số kỹ thuật giảm nhiễu phân bố liều mô MCNP thăng giáng gây nên MCNP5 không hữu hiệu việc mô tính tốn liều xạ trị, mà cịn ứng dụng nhiều lĩnh vực nghiên cứu khác kiểm chứng phép đo chiều dày vật liệu, áp dụng việc tính tốn hiệu suất DET HPGe GEM 15P4,… Trong tương lai, nghiên cứu ứng dụng MCNP có ích cho cơng tác đào tạo nhằm phát triển nguồn nhận lực lĩnh vực ứng dụng kỹ thuật hạt nhân công nghiệp Việt Nam 43 TÀI LIỆU THAM KHẢO Tiếng Việt [1] Phan Sỹ Anh, Trần Xuân Trường, Mai Trọng Khoa, Nguyễn Đắc Nhật, Nguyễn Thị The, Nguyễn Thành Chương, Đào Thị Bích Thủy, Trần Đình Hà (2005), Y học Hạt Nhân, ĐH Y Hà Nội [2] Nguyễn Xuân Cử, Bùi Diệu (2011), Cơ sở Vật Lý tiến kỹ thuật xạ trị, NXB Y học, Hà Nội [3] Hà Văn Hải (2010), Luận văn: Xác định vài thông số đặc trưng chùm electron lượng 6MeV, 9MeV, 15MeV phát từ máy gia tốc PRIMUS dùng xạ trị, TP Hồ Chí Minh [4] Phan Thị Thành Lý (2009), Luận văn: Tính tốn liều gây chùm tia photon từ máy gia tốc dùng Monte Carlo code EGSnrc, TP Hồ Chí Minh [5] Đặng Nguyên Phương (2012), Hướng dẫn sử dụng MCNP cho hệ điều hành Window, TP Hồ chí Minh [6] Lê Thanh Xuân (2010), Luận văn: Mơ máy gai tốc tuyến tính dùng xạ trị phương pháp Monte Carlo, TP Hồ Chí Minh Tiếng Anh [7] Ann Barrett, Jane Dobbs, Stephen Morris, Tom Roques (2008), Paractical Radiotherapy Planning, pp.48-257 [8] Micheal Goitein (2008), Radiation Oncology: A Physicssts – Eye View, pp.111114 [9] Prowess (2003), Prowess Panther User Manual Trang mạng [10] http://en.wikipedia.org/wiki/Cancer 44 [11] http://ungthu.net.vn/cac-phuong-phap-dieu-tri-ung-thu [12] http://benhvienungbuouhungviet.vn/tin-tuc/who-canh-bao-bung-phat-dich-ungthu.aspx [13] http://www.khoahoc.com.vn/doisong/yhoc/suc-khoe/22269_gen-bao-ve-co-thekhoi-ung-thu-phoi.aspx [14] http://medlatec.vn/chi-tiet/can-lam-sang/-phat-hien-som-ung-thu-da-day 224350.aspx [15] http://dieutriungthu.net/dieu-tri-ung-thu/xa-tri-bao-ton-trong-cac-truong-hopung-thu-vu-giai-doan-som-bang-may-xa-tri-trong-suat-lieu-cao.aspx [16] http://www.bvtwhue.com.vn/index.asp?folder=KHOALAMSAN&lang=vn&q= 32 [17] http://ungthubachmai.com.vn/kin-thc-y-khoa/item/1328-x%E1%BA%A1tr%E1%BB%8B-c%C3%A1c-ph%C6%B0%C6%A1ng-ph%C3%A1px%E1%BA%A1-tr%E1%BB%8B.html 45 ... pháp để điều trị ung thư 1.2 Phương pháp xạ trị Xạ trị chia làm hai phương pháp chính: xạ trị xạ trị ngồi 1.2.1 Xạ trị (Brachytherapy) Xạ trị hay gọi xạ trị áp sát kỹ thuật điều trị sử dụng nguồn... TỒNG QUAN VỀ UNG THƯ VÀ PHƯƠNG PHÁP XẠ TRỊ 1.1 Tổng quan ung thư xạ trị 1.1.1 Ung thư 1.1.2 Các phương pháp chẩn đoán điều trị ung thư 1.2 Phương pháp xạ trị ... quan ung thư phương pháp xạ trị Chương 2: Giới thiệu máy gia tốc tuyến tính chương trình MCNP5 Chương 3: Lập kế hoạch xạ trị ung thư phổi dựa phần mềm Prowess Panther Chương 4: Mơ máy gia tốc tuyến