автореферат диссертации по приборостроению, метрологии и информационно-измерительным приборам и системам, 05.11.13, диссертация на тему:Система контроля нефтяной скважины на базе импульсного нейтронного зондирования

Введение.

I. Анализ современного состояния аппаратуры нейтронного каротажа.

1.1. Принципы нейтронного исследования скважин.

1.2. Обзор задач по оптимизации аппаратуры нейтронного каротажа.

1.3. Выводы.

II. Обзор информационности и точности импульсной нейтронометрии.

2.1. Информативность нейтронометрии при исследовании неоднородных пластов.

2.2. Выводы.

III. Кислородный нейтрон-активационный метод (КНАМ).

3.1. Введение.

3.2. Описание известной «равновесной» методики КНАМ.

3.3. Формулировка требовании к аппаратуре для метода КНАМ.

3.4. Предлагаемая «импульсная» модификации КНАМ.

3.5. Физическое моделирование задачи определения заколонных перетоков воды по данным КНАМ.

3.6. Изучение нестабильности во времени параметра г\ в пресной воде.

3.7. Примеры исследований скважин с использованием разработанного макета.

3.8. Экспериментальное исследование «импульсной» методики КНАМ.

3.9. Выводы.

IV. Математическая модель временного распределения нейтронов и алгоритм её обработки.

4.1. Описание математической модели пласта.

4.1.2. Предлагаемый алгоритм обработки временных спадов.

4.2. Разработка программы обработки данных импульсного нейтронного каротажа.

4.2.1. Основные положения.

4.2.2. Формат файла исходных данных.

4.2.3. Входные данные при запуске программы.

4.2.4. Алгоритм математической обработки.

4.2.5. Вывод результатов обработки.

4.3. Описание созданного программного комплекса.

4.3.1. Введение.

4.3.2. Основные требования к первичному материалу и формату файлов исходных данных.

4.3.3. Интерфейс программы.

4.3.4. Последовательность обработки данных.

4.3.5. Ввод стартовых значений перед началом математической обработки.

4.3.6. Вывод результатов обработки.

4.3.7. Правила пользования программным комплексом.

4.4. Проверка точности определения 1 для среды с известными параметрами.

4.5. Выводы.

V. Техническая реализация и внедрение.

5.1. Введение.

5.2.Режимы работы АИНКЗб-ЗЦ.

5.3. Описание работы микроконтроллера.

5.3.1 Задачи, выполняемые микроконтроллером.

5.3.2 Описание программы.

5.3.3 Алгоритм работы микроконтроллера.

5.4. Технические параметры геофизического комплекса АИНКЗб-ЗЦ.

5.5. Описание структурной схемы скважинного прибора.

5.6. Состав и соединение комплекса АИНКЗб-ЗЦ.

5.6.1 Состав АИНКЗб-ЗЦ.

5.6.2 Схема соединения комплекса на основе АИНКЗб-ЗЦ.

5.6.3 Описание блока согласования (панели).

5.7. Результаты опробования импульсного генератора нейтронов.

5.8. Выводы.

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

Методы радиоактивного каротажа*, применяющиеся при контроле за разработкой нефтяных месторождений, можно разделить на пять групп: каротаж - процесс обследования нефтегазовой скважины геофизическими приборами). 1 .Гамма-каротаж.

2.Нейтронный каротаж со стационарным (ампульным) источником нейтронов, включающий нейтронный гамма-каротаж (НТК), нейтронный каротаж по тепловым нейтронам (НКТ).

3 .Нейтронный каротаж с импульсным источником нейтронов (ИНК), включающий импульсный нейтронный каротаж по тепловым нейтронам (ИННК) и импульсный нейтронный гамма-каротаж (ИНГК);

4.Каротаж по наведённой активности кислорода, или кислородный нейтрон -активационный метод (КНАМ).

5.Гамма-гамма-каротаж (ГГК), включающий ГГК по цементному кольцу и колонне — гамма-гамма цементометрия и толщинометрия или дефектометрия, и ГГК скважинной жидкости с установками гамма просвечивания и по рассеянному гамма-излучению.

