Khóa Luận Tốt Nghiệp: Xác Định Mật Độ Gỗ Bằng Phương Pháp Gamma Truyền Qua

Khóa luận tốt nghiệp vật lý nghiên cứu phương pháp gamma truyền qua để xác định mật độ gỗ, mang lại ứng dụng thực tiễn trong ngành gỗ.

• Tác giả: Lê Hồng Thiện
• Người hướng dẫn: PGS. Hoàng Đức Tâm
• Trường học: Trường Đại Học Sư Phạm Thành Phố Hồ Chí Minh
• Chuyên ngành: Vật Lý Học
• Đề tài: Xác Định Mật Độ Gỗ Bằng Phương Pháp Gamma Truyền Qua
• Loại tài liệu: khóa luận tốt nghiệp
• Năm xuất bản: 2021
• Địa điểm: Thành Phố Hồ Chí Minh

Mục lục chi tiết

• LỜI CÁM ƠN
• 1. CHƯƠNG 1: TỔNG QUAN
• 1.1. Tông quan về bức xạ gamma tương tác với vật chất
• 1.2. Phương pháp gamma truyền qua
• 1.3. Phương pháp Monte Carlo và chương trình MCNP6
• 1.4. Kết quả và thảo luận
• 2. CHƯƠNG 2: MÔ HÌNH VÀ MÔ PHỎNG
• 2.1. Mô hình trong mô phỏng MCNP6
• 2.2. Giới thiệu chương trình Colegram
• 2.3. Kết quả và thảo luận
• 3. CHƯƠNG 3: XÁC ĐỊNH MẬT ĐỘ CỦA MỘT SỐ LOẠI GỖ
• KIẾN NGHỊ HƯỚNG NGHIÊN CỨU
• TÀI LIỆU THAM KHẢO
• DANH MỤC BẢNG
• DANH MỤC HÌNH ẢNH
• CÁC TỪ VIẾT TẮT

Tóm tắt

I. Tổng Quan Về Phương Pháp Gamma Truyền Qua Để Xác Định Mật Độ Gỗ

Phương pháp gamma truyền qua là một kỹ thuật hiện đại được sử dụng để xác định mật độ gỗ một cách chính xác. Kỹ thuật này dựa trên nguyên lý tương tác giữa bức xạ gamma và vật chất. Khi bức xạ gamma đi qua một mẫu gỗ, cường độ bức xạ sẽ giảm đi do sự hấp thụ và tán xạ. Việc đo lường sự suy giảm này cho phép xác định mật độ của gỗ một cách hiệu quả. Phương pháp này không chỉ đơn giản mà còn mang lại độ chính xác cao, giúp các nhà nghiên cứu và kỹ sư trong ngành gỗ có thể đánh giá chất lượng gỗ một cách nhanh chóng.

1.1. Nguyên Lý Cơ Bản Của Phương Pháp Gamma

Phương pháp gamma hoạt động dựa trên nguyên lý tương tác giữa bức xạ gamma và các electron trong vật chất. Khi bức xạ gamma đi vào mẫu gỗ, một phần năng lượng sẽ bị hấp thụ, dẫn đến sự suy giảm cường độ bức xạ. Điều này cho phép xác định mật độ gỗ thông qua các phép đo chính xác.

1.2. Lợi Ích Của Việc Sử Dụng Phương Pháp Gamma

Sử dụng phương pháp gamma truyền qua mang lại nhiều lợi ích, bao gồm khả năng đo lường không phá hủy, độ chính xác cao và khả năng áp dụng rộng rãi trong các lĩnh vực khác nhau. Kỹ thuật này giúp tiết kiệm thời gian và chi phí trong việc xác định mật độ gỗ.

II. Thách Thức Trong Việc Xác Định Mật Độ Gỗ Bằng Phương Pháp Gamma

Mặc dù phương pháp gamma truyền qua có nhiều ưu điểm, nhưng cũng tồn tại một số thách thức trong quá trình thực hiện. Một trong những vấn đề chính là ảnh hưởng của độ ẩm trong môi trường đến kết quả đo. Độ ẩm có thể làm thay đổi mật độ gỗ, dẫn đến sai số trong kết quả. Ngoài ra, việc chuẩn bị mẫu cũng cần phải được thực hiện cẩn thận để đảm bảo độ chính xác của phép đo.

2.1. Ảnh Hưởng Của Độ Ẩm Đến Kết Quả Đo

Độ ẩm trong gỗ có thể làm thay đổi mật độ, ảnh hưởng đến kết quả đo. Việc kiểm soát độ ẩm là rất quan trọng để đảm bảo tính chính xác của phương pháp gamma truyền qua.

2.2. Sai Sót Kỹ Thuật Trong Quá Trình Thực Hiện

Sai sót kỹ thuật trong quá trình chuẩn bị mẫu và thực hiện phép đo có thể dẫn đến kết quả không chính xác. Cần có quy trình chuẩn hóa để giảm thiểu sai sót này.

III. Phương Pháp Gamma Truyền Qua Quy Trình Thực Hiện Chi Tiết

Quy trình thực hiện phương pháp gamma truyền qua bao gồm nhiều bước quan trọng. Đầu tiên, cần chuẩn bị mẫu gỗ và thiết lập thiết bị đo. Sau đó, bức xạ gamma sẽ được phát ra và đi qua mẫu gỗ. Cuối cùng, cường độ bức xạ gamma sẽ được ghi nhận và phân tích để xác định mật độ gỗ.

3.1. Chuẩn Bị Mẫu Gỗ Để Đo

Mẫu gỗ cần được chuẩn bị cẩn thận, đảm bảo không có tạp chất và độ ẩm được kiểm soát. Việc này giúp đảm bảo kết quả đo chính xác và đáng tin cậy.

3.2. Thiết Lập Thiết Bị Đo

Thiết bị đo cần được hiệu chuẩn trước khi thực hiện phép đo. Việc này giúp đảm bảo rằng các kết quả thu được là chính xác và có thể so sánh được.

IV. Ứng Dụng Thực Tiễn Của Phương Pháp Gamma Trong Ngành Gỗ

