3.1. Physical–Chemical Characterization of Biomasses
All biomasses have a moisture content of less than 12% (
Table 3), which is considered the maximum moisture limit in the production of class A pellets described in ISO17225-6:2014 [
28]. Cork powder was the sample with the lowest moisture content (2.11%). Regarding volatile matter, biochar had the lowest content (34.10%), as expected due to the result of the pyrolysis process, and sawdust had the higher volatile matter content with 76.02%. Swine waste presented the highest ash content (14.10%) of all studied biomasses; this value hampers the production of standardized pellets solely with this biomass [
28].
Touriga Franca, Malvasia Fina, Tinta Barroca, Tinta Roriz, Manteúdo Branco, and Viosinho are all grape varieties that have similar properties, with moderate moisture content and volatile matter but high ash content. These samples had relatively low nitrogen, and no Sulfur content was detected (
Table 4), making them a good option for pellet production. On the other hand, cork powder, grape seeds and skins, and sawdust had the lowest levels of ash (0.54%, 2.90% and 0.63% respectively). The ash content of the stalks was exceedingly high when compared to sawdust or cork powder, both of which are arboreous biomasses like the stalks. This may be have been due to the presence of elements such as K, Na, Ca, and Mg [
35].
Biochar had the lowest hydrogen content of any biomass (1.45%) since a significant amount of hydrogen was removed during the pyrolysis process.
Because swine waste and biochar contained more than 2% nitrogen, it was impossible to produce a standardized pellet from these materials alone [
28]. All other biomasses had a nitrogen content of less than 1%, the maximum value for class A pellet classification, with cork powder having the lowest nitrogen content (0.61%).
Swine waste had a sulfur content of 0.13%, which is below than what is required for a class A pellet. A small amount of sulfur was detected in the seeds and skins samples.
Swine waste presented the lowest LHV of all biomasses studied (14.90 MJ/kg). On the other hand, biochar presented the highest LHV of all biomasses studied (33.35 MJ/kg). The LHV of cork powder and seeds and skin samples was greater than 20 MJ/kg. Grape stalks presented an LHV of around 15.5 MJ/kg, with Touriga Franca being the grape cultivar with the highest LHV (15.85 MJ/kg). Although swine waste had the lowest LHV, it was higher than the ISO standard requirement for non-woody pellets (ISO17225-6:2014).
These results are consistent with the values found in literature. In one study, swine waste was C (%) 35.59, N (%) 1.72, H (%) 4.84, and S (%) 0.44 [
36]. Small differences in the results obtained with those described in the literature can be explained by animal species and diet, bedding, storage, and handling of the manure [
37]. In another study, an unprocessed grape stalks and pomace substrate was subjected to a proximate and elemental analysis [
35]. The results of this study’s elemental analysis are similar to the ones found on the aforementioned study. The differences found on the proximate analyses could be explained by the higher moisture content measured by the authors (20.2%).
3.2. Using the Combination of Biomasses for Improving Swine Waste Use as Fuel
While maintaining a swine waste base, mixtures of biomasses with various stalk contents were studied. To take advantage of cork powder’s high energy value and low ash content, a mixture containing 20% cork powder was considered. This challenge is likely due to cork powder’s unique properties, including its high suberin and wax content [
38,
39], which hinder particle adhesion, and its low moisture content, which reduces the plasticity and compressibility needed for pellet formation. At higher cork powder concentrations, these factors further disrupted the biomass matrix, making pelletization impossible.
The mixture of two materials did not show a significant improvement in ash content when compared to swine waste, and in some cases, there was an increase in ash content, such as in sample 32, which has the highest ash content verified at 16.65% (
Table 5). Sample 27 had the lowest ash content verified (11.18%). As a result, it would be impossible to create a pellet from swine waste and grape stalk that is suitable for ISO 17225-6:2014 classification.
The hydrogen content of the various mixtures was around 4% (
Table 6). The lowest recorded value was 3.92% for sample 31, and sample 32 had the highest hydrogen content verified (4.75%). Combining swine waste with other materials appeared to reduce the sulfur content, as no sulfur was detected in sample 23. Notably, an increase in the Lower Heating Value (LHV) of the pellets was observed when 30 percent of stalk content from red grape varieties was introduced into swine waste mixtures. In contrast, mixtures containing white grape varieties showed no significant variations. However, for higher stalk contents, there were no substantial enhancements in the LHV.
Mixtures of swine waste and two other biomasses were created to increase LHV while decreasing ash and sulfur content. These samples were made with the Touriga Franca cultivar stalks, which presented the highest LHV.
The sample with the lowest moisture content was sample 41 (6.35%). The mixture with the highest moisture content was sample 39 (
Table 7). In general, the moisture contents of these mixtures were lower than those of two-material mixtures. The content of ashes in the mixtures decreased, maintaining the values between 8.23% and 13.18%. Sample 44 was the mixture with the lowest content of ashes. Four mixtures were obtained with contents of ash lower than 10%, which allowed us to classify them as class B pellets. There was a slight decrease in the amount of ash in the mixtures, which remained between 8.23 and 13.18%. Sample 44 had the lowest ash content.
Four combinations were obtained with ash concentrations less than 10%, allowing them to be classified as class B pellets.
Sample 41 had a hydrogen content of 5.84% (
Table 8). However, this sample revealed a nitrogen content above that required to produce a class A pellet but sufficient to produce one of class B. The LHV of the pellets increased as compared to samples of only two components. The seeds and skins, in general, proved to be a better base than the stalks since they achieved higher LHV values. The maximum LHV was found in sample 45, with 19.10 MJ/kg. The sample with the lowest LHV confirmed was sample 39.
Sample 41 stood out as the best-performing mixture due to its balanced characteristics across key parameters. While other samples exhibited better performance in individual aspects—such as higher LHV or lower ash content—these advantages were often offset by unfavorable levels of nitrogen, sulfur, or other components. Sample 41 achieves the optimal compromise across all parameters, making it the best overall performer for producing class B pellets.
Mixtures of swine waste samples enriched with three biomasses, along with a swine waste sample containing 10% of each of the other biomasses, were investigated. The combination of four materials did not significantly enhance the ash content, which remained within the range from 9.63% and 11.96%, the sample 49 exhibiting the highest ash content (
Table 9). Nevertheless, two additional mixtures with ash levels low enough to be classified as class B pellets were identified.
Hydrogen contents of these samples varied between 3.22% and 4.66% for samples 52 and 47, respectively (
Table 10). Nitrogen contents verified for these mixtures varied between 1.32% for sample 51 and 1.86% for sample 53. LHV values stayed very close to 16 MJ/kg when Touriga Franca stalk was used as a second basis. The LHV of these mixtures ranged from 15.95 MJ/kg for sample 48 to 17.24 MJ/kg for sample 53.
These results show that by combining different biomass residues, it was possible to obtain pellets with enough quality to be commercialized as class B pellets.
3.3. Correlation Between Proximate Analysis, Elemental Composition, and Calorific Values of Biomass
To better understand the relationship between the proximate analysis results, elemental composition, and calorific values obtained for the biomass studied, a Principal Component Analysis (PCA) was performed.
The PCA yielded three principal components (PC) that accounted for 86.17% of the total variation in the original dataset (
Figure S1, Supplementary Material). PC1, which accounted for 44.4% of total original variance, correlated positively with the Moisture content and negatively with the fixed carbon, carbon content, HHV, and LHV values. PC2, which explained 25.6% of the total original variation, correlated positively with volatile matter and hydrogen content and negatively with ash and nitrogen content. None of the variables studied were significant in PC3, which explained 16.2% of total original variability.
These results show that there was a positive relationship between heat values and the carbon and fixed carbon content of the biomass samples studied, and a negative relationship between heat values and the moisture content of the samples (
Figure 1A). Because carbon and hydrogen content have the impact of boosting calorific value in the combustion process [
40], a positive correlation should be expected with the hydrogen content. This was not observed, however, due to the presence of the biochar sample with a high LHV and a very low hydrogen content (1.45%). In fact, if a PCA analysis was performed ignoring this sample, it could a positive correlation could be found between the calorific value and hydrogen as would be expected.
The distribution of the sample scores along PC1 and PC2 is presented in
Figure 1B. Most samples were clustered near the origin, but there were samples with high negative values of PC1, such as biochar (sample 11) and cork powder (sample 9). Cluster analysis was performed on the PC1 and PC2 scores of the biomass samples using Euclidean distance and Ward’s hierarchical agglomerative method. This method evaluates cluster distances, using an analysis of variance methodology seeking to minimize the sum of squares of any two hypothetical clusters that may form at each clustering phase.
Figure 1C shows a formation of 20 clusters (considering the cut-off on the 2.6-region obtained from the plot of Linkage Distances (
Figure S2, Supplementary Material). Apart from sample 36, forming an isolated cluster and sample 39, the first five clusters were entirely made up of samples with two mixed biomasses of swine waste and grape stalks. The sixth cluster was composed of samples 49–52, which were mixtures of four biomasses. Two of the samples from these four biomass mixtures, samples 47 and 48, were left out of this cluster, most likely due to their lower ash content, and they appeared in the eighth cluster, showing some similarities to samples 38, 42, and 44. Samples 20 and 37 formed the seventh cluster. The three following clusters were composed of isolated samples. Sample 10 and sample 8 were samples outlined in the projection factor plan. Because it had a significantly lower ash content, sample 5 was separated from all the other stalks in the 12th cluster. The 13th and 14th clusters were composed of samples 43, 45, and 46. Sample 45 was separated from the other two due to its higher LHV caused by the presence of cork powder. All these samples were made up of four biomass mixtures and contain biochar. The remaining clusters consisted of essentially isolated samples.
A multiple linear regression (MLR) analysis was performed using elemental composition data to understand its relative contribution to the LHV value. The MLR models were generated using the “best subset” method, optimized with Adjusted R
2. The Student
t-test (
p < 0.05) was used to evaluate the parameter estimates for all models. The resultant model was statistically significant (R = 0.84; F = 29.33;
p < 0.001), accounting for 71% of the LHV variance as showed in (
Table S1, Supplementary material).
The final regression model is expressed as follows:
where N, C, H, and O represent the nitrogen, carbon, hydrogen, and oxygen content of the samples, respectively. The intercept and coefficients quantify the relative contribution of each variable to the LHV.
Carbon content exhibited the greatest influence on the LHV, as indicated by its highest beta value and its zero-order regression coefficient (r = 0.712), accounting for 50.66% of its variability. This was followed by oxygen (r = − 0.433; 18.79%), hydrogen (r = − 0.115; 1.33%), and nitrogen (r = 0.053; 0.29%). The squared structural correlation coefficients further emphasize carbon’s predominant role (
; 71%), with smaller contributions from oxygen (
; 26%), hydrogen (
r; 2%), and nitrogen (
; 0.4%). These results align with
Figure 2 and
Table S2 (Supplementary Material). The product β × r makes it possible to calculate the partition of the regression effect into non-overlapping parts based on the interaction of the β coefficients and the zero-order correlation coefficients with the dependent variable, [
41], showing that, in this respect, the carbon content (59%) and oxygen content (14%) account for the large part of the variation in the regression equation, followed by hydrogen content (4%) and, finally, the nitrogen content (2%). These results clearly show that for these biomasses, the variation of the LHV is mostly explained by the variation of the carbon and oxygen content and, to a lesser extent, by the hydrogen and nitrogen content.
3.5. Energetic Valorisation
Three scenarios were considered based on key criteria for biomass utilization, including the use of all available biomass from wine and pig production, optimizing the value of the final product (biofuel), prioritizing mixtures with the highest LHV to maximize financial gains, minimizing ash content for cleaner combustion, and enabling the dilution of dangerous contaminating constituents (
Table 13).
The first scenario involves the utilization of all wine by-products in conjunction with swine waste. Leftover seeds and skins are further combined with swine waste and cork powder, while the remaining swine waste is mixed with cork powder and biochar. This scenario aims to maximize energy production by leveraging all available resources, making it the most resource-intensive approach. In contrast, the second scenario focuses on minimizing ash content, which is achieved by blending wine by-products with swine waste, including stalks, seeds, and skins, alongside sawdust. After utilizing the stalks, the remaining seeds and skins are combined with swine waste and cork powder. Additionally, mixtures of swine waste, cork powder, and sawdust are incorporated. This approach prioritizes cleaner combustion while resulting in slightly lower energy potential compared to the first scenario. The third scenario strikes a balance between multiple criteria by replacing sawdust in the initial mixture with cork powder. The resulting blend of swine waste, stalks, seeds and skins, and cork powder achieves a higher LHV and reduced nitrogen content while maintaining compliance with ISO standards. Although this scenario slightly increases ash and sulfur contents compared to the second scenario, these values remain within acceptable limits, enhancing financial benefits while ensuring the production of standardized pellets.
Using the mixtures in Scenario 1, there are 529,695 tons of oil equivalent (toe) of energy available to be valued in Portugal. Using these biomasses for heat generation and assuming a combustion efficiency of 90%, 5,588,677 MWh of heat might be produced. A Rankine cycle with 30% efficiency may produce 1,676,603 MWh of electricity in power production. With an 85% utilization factor, a cogeneration system for heat and electricity would result in a total of 4,037,819 MWh of combined energy (
Table S3, Supplementary Material). With the 2020 natural gas price of EUR 0.0527/kWh in Portugal, potential savings could reach EUR 294,523,282 (as detailed in Subsection 2.4 of the Materials and Methods section, along with the reference). Furthermore, based on the reference rates for wholesale electricity trading in Portugal for the same year (EUR 0.10/kWh as detailed in Subsection 2.4 of the Materials and Methods section, along with the reference), electricity generation could be valued at EUR 167,660,312, and EUR 268,966,261 in a cogeneration configuration (
Table S4, Supplementary Material).
The selection of these combinations increases profits by selecting mixtures with higher LHV, as well as being the optimum method for utilizing wine biomass. However, this scenario implies the production of non-standardized pellets for having ash contents higher than 10%, which makes this scenario impracticable.
In Scenario 2, Portugal has the potential to obtain approximately 485,900 toe of energy. Due to the utilization of mixtures with lower LHV, the potential energy available for conversion in this scenario is 8.5% less than in the first scenario. Nevertheless, it is crucial to highlight that the fifth criterion is satisfied in this case.
With a substantial heat output of 5,115,775 MWh, these mixtures demonstrate their efficiency in converting biomass into thermal energy. Additionally, they contribute significantly to electricity generation, producing 1,534,732 MWh. The most striking aspect is their role in cogeneration, where these mixtures combine to yield a total of 3,696,147 MWh of versatile energy (
Table S5, Supplementary Material). Portugal could be saving EUR 269,601,334 in natural gas if these mixtures were valued in heat production. In renewable electricity production, it might be worth EUR 153,473,245, and in cogeneration, it could be worth EUR 246,206,895 (
Table S6, Supplementary Material).
Scenario 3 makes 485,463 toe of energy available for conversion in Portugal, which is more 589 toe then the second scenario. The substantial heat generation of 5,121,994 MWh demonstrates the efficiency of these biomass combinations, particularly in meeting heating demands.
Moreover, the electricity production output of 1,536,598 MWh is significant. This clean energy source contributes to reducing the carbon footprint and dependence on fossil fuels in Portugal’s energy mix. It aligns with the country’s goals to increase the share of renewable energy sources in its energy production. The cogeneration results, totaling 3,700,641 MWh, highlight the versatility and multifaceted benefits of these biomass mixtures (
Table S7, Supplementary Material). Employing these mixtures for heat generation could result in substantial savings, potentially reaching EUR 269,929,092. When harnessed for electricity production, they could hold an economic value of approximately EUR 153,659,825, and in a cogeneration setup, their worth might extend to an impressive EUR 246,506,212 (
Table S8, Supplementary Material).
These figures represent gross potentials, meaning that expenses related to the collection, transportation, and processing of these biomasses, as well as equipment costs, need to be considered. However, the anticipated energy valorization values suggest a promising future for utilizing this biomass. Incorporating this portion of energy into Portugal’s total primary energy consumption would result in an 8.07% increase in the share of renewable energy sources (
Figure 3), a considerable amount that would undoubtedly help Portugal reduce its energy dependence. Furthermore, based on the average per capita electricity consumption in Portugal of 2348.4 kWh/year [
29], this amount of electricity would be sufficient to meet the energy needs of 654,482 people, which accounts for approximately 6% of the total Portuguese population. Notably, this production would be enough to supply the entire population of the Algarve region, estimated at 438,864 inhabitants in 2018 [
29].
To assess the practicality of this research, an analysis was conducted on a farm located in Portugal’s Douro region. This farm operates with a livestock production of 2000 pigs and an annual wine grape production of 500 tons. For privacy considerations, this establishment is referred to here as the “Case Study Farm”.
By implementing the third scenario while considering the availability of dry biomass resources, including 473,146.76 kg of swine slurry, 10,032.34 kg of stalks, and 32,322.87 kg of seeds and skins, along with the acquisition of 231.6 tons of cork powder and 199.2 tons of sawdust, resulted in a total biomass resource of 946.29 tons available for pellet production. This amounted to 365.31 toe of available energy at the Case Study Farm.
Using these biomasses for heat production might result in approximately 3854 MWh of heat output. When employed for electricity generation, it has the potential to produce 1156 MWh of clean energy. Considering a selling price of EUR 0.20/kg, pellet sales could generate a yearly financial benefit of EUR 189,258 or a potential gain of EUR 115,627 per year from electricity production.
To assess investment uncertainty, a financial analysis was conducted for pellet production. Some economic indicators were calculated based on an annual inflation rate of 3% and a project lifespan of 20 years. The capital expenditure (CAPEX) for the pellet production plant was estimated at EUR 283,937.80, based on market quotations for equipment and infrastructure, including preprocessing, drying, pelletizing, and storage facilities. The operational expenditure (OPEX) of EUR 137,214.99 per year included costs for energy, maintenance, licensing, consumables, labor, and the acquisition of cork powder and sawdust. These estimates were scaled to the plant’s production capacity of 946 tons per year and adjusted to reflect current market conditions. At the end of 20 years, the net present value (NPV) indicates that not only was the initial investment of EUR 283,937.80 recovered, but also an additional surplus of EUR 689,666.32 was generated. The capital investment at the Case Study Farm was valued at an annual rate of 17.78%, which was represented by the internal rate of return (IRR). Additionally, it was determined that the initial investment of EUR 283,937.80 had a payback period of 6.36 years.