Устройство для измерения естественных потенциалов в скважине

 

Союз Советских

Социалистических

Республик

О П И С А Н И Е „„gyg165

ИЗОБРЕТЕН ИЯ

К ПАТЕНТУ (6!) Дополнительный к патенту (51) М. Кл. (22) Заявлено 26.05.71 (21) 1601490/1656424/

/18-25 (23) Приоритет 03.12 70 (32) 03.12.69

G 01 Ч 3/18

Государственный комитет

СССР не делам изобретений и DTKpblTNN (31) 881793 (331 США

Опубликовано 05,08.79. Бюллетень ¹ 29 (53) УДК

550.832 (088.8) Дата опубликования описания 05.08.79 (72) Автор . изобретения

Иностранец

Ник Шустер (США) Иностранная фирма

"Сосьете де Проспексьен-Электрик Шлюмберже" (Франция) (71) Заявитель (54) УСТРОЙСТВО ДЛЯ ИЗМЕРЕНИЯ ЕСТЕСТВЕННЫХ ПОТЕНЦИАЛОВ

В СКВАЖИНЕ

Изобретение относится к технике, применяемой при геофизическом каротаже, например к устройствам для измерений естественных потенциалов (SP) в скважинах.

Известно устройство для измерения естественных потенциалов, возникающих на изолированном долоте. При измерениях используется мостовая схема с автоматической балансировкой моста (1).

К недостаткам устройства относятся невысокая надежность и небольшая точность измерений.

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

Зонд представляет собой электрод, естественный потенциал которого измеряется относительно потенциала электрода, устанавливаемого на поверхности земли у устья скважины (21.

Естественный потенциал зонда создается в результате диффузии солей из пластовых вод в буровой раствор или из бурового раствора в пласт и адсорбции их ионов поверхноспю

Ь горной породы; фильтрации пластовых вод нли бурового раствора; окислительно-восстановительных реакций, протекающих преимущественно на поверхности соприкосновения горных пород с окружающей их средой (другими по10 родами или буровым раствором).

Величина измеряемого потенциала SP осложняется помехами разнообразного; происхождения . Во-первых, величина измеряемой разности потенциалов поляризации пород искажается

15 поляризацией электродов и ее изменением при движении электродов и во времени. Во-вторых, процессы осаждения твердых частиц в буровом растворе искажают результаты измерений. Гальванокоррозия грузов, а также естественные и

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

679165

10

3 те взаимодействия брони кабеля, бурового раствора и обсадных труб создаются гальванические элементы, ток которых изменяется вследствие нестабильного контакта брони кабеля с трубой. Воздействие помех на измерение потенциала SP сказывается еще более при исследовании морских скважин, так как трудно изолировать поверхностный электрод из-за высокой проводимости морской воды. Морская вода способствует образованиюгальванических элементов между различными металлами, ток которых также непостоянен во времени. Указанные причины делают результаты измерения потенциала SP малодостоверными.

Для устранения воздействия указанных помех предлагалось измерять дифференциал или градиент SP и интегрировать эту разность потенциалов. Однако шум постоянного тока или нестабильность усилителя приводит к значительным ошибкам.

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

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

Дель изобретения — повышение точности измерения потенциала SP путем исключения шумов высокой и низкой частоты.

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

Устройство позволяет измерять потенциалы

$Р на каждом скважинном электроде и устранять шумы, оеложняющие эти измерения.

Зонд, состоящий из двух подвижных либо фиксированных электродов, расположенных на скважинном приборе, обеспечивает измерение потенциалов SP каждого электрода и их градиента, Градиент потенциала затем интегриру-. ется. Высокочастотные составляющие P сигналов, измеренные между электродом зонда и броней кабеля, в сочетании с низкочастотными составляющими сигналов, измеренными между электродом зонда и поверхностным электродом, позволяет получить неискаженное

4 значение естественного потенциала $Р в сква. жине.

На фиг. 1 представлена структурная схема предложенного устройства; на фиг. 2 — частотная характеристика измеряемых сигнала и шумов; на фиг. 3, 4 — суммирующая схема, варианты; на фиг. 5 — характеристика фильтров; на фиг. 6, а и б — эпюры сигналов генератора и хронирующего устройства, на фиг. 7 — датчик в виде тороида для измерения; на фиг. 8— измерительная схема, вариант.

Устройство погружено в скважину, пересе, кающую горные породы и заполненную буровым раствором. Зонд 1 с двумя электродами

2 и 3, раснесенными на расстояние а, проводами 4, 5 подключен к входу дифференциального усилителя 6; Электроды 3 и 7 соединены с входом дифференциального усилителя 8. Выходы усилителей 6, 8 подключены к входам низкочастотного фильтра 9 и высокочастотного фильтра 10, Высокочастотный фильтр содержит фильтр нижних частот 11, (ФНЧ),. операционный усилитель 12, фильтр верхних частот (ФВЧ) 13, интегрирующий конденсатор 14, резистор 15 и выключатель 16.

Выходы фильтров 9, 10 подключены к суммирующей схеме 17, содержащей суммирующие сопротивления 18, 19, суммирующий усилитель 20, резистор 21 и регистратор 22.

Устройство работает следующим образом.

РазноСть потенциалов электродов 2 и 3 измеряется дифференциальным усилителем 6, разность потенциалов электродов 3 и 7 — усилителем 8. Схема 17 обрабатывает полученные сигналы.

Прежде чем рассмотреть работу схемы 17, рассмотрим частотный состав приходящих сигналов. Обычный сигнал SP, измеренный усилителем 8, содержит шумы высокой и низкой частот, Шум высокой частоты создается магнитотеллургическими токами, биметиллическими элементами, шумом буровой установки и т. д.

Постоянный ток или шум низкой частоты обусловлен поляризацией электрода и мало влияет на результаты измерений. Таким образом, выхощчой сигнал усилителя 8 может быть обозначен У + И, где Ч вЂ” истинное значение SP (шум отсутствует), а N — составляющая шума высокой частоты (см. фиг. 2, где по оси абцисс указан диапазон частот, а по оси ординат — уровень сигнала $Р) .

Сигнал ESP усилителя 6 зачастую свободен от шума высокой частоты, поскольку источники этих шумов удалены от электродов 2 и 3 и одинаково.=," влияют на оба электрода, поэтому Ь$Р не содержит этих шумов (они взаимно уничтожаются на усилителе 6), По679165 этому выходной сигнал усилителя 6 (Ь$Р, или градиент) содержит истинный сигнал градиента, обозначенный через G, и шум поляризации, обозначенный через Р, т.е. ESP = G + Р.

Фиг. 2 иллюстрирует зависимость мощности

W от частоты для составляющих измеряемого сигнала. Сигнал SP является суммой Ч + N, а сигнал Ь$Р— суммой G + P, поэтому из сиг нала усилителя 8 можно отфильтровать низкочастотную составляющую сигнала, а из сигнала усилителя 6 — высокочастотную составляющую сигнала SP.

Описываемое устройство позволяет измерить сигналы SP и Ь$Р в схеме 17 и получить сигнал SP, в котором скомпенсировано влияние помех. Другими словами, составляющие N и Р сигналов отфильтровываются от компонент Ч и G в цепи обработки сигнала.

Выходной сигнал градиента Ь$Р на выходе усилителя 6 выражается следующим образом: и

gP =(:,+Р =Ч (2+a)- 4Ю)+у= а — +Ð(1

J7 ) / где Ч(Е) — истинное значение SP на глубине

2 в скважине;

Ч(Е 4 ) — истинное значение $Р на глубине

2 + а в скважине.

В функции от времени уравнение (1) имеет вид:

ASP=O - +p = — (+p ) (2) где V = - - скорость движения электродов.

Записывая уравнение (2) через оператор Лапласа, имеем:

Т. (5 й) = 8 г Ч(й) + Р, (3)

С (g) где г — параметр времени;

S — комплексное число.

Тогда на выходе усилителя 8 имеем:

ggp=V(Z)- Y -И. (5) где Ч)) — потенциал электрода 7.

