Magnitudes of seismic events induced by fluid injections in the Earth’s crust

Full Text

ВВЕДЕНИЕ

Сейсмичность индуцируется не только при строительстве плотин [Rastogi et al., 1998; Lamontagne et al., 2006], подземной разработке месторождений [McGarr, 1992; Shen et al., 2008; Sato et al., 1988], но и при закачке воды для извлечения углеводородов из низкопроницаемых резервуаров, что более подробно изложено в последующих разделах, а также в искусственных геотермальных системах (Enhanced Geothermal Systems, EGS) [Ladner et al., 2009; Cladouhos et al., 2010; Nicolas et al., 2011; Convertito et al., 2012]. Последнее явление известно как тригерная сейсмичность или сейсмичность, вызванная ингжекцией флюидов (Injected Induced Seismisity, IIS). IIS возникает, когда напряженное состояние разлома приближается к критерию разрушения из-за снижения эффективного стресса, вызванного увеличением порового давления [Hubbert et al., 1959; Zoback, 2010]. Сейсмичность при небольших землетрясениях не вызывает проблем, в то время как более сильные землетрясения могут вызвать значительные разрушения.

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

Сверхкритический CO₂ закачивается в подземный массив горных пород в рамках проектов по захоронению CO₂ (Carbon capture and storage, CCS). Используя полученные уравнения, оценивается IIS, вызванный CCS.

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

В статье учтены также последние экспериментальные данные по индуцированной сейсмичности, опубликованные в работах [Foulger et al., 2018; Grigoli et al., 2018; Schultz et al., 2020; Kocharyan, 2021; Li et al., 2022].

ВЗАИМОСВЯЗЬ МЕЖДУ ОБЪЕМОМ ЗАКАЧИВАЕМОЙ ВОДЫ И МАГНИТУДОЙ ЗЕМЛЕТРЯСЕНИЙ

Для обоснования уравнений представления верхней границы максимальной магнитуды триггерной сейсмичности рассмотрены фактические данные триггерной сейсмичности, представленные в табл. 1. Они были выбраны таким образом, чтобы учитывался объем закачанной воды и максимальная магнитуда землетрясений. Зависимость между логарифмом объема закачанной воды, V (в м³), и максимальной магнитудой вызванного сейсмического события, M_max, изображена графически, как показано на рис. 1.

 

