Все выпуски

В последние десятилетия важное значение приобретает разработка методов повышения эффективности извлечения углеводородов в месторождениях с нетрадиционными запасами, содержащими в больших количествах газовый конденсат. Это делает актуальным развитие методов математического моделирования, реалистично описывающих процессы фильтрации газоконденсатной смеси в пористой среде.

В данной работе рассматривается математическая модель, описывающая динамику изменения давления, скорости и концентрации компонент двухкомпонентной двухфазовой смеси, поступающей в лабораторную модель пласта, заполненную пористым веществом с известными физико-химическими свойствами. Математическая модель описывается системой нелинейных пространственно-одномерных дифференциальных уравнений в частных производных с соответствующими начальными и граничными условиями. Лабораторные эксперименты показывают, что в течение конечного времени система стабилизируется, что дает основание перейти к стационарной постановке задачи.

Численное решение сформулированной системы обыкновенных дифференциальных уравнений реализовано в среде Maple на основе метода Рунге–Кутты с автоматическим выбором шага. Показано, что полученные на этой основе физические параметры двухкомпонентной газоконденсатной смеси из метана и н-бутана, характеризующие моделируемую систему в режиме стабилизации, хорошо согласуются с имеющимися экспериментальными данными.

Это подтверждает реалистичность выбранного подхода и обоснованность его дальнейшего развития и применения для компьютерного моделирования неравновесных физических процессов в газоконденсатных смесях в пористой среде с целью выработки в перспективе практических рекомендаций по увеличению извлекаемости углеводородного газоконденсата из природных месторождений. В работе представлена математическая постановка системы нелинейных уравнений в частных производных и соответствующей стационарной задачи, описан метод численного исследования, обсуждаются полученные численные результаты в сравнении с экспериментальными данными.

In the last decades, the development of methods for increasing the efficiency of hydrocarbon extraction in fields with unconventional reserves containing large amounts of gas condensate is of great importance. This makes important the development of methods of mathematical modeling that realistically describe physical processes in a gas-condensate mixture in a porous medium.

In the paper, a mathematical model which describes the dynamics of the pressure, velocity and concentration of the components of a two-component two-phase mixture entering a laboratory model of plast filled with a porous substance with known physicochemical properties is considered. The mathematical model is based on a system of nonlinear spatially one-dimensional partial differential equations with the corresponding initial and boundary conditions. Laboratory experiments show that during a finite time the system stabilizes, what gives a basis to proceed to the stationary formulation of the problem.

The numerical solution of the formulated system of ordinary differential equations is realized in the Maple environment on the basis of the Runge–Kutta procedure. It is shown that the physical parameters of the gascondensate mixture, which characterize the modeled system in the stabilization regime, obtained on this basis, are in good agreement with the available experimental data. This confirms the correctness of the chosen approach and the validity of its further application and development for computer modeling of physical processes in gas-condensate mixtures in a porous medium. The paper presents a mathematical formulation of the system of partial differential equations and of respective system stationary equations, describes the numerical approach, and discusses the numerical results obtained in comparison with experimental data.

В химической технологии для получения очищенного конечного продукта часто используется процесс десублимации. Для этого используются охлаждаемые жидким азотом или холодным воздухом емкости. Смесь газов протекает внутри емкости и охлаждается до температуры конденсации или десублимации некоторых компонентов газовой смеси. Конденсированные компоненты оседают на стенках емкости. В статье представлена математическая модель для расчета охлаждения емкостей для десублимации паров охлажденным воздухом. Математическая модель основана на уравнениях газовой динамики и описывает течение охлажденного воздуха в трубопроводе и воздушном теплообменнике с учетом теплообмена и трения. Теплота фазового перехода учитывается в граничном условии для уравнения теплопроводности путем задания потока тепла. Перенос тепла в теплоизолированных стенках трубопровода и в стенках емкости описывается нестационарными уравнениями теплопроводности. Решение системы уравнений проводится численно. Уравнения газовой динамики решаются методом С. К. Годунова. Уравнения теплопроводности решаются по неявной разностной схеме. В статье приведены результаты расчетов охлаждения двух последовательно установленных емкостей. Начальная температура емкостей равна 298 К. Холодный воздух течет по трубопроводу, через теплообменник первой емкости, затем по трубопроводу в теплообменник второй емкости. За 20 минут емкости остывают до рабочей температуры. Температура стенок емкостей отличается от температуры воздуха на величину не более чем 1 градус. Поток охлажденного воздуха позволяет поддерживать изотермичность стенок емкости в процессе десублимации компонентов из газовой смеси. Приведены результаты аналитической оценки времени охлаждения емкости и разности температуры между стенками емкости и воздухом в режиме десублимации паров. Аналитическая оценка основана на определении времени термической релаксации температуры стенок емкости. Результаты аналитических оценок удовлетворительно совпадают с результатами расчетов по представленной модели. Предложенный подход позволяет проводить расчет охлаждения емкостей потоком холодного воздуха, подаваемого по трубопроводной системе.

