Подать статью
Стать рецензентом
Том 245
Страницы:
531-538
Скачать том:

Особенности проектирования вентиляции угольных шахт, применяющих камерно-столбовую систему разработки

Авторы:
С. С. Кобылкин1
А. Р. Харисов2
Об авторах
  • 1 — д-р техн. наук профессор НИТУ "МИСиС", Горный институт ▪ Elibrary ▪ Scopus ▪ ResearcherID
  • 2 — заместитель командира взвода «Военизированного горноспасательного отряда Восточной Сибири» ФГУП «ВГСЧ» Военизированный горноспасательный отряд Восточной Сибири
Дата отправки:
2020-05-12
Дата принятия:
2020-09-22
Дата публикации:
2020-11-24

Аннотация

Безопасность ведения горных работ на угольных шахтах по аэрологическим факторам зависит от качества принятых и реализованных проектных решений по вентиляции. Действующее на данный момент «Руководство по проектированию вентиляции угольных шахт» не учитывает особенности применяемых в России камерно-столбовых систем разработки. Это приводит к повышению опасности возникновения взрывов, пожаров и загазирований. Подробно изученный зарубежный опыт проектирования вентиляции для рассматриваемых систем разработки угольных месторождений позволил сформулировать рекомендации по организации схемы проветривания угольных шахт, применяющих камерно-столбовую систему разработки, и порядок проветривания при многоштрековой подготовке. Наблюдения показали, что применение существующего в России порядка проветривания добычных участков с камерно-столбовой системой разработки усложняет проведение аварийно-спасательных работ. Низкие скорости и разнонаправленное движение воздуха, затрудненный отток тепла, оставление целиков угля повышают риск возникновения и позднего обнаружения эндогенного пожара. Результаты численного моделирования показали, что установка (параллельно штрекам) вентиляционных сооружений в межкамерных целиках позволит повысить надежность проветривания за счет перевода схемы проветривания из сложной диагональной в сложную параллельную. Это также позволит снизить количество воздуха, необходимое для добычного участка, и общее аэродинамическое сопротивление. Проведенные исследования позволили сформулировать требования к порядку проектирования вентиляции угольных шахт, применяющих камерно-столбовую систему разработки, заключающиеся в порядке отработки блоков в панели, а также дополнительном применении вентиляционных сооружений (легких перемычек или парусов).

Ключевые слова:
безопасность вентиляция угольная шахта проектирование камерно-столбовая система разработки мо-делирование схема проветривания
10.31897/PMI.2020.5.4
Перейти к тому 245

Введение. В настоящее время в России работают несколько угольных шахт, на которых применяется камерно-столбовая система разработки (например, шахты «ДенисовскаяЦентральная» и УК «Межегейуголь»). Данные системы разработки распространены при добыче негорючих полезных ископаемых. При добыче угля возникает целый ряд вопросов по обеспечению аэрологической безопасности в части выбора рациональных схем и режимов проветривания, вентиляционных сооружений и порядка контроля параметров рудничной атмосферы. Данные вопросы тесно связаны с порядком ведения горных работ и отработки месторождения. При этом нормативных требований, методических документов по проектированию и контролю вентиляции для таких систем разработки в России не разработано. Для обеспечения безопасности горных работ, а также горноспасательных работ с учетом предстоящих существенных изменений в нормативно-правовой базе и критериев «регуляторной гильотины» [7] необходимы исследования по параметрам проветривания угольных шахт, применяющих камерно-столбовую систему.

Постановка проблемы. Камерно-столбовая система разработки на угольных шахтах активно применяется в США и Австралии [14, 14]. Шахтное поле разделяется на панели, те, в свою очередь, отрабатываются проходческими комбайнами, откатка осуществляется, как правило, самоходными вагонами и конвейерами. Целики вынимаются обратным ходом с применением самоходных крепей. Проветривание таких выемочных участков осуществляется путем разделения панели на две или три части (рис.1). По одним горным выработкам воздух подается, а по другим выдается исходящая струя воздуха, при этом допускается поддерживать несколько выработок в «нейтральном» состоянии проветривания, без определенного направления движения воздуха, с соблюдением минимальной скорости движения воздуха, равной 250 FPM (1,27 м/с) [9].

