An assessment of ensemble learning approaches and single-based machine learning algorithms for the characterization of undersaturated oil viscosity

Background: Prediction of accurate crude oil viscosity when pressure volume temperature (PVT) experimental results are not readily available has been a major challenge to the petroleum industry. This is due to the substantial impact an inaccurate prediction will have on production planning, reservoir management, enhanced oil recovery processes and choice of design facilities such as tubing, pipeline and pump sizes. In a bid to attain improved accuracy in predictions, recent research has focused on applying various machine learning algorithms and intelligent mechanisms. In this work, an extensive comparative analysis between single-based machine learning techniques such as artificial neural network, support vector machine, decision tree and linear regression, and ensemble learning techniques such as bagging, boosting and voting was performed. The prediction performance of the models was assessed by using five evaluation measures, namely mean absolute error, relative squared error, mean squared error, root mean squared error and root mean squared log error


Background
Among all the properties of a fluid, viscosity plays a key role in influencing and controlling flow behavior either in pipelines hydraulics or in porous media such as a reservoir.The term "viscosity" is described as the internal resistance fluid encounters during flow [1].In the petroleum industry, a reliable determination of oil viscosity is essential for the analysis of pressure drop which results from fluid flow through tubing or pipes, porous media and so on, for ascertaining well productivity through proper tubing, pipelines and pump size selection, for implementation of secondary recovery processes, enhanced oil recovery processes, reservoir simulation, sound design facilities and for optimum reservoir management [2,3,4,5].Physical properties and thermodynamics of reservoir fluids such as temperature, pressure, bubble point pressure (Pb), gas/oil ratio (GOR), gas gravity and oil gravity are all dependent on oil viscosity [6] which is usually measured in the laboratory at constant temperature and different pressures considering its variability under changing operating conditions.Empirical correlations can be utilized to estimate qualities over a greater range of pressure and temperature in such instances.However, there are a lot of correlations in the literature, so it is hard to know which one to utilize in each circumstance.Standing [7] developed the first widely accepted correlation for predicting saturation pressure and crude oil formation volume factor.Four parameters influenced this relationship: the gas-to-oil ratio, temperature and the gravity of both oil and gas.Many more correlations were reported for other crude samples based on the same characteristics, with more experimental data used in general than in Standing's study [8][9][10][11][12][13][14][15][16].
Correlations are being employed in predicting viscosity especially when experimental data are not available at temperatures aside from reservoir temperature.These correlations are gotten at three different conditions, namely at Pb, above Pb and below Pb.The term "saturated oil viscosity" describes the viscosity of the oil when bubble point pressure is reached occurring at reservoir temperature (Tr).Any extended reduction in pressure after Pb is reached results in the release of gases which were initially dissolved in oil.This continues until pressure declines to the pressure of the atmosphere till no more gas initially dissolved in the oil is remaining.The viscosity of such oil is called dead oil viscosity.The undersaturated oil viscosity is the viscosity of the crude oil at a pressure above the bubble point and reservoir temperature [17].If crude oil is undersaturated at the initial reservoir pressure, the viscosity will drop somewhat as reservoir pressure falls.At saturation pressure, the viscosity will be at its lowest.As the reservoir pressure is reduced below the bubble threshold, the evolution of gas from solution increases the density and viscosity of the crude oil.
Generally, the majority of the discovered oil reservoirs in the Niger-Delta deep water environments of Nigeria are often highly "undersaturated." The majority of these reservoirs continue to exist in an undersaturated state during production because of aquifer support or water flooding projects done to sustain the reservoir pressure (Pr) to always exceed Pb.Hence, the need for a reliable predictive model or tool to predict undersaturated oil viscosity has led several authors such as Refs.[18][19][20][21][22][23][24][25][26][27][28][29] that developed empirical models currently in use in the petroleum industry.For instance, Shokir and Ibrahim [30] presented a new undersaturated crude oil viscosity model using multi-gene genetic programming (MGGP).Their results specified that the new MGGP-based model is useful in the prediction of undersaturated oil viscosity.A computer-based model, MLP-NN, was designed by Moghadam et al. [31] for the estimation of the viscosity of dead oil, saturated oil and undersaturated oil.The result outcomes indicated the precision and reliability of the method.Sinha et al. [32] developed two different approaches to estimate the dead oil viscosity form.The first approach was based on a trainable explicit expression that uses a reference viscosity measured at any temperature.Our second method utilized a hybrid machine learning method.Their approach uses a richer dataset with very limited input parameters.Beal [19] proposed the use of a graphical correlation with a mean deviation of 2.7% used to estimate undersaturated oil viscosity been a dependent variable of Pr, Pb and viscosity oil at Pb. Vazquez and Beggs [25] utilized numerous sets of PVTmeasured values for their model which had −7.541% as the mean percentage error.Petrosky and Farshad [13] utilized data from the Gulf of Mexico for their correlation which had an outcome of −0.19% for relative error and 4.22% for standard deviation.Kartoatmodjo and Schmidt [33] correlation is a modification of Beal's correlation with −4.29% as the mean error.Elsharkawy and Alikhan [23] correlation had 1.2% for mean absolute relative error and 0.022 for its deviation [2].These techniques of estimating crude oil viscosity were mostly based on applied mathematical approaches and are able to estimate the viscosity of the oil based on API gravity of oil, oil formation factor (Bo), solution GOR and Tr.In addition, the correlations had an error ranging from 25 to 40% and were mostly constrained to a particular location from which the sample data were collected; hence, their accuracy level was limited as such they could not be generalized [34].
To overcome most of these challenges in the correlation of oil viscosity, soft-computing and machine learning-based computational models were adopted in estimating the viscosity of crude oil.Linear regression, artificial neural networks, support vector machines, decision tree and so many other machine learning methods (see [35][36][37][38][39][40][41][42][43][44]) have been reported to have performed more efficiently than conventional empirical correlations.A study that involved predicting the viscosity of crude oil samples from Nigeria using a neural network with 0.99 as the coefficient of correlation presented an improved performance when compared to already developed empirical correlations.Also, studies on using neural networks to estimate the viscosity of crude oil samples in Iran and Oman have performed better than existing correlations.A report of studies on estimating the viscosity of crude oil in Canada using neural networks and support vector machine was also accurately predicted, therefore affirming the applicability of these algorithms in estimating the viscosity of crude oil is reliable [45].
Consequently, learning techniques that combine different machine learning algorithms that ally to improve predictions (stacking) or decrease variance (bagging and bias boosting) were further adopted in view that the outcome will be better than single-based machine learning classifiers and regressors.These learning techniques were referred to as ensemble learning techniques.A machine learning paradigm in which numerous learners are trained to tackle the same issue is known as ensemble learning [46].Unlike traditional machine learning algorithms that attempt to learn a single hypothesis from training data, ensemble methods attempt to create a number of hypotheses and aggregate them for usage.Ensemble methods were discovered to be more efficiently applied to datasets of points as over-fitting was greatly avoided in the course of applying these techniques [45][46][47][48].Dietterich [49] outlined three key arguments for employing an ensemble-based system in his 2000 review article: i) statistical; ii) computational and iii) representational.Succinctly, ensemble learning techniques have found application in several fields and though has started gaining the attention of researchers in the petroleum engineering field.Santos et al. [50] utilized a neural network ensemble in identifying lithofacies.Gifford and Agah [51] utilized techniques that classify multiagents for lithofacies recognition.Masoudi [52] was able to identify Sarvak productive zone formation by integrating the outputs from fuzzy and Bayesian classifiers.Davronova and Adilovab [53] did a comparative analysis of the ensemble methods for drug design.Some recent works that have deployed ensembled learning could be found in Refs.[54][55][56][57][58][59][60][61].
Researchers in the petroleum industry have also found ensembled learning as a veritable tool to make phenomenal changes in the business.Anifowose et al. [62] characterized reservoir by applying the stacked generalization ensemble in order to improve supervised machine learning algorithms capabilities in prediction.Anifowose et al. [63] recommended that ensemble techniques find relevant applications in the oil and gas industry in his review work on the application of the ensemble method.Bestagini et al. [64] classified Kansas oil-field data lithofacies by applying the random forest ensemble method.Xie et al. [65], and Tewari and Dwivedi [66] presented a study which compared the efficiency of ensemble methods applied in recognizing lithofacies.Tewari and Dwivedi [67] also did a comparison work on lithofacies identification.Bhattacharya et al. [68] utilized ANN and RF learning techniques in the prediction of the daily production of gas from the unconventional reservoir.Tewari [69,70] utilized ensemble methods to estimate the recovery factors of various reservoirs.The accuracies reported in these applications supported the argument that the application of ensemble classifiers or regressors was more precise when compared to the outcome of the discrete classifiers or regressors.Not quite much has been known concerning the comparison of ensemble regressors in predicting the PVT property of reservoir fluids such as oil viscosity.Given these observations, this study seeks to compare single-based machine learning techniques and ensemble learning techniques in predicting undersaturated oil viscosity.
The outcome of this paper brings more clarity to the applications of diverse ensemble techniques in predicting oil viscosity.Again, it helps the petroleum/reservoir engineer with an optimum predictive modeling tool for their tasks as a wide range of methods have been developed and harnessed.It also exposes additional tools and techniques to improve accuracy in the prediction of oil viscosity.Finally, to the best of authors' knowledge, this work has presented a more robust comprehensive comparative analysis between single-based machine learning techniques and ensemble learning techniques than it has been in previous literature.

Single-based machine learning algorithms
This study adopts four (4) baseline machine learning techniques, namely ANN, SVM, DT and LR, based on their superiority and widespread applications in estimating PVT properties of reservoir fluid in research.
(1) Decision tree (DT) algorithm: This algorithm is in the form of a tree-like structure that makes use of a branching methodology in clarifying each likely single result for a particular prediction/decision.Decision tree in recent times has been increasingly applied for various classification tasks because of their simplicity, ease of interpretation, low-slung cost of computation and graphical representation ability.The appropriate property for each generated tree node uses an information gain approach and any attribute that has the maximum information gain is selected as the current node for the test attribute.The operation of a decision tree algorithm on a dataset (DS) is expressed below.Firstly, the entropy value estimate of DS is presented mathematically in Eq. 1 as, where E(S) is the entropy of a DS collection, p i represents the number of proportion instances which belongs to the class i and m represent the number of classes in the system.Secondly, information gain (IG) for a particular attribute, for instance, K in a collection S is presented mathematically in Eq. 2 as, where S u is the set of value v for attribute K instances.(2) Support vector machine (SVM) algorithm: SVM is also a type of supervised machine learning technique which can be applied for both classification and regression tasks.It serves as the linear separator inserted between two data nodes to identify two different classes in environs composed of multiple dimensions.The implementation of SVM is as follows.Assume DS to be the training dataset, DS = (x i , y i , ..., (x n , y n ))ǫX, R .In support vec- tor machine, DS is represented as points in an N − dimensional space and then attempts to develop a hyperplane that will divide the space into specific class labels with a right margin of error.Support vector machine optimization algorithm is represented mathematically in Eq. 3 as, Mapping θ of vectors u i (DS) to a higher- dimensional space.The support vector machine then locates best-fit margin that divides the hyperplane linearly in this dimension.The formulation of the kernel function is presented as K (u i , u j ) ≡ θ(u i ) T (θu j ) .The radial basis function (RBF) kernel adopted in this study can expressed mathematically in Eq. 4 as, (1) where (u i − u j ) is the Euclidean distance between two data points.(3) Artificial neural networks (ANN) algorithm: ANN also referred to as neural network (NN) is a connection of interdependent layers that receives input, initializes it and passes it over to the next available component.The multilayer perceptron (MLP) for the ANN was adopted in this work.Multilayer perceptron maps a function f (.) : R D → R 0 , and dataset (dS) is trained, where (d) stands for dS input dimension, and 0 represents the number of dimensions of required output.Assume U is a set of features, and z expected output, taking U = {u 1 , u 2 , u 3 , ..., u D } , the MLP now learns the nonlinearity of the function approximators for classification and regression tasks.Stochastic gradient descent, limited-memory-Broyden Fletcher-Goldfarb-Shanno (lbfgs) or adaptive moment estimation (ADAM) can be used in training multilayer perceptron but the Tikhonov regularizer and Adam optimizer have been utilized in this work.In each layer, the activation function that was adopted in this work is rectified linear unit (ReLU) function.For each layer, the mapping functions can be expressed as Eq. 5, while the algorithm for backpropagation was utilized to train the MLP.
ω l is the weight matrix, β l is the bias and α is the weighted input sum.

Ensemble learning methods (EMs)
EM is an approach in machine learning in which combined models (usually referred to as "weak learners") that have been trained to solve a particular problem with their results aggregated to yield more improved output [86].Weak models when appropriately combined yields more robust and/or accurate models which is the primary hypothesis upon which ensemble learning is based on.Thus, ensemble methods accentuate the strength and reduce the weakness of the single classifiers or predictor.In ensemble method, diverse single classifiers or regressors are independently trained using same or differing dataset though not having the same variables [87].
The target output is derived by determining the mean of the individual single-based classifier/regressor output.Figure 1 shows an illustration of bias-variance trade-off.
To be able to solve a problem, the degree of freedom of the proposed base learning model should be sufficient to solve any underlying complexity of the dataset been worked upon in addition, the degree of freedom required (5) Z l = ω (l)T × α (l−1) + β (l)  should also avoid high variance so the model could be robust in both classification and regression tasks.This condition is referred to as bias-variance trade-off.The following factors need careful consideration when applying ensemble classifier and regresses models; difficulties in identifying the most suitable classifier or regressor for a particular application domain due to the numerous available methods, the amount of single regressors or classifiers to combined for improved accuracy as well as combination technique that will be most suitable for the different single classifiers and regressors to yield a better outcome.Therefore, presented below are discussions of few techniques for ensemble learning.Ensemble learning techniques can are grouped into diverse powerful methods such as bagging, boosting, voting, stacking and so on.
(1) Bagging approach Bagging stands for bootstrap aggregating which mostly finds its application in regression and classification tasks.Bagging increases models' accuracy by adopting decision trees techniques by reducing the variance to an optimum point.One of the major challenges of most predictive models is over-fitting but by reducing variance which invariably increases accuracy this challenge is eliminated.The advantage of bagging is that most often the learners with a weak base combine to form a more stable single learner while its high computational cost is one of its major limitations.In addition, when the appropriate procedure for bagging is ignored, it could result in more bias in models.There are two classifications of bagging, namely bootstrapping and aggregation.In bootstrapping, samples are gotten from the entire populace (set) by a procedure called the replacement method.The replacement method aids in making the selection process to be randomized.
To complete the procedure, the base learner runs on the samples.In aggregation, all possible prediction outcomes are usually incorporated and randomized.Aggregation is a function of the probability of the bootstrapping procedures or based on all the results of the models used in prediction.(2) Boosting approach Boosting learns the mistakes of the past predictor for improved and enhanced future predictions.This approach combines many base learners that are weak to obtain one strong learner, which improves the model's predictability significantly.In boosting, weak learners are arranged in sequence so that each of them learns from the preceding one in the sequence, thereby resulting in a better predictive model.The different forms of boosting include gradient boosting, adaptive boosting (also referred to as AdaBoost) and extreme gradient boosting (also referred to as XGBoost).AdaBoost uses weak learners in the form of decision trees, with one split that is mostly referred to as decision stumps.This decision stump contains several observations of like weights.In gradient boosting, predictors are added sequentially so that leading predictors correct their upcoming predictors targeted at improving the model's accuracy by canceling out the effect of the predictor's errors in the previous ones, the new predictors are usually fitted for this purpose.The gradient booster identifies the problems in the predictions of the learners through the help of gradient descent, and counters them necessary.In XGBoost, decision trees are used alongside boosted gradient to improve the speed and performance of the learner.XGBoost depends mainly on the speed of computation and the target model's performance.(3) Voting approach A voting ensemble is also called a majority voting ensemble.It is an ensemble technique which combines multiple models' predictions and is used to improve the model performance of any single model.Voting can be applied for both regression and classification tasks.Regression tasks usually require computing the mean of the model's predictions; hence, predictions are the average of contributing models.For classification tasks, each label's prediction is usually presented so that the predicted one is the most voted label.The voting ensemble may be in the form of a meta-model or a model of models.The term "meta-model" describes a set of existing trained models while the existing models would not be aware of their usage in the ensemble.A voting ensemble is better applied if two or more models with good and aggregable predictions.In a nutshell, ensemble approaches combine numerous models rather than using just one in order to increase the accuracy of outcomes in models.The combined models considerably improve the accuracy of the outcomes.This has increased the acceptance of ensemble techniques in machine learning.By integrating numerous models into one highly dependable model, ensemble approaches seek to increase predictability in models.The three most common ensemble techniques are stacking, bagging and boosting.For regression and classification, ensemble approaches work best because they lower bias and variance and increase model accuracy.

Study framework
Figure 2 shows the framework adopted in this research work which implemented single-based machine learning algorithms (ANN, SVM, DT and LR) and ensemble learning algorithms (bagging, boosting and voting) using the ANN, SVM, DT and LR for predicting undersaturated oil viscosity followed by an evaluation of their accuracies and error measurement, respectively.The Fig. 2 Study framework procedure for this study shown in Fig. 2 is classified into three main phases, namely the data preprocessing stage, the model building stage which includes both the singlebased methods and the ensemble methods, respectively, and finally the comparison of the models' accuracies and error metrics.

Research data
Experimental PVT dataset of fifteen (15) bottom-hole samples collected from different Niger-Delta oil fields in Nigeria yielding a total of one hundred and three (103) data points were used in the models.After running the training and testing scenarios on the dataset for both single-based learners and a sample of the ensemble methods, for model training, 70% of the total dataset was utilized while the remaining 30% was utilized to predict undersaturated oil viscosity values.This research involves experimental PVT data of fifteen (15) bottomhole samples taken from different Niger-Delta oil fields in Nigeria.These reservoirs were still in the "undersaturated" condition as some were been produced with the support of water flooding while others were at the early stage of initial production when the sample was collected for analysis in the laboratory.The bottom-hole samples were flash separated to derive the API gravity of oil, relative density of gas and solution GOR values.The viscosity data were obtained by rolling ball viscometer at varying pressures.Data from four reservoirs were used in this study.The total selected data points are one hundred and three (103).The range of values for the data is outlined in Table 1.

Model setup and evaluation
Scikit-learn library [88] and Python programming language was used to build a total of twenty-three (23) different models in this study of predicting undersaturated oil viscosity which comprises of four (4) using singlebased machine learning algorithms, five (5) using bagging ensemble method, six (6) using boosting ensemble method and eight (8) using voting ensemble method.The base learners' parameters were set as follows: For the ANN, multilayer perceptron was adopted with thirteen (13) hidden layers (HL).The maximum iteration was set to 2000, optimizer = lbfgs, activation = ReLU activation function.For SVM, the RBF kernel was utilized.The decision tree method setting, criterion = MSE .These were implemented on an intel corei7 64 bits with 64 GB memory laptop.
The determination of the performance between ensemble methods and single-based machine learning techniques involved using these diverse approaches to estimate the measured dataset for the undersaturated oil viscosity property.Followed by a comparative analysis between the outcome of ensemble methods and single-based machine learning methods to measure their agreement with experimental data to ascertain the efficiency of each of the approaches.The model performances were assessed under the design format of 70% training dataset and 30% testing data.To evaluate the model performance, five (5) widely known evaluation metrics were used, namely five evaluation measures, namely MAE, R 2 , MSE, RMSE and RMSLE and their formulas are shown in Eqs.6-10.z i are predicted values, u i are observed values, ẑ is the mean predicted values, û is the mean observed values, and n is the sum of the instances.A regression model's accuracy will be higher if its MAE, MSE and RMSE values are lower.However, a greater R 2 value is seen as pref- erable.The effectiveness of a linear regression model in fitting a dataset is quantified by all evaluation metrics, including MAE, MSE, RMSE, RMSLE and R 2 .While R 2 indicates how well the predictor variables can account for 2 The ability of a system or model to maintain stability and undergo only little (or no) changes when subjected to noisy or inflated inputs is known as robustness.Therefore, an outlier's impact on a robust system or measure must be reduced.Since the squaring of the errors will place higher importance on outliers in this scenario, it is simple to conclude that some evaluation measures, like MSE, may be less robust than others, like MAE.Then, using RMSE, the MSE error is square-rooted to return it to its original unit while keeping the characteristics of penalizing higher errors.As such, the five measures were deployed to take care of both the robust and non-robust systems.
The entropy value estimate of the dataset (DS) as shown in Eq. 1 essentially informs us how impure a set of data is.Non-homogeneity is described here by the term "impure." In other words, homogeneity is measured by entropy.It provides information on an arbitrary dataset's impureness and non-homogeneity.The entropy of a set of instances or dataset, DS, which includes examples of a target idea in both positive and negative light with respect to this Boolean categorization is given in Eq. 11 as, P ⊕ is the portion of positive examples and P ⊖ is the por- tion of negative examples in DS.Summarily, when the data collection is entirely non-homogeneous, the entropy is highest; when it is homogeneous, it is lowest.

Results
To ascertain an optimal number of hidden layers for multilayer perceptron/ANN to be used, a range of experiments was carried out by testing 1 to 109 hidden layers.Figure 3 shows that the ANN achieved the highest accuracy (R 2 = 0.990) at different hidden layers.However, in this study, thirteen (13) hidden layers were adopted because the error values obtained using thirteen (13) hidden neurons were considerably low (MAE = 0.0009, MSE=1.20 − e06, RMSE=0.001095)when compared to other numerous hidden layers that had the same value for Coefficient correlation with it.
Model proficiency evaluation formed the first part of the analysis.The correlation coefficient (R 2 ) is calculated for each of the models under different train/test conditions.Based on the prediction performance gotten, as well as consideration of the recent application of some specific single-based machine techniques in the petroleum industry in predictive modeling, only LR, SVM, (11) Entropy(DS) = − P ⊕ log 2 P ⊕ + P ⊖ log 2 P ⊖ ANN and DT out of five regressors were chosen to ascertain the optimal training/testing scenario.Altogether, five different training and testing scenarios based on different training and testing dataset ratios were tried and tested for each single-based machine learning algorithm and bagging ensemble method as shown in Figs. 4 and 5.As can be observed in Fig. 4, the representation of MLP in almost all scenarios has been good, with the highest correlation coefficient equal to 0.99946 under 70-30% scenario which no other single-based machine learning technique could achieve.On the other hand, the performance of SVM was the lowest with correlation coefficient equal to − 0.09743 under 90-100% scenario.The correla- tion coefficients obtained by each of the techniques were high except for SVM which was extremely low.However, the consistency of ANN to standout in all the scenarios is quite noticeable followed by LR and then DT which were at a close range to ANN. Figure 5  It is based on this result that the dataset ratio adopted throughout this work is 70% for training and 30% for testing.
After training all the 23 models built in this work using both the single-based machine learning techniques and ensemble method techniques with 70% of the dataset, the Fig. 3 Hidden layer selection using correlation coefficient models were up to be tested and evaluated.The second part of the analysis is the comparison of individual performance of the single-based machine learning methods.This was achieved by utilizing the last data group comprising 30% of the total dataset (29 data points) gathered for this work, which was not seen by any of these developed models during training was used.Table 2 presents the prediction performance of the single-based machine learning algorithm on the 30% test data used in this work.From the information provided in Table 2 and comparing the evaluation metrics of each of the methods, SVM has the worst prediction performance for the dataset used in this work, while the other three algorithms performed significantly well with ANN prediction performance topping every other single-based algorithm having the highest (R 2 ) value (0.99665) which is very close to 1.00 depicting interdependence of variables as well as the lowest RMSE value (0.00042).
For an in-depth comparative analysis, the different ensemble methods (bagging, boosting and voting) were applied using the single-based machine algorithms (ANN, SVM, DT and LR) under study as base learners to further reveal the difference in prediction performance of the single-based machine learning techniques and the ensemble learning techniques.Tables 3, 4 and 5 present the prediction performance of the various ensemble methods under study.From the information provided in Table 3, the application of the bagging approach of ensemble technique resulted in significant improvement in the R 2 value and reduction of MAE, MSE and RMSE values, respectively.However, R 2 value of the bagged-SVM slightly improved when compared to its corresponding value in Table 2 when it is a singlebased algorithm.Still, bagged-SVM performed the poorest when compared to other bagged base learners.This is obvious from the negative value of the R 2 value ( − 0.14575) which depicts that this model could decipher the interdependent relationship of the variables in the dataset.Also under the bagging procedure of ensemble learning, bagged-LR stands out with the highest value of R 2 value (0.99894) which is close to 1.00 and the lowest error evaluation metrics values (MAE = 0.00204, MSE = 5.723E-06 , RMSE = 0.00456) across board; all indicating an excellent prediction.Its comparison with the prediction performance of single-based LR also reveals significant improvement.
From Table 4, the application of the boosting approach of ensemble technique on the base learners resulted in AdaBoost-LR having the best prediction performance with R 2 value (0.99930) and very low errors evalua- tion metric values (MAE = 0.00165, MSE = 3.76101e − 05, RMSE = 0.00193).It is also observed that boosting approach performed better than the bagging approach with the LR, ANN and SVM single-based algorithm as base learners.
From Table 5, the application of the voting approach of ensemble technique on the combination base learners further enhanced improvement in prediction performance.The most conspicuous of them all is SVM, which is a single-based algorithm with very low and negative R 2 values, indicating its weakness in prediction.When combined with other single-based algorithm in voting, significantly rose with R 2 values (0.73399, 0.87540, 0.88078, 0.88015 and 0.92915) for each of the combinations in voting approach.Figures 6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26     The non-fitness of the points on the 450 line reveals a poor agreement between the measured values and predicted values.Hence, for this dataset, we can conclude that those respective models are weak and unreliable.However, this situation was the difference in the voting approach of the ensemble method as all the combinations involving SVM turned out a good fit as can be found in Figs. 20, 23, 24, 25 and 27.The fitting of the data points on the 450 lines reveals the accuracy of the predicted results using these approaches.This finding indicates that ensemble methods with voting in this particular case have the potential of boosting weak single-based machine learning algorithms.This is evident the way the previous weak single-based SVM could produce an accurate prediction of the undersaturated oil viscosity having combined with other techniques using the voting ensemble method.

The predictive models
The error evaluations using the five measuring metrics reveal that the ensemble learners had the lowest MAE, MSE, RMSE and RMSLE values and highest R 2 val- ues when compared with the single-based algorithms.Tables 4 and 5 show that the approach described in this work is suitable for predicting undersaturated oil viscosity, with the bagged ensemble method providing the best solution.The Al-Kafaji [89] and Bennison [90] procedures are ineffective for crude oil gravities less than a particular threshold.The most frequent approach for obtaining oil viscosity data is viscosity correlation, which is extremely helpful and successful in estimating oil viscosity at various temperatures for various oil kinds.The limitation on the parameters from which these correlations have been formed is the essential element in their implementation.As a result, these connections are region-specific and cannot be generalized worldwide.
According to the literature, measuring heavy oil viscosity at low temperatures might be difficult since the predicted values of the viscosity frequently surpass the maximum limit of the equipment.This model is suitable for a wide range of oil viscosities and temperatures.
The suggested oil viscosity correlation does not need any compositional investigation of the oil or asphalt content, which is a significant benefit of the ensembled model.The created model outperforms both a singlelearner method and the leading correlation in terms of predictive potential.The model may be used as a quick tool to validate the quality of experimental data and/or the validity and accuracy of various viscosity models, particularly when there are differences and uncertainties in the datasets.The suggested model, which has increased the accuracy and efficiency of crude oil viscosity, may be used in any reservoir simulator software.This demonstrates the ensemble leaners algorithm's

The sensitivity analysis of the ensemble learners
To show the impact of all input factors, such as GOR, temperature, gas gravity and oil gravity, the effects of the individual independent variables on the ensemble learners were examined.The outcomes of the ensemble learners' sensitivity analysis are shown in Fig. 28.The rank correlation coefficients between the output variable and the samples for each of the input parameters are displayed in this figure.In general, the effect of any input in deciding the value of the output grows as the correlation coefficient between that input and the output variable rises.The chart makes it clear that the main factor affecting oil viscosity is GOR.

Conclusion
A rigorous comparison is presented between two different types of regression schemes viz: single-based machine learning techniques (such as ANN, SVM, DT and LR) and ensemble learning techniques (such as bagging, boosting and voting) to predict undersaturated oil viscosity of Nigeria's crude oil.Four (4) variables were utilized, namely API gravity, gas/oil ratio, bubble-point pressure and reservoir temperature.The sensitivity analysis reveals that the oil viscosity is mostly affected by the gas/oil ratio, GOR.
The response variables in absolute terms were measured via the MAE, MSE, RMSE and RMSLE.When compared to the single-based learner, the values of these error evaluation metrics are much lower in all ensemble approaches.In the same way, the R 2 values are higher in ensemble methods than in the singlebased learner counterpart.This shows that ensemble models are consistent and reliable when compared with their respective individual single-based models.However, the results obtained during the testing and training using the bagging ensemble learners reveal that the developed bagged ensemble method has better capabilities in making known the uncertainties embedded in the data as their correlation coefficient values were significantly improved when compared with the other learners.It also reveals that the created model outperforms both a single-learner method and the leading correlation procedure in terms of predictive potential.
Through the voting ensemble technique, the weak prediction made by single-based SVM had a Finally, this study looked at a method for estimating viscosity regardless of temperature, oil type or other crucial characteristics that are difficult to collect experimentally.In the quest for more rigorous and reliable tools by petroleum and reservoir engineers, future works in this domain should incorporate ensemble learning techniques in predicting other PVT properties such as GOR, Bo and isothermal compressibility of oil.

Fig. 1
Fig. 1 Illustration of bias-variance trade-off reveals that the correlation coefficient values of bagged-LR, bagged-DT and bagged-ANN are considerably close ties with just a minute difference across all the five different training and testing scenarios based on different training and testing dataset ratios.It is also worthy of noting that the correlation coefficient values of DT, LR and ANN of the single-based algorithm were improved upon in the bagging ensemble method with the most conspicuous being the DT.Yet 70-30% of training and testing scenarios gave the highest correlation coefficient equal to 0.99964.

Fig. 4 Fig. 5
Fig. 4 Training and testing scenario for single-based machine learning algorithms and 27 illustrate cross plots of the predicted undersaturated oil viscosity values versus the experimental/measured undersaturated oil viscosity values.The cross plot reveals the level of coherence between the measured and the estimated values.A perfect agreement should have all the points lie on the 450 line on the cross plot.In every instance, all the graphs showed the cloud of tightest points around the 450 lines.This signifies an accepted agreement apart from Figs. 7, 11 and 17 representing the cross plots for

Fig. 10
Fig. 10 Cross plot of bagged-ANN predicted values versus measured values

Table 1
Statistical description of the PVT data the variation in the response variable, MAE, MSE, RMSE and RMSLE indicate how well a regression model can predict the value of a response variable in absolute terms.

Table 2
Prediction performance of different single-based machine learning algorithms

Table 3
Prediction performance of bagging approach type of ensemble technique

Table 4
Prediction performance of boosting techniques of ensemble method

Table 5
Prediction performance of voting approach type of ensemble technique SVM, bagged-SVM and AdaBoost-SVM, respectively.