# In Situ Bioremediation of Crude Petroleum Oil Polluted Soil Using Mathematical Experimentation

Authors: Modupe Elizabeth Ojewumi; Moses Eterigho Emetere; Damilola Elizabeth Babatunde; Joshua Olusegun Okeniyi
Journal: International Journal of Chemical Engineering
 (2017)

---

## 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.


---

## Body

## 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 [7–12]. 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 1–4. 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.Table 1
Typical soil properties.
Soil type
Particle diameter (mm)
Porosity,εp
Permeability,k (m2)
Clay
<0.002
0.4–0.7
10−18–10−15
Silty
0.002–0.05
0.35–0.5
10−18–10−12
Sandy
0.05–2.00
0.25–0.5
10−14–10−10
Gravel
<2
0.25–0.4
10−18–10−10Table 2
Typical first-order biodegradation rate constants from field studies.
Field study location
Shoreline type
Oil type
Treatment
First-order biodegradable rate day−1
Reference
Alkanes
PAHs
Delaware
Sandy beach
Bonny light crude oil
Control
0.026
0.021
[12]
Nutrient
0.056
0.031
Inoculum
0.045
0.3026
Quebec, Canada
Tidal fresh water
Mesa light crude oil
Control
0.0028
0.0028
[19]
Nutrient
0.0023–0.034
0.0016–0.041
Texas
Brackish Wetland
Phase II: Arabian light crude oil
Control
0.019
0.017
[20]
Nutrient
0.042–0.061
0.018–0.027
Phase III: Arabian medium crude oil
Control
0.020
0.015
Nutrient
0.024
0.013
Inoculum
0.019–0.030
0.016–0.017Table 3
Relations for diffusion in porous solids [18].
Mechanism
Equations
Applies to
Bulk diffusion in pores
D
e
f
f
=
ɛ
p
D
τ
Gases or liquids in large pores
Knudsen diffusion
D
k
=
48.5
d
p
o
r
e
T
M
1
/
2in m2/s  Dkeff=εpDkτ  Ni=-DkdC1dz
Dilute (low pressure) gases in small poresNkn=λdpore>10
Combined bulk and Knudsen diffusion
D
e
f
f
=
1
-
α
x
A
D
e
f
f
+
1
D
k
e
f
f
-
1  α=1+NANB  Deff=1Deff+1Dkeff-1
N
A
≠
N
B    NA=NB
Surface diffusion
J
S
i
=
-
D
s
e
f
f
ρ
p
d
q
i
d
z    Dseff=εpDsτ  Dsθ=Dsθ=01-θ  Ds=-Dsqexp⁡-EsRT
Adsorbed gases and vapoursθ = fractional surface coverage ≤0.6Table 4
Maximum recorded growth rates for some organisms measured near their respective optimal temperature in complex media [16].
Organism
Temperature, °C
Doubling Time,td (h)
Vibrio natriegens
37
0.16
Bacillus stearothermophilus
60
0.43
Escherichia coli
40
0.38
Bacillus subtilis
40
0.43
Pseudomonas putida
30
0.75
Vibrio marinus
15
1.35
Rhodobacter sphaeroides
30
2.2Other Useful Data. Approximate predictions of Tortuosity factors for common adsorbents are as follows [18]:(1)τp=2-εp2εpMackieandMeares,τp=1εpWakaoandSmith,τp=εp+1.51-εpSuzukiandSmith.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:34.8×10-3 Pa·s.(c)
Approximately apparent diffusivity of oil sample:1.018×10-3 cm2/s.For this work, the following parameters were used:Soil type: sandy.Porosity,εp: 0.43.Tortuosity factor: 6.0.Growth rate organism:5.324×10-7/s (hydrocarbon utilizers).Soil initial concentration[C]=0.38g/g; 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 is(2)u=kpμΔϕ,where ϕ=∫Popdp/p+gz  (Hubbest).Equation (2) is generated as follows:For a one-dimensional flow of a fluid, it may be stated as(3)Q=AKDhZ.The differential form is shown below:(4)Q=AKdhdz=-AKdP/Pg+Zdz.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,(5)V=QA=-KdP/Pg+Zdz.For liquids other than water, the hydroid conductivity can be replaced by the permeability of the porous media, such that(6)V=KpgμdP/Pg+Zdz.This gives(7)V=-Kμ=dP+Zpgdz,where K=Kpg/μ and K is the permeability. Zpg is a constant and its differential is zero.Therefore (7) becomes (8)V=-Kμdpdz.Thus, for flow in the x, y, and z directions, respectively, Darcy’s law can be experienced.(9a)ux=Kxμdpdx,(9b)uy=-Kyμdpdy,(9c)uz=-Kzμdpdz+ρ,where u is volumetric rate of flows per unit cross-sectional area. μ is viscosity of the fluid. K is permeability (M2). ρ is density of the fluid (kg/m3). Z is height. ux is mass-average velocity in the x-direction (m/s). uy is mass-average velocity in the y-direction (m/s). uz is mass-average velocity in the z-direction (m/s). V 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,DAB, of a constituent A in solution is B, which is a measure of its diffusive mobility. (10)JA=-DAB∂CA∂z=-cDAB∂XA∂z.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 is(11)μk=µmaxSKS+S,where µmax is maximum growth rate achievable when S≫KS. S is substrate concentration. KS 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(12)R=µC,where R is model reaction term. µ is first-order degradation rate coefficient. C is substrate concentration.From Figure1, the continuity equation for flow could be written as(13)RateofoilEnteringelementalVolume-Rateofmassofoilleavingtheelementmass-Consumptionofoilbymicrobes=Themassrateofaccumulation,Rateofoilenteringtheelementalvolume-Consumptionofoilintheelementalvolume=Rateofaccumulationofoilintheelementalvolume.The 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.Figure 1
Elemental volume of oil contaminated soil.This is given by the summation of the net mass gains in theX, Y, and Z-direction.Let the elemental volume =ΔxΔyΔz = [M3].Let the molecular weight of oil =Mwo = [kg/mol].Let the molar flux of oil =No [mol/m2s].The mass rate of flow of componentA (oil) into the three faces with a common corner at A is therefore written as(14)Mw0N0xxΔyΔz+N0,yyΔxΔz+N0zΔxΔy,where N0,x signifies the X-directed flux and (N0,x) its value at location X.Similarly, the mass rate of flow out of the three faces with a common corner atB is(15)Mw0N0xx+ΔxΔyΔz+N0,yy+ΔyΔxΔz+N0zz+ΔzΔxΔy.Since the net rate oil entering the elemental volume is the difference between the entering and leaving oil flows rates (14) and (15) become(16)Mw0N0xx+Δx-N0,kkΔyΔz+N0,yy+Δy-N0yΔxΔy+N0,zz+Δz-N0,zzΔxΔy.If, in addition, the rate of biochemical reaction of oil in the element volume is given as(17)R0=Molvolumetime=MolM2S,therefore its production rate is(18)MwoR0ΔxΔyΔ2masstime.The rate of accumulation with the soil elemental volume is given by the following.If the total oil in the elemental volume =Po = ΔxΔyΔzPo, therefore its rate of accumulation is (19)ΔxΔyΔzdpodt.Combining (16), (18), and (19), dividing by ΔxΔyΔz, and taking the limit as the three distances become zero give(20)Mw0∂N0,x∂xΔxΔyΔz+∂N0,y∂yΔxΔyΔz+∂N0z∂2ΔxΔyΔz-MwoR0ΔxΔyΔz∂po∂t=ΔxΔyΔz∂po∂t,(21)Mw0∂N0,x∂x+∂N0y∂y+∂N0,z∂zΔxΔyΔz-MwoRo∂po∂tΔxΔyΔz=∂po∂tΔxΔyΔz,Mw0∂N0,x∂x+∂N0y∂y+∂N0,z∂zMwoRo∂po∂t=-∂po∂t∴Mw0∂N0,x∂x+∂N0y∂y+∂N0,z∂z+∂po∂t=-MwoR0.The molar flux of the oil N0 is made up of two parts: The adjusted amount or rate from the bulk motion Bo and the fraction XA of N which is Bo and resulting from diffusion JA(22)No=Bo+Jo.Equation (22) expressed in terms of mass in the X-direction gives(23)Mw0Nox=UxPo+Mw0Jox,where Ux is the mass-average velocityJ
0 is the diffusion flux given by Fick’s law(24)J0x=-Do∂∂x.Differentiating (23), we have(25)Mw0∂Noy∂x=Ux∂po∂x+Po∂Ux∂x+Mw0∂J0,x∂x.Similarly written for y and z-direction, we have (26)Mw0∂Noy∂y=Uy∂po∂y+Po∂Uy∂y+Mw0∂J0,y∂y,Mw0∂N0z∂z=Uz∂Po∂z+Po∂Uz∂z+Mwo∂Joz∂z.Substituting (25) and (26)(27)Po∂ux∂x+∂uy∂y+∂uz∂z+Ux∂po∂x+Uy∂po∂y+Uz∂po∂z︸a+Mw0∂Jox∂x+∂Jo,y∂y+∂Jo,z∂z︸b+∂po∂t︸c=-MwoRo_.Let us assume that density and flux are constants; then a, b, c = 0, so we have(28)Po∂ux∂x+∂uy∂y+∂uz∂z=-MwoRo,p0Mw0=C.Substituting for J0x,J0y&J0z in (24) and assuming that the diffusivity Do is equal to the x-, y-, and z-direction therefore (27) gives (29)Ux∂po∂x+Uy∂po∂y+Uz∂po∂z+Po∂uk∂x+∂uy∂y+∂uz∂z-MwoDo∂2C∂x2+∂2C∂y2+∂2C∂z2+∂po∂t=-MwoRo.Dividing by molecular weight, we have(30)Ux∂c∂x+Uy∂c∂y+Uz∂x∂z+C∂ux∂x+∂uy∂y+∂uz∂z-Do∂2C∂x2+∂2C∂y2+∂2C∂z2+∂C∂t=-Ro,(31)∂Ux∂x+∂Uy∂y+∂Uz∂z=0Continuityequation∴Ux∂c∂x+Uy∂c∂y+Uz∂c∂z-Do∂2C∂x2+∂2C∂y2+∂2C∂z2+∂C∂t=-Ro.Assuming that velocity = 0, there is no chemical reactionEquation (31) reduces to Fick’s law:31a∂c∂t=Do∂2C∂x2+∂2C∂y2+∂2C∂z2.Assuming a unidirectional diffusion in a particular dimension is taken as z-direction, (31) becomes(32)Uz∂c∂z-Do∂2C∂z2+∂C∂t=-Ro,∂c∂z-Do∂2C∂z2+Uz∂C∂t-Ro,where Ro 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,31bμ=μmaxCKo+C,where Ko is saturation constant.The obtained mathematical model is31c∂c∂t=D=∂2C∂Z2-Uz∂C∂z-μC.


## 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 is(33)∂C∂t=D∂2C∂z2-uz∂C∂z-μC.The expression above is an initial value problem in terms of time, t, but boundary conditions can be generated for the variable z.Initial conditions are(34)C=C0,z=0,t=0;C=0,z>0,t=0.Boundary conditions are(35)D∂C∂tz=0,t>0=uzC0-uzC+D∂2C∂z2z=L,t>0=0,where C0 is initial oil concentration and L is total height.Defining the dimensionless variables,(36)aa=DL2,ab=-uzL,ac=-μ,so that the model equation (21) becomes(37)dc-dt=aa∂2c-∂z2+ab∂c-∂z+acc-,where c-=C/C0 is relative (dimensionless) concentration.Considering the nature of the boundary condition,D(∂C/∂z)=0, at z=0, the transformation below is used z=z-2(38)dz=2zdz.So that: (39a)∂c-∂z=∂c-∂z∗∂z∂z(39b)=2z∂c-∂z(39c)=2√z∂c-∂z.Also,(40a)∂2c-∂z2=∂∂z∂c-∂z.This, from (39c), becomes(40b)∂∂z2√z∂c-∂z=212√z∂c-∂z+√z∂2c-∂z2=2√z∂2c-∂z2+1√z∂c-∂z.Substituting (39c) and (40b) into (37) yields(41a)dc-dt=2aa√z∂2c-∂z2+aa√z∂c-∂z+2ab√z∂c-∂z+acc-.This can be factorized to give(41b)dc-dt=2aa√z∂2c-∂z2+√zaaz+2ab∂c-∂z+acc-.The coefficients in this expression can be by single terms such that (42)a1=2aaz,a2=zaaz+2ab,a3=ac,to give new ODE (43)dC-dt=a1∂2C-∂z2+a2∂C-∂z+a3C-.Rendering the initial conditions (at t=0) dimensionless yields(44a)C-=1,z-=0,t=0;C-=0,z->0,t=0.Rendering the other boundary conditions (at z=L) dimensionless yields(44b)ad∂C-∂z-=C--1,where ad=-aa/ab.


