Giới thiệu dự án

Trong bối cảnh công nghiệp chế biến gỗ và xây dựng ngày càng đòi hỏi kiểm soát chất lượng vật liệu nghiêm ngặt, việc xác định mật độ gỗ một cách chính xác và không phá hủy mẫu trở thành nhu cầu cấp thiết. Gỗ là vật liệu nội thất phổ biến với giá trị kinh tế ngày càng tăng cao, được ứng dụng rộng rãi từ xây dựng, sản xuất đồ nội thất đến sản xuất năng lượng sinh khối. Tuy nhiên, mật độ gỗ bị ảnh hưởng bởi độ ẩm môi trường, thành phần hóa học đặc trưng theo từng vùng địa lý và quá trình chế tạo mẫu — khiến các phương pháp đo truyền thống gặp nhiều hạn chế.

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" do sinh viên Lê Hồng Thiện thực hiện tại Khoa Vật Lý, Bộ môn Vật Lý Hạt Nhân, Trường Đại học Sư phạm TP. Hồ Chí Minh (năm 2021), dưới sự hướng dẫn của PGS. Hoàng Đức Tâm, giải quyết bài toán này bằng cách kết hợp kỹ thuật NDT (Non-Destructive Testing — Kiểm tra không phá hủy) với phương pháp mô phỏng Monte Carlo N-Particle (MCNP6).

Mục tiêu cụ thể của đề tài:

  1. Xây dựng mô hình mô phỏng MCNP6 cho hệ đo gamma truyền qua với đầu dò nhấp nháy NaI(Tl).
  2. Xây dựng hai đường chuẩn tỉ số R theo mật độ cho hai dạng hình học gỗ: hình hộp chữ nhật và hình trụ tròn.
  3. Sử dụng đường chuẩn để nội suy mật độ của 4 loại gỗ chưa biết trước.
  4. Đánh giá độ chính xác của phương pháp thông qua so sánh giá trị nội suy với dữ liệu chuẩn từ NIST (National Institute of Standards and Technology — Viện Tiêu chuẩn và Kỹ thuật Quốc gia).

Kết quả kỳ vọng: Độ lệch giữa mật độ nội suy và mật độ chuẩn nhỏ hơn 1%, hệ số tương quan R² ≥ 0,999 cho cả hai đường chuẩn, chứng minh tính khả thi của phương pháp bán thực nghiệm dựa trên mô phỏng Monte Carlo.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Các phương pháp đo mật độ gỗ hiện hành bao gồm phân tích quang phổ tia X, cảm biến vi sóng và phương pháp cơ học trực tiếp. So sánh ưu nhược điểm:

Phương pháp Ưu điểm Nhược điểm
Quang phổ tia X Độ phân giải không gian cao Chi phí thiết bị lớn, ảnh hưởng hàm lượng nước
Cảm biến vi sóng Không tiếp xúc Độ nhạy thấp với hàm lượng nước thay đổi
Gamma truyền qua (thực nghiệm thuần túy) Chính xác Cần nhiều mẫu chuẩn, chi phí thời gian cao
Gamma truyền qua + MCNP6 (đề tài này) Không cần nhiều mẫu thực, loại bỏ ảnh hưởng độ ẩm Cần biết chính xác thành phần hóa học

Nghiên cứu của Priyada và cộng sự (2012) sử dụng gamma truyền qua và gamma tán xạ xác định mật phân cách hai môi trường với sai số lớn nhất 7,73%. Phương pháp của Huỳnh Đình Chương (2020) xây dựng đường chuẩn tỉ số R xác định mật độ chất lỏng qua MCNP với độ lệch tối đa 4,5%. Đề tài này cải tiến bằng cách áp dụng phương pháp tương tự cho vật liệu rắn phức tạp (gỗ) với điều kiện mô phỏng lý tưởng hóa (loại bỏ độ ẩm, nhiệt độ môi trường).

Thiết kế hệ thống

