HYDRODYNAMIC MODEL OF UNSTEADY MOVEMENT OF OIL AND GAS IN THE RESERVOIR-WELL SYSTEM
Abstract (English):
A hydrodynamic model of unsteady movement of oil and gas in the reservoir-well system is constructed and solutions to the system of differential equations are given. The obtained analytical formula makes it possible to determine the dynamics of the pressure at the bottom of the well and the productivity of the formation, depending on the parameters of the system. Numerical calculations have been carried out for practical values of the reservoir-well systemparameters.

Keywords:
reservoir-wells, Laplace transform, fluid movement, gas movement, continuity equation

Introduction

Gas lifting of liquid from a well is widely used in oil production. The increase in the efficiency of the hoist is of great practical and scientific importance. Despite the fact that a lot of works [1-6] are devoted to this problem, the issue of increasing the efficiency of gas lift taking into account the movement of fluid in the reservoir-well system is still poorly understood. And the movement of fluid during oil production occurs in the reservoir-well system. Therefore, modeling and studying the movement of gas lift fluid in the reservoir-well system is of great practical and scientific importance.

Statement and solution of the problem

Consider plane-radial filtration of a homogeneous fluid in a circular homogeneous reservoir.

The differential equation of the piezoconductivity of a plane-radial fluid flow has the form [7,8]

(1)

where ;.

Initial and boundary conditions

(2)

(3)

(4)

The solution of equation (1) with the initial (2) and boundary conditions (3) and (4) has the form [7,9]

(5)

wherepressure on the borehole wall, at the initial moment of time,  differential pressure to which the initial differential decreases,  - pressure on the borehole wall at the bottom of the borehole after pressure change

- roots of the transcendental equation

(6)

The fluid flow rate at the moment t through the lateral surface of the well with radius is determined by the formula

(7)

Then from expression (7) taking into account formula (5), we obtain

(8)

where                         ,

Gas movement in annular space.

Now consider the motion of gas in the annular space. The equations of gas motion in annular space and continuity are described by the equations of I.A. Charny [10,11]

,                                (9)

Differentiating the first equation of expression (9) in the , coordinate, and the second in time  , and subtract one from the other and get:

(10)

Having placed the origin of the coordinate axis at the wellhead and directed it downward, for the initial and boundary conditions we will have:

;                                                                      (11)

(12)
(13)

(14)

where- gas injection pressure, - bottomhole pressure, -speed of sound in gas, - drag coefficient, t time x-coordinate,- mass flow rate of gas in the flow area of the annular space, - gas density at a given pressure, P-pressure in any cross-section of the annular space.

The solution of equation (10), taking into account the boundary conditions (13) and (14), will be sought in the form [12,13]

(15)

Substituting expression (15) into equation (10), multiplying both sides of the resulting expression by  and integrating it fromto , we get the equation:

(16)

Applying the Laplace transform [14,15] from expression (16) we obtain

(17)

and are determined from expression (15) and have the form

respectively, the initial value of the gas supply pressureand bottomhole pressure.

Applying the Laplace transform from expression (15), we obtain

(18)

Substituting expression (17) into formula (18), and the resulting expression into the first equation of system (9) after the Laplace transform, we obtain:

(19)

where- initial mass flow rate of gas.

The movement of the gas-liquid mixture in the riser pipe.

The gas supplied through the annular space and the liquid coming from the formation are mixed in the riser string and rises through it. The density of the mixture in this case can be determined by the formula [6]

where- respectively, the density of the mixture, oil and gas at a given pressure, -mass fraction of oil in gas.

In the first approximation, the mixture is assumed to be homogeneous and the interactions between the liquid and gas bubbles are neglected.

Then the equation of motion of the mixture will have the same form as equation (10). By placing the origin of the х1 coordinate axis at the lower end of the pipe string and directing it upwards for the initial and boundary conditions, we will have

(20)

(21)
(22)

(23)

The solution of equation (10), taking into account the boundary conditions (22) and (23), will be sought in the form [9]