### 4.2. Discretization of Model

Orthogonal collocation is applied to boundary value problem. In this case, “z” is bounded but “t” is not. Carrying out discretization, the following transformations are used:41a1∂C-∂z=∑1N+1Ai·jC-j,∂2C-∂z2=∑1N+1Bi.jC-j.Substituting these expressions into (43) yields (45)dC-idt=ai∑1N+1BijC-j+a2∑1N+1AijC-j+a3C-i.This can be written as (46a)dC-idt=+a1∑1NBijC-j+a1Bi,N+1C-N+1+a2∑1NAijC-j+a2Ai,N+1C-N+1+a3C-i.And by grouping like terms gives(46b)dC-idt∑1Na1Bij+a2AijC-j+a1Bi,N+1+a2Ai,N+1C-N+1+a3C-i.Discretizing the other boundary condition (at z=L) yields (47a)a4∑1N+1AN+1,jC-j=C-N+1-1,where a4=2abz (from (39c) and (44a)).As in (46a), (47a) can be expanded as(47b)a4∑1NAN+1,jC-j+a4AN+1,N+1C-N+1=C-N+1-1,so that(48)c-N+1=a4∑1NAN+1,jC-j+11-a4AN+1,N+1,which when substituted into (46b) gives the expression(49)dC-idt=∑1NaiBij+a2AiJC-j+aiBi,N+1+a2Ai,N+1a4∑1NAN+1,jC-j+11-a4AN+1,N+1+a3c-i.In simple terms, the expression above can be written as(50)dC-idt=∑1NfijC-j+ei+a3C-i,where(51)fij=aiBij+a2Aija4a1Bi,N+1+a2Ai,N+1AN+1,j1-a4AN+1,N+1,ei=a1Bi,N+1+a2Ai,N+11-a4AN+1,N+1.As a matrix, (48) could be shown as(52)C1C2MCN=f11f12Kf1Nf21∧∧f2NM∧∧MfN1∧∧fNNC1C2MCN+e1MMeN+a3C1C2MCN.Rewriting this in vector form yields (53)dCdt=f+1a3C+e,dCdt=f∗C+e.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.


