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

GEODYNAMICS & TECTONOPHYSICS

PUBLISHED BY THE INSTITUTE OF THE EARTH'S CRUST SIBERIAN BRANCH OF RUSSIAN ACADEMY OF SCIENCES

ISSN 2078-502X

2019 VOLUME 10 ISSUE 2 PAGES 471-481

https://doi.org/10.5800/GT-2019-10-2-0422

Prediction of compression and extension zones in geological

structures based only on the velocities of longitudinal waves in the geological medium

B. P. Sibiryakov1, 2, E. A. Khogoev1

1 A.A. Trofimuk Institute of Petroleum Geology and Geophysics, Siberian Branch of RAS, Novosibirsk, Russia

2 Novosibirsk State University, Novosibirsk, Russia

Abstract: The article presents accurate solutions for the problem for two elastic half-spaces with an arbitrary curvilinear interface. Our study shows that dilatation solutions (Poisson integrals) are dependent on neither an overall compression modulus nor the Poisson ratio, and depend only on the velocity of longitudinal waves. These specific solutions can be supplemented by general solutions for an incompressible elastic medium, and the boundary conditions of the rigid contact for the sum of the solutions can thus be satisfied. Relatively simple calculations make it possible to determine the divergence of the displacement field and reduce the entire problem solving process to a study of Poisson equations with a known divergence. Furthermore, predictions of volumetric compression or extension are important for geological investigations, since the zones characterized by reduced pressure rates may act as fluid at-tractors.

Key words: dilatation; forecasting; prediction; pressure; velocity of longitudinal waves; fundamental solutions; Poisson integral; random walk method (RWM) for elliptic equations

RESEARCH ARTICLE Received: September 13, 2018

Revised: March 15, 2019 Accepted: April 3, 2019

For citation: Sibiryakov B.P., Khogoev E.A., 2019. Prediction of compression and extension zones in geological structures based only on the velocities of longitudinal waves in the geological medium. Geodynamics & Tectonophysics 10 (2), 471-481. doi:10.5800/GT-2019-10-2-0422.

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

Б. П. Сибиряков1- 2, Е. А. Хогоев1

1 Институт нефтегазовой геологии и геофизики им. А.А. Трофимука СО РАН, Новосибирск, Россия

2 Новосибирский национальный исследовательский государственный университет, Новосибирск, Россия

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

Ключевые слова: дилатация; прогноз; давление; скорость продольных волн; фундаментальные решения;

интеграл Пуассона; метод случайных блужданий для эллиптических уравнений

1. Введение

Методы структурной сейсмологии используются в основном для определения конфигурации сейсмических границ, а также для определения параметров геологической среды, прежде всего скоростей продольных и поперечных волн [Suvorov et al., 2002; Puzirev et al., 1979; Puzyrev, 1993]. Однако возможности этих методов гораздо шире. В принципе, возможен прогноз состояния среды, прежде всего напряженного состояния. Известно, что расчет напряженного состояния по данным о пластовых скоростях продольных и поперечных волн, а также плотности пород и конфигурации структур сводится к интегрированию уравнений равновесия упругого тела с граничными условиями непрерывности векторов перемещений и нагрузок на сейсмических границах [Sibiryakov et al., 2004, 2013, 2015]. Существуют методы, использующие все точки трехмерного пространства, такие как метод конечных элементов, и методы более экономичные, где результирующая линейная система содержит информацию только о точках границ геологической среды. Это метод граничных интегральных уравнений [Kupradze, 1963]. Тем не менее вопрос о прогнозе напряженного состояния нельзя считать полностью решенным. Дело здесь не только в том, что большие линейные системы требуют методов регуляризации, что снижает достоверность расчетов.

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

2. Постановка задачи

Рассмотрим простейшую ситуацию. Пусть верхней границей полупространства является плоскость г = 0. Поверхность г = г2(х,у) есть поверхность структуры, разделяющей упругие среды, а

x

Рис. 1. Схема границ в среде: (а) - в трехмерном виде; (б) - в сечении по координате Y=0. Fig. 1. Schematic boundaries in the medium: (а) - 3D; (б) - cross section along coordinate Y=0.

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