Рис.1. Общий вид схем проходки при многоштрековой подготовке (а) и проветривания добычных участков (б) шахт США и Австралии 1 – перемычка; 2 – свежая струя воздуха; 3 – исходящая струя воздуха; 4 – конвейер; 5 – монорельс; 6 – вентиляционный парус
Рис.2. Схема проветривания призабойной части с применением перегородки (а) и вентилятора местного проветривания (б) 1 – перемычка; 2 – свежая струя воздуха; 3 – исходящая струя воздуха; 4 – конвейер; 5 – монорельс; 6 – вентиляционный парус; 7 – проходческий комбайн; 8 – перегородка; 9 – вентилятор местного проветривания и воздуховод

Особенность таких схем проветривания заключается в обеспечении воздухом призабойной части путем использования продольных перегородок (рис.2, а) или установкой в непосредственной близости вентиляторов местного проветривания (ВМП) (рис.2, б). При этом наблюдается четкое направленное движение свежей и исходящих струй воздуха. По расчетам W.R.Kennedy [14] данная схема позволяет снизить аэродинамическое сопротивление выемочного участка в девять раз с 0,09 до 0,01 km. Аэродинамическая схема представлена сложным параллельным соединением и является возвратноточной. Добычные и проходческие работы совмещены и имеют единый фронт.

На российских угольных шахтах схема проветривания имеет ряд существенных отличий:

1. Вначале одним или двумя спаренными штреками нарезается выемочный участок (выемочный столб или панель).Участок делится на блоки, которые в произвольном порядке отрабатываются проходческими комбайнами, при этом нет единого фронта добычных работ.

2. Воздух в забои подается путем направления основной струи в отрабатываемый блок за счет установок вентиляционных окон в штреках или по воздуховодам вентиляторами местного проветривания, установленных в штреках.

Рис.3. Вентиляционный план шахты УК «Межегейуголь» Условные обозначения см. на рис.1

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

Исследования проводились на шахте УК «Межегейуголь», которая расположена в Республике Тыва. Разрабатываемый пласт имеет наклонное падение, что обуславливает вскрытие его наклонными стволами. Одной из отличительных особенностей данного горного предприятия является то, что его проектировала австралийская компания Snowden, именно ее специалистами была предложена камерно-столбовая система разработки. Месторождение на данном этапе отработки вскрыто четырьмя наклонными стволами, при этом общая схема проветривания – центрально-фланговая. На шахте применятся нагнетательный способ проветривания (рис.3).

По данным геологических изысканий, месторождение имеет несколько зон с выжженным из-за самовозгорания углем. Это связано с наклонным падением пласта, выходом его на поверхность и наличием аэродинамической связи с поверхностью. При полном обрушении количество трещин, связывающих пласты угля с поверхностью, может быть существенно увеличено, что значительно повысит риск возникновения эндогенных пожаров [8, 17]. В работе находится один вентилятор главного проветривания ВО-A22-14ВР, воздух подается по западному вентиляционному стволу. Исходящая струя из шахты выдается через западный конвейерный, западный транспортный и западный фланговый стволы. Метаноносность разрабатываемого пласта с глубиной ведения горных работ увеличивается с 4 до 11 м3/т. При этом проветривание выемочных блоков предусматривается за счет общешахтной депрессии. Дополнительно для вентиляции тупиковых частей горных выработок с общей протяженностью более 10 м применяются вентиляторы местного проветривания. Они устанавливаются на штреках, оконтуривающих добычной участок (рис.3).

По данным геологических изысканий, месторождение имеет несколько зон с выжженным из-за самовозгорания углем. Это связано с наклонным падением пласта, выходом его на поверхность и наличием аэродинамической связи с поверхностью. При полном обрушении количество трещин, связывающих пласты угля с поверхностью, может быть существенно увеличено, что значительно повысит риск возникновения эндогенных пожаров [8, 17]. Составленная условная аэродинамическая схема (рис.4) позволяет наглядно оценить количество зон с трудно прогнозируемым направлением воздуха. Внутри блоков горные выработки проветриваются с низкими скоростями движения воздуха. Проведенные натурные измерения показали, что скорость движения воздуха менее 0,25 м/с. Периодически за короткий промежуток времени в междукамерных целиках происходит опрокидывание вентиляционных струй. Это связано с применяемым вспомогательным транспортом для доставки людей (машины марки Driftrunner) и грузов (погрузочно-доставочные машины CL210, дизелевозный манипулятор DMZ50F). Заезжая в горные выработки выемочного участка и блока, они перекрывают часть сечения, что приводит к увеличению аэродинамического сопротивления за счет лобового сопротивления.

