NUMERICAL SOLUTION TO THE PROBLEM OF TWO–PHASE FILTRATION WITH HETEROGENEOUS COEFFICIENTS BY THE FINITE ELEMENT METHOD
Аннотация
Рассматривается процесс фильтрации двухфазной жидкости в пористой неоднородной среде. Данный процесс описывается связанной системой уравнений для насыщенности, скорости фильтрации и порового давления. Рассмотрены математические модели с учетом и без учета капиллярных сил, при наличии которого для насыщенности имеем нестационарное уравнение конвекции-диффузии. Поскольку данный процесс характеризуется существенным преобладанием конвективного слагаемого в уравнении для насыщенности, используются противопотоковые аппроксимации посредством добавления неоднородной искусственной диффузии. Скорость и давление аппроксимируются с использованием смешанного метода конечных элементов. Представлены результаты численных расчетов для двумерного случая с сильно неоднородными коэффициентами проницаемости пористой среды. Рассмотрены несколько случаев, связанных с линейными и нелинейными коэффициентами относительной проницаемости флюида и наличием капиллярных сил. We consider the process of filtration of a two-phase fluid in a porous, heterogeneous medium. This process is described by a coupled system of equations for saturation, filtration rate, and pore pressure. We consider mathematical models with and without capillary forces, in the presence of which, for saturation, we have a nonstationary convection-diffusion equation. Since this process is characterized by a significant predominance of the convective term in the equation for saturation, countercurrent approximations are used by adding non-uniform artificial diffusion. Speed and pressure are approximated using a mixed finite element method. The results of numerical calculations for a two-dimensional case with strongly heterogeneous permeability coefficients of a porous medium are presented. Several cases of relative fluid permeability associated with linear and nonlinear coefficients and the presence of capillary forces are considered.