CC BY
47
НАУКИ О ЗЕМЛЕ
УДК 622.831.322
ЗАВИСИМОСТЬ ВИДА ГАЗОДИНАМИЧЕСКОГО ЯВЛЕНИЯ ОТ ПРОНИЦАЕМОСТИ МАССИВА С УЧЕТОМ ГАЗОГИДРАТОВ
A CERTAIN TYPE OF GAS-DYNAMIC PHENOMENA DEPENDING ON THE POROSITY OF COAL SEAMS IN THE DECOMPOSITION OF GAS HYDRATES
Ким Татьяна Леонидовна,
канд.техн. наук, e-mail: tanyakim@list.ru, Kim Tatyana L., С.Sc. ( Engineering) Дырдин Валерий Васильевич, доктор техн. наук, профессор, e-mail: wdl941@live.ru Dyrdyn Valeriy V., Dr.Sc. ( Engineering), professor
Кузбасский государственный технический университет имени Т.Ф. Горбачева, 650000, Россия, г. Кемерово, ул. Весенняя, 28
Т. F. Gorbachev Kuzbass State Technical University, 28 street Vesennyaya, Kemerovo, 650000, Russian Federation
Аннотация: Статья посвящена определению вида газодинамических явлений при возможной диссоциации газогидратов с учетом проницаемости угольного массива. Установлена закономерность распределения газового давления в краевой зоне угольного пласта после диссоциации газогидратов в зависимости от гидратонасыщенности и с учетом изменяющегося коэффициента проницаемости краевой зоны угольного пласта. Установлено, что вид газодинамического явления при наличии газогидратов в поровом пространстве каменных углей зависит от критического значения проницаемости краевой зоны.
Abstract: The publication is devoted creation of a method of the forecast gas-dynamic phenomena occurring under the influence of gas hydrates in coal seams with application of a method of mathematical modelling. For the first time the minimum humidity at which formation of gas hydrates in coal seams is is possible is found. During work there was a definition of speed of a dissociation of hydrates of gas depending on a size of corpuscles by a method of mathematical modelling. Also conformity is found between velocity of move of boundary line of a dissociation of gas hydrates in space of pores of coal and velocity of extending of a wave of crushing. The interconnection of allocation of gas pressure in a band of edge of a layer of coal after a dissociation of gas hydrate depending on its saturation capacity and taking into account a changing coefficient of permeability is spotted. Including it has been spotted that the kind gas-dynamic processes in the presence of gas hydrates in space of pores of coals depends on critical value ofpermeability in a band of edge of a coal seam.
Ключевые слова: Газодинамические явления, скорость диссоциации, газогидраты, коэффициент проницаемости, внезапный выброс.
Keywords: gas-dynamic phenomena, the rate of dissociation, gas hydrates, the permeability coefficient, sudden emission.
Экспериментально установлено [1], что при определенных термодинамических условиях (температура и давление) и влажности в поровом пространстве угольных пластов могут образовываться газогидраты метана. В силу малых размеров газо-гидратных частиц скорость их разложения достаточна высока и при 4 нм может составлять 28 мкм/с. Разложение газогидратных частиц происходит в случае, если давление в системе газ-массив при определенной температуре может скачком уменьшиться до значений меньших равновесных. В таких случаях на участке разложения газовое давление может скачком увеличиться. С этой целью, используя методы численного и ма-
тематического моделирования, рассчитаем возможное приращение давления газа впереди забоя подготовительной выработки с учетом гидратонасыщенности угольного пласта и коэффициента проницаемости. Для этого воспользуемся дифференциальным уравнением Буссинеска для одномерного изотермического течения, полученного Л. С. Лейбензоном [2]:
д1 2т/иГ дх где га- коэффициент пористости среды; р -давление газа; цг- динамическая вязкость газа;
4
Т.Л. Ким, В.В. Дырдин
к - коэффициент проницаемости.
При решении задачи принято, что в точке максимума зоны опорного давления имеем минимум проницаемости массива. За зоной влияния выработки коэффициент проницаемости принимает постоянное значение к^. Другие граничные условия для коэффициента проницаемости: а) проницаемость на забое выработки к(х) | _Q = к0,
т. е. начало оси х совпадает с забоем выработки, а х > 0 - в глубину массива;
б) *з(*) Ц = М*) | . = к ' где * = рас"
стояние от забоя выработки до начала зоны, где происходит диссоциация газогидратов, т. е. принято, что диссоциация газогидратов начинается на границе зон неупругих (х3 )и упругих деформаций (х4).
В такой постановке данная задача является автомодельной, и для ее решения используем явные разностные схемы и математический пакет Maple 14, то есть вместо производных в дифференциальном уравнении будем использовать их конечноразностные аппроксимации. Аппроксимируя частные производные уравнения (1) конечными разностями, получаем систему линейных алгебраических уравнений для определения давления.
Предварительно определим число молей газа v, содержащегося в гидратном состоянии в угольном пласте [3]:
^ _ ^угля О гидрата ^
Меоды»
>
W р
где ft =-—— - гидратонасыщенность; плот-
Ph
ность воды pw = 1000 кг/м3; 8=0.126 - природная масса газа в единице массы газогидрата; плотность гидрата природного газа ph = 883 кг/м3; W - природная влажность угля; \1воды - молярная масса воды; п - число молекул воды, приходящихся на одну молекулу газа в элементарной ячейке газогидрата формы KC-I.
В качестве примера рассчитаем распределение давления газа после диссоциации газогидратов для условий пласта Владимировский шахты «Северная». Результаты расчета приведены на рис.1.
Дополнительный объем свободного газа зависит от протяженности зоны, содержащей газогидраты и гидратонасыщенности. Как следует из полученных результатов, добавочное давление газа в зоне разложения газогидратов увеличивается в 2-2.5 раза, что влияет на увеличение газового давления впереди забоя подготовительной выработки и может приводить к внезапным выбросам угля и газа и внезапным газовыделениям.
Для определения вида ГДЯ (внезапный выброс или загазирование) рассчитаем коэффициент
выбросоопасности пласта Я" по методике профессора В. И. Мурашёва [4-6] в зависимости от коэффициента проницаемости массива с учётом прироста давления вследствие диссоциации газогидратов метана, а по точке перехода кривой через ноль определим критическое значение проницаемости ккр для данного пласта, которое и позволит разграничить области, опасные по внезапным выбросам и возможному загазированию выработок.
Р, 106Па
Рис. 1. Распределение давления газа р после начала диссоциации газовых гидратов и изменяющейся проницаемости краевой зоны угольного пласта: 1 - исходное распределение газового давления; 2 - при гидратонасыщенности 1.9%; 3 - при гидратонасыщенности 2.6%; 4 - при гидратонасыщенности 3.9 %; 5 - при гидратонасыщенности 5.2 %, и- давление газа за зоной влияния вы-
1 Шал
работки
Зависимость Rn = f(r), построенная для
условий Кузбасса, представлена на рис. 2.
Для практического использования полученных результатов для разграничения участков угольного пласта, опасных по внезапным выбросам угля и газа, и загазированию выработок нами предлагается следующая методика.
Впереди забоя подготовительных выработок не образуется значительных концентраций напряжений, поэтому примем экспоненциальный закон уменьшения проницаемости угольного пласта вглубь массива:
к(х) = к0е~ЬЫх\ (2)
где k0- коэффициент газопроницаемости угля в состоянии полной разгрузки, м2 (определяется в лабораторных условиях по В. В. Ходоту); Ъ - некоторая постоянная, зависящая от компрессионных свойств угольного пласта; cr(x) - напряжения в рассматриваемой точке.
п s R, 10 Па/м
N
\
\ • вне !зап ные вы£ jpoc ы
\ ►
Ч ч
\ »
о.с 01 о,с 02 0.( D3 0.C оШ 5Ш 06 o.i оЛЯ иш !оШ 4 И о.с ТТЛ! то 140.( Т5Ю i7 k.
не опа< ;но 4 i £ > i 4
3 агаз 1_1 иро 1_ ват L—— ie в Mpal >ОТК и
Рис. 2. Зависимость вида ГДЯ от проницаемости массива на границе раздела упругих и неупругих
Зная к0 и <у(х), по формуле (2) рассчитаем
проницаемость в любой точке впереди забоя подготовительной выработки, в том числе и на границе зон упругих и неупругих деформаций. Напряжения на забое и в максимуме опрного давления определим по известной методике ВНИМИ [7,8]. Закон изменения напряжений в массиве от забоя до максимума примем линейным. Расстояние от забоя до зоны упругих деформаций определим одним из геофизических методов, например, методом скважинной электрометрии по точке перегиба экспериментальной кривой. Приведенная методика позволит определить значения проницаемости угольного пласта на границе раздела неупругих и упругих деформаций массива, т. е. где возможна начальная диссоциация газогидратов.
При значениях коэффициента проницаемости к на границе раздела упругих и неупругих деформаций, превышающих критическое значение для данного угольного пласта, происходит загазиро-вание выработки (так как происходит интенсивная фильтрация газа в сторону забоя выработки), а
деформаций.
при к < кКр - формирование выбросоопасной ситуации.
В качестве примера для пласта «Владимиров-ский» шахты Северная рассчитано, что значение коэффициента проницаемости на расстоянии /=3.5
-15 „2
м от забоя равно 3-10 '"м", а критическое значение для этих условий составляет 10-10"15 м2, т.е. к
(3.5)= 3-10"15м2 < ккр = 10-10~15м2, поэтому данный участок опасен по внезапным выбросам угля и газа (рис. 2).
Таким образом, диссоциация газогидратов может существенно влиять на формирование вида газодинамического явления, поэтому для повышения безопасности подземной отработки угольных пластов нужно на плане ведения горных работ в пределах шахтного поля и горизонта выделять по термодинамическим параметрам и влажности потенциально опасные участки по наличию газогидратов, что позволит заранее разработать горнотехнические мероприятия по переходу забоями выработок данных участков.
СПИСОК ЛИТЕРАТУРЫ
1. Ким, Т. J1. Образование твердых растворов природного газа в угольной матрице //Вестник Кузбасского государственного технического университета. - 2013. - № 6. - С. 9 -12.
2. Максимов, А. М. Оценка возможности выбросов газа при разложении газовых гидратов в пласте / А. М. Максимов, В. С. Якушев, Е. М. Чувилин // Доклады Академии наук. - 1997. - том 352, № 4. - С. 532-534.
3. Ким, Т. JI. Влияние твердых растворов природного газа на газодинамические процессы впереди забоя подготовительной выработки/ Т. JI. Ким, В. В. Дырдин, А. А. Мальшин, С. А. Шепелева // Вестник Кузбасского государственного технического университета.- 2012. - №3. - С. 12-15.
4. Мурашев, В. И. Методическое руководство по прогнозу зон, потенциально опасных по газодинамическим явлениям, при проведении подготовительных выработок с учетом твердых растворов природного газа по типу газогидратов» / В. И. Мурашев, Т. JI. Ким, В. В. Дырдин. - Кемерово, 2014. - 19 с.
5. Инструкция по безопасному ведению горных работ на пластах, опасных по внезапным выбросам угля (породы) и газа (РД 05-350-00). Предупреждение газодинамических явлений в угольных шахтах (Сборник документов) [Текст]. - Москва: Государственное предприятие НТЦ по безопасности в промышленности Госгортехнадзора России, 2011. - С. 120-303
6. Тарасов, Б. Г. Определение потенциально выбросоопасных участков впереди подготовительных
6
Т.Л. Ким, В.В. Дырдин
выработок по электрическим параметрам среды и естественного геоэлектрического поля / Б. Г. Тарасов, В.В. Дырдин, В. И. Мурашев, В. В. Иванов, Г. В. Кроль / Подземная разработка мощных угольных пластов. - Сб. научн. тр. / Кузбасск. полит, ин-т. Вып. 2. - Кемерово. - 1974.
7. Расчет и экспериментальная оценка напряжений в целиках и краевых частях пласта (методическое указание). ВНИМИ.- JL - 1973.
8. Яновская М.Ф., Премыслер Ю.С. Номограммы для расчета газовыделения при разрушении угля. - М.: ИГД им. А.А.Скочинского, 1967. 171 с.
REFERENCES
1. Kim, T. L. Obrazovanie tverdyh rastvorov prirodnogo gaza v ugol'noj matrice // Vestnik Kuzbasskogo gosudarstvennogo tekhnicheskogo universiteta. - 2013. - № 6. - S. 9 -12.
2. Maksimov, A. M. Ocenka vozmozhnosti vybrosov gaza pri razlozhenii gazovyh gidratov v plaste / A. M. Maksimov, V. S. Jakushev, E. M. Chuvilin // Doklady Akademii nauk. - 1997. - torn 352, № 4. - S. 532 -534.
3. Kim, T. L. Vlijanie tverdyh rastvorov prirodnogo gaza na gazodinamicheskie processy vperedi zaboja podgotovitel'noj vyrabotki/ T. L. Kim, V. V. Dyrdin, A. A. Mal'shin, S. A. Shepeleva // Vestnik Kuzbasskogo gosudarstvennogo tekhnicheskogo universiteta.- 2012. - №3. - S. 12-15.
4. Murashev, V. I. Metodicheskoe rukovodstvo po prognozu zon, potencial'no opasnyh po gazodinamich-eskim javlenijam, pri provedenii podgotovitel'nyh vyrabotok s uchetom tverdyh rastvorov prirodnogo gaza po tipu gazogidratov» / V. I. Murashev, T. L. Kim, V. V. Dyrdin. - Kemerovo, 2014. - 19 s.
5. Instrukcija po bezopasnomu vedeniju gornyh rabot na plastah, opasnyh po vnezapnym vybrosam uglja (porody) i gaza (RD 05-350-00). Preduprezhdenie gazodinamicheskih javlenij v ugol'nyh shahtah (Sbornik dokumentov) [Tekst]. - Moskva: Gosudarstvennoe predprijatie NTC po bezopasnosti v promyshlennosti Gos-gortehnadzora Rossii, 2011. - S. 120-303
6. Tarasov, B. G. Opredelenie potencial'no vybrosoopasnyh uchastkov vperedi pod-gotovitel'nyh vyrabotok po jelektricheskim parametram sredy i estestvennogo geojelektricheskogo polja / B. G. Tarasov, V. V. Dyrdin, V. I. Murashev, V. V. Ivanov, G. V. Krol' / Podzemnaja razrabotka moshhnyh ugol'nyh plastov. - Sb. nauchn. tr. / Kuzbassk. polit, in-t. Vyp. 2. - Kemerovo. - 1974.
7. Raschet i jeksperimental'naja ocenka naprjazhenij v celikah i kraevyh chastjah plasta (metodicheskoe ukazanie). VNIMI.-L. - 1973.
8. Janovskaja M.F., Premysler Ju.S. Nomogrammy dlja rascheta gazovydelenija pri razrushenii uglja. -M.: IGD im. A.A.Skochinskogo, 1967. 171 s.
Поступило в редакцию 15.03.2016 Received 15 March 2016
