Tính toán che chắn cho phòng xạ trị dùng máy gia tốc bằng monte carlo code EGSnrc

Tính toán che chắn cho phòng xạ trị dùng máy gia tốc bằng monte carlo code EGSnrc Tính toán che chắn cho phòng xạ trị dùng máy gia tốc bằng monte carlo code EGSnrc Tính toán che chắn cho phòng xạ trị dùng máy gia tốc bằng monte carlo code EGSnrc Tính toán che chắn cho phòng xạ trị dùng máy gia tốc bằng monte carlo code EGSnrc Tính toán che chắn cho phòng xạ trị dùng máy gia tốc bằng monte carlo code EGSnrc Tính toán che chắn cho phòng xạ trị dùng máy gia tốc bằng monte carlo code EGSnrc Tính toán che chắn cho phòng xạ trị dùng máy gia tốc bằng monte carlo code EGSnrc Tính toán che chắn cho phòng xạ trị dùng máy gia tốc bằng monte carlo code EGSnrc Tính toán che chắn cho phòng xạ trị dùng máy gia tốc bằng monte carlo code EGSnrc Tính toán che chắn cho phòng xạ trị dùng máy gia tốc bằng monte carlo code EGSnrc Tính toán che chắn cho phòng xạ trị dùng máy gia tốc bằng monte carlo code EGSnrc Tính toán che chắn cho phòng xạ trị dùng máy gia tốc bằng monte carlo code EGSnrc Tính toán che chắn cho phòng xạ trị dùng máy gia tốc bằng monte carlo code EGSnrc

