Mathematical modeling of two-phase compressible fluid filtration based on modified adaptive method of minimum amendments

The work objective is to build and investigate the modified adaptive method of minimum amendments (MAMMA) which is destined for the numerical simulation of the two-phase compressible fluid filtration in porous media. This approach allows overcoming the known use limitations of other methods of the finite-difference equations solution, such as: crucial differential pressures acting on the oil-and-water bearing formation; and the compressibility of the medium at the considerable gas content in the oil phase. An approximation method - an explicit one for defining the function of water saturation, and an implicit one for the pressure function computation - is selected as the research basis. When setting the initial boundary value problem and its sampling, the process of the two-phase compressible fluid filtration in the space-dimensional domain with the lateral area bounded below by the subface of stratum, and above - by the bed top, is considered. A two-layer iterative method of the variational type - a modified method of minimal amendments adapted for solving finite-difference equations of the two-phase compressible fluid with a non-selfadjoint operator under the most general assumptions on the properties of the grid-problem operator is built. It is shown that a MAMMA has the asymptotic convergence rate characteristic of the “classical” alternate triangular method that does not use the Chebyshev acceleration technique and can be applied to the problems with a self-adjoint operator. Numerical experiments have confirmed the high efficiency of MAMMA. It is established that to achieve the specified accuracy, the number of iterations at the MAMMA reduces to 3-20 times as compared to the method of Seidel and the overrelaxation method.

задачи фильтрации двухфазных сжимаемых жидкостей, сеточные уравнения, несамосопряженный оператор, адаптивный метод минимальных поправок, two-phase compressible fluid filtration problems, finite-difference equations, non-selfadjoint operator, adaptive method of minimum amendments

1. Sukhinov, А.I., Grigoryan, L.A., Sukhinov, A.A. Usovershenstvovannyy poperemenno-treugol'nyy metod chislennogo resheniya prostranstvenno-trekhmernykh zadach fil'tratsii dvukhfaznoy neszhimaemoy zhidkosti. [Advanced alternate triangular method of numerical solution of spatial three-dimensional problems of two-phase filtration of incompressible fluid.] Vestnik of DSTU, 2016, vol. 16, no. 1, pp. 5–18 (in Russian).

2. Konovalov, А.N. Zadachi fil'tratsii mnogofaznoy neszhimaemoy zhidkosti. [Problems of multiphase filtration of incompressible fluid.] Novosibirsk: Nauka, 1988, 166 p. (in Russian).

3. Vabishchevich, P.N. Yavno-neyavnye vychislitel'nye algoritmy dlya zadach mnogofaznoy fil'tratsii. [Explicitimplicit numerical algorithms for multi-phase filtration problems.] Mathematical Models and Computer Simulations, 2010, vol. 22, no. 4, pp. 118–128 (in Russian).

4. Dushin, V.R., et al. Mathematical modeling of flows in porous media. WSEAS Transactions on fluid mechanics, 2014, vol. 9, pp. 130–166.

5. Sukhinov, A.I., Chistyakov, A.E. Adaptive modified alternating triangular iterative method for solving grid equations with a non-self-adjoint operator. Mathematical Models and Computer Simulations, 2012, vol. 4, iss. 4, pp. 398–409.

6. Sukhinov, A.I., Shishenya, A.V. Increasing efficiency of alternating triangular method based on improved spectral estimates. Mathematical Models and Computer Simulations, 2013, vol. 5, iss. 3, pp. 257–265.

7. Sukhinov, A.I., Chistyakov, A.E., Shishenya, A.V. Modifikatsiya metoda minimal'nykh popravok dlya resheniya setochnykh uravneniy s nesamosopryazhennym operatorom. [Modification of minimal residuals iterative method for solving grid equations with non-selfadjoint operators.] Izvestiya SFedU. Engineering Sciences, 2013, no. 4, pp. 194–202 (in Russian).