Phương pháp gamma truyền qua đã được áp dụng rộng rãi trong ngành gỗ để xác định mật độ gỗ. Kỹ thuật này không chỉ giúp đánh giá chất lượng gỗ mà còn hỗ trợ trong việc phân loại và lựa chọn gỗ cho các ứng dụng khác nhau. Nhiều nghiên cứu đã chỉ ra rằng phương pháp này mang lại kết quả chính xác và đáng tin cậy.

4.1. Đánh Giá Chất Lượng Gỗ

Phương pháp gamma giúp đánh giá chất lượng gỗ một cách nhanh chóng và hiệu quả. Điều này rất quan trọng trong ngành sản xuất và chế biến gỗ.

4.2. Phân Loại Gỗ Theo Mật Độ

Việc phân loại gỗ theo mật độ giúp các nhà sản xuất có thể lựa chọn loại gỗ phù hợp cho từng ứng dụng cụ thể, từ đó nâng cao hiệu quả sản xuất.

V. Kết Luận Về Phương Pháp Gamma Truyền Qua Trong Xác Định Mật Độ Gỗ

Phương pháp gamma truyền qua là một công cụ hữu ích trong việc xác định mật độ gỗ. Mặc dù còn một số thách thức cần khắc phục, nhưng với những lợi ích mà nó mang lại, phương pháp này chắc chắn sẽ tiếp tục được phát triển và ứng dụng rộng rãi trong tương lai. Việc cải tiến quy trình và công nghệ sẽ giúp nâng cao độ chính xác và hiệu quả của phương pháp này.

5.1. Tương Lai Của Phương Pháp Gamma Trong Nghiên Cứu

Phương pháp gamma có tiềm năng lớn trong nghiên cứu và ứng dụng trong ngành gỗ. Việc phát triển công nghệ mới sẽ giúp cải thiện độ chính xác và khả năng ứng dụng của phương pháp này.

5.2. Khuyến Nghị Đối Với Nghiên Cứu Tiếp Theo

Cần tiếp tục nghiên cứu và phát triển các phương pháp mới để cải thiện độ chính xác của việc xác định mật độ gỗ. Việc này sẽ giúp nâng cao chất lượng sản phẩm trong ngành gỗ.

Nội dung chính

Bối cảnh và vấn đề nghiên cứu

Trong kỹ thuật công nghiệp và đời sống, gỗ là loại vật liệu tự nhiên được ứng dụng rộng rãi trong sản xuất đồ nội thất, kết cấu xây dựng và chế biến nhiên liệu. Giá trị sử dụng, tính chất cơ học, độ bền uốn, độ ổn định kích thước, khả năng gia công và bảo quản của gỗ phụ thuộc mật thiết vào mật độ khối. Mật độ còn là tiêu chí quan trọng để phân loại phẩm cấp gỗ, đánh giá cấu trúc tế bào và phân biệt gỗ già với gỗ non. Tuy nhiên, việc xác định mật độ gỗ trong thực tế gặp nhiều trở ngại do độ ẩm môi trường biến thiên, thành phần hóa học phức tạp (thay đổi theo loài, tuổi cây và điều kiện thổ nhưỡng từng vùng địa lý), cũng như các sai số cơ học phát sinh trong quá trình gia công, cắt gọt mẫu đo.

Kiểm tra không phá hủy (Non-Destructive Testing - NDT) là phương pháp phân tích cho phép kiểm tra cấu trúc bên trong hoặc bề mặt vật liệu mà không làm tổn hại đến khả năng sử dụng của đối tượng. Trong nhóm kỹ thuật NDT chiều sâu, phương pháp gamma truyền qua và gamma tán xạ sở hữu nhiều ưu thế về tính chuẩn xác và khả năng lặp lại. Các công trình thực nghiệm trước đây thường đòi hỏi chế tạo hàng loạt mẫu vật chuẩn và tiến hành đo đạc nhiều lần, dẫn đến chi phí cao và dễ phát sinh sai số do nhiệt độ, độ ẩm môi trường. Sự phát triển của các phương pháp tính toán số học, điển hình là kỹ thuật Monte Carlo, mở ra giải pháp mô phỏng các quá trình vật lý hạt nhân trong môi trường lý tưởng, giúp loại bỏ các tác nhân gây nhiễu môi trường trước khi ứng dụng vào đo đạc thực tế.

Khóa luận tốt nghiệp của tác giả Lê Hồng Thiện tập trung giải quyết bài toán trên thông qua việc kết hợp mô phỏng số và xử lý phổ hạt nhân. Các mục tiêu và nhiệm vụ nghiên cứu cụ thể trong văn bản bao gồm:

1. Tổng quan cơ sở lý thuyết về các cơ chế tương tác bức xạ gamma với vật chất (hiệu ứng quang điện, tán xạ Compton, hiệu ứng tạo cặp) và quy luật suy giảm của chùm photon truyền qua môi trường hấp thụ.
2. Mô hình hóa cấu hình đo bức xạ gamma truyền qua bằng chương trình mô phỏng Monte Carlo MCNP6, bao gồm thiết lập hệ chuẩn trực chì, nguồn phóng xạ Cs-137, đầu dò nhấp nháy NaI(Tl) và các mẫu gỗ với hai cấu hình hình học (hình hộp chữ nhật và hình trụ tròn).
3. Ứng dụng phần mềm Colegram để làm khớp phổ năng lượng photon ghi nhận được, bóc tách diện tích đỉnh hấp thụ năng lượng toàn phần khỏi nền tán xạ Compton.
4. Xây dựng đường chuẩn tương quan giữa tỉ số truyền qua $R$ và mật độ khối đối với 10 loại gỗ trong khoảng mật độ từ 0,405 g/cm³ đến 1,110 g/cm³, sau đó kiểm chứng độ chính xác nội suy trên 4 loại gỗ thử nghiệm.

Đối tượng nghiên cứu của khóa luận là đặc trưng suy giảm chùm tia gamma đơn năng khi đi qua 14 loại vật liệu gỗ tự nhiên có thành phần hóa học và mật độ khối khác nhau. Phạm vi nghiên cứu về mặt không gian tập trung vào việc mô hình hóa hệ đo tại Phòng thí nghiệm Vật lý hạt nhân, Trường Đại học Sư phạm TP.HCM; phạm vi điều kiện vật lý xét ở trạng thái mẫu gỗ sấy khô lý tưởng (độ ẩm $\theta = 0$) với độ dày vật liệu cố định là 5 cm, dưới năng lượng phát xạ 661,66 keV của đồng vị phóng xạ Cs-137.