N O N N UYỄN Ị ẨM Ú N ÍN OÁN E ẮN O P ÒN X Ị DÙN MÁY Ố ẲN BẰN MON E LO CODE EGSnrc LU N N L – ăm 2013 N N MỤC LỤC DANH MỤC CÁC KÍ HIỆU VÀ CHỮ VIẾT TẮT DANH MỤC CÁC BẢNG DANH MỤC HÌNH VẼ MỞ ĐẦU CHƯƠNG AN TOÀN BỨC XẠ TRONG XẠ TRỊ 12 1.1 Các khuyến cáo tiêu chuẩn an toàn xạ (ATBX) xạ trị 12 1.1.1 Các khuyến cáo tiêu chuẩn ATBX theo ICRP 12 1.1.1.1 Các đại lượng liều áp dụng tính tốn che chắn 13 1.1.1.2 Giới hạn liều người 17 1.1.2 Các qui định ATBX Việt Nam 18 1.1.2.1 Giới hạn liều nghề nghiệp 18 1.1.2.2 Giới hạn liều công chúng 19 1.2 Những biện pháp nhằm hạn chế tiếp xúc với chùm tia xạ 19 1.3 Cơ sở lý thuyết tính tốn che chắn ATBX 22 1.3.1 Các thuật ngữ tính tốn thiết kế che chắn 22 1.3.1.1 Khu vực kiểm sốt khu vực khơng kiểm sốt 22 1.3.1.2 Khối lượng công việc (Workload-W) 22 1.3.1.3 Hệ số sử dụng (Use Factor-U) 23 1.3.1.4 Hệ số chiếm (Occupancy Factor-T) 23 1.3.1.5 Bức xạ sơ cấp xạ thứ cấp 24 1.3.1.6 Rào chắn bảo vệ 24 1.3.2 Lý thuyết tính tốn che chắn sơ cấp 25 1.3.2.1 Hệ số truyền qua rào sơ cấp ( B pri ) 25 1.3.2.2 Số lớp che chắn (n) 26 1.3.2.3 Bề dày tường che chắn (t) 26 1.3.2.4 Hệ số truyền qua (B) biết trước bề dày tường che chắn 27 1.3.2.5 Liều tương đương (Hpri) truyền qua bề dày vật liệu che chắn 27 1.3.3 Lý thuyết tính tốn che chắn thứ cấp 27 1.3.3.1 Hệ số truyền qua rào thứ cấp xạ tán xạ từ bệnh nhân 27 1.3.3.2 Liều tương đương bên ngồi phòng máy tán xạ 28 1.3.3.3 Hệ số truyền qua rào thứ cấp xạ tán xạ bị rò rỉ  BL  28 1.3.3.4 Liều tương đương bên ngồi phòng máy xạ rò rỉ 28 1.3.3.5 Bề dày lớp che chắn chùm thứ cấp 28 CHƯƠNG ĐỐI TƯỢNG NGHIÊN CỨU VÀ PHƯƠNG PHÁP TÍNH 29 2.1 Máy gia tốc tuyến tính HPD Bệnh viện Chợ Rẫy 30 2.1.1 Nguyên tắc hoạt động 30 2.1.2 Cấu tạo đầu máy gia tốc 31 2.1.3 Tương tác xạ phát từ máy gia tốc với vật chất 32 2.1.4 Phương thức xạ trị đầu máy gia tốc HPD BVCR 33 2.2 Chương trình EGSnrc (Electron Gamma Shower) 34 2.2.1 Nguyên tắc tính liều xạ theo phương pháp mô MCNP 34 2.2.2 Nguyên tắc hoạt động code EGSnrc 34 2.2.3 Cấu trúc chung code EGS 35 2.2.4 Hệ thống thư mục code EGSnrc 36 2.2.5 Mối quan hệ BEAMnrc DOSXYZnrc 37 2.3 Vai trò BEAMnrc mô 38 2.3.1 Khai báo cho BEAMnrc 38 2.3.2 File không gian pha *.egsphant 39 2.4 Vai trò DOSXYZnrc tính tốn phân bố liều 41 2.4.1 Khai báo phantom 41 2.4.2 Các định dạng file Output 42 CHƯƠNG ÁP DỤNG CODE EGSnrc TRONG MƠ PHỎNG VÀ TÍNH TỐN SO SÁNH KẾT QUẢ MƠ PHỎNG GIỮA CODE EGSnrc VỚI KẾT QUẢ TÍNH TỐN THEO NCRP 151 VÀ KẾT QUẢ ĐO TẠI BVCR 44 3.1 Tính phân bố liều khu vực công chúng theo NCRP 151 44 3.1.1 Điều kiện làm việc liệu tính tốn 45 3.1.2 Kết tính phân bố liều theo NCRP 151 cho khu vực công chúng 46 3.1.2.1 Liều bên tường sơ cấp  H pri  46 3.1.2.2 Liều bên tường thứ cấp  H sec  46 3.1.2.3 Tổng liều xạ bên ngồi phòng máy gia tốc 47 3.1.3 Kết đo phân bố liều phòng máy gia tốc BVCR 47 3.2 So sánh phương pháp tính giải tích theo NCRP 151 Siemens 48 3.2.1 Phương pháp tính liều giới hạn theo NCRP 151 48 3.2.2 Phương pháp tính liều giới hạn BVCR theo Siemens 51 3.3 Tính phân bố liều phòng máy gia tốc dùng code EGSnrc 53 3.3.1 Mơ hình hố đầu máy gia tốc dùng BEAMnrc 53 3.3.2 Tạo file không gian pha 55 3.3.2.1 Nguồn sử dụng BEAMnrc 55 3.3.2.2 Các thông số mô tả nguồn 57 3.3.3 Q trình tính liều với DOSXYZnrc 58 3.3.3.1 Khai báo phantom 59 3.3.3.2 Nguồn thông số vận chuyển 60 3.3.3.3 Trình bày kết 62 KẾT LUẬN 68 HƯỚNG PHÁT TRIỂN 70 TÀI LIỆU THAM KHẢO 71 PHỤ LỤC 74 DANH MỤC CÁC KÍ HIỆU VÀ CHỮ VIẾT TẮT 2D: Two-Dimensional 3D: Three-Dimensional ATBX: An Toàn Bức Xạ BVCR: Bệnh viện Chợ Rẫy CM: Component Module ECUT: Electron cut-off energy EGS: Electron Gamma Shower ETRAN: Electron TRANsport FS: Field Size GEANT4: Geometry Add Tracking IAEA: International Atomic Energy Agency ICRP: International Commission on Radiological Protection MCNP: Monte Carlo N Particle NCRP: National Council on Radiation Protection and Measurements NRCC: National Research Council of Canada OMEGA: Ottawa Madision Electron Gamma Algorithm PCUT: Photon Cut-Off Energy PENELOPE: Penetration and Energy LOss of Positrons and Electron SAD: Source to Axis Distance SSD: Source to Surface Distance W: Workload U: Use Factor T: Occupancy Factor TVL: Tenth Value Layer TVLe: Equilibrium Tenth Value Layer TVL1: First Tenth Value Layer 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) 15 Bảng 1.2: Các trọng số mô đặc trưng cho mô quan thể 16 Bảng 1.3: Giới hạn liều qua thời kỳ ICRP 18 Bảng 1.4: Giá trị hệ số sử dụng chùm tia sơ cấp 23 Bảng 1.5: Giá trị hệ số chiếm 24 Bảng 3.1: Kết đo liều Bệnh Viện Chợ Rẫy 48 Bảng 3.2: Bề dày che chắn cho khu vực nhân viên tính theo NCRP 151 Siemens 53 Bảng 3.3: Khoảng lượng xác suất tương ứng với photon MV 56 Bảng 3.4: Mảng giá trị liều khu vực nhân viên đặt phantom nước 62 tâm phòng máy gia tốc 62 Bảng 3.5: Mảng giá trị liều khu vực cơng chúng đặt phantom nước tâm phòng máy 63 Bảng 3.6: So sánh giá trị liều (mSv/năm) khu vực nhân viên 64 Bảng 3.7: So sánh giá trị liều (mSv/năm) khu vực công chúng 64 Bảng 3.8: Mảng giá trị liều thay phantom nước khơng khí loại bỏ tường che chắn 65 Bảng 3.9: So sánh suy giảm chùm tia theo lý thuyết (quy luật giảm theo bình phương khoảng cách) mô 66 Bảng 3.10: Mảng giá trị liều thay phantom nước khơng khí phòng máy che chắn với tường dày 130 cm 66 Bảng 3.11: So sánh suy giảm chùm tia theo quy luật hàm logarit mơ 67 DANH MỤC HÌNH VẼ Hình 1.1: Biểu diễn độ suy giảm cường độ chùm tia theo khoảng cách 21 Hình 1.2: Biểu diễn phụ thuộc độ rộng chùm tia vào bề dày lớp che chắn 21 Hình 1.3: Sơ đồ che chắn xạ sơ cấp thứ cấp 25 Hình 1.4: Khoảng cách từ nguồn đến điểm Q bên tường bảo vệ 26 Hình 1.5: Hình biểu diễn khoảng cách d sca , d sec , F d L 27 Hình 2.1: Sơ đồ nguyên tắc hoạt động máy gia tốc 30 Hình 2.2: Thành phần đầu điều trị máy gia tốc tuyến tính 31 Hình 2.3: Sơ đồ tạo thành hiệu ứng xạ tương tác với vật chất 33 Hình 2.4: Cấu trúc hệ thống code EGSnrc sử dụng với user code 35 Hình 2.5: Cấu trúc đường dẫn HEN_HOUSE 36 Hình 2.6: Cấu trúc đường dẫn EGS_HOME 37 Hình 2.7: Mối liên hệ BEAMnrc DOSXYZnrc 38 Hình 2.8: Giao diện CM đầu máy gia tốc HPD 38 Hình 2.9: Giao diện cho khai báo thông số vận chuyển, thơng số vào 39 Hình 2.10: Giao diện khai báo vùng ghi file không gian pha 40 Hình 2.11: Giao diện khai báo Phantom DOSXYZnrc 42 Hình 3.1: Sơ đồ vị trí tính liều bên ngồi phòng máy gia tốc 44 Hình 3.2: Các khoảng cách tính tốn che chắn 45 Hình 3.3: Hướng chiếu chùm tia ứng với góc quay khác 49 Hình 3.4: Mơ hình 2D đầu máy gia tốc HPD 55 Hình 3.5: Giao diện nguồn khai báo cho Beamnrc 58 Hình 3.6: Giao diện khai báo vùng tính liều 60 Hình 3.7: Giao diện khai báo vật liệu cho phantom 60 Hình 3.8: Giao diện khai báo thông số cho nguồn DOSXYZnrc 61 Hình 3.9: Tọa độ File khơng gian pha đến từ nguồn 61 Hình 3.10: Đồ thị biểu diễn suy giảm liều phantom nước đặt phòng máy 62 Hình 3.11: Đồ thị biểu diễn suy giảm liều phòng máy gia tốc khơng che chắn 65 Hình 3.12: Đồ thị biểu diễn suy giảm liều thay phantom nước 67 MỞ ĐẦU Nhằm thực chương trình phòng chống ung thư quốc gia, nhiều bệnh viện Việt Nam xúc tiến việc lắp đặt máy gia tốc cho xạ trị Do suất liều cung cấp máy gia tốc cao, phòng xạ trị phải che chắn thích hợp để bảo đảm mặt an tồn xạ cho nhân viên cơng chúng Việc thiết kế che chắn cho phòng đặt máy tốn vừa mang tính khoa học vừa có yếu tố kinh tế: mặt phải đáp ứng qui định liều giới hạn, mặt khác cần hợp lý mặt chi phí Để giải tốn này, cần có thơng tin suất liều, tần suất phát tia, số lượng bệnh nhân, góc chiếu trung bình, thời gian lưu trú nhân viên hay cơng chúng sở, bố trí khu vực lân cận, v.v Cuối cùng, biết phổ lượng chùm xạ, ta mơ tả suy giảm chùm tia qua bề dày vật chất khác Bên cạnh phương pháp tính tốn giải tích trình bày tài liệu IAEA [10] hay NCRP 151 [14], phương pháp Monte Carlo cho phép tính tốn che chắn Đặc biệt, chương trình Monte Carlo EGSnrc, với code chuyên dụng BEAMnrc DOSXYZnrc, xây dựng cho tính tốn với máy gia tốc, giúp tính tốn che chắn cho phòng máy gia tốc Cơng cụ công nhận rộng rãi giới vật lý y khoa xem tiêu chuẩn để đánh giá tính tốn liều lượng Để mô tả phổ lượng chùm xạ nhờ EGSnrc tính tốn che chắn, cần phải có hiểu biết cấu trúc máy gia tốc Với máy gia tốc cụ thể điều khả thi Luận văn nhằm mục đích nghiên cứu áp dụng EGSnrc vào việc tính tốn che chắn cho phòng chứa máy gia tốc Kết tính tốn so sánh với phương pháp tính khác để đánh giá độ xác rút kết luận cần thiết Từ mục đích nội dung cơng việc trên, luận văn bố cục bao gồm ba chương 154, 130, 30, -398 30, 75, 55, 208, 20, 10, 20, 130, 130, 30, 1, 1, 1, 11, 1, 10, 2, 0.001 2, 2, 1, 1, 1, 10, 2, 0.001 2, 2, 2, 10, 1, 10, 3, 3.2 2, 2, 11, 11, 1, 10, 2, 0.001 3, 3, 1, 1, 1, 10, 2, 0.001 3, 3, 2, 2, 1, 10, 3, 3.2 3, 3, 3, 3, 1, 10, 2, 0.001 3, 3, 4, 8, 1, 10, 3, 3.2 3, 3, 9, 9, 1, 10, 2, 0.001 3, 3, 10, 10, 1, 10, 3, 3.2 3, 3, 11, 11, 1, 10, 2, 0.001 4, 4, 1, 1, 1, 10, 2, 0.001 4, 4, 2, 2, 1, 10, 3, 3.2 4, 4, 3, 9, 1, 10, 2, 0.001 4, 4, 10, 10, 1, 10, 3, 3.2 85 4, 4, 11, 11, 1, 10, 2, 0.001 5, 7, 1, 1, 1, 10, 2, 0.001 5, 7, 2, 2, 1, 10, 3, 3.2 5, 7, 3, 4, 1, 10, 2, 0.001 5, 7, 5, 7, 1, 10, 1, 5, 7, 8, 9, 1, 10, 2, 0.001 5, 7, 10, 10, 1, 10, 3, 3.2 5, 7, 11, 11, 1, 10, 2, 0.001 8, 8, 1, 1, 1, 10, 2, 0.001 8, 8, 2, 2, 1, 10, 3, 3.2 8, 8, 3, 9, 1, 10, 2, 0.001 8, 8, 10, 10, 1, 10, 3, 3.2 8, 8, 11, 11, 1, 10, 2, 0.001 9, 9, 1, 1, 1, 10, 2, 0.001 9, 9, 2, 2, 1, 10, 3, 3.2 9, 9, 3, 3, 1, 10, 2, 0.001 9, 9, 4, 8, 1, 10, 3, 3.2 9, 9, 9, 9, 1, 10, 2, 0.001 9, 9, 10, 10, 1, 10, 3, 3.2 9, 9, 11, 11, 1, 10, 2, 0.001 10, 10, 1, 1, 1, 10, 2, 0.001 10, 10, 2, 10, 1, 10, 3, 3.2 10, 10, 11, 11, 1, 10, 2, 0.001 11, 11, 1, 11, 1, 10, 2, 0.001 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 6, 6, 1, 11, 6, 6, 0, 1, 11, 6, 6, 6, 6, 2, 0, 0, 0, 0, 0, 0, 0, 86 2, 2, 0, 0, 0, 90, 180, 0, 180, 0, 0, 0, 0, 2, 2, 2, 30, 0, 30, 10, 30 C:/egsnrc_mp/BEAM_MGT6MV/6MVCR_Spectrum.egsphsp1 400000000, 0, 300, 33, 97, 40, 0, 1, 0, 0, , 0, 0, 0, 0, ######################### :Start MC Transport Parameter: Global ECUT= 0.7 Global PCUT= 0.01 Global SMAX= 1e10 ESTEPE= 0.25 XIMAX= 0.5 Boundary crossing algorithm= EXACT Skin depth for BCA= Electron-step algorithm= PRESTA-II Spin effects= On Brems angular sampling= Simple Brems cross sections= BH Bound Compton scattering= Off Compton cross sections= default Pair angular sampling= Simple Pair cross sections= BH Photoelectron angular sampling= Off Rayleigh scattering= Off Atomic relaxations= Off Electron impact ionization= Off Photon cross sections= xcom Photon cross-sections output= Off :Stop MC Transport Parameter: ######################### 87 PHỤ LỤC C: file output *.egslst với phantom nước đặt tâm phòng máy gia tốc ******************************************************************* NRCC/UW EGSnrc user-code DOSXYZnrc ($Revision: 1.50 $ last edited $Date: 2013/02/13 21:57:55 $) ON gnu-win32 (gnu_win32) 07:43:53 Aug 24 2013 ******************************************************************* ** DOSXYZnrc ** Z pronounced zed ** ** ** Code developed at the NationalResearch Council of Canada and ** ** University of Wisconsin as part of the OMEGA project ** This is $Revision: 1.50 $ last edited $Date: 2013/02/13 21:57:55** ** ******************************************************************* The following parameters may be adjusted in dosxyz_user_macros.mortran $MXMED: Max number of media: $MXSTACK: Max stack size: 10000 $IMAX,etc: Max dose scoring regions in x,y,z directions: 128 128 128 $MAXDOSE: Max dose scoring regions consistent with above:2097153 $DOSEZERO(=0) 1=> all doses with uncert > 50% are zeroed in 3ddose file The following parameters may be adjusted in srcxyz.macros $INVDIM: number of elements in inverse CPD for input energy spectra = 1000 $NENSRC: number of bins in input energy spectrum = 200 Title: INPUT-DOSXYZnrc-MGT6MV =========================================================== Number of media (min = 1, max = 7, => CT data): Medium 1: H2O700ICRU Medium 2: AIR700ICRU 88 Medium 3: Concrete ECUTIN,PCUTIN,(ESTEPE,SMAX DUMMY INPUTS): 0.700 0.010 0.000 0.000 0.000 0.000 # regions in x (max= 128),y (max= 128),z (max= 128) directions (if end of this input Input groups of regions for which ECUT and PCUT are not defaults NB This option is disabled, just input zeros Dummy values of lower,upper i, j, k, ECUT, PCUT Found blank line => end of this input Enter numbers on one line pairs defininglower,upper x,y,z indicies of dose regions for which results are to be output IZSCAN: non-zero for z-scan/page MAX20: if any one = 1, output summary of max 20 doses end signaled by first pair both zero forno dose printed, MAX20 is still read from first line 6 11 6 0 11 6 6 Found blank line => end of this input The material in the region outside the phantom is:AIR700ICRU The thickness of this region (in x, y & z direction) is: 30.000 cm Particles will be read from file: 92 C:/egsnrc_mp/BEAM_MGT6MV/6MVCR_Spectrum.egsphsp1 Total number of particles in file Total number of photons : : 999165601 997900575 The rest are electrons/positrons Maximum kinetic energyof the particles: 6.985 MeV Minimum kinetic energy ofthe electrons: 0.189 MeV # of particles incident fromoriginal source: 4049999872.0 NCASE,IWATCH,TIMMAX,INSEED1,INSEED2,BEAM_SIZE,ISMOOTH,IRES TART,IDAT,REJECT,ESAVE_GLOBAL,NRCYCL,IPARALLEL,PARNUM,n_sp lit,ihowfarless: 400000000 300.00 33 97 40.00 0 0.00 0 0 Above RNG seeds will be replaced by those from previous run =========================================================== Electron/Photon transport parameter =========================================================== Photon cross sections xcom Compton cross sections default Photon transport cutoff(MeV) 0.1000E-01 Pair angular sampling SIM Pair cross sections BH Triplet production Off Bound Compton scattering OFF Radiative Compton corrections Off Rayleigh scattering OFF Atomic relaxations OFF Photoelectron angular sampling OFF Electron transport cutoff(MeV) 0.7000 Bremsstrahlung cross sections BH 93 Bremsstrahlung angular sampling SIM Spin effects On Electron Impact Ionization Off Maxium electron step in cm (SMAX) 0.1000E+11 Maximum fractional energy loss/step (ESTEPE) 0.2500 Maximum 1st elastic moment/step (XIMAX) 0.5000 Boundary crossing algorithm EXACT Skin-depth for boundary crossing (MFP) Electron-step algorithm 3.000 PRESTA-II =========================================================== ******************************************************************* Summary of source parameters (srcxyznrc $Revision: 1.32 $) ******************************************************************* Full phase space input for each incident particle x-coordinate of the isocenter, 0.0000 cm y-coordinate of the isocenter, 0.0000 cm z-coordinate of the isocenter, 0.0000 cm Polar angle of source plane: 90.0000 degrees Azimuthal angle of origin in source plane: 180.0000 degrees Distance from isocenter to origin in source plane: Source plane rotation angle, 180.0000 degrees Total number of particles in phase space file: Particles to be simulated: all Medium AE AP H2O700ICRU 0.700 0.010 AIR700ICRU 0.700 0.010 Concrete 0.700 0.0000 cm 0.010 Histories finshed in previous run 2000000000 Equivalent no of primary histories 2489781843 94 999165601 Total previou CPU time (hour) 52.134 Random number seeds for restart ixx jxx = 95 31 No range rejection *************************************************************** Histories to be simulated for this run 400000000 Histories to be analyzed after this run 2400000000 *************************************************************** Elapsed wall clock time to this point= CPU time so far for this run = 1.964 s 0.811 s BATCH # TIME-ELAPSED TOTAL CPUTIME RATIO TIME OF DAY RNG pointers 0.0 0.0 0.00 07:43:55 ixx jxx = 95 31 2864.2 2785.9 1.03 08:31:39 ixx jxx = 91 27 6448.5 6227.4 1.04 09:31:23 ixx jxx = 36 69 11040.3 10819.0 1.02 10:47:55 ixx jxx = 14421.3 14172.4 1.02 11:44:16 ixx jxx = 64 97 17286.2 16976.0 1.02 12:32:01 ixx jxx = 83 19 21878.5 21568.0 1.01 13:48:33 ixx jxx = 73 25896.3 25567.1 1.01 14:55:31 ixx jxx = 22 55 28749.9 28356.7 1.01 15:43:05 ixx jxx = 10 32740.1 32322.4 1.01 16:49:35 ixx jxx = 30 63 34 39 Total CPU time for run = 36909.0 s = 10.253 hr => 39014849 hist/hr On gnu-win32 (gnu_win32) CPU time for this+previous run =224590.2 s = 62.386 hr => 38470064 hist/hr On gnu-win32 (gnu_win32) ******************************************************************* # of particles read from ph-sp file (N_read) # of primary (non ph-sp) histories read from ph-sp file 95 = 2400000016 = 2981903490 # of particles discarded dueto charge/LATCH/W/multiple passer = # of particles discarded because beyond BEAM_SIZE = 16 # of photons rejected because beyond DBS splitting radius = # of particles that missed geometry = N_used/N_read = 1.000 # of times each particle in ph-sp file recycled (last particle may be recycled less than this) # of ph-sp particles simulated (N_used) = = 2400000000 # of times ph-sp file restarted in this run = ******************************************************************* Fraction of incident energy deposited in the phantom = 0.9853 Fraction of incident energy deposited in the region surrounding the phantom when incident particles go through it = 0.0001 Number of charged particle steps simulated, N_step = Number of charged particle steps/incident fluence = 98437731947 1.01189E+01 No of PRESTA-II steps/total no of charged particle steps = 0.15609 1INPUT-DOSXYZnrc-MGT6MV Elec/positron planar energy fluence scored in a 1cm**2 circle centered at z-axis on the phantom surface = 1.645E+03(1/cm**2) Photon planar energy fluence scored in a 1cm**2 circle centered at z-axison the phantom surface = 1.258E+06(1/cm**2) Elec/positronplanar fluence scored in a 1cm**2 circle centered at z-axis on the phantom surface = 1.015E+03(1/cm**2) Photon planar fluence scored in a 1cm**2 circle centered at z-axis on the phantom surface = 6.892E+05(1/cm**2) No of particles incident from phase space file = 2400000000 No of particles incident from original source = 9728116736.0 DOSXYZnrc ($Revision: 1.50 $) Dose outputs (Dose/incident particle from original source, Gy) 96 for z= -10.000 to ybounds:******* j= xbounds ( 0.000 -473.500 k= -343.500 -189.500 -25.000 -5.000) 5.000 0.000E+00-99.9% 3.325E-23- 0.5% 2.301E-21- 1.5% 1.200E-200.7% for z= -10.000 to ybounds:-25.000 j= xbounds ( 0.000 -5.000 k= 5.000 25.000 189.500 -5.000) 5.000 1.222E-18- 0.0% 1.096E-17- 0.0% 3.517E-19- 0.2% 7.828E-21- 0.8% for z= -10.000 to ybounds:189.500 j= xbounds ( 0.000 343.500 10 k= 473.500 503.500 11 -5.000) 5.000 1.728E-21- 1.7% 2.220E-23- 0.9% 0.000E+00-99.9% 1INPUT-DOSXYZnrc-MGT6MV Elec/positron planar energy fluence scored in a 1cm**2 circle centered at z-axis on the phantom surface = 1.645E+03(1/cm**2) Photon planar energy fluence scored in a 1cm**2 circle centered at z-axison the phantom surface = 1.258E+06(1/cm**2) Elec/positronplanar fluence scored in a 1cm**2 circle centered at z-axis on the phantom surface = 1.015E+03(1/cm**2) Photon planar fluence scored in a 1cm**2 circle centered at z-axis on the phantom surface = 6.892E+05(1/cm**2) No of particles incident from phase space file = 2400000000 No of particles incident fromoriginal source = 9728116736.0 DOSXYZnrc ($Revision: 1.50 $) Dose outputs (Dose/incident particle from original source, Gy) 97 for z= -10.000 to xbounds:******* i= ybounds ( 0.000 -520.000 k= -445.000 -390.000 -25.000 -5.000) 5.000 0.000E+00-99.9% 0.000E+00-99.9% 1.115E-23- 0.9% 9.405E-221.2% for z= -10.000 to xbounds:-25.000 i= ybounds ( -5.000 0.000 k= 5.000 25.000 390.000 -5.000) 5.000 8.625E-20- 0.0% 1.096E-17- 0.0% 0.493E-18- 0.0% 7.537E-19- 0.1% for z= -10.000 to 0.000 xbounds:390.000 445.000 i= ybounds ( 5.0 10 k= 520.000 550.000 11 -5.000) 1.368E-19- 0.0% 4.448E-22- 0.4% 1.230E-24-59.6% 98 Phụ lục D Bảng: Các điều kiện áp dụng tính tốn mơ Điều kiện tính tốn & mơ Liều giới hạn phép (P) vùng kiểm soát Giá trị/đơn vị 0,1 mSv/tuần (5 mSv/năm) Liều giới hạn phép (P) vùng khơng kiểm sốt 0,02 mSv/tuần (1 mSv/năm) Suất liều Isocenter 2Gy/phút Số lượng bệnh nhân ngày 100 người Liều chiếu cho bệnh nhân Gy Thời gian làm việc máy gia tốc giờ/ngày ngày/tuần Khoảng cách từ nguồn đến Isocenter 100 cm Máy gia tốc phát chùm photon có lượng MV Khoảng cách từ isocenter đến điểm bên tường che d C = 5,50 m chắn cách tường 0.3 m Khoảng cách từ nguồn đến điểm bên tường che d pri = 6,50 m chắn cách tường 0.3 m Workload 1000 Use factor Occupancy factor (khu vực nhân viên) 99 ... đặt máy gia tốc cho xạ trị Do suất liều cung cấp máy gia tốc cao, phòng xạ trị phải che chắn thích hợp để bảo đảm mặt an toàn xạ cho nhân viên cơng chúng Việc thiết kế che chắn cho phòng đặt máy. .. pháp Monte Carlo cho phép tính tốn che chắn Đặc biệt, chương trình Monte Carlo EGSnrc, với code chuyên dụng BEAMnrc DOSXYZnrc, xây dựng cho tính tốn với máy gia tốc, giúp tính tốn che chắn cho phòng. .. trúc máy gia tốc Với máy gia tốc cụ thể điều khả thi Luận văn nhằm mục đích nghiên cứu áp dụng EGSnrc vào việc tính tốn che chắn cho phòng chứa máy gia tốc Kết tính tốn so sánh với phương pháp tính