Next Article in Journal
Energy-Aware Scheduling Based on Marginal Cost and Task Classification in Heterogeneous Data Centers
Next Article in Special Issue
Enhanced Biogas Production by Ligninolytic Strain Enterobacter hormaechei KA3 for Anaerobic Digestion of Corn Straw
Previous Article in Journal
On-Line Optimization of Energy Consumption in Electromagnetic Mill Installation
Previous Article in Special Issue
Enhancement of Food Waste Management and Its Environmental Consequences
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of Biological Efficiency Assessment Methods and Their Application to Full-Scale Biogas Plants

1
State Institute of Agricultural Engineering and Bioenergy, University of Hohenheim, Garbenstraße 9, 70599 Stuttgart, Germany
2
Institute of Agricultural Engineering, Kiel University, Olshausenstraße 40, 24098 Kiel, Germany
3
DBFZ Deutsches Biomasseforschungszentrum Gemeinnützige GmbH, Biochemical Conversion Department, Torgauer Straße 116, 04347 Leipzig, Germany
4
Rytec GmbH, Consulting and Research, Pariser Ring 37, 76532 Baden-Baden, Germany
5
Tropics and Subtropics Group, Institute of Agricultural Engineering, University of Hohenheim, 70599 Stuttgart, Germany
*
Author to whom correspondence should be addressed.
Energies 2021, 14(9), 2381; https://doi.org/10.3390/en14092381
Submission received: 19 March 2021 / Revised: 14 April 2021 / Accepted: 20 April 2021 / Published: 22 April 2021
(This article belongs to the Special Issue Enhanced Biogas Production)

Abstract

:
For calculation of biological efficiency of a biogas plant (BP), it is required to determine the specific methane potential (SMP) of the substrate. A study comparing available methods for determination of SMP and the comparison with data of full-scale BPs is missing but necessary according to the differences in process conditions between both. Firstly, mass and mass associated energy balances of 33 full-scale BPs were calculated and evaluated. The results show plausible data for only 55% of the investigated BPs. Furthermore, conversion and yield efficiencies were calculated according to six different methods for SMP determination. The results show a correlation between the measured on-site specific methane yield and the calculated SMP by methods based on biological degradability. However, these methods underestimate the SMP. Calculated SMPs based on calorific values are higher, but less sensitive. A combination of biochemical and energetical methods is a promising approach to evaluate the efficiency.

1. Introduction

In Germany alone, approximately 8700 biogas plants (BP) were in operation in the year 2020 with an installed capacity of 6.2 GWel in total [1]. However, the installation of new BPs has come to a near standstill in recent years, mainly due to the lower feed-in tariffs and premiums for electricity generation from the amendments to the Renewable Energy Sources Act (EEG) since 2012 [2]. At the same time, BPs built before 2012 will phase out from their guaranteed 20-year time period of tariffs and premiums in the coming years [3]. In order to compensate for falling tariffs, the substrate conversion rate has to be increased. Therefore, methods are needed to assess the biological process efficiency. The efficiency of a conversion process like anaerobic digestion can be defined by the extent of conversion of an initially available potential of the substrate in two ways: (i) the, produced methane related to the methane potential of the substrate called yield efficiency (YE) in the following and (ii) the actually used potential, i.e., the differences between the methane potential of the substrate and the residual methane potential of the digestate related to the methane potential of the substrate, called conversion efficiency (CE) in the following.
While the substrate is weighed and recorded in most full-scale biogas plants, the output in the form of digestate and biogas is hardly ever measured and must therefore be estimated. Schievano et al. [4] estimated the methane production of the full-scale biogas plant based on the electricity generation of the combined head and power plant (CHP). Thus, under steady state conditions the mass of digestate can be calculated from the mass balance of the input (substrate) and output (biogas and digestate).
In addition to the mass flows, the methane potential of substrates and digestates must be estimated. Various approaches to this problem can be found in the literature. Several studies are based on organic dry matter content (ODM) degradation (CEODM) [5,6,7]. Since the range of substrates typically used in practice varies greatly in terms of the proportion of anaerobically degradable organic matter, a substrate-independent comparison of BPs with this method is not possible [8]. To overcome this disadvantage, other methods based on a more precise determination of the anaerobic conversion potential of substrates have been developed. One of the most common methods is the so-called biochemical methane potential test (BMP), which is a laboratory batch test performed according to VDI 4630 [9]. Anaerobic biodegradability and microbial growth are considered by the method and the methane potential can be calculated directly for any kind of substrates and digestates. Several impact factors were reported and show which efforts are necessary for getting reproducible results [10,11]. To create a common database for specific methane yields of commonly used substrates, BMP results of various laboratories were compiled and published by the Association for Technology and Structures in Agriculture (KTBL) [12], which are referred to as “KTBL values” in the following.
However, the KTBL values do not consider the effects of natural fluctuation in dry matter content (DM), ODM and the degradability of ODM. These parameters are considered by a method based on fermentable organic dry matter (FOM), which was developed by Weißbach [13]. The FOM represents a part of ODM, from which the non-degradable organic fraction is subtracted. Weißbach [14,15,16] defined fermentation coefficients based on feeding experiments with sheeps where FOM of the substrates was derived from a correlation between fiber content in substrate and fiber content in the excrements [13,14,15,16,17].
In contrast to the aforementioned specific considerations of the biological anaerobic process, in industrial energy conversion processes the energy output is simply related to the calorific value of the substrate, representing the total energy content (tE). To compare the efficiency of the biogas process with other industrial conversion processes, two methods were developed to consider the non-fermentable fraction of the substrates. Mächtig et al. [18] suggested subtracting the energy of the lignin fraction from the gross energy content, as lignin is the largest non-biodegradable fraction in the feedstock and called this the anaerobically degradable energy (adE). However, the method does not take into account that besides lignin, the encrusted parts of hemicellulose and cellulose are also inaccessible to microorganisms. Moreover, the authors reported a relative high error of prediction of the lignin content by this method. Alternatively, Fischer et al. [19] proposed to correct the calorific value by the fermentation coefficient from the FOM methodology to calculate the energy content of the fermentable organic dry matter (EFOM). Herein, the non-degradable part of the substrate is assumed to be lignin.
Although biological efficiency is very important for optimizing biogas plants, there are only few scientific publications on this topic and there is no comparative study of the presented evaluation methods for calculating the yield and conversion efficiency of the biological process in large-scale BPs, nor is there a systematic analysis on their validity.
A state-supported monitoring program with 33 participating full-scale biogas plants has now provided the first opportunity for this research. The aim of this study was to compare different methods of evaluating efficiency in order to identify the most suitable method. Especially, the differences in methodology and accuracy of data collection as well as the problems in transferring laboratory data to full-scale agricultural BPs are included in the consideration.

2. Materials and Methods

2.1. Study Design

In this study, 33 full-scale agricultural BPs in Germany were investigated in the years 2016–2018 for the efficiency of their digestion processes. The descriptive process data of the BPs such as temperature in the first digester, hydraulic retention time (HRT) in heated and gas tight system, the organic loading rate (OLR) and the feed ratio of manure, are listed in Table 1, the substrate ratio of feeding is listed in Table S1. The 33 BPs are a selection from 61 BPs from the biogas measurement program III [20], for which mass and energy balances could be calculated consistently according to the procedure described below. All examined BPs utilize the produced biogas in CHP units exclusively. The average electrical power varied from 73 kW to 1796 kW, calculated by dividing the accumulated amount of electricity generated during the year under investigation by 8760 h. Different methods were applied to quantify the initial potential for biogas formation of the substrates as well as of the residual potential of the digestates to every BP. All methods are based on a common mass balance. Mass balances were calculated based on monthly values of the material flows. Efficiency indicators were calculated as yearly values based on the sum of the monthly values from mass balances.

2.2. Data Collection and Laboratory Analyses

2.2.1. Process Data

For calculation of the different efficiency indicators, process data of the examined BP were necessary. These were recorded automatically from sensors by process control systems or manually by plant operators.
The input mass of every solid substrate was weighed daily by scales. The input volume of every liquid substrate was measured daily via flowmeters or calculated based on livestock units. The concentration of methane in the produced biogas was measured mostly continuously but at least monthly. Operating hours and produced electricity from CHP units was recorded monthly, as well as the feed-in amount of electricity and the amount of internal electricity consumption of the BP. For pilot injection CHP engines, the amount of consumed ignition oil was recorded monthly.

2.2.2. Sampling