Рис. 1. Теоретические зависимости и экспериментальные данные, связывающие объем закачиваемого флюида V (м3) и максимальную магнитуду Mw(max) триггерного сейсмического события. Легенда: A-I – экспериментальные данные по инжекции воды (см. табл. 1); 7 – теоретическая зависимость для воды, уравнение (7); J-K – экспериментальные данные по инжекции scCO2 (см. табл. 1); 8 – теоретическая зависимость для scCO2, уравнение (8); N-R ‒ данные по объемам извержений вулканов (см. табл. 1); 9 ‒ теоретическая зависимость для магмы, уравнение (9); Model – результаты CFRAC-моделирования инжекции магмы в сдвиговую трещину в основании Мутновского вулкана [Kiryukhin et al., 2022, вариант #10, Table 3].

 

Авторы не ставят задачу провести всеобъемлющий обзор IIS. Ниже приводятся лишь краткие пояснения к конкретным сценариям закачки. Максимальная магнитуда 4.7 была зарегистрирована в Арканзасе в результате закачки воды при добыче нефти и природного газа, когда сточные воды от гидроразрыва пласта закачивались для повышения нефтеотдачи (случай A [Horton et al., 2012]). Загрязненная вода с военного завода закачивалась в подземный резервуар в Денвере, штат Колорадо, с 1962 по 1967 гг., и в результате была зарегистрирована сейсмичность с максимальной магнитудой 5.3 (случай B [Healy, 1968; Zoback, 2010]). Триггерная сейсмичность зарегистрирована при закачке воды в EGS для извлечения геотермальной энергии (случаи C и D [Cladouhos, 2010; Ladner, 2009]). В качестве примера можно привести геотермальное месторождение Гейзерс в Калифорнии, где была зарегистрирована максимальная магнитуда 4.4 (случай С [Cladouhos, 2010]). В уникальном случае рассол из купола каменной соли был закачан в скважину, чтобы избежать загрязнения реки; сейсмическая активность магнитудой 4.4 была зарегистрирована в Парадокс Вэлли, Калифорния (случай J [Block et al., 2014]). Известны случаи, когда вода закачивалась на глубину до 9 км (случай I, немецкая Kontinentales Tiefbohrprogramm der Bundesrepublik Deutschland (KTB) [Zoback et al., 1997]) или даже в разлом Нодзима после землетрясения в Кобе в 1995 г. (случай G [Nishigami et al., 2002]) в научных целях.

A. Nicol с соавторами [Nicol et al., 2011] рассмотрели сценарии закачки воды в горячую сухую породу и для добычи углеводородов и вывели следующее уравнение:

$M_{\max }=0.77\cdot Ig\left(V\right)-0.51$ (1)

G. Zöller и M. Holschneider [Zöller, Holschneider, 2014] и J.H. Dieterich с соавторами [Dieterich et al., 2015] оценили теоретическую максимальную величину магнитуды триггерной сейсмичности следующим образом:

$M_{\max }=0.67\cdot Ig\left(V\right)$ (2)

A. McGarr [2014] также теоретически вывел уравнения, но не для максимальной магнитуды, а для верхней границы сейсмического момента M₀ триггерного землетрясения в виде:

$M_0=G\cdot V$ (3)

где G ‒ модуль сдвига, и подтвердил, что уравнение (3) может представлять верхний предел для сейсмического момента сильнейшего триггерного землетрясения в 18-ти случаях IIS с модулем сдвига 30 ГПа. Используя взаимосвязь между сейсмическим моментом M₀ и величиной M:

$\lg \left(M_o\right)=9.05+1.5\cdot M$ (4)

и принимая модуль сдвига равным 30 ГПа, коэффициент трения μ = 0.6, параметр в законе Гутенберга-Рихтера b = 1, уравнение (3) можно преобразовать к следующему виду:

$M_{\max }=0.67\cdot Ig\left(V\right)+0.95$ (5)

где M_max ‒ верхняя граница максимальной величины триггерной сейсмичности (см. рис. 1).

Уравнения (1) и (2) можно обобщить в виде уравнения:

$M_{\max }=A\cdot Ig\left(V\right)+B$ (6)

где A и B ‒ константы, определяющие верхнюю границу максимальной величины IIS, будет предпочтительнее, чем использование уравнения, подобного уравнению (1) или (2), для представления максимальной величины IIS.

В данной работе коэффициент B скорректирован, поскольку случаи B и H немного превышают оценки по уравнению (5). Однако A был оставлен неизменным, поскольку все три теоретических исследования показали одно и то же значение; B принят равным 1.42 для удовлетворения условий В и Н (см. рис. 1). В итоге уравнение:

$M_{\max }=0.67\cdot Ig\left(V\right)+1.42$ (7)

предлагается использовать для оценки верхней границы M_max по IIS на основе случаев от A до J в табл. 1. Это уравнение соответствует относительно высокому модулю сдвига 151 ГПа, имеющему “эффективный” смысл при принятых выше значениях коэффициента трения, параметра b (в законе Гутенберга-Рихтера), геологических и других неопределенностях. Сейсмическая энергия максимальных сейсмических событий составляет приблизительно от 1/10 до 1/10000 энергии закачки, которая может быть рассчитана как произведение устьевого давления и объема закачиваемой воды. По фактическим данным (см. табл. 1) закачка воды не приводила к высвобождению энергии, превышающей энергию закачки.

 

Таблица 1. Примеры триггерной сейсмичности, вызванные инжекцией воды (A-J), сверхкритического CO₂ (K-M) и магмы (N-R)

	Резервуар	Флюид	Глубина инжекции, м	Объем инжекции, млн м3	Устьевое давление, MPa	Мmax	Ссылки
A	Well 1& 5, North Central AK, USA (2010)	Сточные воды	1800	0.63	15.7	4.7	[Horton et al., 2012]
B	Denver, USA (1962‒1967)	Загрязненные сточные воды	3700	0.62	11.4	5.3	[Healy et al., 1960; Zoback, 2010]
C	The Geysers, USA (1972‒2006)	Вода для EGS	6500‒8500	850.0	0.53	4.4	[Cladouhos, 2010]
D	Basel 1, Switzerland (2006)	Вода для EGS	5000	0.012	30	3.4	[Ladner, 2009]
E	GPK3, Soultz, France (1987‒2007)	Вода для EGS	5000	0.205	7	2.2	[Nicolas, 2011]
F	Nojima fault, Japan (1997)	Вода, эксперимент	540	0.000069	4.6	0.6	[Nishigami et al., 2002]
G	Nojima fault, Japan (2000)	Вода, эксперимент	540	0.000131	4.6	0.3	[Nishigami et al., 2002]
H	Matsushiro, Japan (1970)	Вода, эксперимент	1800	0.000302	5	2.8	[Nishigami et al., 2002]
I	KTB, Germany (1994)	Раствор калия	9030	0.000200	53	1.2	[Zoback, Harjes, 1997]
J	Paradox Valley, US (2013)	Рассол	4300‒4800	9.5	81	4.4	[Block et al., 2014]
K	Decatur, IL, USA (2011‒2014)	CO2 для CCS	2100	2.1	–	1.3	[Kayen et al., 2015]
L	Sleipner, Norway (1996)	CO2 для CCS	800‒1100	28.8	–	4.7	http://www.secarbon.org/wp-content/uploads/2011/05/Hagen.pdf
M	Cogdell, TX, USA (2004‒2012)	CO2 для CCS	<2100	37.2	–	4.4	[Gan, Frohlich, 2013]
N	Барбарунга, Исландия (2014)	магма	–	480.0–510.0	–	> 5	[Sigmundsson et al., 2015]
O	П. Толбачик (2012)	магма	–	550.0	–	4.9	[Belousov et al., 2015], КФ ФИЦ ЕГС РАН
P	Ключевской (2000‒2021)	магма		587.0		4.6	[Coppola et al., 2021] КФ ФИЦ ЕГС РАН
R	Безымянный (2000‒2021)	магма		102.0		3.9	[Coppola et al., 2021] КФ ФИЦ ЕГС РАН

 

СЕЙСМИЧНОСТЬ, ВЫЗВАННАЯ ДОБЫЧЕЙ СЛАНЦЕВОГО ГАЗА

Общемировое количество воды, добываемой при традиционной добыче нефти и газа, составляет приблизительно 300 миллионов баррелей/год (веб-сайт JOGMEC, по данным от 15 октября 2015 г. в Японии). Общемировое производство природного газа составляет 2870 Гм³/год (2006 г.) или 19 млрд баррелей/год нефтяного эквивалента. Мировая добыча сырой нефти составляет приблизительно 30 млрд баррелей/год. Общая добыча составляет 49 млрд баррелей/год, а количество добываемого рассола составляет примерно 0.6% от добычи нефти и газа. Рассол обрабатывается и закачивается в подземные резервуары.

Для добычи сланцевого газа требуется от 8000 до 25000 м³ воды для гидроразрыва пласта на скважину (по данным Массачусетского технологического института, 2011 г., http://mitei.mit.edu/publications/reports-studies/future-natural-gas). Например, количество используемой воды составляет 7.1 млн м³/год в Барнетте (Техас), 5.1 млн м³/год в Файетвилле (Арканзас), 2.7 млн м³/год в Хейнесвилле (Лос-Анджелес), и до 13.5 млн м³/год в Марселлусе. Из воды, полученной при гидроразрыве пласта, 15‒25% уходит обратно в течение 30 дней. Эта вода обрабатывается и закачивается в подземные резервуары. В 2009 г. объем добычи сланцевого газа составил около 85 млн м³ (нефтяной эквивалент). Объем закачиваемой воды составляет примерно 33% от объема добычи газа только при гидроразрыве пласта на пяти основных месторождениях сланцевого газа или, по крайней мере, в 50 раз больше, чем объем закачиваемого рассола при добыче обычных нефти и газа. Следовательно, при добыче сланцевого газа ожидается гораздо более сильная сейсмичность.

Землетрясения М 2.7 и М 1.5 произошли 1 апреля 2011 г. и 27 мая 2011 г. соответственно в 300 м от зоны гидроразрыва пласта Preese Hall Well PH1 компании Cuadrilla Resources в Великобритании (веб-сайт Британской геологической службы http://earthquakes.bgs.ac.uk/research/earthquake_hazard_shale_gas.html по данным от 20 октября 2015 г.). В Оклахоме триггерная сейсмичность с максимальной магнитудой 2.8 возникла в пределах 4 км от зоны гидроразрыва пласта на скважине 4-18 блока B Пикет 17‒18 ноября 2011 г. [Holland, 2013]. A.L. Lleons и A.J. Michael [Lleons, Michael, 2012] показали, что количество землетрясений с магнитудой более 3 в Оклахоме и Арканзасе внезапно увеличилось в 2009 г. Зарегистрировано внезапное увеличение частоты сейсмических событий вокруг сланцевого месторождения Файетвилл в Арканзасе. Это увеличение и событие M 5.7 в штате Оклахома в ноябре 2011 года [Keranen et al., 2012] считаются искусственно вызванными явлениями. Случайная нулевая гипотеза, принятая для последовательности возникновения землетрясений, была отвергнута с p-значением менее 0.1%. [Ellsworth et al., 2012] получили аналогичные результаты. Далее, [Shirezaei et al., 2016] показали наблюдаемый подъем поверхности земли в результате закачки сточных вод в восточном Техасе, их модель позволила сделать вывод, что увеличение порового давления более чем на 1 МПа может вызвать землетрясения, включая самое сильное землетрясение M 4.8 в восточном Техасе 17 мая 2012 г.

Эти наблюдения и анализ показывают, что триггерная сейсмичность может возникать при закачке воды, связанной с добычей сланцевого газа. Приведенное выше уравнение (7) может быть использовано для количественного анализа IIS в контексте оценки возможного объема воды при заданном предельном значении M_max.

Подставляя количество воды для основных месторождений сланцевого газа в уравнение (7), ожидается максимальное значение М 5.7‒6.2 при добыче в течение 1 года на каждом участке. Максимальная величина составит М 6.7‒7.2 при 30-летней добыче. Учитывая, что рост М_max в зависимости от времени довольно неравномерен, а поровое давление уменьшается со временем за счет диффузии, как указано в работе [McGarr, 2014], оценка для 30-летней добычи может быть завышенной. Однако, учитывая, что в 2011 г. уже произошло землетрясение M 5.7 вблизи Праги, Oклахома [Keranen et al., 2012], а затем (3 сентября 2016 г.) произошло землетрясение M 5.8 в 15 км к северо-западу от Пауни, Oклахома (веб-сайт Геологической службы США в поиске архивов землетрясений http://earthquake.usgs.gov/earthquakes/search/ 8 октября 2016 г.), вышеуказанная верхняя граница близка к реальности.

ТРИГГЕРНАЯ СЕЙСМИЧНОСТЬ ПРИ CCS (СARBON СAPTURE AND STORAGE)

При эксплуатации пилотных установок CCS для захоронения суперкритического СО₂ сообщалось о случаях триггерной сейсмичности, включая случай в Декатуре, штат Иллинойс, с М_max = 1.3 [Kaven et al., 2016]. Также сообщалось о случае после закачки CO₂ для повышения нефтеотдачи в Когделле, штат Техас, с М_max = 4.4 [Gan, Frohlich, 2013]. Сейсмические события в пределах 200 км от нагнетательной скважины до и после закачки CO₂ в известном проекте CCS в Слейпнере, Норвегия, были исследованы с помощью поиска архивов землетрясений на сайте USGS (веб-сайт Геологической службы США, http://earthquake.usgs.gov/earthquakes/eventpage/nc72282711#general_summary 20 октября 2015 г.) и проанализированы в зависимости от расхода закачки CO₂ (веб-сайт Statoil http://www.secarbon.org/wp-content/uploads/2011/05/Hagen.pdf 20 октября 2015 г.). Количество сейсмических событий уменьшилось с началом закачки. Однако после 2000 г. сейсмичность активизировалась, и было отмечено несколько более значительных сейсмических событий, одно из которых достигло величины M 4.7. Сейсмичность прекратилась через несколько лет после окончания закачки. Одним из механизмов снижения сейсмичности может быть восстановление порового давления после добычи нефти и газа.

Хотя в работе [Nishimoto et al., 2007] показано, что влияние порового сверхкритического CO₂ на механические свойства породы при трехосном сжатии почти такое же, как влияние поровой воды, данные из рассмотренных выше случаев K-M позиционируются ниже уравнения (7) (см. рис. 1). Поэтому мы уменьшили константу B для землетрясений, вызванных CCS, чтобы оценить верхнюю границу максимальной магнитуды (см. рис. 1):

$M_{\max }=0.67\cdot Ig\left(V\right)-0.30$ (8)

Сравнивая уравнение (8) с уравнением (7), верхняя граница M_max при том же объеме закачки CO₂ меньше на 1.72, чем при закачке воды. Опыта достаточного количества случаев CCS нет, поэтому землетрясения, превышающие предельное уравнение (8), могут произойти в будущем. На данном этапе изученности можно предположить, что разница в сжимаемости между водой и сверхкритическим CO₂ может вызвать эти отклонения.

Выбросы CO₂ в Японии составляют приблизительно 1.2 Гт/год. Это количество можно преобразовать в объем сверхкритического CO₂ в 2.56 Гм³/год, исходя из плотности 469 кг/м³. Глубина закачки должна составлять не менее 740 м, исходя из сверхкритических условий (давление > 7.4 МПа и температура > 304 К), и для упрощения расчетов была принята равной 800 м.

Если предположить, что площадка CCS будет эксплуатироваться в течение 30 лет, то объем закачки должен составлять менее 62000 м³/год для площадки CCS при прогнозном значении М 6.5 (формула 8). Этот объем составляет всего 1/41000 выбросов CO₂ в Японии. Целевое значение сокращения выбросов CO₂ с помощью CCS, предложенное Японской межправительственной группой экспертов по изменению климата (IPCC), составляет 17% в 2050 г. Следовательно, для закачки 17% CO₂ потребуется 7000 объектов CCS при этом прогнозное значение М_max составит М 6.5, землетрясения такой магнитуды происходят в Японии приблизительно раз в три года и наносят существенный ущерб.

Как полагают некоторые исследователи [Cappa, Rutqvist, 2011, 2012; Zoback, 2013; Zoback, Gorelik, 2012], триггерная сейсмичность может уменьшиться в результате утечек CO₂.

ТРИГГЕРНАЯ СЕЙСМИЧНОСТЬ ПРИ МАГМАТИЧЕСКОЙ АКТИВНОСТИ, ПРЕДШЕСТВУЮЩЕЙ ИЗВЕРЖЕНИЯМ ВУЛКАНОВ

Триггерная сейсмичность с максимальной магнитудой более 5 при извержении вулкана Бардарбунга в Исландии в 2014 г. была вызвана инжекцией магмы, при этом объем инжекции магмы оценивается по деформациям поверхности от 0.48 до 0.51 км³, объем изверженного материала на поверхности составил 1.4 км³ (случай N, [Sigmundsson et al., 2015]).

Максимальная магнитуда при извержении вулкана Толбачик в 2012 г. составила 4.9 по данным КФ ФИЦ ЕГС РАН, при объеме изверженного материала 0.55 км³ [Belousov et al., 2015] (случай O).

Извержения вулканов Ключевской и Безымянный фиксируются тепловой дистанционной спутниковой съемкой, которая позволяет оценивать объемы извержений (без учета эксплозивных продуктов) по корреляции между тепловой мощностью и приращениями объемов построек соответствующих вулканов (https://www.mirovaweb.it/?country_id=4) [Copolla et al., 2021]. Информация по триггерной сейсмичности, предшествующей извержениям вулканов, доступна из каталогов ГС РАН. Исходя из этого, табл. 1 дополнена данными, включающими объемы извержений вулкана Ключевской (2003‒2021 гг.) и вулкана Безымянный (2000‒2020 гг.). C учетом вышесказанного, рис. 1 дополнен четырьмя точками N, O, Р и R.

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

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

$M_{\max }=0.67\cdot Ig\left(V\right)-0.82$ (9)

Моделирование инжекции магмы в сдвиговых геомеханических условиях в фундаменте Мутновского вулкана [Kiryukhin et al., 2022] показывает, что зависимость M_max от объема инжекции приближается к теоретической зависимости для воды (см. рис. 1).

ОБСУЖДЕНИЕ

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

Самая низкая сжимаемость из рассматриваемых флюидов у воды 4.9∙10^-10 Па^-1, более высокая сжимаемость у суперкритического CO₂ (плотность CO₂ при 350 К и 8 МПа составляет 164 кг/м³, а при 305 К и 30 МПа ‒ 941 кг/м³ [Pruess, 2005]) – это может завышать объемы, использованные для отображения экспериментальных точек суперкритического CO₂ на рис. 1.

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

Важным следствием из анализа распределения опытных данных, показанных на рис. 1, является то, что значения максимальных магнитуд триггерной сейсмичности не превысили теоретически обоснованного предела для воды (уравнение 7) и экспериментально установили предел для суперкритического СО₂ (уравнение 8) и магмы (уравнение 9).

Заметим также, что модель McGarr (предельная магнитуда землетрясения, которое может быть инициировано закачкой флюида, пропорциональна объему закачанного флюида) нуждается в постоянном обновлении по мере накопления экспериментальных данных. Например, землетрясение 2017 г. M_w 5.5 в Pohang (Корея), которое произошло вблизи участка при проведении гидроразрыва с закачкой ~10⁴ м³ воды [Grigoli et al., 2018], в то время как согласно модели McGarr для инициирования землетрясения такой магнитуды требуется объем закачки на три порядка больше. Приводятся также свидетельства о влиянии на сейсмичность чрезвычайно слабых вариаций гидростатического давления (порядка миллибар) [Faugler et al., 2018].

Предлагаемое нами уравнение (7) оценивает магнитуду землетрясения в Pohang (Корея) величиной 4.1. Наша оценка меньше, чем наблюдаемое 5.5. Однако она гораздо лучше, чем у McGarr. По мере накопления экспериментальных данных возможно возникнут достаточные основания чтобы скорректировать коэффициент “В” в нашем уравнении. Заметим также, что принятый объем закачки 10000 м³ – это всего лишь предположение в работе [Grigoli et al., 2018].

ЗАКЛЮЧЕНИЕ

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

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

БЛАГОДАРНОСТИ

Авторы благодарят аспиранта Кея Такахаши за помощь в расчетах объемов инжекции по литературным данным, А.В. Соломатина и А.В. Сергееву за полезные замечания на стадии подготовки рукописи. Авторы признательны рецензентам за конструктивные замечания, которые позволили усовершенствовать статью.

Работа выполнена с использованием данных, полученных на уникальной научной установке “Сейсмоинфразвуковой комплекс мониторинга арктической криолитозоны и комплекс непрерывного сейсмического мониторинга Российской Федерации, сопредельных территорий и мира” (https://ckp-rf.ru/usu/507436/).

КОНФЛИКТ ИНТЕРЕСОВ

Авторы данной работы заявляют, что у них нет конфликта интересов.

ФИНАНСИРОВАНИЕ РАБОТЫ

Исследование выполнено при финансовой поддержке РФФИ и Научным и технологическим исследовательским советом Японии в рамках научного проекта РФФИ № 21–55–50003 ЯФ_а “Магматический фракинг и флюидные потоки в вулканических структурах”.

About the authors

A. V. Kiryukhin

Institute of Volcanology and Seismology FEB RAS

Author for correspondence.
Email: AVKiryukhin2@mail.ru
Russian Federation, bulvar Piipa, 9, Petropavlovsk-Kamchatsky, 683006

Y. Fujii

Hokkaido University

Email: AVKiryukhin2@mail.ru

Faculty of Engineering

Japan, N13W8, Sapporo, 060-8628

B. A. Alam

Military Institute of Science and Technology

Email: AVKiryukhin2@mail.ru

Department of Petroleum & Mining, Faculty of Civil Engineering

Bangladesh, Dhaka, 1216

E. V. Chernykh

Institute of Volcanology and Seismology FEB RAS

Email: AVKiryukhin2@mail.ru
Japan, bulvar Piipa, 9, Petropavlovsk-Kamchatsky, 683006

References

Supplementary files