Tổng quan nghiên cứu

Nhu cầu nguyên liệu gỗ phục vụ công nghiệp sản xuất giấy tại Việt Nam luôn đối mặt với nhiều biến động phức tạp về thị trường và nguồn cung. Năm 2010, Tổng công ty Giấy Việt Nam đã tiến hành đầu tư nâng công suất sản xuất giấy lên mức 120.000 tấn/năm, đồng thời xây dựng quy hoạch dự án bột giấy tẩy trắng với công suất thiết kế đạt 250.000 tấn/năm. Trước đó, dự án mở rộng giai đoạn I vào năm 2003 đã nâng sản lượng bột giấy từ 48.000 tấn/năm lên trên 50.000 tấn/năm. Tuy nhiên, nhiều dự án mở rộng quy mô ngành giấy thời gian qua bị đình trệ hoặc suy giảm hiệu quả sinh lời do thiếu hụt công tác phân tích, dự báo chính xác nguồn cung gỗ đầu vào.

Đứng trước nguy cơ cạnh tranh gay gắt từ các tập đoàn nước ngoài và sự thay đổi trong thói quen tiêu dùng sản phẩm in ấn, bài toán tự chủ và tối ưu hóa nguồn nguyên liệu trở thành điều kiện sống còn của doanh nghiệp. Nghiên cứu tập trung giải quyết vấn đề khai phá dữ liệu chuỗi thời gian nhằm xây dựng mô hình dự báo chính xác nhu cầu và khả năng đáp ứng nguồn gỗ nguyên liệu cho Nhà máy Giấy Bãi Bằng. Mục tiêu cụ thể là lượng hóa xu hướng biến động của nguồn gỗ tự cung cấp và gỗ mua ngoài, từ đó đưa ra dự báo trung hạn cho giai đoạn 10 năm tiếp theo đến năm 2020.

Phạm vi nghiên cứu được triển khai dựa trên tập dữ liệu lịch sử giai đoạn 2000 - 2012 tại Nhà máy Giấy Bãi Bằng và các vùng nguyên liệu phụ cận khu vực miền Bắc, Bắc Trung Bộ. Ý nghĩa của nghiên cứu thể hiện qua việc cung cấp công cụ định lượng hỗ trợ nhà quản lý tối ưu hóa tồn kho, giảm thiểu rủi ro biến động giá mua ngoài khoảng 15% đến 20%, tạo lợi thế cạnh tranh bền vững cho thương hiệu giấy Bãi Bằng trên thị trường nội địa.

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 quy trình phát hiện tri thức từ cơ sở dữ liệu và lý thuyết phân tích chuỗi thời gian thực. Chuỗi thời gian được hiểu là một tập hợp các giá trị quan sát của biến thực theo các khoảng thời gian bằng nhau, bao gồm bốn thành phần cơ bản: thành phần xu thế dài hạn, thành phần mùa vụ, thành phần chu kỳ và thành phần bất thường ngẫu nhiên.

Mô hình toán học trung tâm là phương pháp Box-Jenkins với mô hình tích hợp tự hồi quy trung bình trượt ARIMA. Mô hình này kết hợp quá trình tự hồi quy bậc p (biểu diễn giá trị hiện tại theo các giá trị trễ quá khứ) và quá trình trung bình trượt bậc q (biểu diễn giá trị theo trọng số của các sai số ngẫu nhiên hiện hành và trễ). Để vận hành ARIMA, chuỗi số liệu phải thỏa mãn tính dừng yếu, nghĩa là giá trị trung bình và phương sai không đổi theo thời gian, đồng phương sai chỉ phụ thuộc vào độ trễ thời gian. Các khái niệm bổ trợ then chốt bao gồm hàm tự tương quan mẫu, hàm tự tương quan từng phần và kiểm định nhiễu trắng thuần túy nhằm đánh giá phần dư của mô hình.

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

Nghiên cứu sử dụng tập dữ liệu chuỗi thời gian thực thường niên kéo dài 13 năm liên tục từ năm 2000 đến năm 2012. Cỡ mẫu nghiên cứu gồm 13 quan sát hàng năm đối với từng biến số chính: sản lượng bột giấy, sản lượng giấy thành phẩm, lượng gỗ nguyên liệu tự cung cấp, lượng gỗ mua ngoài và lượng gỗ xuất khẩu.

