Nghiên cứu thiết kế và chế tạo hệ thống thiết bị thực nghiệm để đo một số đặc trưng vật lý neutron, phân tích kích hoạt và định liều neutron phục vụ công tác đào tạo nhân lực hạt nhân

BỘ KHOA HỌC VÀ CÔNG NGHỆ VIỆN NĂNG LƯỢNG NGUYÊN TỬ VIỆT NAM BÁO CÁO TỔNG KẾT ĐỀ TÀI KHOA HỌC CÔNG NGHỆ CẤP BỘ NĂM 2010-2012 NGHIÊN CỨU, THIẾT KẾ VÀ CHẾ TẠO HỆ THỐNG THIẾT BỊ THỰC NGHIỆM ĐỂ ĐO MỘT SỐ ĐẶC TRƯNG VẬT LÝ NEUTRON, PHÂN TÍCH KÍCH HOẠT VÀ ĐỊNH LIỀU NEUTRON PHỤC VỤ CÔNG TÁC ĐÀO TẠO NHÂN LỰC HẠT NHÂN (Mã số: 03/10/NLNT) Cơ quan chủ trì: Viện Nghiên cứu hạt nhân Chủ nhiệm: TS NCVC Nguyễn Văn Hùng LÂM ĐỒNG, THÁNG 11/2012 DANH SÁCH NHỮNG NGƯỜI CHÍNH THAM GIA THỰC HIỆN ĐỀ TÀI TT Họ tên Học hàm, học vị, ngạch công chức Đơn vị công tác Nguyễn Văn Hùng TS, NCVC Trung tâm đào tạo, Viện NCHN Nguyễn Minh Tuân CN, NCV Trung tâm LPƯ, Viện NCHN Đặng Lành ThS, NCV Phòng Vật lý – Điện tử HN, Viện NCHN Nguyễn Xuân Hải TS, NCV -nt- Cao Đơng Vũ ThS, NCV Trung tâm phân tích, Viện NCHN Trần Trí Viễn CN, NCV Trung tâm LPƯ, Viện NCHN Trần Quốc Dưỡng CN, NCV -nt- Trang Cao Sử CN, NCV -nt- Phạm Ngọc Sơn ThS, NCV Phòng Vật lý – Điện tử HN, Viện NCHN 10 Trần Quang Thiện CN, NCV Trung tâm phân tích, Viện NCHN 11 Hồng Sỹ Minh Phương CN, NCV Trung tâm đào tạo, Viện NCHN i MỤC LỤC Tiêu đề Những ký hiệu chữ viết tắt Tóm tắt Abstract MỞ ĐẦU Tình hình nghiên cứu nước nước Nhu cầu lý thực Đề tài Mục tiêu đề tài Các nội dung nghiên cứu Đề tài Ý nghĩa khoa học thực tiễn Đề tài Cách tiếp cận phương pháp nghiên cứu sử dụng Đề tài Đơn vị thực Thời gian thực Nguồn mức kinh phí cấp Phần I TỔNG QUAN LÝ THUYẾT Chương Đặc trưng làm chậm khuếch tán neutron môi trường nước 1.1 Một số khái niệm 1.2 Đặc trưng làm chậm neutron môi trường nước 1.3 Đặc trưng khuếch tán neutron môi trường nước 1.4 Chiều dài khuếch tán L 1.5 Chiều dài làm chậm Ls 1.6 Diện tích di cư M2 Chương Phương pháp mơ Monte-Carlo chương trình MCNP 2.1 Phương pháp Monte-Carlo 2.2 Chương trình MCNP 2.3 Input file số lệnh thường dùng MCNP 2.3.1 Tiêu đề 2.3.2 Cell cards (Các thẻ ô) 2.3.3 Surface cards (Các thẻ bề mặt) 2.3.4 Data cards 2.4 Output file chương trình MCNP 2.5 Số liệu hạt nhân MCNP 2.5.1 Số liệu hạt nhân 2.5.2 Số liệu tiết diện tương tác Chương Phương pháp phân tích kích hoạt neutron dùng nguồn đồng vị 3.1 Phương pháp phân tích kích hoạt neutron Trang vi viii viii 1 3 4 4 5 10 11 12 13 13 13 15 15 15 16 17 19 20 20 20 21 21 ii 3.1.1 Nguyên lý phương pháp NAA 3.1.2 Phương trình kích hoạt 3.1.3 Các phương pháp chuẩn hóa xác định hàm lượng nguyên tố NAA 3.2 Xác định thông số phổ neutron phục vụ cho phân tích kích hoạt 3.2.1 Phương pháp xác định hệ số lệch phổ  3.2.2 Thực nghiệm xác định tỉ số thơng lượng f 3.2.3 Tính tốn thơng lượng neutron nhiệt 3.2.4 Tính tốn thơng lượng neutron nhiệt Chương Phương pháp đo neutron dựa tương tác hạt tích điện 4.1 Giới thiệu 4.2 Đo neutron dùng ống đếm Bo dựa phản ứng 10B(n, α)7Li 4.3 Ống đếm tỉ lệ 3He dựa phản ứng 3He(n, p)3H Chương Phương pháp đo suất liều neutron dựa phép đo thông lượng 5.1 Giới thiệu 5.2 Đo thơng lượng neutron dùng dị kích hoạt 5.2.1 Phương trình kích hoạt 5.2.2 Xác định thơng lượng neutron nhiệt nhiệt 5.3 Xác định suất liều neutron từ kết đo thông lượng 21 22 24 26 26 30 30 31 32 32 32 36 38 38 38 38 41 41 Phần II THỰC NGHIỆM VÀ KẾT QUẢ Chương Giới thiệu chung hệ Howitzer chứa nước dùng nguồn neutron đồng vị 6.1 Hệ Howitzer chứa nước 6.2 Nguồn neutron đồng vị sử dụng Đề tài 6.2.1 Nguồn 252Cf 6.2.2 Nguồn 241Am-Be Chương Xây dựng cấu hình hệ Howitzer chứa nước dùng nguồn neutron đồng vị mơ MCNP 7.1 Mơ tả hình học vật liệu 7.1.1 Mơ tả hình học Howitzer chứa nước 7.1.2 Vật liệu sử dụng 7.2 Xây dựng input file 7.2.1 Thiết lập cell card 7.2.2 Thiết lập surface card 7.2.3 Thiết lập Data card 7.3 Xây dựng cấu hình Howitzer chứa nước mơ MCNP 7.3.1 Mơ hình dạng kích thước Howitzer nước dùng nguồn neutron 7.3.2 Các hình ảnh Howitzer nước vẽ chương trình MCNP 7.4 Kết thảo luận 7.4.1 Kết tính tốn đặc trưng làm chậm khuếch tán neutron 7.4.2 Đánh giá an toàn xạ sử dụng Howitzer chứa nước mô MCNP 43 43 43 45 45 46 47 47 47 48 49 49 51 52 53 53 54 56 56 62 iii 7.4.3 Kết luận Chương Thiết kế chế tạo hệ thống thiết bị thực nghiệm đo neutron 8.1 Sửa chữa nâng cấp hệ đo neutron Ấn Độ 8.1.1 Tình trạng hệ ghi đo neutron ghép nối máy tính Ấn Độ 8.1.2 Thiết kế, chế tạo khối nguồn nuôi thấp 8.1.3 Thiết kế, chế tạo khối cao kV 8.1.4 Thiết kế, chế tạo khối khuếch đại phổ kế ngõ tựa Gauss 8.1.5 Thiết kế, chế tạo khối MCA 4k 8.1.6 Tổng quan ngôn ngữ LabVIEW 8.1.7 Kết kiểm tra khối cao kV 8.1.8 Kết kiểm tra khối khuếch đại phổ kế 8.1.9 Kiểm tra độ phi tuyến vi phân tích phân hệ thống 8.1.10 Hình ảnh hệ đo neutron sửa chữa nâng cấp 8.2 Chế tạo hệ thống Howitzer chứa nước 8.2.1 Cấu hình kích thước hệ Howitzer chứa nước 8.2.2 Mơ tả bổ sung vẽ Chương Đo thực nghiệm hệ thống Howitzer chứa nước chế tạo 9.1 Đo phân bố thông lượng suất liều neutron 9.1.1 Thiết bị dụng cụ thực nghiệm 9.1.2 Đo phân bố thông lượng suất liều neutron theo chiều bán kính ngang nguồn 252Cf 9.1.3 Đo phân bố thông lượng suất liều neutron theo chiều thẳng đứng nguồn 252Cf 9.1.4 Đo phân bố thông lượng suất liều neutron theo chiều bán kính ngang nguồn 241Am-Be 9.1.5 Kết đo phân bố thông lượng suất liều neutron theo chiều thẳng đứng nguồn 241Am-Be 9.1.6 Khảo sát suất liều neutron gamma nguồn 252Cf máy đo liều xách tay 9.1.7 Nhận xét 9.2 Đo đặc trưng làm chậm khuếch tán neutron môi trường nước 9.2.1 Chiều dài làm chậm khuếch tán neutron phụ thuộc vào nhiệt độ nguồn 241Am-Be 9.2.2 Chiều dài làm chậm khuếch tán neutron nguồn 252Cf 9.2.3 Nhận xét 9.3 Khảo sát thông số phổ neutron khả xác định hàm lượng số nguyên tố mẫu xi măng cà phê 9.3.1 Giới thiệu 9.3.2 Hệ đo mẫu 9.3.3 Khảo sát thông số phổ neutron 9.3.4 Các bước thực nghiệm 9.3.5 Phân tích hàm lượng số nguyên tố mẫu xi măng cà phê 9.3.6 Kết phân tích Mn Na mẫu oxit 64 65 65 65 66 68 70 71 77 77 78 78 80 81 81 88 90 90 90 93 100 102 103 104 108 109 109 117 120 121 121 122 123 125 127 129 iv 9.3.7 Nhận xét Chương 10 Xây dựng thực hành 10.1 Đo thông lượng neutron môi trường nước dùng nguồn đồng vị phương pháp kích hoạt dị 10.2 Đo độ độ dài làm chậm neutron môi trường nước dùng nguồn đồng vị 252Cf 10.3 Đo diện tích di cư (M2) neutron môi trường nước dùng nguồn đồng vị 252Cf 10.4 Xác định hàm lượng nguyên tố phương pháp phân tích kích hoạt dùng nguồn neutron đồng vị 10.5 Đo suất liều neutron môi trường nước dùng nguồn đồng vị phương pháp kích hoạt dị Chương 11 Thiết lập hai toán mẫu (Benchmark) 11.1 Bài tốn 1: Tính tốn mơ MCNP thí nghiệm đo phân bố thông lượng, đặc trưng làm chậm khuếch tán neutron hệ Howitzer chứa nước dùng nguồn đồng vị 241Am-Be 11.2 Bài tốn 2: Tính tốn mơ MCNP thí nghiệm đo phân bố thơng lượng, đặc trưng làm chậm khuếch tán neutron hệ Howitzer chứa nước dùng nguồn đồng vị 252Cf 131 132 132 Phần III KẾT LUẬN Kết thực Đề tài Sản phẩm KHCN ứng dụng Đề tài 2.1 Hệ thống thiết bị thực nghiệm đo neutron 2.2 Các cơng trình khoa học cơng bố 2.3 Các luận văn thạc sĩ hoàn thành 2.4 Hướng dẫn thực hành cho đoàn học viên Những đề nghị nghiên cứu 184 184 185 185 185 185 186 186 TÀI LIỆU THAM KHẢO 187 PHẦN PHỤ LỤC Phụ lục I (Chi tiêu kinh phí) Phụ lục II (Chương 7) Phụ lục III (Chương 8) Phụ lục IV (Chương 9) Phụ lục V (Các báo công bố) 190 140 144 150 159 169 169 176 v NHỮNG KÝ HIỆU VÀ CHỮ VIẾT TẮT - ATBX: An toàn xạ - DNL: Độ phi tuyến vi phân - ĐHBK: Đại học Bách khoa - ĐHKHTN: Đại học Khoa học tự nhiên - ĐHSP: Đại học Sư phạm - EVN (Vietnam Electricity): Tập đoàn Điện lực Việt Nam - FTC (Follow-up Training Course): Khóa huấn luyện (tự tổ chức) - JAEA (Japan Atomic Energy Agency): Cơ quan Năng lượng nguyên tử Nhật Bản - JTC (Joint Training Course): Cùng tổ chức khóa huấn luyện - HTQT: Hợp tác quốc tế - IAEA (International Atomic Energy Agency): Cơ quan Năng lượng nguyên tử quốc tế - INAA (Instrumental Neutron Activation Analysis): Phân tích kích hoạt neutron dụng cụ - INL: Độ phi tuyến tích phân - ITP (Instructor Training Program): Chương trình đào tạo giảng viên - IE (Applicaion of nuclear technique in industry and environment): Ứng dụng Kỹ thuật hạt nhân công nghiệp môi trường - KHCN: Khoa học công nghệ - KHKTHN: Khoa học kỹ thuật hạt nhân - KTHN: Kỹ thuật hạt nhân - L (Diffusion length): Chiều dài khuếch tán - LS (Slowing-down length): Chiều dài làm chậm - LabVIEW (Laboratory Virtual Instrument Engineering Workbench): Ngôn ngữ lập trình dạng sơ đồ hình ảnh - LOD (Limit of detection): Giới hạn xác định - M2 (Migration area): Diện tích di cư vi - LPƯHNĐL: Lị phản ứng hạt nhân Đà Lạt - NAA (Neutron Activation Analysis): Phân tích kích hoạt neutron - NCHN: Nghiên cứu hạt nhân - NLNT: Năng lượng nguyên tử - NLNTVN: Năng lượng nguyên tử Việt Nam - NRPB (Nuclear power and Renewable energy Projects pre-investment Boad): Ban chuẩn bị đầu tư dự án điện hạt nhân lượng tái tạo - NuTEC hay NuHRDeC (Nuclear Human Resourse Development Center): Trung t©m Phỏt trin ngun nhõn lc hạt nhân - PTN: Phũng thí nghiệm - to: Nhiệt độ - TTĐT: Trung tâm đào tạo - TTLPƯ: Trung tâm Lò phản ứng - TTMT: Trung tâm Môi trường - TTPT: Trung tâm Phân tích - VL-ĐTHN: Vật lý – Điện tử hạt nhân - VLHN: Vật lý hạt nhân vii TÓM TẮT Báo cáo trình bày kết nghiên cứu mơ Monte-Carlo để thiết kế cấu hình hệ thống thiết bị thực nghiệm Howitzer chứa nước đa chức dùng nguồn neutron đồng vị 252Cf 241Am-Be Từ chế tạo hệ Howitzer chứa nước, gồm: thùng nhơm lớp chứa nước hình trụ với kích thước 1,2mx1,2m kích thước ngồi 1,3mx1,3m; hệ trao đổi nhiệt bơm nước tuần hoàn; hệ giá đỡ di chuyển nguồn neutron; hệ đo neutron dùng ống đếm 3He BF3 Từ hệ thiết bị chế tạo này, nghiên cứu tiến hành về: đo thơng lượng suất liều neutron dùng dị kích hoạt vàng; đo chiều dài làm chậm khuếch tán diện tích di cư neutron mơi trường nước; khảo sát thông số phổ neutron xác định hàm lượng số nguyên tố mẫu xi măng cà phê; thiết lập toán mẫu; xây dựng thực hành liên quan phục vụ công tác đào tạo Từ kết đề tài, đăng tải báo khoa học, hoàn thành hướng dẫn luận văn thạc sĩ hệ Howitzer phục vụ cho đoàn sinh viên/cán sở khác đến Viện NCHN thực tập Sản phẩm thực tiễn đề tài hệ thiết bị thực nghiệm Howitzer chứa nước đa chức dùng nguồn neutron đồng vị nghiên cứu chế tạo lần Việt Nam, hữu ích cho công tác đào tạo nhân lực hạt nhân ABSTRACT This report presents and research results and Monte-Carlo simulation for designing a configuration of system on experimental equipment including multi-function water-filled Howitzer using isotope neutron sources of 252Cf and 241Am-Be From that, the water-filled Howitzer was made It consists of the blocks as follows: water-filled two-layer aluminum cylinder tank with the inside-sizes of 1.2mx1.2m and the outside-sizes of 1.3mx1.3m; temparature exchanging system and circulation water-ejector; system for supporting and moving neutron sources; neutron measuring system using 3He and BF3 counters From the Howitzer system, researches were carried out for: measuring neutron flux and rate-dose using the activated gold foils; measuring slowing-down and diffusion lengths as well as migration area of neutrons in water medium; surveying spectrum parameters of neutron and determining contents of some elements in cement and coffee samples; establishing two benchmarks; setting up five related experimentally exercises for training From the results of the Project, three scientific articles were issued, five Msc thesies were completed, and the Howitzer system served at once three teams of students/staffs from the other institutions to the Nuclear Research Institute for practice The practical product of the Project has been the multi-function water-filled Howitzer using isotope neutron sources that was researched and made fist time in Vietnam It has been very useful for traning activity in the nuclear field viii MỞ ĐẦU Tình hình nghiên cứu nước ngồi nước 1.1 Tình hình nghiên cứu nước ngồi Trên giới, để phục vụ nghiên cứu đào tạo nguồn nhân lực lĩnh vực vật lý neutron vật lý lò phản ứng, trường đại học (Đại học công nghệ Tokyo, Nhật Bản; Đại học MIT, Mỹ; Đại học Thanh Hoa, Bắc Kinh, Trung Quốc, v.v…), Trung tâm đào tạo thuộc Viện nghiên cứu (như NuHRDeC/JAEA, Nhật Bản; NTC/KAERI, Hàn Quốc, v.v.) có phịng thí nghiệm với đầy đủ trang thiết bị thực hành để để học viên tiến hành đo đạc đặc trưng vật lý neutron như: đo phổ thông lượng neutron, độ dài làm chậm khuếch tán neutron, kích hoạt mẫu nguồn neutron đồng vị, định liều neutron thực tính tốn mơ Ở NuHRDeC/JAEA, Nhật Bản có phịng thí nghiệm với đầy đủ trang thiết bị để nghiên cứu, đo đạc thực nghiệm đặc trưng làm chậm khuyếch tán neutron môi trường graphite môi trường nước sử dụng nguồn neutron đồng vị 252Cf 241Am-Be Vừa qua (năm 2006-2009), số cán Viện NCHN sang NuHRDeC/JAEA, Nhật Bản thực tập nghiên cứu phịng thí nghiệm theo Chương trình đào tạo giảng viên ITP (Instructor Training Program) Với phịng thí nghiệm, tùy theo thiết bị đo nguồn neutron sử dụng mà có tính tốn, thiết kế khơng gian làm việc phép đo đạc thực nghiệm khác nhằm phục vụ tốt công tác huấn luyện, đào tạo nghiên cứu cho đối tượng học viên khác (sinh viên, học viên cao học, thực tập sinh, v.v…) đảm bảo vấn đề an tồn xạ (ATBX) Các phịng thí nghiệm bắt buộc phải có sở đào tạo nhằm đào tạo nguồn nhân lực lĩnh vực hạt nhân nói chung điện hạt nhân nói riêng 1.2 Tình hình nghiên cứu nước Ở nước ta, có số sở nghiên cứu (như Viện KHKTHN Hà Nội, Viện Vật lý thuộc Viện KHCN Việt Nam) có số nghiên cứu lý thuyết thực nghiệm đặc trưng neutron rời rạc/nhỏ lẻ, chưa có hẳn phịng thí nghiệm độc lập để phục vụ công tác huấn luyện đào tạo, mà thực số thí nghiệm đo thơng lượng neutron, phân tích kích hoạt neutron cho khóa luận tốt nghiệp hay thực tập học viên Viện Nghiên cứu hạt nhân (NCHN), Đà Lạt sở nghiên cứu có Lị phản ứng nghiên cứu IVV-9 với công suất danh định 500 kW phục vụ tốt cho công tác huấn luyện, đào tạo, giảng dạy nghiên cứu lĩnh vực vật lý hạt nhân thực nghiệm nói chung, đặc biệt lĩnh vực vật lý neutron vật lý lị phản ứng nói riêng Tuy nhiên, để chuẩn hóa cơng tác huấn luyện đào tạo lĩnh vực sở đào tạo nước ngoài, trước tiến hành thực nghiệm lò phản ứng IVV-9, kênh chiếu xạ ngang lò phản ứng, kênh chiếu xạ khô (kênh No 13-2, No 7-1) học viên cần phải thực thí nghiệm (về đo thơng lượng phân bố thơng lượng neutron, phân tích mẫu kích hoạt dùng nguồn neutron đồng vị, làm chậm khuếch tán neutron môi trường vật lý khác graphite, nước nhẹ, v.v…) hệ thống thiết bị thực nghiệm đa phù hợp để có hiểu biết tốt đặc trưng neutron môi trường chất làm chậm, phương pháp phân tích kích hoạt neutron định liều neutron dùng nguồn đồng vị Nhưng Việt Nam chưa có sở nghiên cứu/đào tạo có Bảng Sự thay đổi chiều dài làm chậm Ls theo nhiệt độ nguồn 241Am-Be Nhiệt độ T(0C) 35 40 45 50 Ls ΔLs 7,883 7,883 7,882 7,883 0,8 0,9 0,8 0,8 Sai số tương đối (%) 10,1 11,4 10,1 10,1 2.3.2.2 Diện tích i cư chiều ài khuếch tán: Để xác định diện tích di cư phải đo tổng neutron nhiệt nhiệt Kết thực nghiệm diện tích di cư M2 chiều dài khuếch tán L phụ thuộc nhiệt độ môi trường nước nguồn 241Am-Be nêu Bảng Bảng Sự thay đổi chiều dài khuếch tán L theo nhiệt độ nguồn 241Am-Be Nhiệt độ T(0C) M2(cm2) L(cm) σ/σ0 Lσ/σ0(cm) 35 70,079 2,817 0,994 2,801 40 70,239 2,846 0,993 2,825 45 70,386 2,871 0,991 2,844 50 70,499 2,891 0,989 2,858 Với: σ0 khối lượng riêng nước 00C; σ khối lượng riêng nước nhiệt độ T 2.3.2.3 Chiều ài làm chậm khuếch tán neutron 35oC với nguồn 252Cf: Đo đạc tính tốn tương tự nguồn 241Am-Be, chiều dài làm chậm Ls nguồn 252 Cf 35oC xác định là: Ls = 6,0 ± 0,9 cm diện tích di cư M2 nêu Bảng Bảng Các giá trị đo M2 nguồn 252Cf nhiệt độ 350C Số lần đo M2(cm2) 44,332 44,331 44,334 Giá trị trung bình 44,3 Sai số 3,9 Sai số tương đối (%) 8,8 2.3.3 Khảo sát thông số phổ neutron khả xác định hàm lượng số nguyên tố mẫu xi măng cà phê 2.3.3.1 Khảo sát thông số phổ neutron: Dùng monitor Au dạng (Au-M2, Au-M3) monitor Mo dạng dây (Mo-1, Mo-2) Cân tính bề dày monitor, chiếu khơng bọc có bọc Cd, sau đo hoạt độ phổ kế gamma HP e Kết tính thông số фth, фe, α, f nêu Bảng Bảng Bảng Kết tính фth, фe, f bỏ qua ảnh hưởng hệ số che chắn neutron hệ số α фth (2,09 ± 0,21).105 (n/cm2.s) фe (1,13 ± 0,11).103 (n/cm2.s) f 184,1 ± 14,2 Bảng Kết tính фth, фe, α, f có tính đến ảnh hưởng hệ số che chắn neutron hệ số α фth фe α f (1,42 ± 0,18).105 (n/cm2.s) (1,44 ± 0,19).103 (n/cm2.s) 0,27 ± 0,04 98,4 ± 6,6 2.3.3.2 Phân tích hàm lượng số nguyên tố mẫu xi măng cà phê Trên sở xác định thông số phổ neutron vị trí chiếu, tiến hành thực nghiệm phân tích nguyên tố mẫu xi măng, cà phê, dùng mẫu chuẩn xi măng (NIST-114p) Kết phân tích phương pháp ko-NAA phương pháp tương đối mẫu xi măng Hà Tiên Sao Mai mẫu cà phê nêu tương ứng Bảng 10, Bảng 11 Bảng 12 Bảng 10 Kết phân tích hàm lượng Al, Na mẫu xi măng Hà Tiên PP k0-INAA PP tương đối Mẫu Al Na Al Na Xi măng Hà Tiên 8520 1920 20260 3994 Xi măng Hà Tiên 16490 1841 29080 2753 Xi măng Hà Tiên 12234 1621 20845 3425 Trung bình 12415 1794 23395 3391 Bảng 11 Kết phân tích hàm lượng Al, Na mẫu xi măng Sao Mai Mẫu Xi măng Sao Mai Xi măng Sao Mai Trung bình PP k0-INAA Al Na 5472 1471 3411 945 4442 1208 PP tương đối Al Na 13320 3222 9939 2278 11630 2750 Bảng 12 Kết phân tích hàm lượng Na mẫu cà phê phương pháp k0-NAA Mẫu Cà Phê Trung Nguyên Cà Phê Trung Nguyên Cà Phê Thanh Tiến Cà Phê Thanh Tiến Na 444 449