Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Богданов, Алексей Иванович. Интерпретация сейсмических годографов

.pdf
Скачиваний:
6
Добавлен:
29.10.2023
Размер:
10.99 Mб
Скачать

определим положение фронта отраженной волны в момент вре­ мени ti — kt. Аналогичным путем найдем положение фронта отраженной волны в любые моменты времени.

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

Размеры отражающей площадки определяются точками ее пересечения с прямыми, проведенными перпендикулярно фронтам

отраженных

воли из крайних то­

 

 

чек интервала профиля, в преде­

 

 

лах которого был наблюден про­

 

 

дольный

годограф

отраженных

 

 

волн.

 

 

 

 

 

 

 

 

 

Недостаток метода полей вре­

 

 

мен заключает&я в громоздких

 

 

построениях

фронтов

главным

 

 

образом

отраженных волн.

Этим

 

 

объясняется весьма

ограниченное

 

 

применение

 

метода

полей

вре­

 

 

мен в практике сейсморазведочных

 

 

работ, несмотря на его универ­

 

 

сальность.

 

 

 

 

 

 

 

Метод мнимой линии

 

 

сводится к следующему. Имея

 

 

линейный

продольный

годограф

 

 

отраженных

волн

и

зная

ско­

Рис. 16. Построение отражающей

рость

гу в

среде,

покрывающей

границы при помощи линейки (по

отражающую

границу

раздела,

В. Д. Завьялову и Ю. В. Тимо­

строим мнимую линию как огибаю­

шину).

 

1 — отражающая граница;

2 — мни

щую

к

семейству

окружностей.

мая линия h.

 

Центры

семейства

расположены

 

 

 

в произвольных точках Р дневной

произведениям из

скорости

поверхности.

Радиусы

РО"

равны

v1 на время

tp прихода отраженной волны в точку Р. Точки, на­

ходящиеся на отражающей границе

R, лежат на перпендикуля­

рах, восстановленных к серединам прямых, соединяющих любую точку мнимой линии с пунктом взрыва (рис. 8). На основании этого положение отражающей границы R может быть найдено как положение огибающей к семейству перпендикуляров, вос­ становленных к серединам прямых 00". Для отыскания отражаю­ щей границы R В. Д. Завьяловым и Ю. В. Тимошиным предло­ жено пользоваться двумя перпендикулярно скрепленными между собой линейками, одна из которых (измерительная) имеет две равномерные шкалы с общим нулем в точке пересечения линеек (рис. 16). Зная положение пункта взрыва О и мнимой линии L, устанавливаем измерительную линейку так, чтобы на точке

39

взрыва

и на мнимой линии отмечались одинаковые значения но

шкале.

Вследствие этого перпендикулярная линейка будет уста­

новлена

по середине отрезка, соединяющего

точку взрыва О

с произвольной точкой О", находящейся на

мнимой линии.

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

Всинклинальных зонах огибающая будет проходить сверху,

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

§6. ОПРЕДЕЛЕНИЕ ЭЛЕМЕНТОВ ЗАЛЕГАНИЯ ПЛОСКОЙ ОТРАЖАЮЩЕЙ ПОВЕРХНОСТИ И СКОРОСТИ РАСПРОСТРАНЕНИЯ ВОЛН В ПОКРЫВАЮЩЕЙ ЕЕ ОДНОРОДНОЙ СРЕДЕ ПО ГОДОГРАФАМ

ОТРАЖЕННЫХ ВОЛН

Но поверхностному годографу. Если отражающая поверх­ ность плоская и известна скорость vY распространения упругих волн в покрывающей ее однородной среде, то элементы залегания ее в пространстве могут быть определены по заданным или опре­ деленным параметрам х, у, t, тх и ту любой одной точки Р по­ верхностного годографа отраженных волн на основании следую­ щих формул: (1.30), (1.31), (1.32), (1.40), (1.41) и (1.42).

Эти формулы имеют наиболее простой вид, если точка Р взята в начале координат (в пункте взрыва), т. е.

*о = — *4% (Тх)0,

lg а =

,

 

 

Пх/о

Уо =

(Ту)0, sin ф =

т0,

ZQ = V1t /1 -

*,(?)„ h = ^.

По тем же параметрам и формулам (I. 34),

(I. 35), (I. 36) можно

рассчитать координаты х',

у' и z

точки, лежащей на отражающей

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

(/гр)

