Muhammad Ahsan
Amjed
*ab,
Filip
Sobic
a,
Matteo C.
Romano
a,
Tiziano
Faravelli
b and
Marco
Binotti
*a
aDepartment of Energy, Politecnico di Milano, via Lambruschini 4, 20156, Milano, Italy. E-mail: muhammadahsan.amjed@polimi.it; marco.binotti@polimi.it
bDepartment of Chemistry, Materials, and Chemical Engineering, “G. Natta”, CRECK Modelling Lab, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133, Milano, Italy
First published on 22nd July 2024
Pyrolysis has become one of the most attractive options for converting carbonaceous biomass into bio-oil or biochar. This study explores a novel solar pyrolysis process intended to produce both bio-oil and biochar, thereby improving carbon efficiency. Aspen Plus and SolarPILOT were used to model a 10 MW biomass pyrolysis plant thermally sustained by hot particles from a falling-particle solar tower receiver. A yearly analysis was carried out for three configurations to estimate the annual production of oil and biochar. The results showed that the hybrid plant, combining solar receiver and biochar backup combustor, leads to the lowest cost of bio-oil (18.7 € per GJ, or 0.29 € per kg) and a carbon efficiency of 83%. Whereas, the plant fully sustained by solar power achieves a carbon efficiency of 90%; however, it results in a significantly higher cost of bio-oil (21.8 € per GJ, or 0.34 € per kg) due to the larger size of particle storage and a lower capacity factor of the pyrolysis plant. In comparison, a conventional pyrolysis plant with no biochar production yielded the most expensive option in terms of the cost of produced bio-oil (27.5 € per GJ) and features the lowest carbon efficiency (74%). Sensitivity analysis shows that the pyrolyzer Capex, operational cost, biochar market price, plant availability and discount rate significantly affect bio-oil production cost.
The pyrolysis process is endothermic and requires an external source of energy to heat up the feedstock and break down the molecular structure of the carbonaceous material. This thermal input can be supplied by various sources, e.g. direct combustion of its by-products, electrical resistive heating or solar thermal heating.6,7 The limit of conventional pyrolysis processes sustained by combustion is the loss of part of the biogenic carbon as CO2. Supplying the process heat via electric or solar heating allows for the reduction of CO2 emissions and improves the carbon efficiency and the amount of bio-based products.
Concentrated solar power (CSP) systems use a series of mirrors that concentrate the solar radiation towards a receiver where solar energy is converted into thermal energy, which can either be used directly or to produce electricity using a thermodynamic cycle. The possibility of storing heat in a low-cost thermal energy storage (TES) system for later use allows for decoupling the availability of solar radiation from thermal power production. In contrast to solar photovoltaics (PV), CSP systems use almost the entire spectrum of solar radiation to produce heat typically in the range of 400–2000 °C, which may be used to drive chemical reactions such as pyrolysis.8,9 The integration of CSP into conventional pyrolysis may help fulfil the energy requirements of the endothermic pyrolysis reactor.6,10,11
Joardder et al.12 conducted a lab scale study, in which solar radiation was directly focused on a pyrolysis reactor. The results show that solar-based pyrolysis can cut 33% of fuel production costs and 32.4% of carbon emissions. Moreover, combining CSP with conventional pyrolysis allows for the storage of solar thermal energy in terms of transportable, high energy density and upgraded fuel.10,13 It also reduces reactor dependency on its by-product combustion, which ultimately reduces its carbon footprint.14,15
Other studies have investigated the integration of CSP in the conventional pyrolysis of biomass using direct radiation on a reactor in a lab-scale environment in which multiple types of biomasses went through a thermochemical conversion process with CSP assistance.6,16–21 In addition to biomass, the solar pyrolysis of other feedstocks was assessed. For instance, Zeaiter et al.22 assessed the pyrolysis of scrap rubber using Fresnel lenses to concentrate the radiation over a tube filled with feedstock, with and without a catalyst. Hosseini et al.8 recently modeled coal drying and pyrolysis using CSP under four different scenarios considering different collectors, pyrolysis temperatures, and heat transfer fluids (HTF). They found that the integration of solar towers at optimized sizes of solar multiple and storage can meet 12.8% of the annual thermal energy demand for the drying and pyrolysis of coal.
In some lab scale studies,6,12,16,19–23 fixed bed or rotating pyrolysis reactors were used with solar radiation directly concentrated onto the reactor surface by Fresnel lenses or parabolic dishes, which attained low and medium temperature ranges from 250 to 600 °C based on design parameters and pyrolysis reactor requirements.
Many types of CSP systems can be used for solar pyrolysis, e.g. dishes, linear concentrators and power towers.10,24 Contrarily to conventional power towers, which typically adopt molten salts as a heat transfer medium, this work considers advanced falling particle tower systems, recently proposed as advanced solutions for next-generation CSP systems.25 The main idea behind the falling particle receiver is to avoid the maximum temperature limit given by the stability of the molten salts (∼565 °C), allowing for direct heating of the particles, which enables increasing the peak flux on the receiver from 1 MW m−2 in the case of conventional receivers to at least 2 MW m−2.26 In the case of pyrolysis, another important advantage is the possibility of using these particles directly in the reactor, thus avoiding the use of an intermediate heat exchanger. A recent comparative study by Jie Ling et al.27 also suggested the above-mentioned claim and proposed a CSP beam-down tower and solid particle receiver with a fluidized reactor for high efficiency and high-temperature solar pyrolysis operations.
Previous studies on CSP-based pyrolysis predominantly focus on small-scale laboratory experiments involving direct irradiation on the reactor. The solar-driven pyrolysis using direct radiation faces several challenges, including uneven heat distribution within the reactor, thermal stress caused by variations in solar flux, and lower oil yields due to the slow pyrolysis of biomass. Furthermore, there is no control over the temperature and heating rates, and the process stops working in the absence of sunlight.7,10 To address these issues while maintaining the concept of fast pyrolysis using a fluidized bed reactor, indirect solar CSP heating was adopted through solid PHCs. This solution also offers thermal storage for longer periods in the absence of sunlight, ensuring the smooth operation of pyrolysis.
There is a notable absence of studies examining CSP-based pyrolysis with solid PHC falling particle receivers on an industrial scale. Moreover, there is a significant gap in the literature regarding the techno-economic analysis of such industrial-scale solar-assisted pyrolysis plants. Our work aims to address this gap and serves as a starting point for future projects in this area. To the best of our knowledge, this work explores for the first time the process integration of an industrial scale biomass pyrolysis plant with a falling particle solar tower system using the same solid particle heat carrier (PHC) for the biomass pyrolysis reactor and as a heat carrier in the solar receiver. The main purpose of this study, which is conducted as part of the EU Pysolo project,28 is to perform a techno-economic assessment of 10 MWth fast pyrolysis plants for bio-oil and biochar production integrated with a falling particle solar tower system, starting from existing industrial scale models for the pyrolysis process and CSP system. The following three cases are assessed and compared through techno-economic indicators:
• Conventional pyrolysis process: reference configuration, where the heat for the pyrolysis process is supplied through the combustion of a fraction of the pyrolysis products (char and pyro-gases).
• Solar-based pyrolysis process: heat for the pyrolysis process is provided only by solar heat produced using the CSP system equipped with a Thermal Energy Storage (TES) system, and all the produced biochars is exported as products.
• Hybrid pyrolysis process: heat for the pyrolysis process is supplied either from the CSP system or, when no solar heat is available, from the combustion of a fraction of the pyrolysis products, resulting in the export of most of the biochar produced.
Woody biomass of poplar tree that grows abundantly in Europe, Canada and South America is considered biomass feedstock. The plant is sized to convert a biomass input of 50 dry t per h (10 MWth on a lower heating value basis). The relatively small plant size is representative of a scenario of multiple pyrolysis plants erected at a biomass point source aimed at generating a high volumetric energy density product to be transported to a central refining plant. Transportation of compact products instead of bulky biomass from point sources helps reduce both transportation costs and carbon emissions. The main assumptions for modeling the pyrolysis plant are as follows:
(1) The plant operates under steady-state conditions.
(2) All the plant equipment is modelled with a zero-dimensional approach.
(3) Particle Heat Carriers (PHC) used in the model are considered chemically inert in the pyrolysis reactor.
A block diagram of the conventional pyrolysis section is depicted in Fig. 1A. Biomass is initially dried using hot flue gases obtained from the combustion of pyrolysis products. Then, it is fed to a fast pyrolysis fluidized bed reactor (FBR) together with hot PHC and used as a heat carrier. The solid products and the PHC exiting the pyrolysis reactor are then sent to a combustor, where the combustion of char increases the PHC temperature before it is recirculated to the FBR. The gaseous products exiting the FBR are then cooled and separated in non-condensable pyro-gases, partly recirculated to the reactor and partly combusted to thermally sustain the process.
A schematic of the solar-driven pyrolysis plant is reported in Fig. 1B. In this configuration, the solid PHC and char exiting the FBR are separated. The PHC particles are sent back to the CSP section of the plant, and the biochar is extracted from the plant as an additional product. The CSP section includes the falling particle receiver on top of the power tower and the TES system. In the receiver, the solid particles acting as PHC are directly irradiated to raise their temperature and then sent back to the FBR. The solar field is modelled using SolarPILOT,31 a software developed by NREL for the design and performance simulation of power tower heliostat fields, while the falling particle receiver is modelled with an in-house code.32
For the solar-based pyrolysis plant, two PHCs are considered (sand and ceramic particles), and two different operational modes are considered (solar only and hybrid mode). For each plant, a yearly product yield is assessed, together with annual emission, CSP plant efficiencies, energy conversion efficiencies, minimum fuel sale price (MFSP) and availability.
The biomass considered in the model is hybrid poplar wood, with thermochemical analysis data presented in Table 1. For ambient air, summer conditions are assumed (32.2 °C, 10132 bar, 75.5% relative humidity).29 The thermochemical properties of all the species used in Aspen Plus to model pyrolysis products are given in ESI Section 1.†
Elemental analysis (% wt on dry basis) | |
C | 50.94 |
H | 6.04 |
O | 41.90 |
N | 0.17 |
S | 0.03 |
Ashes | 0.92 |
Proximate analysis (% wt on dry basis) | |
Volatile matters | 84.88 |
Fixed carbon | 14.2 |
Ashes | 0.92 |
Sulfur analysis (% wt on dry basis) | |
Pyritic sulfur | 0.03 |
Sulphate | 0.00 |
Organic sulfur | 0.00 |
Calorific values (MJ kg−1) | |
HHV | 14.0 |
LHV | 12.3 |
Humidity, dry basis | 30% |
Pyrolysis plant energy conversion efficiency ηpyro plant is calculated by applying eqn (1):
(1) |
The CSP plant solar-to-thermal efficiency ηsol–th is calculated by applying eqn (2) as the product of the solar field optical efficiency, ηopt, and of the falling particle receiver thermal efficiency, ηth,rec. rec and PHC,rec are the solar power incident on the receiver and thermal power delivered to PHC, respectively. DNI is the direct normal irradiance and Ah is the total heliostat area.
(2) |
Carbon efficiency εC is calculated using eqn (3), where yC,prod,i and yC,biom are the carbon content in ith product and biomass, respectively.
(3) |
The thermal energy usage efficiency of the CSP-based pyrolysis plant ηth,use is calculated using eqn (4), where def is the total solar power loss through defocusing.
(4) |
Emission to oil ratio (ETO) is calculated using eqn (5). E and Oil are the CO2 emissions and oil production in kg and GJ, respectively, from the pyrolysis plant.
(5) |
Considering that biogenic CO2 emissions are climate neutral, CO2 emission credits associated with biochar production are also computed as net negative emission to oil ratio (ETOnet) using eqn (6).
(6) |
• Section 1, Biomass pretreatment: The first section comprises the preparation of biomass before the pyrolysis process via drying and grinding.
• Section 2, Pyrolysis: This section consists of a fluidized bed pyrolysis reactor, where dried biomass is converted by contact with hot PHC and fluidizing gases.
• Section 3, Solid removal: It encompasses the cyclone filter, which is responsible for the separation of volatile products from solids entrained from the reactor at high temperatures.
• Section 4, Bio-oil recovery: In this section, quenching columns are used for the condensation and collection of bio-oil along with other auxiliaries and bio-oil filters.
• Section 5, Combustion: This section includes the combustor block, where biochar, a portion of pyrolytic gases, and retentate of the filter are burnt to heat the circulating PHC.
A single-pass direct feed dryer is assumed to consist of a rotating drum where biomass is put in contact with the hot gases, which provide the required heat for drying. The biomass is collected from the outlet as the drum rotates, while the humid gases are vented into the atmosphere. For this type of dryer, the inlet hot gases have a temperature between 230 °C and 1100 °C and an outlet temperature between 70 °C and 110 °C to prevent condensation of moisture in the vapor stream.33
The biomass temperature and residual moisture at the dryer exit (stream 2) were assumed to be equal to 71.7 °C and 10%, respectively. The temperature of the hot gases (stream 16) fed to the dryer is reduced by mixing ambient air (stream 25) with flue gases exiting the combustor (steam 23). ESI Table 3† reports the parameters assumed in the simulation.
An RYield reactor with a fixed outlet composition, operating temperature and pressure obtained from the study by Jones et al.29 is adopted in Aspen Plus to model the reactor. Table 2 describes the product yield from the pyrolysis reactor.
Species | Concentration (% wt) | Species | Concentration (% wt) | Species | Concentration (% wt) | Species | Concentration (% wt) |
---|---|---|---|---|---|---|---|
H2O | 5.89 | C2H6 | 0.64 | C3H6O2 | 1.30 | C12H8O | 0.35 |
H2 | 0.07 | C3H8 | 0.02 | C6H6O2 | 1.09 | C20H26O8 | 1.43 |
O2 | 0.19 | C4H10 | 0.00 | C10H12O2 | 0.78 | C21H26O8 | 0.31 |
N2 | 0.04 | C2H4 | 1.46 | C6H10O5 | 0.93 | C8H11N | 0.01 |
CO | 34.88 | C3H6 | 1.02 | C12H22O11 | 5.29 | S | 0.01 |
CO2 | 32.87 | C4H6O2 | 0.87 | C20H28O2 | 0.47 | C | 2.42 |
CH4 | 4.63 | C8H8O3 | 1.09 | C16H16O2 | 1.71 | ASH | 0.22 |
Downstream of the pyrolyzer, all products flow into the solid separation unit, which consists of a cyclone filter to remove almost all solids from hot vapors. Solid materials, mainly PHC, char and ashes, are transferred to the combustor, while volatile matter is directed to the condensation unit. The PHC separation efficiency is assumed to be 99.9% while the char and ash separation efficiency is assumed to be 92% because of the finer particle sizes. ESI Table 4† depicts the assumptions and parameters used in this plant section.
Before storage, bio-oil passes through the filtration process, where solid particles are removed along with a small portion of oil (stream 8). This is burnt together with other components in the combustor to fulfil the energy needs of the reactor.
Combustion is carried out with ambient air, with 25% excess with respect to stoichiometric combustion. Thermal losses in the combustor are assumed to be equal to 1% of the biomass LHV. Downstream of the combustor, the hot PHC is separated by applying a cyclone and recirculating it to the pyrolyzer. The resulting hot gas is cooled and diluted with air to achieve 345.5 °C at the dryer inlet. Entrained ashes along with a small quantity of PHC are separated by filtration air dilution, and the drying gas temperature is tuned to achieve a target dryer exit temperature of 71.7 °C.29 ESI Table 6† lists the assumptions and parameters used to model the combustion section.
(1) Introduction of CSP block to replace the biochar combustor;
(2) Introduction of an additional combustor for the combustion of excess gases for biomass drying;
(3) Introduction of biochar cooling and recovery unit.
Fig. 3 shows a schematic of the integrated process considering the above-mentioned changes. The pyrolizer outlet stream (stream 4) passes through a cyclone, where solid particles are separated from pyrolysis vapors. Solids are then sent to a char-PHC separator, where char is recovered and PHC is returned to the CSP loop. Char/PHC separation is expected to be a challenging process that may be obtained using different solid/solid separation techniques by exploiting the different particle sizes and densities. In the presented analysis, char separation efficiency from PHC is assumed to be equal to 99.99%. A PHC make up (PHC-MK) is considered to account for losses during the separation of char and PHC. The separated PHC and make-up PHC flows are heated up to the same target temperature as the baseline plant (609 °C) in the CSP unit. ESI Table 7† illustrates the parameters used in CSP-based pyrolysis after changes in the conventional pyrolysis model.
Fig. 3 Schematic diagram of the modelled CSP-based and hybrid pyrolysis plants. Dashed lines represent hybrid system operation when solar heat is not available. |
In the CSP-based process, one combustor block is included to generate hot gases for the dryer. In this case, as char is recovered as an additional product, only light gaseous products are burnt in the combustor. The hot gases exiting the combustor are partly cooled before feeding the dryer via heat exchange and air mixing to avoid the loss of biomass volatiles in the drying process.
The plant is assumed to be located in Seville (Spain). DNI data are taken from SolarPILOT.31 The thermal efficiency under any given condition is obtained using the in-house model described by Pasqualotto et al.32 This model allows estimating the particle receiver thermal efficiency for a given geometry and given particle type and size. It accounts for the drag force effect on the particles, considers a 2D discretization of both solar flux on the receiver and particle properties inside the curtain, and considers the variable reflectivity of the wall behind the curtain. In this work, the sizing procedure is consistent with that used by Pasqualotto et al.,32 while two types of PHCs are considered, namely CARBO ACCUCAST ID 50 (spherical sintered-bauxite particles) and sand, whose properties are shown in ESI Table 9.† ACCUCAST ID 50, used at Sandia Laboratories,34 was identified as the best candidate for particle receivers due to its high solar absorptivity. However, sand is considered a lower-cost alternative35 although featuring a considerably lower solar absorptivity, resulting in lower receiver thermal efficiency.
The design and simulation of the solar field are carried out using SolarPILOT.31 Five solar field designs are considered by varying the solar multiple (SM). The SM characterizes the relative size of the CSP plant with respect to the pyrolysis plant, and it is defined as the ratio between the heat provided by the solar receiver in design conditions (i.e. 21st June, Solar noon) and the heat demand of the pyrolysis unit (1.73 MWth):
(7) |
The tower height for different SM values is obtained by linear interpolation between 25 m for a 1.2 MWth receiver, as reported by Coventry et al.,36 and 55 m for a 20 MWth receiver, as reported by Frantz et al.37 All the other input parameters, presented in Table 3, are common for all the fields. The required incident power on the receiver for each solar field is determined by dividing the product of the pyrolizer power requirement and the theoretical value of the solar multiple (SMT) with a guess value of thermal efficiency (assumed equal to 85% based on the thermal efficiency results obtained by Pasqualotto et al.32). However, it must be stressed that SM used in determining the required incident power on the receiver during the solar field sizing is, as aforementioned, only theoretical. This is because the actual value of SM depends on the actual design thermal efficiency of the receiver, which in turn depends, among others, on the type of PHC adopted. Therefore, the same solar field will have different actual SM values for Carbo ID 50 and sand. Aperture sizes for sizing each solar field are selected to meet a maximum flux of 2 MW m−2 (ref. 26) while keeping the average flux on the receiver equal to the conventional molten salts receiver (0.5 MW m−2 according to Blanco and Santigosa39). For the sake of simplicity, the receiver aspect ratio is assumed to be equal to one, as in the study by Pasqualotto et al.32
Parameter | Value | References |
---|---|---|
Site | Seville | 31 |
Latitude | 37.4° | |
Longitude | −5.9° | |
Design DNI (W m−2) | 900 | |
Design point | 21st June, Solar noon | |
Heliostat area (m2) | 16 (4 × 4 m) | 38 |
Heliostat focusing type | At slant | 38 |
Heliostat total reflected image error (mrad) | 3.07 | 31 |
Heliostat reflectivity (—) | 0.94 | 38 |
Design power to the HTM, SM = 1 (MW) | 1.73 | |
First guess thermal efficiency (—) | 0.85 |
Fig. 4 represents solar fields generated using the previously described approach and their design optical efficiency (not considering the intercept efficiency, i.e. neglecting the “spillage” losses of reflected solar radiation that does not hit the receiver aperture due to manufacturing errors) for the cases with SMT equal to 1, 5, and 9.
For each solar field and particle type, the optimal aperture size is selected to minimize in design conditions the trade-off between the spillage losses (complement of the intercept efficiency), and thermal losses, thus maximizing the design solar-to-thermal efficiency. Fig. 5 shows the values of solar-to-thermal efficiency for all considered sizes, SMT, and particle types as a function of the receiver aperture size. Table 4 illustrates the design parameters of the best performing aperture sizes.
Fig. 5 Solar-to-thermal efficiency of bauxite and sand particles for different SMT, as a function of the receiver aperture size. |
SMT = 1 | SMT = 3 | SMT = 5 | SMT = 7 | SMT = 9 | ||||||
---|---|---|---|---|---|---|---|---|---|---|
Carbo ID 50 | Sand | Carbo ID 50 | Sand | Carbo ID 50 | Sand | Carbo ID 50 | Sand | Carbo ID 50 | Sand | |
Tower height (m) | 25.8 | 31.4 | 36.9 | 42.4 | 47.9 | |||||
Number of heliostats | 191 | 591 | 996 | 1404 | 1801 | |||||
Cosine efficiency (%) | 88.75 | 86.10 | 85.55 | 85.29 | 85.27 | |||||
Optical efficiency (%) | 74.61 | 74.61 | 71.62 | 71.62 | 70.92 | 70.50 | 70.10 | 70.10 | 70.07 | 69.75 |
Thermal efficiency (%) | 82.71 | 72.29 | 85.45 | 74.58 | 85.68 | 75.21 | 86.28 | 75.21 | 86.32 | 75.63 |
Aperture size (m) | 1.8 | 1.8 | 3.0 | 3.0 | 3.9 | 3.8 | 4.5 | 4.5 | 5.1 | 5.0 |
Solar-to-thermal efficiency (%) | 61.71 | 53.93 | 61.20 | 53.42 | 60.76 | 53.02 | 60.48 | 52.72 | 60.48 | 52.75 |
Receiver power (MWth) | 1.702 | 1.488 | 5.262 | 4.601 | 8.858 | 7.743 | 12.460 | 10.906 | 16.072 | 14.062 |
Actual SM (—) | 0.984 | 0.860 | 3.041 | 2.659 | 5.120 | 4.476 | 7.202 | 6.304 | 9.290 | 8.129 |
From the results presented in Table 4, it is possible to draw the following conclusions:
• Optical efficiency decreases with SMT. This occurs because the cosine efficiency decreases as the average angle between the sun rays and the heliostat's surface normally increases more by adding new heliostats than it decreases by increasing the tower height.
• A considerable difference in the thermal efficiency (∼11%) between Carbo ID 50 and sand is obtained due to the lower solar absorptivity of sand with respect to Carbo ID 50 (55% versus 90.6%).
Taking everything into account, solar-to-thermal efficiency decreases slightly with SMT because an optical efficiency decrease has a higher impact than an increase in thermal efficiency. Additionally, Carbo ID 50 has ∼8% higher value of solar-to-thermal efficiency compared to sand because of the differences in thermal efficiency. The actual SM values are obtained by dividing the design thermal power of the receiver by the thermal power required by the pyrolizer. The obtained thermal efficiency values for Carbo ID 50 are slightly higher than those reported in the literature for similar average heat flux on the receiver:40 this is consistent with the lower receiver operating temperatures considered in this work (434–609 °C) compared to particle receivers in the literature, which are mainly designed for electricity generation purposes (operating at 550–750 °C).
After defining the design of the solar field for each SMT value, it is possible to compute the off-design performance of the receivers. The performance of each receiver is simulated by varying the DNI and thus the incident thermal power on the receiver from 0.2 to 1.15 of the nominal value while tuning the PHC flow rate to maintain a fixed temperature, which is equal to the design value. The results of this analysis are shown in Fig. 6a as thermal efficiency curves vs. thermal input to the receiver relative to the nominal one (Qin,rec/Qin,rec,nom). The thermal efficiency of both Carbo ID 50 and sand decreases more quickly as the incident thermal power decreases for smaller SMT. This is because the curtain absorptivity decreases more rapidly due to the lower mass flow rate per unit length for the smaller SMT receivers, which is ultimately because the receiver aperture size grows less than proportionally with the SMT (see Table 4). To better illustrate this, Fig. 6b shows the trends of the average absorptivity of the particle curtain and the mass flow rate per unit length of the receiver.
Solid PHC needs a lift to carry particles to the solar receiver from storage. The electrical power consumption for particle lift is calculated using eqn (8).41Plift is the power required to lift PHC mass flow rate ṁp in kg h−1 of solid particles to height Hlift in meters with an efficiency of ηlift that is considered as 0.8.41 An additional 10% of overall electricity consumption is added to consider the consumption of light and other instruments related to the CSP plant.
(8) |
The dashed lines in Fig. 3 illustrate the configuration when the plant is running in hybrid mode. This operating mode involves a reduction in the biochar output but allows for increasing the capacity factor of the pyrolysis plant, possibly improving the economic KPIs.
The Net Present Value (NPV) is set to zero at the end of plant life with a fixed discount rate i, as shown in eqn (9):
(9) |
• TCI represents the total capital investment for the pyrolysis and CSP plants (€),
• t is the year and includes a construction time of 3 years (t = −2, −1, 0),
• MFSPt represents the minimum selling price of the bio-oil (€ per t),
• Qt represents the quantity of bio-oil produced per year (t per year),
• Pby-prod,t gives the unit selling price of the byproducts produced by the plant/year (€ per t),
• Qy-prod,t represents the quantity of the byproducts produced by the plant/year (t per year),
• COP,VAR,t is the variable operational cost associated with the raw material and electricity per year (€ per year),
• COP,FIX,t is the fixed operational cost per year (€ per year),
• LR,t represents the interest paid on loan per year (€ per year),
• Tt represents the taxes paid for the year (€ per year).
• i is the discount rate (%).
According to data available in the literature,43 the biochar price may vary significantly, ranging between 0.37 and 6.42 € per kg, with average wholesale and retail prices of 1.89 € per kg and 2.83 € per kg, respectively. For the calculation of MFSP, the average wholesale price of 1.89 € per kg is considered, subject to a sensitivity analysis.
(10) |
(11) |
The total installed cost Cinst is estimated by multiplying the component cost by the installation factor (eqn (12)). The installation factor (ESI Table 11†) accounts for additional costs related to component installation, electrical network connections, control systems and civil work.
Cinst = finst × Ccomponent | (12) |
The total direct plant cost is estimated by summing up the total installed cost with the cost of additional materials, costs of additional components, and site development cost, which are assumed to be 4.5%, 4% and 10% of component purchase cost, respectively.29 The cost of land used to install the plant is also considered in the direct costs.
The overall cost of the pyrolizer, including the main reactor, cyclones, combustor, oil recovery system, compressors, and handling systems, is equal to 34.86 M€ for a capacity of 400 dry tons per day of treated biomass.29,45 The cost of biomass pretreatment includes the dryer, entire handling area, handling system, storage system and biomass grinder, which is scaled from the literature42 to our desired reference size and year. Similarly, the cost of utilities and auxiliary components, e.g. waste water treatment, cooling tower, products storage system, ventilation system, fire prevention system, and waste disposal system, is also considered (see ESI Table 11†). Apart from components costs, land cost was also considered assuming a pyrolysis plant footprint of 60000 m2 and an average land cost of 1.75 € per m2.46,47
Indirect costs, accounting for 60% of the total direct costs, are used in this work. These costs include construction charges, legal expenses, engineering start-up and contingencies costs.
Fixed Capital Investment (FCI) of the CSP section is computed by summing up the heliostat field, receiver, tower, TES, particle and particle elevator costs, computed with correlations and reference values adapted from the literature (see ESI Tables 13 and 14†) according to eqn (13):
FCI = Chel,field + Crec + Ctow + CTES,hot + CTES,cold + Cparticles + Celevator. | (13) |
The fixed operating costs include the wages of the employees and other expenses that are linked to the insurances, maintenance and taxes along with any other benefits and overheads. The annual costs related to the wages are estimated based on the data provided in ESI Table 16.†
Streams name/no. | Stream description | Flow rate (kg h−1) | T (°C) | P (bar) | Composition (% wt) | % C yield | |||||
---|---|---|---|---|---|---|---|---|---|---|---|
C | H | O | N | Ashes | PHC | ||||||
a Molar composition of pyro-gas: 52% CO, 31% CO2, 12% CH4, 4% C2+, 1% H2. b Molar composition of flue gases: 63% N2, 8.2% O2, 9% CO2, 19% H2O, 0.8% Ar. | |||||||||||
Input streams | |||||||||||
1 | Biomass (wet) | 2930 | 20 | 1.013 | 35.7 | 7.6 | 56.0 | 0.1 | 0.64 | — | — |
14 | Combustion air | 3829 | 32.2 | 1.013 | 0.01 | 0.2 | 24.8 | 75 | — | — | — |
25 | Air to dryer | 2094 | 32.2 | 1.013 | 0.01 | 0.2 | 24.8 | 75 | 0.00 | — | — |
PHC-MK | PHC make-up | 3.2 | 32.2 | — | — | — | — | — | — | 100 | — |
Output streams | |||||||||||
12 | Pyro-gasa | 135 | 72.03 | 1.59 | 39.6 | 2.3 | 58.1 | 0.0 | — | — | 5 |
Bio-oil | Bio-oil | 1732 | 54.4 | 1.24 | 41.6 | 7.8 | 50.6 | 0.0 | — | — | 69 |
17 | Flue gasesb | 6969 | 71.7 | 1.013 | 3.9 | 1.4 | 30.9 | 63.8 | — | — | 26 |
Ash | Ash to disposal | 22 | 608.4 | — | — | — | — | — | 85.6 | 14.4 | — |
Table 6 shows the energy balance and the main performance indicators of the conventional pyrolysis plant. Along with 10 MWLHV of input biomass, 0.1 MW of electricity is consumed to run compressors and heat rejection auxiliaries for the pyrolysis process. The chemical energies contained in bio-oil and pyro-gas are 7.57 MWLHV and 0.36 MWLHV, respectively. Heat losses from the plant units (e.g. pyrolyzer, combustor, compressor and filters) are equal to 2.18 MW, accounting for almost 22% of total input biomass LHV. Heat dissipation losses also occur during quenching at the dryer outlet and char gas and oil cooling processes. The energy conversion efficiency of 73.9% demonstrates that even a conventional pyrolysis process can retain a substantial portion of the input biomass energy in the form of produced bio-oil and pyro-gas. A carbon efficiency of 74.1% indicates that most of the carbon in the biomass is retained in useful products rather than being emitted as CO2. However, a significant portion is emitted into the environment due to the combustion of by-products. In general, the biomass conversion process is carbon neutral, so the net emission is still considered zero.
Parameter | Value |
---|---|
Biomass input (kWLHV) | 10000 |
Bio-oil production (kWLHV) | 7574 |
Pyro gas production (kWLHV) | 362 |
Heat rejected to ambient (kW) | 1986 |
Heat losses at combustor and pyrolyzer (kW) | 200 |
Compressor electric power (kW) | 118 |
Electric consumption for heat rejection (kW) | 15 |
Electric consumption for biomass pretreatment (kW) | 200 |
Other auxiliaries consumption (kW) | 34 |
Energy conversion efficiency, ηpyro plant (%) | 73.9% |
Carbon efficiency, εC (%) | 74.1% |
Emission to oil ratio, ETO (kgCO2 GJoil−1) | 36.5 |
Net emission to oil ratio, ETOnet (kgCO2 GJoil−1) | 0.0 |
Streams name/no. | Stream details | Flow (kg h−1) | T (°C) | P (bar) | Composition (% mass) | % C yield | |||||
---|---|---|---|---|---|---|---|---|---|---|---|
C | H | O | N | Ashes | PHC | ||||||
a Molar composition of flue gases: 58.8% N2, 12.8% O2, 4.8% CO2, 23% H2O, 0.7% Ar. | |||||||||||
Input streams | |||||||||||
Biomass | Biomass (wet) | 2930 | 20 | 1.013 | 35.7 | 7.6 | 56.0 | 0.1 | 0.64 | — | — |
14 | Combustion air | 3659 | 25 | 1.59 | 0.0 | 0.0 | 23.5 | 76.5 | — | — | — |
PHC-MK | PHC make-up | 3 | 25 | 1.22 | — | — | — | — | — | 100 | — |
Output streams | |||||||||||
Bio-oil | Boi-oil | 1732 | 54.4 | 1.24 | 41.6 | 7.8 | 50.6 | 0.0 | — | — | 69.0 |
Biochar | BioChar | 246 | 25 | 1.22 | 83.1 | 1.7 | 6.6 | 1.4 | 7.1 | — | 19.6 |
Sludge | Sludge | 59 | 54.4 | 1.24 | 38.3 | 7.2 | 46.6 | — | 2.4 | 5.4 | 2.2 |
17 | Flue gasesa | 4555 | 71.7 | 1.013 | 2.1 | 1.7 | 34.7 | 61.5 | — | — | 9.3 |
Table 8 illustrates the energy balance and the main performance indicators of the CSP-based pyrolysis plant. In solar-based pyrolysis, the CSP block provides 1730 kWth to the pyrolyzer by transferring thermal power through the PHC. Electrical power consumption, biomass, oil LHV and heat dissipation are similar to the conventional pyrolysis plant. Char and sludge chemical power are additional outputs that lead to an increase in energy efficiency ranging from 73.9% to 78.6%, indicating a more effective use of the input biomass energy. The additional energy input from concentrated solar power (CSP) helps to achieve this gain in efficiency. Carbon efficiency improves from 74.1% to 90.7%, meaning that a greater proportion of biomass carbon is retained in useful products (bio-oil and biochar). This higher carbon efficiency reflects the lower carbon losses and better utilization of feedstock biomass carbon.
Parameter | Value |
---|---|
Biomass input (kWLHV) | 10000 |
Concentrated solar power (kW) | 1730 |
Bio-oil production (kWLHV) | 7574 |
Char production (kWLHV) | 1998 |
Sludge production (kWLHV) | 238 |
Heat rejected to ambient (kW) | 1847 |
Heat losses at combustor and pyrolyzer (kW) | 200 |
Compressor electric power (kW) | 124 |
Electric consumption for heat rejection (kW) | 15 |
Electric consumption for biomass pretreatment (kW) | 200 |
Other auxiliaries consumption (kW) | 34 |
Electric consumption for PHC lift (kW) [SMT = 1] | 5.1 |
Energy conversion efficiency, ηpyro plant (%) | 78.6% |
Carbon efficiency, εC (%) | 90.7% |
Emission to oil ratio, ETO (kgCO2 GJoil−1) | 13.06 |
Net emission to oil ratio, ETOnet (kgCO2 GJoil−1) | −27.5 |
The emission to oil ratio significantly decreases from 36.5 to 13.1 kgCO2 GJoil−1, demonstrating a reduction in the carbon footprint of the produced bio-oil. Moreover, the net negative emissions of −27.5 kgCO2 GJoil−1 indicate that the process can act as a carbon sink, removing more CO2 from the atmosphere. In terms of energy content, the production of char and sludge contributes 1998 kWLHV and 238 kWLHV, respectively.
In the hybrid pyrolysis plant, the stream flows are the same as in the CSP-based case, provided hot PHC is available from the receiver or form the high-temperature PHC storage. When hot PHC is not available, sludge and almost 75% char are sent to the combustor to provide an alternative heat source for the PHC, as shown in Fig. 3.
Components | Conventional M€ | CSP based M€ | Hybrid M€ |
---|---|---|---|
a The values of CSP plant CAPEX are at the selected size of CSP plant based on MFSP. b This revenue is based on the MFSP of oil at the selected plant size. | |||
Pyrolysis plant CAPEX (TCI) | 20.25 | 19.43 | 21.11 |
Plant (pyrolizer and oil recovery system) | 8.53 | 8.53 | 8.53 |
Solid char combustor | 0.94 | — | 0.94 |
Gas combustor | — | 0.48 | 0.48 |
Biomass pretreatment | 1.24 | 1.24 | 1.24 |
Utilities and auxiliaries | 0.48 | 0.48 | 0.48 |
Material stock up | 0.17 | 0.16 | 0.17 |
Additional components | 0.19 | 0.18 | 0.19 |
Site development | 0.42 | 0.40 | 0.43 |
Land cost for pyrolysis plant | 0.10 | 0.10 | 0.10 |
Indirect cost | 7.23 | 6.94 | 7.54 |
Working capital | 0.96 | 0.93 | 1.01 |
CSP plant CAPEX | — | 6.51 | 4.82 |
Heliostat field | — | 3.32 | 2.35 |
Receiver | — | 0.69 | 0.49 |
Tower | — | 0.20 | 0.15 |
Thermal energy storage (TES) | — | 0.72 | 0.68 |
PHC particle | — | 0.04 | 0.04 |
PHC elevator cost | — | 0.16 | 0.09 |
Land cost for CSP plant | — | 0.3 | 0.2 |
Indirect cost & contingencies | — | 1.07 | 0.79 |
Total CAPEX | 20.25 | 25.94 | 25.93 |
Pyrolysis plant total O&M/year | 2.99 | 2.86 | 3.02 |
Variable operational cost (biomass)/year | 0.34 | 0.28 | 0.34 |
Variable operational cost (other materials & disposal)/year | 0.23 | 0.19 | 0.23 |
Fixed operational cost/year | 2.42 | 2.39 | 2.45 |
CSP plant O&M/year | — | 0.13 | 0.1 |
Total O&M/year | 2.99 | 3.0 | 3.1 |
Biochar revenues/year | 0 | 3.36 | 3.31 |
Bio oil revenues/yearb | 6.58 | 4.29 | 4.46 |
Total revenues/year | 6.58 | 7.65 | 7.78 |
Material flows (kt per year) | Conventional | CSP based | Hybrid |
Biomass consumption (kt per year) | 25.7 | 20.6 | 25.7 |
Bio-oil production (kt per year) | 15.2 | 12.8 | 15.2 |
Sludge production (kt per year) | 0 | 0.44 | 0.38 |
Pyro-gas (kt per year) | 1.2 | 0 | 0 |
Char (kt per year) | 0 | 1.8 | 1.8 |
CO2 emissions (kt per year) | 8.7 | 2.6 | 4.5 |
ETO ratio (kgCO2 GJoil−1) | 36.5 | 13.1 | 19.0 |
Net ETO ratio (kgCO2 GJoil−1) | 0 | −27.5 | −22.3 |
MFSP (€ per GJ) | 27.53 | 21.79 | 18.68 |
Fig. 7 MFSP of the CSP-based pyrolysis plant (a) and the hybrid pyrolysis plant (b) at various SMT and storage sizes. |
The total CAPEX of the CSP unit of the selected plant (SMT 7, sand, 24 h TES) is 6.51 M€ (Table 9), which corresponds to about 25% of the total plant cost. The total operational and maintenance cost, including the electricity required to run the lift and other components, is 0.13 M€ per year (Table 9). The total CAPEX and O&M costs of the pyrolyzer for this case are slightly lower than for the conventional case because the combustor size is significantly smaller than the conventional one. For this case, the CAPEX and O&M costs of pyrolyzer are 19.43 M€ and 2.86 M€ per year, respectively (Table 9), while feedstock, consumable and utility costs (e.g. biomass, water, electricity, and disposal cost) depend on the capacity factor, which is around 82% at the selected size of CSP. The calculated MFSP is 0.34 € per kgoil or 21.79 € per GJ.
Despite lower bio-oil production in CSP-based plants (12.8 kt per year vs. 15.2 kt per year for the conventional plant), the process is economically competitive thanks to the additional revenues from biochar production. The CSP-based and hybrid plants generate significant additional revenue from biochar (1.07 M€ per year and 1.2 M€ per year, respectively), thereby offsetting the higher capital and operational costs. This additional income contributes to the economic feasibility of the solar-based pyrolysis systems.
Fig. 11a shows the capacity factor of the pyrolysis plant and the monthly product yield. The plant capacity factor is limited to below 71% from November to February, while it peaks to above 95% in July and August. Bio-oil and biochar productivity as well as CO2 emissions reflect the plant capacity factor trend.
Fig. 11 Monthly distribution of products, emissions and availability in (a) CSP-based pyrolysis case and (b) hybrid pyrolysis at optimized scenarios. |
Fig. 11b shows the monthly production and emission trends of the hybrid plant. As the hybrid plant runs stably throughout the year, the capacity factor of the pyrolysis unit is always 100%, and bio-oil yield depends only on the number of days in each month. In this plant, solar radiation affects biochar production and CO2 emissions. In the winter months, when the plant relies more on biochar combustion, biochar yield is reduced and CO2 emissions increase. The sizing of the solar field and the TES significantly affects the product yield and the carbon balance. ESI Fig. 3a and b† show the CO2 emission trends with respect to SMT and storage size of CSP-based and hybrid pyrolysis, respectively.
Fig. 12 Sensitivity analysis of MFSP under rate of change of different variables for (a) conventional pyrolysis plant, (b) CSP-based pyrolysis plant and (c) hybrid CSP pyrolysis plant. |
In the CSP-based plant, the availability varies based on the SMT and storage size and its analysis is shown in Fig. 7a. The analysis in Fig. 12b refers to the selected optimized case. The DR and Pyrolyzer CAPEX have a significant impact, as in the conventional pyrolysis case. For example, a 60% decrease in DR, taxes and pyrolyzer Capex can decrease MFSP by almost 40% and vice versa. The fixed operational cost is the second most important factor because a 60% change in the rate can cause almost a 37% change in MFSP. The biochar selling price also significantly affects the MFSP. A variation between −80 and +80% (i.e. from 0.38 to 3.4 € per kg) involves a reduction in the MFSP from about 35 to 8.6 € per GJbio-oil. The other parameters, including the CSP plant cost, have a lesser impact on the MFSP (Fig. 12b).
Fig. 12c shows the same analysis for hybrid cases with similar trends. Decreasing availability by 40% increases the MFSP by almost 106%. Due to the lower production of biochar, its selling price affects the MFSP less than in the CSP-based plant: the variation between −80 and +80% leads to a reduction from 29.4 to 7.93 € per GJbio-oil.
The main additional source of revenues for CSP-based and hybrid cases is biochar. Fig. 13 shows that biochar selling prices of at least 1.24 and 0.64 € per kg (corresponding to carbon credits that are equal to 407 € per tCO2 and 210 € per tCO2, respectively) make CSP-based and hybrid pyrolysis competitive with conventional pyrolysis.
Fig. 13 MFSP variation with respect to biochar selling price for CSP-based and hybrid pyrolysis plants. |
This study relied on a simplified pyrolysis model with a constant product yield. Future research activity should focus on the development and implementation of detailed validated pyrolysis reactor models, including chemical kinetics, to evaluate the impact of the operational parameters (temperature, pressure, feedstock, residence time and heating rate) on the product yield and the solar section design. Additionally, PHC-char separation needs to be experimentally assessed and numerically modeled, as the separation efficiency affects the quality of biochar for end use as well as the heat transfer properties of the PHC in the solar receiver. Moreover, the performance of different PHCs through modelling and experimental activities needs to be assessed.
➢ Solar-based pyrolysis can achieve over 90% carbon efficiency, resulting from about 70% of the inlet biogenic carbon retained in the bio-oil and about 20% of the carbon in the biochar. Based on yearly simulations, the hybrid plant results in a carbon efficiency of 83%. The conventional pyrolysis process achieves 74% carbon efficiency. Owing to the carbon stored in biochar, the solar-based and hybrid plants achieve net negative emissions of −27.5. and −22.3 kgCO2 GJoil−1, respectively.
➢ For the hybrid plant, with an assumed biochar selling price of 1.89 € per kg, a minimum bio-oil selling price of 18.68 € per GJLHV was obtained, vs. 27.53 € per GJLHV of the reference pyrolysis process. For the CSP-based plant, the estimated MFSP is equal to 21.79 € per GJLHV. The cost of bio-oil produced in the CSP-based plant is higher than in the hybrid case despite the higher char yield due to the higher impact of Capex caused by the lower yearly capacity factor (82% vs. 100% of the hybrid case). Therefore, an economic-environmental trade-off exists between CSP-based and hybrid configurations, driven by minimum fuel selling price (i.e. fuel production cost) and overall process carbon efficiency.
➢ The sensitivity analysis showed that the variation of the discount rate, plant availability, biomass cost, biochar selling price and pyrolyzer CAPEX exhibits the highest impact on the MFSP. However, the type of particle heat carrier (sand or Carbo ID50 PHC) has a minor impact on the MFSP.
➢ The breakeven biochar selling prices for CSP-based and hybrid cases that can make them economically competitive with conventional pyrolysis are 1.24 € per kg and 0.64 € per kg, which correspond to carbon credit values of 407 € per tCO2 and 210 € per tCO2, respectively.
Capex | Capital expenditure |
CSP | Concentrated solar power |
DR | Discount rate |
DSR | Dissipated solar radiation |
ETO | Emission to oil ratio |
FBR | Fluidized bed reactor |
FCI | Fixed capital investment |
HHV | Higher heating value |
HTF | Heat transfer fluid |
HTM | Heat transfer material |
LHV | Lower heating value |
MFSP | Minimum fuel sale price |
NREL | National Renewable Energy Laboratories |
PC | Product combustion |
PHC | Particle heat carriers |
SMT | Theoretical solar multiple |
TCI | Total capital investment |
TDC | Total direct cost |
TES | Thermal energy storage |
TIC | Total installed cost |
TPC | Total purchase cost |
VBA | Visual basic analysis |
A h | Total heliostats area [m2] |
C inst | Installation cost [€] |
C component | Component cost [€] |
C o | Cost of reference component size [€] |
C x | Cost of selected component size [€] |
C hel,field | Cost of heliostat & field [€] |
C rec | Cost of receiver [€] |
C tow | Cost of tower [€] |
C TES,hot | Cost of hot thermal energy storage [€] |
C TES,cold | Cost of cold thermal energy storage [€] |
C particles | Cost of PHC particles [€] |
C elevator | Cost of PHC elevator [€] |
C OP,VAR,t | Variable operational cost per year [€ per year] |
C OP,FIX,t | Fixed operational cost per year [€ per year] |
DNI | Direct normal irradiance [W m−2] |
E | CO2 emissions [kg s−1] |
ε C | Carbon efficiency [—] |
f inst | Installation factor [—] |
g | Gravitational acceleration [m s−2] |
H lift | Height of PHC lift [m] |
i | Discount rate [%] |
LHVbiom | LHV of biomass [MJ kg−1] |
LHVprod,i | LHV of the ith product [MJ kg−1] |
L R,t | Loan with the inclusion of interest paid per year [€ per year] |
ṁ biom | Mass flow rate of biomass [kg s−1] |
ṁ prod,i | Mass flow rate of the ith product [kg s−1] |
η pyro plant | Pyrolysis plant energy conversion efficiency - |
m p | PHC mass flow rate [kg s−1] |
n | Scaling factor [—] |
η el,ref | Electric efficiency [—] |
η sol–th | Solar-to-thermal efficiency [—] |
η opt | Optical efficiency [—] |
η th,rec | Thermal efficiency [—] |
η th,use | Thermal energy usage efficiency [—] |
η lift | Lift efficiency [—] |
Oil | Oil production [kg s−1] |
P Aux | Electric power used by auxiliaries [W] |
P by-prod,t | Selling price of by-products per year [€ per year] |
P lift | Electric power used by PHC lift [W] |
def | Total solar power loss by defocusing [W] |
PHC,pyro | Thermal power provided to pyrolyzer via PHC [W] |
PHC,rec | Thermal power delivered to PHC [W] |
rec | Solar power incident on the receiver [W] |
Q t | Quantity of bio-oil produced per year [t per year] |
Q y-prod,t | Quantity of the byproducts produced per year [t per year] |
S 0 | Reference component size [t per day] |
S X | Selected component size [t per day] |
T t | Taxes paid per year [€ per year] |
y C,biom | Carbon content of biomass [%] |
y C,prod,i | Carbon content in the ith product [%] |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4se00450g |
This journal is © The Royal Society of Chemistry 2024 |