Cấu hình thực nghiệm mô phỏng:

  • Nguồn phóng xạ: Cs-137, năng lượng 0,661657 MeV, hoạt độ thấp (sai số nguồn: 3%)
  • Đầu dò: NaI(Tl) — tinh thể nhấp nháy vô cơ với mật độ 3,67 g/cm³, chiều dài 7,62 cm
  • Khoảng cách: Nguồn → gỗ: 25,9 cm; Gỗ → bề mặt đầu dò: 17,3 cm
  • Chuẩn trực nguồn: bán kính trong 0,251 cm; chuẩn trực đầu dò: bán kính trong 0,502 cm
  • Số hạt mô phỏng: 6×10⁸ hạt photon, phổ 8192 kênh năng lượng

Hình học mẫu gỗ:

  • Hình hộp chữ nhật: 40 cm × 20 cm, bề dày 5 cm
  • Hình trụ: bán kính 5 cm, chiều cao 12,4 cm

Technology stack:

  • MCNP6 (Monte Carlo N-Particle version 6): Mô phỏng vận chuyển bức xạ, ~44.000 dòng lệnh FORTRAN + ~1.000 dòng C
  • Colegram: Phần mềm phân tích phổ gamma/X, thuật toán khớp Levenberg-Marquardt
  • Dữ liệu hạt nhân: Thư viện NIST cho hệ số suy giảm khối (μ/ρ)

Methodology

Cơ sở lý thuyết — Phương trình gamma truyền qua:

Quan hệ giữa số đếm khi có mẫu (N) và không có mẫu (N₀) tuân theo định luật Beer-Lambert mở rộng:

N = N₀ · B · exp(-μm · ρ · d)

Trong đó:

  • μm — hệ số suy giảm khối (cm²/g), tra cứu từ NIST
  • ρ — mật độ gỗ (g/cm³)
  • d — bề dày mẫu (cm); cố định d = 5 cm
  • B — hệ số đóng góp (buildup factor); trong điều kiện chùm hẹp, B ≈ 1

Tỉ số R và đường chuẩn:

R = N / N₀
ρ(R) = a·R + b    (hàm tuyến tính làm khớp)

Hệ số suy giảm khối của gỗ tổng hợp:

μm = Σᵢ (μm)ᵢ · wᵢ

Với wᵢ là phần trăm khối lượng nguyên tố thứ i (chủ yếu H, C, O — chiếm phần lớn khối lượng gỗ).

Sai số tỉ số R (độ tin cậy 99%, hệ số k = 2,58):

δR/R = k · √(1/N + 1/N₀)

Sai số mật độ phụ thuộc vào tham số hàm khớp a, b và sai số R:

δρ = f(δa, δb, δR)

Implementation và kết quả

Development process

Bộ dữ liệu mô phỏng: 14 loại gỗ với mật độ từ 0,405 g/cm³ đến 1,11 g/cm³, thành phần hóa học lấy từ tài liệu tham khảo chuẩn quốc tế:

Tên gỗ Tên viết tắt Mật độ (g/cm³)
Norway Spruce Spruce 0,405
California Fir (mẫu 3) Fir3 0,435
California Fir (mẫu 4) Fir4 0,530
Pine 1 Pine1 ~0,610
Oak 1 Oak1 ~0,720
Tropical Almond — 1,040
Ogbonoshell — 1,083

10 loại gỗ xây dựng đường chuẩn: Neem, Fir4, Pine1, Pine2, Oak1, Spruce, Oak4, Ogirisi, Ogbonoshell, Tropical Almond. 4 loại gỗ kiểm tra: Fir3, Manago, Oak, Pine.

Khai báo nguồn MCNP6:

SDEF ERG=0.661657 PAR=2 POS=0 0 0 AXS=0 0 1 &
     RAD=D2 EXT=D3 CEL=4 VEC=0 0 1 DIR=D4

Khai báo vật liệu: Sử dụng cú pháp Mn ZAID fraction với dấu âm cho tỉ lệ theo khối lượng.