Cơ sở lý thuyết và phương pháp

Khung lý thuyết của nghiên cứu dựa trên quy luật suy giảm cường độ bức xạ điện từ theo định luật Beer-Lambert. Khi chùm photon hẹp đi qua môi trường đồng nhất có bề dày $d$ (cm) và khối lượng riêng $\rho$ (g/cm³), mối liên hệ giữa số đếm xung ghi nhận tại đỉnh toàn phần khi có mẫu ($N$) và khi không có mẫu ($N_0$) tuân theo công thức:

$$N = N_0 e^{-\mu_l d} = N_0 e^{-\mu_m \rho d}$$

Trong đó, $\mu_l$ là hệ số suy giảm tuyến tính (cm⁻¹) và $\mu_m = \mu_l / \rho$ là hệ số suy giảm khối (cm²/g). Đối với vật liệu xốp có chứa nước, hệ số suy giảm tuyến tính phụ thuộc vào độ ẩm $\theta$ theo biểu thức $\mu_l = \mu_{mw} \rho_w + \theta \mu_{mw} \rho_w$. Trong mô phỏng này, mẫu gỗ được đặt ở điều kiện sấy khô tuyệt đối nên $\theta = 0$. Thành phần hóa học cấu tạo của gỗ chủ yếu gồm ba nguyên tố hydro, cacbon và oxy (tổng lượng cacbon và oxy chiếm phần lớn khối lượng). Do đó, ở vùng năng lượng photon cao của nguồn Cs-137, hệ số suy giảm khối $\mu_m$ của các loại gỗ biến thiên không đáng kể và được coi như một đại lượng gần như không đổi.

Tác giả thiết lập tỉ số truyền qua $R = N / N_0$. Mối tương quan giữa mật độ $\rho$ và tỉ số $R$ được mô tả qua hàm tuyến tính $\rho(R) = aR + b$ hoặc $R(\rho) = a\rho + b$, trong đó $a$ và $b$ là các tham số làm khớp thực nghiệm. Tham số $b$ đóng vai trò hiệu chỉnh các sai lệch do giả định hệ số suy giảm khối là hằng số và việc giả định hệ số tích lũy $B = 1$. Việc sử dụng tỉ số $R$ giúp loại trừ ảnh hưởng của độ trôi thiết bị, biến thiên nhiệt độ và các yếu tố môi trường trong quá trình đo đạc.

Hệ thống phương pháp và công cụ nghiên cứu bao gồm:

• Phương pháp mô phỏng Monte Carlo trên chương trình MCNP6: MCNP6 là chương trình tính toán chuyên dụng trong vật lý hạt nhân, xử lý bài toán vận chuyển bức xạ thông qua việc gieo số ngẫu nhiên theo các hàm phân bố xác suất tương tác hạt, sử dụng kỹ thuật DXTRAN và thuật toán Russian Roulette để theo dõi dấu vết tương tác của photon trong không gian Descartes ba chiều. Cấu hình mô phỏng gồm:
  • Khai báo loại hạt: Thẻ MODE 2 (vận chuyển photon).
  • Nguồn bức xạ: Thẻ SDEF phát hạt photon đẳng hướng năng lượng $ERG = 0,661657\text{ MeV}$ (tương ứng phát xạ gamma của Cs-137), đặt tại gốc tọa độ bên trong buồng che chắn chì có kích thước dài 18,1 cm, rộng 10,2 cm, cao 10,2 cm, bán kính lỗ chuẩn trực nguồn 1,5 cm. Đường kính trong của khối chuẩn trực nguồn chì là 0,251 cm và chuẩn trực đầu dò là 0,502 cm. Khoảng cách từ nguồn đến mẫu gỗ là 25,9 cm, từ mẫu gỗ đến bề mặt đầu dò là 17,3 cm.
  • Số lượng lịch sử hạt mô phỏng: $6 \times 10^7$ hạt, phổ năng lượng phân giải thành 8192 kênh.
  • Đầu dò nhấp nháy NaI(Tl): Mô phỏng chi tiết cấu trúc đa lớp gồm vỏ nhôm bảo vệ dày 0,1 cm, dài 7,79 cm; lớp không khí mật độ 0,001205 g/cm³; lớp phản xạ nhôm oxit; tinh thể nhấp nháy NaI có bề dày 3,81 cm, chiều dài 7,62 cm và mật độ 3,67 g/cm³.
  • Mẫu vật liệu: Khai báo bằng thẻ vật liệu Mn ZAID fraction kết hợp lệnh chuyển hệ trục tọa độ TRn. Mẫu gỗ được khảo sát ở hai dạng hình học có cùng chiều dày truyền qua 5 cm: hình hộp chữ nhật ($40\text{ cm} \times 20\text{ cm} \times 5\text{ cm}$) và hình trụ tròn (chiều cao 12,4 cm, bán kính 5 cm).
• Phương pháp xử lý phổ năng lượng bằng chương trình Colegram: Ứng dụng thuật toán phi tuyến Levenberg-Marquardt kết hợp giữa phương pháp gradient và ma trận Hessian để cực tiểu hóa độ lệch bình phương trong vùng đỉnh năng lượng quan tâm (ROI). Hàm giải tích làm khớp gồm hàm nền One-step nhằm loại bỏ phông tán xạ Compton liên tục và hàm Gauss để tính toán chính xác diện tích đỉnh hấp thụ năng lượng toàn phần.
• Phương pháp truyền sai số thống kê: Bức xạ phân rã tuân theo phân bố Poisson. Sai số tỉ số truyền qua $\sigma_R$ được xác định với hệ số phủ $k = 2,58$, tương ứng độ tin cậy thống kê 99%. Sai số mật độ nội suy $\sigma_\rho$ được tính toán kết hợp từ sai số thống kê số đếm và sai số của các tham số hàm khớp $a$, $b$.
• Nguồn dữ liệu: Thành phần nguyên tố hóa học và mật độ danh định của 14 loại gỗ được kế thừa từ các công bố thực nghiệm của NIST và các tài liệu chuyên ngành quốc tế.

