О комплексировании сейсмических и электромагнитных активных методов для картирования и мониторинга состояния двумерных неоднородностей в N-слойной среде Текст научной статьи по специальности «Физика»

УДК 550.830+539.3

О КОМПЛЕКСИРОВАНИИ СЕЙСМИЧЕСКИХ И ЭЛЕКТРОМАГНИТНЫХ АКТИВНЫХ МЕТОДОВ ДЛЯ КАРТИРОВАНИЯ И МОНИТОРИНГА СОСТОЯНИЯ ДВУМЕРНЫХ НЕОДНОРОДНОСТЕЙ В N-СЛОЙНОЙ СРЕДЕ

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

ABOUT COMPLEX APPROACH OF USING SEISMIC AND ELECTROMAGNETIC ACTIVE METHODS FOR MAPPING AND STATE MONITORING OF 2D HETEROGENEITIES IN A A/-LAYERED MEDIUM

O.A. Hachay, A.Yu. Hachay

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

Ключевые слова: комплексная сейсмическая и электромагнитная методика, алгоритмы решения прямой динамической задачи сейсмики, совместный анализ прямых двумерных задач для сейсмического и электромагнитного поля, двумерная однородная неоднородность в N-слойной среде.

That paper is devoted to analysis of solution 2D direct problem for seismic field in dynamical variant with use integral and integro-differential equations and common analysis with the solution of electromagnetic field for the case H-polarization for the model homogeneous inclusion in iV-layered medium.

Keywords: complex seismic and electromagnetic method, algorithms of solution of the direct dynamical seismic problem, common analysis of the solutions of seismic and electromagnetic field, 2-d homogeneous heterogeneity in a N-layered medium.

Введение

В связи с усложнением решаемых геологических задач все более актуальной становится проблема создания комплексных методик исследования. При комплексировании геофизических методов необходимо пространственное совмещение систем наблюдения и создание единых интерпретационных алгоритмов. В Институте геофизики УрО РАН разработана единая ЗО-методика наблюдений с использованием электромагнитных и сейсмических полей (в динамическом частотно-геометрическом варианте), опирающаяся на единую концепцию трехэтапной интерпретации [10-13], которая включает в себя единую систему векторных наблюдений сейсмических и электромагнит-

Хачай Ольга Александровна - д-р физ.-мат. наук, ведущий научный сотрудник, Институт геофизики УрО РАН; [email protected]

Хачай Андрей Юрьевич - канд. физ.-мат. наук, доцент математико-механического факультета, Уральский государственный университет имени А.М. Горького, г. Екатеринбург; [email protected]

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

Hachay Olga Alexandrovna - PhD, leading scientific expert of the Institute of geophysics of the Ural Branch of RAS; [email protected]

Hachay Andrey Vurievitch - PhD, assistant professor of the Mathematics and mechanics department of Ural State University of A.M. Gorky, Ekaterinburg; [email protected]

посредством преобразования Фурье (для действительной частоты) [И].

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

В работе [4] изложена идея алгоритма построения критерия подобия для сейсмических и электромагнитных полей от сингулярных источников, эквивалентных проявлению искомого локального объекта в рассматриваемых полях. При этом объект в сейсмическом поле аппроксимируется погруженной точечной горизонтально действующей силой в /'/-слойном упругом изотропном полупространстве, в электромагнитном поле это погруженный горизонтальный магнитный диполь в /У’-слойное изотропное проводящее полупространство. Выбор сингулярных источников определялся подобием морфологий создаваемых ими полей при взаимно перпендикулярном направлении моментов электромагнитного и сейсмического сингулярных источников. Разработан итерационный алгоритм вычисления электромагнитного и сейсмического полей при условии расположения аномальных источников в произвольном слое Л'-слойной среды [5-7].

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

1. Прямая динамическая задача сейсмики. Постановка задачи

Основное дифференциальное уравнение динамики упругого тела имеет вид [2]

\\Ли + (к + n)graddiv и + <jF0 = о

д2й

dt2

(1)

Здесь и(их,иу,и;) - вектор смещения; Р0(Р0х,Р0)„Р0:) -вектор массовых сил; X и ц - постоянные Ламэ; а - плотность среды; / - время.

Мы будем изучать сейсмическую задачу для случая установившихся колебаний [2]:

цАй+ (Х + р)^асИіу и + оР+оз~ аи = 0. (2)