Phương pháp chọn mẫu là thu thập dữ liệu thứ cấp toàn phần từ hệ thống thống kê chính thức của Tổng công ty Giấy Việt Nam và Nhà máy Giấy Bãi Bằng. Lý do lựa chọn mô hình ARIMA và phần mềm kinh tế lượng EViews xuất phát từ đặc thù của chuỗi số liệu kinh tế có dung lượng mẫu vừa phải, không đòi hỏi các biến ngoại sinh phức tạp nhưng vẫn đảm bảo độ tin cậy cao nhờ nguyên lý để dữ liệu tự nói. Quy trình nghiên cứu được triển khai theo 4 bước lặp chuẩn hóa của Box-Jenkins: nhận dạng mô hình thông qua biểu đồ tự tương quan; ước lượng các tham số bằng phương pháp bình phương nhỏ nhất; kiểm định chẩn đoán tính nhiễu trắng của phần dư; và thực hiện dự báo xu hướng đến năm 2020. Toàn bộ quá trình phân tích và ước lượng được hoàn thành trong quý IV năm 2013.

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

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

Quá trình khai phá dữ liệu và ước lượng mô hình ARIMA(1,0,1) trên phần mềm EViews đã đem lại bốn phát hiện trọng tâm có giá trị thực tiễn cao:

Thứ nhất, lượng gỗ nguyên liệu tự cung cấp của đơn vị duy trì xu hướng tăng trưởng tích cực từ 120.350 nghìn tấn vào năm 2000 lên mức đỉnh 177.099 nghìn tấn năm 2008, đạt mức tăng 47,15%. Sau giai đoạn điều chỉnh, nguồn tự cấp ổn định ở mức 145.950 nghìn tấn vào năm 2012, thể hiện khả năng duy trì vùng nguyên liệu tự chủ ở quy mô trên 140.000 nghìn tấn/năm.

Thứ hai, sản lượng gỗ mua ngoài phát sinh đột biến sau dự án mở rộng giai đoạn I. Từ mức 0 nghìn tấn trong giai đoạn 2000 - 2002, lượng mua ngoài đã tăng vọt lên 98.000 nghìn tấn năm 2003 và đạt đỉnh 166.490 nghìn tấn vào năm 2011. Tỷ trọng gỗ mua ngoài chiếm từ 40% đến 57% tổng nhu cầu nguyên liệu của nhà máy trong giai đoạn 2005 - 2012.

Thứ ba, sản lượng bột giấy và giấy thành phẩm có mối tương quan thuận rất chặt chẽ. Sản lượng bột giấy tăng từ 48.250 nghìn tấn năm 2000 lên 75.183 nghìn tấn năm 2008, tương ứng mức tăng 55,82%. Sản lượng giấy thành phẩm tăng tương ứng từ 81.356 nghìn tấn lên 112.031 nghìn tấn, phản ánh hiệu suất vận hành của các dây chuyền mở rộng.

Thứ tư, nguồn cung gỗ tại khu vực phía Bắc ghi nhận sự bùng nổ của sản lượng xuất khẩu dăm gỗ ra thị trường quốc tế. Từ mức 150.565 nghìn tấn năm 2004, lượng xuất khẩu đã tăng phi mã lên 2.350.200 nghìn tấn vào năm 2011, gấp hơn 15,6 lần sau 7 năm.

Thảo luận kết quả

Trên các biểu đồ chuỗi thời gian và đồ thị hồi quy tương quan tuyến tính kết xuất từ EViews, đường xu thế dài hạn của biến tự cấp và mua ngoài phản ánh rõ ràng tính bù trừ lẫn nhau trong cung ứng sản xuất. Đồ thị hàm tự tương quan mẫu cho thấy chuỗi sai phân bậc một đạt tính dừng nhanh chóng, giúp mô hình ARIMA(1,0,1) đạt hệ số ước lượng tin cậy cao, phần dư hoàn toàn là nhiễu trắng ngẫu nhiên.