В данной работе рассматривается метод импульсного нейтронного каротажа (ИНК) и кислородный нейтрон - активационный метод (КНАМ). Первое опробование метода ИНК, проведенное в 1959 г. на Туймазинском и Ромашкинском месторождениях нефти, показало его эффективность при определении местоположения водонефтяного контакта и выделения обводненных интервалов и пластов. В дальнейшем была подтверждена высокая разрешающая способность ИНК, позволяющая решать задачу определения интервалов обводнения пластов в обсаженных скважинах при значительно меньшей минерализации пластовых вод и меньшей пористости коллектора, чем это необходимо для стационарных нейтронных методов. Эти достоинства импульсного нейтронного каротажа обусловили его быстрое и широкое распространение по всей стране. Важным преимуществом метода ИНК является возможность количественного определения среднего времени жизни тепловых нейтронов в горных породах х, которое количественно может быть связано с нефтенасыщенностью пород.

Однако в процессе многолетней эксплуатации выявились недостатки аппаратуры ИНГК. Такие как большой диаметр спускаемого скважинного прибора, малая скорость регистрации, низкая надежность, сложность в эксплуатации и др. В настоящее время возникла необходимость создания геофизической аппаратуры ИНГК на основе современной элементной базы с цифровой регистрацией получаемой информации. Аппаратурный комплекс должен удовлетворять всем метрологическим требованиям, предъявляемым к аппаратуре данного типа. Кроме того, разрабатываемый геофизический комплекс должен иметь возможность для проведения исследований скважин не только методом ИНГК, но и методами ГК и КНАМ. Необходимо, чтобы скважинный прибор имел сверхмалый диаметр внешнего, защитного кожуха (до 36 мм), для проведения исследований скважин через насосно-компрессорные трубы (НКТ), т.к. их диаметр 45^-60мм. Точность определения нейтронного параметра исследуемого пласта (х-времени жизни тепловых нейтронов) должна быть не хуже 2%.

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

Для достижения поставленной цели, в диссертации решаются следующие задачи:

1. Поиск основных оптимальных аппаратурных и методических параметров аппаратуры нейтронного каротажа (НК).

2. Исследование влияния различных параметров пласта на информативность метода импульсного нейтронного каротажа (ИНК).

3. Разработка математической модели и алгоритма обработки данных, получаемых методом ИНК.

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

Научная новизна исследований представлена следующими результатами:

- Предложены и разработаны методика обработки и алгоритм вычисления т-пласта (время жизни тепловых нейтронов) по данным ИНГК, позволяющие количественно определять нефтенасыщенность пласта.

- Разработана новая модификация КНАМ, позволяющая в 2 раза повысить точность определения скорости жидкости в стальной колонне скважины и за ней и обладающая расширенным диапазоном измеряемых

•з о потоков (от 3 м /сут до более 100 м /сут для 128мм колонны).

- Предложен и создан трёхзондовый вариант скважинного прибора, реализующий методы ИНГК, КНАМ и ГК, с улучшенными эксплуатационными характеристиками.

Практическая ценность работы состоит в том, что:

1. Разработан и исследован новый аппаратно-програмный комплекс АИНКЗб-ЗЦ с улучшенными точностными и эксплуатационными характеристиками, позволяющий проводить исследования нефтяных и газовых скважин методами ИНГК, ГК и КНАМ, состоящий из:

- Скважинного прибора сверх малого диаметра (36 мм).

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

- Программного обеспечения для управления скважинным прибором и обработки данных, полученных комплексом АИНКЗб-ЗЦ, с использованием современных вычислительных средств.

2. Практически проверена новая модификация КНАМ.

3. Результаты диссертационной работы внедрены в АО «Татнефтегеофизика». В настоящее время разработанный комплекс эксплуатируется в ОМП-ЗО (Бугульма), АУГР (Альметьевск), Тимано-Печёрской экспедиции (республика Коми), Нефтегорская геофизическая экспедиция (Самарская область), Узеньгеофизика (республика Казахстан), ЧЭГИС (Пермская область).

Основные положения, выносимые на защиту;

1. Результаты исследования влияния различных параметров изучаемого пласта (пористость и минерализация) на информативность и достоверность получаемых результатов.

2. Методика обработки и алгоритм вычисления т-пласта (время жизни тепловых нейтронов) по данным ИНГК. Предложен алгоритм количественного вычисления т-пласта, основанный на приведении временного распределения тепловых нейтронов и соответствующих им у-квантов к экспонентной модели распределения.

3. Разработанный и запущенный в мало серийное производство новый цифровой аппаратно-программный комплекс трёх зондового

11 импульсного нейтрон-гамма исследования нефтяных и газовых скважин сверхмалого диаметра 36мм (АИНКЗб-ЗЦ).

4. Разработанная новая «импульсная» модификация КНАМ, позволяющая вычислять скорость жидкости в стальной колонне скважины и за колонной. Данная модификация КНАМ основана на меточном способе измерения скорости, в котором радиоактивные метки ставятся путём активации ядер кислорода.

Автор выражает глубокую благодарность за всестороннюю поддержку начальнику Научно-технического управления АО «Татнефтегеофизика» Юсупову Р.И., и главному инженеру НТУ Воронкову Л.Н. Получение приведённых в работе результатов было бы не возможно без большой квалифицированной помощи в различное время в выполнении исследований, опробовании и внедрении технологий сотрудников Опытно-методической партии №30 Баженова В.В., Лифантьева В.А., Ведерникова И.Р., Наплавкова Е.И., Хуснуллина P.A., Петрухина A.C., Труханова О.В., Батырева В.А., Крюкова В.П. и др.

В диссертации представлены результаты исследований, выполненных лично автором в Научно-техническом управлении АО «Татнефтегеофизика» в период с 1997 по 2002 гг.

I. Анализ современного состояния аппаратуры нейтронного каротажа.

Общие выводы.

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

В заключение диссертации можно сформулировать следующие выводы:

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

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

3. Импульсный нейтронный каротаж (ИНК) обладает рядом преимуществ по сравнению с непрерывным. Такими как безопасность и технологичность исследования, а так же возможность количественного определения параметров исследуемого пласта.

4. Найдено, что минимальная погрешность определения коэффициента диффузии нейтронов в пласте гшп(Овпл) методом двухзондового цифрового ИНК достигается при расстоянии детекторов от излучателя Ъ\= О 30 см, х2= 60 80 см и временных задержках 14«700-^-900 мкс.

5. Определены два основных условия информативности импульсного нейтронного каротажа: а) Ха ш^Ха с, где Ха пл -макроскопическое сечение захвата нейтронов в пласте, 2а с - макроскопическое сечение захвата нейтронов в скважине, б) 1 > 1а с, где ^временная задержка после вспышки нейтронов, ^ с -время диффузии нейтронов в зоне скважины. Показано, что информативность метода ИНГК выше, чем метода ИННК.

6. Из различных источников найдены зоны информативности (применимости) метода ИНГК в зависимости от пористости и минерализации пластовых вод.

7. Информативность и помехоустойчивость декремента затухания X во всех случаях заметно выше, чем скорость счёта Г

8. Для получения надёжной диаграммы декремента затухания X от глубины, необходима непрерывная цифровая запись декремента затухания с малым шагом квантования (например, на магнитный носитель), а измерения Х(Н) на точках или аналоговая регистрация X непригодны.

9. В результате исследования известной «равновесной» модификации КНАМ подтверждена возможность определения скорости жидкости в колонне и за колонной с помощью разработанного макета, который обладает повышенной чувствительностью (1,5 раза) к изменению скорости жидкости, по сравнению с аналогичными приборами. Полученная точность определения расхода составила 10%.

10. Макет был испытан как при каротаже со скоростью движения каротажного зонда 30+120 м/ч, так и при измерениях на точках. Получена чувствительность к расходу воды в колонне (5 дюймов) - 10 м3/сут и к расходу за колонной 0,7 м /сут при точности 10%. Выявлены недостатки «равновесной» модификации КНАМ: а) небольшая верхняя граница о измерений расхода, составляющая 70 м /сутки; б) искаженные оценки при наличии встречных потоков (например, поток вниз в колонне нагнетательной скважины и заколонный переток вверх).