Введем обозначения:

Х + Ц (У 2 I 2

: А н---------graddiv; —а> -к2

ц (Д

(3)

Тогда уравнение (2) с учетом (3) можно переписать в виде

А*й + к}й = —F.

(4)

Используя векторные преобразования [1], получим (А, + 2ц) graddiv it - \xrotrotii + а2аи = -аР. (5)

Представим и и F в виде:

и = grado? + roty ; F = gradO + rot'Y . (6)

Введем обозначения:

и, = grad<p; Ti2 - rot\\i. (6’)

Тогда

divü = diviii = Дф; rotii = rotu2 =

= graddivij> - Дф (7)

и

graddivu = grad Дер; rot rotii = -rotA\\i. (8)

Внесем в (5) выражения (6), (8):

(X + 2[i)gradAty + \x¡'oiA\\í + ю 2ogradty +

+ Gi2arot\\i = -aigrad^b + rot^V). (9)

Следовательно, чтобы удовлетворить (9), достаточно положить

(X + 2\i)gradA<p+ (ü2agradq> - -agradeE>;

p.rotA\¡i + <o2orot\\i = -aro 14*,

или

(к + 2ц) Дф + ю2аф = -оФ;

цЛф + ©2<j\j> = -a*í/.

Введем дополнительно обозначение:

(10)

со

■к2.

■ <П)

С учетом (11) уравнения (10) можно записать в виде:

Дф + &і ф = -

СТ -Ф; Avfr + к2ц ■

-Щ.

(12)

\ + 2ц ~ ¡а

Будем предполагать, что среда является упругой и выполняется известная линейная зависимость между компонентами тензора напряжений и деформаций [2]:

(7Х = Xdivu + 2¡i

у ^ух И

У*

V

дих ди —- + —-ду дх

ди

ст(, = Xdhm + 2\i——: ду

а_ = Xdivu + 2ц

ди,

~йГ

V = hv = Iа

(,3)

диу ди. dz ду

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

Для системы уравнений (12) в [2] сформулированы условия излучения в виде:

щ = 0(1),

ги2 =0(1), г \ —--¡к2й2

дг

дщ

дг

дщ

&\Щ =о(1);

(14)

: о(1).

2. Задача о дифракции звука на двумерной упругой неоднородности, расположенной в у-м слое 7У-слойной среды

Рассмотрим более общую задачу дифракции звука на двумерной упругой неоднородности, расположенной в у-м слое Л'-слойной среды. Вывод соответствующих уравнений для задачи дифракции звука на двумерной упругой неоднородности, расположенной в упругом пространстве, сделан в работе [2].

Эту задачу будем решать, используя подход, изложенный в работе [3]. Массовые силы будем считать потенциальными и сосредоточенными в первом слое Л'-слойной среды. Плоскость ХОУ совпадает с верхней плоскостью 1-го слоя, г = 0. Ось 02 направлена вертикально вниз. Образующие двумерной неоднородности в виде цилиндра произвольного сечения направлены вдоль оси 07. Положим ц = 0, тогда в каждом из слоев выполняется первое уравнение из (12), преобразованное к виду:

Аф, +^ф, =-2п/,(М)-

ст (15)

й = gradщ /,{М) =—^— Ф,

где / = 1,..., и; Ф,

2 пХ,

= Ф при г

1, при Ф, = 0.

Волновое число в г-м слое равно согласно (11): Введем обозначение:

(16)

к(М) = к]: - волновое число в слоистой среде;

/ =1, ...,N11

Г кх 7, Л/ е 50; К{М) = \.

I ЦМ),М«50;

ку, = ю

2

X,,

(17)

Индекс у7 обозначает свойства среды внутри неоднородности, находящейся в пространстве с поверхностью 50. В общем случае в произвольном слое или внутри неоднородности уравнение (15) с учетом (16) и (17) будет иметь вид:

Дф| +К2(М)^ =-2%/){М) . (18)

Граничные условия в среде без разрывов, заключающиеся в непрерывности вектора смещения и тензора напряжений согласно (6) и (13) на границах раздела имеют вид: дф, 5(Р,+1=0|

дг дг дп)ц I дп