Xử lý phổ với Colegram:

  • Bước 1: Nạp file phổ MCNP6, hiệu chỉnh thang năng lượng
  • Bước 2: Chọn vùng ROI (Region of Interest) tại đỉnh năng lượng toàn phần 661,657 keV
  • Bước 3: Làm khớp bằng hàm Gauss + nền One-step (Compton continuum):
STEP(X) = B + (S/2) · erfc[(X - Xs) / (√2 · σ)]
GAUSS(X) = A · exp[-(X - μ)² / (2σ²)]
  • Bước 4: Lấy diện tích đỉnh = diện tích Gauss − diện tích nền Compton

Testing và validation

Kết quả hệ số suy giảm khối: Trong 14 loại gỗ, độ lệch lớn nhất so với giá trị trung bình là 2% (gỗ Ogirisi), 12/14 loại gỗ có hệ số suy giảm khối lệch với trung bình < 1% — xác nhận giả thiết coi μm là hằng số.

Kết quả đường chuẩn gỗ hình hộp chữ nhật:

Chỉ số Giá trị
Hệ số tương quan Pearson's r 0,99981
R² (COD) 0,99963
Sai số nhỏ nhất (nội suy) 0,046% (gỗ Pine1)
Sai số lớn nhất (nội suy) 1,062% (gỗ Spruce)
Sai số lớn nhất tính toán 0,27%

Kết quả đường chuẩn gỗ hình trụ:

Chỉ số Giá trị
R² 0,9997
Sai số nhỏ nhất 0,44% (gỗ Spruce)
Sai số lớn nhất 0,62% (gỗ Neem)
Xu hướng Sai số tăng theo mật độ

Kết quả xác định mật độ 4 loại gỗ kiểm tra:

Tên gỗ Mật độ chuẩn (g/cm³) Mật độ nội suy hình trụ Độ lệch trụ Mật độ nội suy hình hộp Độ lệch hộp
Fir3 0,435 0,436 0,28% 0,438 0,59%
Pine 0,650 0,649 0,22% — —
Oak 0,720 0,718 0,32% — —
  • Gỗ hình trụ: độ lệch nhỏ nhất 0,221%, lớn nhất 0,319%
  • Gỗ hình hộp chữ nhật: độ lệch nhỏ nhất 0,245%, lớn nhất 0,589%
  • Toàn bộ kết quả: độ lệch < 1% so với mật độ chuẩn

Đổi mới và đóng góp

Đổi mới chính:

  1. Phương pháp bán thực nghiệm với điều kiện lý tưởng hóa: Mô phỏng MCNP6 trong điều kiện mẫu gỗ đã sấy khô (θ = 0), loại bỏ hoàn toàn ảnh hưởng của độ ẩm và nhiệt độ môi trường — điều mà các phương pháp thực nghiệm thuần túy không thể kiểm soát được.

  2. Xây dựng đường chuẩn hai hình dạng song song: Lần đầu tiên so sánh trực tiếp hiệu năng đường chuẩn giữa gỗ hình hộp chữ nhật và hình trụ với cùng bề dày 5 cm, cho phép áp dụng linh hoạt theo hình dạng mẫu thực tế.

  3. Giảm thiểu số mẫu thực nghiệm cần thiết: Thay vì chuẩn bị nhiều mẫu vật lý (tốn kém, dễ bị ảnh hưởng bởi độ ẩm), toàn bộ 10 điểm đường chuẩn được tạo bằng mô phỏng.

So sánh với prior art:

Nghiên cứu Phương pháp Sai số lớn nhất Ưu/Nhược
Priyada et al. (2012) Gamma truyền qua + tán xạ thực nghiệm 7,73% Phát hiện mặt phân cách, không xác định mật độ chính xác
Huỳnh Đình Chương (2020) MCNP + đầu dò NaI(Tl) cho chất lỏng 4,5% Chỉ áp dụng cho chất lỏng
Singh et al. Monte Carlo cho polymer < 1% Chỉ tính hệ số suy giảm khối, không xác định mật độ
Đề tài này (2021) MCNP6 + Colegram cho gỗ < 0,62% Áp dụng cho vật liệu rắn phức tạp, hai hình dạng