Nguyên nhân cốt lõi dẫn đến tình trạng xuất khẩu gỗ dăm tăng vọt là do diện tích rừng trồng nguyên liệu tại các tỉnh miền Bắc và Bắc Trung Bộ phát triển mạnh mẽ vượt xa năng lực chế biến tại chỗ. Kết quả này hoàn toàn tương đồng với các báo cáo lâm nghiệp ngành giấy khi chỉ ra tình trạng dư thừa nguồn cung thô cục bộ nhưng thiếu hụt nguồn cung chế biến sâu. So với các nghiên cứu kinh tế lượng đa biến, mô hình chuỗi thời gian ARIMA tuy chỉ sử dụng độ trễ nội tại nhưng lại phản ánh chính xác xu thế quán tính, dự báo đến năm 2020 nguồn nguyên liệu nội địa hoàn toàn đủ khả năng cung ứng cho việc tăng công suất nếu nhà máy có chính sách thu mua hợp lý.

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

Dựa trên kết quả dự báo và phân tích cung cầu thực tế, bốn nhóm giải pháp chiến lược được đề xuất như sau:

Thứ nhất, tái cấu trúc quy hoạch vùng nguyên liệu và chính sách bao tiêu nhằm nâng cao tỷ lệ tự chủ nguồn cung gỗ lên mức 65% đến 70% tổng sản lượng đầu vào. Ban Lãnh đạo Tổng công ty Giấy Việt Nam cần trực tiếp chỉ đạo triển khai kế hoạch liên kết dài hạn với các công ty lâm nghiệp địa phương trong giai đoạn 2014 - 2018, bảo đảm ổn định tối thiểu 160.000 nghìn tấn gỗ/năm.

Thứ hai, đẩy mạnh đầu tư mở rộng dây chuyền sản xuất bột giấy thương phẩm và chế biến gỗ sâu để tận dụng nguồn nguyên liệu dôi dư đang bị xuất khẩu thô. Khối Kỹ thuật và Ban Quản lý dự án Nhà máy Giấy Bãi Bằng cần hoàn thiện báo cáo khả thi giai đoạn 2015 - 2020, hướng tới mục tiêu gia tăng giá trị gia tăng ngành giấy thêm khoảng 25% đến 30%.

Thứ ba, thiết lập hệ thống cảnh báo và thu mua linh hoạt nhằm kiểm soát biến động giá nguyên liệu mua ngoài. Phòng Thu mua nguyên liệu phối hợp cùng Phòng Kế hoạch thực hiện theo dõi sát sao thị trường xuất khẩu dăm gỗ trong giai đoạn 2014 - 2016, duy trì lượng mua ngoài từ 90.000 đến 110.000 nghìn tấn/năm với mức giá cạnh tranh.

Thứ tư, số hóa công tác dự báo kinh tế lượng thông qua việc chuẩn hóa quy trình chạy mô hình ARIMA trên phần mềm EViews tại cấp quản lý vận hành. Phòng Công nghệ thông tin có trách nhiệm định kỳ cập nhật chuỗi dữ liệu hàng quý từ năm 2014, phấn đấu duy trì sai số dự báo trung bình dưới 5%.

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

Luận văn cung cấp nhiều giá trị khoa học và ứng dụng thực tiễn cho bốn nhóm đối tượng chuyên môn:

Nhóm cán bộ quản trị chiến lược và ban giám đốc tại các nhà máy giấy, bột giấy, doanh nghiệp lâm sản có thể ứng dụng kết quả dự báo nguồn gỗ để lập kế hoạch mở rộng nhà xưởng, quyết định đầu tư máy móc và chuẩn bị ngân sách thu mua nguyên vật liệu.

Nhóm chuyên viên phân tích chuỗi cung ứng và kinh tế lượng có thể tham khảo quy trình thực hành chi tiết 4 bước của phương pháp Box-Jenkins, từ xử lý tính dừng chuỗi số liệu đến ước lượng mô hình ARIMA(1,0,1) trên phần mềm EViews.

Nhóm giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Khoa học máy tính, Hệ thống thông tin quản lý có thể sử dụng đề tài làm bài học thực nghiệm mẫu mực về ứng dụng khai phá dữ liệu vào bài toán kinh tế công nghiệp.

Nhóm cán bộ cơ quan quản lý nhà nước về tài nguyên rừng và các hiệp hội ngành nghề có thêm căn cứ dữ liệu để hoạch định chính sách phát triển rừng trồng nguyên liệu, cân đối giữa chế biến trong nước và hạn chế xuất khẩu thô.

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