[0;(ю2ф, +Ф,)]-[ст,+1((й2фт + Ф,+1)] = о|ге/. ; |^0(со2ф + Ф^ =|^а(со2ф + Ф)^ ; (19)

[а(ш2ф+Ф)]|г=0 =0,

где индекс у/ обозначает значения о, ф, Ф на границе неоднородности с внутренней стороны, ¡а -с внешней стороны границы неоднородности, которая расположена в у-м. слое; Ь - граница раздела слоя, с индексом / - со стороны /-го слоя, с индексом /+1 — со стороны /+1-го слоя. Условия затухания на бесконечности согласно (6) и (14) имеют вид:

rgradцi = 0(1),

Пусть

г(^--/^ф,.) = о(1). ог

(20)

л

ф, =Ф,-Ф-, (21)

где / = 1,...,/, ..., А'; - потенциал нормального

сейсмического поля в слоистой среде в отсутствие 0

-у,=ф/

(22)

неоднородности: ф“( = ф° и

Аф°+*~фг0=-2тс/;(м);

Эф° 5ф®+1

= 0

ге1, .

дг дг

ст,(ш2ф° +Ф;)]-[а,ч1(ю2ф‘)+1 +Ф/+])] = о|

(23)

ТТ о Зф

На контуре неоднородности ф и —— непрерыв-

дп

ны. ф, - потенциал аномального сейсмического поля, который, с учетом (18), (17) и (21) удовлетворяет уравнению

Дф, +К2(М) ф, = -(К2{М)-к2(М))^ и граничным условиям:

дЦ>, 5ф,+1 дг дг

= 0,

б/.; >

■ Ж

дп

дп )

№

[<т,со2ф, ] -[ст,+1а>2ф/+1 ] = оЦ . На контуре неоднородности [аф] -[оф] =(ауо-а/7)ф°.

(24)

(25)

(26)

Функция источника сейсмического поля С/Кр (М, М )

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

+к2С$р,(М,М°) = -2715(М - М°) и граничным условием:

(27)

= 0

]и=о =о;

дг дг

] - [сты®2<%,/+1 ] = о| ■

На контуре неоднородности Ох (М,М°)

(28)

дОцр(М,М )

./«

дп

непрерывны. Применим формулу

Грина [3] для функций ф, и С-Ар(М,М )для каждого слоя /У-слойной среды при / Фу

дп

-С,

5р,і

дп

й1,-

дС

Ач]

дп

-а,

лЧр,і+1

5фм

дп

<а,

/+1

271фДМ0), А/0 є 5,;

Я*.

ц,

дО.

3>.У

дп

-с?»

$р,і

д<? } дп

Щг

дп

'« 90С

фуо—Г—-Ой

-б,

9«

^+1 +

Эи

>> '

дфу0

дп

Ща =

¡27Гф,(Л/0), М° є 5, ;

(29)

(30)

[ о, ми г 5,-.

Помножим выражения (29) и (30) на а, и сложим их с учетом граничных условий (28), в результате чего получим:

Ні

2п

$(Ф.,а-

дс,

Яр,]

Эф

^р,} '

¡О

дп

)с1с ■

