Next Article in Journal
Detection of Growth Change of Young Forest Based on UAV RGB Images at Single-Tree Level
Next Article in Special Issue
Influence of Tree, Stand, and Site Attributes on Hardwood Product Yield: Insights into the Acadian Forests
Previous Article in Journal
Refined Assessment of Economic Loss from Pine Wilt Disease at the Subcompartment Scale
Previous Article in Special Issue
Gradeability of a Forwarder Based on Traction Performance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Benefits and Requirements of Mathematical Optimization in the Allocation of Wood to a Network of Forest Product Mills

1
Bureau de Mise en Marché des Bois, Ministère des Ressources Naturelles et des Forêts, Québec, QC G1H 6R1, Canada
2
FORAC Research Consortium, Université Laval, Québec, QC G1V 0A6, Canada
*
Author to whom correspondence should be addressed.
Forests 2023, 14(1), 140; https://doi.org/10.3390/f14010140
Submission received: 2 December 2022 / Revised: 6 January 2023 / Accepted: 8 January 2023 / Published: 12 January 2023

Abstract

:
Supply planning is a challenge for the forest industry in the context of natural forests characterized by heterogeneity among raw materials. Several mathematical models have been proposed in the literature to support forest planning, though few have been used by companies. The complexity of the natural environment and the expertise required to use these models limit their application. Nevertheless, these tools can significantly improve profitability. Three main elements were analyzed to assess benefits fostered by computer-optimized planning: (i) assessing the potential of implementing mathematical optimization in companies, particularly by pinpointing the additional resources necessary; (ii) determining the benefits of mathematical optimization to support planning decisions in an industrial context; and (iii) analyzing the impact of variation in information precision. LogiLab, an optimization software was used to find the optimal allocation of raw materials to an industrial network of five mills. The plan produced using optimization was compared to the plan generated by company personnel. The optimized plan generated a nearly 20% greater net profit than the current planning method. This difference was in part due to the more efficient allocation of raw materials to mills. It also highlighted numerous benefits, including a 50% reduction in the time required to produce plans. Furthermore, if attributed volume can be distributed among sawmills, it would be possible to save CAD 3.21/m3 using optimized plans, greatly exceeding initial costs related to system implementation.

1. Introduction

Forests cover a productive area of over 400,000 km2 and play an important role in Quebec’s economy [1]. In 2020, there were 53,000 direct jobs in the forest industry and exports totaled CAD 9.4 bn. Forest product companies are aiming for optimal use of available resources at a time when markets are increasingly competitive, and the supply of wood with desired characteristics is limited. The current context of market competition highlights the importance of optimizing a company’s supply system. The growing complexity of factors to consider makes forest planning increasingly difficult. These factors include requirements for raw material, increasingly dispersed harvesting areas, the diversity of final products to be generated, wood harvesting methods to use, variety of transportation methods, several stakeholders to consider on the territory, environmental regulations, and silvicultural prescriptions to apply. In the literature, various studies have investigated wood supply network planning; many optimization models with beneficial impacts on network profitability have been developed.
Morneau-Pereira et al. [2] noted that supply planning in the forest industry was carried out by managers using a trial-and-error method based on the latter’s experience and intuition. Heuristic approaches remain today’s standard. Numerous criteria must be considered, and many decisions must be made to plan forest actions efficiently, leaving little leeway for managers due to environmental requirements, social apprehension, available infrastructure, and raw material characteristics. In addition, managers must choose the type of transportation method and type of harvesting to use, decide the location and moment of harvesting, and decide to which mills harvested wood will go. All these criteria make the planners’ job difficult. Although optimization has been used to solve forest planning problems for several decades, the nature of problems has continued to develop, making the task of optimization increasingly complex and bringing with it a greater number of constraints [3,4]. While the literature has proposed many mathematical planning models, very few of these models are currently used in companies. We speculate that the complexity of the matter, a constantly changing natural environment, and the expertise necessary to use these models limit their application in an industrial context.
This project primarily aimed to analyze mathematical optimization’s potential to improve wood supply planning. To do this, a gap analysis was performed between current planning in companies and planning generated by optimization models. As both types of planning do not necessarily use the same resources and skills (for example, software, workforce, completion time, training, knowledge), costs and conditions necessary for using one or both approaches have been formally recognized and described. Three specific goals were established: (i) assess the potential of implementing mathematical optimization in companies; (ii) compare the benefits and advantages of using mathematical optimization versus the company’s current planning methods to support planning decisions; and (iii) analyze the impact of variation in forest inventory information used to prepare plans depending on one or the other planning methods under study.
This article is divided as follows. Section 2 provides a literature review on procurement network planning, forest planning of the study area, and various models that have been developed in relation to the problem being addressed. Section 3 presents the method used, case study, and decision-making tool used. Section 4 presents the results, including the requirements for the implementation of the optimization approach, the comparison of the two planning methods under study, and the analysis of gaps between a set of planning scenarios. Section 5 provides a conclusion that emphasizes the key components addressed in the research.

2. Literature Review

