Бесплатный автореферат и диссертация по геологии на тему

Методы быстрого 3D-моделирования полей ядерной геофизики

ВАК РФ 04.00.22, Геофизика

Содержание диссертации, кандидата физико-математических наук, Кулешова, Людмила Борисовна

ВВЕДЕНИЕ

Глава 1. МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПОЛЕЙ ЯДЕРНОЙ ГЕОФИЗИКИ: СОСТОЯНИЕ И ПРИКЛАДНЫЕ ВОЗМОЖНОСТИ

1.1. Прикладные возможности и фундаментальные аспекты теории ядерных методов ГИС

1.2. Состояние и проблемы математического моделирования ядерно-геофизических полей

Глава 2. РАЗРАБОТКА БЫСТРЫХ КОМБИНИРОВАННЫХ

МЕТОДОВ И АЛГОРИТМОВ ЗБ-МОДЕЛИРОВАНИЯ ПОЛЕЙ ЯДЕРНОЙ ГЕОФИЗИКИ

2.1. Принципы нового комбинированного вычислительного аппарата ядерной геофизики ПОЛЕ

2.2. Физико-математические модели вычислительного аппарата ПОЛЕ: состав и свойства

2.3. Модели взаимодействия нейтронов и гамма-квантов, модели состава сред и ядерно-петрофизические модели

2.4. Модель переноса нейтронов и постановки краевых задач

2.5. Построение адаптивных сеток и оценки теории сеточных схем для краевых задач ГИС

2.6. Конечно-разностная аппроксимация стационарных и нестационарных трехмерных краевых задач теории ГИС в вертикальных скважинах

2.7. Конечно-разностная аппроксимация стационарных и нестационарных трехмерных краевых задач теории ГИС в горизонтальных и наклонных скважинах

2.8. Групповые экономичные конечно-разностные алгоритмы ЗБ-моделирования нейтронных полей в ядерной геофизике

2.9. Быстрые алгоритмы ЗБ-моделирования физических спектров и интегральных потоков гамма-квантов на основе комбинированной вычислительной схемы

2.10. Быстрые алгоритмы 3Б-моделирования аппаратурных спектров и интегральных потоков гамма-квантов на основе комбинированной вычислительной схемы

2.11. О возможности обобщения вычислительного аппарата ЗБ-моделирования ПОЛЕ на электрические и тепловые поля ГИС в вертикальных, горизонтальных и наклонных скважинах

Глава 3. ПАКЕТ ПРОГРАММ "ПОЛЕ" БЫСТРОГО

ЗБ-МОДЕЛИРОВАНИЯ ПОЛЕЙ И ПОКАЗАНИЙ МЕТОДОВ ЯДЕРНОЙ ГЕОФИЗИКИ

3.1. Принципы программирования и проектирования пакетов прикладных программ для моделирования геофизических полей и решения прямых задач ГИС.

3.2. Состав, структура и функционирование пакета ПОЛЕ

3.3. Апробация пакета ПОЛЕ точными решениями

3.4. Тестирование ядерного пакета ПОЛЕ данными Монте-Карло и физических экспериментов на моделях пластов.

3.5. Сводные научно-прикладные характеристики пакета ПОЛЕ

Глава 4. ОПРОБОВАНИЕ И ПРИМЕНЕНИЕ ПАКЕТА "ПОЛЕ" ДЛЯ РЕШЕНИЯ ЗАДАЧ ГЕОФИЗИЧЕСКОЙ НАУКИ, ПРОИЗВОДСТВА И НИОКР

4.1. Алгоритмы и пакет программ ДЕЛЬТА для оптимизации зондов, режимов измерений и разработки интерпретационного и метрологического обеспечения аппаратуры ядерной геофизики

4.2. Оптимизация зондов и создание систем палеток и поправок для малогабаритных трехзондовых приборов ННКнт и ННКт

4.3. Методическое обеспечение переинтерпретации данных при подсчете запасов в юрских нефтеносных отложениях Зап. Сибири

4.4. Оптимизация зондов и интерпретационное обеспечение аппаратуры углерод-кислородного С/О-каротажа

4.5. Исследование влияния различных типов утяжелителей бурового раствора на показания ядерного каротажа в разрезах баженовской свиты продуктивных отложений Зап. Сибири

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

4.7. Исследование аппаратом ПОЛЕ задач скважинной гамма-спектрометрии, недоступных физическому моделированию

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

Введение Диссертация по геологии, на тему "Методы быстрого 3D-моделирования полей ядерной геофизики"

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

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

Моделирование методом статистических испытаний ограничено из-за его медленной сходимости. Она делает крайне трудным достижение за разумное время нужной точности для наиболее перспективных «многоканальных» методов ядерной геофизики (ЯГ): импульсных, спектрометрических, многозондовых и особенно их комбинаций.

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

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

В итоге интерпретационно-методическое обеспечение многих приборов ЯГ зачастую является весьма неполным, недостаточно достоверным и хронически отстает от потребностей практики, а сроки разработки новых АМК обычно оказываются слишком большими.

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

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

Предпосылками постановки и успешного выполнения такой работы явились: а) большой опыт численного моделирования задач ЯГ и высокий уровень развития теории ядерного каротажа в России вообще; б) резкое увеличение производительности компьютеров и эффективности их системного обеспечения (ОС, компиляторов, ИСР); в) развитие экономичных методов линейной алгебры и теории сеточных схем.

Цель и задачи работы. Цель работы состоит в разработке математического и вычислительного аппарата ЗБ-моделирования ядерно-геофизических полей и показаний ГИС (пакет алгоритмов и программ ПОЛЕ), который удовлетворяет одновременно 4-м взаимно-противоречивым требованиям, выдвигаемым практикой:

