Next Article in Journal
Crude Oil Bioremediation: From Bacteria to Microalgae
Next Article in Special Issue
Anaerobic Co-Digestion of Sewage Sludge and Trade Wastes: Beneficial and Inhibitory Effects of Individual Constituents
Previous Article in Journal
A Technique for Transformer Remnant Cellulose Life Cycle Prediction Using Adaptive Neuro-Fuzzy Inference System
Previous Article in Special Issue
Toward the Transition of Agricultural Anaerobic Digesters into Multiproduct Biorefineries
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Shortening the Standard Testing Time for Residual Biogas Potential (RBP) Tests Using Biogas Yield Models and Substrate Physicochemical Characteristics

1
School of Water, Energy and Environment, Cranfield University, Cranfield MK43 0AL, UK
2
School of Aerospace, Transport and Manufacturing, Cranfield University, Cranfield MK43 0AL, UK
*
Author to whom correspondence should be addressed.
Processes 2023, 11(2), 441; https://doi.org/10.3390/pr11020441
Submission received: 16 December 2022 / Revised: 26 January 2023 / Accepted: 27 January 2023 / Published: 1 February 2023
(This article belongs to the Special Issue New Frontiers in Anaerobic Digestion (AD) Processes)

Abstract

:
The residual biogas potential (RBP) test is a procedure to ensure the anaerobic digestion process performance and digestate stability. Standard protocols for RBP require a significant time for sample preparation, characterisation and testing of the rig setup followed by batch experiments of a minimum of 28 days. To reduce the experimental time to obtain the RBP result, four biogas kinetic models were evaluated for their strength of fit for biogas production data from RBP tests. It was found that the pseudo-parallel first-order model and the first-order autoregressive (AR (1)) model provide a high strength of fit and can predict the RBP result with good accuracy (absolute percentage errors < 10%) using experimental biogas production data of 15 days. Multivariate regression with decision trees (DTs) was adopted in this study to predict model parameters for the AR (1) model from substrate physicochemical parameters. The mean absolute percentage error (MAPE) of the predicted AR (1) model coefficients, the constants and the RBP test results at day 28 across DTs with 20 training set samples are 4.76%, 72.04% and 52.13%, respectively. Using five additional data points to perform the leave-one-out cross-validation method, the MAPEs decreased to 4.31%, 59.29% and 45.62%. This indicates that the prediction accuracy of DTs can be further improved with a larger training dataset. A Gaussian Process Regressor was guided by the DT-predicted AR (1) model to provide probability distribution information for the biogas yield prediction.

1. Introduction

Digestate is a by-product from the anaerobic digestion (AD) process. Due to its high nutrient value, it can be used as soil improver or fertiliser if the digestate is proved to be valorised and can meet relevant quality standards [1]. The digestate stability can be evaluated with a residual biogas potential (RBP) test. The test typically is required to be carried out under mesophilic conditions for at least 28 days with an appropriate inoculum-to-substrate ratio and micro- and macronutrients supplemented to avoid the inhibition of biogas production [2]. The digestate is considered to have consistent quality if the RBP test biogas yield is below 0.25 L/g volatile solids (VS), as recommended in the Publicly Available Specification 110 (PAS110), which is a key element of the UK Government’s anaerobic digestion quality protocol [3].
The 28-day continuous monitoring of biogas production in an RBP test is time-consuming and onerous for commercial AD operators. This limits the adoption of RBP tests and regulated markets for digestate. There have been many attempts to find alternative approaches to RBP that offer rapid tests result. These include assessing acid production after the inhibition of methanogenesis [4,5,6] and assessing the digestion of the organic fraction of the digestate after separation of the microbial cell component [7]. Nevertheless, both approaches are of great complexity and further research is needed [7].
Additionally, some researchers have attempted to relate the RBP test results to digestate physicochemical characteristics. For example, the theoretical biogas potential calculated based on the stoichiometric methane conversion from volatile fatty acid (VFA) concentrations and the soluble chemical oxygen demand (sCOD) in the digestate sample were suggested in PAS110 as preliminary pass/fail indicators for the RBP tests [8]. Although these preliminary indicators provide useful information about digestate stability, they cannot be correlated with RBP values or provide reaction kinetics that can be used to reveal further digestion performance information including inhibition. In a previous work reviewing the application of the RBP test for PAS110 [9], correlations between RBP values and various characteristics, including VFA and sCOD, were investigated. Although some low to moderate levels of correlations were found for total VFA, total solids and volatile solids, which account for 40%, 36% and 29% of variation in RBP values, respectively, none of the indicators are sufficiently reliable to predict the RBP values accurately [9].
Other researchers have evaluated using empirical biogas production models, including first-order kinetic and Gompertz models, and experimental biogas yield data from the initial stage of the RBP tests to fit specific accumulative biogas production data from RBP tests [10,11,12]. This has led to a promising experimental and modelling ‘hybrid’ approach using experimental data collected from a shorter RBP duration (3–7 days instead of 28 days) to calculate model parameters, and then predict the ultimate biogas production. However, the accuracy of prediction is not sufficient to warrant the replacement of RBP with this hybrid approach [12]; therefore, further improvement of the modelling process is required.
In this research, we evaluated the strength of fit for four biogas yield kinetic models including first-order kinetic, modified Gompertz, pseudo-parallel first-order kinetic and autoregressive (AR) time-series to describe the RBP test biogas production process. The models are then calibrated using experimental data collected from shorter RBP tests (5, 10, 15, 20 and 25 days) to calculate model parameters that are then used to predict ultimate biogas production.
In a previous work [9], although using conventional statistical methods, no significant correlations were found between key physicochemical parameters of digestate samples with RBP results. Due to the potential interplay of these parameters, which can influence the RBP results and reaction kinetics, the correlations may be deeply hidden.
Machine learning techniques including multivariate nonlinear regression analysis with decision trees (DT) were applied to predict the parameters of the biogas production model from the physicochemical characteristics of digestate samples. Compared with other multivariate nonlinear regression methods, the DT method is particularly suitable for a training dataset with limited sample size in this study [13]. The uncertainties of the predicted biogas yield were then assessed using a Gaussian process regressor (GPR).
The data processing framework described in this work can potentially have a wider application for other complex biochemical processes that are influenced by multiple physicochemical parameters of the reaction system.

