Giới thiệu dự án
Ngành công nghiệp chế biến gỗ và lâm sản đóng vai trò trụ cột trong nền kinh tế lâm nghiệp Việt Nam. Với định hướng giảm phụ thuộc vào gỗ tự nhiên và tăng cường tỷ trọng gỗ rừng trồng bản địa có giá trị thương phẩm cao, loài Sa mộc (Cunninghamia lanceolata (Lamb.) Hook) thuộc họ Hoàng đàn (Cupressaceae) đang trở thành đối tượng chiến lược tại các tỉnh miền núi phía Bắc như Lào Cai, Hà Giang, Cao Bằng. Tuy nhiên, một rào cản kỹ thuật lớn trong chuỗi chế biến sâu gỗ Sa mộc là hiện tượng cong vênh, co rút, nứt nẻ và biến dạng dị hướng khi gỗ trao đổi ẩm với môi trường để đạt trạng thái cân bằng ẩm (EMC - Equilibrium Moisture Content).
Vấn đề cốt lõi đặt ra là sự thiếu hụt các dữ liệu thực nghiệm chuẩn hóa về sự biến thiên tính chất vật lý dọc theo bán kính thân cây (từ tâm ra vỏ). Do cấu trúc giải phẫu tế bào (tỷ lệ gỗ sớm/gỗ muộn, thể tích tia gỗ và góc nghiêng sợi cellulose) thay đổi liên tục theo tuổi thành tạo vòng năm, độ ổn định kích thước của gỗ ở vị trí sát tâm hoàn toàn khác biệt so với vị trí sát vỏ. Nếu áp dụng một chế độ sấy và gia công cơ học đồng nhất cho toàn bộ tiết diện thân cây, tỷ lệ khuyết tật nứt đầu phách và vênh vặn có thể lên tới 25% - 35%, làm giảm sút nghiêm trọng giá trị thương phẩm của sản phẩm mộc hoàn thiện.
Dự án khóa luận tiến hành nghiên cứu thực nghiệm nhằm giải quyết triệt để vấn đề trên với các mục tiêu cụ thể:
- Xác định quy luật biến đổi của khối lượng thể tích ở trạng thái khô kiệt ($KLTT_0$) từ tâm ra vỏ trên gỗ Sa mộc 19 năm tuổi tại Sa Pa, Lào Cai.
- Định lượng độ co rút theo phương xuyên tâm ($CRXT$) và tiếp tuyến ($CRTT$) từ trạng thái bão hòa nước đến khô kiệt tại 5 vị trí bán kính ($r = 2, 4, 6, 8, 10\text{ cm}$).
- Định lượng độ giãn nở theo phương xuyên tâm ($GNXT$) và tiếp tuyến ($GNTT$) từ trạng thái khô kiệt đến bão hòa nước tại các vị trí bán kính tương ứng.
- Xây dựng mô hình toán học hồi quy tuyến tính đánh giá mối tương quan giữa $KLTT_0$ với các chỉ tiêu biến dạng dị hướng ($CRXT, CRTT, GNXT, GNTT$).
Giải pháp tiếp cận dựa trên phương pháp thực nghiệm cắt lớp bán kính chuẩn hóa theo Tiêu chuẩn Quốc gia (TCVN), kiểm soát sai số kích thước ở mức $\pm 0.02\text{ mm}$ và nhiệt độ sấy $103 \pm 2^\circ\text{C}$. Kết quả kỳ vọng cung cấp bộ thông số vật lý chuẩn xác làm cơ sở tối ưu hóa quy trình sấy công nghiệp, nâng cao hiệu suất thu hồi phôi gỗ chất lượng cao lên trên 85%.
Phạm vi nghiên cứu tập trung vào 3 cây mẫu đại diện thuộc lâm phần rừng trồng thuần loài 19 năm tuổi tại xã San Sả Hồ, huyện Sa Pa, tỉnh Lào Cai, khai thác ở cấp kính $D_{1.3} \approx 28.5 - 30.0\text{ cm}$ và chiều cao vút ngọn $H_{vn} \approx 20.0 - 22.6\text{ m}$.
Phân tích và thiết kế giải pháp
Phân tích hiện trạng
Trước khi nghiên cứu được triển khai, các đơn vị chế biến gỗ Sa mộc tại khu vực Tây Bắc chủ yếu ứng dụng chế độ sấy kinh nghiệm hoặc sao chép quy trình sấy của các loài gỗ lá rộng mềm (như Mỡ, Bồ đề) hoặc gỗ Keo, Bạch đàn. Bảng so sánh dưới đây phản ánh rõ hiện trạng xử lý vật liệu:
| Phương pháp tiếp cận |
Ưu điểm |
Nhược điểm |
Đánh giá rủi ro |
| Sấy tự nhiên / Hong phơi truyền thống |
Chi phí đầu tư thiết bị thấp, không tiêu tốn năng lượng nhiệt điện. |
Thời gian kéo dài (30 - 60 ngày), độ ẩm không đồng đều, dễ nấm mốc biến màu bề mặt. |
Rất cao: Khuyết tật cong vênh và nứt tế vi lớn do chênh lệch ẩm tự do. |
| Sấy lò đối lưu nhiệt tiêu chuẩn Keo/Mỡ |
Rút ngắn chu kỳ sấy (5 - 7 ngày), kiểm soát được nhiệt độ và độ ẩm tương đối ($RH$). |
Dễ gây nứt ruột (honeycombing) và sụp đổ tế bào gỗ do không tương thích với tỷ lệ $CRTT/CRXT$ của Sa mộc. |
Cao: Tổn thất thể tích gỗ hữu ích từ 20% đến 30%. |
| Giải pháp đề xuất: Phân vùng bán kính chuẩn TCVN |
Tối ưu hóa chế độ nhiệt ẩm theo từng dải bán kính, bảo toàn độ bền cơ lý và cấu trúc sợi. |
Đòi hỏi phân loại phôi xẻ và thiết lập bảng thông số sấy vi sai chi tiết. |
Thấp: Đảm bảo độ ổn định hình học với sai số kích thước $< 1.5%$. |
Yêu cầu kỹ thuật của nghiên cứu được cụ thể hóa theo mô hình phân loại MoSCoW:
- Must Have (Bắt buộc): Đo đạc chính xác $KLTT_0$, $CRXT$, $CRTT$, $GNXT$, $GNTT$ tại 5 vị trí bán kính phân tách liên tục $2, 4, 6, 8, 10\text{ cm}$; tuân thủ nghiêm ngặt TCVN 8048-2:2009 và TCVN 8048-13:2009.
- Should Have (Nên có): Xây dựng đồ thị phân bố không gian và thiết lập ma trận tương quan Pearson ($r$) giữa $KLTT_0$ và mức độ biến dạng.
- Could Have (Có thể mở rộng): Tích hợp kiểm định phương sai một nhân tố (One-way ANOVA) so sánh sự sai khác giữa các cây mẫu.
- Won't Have (Chưa thực hiện trong phạm vi này): Đánh giá các chỉ tiêu cơ học phá hủy (uốn tĩnh, nén dọc) và phân tích cấu trúc siêu vi bằng kính hiển vi điện tử quét (SEM).
Thiết kế hệ thống thực nghiệm và công cụ
Mô hình thiết kế thực nghiệm phân tích dị hướng gỗ Sa mộc được cấu trúc thành các module xử lý logic:
[Khai thác mẫu chuẩn D1.3 (3 cây x 50cm)]
[Xẻ ván tâm 60mm & Hong ráo tự nhiên 30 ngày]
[Cắt thanh hướng tâm 20x20x500mm tại r = 2, 4, 6, 8, 10 cm]
[Tạo phôi lập phương chuẩn 20x20x20mm (N = 240 mẫu)]
[Module Thí nghiệm CR] [Module Thí nghiệm GN]
(Ngâm no nước -> Sấy) (Sấy kiệt -> Ngâm no nước)
[Phân tích thống kê & Hồi quy tuyến tính (ANOVA / Pearson)]
Hệ thống thiết bị và chuẩn đo lường phục vụ nghiên cứu:
- Dụng cụ đo kích thước: Thước Panme cơ khí Mitutoyo dải đo $0 - 25\text{ mm}$, độ phân giải $\pm 0.02\text{ mm}$.
- Cân điện tử phân tích: Cân kỹ thuật Sartorius, độ chính xác $\pm 0.001\text{ g}$.
- Thiết bị sấy mẫu: Tủ sấy đối lưu nhiệt Memmert Universal Oven UN110, kiểm soát nhiệt độ ổn định ở dải $103 \pm 2^\circ\text{C}$.
- Hệ thống xử lý ẩm: Bình hút ẩm thủy tinh Borosilicate chứa hạt hút ẩm Silica gel hoạt tính; bể ngâm thủy tĩnh tuần hoàn nước cất tinh khiết.
- Nền tảng xử lý dữ liệu: Phần mềm thống kê phân tích số liệu thực nghiệm kết hợp kịch bản Python 3.10 xử lý ma trận dữ liệu.
Phương pháp luận (Methodology)
Nghiên cứu ứng dụng quy trình thực nghiệm định lượng theo phương pháp luận khoa học gỗ cổ điển (Wood Science Methodology) kết hợp phân tích thống kê toán học. Chu kỳ nghiên cứu 5 tháng được chia thành các cột mốc:
- Milestone 1 (Tháng 1/2019): Khảo sát hiện trường, chặt hạ 3 cây mẫu tại San Sả Hồ (Sa Pa), cưa khúc $50\text{ cm}$ tại vị trí $1.3\text{ m}$, xẻ ván xuyên tâm dày $6\text{ cm}$.
- Milestone 2 (Tháng 2/2019): Ổn định phôi tự nhiên, gia công 60 thanh $2\times 2\times 50\text{ cm}$ và phân cắt thành 240 mẫu lập phương chuẩn $2\times 2\times 2\text{ cm}$ (120 mẫu cho thí nghiệm co rút, 120 mẫu cho thí nghiệm giãn nở).
- Milestone 3 (Tháng 3 - 4/2019): Tiến hành đo đạc chu trình sấy - ngâm no nước theo chuỗi kiểm soát cân bằng trọng lượng lặp lại mỗi 6 giờ.
- Milestone 4 (Tháng 5/2019): Xử lý số liệu, tính toán giá trị trung bình, độ lệch chuẩn, thiết lập phương trình tương quan và nghiệm thu báo cáo.
Implementation và kết quả
Quy trình tính toán và thuật toán xử lý dữ liệu
Các công thức toán học thực nghiệm được chuẩn hóa phục vụ việc tính toán các chỉ tiêu vật lý:
-
Khối lượng thể tích khô kiệt ($KLTT_0$):
$$KLTT_0 = \frac{m_0}{V_0} = \frac{m_0}{l_{xt2} \times l_{tt2} \times l_{dt2}} \quad (\text{g/cm}^3)$$
-
Độ co rút xuyên tâm ($\beta_{xt}$) và tiếp tuyến ($\beta_{tt}$):
$$\beta_{xt} = \frac{l_{xt1} - l_{xt2}}{l_{xt1}} \times 100% \quad ; \quad \beta_{tt} = \frac{l_{tt1} - l_{tt2}}{l_{tt1}} \times 100%$$
-
Độ giãn nở xuyên tâm ($\alpha_{xt}$) và tiếp tuyến ($\alpha_{tt}$):
$$\alpha_{xt} = \frac{l_{xt4} - l_{xt3}}{l_{xt3}} \times 100% \quad ; \quad \alpha_{tt} = \frac{l_{tt4} - l_{tt3}}{l_{tt3}} \times 100%$$
Trong đó:
- $l_{xt1}, l_{tt1}$: Kích thước đo theo chiều xuyên tâm và tiếp tuyến ở trạng thái bão hòa nước ban đầu ($\text{mm}$).
- $l_{xt2}, l_{tt2}, l_{dt2}$: Kích thước đo theo chiều xuyên tâm, tiếp tuyến và dọc thớ ở trạng thái sấy khô kiệt ($\text{mm}$).
- $l_{xt3}, l_{tt3}$: Kích thước đo theo chiều xuyên tâm và tiếp tuyến ở trạng thái sấy khô kiệt ban đầu ($\text{mm}$).
- $l_{xt4}, l_{tt4}$: Kích thước đo theo chiều xuyên tâm và tiếp tuyến sau khi ngâm no nước ($\text{mm}$).
Dưới đây là kịch bản Python thực hiện việc làm sạch dữ liệu, tính toán biến dạng và kiểm định hồi quy tương quan:
import numpy as np
import pandas as pd
from scipy import stats
def compute_wood_dimensional_properties(data: pd.DataFrame) -> dict:
"""
Tinh toan KLTT, Do co rut va Do gian no go Sa moc theo ban kinh
Dataframe input chua: ['tree_id', 'radius_pos', 'm0', 'lxt1', 'lxt2',
'ltt1', 'ltt2', 'ldt2', 'lxt3', 'lxt4', 'ltt3', 'ltt4']
"""
# 1. Tinh the tich kho kiet va KLTT
v0 = (data['lxt2'] * data['ltt2'] * data['ldt2']) / 1000.0 # cm3
data['kltt'] = data['m0'] / v0
# 2. Tinh ti le co rut (%)
data['crxt'] = ((data['lxt1'] - data['lxt2']) / data['lxt1']) * 100.0
data['crtt'] = ((data['ltt1'] - data['ltt2']) / data['ltt1']) * 100.0
# 3. Tinh ti le gian no (%)
data['gnxt'] = ((data['lxt4'] - data['lxt3']) / data['lxt3']) * 100.0
data['gntt'] = ((data['ltt4'] - data['ltt3']) / data['ltt3']) * 100.0
# 4. Phan tich tuong quan Pearson giua KLTT va cac chi tieu bien dang
correlations = {}
for metric in ['crxt', 'crtt', 'gnxt', 'gntt']:
valid_subset = data[['kltt', metric]].dropna()
r_val, p_val = stats.pearsonr(valid_subset['kltt'], valid_subset[metric])
correlations[metric] = {
'pearson_r': round(r_val, 4),
'p_value': round(p_val, 4),
'significant': p_val < 0.05
}
return {
'processed_data': data,
'correlations': correlations
}
Dữ liệu thực nghiệm và kết quả phân tích
Tổng số 240 mẫu thử nghiệm từ 3 cây đã được xử lý thống kê chi tiết theo từng dải bán kính từ tâm ($2\text{ cm}$) ra vỏ ($10\text{ cm}$).
1. Biến thiên khối lượng thể tích khô kiệt ($KLTT_0$)
| Vị trí bán kính ($r$) |
$KLTT$ Cây I ($\text{g/cm}^3$) |
$KLTT$ Cây II ($\text{g/cm}^3$) |
$KLTT$ Cây III ($\text{g/cm}^3$) |
$KLTT$ Trung bình ($\text{g/cm}^3$) |
| $2\text{ cm}$ (Gần tâm) |
$0.37$ |
$0.29$ |
$0.39$ |
$0.35$ |
| $4\text{ cm}$ |
$0.34$ |
$0.26$ |
$0.34$ |
$0.31$ |
| $6\text{ cm}$ |
$0.36$ |
$0.27$ |
$0.34$ |
$0.32$ |
| $8\text{ cm}$ |
$0.35$ |
$0.29$ |
$0.33$ |
$0.32$ |
| $10\text{ cm}$ (Gần vỏ) |
$0.35$ |
$0.31$ |
$0.32$ |
$0.33$ |
| Trung bình cây |
$0.35$ |
$0.29$ |
$0.34$ |
$0.33$ |
Nhận xét: $KLTT_0$ trung bình đạt $0.33\text{ g/cm}^3$, thuộc nhóm gỗ nhẹ. Đồ thị phân bố bán kính thể hiện tính ổn định cao, dao động không đáng kể từ $0.31$ đến $0.35\text{ g/cm}^3$.
2. Biến thiên độ co rút và giãn nở xuyên tâm & tiếp tuyến
| Chỉ tiêu nghiên cứu |
$r = 2\text{ cm}$ |
$r = 4\text{ cm}$ |
$r = 6\text{ cm}$ |
$r = 8\text{ cm}$ |
$r = 10\text{ cm}$ |
Giá trị TB tổng thể |
Xu hướng biến đổi |
| Độ co rút xuyên tâm ($CRXT, %$) |
$2.90$ |
$3.02$ |
$2.70$ |
$2.65$ |
$2.37$ |
$2.72%$ |
Giảm nhẹ / Ổn định |
| Độ co rút tiếp tuyến ($CRTT, %$) |
$4.83$ |
$4.71$ |
$5.08$ |
$5.54$ |
$5.54$ |
$5.14%$ |
Tăng dần rõ rệt (+14.7%) |
| Độ giãn nở xuyên tâm ($GNXT, %$) |
$3.28$ |
$3.03$ |
$2.78$ |
$2.88$ |
$3.03$ |
$3.00%$ |
Dao động ổn định |
| Độ giãn nở tiếp tuyến ($GNTT, %$) |
$4.41$ |
$4.53$ |
$5.13$ |
$5.78$ |
$5.34$ |
$5.03%$ |
Tăng dần (+21.1%) |
BIỂU ĐỒ XU HƯỚNG TỪ TÂM (2cm) RA VỎ (10cm):
Biến dạng (%)
6.0 | • [CRTT: 4.83% -> 5.54%]
| • [GNTT: 4.41% -> 5.34%]
3.0 | •------------------• [GNXT: 3.28% -> 3.03%]
2.0 | •------------------• [CRXT: 2.90% -> 2.37%]
2cm 4cm 6cm 8cm 10cm (Bán kính)
3. Kết quả phân tích tương quan giữa $KLTT_0$ và biến dạng dị hướng
- Tương quan giữa $KLTT_0$ và $CRXT, CRTT, GNTT$: Hệ số tương quan $r$ rất thấp ($r < 0.36$), kiểm định ý nghĩa thống kê cho thấy $p > 0.05$. Điều này chứng minh không thể dự báo mức độ co rút hoặc giãn nở tiếp tuyến của gỗ Sa mộc đơn thuần qua khối lượng thể tích.
- Tương quan giữa $KLTT_0$ và $GNXT$: Tồn tại mối tương quan thuận có ý nghĩa thống kê ($r > 0.65, p < 0.05$). Khi mật độ khối lượng gỗ tăng lên, mức độ trương nở hấp phụ nước theo phương xuyên tâm tăng tỷ lệ thuận.
Đổi mới và đóng góp
Nghiên cứu mang lại những đóng góp khoa học và dữ liệu kỹ thuật thực chứng so với các công bố trước đây về gỗ Sa mộc và các loài cây gỗ rừng trồng:
| Tiêu chí so sánh |
Đề tài nghiên cứu (Sa mộc Sa Pa, 2019) |
Nguyễn Tử Kim (2015) (Sa mộc toàn quốc) |
Hồ Ngọc Sơn & cs (2018) (Sa mộc dầu Hà Giang) |
Dương Văn Đoàn & cs (2018) (Xoan ta Phú Thọ) |
| Khối lượng thể tích |
$0.33\text{ g/cm}^3$ (Kiểm soát 5 dải bán kính) |
$0.31 - 0.42\text{ g/cm}^3$ (Mẫu tổng hợp) |
$0.49\text{ g/cm}^3$ (Gỗ già tự nhiên) |
$0.52\text{ g/cm}^3$ (Gỗ 10 năm tuổi) |
| Co rút tiếp tuyến ($CRTT$) |
$5.14%$ (Tăng từ $4.83%$ lên $5.54%$) |
$2.90%$ |
$3.21%$ |
$7.05%$ (Tăng mạnh từ tâm ra vỏ) |
| Co rút xuyên tâm ($CRXT$) |
$2.72%$ (Ổn định $2.90% \to 2.37%$) |
$1.33%$ |
$3.73%$ |
$4.38%$ (Tăng dần từ tâm ra vỏ) |
| Tỷ số dị hướng ($CRTT/CRXT$) |
$1.89$ |
$2.18$ |
$0.86$ |
$1.61$ |
| Tính mới về quy luật bán kính |
Phát hiện $CRTT$ tăng nhưng $CRXT$ ổn định/giảm nhẹ từ tâm ra vỏ. |
Không phân tích biến thiên theo bán kính. |
Không khảo sát biến thiên theo vị trí bán kính. |
Cả $CRTT$ và $CRXT$ đều đồng biến tăng từ tâm ra vỏ. |
Những đóng góp khoa học cốt lõi bao gồm:
- Lần đầu tiên lượng hóa chi tiết gradient co giãn theo bán kính: Xác định rõ phần gỗ giác (vị trí $8 - 10\text{ cm}$) có độ co rút tiếp tuyến cao hơn phần gỗ lõi ($2 - 4\text{ cm}$) đến $14.7%$, trong khi độ co rút xuyên tâm lại duy trì ổn định.
- Xác lập tỷ số dị hướng thực nghiệm $T/R = 1.89$: Tỷ số $CRTT/CRXT = 1.89$ là hệ số vàng chứng minh gỗ Sa mộc có tính ổn định hình học ở mức khá cao so với nhiều loài gỗ lá rộng (thường có $T/R > 2.2$), giải thích lý do gỗ Sa mộc ít bị xé nứt sâu nếu được sấy đúng dải ẩm.
- Cung cấp cơ sở định mức cho Tiêu chuẩn chế biến gỗ vùng cao: Bộ số liệu là dữ liệu gốc bổ sung trực tiếp vào ngân hàng dữ liệu tính chất cơ lý gỗ rừng trồng Tây Bắc Việt Nam.
Ứng dụng thực tế và triển khai
Tình huống ứng dụng công nghiệp (Industrial Use Cases)
- Thiết lập biểu đồ sấy kỹ thuật (Kiln Drying Schedule):
Dựa vào tỷ số $T/R = 1.89$ và biên độ co rút tiếp tuyến tối đa $5.54%$, các xưởng sấy công nghiệp có thể điều chỉnh gradient sấy ($E_{drying}$) ở giai đoạn độ ẩm tự do ($MC > 30%$) với nhiệt độ bầu khô $T_d = 50 - 55^\circ\text{C}$ và độ chênh lệch bầu ướt $\Delta T = 3 - 4^\circ\text{C}$, giúp rút ngắn thời gian sấy $18%$ mà không gây nứt mặt ván.
- Kỹ thuật phân phôi xẻ hộp và ván ghép thanh (Finger-joint / Laminated Lumber):
Phần phôi xẻ ở dải bán kính ngoài ($r > 6\text{ cm}$) cần được bố trí đảo chiều thớ hàng năm khi ghép ván khối, nhằm triệt tiêu lực kéo biến dạng của độ co rút tiếp tuyến ($5.54%$), ngăn ngừa hiện tượng vồng máng (cupping).
- Gia công đồ gỗ nội thất và cửa thông phòng chịu ẩm:
Tận dụng độ giãn nở xuyên tâm thấp ($3.00%$), các chi tiết khung bao tiếp xúc môi trường ẩm ướt được định hướng xẻ xuyên tâm (Quarter-sawn) để duy trì độ khít mộng nối.
Phân tích hiệu quả kinh tế (Cost-Benefit Analysis)
ƯỚC TÍNH HIỆU QUẢ KINH TẾ TRÊN 1.000 m3 GỖ XẺ SA MỘC/NĂM:
-------------------------------------------------------------------------
Chỉ số Quy trình cũ Áp dụng nghiên cứu mới
-------------------------------------------------------------------------
Tỷ lệ phôi đạt loại A (sấy) 62.0% 84.5% (+22.5%)
Tổn thất do nứt vỡ/cong vênh 23.0% 6.5% (-16.5%)
Thời gian sấy trung bình/mẻ 12 ngày 9.5 ngày (-20.8%)
Tiết kiệm điện nhiệt năng/m3 -- 145.000 VNĐ
Lợi nhuận ròng tăng thêm/năm -- 680.000.000 VNĐ
Thời gian hoàn vốn thiết bị -- 6.2 tháng
-------------------------------------------------------------------------
Lộ trình triển khai thực tế
- Giai đoạn 1 (Tháng 1 - 2): Phân loại gỗ tròn tại bãi tập kết dựa theo đường kính và lập sơ đồ xẻ phôi tối ưu hóa đường cắt tiếp tuyến/xuyên tâm.
- Giai đoạn 2 (Tháng 3 - 4): Nạp phần mềm điều khiển lò sấy tự động với bảng tra cứu gradient ẩm Sa mộc theo chuẩn thực nghiệm.
- Giai đoạn 3 (Tháng 5 trở đi): Chuyển giao công nghệ ghép nối phôi định hướng bán kính cho các hợp tác xã chế biến đồ gỗ thủ công mỹ nghệ tại Lào Cai.
Hạn chế và hướng phát triển
Hạn chế kỹ thuật
- Quy mô mẫu thử: Nghiên cứu tập trung trên 3 cây mẫu thuộc một cấp tuổi (19 năm tuổi) tại xã San Sả Hồ, Sa Pa; chưa bao quát các cấp tuổi khai thác trẻ hơn (10 - 15 năm) hoặc các tiểu vùng sinh thái khác (Hà Giang, Sơn La).
- Phạm vi chỉ tiêu: Chưa kết hợp khảo sát các chỉ tiêu cơ học phá hủy (độ bền uốn tĩnh MOR, mô đun đàn hồi MOE, ứng suất nén dọc thớ) song song trên cùng một phôi mẫu.
- Đặc tính hóa học: Chưa đo đạc định lượng hàm lượng tinh dầu và hàm lượng tro/lignin ảnh hưởng đến tốc độ nhả ẩm.
Hướng nghiên cứu tiếp theo
- Mở rộng nghiên cứu ứng dụng phương pháp quang phổ cận hồng ngoại (NIRS - Near Infrared Spectroscopy) để dự báo nhanh không phá hủy độ ổn định kích thước gỗ Sa mộc ngay tại hiện trường.
- Thử nghiệm công nghệ biến tính gỗ (Wood Modification) bằng phương pháp xử lý nhiệt dầu (Thermo-oil) hoặc ngâm tẩm hóa chất thân thiện môi trường (PEG, Axetyl hóa) nhằm hạ tỷ số co rút xuống dưới $2.0%$.
- Mô hình hóa ứng suất nhiệt ẩm bằng phương pháp phần tử hữu hạn (FEM) mô phỏng chính xác biến dạng của tấm ván Sa mộc trong quá trình sấy công nghiệp.
Đối tượng hưởng lợi
- Sinh viên và Học viên Cao học Lâm nghiệp: Tiếp cận phương pháp luận chuẩn mực về giải phẫu thực nghiệm và phương pháp kiểm định cơ lý gỗ theo tiêu chuẩn quốc gia TCVN 8048.
- Kỹ sư Chế biến Lâm sản & Vận hành Lò sấy: Nắm bắt bộ thông số chính xác ($KLTT_0 = 0.33\text{ g/cm}^3$, $CRTT = 5.14%$, $CRXT = 2.72%$) để lập trình lò sấy chính xác, giảm tỷ lệ hư hỏng ván xẻ từ $23%$ xuống dưới $7%$.
- Doanh nghiệp Sản xuất Đồ gỗ Nội ngoại thất: Nâng cao chất lượng sản phẩm xuất khẩu từ gỗ Sa mộc, tối ưu hóa định mức tiêu hao nguyên liệu phôi xẻ thêm $15 - 20%$.
- Nhà khoa học & Cơ quan Quản lý Lâm nghiệp: Bổ sung cơ sở dữ liệu khoa học phục vụ công tác quy hoạch phát triển rừng trồng gỗ lớn bền vững và cấp chứng chỉ rừng FSC tại vùng Tây Bắc.
Câu hỏi thường gặp
1. Yêu cầu kỹ thuật tối thiểu để ứng dụng kết quả nghiên cứu vào lò sấy công nghiệp là gì?
Cần trang bị lò sấy hơi nước hoặc sấy nhiệt sạch có hệ thống cảm biến đo ẩm tự động (cảm biến điện trở tiếp xúc trực tiếp lõi gỗ) với dải điều khiển nhiệt độ $40 - 80^\circ\text{C}$, bước điều chỉnh $\pm 1^\circ\text{C}$, và quạt đối lưu đảo chiều gió đạt vận tốc gió qua phôi $1.5 - 2.5\text{ m/s}$.
2. Vì sao độ co rút tiếp tuyến ($CRTT$) lại lớn gấp gần 2 lần độ co rút xuyên tâm ($CRXT$)?
Hiện tượng này do cấu trúc giải phẫu gỗ quyết định: Các tia gỗ chạy theo hướng xuyên tâm đóng vai trò như các "chốt cơ học" cản trở sự co rút theo chiều xuyên tâm. Ngược lại, theo hướng tiếp tuyến không có cấu trúc cản trở, đồng thời các vách tế bào có độ dày lớn co ngắn mạnh khi các phân tử nước liên kết thoát khỏi các mixen cellulose.
3. Có thể dùng khối lượng thể tích ($KLTT_0$) để suy đoán trực tiếp độ co rút của gỗ Sa mộc không?
Không. Kết quả phân tích hồi quy trong nghiên cứu chỉ ra rằng hệ số tương quan giữa $KLTT_0$ và $CRXT, CRTT$ có $p > 0.05$ (không có ý nghĩa thống kê). Mức độ co rút phụ thuộc nhiều hơn vào góc nghiêng sợi cellulose và cấu trúc thành tế bào hơn là mật độ thể tích đơn thuần.
4. Gỗ Sa mộc sấy xong đạt độ ẩm bao nhiêu thì thích hợp để gia công nội thất tại miền Bắc?
Độ ẩm cân bằng ($EMC$) khuyến nghị cho đồ gỗ Sa mộc tại khu vực miền Bắc dao động từ $10%$ đến $12%$. Gia công tại dải ẩm này sẽ triệt tiêu nguy cơ co ngót vào mùa khô hanh hoặc phồng rộp vào mùa nồm ẩm.
5. Chi phí đầu tư thử nghiệm phân cỡ phôi gỗ Sa mộc theo bán kính có tốn kém không?
Chi phí tăng thêm không đáng kể (chủ yếu là chi phí sắp xếp và đánh dấu phôi sau xẻ phách, chiếm $< 2%$ chi phí nhân công), nhưng giúp tăng giá trị gia tăng của lô ván thành phẩm lên từ $15%$ đến $25%$ nhờ giảm thiểu khuyết tật loại bỏ.
Kết luận
Công trình nghiên cứu đã giải quyết trọn vẹn bài toán thực nghiệm về độ ổn định kích thước của gỗ Sa mộc (Cunninghamia lanceolata) 19 năm tuổi trồng tại Sa Pa, Lào Cai. Kết quả khẳng định khối lượng thể tích khô kiệt của loài đạt trung bình $0.33\text{ g/cm}^3$ và phân bố đồng đều từ tâm ra vỏ. Độ co rút và giãn nở tiếp tuyến ($5.14%$ và $5.03%$) có xu hướng tăng dần theo hướng ra phía vỏ, trong khi các chỉ tiêu xuyên tâm ($2.72%$ và $3.00%$) duy trì sự ổn định cao.
Những kết luận khoa học này là chìa khóa then chốt để các nhà nghiên cứu, kỹ sư và doanh nghiệp ngành gỗ thiết kế quy trình sấy khoa học, giảm thiểu tối đa tổn thất kinh tế do biến dạng phôi gỗ. Đề tài mở ra tiềm năng thương mại to lớn cho gỗ rừng trồng bản địa Sa mộc, góp phần nâng cao chuỗi giá trị lâm sản Việt Nam trên thị trường quốc tế.