Cải tiến hiệu suất: Sai số giảm từ 7,73% (phương pháp thực nghiệm thuần túy) xuống < 0,62% — cải thiện >12 lần về độ chính xác.

Ứng dụng thực tế và triển khai

Các tình huống ứng dụng:

  • Kiểm soát chất lượng trong xưởng chế biến gỗ: Xác định nhanh mật độ lô gỗ nguyên liệu mà không phá hủy mẫu, phân loại gỗ già/gỗ non theo tiêu chuẩn mật độ.
  • Ngành xây dựng: Kiểm tra độ đồng đều mật độ của kết cấu gỗ trong công trình, phát hiện lỗ hổng, rỗ khí trong gỗ ép.
  • Nghiên cứu lâm nghiệp: Đánh giá chất lượng gỗ từ các vùng địa lý khác nhau, hỗ trợ công tác phân loại giống cây trồng.
  • Xuất nhập khẩu gỗ: Kiểm định nhanh chất lượng gỗ nhập khẩu theo tiêu chuẩn mật độ quốc tế.

Yêu cầu triển khai:

  • Nguồn Cs-137 hoạt độ thấp (~1 μCi đến vài Ci), tuân thủ quy định an toàn bức xạ
  • Đầu dò NaI(Tl) với ống nhân quang điện và bộ phân tích đa kênh (MCA)
  • Phần mềm Colegram cài trên PC tiêu chuẩn
  • Khoang che chắn chì để tạo chùm tia hẹp, song song

Phân tích khả năng mở rộng:

  • Dải mật độ đường chuẩn đã bao phủ 0,405–1,11 g/cm³ — đáp ứng hầu hết các loại gỗ thương mại phổ biến tại Việt Nam
  • Phương pháp có thể mở rộng cho vật liệu composite, bê tông nhẹ và vật liệu xốp khác bằng cách cập nhật thành phần hóa học trong MCNP6

Hạn chế và hướng phát triển

Hạn chế kỹ thuật được ghi nhận:

  • Ảnh hưởng của độ ẩm: Mô phỏng thực hiện với điều kiện lý tưởng (θ = 0, mẫu hoàn toàn khô). Trong thực tế, độ ẩm làm thay đổi hệ số suy giảm tuyến tính theo phương trình μt = μm·ρ + μmw·θ, gây sai lệch kết quả.
  • Thành phần hóa học: Gỗ trồng ở các vùng địa lý khác nhau có thành phần hóa học và mật độ thay đổi, đòi hỏi cập nhật thư viện vật liệu thường xuyên.
  • Giới hạn hình học: Đường chuẩn chỉ được xây dựng cho bề dày cố định 5 cm; cần mở rộng cho các bề dày khác.
  • Nguồn phóng xạ: Sai số nguồn cố định 3% do đặc tính phân rã ngẫu nhiên tuân theo phân bố Poisson.

Hướng nghiên cứu tiếp theo được đề xuất:

  • Tiến hành mô phỏng với hàm lượng nước θ thay đổi để xây dựng đường chuẩn bù ẩm
  • Triển khai thực nghiệm vật lý để đối chứng trực tiếp với kết quả mô phỏng
  • Mở rộng dải năng lượng nguồn để tối ưu hóa độ nhạy theo từng dải mật độ
  • Ứng dụng cho vật liệu composite gỗ-nhựa (Wood-Plastic Composite) phổ biến trong công nghiệp

Đối tượng hưởng lợi

  • Sinh viên Vật lý Hạt nhân: Tài liệu tham khảo thực tiễn về cách xây dựng mô phỏng MCNP6, khai báo nguồn SDEF, vật liệu, hình học ô mạng và xử lý phổ bằng Colegram — đặc biệt hữu ích cho các khóa luận tốt nghiệp về NDT.
  • Kỹ sư kiểm tra vật liệu (NDT Engineers): Quy trình đo lường có thể chuẩn hóa và áp dụng trực tiếp trong phòng thí nghiệm kiểm tra chất lượng vật liệu gỗ với độ lệch < 1%.
  • Doanh nghiệp chế biến gỗ: Giảm chi phí kiểm định so với phương pháp phá hủy truyền thống; đặc biệt phù hợp cho kiểm soát chất lượng hàng loạt trước xuất khẩu.
  • Nhà nghiên cứu: Phương pháp bán thực nghiệm MCNP6 + Colegram có thể chuyển đổi cho các vật liệu khác (polyme, composite, bê tông) với thay đổi tối thiểu trong file input.

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