Mô hình ARIMA(1,0,1) áp dụng trong luận văn có ưu điểm gì nổi bật?
Mô hình ARIMA(1,0,1) kết hợp tối ưu giữa yếu tố tự hồi quy bậc 1 và trung bình trượt bậc 1, giúp nắm bắt trọn vẹn quán tính lịch sử và triệt tiêu sai số ngẫu nhiên của chuỗi số liệu nguyên liệu gỗ giai đoạn 2000 - 2012 mà không cần thu thập quá nhiều biến số kinh tế vĩ mô phức tạp.

Tại sao tác giả lại lựa chọn phần mềm EViews để xử lý chuỗi thời gian?
EViews là phần mềm chuyên dụng hàng đầu trong phân tích chuỗi thời gian kinh tế lượng, cung cấp đầy đủ công cụ vẽ đồ thị kiểm tra tính dừng, đồ thị tương quan mẫu và thuật toán ước lượng bình phương nhỏ nhất giúp thao tác nhanh chóng cùng kết quả trực quan chính xác.

Dữ liệu quan sát 13 năm có đủ độ tin cậy để đưa ra dự báo trung hạn đến năm 2020 không?
Chuỗi dữ liệu 13 năm (2000 - 2012) phản ánh đầy đủ chu kỳ biến động trước và sau khi nhà máy mở rộng công suất. Dù cỡ mẫu vừa phải, quy trình kiểm định phần dư nhiễu trắng cho thấy mô hình đạt độ ổn định vững chắc để dự báo xu thế trung hạn 5 đến 10 năm.

Lượng gỗ xuất khẩu bùng nổ ảnh hưởng thế nào đến bài toán nguyên liệu của Bãi Bằng?
Sự gia tăng xuất khẩu từ 150.565 nghìn tấn năm 2004 lên 2.350.200 nghìn tấn năm 2011 chứng minh nguồn cung gỗ vùng phụ cận rất dồi dào, đồng thời tạo sức ép cạnh tranh thu mua, đòi hỏi Bãi Bằng phải nâng cao năng lực chế biến để giữ chân nguồn gỗ nội địa.

Hạn chế chính của mô hình dự báo chuỗi thời gian trong luận văn là gì?
Hạn chế lớn nhất là mô hình ARIMA đơn biến chỉ dựa vào quy luật quá khứ nội tại của dữ liệu, chưa tích hợp được tác động đột ngột từ biến động chính sách xuất nhập khẩu, thuế quan hoặc các cú sốc kinh tế vĩ mô ngoài dự kiến.

Kết luận

  • Nghiên cứu đã hệ thống hóa toàn diện cơ sở lý thuyết về khai phá dữ liệu chuỗi thời gian và phương pháp luận Box-Jenkins ứng dụng trong công nghiệp sản xuất giấy.
  • Xây dựng thành công mô hình dự báo ARIMA(1,0,1) trên phần mềm EViews dựa trên chuỗi số liệu thực tế giai đoạn 2000 - 2012 của Nhà máy Giấy Bãi Bằng.
  • Khẳng định nguồn cung gỗ nguyên liệu tại khu vực miền Bắc và Bắc Trung Bộ hoàn toàn đủ năng lực đáp ứng ổn định cho nhu cầu sản xuất đến năm 2020.
  • Đề xuất định hướng chiến lược đầu tư mở rộng sang lĩnh vực bột giấy thương phẩm và chế biến gỗ sâu nhằm giảm tỷ lệ xuất khẩu thô dăm gỗ.
  • Mở ra tiền đề phát triển các mô hình dự báo kinh tế lượng đa biến và tự động hóa chuỗi cung ứng trong giai đoạn tiếp theo 2020 - 2025.

Đóng góp lớn nhất của luận văn là giải quyết hài hòa khoảng cách giữa kỹ thuật khai phá dữ liệu trong khoa học máy tính và bài toán quản trị sản xuất thực tế tại doanh nghiệp quy mô lớn. Doanh nghiệp cần khẩn trương số hóa công tác dự báo và triển khai áp dụng các mô hình định lượng vào hoạt động lập kế hoạch kinh doanh hàng năm.