Методическое руководство по интерпретации данных ГИС. Руководство к практическим занятиям по курсу "Автоматизированная интерпретация данных геофизических исследований скважин " для студентов IV курса
 Скачать 1.56 Mb.
|
|
Автоматизированная интерпретация данных геофизических исследований скважин (на примере системы ГИНТЕЛ) Методическое руководство к практическим занятиям по курсу “Автоматизированная интерпретация данных геофизических исследований скважин ” для студентов IV курса КФУ, обучающихся по направлению “Геология” (бакалавр геологии, специализация геофизика) На сегодняшний день большинство производственных операций производится при помощи вычислительных машин. Не исключение и интерпретация данных геофизических исследований скважин (далее ГИС). В данном методическом пособии рассмотрены основные вопросы автоматизированной интерпретации данных ГИС (далее АИГИС). В настоящее время существует большое количество различных программных пакетов для интерпретации данных ГИС: ГИНТЕЛ, ПРАЙМ, TechLog и пр.  Рис.1 Задачи АИГИС. По существу говоря, АИГИС есть интерпретация данных с использованием IT. Т.е. использование компьютерных инструментов для облегчения расчётов и принятия решений человеком (интерпретатором). И здесь нет разницы в алгоритмах и подходах к интерпретации при “ручной” работе и на ЭВМ. Поэтому надо иметь чёткое представление об интерпретации ГИС вообще. Задачи интерпретации данных ГИС можно разделить на задачи общего характера и задачи детального исследования. Задачи общего характера:
Задачи детального исследования - количественное определение параметров, таких как:
Приведенные задачи решаются в 3 этапа: 1. Индивидуальная интерпретация диаграмм каждого геофизического метода в отдельности, с целью выделения в разрезе пород отличной от вмещающих пород геофизической характеристикой данного геофизического поля; 2. Комплексная интерпретация. На этом этапе комплексно (совместно) интерпретируются диаграммы разных геофизических методов по одной и той же скважине, при этом учитываются геологические и другие (промысловые, технические) данные. В результате составляется геологический разрез, и выделяются промышленно интересные объекты (продуктивные горизонты); 3. Обобщающая интерпретация заключается в анализе данных ГИС по площади. В той или иной мере все эти задачи автоматизированы. Однако достичь полной автоматизации процесса интерпретации на современном уровне развития техники и технологий пока не удаётся. Поэтому принято говорить об автоматизированной, а не об автоматической интерпретации. Основная проблема, с которой столкнулись геофизики при разработке автоматизированных систем интерпретации это формализация задач интерпретации, т.е. описание геологических процессов, явлений и объектов языком математики (создание математической модели). Вычислительная машина способна решать только конкретные формализованные задачи в виде алгоритмов. На данный момент компьютеры позволяют производить громоздкие вычисления за короткое время. Это сокращает временные затраты на интерпретацию. Однако, машина не способна контролировать качество получаемого результата. Поэтому задача контроля произведенной машиной интерпретации и ответственность за результат лежит на человеке, проводящем интерпретацию. Исходными данными для автоматизированной интерпретации служат:
 Рис.2. Примеры форматов данных ГИС