The input substrate and digestate material streams had to be characterized for their inherent methane and energy potentials by laboratory analysis. Therefore, sampling of these materials was necessary. The sampling method was standardized. Samples of input substrates were taken monthly. Solid substrates were sampled from silos or piles by taking subsamples at 20 cm depth from different positions and mixing them together. Liquid samples, e.g., manure, were taken after complete mixing of storage tanks. Digestate samples were taken at the last gas-tight digester in flow direction at least every three months. A minimum of 10 L of the digestate were discarded before taking samples from sampling pipes. Afterwards, the samples were cooled down to 2 °C ± 5 °C on field-site and stored at −20 °C ± 1 °C in the laboratory except for digestate samples for BMP test. These samples were cooled down to 4 C ± 1 °C on field-site and stored at the same temperature in the laboratory.

2.2.3. Dry Matter Content (DM) and Organic Dry Matter Content (ODM)

DM and ODM were analyzed for every sample at least in duplicate according to German standards DIN EN 15935 [21]. The DM and ODM of all silage and sugar beet samples were corrected for volatile fatty acids according to Weißbach et al. [13,14,15]. The pH values and acid concentrations were measured in duplicate for the DM/ODM correction of silages and sugar beet. More information about the measurement methods can be found in the literature [22,23,24].

2.2.4. Biochemical Methane Potential Test

The biochemical methane potential (BMP) test was done in three different laboratories, two of them used the Hohenheim biogas yield test (HBT). HBT is a continuous mixed system with a 100 mL-syringe as the reactor. The determination of the methane percentage was done by an infrared-spectrometric methane-sensor Advanced Gasmitter, Pronova Analysetechnik (Berlin, Germany). Between both HBTs only the methane quality measurement was different. BP 5–17 and 24–32 were measured in dry gas and BP 1–4 and 20–23 were measured in wet gas and corrected afterwards. The third laboratory used a Bergedorf fermentation test [9] with a volume of 1.5 L and measured the BP 18–19 and 33. Gas was measured by a tipping cell counter (MilligasCounter, Ritter Apparatebau GmbH, Bochum, Germany). Methane percentage was measured by an infrared sensor (Awite Bioenergie GmbH, Langenbach, Germany).
The BMP test of all substrate samples was done once year for all substrates according to VDI 4630 [9]. Temperature for all tests was 37 °C and the tests stopped, for HBTs, after 35 days and for Bergedorf fermentation after reaching less than 0.5% total methane production in one day. All samples were measured at least in triplicates and with an inoculum to substrate ODM ratio of 2:1. Inocula from a 400 L laboratory reactor was used to investigate samples from BP 5–17 and 24–32, inocula from a wastewater treatment plant to investigate samples from BP 1–4 and 20–23 and inocula from 2500 L reactor to investigate samples from BP 18–19 and 33. Hülsemann et al. [11] describes in detail the systems, the inocula and the measurement accuracy. The specific methane potential of the digestate was measured according to VdLuFa [25] at 37 °C for 60 days.

2.2.5. Fiber Content

The raw fiber content of biomass in general describes the organic, non-fat, acid- and alkali-insoluble fraction and was determined for FOM calculation. Raw fiber was determined according to the protocol published by Dittrich-Zechendorf [26]. The fat free sample material was successively treated in boiling sulphuric acid and caustic potash of a defined concentration. The remainder is separated, dried, weighed and finally incinerated at 500 °C. The loss of mass during incineration corresponds to the raw fiber content of the sample. The fiber content of each of the BP’s substrates was analyzed four times a year and in triplicates for statistical validation.

2.2.6. Gross Calorific Value

The gross calorific value (GCV) was determined according to DIN EN ISO 18125 [27]. A correction for sulfur or nitrogen content was omitted. The samples were dried at 105 °C and milled to a particle size smaller than 1 mm prior to analysis. This way, dry matter specific GCVs were measured. The GCV was corrected for volatile fatty acids and alcohols, like done in DM determination for silage samples. For the substrate and digestate of each BP, the gross calorific value was determined four times in a year in triplicates.

2.3. Mass Balance

The system boundary for mass balances was chosen to enclose all gas-tight digesters (Figure 1) in order to describe the biological process alone and also be independent of the BP’s biogas conversion technologies [28]. For mass balances all material streams crossing the system boundary were determined.
All examined BPs performed a continuous digestion process. Feeding intervals were between 20 min and 1 h. Fluctuating liquid filling levels in covered digestate storage tanks, due to agricultural logistics in digestate application, were handled by choosing an investigation period of 12 months. This way, seasonal variations during the observation period were reduced. Therefore, steady state conditions were assumed. Based on the defined system boundary, three material streams had to be considered in the mass balance: substrate, biogas and digestate. Recycled digestate was not considered in the balance, since it is not leaving the system boundary (Figure 1). In case of solid-liquid-separation prior to recycling of the liquid digestate, both fractions were considered in the mass balance.
The following equation was used for steady-state mass balancing:
Δ m = 0 = m sub m gas m dig
where msub is the mass of the substrate, mgas the mass of the biogas and mdig the mass of the digestate.

2.3.1. Substrate

The masses and/or flow rates of substrate input were measured at every BP. Liquid substrate intake was measured or calculated as flow rate and converted to mass assuming a density of 1000 kg m³ due to the high water contents (more than 90%).

2.3.2. Biogas

Measurement equipment for biogas flow is rarely available at the examined agricultural biogas plants. In order to apply a common procedure, the amount of utilized biogas was assumed to be equal to the produced amount of biogas, which was based on the amount of electricity, which was fed into the grid. This value was measured and recorded by grid operators and is the only reliable value that is available for each of the BP. Losses by biogas leakages in pipes, digester roofs or pressure relief valves were not measured and could therefore not be considered.
For the calculation of produced biogas, the following assumptions were made:
  • the transformation loss between the feed-in point and the CHP unit is 2% of the amount of electricity fed in, which is the average value of all examined BPs for which this transformation loss could be calculated based on alignment of electricity measurement at CHP units and at the grid access point;
  • gas leakages and losses, like gas burned in the emergency flares, were not considered;
  • for BPs using more than one CHP unit: The fed-in electricity from the BP was allocated to single CHP units by their rated power and operating hours in the investigated 12-month period;
  • the electrical efficiency of CHP units was assumed to be equal to the published values by the manufacturers minus 3.1%, which reflects the average efficiency loss as determined by Aschmann and Effenberger [29]. Factors for efficiency loss are engine wear, site of installation above sea level, properties at different loads and engine settings. A higher accuracy in this value was not possible with the available data;
  • the biogas was simplified to consist of methane and carbon dioxide only. The methane concentration was measured at every BP at least one time a month. The residual was assumed to be carbon dioxide only. Justification: Water vapor in produced biogas was condensed by gas cooling at every BP. The condensed water was pumped into digestate storage tanks and accounts to digestate mass in the mass balance. Other trace gas components produce a negligible error.
The gross electrical energy produced by the CHP units of a BP was calculated by:
W gross = 1.02 × W grid + W oc
where Wgross gross electricity, Wgrid electricity feed in grid and WOC is the amount of produced electricity which is used internally by the BP. 31 of the examined biogas plants inject the total produced electricity, so that WOC equals zero in these cases.
For BPs with more than one CHP unit, the ratio ak of the electricity produced by a single CHP k was calculated by:
a k = t k × P k k t k × P k
The energy of utilized biogas in all CHP units of a BP was calculated by:
E gas = k a k ( η k 0.031 ) × W gross
with ηk the electrical efficiency of the CHP units.
The utilized volume of biogas was calculated by:
V gas = E gas W PIO NCV CH 4 × x CH 4
where NCV CH 4 the net calorific value, x CH 4 the methane ratio and WPIO the amount of energy from ignition oil, which was subtracted for BPs, which used pilot injection engines as CHP. Finally, the mass of biogas was calculated by:
m gas = V gas × x CH 4 × ρ CH 4 + 1 x CH 4 × ρ CO 2
With the densities ρ CH 4 = 0.72 kg m−3 and ρ CO 2 = 1.98 kg m−3 at temperature of 0 °C and pressure of 101,325 kPa
Relating the utilized volume of biogas to the amount of ODM of fed substrate in the same timespan leads to the on-site specific methane yield (SMY):
On - site   SMY = V gas m sub , ODM

2.3.3. Digestate

The mass and volume of digestate was not recorded at the examined BPs and therefore for the calculation using the mass balance:
m dig = m sub m gas

2.3.4. Organic Dry Matter Content (ODM) Material Balance

An ODM material balance was used to evaluate the plausibility of the mass balance. Mass of ODM for each substrate, gas and digestate (i) was calculated for each month j = 1, …, 12 and summed up for the year.
Two cases were investigated. First, it was assumed, that biogas was produced from ODM only. Second, water incorporation of 10% of the mass of biogas was considered, originating from cellulose degradation during digestion.
Evaluation was done by calculating the ODM balance residual (ODMBR) in the following equation:
ODMBR   % = sub   i j m sub , i , j × DM i , j × ODM i , j j m gas , j × 1 w H 2 O dig   g j m dig , g , j × DM g , j × ODM g , j sub   i j m sub , i , j × DM i , j × ODM i , j × 100 %
where w H 2 O = 0.1 for 10 % water incorporation or w H 2 O = 0 for neglecting water incorporation.

2.4. Energy Balance

The system boundary for energy and energy balances was chosen equal to the one for mass balance (Figure 1). Energy flows associated with mass flows were examined to evaluate energy-based efficiency indicators as well as to prove the plausibility of the used process data for mass balancing. According to the defined system boundary the only mass associated energy intake is the substrate input. Energy outputs are biogas and digestate. The residual of the input energy leaves the process as a loss of heat and entropy as a result of the conversion process of the organic matter.
The following balance was applied to calculate the part of the substrate energy, which is lost in the digestion process as heat, entropy or leakage of biogas:
E loss = E substrate E digestate E gas
Esubstrate is the potential of the substrate fed in a BP as input. Edigestate is the potential of the digestate as output, Egas is the potential of produced biogas and Eloss is the potential lost during the process (e.g., burned gas by emergency flare, leakage, open pressure relief valve, heat loss) (Equation (10)).
The energy balance residual (EBR), was calculated as:
EBR   % = sub   i j m sub , i , j × DM i , j × GCV i , j j V biogas , j × x CH 4 , j × GCV CH 4 digestate   g j m dig , g , j × DM g , j × GCV g , j sub   i j m sub , i , j × DM i , j × GCV i , j × 100 %

2.5. Specific Methane Potential of Substrate Mixtures

Six different methods were used to calculate the annual average ODM-specific methane potential (SMP) of the fed substrate mixtures of the 33 BPs. The formulas represent the ODM-weighted average of the SMP of the single substrates used. For methods based on energy quantities, the SMP was calculated by simply dividing the energy amount by the GCV of methane. By that, conversion losses are neglected.

2.5.1. Biochemical Methane Potential Test (BMP)

The BMP is a batch test, with which the specific methane potential of the substrate can be measured directly. The mixture from direct BMP measurements was calculated by:
SMP BMP , SM = sub   i j m sub , i , j × DM i , j × ODM i , j × SMP BMP , i sub   i j m sub , i , j × DM i , j × ODM i , j

2.5.2. Values According to Literature (KTBL)

The specific methane potential of the substrate mixture using KTBL values from literature was calculated by:
SMP KTBL ,   SM = sub   i j m sub , i , j × DM i , j × ODM i , j × SMP KTBL , i sub   i j m sub , i , j × DM i , j × ODM i , j
SMPKTBL,i is the SMP of a single substrate i from the KTBL- values, which are based on the results of biochemical methane potential test and are given for the most common agricultural substrates in Germany [12]. For rarely used substrates, for which KTBL values were not defined, SMPBMP,i values were used.

2.5.3. Fermentable Organic Dry Matter (FOM)

The FOM is a part of ODM, according to Equation (14). FOM considers that not all of the ODM can be anaerobically degraded. Calculation is done by subtracting a part of the fiber content from ODM, which is a substrate specific calculation of the so called fermentation coefficient fi based on regression equations investigated by fermentation experiments in sheep stomach [13]. For manure, fermentation coefficients were defined based on BMP results [15,16].
FOM i = f i × ODM i
The specific methane potential of the substrate mixture based on FOM was calculated by:
SMP FOM ,   SM = sub   i j m sub , i , j DM i , j × FOM i , j × SMP avFOM sub   i i m sub , i , j × DM i , j × ODM i , j
SMPavFOM is the average specific methane potential of FOM, which is used as a constant factor according to Weißbach [13] (Equation (16)).
SMP avFOM = 420   L   kg 1
For substrates without a regression equation by Weißbach [13,14,15,16,17], the fermentation coefficient f was calculated. The SMPBMP,i was used to determine the fermentation coefficient for FOM calculation by:
f i = SMP BMP , i SMP avFOM

2.5.4. Energy of Fermentable Organic Dry Matter (EFOM)

With EFOM, the FOM is considered energetically by measuring the gross calorific value of ODM and distinguishing degradable and non-degradable fraction using the fermentation coefficient from FOM-methodology according to Fischer et al. [19]. Herein, the non-degradable fraction is assumed to be lignin.
The specific methane potential of the substrate mixture based on the energy of FOM was calculated by:
SMP EFOM , SM = sub   i j m sub , i , j × DM i , j × ODM i , j × GCV FOM , i , j sub   i j m sub , i , j × DM i , j × ODM i , j × GCV CH 4
GCVFOM,i,j is the gross calorific value of FOM, calculated according to equation 7 of Fischer et al. [19]. The gross calorific value of methane GCVCH4 was assumed to be 39.73 MJ m−3 [30].

2.5.5. Anaerobically Degradable Energy (adE)

The adE is the energy of anaerobically degradable parts of ODM, which is calculated by subtracting the energy of lignin from tE. In contradiction to EFOM methodology, the non-degradable part of ODM is not calculated using fermentation coefficients, but by estimating the lignin content of digestates from measured gross calorific values of the digestates. As lignin is assumed to be non-degradable, the energy of lignin in the substrate mixtures has to be the same amount as in the respective digestates.
The specific methane potential of the substrate mixture using the anaerobically degradable energy content was calculated by:
SMP adE , SM = sub   i j m sub , i , j × DM i , j × GCV i , j dig   g j m dig , g , j × DM g , j × ODM g , j × w lig , g ,   j × GCV lig sub   i j m sub , i , j × DM i , j × ODM i , j × GCV CH 4
The lignin content wlig in digestate was calculated from measured gross calorific values of the digestate according to the regression equation M1 by Mächtig et al. [18]. The authors stated a gross calorific value of lignin of 33.74 MJ kg−1 for this model, which was used here for GCVlig.

2.5.6. Total Energy (tE)

The total energy content is calculated from measured gross calorific values and does not consider the different anaerobic degradability of substrates. The specific methane potential of the substrate mixture using the total energy content was calculated by:
SMP tE , SM = sub   i j m sub , i , j × DM i , j × GCV i , j sub   i j m sub , i , j × DM i , j × ODM i , j × GCV CH 4

2.6. Efficiency Indicators

2.6.1. Yield Efficiency

The yield efficiency (YE) is an efficiency indicator, which describes which part of the initial potential of the input was converted to the desired output of methane.
YE = amount   of   produced   methane methane   potential   of   input   substrate
The numerator can either be the mass of methane or the energy of methane. The denominator, which represents the initial potential, has consequently to be given in the same dimension. YE was calculated using different methods to quantify the methane or energy potential.

Yield Based on BMP

The YE based on biochemical methane potential tests YEBMP was calculated as follows:
YE BMP = month   j V gas , j × x CH 4 ,   j sub   i j m sub , i , j × DM i , j × ODM i , j × SMP BMP , i

Yield Based on KTBL

The YE based on KTBL- values YEKTBL was calculated as follows:
YE KTBL = month   j V gas , j × x CH 4 ,   j sub   i j m sub , i , j × DM i , j × ODM i , j × SMP KTBL , i
SMPKTBL,i is the specific methane potential (SMP) from the KTBL- values [12]. Where KTBL values were not defined, SMPBMP,i values were used.

Yield Based on FOM

The YE based on methane potential of fermentable organic dry matter content YEFOM was calculated as follows:
YE FOM = month   j V gas , j × x CH 4 ,   j sub   i j m sub , i , j × DM i , j × FOM i , j × SMP FOM

Yield Based on EFOM

The YE based on the energy of fermentable organic dry matter YEEFOM was calculated as follows:
YE EFOM = month   j V gas , j × x CH 4 ,   j × GCV CH 4 sub   i j m sub , i , j × DM i , j × ODM i , j × GCV FOM , i , j

Yield Based on adE

The YE based on the energy of anaerobically degradable part of the substrate YEadE was calculated as follows:
YE adE = month   j V gas , j × x CH 4 ,   j × GCV CH 4 sub   i j m sub , i , j × DM i , j × GCV i , j dig   g j m dig , g , j × DM g , j × ODM g , j × w lig , g ,   j × GCV lig

Yield Based on tE

The YE based on the total energy of substrate YEtE was calculated as follows:
YE tE = month   j V gas , j × x CH 4 ,   j × GCV CH 4 sub   i j m sub , i , j × DM i , j × GCV i , j

2.6.2. Conversion Efficiency

The conversion efficiency (CE) is an efficiency indicator, which describes which part of the initial potential of the input was converted.
CE = potential   of   input   substrate potential   of   the   digestate potential   of   input   substrate
Different methods were used for the calculation of CE.

Conversion Based on BMP

