Математическая модель отбора газа из газогидратного пласта с учетом образования льда

Тип работы:
Реферат
Предмет:
Механика


Узнать стоимость

Детальная информация о работе

Выдержка из работы

УдК 532. 546:536. 421
математическая модель отбора газа из газогидратного пласта с учетом образования льда
Хасанов М.к., доровская М.с.
Стерлитамакский филиал Башкирского государственного университета, Стерлитамак, e-mail: hasanovmk@mail. ru
Рассмотрены особенности разложения газогидратов при отборе газа из пористой среды, в исходном состоянии заполненной газом и гидратом. При постановке задачи принято, что в результате отбора газа в пласте возникает три характерные области: ближняя, где поры заполнены газом и водой, промежуточная, в которой газ, вода и гидрат находятся в равновесии, и дальняя, которая заполнена газом и гидратом. Соответственно возникают две подвижные поверхности: между дальней и промежуточной областями, где начинается диссоциация гидрата, и между ближней и промежуточной областями, на которой заканчивается процесс разложения гидрата. В автомодельном приближении получено аналитическое решение данной задачи. Установлена возможность существования решений, согласно которым разложение газогидрата может происходить как на газ и воду, так на газ и лед. Показано, что в случае разложения на газ и лед необходимо вводить третью подвижную межфазную границу, на которой происходит образование льда и соответственно дополнительно рассматривать четвертую область, в которой одновременно сосуществуют в состоянии фазового равновесия газ, лед и гидрат. Определены критические значения давления, определяющего возникновение четвертой области, в которой одновременно сосуществуют в состоянии фазового равновесия газ, лед и гидрат.
ключевые слова: газовый гидрат, диссоциация, автомодельное решение, пористая среда
MATHEMATICAL MODEL OF GAS FROM GASHYDRATE STRATUM TAKING INTO ACCOUNT THE FORMATION OF ICE
Khasanov M.K., Dorovskaya M.S.
Sterlitamak branch of Bashkir State University, Sterlitamak, e-mail: hasanovmk@mail. ru
The features of the decomposition of gas hydrates during gas extraction from a porous medium, in the initial state and filled with gas hydrate. When the formulation of the problem, it is assumed that as a result of selection of gas formation occurs in three distinct areas: proximal, wherein the pores are filled with gas and water, an intermediate in which the gas hydrate and water in equilibrium, and the long, which is filled with gas and hydrate. Accordingly, there are two moving surfaces, between the distal and intermediate regions, where it starts to hydrate dissociation, and between the near and intermediate regions, which ends the process of hydrate decomposition. In similar approximation, an analytical solution of this problem. The possibility of the existence of solutions, according to which the gas hydrate decomposition can occur at both the gas and the water, so the gas and ice. It is shown that in the case of decomposition into gas and ice must enter the third mobile phase boundary at which the formation of ice and accordingly further contemplates a fourth region, which coexist in the phase equilibrium state of gas hydrate and ice. The critical value of the pressure determines the occurrence of the fourth area, which coexist in the phase equilibrium state of gas hydrate and ice.
Keywords: gas hydrate, dissociation, self-similar solution, porous medium
Актуальность исследований процессов образования и разложения газогидратов в пористых средах обусловлена все возрастающим интересом к проблеме газогидратов, связанным, в первую очередь, с признанием того факта, что в перспективе природные газогидраты могут стать новым источником углеводородного газа благодаря значительным ресурсам, неглубокому залеганию и концентрированному состоянию в них газа [1, 3]. В связи с этим за последние годы резко возрос интерес к исследованиям и разработкам технологий, позволяющих использовать газогидраты в виде альтернативного углеводородного сырья.
Общепринятые способы извлечения природного газа из газовых гидратов включают воздействие нагревания и/или пониженного давления на газовые гидраты с целью высвобождения природного газа. Однако данные способы требуют подвода
к системе значительного количества энергии, что ведет к высоким затратам на извлечение. Одним из вариантов решения данной проблемы могла бы быть разработка таких технологий добычи газа из газогидратов, при которых разложение гидрата происходит не на газ и воду, а на газ и лед. Это позволило бы уменьшить энергетические затраты на разработку газогидратных месторождений, поскольку удельная теплота фазового перехода гидрат-лед значительно ниже теплоты перехода гидрат-вода. Очевидно, что данные технологии должны быть подкреплены соответствующими расчетами по теоретическим моделям, адекватно отражающим реальные условия.
Математическая модель процесса разложения газогидрата в пористой среде на газ и воду для случая протекания фазовых переходов на фронтальной границе построена в работах [2, 6, 7]. Некоторые особенно-
¦ ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ ¦
2847
сти процесса образования газогидрата в пористой среде рассмотрены в работах [4,5]. В настоящей работе представлено решение задачи о диссоциации газогидрата в пористой среде при депрессионном воздействии в случае формирования объемной области фазовых переходов.
Постановка задачи и основные уравнения
Для описания процессов тепломассопе-реноса, сопровождающихся диссоциацией газогидрата в пористой среде, примем следующие допущения. Процесс однотемпера-турный, т. е. температуры пористой среды и насыщающего вещества (газа, гидрата или воды) совпадают. Гидрат является двух-компонентной системой с массовой концентрацией газа О. Кроме того, скелет пористой среды, газогидрат и вода несжимаемы и неподвижны, пористость постоянна, газ является калорически совершенным.
В рамках отмеченных допущений запишем для осесимметричной задачи уравнения сохранения массы, закон дарси, уравнения состояния газа и притока тепла (без учета баротермического эффекта): () 1 3
дt
г дг
ы о/
с К ЭР
дг
Р=РЛТ& gt-
д (с 8Т
(1)
д_ г дг
. дг,

Здесь т — пористость- р — давление, Т — температура- р. и — истинная плотность и насыщенность пор. -й фазы (. = к, I, g) — индексы к, I и g относятся к параметрам гидрата, воды и газа соответственно- и, к и ц — скорость, проницаемость и динамическая вязкость для газовой фазы- Ьк — удельная теплота разложения гидрата- с — удельная теплоемкость газа- рс и X -удельная объемная теплоемкость и коэффициент теплопроводности системы.
Зависимость коэффициента проницаемости для газа к от газонасыщенности можно задать на основе формулы Козени

^(^Ъи3).
Значения температуры и давления в области разложения гидрата связаны условием фазового равновесия [1]

Т = Т0 + Т* 1п
Р8~
(2)
о /
где Т0 — исходная температура системы, рк0 -равновесное давление, соответствующее исходной температуре, Т* - эмпирический параметр, зависящий от вида газогидрата.
При диссоциации газогидрата в пористом пласте возникают зоны, в которых газ, вода и гидрат могут находиться в различных состояниях. На поверхностях разрыва между этими зонами, где терпят скачки насыщенности фаз, а также потоки массы и тепла, выполняются соотношения, следующие из условий баланса массы и тепла:
т (5лРл (1 — + $ Р1) г (
(*)
= 0,
т
РА К -Г (,)) -РА^Г (
'-(*)
= 0, (3)
Г, дТ 1
X- =
дг
тР^А г,)
Здесь [у] - скачок параметра у на границе между зонами- г ^) — скорость движения этой границы- нижний индекс в скобках (5) относится к параметрам на границе между зонами. Температура и давление на этих границах полагаются непрерывными.
Будем полагать, что пласт в начальный момент времени насыщен газом и газогидратом, давление р0 и температура Т0 которых в исходном состоянии соответствуют термодинамическим условиям существования их в свободном состоянии (р0 & gt- р50). Пусть давление и температура изначально одинаковы во всем пласте
X = 0: Т = Т0, р = р0 (г & gt- О.
Пусть через скважину, вскрывшую пласт на всю толщину, отбирается газ с постоянным массовым расходом Q (на единицу высоты скважины) при условии отсутствия кондуктивного потока тепла на границе скважины. В результате отбора газа в пласте возникает три характерные области: ближняя (первая), где поры заполнены газом и льдом (водой), промежуточная (вторая), в которой газ, лед (вода) и гидрат находятся в равновесии, и дальняя (третья), которая заполнена газом и гидратом. Во второй зоне происхо-
дит диссоциация гидрата. Соответственно возникают две подвижные поверхности: между дальней и промежуточной областями, где начинается плавление гидрата, и между ближней и промежуточной областями, на которой заканчивается процесс разложения гидрата. Исходная гидратона-сыщенность пласта V равна гидратонасы-щенности третьей области и постоянна во всех точках этой области.
Рассматривая достаточно большие времена после начала отбора газа, когда радиус ближней (первой) области значительно превышает радиус скважины г можно полагать, что размер скважины слабо сказывается на особенностях протекания процесса. Тогда с учетом закона Дарси и уравнения состояния для газа условия на границе скважины имеют вид:
kg л
М7
2 Л
дг
Q,
дт 0 (
— = 0 (г дг v w
¦ 0, t & gt- 0).
Автомодельное решение
Сформулированная задача имеет автомодельное решение. Введем, автомодельную переменную ^) = r/¦^x (¦Г)t, где х (Т = и рс — коэффициент температуропроводности пласта.
В автомодельных переменных решение для распределения давления (после применения метода линеаризации Лейбензона) и температуры в каждой из областей может быть записано в виде:
а) в ближней области (0 & lt-? & lt- ?):
Рс1)
'--Pi)
ем л»)
пк,
(1)
1
Л т
ехр

4ti

Я
T = T
shm = 0-
(4)
б) в промежуточной области (? & lt-? & lt- ?):
2 2, / 2 2 Л P (2)=P (m)+(P (s)-P (m))
J, ехр
v 4fb& gt--

J? exp
«М ^
T (2) = T0 + T 1П (P (2) /Ps0), Sm =V-
V 4iW
т-(Л (2) -l)
я,
2(фАГ,+л (2)АГ)
-In
в) в дальней области (?m) & lt-? & lt-®(m)):
(my
со j
(pU —)J, ехр
2 2 i ^ Аз) =Ро+~
V 4т1(3)У

Г 1
J, еХР
«м q
. X
V 4tW


(Р (т) ~){ к еХР I '-
(Ре 4 2 р
W 1 Р (3) о У

f 1 (Г ^(3) 2


(5)
(6)
Здесь
Sh (3) = v.

К. j) = k0 j), PeW) =
PsqcAjоД& gt-
Фл
& quot- -, ДГ = -, р «- плотность газа при давлении р0.
P, o (l-v)
рс
¦ ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ ¦
2849
Параметры первой, второй и третьей областей снабжены нижними индексами в скобках. = 1, 2 и 3- нижний индекс в скобках (5) и относится к параметрам на границе между первой и второй областями, (т) — между второй и третьей зонами соответственно.
На основе решений (4) — (6) и условий баланса массы и тепла (3) на межфазных границах получена система трансцендентных уравнений, исходя из которой численно определяются автомодельные координаты данных границы и значения параметров на них.
Т0 = 275 К И
Рис. 1. Распределение по переменной? температуры, давления и гидратонасыщенности: Q = 0,035 кг/(мс) (1), Q = 0,045 кг/(мс) (2)
Результаты расчетов
На рис. 1 приведены распределения давления, температуры и гидратонасы-щенности для разных значений массового расхода отбора метана из гидратосодержа-щего пласта. Для параметров, характеризу-
ющих систему, приняты следующие значения: т = 0,1, О = 0,12, V = 0,2, р = 5 МПа,, Т* = 10 К, р80 = 3,3 МПа, 520 Дж/(К-кг), к0 = 10−14 м2, р. = р. = =^00 кг/м3, рс = 2,5−106 Дж/(Юкг),
с^ = 1560 Дж/(Ккг), X = 2 Вт/(мК),
ц = 10−5 кг/(мх), Ь)1 = 5 105 Дж/кг. В соответствии с рисунком, если массовый расход отбора газа достаточно велик (кривая 2), то температура на границе, разделяющей первую и вторую области, опускается ниже температуры замерзания воды, что соответствует «переохлаждению» воды. Следовательно, в этом случае модель с двумя подвижными межфазными границами не позволяет построить физически непротиворечивое решение. Для устранения данного противоречия необходимо вводить третью подвижную межфазную границу, на которой происходит образование льда, и соответственно дополнительно рассматривать четвертую область, в которой одновременно сосуществуют в состоянии фазового равновесия газ, лед и гидрат.
Рис. 2. Зависимость критического значения массового расхода от абсолютной проницаемости пласта: Т0 = 275 К, р0 = 5 МПа
На рис. 2 изображена зависимость критического значения массового расхода отбора газа от абсолютной проницаемости, определяющего возникновение двух протяженных областей разложения гидрата как на газ и воду, так и на газ и лед. Данная кривая разделяет плоскость (к Q) на две области, каждая из которых соответствует различным режимам разложения гидрата. В области под кривой существует решение с разложением гидрата на газ и воду. Область выше кривой отвечает формированию двух протяженных областей фазовых переходов с разложением газогидрата как на газ и воду, так и на газ и лед. При этом в соответствии с рисунком режим с образованием
льда характерен для низкопроницаемых пористых сред. Это обусловлено тем, что снижение проницаемости пласта при фиксированном массовом расходе отбора газа обуславливает значительное снижение давления непосредственно вблизи скважины, тем самым приводя к уменьшению давления и соответственно температуры на границе разложения гидрата ниже температуры замерзания воды.
Выводы
В работе исследованы особенности разложения газогидрата в пористой среде при отборе газа. Показано, что в зависимости от интенсивности отбора газа и исходных параметров системы разложение газогидрата может происходить либо на газ и воду, либо как на газ и воду, так и на газ и лед.
Исследование выполнено при финансовой поддержке Российского фонда фундаментальных исследований в рамках научного проекта 14−01−31 089 мола.
список литературы
1. Бык С. Ш., Макогон Ю. Ф., Фомина В. И. Газовые гидраты. — М.: Химия, 1980. — 296 с.
2. Васильев В. И., Попов В. В., Цыпкин Г. Г. Численное исследование разложения газовых гидратов, сосуществующих с газом в природных пластах // Изв. РАН. Механика жидкости и газа. — 2006. — № 4. — С. 127−134.
3. Истомин В. А., Якушев В. С. Газовые гидраты в природных условиях. — М.: Недра, 1992. — 236 с.
4. Хасанов М. К., Столповский М. В., Кильдибаева С. Р. Численное моделирование образования газогидрата в пористой среде при инжекции газа // Вестник Башкирского университета. — 2013. — Т. 18, № 4. — С. 969−972.
5. Khasanov M.K., Gimaltdinov I.K. and Stolpovsky M.V. Specific features of the formation of gas hydrates during the injection of a cold gas into a porous medium saturated with a gas and water // Theoretical Foundations of Chemical Engineering -2010, Vol. 44, № 4 — P. 424−431.
6. Tsypkin G.G. Mathematical model for dissociation of gas hydrates coexisting with gas in strata // Doklady Physics -2001, Vol. 46, № 11 — P. 806−809.
7. Tsypkin G.G. Analytical solution of the nonlinear problem of gas hydrate dissociation in a formation // Fluid Dynamics — 2007, Vol. 42, № 5 — P. 798−806.
References
1. Byk S. Sh., Makogon Ju.F., Fomina V.I. Gazovye gidraty., Moscow: Khimija, 1980, 296 p.
2. Vasilev V.I., Popov V.V., Cypkin G.G., Chislennoe issle-dovanie razlozhenija gazovyh gidratov, sosushhestvujushhih s gazom v prirodnyh plastah. Izv. RAN. Mehanika zhidkosti i gaza., 2006. no. 4. pp. 127−134.
3. Istomin V.A., Jakushev V.S. Gazovye gidraty v prirodnyh uslovijah, Moscow: Nedra, 1992, 236 p.
4. Khasanov M.K., Stolpovskij M.V., Kil'-dibaeva S.R., Chislennoe modelirovanie obrazovanija gazogidrata v poris-toj srede pri inzhekcii gaza. Vestnik Bashkirskogo universiteta, 2013, Vol. 18, no. 4, pp. 969−972.
5. Khasanov M.K., Gimaltdinov I.K. and Stolpovsky M.V., Specific features of the formation of gas hydrates during the injection of a cold gas into a porous medium saturated with a gas and water, Theoretical Foundations of Chemical Engineering, 2010, Vol. 44, no. 4, pp. 424−431.
6. Tsypkin G.G. Mathematical model for dissociation of gas hydrates coexisting with gas in strata, Doklady Physics, 2001, Vol. 46, no. 11, pp. 806−809.
7. Tsypkin G.G. Analytical solution of the nonlinear problem of gas hydrate dissociation in a formation, Fluid Dynamics, 2007, Vol. 42, no. 5, pp. 798−806.
Рецензенты:
Мустафина С. А., д.ф. -м.н., профессор, декан физико-математического факультета Стерлитамакского филиала Башкирского государственного университета, г. Стерлитамак-
Михайлов П. Н., д.ф. -м.н., профессор, заведующий кафедрой «Алгебра, геометрия и методика обучения математике» Стерли-тамакского филиала Башкирского государственного университета, г. Стерлитамак.
Работа поступила в редакцию 01. 04. 2015.

ПоказатьСвернуть
Заполнить форму текущей работой