Развитие компьютерных технологий статистического моделирования для решения задач стационарного и импульсного нейтронного каротажа скважин тема диссертации и автореферата по ВАК РФ 25.00.10, кандидат технических наук Цымбал, Дмитрий Васильевич

• Цымбал, Дмитрий Васильевич
• кандидат технических наук
• 2005

Оглавление диссертации кандидат технических наук Цымбал, Дмитрий Васильевич

Глава 1. Развитие технологий стандартного ННК (CN 2435/38 Western Atlas) с использованием компьютерного моделирования

1.1 Краткая история вопроса

1.2 Описание прибора CN

1.3 Компьютерная модель

1.4 Тестирование модели

1.5 Краткое описание процесса калибровки

1.6 Обоснование двухшаговых систем расчетов. Сходимости статистических характеристик

1.7 Краткое описание системы статистических тестов

1.8 Создание прецизионных компьютерных моделей (КоМоД2) для приборов НК

1.9 Обоснование компьютерных моделей ННК на расчетах калибратора

1.10 Результаты расчетов и их обсуждение

1.11 Сопоставления расчетов, выполненных по программам MCNP4C и MCNP5. Влияние уточненных ядерных данных

1.12 Результаты расчетов палеток для обсаженной скважины

1.13 Стандартный прибор ННК в условиях скважины с двухколонной обсадкой

1.14 Вывод к главе

Глава 2 Компьютерное моделирование при проектировании новых приборов ННК для каротажа горизонтальных скважин (вторых стволов) и систем интерпретации.

2.1 Краткая история вопроса

2.2 Постановка задачи

2.3 Особенности компьютерных моделей

2.4 Выводы к главе

Глава 3. Компьютерное моделирование задач импульсного нейтронного каротажа. Сравнительный анализ эффективности применения ИННК и ИНГК на единой аппаратурной основе прибора РБК-ЮО

3.1 Краткая история вопроса

3.2 Сопоставительный анализ применения ИННК и ИНГК в Западной Сибири (модельные и скважинные измерения)

3.3 Краткое описание аппаратуры РБК-ЮО

3.4 Метрология ИНК

3.5 Структура сигнала ИНК

3.6 Радиальное зондирование

3.7 Термализация в задачах ИННК/ИНГК

3.8 Задача ИНК о двухкомпонентном разложении сигнала

3.9 Выводы к главе

Введение диссертации (часть автореферата) на тему «Развитие компьютерных технологий статистического моделирования для решения задач стационарного и импульсного нейтронного каротажа скважин»

Актуальность работы

Двухзондовые отечественные приборы стационарного нейтронного каротажа (СРК, РКС) появились 30 (примерно) лет назад на смену однозондовым приборам типа СП-62, ДРСТ. Тогда же была разработана идеология создания системы палеток для новой аппаратуры на единой расчетной основе метода Монте-Карло (Еникеева Ф.Х, Журавлев Б.К. Велижанин В.А. Стариков В.Н. и др). Расчетная палеточная база охватывала следующие факторы: открытый ствол-отклонение прибора от стенки, переменный диаметр, переменная ^ толщина и плотность глинистой корки, переменная литология, переменная минерализация в скважине и пласте, переменная температура и давление (8-9 факторов, которые распадались на 2 группы). В обсаженной скважине учитывали переменный диаметр, толщину цементного кольца, толщину и диаметр обсадки. Объединение всех поправок в единую расчетную базу (с быстрой выборкой) взамен разрозненной системы различных поправок было огромным шагом вперед и обеспечило успех расчетных Монте-Карловских методов в геофизике.

За рубежом эта технология появилась чуть позднее (1Ви11ег, 1.ВпЙ1те1з1ег) однако сразу стала обязательной для всех ведущих геофизических компаний. Расчетная палеточная база облегчила внедрение компьютерной интерпретации нейтронного каротажа. Выходным параметром стала нейтронная пористость, исправленная за все мешающие факторы. Вначале эта технология появилась на Западе, затем стала внедряться в России (хотя и с некоторым трудом). Ш

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

В последние годы бурение горизонтальных скважин (вторых стволов) во всем мире и в России развивается исключительно быстро. Для каротажа вторых стволов в России применяют разные технологии, в частности, малогабаритную автономную аппаратуру на буровом инструменте. В отличие от технологии ГОРИЗОНТ, которая не допускала промывки горизонтального ствола в процессе каротажа, в разрабатываемой автономной аппаратуре ПетроАльянса предполагается использование проточной схемы. В процессе каротажа производится промывка, чтобы избежать заклинивания прибора и не иметь ограничений по времени работы.