или к поверхности наблюдения (z').

1

Изложенный прием является развитием широко известного способа

засечек; он был опубликован в 1955 г., но не получил еще широкого произ-

. водственного опробования.

40

Глубину z' можно найти из уравнения (I. 8), решая его относи­ тельно z':

Zo

Полагая в этом уравнении х = х, а у' = у, заменяя ха, yt) и z0 через их значения из (I. 30), (I. 31), (I. 32), выполнив преобра­ зования, найдем

2 vtt У1 —г>!2 Т2

7

Глубину hv можно определить по координатам х , у' и z осно­ вания перпендикуляра, опущенного из точки Р на отражающую

плоскость. По известным х , у' и z находим hp

из уравнения

hP = V (х' — .г)2 + (у' — у)2 4- (z')2 ,

(1.46)

где х и у — координаты точки Р.

Координаты х', у' и z можно получить из совместного реше­ ния уравнения (1.8) и уравнения прямой, перпендикулярной отражающей плоскости и проходящей через точку Р.

Уравнение прямой в пространстве, проходящей через точку Р и нормальной к отражающей плоскости, можно записать через координаты х и у точки Р и направляющие косинусы этой прямой в пространстве. Направляющие же косинусы равны направляю­ щим косинусам перпендикуляра, опущенного из пункта взрыва на отражающую плоскость (см. § 4). Следовательно, уравнение перпендикуляра, опущенного из точки Р на отражающую пло­ скость, можно представить в виде

 

 

х'~х

= У'~У = z'

(1.47)

 

 

 

хо

 

Уо

 

где х ,

у' и

z' — текущие координаты точек, лежащих на перпен­

дикуляре, опущенном из точки Р

на отражающую

плоскость;

х0, у0

и z0

— координаты мнимого пункта взрыва.

 

Из

уравнения следует,

что

 

 

 

 

 

X

'

 

z'

 

 

 

 

 

= х + х0

— ,

 

 

 

 

 

,

I

z'

 

 

 

 

У

 

= У + Уо~..

 

 

 

 

 

 

 

Z0

 

Подставляя х' и у’ в уравнение (I. 8) и решая полученное урав­

нение относительно z ,

получим

 

 

 

 

z

_

 

2/г2 — (хох + УоУ) „

 

 

 

=

------- —------- z0.

 

 

 

 

 

 

4/г2

 

 

41

После подстановки

г

в

уравнение для

х и у' имеем

1

.

,

2h2 — (xox4-?/„?/)

 

— •'

+ J о

/(Л2

 

У

—У+Уо------/Л2

Подставив х', у' и z в уравнение (I. 46) и произведя преобра­ зования, находим

,2h2 — (x0x + y0y)

 

Пр -

2h

После замены х0, у0

и h через параметры х, у, t, тх и ху точки

Р и преобразований

имеем

 

 

2 |/v12t2-l[-ri—2v12t (х т-с-)-у ху)

Для плоской отражающей поверхности х0, у0 и z0 — величины постоянные. Они полностью определяют положение отражающей плоскости в пространстве и могут быть получены независимо по параметрам любых точек поверхностного годографа отраженных волн. На основании этого представляется возможность по задан­ ным или наблюденным параметрам любой пары точек поверхност­ ного годографа отраженных волн определить независимо не только х0, у0 и z0, но и скорость vr. Записав полученные выражения для х0, у0 и z0 через параметры двух произвольных точек В и А по­ верхностного годографа отраженных волн, получим три уравне­

ния, по которым

можно определить скорость г^:

 

Хь — V-^tb (тж)ь = Ха — l\2ta (Тх)а,

 

уь

Vy2tb (Ау)ь = Уа — V2ta

 

Vrtb V 1

— У].2 Ть2 = Vyta

— Гх2 Та2 .

Решение этих

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

дает

v

= 1/ tb2-ta2 —

1

v tb2Tb2~ta2xn2 ■

Подставив vr в уравнения для х0, у0, z0, можно получить девять независимых выражений для их определения, наиболее удобными из которых являются следующие:

Xgtb (Xx)b Xbtg (Хх)д

lb (Тт)Ь (Хх)а

42

..

*Уа-ъ(Ху)ь

ybta(ty)a

,, r<>\

Уоtb(xv)b-ta(xy)a ’

(1 • 5<))

zo = 7T- ?<az 2 т г V (tb2 - 7a2) (Tb2 - Xa2) .

Zb

Tb ta Xa

 

Если был полностью наблюден поверхностный годограф волн, отраженных от плоской поверхности, изображенной в виде карты изохрон, то по положению центра окружности семейства изохрон можно найти координаты л0 и //0 мнимого пункта взрыва и тем самым направление восстания или падения отражающей пло­ скости. Направление восстания определяется направлением пря­ мой, проведенной из начала координат (пункта взрыва) в центр окружности семейства изохрон.

Непосредственное определение хп и у0 с карты изохрон позво­ ляет найти независимо z0 и г\ по более простым формулам, в ко­ торые не будут входить Хх и ту, определяемые недостаточно точно. Исходные формулы для z0 и vi можно получить из уравне­ ния поверхностного годографа волн, отраженных от плоской границы раздела, если его записать для двух произвольных точек В А поверхностного годографа с координатами хь, Уь, tb и ха, уа, ta. Тогда имеем

v^tb2 = (хь — .г0)2 + (уь — Уо)2 + z02,

гЛа2 = (Ха — А))2 + (Уа — Уо)2 + Z02.

Разделив первое уравнение на второе и разрешив полученное соотношение относительно z0, найдем

- _ I /fa2 [fob —А>)2 + (№ —Уо)2] —гЬ2 [жа—Жр)2 + (Уа —Уо)2] ° Г Zb2 —/а2

Если из первого уравнения

вычесть

второе и решить

полу­

ченное равенство относительно

vi, то

получим

 

1/ (^ь —А)2 —(«а—а;о)2 + (уь —Уо)2 —(ya —Уо)2

,,

V1 = Г

 

 

и

Уравнения для z0 и vi могут быть записаны в более простой форме, а именно,

, _ 1/ ZaW-/b2^a2

°V tb*-t a*

=1/ 7У-Да2

где Вь и Ra — радиусы изохрон, проходящих через точки В и А, равные

Rb = V(хь — А))2 + (уь — Уо)2 > Ra = V (Ха — Т0>)2 + (уа — у0)2.

4.3

Имеется принципиальная возможность [52 ] определить vi без знания х0 и у0 только по параметрам х, у и t трех точек С, В и А поверхностного годографа волн, отраженных от плоскости, если на параметры х и у этих точек наложить следующее ограничение. Все три точки должны лежать па одной прямой линии, как угодно ориентированной относительно осей Ох и Оу, па равном расстоя­ нии S друг от друга. В этом случае

Хс — Хь = Хь — Ха = &х,

Ус— Уь = уь — Уа= А у,

5 = /(Ах)2 + (Ау)2.

Запишем уравнение (I. 51) для точек С и В и В и А:

Ч2 *с( 2 — tb2) = хс2 — хь2 — 2х0 (хс — хь) + Ус2 — Уь2 — 2у0 (ус — уь) =

= А X (хс + хь — 2ж0) + А у (ус + уь — 2у0),

V 2 (tb2 — ta2) = Хь2 — Ха2 — 2.Г0 (хь — Ха) + уЬ2 ~ Уа2 — 2уд(уъ~ Уа) =

= А ж (хь + Ха — 2хд) + А у (уъ + уа — 2у0).

Разность этих двух уравнений будет

г-’l2 (tc2 — 2tb2 + ta2) = к х(хс — Хь + Хь — Ха) +

+ А у (Ус — уь 4- уь — Уа) = 2 А ж2 + 2\у2 = 2S2,

откуда

 

5

 

(1.52)

V1 = ---- ---------------------

рЛ

—2tb2

-\-tg2

 

 

Если в двух точках В и А поверхностного годографа время прихода отраженных волн остается постоянным, т. е. если tb = = ta, то формула для определения vi получает вид:

S

1/«с2 —Zb2

Г2

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

ние

годографа

отраженных волн

на этом профиле имеет вид:

 

 

t = -L.y (х-х0}2 + Н2,

где

Н — длина

перпендикуляра,

опущенного на этот профиль

из мнимого пункта взрыва.

44

В зависимости от того, является ли профиль продольным или непродольпым, и от ориентировки его относительно залегания отражающей плоскости II будет различным (см. § 4).

Линейные годографы волн, отраженных от плоской поверх­ ности, представляют собой гиперболы, симметричные относи­ тельно оси, параллельной оси Ot и смещенной по оси Ох на отре­ зок х0 от начала координат. Следовательно, величина х0 равна абсциссе точки минимума линейного годографа и определяется из наблюденного годографа непосредственно. Ордината же точки

минимума линейного

годографа

 

 

 

t

- 11

 

откуда

£мин ~

 

 

 

 

 

и

,

 

 

11

Н

^1^мин

ИЛИ

-мин .

Подставляя в выражение для t вместо vi его значение через Н и

£мин и решив полученное

уравнение относительно Н, найдем

 

_

(х хо)

^мин

 

гмин

Если же в выражение для t вместо Н подставить его значение через vi и fMIIH и решить полученное уравнение относительно щ, то получим

,, _

х0)

 

 

По приведенным формулам могут быть независимо определе­

ны величины vi и Н по координатам х0 и £мин

точки минимума

линейного годографа и по

координатам х и t

любой его

другой

точки.

 

 

 

 

Эту же задачу можно решить по параметрам х, t и' хх

любой

пары или по параметрам

х,

t трех точек линейного годографа,

расположенных на равных расстояниях S друг от друга,

по фор­

мулам предыдущего раздела.

Для повышения точности определения величин vx и, следова­ тельно, II по приведенным выше формулам следует брать не одну, а несколько пар или троек точек на годографе, так как снимаемые с годографа в этих точках параметры t, х и т0 мо­ гут содержать те или иные погрешности. Вследствие этого vr и Н, определенные по различным парам или тройкам точек, бу­ дут несколько различаться между собой и для нахождения наи­ более достоверных их значений следует брать средние значения из полученного ряда величин. Аналитический расчет значений г^и Н по параметрам нескольких пар или троек точек наблюденного годографа и вычисление их средних значений сложны и трудо­ емки. Для упрощения и ускорения определения vi и Н по наблю­

45

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

Рис. 17. Определение vr способом квадратных координат по линейному го­ дографу отраженных волн.

широкое применение на практике получили следующие способы определения щ и Н по наблюденным линейным годографам отра­ женных волн.

Определение и Н способом квадратичных координат. Ура­ внение (I. 21) после возведения в квадрат можно записать в виде

А(^-^о)2+4-//2- (z-53)

Введем новые переменные

t2 = Т, (х-х0)2 = X.

Тогда уравнение (I. 53) перепишем в виде

Т~±-Х + -1ГГР.

(1.54)

46

Последнее уравнение представляет собой в осях координат X и Т уравнение прямой линии с угловым коэффициентом 1/г/, кото­ рая отсекает на оси времен (при X = 0) отрезок Я2/^2.

Продолжение этой прямой отсечет па оси X (при 7 = 0) отре­ зок Н2 (рис. 17). Следовательно, перестраивая наблюденный годограф отраженных волн в оси координат Т и X, по наклону полученной прямой определяем скорость распространения упру-

Рис. 18. Определение способом постоянной разности по ли­ нейному годографу отраженных волн.

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

Для того чтобы перестроить наблюденный годограф в новые оси координат, нужно знать х0. Ошибки в ее определении приво­ дят к ошибкам в определении щ и Н. В этом заключается недо­ статок способа квадратичных координат.

Определение г>1 способом постоянной разности (по В. А. Бу-

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

47

Возьмем на наблюденном годографе отраженных волн две

точки Л и В с координатами ха,

ta, хь = ха +

(рис. 18).

Напишем для этих двух точек два

уравнения:

 

[(.Га -I- d) — То]2 + //2 = У12/ь2,

(,Г„ — ,Г0)2 4- //2 = V^t,,2

Вычитая из первого уравнения второе и преобразуя, получим

2xad 4- d2 — 2.rnd = Vj2 (tb2 — ta2).

Рис. 19. Определение

способом асимптот по линейному годографу отра­

 

женных волн.

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

tb2 — ta2 = Q.

Тогда можем написать, опуская у х значок а, что

2d

d2—2dx0

(1.55)

Q

— — х 4------------ .

Х

Vi

V]2

 

В системе координат х и

Q равенство (I. 55) будет уравнением

прямой с угловым коэффициентом, обратно

пропорциональным

г?!2.

 

 

 

На основании (1.55) наблюденный годограф отраженных волн может быть перестроен в оси координат Q и х. Для последова­ тельного ряда пар точек годографа с разностью абсцисс, равной d, необходимо вычислить разность квадратов их ординат и в вы­ бранной системе координат Q и х отложить эти разности квадратов ординат против значения наименьшей абсциссы соответствующей пары точек. Полученные точки с координатами Q и х должны лечь на одну прямую линию, тангенс угла наклона которой будет

равен 2d/y12.

способом асимптот (по С. Д. Шушакову).

Определение

Выше указывалось,

что асимптоты гиперболических линейных

48

Соседние файлы в папке книги из ГПНТБ