The Energy Potential, Environmental Impact, and Occupational Health and Safety Potential of Biogas Obtained from Filter Cake in Artisanal Panela Production

1. Introduction

Globally, around 80% of Saccharum officinarum crops are used for sugar production and 20% for ethanol [1]. In 2020, around 1,869.7 million tons of sugarcane were harvested, and Brazil was the main producer, growing 757.1 million tons, followed by India, China, Pakistan, and Thailand [2]. Brazilian productivity can achieve up to 72 t/ha, with the southeast region accounting for 60% of national production [3]. The sugar industry generates approximately 279 million metric tons of waste each year, including bagasse and filter cake [4]. The latter, also known as cachaza, is produced during juice clarification and has a high moisture content (≈75%) as well as compounds such as cellulose, hemicellulose, lignin, sucrose, and minerals [5]. The composition of filter cake varies depending on the type of sugarcane, the soil, the climate, and the efficiency of the milling process [6].

The province of Pastaza, in the Ecuadorian Amazon, has an estimated sugarcane production potential of 4,500 hectares, with an average yield of 55 t/ha, mainly for the production of derivatives such as panela in blocks, granulated panela, aguardiente (a type of sugarcane liquor), and fresh sugarcane sticks [7]. The local sugarcane industry is made up of more than 110 small factories, which, alongside their products, generate a considerable amount of organic waste that accumulates around these facilities [8]. This activity depends on the use of solid biomass wood and bagasse as a source of thermal energy, which leads to a high consumption of forest resources and the emission of incomplete combustion gases. Therefore, the valorization of filter cake as an energy substrate is a strategic alternative to reduce deforestation, decrease the use of traditional biomass, and strengthen the sustainability of artisanal sugarcane systems. However, there is no research available that integrates stoichiometric simulation with TRACI (Tool for the Reduction and Assessment of Chemical and Other Environmental Impacts) and social life cycle assessment (S-LCA) tools to demonstrate the energy potential of biogas obtained from filter cake and its environmental and social benefits in artisanal panel systems.

The processing of sugarcane using rudimentary technologies continues to be common practice within small factories in the province of Pastaza. In several sugar mills (paneleras), filter cake is considered a waste product due to its high organic load, with chemical oxygen demand (COD) values of up to ≈ 90,000 mg/L and biochemical oxygen demand (BOD₅) values of up to ≈ 40,000 mg/L [9]. The final disposal of filter cake is usually carried out without prior treatment, thus intensifying the contamination of soils, water bodies, and areas surrounding the trapiches [10]. [11] reported that 71.4% of sugar mills do not apply any kind of environmental management, while [12] identified that 89.5% of environmental impacts come from processes related to the management of waste, such as filter cake, and other stages of the panela production process and cane cutting.

Anaerobic digestion represents an appropriate technological option for the treatment of agro-industrial waste with high moisture content, such as filter cake. This technology, in addition to being economically accessible, allows for the production of biogas as a renewable energy source and a stabilized by-product (digestate) with potential use as fertilizer [13]. Previous records of sugarcane waste digestion have focused on kinetic estimates or laboratory conditions, without incorporating a comprehensive approach that links biogas production with environmental sustainability indicators. This gap limits the possibility of comparing the energy efficiency of biogas with traditional fuels and makes it difficult to quantify its real environmental benefits.

The yield of anaerobic digestion is directly dependent upon the type of substrate used [14]. One study found that the combination of filter cake with agricultural waste achieved a production of 120 m³/t of methane (CH₄), whereas the exclusive use of filter cake yielded outputs ranging from 73.7 to 129.5 mL/g of CH₄ [15]. Similarly, [16] reported the generation of 0.360 m³ of biogas from 10 kg of cattle manure, and [17] achieved an 81% biogas yield using vinasses as substrate. The quality of biogas is associated with the proportion of methane present, which is the component with the highest calorific value. In simple applications, such as stoves, heaters, or generators, the removal of carbon dioxide (CO₂) is not required. However, for more demanding uses like vehicle fuel, it is necessary to purify the biogas to maximize its energy efficiency [18].