Импульсный нейтронный каротаж (ИНК) был разработан в начале 60гг. Первые успешные испытания были проведены в 1961г. Метод был практически реализован и внедрен усилиями очень многих специалистов, среди которых отметим Алексеева Ф.А., Ерозолимского Б.Г, Школьникова A.C., Беспалова Д.Ф., Шимелевича Ю.С., Поляченко А.Л., Воронкова JI.H. L.Allen, B.Mills, R.Caldwell и многих других. Метод широко применяется во всем мире, однако новые задачи требуют развития и этого метода.

Библиография работ по нейтронному каротажу огромна. Главный вклад в создание современной картины метода внесли Дядькин И.Г., Гулин Ю.А. Кантор С.А. Поляченко A.JI., Кожевников Д.А. Резванов P.A., Еникеева Ф.Х., Журавлев Б.К. Велижанин В.А. Хаматдинов Р.Т. Стариков В.Н. Лухминский Б.Е. Давыдов Ю.Б., Крапивский Е.И., J. Schweitzer, D. Oliver, D.Ellis, R.Gardner, L.Allen, B.Mills, R.Caldwell и очень многие другие.

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

Щ необходим инструмент для быстрого решения этих задач методами математического моделирования (компьютерного эксперимента). Этому посвящена работа.

Цели работы и решаемые задачи:

1. Создание иерархии (семейства) компьютерных моделей на основе метода Монте-Карло для расчета актуальных задач стационарного и импульсного каротажа скважин.

2. Тестирование моделей (Validation & Verification) на стандартных калибраторах и модельных полигонах для получения т погрешностей расчета, адекватных погрешностям калибровки аппаратуры (1-2%).

3. Расширение стандартного набора палеток CN Western Atlas путем расчета палеток и поправок на нетрадиционную литологию коллектора (полимиктовый песчаник, наличие нейтронных поглотителей, нетрадиционный гранитоидный коллектор), на измерения в обсаженной скважине (гель-цементом, двойная колонна и т.д.)

4. Развитие элементов компьютерной метрологии путем сопоставления калибровочных кривых для различных модельных полигонов (например, Хьюстонского и Раменского), а также путем сопоставления чувствительности к пористости различных приборов (CN Western Atlas, автономный прибор ЗННК, двухзондовый С\0-генератор МАРКА).

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

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

Научная новизна:

1. Созданы элементы компьютерной метрологии ННК для серийной аппаратуры CN Western Atlas, сопоставлены калибровочные зависимости для различных модельных полигонов и расширен основной набор палеток за счет нетрадиционных коллекторов.

2. На основе компьютерного моделирования выбраны параметры многозондовой автономной аппаратурой ННК и созданы новые приемы 2D и 3D интерпретации.

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

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

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

2. На основе математического моделирования созданы новые приемы 2D и 3D интерпретации для автономной нейтронной аппаратуры каротажа горизонтальных скважин и оптимизированы ее параметры.

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

Апробация работы

Основные положения диссертации докладывались на V, VI, УП Международных конференциях «Новые идеи в науках о Земле» (2001, 2003, 2005 годы), на Молодежной секции научно-практической конференции «ГЕОМОДЕЛЬ-2002» (2002 год) и на Научно-практической конференции «Ядерная геофизика - 2004» (2004 год, Санкт-Петербург).

Практическая реализация и внедрение результатов работы

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

Результаты, полученные автором, использованы для интерпретации и внесения поправок в стандартный каротаж, выполняемый с помощью аппаратуры CN в ПетроАльянсе. Результаты работы использованы также специалистами ЦГЭ Минэнерго при ревизии запасов месторождения Белый Тигр (СРВ).

Структура и объем

Диссертация состоит из 110 страниц текста с приложениями, 49 рисунка, 3 таблиц и 87 позиций библиографии.

Благодарности

Автор благодарит специалистов ПетроАльянса В.П. Стенина, В.В. Метлина, А.А. Чередниченко, В.В. Костылева, М.А.Ракитина, а также специалистов Когалымского и Усинского подразделений за поддержку; специалистов ГЕОКОНА, ВНИИА, сотрудника ОИЯИ Рогова А.Д, сотрудников кафедры ЯРМ и ГИ и кафедры математического моделирования МГТРУ за полезные советы, своего руководителя Б.Е.Лухминского.

Заключение диссертации по теме «Геофизика, геофизические методы поисков полезных ископаемых», Цымбал, Дмитрий Васильевич

3.9 Выводы к главе 3. Реализована компьютерная модель, адекватно описывающая параметры прибора РОК-100.