Mosconi [5] defined the wood supply network as the “people, organizations and activities from procurement to deliveries at processing mills” [unofficial translation]. Supply serves as the point of departure for creating a value network. According to Lehoux et al. [6], network performance essentially relies on alignment between supply steps and production tasks until final product sales. Planning operations within this network is carried out according to one of two main approaches used in network planning: centralized and distributed planning [7]. The centralized planning approach is based on a global view of the supply chain where management is coordinated by one parent company that considers all circulating flows [8]. Planning decisions are based on the entire supply chain in order to maximize global benefits [9]. In contrast, distributed planning considers independent companies when making decisions. Interactions among various actors in the chain are based on direct client–supplier relationships to ensure that the upstream supplier can meet the demand of the downstream client [10].
In network planning, numerous parameters can influence decision making. Planning supports decision making by developing scenarios the manager will analyze before choosing which option suits them best [11]. The planning process in the forestry sector is usually based on trial and error through planners’ experience and intuition [12]. Mosconi’s study [5] establishes the planning method used by a forest product company operated in Quebec. Mosconi listed various tools used by the company, including Excel spreadsheets, the geographic information system ArcGIS, simulation software FPinterface, and internal software which were generally used to correct and analyze various planning scenarios. LeBel et al. [13] state that the use of advanced technological planning tools would allow the rapid assessment of several scenarios in addition to visualizing impacts related to changes proposed by various stakeholders. However, according to study partners, the aforementioned initiatives were indeed attempted in the past and the use of tools was abandoned because of the models’ complexity, the time required, and the financial and human resources needed to use these tools.
Many studies have addressed computerized tools; in particular, in the 1990s [14], the development of simulation and optimization models greatly expanded. Thereafter, increasingly efficient models were developed to support decision making in planning forest supply activities (see: [2,15,16,17,18,19,20,21]. Table 1 summarizes the characteristics of a selection of optimization models and lays out the potential benefits of their use. As noted previously, many models were developed but not used concretely in companies. Developed models often include variables that make model use complex. In reality, these variables are often inappropriate for the planning to be applied and are not adapted to the industrial context.
Given the study’s context and goals, the LogiOpt planning system first developed by Morneau-Pereira [2] seems a useful foundation to develop new functionalities required for application in companies. This system can be used as a coupled approach with software currently used in companies and an optimization model. It could also promote the integration of supply to process, and process to market as suggested by [22,23]
Unlike previous studies that presented the same subject, this study establishes the conditions necessary to implement mathematical optimization, to improve and simplify the planning process based on planners’ experience and intuition. This study also quantifies the costs and benefits of implementing a planning approach using optimization.

3. Materials and Methods

The first component of our methodological approach is to compare an actual plan (as established by planners in a forest product company), with an optimal plan made by a decision-making tool. To assess the current planning approach, work carried out during the annual programming process (APP) was noted. Each step was observed and described to have an overall view of the method and requirements related to implementation. This documenting of the process used by companies made it possible to easily pinpoint implementation steps that could be optimized.
A case study was carried out featuring a large-scale Canadian forest product company supplying five sawmills in the same general area. These mills mostly transform three tree species into lumber of various sizes: black spruce (Picea mariana M.), balsam fir (Abies balsamea L.), and jack pine (Pinus banksiana L.). Annually, the largest of the five mills can process 1 million m3, and the smallest one, 0.48 million m3. The average mill capacity is 0.8 million m3. By-products such as shavings and sawdust were also considered, though they are less desirable because they reduce generated profits. Information about the network and products were used to model the problem. This case study used the LogiLab optimization platform [24]. The latter is based on a linear programming mathematical optimization model (see Appendix A) and can include many key factors to be considered. It also allows the consideration of transport, harvesting, heterogeneity of raw materials, and factors associated with sawmill performance when planning wood supply from the forest to in-plant processing. The goal is to maximize net value from the sale of finished products for the entire network. The solution may therefore generate a supply plan that is more costly to harvest and transport if it yields revenues that compensate for these extra costs.
A complete understanding of industry needs was required in order to apply an optimization model to the case study. We configured the generic mathematical model used in the LogiLab platform and considered the partner company’s specific context (for example, bucking rules, log assortments, and transportation modes). In so doing, we made changes to the generic model and added constraints to reflect the company’s context.
More specifically, we incorporated tree bucking patterns used depending on the forest harvesting area (two lengths are possible for logs: 500 cm and 560 cm–16′ and 18′, respectively, in local measurement units). Yet only two of five sawmills in the network were capable of transforming 560 cm logs [Equation (1)]. Forest harvest teams sorted harvested stems in three piles, based on species or diameter depending on procurement region; each pile could be sent to a maximum of two sawmills [Equation (3)]. Wood is procured from 66 harvesting units, generating 2.8 million m3. Wood transport was performed by non-standard (oversized) trucks or four-axle trailers. Volume transported by non-standard trucks originated only from certain areas and could only be sent to three of the five sawmills by way of logging roads [Equation (2)]. We added the following components to the model:
New Set
B U Set of harvesting blocks
S 18 p i W Set of the bucking process in 18’ logs (560 cm)
I Set of piles
P i Set of products on the pile i
B i B Set of harvesting blocks on an area to harvest
New Constant
η Total transformation capacities at mills accepting 18’ logs (560 cm)
New variable
L u w Binary variable indicating whether products from u are sent to w
Mathematical formulation
t T p P w S 18 p i u U γ p w Y t u w η (1)
e δ u   δ u + t T p P F e p t M L u w u B , w U , i I
(2)
w U L u w   2 u B i
(3)
In addition, the company must establish a few parameters applicable to the model (for example, fuel cost adjustments, harvested areas positioning, and type of harvesting). Once optimization is completed, a summary of results can be seen in the LogiLab web interface. It is also possible to export a detailed solution file in Excel, which permits the generation of tables and graphs to make results analysis easier. The Figure 1 flowchart lays out the various steps and files needed to create the model and visualize results.
The estimation of revenue generated via completed product sales was made possible using OptitekTM, a software developed by FPInnovations. For each network mill, OptitekTM simulated products generated for each log that entered the mill according to the current configuration of saw lines and log characteristics (for example, curvature, taper, and size). Quebec’s forest area is divided into 191 stumpage zones to reflect the market value of standing forest for each cubic meter (m3). This value reflects market conditions and is adjusted every three months [25]. Harvesting area data for company planning were produced using ten-year inventories from the Ministère des Forêts, de la Faune et des Parcs [government] (MFFP). Data were then processed using the FPInterfaceTM aggregation tool (a complementary FPInterfaceTM module) to create harvesting blocks that work as procurement points to allocate wood to each mill. Figure 2 shows an example of harvesting blocks created by aggregating polygons in shapefile by attribute. A total of 66 blocks were created. Each harvesting block became a grouping of spatially dispersed harvesting areas.
Managers were questioned about their specific tasks so as to integrate the new optimization process into their current company planning process. Based on the answers received, we created a diagram that summarized the planning process, which was validated by the partner company. Planning steps, task descriptions for each staff member, hours spent on each task, software programs used for the process, and the approximate cost of developing a plan were laid out. We pinpointed various alternative and realistic scenarios in collaboration with planners.
Seven scenarios were generated in total: four optimized scenarios, a manual scenario, and two scenarios with adjusted inventory data. The latter two scenarios assessed the impact of imprecision of data on manual and optimal plans. Table 2 shows the list of scenarios and the comparisons performed.
  • The “optimal baseline” (S1) plan was the closest to the current process in the partner’s company. It was the reference point for comparison to other scenarios.
  • The optimal scenario with “modified piles” (S2) compares the impact of sorting logs based on diameter instead of species. Piles tested in this scenario were (i) a pile of fir trees (all sizes), (ii) a pile of spruce and jack pine of small diameter (<16 cm), and (iii) a pile of spruce and jack pine of large diameter (≥16 cm).
  • The optimal scenario with “set bucking” (S3) requires each block to be harvested with a 500 cm (16’) length, regardless of destination. This scenario makes comparison with the reference scenario possible because two bucking lengths (500 and 560 cm) are allowed depending on the destination.
  • The “combined” optimal scenario (S4) uses pile sorting depending on diameter, species, and 500 cm bucking length.
  • The “manual” plan (S6) is where allocation decisions are made by planners without optimization.
  • Both scenarios with “adjusted inventory” (S5 and S7) consider actual data (exact values) of inventories that were measured and cross-referenced at the mill (available for some of the harvested blocks).
Lastly, in order to obtain comparable and realistic scenarios, we constrained the optimization model to at least meet the volume allocation specified on the company timber supply guarantee (TSG). This is to ensure each mill received at least 100% of its contracted wood supply. We also used a constraint to ensure optimal scenarios harvested the exact volume anticipated in the manual planning of the network, which was nearly 130% of all TSGs. The extra volume is usually obtained from an auction system managed by a state agency.

4. Results

4.1. Required Resources

Industry has the authority to act on a limited number of steps in the public forest planning process, as other steps are controlled by other agents. We identified seven main steps to describe the planning process:
  • Timber supply allocations (chief forester)
  • Choice of areas of interest (operational joint committee)
  • Volume available from auctions (Bureau de mise en marché des bois-BMMB)
  • Development of the individual annual programming process (IAPP) by timber supply (TS) beneficiaries
  • Sharing of individual annual programming processes (IAPP)
  • Submission of the IAPP to BMMB
  • Adaptation of harmonization measures if required
Of the seven steps, industry can provide input only on the fourth and seventh; thus, we proposed optimization for these two steps. The cost of the company’s manual plan was calculated using the hours needed per employee category for the development of each process step, and their average pay scale; necessary resources such as material and software were also considered. The process map (Figure 1) and interviews with company personnel allowed us to estimate labor time and cost.
As shown in Figure 3, 12 workers were involved, most of whom had a college degree. In total, between 56 and 128 workdays (eight hours per day) were needed for planning, amounting to an estimated cost between CAD 23,873 and CAD 30,918. The equipment needed for plan development included the use of many computerized tools, including Microsoft Excel, ArcGIS, SyGif_GT/integrated model in ArcMap, Microsoft Word, and Microsoft Access. Costs connected to the use of these software programs were not calculated in the plan production cost because the company uses most of these computerized tools daily.
The initial optimal plan development in LogiLab includes the creation of input files. Once files have been imported into LogiLab, the problem is set, and optimization can be run. Generally, the resolution time for scenarios included in our optimization model was no more than 15 min. The creation and analysis of scenarios via LogiLab required a professional who had a bachelor’s degree in forestry, which corresponds to university-level training. Time required to generate an optimal planning scenario with LogiLab, create an input file, and incorporate a detailed solution in the control panel was estimated at a day and a half by a process expert. This cost was estimated at CAD 1119 for the first year and was cut by half in the following years. Through the incorporation of optimization in the planning process within the company, we estimated that the time currently used for planning would be minimally cut in half, amounting to an average savings of over CAD 13,698.
In addition, this process required the use of many computerized tools and software programs such as FPInterfaceTM and OptitekTM developed by FPInnovations, ArcGis, and CPLEX. According to Volpé [26], the FPInterfaceTM license cost CAD 3500 in addition to the annual subscription of CAD 1150. The CPLEX license was estimated at CAD 3573. The cost of OptitekTM was not included in the planning process because it is already in use by the partner company. Savings of CAD 4355 were possible as of the first year of optimal planning implementation within the company (Table 3). In the following years, savings were CAD 8415.
The current manual planning process requires more total time than the optimized planning process. Thus, we estimate that it would be more profitable to adopt the optimized process with LogiLab as soon as it is deployed. Investments in workforce training would be necessary in order to make the leap from the manual plan to the optimal plan and to master specific aspects of LogiLab. However, these investments should be negligible, given current workers’ qualifications. Since basic software such as ArcGis and Microsoft Excel are already used in the current manual planning process, it is possible to choose the optimal planning process with little effort.

4.2. Using Initial Inventory Data

The optimal reference plan adapted for the case study is the plan with the highest net value, followed closely by the optimal plan with modified pile constraints, as shown in Figure 4. The manual plan has the lowest net value with a 19.26% difference as opposed to the optimal plan of reference with a difference of approximately CAD 9 M. Other monetary values are presented in variation in relation to the optimal scenario (net value was set at 1.000) because of confidentiality.
In the manual plan, forecast revenues for sawing products are the lowest, whereas revenues for by-products are the highest (Table 4). The highest revenues for sawing products are provided by the optimal plan with modified pile constraints; the lowest by-product revenues are related to the optimal plan of reference.
Though the manual method attempts to lessen expenses associated with transportation distances, it has the highest transportation costs, amounting to 1.7% more than transportation costs in the optimal plan of reference. With respect to production costs, the same trend is observed: the highest costs are generated by the manual plan, with 2.7% more than the lowest production costs in the optimal plan of reference (see Table 4).
The analysis shows little variation among the optimal plan scenarios. The optimal plan of reference is well adapted to the company’s current operations; that scenario generates the greatest benefits. This result originates essentially from a drop in transportation and mill production costs. As with the optimal scenario with modified piles, the optimal scenario of reference is less restricted by product types allowed into the piles, enabling a more efficient allocation of products, more in line with transformation processes at the sawmills and travel distances. This is why the scenario with modified piles has a slightly lower net value even though it generated the highest lumber revenues.
When comparing the manual plan to other plans, we see a difference in costs, revenues, and net value generated as opposed to optimal scenarios that present more appealing results. In our analysis, revenues from by-products (e.g., chips, sawdust) have no benefit for the company; because the value of these by-products is lower than costs, they are generally sold at a loss.
The major difference between net values in the optimal reference scenario and the manual scenario is mainly the result of the allocation of harvest areas and volume depending on transformation performance at the sawmills. The more forest products (logs) shipped to sawmills have characteristics akin to the specific capacities of a given mill, the less they produce undesirable by-products and the more they produce high-value products (lumber), which has a direct impact on the net value generated by the mill.
Since the goal of this study is to maximize the network’s net value, the use of the optimization model allowed to send harvested volumes to the best-adapted sawmill, as a priority, to increase overall generated net value, all the while complying with imposed restrictions on minimum deliveries to each mill (Figure 5).
Comparing the allocation of expected harvest blocks to manual and optimal scenarios of reference, very few blocks were chosen and allocated the same among scenarios (Figure 6). Both scenarios allocated only 11 harvest blocks (16%) the same way, i.e., the same blocks and same volume shipped to the same sawmills. In comparing the allocation of harvest blocks between the two scenarios, we see that:
  • 91% of allocated harvest blocks to sawmill 1 were the same, though proportions were different among scenarios; none of the blocks were allocated the same way (shipped volume proportions/different blocks); 9% of allocated blocks were totally different;
  • 26% of harvest blocks to sawmill 2 were the same; these blocks were allocated the same way among scenarios (74% were totally different);
  • 63% of harvest blocks to sawmill 3 were the same, though proportions were different and no blocks were allocated the same way (37% were totally different);
  • 32% of harvest blocks to sawmill 4 were the same, including only 3% that were allocated the same way (68% were totally different);
  • 45% of harvest blocks to sawmill 5 were the same but in different proportions; none were allocated the same way (55% were totally different).
Allocating harvest areas to the sawmilling network is more profitable when optimization is used because the model considers transportation distances, transportation costs by mode, mill performance by imported forest products, capacity of mill processing, characteristics of raw material processed, value of generated finished products, and processing costs per mill and product. Manual planning is limited to matters related to transportation costs and uses approximations to guide mill process performance.
Our results confirmed that decisions made by the company concerning log classes in the forest were appropriate, as benefits were higher when piles no longer had to be sorted by diameter but only by species for all areas and when 560 cm logs were allowed to be produced.

4.3. Using Verified Inventory Data

The impact of forestry inventory data accuracy on planning was assessed using actual volume as measured following harvest. For our case study, we received exact (real) harvested volumes of 18 out of 66 blocks, representing approximately 27% of all harvest blocks. As illustrated in Figure 7, adjustments made to expected inventory data using data measured at the mill correspond to a drop of 8.18% of the total harvest volume and 11.6% of the blocks’ total area. Adjusted volume for 18 harvest blocks was lowered by 22.44% as opposed to the expected volume before harvesting. The proportion of black spruce and jack pine increases, whereas fir tree volume decreases when comparing forecast data with actual data.
Initial inventory data were replaced by actual values in the mathematical model; a new optimized solution was calculated with an adjustment of volumes and areas. Results of the optimal and manual plans with real data were compared with plans using initial inventories (forecast), as shown in Figure 8. Trends similar to those in manual and optimal plans were seen when inventories were adjusted. In adjusted manual and optimal plans (with 27% of adjusted inventories), there was a drop in harvested forest volume and factory products, a drop in costs and revenues, and a drop in net value. Few differences were seen between optimal and manual plans with inventory data adjusted in comparison to initial plans. On the basis of these results, the optimal plan remains the “best plan” even when inventory data change. In comparison to reference plans, both adjusted plans dropped at about the same level (around 14% for the optimal plan and 13.6% for the manual plan). The net value generated by optimization was still higher than the net value in the manual plan (a difference of CAD 7.5 M).
According to these data, the net value generated by the optimal plan of reference with adjusted inventory dropped slightly more than the value generated by the manual plan with adjusted inventory. For each cubic meter that was overestimated, the net value computed by the optimal plan generated losses of around CAD 5.01/m3 when compared to the manual plan with adjusted inventory. However, all inventories—considering that the manual plan’s net value does not drop post-adjustment—would need to be 30% less for the manual plan to become a better option and generate an equivalent or higher net value. An error close to 10%, around 286,013 m3, in harvested volume related to inventory data is a drop of 14% in net value in the optimal plan and 13.6% in the manual plan. In any case, these results point toward the importance of improving the precision of pre-harvest inventories.

5. Conclusions

Planning wood supply for forest product mills is often done based on planners’ experience and intuition. However, wood supply planning requires awareness of numerous criteria. It is generally accepted that the use of optimization helps planners make better decisions. In a rigorous manner, our study quantifies the benefits of a computerized optimization system to plan forest-mill allocations. Further, the effort and costs of setting up an optimization tool in a company are estimated using the real context of a company based in Canada. The literature indicates that previously, many theoretical models have been proposed—with few or no documented examples of application in companies. Thus, our study wanted to verify whether or not the complexity of the problem, the constantly changing natural environment, and the expertise required to use these models were the limiting factors to their application in an industrial context. To verify our initial hypotheses, we (i) assessed the potential for implementing mathematical optimization in a company and (ii) demonstrated the benefits of using mathematical optimization to support planning decisions. The study was conducted with the collaboration of a large-scale forest product company with five sawmills distributed across a large territory.
Results demonstrated that adopting the optimal planning process can lead to savings as early as the first year after implementation. However, implementing the proposed planning process will require some effort and changes within the company. Among these requirements is the purchase of certain software programs to create databases needed to use LogiLab, including FPInterfaceTM, OptitekTM, or CPLEX (a mathematical solver). A minimum of training in mathematical optimization and general programming knowledge is all that is required from staff. The use of an optimal planning process as a complement to the current planning method should minimally decrease the time dedicated to planning by half.
The results of a manual plan and an optimized plan were compared to demonstrate the benefits of using mathematical optimization to support planning decisions. Results showed little variation among optimal plans that incorporated different operational constraints. The effect of these variations revealed that the proposed optimal planning process was robust to operational variations. Results even demonstrated that the net value anticipated by the developed optimal plan was 19.26% higher than the value developed by the current company method. This increase can be attributed to an improved allocation of harvest blocks to specific mills. This more efficient allocation would lower transportation costs and increase average transformation performance in the network sawmills, enabling a greater volume of lumber—the main source of revenue—in addition to lowering production costs. It was considered and validated with company planners that all 66 areas were accessible except during the thaw season. It is however possible that the optimized solution would generate periodic imbalance between the volume to be transported and local trucking capacity. Such a problem is frequently reported in the literature [27]. Attention must therefore be paid to this element when producing the bi-weekly transport schedule at the operational level. Our results also demonstrate that it is crucial to know in advance (by at least one year) which sectors will supply mills. With respect to the impact of varying the quality of inventory data, further research using a larger quantity of verified data would be needed.

Author Contributions

Conceptualization, L.L. and M.A.; methodology, M.A., L.L. and E.B.; software, M.A.; validation, M.A., L.L. and E.B.; formal analysis, M.A., L.L. and E.B.; investigation, M.A. and E.B.; resources, L.L.; data curation, M.A. and E.B.; writing—original draft preparation, M.A.; writing—review and editing, L.L. and E.B.; visualization, M.A.; supervision, L.L.; project administration, L.L.; funding acquisition, L.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by NSERC and the FORAC research consortium, grant number IRCPJ 545469.

Data Availability Statement

By agreement with the partnering company, data pertaining to manufacturing processes cannot be shared. Other data is available on demand.

Acknowledgments

The authors would like to thank the FORAC research consortium and its partners for their financial support of this project.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

LogiLab was developed by the FORAC Research Consortium. It offers a generic supply chain (SC) model based on mathematical programming, to model a forest products network (i.e., business units such as forests, sawmills, paper mills, wood pellet mills, customers, etc., transformation and storage processes as well as the transportation modes used within the SC). All constraints related to supply, production, storage, and transport capacities are considered, as well as customer demand for finished products. LogiLab finds the best possible exploitation of available resources (raw material, production, and storage capacity, etc.), by minimizing total cost, or by maximizing the generated profits. The mathematical formulation (profit maximization model) is presented below.
  • T   Set of time periods
  • U   Set of business units
  • K   Set of capacity types (production and storage)
  • W   Set of processes (production and storage)
  • W t u     W   Set of processes available at business unit u at period t
  • P   Set of products
  • E   Set of existing links (edges) between business units
  • δ u +     E   Set of edges incidents to u
  • δ u     E   Set of edges extending from u
  • Decision Variables
  • Y w u t   Quantity of process w used at business unit u at period t
  • D p u t   Quantity of product p sold at business unit u at period t
  • F p e t   Flow of product p on edge e at period t
  • Parameters
  • q k u t   Quantity of capacity type k available at business unit u at period t
  • f p e t l   Minimum flow of product p on edge e at period t
  • f p e t u   Maximum flow of product p on edge e at period t
  • f e t l   Minimum flow of all products on edge e at period t
  • f e t u   Maximum flow of all products on edge e at period t
  • c w   Cost of using process w
  • c p e t f   Cost of transporting product p on edge e at period t
  • l p e t   Transportation time of product p on edge e at period t
  • s w   Time required for manufacturing product p with process w
  • α p w   Quantity of input product p required by process w
  • γ p w   Quantity of output product p manufactured by process w
  • λ k u w   Quantity of capacity type k at business unit u consumed by process w
  • d p u t   Demand of product p at business unit u at period t
  • ρ p u t   Selling value of product p at business unit u at period t
  • Objective Function
m a x   t T ( w W ( p P d p u t > 0 ρ p u t D p u t w W t u c w Y w u t ) e E ( p P c p e t f F p e t ) )  
  • Constraints
t 1 T ( w W t u t 1 + s w = t γ p w Y w u t 1 ) + e δ u + ( t 2 T t 2 + l p e t 2 = t F p e t 2 ) w W t u α p w Y w u t e δ u F p e t D p u t = 0
  t T ,     u U ,     p P
D p u t   d p u t   t T ,     u U ,     p P
w W t u λ k u w Y w u t   q k u t   t T ,     u U ,     k K
f e t l     p P F p e t     f e t u     t T ,     e E
f p e t l     F p e t     f p e t u     t T ,     e E ,     p P  
F p e t ,   D p u t ,   Y w u t   0   t T ,     e E ,     p P ,     w W

References

  1. Ministère des Forêts, de la Faune et des Parcs. Ressources et Industries Forestières du Québec. Portrait Statistique 2020; Gouvernement du Québec: Québec, QC, Canada, 2021.
  2. Morneau-Pereira, M.; Arabi, M.; Ouhimmou, M.; Gaudreault J et Nourelfath, M. An optimization and simulation framework for integrated tactical planning of wood harvesting operations, wood allocation and lumber production. In Proceedings of the MOSIM14, 10th International Conference on MOdeling, Optimization and SIMlation, Nancy, France, 5–7 November 2014. [Google Scholar]
  3. Rönnqvist, M. Optimization in forestry. Math. Program Ser. 2003, 97, 267–284. [Google Scholar] [CrossRef]
  4. Rönnqvist, M.; D’Amours, S.; Weintraub, A.; Jofre, A.; Gunn, E.; Haight, R.G.; Romero, C. Operations research challenges in forestry: 33 open problems. Ann. Oper. Res. 2015, 232, 11–40. [Google Scholar] [CrossRef]
  5. Mosconi, E.P. L’apport D’une Approche Multi-Dimensionnelle de Gestion des Connaissances à la Performance Décisionnelle des Organisations. Ph.D. Thesis, Université Laval, Quebec City, QC, Canada, 2014; 230p. [Google Scholar]
  6. Lehoux, N.; Marier, P.; D’Amours, S.; Ouellet, D.; Beaulieu, J. Value Creation Network of Canadian Wood Fibre; Centre Interuniversitaire de Recherche sur les Réseaux D’entreprises, la Logistique et le Transport (CIRRELT): Montreal, QC, Canada, 2015; T-2012-34; 210p. [Google Scholar]
  7. Somers, F. Hybrid: Unifying centralised and distributed network management using intelligent agents. In Proceedings of the Network Operations and Management Symposium, Kyoto, Japan, 15–19 April 1996; pp. 34–43. [Google Scholar] [CrossRef]
  8. François, J. Planification des Chaînes Logistiques: Modélisation du Système Décisionnel et Performance. Ph.D. Thesis, Université de Bordeaux, Bordeaux, France, 2007. [Google Scholar]
  9. Sarrazin, F.; LeBel, L.; Lehoux, N. Regional logistics centre: Application to Eastern Canada. INFOR Inf. Syst. Oper. Res. 2021, 59, 439–464. [Google Scholar] [CrossRef]
  10. Eriksson, M.; Lebel, L.; et Lindroos, O. Management of outsourced forest harvesting operations for better customer-contractor alignment. For. Policy Econ. 2015, 53, 45–55. [Google Scholar] [CrossRef]
  11. Boukherroub, T.; D’amours, S.; Rönnqvist, M. Sustainable forest management using decision theaters: Rethinking participatory planning. J. Clean. Prod. 2018, 179, 567–580. [Google Scholar] [CrossRef]
  12. Arabi, M.; Gaudreault, J.; Nourelfath, M.; Favreau, J.; Morneau-Pereira, M. Integrating optimization and simulation for supply chain tactical planning in the forest Products Industry. In Proceedings of the International Conference on Information Systems, Logistics and Supply Chain (ILS), Quebec City, QC, Canada, 26–29 August 2012. [Google Scholar]
  13. LeBel, L.; Laliberté, F.; Simoneau, L.; Morin, F.; et Gharbi, C. Planification Forêt Efficiente (Phase 1), Rapport Final; Consortium de Recherche FOR@C, Groupe Initiative Mauricie (GIM): Mauricie, QC, Canada, 2015; pp. 34–43. [Google Scholar]
  14. Messier, C.; Bernier, P.Y.; Lo E et Coates, D. Le développement des modèles de simulations en foresterie: En pleine expansion. L’Aubelle 2001, 136, 10–13. [Google Scholar]
  15. Karlsson, J.; Rönnqvist, M.; et Bergström, J. An optimization model for annual harvest planning. Can. J. For. Res. 2004, 34, 1747–1754. [Google Scholar] [CrossRef]
  16. Beaudoin, D.; LeBel, L.; et Frayet, J.-M. Tactical supply chain planning in the forest products industry through optimization and scenario-based analysis. Can. J. For. Res. 2007, 37, 128–140. [Google Scholar] [CrossRef] [Green Version]
  17. Marques, A.F.; de Sousa, J.P.; Rönnqvist, M.; et Jafe, R. Combining optimization and simulation tools for short-term planning of forest operations. Scand. J. For. Res. 2014, 29 (Suppl. S1), 166–177. [Google Scholar] [CrossRef]
  18. Gautam, S.; Lebel, L.; Beaudoin, D.; et Simard, M. Modelling hierarchical planning process using a simulation-optimization system to anticipate the long-term impact of operational level silvicultural flexibility. Sci. IFAC Pap. 2015, 48, 616–621. [Google Scholar] [CrossRef]
  19. Eyvindston, K.; et Kangas, A. Stochastic goal programming in forest planning. Can. J. For. Res. 2014, 44, 1274–1280. [Google Scholar] [CrossRef]
  20. Mobtaker, A.; Ouhimmou, M.; Rönnqvist, M.; Paquet, M. Development of an economically sustainable and balanced tactical forest management plan: A case study in Quebec. Can. J. For. Res. 2018, 48, 197–207. [Google Scholar] [CrossRef] [Green Version]
  21. Mobtaker, A.; Ouhimmou, M.; Audy, J.F.; Rönnqvist, M. A review on decision support systems for tactical logistics planning in the context of forest bioeconomy. Renew. Sustain. Energy Rev. 2021, 148, 111250. [Google Scholar] [CrossRef]
  22. Dumetz, L.; Gaudreault, J.; Bril El-Haouzi, H.; Thomas, A.; Lehoux, N.; Marier, P. Tactical-operational coordination of a divergent production system with coproduction: The sawmilling challenge. INFOR Inf. Syst. Oper. Res. 2021, 59, 377–399. [Google Scholar] [CrossRef]
  23. Hoogstra-Klein, M.A.; Meijboom, K. A qualitative exploration of the wood product supply chain—Investigating the possibilities and desirability of an increased demand orientation. For. Policy Econ. 2021, 133, 102606. [Google Scholar] [CrossRef]
  24. Boukherroub, T.; LeBel, L.; Lemieux, S. An integrated wood pellet supply chain development: Selecting among feedstock sources and a range of operating scales. Appl. Energy 2017, 198, 385–400. [Google Scholar] [CrossRef]
  25. Bureau de Mise en Marché du Bois. Grille de Taux de la Valeur Marchande des Bois sur Pied, Disponible en Ligne. 2019. Available online: https://bmmb.gouv.qc.ca/publications-et-reglements/tarification-forestiere/ (accessed on 24 February 2021).
  26. Volpé, S. Modèle D’évaluation de la Biomasse Forestière, FPinnovation. 34p. 2021. Available online: http://partenariat.qc.ca/wp-content/uploads/2019/02/FPI-Biomass-model_webinar_Feb-28_2019_FR.pdf (accessed on 10 December 2020).
  27. Malladi, K.T.; Sowlati, T. Optimization of operational level transportation planning in forestry: A review. Int. J. For. Eng. 2017, 28, 198–210. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the procurement optimization process used in this project.
Figure 1. Flowchart of the procurement optimization process used in this project.
Forests 14 00140 g001
Figure 2. An example of harvesting blocks before (left) and after aggregation by area (right) with the aggregation tool FPInterfaceTM.
Figure 2. An example of harvesting blocks before (left) and after aggregation by area (right) with the aggregation tool FPInterfaceTM.
Forests 14 00140 g002
Figure 3. Education received by employees in charge of planning.
Figure 3. Education received by employees in charge of planning.
Forests 14 00140 g003
Figure 4. Relative variation of net value with scenarios in comparison to the optimal scenario.
Figure 4. Relative variation of net value with scenarios in comparison to the optimal scenario.
Forests 14 00140 g004
Figure 5. Proportion of volumes allocated to factories in the optimal scenario of reference in comparison to the manual scenario.
Figure 5. Proportion of volumes allocated to factories in the optimal scenario of reference in comparison to the manual scenario.
Forests 14 00140 g005
Figure 6. Comparing allocation of harvest blocks in mills between manual and optimal scenarios of reference.
Figure 6. Comparing allocation of harvest blocks in mills between manual and optimal scenarios of reference.
Forests 14 00140 g006
Figure 7. Variation between real-time inventory and initial inventory.
Figure 7. Variation between real-time inventory and initial inventory.
Forests 14 00140 g007
Figure 8. Impact of changing inventory data on the optimal (a) and manual (b) plans.
Figure 8. Impact of changing inventory data on the optimal (a) and manual (b) plans.
Forests 14 00140 g008
Table 1. Elements considered and benefits of relevant models.
Table 1. Elements considered and benefits of relevant models.
AuthorsForest
Management
Harvest
Operations
Road
Network
TransportProcessing (Sawmill)Forest
Inventory
Benefits
Karlsson et al. (2004) [15] ↓ Costs
Beaudoin et al. (2007) [16] ↑ Profits
Marques et al. (2014) [17] ↓ Time
Gautam et al. (2014) [18] ↑ Profits
Eyvindson &
Kangas (2014) [19]
↑ Profits
Morneau-Pereira et al. (2014) [2] ↑ Profits
Mobtaker et al. (2018, 2021) [20,21] ↓ Costs
Table 2. Comparison matrix for the experiment.
Table 2. Comparison matrix for the experiment.
ScenarioS1S2S3S4S5S6S7
S1. Optimal (baseline)
S2. Optimal with modified piles
S3. Optimal with set bucking rules
S4. Optimal combined (modified piles and set bucking rules)
S5. Optimal with adjusted inventory
S6. Manual
S7. Manual with adjusted inventory
Table 3. Comparison between current planning costs and optimized planning costs for the first year of implementation and following years.
Table 3. Comparison between current planning costs and optimized planning costs for the first year of implementation and following years.
Optimized Planning Process
ExpensesYear 1Year 2 and onCurrent Method
Company personnelCAD 14,817 CAD 14,258 CAD 27,395
Software licenseCAD 8223 CAD 4723 -
TotalCAD 23,040 CAD 18,981 CAD 27,395
Table 4. Revenues and costs variations for selected scenarios in comparison to the optimal scenario (set at 1.000).
Table 4. Revenues and costs variations for selected scenarios in comparison to the optimal scenario (set at 1.000).
Modified Piles (S2)Set Bucking (S3)Combined
(S4)
Manual
(S6)
Expected revenues from the sale of lumber1.0000.9930.9930.975
Expected revenues from the sale of by-products1.0001.0001.0011.010
Transportation costs1.0011.0051.0061.017
Production costs1.0001.0001.0001.027
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Auger, M.; LeBel, L.; Brotherton, E. Benefits and Requirements of Mathematical Optimization in the Allocation of Wood to a Network of Forest Product Mills. Forests 2023, 14, 140. https://doi.org/10.3390/f14010140

AMA Style

Auger M, LeBel L, Brotherton E. Benefits and Requirements of Mathematical Optimization in the Allocation of Wood to a Network of Forest Product Mills. Forests. 2023; 14(1):140. https://doi.org/10.3390/f14010140

Chicago/Turabian Style

Auger, Maxime, Luc LeBel, and Edith Brotherton. 2023. "Benefits and Requirements of Mathematical Optimization in the Allocation of Wood to a Network of Forest Product Mills" Forests 14, no. 1: 140. https://doi.org/10.3390/f14010140

APA Style

Auger, M., LeBel, L., & Brotherton, E. (2023). Benefits and Requirements of Mathematical Optimization in the Allocation of Wood to a Network of Forest Product Mills. Forests, 14(1), 140. https://doi.org/10.3390/f14010140

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop