ФАЗОВЫЕ ПЕРЕХОДЫ МЕТАНА ПРИ ИЗМЕНЕНИИ ТЕРМОДИНАМИЧЕСКИХ ПАРАМЕТРОВ ГОРНОГО МАССИВА
Для построения математической модели фазовых превращений метана примем во внимание, что образование метана в период накопления торфяника и постепенного погребения его под наносы последующих отложений происходило при температурах 150-300oС, когда сорбционная способность угля была близка к нулю. В дальнейшем, в процессе инверсии и понижении температуры, часть метана сорбировалась углем, часть оставалась в свободном состоянии как в трещинах и микропорах угля, так и в коллекторах вмещающих пород. Дальнейшее изменение термодинамических параметров угленосной толщи влекло за собой переход свободного газа в гидратированное состояние. Образование гидратов метана происходит либо при низких температурах (t=12-14 oС) при P=10 МПа, либо при высоком гидростатическом давлении, большем чем в современных условиях. Например, для равновесного состояния гидрата метана. [1]
. (1)
Образование гидратов в пористой среде при наличии центров кристаллизации имеет объемно-диффузионный хаpaктер. Одним из основных факторов, определяющих условия существования плотных гидратов в пористой среде, является упругость газа в водном растворе, находящемся в контакте с гидратом ρw и в гидрате ρH при заданной температуре. Величина ρH определяется составом гидрата и его температурой. Параметр ρw зависит от растворимости газа в воде при заданных давлении и температуре [2].
Условием существования гидрата является ρw ≥ ρH. Растворимость газа в воде, находящейся в контакте с гидратом, всегда ниже, чем в отсутствии гидрата [2].
Степень газонасыщенности водного раствора, контактирующего с гидратом, определяется диффузионными потоками газа:
а) потоком рассеяния в вышележащих горизонтах или в омывающие воды (Д1);
б) потоком из нижележащих пластов (Д2).
Гидрат не будет диссоциировать при Д2 ≥ Д1. Из этого условия следует, что чем глубже гидросодержащие породы, чем меньше и однороднее поры, тем ниже Д1, а следовательно, тем выше возможность сохранения гидрата.
При изменении горнотехнических условий происходит понижение давления газа ниже давления разложения гидрата при существующей пластовой температуре, т.е. ρw становится меньше ρH (ρw < ρH), что вызывает диссоциацию гидратной фазы, т.е. газ переходит в свободное состояние и существенно изменяет газодинамику угленосной среды.
Для построения математической модели фазовых переходов метана выделим элементарный объем горного массива V, ограниченный поверхностью G. Будем предполагать, что в этом объеме газ находится в трехфазном состоянии (свободном, сорбированном и гидратном). Тогда масса газа, заключенного в этом объеме, в момент времени t будет равен
(2)
где ρi - плотность газа в каждой фазе, r(x1,x2,x3) - радиус-вектор элемента объема dV.
В процессе движения из объема V через его границу G в единицу времени вытекает количество свободного газа, равное
(3)
где ρ1 - плотность свободного газа,
n - единичный вектор внешней нормали к поверхности G.
Тогда баланс газа за время dt будет
(4)
Для течений, не имеющих сильных разрывов, интегральные уравнения заменим дифференциальными, описывающими процесс фазового перехода метана. Понижение давления до величины разложения гидрата определило подвижную границу диссоциации газа .
В зоне I (0 ≤ x ≤ ) произошла десорбция газа и закончилось разложение гидрата и газ перешел в свободное состояние. В зоне II ( ≤ x ≤ ∞) процесс диссоциации еще не наступил и гидратная фаза сохраняется.
В одномерном случае процесс диссоциации описывается следующей краевой задачей:
Здесь = ρRT соответственно в зонах диссоциации и гидратирования,
D = λU - коэффициент конвективной диффузии, м2/сутки;
λ - параметр дисперсии, M;
γ = ηU - коэффициент гидратирования (диссоциации), сутки-1;
η - константа скорости гидратирования, м-1;
t1 - время окончания первой стадии, сутки.
Введем безразмерные координаты
(8)
1. Рассмотрим стадию диссоциации газа в зоне . Задача (5) в обозначениях (8) принимает вид
(9)
Решение задачи (9) в виде составного разложения по степеням e имеет вид [3]
(10)
Подставляя (10) в (9) и приравнивая члeны при одинаковых степенях e,получим
(11)
(12)
(13)
(14)
Решая (11)-(14) методом хаpaктеристик, найдем
(15)
Если , то При функция определяется из (12), (13) при конкретном задании функции .
Так как θ1 > θ0(h), то на хаpaктеристике функция терпит разрыв, т.е. разложение решения задачи (11) в виде (12) справедливо всюду, за исключением окрестности хаpaктеристики . В работе [3] приведено решение данной задачи в виде
(16)
Общее решение, справедливое во всей области фазовых превращений метана, можно получить методом аддитивного составления [3]
(17)
Формула (17) справедлива при
При и из (17) имеем
(18)
При этих условиях точное решение задачи (11) получено в виде
(19)
Незначительное различие точного (19) и приближенного решения (18) указывает на эффективность приближенных методов для решения данной задачи.
Из анализа полученного решения можно определить подвижную границу диссоциации газа . Граница движется противоположно направлению фильтрации. При большой скорости ведения горных работ граница может приближаться к кромке пласта. Для приблизительной оценки величины предположим, что давление на границе зоны диссоциации совпадает с давлением свободного газа в начальный период разложения гидрата, т.е. , тогда уравнение баланса газа на границе можно представить в виде
(20)
где определяется по формуле (1). Разделяя переменные в уравнении (20), найдем
(21)
Из формулы (21) следует, что гидратное давление снижается довольно резко до некоторого значения, соответствующего разложению гидрата, что свидетельствует о быстром его разложении и незначительных размерах зоны диссоциации ( ).
Расчеты по формуле (21) показывают, что зона диссоциации находится в зоне упругих деформаций угольного пласта. Этот факт говорит о том, что в этой области происходит начительный рост концентрации газа. Общее количество диссоциированного газа будет
(22)
Исходя из формулы состава кристаллогидрата ( ) и с учетом выражения (22) можно определить массу газа и число объемов газа в одном объеме гидрата. Например, для P0=1 МПа, τ=0 оС. Число объемов газа в единице объема гидрата составляет 101,44 кг/м3. Таким образом, один объем гидрата содержит более сотни объемов газа, это свидетельствует о том, что гидраты (особенно метана) отличаются значительными запасами внутренней энергии и высокой концентрацией газа. Следовательно, зоны разложения гидрата являются зонами повышенной газодинамической активности пласта.
СПИСОК ЛИТЕРАТУРЫ:
- Веригин Н.Н. О разложении гидрата газа в плате. //МЖГ.- 1982 - №4.-с.171 -173.
- Макогон Ю.Ф. Гидраты природных газов.- М.: Недра,1974.
- Ван-Дайк М. Методы возмущений в механике жидкости. М., Мир, 1967.
Статья в формате PDF 244 KB...
03 11 2024 18:24:53
Статья в формате PDF 220 KB...
02 11 2024 9:57:57
Статья в формате PDF 118 KB...
01 11 2024 4:16:12
Статья в формате PDF 111 KB...
29 10 2024 17:52:56
28 10 2024 13:54:25
Статья в формате PDF 101 KB...
26 10 2024 6:41:24
Статья в формате PDF 312 KB...
25 10 2024 10:34:32
Статья в формате PDF 293 KB...
24 10 2024 20:55:55
Статья в формате PDF 281 KB...
23 10 2024 0:39:52
Статья в формате PDF 102 KB...
22 10 2024 3:46:33
Статья в формате PDF 114 KB...
21 10 2024 13:52:40
Статья в формате PDF 108 KB...
20 10 2024 2:29:47
Статья в формате PDF 253 KB...
19 10 2024 0:37:31
Статья в формате PDF 277 KB...
18 10 2024 16:49:39
Проведен анализ эффективности курсового гетеросуггестивного воздействия на функциональное состояние ЦНС у женщин репродуктивного возраста. С помощью методов электроэнцефалографии и спектрального анализа вариабельности сердечного ритма получены достоверные данные о положительной динамике на центральном и вегетативном уровнях обеспечения психофизиологической устойчивости обследованных женщин. ...
17 10 2024 1:10:12
Статья в формате PDF 108 KB...
16 10 2024 3:35:22
Представлены результаты поисковых экспериментов по применению неразрушающего способа сертификации резонансной древесины ели на поленьях, основанного на измерении скорости звука поперек полена. ...
15 10 2024 13:29:42
Статья в формате PDF 312 KB...
14 10 2024 23:47:58
Статья в формате PDF 161 KB...
13 10 2024 13:13:15
Статья в формате PDF 117 KB...
12 10 2024 23:31:26
Статья в формате PDF 106 KB...
11 10 2024 14:50:54
Статья в формате PDF 121 KB...
10 10 2024 8:52:38
Статья в формате PDF 376 KB...
09 10 2024 16:31:40
Статья в формате PDF 109 KB...
08 10 2024 23:22:58
Статья в формате PDF 100 KB...
06 10 2024 14:43:20
Статья в формате PDF 123 KB...
05 10 2024 6:42:22
Статья в формате PDF 217 KB...
03 10 2024 22:45:25
Статья в формате PDF 101 KB...
02 10 2024 17:27:50
Статья в формате PDF 125 KB...
30 09 2024 17:18:41
Статья в формате PDF 109 KB...
29 09 2024 12:57:43
Статья в формате PDF 102 KB...
28 09 2024 8:36:12
Статья в формате PDF 135 KB...
27 09 2024 15:29:17
Статья в формате PDF 124 KB...
26 09 2024 19:22:35
Статья в формате PDF 122 KB...
25 09 2024 7:48:57
Еще:
Поддержать себя -1 :: Поддержать себя -2 :: Поддержать себя -3 :: Поддержать себя -4 :: Поддержать себя -5 :: Поддержать себя -6 :: Поддержать себя -7 :: Поддержать себя -8 :: Поддержать себя -9 :: Поддержать себя -10 :: Поддержать себя -11 :: Поддержать себя -12 :: Поддержать себя -13 :: Поддержать себя -14 :: Поддержать себя -15 :: Поддержать себя -16 :: Поддержать себя -17 :: Поддержать себя -18 :: Поддержать себя -19 :: Поддержать себя -20 :: Поддержать себя -21 :: Поддержать себя -22 :: Поддержать себя -23 :: Поддержать себя -24 :: Поддержать себя -25 :: Поддержать себя -26 :: Поддержать себя -27 :: Поддержать себя -28 :: Поддержать себя -29 :: Поддержать себя -30 :: Поддержать себя -31 :: Поддержать себя -32 :: Поддержать себя -33 :: Поддержать себя -34 :: Поддержать себя -35 :: Поддержать себя -36 :: Поддержать себя -37 :: Поддержать себя -38 ::