Next Article in Journal
A Low-Cost Spore Trap Allows Collection and Real-Time PCR Quantification of Airborne Fusarium circinatum Spores
Next Article in Special Issue
Productivity and Oil Content in Relation to Jatropha Fruit Ripening under Tropical Dry-Forest Conditions
Previous Article in Journal
The Concept of Landscape Structure, Forest Continuum and Connectivity as a Support in Urban Forest Management and Landscape Planning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimizing the Location of Biomass Energy Facilities by Integrating Multi-Criteria Analysis (MCA) and Geographical Information Systems (GIS)

1
ARC Centre for Forest Value, Discipline of ICT, College of Science and Engineering, University of Tasmania, Hobart, TAS 7001, Australia
2
Forest Industries Research Centre, University of the Sunshine Coast, Locked Bag 4, Maroochydore DC, Queensland 4558, Australia
3
Private Forests Tasmania, 30 Patrick Street, Hobart, TAS 7000, Australia
*
Author to whom correspondence should be addressed.
Submission received: 4 September 2018 / Revised: 12 September 2018 / Accepted: 18 September 2018 / Published: 20 September 2018
(This article belongs to the Special Issue Supply Chain Optimization for Biomass and Biofuels)

Abstract

:
Internationally forest biomass is considered to be a valuable renewable energy feedstock. However, utilization of forest harvesting residues is challenging because they are highly varied, generally of low quality and usually widely distributed across timber harvesting sites. Factors related to the collection, processing and transport impose constraints on the economic viability of residue utilization operations and impact their supply from dispersed feedstock locations. To optimize decision-making about suitable locations for biomass energy plants intending to use forest residues, it is essential to factor in these supply chain considerations. This study conducted in Tasmania, Australia presents an investigation into the integration of Multi-criteria analysis (MCA) and Geographical Information systems (GIS) to identify optimal locations for prospective biomass power plants. The amount of forest harvesting biomass residues was estimated based on a non-industrial private native resource model in Tasmania (NIPNF). The integration of MCA and a GIS model, including a supply chain cost analysis, allowed the identification and analysis of optimal candidate locations that balanced economic, environmental, and social criteria within the biomass supply. The study results confirm that resource availability, land use and supply chain cost data can be integrated and mapped using GIS to facilitate the determination of different sustainable criteria weightings, and to ultimately generate optimal candidate locations for biomass energy plants. It is anticipated that this paper will make a contribution to current scientific knowledge by presenting innovative approaches for the sustainable utilization of forest harvest residues as a resource for the generation of bioenergy in Tasmania.

1. Introduction

As a response to increased fuel cost and environmental concerns, the use of renewable energy source is being highly considered as an alternative to fossil fuel [1]. Globally, forest biomass is considered one of the main renewable energy sources, and its utilization as an energy source has increased rapidly in the last decades [2]. Policies supported by the public and all political parties has resulted in financial support and subsidies by governments globally, especially in Europe [3]. In the case of the United States, forest residues are also becoming valuable materials such as biomass energy feedstocks and engineered wood products [4]. Wood residues from commercial timber harvesting have been utilized to make various timber products for decades. Nearly all particle boards in the United States are made from timber processing residues, and about one-third of all pallets produced annually are made from recycled wood [5].
Tasmania, Australia’s southernmost state, has a long history of forest utilization for timber production. In the past, the majority of timber harvest concentrated on native forests, but in the future, plantations will become the primarily source of timber products [6]. In contrast to the bioenergy trends observed globally, forest biomass and residue utilization has been insignificant in Tasmania (as in the rest of Australia), and has got little political or public support. Because of the lack of economic incentives, a significant amount of harvesting and processing residues end up being burnt in the forest or dumped in landfills [7]. The harvesting and processing of native forests produce a significant volume of timber residues that could be utilized by biomass plants in Tasmania for the production of electricity from renewable forest biomass resources. Currently, there are no biomass plants in Tasmania; however, a project has been proposed to build and operate an energy plant near the town of Huonville, in southern Tasmania, but it is still under review waiting for its approval by the Tasmanian Government [7].
Forest residues, including unmerchantable trees, small-diameter trees, tops, limbs and chunks produced from mechanical thinning and conventional saw-timber harvesting operations provide an opportunity to produce bioenergy and bio-based forest products. New technologies that are capable of converting forest residues into high quality and sustainable bioenergy and useful bio-based products are emerging [8]. To a large degree, the underutilization of forest residues can be attributed to the high cost of collecting and transporting these residues to end user markets and the low market prices paid for delivered forest residues. Several studies have investigated innovative forest biomass operations that effectively improve access to harvesting sites and economic efficiency [9,10,11]. However, the inherent inefficiency of transporting low-density and high moisture content (MC) biomass feedstock to biomass facility still remains a fundamental economic barrier to its increased utilization.
The location of a bioenergy plant is affected by a wide variety of factors and criteria. For instance, the location of biomass facility is highly influenced by the location of biomass feedstock [12]. An effective and efficient supply chain and optimized logistics system are vital for the success of a biomass-based energy industry. Since transportation costs impact highly total biomass fuel costs, site selection for new biomass facilities in the proximity of currently non-used available biomass resources are the preferred option when designing biomass supply chain networks [13]. In addition, biomass plant design should factor in environmental and socio factors to determine the number, location, and size of biomass facilities within the network of feedstock collection points, plants, and storage units [14]. This design gives form, structure, and shape to the entire supply chain and logistics system [15]. For this reason, the optimal locations of the biomass facilities play a significantly important role in the biomass energy supply chain [16].
A Geographical Informational System (GIS) is a tool broadly used to investigate the potential availability of biomass feedstocks [17], and to minimize transportation costs through logistics analysis and distance calculations [18]. A GIS network analysis, location–allocation analysis tool can simulate the site competitions for biomass resources [19]. The integrated GIS and Multi Criteria Analysis (MCA) is one of the preferred tools for the selection and location of biomass facilities and for assessing the relative importance of the economic, environmental, and social criteria affecting this site selection [20,21,22,23,24,25]. The Analytical Hierarchy Process (AHP), which is based on outranking techniques, is one of the most applied techniques in MCA [26]. The methodology includes a system to get an estimate of weighting to factors that impact decisions. The influencing impact of the factors for the selection and location of biomass energy facilities are not always the same, and in reflection of regional priority demands, it is necessary to assign them a different relative weight [27].
This study aims to investigate the optimal location of prospective biomass power plants in Tasmania, Australia using forest harvesting residue estimation. Integrated GIS and AHP models, along with a supply chain cost analysis, were developed to determine the best candidate locations of biomass plants. Three main criteria (economic, environmental, and social) and sub-criteria were established and factored in to determine the best biomass facility location. Additionally, various moisture contents of biomass feedstock were assumed to determine their impact on biomass transportation cost. It is expected that the results of this study make a contribution to the future Tasmania biomass energy and provide vital information to stakeholders so that they make good planning decisions about the utilization of forest residues for heat and electricity production in the region.

