Abstract

Supply chain optimization for biomass-based power plants is an important research area due to greater emphasis on green energy sources. This paper develops and applies two optimization models to analyze the impacts of biomass competition on cost structures and gross margins for four competing biomass-based power plants in northwestern Ontario. Model scenarios are run to study the impacts of changes in parameters relevant to biomass type and processing technology, and prices of inputs and outputs on procurement costs. Cost minimization model shows that per unit procurement costs are directly proportional to the size of the power plants in all scenarios. Profit maximization model, on the other hand, shows that FMUs that are closer to the power plants make higher gross margins. However, the margins significantly increase for FMUs that are close to the power plants potentially offering higher prices. The variations in costs and gross margin structures under various model scenarios are explained by location of depletion cells relative to power plants, availability of each type of biomass in depletion cells, biomass demands, and differential processing costs for two types of biomass. These results can aid decision makers to make improved decisions related to biomass supply chains for bioenergy production.

1. Introduction

Forest biomass has been recognized as an alternative energy source since it is renewable and CO2 neutral [1]. However, renewable energy production from forest biomass faces many challenges due to uncertainty of its continuous supply [24]. The ever increasing demand of biomass for bioenergy production has enormously increased the transportation distances and costs for forest biomass procurement, which is spread over large geographical locations at varying distances from the power plant [1, 5]. Further, the optimal harvest schedules become complicated if more power plants compete for available biomass feedstock over given space and time. This suggests a need for developing and using optimization models to analyze such problems for improved decision making in bioenergy production.

Canada’s dependence on fossil fuel has changed in recent times, and forest biomass has become an important part of its energy picture, supplying about 4.7% of our primary energy demand, the second largest source of renewable energy after hydroelectricity [6]. A major application of bioenergy is found in the forest products industry of Canada. Beyond the forest industry, several independent power plants generate electricity from forest biomass in Canada. Currently, three major combined heat and power (CHP) plants in northwestern Ontario (NWO) are using forest biomass feedstock. These include Thunder Bay CHP plant (with annual feedstock demand of 730,000 green tonnes (gt)), Fort Frances CHP plant (800,000 gt), and Dryden CHP plant (480,000 gt). Atikokan power-generating station, another power plant in NWO, is currently being converted to use forest biomass feedstock instead of coal requiring 200,000 gt annually [7, 8]. When all these plants become fully operational, the power plant managers (buyers of biomass) will need to have a sustainable and cost effective supply of biomass, and the contractors (suppliers of biomass) would like to increase their gross margins (profits) from biomass supply operations. Therefore, the optimization of forest biomass feedstock supply needs to be addressed both with cost optimization and profit maximization points of views.

Forest harvest residue (FHR), which includes tops and branches and unmerchantable wood left after stand harvesting, and underutilized wood (UW), which includes unharvested tree species that are not commercially important for timber, are the two most common forms of forest biomass used for bioenergy. Further, there are numerous options for comminuting forest biomass and several trucking and loading options having varying costs. Optimizing biomass procurement is, therefore, a complex problem with several supply and demand constraints.

A number of computational techniques from heuristics to advanced optimization models have been developed to model location specific biomass procurement problems [914]. These studies focus on bioenergy plant location problems for widespread forest biomass applications. The supply network normally consists of a single bioenergy production plant with many biomass supply regions [1517]. In this context, forest biomass transportation costs constitute the majority of forest fuel supply chain costs for energy production. Further studies in biomass procurement, therefore, focused on reducing overall transportation costs [1821]. Linear programming (LP) models, which are a common technique in operations research, have been used to optimize transportation distances and costs for forest biomass procurement for energy production at regional levels [22, 23]. Forest biomass feedstock data have also been analyzed using geographic information systems (GISs) techniques [24, 25]. LP models have also been developed using GIS-based forest biomass data both for estimating the feedstock availability and for reducing transportation distances and costs [2629]. However, all these studies have been conducted either in Nordic or in central European countries.