Рис.4. Упрощенная схема проветривания шахты УК «Межегейуголь» 1 – вентилятор главного проветривания; 2 – вентилятор местного проветривания; 3 – ветвь со свежей струей воздуха; 4 – ветвь с исходящей струей воздуха

Также следует отметить, что для увеличения производительности по добыче угля отработка панели на шахте иногда проводится единовременно в нескольких блоках (см. рис.3). Это не исключает возможность загазирования горных выработок за счет попадания метана через конвейерный штрек, где газовый контроль осуществляется в соответствии с ФНиП только на исходящей с участка и у вентиляторов местного проветривания, установленных на вентиляционном штреке. Об изменении порядка контроля сообщалось ранее [3, 12], но для схем с камерно-столбовой системой разработки таких предложений на обсуждение не выносилось.

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

Так 5 августа 2017 г. на шахте УК «Межегейуголь» произошел пожар с задымлением горных выработок. На поверхность было выведено 93 человека. Из-за сложной схемы проветривания первоначально была введена в действие позиция № 17 (пожар в конвейерном штреке 2-2), но после проведения разведки выяснилось, что пожар произошел в ходке 3-15. С момента начала аварии до обнаружения очага пожара прошло 89 минут. Причиной пожара явилось возгорание конвейерной ленты. Дым распространился по всем горным выработкам, что усложнило поиск места очага. Однако горноспасателям удалось локализовать распространение пожара и не допустить возгорания целиков угля в выработанном пространстве.

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

Методология и анализ результатов исследований. Для оценки эффективности проветривания и обеспечения безопасности по аэрологическому фактору для сложных диагональных схем, которыми представлены выемочные участки, были использованы специализированные программные средства «Вентиляция 2.0» и «ВентПЛА», также была выполнена верификация путем проведения натурных измерений. [9].

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

Расчетами установлено, что в ряде ходков и промежуточных штреках внутри блока скорость движения воздуха менее допустимых значений (0,25 м/с) или равна нулю, что подтверждается измерениями в шахте.

Для проветривания одного забоя, по шахтным расчетным данным, необходимо 6,25 м3/с воздуха. Панель, состоящая из блоков (до семи в панели), разделяется на штреки (до семи в блоке) и ходки (до девяти в блоке). Поддерживаемые горные выработки (без перемычек), представленные единой сетью, необходимо проветривать с учетом минимальной скорости движения воздуха. Требуемое количество воздуха составит Qблокблок = n(VS) = (9 + 7) 0,25·21 = 84 м3/c. С учетом одновременной работы четырех забоев расход воздуха будет равен Qпан = 84 + 6,25∙4 = 109 м3/c. В настоящее время воздух в один блок панели подается в меньшем количестве.

Следующим шагом была оценка влияния увеличения общей протяженности горных выработок в панели на общее аэродинамическое сопротивление как основной элемент шахтной вентиляционной сети (ШВС), изменение которого влияет на аэрологическую безопасность. При математическом моделировании проводился расчет величины аэродинамического сопротивления трения (без учета местных и лобовых сопротивлений). При этом принималось, что на каждом этапе полностью отрабатывается один ходок с примыкающими к нему участками внутриблоковых штреков (рис.5). Данные расчетов приведены в таблице с графиком на рис.6.

Рис.5. Порядок компьютерного моделирования воздухораспределения в упрощенной шахтной вентиляционной сети исследуемого объекта: a – 1 этап (2 ходка); б – 2 этап (3 ходка); с – 3 этап ( 4 ходка); г – 4 этап (5 ходков)
Рис.6. Результаты математического моделирования (а) и график зависимости общего аэродинамического сопротивления шахты от общей протяженности горных выработок в панели (б)

Из полученных результатов можно сделать вывод о том, что с увеличением количества выработок общее аэродинамическое сопротивление при данной схеме уменьшается незначительно (в 1,02 раза), при этом существенно (более чем на 1000 м) увеличивается общая протяженность горных выработок. Уменьшение аэродинамического сопротивления происходит по линейной зависимости с весьма близким к нулю угловым коэффициентом (рис.6, б), однако, при использовании существующих специализированных программных средств («Вентиляция 2.0» и «ВентПЛА») не учтены лобовые и местные сопротивления, что не позволяет дать полную оценку влияния применяемой системы разработки на аэрологическую безопасность. Исследования могут быть продолжены только при использовании трехмерного компьютерного моделирования, которое широко применяется в ведущих университетах и научных организациях [2, 10, 15, 16]. Данный вопрос требует проведения дальнейших исследований, а по текущим данным можно сделать вывод, что при увеличении общей протяженности горных выработок в применяемой схеме незначительно увеличивается общее аэродинамическое сопротивление.

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

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

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

