532.546 The efficient method for numerical solution of dual porosity reservoir model
Maykov D. N. (Udmurt Federal Research Center of the Ural Branch of the Russian Academy of Sciences), Makarov S. S. (Udmurt Federal Research Center of the Ural Branch of the Russian Academy of Sciences)
MODEL, DUAL POROSITY RESERVOIR NUMERICAL SOLUTION, MULTILATERAL WELL, DIFFUSIVITY EQUATION, PARAMETRIC ANALYSIS
doi: 10.18698/2309-3684-2024-3-317
This paper presents a method for accelerating the numerical solution of the diffusivity equation based on the Warren-Root model. A general differential equation system describing the filtration model from the matrix to the fracture is written through the complex parameters of the transmissivity ratio, the storability ratio, and the volumetric average permeability of the fracture system. The proposed method for accelerating the numerical solution of a differential equation system describing a double-porosity reservoir model is based on converting the traditional form of a finite-difference approximation of the system for two differential equations into one equation. A stable implicit difference scheme is used to obtain a finite-difference approximation of the parameters. Boundary conditions of the first and second kind are considered: a constant pressure boundary and an impermeable boundary. The results of the test calculations using the proposed method are compared with the analytical solution. The pressure change in the well was compared, calculated by numerical and analytical methods. The pressure in the well was calculated using the Peaceman method with the effective radius for the Voronoi grid cell. A numerical analysis of the parameters of a multilateral well in a double porosity formation model is carried out. A two-dimensional Cartesian unstructured irregular Voronoi grid was used as the calculated grid. The numerical calculations of the matrix equations were carried out by three different methods: the biconju-gate gradient stabilized method with ILU(0) preconditioning, the Gauss-Seidel method with relaxation, and the Newton method. It is shown that the implementation of a differential equation system according to the proposed method significantly reduces the complexity of the numerical solution and reduces the calculation time of the filtration process modeling and pressure transient analysis interpretation.
[1] Walker A.C. Estimating reservoir pressure using the principle of superposition. SPE, 1968. DOI: 10.2118/2324-MS
[2] Cinco-Ley H., Samaniego F. Use and misuse of the superposition time function in well test analysis. SPE-19817,1989. DOI: 10.2118/19817-MS
[3] Maykov D.N., Isupov S.V., Makarov S.S., Anikanov A.S. Metod uskoreniya rascheta davleniya pri izmenyayushchihsya debitah po istorii ekspluatacii skvazhiny [The efficient method for pressure calculation at variable rate]. Neftyanoe hozyajstvo [Oil Industry], № 9, 2021, pp. 105–107. DOI: 10.24887/0028-2448-2021-9-105-107
[4] Maykov D.N., Makarov S.S. Numerical investigation of optimization algorithms for the hydrodynamic model adaptation based on the well test results. Mathematical Models and Computer Simulations, 2022, vol. 34, № 9, pp. 71-82.
[5] Gordonov A.O., Laptev I.V., Sidorenko N.Yu. Mathematical modeling of laminar and turbulent filtration processes of liquid in compressible medium in porous meshmaterials. Mathematical modeling and computational methods, 2023, № 2, pp. 67-89. DOI: 10.18698/2309-3684-2023-2-6789
[6] Kholmurodov A.E., Dilmuradov N. Mathematical modeling of one-dimensional non-linear motion in a fluid-saturated porous medium. Mathematical modeling and computational methods, 2018, № 1, pp. 3-15. DOI: 10.18698/2309-3684-2018-1-315
[7] Blonsky A.V., Mitrushkin D.A., Savenkov E.B. Discrete fracture network modelling: physical and mathematical model. Keldysh Institute PREPRINTS, 2017, № 65, p. 28. DOI 10.20948/prepr-2017-65
[8] Maaref S., Moosavi R., Riazi M. Comparison of numerical methods to solve Warren-Root PDE of naturally fractured reservoirs. The 1st National Conference on Oil and Gas Fields Development, 2015.
[9] Istrafilov M.Y., Morozkin N.N. Solution of the filtration problem in the fissured formation using warren-root model. Vestnik Bashkirskogo universiteta, 2017, vol. 22, no. 1, pp. 15-19.
[10] Afanaskin I.V., Volpin S.G., Lomakina O.V., ShteynbergYu. M. Carbonate reservoirs crosswell survey interpretation by a two-rate test using numerical models. Software & Systems, 2018, no. 3, pp. 500-506.
[11] Vasilyev V.I., Vasilyeva M.V., Grigorev A.V., Prokopiev G.A. Mathematical modeling of the two-phase fluid flow in inhomogeneous fractured porous media using the double porosity model and finite element method. Lobachevskii Journal of Mathematics, 2018, vol. 160, no. 1, pp. 165-182.
[12] Vabishchevich P.N., Grigoriev A.V. Numerical modeling of a fluid flow in anisotropic fractured porous media. Numerical Analysis and Applications, 2016, vol. 19, no. 1, pp. 61-74. DOI: 10.15372/SJNM20160105.
[13] Vabishchevich P.N., Grigoriev A.V. CHislennoe modelirovanie fil'tracii na osnove modeli dvojnoj poristosti [Numerical simulation of filtration based on the dual porosity model]. Superkomp'yuternye tekhnologii matematicheskogo modelirovaniya: trudy vtoroj mezhdunarodnoj konferencii [Second international conference Supercomputer technologies of mathematical modelling], 2013, pp. 50-61.
[14] Grigoriev A.V. Chislennoe modelirovanie filtracii v treshchinovato-poristoj srede [Numerical modeling of filtration in a fractured porous medium]. Matematicheskie zametki YAGU [Mathematical Notes of NEFU], 2013, vol. 20, no. 2, pp. 237-245.
[15] Warren J.E., Root P.J. The behavior of naturally fractured reservoirs. SPE Journal, 1963, pp. 244–255. DOI: 10.2118/426-PA
[16] Stewart G. Well test design and analysis. PennWell Corp., 2011, pp. 549–603.
[17] Abdus Satter, Ghulam M. Iqbal. Reservoir Engineering: The Fundamentals, Simulation, and Management of Conventional and Unconventional Recoveries. Texas, Gulf Professional Publishing, 2016, 486 p.
[18] Fortune S. A sweepline algorithm for Voronoi Diagrams. Algorithmica, 1987, vol. 2, pp. 153-174.
[19] Zhang L., Wenshu Zha, Lu D. Generation and application of adaptive PEBI grid for numerical well testing (NWT). Proceedings 2013 International Conference on Mechatronic Sciences, 2013. DOI: http://dx.doi.org/10.1109/MEC.2013.6885542
[20] Maykov D.N., Makarov S.S. Parametricheskij analiz modeli mnogozabojnoj skvazhiny v plaste s dvojnoj poristost'yu [The multilateral well model with complete vertical opening of naturally fractured reservoir]. Vestnik TyumGU. Fiziko-matematicheskoe modelirovanie. Neft', gaz, energetika [Tyumen State University Herald. Physical and Mathematical Modeling. Oil, Gas, Energy], 2023, vol. 9, no. 3(35), pp. 100-116.
[21] Maykov D.N., Borkhovich S.Yu. Analiticheskaya model' mnogozabojnoj skvazhiny s polnym vertikal'nym vskrytiem plasta [Analytical model of a multilateral well with full vertical penetration of the formation]. Neft'. Gaz. Novacii [Oil. Gas. Innovations], no. 11 (240), 2020, pp. 61–65.
[22] Peaceman D.W. Interpretation of well-block pressures in numerical reservoir simulation with nonsquare grid blocks and anisotropic permeability. Society of Petroleum Engineers Journal, 1983, 23(03), pp. 531–543.
[23] Syihab Z. Simulation on Discrete Fracture Network Using Flexible Voronoi Grid System. Petroleum Engineering, Texas A&M University, 2009, 194 p.
[24] Saad Y. Iterative methods for sparse linear systems. SIAM, 2003, 184 p.
[25] Glazyrina L.L., Karchevskij M.M. Vvedenie v chislennye metody: uchebnoe posobie [Introduction to numerical methods: textbook]. Kazan', Kazanskij universitet, 2017, 122 p.
[26] Bahvalov N.S., Zhidkov N.P., Kobelkov G.M. CHislennye metody [Numerical methods]. Moscow: BINOM. Lab. znanij [BINOM. Lab. knowledge], 2007, 636 c.
[27] Bourdet D. et al. A new set of type curves simplifies well test analysis. World Oil, 1983, pp. 95–106.
Майков Д.Н., Макаров С.С. Метод ускорения численного решения дифференциального уравнения пьезопроводности модели пласта с двойной пористостью. Математическое моделирование и численные методы, 2024, № 3, с. 3–17.
Download article
Количество скачиваний: 455