2. Materials and Methods

The schematic overall flowchart of the approaches used and implemented in the study is presented in Figure 1. Available land areas were investigated using a GIS restriction model. Integrated GIS and the multicriteria assessment method AHP were applied to assign a weight to each main criterion (economic, environmental, and social) and each sub criteria [28]. Spatial analysis required a large GIS dataset, both in raster and vector formats. In this study, most of the GIS data was provided by industrial partners (Private Forests Tasmania—PFT, Esk mapping company and Sustainable Timber Tasmania—STT) and gathered from different institutional websites, such as the Australian Government National Map. The optimal location of biomass facilities comprised three analytical phases including, land availability, land suitability, and location–allocation of biomass facilities. A selection of the available land for the biomass plant construction was initially investigated applying a restriction model (land availability); then, a number of suitable locations were identified applying a land suitability model according to a set of specific criteria; finally, selected candidate optimal locations were found by minimizing the total transportation distance of the feedstock location to each plant (location–allocation analysis). Detailed information of the methods applied is provided in the following sections.

2.1. Estimations of Biomass Availability in Tasmania

Potential forest harvesting and availability of residues have been estimated for public and private managed forests by STT and PFT, respectively. Modelling of public native forests and plantations have been conducted by STT as part of Tasmanian Forests Agreement (TFA) implementation process, while wood flow modelling for private native forests and plantations have been completed by PFT. The definition of two types of forest residues is described in Table 1. STT forest modelling methodology for public eucalypt native forests was described in [29,30,31]. In the private forest sector, the native forest resource is highly dispersed; according to PFT, the forest estate is distributed across 7500 land owners. Due to the difficulties in aligning land owner objectives over a broad ownership range, residue volumes were estimated using a sustainable yield model developed by PFT in 2012 [32]. The model assumes that the available native forest resource has been harvested and reforested repeatedly over a period of 360 years to allow several rotations of slower growing dry eucalypt forests to be included in the wood flow. In addition, softwood residue estimation was concentrated on private ownership. The model used to estimate softwood residue availability was not based on softwood industry projections and may therefore not reflect actual yields. The details and key assumptions associated with these forest residue estimations were described in the Residue Solutions Project 2014 [33]. As part of the forecast residue modelling, STT and PFT estimated the volume of processing residues potentially available based on sawlogs and peeler log volumes. The analysis assumes that 60% of native forest and hardwood plantation sawlogs, 50% of softwood plantation sawlogs, and 11% of native forest and hardwood plantation peeler logs are converted into processing residues. Additionally, to account for supply chain costs, areas with residues under 1000 m3 were not included in the biomass availability analysis. The potential Tasmanian harvest volumes were projected from 2014 to 2019.
This paper explores a range of biomass energy options, with modelling of a biomass combustion plant as the primary focus. As a result, non-woody and wet residues (from aquaculture processing, dairy processing, and waste water treatment) have not been considered. On-going research is examining these wet biomass residues, and this will be published in a subsequent research paper.

2.2. Land Availability Analysis Using a Restriction Model

The location of a biomass plant must satisfy a number of criteria and constraints such as geological and environmental reserve restrictions, which are imposed by several government regulations. Tasmanian forest ecosystems are largely protected within parks and reserves. This and the critical requirements outlined by the regional or national levels, limit the potential location of biomass energy plants. Given that energy plants are inexistent in Tasmania, constraints and restriction factors were established from a prefeasibility study conducted by Dorset [34] and from previous published articles [12,23,35,36,37,38,39,40,41,42]. Several constraints, both geological and environmental, were considered to select the potential location of biomass plants. Thus, protected and reserved areas, including threatened species areas, wildlife habitat areas, water pollution areas, and areas of interest for the communities were excluded for the location of biomass plants. The value related to spatial constraints in the dataset was arranged and converted into geo-referred layers. In the first step, buffer zones were created through a proximity analysis, using the GIS layers. Based on the constraint descriptions, buffer zones were created along the margins of some of these areas to delimit a minimum protection area. In this study, the available biomass plant locations were always established outside of buffer areas. A list of constraints for the selection of biomass plants is presented in Table 2. Each GIS data layer including specific information of constraint values was converted and reclassified into a binary raster data format (20 m × 20 m cell size). Reclassified binary raster data were assigned with a 0 or 1 value, in which cells assessed “0” were considered restricted area and cells assessed “1” were considered available area with respect to that factor. Using reclassified binary maps, a final combined binary map was generated by multiplying all the reclassified binary raster layers using the raster calculation tool in ArcGIS 10.6 (ArcGIS Desktop: Release 10. Redlands, CA: Environmental Systems Research Institute) [43]. The output of the combined map was called “Available map”.

2.3. Land Suitability Analysis with a Weighted Linear Model

