Abstract

Mathematical modelling of in situ (on site) bioremediation of crude petroleum polluted soil was investigated. An unsteady state mathematical model based on bulk flow of oil through the soil and molecular diffusion through the pores of the soil was developed. The parabolic partial differential equation developed was resolved into a system of ordinary differential equations (ODEs) by orthogonal collocation method and the necessary boundary condition was used. The resultant system of ODE was solved using fourth-order Runge-Kutta method. The simulated data gave a good agreement with experimental data.

1. Introduction

One of the major concerns of oil industry today is how to improve the recovery of large percentage of oil remaining unrecovered in the old and new depleted producing fields. Microbial treatment of oil well system and reservoir to control paraffin deposition is one of the efficient enhanced oil recovery methods [1]. Bioremediation project in practical sense is very expensive to execute. The accepted disposal methods of incineration or burial of insecure landfills are now expensive when amounts of contaminants are large. Mechanical and chemical methods generally used to remove hydrocarbons from contaminated sites have limited effectiveness and can be expensive [2]. The aftermath of bioremediation activities may be dangerous to life in the predetermined region. Hence, the first step towards achieving a safe bioremediation process is developing mathematical models that comprise all necessary parameters that determine the entire process.

References [3, 4] reported that, in Taiwan, over 400 gas stations had encountered problems of contamination from petroleum hydrocarbons, due to aged pipelines. Leaks and accidental spills occur regularly during the exploration, production, refining, transport, and storage of petroleum and petroleum products [3]. Release of hydrocarbons into the environment whether accidentally or due to human activities is a main cause of water and soil pollution [4]. This is one of the major environmental problems today resulting from the activities of the petrochemical industry. Hydrocarbon components have been known to belong to the family of carcinogens and potent immunotoxic organic pollutants which are harmful to all living things both terrestrial and aquatic [5]. Soil contamination with hydrocarbons causes extensive damage of local system since accumulation of pollutants in animals and plant tissue may cause death or mutations [6]. Among the many techniques employed to decontaminate the affected sites, in situ bioremediation using indigenous microorganisms is by far the most widely used [712]. This approach to reclaiming contaminated land reduces the threat to groundwater and enhances the rate of biodegradation [13]. The technology commonly used for the soil remediation includes mechanical, burying, evaporation, dispersion, and washing. However, these technologies are expensive and can lead to incomplete decomposition of contaminants [2, 6, 14, 15].

Mathematical model is regarded as a decision tool that assists decision makers in effectively dealing with complex issues such as oil spillage on soil surfaces [1]. When oil is released into the subsurface, it travels downward until it encounters a low permeability layer that it cannot penetrate because of large capillary force. If oil spill is present in sufficient quantity the oil may pool on the low permeability materials or more latterly, following the geological gradient of this materials contact of oil with groundwater results in dissolution of the soluble content of oil and subsequent contamination of the groundwater [16].

Most mathematical models follow a mathematical structure that can be used to describe and study a real solution. Hence, mathematical model is a representation in mathematical terms of certain aspects of nonmathematic system. The arts and crafts of mathematical modelling are exhibited in the construction of models that not only are consistent in themselves and mirror the behaviours of their prototype, but serve some exterior purpose.

2. Data Collection

In this paper it is intended to describe the bioremediation of land on which oil has spilled using a mathematical model such that it is possible to predict the residual concentration of the oil at any time. The model was divided into parts which includes the following physical principles [17], that is,(i)law of conservation of mass (continuity equation),(ii)Darcy’s diffusion in liquids (Fick’s laws),(iii)Monod kinetics.In a porous fluid flow medium, the most significant quantity that must be conserved is the mass. In the case of oil biodegradation, no oil is generated; only consumption by microbes takes place. We shall adopt experimentally verified parameters shown in Tables 14. Table 1 gives the soil properties of the polluted soil that shall be assumed. Table 2 illustrates the first-order biodegradation rate constants from field studies. Table 3 shows the relations for diffusion in porous solids. Table 4 reveals the maximum recorded growth rates for some organisms measured near their respective optimal temperature in complex media.

Other Useful Data. Approximate predictions of Tortuosity factors for common adsorbents are as follows [18]:This means that, for sandy soil, with irregular particle and pores shape and size, the tortuosity value will be much higher than the highest one from these correlations:(a)Specific gravity of crude oil (sour blend): 0.842.(b)Dynamic viscosity:  Pas.(c)Approximately apparent diffusivity of oil sample:  cm2/s.For this work, the following parameters were used:Soil type: sandy.Porosity, : 0.43.Tortuosity factor: 6.0.Growth rate organism: (hydrocarbon utilizers).Soil initial concentration ; soil concentration in g/g soil can be converted to g/m3 by multiplying the former by soil bulk density.

3. Theories of Formulation

