# On the Estimation of the Density of Brine with an Extensive Range of Different Salts Compositions and Concentrations

^{*}

**Corresponding Author(s):**Halali MA, Young Researchers and Elite Club, North Tehran Branch, Islamic Azad University, Tehran, Iran, Email: [email protected] Mohammadi AH, Institut de Recherche en Génie Chimique et Pétrolier (IRGCP), Paris, France, Email: [email protected]

*
Received Date: May 01, 2016 /
Accepted Date: May 20, 2016 /
Published Date: May 26, 2016 *

### Abstract

Accurate determination of brine density is crucial for an efficient designing of various systems through where brine take an undeniable part. The importance of determining brine density is well corroborated by recent ongoing research surge in this scope. Despite existence of several experimental and empirical approaches, a reliable and robust model seems to be requisite for precise determination of brine density. In regard to high performance and great robustness of soft computing approaches for unraveling science and engineering problems, this article proposes LS-SVM and MLP-NN models to determine the brine density. Both models are developed over 1868 data points including both natural and synthetic brines. The proposed models take account an extensive range of input parameters such as temperature, pressure, and concentration. The developed models can significantly estimate the target values with respect to high values of R2, which are 0.999999 and 1.000000 for MLP-NN and LS-SVM models, respectively. Considering high accuracy and swiftness of proposed models, they can be great assets to science and engineering scopes.

**Keywords:** LS-SVM; MLP-NN; Model; Brine density; Natural and synthetic brine.

#### Nomenclature

AARD: Average Absolute Relative Deviation; ANFIS: Adaptive Neuro-Fuzzy Inference System; ARD: Average Relative Deviation; ANN: Artificial Neural Network; LS-SVM: Least Square Support Vector Machine; MLP-NN: Multilayer Perceptron Neural Network; MSE: Mean Squared Error; RMSE: Root Mean Square Error; STD: Standard Deviation Errors; N: the Number of data points

#### Introduction

Nowadays, brines properties play a pivotal role in various engineering and scientific scopes including chemical, physical, geochemical, and geothermal domains. More specifically, these researches are mostly accomplished around fluid inclusion surveys, CO_{2} sequestration, water desalination, studies of fluid-rock interaction, fluid flow simulation and enhanced oil recovery (EOR) [1-3]. Take for instance drilling industry, which is well recognized as a costly operation, most of drilling fluids categorized as water based fluids are whether built on the basis of brine or in evident contact with salt formations, either ways the brine density take an undeniable part during the whole drilling procedure. Conventionally, a drilling operation faces serious challenges, which are in close relation with brine density and solubility through the fluid [4]. In other words, the brine density can be a representative index of well overhead pressure as well as affecting numerous parameters such as pH, contaminations (Ca^{++}, Mg^{++}), shale swelling problem, wash out problem, waste management, etc. Therefore, by determining the brine density precisely not only can we pave the way for a safe and secure operation but also we can manage different problems appropriately. Moreover, brine completion fluid has widespread application in completion and workover operations [5]. Generally, such fluids own high density to provide high fluid column hydrostatic pressure as well as overcoming high formation pore pressure. The aforesaid fluid should be solid free, thus the defined fluid weight should be obtained from saturated salts solutions including CaCl_{2}, CaBr_{2} and ZnBr_{2} [6]. The effect of temperature, pressure, and concentration on the brine properties is undeniable. Thus, precise estimation of brine density is pressure, which cause overbalance condition but also to avoid high cost of expensive completion fluid due to its overproduction in case the well experiences serious loss through the formation. In addition, the phase equilibrium and PVT properties of the CO_{2}-NaCl–H_{2}O system have been investigated broadly since it has extensive applications in geochemical processes [7,8]. A wide research effort in determining PVT properties of CO_{2}-NaCl–H_{2}O geological fluids in industrial processes including the oil and gas production, Enhanced Oil Recovery (EOR), natural gas clathrate engineering, geothermal exploitation, the exhaust gases treatment, supercritical fluid extraction and oxidation, wastewater and waste liquids, fertilizers, and seawater desalination, etc. vouches for the fact that recent interests are developing around PVT behavior of CO_{2}-NaCl–H_{2}O geological fluids [9]. Interestingly, obtaining precise PVT properties of brine is requisite for avoiding errors in prediction of other petroleum parameters [10]. As a matter of the fact, with a recent considerable growth in oil production rate, the amount of brine production increases accordingly. There is no doubt that managing such a great amount of produced brine requires a tremendous cost. More elaborately, aquifers are the available water surround hydrocarbon reservoirs and assist the reservoir to produce robustly in various ways: peripheral water drive, edge water drive, and bottom water drive [11]. Through so-called water drive approach the new wave of brine pushes the hydrocarbon and increases the reservoir pressure in order to produce with a higher flow rate. Unfortunately, brine may also be produced in conjunction with hydrocarbon. The required not only to prevent inaccurate determination of overhead pressure, which cause overbalance condition but also to avoid high cost of expensive completion fluid due to its overproduction in case the well experiences serious loss through the formation. In addition, the phase equilibrium and PVT properties of the CO_{2}-NaCl–H_{2}O system have been investigated broadly since it has extensive applications in geochemical processes [7,8]. A wide research effort in determining PVT properties of CO_{2}-NaCl–H_{2}O geological fluids in industrial processes including the oil and gas production, Enhanced Oil Recovery (EOR), natural gas clathrate engineering, geothermal exploitation, the exhaust gases treatment, supercritical fluid extraction and oxidation, wastewater and waste liquids, fertilizers, and seawater desalination, etc. vouches for the fact that recent interests are developing around PVT behavior of CO_{2}-NaCl–H_{2}O geological fluids [9]. Interestingly, obtaining precise PVT properties of brine is requisite for avoiding errors in prediction of other petroleum parameters [10]. As a matter of the fact, with a recent considerable growth in oil production rate, the amount of brine production increases accordingly. There is no doubt that managing such a great amount of produced brine requires a tremendous cost. More elaborately, aquifers are the available water surround hydrocarbon reservoirs and assist the reservoir to produce robustly in various ways: peripheral water drive, edge water drive, and bottom water drive [11]. Through so-called water drive approach the new wave of brine pushes the hydrocarbon and increases the reservoir pressure in order to produce with a higher flow rate. Unfortunately, brine may also be produced in conjunction with hydrocarbon. The more the reservoir produces, the lower the reservoir pressure becomes and the more the brine will be produced. The production of brine along with hydrocarbon would have an adverse effect on the produced hydrocarbon both qualitatively and quantitatively [12]. This would even cause well shut-in due to lack of appropriate surface treatment. It is important to note that in some cases even in reservoirs enjoying most recent modern approaches the produced fluid contains more than 90% brine in volume [13]. Furthermore, it is contented by many recent research surveys that the hydrocarbon recovery through the water flooding approach can be boosted by means of adjusting the salinity of injected brine [14]. To this end, brine properties should be determined with the utmost accuracy. The most outstanding properties are density, viscosity, and solubility. Getting acquainted with aforementioned properties is crucial for evaluating the effect of volume and movement of water in the reservoir, wellbore, and through the surface facilities [15]. In general, widespread applications of brine properties in other fields induce a vast surge in literature to be engendered elaborately [16-19]. The properties of brine can be obtained through different approaches including laboratory experiments, available models and correlations, and also soft computing methods. Actually, laboratory studies were recognized to be the most solid and precise method. Nevertheless, this approach was found to be expensive and timeconsuming [20]. Thus, in the absence of laboratory experiments, other methods such as implementing empirical models and correlations were preferred to determine brine properties. However, these models are suffered from uncertainties and limitations [21]. All in all, apart from the time and cost issues, the success of different vital operations pivots on implementing the right approaches. Hereupon, the power of economy in associate with high degree of accuracy expected for determination of brine properties force us to seek for new alternatives to fulfill this demand. In following a review of all mentioned methods is provided.

**Determination of brine properties**

Basically, there are three main approaches appertain to determination of brine properties.

1- Methods, which are on the basis of laboratories surveys;

2- Methods, which link the brine properties to the empirical correlations; and

3- Methods, which are based on soft computing analysis.

These three approaches are discussed further below.

**Laboratorial studies**

Thus far, extensive experimental studies have been conducted to determine brine density in a wide span of temperature, pressure, and concentration. Gibson and Loeffler [22] reported the specific volume of brine for temperatures between 29815.15 and 385.15 K, pressures from 0.1 to 100 MPa, and salinity from 5 to 25 wt% NaCl. Hilbert [23] accomplished a massive research work to measure the specific volume of brine for temperatures, pressures, and salinities from 293.15 to 673.15 K, 10 to 400 MPa, and from 0 to 25 wt% NaCl, respectively. Ghafri et al. [24] reported the density of NaCl (aq) in a temperature range between 283 and 472 K and pressures up to 68.5 MPa and molality of 1.06, 3.16, and 6 mol/kg. In a similar attempt, Kumar [25] measured the density of SrCl_{2} (aq) for temperatures of 323.15 to 473.15 K and at pressure of 202.7 MPa and for concentrations up to 2.7 mol/kg. Pitzer et al. [26] conducted an experimental study to measure the densities of MgSO_{4} and Na_{2}SO_{4} up to 470 K, 10 MPa and 1 m ionic strength. Similarly, Obsil et al. [27] examined densities of Na_{2}SO_{4} (aq) and K_{2}SO_{4} (aq) using vibrating-tube densitometry at temperatures from 298.15 K to 572.7 K and pressures up to 30.7 MPa as well as at concentrations of 0.01 mol/ kg to 1.0 and 0.5 mol/kg, respectively. Moreover, the density of NaCl (aq) and other salts were reported by Gates and Wood [28] from 0.1013 to 40 MPa and temperature of 298.15 K and in a concentration range of 0.05 to 5.0 mol/kg. Crovetto et al. [29] measured the vapor pressure and density of NaCl (aq) at temperature of 623 K and at molalities of 0.25, 0.5, 1, and 3 mol/kg. Sharygin et al. [30] investigated the densities of aqueous solutions of Na_{2}CO_{3} and NaHCO_{3} for temperatures in a range of 298 K to 623 K and pressures up to 28 MPa and with molalities from 0.1 mol/kg to 1.0 mol/kg implementing vibrating-tube flow densimeter. Further experimental research efforts can be found elsewhere in open literature [31-35].

**Empirical correlations**

Hass [36] presented a model to estimate the density of vaporsaturated NaCl (aq) applying the empirical Masson’s Rule. This model is applicable to predict density in a temperature range of 348.15 to 598.15 K and up to saturation of 7.3 molal. Similarly, Philips et al. [36] brought forward a model to determine brine density for temperatures, pressures, and malalities of 283.15 to 623.15 K, up to 50 MPa, and 0.25 to 5 mol/kg, respectively. This model can predict the experimental values with maximum uncertainty of ± 2%. In another attempt, Rogers and Pitzer [37] conducted a very elaborated computational research to introduce a semiempirical model in order to explain thermal properties of NaCl (aq) for concentration range of 0.1 to 5 molal as well as determine volumetric data for temperatures from 273.15 to 573.15 K and pressures of 0 to 100 MPa. As a result, they provided comprehensive reports for values of expansivity, specific volume, and compressibility versus temperature, pressure, and molality. Ghafri et al. [1] developed an empirical model to predict the density, apparent molar volume, and isothermal compressibility of various brines employing corresponding experimental data. They were taking into account the whole range of experimental data including temperature, pressure, and malality. Bahadori et al. [38] presented an Arrhenius type function to discern the different characteristics of reservoir brine such as density, vapor pressure, and enthalpy over a wide range of 5-25% salt content by mass, and temperatures above 303.15 K. The aforementioned model enjoys lower roundoff error by excluding massive complexities of mathematics and provided in a pretty user-friendly configuration. More correlations are available in literature [39,40].

**Soft computing techniques**

Nowadays, soft computing approaches are well recognized as beneficial robust tools, which take a significant part in analyzing and unraveling challenging problems in various scopes of science and engineering (for instance [41-46]). The preponderance of these computer based approaches includes determining a target function with high degree of precision compared to many published empirical and analytical solutions [47,48].

Intelligent methods including the artificial neural network (ANN) [49-51], Adaptive Neuro-Fuzzy Inference System (ANFIS) models [52,53], Least Square Support Vector Machine (LS-SVM) [54-56], Multilayer Perceptron Neural Network (MLP-NN) [45,57-59], and Radial Basis Function Neural Networks (RBF-NN) [60-62] are amazingly robust and reliable tools for data analysis and interpretation that can be employed to predict regression and classification problems.

LS-SVM and MLP-NN as intelligent methods have been satisfactorily implemented in many scopes of science and engineering to unravel very complexes computational challenges; a fact evidently acknowledged by recent researches organized around investigating PVT properties of reservoir hydrocarbons [54,55,63,64], pressure gradient estimation of multiphase flow [65,66], and retention and solubility determinations [67,68]. In addition to previous domains, Arabloo et al. [69] proposed a LS-SVM algorithm to estimate properties of reservoir brine including liquid saturation vapor pressure, density and enthalpy. The model results were in acceptable agreement with the experimental data regarding the R^{2} value of 0.999. Moreover, Kamari et al. [70] presented an analogous model to determine crude oil salt content, their model was tested over 63 data points from literature, which exposed satisfactory accuracy thanks to high R^{2} value of 0.9999. Tatar et al. [71] investigated the reservoir brine properties using RBF-NN. The results for density, enthalpy and vapor pressure were comparatively convincing with respect to low corresponding RMSE of 0.270810, 0.455726, and 1.264687, respectively. In a similar attempt, Tatar et al. [72] implemented MLP and GA-RBF models to determine reservoir formation water density. Their models exposed high accuracy since the estimations were in good agreement with experimental data regarding low RMSE of 3.27E-05.

Based on the convincing background of incorporating least square support vector machine and multilayer perceptron neural network to resolve an extensive range of engineering problems, we are pursuing ongoing development and application of LS-SVM and MLP-NN for modeling of brine density of different salts. Therefore, the objective of this study is to develop robust soft computing based models for accurate determination of the density of the brines of different salts for a wide range of concentrations. Underlying the contribution of the article, our developed LS-SVM and MLP-NN models cover a broad range of input flow data including 1868 data sets. Worth noting that the input data cover both natural and synthetic brine densities, which are gathered for the first time.

Indeed, the superior high accuracy in determining density of the brines of different salts is provided implementing the proposed LSSVM and MLP-NN models, which reveals that the aforesaid models would be great assets for science and engineering activities in this domain. The remainder of this article is organized as following; section 3 appertains to background and development of each model. Afterwards, the result of proposed models will be presented in section 4 and both models will be compared with each other. At last, section 5 brings forward a conclusion of the article.

#### Details of Model

**Data acquisition**

In order to develop a reliable model it is essential to incorporate valid and authentic data, which covers a wide range of variables [45,73-75]. A set of 1868 data point are utilized in this study, which is gathered from literature [1,24]. The density of pure brines of CaCl_{2}, MgCl_{2}, KI, NaCl, KCl, AlCl_{3}, SrCl_{2}, Na_{2}SO_{4}, and NaHCO_{3} as well as their different combinations with different salt concentration are investigated in this study. Salt concentration of the samples spans in the range of 0.359-6 mol/kg. The details of the input and output parameters are listed in **Table 1**.

Input parameter | Minimum | Maximum | Average | Standard Deviation |
---|---|---|---|---|

CaCl_{2} (mole fraction of salt) |
0 | 1 | 0.113892 | 0.305974 |

MgCl_{2}(mole fraction of salt) |
0 | 1 | 0.111973 | 0.311183 |

KI (mole fraction of salt) | 0 | 1 | 0.105996 | 0.307915 |

NaCl (mole fraction of salt) | 0 | 1 | 0.274283 | 0.413956 |

KCl (mole fraction of salt) | 0 | 1 | 0.121051 | 0.298603 |

AlCl_{3} (mole fraction of salt) |
0 | 1 | 0.034261 | 0.181948 |

SrCl_{2} (mole fraction of salt) |
0 | 1 | 0.101313 | 0.3016 |

Na_{2}SO_{4} (mole fraction of salt) |
0 | 1 | 0.069138 | 0.250543 |

NaHCO_{3} (mole fraction of salt) |
0 | 1 | 0.068126 | 0.25071 |

Salt Concentration (mol/kg) | 0.359 | 6 | 2.232862 | 1.63979 |

T (k) | 283.15 | 473.02 | 379.5151 | 57.82789 |

p (Mpa) | 0.9 | 68.6 | 35.40075 | 22.17534 |

ρ (kg/m^{3}) |
891.65 | 1406.89 | 1098.936 | 100.0464 |

**Table 1:** Statistical values of the input and target data.

**Backgrounds of modeling**

Artificial neural networks are algorithms capable of learning from experience, improving their performance and adapting themselves to the changes in the environment [76]. ANNs have been used for monitoring, controlling, classification, and simulation of activated sludge processes. ANNs can provide numerous benefits such as ability to process a large amount of data and generalize the results, require much less statistical training process, able to diagnose nonlinear relationships between dependent and independent parameters implicitly as well as capable of diagnosing all possible interactions between predictor parameters.

On the other hand, the “black box” nature of ANNs, the great computational stress, susceptibility to over-fitting, and the empirical nature of model generation are the most evident disadvantages of ANNs [44,77]. The two implemented models are described further below.

The support vector machines: The support vector machine (SVM) is a supervised learning technique from the field of machine learning with capability of both classification and regression analysis [78-81]. On the other side, the necessity to resolve a large-scale quadratic programming problem is one of the main defects of the SVM [82]. To deal with this challenge, a new technique so-called Least-Squares SVM (LS-SVM) is presented, which is a modification of the traditional SVM. This technique solves linear equations (linear programming) rather than quadratic programming problems to attenuate the complexity of optimization process [83-85]. Considering the problem of approximating a given dataset with a nonlinear function:

(1)

Where represent dot product; Φ(x) represents the nonlinear function that applies linear regression; b and w are bias terms and weight vector, respectively. In LS-SVM for function prediction, the optimization problem is formulated as [82,86]:

(2)

(3)

where, e_{k}∈R are error variables; and γ ≥0 is a regularization constant. To solve this optimization problem, Lagrange function is established as [82,86]:

(4)

where, *α _{k}R *are Lagrange multipliers. The solution of Eq. (4) can be determined by partially differentiating with respect to

*w, b, e*and α

_{k}*[82,86]:*

_{k}(5)

By defining *1 _{v} = [1;…1], Y = [y1;…;yN], α = [α_{1};…;α_{N}]* and excluding w and e, the following linear equations are obtained [82]:

(6)

where, *I _{N}* refers to an N×N identity matrix and Ω is the kernel matrix that is introduced as [82]:

(7)

There are several kernel functions including linear, spline, polynomial, and radial basis function [87,88]. On the side, radial basis function, a Gaussian function, and polynomial function are the most widespread functions.

(8)

(9)

where, σ denotes the width of the RBF, which controls the regression capability and *d* is the polynomial degree.

**Multilayer perceptron networks: **Multilayer Perceptron Neural Networks (MLP-NNs) consist of three different kind of layers, namely, input layer, hidden layer(s), and output layer (**Figure 1**). A single MLP may have one or more hidden layer. Each layer comprises some neurons accordingly. The number of neurons in the input and output layers are attributed to the number of input and output data, respectively. The number of hidden layers as well as included neurons in them is optional and can be specified either intelligently or by trial and error to gain the most efficient performance. Minimum Square Error (MSE) demonstrates the performance of the presented network. In such networks, the error is back-propagated through the network and the weights and biases are optimized through some iteration called epochs. The number of epochs should be determined in a manner that the network neither undertrain nor overtrain. In the former, the network does not have sufficient time to complete the learning process. In the latter, the network does not learn but memorizes. This leads to poor performance of network in estimation of test data set. **Figure 1** demonstrates a general form of a fully connected multilayer perceptron with two hidden layers.

**Figure 1: **Structure of a fully connected multilayer perceptron with two hidden layers [90].

There are two different types of signals in MLP network, namely, function signal and error signal. In each node that function signal passes, it is calculated as a function of the inputs values and their corresponding weights and biases. This signal is also called input signal. No computational task is performed in input neurons. In following layers, the activation levels are calculated in each neuron and the results are propagated through the next layers until reaching the output layer. The output of the network is specified in the output layer's neuron(s). Afterwards, this value is compared with the target value and error is calculated. Because the essence of the supervised learning paradigm, error is back-propagated through the network to regulate weights and biases. Here, the other type of signal, namely, error signal is implemented. This signal propagates backward from the output neurons into the network. **Figure 2** indicates a portion of the MLP network demonstrating function and error signals [89].

Hidden layer neurons take a significant part in performance of the MLP as they behave as feature detectors. During training process, hidden neurons start to reveal the features of dedicated data for training. This is accomplished by nonlinear transformation of the input data to feature space. In the new space, it is may be easier to separate data from each other for classification. Formation of feature space makes the difference between multilayer perceptron and Rosenblatt perceptron. Each neuron uses an activation function. Sigmoid functions, which are "S" shaped graph, are the most frequent type of activation functions. Further details about Sigmoid functions is available in literature [90].

**Designing LS-SVM and MLP models representing density of brines**

The first step to start the simulation is data normalization. All the data points including input and output parameters were normalized between -1 and 1 using the following formula:

(10)

At the next step, 80% of the total data were randomly allocated for training and the rest for testing the developed networks. This division is such that there is no local accumulation of train or test data points. Matlab^{®} 2014a was utilized to implement the MLP-NN and LS-SVM codes.

**Development of LS-SVM:** In order to prevent overfitting problems, the data set was separated into two subdata sets. To this end, 80% and 20% of the main data set was randomly selected for the training set and test set, respectively.

During the computation, the extensively implemented kernel function, i.e., radial basis function (RBF), has been utilized. It has the general form as following [63,64,91]:

(11)

in which, σ is a decision parameter that is determined during the optimization calculations [84].

The Mean Square Error (MSE) between the developed model results and corresponding experimental data reported in the literature, as defined by equation (12), was considered as objective function during model computation [92,93].

(12)

in which t and o are target value and estimated value, respectively; and *n* is number of the data points.

**Development of MLP-NN:** Cybenko [94] stated that it is mathematically proved that every function can be estimated adequately by an MLP-NN with only one hidden layer. Therefore, only MLP-NNs with only one hidden layer are investigated in this study. Different MLP-NNs structures with 4-25 neurons in the hidden layer were examined and finally it was proved that the MLP-NN with 24 neurons in the hidden layer leads to the best prediction of brine density. The MSE of the examined neural networks is shown in **Figure 3**.

**Accuracy of the proposed model and validation**

Both graphical and statistical methods are utilized to validate and show the accuracy of the proposed models. Four different statistical parameters including correlation factor (R^{2}), Average Absolute Relative Deviation (AARD), Standard Deviation (STD), and Root Mean Squared Error (RMSE) are employed (Equations (13)-(16)) to investigate the accuracy of the proposed models. The formulation of these parameters is as follows:

. (13)

(14)

(15)

#### Results and Discussion

The LS-SVM has two tuning parameters called, γ and σ^{2}. Utilizing simulated annealing optimization algorithm, the optimal values for the mentioned values were determined to be 3630241399 and 7.58, respectively.

The statistical parameters obtained from the developed models are demonstrated in **Table 2**. The statistical parameters are including R^{2}, AARD, STD, and RMSE. As it can be seen, the statistical parameters are calculated for both train and test data sets. With regard to high values of correlation factor for MLP-NN and LS-SVM models, which are 0.999999 and 1.000000, respectively, both developed MLPNN and LS-SVM models enjoy high level of accuracy in predicting experimental data. Furthermore, low values of AARD, STD and RMSE for proposed MLP-NN model which are 0.006893, 0.000115 and 0.117512, respectively, validate the great authenticity of the MLP-NN model. Moreover, the low values of AARD, STD and RMSE obtained from the proposed LS-SVM models which are 0.003657, 5.11E-05 and 0.054512, respectively well demonstrate that the model estimations are in great accordance with experimental data points.

R^2 | AARD | STD | RMSE | N | ||
---|---|---|---|---|---|---|

MLP-NN | Train Data | 0.999999 | 0.006483 | 0.000101 | 0.104272 | 1494 |

Test Data | 0.999998 | 0.008532 | 0.000159 | 0.159809 | 374 | |

All Data | 0.999999 | 0.006893 | 0.000115 | 0.117512 | 1868 | |

LS-SVM | Train Data | 1.000000 | 0.003246 | 4.46E-05 | 0.047559 | 1494 |

Test Data | 0.999999 | 0.0053 | 7.15E-05 | 0.076201 | 374 | |

All Data | 1.000000 | 0.003657 | 5.11E-05 | 0.054512 | 1868 |

**Table 2: **Statistical parameters of the proposed models.

Accuracy of the developed models is further investigated by means of graphical analysis methods such as cross plot, relative deviation, and error distribution. Cross plots indicate the degree of validation, whereas error distribution explains if the developed models have an error trend.

**Figures 4-7** are indicative of comparison between MLP-NN and LS-SVM models estimated/represented brine density and target (experimental) values for both train and test data sets.

**Figures 4(a)** and **4(b)** display a scatter plot of experimental brine density versus developed MLP-NN and LS-SVM models estimation, respectively. These plots are for both train and test data sets. At it is displayed, presence of almost all data points on the bisector of the first quadrant leads us to a conclusion that a close agreement is obtained between the predictions of both MLP-NN and LS-SVM models and the experimental data points for both train and test data sets.

The degree of correspondence between experimental data points and MLP-NN and LS-SVM models estimated values is further demonstrated in **Figures 5(a)** and **5(b)**, respectively. **Figure 5** displays the relative deviations of the proposed models predictions versus experimental data points for both train and test data sets. According to the **Figure 5**, the predicted brine densities are in a desired agreement with target values since the distribution of relative deviation is close to the horizontal line and almost coincides that. For the sake of better visual comparison, the estimated and target values are displayed simultaneously versus the data point number (i.e., data index) for both MLP-NN and LS-SVM models in **Figures 6(a)** and **6(b)**, respectively. It goes without saying that coincidence of the experimental and predicted values signifies the accuracy of both MLP-NN and LS-SVM models. In order to shed more light on the reliability of developed models, the error distribution plot is demonstrated for both proposed models in **Figure 7**. For both models, as it is displayed, the error distribution is completely symmetric around the center of 0. To round it off, RMSE bar plot of both MLP-NN and LS-SVM models is illustrated for train, test, and all data sets in **Figure 8**. The result explains that the LS-SVM model slightly outperforms MLP-NN model, however, both models enjoy great accuracy with respect to low values of RMSE (0.054512 and 0.117512, respectively) and other statistical parameters including R^{2}, AARD and STD.

#### Conclusion

Various scopes of industry face a serious challenge for precise determining of brine density. These scopes are more around fluid inclusion surveys, CO2 sequestration, studies of fluid-rock interaction, drilling industry and Enhanced Oil Recovery (EOR). In spite of recent breakthrough in determining brine properties including brine density, a swift reliable model is required for an accurate determination of brine density. Hereupon, two soft computing models named, LS-SVM and MLP-NN have been proposed in this paper to address this need. These models cover both natural and synthetic brines as well as an extensive range of input parameters including temperature, pressure, and concentration.

The reliability of designed models was verified in favor of employing statistical parameters. Indeed, the LS-SVM and MLP-NN revealed R^{2} of 1.000000 and 0.999999, AARD of 0.003657 and 0.006893, STD of 5.11E-05 and 0.000115 and RMSE of 0.054512 and 0.117512, respectively. Hereupon, it can be concluded that developed models are reliable enough to predict brine density, which correspondingly they can be implemented to assist delicate industrial designations.

#### References

- Al Ghafri SZ, Maitland GC, Trusler JM (2013) Densities of SrCl2 (aq), Na2SO4 (aq), NaHCO3 (aq), and Two Synthetic Reservoir Brines at Temperatures between (298 and 473) K, Pressures up to 68.5 MPa, and Molalities up to 3 mol• kg–1. J ChemEng Data 58: 402-412.
- Huizenga JM, Jens Gutzmer, David Banks, Lynnette Greyling (2006) The Paleoproterozoic carbonate-hosted Pering Zn–Pb deposit, South Africa. II: fluid inclusion, fluid chemistry and stable isotope constraints. MineraliumDeposita 40: 686-706.
- Kumar A, Myeong HN, Robin CO, Gary AP, Larry WL, et al., (2005) Reservoir Simulation of CO 2 Storage in Aquifers. Spe Journal 10: 336-348.
- Borchardt JK, Rao SP (1985) High density well drilling, completion and workover brines, fluid loss reducing additives therefor and methods of use, Google Patents.
- Zhang H, et al., (2014) High density brines for use in wellbore fluids, Google Patents.
- Krook G, Boyce T (1984)Downhole density of heavy brines. in SPE Formation Damage Control Symposium, Society of Petroleum Engineers.
- Driesner T, Heinrich CA (2007)The system H2O–NaCl. Part I: Correlation formulae for phase relations in temperature–pressure–composition space from 0 to 1000°C, 0 to 5000bar, and 0 to 1 X NaCl. GeochimicaetCosmochimicaActa 71: 4880-4901.
- Driesner T (2007)The system H2O–NaCl. Part II: Correlations for molar volume, enthalpy, and isobaric heat capacity from 0 to 1000°C, 1 to 5000bar, and 0 to 1 X NaCl. GeochimicaetCosmochimicaActa 71: 4902-4919.
- Duan Z,Hu J, Li D, MaoS (2008) Densities of the CO2–H2O and CO2–H2O–NaCl systems up to 647 K and 100 MPa. Energy & Fuels 22:1666-1674.
- Numbere D, Brigham WE, Standing M (1977) Correlations for Physical Properties of petroleum reservoir Brines, Stanford Univ, CA (USA). Petroleum Research Inst.
- Ahmed T (2006) Reservoir engineering handbook, Gulf Professional Publishing.
- Ogunberu A,Ayub M (2007)An analytical model for brine-induced oil recovery. Petroleum Science and Technology 25: 1571-1581.
- Bahadori A,Vuthaluru HB (2010) Prediction of silica carry-over and solubility in steam of boilers using simple correlation. Applied Thermal Engineering 30: 250-253.
- Austad T, RezaeiDoust A, Puntervold T (2010) Chemical mechanism of low salinity water flooding in sandstone reservoirs. in SPE improved oil recovery symposium. Society of Petroleum Engineers.
- Spivey J, McCain Jr W, North R (2004) Estimating density, formation volume factor, compressibility, methane solubility, and viscosity for oilfield brines at temperatures from 0 to 275°C, pressures to 200 MPa, and salinities to 5.7 mole/kg. Journal of Canadian Petroleum Technology 43.
- Garcia JE (2001) Density of aqueous solutions of CO2. Lawrence Berkeley National Laboratory.
- Yan W, Huang S, Stenby EH (2011) Measurement and modeling of CO2 solubility in NaCl brine and CO2–saturated NaCl brine density. International Journal of Greenhouse Gas Control 5:1460-1477.
- Nomeli MA, Tilton N, Riaz A (2014) A new model for the density of saturated solutions of CO2–H2O–NaCl in saline aquifers. International Journal of Greenhouse Gas Control 31:192-204.
- Gupta D, Carman PS, Venugopal R (2012) High-density brine-based fracturing fluid for ultra-deep fracturing stimulation. in SPE Annual Technical Conference and Exhibition.Society of Petroleum Engineers.
- Arabloo M, Shokrollahia A, Gharagheizib F, Mohammadi AH, et al., (2013)Toward a predictive model for estimating dew point pressure in gas condensate systems. Fuel Processing Technology 116: 317-324.
- Halali MA, Azaria V, Arablooa M, Mohammadib AH, Bahadorid A (2016) Application of a radial basis function neural network to estimate pressure gradient in water–oil pipelines. Journal of the Taiwan Institute of Chemical Engineers 58:189-202.
- Gibson RE,Loeffler OH (1949)Pressure -volume-temperature relations in solutions viii. The behavior of some solutions of electrolytes in water. Annals of the New York Academy of Sciences 51: 727-752.
- Anderson JH (1978)Multiple well dual fluid geothermal power cycle, Google Patents.
- Al Ghafri S, Maitland GC, Trusler JM (2012) Densities of Aqueous MgCl2 (aq), CaCl2 (aq), KI (aq), NaCl (aq), KCl (aq), AlCl3 (aq), and (0.964 NaCl+ 0.136 KCl)(aq) at Temperatures Between (283 and 472) K, Pressures up to 68.5 MPa, and Molalities up to 6 mol• kg–1. Journal of Chemical & Engineering Data 57:1288-1304.
- Kumar A (1986) Densities of aqueous strontium chloride solutions up to 200. degree. C and at 20 bar. J Chem Eng Data. 31:347-349.
- Phutela RC,Pitzer KS (1986) Densities and apparent molar volumes of aqueous magnesium sulfate and sodium sulfate to 473 K and 100 bar. J Chem Eng Data 31: 320-327.
- Obšil M, Majer V (1997) Densities and Apparent Molar Volumes of Na2SO4 (aq) and K2SO4 (aq) at Temperatures from 298 K to 573 K and at Pressures up to 30 MPa. Journal of Chemical & Engineering Data 42: 137-142.
- Gates JA, Wood RH (1985) Densities of aqueous solutions of sodium chloride, magnesium chloride, potassium chloride, sodium bromide, lithium chloride, and calcium chloride from 0.05 to 5.0 mol kg-1 and 0.1013 to 40 MPa at 298.15 K. J Chem Eng Data 30: 44-49.
- Crovetto R, Lvov SN, Wood RH (1993)Vapor pressures and densities ofNaCl (aq) and KCl (aq) at the temperature 623 K and CaCl2 (aq) at the temperatures 623 K and 643 K. The Journal of Chemical Thermodynamics 25: 127-138.
- Sharygin AV, Wood RH (1998) Densities of aqueous solutions of sodium carbonate and sodium bicarbonate at temperatures from (298 to 623) K and pressures to 28 MPa. The Journal of Chemical Thermodynamics 30:1555-1570.
- Isdale J, Morris R (1972) Physical properties of sea water solutions: density. Desalination 10: 329-339.
- Emerson W, Jamieson D (1967) Some physical properties of sea water in various concentrations. Desalination 3: 213-224.
- Chen CT, Millero FJ (1976)The specific volume of seawater at high pressures. Deep Sea Research and Oceanographic Abstracts 23: 595-612.
- Guliyev T (2009) The (
*r*,*r*,*T*) properties and apparent molar volumes*V*fof (ZnBr2+ C2H 5OH). The Journal of Chemical Thermodynamics 41: 1162-1169. - Safarov J(2012) (p, r, T) properties of seawater: Extensions to high salinities. Deep Sea Research Part I: Oceanographic Research Papers 65: 146-156.
- Phillips SL (1981) A technical databook for geothermal energy utilization. Lawrence Berkeley National Laboratory.
- Rogers P,Pitzer KS (1982) Volumetric properties of aqueous sodium chloride solutions. Journal of Physical and Chemical Reference Data 11: 15-81.
- Bahadori A, Al-Haddabi M, Vuthaluru H (2013) The estimation of reservoir brine properties during crude oil production using a simple predictive tool. Petroleum Science and Technology31: 691-701.
- Sorelle RR (1982) Mathematical field model predicts downhole density changes in static drilling fluids. SPE Annual Technical Conference and Exhibition, Society of Petroleum Engineers.
- Potter RW, Brown DL (1977) The Volumetric Properties of Aqueous Sodium Chloride Solutions from 0° to 500°C at Pressures Up to 2000 Bars Based on a Regression of Available Data in the Literature: Preliminary Steam Tables for NaCl Solutions, US Government Printing Office.
- Valle A, Utvik OH (1997) Pressure drop, flow pattern and slip for two phase crude oil/water flow: experiments and model predictions.Ichmt Digital Library Online.
- Brinkman H (1952)The viscosity of concentrated suspensions and solutions. The Journal of Chemical Physics20: 571-571.
- Himmelblau DM (2000) Applications of artificial neural networks in chemical engineering. Korean Journal of Chemical Engineering 17: 373-392.
- Chamkalani A (2013) Soft computing method for prediction of CO2 corrosion in flow lines based on neural network approach. Chemical Engineering Communications 200: 731-747.
- Tatar A, Shokrollahi A, Azari V(2015) A hybrid intelligent computational scheme for determination of refractive index of crude oil using SARA fraction analysis. The Canadian Journal of Chemical Engineering 93: 1547-1555.
- Ghiasi MM, Jaghdhan ZE, Lee M, Abbas A, Bahadori A et al., Development of soft computing methods to predict moisture content of natural gases. Journal of the Taiwan Institute of Chemical Engineers 55: 36-41.
- Lee S(1999) Measuring thermal conductivity of fluids containing oxide nanoparticles. Journal of Heat Transfer 121: 280-289.
- Huang Y (2003) Development of an artificial neural network model for predicting minimum miscibility pressure in CO< sub> 2 flooding. Journal of Petroleum Science and Engineering 37: 83-95.
- Zahedi G, Karami Z, Yaghoobi H (2009) Prediction of hydrate formation temperature by both statistical models and artificial neural network approaches. Energy Conversion and Management 50: 2052-2059.
- Elgibaly AA, Elkamel AM (1998) A new correlation for predicting hydrate formation conditions for various gas mixtures and inhibitors. Fluid Phase Equilibria 152: 23-42.
- Ghiasi MM (2012) Initial estimation of hydrate formation temperature of sweet natural gases based on new empirical correlation. Journal of Natural Gas Chemistry 21: 508-512.
- Chang F-J, Chang Y-T (2006) Adaptive neuro-fuzzy inference system for prediction of water level in reservoir. Advances in Water Resources 29: 1-10.
- Khajeh A, Modarress H, Rezaee B (2009) Application of adaptive neuro-fuzzy inference system for solubility prediction of carbon dioxide in polymers. Expert Systems with Applications 36: 5728-5732.
- Fayazi A, ArablooM, Shokrollahi A, ZargariMd-H, et al. (2013) State-of-the-Art Least Square Support Vector Machine Application for Accurate Determination of Natural Gas Viscosity. Industrial & Engineering Chemistry Research 53: 945-958.
- Fayazi A, Arabloo M, Mohammadi AH (2014)Efficient estimation of natural gas compressibility factor using a rigorous method. Journal of Natural Gas Science and Engineering 16: 8-17.
- Shokrollahi A, Tatar A, Safari H (2015)On accurate determination of PVT properties in crude oil systems: Committee machine intelligent system modelling approach. Journal of the Taiwan Institute of Chemical Engineers 55: 17-26.
- Ligtenberg J (2003) Unravelling the petroleum system by enhancing fluid migration paths in seismic data using a neural network based pattern recognition technique. Geofluids 3: 255-261.
- Karimpouli S, FathianpourN, Roohi J(2010) A new approach to improve neural networks' algorithm in permeability prediction of petroleum reservoirs using supervised committee machine neural network (SCMNN). Journal of Petroleum Science and Engineering 73: 227-232.
- Tatar A, Naserib S, BahadoricMd, Hezaved AZ, Kashiwao T, et al. (2016) Prediction of carbon dioxide solubility in ionic liquids using MLP and radial basis function (RBF) neural networks. Journal of the Taiwan Institute of Chemical Engineers 60: 151-164.
- Sayahi T, Tatar A, Bahrami M (2016) A RBF model for predicting the pool boiling behavior of nanofluids over a horizontal rod heater. International Journal of Thermal Sciences 99: 180-194.
- Naseri S, Tatar A, Shokrollahi A (2016) Development an accurate method to prognosticate choke flow coefficients for natural gas flow through nozzle and orifice type chokes. Flow Measurement and Instrumentation 48: 1-7.
- Naseri S, Tatar A, Bahadori M, Rozyn J, Kashiwao T, et al. (2015) Application of Radial Basis Function Neural Network for Prediction of Titration-Based Asphaltene Precipitation. Petroleum Science and Technology 33: 1875-1882.
- Shokrollahi A, Arabloo M, Gharagheizi F, Amir Md-H (2013) Intelligent model for prediction of CO2 – Reservoir oil minimum miscibility pressure. Fuel 112: 375-384.
- Farasat A, Shokrollahi A, Arabloo M, Gharagheizi F, Amir Md-H (2013)Toward an intelligent approach for determination of saturation pressure of crude oil. Fuel Processing Technology 115: 201-214.
- Fazavi M, Hosseini SM, Arabloo M, Shokrollahi A, Taleghani MN, et al. (2013) Applying a Smart Technique for Accurate Determination of Flowing Oil/Water Pressure Gradient in Horizontal Pipelines. Journal of Dispersion Science and Technology 35: 882-888.
- Rosa ES, Salgado RM, Ohishi T, Mastelari N(2010) Performance comparison of artificial neural networks and expert systems applied to flow pattern identification in vertical ascendant gas–liquid flows. International Journal of Multiphase Flow 36: 738-754.
- Kamari A, Gharagheizi F, Bahadori A, Mohammadi AH (2014) Determination of the equilibrated calcium carbonate (calcite) scaling in aqueous phase using a reliable approach. Journal of the Taiwan Institute of Chemical Engineers 45: 1307-1313.
- Safari H, Shokrollahi A, Jamialahmadi M, Ghazanfari MH, Bahadori A, et al. (2014) Prediction of the Aqueous Solubility of BaSO4Using Pitzer Ion Interaction Model and LSSVM Algorithm. Fluid Phase Equilibria 374: 48-62.
- Arabloo M, Ziaee H, Lee M, Bahadori A (2015) Prediction of the properties of brines using least squares support vector machine (LS-SVM) computational strategy. Journal of the Taiwan Institute of Chemical Engineers 50: 123-130.
- Kamari A, Bahadori A, Mohammadi AH(2015) On the determination of crude oil salt content: Application of robust modeling approaches. Journal of the Taiwan Institute of Chemical Engineers 55: 27-35.
- Tatar A,Naseri S, Lee M, Bahadori A,Sirach N (2015) Prediction of reservoir brine properties using radial basis function (RBF) neural network. Petroleum 1: 349-357.
- Tatar A, Naseri S, Bahadori M, Rozyn J, Lee M, et al. (2015) Evaluation of Different Artificial Intelligent Models to Predict Reservoir Formation Water Density. Petroleum Science and Technology 33: 1749-1756.
- Nouri-Taleghani M, Mahmoudifar M, Shokrollahi A, Tatar A, Karimi-Khaledi M, et al. (2015) Fracture density determination using a novel hybrid computational scheme: a case study on an Iranian Marun oil field reservoir. Journal of Geophysics and Engineering 12: 188.
- ShokrollahiA, Tatar A, Safari H(2015) On accurate determination of PVT properties in crude oil systems: Committee machine intelligent system modeling approach. Journal of the Taiwan Institute of Chemical Engineers.
- Tatar A, Shokrollahi A, Lee M, Kashiwao T, Bahadori A (2015) Prediction of supercritical CO2/brine relative permeability in sedimentary basins during carbon dioxide sequestration. Greenhouse Gases: Science and Technology 5: 756-771.
- Abido MA, Abdel-Magid YL(1997) On-line identification of synchronous machines using radial basis function neural networks. Power Systems, IEEE Transactions 12: 1500-1506.
- Majidi SMJ, Shokrollah A, Arabloo M, Soleymanloo RM, Masihi M (2014) Evolving an accurate model based on machine learning approach for prediction of dew-point pressure in gas condensate reservoirs. Chemical Engineering Research and Design 92: 891-902.
- Vapnik V (1995) The nature of statistical learning theory, Springer-Verlag, New York pp: 314.
- Baylar A, Hanbay D, Batan M(2009) Application of least square support vector machines in the prediction of aeration performance of plunging overfall jets from weirs. Expert Systems with Applications 36: 8368-8374.
- Übeyli˙ ED (2010)Least squares support vector machine employing model-based methods coefficients for analysis of EEG signals. Expert Systems with Applications 37: 233-239.
- Chen TS (2011) A Novel Knowledge Protection Technique Base on Support Vector Machine Model for Anti-classification. Electrical Engineering and Control 98: 517-524.
- Suykens JAK (2002) Least Squares Support Vector Machines, Singapore, World Scientific Pub. Co. Pte. Ltd.
- Suykens JAK, Vandewalle J (1999) Least Squares Support Vector Machine Classifiers. Neural Processing Letters9: 293-300.
- Ghiasi MM (2014) Robust Modeling Approach for Estimation of Compressibility Factor in Retrograde Gas Condensate Systems. Industrial & Engineering Chemistry Research 53: 12872-12887.
- Mahmoodi NM, Arabloo M, Abdi J (2014) Laccase immobilized manganese ferrite nanoparticle: Synthesis and LSSVM intelligent modeling of decolorization. Water Research67: 216-226.
- Wang H, Hu D (2005) Comparison of SVM and LS-SVM for regression. In International Conference on Neural Networks and Brain. Beijing, China: IEEE.
- Gunn SR (1998) Support Vector Machines for Classification and Regression,University of Southampton: Faculty of Engineering, Science and Mathematics School of Electronics and Computer Science.
- Müller KR, Mika S, Rätsch G, Tsuda K, Schölkopf B (2001) An introduction to kernel-based learning algorithms. IEEE Trans Neural Netw 12: 181-201.
- Haykin S (2012) Neural Networks and Learning Machines. 209: Pearson Education.
- Ranka S, Mohan CK, Mehrotra K, Menon A (1996) Characterization of a Class of Sigmoid Functions with Applications to Neural Networks. Neural Netw 9: 819-835.
- Rafiee-Taghanaki S (2013) Implementation of SVM framework to estimate PVT properties of reservoir oil. Fluid Phase Equilibria 346: 25-32.
- Ghiasi MM (2013) Novel methods predict equilibrium vapor methanol content during gas hydrate inhibition. Journal of Natural Gas Science and Engineering. 15: 69-75.
- Ghiasi MM, BahadoriA, Zendehboudi S (2014) Estimation of the water content of natural gas dried by solid calcium chloride dehydrator units. Fuel 117: 33-42.
- Cybenko G (1989) Approximation by superpositions of a sigmoidal function. Math Signals Systems 2: 303-314.

Citation: Tatar A, Halali MA, Mohammadi AH (2016) On the Estimation of the Density of Brine with an Extensive Range of Different Salts Compositions and Concentrations. J Thermodyn Catal 7: 167. Doi: 10.4172/2157-7544.1000167

Copyright: © 2016 Tatar A, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.