^{1}

^{*}

^{2}

^{1}

This study aims to develop a mathematical analysis for one-dimensional modeling of a radial flow through a production well drilled in a confined aquifer, in the case of steady-state flow conditions. An analytical solution has derived from that expression for estimation of drawdowns according to different flowrates. Through that process, the evaluation of static pressure, the calculation of hydraulic charge due to the waterflow through the well is evaluated, the drawdowns curves are drawn and at last, the obtained curves are analyzed. The curves obtained for the different flow rates have an asymptotic direction, the axis of the hydraulic charges. The variation of the hydraulic charge depends on the radial distance for different flow rates. The P point, is a common point of all curves obtained for different production flowrates in the well. This point is where the well production flowrate is optimum for the optimal hydraulic charge.

Groundwater flows from the interconnections of aquifers to the producing well by radial flow (waterflow between aquifers and wells) and obeys to the physical phenomenon based on the relevant physical principles: Darcy’s law and mass balance, (Equations (1) and (2) respectively) which are fundamental equations for developing theories in groundwater flow; they are:

Darcy’s law:

v = − K ∇ h (1)

and the mass balance

d i v − ( K ∇ h ) = 0 (2)

where K = k ⋅ ρ w ⋅ g μ and h = z h + p ρ w ⋅ g .

This developed theory of flow modeling process is applicable to several areas of flow mechanics, i.e. geothermal reservoirs, hydrodynamics, and mining process.

During the waterflow through the well, the hydraulic charge in the confined aquifer varies from the hydraulic charge h_{R}, due to the well-range R, to the hydraulic charge h_{W} due to the radius of the well a. The present study develops first, the calculation of the drawdown, therefore, the drawdown curve (cone of depression) in the potentiometric surface around the producing well basing on Equation (1), considering that the flow regime in confined aquifer is steady-state, the hydraulic charge at any point of the aquifer remains constant, i.e. does not change with the time; thus ( ∂ h / ∂ t ) = 0 ; and the fluid velocity is independent of the time. The measured data at each time (the well-range depended variable), is a point of the drawdown curve (cone of depression), in the potentiometric surface created by the hydraulic charge (in the well). The developed mathematical model determines the expression of the static pressure which can be expressed by the pressure in the well and the hydraulic charge in the piezometer (the difference between the static level and the dynamic pumping water level) and defines the variation of the drawdown cone in the potentiometric surface. The model aims to calculate the different values of the hydraulic charge in function of the radial distance (r), with a ≤ r ≤ R , according to different values of flowrate at the wellhead.

Secondly, it develops a case study of the appearance of groundwater drawdown curves and analyzes through the developed mathematical model, their variation according of the wellhead flowrates during the water production process and their different curves.

The objective of the present paper is to determinate the static pressure of the confined aquifer and function of the pressure at the bottom of the drilling, and then simulate the variation of the hydraulic charge of the groundwater in function of the well action radius R, and finally to simulate mathematically the potentiometric surface of the depression cone created by the variation of the radial distance (r) and the variation of the flow rates (Q) at the wellhead, which describes the different levels attended from the static pressure to the hydraulic pressure at the bottom of the well.

