Течение двухфазной жидкости в модельной пористой среде, образованной осесимметричными каналами переменного сечения

Abstract

The micro-inhomogeneities of the porous medium have a decisive influence on its filtrationcapacitive properties, such as porosity, absolute permeability, and relative phase permeability. This fact explains the relevance of a study of multiphase flows in the channels that constitute the porous media. There are several approaches to modeling porous media: periodic structures, filtration theory, percolation theory, and statistical hydrodynamics. The approach based on periodic structures, has several advantages: to determine the reservoir properties of the simulated medium (knowing the geometry of one characteristic pore), it is enough to describe the fluid flow inside, based on analytical estimates or from the results of numerical solution of the system of hydrodynamic equations. This article considers the features of the flow of a two-phase fluid in axisymmetric channels of variable cross section with periodic boundary conditions. The saturation of each phase depends on the initial distribution of these phases in the channel and does not depend on time. Furthermore, the flow of a two-phase incompressible linearly viscous fluid is calculated numerically based on the Navier-Stokes system of equations. At each time point, the reduced pulse of both phases was determined. The time-averaged value of the reduced pulse in the quasi-periodic flow regime is used to calculate the relative phase permeabilities in the dependence of the initial ratio of the phase saturations. All calculations are done in open packages. The channel geometry and computational grid are built in the Salome package, the flow is calculated in the OpenFOAM computation package, the visualization of the calculation results is performed in the paraView application.

Микронеоднородности пористой среды имеют определяющее влияние на ее фильтрационно-емкостные свойства, такие как пористость, абсолютная проницаемость, относительная фазовая проницаемость. Поэтому исследование многофазных течений в каналах, составляющих пористые среды, является актуальным. Существует несколько подходов к моделированию пористых сред: периодические структуры, теория фильтрации, теория перколяции, статистическая гидродинамика. Подход на основе периодических структур обладает рядом преимуществ: для определения фильтрационно-емкостных свойств моделируемой среды достаточно, зная геометрию одной характерной поры, описать течение флюида в ней на основе аналитических оценок либо по результатам численного решения системы уравнений гидродинамики. В работе рассмотрены особенности течения двухфазной жидкости в осесимметричных каналах переменного сечения с периодическими граничными условиями. Насыщенность каждой из фаз зависит от начального распределения этих фаз в канале и не зависит от времени. Далее численно рассчитано течение двухфазной несжимаемой линейно-вязкой жидкости на основе системы уравнений Навье – Стокса. В каждый момент времени определялся приведенный импульс обеих фаз. Осредненное по времени значение приведенного импульса на квазипериодическом режиме течения использовано для расчета относительных фазовых проницаемостей в зависимости от начального соотношения фазовых насыщенностей. Все вычисления проделаны в открытом программном обеспечении: геометрия канала и расчетная сетка построены в пакете Salome, расчет течения – в OpenFOAM, визуализация результатов расчетов – в paraView.