(24)

Substituting expression (24) into equation (10), multiplying both sides of the resulting expression byand integrating it from to, we will have:

(25)

Applying the Laplace transform from expression (24), we obtain

(26)

and are determined from expression (24) taking into account the initial conditions (20) and (21) and have the form

where

Applying the Laplace transform from expression (24), we obtain

(27)

Substituting formula (27) into the first equation, system (9) after the Laplace transform, we obtain

(28)

where -the initial value of the mass flow rate of the mixture.

Continuity conditions

(29)

where  and-respectively, the area of the flow section of the pipe string and the annular section. Applying the Laplace transform from expression (29), we obtain

(30)

Substituting expressions (8), (19), and (28) into the continuity condition (30), taking into account only one term of the series in the first approximation, we obtain

(31)

where,

and are determined from the following formulas

,

where -gas density at atmospheric pressure; -acceleration of gravity, - atmospheric pressure; - oil inflow from the reservoir per unit of time at the initial moment of time.

From expression (28), taking into account expressions (26) and (27), one can determine  .

Conclusion

A hydrodynamic model of unsteady movement of oil and gas in the reservoir-skavazhin system has been built and analytical expressions have been obtained that allow, depending on the parameters of the reservoir-well system, to determine the bottomhole pressure and well productivity.

1. Barashkin R.L., Samarin I.V. "Hydrodynamic modeling of the production process of gas-lift wells" // Bulletin of the Tomsk Polytechnic University. 2006. V. 309. № 6

2. Duc Hien Bui, Rais Yanfurovich Nugaev, Rim Khalitovich Khazipov "On the issue of optimizing the operating modes of gas-lift wells" // Scientific and technical journal "Problems of collection, preparation and transportation of oil and oil products." 2012 №1

3. Kashif Rashid, William Bailey, Benoit Couet "A Survey of Methods for Gas-Lift Optimization"// Modelling and Simulation in Engineering, september, 2012.

4. Aliev F.A., Ilyasov M.Kh, Dzhamalbekov M.A. "Modeling the operation of a gas-lift well" // "OilGasScientificResearchProject" Institute, SOCAR, 2008.

5. Aliev F.A., Aliev N.A., Guliev A.P., Tagiev R.M., Dzhamalbekov M.A. "Method for solving a boundary value problem for a system of hyperbolic equations describing motion in a gas-lift process" // Applied Mathematics and Mechanics. 2018.V 82 issue 4.P. 512-518

6. Abbasov E.M. Determination of fluid accumulation time in periodic gas-lift wells. Engineering Physics Journal 2013, V86, №2, P.310-317

7. Shchelkachev V.N. Foundations and applications of the theory of transient filtration. V 24.M.: Oil and gas, 1995, 585 P.

8. Charny A.I. Underground fluid dynamics. Moscow: Gostommehizdat, 1963.

9. Agaeva N.A. "Hydrodynamics of fluid movement in the conjugate system reservoir-well-pipeline" // Engineering Physics Journal, Volume 94, № 1, 2021, P.62-71

10. Charny A.I. Unsteady fluid movement in pipes. Moscow: Bowels, 1975.

11. Huseynzade M.A., Druchina L.I., Petrova O.N., Stepanova M.F. Hydrodynamic processes in complex pipeline systems. Moscow: Bowels, 1991.

12. Abbasov E.M., Kengerli T.S., Abdullaeva N.R. Modeling the process of filtration of a gas-liquid mixture in a conjugated system of plastisin. Engineering Physics Journal.2020, V.92, № 5, P.1122-1196

13. Abbasov E.M., Kerimova Sh.A., Agaeva N.A. Integral modeling of the pressure recovery process. Engineering Physics Journal, 2019, V.92, №6, P.2475-2481

14. Aramonovich I.G., Lunts G.L., Elsgolts E.E. Complex variable functions. Operational calculus. Stability theory. M.: Science, 1968, 390P.

15. Dec. D Guide to the practical application of the Laplace transform. M.: Science,1965.