11 .Предложена новая «импульсная» модификация КНАМ. В результате доказана возможность использования её для определения расхода жидкости в нагнетательных и добывающих скважинах при расходах от

3 3

5 м /сут для 5" (128мм) колонны до более чем 100 м /сут. Точность определения расхода составила 5%. Данный метод может быть применён для определения заколонных перетоков в нагнетательных скважинах.

12.При малых скоростях за счёт турбулентности жидкости замечен так называемый эффект «размывания» радиоактивной метки, который определяет нижнюю границу измеряемого потока.

13.Создан алгоритм количественного вычисления времени жизни тепловых нейтронов (Хпл), который позволяет существенно повысить достоверность измерения нефтенасыщенности и пористости исследуемого пласта.

14.На основе предложенного алгоритма создан программный комплекс по обработке данных ИНГК, ГК и КНАМ. Результаты работы данной программы представлены в универсальном формате (ЬА8-формат), который поддерживается большинством специализированных программ для построения геофизических диаграмм. Программный комплекс вошёл в состав геофизического комплекса аппаратуры импульсного нейтронного каротажа (АИНКЗб-ЗЦ).

15.С помощью разработанной программы проведена обработка данных, полученных в эталонной среде, в результате которой получена точность определения (тш) порядка 1%.

16.Разработан геофизический комплекс импульсного нейтрон-гамма каротажа (АИНКЗб-ЗЦ) нового поколения, обладающего повышенными точностными и информативными характеристиками, который имеет возможность работать не только в режиме ИНГК, но и в режимах ГК и КНАМ.

17.Благодаря разработанной нами трёх зондовой конструкции скважинного прибора (два зонда снизу, относительно нейтронной трубки, и один зонд сверху) достигнуты следующие результаты: a. Повышена надёжность определения т-пласта. b. Производится дублирование показаний зонда «60», за одно исследование. c. Возможно одновременное определение скорости встречных потоков жидкости (в скважине и в цементном кольце).

18.Проведены испытания АИНКЗб-ЗЦ, создана конструкторская документация и комплекс запущен в мало серийное производство в АО «Татнефтегеофизика». В настоящее время разработанный комплекс эксплуатируется в ОМП-ЗО (Бугульма), АУГР (Альметьевск), Тимано-Печёрской экспедиции (республика Коми), Узеньгеофизика (республика Казахстан).

19.Опробование аппаратуры и методики АИНКЗб-ЗЦ осуществлено в более 170 скважинах различной глубины с различным литологическим составом, в различных регионах России, Татарстана, Республики Коми и Казахстана.

124

Отмечена высокая эффективность использования данных ИНГК и КНАМ в нефтяных скважинах в комплексе с другими методами для решения выше приведённых задач.

1. Алексеев Ф.А., Головацкая И.В., Гулин Ю.А. и др. Ядерная геофизика при исследовании нефтяных месторождений. М., «Недра», 1978. 359 с.

2. Воронков JI.H., М.Х. Хуснуллин, Ю.М. Малыгин и др.Совершенствование и внедрение ядерно-геофизических и гидродинамических методов и аппаратуры исследований скважин. Отчёт по теме 30/85, фонд треста «Татнефтегеофизика», Бугульма, 1986, 208 с.

3. Воронков JI.H., М.Х. Хуснуллин, Ю.М. Малыгин и др. Разработка, своение и внедрение новой техники методики радиоактивных исследований скважин.Отчёт по теме 30/73, фонд треста «Татнефтегеофизика», Бугульма, 1973. 205 с.

4. Гермогенова Т.А. Развитие численных методов решения кинетического уравнения в задачах физики засщиты от излучений. Tp.II Всесоюзн. конф.по защите от ионизирующих излучений ядерно-технических установок (Москва, дек. 1978). М.,изд. МИФИ, 1978.

5. Гост 8.207-76.Прямые измерения с многократными наблюдениями. «Стандарт», М., 1979.10с.

