25018 (686696), страница 2

Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 2)

Рис. 2. Иллюстрация к понятию времени запаздывания.

В литературе описано много способов интерпретации, использующих время запаздывания. Например, такие способы предложены Гарднером, Бартелмсом, Таррантом, Виробеком, Барри. Рассмотрим только три последних. Методы, описанные Виробеком и Таррантом, подходят для одиночных годографов, метод Барри дает наилучшие результаты в случае встречных годографов.

б) Метод Барри. Схема, описанная Барри, подобно многим, основанным на временах запаздывания, требует разложения полного времени запаздывания на составляющие члены _S и _g. На рис. 3 изображен приемник R, который регистрирует колебания от источников A и В. Луч BN отражается под критическим углом; следовательно, Q — первый приемник, который зарегистрирует головную волну, порожденную источником В. Пусть _ам — время запаздывания в пункте взрыва A, _nq и _PR — времена запаздывания в пунктах приема Q и R, a _aq и _ar — полные времена запаздывания для траекторий AMNQ и AMPR..

Тогда

Время запаздывания в пункте взрыва В _BN, в случае если наклон границы мал, приблизительно равно времени запаздывания в пункте приема Q _NQ. Следовательно,

.

Времена запаздывания в пунктах приема теперь можно записать в следующем виде:

(3)

Таким образом, время запаздывания в пункте приема R можно получить, если имеются данные для двух пунктов взрыва

Рис. 3. Определение времен запаздывания в пункте взрыва и пункте приема.

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

(4)

(5)

Принимаем, что время запаздывания _bn равно половине t₀ в точке В; это позволяет рассчитать приближенное значение BQ и определить таким образом времена запаздывания для всех сейсмоприемников вправо от Q, которые зарегистрировали колебания, возбужденные в точках А и В.

Этот способ интерпретации включает следующие шаги, которые можно проследить по рис.4:

а) построение годографов по исправленным временам;

б) расчет и построение графиков полных времен запаздывания для всех положений приемников;

в) расчет «сейсмического сноса для сейсмоприемников» (РР' на рис. 4) с помощью соотношения РР' ≈ V_2PRtg², после чего кривые времен запаздывания в п. (б) сдвигаются по направлению к пункту взрыва на эти величины;

г) смещенные на этапе (в) кривые для встречных профилей должны быть параллельны; любое расхождение обусловлено

Рис. 4. Интерпретация встречных профилей по методу времен запаздывания.

неточным выбором значения V₂; следовательно, значение V₂ необходимо исправить и повторить этапы (б), (в), пока кривые не станут параллельны (на практике уточнение V₂ обычно производится только один раз);

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

в) Метод Тарранта. В этом методе времена запаздывания используются для определения положения точки Q (рис. 5, а), в которой энергия, регистрируемая в пункте R, отходит от границы. Обозначив _g время запаздывания, связанное с траекторией QR, запишем

,

Откуда

. (6)

Мы получили уравнение эллипса в полярной системе координат. Эллипс —это геометрическое место таких точек Q (рис. 5, б),

Рис. 5. Интерпретация по методу Тарранта а — связь между точкой приема К и точкой Q отхода от границы; б —схема, поясняющая, что геометрическим местом точек Q является эллипс, в одном из фокусов которого располагается точка R; в — геометрия эллипса, проходящего через точку Q.

