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

Определение удельного электрического сопротивления пластов в тонкослоистом разрезе по комплексу зондов электрического и индукционного каротажа

ВАК РФ 04.00.12, Геофизические методы поисков и разведки месторождений полезных ископаемых

Автореферат диссертации по теме "Определение удельного электрического сопротивления пластов в тонкослоистом разрезе по комплексу зондов электрического и индукционного каротажа"

КОМИТЕТ ПО ГЕОЛОГИИ И ИСПОЛЬЗОВАНИЮ НЕДР ПРИ ПРАВИТЕЛЬСТВЕ РОССИЙСКОЙ ФЕДЕРАЦИИ АУЧНО-ПРОИЗВОДСТВЕННОЕ ГОСУДАРСТВЕННОЕ ПРЕДПРИЯТИЕ О ГЕОФИЗИЧЕСКИМ РАБОТАМ НА НЕФТЬ И ГАЗ "НЕФТЕГЕОФИЗИКА" ВСЕСОЮЗНЫЙ НАУЧНО-ИССЛЕДОВАТЕЛЬСКИИ ИНСТИТУТ

ГЕОФИЗИЧЕСКИХ МЕТОДОВ РАЗВЕДКИ ' ,

СВНИИГЕОФИЗЖА)

На правах рукописи ПОТАПОВ АЛЕКСАНДР ПЕТРОВИЧ

УДК 550.832.7.058

ОПРЕДЕЛЕНИЕ УДЕЛЬНОГО ЭЛЕКТРИЧЕСКОГО . СОПРОТИВЛЕНИЯ ПЛАСТОВ В ТОНКОСЛОИСТОМ РАЗРЕЗЕ ПО КОМПЛЕКСУ ЗОНДОВ ЭЛЕКТРИЧЕСКОГО И ИНДУКЦИОННОГО КАРОТАЖА

Специальность 04.00.12- геофизические методы поисков и разведки месторождений полезных ископаемых

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

Москва - 1992

Робота выполнена во всесоюзном научно-исследовательском и проектно-конструкторском институте геофизических исследований геологоразведочных скважин (ВШИГИС) научно- производственного государствошюго предприятия по геофизическим работам(НПГП ГЕРС) Научный руководитель: - доктор технических наук,

профессор Н. Н. СОХРАНОВ Официальные оппоненты:- доктор технических наук

Е.В.Чаадаев: - кандидат физико-математических наук Л.а.книжнерман Ведущее предприятие:-тюменская тематическая экспедиция

• концерна"тюменьгеология"

Защита диссертации состоится _" 09_1992 г.в 12 час.

на заседании специализированного совета Д071.06.01 по защите диссертаций на. соискание ученой стопени доктора наук при; Всесоюзном научно- исследовательском институте геофизических методов разведки (ВШИГеофизика) по адресу :ЮЮОО,г.Москва, ул.Чернышевского,22

С диссертацией иокно ознакомиться в библиотеке института.

Автореферат разослан _" 08_1992 г.

Отзыв на автореферат в двух экземплярах, ¡заверенный печатью

учреадэния, просим приолеть ученому оекротирю.

Ученый секретарь специализированного совете,, кандидат технических наук

Н. П.Чихов

Обшая характеристика работы Актуальность проблемы Данные геофизических исследований скввмш (ГИС),среди которых наиболее распространенными являются электрические, дают значительный объем Ш1формации о геологическом разрезе. Удельное электрическое сопротивление(УЭС) несот информацию о таких характеристиках плскта.как структура порового проо странства и характер насыщения.

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

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

На этом пути весьма незначительные резервы повышения точности определения УЭС в тонкослоистом разрезе,так как в принципе модель не списывает реальный геоалектрический разрез.Значителышй прогресс может быть достигнут при использовании современных методов решения прямых задач для модели пачки пластов с проникновением(тонкослоистый разрез) и разработка на их основе методов решения обратных задач, способов и алгоритмов интерпретации данных электрокаротака в условиях тонкослоистого разреза.

Новые возможности открываются в связи с применением при обработка материалов ГИС персональных ЭВМ, оснащенных графическими мониторами,что позволяет аффективно использовать диалоговый режим.

Разработка новых систем обработки данных ГИС ставит новые

задачи по совершенствованию технологии определения УЭС в тонкослоистом разрезе, что повисит эффективность и точность определения УЭС.