2. Materials and Methods

2.1. Digestate Samples and RBP Test

The sampling point of the 25 digestate samples was the outlet of the final tank from which the biogas was collected. The AD plants involved in the study are anonymised and coded as ADP1–25. The inoculum was from the anaerobic digester at Millbrook Wastewater (WW) Treatment Plant at Southampton, UK. The RBP test followed the standard procedure described in the PAS110 [3]. Samples were tested in triplicate against two positive controls and three inoculum-only controls.

2.2. Analytical Methods

Twenty physicochemical characteristics were analysed for each digestate, including VFAs (total VFA and acetate), total ammoniacal nitrogen (TAN), total Kjeldahl nitrogen (TKN), alkalinity (total alkalinity (TA), partial alkalinity (PA) and intermediate alkalinity (IA)), TS, VS, pH, COD (total COD and sCOD), trace element (TE) concentrations (cobalt (Co), iron (Fe), molybdenum (Mo) and nickel (Ni)), calorific value (CV) and elemental compositions (C, H, N). TAN/TKN and IA/PA were calculated as two extra metrics. TAN/TKN represents the relative contents of ammonia nitrogen and organically bonded nitrogen and thus how ready the substrate is for microorganism degradation [14]. IA/PA is an indicator of VFA accumulation in the AD process. Ripley et al. (1986) [15] suggest IA/PA < 0.3 indicates a stable state of the anaerobic process.
Determination of VFAs is based on the SCA (1979) [16] method Determination of Volatile Fatty Acids in Sewage Sludge. Supernatant layer from digestate centrifugation with 10% formic acid were quantified in a Shimadzu GC-2010 gas chromatograph with a flame ionisation detector and a capillary column type SGE BP-21. TAN and TKN were determined using a Kjeltech digestion block and steam distillation unit, according to the manufacturer’s instructions (Foss Ltd., Warrington, UK). Alkalinity was measured by titration with 0.25 N H2SO4 to endpoints of pH 5.75 and 4.3 in order to allow calculation of TA, PA and IA [15]. TS and VS were determined with Standard Method 2540 G (APHA, 2005). Total COD and sCOD were analysed by adapting the closed reflux titrimetric method of 5220C, APHA [17]. TE concentrations were determined using ICP-MS or ICP-OES at a UKAS-accredited commercial laboratory (Severn Trent Services, Coventry, UK) after in-house hydrochloric–nitric acid digestion [18]. CV was measured with a CAL2k-ECO bomb calorimeter (CAL2k, Digital Data Systems, Gauteng, South Africa). Elemental C, H, and N analysis was performed using a Flash EA-1112 elemental analyser (Thermo Finnigan, Cheshire, UK).

2.3. Assessing the Strength of Fit for Biogas Production Kinetic Models

