Файл: Романов В.Г. Обратные задачи для дифференциальных уравнений спецкурс для студентов НГУ.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 19.06.2024

Просмотров: 207

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

тыо

tffx)

, зависящей от точки пространства;

в точках

границы

 

приборы фиксируют время пробега

т(х,х")

волн от источника

сх°б 3

до приемника

х е 5 ;•. требуется по временам

т(х,х?),

где

х, х°-

произвольные

точки поверхности S,

найти скорость

V(x)

распространения волн внутри области *SD . Б дальнейшем мы

увидим, что функция

т/Ьс, ccj

удовлетворяет

внутри области

Ф

уравнению первого порядка

 

 

 

 

 

 

 

 

 

 

 

 

\^uudx

Т(х,х°) |=/г(сс),

n-W^^fr),

 

О )

 

причем

т(х°х°)=о1

поэтому в математическом отношении эту за­

дачу можно сформулировать как задачу

отыскания такого дифференци­

ального

уравнения (9)

(то есть

такой

функции

п(х)},

 

 

что

его

решения,

удовлетворяющие

условию

Т(х°х°)

= 0,

принимают

при

произвольных

х, х°е

S

заданные

значения. Точка

х°

является

здесь параметром прямой задачи для уравнения ( 9 ) .

 

 

 

 

 

Своим происхождением обратная кинематическая задача

обязана

геофизике. В 1905-1907 годах Г.Герглотц и Е.Вихерт

рассмотрели

первую обратную кинематическую задачу для сферически симметрич­ ной модели Земли, в предположении, что скорость распространения

возмущений монотонно растет

с глубиной. Ими было показано, что

скорость передачи возмущений

v(x), являющейся при этом предпо­

ложении функцией расстояния до центра Земли, однозначно определя­

ется функцией

т(х, х°)

при фиксированном

значении х ° е б

и

произвольной точке xeS

 

(в данном случае

3 - поверхность

зем­

ного шара), получены явные

формулы для отыскания скорости ~0(х).

Это позволило,

основываясь

на данных сейсмических наблюдений над

землетрясениями, сделать первые выводы о глубинном строении Земли.

Землетрясения, происходящие внутри Земли, вызывают упругие волны двух типов, (см. [*Э2] ) . Волны, идущие с большими скоростя­ ми - волны сжатия и разрежения; они состоят из колебаний, совпа­

дающих по своему

направлению с направлением их распространения.

Эти волны носят

название

продольных. Волны второго типа связаны

с колебаниями типа сдвига

или кручения, перпендикулярными к на­

правлению распространения

колебаний, поэтому они называются по­

перечными. Скорость распространения продольных и поперечных волн связана с упругими свойствами вещества, через которое они прохо­ дят:



В этих формулах Л,

^и-

параметры Ламе, характеризующие сопротив­

ляемость вещества на

сдавливание и на сдвиг, соответственно, JJ -

плотность вещества, Up ,

ttg - скорости распространения

продольной

и поперечной

волн.

 

 

 

 

Упругое

вещество полностью характеризуется параметрами Я ,^и,

р . Поэтому,

если бы скорости распространения волн

vp

, -г& были

известны, соотношения (10) давали бы уже достаточно

большую инфор­

мацию о веществе Земли. В связи с этим одна из основных задач

сейсмологии

и заключается

в отыскании скоростей распространений

продольной и поперечной волн. Для ее решения используются наблю­ дения над режимом колебаний земной поверхности, проводимые на сейсмостанциях. Энергия некоторых землетрясений настолько велика,что результат их воздействия, может быть зафиксирован в любой точке земного шара. Процесс распространения упругих колебаний внутри Земли описывается системой динамических уравнений теории упругости

 

+ £

Ш +

?

w

 

=

р

 

Здесь

•гГ = (id,

itlt v3)

-

вектор

смещений

точек среды

по

осям

ха > Tj, х3

соответственно, Те

-

единичный

орт оси хе .

Ана­

лиз решении уравнений ( I I )

приводит

к выводу,

что продольные

и

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

до

сейсмостанций (в

сейсмологии

такие

таблицы носят

название г о ­

дографов волн [ 33] ).

 

 

 

 

 

Если обозначить

через Г(х,

х°)

- луч, соединяющий

точки

X, ас", а через ds -

элемент длины его дуги в точке х

, то

время

dt

, затрачиваемое

на прохождение пути d$ волной,

имеющей

ско-

10


рость v(x) , можно подсчитать по формуле

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

(12)

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

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

рых зависит скорость) модели имели очень большое значение,

так

как благодаря им удалось сделать некоторые выводы о

структуре

Зем­

ли. Вместе с тем они стимулировали новые постановки

обратных

за­

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

По современным представлениям земной шар состоит из трех ос­ новных частей: коры, оболочки и ядра. Земная кора представляет из себя сравнительно небольшой по мощности слой. Толщина земной коры в среднем составляет 30-50 км, она меньше под океанами и

