![](/user_photo/_userpic.png)
книги из ГПНТБ / Богданов, Алексей Иванович. Интерпретация сейсмических годографов
.pdf■/' |
= У + 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