The CE based on biochemical methane potential tests CEBMP was calculated as follows:
CE BMP = sub   i j m sub , i , j × DM sub , , i , j × ODM sub ,   i , j × SMP BMP , i dig   g j m dig , g , j × DM g , j × ODM g , j × SMP BMP , dig , g , j sub   i j m sub , i , j × DM i , j × ODM i , j × SMP BMP , i
The SMPBMP,dig describe how much gas can be produced out of the digestate by an additional hydraulic retention time of 60 days and 37 °C.

Conversion Based on FOM

The CE based on fermentable organic dry matter CEFOM was calculated as follows:
CE FOM = sub   i j m sub , i , j × DM i , j × ODM i , j dig   g j m dig , g , j × DM g , j × ODM g , j sub   i j m sub , i , j × DM i , j × FOM i , j
Due to missing correlation equations to calculate the residual FOM in the digestate, the amount of degraded organic matter is set as numerator, as by the theory of the method degradation can only occur from FOM.

Conversion Based on EFOM

The CE based on energy fermentable organic dry matter CEEFOM was calculated as follows:
CE EFOM = sub   i j m sub , i , j × DM i , j × GCV i , j dig   g j = 1 12 m dig , g , j × DM g , j × GCV g , j sub   i j m sub , i , j × DM i , j × ODM i , j × GCV FOM , i , j
Due to missing correlation equations to calculate the residual FOM in the digestate and the inherent energy, the amount of converted total energy is set as numerator, as by the theory of this method the converted energy can only be released from degradation of FOM.

Conversion Based on adE

The CE based on the energy of anaerobically degradable part of the substrate CEadE was calculated as follows:
CE adE = sub   i j m sub , i , j × DM i , j × GCV   i , j dig   g j m dig , g , j × DM g , j × GCV g , j sub   i j m sub , i , j × DM i , j × GCV i , j dig   g j m dig , g , j × DM g , j × ODM g , j × w lig , g ,   j × GCV lig
CEadE describes the conversion with the assumption that lignin cannot be degraded.

Conversion Based on tE

The CE based on total energy CEtE was calculated as follows:
CE tE = sub   i j m sub , i , j × DM i , j × GCV i , j dig   g j m dig , g , j × DM g , j × GCV g , j sub   i j m sub , i , j × DM i , j × GCV i , j
GCVg,j describes the energy content in the digestate without distinguishing the non-degradable part.

Conversion Based on ODM

The CE based on the organic dry matter content of the substrate CEODM was calculated as follows:
CE ODM = sub   i j m sub , i , j × DM sub , , i , j × ODM sub ,   i , j dig   g j m dig , g , j × DM dig , g , j × ODM dig , g , j sub   i j m sub , i , j × DM i , j × ODM i , j

3. Results and Discussion

3.1. Mass and Energy Balance

The efficiency determination was based on mass and energy balances. They are fundamental for subsequent calculation of efficiency indicators. The balance residual is determined for ODM material balance and energy balance of the 33 investigated BPs (Equations (8) and (10)). Based on the balance residuals the credibility of measured mass and energy flows can be evaluated [4].
The results show negative as well as positive balance residuals for ODM material and energy balances, which may be caused by several reasons (Figure 2). Positive values indicate lower mass and energy output than input. Negative values indicate a higher mass and energy output than input. Random residuals cannot be completely avoided based on possible non-steady state conditions and the error in sampling regarding the place and time, e.g., the DM and ODM content of substrates show daily variation [31].
The accuracies of DM and ODM measurements between the laboratories have a small deviation of less than 1% and consequently only have a minor effect on the residuals, according to the internal validation of the participating laboratories. Therefore, the calculated residuals cannot solely be based on the measurement accuracy of the laboratory methods. Hence, a more likely reason is the on-site data accuracy.
Positive ODMBRs were observed for 10 BPs. A possible reason for a positive ODMBR is a too high assumed feeding mass by errors in measurement device as well as unknown and unrecorded substrate inputs. Addition of unrecorded substrate streams like rain water or silage seeping water lead to a dilution of digestate and therefore a lower ODM concentration, which consequently results in an underestimation of the digestate’s potential. Another reason for positive ODMBRs could be the assumption of a too high CHP unit efficiency. This results in a lower calculated biogas yield and a positive ODMBR. A leakage of biogas from pipes and pressure relief valves leads to the same error.
In contrast, negative ODMBRs can result from an assumption of a too low CHP unit efficiency. An error in substrate mass measurement devices, e. g. by omitted calibration, can also lead to a negative ODMBR, if the measured mass is smaller than the real input mass.
However, both cases randomly occur across all plants and cannot explain that the most BPs (70%) show negative ODMBRs. A methodical error is the neglect of water incorporation in ODM during the biogas production process. Water incorporation leads to an overestimation of ODM content via biogas output, as biogas is assumed to be built solely from ODM. Based on stoichiometric calculation for cellulose conversion the necessary water addition is 0.11 g g−1 Cellulose. This results in a fraction of 10% of biogas mass build from water. Pröter et al. [22] calculated a value of 11.25% of water incorporation in biogas for a maize silage. Using this factor of 11.25% water incorporation for the substrate mixtures of the BPs leads to an ODMBR of −3% instead of −7% averaged across all investigated biogas plants. However, for BPs with a positive ODMBR, the residuals increase proportionally (maximum +29% instead of +24%). The extent of water incorporation depends on the substrate´s constituents. The investigated BPs use diverse substrates in substrate mixtures, so that using a fixed value for respecting water incorporation for all BPs is not justifiable. Anyway, by disregarding water incorporation, negative ODMBRs are the expected case for accurate process data. Ultimately, the water incorporation only effects the ODM balance for checking the mass balances plausibility and not the mass balance itself.
Even though EBRs are also affected by measurement errors in mass balance, there are other possible factors responsible for the distribution. A positive EBR is observed at more than 80 % of the BPs, which appears to be caused by two main factors. Energy outputs without a relation to the mass flow like heat and entropy loss based on exothermal reactions cannot be detected by the chosen approach. Furthermore, the calculated biogas yield might be too low due to overestimation of CHP unit´s electrical efficiency or due to gas losses. Biogas loss can happen, amongst other reasons, through leakages and the pressure relief valves, just as Liebetrau et al. [32] determined emission factors of agricultural BPs to be in the range of 1.1 to 13.7%CH4 of the total production. A third reason might be, like already mentioned for positive ODMBR, that unrecorded substrate mass was fed to the biogas plant. Even if this unrecorded mass was only water without energy content, it would lead to a wrong mass balance assuming to low digestate mass and therefore would cause a lower energy output by the digestate mass flow.
Negative EBRs are smaller than 5% among all BPs and are only shown by 20% of the BPs. The negative EBRs are linked with the energy balance in which strongly negative ODMBRs also result in negative EBRs. Therefore, reasons for negative EBRs are the same as for negative ODMBRs, i.e., possibly underestimated input mass or underestimated CHP unit efficiency leading to overestimated biogas yields. A negative EBR can be reported only for BPs with a negative ODMBR of more than 14%. This reveals the large impact of the mass balance, as errors here are reproduced in the energy balance.
The reflection of either ODMBR or EBRs alone does not allow a distinct decision regarding most probable errors and the credibility of the material and energy flows. The coherence of mass and energy balance reveals more useful information. From Figure 2, three cases can be distinguished:
(1)
Positive ODM balance residual and positive energy balance residual
Ten BPs show a positive value of ODMBR and EBR. This positive ODMBR increases the positive EBR. For these few plants a clear positive error in feeding mass measurement device, an overestimation of CHP unit efficiency or biogas leakage is assumed. However, correction of this data is not possible in hindsight. Operators of these BPs should check the accuracy of their measurement devices and the efficiency of the CHP unit.
(2)
Negative ODM balance residual and positive energy balance residual
Eighteen BPs show a negative ODMBR and positive EBR. This is attributed to water incorporation and entropy loss. A high data quality for the investigated BPs can be assumed, but more detailed research into water incorporation and entropy loss are necessary to be more specific. From the presented methodology this is the expected case.
(3)
Negative ODM balance residual and negative energy balance residual
Five BPs show negative EBRs, which is generally impossible. However, the negative residuals are comparably small. In coherence with the negative ODMBRs the most probable errors are negative deviations in feeding mass measurement or an underestimated CHP unit efficiency. In this case, the accuracy of mass measurement should also be checked by plant operators.
The residuals follow reasonable directions only in case 2. Therefore, the evaluation suggest that only the 55% of the examined biogas plants give credible results on mass and energy balances by the chosen approach.
BP 12 is a research biogas plant, where data quality is known to be accurate. This BP is one of the BPs with a negative ODM and a positive EBR. This fact strengthens the assumption that BPs with this combination are credible.
The findings can be further distinguished with regard to the diversity of the process parameters of the examined BPs. For BPs fed with 100% energy crops, in every case a negative ODMBR and a positive EBR has been found. The ODM balance of these BP are credible. Different to that all BP with a power production less than 100 kW and also 4 of 6 BP with a manure content higher than 75% show a positive balance residual for energy and ODM. The reason for the latter is, that the installation of measurement devices, e.g., for manure mass flows, is not economically feasible on small-scale BPs. The balances of these plants are less credible.
The ODM and energy balance reveals the problems in carrying out an efficiency analysis. The analysis is based on both balances. The necessary assumptions and the use of on-site process data can result in low accuracy, mainly due to missing measurement equipment or inadequate calibration. The inability to determine water incorporation, entropy loss and gas loss for full-scale biogas plants prevent the determination of a correction factor.