[о(М°)ф(Л/°), М° є -5С; [ 0, М° £Яг5с.

(31)

Применим формулу Грина для функций ф,7 и

^_____Л

2%

с

а-/'г)фс^/-^-^ = а(М°)ср(М0). (33)

2тс

с;

дп

Воспользуемся равенством [3]:

(?Ю -<*ц)

№______Р

2п

Эн

Эф°(М)Л дп

йс =

(а,а-а,)ф°(М0), АГє5с;

Л

0, ЛГ £8г.

(34)

Тогда выражение (33) с учетом (34) можно переписать в виде:

271

(СТ/о ~°>)

иФ(мх?ф;і(м,м°к

а

тм+-

а,-,

уа о

Ф

'(М0)-

§С5ау.^ = ф(М0), М° є5с;

,а„2тг

дп

(35)

(СТ/0~СТ/')|ОлА/^ = ф(М0), М°«5С.

а(М°)2я с

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

СНр! (М,м ) для внутренности области 5С, с учетом (25), (29), (31), (33) получим:

к2 -к2 2п

\\^муз3р^м,м°)йх

)с1с

М

1 дС8р,і

дп

Щл,

дп

[ф(Л/°), М° є5с;

(32)

[ 0, М"г5с.

Помножим выражение (32) на сг^, и сложим с

(32), с учетом граничных условий на контуре С неоднородности получим:

-*5)

2 71

Яф {М)С8р]{М,М°)ёхм-

2 п

дО,„

+^-^ффО(М)^^С-

дп

3. Сопоставление алгоритмов для сейсмического (дифракция звука) и электромагнитного двумерного моделирования (Н-поляризация)

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

ЯХ(М°) = -

кі-кі

2п

\\нх (М)вт (М, М° )ЛМ

к — к КМ к!

к2 271 :с

§НХ(М)

эс,

]Л

І' 1/0/ я^Ол

сіс + ^Н^М"), М° є Бс\ дп к

Нх(МуЗт(М,А^)с1-

2пкц

(36)

'%/ +

к2(М°Хк2-ф

к2к22п

Нх(М)^-^+Нх(М°), М°£3С. дп

Здесь к2(М°) = 10щ()а(М0), Цд = 4тг-10 7 —, о(М°) -

м

проводимость в точке А/0; / - мнимая единица; НХ(М°) - суммарная составляющая магнитного

поля; Н^(М°) - составляющая магнитного поля слоистой среды в отсутствие неоднородности; А2 = т\10<у^, к2 = /(оц0ст(, о., - проводимость

внутри неоднородности, расположенной в у-м слое; сг, - проводимость /-го слоя Л'-слойной среды. Функция Грина Ст (М, М°) в выражениях

(36) удовлетворяет уравнению аналогичному (27) с ~2 2

& {М) = к1 , введенному выше. Граничные условия состоят в непрерывности самой функции на границах раздела слоев, однако нормальная производная ее терпит разрыв вида

1 дС,

ті+І

7+1

дп

(37)

Сопоставим полученные уравнения для сейсмической и электромагнитной задачи:

^М°) = Л ф(М)С^(М,М°)с1тм +

2л

у +^цР(М0),М0 є 5Г; ау,2тг *. дп Ор

к2-

(38)

дС].

к2 2% с

ЦЬнх{М)-^(к+-4гНйх{М0),М0є8с.

■ 2 ' дп ъ1

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

4. Алгоритм моделирования дифракции линейно поляризованной упругой поперечной волны на двумерной неоднородности, расположенной ву-м слое Л'-слойной среды

Рассмотрим задачу о дифракции линейно поляризованной упругой поперечной волны на дву-

мерной неоднородности, расположенной в у-м слое Л'-слойной среды, используя подход, описанный в работе [3] для электромагнитной двумерной задачи (случай Н-поляризации). Геометрическая модель среды аналогична описанной выше в предыдущей задаче. В этом случае в уравнении (5) Ру, Р. равны нулю,

и. = иу = 0, (39)

а составляющая вдоль оси ОХ не зависит от х. Составляющие тензора напряжения (13):

Т’ХУ ^ д, ’ ’в г* ? 5 -у-

су &

Тогда наша задача с учетом (8), сводится к решению следующей задачи:

Аих1+к1их,=-2п/,(МУ, /,(М)=^-РХГ, (41)

где /= 1, ...,я; FЛ■ = F при /= 1, при іф 1 Рхі = 0. Волновое число в і-м слое равно согласно (3): ст,

х2і

= ю-

(42)

Введем обозначение:

к(М) = к2, - волновое число в слоистой среде, /= 1, ..., пи

( > А/ е З'о;

К(М) = <

[к(М), М £ 80;

(43)

=ш2^

Индекс у/, как и в предыдущей задаче указывает на свойства среды внутри неоднородности. В общем случае в произвольном слое или внутри неоднородности уравнение (41) с учетом (42) и (43) будет иметь вид

Л«,, +К2{М)иХ1 = -271иМ). (44)

Граничные условия, заключающиеся в непрерывности вектора смещения и тензора напряжений на всех границах раздела, сводятся к непрерывности иХ1 и (.1, ——. Выполнены условия затухания дп

на бесконечности.

Пусть

(45)

. .... о

где I - 1, ...,у,у?, ..., п; их, - составляющая нормального сейсмического поля в слоистой среде в отсутствие неоднородности: М®, = Иу и

Д и° +к1,и% =-2тг/(М);

И,"

дг4,

дг

-Ц/+Г

сгіхіі+1

52

= 0

гєД

.0

(46)

(47)

ц о дих

На контуре неоднородности их и —- непре-

дп

рывны. II,

составляющая аномального сеисми-

ческого поля, которая, как легко показать аналогично (24), удовлетворяет уравнению

МХ1 +К\М)йх, =-(К2(М)-к\М))й1

и граничным условиям:

дйг і дйг

2п

і

ас;,.

дп

--о,*

5Й..

дп

&]о =

(48)

\\х(М°)йх(Мй), М° є Я

Ц^=Цу. (54)

—Ц/+Г Йт : -ЇІ

дг

:0

х,і их,і+1 иЬєі,- ’ Чї|г=0

дйх

дп

дп

ди

(49)

'¡а

дп

йх,И~йх,}а =0.

Функцию источника сейсмического поля с;л« (М,Л/°) выберем следующим образом: она определяется как решение следующей краевой задачи:

+РсЯх,(М,М°) = -2п5 (М-М°) (50)

и удовлетворяет граничным условиям:

9С?&,-+ 1 Л|

И/ —т— ~ ^+1 —:-------= 01"

[ 0, Миг5у-5С;

Применим формулу Грина для функций йх и ,М<}) для внутренности области Бс, с учетом (44), (46), (48), (50) получим:

к2 -к2 К2р К2]

2%

\\их(М)(],:лі(М,М°)сіхм

,_Ц[Й^

2п / *р дп

г дй*,Л

(~Гс

дп дп

СІСі

\йх(М°1 М° е5с;

(55)

дг

дг

г&ц

(51)

На контуре неоднородности С&(М,М ) и

[ 0, Ми

Помножим выражение (55) на |1/( и сложим с (54),

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

уи(му?& .(М,м°ухм -

дС8,(М, М ) дп

непрерывны. Применим формулу (и /и-Ц7,) г о 96;

Грина [3] для функций м„ и 0&(М,М ) для каждого слоя |У-слойной среды при ¿Фу

са,-

)Ч______гр

2п

с

и

дп

■'Зя.і

дп

..¿гАхд^ыс1с+

2тс £ дп

дв

дп

С

&,/+!'

"*,/+1 дп

2%

§] -= а(М°)йх(М°). (56)

еВ,

/+1

Воспользуемся равенством [2]

\2пйхі(М°), М° є 5,;

0, М° г 5,-.

(52)

(1-У -ц^)

2я

!#

_ 9и°(М)'

дп

дп

При / =у

дп

9м

9и

(II,

|(^-о-ц7,К0(М0), М°е5с; I 0, М°г5г.

а?с =

(57)

I -сад,-

і

+ ф и

^х,7+1 Л

дп

с!!/+, +

Тогда выражение (56) с учетом (57) можно переписать в виде:

(4> -*22/)

¿}±__У_У 2п

Цих(М)С^(М,^хм +^и°х(М°) Н/'

+

•л,7«

9«

С

9м

& ,у ■

х,;а

дп

Асіа =

= моє5с.

дп

12ямх>у(Л/°), М° є Л’у - 5 с; [ 0, М°ї57.-5с.

(53)

^,-2я с

^/>(^2 н —^2/)

(58)

Помножим выражения (52) и (53) на ц, и сложим их с учетом граничных условий (49) и (51), в результате чего получим:

ЦъЮв*, ДМ, М°Утм+м«(М°)+

+ 0*/г^г £^=

ц(Л/))27Т I дп х

В выражениях (58) выписан алгоритм моделирования сейсмического поля для поперечных упругих колебаний в Л'-слойной среде, содержащей двумерную неоднородность. Первое выражение есть интегральное нагруженное уравнение Фредгольма второго рода, решением которого является распределение составляющей вектора упругих смещений для поперечных колебаний внутри неоднородности. Второе выражение есть интегральное представление для вычисления вектора упругих смещений в любом слое Л'-слойной среды.

5. Сопоставление алгоритмов для сейсмического (линейно поляризованная поперечная упругая волна) и электромагнитного двумерного моделирования (Н-поляризация)

Сопоставляя выражения (58) с соответствующими выражениями для электромагнитного поля (Н-поляризация) (36) мы видим наличие подобия интегральных структур этих выражений. Различие имеется в коэффициентах при соответствующих членах в выражениях (36) и (58), что можно учесть в выборе системы наблюдения того или иного поля. Следует учесть также различие в частотной зависимости отклика среды на сейсмическое и электромагнитное возбуждения. Однако при соблюдении коэффициентов подобия сейсмическое поле, возбуждаемое поперечными колебаниями, и электромагнитное поле будут содержать подобную информацию о структуре неоднородной среды и связанном с ней состоянии. Эти результаты подтверждаются натурными экспериментами, изложенными в работах [9, 12].

Заключение

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

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

Литература

1. Кочин, Н.Е. Векторное исчисление и начала тензорного исчисления / Н.Е. Кочин. — М.: Изд-во АН СССР, 1951.-426 с.

2. Купрадзе, В Д. Граничные задачи теории

колебаний и интегральные уравнения / В.Д. Купрадзе. - М.; Л.: Гос. изд-во технико-

теоретической литературы, 1950. — 280 с.

3. Дмитриев, В.И. Дифракция плоского электромагнитного поля на цилиндрических телах, расположенных в слоистых средах / В.И. Дмитриев // Вычислительные методы и программирование в слоистых средах. — М.: Изд-во МГУ, 1965. -Вып. III. - С. 307-316.

4. Хачай, А.Ю. Изучение критерия подобия

сейсмических и электромагнитных полей от погруженных сингулярных источников для осуществления активного мониторинга нестационарной среды / А.Ю. Хачай // Глубинное строение, геодинамика, мониторинг. Тепловое поле Земли. Интерпретация геофизических полей: материалы

Третьих науч. чтений Ю.П. Булашевича. — Институт геофизики УрО РАН, 2005. — С. 143—145.

5. Хачай, А.Ю. Алгоритм решения прямой задачи электромагнитных исследований при возбуждении горизонтальным магнитным диполем, расположенным в произвольном слое п-слойной изотропной проводящей среды / А.Ю. Хачай // Информатика и математическое моделирование. -Екатеринбург, 2006. - С. 136-168.

6. Хачай, А.Ю. Алгоритм решения прямой динамической задачи сейсмики при возбуждении горизонтальной точечной силой, расположенной в произвольном слое N-слойной упругой изотропной среды / А.Ю. Хачай //Информатика и математическое моделирование. — Екатеринбург, 2006. — С. 170-278.

7. Хачай, А.Ю. Алгоритм решения прямой динамической задачи сейсмики при возбуждении точечным источником вертикальной силы, расположенной в произвольном слое N-слойной упругой изотропной среды / А.Ю. Хачай // Информатика и математическое моделирование. — Екатеринбург, 2006.-С. 279-310.

8. Использование комплексной попланшетной сейсмической и электромагнитной методик для решения задач картирования приповерхностных неоднородностей / O.A. Хачай, B.C. Дружинин, Ю.С. Каретин и др. // Геофизика XXI столетия: 2001 г.: сб. тр. третьих геофизических чтений им.

B.В. Федынского. — М.: Научный Мир, 2001. -

C. 327-337.

9. Хачай, O.A. Изучение критерия подобия для сейсмических и электромагнитных исследований в

частотно-геометрическом варианте / O.A. Хачай, Т.А. Хинкина, В. В. Бодин // Астрономо-геодези-ческие исследования. — Екатеринбург: УрГУ. -2001. - С. 30-35.

10. Трехмерный электромагнитный мониторинг состояния массива горных пород / O.A. Хачай, Н.П. Влох, Е.Н. Новгородова и др. // Физика Земли. - 2000. -№12.- С. 1-8.

11. Хачай, O.A. Предпосылки сейсмоэлектро-магнитного мониторинга нестационарной среды / O.A. Хачай, Т.А. Хинкина, В.В. Бодин // Российский геофизический журнал. - 2000. —№ 17-18. - С. 83-89.

12. Метод картирования зон потенциальной неустойчивости массива горных пород различного

вещественного состава с использованием данных динамической сейсмики и электромагнитных индукционных исследований / O.A. Хачай, В.В. Бодин, Е.Н. Новгородова и др. // Горный информационно-аналитический бюллетень. —2001. —Ms 3. —

С. 10-16.

13. Хачай, O.A. Исследование разрешающей способности попланшетной электромагнитной методики для активного картирования и мониторинга неоднородных геоэлектрических сред / O.A. Хачай, Е.Н. Новгородова, А.Ю. Хачай//Физика Земли. - 2003. —№ 1. — С. 30-41.

Поступила в редакцию 13 февраля 2010 г.