Nội dung chính theo từng chương

Chương 1: Cơ sở lý thuyết

Chương 1 hệ thống hóa toàn diện các cơ chế tương tác vi mô của photon với môi trường vật chất. Ba hiệu ứng chính được phân tích gồm:

• Hiệu ứng quang điện: Chiếm ưu thế ở dải năng lượng thấp, photon tương tác truyền toàn bộ năng lượng cho electron liên kết của nguyên tử bia, làm bứt electron ra khỏi nguyên tử.
• Tán xạ Compton: Tương tác tán xạ đàn hồi giữa photon và electron tự do, photon bị lệch hướng và mất một phần năng lượng cho electron giật lùi; là hiệu ứng chi phối chủ yếu ở dải năng lượng trung bình (quanh mức 662 keV của Cs-137).
• Hiệu ứng tạo cặp: Xảy ra khi năng lượng photon tới vượt quá hai lần năng lượng nghỉ của electron (lớn hơn 1022 keV) trong điện trường hạt nhân, tạo ra một cặp electron - positron (được Patrick Blackett phát hiện năm 1948).

Nội dung chương tiếp tục phân tích sự suy giảm của tia gamma trong vật liệu, chỉ rõ nguyên tử có số khối $Z$ lớn như chì có khả năng cản xạ và tiết diện tán xạ vượt trội so với cacbon. Khi chùm tia đi qua vật liệu dày, tán xạ thứ cấp xuất hiện làm chùm tia bị phân kỳ, đòi hỏi phải đưa thêm hệ số tích lũy $B$ vào định luật hấp thụ. Để triệt tiêu bức xạ thứ cấp đưa hệ số $B$ về xấp xỉ 1, tác giả thiết lập điều kiện hình học chùm tia hẹp thông qua hai ống chuẩn trực chì cho nguồn và đầu dò. Cuối chương, tác giả xây dựng lý thuyết lan truyền sai số thực nghiệm, thiết lập công thức tính toán độ bất định cho tỉ số $R$ và mật độ $\rho$ ở mức độ tin cậy 99%.

Chương 2: Phương pháp Monte Carlo và chương trình mô phỏng MCNP6

Chương 2 mô tả chi tiết phương pháp tính toán ngẫu nhiên Monte Carlo bắt nguồn từ Dự án Manhattan do Stanislaw Ulam phát triển, chuyên giải quyết các bài toán tích phân nhiều chiều và vận chuyển hạt phức tạp. Tác giả trình bày cấu trúc mã nguồn MCNP6 (gồm khoảng 44.000 dòng lệnh FORTRAN, 1.000 dòng lệnh C và 400 chương trình con), nguyên lý mã hóa ba loại hạt (neutron ký hiệu 1, photon ký hiệu 2, electron ký hiệu 3) cùng phương thức khai báo không gian hình học thông qua các toán tử giao, hội, bù giữa các mặt phẳng, mặt cầu và mặt trụ trong hệ tọa độ Descartes.

Quy trình mô phỏng hệ đo gamma truyền qua được thiết lập đồng bộ với số hạt phát xạ $6 \times 10^7$ photon. Mười bốn loại gỗ khảo sát được chia thành hai nhóm: 10 loại dùng để xây dựng đường chuẩn (gồm Neem, Fir4, Pine1, Pine2, Oak1, Spruce, Oak4, Ogirisi, Ogbonoshell, Tropical Almond) và 4 loại dùng làm mẫu kiểm tra đánh giá độc lập (gồm Fir3, Mango, Oak, Pine). Cấu trúc đầu dò NaI(Tl) và buồng chì được mô hình hóa sát với thông số kỹ thuật thực tế. Phần cuối chương trình bày chi tiết quy trình 4 bước xử lý phổ trên phần mềm Colegram: (1) Mở tệp số liệu và hiệu chỉnh thang đo năng lượng; (2) Thiết lập vùng ROI bao quanh đỉnh quang điện 661,66 keV; (3) Khai báo hàm nền One-step để làm khớp phông Compton; (4) Khai báo hàm Gauss để làm khớp đỉnh toàn phần, từ đó bóc tách diện tích đỉnh thực với độ chính xác cao.

Chương 3: Kết quả và thảo luận

Chương 3 trình bày các kết quả tính toán số học thu được từ MCNP6 và Colegram. Phân tích dữ liệu hệ số suy giảm khối $\mu_m$ từ thư viện NIST cho thấy, dù mật độ 14 loại gỗ trải rộng từ 0,405 g/cm³ đến 1,110 g/cm³, hệ số $\mu_m$ hầu như không đổi quanh giá trị trung bình (dao động từ 0,077 cm²/g đến 0,082 cm²/g). Độ lệch lớn nhất của $\mu_m$ so với trung bình là 2,0% (ở gỗ Ogirisi) và nhỏ nhất là 0,54%; có 12/14 loại gỗ có độ lệch dưới 1%. Điều này khẳng định cơ sở lý thuyết xem $\mu_m$ là hằng số là hoàn toàn chính xác.

Tiếp đó, tác giả tiến hành xây dựng và đánh giá hai bộ đường chuẩn:

• Mẫu gỗ hình hộp chữ nhật ($40\text{ cm} \times 20\text{ cm} \times 5\text{ cm}$): Số đếm diện tích đỉnh toàn phần giảm đều đặn từ $51.943.760$ xung (đối với gỗ Spruce mật độ 0,405 g/cm³) xuống $24.856.320$ xung (đối với gỗ Ogbonoshell mật độ 1,083 g/cm³) so với số đếm không mẫu $N_0 = 61.214.660$ xung. Đồ thị tương quan giữa tỉ số $R$ và mật độ $\rho$ đạt hệ số xác định tuyến tính $R^2 = 0,9996$ (hệ số Pearson $r = 0,99981$). Đánh giá nội suy trên 10 loại gỗ làm chuẩn cho thấy độ lệch giữa mật độ tính toán và mật độ danh định dao động từ 0,046% (gỗ Pine1) đến 1,062% (gỗ Spruce).
• Mẫu gỗ hình trụ tròn (chiều cao 12,4 cm, bán kính 5 cm): Số đếm diện tích đỉnh toàn phần của mẫu trụ cũng giảm tuyến tính tương ứng theo mật độ, với $N_0 = 61.195.260$ xung. Hàm khớp đường chuẩn đạt hệ số $R^2 = 0,9996$ (hoặc 0,9997). Sai số tính toán mật độ nội suy của các mẫu chuẩn dao động từ 0,44% (Spruce) đến 0,62% (Neem).
• Kiểm chứng độc lập trên các mẫu gỗ thử nghiệm: Tác giả áp dụng hai đường chuẩn vừa xây dựng để xác định mật độ của các mẫu gỗ chưa biết trước gồm Fir3, Pine và Oak. Kết quả kiểm chứng chứng minh cả hai dạng hình học đều cho kết quả nội suy có độ lệch so với giá trị thực tế dưới 1%.