This study proposes a comprehensive approach applied to artisanal panela production, combining stoichiometric simulation in CHEMCAD [19,20] with the TRACI and S-LCA tools. This integration enables a simultaneous analysis of the energy yield of biogas obtained from filter cake and the environmental impacts associated with its use. It also incorporates the social dimension linked to the operation of the biodigester, offering a more complete view of the production system. In this way, a technical basis is established that contributes to optimizing the use of agro-industrial waste and promoting sustainable practices in artisanal panela systems.

The analysis of environmental impact using the TRACI tool allows for the evaluation of emissions associated with thermochemical and biological processes by considering compounds such as CH₄, CO₂, H₂S, and NOₓ [21]. This methodology classifies impacts into categories such as climate change, soil acidification, tropospheric ozone formation, eutrophication, human toxicity, photochemical smog, and ozone depletion [22], providing a quantitative basis for comparing energy alternatives to the use of fossil fuels or solid biomass. The choice of TRACI 2.2 over methodologies such as life cycle assessment (LCA), GREENSCOPE (Gauging Reaction Effectiveness for the Environmental Sustainability of Chemistries with a Multi-Objective Process Evaluator), or WAR (Waste Reduction Algorithm) is justified by its ability to estimate environmental impacts directly from simulated process flows, without requiring complete life cycle inventories. It sets itself apart from LCA, which requires a database, or GREENSCOPE and WAR, which focus on the efficiency and safety of large-scale chemical plants.

In addition to the environmental component, the sustainability of the sugarcane industry includes social and occupational dimensions. In this study, S-LCA was considered to identify occupational health and safety risks for personnel involved in the operation of the biodigester. The Occupational Health and Safety Potential (OHSP) indicator allowed for the quantification of hours potentially lost due to workplace accidents. Therefore, this study methodologically integrated CHEMCAD simulation, S-LCA, and the TRACI model and applied it to the case of an artisanal sugarcane mill in Pastaza. This enabled the simultaneous estimation of biogas yield, solid biomass savings, and associated environmental impact categories, providing a technical basis for the energy transition of the Ecuadorian sugarcane sector.

Based on the above, this study aimed to determine the energy potential, environmental impact, and OHSP of biogas obtained from the bagasse generated in the artisanal production of panela in Pastaza, using simulations in CHEMCAD, TRACI, and S-LCA software to estimate the sustainable performance of the panela system.

2. Materials and Methods

This research followed a sequential methodological procedure aimed at assessing the energy use of filter cake and its effect on the panela system. The first stage involved collecting primary and secondary information from the plant in Pastaza, including data on milling capacity, filter cake generation, fuel consumption, and operating conditions. Based on this, the filter cake was characterized, and representative compositions of moisture, organic matter, and ash were defined, which served as direct input for modeling anaerobic digestion and defining input and output streams in the diagram of the process (Figure 1).

The final phase integrated modeling and simulation of the process in CHEMCAD with two levels of evaluation. First, the calorific value of the biogas and the amount of wood and bagasse that can be replaced in the boiler were calculated, consistent with actual plant data. Second, TRACI 2.2 was applied together with the OHSP indicator within the S-LCA approach to estimate environmental impacts and occupational risks.

2.1. Data Collection

The information used to characterize the filter cake was obtained from scientific sources indexed in international databases (Scopus, SciELO, Elsevier, Springer, and Wiley), using keywords such as filter cake, biogas, and anaerobic digestion. Recent (≤ 5 years) average compositions were prioritized, including moisture, ash, fiber, protein, lipids, and sugars. The information collected was compared with local data on sugarcane yield and panela production in the province of Pastaza [7,23]. This allowed the average compositions obtained from the literature to be adjusted, and it ensured their representativeness in the Amazonian context.

2.2. Characterization of the Filter Cake

The filter cake composition was analyzed based on information obtained from recent scientific literature, prioritizing experimental reference values. Data on crude protein, lipids, ash, sucrose, total fiber, and their fractions were incorporated, along with profiles of amino acids, fatty acids, minerals, and bioactive compounds reported by [5,24,25,26].

2.3. Modeling of the Process

Different boundary conditions for anaerobic digestion were analyzed by considering the upper and lower limits of the main components of the filter cake (Table 1). The purpose was to establish a representative range of compositions, including a base profile, in order to evaluate the system’s response to possible variations in the raw material.