- точности;

- быстроты;

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

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

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

1. Разработка набора физико-математических моделей методов ЯГ, алгоритмов построения адаптивных сеток и разностной аппроксимации краевых задач ядерных и других полей ГИС в сложных ЗБ-геометриях, включая горизонтальные и наклонные скважины; их реализация в программных модулях пакета ПОЛЕ.

2. Разработка методов и алгоритмов быстрого ЗВ-моделирования полей нейтронов и гамма-квантов в ЯГ на основе комбинированных вычислительных схем (КВС); их реализация в программах-"решателях" пакета ПОЛЕ.

3. Создание современной версии пакета прикладных программ ПОЛЕ (РОЬЕ) в целом: архитектуры, управления, баз данных и графического интерфейса ввода/вывода данных.

4. Тестирование математических моделей ЯГ и алгоритмов ЗВ-моделирования данными Монте-Карло и эксперимента.

5. Обобщение быстрых сеточных методов ЗВ-моделирования на электрические и температурные поля ГИС в сложных геометриях, включая горизонтальные и наклонные скважины.

6. Создание алгоритмов и пакета программ ДЕЛЬТА (ВЕЬТА) для оптимизации зондов и режимов измерений аппаратуры ЯГ;

7. Широкое опробование пакета ПОЛЕ для решения задач геофизической науки, производства и НИОКР.

Научная новизна. Создан новый вычислительный аппарата быстрого ЗВмоделирования полей ядерной геофизики, в т.ч. впервые:

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

2. Разработана комбинированная вычислительная схема (KB С) для быстрого 3D-моделирования полей и показаний методов скважинной гамма-спектрометрии, использующая комбинацию "медленных" вычислительных методов (моментов, Монте-Карло, оптимизации/параметризации многомерных функций) и «быстрых» методов (аналитических, аппроксимации/интерполяции, функций Грина, обобщенного спектрально-лучевого приближения (GSRA)).

3. В теории ядерных методов ГИС созданы методы, алгоритмы и пакет программ ЗБ-моделирования полей нейтронов и фотонов, которые удовлетворяют одновременно 3-м взаимно-противоречивым критериям: точности, высокого быстродействия и пригодности для геометрий каротажа любой сложности. При этом достигнутое быстродействие многократно превышает скорости известных отечественных и зарубежных программ-аналогов того же класса сложности и точности.

4. Получены точные конфигурации каротажных диаграмм ядерных методов ГИС в горизонтальных и наклонных скважинах и и найдены зависимости характерных точек диаграмм ЯГ от нового параметра - угла встречи скважины и пласта.

5. Разработаны метод, алгоритм и пакет программ для теоретической компьютерной оптимизации зондов, технологий измерений и метрологических характеристик аппаратуры ЯГ, основанный на быстром численном моделировании полей в режиме online.

Защищаемое научное положение. Разработанная совокупность комбинированных вычислительных схем, физико-математических моделей ЯГ, адаптивных сеток, вычислительных алгоритмов и баз данных, реализованная в пакете программ ПОЛЕ, позволяет выполнять ЗБ-моделирование полей и показаний ядерных методов ГИС, удовлетворяющее одновременно критериям точности, быстродействия, геометрической и ядерно-физической универсальности.

Методы исследования, личный вклад. Теоретические и математические исследования, программирование на языках С++ и Фортран, пакетное проектирование, разработка баз данных, вычислительный эксперимент, тестирование данными Монте-Карло и физического эксперимента, массовое численное моделирование полей и показаний ЯГ, оптимизационные исследования.

Личное участие автора состоит в разработке математических методов, моделей, алгоритмов и С++-программ быстрого 3 D-моделирования геофизических полей, реализованных в действующей версии пакета программ ПОЛЕ, разработке алгоритмов, программ и интерфейса оптимизационного пакета ДЕЛЬТА, тестировании созданного вычислительного аппарата данными Монте-Карло и эксперимента, во всех опробованиях пакетов ПОЛЕ и ДЕЛЬТА при решении задач геофизической науки, производства и НИОКР.

Практическая значимость работы. Разработанные моделирующий пакет ПОЛЕ и оптимизационный пакет ДЕЛЬТА позволяют решать следующие практические задачи и этапы НИОКР:

- быстрой разработки полных и достоверных систем интерпретационных зависимостей для существующей и проектируемой аппаратуры ЯГ;

- внедрения в практику геофизических НИОКР теоретических методов оптимизации зондов, режимов измерений и метрологических характеристик аппаратуры ЯГ (пакет ДЕЛЬТА), что позволит существенно сократить сроки ее разработки;

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

Решен ряд конкретных практических задач производства и НИОКР, поставленных геофизическими организациями России:

- оптимизации зондов и создания систем палеток и поправок для малогабаритных трехзондовых приборов ННКнт и ННКт;

- методического обеспечения переинтерпретации данных ЯГ при подсчете запасов нефти в юрских отложениях Зап. Сибири;

- оптимизации зондов и интерпретационного обеспечения аппаратуры углерод-кислородного С/О-каротажа;

- исследования влияния различных типов утяжелителей бурового раствора на показания методов ядерного каротажа в разрезах баженовской свиты продуктивных отложений Зап. Сибири;

- оптимизации методических характеристик отечественного АМК многозондового импульсного нейтронного каротажа АИНК-89.

Апробация работы, публикации. Основные результаты диссертации докладывались и опубликованы в трудах:

Международной конференции "NUCLEAR GEOPHYSICS" 97" по ядерной геофизике (Краков, октябрь 1997); Всероссийской научно-технической конференции молодых ученых и специалистов "Геофизика-97" (С-Петербург, июнь 1997); Всероссийской конференции "Современная ядерная геофизика при поисках, разведке и разработке нефтегазовых месторождений" (Бугульма, май 1997); Научно-практического семинара «Проблемы качества ГИС» (Тверь, январь 1997); Международного симпозиума "Новая геофизическая техника для исследования бурящихся и действующих вертикальных, наклонных и горизонтальных скважин" (Уфа, апрель 1997); Международной конференции SPWLA-ЕАГО-РГУНГ по скважинной геофизике "МОСКВА-98" (Москва, сентябрь 1998); Всероссийского научно-практического семинара "Ядерная геофизика: состояние и перспективы" (Москва, май 1999).

Основные положения диссертации опубликованы в 10 научных статьях и тезисах конференций и 5 научно-исследовательских отчетах.

Объем и структура работы. Диссертация состоит из введения, 4 глав, состоящих из 26 параграфов, заключения и библиографии. Она содержит страниц основного текста, рисунков, таблиц, библиографию из

Заключение Диссертация по теме "Геофизика", Кулешова, Людмила Борисовна

Выводы для практической интерпретации: а) помехи для НКт четко ранжированы (за исключением углистости) в порядке убывания эффекта следующим образом: известняк (100%) -» шпаты (70%) —» пирит (5%) шпаты (40%) —» сидерит (5%) уголь (5%, Кп>15%). б) Все помехи значительны: > 2% абс. и достигают 5% абс. в) Все помехи НКт с увеличением Кп возрастают по абсолютной величине, а эффект углистости, наоборот, падает (т.к. в буром угле около половины молекул —водород) и в плотных пластах он наряду с карбонатностью является основной помехой. г) Все помехи НКнт четко ранжированы во всем диапазоне Кп от 0 до 30% следующим образом: уголь (« 4%) —» известняк (2 — 3%) —> остальные (до 0,5 — 1%).

Для выполнения данной работы была подобрана и оттестирована математическая модель метода НКт-50 с аппаратурой типа ДРСТ-3-90 в обсаженных скважинах, для которых надежных палеток практически не существует либо очень мало. Погрешность расчетных палеток по ППП ПОЛЕ оказалась около 0,3-1% в единицах Кп — по сопоставлению с экспериментальными палетками Басина-Тюкаева (модели ВНИИЯГГ). С помощью этой модели пакетом ПОЛЕ было рассчитано семейство палеток для практической интерпретации данных НКт в полевошпатных песчаниках юры Зап.Сибири с целью оценки подсчетных параметров.

4.4. Оптимизация зондов и интерпретационное обеспечение аппаратуры углерод-кислородного С/О-каротажа

Основными задачами и этапами исследования (поставлены Тверьгеофизикой, ВНИИГИС и Сургутнефтегеофизикой /102/) являлись:

I. Математическое моделирование показаний аппаратуры С/О-каротажа, разработанной во ВНИИГИС, в терригенных нефтеводонасыщенных коллекторах для широкого диапазона изменения геолого-технических условий (ГТУ): пористости (Кп), нефтенасыщенности (Кн), литологии, диаметра (Вс) и конструкции скважин, длины зонда (X) и т.п.

II. Расчет нескольких типов функционалов Б; [7] от показаний С/О-каротажа (1), которые могли бы быть использованы в качестве интерпретационных параметров при обработке данных С/О-каротажа, и исследование их зависимостей от варьируемых ГТУ.

III. Сравнительный анализ всех сконструированных интерпретационных параметров Б^] с целью выбора наилучших из них для С/О-каротажа по одному из двух критериев: а) минимума погрешности определения коэффициента нефтенасыщенности Кн (или продуктивности q = Кп Кн); б) постоянства интерпретационного параметра на вмещающих («фоновых») породах: плотных песчано-глинистых интервалах, чисто водонасыщенных пластах с переменной пористостью и т.п.

Рассмотрим постановки задач и результаты их решения для каждого из этапов I,

II, III.

I. Численное моделирование показаний С/О-каротажа

1.1. Исходный массив показаний С/О-каротажа—скорости счета ГИНР 1с и До в «углеродном» АЕс) и «кислородном» АЕо) энергетических окнах для описанного ниже диапазона ГТУ был получен с помощью математического моделирования быстрым пакетом прикладных программ (ППП) ПОЛЕ. Аналогичный массив был получен для ГИРЗ.

1.2. Приведем в табл. 4.2 заложенные в базе данных ППП ПОЛЕ спектры ГИНР, испускаемые основными породообразующими и конструктивными элементами. При этом под выходом (Ь) ГИНР данной энергетической линии Еу данного элемента X понимается величина

Е0

Ьх(Еу) = I ёЕп.ахПу(Еп-»Еу),

Епорог где стхПу(Еп-^Еу) — микроскопическое сечение генерации линии ГИНР энергией Еу на ядре элемента X в зависимости от энергии замедляющихся нейтронов Еп,

Епорог — порог (п,у)-реакции неупругого рассеяния генерации линии Еу, Е0 = 14 МэВ — энергия с1-Т источника нейтронов, все Е; — в МэВ, схПу — в миллибарнах (мб), 1 мб = 10 — 27 см2. Зависимости сечений генерации для каждой из линий, приведенных в табл. , от энергии замедляющихся нейтронов с1-Т источника в базе данных пакета ПОЛЕ взяты согласно библиотекам оцененных ядерных данных / 88 /.

Выходы линий ГИНР рассчитаны только для линий с энергией больше 3.2МэВ, которые попадают в интервалы энергетических окон аппаратуры С/О-каротажа. В табл. 4.2 в последней колонке приведены характеристики ГИНР-континуума: пороговая энергия Епорог и выход Ьконт.

