1. Introduction
Energy is unquestionably one of the most significant driving forces for developing and advancing a nation. With the increasing population and advancement of technology, energy consumption and demand are increasing rapidly. Alternative sources of energy are needed to achieve a better and sustainable system. Biodiesel is a cleaner-burning fuel consisting of alkyl esters produced from the transesterification of triglycerides (TGs) or esterification of free fatty acids (FFAs) with low-molecular-weight alcohols [
1]. Biodiesel has been one of the desired alternative fuels to replace petrodiesel due to its similar properties and lower emission of pollutants [
2]. The main sources for biodiesel production are edible and non-edible oil crops, such as palm oil, soybean, canola, rapeseed and microalgae [
3,
4,
5]. The conversion of edible oil crops into biofuel has created many issues around the ideology of using food as an energy source. In addition, concerns were raised in regards to the sustainability of the oil crop plantation. For example, a report revealed that the palm oil industry causes deforestation, habitat loss, water and air pollutions and forest fire issues [
6]. This raises concerns about the sustainability performance of biodiesel, although it is produced from renewable resources.
A variety of research has been conducted to improve the sustainability performance of biodiesel production. For example, an input–output model was proposed to improve the sustainability of oil palm plantations [
7]. The biodiesel supply chain is simulated in a multi-objective model to analyse the impact of market risk attitudes towards sustainability performance [
8]. A study on an animal-fat-based biodiesel supply chain network integrated with risk management and uncertainty consideration showed that the logistics cost was the second-highest expenditure next to installation cost [
9]. The integration of wastewater sludge has shown great potential in biodiesel production based on a data-driven optimisation model with the inclusion of uncertainty analysis and fuzzy approach [
10]. Biodiesel production derived from Jatropha Curcas is designed in a two-stage stochastic programming approach which shows that the integration of different kinds of biodiesels can result in a more comprehensive and green system [
11]. A study shows that microalgae-based biodiesel is feasible; however, the conversion substantially depends on the lipidic characteristic of the feed [
12]. Various studies show that waste cooking oil could be a promising source for biodiesel production for the more sustainable practice of utilising waste as an energy source [
13,
14,
15]. A comparison between petrodiesel, palm biodiesel and opium poppy biodiesel shows that the latter biodiesel has better energetically sustainable performance in the engine [
16].
From the literature, multiple alternative resources have been investigated to improve sustainability performance in biodiesel production. The majority of the researchers investigated the integration of alternative non-edible and renewable oil sources and compared the sustainability performance with the conventional approach. However, a recent study shows that different types of oil (e.g., palm oil, rapeseed oil, soybean oil, etc.) have different sustainability index profiles due to the variation in plantation practice and generation yield [
17]. For instance, in the comparison of sustainability performance between palm oil, soybean oil, rapeseed oil and sunflower oil, the study suggests that the production of palm oil has the highest sustainability index profile in land usage due to the relatively high oil generation yield, but performs poorly in deforestation and carbon footprint, while soybean oil production requires the most land, but has the least impact from fertiliser usage. There is no single renewable source that is able to perform well in all the sustainability index profiles. Nonetheless, utilisation of multiple oil crops as process feed has shown significant improvement in the overall sustainability index profile [
17]. Additionally, a detailed review of various direct transesterification processes has shown that the majority of the existing technologies are capable of handling a diverse feedstock type [
18]. In view of that, there is a research gap to investigate a biodiesel supply chain network that incorporates diverse feedstocks to create a balanced and optimum feedstock sustainability profile. This is critical to address the concern of sustainability issues from high dependence on specific oil crops in biodiesel production, such as extensive deforestation activity from the palm oil plantation or high land usage from soybean plantation. Additionally, diversifying the feedstock could improve the supply chain security issue [
19].
This paper aims to address the issue of non-optimum sustainability performance based on the usage of oil crops for biodiesel production by exploring the possibility of incorporating multiple oil crops as the process feed to create a balanced feed sustainability index profile. As mentioned above, the majority of the sustainability concern of oil crop utilisation of biodiesel is coming from the upstream processes like plantation and milling; this paper will only focus on the sustainability performance of different oil crops’ upstream processes. A mathematical optimisation model is proposed to evaluate the feasibility of incorporating multiple oil crops as process feed for biodiesel production to achieve a better feed sustainability index profile. In order to avoid significant process modification and redesign that will involve high upfront investment, the feedstock selection is based on the limitation of the existing biodiesel production process to handle feed fluctuations. A previous work on biomass supply chain optimisation has tackled a similar problem by introducing a diverse biomass feed selection approach based on existing process technology configurations [
20]. The work proposed a biomass selection model based on the biomass element characteristics instead of biomass species. In terms of biodiesel production from oil crops, studies have shown that the efficiency and performance of the process are subjected to the type of feed and the feed properties such as free fatty acid, iodine value and oxidation stability [
21,
22,
23]. Therefore, the same concept is adopted in this study to investigate the feasibility of integrating multiple oil crops into existing biodiesel processing plants based on the oil crop properties, as well as optimise the overall sustainability index of the feedstock mixture.
2. Model Formulation
Superstructure diagrams are used to illustrate the proposed concept of the oil crop supply chain optimisation model for biodiesel production in this work.
Figure 1 shows the conventional supply chain of biodiesel production from oil crops. In the conventional approach, the majority of the biodiesel refinery utilises a single source of oil crop as the process feed. This creates the concern of unsustainable sources of certain oil crops based on the oil crops’ plantation and refinery processes. For instance,
Figure 2 shows the comparison of the sustainability index profile of several major oil crops obtained from previous studies to optimise the oil crop usage to achieve the highest sustainability level [
17]. A total of five sustainability aspects were considered, including deforestation, oil yield, fertiliser impact, carbon footprint, and water impact [
24,
25,
26]. The study compared four main oil crops in the global market, which included palm oil, rapeseed oil, soybean oil and sunflower oil [
27,
28,
29,
30,
31]. A comparative sustainability index was used to highlight the differences between the sustainability level of each oil crop in the respective sustainability aspect, where 100% indicates the highest sustainability level and 0% indicates the lowest sustainability index. For example, oil palm has the highest oil yield at 610 gal of oil/acre, followed by rapeseed oil, sunflower oil and soybean oil at 122 gal of oil/acre, 98 gal of oil/acre and 46 gal of oil/acre, respectively. In the comparative sustainability index, oil palm is indexed at 100% due to the highest performance, and soybean oil at 0% due to the lowest oil yield among the oil crops, while rapeseed oil and sunflower oil are indexed at 14% and 9%, respectively, based on the linear interpolation from the maximum and minimum value [
17].
Figure 2 demonstrates the changes in the sustainability index profile between single and multiple oil crops utilisation. For instance, if biodiesel production is solely depending on palm oil, the system will have a significant negative sustainability impact on deforestation and carbon footprint, while production from rapeseed oil only will create high fertiliser impact and land usage due to low oil yield. Alternatively, if the biodiesel production utilised a combination of both oil crops, the produced biodiesel could have a more well-balanced sustainability index profile which subsequently eliminates the possible sustainability problems. The sustainability index profile of the mixed feed is determined based on the weight distribution of the oil crop mixture and their respective sustainability index. Note that in the comparative sustainability index approach, a 0% index does not imply that it is not acceptable, instead, it represents the performance is the lowest among the peers. For instance, soybean oil is still considered a potential resource for biodiesel production even though it has a score of 0% in the oil yield index. The comparative sustainability index is used to understand the position of the sustainability level of the process feed among the competing oil crops and the bottleneck of the sustainability profile of each considered oil crop. This can be used to motivate the effort to improve the sustainable practice in all the oil crop industries, as a new breakthrough in any one oil crop’s sustainability performance may affect the sustainability index of the rest of the competitive oil crops. Ideally, there is no upper limit of the best sustainability achievement as the requirement is constantly evolving and demanding as the technologies improve.
Figure 3 shows the proposed model in this paper. The model introduces a novel oil crop selection approach based on bio-oil properties. Instead of limiting specific oil crops as the feedstock of each refinery, the proposed model provides the flexibility of integrating multiple types of oil crops into the existing biodiesel refineries. In order to avoid significant fluctuations in the operation and performance of the biodiesel refineries, only a specific range of bio-oil properties will be accepted in each refinery. This is defined as the properties acceptance range in the context of this paper. The proposed properties acceptance range will be constructed based on the feed properties that would have significant impacts on the production yield and biodiesel quality. For example, the free fatty acid content is considered in this study as part of the feed selection criteria due to the direct influence on the performance of biodiesel production [
32,
33]. Apart from that, the stability of the oil crop and the produced biodiesel is another critical factor in ensuring consistency of the oil and biodiesel quality within the supply chain network. For instance, a study shows that crude palm oil has a shelf life of approximately six months [
34]. Iodine value is used to determine total unsaturation within the fatty acid which influences the oxidation stability of the oil crop and the produced biodiesel [
35,
36,
37]. Estimation of those properties in an oil crop mixture has been proven to correspond to the mixing ratio of the oil crops used in the biodiesel production [
38]. These properties would be included as part of the feed selection criteria in the proposed supply chain model to ensure adequate stability of the oil crop and produced biodiesel. Incorporating the feed selection criteria based on the oil crop properties increases the flexibility of the supply chain model to determine the best combination of oil crop utilisation to achieve a feasible and optimal system.
The following describes the model formulation proposed in this study. The model considers (i) mass balance of the supply chain, (ii) oil crop selection based on properties acceptance range, (iii) sustainability profile evaluation and (iv) cost consideration. Equations (1)–(5) describe the mass balance of the material logistics in the supply chain system. Equation (1) indicates that the total transported oil crops,
should be less than or equal to the availability of the oil crops,
, from each resource location,
f. In this context, each resource location is restricted to a single type of oil crop availability. Equation (2) determines the total oil crops received as process feed,
, in the respective biodiesel refineries,
r. The production rate of biodiesel is governed by Equation (3), where the produced biodiesel,
, in each biodiesel refineries,
r, is calculated based on a fixed biodiesel conversion yield,
, based on the total oil crop received as the process feed,
. Similar to Equation (1), Equation (4) limits the transportation of produced biodiesel,
, to be not more than the produced biodiesel in the respective biodiesel refineries. Lastly in material balance, Equation (5) indicates that all biodiesel demand,
, is required to be fulfilled in the supply chain network.
The second part of the model is to include the process feed selection based on oil crop properties. The respective oil crop properties,
p, from each resource,
f, is presented as
and the selection of the feed is bounded by the property acceptance range in respective biodiesel refineries. This acceptance range is proposed based on the natural oil crop properties’ fluctuation in a typical supply chain system and the original oil crop species used in the refinery. For example, the fluctuation of palm oil properties would be used as the basis to form the acceptance range for a palm-oil-based biodiesel refinery. Equations (6) and (7) limit the feed selection based on the upper and lower properties range,
and
, respectively.
Apart from the oil crop properties integration, the following equation represents the sustainability consideration of the model. The focus of this paper is to evaluate the possibility of producing a balanced sustainability index profile in each biodiesel refinery by utilising diverse oil crop feed (with each unique sustainability index profile) in order to tackle the sustainable concern of over usage of a specific oil crop. Therefore, the proposed formulation focusses on the sustainability index profile contributions from feed selection, in other words, the supply chain network between each resource and the biodiesel refinery,
.
represents the sustainability index,
s, of the respective oil crop from each resource,
f. Equation (8) governs that the combined sustainability profile of the feed used to produce biodiesel should be higher than a specific sustainability index limit,
. Depending on the regional legislation, the sustainability limit can be modified to ensure the proposed supply chain solution achieves the sustainability standard. Note that the sustainability index proposed in this model is associated with the type and origin of the oil crops. The index only considers the upstream processes of oil crops such as plantation and milling processes. Sustainability contributions from logistics and refineries are not included as part of the study in order to highlight the sustainability impact of biodiesel production from the feed selection.
Lastly, the model also considers the economic aspect of the supply chain. Equations (9)–(11) show the cost calculation for oil cost,
; production cost,
; and transportation cost,
, where
refers to the unit cost for each oil crop;
denotes the production cost per unit of biodiesel produced in each refinery; while
and
refer to the unit transportation cost of oil crop and biodiesel, respectively. The logistics cost in Equation (11) considered the distance between the source and destination and the transportation mode, such as land and sea travels. Equation (12) shows the total cost considered in the model.
Two objective functions (different scenarios) are proposed to be investigated based on the proposed model. First, the biodiesel supply chain problem can be optimised by minimising the
(Equation (12)) to determine the sustainability index profile of each refinery feed at the lowest cost. Apart from that, the model can be used to optimise the overall sustainability score (Equation (14)) of the system. Equations (13) and (14) show the calculation of the respective sustainability score of the refinery feed,
, and the overall sustainability score,
S_Overall_Score, of the system. Sustainability factors,
, are introduced to indicate the weightage distribution of the importance of sustainability aspects,
s, in the respective refinery,
r. Each refinery can have its own weightage distribution based on the regional sustainability requirement and policy. For instance, a weightage of 0.2 can be assigned to each sustainability aspect (total of five) shown in
Figure 2 if all aspects are equally important. Alternatively, higher weightage can be assigned for specific sustainability aspects based on the development direction or legislation requirement to give priority on a specific sustainability aspect.
Since the sustainability score is calculated based on the multiplication of the number of oil crops and their respective sustainability index, the range of sustainability scores presented in Equations (13) and (14) will not be the same as the sustainability index, which is in the range of 0% to 100%. In order to provide a better perspective and comparison of the improved feed sustainability, the optimised sustainability score is converted into a sustainability index format by dividing the score with the total amount of feed used. Equation (15) shows the sustainability index profile for each refinery,
, and Equation (16) shows the overall sustainability index of the system,
. These calculations were conducted outside of the optimisation model to simplify the model from a non-linear problem to a linear problem.