We found no such optimization models dealing with biomass supply competition among various power plants in the Canadian context with large contiguous wood supply areas. The purpose of this study is, therefore, to develop and apply optimization models to analyze the impact of interpower plant competition for the available biomass feedstock on cost structures and gross margins of each power plant in NWO. Specifically, the two models (one based on cost minimization and another on profit maximization) are developed. The cost minimization model is applied from the power plant managers’ (buyers) perspective, and the profit maximization model is applied from the biomass suppliers’ perspective. Several sensitivity analyses are run to study the impact of changes in parameters relevant to forest biomass processing technology, biomass types, and prices of inputs and outputs.

2. Methods

The flow chart in Figure 1 presents the methodology for developing the LP models for supplying forest biomass feedstock to the four power-generating stations. The study area, 324 km × 516 km (167,184 km2), consists of 18 FMUs west of Lake Nipigon in NWO. GIS data are collected from Land Information Ontario, Sustainable Forest Licence (SFL) holders, and consultant companies in the formats of Shapefile and Geodatabase [3033]. The original vector data is first converted to raster and finally to spatial database text files for the entire research area using ArcGIS software. Three main spatial layers (land use, forest depletion, and cost layers) are prepared on a raster grid size of 1 km × 1 km (1 km2), where each cell represents a feature [34, 35]. Different data input methods are followed for different types of features. The dominant data input method is used to create a raster layer of forest land use class [35], where a grid code entity is assigned to each feature (productive forest, water/lake, and other land use), which occupies more than 50% of the cell area (detailed methodology of data input methods is provided by Alam et al. [36]). The percent occurrence method is used to prepare the depletion layer, where a code number is assigned to the grid cell depending on its depletion percentage. Before developing the marginal transport cost layer for each power plant, a raster layer of the road network for the research area is also developed using the presence/absence method of data input [35]. A grid code is assigned depending on the presence of different types of roads in a cell. The raster layer of the road network is converted to ASCII (.txt) format. A network optimization model is used to prepare the minimum cost layer for transporting forest biomass through the road network to the four power plants. The location of each power plant is used as the sink node during the data input and optimization processes. The vehicle considered for transporting forest biomass feedstock is a tractor with a 53-foot semitrailer (with belly). In Ontario the allowable legal payload weight for this type of tractor-trailer is taken as 40.55 tonnes (t). The charge-out rate for biomass truck with operator is assumed to be 85 $ ·h−1. A fixed time for loading, unloading, and delay of 2.5 hours per trip is assumed. In this cost layer, the cost of transporting forest biomass from each road cell (1 km2) of the research area to the four power plants is established. Minimum cost zones (tessellations) for transporting forest biomass from each of the 100,061 forest cells to each of the four power plants are established in order to create biomass catchment areas for each power plant.

Two optimization models are developed for analyzing optimal supply of forest biomass from depletion cells to the four power plants. The objective function in the first model is to minimize the total biomass feedstock procurement costs and in the second model is to maximize total profit (gross margin) of FMUs for supplying forest biomass feedstock from the forest cells to power plants. Even though these two models give the same solutions, the relative impacts of changing biomass prices are not captured by the cost minimization model. Changing biomass prices on part of bigger power plants will change the optimal solutions, which is best handled by the profit maximization model. The results of profit maximization are more relevant to the contractors (biomass suppliers) in order to evaluate the profitability of their operations in the given locations (FMUs). The parameters, indices, and variables used for developing the two models using the General Algebraic Modeling System (GAMS) programming language are listed as follows.

Indices:
𝑖: power plants (𝑖: 1 = Thunder Bay CHP plant, 2 = Fort Frances CHP plant, 3 = Atikokan generating station and 4 = Dryden CHP plant),𝑗: forest cell in the depleted forest area (𝑗=1,2,3,,19,315).

Parameters:
HF: harvesting factor,PR: processing (harvesting and grinding/chipping) cost ($gt−1) of FHR at roadside, PU: processing (harvesting and grinding/chipping) cost ($gt−1) of UW at roadside, DB𝑖: annual forest biomass demand (gt) of each power plant, TR𝑗: theoretical availability (gt) of FHR in a forest cell over 7 years (2002–2009), TU𝑗: theoretical availability (gt) of UW in a forest cell over 7 years (2002–2009),AR𝑗: annual technical availability (gt) of FHR in a forest cell, calculated using (1) AR𝑗=TR𝑗×HF7,(1)AU𝑗: annual technical availability (gt) of UW in a forest cell, calculated using (2) AU𝑗=TU𝑗×HF7,(2)TC𝑖𝑗: biomass transportation cost ($gt−1) from the 𝑗th forest cell to the 𝑖th power plant including loading and unloading costs. 𝑃𝑖: price ($gt−1) of forest biomass at the 𝑖th power plant.

