Xác định mật độ của gỗ bằng phương pháp gamma truyền qua

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM TP.HCM KHOA VẬT LÝ BỘ MÔN VẬT LÝ HẠT NHÂN -   - LÊ HỒNG THIỆN KHÓA LUẬN TỐT NGHIỆP XÁC ĐỊNH MẬT ĐỘ CỦA GỖ BẰNG PHƯƠNG PHÁP GAMMA TRUYỀN QUA THÀNH PHỐ HỒ CHÍ MINH – THÁNG NĂM 2021 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM THÀNH PHỐ HỒ CHÍ MINH KHOA VẬT LÝ BỘ MƠN VẬT LÝ HẠT NHÂN XÁC ĐỊNH MẬT ĐỘ CỦA GỖ BẰNG PHƯƠNG PHÁP GAMMA TRUYỀN QUA SINH VIÊN THỰC HIỆN: LÊ HỒNG THIỆN Cán hướng dẫn: PGS TS HOÀNG ĐỨC TÂM Chuyên ngành: Vật Lý học THÀNH PHỐ HỒ CHÍ MINH – THÁNG NĂM 2021 LỜI CẢM ƠN Để hoàn thành tốt đề tài khóa luận “Xác định mật độ gỗ phương pháp gamma truyền qua” Ngoài cố gắng nỗ lực thân, xin gửi lời cảm ơn đến người bên cạnh, giúp đỡ để tơi hồn thành khóa luận cách tốt Lời xin gửi lời cảm ơn sâu sắc đến Thầy Hoàng Đức Tâm – giảng viên Khoa Vật Lý, Trường Đại học Sư phạm Thành Phố Hồ Chí Minh Thầy hướng dẫn, hỗ trợ tơi q trình làm khóa luận truyền đạt kinh nghiệm quý báu Tôi xin cảm ơn đến anh chị trong phịng thí nghiệm Vật lý hạt nhân góp ý để tơi hồn thiện khóa luận Đồng thời, tơi xin cảm ơn ban chủ nhiệm khoa Vật lý tạo điều kiện thuận lợi sở vật chất để tơi hồn thành khóa luận Cuối cùng, tơi xin cảm ơn gia đình, bạn bè hỗ trợ tơi mặt tinh thần suốt q trình làm khóa luận TP.HCM, Ngày tháng năm 2021 Lê Hồng Thiện i XÁC NHẬN CỦA CHỦ TỊCH HỘI ĐỒNG XÁC NHẬN CỦA GIẢNG VIÊN HƯỚNG DẪN Thành phố Hồ Chí Minh, tháng năm 2021 Sinh viên thực Lê Hồng Thiện ii MỤC LỤC Trang LỜI CẢM ƠN i DANH MỤC BẢNG v DANH MỤC HÌNH ẢNH vi CÁC TỪ VIẾT TẮT vii MỞ ĐẦU CHƯƠNG CƠ SỞ LÝ THUYẾT 1.1 Tổng quan xạ gamma tương tác với vật chất 1.1.1 Hiệu ứng quang điện 1.1.2 Tán xạ Compton 1.1.3 Hiệu ứng tạo cặp .6 1.1.4 Suy giảm tia gamma 1.2 Phương pháp gamma truyền qua 1.2.1 Hệ số đóng góp B 1.2.2 Phương pháp gamma truyền qua .9 1.3 Phương pháp truyền sai số 12 CHƯƠNG PHƯƠNG PHÁP MONTE CARLO 14 VÀ CHƯƠNG TRÌNH MƠ PHỎNG MCNP6 14 2.1 Giới thiệu phương pháp Monte Carlo chương trình MCNP 14 2.1.1 Phương pháp Monte Carlo 14 2.1.2 Chương trình MCNP6 .15 2.1.3 Mơ hình mô MCNP6 16 2.2 Mô MCNP cho phương pháp gamma truyền qua 22 2.3 Phương pháp xử lý phổ 25 2.3.1 Giới thiệu chương trình Colegram 25 2.3.2 Các bước thực xử lý phổ 26 2.3.3 Hàm toán làm khớp 28 iii CHƯƠNG KẾT QUẢ VÀ THẢO LUẬN 29 3.1 Xây dựng đường chuẩn để xác định mật độ loại gỗ 29 3.1.1 Với mẫu gỗ dạng hình hộp chữ nhật .30 3.1.2 Với mẫu gỗ dạng hình trụ .34 3.2 Xác định mật độ số loại gỗ 37 KẾT LUẬN .40 KIẾN NGHỊ HƯỚNG NGHIÊN CỨU 41 TÀI LIỆU THAM KHẢO 42 iv DANH MỤC BẢNG Bảng 2.1 Phương trình số loại mặt định nghĩa MCNP 17 Bảng 2.2 Các thơng số đầu dị NaI(Tl) 21 Bảng 2.3 Bảng liệu thành phần mật độ loại gỗ 22 Bảng 2.4 Bảng liệu thành phần mật độ loại gỗ 23 Bảng 3.1 Hệ số suy giảm khối gỗ 29 Bảng 3.2 Xử lý diện tích đỉnh lượng tồn phần gỗ hình hộp chữ nhật 31 Bảng 3.3 Kết đánh giá mật độ nội suy gỗ hình hộp chữ nhật 33 Bảng 3.4 Kết truyền qua gỗ hình trụ 34 Bảng 3.5 Kết truyền qua gỗ hình trụ 35 Bảng 3.6 Đánh giá kết nội suy từ đường chuẩn hai loại gỗ 37 Bảng 3.7 Bảng nội suy mật độ gỗ chưa biết trước 38 Bảng 3.8 Bảng nội suy mật độ gỗ chưa biết trước 38 v DANH MỤC HÌNH ẢNH Hình 1.1 Hiệu ứng quang điện Hình 1.2 Tán xạ Compton Hình 1.3 Hiệu ứng tạo cặp Hình 1.4 Suy giảm tia gamma vật liệu cacbon chì Hình 1.5 Mơ hình phương pháp gamma truyền qua Hình 2.1 Cấu tạo ô mạng 18 Hình 2.2 Cấu trúc đầu dò NaI(Tl) 22 Hình 2.3 Mơ hình truyền qua cho mẫu gỗ hình trụ 24 Hình 2.4 Mơ hình truyền qua cho mẫu gỗ hình hộp chữ nhật 24 Hình 2.5 Hình ảnh phổ chương trình Colegram .26 Hình 2.6 Chọn vùng liệu cần ROI 27 Hình 2.7 Kết hàm làm khớp 28 Hình 3.1 Hệ số suy giảm khối loại gỗ 30 Hình 3.2 Số đếm theo mật độ hình hộp chữ nhật 32 Hình 3.3 Tỉ số R theo mật độ gỗ có dạng hình hộp chữ nhật 32 Hình 3.4 So sánh mật độ nội suy mật độ chuẩn 34 Hình 3.5 Số đếm theo mật độ gỗ hình trụ 35 Hình 3.6 Đồ thị hàm khớp đường chuẩn gỗ dạng hình trụ 36 Hình 3.7 Đánh giá kết nội suy từ đường chuẩn hai loại gỗ 36 Hình 3.8 So sánh giá trị mật độ gỗ tính tốn mật độ gỗ nội suy 39 vi CÁC TỪ VIẾT TẮT Chữ viết tắt NDT NIST Tiếng Anh Tiếng Việt Non-Destructive Testing Kiểm tra không phá hủy National Institute of Standards Viện Tiêu chuẩn Kỹ and Technology thuật quốc gia FWHM Full Width at Half Maximum MCNP Monte Carlo N-Particle vii Độ rộng nửa chiều cao đỉnh Monte Carlo N-hạt MỞ ĐẦU Kỹ thuật phân tích khơng hủy mẫu (NDT) kỹ thuật kiểm tra phân tích sử dụng dị tìm khuyết tật bên bề mặt vật cần kiểm tra mà không làm ảnh hưởng đến khả sử dụng đối tượng, dùng để phát khuyết tật vết nứt, rỗ khí, ngậm xỉ, mối hàn không thấu chân, kiểm tra độ ăn mịn vật liệu, độ cứng bê tơng cốt thép, đo bề dày vật liệu [1],…Việc dò khuyết tật cơng trình giúp tìm vị trí bị hư hỏng thơng số kỹ thuật chưa đạt chuẩn xảy biến dạng qua nhiều năm sử dụng Kiểm tra không hủy mẫu có hiệu mặt chi phí, hạn chế rủi ro, tăng cường an tồn xây dựng Vì kỹ thuật phát sai hỏng thay mặt hàng trước cố xảy mà không ảnh hưởng đến khả sử dụng Kiểm tra khơng hủy mẫu cách kiểm tra có độ xác cao thí nghiệm lặp lại nhiều lần với đối tượng công nghiệp Phương pháp kiểm tra gồm nhiều phương pháp khác chia thành hai nhóm kiểm tra khuyết tật sâu bên trong: phương pháp chụp ảnh phóng xạ, phương pháp siêu âm,… kiểm tra khuyết tật gần bề mặt phương pháp kiểm tra thẩm thấu lỏng, phương pháp kiểm tra bột từ, phương pháp kiểm tra dịng xốy,… [1] Ngồi ra, phương pháp gamma truyền qua gamma tán xạ đề cập để kiểm tra tính chất, cấu trúc vật liệu Năm 2012, Priyada cộng xác định mặt phân cách hai môi trường mật độ phương pháp gamma truyền qua gamma tán xạ [2] Nghiên cứu tập trung vào khả phát mặt phân cách lớp chất lỏng – chất lỏng, chất khí – chất lỏng cho kết xác với sai số lớn 7,73% Năm 2015, V.P.Singh cộng xác định hệ số suy giảm khối cho số loại polyme phương pháp mô Monte Carlo [3] Nghiên cứu sử dụng tia gamma với mức lượng khác để đánh giá thay đổi hệ số suy giảm khối Nghiên cứu sử dụng hình học mơ khác sử dụng cho thực nghiệm Kết thu so với liệu từ NIST có sai số 1% Năm 2020, Huỳnh Đình Chương cộng xây dựng phương pháp bán thực nghiệm để xác định mật độ chất lỏng đầu dò nhấp nháy NaI(Tl) [4] Trong phương pháp này, mật độ xác định cách xây dựng bình % gỗ Ogirisi, độ lệch nhỏ 0,54 % Trong 14 loại gỗ xét đến, có 12 loại gỗ có hệ số suy giảm khối lệch với giá trị trung bình nhỏ % Theo lý thuyết, mật độ độ phụ thuộc vào tỉ số R hệ số suy giảm khối tính số: m = Const Hình 3.1 Hệ số suy giảm khối loại gỗ Hàm (1.11) làm khớp dạng: ρi ( R i ) = a R i + b (3.1) hệ số a b tham số hàm khớp Trong đó, hệ số a dùng để làm khớp hệ số µ𝑚 d Hệ số b hệ số hiệu chỉnh để phương trình 3.1 phù hợp với phương trình 1.11 số giá trị gần xét đến Bởi vì, cho hệ số đóng góp B = bỏ qua tán xạ thứ cấp đến đầu dị tính gần cho hệ số suy giảm khối số 3.1.1 Với mẫu gỗ dạng hình hộp chữ nhật Kết xử lý thu diện tích đỉnh lượng tồn phần loại gỗ khác có mật độ thấp từ 0,405 g cm-3 đến 1,11 g cm-3 Kết trình bày bảng 3.2 Kết N diện tích đỉnh lượng tồn phần khơng 30 có gỗ Tỉ số R tỉ số diện tích đỉnh lượng tồn phần có gỗ so với khơng có gỗ Bảng 3.2 Xử lý diện tích đỉnh lượng tồn phần gỗ hình hộp chữ nhật Đỉnh hấp STT Tên gỗ Mật độ (g cm-3) 𝐍𝟎 thụ Tỉ số lượng R Sai số R toàn phần Spruce 0,405 61214660 51943760 0,164 0,05% Fir4 0,530 61214660 49408950 0,214 0,05% Pine1 0,610 61214660 47779230 0,248 0,05% Pine2 0,650 61214660 47006950 0,264 0,05% Oak1 0,680 61214660 46455720 0,276 0,05% Oak4 0,760 61214660 44966360 0,308 0,05% 1,040 61214660 39876800 0,429 0,05% 1,050 61214660 39664040 0,434 0,05% Tropical Almond Ogbono Shell Ogirisi 1,060 61214660 39502000 0,438 0,05% 10 Neem 1,110 61214660 38938560 0,452 0,05% Từ việc đánh giá kết quả, việc tăng mật độ lên khả chùm tia xuyên qua giảm nhiều Hai cách đưa để xây dựng đường chuẩn là: xây dựng đường tuyến tính mật độ theo tỉ số R xây dựng đường tuyến tính tỉ số R theo mật độ 31 Hai đường chuẩn xây dựng dùng để nội suy lại mật độ gỗ Để tốn đưa ứng dụng với mơ hình thực nghiệm sử dụng đường chuẩn tỉ số R phụ thuộc vào mật độ Hình 3.2 Số đếm theo mật độ hình hộp chữ nhật Hình 3.3 Tỉ số R theo mật độ gỗ có dạng hình hộp chữ nhật 32 Khi đánh giá tỉ số R theo mật độ đồ thị dạng hàm tuyến tính tốt Từ phương trình (1.11) ta có làm khớp theo dạng R i = a ρ + b với hệ số Hệ số Giá trị Sai số b -0,00645 0,00231 a 0,417 0,00285 Từ đường chuẩn vừa xây dựng được, thay lại giá trị tỉ số R để nội suy lại mật độ đánh độ lệch đường chuẩn so mật độ ban đầu Từ Bảng 3.3, độ lệch nhỏ 0,046 % với loại gỗ Pine1 độ lệch lớn 1,062 % với loại gỗ Spruce Bảng 3.3 Kết đánh giá mật độ nội suy gỗ hình hộp chữ nhật Mật độ Mật độ nội suy Sai số (g cm-3) (g cm-3) mật độ Spruce 0,405 0,409 0,24% 1,06% Fir4 0,530 0,529 0,24% 0,14% Pine1 0,610 0,610 0,24% 0,05% Pine2 0,650 0,649 0,24% 0,19% Oak1 0,680 0,677 0,26% 0,43% Oak4 0,760 0,755 0,25% 0,63% Tropical Almond 1,040 1,043 0,26% 0,32% Ogbonoshell 1,050 1,056 0,26% 0,58% Ogirisi 1,060 1,066 0,26% 0,56% 10 Neem 1,110 1,100 0,27% 0,87% STT Tên gỗ Độ lệch Sai số lớn mật độ tính tốn nhỏ 0,3% Tỉ số R2 0,9996 nên đường chuẩn đáp ứng tốt với liệu mơ 33 Trong hình 3.4 biểu diễn độ chênh lệch mật độ nội suy từ đường chuẩn Đường chuẩn mô xây dựng trùng khớp với liệu chuẩn Hình 3.4 So sánh mật độ nội suy mật độ chuẩn 3.1.2 Với mẫu gỗ dạng hình trụ Bảng 3.4 Kết truyền qua gỗ hình trụ STT Tên gỗ Mật độ (g cm-3) Đỉnh hấp thụ lượng toàn phần 𝐍𝟎 Tỉ số R Sai số R Spruce 0,405 44113020 61195260 0,327 0,05% Fir4 0,530 39879780 61195260 0,428 0,05% Pine1 0,610 37297150 61195260 0,495 0,05% Pine2 0,650 36112830 61195260 0,527 0,05% 34 Bảng 3.5 Kết truyền qua gỗ hình trụ STT Tên gỗ Mật độ (g cm-3) Đỉnh hấp thụ lượng 𝐍𝟎 toàn phần Tỉ số R Sai số R Oak1 0,680 35270130 61195260 0,551 0,06% Oak4 0,760 33067710 61195260 0,616 0,06% 1,040 26061140 61195260 0,854 0,06% Tropical Almond Ogbonoshell 1,050 25792580 61195260 0,864 0,06% Ogirisi 1,060 25578040 61195260 0,872 0,06% 10 Neem 1,110 24853520 61195260 0,901 0,06% Hình 3.5 Số đếm theo mật độ gỗ hình trụ 35 Hình 3.6 Đồ thị hàm khớp đường chuẩn gỗ dạng hình trụ Hệ số Giá trị Sai số b 0,0103 0,00405 a -0,828 0,0051 Hàm làm khớp với hệ số R2=0,9996, hàm làm khớp phù hợp với liệu Hình 3.7 Đánh giá kết nội suy từ đường chuẩn hai loại gỗ 36 Tỉ số R2 0,9997 nên đường chuẩn khớp với liệu mô Sai số nhỏ 0,44% gỗ Spruce, sai số lớn 0,62% gỗ Neem Sai số tăng lên mật độ gỗ lớn Đường chuẩn xây dựng đáp ứng tốt với vùng có mật độ thấp Bảng 3.6 Đánh giá kết nội suy từ đường chuẩn hai loại gỗ Mật độ Mật độ nội suy Sai số (g cm-3) (g cm-3) mật độ Spruce 0,405 0,408 0,44% 0,67% Fir4 0,530 0,530 0,46% 0,09% Pine1 0,610 0,610 0,48% 0,06% Pine2 0,650 0,649 0,49% 0,11% Oak1 0,680 0,678 0,49% 0,32% Oak4 0,760 0,756 0,51% 0,57% Tropical Almond 1,040 1,043 0,60% 0,30% Ogbonoshell 1,050 1,056 0,60% 0,54% Ogirisi 1,060 1,066 0,60% 0,54% 10 Neem 1,110 1,100 0,62% 0,86% STT Tên gỗ Độ lệch 3.2 Xác định mật độ số loại gỗ Chúng dùng hai đường chuẩn hai loại gỗ hình hộp chữ nhật gỗ hình trụ để kiểm tra lại mật độ gỗ kết thu Bảng 3.7, Bảng 3.8 Đối với mẫu gỗ dạng hình trụ cho mật độ nội suy lệch nhỏ 0,221 % lệch lớn 0,319 % với mật độ chuẩn Đối với mẫu gỗ dạng hình hộp chữ nhật cho độ lệch nhỏ 0,245 % độ lệch lớn 0,589 % so với mật độ chuẩn Tuy nhiên, sai số phép đo với mẫu gỗ hình trụ lớn so mẫu gỗ hình hộp chữ 37 nhật Kết kiểm tra cho thấy độ lệch hai đường chuẩn vừa xây dựng cho kết nhỏ % Bảng 3.7 Bảng nội suy mật độ gỗ chưa biết trước Tên gỗ Mật độ (g cm-3) Gỗ hình trụ Mật độ nội suy Sai số (g cm-3) mật độ Độ lệch % Fir3 0,435 0,436 0,45% 0,28% Pine 0,650 0,649 0,49% 0,22% Oak 0,720 0,718 0,51% 0,32% Mango 1,100 1,097 0,61% 0,31% Bảng 3.8 Bảng nội suy mật độ gỗ chưa biết trước Tên gỗ Mật độ (g cm-3) Gỗ hình hộp chữ nhật Mật độ nội suy Sai số (g cm-3) mật độ Độ lệch % Fir3 0,435 0,438 0,24% 0,59% Pine 0,650 0,648 0,24% 0,30% Oak 0,720 0,718 0,25% 0,25% Mango 1,100 1,097 0,27% 0,31% 38 Hình 3.8 So sánh giá trị mật độ gỗ tính tốn mật độ gỗ nội suy 39 KẾT LUẬN Đối với đề tài “ Xác định mật độ gỗ phương pháp gamma truyền qua”, tiến hành mơ tồn q trình thực nghiệm chương trình MCNP6 Từ kết mơ phỏng, tơi dùng hai đường chuẩn xây dựng từ hai hình dạng gỗ gỗ hình hộp chữ nhật gỗ hình trụ để tính tốn mật độ loại gỗ Kết thu từ nghiên cứu: Xây dựng hai đường chuẩn với hai hình dạng gỗ Kết đường chuẩn gỗ hình hộp chữ nhật có tỉ số R2 = 0,9996 Sai số lớn đường chuẩn tính tốn 0,27 % Kết đường chuẩn gỗ hình trụ có tỉ số R2 = 0,9946 Sai số lớn tính tốn 0,62 % Từ hai đường chuẩn trên, tơi tiến hành tính toán mật độ loại gỗ biết để so sánh kết tính với kết thực tế Độ lệch lớn gỗ hình hộp chữ nhật 0,31 %, sai số lớn mật độ 0,27 % Độ lệch lớn gỗ hình hộp trụ 0,30 %, sai số lớn mật độ 0,61 % Việc xây dựng đường chuẩn mơ MCNP6 với 10 loại gỗ có mật độ từ 0,405 (g cm-3) đến 1.110 (g cm-3) để đánh giá kết tuyến tính tỉ số R mật độ Hai đường chuẩn vừa xây dựng đáp ứng tốt với liệu 40 KIẾN NGHỊ HƯỚNG NGHIÊN CỨU Khóa luận xây dựng đường chuẩn mơ cho kết đáp ứng tốt với liệu đưa Trong q trình thực thí nghiệm cần phải ý đến độ ẩm môi trường đo Để mơ phổ cách xác cần phải biết xác thành phần loại gỗ Tuy nhiên, gỗ tự nhiên có mật độ thành phần khác tùy vào đặc điểm vùng Do đó, việc xác định xác thành phần gặp khó khăn Tơi có đề xuất số hướng nghiên cứu tiếp theo: • Tiến hành mơ xác định mật độ gỗ với thành phần nước gỗ thay đổi • Tiến hành thực nghiệm xác định mật độ gỗ 41 TÀI LIỆU THAM KHẢO [1] Vũ Tiến Hà (2018), “Ứng dụng phương pháp kiểm tra khơng phá hủy thơng dụng lĩnh vực khí Việt Nam”, hội nghị khoa học công nghệ tồn quốc khí lần thứ V [2] P Pryada, M Margret, R Ramar, Shivaramu (2012), “Intercomparison of gamma ray scattering and transmission techniques for fluid–fluid and fluid–air interface levels detection and density measurements,” Applied Radiation and Isotopes,volume 70, pp 462 – 469 [3] V P Singh, S P Shirmardi, M E Medhat, N M Badiger (2015), “Determination of mass attenuation coefficient for some polymers using Monte Carlo simulation”, Vacuum, volume 119, pp 284-288 [4] H D Chuong, L T N Trang, H D Tam,V H Nguyen, T T Thanh (2020), “A new approach for determining the thickness of material plate using gamma backscattering method”, NDT&E International, volume 113 • [5] Y Chen, Jinzhou Fu, D Baokang, S Qingfeng, L Huiqiao, and Z Tianyou (2020), “Artificial wooden nacre: a high specific strength engineering material”, Acs Nano, 14, pp 2036 – 2043 [6] K Nishimiya, T Hata, Y Imamura and S Ishihara (1997), “Analysis of chemical structure of wood charcoal by X-ray photoelectron spectroscopy,” Journal of Wood Science, pp 56 – 61 [7] B Goy, P Martin and J.-M Leban (1992), “The measurement of wood density by microwave sensor,” Springer-Verlag, 163 – 166 [8] F.S.Malan, P.G Marais (1991), “Direct gamma ray densitometry of wood,” South African Forestry Journal-no-157, pp – [9] Ngô Quang Huy, ”Cơ sở Vật lý hạt nhân”, NXB khoa học kỹ thuật, 2006 42 [10] J E Fernandez (1991), “Compton and Rayleigh double scattering of unpolarized radiation,” Physical Review A [11] Dr.GeorgSchlieper (2000), “Principles of gamma ray densitometry”, Metal Powder Report, Volume 55, pp 20-23 [12] A Kiyani, A A Karami, M Bahiraee and H Moghadamian (2012), “Calculation of gamma buildup factors for point sources,” Advances in Materials Research, pp 93-98 [13] P R.Bevington, D K Robinson (2003), “Data Reduction And Error Analysics For The Physical Sciences,” Emeritus Professor of Physics Case Western Reserve University [14] Đặng Nguyên Phương (2012), “Hướng dẫn sử dụng MCNP cho hệ điều hành Windows” [15] E C Okoroigwe, Z Li, S Kelkar, C Saffron, S Onyegegbu (2015), “Bio-oil yield potential of some tropical woody biomass”, Journal of Energy in Southern Africa, pp 33-41 [16] R J McConn , C J Gesh, R T Pagh, R A Rucker, R G Williams (2011), “Compendium of material composition data for radiation transport modeling”, Radiation Portal Monior Project [17] M C Lépy (2004), “Prensentation of the colegram software”, Laboratoire National Henri Becquerel Trang web [18] https://pdg.lbl.gov/2020/figures/passage/figures/sigma_both_06.pdfm Truy cập vào web Trang web tham khảo Source Particle Data Group, truy cập vào ngày 13/05/2021 [19]https://akela.mendelu.cz/~xcepl/inobio/nove/Wood_anatomy/WAEF-06micro_softwoods.pdf 43 Truy cập vào ngày 13/05/2021 [20] https://www.physics.nist.gov/PhysRefData/Xcom/html/xcom1.html Truy cập vào ngày 03/05/2021 44