## 4.1. Dimensionless Parameter Determination

The obtained mathematical model is(33)∂C∂t=D∂2C∂z2-uz∂C∂z-μC.The expression above is an initial value problem in terms of time, t, but boundary conditions can be generated for the variable z.Initial conditions are(34)C=C0,z=0,t=0;C=0,z>0,t=0.Boundary conditions are(35)D∂C∂tz=0,t>0=uzC0-uzC+D∂2C∂z2z=L,t>0=0,where C0 is initial oil concentration and L is total height.Defining the dimensionless variables,(36)aa=DL2,ab=-uzL,ac=-μ,so that the model equation (21) becomes(37)dc-dt=aa∂2c-∂z2+ab∂c-∂z+acc-,where c-=C/C0 is relative (dimensionless) concentration.Considering the nature of the boundary condition,D(∂C/∂z)=0, at z=0, the transformation below is used z=z-2(38)dz=2zdz.So that: (39a)∂c-∂z=∂c-∂z∗∂z∂z(39b)=2z∂c-∂z(39c)=2√z∂c-∂z.Also,(40a)∂2c-∂z2=∂∂z∂c-∂z.This, from (39c), becomes(40b)∂∂z2√z∂c-∂z=212√z∂c-∂z+√z∂2c-∂z2=2√z∂2c-∂z2+1√z∂c-∂z.Substituting (39c) and (40b) into (37) yields(41a)dc-dt=2aa√z∂2c-∂z2+aa√z∂c-∂z+2ab√z∂c-∂z+acc-.This can be factorized to give(41b)dc-dt=2aa√z∂2c-∂z2+√zaaz+2ab∂c-∂z+acc-.The coefficients in this expression can be by single terms such that (42)a1=2aaz,a2=zaaz+2ab,a3=ac,to give new ODE (43)dC-dt=a1∂2C-∂z2+a2∂C-∂z+a3C-.Rendering the initial conditions (at t=0) dimensionless yields(44a)C-=1,z-=0,t=0;C-=0,z->0,t=0.Rendering the other boundary conditions (at z=L) dimensionless yields(44b)ad∂C-∂z-=C--1,where ad=-aa/ab.