Variables:
𝑅𝑖𝑗: amount of annual FHR harvested (gt) from the 𝑗th forest cell for the 𝑖th power plant, 𝑈𝑖𝑗: amount of annual UW harvested (gt) from the 𝑗th forest cell for the 𝑖th power plant,TB𝑖: annual forest biomass (gt) brought in the 𝑖th power plant, TCB: total annual cost ($) of biomass procurement, calculated using (3) TCB=4𝑖=119315𝑗=1𝑅𝑖𝑗PR+TC𝑖𝑗+4𝑖=119315𝑗=1𝑈𝑖𝑗PU+TC𝑖𝑗,(3)TRB: total annual revenue ($) from forest biomass supply to power plants, calculated using (4) TRB=4𝑖=119315𝑗=1𝑃𝑖𝑅𝑖𝑗+𝑈𝑖𝑗,(4)TPB: total annual profit ($) from forest biomass supply to power plants, calculated using (5) TRB=4𝑖=119315𝑗=1𝑃𝑖𝑅𝑖𝑗+𝑈𝑖𝑗4𝑖=119315𝑗=1𝑅𝑖𝑗PR+TC𝑖𝑗4𝑖=119315𝑗=1𝑈𝑖𝑗PU+TC𝑖𝑗.(5)

The objective functions of two LP models, which minimize overall procurement costs and maximize total profit, are specified by (6) and (7), respectively. The constraints used in these models are specified by (8) to (12). MinimizeTCB=4𝑖=119315𝑗=1𝑅𝑖𝑗PR+TC𝑖𝑗+4𝑖=119315𝑗=1𝑈𝑖𝑗PU+TC𝑖𝑗.(6)MaximizeTPB=4𝑖=119315𝑗=1𝑃𝑖𝑅𝑖𝑗+𝑈𝑖𝑗4𝑖=119315𝑗=1𝑅𝑖𝑗PR+TC𝑖𝑗4𝑖=119315𝑗=1𝑈𝑖𝑗PU+TC𝑖𝑗,(7) subject to 19315𝑗=1𝑅𝑖𝑗+𝑈𝑖𝑗=TB𝑖,(8)4𝑖=1𝑅𝑖𝑗AR𝑗,(9)4𝑖=1𝑈𝑖𝑗AU𝑗,(10)TB𝑖DB𝑖,𝑅(11)𝑖𝑗,𝑈𝑖𝑗0.(12)

Equation (8) is an intermediate equation summing up the plantwise total biomass being harvested. Equations (9) and (10) are the FHR and UW harvesting constraints respectively, which specify that the total amount of FHR and UW harvested from 𝑗th forest cell is less than or equal to the annual amount of technically available FHR and UW in the 𝑗th forest cell, respectively. Equation (11) is the demand constraint, which ensures that the total amount of forest biomass harvested for 𝑖th power plant should at least be equal to biomass feedstock demand for the power plant. Finally, (12) is nonnegativity constraint on decision variables. The estimates of parameters used in the model are presented in Table 1.

Model Scenarios
Table 2 presents the definitions and descriptions of 14 cost minimization model scenarios with their relevancy in relation to the expected changes in technological and economic parameters in the model. Out of 14 scenarios, six scenarios (DEPC, IHRW1, IHRW2, and IHFR) are cost decreasing and seven scenarios (IHOR, DHRW, INPC, DEPC, INDP, DEDP, and INTC) are cost increasing. Table 3 shows the definitions and descriptions of eight different model scenarios for the profit maximization model. The results of this model are presented in terms of per unit gross margins for each FMU by scenarios.

While selecting the model scenarios, first we focus on impacts of changing economic parameters (transportation and processing costs and prices) on cost and gross margin structures. Secondly, the impacts of changing technological parameters (harvesting factors and biomass demand and availability) are investigated.

3. Results and Discussion

3.1. Cost Zones