Цель работы.Разработка способов и программного обеспечения определения УЭС в тонкослоистом разрезе на основе решений прямых и обратных задач с оценкой погрешности полученных результатов.

Основные задачи исследования:

• -исследовать возможности разработки эффективных алгоритмов определения УЭС в тонкослоистом разрезе на основе решений прямых и обратных задач;

-разработать метода решения прямых и обратных задач БКЗ, ИК.БК, которые возможно использовать в алгоритмах интерпретации влектрокаротажа в тонкослоистом разрезе;

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

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

Научная новизна: -получены аналитические выражения для вычисления матриц производных по искомым параметрам различных методов электрического и индукционного каротажа, необходимых при решении обратных задач;

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

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

Защищаемые полозапия: -аналитические выражения для вычисления матриц производных по искомым параметрам, полученлио на основе линеаризации исходного уравнения Риккати, математическое моделирование измерений различными зопда?«1 элэктричоского и индукционного каротажа, позволяющие существенно уменьшить время счета при решении поставленной задачи на ЭВМ;

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

Практнчэская цо!шесть работы зяклачаэтея о том,что результаты исследований и разработанное программное обеспечение используется при интерпретации материалов электрокаротажа в тонкослоистом разрезе,что позволяет получить болов точное значение УЗС пород, за счет использования в качества интерпретационной модели тонкослоистого разреза, итерационных методов решения обратных задач с оценкой достоверности результатов. Разработанное программное обеспечение. используется на ЭЕМ тала ЕС, 1ЕМ РО п различных системах:АС0ИГИС, "ГИНТЕЛ",

••ингис", «подсчет-.

Апробация работы.Остишо полошим настоящей работы докладывались и обсуждались иа научно-технических конференциях посвященных вопросам поиска и разводки нефтянных,газовых и газокон-. денсвтнцх месторождений йанвдной Сибири 1фи геофизических исследовшшях скважин (г.Тюмень, 19&4Г. ,1$йСг. ,1907г. ,19Шг.), на научно - технической коифорошщи "Повышение эффективности геофизических исследований скважин" (г.Октябрьский, 1Э9ог.),на совещании "Применение УШ для обработки геолого-геофизической информации при исследовании скважин"(г.Тверь,1991г.),на совещании "Персональные ЯВМ в I громы слоиои г<)офи:щко"(г.Тюмонь, 1992г. )на семинаре ^Программно-мотодичоокий комплекс геологичоской 4интпрпротации ГйС "(г.Октябрьский, 1092г.) и рядо других конференций и совещаний.

Публикации. По теме диссертации опубликовано 12 печатных робот.

Объем работы. Диссертация состоит из введения, 4 глав и заключения, содержит 110 стрпииц машинописного токста, 41 рисунок, I таблицу. Список литоратуры включает 94 шшмоноьания.

Диссертационная работа выполнена автором в пприод работы во ВНИИГКС и учебы в заочной аспирантуре ЫШИГеофизики в 1987-19Э1ГГ

Содержание работы

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

Автором были получены окончательные выражения для вычисления матриц производных по искомым параметрам,реализован алгоритм моделирования измерений различными зондами БКЗ,ПК,БК,разработан

конкретный лгосоЗ и алгоритм определения УЭС.

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

Разработка эффективных ьычисл!телышх схем решения прямых задач для моделрЯ,огтсившгапи тонкослоистый разрез,более широкое распространенно ЭВМ ьнсокой производительности позволяют распутывать на создание аффективных алгоритмов определения сопротивления.

Автор глубоко признателен научному руководителе, -доктору технических наук,профессору Оохрчнову Н.Н. за руководство, помещу н оссувдонкр направлений исслэдоьянмй на всех этапах, кандидату технических неук Кноллеру Л.Е., под р/ководством и при непосредственней участии которого выполнена работа и выбраны направления исследований.

Автор выракаот .свои искрешшэ благодарность сотрудникам отдела ' теории кг.рст.а и интерпретации данных ВНИИГИС за участие при разработке программ, за помощь при реализации технологии использования программного обеспечения при комплексной геологической ииторпротгщци, подготовке и оформлении м^ториалов диссертации .

И первоЗ глаиэ проведен анализ возможностей разработки cnocccia определении УЭС в тонкослоистом разрезе. Рассмотрено современное состояние решений прямых и обратных задач.

