1. Introduction
The world is facing significant challenges in meeting its energy needs while minimizing its environmental impact. The finite nature of fossil fuels and the increasing concerns about greenhouse gas emissions have led to the exploration of alternative energy sources [
1]. Energy from biomass is one such alternative that has gained widespread attention due to its potential to reduce dependence on fossil fuels while reducing greenhouse gas emissions. Biomass can be used for various purposes, including heat, power, and transportation fuels, with biofuels being a key area of interest.
Hydrothermal liquefaction (HTL) has emerged as a promising technology for the conversion of biomass into biofuels [
2]. HTL uses high temperatures and pressure to convert wet biomass into a liquid bio-crude that can be refined into transportation fuels. The process has several advantages over other biofuel production technologies, including the ability to process a wide range of feedstocks because of its acceptability of water, high energy efficiency as no prior drying of feedstock is required, and the potential to produce high-quality biofuels. In this study, EtOH and FA are used as part of the reaction medium.
The inclusion of ethanol (EtOH) in the HTL process is a significant factor deserving of consideration. It appears that EtOH exerts a promoting effect on bio-crude yield, as shown in previous studies [
3,
4]. This might be due to EtOH’s ability to de-polymerize the lignin structures and esterify the acids from the de-amination of amino acids [
5]. However, derivatives from EtOH, such as ethyl esters, might also become part of the bio-crude phase, resulting in better yields with higher concentrations. In the reaction medium comprising water, EtOH, and FA, there is a potential to produce ethyl formate, a compound that might be incorporated into the bio-crude. To elucidate this phenomenon and accurately assess the contribution of EtOH to the bio-crude, further investigations are warranted as part of future work, including GC–MS analysis to assess compounds that might be derived from EtOH and FA and how much they contribute to the overall bio-crude yield.
In 2008, a significant advancement in the hydrothermal liquefaction process was introduced by incorporating formic acid (FA) as an in-situ hydrogen donor. This innovative approach, known as the lignin-to-liquid process, initially focused on lignin-rich feedstocks, resulting in impressive bio-crude yields ranging from 60–80%, as shown in the study conducted by Kleinert et al. [
6]. Since its inception, formic acid-assisted HTL has undergone continuous refinement, with adaptations including the utilization of water as the reaction medium and different feedstocks such as Spirulina [
7]. Notably, Hegdahl et al. demonstrated the adaptability of this method to digested sewage sludge despite its original design for lignin-rich feedstocks [
8]. This evolution highlights the versatility and potential of the lignin-to-liquid process in accommodating a broader range of feedstock sources and underscores its significance in the context of bio-crude production.
In recent years, co-HTL has gained attention as a means of enhancing and increasing the production of biofuels from biomass, especially in terms of enlarging the feedstock base for continuous production of biofuels [
9]. Co-HTL involves the simultaneous processing of two or more feedstocks, resulting in a potential interaction that can affect the overall yield and quality of the resulting biofuels. In this study, digested sewage sludge and lignocellulosic waste from spruce are investigated regarding an interaction effect on yield and bio-crude quality.
Digested sewage sludge and lignocellulosic waste are both valuable feedstocks for the production of renewable fuels and chemicals. Digested sewage sludge is a byproduct of wastewater treatment processing, which produces biogas by anaerobic digestion. The solid residue still contains high levels of organic carbon, including recalcitrant components and microbial biomass [
10]. It is to some degree used as fertilizer but is also disposed of in landfills or incinerated, which can lead to negative environmental impacts. However, the use of digested sewage sludge as a feedstock for HTL can convert this waste into a valuable resource. Similarly, lignocellulosic waste, such as sawdust and wood chips, is a byproduct of the forestry and wood industries [
11]. These wastes are often burned for energy or disposed of in landfills. However, they can also be converted into biofuels and high-value chemicals through HTL. The use of these waste streams as feedstocks for HTL has the potential to not only reduce waste and emissions but also provide a sustainable source of biofuels and chemicals. By combining these two feedstocks in the HTL process, there could be a potential to achieve a synergistic effect, leading to higher yields and better-quality biofuels.
The inherent heterogeneity observed in the feedstocks [
8], particularly in the DSS, presents a notable source of uncertainty. The variability in DSS composition can be attributed to several factors, including geographical location, seasonal variations, and the specific sources of wastewater. These dynamic environmental conditions can lead to significant fluctuations in the chemical and physical properties of the DSS feedstock, introducing variability between different batches. Even within the same batch, DSS exhibits heterogeneity due to the complex nature of sewage sludge. Replicate experiments conducted on such heterogeneous feedstock are susceptible to inherent variations, which may impact the reproducibility of results. Consequently, it is imperative to acknowledge and address the uncertainties stemming from feedstock heterogeneity in the context of experimental design and data interpretation, as they can influence the robustness and reliability of the findings.
In this paper, we present the results of co-HTL of digested sewage sludge and lignocellulosic waste (spruce) and demonstrate the potential for this approach to increase and enhance the production of biofuels from biomass. This study investigates the effects of different operating conditions on the co-HTL process and evaluates the properties of the resulting bio-crudes.
To investigate the effects of operational conditions, several multivariate methods are employed, including factorial design, principal component analysis (PCA), and partial least squares regression (PLS). Additional experiments are being carried out to further elucidate the effect of the feedstock ratio. These methods enable a thorough exploration of the complex interactions between various parameters, such as temperature, hydrogen donors, and reaction medium, and their impact on the final product.
2. Materials and Methods
2.1. Digested Sewage Sludge
The digested sewage sludge (DSS) utilized in this study was obtained from the Bergen Biogas Facility in Bergen, Norway, following the completion of anaerobic digestion of sewage sludge. The Bergen Biogas Facility is owned by the municipality and sources its input for biogas production from wastewater treatment plants in the vicinity of Bergen [
12]. DSS typically comprises a complex mixture of organic and inorganic components, including residual organic matter, microorganisms, nutrients, and various contaminants. It commonly contains organic materials in the form of cellulosic and lignocellulosic materials, proteins, lipids, and other organic compounds originating from human waste, food residues, and household products. Additionally, DSS might contain inorganic constituents such as metals, salts, and minerals originating from wastewater. Before its acquisition, the DSS feedstock underwent a dewatering process at the Bergen Biogas Facility. Their analysis showed that this particular batch of DSS contains 74.6% water. A total of 61.0% of the dry matter is organic.
2.2. Lignocellulose
Lignocellulosic biomass (LC) is predominantly composed of cellulose, hemicellulose, and lignin, which form the structural basis of plant cell walls. Cellulose provides strength, hemicellulose offers flexibility, and lignin imparts rigidity [
13]. LC also contains minor components like extractives (oils and resins) and ash. The lignocellulosic component of the total feedstock in this study consisted of samples of cut spruce collected from a logging field located on Osterøy, an island situated on the coast of Western Norway. The collected spruce samples were milled to produce coarsely ground sawdust particles measuring approximately 0.5 cm
3 in size. Subsequently, the sawdust was oven-dried at a temperature of 105 °C until the weight change did not exceed 1%, with measurements taken at 24-h intervals.
2.3. Hydrothermal Liquefaction
The HTL process was carried out as shown in
Figure 1 using the following procedure: 1 g of organic sample (calculated by the organic content of DSS provided by Bergen Vann), 1 mL of FA, and enough distilled water/EtOH to fill to constant loading (6.7 mL) were combined in a non-stirred 25 mL Parr reactor from the 4740 series, manufactured by Parr Instrument Company in Moline, IL, USA. The reactor was tightly sealed and placed inside a preheated Carbolite Laboratory high-temperature oven, where it was heated within a temperature range of 320–360 °C for a duration of 3 h. Subsequently, the reactor was allowed to cool to room temperature.
To separate the gaseous phase, the reactor was vented while the solid phase was filtered using a glass fiber filter (Whatman, GF/A, Whatman plc, Maidstone, Kent, UK). Filtration was performed using a solvent mixture composed of ethyl acetate (EtOAc) and tetrahydrofuran (THF) in a ratio of 9:1 (v/v). The organic phase was separated from the aqueous phase by decanting and then subjected to further drying with sodium sulfate. The solvent was subsequently evaporated using a rotary evaporator at 40 °C and 212 mbar. The weight of the reactor was measured before and after venting to determine the yield of the gaseous phase, while the weight of the filtration equipment was measured before and after filtration and drying to determine the solid yield. Finally, the bio-crude yield was determined by weighing the mass after the evaporation of the solvent.
It is important to note that water recovery was not measured due to the substantial water loss during the evaporation of EtOAc, which made accurate measurements impractical in this particular step of the process.
2.4. GC–MS
For the GC–MS analysis of all bio-crude samples, an Agilent Technologies 7890A GC equipped with an auto-sampler and an Agilent 5977A MSD detector were employed. Injection of the samples was carried out in splitless mode, utilizing an injection temperature of 280 °C. A 30 m HP-5 ms column from Agilent Technologies (Santa Clara, CA, USA), featuring an inner diameter of 250 μm and a thickness of 0.25 μm, was utilized. The temperature program for the column included an initial temperature of 50 °C for 2 min, followed by a ramping rate of 10 °C per minute up to 200 °C, and then a further increase at a rate of 20 °C per minute up to 300 °C, which was maintained for 5 min.
The mass spectrometer detector was equipped with an ion-source temperature of 230 °C and covered a mass range from 25 to 400 μ. Positive ionization was employed at 70 eV. Compound identification was conducted using Agilent MassHunter Qualitative Analysis version 10.0 coupled with the NIST 2.0 library. A solvent delay of 5.50 was implemented to ensure accurate detection and analysis of the bio-crude samples.
A silylation process was employed to enhance the separation of the compounds during GC analysis. Firstly, 0.005 g of the bio-crude sample was added to a vial and dissolved in a mixture of ethyl acetate/tetrahydrofuran (EtOAc/THF) (1.5 mL, 9/1, v/v). A total of 0.01 mg/mL of hexadecane was included as an internal standard. A total of 1 mL was extracted from the mixture before adding 0.150 mL of pyridine and 0.150 mL of N,O-bis(trimethylsilyl)trifluoroacetamide. The vials were securely capped and heated at 70 °C for 30 min, followed by cooling to room temperature.
Subsequently, 0.700 mL of the silylated bio-crude was transferred to a new vial and mixed with 0.700 mL of pentane, resulting in a solution containing approximately 1.3 mg/mL of the silylated oil and 4.2 μg/mL of the internal standard. The samples were then cooled at 5 °C overnight and subsequently filtered through a 0.45 μm Puradisc NYL filter before being subjected to GC–MS analysis. This silylation procedure and subsequent filtration step helped to improve the separation and preparation of the bio-crude samples for analysis.
2.5. Elemental Analysis and Heating Value Calculations
Two replicates were analyzed using a Vario EL III instrument in CHNS mode to conduct elemental analysis on the bio-crude samples. Helium was utilized as the carrier gas during the analysis. The oxygen content was estimated as the difference in elemental composition. The obtained data from the elemental analysis were utilized to characterize the bio-crudes based on the hydrogen-to-carbon (H/C) and oxygen-to-carbon (O/C) ratios. Additionally, these ratios were employed to calculate an estimated higher heating value (HHV), providing valuable knowledge of the energy content of the bio-crudes. This calculation is proposed by Channiwala et al. [
14]. The energy recovery was calculated by (HHV
bio-crude × mass
bio-crude)/(HHV
Feedstock × mass
Feedstock). Formic acid and ethanol are not included in the heating value calculations, so a full energy balance for the reaction system is not attained.
2.6. Experimental Design
The experimental setup followed a 2
4−1 fractional factorial design, encompassing four main factors: the mixture ratio of DSS and LC, temperature, EtOH concentration, and the presence of FA. Two center-point experiments were conducted, one including FA and the other without it. The variable ranges for these factors were established based on prior screening studies involving lignocellulosic feedstocks (Alper et al., 2019 [
15] and Xu et al., 2012 [
16]), ensuring comprehensive coverage of oil conversion across a wide spectrum. To ensure the reliability of the results, replicated experiments were conducted to assess their reproducibility. Complete details regarding the experimental conditions can be found in
Table 1.
Based on the results from the fractional factorial design, a new set of experiments was conducted to delve deeper into the analysis of one specific variable while keeping the remaining variables at the levels that resulted in the highest response (
Table 2). By isolating the variable that explained the mixture ratio of the feedstocks and observing its effects within the context of the previously identified optimal levels, a more detailed understanding of its contribution to the overall response is obtained, and a potential curvature in the response could be elucidated. This approach offers the most effective means to demonstrate the multiplicative effect resulting from the mixture of feedstocks and ascertain whether this effect is synergistic or antagonistic.
2.7. Multivariate Data Analysis
To examine the outcomes of the experimental design, principal component analysis (PCA) was conducted using Sirius 11.0 software, resulting in the generation of a biplot. The biplot visually represents the correlations between different factors. In the biplot, factors that correlate positively are positioned on the same side of the origin, as indicated by a plus sign in the results [
17]. Conversely, factors that correlate negatively are positioned on opposite sides of the origin. Factors positioned at 90° angles from each other through the origin are considered uncorrelated. The distance of a factor from its origin reflects its impact on, or sensitivity to, other factors. The further away a factor is from its origin, the greater its influence or susceptibility to other factors.
In addition to PCA, the partial least squares (PLS) method was employed using Sirius 11.0 software. This approach quantifies the influence of each factor on the target outcome. The statistical significance of a PLS model was assessed through key metrics, including R
2 (coefficient of determination), RMSECV (Root Mean Square Error of Cross-Validation), and the number of PLS components employed. R
2 indicates the proportion of variance in the dependent variable that the model explains, with higher values suggesting a more robust model fit. RMSECV quantifies the accuracy of predictions, representing the average error between actual and predicted values, and lower values indicate better predictive performance. The number of PLS components used in the model influences its complexity. Increasing the number of components may lead to overfitting, where the model performs well on the training data but poorly on new data. Therefore, a balanced approach is required to select an optimal number of components that capture the essential variance without introducing excessive complexity [
18].
A statistically significant PLS model exhibits high R2 values, low RMSECV values, and an appropriate number of PLS components. These factors collectively ensure the model’s ability to effectively represent the underlying relationships in the data, providing reliable predictions and meaningful insights.
4. Conclusions
The best experiment from the fractional factorial design shows that mixtures with more DSS, 350 °C, formic acid present, and 50 wt% EtOH give the highest bio-crude yield of 58.9 wt%. For comparison, pure DSS and LC run at these process conditions resulted in a bio-crude yield of 52.5 wt% and 48.3 wt%, respectively. The organic content in the solid residue follows an opposite trend, where the experiment with the highest bio-crude yield also shows the lowest organic content in the solid residue of 3.7 wt%.
The elemental ratios of hydrogen, oxygen, nitrogen, and carbon do not exhibit significant variations across the experiments, with HHVs ranging from 29.3 to 35.6 MJ/kg. Since the HHV is relatively constant across the experiments, the proportion of energy recovered relative to the total energy content in the bio-crude is primarily determined by the amount of bio-crude obtained. The energy recovery achieved in the bio-crude relative to the feedstock ranges from 41 to 103%.
The PLS model was constructed using two PLS components covering 94.25% of the total variance, resulting in a root mean square error of cross-validation (RMSECV) of 6.58 for the prediction of crude yield. The regression equation reveals a synergistic effect by mixing feedstocks, as the quadratic term of the regression equation for mixture ratio shows a negative coefficient. This effect is also shown in the biplot from the principal component analysis. The experiment conducted with a 50/50 mixture ratio and all other levels at high reaction condition levels exhibited the highest correlation with bio-crude yield, and the correlations diminished as the ratios deviated from the 50/50 mark.
GC–MS data revealed novel peaks within the chromatograms of combined feedstocks that were not discernible in either of the individual feedstock chromatograms. These distinctive peaks are associated with new phenolic compounds and further show the effect of combining different feedstocks.
This research advances the field of co-HTL by demonstrating its viability in converting DSS and LC into bio-crude with an impressive yield that stands in the upper region of reported values in the available literature. In addition, this work not only contributes a valuable method for optimizing co-HTL conditions but also underscores the synergistic effects of co-processing DSS and LC, offering both high yield and consistent energy content. The implications of this study extend to advancing sustainable biofuel production through the efficient utilization of diverse biomass feedstocks, paving the way for environmentally friendly and economically viable energy solutions.