Задача сводится к решению уравнения упругого равновесия для двух сред, разделенных произвольно сложной поверхностью структуры г = г2(х,у). На этой поверхности задано шесть граничных условий, т.е. равенство трех компонент вектора нагрузок и вектора перемещений. В то же время на границе г = 0 задаются лишь три граничных условия, т.е. нулевые значения вектора нагрузок.

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

3. Метод решения

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

3.1. Частные решения уравнения равновесия упругого тела

Уравнения упругого равновесия для поля перемещений и

(Л + ¿и)дгайй1уи + /лАи = pg

имеют частные решения в виде свертки тензора Грина с вектором силы тяжести [Kupradze, 1963] с учетом того, что вектор g имеет только одну - вертикальную - компоненту:

«2М = ^ШГ*

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

rfcz(x,y) = [(1+ y%fc + d-r2)(Xfc-ryg:)-y3)]

1

г2(х,у) ] г(х,у)' V2

где у2 = а 5zfe - символ Кронекера. Интегрирование предполагается по объему структуры. В точке x по-

ур

ле вычисляется, а точка у есть переменная точка объема в процессе интегрирования, поэтому символ dVy означает переменный элемент объема структуры. Для задач геодинамики интеграл Пуассона (1) представляет особый интерес. Прежде всего, ускорение силы тяжести практически постоянно и оно не участвует в процессе интегрирования. Кроме того, поле перемещений (1) зависит от скорости продольных и поперечных волн и от объема структуры, но не зависит от плотности. Дивергенция поля перемещений 0о(х) (дилатация поля перемещений) не зависит ни от плотности, ни от скорости поперечных волн, а только от скорости продольных волн в соответствии с выражением

90(х) = divu0Jx) = ow fcw 4nV2 JJJ dz\r(x,y)J У

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

6om = ibm

4 nVf

r

dxdy.

В интеграле (2) M - внутренняя точка геологической структуры, которая, в частности, может и принадлежать границе структуры, а поверхности г2(х,у), г1(х,у) ограничивают объем структуры сверху и снизу. Интегрирование ведется по плоскости проекции этих поверхностей. Расстояние между точками M и поверхностью г(х,у) обозначено г[М, г(х,у)].

Мы можем рассматривать аномальное значение дилатации, выражаемое разностью между общей объемной деформацией и объемной деформацией любого элемента в условиях отсутствия структуры:

4п \V2 V£/JJ

r

dxdy.

Частное решение уравнений теории упругости (2)-(3) вычисляется достаточно просто, однако оно не обязано удовлетворять граничным условиям непрерывности векторов перемещений и нагрузок на поверхности, разделяющей разные среды.

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

О вычислении дилатации и давления. Если имеется слой на полупространстве без каких-либо структур, то существует элементарное решение для напряжений:

^ = (Л+2/г)е22 = pgz; ахх = ауу = (1-2у2)рдг.

Это связано с отсутствием горизонтальных перемещений, а также всех деформаций, кроме одной, вертикальной деформации е22, так что закон Гука дает вертикальные напряжения в виде а22 = {Л+2/и)егг, горизонтальные напряжения в форме ахх = оуу = а давления даются выражением [Sibiryakov et al., 2004]:

Р =

=

В этом элементарном случае давление существенно меньше веса вышележащих пород рдг.

При наличии структуры в поле дилатаций вносятся интегралы по объему структуры в0(М+) — в0(М~), связанные с различными скоростями продольных волн в слое и полупространстве. Иными словами, убирается влияние материала верхнего слоя и добавляется влияние материала нижнего полупространства на дилатацию. Естественно, что само поле дилатации разрывается на всей поверхности раздела слоя и полупространства, осложненных структурой. Возникает необходимость пополнения частного решения общим для выполнения граничных условий непрерывности векторов перемещений и нагрузок. Так как сделана попытка устранить влияние силы тяжести с помощью частного решения, можно попытаться искать общее решение, соответствующее упругой несжимаемой среде.

3.2. Закон Гука для несжимаемой среды

В несжимаемой среде поле перемещений удовлетворяет уравнению Лапласа

Диг = 0.

Напряжения и деформации в несжимаемой среде связаны соотношением:

,, {дщ дик\

а для несжимаемой потенциальной среды (и = дгайФ~) закон Гука есть

= 2 д2Ф

Вектор нагрузок на границе г = г0(х,у) принимает вид:

П о д2Ф д2Ф дщ

Р1 = а1кпк = 2И^тПк = 2ц—=2м~. С граничными условиями для полных перемещений:

и* ~и1 =Щ0

Поля перемещений и нагрузок соответствуют частным решениям и вычисляются достаточ-

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

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

Представим ситуацию, когда структура разделяет два слоя с криволинейной поверхностью и различными упругими свойствами и плотностью, так что верхняя граница первого слоя есть некоторая плоскость г = к. Тогда нижняя граница второго слоя, соответственно г = Н. Поле перемещений на границе г = г0(х,у) можно определить в приграничной точке М+, отстоящей от точки М границы на расстоянии Д по нормали. Поле И((М+) = А*кик + В*кУк, где матрицы А^к есть вероятности перехода от границы г = г0(х,у) в точку М+ а В*к из любой точки верхней границы, например свободной поверхности г = к в ту же самую точку. Аналогично поле по другую сторону поверхности г = г0(х,у) в точке М~ дается равенством:

щ(М~) = А~1кик +в;кук.

Матрица В'[к есть вероятность попадания из приграничной точки М~ в одну из точек нижней плоскости г = Н.

Нормальная производная на границе г = г0(х,у) может быть определена как отношение разности поля перемещений в граничной и приграничной точках, а именно:

dui _ (Sik-AÎk)uk , ^г>+ D + + A,,+ - -Д--\Bik~ Bik )Vk ,

где Д - шаг дискретизации, В*к+ - вероятность попадания не в приграничную точку М+, а в точку границы М. Эта разность может быть как угодно малой. Граничное условие равенства вектора нагрузок дает соотношение:

Pitofc + 0 =M2

Система уравнений (5) не замкнута, так как содержит два множества неизвестных и^ и ик. Поля Уд" и vk известны, так как предполагаются известными нагрузки и перемещения на плоскостях, ограничивающих рассматриваемую область. Однако существует и условие равенства поля перемещений, которое замыкает линейную систему (5).

Второе граничное условие для системы уравнений (5) гораздо более простое, а именно:

ufc +ико =ик +ико.

Тем самым определяется полное поле перемещений как по одну, так и по другую сторону границы z = z0(x,y), т.е. сумма частного и общего решений

\pi1(sik -А&) -ц2(sik -¿rfc)K"fc +0 = Д(piro -О + -B0T)vû - (Bifc -SiTK.

Таким образом, имеется неоднородная система для определения поля перемещений на границе раздела сред, если заданы поля перемещений на удаленных от структуры границах z = 0 и z = Н. Если пренебречь влиянием удаленных границ, то поле перемещений и остальные элементы напряженного состояния вычисляются явно. Таким образом, граничные условия удовлетворяются, и эвристический прием использования частных решений в виде интегралов Пуассона и несжимаемой среды может считаться оправданным.

4. Численное решение задачи расчета дилатации и давления для геологических структур ПО ДАННЫМ только о скоростях продольных волн

Для оценки реальных значений аномальной дилатации проведем расчет по формуле (3) при типичных значениях перепада скоростей. Граница раздела пластов представляет собой антиклинальное поднятие (рис. 2), описываемое формулой:

где гт = 300 м, а = $ = 400.

Это купол амплитудой 300 м на глубине 3 км, с достаточно крутыми углами подъема, так что поведение деформаций сжатия и растяжения является не элементарным. Точка расчета функции определяется координатами (х0,у0,г0) Интегрирование проводится по квадратной области - 2000, 2000 по каждой из координат х, у:

М (д^Уо^о) = ^ - Й™02000 £у0=0-02000 " Т^Г^^ТЩ^) ^^

Результат расчетов представлен на рис. 3, а, б, при двух вариантах перепада скоростей соответственно, 13=3000, У2=4000 м/с и 13=4000, У2=3000 м/с. Для удобства представления результат масштабируется умножением на 106. Области сжатия обозначены синим цветом, а области растяжения - красным.

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

-2800

-2000

Рис. 2. Поверхность раздела сред с различной скоростью продольных волн, модель - антиклиналь. Fig. 2. Interface between the media that differ in the velocities of longitudinal waves. Anticline model.

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

Другая модель разделяющей границы - квадратный выступ (целик) размерами 1x1 км, с тем же перепадом глубин, что и в предыдущем случае (рис. 4).

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

-2000

Рис. 3. Поле аномальной дилатации Д0 при Y=0: (а) - Vi=3 км/с, V2=4 км/с; (б) - Vi=4 км/с, V2=3 км/с; вертикальная ось - глубина Z, горизонтальная ось - координата X; значения Д0 кодируются цветом, шкала справа.

Fig. 3. Field of anomalous dilation Д0 at Y=0: (а) - V1=3 km/sec, V2=4 km/sec; (б) - V1=4 km/sec, V2=3 km/sec; vertical axis - depth Z; horizontal axis - coordinate X; values Д0 are colour-coded; the scale is shown on the right.

-2700-1

n \

-2900^

-3000 1000 V

Y

-10007

-2000 -i-1-1-1-*

-2000 -1000 0 1000 -2000

X

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

Fig. 4. Interface between the media that differ in the velocities of longitudinal waves. Rectangular rock edge model.

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

На рис. 6 изображено одномерное поведение дилатации для купола (слева) и целика (справа) с глубиной для плоскости, проходящей через центр структуры при различных соотношениях скоростей продольных волн сверху и снизу. Рассмотрено четыре варианта соотношения скоростей в первом (верхнем) и втором слое: 1) 71=3000, 72=4000 м/с; 2) 71=3500, 72=4000; 3) 71=4000, 72=3500 м/с и 4) 71=4000, 72=3000 м/с. Обращает на себя внимание совпадение точек пересечения всех кривых и смена их знака в окрестности