для которых отношение QR/QM остается постоянным (равным эксцентриситету , который для эллипса меньше 1, т. е.

/(h + cos) = ,

а следовательно,

= h/(1 - cos). (7)

Большая ось эллипса 2а = ₌₀ + ₌ = 2h/(1 —²). Малую полуось b можно найти, записав y=sin и определив у_max,; это дает b = h(1 — ²)^-1/2. Расстояние от фокуса R до центра эллипса О равно

|₌₀ — а = h /(1 — ) —h/(1 — ²) = a.

Если принять = sin и h = V_2g, выражение (7) перейдет в (6).

Для горизонтальной преломляющей границы получается эллипс (рис. 5, в) с а = V_2g tgsec, b = V_2gtg и OR = V_2gtg². Подобным же образом RQ = b/cos = а и OQR = arctg(ОR/b) = , OQ = ORctg = V_2gtg.

В окрестности Q эллипс можно аппроксимировать окружностью того же радиуса кривизны. Если записать уравнение эллипса в декартовой системе координат

(x/a)²+(y/b)² = 1,

то радиус кривизны r можно выразить как

r = (1+y’²)^3/2/y’’,

где у' = — (b/а)²(х + у) и у" = —(b/а)²(у — ху')/у²; в точке Q у' = 0 и у" = b/а². Следовательно,

r = a²lb = V_2g/cos³ = V_2g tgsec²

и центр кривизны С лежит в точке (0, r — b), т. е. (0, V_2gtg³). Следовательно, CRO = arctg(CO/RO) = , а значит, CRQ — прямой угол.

Чтобы применить описанный метод, мы должны определить скорости V₁ и V₂ и время запаздывания в пункте взрыва _S, а затем рассчитать _g по формуле

_g = t_R — x/V₂ — _S.

После этого можно вычислить OR, OQ и затем найти положение С, проведя перпендикуляр RC к RQ. Из С проводим дугу окружности, соответствующую преломляющей границе в окрестности точки Q. Если наклон границы отличен от нуля, точкой выхода станет Q' и длина дуги QQ' увеличивается при росте угла наклона границы. Но даже для углов падения умеренной величины дуга эллипса QQ' будет близка к дуге окружности, проходящей через Q, и, таким образом, огибающая дуг окружностей достаточно точно отобразит преломляющую границу.

Метод Тарранта удобен, когда наклон границ умеренный или даже большой, а преломляющая граница криволинейна или имеет неправильную форму. Принципиальным ограничением является точность определения V₂.

г) Метод Виробека. Для иллюстрации метода Виробека обратимся к верхней части рис. 6, где показан ряд одиночных годографов. Последовательные шаги интерпретации таковы:

а) строят исправленные годографы и измеряют времена t₀, отсекаемые годографами на оси t;

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

в) строят график времени t₀, деленного на 2, и сопоставляют его со сводной кривой запаздывания; расхождение между двумя

Рис. 6. Интерпретация профилей, отработанных в одном направлении, по методу Виробека.

кривыми указывает на то, что значение V₂ выбрано неверно (см. ниже), поэтому значение, использованное на шаге (б), следует уточнять, пока две кривые не станут «параллельны», после чего кривую t₀/2 дополняют путем интерполяции и экстраполяции так, чтобы она покрывала тот же диапазон, что и сводная кривая временных задержек;

г) кривую t₀/2 преобразуют в кривую глубин, используя выражение

Успех применения метода Виробека зависит от того, является ли кривая приблизительно параллельной кривой t₀/2 Чтобы применять данный метод, не требуется встречных профилей, так как t₀ не зависит от направления, в котором развернута приемная коса.

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

Метод волновых фронтов (полей времен)

а) Метод Торнберга. Целый ряд способов интерпретации данных МПВ, как правило графических, основан на реконструкции волновых фронтов.

Рис. 7 иллюстрирует принцип метода построения волновых фронтов. Фронт преломленной волны, достигающий точки А в момент времени t = 1,600 с, подходит к В, С, ... в моменты 1,600 + t_B, 1,600 + t_c, .... Построив дуги окружностей с центрами В, С, ... и радиусами V₁ t_B, V₁ t_C,,… , мы можем восстановить волновой фронт для t = 1,600 с (AZ) с требуемой точностью. Подобным же образом можно построить волновые фронты для любого момента времени (см., например, приведенный на том же рисунке волновой фронт для t = 1,400 с) Показаны также волновые фронты прямой волны от источника S, являющиеся окружностями.

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

Две системы волновых фронтов, соответствующие прямой волне и волне, преломленной на первой границе, пересекаются

Характеристики

Список файлов курсовой работы