I I

больше под континентами, особенно увеличиваясь в районе горных массивов. Считается, что кора наиболее резко различается по сво­ им скоростнш» свойствам. Здесь налицо довольно существенное раз­ личие в скоростях распространения волн от одного района к друго­ му. Кора отделяется от оболочки поверхностью Мохо, на которой ста­ рость распространения волн имеет разрыв. Оболочка Земли (употреб­

ляется также термин "мантия Земли") имеет

толщину порядка 3000

км

и состоит из твердого аморфного вещества.

Земное ядро обладает

 

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

В настоящее время рядом геофизиков построены одномерные сфе­ рически-симметричные модели распределения скоростей сейсмических

волн внутри Земли (см.,

например, [26,39,50,53j

) . Поотроены они

на основании различного

сейсмологического материала, при интер­

претации его применялись различные приемы, поэтому естественны некоторые различия в этих моделях. Интересно, однако, отметить, что все модели отличаются друг от друга в среднем на величину по­ рядка 10%, то есть достаточно близки. Вместе с тем в геофизике имеются в настоящее время факты, позволяющие предполагать о су­ ществовании внутри Земли неоднородности! по географическим коор­ динатам. К числу этих фактов относятся систематические отклонения в годографах волн, построенных для различных материков, от усред­ ненного годографа, асимметрия гравитационного и магнитного полей. При этом отклонения от годографов, отвечающих сферически-симмет­ ричному распределению скоростей упругих волн, достаточно малы.На­ пример, при временах пробега сейомических волн от источника воз­ мущения до приемника порядка 10-20 минут отклонения в годографах не превышают Ь-10 секунд. Это позволяет предположить, что откло­ нения скоростного строения Земли от сферически-симметричной моде­ ли также малы. Тем не менее эти малые флуктуации в распределении скоростей волн представляют для геофизики большой интерес, так как, возможно, они помогут объяснить механизм развития земной коры,воп­ рос дрейфа материков и т . д . Таким образом, задача определения трехмерного скоростного строения Земли является одной из актуаль­

ных задач геофизики.. В математическом отношении эта задача

связа­

на с решением обратной задачи для

уравнения (

9 ) , когда tf(x)

яв­

ляется произвольной функцией трех

переменных.

К сожалению, полной

12


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

пока нет, здесь имеются лишь отдельные частные результаты.

Нет

также и достаточно обоснованных методов решения неодномерной

об­

ратной задачи, имеются,однако, попытки численного решения

этой

задачи (см. £ 4,5,I9j) .

 

В последнее время появились довольно обоснованные сомнения в том, что скорость в мантии монотонно растет с глубиной. Рядом гео­ физиков высказано предположение о существовании в мантии Земли слоев с пониненной скоростью или, как их обычно называют, волно­ водов. При наличии последних обратная кинематическая задача с ис­ точниками возмущений, расположенными выше волновода, становится неоднозначной. Характер этой неоднозначности был исследован М.Л.Гервером и В.М.Маркушевичем [45,47] . 3 рамках кинематики един­ ственность задачи можно сохранить только располагая источники воз­ мущений в самом слое пониженной скорости f I I ,477 . Для некоторых районов земного шара постановка обратной кинематической задачи с источниками, расположенными внутри мантии, имеет практический смысл. Так, например, в поясе, проходящем через Курильские остро­ ва и Японию, землетрясения, происшедшие за последнее столетие, заполняют целый объем, простирающийся в глубину до 600 км. В це­ лом по земному шару землетрясения располагаются,в основном, в ко­ ре. Поэтому возможное наличие волноводов существенно ограничивает применимость кинематического подхода при изучении мантии Земли.а связи с этим все большее внимание исследователей привлекают урав­

нения динамической теории упругости ( I I ) . Дело в

том,

что

обрат­

ную задачу, заключающуюся в определении

скоростей

Vp,

V*,

можно

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

( I I ) . Имеющиеся на

поверх­

ности Земли сейсмосганции фиксируют от землетрясений полные сейс­

мограммы, то

есть

измеряют смещения tl(x,x°i)

как функции вре­

мени. Здесь,

как

и раньше,

х " - точка приложения

сосредоточенного

воздействия,

X - точка, в которой находится приемник. Располагая

сейсмограммами от достаточно

большого числа землетрясении

для

различных районов земного шара и достаточно густой сетью сейсмо­ сганции можно пытаться определить по ним функции iff, , как функ­ ции географических координат и глубины. Математическая теория по­ следней обратной задачи, которую естественно назвать обратной ди­ намической задачей теории упругости, только начинает создаваться. Имеются только отдельные результаты, относящиеся к чаотной одно­ мерной модели Земли (см._ [2,3,21,28,29,44,48]).

13.