2. Измерения ИНГК обеспечивают большую скорость счета (в 510 раз) по сравнению с ИННК, что приводит к меньшей погрешности измерений и возможности проводить каротаж с большей скоростью.

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

3. Интервал между импульсами, равный 1000 мкс, достаточен для раздельного определения СИГМА-скважины и СИГМА-пласта с погрешностью не хуже 2%. Результаты компьютерного моделирования показывают, что фон предшествующего импульса не влияет на определение СИГМА-пласта, хотя и несколько влияет на определение СИГМА-скважины.

4. Режим ИННК из-за меньшей скорости счета и особенностей формы нестационарного сигнала позволяет определять СИГМА-пласта с погрешностью 3% и выше. Режим ИННК не требует учета фона, и режим генерации сигнала может быть непрерывным (без пауз для определения фона).

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

6. Расчеты задачи о многозондовом зондировании ИННК и ИНГК показали перспективы этого метода. В частности, при наличии минерализованной зоны проникновения радиусом не более 25см возможна качественная оценка СИГМА-пласта.

7. По заказу Петроальянса составлена и опробована программа СИГМАЕ, использующая новые алгоритмические подходы при обработке экспериментальных и расчетных материалов ИНК в двух компонентном разложении.

Заключение.

В работе получены следующие основные результаты:

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

2. Выполнен расчет поправок не нестандартную литологию коллектора, отсутствующих в стандартном наборе палеток CN Western Atlas (полимиктовый песчаник, наличие нейтронных поглотителей, нетрадиционный гранитоидный коллектор, обсаженная скважина с гель-цементом, двухколонная конструкция обсадки скважины и т.д.).

3. Развиты элементы компьютерной метрологии путем сопоставления калибровочных кривых для различных модельных полигонов (например, Хьюстонского и Раменского), а также путем сопоставления чувствительности к пористости различных приборов (CN Western Atlas, автономный прибор ЗННК, двухзондовый С\0-генератор МАРКА).

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

5. Выполнена компьютерная оптимизация элементов новой многозондовой автономной аппаратуры ЗННК (выбор измерительных параметров, оценка помехоустойчивости, чувствительности) и предложены новые приемы 2D и 3D интерпретации ЗННК, обладающие более высокой чувствительностью к пористости и помехоустойчивостью относительно мешающих факторов.

6. Сопоставлены сравнительные достоинства методов ИННК и ИНГК на базе аппаратуры PDK-100. Показана практическая эффективность ИНГК и одновременно выявлена более высокая чувствительность ИННК к пористости. На основе компьютерного анализа перспектив многозондового импульсного зондирования показано, что методом ИНГК возможна оценка СИГМА-пласта (по совокупности зондовых показаний) при проникновении минерализованного бурового раствора в коллектор на глубину не более 25см.

Список литературы диссертационного исследования кандидат технических наук Цымбал, Дмитрий Васильевич, 2005 год

1. Computational Science and Engineering, sept-oct. 2004 Validation and Verification

2. Лухминский Б.Е., Тепляков A.B., Рогов А.Д. Современный этап компьютерного моделирования Монте-Карло с использованием системы MCNP для модернизации ядерно-геофизических технологий исследования скважин, КАРОТАЖНИК, 93, 101

3. Лухминский Б.Е., Тепляков А.В. Применение методов теории возмущений при решении сложных задач ГИС методом Монте-Карло. КАРОТАЖНИК, 107, 70

4. Czubek J.A. Neutron Tool Calibration by a Scaling Procedure. Nucl. Geophys. 8, #3, 261-279, 1994

5. Butler J, McCracken A.K., Kemshell P.B. A Study of Uncertainties in the Prediction of Neutron Porosity Logging-tool Responses. Nucl. Geophys. 3,#3, 179-191, 1989

6. Proceedings of SPWLA-45, Netherlandjune 2004

7. Hemingway J., Plasek R., Grau J. Introduction of Enhanced Carbon-Oxygen Logging for Multiwell Reservoir Evaluation SPWLA-40, 1999 paper О

8. Schlumberger Log Interpretation chart, 2000

9. Денисик C.A., Лухминский Б.Е., Резванов P.A., Метод статистических испытаний в приложении к расчету распределения нейтронов в задачах нейтронного каротажа, Сб. Портативные генераторы нейтронов в ядерной геофизике, Атомиздат, 1962

10. Ю.Денисик С.А., Дядькин И.Г., Кантор С.А., Кожевников Д.А., Лухминский Б.Е., Поляченко A.JI., Резванов Р.А., Рудык Б.М., К теории нейтронного каротажа