LAS-файл это общепринятый в скважинной геофизике международный стандарт записи данных ГИС. Пример LAS-файла с обозначением его составных частей показан на рис.3. Кривые ГИС записаны в цифровом формате, то есть, представлены в виде записей мгновенного значения геофизического параметра пошагово с глубиной. В России расстояние по вертикали между двумя последовательными записями обычно принимается равным 0.2 м (на зарубежных диаграммах – через фут - 0,33м). Это расстояние называется шагом квантования по глубине.  Рис.3. Пример LAS-файла. Геофизический сигнал можно представить в виде суммы полезного сигнала и случайного шума. В скважинной геофизике выделение в разрезе пласта связано с наличием аномалии на зарегистрированной кривой данного метода ГИС. Для уверенного выделения аномалии, вызванной полезным сигналом на фоне шума, необходимо, чтобы в аномалию было вовлечено не менее 3 последовательных точек. Другими словами, для вычленения в разрезе пласта необходимо, чтобы он выделялся не менее чем на 3 точках подряд, то есть минимальная технологически обоснованная мощность пласта составляет - 0.2м * 3 = 0.6м. Исключение составляют высокоразрещающие методы (акустический сканер, микрозонды и пр.), имеющие меньший шаг квантования по глубине. В процессе “ручной” интерпретации последовательно проводятся следующие шаги (этапы). В общем случае (например, в терригенном разрезе) в интервале детальных исследований выделяются пласты, среди них по принципу плотный - неплотный и глинистый - неглинистый выбираются неплотные (пористые) и неглинистые (проницаемые). Среди этого набора пластов, интересующих нас с точки зрения хороших ФЕС, уже по результатам количественных расчётов мы выбираем только те, которые можно отнести к коллекторам (по критическим значениям для каждого пласта), в которых теперь также количественно определяем и насыщение. Процесс автоматизированной интерпретации, по сути, использует такой же граф. Но специфика применения компьютерных технологий требует некоторых особых приёмов, которые и разбираются далее. АИГИС можно разбить на несколько основных этапов: - создание проекта; - загрузка данных; - ввод априорной информации; - расчёт поправок в результаты измерений (кривые ГИС) за влияние различных факторов; - определение границ пластов (расчленение разреза); - расчёт ФЕС выделенных пластов (возможных коллекторов); - определение литологии, насыщения, типа коллектора; - расчёт модели породы выделенного интервала; - формирование заключения; - вывод графических и текстовых материалов в твёрдой копии; - дальнейшие работы по детальной интерпретации, сводной (площадной) интерпретации, создание постоянно-действующей модели (ПДМ) месторождения; В той или иной последовательности все эти этапы наличествуют во всех автоматизированных системах интерпретации данных ГИС, реализующих обработку открытого ствола. В автоматизированной интерпретации данных ГИС существуют два подхода к её реализации — интерпретация поточечная и попластовая. В поточечной интерпретации фильтрационно ёмкостные свойства (ФЕС) рассчитываются последовательно для каждой точки всего интервала обработки – результатом интерпретации являются непрерывные кривые полученных рассчитанных характеристик, не приуроченные к пластам ввиду не расчленённого заранее разреза. При попластовой интерпретации прежде выделяются пласты, затем внутри них снимаются пластовые значения параметра (среднее значение по пласту) и вычисляются ФЕС - результатом интерпретации являются кусочно-непрерывные кривые полученных характеристик, приуроченные к пластам ввиду расчленённого заранее разреза. Оба подхода имеют равные права на существование. В системе ГИНТЕЛ (основная, базовая версия – рис.4.) заложен поточечный принцип интерпретации, но представленный в данном методическом руководстве вариант интерпретации (так называемая ТАТВЕРСИЯ ГИНТЕЛ или ТатТехнология) использует попластовый вариант, т.к. алгоритм интерпретации, принятый на месторождениях Татарстана (Стандарт предприятия) использует именно такой подход.  Рис.4. Функциональные возможности системы ГИНТЕЛРассмотрим процесс автоматизированной интерпретации данных ГИС в системе ГИНТЕЛ 2005 (на примере скважин РТ).  Рис.5. Вид основного монитора ГИНТЕЛ 2005
Для работы в системе ГИНТЕЛ каждому пользователю необходимо создать свою учетную запись. Делается это для сохранения данных в случае, если за одним компьютером работает несколько пользователей.
Необходимым шагом создания своего проекта является разработка базы данных. База данных содержит в себе сведения о регионе, площади и конкретной скважине, по которой будет проводиться интерпретация. Это необходимо для того, чтобы компьютер в автоматическом режиме находил в своей памяти петрофизические зависимости, выявленные (построенные) и заранее загруженные в неё для выбранного региона, площади и конкретного пласта. Подробнее остановимся на вопросе иерархии базы данных. База данных ГИНТЕЛ иерархична, то есть имеет некую степень подчинения одних элементов другим. Высшим в иерархии стоит регион. Под регионом понимается крупная нефтегазовая провинция или административный регион (например, Западная Сибирь или Татарстан). В базе данных может содержаться несколько регионов.  Рис. 6. Создание базы данных. В каждом регионе содержатся площади. Под площадью понимается небольшое месторождение (или его часть), которое в процессе интерпретации рассматривается как единая система с конкретными петрофизическими зависимостями для всей системы в целом. Наконец, в каждой площади находятся скважины. Скважина в ГИНТЕЛ это низшая ступень иерархии базы данных. В ней содержатся кривые ГИС и прочие данные. Стоит отметить, что разные стволы одной скважины рассматриваются как отдельные скважины.
 Рис.7. Разработка проекта.