The physical model of the problem is constituted by a confined aquifer, a porous medium which thickness is ; and a producing drilling whose flow rate is Q. The well cross-section is simulated to a circular section which radius is a. In the cylindrical coordinates system ( O , r , z ) ; h R is the hydraulic charge at the limit of the action radius R in the confined aquifer, and h W is the hydraulic charge of the well; the reservoir, the porous medium which thickness is e; and the producing well which flowrate is Q and h R are the hydraulic charge at the limit of the well-range R in the reservoir [

The well-range is the locus point for a part of potentiometric surface curve where, the pressure does not influence the potentiometric surface in the reservoir during the groundwater exploitation, comparing to its axis, and the value of the depression cone (pumping-out), i.e. the drawdown is not available; ( h R − h w ) = 0 .

The following assumptions for the governing equations are the following:

- the porous medium is isotropic, homogenous;

- permeability, porosity (ϕ), hydraulic conductivity (K), and productive thickness (e), are supposed to be constant;

- the density (ρ), the dynamic viscosity (μ) and the cinematic viscosity (ϑ) of the groundwater are considered as constants.

The water is produced in accordance to

The boundary conditions for the governing equations are the following.

The governing equations are [

The generalized law of Darcy:

v = − K ∇ h with K = k ⋅ ρ w ⋅ g μ (3)

in the case of an incompressible fluid, h = z h + p ρ w ⋅ g

The continuity equation:

∂ ∂ t ( ϕ ) + d i v ( v ) = 0 (4)

The diffusivity equation:

∂ ( ϕ ) ∂ t + d i v ( − K ∇ h ) = 0 (5)

Considering the porosity ϕ independent of time, the porous medium non-deformable, taking only the hydrodynamic part and the temperature constant, ρ depending only on p, we got:

∂ ( ϕ ) ∂ t = ∂ ( ϕ ) ∂ p ∂ p ∂ t = s ρ ⋅ g ∂ p ∂ t (6)

where s = ρ w g ∂ ( ϕ ) ∂ p is the storage-specific coefficient which shows the ability

of the porous medium to release fluid when increasing the pressure. By introducing Equation (6) in Equation (4), we obtain:

s ρ w ⋅ g ∂ p ∂ t + d i v ( − K ∇ h ) = 0 (7)

where h = p ρ w g + z h and ∂ h ∂ t = ( ∂ p ∂ t ρ w g − p g ∂ ρ w ∂ t ) 1 ( ρ w g ) 2

Since the fluid is incompressible, thus, ∂ ρ w ∂ t = 0 ; and ∂ p ∂ t ≈ ρ w g ∂ h ∂ t ; we got the mass balance equation for an incompressible fluid in a non-deformable medium:

s ∂ h ∂ t + d i v ( − K ∇ h ) = 0 (8)

the diffusivity equation for the confined groundwater aquifer.

In steady flow condition, ∂ h ∂ t = 0 and thus:

d i v ( − K ∇ h ) = 0 (9)

This Equation (9) allows to estimate the (hydraulic charge) potentiometric surface created hydraulic charge expression h ( r ) in function of the hydrogeological parameters. In cylindrical coordinates, the diffusivity equation can be written as follows:

∂ 2 h ∂ r 2 + 1 r ∂ h ∂ r + 1 r 2 ∂ 2 h ∂ θ 2 + ∂ 2 h ∂ z 2 = 0 (10)

As the flow is axisymmetric flow, ∂ 2 h ∂ θ 2 = 0 and there is no horizontal flow so ∂ 2 h ∂ z 2 = 0 .

Thus the previous differential equation is reduced to:

∂ 2 h ∂ r 2 + 1 r ∂ h ∂ r = 0 and becomes d 2 h d r 2 + 1 r d h d r = 0 (11)

Therefore 1 r d d r ( r d h d r ) = 0

The integration yields:

r d h d r = C 1 ; d h = C 1 r d r ; and h = C 1 ln r + C 2

Using the boundary conditions:

r = a ; h w = C 1 ln a + C 2 , r = R ; h R = C 1 ln R + C 2

with C 1 = h R − h w ln ( R / a ) and C 2 = h w − h R − h w ln ( R / a ) ln a

The general solution is then:

h = h R − h w ln ( R / a ) ln r + h w − h R − h w ln ( R / a ) ln a ;

thus

h = h R − h w ln ( R / a ) ln ( r a ) + h w (12)

The flow velocity into the well is calculated by using Darcy’s equation and taking into account the previous boundary conditions:

v = − K d h d r = − K h R − h w ln ( R / a ) 1 r

The flow rate can be determined by evaluating the flow at any radial distance r and integrating over the flow surface. Taking into account r = a , we get:

v = − K h R − h w ln ( R / a ) 1 a

Q = ∫ A ( v ) r = a d A = ∫ θ = 0 2π − K ln h R − h w ln ( R / a ) 1 a d θ

Q = − 2 π K e h R − h w ln ( R / a )

Thus

h R − h w = Q 2 π K e ( ln R a ) (13)

According to Laurent et al. [

relation P h q = μ Q 2 π e k ln ( R a ) , and they supposed that ln ( R a ) = 2 π can be used in this case. By replacing Equation (13) in Equation (12) we obtain:

h = Q 2 π K e ( ln r a ) + h w (14)

The hydraulic charge of the well related to the bottom pressure of well is expressed by

h w = p w ρ w ⋅ g (15)

According to Codo, [

following relationship: P w = P s t − μ Q 2 π e K ln ( R a ) ; with ln ( R a ) = 2 π ; so we get:

P w = P s t − μ Q e K (16)

Replacing Equation (16) into Equation (15), then we get:

h w = p s t ρ w ⋅ g − Q e . K (17)

with K = k ρ g μ ; Substituting Equation (17) into Equation (14) we obtained the expression of the hydraulic charge; h = Q 2 π K e ( ln r a ) + P s t ρ w g − Q e K [

Thus the expression of h becomes:

h = Q 2 π K e [ ( ln r a ) − 2 π ] + P s t ρ w ⋅ g

and at least

h = Q 2 π K e ( ln r a ) − Q e K + P s t ρ w ⋅ g (18)

Let us remember that the hydraulic charge h for an incompressible fluid is:

h = v 2 2 g + p ρ w + z h

In the porous medium, the groundwater flow can be neglected ( v = 0 ) thus, v 2 2 g = 0 and then the hydraulic charge is equivalent to the static charge:

h = p ρ w ⋅ g + z h

The static pressure of the water at any point of the aquifer is defined with the equation of Bernoulli [

- The reservoir is at a constant temperature of 20˚C, and so negligible,

- The aquifer is such that the operating range (R) is very large compared to its thickness (e),

- The porous medium is homogeneous, permeable and isotropic ( k h = k v = k ),

- The fluid is incompressible ( ρ w = c s t e ),

- The flow is radial from the reservoir to the drilling under the following conditions: ∂ 2 h / ∂ θ 2 = 0 et ∂ 2 h / ∂ z 2 = 0

- The flow is two-dimensional flow ( x , y ) and ( v z = 0 ) ,

- The flow is steady ( ∂ h ∂ t = 0 ),

- The flow is conservative ( d i v v = 0 ).

At the wellhead, the pressure can be expressed as [

p 2 = p 1 − ρ w ⋅ g ⋅ H − λ ⋅ ρ w H v 2 2 d (17)

with the bottom pressure p 1 = p s t − μ Q 2 π e k ln ( R a ) ; λ = 0.086 R e 0.2 ; R e = v d ϑ and v = 4 Q π d 2 .

After transformations

p 2 = p s t − ρ w ⋅ g ⋅ H − μ Q 2 π ⋅ e ⋅ K ln ( R a ) − 0.06642 ρ w H d 4.8 Q 1.8 ϑ 0.2 (18)

with the approximation of Laurent et al. [

p s t = p 2 + ρ w ⋅ g ⋅ H + μ Q e ⋅ K + 0.06642 ρ w H d 4.8 Q 1.8 ϑ 0.2 (19)

The calculation of the static pressure assumed the use of the parameters of the well and aquifer. Just before the pumping operation ( t = 0 ), the hydraulic charge is uniformly distributed throughout the aquifer, the initial condition is h ( x , y , z , o ) = h R ; ∀ ( x , y , z ) .

The two boundary conditions associated with the problem are the following:

- At the aquifer: r = R , the charge is h R

- At the well: r = a , the charge is h W

From the test pumping carried out on the aquifer of Monzoungoudo by DGEau, we obtained some values of hydrodynamic parameters whose are the following:

- Temperature of the well is constant at 20˚C,

- Water density ρ W = 1000 kg / m 3 .

- Dynamic viscosity μ = 10 − 3 Pa ⋅ s .

- Kinematic viscosity ϑ = 10 − 6 m / s .

Geometric and hydrodynamics data of Monzoungoudo well, are provided and collected in

From the pumping test of the aquifer, the hydrodynamic parameters obtained are the followings:

- Permeability: 2.28 × 10^{−4} m/s.

- Average porosity: 30%.

And the following table shows the geometric and the hydrodynamic characteristics of the Monzoungoudo well.

The following

The flowrates vary from 0.002 m^{3}/s to 0.009 m^{3}/s and the representative curves are simulated to a logarithmic function. The common point of the curves are determined by transforming the Equation (13), and then, we obtain the following relation:

Symbols | Description | Values | Units |
---|---|---|---|

H | Depth of the well | 244.18 | m |

D | Diameter of the well | 0.126 | m |

Q | Flow rate | 2000 | cm3/s |

g | Acceleration of the gravity | 9.81 | m/s2 |

p 2 | Pressure at the head of the well | 4.16 | bar |

ε | Absolute roughness of the pipe | 0.12 | mm |

[ h − p s t ρ w ⋅ g ] − Q 2 π e K [ ln ( r a ) − 2 π ] = 0 , ∀ Q ∈ ℝ (20)

This equation is equal to zero, if simultaneously, we have:

{ h − P s t ρ w g = 0 ln ( r a ) − 2 π = 0 (21)

thus: { h = P s t ρ w g r = a e 2 π .

Using the following data:

P s t = 28.12 bars ; ρ w = 1000 kg / m 3 ; g = 9.81 N / kg ; a = D / 2 = 0.063 m

We get: { h = 28.12 × 10 5 1000 × 9.81 r = 0.063 e 2 π .

At least, we obtain a P point { h = 286.65 m r = 33.74 m .

All the curves are logarithmic curves and all go across, at a common point P ( h = 286.65 m ; r = 33.74 m ). The curves appear as two parts of the curve in two intervals created by this P point; thus, there are the first interval 0.0063 W r W 33.74 and the second one 33.74 W r WR where R is the well-range.

In the first interval, no matter the flow rates, all curves get a common part where all the curves are identic and mixed; and in the second one, the curves are different, in accordance to water flowrates.

Observing curves, we note that, the drawdown (s) tends to zero (by adopting as asymptotic direction, the horizontal direction i.e. the static pressure piezometric surface level) while the radial distance tends to well-range.

- In the case of 0.063 m ≤ r < 33.74 m , the vertical asymptotic direction, is the hydraulic charges axis. In this zone, the drawdown is appreciable until the radial distance r = 33.74 m .

- At the P point (33.74 m; 286.65m), the drawdown is identic no matter the value of flowrate; when r = 33.74 m , all hydraulic charge curves intersect regardless of the value of the flowrate; thus; at this point, the common optimum value of drawdown s o p t is available.

- When r > 33.74 m , for higher production flow rate, the closer representative hydraulic charge curve is in the undisturbed aquifer charge (initial piezometric surface). The water production does not influence the initial piezometric surface. When the production rate increases, then the corresponding drawdowns are low compared to those of previous flow rates. In the same way, the different drawdowns observed for the different hydraulic charge curves (according to the different flow rates) decrease when the radial distance (r) increases and are cancelled while the radius of action (R) of the well is reached.

This paper attended to the mathematical steady-state radial modelling simulation of flow around the production well in the confined aquifer of Monzoungoudo. All the curves, obtained for the different flow rates, have as an asymptotic direction, the axis of the hydraulic charges. The variation of the hydraulic charge is a function of the radial distance for different flow rates. The result shows that, on one hand, the drawdown occurs significantly in the close proximity of the well, for variable values of flow; on the other hand, the obtained point P ( h = 286.65 m ; r = 33.74 m ) is the crossing point of all curves for different production flowrates in the well. It is where we get the optimum of production flowrate for the optimal hydraulic charge in a steady-state flow condition in a confined water table.

This work was supported by ICMPA-UNESCO Chair and LSTE.

We, authors declare that we have no competing interests.

de P. Codo, F., Hountondji, B. and Aina, M.P. (2019) Steady-State Radial Flow Modeling through the Production Well in the Confined Aquifer of Monzoungoudo, Benin. Open Journal of Fluid Dynamics, 9, 107-118. https://doi.org/10.4236/ojfd.2019.92008

a: radius of the well

R: well-range

r: radial distance ( a ≤ r ≤ R )

h_{R}: hydraulic charge in the aquifer

h_{W}: hydraulic charge in the well

p: pressure

p_{st}: static pressure

P_{w}: bottom pressure of the well

R_{e}: Reynolds number

h: hydraulic charge

ρ_{w}: density of the water

g: acceleration of the gravity

z_{h}: altitude

K: hydraulic conductivity of the medium

μ: dynamic viscosity of the fluid

ϑ: kinematic viscosity of the fluid

k: intrinsic permeability of the medium

ϕ: porosity of the medium

v : darcy velocity

P_{hq}: hydrodynamic flow pressure,

s: drawdown

S: storage specific coefficient

H: depth of the well

e: productive thickness of sheet

Q: water flowrate

λ: coefficient of discharge