Four biogas-production kinetic models commonly used to estimate the kinetic constants of the AD process were compared for their abilities to fit the biogas production of RBP tests. These models include three empirical models (first-order kinetic, modified Gomperz and pseudo-parallel first-order models) and a time-series model (first-order autoregressive). The strength of a kinetic model to accurately fit the experimental biogas production data was measured using R2 values, which indicate the percentage of the variance in the responses explained by a model.
Additionally, these models were used to predict the 28-day RBP test results by fitting the models with initial 5, 10, 15, 20 and 25 days’ experimental data using the Matlab R2021b Curve Fitting Toolbox. Based on the absolute percentage error (APE) between the experimental data and model-predicted results, the accuracy of prediction and duration of experimental data required to obtain a sufficiently accurate prediction were investigated for the following four models:
(1) First-order model (FO): The FO model (Equation (1)) is derived from the assumptions that the substrate degradation is a first-order reaction with hydrolysis as the speed-limiting step and the cumulative biogas yield is proportional to the amount of substrate degraded (Equation (2)) [19,20].
y t = y m 1 e k t  
d c d t = k c   c 0 c c 0 = y y m  
where y(t) is the cumulative biogas yield at time t, k is the first-order rate constant, ym is the maximum cumulative gas production, c is the concentrate of the substrate and c0 is the initial substrate concentration.
(2) Modified Gomperz model (MG): The modified Gomperz model (Equation (3)) is derived from the Gomperz model, which is used to describe the microbial activity and has a signature sigmoid shape [21,22]. It describes the biogas production in terms of the exponential growth rates and lag phase duration of anaerobic degradation microorganisms [11].
y t = y m × exp exp R y m × e × λ t + 1
where y(t) is the cumulative gas production at time t, ym is the maximum cumulative gas production (mL CH4/gVS), R is the maximum gas production rate (mL CH4/gVS/d) and λ is the lag phase period or minimum time to produce biogas (days).
(3) Pseudo-parallel first-order model (PP): The pseudo-parallel first-order model (Equation (4)) is considered to be more suitable for describing the biogas yield of mixtures of substrates with different kinetic rates (rapid and slow) [23].
y t = y m 1 P e k 1 t 1 P e k 2 t
where y(t) is the cumulative gas production at time t, ym is the maximum cumulative gas production (mL CH4/g VS), P is the the proportion of the readily degradable material, k1 is the first-order rate constant for readily degradable material, and k2 is the first-order rate constant for less readily degradable material.
(4) First-order autoregressive model (AR (1)): AR (1) is a time-series model that predicts the present timestep based on the observations from previous timesteps. The autocorrelation function (ACF) (autocorrelation between timesteps) plots for all the RBP test biogas yield samples gradually trail off (Figure 1a, using ADP20 as an example). Therefore, the biogas production process is an AR process. Many time-series models that essentially model the randomness of the time series data need the trends in the data to be removed, in other words, to ensure the stationarity of the data. However, the application of the AR model does not intrinsically require transforming the data into stationary data. Thus, the biogas yield data were not converted to a stationary process in this study.
For most of the RBP test data, the cumulative biogas production curves typically exhibit a rapid gas-production stage in the initial three days followed by a noticeable reduction in biogas production rate. This means one set of AR model parameters generally will result in a poor fitting for both the rapid-production stage and the later stage after the initial three days. In this study, to simplify the modelling process, the AR (1) model was only used to model biogas production from day 4 of the RBP tests. The partial autocorrelation function (PACF) (partial autocorrelation between timesteps) plots of data from day four for most of the biogas time series samples were cut off after one lag (Figure 1b, taking ADP20 as an example). Therefore, the biogas production from day four was an AR (1) process (Equation (5)).
The upper asymptote of an AR (1) process is determined by its unconditional mean μ = c/(1 − β), which corresponds to the ultimate biogas yield of a substrate.
X t = β X t 1 + c + ε t
where Xt is the response of the present timestep, the Xt−1 is the response of the previous one timestep, β is the coefficient, c is the constant and εt is the white noise.

2.4. Prediction of Biogas Kinetic Model Coefficients and Constants Using Decision Tree Multivariate Regression Method

2.4.1. Decision Tree Multivariate Regression Method

The best-performing model out of the four kinetic models was further studied using decision tree (DT) multivariate regression analysis. DT is a non-parametric supervised learning technique that can be applied for both regression analysis and classification, providing piecewise constant approximations as prediction results [24].
DTs were trained using all physicochemical characteristics listed in Table 1 and Table 2 as predictors to predict the coefficients and constants of the best-performing biogas production kinetic model. First, RBP results of 20 digestate samples (out of the 25-sample dataset) were used to train the DT model with the five remaining samples as the test dataset for model validation. The training process generated 53,130 groups of different combinations of 20 samples when choosing out of the 25. For practical reasons, a subset of 5000 groups were selected randomly to evaluate the absolute percentage error (APE) and mean absolute percentage error (MAPE) (Equation (6)) of the predicted biogas production model parameters and the APE of the calculated ultimate biogas yield. Then, the DT model was trained and validated using the leave-one-out (LOO) cross-validation method, which apply four additional set of data to the training data to verify if the prediction accuracy would improve.
  A P E i = y i y i ^ y i   M A P E = 1 n i = 1 n y i y i ^ y i
where yi is the model parameter derived by fitting the model to the entire 28-day RBP test data, y i ^ is the model parameter inferred from the physicochemical characteristics by DT and n is the number of training set samples.

2.4.2. Assessing Prediction Uncertainty Using a Gaussian Process Regressor (GPR)

To quantify the uncertainty in the calculated biogas yield using the kinetic model predicted by DT, a Gaussian process regressor (GPR) method was applied in this study. GPR is a kernel-based Bayesian tool to perform nonlinear regression. The process is specified by a mean function m(x) and a covariance function k (x, x’|θ) which defines the covariances between the responses at any two input locations x and x’. θ represents the hyperparameters of the covariance function and their values are learned from the training data by maximising the log marginal likelihood. Once the hyperparameters are decided, the prediction for new input is performed by computing the marginal posterior distribution conditioning on the dimensions with known inputs [25].
The principle of GPR is to predict one timestep further each time by fitting one more datum provided by the biogas-yield kinetic model and the prediction uncertainty bands were also returned. The predictions of the GPR with zero mean and squared exponential kernel and the GPR with linear basis function and squared exponential kernel were compared.

3. Results and Discussion

3.1. Digestate Characterisation and RBP Test Results

The 28-day RBP test results for each set of RBP samples together with the physico-chemical characteristics of each digestate are shown in Table 1 and Table 2. Specific cumulative biogas yield data collected during the 28-day testing period were reported elsewhere [9] and were used in the model fitting and training process in this study.

3.2. Assessing Strength of Fit for Biogas Production Models