Обычно Чд равно нулю без учета составляющих шума N, и уравнение (5) записывается. следующим образом фР = V(Z) +N, (6) Преобразование Лапласа для уравнения (6) дает

1)

Для простоты вместо Ч(2) будем писать Ч.

Из фиг. 2 следует, что оба члена уравнення (7) отличаются по частоте. Поэтому выражение для выходного сигнала схемы 17 должно иметь следующий вид; ,. e,= 1-Ю+Я Н(Ят.Ч+Р) (a) или л ) Г"

Е, =g. + g т Н) V+ .t.М+ Í р, (9) где 1 и Н вЂ” соответственно низко- и высокочастотные функции.

1О

Усилитель 8 соединен с ФНЧ 9, имеющим переходную характеристику 1 уравнения (9).

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

21 в цепи обратной связи.

Аналоп4чно сигнал DSP усилителя 6 подается на вход операционного усилителя 12, снабженного резистором 15 в цепи обратной связи, т.е. сигнал ЬЯР поступает в высокочастот1ный фильтр 10, переходная характеристика которого описывается функцией Н иэ уравнения (9). ФНЧ 11 имеет характеристику 1 и пропускает только низкочастотную часть сигнала

Ь$Р к вычитающему входу усилителя 12, к положительному входу которого подается сигнал Ь$Р or усилителя 6. Усилитель 12 с конденсатором 14 образует интегратор, выход которого соединен с ФВЧ 13, имеющим переходную характеристику Н. Выходной сигнал ФВЧ .13 через усилитель 20 записывается регйстратором 22.