2. Установка в ходках быстровозводимых вентиляционных сооружений типа парус (рис.7, б и 8, а), а также изменение технологии подготовки выемочной панели позволит перейти от сложной диагональной к сложной параллельной схеме проветривания участка. Такое решение обеспечит выполнение правил безопасности в части проветривания за счет диффузии тупиков (образованных из-за установки перемычки в сбойках между штреками). Также это позволит сократить количество воздуха, необходимого для проветривания поддерживаемых горных выработок. Внедрение предлагаемой параллельной схемы проветривания позволит исключить вероятность опрокидывания вентиляционной струи воздуха, что, по наблюдениям, происходит в сложной диагональной схеме. Регулирование воздуха между штреками можно проводить с помощью парусных перемычек. Количество воздуха, необходимое для проветривания одного блока (семь штреков, с учетом исключенных из ШВС девяти ходков), составит Q = n(VS) = 7·0,25·21 = 36,75 м3/c (без учета внедрения параллельной схемы проветривания расход составляет 84 м3/с).

Рис.7. Общая аэродинамическая схема участка (а) и схема с предлагаемыми изменениями – использованием легких металлических перемычек (парусов) и измененной последовательностью в проходке конвейерного штрека (б) 1 – свежая струя воздуха; 2 – исходящая струя воздуха; 3 – направление движения проходческого (добычного) забоя; 4 – действующая горная выработка; 5 – планируемая к проходке горная выработка; 6 – предлагаемая легкая металлическая перемычка или парус

3. Установка вентиляторов местного проветривания должна осуществляться на свежей струе воздуха перед первым забоем или перед основным фронтом ведения горных работ (рис.8, б). Это мероприятие позволит не допустить попадания метана в действующие горные выработки из-за возможного опрокидывания вентиляционной струи.

4. Для повышения уровня аэрологической безопасности также необходимо рассматривать комплекс работ по вентиляции и дегазации, что подробно описано в работах [5, 6].

Рис.8. Графическое представление предлагаемых рекомендаций по установке парусов (а) и размещению вентиляторов местного проветривания (б).

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

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

Анализ результатов компьютерного моделирования показал, что изменение сопротивления в любом из определяющих штреков может привести к опрокидыванию вентиляционной струи внутри панели. Расчетами и натурными измерениями установлено, что в ряде ходков внутри панели скорость движения воздуха менее допустимых 0,25 м/с или вообще равна нулю, что приводит к нарушению ФНиП.

Проведенные исследования могут быть использованы при разработке нормативно-правовой базы с учетом критериев «регуляторной гильотины» [7], а также риск-ориентированного подхода при разработке методики прогноза взрывобезопасности выемочных участков угольных шахт [1] и при разработке новых технологий добычи газоносных угольных пластов [13].

Литература

1. Баловцев С.В. К методике прогноза взрывобезопасности выемочных участков угольных шахт // Горный информационно-аналитический бюллетень (научно-технический журнал). 2018. № 11. С. 218-226.

2. Гендлер С.Г. Методика определения параметров системы подогрева воздуха в железнодорожных тоннелях, расположенных в суровых климатических условиях / С.Г.Гендлер, С.В.Синявина // Записки Горного института. 2017. Т. 224. С. 215-222. DOI: 10.18454/PMI.2017.2.215

3. Каледина Н.О. О выборе способа проветривания тупиковых горных выработок газообильных угольных шахт / Н.О.Каледина, С.С.Кобылкин // Горный журнал. 2014. № 12. С. 99-104.

4. Кобылкин С.С. Трехмерное моделирование при проведении инженерных расчетов по тактике горноспасательных работ / С.С.Кобылкин, А.С.Кобылкин // Горный журнал. 2018. № 5. С. 82-85. DOI: 10.17580/gzh.2018.05.13

5. Малашкина В.А. Направления повышения эффективности использования систем дегазации угольных шахт // Горный информационно-аналитический бюллетень (научно-технический журнал). 2019. № 6. С. 206-214. DOI: 10.25018/0236-1493-2019-06-0-206-214

6. Малашкина В.А.Направления повышения эффективности подземной дегазации для улучшения условий труда шахтеров угольных шахт // Горный информационно-аналитический бюллетень (научно-технический журнал). 2018. № 7. С. 69-75. DOI: 10.25018/0236-1493-2018-7-0-69-75