6. Гутников B.C. Интегральная электроника в измерительных устройствах. -2-е изд., перераб. и доп.- Сан.-Петр.: Энергоатомиздат. 1998.-304 с. ил.

7. Дворянин В.В. Радиотехнические измерения. «Энергия», М.,1990.359с.

8. Зверев Г.Н. Основы теории машинной интерпретации и анализа эффективности ГИС. Докт.дис.М.,МИНХиГП, 1989.

9. Иванкин В.П. К методике исследования скважин нейтронным гамма-методом Уч.зап.Сарат.ГУ.Вып.геол.,1970, т.74, с.311-317.

10. Изучение элементного состава горных пород и насыщающих их флюидов нейтронными методами ядерной геофизики/ К.И.Якубсон, В .В .Стрельченко, В.В.Муравьёв и др. Тр. МИНХиГП, 1974, вып. III, с. 136-204.

11. Казарновский М.В. Теория нестационарной термализации нейтронов.- В кн.:Ядерные реакции и взаимодействие нейтронов свеществом. М.,Наука, 1972, с. 3-98. (Труды ФИАН СССР, т.63)

12. Кантор С.А. , Кожевников Д.А. , Поляченко A.JI. , Шимелевич Ю.С. Теория нейтронных методов исследования скважин. М., Недра, 1985. -224с.

13. Кантор С.А Теоретические основы нейтронных методов исследования горных пород, пересечённых скважиной. Док. Дис. М.,ВНИИЯГ, 1980г.

14. Киргизов Д.И., Ильин Г.И., Оптимизация режимов измерений и аппаратуры геофизических исследований в комплексе импульсного нейтронного каротажа, //Всеросийская студенческая научная конференция, Тез. докл, Самара, СГАУ им. Королёва С.П, 1999.

15. Киргизов Д.И., Ильин Г.И., "Новый аппаратно-программный комплекс скважинного импульсного генератора нейтронов АИНКЗб-ЗЦ и егоприменение в нефтяной геологии"// журнал «Электронное приборостроение» №26, 2002г., Казань, КГТУ им. Туполева А.Н.

16. Кулинкович А.Е. , Басин Я.Н. , Применение теории информации к вопросам радиоактивного каротажа. В кн.: Ядерная геофизика, М.,Недра, 1988, с. 81-98.

17. Клинт Хинт. С (серия «Без проблем!»)/ Пер. с англ. М.: Восточная Книжная Компания, 1997. - 448 е.: ил.

18. Лебедев В.Е. Результаты расчёта комплексом численных методов процесса нестационарной диффузии нейтронов в многопластовых средах со скважиной. Кан. Дис. М., ВНИИЯГ,1968.

19. Микроконтроллеры.Выпуск 1. Перевод с англиского -Б.Я. Прокопенко, М.ДОДЭКА, 1998, 384с.

20. Методы Монте-Карло в физике и геофизике. Уфа, изд. БашГУ,1983.

21. Новые идеи поиска, разведки иразработки нефтяных месторождений: Труды научно практической конференции VII международной выставки «Нефть, газ -2000» (Казань, 5-7 сентября 2000 года). — Казань: Экоцентр, 2000.-468с.

22. Немнюгин С.А. Turbo Pascal. СПб: Издательство «Питер»,2000ю -496с. :ил.

23. Поляченко AJI. Метод теоретико-экспериментальных палеток радиоактивного каротажа (ТЭП PK). В Кн.: Ядерная геофизика при подсчёте запасов нефти и газа. М., ВНИИЯГГ, 1980, с. 51-56.

24. Поляченко A.JI. , Теория обратных задач нефтегазовой скважинной ядерной геофизики. Докт. Дис. М., ВНИИЯГГ , 1982.

25. Попов Н.В. О применении методов неупругого рассеяния нейтронов и наведённой радиоактивности для выделения углесодержащих пластов и оценки их качества. В Кн.:Портативные генераторы нейтронов в ядерной геофизике. М., Госатомиздат, 1962, с 127-144.

26. Путкадзе JI.A. Комплексное использование материалов электрического и импульсного нейтронного каротажа для оценки изменения нефтенасыщенности коллекторов разрабатываемых залежей. Баку, издАзВНИИГеофизика, 1978.