A weighted overlay function was used in ArcGIS 10.6 (suitability analysis), which included a multi-criteria evaluation. Through the suitability analysis, additional geographical information and other layers can be considered by decision-makers as factors that can positively affect the biomass plant location decision based on detected areas already classified as “Available” within the geographical areas [12]. The impact of each of these additional criteria varies, and, therefore, this step is required to correctly assign them a different relative weight [27,28]. This was done by applying an AHP algorithm, and the procedure is presented in detail in [28].
In this study, eight criteria were considered in the suitability analysis. The selected criteria were classified according to three different groups with different objectives, called “Economic”, “Environmental”, and “Social”. In the Economic group, there were four sub criteria layers that had a high correlation with the economic influence of the biomass plants. It was assumed that biomass plants located closer to an industrial area would have a higher weighted value compared to further biomass plants. The distance from road and accessibility of the biomass plant also determined a high weighted value as short distances from existing roads can reduce biomass feedstock transportation costs substantially.
Regarding the second objective group (environmental), it was preferred to locate a power plant in a relatively flat area (slope less or equal to 15%) and in a middle elevation (elevation from 0 to 500 m) to avoid potential damage that could occur on steep sites and sites with sensitive geophysical conditions. In addition, to reduce water contamination risks, the locations of biomass plant were only possible in places located far away from main river and water bodies (interval from 0 to 500 m).
Finally, the criteria assigned to the social group included the local employment rate and the population density in the local government area of the identified biomass plant locations. To support local government economic impact, areas with a low and a middle employment rate were selected as preferred locations. In addition, middle population density areas (from 50 to 300 people/km2) were prioritized for the location of the biomass plants. Weights obtained from an expert consultation survey were assigned to each objective group (economic, environmental, and social) and then to each sub criterion. The sum of the weights in different hierarchical level was equal to one.
Once the weights were assigned to each criterion, a weighted reclassified raster format layer (20 m × 20 m cell size), and corresponding final suitability map, were created by a weighted linear combination of the raster layers. The applied equation was as follows (Equation (1)):
S i   =   p = 1 m w p   C i p ,
where S i is the suitability value of the ith cell in the final grid, w p is the weight allocated to the pth criterion from the AHP analysis, C i p is the ith cell value in the grid of the pth suitable criterion layer, and m is the total number of criteria in this analysis [19]. The values of all the criteria were standardized before including them in the equation. The calculated values Equation (1) in final raster were reclassified into seven level classes, where higher values represented more suitable places for the location of a biomass plant. After reclassified raster values were computed, cells that were assigned to levels 6 and 5 were identified as the most suitable locations for a potential biomass plant. A potential site was located in each cell by overlapping the areas under classes 6 and 5 within a grid mesh of 200 ha.

2.4. Biomass Energy Feedstock Supply Chain Costs

The Network Analysis tool in ArcGIS 10.6, which calculates transportation distances using a road network data set, rather than straight line or radius-based distances, was selected for the analysis. The location–allocation analysis was used to select potential bioenergy facility sites using the p-median problem solver included in the Network Analyst tool [44].
The aim of the location–allocation procedure (p-median problem solver) is to locate p candidate bioenergy facilities among n potential selected suitable bioenergy facility candidates (n > p) while satisfying a set of demands so that the total sum of weighted distances between each demand and facility location is minimized [45]. Additionally, the location–allocation procedure automatically selects the shortest travel or least cost distance so that the sum of the distance between the bioenergy facility candidates locations l and a set of biomass resource location points p is minimized (Equation (2)) [19]:
m i n l p   w p   d l p   x l p ,
where wp = weight associated with each biomass available point p, dlp = distance between available biomass point p and potential bioenergy facility l. xlp = 1 if biomass available point p is assigned to bioenergy facility l, and zero otherwise.
Using Dijkstra’s algorithm, the location–allocation solver generates an origin-destination matrix with the shortest-path costs between bioenergy plant candidates and available bioenergy resource points. It combines several techniques including a vertex substitution heuristic and a refining metaheuristic to achieve a near optimal solution [44]. Transportation costs (Equation (3)) were calculated using the travel distance along the optimal route derived from the location and allocation analysis. Equation (4) presents the equation to calculate the weighted transportation cost for bioenergy plant j ( WTC j ):
C i j   =   F C   +   2 V C k = 1 N d k ,
W T C j   =   ( i = 1 m B i   C i j ) / B t ,
where C i j is the total transportation cost ($ per dry tonne) for the optimal route between the bioenergy plant candidate j and available bioenergy resource location i; B i is the total biomass at source i; m is the total number of bioenergy resource locations assigned to plant j; B t is the total bioenergy resource assigned to plant j; FC is the fixed cost related to collecting the biomass energy feedstock ($ per dry tonne); VC is the variable cost related to travel distance ($ per tonne-km); d k and N are the length of the traveled road segment (km) and number of total segments along the optimal route between the plant j and the bioenergy resource location i, respectively [19].
This study also investigated the effect of MC on transport cost. For an MC of 55%, 50%, 45%, and <40%, transport costs (VC) were assumed to be $0.11/tonne/km, $0.10/tonne/km, $0.09/tonne/km, and $0.09/tonne/km, respectively [46]. This assumes a chip van truck with a volumetric capacity of 83.4 m3 and net payload of 29.2 tonnes, which transports chips from residues with a basic density 450 kg/m3. A value of $9.5 per tonne was assumed for FC, taking an in-field delimber, debarker, chipper machine (Peterson an Astec Inc., Eugene, OR, USA) as the basis for the calculation [47].
Several scenarios were simulated, assuming one to four potential biomass facilities. Through the location–allocation analysis, the locations of the most suitable bioenergy facilities were identified using available harvesting biomass sources in the Tasmania region. The transportation distances from every biomass availability location point to potential bioenergy facility locations were estimated with GIS software and using the actual road network date set. This calculation was applied for every candidate bioenergy facility location, ensuring that each biomass available location point was allocated to only one bioenergy facility.
The size of each facility is highly impacted by the number of bioenergy facility candidates and the extension of the associated biomass available supply zones. Greater and more dispersed feedstocks affect the size of the biomass facility and increase transportation distances and costs [48]. To determine the location and size of the biomass plants, the actual available biomass feedstock was re-estimated based on the results obtained from the location–allocation analysis. Subsequently, the selected supply points that resulted from the location–allocation analysis were extracted and clipped with the original biomass availability GIS data. Thus, the logistically available biomass feedstocks were calculated using the field calculator in ArcGIS 10.6.

3. Results and Discussion

3.1. Estimations of Biomass Availability in Tasmania