3.2. Specific Methane Potential

The specific methane potentials (SMP) from substrate mixtures of 33 biogas plants (BP) in Figure 3 were determined by literature values of Association for Technology and Structures in Agriculture (KTBL), biochemical methane potential test (BMP), fermentable organic matter (FOM), energy of fermentable organic matter (EFOM), anaerobically degradable energy (adE), total energy (tE) and the on-site measured methane yield (on-site SMY). SMPKTBL,SM, SMPBMP,SM, SMPFOM,SM and SMPEFOM,SM show similar tendency for the specific methane potential between the BPs (Figure 3). The tendency is also similar to the on-site SMY, which is the specific methane yield of the fed substrate mixtures. This indicates a good sensitivity of the methods. However, SMPKTBL,SM, SMPBMP,SM, SMPFOM,SM and SMPEFOM,SM are mostly lower compared to the on-site SMY. The main reason is an underestimation of the potential by the methods, as discussed later. Another possible reason, like false process data, is only a possibility for 5 of 33 BPs (Case 3, Figure 3), by underestimation of the fed substrate mass and/or underestimation of CHP unit efficiency.
SMPAdE,SM and SMPtE,SM are much higher than SMP of the other methods (Figure 3). As the methods for adE and tE give energy values as results, they had to be converted to methane yields for comparison. This was done by dividing the energy values by the gross calorific value of methane (Equations (19) and (20)). By that, losses by chemical conversion of substrates to methane are neglected, leading to higher SMP values. Within the other methods, these losses are methodologically included.
SMPadE,SM are lower than SMPtE,SM, as anaerobically non-degradable lignin in the substrates is subtracted from the total energy potential by using adE. The sensitivity of the calculated SMPtE,SM and SMPadE,SM is lower than for the other methods. This is expected for SMPtE,SM, as this value does not reflect the anaerobic degradability of the substrate mixture. In contrast, SMPadE,SM respects the composition of the substrate mixture but does not respect the accessibility of the degradable components for microbial degradation. The degradable components hemicellulose and cellulose are typically incrusted by lignin, which makes them not accessible for microbial degradation without physical pretreatment.
Figure 2 and Figure 3 indicate an effect between the SMY of a BP and the ODMBR. The SMY correlates also according to the fed manure (Table 1), because of less SMY of manure compared to energy crops.
Figure 3 also reveals the tendency that a negative ODM and energy balance follow in lower estimated SMP compared to the on-site SMY, during positive ODM and energy balance follow in a higher estimated SMP compared to the on-site SMY, according to an aforementioned probable negative/positive feeding mass error. These findings support the assumption, that the data quality is low for 45% of the BPs.
Besides the SMP of the substrate, the methane potentials of the digestate were used to determine the conversion efficiency (CE). Only the methods BMP, adE and tE determine the potential of digestate (Figure 4). For FOM and EFOM a required fermentation coefficient for the digestate is not determinable.
In contrast to Figure 3 the SMP of digestates in Figure 4 were determined by biochemical methane potential test (BMP), anaerobically degradable energy (adE) and total energy (tE), based on the ODM in digestates. Again, the energy of adE and tE was converted to SMP by dividing with the gross calorific value of methane and by that neglecting losses of conversion. The potentials show major differences based on different assumptions in the three methods.
SMPtE,dig was the highest calculated SMP and the values are even higher, than the SMP of the substrates calculated by this method, which is caused by using ODM-specific values. In accordance with the lignin ratio in ODM, the ODM specific energy in digestate is higher than in substrate. Since low calorific components of ODM are degraded (mainly carbohydrates with around 17–18 MJ kg−1), high calorific lignin (around 25–26 MJ kg−1 [23]) remains undigested. This is also the reason why the pattern of tE and adE of the 33 BPs is symmetrical (Figure 4). Higher SMPtE,dig means higher lignin contents, resulting in lower SMPadE,dig.
The SMPadE,dig in Figure 4 is around 3 times higher than SMPBMP,dig. Most probable reason for the difference is that hemicellulose and cellulose in the digestates are incrusted by lignin and not accessible for microbial degradation in BMP test. In median, based on the results of this study, 39% of the energy potential of the digestates are lignin. Additional 48% of the energy potential of the digestates are not accessible to microbial degradation, as calculated by the difference of adE and BMP. This should be the maximum possible potential, which is achievable with further treatment. Further research on this topic could lead to a useful parameter set in future. However it has to be mention, that the temperature effects the SMPBMP,dig strongly [33].

3.3. Conversion and Yield Efficiency

The yield efficiencies (YE) were calculated for all BPs according to six different methods based on Equation (21). The results are listed in Table 2 grouped in the three described cases. The calculated YEs show a similar distribution like the SMPs calculated for different methods in Figure 3. Table 2 clearly shows that the mass balance also effects the results of YE for every method. BPs with a positive ODM and energy balance residual (case 1) show low YE for every method according to a positive feeding mass error.
Different from that, YE of BP with negative ODMBRs and EBRs (case 3) show a tendency for high YE according to a negative feeding mass error, but the difference to case 2 is not that obvious. Similar to YE in Table 2, CEs of all BPs are listed in Table 3. The CE was determined based on the potentials of the substrate and the digestate determined by BMP, FOM, EFOM, adE, tE and ODM (Table 3). Dependencies between process data, e.g., process temperature, HRT or OLR, and efficiency indicators were not observed and are therefore not presented.
The KTBL values and the BMP values are based on measured values in a BMP test. Both show a YE higher than 100%. It appears that the methods underestimate the real SMP. The values of the research BP (BP 12) support these results even if the values are near to 100% (BMP 104% and KTBL 101%). Efficiency values higher than 100% was also reported for both in literature [34,35]. Underestimation of SMPBMP,SM and SMPKTBL,SM can result, due to the effects of co-digestion in a full-scale BP instead of mono-digestion in the laboratory batch test [36,37]. Pöschel et al. [38] estimated 10% higher values. Other possible reasons could be methodological differences between lab and full-scale, like the higher temperature in a full-scale biogas plant, a higher consumption for bacterial growth based on a missing steady state condition in laboratory tests as well as unacclimated inoculum or a too short HRT in laboratory tests [39].
CEBMP is the state-of-the-art method to determine the CE. A CE higher than 100 % is not possible for this method. The CE reveals the achievable extra potential by extending the HRT of gas-tight digesters for 60 days in a 37 °C heated digester. In Germany, the total HRT of gas-tight digesters built after 2016 and the HRT of a secondary digester built after 2011 are already higher than 150 days according to EEG 2017. This leads to a high CEBMP found in this study and also to a small spread in the values, which corresponds to results by Ruile et al. of 21 German biogas plants [5]. The CEBMP reveals a small potential in increasing the digestion duration nowadays compared to lower CEBMP in other former studies [5,33,40], which is supported by a comparison of the SMPBMP to the SMPBMP of former studies [34,41]. Further optimization cannot be done for most BPs based on this. However, it is the only CE method, which can show the reachable potential by a longer HRT. This information is of interest for BPs operation in different countries today. For three outliers with a CEBMP lower than 90%, clear reasons can be found, e.g., low retention time in heated system by feeding substrate with a high lignin content.
CEFOM and CEEFOM show a higher spread for CE, but CE and YE higher than 100% can be reported for both. Values higher than 100% for YEEFOM and YEFOM are reasoned by underestimated SMP of the substrate. Lower residual ODM masses in the digestates than are expected by the fermentation coefficients used in these methods are responsible for CEEFOM and CEFOM values above 100%. Based on the results, the applied substrate specific fermentation coefficients are too low. In order to utilize FOM and EFOM for CE or YE calculations, more research on fermentation coefficients is necessary. Especially, the variability of manure cannot be considered adequately by using a fixed fermentation coefficient, which only depends on the type of manure. For the research BP, both methods show a value lower than 100% (FOM 95% and EFOM 86%). Weißbach [42] analyzed the efficiency of three different biogas fermenters of the same biogas plant based on FOM and found CEs between 96.4 and 99.5% and YEs between 99.7 and 103.7%.
YEadE and YEtE are throughout lower than 100%, which is lower than found with other methods. For the research plant BP 12 adE reach a value of 71% and tE of 66%. As YEadE is applied for the first time here, comparative literature values are missing. However, it is possible to calculate expected values by stoichiometry with the so called Buswell’s formula [43]. For example, for cellulose a maximum YEadE of 95% can be expected based on a GCV of cellulose of 17.3 MJ kg−1. Based on the reference value for the specific methane yield of cellulose from VDI 4630 [9], a YEadE of 86% is expected. For cellulose, YEadE and YEtE are equal, as cellulose is completely degradable under anaerobic conditions. The example shows, that the calculated YEadE of the examined BPs are in a reasonable range.
CEadE and CEtE are lower than CE with other methods and show a higher spread compared to CEBMP. The CEadE is lower than CEBMP in every case as expected by the different residual SMP in the digestates. Comparison of CEadE and CEtE shows the energy potential in lignin. The study reveals that in median 9% of the total energy of the substrate is stored in lignin. CEtE and YEtE can be used as a reference to compare the biogas process with other biomass conversion processes, like ethanol production, pyrolysis or biomass combustion. Moreover, CEtE is revealing how much additional energy can be utilized from the digestate. 11–7% of the total energy potential is not used by the BPs according to the results of this work. Bio-economy concepts could increase the utilization of the potential. However, CEtE and YEtE are also strongly affected by the used substrate similar to CEODM and thereby hinder the comparison of BPs, which use substrates with different degradability.
CEODM is used several times for determining the CE and a clear divergence between the BPs can be reported. It is clearly impacted by the fed substrate. The ODM of maize is nearly completely anaerobically degradable, because of the low lignin content. Compared to that, cattle manure has a higher lignin content and a lower CEODM. Therefore a comparison of two BPs with different feedstock like done here would not be useful []. CEtE and CEODM are strongly dependent to each other. Only for one BP the two parameters show 11% difference to each other. In every other case CEODM is 1–4% higher than CEtE as expected, due to the linkage of ODM content and gross calorific value.
As shown, CE or YE are strongly dependent on the SMP of the substrate. KTBL, BMP, FOM and EFOM underestimate the substrate potential. AdE and tE are more coherent, but the sensitivity to the on-site methane yields of the BPs is low and the substrate´s degradability, as an interaction between biomass composition and accessibility for degradation, is not fully considered. Hence, more research is necessary to investigate and implement influencing factors on SMP, especially to describe the differences between laboratory scale, where most of the determination methods were developed and commercial full-scale digestion, where the methods are to be applied. Typically, the accuracy of measurements is higher in laboratory experiments, but process conditions differ to those in full-scale biogas plants. From this point of view, it was a necessity to compare the methods with data from full-scale biogas plants, even though high measurement accuracy is difficult to achieve, as shown with presented mass and energy balances.

