(Luận văn thạc sĩ) xác định hoạt độ mẫu môi trường bằng phương pháp trực tiếp

1 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM TP HỒ CHÍ MINH KHOA VẬT LÝ PHẠM THỊ THÙY DƯƠNG KHÓA LUẬN TỐT NGHIỆP XÁC ĐỊNH HOẠT ĐỘ MẪU MÔI TRƯỜNG BẰNG PHƯƠNG PHÁP TRỰC TIẾP Chuyên ngành: VẬT LÍ HẠT NHÂN NĂM 2019 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM TP HỒ CHÍ MINH KHOA VẬT LÝ XÁC ĐỊNH HOẠT ĐỘ MẪU MÔI TRƯỜNG BẰNG PHƯƠNG PHÁP TRỰC TIẾP Người thực hiện: PHẠM THỊ THÙY DƯƠNG Người hướng dẫn khoa học:ThS LÊ QUANG VƯƠNG LỜI CẢM ƠN Trong suốt trình học tập thực khóa luận này, sinh viên nhận nhiều giúp đỡ từ quý Thầy/Cơ, Giờ đây, khóa luận hồn thành, sinh viên xin gửi lời cảm ơn chân thành đến:  ThS Lê Quang Vương – người hướng dẫn khoa học, tận tình bảo tạo điều kiện thuận lợi để sinh viên hồn thành khóa luận  Quý Thầy-Cô Bộ môn Vật lý Hạt nhân, Khoa Vật lý, Trường Đại học Sư Phạm Tp HCM giảng dạy suốt năm qua Những kiến thức mà sinh viên thu nhận qua giảng, môn học Thầy-Cô tảng để sinh viên tiếp thu giải vấn đề khóa luận Hơn nữa, kiến thức cịn hành trang cho sinh viên tiếp công việc sau  Quý Thầy-Cô Hội đồng bảo vệ khóa luận đóng góp ý kiến để khóa luận hồn thiện  Phịng thí nghiệm Vật lý hạt nhân Trường Đại học Sư Phạm Tp.HCM đáp ứng điều kiện sở vật chất, trang thiết bị cần thiết để học viên thực khóa luận  Viện y tế cơng cộng Thành phố Hồ Chí Minh hỗ trợ đo đạc mẫu phân tích hệ phổ kế gamma  Các thành viên gia đình bên động viên chia sẻ lúc khó khăn để an tâm có thêm động lực thực khóa luận TP Hồ Chí Minh, ngày 07 tháng 05 năm 2019 Sinh viên Phạm Thị Thùy Dương Mục lục LỜI CẢM ƠN DANH MỤC CÁC KÍ HIỆU VÀ CHỮ VIẾT TẮT DANH MỤC CÁC BẢNG DANH MỤC CÁC HÌNH VẼ VÀ ĐỒ THỊ MỞ ĐẦU CHƯƠNG I-TỔNG QUAN 1.1 Tình hình nghiên cứu vấn đề liên quan đến khóa luận 1.1.1 Tình hình nghiên cứu giới 1.1.2 Tình hình nghiên cứu Việt Nam 10 1.1.3 Mục tiêu nội dung khóa luận 13 1.2 Cơ sở lí thuyết 13 1.2.1 Chuỗi phóng xạ tự nhiên 13 1.2.2 Qui luật phân rã phóng xạ 16 1.2.3 Cân phóng xạ 17 1.2.4 Tương tác xạ gamma với vật chất 18 1.2.5 Hiệu suất đỉnh lượng toàn phần 20 1.2.6 Hệ số tự hấp thụ 21 1.3 Kết luận chương 23 CHƯƠNG II-ĐỐI TƯỢNG VÀ PHƯƠNG PHÁP NGHIÊN CỨU 24 2.1 Đối tượng nghiên cứu 24 2.1.1 Hệ phổ kế gamma sử dụng đầu dò HPGe 24 2.1.2 Mẫu chuẩn mẫu phân tích 24 2.2 Phương pháp nghiên cứu 26 2.2.1 Các phương pháp xác định hoạt độ 26 2.2.2 Công thức xác định hoạt độ 27 2.3 Kết luận chương 28 CHƯƠNG III-KẾT QUẢ VÀ THẢO LUẬN 29 3.1 Xây dựng đường cong hiệu suất đỉnh lượng toàn phần 29 3.2 Xác định hoạt độ mẫu phân tích 32 3.3 Kết luận chương 37 KẾT LUẬN VÀ KIẾN NGHỊ 38 PHỤ LỤC 42 DANH MỤC CÁC KÍ HIỆU VÀ CHỮ VIẾT TẮT  Các ký hiệu: Ký hiệu Ý nghĩa λ Hằng số phân rã nguồn đồng vị phóng xạ (1/giây) ρ Mật độ khối lượng (g/cm3) ε Hiệu suất đỉnh lượng toàn phần I Xác suất phát xạ gamma nguồn đồng vị phóng xạ A Hoạt độ riêng đồng vị phóng xạ mẫu (Bq/kg)  Các chữ viết tắt: Chữ viết tắt Tiếng Anh Tiếng Việt ETNA Efficiency Transfer for Nuclide Chương trình chuyển hiệu suất Activity measurements cho phép đo phóng xạ FWHM Full Width at Half Maximum Bề rộng nửa giá trị cực đại HPGe High Purity Germanium Đầu dò bán dẫn Ge siêu tinh khiết IAEA International Atomic Energy Cơ quan lượng nguyên tử Agency quốc tế Monte Carlo N Particles Chương trình mơ Monte MCNP Carlo MDA Minimum Detectable Activity Hoạt độ phát tối thiểu WDXRF Wavelength Dispersive X-Ray Phương pháp tán sắc bước sóng Fluorescence huỳnh quang tia X DANH MỤC CÁC BẢNG STT Chỉ số bảng Nội dung Trang 2.1 Thông số mẫu chuẩn RGU-1 mẫu 25 phân tích 3.1 Kết xác định hiệu suất ghi đỉnh 28 lượng toàn phần mẫu RGU-1 3.2 Độ sai biệt hiệu suất ghi đỉnh lượng 30 toàn phần mẫu chuẩn RGU 3.3 Thành phần nguyên tố hóa học có 32 mẫu chuẩn mẫu phân tích 3.4 Hệ số tự hấp thụ mẫu 33 3.5 So sánh hoạt độ trực tiếp hoạt độ cân 34 mẫu 3.6 Độ sai biệt đồng vị 226Ra 35 đồng vị lại đo trực tiếp mẫu 3.7 Độ sai biệt đồng vị 226Ra 35 đồng vị lại đo cân mẫu 12 3.8 So sánh hoạt độ 226Ra xác định phương pháp trực tiếp phương pháp đo cân hoạt độ 35 DANH MỤC CÁC HÌNH VẼ VÀ ĐỒ THỊ STT Chỉ số hình Nội dung Trang 1.1 Chuỗi phóng xạ tự nhiên Uranium-238 14 1.2 Chuỗi phóng xạ tự nhiên Uranium-235 15 1.3 Hiệu ứng quang điện 18 1.4 Tán xạ Compton 19 1.5 Sự tạo cặp 19 1.6 Quá trình hấp thụ xạ qua vật 21 liệu 2.1 Hệ phổ kế gamma phông thấp Viện y 24 tế công cộng Tp HCM 2.2 Mẫu chuẩn RGU-1 (a) mẫu đá phân 25 tích (b) 3.1 Phổ mẫu chuẩn RGU-1 ghi nhận 28 đầu dò HPGe 10 3.2 Đường chuẩn hiệu suất ghi đỉnh 30 lượng toàn phần mẫu chuẩn RGU-1 11 3.3 Phổ mẫu phân tích M1 ghi nhận 31 đầu dò HPGe GEM P4-83 phương pháp trực tiếp 12 3.4 Phổ mẫu phân tích M1 ghi nhận 32 đầu dị HPGe GEM P4-83 phương pháp cân hoạt độ 14 3.5 Đồ thị so sánh hoạt độ đo trực tiếp 34 hoạt độ đo cân mẫu số 15 3.6 So sánh hoạt độ đồng vị 226Ra mẫu phân tích 36 MỞ ĐẦU Xác định hoạt độ phương pháp đo cân phương pháp sử dụng phổ biến Tuy nhiên, với phương pháp cần phải xử lí mẫu cẩn thận có thời gian nhốt mẫu khoảng 20 đến 30 ngày để trình cân xảy Dẫn đến nhiều nhược điểm mẫu sinh học dễ hư hỏng, mẫu đá xảy tượng thoát khí Rn mẫu khơng đóng kín dẫn đến sai số q trình phân tích Từ nhược điểm phương pháp xác định trực tiếp hoạt độ đề xuất với khả khắc phục nhược điểm phương pháp đo cân Cùng với đầu dò bán dẫn Ge siêu tinh khiết (HPGe) với ưu điểm độ phân giải lượng cao sử dụng phổ biến phép phân tích hoạt độ đồng vị phóng xạ phát gamma Khóa luận thực với mục tiêu xác định hoạt độ số mẫu môi trường phổ gamma ghi nhận đầu dò HPGe phương pháp trực tiếp Để xác định hoạt độ đồng vị phóng xạ hệ phổ kế gamma phương pháp đo trực tiếp, có hai vấn đề chủ yếu cần phải giải là: (1) chuẩn hiệu suất ghi nhận đỉnh lượng tồn phần, (2) tính tốn hệ số hấp thụ vật liệu Hiệu suất đỉnh phụ thuộc vào nhiều yếu tố như: đặc trưng phân rã đồng vị phóng xạ diện mẫu, lượng xạ gamma tới, hình học đầu dị, khoảng cách từ mẫu đến đầu dị, kích thước thành phần hóa học mẫu Hệ số hấp thụ phụ thuộc vào thành phần hóa học mẫu Do đó, việc hiệu chỉnh cho vấn đề cần phải thực phép đo cụ thể, khó khăn gặp phải tiến hành phương pháp trực tiếp Khóa luận tập trung trình bày phương pháp xác định hoạt độ số mẫu môi trường (ở mẫu đá qua xử lí) hệ phổ kế gamma sử dụng đầu dị HPGe Trong q trình xác định, nhận thấy có đóng góp tỷ lệ 226 Ra hệ gamma Sử dụng kết từ nghiên cứu khác tổng số đếm đỉnh 186,21 keV có 57,2% 226 Ra với tồn trạng thái cân [9],[11] Các tính tốn sử dụng để xác định hoạt độ 226 Ra cho phép đo trực tiếp nhanh chóng cách sử dụng đỉnh 186 keV Nội dung khóa luận bao gồm:  Trình bày tổng quan tình hình nghiên cứu chủ đề liên quan đến khóa luận; mục tiêu nội dung nghiên cứu; sở lý thuyết sử dụng khóa luận  Mô tả đối tượng nghiên cứu phương pháp sử dụng khóa luận  Trình bày kết xác định hoạt độ mẫu phân tích sau tính tốn hệ số hấp thụ vật liệu, so sánh phương pháp đánh giá độ sai biệt CHƯƠNG I TỔNG QUAN 1.1 Tình hình nghiên cứu vấn đề liên quan đến khóa luận 1.1.1 Tình hình nghiên cứu giới Hiệu suất ghi đỉnh lượng toàn phần đại lượng quan trọng, ảnh hưởng đến kết phân tích hoạt độ đồng vị phóng xạ hệ phổ kế gamma Trong phép đo mẫu mơi trường, hàm lượng đồng vị phóng xạ diện mẫu thường thấp, nên yêu cầu mẫu đo phải chuẩn bị dạng thể tích với kích thước lớn Khi đó, hiệu suất đỉnh phép đo bị ảnh hưởng hiệu ứng tự hấp thụ mẫu Trong điều kiện hình học đo, hiệu ứng tự hấp thụ phụ thuộc vào thành phần nguyên tố mật độ khối lượng mẫu, mà khác biệt mẫu đo khác Do đó, việc tính tốn hiệu suất đỉnh hệ số hiệu chỉnh tự hấp thụ cần phải thực mẫu đo Đường chuẩn hiệu suất đỉnh cho phép đo quan tâm xác định thực nghiệm nghiên cứu Vukanac cộng [19] Phương pháp địi hỏi phải có mẫu tham khảo với hình học, thành phần nguyên tố mật độ khối lượng tương tự mẫu phân tích Khi đó, hoạt độ đồng vị phóng xạ mẫu phân tích xác định dựa đường chuẩn hiệu suất thực nghiệm Tuy nhiên, phương pháp thực nghiệm khó ứng dụng nhiều phịng thí nghiệm, khơng đủ điều kiện để chế tạo mẫu tham khảo tốn thời gian tiền bạc Do đó, phương pháp tính tốn phát triển để giải vấn đề Một số chương trình mô đxa phát triển sử dụng rộng rãi MCNP, GEANT, EGS chương trình tính toán ETNA, EFTRAN, TRUECOIN Độ tin cậy chương trình tính tốn đánh giá nhiều cơng trình nghiên cứu Trong số nghiên cứu tiêu biểu sau: Radu cộng [18] đánh giá khả phương pháp chuyển hiệu suất, sử dụng chương trình ETNA, để tính tốn hiệu suất đỉnh cho vài hình học đo khác dựa đường chuẩn hiệu suất đỉnh thu từ hình học đo tham khảo nguồn dạng điểm Trong nghiên cứu này, hệ phổ kế gamma sử dụng đầu dò HPGe đồng trục loại n sử dụng để ghi nhận phổ thực nghiệm Các nguồn phóng xạ chuẩn 33 Sau thời gian nhốt mẫu 30 ngày để đạt trạng thái cân phóng xạ đồng vị 238 U sản phẩm cháu (214Pb, 214Bi), mẫu đo với thời gian 86400s Hình 3.4 Phổ gamma mẫu M1 phông đo cân Mẫu chuẩn sử dụng có thành phần hóa học xác định [22], bên cạnh nhận thấy mẫu phân tích có thành phần hóa học giống với mẫu nghiên cứu Kazem cộng [16] Nồng độ nguyên tố mẫu phân tích (bảng 3.3) tính theo trung bình kết [16] sử dụng để xác định hệ số suy giảm tuyến tính Bảng 3.3 Thành phần ngun tố hóa học có mẫu chuẩn mẫu phân tích Nồng độ nguyên tố STT Nguyên tố Z Mẫu chuẩn RGU-1 Mẫu phân tích [22] [16] O 0,533399 0,479129 Al 13 0,001 0,130915 Si 14 0,464 0,274131 34 S 16 0,00002 K 19 0,00002 0,003955 Ca 20 0,0003 0,057517 Ti 22 0,00008 0,000192 Fe 26 0,0003 0,001996 12 Pb 82 0,00008 13 U 92 0,0004 14 Na 11 0,0002 0,051790 15 C 0,0001 16 Mg 12 0,0001 0,00001 17 Th 90 0,000001 18 Mn 14 0,000240 19 Ba 56 0,000117301 Khi xác định thành phần nguyên tố hóa học có mẫu ta sử dụng phần mềm XCOM [21] để xác định hệ số suy giảm khối đỉnh lượng Dựa vào mật độ khối lượng mẫu ta tính hệ số suy giảm tuyến tính mẫu đỉnh lượng có mẫu Sử dụng cơng thức (1.17) ta có hệ số tự hấp thụ mẫu bảng 3.4 Bảng 3.4 Hệ số tự hấp thụ mẫu 63,3 186,2 351,9 609,3 1764 keV keV keV keV keV M1 0,98 1,01 1,01 1,01 1,00 M2 0,98 1,00 1,01 1,01 1,00 M3 0,98 1,01 1,01 1,01 1,00 M4 0,97 1,00 1,01 1,01 1,00 M5 0,96 0,99 1,00 1,00 1,00 Mẫu Sau xác định hệ số tự hấp thụ mẫu sử dụng công thức (2.1) xác định định hoạt độ sau tính toán hệ số tự hấp thụ mẫu Bảng 3.5 So sánh hoạt độ trực tiếp hoạt độ cân mẫu M1 35 Mẫu M1 Đồng vị Năng Hoạt độ Hoạt độ lượng trực tiếp cân (keV) (Bq.kg 1 ) (Bq.kg 1 ) Độ sai biệt (%) 234 Th 63,30 55,14 ± 3,12 57,56 ± 3,21 4,22 226 Ra 186,21 62,77 ± 1,32 68,73 ± 1,39 8,67 351,90 70,24 ± 0,77 70,88 ± 0,78 0,91 214 Pb 214 Bi 609,30 70,00 ± 0,90 69,71 ± 0,90 0,41 214 Bi 1764,30 63,12 ± 1,61 67,54 ± 1,67 6,54 Hình 3.5 Đồ thị so sánh hoạt độ đo trực tiếp hoạt độ đo cân mẫu M1 Hoạt độ xác định có giá trị từ (55,14 ± 3,12) Bq.kg-1 (70,24 ± 0,77) Bq.kg-1 đo trực tiếp (57,56 ± 3,21) Bq.kg-1 (70,88 ± 0,78) Bq.kg-1 đo cân với độ sai biệt nhỏ 9% Bảng 3.6 Độ sai biệt tương đối 226Ra đồng vị lại đo trực tiếp mẫu 36 Hoạt độ 226Ra Hoạt độ đo trực tiếp 186,21 keV (Bq.kg-1) 234 214 214 (%) (Bq.kg-1) Th (63,3 keV) 55,14 ± 3,12 Pb (351,9 keV) 70,24 ± 0,77 214 Độ sai biệt Bi (609,3 keV) 70,00 ± 0,90 Bi (1764,3 keV) 63,12 ± 1,61 13,85 10,62 62,77 ± 1,32 10,32 0,55 Bảng 3.7 Độ sai biệt tương đối đồng vị 226Ra đồng vị lại đo cân mẫu Hoạt độ đo cân (Bq.kg-1) 234 214 214 186,21 keV Độ sai biệt (%) (Bq.kg-1) Th (63,3 keV) 57,56 ± 3,21 Pb (351,9 keV) 70,88 ± 0,78 214 Hoạt độ 226Ra Bi (609,3 keV) 69,71 ± 0,90 Bi (1764,3 keV) 67,54 ± 1,67 16,25 3,13 68,73 ± 1,39 1,42 1,73 Tương tự ta tính hoạt độ mẫu cịn lại khoảng (31,05 ± 2,00) Bq.kg-1 đến (99,26 ± 1,00) Bq.kg-1 đo trực tiếp (57,56 ± 3,21) Bq.kg-1 (106,19 ± 1,87) Bq.kg-1 đo cân hoạt độ (Phụ lục B) Bảng 3.8 Độ sai biệt hoạt độ 226Ra đo trực tiếp đo cân Hoạt độ trực tiếp Hoạt độ cân Độ sai biệt (Bq.kg 1 ) (Bq.kg 1 ) (%) M1 62,77 ± 1,32 68,73 ± 1,39 8,67 M2 65,70 ± 1,38 81,15 ± 1,57 19,03 M3 52,45 ± 1,16 64,70 ± 1,32 18,94 M4 93,19 ± 1,71 106,19 ± 1,87 12,24 M5 61,18 ± 1,20 75,78 ± 1,38 19,28 Mẫu Các giá trị hoạt độ 226Ra trình bày bảng có giá trị từ (65,70 ±1,38) Bq.kg-1 (93,19 ± 1,71) Bq.kg-1 đo trực tiếp từ (64,70 ± 1,32) 37 Bq.kg-1 (106,19 ± 1,87) Bq.kg_1 đo cân với độ sai biệt hai phương pháp 20% Hình 3.6 So sánh hoạt độ đồng vị 226Ra mẫu phân tích 3.3 Kết luận chương Trong chương 3, khóa luận trình bày nội dung: xây dựng đường cơng hiệu suất ghi đỉnh lượng toàn phần mẫu chuẩn RGU-1 với độ sai biệt nhỏ 8% Tính tốn hoạt độ mẫu phân tích có giá trị khoảng từ (55,14 ± 3,12) Bq.kg-1 (70,24 ± 0,77) Bq.kg-1 đo trực tiếp (57,56 ± 3,21) Bq.kg-1 (70,88 ± 0,78) Bq.kg-1 đo cân với độ sai biệt nhỏ 9% Các giá trị hoạt độ 226Ra trình bày bảng có giá trị từ (65,70 ±1,38) Bq.kg-1 (93,19 ± 1,71) Bq.kg-1 đo trực tiếp từ (64,70 ± 1,32) Bq.kg-1 (106,19 ± 1,87) Bq.kg1 đo cân với độ sai biệt hai phương pháp 20% 38 KẾT LUẬN VÀ KIẾN NGHỊ Khóa luận đạt kết sau:  Xây dựng đường cong hiệu suất đỉnh lượng tồn phần với độ sai biệt khơng q 8%  Xác định hệ số tự hấp thụ mẫu  Xác định hoạt độ 226Ra, 234Th, 214Pb, 214Bi trong mẫu phân tích phương pháp trực tiếp phương pháp cân bằng, độ sai biệt hai phương pháp nhỏ 20%  Các kết đạt cho thấy phương pháp đo trực tiếp đáng tin cậy áp dụng để xác định hoạt độ 226Ra sau chuẩn bị mẫu mà không cần thời gian chờ đợi Như vậy, khóa luận hồn thành mục tiêu nghiên cứu đề Đồng thời, để đạt kết tốt việc phân tích hoạt độ đồng vị phóng xạ mẫu mơi trường, luận văn kiến nghị công việc nghiên cứu nên thực sau:  Xác định nồng độ nguyên tố hóa học có mẫu phân tích tính tốn xác hệ số tự hấp thụ  Đánh giá ảnh hưởng hiệu ứng trùng phùng đồng vị phát nhiều lượng gamma 214Pb, 214Bi đến kết xác định hoạt độ  Tiến hành đo đạc nhiều mẫu phân tích với bề dày, thành phần khác để đánh giá độ tin cậy phương pháp 39 TÀI LIỆU THAM KHẢO Tài liệu tiếng Việt [1] Châu Văn Tạo (2013), Vật lí hạt nhận đại cương, Nhà xuất Đại học Quốc gia Thành phố Hồ Chí Minh [2] Hồng Đức Tâm (2019), Phân tích sai số liệu thực nghiệm, Nhà xuất Trường Đại học Sư Phạm Thành phố Hồ Chí Minh [3] Ronald G-William S (1983), Vật lí đại, Nhà xuất Giáo dục Tài liệu tiếng Anh [4] Huy N.Q, Luyen T.V (2004), A method to determine 238U activity in environmental soil samples by using 63.3-keV-photopeak-gamma HPGe spectromete, Applied Radiation and Isotopes, 61, 1419–1424 [5] Huy N.Q (2010), The influence of dead layer thickness increase on efficiency decrease for a coaxial HPGe p-type detector, Nuclear Instruments and Methods in Physics Research A, 621, 390–394 [6] Huy N.Q, Binh D.Q (2014), A semi-empirical approach to analyze the activities of cylindrical radioactive samples using gamma energies from 185 to 1764 keV, Applied Radiation and Isotopes, 94, 82–88 [7] Loan T.T.H, Thanh T.T, Phuong D.N (2007), Gamma spectrum simulation and coincidence summing factor calculation for point sources using MCNP code, Communication in Physics, 2, 110–116 [8] Loan T.T.H, Phuong D.N, Phong D.P.H (2009), Investigating the effect of matrices and densities on the efficiency of HPGe spectroscopy using MCNP, Communications in Physics, 19, 45-52 [9] Thanh.T.T, Ferreux L., Lépy M.C., Tao C.V (2010), Determination activity of radionuclides in marine sediment by gamma spectrometer with anti cosmic shielding, Journal of Environmental Radioactivity, 101, 780-783 40 [10] Thanh T.T, Tao C.V, Loan T.T.H, Au B.H (2012), A detail investigation of interactions within the shielding to HPGe detector response using MCNP code, Kerntechnik, 77, 458-461 [11] Thanh T.T, Vuong L.Q (2018), Validation of an advanced analytical procedure applied to the measurement of environmental radioactivity, Journal of Environmental Radioactivity 184-185, 10 [12] Vuong L.Q, Tao C.V, Thanh T.T (2017), Revision of nuclear data of U-235 and Ra-226 for the 186 keV gamma-ray peak for the determination of activity in environmental samples, J Radioanal Nucl Chem 314:1273-1277 [13] Corte F.De (2005), Direct gamma-spectrometic measurement of the Ra-226 186,2 keV line for detecting U-238/Ra-226 disequilibrium in determining thi environmental dose rate for the luminescence dating of sediments, Aplplied Radiation and Isotopes 63:509-598 [14] Ebaid Y.Y (2009), Use of gamma ray spectrometry for Uranium isotopic analysis in enviromental samples, Physics Department Faculty of Science, Fayoum University, Fayoum, 63514 Egypt [15] Ebaid Y.Y (2005), U-235 -  emission contribution to the 186 keV energy transition of Ra-226 in environmental samples activity calculations, International Congress Seres 1276: 409-411 [16] Kazem K, Ali K, Yilin X, Fatemeh S (2018), Chemical composition of rockforming minerals and crystallization physicochemical conditions of the Middle Eocene I-type Haji Abad pluton, SW Buin-Zahra, Iran, Arabian Journal of Geosciences 11:717 226 [17] Völgyesi P 2014), Using the 186-keV peak for Ra activity concentration determination in Hungarian coal-slag samples by gamma-ray spectroscopy, J Radioanal Nucl Chem 302:375-383 [18] Radu D., Stanga D., Sima O (2009), ETNA software used for efficiency transfer from a point source to other geometries, Applied Radiation and Isotopes,67,1686–1690 41 [19] Vukanac I., Djurasevic M., Kandic A., Novkovic D., Nadjerdj L., Milosevic Z (2008), Experimental determination of the HPGe spectrometer efficiency curve, Applied Radiation and Isotopes, 66 Tài liệu tham khảo website: [20] http://www.nucleide.org/Laraweb/index.php [21] https://physics.nist.gov/PhysRefData/Xcom/html/xcom1.html [22]https://nucleus.iaea.org/rpst/ReferenceProducts/ReferenceMaterials/Radionuclides /IAEA-RGU-1.htm 42 PHỤ LỤC Phụ lục A: Hoạt độ động vị phóng xạ có mẫu phân tích xác định phương pháp trực tiếp phương pháp cân phóng xạ chưa hiệu chỉnh hệ số tự hấp thụ Năng Mẫu Đồng vị lượng (keV) trực tiếp cân Độ sai biệt (%) 1 1 (Bq.kg ) (Bq.kg ) Th 63,30 72,927 ± 3,847 89,065 ± 4,425 18,12 226 Ra 186,21 65,290 ± 1,369 80,639 ± 1,563 19,03 351,90 72,850 ± 0,801 74,349 ± 0,812 2,02 Pb 214 Bi 609,30 70,403 ± 0,913 79,223 ± 0,992 11,13 214 Bi 1764,30 64,794 ± 1,658 72,290 ± 1,762 10,37 234 Th 63,30 47,035 ± 2,793 77,500 ± 3,914 39,31 226 Ra 186,21 52,362 ± 1,158 64,595 ± 1,318 18,94 351,90 56,127 ± 0,652 56,991 ± 0,659 1,52 214 Hoạt độ 234 214 Hoạt độ Pb 214 Bi 609,30 55,359 ± 0,755 57,672 ± 0,776 4,01 214 Bi 1764,30 49,456 ± 1,378 57,720 ± 1,499 14,32 234 Th 63,30 95,617 ± 4,638 77,041 ± 3,980 24,11 226 Ra 186,21 92,674 ± 1,702 105,600 ± 1,856 12,24 351,90 97,892 ± 0,985 101,509 ± 1,011 3,56 214 Pb 214 Bi 609,30 96,094 ± 1,136 99,693 ± 1,168 3,61 214 Bi 1764,30 91,572 ± 2,000 97,243 ± 2,069 5,83 43 234 Th 63,30 56,281 ± 3,187 58,758 ± 3,279 4,22 226 Ra 186,21 62,523 ± 1,312 68,456 ± 1,388 8,67 351,90 69,351 ± 0,764 69,988 ± 0,769 0,91 214 Pb 214 Bi 609,30 69,430 ± 0,893 69,144 ± 0,891 0,41 214 Bi 1764,30 62,940 ± 1,602 67,347 ± 1,663 6,54 234 Th 63,30 32,356 ± 2,089 69,964 ± 3,491 53,75 226 Ra 186,21 61,634 ± 1,213 76,353 ± 1,393 19,28 351,90 68,755 ± 0,719 70,282 ± 0,730 2,17 214 Pb 214 Bi 609,30 70,540 ± 0,856 72,320 ± 0,872 2,46 214 Bi 1764,30 66,746 ± 1,524 70,534 ± 1,572 5,37 Phụ lục B : Hoạt độ động vị phóng xạ có mẫu phân tích xác định phương pháp trực tiếp phương pháp cân phóng xạ hiệu chỉnh hệ số tự hấp thụ Mẫu Đồng vị Năng Hoạt độ Hoạt độ Độ sai lượng trực tiếp cân biệt (keV) (Bq.kg 1 ) (Bq.kg 1 ) (%) 234 Th 63,30 71,74 ± 3,78 87,61 ± 4,35 18,12 226 Ra 186,21 65,70 ± 1,38 81,15 ± 1,57 19,03 351,90 73,91 ± 0,81 75,43 ± 0,82 2,02 214 Pb 214 Bi 609,30 71,08 ± 0,92 79,99 ± 1,00 11,13 214 Bi 1764,30 65,04 ± 1,66 72,56 ± 1,77 10,37 44 Mẫu Đồng vị Mẫu Hoạt độ Hoạt độ Độ sai lượng trực tiếp cân biệt (keV) (Bq.kg 1 ) (Bq.kg 1 ) (%) 234 Th 63,30 45,88 ± 2,72 75,60 ± 3,82 39,31 226 Ra 186,21 52,45 ± 1,16 64,70 ± 1,32 18,94 351,90 56,73 ± 0,66 57,61 ± 0,67 1,52 214 Mẫu Năng Pb 214 Bi 609,30 55,73 ± 0,76 58,06 ± 0,78 4,01 214 Bi 1764,30 49,55 ± 1,38 57,83 ± 1,50 14,32 Năng Hoạt độ Hoạt độ Độ sai lượng trực tiếp cân biệt (keV) (Bq.kg 1 ) (Bq.kg 1 ) (%) Đồng vị 234 Th 63,30 93,94 ± 4,56 75,69 ± 3,91 24,11 226 Ra 186,21 93,19 ± 1,71 106,19 ± 1,87 12,24 351,90 99,26 ± 1,00 102,93 ± 1,03 3,56 214 Pb 214 Bi 609,30 96,98 ± 1,15 100,61 ± 1,18 3,61 214 Bi 1764,30 91,89 ± 2,01 97,58 ± 2,08 5,83 Năng Hoạt độ Hoạt độ lượng trực tiếp cân (keV) (Bq.kg 1 ) (Bq.kg 1 ) Đồng vị Độ sai biệt (%) 234 Th 63,30 31,05 ± 2.00 67.15 ± 3,35 53,75 226 Ra 186,21 61,18 ± 1.20 75.78 ± 1,38 19,28 351,90 69,00 ± 0.72 70.53 ± 0,73 2,17 214 Pb 214 Bi 609,30 70,60 ± 0.86 72.38 ± 0,87 2,46 214 Bi 1764,30 66,65 ± 1.52 70.43 ± 1,57 5,37 45 PHỤ LỤC C: Độ sai biệt đỉnh 226Ra đỉnh lại đo trực tiếp Hoạt độ 226Ra Hoạt độ đo trực tiếp 186,21 keV (Bq.kg-1) 234 214 M2 Th (63,3 keV) 71,74 ± 3,78 Pb (351,9 keV) 73,91 ± 0,81 214 214 1,02 Th (63,3 keV) 45,88 ± 2,72 14,31 Pb (351,9 keV) 56,73 ± 0,66 7,56 55,73 ± 0,76 Bi (1764,3 keV) 49,55 ± 1,38 5,84 Th (63,3 keV) 93,94 ± 4,56 0,79 Pb (351,9 keV) 99,26 ± 1,00 214 93,19 ± 1,71 5,88 6,12 Bi (609,3 keV) 96,98 ± 1,15 Bi (1764,3 keV) 91,89 ± 2,01 1,42 Th (63,3 keV) 31,05 ± 2,00 97,00 Pb (351,9 keV) 69,00 ± 0,72 234 214 214 214 52,45 ± 1,16 7,57 Bi (609,3 keV) 214 M5 11,11 65,04 ± 1,66 214 214 65,70 ± 1,38 Bi (1764,3 keV) 234 M4 8,41 71,08 ± 0,92 214 214 (%) Bi (609,3 keV) 234 M3 (Bq.kg-1) Độ sai biệt Bi (609,3 keV) 70,60 ± 0,86 Bi (1764,3 keV) 66,65 ± 1,52 61,18 ± 1,20 3,90 11,34 13,35 8,21 PHỤ LỤC D: Độ sai biệt đỉnh 226Ra đỉnh lại đo cân Hoạt độ đo trực tiếp (Bq.kg-1) 234 M2 214 87,61 ±4,35 Pb (351,9 keV) 75,43±0,82 Bi (609,3 keV) 186,21 keV (Bq.kg-1) Th (63,3 keV) 214 Hoạt độ 226Ra 79,99±1,00 Độ sai biệt (%) 7,97 81,15±1,57 7,04 1,43 46 214 Bi (1764,3 keV) 72,56±1,77 10,58 Th (63,3 keV) 75,60±3,82 16,85 Pb (351,9 keV) 57,61±0,67 234 214 M3 214 214 58,06±0,78 Bi (1764,3 keV) 57,83±1,50 10,61 Th (63,3 keV) 75,69±3,91 28,73 Pb (351,9 keV) 102,93±1,03 214 214 214 3,07 100,61±1,18 Bi (1764,3 keV) 97,58±2,08 8,11 Th (63,3 keV) 67,15±3,35 11,40 Pb (351,9 keV) 70,53±0,73 214 214 214 106,19±1,87 10,27 Bi (609,3 keV) 234 M5 10,96 Bi (609,3 keV) 234 M4 64,70±1,32 Bi (609,3 keV) 72,38±0,87 Bi (1764,3 keV) 70,43±1,57 75,78±1,38 5,26 6,93 4,49 7,07 47 ... dụng phương pháp xác định hoạt độ trực tiếp hệ phổ kế gamma ghi nhận đầu dò HPGe số mẫu mơi trường sau so sánh với kết đo từ phương pháp xác định cân hoạt độ Từ giải thích xác định hoạt độ số mẫu. .. 2.2 Mẫu chuẩn RGU-1 (a) mẫu đá phân tích (b) 2.2 Phương pháp nghiên cứu 2.2.1 Các phương pháp xác định hoạt độ Phương pháp cân phương pháp xác định hoạt độ phổ gamma thông qua việc xác định hoạt. .. độ đồng vị có mẫu phân tích với phương pháp đo trực tiếp Đồng thời xác định hệ số tự hấp thụ mẫu định xử lí kết hoạt độ với hệ số tự hấp thụ xác định So sánh kết với hoạt độ xác định phương pháp