Функции схемы 17 выражены уравнениями (8) и (9). Иэ уравнения (9) условие точного воспроизведения истинной части Ч сигнала SP выражается как:

L,=ßñÍ =Х. (10)

Переписывая уравнение (10) в функции Н поH ™i. («) зо

Для фильтра 10 переходная характеристика Н выражается следующим образом: (<).)

$г

Иэ уравнения (12) видно, что если L = 1 и

Н) = 1, то уравнение (12) аналогично уравнению (11), а функция Н фильтра 10 будет удовлетворять уравнениям (10) и (11). Если Н=1, 4 то ФВЧ 13 не нужен. Однако несовершенство интегратора, уход нуля усилителя 6 и тд. приводят к тому, что в фильтре 10 практически нужен фильтр 13, который компенсирует лю-. бые помехи в фильтре 10. С друтой стороны, 45 интегратор не обеспечивает точного воспроизведения составляющих сигнала постоянного тока из-за включения в фильтр 10 ФВЧ 13.

Включение переменного резистора 15 параллельно конденсатору 14 с помощью выключателя

16 обеспечивает установку постоянной времени интегратора. Если допустить, что уравнение (10) не выполняется, то и в этом случае устройство обеспечивает измерение истинного сигнала Ч.

Схема 17 осуществляет измерения в функции временной, если зонд перемещается в скважине с постоянной скоростью. В случае нарушения этого условия регистрируемый сигнал осложняется шумом высокой частоть1.

679165

