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

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

УДК 622:51.001.57

Т. Л. Ким, В. В. Дырдин, В. А. Белков

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССА ДИССОЦИАЦИИ КРИСТАЛЛОГИДРАТА В КРАЕВОЙ ЗОНЕ УГОЛЬНОГО ПЛАСТА

Газы при определенном давлении и температуре могут образовывать кристаллогидраты, структура которых может быть различной. Гидраты могут стабильно существовать в широком диапазоне давлений и температур. Например, гидрат

— 8

метана существует при давлениях от 2 • 10 Па

до 2 • 103 MПа и температурах от 70 до 350 К.

Газовый гидрат - образование, в котором молекулы газа заключены в кристаллические ячейки, состоящие из молекул воды, удерживаемые водородной связью.

Для образования гидрата необходимыми условиями являются: наличие газа, воды и давление, превышающее равновесное значение при определенной температуре. Химические связи между молекулами отсутствуют. Молекулы воды объединены водородной связью, легко разрушающейся при понижении давления или повышении температуры.

В настоящее время известно более десяти структур газогидратов, существующих при различных давлениях и температурах. Большая часть новых структур выявлена группой ученых Сибирского отделения РАН[1].

Газогидраты также могут образовываться и находиться в угольных пластах при температурах, больших 0 градусов по Цельсию, но давление при этом должно быть равно или превышать равновесное значение.

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

При разработке некоторых угольных пластов происходят внезапные газовыделения, внезапные выбросы угля и газа и другие подобные явления, которые носят характер катастроф, приводят к

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

Диссоциация кристаллогидратов в уголь-

ном пласте может происходить при резком понижении давления. На рис.1. представлено распределение газового давления впереди забоя подготовительной выработки в естественных условиях до процесса диссоциации кристаллогидратов. Зона 1

- зона ламинарной фильтрации, которая практически совпадает с зоной влияния выработки. Зона 2

- зона, в которой давление газа практически постоянно; 3 - исходное распределение газового давления; 4 - распределение давления газа при увеличении трещиноватости вследствие разрушения части массива; 5 - зона разрушения.

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

Р.Па

Угольный пласт Х,М

V ©

5

У///////////ШШ////////А

Рис.1. Распределение газового давления впереди забоя подготовительной выработки

Уравнение движения метана в пласте угля одномерного изотермического течения описывается уравнением [ 3]:

Ф„

кп

, п=1,2 ,

Э/ 2тп^Г Эх2

где т- коэффициент динамической пористости среды; р -давление; /аГ - динамическая вязкость газа; к-коэффициент проницаемости.

20

Т. Л. Ким, В. В. Дырдин, В. А. Белков

Найдем закон изменения давления газа в областях 1,2 (рис.1). Примем следующие начальное и граничные условия:

р(Х,0) = р0 , при х<0 - начальное давление в выработке;

р 1

[Р,+р/ - 2 Рі/ + Рі-1,/ ] '

Рі

І +1

+ Р,

.

Формула (2) позволяет найти решение в данный момент времени через решение в предыдущий момент.

В результате расчетов при Т=278К получаем распределение газового давления в зонах 1 и 2 для момента времени 1=10с после начала процесса диссоциации кристаллогидратов. Расчет произведен при следующих значениях: рр = 6 -105 Па;

рГ = 12 -105 Па; к1 = 0,01дарси ; к2 = 0,001дарси ; р0 = 105 Па ; тх = 0,10; т2 = 0,05 . 1ЛГ = 0,001 Па • с ;

Из полученных результатов следует, что давление в первой зоне за это время возрастает до величины давления рр, которое составляет для данных условий 6 105 Па на расстоянии 2,5 м от забоя выработки по простиранию пласта, вследствие чего и прекращается процесс диссоциации кристаллогидратов.

Во второй зоне давление от рр возрастает до

давления рГ, соответствующего начальному состоянию. р,105Па-12-

где , -граница разложения кристаллогидратов;

Р 1( Х ,0) = Р 2 (х ,0) = Р(^ , г) = Р г , при х>а, а- протяженность зоны влияния выработки.

Для решения данной задачи используем явные разностные схемы и математический пакет Мар1е14.

Первоначально строим прямоугольную сетку методом конечных разностей, где по оси х будем откладывать значения давлений в зависимости от координаты точки, а по оси у -значения давления в зависимости от времени.

х- = ]Н, у = 0,1,2..п ;г, = ,к, , = 0,1,2..т ;

Заменим частную производную по времени в уравнении (1) конечно-разностной аппроксимацией:

= 1 [р(х, 1 + к) - Р(х, 1)]= 1 [р<,[+1 - Р,-]

Вторую производную по координате в уравнении (1) заменим выражением:

д 2 2 1

-~РГ = —г [р (х + Ъ, г) - 2 р (х, г) + р (х - Ъ, г) ]2 дх Ъ

= ]_

Ъ2

Разрешим полученное выражение относительно значений функции на верхнем временном слое:

1 [р,,у+1 - р,- ]= ЪТ [р,+1,у - 2 р,,у + р,-1,у ]

0 1 2.5 10 х,м

Рис.2. Распределение давления газа в краевой зоне пласта после начала диссоциации кристаллогидратов

Граница разложения кристаллогидрата в первоначальный момент расположена вблизи забоя и с течением времени перемещается вглубь массива. Скорость ее перемещения найдем из выражения:

V

где ' - содержание газа в гидратном состоянии в единице объема пористой среды, зависящее от пористости и влажности угольного пласта; р -

давление разложения гидрата при пластовой температуре; к 1, к2- коэффициенты проницаемости,

соответственно, в областях 1 и 2.

Зависимость скорости перемещения границы разложения гидрата от времени представлен на рис.З.

гдР2

Р0

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

Таким образом, путем решения дифференциальных уравнений второго порядка получены за-

р

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

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

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

1. Макогон, Ю. Ф. Гидраты природных газов /Ю. Ф. Макогон. - М.: Недра, 1974. - 208 с.

2. Шепелева С. А., Дырдин В. В. Расчет газовыделений в скважину при предварительной дегазации угольных пластов / С. А. Шепелева, В. В. Дырдин // Нетрадиционные и интенсивные технологии разработки месторождений полезных ископаемых: сб. науч. статей под общей ред. В. Н. Фрянова. - Вып. 1. -Новокузнецк: СибГИУ, 2008. - С. 59 - 64.

3.Веригин Н. Н., Хабибуллин И. Л., Халиков Г. А. Линейная задача о разложении гидратов газа в пористой среде. Изв. АН СССР, МЖГ,№1, 1980. С.174-177.

Авторы статьи:

Ким

Татьяна Леонидовна, ассистент каф. физики КузГТУ, e-mail: tanyakim@list.ru.

Дырдин Валерий Васильевич, докт.техн.наук, проф., каф. физики КузГТУ, тел.: (+7-384-2)39-63-71.

Белков

Виталий Александрович, студент гр.МА-102 КузГТУ, тел.: +7-908-94б-07-57.