Simulation scenarios were defined based on different feed flow levels (Table 2) by adjusting the percentage of available bagasse according to the theoretical value for the province of Pastaza. These combinations allow simulations to be run that show the behavior of the system under variable load conditions, with an emphasis on the amount of biogas in terms of CH₄, CO₂, and H₂S.

2.4. Simulation of the Process

The simulation was performed in CHEMCAD 6.3.1 by modeling the anaerobic digestion of filter cake using mass balances. Daily capacity was estimated based on sugarcane data in Pastaza [7] and by considering usable area and land set aside for panela. The operating conditions were a mesophilic temperature of 25–35 °C and pH of 6.8–7.2 [27]. These parameters ensure microbiological stability in anaerobic processes.

2.4.1. Estimation of Filter Cake Production in Pastaza

The estimate was made according to the availability of sugarcane and the percentage of filter cake extracted. Sugarcane production and the percentage of usable sugarcane crops were established according to [7], while the estimated cane yield and percentage used for panela production were set according to [28], and the percentage of filter cake was determined as proposed by [29]. Based on these parameters and considering an annual operating period of 365 days, a tentative daily filter cake production of 9,760 kg d⁻¹ was estimated, and this value was used as a reference for defining feed flows and process simulation scenarios.

2.5. Technical Evaluation of Biogas

The Low Heating Value (LHV) was obtained by extrapolating the values published in [30]. The LHV of wood (Eucalyptus globulus) was obtained from [31] and that of bagasse from Equation 1 proposed by [32].

