Magazine of Civil Engineering issn
2712-8172
journal homepage: http://engstroy.spbstu.ru/
Research article UDC 66.067.1 DOI: 10.34910/MCE.119.2
Numerical study of the process of unsteady flow in a three-layer
porous medium
R. Ravshanov1, Z.S. Abdullaev2©, E.V. Kotov313, Sh.N. Turkmanova 2©
1 Research Institute for the Development of Digital Technologies and Artificial Intelligence, Tashkent, Uzbekistan
2 Tashkent Institute of Irrigation and Agricultural Mechanization Engineers, Tashkent, Uzbekistan
3 Peter the Great St. Petersburg Polytechnic University, St. Petersburg, Russia M ekotov. cfd@gmail. com
Keywords: mathematical model, groundwater, groundwater, filtration, analytical solution, numerical algorithm
Abstract. The unsteady fluid flow in a three-layer porous medium is numerically investigated and is an important and topical problem. An analytical solution of the equation for the pressure fluid layer is obtained on the basis of the theory of elastic regime, taking into account the overflow from the coating and the low-permeability layer into the low-permeability bulkhead and external sources that greatly affect the liquid level change. In the main bounded aquifer only horizontal liquid migrations prevail, and in the cover and low-conductive layers only vertical migrations prevail, allowing for horizontal components of the flow rate to be omitted here. Evaporation from the surface of a liquid in a porous medium is considered. Evaporation from the surface of a liquid in a porous medium significantly affects the distribution of overflows in the layers and the distribution of the liquid level in a porous medium. Therefore, when designing vertical drains for enhanced oil recovery in multilayer reservoirs and designing fluid flows in reservoirs, evaporation must be taken into account. The computational experiments have established that the dynamics of changes in the liquid level in a porous medium depends significantly on the evaporation parameter. With an increase in its value, the liquid level in a porous medium decreases proportionally over time. The accuracy of the numerical solution using the balance equation showed that the error does not exceed 1.3%.
Funding: This work is supported by the Russian Science Foundation under grant 21-79-10283, date 29 July 2021.
Citation: Ravshanov, R. Abdullaev, Z.S., Kotov, E.V., Turkmanova, Sh.N. Numerical study of the process of unsteady flow in a three-layer porous medium. Magazine of Civil Engineering. 2023. 119(3). Article no. 11902. DOI: 10.34910/MCE.119.2
1. Introduction
The pressureless flow of fluid through a multilayer porous medium is found in many technical applications. Historically, such flows were considered concerning hydraulic structures. Later, the obtained results and approaches were applied to filter technology, chemical processes and apparatuses, oil production, climate technology, and the analysis of fluid movement in the elements of buildings and structures.
In particular, the paper [1] proposes a mathematical model and a numerical algorithm for monitoring and forecasting groundwater and surface water migrations using geofiltration models.
© Ravshanov, R. Abdullaev, Z.S., Kotov, E.V., Turkmanova, Sh.N., 2023. Published by Peter the Great St. Petersburg Polytechnic University.
The problem analytical solution aims to study the interaction of surface and subsurface water flows. The paper [2] considers two-dimensional steady subsurface water flow in the vertical plane. In that article, the aquifer is idealized as an infinite band, and the channel is modeled as a horizontal equipotential function.
The article [3] proposes a mathematical model and a numerical algorithm for solving the problem of unsteady free flow of groundwater filtration considering well galleries in heterogeneous porous media where the wells differ by their hydrogeological characteristics. The authors of article have compiled an analytical solution to the specified problem to get a linearized system describing groundwater free filtration using the Laplace transform with variable t. To create isolated areas preventing the spread of harmful liquids and to protect groundwater in the interlayers, ratios to support the water heads there have been derived, as well as formulas for determining interlayer water tables in the corresponding zones of the groundwater filtration area [4].
In the article [5], the authors have developed a general mass transfer model based on the infiltration model of shore intake considering water exchange between groundwater and surface water. The mass transfer model describes the salt transfer in groundwater and the kinetics of salt exchange in dry soil. Some attention in modeling the process is paid to developing a methodology for numerical modeling of the groundwater level regime in the zones of shore intake influence. This will ensure reliable calculations when forecasting the operation of wells and estimating groundwater reserves of water intakes in difficult hydrogeological conditions. Numerical algorithms and software are developed to perform calculations on computer systems to solve inverse problems of aquifers' hydrogeological parameters and identify water exchange parameters. The general mass transfer model proposed by the authors of the articles [6, 7] will make it possible to consider the main processes comprehensively and to assess the groundwater mineralization degree and the transfer of contaminants, including hydrocarbon ones, by interacting flows of surface water and groundwater.
The article [8] deals with filtration processes: with a stationary filtration mode on an interfluvial soil body with constant levels; with a nonstationary filtration mode at backwater on an interfluvial soil body; influenced by a rise in the water level on the right border. Depression curves are made for stationary and nonstationary modes, and the correctness of computer program calculations has been proved using the Dupuit formulation. As noted in the paper, the convergence of the actual level values and the model ones has been found unsatisfactory at levels of 0.40 m and 0.25 m and when applying the stationary mode. However, should the level increase to 0.35 m, there are almost no or small discrepancies. The numerical model demonstrates how the level change rate depends on the transmissivity level: the higher the transmissivity level, the faster the stationary filtration mode starts running.
The article [9] provides the results of mathematical modeling of the water table distribution in the underflow talik depending on the intensity of water intake from wells and pits. As concluded by the article's authors on the calculation results, the selected layout of wells along the talik zone of the river and the water intake regime will provide the water volumes required for the development of the mining and processing plant. it is necessary to operate five additional production wells To ensure the estimated demand for drinking water during 2013-2017. The hydrodynamic impact of the pit on the designed water intake wells is also insignificant. By the 31st (2044) year of development, an additional drop in the water level from 0.8 m (well 42) to 8.9 m (well 127) is forecasted, which does not exceed the permissible 65 m.
A mathematical model and a numerical algorithm are developed to forecast the groundwater table on the slope of the river valley [10] To provide hydrogeological forecasts. The article's authors note that it characterizes soils' geological structure, underground flow boundary conditions, and water-conductivity parameters. Hydrogeological parameters and boundary conditions intensity degree are identified by the multivariate numerical modeling. Using the results obtained from it, water flows to drainage structures, and their dependence on technogenic infiltration have been estimated. Parameters of shallow red-brown clays are determined at which the clays promote the local rise of the groundwater level. The article considers the hydrodynamic processes in the mouth of the Temernik River on the right bank of the Don River to give detailed data on the soil body's hydrogeological structure and assess the filtration parameters and boundary conditions. The authors used the results to determine the water flow to the in-situ drainage structures in the seepage zone. In the studies performed by the authors, the method used is based on numerical hydrogeological modeling that systematically includes the interrelated geofiltration parameters and makes it possible to cover the geoecological risk factors to develop effective solutions to combat flooding.
In the paper [11], the author modeled subsurface water filtration through a homogeneous earth dam with vertical slopes on a non-conductive base. The GEO-SLOPE GeoStudio hydrodynamic calculation software package has identified the main characteristics of the filtration flow and constructed the depression curves. To confirm the adequacy of the results of the numerical calculations, the author of the work has compared them with the results calculated using a proven method; the conclusion provides recommendations for using boundary conditions to make modeling more accurate.
The paper [12] considers steady plain filtration in a rectangular dam with a partially non-conductive vertical wall during evaporation from the water table. A mixed multiparameter boundary equation has been formulated for the analytical function theory to study the evaporation effect; it was solved using the method of P.Ya. Polubarinova-Cochina. The proposed model provides the basis for an algorithm for calculating the filtration characteristics of the flow. The results of a hydrodynamic analysis of the flow rate dependencies are presented, as well as the ordinate of the point where the depression curve starts concerning all physical parameters of the scheme. The exact values of the specified characteristics are compared with the known approximate values obtained by other authors out of evaporation process conditions. The results of the study give an idea (at least a qualitative one) on the possible dependence of the migration characteristics when considering the problem of filtration to an imperfect pump well.
The article [13-14] describes a numerical study of the problems of free nonlinear filtration in a trapezoidal and rectangular soil dam with the horizontal layers of different coefficients: filtration, partial saturation zone, vertical dam core, and horizontal drain canals. A generalization to the three-dimensional case is given as well. The proposed grid calculation method can be used to conduct multifactorial studies of the environmental aspects affecting water transfer in in-situ multilayer porous media.
The article [15] compares the solutions for two methods of calculating anisotropic dam filtration; numerical experiments have been performed for several profiles of soil dams, their elements, and antifiltration devices (downstream or upstream shell, screen, and core). Further, the following methods of solving anisotropic problems have been compared: through an imitated hydrodynamic filtration grid made by stretching an orthogonal hydrodynamic grid previously constructed by the EGDA method for a distorted isotropic model of the dam; and through a finite element numerical method supported by the local variation method.
The paper [16] considers the numerical solution of the anisotropic filtration problem for the soil dam. The results of computational experiments are given that were performed on a computer when solving the filtration problem and calculating the stability of the soil dam slope taking into account the anisotropy.
In the dissertation [17], the implementation methods are developed and compared concerning onedimensional models of water runoff along river slopes and the channel network based on the finite difference methods and the finite element method; the behavior patterns are proposed, implemented, and studied for the difference schemes of numerical integration of two-dimensional models of water runoff along the surface of slopes with a topography of different complexity; the numerical integration method has been developed and verified for equations of vertical water transfer in soils based on a four-point implicit scheme; the model of water erosion during rainfall floods is proposed and implemented, which describes the processes of drip and plane erosion and the transfer of soil particles by water flow along the surface of river slopes and the channel network. The advantages of applying finite element schemes are also demonstrated concerning real catchments; effective algorithms for its application have been developed, various methods of combining models of rainwater runoff formation at various catchment schemes are proposed.
In the dissertation [18], a steam-assisted gravity drainage model is developed for the first time taking into account the law of filtration with the ultimate pressure gradient, which allows for describing the main development stages of the steam chamber, namely, its growth to the top of the stratum, horizontal expansion, and expansion of the steam chamber towards the stratum bottom; a method has been developed for determining the anisotropic stratum filtration parameters according to vertical interference test; a mathematical model has been developed for studying the stationary fluid flow to the radial system of horizontal wells in the anisotropic stratum taking into account the influence of hydraulic pressure losses on friction in wellbores; a semi-analytical model is created to describe the process of nonstationary fluid flow to a multisectional horizontal well equipped with inflow control valves and pressure sensors in isolated sections; short-time tests are designed for vertical wells with hydraulic fracture, imperfect wells, and horizontal wells; complex transfer functions are expressed, and amplitude-frequency and phase-frequency specifications are determined for the following systems: 'porous-fractured vertical well' and 'hydraulic fracture in layer of finite conductivity.'
According to the results obtained, many conceptual mathematical and computational models have been developed to forecast the process of groundwater filtration and migration of variable-saturated liquids in multilayer porous media. In contrast to the above works, in this research, the unsteady liquid flow in a three-layer porous media considers external sources that greatly affect the liquid level change.
2. Methods
For a numerical study, specific conditions have been considered. The unsteady flow to the vertical drain canals in a three-layer boundary stratum considers evaporation from the liquid stream's upper surface and the elastic regime in a low-conductive layer. The main bounded aquifer lies under a low-conductive cover stratum and has a low-conductive layer below, facilitating its connection with the bed rock. It is
assumed that the A.N. Myatiyev-G.N. Girinskii theory is true for these conditions. It should be noted that in the main bounded aquifer only horizontal liquid migrations prevail, and in the cover and low-conductive layers only vertical migrations prevail, allowing for horizontal components of the flow rate to be omitted here.
Considering the above, the continuity equation for the cover layer has the following form:
dV
dz
— 0.
(1)
For this purpose, V is the vertical component of the flowrate in the porous media.
Further, applying the Darcy Law and considering the condition of the continuity of the liquid heads on the cover layer bottom, gives the following:
Me
dH1
dt
K„
H - H1 H1
+ q + F;
(2)
and the following equation for the bounded aquifer based on the elastic regime theory including liquid migration from the cover and low-conductive layer
dH df^dH Л df^dH}
dy
M---— —I T--1+---
dt dx V dx) dy
T
V
+K
H1 - H
Hi
- K
dH2 (x, y, -m, t)
dz
(3)
with a low-conductive dam taking into account the elastic filtration mode will be written as follows:
dH 2
a
d 2 H 2
dt
dz2
(4)
Here, Hi is a liquid head in the top layer; це is a free liquid loss or lack of saturation; Ke is a filtration coefficient in the cover layer; H is a liquid head in the intralayer; q is the total infiltration characterizing the actual infiltration and evaporation from the liquid level in a porous medium; x and y are coordinates of the horizontal plane, л is an elastic liquid loss coefficient; F is an external source; T = mk is filtration conductivity; K is a filtration coefficient; m is the thickness of the middle aquifer layer; Kn is a filtration
.. , . Knmn
coefficient; H2(x,y,z,t) is a liquid head in the lower layer; z is a vertical coordinate; an ~ ~ is a
Mn
piezoelectric conductivity coefficient; Цп is an elastic liquid loss coefficient; mn is the thickness of the low-conductive layer (Fig. 1).
3. Results and Discussion
In contrast with the study performed by F.B. Abutaliyev [19] and other authors, the external F source on the cover layer of the porus media was considered.
Figure 1. Schematic model of the flow in the three-layer porous medium.
Thus, the system of differential equations in partial derivatives (2)-(4) describes the flow process in a three-layer porous medium.
It should be noted that in the article [16], an analytical solution to the problem was given, with the linearization of equations (2) and (3) and condition q = const. However, it is necessary to emphasize that q(x,y,Hi,t) is a function of coordinates x, y of the level (Hi) and time (t).
This function depends on Hi as follows. If the liquid level in a porous medium is on the upper horizontal surface of a porous medium, then the function takes maximum values. If the liquid level in a porous medium drops below a critical depth, then the function equals zero.
This monotonically decreasing function reaches its maximum value when Hi coincides with the upper horizontal surface of a porous medium and asymptotically tends to zero at Hi ^Hkp, where Hkp is a critical value of the liquid level in a porous medium. Thus, q = 0 below this critical depth. The value q in such assumptions denotes evaporation from the liquid level in a porous medium.
In the paper [20], the following relation is used to calculate the value q:
q = q0
V mb - H1 ^
v mb - HkpJ
(5)
Here, qo is the intensity of evaporation on the upper horizontal surface of a porous medium; Hkp is the critical subsurface liquid depth; n is an exponent that depends on external factors and the subsurface liquid depth.
It should be noted that, in the article [21], evaporation from the liquid level in a porous medium for a single-layer stratum model was considered at n = 1.
According to the study, the general case of the evaporation task in form (4) for the multilayer stratum model has no solution. It should be noted that, in general, it is difficult and rather impossible to formulate an analytical solution to the problem of unsteady filtration for the system (2)-(3). Therefore, it is advisable to apply the finite difference method to integrate the nonlinear system of equations (2)-(4), taking into account the evaporation in the form (5), and the external source F.
Let us consider the liquid flow to the vertical drainage well drilled into the main aquifer in a limited circular stratum with r <Rk, considering the evaporation in the form (5)
It is assumed that the well is located in the stratum center. Then, due to the flow symmetry in (2)-(4), the system of equations considering the liquid intake to the upper surface of the liquid flow:
F + q0
f TJ Y
1- mb - H1
kp J
V
mb -H>
■Mb
H = K H1
dt '
H - H
1 dH = 1 d_(r H1 + Kb H1 -H Kb dH2 (r, -m, t)
a dt r dr V dr J T H1
1 dH2 _ d2 H
T
dz
an dt dz
With initial and boundary conditions based on the following:
H1 (r,0) = H (r,0) = H2 (r, z,0) = H0;
SH (R, t)_ Qc ,
dr
2nTRc
dH (Rk, t)
dr
= 0;
H 2 (r, -m, -mn, t ) = H 0; H 2 (r, -m, t ) = H (z, t).
(6)
(7)
(8)
(9)
(10)
(11)
(12)
(13)
To integrate the system (6)-(8) numerically into the conditions (9)-(13), the nondimensional variables by the following formulas proceeded:
U = -H^; V = -^~; W = -H2-; s = l n—; z = mbZ;
H_ p H™ p H„ p Rk
- xap
xap 2
- xap
R
t = -^~ t; Q = 2nTHxflpQ .
a
Together with the equations (6)-(13), the equations became as following:
dU U-V
F + P
f \n
(1- 1 -ш
V
1 -U
kp J
dV=i-2s FL+aUzL-в
дт ds2 U
-Y--
дт
dW
U
(14)
f
s, -
m Л
— ,t
mb J
dZ
5
dW d2W
дт dZ ’
U ( s,0) = V ( s,0) = W ( s,Z,0) = W0;
dV ( sc ,t)
ds
dV ( 0,t)
= Q ;
= 0;
(15)
(16)
(17)
(18)
ds
(19)
W
f m + J
s,-------
v mb
_ W0;
(20)
where
W
f m >
s,-----,t
v mb
_ V ( s,t) ;
a _
KbR2.
r\ SJ
R_ Kn Rk . .._ a^b H . s_ amb . Hxap
P m ; Y __ _ 1 Hkar; О 2 ; ^
anRk
THkap" T mb’- KR2 kr anRk* mb
(21)
p0 _ *L; W0 _Jb0-; sc _ &Д ; O' _-0
0 Kb 0 Hxap c Rk 2nTH
xap
According to the analysis of material balance numerical solutions, it is found that the nondimensional spatial variable s must be taken in the form indicated above since only in this case, the solution outcomes
n r
in the vicinity of the well Rc are correctly considered. For example, if s _ —, then the error in the balance
Rk
ratio can be about 20 %.
An implicit finite difference scheme is applied to numerically solve the system of nonlinear equations (14)-(16) with additional conditions (16)-(21).
The segment (so, 0) was split into m equal parts with an increment of As. Then Si = so + iAs, i = 0, 1,2, ..., m.
Due to condition (21), equation (16) must be solved m - 1 time along the lines parallel to the z-axis passing through the points si, I = 0, 1, 2, ..., m-1. To solve this equation with the finite difference method, the segment
f
V
m m + m,
Л
mh
m
b J
is to be divided into l equal segments with the points
m + m
Zj _-
m
+ jAZ; j = 0, 1, 2, ..., l, AZ = const.
’b
The uniform time increment At was introduced. As for the system of equations (14)-(16) with any of the spatial and temporal points {si,Zj,kAT}, i = 0, 1, 2, ..., m-1, j = 0, 1, 2, ..., l - 1, k = 1, 2, ..., a stable implicit finite difference scheme with accuracy was compiled as follows
(As )2 +(AZ)2 + At
V,k -V,k-1 _0l~2Si (V+1,k -2V,k + V-1,k) + aArUiU V,k +eA-(-3Wt
Uj ь 2AZ
l,k + 4Wi,l-1,k - Wi,l-2,k ); (22)
W, jk -Wi, jk-1 _ x(Wi, j+1,k -2W-, j,k +W-, H,k);
(23)
AZa-Uk-1 )_P0
f1-1 -ли1Л Y _ U , k -V, k
1 -Ukp J U,
(24)
i,k
Here
в
Ат
(А?
; Ui,k = U (s,кАт); x
Ат
M2
У,к = V (s, кАт); ^v,k = W (s ,^, кАг);
j = 1,2,..., t -1 i = 1,2,...,m -1 к = 1,2,...
Considering (20), the equation (23) was written as follows
Wi, j, к = Л', j+1,kWi, j+1, к + Bi, j+1,к,
(25)
where the sweep coefficients are defined as follows:
Ai, j+1, к
x
1 + (2 - 4, j к )x
i, j+1, к
W', j к-1 + xBJ, j к .
1 + (2 - Ai, j к )x
4,1,к = 0; Bi,1,к = W0; j = 1,2,..., t -1.
In the equation (22), W was substituted with its value from (25), taking into account the conditions (18) and (22), assuming Uj,к known. The equation will be rewritten in the following form
Vi к = Ci+1,к4+1, к + Di+1,к. (26)
Where
Ci+1, к
et~2SJ
1 + (2 - С1к )et-2 Sj +^1-^1 1 ,к) Uty 2AZ
(4 - Ai,t-1,к) Ai,t,к - 3
Vj к -1 + аАт-
D,
i +1, к
в Ат
2AZ
[( 4 - Ai,t- 1,к ) Bi,t,к - Bi,t-1,к
+ et~2SjDi, к
1 + (2 - Cj к )er2S +аАт_вАт[(4 - A,t-1, к ) A t, к - 3
Ui к 2АЪ
„2S,
Сг, к = 1 -
2в
1 +
аАт вАт
U1 2АС
[( 4 - Ai,t—1,к ) Ai,t,к - 3
Ак = -А$(2 +
t2^|V,к-1 +аАт + 2^[(4- Ai,t-1,к)Bi,t,к -Bi,t—1,к
i = 1,2,..., m -1.
Out of equation (26), considering condition (19), the formula got like this:
Vm, к ^m- 1,к
D
т,к
1 - C
т,к
Thus, once С^,к and Dj,к are calculated according to the recurrent ratios above and Vm-2,h is sequently determined, it will be possible to finally solve the following based on the formula (25)
W, j,к : Witt,к,..., WjJ,к.
(27)
Vo,к
It is assumed that Uj,h Then, for example, assuming that Uj,к = Щк-1, the first approximation Vi/1 = Wjjk(1'> is known. The UjJ1 is found by substituting VmJ1 into equation (14) and integrating it, for
example, by the Adams-Stormer method over [(k - 1)Ar, kAr\. The second approximation Vt,k(2) is determined by it substituted into (26).
The process will be finished when the following iterative process reproducibility condition is met
max
-V
(p)
ik
<S,
where e > 0 is a small value of the calculation error.
Analysis of the numerical calculations performed according to the above algorithm demonstrates that it is advisable to apply damping according to the sequence Vi,kpP with regard to the formula
Vi- k ^ = ( k(p) +(1 -u)V;k(P-1), (28)
where Vi,k(p) is the solution of (2.33). That is, if the sequence Vi,k(p) tends to be 'inconsistent,' the equation (28) stabilizes it.
To assess the accuracy of the numerical solution, the method of balance equations is used. This equation is deriveted as follows. The equation (7) can be written in nondimensional form
sdV d fsdV)
£^- = ^7 £77 + aPo£
дт d£V d£)
f
1 -
1-XU
1 - U
kp )
* dU
-ay£- p£-
дт
dW
f „ m ^
£,-----,т
mb
д£
,(29)
where £ = —
Rk
The equation (29) is integrated throughout the pore space of the aquifer stratum. Then, due to m = const, the integrated form is as follows
1
\£dTd£=-Q‘+ap> f£1 - 7
1
1 -XU
£ £
The notation was introduced:
- U
kP )
Y 1 dW
d£-ayf£---
g
( m ^
£,----,т
v mb )
d£
d£.
(30)
Vcp(T)=f£V (£,т) d£;
£c
Ucp(T)=f£U (£,т) d£;
£c
Then equation (30) can be written as follows
dV
cP
dUc
dt
= -Q + qucn(т)-ау ~tPpq
dt
For this purpose,
(31)
(32)
(33)
qucn (т) = аРо f£
£c
1-
1 -XU
n
1-U
d£;
f
dW
qn (т) = Pf£-
£c
kP ) \
„ m
£,-----,т
v mb )
д£
d£.
(34)
(35)
Integrating (33) by т gives the follows
Vcp (7) = Vcp (0) -Q*t + ]qucn (r)dr-ay[Ucp (r) -Ucp (°)]- Jqn J)J
(36)
That the ratio is called the balance equation. The following is proceed to assess the accuracy of the numerical solution to this equation. The numerical solution corresponding to the moment of time (т) is averaged by formulas (31) and (32). Then the evaporation and liquid migration is calculated by formulas (34) and (35) further substituting these values in (36). Then an approximate equation denoting the accuracy of the approximate solution averaged over the entire porous space is as follows
Vcp (T)+Q*t+arUcp (T)- J qucn (T) dT+Jqn (T)dT *Vcp (°)- arUcp (°) •
(37)
The problem is considered with the following data to illustrate the abovementioned algorithm:
Ro = 10 cm, Rk = 500 m, mb = 40 m, m = 100 m, mn = 10 m, H° = 39 m, n = 0, 1,2, 3;
q° = 0, 0.0036, 0.036 m/day T = 103 m2/day, ^b = 0.1,
Hkp = 37 m
Qc = 1,250 m3/day, a = 106 m2/day, an = 102 m2/day,
K = 0 m/day,
Kn = 0.01 m/day Kb = 1 m/day
In this example, the vertical filtration coefficients vary greatly. The ratio of the filtration coefficients in
k
the cover and high-conductive layers is — = 1° , while the ratio of these coefficients in the high-conductive
K6
k
layer and the low-conductive dam is — = 1, °°° . Until recently, this fact suggested that low-conductive
K6
layers containing liquid are incompressible or, at best, have elastic reserves. However, it is assumed them to be completely negligible. Studies [22] have shown that layers can give a significant amount of liquid, even the low-conductive ones. The liquid saved against other reserves reduces the efficiency of vertical drain canals (wells), increasing the running time at the liquid level in a porous medium.
Figures 2-4 show the results of calculations regarding the liquid level in the cover layer for different parameters n and evaporation qo values. These figures illustrate how the parameter n and evaporation qo affect the distribution of the liquid level in a porous medium.
Figure 2. Change in the value of the liquid level in a porous medium at the cover layer at n = 1
and different evaporation values qo.
According to the curves in Fig. 2, the change in the liquid table value significantly depends on the evaporation value. As the evaporation value increases, the liquid level in a porous medium drops proportionally over time.
Figure 3. Change in the value of the liquid level in a porous medium in the cover layer at n = 2
and different evaporation values qo.
According to the numerical calculations obtained and the curves in Fig. 3 and 4, it can be seen that the change in the value of the liquid level in a porous medium also depends on the change in the parameter n. As it increases, the liquid level in a porous medium drop at different evaporation values [23].
Figure 4. Change in the value of the liquid level in a porous medium in the cover layer at n = 3
and different evaporation values qo.
Fig. 5-7 represent flow curves reduced to an area unit with a cover layer and a low-conductive layer for different evaporation parameters qo and n values. Liquid migrations depend on evaporation parameters.
Figure 5. Flow curves reduced to an area unit with a cover and low-conductive layers at qo = 0.036.
Fig. 2-7 conclude that evaporation from the liquid level in a porous medium significantly affects the liquid migration distribution in the layers and the liquid level in a porous medium distribution. Thus, one must consider evaporation when designing pumped-well drain canals in multilayer medium to improve lands and study liquid flows in layers.
The paper verified the accuracy of the numerical solution according to the balance equation. It appears that the error margin is limited by 1.3 %.
Figure 6. Flow curves reduced to an area unit with a cover and low-conductive layers at
qo = 0.0036.
4. Conclusions
The computational experiments have established that the dynamics of changes in the liquid level in a porous medium depends significantly on the evaporation parameter. With an increase in its value, the liquid level in a porous medium decreases proportionally over time.
An analysis of the obtained numerical calculations showed that the change in the liquid level in a porous medium depends significantly on the parameter n. As its value increases, the liquid level in the porous medium decreases at different values of the evaporation parameter.
The analyzing results of numerical calculations show that evaporation from the surface of a liquid in a porous medium significantly affects the distribution of overflows in the layers and the distribution of the liquid level in a porous medium. Therefore, when designing vertical drains to improve fluid selection in multilayer reservoirs and designing liquid flows in layers, it is essential to take into account evaporation.
Checking the accuracy of the numerical solution using the balance equation showed that the error does not exceed 1.3 %.
References
1. Tskhay A.A., Koshelev, K.B., Kim, N.Yu. Model vzaimodeystviya podzemnykh i poverkhnostnykh vod dlya sistemy podderzhki prinyatiya resheniy [Model of interaction between ground and surface waters for a decision support system]. Informatsionnyye sistemy v ekonomike, ekologii i obrazovanii. Barnaul: Izd-vo AltGTU, 2002. Pp. 39-41. (rus)
2. Anderson E.I. An analytical solution representing groundwater-surface water interaction. Water Resource. Res. 2003. 39(3). P. 1071. D0I:10.1029/2002WR001536
3. Ravshanov, N., Kodirov, K. Modeling of the process of free water filtration of groundwater based on the availability of gallery of wells. Problems of Computational and Applied Mathematics. 2016. 2. Pp. 33-46
4. Ravshanov, N. Abdullaev, Z., Khafizov, O. Modeling the Filtration of Groundwater in Multilayer Porous Media. Construction of Unique Buildings and Structures. 2020. 92. Article No. 9206/ISSN 2304-6295. DOI: 10.18720/CUBS.92.6
5. Kashevarov, A.A. Matematicheskoye modelirovaniye kachestva podzemnykh vod v zonakh vliyaniya beregovykh infiltratsionnykh vodozaborov [Mathematical modeling of groundwater quality in the zones of influence of coastal infiltration water intakes] [Online]. URL: https://www.rfbr.ru/rffi/ru/project_search/o_238194 (date of application: 03.12.2022).
6. Ravshanov, N., Daliev, S., Abdullaev, Z., Khafizov, O. Ground and confined underground waters and their salt content. International Conference on Information Science and Communications Technologies, ICISCT 2020. 2020. 9351467. DOI: 10.1109/ICISCT50599.2020.9351467
7. Ravshanov, N., Daliev, Sh., Abdullaev, Z., Khafizov, O. Ground and confined underground waters and their salt content. IOP Conference Series: Materials Science and Engineering. 2020. 896(1). 012047. DOI 10.1088/1757-899X/896/1/012047
8. Belov, K.V., Lisenkov, A.B., Ponomarev, A.D., Gorbatenko, N.S. Study of fluid filtration in a porous medium using physical and numerical modeling. Bulletin of the Tomsk Polytechnic University. Geo Assets Engineering. 2017. 328(8). Pp. 64-74
9. Buyskikh, A.A., Basistyy, V.A. Evaluation by the method of mathematical modelinof the exploitation resources of underground waters of the gold ore deposit. Tekhnologii tekhnosfernoy bezopasnosti. 2010. 6(34). (rus)
10. Gridnevskiy, A.V. Numerical simulation of the filtration process in the right-bank of the river Don for protection the buildings against rise of groundwater in the city of Rostov-on-Don. Geology and Geophysics of Russian South. 2019. 9(1). Pp. 150-163. (rus)
11. Shalanin, V.A., Patlai, K.I. Numerical simulation of the groundwater filtration process through a rectangular dam of homogeneous soil on a waterproof base. FEFU: School of engineering bulletin. 2019. 2/39. DOI: 10.24866/2227-6858/2019-2-13 (rus)
12. Bereslavskiy, E.N., Dudina, L.M. Problem of filtration in a rectangular web with a partially impenetrable vertical wall. Vestnik of
Saint Petersburg University. Mathematics. Mechanics. Astronomy. 2019. 6(2). Pp. 288-297.
DOI: 10.21638/11701/spbu01.2019.211 (rus)
13. Sheshukov, Ye.G., Kurtseva, K.P. Chislennoye issledovaniye krayevykh zadach nelineynoy filtratsii [Numerical investigation of boundary problems of nonlinear filtration]. Power engineering: research, equipment, technology. 2012. 9-10. Pp. 158-166. (rus)
14. Kurtseva, K.P., Lapin, A.V., Sheshukov, Ye.G. Resheniye setochnymi metodami zadachi filtratsii zhidkosti v plotine pri nelineynom zakone filtratsii [Solution by grid methods of the problem of fluid filtration in a dam with a nonlinear filtration law]. Izvestiya vuzov. Matematika. 1995. 2. Pp. 47-52. (rus)
15. Aniskin, A.N., Memarianfard, M.Ye. Chislennoye modelirovaniye anizotropnoy filtratsii v gruntovykh plotinakh [Numerical simulation of anisotropic filtration in ground dams]. Vestnik MGSU. 2009. 4. Pp. 219-224. (rus)
16. Aniskin, N.A., Memarianfard, M.Ye. Uchet anizotropii v filtratsionnykh raschetakh i raschetakh ustoychivosti otkosov gruntovykh plotin [Accounting for anisotropy in filtration calculations and calculations of slope stability of soil dams]. Vestnik MGSU. 2010. 4. Pp. 388-398. (rus)
17. Demidov V. N Chislennoye modelirovaniye protsessov formirovaniya dozhdevogo stoka [Numerical modeling of rain runoff formation processes]. Doctor thesis. Moscow, 2007. (rus)
18. Morozov P Ye Metody resheniya pryamykh i obratnykh zadach podzemnoy termogidrodinamiki [Methods for solving direct and inverse problems of underground thermohydrodynamics].Moscow, 2022. (rus)
19. Abutaliyev, F.B., Abutaliyev, E.B. Metody resheniya zadach podzemnoy gidromekhaniki na EVM [Methods for solving problems of underground hydromechanics on a computer]. Tashkent: Izd-vo “Fan”, 1968. 196 p.
20. Averyanov, S.F. Fil'tratsiya iz kanalov i yeye vliyaniye na rezhim gruntovykh vod [Filtration from canals and its impact on the groundwater regime]. Moscow: Kolos, 1956. 237 p.
21. Porubarinova-Kochinav, P.Ya. Teoriya dvizheniya gruntovykh vod [Theory of groundwater movement]. Moscow: Nauka, 1977. 664 p. (rus)
22. Khantush, M.S. Novoye o teorii peretikaniya [New about the theory of poking]. Voprosy gidrogeologicheskikh raschetov. Moscow: Mir, 1964. Pp. 43-60. (rus)
23. Ravshanov, N., Abdullaev, Z., Aminov, S., Khafizov, O. Numerical study of fluid filtration in three-layer interacting pressure porous formations. E3S Web of Conferences. 2021.264. 01018. DOI: 10.1051/e3sconf/202126401018
Information about authors:
Normahmad Ravshanov,
E-mail: ravshanxade-09@mail.ru
Zafar Abdullaev, PhD in Physics and Mathematics, ORCID: https://orcid.org/0000-0003-1351-7861 E-mail: abdullaevv. zafar@gmail. com
Evgeny Kotov,
E-mail: ekotov. cfd@gmail. com
Shodiya Turkmanova,
ORCID: https://orcid. org/0000-0003-0521-6882 E-mail: turkmanovashodiya@gmail. com
Received 04.12.2022. Approved after reviewing 26.12.2022. Accepted 26.02.2023.