Luận văn thạc sĩ áp dụng EGSnrc trong việc tính phân bổ liều quanh nguồn phóng xạ dùng trong xạ trị áp sát

THƯ VIỆN BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM TP HỒ CHÍ MINH TRẦN SỸ HUY ÁP DỤNG EGSnrc TRONG VIỆC TÍNH PHÂN BỐ LIỀU QUANH NGUỒN PHÓNG XẠ DÙNG TRONG XẠ TRỊ ÁP SÁT Chuyên ngành: VẬT LÝ NGUYÊN TỬ, HẠT NHÂN VÀ NĂNG LƯỢNG CAO Mã số: 60 44 05 LUẬN VĂN THẠC SĨ VẬT LÝ NGƯỜI HƯỚNG DẪN KHOA HỌC: TS NGUYỄN ĐÔNG SƠN Thành phố Hồ Chí Minh - 2010 DANH MỤC CÁC KÝ HIỆU VÀ CÁC CHỮ VIẾT TẮT Các ký hiệu K a (d ref ) a : Suất kerma khơng khí đo khơng khí khoảng cách 1m tính từ nguồn Sk (air kerma strength): cường độ kerma khơng khí để đặc trưng cho độ mạnh nguồn ΓX: số suất liều chiếu nguồn không bị lọc D : suất liều hấp thụ K (d ) : suất kerma khơng khí X (d ref ) : Suất liều chiếu f med : hệ số chuyển đổi    x : số suất liều đồng vị nguồn an : số dị hướng Λ: số suất liều μ: hệ số suy giảm tuyến tính A (activity): Hoạt độ Aapp (apparent activity): Hoạt độ biểu kiến D: Liều hấp thụ F(r,θ): hàm dị hướng g(r): hàm liều xuyên tâm G(r,θ): hệ số hình học T(r) : hệ số hấp thụ mơ Các chữ viết tắt CT: Computed Tomography EGS: Electron Gamma Shower HDR: High dose rate HVL: Half Value Layer IAEA: International Atomic Energy Agency ICRU: International Commission on Radiation Units and Measurement LDR: Low dose rate MCNP: Monte Carlo Neutron-Photon MRI: Magnetic resonance imaging SD: Standard deviation TG-AAPM: Task Group- The American Association of Physicists in Medicine WHO: World Health Organization XTAS: Xạ trị áp sát MỞ ĐẦU Xạ trị kỹ thuật áp dụng xạ ion hóa điều trị bệnh ung thư Cùng với phẫu thuật hóa trị, xạ trị phương pháp điều trị ung thư quan trọng hiệu Với phát triển khoa học kỹ thuật với đầu tư trang thiết bị đại, ngành xạ trị với ưu điểm riêng áp dụng rộng rãi việc điều trị bệnh ung thư Xạ trị áp sát kỹ thuật điều trị dùng nguồn đồng vị phóng xạ áp sát vị trí cần điều trị Hiệu điều trị phụ thuộc vào việc tính tốn xác phân bố liều thể bệnh nhân Kết tính tốn phân bố liều phụ thuộc phần quan trọng vào việc đánh giá xác độ mạnh nguồn phát xạ Theo khuyến cáo IAEA, nguồn phải đo đạc trước sử dụng để đánh giá độ mạnh, dù có số liệu cung cấp nhà sản xuất [13] Việc đo đạc nguồn không đơn giản, đặc biệt điều kiện thiếu thốn trang thiết bị thích hợp Việt Nam Bên cạnh phép đo, độ mạnh nguồn ước lượng dựa tính tốn Monte Carlo Đề tài luận văn nhằm mục đích tìm hiểu khả Trong số chương trình Monte Carlo sử dụng nay, PENELOPE, MCNP, GEANT4, EGSnrc,… EGSnrc áp dụng cho vùng lượng phù hợp với việc tính liều thừa nhận rộng rãi tiêu chuẩn để tính liều xạ trị Đối với Việt Nam chương trình mẻ Chính mục đích luận văn tìm hiểu code EGSnrc với tinh thần học hỏi cách sử dụng để áp dụng tính liều xạ trị áp sát Nội dung chủ yếu luận văn tìm hiểu áp dụng EGSnrc việc tính phân bố liều số nguồn sử dụng xạ trị áp sát MỤC LỤC TỔNG QUAN .7 CHƯƠNG 1: XẠ TRỊ ÁP SÁT VÀ PHƯƠNG PHÁP TÍNH LIỀU 10 1.1 Giới thiệu tổng quan kỹ thuật xạ trị áp sát 10 1.1.1.Tổng quan 10 1.1.2 Tình hình phát triển xạ trị áp sát Việt Nam 10 1.2 Đặc trưng nguồn xạ 12 1.2.1 Các đại lượng đặc trưng cho độ mạnh nguồn 12 1.2.2 Yêu cầu chung nguồn xạ 13 1.2.3 Nguồn Ir-192 14 1.2.4 Nguồn Cs-137 15 1.3 Phương pháp tính suất liều nguồn xạ dùng xạ trị áp sát 17 1.3.1 Liều hấp thụ D suất liều hấp thụ D 17 1.3.2 Cơng thức tính suất liều hấp thụ 17 1.3.3 Yêu cầu độ xác việc cấp liều cho bệnh nhân 22 CHƯƠNG 2: PHƯƠNG PHÁP MONTE CARLO TRONG VIỆC TÍNH LIỀU CHƯƠNG TRÌNH EGSnrc 24 2.1 Phương pháp mô Monte Carlo 24 2.1.1 Giới thiệu tổng quan phương pháp Monte Carlo 24 2.1.2 Sự tạo số ngẫu nhiên 25 2.1.3 Quá trình tương tác photon 26 2.2 Mô Monte Carlo vận chuyển Photon .27 2.2.1 Mô phóng xạ sơ cấp [4], [18] .27 2.2.2 Chọn loại tương tác [4] 28 2.3 Giới thiệu chương trình EGSnrc 29 2.3.1 Giới thiệu chung 29 2.3.2 DOSRZnrc 31 CHƯƠNG 3: ÁP DỤNG CHƯƠNG TRÌNH EGSnrc CHO CÁC NGUỒN Ir-192 VÀ Cs-137 DÙNG TRONG XẠ TRỊ ÁP SÁT 38 3.1 Thiết lập thông số cấu trúc hình học cho nguồn 38 3.1.1 Thành phần cấu tạo cấu trúc nguồn .38 3.1.2 Khai báo thơng số cho chương trình EGSnrc 40 3.2 Tính giá trị g(r) DF(r, θ) nguồn 46 3.2.1 Kết hàm g(r) .46 3.2.2 Kết hàm DF(r, θ) 51 3.2.3 Kết luận 55 KẾT LUẬN .57 HƯỚNG PHÁT TRIỂN 58 PHỤ LỤC 59 TÀI LIỆU THAM KHẢO .69 TỔNG QUAN Cho tới nay, bệnh ung thư trở thành nguyên nhân đứng thứ hai gây tử vong nhóm bệnh khơng lây nhiễm Tại Việt Nam, theo số liệu thống kê tình trạng bệnh ung thư bệnh viện K, bệnh viện Ung Bướu TP Hồ Chí Minh hai trung tâm chuẩn đoán điều trị ung thư xạ số tỉnh thành, ước tính năm nước ta có khoảng 150.000 trường hợp mắc 75.000 người chết ung thư số có xu hướng ngày gia tăng Dự kiến ngành y tế, đến năm 2020 năm Việt Nam có khoảng 200.000 trường hợp mắc 100.000 trường hợp chết ung thư [24] Tuy nhiên, theo chuyên gia, ung thư bệnh vô phương cứu chữa, phát sớm điều trị phác đồ, bệnh nhân ung thư hoàn toàn chữa khỏi tăng thêm thời gian chất lượng sống cho người bệnh Hiện có phương pháp để điều trị ung thư là: phẫu thuật, hoá trị xạ trị Việc lựa chọn phương pháp điều trị phụ thuộc vào nhiều yếu tố điều kiện điều trị bệnh viện, vị trí khối u, giai đoạn bệnh tình trạng bệnh nhân Trong năm gần đây, xạ trị áp sát phát triển mạnh mẽ phải kể đến đời máy điều trị xạ trị áp sát suất liều cao (high dose rate, HDR) sử dụng nước phát triển giới Ở nước ta, Bệnh viện Ung Bướu TP.HCM đơn vị đưa vào hoạt động máy HDR cho việc điều trị ung thư cổ tử cung (từ năm 2000) Cho đến nay, hàng ngàn ca bệnh điều trị kết thu qua báo cáo tổng kết đáng phấn khởi [25] Kể từ đời vào đầu kỷ 20, xạ trị áp sát tiến hóa có mối liên kết chặt chẽ với vật lý y học Lịch sử 50 năm Hiệp hội Vật lý Y học Mỹ (AAPM) không với xuất vật lý y học nghề trưởng thành, với đổi cách mạng thật vật lý xạ, bao gồm lò phản ứng hạt nhân, máy gia tốc hạt mới, hình ảnh 3D, máy tính hỗ trợ điều trị việc lập kế hoạch, với tiến việc hiểu biết trình vận chuyển xạ điều chế phản ứng lâm sàng, … làm thay đổi đáng kể việc áp dụng kỹ thuật xạ trị áp sát [1] Có lẽ tác động cao tiến công nghệ nửa cuối kỷ qua điều chế đồng vị phóng xạ nhân tạo hệ thống nạp nguồn sau vào xạ trị áp sát Tiến làm giảm đáng kể chi phí, giảm tiếp xúc, tăng tính linh hoạt kỹ thuật [2] Xạ trị áp sát phương thức điều trị nguồn phát xạ (nguồn đồng vị phóng xạ đóng gói) đặt áp sát hay bên khối u Ưu điểm cung cấp liều tập trung vào khối u ảnh hưởng đến mơ lành, nhờ quy luật giảm theo bình phương khoảng cách Tuy nhiên, liều cung cấp cho bệnh nhân lớn, sai sót việc cấp liều cho bệnh nhân dẫn đến nguy hiểm cho mô lành, việc xác định xác liều xạ trị áp sát việc vô quan trọng [4] Phương pháp tính liều hành xạ trị áp sát dựa hình thức TG AAPM-43, dựa vào chồng chập nguồn đơn thu phantom nước với thể tích thích hợp cho tán xạ xạ [20] Cách tiếp cận thơng qua hệ thống điều trị máy tính thay cho phương pháp cổ điển hệ thống Manchester Paris [19] Một bước quan trọng q trình xạ trị tính phân bố liều thể bệnh nhân, hiệu điều trị phụ thuộc nhiều vào khâu Có nhiều phương pháp để tính liều phương pháp Monte Carlo cho phép tính phân bố liều với độ xác cao, hạn chế thời gian tính toán lâu Trong nửa đầu kỉ 20 đến năm 1960, việc tính phân bố liều xạ trị áp sát lâm sàng sử dụng cách tra bảng dựa việc đo buồng ion hóa tích phân Sievert Vào năm 1971, dựa mơ hình phân bố liều nguồn xạ trị áp sát đối xứng trụ môi trường chất lỏng hay rắn, Krishnaswamy thiết lập móng cho việc mơ việc tính liều xạ trị áp sát dùng phương pháp Monte Carlo [1] Qua thập kỷ áp dụng kỹ thuật mô Monte Carlo, việc mô phân bố liều quanh nguồn xạ trị áp sát ngày đạt mức độ chi tiết hơn, chẳng hạn việc tính hệ số chồng chập (build-up) nguồn điểm đẳng hướng cho phân bố liều xạ trị áp sát Berger, webb Fox, báo tiếng Meisberger [17] Nhiều cách tiếp cận tinh vi 3D theo đuổi sau đó, đầu Williamson [1] Với việc tăng tốc độ xử lý máy tính, khả tiếp cận hoạt động hệ thống sử dụng rộng rãi qua nghiên cứu y tế, tính phân bố liều phương pháp mô Monte Carlo đạt đến ngưỡng quan trọng năm 1990, có báo cáo quan trọng AAPM TG-43 năm 1995 [20] Báo cáo thiết lập hình thức luận tiêu chuẩn việc tính liều quanh nguồn xạ trị áp sát, đưa tham số cho dạng nguồn khác đồng vị phóng xạ, thúc đẩy qn sử dụng thơng số tính liều tổ chức riêng biệt Thơng qua đó, nâng cao khả sử dụng rộng rãi phương pháp Monte Carlo việc tính phân bố liều nguồn dùng xạ trị áp sát Phương pháp Monte Carlo mở rộng mã PTRAN Williamson loại nguồn khía cạnh khác nguồn kiểm tra Kể từ sau Krishnaswamy, kỹ thuật mơ Monte Carlo trở thành phương pháp việc tính tốn phân bố liều cho nguồn xạ trị áp sát Tuy nhiên, việc mô cần phải cải tiến để ngày xác hơn, kết cần phải luôn so sánh với thực nghiệm Kết mơ Monte Carlo khơng xác sai khác nguồn mơ tả nguồn thực tế, đơn giản hóa mơ tả phổ xạ, v.v Do đó, việc so sánh kết Monte Carlo kết thực nghiệm cần thiết để phát khác biệt lý thuyết thực tế hệ thống trình đo [1] Một chương trình Monte Carlo sử dụng phổ biến nghiên cứu xạ trị code EGSnrc Mục đích luận văn tìm hiểu chương trình EGSnrc, cụ thể code DOSRZnrc áp dụng để tính phân bố liều quanh nguồn sử dụng xạ trị áp sát, so sánh với thực nghiệm rút kết luận cách sử dụng code phân bố liều qua lớp vỏ nguồn khác Quá trình thực luận văn việc tìm hiểu tổng quan kiến thức liên quan đến xạ trị kỹ thuật tính liều, giới thiệu phương pháp Monte Carlo xạ trị; tìm hiểu đặc điểm, thành phần, sở vật lý code EGSnrc, nhiệm vụ code cách thức sử dụng code DOSRZnrc; sau áp dụng để khảo sát tốn cụ thể mơ tả chi tiết cơng trình trước phân bố liều quanh nguồn Ir-192 Cs-137 [3], xem xét phù hợp kết thu với kết có từ cơng trình, chứng minh tính hữu dụng chương trình để rút học cách sử dụng code Từ mục đích nội dung cơng việc trên, luận văn có bố cục sau: Chương 1- Lý thuyết sở Chương trình bày vấn đề liên quan đến xạ trị áp sát, tình hình xạ trị áp sát Việt Nam, đại lượng đặc trưng cho độ mạnh nguồn, nêu lên yêu cầu độ xác việc tính liều hấp thụ thể bệnh nhân Chương giới thiệu công thức tính liều hấp thụ dựa hình thức luận AAPM-TG43 sử dụng phổ biến số nguồn phóng xạ dùng xạ trị áp sát, đặc biệt hai nguồn Ir-192 Cs-137 mà khảo sát chương sau Chương 2- Lý thuyết chuyên sâu Chương trình bày vận chuyển photon vật chất theo phương pháp mô Monte Carlo, phương pháp Monte Carlo xạ trị áp sát, chương trình EGSnrc code DOSRZnrc chương trình Chương đóng vai trị quan trọng nội dung luận văn Các kiến thức trình bày chương cần thiết để áp dụng EGSnrc DOSRZnrc, mô tả chương Chương 3- Áp dụng chương trình EGSnrc để khảo sát phân bố liều quanh nguồn Ir-192 Cs-137 Chương trình bày nội dung công việc bước việc sử dụng mã DOSRZnrc chương trình EGSnrc để tính liều hấp thụ từ rút kết luận phân bố liều quanh nguồn Ir-192 Cs-137 So sánh với thực nghiệm cơng trình công bố C.Thomason, T.R.Mackie, M.J.Lindstrom P.D.Higgins [3] nhằm chứng minh tính hữu dụng chương trình CHƯƠNG 1: XẠ TRỊ ÁP SÁT VÀ PHƯƠNG PHÁP TÍNH LIỀU Chương trình bày khái niệm kỹ thuật xạ trị áp sát, tình hình phát triển xạ trị áp sát Việt Nam, ý việc sử dụng xạ trị áp sát nguồn thường sử dụng xạ trị áp sát Một số đại lượng công thức quan trọng có liên quan đến việc tính liều cần cho phần sau, giới thiệu 1.1 Giới thiệu tổng quan kỹ thuật xạ trị áp sát 1.1.1.Tổng quan Xạ trị áp sát (XTAS) phương thức điều trị nguồn phát xạ (nguồn đồng vị phóng xạ đóng gói) đặt áp sát hay bên khối u XTAS áp dụng cho nhiều trường hợp ung thư, thành công cho phụ khoa ung thư đầu cổ Ban đầu XTAS phát triển để điều trị khối u nằm sâu mà kỹ thuật xạ trị thời kỳ đầu không mang lại hiệu XTAS hốc thích hợp điều trị phụ khoa, đưa nguồn vào qua âm đạo XTAS kẽ thích hợp cho ung thư đầu cổ dễ tiếp cận qua đường miệng mũi Kỹ thuật cấy nguồn vĩnh viễn để điều trị ung thư tiền liệt tuyến ngày phổ biến [21] XTAS áp dụng độc lập (ung thư tuyến tiền liệt ung thư vú giai đoạn đầu) hay kết hợp với xạ trị (ung thư phụ khoa, ung thư tuyến tiền liệt giai đoạn trễ, ung thư đầu cổ) Cũng áp dụng sau phẫu thuật để diệt phần cịn sót lại mơ ung thư [21] Trong khoa xạ trị, thường khoảng từ 10% đến 20% bệnh nhân điều trị xạ trị áp sát [6] XTAS phát triển mạnh Châu Âu (Paris, Manchester, Stockholm) Trong nửa đầu kỷ 20, nguồn xạ dùng radium Nguồn Ra-226 phát photon lượng cao, chịu tương tác quang điện xương, thích hợp để điều trị mô ung thư nằm gần xương mà không sợ bị hoại tử xương Một ưu điểm khác nguồn Ra-226 có chu kỳ bán rã lớn, nên khơng phải hiệu chỉnh thay nguồn thời gian sử dụng Ưu sau khơng cịn quan trọng xạ trị ngồi có lượng cao phát triển năm 1950-1960 Radium sinh khí radon phóng xạ, làm chúng bị thay nguồn cesium Cs-137, có hạt nhân chất rắn [21 ] Ưu điểm: so với xạ trị từ xa XTAS cung cấp liều tập trung vào khối u ảnh hưởng đến mô lành, nhờ qui luật giảm bình phương theo khoảng cách cường độ Nhược điểm: Chỉ áp dụng cho khối u tập trung nhỏ, ngồi cịn cần phải can thiệp vào thể bệnh nhân, cần nhiều thời gian công sức trình điều trị [6] 1.1.2 Tình hình phát triển xạ trị áp sát Việt Nam Hiện nay, với phát triển công nghệ cơ-điện tử máy tính, xạ trị áp sát đạt bước đột phá phương tiện kỹ thuật điều trị Các máy xạ trị suất liều cao (High Dose Rate - HDR) sử dụng rộng rãi hầu hết sở xạ trị giới Quy trình xạ trị hiệu ứng định luật bình phương tối thiểu Sự giảm tính dị hướng khoảng 3cm tăng tán xạ liều, tán xạ ngày tăng bù lại cho biến đổi liều ban đầu gây suy giảm khác nguồn lớp vỏ Tính dị hướng thể rõ nguồn Ir với lớp vỏ platinum Sự khác biệt tính dị hướng nguồn Ir phù hợp với khác biệt tính chất vật lý lớp vỏ chúng lõi nguồn có lớp vỏ platinum có khối lượng lớn nguồn có lớp vỏ stainless steel Vì platinum có khối lượng nguyên tử mật độ lớn stainless steel suy giảm nhanh hơn, phân bố liều từ nguồn với lớp vỏ platinum mong đợi có tính dị hướng lớn có suy giảm lớn Sự giảm tính dị hướng từ 3cm đến 12,5cm tính từ nguồn có lớp vỏ platinum phù hợp với suy giảm lớn hơn, có tán xạ lớn platinum g(r) thu phù hợp với DF(r, 90o) Từ nhận xét rút số kết luận sau: Ảnh hưởng lớp vỏ nguồn phân bố liều dọc theo đường trung trực nguồn nhỏ Phân bố liều xung quanh nguồn Ir với lớp vỏ platinum tính dị hướng rõ ràng nguồn Ir với lớp vỏ stainless steel Xác minh tính hữu dụng mã Monte Carlo EGSnrc việc mơ hình nguồn Ir Cs để thu phân bố liều Điều có vai trị quan trọng đặc biệt lĩnh vực xạ trị áp sát so với kỹ thuật tính tốn khác chúng khó dễ mắc nhiều lỗi KẾT LUẬN Việc tính liều xạ trị 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 Những nghiên cứu gần cho thấy phương pháp Monte Carlo phương tiện hàng đầu cho việc tính liều xạ trị Tuy nhiên hạn chế thời gian tính tốn nên code chưa đưa vào tính liều xạ trị cách trọn vẹn mà có kết hợp với phương pháp giải tích Trên giới, code Monte Carlo kiểm tra tính xác khơng ngừng cải tiến, nâng cao tốc độ tính tốn để bước áp dụng vào tính liều xạ trị Thành cơng tốn mơ phụ thuộc vào nhiều yếu tố Trong đó, người thực khơng tìm hiểu nguyên tắc hoạt động cách sử dụng code mà phải hiểu biết thật chi tiết đặc trưng chùm xạ chiếu đến, tức hiểu rõ thông tin cấu tạo nguồn xạ phát Với thành cơng từ tốn mơ phỏng, tin tưởng EGS code mạnh việc mô vận chuyển xạ Kỹ thuật Monte Carlo tiếp tục phát triển mạnh vật lý y khoa tốc độ máy tính ngày cao dự án code Monte Carlo tiếp tục phát triển Chúng tơi hồn tồn tin tưởng rằng, kỹ thuật Monte Carlo sử dụng tính liều xạ trị để điều trị bệnh ung thư tương lai không xa HƯỚNG PHÁT TRIỂN Ngồi việc phải nắm rõ vai trị cách sử dụng code chương trình EGSnrc độ xác kết cịn phụ thuộc nhiều vào số lịch sử N Việc cho số lịch sử lớn cho kết với độ xác cao lại địi hỏi thời gian tính tốn lâu Để khắc phục tình trạng người ta đưa thuật toán làm giảm thăng giáng đồng thời đưa chương trình hỗ trợ việc chạy song song máy tính (một chương trình chạy nhiều máy tính kết máy tính tổng hợp lại xem kết việc chạy máy nhất) Ngoài ra, phát triển cơng nghệ máy tính mạnh, việc tạo CPU mạnh có cơng suất lớn góp phần làm thời gian tính tốn giảm xuống đáng kể Trong luận văn cho chạy chương trình EGSnrc Window nên khơng thể tận dụng việc chạy song song máy tính (việc chạy song song hóa thực Linux) Ngồi ra, để tăng độ xác việc xác định liều hấp thụ việc khai báo nguồn xạ xác giống thật tốt Trong luận văn này, mô cách gần cho môi trường chứa Stainless Steel balls nguồn Cs, chúng có dạng hình cầu (mơ gần chúng có dạng hình trụ đồng chất) Thực tế có cách tiếp cận việc mơ cấu trúc hình học nguồn có hình dạng phức tạp [7] Việc làm giảm thời gian tính tốn cách sử dụng thuật toán làm giảm thăng giáng việc chạy song song máy tính việc khai báo cấu trúc hình học cách xác công việc cần tiếp tục thực tương lai PHỤ LỤC Bảng A1 Giá trị DF(r, θ) nguồn Ir bao plantinum, sử dụng phổ Ir192bare.spectrum theo cơng thức báo Góc (độ) 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 0,682 0,703 0,778 0,831 0,867 0,914 0,933 0,953 0,964 0,977 0,986 0,996 0,991 0,998 1,002 0,999 1,006 1,001 1,000 1,5 0,684 0,716 0,760 0,823 0,860 0,904 0,926 0,951 0,969 0,968 0,988 1,000 1,003 0,998 1,011 1,037 1,009 1,010 1,009 0,685 0,703 0,773 0,825 0,873 0,921 0,935 0,946 0,970 0,973 0,983 0,996 1,005 1,006 1,007 1,010 1,012 1,016 1,015 2,5 0,709 0,699 0,760 0,820 0,867 0,921 0,938 0,953 0,966 0,985 0,989 0,998 1,002 1,011 1,009 1,007 1,013 1,018 1,015 0,714 0,731 0,791 0,844 0,878 0,919 0,935 0,955 0,968 0,982 0,991 0,997 0,999 1,008 1,020 1,017 1,011 1,019 1,007 0,740 0,749 0,797 0,858 0,880 0,917 0,940 0,957 0,975 0,979 0,995 1,005 1,003 1,008 1,020 1,006 1,008 1,019 1,017 Khoảng cách (cm) 0,734 0,764 0,723 0,751 0,742 0,757 0,803 0,813 0,820 0,867 0,873 0,871 0,887 0,904 0,892 0,914 0,923 0,934 0,940 0,933 0,937 0,967 0,954 0,935 0,972 0,956 0,961 0,984 0,965 0,961 0,988 0,984 0,985 0,996 0,986 0,983 0,992 1,001 0,982 1,001 0,998 0,988 1,005 1,009 0,992 1,006 0,993 0,996 1,004 1,003 0,999 1,012 1,005 0,992 1,016 1,012 1,006 0,738 0,775 0,801 0,853 0,880 0,894 0,927 0,938 0,943 0,963 0,964 0,967 0,978 0,978 0,974 0,979 0,976 0,987 0,969 0,678 0,769 0,823 0,851 0,867 0,883 0,913 0,915 0,935 0,937 0,946 0,955 0,956 0,957 0,960 0,955 0,972 0,969 0,960 10 0,819 0,764 0,801 0,832 0,845 0,856 0,900 0,898 0,928 0,926 0,936 0,942 0,938 0,937 0,952 0,948 0,943 0,955 0,961 11 0,754 0,775 0,790 0,830 0,833 0,866 0,883 0,890 0,894 0,883 0,907 0,917 0,926 0,918 0,919 0,927 0,942 0,926 0,933 12 0,756 0,755 0,758 0,831 0,827 0,860 0,835 0,859 0,873 0,886 0,895 0,893 0,905 0,910 0,909 0,895 0,902 0,908 0,907 12,5 0,778 0,746 0,772 0,792 0,796 0,830 0,841 0,870 0,867 0,872 0,884 0,872 0,888 0,905 0,901 0,900 0,902 0,895 0,889 Bảng A2 Giá trị DF(r, θ) nguồn Ir bao plantinum, sử dụng phổ Ir192bare.spectrum theo AAPM TG-43 Góc (độ) 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 0,662 0,682 0,755 0,808 0,845 0,892 0,913 0,934 0,947 0,962 0,975 0,986 0,984 0,993 0,999 0,997 1,005 1,001 1,000 1,5 0,669 0,700 0,743 0,806 0,842 0,886 0,908 0,934 0,952 0,953 0,974 0,988 0,991 0,988 1,001 1,027 1,000 1,001 1,000 0,670 0,688 0,755 0,807 0,853 0,901 0,915 0,926 0,951 0,954 0,965 0,979 0,987 0,989 0,990 0,995 0,996 1,001 1,000 2,5 0,695 0,684 0,745 0,804 0,850 0,904 0,920 0,936 0,949 0,969 0,972 0,981 0,986 0,995 0,993 0,992 0,998 1,003 1,000 0,707 0,724 0,783 0,836 0,869 0,910 0,927 0,947 0,959 0,973 0,984 0,989 0,992 1,001 1,013 1,010 1,004 1,012 1,000 0,726 0,735 0,782 0,842 0,863 0,900 0,923 0,940 0,958 0,962 0,977 0,987 0,986 0,991 1,003 0,989 0,991 1,002 1,000 Khoảng cách (cm) 0,721 0,755 0,718 0,739 0,733 0,752 0,789 0,803 0,815 0,853 0,862 0,865 0,872 0,893 0,886 0,899 0,912 0,927 0,925 0,921 0,931 0,950 0,942 0,929 0,955 0,945 0,956 0,967 0,954 0,955 0,972 0,972 0,980 0,980 0,975 0,977 0,976 0,990 0,976 0,985 0,987 0,982 0,989 0,998 0,987 0,990 0,982 0,990 0,988 0,991 0,993 0,996 0,994 0,986 1,000 1,000 1,000 0,761 0,799 0,826 0,880 0,908 0,923 0,957 0,968 0,973 0,994 0,994 0,998 1,009 1,009 1,005 1,010 1,007 1,018 1,000 0,706 0,801 0,857 0,886 0,902 0,919 0,951 0,953 0,973 0,976 0,985 0,995 0,995 0,997 1,000 0,995 1,012 1,009 1,000 10 0,852 0,795 0,834 0,866 0,879 0,891 0,937 0,935 0,966 0,964 0,975 0,980 0,976 0,975 0,991 0,987 0,981 0,994 1,000 11 0,808 0,831 0,847 0,889 0,892 0,928 0,947 0,954 0,958 0,947 0,972 0,983 0,993 0,984 0,985 0,993 1,009 0,993 1,000 12 0,833 0,832 0,836 0,916 0,911 0,947 0,920 0,946 0,962 0,977 0,987 0,984 0,998 1,003 1,002 0,987 0,994 1,000 1,000 12,5 0,875 0,839 0,868 0,890 0,895 0,934 0,945 0,978 0,975 0,980 0,994 0,980 0,998 1,017 1,013 1,011 1,014 1,006 1,000 Bảng A3 Giá trị DF(r, θ) nguồn Ir bao plantinum, sử dụng phổ Ir_alpha_omega.spectrum theo công thức báo Góc (độ) 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 0,698 0,699 0,776 0,823 0,861 0,908 0,932 0,954 0,971 0,974 0,977 0,993 0,997 0,991 0,998 0,998 0,993 0,997 1,000 1,5 0,681 0,685 0,762 0,818 0,862 0,904 0,926 0,947 0,965 0,897 0,985 0,993 0,991 1,001 1,003 1,035 1,004 1,005 1,003 0,658 0,696 0,778 0,832 0,867 0,911 0,936 0,945 0,966 0,981 0,986 0,990 0,997 0,999 0,998 1,011 1,007 1,013 1,010 2,5 0,718 0,711 0,775 0,826 0,877 0,910 0,932 0,951 0,971 0,979 0,983 0,991 0,998 0,999 1,013 1,009 1,009 1,009 1,017 0,704 0,723 0,784 0,829 0,876 0,907 0,933 0,955 0,965 0,988 0,993 0,991 1,004 1,002 1,010 1,009 1,014 1,005 1,014 0,715 0,730 0,777 0,852 0,882 0,914 0,950 0,955 0,974 0,984 0,994 0,996 1,003 1,005 1,016 1,006 1,011 1,004 1,005 Khoảng cách (cm) 0,795 0,789 0,808 0,740 0,757 0,777 0,803 0,785 0,806 0,851 0,864 0,845 0,877 0,876 0,873 0,925 0,915 0,910 0,936 0,928 0,936 0,959 0,959 0,953 0,975 0,959 0,956 0,978 0,968 0,966 0,996 0,978 0,960 0,994 0,984 0,971 0,996 0,986 0,973 0,999 0,998 0,992 1,009 0,996 0,986 1,014 0,996 0,985 1,000 1,000 0,993 1,006 1,003 0,991 1,005 1,003 0,997 0,769 0,765 0,808 0,849 0,870 0,903 0,927 0,929 0,941 0,947 0,961 0,978 0,968 0,973 0,983 0,973 0,978 0,969 0,980 0,733 0,758 0,810 0,841 0,866 0,895 0,903 0,912 0,916 0,921 0,940 0,942 0,955 0,955 0,963 0,954 0,967 0,966 0,974 10 0,807 0,771 0,796 0,836 0,864 0,881 0,889 0,906 0,908 0,915 0,919 0,918 0,938 0,940 0,953 0,946 0,929 0,944 0,943 11 0,695 0,753 0,782 0,816 0,839 0,856 0,874 0,878 0,894 0,907 0,900 0,907 0,917 0,922 0,926 0,929 0,921 0,934 0,932 12 0,770 0,769 0,772 0,802 0,824 0,845 0,849 0,865 0,885 0,868 0,882 0,887 0,899 0,900 0,903 0,901 0,903 0,902 0,895 12,5 0,701 0,753 0,765 0,784 0,809 0,821 0,846 0,863 0,869 0,876 0,870 0,879 0,896 0,880 0,894 0,888 0,896 0,884 0,894 Bảng A4 Giá trị DF(r, θ) nguồn Ir bao plantinum, sử dụng phổ Ir_alpha_omega.spectrum theo AAPM TG-43 Góc (độ) 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 0,677 0,678 0,753 0,800 0,838 0,886 0,912 0,935 0,954 0,959 0,966 0,983 0,990 0,986 0,995 0,996 0,992 0,997 1,000 1,5 0,670 0,674 0,749 0,806 0,849 0,891 0,914 0,936 0,954 0,888 0,977 0,986 0,985 0,995 0,998 1,031 1,001 1,002 1,000 0,646 0,684 0,765 0,817 0,852 0,896 0,921 0,930 0,952 0,967 0,973 0,978 0,984 0,988 0,987 1,001 0,997 1,002 1,000 2,5 0,702 0,696 0,758 0,809 0,859 0,891 0,913 0,932 0,952 0,961 0,965 0,972 0,980 0,981 0,995 0,992 0,992 0,992 1,000 0,692 0,711 0,770 0,814 0,861 0,892 0,918 0,939 0,950 0,972 0,978 0,976 0,989 0,987 0,996 0,995 1,000 0,991 1,000 0,711 0,725 0,771 0,846 0,876 0,908 0,945 0,950 0,968 0,978 0,988 0,990 0,998 1,000 1,011 1,001 1,007 1,000 1,000 Khoảng cách (cm) 0,790 0,786 0,810 0,736 0,754 0,779 0,798 0,782 0,808 0,846 0,861 0,847 0,871 0,873 0,875 0,919 0,912 0,912 0,930 0,925 0,939 0,953 0,956 0,956 0,969 0,957 0,959 0,972 0,965 0,969 0,991 0,976 0,963 0,989 0,982 0,974 0,991 0,983 0,976 0,994 0,996 0,995 1,004 0,994 0,989 1,009 0,993 0,988 0,995 0,997 0,996 1,001 1,000 0,995 1,000 1,000 1,000 0,785 0,781 0,825 0,866 0,888 0,922 0,946 0,949 0,961 0,967 0,981 0,998 0,988 0,993 1,003 0,994 0,999 0,989 1,000 0,752 0,778 0,831 0,864 0,889 0,919 0,927 0,936 0,941 0,946 0,965 0,967 0,981 0,980 0,989 0,979 0,993 0,992 1,000 10 0,856 0,818 0,844 0,886 0,916 0,934 0,942 0,961 0,963 0,970 0,974 0,973 0,995 0,997 1,011 1,004 0,986 1,001 1,000 11 0,746 0,808 0,839 0,875 0,900 0,919 0,938 0,942 0,959 0,973 0,966 0,973 0,984 0,989 0,994 0,997 0,988 1,002 1,000 12 0,860 0,859 0,862 0,896 0,920 0,944 0,948 0,967 0,989 0,970 0,985 0,991 1,004 1,006 1,008 1,006 1,009 1,008 1,000 12,5 0,784 0,843 0,856 0,878 0,905 0,918 0,947 0,966 0,972 0,980 0,974 0,983 1,002 0,985 1,000 0,994 1,002 0,989 1,000 Bảng A5 Giá trị DF(r, θ) nguồn Ir bao stainless steel, sử dụng phổ Ir192bare.spectrum theo cơng thức báo Góc (độ) 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 0,776 0,836 0,934 0,952 0,977 0,987 0,997 1,001 1,000 1,002 1,007 1,004 1,003 1,007 1,002 1,002 1,002 1,001 1,000 1,5 0,743 0,839 0,912 0,939 0,964 0,987 0,991 0,998 1,002 1,006 1,004 1,008 1,010 1,007 1,006 1,012 1,005 1,011 1,011 0,728 0,844 0,921 0,947 0,963 0,987 0,992 1,005 1,004 1,001 1,018 1,005 1,013 1,010 1,014 1,015 1,010 1,017 1,010 2,5 0,733 0,832 0,924 0,949 0,973 0,990 0,997 1,002 1,010 1,005 1,010 1,019 1,008 1,018 1,020 1,016 1,017 1,016 1,016 0,722 0,875 0,933 0,948 0,979 0,985 0,998 1,006 1,002 1,010 1,010 1,018 1,019 1,015 1,025 1,015 1,021 1,024 1,023 0,732 0,863 0,928 0,965 0,966 0,991 0,998 0,999 1,006 1,011 1,013 1,015 1,020 1,021 1,021 1,023 1,019 1,016 1,024 Khoảng cách (cm) 0,783 0,761 0,715 0,876 0,874 0,858 0,926 0,935 0,934 0,961 0,951 0,951 0,972 0,973 0,959 0,985 0,984 0,980 0,996 0,982 0,977 0,997 0,992 0,978 1,018 1,000 0,983 1,009 1,006 1,006 1,015 1,013 1,000 1,010 1,003 0,994 1,015 1,006 0,999 1,020 1,010 1,001 1,017 1,017 1,001 1,017 1,011 0,997 1,016 1,011 1,000 1,015 1,010 1,006 1,016 1,010 1,008 0,789 0,885 0,915 0,953 0,966 0,971 0,969 0,985 0,968 0,988 0,983 0,977 0,984 0,987 0,992 0,992 0,997 0,984 0,982 0,773 0,882 0,911 0,931 0,936 0,942 0,958 0,953 0,981 0,959 0,973 0,977 0,977 0,974 0,977 0,978 0,969 0,981 0,972 10 0,775 0,861 0,897 0,912 0,931 0,925 0,949 0,955 0,946 0,951 0,961 0,956 0,967 0,960 0,957 0,957 0,963 0,951 0,969 11 0,843 0,858 0,892 0,909 0,897 0,914 0,924 0,927 0,923 0,936 0,928 0,936 0,946 0,936 0,939 0,938 0,940 0,933 0,944 12 0,812 0,831 0,862 0,872 0,886 0,900 0,906 0,907 0,897 0,911 0,900 0,913 0,916 0,910 0,915 0,923 0,925 0,910 0,930 12,5 0,784 0,833 0,846 0,852 0,887 0,886 0,877 0,889 0,900 0,904 0,908 0,903 0,900 0,903 0,907 0,918 0,905 0,911 0,905 Bảng A6 Giá trị DF(r, θ) nguồn Ir bao stainless steel, sử dụng phổ Ir192bare.spectrum theo AAPM TG-43 Góc (độ) 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 0,753 0,811 0,907 0,926 0,951 0,963 0,975 0,981 0,983 0,987 0,995 0,994 0,996 1,002 0,999 1,000 1,001 1,001 1,000 1,5 0,726 0,819 0,890 0,917 0,942 0,966 0,971 0,978 0,983 0,988 0,988 0,993 0,996 0,994 0,994 1,000 0,994 1,000 1,000 0,715 0,830 0,906 0,932 0,947 0,971 0,976 0,990 0,990 0,987 1,005 0,993 1,001 0,999 1,003 1,005 1,001 1,007 1,000 2,5 0,718 0,815 0,905 0,930 0,954 0,970 0,978 0,984 0,992 0,988 0,992 1,001 0,991 1,002 1,003 1,000 1,002 1,000 1,000 0,704 0,853 0,909 0,924 0,954 0,960 0,973 0,981 0,977 0,985 0,986 0,994 0,995 0,991 1,002 0,992 0,998 1,001 1,000 0,713 0,841 0,905 0,941 0,941 0,967 0,974 0,975 0,982 0,987 0,989 0,991 0,996 0,997 0,997 1,000 0,995 0,992 1,000 Khoảng cách (cm) 0,770 0,753 0,708 0,862 0,864 0,851 0,911 0,925 0,926 0,945 0,941 0,943 0,956 0,962 0,951 0,969 0,973 0,971 0,979 0,972 0,969 0,981 0,982 0,970 1,000 0,990 0,975 0,992 0,997 0,998 0,999 1,003 0,992 0,994 0,993 0,986 0,999 0,997 0,991 1,004 1,001 0,993 1,001 1,008 0,993 1,001 1,002 0,989 1,000 1,001 0,993 0,999 1,000 0,998 1,000 1,000 1,000 0,803 0,901 0,932 0,970 0,983 0,989 0,987 1,004 0,986 1,006 1,001 0,995 1,003 1,005 1,010 1,010 1,015 1,002 1,000 0,795 0,907 0,937 0,957 0,962 0,969 0,985 0,980 1,009 0,986 1,001 1,004 1,004 1,002 1,005 1,005 0,997 1,009 1,000 10 0,800 0,889 0,926 0,941 0,961 0,955 0,980 0,986 0,976 0,982 0,992 0,987 0,998 0,991 0,988 0,988 0,994 0,982 1,000 11 0,893 0,909 0,945 0,963 0,949 0,968 0,979 0,982 0,977 0,991 0,983 0,991 1,002 0,991 0,994 0,994 0,996 0,989 1,000 12 0,873 0,894 0,926 0,938 0,953 0,967 0,974 0,975 0,964 0,979 0,967 0,982 0,985 0,979 0,984 0,992 0,994 0,978 1,000 12,5 0,866 0,919 0,934 0,941 0,980 0,978 0,968 0,982 0,994 0,999 1,003 0,997 0,994 0,998 1,002 1,014 0,999 1,006 1,000 Bảng A7 Giá trị DF(r, θ) nguồn Ir bao stainless steel, sử dụng phổ Ir192_best_Industries.spectrum theo công thức báo Góc (độ) 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 0,782 0,844 0,937 0,954 0,978 0,989 1,000 1,006 1,003 1,004 1,005 1,008 1,003 1,006 1,010 1,006 1,003 1,004 1,000 1,5 0,746 0,844 0,909 0,954 0,960 0,990 0,990 1,000 1,003 1,008 1,012 1,008 1,017 1,017 1,010 1,042 1,013 1,010 1,012 0,732 0,839 0,921 0,955 0,978 0,986 0,990 0,997 1,013 1,009 1,011 1,011 1,011 1,024 1,018 1,018 1,021 1,020 1,020 2,5 0,726 0,854 0,922 0,945 0,975 0,993 0,997 1,000 1,003 1,012 1,016 1,022 1,011 1,026 1,023 1,025 1,023 1,018 1,023 0,737 0,852 0,916 0,951 0,972 0,989 0,997 1,012 1,008 1,015 1,016 1,020 1,019 1,024 1,022 1,026 1,023 1,016 1,020 0,757 0,883 0,949 0,957 0,979 0,996 0,992 1,008 1,011 1,021 1,025 1,014 1,022 1,025 1,023 1,023 1,026 1,029 1,021 Khoảng cách (cm) 0,772 0,717 0,776 0,865 0,888 0,899 0,926 0,941 0,934 0,969 0,962 0,966 0,978 1,002 0,973 0,983 0,986 0,977 1,002 0,995 0,986 1,004 1,017 1,000 1,013 1,004 0,986 1,011 1,019 1,006 1,019 1,013 1,003 1,016 1,008 0,999 1,011 1,017 1,000 1,021 1,016 1,000 1,035 1,017 0,997 1,014 1,019 1,005 1,030 1,014 1,007 1,017 1,014 1,007 1,025 1,022 1,018 0,768 0,884 0,924 0,943 0,963 0,967 0,976 0,990 0,981 1,003 0,991 1,005 0,982 0,991 0,997 0,992 1,008 0,982 0,998 0,764 0,878 0,921 0,951 0,942 0,953 0,975 0,963 0,977 0,965 0,980 0,975 0,973 0,972 0,975 0,979 0,987 0,974 0,983 10 0,772 0,844 0,891 0,922 0,938 0,927 0,955 0,958 0,958 0,952 0,948 0,961 0,960 0,956 0,962 0,969 0,961 0,965 0,964 11 0,759 0,844 0,888 0,880 0,923 0,921 0,923 0,947 0,937 0,939 0,945 0,928 0,939 0,943 0,943 0,940 0,940 0,945 0,939 12 0,793 0,850 0,865 0,882 0,898 0,909 0,890 0,916 0,913 0,911 0,912 0,919 0,925 0,923 0,914 0,928 0,918 0,921 0,918 12,5 0,752 0,831 0,859 0,857 0,889 0,872 0,897 0,903 0,909 0,899 0,905 0,907 0,899 0,908 0,908 0,903 0,921 0,904 0,902 Bảng A8 Giá trị DF(r, θ) nguồn Ir bao stainless steel, sử dụng phổ Ir192_best_Industries.spectrum theo AAPM TG-43 Góc (độ) 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 0,759 0,819 0,909 0,928 0,953 0,965 0,978 0,985 0,986 0,989 0,993 0,998 0,996 1,001 1,007 1,004 1,002 1,004 1,000 1,5 0,727 0,823 0,886 0,931 0,937 0,967 0,968 0,979 0,983 0,989 0,995 0,992 1,001 1,002 0,995 1,028 1,000 0,998 1,000 0,713 0,817 0,897 0,930 0,952 0,961 0,965 0,973 0,989 0,986 0,988 0,989 0,989 1,003 0,997 0,998 1,001 1,000 1,000 2,5 0,706 0,830 0,897 0,920 0,950 0,967 0,971 0,974 0,978 0,987 0,992 0,997 0,988 1,002 0,999 1,002 1,000 0,995 1,000 0,721 0,833 0,896 0,930 0,950 0,966 0,975 0,990 0,986 0,994 0,995 0,999 0,998 1,003 1,002 1,006 1,003 0,996 1,000 0,740 0,863 0,928 0,935 0,957 0,974 0,970 0,986 0,989 0,999 1,003 0,992 1,001 1,004 1,002 1,002 1,005 1,008 1,000 Khoảng cách (cm) 0,753 0,701 0,762 0,843 0,868 0,882 0,903 0,921 0,917 0,945 0,941 0,948 0,954 0,980 0,955 0,958 0,964 0,959 0,977 0,973 0,969 0,979 0,995 0,982 0,987 0,982 0,969 0,986 0,997 0,989 0,994 0,992 0,986 0,992 0,987 0,982 0,986 0,996 0,983 0,997 0,995 0,983 1,010 0,996 0,980 0,990 0,997 0,987 1,005 0,993 0,989 0,992 0,993 0,989 1,000 1,000 1,000 0,770 0,886 0,926 0,945 0,965 0,970 0,978 0,992 0,983 1,006 0,993 1,008 0,985 0,993 0,999 0,994 1,010 0,985 1,000 0,778 0,893 0,937 0,968 0,958 0,969 0,992 0,979 0,994 0,982 0,997 0,991 0,990 0,989 0,992 0,996 1,004 0,990 1,000 10 0,801 0,875 0,924 0,956 0,973 0,961 0,990 0,994 0,994 0,988 0,983 0,996 0,995 0,991 0,998 1,005 0,997 1,001 1,000 11 0,808 0,899 0,945 0,937 0,983 0,980 0,983 1,009 0,998 1,000 1,006 0,988 1,000 1,004 1,004 1,001 1,001 1,006 1,000 12 0,863 0,926 0,941 0,960 0,977 0,990 0,969 0,997 0,994 0,991 0,993 1,001 1,007 1,005 0,995 1,010 0,999 1,002 1,000 12,5 0,834 0,922 0,952 0,951 0,985 0,967 0,995 1,002 1,008 0,997 1,003 1,005 0,997 1,007 1,007 1,001 1,021 1,002 1,000 Bảng A9 Giá trị DF(r, θ) nguồn Cs theo công thức báo Góc (độ) 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 0,796 0,815 0,891 0,917 0,956 0,966 0,981 0,987 0,994 0,996 1,001 1,000 0,998 0,998 1,002 0,998 0,999 0,999 1,000 1,5 0,777 0,817 0,857 0,909 0,933 0,960 0,969 0,974 0,984 0,906 0,987 0,993 0,997 0,997 0,998 1,002 0,997 0,999 0,998 0,788 0,799 0,857 0,898 0,937 0,940 0,965 0,972 0,980 0,981 0,985 0,987 0,990 0,994 0,996 0,994 0,998 0,999 0,999 2,5 0,792 0,804 0,849 0,890 0,924 0,950 0,958 0,968 0,973 0,983 0,980 0,989 0,987 0,990 0,987 0,989 0,993 0,986 0,987 0,801 0,811 0,861 0,894 0,919 0,938 0,945 0,964 0,968 0,971 0,966 0,984 0,981 0,980 0,980 0,983 0,990 0,990 0,988 0,763 0,819 0,842 0,882 0,915 0,946 0,939 0,955 0,964 0,963 0,971 0,974 0,974 0,969 0,972 0,976 0,973 0,978 0,980 Khoảng cách (cm) 0,832 0,798 0,841 0,796 0,841 0,825 0,878 0,861 0,898 0,897 0,921 0,903 0,931 0,919 0,945 0,923 0,943 0,933 0,945 0,933 0,953 0,943 0,958 0,941 0,955 0,946 0,956 0,948 0,958 0,949 0,965 0,952 0,969 0,949 0,966 0,950 0,968 0,956 0,818 0,809 0,836 0,853 0,876 0,895 0,909 0,911 0,915 0,928 0,920 0,931 0,932 0,929 0,932 0,933 0,935 0,931 0,933 0,864 0,751 0,812 0,844 0,857 0,875 0,896 0,897 0,907 0,911 0,912 0,919 0,918 0,923 0,923 0,927 0,924 0,918 0,926 0,765 0,763 0,801 0,838 0,837 0,858 0,883 0,890 0,890 0,891 0,894 0,895 0,900 0,904 0,901 0,901 0,906 0,911 0,913 10 0,696 0,744 0,815 0,814 0,835 0,852 0,853 0,861 0,856 0,864 0,872 0,880 0,876 0,874 0,877 0,885 0,894 0,883 0,886 11 0,722 0,759 0,781 0,806 0,806 0,835 0,836 0,855 0,849 0,844 0,852 0,853 0,868 0,865 0,873 0,884 0,871 0,858 0,868 12 0,773 0,738 0,761 0,790 0,803 0,821 0,831 0,835 0,828 0,841 0,845 0,836 0,841 0,847 0,837 0,851 0,856 0,851 0,849 12,5 0,725 0,738 0,760 0,775 0,783 0,783 0,807 0,811 0,820 0,812 0,832 0,818 0,831 0,828 0,838 0,840 0,838 0,835 0,837 Bảng A10 Giá trị DF(r, θ) nguồn Cs theo AAPM TG-43 Góc (độ) 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 0,772 0,791 0,865 0,891 0,931 0,943 0,960 0,967 0,977 0,981 0,989 0,990 0,991 0,993 0,999 0,996 0,998 0,999 1,000 1,5 0,768 0,808 0,847 0,900 0,924 0,951 0,961 0,967 0,978 0,901 0,984 0,990 0,996 0,997 0,997 1,003 0,999 1,000 1,000 0,783 0,794 0,852 0,893 0,931 0,936 0,961 0,968 0,977 0,978 0,984 0,987 0,989 0,994 0,996 0,995 0,999 1,000 1,000 2,5 0,798 0,811 0,856 0,898 0,932 0,958 0,967 0,977 0,982 0,994 0,991 1,000 0,999 1,002 0,999 1,002 1,006 0,999 1,000 0,809 0,819 0,869 0,902 0,928 0,947 0,955 0,974 0,978 0,982 0,977 0,995 0,992 0,992 0,993 0,996 1,002 1,003 1,000 0,777 0,834 0,858 0,898 0,931 0,964 0,957 0,973 0,982 0,982 0,989 0,992 0,994 0,989 0,992 0,996 0,993 0,998 1,000 Khoảng cách (cm) 0,859 0,834 0,876 0,867 0,832 0,866 0,868 0,862 0,895 0,906 0,900 0,914 0,926 0,937 0,938 0,950 0,944 0,959 0,961 0,960 0,975 0,975 0,965 0,977 0,974 0,977 0,981 0,975 0,976 0,995 0,985 0,986 0,987 0,990 0,985 0,998 0,986 0,989 1,000 0,988 0,992 0,996 0,990 0,993 0,999 0,997 0,996 1,001 1,001 0,993 1,003 0,998 0,994 0,998 1,000 1,000 1,000 0,933 0,810 0,876 0,911 0,925 0,945 0,968 0,969 0,979 0,983 0,985 0,992 0,991 0,997 0,997 1,000 0,998 0,991 1,000 0,838 0,836 0,877 0,918 0,917 0,941 0,967 0,975 0,975 0,976 0,980 0,981 0,986 0,990 0,987 0,988 0,993 0,999 1,000 10 0,785 0,840 0,920 0,919 0,943 0,962 0,963 0,972 0,967 0,976 0,984 0,993 0,988 0,986 0,990 0,999 1,009 0,997 1,000 11 0,832 0,874 0,900 0,928 0,929 0,962 0,964 0,985 0,978 0,972 0,982 0,983 1,000 0,996 1,006 1,019 1,004 0,988 1,000 12 0,910 0,869 0,896 0,930 0,945 0,967 0,979 0,983 0,975 0,989 0,994 0,984 0,990 0,997 0,986 1,001 1,008 1,002 1,000 12,5 0,867 0,882 0,908 0,926 0,935 0,935 0,964 0,969 0,980 0,970 0,994 0,977 0,992 0,990 1,001 1,003 1,002 0,998 1,000 TÀI LIỆU THAM KHẢO [1] Bruce R Thomadsen, Jeffrey F Williamson, Mark J Rivard, Ali S Meigooni (2008), Anniversary Paper: Past and current issues, and trends in brachytherapy physics, American Association of Physicists in Medicine [2] C A F Joslin, A Flynn, E J Hall (2001), Principles and Practice of Brachytherapy - Using Afterloading systems, Arnold [3] C Thomason, T R Mackie, M J Lindstrom and P D Higgins (1991), The dose distribution surrounding 192Ir and 137Cs seed sources, Phys Med Biod., 36(4), IOP Publishing Ltd [4] D Baltas, L Sakelliou, N Zamboglou (2007), The Physics of Modern Brachytherapy for Oncology - Series in Medical Physics and Biomedical Engineering, Taylor & Francis Group, LLC [5] D W O Rogers, I Kawrakow, J P Seuntjens, B R B Walters, E Mainegra-Hing (2003), NRCC Report PIRS-702, NRC Canada [6] Ervin B Podgorsak (2005): Radiation Oncology Physics: A handbook for Teachers and Students, IAEA, Vienna [7] Gultekin Yegin (2003), A new approach to geometry modeling for Monte carlo particle transport: An application to the EGS code system, Elsevier B.V [8] Habib Zaidi, George Sgouros (2003), Therapeutic Applications of Monte Carlo in Nuclear Medicine, IOP Publishing Ltd [9] Halperin, Edward C et al (2008), Perez and Brady's Principles and Practice of Radiation Oncology, 5th Edition, Lippincott Williams & Wilkins [10] IAEA pub 1296 (2008), Setting up a Radiotherapy Programme, IAEA, Vienna [11] IAEA (2000) DMRP/ARBR-2000-01, DIRAC - Directory of Radiotherapy Centres, IAEA, Vienna, [12] IAEA (2006) Safety Reports Series No 47, Radiation protection in the design of radiotherapy facilities, IAEA, Vienna [13] IAEA IAEA-TECDOC-1274 (2002), Calibration of photon and beta ray sources used in brachytherapy, IAEA, Vienna [14] ICRU (2004) Report No.72, Dosimetry of beta rays and low-energy photons for brachytherapy with sealed sources, International Commission on Radiation Units and Measurements [15] I Kawrakow, D W O Rogers (2001–2006), NRCC Report PIRS-701, NRC Canada [16] Jose Perez-Calatayud, Domingo Granero Cabañero, Facundo Ballester Pallarés (2009), Monte Carlo Application in Brachytherapy Dosimetry, Springer Science + Business Media B.V [17] L L Meisberger, R J Keller, and R J Shalek (1968), The effective attenuation in water of the gamma rays of gold-198, iridium-192, cesium-137, radium-226, and cobalt-60 Radiology 90(5), 953–957 [18] Malvin H Kalos, Paula A Whitlock (2004), Monte Carlo methods, WiILEY-VCH Verlag GmbH & Co KGaA, Weinheim [19] Mark J Rivard, Jack L M Venselaar, Luc Beaulieu (2009), The evolution of brachytherapy treatment planning, American Association of Physicists in Medicine [20] Nath et al (1995), Dosimetry of interstitial brachytherapy sources, Med Phys 22 (2), Feb 1995 [21] P Mayles, A Nahum, J C Rosenwald (2007), Handbook of Radiotherapy Physics - Theory and Practice, Taylor & Francis Group, LLC [22] R E P Taylor, G Yegin and D.W.O.Rogers (2006), Benchmarking BrachyDose: Voxel based EGSnrc Monte Carlo calculations of TG-43 dosimetry parameters, American Association of Physicists in Medicine Trang web liên quan [23] http://www.irs.inms.nrc.ca/papers/PIRS629r/node9.html [24] http://www.baomoi.com/Info/Nguoi-benh-se-duoc-tiep-can-voi-ky-thuat-dieu-tri-ung-thu-tothon/82/4550629.epi [25] http://www.vietduchospital.edu.vn/news_detail.asp?ID=2&CID=2&IDN=4541 ... luận văn tìm hiểu code EGSnrc với tinh thần học hỏi cách sử dụng để áp dụng tính liều xạ trị áp sát Nội dung chủ yếu luận văn tìm hiểu áp dụng EGSnrc việc tính phân bố liều số nguồn sử dụng xạ. .. ngành xạ trị với ưu điểm riêng áp dụng rộng rãi việc điều trị bệnh ung thư Xạ trị áp sát kỹ thuật điều trị dùng nguồn đồng vị phóng xạ áp sát vị trí cần điều trị Hiệu điều trị phụ thuộc vào việc tính. .. tính liều xạ trị áp sát dùng phương pháp Monte Carlo [1] Qua thập kỷ áp dụng kỹ thuật mô Monte Carlo, việc mô phân bố liều quanh nguồn xạ trị áp sát ngày đạt mức độ chi tiết hơn, chẳng hạn việc tính