Численное решение задачи об образовании гидрата двуокиси углерода в пористом пласте, изначально насыщенном гидратом метана

Abstract

Актуальность работы обусловлена теоретическим исследованием процессов образования газовых гидратов в пористых средах применительно к проблемам добычи, транспортировки, консервации и хранения углеводородного сырья. Цель исследования: на основе математической модели рассмотреть в осесимметричной постановке задачу о замещении гидрата метана, изначально содержащегося в пористом пласте, гидратом углекислого газа и выявить закономерности динамики образования газовых гидратов в пористой среде конечной протяженности при нагнетании газа через его левую границу, когда другая внешняя граница является открытой для потока газа. Методы исследования: методы механики многофазных систем совместно с численными методами решения обобщенной задачи Стефана. Результаты. Представлена математическая модель, и проведено численное исследование закачки диоксида углерода в пористый пласт, изначально насыщенный метаном и его гидратом, в осесимметричной постановке. Для решения поставленной задачи, заключающейся в нахождении распределения полей давления и температуры, а также в определении положения границы фазового перехода, был применен численный метод ловли фронтов в узел пространственной сетки. Установлено, что процесс замещения газогидрата сопровождается незначительным понижением давления и повышением температуры среды. Показано, что температура на границе фазового перехода зависит от давления закачиваемого газа. При высоких значениях давления нагнетаемого газа она стремится к своему предельному значению, зависящему от исходной гидратонасыщенности пласта. Исследованы зависимости времени полного вытеснения метана из пласта в зависимости от его проницаемости и давления нагнетаемого газа.

The relevance of the discussed issue is caused by theoretical investigation of gas hydrate generation in porous medium with regard to production, transportation and conservation of hydrocarbon raw material. The main aim of the study is, to consider in axisymmetric formulation the replacement of methane hydrate, initially contained in a porous formation, by carbon dioxide hydrate, on the basis of the mathematical model, and to reveal the regularities of dynamics of gas hydrate generation in porous media of finite extent while injecting gas through its left boundary, with the other boundary being open to gas flow. The methods used in the study: multiphase system mechanics together with computational method of Stefenґs generalized problem. Results of the study. The paper introduces the mathematical model and numerical calculation of carbon dioxide injection into porous formation, initially saturated with methane and its hydrate in the axisymmetric formulation. The authors have used the numerical method of tracking the front node in the spatial grid to solve the problem of finding pressure and temperature distribution, and the boundary of phase transmission. The acquired numerical solutions were compared with analytical self-similar solutions. It was ascertained that gas hydrate displacement is followed by little pressure drop and reservoir temperature increase. It was shown that the temperature at the phase border depends on injected gas pressure. If the pressure is high it tends toward its cutoff value which depends on initial gas hydrate saturation in the reservoir. The authors studied as well the relation of methane total displacement from the reservoir as a function of its permeability and injected gas pressure.