Tổng quan nghiên cứu

Ngành công nghiệp sản xuất bột giấy và giấy giữ vai trò huyết mạch trong nền kinh tế quốc dân, trong đó Tổng công ty Giấy Việt Nam và Nhà máy Giấy Bãi Bằng là đơn vị sản xuất chủ lực với bề dày hoạt động lâu đời. Năm 2003, doanh nghiệp hoàn thành dự án mở rộng giai đoạn 1, nâng công suất bột giấy lên 48.000 tấn/năm và giấy in, viết đạt 50.000 tấn/năm, tiếp tục triển khai các dự án quy mô lớn như dây chuyền bột giấy tẩy trắng 250.000 tấn/năm và mở rộng công suất sản xuất giấy đạt mức 120.000 tấn/năm vào năm 2010. Tuy nhiên, bối cảnh kinh tế mở cửa tạo ra làn sóng đầu tư mạnh mẽ từ các tập đoàn quốc tế của Nhật Bản, Trung Quốc, Thái Lan cùng sự gia tăng cạnh tranh khốc liệt từ nguồn giấy nhập khẩu và các loại hình truyền thông kỹ thuật số. Thực trạng này đặt doanh nghiệp trước yêu cầu sống còn phải tối ưu hóa chi phí và chủ động hoàn toàn nguồn nguyên liệu đầu vào.

Vấn đề then chốt đặt ra cho Nhà máy Giấy Bãi Bằng là nguồn gỗ nguyên liệu phục vụ sản xuất chiếm hơn 60% cơ cấu giá thành nhưng lại chịu nhiều biến động phức tạp về chu kỳ sinh trưởng, thời vụ khai thác và sự dao động của thị trường thu mua tự do. Mục tiêu nghiên cứu của luận văn là khảo sát các phương pháp khai phá dữ liệu hiện đại, phân tích đặc tính chuỗi thời gian thực tế và xây dựng mô hình dự báo định lượng chính xác nhu cầu nguồn gỗ nguyên liệu phục vụ kế hoạch sản xuất ổn định lâu dài. Phạm vi nghiên cứu tập trung vào chuỗi dữ liệu vận hành thực tế tại Nhà máy Giấy Bãi Bằng trong giai đoạn 10 năm từ năm 2003 đến năm 2012, hướng đến xây dựng các phương án dự báo ngắn hạn và trung hạn từ 3 tháng đến 3 năm. Nghiên cứu mang ý nghĩa thực tiễn to lớn khi giúp nâng cao độ chính xác dự báo nguồn cung lên trên 85%, giảm thiểu 15% đến 20% rủi ro đứt gãy nguyên liệu hoặc tồn trữ dư thừa gây ứ đọng vốn lưu động.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng lý thuyết phát hiện tri thức từ cơ sở dữ liệu (Knowledge Discovery in Databases - KDD) và phương pháp luận phân tích chuỗi thời gian kinh tế lượng Box-Jenkins. Khái niệm cốt lõi đầu tiên là khai phá dữ liệu, được hiểu như một quy trình kỹ thuật tự động nhằm trích xuất các thông tin ẩn, hữu ích và có tính dự báo từ các tập dữ liệu có quy mô lớn.

Khái niệm thứ hai là chuỗi thời gian thực (Time Series Data), biểu diễn tập hợp các giá trị của biến quan sát theo trình tự thời gian đều đặn với bốn thành phần cấu thành: thành phần xu thế dài hạn (T), thành phần mùa vụ tuần hoàn (S), thành phần chu kỳ (C) và thành phần sai số ngẫu nhiên bất thường (I).

Khái niệm thứ ba là mô hình tự hồi quy tích hợp trung bình trượt ARIMA(p,d,q). Mô hình này hợp nhất ba quá trình toán học: quá trình tự hồi quy AR(p) mô tả mối quan hệ tuyến tính giữa giá trị hiện tại với các giá trị trễ quá khứ; bậc tích hợp sai phân d(d) biến đổi chuỗi không dừng về trạng thái dừng có kỳ vọng và phương sai không đổi; và quá trình trung bình trượt MA(q) biểu diễn tác động của các cú sốc sai số ngẫu nhiên hiện tại cùng các giá trị trễ của nhiễu trắng.

Phương pháp nghiên cứu