11. М.А. Arnautova, Ya. Z. Kandiev, B.E. Lukhminsky, G.N. Malishkin, Monte Carlo Simulation in Nuclear Geophysics/ Intercomparison of the PRIZMA Monte Carlo Program and Benchmark Experiments, Nucl. Geoph, vol.7, No 3, pp 407-418, 1993

12. А.А. Соколов, Ю.М. Лоскутов, И.М. Тернов. Квантовая механика, М., 1962

13. А.М.Кольчужкин, В.В. Учайкин. Введение в теорию прохождения частиц через вещество, М.,1978

14. Б.Е.Лухминский. Вычислительный эксперимент Монте-Карло для решения задач ядерно-геофизического каротажа рудных скважин, диссертация на соискание ученой степени, Москва 1986

15. В.А.Арцибашев. Ядерно-Геофизическая разведка, Москва АТОМИЗДАТ 1980

16. В.Д. Аллен. Регистрация нейтронов, Государственное издательство литературы в области атомной науки и техники, Москва 1962

17. В.М. Добрынин, Б.Ю. Венделыптейн, Р.А. Резванов, А.Н. Африкян. Промысловая геофизика

18. В.Ф. Турчин. Медленные нейтроны, Государственное издательство литературы в области атомной науки и техники, Москва 1963

19. Г.И. Марчук, В.И. Лебедев. Численные методы в теории переноса нейтронов, Атомиздат, Москва 1971

20. Г.И. Марчук. Методы расчета ядерных реакторов, Государственное издательство литературы в области атомной науки и техники, Москва 1961

21. Д. Каханер, К. Моулер, С. Нэш. Численные методы и программное обеспечение, МИР, Москва, 19 98

22. Д.И. Дьяконов, Е.И. Леонтьев, Г.С. Кузнецов. Общий курс геофизических исследований скважин

23. Дж. Льюинс. Ценность. Сопряженная функция.

24. И.В. Савельев. Курс общей физики. Книга 5, М., «Физматлит», 1998

25. И.Г.Дядькин. Лекции о применении методов Монте-Карло в физике, 1967

26. И.М.Соболь. Метод Монте-Карло, М.,1985

27. И.М.Соболь. Численные методы Монте-Карло, М.,1973

28. К. Бекурц, К. Виртц. Нейтронная физика, Атомиздат, Москва, 1968

29. Каротажник, Выпуск №60, Тверь 1999

30. Л.Г. Петросян. Геофизические исследования в скважинах, крепленных трубами, при изучении разрезов нефтегазовых месторождений, Москва, Недра, 1977

31. Л.П. Абагян, Н.О. Базарянц, И.И. Бондаренко, М.Н. Николаев. Груповые константы для расчета ядерных реакторов. Москва. "Атомиздат". 1964

32. Методические рекомендации по проведению исследованний и интерпретации данных нейтронного каротажа с серийной аппаратурой РК (комплект палеток), НПО "Союзпромгеофизика", Всесоюзный НИИ ядерной геофизики и геохимии, Москва, 1979.

33. Методические рекомендации по проведению исследованний и интерпретации данных нейтронного каротажа с серийной аппаратурой РК, НПО "Союзпромгеофизика", Всесоюзный НИИ ядерной геофизики и геохимии, Москва, 1979.

34. Методические указания по проведению нейтронного и гамма-каротажа в нефтяных и газовых скважинах аппаратурой СРК и обработке результатов измерений, НПО "Союзпромгеофизика", Калинин, 1989.

35. Многогруповая библиотека ядерных констант для ядерно-геофизических расчетов. Версия Б-2. Москва. 1976

36. Многозондовый нейтронный каротаж с аппаратурой К-7. Методическая инструкция по проведению измерений и интерпретации данных (комплект палеток), НПО "Нефтегеофизика", Всесоюзный НИИ ядерной геофизики и геохимии, Москва, 1980.

37. Н.А.Власов. Нейтроны, М.,1971

38. Н.А.Власов. Нейтроны, Государственное Издательство Технико-Теоретической Литературы Москва 1955

39. Н.П.Бусленко, Д.И. Голенко, И.М.Соболь, В.Г. Срагович, Ю.А. Шрейдер. Метод статистических испытаний, М.,1962

40. И.Лейпунский Б.В.Новожилов,В.Н.Сахаров. Распространение гамма-квантов в веществе, М.,1960

41. П. Гриблер, Э. Хенли. Импульсный метод в нейтронной физике, Атомиздат, Москва 1969

42. Под редакцией И. Марчука. Метод Монте-Карло в проблеме переноса излучения, Атомиздат, Москва 1967