Е„ конт | ¿Еп.ахК0НТ(Еп),

Епорог где ахК0НТ(Еп) — микроскопическое сечение генерации ГИНР-континуума на ядре элемента X в зависимости от энергии замедляющихся нейтронов Еп, Епорог — порог ГИНР-континуума, Ео = 14 МэВ — энергия ё-Т источника нейтронов, все Е; — в МэВ, ахконт — в миллибарнах (мб).

Спектр ГИНР-континуума принимался согласно испарительной модели ядра и имеет форму, близкую к распределению Максвелла, с т.н. температурой ядра Тх(Еп). Данные по зависимостям Тх(Еп) и сгхК0НТ(Еп) взяты из / 88 /.

Спектр ГИРЗ элементов взят из справочника / 89 /.

ЗАКЛЮЧЕНИЕ

Разработано новое поколение быстрых комбинированных методов и алгоритмов вычислительной ядерной геофизики, которое реализовано в пакете прикладных программ ПОЛЕ. Он выполняет 313-моделирование полей и показаний ядерных методов ГИС, которое удовлетворяет одновременно 4-м трудным требованиям, выдвигаемым практикой: 1)точности; 2) быстродействия; 3) ядерно-физической универсальности, т.е. включения важнейших методов ЯГ во всех модификациях (стандартных, спектрометрических, импульсных, многозондовых и любых их сочетаниях); 4) геометрической универсальности, т.е. пригодности для любых геометрий измерений при каротаже, включая горизон-тальные и наклонные скважины.

При этом решены следующие задачи (из них ряд впервые):

- разработаны физико-математические и петрофизические модели численной теории методов ЯГ; построены алгоритмы формирования адаптивных сеток и разностной аппроксимации краевых задач ЗБ-моделирования ядерных и других полей ГИС в сложных геометриях, включая горизонтальные и наклонные скважины;

- разработаны быстрые комбинированные методы и алгоритмы ЗБ-моделирования полей нейтронов и гамма-квантов в ядерной геофизике. В ПОЛЕ быстрые вычислительные методы: аналитические, конечные разности, многомерная аппроксимация, функций Грина - входят в пакет непосредственно в виде алгоритмов моделирования, а медленные: Монте-Карло, метод моментов и эксперимент - в виде результатов моделирования и априорной информации.

- создан графический интерфейс ввода/вывода данных для ППП ПОЛЕ;

- математические модели ЯГ и алгоритмы ЗБ-моделирования полей тестированы данными Монте-Карло и физических экспериментов.

- разработанные быстрые сеточные алгоритмы 3 D-моделирования ядерных полей обобщены на электрические и температурные поля ГИС в классе геометрий, включающих горизонтальные и наклонные скважины.

Для оптимизации зондов и режимов измерений, а также оперативной разработки интерпретационного и метрологического обеспечения аппаратуры ядерного каротажа на стадиях НИОКР и производственного применения разработан пакет программ ДЕЛЬТА, в котором роль интеллектуального ядра играет быстрое ЗБ-моделирование полей в режиме on-line пакетом ПОЛЕ. Выпол-нено широкое опробование пакетов ПОЛЕ и ДЕЛЬТА: они примене-ны для решения ряда задач геофизической науки, производства и НИОКР, поставленных геофизическими организациями России. В т.ч. оптимизированы характеристики аппаратурно-методических комплексов многозондового импульсного нейтронного каротажа и С/О-каротажа.

Впервые рассчитаны точные конфигурации каротажных диаграмм ядерных методов ГИС в горизонтальных и наклонных скважинах, позволившие исследовать физические основы и построить важные элементы интерпретации данных ЯГ в таких системах, в т.ч. зависимости характерных точек диаграмм ЯГ от нового параметра - угла встречи скважины и пласта.

Библиография Диссертация по геологии, кандидата физико-математических наук, Кулешова, Людмила Борисовна, Москва

1. Страхов В.Н. Методологические проблемы теории и практики интерпретации данных в прикладной геофизике. Тр. конф. "Вопросы методологии интерпретации геофизических данных в прикладной геофизике" (Москва, 7-8 февраля 1996). М., изд. РАН-ЕАГО, 1996, с.4-20.

2. Страхов В.Н. Научное мировоззрение. 1,2. Геофизика N1, 1993; N2, 1994.

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

4. Поляченко А.Л. Численные методы в ядерной геофизике. М., Энергоатомиздат, 1987, 152 с.

5. Дядькин И.Г. Методы Монте-Карло в физике. В кн."Методы Монте-Карло в физике и геофизике". Уфа, БашГУ, 1973, с.7.

6. Дядькин И.Г. Вычислительный эксперимент Монте-Карло, моделирующий перенос нейтронов и гамма-квантов. В кн."Численные методы в ядерной геофизике. М., Энергоатомиздат, 1987, с.67-110.

7. Хисамутдинов А.И. и др. Алгоритмы Монте-Карло в ядерной геофизике. Новосибирск, Наука СО АН СССР, 1985.

8. Velizhanin V.A., Dyadkin I.G., Enikeeva F.H. et.al. Monte Carlo simulation in nuclear geophysics. I, II. Nuclear Geophysics, v.3, No.4, 5, 1991.

9. Кожевников Д.А. Нейтронные характеристики горных пород. М., Недра, 1982.

10. Давыдов Ю.Б., Кузин В.Ф. Теоретические предпосылки каротажа нейтронов деления. Новосибирск, ВО Наука, 1994, 352 с.