LHV= [(18260 – 207.63 (moisture % sample) – 182.6 (ash % sample) – 31.14 (brix % sample)] kJ.kg^{- 1}

(1)

The energy balance equation $Q_{biogas}=Q_{wood}$ and $Q_{biogas}=Q_{bagasse}$ was used to estimate the amount of wood and bagasse that can be replaced by the biogas generated, considering the heating values of the three fuels evaluated.

The energy savings obtained by replacing solid biomass with biogas were compared with the average aboveground biomass of the Amazon rainforest to establish an ecological equivalence. To do this, a value of 174 Mg ha⁻¹ was used as the average dry biomass density reported by [33], and an average individual biomass of 0.70 t tree⁻¹. Using these references, the proportion of wood replaced that would correspond to the preserved forest area and number of equivalent trees was calculated, thus linking the use of biogas to the potential for rainforest preservation in the Amazon.

The results were compared with recent studies on the anaerobic digestion of filter cake and sugarcane residues carried out by [15,34,35]. The reference datasets included the organic matter composition, moisture content, and lignin content of filter cake, along with the experimental biogas and methane yields reported in those studies. The comparison covered total biogas volume, methane (CH₄) fraction, and hydrogen sulfide (H₂S) concentration, and there was a general alignment between local results and international experimental ranges.

2.6. Environmental Assessment with TRACI 2.2

In estimating the environmental impact using TRACI 2.2, only the gaseous flows generated during anaerobic digestion were considered, as they constitute the main source of direct emissions associated with the energy use of biogas. Liquid and solid flows—corresponding to digestate and process sediments—were not included because their impact requires further treatment depending on the type of disposal or agricultural use, which exceeds the system boundaries defined for this assessment. The quantitative analysis focused on emissions of CH₄, CO₂, H₂S, and NOₓ, selected for their significant contribution to the categories of the impact of climate change (kg CO₂ eq), acidification (kg SO₂ eq), and tropospheric ozone formation (kg O₃ eq), in accordance with the methodological guidelines of the TRACI 2.2 model.

The substances considered for each category were CH₄ and CO₂ for climate change, NOₓ and H₂S for acidification, and CH₄ and NOₓ for photochemical smog. The composition of the biogas was obtained from experimental results. For natural gas (≈98% CH₄), liquefied petroleum gas (LPG), and propane, estimated emission values for complete combustion were used based on technical data [36,37,38]. The impact calculation consisted of multiplying the estimated mass of each substance by its corresponding TRACI characterization factor.

2.7. Occupational Health and Safety Potential (OHSP)

The study incorporated the occupational health and safety component through the use of the OHSP indicator, developed by [39], as a midpoint indicator within the S-LCA. This facilitated the quantification of hours potentially lost due to occupational accidents in each phase of the biogas production process (Equation 2). The reference values came from the [40] database in conjunction with the values from the study panel. The mathematical model is described by the following relationship:

$$OHS{P}_{n,a}=C{F}_{OHS,n,a}\times W{H}_n$$

(2)

where:

$C{F}_{OHS,n,a}$ is the occupational characterization factor, calculated as:

$$\mathrm{C}{\mathrm{F}}_{\mathrm{O}\mathrm{H}\mathrm{S},\mathrm{n},\mathrm{a}}={\displaystyle \frac{24\times \left(\mathrm{Accidents}\mathrm{resulting}\mathrm{in}\mathrm{lost}\mathrm{days}\right)}{\mathrm{Hours}\mathrm{worked}\mathrm{per}\mathrm{year}}}$$

(3)

$W{H}_n$ corresponds to the total work hours in stage n of the production process.

The OHSP indicators module reports that the agro-industrial sector has a global average rate of 2.5 accidents resulting in lost days per 100 workers, with an average annual working time of 2,000 hours per employee and an average duration of 4 days lost per accident. These values were applied as reference parameters, adjusted proportionally to the actual hours of the simulated process.

The filter cake-based biogas system was divided into four operational stages: receiving and preparing the filter cake, loading and feeding the biodigester, operating and controlling anaerobic digestion, and cleaning and discharging the digestate. At each stage, the main exposure factors were identified and the OHSP indicator was estimated, taking into account the characteristics of the work, the annual duration, and the nature of the existing risks.

3. Results

3.1. Bibliographic Characterization

Table 3 shows the compounds present in filter cake, including macronutrients, essential and non-essential amino acids, saturated and unsaturated fatty acids, minerals, and phytochemicals with functional potential. This characterization allows the identification of components suitable for agro-industrial use, such as structural fiber, phytosterols, and residual sugars.

3.2. Mathematical Modeling

Mathematical models for predicting biogas production were developed based on the general mass balance equation and the chemical composition of the waste. The modeling considered the main organic fractions present in filter cake, allowing the formulation of stoichiometric reactions that describe the conversion of organic matter during anaerobic digestion (Table 4).

3.3. Modeling of the Process

3.3.1. Simulation of the Anaerobic Digestion Process

The simulated process starts with the storage of filter cake (Unit 1 in Figure 2), which is subsequently fed into the anaerobic digester operated under mesophilic conditions at approximately 25 °C. This operating temperature was selected to maintain stable microbial activity during anaerobic digestion (Singh et al., 2020). After digestion, two main output streams are generated: digestate and biogas. The digestate is discharged through Stream 4, while the biogas is conveyed through Stream 5 toward a purification unit.

In the purification stage, the biogas comes into contact with iron(III) oxide supplied as nanoparticles or granular hematite. This interaction promotes the removal of hydrogen sulfide from the gas phase. The solid products formed during this stage are discharged through Stream 9, while the purified biogas exits the system as the final gaseous product (Figure 2).

3.3.2. Estimation of Bagasse Production in the Sugar Mill

The estimated amount of bagasse generated in the province of Pastaza is approximately 9,760 kg·d⁻¹, considering an extraction rate of 4% per ton of processed sugarcane [28]. Based on a processing capacity of 10 t·d⁻¹, a total of 400 kg·d⁻¹ of filter cake was estimated as available substrate for biogas production in the simulation. To evaluate system performance under variable operating conditions, the simulation was performed using different feedstock capacities (Table 5).

3.4. Process Simulation

3.4.1. Base Production Case According to the Literature Review

According to [35], 120–155 m³ of biogas can be obtained from one ton of processed filter cake, i.e., a yield of 12–15.5%. Under the base production scenario, a feed flow of 1,000 kg·d⁻¹ (1 t·d⁻¹) of filter cake generated 172.98 m³·d⁻¹ of biogas, corresponding to a yield of 17.2%. The simulated biogas consisted mainly of methane and carbon dioxide, while ammonia and hydrogen sulfide were produced in lower proportions (Figure 3 A–B).

3.4.2. Production Flow Case in Pastaza

Figure 3A shows that biogas production increased as the feed flow increased, with the gas mixture mainly composed of CH₄ and CO₂. Figure 3B indicates the presence of NH₃ and SH₄ in lower proportions. The average methane content obtained was 39.6%, while carbon dioxide reached 58.3%. Ammonia represented 2.03% of the gas mixture, whereas hydrogen sulfide remained at a low level of 0.04%.

3.4.3. Boundary Conditions for the Proposed Model

Figure 4 shows the relationship between filter cake composition, feed flow, and biogas production. Biogas output varied according to both composition and production flow. Among the evaluated scenarios, Composition 6 achieved the highest biogas production, reaching 1,736.40 m³·d⁻¹ at a feed flow of 9,760 kg·d⁻¹. An overall increase in biogas production was observed as the feed flow increased, particularly for Compositions 6, 7, and 8.

Hydrogen sulfide remained stable and below 0.12% in all formulations, as shown in Figure 5. This indicates a low concentration of sulfur compounds in the gas phase and a consistent behavior of H₂S across the evaluated conditions.

3.5. Technical Evaluation of Biogas

The biogas obtained from filter cake showed an average calorific value of 13,900 kJ·m⁻³. In the highest production scenario evaluated, with a biogas flow of 1,288.62 m³·day⁻¹, an estimated saving of 1,242.88 kg·day⁻¹ of wood was obtained, equivalent to 1.24 t·day⁻¹ of solid biomass. In Composition 6, with a methane content of 40.69%, biogas production values of 71.16 and 1,736.4 m³·day⁻¹ were obtained at filter cake feed flows of 400 and 9,760 kg·day⁻¹, respectively. These production levels corresponded to wood savings of 70.86 and 1,728.99 kg·day⁻¹, as well as bagasse savings of 135.91 and 3,316.43 kg·day⁻¹ (Table 6).

3.6. Comparison of the Environmental Impact of Gaseous Fuels Using TRACI

The environmental impact assessment using TRACI factors compared four gaseous fuels: biogas, methane (natural gas), LPG, and propane, considering the stages from production to energy use. Biogas was produced through anaerobic digestion of filter cake under mesophilic conditions and purified using iron(III) oxide, resulting in a gas with a methane content of 40.7% (Table 7). Methane, LPG, and propane were considered based on their conventional industrial production routes from fossil resources.

The results show differences among fuels in the impact categories of climate change, acidification, and photochemical ozone formation. Biogas presented the highest climate change impact, reaching 17.40 kg CO₂ eq·m⁻³. Methane exhibited a lower value of 2.875 kg CO₂ eq·m⁻³, while LPG and propane showed intermediate values of 3.07 and 3.05 kg CO₂ eq·m⁻³, respectively.

In terms of acidification, biogas reached 0.00516 kg SO₂ eq·m⁻³, followed by LPG (0.00280 kg SO₂ eq·m⁻³), propane (0.00210 kg SO₂ eq·m⁻³), and methane (0.00140 kg SO₂ eq·m⁻³). For photochemical ozone formation, LPG and propane presented the highest impacts, with values of 0.09919 and 0.07440 kg O₃ eq·m⁻³, respectively, while biogas and methane reached 0.05894 and 0.04965 kg O₃ eq·m⁻³ (Table 7).

3.7. Occupational Assessment and Risk Analysis for Biogas Production

The results of the OHSP indicator calculation are summarized in Table 8. The stages of receiving and grinding the filter cake and cleaning and discharging the digestate showed the highest OHSP values, reaching 0.08 and 0.096, respectively. These stages also recorded higher numbers of accidents and lost workdays, with two accidents per year and three lost days in both cases.

Lower OHSP values were obtained for loading the biodigester (0.075) and operating the biodigester (0.048). These stages involved fewer accidents and were characterized by more controlled working conditions, with reduced direct exposure to waste and gaseous emissions. The numerical values for hours worked, accidents, days lost, and OHSP indicators for each process stage are presented in Table 8.

4. Discussion

The diversity of compounds identified in filter cake indicates the complexity of this by-product and its potential for agro-industrial use, including structural fiber, phytosterols, and residual sugars. This variability depends on crop management, sugarcane variety, and technological practices applied during processing [6].

The waste contains a wide variety of organic compounds, including polysaccharides and oligosaccharides, amino acids, long-chain fatty acids, alcohols, and phytosterols. These components determine the structure of the stoichiometric equations and lead mainly to the formation of CH₄ and CO₂ as final products, in agreement with previously reported mass balance-based models for anaerobic digestion [42].

The simulated configuration enables the optimization of substrate proportions and operational conditions to enhance methane production during anaerobic digestion [43]. Mesophilic operation supports process stability and biogas yield by favoring the activity of anaerobic microorganisms [27]. The separation of digestate and biogas streams improves operational control and allows independent handling of solid and gaseous products.

The incorporation of iron(III) oxide in the purification stage allows hydrogen sulfide to be converted into stable sulfide compounds such as Fe₂S₃, reducing corrosive effects and improving biogas quality [44]. This mechanism aligns with adsorption-based purification strategies described for biogas treatment systems [45], ensuring safer handling and greater suitability of the biogas for energy applications.

The use of variable feedstock flows allows the assessment of system behavior under realistic production fluctuations commonly observed in artisanal sugar mills. The application of a ±10 % variation in filter cake availability provides a basis for identifying break-even conditions and the most viable production scenarios, considering seasonal and operational variability in the region [46]. This approach improves the robustness of the simulation by accounting for uncertainties in raw material supply.

The biogas yield obtained in this study is higher than the range of 120–155 m³·t⁻¹, equivalent to 12–15.5%, reported for filter cake digestion in previous studies [35]. This difference can be explained by variations in substrate composition and process assumptions applied during simulation. The chemical characteristics of filter cake, influenced by agroecological conditions, crop type, operational efficiency of the sugar mill, and clarification method, play a determining role in biogas production potential [34]. These factors contribute to the variability observed among reported yields and support the use of simulation approaches to evaluate optimized operating scenarios.

The methane fraction obtained is slightly below the typical range of 40–70% reported for biogas systems [47] but close to the value of 42.47% reported for filter cake digestion [48]. The relatively high CO₂ content is associated with the presence of phytosterols in filter cake, whose anaerobic degradation favors the formation of this gas [49]. The detection of NH₃ indicates protein and amino acid degradation during digestion, while the very low SH₄ concentration is consistent with efficient sulfur removal using metal oxides such as Fe₂O₃, which can remove 70–90% of H₂S [50].

From an energy perspective, lower CO₂ concentrations are associated with improved combustion performance and higher calorific values [51]. Under the simulated conditions, the CH₄ fraction obtained from filter cake was sufficiently high and the CO₂ content remained at moderate levels, allowing the direct use of biogas for heating, heat generation, or combined power applications without requiring an additional CO₂ removal stage.

The higher biogas production observed in Compositions 6, 7, and 8 is influenced by differences in raw material composition. For instance, Composition 8 included a mixture with higher alcohol content, which favors microbial activity during the acidogenic phase and contributes to increased biogas volumes. However, despite the volumetric increase, Figure 5 shows a slight decrease in the molar fraction of methane, while carbon dioxide remained the dominant component, reaching 57.2% in Composition 6. This behavior is associated with the anaerobic degradation of phytosterols and long-chain alcohols, metabolic pathways that generate a higher proportion of CO₂ relative to CH₄ [52].

Although present at low concentrations, hydrogen sulfide affects biogas quality due to its corrosive nature. Its limited presence may be associated with the reduced degradation of sulfur-containing amino acids such as cysteine and methionine in the protein fraction of the substrate. Effective H₂S control in anaerobic processes can be achieved through redox condition adjustment and the use of absorbent agents, improving system stability and reducing operational risks [53].

The calorific value obtained for biogas produced from filter cake supports its technical viability as an alternative energy source for panela production systems. Its energy performance is comparable to that of traditional solid fuels commonly used in sugar mills, such as bagasse and Eucalyptus globulus wood, whose heating values have been widely reported in the literature [31,32]. Based on this energy equivalence, the use of 1 m³ of biogas is equivalent to replacing approximately 0.96 kg of wood or 1.85–1.91 kg of bagasse, depending on biogas composition (Table 6).

When extrapolated to continuous annual operation, the estimated biomass savings reach 631.08 t·year⁻¹, corresponding to the conservation of 3.63 ha·year⁻¹ of Amazon rainforest or approximately 902 trees·year⁻¹. This ecological equivalence was established using reference values for aboveground biomass density and average individual tree biomass reported for Amazonian forests [33], providing an environmental interpretation of the energy substitution achieved.

The magnitude of biomass displacement obtained in this study exceeds that reported for domestic biodigester systems, where average savings of 1,856.78 kg of wood·year⁻¹ per installed unit have been documented, equivalent to 0.011 ha·year⁻¹ of preserved forest [54]. This difference is primarily associated with the operational scale of the system, as household biodigesters are designed to supply individual families, whereas the present evaluation considers higher volumes of filter cake generated in artisanal sugar mills.

The estimated wood savings are also consistent with findings from high-altitude rural communities, where biogas systems with methane contents ranging from 35 to 60% have been shown to replace up to 75% of daily firewood consumption, with reported savings of 0.1–0.4 kg per person per day, equivalent to approximately 1 m³ of biogas·person⁻¹·day⁻¹ [55]. This consistency supports the robustness of the simulated results under different geographic and socio-productive contexts.

In addition to energy substitution, replacing the direct combustion of wood and bagasse with biogas contributes to a reduction in emissions of particulate matter, carbon monoxide, and volatile organic compounds, thereby decreasing exposure to pollutants in artisanal heating environments [56]. Furthermore, the diversion of bagasse traditionally used as fuel enables its use as a low-impact agro-industrial input for applications such as composting, artisanal paper production, and lignocellulosic board manufacturing, reinforcing circular resource use and reducing pressure on forest biomass [57,58].

The higher climate change impact observed for biogas is associated with its lower methane concentration and the presence of unburned CH₄ and biogenic CO₂ in its composition. Methane leakage during production and handling plays a critical role, as leakage rates exceeding 2% can offset the climate benefits of biogas when compared with natural gas [59].

Natural gas shows a lower climate impact due to its high purity, with methane concentrations close to 98%, which favors more efficient combustion and lower emissions per unit of useful energy delivered [60]. In contrast, LPG and propane exhibit intermediate climate impacts, reflecting their fossil origin and combustion characteristics.

The acidification potential of biogas is mainly related to the presence of hydrogen sulfide. During anaerobic digestion, a significant fraction of sulfur contained in the substrate can be converted into H₂S, which, if not adequately removed, dominates the acidification profile of biogas systems [61]. This highlights the importance of effective gas purification prior to utilization.

The higher contribution of LPG and propane to photochemical ozone formation is linked to emissions of nitrogen oxides and volatile organic compounds during combustion and handling. These emissions favor tropospheric ozone formation and explain the higher smog-related impacts observed for these fuels [62]. Overall, the comparison shows that although biogas presents higher impacts in certain categories, its renewable origin and potential for emission mitigation through improved control strategies support its role as a sustainable alternative to fossil gaseous fuels.

The higher OHSP values observed during filter cake reception and digestate cleaning are associated with increased exposure to physical, chemical, and biological hazards, including wet waste handling, hydrogen sulfide and ammonia emissions, noise, and elevated temperatures. These stages represent critical points for occupational risk management in biogas production systems. Incorporating occupational indicators into S-LCA frameworks allows environmental performance to be linked with worker health and safety, identifying biomass handling and waste gas management as priority stages [63].

The identification of ATEX hazardous areas is essential for reducing accident risks related to methane and hydrogen sulfide leaks or accumulations in biogas facilities [64]. In this context, the digestate cleaning and discharge stage presents elevated risk due to the simultaneous presence of toxic vapors and thermal stress. Effective risk reduction requires integrated preventive strategies that combine regulatory compliance, organizational measures, and worker training, supported by appropriate personal protective equipment, adequate ventilation, and clear safety signage [65]. The use of respirators with acid gas filters, waterproof gloves, and insulated protective clothing contributes to minimizing exposure and improving occupational safety during high-risk operations.

5. Conclusions

The energy use of filter cake through anaerobic digestion is a viable option for the artisanal sugarcane industry in the province of Pastaza. Composition 6 performed best in the process, with biogas production of up to 1,736.40 m³·day⁻¹ and methane fractions close to 40.7%. This level enables the partial replacement of wood and bagasse used as fuels and is associated with an approximate annual reduction of 631.08 t·year⁻¹ of solid biomass, equivalent to the conservation of around 3.63 ha·year⁻¹ of Amazon rainforest or about 902 trees per year. The environmental analysis showed that the sustainable performance of the system depends on controlling methane leaks and removing H₂S to reduce contributions to climate change and acidification, while the social assessment identified greater occupational risks in the stages of receiving the filter cake and cleaning the digestate.