7. Проект постановления Правительства Российской Федерации «О признании утратившими силу нормативных правовых актов и отдельных положений нормативных правовых актов Правительства Российской Федерации, об отмене некоторых актов федеральных органов исполнительной власти и признании не подлежащими применению некоторых нормативных правовых актов, содержащих обязательные требования, соблюдение которых оценивается при проведении мероприятий по контролю при осуществлении федерального государственного надзора в области промышленной безопасности» подготовлен во исполнение пункта 5 плана мероприятий («дорожной карты») по реализации механизма «регуляторной гильотины», утвержденного Председателем Правительства Российской Федерации Д.А.Медведевым 29 мая 2019 г. № 4714п-П36. URL: regulation.gov.ru/projects#npa=94348 (дата обращения: 18.04.2020).

8. Effect of geodynamic setting on spontaneous combustion of coal waste dumps / A.S.Batugin, A.S. Kobylkin, V.R.Musina // Eurasian Mining. 2019. № 2. P. 64-69. DOI: 10.17580/em.2019.02.14

9. Electronic Code of Federal Regulations (e-CFR) Title 30. Mineral Resources. Cornell Law School. URL: www.law.cornell.edu/cfr/text/30 (дата обращения 18.04.2020).

10. Gendler S.G. Evaluation principles of the dust influence of mining enterprises on the / S.G.Gendler, M.L.Rudakov, V.S.Kuznetsov// Latvian Journal of Physics and Technical Sciences. 2019. Vol. 56. № 3. P. 62-69. DOI: 10.2478/lpts-2019-0020

11. Hartman H.L. Mine ventilation and air conditioning. Malabar, Florida: Krieger publishing company, 1991. 791 p.

12. Kaledina N.O. Ventilation of blind roadways in coal mines: Problems and solutions / N.O.Kaledina, S.S.Kobylkin // Eurasian Mining. 2015. № 2. P. 26-30. DOI: 10.17580/em.2015.02.07

13. Kazanin O.I. Numerical study of the air-gas dynamic processes when working out the mosshny seam with longwall faces / O.I.Kazanin, A.A.Sidorenko, Y.G. Sirenko// ARPN Journal of Engineering and Applied Sciences. 2018. Vol. 13. № 4. P. 1534-1538.

14. Kennedy W.R. Practical mine ventilation. Chicago: Interec Publishing Corporation, 1999. 612 p.

15. Semin M.A.Stability of air flows in mine ventilation networks / M.A.Semin, L.Y. Levin // Process Safety and Environmental Protection. 2019. Vol. 24. P. 167-171. DOI:10.1016/j.psep.2019.02.006

16. Simplifying CFD modeling of longwall gobs with a modular meshing approach /R.C.Gilmore, J.A.Marts, J.F.Brune, S.Saki, G.E.Bogin, J.W.Grubb // Mining Engineering. 2015. Vol. 67. Iss. 3. P. 68-72.

17. Validation of the geometrical model and boundary conditions for modeling the process of air intake into the body of a coal waste dump taking into account area geodynamics / A.Batugin, A.Kobylkin, V.Musina, K.Daniil // 18th International Multidisciplinary Scientific GeoConference Surveying Geology and Mining Ecology Management, SGEM 2018, Sofia, Bulgaria. Papers SGEM 2018. 2018. Vol. 18. P. 1111-1118. DOI:10.5593/sgem2018/1.3/S03.140

Похожие статьи

Численное моделирование двустенного шарового резервуара
2020 М. Г. Каравайченко, Л. И. Газалеев
Фальбанды Керетского архипелага Белого моря: характеристика состава пород и минералов, рудная минерализация
2020 Л. И. Салимгараева, С. Г. Скублов, А. В. Березин, О. Л. Галанкина
Оценка смещений рудных контуров после взрыва с применением BMM-системы
2020 Р. А. Рахманов, Д. Лоеб, Н. И. Косухин
Изучение влияния алюминия на адгезию нержавеющей стали при газопламенном напылении
2020 С. Н. Степанов, Т. А. Ларионова, С. С. Степанов
Практика применения аппаратов магнитной обработки для интенсификации процессов первичной подготовки нефти
2020 И. А. Голубев, А. В. Голубев, А. Б. Лаптев
Оценка влияния градиентов водонасыщенности и капиллярного давления на формирование размера зоны двухфазной фильтрации в сжимаемом низкопроницаемом коллекторе
2020 В. А. Коротенко, С. И. Грачёв, Н. П. Кушакова, С. Ф. Мулявин