The results showed that a large base of biomass resources could be used as feedstock to supply biomass energy plants in Tasmania. The potential biomass availability is shown in Figure 2. The total forecasted biomass availability was 13,803,000 green metric tonnes per year (gmt/year) in the 5-year period from 2014 to 2019. Most of the biomass feedstock corresponded to residues from pulpwood and another stem wood harvesting operation. About 79% of the available biomass feedstock was produced from low quality logs not suitable for wood chipping to make pulp and paper. In addition, 95% of the total available biomass feedstocks were sourced from harvesting residue. Only 5% of total biomass availability was generated from processing residues. In the supply zone, Murchison East (ME) had the largest amount (3,125,000 gmt/year) of potential biomass availability in Tasmania. Bass North was the second largest biomass available zone with an availability of 2,499,000 gmt/year. Results obtained from the biomass availability analysis indicate that most of the available biomass feedstocks were produced from the northern side of Tasmania including Murchison West, East, Mersey, and Bass North supply zone. The distribution of biomass harvesting residue is shown in Figure 3.

3.2. Location of Bioenergy Plants

Figure 4 shows the availability map with the available areas for the location of bioenergy plants. About 49.6% of the Tasmanian land was accounted for as available area, while the non-available areas corresponded primarily to national conservation areas and areas prone to geological risk such as floods and natural disaster. White cells in Figure 4 represent unavailable areas where biomass energy plant cannot be located.
Following the AHP analysis, the first hierarchical level, main criteria, and the second one, sub-criteria, were weighted on the basis of a pairwise comparison from an expert surveying and through a matrix computation [12]. The results of the AHP analysis, including weights and corresponding consistency ratios, are presented in Table 3. The results indicated that the highest weight was given to the economic criteria. This was emphasized at the second hierarchical level (sub-criteria), where feedstock availability scored the highest weight compared to all the second level sub-criteria. In the social criteria (main criteria level), the population of the area was given the highest weight. Finally, environmental factors were considered as the third main criteria, with slope and flat areas located within a distance range between 100 m and 500 m from the water bodies as the preferred sites for the location of bioenergy power plants. The value of consistency ratio (CR) showed a high degree of consistency, confirming the adequacy of the weight values (Table 3). The results of the AHP analysis, and the corresponding weights, were included in the suitability model, where the reclassified weights were used to generate a final suitability map for the potential location of the bioenergy plants (Figure 4). Combined availability and suitability maps generated from raster calculations allowed the identification of the final possible biomass plant locations, with as many as 125 possible locations with a suitability value ≥ 5 being selected (Figure 4). Most of the highly suitable bioenergy plant candidates (suitability value = 6) were located in the Murchison East, Mersey, Bass North, Huon, and Derwent East supply zones. High suitability was attributed to sites with a large amount of biomass residue availability (Figure 2).
To secure enough feedstock for each final candidate location area, only selected locations larger than 10 ha were chosen as potential candidates for the location of the bioenergy plants. The candidates also required to be located near the road network and major roads to allow shorter transportation distances. The most suitable candidates were located close to the Tasmanian road network, indicating that the transport logistics criteria were satisfied. Figure 4 shows the 125 potential bioenergy facility locations that resulted from the suitability analysis.

3.3. Biomass Supply Chain Costs in Tasmania

Several scenarios were simulated with the location–allocation analysis to assess the impact of an increasing number of bioenergy power plants. From the location–allocation analysis, weighted unit transportation costs were calculated for an increasing number of bioenergy power plants. Figure 5 shows the four sites that were selected from the candidate locations (125 candidates) by solving the p-median problem using 15,322 biomass availability sources.
The estimated supply chain costs for all identified locations are presented in Table 4. Harvesting, processing, and loading costs were not included in the supply chain cost calculations, as the primary focus of the study was on transportation costs. As expected, the total transportation cost was highly influenced by the distance between the biomass facility and the biomass feedstock available locations. Scenario 1, in which one biomass facility within a 200-km radius was selected, has the highest coverage rate (83.93%) among the identified biomass facility locations. However, the total transportation cost was very high due to the long distance between the selected facility and the biomass feedstock locations. When increasing the number of facilities, the transportation distance and total cost decreased rapidly compared with a single and long-distance biomass facility scenario. Scenario 2, in which two biomass facilities within a 100-km radius were selected, had a coverage rate that was approximately 74% of total biomass available locations. In Scenario 2, the average round distance was about 50% (97.22 km) of that calculated in Scenario 1 (189.66 km). The total transportation cost decreased substantially by increasing the number of biomass facilities. The results shown in Table 4 highlight the importance of optimizing the location of potential biomass facilities to reduce transportation costs. Based on the truck transportation cost scenario, the optimal number of biomass plants was 3 in the Tasmania region. As seen in Figure 6, there is an indication that transport costs start increasing when four or more biomass plants are built.
Different MC levels (55%, 50%, 45%, and <40%) were used to simulate the impact of this parameter on biomass transportation costs. The assumed transport costs per tonne and kilometer were $0.11/tonne-km (MC = 55%), $0.10/tonne-km (MC = 50%), $0.09/tonne-km (MC = 45%), and $0.08/tonne-km (MC < 45%), respectively. These costs were assumed based on a chip van truck with a volumetric capacity of 83.4 m3 and 29.2 gmt of net payload that transports chips from residues with a basic density = 450 kg/m3. Figure 7 shows the impact of MC on biomass transport costs. A fall in MC from 55% to 40% resulted in lower total transportation costs in all the biomass facility location scenarios analyzed. This confirms the results from previous studies, which have identified MC as one of the most important factors affecting the calorific value and transportation costs in biomass supply chains [49,50,51,52,53,54,55]. Thus, MC and transport distance constitute the factors with the major impacts on the biomass transportation efficiency from harvesting sites to biomass facilities [46,56].

4. Conclusions

This study investigated the optimal location of potential biomass facilities in Tasmania, Australia, using a forest harvesting residue estimation procedure. In addition, an integrated multi-criteria analysis and Geographical Information System (GIS-AHP) model in conjunction with a supply chain cost analysis were developed to determine the best candidate locations. The estimated biomass availability indicated that Tasmania has an abundant biomass feedstock to operate biomass and energy facilities. After running several simulations on transportation cost scenarios, it was determined that three biomass plants within a radius of 80 km was the best option for Tasmania’s future biomass energy plan. Considerable efforts in biomass supply chain research and life cycle analyses are still required to widely investigate the efficiency and total operation costs for the Tasmanian conditions [57,58]. From the supply chain costs analysis conducted in this study, the distance between biomass facilities and biomass feedstock locations, in addition to the MC of biomass feedstock, were the factors with the greatest impacts on transport costs. Previous biomass transportation studies in Australia have identified impact factors that have proved to be the most challenging issues in the forest biomass supply chain [52,55,59]. In future studies, other biomass supply chain costs including harvesting, processing, loading, and storage should be included so that the impacts of the location and size of biomass facilities in Tasmania on overall supply chain costs are properly quantified.