Minimum transportation cost zones (tessellations) for transporting forest biomass from each forest cell of depleted forest area to the four power plants are shown in Figure 2. These minimum cost zones show the forest biomass catchment areas for each power plant. The minimum cost zone for the Dryden CHP plant is the largest, while the Thunder Bay CHP plant catchment area is the second largest among the four cost zones. The Dryden CHP plant is located near the middle of the research area, and there are many forest cells closer to this power plant along with a denser network of higher class and straighter roads. There is no competition for the Dryden CHP plant for forest biomass supply from any other power plant in the research area on the northern side. Similarly, Thunder Bay CHP plant is located in the southeastern part of research area. There is no competing power plant on its northern, eastern, and southern sides in the research area. The other two power plants, Atikokan and Fort Frances, which are located close to each other in the research area, might compete for the same forest biomass from some depletion cells depending upon their requirements. The cost zone map in Figure 2 shows the relative cost competitive areas for each power plant based on minimum transportation cost for each power plant from each forest cell. As more attributes (biomass density and economic) are introduced in the forest cells, the dynamics of competition among the four power plants would change.

3.2. Cost Analysis

The results of the cost minimization model provide an optimal solution for supplying forest biomass feedstock from forest cells to the four power plants on an annual basis by minimizing the total annual procurement (harvesting, processing, and transportation) costs of biomass, subject to the availability of forest biomass in each depleted forest cell and demand of each power plant. The model selects 11,790 and 2,991 cells in total for supplying FHR and UW, respectively, for all four plants. The more FHR cells are selected because the processing cost for FHR is less as compared to UW forest biomass. Figure 3 represents the distribution of optimal depletion forest cells that supply FHR and UW biomass to each power plant in the BASE scenario. Similar distributions of optimal depletion cells for different model scenarios can be produced; these are not reproduced here because of space limitations.

First, we discuss the results relating to the impacts of changing economic parameters (transportation and processing costs and prices) on cost structures. Results of the impacts of changing technological parameters (harvesting factors and biomass demand and availability) are discussed next. The number of cells selected for supplying FHR and UW forest biomass for each power plant, in each model scenario, is presented in Table 4. While observing number of selected cells in Table 4, we find more depletion cells for FHR in all scenarios for all plants, but lesser amount of biomass from FHR type is being harvested (Table 5) in most of the scenarios. This variation is explained by the fact that, on an average, the availability of FHR is about one-fourth of UW biomass type. The source of the information is the input dataset of the model.

Table 5 describes the distribution of biomass (both FHR and UW) amount being harvested from the optimal depletion cells for each plant under different model scenarios. The variations in number of depletion cells (Table 4) and biomass amount (Table 5) being selected by power plants are explained by amount of biomass requirement for each plant, availability of FHR and UW in each cell, relative processing costs of each biomass type, and the competition among the power plant for the biomass in the given cell. Further, these variations in distribution of optimal depletion cells and amount of biomass for both types explain the per unit biomass procurement cost structures for each power plant under different model scenarios (Table 6).

The 13 cost minimization model scenarios analyze the impacts of changing relevant economic (three scenarios) and technological (10 scenarios) parameters (as described in Table 2) on biomass procurement costs structures of the four power plants as shown in Table 6. The results of these scenarios are compared with that of the BASE scenario. The results in Table 6 show that per unit costs for Fort Frances, Thunder Bay and Dryden CHP plants, and Atikokan power plant are in descending order for all model scenarios, which indicate that the higher the biomass demand, the higher the per unit procurement cost. The reason for this is that power plants with more biomass demand have to collect biomass from much longer distances. The range of unit procurement cost ($ODt−1) for Fort Frances, Thunder Bay and Dryden CHP plants, and Atikokan power plant is 60.81–76.81, 59.72–75.37, 56.57–71.31, and 56.03–70.39, respectively. This result shows that the highest per unit procurement cost (for all power plants on average a 16.77% increase from the BASE scenario) occurs with a 20% increase in transportation cost (INTC scenario). The lowest per unit procurement costs occurs with a 10% decrease in processing cost scenario (DEPC scenario). A decrease in processing cost by 10% resulted in a 7.44% decrease in per unit procurement costs for all power plants on average. Since transportation and processing costs form a major component of biomass supply chain, changes in these parameters have major impacts on per unit procurement cost for all power plants. These sensitivity results imply that policy should focus first on controlling transportation costs and then processing costs. The processing cost includes both productive machine time and idle time costs. The productive machine time can be improved by investing in improving the processing technology, and the idle time can be reduced by better scheduling of operations. Improving road network and vehicle configuration for hauling biomass can reduce the transportation costs. Consequently, one can achieve reduced procurements costs in biomass value chains with these improved measures.

