CC BY
3УДК 550.344.6+ 550. 834
Сайипбекова А.М.1, Абдрахматов К.Е.2, Степаненко Н.П.3, Макамбаева Ж.А.1 , Саипбекова С.Э.1, Жыргалбекова Э.Ж.1
1Ошский государственный университет, г. Ош, Кыргызстан 2Институт сейсмологии НАН КР, г. Бишкек Кыргызстан 3ТОО «Национальный научный центр сейсмологических наблюдений
и исследований», г.Алматы, Казахстан
АНАЛИЗ ФИЗИКО - МАТЕМАТИЧЕСКИХ ОСНОВ МЕТОДОВ ОПРЕДЕЛЕНИЯ ДВУМЕРНЫХ СКОРОСТНЫХ МОДЕЛЕЙ ПО
ГЕОТРАВЕРСАМ
Аннотация. В статье проанализированы математические основы профильной сейсмотомографии. Изложены физико-математические основы распространения рефрагированных волн. Лучевое трассирование велось с заданием граничных условий. Определено, что двумерные сейсмические разрезы, построенные с применением детальных наблюдений по системам встречных и нагоняющих годографов, по сравнению со значительно осреднёнными скоростными значениями Vp наиболее информативны. При анализе волновой картины землетрясений до расстояния 350-780 км отмечена весьма сложная картина волнового поля. Определена двумерная скоростная модель по данным ГСЗ и региональных сейсмических станций Института сейсмологии НАН КР, ТОО « Национального научного центра сейсмологических наблюдений и исследований» МЧС РК и КНР.
Ключевые слова: физико-математические основы, анализ волновой картины, измерение волн, время пробега, скорость волн, профильная томография.
ПРОФИЛДЕРДИН ЭКИ eЛЧeМДYY ЫЛДАМДЫКТАРЫН АНЫКТОО МЕТОДДОРУНУН ФИЗИКА-МАТЕМАТИКАЛЫК НЕГИЗИН ТАЛДОО
Кыскача мазмуну. Макалада профилдик сейсмикалык томографиянын математикалык негиздери талданат. Сынып кийинки талаада таралган толкундардын физикалык жана математикалык негиздери кeрсeтYлгeн. Толкун таралуу жолундагы берилген чек ара шарттары кeрсeтYлгeн. Vp ылдамдыктын олуттуу орточо маанилерине салыштырмалуу, Карама каршы, жана орточо кармап туруучу негизги годографтар системаларынын жардамы менен эки тараптан байкоолордун жардамы менен курулган эки eлчeмдYY сейсмикалык участоктор эц так экендиги аныкталган. 350-780 км аралыкка чейинки жер титирeeлeрдYн толкундундардын тYЗYЛYШYн талдоодо толкун талаасынын eтe татаал экендиги белгиленген. Литосфердагы сейсмикалык толкундун ылдамдыгынын эки eлчeмдYY модели тундук Казахстан аймагында терен сейсмикалык изилдоо методу менен ал эми Тянь-Шань, Памир тарабы болсо Кыргызстандын сейсмология институтунун, Казакстандын сейсмологиялык байкоолор жана изилдeeлeр улуттук илимий борборунун, жана Кытай Эл Республикасынын аймактык сейсмикалык станцияларынын маалыматтарынын негизинде аныкталган.
Негизги сездер: физикалык-математикалык негиздери, толкундун тYPYн аныктоо, толкундун келген убакытын eлчee, толкундун ылдамдыгы, профилдик томография.
ANALYSIS OF PHYSICAL AND MATHEMATICAL PRINCIPLES OF METHODS FOR DETERMINITION OF TWO-DIMENSIONAL VELOCITY MODELS BY
GEOTRAVERS
Abstract. The paper analyzes the mathematical principles of profile seismic tomography. The physical and mathematical principles of the propagation of refracted waves are presented. Ray tracing was carried out with setting boundary conditions. It has been determined that two-dimensional seismic profiles created using detailed observations of systems of oncoming and catch-up time-travel curves in comparison with significantly averaged velocity Vp values are the most informative. Under analysis of the wave pattern of earthquakes up to a distance of 350-780 km, a very complex pattern of the wave field was noted. A two-dimensional velocity model was determined based on data from the State Survey and regional seismic stations of the Institute of Seismology NAS KR, the National Scientific Center for Seismological Observations and Research of MES RK and China.
Keywords: physical and mathematical principles, wave pattern analysis, wave measurement, travel time, wave velocity, profile tomography.
В определении и уточнении контуров сейсмоопасных зон территории Кыргызского Тянь -Шаня важную роль играют надёжные количественные сведения о неоднородностях литосферы, точнее сведения о горизонтальных и вертикальных неоднородностях земной коры и верхней мантии, потому что именно они являются индикаторами движений или же «изменения особенностей» глубинных процессов. Особенно, если учесть что Тянь-Шанский ороген -это система сложнопостроенных горстков и гребенов. В рельефе всё это отражается в виде высокогорных хребтов и глубоких межгорных впадин. По данным известного узбекстанского геофизика Борис Борисовича Тальвирского [12], общая амплитуда движений, приведших к образованию этих структур, составляет до 10 -12 километров. Всем известно, что при изучении неоднородностей ведущее место занимают сейсмические методы. В развитии представлений о скоростной модели земной коры и верхней мантии, безусловно, большую роль играют точные экспериментальные сейсмические материалы по отдельным отработанным профилям. Например, при повторной интерпретации отдельных материалов по профилю Балхаш (Казахстан)- Бачу (Китай) -продолжение Иссыкского профиля ГСЗ, рассчитанные математические модели показали, что рефрагированные волны дают богатую информацию о вертикальных градиентах (с отдельных градиентах скоростей в отдельных слоях) [11].
Прежде, чем проанализировать математические основы профильной сейсмотомографии, более досконально остановимся на волновом поле, которое при обработке землетрясений вдоль трансектов по региональной сети имеет большое значение. При обработке профильных наблюдений анализ волнового поля, выделение опорных, доминирующих волн является важным этапом. Для построения годографов представляют интерес участвующие в практике сейсмологической обработки местных землетрясений волны P, Pn, Pg, S, Sn и Lg в интервале эпицентральных расстояний 0- 400 км, а также закритические отражённые волны от границы Мохоровичича (М) (подошвы коры) Pотр., Sотр. в интервале эпицентральных расстояний 80-200 км, не используемые в рутинной сейсмологической обработке. Основной информацией являются вступления прямых P- и S- волн, взятые с бюллетений землетрясений в пределах эпицентральных расстояний 0-150 км.
Предпосылки, обеспечивающие решение задачи:
1). Высокая плотность очагов местных землетрясений.
2). Достаточно густая сеть региональных сейсмических станций.
Решения двумерной инверсии по данным рефрагированных волн. Двумерные сейсмические разрезы, построенные с применением детальных наблюдений по системам встречных и нагоняющих годографов, по сравнению со значительно осреднёнными скоростными значениями Vp наиболее информативны.
Линеаризованная постановка задачи отыскания распределения скорости по годографам рефрагированных волн предложена в работах [1] для сферической системы координат. Для двумерного случая предполагается, что распределение скорости может быть представлено суммой двух функций, одна из которых зависит только от глубины Z, вторая - от координат (X, Z):
кЪ(1)
Если время пробега рефрагированной волны в среде со скоростью V(X, 2) записать в виде
I (х,1) = Т (/ )+(х, /), (2)
где Т (1) - время пробега волны в среде со скоростью V и предположить, что величины и gradV малы, так что квадратами их величин можно пренебречь, то для времени справедливо выражение:
т(х, /) = г . (3)
Интегрирование ведётся вдоль траектории сейсмического луча, распространяющегося в среде с одномерным распределением скорости. Годографы задаются на дуге большого круга Земли и используется сферическая система координат. Решение даётся для поля времен t (х, 1). Важной особенностью принятого способа решения этой задачи является возможность одновременного использования всей системы встречных и нагоняющих годографов при получении двумерного поля значений истинной скорости в плоскости разреза. Для этого наблюдённые годографы преобразуются в двумерное поле времён в функции координат (х/2, Т). При трансформации наблюдённых годографов рефрагированных волн в поле значений истинных скоростей в плоскости разреза, необходимо учитывать следующие факты, например, алгоритм и программа рассчитаны на среды, не содержащие скачкообразного изменения скорости (граница первого рода), что связано с неизбежным искажением наблюдённых годографов в области выхода в первые вступления преломлённых волн (фактически происходит не учёт части петли рефрагированной волны, связанный с наличием в среде границы I и II рода и прослеживаемой во вторых вступлениях). Это всегда приводит к сглаживанию значений скорости вблизи сейсмических границ, что необходимо помнить при подготовке исходной структурно-скоростной модели к её корректировке через итеративное решение прямой задачи.
Для построения разрезов земной коры наблюдаемое волновое поле представляется в виде непрерывного годографа рефрагированных волн с петлями, сопровождающими области смен волн в первых вступлениях. Петли представлены элементами обычно выделяемые при ГСЗ волн: ветвями сменяемой и сменяющейся в первых вступлениях волн, а также элементом, кинематически совпадающим с закритически отражённой волной на всём интервале её прослеживания от «начальной точки» до асимптотического сближения с «прямой» волной.
Зададим годограф в двух видах и приведём основные формулы: Первый вариант. Допустим, что годограф является непрерывной функцией t (х), тогда для каждой точки годографа можно вычислить:
й 1 / ч
— = — = а(х) йх V
Одновременно считается известная обратная функция X(а). Тогда глубина проникания луча вычисляется по формуле Вихерта-Чибисова:
г=-1
тг *
V - X(а)йа
п [^(а-— )2 -1
Если считать годограф рефрагированной волны параболой второй степени, тогда: , = dx - Ьх2 (4)
Интервал АХ осреднённой параболы в точке с абсциссой Х из условия, чтобы «стрела перегиба» параболы не превосходила среднеквадратичную ошибку измерения времени 1 тогда:
в,=^х1 Ь,
4
где Ь - коэффициент. Его можно определить из уравнения (4) путём дифференцирования: dt 1
— = — = а - 2Ьх, dx V
где V - скорость в точке с абсциссой Х. Найдем - Ь: 1
V х V ,
Таким образом, расстояния между точками будет имеет вид:
ь = 1
х
Ах = 2 х
—5,
\5, - х
Ошибка в определении скорости А— связана с Ах, 5, и плотностью наблюдений д.
2-2 I 3(—-Ах -1) А— =- -————-^5,,
АхшП V —■Ах—ш (—-Ах—ш + 1)
где — - количество наблюдений на единице длины годографа.
При очень большой плотности наблюдений можно приближённо написать:
2—
Ах ■ = ,
1—1 А—\
3—2 (в,)2 - А—
(5)
2—- —
При интерпретации полей времени t (х, 1) отсчётные точки автоматически задаются сечением поля А1, при выборе которого по возможности должны быть учтены условия (4) и (5).
Второй вариант. Годографы задаются на дуге большого круга Земли и используется сферическая система координат. Решение ищется для поля времен t (х, 1) рефрагированных волн.
Поле времен t (х, 1) фактически представляет собой способ изображения наблюдаемых времён на плоскости (х, 1). Параметром поля является расстояние (база зондирования) точки регистрации от источника. Время относится к центру соответствующей базы. Соединяя на плоскости (х, 1) времена й, отвечающие фиксированным значениям баз 1к, непрерывными линиями, получим семейство изолиний 1= сом^ которые образуют поле времен 1 (х, 1).
Поле времён характеризуется вертикальными
ГТ1
V dl )
V У х=еот,
Я =
и горизонтальными е (dt
dY
V ^ ) I =оот,
градиентами. Для модели изменения скорости только с глубиной V (Z), обратный вертикальный градиент поля времён рефрагированной волны равен скорости на максимальной глубине проникновения соответствующего луча.
Допустим T (l) - нормальный годограф, V (Z) - нормальный разрез. За нормальный годограф примем осреднённую на всём рассматриваемом профиле функцию t (l)x=const (вертикальное сечение поля времён, аналог годографа ОГТ). Аппроксимируем годограф T (l) отрезками прямых линий. В данном случае для T(lk ) мы имеем:
mi, \ 2k 2k 2k, , l, „ ч
т{1к) = lk + —2cosilk +... + cosik_l + . (6)
V V2 Vk-1 Vk Найдём функцию V (Z) известным способом (4) в виде ступенчатой функции:
V (z >
V при 0 < zkj
V2 при k < zk1 + k2 . (7)
к -i к
Vk -i пРи Z kk <z <Z kk 1 1
Для получения функции V (X, Z) необходимо решить уравнение (3), которое приближённо представим в виде:
п
т(х, I) = £& -Мк (8)
к=1
Здесь «к - аномалия обратной скорости на интервале луча в к- м слое. Аномалия обратной скорости
«к - 1 1
Г(х.г) Г(г)
принимается постоянной величиной в пределах каждого прямого отрезка луча, «к (х)
относится к середине Х соответствующего отрезка луча, индекс к обозначает слой, к которому относится аномалия.
Функция т(х, I ) = £ (х, I) — Т (х, I) описывает аномалии поля времён £ (х, I)
относительно нормального поля времен Т (х, I), соответствующего нормальному разрезу
V Для удобства упростим обозначения, придадим т индекс, соответствующий номеру изолиний поля времен:
т(х ¡1)=т (4 т(х I)=т2 (x),..., т(x, 1к)=тк(х) (9)
Если обозначить через:
^к—1,к = к^ЕРк—1,к
горизонтальную проекцию луча в слое (к-1), преломившегося на кровле слоя Ц, уравнение (9) для произвольной к-той изолинии поля времен можно записать в следующем виде:
008 11
#1
( 1к 1 т
х- — — Ь
2 2
008 1
2к
'к-1
( I
к
х - "Т" + Т Ь2,к + Ь1к
V 22 )
1
1к
Л
г ( 1к 1 г
+#, х+у - ^ Ь1к
+
+ #2
( I
х + к О Ь2,к Ь1,к
V 22 )
1
V
008 1
к-1,к
( I
к-2
( I,
+ ... +
к-I
#к-1 х +1 Ьк-1 | + Х Ьк +#к-1 х + ^ -1 Ьк-1 Ь
22
к-1
+ #к(х-2)ХЬ
(10)
Зная нормальный разрез и аномальное поле времен г(х, I), по уравнению (10) последовательно от #1 (х) до #к (х) найдём значения #к (х) для к - слойного разреза. По изолинии Г)(х) устанавливается аномалия обратной скорости в первом верхнем слое:
# (х) =
(х )
Т\\х
11
(11)
По линии г2(х), зная аномалии #(х) для первого слоя, определяются значения #2 (х) для второго слоя:
Ь(х )-
#2 (х ) =
к 0081
#1
1к К
х - — + - Ь 2 2
1,2
+ #1
х +
— -11, 22
1,к
(12 - 2Ьи2 )
По линии тк(х), зная функции для к-1 слоёв, имеем
#2 (х)= К (х
к1
008/,,
( к
1
(
к
0081
2к
#2
( 1к х-—+ -Ь 2 2
х --+2 !+#,
1
Ь 1 г
х+ —--Ь 22
1,к
к-2
( I 1
г=1
к-1 V
Ь2к + Ь1к |+#2(х)+1 Ьк-1,кЬк +#к-1 х+-1 Ьк-1,к-ЕЬ
22
1=1 )_
} -
к-1
1 к - 2Е Ьк-1,к
г=1
Зная нормальный разрез V (2) и функции обратной скорости #к(х), распределение скорости в каждом слое к определяется как:
(л у1
—(х, 2к ) = ^ - + #к (х
(12)
Во-первых, главное условие при линеаризованной постановке # < 1/У(2). Для консолидированной части земной коры, где преобладают скорости 6.0 - 6.5 км/сек, ? 1
соотношения между # и —-¡-\ принимают примерно следующие значения: для разницы
—(г)
АУ = 0.1 км/сек между значениями истинной V (х,2) и нормальной V (2) скоростей # примерно в 60-70 раз меньше 1АУ(2); для 0.2 км/сек # - в 30 раз меньше 1АУ(2), для АУ = 0.3 км/сек - в 20 раз. Таким образом, нормальный годограф Т(х, 1) или нормальное поле времен Т (х, 1) надо выбирать по полю времен таким образом, чтобы отклонения т(х, 1) были минимальными.
Для получения возможно более полной информации по нижней части коры, при решении обратной задачи, необходимо использовать продолжение, во-вторых, вступлениях «прямой» волны после выхода в первые вступления волны, преломлённой
V
)
к
2
1=1
1=1
1=1
V
V
)
)
на границе Мохоровичича (М). При решении задачи отыскиваются глубины точек Нк, которым соответствуют истинные значения скорости V численно равные кажущейся скорости V на участке годографа, отстоящего от пункта взрыва на расстоянии Хк. Абсциссы точек со значениями Vист находятся путём деления X пополам. Помимо первых вступлений учитывались петли на годографах, связанные с границами раздела. Как мы видим, при обработке профильных наблюдений анализ волнового поля, выделение опорных, доминирующих волн является важным этапом.
В данной работе внимание будет уделено способам обработки рефрагированных волн и закритических отражений. Интерпретация систем годографов рефрагированных волн, расчёты сделаны по программе Мишенькиной З.Р. по модифицированному алгоритму [6]. Предпочтение этому способу решения задачи отдано в связи с тем, что он допускает возможность выделения в разрезе инверсионнных включений в отличие от других программ. Здесь мы остановимся на основных моментах методики.
Скоростные разрезы, полученные описанным способом, рассматриваются как схематическая модель, требующая уточнения. Это связано с выбранным решением одномерной обратной задачи, исключающим возможность выделения слоёв разреза с отрицательным вертикальным градиентом скорости. Все указанные выше свойства лучей и годографов рефрагированных волн в случае линейного изменения скорости с глубиной делают эту скоростную модель удобной для расчётов и общей интерпретации волновых полей ГСЗ. Большое значение имеют и модели с уменьшающимся вертикальным градиентом скорости. Для такой монотонно возрастающей функции V(z) годограф рефрагированной волны однозначен и непрерывен относительно координаты -х. Такие модели мы будем называть нормальными.
Исследования известных сейсмологов с многолетним опытом работы в области сейсмологии и сейсморазведки, докторов наук из Центральной Азии, в частности, в серии опубликованных работах докторов геолого - минералогических наук Ф.Х. Зуннунова в 1976 году [4], 1985 году [5] , Б.Б. Тальвирского в 1982 году [12], доктора технических наук В.И. Шацилова в 1994, 1995,1996 годы [13,14,15], доктора физико-математических наук А.М. Сайипбековой в 1994,1995, 1996, 1998, 2003 годы [13,14,15, 10,11] и в разных периодах на материалах ГСЗ и региональной сети в Средней Азии показали основную роль рефрагированных волн в формировании поля упругих колебаний при больших взрывах.
В наших последних работах при решении прямой задачи был применён лучевой метод определения годографов, а что касается выбора модели, использовалась градиентная модель среды, разработанная в работах профессоров ОшГУ и ОшТУ [10]. Наш многолетний опыт показывает, модели не имеют принципиальной разницы в условиях получения решения, так как слоистую среду можно аппроксимировать непрерывной и, наоборот, градиентную модель представить как предельный случай тонкослоистой. Многие сейсмологи из Москвы на материалах ГСЗ Средней Азии показали основную роль рефрагированных волн в формировании поля упругих колебаний при взрывах [7, 8, 9].
а
6
Т
к
О
1 /
А- о
2 <*,*
к
Рисунок 1. Схема показывает распространение и определение годографа рефрагированных волн. 1- точки разреза, для которых определена скорость, 2 - значение скорости, по которым показаны изолинии скорости. В итоге мы имеем (а) годографы рефрагированных волн и (б) скоростные разрезы (скоростные разрезы построены путём сложения частей годографов).
После формирования исходной скоростной модели коры и верхней мантии проводится итеративное решение прямой двумерной задачи. При решении прямой задачи входными данными являются скелет полного скоростного разреза, построенного в рамках модели с произвольно заданным распределением значений скорости, и система встречных и нагоняющих годографов.
В методическом отношении успешное решение прямой задачи обеспечено возможностью выполнения следующих условий:
1) произвольное перераспределение значений скорости в плоскости разреза;
2) изменение положения сейсмических границ и линии равных скоростей по наклону и глубине;
3) варьирование мощности переходного слоя (сейсмической границы II рода);
4) визуальная оценка траектории распространения луча.
При построении двумерной скоростной модели земной коры и верхней мантии сначала необходимо определить наиболее устойчивую осреднённую модель коры выделить главные слои и крупные блоки, вычислить их скоростные параметры и только по ним обосновывать серьезные геологические выводы.
Двумерная модель среды в виде функции V (х,2) должна удовлетворять следующим требованиям:
а) Распределение скоростей на плоскости Z, X должно соответствовать количественно кинематическим и качественно динамическим характеристикам основных компонентов наблюденного поля.
б) Решение прямой задачи для предполагаемого скоростного разреза должно давать сходимость в допускаемых пределах со всей системой наблюдённых годографов основных волн, а не только с одним каким - либо осреднённым годографом, как это обычно делается при проверке одномерной модели. Всем известно, что реальная геологическая среда сложнее любой модели.
Двумерная скоростная модель литосферы по геотраверсу Манас — Суфикурган
Из всех проведённых специальных сейсмологических работ за инструментальный период наиболее важным являются результаты Международного профиля Тянь-Шань-Гималаи. Северный конец Международного профиля Токтогул - Сринагар пересекает Фергано-Кокшаальскую герцинскую складчатую область в пределах Кыргызстана. Этот профиль ГСЗ представляет особый интерес для теоретико-методических разработок в области оценки сейсмической опасности, так как является единственным пересечением такого рода центральной части зоны действия горизонтальных сил Индийской плиты на Евроазиатскую зону сочленения Памира и Тянь-Шаня, по-видимому, являющихся основной причиной высокой сейсмичности в Центральной Азии.
Профиль известен под разными названиями (Тянь-Шань-Памир-Гималаи, Токтогул - Сринагар, Токтогул - Нанга - Парбат) [3, 4, 5]. Он отработан в рамках работ Международного геодинамического проекта Памир-Гималаи. В этом проекте были задействованы многие геофизические организации мира. Кусочно-непрерывным продольным профилированием с использованием одиночных взрывов был отработан профиль Сринагар (Индия) на юге - Зоркуль - Каракуль - Ош - Токтогул (Кыргызская Республика) на севере и были пересечены высокие горы Гималаев, Каракорума, Гиндукуша, Памира и Южного Тянь-Шаня [6].
Территориально расклад исполнителей был таков: на юге работали специалисты Индии, на территории Пакистана - специалисты Триестского университета (Италия) во главе Бг. Антонио Марусси, в Афганистане был пропуск, а в Таджикистане, Кыргызстане, Узбекистане - работали в основном специалисты Узбекистана под руководством профессора, доктором геолого-минералогических наук Ф. Х. Зуннуновым [3, 4, 5].
Дополнительные исследования ГСЗ на участке Мургаб - оз. Зоркуль в 1978 г. проводили: Институт геологии и геофизики АН УзССР, Институт геологических наук АН КазССР, Южно-Таджикская геофизическая экспедиция Управления геологии Тадж ССР.
I!
5Щ]
б сш
Рисунок 2. Расположение двух профилей I - Венцюань-Терексай (западная часть), II - Манас-Сопукоргон (по направлению международного профиля Тянь-Шань-Памир -Гималаи). 1 - профили, 2-сейсмостанции, 3 - некоторые населенные пункты, 4-граница стран, 5 - значение Vp, 6 - номера отдельных сейсмогенерирующих зон, с севера на восток в частности, (4) - Южно-Таласская,(12)-Таласо-Ферганская, (13)-Сандалашская, (14)-Атойнокская, (15) Карасуйско - Арсланбобская, (16) -Кугарто- Талдыкская, (17) -Южно - Ферганская, (18)- Хайдарканская, (19) - Туркестано -Алайская, (20) -Гиссаро- Кокшаалская, 7 - граница сейсмогенных зон.
Зачастую недостаточная детальность в отдельных районах исследуемой территории при построении трёхмерной скоростной модели земной коры и верхней мантии исключает возможность использования её результатов в комплексе с ГСЗ на сопоставимом уровне. Поэтому разработан геотраверс, в отдельных участках стыкующихся с профилями ГСЗ, и двух связующих вдоль простирания основных структур орогена, а методика профильной сейсмической томографии изложена в работах [9, 10, 11].
В обработку включены записи объёмных волн местных землетрясений с & > 8.1, попадающие в полосу 10 км от осевой линии за период наблюдения 1970 - 1994 гг. Материал за 25-летний период позволил отобрать записи только с чёткими вступлениями фаз. Используя встречные и нагоняющие годографы рефрагированных волн по
региональной сети, совместно с опубликованными годографами ГСЗ рассчитан скоростной разрез по профилю.
В работе использованы записи сейсмических станций Института сейсмологии НАН КР Манас - Аркит - Токтогул - Аккия - Арсланбоб - Суфи-Курган (рисунок 1), положение профиля - II. Профиль пересекает Таласскую антиклиналь, Баубашитинский горный массив, Караунгурскую и Кугартскую синклинали, Узгенскую мульду, Алдиярский массив и северовосточные отроги Алайского хребта.
При построении скоростного разреза вдоль профиля на участке Токтогул - Узген-Ош дополнительно использованы аналоговые записи промышленных взрывов отряда Токтогул. Кроме того, при целевом наборе материалов для системы годографов наряду с данными стандартной обработки, нам пришлось работать непосредственно с сейсмограммами. Последнее относится и к случаю работы в области смены волн в первых вступлениях Р- и Рп и S- на Sn- (эпицентральные расстояния 180-200 км в зависимости от глубины залегания поверхности М), где идентификация классов волн по индивидуальным записям затруднительна по объективным причинам даже для опытных сейсмологов и материалы обработки необходимо критически анализировать с учётом конкретных геологических условий вдоль профиля.
Преломлённые от верхов мантии Рп-волны обладают длинными периодами, за ними следуют более интенсивные и чёткие волны с меньшими периодами, называемые Р и распространяющиеся в коре. Скорости их для сглаженных годографов выдерживаются при значительных перекрытиях в довольно узких пределах 7.8-8.2 км/с; нагоняющие годографы приблизительно параллельны.
Группы первых волн Рg в первых вступлениях очень сильно зависят от мощности консолидированной коры и может 10 до 250-300 км в горных районах. Граничные скорости обычно колеблются в пределах 5.5 - 6.5 км/сек. С выходом в первые вступления волн от границы М наблюдается резкий излом годографа со сменой кажущихся скоростей с 6.5 на 7.8 - 8.2 км/сек. В сейсмологии годографы на близких расстояниях от сейсмических источников к коре обычно выглядят как две пересекающиеся прямые. Одна из них, для меньших расстояний, отвечает скорости 6 км/сек, а другая, для больших расстояний - скорости 8 км/с. Первую относят к прямым волнам Р, распространяющимся в коре и обозначают Р и Pg, что символизирует гранитный слой.
При анализе волновой картины при зондировании литосферы Тянь-Шаня с применением записей землетрясений наблюдается весьма сложная структура волнового поля. Ранее нами отмечалось [11], при энергетическом классе 8 и 10 для глубины очага 5 -10 чёткие записи продольных и поперечных волн наблюдаются до расстояний 250 км. Это на региональном уровне, по региональной сейсмологической сети нормально, чем выше класс тем больше глибина зондирования ( глубина проникновения литосферу), это очевидно.
Фактически при сопоставлении записей взрывов и землетрясений при одинаковых эпицентральных расстояниях, энергия коровых землетрясений 8 класса эквивалентна энергии взрыва в озере 5 тонн тротила при глубине погружения 50 - 65 м. А также мы проанализировали аналоговые записи отдельных взрывов и землетрясений вдоль профиля, определили видимые частоты, количество фаз, интенсивность и кажущиеся скорости сопоставимы с параметрами продольных волн, зарегистрированных от специальных взрывов. Следует отметить, что поперечные волны на Z составляющей приборов до расстояний 80 км сильно забивают область регистрации отраженных волн.
Нам удалось построить двумерную скоростную модель путём комбинирования сейсморазведочных и региональных сейсмологических данных. На рисунке 2 показана только часть профиля II (где использованы данные 7 региональных сейсмостанций
Института сейсмологии НАН КР). В Кыргызстане профиль известен названием -Тянь-Шань-Памир-Гималаи [3, 4, 5]. Выше было отмечено, что профиль отработан в рамках Международного исследовательского геодинамического проекта Памир-Гималаи, поэтому здесь мы будем только «соисполнителями» по нашей территории.
г ;
Ч 9
Рисунок 3. Геологические и скоростные разрезы по профилю Манас-Сопукоргон:
1-3 Обобщенные стратифицированные комплексы;
1- мезокайназойский платформенный и орогенный чехол,
2- террогенно-вулканогенные образования палеозоя,
3- допалеозойские образования.
4-5 Интрузивные комплексы; 4- кислые; 5- основные.
6- Тектонические нарушения.
7- Изолинии скорости в км /с, истинные скорости.
8- Изолинии Vp/Vs (в числителе) и значения коэффициента Пуассона (в знаменателе).
9- Обобщенное положение подошвы коры.
10- Графики в точках пересечения с профилями ГСЗ: пунктирно-профильная томографии.
Наличие в архивах двух сейсмологических Институтов Кыргызской и Казахской Республик бюллетений землетрясений и сейсмограммного материала за 35 летний период наблюдений, позволили построить двумерные скоростные модели по отдельным геотраверсам, которые будут изложены после тщательной проверки с точки зрения статистики совместно с ведущими геологами Института сейсмологии НАН КР в следующих номерах журнала «Вестник Института сейсмологии НАН КР».
В целом, при зондировании литосферы Тянь-Шаня с применением записей землетрясений наблюдается весьма сложная структура волнового поля.
Получена высокая степень соответствия экспериментальных и расчётных годографов для верхней и нижней части коры. Из всех пяти профилей материал данного профиля по количеству и по качеству превосходит, одним словом, по части обеспеченности кондиционным материалам, особый. Подчеркнута четкая смена с образованием петли за счет резкой границы фундамент - осадки. Получены петли, подтверждающие смену волн на внутри коровых границах.
В итоге построена двумерная скоростная модель земной коры и верхней мантии по направлению Международного профиля ГСЗ Тянь-Шан -Памир-Гималаи. Разрез представлен в масштабе 1:500 000 (рисунок 3). Здесь мы представили чёрно -белый вариант, со временем намерены изменить масштаб разреза, для отдельных геологических структур будем еще уточнять с помощью наблюдённого поля цифрой региональной сейсмологической сети Казахстана, Кыргызстана и Таджикстана.
На данный период мы можем отметить, что характерной особенностью разреза является мозаичная форма изолиний скорости в верхней части земной коры, с их сгущением в районе Южно-Ферганского разлома. Обычно зоны аномального сгущения изолиний соответствуют сейсмическим границам, ответственным за наличие петель на годографах (например, район Токтогула). В районе Северного Памира в интервале глубин 80-160 км выявлена низкоскоростная зона.
ЛИТЕРАТУРА
1. Алексеев А.С., Лаврентьев М.М., Нерсесов И. Л., Романов В.Г. Результаты изучения горизонтальных неоднородностей мантии по профилю Памиро - Байкал // Математические проблемы геофизики. Новосибирск: Наука, 1971.- вып. 2, С.143-165.
2. Бессонова Э.Р., Джонсон Л.П., Ситникова Г.А., Фишман В.М. Решение обратной задачи сейсмологии методом 1ау (р). // Вычислительная сейсмология. Москва: Наука,1974. вып.7. С.82-99.
3. Белоусов В.В., Беляевский Н.А., Борисов А.А. Строение литосферы по профилю глубинного сейсмического зондирования Тянь-Шаня -Памир- Каракорум- Гималаи. Советская геология, 1979, № 1, с. 11- 28.
4. Зуннунов Ф.Х. Земная кора Южного и Восточного Узбекистана -// Геофизический сборник АН УкрССР, вып. 72, Киев:Наукова Думка, 1976, С 65-72.
5. Зуннунов Ф.Х. Литосфера Средней Азии по сейсмическим данным. Ташкент: Фан, 1985.-108 C.
6. Мишенькина З.Р., Шелудько И.Ф., Крылов С.В. Использование линеаризованной
постановки обратной кинематической задачи для двумерных полей времени t (x, l) рефрагированных волн // Численные методы в сейсмических исследованиях, Новосибирск: Наука, 1983. С.140-152.
7. Никитин А.А. Теоретические основы обработки геофизической информации. Москва: Недра, 1986. - 342 С.
8. Павленкова Н.И. Волновые поля и модель земной коры. Киев: Наукова думка, 1973. - 19 C.
9. Павленкова Н.И. Метод глубинного сейсмического зондирования основные этапы развития достижения и проблемы. // Физика Земли. 1999. Июль-Август. №7-8. С.3-29.
10. Сайипбекова А.М., Сатыбаев А. Дж. Методы решения обратной кинематической задачи сейсмики. //Сбор. науч. трудов Кыргызско-Узбекского университета, Ош: 1998. С.45-49.
11. Сайипбекова А.М., Абдрахматов К.Е., Степаненко Н.П., Саипбекова С.Э. Скоростная модель литосферы Тянь-Шаня по данным геотраверса Бачу - Балхаш. Известия НАН КР, 2023. №8, С.138-144.
12. Таль-Вирский Б.Б. Геофизические поля и тектоника Средней Азии. М.: Недра, 1982. -271 С.
13. Шацилов В.И., Сайипбекова А.М. Профильная томография литосферы Тянь-Шаня по материалам региональной сейсмологии. Докл. НАН Республики Казахстан, 1994. №6, С. 47-52.
14. Шацилов В.И., Сайипбекова А.М., Грибанов Ю.Е. Изучение глубинного строения Тянь-Шаня по материалам региональной сейсмологии. // Journal Inland Earthquake. Urumqi: 1995. Vol.9, №4, С.374-381.
15. Шацилов В.И., Сайипбекова А.М., Кушчубекова А.А., Искандерова А.М. и др. Комплексная обработка данных региональной сейсмологии и ГСЗ с целью получения детальных скоростных разрезов. // «Проблемы механики и прикладной математики»: Матер. межд. научно- методической конференции посвященной памяти профессора Ф.И. Франкля. - Бишкек: 1996. Том 1, Механика. С.123-125.
Рецензент: д. ф.-м.н., Токтосопиев А.М.
