Áp Dụng Phương Pháp Quét Gamma Cắt Lớp Để Xác Định Vị Trí Nguồn Phóng Xạ Bên Trong Thùng Thải

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM THÀNH PHỐ HỒ CHÍ MINH Vũ Lan Anh ÁP DỤNG PHƯƠNG PHÁP QUÉT GAMMA CẮT LỚP ĐỂ XÁC ĐỊNH VỊ TRÍ NGUỒN PHĨNG XẠ BÊN TRONG THÙNG THẢI LUẬN VĂN THẠC SĨ KHOA HỌC VẬT CHẤT Thành phố Hồ Chí Minh – 2019 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM THÀNH PHỐ HỒ CHÍ MINH Vũ Lan Anh ÁP DỤNG PHƯƠNG PHÁP QUÉT GAMMA CẮT LỚP ĐỂ XÁC ĐỊNH VỊ TRÍ NGUỒN PHĨNG XẠ BÊN TRONG THÙNG THẢI Chuyên ngành: Vật lý nguyên tử Mã số : 60 44 01 06 LUẬN VĂN THẠC SĨ KHOA HỌC VẬT CHẤT NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS TS TRẦN THIỆN THANH Thành phố Hồ Chí Minh – 2019 BỘ GIÁO DỤC VÀ ĐÀO TẠO LỜI CAM ĐOAN Tơi xin cam đoan cơng trình nghiên cứu riêng thân Các số liệu luận văn trung thực Kết luận văn chưa cơng bố cơng trình Tác giả luận văn Vũ Lan Anh LỜI CÁM ƠN Trong trình học tập Trường Đại học Sư Phạm thời gian thực luận văn phịng thí nghiệm Vật lý Hạt nhân, Khoa Vật lý – Vật lý Kỹ thuật, Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Thành phố Hồ Chí Minh, em nhận giúp đỡ tận tình, q giá thầy, cơ, cán anh, chị bạn học viên Thông qua luận văn này, em xin gửi lời cảm ơn sâu sắc đến:  PGS TS Trần Thiện Thanh, người hướng dẫn trực tiếp, tạo điều kiện sở vật chất yếu tố khác để hoàn thành luận văn cách thuận lợi Thầy tận tình bảo, hướng dẫn cặn kẽ giúp đỡ, động viên tháo gỡ khó khăn vướng mắc q trình thực luận văn  ThS Huỳnh Đình Chương, người hướng dẫn thứ hai người anh, ân cần bảo, hướng dẫn nhiệt tình, chi tiết, khơng ngại khó khăn em giải điều cịn vướng mắc q trình thực luận văn bên cạnh bảo thầy Thanh  TS Trần Nhân Giang, giúp đỡ em nhiều ngơn ngữ lập trình q trình dựng ảnh  Bộ môn Vật lý Hạt nhân đáp ứng điều kiện sở vật chất, trang thiết bị cần thiết để em thực luận văn  Các Thầy Cô hội đồng bảo vệ luận văn đọc, nhận xét đóng góp ý kiến giúp luận văn hồn thiện  Các thành viên gia đình dành tất tình yêu thương, hy sinh, lúc bên cạnh, giúp vượt qua khó khăn học tập sống  Phịng sau đại học tạo điều kiện giúp đỡ em để hoàn thành luận văn cách tốt Thành phố Hồ Chí Minh, tháng năm 2019 VŨ LAN ANH MỤC LỤC LỜI CAM ĐOAN LỜI CÁM ƠN DANH MỤC CÁC BẢNG DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ MỞ ĐẦU CHƯƠNG TỔNG QUAN 1.1 Tổng quan đề tài nghiên cứu 1.1.1 Tình hình nghiên cứu giới 1.1.2 Tình hình nghiên cứu nước 1.1.3 Nhận xét 1.1.4 Mục tiêu nội dung nghiên cứu luận văn 1.2 Cơ sở lý thuyết 1.2.1 Phép biến đổi Radon 1.2.2 Kỹ thuật chiếu ngược 11 1.2.3 Kỹ thuật lọc ảnh 13 1.2.4 Kỹ thuật lọc ảnh miền tần số 14 1.3 Tổng kết chương 16 CHƯƠNG THỰC NGHIỆM 17 2.1 Thiết lập thực nghiệm 17 2.2 Quy trình tái tạo ảnh gamma cắt lớp thùng thải 23 2.3 Tổng kết chương 26 CHƯƠNG KẾT QUẢ VÀ THẢO LUẬN 27 3.1 Dữ liệu thu hình chiếu 27 3.2 Kết tái tạo ảnh cắt lớp 28 3.3 Xác định vị trí nguồn phóng xạ 32 3.3 Tổng kết chương 33 KẾT LUẬN VÀ KIẾN NGHỊ 34 TÀI LIỆU THAM KHẢO 35 PHỤ LỤC PL1 DANH MỤC CÁC BẢNG Bảng 2.1 Thông tin nguồn phóng xạ 137Cs sử dụng luận văn 19 DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ Hình 1.1 Mơ tả phép biến đổi Radon 10 Hình 1.2 Ma trận hình chiếu ảnh cắt lớp 11 Hình 1.3 Quá trình ghi nhận hình chiếu tái tạo ảnh 12 Hình 1.4 Sự nhịe ảnh phép chiếu ngược 13 Hình 2.1 Đầu dị nhấp nháy NaI(Tl) sử dụng luận văn 17 Hình 2.2 Cấu hình ống chuẩn trực đầu dị 18 Hình 2.3 Thùng thải sử dụng luận văn 19 Hình 2.4 Vị trí nguồn phóng xạ bên phân đoạn thùng thải 20 Hình 2.5 Hệ đo quét gamma cắt lớp cho thùng thải phóng xạ 21 Hình 2.6 Giao diện phần mềm March điều khiển hoạt động hệ khí 22 Hình 2.7 Sơ đồ thực quét cắt lớp thùng thải 24 Hình 2.8 Phổ gamma ghi nhận diện tích đỉnh 25 Hình 2.9 Sơ đồ mô tả bước xử lý chương trình tái tạo ảnh cắt lớp 26 Hình 3.1 Sự phân bố diện tích đỉnh theo tọa độ đầu dị trục Ox 27 Hình 3.2 Sinogram tương ứng với liệu 28 Hình 3.3 Ma trận hình chiếu tương ứng với liệu 19 hình chiếu 29 Hình 3.4 Ma trận hình chiếu tương ứng với liệu 37 hình chiếu 30 Hình 3.5 Ảnh chụp cắt lớp tương ứng với liệu 19 hình chiếu 31 Hình 3.6 Ảnh chụp cắt lớp tương ứng với liệu 37 hình chiếu 31 Hình 3.7 Ảnh cắt lớp tái tạo kỹ thuật chiếu ngược có lọc 32 MỞ ĐẦU Các hoạt động nghiên cứu ứng dụng công nghệ hạt nhân như: sản xuất điện hạt nhân, điều chế dược phẩm phóng xạ, xạ trị y học chiếu xạ công nghiệp v.v tạo lượng lớn rác thải phóng xạ Những rác thải chứa đồng vị phóng xạ có hoạt độ khác chu kì bán rã từ vài chục năm đến hàng triệu năm Công tác quản lý rác thải phóng xạ ln mối quan tâm quốc gia sử dụng phát triển công nghệ hạt nhân Thơng thường, rác thải phóng xạ chứa đựng bên thùng kín lớn để đảm bảo chúng cách ly hồn tồn với mơi trường bên ngồi Sau đó, hoạt độ đồng vị phóng xạ bên thùng chứa rác thải phóng xạ (sau gọi tắt thùng thải) xác định để phân loại cho phù hợp với quy tắc xử lý quốc gia Việc phân tích hoạt độ đồng vị phóng xạ bên thùng thải thường thực hệ phổ kế gamma Để đạt kết phân tích xác, phân bố đồng vị phóng xạ vật liệu (matrix) bên thùng thải cần phải biết rõ Tuy nhiên, đồng vị phóng xạ phân bố vị trí bên thể tích thùng thải, đồng thời thể tích thùng thải lấp đầy nhiều loại vật liệu khác (đối với trường hợp matrix khơng đồng nhất) Do đó, việc tìm kiếm phương pháp để xác định vị trí đồng vị phóng xạ matrix bên thùng thải cần thiết Phương pháp chụp ảnh gamma cắt lớp ứng dụng rộng rãi cho việc chẩn đốn hình ảnh y tế Phương pháp phân làm hai mơ hình chụp khác Mơ hình thứ có nguồn phóng xạ đầu dị nằm bên ngồi đối tượng cần chụp ảnh để ghi nhận tia gamma truyền qua Mơ hình cho biết thơng tin phân bố vật liệu bên đối tượng Mơ hình cịn lại nguồn phóng xạ đưa vào bên đối tượng, đầu dò đặt bên ngồi để ghi nhận tín hiệu Ảnh từ mơ hình cho thấy phân bố nguồn phóng xạ bên đối tượng Về nguyên lý, phương pháp chụp ảnh gamma cắt lớp ứng dụng cho tốn phân tích thùng thải phóng xạ Do đó, đề tài nghiên cứu xác định phân bố vị trí nguồn phóng xạ vật liệu bên thùng thải phương pháp chụp ảnh gamma cắt lớp cơng việc có ý nghĩa Một điểm luận văn tiến hành khảo sát cho trường hợp có hai nguồn phóng xạ với hoạt độ khác bố trí hai vị trí khác bên thùng thải Luận văn thực với mục tiêu ứng dụng phương pháp chụp ảnh gamma cắt lớp để xác định phân bố nguồn phóng xạ bên thùng thải Trong đó, luận văn tập trung vào việc khai thác thuật tốn chiếu ngược có lọc (filtered back projection algorithm) để dựng ảnh từ liệu ghi nhận hệ đo Vị trí nguồn phóng xạ xác định từ hình ảnh so sánh với vị trí thực tế bên thùng thải để kiểm chứng phương pháp Nội dung nghiên cứu luận văn phần đề tài khoa học công nghệ cấp Đại học Quốc Gia – Hồ Chí Minh mã số C2018-18-04, với tiêu đề “Nghiên cứu phương pháp phân tích khơng phá hủy để xác định phân bố hoạt độ đồng vị phóng xạ bên thùng chứa rác thải hạt nhân có matrix khơng đồng nhất” Nội dung luận văn trình bày 03 chương Trong đó: Chương trình bày tổng quan tình hình nghiên cứu chủ đề liên quan đến luận văn; sở lý thuyết dựng ảnh Chương mô tả chi tiết điều kiện thực nghiệm phép đo hệ phổ kế gamma dùng đầu dị NaI(Tl) thực luận văn; quy trình đo phân tích liệu để thu ảnh chụp cắt lớp Chương trình bày ảnh chụp cắt lớp đạt từ liệu thực nghiệm chương trình dựng ảnh xây dựng; đánh giá ảnh hưởng vài yếu tố lên chất lượng ảnh; so sánh kết dựng ảnh giá trị thực tế CHƯƠNG TỔNG QUAN 1.1 Tổng quan đề tài nghiên cứu 1.1.1 Tình hình nghiên cứu giới Stanga cộng [1] (năm 2012) đưa phương pháp để xác định hoạt độ đồng vị phóng xạ thùng thải kỹ thuật quét gamma Về nguyên lý hoạt động, thùng thải phóng xạ quay liên tục trình phân tích, đầu dị dịch chuyển đến vị trí tương ứng với độ cao khoảng cách với thùng thải khác để ghi nhận phổ gamma cho vị trí đo Về mặt ý tưởng phương pháp, thể tích thùng thải chia thành hệ thống voxel tạo nên hình trụ bên vịng trụ bên Cho voxel đồng vị phóng xạ phân bố đồng có hoạt độ chưa biết Đối với vị trí đo đầu dò, tốc độ đếm đỉnh lượng quan tâm từ phép đo thùng thải tính tổng tốc độ đếm voxel Trong đó, tốc độ đếm cho voxel xác định tích hoạt độ, hiệu suất ghi nhận hệ số hiệu chỉnh voxel Trong trường hợp matrix thùng thải gần đồng hệ số hiệu chỉnh cho Từ phép đo nhiều vị trí khác đầu dị mức lượng khác xạ gamma phát từ đồng vị hệ phương trình đa biến thiết lập Trong đó, hoạt độ phóng xạ voxel biến số; tốc độ đếm phép đo, hiệu suất ghi nhận voxel hệ số phương trình Sau đó, hệ phương trình chuẩn hóa dạng ma trận phương pháp Tikhonov sử dụng để đưa lời giải Từ đó, hoạt độ phóng xạ tổng thùng thải xác định Phương pháp ứng dụng để phân tích thùng thải lấp đầy xi măng Portland (mật độ 2,1 g/cm3) chứa nguồn Eu dạng dây Kết cho thấy độ sai biệt tương đối hoạt độ 152 thực giá trị tính tốn nhỏ 16% Krings cộng [2] (năm 2012) đề xuất phương pháp để xác định hoạt độ đồng vị phóng xạ bên thùng thải không đồng Trong nghiên cứu này, tính khơng đồng matrix xem phân lớp vật liệu rác thải cấu trúc che chắn bên thùng chứa, mà hình học chúng 31 ảnh có lọc cao so với ảnh không lọc Điều giúp cho việc nhận diện xác định vị trí nguồn phóng xạ ảnh cắt lớp xác Hình 3.5 Ảnh chụp cắt lớp tương ứng với dữ liệu 19 hình chiếu (a) khơng có lọc, (b) có lọc Hình 3.6 Ảnh chụp cắt lớp tương ứng với dữ liệu 37 hình chiếu (a) khơng có lọc, (b) có lọc Như đặc trưng kỹ thuật chiếu ngược, ta quan sát thấy vệt sáng hình xuất ảnh cắt lớp Các vệt sáng tương ứng với hình chiếu liệu đầu vào sử dụng để tái tạo ảnh Vì thế, số lượng vệt sáng ảnh tái tạo từ 37 hình chiếu nhiều so với ảnh tái tạo từ 19 hình 32 chiếu Thực tế vệt sáng làm giảm khả nhận diện nguồn phóng xạ có hoạt độ thấp giảm độ xác xác định vị trí nguồn phóng xạ ảnh Về bản, số lượng hình chiếu sử dụng để tái tạo ảnh đủ nhiều vệt sáng ảnh Đồng thời, việc sử dụng nhiều hình chiếu tái tạo ảnh làm tăng lên độ tương phản nguồn khu vực xung quanh hình 3.5 hình 3.6 3.3 Xác định vị trí nguồn phóng xạ Để xác định vị trí nguồn phóng xạ bên thùng thải, trước tiên ta phải tiến hành đồng hệ trục tọa độ thùng thải với hệ trục tọa độ ảnh tái tạo Có thể thấy gốc tọa độ đặt góc phía bên trái ảnh cắt lớp, gốc tọa độ thùng thải đặt tâm thùng Hình 3.7 cho thấy ảnh cắt lớp tái tạo kỹ thuật chiếu ngược có lọc với 37 hình chiếu chuyển gốc tọa độ tâm thùng Vị trí nguồn phóng xạ xác định từ liệu ảnh cắt lớp Hình 3.7 Ảnh cắt lớp tái tạo kỹ thuật chiếu ngược có lọc với 37 hình chiếu chủn gốc tọa độ tâm thùng Khi đó, vị trí nguồn phóng xạ xác định dựa vị trí pixel có cường độ lớn khu vực xung quanh điểm sáng ảnh cắt lớp Cụ 33 thể, vị trí pixel có cường độ lớn (19, 8) điểm ảnh cường độ cao (7, -19) điểm ảnh cường độ thấp Độ dài pixel 10 mm nên tọa độ nguồn mặt phẳng Oxy (190, 80) mm (70, -190) mm Trong thực tế, tọa độ mặt phẳng Oxy nguồn phóng xạ bên thùng thải (190, 77) nguồn Cs-A (77, -190) nguồn 137 137 Cs-B mô tả mục 2.1 Như vậy, có phù hợp tốt vị trí nguồn phóng xạ xác định phương pháp quét gamma cắt lớp vị trí thực tế nguồn bên thùng thải 3.3 Tổng kết chương Trong chương 3, luận văn trình bày nội dung sau: - Trình bày kết việc tái tạo ảnh cắt lớp kỹ thuật chiếu ngược có lọc chiếu ngược khơng lọc với số lượng hình chiếu khác Kết đạt cho thấy kỹ thuật lọc tín hiệu miền tần số khơng gian giúp làm giảm độ nhịe ảnh gia tăng độ tương phản thùng nền, nguồn khu vực xung quanh Đồng thời, việc sử dụng liệu đầu vào với nhiều hình chiếu để tái tạo ảnh cắt lớp giúp đạt chất lượng ảnh tốt - Vị trí hai nguồn phóng xạ 137Cs xác định từ ảnh cắt lớp cho thấy phù hợp tốt với vị trí thực tế chúng bên thùng thải Tuy nhiên, chương trình tái tạo ảnh xây dựng cịn hạn chế việc nhận diện nguồn phóng xạ có hoạt độ thấp có nhiều nguồn phóng xạ nằm phân đoạn 34 KẾT LUẬN VÀ KIẾN NGHỊ Luận văn thực nội dung nghiên cứu sau: - Thiết lập hệ đo quét gamma cắt lớp cho thùng thải phóng xạ sử dụng đầu dò NaI(Tl) Đồng thời, tiến hành thực nghiệm quét thùng thải hệ đo để thu thập liệu cho việc tái tạo ảnh cắt lớp - Một điểm luận văn tiến hành khảo sát cho trường hợp có hai nguồn phóng xạ 137Cs với hoạt độ khác bố trí hai vị trí khác bên phân đoạn thùng thải Bên cạnh đó, ý tưởng việc chế tạo thùng thải phóng xạ gồm nhiều phân đoạn nhỏ giúp việc bố trí thực nghiệm xác thuận tiện cho việc thay mẫu - Xây dựng chương trình tái tạo ảnh cắt lớp kỹ thuật chiếu ngược có lọc dựa ngơn ngữ lập trình Python 3.7.4 - Đánh giá ảnh hưởng việc lọc tín hiệu số lượng hình chiếu lên chất lượng ảnh cắt lớp tái tạo kỹ thuật chiếu ngược - Nhận diện xác định vị trí nguồn phóng xạ 137Cs dựa ảnh chụp cắt lớp tái tạo từ liệu quét thực nghiệm chương trình xây dựng Kết cho thấy có phù hợp tốt vị trí nguồn phóng xạ xác định phương pháp quét gamma cắt lớp vị trí thực tế nguồn bên thùng thải Tuy nhiên, chương trình tái tạo ảnh xây dựng cịn hạn chế việc nhận diện nguồn phóng xạ có hoạt độ thấp có nhiều nguồn phóng xạ nằm phân đoạn Như vậy, luận văn hoàn thành mục tiêu nghiên cứu đề Tuy nhiên, luận văn số hạn chế định từ kiến nghị công việc nghiên cứu nên thực sau: - Tiếp tục nghiên cứu cải tiến thuật toán tái tạo ảnh cắt lớp kỹ thuật chiếu ngược có lọc để nhận diện tốt nguồn phóng xạ cho trường hợp có nhiều nguồn với hoạt độ khác nằm phân đoạn - Tiến hành khảo sát cho thùng thải có matrix khác để đánh giá ảnh hưởng matrix lên kết dựng ảnh 35 TÀI LIỆU THAM KHẢO [1] D Stanga, D Gurau, “A new approach in gamma ray scanning of rotating drums containing radioactive waste”, Applied Radiation and Isotopes 70, pp 2149-2153, 2012 [2] T Krings, E Mauerhofer, “Reconstruction of the isotope activity content of heterogeneous nuclear waste drums”, Applied Radiation and Isotopes 70, pp 1100-1103, 2012 [3] P Filß, “Relation between the activity of a high density waste drums and its gamma count rate measured with an unshielded Ge detector”, Applied Radiation and Isotopes 46, pp 805-812, 1995 [4] T Krings, C Genreith, E Mauerhofer, M Rossbach, “A numerical method to improve the reconstruction of the activity content in homogeneous radioactive waste drums”, Nuclear Instruments and Methods in Physics Research A 701, pp 262-267, 2013 [5] R Tushar, M.R More, R Jilju, S Amar, “Active and passive CT for waste assay using LaBr3(Ce) detector”, Radiation Physics and Chemistry 130, pp 29-34, 2017 [6] T.Q Dung, “New measuring technique for assay of radioactive materials in waste drums”, Progress in Nuclear Energy 33, pp 403-420, 1998 [7] T.H Anh, T.Q Dung, “Evaluation of performance of gamma tomographic technique for correcting lump effect in radioactive waste assay”, Annals of Nuclear Energy 28, pp 265-273, 2001 [8] T.T Thanh, H.T.K Trang, H.D Chuong, V.H Nguyen, L.B Tran, H.D Tam, C.V Tao, “A prototype of radioactive waste drum monitor by nondestructive assays using gamma spectrometry”, Applied Radiation and Isotopes 109, pp 544-546, 2016 [9] L.P Jerry and M.L Jonathan, Medical Imaging Signals and Systems, Pearson Education, New Jersey USA, 2015 36 [10] Mirion Technologies Inc (2017) 802 scintillation detectors Available: https://www.mirion.com/products/802-scintillation-detectors (đăng nhập vào ngày 30/08/2019) [11] Mirion Technologies Inc (2017) OspreyTM – Universal digital MCA tube base for scintillation spectrometry Available: http://www.mirion.com/products/osprey-universal-digital-mca-Tube-basefor-scintillation-spectrometry (đăng nhập vào ngày 30/08/2019) [12] Eckert & Ziegler Catalogue of industrial sources Available: https://www.ezag.com/fileadmin/ezag/useruploads/pdf/isotope/5_industrial_s ources.pdf (trang 34, đăng nhập vào ngày 30/08/2018) PL1 PHỤ LỤC Phụ lục 1: Code chương trình Phyton phiên 3.7.4 import numpy as np import matplotlib.pyplot as plt from PIL import Image, ImageChops from scipy.fftpack import fft, fftshift, ifft arr = [] inp = open ("scan_1.txt","r") #read line into array for line in inp.readlines(): # add a new sublist arr.append([]) # loop over the elemets, split by whitespace for i in line.split(): # convert to integer and append to the last # element of the list arr[-1].append(float(i)) sinogram = np.array(arr) #sinogram = np.transpose(sinogram) print(sinogram) print(np.size(sinogram)) plt.imshow(sinogram) imgplot = plt.imshow(sinogram) plt.colorbar() plt.show() print('Performing backprojection') from skimage.transform import iradon from skimage.io import imread PL2 from skimage import data_dir from skimage.transform import radon, rescale theta = np.linspace(0., 360., min(sinogram.shape), endpoint=False) print('Gia tri theta: ') print(len(theta)) reconstruction_fbp = iradon(sinogram, theta=theta, filter="ramp", interpolation="linear", circle=True) imgplot = plt.imshow(reconstruction_fbp) plt.colorbar() #plt.colorbar() #plt.show() #error = reconstruction_fbp - image #print('FBP rms reconstruction error: %.3g' % np.sqrt(np.mean(error**2))) #imkwargs = dict(vmin=-0.2, vmax=0.2) fig, (ax1,ax2) = plt.subplots(1, 2, figsize=(8, 4.5), sharex=True, sharey=True) #ax1.set_title("Reconstruction\nFiltered back projection") ax1.imshow(reconstruction_fbp) #ax2.set_title("Reconstruction error\nFiltered back projection") #ax2.imshow(reconstruction_fbp - image, cmap=plt.cm.Greys_r, **imkwargs) plt.show() PL3 Bảng P1 Số liệu thực nghiệm ứng với góc quay từ 00 – 1100 thùng thải Tọa độ trục Ox (mm) Góc quay thùng thải (0) 10 20 30 40 50 60 70 80 90 100 110 4025 2822 1966 1557 1179 255 15 286 22 240 10 3780 3006 1113 1428 613 226 377 328 315 20 2790 2579 1253 2055 410 446 0 518 475 189 30 2384 3217 1656 536 951 124 204 51 348 932 715 532 40 2847 2579 1095 700 294 494 537 546 892 1721 2574 50 2377 1805 860 1161 96 340 30 963 2147 4947 7194 60 2056 2011 1217 1055 553 474 435 414 1135 5320 9360 14129 70 1756 1652 828 554 496 134 67 497 3288 9615 16443 22754 80 1972 1502 605 891 222 317 407 7092 16054 23271 30743 90 1603 1364 656 944 154 303 386 1063 11703 20528 26138 29667 100 1322 1364 910 296 269 1114 3696 16378 23352 21192 20443 110 1014 1351 51 778 53 277 2154 6665 20419 19766 14415 12062 120 547 837 386 551 368 252 4987 10485 20627 13960 8498 5520 130 696 823 338 345 276 1249 7919 14573 16561 9034 3747 2597 140 736 1005 29 526 523 2453 11787 19744 11161 4526 1690 1552 150 559 794 460 590 93 5057 15804 16279 7198 2298 1368 1944 160 786 592 281 685 1005 8694 18439 11551 3072 1560 826 1880 170 591 417 317 47 2486 12262 17226 7666 1406 1445 1159 1857 180 481 247 226 644 4485 15079 13565 4507 984 1331 1050 2190 190 201 326 283 959 7124 17758 9179 1664 835 843 1577 2772 200 411 201 68 1913 10373 16304 5919 1304 877 1002 1101 3398 210 318 275 4032 14584 12765 3385 900 1105 1022 1571 5104 220 256 219 809 6605 17648 8802 1528 905 741 1157 2165 9040 230 51 543 1378 9644 16764 5721 788 1003 1026 1431 2480 17265 240 351 84 2901 12440 12742 3769 1036 1137 971 1278 3582 35438 250 689 5301 16087 9704 1704 975 993 1253 1965 5654 60180 260 217 983 8201 16959 5797 824 762 954 1045 1927 10684 86666 270 132 1805 10468 14375 3825 1249 931 1038 1155 2736 22690 115467 PL4 280 28 3886 14314 10764 1870 958 1113 1268 1405 4499 44778 145392 290 752 6672 16977 8067 1139 1227 555 1505 1759 7765 69371 158079 300 1601 9532 15379 5114 514 1090 1374 1535 2552 14919 98990 135401 310 2771 12569 13014 2759 415 932 1181 1478 3245 32683 126977 106603 320 5493 16056 9027 1530 318 1360 1431 2132 6215 56300 155788 76476 330 8809 16619 6719 1655 931 1448 1431 2611 10729 85543 152769 50776 340 11579 14572 3559 1376 751 1624 1642 3697 24554 115110 125684 28067 350 14830 11795 2501 1802 1077 1492 2501 5438 48004 146227 96701 13735 360 17382 9014 1173 1434 1404 1785 3125 10241 76356 162844 67625 7100 370 15800 5917 1370 1807 1857 2498 3915 23693 108944 145191 41539 5311 380 12932 4285 1935 1670 1827 2970 5990 47562 141538 112275 20828 3357 390 10886 3085 1856 2183 2197 3653 12107 78639 168278 82109 10313 2644 400 7949 3034 2011 2491 2234 5453 28708 114601 158941 54760 6076 1897 410 5583 3457 2819 2896 3709 10320 58403 150767 125644 30463 3953 1712 420 6170 4218 3496 3798 5953 24829 94737 178589 92072 14305 3170 1958 430 7238 5457 4389 5361 11965 54429 136075 164362 61588 7954 2292 1549 440 12848 8367 6578 10682 29518 95387 174552 129845 34202 5270 1967 1277 450 37030 17911 15121 27699 68801 141888 188905 93237 18927 3463 1400 1327 460 129738 54471 44375 70538 119415 186187 158587 59065 8952 2692 1285 560 470 270825 140097 107479 129383 172856 204818 117604 31584 5295 2576 836 1131 480 414763 248860 187363 196721 222673 171727 77714 14846 3460 2257 1025 610 490 407796 348444 272885 252404 212477 124862 47571 7836 3878 1686 1011 555 500 262831 333460 290587 238375 157283 79221 18755 5042 2510 1602 683 317 510 120578 220887 220346 171305 107073 41826 8749 4246 2012 1145 394 405 520 34983 116945 135283 106130 58382 18227 6279 2671 1965 1302 301 530 12398 40642 64099 53121 24860 9224 4446 2404 1467 988 251 494 540 5306 16542 21920 20734 12287 6546 3641 2261 884 894 134 550 2708 7350 9832 11043 7261 4195 1978 1594 437 791 207 560 1880 4290 3795 5973 3862 3351 1209 893 468 86 189 22 570 2188 2302 2680 3785 2003 2235 640 449 176 127 484 258 580 1632 2747 1947 1341 821 1357 200 373 197 404 76 PL5 Bảng P2 Số liệu thực nghiệm ứng với góc quay từ 1200 – 2300 thùng thải Tọa độ trục Ox (mm) Góc quay thùng thải (0) 120 130 140 150 160 170 180 190 200 210 220 230 0 180 832 376 1068 1325 2931 2072 2181 2622 2109 10 338 440 862 1499 2688 4535 3024 3607 2355 2589 20 402 1201 1573 2840 3782 7422 6362 5741 3308 3369 30 523 437 755 1392 1892 4882 8899 16026 15460 12732 3492 3423 40 854 1032 789 2080 2550 6933 16904 40199 48040 33760 8765 3774 50 5531 2741 1437 2844 4538 12701 43497 88470 111334 99740 30236 6609 60 13415 7238 2238 3064 6152 29411 86988 150059 193278 199386 110566 17119 70 24740 18637 4213 4318 10678 60967 139502 217637 274060 307264 242907 56536 80 34915 33135 14000 6327 23967 103563 194095 256029 288302 339824 383834 197139 90 35795 44855 32151 11509 56507 149500 227350 222246 212032 260598 414081 398494 100 25466 40481 53386 29607 92849 194994 197987 151992 129113 150281 281017 518472 110 14553 25366 55944 81155 132749 205361 142278 89372 58545 62684 138849 414814 120 6087 12240 40349 135699 179692 163351 89001 37950 19736 20833 40969 212362 130 3163 5855 35012 167774 219140 116916 44089 15227 10022 10397 15918 62893 140 2538 6275 46082 175938 227447 64184 16964 8609 6273 6313 9262 20623 150 2746 9062 72827 183089 194895 38509 8975 5818 4024 5207 5775 11600 160 2491 18955 103126 173000 117203 56060 5809 4622 2943 3711 3440 6894 170 4457 37875 137237 137038 46413 96704 5726 3186 2168 2964 2638 5198 180 6752 62630 165270 100464 22829 93603 8635 3191 1953 2240 2063 3867 190 12312 92596 162436 66040 11069 44406 41374 1670 1418 2040 1049 3332 200 26737 122770 129671 36197 6084 10807 100351 3025 692 1690 72 1759 210 49079 153570 96484 17039 4668 5246 103143 5233 921 1063 189 971 220 75346 161322 65520 8071 3101 3457 45813 23061 233 682 449 230 104928 136277 38038 5551 2548 2514 8438 87307 940 428 226 609 240 134450 105985 17601 4116 2129 1820 2816 117022 1592 143 134 250 158437 73969 8407 2884 1547 1662 1719 69449 7761 1028 0 260 148185 46496 5952 2246 1722 1097 1168 13824 48946 984 245 PL6 270 119023 24989 3473 1692 843 1179 643 3443 94427 920 50 280 89096 11991 2978 1582 1151 767 808 1671 87890 3341 0 290 61056 6734 2386 1219 1150 695 1132 31559 23344 148 300 35940 4836 1676 1201 779 278 0 5176 90566 1052 310 17482 3682 1825 741 754 141 233 2093 121232 2213 291 320 8100 2467 1349 1205 152 405 0 495 65824 9939 751 330 4998 2253 1330 789 374 220 0 482 11329 57324 521 340 3437 1416 1224 894 87 136 0 2448 111969 803 350 2884 1369 704 690 329 0 121 1331 93807 2126 360 2182 1331 256 749 55 0 96 462 30540 38731 370 2043 980 296 546 85 75 0 447 5436 94770 380 1391 495 172 283 0 147 102 544 1833 100431 390 856 919 222 0 120 0 45571 400 1253 987 477 140 0 194 386 290 9421 410 716 673 0 42 50 167 605 487 3483 420 752 626 173 270 0 442 435 171 1376 430 518 459 316 0 49 28 52 716 1520 440 844 665 135 0 73 217 84 107 554 623 1660 450 634 641 145 148 45 438 186 660 769 1725 460 759 454 86 93 528 488 250 313 749 838 1083 2300 470 546 347 138 279 157 545 188 703 906 1075 2953 480 259 312 634 38 99 98 67 270 867 1522 2657 490 185 350 194 681 293 473 293 515 854 1362 2916 500 146 223 384 670 830 143 1600 1955 2454 510 93 129 149 335 712 338 881 1197 1461 1375 3420 520 69 467 157 408 255 765 562 1237 1157 2182 2431 3313 530 135 201 138 509 37 634 823 1823 1850 1981 3212 540 198 403 350 88 344 556 876 1300 1146 2227 2446 4524 550 262 508 344 478 301 677 1615 1200 1965 2309 4470 560 287 357 297 235 367 949 980 1825 1360 2131 2628 5156 570 125 164 97 874 278 1052 408 1569 1553 2493 3167 4577 580 212 433 327 642 221 1093 830 1617 1148 3334 2640 4754 PL7 Bảng P3 Số liệu thực nghiệm ứng với góc quay từ 2400 – 3500 thùng thải Tọa độ trục Ox (mm) Góc quay thùng thải (0) 240 250 260 270 280 290 300 310 320 330 340 350 872 877 420 1043 3134 5551 6851 9291 9816 7468 6215 6779 10 839 1014 771 479 1838 5494 6975 7199 8365 7037 6702 5587 20 1363 1400 746 298 1147 4473 6614 8175 8154 7310 5448 4743 30 1539 1119 300 462 738 3114 4786 7225 7351 6376 5308 4682 40 2045 1519 521 1158 3819 5350 6414 5841 6194 5386 3804 50 3069 1745 1019 468 185 2063 3283 5379 5493 5726 4256 3847 60 1731 1879 682 297 375 1700 2268 5490 4465 4909 4454 3912 70 2980 2335 541 14 66 1469 2955 4629 4062 3518 3833 3399 80 12072 3438 1336 334 123 488 2175 3854 3395 3293 4191 2636 90 40845 2973 1038 145 313 1804 2164 2990 3092 3432 2557 100 198352 6263 2324 820 219 848 2872 3085 3498 2713 2201 110 467613 19533 2318 1274 117 472 1282 3511 2518 1968 2259 120 637315 80499 2283 2065 303 339 146 1031 2907 2492 2393 1793 130 491897 363940 4270 2256 400 523 690 2313 1639 1999 1263 140 215941 743734 16468 3012 389 593 1485 1557 1722 1504 150 47855 760160 82884 2349 871 577 358 550 574 1150 1562 682 160 18791 380620 436029 5546 1097 203 357 338 708 1091 1025 170 10714 89673 915512 14737 1663 438 194 399 819 1015 1353 180 5886 27448 848151 47971 1236 45 134 18 548 89 190 3212 14332 338880 322897 864 1184 153 386 48 218 527 205 200 727 9082 56991 951787 4879 1506 96 360 245 223 304 210 46 4863 21098 1078063 22014 1135 46 79 23 382 456 220 2127 10283 491179 114830 533 1054 117 153 373 527 230 422 3870 80924 684559 4521 730 263 0 176 240 175 25 890 27348 1191402 14782 286 728 127 27 250 382 553 13631 838740 40290 456 848 341 200 388 PL8 260 420 310 187 5387 176925 266426 1908 449 0 48 102 270 465 2614 36962 1021121 6732 48 526 124 236 280 127 403 79 378 13148 1248974 19723 411 802 76 20 290 0 186 1069 4336 562561 69616 410 894 185 309 363 300 96 48 1716 1741 80749 517074 4319 230 437 194 164 310 0 931 988 27109 1216638 13539 530 289 320 159 24 603 913 11623 1042040 32430 1693 223 1020 330 21 244 177 3577 295630 197616 4670 184 547 1011 340 359 406 263 259 1661 49532 824248 9562 18 567 350 198 391 0 412 2780 15607 1143745 23970 1268 678 1705 360 988 386 116 439 1602 3723 663379 102570 5319 1427 4922 370 2069 542 294 468 124 756 1660 123428 530910 11652 2058 8468 380 7655 1515 237 610 323 685 2757 30244 1005147 25262 6821 13682 390 39258 1722 584 342 371 550 798 10689 837944 106806 12914 18002 400 81212 4007 947 368 592 417 533 3462 283952 455702 22332 22297 410 85328 18159 1366 1374 453 1246 2944 48906 857453 52677 24656 420 43530 48064 4066 1568 673 498 251 2815 16369 769265 215959 26989 430 11066 72750 15081 3220 1270 561 494 3061 9379 334181 528342 47041 440 3599 64175 36824 7159 2611 1299 1523 3793 12345 78210 724622 161958 450 1557 31228 60446 20682 6279 3714 4290 8546 17136 39721 544343 380419 460 2520 9084 58353 37106 16439 10373 10779 15104 22671 25819 224740 574047 470 2335 5051 33551 50239 29605 20118 18246 20814 23534 19623 50627 503342 480 2497 4813 14422 44039 41375 31080 25893 25047 20518 12761 15619 265579 490 2698 3397 6118 27247 39494 35186 29319 23919 14875 8180 6127 79488 500 2696 4608 5534 13258 26354 27890 23773 16395 9038 4201 4147 23039 510 3126 4793 4585 7563 13367 19121 15947 10413 4281 2069 2868 8342 520 3980 5659 5106 6196 6813 10440 8481 4858 2263 1160 2605 4712 530 4043 6070 4921 6042 5665 5050 3939 1739 1128 1042 1711 3260 540 3647 5586 5800 5842 5205 4907 2429 962 44 565 1030 2079 550 4471 6596 6856 7505 6642 4571 1398 911 419 297 1431 560 5270 6726 6629 7683 6151 5557 2266 950 518 532 1158 570 4593 7929 6810 8085 6521 5435 3406 907 434 98 775 PL9 580 4513 8226 6990 9020 6981 6850 3339 1469 50 135 124 756 ... xuất để xác định hoạt độ đồng vị phóng xạ bên thùng thải Nhìn chung, kỹ thuật chủ yếu dựa hai phương pháp quét gamma phân đoạn quét gamma cắt lớp Trong đó, kỹ thuật đo dựa phương pháp quét gamma. .. Oxy nguồn phóng xạ bên thùng thải (190, 77) nguồn Cs-A (77, -190) nguồn 137 137 Cs-B mô tả mục 2.1 Như vậy, có phù hợp tốt vị trí nguồn phóng xạ xác định phương pháp quét gamma cắt lớp vị trí. .. hợp thùng thải có matrix nguồn phóng xạ phân bố khơng đồng Bên cạnh đó, phương pháp quét gamma cắt lớp [8] chứng minh cách thức hiệu để xác định phân bố nguồn phóng xạ matrix bên thùng thải Phương