Рис. 5. Поле аномальной дилатации ДО при Y=0; вертикальная ось - глубина Z, горизонтальная ось - координата X, значения ДО кодируются цветом, шкала справа. (а) - Vi=3 км/с, (б) - V2=4 км/с.

Fig. 5. Field of anomalous dilation ДО at Y=0: (а) - Vi=3 km/sec; (б) - Vi=4 km/sec; vertical axis - depth Z; horizontal axis -coordinate X; values Д0 are colour-coded; the scale is shown on the right.

Рис. 6. Аномальная дилатация Д0по вертикали при X=Y=0 при различных вариантах перепада скорости продольных волн V1,V2 на границе, указанных в центре рисунка: (а) - модель - антиклиналь, (б) - модель - прямоугольный выступ.

Fig. 6. Anomalous dilation Ав at X=Y=0 in case of different variants of the differential velocity of longitudinal waves V1 and V2 at the boundary, as shown in the center of the figure: (a) - anticline model, (б) - rectangular rock edge model.

горизонтальной плоскости. Уменьшение перепада скорости на 500 м/с приводит к уменьшению перепада аномальной дилатации приблизительно в два раза.

Давление P в среде связано с дилатацией формулой Р = ^Л + ^ß^O, или иначе Р = (Л + 2/1) (l — ^У2) 0 =