4. Conclusions

Overall, the results of the current research show that benchmarking agricultural biogas plants is challenging regardless of the method chosen. It is worthwhile to check the reliability of the recorded operating parameters by evaluating mass and energy balances, as shown by the fact that plausible results were obtained for only 55% of the investigated biogas plants. Regular calibration of the measuring instruments is therefore highly recommended for determining the biological efficiency but is repeatedly neglected under practical conditions. However, this had no influence on the comparison of the different methods, as the statements were similar for all BPs.
A conclusive recommendation for a specific method to evaluate efficiency is not possible. YEKTBL, YEBMP, YEFOM and YEEFOM show a good relation to the field measurements, but all underestimate the methane potentials of the substrate. The probable reason for this is an underestimation of the fermentable fraction of the substrate. These conclusions are supported by the fact that CEFOM and CEEFOM being higher than 100%. To avoid this, further research is necessary, e.g., towards substrate characterization. The results of YEadE and YEtE, as well as CEadE and CEtE, are more coherent, but the sensitivity is low and the degradability of the substrate as an interaction between biomass composition and degradability is not fully considered. Based on the YE data, an evaluation of the BPs is not possible and not done here. A combination of CEBMP and CEadE is a promising approach, as it can reveal the achievable potential of further substrate pretreatments. This combination should be investigated in further research. Besides the values of CEBMP, an evaluation of biogas plant operation is difficult as reference values are missing.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/en14092381/s1, Table S1: Substrates ratio of feeding of the 33 investigated biogas plants (BP), MS = maize silage, GS = gras silage, WPS = whole plant silage, SB = sugar beet, CE = cereals, PO = Potato, LCM = liquid cow manure, SCM = solid cow manure, HM = horse manure, CM = chicken manure, LPM = liquid pig manure, SM = sheep manure, ML = molasses of lactose production, BT = brewer’s spent grains, FO = Fodder residuals, OT = others.

Author Contributions

Conceptualization, B.H., T.M. and M.P.; methodology, B.H., T.M. and M.P.; software, B.H. and T.M.; validation, B.H., T.M. and M.P.; formal analysis, B.H., T.M. and M.P.; investigation, B.H., T.M. and M.P.; resources, B.H., T.M., M.P., H.O., J.L., E.H. and J.M.; data curation, B.H., T.M. and M.P.; writing—original draft preparation, B.H., T.M. and M.P.; writing—review and editing, H.O., J.L., E.H. and J.M.; visualization, B.H., T.M. and M.P.; supervision, H.O., J.L., E.H. and J.M.; project administration, H.O., J.L., E.H. and J.M.; funding acquisition, H.O., J.L., E.H. and J.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Fachagentur Nachwachsende Rohstoffe (FNR), Gülzow, Germany; and the Federal Ministry of Food and Agriculture (BMEL), Berlin, Germany Germany [grant number 22403715]; This work was done in the project biogas measurement program III.

Institutional Review Board Statement

Not Applicable.

Informed Consent Statement

Not Applicable.

Data Availability Statement

Not Applicable.

Acknowledgments