11. Проблемы и перспективы ядерно-геофизических методов в изучении разрезов скважин. Тезисы докл. Всесоюзной научно-техн. конф. (Обнинск, 19-24 июня 1989). МНТК "ГЕОС", 1989.

12. Поляченко А.Л., Кулешова Л.Б. и др. Развитие оптимизационного "беспалеточного" метода интерпретации данных ядерной геофизики на базе объектно-ориентированного проектирования. В сб. /12/, докл. F3.6.

13. A.L.Polyachenko, I.V.Kozyrev,L.B.Kuleshova, A.G.Timofeev. Optimization of <Chart-Free> interpretation of nuclear logs based on the object-oriented designing.

14. Кожевников Д.А. Спектрометрическая и интегральная модификации гамма-метода: алгоритмическая интерпретация. В сб. /9/, с.271.

15. Кожевников Д.А. Адаптивные и поправочные методики интерпретации данных радиометрии скважин. В сб. /12/, докл. F2.4.

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

17. Briemlister J.F. /Ed./ MCNP a general Monte-Carlo code for neutron and photon transport, version ЗА. Los Alamos Nat. Lab., report LA-7396-M, rev.2, 1986.

18. Shuttleworth Е., Chucas S.J. Linked Monte-Carlo and finite-element diffusion methods for reactor shield design (McBEND). Proc. 6 Int. Conf. on radiation shielding. Tokyo, v.I, 1983, p.180.

19. Butler J., Clayton C.J. A new philosophy for calibration oil well logging tools based on neutrons trasport codes. Trans. 25-th SPWLA annual logging symp., New Orleans, 1984, v. II.

20. K.-J Dunn. A diffusion model for pulsed neutron logging. Geophysics, v.54, No.l, 1989, 100-113.

21. Mickael M., Gardner R.P., Verghese K. McDNL: a new specific purpose Monte Carlo code for simulation of dual-spaced neutron porosity logs. SPWLA-29, 1988, v. II.

22. Gardner R.P. et.al. Complite composition and density correlated sampling in the specific purpose Monte Carlo codes McPNL and McDNL for simultating pulsed neutron and neutron porosity logging tools. Nucl.Geophys., v.3, No.3, 1989, p. 157-165.

23. Еникеева Ф.Х., Журавлев Б.К., Гулин Ю.А. Решение задач нейтронного каротажа нефтяных скважин. В сб. /22/, с.34-56.

24. Журавлев Б.К. Диссертация к.т.н., МИНХиГП, 1985.