27. Резванов Р.А Нейтронные методы исследования нефтяных и газовых скважин (при разведке и разработке залежей, приуроченных к пластам сложного состава и строения). Док. Дис. М. МИНХиГП, 1992.

28. Руководство по применению промыслово-геофизических методов для контроля за разработкой нефтяных месторождений. М., Недра, 1978. 256с.

29. Скважинная ядерная геофизика. Справочник геофизика ( под ред. В.М. Запорожца). М., Недра, 1978.

30. Степанов A.B. К теории переноса нейтронов в неоднородных средах. Тр. ФИАН СССР, 1969, т.44,с. 204-257.

31. Султанов С.А. Контроль за обводнением нефтяных пластов. М., Недра, 1984.

32. Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач. М., Наука, 1979.

33. Фаронов В.В. Турбо Паскаль 7.0. Начальный курс. Учебное пособие. Издание 7-е, переработанное.-М.:Нолидж,2ООО.-576 е., ил.(ISBN 5-89259054-9)

34. Филиппов Е.М. Ядерная геофизика. Новосибирск, Наука, 1973.

35. Физические основы импульсных нейтронных методов исследования скважин/ Ю. С. Шимелевич, С. А. Кантор, А. С. Школьников и др. М., Недра, 1976.

36. Франк Г. М. Импульсный метод исследования свойств медленных нейтронов. Тр. ФИАН СССР, 1962, т.14 с.117-146.

37. Хуснулин М.А. Применение гамма-метода для определения заводнённых пластов. Геология нефти и газа, №12, 1973.

38. Хоровиц П., Хилл У., Искусство схемотехники. М., «Мир», 1998г.

39. Цейтлин В. Г. Разработка теоретических основ количественной интерпретации данных ИНК нефтяных скважин. Канд. Дис. М., ВНИИЯГГ, 1985.

40. Шапошникова Т. А. Создание конечно-разностного аппарата для расчёта полей излучения при нейтронометрии скважин. Кан. Дис. М., ВНИИЯГ, 1982.

41. Шимелевич Ю.С. Физические основы нейтронных методов исследования скважин с использованием временного анализа излучений. Док.дис. М.,ВНИИЯГГ,1968.

42. Широков Ю.М., Юдин Ю.П. Ядерная физика. М., Наука, 1982.51 .Юров В. Assembler. СПб. Литер, 2001. - 624 е.: ил.

43. Ядерная геофизика при исследовании нефтяных месторождений/ Ф.А. Алексеев, И.В. Головацкая, Ю.А. Гулин и др. М., Недра, 1988.

44. Allen L.S., Mills W. R., Caldwell R. L. The effects of fluid invasion in pulsed neitron logging. Geophysocs, 1995, vol. 30, №3, pp 115-121.

45. Karus E.V., Shimelevich Yu. S. Nuclean Geophysics in Prospecting, Exploration and Developmen in Oil and Gas Fields. Int. J. Appl. Isot. 1983, vol. 34, №1, pp. 95-117.

46. Mills W.R., Allen L.S. CUNLAR (Computer usage neutron log analyzer program, 655). -Nucl. Sci. Eng., 1986, vol. 21, № 3, p 441.

47. Pontecorvo B. Neutron well logging.- Oil and Gas J.,1981, vol.40, No 18, pp. 32-33.

48. Tittle C.W., Allen L.S. Theory of neutron logging. Geophysics, 1991, vol. 26, No l,pp. 27-39; 1996, vol. 31, No. l,p.214.

49. Transaction of 96 Internetional Symposium on Well Logging Techniquest for Oilfield Development under Waterflooding. Copyrite 1996 by Petroleum Industry Press.

50. Wyman R.E. How should we measure residual-oil saturation? Bull. Canadian Petr. Geol., 1987, vol 25, №2, pp. 234 - 270.

51. Youmans A.H., Bishop W., Wichman P. A. Application of the neutron lifetime log in new wells. Log Analyst, 1981, vol. 12, № 2, pp 3-13.