Within the 28 days of the RBP test, it was noticeable that biogas yields of some samples in this study had reached a plateau, whilst others still were showing an upward trend close to the end of the 28-day test. This is clearly due to the different concentrations of readily degradable materials in the digestate samples. To distinguish these two types of digestate samples, the 25 digestate samples were classified into two types based on the absolute average of the daily biogas production change percentages in the last four days: (1) Type I: less than 0.5%; (2) Type II: more than 0.5% (Table 3).
Table 4 shows the R2 values for the fits of the three empirical biogas production models to RBP test data. Overall, the modified Gomperz model achieved lower R2 values across the majority of the samples, indicating poorer fitting performance. In contrast, the PP first-order model could describe the biogas production of almost all samples with an R2 value between 97–99%. The FO model performed better at describing Type I samples, whereas the PP first-order model was more suitable for substrates mixed with materials with different reaction-rate constants, and therefore performed better with Type II samples. However, it is worth noting that there are usually multiple sets of optimal solutions of the estimates of the PP first-order model’s parameters (Ym, P, k1 and k2). This is because the nonlinear least-square error function of this model when fitting a particular set of biogas yield data is not always convex. Therefore, it was not chosen for the study of training DTs to predict the model parameters from the digestate physicochemical characteristics in the following section.
In addition, when fitting the PP first-order model in Matlab using the least-square algorithm, the initial value set for the parameter P should avoid 0.5 and k1 and k2 should not be the same. Otherwise, the partial derivatives of the error function with respect to k1 and k2 are the same, which means the moving direction of these two dimensions are the same and the nonlinear search for the minimum value of the error function value will settle at a local minimum point.
The fitting of the AR (1) model after the initial 3-day rapid biogas production stage was comparable to the PP first-order model, with R2 values of 99% for the majority of samples (Table 5), regardless of Type I or Type II data. Figure 2 shows the fits of three empirical models and the AR (1) model to the cumulative biogas yield in RBP tests.
When using these four models to predict the final RBP results based on experimental data collected from a shorter test duration (5, 10, 15, 20 and 25 days), none of the three empirical models and the AR (1) model could achieve a sufficient level of accuracy to replace full-length RBP tests. This conclusion is in agreement with previous studies [12,26]. When using the first-order model, typically the level of prediction accuracy increases when more experimental data points are provided. FO performed moderately better than the MG model at predicting the RBP test result. This result was in agreement with the result of Nielfa et al. (2015) [26] for biomethane potential test biogas production prediction for the organic fraction of municipal solid waste and biological sludge co-digestion. However, both FO and MG models were proved to be unsuitable to predict the RBP result of Type II data.
Using 15 days of experimental data, the PP model could predict the RBP result with APE < 10% in nearly all the samples (both Type I and II). Therefore, the PP model is a preferred option for RBP test-result-prediction when the experiment is half-way through. For the AR (1) model, RBP data fitting and modelling starts from day 4, thus 10, 15, 20 and 25 days of experimental data were used for model fitting. The prediction ability of the AR (1) model was comparable to that of the PP model. Table 6 shows the predictions of RBP test results when fitting an increasing amount of experimental data to four biogas production models (due to the large sample numbers, only seven samples were randomly selected from each group of Type I and II data).

3.3. Biogas Yield Prediction from Digestate Physicochemical Characteristics by DT

DTs were first trained with the physicochemical characteristics of 20 digestate samples as predictors and the fitted coefficients and constants of AR (1) models as responses. A total of 5000 groups of splits between the training set and test set were randomly chosen to evaluate the prediction accuracy. The average MAPEs of the predicted AR (1) model coefficients and constants for among 5000 groups of test sets were 4.58% and 72.04%, respectively. The MAPE of the calculated RBP test result at day 28 from the predicted AR (1) model parameters was 52.125%.
With four more samples provided for the training set, the MAPEs of the predictions for the AR (1) model coefficient and constant with the LOO cross-validation method were 4.31% and 59.29%, respectively (Table 7). The MAPE of the calculated RBP test result at day 28 was 45.620%. With four more data provided, the prediction accuracy of the RBP test result of the DTs improved 12.48%. Additionally, among 25 LOO cross-validation groups, the APEs of predicted RBP test results of five groups were smaller than 10%.

3.4. AR (1) Model Prediction Guide GPR

When training the DTs with LOO cross-validation and ADP3 as the test set, the APE of the predicted biogas yield on day 28 using the inferred AR (1) model was 6.87%. This is used as an example for illustration in Figure 3. By accepting one more input each time from the inferred AR (1) model together with the experimental observations of the first three days, the predictions of every one timestep further of the GPR with zero mean and squared exponential kernel were smaller than AR (1)’s predictions for the first few timesteps and then gradually closer to the predictions of the AR (1) model. This was explained by the increased fitted GPR model length scale and vertical scale when receiving more training data. The smaller the length scale, the curvier the underlying function is, and the smaller the vertical scale, the more concentrated the underlying function is around the mean. In contrast, for the GPR with a linear basis function, the predictions surged for the first few timesteps and then approached the AR (1) model predictions. This corresponded to the decreased slope of the fitted linear mean function. In general, given that only one more timestep is predicted at a time, the selection of zero or linear mean function is not of much concern. The 95% confidence interval of the prediction of GPR narrowed when more data were provided.

4. Conclusions

This study has demonstrated that it is possible to use RBP experimental data collected in the initial stage of the test to predict the 28-day RBP result in a kinetic model fitting exercise. By fitting 15 days of experimental data from RBP tests to kinetic models, the PP first-order model and the AR (1) model achieved a promising accuracy with APE < 10%.
Further study demonstrated using the decision tree (DT) method that AR (1) model parameters can be predicted from the physicochemical characteristics of the digestate samples. This provides potential to further reduce the data requirement to four days of RBP experimental data and thereby significantly reduce the resting time of a standard 28-day RBP test to around four days. It was observed that when more training data were included in the DT machine learning model (from 20 to 24 samples), the prediction accuracy of the RBP result increased by 12.48%. This indicates that collecting more data to include in the model-training process can further improve the prediction outcome.
The framework of predicting kinetic model parameters from the physicochemical characteristics of the substrate can potentially be applied to the yield prediction of the product from other biochemical reaction processes.