Thanks goes to the project partners Bavarian State Research institute for Agriculture (LfL) and to person Robert Kliche and Matthias Effenberger. Furthermore: we would like to thank Jan Postel, Erik Fischer and Markus Lenhart for their contributions.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. FNR. Entwickelung der Biogasanlagenstandorte in Deutschland. Available online: https://mediathek.fnr.de/grafiken/daten-und-fakten/bioenergie/biogas/entwicklung-der-biogasanlagenstandorte-in-deutschland.html (accessed on 14 April 2021).
  2. Torrijos, M. State of Development of Biogas Production in Europe. Procedia Environ. Sci. 2016, 35, 881–889. [Google Scholar] [CrossRef]
  3. Daniel-Gromke, J.; Rensberg, N.; Denysenko, V.; Stinner, W.; Schmalfuß, T.; Scheftelowitz, M.; Nelles, M.; Liebetrau, J. Current Developments in Production and Utilization of Biogas and Biomethane in Germany. Chem. Ing. Tech. 2008, 90, 17–35. [Google Scholar] [CrossRef]
  4. Schievano, A.; D’Imporzano, G.; Salati, S.; Adani, F. On-field study of anaerobic digestion full-scale plants (Part I): An on-field methodology to determine mass, carbon and nutrients balance. Bioresour. Technol. 2011, 102, 7737–7744. [Google Scholar] [CrossRef] [PubMed]
  5. Ruile, S.; Schmitz, S.; Mönch-Tegeder, M.; Oechsner, H. Degradation efficiency of agricultural biogas plants—A full-scale study. Bioresour. Technol. 2015, 178, 341–349. [Google Scholar] [CrossRef] [PubMed]
  6. Hartmann, H.; Ahring, B.K. Anaerobic digestion of the organic fraction of municipal solid waste: Influence of co-digestion with manure. Water Res. 2005, 39, 1543–1552. [Google Scholar] [CrossRef] [Green Version]
  7. Demirer, G.N.; Chen, S. Effect of retention time and organic loading rate on anaerobic acidification and biogasification of dairy manure. J. Chem. Technol. Biotechnol. 2014, 79, 1381–1387. [Google Scholar] [CrossRef]
  8. Schievano, A.; D’Imporzano, G.; Orzi, V.; Adani, F. On-field study of anaerobic digestion full-scale plants (Part II): New approaches in monitoring and evaluating process efficiency. Bioresour. Technol. 2011, 102, 8814–8819. [Google Scholar] [CrossRef]
  9. VDI 4630. Fermentation of organic materials-Characterisation of the Substrate, Sampling, Collection of Material Data, Fermentation Tests; Verlag des Vereins Deutscher Ingenieure: Düsseldorf, Germany, 2016. [Google Scholar]
  10. Weinrich, S.; Schäfer, F.; Bochmann, G.; Liebetrau, J. Value of batch tests for biogas potential analysis. Method comparison and challenges of substrate and efficiency evaluation of biogas plants. IEA Bioenergy Task 2018, 37, 13–28. [Google Scholar]
  11. Hülsemann, B.; Zhou, L.; Merkle, W.; Hassa, J.; Müller, J.; Oechsner, H. Biomethane Potential Test: Influence of Inoculum and the Digestion System. Appl. Sci. 2020, 10, 2589. [Google Scholar] [CrossRef] [Green Version]
  12. KTBL. Faustzahlen Biogas, 3rd ed.; KTBL: Darmstadt, Germany, 2013. [Google Scholar]
  13. Weißbach, F. On Assessing the Gas Production Potential of Renewable Primary Products. Landtechnik 2008, 63, 356–358. [Google Scholar] [CrossRef]
  14. Weißbach, F. Gas production potential of fresh and ensiled sugar beets in biogas production. Landtechnik 2009, 64, 394–397. [Google Scholar] [CrossRef]
  15. Weißbach, F. The gas forming potential of pig slurry in biogas production. Landtechnik 2011, 66, 460–464. [Google Scholar] [CrossRef]
  16. Weißbach, F. The gas forming potential of dry chicken dung in biogas production. Landtechnik 2012, 67, 299–304. [Google Scholar] [CrossRef]
  17. Weißbach, F. Gas production potential of forage and cereal crops in biogas production. Landtechnik 2009, 64, 317–321. [Google Scholar] [CrossRef]
  18. Mächtig, T.; Moschner, C.R.; Hartung, E. Monitoring the efficiency of biogas plants—Correlation between gross calorific value and anaerobically non-degradable organic matter of digestates. Biomass Bioenergy 2019, 130, 105389. [Google Scholar] [CrossRef]
  19. Fischer, E.; Postel, J.; Ehrendreich, F.; Nelles, M. Using the mean fuel efficiency to energetically assess agricultural biogas plants. Agric. Eng. 2016, 71, 139–154. [Google Scholar] [CrossRef]
  20. Biogas-Messprogramm III, 1st ed.; Fachagentur Nachwachsende Rohstoffe: Gülzow, Germany, 2021.
  21. DIN EN 15935:2012-11, Sludge, Treated Biowaste, Soil and Waste—Determination of Loss on Ignition; Beuth Verlag GmbH: Berlin, Germany, 2012.
  22. Pröter, J.; Weinrich, S.; Hofman, J.; Kube, J. Mass balancing of biogas plants. In Collection of Methods for Biogas, 2nd ed.; Liebetrau, J., Pfeiffer, D., Eds.; DBFZ: Leipzig, Germany, 2020; pp. 381–392. [Google Scholar]
  23. Apelt, M. Determination of aliphatic, organic acids and benzaldehyde with headspace GC. In Collection of Methods for Biogas, 2nd ed.; Liebetrau, J., Pfeiffer, D., Eds.; DBFZ: Leipzig, Germany, 2020; pp. 64–68. [Google Scholar]
  24. Steinbrenner, J.; Nägele, H.-J.; Buschmann, A.; Hülsemann, B.; Oechsner, H. Testing different ensiling parameters to increase butyric acid concentration for maize silage, followed by silage separation and methane yield potential of separated solids residues. Bioresour. Technol. Rep. 2019, 7, 100193. [Google Scholar] [CrossRef]
  25. VDLUFA. Federation of German Agricultural Investigation and Research Institutes; VDLUFA Verlag: Darmstadt, Germany, 2007. [Google Scholar]
  26. Dittrich-Zechendorf, M. Determination of total Kjeldahl nitrogen and crude protein. In Collection of Methods for Biogas, 2nd ed.; Liebetrau, J., Pfeiffer, D., Eds.; DBFZ: Leipzig, Germany, 2020; pp. 90–100. [Google Scholar]
  27. DIN EN ISO 18125:2017-08, Solid Biofuels—Determination of Calorific Value; Beuth Verlag GmbH: Berlin, Germany, 2017.
  28. Havukainen, J.; Uusitalo, V.; Niskanen, A.; Kapustina, V.; Horttanainen, M. Evaluation of methods for estimating energy performance of biogas production. Renew. Energy 2014, 66, 232–240. [Google Scholar] [CrossRef]
  29. Aschmann, V.; Effenberger, M. The electrical efficiency of biogas-driven combined heat and power plants. Agric. Eng. 2012, 67, 118–121. [Google Scholar] [CrossRef]
  30. DIN EN ISO 6976:2016-12, Natural Gas—Calculation of Calorific Values, Density, Relative Density and Wobbe Indices from Composition; Beuth Verlag GmbH: Berlin, Germany, 2016.
  31. Banks, C.J.; Chesshire, M.; Heaven, S.; Arnold, R. Anaerobic digestion source-segregated domestic food waste: Performance assessment by mass and energy balance. Bioresour. Technol. 2011, 102, 612–620. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Liebetrau, J.; Reinelt, T.; Clemens, J.; Hafermann, C.; Friehe, J.; Weiland, P. Analysis of greenhouse gas emissions from 10 biogas plants within the agricultural sector. Water Sci. Technol. 2013, 67, 1370–1379. [Google Scholar] [CrossRef] [PubMed]
  33. Angelidaki, I.; Boe, K.; Ellegaard, L. Effect of operating conditions and reactor configuration on efficiency of full-scale biogas plants. Water Sci. Technol. 2005, 52, 189–194. [Google Scholar] [CrossRef]
  34. Lansing, S.; Hülsemann, B.; Choudhury, A.; Schueler, J.; Lisboa, M.S.; Oechsner, H. Food waste co-digestion in Germany and the United States: From lab to full-scale systems. Resour. Conserv. Recycl. 2019, 148, 104–113. [Google Scholar] [CrossRef]
  35. KTBL. Gasausbeute in Landwirtschaftlichen Biogasanlagen. KTBL-Heft 107, 3rd ed.; KTBL: Darmstadt, Germany, 2015; p. 27. [Google Scholar]
  36. Mata-Alvarez, J.; Macé, S.; Llabrés, P. Anaerobic digestion of organic solid wastes. An overview of research achievements and perspectives. Bioresour. Technol. 2000, 74, 3–16. [Google Scholar] [CrossRef]
  37. Astals, S.; Batstone, D.J.; Mata-Alvarez, J.; Jensen, P.D. Identification of synergistic impacts during anaerobic co-digestion of organic wastes. Bioresour. Technol. 2014, 169, 421–427. [Google Scholar] [CrossRef] [Green Version]
  38. Pöschel, M.; Ward, S.; Owende, P. Evaluation of energy efficiency of various biogas production and utilization pathways. Appl. Energy 2010, 87, 3305–3321. [Google Scholar] [CrossRef]
  39. Browne, J.D.; Murphy, J.D. Assessment of the resource associated with biomethane from food waste. Appl. Energy 2013, 104, 170–177. [Google Scholar] [CrossRef]
  40. Weiland, P.; Gemmeke, B.; Rieger, C. Biogas-Messprogramm II: 61 Biogasanlagen im Vergleich; Gülzow-Prüzen, Germany, 2009; ISBN 978-3-9803927-8-5. Available online: https://www.fnr-server.de/ftp/pdf/literatur/pdf_385-messprogramm_ii.html (accessed on 20 April 2021).
  41. Zhou, L.; Hülsemann, B.; Cui, Z.; Merkle, W.; Sponagel, C.; Zhou, Y.; Guo, J.; Dong, R.; Müller, J.; Oechsner, H. Operating Perfomance of Full-Scale Agricultural Biogas Plants in Germany and China: Results of a Year-Round Monitoring Program. Appl. Sci. 2021, 11, 1271. [Google Scholar] [CrossRef]
  42. Weißbach, F. Degree of utilization of primary renewable products in biogas production. Landtechnik 2009, 64, 18–21. [Google Scholar] [CrossRef]
  43. Symons, G.E.; Buswell, A.M. The Methane Fermentation of Carbohydrates. J. Am. Chem. Soc. 1933, 55, 2028–2036. [Google Scholar] [CrossRef]
