CC BY
27
УДК 550.34
МОДЕЛИРОВАНИЕ ТОЧЕЧНЫХ ИСТОЧНИКОВ МИКРОСЕЙСМИЧЕСКОГО ИЗЛУЧЕНИЯ В АНИЗОТРОПНЫХ ГОРНЫХ ПОРОДАХ
Антон Витальевич Азаров
Институт горного дела им. Н. А. Чинакала СО РАН, 630091, Россия, г. Новосибирск, Красный проспект, 54, аспирант, тел. (383)335-96-42, e-mail: [email protected]
Александр Сергеевич Сердюков
Институт горного дела им. Н. А. Чинакала СО РАН, 630091, Россия, г. Новосибирск, Красный проспект, 54, кандидат физико-математических наук, младший научный сотрудник, тел. (383)335-96-42, e-mail: [email protected]
В работе рассматривается задача моделирования микросейсмического излучения от источников точечного типа в анизотропных горных породах. В качестве примера рассматривается анизотропия, связанная с трещиноватостью горных пород. Для описания распространения упругих волн в таких породах предлагается использовать систему уравнений упругости для трансверсально-изотропных сред. Также в работе описана схема получения формул компонент тензора сейсмического момента описывающего точечные источники в случае транс-версально-изотропной среды с горизонтальной осью симметрии. Приведены численные результаты моделирования упругих волн для трехмерной неоднородной модели среды с анизотропным слоем.
Ключевые слова: уравнения упругости, упругие волны, микросейсмический источник, анизотропия горных пород, численное моделирование.
SIMULATION OF MICRO-SEISMICITY OF POINT SOURCES OF RADIATION IN ANISOTROPIC ROCKS
Anton V. Azarov
Chinakal Institute of Mining, Siberian Branch, Russian Academy of Sciences, 630091, Russia, Novosibirsk, 54 Krasny prospect, Postgraduate student, tel. (383)335-96-42, e-mail: [email protected]
Alexander S. Serdyukov
Chinakal Institute of Mining, Siberian Branch, Russian Academy of Sciences, 630091, Russia, Novosibirsk, 54 Krasny prospect, Candidate of Physico-Mathematical Sciences, Junior Researcher, tel. (383)335-96-42, e-mail: [email protected]
In the article we study the problem of modelling of microseismic propagation from the point sources in anisotropic rock mass. As an example, we study the anisotropy which appears due to rock mass fracturing. To describe the elastic waves' propagation in the rock mass we suggest using the system of elastic equations for transversely isotropic mediums. Also, scheme of obtaining seismic moment tensor component describing point sources in the case of a transversely isotropic medium with the horizontal axis of symmetry is given in the article. The numerical results of the elastic waves' modelling for three-dimensional inhomogeneous medium model with an anisotropic layer are presented in the work.
Key words: elastic equation, elastic waves, microseismic source, rock anisotropy, numerical modeling.
Использование простых моделей вмещающих сред является одним из недостатков современных алгоритмов обработки данных микросейсмического мониторинга. По причине того, что горные породы обладают анизотропией, то учет данного свойства может повысить точность определения параметров сейсмических источников.
Известно, что порода в угольных пластах содержит множество трещин. При этом размеры данных неоднородностей много меньше длин волн, соответствующих регистрируемому частотному диапазону. При разработке угольных пластов типична ситуация наличия системы трещин вертикального падения, ориентированных вдоль простирания пластов. Не ограничивая общности, рассмотрим систему трещин, параллельных плоскости (у; z) в декартовой системе координат х, у, z (где z - глубина). Процесс распространение сейсмических волн в такой ситуации хорошо описывается в рамках модели эффективной трансверсально-изотропной модели (Т1) [1], параметры которой зависят от объемной плотности трещин и инвариантны относительно поворота вокруг оси х (ось симметрии Т1). Закон Гука в матричной форме в такой среде имеет вид:
ч/ (с с Мз с Мз 0 0 0 ^
°"22 с с С -2 С ^33 44 0 0 0 ^22
°"зз с С -2 С ^33 44 с 0 0 0 £"33
°"23 0 0 0 с 44 0 0 2 £
°13 0 0 0 0 С66 0 2е1Ъ
Ч°"12 V0 0 0 0 0 с ^66 у V 12 /
где сг , £г] - компоненты тензоров напряжений и деформаций, а ненулевые компоненты матрицы жесткости могут быть выражены через параметры Томсона, широко используемые для описания упругих свойств трансверсально-изотропных сред:
С33=рКр2(2£ + 1); С44 = рК2(2у + 1); С66 = рК2; Сп=рУр2-Си=р(^1(Гр2- У52)((2^ + \)Ур2 - К2)- V2);
Для построения зависимости элементов матрицы жесткости С1} трансвер-
сально-изотропной среды (1) от объемной плотности трещин можно использовать дифференциального приближения А.С.Вавакина и Р.Л.Салганика [2]. Учет трещиноватой структуры породы производится внесением поправок в тензор податливости £ = С-1. Предполагается, что влияние трещин проявляется только в компонентах тензора податливости, связывающих нормальные напряжения и деформации, направленные перпендикулярно трещинам, а также касательные напряжения и деформации в их плоскости. Для рассматриваемой задачи это предположение формально записывается следующим образом:
^¿п-есрд зг(з)
где , Б - соответственно, компоненты тензора податливости эффективной
модели и вмещающей среды, Q и Р- поправочные коэффициенты, задаваемые системой дифференциальных уравнений:
3 3 л/го + г^р+^ро-г2) (4)
В = (Р--у)4(<2-У2)( 1 -V2), 0(0) = 1, Р(0) = 1,
где V - коэффициент Пуассона вмещающей сред, м? - переменная, изменяющаяся в пределах от нуля до заданного значения объемной плотности трещин р/г.
Для создания эффективной модели трещиноватой среды необходимо построить решение системы (4) при и> = р/г. В промежуточных точках ч'е[0,/?/( ]
оно соответствует моделям с плотностью трещин равной w. Это означает, что поправочные коэффициенты, компоненты матрицы жесткости Су. и параметры
анизотропии могут быть вычислены, как функции плотности трещин путем решения системы (4).
Микросейсмические источники любой природы проще всего описывать как разрывы поля смещений в сплошной упругой среде [3]. Рассмотрим разрыв поля смещений [и] вдоль поверхности £ с нормалью V, расположенной в области Б, заполненной упругой средой с тензором жесткости с^ .Компоненты поля смещений в точке x области Б дается формулой:
ми(х,0 {{ //, у}.ст (5)
где ^ компоненты тензора Грина.
Полю смещений, определяемому (5), соответствует набор моментов сил:
= Щ у]С1Ш<1^ (6)
Заметим, что формулы (5)-(6) не предполагают изотропности среды. В частности они остаются справедливыми и в случае трансверсально-изотропной среды, в которой закон Гука имеет вид (1). Компоненты тензора жесткости с
преобразуются в элементы матрицы жесткости Си путем замены пар индексов
(/,/)->•£, по формуле:
к = 1дуа I = р8рч + (9——— 8рч), (7)
где 51} - символ Кронекера.
В случае трансверсально-изотропной среды с горизонтальной осью симметрии осью симметрии, используя формулы (5), (6), (7) получаем следующие
выражения для ненулевых компонент тензора жесткости с использованием параметров Томпсена [1]:
С1Ш = РУР2> С2222 = С3333 = РУр\2б + 1),
С2233 C3322 = pVp2(2s + \)-2pVs2(2y + \),
С1122 ~~ C1133 = сзз11=C2211=Цфр2 -Vs 2)((2¿+1)vp2 -Vs2)-Vs2),
С2323 — C3223 = C3232 = C2332 =pVS2(2y + l),
С2112 — C2121 ~ C1221 — C1212 — C1331 — C1313 — C3113 — C3131 — pVs2.
Для получения моментов эквивалентных пар сил для микросейсмических источников в трансверсально-изотропной среде необходимо подставить выражения в выражения (5), (6).
В качестве примера, рассмотрим модель точечного источника типа трещина отрыва. Этот источник может подходить для описания сейсмический эмиссии, связанной с фильтрационными процессами. Пусть подвижка при трещине отрыва проходит в направлении нормали к поверхности отрыва. Углы ф, 5, Л определяющие ориентацию плоскости разрыва и направление вектора смещения показаны на рис. 1.
Компоненты вектора нормали v к поверхности разрыва и вектора смещения [u], определяются для рассматриваемого типа источника по формулам:
\их ] = -и sin 8 sin ф, [и2 ] = и sin 8 cos ф, [w3 ] = -и cos 8,
(7)
Vj = - sin 8 sin ф ^ v2 = sin¿> COS<^, V3=-COS£>, v '
где, индекс 1 соответствует оси x, 2 - y, 3 - z, a u есть модуль вектора смещения. Подставляя (7) в (6), для трансверсально-изотропной среды получаем:
Мп = иАСп sin2 8 sin2 ф + iiAC,l3 sin2 8 eos2 фс22U + nAC,l3 eos2 8 Mn = -uAC66 sin2 8&\п2ф =Mn Ml3 = uAC66 sin 28 sin ф = M3l
o o o o o
M22 = uACX3 sin 8 sin ф + uAC33 sin 8 eos ф + nA(C33 - 2C44) eos ¿> M23 = -м^4С44 sin 2¿> eos ф = M32
• o • o ✓ o o
M33 = м^4С13 sin ¿> sin ф + uA(C33 - 2C44) sin 8 eos ф + uAC33 eos 8
где А есть площадь поверхности разрыва ¿/I, а и модуль вектора разрыва смещений, а элементы. По аналогии выводятся компонент тензора сейсмического для чистого сдвига, но из-за громоздкости здесь не приводятся.
Далее в работе проводилось моделирование волновых полей от точечных источников расположенных в анизотропных горных породах. Для численных расчетов использовалась трехмерная скоростная модель среды изображенная на рис. 2а). Анизотропия вводилась только в слое, выделенного на рис. 2а) черным цветом.
у, м
X, м
а) б)
Рис. 2. а) скоростная модель среды; б) снимок волнового поля в момент времени 1=0.045 с
Решения уравнений упругости проводилось конечно-разностным методом, используя схему Вирьё на сдвинутых сетках[4]. За основу при моделировании был взят программный пакет 8орЫ3ё[5]. В качестве примера, на рис. 2б) приводится снимок волнового поля, от источника сдвигового типа, который был расположен в анизотропном слое в момент времени 0.045 секунд.
Выводы.
Предложена схема моделирования сейсмических полей от точечных источников в анизотропных горных породах. Выведены формулы компонент тензора сейсмического момента, в случае трансверсально-изотропной среды. Проведено численное моделирование сейсмических полей от точечных источников в трехмерной неоднородной среде с анизотропным слоем. Полученные результаты могут быть использованы для построения синтетических сейсмограмм.
Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 16-35-00513 мол а.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Thomsen L. Weak elastic anisotropy // Geophysics. - Т.51. - № 10. - С. 1954-1966.
2. Салганик А.С., Вавакин Р.Л. Эффективные упругие характеристики тел с изолированными трещинами, полостями и жесткими неоднородностями // Известия АН СССР, Механика твердого тела. - 1978. - С. 95-107.
3. Аки К., Ричардсон П. Количественная сейсмология. Теория и методы: в 2-х томах. -М., 1983.
4. Virieux J. P-SV wave propagation in heterogeneous media: Velocity-stress finite-difference method //Geophysics. - 1986. - Т. 51. - №. 4. - С. 889-901.
5. Bohlen T. Parallel 3-D Viscoelastic finite difference seismic modelling //Computers & Geosciences. - 2002. - Т. 28. - №. 8. - С. 887-899.
© А. В. Азаров, А. С. Сердюков, 2016