After analyzing the sensitivities of changing economic parameters to the model solutions, we investigate further the impacts of changing relevant technological parameters in order to evaluate the resulting cost structures due to potential changing conditions of these parameters in future. In model scenarios INBA, IHRW1, IHRW2, IHFR, and DEDP, the per unit procurement costs decreased from the BASE scenario for all power plants with the highest decrease at the Atikokan power plant (Table 6). On the other hand, under the model scenarios DEBA, DHRW, DHFR, IHOR, and INDP the per unit procurement costs increased from the BASE scenario for all power plants with the highest increase at the Atikokan power plant. The explanation for these higher changes in per unit procurement costs for the Atikokan power plant is that it has the lowest biomass demand compared to other power plants; it is located near the favorable depletion cells and biomass availability (Figures 2 and 3). However, in IHOR (only FHR biomass with 20% more harvesting factor) scenario, the increase in per unit cost for Atikokan power plant is the least compared to others. Although more FHR biomass is available in the IHOR scenario, there is higher competition for this type of biomass because UW is not being collected, and hence there is an increase in per unit cost compared to the BASE scenario for all power plants. In this scenario, the Atikokan power plant has least increment in per unit cost due to its lower demand for biomass relative to the other power plants. After the Atikokan plant, Fort Frances, the biggest power plant, shows higher increases and decreases in per unit procurement costs for most of the scenarios (Table 6), albeit with small variations among the plants by scenarios. The more interesting results are found in the IHOR scenario, where Dryden CHP plant has the highest increase (9.97%) in per unit cost compared to the BASE scenario. Further, the aggregate impacts of the economics and technological model scenarios are reflected in Figure 4, where we present the power plantwise total biomass procurement costs by scenarios. The bars in Figure 4 contrast with the results of per unit costs by scenarios. For example, the highest total costs for each power plant are shown in the INDP scenario, whereas the highest per unit cost is found in the INTC scenario.

Increasing harvesting factors causes more declines in per unit procurement costs suggesting that the present harvesting factor level needs to be increased for reduced per unit procurement cost. This is possible given that the depletion cells being used in this study are only from productive forest area, and it excludes the ecologically sensitive areas in NWO. Further, IHOR scenario showing us the highest per unit procurement cost, which suggests current emphasis on collecting only FHR for bioenergy purpose, would lead to higher per unit procurement cost. This, therefore, suggests for harvesting standing underutilized tree species to reduce per unit procurement cost.

3.3. Profit Analysis

The results of the profit maximization model provide an optimal solution for supplying forest biomass feedstock from forest cells to the four power plants on an annual basis by maximizing the total gross margin, subject to the availability of forest biomass in each depleted forest cell and meeting the biomass demand of each power plant. These results of gross margin structures are relevant from the biomass supplier’s (contractors) perspective as these are based on biomass supplying FMUs (Table 7), where the contractors operate their business. The results of per unit profit ($ODt−1) for each FMU that supplies biomass to the four power plants in the BASE and seven other scenarios are shown in Table 7. The profit for any FMU depends on the price of biomass being offered by the power plants. In the BASE scenario, most of the FMUs make positive profit.

