Комплексирование акустических, гравитационных и геомеханических полей в иерархических средах

Хачай Ольга Александровна; Хачай Олег Юрьевич; Хачай Андрей Юрьевич

О.А. Хачай, О.Ю. Хачай, А.Ю. Хачай

комплектирован™

АКУСТИЧЕСКИХ, ГРАВИТАЦИОННЫХ И ГЕОМЕХАНИЧЕСКИХ ПОЛЕЙ В ИЕРАРХИЧЕСКИХ СРЕДАХ*

Построены итерационные алгоритмы 2-D моделирования для дифракции звука и линейно поляризованной поперечной упругой волны на включении с иерархической упругой структурой, расположенной в J-ом слое N-слойной упругой среды. Рассмотрен случай, когда плотность включения каждого ранга отличается от плотности вмещающей среды, а упругие параметры совпадают с упругими параметрами вмещающего слоя, а также, когда плотность включения каждого ранга совпадает с плотностью вмещающей среды, а упругие параметры включения каждого ранга отличаются от упругих параметров вмещающего слоя. Использован метод интегральных и интегро-дифференциальных уравнений для пространственно-ча-статного представления распределения волновых полей. Из построенной теории следует, что при комплексировании акустических, геомеханических и гравитационных полей необходимо использовать такие данные, которые получены в рамках систем наблюдения, настроенных на исследование иерархической структуры среды. Ключевые слова: иерархическая среда, гравитационное поле, акустическое поле, геомеханическое поле, итерационный алгоритм 2D моделирования.

Введение

Последние десятилетия характеризуются активным развитием наук о Земле. Мы воспользуемся материалами, изложенными в [1], написанными академиком А.Н. Дмитриевским. «Геофизика XX века — это понимание того, что:

• геофизические поля являются индикаторами процессов, протекающих в литосфере; геофизические параметры, регистрируемые дистанционно, могут иметь функциональную или тесно-корреляционную связь с вещественно-структурными характеристиками геологической среды (на макро-, мезо- и микро уровнях);

* Работа выполнена при поддержке Программы государственной поддержки ведущих научных школ (НШ-9356.2016.1).

УДК 550.830, 539.3

• анализ пространственно-временного и энергетического распределения геофизического поля может дать информацию о пространственно-временном распределении свойств геологической среды;

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

Геофизика позволила ответить на следующие принципиальные для геолога вопросы:

• какова глубина залегания и геометрия (иногда морфология) исследуемого геообъекта;

• каков вещественно-структурный состав геообъекта;

• где и как расположены субвертикальные и субгоризонтальные ноднородности и прежде всего зоны макротрещиноватости;

• каковы условия фильтрации для флюидов-порозаполни-телей;

• каковы термодинамические условия залегания искомого объекта.

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

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

В работах [2—4] построены алгоритмы моделирования в электромагнитном случае для 3^ неоднородности, в сейсмическом случае для 2^ неоднородности для произвольного типа источника возбуждения ^-слойной среды с иерархическим проводящим или упругим включением, расположенным в /-ом слое. Построен итерационный алгоритм решения обратной задачи для 2^ дифракции линейно поляризованной упругой поперечной волны для ^-слойной упругой среды с упругим аномально напряженным включением иерархического типа на основе уравнения теоретической обратной задачи [5]. В настоящей работе будет исследован вопрос о комплексировании разных типов геофизических полей, распространяющихся в иерархических структурных средах.

Моделирование дифракции звука

на двумерной аномально плотностной неоднородности

иерархического типа, расположенной в Л-слойной

упругой среде

В работе [3] описан алгоритм моделирования дифракции звука на двумерном упругом иерархическом включении, расположенном в J-ом слое ^-слойной среды. С5р:]- (М, М0) — функция источника сейсмического поля, краевая задача для которой сформулирована в работе [2], Ь/^ч = га2(ст^ / X) — волновое число для продольной волны, в приведенном выражении индекс ji обозначает принадлежность свойств среды внутри неоднородности, ja — вне неоднородности, X — постоянная Ламэ, а — плотность среды, га — круговая частота, и = дтайф — вектор смещений, ф0 — потенциал нормального сейсмического поля в слоистой среде в отсутствие неоднородности: ф0 = ф°а. Будем считать, что упругие параметры иерархического включения для всех рангов I и вмещающего слоя одинаковы, а плотность иерархического включения для всех рангов отличается от плотности вмещающей среды, тогда система уравнений [3] перепишется в виде:

) ДОф (М)О^ (М, м 0) -