Широкое практическое применений получили аналитические методы (Антонов Ю.Н..Дмитриев В.И..Свдорчук А.И.,Чаадаев Е.В.), ка осново которых рассчитаны и получили признание палетки БКЗ,11К,БК (Аксельрод G.M. .Альпин Л.М. .Бондяренко И.Т., Бричепко

И.П./Зефиров H.H.. Йурйшюа D.H.,Кучеров P.A..Сохранов H.H., Сидорчук А.И., Чаадаев Е.Ь. и др.).

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

Созданные к настоящему ьромини алгоритмы . рошения прямых задач для срсд с разнотипными границами раздела на основе интегральных уравнений, коночных разностей и конечных оломеьтов (Альшш Л.М., Дашовский Ю.Л., Вильго U.U..Захаров Е.Ь..Колосов

A.А.,Морик В.Р., Иантьхии Fi.А. /Гамровский Л.А.,Зпов М.И., Ярмахов И.Г.и др.) в основном используются при моделирований» измерений для целей изучения влияния различных факторов на показания мотодов и построения ипториротациошшх палеток.

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

B.Л., Иванов И. Т. .Кшшюрман Л .Л. /Гамярчсшко Т.В. и др. ),и могут быть использованы при автоматической интерпретации. Они особенно аффективны для моделей с малым числом границ, для многослойных моделей (тонкослоистый ризроз) время счета увеличивается пропорционально числу пластов. аффективное их использование требует разработки специальных процедур интерпретации.

При использовании ЭВМ для решения обратных задач ЭК сформировались два различных подхода к peujeioco отой проблемы, и, ооотватствэнно резработанзше алгоритмы и программы уоловно моаао разделить на две группы.

К первой группе относятся алгоритмы, модолирущие процадуру

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