Kết quả và đóng góp

Nghiên cứu đã đạt được các kết quả định lượng cụ thể thông qua số liệu trích xuất từ mô phỏng MCNP6 và xử lý đỉnh phổ Colegram:

• Xác nhận tính bất biến tương đối của hệ số suy giảm khối $\mu_m$ đối với chùm tia photon năng lượng 661,66 keV đi qua các loại gỗ tự nhiên: độ lệch của 12/14 loại gỗ so với giá trị trung bình là dưới 1%.
• Thiết lập thành công hai phương trình đường chuẩn tuyến tính mô tả mối quan hệ giữa tỉ số truyền qua $R$ và mật độ gỗ cho hai dạng hình học với độ tương quan cao ($R^2 = 0,9996$).
• Kiểm chứng khả năng xác định mật độ trên các mẫu gỗ thử nghiệm độc lập. Bảng số liệu dưới đây tổng hợp kết quả đối sánh giữa mật độ thực tế và mật độ nội suy:

Dạng hình học mẫu	Tên loại gỗ	Mật độ chuẩn (g/cm³)	Mật độ nội suy (g/cm³)	Sai số mật độ (%)	Độ lệch tương đối (%)
Hình trụ tròn	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%
Hình hộp chữ nhật	Fir3	0,435	0,438	-	0,589% (hoặc 0,31%)
	Pine	0,650	-	-	-
	Oak	0,720	-	-	0,245%

Các kết quả trên cho thấy:

• Đối với mẫu gỗ dạng hình trụ, độ lệch giữa mật độ nội suy và mật độ chuẩn dao động trong khoảng từ 0,221% đến 0,319%, sai số phép đo nằm trong khoảng 0,45% - 0,51%.
• Đối với mẫu gỗ dạng hình hộp chữ nhật, độ lệch giữa mật độ nội suy và mật độ chuẩn dao động trong khoảng từ 0,245% đến 0,589% (tổng kết chung toàn bộ các mẫu kiểm tra ghi nhận độ lệch cực đại dưới 0,31% - 0,59%).
• Mẫu hình hộp chữ nhật cho sai số thống kê nhỏ hơn so với mẫu hình trụ tròn, song cả hai cấu hình đều đạt độ chính xác cao với sai số tổng thể nhỏ hơn 1%.

Về mặt đóng góp khoa học và giải pháp kỹ thuật, nghiên cứu đã đề xuất một quy trình bán thực nghiệm tiên tiến: sử dụng chương trình MCNP6 để thiết lập đường chuẩn phân tích thay cho việc phải chế tạo nhiều mẫu chuẩn cơ học phức tạp. Cách tiếp cận này loại trừ triệt để sự ảnh hưởng của độ ẩm môi trường và sai sót gia công mẫu, giúp tối ưu hóa thời gian và chi phí thực nghiệm nhưng vẫn bảo đảm độ tin cậy học thuật cao.

Hạn chế và hướng nghiên cứu tiếp

Khóa luận thẳng thắn chỉ ra các hạn chế xuất phát từ phạm vi nghiên cứu lý thuyết và điều kiện mô phỏng:

• Mô hình MCNP6 được thiết lập trong điều kiện lý tưởng với mẫu gỗ sấy khô hoàn toàn ($\theta = 0$), chưa tính đến sự tác động trực tiếp của độ ẩm môi trường tự nhiên lên hệ số suy giảm khối và mật độ thực tế của gỗ trong môi trường làm việc thông thường.
• Thành phần hóa học của gỗ trong tự nhiên biến động phụ thuộc vào thổ nhưỡng và khí hậu từng vùng địa lý, gây khó khăn cho việc khai báo chính xác tuyệt đối tỷ lệ các nguyên tố trong thẻ vật liệu của MCNP6.
• Toàn bộ kết quả mới dừng lại ở mức độ mô phỏng số học trên máy tính kết hợp xử lý phổ giải tích, chưa được đối chứng trực tiếp với hệ thực nghiệm đo đạc thực tế tại phòng thí nghiệm.

Từ các hạn chế trên, tác giả đề xuất hai hướng nghiên cứu tiếp theo:

1. Tiếp tục tiến hành mô phỏng bài toán truyền qua gamma trên MCNP6 với hàm lượng nước trong gỗ thay đổi ($\theta > 0$) để khảo sát định lượng ảnh hưởng của độ ẩm môi trường đến đường chuẩn mật độ.
2. Thiết lập hệ thực nghiệm đo đạc thực tế tại phòng thí nghiệm Vật lý hạt nhân với nguồn Cs-137 và đầu dò NaI(Tl) để kiểm chứng, hiệu chỉnh các số liệu thu được từ mô phỏng.

Giá trị tham khảo

Khóa luận là tài liệu tham khảo có giá trị học thuật và thực tiễn cho nhiều đối tượng:

• Sinh viên, học viên chuyên ngành Vật lý hạt nhân, Vật lý kỹ thuật: Cung cấp tài liệu tham khảo chi tiết về quy trình lập trình mô phỏng MCNP6 (khai báo nguồn tổng quát SDEF, thẻ vật liệu hợp chất, thẻ chuyển hệ tọa độ TRn, thiết lập hình học đầu dò nhấp nháy NaI(Tl) và hệ chuẩn trực chì).
• Người làm nghiên cứu trong lĩnh vực NDT và ứng dụng bức xạ: Cung cấp phương pháp luận chuẩn xác về kỹ thuật xử lý phổ bức xạ gamma trên phần mềm Colegram, phương pháp bóc tách đỉnh quang điện bằng hàm Gauss và hàm nền One-step theo thuật toán Levenberg-Marquardt.
• Kỹ sư và nhà nghiên cứu công nghệ vật liệu, lâm sản: Cung cấp cơ sở dữ liệu về thành phần nguyên tố, mật độ khối và hệ số suy giảm photon của 14 loại gỗ tự nhiên, cùng giải pháp xác định nhanh khối lượng riêng của gỗ mà không phá hủy mẫu vật.

Câu hỏi thường gặp

Tại sao hệ số tích lũy bức xạ $B$ trong nghiên cứu này được xem xấp xỉ bằng 1?

Trong điều kiện hình học chùm tia rộng, các photon tán xạ Compton và tán xạ Rayleigh thứ cấp từ vật liệu bia, không khí và cấu kiện che chắn có thể đi vào đầu dò, làm sai lệch định luật suy giảm hàm mũ thông thường. Để khắc phục, tác giả đã bố trí hệ hai khối chuẩn trực bằng chì có đường kính trong rất nhỏ (0,251 cm ở nguồn và 0,502 cm ở đầu dò) với bề dày đủ lớn. Cấu hình này tạo ra chùm tia song song, hẹp và ngăn chặn hầu hết các bức xạ tán xạ nhiều lần đến đầu dò, khiến phần đóng góp của bức xạ thứ cấp là không đáng kể, do đó hệ số tích lũy $B$ được coi xấp xỉ bằng 1.

Nguồn phóng xạ và đầu dò NaI(Tl) được cấu hình trong MCNP6 với những thông số nào?

Nguồn phóng xạ được khai báo qua thẻ SDEF phát hạt photon năng lượng 0,661657 MeV (tương ứng đồng vị Cs-137), đặt tại gốc tọa độ trong buồng chì kích thước $18,1\text{ cm} \times 10,2\text{ cm} \times 10,2\text{ cm}$ với bán kính chuẩn trực 1,5 cm. Đầu dò NaI(Tl) được mô phỏng gồm lớp vỏ nhôm dày 0,1 cm, dài 7,79 cm; lớp không khí mật độ 0,001205 g/cm³; lớp phản xạ nhôm oxit; và tinh thể nhấp nháy NaI có bề dày 3,81 cm, chiều dài 7,62 cm, mật độ 3,67 g/cm³. Khoảng cách từ nguồn đến mẫu gỗ là 25,9 cm và từ mẫu gỗ đến mặt đầu dò là 17,3 cm; số lượng lịch sử hạt mô phỏng là $6 \times 10^7$ hạt trên 8192 kênh năng lượng.

Hai hàm toán học nào được chương trình Colegram sử dụng để xác định diện tích đỉnh hấp thụ năng lượng toàn phần?

Chương trình Colegram sử dụng thuật toán Levenberg-Marquardt để khớp hai hàm giải tích: hàm nền One-step (STEP) dùng để mô phỏng và loại bỏ phông tán xạ Compton liên tục phía dưới đỉnh phổ, và hàm phân bố chuẩn Gauss dùng để làm khớp hình dạng đỉnh hấp thụ năng lượng toàn phần. Diện tích thực của đỉnh quang điện được tính bằng diện tích dưới đường Gauss trừ đi phần diện tích phông nền One-step tương ứng.

Vì sao hệ số suy giảm khối của các loại gỗ khác nhau có thể coi như một hằng số ở mức năng lượng 662 keV?

Thành phần hóa học chủ yếu của các loại gỗ tự nhiên là hợp chất hữu cơ cấu tạo từ ba nguyên tố hydro, cacbon và oxy (ký hiệu dạng $H_x C_y O_z$), trong đó tỷ lệ cacbon và oxy chiếm phần lớn khối lượng. Ở dải năng lượng photon cao (661,66 keV của nguồn Cs-137), cơ chế tương tác chiếm ưu thế là tán xạ Compton trên electron tự do, tiết diện tương tác tỷ dụng chủ yếu với tỷ số $Z/A$ (đại lượng này đối với cacbon và oxy là gần như tương đương nhau). Dữ liệu NIST cho thấy hệ số suy giảm khối của 14 loại gỗ dao động rất hẹp quanh 0,077 - 0,082 cm²/g, với độ lệch so với trung bình của 12/14 loại gỗ đều dưới 1%, cho phép xem hệ số này là một hằng số trong phương trình nội suy.

Độ chính xác nội suy mật độ trên các mẫu gỗ thử nghiệm đạt mức độ nào?

Kết quả kiểm tra trên các mẫu gỗ kiểm chứng độc lập (Fir3, Pine, Oak) cho thấy độ lệch tương đối giữa mật độ nội suy từ đường chuẩn và mật độ chuẩn thực tế đạt mức rất cao: đối với mẫu hình trụ dao động từ 0,221% đến 0,319% (sai số phép đo từ 0,45% đến 0,51%); đối với mẫu hình hộp chữ nhật dao động từ 0,245% đến 0,589%. Toàn bộ các kết quả kiểm chứng đều có độ lệch tương đối nhỏ hơn 1%, chứng minh tính hiệu quả và độ tin cậy vượt trội của đường chuẩn xây dựng từ mô phỏng MCNP6.

Kết luận

Khóa luận tốt nghiệp của tác giả Lê Hồng Thiện đã ứng dụng thành công chương trình Monte Carlo MCNP6 kết hợp phần mềm xử lý phổ Colegram để thiết lập đường chuẩn xác định mật độ gỗ bằng phương pháp gamma truyền qua với nguồn Cs-137 và đầu dò NaI(Tl). Nghiên cứu đã xây dựng hoàn chỉnh hai đường chuẩn có hệ số tương quan tuyến tính cao ($R^2 = 0,9996$) cho hai dạng hình học mẫu hộp chữ nhật và mẫu trụ tròn trong dải mật độ từ 0,405 g/cm³ đến 1,110 g/cm³. Kết quả kiểm chứng trên các mẫu gỗ độc lập khẳng định độ chính xác của phương pháp với độ lệch so với mật độ chuẩn dưới 1% và sai số thống kê cực đại chỉ 0,62%. Công trình cung cấp một quy trình tính toán bán thực nghiệm chặt chẽ, mở ra hướng ứng dụng hiệu quả kỹ thuật kiểm tra không phá hủy trong nghiên cứu vật lý hạt nhân và công nghệ vật liệu gỗ.