Следующим шагом является загрузка исходных LAS-файлов в систему (Рис. 8). LAS-файлы содержат информацию о типе кривых, их интервале, шаге квантования и собственно скважинные данные. Программа автоматически определяет тип кривой по её мнемонике. (Под мнемоникой мы понимаем уникальное имя каждой кривой, например SP, GR, Dt и пр). Так же из LAS файлов считывается интервал кривой. В дальнейшем интерпретатор, основываясь на этой информации, устанавливает рабочий интервал проекта, то есть интервал, в котором будет производиться интерпретация. После загрузки LAS-файлы преобразуются во внутренний формат ГИНТЕЛа с расширением *.well. В дальнейшем все операции производятся именно над этими данными.  Рис.8. Загрузка Las файлов. 5. Ввод данных инклинометрии  Рис.9. Ввод данных инклинометрии. Данные инклинометрии – это информация о пространственном положении скважины. В автоматизированных системах, к сожалению, нет единого формата представления данных инклинометрии, но мы рассмотрим наиболее существенные особенности её. Инклинометрия представляется в виде таблицы с некоторым набором столбцов. Число их может изменяться от 3 до 11, однако значимыми являются только 5. Важно отметить, что имена столбцов не имеют единого стандарта. Рассмотрим наиболее важные из них (далее в скобках будут приводиться возможные имена столбца). Глубина (измеренная глубина, относительная глубина, кабельная глубина, MD) – столбец указывает на длину кабеля, то есть на расстояние от устья скважины до текущей точки в метрах. Как правило, шаг глубины 10 или 20 м. Угол (Альфа, наклон скважины, α, Incl) – столбец указывает угол наклона скважины. Указывается в градусах или в радианах. Магнитный азимут (Азимут, Az) - столбец указывает угол между направлением скважины и направлением на магнитный полюс в горизонтальной плоскости. Указывается в градусах или в радианах. ΔX (+Восток – Запад, dX) – столбец указывает на смещение ствола скважины в направлении X. Положительное значение указывает на смещение в восточном направлении, отрицательное – на смещение в западном. ΔY (+Север – Юг, dY) – столбец указывает на смещение ствола скважины в направлении Y. Положительное значение указывает на смещение в северном направлении, отрицательное – на смещение в южном. Для построения ствола скважины достаточно 3 столбцов (Глубина, Угол и Азимут или Глубина, ΔX и ΔY). 6. Анализ базы данных по скважине. Следующим шагом является проверка корректности загруженных данных в табличном (Рис.10) и графическом (Рис.11) виде. Главным окном доступа в базу данных по скважине является диспетчер базы данных (Рис.10). Здесь отображается вся информация по конкретной скважине в виде таблицы. Столбцы указывают Имя, Префикс, Интервал, Шаг и направление данных. С именем WELL хранится информация о скважине. С именем ListZones хранится информация о списке зон проекта (объектов интерпретации). С именем INKL хранится таблица данных инклинометрии. С именем CORE хранится таблица данных керна. Далее с именами мнемоник хранятся кривые ГИС.  Рис.10. Проверка корректности загруженных данных в табличном виде. Диспетчер данных по скважине. Проверка данных в графическом формате проводится на планшете данных по скважине. В системе ГИНТЕЛ есть набор шаблонов планшетов, также существует возможность создавать свои собственные планшеты и шаблоны. На планшете отображаются кривые ГИС, Измеренные и Абсолютные глубины. Существует возможность отображать кривые в линейном и логарифмическом масштабе параметра. На планшете можно производить различные операции для более наглядного отображения кривых – добавление и удаление колонок или отдельных кривых, изменение масштаба геофизического параметра, перемещение кривых из одной колонки в другую.  Рис.11. Проверка корректности загруженных данных в графическом виде. Планшет данных по скважине. 7. Создание документа TDoc. Далее рассмотрим ТатТехнологию – алгоритм интерпретации данных ГИС, разработанный в ООО «ТНГ-ГРУПП» и широко применяющийся в Татарстане. Первым шагом является создание документа TDoc, содержащего априорные данные, необходимые для интерпретации. К ним относятся данные: – по скважине: код площади, номер скважины, кровля, подошва интервала обработки, координаты устья, УЭС и плотность раствора, УЭС фильтрата и т.д.; – о методах исследования: скорость их регистрации, типы приборов, длины зондов, цена условной единицы и одного микрорентгена в час для НГК и ГК, время интегрирования (интегрирующая ячейка); – об опорных пластах (или их отсутствии); – о границах стратиграфических горизонтов. Рассмотрим структуру файла по частям. 1) СКВАЖИНА: 1 ПЛОЩАДЬ: Примерная КРОВЛЯ: 800.00 ПОДОШВА: 1178.00 $FIELD 214 ; вставить код площади В данной части файла указываются: номер скважины, название площади, отметки кровли и подошвы интервала обработки, а также указывается код площади, по которому программа в памяти компьютера находит заранее записанные петрофизические зависимости для этой прощади. 2) Параметры ПЖ $DM 1.12 ; Плотность ПЖ,г/см3 $RM 0.4 ; УЭС ПЖ,Омм $VIS 28 ; Вязкость ПЖ,с $LOS 8 ; Водоотдача,см3/с $DD 4 ; Дзп/Дс $DN 0.216 ; Дс номинальный $ZP 0 ; Проникновение: 0 - повышающ., 1 - понижающ. $TR глинистый ; Тип раствора: глинистый, известковый и т.д. $TZ 40 ; Температура на забое, град.Ц $DT 14.09.08 ; Дата каротажа (чч.мм.гг) Определение характера проникновения: подошва тип: 0 - повышающ., 1 - понижающ. $CP ; CharacterPenetration (характер проникновения) 1178.0 0 $END_CP Здесь указываются параметры промывочной жидкости, такие как: плотность, удельное электрическое сопротивление, вязкость, водоотдача, отношение диаметра зоны проникновения к диаметру скважины, номинальный диаметр скважины, тип проникновения, температура на забое и дата каротажа. 3) Параметры аппаратуры Метод Ед.изм Шаг_кв Скор_рег Тип_приб Тип_зонд Пост_инт Ц_усл Отн_кнл $TOOL NGL усл.ед. 0.2 400 РК-3-76 - 3 15400 0.11 $TOOL GR мкр/ч 0.2 400 РК-3-76 - 3 1972 0.11 $TOOL ILD Омм 0.2 2500 4И1 4F0.75 - - - $TOOL LLD Омм 0.2 2500 БК-3 БК-3 - - - $TOOL SP мВ 0.2 2500 K-1 - - - - $TOOL PZ Омм 0.2 2500 K-1 М11N0.5A - - - $TOOL CALI м 0.2 2500 K-1 - - - - В данной части указываются параметры аппаратуры: используемый метод (в приведённом списке методов перечислены те, которые входят в обязательный комплекс ГИС на месторождениях РТ), единица измерения регистрируемого геофизического параметра, шаг квантования, скорость регистрации, тип используемого прибора, тип зонда, постоянная интегрирования, цена условной единицы (для РК) и отношения каналов (для РК). С учетом этих параметров программа выбирает необходимую информацию из памяти компьютера и вводит нужные поправки в геофизические кривые для интерпретации. 4) Стратиграфия Вместо * вставить подошву горизонта $STRAT p 1028.0 Подольский k 1092.8 Каширский vrt 1114.6 Верейский(терриген.) vrk 1143.2 Верейский(карбонат.) bash 1178.0 Башкирский prt * Протвинский mich * Тар-Стеш_Мих-Венев alex * Алексинский al * Алексинский и Тульский tula_k * Тульский(карбонат.) tula_t * Тульский(терриген.) аbobr * Бобриковский и Тульский bobrik * Бобриковский rd * Радаевский t * Верхнетурнейский подъярус kzl * Кизеловский cher * Черепетский up_mal * Упино-Малевский zvl * Заволжский d_l * Данково-Лебедянский el * Елецкий zd * Задонский vr_evl * Воронеж_Евлано-Ливен br * Бурегский sm * Семилукский srg * Саргаевский d * Терригенные отложения Девона d0 * Кыновский(Д0) d1 * Пашийский(Д1) d2 * Муллинский(Д2) d3 * Ардатовский(Д3) d4 * Воробьевский(Д4) d5 * Эйфельский ярус(Д5) $END_STRAT В разделе «стратиграфия» указываются горизонты исследуемой территории. Вместо звездочек записывается их подошва, для того, чтобы программа использовала известные петрофизические связи и зависимости для данных горизонтов, заранее записанные в памяти компьютера. Очевидно, что разбивка разреза на стратиграфические горизонты требует предварительной работы над диаграммным материалом комплекса методов ГИС. 5) Таблица опорных для РК $TABL_OP НГК_ДЕВ -9999 ;NGKmin НГК_УГЛ -9999 ;NGKmin НГК_ВЕР 1.1 ;NGKmin НГК_ФАМ -9999 ;NGKmax НГК_АЛ 3.8 ;NGKmax ГК_НЕГЛ 1.5 ;GKmin ГК_ДЕВ -9999 ;GKmax ГК_СЕМ -9999 ;GKmax ГК_УГЛ -9999 ;GKmax ГК_ВЕР 17.0 ;GKmax $END_TABL_OP Для горизонтов выделяют опорные значения максимума и минимума по замерам РК. Программа выбирает их автоматически и записывает в данную таблицу файла TDoc. Интерпретатор имеет возможность редактировать значения по опорным горизонтам. В дальнейшем эти данные используются при расчёте двойного разностного параметра (Ag и Ang). Оставшаяся часть файла относится к описанию конструкции скважины. |