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

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

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

■/'

= У + 4(-/»--/)

 

 

.

(1-28)

 

2' = 4- V(lyt)2— (х — rz0)2 — (?/ — ?/0)2 X

 

 

 

(гц/)2 — х2 — у2

 

 

 

(1.29)

 

(Vil)2 — х(х — хп) — у (у — у„)

 

 

 

 

 

 

Таким

образом, координаты х , у' и z

точки

отражения 6

выражены нами через х и у точки Р и х0, у0

и

z0 точки О". Ко­

ординаты х и у точки Р нам известны, а х0,

у0

и z0

можно опре­

делить через параметры поверхностного годографа отраженных

волн в

точке Р следующим образом.

t

раздельно по х

 

у,

Дифференцируя выражения (Т. 7) для

и

получим

_

_

1 _

 

2 (х — хп)

 

_ х — хп

 

 

 

 

c)t.

 

 

 

 

 

 

 

~

 

vx*~

2(4 У(х - хи)2~+7у - у0)2

zu2 ~ ^21

'

 

 

 

дУ

 

х _ Д_

 

2 (у — у0)У —Уо

 

 

 

 

 

У

1*'«

2У1/(^-^)2 + (?/-?/0)2 +V

'

 

откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X — ,г0 — v^t хх или ./■() = ./■ — vft хх,

(1.

30)

 

 

 

 

У — Уо = Vi2t ту или ?/о

= У —

 

ty,

(Ь 31)

где

хх

и

т;/

— градиенты

времени

поверхностного годографа

в точке

Р, равные величинам, обратным

составляющим

г’х* и

Vy*

кажущихся скоростей

распространения

 

фронта отраженных

волн через эту точку в направлениях, параллельных осям ко­ ординат (Ох и Оу). После подстановки полученных значений х0 и у0 в выражение для t найдем значение z0:

z0 = vit ]/1 — г^Тх2 — г?!2 Ту2 = Vyt Kt — г’,2 г2,

(1.32)

где

 

X = |/тх2+ Ту2".

(1.33)

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

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

в

общем виде значений х0>

уо и z0 в выражения для х’, у' и

z'

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

жены только через параметры любой точки поверхностного го­ дографа отраженных волн и скорость следующим образом:

(1.34)

29

I

1

\ (гц/)2— (ж2 + у2)

У — у

2

V1t — (»J хх)х — (г\ xy)y ~

г

'

1 1'Т------

7----- v>—~7

-------7? (г>102 — (ж2

4- у2)

 

= _ I 1

- (г,т,)—

(!>1 т ,)2

 

 

 

= 4 /1 - (»,г.)2 - (-.т,)2 ,., _

 

ИЛИ

 

 

 

 

 

 

 

 

1 2xt-[(v1l)2-{-x2 — y2]xx — 2xyxy

 

 

~ 2

t - zr^yry

,1 2yt-[(i\t)2-\-y2—x2]xv — 2xyxx

У

-

~-------------------

'1

t — хт

— „т----------------

 

 

— х тх — У

''у

<'-35)

=

. (I. 36>

.. „ VbdQ

п,

<'• 39>

где I — расстояние от пункта взрыва О до точки приема Р отра­ женной волны. Первые три формулы в несколько иной форме и другим путем были впервые выведены И. С. Берзон [6].

Азимут падения отражающей поверхности в точке отражения S может быть определен по формуле

tg а = у± =

.

(1. 40)

Хо

X — v±2t хх

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

:_______ /V + Уо2

г2 — 2г/!Ч (ж тх 4- у Ху) + v^l2x2

,т ,

P'V + Уо2 + zor

/2 —2г/12г(жТх+уту)4-г>12г2

 

где

 

 

I = Ух2 4- у2 ,

т = Ухх2 4- ту2.

 

На основании приведенных выше формул для любой

точки

Р(х, у) можно рассчитать координаты точек отражения S, а также угол и азимут падения отражающей поверхности в точке отраже­ ния и сразу же использовать их для построения структурной карты. Если скорость в среде, покрывающей отражающую по­ верхность, функция I (см. гл. III, § 2), то в формулах (1.34), (1.35), (1.36) для каждого значения I должна быть взята соответ­ ствующая величина эквивалентной скорости.

Формулы (1.34), (1.35), (1.36) и расчет координат х , у' и z могут быть упрощены, если на карте изохрон выбираются не про­ извольные точки, а точки, в которых тх = 0 или ту = 0, или же ТОЧКИ, В которых Тх = Ту .

30

Очевидно, что хх = 0 в тех точках карты, в которых касатель­ ные к изохронам параллельны оси Ох. т(/ = 0 в тех точках, в ко­ торых касательные к изохронам параллельны оси Оу. И, наконец, точками карты пзохроп, в которых хх = хи по абсолютной ве­ личине, будут точки, в которых касательные к изохронам с осями координат составляют угол, равный 45°.

Рис. 12. Карта изохрон поверхностного годографа волн, отраженных от не­ плоской границы раздела, покрытой однородной средой.

О — пункт взрыва; А, В, С — точки нарты изохрон.

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

волн возрастает в сторону положительного направления

осей

Ох и Оу, и, наоборот, отрицательные значения градиента

будут

в тех точках карты изохроп, в которых время нарастает в напра­ влении, противоположном положительным направлениям осей

Ох и Оу.

31

На карте изохрон, кроме областей, в которых хх и ту одинако­ вого знака, будут области, где тх и ту разных знаков (тх поло­ жительно, Ту отрицательно и наоборот).

Карта изохрон волн, отраженных от неплоской отражающей границы раздела, изображена на рис. 12. В точках, обозначенных буквой А, тх = 0. В точках, обозначенных буквой С, тх = Ту по абсолютной величине и имеют одинаковые знаки в верхнем левом и нижнем правом углах карты (в верхнем левом углу они

отрицательны, в нижнем правом — положительны). В

этих

же

точках в верхнем правом и нижнем левом углах карты

тх и

ту

имеют разные знаки (в верхнем правом углу тх положительно,

а Ту отрицательно; в нижнем левом углу тх отрицательно, а

ту

положительно). В точках А, расположенных в верхней части рис. 12, Ту отрицательно, а в точках А, расположенных в нижней части карты, Ту положительно. В точках В, расположенных

справа от

центра карты, тх положительно, а слева тх

отрица­

тельно; во

всех точках В ту = 0.

 

В рассмотренных выше частных случаях расчетные

формулы

несколько упрощаются. Для точек,

в

которых тх

= 0,

выраже­

ния (1.34), (1.35)

 

и

(1.36)

приобретают

вид:

 

 

 

 

 

 

 

 

х

= х,

 

 

 

 

 

 

 

1

(М2_/2

_ ?#-4’[(М2 + у2-а;21Т

 

 

 

2

t — ух

 

 

 

 

 

t—ух

 

 

z

,

=

1

i 1

 

 

 

7

(УВ)2 —12

 

 

 

-7

У 1

г

0’1-Ч-2—,

 

 

 

 

 

2сг

 

 

'

t — ух

 

 

так как в этом случае ту

— т;

т может быть как положительно,

так и отрицательно.

 

ху

 

= 0, будем иметь

 

 

Для точек, в которых

 

 

 

 

1

(М2_/2

Xt~ 4-[(^02+^-У2]

 

X = X —■ -7Г- т -

-----

= ---------------------------------

,

у' --= У,

?'

1

 

 

177------ 7------\2

(y3Z) ---

 

 

-Т---- 1/

1

(i’i

Т 2

--------

 

 

 

 

 

 

2.v,

 

г

 

' 1

'

t—хх

 

 

так как в этом случае

тх = т; т по-прежнему может быть как по­

ложительно, так и отрицательно.

 

 

 

 

 

Для точек, в которых

тх

 

 

и

ту

одинаковы по величине и по

знаку, получим следующие выражения для х , у'

и z':

 

 

 

 

1

 

/

 

ч

(г^)2—г2

 

 

 

 

 

-----7=

(V1

т)

 

 

 

 

 

 

 

 

2/2

 

 

 

 

 

 

(®+у)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ч^т) ... ^2-12 2 ^2

32

'

1

Z

------ 7------v?

 

(г,^)2—^2

 

 

= —

у 1 — (Г! т)2----- ,

 

 

 

 

yiz~ -7=- (»+?/)

где

 

 

 

 

 

V 4

 

 

 

 

т

 

 

 

 

 

 

 

 

 

 

Тх — Ху — —.

 

 

 

 

 

У2

При этом т может быть и положительным и отрицательным.

Для точек, в которых

хх и

ху одинаковы по величине, но раз­

личны по знаку, формулы для х', у'

и

z' целесообразнее оставить

в форме уравнений

 

(1.34),

(1.35)

и

(1.36).

В точке минимума поверхностного годографа отраженных волн или в точке карты изохрон с минимальным временем при­ хода отраженной волны, когда хх = ху = 0:

Для точек карты изохрон, лежащих на оси Ох, формулы (I. 34), (1. 35) и (I. 36) приобретают вид (при г/=0);

,

, =

 

 

 

_ ^-4-[(^)г+^2]Тх

 

Т

2 У1

V1t— (v1Xx)x

 

t—х Хх

 

v' =____L (г, т )

— (г?! xx) x

=

t—x xx

 

У

2 ' 1 y'

 

 

z' = -у /1 -УхТ^-УхТ,)2

V{^~X\X =

 

 

z

 

 

u^i— t,x) x

 

 

 

=-4- /i-yi^-^T,)2 m2-* 2.

 

 

 

4.0'Д

 

 

 

t 11 1 «v X*

 

Для точек карты изохрон,

лежащих на оси Оу, получим (при

х = 0)

 

 

 

 

 

 

 

 

Ж' = _ J- (г; т )

^2~У2

=

t — yxy

 

 

2

1 х' v1t—(v1xv)y

 

 

1

(М2-У _

yt~~

 

У

У

2

V1t — (V1 Ху)у

 

t — yXy

 

 

z' = 4 У 1 — (г?! тж)2 — (г?х Ту)2

—4"/—4~

 

 

2 г

х 1 '

\

у/

Vlt — х \у

 

 

 

= -^Г /1 - Ух Тх)2 - Ух ту)2

 

 

 

^^1

 

 

 

1 — У 'у

 

Величины х, у, I и

t могут быть сняты с карты изохрон с до­

статочной

точностью.

Если известна

величина ~vr,

то точность

3 заказ 1947.

33

определения координат х', у', z точек отражения S будет зависеть только от точности определения величины тх и ту.

Необходимо отметить, что практически по приведенным фор­ мулам следует вычислять только две координаты точек отражения, а именно z и х' или z и у'. Третью координату можно определить графическим путем, так как проекция точек отражения на днев­ ную поверхность будет лежать на направлении максимального

Рис. 13. Структурная карта залегания отражающего горизонта для пло­ щади, в пределах которой наблюдена карта изохрон, изображенная на рис. 12.

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

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

Для примера на рис. 13 приведена схематическая структурная карта, характеризующая залегание отражающего горизонта. Она построена по карте изохрон, показанной на рис. 12.

На рис. 13 ориентировочно заштрихована та площадь, строе­ ние которой может быть изучено при помощи карты изохрон, изображенной на рис. 12. Эта площадь примерно в 4 раза меньше той площади, на которой была наблюдена карта изохрон.

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

84

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

(1.27),

(1.28),

(1.29),

а

именно,

 

 

х' =

,

у'

=

,

z' = -|- |/(z?jt)2 — х2 — у2 .

Допускаемые в этом случае погрешности в определении ве­

личин

х', у' и

z

имеют следующий вид

(по Е. Н. Рязановой):

 

 

 

 

1

 

(e^-j-xy

 

 

Дж — 2

 

(х~Хо} — у(у~уо)'

 

 

 

л

_ _L,,

(г^-фаа/

х*

 

 

 

________ у«_________

У2 Уо (г^102— х (х—ж0) — у (у — у0) ’

 

Дг' =

V(vlt)2 — (х — ж0)2 — (у — уо)2 X

 

 

(г>1£)2 — х2— у2

- Х‘ - »’

 

 

 

~

 

При

глубине

залегания отражающей

поверхности

около

1000 м,

скорости

в покрывающей среде vr

2000 м/сек, углах

наклона 1,5° погрешности в определении х

и у' измеряются де­

сятками

метров,

а в определении z — единицами метров.

Определение конфигурации неплоской

отражающей

поверх­

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

вближайшем удалении влево и вправо от этой плоскости. При этом предположении можно считать, что точки отражения

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

получить из

уравнений

(1.34), (1.35), (1.36) и (1.41), положив

в них у = 0

и

Ту = 0.

Тогда

 

 

 

 

(у102—%2

Xt----i-[(2M)2 + a;2)]'t*

= х — ~ (г?! тж)

(1.42)

 

 

 

—(Vj тх)х

t—x xx

 

,

1

лГ»------ "l-------ST

W)2 —

*3

 

 

 

 

35

1

(v^ — x1

(1.43)

~ 2vt

t—x xx

 

tg(p =

x—v^t xx

(1.44)

 

vrt 1 — У2 Tx2

Конфигурацию отражающей поверхности можно определить по продольному годографу отраженных волн не только аналити-

Рис. 14. Построение отражающей площадки по ли­ нейному продольному годографу отраженных волн методом эллипсов.

чески по приведенным выше формулам, но и графически. Для ре­ шения этой задачи используются метод эллипсов (П. П. Лазарева),

метод окружностей

(Н. Н.

Пузырева), метод полей времен

(10. В. Ризниченко)

и метод

мнимой линии (В. Д. Завьялова,

Ю. В. Тимошина).

 

 

Метод эллипсов сводится к следующему. Имея годо­ граф отраженных волн и зная тем самым время ta прихода отра­ женной волны в некоторую точку А, лежащую на линейном про­ дольном профиле, а также скорость гх в среде, покрывающей отражающую границу раздела, находим длину пути пробега Sa отраженной волны из пункта взрыва О до отражающей поверх ности и обратно до точки A: Sa =

36

Взяв точки О и А за фокусы, строим графически эллипс с дли­ ной большой оси, равной величине Sa (рис. 14). Сумма расстояний от любой точки, лежащей на этом эллипсе, до точек О и А будет равна Sa. Таким образом, можно лишь утверждать, что где-то на этом эллипсе находится точка, лежащая на отражающей по­ верхности, но ее положение пока остается неопределенным.

Взяв на линейном профиле еще ряд точек В, С, D и т. д., удаленных на разные расстояния и в разных направлениях отно­ сительно пункта взрыва О, и зная времена прихода tb, tc, td и т. д. в эти точки волны, отраженной от одной и той же отражающей поверхности, определяем длины путей пробега Sb, Sc, Sd и т. д. отраженных волн до этих точек. Беря последовательно в качестве фокусов точки О и В, О и С, О и D и т. д., строим эллипсы с боль­ шими осями, равными соответственно Sb, Sc и Sd. Когда за точку приема упругих волн взята точка взрыва О, эллипс превращается в окружность с радиусом В = So/2 = vit0/2. На каждом из построенных эллипсов будет находиться по одной точке, лежащей на отражающей поверхности, от которых волны после отражения пришли в В, С, D и другие точки профиля. Положение криволи­ нейной отражающей площадки определится как огибающая к по­ строенному семейству эллипсов. Размер отражающей площадки, соответствующей наблюденному годографу, будет ограничен точ­ ками касания огибающей к эллипсам, за переменные фокусы которого были взяты крайние точки участка профиля, в пределах которого из пункта взрыва О был наблюден продольный годограф отраженных волн. Недостатком этого способа является трудо­ емкость построений семейства эллипсов. При определении близких к горизонтальным границ раздела вместо эллипсов можно строить дуги окружностей с центром, расположенным посредине между пунктом взрыва и точкой профиля (Л), которую в методе эллипсов брали за второй фокус, и радиусом

Ва = -1- /Л2 - (ОД? .

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

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

представлять собой

дуги

окружностей с радиусами Bi = viti.

Если между моментами времени

ti, выбранными для построения

фронтов

падающих

волн,

интервалы времени

A t = const, то

фронты

падающих

волн

будут

представлять

собой семейство

37

концентрических окружностей с равными расстояниями друг от друга.

Положение фронтов отраженных волн в различные моменты

времени Ц

может быть определено по заданному годографу на

основании принципа Гюйгенса. В момент времени

меныний ми­

нимального

времени £мнн прихода отраженной волны к какой-

либо точке

продольного профиля, фронт отраженной волны на-

Рис. 15. Построение отражающей площадки по ли­ нейному продольному годографу отраженных волн ме­ тодом полей времен.

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

ных из нескольких точек дневной

поверхности с радиусами

Rn = v^tn— ti), где

tn— время

прихода отраженной волны

в ту или иную точку

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

Установив положение фронта отраженной волны в момент

времени ti,

находим его положение в моменты времени ti— \t,

k — 2\t и

т. д.

Для этого из точек, лежащих па фронте волны

в момент времени

ti, следует провести дуги окружностей с ра­

диусами г = A tv±. По огибающей к новому семейству окружностей

38

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