Under this section, we considered the full formulations of Fick’s and Darcy’s law. Darcy’s law expresses the fact that the volumetric flow rates per unit cross-sectional area at any point in a uniform porous medium are proportional to the gradient in potential in the direction of flow at that point [17]. The law is valid for laminar flow at low Reynold’s number and its mathematical expression iswhere   (Hubbest).

Equation (2) is generated as follows:

For a one-dimensional flow of a fluid, it may be stated asThe differential form is shown below:The negative sign shows that the head decreases in the direction of flow.

Dividing (4) by the flow area will give Darcy’s law; thus,For liquids other than water, the hydroid conductivity can be replaced by the permeability of the porous media, such thatThis giveswhere and is the permeability. is a constant and its differential is zero.

Therefore (7) becomes Thus, for flow in the , , and directions, respectively, Darcy’s law can be experienced.where is volumetric rate of flows per unit cross-sectional area. is viscosity of the fluid. is permeability (M2). is density of the fluid (kg/m3). is height. is mass-average velocity in the -direction (m/s). is mass-average velocity in the -direction (m/s). is mass-average velocity in the -direction (m/s). is Darcy’s flux (m/s).

Fick’s law deals with movement of individual molecules through a substance by virtue of their thermal energy. Diffusion occurs in substances as a result of concentration gradient; materials move from regions of higher concentration to regions of lower concentration to attain uniformity in substances. The rate at which a solute moves at any point in any direction must therefore depend on this concentration gradient at that point and in that direction. This rate is measured in a direction normal to the diffusion.

The diffusivity or diffusion coefficient, , of a constituent in solution is , which is a measure of its diffusive mobility. The negative signs emphasize that diffusions occur in the direction of a drop in concentration [17].

Monod Kinetic. Monod kinetics describe(a)Enzymes-substrate activity(b)Adsorption isotherm behaviours(c)Microbial growth rate(d)Microbial uptake of substrateThe Monod equation iswhere is maximum growth rate achievable when . is substrate concentration. is value of the concentration of the concentration at which the specific growth rate has half its value.

Simple first-order degradation models could be written as where is model reaction term. is first-order degradation rate coefficient. is substrate concentration.

From Figure 1, the continuity equation for flow could be written asThe net rate oil entering the elemental volume equals difference between the rate of oil entering the elements and the rate of mass of oil leaving the elemental volume.

This is given by the summation of the net mass gains in the , , and -direction.

Let the elemental volume = = [M3].

Let the molecular weight of oil = = [kg/mol].

Let the molar flux of oil = [mol/m2].

The mass rate of flow of component (oil) into the three faces with a common corner at is therefore written aswhere signifies the -directed flux and () its value at location .

Similarly, the mass rate of flow out of the three faces with a common corner at isSince the net rate oil entering the elemental volume is the difference between the entering and leaving oil flows rates (14) and (15) becomeIf, in addition, the rate of biochemical reaction of oil in the element volume is given astherefore its production rate isThe rate of accumulation with the soil elemental volume is given by the following.

If the total oil in the elemental volume = = , therefore its rate of accumulation is Combining (16), (18), and (19), dividing by , and taking the limit as the three distances become zero giveThe molar flux of the oil is made up of two parts: The adjusted amount or rate from the bulk motion and the fraction of which is and resulting from diffusion Equation (22) expressed in terms of mass in the -direction giveswhere is the mass-average velocity

is the diffusion flux given by Fick’s lawDifferentiating (23), we haveSimilarly written for and -direction, we have Substituting (25) and (26)Let us assume that density and flux are constants; then , , = 0, so we haveSubstituting for in (24) and assuming that the diffusivity is equal to the -, -, and -direction therefore (27) gives Dividing by molecular weight, we haveAssuming that velocity = 0, there is no chemical reaction

Equation (31) reduces to Fick’s law:Assuming a unidirectional diffusion in a particular dimension is taken as -direction, (31) becomeswhere is the rate of consumption of oil by microorganism which can be equivalent to “.” is the growth rate of the microorganism.

Monod kinetics describes the growth rate of organism; that is,where is saturation constant.

The obtained mathematical model is

4. Model Decomposition

The methods used for this study are listed below:(i)Conversion of variables to dimensionless form.(ii)Discretion of model.

4.1. Dimensionless Parameter Determination

The obtained mathematical model isThe expression above is an initial value problem in terms of time, , but boundary conditions can be generated for the variable .

Initial conditions areBoundary conditions arewhere is initial oil concentration and is total height.

Defining the dimensionless variables,so that the model equation (21) becomeswhere is relative (dimensionless) concentration.

Considering the nature of the boundary condition, , at , the transformation below is used So that: Also,This, from (39c), becomesSubstituting (39c) and (40b) into (37) yieldsThis can be factorized to giveThe coefficients in this expression can be by single terms such that to give new ODE Rendering the initial conditions (at ) dimensionless yieldsRendering the other boundary conditions (at ) dimensionless yieldswhere .

4.2. Discretization of Model