## 4.2. Discretization of Model

Orthogonal collocation is applied to boundary value problem. In this case, “z” is bounded but “t” is not. Carrying out discretization, the following transformations are used:41a1∂C-∂z=∑1N+1Ai·jC-j,∂2C-∂z2=∑1N+1Bi.jC-j.Substituting these expressions into (43) yields (45)dC-idt=ai∑1N+1BijC-j+a2∑1N+1AijC-j+a3C-i.This can be written as (46a)dC-idt=+a1∑1NBijC-j+a1Bi,N+1C-N+1+a2∑1NAijC-j+a2Ai,N+1C-N+1+a3C-i.And by grouping like terms gives(46b)dC-idt∑1Na1Bij+a2AijC-j+a1Bi,N+1+a2Ai,N+1C-N+1+a3C-i.Discretizing the other boundary condition (at z=L) yields (47a)a4∑1N+1AN+1,jC-j=C-N+1-1,where a4=2abz (from (39c) and (44a)).As in (46a), (47a) can be expanded as(47b)a4∑1NAN+1,jC-j+a4AN+1,N+1C-N+1=C-N+1-1,so that(48)c-N+1=a4∑1NAN+1,jC-j+11-a4AN+1,N+1,which when substituted into (46b) gives the expression(49)dC-idt=∑1NaiBij+a2AiJC-j+aiBi,N+1+a2Ai,N+1a4∑1NAN+1,jC-j+11-a4AN+1,N+1+a3c-i.In simple terms, the expression above can be written as(50)dC-idt=∑1NfijC-j+ei+a3C-i,where(51)fij=aiBij+a2Aija4a1Bi,N+1+a2Ai,N+1AN+1,j1-a4AN+1,N+1,ei=a1Bi,N+1+a2Ai,N+11-a4AN+1,N+1.As a matrix, (48) could be shown as(52)C1C2MCN=f11f12Kf1Nf21∧∧f2NM∧∧MfN1∧∧fNNC1C2MCN+e1MMeN+a3C1C2MCN.Rewriting this in vector form yields (53)dCdt=f+1a3C+e,dCdt=f∗C+e.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 z-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.Figure 2
Possibilities of the bioremediation process: (a) Solution 1 atD=0.1, U=0.25; (b) Solution 1 at D=0.1, U=0.45; (c) Solution 2 at D=0.1, U=0.25; (d) Solution 2 at D=0.1, U=0.45; (e) Solution 3 at D=0.1, U=0.25; (f) Solution 3 at D=0.1, U=0.45.
(a)
(b)
(c)
(d)
(e)
(f)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 (Figures2-3).Figure 3
Comparative analysis of experimental and simulated results.Also, noticeable from Figure4, 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.Figure 4
Comparative analysis of experimental and simulated results at Ce/0.38 (dimensionless).Figure 5
Comparative analysis of experimental and simulated results on soil depth.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. Figure6 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.Figure 6
Simulated results of varying soil depths.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.



---

*Source: 5184760-2017-07-04.xml*