Figure 1. System boundary for efficiency assessment; starting with the first digester and ending with the last covered digester (CHP—combined heat and power plant).
Figure 1. System boundary for efficiency assessment; starting with the first digester and ending with the last covered digester (CHP—combined heat and power plant).
Energies 14 02381 g001
Figure 2. Balance residual of the organic dry matter content (ODM) material balance without and with considering a 10% water incorporation and the energy balance for 33 biogas plants (BP). A positive value means missing mass/energy in the output. BPs are grouped in three categories (Case 1: positive ODM and energy balance, Case 2: negative ODM balance and positive energy balance and Case 3: negative ODM and energy balance), sorted by the three cases at first and the ODM balance residual without consideration of water incorporation as second.
Figure 2. Balance residual of the organic dry matter content (ODM) material balance without and with considering a 10% water incorporation and the energy balance for 33 biogas plants (BP). A positive value means missing mass/energy in the output. BPs are grouped in three categories (Case 1: positive ODM and energy balance, Case 2: negative ODM balance and positive energy balance and Case 3: negative ODM and energy balance), sorted by the three cases at first and the ODM balance residual without consideration of water incorporation as second.
Energies 14 02381 g002
Figure 3. The specific methane potential (SMP) from substrate mixtures of 33 biogas plants (BP) determined by literature values of Association for Technology and Structures in Agriculture (KTBL), biochemical methane potential test (BMP), fermentable organic matter (FOM), energy of fermentable organic matter (EFOM), anaerobically degradable energy (adE), total energy (tE) and the on-site measured methane yield (on-site SMY). BPs are separated according to three cases (Case 1: positive ODM and energy balance, Case 2: negative ODM balance and positive energy balance and Case 3: negative ODM and energy balance).
Figure 3. The specific methane potential (SMP) from substrate mixtures of 33 biogas plants (BP) determined by literature values of Association for Technology and Structures in Agriculture (KTBL), biochemical methane potential test (BMP), fermentable organic matter (FOM), energy of fermentable organic matter (EFOM), anaerobically degradable energy (adE), total energy (tE) and the on-site measured methane yield (on-site SMY). BPs are separated according to three cases (Case 1: positive ODM and energy balance, Case 2: negative ODM balance and positive energy balance and Case 3: negative ODM and energy balance).
Energies 14 02381 g003
Figure 4. The spec. methane potential of digestate (SMP) of 33 biogas plants (BP) determined by biochemical methane potential test (BMP), anaerobically degradable energy (adE) and total energy (tE), based on the ODM in digestates.
Figure 4. The spec. methane potential of digestate (SMP) of 33 biogas plants (BP) determined by biochemical methane potential test (BMP), anaerobically degradable energy (adE) and total energy (tE), based on the ODM in digestates.
Energies 14 02381 g004
Table 1. Average electrical power production, temperature of the first digester, number of hydrolysis stages (HY), primary digesters (FD), secondary digesters (SD) and storage tanks (ST), hydraulic retention time (HRT), organic loading rate (OLR), manure share of the 33 investigated biogas plants (BP). BP 12 is a research biogas plant with high measurement accuracy.
Table 1. Average electrical power production, temperature of the first digester, number of hydrolysis stages (HY), primary digesters (FD), secondary digesters (SD) and storage tanks (ST), hydraulic retention time (HRT), organic loading rate (OLR), manure share of the 33 investigated biogas plants (BP). BP 12 is a research biogas plant with high measurement accuracy.
BPAverage Electrical PowerTemperature First DigesterStagesHRT HeatedHRT Gas-TightOLRManure Share
kW°CHY/FD/SD/STddkgODM d−1 m−3%mass
173420/1/0/1501533.091.4
253227–331/1/1/473732.275.5
374440/1/1/11481481.181.9
467136–400/1/2/12312751.036.4
51229450/1/1/2722133.358.0
677420/1/0/1732211.982.1
7498440/1/1/21273463.432.9
8209430/1/0/266664.910.5
931638–470/1/1/21041563.137.7
10358420/1/1/12252251.67.7
1150847–530/1/1/1591164.834.9
12207450/2/1/01421421.843.6
13512400/2/1/261612.667.9
14451440/1/0/1421684.950.9
15942410/1/1/181813.80.3
16469400/1/1/11131132.332.5
17170643–450/2/1/1721184.30.0
18649450/1/1/11331892.50.0
19571431/2/0/1631293.556.3
2019939–450/1/0/1651342.552.9
211796430/2/1/2731564.10.0
2263543–490/1/1/1711683.251.5
2345952–590/1/1/51011012.934.6
2438135–430/1/0/178782.251.6
25560430/1/1/2892181.955.9
26712420/1/0/1541043.175.6
27739430/1/1/11241922.50.0
28557440/2/0/645453.673.0
29371430/1/1/1872082,342.7
30515420/1/0/1812263.632.0
31511421/2/1/11201201.462.1
32512440/1/1/161961.684.9
33975500/1/2/2592724.351.7
Table 2. Yield efficiency (YE) of 33 biogas plants (BP) determined by literature values of Association for Technology and Structures in Agriculture (KTBL), biochemical methane potential test (BMP), fermentable organic matter (FOM), energy fermentable organic matter (EFOM), anaerobically degradable energy (adE) and total energy (tE). The minimum (Min), maximum (Max) and median values were calculated for each parameter. The median of BP fitting Case 2 were also calculated. BP 12 is a research biogas plant with high measurement accuracy.
Table 2. Yield efficiency (YE) of 33 biogas plants (BP) determined by literature values of Association for Technology and Structures in Agriculture (KTBL), biochemical methane potential test (BMP), fermentable organic matter (FOM), energy fermentable organic matter (EFOM), anaerobically degradable energy (adE) and total energy (tE). The minimum (Min), maximum (Max) and median values were calculated for each parameter. The median of BP fitting Case 2 were also calculated. BP 12 is a research biogas plant with high measurement accuracy.
Yield Efficiency [%]
BPBMPKTBLFOMEFOMadEtE
Case 11987673674236
21019798926153
3927977765045
486761041075449
5928796956759
6889887757550
71089998916963
8778285856055
9106971101106862
10909992857572
Case 2111059799956962
1210410195867166
131141031091016961
141141041031008072
151011171131078783
16961131061028581
171101161101038683
181221161321328982
191071091121067770
201101021151076959
21981161101038882
221151131101058275
231141101281228476
24961091091037570
25991131081028274
261381151241097059
27941161111048782
281181111151087870
Case 3291091131091028576
301001161161118274
311221241281228880
321321291261167968
331231161161098679
Min777673674236
Max1381291321328983
Median1051091091037770
Median Case 21091101121057973
Table 3. Conversion efficiency (CE) of 33 biogas plants (BP) determined by biochemical methane potential test (BMP), fermentable organic matter (FOM), energy fermentable organic matter (EFOM), anaerobically degradable energy (adE), total energy (tE) and organic dry matter content (ODM). The minimum (Min), maximum (Max) and median value were calculated for each parameter. The median of BP fitting Case 2 were also calculated. BP 12 is a research biogas plant with high measurement accuracy.
Table 3. Conversion efficiency (CE) of 33 biogas plants (BP) determined by biochemical methane potential test (BMP), fermentable organic matter (FOM), energy fermentable organic matter (EFOM), anaerobically degradable energy (adE), total energy (tE) and organic dry matter content (ODM). The minimum (Min), maximum (Max) and median value were calculated for each parameter. The median of BP fitting Case 2 were also calculated. BP 12 is a research biogas plant with high measurement accuracy.
Conversion Efficiency [%]
BPBMPFOMEFOMadEtEODM
Case 1189101123776669
290104121807073
39580104686164
496133159807275
597123120847579
6948680805364
79691117888284
89399107756972
998121144898285
1099100104928890
Case 21197101111817376
129491101847780
1393111120817275
149695105847680
1597111112918789
1697106105878385
1798106109918990
1898126136928588
1985105112827478
2095106112726165
2199108106918588
2297108112878083
2396116125867881
2489106108797375
2596100104847680
269196120776467
2798108112938991
2893106109787073
Case 32997104100837579
309397105776973
3198115117857780
3295104112766669
3395101106847780
Min858080685364
Max99133159938991
Median96105112847579
Median Case 295106112857780
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hülsemann, B.; Mächtig, T.; Pohl, M.; Liebetrau, J.; Müller, J.; Hartung, E.; Oechsner, H. Comparison of Biological Efficiency Assessment Methods and Their Application to Full-Scale Biogas Plants. Energies 2021, 14, 2381. https://doi.org/10.3390/en14092381

AMA Style

Hülsemann B, Mächtig T, Pohl M, Liebetrau J, Müller J, Hartung E, Oechsner H. Comparison of Biological Efficiency Assessment Methods and Their Application to Full-Scale Biogas Plants. Energies. 2021; 14(9):2381. https://doi.org/10.3390/en14092381

Chicago/Turabian Style

Hülsemann, Benedikt, Torsten Mächtig, Marcel Pohl, Jan Liebetrau, Joachim Müller, Eberhard Hartung, and Hans Oechsner. 2021. "Comparison of Biological Efficiency Assessment Methods and Their Application to Full-Scale Biogas Plants" Energies 14, no. 9: 2381. https://doi.org/10.3390/en14092381

APA Style

Hülsemann, B., Mächtig, T., Pohl, M., Liebetrau, J., Müller, J., Hartung, E., & Oechsner, H. (2021). Comparison of Biological Efficiency Assessment Methods and Their Application to Full-Scale Biogas Plants. Energies, 14(9), 2381. https://doi.org/10.3390/en14092381

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop