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

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

.pdf
Скачиваний:
2
Добавлен:
24.10.2023
Размер:
11.24 Mб
Скачать

 

 

 

 

 

 

 

02

она

определяет величину

пара­

 

 

 

 

 

 

 

метров

pi=\kir\

 

и

рг=\Ыг\.

В

слу­

 

 

 

 

 

 

 

чае

однородного полупространства по­

 

 

 

 

 

 

 

грешность формул для зоны малых па­

 

 

 

 

 

 

 

раметров не превышает 20% при усло­

 

 

 

 

 

 

 

вии /?1=£П,5, над двухслойным

разрезом

 

 

 

 

 

 

 

при

О2ф0

и r>h\

та же

точность ап­

 

 

 

 

 

 

 

проксимации

реализуется

при

 

р2^1,о.

 

 

 

 

 

 

 

При

непроводящем

основании

послед­

 

 

 

 

 

 

 

нее

ограничение следует

заменить на

 

 

 

 

 

 

 

P s =

I k I

 

 

1,5.

Экспериментальным

 

 

 

 

 

 

 

критерием наблюдения поля в области

 

 

 

 

 

 

 

малых

параметров

является

пропор­

 

 

 

 

 

 

 

циональное изменение измеряемой ха­

 

 

 

 

 

 

 

рактеристики поля с частотой или, что

 

 

 

 

 

 

 

то же самое, неизменность значений о

 

 

 

 

 

 

 

при изменении

частоты.

 

 

 

 

 

 

 

 

 

 

Аналогичное

рассмотрение

особен­

Рис. 5. Двухслойные кривые

 

 

оЬ/а

ностей

электромагнитного

поля

над

 

 

двухслойной средой можно провести и

a — 0 2 / O i = 4;

б — a2 /<Ji = '/4.

Шифр

в дальней зоне (области больших па­

кривых — q=

I kihi I

 

 

 

 

 

 

 

 

 

 

раметров).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На основе

формулы

(III . 58)

при z =

/ i = 0

и с учетом

значений

первичного поля получим:

 

 

 

 

*=—h-k\r^2

 

 

 

 

 

k\r2 я

 

;

 

r=-i£rK

 

 

 

 

-67)

 

 

1 8

2

 

h

 

 

 

 

е

 

 

 

 

 

(1И

R определено в формулах ( I I I . 5 3 ' ) — ( I I I . 5 4 ' ) . Вводя в соответствии

с выражением

(III . 64)

значения кажущегося

сопротивления

и ис­

ключая таким образом из формулы

(111.67)

множители, не несущие

• полезной информации, будем иметь:

 

 

 

 

 

 

 

 

 

 

 

 

 

Р г = Р г = Р 9 = Р 1 | # 1 2 .

 

 

 

 

(Ш.67')

Из полученных выражений следует, что в рассматриваемой об­

ласти изменение разноса установки г не дает полезной

информации

о геоэлектрическом разрезе. Изучение его в дальней зоне может базироваться только на изменении частоты поля и на использова­ нии явления скин-эффекта:

при

 

 

« o - ~ 0 ( | A , A i l ^ 0 ) / ? — и

Р —

при

 

 

со—>• оо (| kxhx | —*• оо)

>1

и р—»"Pl.

Из формулы ( I I I . 6 7 ' ) следует, что амплитудные характеристики

любых компонент электромагнитного

поля

вертикального магнит­

но

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

юз Vю 1Lf

hi

hi

для ряда значений

рг/pi, приведе­

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

скольку

последние

определяются не

модулем

функции

I R \, а ее дейст­

вительной и мнимой частями. На

рис.

6 приводятся также

кривые

^ а / ь ,

заимствованные~ u l l „ „ . v - «изw

работы

Г. В. Приса [47]. Из сопоставления их с кривыми pz следует, что они бо­ лее быстро (при меньших значениях

Рис. 6. Двухслойная палетка частотных дипольных индукци­ онных зондирований (область больших параметров)

1 - 1

иЫа'

Шифр кривых —

Р2/Р1

 

 

 

-77-) выходят

на асимптоту второго слоя

(особенно при

p 2 > p i ) .

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

при измерении

~ ~ сокращается диапазон

частот,

необходимый для исследования

па

 

 

слоистого разреза.

 

При уменьшении разносов или частоты

условия дальней зоны

могут нарушаться. В этом случае кривые р начинают отличаться от кривых дальней зоны и не выходят на асимптоту второго слоя. Это

видно, например, из рис. 7, на котором изображены кривые р для

Рис. 7. Двухслойные кри­ вые рь/а

Шифр кривых — r / A i

8 76 32 64-

61

дальней

зоны и для

некоторых

конкретных отношений

Изме­

нение

hi

при заданных значениях

hi

определяет величину

пара-

 

r

 

 

 

 

 

метра р\ = I kir I и косвенным

образом

сказывается на

выполнении

условий

 

\kir\ » 1 и

\k2r\ ^>1 . Экспериментальным критерием

того,

что наблюдения поля проводятся в дальней зоне, является неизмен­ ность значений р, вычисленных по формулам ( I I I . 6 4 ) , при измене­ нии разноса установки.

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

структуры,

может

быть

выражено в

табулированных

функциях

в случае S-проводимости

при дополнительном

предположении о"г =

= 0 («зона

S » ) . На

основе формулы

(III . 52)

функция

отражения

в зоне 5 может быть записана в виде

 

 

 

 

 

 

 

 

 

(111.52')

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

 

 

оо

 

 

 

?(Ps)=-M-$

 

Л(Хр)

~ — - ( / > I = |

 

k\\h,r^Sr),

исследованную и

табулированную в работе

Л. Б.

Гасаненко и

Г. П. Шолпо [17]:

 

 

6

 

4

р%

F(ps),

 

ev =

^

hr=~i-^-F'{ps),

--?fF'(ps).

 

 

 

 

 

0И.68)

Принятые при выводе этих формул предположения практически обозначают, что разнос применяемой установки г значительно пре­ вышает мощность проводящей толщи hi(r^>hi) и что в этой толще не происходит поглощения поля (|kihi|->•0, но \k\\hir конечно).

При малых p s (низкие частоты)

Ps

 

 

 

 

 

J - g - ( I n / ? | - C ) ,

 

 

 

 

 

p i

Ps

А

 

•I

Pi

1

Ps 1

-

4

2

 

 

 

1 "

 

 

 

при больших p s

(высокие

частоты)

(1пр| — С — 1 )

и

36

 

,

 

. 6

12

Здесь С=0,5772-

Ps

 

 

 

Ps

Ps

постоянная

Эйлера.

 

(111.68')

(Ш.68")

62

Рис. 8. Графики магнитного поля

вертикального магнитного диполя

над двухслойной средой (5 — про­

водимость):

 

 

/ — R e f t z + 1 ;

2 — Imhz;

3 — Re hr\ 4 —

lmhr;

5-\ha\;

6-\hb\;

7-\hb/ha\;

Я — Ф„ .

град;

9 — ya,

град

 

Графики

магнитного

по­

 

 

 

 

ля

(ортогональные

компо­

 

 

 

 

ненты

и

эллиптические

ха­

 

 

5'

 

рактеристики)

для

рассмат­

90°-ho

 

риваемого случая

приведены

80

7

 

 

на

рис.

8. На

основе

этих

5

7

 

графиков или непосредствен­

70 '

3

 

 

но

 

на

основе

 

формул

ВО- 2

 

 

(111.68) — (111.68")

может

50 ю-'

97

8

быть

найдена

 

продольная

 

 

проводимость

 

5

толщи

по­

• 7

 

 

верхностных

отложений,

пе­

30 -5

 

 

 

рекрывающих

непроводящее

20 . 3

 

 

основание. Из

сопоставления

10 -

2

 

 

формул

(111.68')

и

(111.65)

0 Чо~г

 

 

следует,

что

основное отли­

 

W

2 3 5 7 10 2 3 S

7 70 2 Зр<

чие

случая

S-проводимости

 

 

 

 

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

Из сопоставления формул (111.68") с формулами (111.65), ха­ рактеризующими поведение поля над двухслойной средой в дальней зоне, следует, что в случае 5-проводимости при больших парамет­ рах p s поле убывает с частотой на порядок быстрее. В этом сказы­ вается дополнительная идеализация рассматриваемой электриче­ ской ситуации (| 1-ИЗ, | &i | г - » - оо) . При дальнейшем увеличе­ нии частоты, когда предположение | kih\ | ~ 0 станет неправомочным, поведение поля будет определяться формулами для однородного полупространства с проводимостью Oi в дальней зоне.

Магнитное поле диполя, приподнятого над слоистой средой

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

63

а

определяется в

основном

высотой

установки z + h. Поэтому

при

z + h>p становится целесообразным

вводить магнитные числа

вто­

ричного поля соотношениями

 

 

 

 

 

 

 

"

1 ~

М

/г)3

и*

 

 

 

 

4тс(г +

1 1 1

 

 

и

рассматривать

их как

функции

параметров q2

=0JU.CD(z + h)2 и

q2

= oi\xahi(z + h).

Такая система обозначений формально отражает

то

обстоятельство,

что в

аэроварианте

дипольных

индукционных

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

При выполнении условия z + h>p величина наблюдаемого вто­ ричного поля резко уменьшается и доступными для измерения ста­ новятся главным образом его мнимые части, которые легче выде­ лить на фоне большого первичного поля. В связи с этим основное внимание ниже будет уделено анализу мнимых частей ортогональ­ ных компонент вторичного поля.

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

предельных

случаях.

При достаточно

большой высоте

установки

(z + h^$>p)

 

на

основе

формулы

(3.68), приведенной

в

работе

И. С. Градштейна и И. М. Рыжика [19], могут быть найдены

выра­

жения для поля диполя, приподнятого

над однородным

полупро­

странством:

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

В этих формулах и = — t f i q u

Hi(«)

и Ni(u) —функции

Струве

и Неймана.

 

 

 

 

 

 

 

На рис. 9 представлены графики мнимой части магнитных чи­

сел h* и h* для

предельного случая ( I I 1.69) и для ряда значений

Z

р

 

 

 

 

 

 

 

z + h

 

 

 

 

 

 

 

 

d— , найденные численным интегрированием формул (111.56). тт Р , z+h . .

Из сопоставления этих графиков следует, что при -s>4 по­ грешность аппроксимации нормального поля приподнятого диполя асимптотическими формулами (111.69) не превышает 20% . При­ мерно такая же точность аппроксимации имеет место и в общем случае слоистого разреза, поэтому в дальнейшем будем рассматри­ вать только предельный случай:

z + h

>- оо.

Р

В области достаточно низких частот (малых параметров) на ос­ нове формул (111.57) нетрудно записать обобщенное выражение для

64

Рис. 9. Поле вертикального маг­ нитного диполя, приподнятого над однородным полупространством

 

 

г + А

 

 

Шифр кривых — d =

 

 

 

вторичного поля:

 

 

 

Im

h*z=

4fz

 

 

z + h

Im hu

4jt

 

(111.70)

 

 

 

 

где

 

л f

 

 

 

Я/г

 

1

+

 

 

 

 

 

 

 

1

 

 

 

 

 

2?

 

 

(111.70')

 

1 + 2/£

 

 

 

 

 

4\

] /

\ +

g

+ 4(5

 

V2

 

1 +

£/2

 

Yi+ijg-

 

 

i

 

 

 

4(5

 

 

 

 

 

~ q s '

1+2/g

 

*

В этих формулах

(Ш.70")

 

z + h

Ai

Dsxhx

с2 (г + /г)

В

соответствии с

( I I I . 7 0 ' ) — (111.70")

обобщенные

параметры

Яft и

ЯТР обладают

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

в

параметры

qi/i2{hi-+oo),

q2/i2(hi^0)

и q8(hi-+0,

a i - > o o ,

5 =

const, a2 ->-

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

При увеличении высоты

установки z + h вертикальная составляю­

щая вторичного поля I m Hz

от однородного полупространства убы-

1

 

 

 

 

 

 

вает как —— г - >

а

от

 

пласта

конечной

мощности — как

z+h

 

 

 

 

 

 

(z +1 hyT' Убывание

I m # p

в

обоих

случаях происходит на поря­

док быстрее. Разное

убывание

с высотой поля

от ограниченных и

5 Заказ № 271

65

неограниченных пластов является предпосылкой для проведения

высотных зондирований в этой

области.

 

 

 

 

 

 

В области

больших

параметров (достаточно

больших

высот —

z + h^hi) на основе формулы

(III . 58) также можно записать обоб­

щенные выражения для мнимых частей вторичного поля:

 

 

 

 

т

ь*

1

2

 

z + h , , *

24

 

/ I I T _ 1 Ч

 

 

 

\mhz=—j-;

 

 

— I m / z p = — т

;

 

 

(III.71)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

1 + m2e~2

V l q

2OTCOS Y'2q • e / 2 <

7

 

 

 

 

 

1 _ O T 2

-2 Y-2Q + 2 M

S I N

У2д-е~У2"

 

 

 

 

 

 

 

 

2

1 +

V2* +

 

 

 

 

 

 

 

 

- w - V b

+ v b •

 

 

( 1 1 , - 7 1 ' )

Здесь m = - ~ ^ : — < 7 = U i f t i | ,

a =

a i A i l / - ^ - .

 

 

 

 

УР2+УР1

 

 

 

 

 

 

 

0"2

 

 

 

 

 

Из приведенных

формул следует,

что

обобщенные параметры

q 2 h

при q-^oo

переходят в

i2qit

при ^ - > 0 в У2 ?2 и при

а2 ->-0,

но

OdAi = 5 = const — в

<7|.

Параметр

д 2

строго

 

пропорционален

высоте установки \z + h\; это приводит к тому, что мнимая часть вертикальной составляющей поля независимо от строения среды убывает с высотой в рассматриваемой области как 1/(г+/г)4 , а мни­ мая часть горизонтальной составляющей над любым слоистым раз­ резом — как l/(z + h)5. Отсюда следует, что высотные зондирования в области больших параметров невозможны. От частоты параметр <72 , а значит и вторичное поле, зависит сложным образом. При

w ->• 0 и со оо поле становится таким же, как над однородным по­ лупространством с сопротивлением, соответственно рг и pi. Это предопределяет возможность осуществления частотных зондирова­ ний в этой области.

Таким образом, при дипольных исследованиях в аэроварианте дело обстоит так же, как и в наземном случае: в области достаточно низких частот можно проводить высотные зондирования, а в об­ ласти достаточно больших высот —частотные. Малость частоты определяется требованием <7г<С1, а большая высота — условием

На основе формул (III . 70) и ( I I I . 7 1 ) , полагая в них qfz — q f P = = qil~}j2 или q 2 h = ~tf 2 q i y можно ввести понятия кажущейся прово­ димости а или сопротивления р в аэроэлектроразведке:

~

4

*

81m h*

 

° = - ^ u r w l m h * = w e U ) ;

( ш - 7 2 )

V 7 = V

^ ^ + h )

Im Н1=^Щ±^-

 

Im h*P. (Ш.73)

66

В соответствии с этими формулами а и р могут быть определены из результатов измерения поля. Палеточные кривые о в зависимо­ сти от высоты установки (z + h/hi) и р — в зависимости от частоты (параметра q = \kihi\) могут быть рассчитаны по формулам

-2*-=2^Ц-,

^ = 4 - ^ ;

(Ш.74)

^=У2-ЩГ. (111.75)

Результаты этих расчетов изображены на рис. 10 и 11, предста­ вляющих собой двухслойные палетки высотных или частотных зон­ дирований (рис. 10 — только для о г ) . Их вид и способы использо­ вания при интерпретации принципиально не отличаются от соот­ ветствующих палеток для наземных зондирований. Эксперимен­ тальными критериями правомочности применения этих палеток для

интерпретации являются: постоянность величины о при изменении

частоты и р—при изменении высоты установки.

 

Если не задаваться ограничивающими требованиями {qi^$>l

или

z + h^>hi), то достаточно компактные палетки для поля диполя

над

двухслойной средой могут быть получены лишь для некоторых ча­ стных случаев. Наиболее важными среди них являются: 1) тонкий пласт с заданной продольной проводимостью S, подстилаемый сре­ дой конечного сопротивления; 2) пласт конечных мощности и про­

водимости,

лежащий на

непроводящем

основании. На

рис. 12, с и

13, а приведены графики Im/г* , построенные для обоих

упомянутых

случаев в

зависимости

от параметра

<7/. Поскольку частота поля

входит в этот параметр в виде отдельного множителя )/(о, эти кри­ вые хорошо приспособлены для интерпретации измерений на раз­

ных частотах. При продольной

проводимости 5

кривые

индексиро-

 

D

 

aiht

——,

 

а параметр qfZ =

 

if

,

,

1

п

ваны значением 6 = — ;

 

*

qsv

1 + - ^ — -

Д л я

 

 

o2(z

+ h)

r

г

, 0 Г

Г

 

26

 

 

 

(Г9

2 + П

 

 

 

 

 

26

 

пласта

ной

м

 

 

 

 

 

 

z+h

 

,

a

<7/z

конечной

мощности шифр

кривых — g =

 

<7i

1

 

 

 

 

 

 

 

hi

 

 

 

У2

У 1 + £ / 2

. Левые ветви всех кривых на этих палетках сли-

 

 

 

 

 

 

 

 

 

 

 

 

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

5*

67

0,07

0,01 2 3 5 70,1 Z 3 5 71,0 Z 3 5 7 70 2 3 5 7100 Z

Рис. 10. Двухслойная палетка высотных зон­ дирований (z+h>p)

Шифр кривых — y=a2faj

J 5 70,07 2 3 5 70,1 Z 3 5 77,0 2 3 5 7

Рис. 11. Двухслойная палетка частотных зонди­ рований ( z + / i > р)

Шифр кривых — p2 /Pi

Jmh-,

1,0

 

 

 

 

 

 

 

 

 

 

 

 

 

7

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

16

 

 

3

 

 

 

//

 

 

 

 

 

 

 

 

 

 

•1/8

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

,0,1

 

 

 

1/

 

 

 

 

 

 

 

 

 

 

 

 

Т/г

 

 

 

 

 

 

 

 

 

 

 

i-

f/4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

fii

 

 

 

 

 

 

 

 

 

5

 

 

 

Hi-///6

 

 

 

 

 

 

 

 

3

 

 

'kt

 

 

 

 

 

 

 

 

 

 

41 h

 

 

 

 

 

\

 

 

2

 

 

 

 

 

 

 

 

R01

 

 

f

 

 

 

 

 

 

 

\

 

 

0,1

2

3

5

71,0

2

3

5

7 10

2

3

5

7 ^

Рис.

12.

Графики поля вертикального диполя

 

Im hz*

 

над тонким проводящим

пластом

(02 ф 0)

 

а — в

зависимости от

параметра

 

qfz

(шифр

кривых —

Р= Oj"i (z^'+и\п) j\ в зависимости от параметра Яи (шифр

кривых — a = ( T i f t i l /

У""

У

о2

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