Author Contributions

Conceptualization, W.G., P.L. and Y.J.; Methodology, Y.L., W.G. and Y.J.; Formal analysis, Y.L.; Writing—original draft, Y.L.; Writing—review & editing, W.G., P.L. and Y.J.; Supervision, W.G., P.L. and Y.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analysed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lee, M.E.; Steiman, M.W.; St Angelo, S.K. Biogas digestate as a renewable fertilizer: Effects of digestate application on crop growth and nutrient composition. Renew. Agric. Food Syst. 2021, 36, 173–181. [Google Scholar] [CrossRef]
  2. Walker, M.; Banks, C.; Heaven, S.; Frederickson, J. OFW004-005. Waste and Resources Action Programme, Residual Biogas Potential Test for Digestates; University of Southampton & Open University: Southampton, UK, 2010. [Google Scholar]
  3. WRAP. PAS 110:2014. Specification for Whole Digestate, Separated Liquor and separated Fibre Derived from the Anaerobic Digestion of Source-Segregated Biodegradable Materials. 2014. Available online: https://wrap.org.uk/sites/default/files/2021-03/PAS110_2014.pdf (accessed on 26 June 2022).
  4. Lie, E.; Welander, T. A method for determination of the readily fermentable organic fraction in municipal wastewater. Water Res. 1997, 31, 1269–1274. [Google Scholar] [CrossRef]
  5. Ruel, S.M.; Comeau, Y.; Héduit, A.; Deronzier, G.; Ginestet, P.; Audic, J.M. Operating conditions for the determination of the biochemical acidogenic potential of wastewater. Water Res. 2002, 36, 2337–2341. [Google Scholar] [CrossRef] [PubMed]
  6. Mohan, S.V.; Babu, V.L.; Sarma, P.N. Effect of various pretreatment methods on anaerobic mixed microflora to enhance biohydrogen production utilizing dairy wastewater as substrate. Bioresour. Technol. 2008, 99, 59–67. [Google Scholar] [CrossRef] [PubMed]
  7. Banks, C.J.; Heaven, S.; Zhang, Y.; Sapp, M.; Thwaites, R. OMK002-014. Waste and Resources Action Programme. In Review of the Application of the Residual Biogas Potential Test; University of Southampton: Southampton, UK, 2013. [Google Scholar]
  8. WRAP. OFW004-005. Residual Biogas Potential Test for Digestates, Development and Evaluation of a Method for Testing the Residual Biogas Potential of Digestates (Section 4, Full Description of the RBP Test). 2010. Available online: https://www.ktbl.de/fileadmin/user_upload/Allgemeines/Download/Ringversuch-Biogas/Residual-Biogas-Potential.pdf (accessed on 20 July 2022).
  9. WRAP. OMK002-014. Review of the Application of the Residual Biogas Potential Test. 2013. Available online: http://www.organics-recycling.org.uk/uploads/article2652/PAS110%20digestate%20stability%20review.pdf (accessed on 13 July 2022).
  10. Strömberg, S.; Nistor, M.; Liu, J. Early prediction of biochemical methane potential through statistical and kinetic modelling of initial gas production. Bioresour. Technol. 2015, 176, 233–241. [Google Scholar] [CrossRef] [PubMed]
  11. Pramanik, S.K.; Suja, F.B.; Porhemmat, M.; Pramanik, B.K. Performance and kinetic model of a single-stage anaerobic digestion system operated at different successive operating stages for the treatment of food waste. Processes 2019, 7, 600. [Google Scholar] [CrossRef]
  12. Calabrò, P.S.; Folino, A.; Maesano, M.; Pangallo, D.; Zema, D.A. Exploring the Possibility to Shorten the Duration and Reduce the Number of Replicates in Biomethane Potential Tests (BMP). Waste Biomass Valorization 2022, 2022, 939. [Google Scholar] [CrossRef]
  13. Hastie, T.; Tibshirani, R.; Friedman, J.H. Chapter 9: Additive Models, Trees, and Related Methods. In The Elements of Statistical Learning: Data Mining, Inference, and Prediction; Springer: New York, NY, USA, 2009; Volume 2, pp. 1–758. [Google Scholar]
  14. Deng, Z.; van Linden, N.; Guillen, E.; Spanjers, H.; van Lier, J.B. Recovery and applications of ammoniacal nitrogen from nitrogen-loaded residual streams: A review. J. Environ. Manag. 2021, 295, 113096. [Google Scholar] [CrossRef] [PubMed]
  15. Ripley, L.E.; Boyle, W.C.; Converse, J.C. Improved alkalimetric monitoring for anaerobic digestion of high strength wastes. J. Water Pollut. Control Fed. 1986, 58, 406–411. [Google Scholar]
  16. SCA. Determination of volatile fatty acids in sewage sludge. In Methods for the Examination of Waters and Associated Materials; Standing Committee of Analysts; HMSO: London, UK, 1979. [Google Scholar]
  17. APHA. Standard Method for Examination of Water and Wastewater, 21st ed.; AWWA, WPCF: Washington, DC, USA, 2005. [Google Scholar]
  18. SCA. Methods for the determination of metals in soils, sediments and sewage sludge and plants by hydrochloric-nitric acid digestion. In Methods for the Examination of Waters and Associated Materials; Standing Committee of Analysts; HMSO: London, UK, 1986. [Google Scholar]
  19. Zahan, Z.; Othman, M.Z.; Muster, T.H. Anaerobic digestion/co-digestion kinetic potentials of different agro-industrial wastes: A comparative batch study for C/N optimisation. Waste Manag. 2018, 71, 663–674. [Google Scholar] [CrossRef] [PubMed]
  20. Kafle, G.K.; Chen, L. Comparison on batch anaerobic digestion of five different livestock manures and prediction of biochemical methane potential (BMP) using different statistical models. Waste Manag. 2016, 48, 492–502. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Gompertz, B. XXIV. On the nature of the function expressive of the law of human mortality, and on a new mode of determining the value of life contingencies. In a letter to Francis Baily, Esq. F. R. S. &c. Philos. Trans. R. Soc. Lond. 1825, 115, 513–583. [Google Scholar]
  22. Zhu, H.; Yang, J.; Cheng, X. Application of Modified Gompertz Model to Study on Biogas production from middle temperature co-digestion of pig manure and dead pigs. E3S Web Conf. 2019, 118, 1–7. Available online: https://www.e3s-conferences.org/articles/e3sconf/pdf/2019/44/e3sconf_icaeer18_03022.pdf (accessed on 10 October 2021). [CrossRef]
  23. Pererva, Y.; Miller, C.D.; Sims, R.C. Existing empirical kinetic models in biochemical methane potential (BMP) testing, their selection and numerical solution. Water 2020, 12, 1831. [Google Scholar] [CrossRef]
  24. Breiman, L.; Friedman, J.H.; Olshen, R.A.; Stone, C.J. Regression. In Classification and Regression Trees, 1st ed.; Routledge & CSC Press: Boca Raton, FL, USA, 2017; pp. 216–264. [Google Scholar]
  25. Rasmussen, C.E.; Williams, C.K.I. Regression. In Gaussian Processes for Machine Learning; Dietterich, T., Bishop, C., Heckerman, D., Jordan, M., Kearns, M., Eds.; MIT Press: Cambridge, MA, USA, 2006; pp. 7–29. [Google Scholar]
  26. Nielfa, A.; Cano, R.; Fdz-Polanco, M. Theoretical methane production generated by the co-digestion of organic fraction municipal solid waste and biological sludge. Biotechnol. Rep. 2015, 5, 14–21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. The ACF (a) and PACF (b) of sample ADP20 biogas production time series from day 4. The blue lines are the confidence bounds for a significant correlation. The red lines are the autocorrelation and partial autocorrelation between the current time step and the time steps of different lags.