Trích đoạn nội dung tài liệu

Chương 1. Trình bày về cơ chế tương tác bức xạ gamma trong vật chất là hiệu ứng quang điện, hiệu ứng Compton và hiệu ứng tạo cặp. Lý thuyết về phương pháp gamma truyền qua và tính sai số của phép đo. Trình bày về chương trình mô phỏng Monte Carlo.

Trong chương này, trình bày về cách khai báo nguồn. chương trình Colegram. xử lý mô phỏng. Kết quả và thảo luận.

Trình bày kết quả xử lý với hai kích thước là gỗ dạng hình hộp chữ nhật và gỗ hình trụ tròn có cùng bê dày là 5 em. Đánh giá vả so sánh kết quả, xây dựng đường chuẩn mô phỏng với hai mẫu gỗ, độ tin cậy phép đo. Ứng dụng của phép đo trong thực tế. Tông quan về bức xạ gamma tương tác với vật chất Trong chương này, trình bày các cơ chế tương tác bức xạ gamma với vật chất dé khảo sát các bài toán nghiên cứu và ứng dụng vật lý hạt nhân trong công nghiệp.

môi trường và lĩnh vực kinh tế. Tương tác của bức xạ gamma với môi trường vật chất sẽ khác nhau với mỗi mức năng lượng. Tia gamma tương tác với môi trường chủ yếu là tương tác điện từ. Bức xạ gamma không mang điện nên không chịu ảnh hưởng của trưởng Coulomb và nó chỉ tương tác với electron trong vùng có bán kính khoảng 10 !! m.

Cường độ của chùm tia gamma giảm dan do hap thụ hay tán xạ với vật chất. Bức xạ gamma là bức xạ điện từ có bước sóng cực ngắn, với bước sóng thấp nhật khoảng 10 -!? m và có khả năng ion hóa mạnh trong môi trường vật chất. Trên Trái Dat bức xạ gamma thường được sinh ra bởi sự phân rã của các hạt nhân không bên và các bức xạ thứ cấp từ các tương tác của các hạt trong tia vũ trụ. Các hiệu ứng của bức xạ đi qua môi trường vật chất là tán xạ Compton, hiệu ứng tạo cặp.

hiệu ứng quang điện. Ngoài ra, chương | trình bày quy luật suy giảm cường độ bức xạ gamma trong môi trường vật chất. Bức xạ ion hóa đi vào trong vật chất sẽ truyền năng lượng, kích thích và ion nguyên tử xung quanh. Kích thích là quá trình nguyên tử hấp thụ năng lượng và chuyển lên trạng thái kích thích.

O trạng thái nguyên tử, phân tử phát ra photon năng lượng, bức xạ nhiệt hay phan ứng hóa học đề trở lại trạng thái cơ bản. Quá trình nguyên tử, phân tử bị kích thích thì không có electron nào bị bức ra khỏi nguyên tử, phân tử. lon hóa là quá trình electron bị bức ra khỏi nguyên tử, phân tử do nó nhận được năng lượng lớn hơn năng lượng liên kết với nguyên tử, phân tử. Quá trình ion hóa tạo ra một cặp ion là ion âm (electron thoát ra) và ion đương (phan còn lại của nguyên tử, phân tr).

Hiệu ứng quang điện Hiệu ứng quang điện là dang tương tác phô biến khi năng lượng của tia gamma có năng lượng lớn hơn năng lượng liên kết clectron nguyên tử với hạt nhân. Tia gamma tương tác với electron phá vỡ liên kết của nó và nguyên tử, năng lượng của nó truyền cho electron và nguyên tử bị kích thích. Các electron sẽ hap thụ hoàn toàn năng lượng photon đến đẻ phá vỡ năng lượng liên kết và thoát ra khỏi về bề mặt kim loại với động năng cực đại. Albert Einstein dựa vào lý thuyết lượng tử ánh sáng của Max Planck đã giải thích thành công hiệu ứng quang điện.

Năng lượng của photon tới ít nhất phải bằng năng lượng lượng liên kết electron thì hiệu ứng quang điện mới xảy ra. Hiệu ứng quang điện không xảy ra với electron tự đo. © Electron thoát ra hd Hình 1. Hiệu ứng quang điện 1.

Tán xạ Compton Electron tin xg ® .2 Tan xa Compton Tan xạ Compton là hiện tượng tin xạ đản hoi giữa photon va electron tự đo. Khi xảy ra tán xạ Compton thì photon sẽ mắt đi một phần năng lượng và sẽ bị lệch phương so với ban đầu. Hiệu ứng tạo cặp Hiệu ứng tạo cặp xảy ra khi năng lượng của tia gamma lớn hơn hai lần năng lượng nghỉ của electron (1022 keV) lan truyền trong môi trường vật chất, vùng lân cận của hạt nhân nguyên tử thì có thê xảy ra hiệu ứng tạo cặp electron và positron. Dé xảy ra quá trình tạo cặp thì năng lượng tới của photon phải lớn hon năng lượng nghỉ của hai hạt và tuần theo định luật bảo toàn năng lượng và định luật bảo toàn động lượng.

Đôi với các photon có năng lượng lớn hơn 1022 keV, quá trình tạo cặp là phương thức chiếm ưu thể của tương tác photon với vật chất. Những tương tác này lần đầu tiên được quan sát thay trong buông mây được tim thay bởi Patrick Blackett, ông được trao giải Nobel Vật lý năm 1948. Suy giảm tia gamma Sự suy giảm của một chùm tia gamma hẹp đi qua một lớp vật liệu đồng nhất được đặc trưng bởi hệ số suy giảm. Hệ số này càng cao thì suy giảm càng nhanh.

