Пермь, Пермский край, Россия
сотрудник с 01.01.1998 по 01.01.2025
Пермь, Пермский край, Россия
сотрудник
Россия
сотрудник
Россия
ВАК 1.6.10 Геология, поиски и разведка твердых полезных ископаемых, минерагения
ВАК 1.6 Науки о Земле и окружающей среде
ВАК 1.2.1 Искусственный интеллект и машинное обучение
ВАК 1.2.2 Математическое моделирование, численные методы и комплексы программ
УДК 55 Геология. Геологические и геофизические науки
Рассмотрены возможности современных технологий основных методов разведочной геофизики для решения задач, возникающих при эксплуатации калийного рудника. Представлены конкретные примеры применения комплексов геофизических методов для прогноза особенностей геологического строения на новых площадях отработки с целью информационного обеспечения геомеханических расчетов параметров ведения горных работ, мониторинга строения и физических свойств водозащитной толщи над ранее отработанными площадями и на аварийных участках, контроля состояния разного рода элементов подземных горнотехнических конструктивов, изучения устойчивости элементов зданий и сооружений на подработанной территории и оценки экологических рисков.
Калийный рудник, водозащитная толща, сейсморазведка, электроразведка, гравиразведка, геомеханический прогноз, шахтный ствол, межкамерный целик
Введение
Проблема обеспечения продовольствием постоянно растущего населения планеты характеризуется абсолютной актуальностью. Ее решение неразрывно связано с объемами и качеством минеральных удобрений, основным источником которых являются месторождения калийных и магниевых солей.
Россия занимает второе место в мире после Канады по количеству извлекаемых запасов и по объему прогнозных ресурсов калийных солей. В стране сосредоточена почти четверть мировых извлекаемых запасов, а объем ресурсов превышает 15 % мировых. Более 80 % разведанных запасов сосредоточено в крупнейшем в Европе Верхнекамском месторождении калийных и магниевых солей (ВКМКС) в Пермском крае. Разведано и эксплуатируется Гремячинское месторождение в Волгоградской области, открыты Нивенское в Калининградской области, Западно- и Восточно-Петриковское месторождения в Саратовской области, Непское в Иркутской области. Все перечисленные месторождения на территории России, а также за ее пределами, например Саскачеванское в Канаде, Старобинское в Белоруссии, Тюбегатанское в Узбекистане, различаются горно-геологическими условиями разработки. Независимо от этого, главным условием обеспечения безопасной подземной разработки месторождений водорастворимых полезных ископаемых является сохранение полной водонепроницаемости пород водозащитной толщи (ВЗТ), перекрывающей доступ подземных вод в горные выработки. Нарушение сохранности ВЗТ привело к затоплению десятков калийных и каменно-соляных рудников в разных странах мира.
Например, для Саскачеванского месторождения главный камень преткновения в обеспечении безопасной разработки – зоны обрушения (collapse zone ) [1]. Подобные неоднородности большей частью приурочены к краевым частям рифогенных структур девонского возраста [2]. Для Тюбегатанского месторождения в рамках решения той же задачи основную проблему представляют элементы разрывной тектоники [3, 4].
С целью выявления аномальных геологических особенностей ВЗТ и определения степени их опасности для горных работ в основном применяются геофизические методы. Обоснование применимости того или иного метода базируется на физико-геологической модели исследуемого интервала геологического разреза. В условиях пластового залегания контролируемых отложений, которое характерно для большинства калийных месторождений, наибольшее распространение получили сейсморазведочные исследования.
На Саскачеванском месторождении еще в 80-х годах прошлого века достаточно широко внедрялась наземная 2D- и 3D-сейсморазведка [5] в комплексе с шахтными сейсмоакустическими исследованиями [6]. Аналогичный комплекс методов развивается на Гремячинском месторождении [7]. Для Тюбегатанского месторождения уже в этом веке выполнен целый комплекс геофизических методов, включающий наземные: сейсмо-, электро- и гравиразведку – и скважинные сейсмические исследования [4].
На сегодняшний день наиболее широкое применение геофизических методов для решения горно-геологических задач, возникающих на различных стадиях горного производства, отмечается в пределах ВКМКС.
Верхнекамское месторождение эксплуатируется с 30-х годов прошлого столетия, а геофизические работы проводятся на его территории с 1920-х годов. Применялись в основном гравиметрические, магнитометрические, сейсмические (МПВ) и электрометрические исследования. На сегодняшний день продолжается применение сейсмоакустических, гравиметрических и электрометрических методов разведочной геофизики. Подобное содержание геофизического комплекса базируется на физико-геологической модели исследуемого интервала геологического разреза [8] (рис. 1) с учетом современного уровня развития данных методов.
Рисунок 1. Обобщенная физико-геологическая модель для ВКМКС
Метод и теория
В перечень важнейших задач, возникающих в пределах современного эксплуатируемого калийного месторождения, входят:
-
прогноз особенностей геологического строения на новых площадях отработки с целью информационного обеспечения геомеханических расчетов параметров ведения горных работ;
-
мониторинг строения и физических свойств водозащитной толщи над ранее отработанными площадями и на аварийных участках;
-
контроль состояния разного рода элементов подземных горнотехнических конструктивов;
-
изучение устойчивости элементов зданий и сооружений на подработанной территории, оценка экологических рисков.
Для их решения применяются геофизические технологии различной разрешающей способности и информационной составляющей для наблюдений с поверхности и во внутренних точках среды. Основные особенности данных технологий обусловлены как природными, так и техногенными факторами. К природным факторам относятся незначительные глубина залегания и размеры объектов исследований, значительная латеральная изменчивость их строения и свойств, нефтеносность подстилающих отложений. Техногенные включают значительные объемы выработанного пространства, территориальную совмещенность с градопромышленными агломерациями.
1. Прогноз особенностей геологического строения
Здесь ведущее место занимают сейсмо- и электроразведочные методы исследования. Из сейсморазведочных применяется методика многократных перекрытий (ММП) с использованием интерференционных систем регистрации [9]. Параметры системы регистрации обеспечивают уверенную регистрацию отраженных волн в целевом интервале глубин (до 500 м) в частотном диапазоне до 120 Гц, что соответствует требованиям высокоразрешающей сейсморазведки (рис. 2). Цифровая обработка направлена на выделение опорных отражающих горизонтов в интервале глубин от 100 до 400 м: ТКТ – кровля терригенно-карбонатной толщи, СМТ – кровля соляно-мергельной толщи, ПКС – кровля покровной каменной соли или ВС (верхняя соль), Карн. – кровля продуктивной толщи, Сил. – граница между карналлитовой и сильвинитовой пачками, МГ – интервал маркирующей глины, С – подошва подстилающей каменной соли. При этом не допускается их повышенная регуляризация, традиционная при решении типовых структурных задач.
При выполнении площадных (квази-3D) и пространственных (3D) наблюдений информационное обеспечение интерпретационного этапа расширяется за счет возможности анализа пространственных распределений сейсморазведочных параметров. В основном решаются задачи, связанные с локализацией объектов поиска.
Квазипространственные наблюдения организуются при отсутствии прямого доступа к изучаемым участкам: под зданиями и сооружениями, на аварийных участках. Линии приема и возбуждения упругих колебаний разносятся по границам изучаемой территории. Суммирование пространственных выборок по пунктам приема (ОР), пунктам возбуждения (SP) и точкам отражения (ОГТ) позволяет конкретизировать положение аномалии образующих объектов относительно линии точек отражения [10].
Наблюдения 3D выполняются на участках, требующих детальных исследований. Получаемые кубы сейсмической информации характеризуются повышенной пространственной разрешенностью относительно профильных и площадных систем наблюдений, что позволяет прослеживать в амплитудных представлениях геометрические трансформации локальных литолого-структурных неоднородностей ВЗТ (рис. 3).
Рисунок 2. Динамический временной разрез по результатам наземных сейсморазведочных наблюдений
Рисунок 3. Результаты обработки данных малоглубинных исследований 3D: вертикальные срезы куба сейсмических данных меридионального (а) и широтного простирания (б)
Из электроразведочных в основном применяются методы на постоянном токе в комплексе с оценкой изменчивости естественной поляризации (ЕП): вертикальное электрическое зондирование (ВЭЗ), симметричное электропрофилирование (СЭП) и ЕП (рис. 4).
Рисунок 4. Геоэлектрический разрез по данным ВЭЗ
Комплексная интерпретация сейсмо- и электроразведочных данных позволяет локализовать потенциально опасные участки, характеризующиеся пониженными упругими свойствами исследуемых отложений и повышенной интенсивностью гидрогеологических процессов в надсоляной толще, что учитывается при формировании исходной физико-геологической модели для геомеханических расчетов [4] (рис. 5).
Для детализации геологического строения и определения свойств породного массива на уровне продуктивных пластов и прилегающих к ним пластов ВЗТ выполняются сейсморазведочные исследования из горных выработок (рис. 6). Регистрация сейсмических данных во внутренних точках среды обеспечивает разрешающую способность до 1 м в направлении распространения упругих волн. При этом положение нижней границы выделяемых участков определяется наиболее достоверно.
При проведении работ во внутренних точках среды за счет изменения частотного диапазона возможно появление дополнительных ОГ (снизу вверх по разрезу): Гк-, Дк-, Ек-, Ик-кровля карналлитовых пластов Г, Д, Е, И; ПКСп-, ПКСк-подошва и кровля покровной каменной соли. Целевой интервал при наблюдениях в верхнее полупространство, как правило, не превышает 150 м [11] (рис. 6).
Рисунок 6. Динамический временной разрез по результатам шахтных сейсморазведочных наблюдений
В случае реализации наземно-подземного комплекса сейсмических исследований выполняется построение разноуровневых моделей реальных сред, согласованных по эффективным кинематическим параметрам. За счет разрешающей способности шахтных сейсмических исследований повышается точность вертикальной локализации геологических неоднородностей, определения их строения и свойств. Физико-геологические модели повышенной детальности, построенные по результатам комплекса сейсморазведочных исследований, используются в текущих и в прогнозных геомеханических оценках состояния ВЗТ и на потенциально опасных участках шахтных полей действующих калийных рудников [12].
В случае отсутствия горнотехнических условий проведения сейсмоакустических наблюдений из горных выработок повышение разрешающей способности интерпретационных заключений возможно за счет включения в комплекс сейсморазведочного контроля горнотехнических ситуаций блока скважинных модификаций. Содержание данного блока определяется необходимой разрешающей способностью и качеством параметрического обеспечения, как на поверхности, так и во внутренних точках среды [13].
С целью определения скоростных параметров породного массива наиболее информативными для интервала малых глубин зарекомендовали себя продольное, непродольное и обращенное вертикальное сейсмопрофилирование (ВСП), а также межскважинное просвечивание [14, 15].
В целях картирования малоразмерных природно-техногенных объектов предложено скважинное сейсмическое профилирование в рамках методики многократных перекрытий (рис. 7) [16]. Предполагается наличие вертикально ориентированных отражающих границ, удаление которых от скважины оценивается в горизонтальном направлении.
Рисунок 7. Динамический временной разрез по результатам шахтных сейсморазведочных наблюдений. Результаты обработки скважинного ОГТ, разрезы: а) временной, б) динамический, в) скоростной, г) частотный
По гравиметрическим данным выявляются зоны распространения и вероятный интервал глубин разуплотнения горных пород. Например, на рис. 8 приведено вертикальное сечение гравитационного поля через скважины 1 и 2. В районе скважин наблюдаются две изолированные отрицательные аномалии на эффективных глубинах порядка 100 м, отражающие зоны разуплотнения в надсоляной толще.
Рисунок 8. Вертикальное сечение трехмерной трансформанты гравитационного поля
2. Мониторинг строения и физических свойств
При мониторинге территории действующих шахтных полей с целью определения участков, требующих повышенной детальности сейсмо- и электроразведочных исследований, выполняются гравиметрические наблюдения. Среднеквадратическая погрешность определения силы тяжести не должна превышать ± 0,02 мГал. Основная задача – картирование зон разуплотнения в исследуемом интервале породного массива. Количественная оценка глубины залегания источников локальных аномалий определяется на основе решения трехмерной обратной задачи гравиразведки с использованием монтажных технологий (рис. 9).
Рисунок 9. Пространственное распределение функции локализации источников локальных гравитационных аномалий
Сочетание полученных результатов в пределах одной территории работ позволяет конкретизировать природу интерпретационных заключений. Например, по совокупности негативных изменений геофизических параметров в пределах горного отвода калийного месторождения в рамках геофизического комплекса, включающего сейсмо-, электро- и гравиразведку, выделено три участка (рис. 10): 1) на северо-востоке (зоны 1 и 2 сейсморазведочных аномалий) характеризуется пониженными значениями упругих параметров – скорости и плотности, а также удельного электрического сопротивления; 2) в центре (зона 4 сейсморазведочных аномалий) аналогичные признаки по негативным изменениям всех трех составляющих; 3) на юго-западе площади работ отмечается снижение значений упругих параметров при противоречивом поведении электрометрических параметров от незначительного понижения к преобладающему повышению значений.
Подобные сочетания рассматриваемых геофизических параметров позволяют предположить понижение прочностных свойств в целом для надсоляной и соляной толщ на всех трех участках, при этом на первых двух имеет место реализация негативных гидрогеологических процессов.
Рисунок 10. Результаты комплексной интерпретации сейсмо-, электро- и гравиразведочных исследований
В пределах действующих шахтных полей на месторождениях водорастворимых полезных ископаемых возникновение так называемых аварийных участков не является чем-то из ряда вон выходящим. На сегодня известно около 100 рудников, вышедших из строя вследствие прорыва в выработанное пространство надсоляных вод. В пределах ВКМКС таковых два: на Балахонцевском (1986) и Березниковском (2006) участках. При эксплуатации подобных месторождений вероятно появление участков с потенциальными рисками нарушения сплошности водозащитной толщи (ВЗТ) за счет недоучета локальных особенностей геологического строения и латеральной изменчивости физико-механических свойств породного массива.
Основной критерий выделения – увеличение скоростей оседания, возможна микросейсмическая активность. На первом этапе их изучения выполняются: визуальное обследование горных выработок, маркшейдерские и геофизические наблюдения с последующей геомеханической оценкой [12].
Далее по результатам геомеханических расчетов определяется необходимость организации геофизических мониторинговых наблюдений, их содержание и периодичность. Периодическое получение информации об изменении упругих свойств позволяет спроектировать и реализовать необходимые горнотехнические меры для поддержания устойчивости породного массива. Данная возможность крайне важна в аварийных ситуациях для оценки эффективности мер, принимаемых по их купированию (рис. 11).
Рисунок 11. Инъекция тампонажного раствора в зону максимальных негативных изменений значений упругих параметров
Оценка эффективности подобной меры выполняется по изменчивости значений интервальных скоростей распространения упругих волн в зоне воздействия (рис. 12). В приведенном примере на забое инъекционной скважины в интервале ОГ ВС – Ек отмечается рост значений интервальной скорости до 200 м/с. Подобное изменение упругих параметров позволяет предположить укрепление дезинтегрированной зоны породного массива.
Рисунок 12. Сейсмогеологические разрезы: а) сентябрь, б) ноябрь, в) разность разрезов
В случае необходимости организации непрерывного онлайн-мониторинга упругих свойств породного массива применяется технология распределенных оптоволоконных акустических датчиков [17].
На одном из аварийных участков соляного рудника смонтирована оптоволоконная система, представляющая собой единую оптическую линию, позволяющую регистрировать разнотипные прямые и преломленные волны одновременно во всех скважинах. Достоверность информации подтверждается сравнением данных, зарегистрированных традиционными датчиками и оптоволоконной системой (рис. 13).
Основным недостатком системы DAS на данный момент является низкое отношение сигнал/шум. В этой связи в условиях мониторинга с постоянным источником упругих колебаний применяется повышенное количество возбуждений [18].
Рисунок 13. Сравнение волновых картин, регистрируемых гидрофонами (а, б) и с помощью распределенной оптоволоконной системы (в, г)
3. Контроль состояния
Геометрия и расположение горных выработок в шахтном поле ограничивает пространственную приуроченность получаемой информации о строении и свойствах породного массива. С другой стороны, за счет особенностей конфигурации выработанного пространства возможны исследования с различной пространственной ориентацией, комплексирование которых формирует объемные представления об особенностях строения и свойствах породного массива в окрестностях горной выработки.
На большинстве калийных месторождений при подземной разработке применяется камерно-столбовая система. Контроль состояния обеспечивающих ее устойчивость горнотехнических элементов – обязательная и необходимая мера. Известны методические решения по прогнозу остаточного срока службы межкамерных целиков, основанные на информации об ослаблении прочностных свойств пород в приконтурной части целика за счет применения сейсмоакустического просвечивания на отраженных волнах [19] (рис. 14).
Рисунок 14. Временной разрез межкамерного целика: I – граница ближней зоны влияния выработки, II – граница дальней зоны влияния выработки, III – дальняя стенка выработки
Освоение подземного пространства любого калийного рудника начинается с проходки шахтных стволов. Наиболее распространен способ заморозки, направленный на купирование пересекаемых водоносных пластов. При его реализации обязателен мониторинг формирования ледопородного массива, где применяются сейсмоакустические исследования в рамках методических решений скважинного ОГТ [16]. При последующей эксплуатации ствола за счет прекращения замораживающих процедур вероятно появление локальных водопритоков в затюбинговом пространстве. Здесь для изучения малоразмерных объектов (< 1 м) породного массива, залегающих на незначительной удаленности от выработанного пространства (< 10 м), по очевидным причинам целесообразно проводить исследования в акустическом диапазоне частот (> 3 кГц) [11].
Возможности применения сейсморазведочных подходов в акустическом диапазоне частот, как и мониторинге ледопородного массива, обуславливаются наличием вертикально ориентированных отражающих границ – тюбинг, цементное кольцо, породный массив. Данная модель характеризуется изменчивостью упругих и геометрических параметров основных структурных элементов. Подобные исследования реализованы на шахтных стволах различных рудников. На рис. 15 представлены результаты одного из таких исследований. Работы проводились с целью выявления зон повышенной циркуляции воды за стволовой крепью.
Рисунок 15. Результаты цифровой обработки акустических данных ОГТ в диапазонах частот: а) 10–20 кГц; б) 3–10 кГц
Обработка проведена отдельно в двух неперекрывающихся частотных диапазонах: 1) 3–10 кГц, 2) 10–20 кГц. Разделение волнового поля на условно низкочастотную и высокочастотную составляющие акустических колебаний позволило повысить разрешающую способность в ближайшем интервале исследований (до 3–5 м) и увеличить глубинность наблюдений до 10 м.
4. Изучение устойчивости
Территории шахтных полей часто совмещены с урбанизированными территориями, что представляет потенциальную угрозу для их инфраструктуры. Для подобных территорий сформированы две системы мониторинга породного массива. Первая система мониторинга, представленная в разделе 2, направлена на контроль изменений строения и свойств разрабатываемого интервала глубин породного массива, отвечающего за устойчивость подземных и поверхностных сооружений. Подобная организация системы мониторинга и легла в основу контроля строения и свойств верхней части разреза (ВЧР) урбанизированных территорий над горными выработками. Геофизические исследования представлены в основном наземными и скважинными сейсморазведочными исследованиями, как ВЗТ, так и ВЧР в комплексе с наземной электроразведкой на постоянном токе [20].
Детальность, периодичность наблюдений и необходимость дополнительных видов мониторинга определяется с учетом получаемых результатов. По результатам геолого-геофизических исследований выполняется количественная геомеханическая оценка деформаций породного массива с учетом оседаний земной поверхности. На ее основе осуществляется прогноз их воздействия на здания и сооружения с рекомендациями по конструктивным мерам защиты и организации онлайн-мониторинга зданий и упругих свойств верхней части разреза.
Особое значение в пределах подобных территорий приобретают вопросы контроля и экологической обстановки, влияющей на повседневную жизнедеятельность населения. Калийное производство, как и любое другое, связано с появлением определенного вида отходов, которые утилизируются в рамках согласованных полигонов. Наиболее перспективным способом экологической диагностики загрязнения подземных вод территорий, подвергающихся техногенным воздействиям от солеотвалов и шламохранилищ, является применение электроразведки методом сопротивлений – СЭП, ВЭЗ, электротомография. На основе анализа данных качественной и количественной интерпретации определяется пространственная закономерность аномальных участков пониженных значений электрических сопротивлений, связанных с повышенной минерализацией подземных вод.
Рисунок 16. Результаты электроразведки при выделении зоны пониженных сопротивлений
Выводы
Применение современных геофизических технологий позволяет снизить горно-геологические риски, возникающие при эксплуатации любого калийного месторождения. Результаты их применения могут использоваться как напрямую при прогнозировании разного рода деструктивных процессов породного массива в интервале техногенного воздействия, так и для параметрического обеспечения геомеханических оценок степени данного воздействия.
1. Халабура С. П. (Стив) Обзор геологии добычи калийных солей методом растворения в Саскачеване / Халабура С. П. (Стив), Харди М. П. // Solution Mining Research Institute: Technical Meeting. – 2007
2. Гендзвилл Д. Дж. Характеристика скального массива вокруг калийного рудника в Саскачеване с использованием геофизических методов: обзор / Gendzwill D. J., Stead D. // Canadian Geotechnical Journal. – Т. 29. – № 4. – 1992. – С. 666–674.
3. Гагала Л. Таджикский бассейн и Юго-Западный Тянь-Шань, Северо-Западная Индо-Азиатская зона столкновения: 1. Структура, кинематика и соляная тектоника в Таджикском складчато-надвиговом поясе Западного предгорья Памира / Гагала Л., Ратшбахер Л., Рингенбах Дж.-К., Куфнер С.-К., Шурр Б., Дедов Р. [и др.] // Tectonics. – 2020.
4. Барях А. А. Информационное обеспечение геомеханических расчетов устойчивости подработанного породного массива со сложным тектоническим строением / Барях А. А., Санфиров И. А., Дьяконов М. В., Лобанов С. Ю., Никифорова А. И. // Russian Journal of Earth Sciences. – 2024. – Т. 24. – ES1012. – DOI:https://doi.org/10.2205/2024es000894.
5. Гендзвилл Д. Дж. Холм Виннипегосис и прерийные эвапориты Саскачевана – сейсмическое исследование / Гендзвилл Д. Дж. // Бюллетень Американской ассоциации геологов-нефтяников. – 1978. – Т. 62. – С. 73–86
6. Гендзвилл Д. Дж., Рэнди Б. Высокоразрешающие сейсмические отражения в калийном руднике / Гендзвилл Д. Дж., Рэнди Б. // Geophysics. – 1993. – Т. 58. – № 5. – С. 741–748.
7. Кулагов Е. В. Сейсморазведка методом отраженных волн в подземных выработках Гремячинского месторождения калийных солей / Кулагов Е. В., Хвесеня С. С., Иванов С. В., Богдан С. И. // Инженерная и рудная геофизика – 2023: сборник материалов конференции. – С-Петербург, 2023. – С. 4.
8. Глебов С. В. Геофизическое обеспечение разработки Верхнекамского месторождения солей / Глебов С. В. // Горный информационно-аналитический бюллетень. – 2004. – №9. – С. 89–92.
9. Санфиров И. А. Рудничные задачи сейсморазведки МОГТ / Санфиров И. А. – Екатеринбург: УрО РАН, 1996.
10. Санфиров И. А. Развитие интерпретационных возможностей сейсморазведочных исследований при освоении месторождений водорастворимых полезных ископаемых / Санфиров И. А., Жикин А. А., Богданов Р. А., Фатькин К. Б. // Горный журнал. – 2021. – №4. – С. 28–33.
11. Санфиров И. А. Комплексные методические решения в шахтной сейсмоакустике / Санфиров И. А., Бабкин А. И., Ярославцев А. Г. // Геофизика. – 2014. – № 5. – С. 10–15.
12. Барях А. А. Сейсмо-геомеханический прогноз состояния водозащитной толщи на калийных рудниках / Барях А. А., Санфиров И. А., Федосеев А. К., Бабкин А. И., Цаюков А. А. // Физико-технические проблемы разработки полезных ископаемых. – 2017. – № 6. – С. 10–22.
13. Санфиров И. А. Горнотехнические приложения малоглубинной скважинной сейсморазведки / Санфиров И. А., Чугаев А. В., Бабкин А. И., Лисин В. П., Бобров В. Ю. // Геофизика. – 2018. – № 5. – С. 24–30.
14. Лисин В. П. Исследование потенциально опасных зон ускоренного оседания грунтовой толщи шахтных полей методами наземной и скважиной сейсморазведки / Лисин В. П., Чугаев А. В., Санфиров И. А. // Инженерная геофизика 2018 (EAGE): сборник материалов конференции. – Алма-Ата, 2018.
15. Чугаев А. В. Определение скоростных параметров разреза по сейсмическим скважинным исследованиям / Чугаев А. В., Лисин В. П., Санфиров И. А., Никифоров В. В., Зеленин В. П. // Инженерная геофизика 2017 (EAGE): сборник материалов конференции. – Кисловодск, 2017.
16. Патент № 2706910 C1 Российская Федерация, МПК G01V 1/00, G01V 1/28, G01V 1/30. Способ контроля толщины ледопородного ограждения при строительстве шахтных стволов: № 2019109327: заявл. 29.03.2019: опубл. 21.11.2019 / Санфиров И. А. , Ярославцев А. Г., Бабкин А. И., Чугаев А. В.
17. Санфиров И. А. Направления развития сейсморазведочного мониторинга водозащитной толщи действующих соляных рудников / Санфиров И. А., Чугаев А. В., Трапезникова А. Б., Тезиков А. Д. // Горный журнал. – 2023. – № 11. – С. 27–31.
18. Патент № 2809469 C1 Российская Федерация, МПК G01V 1/42, G01D 5/26. Способ и система сейсмоакустического контроля массива горных пород: № 2023109620: заявл. 14.04.2023: опубл. 12.12.2023 / Чугаев А. В., Санфиров И. А., Кузнецов А. И., Богданов Р. А.
19. Барях А. Инструментальный контроль состояния соляных целиков для прогнозирования их остаточного ресурса работы / Барях А., Санфиров И., Асанов В., Бабкин А., Гегин А., Бруев А. // Горная геотехника и подземное строительство в начале XXI века (XXIV зимняя школа механики горных пород) 12–16 марта 2001 г. – Вроцлав, 2001. – С. 15–24.
20. Санфиров И. А. Малоглубинные геофизические исследования на Верхнекамском месторождении калийных солей / Санфиров И. А., Степанов Ю. И., Фатькин К. Б., Герасимова И. Ю., Никифорова А. И. // Физико-технические проблемы разработки полезных ископаемых. – 2013. – № 6. – С. 71–77.