43. Р.А.Резванов. Радиоактивные и другие неэлектрические методы исследованния скважин, Москва. "Недра". 1982

44. С.М.Ермаков. Метод Монте-Карло и смежные вопросы.

45. С.М.Ермаков, Г.А.Михайлов. Курс статистического моделирования, М.,1976

46. С.С. Итенберг, Т.Д. Дахкильгов. Геофизические исследования в скважинах. Москва. "Недра". 1982

47. Ю.А. Цирлин, М.Е. Глобус, Е.П. Сысоева. Оптимизация детектирования гамма-излучения сцинтилляционными кристаллами, М., «Энергоатомиздат», 1991

48. Bakhtejarov S.D., Kamenshtchikov A.F., Lukhminsky В.Е. Monte Carlo and parallel computations on transputer nets

49. C. D. Zerby. A Monte Carlo Calculation of the Response of Gamma-Ray Scintillation Counters, Oak Ridge National Laboratory, Oak Ridge, Tennessee

50. D.E. Peplow, R.P. Gardner, K. Verghese. Sodium Iodide Detector Response Functions Using Simplified Monte Carlo Simulations and Principal Components, Nuclear Geophysics, vol.8, no.3, 1994

51. G.J. Roach, J.S. Watt, H.W. Zastawny, P.E. Hartley, W.K. El???. Multiphase Flowmeter for Oil, Water and Gas in Pipelines Based on Gamma-Ray Transmission Techniques, Nuclear Geophysics, vol.8, no.3, 1994

52. H. Zheng. How the Depth of Investigation Changes with Spacing: A Computer Model for Carbon/Oxygen Logging, Petrophysics, Vol.42, No.l, Jan.-Feb. 2001

53. I.M. Sobol. An efficient approach to multicriteria optimum design problems, Surveys on Mathematics for Industry 1 (1992): 259-281

54. I.M. Sobol, Yu.L. Levitan. A Pseudo-Random Number Generator for Personal Computers, An International Journal Computers and Mathematics with Applications 37 (1999): 33-40

55. Judith F. Briesmeister, Editor. MCNP™ A General Monte Carlo N-Particle Transport Code Version 4C, Manual, March, 2000

56. Randy A. Schwarz, Lee L. Carter, Wendi Brown. Particle Track Visualization Using the MCNP Visual Editor

57. Thomas E. Booth. A Sample Problem for Variance Reduction in MCNP, Los Alamos National Laboratory, LA-10363-MS, october 1985

58. W.E. Sharp, Carter Bays. A review of portable random number generators, Department of Geological Sciences and Department of Computer Science, University of South Carolina, Columbia, U.S.A., 1991

59. Ядерная геофизика, труды ВНИИЯГГ, 1963г

60. Ядерная геофизика, труды ВНИИЯГТ, т7, 1967г

61. Теленков В.М. в Сб Современная ядерная геофизика при поисках, разведке и разработке нефтегазовых месторождений Москва 2004 стр 402-426

62. Лухминский Б.Е. Обзор по зарубежным генераторам нейтронов. 1998-2002

63. А.А.Чередниченко, Б.Е.Лухминский. Первая редакция. Генераторы нейтронов нового поколения (краткий обзор).

64. Интерактивная ЗД визуализация при проектировании скважин

65. Лухминский Б.Е. Генераторы нейтронов (обзор). (90-95гг)1. Петроальянс, 2000г

66. А.А.Чередниченко, Б.Е.Лухминский. Первая редакция. Развитие метрологической системы для спектрометрических генераторов типа MSI-CO.

67. Расчетное сопровождение современного ИНК в ПетроАльянсе, петроальянс 2002

68. M.W.Mickael, W.A.Gilchrist, Jr, R.J.Mirzwinski, G.N.Salaita, R.T. Rajasingam., Interpretation of the response of a new through-tubing carbon\oxygen instrument using numerical modeling techniques, SPWLA-38, Jl-Jl 1,1997.

69. Morris M, Hemingway J. Continuous Oil, Gas and Water Holdup Using Pulsed-Neutron Spectroscopy Techniques, SPWLA-40, N1-N13

70. POLARIS рекламное описание.

71. D.C.Stromswold, W.R.Mills, R.D.Wilson, T.K.Cook. Formation porosity measurement using epithermal neutron lifetime. IEEE transactions on Nuclear Science, Vol. 36. No. 1, February 1989

72. Лухминский Б.Е. Обзор зарубежных публикаций по горизонтальному бурению и геофизическому сопровождению, 1994, 98, 2000гг

73. Лукьянов Э.Е. Состояние и перспективы развития геофизических исследований в горизонтальных скважинах. Тверь 1994