7 ( 4 о\ о Л+2/

= р Vp -y^J 0, так как V^ = Для оценки измененного давления в среде, связанного с рассчитанной

дилатацией, установим следующие (типичные) значения параметров: у2=1/3; р=3000 кг/м3, Vp=3000 м/с.

Изменение давления, обусловленного аномальной дисперсией, которое в наших примерах составляет от A0min=-5O.O-1O-6 до Д#тах=50.0-10-6, или |Д#|=10-4, может быть оценено как

АР = 3 • 103 • (3-103)2 (l- |Д0| = 27 • 109Ц |Д0| = 15 • 109 -10"4 = 15 • 105 Па.

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

Р = 3-103 -9.81-3-103 (l- = 9.81 •lO6 -5 = 49-106 = 49 • 106 Па.

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

5. Обсуждение

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

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

6. Заключение

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

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

7. Литература / References

Kupradze V.D., 1963. Methods of Potential in the Theory of Elasticity. Fizmatgiz, Moscow, 472 p. (in Russian) [Купра-дзе В.Д. Методы потенциала в теории упругости. М.: Физматгиз, 1963. 472 с.].

Puzirev N.N., Mandelbaum M.M., Krylov S.V., Mishenkin B.P., Mishenkina Z.R., Petrick G.V., Seleznev V., 1979. New data from explosion seismology in the Baikalian rift zone. Tectonophysics 56 (1-2), 128. https://doi.org/10.1016/ 0040-1951(79)90023-4.

Puzyrev N.N. (Ed.), 1993. Detailed Seismic Studies of the Lithosphere Using P- and S-waves. Nauka, Novosibirsk, 199 p. (in Russian) [Детальные сейсмические исследования литосферы на P- и S-волнах / Ред. Н.Н. Пузырев. Новосибирск: Наука, 1993. 199 с.].

Sibiryakov B.P., Leite L.W.B., Sibiryakov E.P., Vieira W.W.S., 2015. Local low pressure areas in anticline structures. Revista Brasileira de Geofísica 33 (2), 225-236. https://doi.org/10.22564/rbgf.v33i2.716.

Sibiryakov B.P., Leite L.W.B., Vieira W.W.S., 2013. Model of the structured continuum, and the relation between specific surface area, porosity and permeability. Revista Brasileira de Geofísica 31 (4), С. 559-568. https://doi.org/ 10.22564/rbgf.v31i4.337.

Sibiryakov B.P., Sibiryakov E.B., Glebov A.F., Nesterov V.N., Sokolov E.P., 2004. Prediction of stress and hydrodynamic behavior of fluids from multicomponent seismic data. Geologiya i Geofizika (Russian Geology and Geophysics) 45 (6), 709-715.

Suvorov V.D., Mishenkina Z.M., Petrick G.V., Sheludko I.F., Seleznev V.S., Solovyov V.M., 2002. Structure of the crust in the Baikal rift zone and adjacent areas from deep seismic sounding data. Tectonophysics 351 (1-2), 61-74. https:// doi.org/10.1016/S0040-1951(02)00125-7.

СВЕДЕНИЯ ОБ АВТОРАХ | INFORMATION ABOUT AUTHORS

Борис Петрович Сибиряков

докт. физ.-мат. наук, профессор, в.н.с.

Институт нефтегазовой геологии и геофизики им. А.А. Трофимука СО РАН

630090, Новосибирск, пр. Академика Коптюга, 3, Россия

Новосибирский национальный исследовательский

государственный университет

630090, Новосибирск, ул. Пирогова, 2, Россия

И e-mail: SibiryakovBP@ipgg.sbras.ru © https://orcid.org/0000-0002-4634-7769

Boris P. Sibiryakov

Doctor of Physics and Mathematics, Professor, Lead Researcher

A.A. Trofimuk Institute of Petroleum Geology and Geophysics, Siberian Branch of RAS

3 Academician Koptug ave., Novosibirsk 630090, Russia

Novosibirsk State University 2 Pirogov street, Novosibirsk 630090, Russia

Евгений Андреевич Хогоев

канд. физ.-мат. наук, с.н.с.

Институт нефтегазовой геологии и геофизики им. А.А. Трофимука СО РАН

630090, Новосибирск, пр. Академика Коптюга, 3, Россия

И e-mail: HogoevEA@ipgg.sbras.ru © https://orcid.org/0000-0002-8941-6840

Evgenii A. Khogoev

Candidate of Physics and Mathematics, Senior Researcher

A.A. Trofimuk Institute of Petroleum Geology and Geophysics, Siberian Branch of RAS

3 Academician Koptug ave., Novosibirsk 630090, Russia