Although 19,315 depletion cells used in this study are from 18 FMUs, the optimal depletion cells selected by the model fall only in 13 FMUs for all the scenarios. In the BASE scenario, only the Kenora and Lake Nipigon FMUs show negative gross margins (Table 7) as the optimal depletion cells falling in these FMUs are relatively farther away from the power plants (Figure 2). While testing the sensitivity of increasing truck charge rate by 20% (INTC scenario), only four FMUs (Dog River-Matawin, Dryden, English River, and Lac Seul FMUs) make positive profit. These results suggest that, in future, contractors operating in most of the FMUs will face negative gross margins, thereby halting the biomass supply hindering smooth operation of the power plants. We, therefore, tested the sensitivity of increased per unit price on part of the larger power plants to the resultant gross margin structures of biomass supplying FMUs. For this purpose, a combination of price increase scenarios at the two larger power plants (Fort Frances and Thunder Bay) is investigated. For an increase in price by 10% (IRTB1) and 20% (IRTB2) at the Thunder Bay CHP plant, the per unit profit increases by 148% and 296% on average for all FMUs, respectively. The FMUs, which are located near the Thunder Bay CHP plant (Black-Sturgeon Forest, Dog River-Matawin Forest, English River, Lakehead Forest, Spruce River Forest, and Lake Nipigon Forest), make more per unit profit from an increase in price at the Thunder Bay CHP plant (ranging from 3.82 to 7.32 $ODt−1 in IRTB1 and from 10.32 to 13.82 $ODt−1 in IRTB2) because the transportation costs from these FMUs to this plant are low. Similarly, for an increase in price by 10% (IRFF1) and 20% (IRFF2) at the Fort Frances CHP plant, the average per unit profit of FMUs increases by 54% and 108% on average for all FMUs, respectively. Finally, if there is an increase in price of forest biomass at both the Thunder Bay and Fort Frances CHP plants by 10% (IRTF1) and 20% (IRTF2), the average per unit profit of FMUs increases by 202% and 404%, respectively. Therefore, an increase in price of forest biomass at a power plant greatly helps to increase the profits of close by FMUs but only marginally helps to increase the profits of faraway FMUs. It is, therefore, advisable to put smaller power-plants in a wider range of area in order to make biomass collection a profitable business for FMUs. This will also help in reducing the distribution cost of power-to-small communities and increasing the rural employment through bioenergy production systems in the region.

4. Conclusions

This study investigates an optimal biomass supply chain for four large-scale biomass-based power plants in NWO, which has been a priority research area recently due to greater emphasis put on green energy sources in Canada and around the globe. Two optimization models have been developed and applied to analyze the impacts of biomass competition, likely to be created due to increasing demand for biomass feedstock on part of the four power plants in this region, on cost and gross margin structures for biomass based power plants. Twenty-two model scenarios (14 for cost minimization and 8 for profit maximization model) are used to study the impacts of changes relevant to economic and technological parameters of biomass procurement for each power plant and FMUs. The results show that per unit procurement costs are directly proportional to the size of the power plants in all scenarios. The highest percent increase in per unit biomass procurement costs from the base scenario is found in the increased truck charge rate by 20% scenario, followed by the scenario where only forest harvest residue biomass is used by all power plants.

The results from increased truck charge rate and FHR only scenarios point to some of the important policy implications for biomass supply chain management in NWO. The policy should focus more on controlling transportation costs by investing in transport infrastructure and logistics and rethinking of current emphasis on collecting only FHR biomass for bioenergy purpose in order to reduce the biomass procurement costs. Further, harvesting standing underutilized tree species reduces per unit procurement cost. Results from other scenarios are also important and interesting to evaluate the impacts of changing economic and technological parameters on cost structures for designing the improved decision-making strategies for power plants.

On the other hand, the results from profit maximization model are more useful to the biomass suppliers (contractors), who are interested in making maximum profits in supplying forest biomass feedstock to power plants from given FMUs. They can maximize their profits from FMUs that are closer to the power plants. However, their profits significantly increase if the power plants offer higher prices. This is possible only if the suppliers meet the quality standards and lead time requirements of the buyers. The variations in costs and gross margin structures under various model scenarios are explained by location of depletion cells relative to power plants, availability of each type of biomass in depletion cells, biomass demands, and differential processing costs for two types of biomass. This modeling framework can be applied elsewhere to study the similar problem of biomass supply chain. The results of such modeling can help decision makers make improved decisions relating to biomass supply chains for bioenergy production.

Acknowledgments

The authors are grateful that the Ontario Ministry of Energy and Infrastructure through the Ontario Centres of Excellence, NSERC Strategic Network on VCO, Abitibi-Bowater Inc., and Thunder Bay Community Economic Development Commission funded this research. They are thankful to Lakehead University, Confederation College, FPInnovations FERIC Division, Lac Des Mille Lac First Nation, Buchanan Group, and Ontario Ministry of Natural Resources for further support.