Nguồn dữ liệu của luận văn là tập dữ liệu thứ cấp được trích xuất từ hệ thống thống kê kế hoạch, nhật ký sản xuất và phòng vật tư của Nhà máy Giấy Bãi Bằng, trải dài liên tục qua 120 tháng quan sát từ tháng 1 năm 2003 đến tháng 12 năm 2012.

Phương pháp chọn mẫu là chọn mẫu toàn bộ theo chuỗi thời gian định kỳ hàng tháng với 12 kỳ quan sát mỗi năm, tập trung vào hai chuỗi biến số định lượng then chốt: sản lượng gỗ tự cấp từ các vùng nguyên liệu lâm nghiệp liên kết và sản lượng gỗ thu mua ngoài thị trường tự do.

Lý do lựa chọn phương pháp phân tích chuỗi thời gian Box-Jenkins kết hợp phần mềm chuyên dụng EViews là nhờ khả năng xử lý chuyên sâu các biến động chuỗi kinh tế. Phương pháp này vượt trội hơn các phương pháp phán đoán chuyên gia Delphi hay mạng nơ-ron đa lớp vốn dễ rơi vào tình trạng quá khớp (overfitting) với số lượng tham số lớn. Thông qua quy trình 4 bước chuẩn tắc gồm nhận dạng mô hình qua đồ thị tự tương quan ACF và tự tương quan riêng phần PACF, ước lượng tham số bằng phương pháp bình phương nhỏ nhất (OLS), kiểm định chẩn đoán tính nhiễu trắng của phần dư và tiến hành dự báo điểm với khoảng tin cậy xác thực.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Phát hiện thứ nhất chỉ ra cơ cấu nguồn gỗ nguyên liệu của Nhà máy Giấy Bãi Bằng tồn tại sự mất cân đối lớn giữa nguồn lực tự chủ và thị trường ngoài. Lượng gỗ tự cấp chỉ đáp ứng được khoảng 35% đến 40% tổng nhu cầu sản xuất thực tế, trong khi nhà máy phải phụ thuộc từ 60% đến 65% vào nguồn gỗ thu mua bên ngoài từ các đơn vị và hộ dân trồng rừng.

Phát hiện thứ hai làm sáng tỏ tính biến động mùa vụ rõ nét của nguồn nguyên liệu đầu vào. Sản lượng tiêu thụ và nhập kho gỗ nguyên liệu tăng mạnh khoảng 25% trong giai đoạn quý 2 từ tháng 3 đến tháng 6 hàng năm nhằm đẩy mạnh sản xuất giấy in, giấy viết phục vụ nhu cầu mùa tựu trường, sau đó giảm sâu khoảng 18% trong các tháng từ tháng 8 đến tháng 11.

Phát hiện thứ ba khẳng định mô hình ARIMA(1,0,1) là mô hình kinh tế lượng tối ưu nhất để mô phỏng tương quan giữa sản lượng gỗ tự cấp và gỗ mua ngoài. Mô hình đạt hệ số xác định cao, mức sai số dự báo trung bình được kiểm soát dưới 5%, đồng thời phản ánh mối quan hệ hồi quy tuyến tính chặt chẽ giữa lưu lượng bột giấy sản xuất và công suất 120.000 tấn giấy thành phẩm mỗi năm.

Thảo luận kết quả

Sự biến động mạnh của nguồn gỗ mua ngoài bắt nguồn trực tiếp từ việc thiếu quy hoạch liên kết vùng nguyên liệu bền vững và sự cạnh tranh thu mua gay gắt từ các xưởng sản xuất dăm gỗ xuất khẩu tư nhân. So sánh với các nghiên cứu kinh tế lượng của Armstrong và Makridakis, kết quả này hoàn toàn thống nhất với nhận định rằng khi dự báo các chuỗi dữ liệu sản xuất công nghiệp, việc phân rã yếu tố mùa vụ kết hợp xử lý sai phân chuỗi dừng mang lại độ tin cậy cao hơn hẳn các mô hình ngoại suy xu hướng đơn thuần.

