![](/user_photo/_userpic.png)
книги из ГПНТБ / Светов Б.С. Теория, методика и интерпретация материалов низкочастотной индуктивной электроразведки
.pdf
|
|
|
|
|
|
|
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 |