2% Бс1

-$ о ^¿с = ф(М°), М0 е 3С1 щ2% а дп

¡и (Ки - К')

Ц ф, (м)о5р ] (М, М0 )йтм + ф0-! (м0) -

т(М °)2 П м

г0\ лжО

¡а У« К- 5ф;

= ф(М0),М0 е 5с

ст(М0)2п 5р,У дп

0Бру (м, м0) — функция источника акустического поля, она совпадает с функцией [3], к2]-и = ю2(ст¡и / (X¡и ); X¡и = X¡а; - волновое число для продольной волны, в приведенном выражении индекс ji обозначает принадлежность свойств среды внутри неоднородности, ja — вне неоднородности, I = 1..Х — номер иерархического уровня, ф0 — потенциал нормального акустического поля в слоистой среде в отсутствие неоднородности предыдущего ранга, если I = 2..Х ф0 = ф;-1, если I = 1, ф0 = ф0, что совпадает с соответствующим выражением [3].

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

Моделирование дифракции упругой поперечной волны на аномально плотностной неоднородности иерархического типа, расположенной в ^-алойной упругой среде

Аналогично (1) выписывается такой же процесс для моделирования распространения упругой поперечной волны в ^-слой-

ной среде с двумерной иерархической структурой произвольной морфологии сечения с использованием интегральных соотношений, выписанных в работе [3].

(/ш - К )

2п

И их1 (М)О^(М, М

(2)

+и°х(1 _0(м°) = их1 (Мх),Мх е 3С1

Д и-(м)^(м,М +

_)(Мх) = их1 (Мх),Мх й

где 05а(М, М0) — функция источника сейсмического поля рассматриваемой задачи, она совпадает с функцией Грина, выписанной в работе [3] для соответствующей задачи, = га2 (а/ ) — волновое число для поперечной волны, ^ — постоянная Ламэ,

и , — составляющая вектора смещений, I = 1...Х — номер иерарх «

хического уровня, их1 — составляющая вектора смещений сейсмического поля в слоистой среде в отсутствие неоднородности предыдущего ранга, если I = 2..Х и0х1 = их(1 -1), если I = 1, и0х1 = и- , что совпадает с соответствующим выражением для нормального поля в работе [3]. Следует отметить что структура уравнений (1) совпадает с общим случаем, когда иерархическая неоднородность имеет не только плотностные параметры отличные от параметров вмещающей среды, но и упругие параметры на всех рангах отличаются от упругих параметров вмещающего слоя. Отличие этой задачи заключается только в значениях волнового числа. Таким образом, более чувствительным к области плотностных неоднородностей в массиве является отклик среды, связанный с поперечной волной. Это следует учитывать при картировании сложно организованной геологической среды.

Моделирование дифракции звука на двумерной аномально напряженной неоднородности иерархического типа, расположенной в Л-слойной упругой среде

В работе [3] описан алгоритм моделирования дифракции звука на двумерном упругом иерархическом включении, расположенном в /-ом слое ^-слойной среды. С5р}- (М, М0) — функция источника сейсмического поля, краевая задача для которой сформулирована в работе [2], = га2 (а ~ / X) — волновое число для продольной волны, в приведенном выраже-

нии индекс^ обозначает принадлежность свойств среды внутри неоднородности, ja — вне неоднородности, X — постоянная Ламэ, а — плотность среды, га — круговая частота, и = дтайф —

вектор смещений, ф0 — потенциал нормального сейсмического

0 0

поля в слоистой среде в отсутствие неоднородности: ф/; = ф;-а . Будем считать, что плотность иерархического включения для всех рангов I и вмещающего слоя одинаковы, а упругие параметры иерархического включения для всех рангов отличаются от упругих параметров вмещающей среды, тогда система уравнений [3] перепишется в виде:

(Ъ2 - Ъ2 )

- §§ф (М)О^ (М, М°)йтм +

2п ^^

+ф0-1(М0) = ф(М0),М0 е БС1 (к1'" - ^) Ц Ф, (М)С5р ] (М, М0 )йтм

(3)

а(М°)2п м

+ф0-1(м0) = ф,(м0),м0 Й Scl Обозначения те же, что и для системы уравнений (1).

Моделирование дифракции упругой поперечной волны на аномально напряженной неоднородности иерархического типа, расположенной в Л-слойной упругой среде

Аналогично (3) выписывается такой же процесс для моделирования распространения упругой поперечной волны в ^-слойной среде с двумерной иерархической структурой произвольной морфологии сечения с использованием интегральных соотношений, выписанных в работе [3].

^ -^ Ц их1 (М)О^(М,М°Мтм +^Ч0(1-1)(М°) +

2п БСЦ ^¡и

§"" (м) дП

Рщ2п С, дп (4)

* (М0Ц и,(М)О^.(М, М°)йтм + и^ш0) +

+ § и, (М) ^ йс = и, (М0), М0 е

ц(М°)2п 5Ск

\и, (М)^^-йс = и, (М0), М0 е 5С

ja М^ )Г .4 , _ /,.0 4 д.0

ц(М°)2п С, х1 дп

где (М, М0) — функция источника сейсмического поля рассматриваемой задачи, она совпадает с функцией Грина, выписанной в работе [3] для соответствующей задачи, = га2 (а/ ) ф Цja — волновое число для поперечной волны, а= а]а, ^ — постоянная Ламэ, их1 — составляющая вектора смещений, I = 1..^ — номер иерархического уровня, и^ — составляющая вектора смещений сейсмического поля в слоистой среде в отсутствие неоднородности предыдущего ранга, если I = 2..Х и— = их(1 _0 , если I = 1, их1 = их , что совпадает с соответствующим выражением для нормального поля в работе [3]. Следует отметить, что структура уравнений (4) совпадает с общим случаем, когда иерархическая неоднородность имеет не только упругие параметры отличные от параметров вмещающей среды, но и плот-ностные параметры на всех рангах отличаются от плотностных параметров вмещающего слоя. Отличие этой задачи заключается только в значениях волнового числа. Таким образом, более чувствительным к области упругих неоднородностей в массиве является отклик среды, связанный с продольной волной. Это следует учитывать при картировании сложно организованной геологической среды.

Заключение

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

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

СПИСОК ЛИТЕРАТУРЫ

1. Дмитриевский А. Н. Избранные труды. Т. 2. — М.: Наука, 2009. — 435 с.

2. Хачай О. А., Хачай А.Ю. О комплексировании сейсмических и электромагнитных активных методов для картирования и мониторинга состояния двумерных неоднородностей в Ж-слойной среде // Вестник ЮУРГУ. Серия «Компьютерные технологии, управление, радиоэлектроника». - 2011. - № 2(219). - С. 49-56.

3. Хачай О. А., Хачай А.Ю. Моделирование электромагнитного и сейсмического поля в иерархически неоднородных средах // Вестник ЮУРГУ. Серия «Вычислительная математика и информатика». -2013. - Т. 2. - № 2. - С. 49-56.

4. Хачай О. А., Хачай О. Ю., Хачай А. Ю. Новые методы геоинформатики мониторинга волновых полей в иерархических средах // Геоинформатика. - 2015. - № 3. - С. 45-51.

5. Хачай О.А., Хачай О.Ю., Хачай А.Ю. Определение поверхности аномально напряженного включения в иерархической слоисто-блоковой среде по данным акустического мониторинга // Горный информационно-аналитический бюллетень. - 2016. - № 4. - С. 354-366.

КОРОТКО ОБ АВТОРАХ

Хачай Ольга Александровна - доктор физико-математических наук,

ведущий научный сотрудник, старший научный сотрудник,

Институт геофизики им. Ю.П. Булашевича

Уральского отделения РАН, e-mail: olgakhachay@yandex.ru,

Хачай Олег Юрьевич1 - кандидат физико-математических наук, доцент,

Хачай Андрей Юрьевич1 - кандидат физико-математических наук, доцент,

1 Уральский федеральный университет

им. первого Президента России Б.Н. Ельцина.

Gornyy informatsionno-analiticheskiy byulleten'. 2017. No. 4, pp. 328-336. O.A. Khachay, O.Yu. Khachay, A.Yu. Khachay COMPLEXING ACOUSTIC, GRAVITATIONAL AND GEOMECHANICAL FIELDS IN HIERARCHICAL MEDIA

A very significant result of previous century is a conclusion of fundamental significance of the role of block hierarchic structure of rocks and massive for explaining the existence of

UDC 550.830, 539.3

ззБ

a wide set of nonlinear geomechanical effects and occurring complicated self organizing geo systems by analyze of forming large and super large deposits. The hierarchic structure is characteristic of many systems. Especially for the Earth's lithosphere, for which more than 30 hierarchic levels had been divided using geophysical information. Thus it is needed to take that into account by constructing new complex systems for geophysical investigation of the Earth's lithosphere. There are constructed iteration algorithms for 2-D modeling of sound wave diffraction and linear polarized transversal elastic wave on an inclusion with hierarchic elastic, located in the J-th layer of the N-layered elastic medium. It is considered a case, when the density of the inclusion of each rank differs from the density of the surrounded layer, and the elastic parameters are equal to the elastic parameters of the surrounded layer. Also we considered the case when the density parameters of the hierarchic inclusions are equal to the density of the layer, where they are located and their elastic parameters differ from the parameters of the layer. We used the method of integral and integral-differential equations for space frequency representation of the wave field distribution. From the constructed theory it is evident that by integrating seismic and gravitation fields it is needed to use such data, which are obtained by the systems of observation, that are set to observe the medium hierarchical structure. The use the values of density, derived from the correlation dependence between the velocity values of the longitudinal wave, derived from the kinematic interpretation of seismic data, and density, for constructing the density model using gravitational data can lead to the discrepancy of the obtained model and matter content of the investigated geological medium. These results are the base of constructing new mapping systems of geological systems. Especially it is significant for mapping oil and gas deposits and for forecasting of their effectively outworking.

Key words: hierarchical medium, gravitational field, acoustic field, geomechanical field, 2D iteration modeling algorithm.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

AUTHORS

Khachay O.A, Doctor of Physical and Mathematical Sciences, Leading Researcher, Senior Researcher, e-mail: olgakhachay@yandex.ru, Institute of Geophysics of Ural Branch of Russian Academy of Sciences, 620016, Ekaterinburg, Russia,

Khachay O.Yu.1, Candidate of Physical and Mathematical Sciences, Assistant Professor, KhachayA.Yu.1, Candidate of Physical and Mathematical Sciences, Assistant Professor, 1 Ural Federal University, 620002, Ekaterinburg, Russia.

ACKNOWLEDGEMENTS

This study has been supported by the Program of Governmental Support for Leading Schools, Project No. NSh-9356.2016.1.

REFERENCES

1. Dmitrievskiy A. N. Izbrannye trudy. T. 2 (Selected works, vol. 2), Moscow, Nauka, 2009, 435 p.

2. Khachay O. A., Khachay A. Yu. Vestnik YuURGU. Seriya «Komp'yuternye tekhnologii, upravlenie, radioelektronika». 2011, no 2(219), pp. 49—56.

3. Khachay O. A., Khachay A. Yu. Vestnik YuURGU. Ceriya «Vychislitel'naya matema-tika i informatika». 2013, vol. 2, no 2, pp. 49—56.

4. Khachay O. A., Khachay O. Yu., Khachay A. Yu. Geoinformatika. 2015, no 3, pp. 45-51.

5. Khachay O. A., Khachay O. Yu., Khachay A. Yu. Gornyy informatsionno-analitiches-kiy byulleten'. 2016, no 4, pp. 354-366.

COMPLEXING ACOUSTIC, GRAVITATIONAL AND GEOMECHANICAL FIELDS IN HIERARCHICAL MEDIA

Текст научной работы на тему «Комплексирование акустических, гравитационных и геомеханических полей в иерархических средах»