Author Contributions

Conceptualization, H.W. and P.T.; Methodology, H.W.; Software, H.W.; Validation, M.A., P.T. and M.M.; Formal Analysis, H.W.; Investigation, H.W.; Resources, M.M. and M.A.; Data Curation, H.W. and M.M.; Writing-Original Draft Preparation, H.W.; Writing-Review & Editing, P.T., M.S.T., and M.A.; Visualization, H.W.; Supervision, P.T. and M.S.T.; Project Administration, P.T.; Funding Acquisition, P.T. Authorship must be limited to those who have contributed substantially to the work reported.

Funding

This research was funded by Australian Research Council, Industrial Transformation Training Centres scheme, grant number (IC150100004).

Acknowledgments

This research used in this paper is funded through a major Australian Research Council (ARC), Industrial Transformation Training Centres scheme, research project called: The ARC Centre for Forest Value <www.utas.edu.au/arc-forest-value>.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Solomon, S.; Qin, D.; Manning, M.; Averyt, K.; Marquis, M. Climate Change 2007-the Physical Science Basis: Working Group I Contribution to the Fourth Assessment Report of the IPCC; Cambridge University Press: Cambridge, UK, 2007; Volume 4. [Google Scholar]
  2. Ferranti, F. Energy Wood: A Challenge for European Forests Potentials, Environmental Implications, Policy Integration and Related Conflicts; EFI Technical Report 95; European Forest Institute: Joensuu, Finland, 2014; p. 157. [Google Scholar]
  3. Rothe, A. Forest Biomass for Energy: Current and Potential Use in Tasmania and a Comparison with European Experience; Zentrum Wald Forst Holz Weihenstephan, University of Applied Sciences, Weihenstephan-Triesdorf: Hobart, Australia, July 2013. [Google Scholar]
  4. Han, S.-K.; Han, H.-S.; Bisson, J.A. Effects of grate size on grinding productivity, fuel consumption, and particle size distribution. For. Prod. J. 2015, 65, 209–216. [Google Scholar] [CrossRef]
  5. McKeever, D.B. Inventories of woody residues and solid wood waste in the United States, 2002. In Proceedings of the Ninth International Conference, Inorganic-Bonded Composite Materials, Vancouver, BC, Canada, 10–13 October 2004; p. 2004. [Google Scholar]
  6. URS-report. Stage 1-Residue Options Identification and Analysis; URS Australia Pty Ltd.: Southbank, VIC, Australia, 2015. [Google Scholar]
  7. Rothe, A.; Moroni, M.; Neyland, M.; Wilnhammer, M. Current and potential use of forest biomass for energy in Tasmania. Biomass Bioenergy 2015, 80, 162–172. [Google Scholar] [CrossRef]
  8. Han, H.-S.; Jacobson, A.; Bilek, E.T.; Sessions, J. Waste to wisdom: Utilizing forest residues for the production of bioenergy and biobased products. Am. Soc. Agric. Biol. Eng. 2018, 5–10. [Google Scholar]
  9. Rawlings, C.; Rummer, B.; Seeley, C.; Thomas, C.; Morrison, D.; Han, H.; Cheff, L.; Atkins, D.; Graham, D.; Windell, K. A Study of How to Decrease the Costs of Collecting, Processing and Transporting Slash; Studies for Three Sites in Montana; Montana Development Council: Montana, MT, USA, 2004. [Google Scholar]
  10. Harrill, H.; Han, H.-S. Application of hook-lift trucks in centralized logging slash grinding operations. Biofuels 2010, 1, 399–408. [Google Scholar] [CrossRef]
  11. PI, D.W.P.; Harrod, R.J.; Dodson, E.K.; Lehmkuhl, J.F.; Ohlson, P. Fuel Succession, Post-Fire Logging, and Future Fire Behavior: Addressing the “Reburn Problem”; 06-3-4-16; Final Report to the Joint Fire Sciences Program; Wenatchee Forestry Sciences Laboratory: Wenatchee, WA, USA, May 2010. [Google Scholar]
  12. Delivand, M.K.; Cammerino, A.R.B.; Garofalo, P.; Monteleone, M. Optimal locations of bioenergy facilities, biomass spatial availability, logistics costs and ghg (greenhouse gas) emissions: A case study on electricity productions in south Italy. J. Clean. Prod. 2015, 99, 129–139. [Google Scholar] [CrossRef]
  13. Noon, C.E.; Daly, M.J. Gis-based biomass resource assessment with bravo. Biomass Bioenergy 1996, 10, 101–109. [Google Scholar] [CrossRef]
  14. Burak, A.; Harry, C.; David, W.; Kevin, G.; Sujith, S.; Mario, E.; Norman, S. Woody biomass and mill waste utilization opportunities in alabama: Transportation cost minimization, optimum facility location, economic feasibility, and impact. Environ. Prog. Sustain. Energy 2011, 30, 720–732. [Google Scholar]
  15. Bowersox, D.J.; Closs, D.J.; Cooper, M.B. Supply Chain Logistics Management; McGraw-Hill: New York, NY, USA, 2002; Volume 2. [Google Scholar]
  16. Gold, S.; Seuring, S. Supply chain and logistics issues of bio-energy production. J. Clean. Prod. 2011, 19, 32–42. [Google Scholar] [CrossRef]
  17. Beccali, M.; Columba, P.; D’Alberti, V.; Franzitta, V. Assessment of bioenergy potential in sicily: A gis-based support methodology. Biomass Bioenergy 2009, 33, 79–87. [Google Scholar] [CrossRef]
  18. Montgomery, T.D.; Han, H.-S.; Kizha, A.R. Modeling work plan logistics for centralized biomass recovery operations in mountainous terrain. Biomass Bioenergy 2016, 85, 262–270. [Google Scholar] [CrossRef]
  19. Sultana, A.; Kumar, A. Optimal siting and size of bioenergy facilities using geographic information system. Appl. Energy 2012, 94, 192–201. [Google Scholar] [CrossRef]
  20. Yue, C.-D.; Yang, G.G.-L. Decision support system for exploiting local renewable energy sources: A case study of the chigu area of southwestern Taiwan. Energy Policy 2007, 35, 383–394. [Google Scholar] [CrossRef]
  21. Buchholz, T.; Rametsteiner, E.; Volk, T.A.; Luzadis, V.A. Multi criteria analysis for bioenergy systems assessments. Energy Policy 2009, 37, 484–495. [Google Scholar] [CrossRef]
  22. Van Dael, M.; Van Passel, S.; Pelkmans, L.; Guisson, R.; Swinnen, G.; Schreurs, E. Determining potential locations for biomass valorization using a macro screening approach. Biomass Bioenergy 2012, 45, 175–186. [Google Scholar] [CrossRef]
  23. Perpiña, C.; Martínez-Llario, J.C.; Pérez-Navarro, Á. Multicriteria assessment in gis environments for siting biomass plants. Land Use Policy 2013, 31, 326–335. [Google Scholar] [CrossRef]
  24. Mourmouris, J.; Potolias, C. A multi-criteria methodology for energy planning and developing renewable energy sources at a regional level: A case study Thassos, Greece. Energy Policy 2013, 52, 522–530. [Google Scholar] [CrossRef]
  25. Boggia, A.; Cortina, C. Measuring sustainable development using a multi-criteria model: A. case study. J. Environ. Manag. 2010, 91, 2301–2306. [Google Scholar] [CrossRef] [PubMed]
  26. Wang, J.-J.; Jing, Y.-Y.; Zhang, C.-F.; Zhao, J.-H. Review on multi-criteria decision analysis aid in sustainable energy decision-making. Renew. Sustain. Energy Rev. 2009, 13, 2263–2278. [Google Scholar] [CrossRef]
  27. Massimo, M. Optimal Plant Size and Feedstock Supply Radius: Minimize the Production or Maximize the Profits? In Proceedings of the 20th European Biomass Conference and Exhibition, Milan, Italy, 18–22 June 2012; pp. 1–7. [Google Scholar]
  28. Saaty, T.L. Decision making with the analytic hierarchy process. Int. J. Serv. Sci. 2008, 1, 83–98. [Google Scholar] [CrossRef]
  29. Whiteley, S. Calculating the sustainable yield of tasmania’s state forests. TASFORESTS-HOBART- 1999, 11, 23–34. [Google Scholar]
  30. Forestry Tasmania. Sustainable High Quality Eucalypt Sawlog Supply from Tasmanian State Forest; Forestry Tasmania: Tasmanian, Australia, 2007. [Google Scholar]
  31. Burgman, M.; Robinson, A. Review of Tasmanian Forest Estate Wood Supply Scenarios; Final Report to the Independent Verification Group, Intergovernmental Agreement; University of Melbourne: Melbourne, Australia, March 2012. [Google Scholar]
  32. Wilson, J. Private Forest Resourcing Model 2012: Where Are the Best Locations for a Hypothetical Tasmanian Wood Processing Plant; Unpublished Report of Esk Mapping & GIS to Private Forests; Private Forests Tasmania: Tasmania, Australia, 2012. [Google Scholar]
  33. PFT. Residue Solution Project (Private Forests Tasmania); Private Forests Tasmania: Tasmania, Australia, 2014. [Google Scholar]
  34. Wilson, J. Dorset Woody Biomass Pre-Feasibility Study 2013: A Regional Inventory of Potential Woody Biomass Resources Surrounding Scottsdale; Private Forests Tasmania: Tasmania, Australia, October 2013. [Google Scholar]
  35. Viana, H.; Cohen, W.B.; Lopes, D.; Aranha, J. Assessment of forest biomass for use as energy. Gis-based analysis of geographical availability and locations of wood-fired power plants in portugal. Appl. Energy 2010, 87, 2551–2560. [Google Scholar] [CrossRef]
  36. Shi, X.; Elmore, A.; Li, X.; Gorence, N.J.; Jin, H.; Zhang, X.; Wang, F. Using spatial information technologies to select sites for biomass power plants: A case study in guangdong province, china. Biomass Bioenergy 2008, 32, 35–43. [Google Scholar] [CrossRef]
  37. Herrera-Seara, M.; Aznar Dols, F.; Zamorano, M.; Alameda-Hernández, E. Optimal location of a biomass power plant in the province of granada analyzed by multi-criteria evaluation using appropriate geographic information system according to the analytic hierarchy process. In Proceedings of the International Conference on Renewable Energies and Power Quality (ICREPQ’10), Granada, Spain, 23–25 March 2010; University of Granada: Granada, Spain, 2010. [Google Scholar]
  38. Zhang, F.; Johnson, D.M.; Sutherland, J.W. A gis-based method for identifying the optimal location for a facility to convert forest biomass to biofuel. Biomass Bioenergy 2011, 35, 3951–3961. [Google Scholar] [CrossRef]
  39. Perpina, C.; Alfonso, D.; Pérez-Navarro, A.; Penalvo, E.; Vargas, C.; Cárdenas, R. Methodology based on geographic information systems for biomass logistics and transport optimisation. Renew. Energy 2009, 34, 555–565. [Google Scholar] [CrossRef]
  40. Rodríguez, R.; Gauthier-Maradei, P.; Escalante, H. Fuzzy spatial decision tool to rank suitable sites for allocation of bioenergy plants based on crop residue. Biomass Bioenergy 2017, 100, 17–30. [Google Scholar] [CrossRef]
  41. Comber, A.; Dickie, J.; Jarvis, C.; Phillips, M.; Tansey, K. Locating bioenergy facilities using a modified gis-based location–allocation-algorithm: Considering the spatial distribution of resource supply. Appl. Energy 2015, 154, 309–316. [Google Scholar] [CrossRef]
  42. Zubaryeva, A.; Zaccarelli, N.; del Giudice, C.; Zurlini, G. Spatially explicit assessment of local biomass availability for distributed biogas production via anaerobic co-digestion–mediterranean case study. Renew. Energy 2012, 39, 261–270. [Google Scholar] [CrossRef]
  43. ESRI. Advanced GIS Spatial Analysis Using Raster and Vector Data; Environmental System Research Institute: New York, NY, USA, 2001. [Google Scholar]
  44. ESRI. Algorithm Used by Network Analyst; ESRI: Redlands, CA, USA, 2011. [Google Scholar]
  45. Daskin, M.S. Network and Discrete Location: Models, Algorithms, and Applications; John Wiley & Sons: Queens Land, Australia, February 2011. [Google Scholar]
  46. Ghaffariyan, M.R.; Acuna, M.; Brown, M. Analysing the effect of five operational factors on forest residue supply chain costs: A case study in western australia. Biomass Bioenergy 2013, 59, 486–493. [Google Scholar] [CrossRef]
  47. Ghaffariyan, M.R.; Sessions, J. Comparing the Efficiency of Four Harvesting Methods in a Blue Gum Plantation in South-West Western Australia; CRC for Forestry: Tasmania, Australia, 2012. [Google Scholar]
  48. Monteleone, M.; Cammerino, A.; Karaoglanoglou, L.; Koullas, D.; Koukios, E. Agro-energy chains and potential from dedicated crops, agricultural and agro-industrial residues: The comparative study of the two south european regions in a 2020 perspective. In Proceedings of the 20th European Biomass Conference and Exhibition, Millan, Italy, 18–22 June 2012; pp. 371–376. [Google Scholar]
  49. Röser, D.; Mola-Yudego, B.; Sikanen, L.; Prinz, R.; Gritten, D.; Emer, B.; Väätäinen, K.; Erkkilä, A. Natural drying treatments during seasonal storage of wood for bioenergy in different European locations. Biomass Bioenergy 2011, 35, 4238–4247. [Google Scholar] [CrossRef]
  50. Ghaffariyan, M.R.; Acuna, M.; Brown, M. Natural drying and optimising a forest residue supply chain to reduce the total operating costs: A case study in western Australia. Ind. Bull.-Aust. For. Oper. Res. Alliance (AFORA) 2014, 8, 1–5. [Google Scholar]
  51. Murphy, G.; Kent, T.; Kofman, P.D. Modeling air drying of sitka spruce (picea sitchensis) biomass in off-forest storage yards in Ireland. For. Prod. J. 2012, 62, 443–449. [Google Scholar] [CrossRef]
  52. Acuna, M. Timber and biomass transport optimization: A review of planning issues, solution techniques and decision support tools. Croat. J. For. Eng. J. Theor. Appl. For. Eng. 2017, 38, 279–290. [Google Scholar]
  53. Sosa, A.; Acuna, M.; McDonnell, K.; Devlin, G. Controlling moisture content and truck configurations to model and optimise biomass supply chain logistics in Ireland. Appl. Energy 2015, 137, 338–351. [Google Scholar] [CrossRef] [Green Version]
  54. Sosa, A.; Acuna, M.; McDonnell, K.; Devlin, G. Managing the moisture content of wood biomass for the optimisation of ireland’s transport supply strategy to bioenergy markets and competing industries. Energy 2015, 86, 354–368. [Google Scholar] [CrossRef]
  55. Acuna, M.; Anttila, P.; Sikanen, L.; Prinz, R.; Asikainen, A. Predicting and controlling moisture content to optimise forest biomass logistics. Croat. J. For. Eng. J. Theor. Appl. For. Eng. 2012, 33, 225–238. [Google Scholar]
  56. Möglichkeiten und Probleme Beim Holztransport. Available online: https://zidapps.boku.ac.at/abstracts/oe_list.php?paID=3&paSID=5784&paSF=&paCF=0&paLIST=02006: na (accessed on 10 September 2018).
  57. Cherubini, F.; Bird, N.D.; Cowie, A.; Jungmeier, G.; Schlamadinger, B.; Woess-Gallasch, S. Energy-and greenhouse gas-based LCA of biofuel and bioenergy systems: Key issues, ranges and recommendations. Resour. Conserv. Recyc. 2009, 53, 434–447. [Google Scholar] [CrossRef]
  58. Ghaffariyan, M.R. Review of European biomass harvesting technologies. Silva Balc. 2010, 11, 5–20. [Google Scholar]
  59. Acuna, M.; Mirowski, L.; Ghaffariyan, M.R.; Brown, M. Optimising transport efficiency and costs in Australian wood chipping operations. Biomass Bioenergy 2012, 46, 291–300. [Google Scholar] [CrossRef]