Dữ liệu nghiên cứu được trình bày trực quan hóa thông qua hệ thống đồ thị và bảng biểu chuyên sâu trên EViews. Cụ thể, biểu đồ tương quan phân tán (scatter plot) kết hợp đường xu hướng tuyến tính biểu thị rõ nét mối quan hệ tỷ lệ thuận giữa lượng gỗ tự cấp và mua ngoài. Đồ thị hàm tự tương quan mẫu SAC và tự tương quan riêng phần SPAC minh họa trực quan quá trình tắt lịm dần của chuỗi dừng, giúp xác định chuẩn xác bậc tự hồi quy p = 1 và bậc trung bình trượt q = 1. Đồng thời, biểu đồ đường so sánh giữa số liệu thực tế và số liệu dự báo qua 120 chu kỳ quan sát cho thấy đường dự báo bám sát chặt chẽ đường thực nghiệm với độ lệch dư phân tán ngẫu nhiên dạng nhiễu trắng thuần túy.

Đề xuất và khuyến nghị

  • Chuẩn hóa quy trình số hóa và làm sạch dữ liệu đầu vào: Tiến hành xây dựng hệ thống phần mềm quản lý cơ sở dữ liệu chuỗi thời gian tập trung tại Nhà máy Giấy Bãi Bằng, phấn đấu đạt tỷ lệ tự động hóa 100% việc cập nhật số liệu nhập gỗ hàng tuần trong vòng 6 tháng tới, do Phòng Công nghệ Thông tin phối hợp cùng Phòng Kế hoạch Vật tư trực tiếp triển khai.
  • Ứng dụng mô hình ARIMA(1,0,1) vào lập kế hoạch tác nghiệp: Ban Giám đốc nhà máy cần tích hợp mô hình ARIMA(1,0,1) trên phần mềm EViews vào quy trình lập kế hoạch điều hành sản xuất định kỳ 3 tháng và 12 tháng, hướng tới mục tiêu duy trì sai số dự báo nguồn gỗ dưới ngưỡng 6% trong suốt lộ trình 12 tháng áp dụng.
  • Tái cấu trúc và mở rộng diện tích vùng nguyên liệu tự cấp: Nâng tỷ trọng cung ứng gỗ tự cấp từ mức 35% hiện tại lên tối thiểu 50% tổng nhu cầu sản xuất trong giai đoạn 3 đến 5 năm tới, do Phòng Nguyên liệu chủ trì kết hợp cùng các công ty lâm nghiệp địa phương triển khai mở rộng diện tích trồng rừng gỗ lớn.
  • Thiết lập chuỗi liên kết bao tiêu giá ổn định với nông dân: Ký kết hợp đồng bao tiêu dài hạn với trên 500 hộ gia đình và hợp tác xã trồng rừng trong bán kính 100 km, nhằm kiểm soát biên độ biến động giá mua ngoài không vượt quá 10%/năm, hoàn thành triển khai trong vòng 2 năm do Phòng Vật tư Thương mại phụ trách.

Đối tượng nên tham khảo luận văn

  • Ban giám đốc và nhà quản trị chuỗi cung ứng ngành giấy: Nhận được giải pháp định lượng khoa học để cân đối nguồn nguyên liệu đầu vào, giảm áp lực biến động giá thị trường, tối ưu hóa quy mô vận hành từ 50.000 tấn đến 120.000 tấn giấy/năm và hạ giá thành sản phẩm.
  • Chuyên viên phân tích dữ liệu và chuyên gia kinh tế lượng: Nắm vững quy trình ứng dụng thực tế của mô hình Box-Jenkins ARIMA trên phần mềm EViews đối với các bài toán chuỗi thời gian công nghiệp thực tế có tính mùa vụ cao.
  • Giảng viên, học viên cao học và sinh viên ngành Khoa học máy tính: Sử dụng làm tài liệu tham khảo chuẩn mực về phương pháp luận nghiên cứu phát hiện tri thức từ cơ sở dữ liệu (KDD), kỹ thuật xử lý dữ liệu chuỗi thời gian và cấu trúc trình bày luận văn thạc sĩ.
  • Các doanh nghiệp sản xuất chế biến lâm sản và bột giấy tại Việt Nam: Tham khảo kinh nghiệm thực tế về mô hình hóa nguồn cung, xác định điểm cân bằng giữa nguồn lực nội bộ và thu mua thương mại để xây dựng chiến lược phát triển bền vững.

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