50 к (+gT f+gT 1+«Мт

L 2 или

Н- <р тг т1(рт ) 55

Уравнение (14) означает преобразование Jlanнаса для фильтра нижних частот, а уравнение (l6) — для полосового фильтра, 7

Компенсация этого шума обеспечивается интегратором; интегрирующим сигнал в функции глубин. Кроме того, изменяют фильтры 9, 11, 1 3, .

На фиг. 3 показана суммирующая схема, реализующая коррекцию сигнала при изменении скорости движения зонда в скважине, Сигнал ASP подается на вход усилителя 23, коэффициент усиления К которого изменяется в зависимости от изменения скорости движения зонда в скважине. Вращающийся ролик

24 фиксирует это изменение, Ролик 24 передает вращение валу 25, который управляет коэффициентом К усилителя 23 в функции скорости передвижения кабеля.

Усилитель 23 содержит входное сопротивление 26, соединенное со входом операционного усилителя 27, другой вход которого соединен с поверхностным электродом 8. В цепи обратной связи усилителя 27 включен потенциометр

28, управляемый валом 29 измерителя скорос.ти 30. Измеритель скорости работает в функции от скорости вращения вала 25. Таким образом, изменение коэффициента К определяется изменением скорости перемещения кабеля.

На вход интегратора 31 поступают сигналы от усилителя 8 и 23. На выходе интегратора

31 включен ФНЧ 32, имеющий переходную функцию -- †- —, где T> — постоянная времени фильтра 32. Обратная величина выходного сигнала ФНЧ 32 и входной сигнал PS суммируются интегратором 33, вйход которого соединен со вторым ФНЧ 34, имеющим переходнУю хаРактеРистикУ -т, ale T + Га постоянная времени фильтра 34. Выходные сигналы фильтров 32 и 34 суммируются интегратором 35 и регистрируются регистратором 36, который приводится во вращение валом 25 в функции перемещения кабеля. Поэтому компенсированный сигнал PS записывается как функция глубины скважины.

Для переходных характеристик L и Н схемы на фиг. 3 определены значения: (g — — ) + ((33 4s

1 Х 1

=, 7 + т + 8T2 или

3+8 (Т1+Т2)

У

У+От)(572) тогда

Из уравнений (14) и (16) выражение L+StH согласно (9) можно записать как: +о(т,+ У $ К т2 (+ $ т,) (+ T )

На фит. 5 показана зависимость амплитуд сигналов от частоты для переходных характеристик 1, Н и 1 + 8TH — см. уравнения (14), (16) и (17) соответственно. Эти характеристики описывают фильтр нижней частоты, полосовой фильтр и фильтр более высокой частоты, чем частота предшествующих фильтров.

Иэ уравнения (9), фиг. 2 и 5 следует, что путем выбора постоянной времени устройства, изображенного на фиг. 3, можно исключить все составляющие N шума высокой частоты.

Переходная характеристика 1 на фиг. 5 может быть такой, чтобы снизить до нуля этот шум. Составляющая P поляризации также может быть исключена, так как она уменьшается до нуля при частоте, меньшей частоты диапазона, который охватывает переходная характеристика Н. Следует заметить, что составляющая Ч сигнала PS лежит в области частот, где может быть выбрана характеристика 1 + ЯсН", поэтому выбор характеристики 1 + В И должен быть таков, чтобы ее частота и частота сигнала Ч не совпадали. По уравнению (10) для точного воспроизведения составляющей Ч сигнала PS необходимо выполнение этого условия во scar диапазоне частот измерения составляющей. Таким образом, устройство на фиг. 3 должно работать в диапазоне частот изменения сигнала.

Если зонд перемещается в скважине с постоянной скоростью, то устройство на фиг, 3 обеспечивает точное измерение сигнала PS. Если же скорость перемещения зонда изменяется, то это устройство позволяет компенсировать влияние этих изменений на результаты измерений. Для этого сопротивление потенциометра 28 изменяется в функции от скорости зонда таким образом, чтобы член Kt в уравнении (17) оставался постоянным. Для этого необходимо, чтобы К.г =Т . Тогда уравнение (17) преобразуется в следующее:

g- 8(т +7 )+pат т

f„+ «1g f- — . (18) (r S ) (1+ Т )

Иэ уравнения (18) видно, что числитель равен знаменателю, и Ь+бтН =М в соответствии с уравнением (10).

При выборе разноса электродов зонда (расстояния а) и постоянных времени Т1 и Т, необходимо учесть ряд факторов.

Необходимо обеспечить замер градиента и получить сигнал, достаточный для регистрации ) 679165

9 дифференциальным усилителем 6. Компромиссное решение достигается при а = 0,3 г 0,6 ь1

Для выбора Т, и Т (Т, = Тр) необходимо обеспечить такое время, чтобы характеристика достигла нуля до того, как сигнал М будет существенно влиять на результаты измерений. Величина времени должна быть такая, чтобы переходная характеристика Н была бы равна нулю при частотах, на которых состав о ляющая P сигнала PS существует. Кроме того, скорость перемещения электродов определяет выбор Т, и Т . Для скорости 0,3 м/с приемлемые результаты получают при постоянных времени от 5 до 40 мкс.

На фиг. 4 дана электрическая схема устройства для изменения коэффициента усиления К как функции скорости перемещения зонда.

Провода 4, 5 соединяют электроды 2 и 3 соответственно через сопротивления 37 и 38 со входами операционного усилителя 39, на выходе которого сигнал пропорционален BSP.

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

39 соединен с электродом 7 сравнения. Усилитель 39 охвачен обратной связью по сопротивлению 40. Через резисторы 41 и 42 выходной сигнал усилителя 39 подается на вход усилителя 43.

Потенциал нижнего электрода 3 по проводу .4 подается через входное сопротивление 44 на вход усилителя 43. Цепь обратной связи этого усилителя содержит емкость 45 и резистор 46. Постоянная времени этой цепи рав-. на Т,, Выходной сигнал усилителя 43 и потенциал SP, подаваемый по проводу 4, суммируются на входе усилителя 47 через два входных резистора 48 и 49 соответственно. Оба других входа усилителей 43 и 47 соединены с электродом 7. Усилитель 47 содерх.ит цепь обратной связи, составленную резистором 50 и ем- 4, костью 51, постоянная времени которой Т .

Выходные сигналы усилителей 43 и 47 суммируются через два сопротивления 52 и 53 иа входе операционного усилителя 54, на выходе которого подключен регистратор 55 для записи 4 компенсированного сигнала SP в функции глубины скважины.

Переключающий триод 56 соединяет точку, общую для резисторов 41 и 42, через периодически повторяющиеся промежутки времени, 50 определяющие скорость перемещения кабеля, с землей, Для определения этой скорости импульсный генератор 57, содержащий вращающийся ролик 58„взаимодействующий с кабелем, вырабатывает импульсы, соответствующие изменению величины перемещения кабеля.

Эти импульсы возбуждают кипп-реле 59 с постоянным временным периодом, которое управляет работой триода 56.

Н Ри g (г1) к= „ иТ

Учитывая, что R Ть и r п о сeт о я н н ы, коэффициент усиления усилителя К будет функцией только и

Необходимо также устранить влияние составляющей P измеряемого сигнала. Согласно уравнению (9); член P поляризации увеличивается переходной характеристикой Н, которая для устройства на фиг. 4 представлена урави нением (16). Заменяя К на Т ° o поскольа ку КГ=Т и Т = —,, выражение для профильтрованного сигнала P приводят к виду; (.8, т,3

Я+ 6 Тт ) (f+ 3 Г )

Отсюда легко определить член Ри. Таким образом, если изменяется скорость и и P имеет уровень постоянного тока, отличный от нуля, то измеряемый сигнал 1О будет искажен в результате действия переходной характеристики высокой частоты (2Z) 5"Т, Т /C(+ Sr,)(+ят ц

Во избежание этого во, входном канале, по которому поступает сигйал SzV+P цепи обработки сигнала (фиг. 4 или 3), устанав ливают фильтр верхних частот (тиристор 60), устраняющий влияние постоянной составляющей P поляризации сигнала SP. Однако этот

На фнг. б, а и б, показано изменение импульсов генератора 57 и хронируюших импульсов кипп-реле 59.

Импульсы генератора 57 могут изменяться по частоте, зависяшей от скорости движения кабеля. Кипп-реле 59 вырабатывает на каждый из импульсов генератора 57 сигнал определенной длительности Т . Эти хронирующие импульсы выключают триод 56, так что величина сигнала на входе усилителя 43 изменяется пропорционально скорости движения кабеля.

Если генератор вырабатывает импульсы через каждый отрезок "в" перемещения кабеля, то временные интервалы между последовательными импульсами будут равны —,„. Учитывая, что общее сопротивление резисторов 41 и 42 равно т, величина среднего тока на входе усилителя 43 будет равна:

ygP Гд (19)

Icpeg =, г и д ри. т.ч

+cpeg (20)

Sr

Пусть сопротивление резистора 46 равно R, тогда коэффициент К усиления усилителя 43 равен:

679165

11 фильтр не нуЖен, если скорость движения скважинного снаряда стабилизирована.

Если переходная характеристика фильтра 60 равна - -„+, то коэффициент усиления К

7 имеет следующее выражение: ь: = И (u,T) (g3), 3-ST„a i

Введение фильтра 60 в устройство на фиг. 4 требует уточнения переходных характеристик

° 1 и Н. Однако, если. Тз>> Т,, Т, то принятые значения 1. и Н могут быть сохранены, при этом будет наблюдаться незначительное отклонение коэффициента усиления для L +бай, приведенное в уравнении (10).

Следует заметить, что измерение тока сигнала SP можно использовать для получения гради ента сигнала вместо-применения двух расположенных рядом электродов. Для этого случая ток сигнала SP равен:

1 дч 6 Z (2Ч)

$r — н где G удельная проводимость раствора в скважине, А — объем жидкости, проходящей через датчик, выполненный в виде тороида; (3 А обычно постоянна, следовательно, результат измерения пропорционален, т.е.

dv

<ГЯ этот случаи аналогичен использованию двух электродов.

На фиг. 7 показано применение датчика, выполненного в виде тороида для измерения

SP.

Обычно SP измеряется электродами 3 и 7, включенными на вход усилителя 8, Градиент измеряется, например, индукционным датчиком 6l, подключенным к усилителю 6.

Выходные сигналы усилителей б и 8 обрабатываются затем по описанной схеме измерения компенсированного сигнала SP

Вместо градиента .сигнала SP для измерения высокочастотного сигнала может применяться сигнал, измеренный обычным образом электродом 3 и броней кабеля.

Указанную возможность иллюстрирует фиг.8.

Зонд 62, укрепленный на кабеле 63, содержит электрод 64. Нижняя часть кабеля покрыта изолирующим материалом 65. Провода бб и 67 соединяют электрод 64 и броню кабеля бЗ с измерительной схемой.

Усилитель 68 регистрирует разность потенциалов брони кабеля и электрода 7 сравнения.

12 Усилитель 69 аналогично регистрирует сигнал сравнения скважинного и поверхностного электродов. Высокочастотный сигнал SP с брони кабеля выделяется ФВЧ 70, Устройство 71 интегрирует эти сигналы, после чего они регистрируются самописцем 72 в виде функции глу-, бины скважины.

Усилитель 69 содержит низкочастотную составляющую сигнала, но искажен шумом пере1о менного тока. Наоборот, сигнал с выхода усилителя 68 содержит высокочастотную информацию SP, но искажен низкочастотным шумом.

Выделяя низкочастотную составляющую из первого измерения (ФНЧ 73) и высокочастотную (ФВЧ 70) — иэ второго, можно получить неискаженное значение сигнала $Р.

Следует указать, что реализация устройства возможна в виде поверхностной или скважинной схемы.

Вместо выделения- частотных составляющих сигнала SP с помощью описанного устройства возможна регистрация указанных составляющих SP с целью последующей обработки информации вычислительными машинами. В

25 этом случае применение описанного устройства не рекомендуется.

Формула изобретения зо

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

Источники информации, принятые во внимание при экспертизе

1. Патент США ND 2973505, кл. 340 — 18, 28.02.1961.

2. Дахнов В. H. Промысловая геофизика.

Гостоптехиздат, 1959, с. 527 — 529.