Figure 1. The ACF (a) and PACF (b) of sample ADP20 biogas production time series from day 4. The blue lines are the confidence bounds for a significant correlation. The red lines are the autocorrelation and partial autocorrelation between the current time step and the time steps of different lags.
Processes 11 00441 g001
Figure 2. The fits of the three empirical models and the AR (1) model to the entire RBP test process for Type I (a) and Type II (b) data.
Figure 2. The fits of the three empirical models and the AR (1) model to the entire RBP test process for Type I (a) and Type II (b) data.
Processes 11 00441 g002
Figure 3. The prediction and its probability distribution of the GRP guided by the prediction of the inferred AR (1) model. (a) Zero mean and squared exponential kernel. (b) Linear basis function and squared exponential kernel.
Figure 3. The prediction and its probability distribution of the GRP guided by the prediction of the inferred AR (1) model. (a) Zero mean and squared exponential kernel. (b) Linear basis function and squared exponential kernel.
Processes 11 00441 g003
Table 1. RBP test results on day 28 (d-28) and the physicochemical characteristics of 25 samples (Part 1).
Table 1. RBP test results on day 28 (d-28) and the physicochemical characteristics of 25 samples (Part 1).
ADP No.d-28AcetateTotal-VFATANTKNTAN/TKNTAPAIAIA/PATSVS
L/g VSmg/Lg N/kg kg/kg CaCO3g/kg
10.13178.18193.571.703.530.4812.529.033.490.3949.8437.12
20.1829.0837.350.623.280.197.976.421.550.2449.1130.39
30.13846.781016.417.9812.370.6433.3825.877.510.2993.2867.65
40.24684.05775.994.046.480.6217.5212.425.100.4146.4332.59
50.12247.38332.575.127.480.6924.5016.627.880.4756.5737.49
60.06111.56115.642.804.910.5746.4328.1218.310.65178.3450.05
70.07183.61187.392.694.440.6158.2020.9237.271.78137.4340.44
80.2614.42324.320.401.290.313.312.201.110.5017.1612.11
90.362633.799262.766.549.230.7123.9915.628.360.5447.7934.51
100.17134.49204.903.325.210.6417.6413.104.540.3546.5329.51
110.302662.153963.272.223.140.719.455.224.230.8120.3812.74
120.1723.2936.682.714.640.5812.228.983.240.3636.9126.50
130.16335.54364.462.805.290.5315.7312.143.590.3043.5929.98
140.1319.2519.250.442.250.203.971.852.111.1435.0421.47
150.09250.28259.053.586.120.5817.2013.663.540.2658.0043.56
160.382706.803871.623.004.870.6211.917.764.140.5336.6526.88
170.2636.6950.291.492.850.527.715.791.910.3329.8216.42
180.25832.241440.642.424.270.5713.709.983.720.3759.0444.97
190.1439.8253.934.626.850.6821.1316.314.820.3052.9238.76
200.19184.69218.065.768.550.6723.5618.345.230.2964.6444.93
210.3389.69260.153.005.610.5424.8514.7010.140.69211.53106.18
220.22299.39353.292.253.110.7212.7810.072.710.2720.869.85
230.29270.61411.154.246.560.6517.2313.243.990.3052.4932.30
240.28215.32250.624.476.980.6418.1913.874.330.3150.2534.87
250.29241.71272.534.777.010.6819.5215.154.380.2946.8431.96
Table 2. RBP test results on day 28 (d-28) and the physicochemical characteristics of 25 samples (Part 2).
Table 2. RBP test results on day 28 (d-28) and the physicochemical characteristics of 25 samples (Part 2).
ADP No.d-28pHCVCHNCoFeMoNiTotal-CODsCOD
L/g VSMJ/kg TS%mg/Lg O2/L
10.138.3016.3839.204.485.610.85130.530.291.2042.776.67
20.187.4114.8937.444.575.530.081144.150.300.5040.352.73
30.138.3518.6134.804.527.340.241020.820.431.7588.6421.56
40.248.1717.5340.134.896.610.282031.830.190.6044.0412.14
50.128.4516.1338.324.206.050.2190.760.180.3940.6915.16
60.068.143.2917.081.712.201.082326.910.947.7457.2911.19
70.078.164.4112.511.592.250.981923.480.695.0639.328.45
80.267.3320.2646.255.586.110.0387.360.070.2428.623.16
90.368.3522.2647.755.328.510.09539.070.140.5075.3935.94
100.178.0414.7836.854.005.740.42174.880.300.7136.127.82
110.307.6217.0238.874.436.300.04210.670.080.1928.6910.00
120.178.1018.2642.434.886.130.25234.070.080.2831.6810.17
130.168.3716.2138.474.476.651.04111.120.301.0246.889.53
140.137.5015.2532.944.505.780.451166.461.461.6719.634.50
150.098.9019.8644.135.285.220.09151.820.251.0865.6311.87
160.387.9219.2243.495.626.540.0964.010.180.2639.1818.84
170.267.9213.3530.673.825.670.19348.020.280.5626.516.39
180.258.1617.7642.164.764.349.06227.305.6830.3755.7417.00
190.148.4217.9044.564.676.871.42293.850.501.3358.5814.74
200.198.1516.4138.554.556.401.211056.090.451.5568.5215.57
210.338.0712.2433.212.642.352.584249.891.1713.7982.9212.59
220.228.3210.7424.932.885.270.8936.940.260.8217.769.52
230.298.5413.8432.934.065.590.18649.780.160.4654.3316.54
240.288.5416.8238.594.826.650.38679.620.160.5867.3915.86
250.298.6716.5036.824.715.900.37797.930.150.5361.5616.83
Table 3. The classification of two types of RBP test biogas-production time series according to the absolute average of the daily biogas production change percentages in the last four days.
Table 3. The classification of two types of RBP test biogas-production time series according to the absolute average of the daily biogas production change percentages in the last four days.
Type IADP34511121316192021
Increase (%)0.380.040.260.270.360.310.180.360.40.50
Type IIADP126789101415171822232425
Increase (%)1.071.131.231.560.700.580.580.940.951.131.000.920.770.860.80
Table 4. R2 values of the fits of three empirical models (FO, MP and PP) to the Type I and Type II RBP test biogas-production time series (R2 values larger than 97% are in bold).
Table 4. R2 values of the fits of three empirical models (FO, MP and PP) to the Type I and Type II RBP test biogas-production time series (R2 values larger than 97% are in bold).
Type I Fitting R2 Values (%)
ADP34511121316192021
FO98.598.497.986.899.098.299.597.796.497.2
MG97.694.096.680.796.191.896.094.988.361.6
PP98.598.797.994.699.199.599.697.898.797.6
Type II Fitting R2 Values (%)
ADP126789101415171822232425
FO97.198.686.577.283.297.398.497.697.597.595.596.098.899.199.5
MG88.893.978.685.868.299.193.090.874.988.485.685.996.797.399.0
PP99.999.692.892.698.097.498.999.598.997.899.198.998.899.199.5
Table 5. R2 values of the fits of the AR (1) model to the Type I and Type II RBP test biogas-production time series from day 4 (R2 values larger than 97% are in bold).
Table 5. R2 values of the fits of the AR (1) model to the Type I and Type II RBP test biogas-production time series from day 4 (R2 values larger than 97% are in bold).
Type I Fitting R2 Values (%)
ADP34511121316192021
AR (1)98.299.497.497.098.699.399.899.798.990.1
Type II Fitting R2 Values (%)
ADP126789101415171822232425
AR (1)99.299.885.480.995.999.199.599.798.999.199.699.598.799.199.9
Table 6. The level of accuracy (expressed in APE%) of predicted 28-day RBP test result from fitting the mathematical models to an increasing number of experimental observations from 5 days to 25 days. APE less than 10% indicates accurate predication (in bold).
Table 6. The level of accuracy (expressed in APE%) of predicted 28-day RBP test result from fitting the mathematical models to an increasing number of experimental observations from 5 days to 25 days. APE less than 10% indicates accurate predication (in bold).
Type IDaysFOMGPPAR (1)Type IIDaysFOMGPPAR (1)
ADP 3577.1624.2335.08/ADP 2546.4651.9428.00/
104.989.627.261.151024.3437.7724.6512.20
151.015.142.700.471517.3424.486.864.52
202.101.550.881.712010.9814.872.021.63
252.290.073.231.81256.406.250.430.27
ADP 4557.9926.6552.67/ADP 853.6527.783.8080/
107.5015.327.193.971015.2417.5615.1312.57
155.005.534.656.001512.0312.028.318.62
201.841.714.262.63209.629.625.238.20
250.400.221.971.91253.483.442.765.74
ADP 55159.2517.17105.36/ADP 105163.3937.6229.58/
101.2111.801.048.22107.6325.7419.182.99
151.491.721.541.83157.7316.4112.750.03
200.800.861.512.33205.4010.2021.580.83
251.301.291.692.08253.496.224.851.02
ADP 13523.1243.0825.98/ADP 1450.1644.205.68/
1021.7427.407.710.561027.0432.5727.0644.65
1511.8512.9622.508.861517.0419.331.208.78
204.954.1910.033.702010.6710.335.623.70
252.501.872.912.58256.584.851.922.04
ADP 16519.9037.7117.89/ADP 1751.646445.2011.25/
104.6618.194.765.691011.9926.915.0247.42
153.388.203.050.571513.8617.950.840.25
201.092.653.701.372010.9011.124.123.06
250.610.951.000.89256.004.611.531.94
ADP 195183.8924.27146.92/ADP 1852.9845.763.74/
102.3313.072.338.621030.9434.0530.935.6969
153.137.293.411.311522.3622.197.510.57
202.783.232.790.362012.6411.804.161.37
251.961.791.910.66256.294.442.960.89
ADP 20524.1131.6227.78 ADP 25581.3654.3328.59/
1011.1621.668.6448.841036.3323.982.2414.84
158.9313.724.540.62158.5214.081.091.92
207.159.660.061.15203.267.757.050.48
255.326.691.120.51251.064.501.340.11
Table 7. The APEs for the AR (1) model coefficient and constant prediction and the biogas yield at 28th day prediction from the LOO cross-validation method.
Table 7. The APEs for the AR (1) model coefficient and constant prediction and the biogas yield at 28th day prediction from the LOO cross-validation method.
Test Sample No.Fitted AR (1) CoefficientPredicted AR (1) CoefficientAR (1) Coefficient APEFitted AR (1) ConstantPredicted AR (1) ConstantAR (1) Constant APEReal RBP Test ResultPredicted RBP Test ResultRBP Test Result APE
ADP10.9320.9370.5%0.0100.0101.1%0.1320.1374%
ADP20.9490.9341.6%0.0110.01965.6%0.1820.24333.5%
ADP30.7890.88412%0.0280.01740.2%0.1320.1416.9%
ADP40.8760.8962.2%0.0300.03517.5%0.2360.32136%
ADP50.8340.8501.9%0.0200.03684%0.1170.236101.9%
ADP60.9460.81314%0.0040.011174.7%0.0620.0577.9%
ADP70.9690.9323.8%0.0030.011293%0.0660.137107.4%
ADP80.9290.9370.9%0.0190.03684.6%0.2610.49088%
ADP90.8530.8934.8%0.0550.03241.2%0.3610.29119.5%
ADP100.8900.8940.4%0.0200.02316.1%0.1740.20618.6%
ADP110.9000.9414.6%0.0300.03514.6%0.3010.50668.4%
ADP120.8320.95014.1%0.0280.03523.5%0.1710.523205.7%
ADP130.9060.8921.5%0.0160.00943%0.1640.08647.6%
ADP140.9410.9360.6%0.0090.01014.6%0.1290.1387.1%
ADP150.9010.8476.1%0.0090.020123%0.0860.12646.5%
ADP160.8640.8422.5%0.0530.03238.7%0.3790.20546%
ADP170.9340.9360.2%0.0190.00952.8%0.2570.13945.9%
ADP180.9610.8709.4%0.0130.041211.6%0.2540.30721%
ADP190.8610.8984.2%0.0200.02316.5%0.1390.21252.5%
ADP200.8870.8940.8%0.0220.0221.8%0.1940.2024.3%
ADP210.8180.8929.2%0.0060.01075.2%0.3250.09271.6%
ADP220.9100.9393.2%0.0210.0227%0.2220.31541.9%
ADP230.8500.8460.5%0.0430.03321.4%0.2860.21524.7%
ADP240.8720.8962.7%0.0360.0345.2%0.2800.31312%
ADP250.9010.8476.0%0.0300.03515.3%0.2870.22421.7%
MAPE 4.3% 59.3% 45.6%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, Y.; Guo, W.; Longhurst, P.; Jiang, Y. Shortening the Standard Testing Time for Residual Biogas Potential (RBP) Tests Using Biogas Yield Models and Substrate Physicochemical Characteristics. Processes 2023, 11, 441. https://doi.org/10.3390/pr11020441

AMA Style

Liu Y, Guo W, Longhurst P, Jiang Y. Shortening the Standard Testing Time for Residual Biogas Potential (RBP) Tests Using Biogas Yield Models and Substrate Physicochemical Characteristics. Processes. 2023; 11(2):441. https://doi.org/10.3390/pr11020441

Chicago/Turabian Style

Liu, Yanxin, Weisi Guo, Philip Longhurst, and Ying Jiang. 2023. "Shortening the Standard Testing Time for Residual Biogas Potential (RBP) Tests Using Biogas Yield Models and Substrate Physicochemical Characteristics" Processes 11, no. 2: 441. https://doi.org/10.3390/pr11020441

APA Style

Liu, Y., Guo, W., Longhurst, P., & Jiang, Y. (2023). Shortening the Standard Testing Time for Residual Biogas Potential (RBP) Tests Using Biogas Yield Models and Substrate Physicochemical Characteristics. Processes, 11(2), 441. https://doi.org/10.3390/pr11020441

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