Orthogonal collocation is applied to boundary value problem. In this case, “” is bounded but “” is not. Carrying out discretization, the following transformations are used:Substituting these expressions into (43) yields This can be written as And by grouping like terms givesDiscretizing the other boundary condition (at ) yields where (from (39c) and (44a)).

As in (46a), (47a) can be expanded asso thatwhich when substituted into (46b) gives the expressionIn simple terms, the expression above can be written aswhereAs a matrix, (48) could be shown asRewriting this in vector form yields The resultant equation above representing a system of N-coupled ordinary differential equations with initial conditions stated in (34) can be solved using a numerical integration method such as the fourth-order Runge-Kutta method.

5. Results and Discussion

Using (33) we obtained three different type solutions as shown in Figures 2(a)2(f). The solution represents the various possibilities in the bioremediation process. Figure 2(a) shows how the residual oil concentration is expected to reduce with time. This solution is trivial because it expresses the normal scenario at the completion of a successful field work. Figure 2(b) shows what happens when the mass-average velocity in the -direction is increased from 0.25 m/s to 0.45 m/s. At shorter time, the residual oil concentration is expected to reduce. However, at different scenario, for example, considering a heterogeneous soil sample, the concentration may increase with time at a defined soil depth (Figures 2(c) and 2(d)). However, at varying or selective batch treatment (Figures 2(e) and 2(f)), there could be a more effective treatment. Hence we applied the model to practical or experimental results.

The residual substrate concentration reduces with time at any depth in the soil clearly indicating the processes considered in the model development (bulk flow, molecular diffusion, and microbial degradation). While experimental values show reduction from 0.38 g/g soils to 0.14/g soils after 36 days, simulated values drop from 0.38 g/g soil to 0.08 g/g soil. These values yield comparable overall disappearance rates of 0.006 g/g soil per day and 0.0083 g/g soil per day, respectively (Figures 2-3).

Also, noticeable from Figure 4, the concentration of residual oil at soil depths closer to the surface rises rapidly initially before displaying the trend of concentration drop at the surface. This is as a result of more rapid accumulation in these layers than in layers farther from the surface. This trend reflects the effect of the soil properties and distance on the oil flow as aptly described by the model’s boundary conditions. Reasonably, Figure 5 indicates a very small amount of oil in deeper layers of soil at the earlier time that gradually builds up over time with increasing bulk flow.

For the soil type considered in the experiment used as comparison (sandy soil), the model displays similar trend in the range of oil velocity values permissible for the soil type (varying only soil permeability). This further confirms the model developed and the computer program package for its solution, accurate tools for its solution, and accurate tools for predicting the fate of crude spilled on soil. Figure 6 shows that the concentration of residual oil at the soil depths closer to the surface rises rapidly initially before displaying the trend of concentration drop at the surface. This is as a result of more rapid accumulation in these layers than layers farther from the surface.

This trend reflects the effect of the soil properties and distance on the oil flow as described by the model’s boundary conditions.

6. Conclusion

Bioremediation provides the best environment-friendly approach to treating oil pollution by spillage. Requirements for success in carrying out this technique of remediation include adequate knowledge of the soil characteristics nature of microorganisms that degrade crude oil (whether present in the soil or not) and ample time. Because soil types and microbial growth varies with pollutant concentration and time, bioremediation of polluted site is usually specific to the site.

With adequate data on the soil and required microbe characteristics, the right mathematical description of the process can be developed, as achieved in this work. This can be used to predict the fate of the oil with time and the extent of pollution; to predict time when complete remediation would have taken place; and to investigate possible methods of increasing the rate of bioremediation. Ultimately, the cost of bioremediation is reduced when the above roles of model development and solution are carried out.

The model developed in the course of this work suits only a simple case where no temperature effects come to play and underground water is not affected by the oil spill. The availability of all relevant experimental data plays a crucial part in comparing trend by simulation with that obtained from experiment. In the case of bioremediation, the characteristics of the particular soil type being remediated, microbial growth rates at prevailing environmental conditions, and relevant properties of the particular oil type spilled are necessary to prove the reliability and robustness of the model and its solution method.

Variables

:Permeability (m2)
:Diffusion coefficient (length2/time or cm2/sec.)
:Knudsen diffusion (in m2/s)
:Density of the fluid (kg/m3)
:Velocity of the fluid (m/s)
:Characteristic linear dimension (m)
:Dynamic viscosity of the fluid with respect to object (Pas or Ns/m2 or kg/m·s)
:Volumetric flow rate (m3/s)
:Pipe’s cross-sectional area (m2)
:Doubling Time (h)
:Mass-average velocity in the -direction (m/s)
:Mass-average velocity in the -direction (m/s)
:Mass-average velocity in the -direction (m/s)
:Darcy’s flux (m/s)
:Maximum growth rate (time−1).

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors appreciate the partial sponsorship of Covenant University, Ota, Nigeria.