Fi наиболее общей постановке задача интерпретации данных комплекса ЕКЗ.ИК.БК может ставиться как обратная задача восстановления сопротивления сроды p(r,z,(p) по результатам измерений на осу. сквакшш. Известно (Друскин В .Л.), что при такой постановке задача имозт единственность решения, но при ее практическом решении встречаются определенные трудности. Зто связано, в первую очередь, с некорректностью самой задачи; с другой стороны, существушие решешия прямых задач для такого класса моделой недостаточно эффективны и требуют значительных ресурсов ЭВМ по времени счета и по' памяти.

Известен подход с использованием методов поверхноетних зарядов при рошетп! обратной Еадачи для сред сложной геометрии (Кулишсович А.Е.).

Получено численное решение двумерной обратной задачи для моделей пачки пластов с прошкновениом, когда известно нологениэ гр аниц (Друскин В.Л. ,К:шжнирмак Л.А.). П этом случае двумерная обратная задача рззбивзэтея на (1 обратных одномерных задач (гдв .V- число пластов). Но в практики обработки данных ЭК наибольшее применение получили iumcc моделей пласта бесконечной толщины и одиночного пласта с проникновением (Зупдэлевич С.М.,Комаров С.К.,Кнеллер Л.Е., Краспогюн М.Д..Пантюхин В.А.,Сохранов H.H. п Др.).

Широко используются при автоматической интерпретации различные аппроксимации решений прямых и обратных задач:обратные фильтры, мультипликативные поправки,уравнения с эмпирическими коэффициентами ( Аксельрод С.М.,Лфонша Н.Ы., Ваксман К.Г., Вагина E.H..Зверев Г.Н..Ингерман В.Г.,Кашик А.С.,Кнеллер Л.Е., Плюснин М.И.,Чуринова И.М.,и др.).

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

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

Вторая глава посвящена разработке способа определения удельного электрического сопротивления в тонкослоистом разрезе.

Эффективно решать подобные многомерные задачи позволяет использование регуляризовашшх методов оптимизации для решения обратных задач и прямых задач для сложных моделей,описывающих тонкослоистый разрез(Кнеллэр Л.Е.;.Формально необходимо подобрать такие параметры геоэлектрического разреза ( pnt, pant -соответственно сопротивления пласта и зоны проникновения ;I>3ni-диаметр зоны проникновения ; t=i,...,N.. где if-чиоло пластов),чтобы теоретически рассчитанные по прямым задачам кривые кажущегося сопротивления (КС) различных методов (БКЗ.ИК.БК). минимально

отличались от фактических кривых КС, с учетом погрешностей в измерениях.

Непосредствешга процедуре подбора параметров геовлек-трического разреза потребует многократного обращония к двумерной прямой задаче (ДПЗ), что исключает использование такого подхода при практической интерпретации материалов ГИС из-за значительных* затрат времени ЭВМ. Так, на одну итерацию необходимо (Я + 1) обращение к ДПЗ для вычисления,на основе конечных прлгращений, матриц производных (МП) по искомым параметрам, которые являются основой в итерационном процессе, и непосредственно для определения кривой КС для конкретного метода.

В работе приведен следующий подход. На первых шагах итерационного процесса,когда,как правило,задача интерпретации решается при больших значениях дисперсии (разница мввду теоретически рассчиташгами и фактически?.™ кривыми КС),используются быстрые алгоритмы решения прямых задач (БПЗ) на основе комбинации Фундаментальных решений для сред с вертикальной и радиальной неоднородностью, что позволяет с достаточной для практических целей точностью моделировать измерения различных зондов электрокаротажа. Такой подход позволил получить аналитстаскив выражения для вычисления матриц производных за сопоставимое с прямой задачей время очотп. По доо'лшвяшо некоторого минимума дисперсии БПЗ настраивается по ДПЗ, и процесс продолжается , пока не будет достигнут неснихаемый минимум дисперсии.

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

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

Для реализации приведенного спосооа решены прямые задачи различных методов электрокаротажа (ЬКЗ.ИК.БК) для сред с непрерывным распроделением сопротивления по глубине,для радиалыю-неоднородных сред и , как частный случай, для срэд с плоско-параллельными и коаксиилыючушшдрическими границами раздела, на основы которых возможно получить аналитические выражения для вычисления матриц производных по искомым параметрам.

Решение исходных уравнений непрерывности (ЭК) и Гельмгольца (МК)сведены к решению задачи Коши для уравнений Ршгкати. Для вычисления МП используется метод локальной линеаризации, который заключается в следующем. Если внести некоторые малые возмущения (йР ) в параметры модели - . Рзп1.дзп4)|Т0 8Т0 приведет к малому возмущению в сигнале { брк ), а уравнение Риккати сведется к линейному уравнению, которое имзет аналитическое решение, что позволило получить линейную зависимость мваду( 8р ) м ( 2? ), которая по существу является матрицей производной.

Полученные решония прямых и обратных задач потребовали разработки алгоритмов интегрирования уравнения Риккати на комплексной плоскости, сысс!;озф|«ктивных процедур впчкелонял модифицкропашшх фугашкй Бесселя от комплексного аргумента, квадратур интегрирования осцилируюад функ1^й. •

Б третьей главе описывается математическое обеспечение дпь определения УЗС по' комплексу зондов апэктрокаротака на основа

мщений прямых и обратных задач.

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

В качества двумерней _ прямой задачи в математическое >беспвчещ1е включены программы на основе метода прогстеси коэффициентов Фурье (В.Л.Друскин), математическая постановка и шеленная реализация проведена В.А.Пантюхиным (ВНИГИК).

Для расчета двумерной матрицы производных и моделирования сриьых КС градиент и потенциал-зондов на начальных итерациях разработаны БПЗ на основа фувдамэнталышх решений для моделей с 1эртакалькой неоднородностью и рядизльно неоднородной среда. Суть шгоритаа заключается в следуккцем. На первом этапе решается фямпя задача и вычисляются матрицы производных для :'здиплыю-неоднородней срода, а затем на основе полученных фиведешшх значений сопротивления пластов решается задача для «оделэй с вертикальной неоднородностью. В работе получены »звисимости, рассчитанные на основе ДПЗ, позволяющие корректировать кривыо КС и МП,полученные на основе БПЗ.Приведены примеры соррокции кривых КС для градизнт-зовдов в тонкослоистом разрэзз и зопоставченив о расчетами по ДПЗ,показано, что алгоритм юделирования измерений и расчета МП позволяет уверенно решать даумарную обратную задачу.

Использование низких частот тока питания зондов серийной яшаратури ИК позволяет разработать более простой, чем для зондов

БКЗ алгоритм вычисления измерений и двумерной матрицы производных. На основе формул, учитывающих скин-эффект в дальней зоне (Кауфман A.A..Плюснин М.И..Соколов В.П.) получено выражение для модели с изменением проводимое™ в радиальном и вертикальном направлениях, причем в континуальной постановке с использованием прямых задач для моделей с произвольным изменением проводимости по глубине.полученных во второй главе.

Двумерная матрица производных определяется как произведение МП* для вертикально-неоднородной среда на МП для модели с радиальной неоднородностью. По точности и времени счета алгоритм удовлетворяет требованиям практики.

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

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

Расчеты по программе моделирования зондов ЕК сопоставлялись с известными палаточными зависимостями для трехэлектродного и двухзондового ыногоэлоктродного БК (БКС-2), а также с кривыми профилирования. Во всех ситуациях он обеспечивает достаточную для практических целей точность.

Таким образом, разработано математическое обеспечение

моделирования измерений БКЗ, ИК, БК,ио точности удовлетворяющее требованиям необходимым при итерационном подходе к определении сопротивлештя.Получешшй алгоритм счисления двумерной МП обеспечивает уверенную сходимость итерацкошгого процесса и позволяет значительно сократить время счета.

В четвертой глазе приводятся результаты опробования и использования программного обеспечения определения удельного электрического сопротивления пс комплексу разнотипных зондов олектрокоротака в тонкослоистом рэзрпзо (ЕЬЕЗТП).

Математическое обеспечение позволяет проводить интерпретацию в слодуыцих технологических режимах:

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

-Задаются границы пластов, признаю! коллекторов и определяются параметры пластов ( Р31П, Рп{. 1>:зп1 ).

-Задаются границы пластов, "признаки коллекторов, и значение диаметра зоны проникновения в каждом пласте, определенного по косвенным методам,определяются сопротивления пластов и зон проникновения ( рп4, р зп1)-

-Предусматриваются задания сопротивлений зон проникновения (Рзп1), определенных по коротким градиент-зондам, либо МЕК, либо вычисленные на основе пэчрофизичвсюи зависимостей,определяются

Рп» • !>зт'

-В комплексе программ предусмотрен режим задания УЭО про-пластков, если оно известно по друг™ оценкам. Например, сопротивление глин, плотных высокоомных прослоев. В такой постановке решается задача меньшей размерности, повышается достоверность определения УЭС коллекторов, что особенно вакнс в тонкослоистом

разрезе.

Моделирование кривых КС для заданной модели геоэлектрического разреза позволяй'.' оценить качество интерпретации по более простым алгоритмам и г.рп необходимости уточнить УЭС но программе ELKSTR.

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

Во всех случаях подтверждена достаточно быстрая сходимость алгоритма и приемлемое время счета. Для точных данных достигается дисперсия порядка Iü'^-lü"4- за Ьь итераций и 2-4 обращения к ДОЗ, при этом параметры цласгор восстанавливаются точно.

Для данных с • помехами дисперсия, как правило, достигает величины того же порядка, как и сами помехи. В этом случае необходимо от 1-го до 3-х обращений к. ДПЗ.

■ Комплекс программ KLESTK прошел широкое опробование при интерпретации на реальном скьажгапюм материале. Исследовались различные возможности:

-интерпретация при отбивке границ пластов по БК2; -отбивка границ пластов с учетом других методов (в первую очередь БК);

-интерпретация материалов для произвольной модели (поточечная обработка);

-использование полного и ограниченного комплекса аондов; -использование характерных точек, снятых с кривых КС; -априорное задание диаметра зоны проникновения, поиск

параметра пластов для случая определения общего диаметра зоны проникновения для пачки пластоп;

-комплексное использование с более простыми алгоритмами (для мидоли пласта бесконечной толщины).

Результата опробовании по,, липли, что при качественном материпло дисперсия, кик правило, достигает 0,02-0,СО за 3-5 иторпций. Ирак*» теки во псих случаях достаточно одного обращения к Д113. Ислюльпаипнко модели с произвольным распределенном сопротивления по глубине нояшллот достигнуть оолоо меньших диспорсий, чпк для модолай о границами пластов, но в такой постановки требуются значительные затраты времени ЭВМ.

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

Важным резервом экономии вромони счета, является использование не всех точок с кривой КС, а некоторых характерных значений с изо, которые могут бить выбршш программой, либо заданы интерпретатором. На последней итерации НПЗМОИЮ уточнить результаты интерпретации по всей кривой НС. Преподанные исследования показали, что использование характерных точек позволяет значительно уменьшить вромя счета без потери точности в онрадолоши УЭО пластов.

Наиболее оЭДйчтивмо использовать комплекс программ возможно в следующем технологическом режиме. На первом этапе проводится интерпретация материалов ЭК по программам для пласта бесконечной толщины, требующих незначительных затрат времени ЭЕ.Ч, затем по программе ШБТП оценивается, качество интерпретации и при

необходимости уточняется сопротивление тох интервалов, где качество неудовлетворительно (большее расховдение фактических и теоретически расчитанных кривых КС, либо значительнее погрешности в УЭС).

' По. результатам опробовшшя показана работоспособность комплекса программ и ого преимущество при определении УЭС » тонкослоистом разрезе за счет использования интерпретационной модели, описывающей реальный геоэлектрический разроз, возможности оценки качества интерпретации и вычисления погрешностей ояро-деления УЭС.

Заключение

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

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

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

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

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

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

-Исследованы вопросы интерпретации для модели пачки пластов с проникновением (тонкослоистый разрез).

-Прадложона технологическая схема определения УЭО в тонкослоистом разрезе, которая позволяет повысить дбстогсрность определения сопротивления и уменьшить тробуомое время счота за счет 1гредварительной оцонки £>зп, поиска общего для пачки плаотов диаметра зоны проникновения, задания сопротивления пропластков, использования характерных точек.

-Проведена адаптация программного обеспечения для персональных компьютеров типа И5М-РС и их подключение к системам АСОИГИС, ('ИНТЕЛ, ПОДСЧЕТ, ИНГИС, отработала методика его использования на материалах различных регионов.

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

I.Численное решение прямой задачи теории индукционного каротаяа при радиолою - неоднородной зоне проникновения.// Прикладная геофизика, Вып. 108, М., Недра, 19Й4, о. ПУ-123. (соавтор Клеллер Л.Е.)

2.Быстрый алгоритм решения прямых задач бокового каротака. //Прикладная геофизика Вып. 112, М., Недра, 1984, с. 141-144. (соавтор Кнеллер Л.Е.)

3.Решение прямой и обратной задач индукционного каротажа для повышения точности определения сопротивления пластов при наличии вертикальной неоднородности вмещающих пород. - сб.тезисов"Повышение эффективности геофизических методов исследования скважин".Октябрьский,1988,с.58-59

4.Опыт разработки программного обеспечения длл интерпретации данных элоктрокаротажа на основе прямых задач и итерационных мотодов оптимизации.-в сб.тезисов "Повышение эффективности геофизических методов исследования скважин". Октябрьски!!, 1930, с. Ь0-Ь7.

5.Метода ■решения прямых и обратных задач элоктрокаротажа. М., 1909, с.53(-Мат. методы и автоматиз. системы в геологии:Обзор. ВНИИ экон. минер, сырья и геол.-развед. работ ВИЭМС.) (соавторы Кнеллер Л.Е. .Пантнхин В.А)

6.Решение прямой и обратной задач электрокаротажа для радиально-неоднородных сред.// Геология и геофизика, 1989,

& I, с. 88-96. (соавтор Кноллер Л.Е.)

7.Обратная задача электромагнитного каротажа для сред с радиальным градиентом электропроводности.//Передовой научно-производственный опыт.-Информационный сборник ВИЭМС,выл.8,с.11-19. (соавтор Кнеллер Л.Е.)

8. Решениэ прямой и обратной задач электрокаротажа для сред с радиальным градиентом электрического сопротивления//Г1е ре до-вой научно-производственный опыт.-Информационный сборник, ВИЭМС.еып,8,с.З-И. (соавтор Кноллер Л.Е.)

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

кого разреза'.// Нефтегазовая геология и геофизика. ЭИ. ВНШОЭНГ М.1990,ВЫЛ.7;С.34-40..

( соавторы Кноллер Л.Е..Кнеллер О.М.)

10.Решение прямой и обратной задач индукционного каротажа о учетом вертикальной и радиальной неоднородности геовлектрического разреза.//Изв.Вузов.Геология и разведка * 9 1990,0.95-102.

(соавтор Кнеллер Л.Е.)

11.Решение прямой и обратной задач ндукционного каротажа для сред с произвольным и дискретным распределением проводимости по глубина.//Теология и геофизика * Б, 1990, с. 112-132.

(соавтор Кнеллер Л.Е.)

12.Определение сопротивления в тонкослоистом разрезе по данным электрокаротажа на основе решений прямых и обратных задач(метода,алгоритмы,ре зультаты)М.1991,Деп.ВНИТИ, * 2255-В91 с. 91 (соавторы Кнеллер Л.Е.,Кнеллер О.М.)

Подписано к печати 4.08.92.Тира* 100. Ротопринт ТатПИПИнвфть,Бугульма,М.Джалиля,32.