ФАЗОВЫЕ ПЕРЕХОДЫ МЕТАНА ПРИ ИЗМЕНЕНИИ ТЕРМОДИНАМИЧЕСКИХ ПАРАМЕТРОВ ГОРНОГО МАССИВА

1

• Авторы
• Файлы

Беспятов Г.А. Статья в формате PDF 127 KB

Для построения математической модели фазовых превращений метана примем во внимание, что образование метана в период накопления торфяника и постепенного погребения его под наносы последующих отложений происходило при температурах 150-300^oС, когда сорбционная способность угля была близка к нулю. В дальнейшем, в процессе инверсии и понижении температуры, часть метана сорбировалась углем, часть оставалась в свободном состоянии как в трещинах и микропорах угля, так и в коллекторах вмещающих пород. Дальнейшее изменение термодинамических параметров угленосной толщи влекло за собой переход свободного газа в гидратированное состояние. Образование гидратов метана происходит либо при низких температурах (t=12-14 ^oС) при P=10 МПа, либо при высоком гидростатическом давлении, большем чем в современных условиях. Например, для равновесного состояния гидрата метана. [1]

  .                                       (1)

Образование гидратов в пористой среде при наличии центров кристаллизации имеет объемно-диффузионный хаpaктер. Одним из основных факторов, определяющих условия существования плотных гидратов в пористой среде, является упругость газа в водном растворе, находящемся в контакте с гидратом ρ_w и в гидрате ρ_H при заданной температуре. Величина ρ_H определяется составом гидрата и его температурой. Параметр ρ_w зависит от растворимости газа в воде при заданных давлении и температуре [2].

Условием существования гидрата является ρ_w ≥ ρ_H. Растворимость газа в воде, находящейся в контакте с гидратом, всегда ниже, чем в отсутствии гидрата [2].

Степень газонасыщенности водного раствора, контактирующего с гидратом, определяется диффузионными потоками газа:

а) потоком рассеяния в вышележащих горизонтах или в омывающие воды (Д₁);

б) потоком из нижележащих пластов (Д₂).

Гидрат не будет диссоциировать при Д₂ ≥ Д₁. Из этого условия следует, что чем глубже гидросодержащие породы, чем меньше и однороднее поры, тем ниже Д_1, а следовательно, тем выше возможность сохранения гидрата.

При изменении горнотехнических условий происходит понижение давления газа ниже давления разложения гидрата при существующей пластовой температуре, т.е. ρ_w становится меньше ρ_H (ρ_w < ρ_H), что вызывает диссоциацию гидратной фазы, т.е. газ переходит в свободное состояние и существенно изменяет газодинамику угленосной среды.

Для построения математической модели фазовых переходов метана выделим элементарный объем горного массива V, ограниченный поверхностью G. Будем предполагать, что в этом объеме газ находится в трехфазном состоянии (свободном, сорбированном и гидратном). Тогда масса газа, заключенного в этом объеме, в момент времени t будет равен

                     (2)

где ρ_i - плотность газа в каждой фазе, r(x₁,x₂,x₃) - радиус-вектор элемента объема dV.

В процессе движения из объема V через его границу G в единицу времени вытекает количество свободного газа, равное

               (3)

где ρ₁ - плотность свободного газа,

n - единичный вектор внешней нормали к поверхности G.

Тогда баланс газа за время dt будет

                                         (4)

Для течений, не имеющих сильных разрывов, интегральные уравнения заменим дифференциальными, описывающими процесс фазового перехода метана. Понижение давления до величины разложения гидрата определило подвижную границу диссоциации газа .

В зоне I (0 ≤ x ≤ ) произошла десорбция газа и закончилось разложение гидрата и газ перешел в свободное состояние. В зоне II ( ≤ x ≤ ∞) процесс диссоциации еще не наступил и гидратная фаза сохраняется.

В одномерном случае процесс диссоциации описывается следующей краевой задачей:

Здесь  = ρRT соответственно в зонах диссоциации и гидратирования,

D = λU - коэффициент конвективной диффузии, м²/сутки;

λ - параметр дисперсии, M;

γ = ηU - коэффициент гидратирования (диссоциации), сутки^-1;

η - константа скорости гидратирования, м^-1;

t₁ - время окончания первой стадии, сутки.

Введем безразмерные координаты

                               (8)

1. Рассмотрим стадию диссоциации газа в зоне . Задача (5) в обозначениях (8) принимает вид

                              (9)

Решение задачи (9) в виде составного разложения по степеням e имеет вид [3]

                  (10)

Подставляя (10) в (9) и приравнивая члeны при одинаковых степенях e,получим

                               (11)

                (12)

                                   (13)

                             (14)

Решая (11)-(14) методом хаpaктеристик, найдем

                           (15)

Если , то  При  функция  определяется из (12), (13) при конкретном задании функции .

Так как θ₁ > θ₀(h), то на хаpaктеристике  функция  терпит разрыв, т.е. разложение решения задачи (11) в виде (12) справедливо всюду, за исключением окрестности хаpaктеристики . В работе [3] приведено решение данной задачи в виде

                 (16)

Общее решение, справедливое во всей области фазовых превращений метана, можно получить методом аддитивного составления [3]

                      (17)

Формула (17) справедлива при

При  и  из (17) имеем

                           (18)

При этих условиях точное решение задачи (11) получено в виде

               (19)

Незначительное различие точного (19) и приближенного решения (18) указывает на эффективность приближенных методов для решения данной задачи.

Из анализа полученного решения можно определить подвижную границу диссоциации газа . Граница  движется противоположно направлению фильтрации. При большой скорости ведения горных работ граница  может приближаться к кромке пласта. Для приблизительной оценки величины  предположим, что давление на границе зоны диссоциации совпадает с давлением свободного газа в начальный период разложения гидрата, т.е. , тогда уравнение баланса газа на границе  можно представить в виде

                            (20)

где  определяется по формуле (1). Разделяя переменные в уравнении (20), найдем

                        (21)

Из формулы (21) следует, что гидратное давление  снижается довольно резко до некоторого значения, соответствующего разложению гидрата, что свидетельствует о быстром его разложении и незначительных размерах зоны диссоциации ( ).

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

                  (22)

Исходя из формулы состава кристаллогидрата ( ) и с учетом выражения (22) можно определить массу газа и число объемов газа в одном объеме гидрата. Например, для P₀=1 МПа, τ=0 ^оС. Число объемов газа в единице объема гидрата составляет 101,44 кг/м³. Таким образом, один объем гидрата содержит более сотни объемов газа, это свидетельствует о том, что гидраты (особенно метана) отличаются значительными запасами внутренней энергии и высокой концентрацией газа. Следовательно, зоны разложения гидрата являются зонами повышенной газодинамической активности пласта.

СПИСОК ЛИТЕРАТУРЫ:

1. Веригин Н.Н. О разложении гидрата газа в плате. //МЖГ.- 1982 - №4.-с.171 -173.
2. Макогон Ю.Ф. Гидраты природных газов.- М.: Недра,1974.
3. Ван-Дайк М. Методы возмущений в механике жидкости. М., Мир, 1967.

---