Hệ số suy giảm là thước đo xác suất tương tác của tia gamma với vật chất, tương tác với vật chất càng nhiều thì năng lượng mat mát càng lớn. Các hiệu ứng như hiệu ứng quang điện, hiệu ứng Compton va hiệu ứng tạo cặp đồng thời cùng đóng góp vào hệ số suy giảm. Thường có một xác suất đáng kẻ là tia gamma truyền đi một khoảng cách tương đối lớn mà không tham gia vào bat kỳ tương tác nào. Kết qua là bức xạ gamma vừa xuyên thấu hơn vừa khó dừng hơn so với bức xa alpha hoặc beta.

Xác suất tương tác của tia gamma với hạt nhân có số Z lớn thì cảng cao. LkeV 1 MeV 1keV 1 MeV Nang lượng Photon Hình 1. Suy giảm tia gamma trong vật liệu cacbon va chi [18] Trong Hình 1.4, đối với vật liệu có số nguyên tứ Z lớn như chi thì suy giảm chùm tia gamma sẽ lớn hơn là nguyên tử có số Z nhỏ như cacbon. Bởi vì những nguyên tử có Z càng lớn thì lớp điện tử xung quanh hạt nhan sẽ dày đặc hơn so với những nguyên tử có Z nhỏ.

Tiết diện tán xa Rayleigh và tiết diện hap thụ quang điện trong cacbon chỉ xảy ra ở mức năng lượng 20 keV. Trong khi đó, đối với hạt nhân chì có tiết điện tán xạ Rayleigh xảy ra ở vùng năng lượng 600 keV. Khi tia gamma đi vào môi trường vật chất. khả năng tương tác với lớp điện tử xung quanh hạt nhân có số Z lớn sẽ cao hơn nên chùm tia gamma khi đi vào trong vật liệu có số Z lớn sẽ bị hap thụ nhiều hơn.

Phương pháp gamma truyền qua 1. Hệ số đóng góp B Khi chùm bức xạ gamma đi qua vật chất thì bức xạ gamma sẽ giảm theo công thức (1.3) với điều kiện hình học tốt như bẻ dày của chất hap thụ có bé day mỏng và các chùm bức xạ song song. I = lạ e~t#mPX (13) Trong một số điều kiện hình học có bê dày vật liệu hap thụ tương đối lớn và các bức xạ gamma phân kì thì công thức (1.3) không còn phù hợp nữa, do đó cần phải chỉnh sửa bang cách thêm vào công thức trên một hệ số đóng góp B [12]: = In Be tm?P* (1. Hé sé đóng góp phụ thuộc vào năng lượng của bức xa gamma và vật liệu.

Khi bức xạ gamma đến vật liệu bia sẽ gây ra các tương tác như tán xạ compton và tán xạ Rayleigh và sinh ra các bức xạ gamma thứ cấp. Các bức xạ thứ cấp này đi đến đầu dd và không thê phân biệt so với bức xạ gamma bị suy giảm. Các nguyên nhân sinh ra tán xạ thứ cấp: © Tan xạ trong vật liệu bia; © Tan xạ trong không khí; 3° Tán xạ với vật liệu chế tạo đầu dò; © Tan xạ với vat liệu chê tạo nguồn và vật liệu che chan. Vì hệ số đóng góp nảy phức tạp nên trong đo thực nghiệm thường ưu tiên việc sử dụng chùm tia song song và hẹp.

Dé tạo ra chùm tia song song và hẹp thì cần phải chuẩn trực bằng chì có đường kính trong nhỏ và bê day đủ lớn. Ngoài ra, để có thể hạn chế thành phan tán xạ nhiều lần đến đầu dò thì ta sử dụng ống chuẩn trực dau dò. Trong bài toán mô phóng, chúng tôi sử đụng hai khối chuẩn trực nguồn và chuẩn trực đầu đò có đường kính trong lần lượt 0,251 em và 0,502 em. Do đóng góp của thành phan thứ cấp nhỏ hơn nhiều so với trước khi chưa đóng góp nên xem hệ số tích lũy bằng 1.

Phương pháp gamma truyền qua Khi một chùm tia photon được phát ra từ đồng vị phóng xạ đi vào đầu dò, bức xạ gamma tương tac với các clectron bên trong vật liệu, một số photon có khả năng bị hap thụ hoặc bị tán xạ với các electron trong môi trường. Dựa vào sự suy giảm chùm tia photon khi đi qua mẫu, phương pháp gamma truyền qua được sứ dụng dé tim môi liên hệ cường độ chùm tia với mật độ. Mỗi liên hệ giữa số đếm được ghi nhận từ dau đò được bỗ trí thí nghiệm ở Hình 1.5 với khi không có mẫu là No và khi có mẫu là N, tuân theo công thức 1.5 [11]: Chuan tri dav db Tinh thé Nal(Th) ¥ Vật liệu Ông nhân qaang điện Võ ahéen Hình 1. Mô hình phương pháp gamma truyền qua N = Np ented (1.5) trong đó; d là bề dày của gỗ (cm); Ht, là hệ số hap thụ tuyến tính (em? g}), Hệ số suy giảm tuyến tính thay đôi theo thành phan vật chất và mật độ.

Sự phụ thuộc của nguyên tử Z chiếm ưu thế va sự thay đối thành phần hóa học của nguyên liệu gây ra sự thay đôi của hệ số suy giảm tuyến tính ngay cả đối với các loại vật liệu có cùng khối lượng riêng. Ở năng lượng cao thì sự phụ thuộc của hệ số suy giảm tuyến tính vào mật độ là chiếm ưu thé hơn. Hệ số suy giảm tuyến tính được tính như sau: H, 2Í X.6 Đặc tính của hệ số suy giảm tuyến tính phải được xác định độc lập với tỉ trọng của nó. Đối với môi trường xốp như gỗ, hệ số suy giảm tuyến tính là một hàm của khối lượng riêng và hàm lượng nước trong vật liệu đó theo phương trình [S]: Hy = Ha Py FO bine (1.7) trong do: It,, là hệ số suy giảm khối của gỗ: Umw hệ số suy giảm khối của nước; p, (gem”) là mật độ tương ứng mỗi loại gỗ: @ là hàm lượng nước có trong gỗ.

10 Trong bài toán mô phỏng, mẫu gỗ đang xét được say khô và không lẫn nước nên 0=0. Phương trình (1.