Figure 1. The overall approach to determine optimal locations of biomass facilities.
Figure 1. The overall approach to determine optimal locations of biomass facilities.
Forests 09 00585 g001
Figure 2. The forecast Tasmania state-wide residue flow estimates annually in thousand green metric tonnes per year in the 2014 to 2019 period.
Figure 2. The forecast Tasmania state-wide residue flow estimates annually in thousand green metric tonnes per year in the 2014 to 2019 period.
Forests 09 00585 g002
Figure 3. Suitability map with the final potential bioenergy plant locations in Tasmania.
Figure 3. Suitability map with the final potential bioenergy plant locations in Tasmania.
Forests 09 00585 g003
Figure 4. Map of restricted and available areas, and suitability map for bioenergy plant locations in Tasmania.
Figure 4. Map of restricted and available areas, and suitability map for bioenergy plant locations in Tasmania.
Forests 09 00585 g004
Figure 5. Optimal number and location of potential bioenergy plants. Lines show connections between origin and destinations; Tasmania road network was considered for the analysis; four scenarios were evaluated: (a) Scenario 1, 200 km radius selected one biomass facility (location code: Select 1), (b) Scenario 2, 100 km radius selected two biomass facilities (location code: S2_A_91 and S2_B_100), (c) Scenario 3, 80 km radius selected three biomass facilities (location code: S3_A_90, S3_B_92, and S3_C_100), and (d) Scenario 4, 80 km radius selected four biomass facilities (location code: S4_A_15, S4_B_90, S4_C_92, and S4_D_100).
Figure 5. Optimal number and location of potential bioenergy plants. Lines show connections between origin and destinations; Tasmania road network was considered for the analysis; four scenarios were evaluated: (a) Scenario 1, 200 km radius selected one biomass facility (location code: Select 1), (b) Scenario 2, 100 km radius selected two biomass facilities (location code: S2_A_91 and S2_B_100), (c) Scenario 3, 80 km radius selected three biomass facilities (location code: S3_A_90, S3_B_92, and S3_C_100), and (d) Scenario 4, 80 km radius selected four biomass facilities (location code: S4_A_15, S4_B_90, S4_C_92, and S4_D_100).
Forests 09 00585 g005
Figure 6. Variation of total transportation cost with number of plants.
Figure 6. Variation of total transportation cost with number of plants.
Forests 09 00585 g006
Figure 7. Impact of moisture content (MC) on transport costs (for an explanation of the facility location codes, see Figure 5).
Figure 7. Impact of moisture content (MC) on transport costs (for an explanation of the facility location codes, see Figure 5).
Forests 09 00585 g007
Table 1. Definition of forest biomass availability.
Table 1. Definition of forest biomass availability.
Biomass TypeResidue SourceDefinition
Harvesting residuesPulpwoodLogs of any size that are not suitable for solid wood processing to a small-end-diameter of 8 cm
Other stem wood materialNon-merchantable stem wood that is usually left in the forest. It also includes stumps but excludes limbs, foliage and roots
Processing residuesSawlogs or peeler logOffcuts, woodchips and sawdust from timber processing
Table 2. List of the constraints to build the restriction model developed in this study.
Table 2. List of the constraints to build the restriction model developed in this study.
ConstraintBuffer (Meters)
Game reserve-*
Historic site-*
National park-*
Nature reserve-*
State reserve-*
Wellington park-*
Public authority land within world heritage area (WHA)-*
Conservation area-*
Nature recreation area-*
Regional reserve-*
Informal reserve on sustainable timber Tasmania managed land-*
Future potential production forest-*
Informal reserve on other public land-*
Private land within WHA-*
Private sanctuary-*
Private sanctuary and conservation covenant-*
Private nature reserve and conservation covenant-*
Conservation covenant perpetual -*
Indigenous protected area-*
Conservation covenant-*
Management agreement-*
Management agreement and stewardship agreement-*
Stewardship agreement-*
Other private reserve-*
Tasmania point of interest-*
Tasmania threatened species-*
Lake bordering areaAround lake bordering area 300 m buffer zones
River, streams and waterwaysAround main river, river, streams area 150 m buffer zones
-* Constraint values were converted into a binary raster format (0 = restricted area, 1 = available area).
Table 3. Weights and preference factors in the Analytical hierarchy process (AHP) procedure.
Table 3. Weights and preference factors in the Analytical hierarchy process (AHP) procedure.
Main CriteriaWeights (Criteria)Second Level: Sub-CriteriaWeights (Sub-Criteria)Total WeightConsistency Ratio
Economic0.540Feedstock availability0.7470.4030.04
Industrial area0.1060.057
Transport logistics0.1470.079
Environmental0.163Elevation0.2400.0390.056
Slope0.4000.065
Water bodies0.3600.059
Social0.297Local employment rate0.4000.1190.061
Population0.6000.178
Sum of the weights1.0 3.01.0
Table 4. Estimated transportation cost of biomass feedstock.
Table 4. Estimated transportation cost of biomass feedstock.
Facility Location aRange (km) bDemand (No.) cCoverage Rate (%) dTotal AreaBiomass Availability (Green Metric Tons) eAve Distance (km) f$/tonne/km gTotal Cost ($) h
S120012,85983.93463,54031,757,500189.960.11663,586,379
S2_A_91100538035.11188,67413,375,50097.220.11143,040,272
S2_B_100100597338.98222,56814,694,00090.140.11145,696,888
S3_A_9080177211.5766,9585,148,00061.820.1135,007,430
S3_B_9280368024.02120,8958,402,00080.840.1174,713,945
S3_C_10080464430.31183,70311,636,00079.100.11101,244,836
S4_A_1580166410.8628,6642,743,50064.900.1119,585,847
S4_B_9080177211.5766,9585,148,00061.820.1135,007,430
S4_C_9280345422.54115,9507,994,00078.500.1169,028,190
S4_D_10080464430.31183,70311,636,00079.100.11101,244,836
a selected number of biomass facilities and identified location ID (location information seen above Figure 5). b coverage radius distance from selected facilities. c total number of biomass availability locations to selected biomass facility. d coverage rate of selected biomass facilities (selected facility covered biomass demand location/total number of biomass availability locations) and sum of coverage rate in selected facility number is total coverage rate in each scenario. e sum of total biomass availability in allocated biomass demand locations. f average round distance between identified facility locations and biomass demand locations. g truck transportation cost assumed that MC (moisture content) = 55%, the chip van truck has a volume capacity of 83.4 m3, a net payload of 29.2 tonnes, and transport chips with a basic density = 450 kg/m3. h The cost estimation was based on average distance between identified biomass facilities and demand biomass availability locations.

Share and Cite

MDPI and ACS Style

Woo, H.; Acuna, M.; Moroni, M.; Taskhiri, M.S.; Turner, P. Optimizing the Location of Biomass Energy Facilities by Integrating Multi-Criteria Analysis (MCA) and Geographical Information Systems (GIS). Forests 2018, 9, 585. https://0-doi-org.brum.beds.ac.uk/10.3390/f9100585

AMA Style

Woo H, Acuna M, Moroni M, Taskhiri MS, Turner P. Optimizing the Location of Biomass Energy Facilities by Integrating Multi-Criteria Analysis (MCA) and Geographical Information Systems (GIS). Forests. 2018; 9(10):585. https://0-doi-org.brum.beds.ac.uk/10.3390/f9100585

Chicago/Turabian Style

Woo, Heesung, Mauricio Acuna, Martin Moroni, Mohammad Sadegh Taskhiri, and Paul Turner. 2018. "Optimizing the Location of Biomass Energy Facilities by Integrating Multi-Criteria Analysis (MCA) and Geographical Information Systems (GIS)" Forests 9, no. 10: 585. https://0-doi-org.brum.beds.ac.uk/10.3390/f9100585

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