Vì sao luận văn lựa chọn mô hình ARIMA thay vì áp dụng mạng Nơ-ron nhân tạo?
Mạng Nơ-ron thường đòi hỏi tập dữ liệu rất lớn và tiềm ẩn rủi ro quá khớp nghiêm trọng khi cấu trúc mạng phức tạp, đồng thời khó giải thích ý nghĩa kinh tế. Với chuỗi dữ liệu 120 tháng quan sát, mô hình ARIMA(1,0,1) mang lại tính tường minh cao, cấu trúc gọn nhẹ và kiểm soát sai số ổn định dưới 5%.

Quy trình phát hiện tri thức KDD trong luận văn bao gồm những bước cơ bản nào?
Quy trình KDD gồm 6 bước liên hoàn: thu thập và lựa chọn dữ liệu gỗ nguyên liệu, làm sạch và tiền xử lý dữ liệu thô, chuyển đổi chuỗi thời gian về dạng dừng, trích xuất mẫu dự báo qua mô hình ARIMA, kiểm định chẩn đoán phần dư và ứng dụng mô hình vào lập kế hoạch thực tế tại nhà máy.

Cơ cấu nguồn gỗ tự cấp và mua ngoài tác động như thế nào đến rủi ro sản xuất?
Nguồn gỗ tự cấp chỉ đáp ứng khoảng 35% đến 40% khiến nhà máy phải phụ thuộc hơn 60% vào thị trường tự do. Sự mất cân đối này làm tăng độ nhạy cảm trước các đợt tăng giá nguyên liệu từ 8% đến 12% trong mùa cao điểm, đe dọa trực tiếp đến tính liên tục của dây chuyền sản xuất.

Yếu tố tính mùa vụ của chuỗi thời gian nguyên liệu gỗ thể hiện rõ nhất vào giai đoạn nào?
Yếu tố mùa vụ xuất hiện rõ rệt nhất vào quý 2 từ tháng 3 đến tháng 6 khi nhu cầu sản xuất giấy in, giấy viết phục vụ niên khóa học mới tăng trên 25%. Ngược lại, nhu cầu tiêu thụ nguyên liệu giảm bình quân 18% từ tháng 8 đến tháng 11 hàng năm.

Phần mềm EViews hỗ trợ những tính năng trọng yếu nào trong việc xây dựng mô hình?
EViews đóng vai trò công cụ phân tích kinh tế lượng giúp nhận dạng tính dừng qua đồ thị tự tương quan ACF và PACF, thực hiện sai phân bậc d, ước lượng các tham số mô hình bằng thuật toán OLS, kiểm tra tính phân phối chuẩn của phần dư và kết xuất dự báo điểm với độ tin cậy 95%.

Kết luận

  • Hệ thống hóa toàn diện nền tảng lý thuyết khai phá dữ liệu chuỗi thời gian và quy trình phát hiện tri thức KDD ứng dụng trong sản xuất công nghiệp.
  • Khảo sát và phân tích sâu sắc thực trạng biến động nguồn gỗ nguyên liệu tại Nhà máy Giấy Bãi Bằng qua chuỗi quan sát 120 tháng liên tục.
  • Xây dựng và kiểm định thành công mô hình dự báo ARIMA(1,0,1) trên phần mềm EViews với độ chính xác cao và sai số dưới 5%.
  • Làm rõ cấu trúc phụ thuộc 65% vào nguồn gỗ mua ngoài cùng quy luật mùa vụ tăng 25% trong quý 2 hàng năm của nhà máy.
  • Đề xuất hệ thống giải pháp quản lý đồng bộ từ nâng cấp phần mềm dữ liệu đến mở rộng vùng rừng tự cấp lên 50% tổng nhu cầu.

Đóng góp quan trọng nhất của luận văn là việc kết nối thành công phương pháp luận khoa học máy tính với bài toán thực tiễn quản trị chuỗi cung ứng lâm sản. Kế hoạch tiếp theo cần tập trung tích hợp mô hình dự báo vào hệ thống ERP của doanh nghiệp trong 6 tháng tới và mở rộng chu kỳ dự báo lên 36 tháng. Các nhà quản lý và nhà nghiên cứu quan tâm có thể tiếp tục phát triển hướng đề tài này bằng cách tích hợp thêm các yếu tố ngoại sinh như biến động giá xăng dầu và chính sách xuất khẩu dăm gỗ để gia tăng độ chuẩn xác cho mô hình.