Файл: Романов В.Г. Обратные задачи для дифференциальных уравнений спецкурс для студентов НГУ.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.