25. Головацкий С.Ю. Особенности вычислительного эксперимента для задач гамма-каротажа. В кн. Методы Монте-Карло в вычислительной математике и математической физике (докл. 7 Всесоюз. совещ. Новосибирск, ВЦ СО АН СССР, 1985, с.340-343.

26. Математическое моделирование в ядерной геофизике. Уфа, изд. БашФАН СССР, 1979, 183 с.

27. Лухминский Б.Е. Диссертация д.ф.м.н., МГРИ, 1985.

28. Лебедев В.Е., Денисик С.А. и др. Об интегрировании уравнения диффузии нейтронов методом сеток применительно к задачам ядерной геофизики. Труды ВНИИЯГГ, вып. 7, 1969, с.20-25.

29. Поляченко А.Л. Быстрое математическое моделирование задач ядерной геофизики: пакет ПОЛЕ. НТВ "Каротажник", Тверь, изд. АИС, 1996, вып.28, с.72-85.

30. Поляченко А.Л., Кулешова Л.Б. Математические модели импульсного и стационарного нейтронного каротажа при решении прямых и обратных задач пакетом ПОЛЕ. "Каротажник", Тверь, изд.АИС, 1996, вып.29, с.91-104.

31. А.Л.Поляченко, Л.Б.Кулешова. Быстрое математическое моделирование прямых задач ядерного каротажа и других методов ГИС.

32. A.L.Polyachenko, L.B.Kuleshova. Fast mathematical modelling of direct problems for Nuclear and other Well Logging methods. В сб. /33/, с.

33. Басс Л.П., Гермогенова Т.А., Хмылев А.Н. Модульная структура программ в осесимметричных задачах теории переноса. Система "Радуга". Препринт ИПМ АН СССР N97. М.; 1973

34. Басс Л.П. и др. Пакет ЗАЩИТА. Версия 3.0. Основные возможноти функционального и системного наполнения. Препринт ИПМ АН СССР N27. М.; 1987, 24с.

35. Башмачникова А.Н., Зизин М.Н. и др. ФИХАР модульная система программ для реакторных расчетов. Труды III семинара по комплексам программ математической физики. Новосибирск, 1973

36. Марков Ю.В., Зизин М.Н. и др. Структура комплекса НФ-6 для реакторных расчетов на БЭСМ-6. Комплексы программ математической физики. Новосибирск, 1978.

37. Горбенко Н.И., Ильин В.П. и др. Пакет программ ЭРА для автоматизации электрооптических расчетов. Численные методы решения задач электронной оптики. Новосибирск, ВЦ СО АН СССР, 1979, с.34-60.

38. ППП ЭФЕС-3 для решения трехмерных краевых задач. Препринт N 561 ВЦ СО АН СССР. Новосибирск, 1985.

39. Борисов В.М. Разработка пакетов программ вычислительного типа. -М: Изд. МГУ, 1990, 124с.

40. Пакеты прикладных программ. Проблемы и перспективы (Сер. Алгоритмы и алгоритмические языки). М., Наука, 1982.

41. Абагян Л.П. и др. Групповые константы для расчета реакторов и защиты (БНАБ-78). М, Атомиздат, 1982, 311 с.

42. Абагян Л.П. и др. Групповые константы для расчета ядерных реакторов (БНАБ-64). М., Атомиздат, 1964.

43. ENDF/B-IV cross-sections measurements standarts/ Garber D. et al. BNL-NCS-50464, ENDF-225, N.-Y., 1975.

44. Garber D.I., Kinsey R.R. Neutron cross section curves. BNL-325, 3-rd ed., v.2, 1976.

45. Internatinal conference "NUCLEAR GEOPHYSICS'97" (Krakow, 20-23 Oct. 1997). Abstracts, Krakow, 1997.

46. Годунов С.К., Рябенький B.C. Разностные схемы (введение в теорию). М., Наука, 1977, 440 с.

47. Самарский A.A., Николаев Е.С. Методы решения сеточных уравнений. М., Наука, 1978,-590 с.

48. Вычислительная математика и техника в разведочной геофизике. Справочник геофизика, изд.2 /Под ред. В.И.Дмитриева. М., Недра, 1990, 498 с.

49. Качанов И.Н. Оценка групповых выходов гамма-излучения при радиационном захвате тепловых нейтронов. Вопросы дозиметрии и защиты от излучений. Тр. МИФИ, вып. 20, М., 1981, с.43-48.

50. Тихонов А.Н., Самарский A.A. О разностных схемах для уравнений с разрывными коэффициентами. ДАН СССР, т. 108, вып.З, 1956; ДАН СССР, т. 124, N 4, 1959.

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

52. Гермогенова Т.А., Дегтярев С.Ф., Орлов В.В. и др. Перенос быстрых нейтронов в плоских защитах. М., Атомиздат, 1971.

53. Поляченко А.Л. Теория скважинной ядерной геофизики. В справочнике "Скважинная ядерная геофизика", гл.1, Недра, изд.2, 1990, с.6-32.

54. Поляченко А.Л. Диссертация д.ф.м.н. "Теория обратных задач нефтегазовой скважинной ядерной геофизики". ВНИИЯГГ, 1983.

55. Ильин В.П. Численные методы решения задач электрооптики. Новосибирск, Наука, 1974.

56. Ильин В.П. Численные методы электрофизики: Наука, 1986.

57. Городков С.С. Сравнение сеточных алгоритмов расчета нейтронных полей в больших реакторах. Атомная энергия, т.85, вып.З, 1988, с. 184.

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

59. Булеев Н.И. Пространственная модель турбулентного обмена. М.,Наука, 1989.

60. Марчук Г.И. Методы вычислительной математики. М., Наука, 1977, 456с. Марчук Г.И., Лебедев В.И. Численные методы в теории переноса нейтронов.1. S3 21. М., Атомиздат, 1971.

61. Вабищевич П.Н., Макаров М.М. Институт математического моделирования РАН. Препринт N 19, 1993.

62. Поляченко A.JL, Гаврина Т.Е., Шапошникова Т.А. Расчет нейтронных групповых констант горных пород (программа КОНГ). М., ВНТИЦентр-ГФАП, "Алгоритмы и программы", N 3, П001377, 1976.

63. Николаев М.Н., Савоськин М.М. и др. Многогрупповое приближение в теории переноса нейтронов. М., Энергоатомиздат, 1984, 256 с.

64. Вычислительные методы в физике реакторов/Под ред. Х.Гринспена. М., Атомиздат, 1972, 372 с.

65. Физические проблемы защиты реакторов. М., Атомиздат, 1971.

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

67. Вигнер Е., Вейнберг А. Физическая теория ядерных реакторов. М., ИЛ, 1961.

68. Бергельсон Б.Р., Суворов А.П., Торлин Б.З. Многогрупповые методы расчета защиты от нейтронов. М., Атомиздат, 1970.

69. Барсуков O.A., Авзянов B.C. Пространственно-энергетическое распределение нейтронов в системе скважина-пласт. Атомная энергия, т.10, вып. 5, 1961, с.478-486.

70. Поляченко А.Л., Шапошникова Т.А. и др. Альбом расчетных палеток нейтронного каротажа. ВНИИЯГГ, 1977, 73 с.

71. Зендриков Ю.М., Поляченко А.Л. и др. Расчет показаний нейтронных методов при изучении слоистого разреза. Сб. Ядерно-геофизические и геоакустические методы выделения продуктивных пластов в обсаженных скважинах. Тр.ВНИИЯГГ, вып.11, Недра, 1972, с.23-32.

72. McCallen C.W.J. SNAP-3D three-dimensional neutron diffusion code. TRG report 2677R/U.K.A.E.A., Risley, 1975.

73. Badruzzman A., Chiarmonte J. A comparison of Monte Carlo and discrete ordinates methods in three-dimensional well-logging problem. Trans. Amer. Nucl. Soc., v. 50, 1985, p.256-267.

74. Кулешова Л.Б., Поляченко А.Л. Быстрое математическое моделирование прямых задач ядерной геофизики. "Геоинформатика", N 2, 1998, с.49-53.

75. Кулешова Л.Б. Быстрое ЗО-моделирование нейтронных полей в ядерной геофизике. Всеросс. научно-техн. конф. молодых ученых и специалистов "Геофизика-97" (С-Петербург, 3 -6.06.97). Тезисы докладов. С.-ПГУ-ЕАГО, С-Петербург, 1997, с. 19.

76. Самарский A.A., Гулин A.B. Численные методы. М., Наука, 1989.

77. Самарский A.A., Андреев В.Б. Разностные методы для эллиптических уравнений. М., Наука, 1976.

78. Самарский A.A. Теория разностных схем. М., Наука, 1983.

79. Тихонов А.Н., Самарский A.A. Однородные разностные схемы на неравномерных сетках. ЖВМ и МФ, 1962, 2, вып.5.

80. Тихонов А.Н., Самарский A.A. Уравнения математической физики. М., Наука, 1972.

81. Булеев Н.И. и др., Разностные уравнения диффузии на границах раздела двух сред. "Вопросы атомной науки и техники". Сер. Физика и техника ядерных реакторов, Вып. 7 (29), М., 1982, с.14-18.

82. Богданова М.С., Капорин И.Е., Николаев Е.С. и др. Комплекс программ для решения разностных эллиптических краевых задач. Пакеты прикладных программ: проблемы и перспективы. М.; Наука, 1982.

83. Ильин В.П. Об одном варианте системы модульного программирования. Препринт ВЦ СО АН СССР, N9. Новосибирск, 1976

84. Goldstein Н., Wilkins J.E. Jr. calculations of the penetration of gamma-rays. NYO-3075, June 30, 1954.

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

86. Гольдштейн Г. Основы защиты реакторов. М., Атомиздат, 1961, 343с.

87. Бродер Д.Л., Попков К.К., Рубанов С.М. Малогабаритная защита реакторов. М., Атомиздат, 1967, 365с.

88. Поляченко А.Л. Комбинированный вычислительный аппарат решения прямых задач скважинной гамма-спектрометрии. В / 9 /, с.66.

89. Поляченко А.Л. и др. Отчеты НИР по темам NN 69, 222, 1/45 контрактов с Мингео СССР. М., ВНИИГеосистем, 1989 93 .

90. Поляченко А.Л. Состояние и задачи развития теории и вычислительного аппарата ядерного каротажа. В 191, с. 34.

91. Калиткин H.H. Численные методы. М., Наука, 1976.

92. Головацкий С.Ю. Диссертация к.т.н. "Разработка системы метрологического обеспечения интегрального гамма-каротажа нефтегазовых скважин. М., ВНИИГеоинформсистем, 1989.

93. А.Л.Поляченко, С.А.Кантор, К.С.Михалев, А.Г.Шпаер. Комбинированные методы решения прямых задач скважинной гамма-спектрометрии.

94. A.L.Polyachenko, S.A.Kantor, K.S.Mihalev, A.G.Shpaer. Combinated methods for solution of well gamma-spectrometry direct problems. (Краков, 20-23 октября 1997). Тез. докл.

95. Злоказов В.Б. Малопараметрический способ аппроксимации спектров. Препринт ОИЯИ PI 1-86-135, Дубна, 1986, Юс.

96. Кантор С.А. Аппроксимация функции отклика сцинтилляционного гамма-спектрометра по характерным точкам. "Вопросы атомной науки и техники". Сер. Радиационная техника, вып.З (43), 1990, с. 15-20.

97. Дахнов В.Н. Интерпретация результатов геофизических исследований разрезов скважин. М., Недра, 1972, 365с.

98. Поляченко А.Л., Кулешова Л.Б., Ломтев И.Л. и др. Численное решение прямых задач ядерного, электрического и температурного каротажа для горизонтальных и наклонных скважин (пакет ПОЛЕ).

99. A.L.Polyachenko, L.B.Kuleshova, I.L.Lomtev et al. В /12/, докл. 11.6. Numerical Solution of Direct Problems of Nuclear, Electrical, and Temperature Logging in Horizontal and Deviated Wells: POLE Software Packages. В /12/, доклад 11.6.

100. Чекалюк Э.Б. Термодинамика нефтяного пласта. М.,Недра, 1965, 283 с.

101. Филиппов А.И. Скважинная термометрия переходных процессов. Диссертация д.ф.м.н., Уфа, 1990.

102. Валиуллин P.A. Дисс. д.т.н. "Термические методы диагностики нефтяных пластов и скважин". Тверь, НТП ГЕРС, 1996.

103. Кулинкович А.Е. К решению задачи теории нейтронного каротажа методом «групп». Ядерная геофизика. М., Недра, 1963.

104. Самарский A.A., Карпов В.Я., Корягин Д.А. Принципы разработки пакетов прикладных программ для задач математической физики. ЖВМиМФ, 1978, т. 18, N2.

105. Ершов А.П., Ильин В.П. Пакеты программ технология решения прикладных задач. Препринт ВЦ СО АН СССР, N121. Новосибирск, 1978

106. Горбунов-Посадов М.М. и др. Системное обеспечение пакетов прикладных программ / под ред. А.А.Самарского. М., Наука, 1990,- 208 с.

107. Горбунов-Посадов М.М., Карпов В.Я. и др. Пакет САФРА: программное обеспечение вычислительного эксперимента. Сб. "Пакеты прикладных программ". Вычислительный эксперимент. М., Наука, 1983.

108. Рвачев B.JL, Манько Г.П. и др. Организация функционального наполнения в генераторе программ "ПОЛЕ-ЗМ". В кн. Пакеты прикладных программ. Системное наполнение. М., Наука, 1984, с.89-93.

109. Рвачев В.Л., Манько Г.П. Генераторы программ серии ПОЛЕ для решения задач математической физики. В кн. Численные методы механики сплошной среды. Новосибирск, ИТПМ СО АН СССР, т.12, N3, 1981.

110. Сторм Э., Исраэль X. Сечения взаимодействия гамма-излучения. Справочник. М., Атомиздат, 1973, 265 с.

111. Труханов Г.Я. и др. Ядерно-физические константы взаимодействия нейтронов с элементами, входящими в состав атмосферы и земной коры. М., Энергоатомиздат, 1981,- 304 с.

112. Бычков В.М. и др. Сечения пороговых реакций, вызываемых нейтронами. М., Энергоатомиздат, 1982, 216 с.

113. Разведочная ядерная геофизика: Справочник геофизика/ Под ред. О.Л. Кузнецова и А.Л. Поляченко. 2-е изд., М., Недра, 1986, - 432 с.

114. Абагян A.A. и др. Вторичное гамма-излучение в радиационной защите. М., Энергоатомиздат, 1984, 176 с.

115. Барышев Ю.Ф., Миллер В.В. Метод определения среднего времени жизни нейтронов. В кн. Новые методы и аппаратура ядерной геофизики и геохимии. М., ВНИИЯГГ, 1970, с.63-71.

116. Кожевников Д.А., Марьенко H.H. и др. Экспериментальное изучение распределения нейтронов в однородных средах. В кн. "Ядерно-физические методы изучения горных пород и насыщающих их флюидов". Труды МИНХ и ГП, вып. III, М., Недра, 1974, с.40-63.

117. Гулин Ю.А. О характере зависимости показаний нейтронного каротажа от пористости пород. Прикладная геофизика, вып.77,1976, с.204-211.

118. Иванкин В.П. Тр. HB НИИГГ, вып.4, Саратов, 1966, с.301-321.

119. Амурский А.Г., Блюменцев A.M. и др. Разработка аппаратуры многозондового импульсного нейтронного каротажа. Доклад на /109 /.1. S3?

120. Гулин Ю.А. Гамма-гамма методы исследования нефтяных скважин. М., Недра, 1976,- 160 с.

121. Поляченко А.Л., Цейтлин В.Г. Оптимизация режимов измерений и параметров аппаратуры импульсного нейтронного каротажа. В кн. Ядерно-геофизические и гоакустические методы поисков и разведки месторождений полезных ископаемых. М., ВНИИЯГГ, 1975, с.143-160.

122. Резванов P.A. Радиоактивные и другие неэлектрические методы исследования скважин. М., Недра, 1982, 367 с.

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

124. Кучурин Е.С., Кулешова Л.Б., Поляченко А.Л. и др. Анализ интерпретационных параметров углерод-кислородного каротажа.

125. E.S.Kuchurin, L.B.Kuleshova, A.L.Polyachenko, V.V.Vinokurov. Analysis of С/О Log Interpretation Parameters. В /12/, доклад ВЗ.З.

126. Кулешова Л.Б., Поляченко А.Л. Отчет о НИР по договору с ВНИИГИС «Разработка средств и способов интерпретации данных аппаратуры С/О-каротажа втерригенных нефтеводонасыщенных коллекторах». М., ВНИИГеосистем, 1997.

127. Кулешова Л.Б., Поляченко А.Л. Отчет о НИР по договору с трестом Сургутнефтегеофизика «Исследование влияния различных типов утяжелителей бурового раствора на показания приборов 2ННК, ГК, ГГК». М., 1998.

128. Поляченко А.Л., Бородин С.Г., Кулешова Л.Б. Отчет о НИР по договору с ВНИИГИС "Теоретическое и интерпретационное обеспечение новой многозондовой аппаратуры импульсного нейтронного каротажа (ИННК-М)".1. М., ВНИИГеосистем, 1995.

129. Отчет о НИР по договору с ВНИИЯГГ «Разработка палеточного обеспечения обработки данных ИНК при двухэкспоненциальном представлении сигнала ИНК». М., ВНИИГеосистем, 1994.

130. Поляченко А.Л., Кулешова Л.Б. Ядерно-каротажные диаграммы в горизонтальных и наклонных скважинах.

131. A.L.Polyachenko, L.B.Kuleshova. Nuclear Logs from Horizontal and Deviated Wells. В / 12/, доклад 11.7.

132. Кулешова Л.Б., Поляченко А.Л. Отчет о НИР по договору с НПФ ГОРИЗОНТ N 39-97 «Разработка интерпретационного обеспечения зондов нейтронного каротажа аппаратуры АМК «ГОРИЗОНТ» для исследования горизонтальных скважин». М., ВНИИГеосистем, 1997.

133. Всероссийский научно-практический семинар "Ядерная геофизика. Современное состояние и перспективы развития". Москва, ВНИИГеосистем , 18-20 мая 1999 .

134. Кухаренко Н.К., Басин Я.Н., Тюкаев Ю.В. Методика определения пористости карбонатных пластов по данным нейтронного каротажа (временное наставление). М., ВНИИЯГГ, 1968, 110 с.

135. Басин Я.Н., Тюкаев Ю.В. Методические рекомендации по проведению исследований и интерпретации данных нейтронного каротажа с серийной аппаратурой PK. М,. ВНИИЯГГ, 1979.

136. Дрожжинов Ю.Н., Москалев О.Б. Поток однократно рассеянных нейтронов от точечного источника в бесконечной среде. Доклады АН СССР, 1970, т. 193, N6, с.1271.

137. Дахнов В.Н. Электрические и магнитные методы исследования скважин. М., Недра, 1969.

138. Абрамович М., Стиган И. Справочник по специальным функциям (пер. с англ.). М., Наука, 1979.

139. Люк Ю. Специальные математические функции и их аппроксимации. М., Мир, 1980.

140. Калмыков Г.А., Миллер В.В. и др. Многоканальный спектрометрический гамма-каротаж ( аппаратура, методика, программное обеспечение, метрология). Доклад на /109 /.

141. Кулешова Л.Б. Оптимизация зондов и режимов измерений аппаратуры ядерного каротажа. Доклад на Всерос. науч-практ. семинаре "Ядерная геофизика Современное состояние и перспективы развития" (Москва, 18-20 мая 1999).