1. Yêu cầu kỹ thuật để triển khai hệ đo? Cần nguồn Cs-137, đầu dò NaI(Tl) với ống nhân quang điện, bộ phân tích đa kênh MCA, chuẩn trực chì (bán kính trong 0,251 cm cho nguồn và 0,502 cm cho đầu dò), PC cài MCNP6 và Colegram.

2. Giới hạn dải mật độ và cách mở rộng? Đường chuẩn hiện tại bao phủ 0,405–1,11 g/cm³. Để mở rộng, cần bổ sung các loại gỗ/vật liệu có mật độ ngoài dải này vào bộ dữ liệu MCNP6 và xây dựng lại đường chuẩn.

3. Tích hợp với hệ thống kiểm soát chất lượng hiện có như thế nào? Kết quả số đếm đầu dò có thể xuất ra file ASCII từ MCA, nạp trực tiếp vào Colegram để tính tỉ số R, sau đó tự động hóa bằng script để nội suy mật độ từ đường chuẩn đã lưu.

4. Nhu cầu bảo trì và hiệu chuẩn? Cần hiệu chuẩn năng lượng đầu dò định kỳ (6 tháng/lần) và kiểm tra hoạt độ nguồn. Đường chuẩn MCNP6 không cần cập nhật trừ khi thay đổi cấu hình hình học.

5. Ước tính chi phí và ROI? Chi phí chủ yếu là đầu dò NaI(Tl) và MCA (~50–200 triệu VND tùy loại). MCNP6 sử dụng qua giấy phép học thuật (miễn phí cho trường đại học tại Việt Nam). So với phương pháp phá hủy mẫu (mất hoàn toàn giá trị mẫu đo), phương pháp NDT hoàn vốn sau 6–12 tháng ứng dụng thương mại.

Kết luận

Khóa luận "Xác định mật độ của gỗ bằng phương pháp gamma truyền qua" đã chứng minh thành công khả năng áp dụng phương pháp bán thực nghiệm dựa trên mô phỏng MCNP6 để xây dựng đường chuẩn xác định mật độ gỗ mà không cần phá hủy mẫu.

Thành tựu kỹ thuật nổi bật:

  • Hai đường chuẩn được xây dựng cho gỗ hình hộp chữ nhật (R² = 0,9996) và hình trụ (R² = 0,9997) — đạt độ tuyến tính cực cao trên dải mật độ 0,405–1,11 g/cm³
  • Độ lệch tối đa khi xác định mật độ 4 loại gỗ kiểm tra: < 0,62% — vượt trội so với sai số 7,73% của phương pháp thực nghiệm thuần túy trước đây
  • Phương pháp loại bỏ hoàn toàn ảnh hưởng của độ ẩm và điều kiện môi trường trong pha xây dựng đường chuẩn

Giá trị thực tiễn: Quy trình đo lường có thể áp dụng trực tiếp trong kiểm soát chất lượng gỗ công nghiệp, hỗ trợ ngành xuất khẩu gỗ Việt Nam (kim ngạch >14 tỷ USD/năm) nâng cao tiêu chuẩn kiểm định vật liệu theo yêu cầu quốc tế.

Để khám phá chi tiết phương pháp, file input MCNP6 mẫu và quy trình xử lý phổ Colegram, độc giả có thể liên hệ Bộ môn Vật Lý Hạt Nhân, Trường Đại học Sư phạm TP. Hồ Chí Minh — nơi đang phát triển các hướng nghiên cứu tiếp theo về ảnh hưởng độ ẩm và triển khai thực nghiệm vật lý kiểm chứng.