For the production of purified final product in chemical engineering used the process of desublimation. For this purpose, the tank is cooled by liquid nitrogen or cold air. The mixture of gases flows inside the tank and is cooled to the condensation or desublimation temperature some components of the gas mixture. The condensed components are deposited on the walls of the tank. The article presents a mathematical model to calculate the cooling air tanks for desublimation of vapours. A mathematical model based on equations of gas dynamics and describes the movement of cooled air in the duct and the heat exchanger with heat exchange and friction. The heat of the phase transition is taken into account in the boundary condition for the heat equation by setting the heat flux. Heat transfer in the walls of the pipe and in the tank wall is described by the nonstationary heat conduction equations. The solution of the system of equations is carried out numerically. The equations of gas dynamics are solved by the method of S. K. Godunov. The heat equation are solved by an implicit finite difference scheme. The article presents the results of calculations of the cooling of two successively installed tanks. The initial temperature of the tanks is equal to 298 K. Cold air flows through the tubing, through the heat exchanger of the first tank, then through conduit to the heat exchanger second tank. During the 20 minutes of tank cool down to operating temperature. The temperature of the walls of the tanks differs from the air temperature not more than 1 degree. The flow of cooling air allows to maintain constant temperature of the walls of the tank in the process of desublimation components from a gas mixture. The results of analytical evaluation of the time of cooling tank and temperature difference between the tank walls and air with the vapor desublimation. Analytical assessment is based on determining the time of heat relaxation temperature of the tank walls. The results of evaluations are satisfactorily coincide with the results of calculations by the present model. The proposed approach allows calculating the cooling tanks with a flow of cold air supplied via the pipeline system.

Одна из трудностей разработки газоконденсатных месторождений обусловлена тем, что часть углеводородов газоносного слоя присутствует в немв виде конденсата, который застревает в порах пласта и извлечению не подлежит. В этой связи активно ведутся исследования, направленные на повышение извлекаемости углеводородов в подобных месторождениях. В том числе значительное количество публикаций посвящено развитию методов математического моделирования прохождения многокомпонентных газоконденсатных смесей через пористую среду в различных условиях.

В настоящей работе в рамках классического подхода, основанного на законе Дарси и законе неразрывности потоков, сформулирована математическая постановка начально-граничной задачи для системы нелинейных дифференциальных уравнений, описывающая прохождение многокомпонентной газоконденсатной смеси через пористую среду в режиме истощения. Разработанная обобщенная вычислительная схема на основе конечно-разностной аппроксимации и метода Рунге – Кутты четвертого порядка может использоваться для расчетов как в пространственно одномерном случае, соответствующемусловиям лабораторного эксперимента, так и в двумерном случае, когда речь идет о моделировании плоского газоносного пласта с круговой симметрией.

Численное решение упомянутой системы уравнений реализовано на основе комбинированного использования C++ и Maple с применением технологии параллельного программирования MPI для ускорения вычислений. Расчеты выполнены на кластере HybriLIT Многофункционального информационно-вычислительного комплекса Лаборатории информационных технологий Объединенного института ядерных исследований.

Численные результаты сопоставлены с данными о динамике выхода девятикомпонентной углеводородной смеси в зависимости от давления, полученными на лабораторной установке (ВНИИГАЗ, Ухта). Расчеты проводились для двух типов пористого наполнителя в лабораторной модели пласта: терригенного (при 25 ^◦С) и карбонатного (при 60 ^◦С). Показано, что используемый подход обеспечивает согласие полученных численных результатов с экспериментальными данными. Путем подгонки к экспериментальным данным по истощению лабораторной модели пласта получены значения параметров, определяющих коэффициент межфазного перехода для моделируемой системы. С использованием тех же параметров было проведено компьютерное моделирование истощения тонкого газоносного слоя в приближении круговой симметрии.

One of problems in developing the gas condensate fields lies on the fact that the condensed hydrocarbons in the gas-bearing layer can get stuck in the pores of the formation and hence cannot be extracted. In this regard, research is underway to increase the recoverability of hydrocarbons in such fields. This research includes a wide range of studies on mathematical simulations of the passage of gas condensate mixtures through a porous medium under various conditions.

In the present work, within the classical approach based on the Darcy law and the law of continuity of flows, we formulate an initial-boundary value problem for a system of nonlinear differential equations that describes a depletion of a multicomponent gas-condensate mixture in porous reservoir. A computational scheme is developed on the basis of the finite-difference approximation and the fourth order Runge .Kutta method. The scheme can be used for simulations both in the spatially one-dimensional case, corresponding to the conditions of the laboratory experiment, and in the two-dimensional case, when it comes to modeling a flat gas-bearing formation with circular symmetry.

The computer implementation is based on the combination of C++ and Maple tools, using the MPI parallel programming technique to speed up the calculations. The calculations were performed on the HybriLIT cluster of the Multifunctional Information and Computing Complex of the Laboratory of Information Technologies of the Joint Institute for Nuclear Research.

Numerical results are compared with the experimental data on the pressure dependence of output of a ninecomponent hydrocarbon mixture obtained at a laboratory facility (VNIIGAZ, Ukhta). The calculations were performed for two types of porous filler in the laboratory model of the formation: terrigenous filler at 25 .„R and carbonate one at 60 .„R. It is shown that the approach developed ensures an agreement of the numerical results with experimental data. By fitting of numerical results to experimental data on the depletion of the laboratory reservoir, we obtained the values of the parameters that determine the inter-phase transition coefficient for the simulated system. Using the same parameters, a computer simulation of the depletion of a thin gas-bearing layer in the circular symmetry approximation was carried out.